搞配网故障分析的人,十有八九都被小电流接地系统的建模恶心过。单相接地故障选线这个东西,本身逻辑不复杂,真正麻烦的是仿真模型要覆盖多种运行方式——今天中性点不接地,明天经消弧线圈接地,后天又改成经电阻接地,每一种方式下的故障特征、零序分量分布、选线判据完全不是一个套路。你手动改一次模型结构还能忍,改十次就真的想摔键盘了。
我最近在 MATLAB 2020a 上把这套模型重新搭了一遍,把三种中性点接地方式做成了运行时可切换的结构,算是把这个老大难问题彻底解决了。模型不是说多高级,但胜在架构清爽:一个主体系统,三套中性点支路,通过开关组自由切换,变电站出线的对地电容、消弧线圈补偿度、接地电阻阻值全部参数化。这篇文章就把我的建模思路、关键参数选取和踩过的坑完整记录下来,给同样在做小电流选线仿真的朋友一个能直接抄作业的参考。
1. 小电流系统为何让建模人员头疼
1.1 一个故障场景引出的建模需求
先还原一下实际场景。配电网10kV线路发生单相接地故障时,系统允许带故障运行一到两个小时,但要尽快找到故障线路并隔离,否则容易扩大成相间短路。故障选线装置的研发和验证,靠现场做实验成本太高、风险大,只能靠仿真。可问题是,不同中性点接地方式下,故障电流的幅值和路径差异极大,选线算法要适配所有场景,建模就必须把三种方式都覆盖。
一开始我的做法很粗暴:搭三套独立模型。一套不接地,一套带消弧线圈,一套带电阻,需要验证哪种方式就打开哪个模型。结果就是项目代码量翻了三倍,每次修改对地电容参数要在三个模型里挨个改,稍不注意还会出现三套模型参数不一致的问题,仿真结果互相矛盾时根本没法排查。
调参调了半个月之后我终于想明白:选线算法面对的应该是同一套系统,只是中性点处理方式不同。建模的逻辑应该跟实际操作保持一致——用一个断路器去投切中性点设备,而不是从一个模型复制粘贴到另一个模型。
1.2 三种中性点接地方式的本质差异
从电气原理上讲,这三种接地方式本质上是调控故障电流路径和幅值的手段。
中性点不接地(IT系统)是最原始的方式。单相接地时故障电流只能通过对地电容形成回路,电流不大,一般10kV线路电容电流在几安到几十安之间,电弧往往能自熄灭。但非故障相的电压会升高到线电压水平,长期运行对绝缘是考验。
经消弧线圈接地(谐振接地系统)是在中性点和地之间串入一个电感量很大的线圈。单相接地时,电感电流与电容电流方向相反,互相抵消,故障点的残余电流可以被压缩到5A以内,电弧基本瞬间熄灭,这是目前中压配网用得最多的方式。代价是接地故障时故障线路的特征量被严重削弱,选线难度反而增加了。
经电阻接地(通常是中电阻或低电阻)的思路是反过来的:人为增大故障电流到几十安甚至几百安,通过零序过流保护和方向判别就能快速、精准地切除故障线路。代价是故障点能量大,系统不能带故障运行。
三种方式的故障特征差别,我整理了一个对照表方便大家理解:
| 接地方式 | 单相接地故障电流 | 非故障相电压 | 故障线路零序特征 | 选线难度 |
|---|---|---|---|---|
| 不接地 | 数A~数十A(电容电流) | 升高√3倍 | 零序电流最大,方向与非故障线路相反 | 中等 |
| 经消弧线圈 | 残流小于5A | 升高√3倍 | 零序电流被补偿,方向可能反向 | 高 |
| 经电阻 | 数十A~数百A | 被限制 | 零序电流幅值大、方向明确 | 低 |
看到这个表就明白了,建模的核心任务不只是把三条支路画出来,而是要让模型在切换后自动呈现出各自的特征差异,选线算法才能有的放矢。
2. 模型架构设计:把三种模式装进同一个模型
2.1 从“改参数”到“切开关”的方案演进
先说结论:在 Simulink 里实现三种接地方式切换,最好用的不是什么高级模块,而是一组三相断路器。把消弧线圈和接地电阻做成两条可投切的支路,再配合一个直接接地支路,通过断路器的开合状态组合来选择当前的接地模式。
第一次尝试的方案是在模型里面放一个手动开关,直接把中性点引线切到不同支路上去。听起来直接,但 Simulink 的手动开关在仿真过程中无法动态切换,而且切换瞬间容易产生数值突变,步长直接崩掉。后来换成了可控断路器(Three-Phase Breaker),用 Step 信号或者手动置位来控制开合,问题就解决了。2020a 的 Simscape Electrical 库里面断路器模块支持外部控制信号,可以在仿真运行中实时投切,这样就能做动态场景测试——比如从谐振接地切换到电阻接地,观察选线装置在过渡过程中的响应。
整体架构上,我把它分成了三个层级。顶层是电网电源与主变压器,中间是出线线路及负荷,底层是中性点设备组合。每一层都是独立的子系统,层与层之间用三相电气端口连接,方便之后替换参数或增删线路。
2.2 主变压器与中性点引出的关键细节
一个小细节特别容易坑人:10kV 配电网的主变压器通常是 Y/Δ 接线,低压侧(10kV 侧)是 Δ 接法,天然没有中性点引出。如果要搭建经消弧线圈或经电阻接地系统,必须先加一台接地变压器(通常采用 Z 形接线的接地变)人为制造中性点。
仿真模型里我用的是 Simscape Electrical 的 Three-Phase Transformer (Zigzag) 模块,如果组件库里没有现成的 Z 形接法,也可以用三个单相变压器按 Z 形绕组方式自己搭,但绕组极性一定要核对清楚,接反了零序阻抗就完全不对了。
中性点引出后,再接一个三绕组测量模块,方便零序电压的获取。主变压器的连接组别我设置的是 Dyn11,这是国内配网最常见的接线方式,模拟结果更贴近工程实际。
2.3 切换逻辑与互锁设计
三种接地方式的切换用三组断路器来实现。我管它们叫 BRK_Direct(直接接地)、BRK_Arc(消弧线圈投入)、BRK_Res(电阻接地)。逻辑说起来很简单:
- 不接地:三组断路器全部分闸,中性点对地悬空
- 经消弧线圈:BRK_Arc 合闸,其余分闸
- 经电阻:BRK_Res 合闸,其余分闸
为了防止人为误操作把两组断路器同时合上(比如消弧线圈和电阻同时挂在中性点上),我在控制逻辑后面加了一个互锁判断。实现方式不复杂,用几个 Logical Operator 模块做条件判断,保证任一时刻只有一组断路器能够合闸。这个保护机制在做自动化批量仿真时尤其重要——批处理跑几百组工况,一旦哪组参数把两个开关同时置1,仿真结果直接乱套,浪费时间不说,还容易得出错误结论。
运行时的切换我用了两种控制源:一种是 Step 信号,可以设定在仿真到特定时间时动作;另一种是手动开关,方便交互式验证。实际测试下来,断路器在合闸瞬间会产生一个过渡过程,如果系统中消弧线圈或电容储存的能量较大,这个暂态过程会长一点,所以在切换后最好预留 0.1~0.2s 的稳定时间再开始分析波形,这块我在后面故障注入部分还会细说。
3. 核心参数与建模实操要点
3.1 线路对地电容与消弧线圈补偿度的计算
仿真参数不是随便填的,建模之前需要先做一轮手算。以一条典型的 10kV 电缆出线为例,单公里电缆对地电容大约在 0.2~0.5μF 之间,架空线小一些,约 0.01~0.05μF/km。假设系统有五条出线,总长度加起来约 30km,平均对地电容按 0.15μF/km 估算,总对地电容 C 大约 4.5μF。
在 10kV 系统(相电压 Uφ = 5.77kV)、50Hz 工况下,总的零序电容电流可以这样估算:
I_C = 3 × Uφ × ωC = 3 × 5770 × 314 × 4.5×10⁻⁶ ≈ 24.5A
这个数值是设计消弧线圈补偿容量的直接依据。按过补偿度 8% 设计,消弧线圈的感性电流需要做到 26.5A 左右。对应的电感值可以反推:
I_L = Uφ / (ωL) → L = Uφ / (ω × I_L) ≈ 5770 / (314 × 26.5) ≈ 0.694H
建模时消弧线圈模块我直接填了电感值,注意 Simscape 里面电感串了个等效电阻(品质因数按 30 考虑),这样零序回路才有一个真实的阻尼,波形不会出现单调衰减的假象。
3.2 中性点电阻阻值怎么选才合理
经电阻接地系统的电阻参数也有讲究。阻值越小,故障电流越大,选线越容易,但对系统和设备冲击也大;阻值太大,故障电流与电容电流混在一起,选线的灵敏度又不够。
工程上有一个经验做法:电阻电流 IR 取电容电流 IC 的 1.5~2 倍,这样故障点的有功分量明显大于容性分量,方向判别才好做。沿用上面 IC ≈ 24.5A 的例子,IR 取 40A,电阻值就是:
R = Uφ / IR = 5770 / 40 ≈ 144Ω
我模型中最终设置的是 100Ω 固定电阻,对应 IR ≈ 57.7A,故障特征非常明显,而且对 10kV 系统来说完全在设备和电缆承受范围内。
3.3 线路模块与测量环节的选择
线路建模我用了 Simscape Electrical 的 Pi Section Line(π 型等值电路),每条出线有三段 π 型级联,而不是直接把整条线路用一个 π 型代替。原因是单相接地时零序电压电流在整条线路上分布并不均匀,多段 π 型能稍微逼真一点,仿真速度损失也不大。输电线路参数我按正序与零序分别设定,这里特别提醒一下:不少新手图省事只设置了正序参数,零序参数把它留空,结果出来的零序电流完全不对。零序阻抗尤其是零序电容,决定了对地故障路径,必须单独设置。
测量环节用的是 Three-Phase V-I Measurement 模块,接在母线出口和每条线路的首端。要测零序分量的话,我习惯接一个 Sequence Analyzer(序分量分析模块),直接输出零序电压和零序电流的幅值与相位。2020a 版本里 Sequence Analyzer 依然是通过相量方式计算序分量的,注意它需要一定时间的数据窗口才能稳定输出,所以我通常等仿真跑到稳定后再取值,而不是一上来就看第一个周期的分析结果。
3.4 求解器设置与仿真速度优化
这一步是我花时间最多的地方。小电流系统里同时存在小时间常数的暂态(断路器动作)、大时间常数的工频稳态(50Hz)以及非线性元件(变压器饱和、消弧线圈电感),求解器设置不当会非常痛苦。
我最终采用的是变步长 ode23tb 求解器,相对容差设为 1e-3,最大步长为 1e-4s。这个组合在精度和速度之间找到了一个不错的平衡点。之前试过 ode45,仿真速度是上去了,但断路器切换瞬间高频振荡捕捉不到,选线算法跑出来的结果跟实际偏差很大。
还有一个提速的技巧:如果只关心第二种方式的稳态和选线动作过程,可以把变压器饱和特性先关掉(设置 Simply 饱和模型),励磁涌流和铁芯饱和非线性会影响仿真初始阶段的收敛,而对故障稳态特性的影响可以忽略。测试完成后需要精确模拟故障暂态的时候再把它打开。用这个方式,原来一组 2 秒仿真时长跑 15 分钟,优化后大约 4 分钟就能跑完。
4. 故障注入与选线算法验证
4.1 故障模块的设置与触发逻辑
故障注入我用的是 Three-Phase Fault 模块,放置在某一回出线线路的 30% 处(从母线侧算起)。故障类型设置为单相接地(Phase A 对地),故障电阻我给了两个典型值:金属性接地用 0.01Ω,弧光接地用 10Ω 左右模拟电弧电阻。不同接地电阻下选线装置的灵敏度表现会有不小差异,多测几个阻值很有必要。
故障触发时间设定为 t = 0.5s,持续 0.5s 后自动切除。这时我会提前把中性点接地方式的切换设置在 t = 0.2s 完成,这样故障发生时系统已经稳定在目标接地模式,不会把切换暂态和故障暂态叠加在一起,波形分析会更干净。
4.2 三种模式下的波形特征对照
模型跑出来后,在 Scope 里面观察母线的零序电压和故障线路与非故障线路的零序电流,特征非常清楚。
中性点不接地模式下,故障线路零序电流明显大于非故障线路,相位大约相差 180°,但不接地系统的容性电流大小跟线路总长度强相关,短线系统故障电流可能只有几安培,对测量精度要求很高。
经消弧线圈接地模式下,故障线路的零序电流变得很小,被消弧线圈的感性电流抵消了一部分,在过补偿状态下故障线路零序电流的相位会与非故障线路趋于一致,仅靠幅值和方向就没法选线了。这时候需要依赖五次谐波法或暂态量分析法,模型中我用到了 FFT 分析模块提取五次谐波分量。
经电阻接地模式下,故障线路零序电流和零序电压的关系就很清楚了,零序功率方向具有唯一性,直接用传统的零序功率方向保护就能选出来。
三种模式下故障点的零序电压波形都近似等于相电压幅值(约 5.77kV),这是单相接地故障的共性特征,可以用来作为选线启动判据。
4.3 选线判据在模型中的落地
我在这套模型上验证了三种典型的选线算法。
第一种是零序电流幅值法,适用于中性点不接地系统。仿真实测零序电流:故障线路约为 60A,非故障线路约为 12~20A,幅值差异很显著,取最大者即可选线。用 MATLAB 的max()和find()组合实现很直接。
第二种是五次谐波法,用于消弧线圈接地系统。因为消弧线圈只在基波频率下完全补偿,五次谐波频率下它的感抗是基波的五倍,补偿效果大幅削弱,所以故障线路的五次谐波零序电流依然大于非故障线路。模型的输出经过 FFT 后取 250Hz 分量,可靠区分。
第三种是零序功率方向法,用于电阻接地系统。在 Simulink 中可以用瞬时功率计算模块对零序电压和零序电流的瞬时值做乘积再取平均,功率为负值的即为故障线路。
三种算法在一个统一的选线判断脚本中运行,数据直接从仿真工作区导入。之前手动改模型的结构,现在切换参数后直接运行统一脚本,验证效率高出一大截。
5. 常见问题与排查记录
这个模型从零搭起来到基本稳定,踩的坑不算少。我把典型的几个问题整理成一个排查表,大家遇到类似问题能少走弯路。
| 现象 | 可能原因 | 排查思路及解决方式 |
|---|---|---|
| 断路器切换瞬间仿真步长崩掉 | 切换时系统储能元件状态突变 | 使用 Breaker 模块的 Snubber 参数(默认 RC 吸收回路),或降低最大步长 |
| 零序电压长时间不恢复 | 消弧线圈品质因数太高、阻尼不足 | 给电感串联合适的等效电阻(Q=20~30) |
| 中性点悬空时电压悬浮 | 不接地方式的断路器和系统中存在极小的对地泄漏路径 | 在不接地支路上串一个 1MΩ 级高阻,近似开路但避免数值奇异 |
| 变压器励磁涌流影响故障波形 | 变压器模块饱和参数设置不当 | 稳态验证阶段先关闭饱和,故障暂态分析再开启 |
| 零序分量幅值随仿真时间漂移 | 变步长求解器在稳态下步长过大 | 限制最大步长不超过 1e-4s,或改用固定步长仿真 |
| Sequence Analyzer 输出滞后 | 序分量分析需要数据窗口时间 | 波形稳定后再记录,至少等待两个工频周期 |
另外还有一个特别坑的地方,是 Three-Phase Fault 模块设置成 A 相接地后,其他两相的对地电容依然会通过电源内阻抗形成回路,这在小电流建模里是正常的,但如果你变压器和线路的零序阻抗没有设置准确,得到的三相不平衡电压相位就可能跟你预期的不一致。遇到这种问题检查零序回路,比检查故障模块参数更有用。
模型验证完毕之后我又做了一次批处理测试:写了一个循环,把消弧线圈的补偿度从 5% 扫描到 12%,每次仿真完自动从工作区提取故障线路的特征量并绘制曲线。这个过程中互锁逻辑发挥了很大作用,因为批处理脚本不会像人一样思考,很多参数组合会产生冲突,互锁开关保证了三组断路器不会同时处于合闸状态,批处理才能安全跑完。
6. 关于这套模型的一些体会
搭完这套模型之后,我最大的感受是:小电流接地系统建模的难点从来不在某个元件的参数有多难设,而在于系统方案的灵活性和可复用性。当你把三种接地方式做成可切换的结构,相当于把研究对象从一条具体的配网线路,扩展成了面向整个中压配电网故障特征的通用验证平台。
在实际操作中,我建议大家搭模型的时候时刻记住这条原则:仿真结构要尽量接近实际系统的操作方式。在地面系统里,中性点设备的投切靠的是断路器;在仿真里,就也应该用断路器来模拟。这样不仅模型行为更真实,而且当你需要扩展第四种接地方式的时候,只需要再并一条支路、增加一组断路器就行,架构上不需要做任何推翻重来的工作。
后续我还打算在这套模型上做两件事:一是把电缆和架空线的混合出线结构加上去,让线路模型更贴近真实配网;二是接入一个实际的选线装置嵌入式代码,做闭环仿真测试。等有结果了再来分享。