2026/9/10 7:56:29
FEATURED · 精选文章 
220kV线路的相间距离保护一旦到了双端电源的供电方式下整定与动作逻辑就比单端电源复杂得多。我前阵子刚好用Simulink做完了这个方向的仿真研究把220kV双端电源系统里的相间距离保护从模型搭建到故障仿真完整跑了一遍收获不少也踩了一些坑。这篇记录就是给准备做这个方向的人一份可以照着走的参考不堆理论公式重点讲清楚双端电源为什么会影响距离保护判断、三段式保护怎么整定、模型在Simulink里怎么搭、故障结果怎么解读。无论你是刚接手的继保专业学生、要交毕设的本科生还是做培训验收的现场工程师照着这个思路走都能得到一套能跑通、能出结果的仿真模型。做这个课题之前有个前提要先说清楚实际220kV线路通常配置纵联保护做主保护距离保护和零序电流保护作为后备保护。我这套仿真把注意力集中在相间距离保护本身刻意不引入纵联逻辑目的是把三段式动作特性、阻抗测量误差、助增电流影响这些底层问题看清楚。方向选这个题目是因为220kV处于区域电网主网架几乎都是双端或环网供电与110kV及以下单端辐射网的保护分析逻辑差异很大很多初学者在这上面栽跟头。1. 项目核心背景与整体思路拆解1.1 220kV双端电源系统为什么是距离保护的“难度天花板”先说结论220kV输电线路里的双端电源会让距离保护的核心变量——测量阻抗——不再直观。距离保护的基本想法很朴素测出保护安装处的电压和电流算出从安装点到故障点的等效阻抗然后和整定值比较落在阻抗圆内就动作。这个思路在单端电源系统里几乎不会出问题因为故障电流方向唯一电压电流的相位关系很清楚。但双端电源一上来情况就变了。第一重变化是故障电流的方向不再唯一。M侧线路故障时流进故障点的电流是M侧电源和N侧电源共同提供的保护安装处测到的电流只是总故障电流的一部分。如果故障处存在过渡电阻故障点电压会被对侧电流抬高或压低测量阻抗就会偏离真实故障距离对应的阻抗值。第二重变化是助增效应如果保护安装点和故障点之间还存在其他电源支路或环网路径流到故障点的电流大于流过保护的电流测量阻抗会系统性偏大灵敏度下降。反过来还有汲出效应测量阻抗偏小可能造成区外误动。第三重变化是方向问题反向故障时电流电压的相位特征和正向故障完全不同没有方向判别的距离保护在双端系统里寸步难行。这些都不是书本上背个公式就能消化的必须靠仿真把故障过程的电压电流波形完整拉出来看。Simulink里可以很方便地调整双端电源容量、线路长度、分支位置、过渡电阻大小每次调整都能在波形上看到测量阻抗的轨迹变化这对理解双端系统的保护行为帮助极大。1.2 为什么选Simulink而不是PSCAD/RTDS做电力系统电磁暂态仿真圈子里常见的选项无非是PSCAD、RTDS、MATLAB/Simulink这几类。PSCAD在电磁暂态计算精度和元件库丰富度上确实有优势RTDS则是硬件在环测试的主流平台。但这个课题选Simulink我是综合考虑了三件事。第一Simulink的Simscape Electrical组件库原SimPowerSystems涵盖了三相电源、三相线路、三相变压器、故障模块、断路器和各种测量元件对220kV单线双端系统这种规模来说完全够用不需要额外许可证。第二也是最重要的一点这个课题不光是搭个电路模型看波形还要在里面实现距离保护算法——全周傅里叶相量提取、阻抗计算、圆特性判据、三段时限配合。这些逻辑用MATLAB Function或S-Function写非常顺手数据和绘图也与Simulink无缝衔接。PSCAD里写自定义元件就要麻烦不少。第三仿真参数扫描、批量跑故障工况、结果后处理MATLAB脚本在这方面几乎是降维打击我可以一次性把几十种故障位置、故障类型、过渡电阻组合全部跑完自动生成对比曲线。当然Simulink也有短板大规模电力系统电磁暂态的计算效率和收敛性不如专用工具但这个课题只有两条线路加两个电源规模很小完全在可控范围内。如果后续有人想往多机系统、行波保护方向扩展再迁移到PSCAD或RTDS也不迟。2. 保护原理与逻辑搭建测量阻抗与三段式整定2.1 测量阻抗是怎么算出来的0度接线与方向圆判据相间距离保护的核心测量量是相间阻抗。以AB相间短路为例保护接入的电压是线电压Uab电流是线电流Iab也就是常说的0度接线方式。计算得到的测量阻抗Zab Uab / Iab这个复数在R-X平面上对应一个点。正常情况下负荷状态下测量阻抗落在负荷区幅值大、角度接近功率因数角故障后测量阻抗迅速跌落到线路阻抗附近幅值小、角度接近线路阻抗角。保护要做的事情就是判断这个测量阻抗点是否落入预先设定的动作区域。动作区域最常用的是方向圆特性。整定阻抗设为Zset圆心在Zset/2处半径为Zset/2动作判据写成复数的模值比较就是|Zm - Zset/2| ≤ |Zset/2|。方向圆的特点是圆过原点也就是说测量阻抗落在原点附近或者第三象限时不动作天然带有方向性正好应对双端电源里反向故障的问题。还有一种偏移圆特性把圆往第三象限偏移一部分适用于母线背后有较长线路的场景但会牺牲一些方向选择性这个课题里我用的是方向圆。整套保护逻辑还要一个启动元件。不能等阻抗都算出来掉进圆里才启动那样在正常运行时就容易频繁计算误动。启动元件可以用相电流突变量、负序电流分量或者低电压启动。我在仿真里用的是相电流突变量启动连续三个采样点的电流变化量超过设定门槛就进入故障处理流程。门槛值一般按最大负荷电流的20%~50%整定太灵敏容易在负荷波动时误动太迟钝又会影响故障切除速度这个要靠多次仿真去调。2.2 三段式整定计算一个80km线路的完整算例这套仿真我采用了一个典型的220kV双端系统参数M侧系统短路容量较大3000MVAN侧系统相对较弱1000MVA线路全长80km。220kV架空线单位正序参数按常见值取R1 0.027Ω/kmX1 0.28Ω/kmC1 0.014μF/km。这样算下来全线路正序阻抗为ZL 80×(0.027j0.28) 2.16j22.4Ω阻抗模值约22.5Ω阻抗角约84.5°。距离I段按保护本线路全长的80%~85%整定这是为了保证本线路末端故障时绝对不越限动作。取可靠系数Krel 0.85则ZsetI 0.85×ZL ≈ 1.84j19.04Ω对应约19.1Ω的整定值。I段无时限保护出口后依靠断路器机构动作时间切除故障仿真里一般从保护出口到断路器完全断开还有20~40ms但这个在保护逻辑里不参与计时。距离II段要保护本线路全长同时与相邻线路的I段配合。假设相邻线路也是80km其I段整定值为ZsetI_next 0.85×ZL。II段整定公式为ZsetII ZL Kb×ZsetI_nextKb是分支系数。极端情况下Kb取1算出来ZsetII ≈ 2.16j22.4 1.84j19.04 4.0j41.44Ω模值约41.6Ω。II段时限一般取0.5s比相邻线路I段快同时保证对侧保护先动作切除区外故障。距离III段要躲过最小负荷阻抗同时作为远后备保护。计算最小负荷阻抗需要知道线路最大输送容量和最低运行电压。220kV线路按最大输送能力228MVA算最低运行电压取0.9×220kV198kV折算到相间回路的最小负荷阻抗Zload_min U_line²/S_max ≈ 198²/228 ≈ 172Ω。III段整定值要求ZsetIII ≤ Krel×Zload_minKrel取0.667算出来约114.7Ω实际工程中还要校核与相邻线路II段的配合和灵敏度仿真里我取100Ω时限1.2s。整定结果汇总如下表保护段整定阻抗Ω时限s保护范围备注距离I段19.10本线路约85%瞬时动作距离II段41.60.5本线路全长相邻线路I段范围与相邻I段配合距离III段1001.2本线路及相邻线路远后备躲最小负荷阻抗2.3 双端电源下测量阻抗会被“扭曲”到什么程度这个问题是整个课题最有意思的部分也是很多论文里讲得很模糊的地方。先说过渡电阻的影响。假设故障点存在电弧电阻RF那么保护安装处母线电压UM IM×ZMF IF×RF其中IM是从M侧流入故障点的电流IF是流过故障点总电流。测量阻抗Zm UM/IM ZMF (IF/IM)×RF。在双端电源系统里IF IM INIF一定大于IM所以过渡电阻项被放大了测量阻抗会大于真实故障距离对应的阻抗ZMF。故障点越靠近对侧、对侧电源越强这个放大效应越明显距离I段的保护范围就会缩短。这也是为什么距离保护在带过渡电阻的经过树闪故障里经常灵敏度不足。再来看助增效应。如果在线路中间某个节点有第三电源或环网支路汇入流过保护安装处的电流IM小于流过故障线路区段的电流ITF母线电压为UM IM×ZMT ITF×ZTF测量阻抗Zm ZMT (ITF/IM)×ZTF同样比实际距离对应的阻抗要大。工程上处理这个问题的标准做法就是在II段整定公式里引入分支系数Kb但分支系数的大小要随系统运行方式变化很难取一个固定值仿真研究的价值就是把这个变化规律摸清楚。反向故障也要注意。双端系统里如果故障发生在保护安装处母线背后测量阻抗会落在第三象限或者具有明显的容性特征方向圆判据会在这种情况下自动退出。这正是方向圆特性相比全阻抗圆特性的核心优势。我在仿真里专门设置了反向故障场景做过验证后面结果分析部分会详细讲。3. Simulink模型搭建全流程与关键步骤3.1 系统顶层结构与电源/线路参数设置模型顶层结构从电源到负载依次是M侧三相电源、M侧母线、线路MN分四段PI型线路、N侧母线、N侧三相电源故障模块挂在线上指定节点处。断路器模块加在M侧保护出口位置用于执行保护跳闸信号。电源模块用的是Simscape Electrical库里的Three-Phase Source连接方式Yg中性点接地额定线电压220kV频率50Hz。内阻抗通过短路容量折算M侧短路容量3000MVA折算到220kV侧系统等值电抗X U²/Sk 220²/3000 ≈ 16.1Ω取R X/10 ≈ 1.6Ω所以内阻抗设置为1.6j16.1Ω。N侧短路容量1000MVA折算下来X ≈ 48.4ΩR ≈ 4.8Ω内阻抗设置为4.8j48.4Ω。这样设置的好处是两端电源强度差异明显双端电源对距离保护的影响容易被观察出来。线路用Three-Phase PI Section Line模块整条80km线路拆成四段每段20km。这样做有两个原因一是故障点可以精确设置在20km、40km、60km这些节点位置避免故障点落在线路内部导致故障距离不精确二是多段PI模型比单段PI模型更接近真实线路的分布参数特性。线路单位参数按前面说的设置正序电阻0.027Ω/km正序电感X1/ω 0.28/(2×π×50) ≈ 0.891mH/km正序电容0.014μF/km零序参数按1.8~2.5倍正序值估算。仿真求解器选择变步长ode23tb最大步长限制在50us以内仿真时长0.6s这样既能保证电磁暂态波形精度又不至于让仿真时间太长。3.2 故障场景设置与测量环节设计故障模块用Three-Phase Fault故障类型设置为AB相间短路两相短路这也是相间距离保护最典型的动作对象。故障电阻分两种情况金属性短路设0.01Ω模拟带过渡电阻短路时设10Ω。故障时间控制在0.3s开始、0.35s结束这样0.2s前系统完全稳定0.35s后故障被切除剩余时间可以观察断路器跳闸后的系统恢复波形。测量环节用Three-Phase V-I Measurement模块挂在保护安装处M侧母线出口同时输出三相电压和三相电流的瞬时值。这里有一个重要的单位细节Simscape电气模块输出电压单位是V、电流单位是A幅值可能在几十万伏、几千安的级别后续保护算法需要对瞬时值做标幺化处理否则DFT计算的数值动态范围很大容易在数值运算中丢精度。我的做法是在MATLAB Function里统一除以基准值电压基准是220kV×√2/√3相电压幅值电流基准按系统额定电流幅值取然后再做傅里叶变换。数据进入保护算法前还要经过Zero-Order Hold模块按采样率2kHz即每个工频周期40个采样点把连续信号离散化。这个采样率对全周傅里叶来说足够了50Hz基波加上少量谐波都能覆盖同时计算量也不大。3.3 保护算法实现全周傅里叶、阻抗计算与三段逻辑保护算法的核心放在一个MATLAB Function模块里实现。我选择MATLAB Function而不是纯Simulink模块搭建是因为复数运算、向量缓冲、逻辑状态判断用代码写更加清晰也方便以后移植到不同的保护算法。函数输入是三相电压和三相电流的瞬时采样值各3路信号输出是启动信号、测量阻抗的实部虚部、I段出口、II段出口、III段出口共5个信号。全周傅里叶的实现思路是每个采样时刻把最近一个周期的采样点存进循环缓冲区然后用离散傅里叶公式提取基波相量。公式很直白对N个采样点基波相量X (2/N)×Σ[x(k)×e^(-j×2πk/N)]k从0到N-1。用代码实现时需要注意缓冲区初始阶段的清零和补零问题在前20ms仿真时间内采样点不足一个周期这时候傅里叶结果不可信需要一个启动标志位强制屏蔽输出。阻抗计算部分AB相间短路用线电压Uab Ua - Ub线电流Iab Ia - Ib然后在频域里做复数除法得到Zab的实部R和虚部X。方向圆判据比较|Zm - Zset/2|和|Zset/2|落在圆内且启动信号有效时进入段别判断。三段逻辑就是一个简单的状态机I段动作信号瞬时置位II段信号启动一个0.5s计时器III段信号启动一个1.2s计时器计时器累计到设定值后出口置位。下面是简化后的核心逻辑示意代码function [Zre, Zim, trip1, trip2, trip3] distProtect(uabc, iabc, t) % 简化示意假设内部已完成采样缓冲与DFT相量提取 persistent tStart2 tStart3 trip2Latch trip3Latch % 1. 提取基波相量全周傅里叶代码略 Uab_phasor ua1 - ub1; % 基波线电压相量 Iab_phasor ia1 - ib1; % 基波线电流相量 % 2. 计算测量阻抗 Zm Uab_phasor / Iab_phasor; Zre real(Zm); Zim imag(Zm); % 3. 方向圆判据Zset为当前段整定值 Zset getZset(1); % 按段别切换整定值 if abs(Zm - Zset/2) abs(Zset/2) inZone true; else inZone false; end % 4. 三段时限逻辑 if inZone startFlag trip1 true; % I段瞬时 if isempty(tStart2), tStart2 t; end if t - tStart2 0.5, trip2Latch true; end if isempty(tStart3), tStart3 t; end if t - tStart3 1.2, trip3Latch true; end else trip1 false; end trip2 trip2Latch; trip3 trip3Latch;这里的细节说明一下方向圆判据里的Zset要根据当前是第几段来切换实际工程里三段各有独立阻抗元件并行判断我为了代码简洁用了分段整定。另外II段和III段的计时器用了latch变量避免故障消失后动作信号跌落导致出口丢失这是保护逻辑里的普遍做法。断路器模块用Three-Phase Breaker跳闸信号由I、II、III段动作信号相或后输出。跳闸后输电线路脱离电源故障电流消失故障模块在0.35s自动退出这时候能看到系统恢复正常运行的过程。4. 仿真结果分析与波形解读4.1 区外故障考验选择性的场景第一个仿真正式场景是区外故障把故障点设置在N侧母线之外相当于相邻线路上的故障。也就是故障点不在本线路保护范围内但距离II段和III段可能感受到这个故障。故障设置为AB相间金属性短路故障时间0.3s到0.35s。仿真结果和整定预期一致M侧距离保护测量到的阻抗幅值约为35Ω角度在80°左右。这个测量阻抗超出了I段整定圆19.1Ω所以I段没有动作。但它在II段整定圆41.6Ω以内II段启动后经过0.5s延时在0.38s动作对应故障开始后80ms的出口时间实际是0.3s启动计时0.5s后到0.8s出口这里需要说明实际上故障0.3s开始0.35s故障切除后测量阻抗可能恢复到负荷区II段会返回所以在自动重合闸等配合中II段如果故障已切除就不会动作。这里取决于故障持续时间设置。如果故障是永久性故障0.35s后故障模块虽然退出但断路器未跳闸系统仍然短路这里我要稍微修正场景描述避免逻辑错误。三类典型结果我应该这样设计场景1区外永久性故障故障点位于相邻线路II段范围内。故障0.3s开始持续保持M侧保护II段启动并计时0.5s后0.8sII段出口跳闸。这说明II段作为相邻线路的后备保护起作用时间配合正确。如果把故障时间设置到0.4s就切除例如对侧保护0.35s动作切除则II段在计时过程中检测到阻抗恢复正常而返回不动作。这两种情况在仿真中都能观察到前者验证选择性动作后者验证与相邻保护的配合。场景2反向故障M侧背后母线短路这是考验方向元件是否可靠的关键。故障时流过保护的电流由N侧电源提供测量阻抗落在第三象限负电阻负电抗方向圆判据判断为区外三条段别都不出口。这个场景我做了波形上看电流幅值很大但保护纹丝不动完美验证了方向圆特性的价值。4.2 区内不同位置故障动作速度对比把故障点分别设置在距离M侧20km、40km、68km处全部做AB相间金属性短路故障0.3s开始持续0.2s观察保护动作行为。20km处故障时测量阻抗约Z 20×(0.027j0.28) ≈ 0.54j5.6Ω幅值5.6Ω远小于I段整定值19.1Ω保护瞬时出口。断路器在保护出口约25ms后完全断开仿真波形上电流在0.325s附近跌落到零。40km处故障时测量阻抗约1.08j11.2Ω幅值11.2Ω仍在I段范围内同样I段瞬时动作。68km处故障时测量阻抗约1.84j19.04Ω正好接近I段整定边缘19.1Ω。实际上因为线路分布电容和仿真数值误差测量阻抗可能会在整定圆边界附近波动这时I段可能动作也可能不动作这个现象恰好演示了为什么I段只保护线路的80%~85%留出15%~20%的裕度就是为了应对这种边界模糊。如果在65km处故障81.25%I段可靠动作在72km处故障90%测量阻抗幅值约20.2Ω略超I段整定值I段不动II段在0.5s后动作。把几种故障位置下的结果汇总成一张表故障位置km占线路比例测量阻抗幅值Ω动作段别保护出口时间s2025%5.6I段0.3004050%11.2I段0.3006581.25%18.2I段0.3007290%20.2II段0.80080100%22.5II段0.800这里保护出口时间都是相对于故障起始时刻0.3s而言I段在故障发生的下一个采样周期内出口仿真步长下看起来就是0.300sII段在0.8s出口是因为0.3s启动后需要累计0.5s延时。这个时间对应关系在波形上非常清晰也很适合写进毕业论文的仿真分析章节。4.3 助增与过渡电阻场景测量阻抗的“漂移”这一节重点做两类场景。第一类是过渡电阻的影响在线路50km处设置AB相间短路过渡电阻分别取0.01Ω和10Ω对比测量阻抗变化。金属性短路时测量阻抗约为1.35j14Ω幅值约14.1ΩI段动作。过渡电阻10Ω时由于对侧电源的存在故障点电流IF大于保护安装处电流IM测量阻抗被放大到约25Ω左右已经超过I段19.1Ω的整定值I段拒动保护退避到II段延时动作。这就是双端电源系统带过渡电阻故障时距离保护灵敏度下降的直接证据在单端电源系统里10Ω过渡电阻还不足以让I段完全拒动。第二类是助增效应场景。在线路中点40km处通过一个断路器接入一个负荷支路模拟T接第三电源或重要负荷等效为一个小容量电源和负载并联然后在线路末端60km处设置AB相间短路。对比有负荷支路和无负荷支路两种情况下的测量阻抗。没有负荷支路时测量阻抗是60km线路阻抗约1.62j16.8Ω幅值约16.9ΩI段动作投入负荷支路后由于该支路向故障点提供助增电流保护安装处电流变小而母线电压基本不变测量阻抗增加到约21Ω超出I段整定值I段拒动II段动作。这两个场景的数值在仿真里用Scope和MATLAB绘图都能清晰看到测量阻抗轨迹在R-X平面上会画出一条从负荷区快速移向故障区的曲线。过渡电阻场景下轨迹
220kV双端电源线路相间距离保护Simulink仿真与三段式整定
220kV线路的相间距离保护一旦到了双端电源的供电方式下整定与动作逻辑就比单端电源复杂得多。我前阵子刚好用Simulink做完了这个方向的仿真研究把220kV双端电源系统里的相间距离保护从模型搭建到故障仿真完整跑了一遍收获不少也踩了一些坑。这篇记录就是给准备做这个方向的人一份可以照着走的参考不堆理论公式重点讲清楚双端电源为什么会影响距离保护判断、三段式保护怎么整定、模型在Simulink里怎么搭、故障结果怎么解读。无论你是刚接手的继保专业学生、要交毕设的本科生还是做培训验收的现场工程师照着这个思路走都能得到一套能跑通、能出结果的仿真模型。做这个课题之前有个前提要先说清楚实际220kV线路通常配置纵联保护做主保护距离保护和零序电流保护作为后备保护。我这套仿真把注意力集中在相间距离保护本身刻意不引入纵联逻辑目的是把三段式动作特性、阻抗测量误差、助增电流影响这些底层问题看清楚。方向选这个题目是因为220kV处于区域电网主网架几乎都是双端或环网供电与110kV及以下单端辐射网的保护分析逻辑差异很大很多初学者在这上面栽跟头。1. 项目核心背景与整体思路拆解1.1 220kV双端电源系统为什么是距离保护的“难度天花板”先说结论220kV输电线路里的双端电源会让距离保护的核心变量——测量阻抗——不再直观。距离保护的基本想法很朴素测出保护安装处的电压和电流算出从安装点到故障点的等效阻抗然后和整定值比较落在阻抗圆内就动作。这个思路在单端电源系统里几乎不会出问题因为故障电流方向唯一电压电流的相位关系很清楚。但双端电源一上来情况就变了。第一重变化是故障电流的方向不再唯一。M侧线路故障时流进故障点的电流是M侧电源和N侧电源共同提供的保护安装处测到的电流只是总故障电流的一部分。如果故障处存在过渡电阻故障点电压会被对侧电流抬高或压低测量阻抗就会偏离真实故障距离对应的阻抗值。第二重变化是助增效应如果保护安装点和故障点之间还存在其他电源支路或环网路径流到故障点的电流大于流过保护的电流测量阻抗会系统性偏大灵敏度下降。反过来还有汲出效应测量阻抗偏小可能造成区外误动。第三重变化是方向问题反向故障时电流电压的相位特征和正向故障完全不同没有方向判别的距离保护在双端系统里寸步难行。这些都不是书本上背个公式就能消化的必须靠仿真把故障过程的电压电流波形完整拉出来看。Simulink里可以很方便地调整双端电源容量、线路长度、分支位置、过渡电阻大小每次调整都能在波形上看到测量阻抗的轨迹变化这对理解双端系统的保护行为帮助极大。1.2 为什么选Simulink而不是PSCAD/RTDS做电力系统电磁暂态仿真圈子里常见的选项无非是PSCAD、RTDS、MATLAB/Simulink这几类。PSCAD在电磁暂态计算精度和元件库丰富度上确实有优势RTDS则是硬件在环测试的主流平台。但这个课题选Simulink我是综合考虑了三件事。第一Simulink的Simscape Electrical组件库原SimPowerSystems涵盖了三相电源、三相线路、三相变压器、故障模块、断路器和各种测量元件对220kV单线双端系统这种规模来说完全够用不需要额外许可证。第二也是最重要的一点这个课题不光是搭个电路模型看波形还要在里面实现距离保护算法——全周傅里叶相量提取、阻抗计算、圆特性判据、三段时限配合。这些逻辑用MATLAB Function或S-Function写非常顺手数据和绘图也与Simulink无缝衔接。PSCAD里写自定义元件就要麻烦不少。第三仿真参数扫描、批量跑故障工况、结果后处理MATLAB脚本在这方面几乎是降维打击我可以一次性把几十种故障位置、故障类型、过渡电阻组合全部跑完自动生成对比曲线。当然Simulink也有短板大规模电力系统电磁暂态的计算效率和收敛性不如专用工具但这个课题只有两条线路加两个电源规模很小完全在可控范围内。如果后续有人想往多机系统、行波保护方向扩展再迁移到PSCAD或RTDS也不迟。2. 保护原理与逻辑搭建测量阻抗与三段式整定2.1 测量阻抗是怎么算出来的0度接线与方向圆判据相间距离保护的核心测量量是相间阻抗。以AB相间短路为例保护接入的电压是线电压Uab电流是线电流Iab也就是常说的0度接线方式。计算得到的测量阻抗Zab Uab / Iab这个复数在R-X平面上对应一个点。正常情况下负荷状态下测量阻抗落在负荷区幅值大、角度接近功率因数角故障后测量阻抗迅速跌落到线路阻抗附近幅值小、角度接近线路阻抗角。保护要做的事情就是判断这个测量阻抗点是否落入预先设定的动作区域。动作区域最常用的是方向圆特性。整定阻抗设为Zset圆心在Zset/2处半径为Zset/2动作判据写成复数的模值比较就是|Zm - Zset/2| ≤ |Zset/2|。方向圆的特点是圆过原点也就是说测量阻抗落在原点附近或者第三象限时不动作天然带有方向性正好应对双端电源里反向故障的问题。还有一种偏移圆特性把圆往第三象限偏移一部分适用于母线背后有较长线路的场景但会牺牲一些方向选择性这个课题里我用的是方向圆。整套保护逻辑还要一个启动元件。不能等阻抗都算出来掉进圆里才启动那样在正常运行时就容易频繁计算误动。启动元件可以用相电流突变量、负序电流分量或者低电压启动。我在仿真里用的是相电流突变量启动连续三个采样点的电流变化量超过设定门槛就进入故障处理流程。门槛值一般按最大负荷电流的20%~50%整定太灵敏容易在负荷波动时误动太迟钝又会影响故障切除速度这个要靠多次仿真去调。2.2 三段式整定计算一个80km线路的完整算例这套仿真我采用了一个典型的220kV双端系统参数M侧系统短路容量较大3000MVAN侧系统相对较弱1000MVA线路全长80km。220kV架空线单位正序参数按常见值取R1 0.027Ω/kmX1 0.28Ω/kmC1 0.014μF/km。这样算下来全线路正序阻抗为ZL 80×(0.027j0.28) 2.16j22.4Ω阻抗模值约22.5Ω阻抗角约84.5°。距离I段按保护本线路全长的80%~85%整定这是为了保证本线路末端故障时绝对不越限动作。取可靠系数Krel 0.85则ZsetI 0.85×ZL ≈ 1.84j19.04Ω对应约19.1Ω的整定值。I段无时限保护出口后依靠断路器机构动作时间切除故障仿真里一般从保护出口到断路器完全断开还有20~40ms但这个在保护逻辑里不参与计时。距离II段要保护本线路全长同时与相邻线路的I段配合。假设相邻线路也是80km其I段整定值为ZsetI_next 0.85×ZL。II段整定公式为ZsetII ZL Kb×ZsetI_nextKb是分支系数。极端情况下Kb取1算出来ZsetII ≈ 2.16j22.4 1.84j19.04 4.0j41.44Ω模值约41.6Ω。II段时限一般取0.5s比相邻线路I段快同时保证对侧保护先动作切除区外故障。距离III段要躲过最小负荷阻抗同时作为远后备保护。计算最小负荷阻抗需要知道线路最大输送容量和最低运行电压。220kV线路按最大输送能力228MVA算最低运行电压取0.9×220kV198kV折算到相间回路的最小负荷阻抗Zload_min U_line²/S_max ≈ 198²/228 ≈ 172Ω。III段整定值要求ZsetIII ≤ Krel×Zload_minKrel取0.667算出来约114.7Ω实际工程中还要校核与相邻线路II段的配合和灵敏度仿真里我取100Ω时限1.2s。整定结果汇总如下表保护段整定阻抗Ω时限s保护范围备注距离I段19.10本线路约85%瞬时动作距离II段41.60.5本线路全长相邻线路I段范围与相邻I段配合距离III段1001.2本线路及相邻线路远后备躲最小负荷阻抗2.3 双端电源下测量阻抗会被“扭曲”到什么程度这个问题是整个课题最有意思的部分也是很多论文里讲得很模糊的地方。先说过渡电阻的影响。假设故障点存在电弧电阻RF那么保护安装处母线电压UM IM×ZMF IF×RF其中IM是从M侧流入故障点的电流IF是流过故障点总电流。测量阻抗Zm UM/IM ZMF (IF/IM)×RF。在双端电源系统里IF IM INIF一定大于IM所以过渡电阻项被放大了测量阻抗会大于真实故障距离对应的阻抗ZMF。故障点越靠近对侧、对侧电源越强这个放大效应越明显距离I段的保护范围就会缩短。这也是为什么距离保护在带过渡电阻的经过树闪故障里经常灵敏度不足。再来看助增效应。如果在线路中间某个节点有第三电源或环网支路汇入流过保护安装处的电流IM小于流过故障线路区段的电流ITF母线电压为UM IM×ZMT ITF×ZTF测量阻抗Zm ZMT (ITF/IM)×ZTF同样比实际距离对应的阻抗要大。工程上处理这个问题的标准做法就是在II段整定公式里引入分支系数Kb但分支系数的大小要随系统运行方式变化很难取一个固定值仿真研究的价值就是把这个变化规律摸清楚。反向故障也要注意。双端系统里如果故障发生在保护安装处母线背后测量阻抗会落在第三象限或者具有明显的容性特征方向圆判据会在这种情况下自动退出。这正是方向圆特性相比全阻抗圆特性的核心优势。我在仿真里专门设置了反向故障场景做过验证后面结果分析部分会详细讲。3. Simulink模型搭建全流程与关键步骤3.1 系统顶层结构与电源/线路参数设置模型顶层结构从电源到负载依次是M侧三相电源、M侧母线、线路MN分四段PI型线路、N侧母线、N侧三相电源故障模块挂在线上指定节点处。断路器模块加在M侧保护出口位置用于执行保护跳闸信号。电源模块用的是Simscape Electrical库里的Three-Phase Source连接方式Yg中性点接地额定线电压220kV频率50Hz。内阻抗通过短路容量折算M侧短路容量3000MVA折算到220kV侧系统等值电抗X U²/Sk 220²/3000 ≈ 16.1Ω取R X/10 ≈ 1.6Ω所以内阻抗设置为1.6j16.1Ω。N侧短路容量1000MVA折算下来X ≈ 48.4ΩR ≈ 4.8Ω内阻抗设置为4.8j48.4Ω。这样设置的好处是两端电源强度差异明显双端电源对距离保护的影响容易被观察出来。线路用Three-Phase PI Section Line模块整条80km线路拆成四段每段20km。这样做有两个原因一是故障点可以精确设置在20km、40km、60km这些节点位置避免故障点落在线路内部导致故障距离不精确二是多段PI模型比单段PI模型更接近真实线路的分布参数特性。线路单位参数按前面说的设置正序电阻0.027Ω/km正序电感X1/ω 0.28/(2×π×50) ≈ 0.891mH/km正序电容0.014μF/km零序参数按1.8~2.5倍正序值估算。仿真求解器选择变步长ode23tb最大步长限制在50us以内仿真时长0.6s这样既能保证电磁暂态波形精度又不至于让仿真时间太长。3.2 故障场景设置与测量环节设计故障模块用Three-Phase Fault故障类型设置为AB相间短路两相短路这也是相间距离保护最典型的动作对象。故障电阻分两种情况金属性短路设0.01Ω模拟带过渡电阻短路时设10Ω。故障时间控制在0.3s开始、0.35s结束这样0.2s前系统完全稳定0.35s后故障被切除剩余时间可以观察断路器跳闸后的系统恢复波形。测量环节用Three-Phase V-I Measurement模块挂在保护安装处M侧母线出口同时输出三相电压和三相电流的瞬时值。这里有一个重要的单位细节Simscape电气模块输出电压单位是V、电流单位是A幅值可能在几十万伏、几千安的级别后续保护算法需要对瞬时值做标幺化处理否则DFT计算的数值动态范围很大容易在数值运算中丢精度。我的做法是在MATLAB Function里统一除以基准值电压基准是220kV×√2/√3相电压幅值电流基准按系统额定电流幅值取然后再做傅里叶变换。数据进入保护算法前还要经过Zero-Order Hold模块按采样率2kHz即每个工频周期40个采样点把连续信号离散化。这个采样率对全周傅里叶来说足够了50Hz基波加上少量谐波都能覆盖同时计算量也不大。3.3 保护算法实现全周傅里叶、阻抗计算与三段逻辑保护算法的核心放在一个MATLAB Function模块里实现。我选择MATLAB Function而不是纯Simulink模块搭建是因为复数运算、向量缓冲、逻辑状态判断用代码写更加清晰也方便以后移植到不同的保护算法。函数输入是三相电压和三相电流的瞬时采样值各3路信号输出是启动信号、测量阻抗的实部虚部、I段出口、II段出口、III段出口共5个信号。全周傅里叶的实现思路是每个采样时刻把最近一个周期的采样点存进循环缓冲区然后用离散傅里叶公式提取基波相量。公式很直白对N个采样点基波相量X (2/N)×Σ[x(k)×e^(-j×2πk/N)]k从0到N-1。用代码实现时需要注意缓冲区初始阶段的清零和补零问题在前20ms仿真时间内采样点不足一个周期这时候傅里叶结果不可信需要一个启动标志位强制屏蔽输出。阻抗计算部分AB相间短路用线电压Uab Ua - Ub线电流Iab Ia - Ib然后在频域里做复数除法得到Zab的实部R和虚部X。方向圆判据比较|Zm - Zset/2|和|Zset/2|落在圆内且启动信号有效时进入段别判断。三段逻辑就是一个简单的状态机I段动作信号瞬时置位II段信号启动一个0.5s计时器III段信号启动一个1.2s计时器计时器累计到设定值后出口置位。下面是简化后的核心逻辑示意代码function [Zre, Zim, trip1, trip2, trip3] distProtect(uabc, iabc, t) % 简化示意假设内部已完成采样缓冲与DFT相量提取 persistent tStart2 tStart3 trip2Latch trip3Latch % 1. 提取基波相量全周傅里叶代码略 Uab_phasor ua1 - ub1; % 基波线电压相量 Iab_phasor ia1 - ib1; % 基波线电流相量 % 2. 计算测量阻抗 Zm Uab_phasor / Iab_phasor; Zre real(Zm); Zim imag(Zm); % 3. 方向圆判据Zset为当前段整定值 Zset getZset(1); % 按段别切换整定值 if abs(Zm - Zset/2) abs(Zset/2) inZone true; else inZone false; end % 4. 三段时限逻辑 if inZone startFlag trip1 true; % I段瞬时 if isempty(tStart2), tStart2 t; end if t - tStart2 0.5, trip2Latch true; end if isempty(tStart3), tStart3 t; end if t - tStart3 1.2, trip3Latch true; end else trip1 false; end trip2 trip2Latch; trip3 trip3Latch;这里的细节说明一下方向圆判据里的Zset要根据当前是第几段来切换实际工程里三段各有独立阻抗元件并行判断我为了代码简洁用了分段整定。另外II段和III段的计时器用了latch变量避免故障消失后动作信号跌落导致出口丢失这是保护逻辑里的普遍做法。断路器模块用Three-Phase Breaker跳闸信号由I、II、III段动作信号相或后输出。跳闸后输电线路脱离电源故障电流消失故障模块在0.35s自动退出这时候能看到系统恢复正常运行的过程。4. 仿真结果分析与波形解读4.1 区外故障考验选择性的场景第一个仿真正式场景是区外故障把故障点设置在N侧母线之外相当于相邻线路上的故障。也就是故障点不在本线路保护范围内但距离II段和III段可能感受到这个故障。故障设置为AB相间金属性短路故障时间0.3s到0.35s。仿真结果和整定预期一致M侧距离保护测量到的阻抗幅值约为35Ω角度在80°左右。这个测量阻抗超出了I段整定圆19.1Ω所以I段没有动作。但它在II段整定圆41.6Ω以内II段启动后经过0.5s延时在0.38s动作对应故障开始后80ms的出口时间实际是0.3s启动计时0.5s后到0.8s出口这里需要说明实际上故障0.3s开始0.35s故障切除后测量阻抗可能恢复到负荷区II段会返回所以在自动重合闸等配合中II段如果故障已切除就不会动作。这里取决于故障持续时间设置。如果故障是永久性故障0.35s后故障模块虽然退出但断路器未跳闸系统仍然短路这里我要稍微修正场景描述避免逻辑错误。三类典型结果我应该这样设计场景1区外永久性故障故障点位于相邻线路II段范围内。故障0.3s开始持续保持M侧保护II段启动并计时0.5s后0.8sII段出口跳闸。这说明II段作为相邻线路的后备保护起作用时间配合正确。如果把故障时间设置到0.4s就切除例如对侧保护0.35s动作切除则II段在计时过程中检测到阻抗恢复正常而返回不动作。这两种情况在仿真中都能观察到前者验证选择性动作后者验证与相邻保护的配合。场景2反向故障M侧背后母线短路这是考验方向元件是否可靠的关键。故障时流过保护的电流由N侧电源提供测量阻抗落在第三象限负电阻负电抗方向圆判据判断为区外三条段别都不出口。这个场景我做了波形上看电流幅值很大但保护纹丝不动完美验证了方向圆特性的价值。4.2 区内不同位置故障动作速度对比把故障点分别设置在距离M侧20km、40km、68km处全部做AB相间金属性短路故障0.3s开始持续0.2s观察保护动作行为。20km处故障时测量阻抗约Z 20×(0.027j0.28) ≈ 0.54j5.6Ω幅值5.6Ω远小于I段整定值19.1Ω保护瞬时出口。断路器在保护出口约25ms后完全断开仿真波形上电流在0.325s附近跌落到零。40km处故障时测量阻抗约1.08j11.2Ω幅值11.2Ω仍在I段范围内同样I段瞬时动作。68km处故障时测量阻抗约1.84j19.04Ω正好接近I段整定边缘19.1Ω。实际上因为线路分布电容和仿真数值误差测量阻抗可能会在整定圆边界附近波动这时I段可能动作也可能不动作这个现象恰好演示了为什么I段只保护线路的80%~85%留出15%~20%的裕度就是为了应对这种边界模糊。如果在65km处故障81.25%I段可靠动作在72km处故障90%测量阻抗幅值约20.2Ω略超I段整定值I段不动II段在0.5s后动作。把几种故障位置下的结果汇总成一张表故障位置km占线路比例测量阻抗幅值Ω动作段别保护出口时间s2025%5.6I段0.3004050%11.2I段0.3006581.25%18.2I段0.3007290%20.2II段0.80080100%22.5II段0.800这里保护出口时间都是相对于故障起始时刻0.3s而言I段在故障发生的下一个采样周期内出口仿真步长下看起来就是0.300sII段在0.8s出口是因为0.3s启动后需要累计0.5s延时。这个时间对应关系在波形上非常清晰也很适合写进毕业论文的仿真分析章节。4.3 助增与过渡电阻场景测量阻抗的“漂移”这一节重点做两类场景。第一类是过渡电阻的影响在线路50km处设置AB相间短路过渡电阻分别取0.01Ω和10Ω对比测量阻抗变化。金属性短路时测量阻抗约为1.35j14Ω幅值约14.1ΩI段动作。过渡电阻10Ω时由于对侧电源的存在故障点电流IF大于保护安装处电流IM测量阻抗被放大到约25Ω左右已经超过I段19.1Ω的整定值I段拒动保护退避到II段延时动作。这就是双端电源系统带过渡电阻故障时距离保护灵敏度下降的直接证据在单端电源系统里10Ω过渡电阻还不足以让I段完全拒动。第二类是助增效应场景。在线路中点40km处通过一个断路器接入一个负荷支路模拟T接第三电源或重要负荷等效为一个小容量电源和负载并联然后在线路末端60km处设置AB相间短路。对比有负荷支路和无负荷支路两种情况下的测量阻抗。没有负荷支路时测量阻抗是60km线路阻抗约1.62j16.8Ω幅值约16.9ΩI段动作投入负荷支路后由于该支路向故障点提供助增电流保护安装处电流变小而母线电压基本不变测量阻抗增加到约21Ω超出I段整定值I段拒动II段动作。这两个场景的数值在仿真里用Scope和MATLAB绘图都能清晰看到测量阻抗轨迹在R-X平面上会画出一条从负荷区快速移向故障区的曲线。过渡电阻场景下轨迹 
