小电流接地系统单相接地故障选线：Simulink仿真与零序电流分析

小电流单相接地系统故障选线：MATLAB Simulink模型仿真与理论支撑下的研究

小电流接地系统的单相接地故障选线，可能是配电自动化方向里最磨人的课题之一。我接手第一个用MATLAB Simulink做故障选线仿真的任务时，以为搭几条出线、挂一个三相短路故障模块就能轻松复现现场录波，结果被零序电流的方向问题折腾了整整一周。后来才明白，这类仿真难的不是建模操作，而是模型能不能把工程实际中的故障物理过程忠实地承载出来。

这个方向的专业跨度其实不小：既要懂中性点接地方式与零序网络的电气原理，又要把暂态、稳态的故障特征量提取出来，最后还要在Simulink里把线缆参数、消弧线圈、故障时刻、采样率这些细节一件件对清楚。这篇文章把我从理论判据到模型搭建、从波形分析到参数整定的完整思路整理成一份可复现的记录，适合正在做配网故障选线仿真、写相关论文，或者准备用Simulink做继电保护算法验证的同学直接参考。

1. 为什么单相接地选线一直是个老大难

1.1 三种接地方式决定了故障电流的量级

小电流接地系统，业内通常指中性点不接地系统和中性点经消弧线圈接地系统。国内配电网大量采用这两种方式，核心原因是单相接地故障时可以带故障继续运行一段时间——中性点不接地系统允许运行不超过两小时，消弧线圈系统因为故障点残流被补偿得更小，电弧容易自熄，很多时候故障会自动消失。

单相接地时，接地相电压下降，非故障相电压升高到线电压，但三相之间的线电压依然对称，所以三相负荷还能勉强维持运行。代价是系统对地绝缘长期承受过电压，一旦另一个点再发生接地，就会发展成两点接地短路，甚至变成相间故障。因此运行规程要求在带故障运行的窗口期里尽快找到故障线路并处理。

如果中性点改为经小电阻接地，故障电流可以达到几百安培到上千安培，保护很容易检测，直接跳闸切除。但代价是任何一次瞬时性接地都会引起用户停电。对供电可靠性要求高的区域，大量用户不能接受这种“一接地就跳”的方式，所以小电流接地系统至今仍是配网主流，选线问题也就一直存在。

1.2 故障特征量为什么那么微弱

单相接地故障的短路电流并不来自电源侧，而是来自线路对地电容的充放电。中性点不接地系统里，非故障线路的零序电流就是本线路对地电容电流，通常只有几安到几十安；故障线路的零序电流等于所有非故障线路电容电流的矢量和，虽然相对大一些，但和负荷电流、短路电流相比依然是“小信号”。

更麻烦的是消弧线圈的介入。消弧线圈串在变压器中性点和地之间，接地故障时产生的电感电流会去补偿线路电容电流。完全补偿状态下，故障点残流很小，弧光可以熄灭，但故障线路的零序电流幅值也会被压到极低，甚至相位发生反转。传统选线判据大多建立在“故障线路零序电流最大且方向相反”这个基础假设上，消弧线圈一补偿，这个根基就没了。这也是为什么工程上把这类系统的选线称为“老大难”：故障量小、方向不定、工况一多变判据就失效。

1.3 选线问题在实际工程中到底有多痛

实际运维里，单相接地发生后变电站的接地选线装置会报警，但经常出现选线结果和真实故障线路不一致的情况。运行人员只能采用“拉路法”——一条线一条线地试拉，看接地信号是否消失。这个方法原理简单，但每拉一条线路就意味着一次停电，对用户的影响特别大，尤其是拉到了带重要负荷的线路，投诉和考核都会立刻找上门。

即便选线装置给出了结果，现场还经常要面对高阻接地、弧光接地、消弧线圈过补偿不足等复杂情况。所以无论是做保护装置研发、故障指示器算法，还是做配电自动化主站的选线策略，第一步都是先把各种故障场景下的信号特征用仿真工具验证清楚。Simulink的优势在于它能把电力系统模型、故障触发逻辑、信号采集和后处理算法放在同一个环境里跑，改一个参数就能对比一组工况，这对选线算法研究来说非常合适。

2. 理论支撑：选线判据与特征量提取

2.1 稳态判据：幅值比较、方向比较与五次谐波

稳态选线的核心是零序电流幅值和相位的比较。中性点不接地系统发生单相接地后，对于非故障线路，零序电流是本线路对地电容电流，方向从母线流向线路；对故障线路，零序电流是所有非故障线路电容电流之和，方向从线路流向母线。所以理论上存在两个明显特征：故障线路零序电流幅值最大，相位与健全线路相反。

在Simulink里验证这个方法，最直接的做法是提取各条线路首端的零序电流，用FFT取基波幅值和相位，记录各线路的幅值和相位关系。低阻接地、故障时刻合理的情况下，这个判据非常可靠。但它的前提是流过每条线路的零序电流主要由对地电容决定，如果系统里有较大的不对称负荷，或者线路参数差异过大，幅值和相位的区分度都会下降。

五次谐波法是为了绕开消弧线圈的基波补偿而提出的。消弧线圈按基波频率整定，对150Hz的五次谐波呈现很大阻抗，补偿作用很弱。因此五次谐波分量里，故障线路和非故障线路之间的幅值、方向关系仍然近似满足稳态判据。这个方法的软肋在于系统中五次谐波源不稳定，如果上游没有足够的谐波背景，故障和非故障线路的五次谐波都不明显，判据灵敏度和可靠性都受影响。

2.2 暂态判据：首半波法与小波分析

单相接地的瞬间，故障相电压发生突变，对地电容通过故障点快速放电，同时健全线路的对地电容被充电，这个过程会产生幅值比稳态大数倍甚至十几倍的暂态电流，而且持续时间很短。首半波法就是利用故障后第一个半波内暂态零序电流的极性和幅值来区分故障线路和健全线路：故障线路的暂态电流首半波极性，与健全线路明显相反。

暂态信号需要足够高的采样率和带宽，否则首半波被滤没了。仿真时要让求解器步长足够小，后处理时还要加带通或高通滤波，保留暂态主导频段而滤掉稳态工频分量。小波变换是另一种常用的暂态分析工具，它能把信号按频带展开，在特定尺度上观察故障起始时刻的奇异性和能量分布。MATLAB里可以用小波工具箱的cwt或wavedec函数，也可以用dwt2之类的基础函数自己搭特征提取流程。实际工程落地时，硬件采样率往往有限，小波分析的频带选择要和采样率匹配，这一点在仿真阶段就要提前考虑，否则现场复现时算法性能会大打折扣。

2.3 信号注入法与综合判据

还有一种思路是主动注入信号。接地故障发生后，通过电压互感器或者专用信号源向系统注入一个区别于工频的特定频率信号，用信号探测器在馈线出口和沿线查找。这种方法不受消弧线圈补偿和过渡电阻的影响，但需要额外硬件和注入装置，适合有条件的配电站做辅助定位。

工程上真正实用的选线装置，几乎没有只靠单一种判据的。常见做法是同时计算稳态幅值、稳态相位、五次谐波、暂态能量等多个特征量，每个判据输出一个候选结果，再通过加权投票或置信度融合给出最终选线结论。Simulink里做综合判据的优势是可以搭建多路特征提取模块并行输出，最后用MATLAB脚本做融合决策。这样既能对比各判据在不同故障场景下的表现，也能提前暴露“某一判据在特定工况下失效”的问题。

2.4 判据与Simulink实现手段的对应关系

3. Simulink模型搭建思路与实操细节

3.1 系统拓扑与元件选型

Simulink里做配电网故障仿真，需要用Simscape Electrical的Specialized Power Systems库来搭建三相电路。电源用Three-Phase Source，变压器用Three-Phase Transformer（Two Windings或Three Windings），线路用PI Section Line，负荷用Three-Phase Parallel RLC Load，故障用Three-Phase Fault。

搭模型前先想清楚系统电压等级和接地方式。以最常见的10kV中性点经消弧线圈接地系统为例：10kV母线由上级110/10kV变压器供电，变压器低压侧通常是Y接或Z接，中性点引出来接消弧线圈。馈线三条起步，太多线路模型会变得很重，三条已经足够展示健全线路、故障线路和拓扑差异。每条线路用两个或三个PI段串联，模拟分布电容的集中等效。电源内阻抗不能设置为0，否则开关动作瞬间会产生数值震荡，这一点后面还会专门说。

故障模块的设置也要注意。Three-Phase Fault模块支持指定故障相、接地电阻和动作时间，可以设置成单相接地。如果要做高阻接地分析，把接地电阻设为几百欧到几千欧。动作时刻建议设在系统稳定运行之后，比如0.05秒，故障持续时间按需设定，之后可以自动切除。

3.2 零序电压与零序电流的提取方法

零序分量的定义是三相瞬时值相加的三分之一，这在Simulink里可以用三个电压表或电流表接到求和模块来实现。比如零序电压U0等于(A相电压+B相电压+C相电压)/3，零序电流的3倍3I0等于(A相电流+B相电流+C相电流)，直接用Add模块把三相电流加起来就是3I0，不需要再除以3。实际工程里电流互感器开口三角测到的就是3I0，直接提这个量更方便后续和实测数据对拍。

更推荐的做法是：把测量到的三相电压、三相电流通过To Workspace模块导出到MATLAB工作区，再写脚本做Fortescue对称分量变换。这个做法的好处是灵活，仿真完成后可以随时调整滤波频率、FFT窗长、判据阈值，不需要重新跑一遍模型。尤其是在做多次工况对比时，把所有原始波形批量导出、用脚本统一处理，比在Simulink里挂一屏示波器高效得多。

提零序量时必须注意参考方向的一致性。电流测量的参考方向如果是“从母线流向线路”，那么所有线路都用同一个参考方向，否则后面做相位比较时会出现系统性反转，让你误以为算法有问题。这个坑我踩过：三条线路中有一条的电流测量模块方向插反了，波形上看起来故障线路和非故障线路的方向关系完全乱掉，排查了半天才发现是测量方向的问题。

3.3 消弧线圈与补偿度建模

消弧线圈本质上是一个带铁芯的电感，仿真模型里既可以用Series RLC Branch直接放一个电感，也可以用专门的可控电抗器模型。但纯电感在仿真里很容易引起持续振荡，实际消弧线圈是有功损耗的，所以模型里必须串联一个电阻，阻值按消弧线圈品质因数折算，一般取感抗的5%左右。

补偿度的计算要先知道系统对地电容电流。假设三条馈线对地电容总和CΣ已知，那么不接地系统单相接地时的电容电流约为 (I_C = 3\omega C_\Sigma U_\phi)。消弧线圈完全补偿时满足 (3\omega C_\Sigma = 1/(\omega L))，即 (L = 1/(3\omega^2 C_\Sigma))。实际运行采用过补偿，补偿度p取5%到10%，此时电感值要略小于完全补偿值。仿真中可以通过调整电感值对比过补偿、欠补偿、完全补偿三种工况下的选线效果。

我常用的做法是把补偿度写成一个可调参数放在MATLAB脚本里，用assignin把数值写进Simulink工作区，这样批量跑工况时只需要改脚本参数，模型不用动。建模的时候还要注意，消弧线圈接在中性点和地之间，如果变压器没有引出中性点，需要先加一台Z形接地变压器来构造中性点。很多初学者直接在三角形侧接消弧线圈，结果完全不对，这个细节非常关键。

3.4 采样设置、求解器选择与信号滤波

仿真时长一般取0.1到0.3秒，故障发生在0.02秒之后，这样既能观察故障前稳态，又能记录完整的暂态过程。求解器推荐变步长ode23tb，因为它对电力电子开关和故障这类刚性事件处理得比较好。如果为了后处理方便要固定采样率，那就用定步长，步长一般取1e-5秒到1e-4秒。对50Hz工频系统，一个周期20毫秒，定步长1e-4时每周期有200个点，做基波FFT够用；但做暂态首半波分析，一个周期至少要有1280点才比较稳，建议步长不超过2e-5秒。

Simulink里的Scope模块显示的波形肉眼直观，但波形数据如果要用于定量分析，要设置成“Save data to workspace”形式。To Workspace模块建议用“Structure With Time”格式输出，这样时间向量和信号值一一对应，在MATLAB里处理起来最方便。滤波环节可以用Analog Filter Design模块实现，也可以用脚本在数据导出后处理，例如对零序电压和零序电流做带通滤波提取五次谐波，或者用高通滤波器剥离工频保留暂态。

4. 仿真流程与关键参数整定

4.1 典型参数选取：一个10kV三馈线案例

为让仿真结果有参考意义，给出一个典型案例参数。系统是10kV中性点经消弧线圈接地，母线侧变压器额定容量按10MVA考虑，阻抗百分比大概10%，这些细节对故障电流的计算有一定影响，不宜随意省略。三条架空馈线长度分别取8km、12km和15km，单位长度参数可以参考典型架空线型号取值：正序电阻0.27Ω/km，正序电感1.3mH/km，正序电容0.0095μF/km，零序电容可以取0.006μF/km的量级。电缆线路参数则差异很大，对地电容比架空线大将近一个数量级，如果仿真的是电缆网络，需要单独查厂家型号参数。

三条线路的走向和长度不要设成完全一样，因为现实中线路参数总是有差异，完全对称的模型会让选线仿真过于“理想化”，掩盖掉幅值比较法在参数不平衡时的潜在问题。负荷方面每条线路接一个恒阻抗负荷，功率因数按0.85左右设置，这样故障前电流和电压的相位关系比较贴近实际。

4.2 故障工况的设计方法

故障工况设计要覆盖故障相、故障时刻、故障电阻、补偿度这几个维度。故障相可以分别做A相、B相、C相接地，验证算法对相别的对称性。故障时刻通过设定故障模块的动作瞬间来改变故障初相角，电压过零附近和电压峰值附近故障的暂态特征差别非常明显。故障电阻建议做四档：0.01Ω模拟金属性接地，20Ω模拟弧光接地，1000Ω模拟经过树枝或潮湿表面的高阻接地，10kΩ模拟绝缘严重劣化。

消弧线圈补偿度设置为过补偿5%、过补偿10%、欠补偿5%、完全补偿四种。每种工况下记录各线路零序电流的波形、零序电压波形和处理后的特征量。把这批工况全部导出保存后，可以绘制一张工况对比表，清晰看出哪些判据在哪些工况下失真。

4.3 Simulink仿真结果的导出与自动选线实现

仿真跑完后，波形不会自己变成选线结论，需要写脚本处理。下面是一个典型的后处理流程，按采样率和FFT窗长提取各线路零序电流基波幅值和相位。

% 假设Simulink导出的结构体数组 data_i0 中保存了各线路3I0 Fs = 10000; % 采样率 f0 = 50; % 基波频率 T = 1/Fs; winLen = 0.04; % 两个工频周期窗 nWin = round(winLen * Fs); for k = 1:3 x = data_i0{k}.signals.values; t = data_i0{k}.time; % 取故障后稳定区间 idx = find(t > 0.1, 1); seg = x(idx:idx+nWin-1); % FFT取基波 X = fft(seg .* hann(nWin)'); fvec = (0:nWin-1) * Fs / nWin; [~, pos] = min(abs(fvec - f0)); amp(k) = 2 * abs(X(pos)) / nWin; pha(k) = angle(X(pos)); end % 幅值最大者判为故障线路 [~, faultLine] = max(amp);

这里有一个细节：FFT之前要加窗，直接用原始截断数据会产生频谱泄漏，幅值误差可能达到百分之十几。汉宁窗是比较稳妥的选择，它虽然会降低频率分辨率，但基波附近没有其他强分量时影响很小。如果要做五次谐波分析，窗长至少取0.1秒，保证150Hz和50Hz、100Hz各频率点能区分开，缩短窗长会导致谱线间隔变大，谐波分量混叠。

自动选线时，如果只用幅值最大原则，消弧线圈完全补偿或过补偿时，故障线路的零序电流幅值可能不是最大的，这时幅值法会直接选错。这就是为什么后面要加暂态能量、首半波极性等辅助判据。用脚本把多判据结果合并，再输出一个综合推荐结论，比单判据仿真有更好的参考价值。

4.4 多判据横向对比与阈值整定

阈值整定是选线算法落地的关键。稳态幅值法的启动阈值可以按照正常运行时最大不平衡零序电流的3到5倍来设，仿真中则可以用非故障线路零序电流最大幅值的1.2到1.5倍作为故障线路判定门槛。这个系数不能拍脑袋，而是要根据仿真中过渡电阻的变化趋势来定：金属性接地时故障线路幅值是健全线路的3倍以上，高阻接地时可能只有1.5倍，阈值设太高会漏选高阻接地，设太低又会误选健全线路。

首半波法主要看第一个半波的极性，故障线路和健全线路的首半波极性相反。仿真中判断极性指的是故障发生后的暂态电流从负到正还是从正到负，建议加一个很小的高频带通滤波，把工频和直流分量都滤掉，只观察暂态振荡的第一波。故障发生在电压过零附近时，暂态能量非常小，首半波极性很容易被噪声干扰，这种情况要和稳态法结果取“或”的关系，而不是“与”的关系。

能量法判据中，直接在MATLAB里把U0和I0信号相乘，再用trapz函数做数值积分，得到的是故障路径上零序能量函数的极性。健全线路的极性为负，故障线路的极性为正，且幅值明显大于健全线路。需要留意的是零序电压互感器开口三角的极性相位关系，仿真模型阶段如果不一致，能量法极性会整个反转。

5. 常见问题与排查技巧实录

5.1 零序电流方向“不认账”的原因

仿真结果里最常出现的诡异现象是：故障线路和非故障线路的零序电流方向看起来差不多，或者时正时负，完全没有理论上的“反相”关系。排查时先确认所有电流测量模块的参考方向是否统一，再确认故障线路的电流是否存在电压互感器极性反转。另一个常见原因是过渡电阻设置过大，金属性接地变成几百上千欧的高阻接地后，零序电流幅值很小，相位受负荷和线路不对称影响很大，方向关系不再稳定。

还有一个很隐蔽的因素：如果故障点离母线比较远，线路的对地电容分布在多个PI段上，实际测到的是线路首端到故障点之间这一段的对地电容电流，和理论推导里“整条线路电容电流”并不完全一致。所以仿真中线路用多少个PI段，对相位精度是有影响的。把每条线路拆成三个PI段以上，相位误差会明显降低。

5.2 模型数值震荡与波形毛刺的处理

故障瞬间断路器动作或故障模块切入，电路拓扑突变，数值求解器会产生高频数值震荡。如果示波器上看到尖峰毛刺像“狗啃过”一样，首先看电源内阻抗是否设成了0，然后看步长是否过大。把求解器切换到ode23tb、减小最大步长，通常能解决大部分问题。还可以在母线处加一个小的对地电容，模拟实际的母线杂散电容，这能让电压突变变得平缓一点，但要注意加电容不能太大，否则会改变系统的电容电流分布。

消弧线圈模型里如果只用纯电感，切除故障或故障点熄弧后会出现持续性的低频振荡波形，这会严重干扰首半波法和能量法的判据。给消弧线圈并上阻尼电阻，或者在电感支路串联一个小电阻，振荡衰减速度就会正常。

5.3 中性点引出错误与接地变压器缺失

经消弧线圈接地系统里，如果系统没有可用的中性点（例如全部是三角形接线），必须先加接地变压器。仿真中很多人直接把消弧线圈挂在一个不相干的节点上，结果系统根本不存在零序通路，仿真结果自然乱七八糟。接地变压器在Simulink里可以用Zigzag Transformer模块建模，也可以用一个Yd接线的变压器来构造，但要注意容量选择，必须能够承受消弧线圈的补偿电流。

如果系统本身是Yg接线的变压器，也要确认中性点是否允许引出。继电保护仿真的严谨性在于一次系统零序网络要闭合，否则零序电压、零序电流都失去物理意义。出现“零序电压很大但每条线路零序电流都是0”这种异常结果时，先检查零序网络是否闭合。

5.4 电压互感器铁磁谐振导致的误判

实际系统中电压互感器的励磁电感会和线路对地电容构成谐振回路，故障后可能发生铁磁谐振，表现为零序电压异常升高、三相电压数值紊乱。仿真里如果挂接了带非线性磁化特性的PT模型，也可能复现这种现象。遇到这个情况，选线算法会收到严重干扰，因为零序电压这个基准量本身已经失真。

仿真阶段建议用线性电压互感器建模，避免非线性谐振干扰掩盖选线算法本身的问题。但实际上现场确实会发生铁磁谐振，所以工程判据里通常会加延时或电压互感器消谐措施，这个复杂工况我在仿真里会单独作为“骚扰场景”来测试，而不是混在主故障工况里。

5.5 常见问题速查表

最后，聊几点实操中的体会

模型搭得再漂亮，最终还是要回到现场录波数据上去验证。我自己做下来的最大一个感受是：别一上来就追求小波神经网络之类的高级算法，先把稳态幅值、相位方向、五次谐波这几套基础判据在Simulink里“跑顺”，再叠加暂态能量和首半波分析，工业项目里的综合选线策略基本都是这个递进套路。另一个经验是把每个工况的仿真模型参数、故障设置、导出数据和后处理脚本打包保存，因为给别人复现“为什么某个判据在这个场景下选错了”时，这些原始记录比任何结论都更有说服力。做仿真研究的人最容易忽略的就是“可复现性”，而选线研究恰恰是一场需要反复对照工况的持久战。如果你也在搭小电流接地系统的Simulink选线模型，希望这份记录能帮你少踩几个我踩过的坑。