1. 输电线路故障诊断的行波理论基础
行波理论在电力系统故障诊断中的应用可以追溯到上世纪50年代，其核心思想是利用故障发生时产生的暂态行波信号进行故障定位和类型识别。当输电线路发生故障时，故障点会瞬间产生向线路两侧传播的电压和电流行波，这些行波在传播过程中遇到阻抗不连续点（如母线、分支点等）会发生反射和折射现象。
行波传播具有几个重要特性：首先，行波传播速度接近光速（约3×10^8m/s），这使得基于行波的故障定位可以达到很高的时间分辨率；其次，不同频率的行波分量在传播过程中衰减特性不同，高频分量衰减更快；最后，三相线路中的行波可以通过相模变换解耦为独立的模量，这大大简化了分析过程。
在工程实践中，行波法相比传统阻抗法具有明显优势。阻抗法受系统运行方式、线路参数等因素影响较大，而行波法主要依赖行波到达时间差，受系统运行状态影响小。此外，行波法可以实现更高精度的故障定位，理论上误差可以控制在几百米以内。
2. 仿真模型构建与参数设置
2.1 Simulink模型架构设计
基于文献复现的仿真模型采用模块化设计思路，主要包含以下几个关键部分：
-
电源系统模块：模拟220kV三相电源，包含正序、负序和零序阻抗参数设置。电源频率设为50Hz，内阻抗根据典型220kV系统参数配置。
-
输电线路模块：采用分布参数模型，线路总长度200km。线路参数包括：
- 正序电阻：0.021Ω/km
- 正序电感：1.154mH/km
- 正序电容：0.009μF/km
- 零序电阻：0.115Ω/km
- 零序电感：3.462mH/km
- 零序电容：0.006μF/km
-
故障模拟模块：可设置多种故障类型（单相接地、相间短路、三相短路等），故障起始时间、过渡电阻等参数可调。本案例设置AC相间短路故障，故障点距L端90km。
-
测量模块：在线路两端安装电流互感器和电压互感器，采样率设为1MHz以满足行波信号采集需求。
2.2 关键参数选择依据
线路模型选择分布参数而非集中参数，是因为行波分析涉及高频暂态过程，分布参数模型能更准确反映行波传播特性。采样率设置为1MHz是基于以下考虑：
-
行波波头包含丰富的高频分量，根据Nyquist定理，采样率至少应大于最高有效频率的2倍。
-
实际行波采集装置的采样率通常在1-10MHz范围，1MHz是兼顾精度和计算复杂度的折中选择。
-
对于200km线路，行波传播时间约0.67ms，高采样率可确保准确捕捉行波到达时刻。
3. 故障信号处理与特征提取
3.1 相模变换实现故障类型识别
相模变换是将三相系统解耦为独立模量的数学工具，常用的变换包括Karenbauer变换和Clarke变换。本模型采用Karenbauer变换：
code复制[I0 Iα Iβ] = 1/3 * [1 1 1; 2 -1 -1; 0 √3 -√3] * [Ia Ib Ic]
其中I0为零序分量，Iα和Iβ为线模分量。不同故障类型的模量特征如下：
- 单相接地：零序分量显著
- 相间短路：线模分量显著，零序分量小
- 三相短路：仅线模分量
通过监测各模量幅值比例，可准确判断故障类型。本案例中AC相短路时，Iα和Iβ分量明显增大而I0变化很小，符合相间短路特征。
3.2 小波变换用于波头检测
小波变换非常适合分析非平稳的行波信号，其优势在于：
- 多分辨率分析能力：可以同时在时域和频域局部化信号特征
- 对奇异点敏感：能够有效检测行波波头的突变点
本模型选用db4小波进行6层分解，原因如下：
- db4小波具有适中的支撑长度和消失矩，在计算复杂度和时频分辨率间取得平衡
- 6层分解可将信号分解到合适的频带，D1层（最高频）最适合捕捉波头突变
波头检测的具体步骤：
- 对α模量电压信号进行小波分解
- 提取D1层细节系数
- 寻找细节系数的模极大值点
- 通过模极大值极性判断初始波头与反射波头
4. 故障定位算法实现
4.1 双端行波定位原理
双端法利用线路两端检测到的初始行波到达时间差进行定位，计算公式为：
code复制x = (L - v*Δt)/2
其中：
- x：故障点距L端距离
- L：线路总长度（200km）
- v：行波传播速度（约光速的98%）
- Δt：行波到达两端的时间差
本案例中，测得Δt=0.000133s，计算得到x=89.775km，与预设值90km非常接近。
4.2 提高定位精度的措施
-
波头精确标定：采用小波变换模极大值法，相比传统阈值法可减少1-2个采样点的误差。
-
速度校正：考虑频率相关的地模速度衰减，对零模分量采用较低传播速度（约光速的70%）。
-
数据同步：两端数据采用GPS同步时钟，时间同步误差控制在100ns以内。
-
抗干扰处理：在信号预处理阶段加入数字滤波，抑制噪声对波头检测的影响。
5. 仿真结果分析与验证
5.1 典型波形分析
-
正常状态波形：三相电流呈平衡正弦波，幅值相等相位互差120度。
-
故障后波形：
- AC相电流突然增大，相位趋于相同
- B相电流略有变化但不显著
- 行波分量在故障后0.3ms内达到峰值
-
模量分量：
- α分量反映AC相间故障特征
- β分量变化较小
- 零序分量几乎为零，验证了不接地故障特性
5.2 定位误差分析
本案例定位误差仅0.112%，主要误差来源包括：
- 行波速度假设为恒定值，实际受频率影响
- 波头检测存在±1个采样点的理论误差
- 线路参数的不确定性
- 仿真模型的理想化假设
通过多次仿真测试，该方法在不同故障类型和位置下的平均定位误差小于0.2%，满足工程精度要求。
6. 工程应用中的注意事项
-
采样系统要求：
- 采样率不低于500kHz
- ADC分辨率建议12位以上
- 模拟带宽需覆盖10kHz-1MHz
-
现场安装要点：
- 电流互感器应具有宽频带特性
- 信号传输采用屏蔽电缆
- 确保GPS天线良好接收
-
算法实现建议：
- 实时处理采用FPGA实现小波变换
- 保留原始数据用于事后分析
- 建立典型故障波形库辅助诊断
-
系统维护提示：
- 定期校验时间同步精度
- 监测CT/PT的频率响应特性
- 更新线路参数变化
7. 模型扩展与改进方向
-
加入噪声测试：在信号中加入实测噪声，验证算法鲁棒性。
-
考虑弧光电阻：模拟不同过渡电阻下的故障特性。
-
复杂网络扩展：研究T接线路、环网等拓扑下的行波传播特性。
-
AI辅助诊断：结合深度学习算法提高故障类型识别准确率。
-
硬件在环测试：连接实际采集装置进行闭环验证。
这个仿真模型为深入研究行波故障诊断提供了良好平台，通过调整参数和拓扑可以研究各种复杂场景下的行波特性。在实际应用中，还需要考虑通信延迟、数据同步等工程因素，这些都可以在现有模型基础上进一步扩展研究。