输电线路行波测距原理与Simulink仿真复现全流程解析

搞电力系统故障分析的人，十有八九都绕不开“行波测距”这四个字。最近我把一篇题目为《基于行波理论的输电线路故障诊断方法研究》的论文在 Simulink 里完整复现了一遍。整个流程走下来，从双端电源系统建模、故障模块配置，到小波变换提取行波波头、双端测距公式验证，中间踩了不少坑，也把论文里的关键图表一一对应上了。如果你正在做毕业设计、论文复现，或者想从零了解行波测距在仿真里到底怎么落地，这篇内容应该能帮你把主线理清楚。

先说一下这个内容解决什么问题。输电线路是电网里故障率最高的设备，故障后能不能快速、准确地找到故障点，直接影响检修效率和供电恢复。行波法不像传统阻抗法那样依赖工频电压电流的稳态关系，而是直接检测故障瞬间产生的行波波头到达测量点的时间差，原理上几乎不受过渡电阻、负荷电流和系统两侧电源结构的影响。因此它的测距精度通常高于传统方法。对研究生来说，这既是一个经典的论文方向，也是 Simulink 暂态仿真、信号处理、小波分析等多种技能的集中训练场。

下面我按实际复现顺序来写，先讲原理和模型选型，再拆解 Simulink 搭建细节，然后给出小波波头检测和测距算法的完整实现，最后分享常见问题排查思路。每个环节我都尽量把参数和单位写清楚，免得你照着做的时候还要自己去猜。

1. 行波故障诊断的核心逻辑与仿真选型

1.1 为什么故障测距要选行波法

传统阻抗测距的基本思路，是利用故障后保护安装处感受到的电压、电流计算出故障回路阻抗，再用单位长度阻抗换算距离。听起来很直接，但实际工况里问题很多。过渡电阻会显著抬高测量阻抗中的电阻分量，让测距结果偏大；双侧电源线路中，对侧电源提供的短路电流会影响电压电流相位关系，导致阻抗角失真；负荷电流、线路分布电容也会介入。误差百分之几到十几都很常见。

行波法完全换了一个思路。故障发生的一瞬间，故障点等效于突然施加了一个阶跃电压源，会产生向线路两端传播的暂态行波。这个行波以接近光速的速度沿输电线路传播，只要我们能准确捕捉到它到达线路两端母线的时间，就能用非常简单的时间差公式计算出故障距离。整个过程不依赖故障后的稳态分量，过渡电阻只影响行波幅值的大小，几乎不影响到达时刻，所以高阻接地故障场景下优势格外明显。

打个生活比方：你站在长走廊一头敲一下水管，水管里的冲击波会传到另一头，根据冲击波到达的时间和波速，就能估算出敲击点离你多远。行波测距就是这个思路的电力系统版本，只不过把“声音”换成电磁波，“走廊”换成输电线路。

1.2 行波的产生、反射与传播规律

故障点产生的初始行波包含两个方向：一个向 M 端母线传播，一个向 N 端母线传播。行波沿线传播时，如果遇到波阻抗不连续的地方，比如母线、故障点、线路分支、电缆与架空线连接处，就会发生反射和透射。对于双端测距来说，我们只关心两端测量点第一次接收到初始行波波头的时刻，不需要去分辨那些复杂反射波。这也是双端法比单端法更容易实现的原因。

双端测距公式是：

x = [L + v × (tM - tN)] / 2

其中 x 是故障点距 M 端的距离，L 是线路全长，v 是行波速度，tM 和 tN 分别是初始行波波头到达 M 端和 N 端的绝对时间。推导也很简单：设故障发生时刻为 t0，则 tM = t0 + x/v，tN = t0 + (L-x)/v，两式相减消掉 t0 就能得到上面的公式。公式里不需要知道故障发生时刻，这是双端法的最大优势。

如果是单端测距，则需要利用故障点反射波与初始波头的时间差，公式是 x = v × Δt / 2。单端法只需要一端安装测量装置，更省钱，但波头识别难度高很多，因为需要区分故障点反射波和母线反射波。很多论文会先做双端，再延伸到单端；我的复现也以双端为主。

1.3 为什么必须用分布参数线路模型

Simulink 里输电线路有两种常用模型：PI 型集中参数线路和分布式参数线路。很多人第一次做行波仿真时不注意，直接拖了一个 PI Section Line 模块，结果发现测距怎么都不对。原因很简单：PI 模型把线路的电阻、电感、电容集中到几个元件里，纯粹为了工频稳态和短路电流计算设计，它无法体现行波沿线路逐米传播的时间延迟。

行波测距的物理本质就是时间测量，所以必须使用 Distributed Parameters Line 模块。这个模块内部采用特征线法或其他行波传播模型求解，能够真实反映行波在沿线各个位置的时延和反射透射过程。可以这么理解：PI 模型像一根一次性把水灌满的水管，分布式参数线路则像沿整根水管连续分布的水锤波，你敲一下末端，波头要花真实的时间传到另一端。选错模型，后面所有工作都白做。

2. Simulink 模型搭建与核心模块配置

2.1 系统拓扑与模块选型

我的复现模型采用双端 220kV 输电系统结构：M 端和 N 端各有一个等效电源，中间通过一条 100km 架空线路相连，故障点放在线上指定位置。两侧电源都带等效内阻抗，模拟实际电网的戴维南等值。线路两端分别安装三相电压电流测量模块，把行波信号导出到 MATLAB 工作空间做后续处理。

具体用到的主要模块如下：

电源模块的线电压我设为 220kV，频率 50Hz，内阻抗按短路容量等值出来的一组小电阻电感即可。重点是线路模块必须选 Distributed Parameters Line，如果界面上看到 PI Section Line 选项，不要用。

2.2 线路参数与故障模块设置

分布参数线路需要填入正序和零序的电阻、电感、电容每公里参数。下面这组参数是我从一篇常见 220kV 架空线论文里取的典型值，仿真完全够用：

这组参数算出来的线模波速大概是：

v = 1 / sqrt(L1 × C1) = 1 / sqrt(0.9e-3 × 12.7e-9) ≈ 2.96e8 m/s

也就是约 0.296 km/μs。100km 线路单程传播时间约 338μs，所以在故障后很短的时间窗内就能看到行波到达两端。记得把线路长度填 100km，故障点位置可以放在 20km、40km、80km 等不同位置来反复验证。

故障模块我用 Three-Phase Fault。最常见的验证场景是 A 相单相接地，故障类型勾选 A 相接地；也可以做 AB 两相接地、ABC 三相短路。过渡电阻设为 10Ω，故障时间写为 [0.025 0.05]，也就是 0.025s 时刻发生故障，0.05s 切除。选 0.025s 是因为在这个时刻附近 A 相电压接近峰值，故障产生的行波幅值比较大，波头容易观察。如果你用的电源相位不是从 0 开始，就先把模型跑一遍看故障相电压波形，把故障触发点设在电压峰值附近。行波法虽然对过渡电阻不敏感，但如果故障发生在电压过零点附近，行波幅值会非常小，给后续波头识别增加难度。

2.3 离散仿真环境与信号采集

行波仿真最关键的一个设置是仿真模式。我强烈建议在 Powergui 里选择 Discrete 离散化模式，采样时间设为 1e-6 秒，也就是 1MHz 采样率。Simulink 求解器也要选固定步长，步长与 Powergui 的采样时间保持一致，否则会出现波头被平滑或者仿真极慢的问题。

1μs 的采样分辨率对应多少测距误差呢？行波波速 0.296km/μs，也就是说 1μs 的时间量化大约对应 296m 的距离误差。双端公式里时间差除以 2，所以最终测距误差通常在几百米量级。如果你希望精度更高，可以把采样时间压缩到 2e-7 秒，但仿真耗时会明显增加。

测量模块用 Three-Phase V-I Measurement 分别挂在 M 端和 N端，紧接着接 To Workspace 导出电流信号。我一般导出六路信号：iAM、iBM、iCM、iAN、iBN、iCN。电压也可以同时导出，方便后面用电压行波做交叉验证，但只需要线路电流时没必要把全部节点电压都导出，否则数据量会很大。To Workspace 的 Save format 我习惯用 Timeseries，便于在脚本里直接读取时间和数据。

3. 行波波头检测与测距算法实现

3.1 相模变换：不要把三相电流直接拿来用

三相输电线路各相之间存在电磁耦合，A、B、C 相电流不是相互独立的。如果直接把某一相电流当作单相行波来处理，会混入其他相的耦合分量，波头位置可能会偏移。工程和论文里通常都会先做相模变换，把三相电流变换成线模分量和零模分量。

零模分量通过大地和避雷线回路传播，波速受大地电阻率、土壤参数影响很大，而且频变特性明显，不适宜用来测距。线模分量基本沿架空导线回路传播，波速稳定，接近光速，是测距的首选。所以实际操作中，我从三相电流里构造线模分量，再送到小波变换模块里检测波头。

简单实现可以这样取：

% ia, ib, ic 是从 Simulink 导出的三相电流数组 i_line = (ia - ib) / sqrt(2); % 近似线模分量，具体系数以你的变换矩阵为准

如果你的论文里有完整的 Karenbauer 变换矩阵，直接用矩阵乘法更好。但不管用哪种变换，核心原则是一样的：用线模分量测距，不要用零模分量。零模波速慢，如果错误地把零模波头当成线模波头，测距结果会偏大很多。

3.2 小波变换模极大值定位波头

行波波头本质上是信号里的一个突变成分，落在高频段。工频分量能量很低频，用高速滤波器就能把突变分离出来，但简单滤波器会带来相位延迟，定位精度不够。小波变换的好处是能够在不同尺度下分析信号，突变点会在细节系数上表现出模极大值，而且定位能力比普通滤波好得多。

我使用的是 db4 小波，做 6 层多分辨率分解。Matlab 脚本大致如下：

fs = 1e6; % 采样率 1MHz，对应采样间隔 1us t = out.tout; % 截取故障前 0.0005s 到故障后 0.003s 的区间，减少无关数据量 idx = t >= 0.0245 & t <= 0.028; i_line = i_line(idx); t_x = t(idx); % 多分辨率分解和细节信号重构 [C, L] = wavedec(i_line, 6, 'db4'); D1 = wrcoef('d', C, L, 'db4', 1); % 第1层细节，最高频 D2 = wrcoef('d', C, L, 'db4', 2); % 第2层细节 % 寻找第二个尺度细节系数的模极大值位置 absD = abs(D2); th = 0.2 * max(absD); pos = find(absD > th, 1, 'first'); t_arrive = t_x(pos);

wrcoef 重构出来的细节信号长度和原信号一致，索引可以直接对应到时间轴，比单层 dwt 输出需要换算索引更省事。阈值取 0.2 倍最大值是我试出来比较稳定的经验值，因为初始行波波头通常是整段数据里幅值最大的突变。如果你的模型里存在明显干扰，阈值可以适当提高。

需要留意的是，D1 层级太高频，容易把数值振荡和开关噪声当成波头；D3 或 D4 虽然更平稳，但定位时间会有几个采样点的偏移。我通常以 D2 或 D3 为准，两个尺度都看一下，如果结果一致，就说明波头定位是可靠的。

3.3 双端测距公式与时间同步处理

从 M、N 两端分别得到初始波头到达时刻后，代入双端公式即可：

L = 100; % 线路长度 km v = 0.296; % 波速 km/us tM_us = tM_arrive * 1e6; % 把仿真时间换算成微秒 tN_us = tN_arrive * 1e6; x_cal = (L + v * (tM_us - tN_us)) / 2;

注意单位统一。tM 和 tN 是绝对到达时刻，不是时间差。由于故障点更靠近 M 端时，M 端先收到波头，tM - tN 会是负数，公式里 L 加上负数再除以 2，依然能算出正确距离。我一开始看到这个负号还怀疑公式写错了，其实把故障距离代进去验证一次就知道没问题。

真实系统里双端测距必须解决两端时钟同步问题。仿真模型天然共享同一个仿真时钟，所以不存在对时误差。但如果你在论文里讨论工程实现，一定要提双端对时装置，否则审稿人很容易揪住这个问题。

3.4 波速校验技巧

行波测距的误差很大程度上取决于波速取得到底准不准。理论波速可以直接用 L1 和 C1 算，但仿真里线路模型、导线参数、大地回路都会对实际传播速度产生影响。更稳妥的做法是用一个已知故障距离来反推波速。

方法很简单：先在线路某个已知位置设一个故障，比如 x = 40km，记录两端波头到达时间 tM 和 tN，然后用公式：

v_cal = abs((2x - L) / (tM - tN))

反推出当前模型里的实际波速。我反推出来的值通常在 0.294 到 0.298 km/μs 之间，和理论值很接近。把这个反推波速再放回批量测距脚本里，比直接拿理论光速去算要稳得多。这个“用已知答案标定参数”的做法在复现论文时非常实用。

4. 仿真结果分析与论文图表复现

4.1 典型故障场景的测距结果

仿真模型跑通后，我做了三组典型故障场景：A 相单相接地、AB 两相短路接地、ABC 三相短路。故障距离分别取 20km、40km、80km，采样率 1MHz，过渡电阻 10Ω。一次完整仿真时长 0.05s，结果是稳定的。

从示波器波形能明显看到：当故障发生在 20km 处，M 端电流行波在故障后约 68μs 出现第一个尖峰，N端则在约 270μs 后出现。这个时间差代入双端公式正好得到接近 20km 的结果。小波 D2 系数在对应位置也会出现一个非常突出的模极大值，位置和小波包络的视觉观测完全一致。

我其中一组复现结果如下：

绝对值都在 0.5km 以内，满足论文里常见的测距精度指标。剩下的几十米误差主要来自 1μs 采样量化和波头定位偏差，属于正常现象。如果采用 0.1μs 采样率或者对模极大值做插值，误差还能进一步缩小，但仿真时间会成倍增加。

4.2 误差来源与过渡电阻影响

复现过程中，我特意分析了过渡电阻对测距结果的影响。分别在 10Ω、100Ω、300Ω 过渡电阻下做了 A 相接地故障测试，结果发现：过渡电阻增大后，行波波头幅值明显变小，但到达时刻几乎没有变化，测距误差始终在 0.5km 以内。这说明行波法在高阻接地场景下确实比阻抗法稳定得多，这也正是论文里反复强调的核心优势。

影响测距误差的主要因素有三个：

第一个是采样率。1μs 采样率带来约 0.3km 的量化误差，提高采样率是最直接的改善手段。第二个是波速取值。线路参数会随频率变化，如果取理论光速而没有用实际模型参数校准，误差可能达到几公里。第三个是小波波头定位的边界效应。靠近数据截断位置时，小波系数边缘会出现异常，所以截取时间窗时要在故障点前后留足余量。

母线反射波的干扰也会影响定位。双端测距只认第一个波头，如果程序不小心把第二个波头识别成初始波头，测距结果会偏移好几个公里。这也是为什么我在自动检测脚本里严格要求只取第一个模极大值，而且要查看两端的检测结果是否与故障距离相匹配。

4.3 批量仿真与误差曲线绘制

论文里通常都会给一张“不同故障距离下的测距误差”曲线，这种图需要批量跑仿真才能画出来。我的做法是把故障距离设成变量，在 MATLAB 脚本里循环修改模型参数，每次仿真后自动计算测距结果并保存。

distances = [5 10 20 30 40 60 80 90]; results = []; for x = distances % 修改模型中故障位置的参数，具体信号名以实际模型为准 % set_param('line_model/Fault', 'Distance', num2str(x)); out = sim('line_model', 'StopTime', '0.05'); % 从 out 中提取 M、N端电流波头到达时刻 % tM_arrive = detect_wavehead(out, 'M'); % tN_arrive = detect_wavehead(out, 'N'); v = 0.296; x_cal = (L + v * (tM_arrive - tN_arrive) * 1e6) / 2; results = [results; x, x_cal, x_cal - x]; end figure; plot(distances, results(:,3), 'o-'); xlabel('实际故障距离 (km)'); ylabel('测距误差 (km)');

这个脚本跑完后，得到的误差曲线非常平滑，方便放在论文结果分析部分。需要注意的是，Simulink 模型里通过 set_param 修改参数时，参数名和路径要格外小心。我更推荐把关键参数做成 MATLAB 工作区变量，在模型对话框里直接引用变量名，这样循环仿真时只需要给变量重新赋值，代码会干净很多。

5. 常见问题与排查技巧实录

5.1 波头找不到或者时有时无

这是复现行波测距时最多人遇到的问题。我第一版模型就栽在这里，明明故障加了，电流波形也有变化，但小波模极大值的位置每次都不一样。后来逐个排查才发现，是故障时刻刚好落在电压过零点附近，行波幅值小到几乎被工频分量淹没。

解决办法是把故障时刻调整到故障相电压峰值附近。电源是 50Hz，电压周期 20ms，峰值出现在 5ms、15ms、25ms 这些位置。如果你的电源初始相位不是 0，先看一眼波形再设故障时刻。另外，不要只取第 1 层小波细节，D1 包含太多高频噪声，把定位结果搞得忽前忽后。我建议至少看 D2 和 D3 两个尺度，以稳定出现的极大值位置为准。

还有一个经验是截取时间窗。如果从仿真 0 时刻开始截取，故障前几十毫秒的数据里会有很多工频和启动暂态，小波最大极值往往落在故障前区域。正确做法是只关注故障前 0.5ms 到故障后 2ms 的时间窗口，这样波头就是窗口内最突出的突变。

5.2 测距误差偏大时的排查方向

如果仿真跑出来的测距结果与实际故障位置差得很远，别急着怀疑算法，先按下面这个思路逐项排查：

最常见的问题还是波速。我见过有人直接用光速 3e8 m/s 去算，误差一下就到 2km 以上。架空线的线模波速确实接近光速，但通常只有 0.29 到 0.30 km/μs 左右，1%~2% 的偏差在线路全长上就会放大成可观的测距误差。所以务必用你模型里的线路参数重新计算波速。

5.3 仿真慢、不收敛、数值振荡怎么办

用连续求解器跑带分布参数线路的行波仿真，可能会慢到怀疑人生。尤其是步长默认通过误差控制自动调整，遇到暂态过程会疯狂加密计算，一个 0.05s 的仿真可能跑十几分钟。解决方法是把 Powergui 设成 Discrete 模式，求解器用固定步长的 discrete 求解器，步长与 Powergui 的采样时间一致。这样仿真速度通常会快一个数量级以上。

分布参数线路在离散仿真里有时会出现高频数值振荡，表现为波头上叠加很多毛刺，影响波头定位。这类振荡常见原因是电路里的电感电容和求解器步长相互作用。处理方法是在母线测量点附近适当加入 snubber 支路，或者在 fault 模块里调节缓冲电阻参数。多数情况下，只要 Powergui 采样时间不小于仿真步长，且步长足够小，就不会出现严重振荡。

数据存储也是一个容易忽略的点。如果把所有节点电压电流都加到 To Workspace，仿真结束时会导出巨大的矩阵，内存占用高，处理也慢。我只保存两端线路电流，最多再加两端各相电压，其他中间量一律不导。仿真时长也不需要做满 0.05s，故障后 2ms 行波过程已经全部结束，如果只是为了验证测距，甚至可以把 StopTime 设成 0.028s。

6. 复现经验总结与后续扩展建议

6.1 复现这类论文的正确顺序

复现行波测距论文，最大的坑不是不会用某个模块，而是没有先建立对结果的预期就直接开跑。我常用的顺序是这样的：先把论文里的公式和参数抄一遍，确认每一个变量的单位和取值范围；然后建一个最简单的单相分布参数线路模型，用一个电压阶跃源验证行波能在线上以正常速度传播；接着把三相系统、故障模块加进去，确认能跑出清晰的故障暂态波形；最后才加小波分析和自动测距脚本。

这个顺序看起来多了一步，实际上最省时间。因为一旦后面测距结果不对，你可以分出来到底是模型物理过程错了，还是信号处理定位错了。如果直接搭一个大模型再排查，变量太多，很难定位问题。

“带着答案复现”也是很重要的一招。论文里如果给出了故障距离和测距结果，你先把理论值算一遍。比如故障在 20km，波速 0.296km/μs，M 端和 N 端的时间差理论值大约是 -203μs。看到这个数量级之后再去跑仿真，只要数量级对不上，马上就能发现参数设置错误，不需要等所有后期处理做完才意识到出了问题。

6.2 从双端测距走向更完整的故障诊断

这次复现完成后，我认为这个模型很容易扩展成更完整的故障诊断方案。你可以把双端测距改成单端测距，核心难点从“识别第一个波头”变成“区分故障点反射波和母线反射波”，这需要补充故障点反射和母线反射的极性判断规则。也可以在线路中间设置功率分支或电缆连接点，研究行波在波阻抗不连续处的透射反射特征，这对混合线路很有意义。

故障类型判别也是自然延伸。零模分量和线模分量的波头幅值、极性关系在不同故障类型下是有规律的，比如接地故障会激发出明显零模分量，而相间短路几乎没有零模。利用这些特征可以做一个简单的故障类型分类器。再往上走，还可以把多层小波能量、模极大值组合等特征喂给分类模型，实现行波特征驱动的智能故障识别。

最后说一个我在实际复现中的体会：行波测距这类论文，最容易让你翻车的往往不是算法本身，而是仿真里的时序和采样口径。我把两个测量端的采样起始时刻统一从 0 开始后，之前百思不得其解的随机误差一下子消失了。这个坑很基础，但真踩过一次之后，你会对“时间同步”这三个字有完全不一样的理解。希望这篇内容能帮你少走一次这样的弯路。