简介面向电力系统研究、开发与运维人员聚焦基于电磁时间反转EMTR理论的电网故障定位问题资源将论文复现与工程实践相结合适合具备电磁学、信号处理和编程基础的研发者学习。虽然包体精简但内容并未缩水全部资料集中为1份PDF文档约680KB既包含EMTR故障定位原理与频域传输线方程推导也提供完整的Python代码覆盖传输线模型定义、故障信号生成、EMTR故障定位算法和模拟实验设计。目前已有136人学习下载。文档还讨论了观测点数量对定位精度的影响并通过小型同轴电缆硬件模拟实验与电磁暂态仿真案例验证了方法在不同场景下的有效性针对时间反转不变性、反射系数影响等关键理论做了验证同时提出多观测点优化布置、时频域混合处理等工程应用建议便于读者结合代码逐段调试将算法落地到实际电力系统项目中。1. EMTR 电网故障定位为什么值得自己复现一遍一条 110 kV 双端输电线路发生单相接地后保护动作只是第一步——运维真正需要回答的问题是故障点在几十公里走廊里的哪一段。传统阻抗法受过渡电阻干扰单端行波法在波头识别上有明显盲区。基于电磁时间反转EMTR的电网故障定位方法把故障定位从“找波头时刻”换成了“找能量聚焦点”记录故障暂态波形把波在时间上倒过来重新注入线路能量会在故障位置附近形成一个峰值。这个方法不需要训练样本不依赖精确的阈值标定适合作为现有行波定位装置的算法替代。下面从物理原理、MATLAB 实现到实验验证给出完整可复现的一条路径看完你可以直接拿自己的录波数据跑一轮定位。2. EMTR 定位原理时间反转在输电线路里如何聚焦故障点时间反转最早在声呐和超声成像里被反复验证先用阵列记录水下目标反射的回波再把采样信号做时间反转后重新发射声波会沿原路径返回并在目标位置重新聚焦。这个现象被称为时间反转镜。把同样的逻辑搬到输电线路只需要记住一个关键事实故障点不是普通位置而是线路阻抗的不连续点故障瞬间产生的暂态行波从这里出发分别向线路两端传播。2.1 时间反转为什么能让故障点聚焦这里起作用的是电磁场的互易性。在无损线性近似下从故障点 xf 传播到 A 端的波与从 A 端注入、传到 xf 的波经历的幅度和相位变化完全相同反向传播不过是把路径倒过来走一遍。更直接的是频域视角时间反转等价于在频域对信号取相位共轭。一个正向传播的波相位是-ωx/v反转后相位变成ωx/v这正是沿原路径反向传播的波。两端记录到的信号各自做时间反转再作为激励源重新注入两个波前会在故障点同时到达——注意是“同时”——叠加出的干涉峰就是定位结果。2.2 从双端行波模型到聚焦函数设线路总长 L故障点距 A 端 xf故障暂态源为 s(t)。理想无损线路下两端收到的信号可以写成UA(t) s(t - xf/v) 反向反射分量UB(t) s(t - (L - xf)/v) 反向反射分量先忽略反射只看第一次到达的波前。变换到频域UA(w) S(w) · exp(-j·w·xf/v)UB(w) S(w) · exp(-j·w·(L-xf)/v)对每一个假设的故障位置 x构造聚焦谱G(w, x) UA(w) · conj(HA(w, x)) UB(w) · conj(HB(w, x))其中HA exp(-j·w·x/v)是假设位置到 A 端的传播算子HB exp(-j·w·(L-x)/v)同理。取共轭正是在频域里做时间反转。然后计算聚焦能量E(x) ∫ |G(w, x)|² dw当x xf时两个指数项的相位正好被共轭项抵消G(w, xf) 2·S(w)能量取得极大值。x偏离故障位置时两个指数项之间出现2w(x - xf)/v的相位差积分后能量下降。整条线路上E(x)的峰值位置就是 EMTR 估计的故障距离。这段推导是论文复现的核心。所谓的“代码”在上面的公式里实际上只剩两件事把两端波形变换到频域再对每个候选位置算一组带相位共轭的内积。2.3 EMTR 与传统行波法在定位机制上的差异对比维度传统双端行波法基于 EMTR 的电网故障定位核心输入波头到达时刻整段故障暂态波形对波头阈值依赖高阈值设置直接影响判定结果低积分运算对单点噪声不敏感波速误差影响直接乘进距离公式产生等比例偏差使聚焦峰钝化位置偏移相对较小反射波处理通常作为干扰剔除反射本身参与聚焦可用短窗隔离实现代价一次时间差计算极低需遍历候选位置做频域相关计算量随网格数线性增长但这不是说 EMTR 全面优于行波法。当线路损耗和频变参数明显时比如长距离电缆EMTR 的反演算子必须包含幅度补偿项否则聚焦峰被拉宽。复现论文时最常遇到的坑也在这里解析算法里用恒定波速仿真软件里用的是频率相关线路模型两边的聚焦位置会出现系统性偏移。常见做法是回注步骤直接把两端记录电压作为源注入 PSCAD 或 Simulink 模型让仿真器自带的线路模型承担补偿代价是每个候选点都要完整跑一次电磁暂态仿真。3. 论文复现第一步MATLAB 实现双端 EMTR 定位代码搜“论文复现”时大家默认先找 GitHub 源码如果找不到就开始焦虑这基本是深度学习社区带起来的习惯。EMTR 论文的复现路径完全不同它没有训练、没有权重核心是物理方程和传播算子。代码写得再笨只要 UA、UB 和传播算子对得上结果就不会差。反向的教训是对不上的时候哪怕有源码也照样白跑。所以我的复现习惯是三步先用仿真生成一组“已知答案”的波形把 EMTR 算法写成独立函数最后换参数做交叉验证。3.1 EMTR 复现的通用数据流程获取故障暂态数据真实录波、PSCAD 输出或实验室采集采样率建议不低于 1 MSa/s双端对齐保证 A、B 两端时间起点一致预处理去直流、带通滤波消除工频分量和高频噪声算法处理沿线网格逐点计算聚焦能量输出聚焦能量峰值对应的距离以及一条便于人工复核的聚焦曲线。第 4 步有两种实现方式。解析算子用 FFT 和相关运算速度快适合论文复现和在线部署仿真回注把反转信号作为激励重新塞回电磁暂态模型精度高但慢。下面的代码采用解析算子如果你手里已经有 PSCAD 或 Simulink 的仿真波形直接替换 UA 和 UB 即可。3.2 正向数据生成先制造一个“已知答案”的输入% 模拟一条双端无损架空线故障点产生宽带暂态脉冲 L 100; % 线路长度 km vp 2.98e5; % 波速 km/s约 0.98 倍光速 xf 67.3; % 故障点距 A 端 km fs 1e6; % 采样率 Hz t (0:2/fs:6e-3).; % 6 ms 观测窗间隔 2 us % 宽带故障源高斯包络调制正弦中心频率 200 kHz f0 200e3; tau 1e-5; s exp(-((t - 1e-3)/tau).^2) .* exp(1j*2*pi*f0*(t - 1e-3)); s real(s); s(t 1e-3) 0; % 故障发生前的时刻置零 % 行波从故障点传播到两端 delayA xf / vp; delayB (L - xf) / vp; UA 0.5 * interp1(t, s, t - delayA, linear, 0); UB 0.5 * interp1(t, s, t - delayB, linear, 0);代码里的tau取 10 us是为了让信号具备足够宽的频带。EMTR 聚焦峰的尖锐程度直接依赖信号带宽窄带源在后续的聚焦函数里只会得到一片平坦的能量分布。实际故障暂态不是这种高斯脉冲而是叠加在工频上的衰减振荡这里生成标准数据的目的只是验证算法本身不追求精确模拟某条线路。3.3 EMTR 频域定位函数相位共轭与聚焦能量function [dist, E, x_grid] emtr_fault_loc(UA, UB, fs, L, vp) % EMTR 双端故障定位频域聚焦能量搜索 % 输入: % UA, UB : 线路两端电压时域序列列向量长度一致 % fs : 采样率 Hz % L, vp : 线路总长 km参考波速 km/s % 输出: % dist : 定位到的故障距离 km从 A 端算起 % E : 聚焦能量数组 % x_grid : 对应候选距离网格 km N length(UA); % 去直流并按 Hann 窗加权抑制频域泄漏 w 0.5 * (1 - cos(2*pi*(0:N-1)./(N-1))); UA (UA - mean(UA)) .* w; UB (UB - mean(UB)) .* w; UA_f fft(UA); UB_f fft(UB); freqs (0:N-1). * fs / N; dx 0.5; % 距离搜索步长 km x_grid (dx:dx:(L-dx)).; E zeros(size(x_grid)); for k 1:length(x_grid) x x_grid(k); % 从假设位置 x 到 A/B 端的传播算子指数取共轭 HA exp(1j * 2 * pi * freqs * x / vp); HB exp(1j * 2 * pi * freqs * (L - x) / vp); % 聚焦谱时间反转在频域等价于相位共轭 G UA_f .* conj(HA) UB_f .* conj(HB); % 忽略直流只累计正频率下的能量 idx 2:ceil(N/2); E(k) trapz(freqs(idx), abs(G(idx)).^2); end [~, imax] max(E); dist x_grid(imax); end调用方式fs 1e6; vp 2.98e5; L 100; dist emtr_fault_loc(UA, UB, fs, L, vp); fprintf(A 端到故障点距离%.2f km\n, dist); % 真实故障距离 67.3 km输出应落在 67.0 或 67.5 km 附近代码里的逻辑关键在三个地方。其一conj(HA)是频域的时间反转如果用时域写法就是flipud(UA)再做相关两者等价其二G把两端信号按假设位置做相位对齐对齐得越好|G|^2的积分越大其三trapz没有乘频率间隔df因为聚焦函数的相对大小不受影响峰值位置不变省一次乘法。提示真实录波数据通常带通道固有延时直接定位之前先做双端互相关对齐否则聚焦峰会整体偏移几百米甚至几公里。3.4 换用自己数据的接入方式data readtable(recording.csv); % 期望列: t, uA, uB t data.t; fs 1 / mean(diff(t)); % 实际采样率 UA data.uA; UB data.uB; dist emtr_fault_loc(UA, UB, fs, L, vp);PSCAD 导出的数据通常等间隔但文件名和列名各版本不一读进来后先确认时间列不要直接把第二三列当 UA、UB。Simulink 输出到工作区的timeseries需要先.Data取数据再把时间向量换算成 fs。只要最后 UA、UB 是等间隔采样的列向量函数内部不关心数据来自哪里。4. EMTR 定位精度受什么影响采样率、带宽、波速与窗口参数EMTR 的定位效果不是一个参数决定的。复现论文时最常见的问题是换了自己的线路模型以后聚焦峰变钝峰值位置偏移于是怀疑算法写错了。大部分情况下算法没错是带宽、波速和窗口没配合好。4.1 带宽决定分辨率和旁瓣先算 v/(2B)理想无损耗模型下聚焦能量峰的主瓣宽度大约是v / (2B)B 是信号的有效单边带宽。架空线波速约2.98e5 km/s如果故障暂态的带宽是 100 kHz主瓣宽度就是1.49 km。这时候搜索网格取 0.5 km 是合理的。如果信号是窄带的比如主要能量落在几 kHz 附近主瓣宽度十几公里整条线上的 E(x) 拉不开对比度算法自然失效。EMTR 的工作对象是故障行波频率从几 kHz 到几百 kHz而不是稳态工频。对录波数据先做频谱分析确认有效带宽再决定网格密度。4.2 波速设定架空线和电缆差一个量级架空线的波速接近光速约2.98e5 km/s电力电缆受绝缘材料介电常数影响波速通常在1.5e5 ~ 2.0e5 km/s。如果波速设错 5%双端 EMTR 的两个传播算子会系统性失配聚焦峰不会消失但位置会偏。工程技巧不要用资料里的经验值硬套用一条已知长度线路的开关脉冲或故障记录做一次校准反算出等效波速再投入定位。这一点对架空线混合电缆线路尤其重要混合线路的“等效波速”依赖各段长度比例没有固定的参考值。参数影响方向取值的典型做法采样率 fs越高保留的行波高频越完整但数据量增大1 MSa/s 起步按最高频率的四倍以上取信号带宽 BB 越大聚焦峰越尖锐分辨率越高先做频谱分析确认有效频段波速 vp估错导致系统性定位偏移用已知长度线路实测校准窗口长度太短丢低频太长引入母线反射覆盖 L/vp 或略长候选网格 dx不影响峰值位置的理论极限只影响搜索粒度取 v/(2B)/3 以下4.3 采样率、滤波器和窗函数怎么配合fc_lo 5e3; % 高通截止滤掉 50/60 Hz 工频 fc_hi 300e3; % 低通截止远离噪声频带 [b, a] butter(4, [fc_lo fc_hi]/(fs/2), bandpass); UA filtfilt(b, a, UA); UB filtfilt(b, a, UB);滤波器必须用filtfilt做零相位滤波。普通filter会引入与频率相关的相位延迟等于额外加了一段不确定的传播延时聚焦峰会因此偏移。窗函数的作用是抑制频谱泄漏但也会削弱数据两端的幅度窗口长度至少要覆盖故障点到最远端的往返时间否则低频成分被截断E(x) 曲线出现起伏。4.4 高阻故障、母线反射与三端线路的扩展高阻接地故障的暂态幅值小但 EMTR 用的是积分能量只要信噪比不是太低聚焦峰依然存在。这时更需要窄带滤波把能量集中在行波的频带内避免工频背景淹没有效信号。母线反射会造成虚假聚焦峰常见做法是缩短窗口只保留第一次到达波前的数据或者在线路模型里显式加入反射系数。三端线路不建议直接写三端口反演算子而是对两对端口分别做双端 EMTR得到两条聚焦曲线再取两条曲线峰值的空间交集。这个做法对短线路和中等线路都稳定定位结果还可以人工核对两条曲线的直观一致性。5. 用实验数据验证 EMTR 定位精度的三个收尾技巧5.1 缩比线路实验用同轴电缆模拟双端故障实验室没有必要搭建真实的输电线路。用一段百米级同轴电缆就能验证 EMTR 定位流程电缆末端接失配负载模拟故障点信号发生器注入宽带脉冲示波器两个通道同时采集电缆两端电压。同轴电缆的波速大约是0.66c ≈ 1.98e5 km/s代入代码里的vp。电缆很短时注意fs要足够高否则一个采样间隔就跨过好几个候选网格聚焦峰根本不可能落在正确位置。5.2 从录波文件到定位函数的完整预处理链data readtable(fault.csv); % 期望列: t, uA, uB t data.t; fs 1 / mean(diff(t)); UA data.uA - mean(data.uA); % 去直流 UB data.uB - mean(data.uB); [b, a] butter(4, [1e4 3e5]/(fs/2), bandpass); UA filtfilt(b, a, UA); UB filtfilt(b, a, UB); dist emtr_fault_loc(UA, UB, fs, L, vp);值得注意的是采样率这里不能用 CSV 文件里标注的名义值应该用mean(diff(t))反推。不少录波装置标称 1 MHz实际文件里带有微小的时钟抖动直接套标称值会引入系统误差。5.3 极性反转与延时同步用互相关给双端信号“对表”双端同步通常靠 GPS 时钟但实际记录依然可能出现固定延时。把其中一端信号相对另一端做互相关取最大相关对应的 lag 作为采样点偏移再移位对齐[c, lags] xcorr(UA, UB, coeff); [~, im] max(c); shift lags(im); % 采样点数偏移 if shift 0 UB circshift(UB, shift); else UA circshift(UA, -shift); end互相关归一化到coeff后相关峰接近 1 说明两端极性一致如果相关峰接近 -1说明两端电压极性相反需要在输入定位函数前对其中一个信号乘 -1 对齐。极性错误会让聚焦函数找错侧甚至把峰压在线路端点附近。这个检查必须在滤波器之前做还是之后做建议先滤波再互相关滤掉工频后行波波前的相关峰更干净延时估计也更稳定。本文还有配套的精品资源点击获取