单相接地故障是配电网里最头疼的故障类型之一，尤其是中性点不接地或者经消弧线圈接地的系统里，故障电流小、特征不明显，现场选线一直是老大难问题。用MATLAB 2022a把它建模仿出来，等于把故障过程搬到了电脑里，波形成什么样子、电压电流怎么变化、零序分量怎么提取，都能看得清清楚楚。这篇博文就围绕这个配电网单相接地故障仿真模型，把设计思路、核心建模细节、实操步骤和踩坑经验一次性讲透，适合电力专业学生、继电保护工程师和对配电网仿真感兴趣的入门者。
1. 模型整体设计思路：先搞懂要仿真的对象
1.1 仿真模型解决什么问题
很多刚接触配电网仿真的人一上来就找模型、跑仿真，却忽略了一个最基本的问题：这个模型到底要回答什么。这套模型表面上是在做一个10kV配电网的单相接地故障仿真，再附带两相接地对照，但实际上它解决的是两个层面的需求。
第一层是认知需求。单相接地发生后，非故障相的电压会升高到原来的根号3倍，零序电压出现，故障相电流有变化但变化幅度不大，这些结论在书本上只是几个公式，但仿真可以让你直观看到波形变化的过程。两相接地对比就更直观了，一个故障电流很小、一个故障电流很大，放在一起看才能对配电网故障特性有真正的体感。
第二层是工程需求。配电网单相接地故障的选线、定位、保护整定，都依赖对故障暂态过程的准确理解。模型跑出来的零序电流波形和零序电压波形，就是后续算法研究的基础数据。所以这个模型不是交作业用的演示模型，而是能够支撑故障特征分析、保护策略验证的一个基础平台。
1.2 为什么选择MATLAB Simulink来搭
电力系统仿真工具有很多，像PSCAD、EMTP、BPA等，都是老牌电力系统仿真软件。但这个场景我更推荐MATLAB Simulink里的Simscape Electrical库（老版本叫SimPowerSystems），原因有三个。
第一个原因是元件库足够全面。三相电源、变压器、分布参数线路、三相故障模块、电压电流测量模块都是图形化的，直接拖进画布连线就行，不需要像EMTP那样手写卡片，对入门者非常友好。
第二个原因是波形查看和数据后处理方便。Simulink里的Scope示波器可以直接看波形，To Workspace模块可以把数据导到MATLAB工作区继续做FFT分析、特征提取。做配电网故障研究的人，终究要写代码处理数据，MATLAB生态天然衔接这一步。
第三个原因是版本适配宽松。MATLAB 2022a的实际运行环境对主流的64位Windows系统支持很稳定，Simscape Electrical在这个版本上无论是元件库完整性还是求解器稳定性都处于比较好用的状态。只要电脑配置不是太老，跑这个级别的配电网模型完全不吃力。
1.3 单相接地和两相接地为什么要一起做
做故障仿真的很多人容易忽略模型边界条件的完整性。只做单相接地，你看到的只是小电流接地系统中最典型的故障形态；只做两相接地，又缺少和单相接地对比的尺度。这套模型把两种故障都放进来了，其实是用同一个平台、同样的一套参数，通过切换故障模块的设置来完成对照实验。
这样做的好处非常明显。第一，模型参数完全一致，排除了因为系统参数不一致导致的波形差异，对比结论更可信。第二，一次搭好的模型，改故障类型只需要改故障模块里的选项，工作量小。第三，从工程角度看，两相接地故障往往意味着保护需要更果断地切除，放在一个仿真模型里做对比，能帮助你建立对“故障严重程度”的量化感知。
2. 核心建模细节与参数设置
2.1 系统结构：从电源到负荷的完整链条
我这套模型的整体拓扑是标准的10kV配电网结构：一个无穷大系统等效电源、一台升压变压器（系统侧110kV到10kV）、一条10kV母线，母线上接了两条馈线，一条长12km，一条长8km，末端分别带恒阻抗负荷。故障点设置在长馈线距离母线5km的位置，这样故障点前后都有线路段，能模拟出完整的电压跌落和零序分布。
这里有一个细节很多人会忽略：系统的中性点是悬空的。10kV配电网大多是中性点不接地或经消弧线圈接地，这是小电流接地系统的根本特征，也是单相接地故障电流很小的原因所在。如果误把变压器中性点直接接地，仿真出来的单相接地故障就是大电流接地系统的特性，波形特征完全不对。所以在2022a版本的Simulink里搭变压器时，要检查绕组连接组别，确保Yg绕组在低压侧而不是高压侧，或者干脆高压侧用Y（不引出中性点）。
2.2 电源、变压器和线路的参数选取
电源模块参数上，我选择Three-Phase Source，线电压设置为110kV，频率50Hz，短路容量设置为1000MVA，这个数值对应的是系统侧的一次变电所容量，足够模拟一个相对刚性的电源。变压器额定容量取10MVA，变比110/10.5kV，连接组别用的是常见Yd11，这样10kV侧就没有直接接地的中性点，符合中性点不接地系统的前提。
线路参数是单相接地仿真的关键。我用分布参数线路模型（Distributed Parameters Line），三分之二的架空线共性参数按照典型10kV架空线的单位参数设置：正序电阻0.27/km，正序电感1.09mH/km，正序电容0.0095F/km，零序电阻2.7/km，零序电感4.2mH/km，零序电容0.006F/km。如果你手头没有具体的线路参数，用这些典型值作为起步是完全可以的，后续再根据实际工程数据修正。
负荷我用的是Three-Phase Series RLC Load，每条馈线末端负荷设为2MW、0.4Mvar。恒阻抗负荷有个好处，就是仿真收敛性比恒功率负荷好，对故障暂态过程的影响也更直观。刚开始搭模型不建议上感应电动机这类动态负荷，会引入额外的机电暂态过程，干扰你对故障特征的判断。
2.3 故障模块设置：把单相接地和两相接地区分开
Simulink的Three-Phase Fault模块是这个模型的核心开关。模块支持按时间和外部信号两种方式控制故障切换，我建议用默认的Timed（定时控制）方式，因为操作简单直观。
单相接地故障的设置方法是：在Fault Phase的选项卡里勾选A相，然后勾选Ground（接地），把故障电阻设成小阻值。这个故障电阻很讲究，工程上最典型的是金属性接地（0.01~0.5），还有一种带过渡电阻的情况（比如经树枝接地，阻值10~100）。仿真里想先看最清晰的波形，建议从金属性接地开始，阻值设0.01，后续再改成10、50、100做对比。Transition Time设置成故障发生和切除的时刻，我习惯设成[0.1 0.3]，意思是0.1秒发生故障、0.3秒切除，仿真到0.5秒能完整看到故障前、故障中、切除后三个阶段。
两相接地的设置更简单：在Fault Phase勾选A相和B相，同时勾选Ground，其他参数不用改动。同样是0.1秒发生、0.3秒切除，故障电流会显著上升，波形特征和单相接地放在一起有非常强烈的视觉对比。
故障模块还有一个容易被忽略的专项参数：Transitional Impedance（过渡阻抗）。故障导通时过渡阻抗应该设得很小，比如0.01，模拟可靠导通；故障断开时过渡阻抗设置为很大，比如1e6，模拟完全开路。这个参数如果设反了，会出现故障模块关不断或者漏电流的情况，波形怎么调都不对。
2.4 零序电压和零序电流的提取方法
单相接地故障分析里有三个核心观测量：三相电压、三相电流、零序电压和零序电流。前两个直接用Three-Phase V-I Measurement模块就能测量，但零序分量需要额外处理。
我用的方法是先把三相电压电流通过测量模块变成Simulink信号，然后用Demux把A、B、C三相分开，再用三个加法器（Add模块）相加。三相电压之和除以3就得到零序电压，三相电流之和除以3得到零序电流。也可以用Simscape Electrical里面自带的Zero-Order Sequence Transformer模块，放在信号线上直接出零序分量，效果一样，但新手容易在信号类型转换上出问题，所以我建议老老实实用加法器+增益模块，原理清晰，排错方便。
零序电压和零序电流的比例关系也有讲究。小电流接地系统里，零序电压基本是系统对地电容电压的重分布结果，幅值约为相电压水平；而零序电流是线路对地电容电流，幅值往往只有几安培到几十安培。这个数量级的差异会在后续判断故障线路时非常重要，看波形时可以把零序电流通道适当放大，不要因为显示尺度问题误判为零序电流不存在。
3. 实操过程：从搭模型到拿到正确波形
3.1 仿真参数与求解器设置
模型搭好之后，仿真参数的设置直接影响能不能跑通、波形准不准。Simulink默认的求解器配置跑简单电路模型问题不大，但配电网含分布参数线路和电力电子开关，仿真步长的选择需要有一点讲究。
我实测下来，采用变步长求解器（Variable-step）、求解器选用ode23tb是比较稳的组合。ode23tb对刚性系统有优化，而配电网线路模型里的电感和电容时间常数差异很大，刚性问题的可能性很高，用ode23tb能减少很多不收敛的报错。最大步长设置为1e-4秒（也就是0.1ms），仿真停止时间0.5秒。这样既保证了暂态过程的精度，又不至于让仿真速度慢到无法接受。
还有几个仿真选项建议按下面这个表设置，排错的时候能省很多时间：
| 参数项 | 推荐设置 | 说明 |
|---|---|---|
| 仿真时长 | 0.5 s | 覆盖故障前0.1s、故障中0.2s、切除后0.2s |
| 求解器 | ode23tb | 对含线路模型和刚性系统更友好 |
| 最大步长 | 1e-4 s | 保证暂态波形细节不丢失 |
| 零穿越检测 | 开启 | 避免故障切换时刻的检测问题 |
| 信号采样时间 | 连续 | 后续导出数据到工作区更方便 |
3.2 单相接地故障的波形特征与判据
仿真跑通后，第一个应该打开的示波器是三相电压。正常情况下三个相电压是对称的，A相接地后，A相电压幅值迅速压向零（考虑故障电阻很小的情况，A相电压会跌落到一个接近零的数值），而B相和C相的电压幅值会明显升高，达到原来的1.7倍左右，也就是接近线电压的水平。这是单相接地最直观的波形特征。
第二个必看的是零序电压。正常运行时零序电压接近零，故障发生后零序电压在短时间建立起来，幅值大约等于相电压。如果看到零序电压升高的同时，零序电流跟着出现，而且故障线路的零序电流比非故障线路明显大，这个故障特征就闭环了。
我还建议打开三相电流的波形看看。金属性单相接地时，故障相的电流变化并不像一些人想象的那么大，这是因为中性点不接地系统的单相接地电流本质上是全系统对地电容电流，数值有限。很多人第一次看到这个波形时很困惑：为什么故障相电流看起来没什么特别大的变化？这正是小电流接地系统的特点，也是为什么单相接地保护比短路保护困难的原因。
3.3 两相接地故障的波形对比
把故障模块改成A、B两相接地，重新跑一遍仿真，对比就出来了。B相电压在故障期间明显降低，A相也是接近零，C相电压升高，但更重要的是线电流会大幅上升，达到几千安的级别，这跟单相接地只有几十安电容电流有本质区别。在人工排查或者保护配合时，两相接地故障属于必须快速切除的故障类型，原因就在这里。
波形对比的时候我有一个习惯：把单相接地和两相接地的波形截图放在同一个时间轴下对照，标注出电压幅值变化比例、电流幅值数量级、零序分量是否稳定这几个关键维度。这样无论自己后续写分析报告还是给别人讲原理，都会非常直观。
3.4 数据导出与后处理
Scope里面看波形只是第一步，真正做研究或者写报告时，需要把仿真数据导出来分析。在Three-Phase V-I Measurement模块的输出端口接一个To Workspace模块，变量名设成V_abc和I_abc，保存格式选Array。仿真结束后这些数据会自动存到MATLAB工作区。
用MATLAB脚本可以快速算零序分量并绘制曲线。我个人习惯是先跑一次0.5秒的仿真，然后把数据导出来，用脚本把零序电压和零序电流画在同一张图里，同时标注出故障开始和切除的时间刻度。这样生成的图既清晰又能直接用于论文或者技术报告。
4. 常见问题与排查技巧实录
4.1 仿真报错不收敛，先查这几处
配电网模型仿真失败的案例我见过不少，最典型的报错就是Simulink提示无法收敛或者求导失败。这时候不要慌，排查顺序我建议：先看系统是否存在代数环，再看采样时间和步长设置，最后看故障模块的过渡阻抗有没有设成0值，尽量避免把过渡阻抗设成精确的0，用0.001之类的极小值代替会更稳妥。
还有一个常被忽略的点：如果变压器中性点没有处理干净，仿真时会直接报错或者算出完全不合理的波形。检查一下变压器模块的绕组连接方式，确认10kV侧是Y型连接但没有接地引线，这样才能保持中性点不接地系统的性质。
4.2 波形毛刺多、有高频振荡怎么处理
波形出现大量毛刺和振荡，多半不是模型本身的错，而是数值计算在高速开关瞬间产生的。故障模块在0.1秒动作切换时，线路电感电流需要突变，数值上就会产生尖峰。解决方法是适当调小最大步长，同时开启零穿越检测。还有一种稳妥的处理方式，就是在仿真结果输出后用MATLAB的滤波函数做后处理，但不建议一上来就滤波，因为故障暂态本身包含有用高频信息，盲目的平滑可能丢掉关键特征。
如果高频振荡非常严重且不衰减，那要检查是不是线路模型选成了集中参数且分段数太少。真实配电网的线路有分布式特性，可以用分布参数模型，或者把集中参数改成多段pi等效，再跑一遍对比波形。
4.3 零序电压测出来一直是0
这个问题很常见。原因基本只有一个：系统对称性太好，或者测量方法不对。如果系统完全对称并且负荷也完全平衡，单相接地时零序电压仍然会出现，因为故障破坏了系统对称性，所以正常来说零序电压一定会出现。
测出来一直是0，要先检查加法器部分接线。用Demux分离三相之后，要注意别把测量信号里的下标顺序搞错，A、B、C分别对应哪个信号，错了的话相加结果自然不对。还有一个检查点：电压测量模块的输出单位是V还是kV，如果单位不统一，算出来的零序分量会差一个数量级，看起来像是接近0。
4.4 故障模块切换时刻的坑
最后说一个我在实操中踩过的坑。Three-Phase Fault的Transition Time如果设置成[0.1]，这个值会被当成两个状态切换时间，也就是0.1秒发生故障，之后永远保持在故障状态，直到仿真结束。我一开始参照某些教程只写了一个时间点，导致0.3秒后系统并没有恢复，分析切除过程时完全对不上。
所以如果你需要看故障切除后的暂态过程，一定要写[0.1 0.3]两个时间点；如果只想看故障持续状态，用[0.1]也没问题，但心里要清楚这个含义。
5. 模型可以做哪些后续扩展
5.1 从故障模型走向选线算法验证
模型跑通了波形，下一步自然就是干活。配电网单相接地故障研究里最大的工程问题是故障选线，也就是在母线有多条出线的情况下，怎么判断是哪条线路发生了接地。现有的选线方法，比如零序电流幅值法、零序电流相位法、五次谐波法、暂态首半波法等等，本质上都需要先有一个准确的故障模型来产生测试信号。
这套模型里在母线上加了几条馈线，每条馈线上装一套零序电流测量，就能拿出来当不同选线算法的测试平台。把算法用MATLAB脚本实现，跑完仿真后读取各路零序电流数据，在指令窗口里看看算法的判断结果和实际设定的故障位置是否一致。这样比对着理论公式空想要有价值得多。
5.2 加入消弧线圈和不同中性点接地方式
我国配电网里有大量经消弧线圈接地系统，故障特征和中性点不接地系统不尽相同。要扩展这个方向并不难：在变压器中性点位置加一个可控的消弧线圈模块，设置成过补偿状态（补偿度一般在8%~10%左右），再跑同样的接地故障，观察零序电流幅值和相位的变化。
消弧线圈的引入会让故障线路的零序电流方向和幅值特点发生变化，很多选线方法在这种场景下会失效，这个话题很值得研究。Simulink模型里加入消弧线圈就是多一个RLC串联支路的事，核心参数是电感值的整定，可以根据系统对地电容电流计算出来。
5.3 故障定位和更复杂的场景
再往后，还可以往这个模型里加故障定位算法的验证场景。配电网单相接地的故障定位比输电网难得多，因为故障电流不明显，传统的阻抗法很难适用。可以考虑用故障暂态行波法，模型里的分布参数线路能产生比较真实的暂态行波过程，在Simulink里观察行波在故障点和母线之间的反射，从而估算故障距离。这个方向需要更高的采样率，对仿真步长设置要求也更苛刻，但模型的核心结构完全不用动。
还有一种常见的扩展思路是把故障模型嵌入到更完整的配电自动化仿真环境里，比如加入馈线自动化逻辑（FA）。把Simulink模型跑出来的故障特征信号导出后，接入用Stateflow编写的馈线自动化控制逻辑，形成带保护逻辑的闭环验证环境。这样整个模型就不只是一个波形演示器，而是一个可以反复验证控制策略的仿真平台。
我个人在实际操作中的体会是：能调通一个经典的配电网单相接地故障模型，比单纯看十篇论文都有用。波形上每一个特征的变化，背后都是系统结构、中性点方式、线路参数、故障条件的综合体现。把模型当成工具，以它的参数调整作为实验变量，才能真正理解配电网接地故障的本质。
