简介：这份资源围绕小电流接地系统单相接地故障选线主题，面向电力系统相关专业的毕业设计者、研究故障选线算法的学生以及从事继电保护工作的工程技术人员。内容从故障特征分析出发，阐述利用稳态特征与暂态特征结合进行选线的原理，并基于零序电流大小实现故障定位，同时还给出MATLAB仿真模型的构建思路与验证过程，可帮助读者快速理解选线方法、掌握建模仿真流程。压缩包内为1个docx文档，共270KB，虽然文件数量精简，但文档包含摘要、关键词、理论分析、仿真结果等完整章节，能够系统呈现研究脉络。目前已有1112人浏览学习，适合需要撰写毕业设计论文或进行故障选线仿真研究的读者参考。文档将理论分析与仿真验证相结合，既能作为方案设计的依据，也可为后续扩展研究提供思路与模板。
1. 小电流接地系统单相接地故障选线：故障电流只有几安培，为什么保护装置集体失灵？
10kV配网出线柜突然报出接地信号，零序电压越限，选线装置却把结果指向母线段。值班员拉开一条出线，接地信号还在；再拉一条，负荷已经损失了，故障却依然存在。这不是某个特定现场的事故，而是小电流接地系统单相接地故障选线这个老问题每天都在上演的剧本。中性点不接地或经消弧线圈接地的系统里，单相接地故障电流通常只有几安培到几十安培，比负荷电流小一两个数量级，传统过流保护根本启动不了。运维人员能拿到的信号只有零序电压抬升和弥散的零序电流，而这两者在故障前后、在消弧线圈补偿前后、在电弧间歇性燃灭过程中的表现完全不同。
这篇文章把选线这个黑匣子拆开：先讲故障后电气量的真实分布，再给出可复现的仿真模型和算法代码，最后把现场最常见也最容易被忽略的五个坑逐一排掉。无论你是做配电自动化调试的工程师、变电站后台维护人员，还是刚接触继电保护的学生，按这套思路都能从零跑通选线流程，并判断哪些方法在自己的网架结构里真正可靠。
2. 选线为什么难：故障特征量的分布规律与选线方法选型
2.1 单相接地后的电压电流分布：零序电流的真实走向
中性点不接地系统发生单相接地时，故障相电压被拉低，非故障相电压升高到线电压，系统中性点电位发生偏移。每一条非故障线路的对地电容电流从母线流向线路，而故障线路上的零序电流等于全系统所有非故障线路电容电流之和，方向从线路流向母线。这个方向关系是稳态量选线方法的基础，也是很多现场人员理解上最模糊的地方：零序电流不是故障点流出来的“故障电流”，而是系统电容电流在故障相上的回流，必须在开关场把所有出线的零序CT极性统一，并让母线TV开口三角正确反映零序电压相位，方向判据才有效。
经消弧线圈接地的系统里，情况发生根本变化。消弧线圈在接地瞬间向故障点注入感性电流，用于补偿系统的容性电流。当过补偿度取8%到10%时，故障线路的零序电流从原来的容性主导变成感性主导，幅值变小，方向与非故障线路趋于一致。这一个变化直接影响以“零序电流最大”“方向相反”为核心的稳态判据，也是许多老选线装置在消弧线圈投入后正确率暴跌的根本原因。
2.2 稳态量法与暂态量法的适用边界：过渡电阻和电弧稳定性的影响
稳态量法包括零序电流幅值比较、相位比较、零序有功分量法和五次谐波法。幅值比较法在出线数量多、线路长度差异大的情况下容易出现误判，因为长线路的对地电容电流可能比故障线路的零序电流还要大；有功分量法利用故障线路零序有功功率远大于非故障线路这一特征，在中性点经电阻接地系统中效果不错，但在纯不接地系统里没有有功通路，接地点的弧光电阻又极不稳定，测出来的有功值往往只有几瓦到几十瓦，抗干扰能力有限。
暂态量法利用接地瞬间的暂态零序电流。这个暂态过程的频率通常在几百赫兹到几千赫兹之间，幅值可以比稳态量大一个数量级以上，而且不受消弧线圈补偿的遮蔽影响，因为消弧线圈的电感在暂态高频下呈高阻，相当于被自然旁路。首半波法比较的就是故障线路与非故障线路暂态零序电流第一个周波的幅值和方向。它的前提是故障发生在相电压接近峰值时刻，如果接地瞬间相电压接近过零点，暂态能量弱，首半波方向可能被不平衡电流掩盖，这就是首半波法“时灵时不灵”的根源。
2.3 选线方法选型决策表
| 系统接地方式 | 推荐方法 | 方法依赖的关键量 | 可靠性的最大威胁 |
|---|---|---|---|
| 中性点不接地 | 零序电流幅值比较+首半波法融合 | 各出线零序CT幅值、暂态方向 | 长线路电容电流过大、CT变比不一致 |
| 经消弧线圈接地 | 暂态量选线为主，五次谐波辅助 | 暂态零序电流首半波、5次谐波分量 | 过补偿度设定不合理、电弧不稳定 |
| 经电阻接地 | 零序有功分量方向法 | 零序电压与零序电流的相位差 | CT极性接反、消弧电阻阻值漂移 |
选型时的关键判断依据是消弧线圈欠补偿还是过补偿。欠补偿状态下谐振过电压风险高，现场一般不用；过补偿状态下故障线路零序电流方向不再具有区分度，相位比较法直接失效。我一般会先确认系统中性点接地方式和消弧线圈脱谐度，再决定算法核心是走稳态方向还是暂态方向——顺序反了，装置投上去第一天就可能在雷击故障时误动。
3. 把故障特征复现出来：搭建10kV配电网单相接地仿真模型的完整步骤
3.1 仿真模型拓扑与元件参数
没有真实故障录波数据时，验证选线算法最可靠的方式是电磁暂态仿真。常见做法是在仿真环境中搭建一个简化10kV配电网，包含一条母线、三到五条出线、一个中性点接地装置（可切消弧线圈）和一组可调过渡电阻的故障模块。架空线路用集中参数模型即可满足选线算法验证需求，但要注意线路长度、对地电容和零序阻抗的真实比例，不能随便填数字。线路参数按典型钢芯铝绞线取值，单位长度的正序、零序阻抗和对地电容必须分开设置，零序参数误差会直接扭曲仿真中的稳态零序电流分布。
母线侧接三相电压源模拟无穷大系统，出线侧不接负荷或只接少量平衡负荷，因为选线关心的是电容电流分布，负荷电流在这个场景里属于干扰项。消弧线圈建模时用可调电感串联阻尼电阻，电感值根据系统电容电流和目标脱谐度换算。故障模块设置在一条出线中段，用受控开关控制接地投入时刻，过渡电阻从0.1Ω到10kΩ可调，这样就能覆盖金属性接地、经过渡电阻接地和电弧接地三种典型工况。
3.2 设置单相接地故障与数据采集
仿真里最容易被忽略的是采样率和数据窗长度。选线算法要分析暂态首半波，采样率低于20kHz时首半波的波形细节丢失严重，方向判据的可靠性没有说服力。数据窗取故障前30ms到故障后100ms，既能看暂态过程，也能覆盖稳态段，方便同时验证稳态量法和暂态量法。故障合闸时刻选择相电压峰值附近，否则暂态能量不足，首半波特征不明显，这个时间参数要在仿真脚本里明确写出来。
故障相、故障时刻、过渡电阻和故障线路编号每个参数组合都要跑一遍，并把各出线零序电流、零序电压和中性点电流存成独立文件，文件名里带参数标签。我一般按“线路编号+过渡电阻+合闸角度”组织数据文件，比如L2_R500_90代表2号线路在电压峰值时刻经500Ω电阻接地，这样后期做算法对比时不用翻仿真模型就能定位数据来源。
3.3 为什么必须做消弧线圈补偿度扫描
消弧线圈过补偿度从5%到15%之间变化时，故障线路稳态零序电流的幅值和相位可能完全翻转，这是选线算法在谐振接地系统里反复踩坑的根源。仿真中把脱谐度扫描作为标准步骤，每个脱谐度下重新记录零序电流的幅值、相位和首半波波形，你会清楚看到：在过补偿度大于10%后，稳态相位方向翻转，而暂态首半波方向不受影响。这组数据直接决定算法选型，而不是靠记忆中的“消弧线圈投入后要用暂态法”这种笼统结论。
MATLAB或Python环境下的仿真参数脚本可以用一段配置代码概括，下面给一个典型的中性点经消弧线圈接地系统仿真参数设置示例：
# 仿真参数配置示例（可直接用于电磁暂态仿真初始化） grid = { 'voltage_ll': 10e3, # 线电压有效值，单位V 'frequency': 50, # 工频50Hz 'sample_rate': 200e3, # 采样率200kHz，确保暂态高频不混叠 'lines': [ # 每条出线的长度和单位零序电容 {'length_km': 6.3, 'C0_per_km': 5.2e-9}, {'length_km': 4.8, 'C0_per_km': 5.2e-9}, {'length_km': 8.5, 'C0_per_km': 4.8e-9}, # 电缆段较多的出线 {'length_km': 3.2, 'C0_per_km': 5.2e-9}, ], 'arc_suppression_coil': { 'active': True, # True=经消弧线圈接地，False=中性点不接地 'over_compensation': 0.10, # 过补偿度10% 'damping_resistor': 15.0, # 阻尼电阻，单位Ω }, 'fault': { 'faulted_line': 2, # 故障线路索引，从0开始，2代表第3条出线 'fault_phase': 'A', # A相接地 'fault_resistance': [0.1, 100, 1000, 5000], # 扫描过渡电阻 'fault_time': 0.08, # 故障时刻，单位s 'close_angle_deg': 90, # 合闸角90度=相电压峰值处 } }这段配置把仿真中直接影响选线结果的关键参数全部显式化。sample_rate取200kHz是为了给首半波分析留足裕度，实际录波装置采样率往往只有4kHz到10kHz，仿真时用高采样率得到的结果再降采样，模拟真实装置的处理链路，这样算法在真机上运行时才不会有“仿真效果好、现场全不对”的落差。over_compensation设为0.10是过补偿的典型取值，修改这个值就能自动扫描谐振接地的最不利工况。close_angle_deg填90度表示相电压峰值时合闸，首半波暂态能量最大，适合验证算法上限；要测试算法下限，把这个值调到0或180度附近。
4. 从波形到判据：零序电流特征提取与选线算法实现
4.1 零序电流的提取与滤波预处理
现场装置输入的通常是三相电流，零序电流由软件计算得到。工程上常见做法是直接用三相电流之和i_a + i_b + i_c，但CT饱和、测量噪声和三相不平衡会在零序通道里引入很大的虚假分量。更稳妥的做法是使用零序CT的二次侧信号，同时在进入判据前做带通滤波。带通滤波器的下限频率设为50Hz附近以保留稳态基波分量，上限频率设为3000Hz以保留暂态高频信息，带宽设计要同时满足稳态量法和暂态量法的需求，不能为了滤噪声把暂态分量一起滤掉。
滤波器的相频特性对首半波方向判据影响极大。普通二阶巴特沃斯滤波器在截止频率附近有约90度的相位延迟，如果只做一次滤波，首半波的方向可能被延迟到错误的位置。我一般会用线性相位特性的FIR滤波器，或者对所有参与方向判据的通道使用完全相同的滤波器并做群延迟补偿，保证各出线之间的相位关系不被滤波器破坏。滤波器阶数别追求太高，64阶FIR在4kHz采样率下已经足够，阶数过高会带来明显的计算延迟，密集出线的配电终端装不下这么大计算量。
4.2 稳态幅值比较法与相位判据的代码实现
中性点不接地系统中，稳态判据的核心是找出零序电流最大且方向与其他线路相反的那条出线。实现时先计算故障后稳态段的零序电流有效值和相位，再对多条出线做比较。相位要以零序电压为基准，因为各出线零序电流的相位差本身就是相对量，没有统一基准就无法得出“方向相反”这个结论。代码里可以按零序电压的过零时刻建立相位参考系，把零序电流投影到参考系上，用实部符号表示方向。
def steady_state_selection(i0_matrix, u0, sample_rate, fault_time, steady_start=0.04, steady_len=0.02): """ 稳态零序电流幅值与相位联合判据 i0_matrix: 二维数组，shape=(line_num, sample_len)，各出线零序电流 u0: 一维数组，母线零序电压 sample_rate: 采样率，单位Hz fault_time: 故障起始时刻，单位s """ start_idx = int((fault_time + steady_start) * sample_rate) end_idx = int(start_idx + steady_len * sample_rate) # 取故障后稳态数据窗 seg_u0 = u0[start_idx:end_idx] seg_i0 = i0_matrix[:, start_idx:end_idx] # 以零序电压为基准，用DFT提取基波相量 import numpy as np ph_u0 = np.exp(1j * np.angle(np.sum(seg_u0 * np.exp(-1j * 2 * np.pi * 50 * np.arange(len(seg_u0)) / sample_rate)))) ph_i0_list = [] for idx, seg in enumerate(seg_i0): ph_i = np.sum(seg * np.exp(-1j * 2 * np.pi * 50 * np.arange(len(seg)) / sample_rate)) # 归算到零序电压基准下 ph_i0_list.append(ph_i * np.conj(ph_u0) / (len(seg) / sample_rate)) # 幅值最大且实部为负的线路为故障线路 mag_list = [abs(val) for val in ph_i0_list] real_list = [val.real for val in ph_i0_list] max_idx = int(np.argmax(mag_list)) if real_list[max_idx] < 0: return max_idx, {'magnitudes': mag_list, 'real_parts': real_list} return -1, {'magnitudes': mag_list, 'real_parts': real_list} # -1 表示母线故障这段代码里用DFT提取基波相量并投影到零序电压基准方向，ph_u0是零序电压的归一化相量，ph_i0_list中各元素乘上np.conj(ph_u0)等效于把电流相量旋转到电压基准坐标系下。故障线路特征表现为幅值最大、实部为负（方向相反），与理论分析一致。steady_start=0.04表示故障后40ms开始取数据窗，避开故障瞬间的暂态非周期分量，steady_len=0.02对应一个周波的稳态数据，足以稳定提取基波幅值。实际使用中如果发现长线路频繁误判，优先看real_parts列表里是否有两条线路同时接近负向最大，这种情况需要切换到暂态判据或缩小参与比较的线路范围。
4.3 暂态首半波法的实现与灵敏度对比
暂态首半波法选取故障后第一个周波内的暂态电流进行方向判别，其核心是用相位电压参考信号对齐各线路的暂态起始时刻，再比较各线路暂态电流初始极性和峰值大小。故障线路的暂态零序电流极性与非故障线路相反，且幅值通常最大。首半波法的天然优势是不受消弧线圈过补偿影响，劣势是依赖故障初相角，合闸角接近零时暂态能量不足，容易漏判。
def first_half_wave_selection(i0_matrix, u0, sample_rate, fault_time): """ 暂态首半波方向判据 返回故障线路索引，-1表示母线故障或无法判定 """ start_idx = int(fault_time * sample_rate) half_wave_len = int(0.01 * sample_rate) # 10ms，对应50Hz工频半波 # 首半波积分极性：电流在时间轴上的有向面积 area_list = [] for idx, seg in enumerate(i0_matrix): seg_window = seg[start_idx:start_idx + half_wave_len] # 积分前先做零均值处理，去除直流偏置 seg_window = seg_window - np.mean(seg_window) area = np.sum(seg_window) / sample_rate area_list.append(area) # 故障线路首半波积分面积最大且与其他线路符号相反 abs_areas = [abs(val) for val in area_list] max_idx = int(np.argmax(abs_areas)) # 判断是否有其他线路与最大者符号相反 sign_max = np.sign(area_list[max_idx]) opposite_exists = any(np.sign(val) == -sign_max and abs(val) > 0.2 * abs_areas[max_idx] for idx, val in enumerate(area_list) if idx != max_idx) if opposite_exists: return max_idx, {'areas': area_list} return -1, {'areas': area_list}这段代码以故障发生时刻后10ms为数据窗，计算零序电流有向面积并比较极性。opposite_exists的判断条件里加了一个0.2倍幅值阈值，避免某条非故障线路的微弱暂态残留被误认为方向相反的判据。首半波法的关键参数是数据窗长度：窗太短会截断暂态主频分量，窗太长则容易混入稳态工频成分。完整周波内的暂态分量通常集中在故障后3到5ms内，取10ms是为了兼容不同暂态主频的线路，如果发现多次误判，可以把这个窗缩短到5ms再对比结果。
5. 选线装置现场投运避坑指南：5个必踩的现场问题
5.1 问题一：选线结果指向母线，短时复归后又再次动作
现象是装置报“母线接地”，但拉开母线上所有出线后零序电压仍然存在，后台显示接地信号反复动作、复归、再动作。原因多半是间歇性电弧接地，电弧在故障相电压过零附近熄灭、重新燃弧，每次燃弧都会产生新的暂态过程，但故障点始终在母线上。若故障发生在TV一次侧或母线避雷器，也是一个偏门但真实存在的情况，选线装置自身逻辑正确，但排查方向被错误引向出线侧。解决方法是区分“母线TV故障”和“用户侧间歇接地”：查看故障时段各出线零序电流是否同时为零，如果所有出线都没有零序电流而零序电压有规律波动，重点检查TV二次侧和母线避雷器状态。
5.2 问题二：相同型号电缆线路之间反复误选
同一批电缆、同样的长度和截面，两条出线的对地电容电流几乎相等，幅值比较法的裕度不足，噪声稍有波动就会翻转选线结果。解决方法是引入相位信息和暂态量辅助判据，不要只依赖幅值。电缆线路的暂态零序电流主频比架空线高，首半波特征更明显，故障线路在暂态窗内的积分面积通常比非故障线路大一个数量级。工程现场我一般会把幅值比较结果和首半波结果做“与”逻辑：只有当两条判据指向同一条线路时才出口，否则报告“无法确定，待人工复核”。
5.3 问题三：雷雨季节装置频繁启动，却找不到对应故障
雷击引起的感应过电压和绝缘子闪络会产生短暂的高频零序电流脉冲，选线装置灵敏度设得过高时会把这些脉冲当成接地故障启动。原因不是装置坏了，而是启动判据里缺少“持续确认”机制。解决方法是把启动条件从瞬时越限改为“零序电压越限持续50ms以上且零序电流幅值稳定超过阈值”，同时增加三相电压不平衡闭锁逻辑。雷击脉冲持续时间短、幅值衰减快，用持续确认机制可以有效过滤。若故障是绝缘子闪络，第一次脉冲可能已经造成绝缘损伤但未完全接地，后续雨水渗入后才发展为稳定接地，此时录波数据里能看到明显的两次暂态过程，建议保留完整录波文件以便事后溯源。
5.4 问题四：消弧线圈投入后，方向判据全部反转
消弧线圈投入并整定到过补偿状态后，原本可靠的“零序电流最大+方向相反”判据开始把故障线路判断为母线故障。原因前面已经分析过：过补偿状态下故障线路的零序电流中感性分量抵消电容性分量，剩余电流的相位与非故障线路趋于一致。解决方法是针对谐振接地系统改用暂态首半波方向判据，或者用零序有功分量法在消弧线圈的阻尼电阻上构造有功通路。现场的后悔药是：在消弧线圈参数变更后，必须重新做一次短路接地试验验证选线方向，不能直接沿用旧定值。
5.5 问题五：CT极性接反导致所有相位判据失灵
新建间隔送电后选线装置频频误选，相位数据显示多条线路的零序电流方向同时“相反”。这是典型的CT极性接反问题。零序CT一般规定一次侧P1朝向线路、P2朝向母线，二次侧同名端接反后电流相位翻转180度。解决方法是送电前用一次通流试验验证极性，或者利用正常运行时各出线零序电流的相位分布做宏观校验：无故障时各线路零序电流由不平衡电压产生，相位应大体一致，若某一条线路明显反向，优先检查极性。
一条可以抄进作业的排查口诀：幅值最大的不一定是故障线路，方向颠倒的才是。现场排查时先做零序电压相位基准验证，再逐条核对零序CT极性，之后才允许看选线结果，顺序反了，查一晚上也找不到问题。
6. 把选线结果往前推一步：故障录波数据驱动的零序电流突变能量检测
首半波法和稳态比较法解决的是“哪条线故障”，而现场真正难处理的是“故障特征处于两种方法的盲区”——高阻接地且初相角接近零度。此时暂态能量弱、稳态幅值小，两条判据可能同时失效。一个实际有效的技巧是做零序电流突变能量检测：对故障前一个周波的电流波形做滑动窗口能量计算，故障瞬间能量突变最大的线路作为候选故障线路。
def transient_energy_detection(i0_matrix, sample_rate, fault_time): """ 计算每路零序电流在故障前后突变能量，能量突变最大者判为故障线路 用短时能量窗口，窗口长度取2ms，滑动步长0.5ms """ window_len = int(0.002 * sample_rate) step_len = int(0.0005 * sample_rate) pre_len = int(0.02 * sample_rate) # 故障前20ms results = [] for idx, seg in enumerate(i0_matrix): energy_list = [] start_pre = int(fault_time * sample_rate) - pre_len start_post = int(fault_time * sample_rate) for win_start in range(start_pre, start_post, step_len): window = seg[win_start:win_start + window_len] energy_list.append(np.sum(window ** 2) / sample_rate) # 突变能量 = 故障后首个窗口能量 / 故障前窗口平均能量 pre_energy = np.mean(energy_list[:-10]) if len(energy_list) > 10 else np.mean(energy_list) post_energy = energy_list[-3] if len(energy_list) >= 3 else energy_list[-1] mutation_ratio = post_energy / max(pre_energy, 1e-12) results.append(mutation_ratio) # 突变能量比最大的线路判为故障线路 candidate = int(np.argmax(results)) return candidate, results突变能量比阈值的经验值一般取20倍以上，小于这个值说明暂态过程不明显，判据不可靠。配合录波分析时，我会把能量突变曲线和零序电压曲线叠在同一个时间轴上看，确认突变时刻与零序电压起动时刻是否一致，如果不一致说明算法里用的故障时刻有偏差，需要返回到启动判据里修正时标。这个方法不能单独作为选线出口判据，但用来缩小范围、验证前两种方法的结论，效果非常直接。
这个技巧本质上是给选线结果做“第二重保险”。我习惯在每次故障处置完成后，把装置的动作报告、录波数据和最终人工排查结果归档成一个带时间戳的数据包，三套数据对照复盘。坚持半年，你会发现自己对选线装置的信任程度和使用水平一起涨上来。希望帮到你，也让下次那个黑匣子不再给你添堵。
本文还有配套的精品资源，点击获取
