☰
10kV配电网单相接地故障仿真:中性点不接地与消弧线圈补偿对比
2026/10/6 4:59:03 网站建设 项目流程

10kV配电网的单相接地故障,占了整个配电系统故障的七八成。为什么大家都盯着这个故障类型不放?因为它牵着一个非常关键的设计决策:系统中性点到底怎么接地。同一个故障,放在中性点不接地系统里,可能只是绝缘报警,系统带着故障继续跑;放在经消弧线圈接地的系统里,故障点电弧很快熄灭,线路不用停电;但如果放任不管,电弧重燃、过电压升高,最终会从单相接地发展成相间短路,直接扩大事故范围。

我这次在Simulink里把这两种接地方式的配电系统都搭了出来,用同一个单相金属性接地故障场景做对比仿真。从线路电容怎么折算、消弧线圈电感怎么算,到模块参数在哪里填、波形怎么判读,都按实际操作顺序记了下来。中间踩了几个坑,比如变压器中性点引不出来、故障相电压相位不对、消弧线圈补偿前后波形看不出差别这一类问题,也一并写了排查思路。如果你正在做配网仿真,或者刚接触Simulink里的电力系统电磁暂态仿真,这篇可以直接当操作参考。

1. 仿真研究的落脚点:中性点接地方式到底在解决什么问题

1.1 先从配电网的故障特征说起

配电系统的故障类型里,单相接地占比极高,原因也不难理解:架空线裸露在外面,雷击、树枝、鸟害、绝缘子污闪,绝大多数故障都发生在相对地之间。电缆线路虽然封闭,但绝缘受潮、外力破坏造成单相接地的情况同样常见。

这类故障有两个显著特点:第一,故障电流往往不算大,特别是中性点不接地系统里,故障点流过的不是金属性短路电流,而是由线路对地电容决定的容性电流,通常只有几安到几十安;第二,很多单相接地属于瞬时性故障,电弧本身不稳定,一会儿燃一会儿灭。所以工程上一直有个矛盾——如果一发生单相接地就跳闸,瞬时性故障也会变成永久性停电,供电可靠性指标会非常难看。但如果不处理,故障点电弧可能持续燃烧,造成设备烧损,甚至发展为两相短路。

中性点接地方式就是用来平衡这个矛盾的。不接地也好,经消弧线圈接地也好,它们的共同思路是:允许系统在单相接地故障下继续运行一段时间,让运行人员有时间排查和转移负荷;区别在于,消弧线圈接地系统会主动抑制故障点电弧,而不是单纯“硬扛”。

1.2 两条完全不同的技术路线

中性点不接地系统,从名字就能看出来,中性点和地之间没有任何人为的低阻通路。正常运行的时候,三相对地电容电流对称,中性点电位近似为零。发生单相接地以后,故障相对地电压被压到接近零,另外两相对地电压升高到线电压水平,中性点电位则被迫位移到接近故障前相电压的数值。此时故障点流过的是电容电流,只要这个电流足够小,电弧就可能自行熄灭,系统也能继续运行。问题是,如果线路电容电流偏大,比如以电缆为主的城市配网,故障电流可能达到几十安培,电弧很难熄灭,绝缘薄弱点会一个接一个被打穿。

消弧线圈接地系统针对的正是这个问题。它的做法是在变压器中性点接入一个电感线圈,发生单相接地时,线圈在相电压作用下产生一个感性电流,方向和容性电流相反,在故障点相互抵消,使残流降到很小。电弧在电流过零时熄灭后,由于恢复电压上升速度变慢,不容易发生重燃。

这条技术路线在工程上极其成熟,城市配网、矿区电网里大量采用。仿真的价值就在于:你不需要真实搭一套10kV系统,就能把这两种接地方式在故障瞬间的电压电流行为完整复现出来,用来验证规律、优化参数甚至做培训演示。

2. 故障瞬间的物理图景:电容电流和补偿电流是怎么算出来的

2.1 中性点不接地:故障电流是电容给的

建模之前,建议先把物理过程吃透。对中性点不接地系统,发生A相金属性接地时,A相电位被强制拉低到地电位。如果忽略接地电阻,可以认为A相电压瞬间变成0,B相和C相的对地电压从原来的相电压上升到线电压。中性点也不再停在原来的电位上,而是发生了位移,位移电压的大小约等于故障前的相电压,方向与故障相电压相反。

这时候故障点流过的电流,是B相和C相对地电容放电、充电形成的电流,而不是电源提供的短路电流。单相接地电容电流的标准工程计算式是:

IC = 3 × ω × C0 × U_ph

其中C0是系统的每相对地电容,U_ph是相电压有效值,ω是工频角频率。注意这里用的是每相对地电容,不要跟三相总电容搞混。

举个具体例子。10kV系统,相电压有效值约5.77kV,假设两条出线折算下来的每相对地电容C0=1.6μF,那么:

IC = 3 × 314 × 1.6×10⁻⁶ × 5770 ≈ 8.7A

这个数值非常典型。架空线为主的系统,电容电流也就几安;电缆线路多了以后,电容电流轻松超过20A,这时候再用不接地方式就很不合理。

如果你在仿真里发现故障电流和手算值对不上,不要急着怀疑公式,先检查线路模型的单位长度电容参数是不是设置成了三相总和,或者线路长度填错了。后面会专门讲这个排查点。

2.2 消弧线圈:用反向电流把故障点电流顶回去

消弧线圈接在中性点和地之间。单相接地时,中性点电位变成接近相电压,消弧线圈两端承受的就是这个电压。于是线圈中流过感性电流:

IL = U_ph / (ωL)

这个电流滞后电压90°,而电容电流超前电压90°,两者在故障点的相位正好相反。假设消弧线圈电流大小刚好等于系统电容电流,故障点的总电流就只剩一个很小的差值,电弧自然容易熄灭。

但这里有个工程上必须避开的坑:不能让消弧线圈刚好工作在“全补偿”状态,也就是IL=IC。全补偿时,系统恰好处于工频谐振点附近,如果中性点存在不对称电压,会产生很高的串联谐振过电压,危及绝缘。所以实际运行都采用过补偿,让IL比IC大一点点,残流呈感性。衡量这个“大一点点”的指标是脱谐度:

v = (IC - IL) / IC

过补偿时IL>IC,v是负值,但大家习惯直接说“过补偿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 参数表:这些数可以直接抄

我用的一组典型参数如下,线路参数是参考配电线路常用典型值取的,你完全可以根据自己的场景改:

电源与变压器:

参数设置值
电源电压110kV,50Hz,三相对称
变压器额定容量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.7A,IL≈9.4A,L≈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 求解器设置:为什么我推荐ode23tb

Simulink默认的求解器对纯控制类模型很友好,但电力电子和电力系统电磁暂态模型里存在大量刚性环节,比如变压器的漏感、线路的分布电容、故障切换瞬间的快速暂态,用默认变步长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相电压迅速跌到接近0,B相和C相电压升高到接近14.1kV,这个数值正好是线电压的峰值。这个特征非常清晰:单相接地后,非故障相电压必然升高到线电压水平。

消弧线圈接地系统的电压波形与不接地系统几乎一样。原因很容易理解,消弧线圈只改变了故障点的电流分布,并没有改变“故障相被拉到地电位、非故障相升高√3倍”这个基本物理事实。所以如果你看到两种系统的电压波形长得差不多,不要觉得模型错了,这是正常的。

4.2 故障点电流:补偿效果就看这一条波形

不接地系统的故障点电流波形是一个稳定的容性电流。有效值大概在8到9A,峰值12A左右,和前面手算的IC=8.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这个模型虽然简化了不少工程细节,但已经把两种接地方式最核心的行为差异还原得很清楚。你如果把它改造成带多条电缆出线的模型,电容电流变大,消弧线圈的补偿效果会更加明显,值得一试。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询