配电网单相接地选线仿真：基于MATLAB/Simulink的建模与验证

干配电网继电保护的同仁，应该都有过这种体验：母线发单相接地报警，选线装置给出一个“疑似故障线路”，结果现场试拉了一轮，真正故障的却是另一条线。尤其是在中性点经消弧线圈接地的小电流系统里，故障残流小、相位关系飘忽不定，选线这个老问题硬是坑了一批又一批人。我的做法是，与其靠经验去猜，不如直接用MATLAB/Simulink把不同中性点接地方式下的单相接地故障仿真模型搭出来，用波形图和零序分量数据说话。这套模型不仅能复现“故障线路零序电流最大”“方向相反”这些理论结论，还能帮我们验证选线判据在暂态和稳态情形下的实际表现，算得上选线方案开发前最划算的一笔投入。本文就围绕这套模型展开，适合正在做配电网故障分析、选线装置算法验证或电力系统课程设计的读者参考。

1. 单相接地选线这么难，仿真模型能帮上什么忙

1.1 选线难题的根源：故障量实在太小

小电流接地系统（中性点不接地、经消弧线圈接地）发生单相接地后，故障点流过的仅是系统对地电容电流，最小的几安培，大的也就几十安培，比负荷电流小一两个数量级。而且三相线电压仍保持对称，负荷照常运行，继电保护不像相间短路那样立即跳闸。规程允许带故障运行一段时间，可值班员心里都清楚，非故障相电压已经升高到线电压，绝缘薄弱点随时可能发展成两相短路，电弧还会引起过电压。选线的目的，就是在众多出线里找到发生接地的那一条，然后安排转移负荷或者停电处理。

难点在哪儿呢？稳态情况下，零序电流的幅值和相位在消弧线圈补偿后变得很不明显。特别是过补偿运行时，故障线路的零序电流方向跟不接地系统相反，如果选线装置还沿用稳态幅值或方向判据，就容易选错。暂态过程虽然信息丰富，但衰减快、频率高，传统电磁式PT/CT难以准确传变，也给工程实现增加了不少难度。说到底，选线本质上是一个“信号弱、干扰大、判据适应性差”的问题，需要把故障机理吃透才能解决好。

1.2 仿真模型能提供什么价值

搭建仿真模型，不是简单去验证课本结论，而是给自己提供一个可以反复折腾的“试验台”。比如说我想研究消弧线圈补偿度从欠补偿到过补偿，零序电流相位到底怎么翻转；想看看故障发生时刻在电压过零和电压峰值时，暂态分量有多大区别；想验证五次谐波判据在长线路占比很高的电网上还灵不灵——这些在真实电网里动一次刀很贵，但在Simulink里改一个参数、点一下运行，波形和录波数据立刻就能出来。

另一个关键价值在于算法验证。选线算法的开发往往需要大量样本数据来整定阈值和评估鲁棒性，单靠现场录波肯定不够。仿真模型配合Matlab脚本，可以批量修改故障线路、故障位置、故障电阻、故障时刻、系统运行方式等参数，一次跑出几百个案例，导出的零序电流零序电压数据直接喂给选线判据做统计分析。这种“仿真+算法验证”的工作流，是目前我做过的最省力的选线研究方式。

2. 中性点三种接地方式在Simulink里的建模区别

2.1 中性点不接地：最容易被忽略的对地电容

建模前先把系统参数定下来。以最常见的10kV配电网为例：三相电源10kV/50Hz，变压器容量选6.3MVA或10MVA，高压侧Y接线，中性点引出。不接地就是变压器中性点悬空，不做任何连接，Simulink里对应的是三线制或者带中性点但不连接。

不接地系统最容易忽略的是线路对地电容。单相接地电流完全由线路和负载侧对地电容形成，所以线路模型里的零序电容参数必须给准。我一般用集中参数π型线路，每条出线取典型单位参数：正序电阻0.27Ω/km，正序电感1.33mH/km，正序电容0.339μF/km；零序电阻2.7Ω/km，零序电感4.13mH/km，零序电容0.28μF/km（架空线和电缆差别很大，具体以实际线路手册为准）。线路长度按5km、8km、12km、20km设置几回线，这样各线路的对地电容电流大小就有明显区分度，后面波形对比才好看。

2.2 经消弧线圈接地：补偿度是建模的灵魂

消弧线圈接在变压器中性点和地之间，模型里就是一个电感串联一个阻尼电阻，通常还会带一个小电阻模拟零序阻尼。电感值要根据系统对地电容电流来算：需要补偿的容性电流I_C约等于所有线路对地电容电流之和，消弧线圈在工频下的感抗X_L=U_ph/I_C。如果按过补偿10%设计，实际补偿电流I_L=1.1×I_C，感抗X_L=U_ph/(1.1×I_C)，电感L=X_L/(2πf)。

我举个例子：假设系统对地电容电流总和为30A，10kV系统相电压约5774V，过补偿10%后补偿电流33A，X_L≈175Ω，L≈0.557H。阻尼电阻通常取感抗的10%~15%，约17.5Ω到26Ω。仿真时把这个RLC串联支路放到中性点与地之间即可，Simulink里用Series RLC Branch或者Three-Phase Series RLC Branch，注意电感后必须可靠接地。如果要研究五次谐波选线，最好给变压器或电感做饱和特性，不然谐波源缺失，后面验证判据会落空。

2.3 经电阻接地：选线判据反而最简单

中性点经电阻接地分为小电阻接地（几十欧以下）和中电阻接地（100~1000Ω）。小电阻接地时，单相接地故障电流可达数百安培，故障特征非常明显，零序电流保护可以直接动作跳闸，选线压力不大。这个时候中性点与地之间只接一个三相电阻支路，电阻值按工程实际取6Ω、10Ω或20Ω都行。中电阻接地则介于两者之间，主要用于限制过电压，同时提供一定的零序电流。

我把三种接地方式搭成同一个系统的变子部分，也就是用三个不同的中性点支路模块，通过手动开关或者变子系统（Variant Subsystem）切换。这样做的好处是：线路和故障参数完全一致，只改中性点支路，就能在同一个Simulink环境下对比三种接地方式的零序特征，不需要复制三个大模型，后面做参数扫描也清爽得多。

3. 从零开始搭一套选线仿真模型：关键模块与参数

3.1 系统级拓扑：电源、变压器、线路、负载

我习惯用Simscape Electrical里的Specialized Power Systems库（老版本叫SimPowerSystems）。搭建顺序：三相电源→变压器→母线→多条π型线路→负载；故障模块放在需要考察的出线上；测量模块放在母线出线断路器处。

电源参数比较容易踩坑。Three-Phase Source里填的是线电压有效值，比如10kV系统填10e3，频率50Hz，接线方式选Y接，中性点不接地。变压器用Three-Phase Transformer（Two Windings），高压侧Yg，低压侧Δ；高压侧中性点引出后接入中性点支路。如果变压器模型选Yg，就会自带一个中性点端子，方便接消弧线圈或电阻。线路我用Three-Phase Pi Section Line，每段按5km来串联，两段就是10km，这样既能近似长线的分布效应，又不至于让仿真慢到没法忍。线路末端接三相并联RL负载，参数按每条线路所带负荷设置，比如5MW+j0.5Mvar。

3.2 故障与测量模块的放置

故障用Three-Phase Fault模块，选择A相接地，在Phase-A Fault那栏打勾。故障时刻我习惯用外部控制信号，接一个Step模块：0.2s之前为0，0.2s时跳变到1，触发故障。故障电阻设成0.1Ω（金属性接地），也可以做成变量，在Matlab脚本里循环赋值，专门跑高阻接地场景。

测量方面，母线处接Three-Phase V-I Measurement，把三相电压电流引出来。为了得到零序电压和零序电流，最干净的做法是在Matlab脚本里合成：i0=(ia+ib+ic)/3，u0=(ua+ub+uc)/3。注意三相信号必须时间对齐，离散仿真里如果不同数据点没有对齐会有误差，所以我把To Workspace的采样时间设成跟模型步长一致，并保存为Structure With Time格式，后面处理比较方便。Simulink里也有零序分量提取模块（Sequence Analyzer），但我觉得脚本合成更灵活、不占模型资源。

3.3 求解器与仿真参数的设置

这是很多模型跑不出理想波形的关键。我推荐用离散仿真，Powergui的仿真类型选Discrete，采样时间设成10μs或20μs，这样才能捕捉到几百赫兹到几千赫兹的暂态高频分量；如果只关心稳态基波波形，50μs也够。求解器选ode23tb或ode15s，容差保持默认。仿真时长设置0.4s左右：前0.2s建立稳态，0.2s发生故障，后面0.2s用于观察暂态和稳态。如果要做暂态能量判据，建议把波形记录窗口覆盖故障前至少一个周波，方便对齐零序电压过零点。

核心参数参考如下表：

4. 故障波形里藏着哪些选线特征（重点看零序）

4.1 不接地系统的波形：幅值最大+方向相反

以4条出线为例，线路1长20km（电容电流最大），线路2长12km，线路3长8km，线路4长5km。假设线路4在0.2s发生A相金属性接地。母线零序电压从0跃变到接近相电压峰值（约5774V的幅值）。观察各线路首端零序电流：非故障线路1/2/3的零序电流是各自线路的对地电容电流，幅值约几安到十几安，相位超前零序电压90°（容性方向）；故障线路4的零序电流等于其余三条线路零序电流之和，幅值明显最大，相位滞后零序电压90°（从线路流向母线）。示波器上看起来，故障线路4的零序电流和线路1的零序电流波形几乎反相，这就是稳态幅值加方向选线的基础。

这里要提醒一点：如果故障线路恰好是短线路，而某条正常线路很长，故障线路电容电流相对非故障线路优势就不大。比如线路4只有5km，线路1有20km，故障线路零序电流幅值理论上是线路1+2+3之和，仍然大于任意单条非故障线路，所以幅值比较在理论上成立，但现场CT误差大，实际还是倾向用方向判据或者暂态判据。

4.2 经消弧线圈接地后的相位翻转现象

同样的故障场景，把中性点换成消弧线圈，过补偿10%后，故障线路零序电流变成补偿后的残流。因为过补偿时电感电流大于系统总电容电流，故障线路的零序电流方向变成与不接地系统相反，超前零序电压90°（感性方向）。从波形上看，故障线路4的零序电流相位与正常线路1-3的相位相同，都是超前，如果还套用“方向相反”的判据，必然误判。幅值上，故障线残流可能很小，有时候甚至小于某条长正常线路的电容电流，稳态幅值判据也不可靠。

如果改成欠补偿，故障线路零序电流方向又会变回容性方向，但系统运行方式一变，补偿度就会漂移，稳态特征非常不稳。这时候必须看暂态：故障开始的最初几个毫秒，零序电压突变量和零序电流突变量里包含丰富的暂态分量，其中故障线路的暂态零序电流首半波极性与零序电压相反，非故障线路相同。在Simulink波形里能看到故障线路的高频毛刺明显更大，尤其在相电压峰值附近接地时。

4.3 经小电阻接地时波形清晰得令人感动

中性点经10Ω电阻接地，单相接地故障电流主要由中性点电阻提供的零序电流构成，数值可达几百安。各线路零序电流中，故障线路的零序电流幅值显著大于非故障线路，方向相反，而且基波分量非常稳定。波形上几乎不用滤波，直接用基波有效值幅值比较就能选出故障线。这也是为什么很多城市配电网改成小电阻接地——代价是接地故障必须立即跳闸，供电可靠性确实降低了，但保护的选择性和灵敏性都大幅改善。

三种接地方式波形特征对比如下表：

5. 选线判据在Simulink/Matlab里的实际验证手法

5.1 零序电流幅值/方向判据的脚本化实现

模型跑完后，我把To Workspace导出的三相电流数据放到Matlab脚本里合成零序电流，然后计算基波幅值和相位。幅值判据就是比较各线路零序电流的基波有效值；方向判据需要结合零序电压相位，看零序电流是超前还是滞后。代码框架大致这样：

% 假设 ts 为采样时间，ia1,ib1,ic1 为线路1的三相电流信号 i0_1 = (ia1.data + ib1.data + ic1.data) / 3; Fs = 1 / ts; N = length(i0_1); Y = fft(i0_1); % 基波对应的频率分量索引（忽略直流） f = (0:N-1) * Fs / N; idx50 = find(abs(f - 50) < 1); amp1 = abs(Y(idx50)) / N * 2; % 基波幅值 phase1 = angle(Y(idx50)); % 基波相位

不过实际我更推荐用Powergui的FFT Analysis工具，它可以直接看选定信号的基波幅值相位，不用自己在脚本里算。批量验证判据时，再写循环脚本批量提取数据，统一计算。核心是数据要干净，负序、零序分量没有混叠，所以采样时间务必要小。

5.2 五次谐波判据的仿真条件

五次谐波判据针对消弧线圈系统，原理是消弧线圈只在基波下补偿，对五次谐波感抗变成基波的5倍，所以五次谐波电流基本不补偿，仍按容性分布；故障线路的五次谐波电流方向与非故障线路相反。要在仿真中造出五次谐波，需要变压器饱和或者在线路中注入谐波源。我用的是让变压器模型勾选Enable saturation，设置饱和特性参数，让励磁电流在过电压时产生畸变，自然带出5次谐波。

我实测发现一个问题：只靠变压器饱和产生的5次谐波非常弱，只有基波百分之几，加上测量噪声后不太容易分析。更可控的做法是用一个并联的谐波电流源（Simulink里用Three-Phase Harmonic Suppressor或者直接加电流源），在母线上注入1-2A的5次谐波零序电流，模拟实际系统里的背景谐波。这样各线路5次谐波零序电流分布就很清晰，故障线路方向相反的特征也明显。

5.3 暂态首半波与暂态能量判据

首半波原理：接地瞬间的暂态电容电流远大于稳态电容电流，且频率高，故障线路与零序电压的极性关系在首半波时区别于正常线路。Simulink里要看到准确的暂态过程，采样时间必须小于20μs。从波形上手动找首半波极性不难，但批量验证就要写检测脚本：先找零序电压故障后第一个过零点，再看零序电流在这个区间里是正还是负。

暂态能量法更抗干扰，定义暂态零序能量为W = ∫v0(t)·i0(t)dt，在故障后一个或半个周波内积分。正常线路W为正或很小，故障线路W显著偏大且符号相反。我通常用梯形积分函数trapz实现，然后对大量故障案例统计W的分布，整定选线阈值。仿真模型的优势在这里体现得很充分：故障位置、故障电阻、故障起始角全都能扫，样本量比现场录波不知道多到哪里去了。

6. 我在这套模型上踩过的坑与调参经验

6.1 线路测量极性搞反，选线方向全反

头一回搭模型时，我把线路零序电流互感器极性接反，导致正常线路和故障线路的方向特征完全颠倒，选线判据跟着错位。后来逐个检查测量模块的正方向定义：Three-Phase V-I Measurement在母线出线侧测线路电流，正方向默认是母线流向线路；故障点设置在线路上，短路电流从两侧流向故障点。实际工程中零序CT的极性也有统一要求，仿真阶段如果不把方向理清楚，后面所有分析都白搭。

检查方法很简单：不接地系统单相接地时，非故障线路零序电流是容性的，超前零序电压90°；如果画出来是滞后，大概率极性反了。另外可以对比两条不同长度非故障线路的零序电流幅值，应该大致跟线路长度（电容）成正比，这一点也能帮助确认接线正确。

6.2 消弧线圈电感整定不能只算线路电容

我一开始计算消弧线圈补偿电流时，只把各条线路线路电容电流加起来，忽略了配电变压器低压侧对地电容以及母线侧杂散电容，结果补偿度偏小，故障残流巨大，波形完全失真。后来用Powergui的Impedance vs Frequency功能扫了母线零序阻抗，找到工频附近的谐振点，从谐振点反推出系统实际总电容电流，再回头调整电感参数，模型才靠谱。

具体做法是：在母线处加一个很小的扰动电流源，频率从1Hz到1000Hz扫频，观察母线零序电压的幅值响应。谐振峰对应的频率和系统电容、电感参数直接关联。这个功能不复杂，但非常实用，建议搭完模型先跑一次扫频再设置消弧线圈参数。

6.3 故障时刻和故障电阻对暂态波形的影响极大

故障时刻选在相电压峰值与过零，暂态分量差出一个数量级。峰值附近接地时du/dt最大，暂态电容电流可达数十安培；过零附近接地时暂态很小甚至看不清。故障电阻提升到1000Ω时，暂态分量被明显阻尼，首半波极性几乎淹没在噪声里。很多论文里给出的漂亮波形都是金属性接地加电压峰值时刻的极端工况，实做时不能只看这种理想结果。

所以我做验证时会用parfor批量扫描：故障时刻从0.18s到0.22s步长0.002s，故障电阻从0.1Ω到1000Ω按对数均匀分布，每条线路都轮流做故障线。跑完一批数据，统计各判据的准确率和失效工况，再针对失效区域决定要不要加辅助判据。这套流程我已经重复了很多次，是选线研究中性价比最高的一步。

6.4 数据保存与内存优化的两个小技巧

To Workspace模块如果同时保存十几路三相波形，采样时间10μs，仿真0.4s，就是十几条×3×40000点，内存压力不小。我的习惯是只保存需要分析的三相电流和母线三相电压，其他中间量用Terminator或者干脆不接入记录模块。其次，To Workspace输出格式建议选Structure With Time，这样在脚本里直接取.data和.time字段，不用再处理bus数据结构。

另外，离散仿真数据默认是均匀采样，但要注意故障时刻0.2s和采样步长10μs不一定正好对齐，导致暂态起始点有小偏差。影响不大，但如果做首半波算法验证，建议把故障触发信号也记录到WorkSpace，确保取数据窗口时能精确对齐故障瞬间。

最后分享一个让我效率翻倍的办法：把三种中性点接地支路封装成Variant Subsystem，用SS中性点选择变量控制到底接入消弧线圈、电阻还是悬空。我把这个变量设成Simulink参数，配合脚本循环，一套模型就能把三种接地方式、几十条线路参数、上百组故障工况全部跑完。选线仿真这件事，模型本身不是核心，核心是你能不能快速批量制造出足够多样的故障样本数据来检验你的判据。这套模型的价值就在于，把原本需要现场蹲守才能积累的故障录波，变成了按一下运行就能得到的东西。