配电网行波测距实战：从Comtrade解析到小波波头检测的完整代码实现

1. 配电网行波测距技术基础行波测距技术是电力系统故障定位的重要手段它利用故障发生时产生的暂态行波信号特性进行定位。当输电线路发生故障时会在故障点产生向线路两端传播的行波信号通过检测行波到达测量点的时间差结合线路参数就能计算出故障距离。与传统阻抗法相比行波测距具有以下优势定位精度高通常误差在300米以内不受系统运行方式影响能适应高阻接地故障响应速度快毫秒级在实际工程中完整的行波测距系统包含以下几个关键环节故障录波数据采集Comtrade格式信号预处理与暂态提取行波合成前向波与反向波波头精确检测小波变换故障类型判断距离计算与结果可视化2. Comtrade数据解析实战2.1 Comtrade格式解析Comtrade是电力系统标准的故障录波格式包含.cfg配置文件和.dat数据文件。我们首先需要解析这些文件获取电压电流信号class ComtradeConfig: 解析Comtrade配置文件.cfg def __init__(self, cfg_path): if not op.exists(cfg_path): raise FileNotFoundError(f配置文件不存在: {cfg_path}) with open(cfg_path, r) as f: lines [line.strip() for line in f.readlines() if line.strip()] # 解析文件基本信息 self.file_info FileInfo(lines[0]) self.channel_info ChannelInfo(lines[1]) # 解析模拟通道 self.analog_channels [] for i in range(self.channel_info.analog): self.analog_channels.append(AnalogInfo(lines[2 i])) # 解析数字通道 start_idx 2 self.channel_info.analog self.digital_channels [] for i in range(self.channel_info.digital): self.digital_channels.append(DigitalInfo(lines[start_idx i]))2.2 数据读取优化技巧在实际工程中我们经常遇到大容量录波文件超过1GB。这时直接读取整个文件会消耗大量内存推荐使用分块读取策略def read_large_comtrade(file_path, chunk_size100000): 分块读取大容量Comtrade文件 parser ComtradeParser(file_path) total_samples parser.config.sample_info[0].end for i in range(0, total_samples, chunk_size): end_idx min(i chunk_size, total_samples) chunk_data parser.data.raw_data[i*parser.data.sample_size : end_idx*parser.data.sample_size] yield process_chunk(chunk_data) # 处理数据块3. 信号预处理与行波合成3.1 克拉克变换实现克拉克变换将三相量转换为模量是行波分析的基础def clarke_transform(three_phase_data): 克拉克变换三相量转模量 clarke_matrix np.array([ [1/3, 1/3, 1/3], # 0模 [1, -1/2, -1/2], # 1模 [0, np.sqrt(3)/2, -np.sqrt(3)/2] # 2模 ]) return np.matmul(clarke_matrix, three_phase_data)3.2 行波合成算法根据电压电流模量计算前行波和反行波def synthesize_traveling_waves(voltage_modal, current_modal, z0, z1): 合成行波0模和1模 uf np.zeros_like(voltage_modal) ur np.zeros_like(voltage_modal) # 0模行波 uf[0] (voltage_modal[0] current_modal[0] * z0) / 2 ur[0] (voltage_modal[0] - current_modal[0] * z0) / 2 # 1模行波波阻抗较小幅值更明显 uf[1] (voltage_modal[1] current_modal[1] * z1) / 2 ur[1] (voltage_modal[1] - current_modal[1] * z1) / 2 return uf, ur4. 小波变换波头检测4.1 小波去噪处理实际信号常含有噪声需要先进行小波去噪def wavelet_denoise(signal, waveletdb8, level3): 基于软阈值的小波去噪 coeffs pywt.wavedec(signal, wavelet, levellevel) # 计算阈值基于Donoho-Johnstone方法 sigma np.median(np.abs(coeffs[-1])) / 0.6745 threshold sigma * np.sqrt(2 * np.log(len(signal))) # 应用软阈值 coeffs[1:] [pywt.threshold(c, threshold, modesoft) for c in coeffs[1:]] return pywt.waverec(coeffs, wavelet)4.2 波头检测算法通过小波分解后的细节系数检测波头位置def detect_wavefront(detail_coeffs, threshold_ratio0.6): 基于细节系数的波头检测 abs_coeffs np.abs(detail_coeffs) max_val np.max(abs_coeffs) if max_val 0: return None threshold max_val * threshold_ratio for i in range(len(abs_coeffs)): if abs_coeffs[i] threshold: # 检查是否为局部最大值 window abs_coeffs[max(0,i-3):min(len(abs_coeffs),i4)] if abs_coeffs[i] np.max(window): return i return None5. 完整工程实现5.1 故障定位主流程def fault_location_analysis(): # 1. 加载两端数据 data_m load_measurement_data(M_CFG_PATH) data_n load_measurement_data(N_CFG_PATH) # 2. 提取暂态信号 transient_voltage_m [extract_transient_signal(data_m[voltage][i], data_m[fault_time], data_m[fs]) for i in range(3)] # 3. 克拉克变换 voltage_modal_m clarke_transform(np.array(transient_voltage_m)) # 4. 行波合成 uf_m, ur_m synthesize_traveling_waves(voltage_modal_m, ...) # 5. 小波波头检测 coeffs wavelet_decompose(uf_m[1]) detail extract_detail_coeffs(coeffs, 3) wavefront_idx detect_wavefront(detail) # 6. 距离计算 distance distance_measurement(t_m, t_n, line_params)5.2 结果可视化技巧使用Matplotlib创建专业级可视化报告def create_visual_report(): fig plt.figure(figsize(10, 8)) # 1. 行波信号对比 ax1 fig.add_subplot(311) ax1.plot(t_m, uf_m[1], labelM端前行波) ax1.plot(t_n, uf_n[1], labelN端前行波) # 2. 小波细节系数 ax2 fig.add_subplot(312) ax2.plot(detail_m, labelM端细节系数) # 3. 故障位置示意图 ax3 fig.add_subplot(313) ax3.plot([0, line_length], [0, 0], b-, linewidth3) ax3.plot([fault_pos, fault_pos], [-1, 1], r--)6. 工程实践中的注意事项在实际项目中有几个关键点需要特别注意波阻抗选择不同电压等级线路的波阻抗差异较大500kV线路典型值为260Ω220kV线路为380Ω。错误的值会导致行波幅值计算偏差。采样率要求行波测距需要高采样率通常≥1MHz低采样率会导致波头检测误差。如果原始数据采样率不足可以考虑插值处理。接地方式影响小电流接地系统的行波信号较弱需要更高的检测灵敏度。这时可以适当降低波头检测阈值。线路参数校正实际波速受线路参数影响建议定期进行波速实测校正。可以通过已知距离的人工短路试验来校准。多端数据同步对于T接线路等复杂拓扑需要确保各端数据时间同步必要时使用GPS对时信号。我在某220kV线路项目中曾遇到波头检测不稳定的问题后来发现是线路参数中的波速取值不当。通过现场实测将波速从光速的98%调整为95%后定位精度从500米提升到了150米以内。