张源源 , 庞华基 , 刘 钊 , 宋 琳 , 邓 猛

(1. 青岛市气象灾害防御技术中心, 山东 青岛 266003; 2. 青岛市气象灾害防御工程技术研究中心, 山东 青岛 266003; 3. 山东省气象局气象灾害防御技术中心, 山东 济南 250000)

目前, 对闪电的定位主要采用时差法, 即利用闪电电磁场波峰到达多个测站的时间来计算闪电发生的位置, 我国气象和电力系统建设的闪电定位网均是采用这种定位方法[1-3]。多人[4-10]研究揭示, 当雷电电磁场沿地表传播时, 地表的起伏和下垫面电导率的不同, 均会对雷电波传输产生影响。Weidman等[11-12]和 Willett[13]等观测表明, 随着观测距离的增大, 海面对雷电电磁波谱具有低通滤波功能, 超过1 MHz的高频分量按照频率平方的倒数迅速衰减, 超过20 MHz的则几乎观测不到。Ming和Cooray等[14]和张其林等[15]利用数值模式从理论上验证了海面对雷电电磁场传播的影响。李东帅等[16-18]利用时域有限差分算法(FDTD)分析起伏地表的雷电电磁场传播特性及其对闪电定位系统的影响, 指出由起伏地表引起的雷电电磁波上升沿时间的延迟, 会对基于时差法的闪电定位系统产生影响。

综上, 前人对陆地对雷电电磁场的传输影响以及其对闪电定位误差的影响进行了卓有成效的研究, 但对海面对雷电电磁场的传输影响及其对闪电定位误差的影响的研究较少, 本文拟在前人的基础上, 利用Wait近似算法分析海浪对雷电波传输的影响, 并进一步分析这种影响对基于时差法定位闪电带来的误差。

近海的闪电探测和定位对港口作业、船舶的航行及经济活动十分重要, 因此研究近海起伏海面对闪电定位系统的影响具有重要意义。本文选取一次雷暴天气过程, 将实测的浮标站风场数据输入模式算法中, 有针对性地研究海面对雷电电磁传播的影响, 为近海雷电探测及定位误差修订提供参考。

1 起伏海面的雷电垂直电场近似算法

1.1 计算起伏海面的雷电垂直电场

假定地面的电导率无限大, 闪电通道垂直于地面,则地表任意一点处的垂直电场的时域表达式为[19-20]:其中,0ε为自由空间的介电常数, z′是地面观测到闪电回击电流脉冲波头的位置, tb(z′)是地面观测到闪电回击电流脉冲从地面到目前波头位置的时间,tb(z′)=(z′2+r2)1/2/c, r是观测点距雷电回击通道的水平距离, c是光速。L′(t)是t时刻观测到的回击通道长度,v是闪电回击速度, i(z′,t-R(z′)/c)是回击电流沿通道的分布, ()zθ′是z′与观测点之间的角度, R(z′)是z′与观测点之间的距离。雷电回击过程如图1所示。

图1 雷电回击过程示意图Fig. 1 Geometry of lightning return stroke

根据 Wait算法, 雷电回击过程产生的垂直电场沿起伏海面传播时, 距雷电回击通道水平距离r处的垂直电场为[14]:

为了研究Wait近似算法的准确性, Shoory等人[24]分析了混合传输路径雷电电磁传播的 Wait算法, 指出在距离闪电通道10 km时, Wait近似算法的计算得到的峰值电场最大观测误差为 4.8%, 峰值上升时间的最大误差为 18%。张其林等人[25]又进一步验证了距离闪电通道200～1 000 m范围内Wait近似算法的精确度, 认为在距离闪电通道几百米的近距离范围内, Wait近似算法的精确度也可以接受。

按照Barrick等效表面阻抗理论[26-27], 起伏海面的等效表面阻抗Δeff可表示为:

其中, Δ为平静海面的等效表面阻抗, Δ′为起伏海面引起的表面阻抗的增量。

其中, γ和η分别为x和y方向的波数(或空间频数)。V(γ, η)是不规则起伏地表的高度谱密度函数。如果已知 V(γ, η), 根据公式(6)—(12)可以计算起伏海面的归一化等效表面阻抗Δeff。

1.2 模拟起伏海面的模型

在模拟海面的雷电电磁传播过程中, 选取合适的海面模型非常关键, 近年来人们发现海浪既具有规律性又具有随机性, 利用二维分形模型能够更为贴切地描述真实海面。本文采用王运华[28]提出的一种改进二维分形海面模型, 这种改进模型同时考虑了空间波数大于和小于基波波数的两种波浪对海面表面轮廓的影响, 海面谱的表达式为:

其中 S(γ, η)为谱密度, D(γ, η, φ)为方向分布函数。

其中φ为海浪和风向之间的角度,β0为雷电电磁场和风向之间的角度。

2 个例分析

2.1 长门岩浮标站风速数据的应用

本文选取2015年4月2日这次雷暴过程, 将青岛长门岩浮标站的风速数据输入到以上的模式中,分析此次雷暴过程青岛近海起伏海浪对雷电电磁场的影响。

青岛长门岩浮标站位于距青岛海岸线 23.3 km的长门岩岛。这里采用的是 10分钟的风速平均值,也就是十分风速。图2给出了2015年4月2日雷暴天气过程时的十分风速, 可以看出, 此次雷暴过程不同时刻风速的变化非常大, 十分风速的最大值达到13.4 m/s, 最小值为2.5 m/s。

2.2 雷电回击通道基电流的选取

回击通道底部电流即基电流, 基电流中的击穿电流和电晕电流都采用双Heilder指数表达式计算[29],其表达式如下

其中,I01、I02分别表示击穿电流和电晕电流最大值,η1、η2为修正因子,τ11、τ12、τ21、τ22分别表示用来确定电流上升时间、衰减时间以及峰值的时间常量,n1、n2为电流陡度因子。这些参数具体取值见表1, 则闪电回击通道底部基电流的波形如图3所示。

表1 继后回击各参数的取值[30]Tab. 1 Parameters of the subsequent return strokes

图3 闪电通道底部采用的基电流波形Fig. 3 Channel-based current waveform corresponding to a typical subsequent return stroke

闪电回击模型选取MTLL(the modified transmissionline model with linear current decay with height)模型[31],即回击电流波形随高度呈线性衰减, 其表达式为:

2.3 海浪对雷电电磁场传播的影响

图 4给出了起伏海面对地闪回击垂直电场的影响,观测点距闪电通道的水平观测距离d分别为1、30、100、200 km。假定电磁场传播方向与风向垂直。曲线1表示电导率无限大的理想海面, 曲线2、3、4、5分别表示风速为2.5、5.2、8.7、13.4 m/s的海面垂直电场强度。

图4 起伏海面对时域垂直电场的影响Fig. 4 Propagation effect of sea surface on the lightning vertical electric field in the time domain at distances of 1 km (a, b), 30 km(c, d), 100 km (e, f), and 200 km (g, h) from the lightning channel

由图4可以看出, 观测距离相同时, 海面的起伏对垂直电场峰值的影响不显著, 但对波形的上升时间的影响较为明显。风速越大, 即海面的起伏程度越大, 上升时间越长; 当观测距离变化时, 垂直电场达到半峰值的时间随着距离而增大。

为了具体分析风速对上升沿时间影响的原因,图5给出了图4对应的垂直电场频谱。可以看出, 随着观测距离的增加, 地闪回击电磁波谱的高频成分优先衰减, 观测距离小于100 km时, 频率超过10 MHz以上的部分按照 1/f2的速度快速衰减(f为频率)。超过100 km时, 电磁波谱的中低频段也逐渐开始衰减。

同样的, 我们对水平磁场进行分析研究, 发现海面的起伏对水平磁场的影响与垂直电场类似。差别在于, 二者到达半峰值的时间不同, 相同的海表情况和相同的观测距离, 水平磁场到达半峰值的时间滞后要短一些。

因此, 当观测距离很远且风速较大时, 无论是利用垂直电场还是水平磁场, 进行雷电放电参数的反演, 海面起伏的影响都不容忽视。

2.4 海浪对闪电定位精度的影响

根据以上对雷电电磁波沿起伏海面传播特性的研究揭示, 海面的起伏会导致地闪回击电磁波上升沿时间的有不同程度的滞后, 这会引起基于时差法的闪电定位系统在定位时产生误差, 从而影响闪电定位的精度。

2.4.1 山东省闪电定位系统介绍

图5 起伏海面对垂直电场频谱的影响Fig. 5 Propagation characteristics of sea surface on the lightning vertical electric field in the frequency domain at distances of 1 km (a, b), 30 km (c, d), 100 km (e, f), and 200 km (g, h) from the lightning channel

山东省闪电定位系统于2006年6月布设完成, 由中国华云技术开发公司研制生产。系统主要由采用了时差法(TOA)进行闪电定位, 由 13个闪电定位仪组成, 分别在章丘、龙口、荣成、即墨、日照、东明、东平、沾化、夏津、鱼台、蒙阴、郯城和昌邑, 位置如图 6所示。山东省内大部分地区闪电探测效率理论值为 95%(图 6)。

图6 山东省LD-II闪电定位系统探测效率图Fig. 6 Detection efficiency of LLS in Shandong Province

2.4.2 计算结果及分析

对于时差法的闪电定位系统, 为了减小下垫面的影响, Cooray、Honma等人[32-33]提出了几种确定与测量波形相关的到达时间的方法。在本文中, 我们将考虑两种方法: (1)对应于电场峰值的波头时间, (2)对应于电场的半峰值时间。

表2和表4分别给出了此次天气过程中不同风速下电场峰值对应的波头时间和半峰值时间。表3和表5分别对应于表2和表4的延迟时间。值得注意的是,当观测距离较近时, 风速对电场峰值的波头时间影响很小。电场的延迟时间随风速以及观测距离的增大而增大。当距离可达100 km以上时, 风速对峰值场的时间有显著影响。当风速达到13.4 m/s, 观测距离为200 km时, 波头时间为0.575 μs, 对应的延迟时间为0.075 μs, 半峰值时间为0.35 μs, 对应的延迟时间为 0.165 μs。

表2 电场峰值的波头时间(单位: μs)Tab. 2 Time corresponding to the peak lightning electric field

表3 表2对应的不同风速下起伏海浪相对于平静海面(v=0 m/s)时的延迟时间(单位: μs)Tab. 3 Time delay corresponding to Table 2

表4 电场的半峰值时间(单位: μs)Tab. 4 Time corresponding to 50% of the peak lightning electric field

表5 表4对应的不同风速下起伏海浪相对于平静海面(v=0 m/s)时的延迟时间(单位: μs)Tab. 5 Time delay corresponding to Table 4

为了进一步讨论海面起伏对时差法闪电定位探测精度的影响, 我们选取了2015年4月2日天气过程中发生的两次闪电(图 7)。这两次雷击共利用了四个测站进行定位, 分别为即墨站(36.392°N, 120.469°E)、日 照 站 (35.432°N, 119.539°E)、 荣 成 站 (37.17°N,122.493°E)和龙口站(37.64°N, 120.341°E)。表 6 给出了输入参数和结果。计算得出, 当风速达到13.4 km/s时, 海面起伏对闪电定位精度有显著影响, 由海面的起伏引起的定位误差分别为12 662 m和3 944 m。第一次闪电定位误差较大的原因是由于第一次闪击点更靠近闪电定位系统组网探测的外围区域, 产生了较大的布站误差, 该误差是由站点之间布放的位置形状产生的。

图7 海面起伏对闪电定位的影响Fig. 7 Effect of rough sea surface on lightning location

表6 图7对应的计算参数和结果Tab. 6 Input parameters and results corresponding to Figure 7

3 结论

本文分析了海面起伏对雷电垂直电场的影响,并选取2015年4月2日雷暴天气过程, 将青岛浮标站测量的风速数据输入计算模型。结果表明, 海面起伏对垂直电场的峰值影响较小, 但会引起波形的上升时间的延迟, 风速越大, 即海面的起伏程度越大,影响越明显。随着观测距离的增加, 雷电垂直电场波形的上升时间逐渐变长。分析垂直电场频谱发现, 随着观测距离的增加, 地闪回击电磁波谱的高频成分优先衰减, 观测距离小于100 km时, 频率超过10 MHz以上的部分快速衰减。当观测距离超过 100 km时,电磁波谱也逐渐由中低频段开始衰减。

雷电电磁场沿起伏海面引起的上升时间的延迟,将会影响基于时差法闪电定位系统的定位精度, 通过个例验证, 定位误差可达几至十几公里。由于闪击点位置的不同, 定位误差的差异可能很大。这其中也包括了闪电定位仪组网形状产生的系统误差。