换流变压器空负载损耗仿真:Simulink建模与参数计算全流程
2026/9/8 21:14:11 网站建设 项目流程

作为一个常年和电力电子仿真打交道的人,我最近刚好把换流变压器的空负载损耗仿真模型完整搭了一遍,从参数计算到Simulink建模再到结果分析,整个过程踩了不少坑。这篇文章我把当时的思路、参数计算、建模步骤和排查过程整理出来,希望能给正在做直流输电或变压器仿真的朋友一些参考。

先说清楚这个东西到底解决什么问题。换流变压器和普通电力变压器不一样,它在直流输电系统里承担着换流阀与交流系统之间的电压变换和电气隔离作用,运行工况非常复杂:阀侧电压电流含有大量谐波,直流偏磁风险高,绝缘结构特殊。损耗仿真如果做得糙,后面热设计、效率评估、经济性分析全部要返工。所以损耗仿真这件事,看起来只是算两个数(空载损耗P0和负载损耗Pk),实际上一旦做细,牵扯到谐波、饱和、频率特性,工程量并不小。

本文适合刚接触Simulink电力系统仿真的学生,也适合已经在做直流输电相关研究的工程师。我会从模型思路、参数计算、仿真搭建到问题排查完整走一遍,过程中所有关键参数都给出计算过程,方便你直接套用。

1. 模型整体思路与方案选型

1.1 为什么用Simulink而不是自编程序

换流变压器损耗仿真听起来可以写个解析公式算一算,但实际工程场景里没这么简单。空载损耗取决于铁芯的磁化特性和铁耗特性,负载损耗取决于绕组电阻和漏抗,这两部分在谐波环境下还要考虑频率的影响。如果单纯用公式估算,谐波电流作用下的附加损耗很难算准。而Simulink/Simscape的优势在于,它可以构建完整的电路拓扑,用电磁暂态的视角去看变压器的工作状态,谐波的影响是自动被包含在波形计算里的。

另外,Simulink生态里有现成的变压器模型。Simscape Electrical库中的Specialized Power Systems(SPS)提供了单相和三相变压器模型,包括线性变压器、饱和变压器等,这些模型内置了磁路饱和、铁耗计算能力,对于损耗仿真来说已经非常够用。比起自己写差分方程求解磁链电流关系,用现成模型至少能节省一周时间。

我的总体思路是:用SPS库中的三相双绕组饱和变压器作为核心,在原方和副方分别接电源和负载,在关键节点接入电压电流测量模块,最后用有功功率测量模块分别测出空载和负载工况下的损耗值。整个模型结构清晰,后续接整流器、做谐波分析也非常方便。

1.2 空载损耗和负载损耗的物理本质

这两个概念我先用大白话说清楚,因为很多初学者一上来就被公式绕晕了。

空载损耗,又叫铁耗,是变压器一次侧加额定电压、二次侧开路时消耗的有功功率。它主要由磁滞损耗和涡流损耗组成,取决于铁芯材料的磁化特性和运行频率。磁滞损耗正比于频率和磁滞回线面积,涡流损耗正比于频率的平方和磁密的平方。这两个分量加起来,最常用的工程近似就是Steinmetz公式:

P_fe = k × f^α × B_max^β

其中k是材料常数,α约为1~1.6,β约为1.6~2.0,不同硅钢片差异很大。Simulink的变压器模型里,铁耗就是按类似逻辑计算的,所以参数能不能填准,直接决定空载损耗仿真准不准。

负载损耗,又叫铜耗,是变压器二次侧短路(或通过额定负载)时、额定电流流过绕组产生的有功损耗,主要由绕组电阻的焦耳热决定。基本公式就是P = 3I²R(三相)。谐波存在时,绕组还会出现额外的涡流损耗,这部分和频率的平方相关,在换流变压器场景里尤其不能忽略。

我在仿真里先只考虑基波下的空载和额定负载损耗,把谐波附加损耗作为后续扩展项,这一点在后面5.1节会再展开。

1.3 方案选型:连续还是离散仿真

用Simulink做电力电子仿真,首先要决定用连续求解器还是离散仿真。Simscape的Electrical库推荐用离散求解器,尤其是模型中带有电力电子开关器件时。离散仿真用固定的采样步长,避免了连续求解器在开关切换时刻的迭代收敛问题,速度也快很多。虽然我现在做的是纯变压器空负载仿真,没有电力电子开关,但我还是建议从第一天就用离散模式,因为换流变压器的应用场景终究是要接换流阀的,早点定下离散仿真的框架,后面不用重新搭模型。

采样时间的选择也讲究。我的经验是先粗后细:先用Ts=50μs跑通流程,确认拓扑正确、数值稳定,再加密到5μs对比结果。对50Hz工频系统来说,50μs的离散步长代表每个工频周期采样400点,对于损耗计算够用;如果涉及更高次谐波(比如关注11、13次谐波),就适当加密。

2. 核心细节解析与实操要点

2.1 Simulink变压器模型选型:线性还是饱和

SPS库里的变压器模型分几种:单相变压器、三相双绕组变压器、三相三绕组变压器,每种又分线性(Linear)和饱和(Saturable)选项。损耗仿真想准确,理论上必须用饱和模型,因为铁耗的模拟依赖于磁路的非线性特性。但饱和模型的参数多一些,多了一个励磁支路的饱和特性曲线需要设置。

我实际搭建时是分两步走的:第一步先用线性变压器模型,把整个仿真链路跑通,确认电压电流关系正确,损耗的大致数值量级对得上;第二步再把线性模型换成饱和模型,填上铁耗系数和磁化曲线,仿真结果下来差异其实不小,尤其是空载电流波形,饱和模型下是明显的尖顶波,线性模型下是正弦波,实际变压器确实是尖顶波。这一步做完,我对仿真结果的信服度才上来。

2.2 变压器模型内部如何计算损耗

这里要稍微深入一下模型内部,因为如果不理解它怎么算的,参数填错了你根本发现不了。饱和变压器模型内部把变压器等效为电感-电阻-理想变压器组合:原副方绕组分别有电阻R1、R2和漏感L1、L2,励磁支路由励磁电感Lm和铁耗电阻Rm并联组成。

空载损耗就落在铁耗电阻Rm上,这个电阻的取值是根据你填的铁耗特性反推的。负载损耗则主要落在R1和R2上,谐波情况下还要考虑频率对电阻的影响,但SPS基础变压器模型默认是恒定电阻,所以如果你想模拟谐波附加损耗,需要自己扩展,后面会讲。

我在仿真时习惯在变压器输出端接一个三相有功无功测量模块,同时用电压电流测量模块输出瞬时波形,两头互相印证。只靠有功功率模块读数有时候会被数值噪声干扰,配合瞬时波形确认稳态情况,结果更可靠。

2.3 测量模块与数据读取的要点

损耗仿真最怕的就是“算出来的数不敢信”。为了拿到可信数据,我一般会同时用两种方式读取损耗。

第一种是用SPS库里的“Active & Reactive Power”模块,它直接输出三相有功功率,把空载工况下读到的有功就是P0,负载工况减去空载损耗就是Pk。第二种是用电压电流有效值自己算:P = 3 × V × I × cosφ,三相平衡时用单相数据乘3。这两种方法相互校验,基本能确定模型和测量环节没出错。

读取数据时的另一个坑是测量窗口。仿真启动瞬间会出现明显的暂态过程,如果从0时刻开始取平均,损耗值会被暂态能量“污染”,数值失真。我习惯在仿真时间轴中段选一个稳定的时间窗口,比如0.1秒到0.3秒,只统计这个窗口内的平均功率。这个细节看着小,其实影响很大,能产生百分之几到百分之十几的偏差。

2.4 仿真求解器与步长设置建议

求解器配置这块,虽然听起来老生常谈,但对损耗仿真精度的影响是最直接的。Simscape Electrical模型在连续仿真下,推荐使用ode23tb或ode15s这类刚性求解器;如果使用离散模式,需要设置采样时间Ts。我本人更偏向离散模式,原因上面说过——为了以后接换流阀模型。

如果你用的是离散模式,那求解器选“discrete”(离散),不需要像连续系统那样选具体算法。唯一需要认真设的就是采样时间。我实测过:同一个空载损耗模型,Ts=50μs时仿真结果和Ts=5μs相比,差值约1%左右,对于工程预研完全够用。如果Ts继续降低到1μs,性能收益已经很微弱,但仿真时间会翻倍。建议你根据模型复杂度和精度需求折中选取,不是越小越好。

3. 实操过程与核心环节实现

3.1 从零搭建换流变压器损耗仿真线路

我以一台三相双绕组换流变压器为例,参数如下:

  • 额定容量:100MVA
  • 电压比:220kV/35kV
  • 频率:50Hz
  • 联结组别:Ynd11(高压侧星形中性点引出,低压侧三角形)
  • 空载损耗:80kW(约定值)
  • 负载损耗:400kW(约定值)

完整的搭建步骤如下。

第一步,打开MATLAB,在命令行窗口输入simulink,新建一个空白模型。从Simscape → Electrical → Specialized Power Systems库中拖入“Three-Phase Transformer (Two Windings)”。这个模块是三相双绕组变压器,支持自定义联结组别,双击打开参数设置界面。

第二步,配置变压器参数。在“Configuration”选项卡里选联结方式Ynd11,勾选“Saturable core”启用饱和模型,把“Simulate core losses”选项打开。然后在“Parameters”选项卡里,填入额定容量、额定频率、原副方额定电压、绕组电阻和漏抗,单位用标幺值更省事。

第三步,搭电源和负载侧电路。原方接一个“Three-Phase Source”,电压设220kV、50Hz,内阻设很小(比如0.001Ω)。副方接“Three-Phase RLC Load”,做空载仿真时把负载设成开路(电阻设无穷大);做额定负载仿真时,把负载功率设为变压器额定容量100MVA。

第四步,加入测量环节。在变压器原方和副方分别插入“Three-Phase V-I Measurement”模块,连接Scope看波形;再在电源和变压器之间加一个“Active & Reactive Power”模块,用来读有功功率。

第五步,配置Powergui模块。模型里必须放一个“Powergui”模块,双击它选择离散仿真模式,设置采样时间Ts。我第一版用Ts=50μs,后续对比再加密到5μs。

3.2 变压器模型参数计算全流程

变压器模型参数是仿真的灵魂。我直接说计算过程,以标幺值为主:

已知额定容量S = 100MVA,高压侧额定电压U1 = 220kV,低压侧U2 = 35kV,那么:

  • 高压侧额定电流:I1 = S / (√3 × U1) = 100×10⁶ / (1.732×220×10³) ≈ 262.4A
  • 低压侧额定电流:I2 = S / (√3 × U2) = 100×10⁶ / (1.732×35×10³) ≈ 1649.6A

如果已知负载损耗Pk = 400kW(通常指75℃或145℃下的值),且认为高低压绕组电阻消耗的损耗等量(实际上可以按损耗分配系数,这里简化),那么:

  • 高压绕组电阻:R1 = Pk / (3 × I1²) = 400000 / (3 × 262.4²) ≈ 1.936Ω
  • 高压侧基准阻抗:Z1_base = U1² / S = (220×10³)² / 100×10⁶ = 484Ω
  • 高压绕组电阻标幺值:R1_pu = R1 / Z1_base ≈ 0.004

这个0.004左右的标幺值非常典型,实际变压器的绕组电阻标幺值大概就在0.002~0.01之间,如果你算出来是0.1,那大概率是单位换算出了问题。

短路阻抗(漏抗)的标幺值,一般变压器铭牌上会给短路电压百分比Uk%,一般换流变压器在12%~18%之间(为了限制短路电流和阀侧故障电流,比普通电力变压器略高)。假设Uk% = 15%,则:

  • 漏抗标幺值:Xk_pu ≈ Uk% / 100 = 0.15
  • 高压侧漏抗实际值:X1 = 0.15 × 484 = 72.6Ω,对应电感L1 = X1 / (2πf) ≈ 0.231H(这个值是高低压绕组漏抗折算到高压侧的总值,模型里如果允许分别填高低压漏抗,我一般各取一半0.075pu)

空载损耗P0 = 80kW,对应励磁支路。空载电流I0%一般在0.2%~0.5%,取0.3%:

  • 励磁电流标幺值:I0_pu = 0.003
  • 励磁支路导纳标幺值Ym_pu = I0_pu ≈ 0.003(忽略相位角时近似)
  • 励磁电导标幺值:Gm_pu = P0 / S = 80×10³ / 100×10⁶ = 0.0008
  • 励磁电纳标幺值:Bm_pu = √(Ym² - Gm²) = √(0.003² - 0.0008²) ≈ 0.00289

到这里,变压器模型参数就齐了。我在Simulink里填参数时,全部用标幺值填入,因为Specialized Power Systems的变压器模块默认支持SI和pu两种单位,用pu可以避免折算变比时的麻烦。

再说铁耗系数。饱和变压器模型里如果打开了“Simulate core losses”,会要求填Steinmetz系数:k、α、β。这个系数怎么定?最稳的做法是用铭牌空载损耗反推。比如额定频率50Hz、额定电压下的空载损耗是80kW,你可以先用k的初值凑,再微调让仿真结果贴近80kW。我第一版直接用默认系数算出来只有40kW,后来根据铭牌数据反推才修正过来,这一步不能省。

3.3 空载损耗与负载损耗仿真操作实录

空载损耗仿真

  1. 副方三相RLC负载设为开路,等效于不接负载。
  2. 原方接额定电压220kV的50Hz三相电源。
  3. 仿真时长设0.5s,求解器用离散模式,Ts=50μs。
  4. 仿真结束后,查看原方电压电流波形。正常情况:空载电流幅值很小(额定电流0.3%左右),波形是尖顶波,不是正弦波。这是铁芯饱和的最直观证据。
  5. 读取“Active & Reactive Power”模块的有功功率,用0.2s到0.5s之间的平均数值作为空载损耗。

按我们的参数,空载电流有效值应该在0.78A左右(262.4×0.003),实际上因为励磁电流非线性,可能稍微偏高一点,只要在0.5~1.2A范围内都合理。

负载损耗仿真

  1. 副方接额定负载,负载功率设为100MVA,功率因数可以取0.8或0.9。
  2. 原方仍然接220kV额定电压,仿真时间0.5s。
  3. 读取原方有功功率P_input,减去空载损耗P0,差值就是负载损耗Pk。

这里有个细节:我说的减法是基于“损耗分离”的近似——认为空载损耗只和电压相关,负载损耗只和电流相关。变压器实际运行中这两者不是完全解耦的,但在工程精度要求内这样处理完全够用。

我实测的结果:空载损耗仿真读数为81.2kW,和设定值80kW只差1.5%,误差主要来自励磁支路的非线性特性和测量窗口的暂态混入。负载损耗读数398kW,误差0.5%。对于工程前期评估,这个精度完够用。

3.4 损耗分离与结果可视化

仿真跑完之后,数据一定要整理成能“给别人讲清楚”的形式。我常用的方法是:

  • 在MATLAB工作区里用脚本读取仿真输出:
% 读取仿真输出数据 load('simout.mat'); active_power = simout.active_power.Data; % 有功功率数组 time = simout.active_power.Time; % 时间数组 % 只取稳态段(比如0.1s之后) steady_idx = time > 0.1; p0_measured = mean(active_power(steady_idx)); fprintf('空载损耗测量值: %.2f kW\n', p0_measured / 1e3);
  • 分别跑空载、额定负载两种工况,得到P0和Pk。
  • 再用损耗随负载率变化的公式:P_total = P0 + β²Pk,其中β是负载率。
  • 画一条“损耗-负载率”曲线,用Matlab的plot,既直观又专业。

换流变压器损耗仿真的最终目的通常有两个:一是给设计阶段做热校核(损耗决定温升),二是给运行阶段做效率评估(损耗决定经济性)。所以最后我会把不同负载率、不同功率因数下的损耗表导成Excel,方便后续接入经济性分析。

4. 常见问题与排查技巧实录

Simulink仿真变压器的坑,我基本都踩过。这里挑几个最高频的,写成速查表,你照着排查能省大量时间。

4.1 仿真速度慢得像蜗牛

现象:仿真跑起来极其慢,几乎每步都要卡一下,0.5s的仿真时长跑了20分钟还没结束。

原因

  1. 求解器设置不合理,用了ode45去解刚性系统(带电力电子或不连续环节)。
  2. 模型里有无穷大阻抗/导纳的悬空节点,导致代数环或数值振荡。
  3. 变压器模型饱和特性给的饱和点太密,导致插值计算慢。

解决办法

  • 改用ode23tb或ode15s,或者直接切离散模式。
  • 把Powergui切换为离散模式,采样时间从50μs起步,如果还慢就调到100μs看趋势。
  • 检查有没有悬空端口,尤其是变压器中性点、三角绕组不接地的那一侧,Simulink里悬空节点会引发数值问题,可以接地(Ground模块)或者加高阻接地。

4.2 仿真直接发散,波形飞到天际

现象:电压电流瞬间变成1e30这种离谱数值,波形直接拉满纵轴。

原因

  1. 变压器参数单位填错,常见的是把标幺值当国际单位填,或反之。
  2. 负载或电源初始相位设置冲突。
  3. 变压器励磁电感太小,导致励磁电流极大。

解决办法

  • 先检查变压器模块的所有参数是不是按铭牌或按上面计算过程填的,特别是短路阻抗百分比,填成0.15还是15,区别巨大。
  • 把仿真初始时刻的电压源相位设为0,负载初始电压不要给太大。
  • 把变压器饱和特性曲线砍掉,先按线性模型跑通,再加非线性。
  • 用“Powergui → Tools → Compute Initial States”让软件自动计算初值,这是一个救命功能。

4.3 空载损耗仿真结果比设计值差一半以上

现象:标称80kW的空载损耗,仿真出来只有30kW或者变成150kW。

原因

  1. 铁耗系数(Steinmetz系数)没填,或者填的是默认值,和实际硅钢片损耗特性不匹配。
  2. 空载电流取的百分数不对。
  3. 测量窗口覆盖了暂态启动过程,平均值被拉偏。

解决办法

  • 在变压器模块里仔细填铁耗三个系数:k、α、β(常数、频率指数、磁密指数)。如果只有一点实测数据,可以在额定频率、额定电压下用单点拟合:先按50Hz、1.7T工况校准k。
  • 测量时取稳态窗,比如0.1s之后,而不是全程平均。
  • 对比励磁导纳计算值和模块默认值,确保Y、G、B三者的量级一致。

4.4 谐波分析时FFT结果一团糟

现象:用Powergui的FFT Analysis看电流频谱,发现50Hz基波分量几乎为0,全是乱七八糟的间谐波。

原因

  • FFT窗口没对准整周期,频谱泄漏严重。
  • 采样时间不满足整周期采样条件。

解决办法

  • FFT分析前,把信号保存为“Structure with time”格式,在Powergui里设置好基波频率50Hz和分析窗口周期数(比如5个周期)。
  • 确保数据长度是周期的整数倍,比如0.1s(5个周期),而不是0.13s这种非整数周期。
  • 如果信号里有大量开关谐波,先经过一个低通滤波器,或者把采样率提高,避免混叠。

4.5 换流变压器特有的阀侧谐波问题

现象:阀侧电流波形明显畸变,谐波含量远超预期,损耗数值不稳定。

原因

  1. 阀侧连接的是整流器,电流含大量特征谐波(6k±1次,12脉波系统主要是11、13次),这些谐波在绕组里产生附加损耗(涡流损耗),模型如果不考虑频率特性,损耗会偏低。
  2. 换流变压器的漏抗和谐波阻抗特性与普通变压器不同,简单模型无法准确反映。

解决办法

  • 如果研究的是谐波损耗,可以考虑用Simulink里的阻抗-频率特性分析(Impedance vs Frequency Measurement),给变压器加扫频激励,得到频率响应后拟合频率相关的绕组电阻模型。
  • 简单做法:在绕组电阻中加一个频率修正系数,R(f) = R_dc × (1 + k×f²),k根据经验取1e-7左右(针对大型油浸式变压器涡流损耗比例),这样估算谐波附加损耗比较快。
  • 追求精确的话,直接上电磁暂态级建模工具,Simulink侧重系统级,变压器内部的集肤效应和邻近效应用专门软件更合适。

5. 这个模型还能怎么扩展

换流变压器损耗仿真模型最怕“搭完就扔”。站在我的经验角度,这个模型有很高的扩展价值,至少可以往以下几个方向延伸:

5.1 完整直流输电系统中的损耗评估

把换流变压器模型放进一个完整的直流输电系统里,整流器、逆变器、直流线路、无功补偿、滤波器都加上。这时变压器的损耗就不再是孤立评估,而是整个换流站损耗的一部分。换流站损耗评估里,变压器损耗占比通常在30%~40%,其余主要是阀损耗,所以精确的变压器损耗模型直接影响整个换流站效率计算。

我实际做过的做法是:保留变压器损耗计算模块,输入侧接整流桥,输出侧测量有功功率的变化,然后在Simulink里做一个“损耗在线计算”子系统,自动把空载损耗和负载损耗分离出来。做这个扩展的时候,你会发现变压器的谐波损耗因为整流器而显著增加,这时就要回到4.5节里的频率修正方法。

5.2 与热模型耦合,做温升预估

变压器损耗算出来不是终点,最终一定有人问你:“损耗这么大,温升怎么样?冷却器选多大?”所以我会把损耗计算结果作为输入,接到变压器热模型中。热模型用的是IEEE或IEC标准的油浸式变压器热路模型:顶层油温升、绕组热点温升、时间常数。Simulink里搭一个热路模型其实不复杂,损耗作为热源,热容和热阻用集总参数。

这样扩展后,整个模型就从“电-磁”扩展到了“电-磁-热”,在工程上的说服力完全不一样。我在评估换流变的负载能力(比如过负荷1.2倍能撑多久)时,这个耦合模型非常有用。

5.3 参数化扫描与优化

Simulink原生的参数扫描体验一般,推荐搭配Simulink Design Optimization,或者直接用MATLAB脚本循环跑不同参数组合。比如固定电压等级,扫描不同短路阻抗(12%~18%)、不同空载电流(0.2%~0.6%)下的损耗组合,生成损耗等高线图,这样在变压器招投标阶段可以快速给出不同方案的损耗对比。

脚本控制仿真、自动改参数、批量保存结果,我建议用sim命令写一个批量仿真脚本,这是把Simulink模型变成生产力工具的关键一步,避免每次手动改参数、手动截图、手动记录,效率至少提升一个量级。

6. 一点个人经验

做变压器损耗仿真这几年,我有两点体会最深。

第一,仿真模型的精度上限,取决于你对变压器物理的理解,而不是Simulink的熟练度。模块谁都会拖,但明白每个参数背后是什么物理量、量级应该在什么范围,才能判断仿真结果是对是错。我见过很多初学者拿着仿真结果直接写报告,一看空载电流5%,明显不合理,还分析得头头是道。所以跑完仿真一定要用理论值做一次数量级校验,这是从业者的基本素养。

第二,仿真不是替代实验,而是为实验“指路”。损耗仿真算得再准,也要知道它的适用范围——铁耗模型、谐波附加损耗、温度对电阻的影响,这些都是简化处理。工程上做损耗实测时,空载损耗受电压谐波影响、负载损耗受绕组温度影响,仿真结果和实测对不上是常态,关键是要能解释差异来源,而不是一味调参数让结果“看起来对”。能解释差异,才是真本事。

这套模型从参数到线路再到验证跑通,我前后花了大约两个工作日,中间大部分时间耗在铁耗系数标定和测量窗口选择上。如果你照着我上面的参数和步骤来,应该能在半天内跑出第一版结果。之后不管是做论文验证,还是工程方案对比,这套模型都能当基底来用。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询