企业官网建设流程全解析

继电保护这块，小电流接地系统的单相接地故障仿真，几乎是每个做配电系统仿真的人都会碰到的课题。我自己最初做这个方向的时候，最难的不是保护逻辑本身，而是“模型到底该怎么搭”——中性点不接地和经消弧线圈接地，这两种方式在Simulink里的接法只差一个支路，但波形特征和保护判据完全不同。这篇内容我把整套仿真模型的搭建思路、消弧线圈参数计算、故障设置细节以及常见坑都整理出来，正在做课程设计、毕业设计或者工程预研的朋友可以直接照着上手，省掉自己反复试错的时间。

1. 仿真之前，先搞懂小电流接地系统在仿什么

1.1 两种接地方式：原理差异直接决定模型结构

小电流接地系统，通俗讲就是中性点不直接接地的系统，常见的有中性点不接地和经消弧线圈接地两种。所谓“小电流”，指的是单相接地故障时，流过故障点的电流不是靠电源短路电流维持的，而是靠线路对地电容形成的容性电流，数值通常只有几安到几十安。

中性点不接地系统（工程上也叫IT系统），电网和大地之间只有三相对地电容这条“弱联系”。发生单相金属性接地时，A相电压被拉到接近0，B、C两相电压升高到线电压，也就是从原来的5.77kV涨到10kV左右（以10kV系统为例），三相线电压仍然对称，所以系统还能带故障运行一到两个小时。但非故障相电压抬高，对绝缘是个考验。

经消弧线圈接地系统，就是在中性点与大地之间串入一个电感线圈。单相接地时，消弧线圈提供一个感性电流，去抵消三相对地电容产生的容性电流，故障点残流变小，电弧更容易自行熄灭。这样做最大的好处是减少弧光接地过电压，代价就是故障特征变得不明显，继电保护选线变得很困难。

从仿真建模的角度看，两种系统的区别本质上就在一根线上：中性点对地之间接不接消弧线圈。这给了我们一个非常好的便利——搭好一套主电路，通过切换中性点支路就能跑出两套完整的结果。

1.2 消弧线圈补偿的参数计算

消弧线圈的补偿效果，核心是看电感L选多大。推导并不复杂：单相接地后中性点电压偏移到相电压Uφ，零序等效网络里三相对地电容3C0和消弧线圈L是并联的。完全补偿时，要求电感电流等于电容电流：

Uφ / (ωL) = 3ωC0Uφ

两边消掉Uφ，得到：

L = 1 / (3ω²C0)

我以一个10kV、50Hz、出线10km的简单系统为例。架空线对地电容取经验值C0=1μF每相（大约0.1μF/km），那么：

U_ll = 10e3; % 线电压 10kV U_phase = U_ll / sqrt(3); f = 50; omega = 2 * pi * f; C0 = 1e-6; % 每相对地电容 I_cap = 3 * omega * C0 * U_phase; % 单相接地时的总电容电流 L_comp = U_phase / (omega * I_cap); % 完全补偿电感 L_over = U_phase / (omega * 1.1 * I_cap); % 过补偿10% fprintf('总电容电流: %.2f A\n', I_cap); fprintf('完全补偿电感: %.3f H\n', L_comp); fprintf('过补偿10%%电感: %.3f H\n', L_over);

算出来大约是：总电容电流5.44A，完全补偿电感3.38H。工程上不会用完全补偿，因为三相电容并不完全对称，完全补偿容易引起串联谐振过电压，通常采用过补偿运行，让残流呈感性。取过补偿度10%时，IL=1.1×5.44≈6.0A，对应的电感约为3.07H。

注意：消弧线圈实际不是纯电感，它带有阻尼电阻（有功损耗）。仿真时最好在电感旁边并联一个几百欧到几千欧的电阻，否则LC谐振回路的阻尼接近零，波形会出现明显的高频振荡，数值仿真也容易发散。

2. Simulink模型完整搭建：模块清单与接线要点

2.1 主电路拓扑与模块选型清单

搭建平台是MATLAB/Simulink里的Simscape Electrical Specialized Power Systems库，早一点的版本叫SimPowerSystems。核心思路是用三个单相AC Voltage Source搭出一个浮动星形电源，把中性点引出来，这样切换到消弧线圈接地时非常方便。

为什么不用现成的Three-Phase Source模块？因为Three-Phase Source的内部接法通常是Yg，中性点在模块内部直接接地，引不出来。虽然可以靠一台Yd接线的主变压器把高压侧中性点引出来接消弧线圈，更贴近实际工程，但模型会多一层变压器，对初学者不友好。用三个单相交流源构造虚拟中性点，是最快、最容易理解的做法，切换对比也很直观。

模块清单如下：

三个单相AC Voltage Source的幅值必须填相电压峰值，不是线电压有效值。10kV系统的相电压有效值是10k/√3=5774V，对应峰值是8165V。如果直接填10000V，整条母线电压会被抬高√3倍，后面的零序电压和电容电流全都不对。这是仿真里很容易犯的第一个错。

线路阻抗方面，10km架空线路的正序电阻大概在1Ω左右，正序电感大概0.01H，对单相接地故障影响不大，但建议保留，避免模型过于理想化导致电流波形失真。对地电容才是关键，每相1μF对应了大约5.4A的接地电容电流，已经能明显看到故障特征。

2.2 零序电压和零序电流的测量方法

小电流接地系统的故障分析，核心看两个量：零序电压3U0和零序电流3I0。在Simulink里提取这两个量有几种方式，最直观、最不容易出问题的就是直接做三相瞬时值的代数和。对称分量法里，三相电压相加就等于3倍零序分量，瞬时值同样满足这个关系。

具体做法：用Three-Phase V-I Measurement测出三相电压和电流后，拉一个Simulink的Sum模块，把Va、Vb、Vc加在一起就是3U0，把Ia、Ib、Ic加在一起就是3I0。后续可以再用RMS模块看有效值，或者接Fourier模块提取基波幅值、相位。

如果想要直接看到有效值，我习惯在Sum后面再接一个RMS模块，设置基波频率50Hz。这样观察稳态幅值非常方便，尤其是比较不同线路的零序电流大小，比在Scope里肉眼对比波形可靠得多。

2.3 故障环节与保护动作逻辑的设置

Three-Phase Fault模块的设置有几个关键点。在“Fault between Phase A and Ground”前打勾，故障电阻Ron填0.001Ω模拟金属性接地；如果想研究经过渡电阻接地，改成10Ω、100Ω甚至1000Ω就可以，但波形幅度会明显变小。

故障时刻设置要注意状态向量和切换时间是一一对应的：Transition statuses填[0 1 0]，Transition times填[0.1 0.3]。意思就是0到0.1s正常运行，0.1s时故障投入，0.3s时故障切除。很多新手在这里把状态向量写反，导致模型从头到尾都在故障状态，或者永远不故障，需要花点时间排查。

如果要演示保护动作，可以在故障线路首端串一个Three-Phase Breaker，控制信号做成“零序电压启动”逻辑。比如用3U0的RMS值和定值比较，输出经过一个滞环Relay，再接一个RS触发器锁存跳闸信号。为什么必须用RS触发器锁存？因为故障特征在断路器跳开后会立即消失，如果不锁存，跳闸信号跟着复位，断路器可能马上又合上，形成反复分合。

3. 仿真结果对照解读：零序量在两种系统里差别有多大

3.1 中性点不接地系统：故障电压抬高和零序电流分布

仿真跑起来以后，最直观的现象是：0.1s之前三相电压对称，每相相电压有效值5774V；0.1s A相接地后，A相电压瞬间塌到接近0，B、C相电压升高到接近10kV，零序电压3U0的有效值大约17.3kV，也就是3倍相电压。

这个结果和理论完全吻合。需要注意的是，B、C相电压升高后会有短暂的高频振荡，这是对地电容与线路电感的自然频率响应，不是模型错误。后面会平稳到10kV左右。

零序电流的分布更有意思。非故障线路首端测到的3I0，就是本线路三相对地电容电流的矢量和，数值大约等于本线路的电容电流；故障线路首端测到的3I0，是所有非故障线路对地电容电流的叠加，因为故障点必须由全系统电容电流汇集提供，而中性点又不接地，这个电流路径只能通过故障线路构成回路。所以故障线路的零序电流明显大于任何一条非故障线路，方向是线路流向母线（以母线指向线路为正方向看是反的）。这正是传统零序电流比幅比相选线能够成立的基础。

3.2 经消弧线圈接地系统：残流被补偿后发生了啥

把中性点对地的消弧线圈接入模型，电感按前面算的3.07H设置，并保留阻尼电阻，重新跑一遍。故障相电压、非故障相电压、3U0的波形和不接地系统几乎一样，非故障相依然升高到线电压，零序电压依然在17.3kV左右。原因是消弧线圈只补偿零序回路中的无功电流，对电压分布基本没有影响。

差别集中在零序电流上。不接地系统里故障线路和非故障线路的零序电流幅值差异明显，经消弧线圈接地后，故障线路上的容性零序电流被消弧线圈提供的感性电流大幅度抵消，故障点残流从5.4A左右降到不足1A，这就是消弧线圈“消弧”效果的来源。从两条出线首端看，故障线路与非故障线路的零序电流幅值变得接近，相位关系也变得模糊，此时再靠“哪条线路零序电流大就选哪条”的逻辑，大概率会误选。

3.3 从波形反推保护逻辑该怎么设计

两套模型对比下来，能很清楚看到一个问题：零序电压启动判据在两种系统下都成立，故障前后3U0从接近0跳到17.3kV，是一个非常可靠的启动量。但零序电流幅值保护在经消弧线圈接地系统里就不太可靠了。

所以实际工程里，中性点不接地系统可以尝试零序电流比幅比相，经消弧线圈接地系统更多依靠零序电压启动，再配合暂态量选线、有功分量法或者注入信号法。做仿真验证时，可以在模型中同时输出故障线路和非故障线路的零序电流相位，自己画一下矢量关系，理解会深很多。

两种系统主要电气量对比如下：

4. 常见问题与排查技巧实录

4.1 断路器、故障模块经常被忽略的设置

用Three-Phase Breaker做保护跳闸时，第一坑就是控制方式。模块内部有Switching times这个参数，一旦填了时间，外部控制信号就失效了。如果发现给了跳闸信号但开关纹丝不动，先去看这个参数是不是留空了。其次看Initial status，初始状态0表示断开，1表示闭合。很多模型跑出来一开始就没电，就是这里设成了0。

另一个典型问题是断路器断开瞬间仿真发散。断开感性电流会产生很高的电压尖峰，如果没有吸收回路，数值求解器容易“爆掉”。解决办法是给断路器合理设置Snubber resistance，或者在三相断路器旁边并联RC吸收支路。对纯仿真来说，如果只是想模拟故障切除，直接用Three-Phase Fault模块的第二次切换时间更省事，没必要强制走断路器。

4.2 Powergui与仿真步长选不好就发散

Specialized Power Systems库里所有模型都必须至少放一个Powergui模块，这是仿真的前置条件。很多新手的模型报错就是因为忘了拖这个模块进来，或者放了一个但没有设置仿真模式。

我建议直接用Discrete离散模式，Ts设50e-6s，也就是50微秒。对50Hz工频系统来说，每个周期采样400个点，波形足够平滑。连续模式下，消弧线圈和对地电容构成的谐振回路是典型的刚性系统，ode23tb有时候会特别慢，甚至卡死。如果必须要用连续模式，优先试ode15s，并且把相对容差放到1e-4左右。

4.3 消弧线圈参数和数值稳定性的坑

消弧线圈不是光给一个电感就行。三相对地电容加消弧线圈电感，本质上是一个几乎无阻尼的LC谐振回路，如果仿真里只有L和C，功率振荡衰减极慢，Scope里看到的波形会叠加明显的高频分量，影响结果判断。常见做法是给线圈并联一个阻尼电阻，模拟实际消弧线圈的有功损耗和阻尼率，取值在数百欧到数千欧之间。

还有一个非常隐蔽的浮地问题：用三个单相AC电源搭浮星形中性点时，如果不接地、不接消弧线圈，整个对地网络处在悬浮状态，某些版本会报矩阵奇异或者零序求解失败。解决方法是给中性点对地接一个大电阻，比如1MΩ，作为数值参考点，对工频计算结果的影响小到可以忽略。这个电阻只起数值稳定作用，不要把它当成真实接地通路。

另外，故障电阻设置也经常出问题。过渡电阻一旦超过几百千欧，故障特征被削弱得几乎看不见，3U0抬升不明显，波形上看起来像“没故障”。做故障仿真时，金属性接地用0.001Ω，模拟过渡电阻先从1Ω、10Ω、100Ω这几个量级试，对比零序量的变化趋势，不要一上来就填一个非常大的值。

5. 最后再分享一点我的实操体会

这类仿真的价值不只是把波形跑出来，更在于模型给了你一个可以反复“做实验”的平台。建议你不要只盯着Scope里的波形看，而是把关键数据导出来，比如用To Workspace把3U0、故障线路3I0、非故障线路3I0存成MATLAB变量，然后自己画幅值对比和相位差。当你亲眼看到同一套参数下，仅仅因为中性点接了一个电感，零序电流分布就从“对比明显”变成“几乎没差别”，比看理论书理解得透彻得多。

做完基础仿真之后，可以继续往几个方向扩展：一是验证不同的选线算法，比如暂态首半波法、五次谐波法、有功分量法，看看各自在消弧线圈接地系统里的表现；二是把架空线换成电缆，电容电流从几安变成几十安甚至上百安，消弧线圈容量和补偿效果又会不一样；三是加上电弧模型，模拟接地电弧的间歇性重燃，研究弧光接地过电压。这些内容都是在现在这个模型骨架上加细节就能实现的。做仿真这件事，最怕的就是只会按按钮、不会分析结果，所以就算模型再小，也建议每个关键波形都拿理论推导验证一遍再往下走。