做电力系统继电保护或者故障分析这块，两相短路故障的仿真模型是一个绕不开的基础工具。不管是做保护定值整定、故障反演，还是给课程设计交差，大家都得在Matlab/Simulink里把这套东西搭起来过一遍。两相短路属于不对称故障，它的特征、计算方法、保护算法的逻辑，都和三相短路、单相接地完全不同，如果只是拿现成模型跑个结果，很难吃透里面的门道，一旦自己从零搭就抓瞎。
这篇博文我用自己的实操经验，把两相短路故障仿真模型从原理、模块搭建、参数设置到波形解读、问题排查整个流程重新走一遍。模型采用电力系统常用的分布式参数输电线路结构，故障模块模拟A、B相金属性短路，配套电压电流测量和序分量提取。适合正在做电力系统仿真课设、搞继电保护算法验证、或者想把不对称故障机理弄明白的同行参考，新手也可以按步骤复现。
1. 两相短路故障的原理与建模前准备
1.1 两相短路时电气量变化的三个关键特征
两相短路（通常写作两相相间短路）指的是电力系统中两个不同相别的导体发生短接，最常见的是A相与B相短路、B相与C相短路这类形式。系统发生两相短路后，电气量变化有几个非常明显的特征，建模前先把这个底子打扎实，后面看波形才心里有数。
第一个特征是故障相电流显著增大且相位相反。假设A相和B相发生金属性短路，两个故障相流过的短路电流幅值相等，相位相差180度，C相电流基本不受影响。这里的“金属性短路”是指故障点过渡电阻近似为零，理想情况下短路点的电压为零，实际仿真中过渡电阻可以设成0.001欧姆这种特别小的值来模拟。
IA = -IB IC = 0（近似）第二个特征是没有零序分量。两相短路不接地，短路回路中没有大地参与，所以零序电流只有在变压器中性点接地、并且存在接地通路时才可能出现。而两相短路发生在相与相之间，完全没有和地的连接，零序分量在理论上是零。这一点非常重要，直接在保护方案上和单相接地故障区分开——接地故障保护要大量依赖零序分量，而相间短路保护只能靠相电流变化和相间电压变化来识别。
第三个特征是故障相的相电压降低。以A、B两相短路为例，故障点处A相和B相的电压近似相等（如果短路阻抗为零则完全相等），非故障相C相电压变化相对较小。在保护安装处测量时，线电压UAB下降明显，UBC和UCA也都会有不同程度的变化，具体大小取决于故障位置距离电源的远近。离电源越近，故障相电压跌落越厉害，短路电流也越大。
1.2 对称分量法视角下的两相短路
分析不对称故障的数学工具是对称分量法，把三相不对称的电流电压分解成正序、负序、零序三组对称分量。两相短路的序分量特征和单相接地短路有本质区别。
以A、B两相短路为例，选基准相为A相（特殊相），边界条件可以写成：
IA + IB + IC = 0（无零序电流） IA = -IB（两故障相电流等量反向）对应到序分量上：
I0 = 0 I1 = -I2（正序电流和负序电流等量反向）也就是说，两相短路时零序电流严格为零，正序电流和负序电流大小相等、方向相反。故障电流的幅值用序分量网络计算时，等于正序电流分量的根号3倍。
在Simulink模型里，虽然可视化界面上看到的是三相瞬时值，但底层的计算逻辑和波形分析的逻辑都离不开序分量的概念。所以在模型里我通常同时搭建一个序分量提取模块，这样不仅可以看三相波形，还能直接看到正序、负序、零序的数值变化，在做保护逻辑验证时特别有用。
2. 仿真模型整体架构与关键模块实现
2.1 模型拓扑结构设计
我搭建的两相短路仿真模型采用典型的“电源-线路-故障点-负载”结构。这个拓扑是电力系统仿真最基础的结构，既能覆盖故障分析的主要场景，又不会因为系统过于复杂导致计算缓慢或者难以调试。
主电路从三相电压源开始，经过一段输电线路到达故障点，故障点之后再接一段线路和三相负载。电源侧采用等效电压源，内部设置三相电压幅值和相位；故障点用Simulink库里的Three-Phase Fault模块实现，这个模块是Simscape Electrical（原SimPowerSystems）库里的核心故障元件，可以灵活配置多种故障类型和故障时序。
模型整体连接关系如下：
三相电源 -> 三相电压电流测量模块 -> 输电线路（分布式参数） -> 故障模块 -> 输电线路 -> 三相负载这里有一点需要专门提醒：电压电流测量模块一定要放在故障模块之前，这样采集到的才是故障发生前的正常运行量和故障发生后的短路电流。如果把测量模块放在故障模块之后，你看到的就只是一部分电流，而不是保护安装处感受到的完整故障电流。
2.2 输电线路模型：分布参数和集中参数的取舍
Simulink的输电线路模型分为集中参数（pi型）和分布参数（Bergeron型）两种，初学时经常会忽略两者区别，直接用了默认的pi型模型，这样做小型模型问题不大，但对于线路长度超过几十公里，或者需要研究行波特性、故障暂态过程的时候，pi型模型的精度就明显不够了。
我做两相短路模型时选择了分布参数模型。分布参数模型考虑线路电感和电容沿线路均匀分布的特性，用波过程的观点来处理，能够比较真实地反映故障瞬间电压电流行波在线路上的传播和折反射现象。更关键的是，在继电保护算法验证中，提取故障特征量（比如暂态分量、谐波分量）时，分布参数模型的波形更接近实际物理过程。
分布参数模型的参数需要设置单位长度电阻R、电感L和电容C，以及线路长度Length。我在模型里使用了一个常用的110kV架空线路参数：
单位长度正序电阻： 0.01273 Ω/km 单位长度正序电感： 0.9337e-3 H/km（约0.293 Ω/km） 单位长度正序电容： 12.74e-9 F/km（约0.401e-6 S/km） 线路长度： 100 km正序阻抗Z1和零序阻抗Z0如果有实测数据就分别设置，没有实测数据时，工程上经常近似取Z0 = 3×Z1来估算。对于我的模型来说，因为验证的是两相短路这种故障类型，零序阻抗不影响故障电流计算结果（零序电流为零），所以零序参数只需要随便填一个合理的值即可，不会影响波形结果。
2.3 电源和负载参数整定思路
三相电压源的电压等级我选了110kV，频率50Hz。这个电压等级是国内高压电网最常见的电压等级，用这个参数做出来的仿真结果有较好的参考意义。电源内部的电阻和电感采用默认的0.002pu和0.08pu（标幺值），这些值模拟的是电源内部的等效阻抗，如果偏离实际太远，会直接影响短路电流的大小。
负载参数需要结合线路传输功率来设置。我设置了额定有功功率100MW、无功功率10Mvar，功率因数约为0.995，模拟高压输电线路末端的重载工况。这里要注意的是，负载阻抗的大小直接决定了故障前的稳态电流水平，如果负载阻抗太小、电流太大，故障电流的相对变化就不明显；如果负载阻抗太大、正常电流太小，故障电流的绝对数值可能很小，看起来不明显。合理的负载设置是让额定电流在几百安培的级别，这样故障电流能到几千安培，波形变化就非常直观了。
3. 核心参数计算与仿真环境配置
3.1 短路电流的理论计算
动手搭模型之前，先把理论值算出来，后面和仿真结果对比的时候才有依据。以A、B相两相短路为例，故障点远离电源时，忽略线路电容的影响，短路电流有效值的简化计算公式为：
I_k2 = (E / (Z1 + Z2)) × √3其中E是故障点正常时的相电压（110kV系统取相电压约63.5kV），Z1和Z2分别是系统到故障点的正序和负序等值阻抗。
把输电线路阻抗折算一下。线路长100km，单位长度正序阻抗0.01273 + j0.293欧姆每公里，那么总正序阻抗为：
Z1 = 0.01273 × 100 + j0.293 × 100 = 1.273 + j29.3 Ω ≈ 29.33∠87.5° Ω两相短路计算时，负序阻抗和正序阻抗近似相等（静止元件满足Z2 = Z1），所以：
I_k2 = 63.5kV × √3 / (2 × 29.33) = 110kV / 58.66 ≈ 1.875kA额定电流大约是100MW / (√3 × 110kV) ≈ 525A，那么两相短路电流大约是额定电流的3.57倍。这个倍数关系很重要，后面看波形的时候可以直接验证：故障相电流峰值大约从正常的740A左右（峰值，1.414×525A）跳到2.65kA左右（峰值，1.414×1.875kA）。注意这里没有考虑电源内阻抗，实际模型里如果电源内阻设置得比较大，短路电流会略小一点。
3.2 故障模块配置方法与参数详解
Three-Phase Fault模块是整个模型的核心。双击模块进入参数配置界面后，有几个地方容易被忽略：
故障类型选择。在Fault between选项中，要选择Phase A和Phase B，对应的就是A相和B相两相短路。这里要特别小心，如果误选了Phase A、Phase B和Phase C三个选项，那就变成三相短路了，整个模型的性质就完全变了。
过渡电阻设置。在Fault Resistance里填0.001欧姆，模拟金属性短路。不建议填0，因为数值0在仿真算法中有时会造成数值计算奇异，而0.001欧姆对结果的影响微乎其微，却能显著提升仿真收敛性。
故障时序设置。在External control of fault timing这个选项里，如果你勾选了外部控制，故障的触发时间需要额外的逻辑信号来给定；如果不勾选，直接在Switch times输入框中写时间数组即可。我习惯用内部时序开关，写为[0.1 0.3]，表示0.1秒时故障开始，0.3秒时故障切除，故障持续0.2秒。
这里解释一下为什么故障持续时间选0.2秒。0.1秒对应的是系统稳定运行了5个周波（50Hz系统一个周波是20ms，5个周波正好0.1秒），此时暂态过程已经基本衰减完毕，故障后的稳态波形清晰稳定；0.2秒的故障持续时间足够观察暂态冲击和稳态短路电流，又不至于让模型运行时间过长。
3.3 仿真求解器配置与采样参数选择
Simulink模型的仿真设置在这个模型里体现为精度和速度的平衡。我的配置如下：
求解器类型：变步长（Variable-step） 求解器算法：ode23tb 最大步长： 1e-4 秒（0.1ms） 最小步长： 自动 相对误差： 1e-3 仿真时间： 0.5 秒这里每一步都是调试出来的经验值。变步长配合ode23tb算法是Simscape Electrical电气系统仿真的常用组合，因为电气系统存在各种时间常数差异极大的元件——电感的暂态时间常数可能只有几毫秒，而仿真总时间要看到0.5秒（25个周波），固定步长要么太慢要么精度不够。ode23tb对刚性系统处理能力强，适合这种含有电感和高次谐波的电路模型。
最大步长设置0.1ms是为了避免步长跨越故障切换时刻。如果最大步长太大，比如默认的0.01秒，仿真步长在0.1秒附近跨过故障发生的边界时，有可能出现严重的数值误差，甚至导致波形出现奇怪的毛刺。这一点也是我在反复调试中明确感受到的。
采样和录波环节，我用了Three-Phase V-I Measurement模块，它可以直接输出三相电压电流信号到Scope或者工作空间。这里有一个使用技巧：测量模块输出的信号格式是Simulink的复合信号（Complex），如果直接接入普通的Scope端口显示，波形会是虚部和实部混合的复杂数据。我一般在测量模块输出的位置加一个Simulink信号类型转换，或者利用Three-Phase V-I Measurement模块内置的输出模式，选择输出相电压的幅值实际值和相位，这取决于你后续是看瞬时波形还是看相量结果。
4. 仿真结果深度解读与对比分析
4.1 故障波形特征解读
把故障前、故障中、故障切除后的波形完整跑一遍后，Scope上能看到几个非常典型的阶段。
0到0.1秒是正常运行阶段，三相电流为标准正弦波，幅值一致、相位相差120度。A相电流峰值约740A，与计算值吻合。
0.1秒到0.3秒是故障阶段，这是分析的核心。A相电流和B相电流几乎瞬间增大到两千安培级别的峰值，两个波形呈反相关系——A相电流达到正峰值时，B相电流达到负峰值；C相电流基本维持故障前的正弦波形。这就是两相短路最直观的特征波形图。
电压波形方面，故障瞬间A相电压和B相电压同时大幅跌落，但两者并不完全为零，因为故障点距离电源还有100km的线路阻抗，线路上存在压降，保护安装处的电压是故障点电压加上线路压降的结果。这一点和“故障点电压为零”的教科书假设有差异，实际仿真波形中表现出来的是电压大幅下降但残压不低。残压的具体数值取决于系统阻抗和线路阻抗的比值，这也解释了为什么距离保护可以靠测量阻抗的变化来判断故障位置。
0.3秒故障切除后，电流和电压经过一个短暂的暂态振荡过程恢复到正常运行值。这个暂态振荡有几个周波的振荡过程，幅值逐渐衰减，是线路储能元件（电感和电容）释放能量的自然响应，属于正常现象。如果振铃时间太长或者不衰减，就要检查是不是线路参数设置过于理想（零阻尼）导致的。
4.2 两相短路vs三相短路vs单相接地：仿真结果对比
同一个模型基础上，把故障模块的选相改成三相或者单相接地，对比三种故障的波形特征，理解会更立体。
| 对比维度 | 两相短路（A-B） | 三相短路 | 单相接地（A相） |
|---|---|---|---|
| 故障电流大小 | I_k2 = (√3/2) × I_k3，约为三相短路的0.866倍 | 最大，作为短路电流计算的基准 | 取决于接地方式，中性点直接接地系统中单独一相接地电流可达很大值 |
| 零序电流 | 严格为零 | 严格为零 | 很大，是接地保护的核心判据 |
| 故障相电压 | 两故障相电压大幅跌落，非故障相变化较小（不接地系统） | 三相电压不同程度跌落 | 故障相电压降为0，非故障相电压升高（不接地系统升高至线电压） |
| 对称性 | 不对称 | 三相对称 | 不对称 |
| 保护应用 | 相间短路保护、距离保护 | 最严重故障类型，用于整定 | 零序保护、接地距离保护 |
这个对比表格是我的参考框架。实际做仿真时，如果把三种故障都跑一遍放在同一个报告里，评审或者答辩时会直观很多，因为它直接体现仿真模型的可扩展性和对故障机理的理解深度。
4.3 序分量提取与理论相互验证
为了定量验证仿真模型的正确性，我在模型中增加了序分量提取环节。利用Simscape库中的Sequence Analyzer模块（傅里叶变换方式提取基波正序、负序、零序分量），在故障稳定运行区间（比如0.2到0.3秒之间）读取数值。
实测结果：
正序电流：约1.09kA 负序电流：约1.08kA 零序电流：接近0和理论计算对比：两相短路时正序电流等于负序电流，零序为零。1.09kA和1.08kA的微小差值是傅里叶变换窗长和系统衰减带来的误差，属于正常范围。
这个验证步骤的意义在于：你不仅看到了“波形像两相短路”，还从序分量维度确认了“模型确实是两相短路”。对于做保护算法研究的同行来说，这个环节几乎是强制性的，因为序分量直接对接保护算法的输入量。
5. 常见问题排查与实操经验记录
5.1 仿真报错的三种常见原因及处理
错误一：求解器不收敛或仿真速度极慢
现象：仿真运行到0.1秒附近时报错“Simulink cannot solve the algebraic loop”或者每一步计算时间极长。
原因分析：最常见的是模型中存在代数环——即某个模块的输出直接又反馈到该模块的输入，而且中间没有延时或者状态来打破这个环。在我的模型里，如果电压电流测量模块和故障模块的连接顺序不正确，或者测量信号被直接用于控制故障模块本身的触发时序，就会形成代数环。
处理办法：在故障模块的时序控制端加一个Memory模块或者Unit Delay模块，打破代数环；或者把故障时序改成内部设置不用外部信号控制。
错误二：波形出现严重的数值毛刺
现象：故障发生的瞬间，电流波形出现来回跳变的尖刺，甚至跳到几千千安级别，明显不符合物理规律。
原因分析：最大步长设置过大，仿真步长跨越了故障切换时刻，导致数值求解器无法精确捕捉开关动作瞬间的状态跳变。
处理办法：把最大步长减小到1e-4秒；同时如果使用ode15s或ode23tb算法，可以把相对误差从默认的1e-3调到1e-4或者1e-5，牺牲少量计算速度换取波形质量。
错误三：零序电流不是零，出现明显的非零数值
现象：理论上两相短路零序电流为零，但仿真结果显示零序电流有几百安培。
原因分析：电源侧的接地方式设置有问题。如果三相电源内部的中性点接地阻抗设置不为无穷大，那么系统的零序网络就有了通路，即使故障本身不接地，三相不平衡也会在零序回路中感应出电流。另外一个可能原因是在变压器（如果有）的连接组别设置中，中性点接地方式配置不对。
处理办法：检查Three-Phase Source模块中的Neutral parameter设置，高压电源中性点应该设置为Inf（无穷大）或Open，确保零序无通路。
5.2 从初学到精通的几条实操心得
心得一：先用固定参数跑通，再折腾变量。
最开始的模型一定要用最简单的参数把它跑通，不要一上来就搞什么故障电阻渐变、故障时间序列变化。我自己的经验是先固定一处故障位置、固定故障电阻0.001欧姆、固定故障时间0.1s-0.3s，确认波形正确后，再批量做不同故障条件的对比。如果一开始就加了一堆变量，万一结果异常，你都不知道是故障模型问题还是参数扫描设置的问题。
心得二：波形显示要合理设置坐标范围。
Scope窗口默认的坐标范围是自适应调整的，故障电流从几百安培跳到几千安培之后，故障前那段波形会被压扁，看不出正常工况下的幅值。我在观察波形时，经常将两个波形窗口并列——一个看全貌（0到0.5秒的总览），一个局部放大（0.08到0.35秒的故障区间），这样两方面信息都能捕捉到。
心得三：数据导出到工作空间方便后处理。
Scope波形虽然直观，但做定量分析还是要靠数据。我一般在模型里加To Workspace模块，把三相电流电压以结构体形式存入工作空间，然后用MATLAB脚本画出更美观的图、计算序分量、或者和理论值做误差分析。实测下来，这种流程比在Scope里肉眼观察靠谱得多，也方便生成报告图片。
6. 模型扩展思路与工程应用建议
6.1 从仿真模型到保护算法验证
两相短路仿真模型的真正价值，除了让你理解故障机理，更在于作为保护算法验证的测试平台。
我做过的几个典型的模型扩展方向：
方向一：距离保护动作逻辑验证。利用仿真得到的三相电压电流波形，计算测量阻抗，模拟距离保护在故障期间是否能正确测量到故障距离对应的阻抗值。这个可以配合Simulink里的继电器逻辑模块实现。
方向二：故障选相逻辑验证。基于两相短路时没有零序电流、负序电流不等于零的特征，写一个简单的选相判据，看看仿真数据能不能准确识别出这是AB相间故障而不是AC相或BC相故障。拿这套模型出来测试，选相逻辑的正确率一眼就能看出来。
方向三：并补电容器投入对故障特征的影响。在线路中点加入并联补偿电容器后重新仿真两相短路，观察故障电流的暂态特征变化，这对做高压输电线路保护研究很有参考意义。
6.2 给刚入门同行的一个参数速查建议
如果你刚开始做电力系统仿真，对参数选取没有概念，可以先用下面这组保险参数作为起点，跑通后再根据实际系统逐步调整：
系统电压： 110kV，50Hz 电源内阻抗： 0.002pu + j0.08pu 线路长度： 100km 线路单位阻抗： R = 0.01273Ω/km，L = 0.9337mH/km，C = 12.74nF/km 负载功率： 100MW + j10Mvar 故障电阻： 0.001Ω 故障时序： [0.1, 0.3] 仿真时长： 0.5s 求解器： ode23tb，最大步长 1e-4s，相对误差 1e-3这组参数跑出来，理论短路电流约1.875kA，仿真结果偏差通常在2%以内。等你把这套模型吃透了，再往不同电压等级（10kV配电网、220kV输电网）、不同故障类型（单相接地、两相接地、三相短路）、不同系统结构（双端电源、含变压器的多级网络）扩展，就有了一套可以反复使用的仿真基础框架，而不需要每次都从零开始搭。
我个人的体会是，Simulink仿真最费时间的反而不是搭模型本身，而是调参数和排查那些莫名其妙的报错。把一套标准模型的参数配置和常见问题经验沉淀下来，后面做任何故障分析都能快速进入状态。上面提到的这份参数速查和排查方法，是我踩过不少坑之后总结下来的，希望对同行有所帮助。
