简介：面向电力系统专业学生与工程技术人员的MATLAB/SIMULINK输电线路故障仿真分析文档，重点讲解短路故障的建模与仿真方法。内容从短路故障类型、MATLAB/SIMULINK基础入手，系统梳理短路计算原理与步骤，并给出三相短路系统仿真模型中的三相电压电流测量、变压器、三相线路、三相电源、线路故障及并联RLC负荷等模块的参数设置方法，附带算例与仿真结果分析。资料为docx格式，共1个文件，打包大小1.27MB，已有227人学习。读者可借助该文档快速掌握输电线路故障仿真流程，理解故障暂态过程，并用于故障定位、保护性能评估与系统稳定性分析等场景。
1. 项目整体思路与方案选型
1.1 为什么用MATLAB做输电线路故障仿真
先说个背景。我在电力系统方向摸爬滚打这些年，接触过的仿真工具不下五种，PSCAD、ATP-EMTP、PSASP都有涉猎，但遇到需要快速验证故障特性、分析保护算法、跑批量工况对比的场景，我第一反应还是开MATLAB。原因很直接：建模灵活、改参数方便、后处理能力强。
输电线路故障仿真这件事，本质上是在一个可控的数字环境里复现短路、接地、断线等异常工况下电压电流的暂态过程。现场不可能真拿一条110kV线路去反复做短路试验，成本和安全都不允许，所以仿真成了唯一可行的路径。MATLAB的优势在于，它既有Simulink这种图形化建模环境，又能用M脚本精确控制仿真流程、批量扫描故障条件、自定义故障类型和测量逻辑，这两者结合起来，做故障分析要比纯图形化软件灵活得多。
从项目实操角度说，这个标题里的“仿真分析”四个字很关键——它意味着交付物不只是跑出一条波形图，而是要能从波形里反推出故障类型、故障距离、暂态特征这些有工程价值的结论。我习惯的流程是：先在Simulink里建立输电线路模型和故障模块，搭好测量环节，再通过脚本批量修改故障发生时间、故障类型、过渡电阻等参数，最后对仿真数据进行FFT、对称分量法等分析，形成一套“建模—仿真—分析—结论”的完整链路。
1.2 方案选型：Simulink为主还是M文件为主
很多初学者纠结一个问题：Simulink和M脚本到底该用哪个？我给出的建议是：两者都用，但分工明确。
Simulink负责电路拓扑的搭建。输电线路、电源、故障模块、测量模块这些都是物理元件，在Simulink里拖拽连线，直观而且不容易出错。尤其是Three-Phase Fault模块，可以直接选故障类型（A相接地、AB相间短路、ABC三相短路等），还能设置故障时序，比手动搭故障网络省事得多。
M文件负责仿真控制和数据分析。比如你要研究不同过渡电阻下故障电流的变化规律，一个个手动改Simulink里的参数再点运行，效率太低。用set_param命令去修改故障模块的过渡电阻和故障时间，循环跑10次、20次仿真，最后一次性汇总波形和数据，这才是MATLAB真正的价值所在。
我实测下来，这样的组合方式在开发效率上比纯Simulink操作快三到四倍，而且后期出报告、画对比图也会顺手很多。如果只在Simulink里点点鼠标，仿真完了再截图导出，你会被大量的重复劳动淹没。
1.3 输电线路故障类型概述
做故障仿真之前，脑子里必须有一张清晰的故障分类表。电力系统里输电线路故障按相别和接地情况进行组合，最常见的是以下几类：
| 故障类型 | 故障相别 | 接地情况 | 工程中常见程度 |
|---|---|---|---|
| 单相接地短路 | A相（或B、C相） | 经电弧或过渡电阻接地 | 最高，约占80%以上 |
| 两相相间短路 | AB、BC、CA | 不接地 | 较少 |
| 两相接地短路 | AB、BC、CA | 经电弧或过渡电阻接地 | 时有发生 |
| 三相短路 | ABC | 一般不接地 | 很少但最严重 |
| 断线故障 | 单相或两相断线 | 可能伴随接地 | 偶发 |
其中单相接地故障是需要重点研究的对象。原因很简单：它在实际电网中占比最高，而且受过渡电阻影响大，故障特征复杂。我的仿真项目通常以单相接地作为默认工况，再根据研究需要扩展到其他故障类型，这样循序渐进，不容易乱。
2. 输电线路故障建模的核心原理与参数计算
2.1 线路模型的选择：集中参数还是分布参数
在Simulink里搭建输电线路，第一步是选线路模型。Simulink的Three-Phase Distributed Parameters Line和Three-Phase PI Section Line是两种常用模块，它俩的适用场景差别很大。
PI型线路模块用集中参数等效，原理是把整条线路的阻抗和对地导纳集中成几个π形电路。优点是结构简单、仿真速度快，适合长度在几十公里以内的短线路，或者对精度要求不高的教学演示。但线路长了之后，分布电容、波过程这些效应会变得明显，PI模型就不再准确。
分布式参数线路模块基于Bergeron模型，考虑了行波在输电线路上的传播时间，对长线路的暂态过程描述更精确。代价是仿真步长要取小，计算量明显增加。
我做这个项目时选用了分布式参数模型，线路长度取100公里。为什么？因为我要分析的故障暂态过程本身就包含高频分量，如果模型本身就不准，后续的波形谐波分析和行波故障定位研究就没有意义了。仿真步长我用的是固定步长，设成20微秒，这是在计算速度与波形分辨率之间反复权衡后确定的。20微秒对应50kHz的采样率，能覆盖到暂态中最关心的几次谐波。
2.2 故障模块的参数设计：过渡电阻与故障时序
Three-Phase Fault模块是Simulink电力系统库里的核心故障元件，它的参数设置直接影响仿真结果的物理合理性。这里面最容易踩坑的就是过渡电阻。
过渡电阻是指故障点电弧电阻加上杆塔接地电阻等综合等效电阻。金属性短路的过渡电阻接近0，但实际中的故障多是非金属性的。比如导线经过树枝放电，过渡电阻可能有几十欧姆；经过绝缘子表面污闪，可能到几百欧姆；经过高阻地面接地，甚至能达到1000欧姆以上。我在研究单相接地故障时，通常会做一组过渡电阻扫描：0.01Ω、10Ω、50Ω、100Ω、500Ω，观察故障电流幅值和相位的变化趋势。
故障时序设置也要讲究。Three-Phase Fault模块里有Fault Timing参数，可以设定故障起始时间和结束时间。这里有一个经验：故障发生时刻最好设置在电压峰值附近，这样得到的故障电流初始值最大，暂态特征最明显，便于分析。另外，故障持续时间的设置取决于你要研究的是暂态过程还是稳态过程。做故障行波分析，持续几个毫秒就够；做保护动作特性分析，可能要持续到保护动作逻辑完成，一般设置0.1秒以上。
我一般把仿真总时长设为0.5秒，故障在0.15秒时发生，持续0.2秒，这样既有足够的故障前稳态数据做基准，也留给故障后暂态和稳态足够的分析窗口。
2.3 Three-Phase Fault模块的进阶配置技巧
默认配置下，Three-Phase Fault模块是通过开关逻辑把故障相连接到地或者相间短接，但工程实际中的故障不是那么理想的。这里分享一个进阶配置：
故障相的选择。在模块的Fault Phases选项里勾选A、B、C相和Ground。做单相接地就只勾A相和Ground；做两相相间短路就勾AB相但不勾Ground；做三相短路就ABC三相连通但不接地。这些组合在模块内部会自动生成对应的短路网络，不需要你手动搭复杂的开关电路。
故障电阻的时变性。如果你想模拟电弧故障的动态过程，可以把故障电阻设成随时间变化的数组，通过信号输入的方式动态切换。虽然Simulink模块默认只支持固定电阻值，但可以通过外接一个受控电压源或者自定义故障子系统来模拟电弧的非线性特性。这个做法我在项目后期研究电弧重燃现象时用过，效果相当逼真。
3. 仿真系统的完整搭建与运行设置
3.1 Simulink整体拓扑结构
完整的输电线路故障仿真模型包括四个部分：电源模块、线路模块、故障模块、测量与记录模块。
电源模块我用的是Three-Phase Source，线电压有效值110kV，频率50Hz。这个电压等级对应的是高压输电线路，也是电力系统里最具代表性的仿真对象。电源的内阻我设得比较小，X/R比按典型发电厂升压变压器出口的特性来取，约为10。
线路模块就是前面说的分布式参数线路，接在电源和负载之间，故障点设置在距离电源端30公里的位置。这里注意，Three-Phase Distributed Parameters Line通常需要串联两个线路模块来精确控制故障点的位置：一段从电源到故障点，另一段从故障点到负荷末端。如果你只用一个线路模块，故障点只能固定在首端或末端，无法模拟线路中间故障。
测量模块使用Three-Phase V-I Measurement，同时测量电压和电流，输出给示波器和工作区变量。这个模块要放在故障点的电源侧，用来记录保护安装处感受到的电气量——这是故障分析最关注的位置。
系统末端接一个负载模块（Three-Phase Series RLC Load），有功功率100MW，无功功率30Mvar，模拟实际负荷。为了简化分析，负载保持三相对称。
3.2 仿真参数的合理设定
模型的仿真参数在Simulink的Configuration Parameters里设置。这里不建议用默认参数，我踩过这个坑。
求解器选择。我选用ode23tb，这是一款适合刚性问题的求解器，电力电子和电力系统中的非线性开关过程用它能有效避免数值振荡。固定步长模式下，步长20微秒；仿真时长0.5秒。固定步长比变步长更容易进行后期数据的等间隔采样，对FFT分析有好处。
容差设置。相对容差和绝对容差都设成1e-4，这个精度对于工频范围内的波形分析已经足够。如果再往小了设，计算时间会明显增加，但波形质量提升有限，性价比不高。
3.3 从Simulink高效导出数据
结束仿真之后，数据导出是下一步分析的关键环节。我习惯用sim命令从脚本里跑仿真，然后把Measurement模块的记录直接写到工作区。示例代码如下：
% 仿真控制与数据导出 mdl = 'transmission_line_fault'; load_system(mdl); % 设置故障类型: 1-A相接地, 2-AB相间短路, 3-ABC三相短路 set_param([mdl '/Fault'], 'FaultType', '1'); % 设置过渡电阻 set_param([mdl '/Fault'], 'Rfault', '10'); % 运行仿真 out = sim(mdl); % 提取三相电压电流波形 t = out.tout; Va = out.Va; Vb = out.Vb; Vc = out.Vc; Ia = out.Ia; Ib = out.Ib; Ic = out.Ic; % 保存数据到工作区，便于后续多种分析 save('fault_data.mat', 't', 'Va', 'Vb', 'Vc', 'Ia', 'Ib', 'Ic');这段代码的核心逻辑很直接：通过set_param修改模型参数，然后循环跑仿真。我在实际项目中，把这一层封装成了批量扫描函数，输入是故障类型数组和过渡电阻数组，输出是汇总好的所有波形数据。批量跑完20组工况，大约需要十几分钟，这比手动逐次点击仿真节省了大量时间。
有一点需要注意：Simulink模块的Simscape Electrical（即老的SimPowerSystems）接口在不同版本的MATLAB里参数名可能会有差异。我在R2020a和R2023b两个版本下都跑过，Rfault这个参数路径是一致的，但如果你用更早的版本，建议先在命令行里用get_param查一下参数列表，避免路径报错。
3.4 波形输出与初步判读
仿真跑完之后，首先观察的是三相电流电压波形。A相接地故障发生后，最显著的特征是故障相电流急剧增大、电压骤降，非故障相的电压电流则维持相对正常的水平。
用FFT工具分析故障电流分量，可以看到除了基波50Hz之外，还会出现明显的2次、3次、5次谐波。这些谐波成分与故障点的电弧特性、线路分布电容都有关系。对于刚入门的朋友，可以从基波含量的变化先判断故障的性质：金属性短路时故障电流基波分量最大，随着过渡电阻增大，基波分量呈下降趋势。
4. 常见问题与排查技巧实录
4.1 仿真不收敛或数值振荡
这是Simulink电力系统仿真中最常见的问题，具体表现是仿真中途报错、波形发散成一条直线或者剧烈毛刺抖动。多数情况是求解器选择和步长设置不合理。
排查思路是这样的：如果仿真一开始就发散了，先检查线路模块的分布式参数设置是否合理，电感电容值是否因单位换算错误导致数量级异常；如果是在故障切换瞬间发散的，优先尝试减小仿真步长，或者把求解器改回变步长的ode23tb、ode15s。固定步长模式下，步长不要超过系统最小时间常数的十分之一。
关于步长选择，我自己的经验是先用变步长跑通模型，观察仿真输出的自然采样周期，再据此设定固定步长，这样能兼顾准确性和效率。上来就拍脑袋设一个固定步长，很容易因为步长过大丢失暂态峰值。
4.2 波形毛刺过多、高频噪声明显
故障发生瞬间，电流电压波形上会出现高频振荡，这是正常的暂态现象。但如果振荡持续不衰减，或者毛刺幅度大到淹没了主要波形，就有问题了。
一个主要原因是测量模块的采样率不足。如果你想分析暂态高频分量，仿真步长就要相应缩小。我之前做过一次对比，步长从50微秒改到10微秒之后，故障行波附近的波形细节丰富了很多，这对后续基于行波原理的故障测距研究至关重要。
另一个容易被忽视的问题是模型里缺少必要的阻尼元件。在故障点或线路末端并联合适的电阻可以加速暂态衰减，让稳态波形更快呈现出来。不过加阻尼电阻要谨慎，不要影响正常工况下的潮流分布。
4.3 故障模块的触发时刻不准确
我在初学阶段遇到过一种情况：设置了故障0.15秒发生，但波形显示故障时刻明显不对。后来发现，问题出在故障模块的Switch Time设置上。Three-Phase Fault模块属于Breaker类元件，其动作时间必须严格按照Transition Time的设置值，且需要确保外部控制信号的逻辑与之匹配。如果你在外部信号和内部定时的两套控制逻辑中混用，就会产生冲突。
所以我的习惯是：统一使用模块内部的Fault Timing，不启用外部控制端。这样只需要管理一套时间参数，排查问题也更简单。
4.4 常见问题速查表
| 现象 | 可能原因 | 解决方案 |
|---|---|---|
| 仿真直接报错 | 参数单位错误、连接错误 | 检查线路和故障模块参数，确认三相连接 |
| 波形发散无规律 | 步长过大、求解器不匹配 | 换用ode23tb，减小固定步长 |
| 故障相电流异常偏小 | 过渡电阻设置过大 | 检查Rfault值，按工况合理设定 |
| 电压跌落不明显 | 系统容量过大或负载过小 | 调整电源短路容量和负载功率 |
| 频谱分析出现混叠 | 采样率不足 | 减小步长，提高等效采样频率 |
| 批量仿真时参数没更新 | set_param路径错误 | 使用get_param确认路径后修改 |
4.5 独家避坑技巧
根据我个人经验，有几个坑是反复遇到的，这里集中分享出来。
第一，保存模型之前一定要先跑通一次完整仿真。我经历过太多次改了几行代码、保存模型再打开后，Simulink模块之间的连接出现异常的问题。运行前用Simulink的诊断检查功能扫描一遍，能省掉大量排查时间。
第二，批量仿真一定要在脚本里设置warning off，把模型编译过程中的冗余警告压制掉，否则循环20次仿真，光刷屏警告就够你烦躁的。但注意，真的遇到了报错信息，一定要把警告重新打开看清楚，不能一刀切。
第三，FFT分析时，Data from workspace的采样时间要与你仿真步长严格对应，否则频谱会多出很多不存在的频率分量。我习惯在仿真结束之后，对数据先做一次重采样，统一到1kHz或者10kHz的采样率，再做频谱分析。
5. 从仿真数据到工程结论的进阶分析
5.1 对称分量法提取故障特征
拿到三相电压电流数据之后，最核心的分析方法是对称分量法。把三相不对称的故障量分解成正序、负序、零序三个对称系统，可以直接得到故障类型判别的判据。
MATLAB里做对称分量法有三种常用途径：自建变换矩阵、调用phase2seq函数、或者用Simscape Electrical里的Sequence Analyzer模块。前两种适用于脚本分析，第三种适合在线监测。我的项目里用的是自建变换法，代码如下：
function [seq] = sym_components(abc) % 对称分量变换矩阵 a = exp(1j * 2 * pi / 3); A = [1 1 1; 1 a^2 a; 1 a a^2]; seq = A \ abc; end用这个函数把A相故障时的三相电流做一次变换，你会发现正序电流幅值最大，负序电流有一定数值，零序电流则取决于系统接地方式。如果是中性点直接接地系统，单相接地时零序电流会非常明显，这也是零序保护能够正确动作的物理基础。
5.2 故障测距的基本逻辑
故障测距是输电线路故障仿真分析里最有实际价值的延伸方向。简单说，它的原理是利用故障时电压电流的工频分量，通过线路的阻抗参数计算出故障点到测量端的距离。
在Simulink模型里，故障点设在30公里处，利用仿真得到的电压电流量反算故障距离，就是验证测距算法的最好方式。具体做法是：提取故障后稳态的电压和电流相量，按线路的正序阻抗参数估算距离。
这里要注意的是，过渡电阻会让测距结果出现偏差。金属性短路时测距误差较小，但过渡电阻增大到100欧姆以上时，单端测距结果会有明显偏移，这是实际工程中的行业性难题。一个缓解思路是通过两侧数据同步测量或者引入故障分量法进行修正。这些内容可以带着仿真数据作为支撑单独展开研究。
5.3 批量工况对比的价值
仿真分析做得多了之后，我最大的体会是：单次仿真的意义极其有限，批量工况对比才能体现建模工作的价值。
把故障类型、过渡电阻、故障位置、负荷水平这些变量排列组合之后批量仿真，用surf或contour图展示故障电流随不同参数的变化关系，会得到很直观的工程结论。比如你可能会发现：过渡电阻从10欧姆升到500欧姆时，故障电流下降不是线性的，而是呈现出明显的非线性衰减特征，并且故障位置越靠近末端，这种衰减越明显。
这些规律性的结论，无论在学术论文写作里还是在实际保护整定计算中，都是非常有说服力的依据。建议各位在做仿真的时候，不要只满足于“跑一次仿真相当于做一次实验”的层次，而要把仿真平台当成一个能够遍历各种工况的“数字试验场”，这才能发挥出MATLAB的全部潜力。
本文还有配套的精品资源，点击获取
