简介这份文档面向电气工程及其自动化专业学生、配电网运维与继电保护技术人员围绕366kV中低压配电网广泛采用的小电流接地系统梳理单相接地的故障机理与分析思路。内容先比较中性点不接地与经消弧线圈接地两种接线方式的基本原理、运行特点及参数差异再分析故障时稳态与暂态电气量的变化包括故障相电压降为零、非故障相电压升为线电压、各相对地电容电流与系统零序电压情况并借助Matlab/Simulink搭建仿真模型分别对两种接地方式给出仿真结果与线路主要参数波形。文件包仅含1个doc文档压缩包约4.58MB属于图文与公式并重的论文式资料便于按章节阅读和摘录。目前已有229人学习适合用来理解80%以上配网故障的成因、掌握仿真建模步骤并整理故障选线困难的原因与后续注意问题。1. 366kV 配电网为什么允许带故障运行两小时调度台最常接到的一类电话是某条 10kV 出线 A 相掉在地上了用户那头灯都没闪一下。值班员切到母线电压画面看三相对地电压里 A 相归零B 相和 C 相各涨到 10kV 线电压值但三块线电压表纹丝不动还是 10kV 且互差 120°。这就是小电流接地系统单相接地故障最反直觉的地方——故障相电压消失、非故障相电压升高到原来的 √3 倍线电压却保持完全对称所以规程允许继续运行 12 小时让运维去找那条故障出线。我国 366kV 中低压配电网绝大多数采用中性点非有效接地方式单相接地故障占全网故障的 80% 以上故障电流却只有几安培到几十安培混杂在上百安培的负荷电流里。要在这点信噪比下判断哪条出线接地难点不在保护整定而在故障特征的提取。这份资源要解决的就是这条链路先把中性点不接地和经消弧线圈接地两种方式下的稳态、暂态电气量算清楚再用 Matlab/Simulink 把波形复现出来最后从波形里提炼可用的选线判据。做配网自动化、继电保护选线装置或者正在做电气工程毕设的人都能直接拿去复现。2. 中性点接地方式选型不接地与经消弧线圈接地的参数边界2.1 有效接地与非有效接地的判据工程上区分大电流接地和小电流接地不看中性点接了什么东西而看零序阻抗与正序阻抗的比值。零序电抗与正序电抗之比 X0/X1 ≤ 3且零序电阻与正序电抗之比 R0/X1 ≤ 1属于中性点有效接地单相接地就构成低阻短路回路保护必须立即跳闸反过来 X0/X1 3、R0/X1 1就是中性点非有效接地也就是常说的小电流接地系统接地电流小允许带故障运行。这四种方式在参数上的差异决定了后面选线和建模仿真时的边界条件接地方式接地电流量级故障后是否跳闸主要限制中性点直接接地数千安至数十千安立即跳闸供电可靠性低跳闸次数多中性点经小电阻接地1001000A立即跳闸需配零序保护电缆网常用中性点不接地几安至几十安可运行 12h弧光过电压选线困难中性点经消弧线圈接地补偿后残流 110A可运行 12h补偿度失配会谐振注意中性点不接地并不是真的什么都不接而是经过变压器中性点的等效对地电容接地这个电容很小所以中性点对地阻抗极大零序阻抗才大。2.2 中性点不接地系统的电容电流估算与规程门槛选不选消弧线圈第一步是把系统的单相接地电容电流算出来。架空线和电缆的估算公式差别很大架空线每相对地电容小10kV 架空线大约每 30km 产生 1A 的电容电流电缆因为相间和相对地距离近同样的长度电容电流要大一个数量级。常见做法是按线路类型分别估算再求和% 中性点不接地系统单相接地电容电流估算 % 输入线电压等级、各类线路总长度 Ue 10; % 线电压单位 kV L_oh 45; % 架空线总长度单位 km L_cb 6; % 电缆总长度单位 km k_oh 3.0; % 架空线经验系数一般取 2.7~3.3 k_cb 0.1; % 电缆线路系数单位 A/(kV*km) Ic_oh k_oh * Ue * L_oh * 1e-3; % 架空线电容电流A Ic_cb k_cb * Ue * L_cb; % 电缆电容电流A Ic Ic_oh Ic_cb; % 全系统电容电流A fprintf(架空线电容电流 %.2f A\n, Ic_oh); fprintf(电缆电容电流 %.2f A\n, Ic_cb); fprintf(系统总电容电流 %.2f A\n, Ic);这段脚本里的 k_oh 是关键参数它随导线半径、几何均距、对地高度变化架空线取 2.7 偏保守、取 3.3 偏激进实际工程按同区域已投运线路的实测值反推更靠谱。规程门槛是310kV 配电网单相接地电容电流大于 10A、35kV 大于 10A部分区域按 30A 执行时应装设消弧线圈。按上面这组参数算出 Ic 大约 27A已经远超门槛中性点不接地方案在这个规模下就不合适了——一旦发生弧光接地故障点电流不足以自行熄弧电弧会间歇性重燃。2.3 消弧线圈补偿度、脱谐度与电感值计算消弧线圈的本质是一个可调电感接在变压器中性点与地之间用电感电流去抵消系统的电容电流。补偿度 p 定义为 (I_L - I_C)/I_C脱谐度 v 定义为 (I_C - I_L)/I_C。全补偿 p0 会让电感与对地电容在工频下谐振中性点位移电压急剧升高实际运行中禁止欠补偿 p0 在系统线路退出、电容电流下降时会滑向全补偿也危险工程上普遍采用过补偿p 取 5%10%。% 消弧线圈电感值与脱谐度计算 Ue 10e3; f 50; w 2*pi*f; Ic 27.0; % 上一步算出的系统电容电流A p 0.08; % 过补偿度工程常用 5%~10% IL (1 p) * Ic; % 消弧线圈应提供的补偿电流A Up Ue / sqrt(3); % 相电压有效值V XL Up / IL; % 消弧线圈感抗ohm L XL / w; % 消弧线圈电感H v (Ic - IL) / Ic; % 脱谐度过补偿时为负值 fprintf(补偿电流 %.2f A感抗 %.2f ohm电感 %.4f H脱谐度 %.2f%%\n, ... IL, XL, L, v*100);代入 Ic27A、p0.08 算下来补偿电流约 29.2A感抗约 197.8Ω电感约 0.63H脱谐度约 -8%。这三个数直接决定 Simulink 模型里消弧线圈支路的参数千万不要在仿真里随手填一个电感值——感抗偏差 10%稳态残流的方向和幅值就会跟着偏后面首半波判据的验证结果全不可信。另外要注意消弧线圈通常是分接头可调的仿真时最好把当前档位对应的电感和系统实际电容电流分别设为变量方便扫参数。3. 单相接地故障稳态与暂态电气量的定量特征3.1 稳态零序电压与各相对地电压的相量关系中性点不接地系统正常运行时三相对地电容对称中性点对地电压近似为零三相对地电压就是三相相电压零序电压为零。A 相 d 点金属性接地后该点电位被强行拉到地电位各相对地电压变成相电动势减去 E_AU_A 0U_B E_B − E_AU_C E_C − E_A。按相量算E_B − E_A 的模是 √3 倍相电压相位滞后 E_A 150°E_C − E_A 同理超前 150°两个非故障相对地电压都升高到线电压。零序电压按定义取三相对地电压的平均值U_0 (0 E_B − E_A E_C − E_A)/3 −E_A。幅值等于相电压方向与故障相电动势相反。母线电压互感器开口三角绕组输出的是 3U_0二次侧额定相电压 57.7V 时开口三角满量程就是 100V。所以现场看到的现象是接地相电压表归零另外两相涨到线电压开口三角从几伏跳到接近 100V。这三个量同步出现基本可以判定发生了单相接地剩下的才是选线。3.2 零序电流分布故障线路与非故障线路的极性差异区分故障线路和非故障线路靠的是零序电流的分布规律。母线上每一条非故障出线的零序电流只包含本线路自身的对地电容电流方向由母线流向线路故障出线的零序电流则是全系统所有非故障元件对地电容电流的总和方向由线路流向母线。两者幅值上差一个本线路电容电流 vs 系统总电容电流的量级相位上差约 180°。判据维度非故障线路故障线路中性点不接地故障线路过补偿消弧线圈零序电流幅值自身电容电流小全系统电容电流之和最大残流可能不是最大零序电流方向母线→线路线路→母线过补偿时方向可能反转基波极性法基准与非故障线路反相失效需换判据五次谐波极性法基准反相且幅值最大仍可用消弧线圈对 250Hz 补偿弱暂态首半波法基准首半波极性相反、幅值大仍可用这张表里最关键的一行是过补偿时方向可能反转。中性点不接地系统里故障线路零序电流幅值最大且极性相反这两条同时成立一旦中性点经消弧线圈且过补偿运行电感电流把电容电流盖过去故障线路的稳态基波零序电流方向翻转幅值也不再最大。很多选线装置在纯不接地系统上调好了接到谐振接地站就频繁误判根因就在这里。3.3 暂态过程电容电流自由振荡与电感电流的频段分离稳态信号只有几安培、还受补偿方式影响所以工程上更看重故障发生后的暂态量。暂态过程包含两个分量一个是暂态电容电流由线路等效电感和对地电容构成串联谐振回路自由振荡频率 f 1/(2π√(L_T·C_0))架空线一般在 3003000Hz电缆线路因为 C_0 大频率偏低另一个是暂态电感电流由消弧线圈电感和对地电容构成回路振荡频率接近工频衰减时间常数大得多。两个分量的频率差了将近两个数量级这是暂态选线能成立的根本原因。故障发生后最初半个周期内电容电流分量占主导故障线路的首半波零序电流幅值最大且与非故障线路的首半波极性相反——这就是首半波法的物理依据。等电容分量衰减完电感分量还在稳态基波判据才接手。仿真建模时如果把求解器步长设得太大3000Hz 的自由振荡分量直接就被数值积分抹掉了波形上只剩工频包络首半波判据自然验证不出来。4. MATLAB/Simulink 建模从线路参数到零序电流波形4.1 SimPowerSystems 元件选型与模型拓扑Simulink 里搭小电流接地系统模型元件选择直接决定仿真结果能不能用。常见做法是用 Simscape Electrical旧版叫 SimPowerSystems里的这几个模块模型环节推荐模块关键参数电源Three-Phase Source线电压、频率、内阻抗主变Three-Phase Transformer (Two Windings)变比、漏抗、接线组别出线Distributed Parameters Line 或 PI Section Line正序/零序 R、L、C线路长度故障点Three-Phase Fault故障相、过渡电阻、投切时间负荷Three-Phase Series RLC Load有功、感性无功消弧线圈Series RLC Branch单相电感值、等效电阻测量Three-Phase V-I Measurement各出线三相电流线路模型的选择是第一个分叉点。PI Section Line 用集中参数 π 型等值计算快但一条出线如果超过几十公里、又要看高频暂态集中参数会把行波过程抹平Distributed Parameters Line 用分布参数能保留暂态高频分量代价是仿真步长必须足够小。我一般先用 PI Section Line 把稳态量对一遍确认电容电流和零序电压对得上再换成分布参数模型看暂态首半波。4.2 用脚本批量配置模型参数手工拖滑块改参数跑两三组就乱了尤其是要扫过渡电阻、故障相、补偿度这些变量的时候。用 set_param 写脚本配置每次仿真前统一赋值结果才可对比% 小电流接地系统 Simulink 模型参数批量配置 model small_current_gnd; % 已建好的模型名 load_system(model); set_param([model /Three-Phase Source], ... PhaseVoltage, 10e3, ... % 线电压有效值 V Frequency, 50); set_param([model /Three-Phase Fault], ... FaultPhase, A, ... % 故障相设为 A 相 FaultResistance, 0.01, ... % 过渡电阻0.01 近似金属性接地 SwitchTimes, [0.06 0.16]); % 0.06s 发生0.16s 切除 set_param([model /Arc_Suppression_Coil], ... Inductance, 0.63, ... % 消弧线圈电感 H按 2.3 节算出 Resistance, 2); % 线圈等效电阻 ohm set_param(model, StopTime, 0.2); % 仿真总时长 0.2s save_system(model);这段脚本的要点在 FaultResistance 和 SwitchTimes 两个参数。过渡电阻设成 0.01Ω 相当于金属性接地零序电压能冲到接近相电压实际配网里架空线掉地、树枝搭线这类故障过渡电阻常在几十欧到几千欧零序电压会明显降低选线难度成倍上升所以扫参数时建议至少跑 0.01Ω、100Ω、1000Ω 三档。SwitchTimes 设成两元素向量前半段是故障发生时刻后半段是切除时刻切除时刻留出 100ms 是为了让暂态分量衰减完、稳态波形稳定下来。4.3 求解器与仿真步长的取值Simscape Electrical 是刚性系统不能拿 ode45 硬跑。常见做法是选 ode23tb 或 ode15s相对误差容限设在 1e-4 到 1e-5 之间。最大步长是这里最容易踩坑的参数要捕捉 3000Hz 的暂态分量按照每个周期至少 20 个采样点算步长必须小于 1/(20×3000) ≈ 16.7µs如果线路模型换成分布参数、对地电容更大自由振荡频率还要往上走步长得再降一个量级一般设 1e-5 s 比较稳妥。Powergui 模块的采样时间要和最大步长对齐否则输出的波形会出现混叠看起来像是故障后电流出现了莫名其妙的低频拍频其实是数值误差。4.4 故障波形后处理零序电流提取与首半波极性比较仿真跑完从示波器里肉眼看波形只能定性判断要做选线判据验证必须把数据导出来算。零序电流的定义是三相电流之和即 3I_0 i_A i_B i_C按出线分别计算再比较故障后首半波的极性% 零序电流提取与首半波极性比较 t out.iabc.time; iabc squeeze(out.iabc.signals.values); % 维度 [N x 3 x nLine] Ts mean(diff(t)); % 仿真步长 s N size(iabc, 1); nL size(iabc, 3); i0 zeros(N, nL); for k 1:nL i0(:,k) iabc(:,1,k) iabc(:,2,k) iabc(:,3,k); % 逐出线求 3I0 end f0 find(t 0.06, 1, first); % 故障发生时刻对应索引 w round(0.01 / Ts); % 取故障后 10ms 作为观察窗 seg i0(f0:f0w, :); polar zeros(1, nL); for k 1:nL [~, im] max(abs(seg(:,k))); % 窗内绝对值最大的点 polar(k) sign(seg(im, k)); % 该点极性即首半波极性 end fprintf(各出线首半波极性); fprintf(%d , polar); fprintf(\n);观察窗取 10ms是因为首半波现象集中在故障后半个工频周期内窗口拉太长会把衰减后的稳态分量混进来极性判断被干扰。另外这段代码取的是窗内绝对值最大点而不是第一个过零点原因是实测和仿真波形都带噪声过零点位置容易抖动取幅值峰值点更稳。如果所有出线极性都一样先别急着改判据优先检查三件事仿真步长是不是太大、故障时刻索引是不是取错了、三相电流测量方向是不是有出线接反了。5. 从仿真波形到选线判据几个容易翻车的细节把波形跑出来只是第一步真正决定这套模型有没有价值的是判据在参数扫动下稳不稳。中性点不接地模型可以固定用基波零序电流幅值加极性判断但谐振接地模型必须准备两套后备判据。第一套是五次谐波法。消弧线圈的感抗与频率成正比在 250Hz 下的感抗是工频的 5 倍但它要补偿的电容电流在 250Hz 下只增大 5 倍——看似刚好抵消实际上消弧线圈的补偿度是按工频整定的五次谐波下补偿严重不足所以五次谐波零序电流的分布规律基本退化成中性点不接地系统的样子故障线路幅值最大、极性相反。提取方式可以用 FFT也可以直接上带通滤波器代码里常用的做法是对故障后稳态段做单频点 DFT。第二套是暂态首半波法不依赖稳态量适合间歇性弧光接地——这种故障的稳态量根本不存在每隔几个周期才重燃一次基波判据完全失效但每次重燃的首半波极性是稳定的。首半波法的窗口别贪长510ms 足够取长了反而被暂态电感分量污染。验证判据时建议固定一张参数扫描表每改一个变量跑一轮观察判据输出是否随之跳变扫描变量取值范围观察指标期望结果过渡电阻0.01 / 100 / 1000 Ω零序电压幅值、故障线路 3I0幅值单调下降极性不变故障相A / B / C开口三角相位随故障相旋转 120°补偿度 p-10% / 0 / 8%故障线路基波 3I0 方向p0 时方向反转故障位置线路首端 / 中段 / 末端首半波幅值首端最大末端最小线路模型PI 型 / 分布参数暂态振荡频率分布参数下频率更高有一点容易被忽略仿真里三相线路参数是对称的实际配网三相换位不完全、沿线负荷不平衡稳态下本来就有几安培的不平衡零序电流。如果仿真中直接用 3I_0 的绝对值做阈值判断到了现场阈值一压就误动。常见做法是在仿真里人为注入 2%5% 的三相不对称看判据还能不能扛住这比把模型调得完美对称再宣布判据有效要靠得住得多。本文还有配套的精品资源点击获取