搞配网故障分析的人十有八九都被小电流接地系统的建模恶心过。单相接地故障选线这个东西本身逻辑不复杂真正麻烦的是仿真模型要覆盖多种运行方式——今天中性点不接地明天经消弧线圈接地后天又改成经电阻接地每一种方式下的故障特征、零序分量分布、选线判据完全不是一个套路。你手动改一次模型结构还能忍改十次就真的想摔键盘了。我最近在 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 吸收回路或降低最大步长零序电压长时间不恢复消弧线圈品质因数太高、阻尼不足给电感串联合适的等效电阻Q20~30中性点悬空时电压悬浮不接地方式的断路器和系统中存在极小的对地泄漏路径在不接地支路上串一个 1MΩ 级高阻近似开路但避免数值奇异变压器励磁涌流影响故障波形变压器模块饱和参数设置不当稳态验证阶段先关闭饱和故障暂态分析再开启零序分量幅值随仿真时间漂移变步长求解器在稳态下步长过大限制最大步长不超过 1e-4s或改用固定步长仿真Sequence Analyzer 输出滞后序分量分析需要数据窗口时间波形稳定后再记录至少等待两个工频周期另外还有一个特别坑的地方是 Three-Phase Fault 模块设置成 A 相接地后其他两相的对地电容依然会通过电源内阻抗形成回路这在小电流建模里是正常的但如果你变压器和线路的零序阻抗没有设置准确得到的三相不平衡电压相位就可能跟你预期的不一致。遇到这种问题检查零序回路比检查故障模块参数更有用。模型验证完毕之后我又做了一次批处理测试写了一个循环把消弧线圈的补偿度从 5% 扫描到 12%每次仿真完自动从工作区提取故障线路的特征量并绘制曲线。这个过程中互锁逻辑发挥了很大作用因为批处理脚本不会像人一样思考很多参数组合会产生冲突互锁开关保证了三组断路器不会同时处于合闸状态批处理才能安全跑完。6. 关于这套模型的一些体会搭完这套模型之后我最大的感受是小电流接地系统建模的难点从来不在某个元件的参数有多难设而在于系统方案的灵活性和可复用性。当你把三种接地方式做成可切换的结构相当于把研究对象从一条具体的配网线路扩展成了面向整个中压配电网故障特征的通用验证平台。在实际操作中我建议大家搭模型的时候时刻记住这条原则仿真结构要尽量接近实际系统的操作方式。在地面系统里中性点设备的投切靠的是断路器在仿真里就也应该用断路器来模拟。这样不仅模型行为更真实而且当你需要扩展第四种接地方式的时候只需要再并一条支路、增加一组断路器就行架构上不需要做任何推翻重来的工作。后续我还打算在这套模型上做两件事一是把电缆和架空线的混合出线结构加上去让线路模型更贴近真实配网二是接入一个实际的选线装置嵌入式代码做闭环仿真测试。等有结果了再来分享。