继电保护这块小电流接地系统的单相接地故障仿真几乎是每个做配电系统仿真的人都会碰到的课题。我自己最初做这个方向的时候最难的不是保护逻辑本身而是“模型到底该怎么搭”——中性点不接地和经消弧线圈接地这两种方式在Simulink里的接法只差一个支路但波形特征和保护判据完全不同。这篇内容我把整套仿真模型的搭建思路、消弧线圈参数计算、故障设置细节以及常见坑都整理出来正在做课程设计、毕业设计或者工程预研的朋友可以直接照着上手省掉自己反复试错的时间。1. 仿真之前先搞懂小电流接地系统在仿什么1.1 两种接地方式原理差异直接决定模型结构小电流接地系统通俗讲就是中性点不直接接地的系统常见的有中性点不接地和经消弧线圈接地两种。所谓“小电流”指的是单相接地故障时流过故障点的电流不是靠电源短路电流维持的而是靠线路对地电容形成的容性电流数值通常只有几安到几十安。中性点不接地系统工程上也叫IT系统电网和大地之间只有三相对地电容这条“弱联系”。发生单相金属性接地时A相电压被拉到接近0B、C两相电压升高到线电压也就是从原来的5.77kV涨到10kV左右以10kV系统为例三相线电压仍然对称所以系统还能带故障运行一到两个小时。但非故障相电压抬高对绝缘是个考验。经消弧线圈接地系统就是在中性点与大地之间串入一个电感线圈。单相接地时消弧线圈提供一个感性电流去抵消三相对地电容产生的容性电流故障点残流变小电弧更容易自行熄灭。这样做最大的好处是减少弧光接地过电压代价就是故障特征变得不明显继电保护选线变得很困难。从仿真建模的角度看两种系统的区别本质上就在一根线上中性点对地之间接不接消弧线圈。这给了我们一个非常好的便利——搭好一套主电路通过切换中性点支路就能跑出两套完整的结果。1.2 消弧线圈补偿的参数计算消弧线圈的补偿效果核心是看电感L选多大。推导并不复杂单相接地后中性点电压偏移到相电压Uφ零序等效网络里三相对地电容3C0和消弧线圈L是并联的。完全补偿时要求电感电流等于电容电流Uφ / (ωL) 3ωC0Uφ两边消掉Uφ得到L 1 / (3ω²C0)我以一个10kV、50Hz、出线10km的简单系统为例。架空线对地电容取经验值C01μ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%时IL1.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 × 3峰值8165V相电压峰值50Hz相位0°/-120°/120°线路阻抗Three-Phase Series RLC Branch或三个单相Series RLC BranchR1ΩL0.01H对地电容三个单相Parallel RLC Branch分别接A/B/C相与地之间C1μFR、L按需故障模块Three-Phase FaultA相接地故障电阻0.001Ω测量模块Three-Phase V-I Measurement输出Vabc、Iabc消弧线圈Series RLC Branch接在中性点与地之间L3.07H阻尼R≈200Ω~2kΩ数值参考电阻单相Resistor1MΩ不接地系统中性点对地求解器设置PowerguiDiscreteTs50e-6s三个单相AC Voltage Source的幅值必须填相电压峰值不是线电压有效值。10kV系统的相电压有效值是10k/√35774V对应峰值是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之前三相电压对称每相相电压有效值5774V0.1s A相接地后A相电压瞬间塌到接近0B、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是一个非常可靠的启动量。但零序电流幅值保护在经消弧线圈接地系统里就不太可靠了。所以实际工程里中性点不接地系统可以尝试零序电流比幅比相经消弧线圈接地系统更多依靠零序电压启动再配合暂态量选线、有功分量法或者注入信号法。做仿真验证时可以在模型中同时输出故障线路和非故障线路的零序电流相位自己画一下矢量关系理解会深很多。两种系统主要电气量对比如下对比项中性点不接地经消弧线圈接地故障相电压接近0接近0非故障相电压升高到线电压升高到线电压3U0有效值约3倍相电压约3倍相电压故障点残流约5.4A容性小于1A且呈感性故障线路零序电流明显大于非故障线路与非故障线路差别很小常规比幅选线可行基本失效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变量然后自己画幅值对比和相位差。当你亲眼看到同一套参数下仅仅因为中性点接了一个电感零序电流分布就从“对比明显”变成“几乎没差别”比看理论书理解得透彻得多。做完基础仿真之后可以继续往几个方向扩展一是验证不同的选线算法比如暂态首半波法、五次谐波法、有功分量法看看各自在消弧线圈接地系统里的表现二是把架空线换成电缆电容电流从几安变成几十安甚至上百安消弧线圈容量和补偿效果又会不一样三是加上电弧模型模拟接地电弧的间歇性重燃研究弧光接地过电压。这些内容都是在现在这个模型骨架上加细节就能实现的。做仿真这件事最怕的就是只会按按钮、不会分析结果所以就算模型再小也建议每个关键波形都拿理论推导验证一遍再往下走。