尧图网站建设 
1. 纵联差动保护为什么非得用“差动制动”两条腿走路群里经常有人问线路纵联差动保护教材上几句话就讲完了到了Simulink仿真却怎么都调不对要么区外故障乱跳要么区内故障拒动。我自己当年做这个方向的时候也一样把两侧电流一相减就以为完事了结果波形出来完全不是那么回事。这篇内容就围绕线路纵联差动保护的Simulink仿真展开把原理、建模、波形分析和影响因素串到一起给正在做课程设计、毕业设计或者打算发小论文的同学一个可以直接上手的参考。1.1 从基尔霍夫定律说起内部故障与外部故障的本质区别先聊原理。纵联差动保护的核心依据是基尔霍夫电流定律对于一个无分支的线路段正常运行时流入电流等于流出电流。保护装在被保护线路的两端M侧和N侧各有一套电流互感器两侧的电流信号通过光纤通道或导引线送到同一个保护装置里做比较。这里有一个关键约定规定电流的正方向。通常规定从母线流向被保护线路为正方向。两侧电流互感器的极性端朝向母线或者说一次电流从母线侧流入互感器。在这个约定下正常运行或区外故障时M侧电流从母线流向线路N侧电流从线路流向母线两侧电流方向相反相量和近似为零。区内故障时短路点从两端都获得电流两侧电流方向相同相量和就是故障电流。所以最朴素的做法就是构造差动电流 ( I_d |\dot{I}_M \dot{I}_N| )。区内故障时 ( I_d ) 很大区外故障和正常运行时 ( I_d ) 接近零。这就是“差动”二字的由来。但问题紧接着就来了。区外故障时穿越电流可能达到几千安培如果两侧电流互感器特性不完全一致哪怕差个0.5%差动电流也有几十安培。要是保护门槛设得太低区外故障就会误动设得太高区内轻微故障又可能拒动。这正是“制动”量介入的原因。1.2 比率制动特性不光是躲穿越电流更是躲CT饱和教科书里常见的比率制动判据可以写成[ I_d I_{set0}, \quad I_d K_r \cdot I_{res} ]其中制动电流一般取 ( I_{res} |\dot{I}_M - \dot{I}_N| / 2 )有的厂家也取最大值或者和值的一半但原理一样。简单理解区外故障时穿越电流大两侧电流几乎相等但方向相反( \dot{I}_M - \dot{I}N ) 的模接近两倍穿越电流所以制动电流 ( I{res} ) 自然会变大而差动电流只有两侧不平衡带来的很小分量。这样一来工作点落在制动区内保护不动作。而区内故障时两侧电流方向一致差动电流是故障电流制动电流反而比较小工作点落在动作区保护灵敏动作。顺带一提纵联差动保护比距离保护、过流保护的优势在于它是“全线速动”的——不需要等待对端信号也不存在距离保护的末端灵敏度问题线路内部任意一点故障都能瞬时切除。这也是为什么高压输电线路普遍把它作为主保护而把距离保护作为后备保护。明白了这一层再看Simulink建模思路就清晰多了。2. Simulink模型搭建先把“被保护线路”这个舞台搭起来说句实话原理再熟一到Simulink里找不到模块、参数不会设照样卡壳。我用的环境是MATLAB R2022b专用的电力系统模块集在Simscape Electrical Specialized Power Systems下面旧版本叫SimPowerSystems操作方法大同小异。2.1 主电路参数设置电源、线路、CT怎么选先搭一个220kV双端供电系统。打开Simulink库浏览器依次拿以下模块两端电源Two Three-Phase Source线电压有效值设为220kV频率50Hz。为了让区外故障和区内故障场景都能看清楚两端电源的相位可以差一点比如M侧相位0度N侧相位滞后10度模拟输送功率的工况。线路用Three-Phase Distributed Parameters Line长度设100km。这里不推荐用PI型集中参数模型做差动保护的仿真因为长线路的分布电容对差动电流有明显影响而PI模型只能近似模拟分布效应做影响因素分析时误差会误导结论。电流测量每一侧用Three-Phase V-I Measurement接在断路器出口和线路之间电压电流一起测出来后面保护逻辑要从这里取电流信号。电源和线路的参数我直接给一组实测好用的M侧电源220kV50Hz三相短路容量5000MVAX/R比10N侧电源220kV50Hz三相短路容量3000MVAX/R比10。线路正序参数R10.01273Ω/kmL10.9337mH/kmC112.74nF/km。注意一个细节三个电压源中性点接地方式要选对。220kV系统中性点直接接地所以要打开电源模块的“neutral”设置否则做单相接地短路时故障回路不闭合波形会非常怪。2.2 故障与测量环节短路点怎么放、电流方向怎么标故障模块用Three-Phase Fault可以设置接地电阻和相间电阻还能通过外部信号控制故障时刻。我这里习惯把故障时间设置在0.15s持续0.1s0.3s清除。这个时间段足够保护判断动作了。关键点来了故障点位置决定了是区内还是区外。我的做法是区内故障故障点放在线路上的0.6倍长度处这是保护范围正中间用来验证灵敏动作。区外故障故障点放在M侧母线与M侧断路器之间模拟M侧区外短路此时经过线路的是穿越性故障电流。电流极性问题绝对要重视。差动保护的核心是电流方向比较所以测量的电流方向必须与咱们定义的“母线流向线路为正”一致。在Three-Phase V-I Measurement模块中电流的正方向是“从正端子流入从负端子流出”。搭模型时让母线侧连接测量模块的正端子线路侧连接负端子这样测出来的电流就是线路正方向电流。两侧都要保持一致否则你会发现区外故障时差动电流比区内故障还大——这几乎是我见过最多人踩的坑。2.3 保护逻辑模型从瞬时值到比率差动判据的完整链路电流信号从测量模块出来先送到Three-Phase RMS模块取有效值再进入保护逻辑。但RMS模块会有延迟对瞬时动作的纵联保护仿真不够精准。更专业的做法是用Discrete Fourier模块提取工频基波分量的幅值和相位然后按相计算。保护逻辑部分的搭建可以复用Simulink基本数学模块用减法模块( I_M - I_N )再用绝对值模块。用加法模块( I_M I_N )再用绝对值模块。差动电流 ( I_d |I_M I_N| )。制动电流 ( I_{res} |I_M - I_N| / 2 )这里加个Gain增益0.5。判断部分我建议用比较模块加上与门。判断条件拆成两路一路判断 ( I_d ) 是否大于最小动作门槛 ( I_{set0} )比如200A按额定电流1500A、CT变比1000/5折算过来的话大概是整定值的0.3倍另一路判断 ( I_d ) 是否大于 ( K_r \times I_{res} )其中 ( K_r ) 取0.5~0.7。两路结果相与再经过一个几十毫秒的延时确认输出跳闸信号。为了波形便于观察跳闸信号接一个断路器模块切除线路两侧的断路器。3. 仿真波形解读内部故障、区外故障、正常负荷三种场景仿真模型搭好之后最重要的是会读波形、会验证保护逻辑是否真的按预期动作。这里我分三种典型工况来解读顺便说说怎么判断结果对不对。3.1 区内故障可靠动作的典型波形特征设置故障点在60km处A相接地短路故障电阻0.01Ω故障时间0.15s开始。打开Scope可以看到故障发生前线路潮流不大两侧电流相位差10度左右差动电流很小只有几十安培的不平衡分量。故障瞬间两侧电流都显著增大M侧和N侧电流波形相位接近同向差动电流 ( I_d ) 从几十安培跃升到几千安培制动电流 ( I_{res} ) 增大但不及 ( I_d ) 的增幅工作点越过比率制动特性曲线进入动作区。保护逻辑输出的跳闸信号在故障后大约10~20ms内变为1随后断路器断开两侧电流归零。这里要注意从故障发生到保护出口总会有延时仿真里这个时间取决于你设置的延时模块和数值积分步长。实际工程中纵联差动保护的整组动作时间一般在20~40ms仿真结果在这个范围内就是合理的。3.2 区外故障比率制动如何把误动压下去把故障点改到M侧母线附近M侧断路器与电源之间还是A相接地故障时间0.15s。这个场景是检验比率制动特性的试金石。故障发生后穿越电流很大M侧和N侧电流的幅值都迅速增大。但注意方向两侧电流相位相反差动电流 ( I_d ) 很小而制动电流 ( I_{res} ) 变得非常大工作点落在制动特性曲线的制动区内保护不动作。我实测的时候发现一个很有意思的现象如果两侧CT模型参数完全一致、电流测量也没有误差区外故障的差动电流几乎接近零这其实是理想化的结果。现实中总会有点不平衡但正因为有了制动电流这个分母工作点依然离动作边界很远。这就是比率制动“以穿越电流的大小动态调整动作门槛”的精髓。3.3 正常负荷与轻微不平衡保护不应动作的底线正常负荷工况下两侧电流相位差来源于两端电源功角差电流幅值在几百安培量级。由于两侧测量的电流是同一个电流只是一端流入、一端流出差动电流理论上为零实际上只有数值计算引入的微小误差。保护逻辑输出保持为0波形稳定没有误动风险。需要特别留意的还有一种“微弱故障”工况比如故障点经100Ω高阻接地。这种故障电流很小可能只有几百安培差动电流也很小如果最小动作门槛 ( I_{set0} ) 设得过高就直接拒动了。仿真时建议专门做一组高阻接地试验用来体会“灵敏度”与“可靠性”之间的权衡这也是后面讲影响因素时的重要伏笔。4. 相关影响因素探究为什么同一套方案换个工况就翻车标题里“相关影响因素”这五个字其实是整个项目里最值得写进论文的部分。很多同学做完基础仿真就停了但评阅老师真正想看的是你能不能分析出哪些因素会改变保护的行为并在模型里定量复现。我逐个说。4.1 分布电容电流高压长线路里的“隐形差流”输电线路对地电容是客观存在的电压等级越高、线路越长电容电流越大。正常工作电压下线路两端都会流出一部分“充电电流”这部分电流流进线路后并不从另一端流出而是“漏”到大地里去了。问题就出在这差动保护的两端电流相加并不完全为零因为有一小部分电流被电容吃掉了形成没有故障时的差动电流 ( I_{cap} )。电压越高、线路越长( I_{cap} ) 越大。对330kV及以上级别的长线路这个值可能达到几十安培甚至上百安培。如果不做补偿固定门槛就会误动。仿真里验证这个影响很直接把线路长度从60km放宽到300km正常负荷工况下观察差动电流的稳态值。你会清楚地看到它随线路长度上升。要消除它可以加电容电流补偿回路也就是用电压测量值乘以电容参数算出补偿电流再从差动电流里扣掉。不过补偿精度总有限度所以实际工程对超长线路还会提高最小动作门槛牺牲一点点灵敏度换可靠性。4.2 CT饱和与变比不一致制动特性的设计初衷前面说过只要两侧CT特性不一致区外故障的大穿越电流就会在二次侧产生虚假差流。最极端的场景是大穿越电流引起一侧CT饱和饱和一侧的二次电流波形出现削顶变形幅值明显下降于是两侧二次电流不再对称差动电流变大。在Simulink里复现这个现象有两种办法。一种是给两侧电流测量模块分别乘一个不同的变比系数模拟变比不一致。比如M侧增益1.0N侧增益1.02区外故障时差动电流立刻就从零变成穿越电流的2%效果非常直观。另一种是用带饱和特性的CT模型库里有专门的Current Transformer模块封装可以设置饱和拐点电压和剩磁系数更贴近物理实际。比率制动正是干这个事的。穿越电流大制动电流跟着大即使差动电流因为CT饱和涨到几百安培工作点也还落在制动区。做论文的同学这个场景一定要做它最能体现你对制动特性的理解深度。4.3 两端数据同步误差物理仿真里藏得最深的坑纵联差动保护需要比较两端电流的相位所以两端采样数据必须严格同步。实际工程里靠的是GPS对时或者光纤通道的同步算法误差通常控制在几十微秒以内。但别小看这几十微秒折算到工频50Hz那就是1~2度的相位误差。为什么相位误差会转化成差动电流把区外故障看成两端电流相位相反如果一端的数据延迟了1ms算出来的 ( \dot{I}_M \dot{I}_N ) 就不再是零而是一个与穿越电流成正比的虚假差流。穿越电流越大这个误差越明显。仿真里可以用一个Variable Time Delay或者Transport Delay模块模拟这个误差设1ms延时异步端电流然后观察区外故障差动电流的变化。如果你的保护门槛整定得很低这个误差就有可能造成误动。这也是为什么工程整定计算里要预留“可靠性系数”——所有误差源产生的总不平衡电流都被折算进整定值里。4.4 故障电阻与故障类型的影响故障电阻的影响前面简单提到过。金属性短路时故障电流大、差动电流大保护动作可靠。但高阻接地时故障电流可能不足金属性短路的十分之一差动电流也随之缩水灵敏度压力全在最小动作门槛 ( I_{set0} ) 的整定上。故障类型的影响也要仿真验证一下。三相短路时三相差动电流都很大ABC三相保护全部动作逻辑上没有任何区分的必要。但单相接地故障时只有故障相的动作量显著变化另外两相差动电流基本不变这正好体现了分相式差动保护“选相跳闸”的优势。Simulink的波形上你可以清楚看到非故障相的差动电流在故障前后几乎一模一样这个特征可以做进保护逻辑里做选相也是论文里一个不错的创新点。影响因素这块内容比较多我把要点整理成一张表方便写论文时直接引用影响因素作用机制典型表现Simulink模拟方法工程应对分布电容电流充电电流流入大地两侧电流不平衡线路越长差流越大正常时就有稳态差流加长线路长度电容电流补偿提高门槛CT饱和/变比不一致二次侧电流畸变两侧不平衡区外故障时出现大的虚假差流修改增益系数或使用饱和CT模型比率制动合理整定可靠性系数数据同步误差两端电流相位比较失真穿越电流大时差流显著增大Transport Delay加延时高精度同步通道预留整定裕度故障电阻故障电流被电阻限制量值减小高阻故障时差流小可能低于门槛加大Fault模块的接地电阻降低最小门槛提高灵敏度故障类型不同相别的故障特征不同单相故障仅有对应相差流变化切换Fault模块故障类型分相差动选相跳闸5. 调试心得仿真发散、极性错误、逻辑不动作的排错思路模型搭好了不等于就能跑出理想波形。这一章把我这些年仿真中反复踩过、反复排掉的坑集中说一下。5.1 仿真发散怎么办求解器与步长的组合Simulink做电力系统电磁暂态仿真最常用的求解器是ode23tb或ode15s。原因很简单电力电子开关、故障瞬间电流突变这些事件让系统呈现明显的刚性特征显式求解器极易发散或者步长小到跑不动。如果模型一运行就发散我的排查顺序是先看有没有代数环。比如直接用电流信号去控制断路器而断路器又影响电流这是典型的代数环可以在模型里用Memory模块或者把信号延迟一拍打断它。再看求解器类型。把solver设为ode23tbMax Step Size设成1e-4基本能解决90%的解算问题。最后看线路参数。分布式参数线路如果和电源的短路容量不匹配初始化时会给出警告要检查电源的串联阻抗是不是合理绝不能是零。5.2 保护不动作的第一排查顺序极性→变比→判据方向保护逻辑明明加了可故障发生后跳闸信号就是不出来。这种情况我最先查的是电流极性不是判据参数。把两侧电流的波形放在同一个Scope里叠加显示区外故障时两条波形应该相位相反区内故障时相位相同。如果发现反了把测量模块的正负端对调一下就好。第二步查变比一致性。M侧和N侧各有一个Gain模块用于模拟CT变比数值必须一致否则区外故障就会触发误动而区内故障时判据也可能因为不平衡电流的方向问题被制动量压住。第三步查判据方向。检查逻辑表达式里 ( |I_M I_N| ) 是不是真的连接成加法关系。很多人画图时不小心把加法模块的符号搞反变成了 ( |I_M - I_N| )那差动和制动就完全颠倒了保护逻辑就成了一个“穿越电流越大越爱动”的错误装置。5.3 用MATLAB脚本验证保护逻辑比看Scope更靠谱Scope看图终究有主观性而且波形细节看不清。我建议把关键信号通过To Workspace模块导出到MATLAB工作区比如导出差动电流 ( I_d )、制动电流 ( I_{res} )、跳闸信号 ( trip )然后在MATLAB脚本里做后处理分析。比如可以在脚本里画出差动电流随时间变化的曲线同时用一条水平虚线画出门槛值很直观就能看出故障后的差流是不是超过门槛。更进一步还可以把所有采样点的 ( (I_{res}, I_d) ) 画在比率制动特性平面上把动作边界和实际轨迹一起画出来观察工作点是进入动作区还是制动区。这一步写进论文里比贴十张Scope截图更有说服力。根据我个人经验来说纵联差动保护的Simulink仿真想做好功夫花在“原理定方向、建模抠细节、分析定深度”这三层上。方向层面搞清差动和制动的配合逻辑建模层面盯住CT极性和测量方向这类看似不起眼却致命的小细节分析层面用影响因素和参数扫描来体现对保护行为的整体把握一篇像样的毕业论文或技术报告就出来了。最后再分享一个小技巧把故障电阻、线路长度、CT变比误差这几个参数分别做成变量用MATLAB循环批量跑仿真把所有结果汇总成一张趋势表你会看到保护性能随参数变化的全貌这部分数据写进论文里非常出彩。 线路纵联差动保护Simulink仿真：原理、建模与影响因素
1. 纵联差动保护为什么非得用“差动制动”两条腿走路群里经常有人问线路纵联差动保护教材上几句话就讲完了到了Simulink仿真却怎么都调不对要么区外故障乱跳要么区内故障拒动。我自己当年做这个方向的时候也一样把两侧电流一相减就以为完事了结果波形出来完全不是那么回事。这篇内容就围绕线路纵联差动保护的Simulink仿真展开把原理、建模、波形分析和影响因素串到一起给正在做课程设计、毕业设计或者打算发小论文的同学一个可以直接上手的参考。1.1 从基尔霍夫定律说起内部故障与外部故障的本质区别先聊原理。纵联差动保护的核心依据是基尔霍夫电流定律对于一个无分支的线路段正常运行时流入电流等于流出电流。保护装在被保护线路的两端M侧和N侧各有一套电流互感器两侧的电流信号通过光纤通道或导引线送到同一个保护装置里做比较。这里有一个关键约定规定电流的正方向。通常规定从母线流向被保护线路为正方向。两侧电流互感器的极性端朝向母线或者说一次电流从母线侧流入互感器。在这个约定下正常运行或区外故障时M侧电流从母线流向线路N侧电流从线路流向母线两侧电流方向相反相量和近似为零。区内故障时短路点从两端都获得电流两侧电流方向相同相量和就是故障电流。所以最朴素的做法就是构造差动电流 ( I_d |\dot{I}_M \dot{I}_N| )。区内故障时 ( I_d ) 很大区外故障和正常运行时 ( I_d ) 接近零。这就是“差动”二字的由来。但问题紧接着就来了。区外故障时穿越电流可能达到几千安培如果两侧电流互感器特性不完全一致哪怕差个0.5%差动电流也有几十安培。要是保护门槛设得太低区外故障就会误动设得太高区内轻微故障又可能拒动。这正是“制动”量介入的原因。1.2 比率制动特性不光是躲穿越电流更是躲CT饱和教科书里常见的比率制动判据可以写成[ I_d I_{set0}, \quad I_d K_r \cdot I_{res} ]其中制动电流一般取 ( I_{res} |\dot{I}_M - \dot{I}_N| / 2 )有的厂家也取最大值或者和值的一半但原理一样。简单理解区外故障时穿越电流大两侧电流几乎相等但方向相反( \dot{I}_M - \dot{I}N ) 的模接近两倍穿越电流所以制动电流 ( I{res} ) 自然会变大而差动电流只有两侧不平衡带来的很小分量。这样一来工作点落在制动区内保护不动作。而区内故障时两侧电流方向一致差动电流是故障电流制动电流反而比较小工作点落在动作区保护灵敏动作。顺带一提纵联差动保护比距离保护、过流保护的优势在于它是“全线速动”的——不需要等待对端信号也不存在距离保护的末端灵敏度问题线路内部任意一点故障都能瞬时切除。这也是为什么高压输电线路普遍把它作为主保护而把距离保护作为后备保护。明白了这一层再看Simulink建模思路就清晰多了。2. Simulink模型搭建先把“被保护线路”这个舞台搭起来说句实话原理再熟一到Simulink里找不到模块、参数不会设照样卡壳。我用的环境是MATLAB R2022b专用的电力系统模块集在Simscape Electrical Specialized Power Systems下面旧版本叫SimPowerSystems操作方法大同小异。2.1 主电路参数设置电源、线路、CT怎么选先搭一个220kV双端供电系统。打开Simulink库浏览器依次拿以下模块两端电源Two Three-Phase Source线电压有效值设为220kV频率50Hz。为了让区外故障和区内故障场景都能看清楚两端电源的相位可以差一点比如M侧相位0度N侧相位滞后10度模拟输送功率的工况。线路用Three-Phase Distributed Parameters Line长度设100km。这里不推荐用PI型集中参数模型做差动保护的仿真因为长线路的分布电容对差动电流有明显影响而PI模型只能近似模拟分布效应做影响因素分析时误差会误导结论。电流测量每一侧用Three-Phase V-I Measurement接在断路器出口和线路之间电压电流一起测出来后面保护逻辑要从这里取电流信号。电源和线路的参数我直接给一组实测好用的M侧电源220kV50Hz三相短路容量5000MVAX/R比10N侧电源220kV50Hz三相短路容量3000MVAX/R比10。线路正序参数R10.01273Ω/kmL10.9337mH/kmC112.74nF/km。注意一个细节三个电压源中性点接地方式要选对。220kV系统中性点直接接地所以要打开电源模块的“neutral”设置否则做单相接地短路时故障回路不闭合波形会非常怪。2.2 故障与测量环节短路点怎么放、电流方向怎么标故障模块用Three-Phase Fault可以设置接地电阻和相间电阻还能通过外部信号控制故障时刻。我这里习惯把故障时间设置在0.15s持续0.1s0.3s清除。这个时间段足够保护判断动作了。关键点来了故障点位置决定了是区内还是区外。我的做法是区内故障故障点放在线路上的0.6倍长度处这是保护范围正中间用来验证灵敏动作。区外故障故障点放在M侧母线与M侧断路器之间模拟M侧区外短路此时经过线路的是穿越性故障电流。电流极性问题绝对要重视。差动保护的核心是电流方向比较所以测量的电流方向必须与咱们定义的“母线流向线路为正”一致。在Three-Phase V-I Measurement模块中电流的正方向是“从正端子流入从负端子流出”。搭模型时让母线侧连接测量模块的正端子线路侧连接负端子这样测出来的电流就是线路正方向电流。两侧都要保持一致否则你会发现区外故障时差动电流比区内故障还大——这几乎是我见过最多人踩的坑。2.3 保护逻辑模型从瞬时值到比率差动判据的完整链路电流信号从测量模块出来先送到Three-Phase RMS模块取有效值再进入保护逻辑。但RMS模块会有延迟对瞬时动作的纵联保护仿真不够精准。更专业的做法是用Discrete Fourier模块提取工频基波分量的幅值和相位然后按相计算。保护逻辑部分的搭建可以复用Simulink基本数学模块用减法模块( I_M - I_N )再用绝对值模块。用加法模块( I_M I_N )再用绝对值模块。差动电流 ( I_d |I_M I_N| )。制动电流 ( I_{res} |I_M - I_N| / 2 )这里加个Gain增益0.5。判断部分我建议用比较模块加上与门。判断条件拆成两路一路判断 ( I_d ) 是否大于最小动作门槛 ( I_{set0} )比如200A按额定电流1500A、CT变比1000/5折算过来的话大概是整定值的0.3倍另一路判断 ( I_d ) 是否大于 ( K_r \times I_{res} )其中 ( K_r ) 取0.5~0.7。两路结果相与再经过一个几十毫秒的延时确认输出跳闸信号。为了波形便于观察跳闸信号接一个断路器模块切除线路两侧的断路器。3. 仿真波形解读内部故障、区外故障、正常负荷三种场景仿真模型搭好之后最重要的是会读波形、会验证保护逻辑是否真的按预期动作。这里我分三种典型工况来解读顺便说说怎么判断结果对不对。3.1 区内故障可靠动作的典型波形特征设置故障点在60km处A相接地短路故障电阻0.01Ω故障时间0.15s开始。打开Scope可以看到故障发生前线路潮流不大两侧电流相位差10度左右差动电流很小只有几十安培的不平衡分量。故障瞬间两侧电流都显著增大M侧和N侧电流波形相位接近同向差动电流 ( I_d ) 从几十安培跃升到几千安培制动电流 ( I_{res} ) 增大但不及 ( I_d ) 的增幅工作点越过比率制动特性曲线进入动作区。保护逻辑输出的跳闸信号在故障后大约10~20ms内变为1随后断路器断开两侧电流归零。这里要注意从故障发生到保护出口总会有延时仿真里这个时间取决于你设置的延时模块和数值积分步长。实际工程中纵联差动保护的整组动作时间一般在20~40ms仿真结果在这个范围内就是合理的。3.2 区外故障比率制动如何把误动压下去把故障点改到M侧母线附近M侧断路器与电源之间还是A相接地故障时间0.15s。这个场景是检验比率制动特性的试金石。故障发生后穿越电流很大M侧和N侧电流的幅值都迅速增大。但注意方向两侧电流相位相反差动电流 ( I_d ) 很小而制动电流 ( I_{res} ) 变得非常大工作点落在制动特性曲线的制动区内保护不动作。我实测的时候发现一个很有意思的现象如果两侧CT模型参数完全一致、电流测量也没有误差区外故障的差动电流几乎接近零这其实是理想化的结果。现实中总会有点不平衡但正因为有了制动电流这个分母工作点依然离动作边界很远。这就是比率制动“以穿越电流的大小动态调整动作门槛”的精髓。3.3 正常负荷与轻微不平衡保护不应动作的底线正常负荷工况下两侧电流相位差来源于两端电源功角差电流幅值在几百安培量级。由于两侧测量的电流是同一个电流只是一端流入、一端流出差动电流理论上为零实际上只有数值计算引入的微小误差。保护逻辑输出保持为0波形稳定没有误动风险。需要特别留意的还有一种“微弱故障”工况比如故障点经100Ω高阻接地。这种故障电流很小可能只有几百安培差动电流也很小如果最小动作门槛 ( I_{set0} ) 设得过高就直接拒动了。仿真时建议专门做一组高阻接地试验用来体会“灵敏度”与“可靠性”之间的权衡这也是后面讲影响因素时的重要伏笔。4. 相关影响因素探究为什么同一套方案换个工况就翻车标题里“相关影响因素”这五个字其实是整个项目里最值得写进论文的部分。很多同学做完基础仿真就停了但评阅老师真正想看的是你能不能分析出哪些因素会改变保护的行为并在模型里定量复现。我逐个说。4.1 分布电容电流高压长线路里的“隐形差流”输电线路对地电容是客观存在的电压等级越高、线路越长电容电流越大。正常工作电压下线路两端都会流出一部分“充电电流”这部分电流流进线路后并不从另一端流出而是“漏”到大地里去了。问题就出在这差动保护的两端电流相加并不完全为零因为有一小部分电流被电容吃掉了形成没有故障时的差动电流 ( I_{cap} )。电压越高、线路越长( I_{cap} ) 越大。对330kV及以上级别的长线路这个值可能达到几十安培甚至上百安培。如果不做补偿固定门槛就会误动。仿真里验证这个影响很直接把线路长度从60km放宽到300km正常负荷工况下观察差动电流的稳态值。你会清楚地看到它随线路长度上升。要消除它可以加电容电流补偿回路也就是用电压测量值乘以电容参数算出补偿电流再从差动电流里扣掉。不过补偿精度总有限度所以实际工程对超长线路还会提高最小动作门槛牺牲一点点灵敏度换可靠性。4.2 CT饱和与变比不一致制动特性的设计初衷前面说过只要两侧CT特性不一致区外故障的大穿越电流就会在二次侧产生虚假差流。最极端的场景是大穿越电流引起一侧CT饱和饱和一侧的二次电流波形出现削顶变形幅值明显下降于是两侧二次电流不再对称差动电流变大。在Simulink里复现这个现象有两种办法。一种是给两侧电流测量模块分别乘一个不同的变比系数模拟变比不一致。比如M侧增益1.0N侧增益1.02区外故障时差动电流立刻就从零变成穿越电流的2%效果非常直观。另一种是用带饱和特性的CT模型库里有专门的Current Transformer模块封装可以设置饱和拐点电压和剩磁系数更贴近物理实际。比率制动正是干这个事的。穿越电流大制动电流跟着大即使差动电流因为CT饱和涨到几百安培工作点也还落在制动区。做论文的同学这个场景一定要做它最能体现你对制动特性的理解深度。4.3 两端数据同步误差物理仿真里藏得最深的坑纵联差动保护需要比较两端电流的相位所以两端采样数据必须严格同步。实际工程里靠的是GPS对时或者光纤通道的同步算法误差通常控制在几十微秒以内。但别小看这几十微秒折算到工频50Hz那就是1~2度的相位误差。为什么相位误差会转化成差动电流把区外故障看成两端电流相位相反如果一端的数据延迟了1ms算出来的 ( \dot{I}_M \dot{I}_N ) 就不再是零而是一个与穿越电流成正比的虚假差流。穿越电流越大这个误差越明显。仿真里可以用一个Variable Time Delay或者Transport Delay模块模拟这个误差设1ms延时异步端电流然后观察区外故障差动电流的变化。如果你的保护门槛整定得很低这个误差就有可能造成误动。这也是为什么工程整定计算里要预留“可靠性系数”——所有误差源产生的总不平衡电流都被折算进整定值里。4.4 故障电阻与故障类型的影响故障电阻的影响前面简单提到过。金属性短路时故障电流大、差动电流大保护动作可靠。但高阻接地时故障电流可能不足金属性短路的十分之一差动电流也随之缩水灵敏度压力全在最小动作门槛 ( I_{set0} ) 的整定上。故障类型的影响也要仿真验证一下。三相短路时三相差动电流都很大ABC三相保护全部动作逻辑上没有任何区分的必要。但单相接地故障时只有故障相的动作量显著变化另外两相差动电流基本不变这正好体现了分相式差动保护“选相跳闸”的优势。Simulink的波形上你可以清楚看到非故障相的差动电流在故障前后几乎一模一样这个特征可以做进保护逻辑里做选相也是论文里一个不错的创新点。影响因素这块内容比较多我把要点整理成一张表方便写论文时直接引用影响因素作用机制典型表现Simulink模拟方法工程应对分布电容电流充电电流流入大地两侧电流不平衡线路越长差流越大正常时就有稳态差流加长线路长度电容电流补偿提高门槛CT饱和/变比不一致二次侧电流畸变两侧不平衡区外故障时出现大的虚假差流修改增益系数或使用饱和CT模型比率制动合理整定可靠性系数数据同步误差两端电流相位比较失真穿越电流大时差流显著增大Transport Delay加延时高精度同步通道预留整定裕度故障电阻故障电流被电阻限制量值减小高阻故障时差流小可能低于门槛加大Fault模块的接地电阻降低最小门槛提高灵敏度故障类型不同相别的故障特征不同单相故障仅有对应相差流变化切换Fault模块故障类型分相差动选相跳闸5. 调试心得仿真发散、极性错误、逻辑不动作的排错思路模型搭好了不等于就能跑出理想波形。这一章把我这些年仿真中反复踩过、反复排掉的坑集中说一下。5.1 仿真发散怎么办求解器与步长的组合Simulink做电力系统电磁暂态仿真最常用的求解器是ode23tb或ode15s。原因很简单电力电子开关、故障瞬间电流突变这些事件让系统呈现明显的刚性特征显式求解器极易发散或者步长小到跑不动。如果模型一运行就发散我的排查顺序是先看有没有代数环。比如直接用电流信号去控制断路器而断路器又影响电流这是典型的代数环可以在模型里用Memory模块或者把信号延迟一拍打断它。再看求解器类型。把solver设为ode23tbMax Step Size设成1e-4基本能解决90%的解算问题。最后看线路参数。分布式参数线路如果和电源的短路容量不匹配初始化时会给出警告要检查电源的串联阻抗是不是合理绝不能是零。5.2 保护不动作的第一排查顺序极性→变比→判据方向保护逻辑明明加了可故障发生后跳闸信号就是不出来。这种情况我最先查的是电流极性不是判据参数。把两侧电流的波形放在同一个Scope里叠加显示区外故障时两条波形应该相位相反区内故障时相位相同。如果发现反了把测量模块的正负端对调一下就好。第二步查变比一致性。M侧和N侧各有一个Gain模块用于模拟CT变比数值必须一致否则区外故障就会触发误动而区内故障时判据也可能因为不平衡电流的方向问题被制动量压住。第三步查判据方向。检查逻辑表达式里 ( |I_M I_N| ) 是不是真的连接成加法关系。很多人画图时不小心把加法模块的符号搞反变成了 ( |I_M - I_N| )那差动和制动就完全颠倒了保护逻辑就成了一个“穿越电流越大越爱动”的错误装置。5.3 用MATLAB脚本验证保护逻辑比看Scope更靠谱Scope看图终究有主观性而且波形细节看不清。我建议把关键信号通过To Workspace模块导出到MATLAB工作区比如导出差动电流 ( I_d )、制动电流 ( I_{res} )、跳闸信号 ( trip )然后在MATLAB脚本里做后处理分析。比如可以在脚本里画出差动电流随时间变化的曲线同时用一条水平虚线画出门槛值很直观就能看出故障后的差流是不是超过门槛。更进一步还可以把所有采样点的 ( (I_{res}, I_d) ) 画在比率制动特性平面上把动作边界和实际轨迹一起画出来观察工作点是进入动作区还是制动区。这一步写进论文里比贴十张Scope截图更有说服力。根据我个人经验来说纵联差动保护的Simulink仿真想做好功夫花在“原理定方向、建模抠细节、分析定深度”这三层上。方向层面搞清差动和制动的配合逻辑建模层面盯住CT极性和测量方向这类看似不起眼却致命的小细节分析层面用影响因素和参数扫描来体现对保护行为的整体把握一篇像样的毕业论文或技术报告就出来了。最后再分享一个小技巧把故障电阻、线路长度、CT变比误差这几个参数分别做成变量用MATLAB循环批量跑仿真把所有结果汇总成一张趋势表你会看到保护性能随参数变化的全貌这部分数据写进论文里非常出彩。