地形对负地闪回击电场时频特性的影响：一种路径增量法

《Remote Sensing》：Terrain Effects on Time-Frequency Characteristics of Negative Return Stroke Electric Fields: A Path-Incremental Method

　　

闪电定位网络记录的波形混合了源变异性与传播效应，使地形相关畸变难以分离。研究人员开发了一种路径增量法，其中同一次回击被两个台站观测时，两站之间的附加路径段承载传播信号。从554次负云对地回击的4466条记录中，筛选出469对近共线台站对，按附加路径长度分层，并在各层内以数字高程模型(DEM)地形复杂度指数中位数分为平滑与崎岖两组。崎岖路径在近–远站对比中表现出上升时间(Tr)与半峰值宽度(Tw)的最清晰增大，而下降时间(Tf)与过零时间(T0)无稳定排序，表明地形主要作用于脉冲前沿与主峰。频域中，崎岖路径在50–150 kHz与100–300 kHz保留更低的归一化频带能量分数，后续回击在100–300 kHz最系统；200–500 kHz响应不稳定。独立混合效应模型将路径伸长、视线遮挡与衍射几何分别关联到后续回击100–300 kHz比值降低39.7–42.5%。12条代表性路径轴对称有限差分时域(FDTD)仿真在匹配源电流下复现了两频带的三组地电导率排序。上升时间、半峰值宽度与归一化100–300 kHz比是最一致的地形相关指标，支持地形感知波形解译。

研究背景方面，负云对地回击产生的闪电电磁脉冲(LEMP)是连接源过程、传播环境与工程应用的关键可观测量，甚低频/低频(VLF/LF)闪电定位系统依赖多站电场记录反演位置与峰值电流，但观测波形受源电流、传播距离、大地电参、地形起伏与电离层边界共同调制。已有观测与仿真多基于平地或水平分层假设，理想构型难以刻画真实山地路径的路径伸长、多径重组、视线遮挡与衍射；大样本、多站、同事件自然负地闪数据稀缺，地形对T_r、T_w、T_f、T₀的约束及其频依性与首/后续回击差异尚未明确。该研究即在此背景下开展，旨在隔离源变异性、提取地形传播信号。

关键技术方法上，研究人员采用国网VLF/LF闪电观测系统2020–2025年65站自然负地闪电场波形（采样间隔400 ns，GPS授时精度0.1 μs），筛选同事件近共线（方位差<0.001 rad）台站对469对；基于DEM计算路径几何伸长系数η、视线遮挡积分V_los、菲涅耳–基尔霍夫参数ν，PCA合成地形复杂度指数(TCI)解释85.54%方差，按附加距离三分位分层后组内以TCI中位数分平滑/崎岖亚组；时域比对近–远ΔT_r、ΔT_w、ΔT_f、ΔT₀，频域做基线扣除单边FFT并按总正频能量归一得远/近比R_b，对50–150、100–300、200–500 kHz做分层与混合效应模型（含随机截距）；另建二维轴对称球坐标FDTD，15 m网格、子网格土体体积分数映射、CPML边界，Heidler函数加MTLE通道模型，在0.001/0.005/0.01 S/m三地电导率下对两条代表路径（平滑/崎岖）的首与后续回击共12例仿真。

研究结果部分，3.1节距离依赖性分析显示50–600 km内T_r、T_f、T_w、T₀、E_os/E_peak几何均值随距增大，但单站跨事件无法分离地形，故转入同事件增量分析。3.2节时域对比表明，崎岖减平滑的中位差在首击中等距离段ΔT_r为0.101 μs（p=0.003）、ΔT_w为0.148 μs（p=0.007），后续回击长距离段ΔT_r为0.132 μs（p=0.007）；ΔT_f与ΔT₀跨层变号且CI含零，证明地形集中于前沿与主峰加宽。3.3节频带替换分析显示后续回击近站100–300 kHz能量分数约为首击2.89倍；崎岖组在100–300 kHz对后续回击三距离层R_b中位差均负（?0.340/~?0.530/~?0.418，p≤0.0303），混合效应模型给出后续回击该带η_z、V_los,z、ν_z每1σ主效应降39.7–42.5%，距离×ν_z交互项显著；200–500 kHz无一致排序。3.4节代表案例显示崎岖路径远站50–150与100–300 kHz远/近比（0.647、0.776）低于平滑路径（1.084、0.812），与全样本同向。3.5节FDTD在0.005与0.01 S/m下12种情形中崎岖路径R_b均低于平滑，100–300 kHz首击降25.6–75.8%、后续降17.5–82.1%，复现观测排序。

讨论部分指出，T_r与T_w对几何相关反射/衍射敏感而尾部参数受源与干涉混淆；100–300 kHz是能量可测性与快变成分敏感性的最优平衡带，但属样本经验界非普适边界；后续回击因前沿陡、初始高频占比大，地形对比更易析出；FDTD证地形引起频率选择性重分布而非纯宽带衰减，200–500 kHz受地电导率调制。地形致前沿变宽使恒比阈值越靠峰到达时偏差越大（5%→50%阈值中位偏差0.460→1.072 μs，等效测距误138→321 m），早取上升沿前段可抑定位偏；峰值电流反演因场–流映射畸变需按地形与地电导率校准，绝对0–50 kHz能比R_abs崎岖廓线0.799–0.899高于平滑0.555–0.583，每百公里附加路径首击峰流偏倚多6.4–21.8百分点、后续多12.8–30.4百分点。

结论部分翻译：研究人员利用554次负云对地回击4466条电场记录，以同事件近共线台站对隔离地形效应，使每次比对共享单一源。崎岖附加路径引起上升时间与半峰值宽度更大增量，下降时间与过零时间无稳定地形排序，地形特征集中于脉冲前沿与主峰。频谱中崎岖路径在50–150与100–300 kHz保留更小归一化能量份额，后续回击在100–300 kHz最一致，200–500 kHz随地电导率变化。上升时间、半峰值宽度与100–300 kHz能量保留比因此构成传播相关波形修正的实用指标，尤以后续回击为甚。因修正集中于前沿与主峰，其渗入两类派生量：恒比标记近峰时到达时偏差增大，故上升沿早取标记减小定位误；峰值电流反演经场–流映射偏倚，需地形与地电导率知情校准。该证据未分辨单一传播机制或定义普适地形阈值，仿真为无噪基线，校准多站观测、空间变地参及三维建模仍待补入业务修正前验证。论文发表于《Remote Sensing》。