简介这份文档面向电气工程及其自动化专业的本科生、研究生及配电网运维技术人员围绕366kV中低压配电网广泛采用的小电流接地系统系统讲解单相接地故障的稳态与暂态特征、中性点不接地与经消弧线圈接地两种方式的原理差异以及基于Matlab的仿真建模与结果分析方法。资源包内含1个doc文件约4.58MB为完整的毕业设计论文涵盖摘要、中性点接线方式比较、故障电气量分析、仿真模型搭建、波形图展示与结论总结等章节可直接用于课程设计、毕业设计选题参考或故障选线研究入门。目前已有229人学习下载。读者可从中获取故障相电压降为零、非故障相电压升为线电压等关键规律的推导过程掌握利用Simulink设置线路参数、复现两种接地系统仿真波形的具体思路并借鉴论文的章节组织与分析方法为配电网单相接地故障定位与选线研究提供可复用的参考框架。1. 小电流接地系统单相接地故障分析为什么它是配网里最容易被低估的“灰犀牛”配电网里有一种故障保护不跳闸、设备不冒烟、值班员盯着屏幕只能看到母线电压轻微晃动但线路下方可能已经有电缆在慢慢烧蚀、绝缘子在持续爬电。这就是小电流接地系统中性点不接地或经消弧线圈/小电阻接地发生单相接地后的典型场景。它不像相间短路那样“炸裂”却占了配网故障总数的六成以上而且选线不准、定位不到、拖成永久故障的案例比比皆是。这篇笔记面向做配网自动化、继电保护整定、故障选线算法或运维排故的工程师从建模、特征提取、选线判据到现场排查把“小电流接地系统单相接地故障分析”这件事拆成能复现的步骤。新手能照着搭模型跑波形熟手能看到判据边界和参数陷阱。2. 先搞清楚故障特征从哪来中性点接地方式决定了你看到什么2.1 三种接地方式下单相接地到底发生了什么小电流接地系统的核心特征是中性点非有效接地单相接地后故障电流被限制在很小的量级线电压三角形基本不变所以系统可以带故障运行一到两小时。但“能带故障运行”不等于“可以不管”因为非故障相电压会升高到线电压绝缘薄弱点随时可能被击穿成两相短路。中性点不接地时单相接地电流等于全网对地电容电流的矢量和数值通常在几安到几十安。中性点经消弧线圈接地时电感电流补偿电容电流故障点残流可以压到5安以下但补偿方式过补偿、欠补偿、全补偿直接决定了零序电压和零序电流的相位关系。中性点经小电阻接地时故障电流被限制在几百安零序保护可以直接动作跳闸但已经不算是“小电流”了很多选线算法在这个场景下会失效。理解这一点很关键你拿到的零序电流波形是全网对地电容分布、消弧线圈补偿度、故障点过渡电阻共同作用的结果不是一条“标准正弦”。很多选线算法在仿真里准确率95%到了现场掉到60%根子就在这里。2.2 零序电压和零序电流的相位关系选线判据的物理基础故障选线最常用的判据是“零序电流幅值最大”和“零序电流方向与零序电压的关系”。前者在金属性接地时好用后者在经消弧线圈接地时更可靠。原理是这样的非故障线路的零序电流是本线路对地电容电流方向从母线流向线路故障线路的零序电流是所有非故障线路电容电流之和方向从线路流向母线。两者相位差接近180度。但注意这个关系只在基波分量上成立。现场波形里混着大量5次、7次谐波还有消弧线圈产生的感性电流直接做FFT取基波幅值比较有时候故障线路的基波幅值反而不是最大的。我一般会同时看基波幅值、5次谐波幅值和首半波极性三个判据投票比单一判据稳得多。2.3 用Python搭一个可复现的仿真模型下面这段代码用简单的集中参数模型模拟一条母线上三条出线其中线路2发生单相接地生成零序电压和三条线路的零序电流波形。目的是让你能快速验证选线判据不用等现场录波。import numpy as np import matplotlib.pyplot as plt # 参数设置50Hz系统采样率10kHz仿真0.2秒 fs 10000 t np.arange(0, 0.2, 1/fs) f0 50 w0 2 * np.pi * f0 # 三条线路的对地电容单位微法线路2为故障线路 C np.array([0.8e-6, 1.2e-6, 0.6e-6]) # 消弧线圈电感按过补偿8%设置 L 1.0 / (1.08 * w0**2 * np.sum(C)) # 故障发生时刻0.06秒 t_fault 0.06 fault_idx np.searchsorted(t, t_fault) # 零序电压故障后非故障相电压升高零序电压为相电压 U_phase 220 / np.sqrt(3) # 相电压有效值 u0 np.zeros_like(t) u0[fault_idx:] U_phase * np.sqrt(2) * np.sin(w0 * t[fault_idx:]) # 各线路零序电流非故障线路为自身电容电流故障线路为其余线路电容电流之和 i0 np.zeros((3, len(t))) for k in range(3): if k 1: # 故障线路 i0[k, fault_idx:] -np.sum(np.delete(C, k)) * np.gradient(u0[fault_idx:], 1/fs) else: i0[k, fault_idx:] C[k] * np.gradient(u0[fault_idx:], 1/fs) # 叠加5次谐波模拟现场波形畸变 for k in range(3): i0[k, fault_idx:] 0.05 * np.max(np.abs(i0[k, fault_idx:])) * np.sin(5 * w0 * t[fault_idx:]) # 绘图 fig, axes plt.subplots(2, 1, figsize(10, 6)) axes[0].plot(t, u0) axes[0].set_ylabel(零序电压 (V)) axes[0].set_title(单相接地零序电压波形) for k in range(3): axes[1].plot(t, i0[k], labelf线路{k1}) axes[1].set_ylabel(零序电流 (A)) axes[1].set_xlabel(时间 (s)) axes[1].legend() axes[1].set_title(三条线路零序电流对比) plt.tight_layout() plt.show() # 简单选线判据故障后取一个周波比较基波幅值 start fault_idx int(0.02 * fs) end start int(0.02 * fs) amp [np.max(np.abs(i0[k, start:end])) for k in range(3)] print(f各线路零序电流幅值: {amp}) print(f幅值最大线路: {np.argmax(amp) 1})这段代码的逻辑说明用集中电容模型模拟对地电容电流消弧线圈按过补偿设置故障线路的零序电流方向取反。参数说明C数组控制各线路电容改变它就能模拟不同线路长度L由补偿度决定过补偿8%意味着电感电流比电容电流大8%t_fault控制故障时刻。运行后你会看到故障线路的零序电流幅值明显大于非故障线路但叠加5次谐波后如果只看峰值可能误判所以代码最后用基波幅值比较。实际现场录波里谐波含量更高建议先带通滤波再比较。3. 选线判据怎么选、参数怎么设从幅值法到群体比幅比相3.1 幅值法、方向法、谐波法的适用边界幅值法最简单比较各线路零序电流基波幅值最大者判为故障线路。它在中性点不接地且线路长度差异不大时准确率不错但遇到以下情况会翻车故障线路很短、非故障线路很长非故障线路的电容电流反而更大经消弧线圈接地且接近全补偿时故障线路基波幅值可能小于非故障线路。方向法比较零序电流与零序电压的相位故障线路零序电流滞后零序电压约90度容性或超前90度感性非故障线路相反。方向法在消弧线圈接地时更稳但需要精确的相位测量现场CT极性接反、采样不同步都会导致误判。谐波法利用故障线路5次谐波电流幅值最大的特征因为消弧线圈对谐波补偿效果差谐波电流主要走故障线路。谐波法在经消弧线圈接地时表现好但要求采样率足够高且谐波源如变频器多的线路容易受干扰。我一般会做群体比幅比相先按幅值排序取前三条线路再比较它们的相位关系故障线路的相位应该与其他线路相反。这样比单一判据稳计算量也不大。3.2 用零序导纳法做选线一个可落地的判据实现零序导纳法把每条线路的零序电流除以零序电压得到导纳故障线路的导纳实部为负电流流出母线非故障线路为正。这个方法对CT极性不敏感也不需要知道系统电容电流总量适合现场改造项目。import numpy as np def select_fault_line(u0, i0_list, fs10000, f050): 零序导纳法选线 u0: 零序电压数组 i0_list: 各线路零序电流数组列表 fs: 采样率 f0: 工频 返回故障线路索引从0开始 # 取故障后一个周波的数据 n len(u0) start int(0.02 * fs) end start int(0.02 * fs) u u0[start:end] # 用FFT提取基波相量 def phasor(x): X np.fft.fft(x) k int(f0 * len(x) / fs) return X[k] * 2 / len(x) U phasor(u) Y_real [] for i0 in i0_list: I phasor(i0[start:end]) Y I / U if abs(U) 1e-6 else 0 Y_real.append(np.real(Y)) # 导纳实部为负的线路判为故障线路 fault_idx int(np.argmin(Y_real)) return fault_idx, Y_real # 用上一节的仿真数据测试 # 假设 u0 和 i0 已经生成 # fault_idx, Y_real select_fault_line(u0, [i0[0], i0[1], i0[2]]) # print(f故障线路: {fault_idx 1}, 各线路导纳实部: {Y_real})逻辑说明先取故障后一个周波的稳态数据用FFT提取基波相量然后计算每条线路的零序导纳。参数说明start和end控制取窗位置必须避开故障暂态过程phasor函数里的k是基波对应的FFT bin采样率和窗长变了要重新算。注意如果零序电压很小比如金属性接地时接近零导纳法会失效这时候要切回幅值法或谐波法。3.3 消弧线圈补偿度对选线的影响过补偿、欠补偿、全补偿消弧线圈的补偿度直接改变零序电流的分布。过补偿时电感电流大于电容电流故障线路的零序电流呈感性方向与非故障线路相同幅值法完全失效必须用方向法或谐波法。欠补偿时电感电流小于电容电流故障线路零序电流仍呈容性幅值法还能用但残流较大电弧不易熄灭。全补偿时电感电流等于电容电流故障点残流接近零零序电压最高但选线最困难因为故障线路的零序电流几乎为零。现场整定一般选过补偿补偿度5%到10%目的是保证故障点残流呈感性电弧容易熄灭同时避免谐振。但过补偿度不能太大否则故障线路零序电流太小选线灵敏度不够。我见过一个站消弧线圈过补偿度设到20%结果单相接地时故障线路零序电流只有0.3安选线装置直接拒动。后来调到8%问题解决。4. 现场排查与录波分析从选线装置报警到找到故障点4.1 选线装置报了线路但找不到故障点先查这四处选线装置报警只是第一步真正找到故障点往往要花几个小时。常见情况是选线装置报了某条线路但巡线人员走了一遍没发现明显故障。这时候先查四处第一选线装置的零序电压门槛是不是设得太低导致把不平衡电流误判为接地第二故障线路的零序CT极性是不是接反了方向法会因此选错第三故障点是不是在电缆沟里或者架空线绝缘子表面肉眼看不到明显烧蚀第四是不是间歇性弧光接地故障点已经暂时恢复绝缘。我一般会先调录波看零序电压的波形。如果是间歇性弧光零序电压会呈现明显的“熄弧-重燃”特征电压波形有平顶或尖峰。如果是稳定接地零序电压是标准正弦。然后看故障线路的零序电流如果幅值很小但5次谐波明显说明消弧线圈补偿度较高故障点残流小故障点可能在电缆中间接头或绝缘子内部。4.2 用录波数据反推故障点位置行波法在小电流接地里的局限行波法测距在输电线路里很成熟但在小电流接地系统里用起来很别扭。原因是单相接地故障电流小行波信号弱而且配网线路分支多、波阻抗不连续行波在分支点会多次反射波形识别困难。我试过用行波法在一条10kV线路上测距误差超过500米还不如人工巡线。更实用的方法是“零序电流突变量法”在变电站出线处和线路分段开关处同时录波比较零序电流的突变时刻故障点就在突变时刻差异最大的两个测点之间。这个方法要求各测点时间同步误差不超过1微秒。如果没有同步时钟可以用GPS授时模块成本不高。4.3 一个现场录波分析脚本从COMTRADE文件到选线结论现场录波装置一般输出COMTRADE格式文件包含通道配置和数据。下面这段代码读取COMTRADE的CFG和DAT文件提取零序电压和零序电流通道做选线判断。import numpy as np def read_comtrade(cfg_file, dat_file): 读取COMTRADE文件返回采样率、通道名、数据矩阵 简化版只处理ASCII格式 with open(cfg_file, r) as f: lines f.readlines() # 第一行站名,装置名,版本 # 第二行通道总数,模拟通道数,数字通道数 parts lines[1].strip().split(,) n_channels int(parts[0]) n_analog int(parts[1]) # 读取通道名 ch_names [] for i in range(2, 2 n_analog): ch_names.append(lines[i].strip().split(,)[1]) # 读取采样率 # 简化处理假设第一行采样率信息在后面的行 fs 10000 # 默认值实际应从CFG读取 # 读取DAT文件 data np.loadtxt(dat_file, delimiter,) # 第一列是采样序号后面是通道数据 t data[:, 0] values data[:, 1:1n_analog] return fs, ch_names, t, values def analyze_fault(fs, ch_names, t, values): 从录波数据中找零序电压和零序电流通道做选线 # 找零序电压通道通道名含U0或3U0 u0_idx None i0_indices [] for i, name in enumerate(ch_names): if U0 in name or 3U0 in name: u0_idx i if I0 in name or 3I0 in name: i0_indices.append(i) if u0_idx is None or len(i0_indices) 0: print(未找到零序电压或零序电流通道) return u0 values[:, u0_idx] i0_list [values[:, i] for i in i0_indices] # 找故障起始点零序电压突变 threshold 0.1 * np.max(np.abs(u0)) fault_start np.where(np.abs(u0) threshold)[0] if len(fault_start) 0: print(未检测到故障) return start fault_start[0] # 取故障后一个周波 end min(start int(0.02 * fs), len(u0)) u_seg u0[start:end] # 比较各线路零序电流幅值 amps [] for i0 in i0_list: amps.append(np.max(np.abs(i0[start:end]))) fault_line np.argmax(amps) print(f故障起始点: {t[start]:.4f}s) print(f各线路零序电流幅值: {amps}) print(f选线结果: 线路{fault_line 1}) # 使用示例 # fs, ch_names, t, values read_comtrade(fault.cfg, fault.dat) # analyze_fault(fs, ch_names, t, values)逻辑说明先读CFG文件拿通道名再读DAT文件拿数据然后按通道名找零序电压和零序电流。参数说明threshold控制故障起始点检测灵敏度设大了会漏掉高阻接地设小了会把不平衡电流当故障。实际使用时CFG文件里的采样率信息要正确解析不同厂家的COMTRADE格式有差异可能需要调整解析逻辑。5. 避坑指南小电流接地选线里最容易翻车的五个地方5.1 零序CT极性接反方向法全盘皆输现象选线装置频繁报某条线路接地但现场检查该线路绝缘正常其他线路反而有放电痕迹。原因零序CT的P1、P2方向接反导致该线路零序电流相位翻转180度方向法把它判为故障线路。解决在停电状态下用一次电流法校验零序CT极性或者用选线装置的“极性自检”功能。如果已经送电可以临时把方向法切换为幅值法先保证不误跳。5.2 消弧线圈补偿度设得太高故障线路零序电流被“淹没”现象单相接地后选线装置不报警或者报警线路与实际故障线路不符。原因消弧线圈过补偿度太大故障点残流很小故障线路零序电流幅值低于非故障线路。解决把补偿度调到5%到10%同时启用谐波法或导纳法作为辅助判据。如果消弧线圈是自动调谐的检查调谐控制器是不是按“残流最小”原则工作应该改成“残流感性且选线灵敏度足够”。5.3 采样不同步导致相位比较失效现象方向法选线结果随机跳变同一故障多次录波选出不同线路。原因各线路零序电流采样不同步相位差引入额外误差。解决选线装置必须用同一套AD同步采样或者用GPS对时。如果现场是多台装置分布式录波时间同步误差要控制在1微秒以内否则相位比较没有意义。5.4 间歇性弧光接地被当成稳定接地处理现象选线装置报了故障线路但巡线人员到现场时故障已经消失线路恢复送电后再次接地。原因间歇性弧光接地的零序电压和电流是断续的选线装置如果只取一个周波做判断可能刚好取到熄弧期间的数据。解决选线算法要连续判断多个周波或者用“零序电压突变量启动多次投票”策略。现场排查时重点检查电缆中间接头、绝缘子污秽区、避雷器计数器动作情况。5.5 高阻接地时零序电压低于门槛选线装置拒动现象线路已经接地但选线装置不启动值班员通过其他现象如电压不平衡才发现。原因故障点过渡电阻大比如树线放电、水泥杆塔放电零序电压只有几百伏低于装置启动门槛。解决降低零序电压启动门槛但要注意躲过正常不平衡电压。我一般会把门槛设在15%相电压左右同时增加零序电流突变量启动作为补充。如果还是不行只能靠人工巡视和红外测温。6. 进阶技巧用暂态首半波极性做快速选线以及一个我踩过的坑暂态首半波极性法利用故障发生瞬间的暂态零序电流首半波方向来判断故障线路。故障线路的首半波零序电流方向与非故障线路相反而且暂态分量幅值比稳态大几倍信噪比高。这个方法在消弧线圈接地系统里特别有用因为消弧线圈对暂态高频分量几乎不补偿。实现上取故障后1到2毫秒的数据做小波变换或者带通滤波500Hz到3kHz提取首半波极性。下面是一个简化实现import numpy as np from scipy.signal import butter, filtfilt def first_half_wave_polarity(u0, i0_list, fs10000): 暂态首半波极性法选线 u0: 零序电压 i0_list: 各线路零序电流 fs: 采样率 # 带通滤波 500-3000Hz b, a butter(4, [500/(fs/2), 3000/(fs/2)], btypeband) # 找故障起始点 threshold 0.1 * np.max(np.abs(u0)) start np.where(np.abs(u0) threshold)[0][0] # 取故障后2ms数据 end min(start int(0.002 * fs), len(u0)) polarities [] for i0 in i0_list: i_filt filtfilt(b, a, i0[start:end]) # 首半波极性第一个过零点的方向 sign np.sign(i_filt[np.argmax(np.abs(i_filt))]) polarities.append(sign) # 极性与其他线路相反的判为故障线路 # 简化取极性与众数不同的线路 from collections import Counter count Counter(polarities) majority count.most_common(1)[0][0] fault_idx [i for i, p in enumerate(polarities) if p ! majority] return fault_idx, polarities # 使用示例 # fault_idx, polarities first_half_wave_polarity(u0, [i0[0], i0[1], i0[2]]) # print(f故障线路: {[i1 for i in fault_idx]}, 各线路首半波极性: {polarities})逻辑说明带通滤波提取暂态高频分量取故障后2毫秒内幅值最大点的符号作为首半波极性。参数说明butter的截止频率根据系统暂态频率范围调整配网电缆线路暂态频率通常在1kHz到5kHz架空线略低。filtfilt是零相位滤波避免引入额外相移。注意如果故障发生在电压过零点附近暂态分量很小首半波法可能失效这时候要切回稳态判据。我踩过的一个坑有一次在一个经消弧线圈接地的变电站首半波法选线连续三次都选错了。后来查录波发现故障线路的零序CT安装在电缆屏蔽层接地线上而屏蔽层电流里混入了大量高频干扰首半波极性被干扰淹没。解决办法是把零序CT移到电缆三相合并后的接地线上或者增加模拟滤波电路。这个坑让我明白再好的算法也架不住CT安装位置不对。现场做选线我的习惯是先看零序电压波形判断故障类型稳定/间歇/弧光再看零序电流幅值和相位做初步选线最后用首半波法或谐波法交叉验证。三个方法结论一致才下结论不一致就调录波仔细分析。这套流程帮我避免了好几次误判。希望帮到你。本文还有配套的精品资源点击获取