10kV配电网的单相接地故障占了整个配电系统故障的七八成。为什么大家都盯着这个故障类型不放因为它牵着一个非常关键的设计决策系统中性点到底怎么接地。同一个故障放在中性点不接地系统里可能只是绝缘报警系统带着故障继续跑放在经消弧线圈接地的系统里故障点电弧很快熄灭线路不用停电但如果放任不管电弧重燃、过电压升高最终会从单相接地发展成相间短路直接扩大事故范围。我这次在Simulink里把这两种接地方式的配电系统都搭了出来用同一个单相金属性接地故障场景做对比仿真。从线路电容怎么折算、消弧线圈电感怎么算到模块参数在哪里填、波形怎么判读都按实际操作顺序记了下来。中间踩了几个坑比如变压器中性点引不出来、故障相电压相位不对、消弧线圈补偿前后波形看不出差别这一类问题也一并写了排查思路。如果你正在做配网仿真或者刚接触Simulink里的电力系统电磁暂态仿真这篇可以直接当操作参考。1. 仿真研究的落脚点中性点接地方式到底在解决什么问题1.1 先从配电网的故障特征说起配电系统的故障类型里单相接地占比极高原因也不难理解架空线裸露在外面雷击、树枝、鸟害、绝缘子污闪绝大多数故障都发生在相对地之间。电缆线路虽然封闭但绝缘受潮、外力破坏造成单相接地的情况同样常见。这类故障有两个显著特点第一故障电流往往不算大特别是中性点不接地系统里故障点流过的不是金属性短路电流而是由线路对地电容决定的容性电流通常只有几安到几十安第二很多单相接地属于瞬时性故障电弧本身不稳定一会儿燃一会儿灭。所以工程上一直有个矛盾——如果一发生单相接地就跳闸瞬时性故障也会变成永久性停电供电可靠性指标会非常难看。但如果不处理故障点电弧可能持续燃烧造成设备烧损甚至发展为两相短路。中性点接地方式就是用来平衡这个矛盾的。不接地也好经消弧线圈接地也好它们的共同思路是允许系统在单相接地故障下继续运行一段时间让运行人员有时间排查和转移负荷区别在于消弧线圈接地系统会主动抑制故障点电弧而不是单纯“硬扛”。1.2 两条完全不同的技术路线中性点不接地系统从名字就能看出来中性点和地之间没有任何人为的低阻通路。正常运行的时候三相对地电容电流对称中性点电位近似为零。发生单相接地以后故障相对地电压被压到接近零另外两相对地电压升高到线电压水平中性点电位则被迫位移到接近故障前相电压的数值。此时故障点流过的是电容电流只要这个电流足够小电弧就可能自行熄灭系统也能继续运行。问题是如果线路电容电流偏大比如以电缆为主的城市配网故障电流可能达到几十安培电弧很难熄灭绝缘薄弱点会一个接一个被打穿。消弧线圈接地系统针对的正是这个问题。它的做法是在变压器中性点接入一个电感线圈发生单相接地时线圈在相电压作用下产生一个感性电流方向和容性电流相反在故障点相互抵消使残流降到很小。电弧在电流过零时熄灭后由于恢复电压上升速度变慢不容易发生重燃。这条技术路线在工程上极其成熟城市配网、矿区电网里大量采用。仿真的价值就在于你不需要真实搭一套10kV系统就能把这两种接地方式在故障瞬间的电压电流行为完整复现出来用来验证规律、优化参数甚至做培训演示。2. 故障瞬间的物理图景电容电流和补偿电流是怎么算出来的2.1 中性点不接地故障电流是电容给的建模之前建议先把物理过程吃透。对中性点不接地系统发生A相金属性接地时A相电位被强制拉低到地电位。如果忽略接地电阻可以认为A相电压瞬间变成0B相和C相的对地电压从原来的相电压上升到线电压。中性点也不再停在原来的电位上而是发生了位移位移电压的大小约等于故障前的相电压方向与故障相电压相反。这时候故障点流过的电流是B相和C相对地电容放电、充电形成的电流而不是电源提供的短路电流。单相接地电容电流的标准工程计算式是IC 3 × ω × C0 × U_ph其中C0是系统的每相对地电容U_ph是相电压有效值ω是工频角频率。注意这里用的是每相对地电容不要跟三相总电容搞混。举个具体例子。10kV系统相电压有效值约5.77kV假设两条出线折算下来的每相对地电容C01.6μF那么IC 3 × 314 × 1.6×10⁻⁶ × 5770 ≈ 8.7A这个数值非常典型。架空线为主的系统电容电流也就几安电缆线路多了以后电容电流轻松超过20A这时候再用不接地方式就很不合理。如果你在仿真里发现故障电流和手算值对不上不要急着怀疑公式先检查线路模型的单位长度电容参数是不是设置成了三相总和或者线路长度填错了。后面会专门讲这个排查点。2.2 消弧线圈用反向电流把故障点电流顶回去消弧线圈接在中性点和地之间。单相接地时中性点电位变成接近相电压消弧线圈两端承受的就是这个电压。于是线圈中流过感性电流IL U_ph / (ωL)这个电流滞后电压90°而电容电流超前电压90°两者在故障点的相位正好相反。假设消弧线圈电流大小刚好等于系统电容电流故障点的总电流就只剩一个很小的差值电弧自然容易熄灭。但这里有个工程上必须避开的坑不能让消弧线圈刚好工作在“全补偿”状态也就是ILIC。全补偿时系统恰好处于工频谐振点附近如果中性点存在不对称电压会产生很高的串联谐振过电压危及绝缘。所以实际运行都采用过补偿让IL比IC大一点点残流呈感性。衡量这个“大一点点”的指标是脱谐度v (IC - IL) / IC过补偿时ILICv是负值但大家习惯直接说“过补偿8%”。一般过补偿脱谐度取5%到10%既能让残流足够小又远离谐振点。对应上面的算例如果取过补偿8%消弧线圈电流大约要做到9.4A电感值代入公式L U_ph / (ω × IL) 5770 / (314 × 9.4) ≈ 1.95H这个电感值对应的线圈容量也就在100kVA以内对一台10kV配电系统来说很常见。2.3 过渡电阻的影响先想清楚实际接地故障很少是纯金属性的很多情况下会有过渡电阻比如树枝接触、杆塔经土壤电阻接地。过渡电阻的存在会让故障相电压降不到零中性点位移电压也相应减小零序电压判据的灵敏度会下降。仿真建模时建议先把金属性接地接地电阻设0.001Ω跑通再换10Ω、100Ω甚至1000Ω做对比。这样你对“故障电压和电流如何随过渡电阻变化”会有一个很直观的认识而不是只看一堆标幺值。3. Simulink建模全流程从模块拖拽到参数设置3.1 模型拓扑五大模块块一条主线我搭的模型结构非常简单你不必一开始就追求高精度先把主回路搭出来、跑出能解释物理规律的波形再逐步细化。整条链路的顺序是110kV电源等值系统 → 降压变压器 → 10kV母线 → 两条出线 → 末端负荷故障点放在其中一条出线的末端测量点放在母线侧和故障支路。模块选择上用到的都是Simscape Electrical里Specialized Power Systems库的标准件功能推荐模块备注电源Three-Phase Source设为110kV等值电源短路容量可以默认取大值变压器Three-Phase Transformer (Two Windings)Dyn11接线低压侧中性点引出线路Three-Phase PI Section Line用π型集中参数模型单位长度参数填入负荷Three-Phase RLC Load每条出线各带一组故障Three-Phase Fault选择A相接地可设置过渡电阻消弧线圈Single-Phase Series RLC Branch接在变压器低压侧中性点与地之间测量Three-Phase V-I Measurement配合Scope观察母线电压、线路电流仿真环境Powergui没有这个模块模型根本跑不起来必须强调一下Powergui。很多新手第一次搭Simscape Electrical模型拖了一堆三相模块点运行直接报错原因就是模型里没有放Powergui。它不是可有可无的“工具箱入口”而是整个电气仿真系统的求解环境必须先拖一个进去。3.2 参数表这些数可以直接抄我用的一组典型参数如下线路参数是参考配电线路常用典型值取的你完全可以根据自己的场景改电源与变压器参数设置值电源电压110kV50Hz三相对称变压器额定容量20MVA变压器变比110kV / 10.5kV变压器接线Dyn11低压侧为Y接且中性点引出变压器短路电阻0.002 pu变压器短路电感0.08 pu线路参数每公里参数线路1架空线线路2电缆长度10km5km正序电阻R10.27 Ω/km0.1 Ω/km正序电感L11.02 mH/km0.25 mH/km正序电容C10.009 μF/km0.3 μF/km零序电阻R00.35 Ω/km0.15 Ω/km零序电感L04.13 mH/km0.35 mH/km零序电容C00.006 μF/km0.2 μF/km两条线路并联后折算出来的每相对地电容大约是1.6μF这个数值就是前面手算消弧线圈参数的来源。你把同样的参数填进自己的模型跑出来的故障电流就会落在8到9A这个量级和理论值对得上。负荷每条出线带一组恒阻抗负荷有功1MW、无功0.2Mvar电压设为10kV。故障发生在负荷侧故障期间这部分负荷会受到严重扰动这是正常的。3.3 消弧线圈参数计算用一个脚本直接算上面手算过一遍实际建模时我用MATLAB脚本算既方便改参数也能避免手误C0 1.6e-6; % 每相对地电容, F U_ph 10e3 / sqrt(3); % 相电压有效值, V omega 2 * pi * 50; % 工频角频率 IC 3 * omega * C0 * U_ph; % 单相接地电容电流, A v 0.08; % 过补偿脱谐度 8% IL (1 v) * IC; % 消弧线圈电流, A L U_ph / (omega * IL); % 消弧线圈电感, H fprintf(IC %.2f A\n, IC); fprintf(IL %.2f A\n, IL); fprintf(L %.3f H\n, L);运行后得到IC≈8.7AIL≈9.4AL≈1.95H。消弧线圈的电阻不能设成0实际线圈总有有功损耗我取了15Ω大约相当于线圈品质因数Q在40左右比较接近真实设备水平。3.4 故障模块设置和两种工况切换故障模块双击后在Fault selection里勾选Phase A故障类型选Phase-to-ground。Transition time有两种方式内部定时或外部信号触发。我建议用内部定时设0.1s开始故障0.25s清除故障整个仿真时长0.3s。这样你能同时看到故障前稳态、故障期间暂态、故障切除后的恢复过程三段波形。故障电阻Ron我设0.001Ω模拟金属性接地。Snubber resistance保持默认的1e6Ω即可它只影响数值求解的收敛性不影响稳态结果。切换两种接地方式的关键点在变压器低压侧绕组中性点不接地系统变压器低压侧设为Y接法中性点不引出、不接地。此时消弧线圈支路直接断开即可。消弧线圈接地系统变压器低压侧改为Yg接法把消弧线圈接到中性点端子和地之间。这里有一个容易卡住的地方。有些版本的变压器模块选Y接法时根本没有中性点端子可以外接设备而选Yg接法时中性点可能是内部直接接地的。不同版本UI不一样如果你找不到外接中性点端子的选项一个笨但可靠的办法是用三个单相变压器搭成三相变压器组低压侧绕组按Yg连接中性点自然就暴露出来了。反正仿真是为了看规律变压器铁芯饱和、磁路耦合这些细节在这个问题里影响不大。3.5 求解器设置为什么我推荐ode23tbSimulink默认的求解器对纯控制类模型很友好但电力电子和电力系统电磁暂态模型里存在大量刚性环节比如变压器的漏感、线路的分布电容、故障切换瞬间的快速暂态用默认变步长ode45很容易出现步长不受控、仿真发散的情况。我建议在求解器设置里做三件事第一求解器选ode23tb这个算法对刚性系统非常稳收敛性比ode45好得多。如果你后面把模型扩展成含逆变器、电力电子器件的系统大概率还是要回到ode23tb。第二最大步长设1e-4秒。不要小看这个限制电气暂态的典型时间常数在毫秒级步长太大不仅波形难看还可能漏掉关键尖峰。第三仿真时长设0.3秒故障从0.1秒开始。前0.1秒是留给系统进入稳态的这个时间足够变压器励磁涌流衰减、线路电容充放电完成。另外Powergui里有一个Load Flow工具可以在仿真前对系统进行潮流初始化。我习惯先把所有母线电压设成额定值运行一次Load Flow再接着跑暂态仿真。这样启动冲击会小很多波形前段也不会因为初始状态没匹配而出现诡异的尖峰。4. 仿真结果解读波形里能看到什么4.1 母线三相电压一眼判断故障类型先看不接地系统的母线电压波形。故障前三相电压对称幅值峰值在8.16kV左右对应10kV系统的相电压幅值10/√3×√2。故障发生后A相电压迅速跌到接近0B相和C相电压升高到接近14.1kV这个数值正好是线电压的峰值。这个特征非常清晰单相接地后非故障相电压必然升高到线电压水平。消弧线圈接地系统的电压波形与不接地系统几乎一样。原因很容易理解消弧线圈只改变了故障点的电流分布并没有改变“故障相被拉到地电位、非故障相升高√3倍”这个基本物理事实。所以如果你看到两种系统的电压波形长得差不多不要觉得模型错了这是正常的。4.2 故障点电流补偿效果就看这一条波形不接地系统的故障点电流波形是一个稳定的容性电流。有效值大概在8到9A峰值12A左右和前面手算的IC8.7A完全吻合。它的相位超前电压90度示波器上看起来就是“电压怎么变电流提前四分之一周期跟着变”。切换到消弧线圈接地系统后同样是0.1秒发生A相接地故障点电流波形会发生肉眼可见的变化波形整体幅度大幅下降只剩下一个很小的残流基本就是过补偿脱谐度对应的那部分电流也就是IL与IC的差值大约0.7A左右。除此之外还会带一点线圈电阻产生的有功分量和谐波分量。这才是消弧线圈真正的作用现场它没有改变系统电压但把故障点的能量输入降到了几乎可以忽略的程度。电弧失去维持电流自然在过零时熄灭不会重燃。4.3 中性点位移电压和3U0两种系统中性点的对地电压在故障期间都大约等于相电压。这个量在工程上对应的是零序电压判据做法一般是取三相电压矢量和的三分之一。仿真里你可以直接测量中性点对地电压也可以在母线电压测量模块后面用算术模块把三相电压加在一起得到3U0波形。故障前3U0接近0故障后3U0跳到约3倍相电压也就是相电压有效值的3倍左右对应到电压互感器二次开口三角绕组就是100V定值。这就是继电保护里“单相接地报警”的核心判据仿真波形里看得一清二楚。4.4 两者放到一张表里对比观测量中性点不接地经消弧线圈接地故障相电压降到约0降到约0非故障相电压升到线电压升到线电压中性点位移电压约等于相电压约等于相电压故障点电流约IC容性主导残流很小过补偿时呈感性电弧自熄能力取决于IC大小IC大时难自熄残流被补偿电弧容易熄灭故障后处理策略可短时继续运行加强监视可短时继续运行残流小更安全这张表基本把一个问题的核心说透了两者的电压行为是相同的差别集中在故障点电流和灭弧能力上。仿真里如果你只想看一个波形就盯住故障点电流那条线它最能说明补偿效果。5. 跑仿真必踩的坑与排查清单5.1 模型跑不通、仿真直接发散电气模型仿发散最常见的原因有三个没放Powergui、求解器不合适、步长约束太松。Powergui缺失的问题前面提过解法就是拖一个进模型。求解器方面不要用ode45换成ode23tb。步长方面最大步长不要超过1e-4如果还是发散继续降到5e-5。还有一种情况是变压器或线路模块里的初始参数填得自相矛盾比如容量的单位、二次侧线电压和相电压的关系弄错。建议先把电源、变压器单独跑一遍确认空载电压正常再往下加线路和负荷层层排查比一次性搭完全部再找问题效率高得多。5.2 故障相电压没有降到0出现这个现象先看故障模块里勾选的故障相是不是A相很多人默认只勾了Phase A但实际模型中测量的却是B相电压这种低级错误最容易漏掉。再看故障电阻Ron。金属性接地时Ron应该设得很小0.001Ω合适如果你设了10Ω甚至更大故障相电压自然降不到0非故障相电压也升不到线电压整个波形看起来会被“削掉一截”。还要留意故障时间设置。如果故障开始时间晚于仿真结束时间或者故障持续窗口太短波形上可能看不到完整的故障段。我一般把故障开始时间设为0.1s仿真时长0.3s留足观察窗口。5.3 消弧线圈看不出来补偿效果如果消弧线圈接地系统的故障点电流和不接地系统差不多大先怀疑电感值有没有填对。用脚本算出L之后直接填进Single-Phase Series RLC Branch的Inductance字段单位是H不要填成mH差了1000倍波形完全没法看。再检查补偿方向。我强调过要过补偿也就是IL大于IC。如果你设置的是欠补偿残流虽然也小了但电容电流的主导地位更强波形相位仍然是容性的。从波形上是能区分出来的残流滞后电压90度说明是过补偿残流超前电压90度说明是欠补偿。还有一个容易忽略的点消弧线圈的有功电阻。真实线圈有损耗电阻对应的是有功电流分量这个分量不能补偿会一直残留在故障点电流里。如果电阻设得太大残流会被有功分量主导看起来“补不干净”。我取15Ω是合理的但如果你发现残流超过好几安培检查一下是不是把线圈电阻设成了几百欧。5.4 波形前段有巨大尖峰或者启动浪涌启动尖峰几乎都是初始状态没设置好造成的。变压器在第一个周波里会经历励磁涌流线路电容也会和电源发生充电暂态。这些暂态不是故障引起的却会干扰你对故障波形的判断。解决办法是先运行Powergui里的Load Flow工具完成初始化然后再跑仿真。如果还是觉得前段波形太乱可以把故障时间往后推到0.2s前0.2秒完全给系统做稳态过渡。波形依然不干净的话把最大步长再缩小一个数量级很多时候“看起来奇怪”只是步长太大导致的高频振荡丢失。5.5 零序电压和理论值对不上零序电压的测量在Simulink里可以直接把三相电压通过一个算术模块相加得到的和除以3就是零序电压不除以3就是3U0。如果你用了序列分析模块注意输出的是对称分量单位可能是标幺值或V要确认清楚。故障期间3U0接近3倍相电压这是单相接地的特征。如果测出来的值只有相电压那么小大概率是你加的是线电压而不是相电压或者测量点放在了负载后面而不是母线侧。6. 写在最后几个让模型更“抗造”的小习惯仿真模型做完了我再分享几个自己用下来的习惯。第一把消弧线圈的电感L、电阻R做成模型里的可调参数不要每次直接改模块内部数值。我习惯用MATLAB脚本统一计算再通过set_param批量写回模型。这样你切换脱谐度、对比不同补偿效果时只需要改脚本里的v值一键完成不用手动翻模块参数。第二波形不要只盯着Scope看。Scope适合快速检查但如果要写报告或论文我会把关键量通过To Workspace模块导到MATLAB工作区再统一用plot函数整理成同一坐标系下的对比图。这样图面干净、线型统一也方便做局部放大。第三做对比仿真时务必保证除了“中性点是否接消弧线圈”之外所有其他参数完全一致。这个要求看似废话但实际操作中很容易因为改了线路长度、负荷大小而不自知最后得出的结论失真。我一般会建一个参数表文件每次改动都登记跑完一组对比之后回头检查一遍确认变量确实只有一个。10kV系统的中性点接地问题说到底是“故障电流控制”和“供电可靠性”之间的平衡。Simulink这个模型虽然简化了不少工程细节但已经把两种接地方式最核心的行为差异还原得很清楚。你如果把它改造成带多条电缆出线的模型电容电流变大消弧线圈的补偿效果会更加明显值得一试。