简介这份文档围绕小电流接地系统单相接地故障展开仿真与分析面向电气工程及其自动化专业的本科生、研究生及配电网运维技术人员帮助解决故障特征理解不清、故障选线困难等实际问题。资源为单个doc文件压缩包约4.58MB内容以毕业论文形式呈现结构完整便于系统学习与参考。文档首先介绍电力系统中性点接线方式对比中性点不接地系统与经消弧线圈接地系统的原理及运行特点随后深入分析单相接地故障时的稳态与暂态电气量并借助Matlab仿真软件搭建故障模型分别对两种接地方式进行仿真给出线路各主要参数的波形图与结论。读者可从中掌握故障相电压降为零、非故障相电压升为线电压等典型特征理解系统可带故障运行1至2小时的原因并获得故障检测与选线的分析思路。目前已有229人学习适合作为课程设计、毕业设计及配电网故障分析的学习素材。1. 小电流接地系统单相接地故障分析为什么中性点不接地反而更难查10kV 配电房里最让人头疼的一幕不是断路器跳闸而是后台报出“母线接地”却没有任何开关动作。小电流接地系统单相接地故障分析说的就是这件事中性点不接地或经消弧线圈接地的中压系统发生单相接地后故障电流只有几安到几十安线电压依然对称设备还能继续跑一到两个小时可故障点在哪条线路、哪一相、离母线多远全靠人去猜。它解决的不是“跳不跳闸”的问题而是“不跳闸的情况下怎么把故障揪出来”。适合配网运维、继保调试、二次设计方向的从业者尤其是手里只有录波文件和零序电流曲线、却要给出结论的人。这篇笔记按“先立理论、再动手复现、最后讲坑”的顺序走所有算例都可以用 Python 直接跑。2. 先搞懂零序电压和零序电流到底从哪来2.1 中性点不接地系统的零序通路正常运行时三相对地电容基本对称流过中性点的电流近似为零。一旦 A 相某处绝缘击穿接地A 相对地电压被拉低到接近零B、C 两相对地电压升高到线电压三相对地电容电流不再平衡这个不平衡电流必须找一条回路流回去。对中性点不接地系统这条回路就是B、C 相电容电流经接地点、经大地、经非故障线路的对地电容返回电源。注意故障线路流过的零序电流是所有非故障线路零序电流之和方向从线路流向母线非故障线路流过的只是自身电容电流方向从母线流向线路。这个方向差异就是后面选线判据的物理基础。零序电压的数值也有讲究。金属性接地时开口三角电压接近 100V对应相电压经过渡电阻接地时会明显下降。工程上常用一个经验关系过渡电阻越大零序电压越低零序电流也越小这也是高阻接地最难查的根本原因。2.2 用一段最小代码把零序量算出来理论讲完必须能算。下面这段代码模拟一个三相对称电源加单相接地用相量法算出零序电压和零序电流参数可以自己改。import numpy as np # 系统参数10kV 系统相电压 5773V U_phase 5773.0 # 相电压有效值(V) f 50.0 # 工频(Hz) omega 2 * np.pi * f # 三相对地电容(单位:F)模拟三条出线 C_lines np.array([1.2e-6, 0.8e-6, 1.5e-6]) # 线路1/2/3 C_total C_lines.sum() # 故障条件A相经过渡电阻 Rf 接地 Rf 500.0 # 过渡电阻(Ω) # 正常三相对地电压相量(以A相为参考) Ua U_phase * np.exp(1j * 0) Ub U_phase * np.exp(-1j * 2 * np.pi / 3) Uc U_phase * np.exp(1j * 2 * np.pi / 3) # 故障后中性点位移电压 Un # 由节点电压法Un -(Ua*Ya Ub*Yb Uc*Yc) / (YaYbYcYf) Ya 1j * omega * C_total Yb 1j * omega * C_total Yc 1j * omega * C_total Yf 1.0 / Rf Un -(Ua * Ya Ub * Yb Uc * Yc) / (Ya Yb Yc Yf) # 故障后各相对地电压 Ua_f Ua Un Ub_f Ub Un Uc_f Uc Un # 零序电压(开口三角)近似为 3*Un 的模 U0 abs(3 * Un) # 各线路零序电流非故障线路为自身电容电流故障线路为其余线路之和 I0_lines [] for i, C in enumerate(C_lines): if i 0: # 假设线路1为故障线路 I0 abs(1j * omega * (C_total - C) * U_phase) else: I0 abs(1j * omega * C * U_phase) I0_lines.append(I0) print(f中性点位移电压 Un {abs(Un):.1f} V) print(f零序电压 3U0 {U0:.1f} V) for i, I in enumerate(I0_lines): print(f线路{i1} 零序电流 {I:.3f} A)这段代码的逻辑说明先按节点电压法求中性点位移电压Un再叠加到各相得到故障后对地电压最后按“故障线路零序电流等于其余线路电容电流之和”的结论算各线路零序电流。参数说明C_lines是各出线对地电容实际工程中可用电缆长度估算电缆约 0.20.3 μF/km架空线约 0.0050.01 μF/kmRf是过渡电阻金属性接地取 110 Ω弧光接地取几百欧高阻接地可到几千欧。把Rf从 10 改到 5000你会看到3U0从接近 100V 掉到十几伏这就是现场“有接地但电压不高”的来源。2.3 消弧线圈接入后判据为什么会失效很多变电站装了消弧线圈目的是补偿电容电流、让电弧自熄。但它一接入故障线路的零序电流里多了一个感性分量和原来的容性分量方向相反。当补偿度接近全补偿时故障线路零序电流可能小到和正常不平衡电流一个量级传统的“零序电流最大”判据直接翻车。常见做法是改用五次谐波分量或首半波极性判据因为消弧线圈对五次谐波的感抗是基波的五倍补偿作用弱故障特征保留得更好。选型时如果现场以电缆为主、电容电流大优先考虑消弧线圈如果以架空线为主、电容电流小中性点不接地加小电流选线装置更经济。3. 用录波数据做单相接地选线的完整流程3.1 数据准备与零序通道对齐现场拿到的一般是 COMTRADE 格式录波包含三相电压、三相电流和零序电压、零序电流通道。第一步不是急着算而是对齐采样率和时标。常见坑是零序电流取自电缆屏蔽层或专用零序 CT和相电流的变比、极性可能不一致。我一般先做两件事确认所有通道采样率一致典型 4kHz 或 12.8kHz再用故障前 5 个周波的数据检查零序通道本底噪声。如果本底噪声超过 0.1A说明零序 CT 接线或屏蔽有问题先处理硬件再谈算法。import numpy as np def load_comtrade_like(data_dict, fs): data_dict: {通道名: 数组}, fs: 采样率(Hz) # 统一转成 numpy 数组 for k in data_dict: data_dict[k] np.asarray(data_dict[k], dtypefloat) return data_dict, fs def preprocess_zero_seq(data, fs, fault_start): 提取故障前后各5周波去直流 n_cycle int(fs / 50) # 一个周波点数 seg data[fault_start - 5*n_cycle : fault_start 5*n_cycle] seg seg - seg.mean() # 去直流分量 return seg逻辑说明load_comtrade_like把各通道统一成浮点数组避免整型溢出preprocess_zero_seq截取故障前后各 5 个周波并去直流去直流是为了防止 CT 饱和或采样偏置导致零序电流虚高。参数说明fault_start是故障起始点可用零序电压突变量超过阈值如 5V来定位fs必须和录波文件一致用错采样率会导致频率计算全错。3.2 基波与五次谐波零序电流提取选线判据的核心是提取故障后一个周波内的零序电流基波幅值和五次谐波幅值。用全波傅里叶或最小二乘拟合都行工程上更常用全波傅里叶计算量小。def dft_amplitude(seg, fs, freq): 计算指定频率的幅值 N len(seg) t np.arange(N) / fs ref np.exp(-1j * 2 * np.pi * freq * t) X np.sum(seg * ref) * 2 / N return abs(X) # 假设已有一组线路的零序电流段 fs 4000 fault_start 2000 I0_line1 np.sin(2*np.pi*50*np.arange(4000)/fs) * 3.0 0.2*np.random.randn(4000) I0_line2 np.sin(2*np.pi*50*np.arange(4000)/fs) * 0.8 0.2*np.random.randn(4000) seg1 preprocess_zero_seq(I0_line1, fs, fault_start) seg2 preprocess_zero_seq(I0_line2, fs, fault_start) I1_base dft_amplitude(seg1, fs, 50) I2_base dft_amplitude(seg2, fs, 50) I1_5th dft_amplitude(seg1, fs, 250) I2_5th dft_amplitude(seg2, fs, 250) print(f线路1 基波{I1_base:.3f}A 五次{I1_5th:.3f}A) print(f线路2 基波{I2_base:.3f}A 五次{I2_5th:.3f}A)逻辑说明dft_amplitude用单频点 DFT 提取指定频率幅值比 FFT 更直接适合只关心基波和五次谐波的场景。参数说明freq取 50 和 250seg长度建议至少一个完整周波太短会频谱泄漏。实际判据可以组合基波幅值最大且五次谐波也最大基本可以锁定故障线路如果基波被消弧线圈补偿掉就看五次谐波谁最大。3.3 用零序电流方向做辅助确认幅值判据在极端情况下会误判比如两条线路电容电流接近。这时加方向判据故障线路零序电流相位与非故障线路相反。实现上比较故障线路和非故障线路零序电流的相位差接近 180° 的为故障线路。def phase_diff(seg1, seg2, fs): 计算两段信号基波相位差(度) N len(seg1) t np.arange(N) / fs ref np.exp(-1j * 2 * np.pi * 50 * t) X1 np.sum(seg1 * ref) X2 np.sum(seg2 * ref) return np.angle(X1 * np.conj(X2), degTrue) pd phase_diff(seg1, seg2, fs) print(f线路1与线路2零序电流相位差 {pd:.1f} 度)逻辑说明用互相关相量法求相位差X1 * conj(X2)的相角就是两信号的相位差。参数说明相位差接近 180° 说明两者方向相反其中幅值大的那条通常是故障线路接近 0° 说明同向一般都不是故障线路。注意相位判据依赖零序 CT 极性统一现场接线不统一时先做极性校正。4. 避坑与排查单相接地分析里最容易翻车的五件事4.1 零序电压有但零序电流全都很小现象后台报接地开口三角电压 60V 以上但每条线路零序电流都不到 0.5A。原因多数是零序 CT 变比选大了或者零序电流回路被短接、开路。解决先查零序 CT 变比和实际一次电流是否匹配10kV 出线零序 CT 常用 50/1 或 100/1如果一次电容电流只有 2A二次侧确实很小需要换小变比或改用高灵敏度选线装置。4.2 消弧线圈一投选线装置就乱报现象消弧线圈投运后原本能正确选线的装置开始频繁误报。原因全补偿状态下故障线路基波零序电流被压到接近零装置按幅值排序必然出错。解决把判据切到五次谐波或首半波同时检查消弧线圈补偿度避免长期运行在全补偿附近通常欠补偿 5%10% 更利于选线。4.3 过渡电阻一大所有判据都失灵现象高阻接地时零序电压只有十几伏零序电流和正常不平衡电流混在一起。原因故障特征随过渡电阻增大而衰减信噪比太低。解决不要只依赖稳态量改用暂态量判据比如首半波极性或暂态零序电流小波变换模极大值暂态量在高阻接地时仍有可检测分量。现场可加装暂态选线装置作为补充。4.4 录波文件通道名对不上算出来的全是错的现象代码跑通但结果和现场结论相反。原因COMTRADE 配置文件里通道顺序和实际接线不一致或者零序电流通道实际接的是 A 相电流。解决每次分析前用故障前稳态数据核对正常运行时零序电流应接近零如果某个“零序”通道有明显 50Hz 分量说明通道接错。4.5 只看单条线路忽略母线接地可能现象所有出线零序电流都不大但零序电压很高。原因故障可能发生在母线上比如母线 PT 或避雷器绝缘击穿。解决先排除母线故障再查线路。方法是对比各出线零序电流之和与母线零序电流如果母线零序电流明显大于出线之和故障在母线。5. 进阶技巧用暂态首半波做高阻接地选线稳态判据在高阻接地时基本没戏真正能救场的是暂态首半波。故障发生瞬间线路对地电容通过故障点放电产生一个频率几百赫兹到几千赫兹的暂态零序电流这个电流不受消弧线圈补偿影响幅值也比稳态大。具体做法取故障后 12ms 数据比较各线路暂态零序电流的首半波极性和幅值故障线路的首半波幅值最大且极性与非故障线路相反。def first_half_wave(seg, fs, fault_start): 提取首半波特征第一个极值点和极性 n_1ms int(fs * 0.001) seg_short seg[fault_start : fault_start 2*n_1ms] idx np.argmax(np.abs(seg_short)) return seg_short[idx], np.sign(seg_short[idx]) # 对多条线路分别提取 for name, seg in [(线路1, seg1), (线路2, seg2)]: amp, pol first_half_wave(seg, fs, fault_start) print(f{name} 首半波幅值{abs(amp):.3f}A 极性{pol})逻辑说明first_half_wave取故障后 2ms 数据找绝对值最大点作为首半波特征。参数说明fault_start要精确到采样点误差超过 1ms 会错过首半波极性用np.sign判断故障线路和非故障线路应相反。实际使用时建议配合小波变换或带通滤波先滤掉工频分量再提取暂态。验证方法拿一段已知故障线路的录波把稳态判据和暂态判据都跑一遍对比选线结果。如果暂态判据在高阻接地案例中正确率明显高于稳态就值得在现场加装暂态选线模块。我自己的习惯是每次分析完都把录波、判据结果和现场巡线结论存成一个案例库下次遇到类似波形先翻案例库比从头算快得多。希望帮到你。本文还有配套的精品资源点击获取