ARTICLE · INTELLIGENCE

战地情报 · 详情页

来自尧图项目组的一线实战观察与深度解析

基于Simulink的换流变压器空负载损耗仿真建模与参数换算

基于Simulink的换流变压器空负载损耗仿真建模与参数换算 做直流输电工程或换流站设备评估时我经常要在只有铭牌参数、没有设计图纸的情况下回答一个问题这台换流变压器到底要承受多少损耗。空载损耗和负载损耗是绕不开的两个基础指标——前者反映铁心材料与工艺水平后者决定绕组载流能力和散热设计。基于Simulink搭建换流变压器空负载损耗仿真线路本质上不是在做电磁场精确求解而是用铭牌参数反推变压器等效电路在仿真环境里把出厂试验重做一遍。这篇文章记录我常用的搭建思路、参数换算和调试要点适合正在用Simulink做电力设备仿真的工程师也适合刚接触变压器模型的研究生照着操作。1. 换流变压器的损耗从哪里来空载与负载的物理本质1.1 空载损耗铁心里的磁滞与涡流空载损耗在行业内通常直接叫铁耗。变压器空载运行时一次侧施加额定电压二次侧开路此时变压器从电网吸收的功率几乎全部消耗在铁心里。铁耗由两部分组成磁滞损耗和涡流损耗。磁滞损耗来自铁心材料在交变磁化过程中磁畴反复翻转所消耗的能量它和磁通密度最大值、频率以及铁心材料的磁滞回线面积直接相关。涡流损耗则是交变磁通在铁心硅钢片内部感应出的环流所产生的焦耳热为了抑制这部分损耗硅钢片被轧得很薄、片间做绝缘处理。工程计算公式里这两部分通常合并成经典表达式PFe Kh × f × Bm^a Ke × f^2 × Bm^2其中Kh和Ke跟材料有关指数a通常在1.6到2.0之间Bm是铁心最大工作磁密。换流变压器的铁心工作磁密往往比普通变压器取得略低一点因为要留出谐波磁通和直流偏磁的裕度。这些微观机理在Simulink仿真里并不会被直接计算。我们搭模型时用的是一个等效铁耗电阻Rc并联在励磁支路上用电阻消耗的有功功率来等效铁耗。这种做法牺牲了对铁心材料物理过程的还原但结果是可用的——因为出厂试验给我们的就是额定电压下的总铁耗P0我们只需要保证仿真中励磁支路在该电压下消耗的功率等于P0就达到了工程精度的要求。1.2 负载损耗绕组电阻加杂散涡流负载损耗又叫铜耗严格说包含两部分绕组导线电阻产生的直流电阻损耗以及漏磁场在绕组、夹件、油箱壁等金属结构件中感应出的杂散涡流损耗。铜耗遵循焦耳定律和电流的平方成正比所以负载损耗随负载率变化非常敏感。杂散损耗则跟漏磁场的分布有很大关系换流变压器因为绕组结构特殊、漏磁较大杂散损耗占比通常比普通双绕组变压器略高。测量负载损耗的标准方法是短路试验把二次侧三相短路在一次侧施加一个较低的电压使绕组电流达到额定值此时变压器输入的有功功率就是额定电流下的负载损耗Pk。有人不理解为什么短路试验能测出负载损耗——因为变压器带负载时负载电流在绕组里产生的损耗与短路试验中额定电流产生的损耗在数值上是等价的只是运行工况不同。1.3 铭牌参数与试验测量仿真要复现的目标搭仿真之前先要清楚手里有哪些参数。一台典型的换流变压器铭牌上会有额定容量、额定电压、频率、联结组别、空载电流百分比、空载损耗、负载损耗和短路阻抗百分比。这些数据就是仿真的输入和校验基准。我用一组示例参数贯穿全文方便读者对照换算参数数值额定容量240 MVA一次侧额定电压330 kV二次侧额定电压110 kV额定频率50 Hz联结组别Yd11短路阻抗uk%18%空载电流i0%0.25%空载损耗P0180 kW负载损耗Pk75℃720 kW一次侧额定电流I1n420 A240 MVA / √3 / 330 kV仿真目标非常明确空载试验模型加额定电压测量得到的有功功率应该约等于180 kW短路试验模型让一次侧电流达到420 A测量得到的有功功率应该约等于720 kW未经温度折算时以设置温度为准。如果偏差在1%以内模型参数就算标定到位。2. 仿真路线怎么选搭等效电路比硬算电磁场更聪明2.1 为什么选Simulink/Simscape而不是有限元换流变压器的损耗理论上可以用有限元软件做二维或三维涡流场计算但那要求完整的铁心结构、绕组尺寸、材料B-H曲线甚至叠片系数这些数据通常只有变压器厂自己掌握。对于工程评估方或者系统研究方来说能拿到的只有铭牌参数和试验报告有限元建模成本太高、周期太长。Simulink加Simscape Electrical原SimPowerSystems的优势在于它不需要三维几何模型只要把变压器T型等效电路中的几个电阻和电感参数设置正确就能在毫秒级仿真时间内复现损耗。对于系统级研究比如换流站损耗统计、温升估算、效率曲线绘制这种精度完全够用。2.2 变压器T型等效电路与模块选型对比变压器线性等效电路的核心是T型结构一次侧漏阻抗Z1、励磁支路Zm、二次侧折算漏阻抗Z2串联。励磁支路由铁耗电阻Rc和磁化电感Lm并联组成Z1和Z2分别由绕组电阻和漏电感构成。在Simulink里实现这个电路有两类路线路线实现方式优点缺点ASimscape Electrical的Three-Phase Transformer (Two Windings)模块参数页直接填额定容量、电压、短路阻抗百分比、空载电流百分比搭起来快自带三相连接内部结构黑盒加谐波源或直流偏置时不如自己搭回路灵活BSimulink基础模块加Three-Phase Series RLC Branch自己搭T型等效每个元件都可独立控制方便注入谐波、设置非线性电感、观察支路内部电压电流需要自己处理三相回路和理想变压器变比连接出错概率高我的习惯是第一次做标准空载和负载损耗验证时用路线A速度快参数也好设置后面要做谐波影响或直流偏磁分析时改成路线B因为这时需要在特定支路上注入谐波电流或者串入直流源黑盒模块反而碍事。2.3 参数换算从铭牌到标幺值的完整计算不管用哪条路线都要先把铭牌参数换算成等效电路参数。以我给出的示例变压器为例逐步算一遍基准阻抗Zb U1n^2 / Sn (330 kV)^2 / 240 MVA 453.75 Ω短路阻抗和负载损耗Zk* uk% / 100 0.18即Zk 0.18 × 453.75 81.68 ΩRk* Pk / Sn 0.72 MW / 240 MVA 0.003即Rk 0.003 × 453.75 1.36 ΩXk* √(Zk*^2 - Rk*^2) √(0.18^2 - 0.003^2) ≈ 0.18即Xk ≈ 81.67 Ω一次侧、二次侧平均分配R1 R2 0.68 ΩX1 X2 40.84 Ω折算到一次侧励磁支路参数Rc* Sn / P0 240 / 0.18 1333.3即Rc 1333.3 × 453.75 ≈ 605 kΩXm* ≈ 1 / (i0% / 100) 1 / 0.0025 400即Xm 400 × 453.75 ≈ 181.5 kΩLm Xm / (2πf) 181500 / 314.16 ≈ 577.7 H这里有个容易算错的地方Rc*是Sn/P0而不是P0/Sn因为电阻标幺值是实际电阻除以基准阻抗而基准阻抗本身就是由U1n^2/Sn定义的代入一整理就得出Sn/P0。我第一次换算时把这个比值写反了仿真空载损耗差了上千倍后来对着功率读数才反应过来。如果使用Simscape的变压器模块可以在参数页直接填标幺值也可以在SI单位页填欧姆和亨利。我建议全程用标幺值因为出厂试验报告里给出的短路阻抗和空载电流本身就是百分数直接换算最不容易出错而且后续要改容量或电压等级时只需要重算基准值。3. 空载损耗仿真线路搭建励磁支路与铁耗电阻的落位3.1 线路结构电压源、测量环节和开路副边用路线A搭建空载损耗仿真时模型配置如下Three-Phase Programmable Voltage Source线电压有效值330 kV频率50 HzThree-Phase V-I Measurement串在一次侧电源和变压器之间用于测量电压电流Three-Phase Transformer (Two Windings)参数按上一节的换算结果填入Three-Phase Active Reactive Power接在测量模块后面读取有功功率Scope观测三相电压、电流和有功功率波形变压器二次侧什么负载都不接直接悬空。空载试验的本质就是副边开路一次侧加额定电压。这里要注意Simscape的电压源默认填的是相电压峰值。330 kV是线电压有效值换算成相电压峰值是330 kV × √2 / √3 ≈ 269.4 kV。很多新手直接填330 kV结果空载损耗是额定值的3倍——因为电压加高了。这个细节在Simulink文档里写得不显眼但在实际使用中几乎每个人都会踩一次。3.2 一次侧电压源的幅值与相位设置避开合闸涌流空载损耗仿真看似简单直接给电压源加额定电压就行但如果参数设置不对波形前几个周期会出现很大的暂态电流导致有功功率读数跳动。问题出在合闸相位。如果在A相电压过零点时投入电源变压器铁心磁通会面临最严重的暂态偏置相当于出现一个直流磁通分量铁心在一段时间内进入饱和区励磁电流可能冲到稳态值的几十倍。这个现象叫合闸涌流和真实变压器空载合闸时一模一样。解决办法有两个第一把Programmable Voltage Source的初始相位设置在电压峰值处也就是Phase A设为90°这样合闸瞬间磁通正好从零开始建立没有直流偏置分量涌流基本可以忽略。第二仿真时间设置得长一些等暂态衰减完再读数。线性电感模型的暂态衰减其实很快但如果是后面加了非线性饱和电感涌流衰减可能需要几百毫秒甚至几秒那时就必须用第一种方法或者加合闸控制逻辑。我通常两种方法同时用初始相位90°加仿真时间0.3秒取最后0.1秒的功率平均值。3.3 空载损耗读数方法与稳定性判断有功功率读数直接看Three-Phase Active Reactive Power模块的输出就行但不要只读某一个个瞬间的值。Simulink的功率测量模块对三相瞬时功率做了低通滤波和平均值处理输出已经接近稳态平均功率不过受数值求解器误差影响波形上会有微小纹波。我的做法是把功率信号接到Scope里再并联一个Mean模块或者直接右键信号选择Mean统计方式设置基波周期为0.02秒对整个仿真时长求平均。更稳妥的是只对最后两个周波求平均因为暂态过程已经衰减完毕。按照示例参数跑完空载损耗读数应该在175 kW到185 kW之间。如果严重偏大先检查电压源的幅值是否按相电压峰值填写如果严重偏小检查励磁支路的Rc是否填入正确。这里有个调试技巧把Scope里同时显示一次侧电压有效值和电流有效值如果电压正确、电流远小于0.25%额定电流那就是Rc或者Lm填反了。4. 负载损耗仿真线路搭建副边短路试验的仿真再现4.1 线路结构副边三相短路与一次侧降压加压负载损耗仿真复现的是短路试验模型配置和空载试验很像只有两处改变一是变压器副边从悬空改为三相短路二是一次侧电压不再是额定电压而是大幅降低后的电压。副边短路在Simscape里可以用Three-Phase Fault模块实现设置三相同时接地短路直接用一个三相Y接的阻抗为零的RLC支路也可以但Three-Phase Fault模块更直观还能设置短路发生时间。一次侧电压源的幅值先按估算值设置Ushort ≈ U1n × uk% 330 kV × 18% 59.4 kV线电压有效值。换算成相电压峰值就是59.4 × √2 / √3 ≈ 48.5 kV。注意短路试验中一次侧电流要达到额定值420 A但二次侧绕组实际流过的电流非常大因为短路状态下变压器相当于一个带漏阻抗的纯电抗回路一次侧施加很小的电压就能推动额定电流。Yd11联结的副边三角形绕组内部环流会让三相短路的电流分布看起来与Y接不同但这个不影响损耗测量功率测量只需要关注一次侧。4.2 电压怎么找阻抗电压百分数的粗调与线性缩放短路试验的电压需要精确调整到让电流等于额定值。理论值是额定电压的18%但因为绕组电阻和漏感的比例、温度换算等原因第一次粗调后电流可能不等于420 A。这里有一个非常实用的方法短路试验回路近似线性电流和电压成正比。如果第一次仿真发现电流只有380 A那就把电压按比例放大420/380倍再跑一次第二次基本就能准确落在额定电流附近。例如第一次设置59.4 kV读得电流380 A则调整后电压为59.4 × (420 / 380) ≈ 65.6 kV。如果读得电流450 A则下调到59.4 × (420 / 450) ≈ 55.4 kV。这种线性缩放调整一般只需要一到两次迭代就能收敛到额定电流比手动试凑效率高很多。读电流时同样要等稳态取最后一个周波的均方根值。当一次侧电流稳定在420 A左右时Three-Phase Active Reactive Power读出的有功功率就是负载损耗仿真值。由于示例参数里绕组电阻是按标幺值Rk* 0.003填入的而这个标幺值对应的正是75℃下的负载损耗所以仿真读数应接近720 kW。4.3 温度折算与75℃参考温度的处理工程上负载损耗必须折算到参考温度最常用的是75℃。铜绕组的温度折算公式是Pk75 Pk_meas × (235 75) / (235 T_meas)其中T_meas是试验时的绕组温度235是铜绕组的温度系数铝绕组用228。仿真模型里的绕组电阻不会随温度自动变化所以要把这一步放到仿真外面处理。我的习惯是先按75℃的目标损耗反推电阻标幺值填入模型这样仿真结果直接就是75℃下的损耗不用再折算。如果研究负载率变化对温升的影响再把电阻设置成对应温度下的值分别跑几个工况。还有一个容易忽略的点杂散损耗中对温度的响应和直流电阻不同。铜绕组的直流电阻随温度升高而增大但涡流杂散损耗随温度升高略有下降因为金属结构件的电阻率升高导致涡流减小。完整的温度折算应该把直流电阻损耗和杂散损耗分开处理分别按正、负温度系数折算后相加。一般工程评估用简化折算就能满足精度但做精细化损耗分离时要注意这个方向。5. 仿真结果校验与数值收敛的坑别让初始条件带偏读数5.1 为什么第一周期读数不能信不管是空载还是负载损耗仿真第一个周期甚至前几个周期的功率波形都不可信。原因有两个层面。第一是数值求解器的初始化。Simscape电网络在仿真启动时需要确定状态变量的初始值变压器电感电流和电容电压的初始条件如果没有显式指定求解器默认从零开始。这本身没有问题但电压源从零相位开始施加时电感的暂态响应会导致电流波形带有衰减的直流分量。第二是代数环问题。如果模型里存在瞬时反馈回路比如功率测量模块的输出参与控制、控制输出又反过来改变电压源求解器在首步迭代中可能还没收敛造成开头几个采样点的数值毛刺。所以无论是读数还是波形观察都要跳过前两个周波。我习惯把仿真时间设到0.3秒以上用Mean模块统计最后0.1秒的平均功率。如果波形显示有功功率还在缓慢爬升说明暂态还没衰减完加大仿真时间。5.2 浮空回路与代数环求解器报错的常见原因Simscape的变压器模型在副边完全悬空时理论上是可以直接仿真的但有时会报出Simulink cannot solve the algebraic loop或者提示某节点没有参考电位。这不是模型原理错而是数值求解器对完全浮空的电位没有参考点。我遇到这个问题时的标准解法在变压器原边或副边的中性点接一个阻值很大的电阻比如1 MΩ到地给网络提供一个高阻参考路径。这个电阻对工频运行的影响微乎其微却能稳定求解器的数值状态。另一个更隐蔽的问题是Three-Phase V-I Measurement的接法。电压测量端必须并联在被测支路两端电流测量端必须串联在被测支路中接反了不仅读数错误还可能破坏电流回路的连续性导致求解器报错。这个在Simscape里是模块画法决定的拖动模块时要注意端口的标注。5.3 误差项分析线性电感、谐波与温度的隐藏偏差即使模型参数完全正确仿真读数和铭牌值之间也会存在微小的系统性偏差我把它拆成几个来源线性磁化电感假设真实变压器的励磁电感随电压变化饱和点附近电感量下降空载电流会略大于线性模型的计算值。但由于空载电流本身很小0.25%对铁耗的读数影响不大。谐波含量Simulink的电压源默认是纯正弦波而实际电网和换流站母线电压总有谐波。谐波电压作用在变压器上会额外增加铁耗和绕组涡流损耗这会导致实际运行损耗高于仿真值。参考温度不一致负载损耗随温度变化如果仿真电阻设定温度和铭牌参考温度不同读数直接对比必然有偏差。把这些误差控制在1%以内是合理的。如果偏差超过2%优先检查参数换算而不是怀疑模块精度。我见过不少人花大量时间调求解器精度结果问题出在Rc和Lm填反了——这种低级错误在示波器上看功率数值时其实很显眼只要对照铭牌值算一下就发现了。6. 扩展方向谐波损耗、直流偏磁与模型复用6.1 特征谐波电流下的损耗修正程序化电压源的妙用换流变压器和普通变压器最大的区别在于它流过大量特征谐波电流。六脉动换流器产生的特征谐波次数是6k±1也就是5、7、11、13次等。谐波电流会造成附加损耗这部分损耗在实际运行中不可忽略。在Simulink里做谐波损耗扩展只需要把Three-Phase Programmable Voltage Source的Harmonics选项卡打开设置各次谐波的幅值比例和相位。例如5次谐波设为基波幅值的15%7次设为10%11次设为6%更精确的比例可以根据直流系统典型潮流反推或用电能质量测试仪实测。谐波损耗的分析思路有两种一种直接看总损耗变化谐波叠加后总有功会上升差值就是谐波附加损耗另一种用FFT把电压电流分解对每次谐波单独计算损耗再求和。前者简单适合工程估算后者能看清各次谐波的贡献适合写分析报告。绕组谐波损耗有个近似规律高频电流的集肤效应使有效电阻随频率近似按f^1.5到f^2增长。所以即使11次谐波电流幅值只有基波的6%其产生的绕组损耗却可能达到基波损耗的几个百分点。6.2 直流偏磁与铁心饱和从线性模型到非线性电感的升级换流变压器在单极运行或地磁感应电流作用下绕组中会流过直流分量。直流偏磁会使铁心工作点偏移严重时半周波饱和励磁电流畸变、损耗增大、噪声加剧这是换流变压器区别于普通变压器的重要故障模式。线性Lm模型无法仿真这种现象因为线性电感的电感量恒定不会饱和。要升级就得把励磁支路的线性电感换成非线性电感Simscape里有Nonlinear Inductor模块可以输入B-H曲线或分段线性化后的磁通-电流特性。具体做法把磁化支路从变压器模块中拆出来用非线性电感替代同时保留Rc铁耗电阻。铁耗电阻此时可以设为常数也可以引入频率和磁密修正看研究深度。注入直流分量的方法是在一次侧串联一个可调直流电压源或者用电流源在变压器带载状态下叠加一个直流偏置电流。非线性电感的参数整定比较麻烦需要知道铁心截面积、磁路长度和B-H曲线。如果没有厂家数据可以用V-I特性拟合法根据空载试验的电压电流数据反推磁化曲线。我通常用饱和点与膝点两个特征点做分段线性化仿真结果已经能定性反映偏磁对励磁电流波形和空载损耗的影响。6.3 模型复用与联合仿真C代码生成、实时仿真与设计校核损耗模型标定好以后它不只是算两个数字还可以作为一个可复用的虚拟变压器嵌到更大的系统里。比如在Simulink里搭一个完整的高压直流输电整流逆变回路把换流变压器模型放在阀厅变压器位置上就能实时观察不同运行工况下变压器的损耗分布配合同步读数的功率测量模块实际上就是一个虚拟功率表。Simulink模型还能通过Embedded Coder生成C代码部署到实时仿真机或控制器原型里。损耗仿真模型生成C代码后可以作为软件在环测试的参考模型也可以集成到保护装置的闭锁逻辑验证中。我自己在工程中把这类模型用于培训仿真平台——学员操作换流站主接线画面上的变压器损耗实时刷新这对理解运行工况很有帮助。模型复用的另一条路是和外部程序联合仿真。Simulink可以被MATLAB脚本批量调用也可以和Carsim、Amesim等专业软件做联合仿真。对换流变压器损耗模型来说联合仿真的典型场景是和多物理场温度模型对接损耗仿真结果作为热源传递给温度场模型温度场模型再把绕组温度反馈回来修正电阻值形成一个损耗-温度闭环。Simulink的S函数或FMI/FMU接口都能实现这类耦合虽然搭建成本高一些但对长期运行的损耗评估非常有价值。我个人始终觉得用Simulink搭换流变压器损耗模型的价值不在于算出那两个额定工况数据而在于留下一个能反复修改、能嵌入系统、能应对谐波和偏磁等非标准工况的分析工具。铭牌参数会变、运行工况会变但T型等效电路和试验复现这套方法论是稳定的。把这条线路摸熟以后处理其他变压器或者电抗器的损耗问题思路基本都是相通的。
RELATED READING

延伸阅读

更多一线实战笔记与深度复盘,助您持续精进