作为一个常年和电力电子仿真打交道的人我最近刚好把换流变压器的空负载损耗仿真模型完整搭了一遍从参数计算到Simulink建模再到结果分析整个过程踩了不少坑。这篇文章我把当时的思路、参数计算、建模步骤和排查过程整理出来希望能给正在做直流输电或变压器仿真的朋友一些参考。先说清楚这个东西到底解决什么问题。换流变压器和普通电力变压器不一样它在直流输电系统里承担着换流阀与交流系统之间的电压变换和电气隔离作用运行工况非常复杂阀侧电压电流含有大量谐波直流偏磁风险高绝缘结构特殊。损耗仿真如果做得糙后面热设计、效率评估、经济性分析全部要返工。所以损耗仿真这件事看起来只是算两个数空载损耗P0和负载损耗Pk实际上一旦做细牵扯到谐波、饱和、频率特性工程量并不小。本文适合刚接触Simulink电力系统仿真的学生也适合已经在做直流输电相关研究的工程师。我会从模型思路、参数计算、仿真搭建到问题排查完整走一遍过程中所有关键参数都给出计算过程方便你直接套用。1. 模型整体思路与方案选型1.1 为什么用Simulink而不是自编程序换流变压器损耗仿真听起来可以写个解析公式算一算但实际工程场景里没这么简单。空载损耗取决于铁芯的磁化特性和铁耗特性负载损耗取决于绕组电阻和漏抗这两部分在谐波环境下还要考虑频率的影响。如果单纯用公式估算谐波电流作用下的附加损耗很难算准。而Simulink/Simscape的优势在于它可以构建完整的电路拓扑用电磁暂态的视角去看变压器的工作状态谐波的影响是自动被包含在波形计算里的。另外Simulink生态里有现成的变压器模型。Simscape Electrical库中的Specialized Power SystemsSPS提供了单相和三相变压器模型包括线性变压器、饱和变压器等这些模型内置了磁路饱和、铁耗计算能力对于损耗仿真来说已经非常够用。比起自己写差分方程求解磁链电流关系用现成模型至少能节省一周时间。我的总体思路是用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库推荐用离散求解器尤其是模型中带有电力电子开关器件时。离散仿真用固定的采样步长避免了连续求解器在开关切换时刻的迭代收敛问题速度也快很多。虽然我现在做的是纯变压器空负载仿真没有电力电子开关但我还是建议从第一天就用离散模式因为换流变压器的应用场景终究是要接换流阀的早点定下离散仿真的框架后面不用重新搭模型。采样时间的选择也讲究。我的经验是先粗后细先用Ts50μ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”离散不需要像连续系统那样选具体算法。唯一需要认真设的就是采样时间。我实测过同一个空载损耗模型Ts50μs时仿真结果和Ts5μ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。我第一版用Ts50μ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 空载损耗与负载损耗仿真操作实录空载损耗仿真副方三相RLC负载设为开路等效于不接负载。原方接额定电压220kV的50Hz三相电源。仿真时长设0.5s求解器用离散模式Ts50μs。仿真结束后查看原方电压电流波形。正常情况空载电流幅值很小额定电流0.3%左右波形是尖顶波不是正弦波。这是铁芯饱和的最直观证据。读取“Active Reactive Power”模块的有功功率用0.2s到0.5s之间的平均数值作为空载损耗。按我们的参数空载电流有效值应该在0.78A左右262.4×0.003实际上因为励磁电流非线性可能稍微偏高一点只要在0.5~1.2A范围内都合理。负载损耗仿真副方接额定负载负载功率设为100MVA功率因数可以取0.8或0.9。原方仍然接220kV额定电压仿真时间0.5s。读取原方有功功率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分钟还没结束。原因求解器设置不合理用了ode45去解刚性系统带电力电子或不连续环节。模型里有无穷大阻抗/导纳的悬空节点导致代数环或数值振荡。变压器模型饱和特性给的饱和点太密导致插值计算慢。解决办法改用ode23tb或ode15s或者直接切离散模式。把Powergui切换为离散模式采样时间从50μs起步如果还慢就调到100μs看趋势。检查有没有悬空端口尤其是变压器中性点、三角绕组不接地的那一侧Simulink里悬空节点会引发数值问题可以接地Ground模块或者加高阻接地。4.2 仿真直接发散波形飞到天际现象电压电流瞬间变成1e30这种离谱数值波形直接拉满纵轴。原因变压器参数单位填错常见的是把标幺值当国际单位填或反之。负载或电源初始相位设置冲突。变压器励磁电感太小导致励磁电流极大。解决办法先检查变压器模块的所有参数是不是按铭牌或按上面计算过程填的特别是短路阻抗百分比填成0.15还是15区别巨大。把仿真初始时刻的电压源相位设为0负载初始电压不要给太大。把变压器饱和特性曲线砍掉先按线性模型跑通再加非线性。用“Powergui → Tools → Compute Initial States”让软件自动计算初值这是一个救命功能。4.3 空载损耗仿真结果比设计值差一半以上现象标称80kW的空载损耗仿真出来只有30kW或者变成150kW。原因铁耗系数Steinmetz系数没填或者填的是默认值和实际硅钢片损耗特性不匹配。空载电流取的百分数不对。测量窗口覆盖了暂态启动过程平均值被拉偏。解决办法在变压器模块里仔细填铁耗三个系数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.1s5个周期而不是0.13s这种非整数周期。如果信号里有大量开关谐波先经过一个低通滤波器或者把采样率提高避免混叠。4.5 换流变压器特有的阀侧谐波问题现象阀侧电流波形明显畸变谐波含量远超预期损耗数值不稳定。原因阀侧连接的是整流器电流含大量特征谐波6k±1次12脉波系统主要是11、13次这些谐波在绕组里产生附加损耗涡流损耗模型如果不考虑频率特性损耗会偏低。换流变压器的漏抗和谐波阻抗特性与普通变压器不同简单模型无法准确反映。解决办法如果研究的是谐波损耗可以考虑用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%明显不合理还分析得头头是道。所以跑完仿真一定要用理论值做一次数量级校验这是从业者的基本素养。第二仿真不是替代实验而是为实验“指路”。损耗仿真算得再准也要知道它的适用范围——铁耗模型、谐波附加损耗、温度对电阻的影响这些都是简化处理。工程上做损耗实测时空载损耗受电压谐波影响、负载损耗受绕组温度影响仿真结果和实测对不上是常态关键是要能解释差异来源而不是一味调参数让结果“看起来对”。能解释差异才是真本事。这套模型从参数到线路再到验证跑通我前后花了大约两个工作日中间大部分时间耗在铁耗系数标定和测量窗口选择上。如果你照着我上面的参数和步骤来应该能在半天内跑出第一版结果。之后不管是做论文验证还是工程方案对比这套模型都能当基底来用。
换流变压器空负载损耗仿真：Simulink建模与参数计算全流程
作为一个常年和电力电子仿真打交道的人我最近刚好把换流变压器的空负载损耗仿真模型完整搭了一遍从参数计算到Simulink建模再到结果分析整个过程踩了不少坑。这篇文章我把当时的思路、参数计算、建模步骤和排查过程整理出来希望能给正在做直流输电或变压器仿真的朋友一些参考。先说清楚这个东西到底解决什么问题。换流变压器和普通电力变压器不一样它在直流输电系统里承担着换流阀与交流系统之间的电压变换和电气隔离作用运行工况非常复杂阀侧电压电流含有大量谐波直流偏磁风险高绝缘结构特殊。损耗仿真如果做得糙后面热设计、效率评估、经济性分析全部要返工。所以损耗仿真这件事看起来只是算两个数空载损耗P0和负载损耗Pk实际上一旦做细牵扯到谐波、饱和、频率特性工程量并不小。本文适合刚接触Simulink电力系统仿真的学生也适合已经在做直流输电相关研究的工程师。我会从模型思路、参数计算、仿真搭建到问题排查完整走一遍过程中所有关键参数都给出计算过程方便你直接套用。1. 模型整体思路与方案选型1.1 为什么用Simulink而不是自编程序换流变压器损耗仿真听起来可以写个解析公式算一算但实际工程场景里没这么简单。空载损耗取决于铁芯的磁化特性和铁耗特性负载损耗取决于绕组电阻和漏抗这两部分在谐波环境下还要考虑频率的影响。如果单纯用公式估算谐波电流作用下的附加损耗很难算准。而Simulink/Simscape的优势在于它可以构建完整的电路拓扑用电磁暂态的视角去看变压器的工作状态谐波的影响是自动被包含在波形计算里的。另外Simulink生态里有现成的变压器模型。Simscape Electrical库中的Specialized Power SystemsSPS提供了单相和三相变压器模型包括线性变压器、饱和变压器等这些模型内置了磁路饱和、铁耗计算能力对于损耗仿真来说已经非常够用。比起自己写差分方程求解磁链电流关系用现成模型至少能节省一周时间。我的总体思路是用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库推荐用离散求解器尤其是模型中带有电力电子开关器件时。离散仿真用固定的采样步长避免了连续求解器在开关切换时刻的迭代收敛问题速度也快很多。虽然我现在做的是纯变压器空负载仿真没有电力电子开关但我还是建议从第一天就用离散模式因为换流变压器的应用场景终究是要接换流阀的早点定下离散仿真的框架后面不用重新搭模型。采样时间的选择也讲究。我的经验是先粗后细先用Ts50μ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”离散不需要像连续系统那样选具体算法。唯一需要认真设的就是采样时间。我实测过同一个空载损耗模型Ts50μs时仿真结果和Ts5μ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。我第一版用Ts50μ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 空载损耗与负载损耗仿真操作实录空载损耗仿真副方三相RLC负载设为开路等效于不接负载。原方接额定电压220kV的50Hz三相电源。仿真时长设0.5s求解器用离散模式Ts50μs。仿真结束后查看原方电压电流波形。正常情况空载电流幅值很小额定电流0.3%左右波形是尖顶波不是正弦波。这是铁芯饱和的最直观证据。读取“Active Reactive Power”模块的有功功率用0.2s到0.5s之间的平均数值作为空载损耗。按我们的参数空载电流有效值应该在0.78A左右262.4×0.003实际上因为励磁电流非线性可能稍微偏高一点只要在0.5~1.2A范围内都合理。负载损耗仿真副方接额定负载负载功率设为100MVA功率因数可以取0.8或0.9。原方仍然接220kV额定电压仿真时间0.5s。读取原方有功功率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分钟还没结束。原因求解器设置不合理用了ode45去解刚性系统带电力电子或不连续环节。模型里有无穷大阻抗/导纳的悬空节点导致代数环或数值振荡。变压器模型饱和特性给的饱和点太密导致插值计算慢。解决办法改用ode23tb或ode15s或者直接切离散模式。把Powergui切换为离散模式采样时间从50μs起步如果还慢就调到100μs看趋势。检查有没有悬空端口尤其是变压器中性点、三角绕组不接地的那一侧Simulink里悬空节点会引发数值问题可以接地Ground模块或者加高阻接地。4.2 仿真直接发散波形飞到天际现象电压电流瞬间变成1e30这种离谱数值波形直接拉满纵轴。原因变压器参数单位填错常见的是把标幺值当国际单位填或反之。负载或电源初始相位设置冲突。变压器励磁电感太小导致励磁电流极大。解决办法先检查变压器模块的所有参数是不是按铭牌或按上面计算过程填的特别是短路阻抗百分比填成0.15还是15区别巨大。把仿真初始时刻的电压源相位设为0负载初始电压不要给太大。把变压器饱和特性曲线砍掉先按线性模型跑通再加非线性。用“Powergui → Tools → Compute Initial States”让软件自动计算初值这是一个救命功能。4.3 空载损耗仿真结果比设计值差一半以上现象标称80kW的空载损耗仿真出来只有30kW或者变成150kW。原因铁耗系数Steinmetz系数没填或者填的是默认值和实际硅钢片损耗特性不匹配。空载电流取的百分数不对。测量窗口覆盖了暂态启动过程平均值被拉偏。解决办法在变压器模块里仔细填铁耗三个系数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.1s5个周期而不是0.13s这种非整数周期。如果信号里有大量开关谐波先经过一个低通滤波器或者把采样率提高避免混叠。4.5 换流变压器特有的阀侧谐波问题现象阀侧电流波形明显畸变谐波含量远超预期损耗数值不稳定。原因阀侧连接的是整流器电流含大量特征谐波6k±1次12脉波系统主要是11、13次这些谐波在绕组里产生附加损耗涡流损耗模型如果不考虑频率特性损耗会偏低。换流变压器的漏抗和谐波阻抗特性与普通变压器不同简单模型无法准确反映。解决办法如果研究的是谐波损耗可以考虑用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%明显不合理还分析得头头是道。所以跑完仿真一定要用理论值做一次数量级校验这是从业者的基本素养。第二仿真不是替代实验而是为实验“指路”。损耗仿真算得再准也要知道它的适用范围——铁耗模型、谐波附加损耗、温度对电阻的影响这些都是简化处理。工程上做损耗实测时空载损耗受电压谐波影响、负载损耗受绕组温度影响仿真结果和实测对不上是常态关键是要能解释差异来源而不是一味调参数让结果“看起来对”。能解释差异才是真本事。这套模型从参数到线路再到验证跑通我前后花了大约两个工作日中间大部分时间耗在铁耗系数标定和测量窗口选择上。如果你照着我上面的参数和步骤来应该能在半天内跑出第一版结果。之后不管是做论文验证还是工程方案对比这套模型都能当基底来用。
