1. 先搞清楚：你的模型到底在“保护”什么
做变压器纵联差动保护仿真的人，十个里有八个踩过同一个坑——模型跑起来了，波形看着也挺像那么回事，但一到区内故障就拒动，区外故障反而误动，调来调去不知道问题出在哪儿。我带过几个做继电保护方向的同学做这类仿真，几乎每个人第一版模型都会在这两个极端之间反复横跳。问题不在于MATLAB用得不熟，而在于对差动保护底层逻辑的理解有断层——仿真只是把理解偏差放大了而已。
变压器纵联差动保护的核心思路其实不复杂：比较变压器各侧电流的相量和，正常运行时（包括区外故障）差流接近零，区内故障时差流急剧增大，保护动作。听起来简单，但真正落到仿真模型里，从CT建模、相位补偿、幅值归算到制动特性整定，每一步都有能让你的模型“看起来对但实际错”的细节。这些细节在教科书上往往一笔带过，在实际工程中却是调试报告里反复出现的老面孔。
这篇内容面向的是正在用MATLAB/Simulink搭建变压器差动保护仿真模型、但被误动和拒动问题卡住的人。无论你是刚接触继电保护仿真的新手，还是已经能跑通模型但对结果不够有把握的老手，下面拆解的这些环节都值得对照检查一遍。我会把每个关键步骤背后的“为什么”讲清楚，把参数选择的计算过程摆出来，把踩过的坑和排查思路整理成可以直接抄作业的清单。全文基于常见的工程实践和仿真经验展开，具体参数需要根据你手头的变压器铭牌和CT配置做调整。
2. 整体方案设计：为什么你的仿真框架从一开始就可能偏了
2.1 差动保护的逻辑骨架与仿真映射
纵联差动保护在仿真里的实现，本质上是对一个电流比较方程的离散化求解。最朴素的表达是：
Id = |I1 + I2 + ... + In|
其中各侧电流需要先归算到同一侧（通常归算到高压侧或额定容量侧），再按变压器接线组别进行相位校正。保护判据通常写成：
Id > Icd（差动速断）或Id > Iop0 + K * Ir（比率制动）
这里Id是差流，Ir是制动电流，Iop0是启动门槛，K是制动斜率。仿真模型要做的，就是把一次侧电流经过CT变换、相位补偿、幅值归算后，代入这个判据，看保护是否动作。
问题在于，很多人在Simulink里搭模型时，直接把三相电流源接到一个“Sum”模块上取差值，然后接一个Relay模块做阈值判断，就认为差动保护搭好了。这种简化模型在理想条件下能跑出“正确”的波形，但一旦加入CT饱和、励磁涌流、区外故障穿越等实际因素，判据的边界条件就会完全失控。原因很简单：你跳过了CT传变特性、相位补偿矩阵、制动电流定义方式这三个决定保护行为的关键环节。
2.2 仿真平台选型：Simulink电力系统模块还是脚本化建模
在MATLAB生态里做变压器差动保护仿真，通常有两条路：一是用Simulink的Simscape Electrical（原SimPowerSystems）工具箱搭电路级模型，二是用MATLAB脚本做信号级仿真。两者各有适用场景，选错了会让后续调试事倍功半。
电路级模型的优势在于能自然引入CT饱和、铁芯非线性、断路器动作等物理过程，适合研究保护与一次系统的交互行为。缺点是仿真速度慢，参数调整不灵活，尤其是CT模型和变压器模型的参数耦合度高，改一个变比可能要重新整定一串参数。信号级模型则相反，你可以用传递函数或分段线性化模型替代CT，用相量或瞬时值信号直接驱动保护判据，调试效率高，但需要你自己保证信号的真实性。
我的建议是：如果你关注的是保护判据本身的正确性（比如比率制动特性、二次谐波闭锁逻辑），优先用信号级建模；如果你要研究CT饱和对保护的影响或涌流工况下的动作行为，必须用电路级模型。两者也可以混合使用——用电路级模型生成故障电流波形，导出到工作空间后用脚本做保护算法验证，这样兼顾了物理真实性和调试灵活性。
2.3 一次系统建模的简化边界在哪里
变压器模型的选择直接影响差流波形。Simulink里常用的三相变压器模型有线性模型、饱和模型和考虑铁芯剩磁的模型。很多人在做差动保护仿真时随手选一个默认的线性变压器，结果涌流工况下差流波形完全不对——因为线性模型根本不产生励磁涌流。
如果你要验证二次谐波闭锁逻辑，变压器模型必须包含饱和特性，并且要设置合理的铁芯饱和参数（饱和磁通密度、剩磁系数）。如果你只关心区内金属性短路和区外故障的稳态差流，线性模型够用，但要注意线性模型在故障暂态过程中不会产生真实的衰减直流分量，这会影响CT饱和的仿真结果。
另一个常被忽略的点是系统阻抗。变压器两侧的等效系统阻抗决定了故障电流的幅值和衰减时间常数。有些人在仿真里把电源设为理想电压源（内阻抗为零），结果区外故障时差流理论上为零，保护当然不动作——但这掩盖了实际系统中CT误差和暂态过程导致的虚假差流。合理的做法是根据短路容量反推系统阻抗，让仿真条件尽量接近真实工况。
3. 核心细节拆解：相位补偿、幅值归算与CT建模
3.1 相位补偿：Y/Δ接线带来的第一道坎
变压器常见的接线组别是Y/Δ-11，两侧电流存在30度相位差。如果不做补偿，正常运行时差流就会等于额定电流的0.5倍左右（因为两个幅值相等、相位差30度的相量相加，其模值为2Icos(75°)≈0.517I），保护直接误动。仿真里必须对两侧电流做相位校正。
常规做法有两种：一是将Y侧电流转换为Δ侧等效电流（Δ侧补偿），二是将Δ侧电流转换为Y侧等效电流（Y侧补偿）。在Simulink里实现时，前者需要对Y侧三相电流做如下运算：
IA_comp = (IA - IB) / √3
IB_comp = (IB - IC) / √3
IC_comp = (IC - IA) / √3
注意这里的√3系数——很多人漏掉它，导致补偿后幅值偏大√3倍，差流在正常运行时就不为零。另一种常见错误是补偿方向搞反了：应该把Y侧补偿到Δ侧，结果把Δ侧补偿到了Y侧，相位差不但没消除反而翻倍。
提示：在Simulink里实现相位补偿时，建议用“MATLAB Function”模块写补偿逻辑，而不是用Sum和Gain模块硬连。原因是补偿矩阵涉及三相之间的交叉运算，用图形化连线容易接错相序，而代码形式一目了然，改起来也方便。
补偿后的电流还需要做幅值归算。如果两侧CT变比不同，必须把二次电流归算到同一基准。归算公式为：
I1_base = I1_secondary * (nCT1 / nCT_base)
其中nCT_base通常取高压侧CT变比。这一步在仿真里容易被跳过，因为很多人直接在一次侧取电流信号做比较，忽略了CT变比的影响。如果你是在一次侧取信号，那确实不需要CT变比归算，但你就失去了仿真CT饱和的能力——这是一个需要提前想清楚的取舍。
3.2 制动电流的定义方式：选错公式，比率制动特性全歪
比率制动判据里，制动电流Ir的定义直接决定了保护的动作边界。常见的定义方式有三种：
| 制动电流定义 | 表达式 | 适用场景 | 特点 |
|---|---|---|---|
| 最大电流法 | Ir = max(|I1|, |I2|) | 两绕组变压器 | 制动量偏大，安全性高，灵敏度略低 |
| 算数和法 | Ir = (|I1| + |I2|) / 2 | 两绕组变压器 | 制动量适中，工程上常用 |
| 相量和法 | Ir = |I1 - I2| / 2 | 多绕组或特殊场景 | 制动量偏小，灵敏度高，安全性需校验 |
仿真里选哪种，取决于你模拟的保护装置采用的实际算法。很多国产保护装置采用算数和法或最大电流法，而一些文献里的仿真为了简化用相量和法。如果你选的制动电流定义和判据里的斜率K不匹配，动作区会明显偏移——比如用最大电流法整定的K值套到相量和法上，区外故障时制动量不足，保护就会误动。
我在早期做仿真时就犯过这个错：判据里写的是Ir = (|I1|+|I2|)/2，但制动斜率K是按最大电流法整定的0.5，结果区外故障穿越时差流还没到动作区，制动量先饱和了，保护拒动。后来把K调到0.3才匹配上。这个教训说明：制动电流定义和斜率整定必须来自同一套逻辑，不能混搭。
3.3 CT建模：理想CT和饱和CT的仿真差异有多大
CT模型是差动保护仿真里最容易被低估的环节。用理想CT（变比线性变换）做仿真，区外故障时差流恒为零，保护永远不会误动——这显然不符合实际。真实CT在故障暂态过程中可能饱和，导致二次电流波形畸变，两侧CT饱和程度不一致时就会产生虚假差流。
在Simulink里建饱和CT模型，通常有两种方式：一是用Simscape Electrical自带的饱和变压器模型配置成CT，二是用MATLAB Function写一个分段线性的励磁特性。前者物理意义清晰但参数多，后者灵活但需要自己定义饱和拐点和励磁电感。
关键参数包括：额定二次电流（通常5A或1A）、准确限值系数（ALF）、二次回路负担（包括引线电阻和继电器内阻）、剩磁系数。其中ALF决定了CT在多少倍额定电流下保持不饱和。比如ALF=20的CT，在一次电流达到20倍额定值时开始进入饱和区。仿真时如果ALF设得过大（比如50），CT几乎不会饱和，你就看不到虚假差流；设得过小（比如5），CT在正常故障电流下就饱和，差流大得离谱。
我的经验值是：做区外故障误动分析时，ALF取10~20比较接近实际中压系统的CT配置；做涌流工况分析时，CT饱和不是主要矛盾，可以先用理想CT跑通逻辑，再逐步加入饱和特性。另外，两侧CT的剩磁方向如果设为相同，饱和程度趋于一致，虚假差流反而小；设为相反或随机，差流会明显增大——这也是实际工程中CT剩磁导致误动的一个典型机理。
4. 实操过程：从零搭一个能跑对的差动保护仿真
4.1 一次系统参数计算与Simulink模型搭建
假设我们有一台三相双绕组变压器，铭牌参数为：额定容量SN = 50 MVA，电压比121/10.5 kV，接线组别Y/Δ-11，短路阻抗Uk = 10.5%。两侧CT变比分别为600/5和3000/5。
第一步，计算各侧额定电流：
高压侧：I1N = SN / (√3 * V1N) = 50e6 / (1.732 * 121e3) ≈ 238.6 A
低压侧：I2N = SN / (√3 * V2N) = 50e6 / (1.732 * 10.5e3) ≈ 2749.3 A
第二步，计算CT二次侧额定电流：
高压侧：I1N_sec = 238.6 / (600/5) ≈ 1.99 A
低压侧：I2N_sec = 2749.3 / (3000/5) ≈ 4.58 A
可以看到两侧二次额定电流不相等（1.99 A vs 4.58 A），必须做幅值归算。归算到高压侧CT二次基准：
低压侧归算值：I2N_sec_comp = 4.58 * (600/5) / (3000/5) = 4.58 * 0.2 = 0.916 A
等等，这里算出来0.916 A和高压侧的1.99 A还是不相等——说明我选的CT变比不匹配。正确的CT选型应该让两侧二次额定电流接近相等。重新选：高压侧CT取300/5（二次额定1.99300/600≈0.995A），低压侧CT取3000/5（二次额定4.58A），归算后4.58(300/5)/(3000/5)=0.458A，仍不匹配。实际上应该让低压侧CT变比更大或高压侧更小。这个计算过程本身就说明了一个常见错误：CT变比选型不当会导致归算后两侧电流不平衡，正常运行时就有差流。在仿真里，你可以直接用理想电流源设定归算后的二次电流，但如果你要模拟真实CT，变比必须选对。
在Simulink里搭建时，建议按以下顺序操作：
- 用“Three-Phase Source”模块建立电源，根据短路容量设定内阻抗（比如短路容量500 MVA，则Zsys = V²/Ssc）。
- 用“Three-Phase Transformer (Two Windings)”模块建立变压器，接线组别设为Yg/Delta-11，短路阻抗设为0.105 pu。
- 两侧分别接“Three-Phase V-I Measurement”模块采集电流。
- 用“Powergui”模块设为离散仿真模式，采样步长建议取50 μs（对应20 kHz采样率），这个步长对工频50 Hz信号来说每周期1000点，足够捕捉暂态细节。
- 故障用“Three-Phase Fault”模块模拟，设置故障电阻和故障时间。
注意：变压器模块的“Saturable”选项如果勾选，需要额外设置饱和特性参数；如果不勾选，涌流工况下差流波形会缺失二次谐波分量，二次谐波闭锁逻辑就无法验证。
4.2 保护判据的Simulink实现与参数整定
保护判据的实现建议用“MATLAB Function”模块，输入各侧归算后的电流瞬时值，输出跳闸信号。核心代码逻辑如下：
function trip = differential_protection(i1a, i1b, i1c, i2a, i2b, i2c) % 输入为两侧三相瞬时电流（已归算到同一基准） % 输出为跳闸信号（1=动作，0=不动作） persistent Id_rms Ir_rms; if isempty(Id_rms) Id_rms = 0; Ir_rms = 0; end % 计算差流和制动电流（瞬时值） id_a = i1a + i2a; id_b = i1b + i2b; id_c = i1c + i2c; ir_a = (abs(i1a) + abs(i2a)) / 2; ir_b = (abs(i1b) + abs(i2b)) / 2; ir_c = (abs(i1c) + abs(i2c)) / 2; % 取三相最大差流和对应制动电流 Id = max([abs(id_a), abs(id_b), abs(id_c)]); Ir = max([ir_a, ir_b, ir_c]); % 比率制动判据 Iop0 = 0.3; % 启动门槛（标幺值，基于CT二次额定电流） K = 0.5; % 制动斜率 Iop = Iop0 + K * Ir; if Id > Iop trip = 1; else trip = 0; end end这段代码里，Iop0和K的整定需要根据保护装置的说明书来。典型值：Iop0取0.3~0.5倍额定电流，K取0.3~0.5。如果K取0.5，意味着制动电流每增加1倍额定值，动作门槛提高0.5倍——这个斜率对应的是区外故障时差流随制动电流增长的速率上限。
参数整定后，必须做三组验证：正常运行（差流应接近零）、区内故障（差流应超过门槛）、区外故障（差流应低于门槛）。如果任何一组不通过，先检查相位补偿和幅值归算，再检查制动电流定义是否与K值匹配。
4.3 二次谐波闭锁与涌流工况的仿真验证
励磁涌流是变压器差动保护误动的头号元凶。涌流波形含有大量二次谐波（占比通常15%~50%），而故障电流的二次谐波含量一般低于10%。因此保护装置通常设置二次谐波闭锁判据：当差流中二次谐波占比超过阈值（典型值15%~20%）时，闭锁差动保护。
在Simulink里实现二次谐波闭锁，需要对差流做傅里叶分析提取二次谐波分量。可以用“Fourier”模块，也可以用MATLAB Function里调用FFT。建议用滑动窗FFT，窗长取一个工频周期（20 ms），这样每个采样点都能输出最新的谐波分析结果。
仿真涌流工况时，变压器的合闸角很关键。合闸角为0度（电压过零点合闸）时涌流最大，合闸角为90度（电压峰值合闸）时涌流最小。建议在仿真里设置合闸角为0度，同时给铁芯设置一定的剩磁（比如0.5 pu），这样涌流波形最恶劣，最能检验闭锁逻辑是否可靠。
我实测下来，如果变压器模型的饱和参数设置合理，合闸角0度、剩磁0.5 pu条件下，差流二次谐波占比可以到30%以上，闭锁逻辑应该可靠动作。但如果你的变压器模型用的是线性模型，二次谐波占比几乎为零，闭锁逻辑永远不会触发——这就是为什么前面强调涌流仿真必须用饱和变压器模型。
5. 常见问题与排查技巧实录
5.1 误动与拒动的典型原因速查表
| 现象 | 可能原因 | 排查方法 | 解决措施 |
|---|---|---|---|
| 正常运行时差流不为零 | 相位补偿方向错误或漏乘√3系数 | 检查补偿前后相量关系，正常运行时差流应小于0.05倍额定 | 修正补偿矩阵，确认Y侧补偿到Δ侧 |
| 正常运行时差流偏大但波形对称 | CT变比归算错误 | 核对两侧二次额定电流是否归算到同一基准 | 重新计算归算系数 |
| 区内故障拒动 | 制动电流定义与K值不匹配，制动量过大 | 打印Ir和Iop曲线，看故障时Iop是否超过Id | 调整K值或更换制动电流定义方式 |
| 区外故障误动 | CT饱和导致虚假差流，或制动斜率K偏小 | 观察区外故障时两侧CT二次电流波形是否畸变 | 增大K值，或加入CT饱和补偿逻辑 |
| 涌流工况误动 | 二次谐波闭锁未生效或阈值过高 | 分析差流二次谐波占比 | 降低闭锁阈值至15%，确认变压器模型含饱和特性 |
| 仿真波形振荡异常 | 仿真步长过大或Powergui离散模式设置不当 | 检查步长是否小于1/(20*f) | 步长取50 μs以下，必要时用变步长求解器 |
5.2 几个让我调了半天的“隐形坑”
第一个坑：Simulink的电流测量模块方向。“Three-Phase V-I Measurement”模块输出的电流方向默认是从端口流向线路，如果你两侧都按同一方向接线，差流计算时符号会反——一侧电流为正，另一侧也为正，加起来是两倍而不是零。解决方法是把一侧的测量模块方向翻转，或者在代码里对一侧电流取负。这个坑隐蔽性极强，因为波形看起来“有差流”但不知道是符号问题。
第二个坑：离散仿真步长与采样率混淆。Powergui里设的“Sample time”是仿真步长，而保护算法里的采样率是另一回事。如果你在MATLAB Function里用persistent变量做滑动窗，窗长按采样点数算，必须确保仿真步长和算法采样率一致。我见过有人在Powergui里设50 μs步长，但算法里按1 ms采样做FFT，结果谐波分析完全对不上。
第三个坑：故障初相角的影响。三相故障模块的故障初相角如果设为默认值，故障电流中可能不含衰减直流分量，CT饱和仿真就不真实。建议把故障初相角设为电压过零点附近（比如0度），这样故障电流有最大直流偏移，最能考验CT模型。
第四个坑：变压器接线组别在Simulink里的表示方式。Simulink变压器模块的接线组别用“Yg”、“Delta1”等字符串表示，但不同版本对组别编号的定义可能有差异。Y/Δ-11在有些版本里要写成“Yg/Delta11”，有些写成“Yg/D11”。如果组别设错，相位补偿的30度方向就反了，差流永远不为零。建议搭好模型后先用额定电流跑一遍，确认差流接近零再往下做。
5.3 仿真结果的可信度校验方法
模型跑通不等于结果可信。我通常用三个方法交叉验证：
方法一：稳态差流校验。正常运行时，差流有效值应小于0.05倍额定电流（考虑CT误差和计算误差）。如果超过0.1倍，说明补偿或归算有问题。
方法二：区外故障制动量校验。在变压器出口处设三相短路，观察制动电流Ir是否随故障电流线性增长，差流Id是否保持在门槛以下。如果Id随Ir同步增长，说明制动斜率K不够大。
方法三：与保护装置录波对比。如果你手头有实际保护装置的故障录波数据（哪怕是历史案例），把同样的电流波形导入仿真模型，看保护判据的输出是否与录波中的动作行为一致。这是最硬核的校验方式，能暴露模型里所有理想化假设带来的偏差。
6. 仿真之外的几点经验
做变压器差动保护仿真，最耗时的部分从来不是搭模型，而是搞清楚“为什么波形不对”。我自己的习惯是：每加一个功能模块（CT饱和、谐波闭锁、制动特性），都先单独验证这个模块的输入输出关系，确认无误后再接入主模型。这样出问题时排查范围小，不至于牵一发动全身。
另外，仿真参数不要追求“一次到位”。先用理想条件跑通逻辑，再逐步加入非理想因素，每加一个就记录差流和制动电流的变化。这个渐进过程本身就是对保护原理的再理解。很多人一上来就把CT饱和、涌流、故障电阻全加上，结果波形一团乱，根本分不清哪个因素在起作用。
最后分享一个我常用的调试技巧：在Simulink里把差流Id、制动电流Ir、动作门槛Iop三条曲线画在同一个示波器上，用不同颜色区分。正常运行时Id贴近零轴，Iop在门槛线上方，Ir在中间；区内故障时Id迅速穿过Iop；区外故障时Id和Ir同步上升但Id始终在Iop下方。这三条线的相对关系一眼就能看出保护行为对不对，比看单条波形直观得多。
