简介一份聚焦GIS局部放电超高频电磁波传播特性的学术文献面向电力设备状态检测工程师、高电压专业研究者及电气工程学生旨在帮助理解UHF法检测局部放电的原理与仿真路径。文章基于时域有限差分法FDTD开展仿真系统阐述了UHF信号在GIS中传播的延迟效应以及经过绝缘子和T型接头时的衰减规律具体给出平均衰减3dB与支路衰减10dB等关键数据并讨论谐振作用对信号传播的减缓影响可作为局部放电定位、传感器优化布置与监测系统设计的参考依据。整包为单个PDF文档大小约229KB内容清晰完整下载后可直接阅读。目前已有110人学习使用对正在开展GIS局部放电研究或工程应用的读者具有较高参考价值。1. 从FDTD仿真看GIS局部放电UHF信号的三维传播GIS局部放电检测里最闹心的问题不是传感器不够灵敏而是你根本不知道测到的信号在GIS腔体里绕了多少路、衰减了多少、在哪个接头又反弹了一次。这篇2006年发表在《高电压技术》的论文用FDTD时域有限差分法把GIS内部UHF电磁波的传播路径、延迟、谐振和衰减做了完整的仿真计算给出了几个可以直接用于工程判断的量化结论信号通过绝缘子平均衰减3dB、通过T型接头拐弯支路平均衰减10dB、导体损耗约3~5dB/km。这些数字对今天做UHF传感器布点、放电源定位和阈值整定仍然有参考价值。文章适合做GIS在线监测、局部放电模式识别和电磁场仿真的工程技术人员阅读尤其是需要把仿真结果换算成现场检测方案的人。2. FDTD建模与仿真参数设置从同轴腔体到高斯脉冲源的逼近2.1 为什么选时域有限差分法而不是频域方法GIS腔体本质上是同轴波导结构内导体是高压导体外导体是接地外壳内部填充SF6气体。局部放电产生的电磁波频带极宽论文中给出现场检测到上升时间小到70ps的放电脉冲对应的频谱分量可以激励到1.5GHz以上UHF检测频段集中在500~2000MHz。这个频带覆盖了TEM模、TE模和TM模多个模式频域方法对每个模式都要单独求解、再叠加处理起来很麻烦而且放电本身是瞬态脉冲用频域方法描述激励源也不自然。FDTD直接在时域上对Maxwell方程组做中心差分离散一次计算就能得到整个频带的响应对瞬态信号天然适配。论文引用Taflove的经典教材作为方法依据剖分采用立方体胞元、边长10mm实际上这个尺寸在自由空间对应约30GHz的频率上限覆盖UHF频段绰绰有余。2.2 同轴GIS腔体的几何建模与网格剖分论文第一个仿真模型是时间延迟模型关键几何参数如下参数数值说明内导体外径 a40 mm高压导杆半径外导体内径 b300 mm外壳内径腔体全长 L4 m模型总长度放电点位置距左侧绝缘子 2 m金属突出物放电突出物长度10 mm模拟尖端放电胞元边长10 mm立方体剖分网格网格规模600×60×60总计算场区这个剖分规模在2006年属于中等偏大的计算量60×60的截面网格意味着同轴截面内径方向30个网格、周向60度分60格空间分辨率10mm对300mm外径的腔体来说曲率建模基本够用。现在个人工作站跑同样的规模用几分钟到十几分钟。网格剖分的关键是内导体和外导体的阶梯近似。FDTD用正交网格逼近圆柱面时会产生台阶误差论文的作法是在剖分截面图上直接把同轴结构映射到网格上没有做共形网格处理。这种简化对UHF频段的工程仿真影响不大因为感兴趣的信号波长在150~600mm之间500~2000MHz10mm的台阶远小于波长。2.3 高斯脉冲源与FDTD计算流程放电源用高斯电流脉冲激励这是局部放电仿真里最常见的激励模型。脉冲幅值10mA、宽度0.6ns、峰值出现在1ns时刻。当放电脉冲为高斯脉冲时辐射电场与电流脉冲呈微分关系所以接收点看到的电场波形是正负对称的两个尖峰。后处理时从FDTD算出的电场分量里提取有效分量E_φ因为接收点置于GIS腔体顶部它对输出电压信号起作用的是切向电场分量。# 高斯脉冲源定义Python示例用于FDTD前置处理器 import numpy as np def gaussian_pulse(t, A0.01, t01e-9, sigma0.25e-9): 高斯电流脉冲 A: 峰值电流 0.01 A (10 mA) t0: 峰值出现时刻 1 ns sigma: 脉冲宽度参数0.6ns半宽对应 sigma ≈ 0.25ns return A * np.exp(-((t - t0) ** 2) / (2 * sigma ** 2)) # 验证采样率 100 THz对应时间步长 0.01 ps t np.linspace(0, 5e-9, 5000) i_pulse gaussian_pulse(t) # 幅值归一化检查 assert np.isclose(np.max(i_pulse), 0.01), 峰值应为 10 mA代码里的sigma参数不是脉冲半宽而是高斯函数的标准差。论文说脉冲宽度0.6ns对应的是从起始到结束的主要能量区间高斯脉冲的工程习惯是取±2.4σ约98%能量作为有效宽度0.25ns的σ对应约1.2ns的完整包络。实际建模时反过来先用电荷转移量估算脉冲面积再反推幅值和宽度。FDTD主循环的稳定性受CFL条件约束时间步长不能超过网格对角穿越时间。10mm网格下最大时间步长约19.2ps论文的仿真时长覆盖100ns对应约5200个时间步这是早年做一次仿真需要等待的原因之一。2.4 吸收边界与计算区域处理FDTD仿真区域必须截断论文的模型两端是绝缘子属于介质不连续面不能简单用PML吸收。实际做法是轴向两端设置PML吸收层论文没有明确说明但按当时主流做法推断同时在被绝缘子封闭的腔体内部波的多次反射是真实物理过程不能吸收掉。这里有个容易踩的坑如果两端直接设PML谐振腔效应就没了算出来的衰减会偏大如果两端用理想导体PEC封闭谐振频率和实际GIS又对不上。论文的做法是把绝缘子建模成介质层让绝缘子成为反射和透射的真实边界腔体内部自然形成驻波。这种处理比整体PML截断更贴近物理场景代价是计算区域必须包含完整的绝缘子结构和足够长的空气段。3. 时间延迟效应UHF信号传播路径与局部放电源定位3.1 延迟效应的物理机制与仿真验证时间延迟的机理不复杂电磁波从放电源到接收点要走有限距离速度近似光速所以不同位置的传感器收到信号的时间不同。论文用两接收点验证了这个效应——接收点A距放电点1m接收点B距A点1m都安装在腔体顶部。关键参数校验放电电流脉冲峰值出现在1ns电磁波从放电点传到A点需要约3.3ns1米除以光速因此A点应该在3.8ns附近收到第一个脉冲。论文的FDTD计算结果与这个理论值完全一致B点延迟类似。这是个必要的正确性校验——如果仿真结果跟理论传播时间差得太多说明剖分、激励或边界设置有问题。这段校验逻辑是全文方法论的骨架后续做任何FDTD仿真都应该先做这类解析对比。3.2 多路径传播与信号持续时间论文指出上述分析只考虑了直线传播路径实际上放电点和接收点之间有无限多条传播路径电磁波可以通过腔体壁反射、绝缘子透射、T型接头拐弯等不同路径到达接收点。这导致接收信号比原始放电脉冲持续的时间长得多论文提到能持续数μs而放电脉冲本身只有ns级。这个特征解释了为什么UHF信号波形看起来拖尾很重。不同路径的信号在不同时间、以不同幅值到达纵向上在腔体端部反射发生之前横向截面内就已出现论文引用的二维谐振现象。这在实际检测中的意义是UHF波形不能简单按到达时间猜放电源位置必须先排除多路径干扰。3.3 基于时间差的定位算法思路延迟效应的工程价值是支撑局部放电源定位。做法是在GIS不同位置安装多个UHF传感器利用信号到达时间差TDOA确定位置。设两个传感器位置分别为S1、S2放电源位置为P到达时间差Δt与位置关系满足d1 c * t1 # 放电源到传感器1的距离 d2 c * t2 # 放电源到传感器2的距离 Δt t1 - t2 # 到达时间差 Δd c * Δt d1 - d2 # 距离差放电源位置在双曲面上双曲面的焦点就是两个传感器位置。工程上布4个以上传感器用非线性最小二乘求解# TDOA定位求解简化实现 from scipy.optimize import least_squares def tdoa_residual(p, sensor_positions, measured_tdoa): sensor_positions: 4个传感器坐标形状 (4,3) measured_tdoa: 相对第一个传感器的到达时间差数组 distances np.linalg.norm(sensor_positions - p, axis1) predicted_tdoa distances - distances[0] # 参考传感器为0号 return predicted_tdoa[1:] / 3e8 - measured_tdoa # 传感器坐标示例沿腔体轴向布4个点 sensors np.array([ [0.0, 0.15, 0.0], [1.5, 0.15, 0.0], [3.0, 0.15, 0.0], [4.0, 0.15, 0.0] ]) # 求解放电源位置 result least_squares(tdoa_residual, x0[2.0, 0.1, 0.0], args(sensors, measured_tdoa))这段代码里measured_tdoa是从各传感器采集信号里提取的首脉冲到达时刻差。但注意论文的仿真结论提示我们通道间信号可能来自不同传播路径所以后续要加上首波识别。实际工程里TDOA定位的精度受限于多路径干扰。放电源到传感器间存在反射路径反射波可能比直射波更强如果触发点设在不合适的阈值上首波识别会出错。这时候要做的是在多个传感器上先做信号到达时间标注去掉明显不符合物理传播时间的野值再用最小二乘平滑。3.4 延迟精度的影响参数论文摘要里有一句容易被忽略的话UHF信号受放电源径向位置、放电通道长度、放电源形状以及放电源与传感器夹角的影响。这意味着时间延迟不是单纯的距离除以光速放电源的辐射方向图会改变不同方向上的有效信号强度。放电通道越长激励出的电磁波模式越多第一个到达的波前可能不是能量最强的波前放电源形状决定了辐射方向图贴着内导体表面放电和悬浮电位放电的辐射特性完全不同。这对定位系统的要求是不能只靠单次触发的时间差应该在多个放电脉冲上做统计平均用一致性好、重复性高的首波时刻作为TDOA输入。4. 反射-谐振与衰减特性波形畸变的物理机制4.1 波阻抗不连续点的反射与透射GIS里的绝缘子、T型接头、弯角、TA电流互感器/TV电压互感器连接处都是波阻抗不连续点。电磁波在这些位置发生部分反射、部分透射入射波与反射波叠加形成谐振当不连续点间的距离等于半波长的整数倍时形成驻波。这个现象对传感器有两个直接影响一是反射波叠加可能增强或抵消特定频率的信号导致传感器收到的信号幅度与放电源真实强度不对应二是谐振效应让信号衰减变慢波形拖尾明显这对时域特征提取有干扰。4.2 谐振频率与GIS腔体尺寸的关系论文给出同轴谐振腔的谐振频率公式$f_{mnp} \sqrt{f_{c,mn}^2 (pv/2L)^2}$其中$f_{c,mn}$是特定模式的截止频率m、n是横向模式指数p是纵向谐振阶数v是波速L是谐振腔长度。可以看出谐振频率由横向截止频率和纵向谐振条件共同决定。谐振仿真模型参数a40mmb380mm两绝缘子间距50cm总长L1.5m。FDTD结果显示信号包络衰减常数是μs量级谐振过程让电场波形变得非常复杂。频谱图显示某些频点能量集中在很窄的带宽内这些频点对应的就是谐振频率。对检测来说谐振频带内的信号是由多路径叠加形成的稳态分布不能直接把幅度当放电强度用但是谐振频率本身又包含了几何结构信息可以利用它对放电位置做粗定位。4.3 均匀腔体中的模式衰减论文把衰减分成两类均匀腔体内的固有损耗和不连续点处产生的损耗。前者包括金属导体有限电导率引起的损耗、填充介质不纯净引起的损耗、以及传播模式低于截止频率时的截止衰减。SF6气体介质损耗很小施工中通常忽略只算导体损耗。FDTD计算的长直腔体单位长度衰减曲线有几个工程要点模式衰减特征TEM衰减最小传输主模式TE11截止频率附近衰减峰值明显TE12、TE21峰值频率接近TE11带宽很窄TM01、TM11高频段衰减较大各模式衰减曲线的峰值对应频率接近原因是在截止频率附近信号在腔体中传播一定距离要经历大量反射单位长度内反射次数增加等效衰减变大。有限电导率金属导体产生的信号衰减系数约3~5dB/km这个数值对短距离GIS几十米来说不是主要矛盾真正要关注的是绝缘子和T型接头。4.4 绝缘子与T型接头的定量衰减结果绝缘子衰减模型采用平板绝缘子简化材料相对介电常数r6厚度5cm。仿真结果显示信号经过绝缘子后幅值衰减且A处第一对脉冲之后的波动是信号在绝缘子处的反射频域上某些频率增益大于0dB是谐振增强的结果。信号经绝缘子后平均衰减3dB。T型接头用三探头模型验证A探头在直通臂上B、C探头在拐弯支路上。结果非常直观拐弯路径比直线路径的衰减严重得多超高频范围内拐弯路径平均衰减10dB、直线路径平均衰减2dB。传播路径平均衰减说明直通路径过T型接头2 dB与绝缘子衰减同量级拐弯路径过T型接头10 dB分支支路信号损失严重绝缘子εr6厚度5cm3 dB反射介质损耗均匀腔体导体损耗3~5 dB/km长距离才需关注10dB是功率比10倍的差距这解释了为什么在T型接头分支末端装传感器经常收不到信号。实际GIS母线间隔通常有好几个盆式绝缘子和至少两三个T型接头信号从放电源传到传感器可能要穿越5~10个不连续点累积衰减可能达到30~50dB这对传感器灵敏度和信号调理电路的动态范围提出了很高要求。5. 传感器布点与放电源定位的实用技巧5.1 时间差定位里的首波识别时间差计算的核心是准确标注每个通道首脉冲的到达时刻。受多路径和噪声影响直接用固定阈值触发不可靠。常见做法是先用小波变换或带通滤波增强信噪比再用能量累积曲线的前沿作为到达时刻估计。论文中放电脉冲上升时间可达70ps但UHF传感器带宽是500~2000MHz实测前沿受传感器响应限制通常在几百ps到1ns之间时间差分辨率约0.1~0.3ns对应定位误差3~9cm。5.2 布点原则与衰减预算根据论文提供的衰减数据可以初步估算传感器的覆盖半径。设传感器灵敏度要求信号到达时信噪比不低于6dB路径上需要穿过2个绝缘子和1个T型接头拐弯那么检测链路总衰减预算为绝缘子2 × 3 dB 6 dB T型接头拐弯1 × 10 dB 10 dB 导体损耗假设30m路径约0.1~0.15 dB 其他裕量6 dB 总衰减预算 ≈ 22 dB这意味着传感器系统至少要有22dB的动态裕量。如果信号源放电量很小pC级这个裕量可能不够需要在T型接头分支上额外加装传感器。论文数据直接支撑了拐弯处必须有独立传感器这个布点建议。5.3 利用谐振效应增强特定频段检测谐振让信号衰减变慢这是劣势也是优势。在GIS工程上可以故意利用腔体谐振特性选择传感器的安装位置和检测频带。如果已知某段GIS的几何尺寸可以先算谐振频率然后把UHF传感器的检测频段对准谐振频率可以获得比非谐振频点更高的等效灵敏度。这相当于用腔体做信号放大器代价是幅度与放电量的线性关系会失真。5.4 数据后处理的最小实现# 从FDTD或实测数据提取衰减量的最小实现 import numpy as np def estimate_attenuation_db(signal_ref, signal_target, fs): 计算目标信号相对参考信号的能量衰减 signal_ref: 参考点信号如A探头 signal_target: 目标点信号如B探头 fs: 采样率Hz def energy(sig): return np.sum(np.abs(sig) ** 2) / fs e_ref energy(signal_ref) e_tgt energy(signal_target) # 保护避免除零 if e_ref 1e-12: return np.inf return 10 * np.log10(e_tgt / e_ref)这段代码可以直接套用于FDTD输出的时域波形计算出两点之间的综合衰减。注意频带要限定在500~2000MHz之间否则带外噪声会干扰能量估算。对实测信号采样率至少要4GSa/s才能覆盖2000MHz的奈奎斯特频率。论文给的核心结论在工程布点里怎么用总结成一句话延迟定位定大致区间、衰减数据定传感器个数、谐振频谱定检测频带。这套思路放到今天依然是GIS-UHF检测系统设计的基本框架。本文还有配套的精品资源点击获取 FDTD仿真揭示GIS局部放电UHF信号传播衰减与定位布点
简介一份聚焦GIS局部放电超高频电磁波传播特性的学术文献面向电力设备状态检测工程师、高电压专业研究者及电气工程学生旨在帮助理解UHF法检测局部放电的原理与仿真路径。文章基于时域有限差分法FDTD开展仿真系统阐述了UHF信号在GIS中传播的延迟效应以及经过绝缘子和T型接头时的衰减规律具体给出平均衰减3dB与支路衰减10dB等关键数据并讨论谐振作用对信号传播的减缓影响可作为局部放电定位、传感器优化布置与监测系统设计的参考依据。整包为单个PDF文档大小约229KB内容清晰完整下载后可直接阅读。目前已有110人学习使用对正在开展GIS局部放电研究或工程应用的读者具有较高参考价值。1. 从FDTD仿真看GIS局部放电UHF信号的三维传播GIS局部放电检测里最闹心的问题不是传感器不够灵敏而是你根本不知道测到的信号在GIS腔体里绕了多少路、衰减了多少、在哪个接头又反弹了一次。这篇2006年发表在《高电压技术》的论文用FDTD时域有限差分法把GIS内部UHF电磁波的传播路径、延迟、谐振和衰减做了完整的仿真计算给出了几个可以直接用于工程判断的量化结论信号通过绝缘子平均衰减3dB、通过T型接头拐弯支路平均衰减10dB、导体损耗约3~5dB/km。这些数字对今天做UHF传感器布点、放电源定位和阈值整定仍然有参考价值。文章适合做GIS在线监测、局部放电模式识别和电磁场仿真的工程技术人员阅读尤其是需要把仿真结果换算成现场检测方案的人。2. FDTD建模与仿真参数设置从同轴腔体到高斯脉冲源的逼近2.1 为什么选时域有限差分法而不是频域方法GIS腔体本质上是同轴波导结构内导体是高压导体外导体是接地外壳内部填充SF6气体。局部放电产生的电磁波频带极宽论文中给出现场检测到上升时间小到70ps的放电脉冲对应的频谱分量可以激励到1.5GHz以上UHF检测频段集中在500~2000MHz。这个频带覆盖了TEM模、TE模和TM模多个模式频域方法对每个模式都要单独求解、再叠加处理起来很麻烦而且放电本身是瞬态脉冲用频域方法描述激励源也不自然。FDTD直接在时域上对Maxwell方程组做中心差分离散一次计算就能得到整个频带的响应对瞬态信号天然适配。论文引用Taflove的经典教材作为方法依据剖分采用立方体胞元、边长10mm实际上这个尺寸在自由空间对应约30GHz的频率上限覆盖UHF频段绰绰有余。2.2 同轴GIS腔体的几何建模与网格剖分论文第一个仿真模型是时间延迟模型关键几何参数如下参数数值说明内导体外径 a40 mm高压导杆半径外导体内径 b300 mm外壳内径腔体全长 L4 m模型总长度放电点位置距左侧绝缘子 2 m金属突出物放电突出物长度10 mm模拟尖端放电胞元边长10 mm立方体剖分网格网格规模600×60×60总计算场区这个剖分规模在2006年属于中等偏大的计算量60×60的截面网格意味着同轴截面内径方向30个网格、周向60度分60格空间分辨率10mm对300mm外径的腔体来说曲率建模基本够用。现在个人工作站跑同样的规模用几分钟到十几分钟。网格剖分的关键是内导体和外导体的阶梯近似。FDTD用正交网格逼近圆柱面时会产生台阶误差论文的作法是在剖分截面图上直接把同轴结构映射到网格上没有做共形网格处理。这种简化对UHF频段的工程仿真影响不大因为感兴趣的信号波长在150~600mm之间500~2000MHz10mm的台阶远小于波长。2.3 高斯脉冲源与FDTD计算流程放电源用高斯电流脉冲激励这是局部放电仿真里最常见的激励模型。脉冲幅值10mA、宽度0.6ns、峰值出现在1ns时刻。当放电脉冲为高斯脉冲时辐射电场与电流脉冲呈微分关系所以接收点看到的电场波形是正负对称的两个尖峰。后处理时从FDTD算出的电场分量里提取有效分量E_φ因为接收点置于GIS腔体顶部它对输出电压信号起作用的是切向电场分量。# 高斯脉冲源定义Python示例用于FDTD前置处理器 import numpy as np def gaussian_pulse(t, A0.01, t01e-9, sigma0.25e-9): 高斯电流脉冲 A: 峰值电流 0.01 A (10 mA) t0: 峰值出现时刻 1 ns sigma: 脉冲宽度参数0.6ns半宽对应 sigma ≈ 0.25ns return A * np.exp(-((t - t0) ** 2) / (2 * sigma ** 2)) # 验证采样率 100 THz对应时间步长 0.01 ps t np.linspace(0, 5e-9, 5000) i_pulse gaussian_pulse(t) # 幅值归一化检查 assert np.isclose(np.max(i_pulse), 0.01), 峰值应为 10 mA代码里的sigma参数不是脉冲半宽而是高斯函数的标准差。论文说脉冲宽度0.6ns对应的是从起始到结束的主要能量区间高斯脉冲的工程习惯是取±2.4σ约98%能量作为有效宽度0.25ns的σ对应约1.2ns的完整包络。实际建模时反过来先用电荷转移量估算脉冲面积再反推幅值和宽度。FDTD主循环的稳定性受CFL条件约束时间步长不能超过网格对角穿越时间。10mm网格下最大时间步长约19.2ps论文的仿真时长覆盖100ns对应约5200个时间步这是早年做一次仿真需要等待的原因之一。2.4 吸收边界与计算区域处理FDTD仿真区域必须截断论文的模型两端是绝缘子属于介质不连续面不能简单用PML吸收。实际做法是轴向两端设置PML吸收层论文没有明确说明但按当时主流做法推断同时在被绝缘子封闭的腔体内部波的多次反射是真实物理过程不能吸收掉。这里有个容易踩的坑如果两端直接设PML谐振腔效应就没了算出来的衰减会偏大如果两端用理想导体PEC封闭谐振频率和实际GIS又对不上。论文的做法是把绝缘子建模成介质层让绝缘子成为反射和透射的真实边界腔体内部自然形成驻波。这种处理比整体PML截断更贴近物理场景代价是计算区域必须包含完整的绝缘子结构和足够长的空气段。3. 时间延迟效应UHF信号传播路径与局部放电源定位3.1 延迟效应的物理机制与仿真验证时间延迟的机理不复杂电磁波从放电源到接收点要走有限距离速度近似光速所以不同位置的传感器收到信号的时间不同。论文用两接收点验证了这个效应——接收点A距放电点1m接收点B距A点1m都安装在腔体顶部。关键参数校验放电电流脉冲峰值出现在1ns电磁波从放电点传到A点需要约3.3ns1米除以光速因此A点应该在3.8ns附近收到第一个脉冲。论文的FDTD计算结果与这个理论值完全一致B点延迟类似。这是个必要的正确性校验——如果仿真结果跟理论传播时间差得太多说明剖分、激励或边界设置有问题。这段校验逻辑是全文方法论的骨架后续做任何FDTD仿真都应该先做这类解析对比。3.2 多路径传播与信号持续时间论文指出上述分析只考虑了直线传播路径实际上放电点和接收点之间有无限多条传播路径电磁波可以通过腔体壁反射、绝缘子透射、T型接头拐弯等不同路径到达接收点。这导致接收信号比原始放电脉冲持续的时间长得多论文提到能持续数μs而放电脉冲本身只有ns级。这个特征解释了为什么UHF信号波形看起来拖尾很重。不同路径的信号在不同时间、以不同幅值到达纵向上在腔体端部反射发生之前横向截面内就已出现论文引用的二维谐振现象。这在实际检测中的意义是UHF波形不能简单按到达时间猜放电源位置必须先排除多路径干扰。3.3 基于时间差的定位算法思路延迟效应的工程价值是支撑局部放电源定位。做法是在GIS不同位置安装多个UHF传感器利用信号到达时间差TDOA确定位置。设两个传感器位置分别为S1、S2放电源位置为P到达时间差Δt与位置关系满足d1 c * t1 # 放电源到传感器1的距离 d2 c * t2 # 放电源到传感器2的距离 Δt t1 - t2 # 到达时间差 Δd c * Δt d1 - d2 # 距离差放电源位置在双曲面上双曲面的焦点就是两个传感器位置。工程上布4个以上传感器用非线性最小二乘求解# TDOA定位求解简化实现 from scipy.optimize import least_squares def tdoa_residual(p, sensor_positions, measured_tdoa): sensor_positions: 4个传感器坐标形状 (4,3) measured_tdoa: 相对第一个传感器的到达时间差数组 distances np.linalg.norm(sensor_positions - p, axis1) predicted_tdoa distances - distances[0] # 参考传感器为0号 return predicted_tdoa[1:] / 3e8 - measured_tdoa # 传感器坐标示例沿腔体轴向布4个点 sensors np.array([ [0.0, 0.15, 0.0], [1.5, 0.15, 0.0], [3.0, 0.15, 0.0], [4.0, 0.15, 0.0] ]) # 求解放电源位置 result least_squares(tdoa_residual, x0[2.0, 0.1, 0.0], args(sensors, measured_tdoa))这段代码里measured_tdoa是从各传感器采集信号里提取的首脉冲到达时刻差。但注意论文的仿真结论提示我们通道间信号可能来自不同传播路径所以后续要加上首波识别。实际工程里TDOA定位的精度受限于多路径干扰。放电源到传感器间存在反射路径反射波可能比直射波更强如果触发点设在不合适的阈值上首波识别会出错。这时候要做的是在多个传感器上先做信号到达时间标注去掉明显不符合物理传播时间的野值再用最小二乘平滑。3.4 延迟精度的影响参数论文摘要里有一句容易被忽略的话UHF信号受放电源径向位置、放电通道长度、放电源形状以及放电源与传感器夹角的影响。这意味着时间延迟不是单纯的距离除以光速放电源的辐射方向图会改变不同方向上的有效信号强度。放电通道越长激励出的电磁波模式越多第一个到达的波前可能不是能量最强的波前放电源形状决定了辐射方向图贴着内导体表面放电和悬浮电位放电的辐射特性完全不同。这对定位系统的要求是不能只靠单次触发的时间差应该在多个放电脉冲上做统计平均用一致性好、重复性高的首波时刻作为TDOA输入。4. 反射-谐振与衰减特性波形畸变的物理机制4.1 波阻抗不连续点的反射与透射GIS里的绝缘子、T型接头、弯角、TA电流互感器/TV电压互感器连接处都是波阻抗不连续点。电磁波在这些位置发生部分反射、部分透射入射波与反射波叠加形成谐振当不连续点间的距离等于半波长的整数倍时形成驻波。这个现象对传感器有两个直接影响一是反射波叠加可能增强或抵消特定频率的信号导致传感器收到的信号幅度与放电源真实强度不对应二是谐振效应让信号衰减变慢波形拖尾明显这对时域特征提取有干扰。4.2 谐振频率与GIS腔体尺寸的关系论文给出同轴谐振腔的谐振频率公式$f_{mnp} \sqrt{f_{c,mn}^2 (pv/2L)^2}$其中$f_{c,mn}$是特定模式的截止频率m、n是横向模式指数p是纵向谐振阶数v是波速L是谐振腔长度。可以看出谐振频率由横向截止频率和纵向谐振条件共同决定。谐振仿真模型参数a40mmb380mm两绝缘子间距50cm总长L1.5m。FDTD结果显示信号包络衰减常数是μs量级谐振过程让电场波形变得非常复杂。频谱图显示某些频点能量集中在很窄的带宽内这些频点对应的就是谐振频率。对检测来说谐振频带内的信号是由多路径叠加形成的稳态分布不能直接把幅度当放电强度用但是谐振频率本身又包含了几何结构信息可以利用它对放电位置做粗定位。4.3 均匀腔体中的模式衰减论文把衰减分成两类均匀腔体内的固有损耗和不连续点处产生的损耗。前者包括金属导体有限电导率引起的损耗、填充介质不纯净引起的损耗、以及传播模式低于截止频率时的截止衰减。SF6气体介质损耗很小施工中通常忽略只算导体损耗。FDTD计算的长直腔体单位长度衰减曲线有几个工程要点模式衰减特征TEM衰减最小传输主模式TE11截止频率附近衰减峰值明显TE12、TE21峰值频率接近TE11带宽很窄TM01、TM11高频段衰减较大各模式衰减曲线的峰值对应频率接近原因是在截止频率附近信号在腔体中传播一定距离要经历大量反射单位长度内反射次数增加等效衰减变大。有限电导率金属导体产生的信号衰减系数约3~5dB/km这个数值对短距离GIS几十米来说不是主要矛盾真正要关注的是绝缘子和T型接头。4.4 绝缘子与T型接头的定量衰减结果绝缘子衰减模型采用平板绝缘子简化材料相对介电常数r6厚度5cm。仿真结果显示信号经过绝缘子后幅值衰减且A处第一对脉冲之后的波动是信号在绝缘子处的反射频域上某些频率增益大于0dB是谐振增强的结果。信号经绝缘子后平均衰减3dB。T型接头用三探头模型验证A探头在直通臂上B、C探头在拐弯支路上。结果非常直观拐弯路径比直线路径的衰减严重得多超高频范围内拐弯路径平均衰减10dB、直线路径平均衰减2dB。传播路径平均衰减说明直通路径过T型接头2 dB与绝缘子衰减同量级拐弯路径过T型接头10 dB分支支路信号损失严重绝缘子εr6厚度5cm3 dB反射介质损耗均匀腔体导体损耗3~5 dB/km长距离才需关注10dB是功率比10倍的差距这解释了为什么在T型接头分支末端装传感器经常收不到信号。实际GIS母线间隔通常有好几个盆式绝缘子和至少两三个T型接头信号从放电源传到传感器可能要穿越5~10个不连续点累积衰减可能达到30~50dB这对传感器灵敏度和信号调理电路的动态范围提出了很高要求。5. 传感器布点与放电源定位的实用技巧5.1 时间差定位里的首波识别时间差计算的核心是准确标注每个通道首脉冲的到达时刻。受多路径和噪声影响直接用固定阈值触发不可靠。常见做法是先用小波变换或带通滤波增强信噪比再用能量累积曲线的前沿作为到达时刻估计。论文中放电脉冲上升时间可达70ps但UHF传感器带宽是500~2000MHz实测前沿受传感器响应限制通常在几百ps到1ns之间时间差分辨率约0.1~0.3ns对应定位误差3~9cm。5.2 布点原则与衰减预算根据论文提供的衰减数据可以初步估算传感器的覆盖半径。设传感器灵敏度要求信号到达时信噪比不低于6dB路径上需要穿过2个绝缘子和1个T型接头拐弯那么检测链路总衰减预算为绝缘子2 × 3 dB 6 dB T型接头拐弯1 × 10 dB 10 dB 导体损耗假设30m路径约0.1~0.15 dB 其他裕量6 dB 总衰减预算 ≈ 22 dB这意味着传感器系统至少要有22dB的动态裕量。如果信号源放电量很小pC级这个裕量可能不够需要在T型接头分支上额外加装传感器。论文数据直接支撑了拐弯处必须有独立传感器这个布点建议。5.3 利用谐振效应增强特定频段检测谐振让信号衰减变慢这是劣势也是优势。在GIS工程上可以故意利用腔体谐振特性选择传感器的安装位置和检测频带。如果已知某段GIS的几何尺寸可以先算谐振频率然后把UHF传感器的检测频段对准谐振频率可以获得比非谐振频点更高的等效灵敏度。这相当于用腔体做信号放大器代价是幅度与放电量的线性关系会失真。5.4 数据后处理的最小实现# 从FDTD或实测数据提取衰减量的最小实现 import numpy as np def estimate_attenuation_db(signal_ref, signal_target, fs): 计算目标信号相对参考信号的能量衰减 signal_ref: 参考点信号如A探头 signal_target: 目标点信号如B探头 fs: 采样率Hz def energy(sig): return np.sum(np.abs(sig) ** 2) / fs e_ref energy(signal_ref) e_tgt energy(signal_target) # 保护避免除零 if e_ref 1e-12: return np.inf return 10 * np.log10(e_tgt / e_ref)这段代码可以直接套用于FDTD输出的时域波形计算出两点之间的综合衰减。注意频带要限定在500~2000MHz之间否则带外噪声会干扰能量估算。对实测信号采样率至少要4GSa/s才能覆盖2000MHz的奈奎斯特频率。论文给的核心结论在工程布点里怎么用总结成一句话延迟定位定大致区间、衰减数据定传感器个数、谐振频谱定检测频带。这套思路放到今天依然是GIS-UHF检测系统设计的基本框架。本文还有配套的精品资源点击获取