企业官网建设流程全解析

最近后台不少兄弟在问220kV电力系统双端电源相间距离保护的Simulink仿真怎么做，尤其是毕业论文做继电保护方向、或者刚进设计院/调试单位想快速上手仿真的朋友。这个题目看起来很“教材”，但真动手搭模型的时候，双端电源的助增效应、测量阻抗计算、三段式配合这些点，每一步都能把人卡到怀疑人生。这篇文章我就以“从零搭出一个能跑、能出波形、能讲清楚动作逻辑”的模型为主线，把参数整定、故障场景设计和结果分析一次讲透。

先明确这个仿真的价值：220kV输电线路是电网骨干网架的主力电压等级，双端电源结构（两侧都有系统电源，负荷潮流可能双向流动）是最常见也最考验保护配合的运行方式。相间距离保护通过测量电压电流算出“测量阻抗”，跟整定阻抗比较来决定动不动——这套逻辑用Simulink做出来之后，你不光能看到保护三段式各自的动作范围，还能直观看到故障位置、过渡电阻、两侧电势相位差对动作边界的影响。适合继保方向学生、青年教师、调试工程师，以及任何想把距离保护从“背公式”变成“看波形”的人。

1. 双端电源系统与相间距离保护的基础认知

1.1 为什么拿双端电源当仿真对象

很多教程或者课程设计都拿单端电源系统来讲距离保护，那个确实简单：短路电流只从本侧电源流出，测量阻抗跟故障距离基本是理想比例关系，整定计算一步到位。但实际220kV线路几乎不可能这么“干净”，绝大多数输电线路两侧都有电源或者电网联络线，短路时故障点电流等于两侧系统共同提供的短路电流之和。

这就带来一个很关键的问题：保护安装处感受到的测量阻抗，不再单纯等于故障距离对应的线路正序阻抗，还会叠加上对侧电源注入电流带来的附加阻抗。举个例子，线路中点发生相间短路时，故障电流一半来自本侧、一半来自对侧，如果两侧系统阻抗参数不一样，两侧电流幅值相位都不一致，测量阻抗就会偏离实际故障距离对应的值——这就是继电保护教材里常说的“助增效应”。单端模型永远复现不了这个现象，而双端电源模型一搭起来，这个偏差会直接反映在动作时间、动作区上，非常有教学和工程参考价值。

另外，双端电源系统还要考虑潮流方向问题。正常运行时有功可能由M侧流向N侧，也可能反过来。相位差变化会改变两侧电流的分配关系，进而影响距离保护的超越和拒动。这些内容放在单端模型里全是“超标内容”，但在双端模型里是绕不开的现实约束。所以我的建议很明确：做220kV距离保护仿真，除非你的课题有特殊限定，否则直接按双端电源来做，信息量和说服力完全不是一个量级。

1.2 相间距离保护三段式到底在保护什么

相间距离保护针对的故障类型是相间短路，包括三相短路、两相短路和两相接地短路。它的基本原理就是：利用线路电压和电流计算出故障点到保护安装处的正序阻抗，再跟预先整定好的阻抗定值比较。保护范围跟电压等级和线路长度没有直接关系，只跟阻抗定值有关，所以天然比电流保护更适合长距离输电线路。

按照常规配置，距离保护分三段（部分场合还有四段）。I段是瞬时动作段，一般按可靠躲过本线路全长的80%到85%整定，因为考虑互感器误差、参数误差和过渡电阻，如果整定到100%，区外故障可能误动，所以留出保护范围之外的“死区”。II段通常按本线路全长加相邻线路的一部分整定，动作时间带0.3到0.6秒延时，用来覆盖I段范围之外的近区故障。III段作为后备保护，整定范围更长，延时也更长，还要保证灵敏度。仿真里如果不把这三段逻辑做全，后面的分析基本等于白做。

这里特别提醒一下：Simulink里的“距离保护模块”如果直接用电力系统库里的Distance Relay（有的版本叫Distance Protection），它内部集成了方向判别和多段逻辑，用起来方便，但不利于讲清楚测量阻抗怎么算。我建议自己搭一个测量阻抗计算子系统，然后跟整定值比较，最后接多段延时逻辑，这样每个环节都能看到中间量，调试起来也直观。

1.3 仿真模型整体架构

整个模型从拓扑上分五块：M侧电源、N侧电源、输电线路（用两段PI模型串联，方便设置不同故障位置）、三相故障模块、保护与断路器闭环。保护部分又细分为三相电压电流测量、相间阻抗计算、方向判别、三段距离逻辑、断路器跳闸信号。

串起来就是：系统正常稳态运行到0.1秒时，故障模块投入，线路某一点发生AB相间短路，保护安装处的电压电流发生突变，阻抗计算模块实时算出测量阻抗Zmeas，三段逻辑模块判断Zmeas是否落在动作区内，满足条件后输出跳闸信号给三相断路器，断路器在下一个过零点（或者约一个半周波后）分闸，切除故障线路。这样一个闭环模型跑完之后，你既能看故障电流电压波形，也能拉出测量阻抗轨迹，还能对比不同保护段动作的先后时间。

2. 220kV系统建模：先把“电网底子”搭对

2.1 双侧电源等值参数

220kV系统的电源在Simulink里一般用Three-Phase Source模块来等效，内部是理想电压源串联等值阻抗。等值阻抗怎么给？工程上最常用的方式是通过母线短路容量反推。比如M侧母线三相短路容量为3000MVA，那么等值阻抗模值为：

|ZsM| = Un² / Ss = 220² / 3000 ≈ 16.13Ω

220kV母线的X/R比通常按10到15考虑，我取X/R=9，这样实部虚部都有，更接近工程实际。于是M侧等值阻抗可以取：

ZsM = 1.731 + j15.55Ω

N侧短路容量较小，取1800MVA，算出来模值约26.89Ω，取X/R=9，得到：

ZsN = 2.885 + j26.34Ω

两侧电源电压幅值都设为220kV（线电压有效值），频率50Hz，M侧初始相位设0度，N侧相位根据潮流方向设定。正常情况下可设N侧滞后15度，模拟功率从M侧送往N侧的运行方式。如果你想研究潮流反向的情况，直接把N侧相位改成超前M侧15度，其他参数不用动，就可以对比同一套保护逻辑在双向潮流下的动作差异——这是双端模型独有的“加分项”。

2.2 线路参数与模型选择

输电线路模型在Simulink里主要有两个选择：Three-Phase PI Section Line（集中参数PI模型）和Distributed Parameters Line（分布参数模型）。220kV线路长度100km以内时，PI模型精度足够了，而且建模简单、仿真速度快，我建议直接用PI模型。不过要注意一点：100km线路用一个PI段模拟，在高频暂态下可能不够精准，稳妥做法是把线路拆成两段或三段，每段50km左右。这样还有一个额外好处——你可以在任何一段的末端设置故障，方便模拟不同的故障发生位置。

线路参数按典型220kV架空线取值，正序电阻每公里0.027Ω，正序电感每公里1.02mH，正序电容每公里0.0113μF。那么对应到串联阻抗和并联电容：

• 正序电阻：0.027Ω/km
• 正序电抗：ωL = 2π×50×1.02e-3 ≈ 0.32Ω/km
• 正序容抗：1/(ωC) = 1/(2π×50×1.13e-8) ≈ 2817Ω·km（换算成并联电容导纳时要小心单位，Simulink面板直接填电容值μF/km即可）

零序参数由于架空地线和大地返回效应的存在，电阻和电抗都会大一些，典型值可取正序的3倍左右，例如零序电阻0.1946Ω/km、零序电感3.06mH/km、零序电容0.0078μF/km。我这里给的都是有名值，在Simulink里直接填就行，不用换算标幺值。

线路模型面板中Sample time建议填0，表示连续采样。系统频率填50Hz。如果后续做暂态分析，建议把“Measurements”选项里的电压电流测量打勾，方便直接把线路端口的瞬时值引出到示波器。

2.3 故障与测量模块配置

故障模块用Three-Phase Fault，这是仿真里最关键的动态事件源。针对相间距离保护的研究对象，我建议做AB相间短路（不接地）和三相短路两个典型场景。在故障模块里，Fault Type选Phase-to-Phase，对应故障相选A和B，可选的“Ground Fault”不要勾选，这样就是纯相间短路。

故障时刻的设置要看系统是否进入稳态。电源和线路都是线性元件，启动后大概0.02到0.04秒就能进入稳定正弦状态，所以故障投入时刻设在0.1秒完全够稳。Transitional State设置成“Fault then Open”，第一转换时间0.1s，第二转换时间0.2s，这样故障持续0.1秒，保护如果动作，断路器在0.2秒前应该完成分闸。

故障电阻默认是0.001Ω，模拟金属性短路。如果你想深入研究过渡电阻对保护的影响，把这个值改成5Ω、10Ω、25Ω分别跑几次就能得到一组对比数据。测量环节用Three-Phase V-I Measurement模块接在保护安装处（M侧断路器外侧、线路入口母线），注意它的电压测量模式选“Phase-to-ground”，电流测量直接变送到保护子系统。电压等级和变比问题在仿真里可以先不折算，直接按一次侧有名值参与计算，后面讲到测量阻抗时统一处理。

3. 距离保护逻辑实现与整定计算

3.1 测量阻抗计算模块

相间距离保护的核心是测量阻抗，ABB、南瑞等厂家装置的通用做法是用故障相的线电压和线电流差来计算。以AB相间故障为例，保护感受到的相间阻抗为：

Zab = Uab / (Ia - Ib)

这里Uab是A相母线电压减去B相母线电压，Iab = Ia - Ib是两相电流之差。对AB相间短路来说，这个比值约等于从保护安装处到故障点的正序阻抗；对其他类型故障，只要不是本相间故障，测出来的阻抗不会落入整定区，保护天然有选相能力。这是相间距离保护跟接地距离保护最大的区别——接地距离还需要零序电流补偿，相间距离不需要。

Simulink里如何实现？从V-I Measurement拉出三相电压电流后，把A相和B相电压做减法得到Uab，A相和B相电流做减法得到Iab，然后送入一个Complex to Magnitude-Angle模块之前，先确认数据格式。更直接的办法是用“Phasor”类型信号路径，在Powergui里把仿真类型设为“Phasor”，然后直接用复数除。但Phasor模式做不了故障暂态细节，我建议还是用连续模式，把瞬时值通过PLL或傅里叶变换提取基波相量后再算阻抗。最简单的做法是先用“Fourier”模块提取50Hz基波的有名值幅值和相位，再把Uab和Iab的复数比值算出来。

如果你不想搭太复杂的相量提取，可以用一个折中方案：把三相电压电流瞬时值分别做增益处理，然后进“Three-Phase Sequence Analyzer”模块，输出正序和相间分量。不过要提醒，Sequence Analyzer在故障瞬间会有一个周波左右的动态响应，结果虽然精确但波形前沿有一点延迟，分析稳态测量阻抗是够用的。我自己习惯的做法是保留原始瞬时值波形，同时另拉一路经过Fourier模块的基波量，两组对照看，既能看清暂态冲击，又能读稳态阻抗。

3.2 三段式整定值计算

整定计算是整个仿真里数字最好看也最容易出错的部分。以100km线路为例，已知线路单位正序阻抗z1 = 0.027 + j0.32 Ω/km，全线路总阻抗：

ZL = z1 × 100 = 2.7 + j32.0Ω

模值约32.11Ω，阻抗角约85.2度，这个角度决定了方向元件最大灵敏角。

距离I段按可靠躲过本线路末端，系数取0.8到0.85，这里用0.85：

ZI = 0.85 × (2.7 + j32.0) = 2.295 + j27.2Ω

也就是说，测量阻抗落入以该值末端的圆内或四边形内时，I段瞬时动作，延时设为0.02s（断路器固有分闸时间约两个周波）。

距离II段要覆盖本线路全长，还要延伸到相邻线路一部分。考虑对侧N母线到相邻线路的配合，若相邻线路阻抗也按同样参数取50km，则其阻抗为1.35 + j16.0Ω。II段定值按本线路全长加上相邻线路的50%：

ZII = ZL + 0.5 × (1.35 + j16.0) = 2.7 + j32.0 + 0.675 + j8.0 = 3.375 + j40.0Ω

延时一般取0.5s，用来配合相邻线路I段，避免区外故障误动。

距离III段按躲过最小负荷阻抗整定。假设线路最大输送功率300MW，功率因数0.9，母线最低运行电压0.95Un，那么最小负荷阻抗模值约为：

|Zload| = (0.95 × 220)² / 300 ≈ 145.6Ω

负荷的阻抗角对应功率因数角arccos0.9≈25.8度。III段定值需要躲过负荷阻抗，并按1.2到1.3的可靠系数缩小，同时满足灵敏度系数≥1.5，所以可以取：

ZIII ≈ 100Ω

方向性不进入III段反向区，仍然按正方向灵敏角85度来整定。

看明白没有？整定过程本质上是“定值躲故障、躲负荷、留灵敏度”的三角平衡。I段要短，II段要跟相邻配合，III段要长但不能碰到正常负荷。仿真里把这些定值分别做成常数模块，通过比较器跟测量阻抗幅值做比较，再叠加上相应的时间延迟，就能得到三段跳闸信号。

3.3 双端电源下的助增效应与超越问题

这一步是双端电源建模的精髓，也是整定计算里最值得展开讲的内容。当两侧都有电源时，故障点的总短路电流I_f = I_M + I_N，其中I_M是本侧保护感受到的电流，I_N是对侧注入的电流。考虑一个很小的过渡电阻Rg，相间故障下保护安装处的测量阻抗可以写成：

Z_meas = m × ZL + (1 + I_N / I_M) × Rg

后面这个(1 + I_N / I_M)就是助增系数的简化表达。如果I_N和I_M同相，助增系数大于1，过渡电阻造成的附加阻抗被放大，测量阻抗比实际故障点对应的线路阻抗更大，保护倾向于“拒动”，容易缩保护范围。如果两侧电流相位差拉开，I_N / I_M出现负实部，附加阻抗可能为复数甚至减小，测量阻抗有可能落入整定区以外更深的区域——这就是“超越”。

在仿真里怎么观察？把过渡电阻从0.001Ω改成10Ω和25Ω分别跑故障，同一故障位置下测量阻抗的末端会沿着“附加阻抗方向”偏移。这个偏移方向完全由对侧电流与本侧电流的相位比决定。如果你把N侧电源相位设成滞后M侧25度，两侧短路电流基本同相，助增系数大约在1.3到1.6之间；如果你把N侧调成超前M侧20度，对侧电流在某些故障点下可能出现相位反转，测量阻抗末端就往反方向偏移。

这个现象对保护整定的启发是：双端电源系统的II段和III段整定不能简单沿用单端比例关系，必要时要根据系统最大/最小运行方式验算助增系数，并校核末端故障时的灵敏度。仿真研究到这一步，才算真正把“双端电源”这个题做出来了，而不是仅仅把两个电源接在线路两端当摆设。

3.4 断路器跳闸与逻辑闭锁

保护动作信号要真正切除故障，还需要跟断路器模块配合。Simulink里三相断路器用Three-Phase Breaker，控制端输入一个外部逻辑信号。把三段保护的输出经过“逻辑或”以后，送到断路器的控制端口即可。要注意断路器模块本身的内部电阻在断开前后是阶跃变化的，默认断态电阻1e6Ω，合态电阻0.001Ω，这个设置影响不大，保持默认就行。如果想模拟实际断路器分闸时刻，可以加一个0.02到0.04秒的延迟，或者用一个PLL检测电流过零点再分闸。

还需要加一个“跳闸保持”逻辑。实际保护装置一旦发出跳闸命令，这个命令是自保持的，不能因为测量阻抗恢复就自动复归断路器，否则故障未切除时测量阻抗波动会导致断路器反复分合。Simulink里可以用一个RS触发器，把三段动作信号作为Set端，故障消失后用一个延迟脉冲作为Reset端，或者干脆用一个Memory模块保持最后的跳闸状态，仿真结果更贴近实际。

4. 仿真场景设置与结果分析

4.1 不同故障位置的动作行为分析

模型搭好之后，最重要的实验就是“故障位置扫描”。把线路分成两段PI模型之后，故障点可以方便设置在20%、50%、80%、以及超过本线路范围的相邻线路段。比如在两段50km线路之间设置故障，对应50%位置；在第二段线路末端设置故障，对应100%位置；再往后加一段模拟相邻线路的PI段，就可以设置120%位置的区外故障。

我实际跑下来的一组典型结果：AB相间金属性短路，故障点在20%位置时，测量阻抗约0.54 + j6.4Ω，远小于I段定值2.295 + j27.2Ω，I段瞬时动作，0.02s后断路器跳闸；故障点在50%位置时，测量阻抗约1.35 + j16.0Ω，依然在I段区内，动作时间基本一致；故障点在85%位置时，测量阻抗2.295 + j27.2Ω落在I段定值边界附近，由于整定延时和阻抗轨迹暂态，可能出现I段拒动，但II段在0.5s延时后可靠动作；故障点在120%相邻线路时，I段、II段都不会动作，III段经过1.0s延时动作切除——当然这里的前提是III段范围覆盖了相邻线路。

把这一组结果整理成表格非常直观，论文里或者技术报告里直接能用：

看这个表格，你会发现距离保护的分段逻辑天然形成了“近处快、远处慢、区外还能后备”的时间阶梯，这个特点就是电网继电保护配合设计的核心框架。

4.2 过渡电阻影响量化对比

金属性短路距离保护表现几乎完美，换上有过渡电阻的真实故障（比如电弧电阻、杆塔接地电阻派生的相间电阻），表现就更有意思了。这里强调一下：相间短路过渡电阻通常比接地短路小，但雷击闪络、树枝碰线、弧光电阻都可能达到几欧到几十欧。

把同一个50%位置的AB相间故障，故障电阻分别设为0Ω、10Ω、25Ω，三次仿真下来观察M侧保护测量阻抗。0Ω时测量阻抗为1.35 + j16.0Ω，刚好在I段区内；10Ω时，因为助增系数约1.4，附加阻抗约14Ω，但注意这个附加阻抗本身带方向，跟线路阻抗角不一致，所以测量阻抗模值增加到约24Ω，角度也可能偏移到60多度，结果仍然在I段内但裕度变小；25Ω时，附加阻抗约35Ω，测量阻抗模值可能超过42Ω，I段直接判为区外，只能等II段0.5s延时动作。

过渡电阻对距离保护的影响本质是“把阻抗末端往某个方向推”，而不是简单的“变大变小”，方向由助增电流相位决定。如果两侧电流相位接近，阻抗扩大，保护范围收缩；如果两侧电流相位反向，阻抗缩小，保护范围扩大甚至超越。这个结论用仿真结果画成阻抗平面图会非常惊艳，我用MATLAB脚本把测量阻抗复数画到R-X平面上，配合三段动作边界圆，一眼就能看出过渡电阻让动作点怎么跑。

4.3 双端电源相位差敏感性分析

除了过渡电阻，双端电源两侧电势相位差也是影响测量阻抗的重要变量。用三个工况对比：N侧滞后15度、同相位、N侧超前15度。实际工程中这两个角度对应的是系统潮流方向：滞后对应M侧送功率出去，超前对应N侧倒送。

仿真结果很典型：N侧滞后15度时，两侧故障电流接近同相，助增系数大约1.2到1.5，测量阻抗略微偏大，保护范围略有收缩，动作时间基本不受影响；N侧超前15度时，对侧电流相位反转，测量阻抗里的附加阻抗变成负方向，极限情况下在末端故障时测量阻抗可能比实际距离还小，II段动作范围外扩，可能出现I段误动风险。这个现象在教科书里叫“超越”或“外伸”，双端仿真里还原得特别清楚。

如果你想把结果做得再深入一点，可以扫多个相位差，比如从-30度到+30度每隔5度取一个值，记录故障点在线路末端时的测量阻抗模值，画一条“阻抗-相角差”曲线。这样能直接看出哪个运行方式下保护最危险，也能验证整定值在最恶劣工况下是否还有选择性。

5. 踩坑实录：常见问题与排查思路

5.1 仿真发散与振荡

Simulink电力系统仿真最常见的错误就是“发散”，表现是波形变成一片噪声或者数值直接爆表到1e20以上。220kV这个量级比较容易踩这个雷，主要原因有三个：一是电力电子开关模块的阻抗阶跃太大，二是仿真算法和步长不匹配，三是Powergui的求解器配置错误。

解决思路按优先级来：在Simulink配置里把求解器选成ode23tb或者ode15s（刚性求解器），最大步长限制在0.0001秒左右，别默认auto。第二，所有断路器和故障模块如果要模拟突然投切，建议把Transitional State用起来，而不是靠断开开关瞬间的高阻变化。第三，如果是Phasor模式下发散，检查你有没有给电源设置Internal Resistance，纯电压源直接接线路确实容易产生数值振荡。

5.2 测量阻抗突变到负值或者跳变

保护逻辑里如果发现测量阻抗在故障瞬间突然变成负值或者在地图上乱跳，十有八九是电流相位提取环节的问题。相间阻抗的计算依赖精确的基波相量，在故障瞬间电流中含有大量衰减直流分量和高次谐波，直接用理想变压器模块的瞬时值做除法，噪声无限大。

我推荐的做法是加一个二阶低通滤波器或者傅里叶变换模块，把基波分量提取出来之后再算幅值和角度。滤波时间常数大约一个周波，故障后0.02到0.05秒内测量阻抗会逐渐稳定到正确值，这个过渡过程本身就是保护动作延时的物理原因，不算问题。另外一个细节：AB相间故障时，要严格用Uab和Ia-Ib计算，不要用三相平均电流；如果用相电流，会有选相元件的问题，阻抗算出来也是错的。

5.3 保护动作了但断路器不跳闸

这个现象我调试那会儿被坑过好几回。逻辑上三段保护输出正常，波形该跳都跳了，但断路器就是不动作，或者动作特别慢。排查顺序先看Three-Phase Breaker模块的“External control”勾没勾，没勾的话外部信号根本进不去；再看断路器控制端的信号类型，Simulink的断路器控制端需要的是开关信号，不是脉冲，最好用0/1常数经过Goto/From或者Data Type Conversion模块再进去。

还有一个容易忽略的问题：断路器模块在断开瞬间，如果仿真算法是变步长，可能因为开关点检测粒度太粗，导致分闸时刻偏差十几毫秒。解决办法是打开Powergui里的“Enable synchronization of Breaker and Fault”选项，让开关动作跟电流过零点对齐，这样断路器电流波形不会出现数值毛刺，跳闸时序也更真实。

5.4 初始化与稳态收敛问题

很多同学搭好模型后一运行，0秒到0.1秒的“正常运行时”波形还算正常，但故障发生瞬间整个系统电流电压突然跳到离谱值。这个可能不是保护逻辑的问题，而是系统在故障前就没有完全进入稳态。220kV输电线有非常大的分布电容，启动瞬间会有一个充电冲击过程，如果仿真时间不够长，到0.1秒故障发生时电容电荷还没稳定。

解决方法是把故障发生时间往后推，一般设到0.15甚至0.2秒，让线路电容充电过程至少持续5个周波以上。也可以先用Powergui的“Load Flow”功能初始化线路和电源，再做动态仿真，这样初始条件更接近稳定运行点。我自己习惯是最开始设0.05秒过渡，0.1秒故障，观察前0.1秒波形确认稳态之后再去分析故障段，而不是一上来就读取全部数据。

整套仿真做下来，我个人体会最深的一点是：Simulink模型跑通只是第一步，真正有价值的是后面那组“位置-阻抗-动作时间”的对照实验。距离保护不是简单的一个阻抗定值比较，而是系统阻抗、过渡电阻、对侧助增、时间配合交织在一起形成的复杂行为。建议你把每个工况下的测量阻抗轨迹都导出来，画在R-X平面上，再叠加上三段动作区边界，这个图一出来，整篇论文的技术含量直接拉满。后续想扩展的话，可以在这个模型上加振荡闭锁、重合闸，甚至把线路换成串补线路，研究串补电容对距离保护的影响，每一步都是可以继续做下去的课题方向。