小电流接地系统单相接地故障Simulink仿真建模与参数计算
2026/9/15 2:00:14 网站建设 项目流程

写这个仿真之前,我先把丑话说在前面:小电流接地系统的故障仿真,难点从来不在Simulink操作本身,而在于你是否真正理解“中性点不接地”和“经消弧线圈接地”这两种方式在故障瞬间到底发生了什么。模型只是把你的理解翻译成波形而已。这个仿真模型我前前后后搭过好几版,从最开始的三相电压直接飙到线电压、消弧线圈参数怎么调都不对,到最后能稳定复现单相接地故障的暂态过程,中间踩了不少坑。这篇博文就把完整的建模思路、参数计算过程、常见问题排查一次说清楚,给后面要做这类仿真的人省点时间。不管是做毕业设计、继电保护整定计算,还是想验证选线算法,这套模型都值得你从头搭一遍。

1. 小电流接地系统仿真:先弄懂这两种接地方式的“脾气”

1.1 为什么要单独研究小电流接地系统

电力系统的中性点接地方式,说白了就是变压器星形侧中性点和大地之间怎么处理。在6kV到66kV这个电压等级的配电网里,大量采用中性点不接地或经消弧线圈接地的方式,统称为小电流接地系统。和110kV及以上大电流接地系统最大的区别在于:发生单相接地故障时,故障电流不是短路电流那种几千安培的“暴力”级别,而是由线路对地电容决定的电容电流,通常只有几安培到几十安培。

这个特点导致了一个非常特殊的运行规则:单相接地故障发生后,三相之间的线电压仍然保持对称,负荷可以继续运行,系统允许带故障运行1到2小时。这对供电可靠性是很大的优势,但代价是保护装置很难通过电流大小来判断故障线路。做仿真模型的核心目的,就是要把这种“小电流”故障下电压电流的细微变化提取出来,为选线、消弧、定值整定提供依据。

1.2 中性点不接地与经消弧线圈接地的本质差异

很多初学者容易把这两种方式混为一谈,其实它们的故障特征有本质区别。中性点不接地系统发生单相金属性接地时,故障相电压降为零,非故障相电压升高到线电压,零序电压升高到相电压水平。这个特征可以用来看故障是否发生、故障相是哪一相,但无法区分故障线路。

经消弧线圈接地系统则是在中性点和大地之间串入一个电感线圈。单相接地时,消弧线圈在零序电压作用下产生感性电流,这个电流恰好和系统对地电容电流方向相反,在接地点相互抵消,使得故障点残流大大减小,电弧更容易自行熄灭。两种方式的本质区别在于:一个是“硬扛”电容电流,一个是“主动补偿”电容电流。仿真模型里最核心的差异也就在这一处——有没有那个消弧线圈,以及它的参数选得对不对。

对比项中性点不接地经消弧线圈接地
故障相电压金属接地时降为零降为零,但受补偿状态影响有偏移
非故障相电压升高到线电压升高到线电压附近
接地故障电流系统对地电容电流补偿后的残流,明显减小
电弧自熄能力弱,易产生间歇性电弧强,过补偿方式下电弧快速熄灭
保护选线难度高,需利用零序电流方向更高,残流太小,方向特征不明显

2. Simulink建模仿真:从模块选型到参数设定的完整思路

2.1 模型整体结构怎么搭

在Simulink里建电力系统仿真模型,我强烈建议你从Simscape Electrical库(早年叫SimPowerSystems)里拖模块,不要自己用数学运算模块拼电路。你需要的核心模块包括:三相电源、变压器、三相PI型输电线路、三相故障模块、电压电流测量单元,以及Powergui模块。Powergui是仿真环境配置的核心,选择离散还是连续求解方式、采样步长多少,都靠它来控制。

整个模型的结构其实就是一条链路:三相电源输出经过升压变压器进入输电线路,线路末端接负荷,故障模块并联在某一相或某两相与地之间。中性点不接地和经消弧线圈接地两种方案,差别就在变压器中性点的处理上。不接地系统直接把中性点悬空或通过高阻抗接地;经消弧线圈系统则在中性点和地之间串入一个电感支路。模型搭好后,通过一个手动开关或参数切换,就能在两种接地方式间自由切换对比。

2.2 选型时最容易出错的几个模块细节

第一个坑是线路模型。很多人用理想的Three-Phase Source直接带电阻负荷,这样仿出来电容电流基本为零,根本看不到小电流接地系统的故障特征。正确做法是用Three-Phase PI Section Line模块,并且在参数里把每相对地电容、相间电容都设进去。线路对地电容是这个仿真的灵魂,没有它,消弧线圈就是空中楼阁。

第二个坑是变压器中性点的引出。Three-Phase Transformer (Two Windings)模块默认是Y接法,但中性点不一定引出来了。你要在绕组连接方式里选择“YN”或“Yg”,这样才能把中性点端子引出到外部电路。另外,变压器模型的选择也很重要,我建议用“Three-phase transformer (two windings)”里的饱和模型还是线性模型要看研究目标。如果只是稳态故障分析,线性模型够用;如果要研究铁芯饱和引起的谐波问题,就得勾选饱和特性,输入饱和曲线数据。

第三个坑是故障模块。Three-Phase Fault模块默认可以设置单相、两相、三相故障,还可以设置故障发生时间和故障持续时间。这里注意设置好接地电阻值,金属性接地就填0.001欧姆左右,弧光接地就填几十到几百欧姆。故障时序用时间序列数组控制,比如[1/50, 1.1/50]表示在0.02秒发生故障、持续0.02秒后切除,这个时间设置要和仿真时长匹配。

2.3 消弧线圈模型怎么选

消弧线圈在Simulink里有几种建模方式:最简方式是用一个Series RLC Branch,在电感支路上设置合适的电感值和电阻值,电阻用来等效线圈的有功损耗,这个方案适用于研究基波补偿特性。进阶方式是用Single-Phase Transformer来模拟带铁芯电感的非线性特征,通过设置饱和特性让它在过电压时呈现非线性,这样的模型可以研究铁芯饱和、谐振等复杂问题。如果你只是为了验证消弧线圈对电容电流的补偿作用,RLC支路完全够用。

消弧线圈的补偿方式也要理解清楚。欠补偿是电感电流小于电容电流,系统在某种运行方式下可能接近全补偿,引发串联谐振过电压,工程上一般不采用。全补偿在理论上让接地残流为零,但实际上因为系统三相对地电容不对称,中性点位移电压会非常高,同样不可取。过补偿是让电感电流略大于电容电流,脱谐度控制在5%到10%之间,这是工程上最常用的方案。仿真模型里就用过补偿参数,我可以给你一组可直接用的参考值。

3. 以10kV系统为参考:故障仿真参数的详细计算过程

3.1 系统基础参数设定

为了让参数计算过程有具体的数字可以依据,我以一个典型的10kV中压配电网为例来演示。三相电源用10.5kV线电压有效值,频率50Hz,经过变压器升压到35kV或110kV再接线路,这里我们简化处理,直接让10.5kV电源经过一个变比1:1的隔离变压器进入10kV线路网络。变压器容量设10MVA,短路阻抗百分比设为6.5%,这些参数影响故障时的电压跌落程度。

线路模型用三段π型线路级联,每段长度10公里,单位长度正序电阻0.27Ω/km、正序电感1.2mH/km、正序电容0.009μF/km,零序参数适当调整。这里的电容值决定了系统总的接地电容电流,是消弧线圈参数设计的原始依据。负荷用三相并联RLC负荷模块,有功功率2MW,无功功率0.2Mvar。故障模块并联在A相和地之间,故障时间设置在0.3秒发生,持续0.1秒后切除。仿真总时长0.6秒,用离散求解器,步长50微秒,这样既能看清故障暂态过程,又不至于算得太慢。

3.2 对地电容电流估算与消弧线圈容量计算

在确定消弧线圈参数之前,必须先估算系统对地电容电流。每一相对地电容约为C0 = 0.009μF/km × 30km = 0.27μF。正常运行时每相电容电流为:

I_C = ω × C0 × U_φ = 2π × 50 × 0.27×10⁻⁶ × (10.5×10³/√3) ≈ 0.52A

单相接地时,流过故障点的总电容电流为:

3I_C = 3 × 0.52 ≈ 1.56A

这个数值看起来不大,但在真实配电网里线路众多,总电容电流很容易达到几十安培,电弧就不容易熄灭了。消弧线圈的作用就是提供一个大小相当但方向相反的感性电流来抵消它。消弧线圈电感值按全补偿计算:

L = 1 / (3ω²C0) = 1 / (3 × (2π×50)² × 0.27×10⁻⁶) ≈ 12.5H

但全补偿会引发串联谐振,实际采用过补偿,脱谐度取-8%左右,这时消弧线圈补偿电流为:

I_L = 3I_C / (1 - 0.08) ≈ 1.56 / 0.92 ≈ 1.70A

对应电感量需要向下调整,同时线圈电阻按感抗的3%到5%估算。例如感抗X_L = 2π×50×L,在过补偿条件下L约11.5H,感抗约3612Ω,电阻取100Ω左右。这样设置后故障点残流大约只有0.14A,消弧效果非常明显,仿真波形上能看到故障电流明显减小。

3.3 仿真时间和求解器选择

仿真时间和步长的选择直接影响波形质量和仿真速度。离散求解器步长我建议设在20到100微秒之间,步长太大,故障瞬间的高频暂态过程会被抹掉;步长太小,0.6秒的仿真跑起来会让人想砸电脑。如果你用的是连续求解器,ode23tb或ode15s这类刚性求解器比较适合电力电子和电力系统故障仿真,但连续模式下波形采样点也可能不够均匀,不利于后续数据分析,所以我个人更推荐离散模式。

4. 单相接地故障仿真结果怎么看:波形特征与数据提取

4.1 中性点不接地系统的故障波形特征

仿真跑完后,重点观察几个量的波形:三相相电压、零序电压3U0、故障点电流、各线路零序电流。在中性点不接地系统里,0.3秒A相金属性接地后,A相电压迅速跌落到接近零,B相和C相电压升高到线电压水平,近似10.5kV。零序电压同时跃升到相电压水平,3U0接近18.2kV的峰值。

这里你会看到故障相电压并非瞬间变成零,而是一个快速的暂态过程,大概持续几毫秒。原因是线路对地电容里的储能要释放出来,和系统电感形成高频振荡。这个高频暂态分量在真实的选线装置里很重要,很多暂态选线算法就是利用故障初始时刻的暂态电流和电压极性来判断故障线路的,所以仿真时千万不要只盯着稳态波形看,把故障发生时刻附近放大再仔细看。

4.2 经消弧线圈接地系统的补偿效果对比

切换到经消弧线圈接地模型,同样设置A相0.3秒接地,你会发现故障点电流大幅减小。以刚才那个参数为例,不接地时故障点电流约1.56A,接入过补偿消弧线圈后,残流降到0.1到0.2A左右。在波形图上最直观的感受是:故障电流的幅值几乎贴着横轴走,和正常工况没有明显区别。

同时你还会注意到一个细节,消弧线圈的接入改变了故障暂态过程的振荡频率和衰减速度。因为电感电流不能突变,故障初始阶段消弧线圈相当于开路,高频电容电流还在;几个毫秒后电感电流建立起来,才开始发挥补偿作用。所以暂态阶段和稳态阶段的故障特征并不一样,这对理解消弧线圈的实际动作过程非常有帮助。

4.3 用FFT分析故障特征分量

Simulink的Powergui里自带FFT分析工具,非常适合对波形做频谱分析。把故障后的电压电流信号导入FFT窗口,可以看到基波分量的幅值以及各次谐波含量。中性点不接地系统发生单相接地时,零序电流的基波幅值等于系统电容电流;而经消弧线圈接地时,零序电流基波被大幅抵消,但5次谐波等高频分量仍然存在。这个特征为谐波选线法提供了依据,也是我为什么建议你把FFT工具用起来的原因——只看时域波形是不够的,频域信息才是保护原理设计的“富矿”。

用Intensity或RMS测量模块也可以直接读出有效值变化趋势。在Simulink里右键信号线,选择“Signal Statistics”就能看到信号的峰值、有效值等统计量,做报告时非常好用。

5. 建模踩坑实录:常见问题与排查技巧

5.1 仿真报错“Numeric Singularity”怎么办

这个报错几乎每个做电力系统仿真的人都遇到过。原因一般是电路里出现了闭环的纯电感回路、电压源直接并联、或者变压器绕组连接方式导致寄生回路。小电流接地系统仿真的重点排查对象是消弧线圈支路:如果你在中性点和地之间直接串入理想电感,而系统中性点没有其他通路,就可能形成数值奇异。解决办法是在电感支路串入一个小电阻(比如0.1欧姆),或者在电感两端并联一个很大的电阻(比如1MΩ),给数值求解提供一个稳定的参考路径。

变压器中性点你没有正确引出也会出问题。检查变压器的绕组连接是否选择了Yg或YN,如果选成了Y(不引出中性点),那么你后面接消弧线圈的端子根本不存在,模型自然就报错了。

5.2 故障后非故障相电压升高不到线电压

有朋友问我,为什么B、C相电压升高后只有8kV而不是理论上的10.5kV?这个问题的根源多半在于系统阻抗太大,故障电流流过电源内阻和线路阻抗产生压降,导致母线电压被拉低。解决办法是检查电源模块的内部阻抗设置,以及变压器短路阻抗百分比,这两个值决定了系统在故障状态下的电压支撑能力。把电源内阻调小、短路阻抗调低,电压升高幅度会趋于理论值。另外,故障模块的接地电阻如果太大,也会让故障相电压无法真正降到接近零。

5.3 消弧线圈投入后零序电流不减少

消弧线圈投入后零序电流没有明显变化,通常是消弧线圈的电感值计算错误,没和系统对地电容形成并联谐振条件。用我前面给的公式重新算一遍,注意每相对地电容C0必须是线路电容值,不是总电容。如果你的是多条线路并联,要把所有线路的对地电容加起来再算。还有一个容易忽略的地方:零序电流测量点要选在消弧线圈支路和线路支路的公共节点上,才能看到补偿后的总电流,选错测量位置会得出完全错误的结论。

5.4 波形出现严重的数值振荡

故障瞬间的高频振荡如果是那种密密麻麻的毛刺,不是真实的物理过程,而是数值振荡。解决办法比较直接:减小仿真步长,或者在测量信号后面加一个滤波环节,或者对故障模块设置较小的过渡电阻。有时候把线路模型从π型改为分布参数模型也能改善高频振荡,但计算量会增大不少。对一般研究目的,π型线路已经把模型的复杂度和精度平衡得很好了。

常见问题核心原因快速处理方案
Numeric Singularity报错电感回路无阻尼通路串联小电阻或并联大电阻
非故障相电压偏低系统阻抗过大调小电源内阻/变压器短路阻抗
零序电流不减小消弧线圈参数不对重新核算对地电容与谐振电感
波形毛刺严重数值振荡减小步长、加滤波器或调小接地电阻
电压瞬间被拉低无法恢复故障电阻设置过小或故障未切除检查故障时序和过渡电阻值

6. 让模型发挥更大价值:扩展方向与最后心得

6.1 把模型接到选线算法里做闭环验证

很多做继电保护的人拿到模型不只是为了看波形,而是要验证自己写的选线算法。我的建议是:在仿真模型里加一个To Workspace模块,把故障期间的零序电压、零序电流数据存到MATLAB工作空间,再通过MATLAB脚本读取并运行你写的算法。这样仿真和算法之间就能形成闭环,选线算法的正确率和抗干扰能力就能批量测试。

现在许多研究都在探索暂态量选线、行波选线、深度学习选线等思路,而这些方法都需要海量的故障波形样本作为支撑。只要把仿真模型里的故障类型、故障角、过渡电阻、线路参数改成变量,然后用脚本批量循环修改参数自动运行仿真,你就能积累足够多的样本数据,这个价值远远超过看几张波形图。

6.2 从单相故障扩展到其他故障类型

这个模型只演示了最常见的单相接地故障,实际上你只需要在Three-Phase Fault模块里切换故障类型,就能仿真两相短路、两相接地短路、三相短路等不同类型的故障。不同故障类型下,零序电压电流的特征完全不同,通过统一对比,你能对小电流接地系统的各种故障形态建立起直觉。我强烈建议把五种典型故障都跑一遍,并把结果整理成表格存档,后面做故障诊断研究和写论文时都是现成的素材。

6.3 我自己反复调整后发现的最有价值的一件事

最后说一个我自己的经验:建模初期不要急着做参数优化的花活,先老老实实把仿真结果和理论值对一遍。比如中性点不接地系统发生单相金属接地时,非故障相电压就应该是线电压,3U0就应该是相电压的三倍,这些数值对不对一验便知。只有基波特征都验证过了,再去看暂态谐波、弧光接地这些进阶内容,否则你在错误模型的基础上做再多分析都是空中楼阁。这套模型我后来改造成标准化的教学用例,每年给学生上手练一遍,凡是肯动手搭模型的人,对继电保护和小电流接地系统的理解深度,都远超过只看课本的人。这个东西值得你花几天时间认真吃透。

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

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

立即咨询