1. 引言
雷电是自然界中的一种长距离瞬时放电现象,也是大气科学的重点研究对象。长期以来,伴随闪电探测手段的提高,相关研究人员对雷暴云电特性以及雷电临近预报方面做了大量深入的研究工作并取得了丰硕的研究成果。Murphy等 [1] 利用大气电场仪、雷达等多种探测资料开发了雷电预警系统。Harnwell [2]、罗福山、张腾飞等人 [3] [4] 对电场仪的标定方面做了进一步的改进工作。张义军等通过对青藏高原东部地区大气电特征的分析发现,雷暴云电结构的差异所呈现出的不同大气电场特征 [5]。吴健等认为大气电场预警指标应以电场强度的阈值、以及电场特征曲线有无快变抖动为预判指标 [6]。王凯、朱浩等统计分析了黄山风景区雷暴过程近地表静电场特征,并建立本地雷电预报方程进而发现本地雷暴地面电场特征主要以负地闪为导向 [7]。相关文献 [8] [9] [10] [11] [12] 将大气电场数据等多种资料融合来提高雷电监测和预警能力。
深入了解大气电场的变化特征对于本地雷电的认识以及雷电警报阈值的确定至关重要,此外也有助于不同强对流天气的过程分析以及气候特征的研究。本文利用2016~2017年南京本地各站夏季雷暴大气电场资料并结合闪电定位、雷达资料重点分析雷暴过境时大气电场特征,在此基础上根据首次地闪发生时刻分析南京本地化各项预警因子阈值,以期尝试探索大气电场资料在改进雷电预报预警中的应用前景。
2. 资料来源
1) 本文使用的数据来源于南京地区大气电场仪(Pre-storm2.0型)观测。在对数据进行质量控制的基础上,选取了南京市区及周各站观测数据资料,其大气电场仪位置分布见图1。
Figure 1. Distribution of atmospheric electric field instrument stations
图1. 大气电场仪站点分布
本文地面电场特征的研究范围定义为31.3328˚N~32.0681˚N、118.5947˚E~119. 0343˚E (见图1)。大气电场仪(Pre-storm2.0型)的性能参数见表1。
Table 1. Main performance index of 1pre-storm2.0 Atmospheric Electric Field Instrument
表1. Pre-storm2.0型大气电场仪主要性能指标
2) 本文所使用的闪电定位资料来源于江苏省气象局闪电定位仪。
3) 本文雷达数据来自于南京多普勒天气雷达站,其探测范围为230 km。
3. 南京地区雷暴天气大气电特征分析
大气电场是大气电学的一个基本的参数,通过对地面大气电场变化特征的研究有助于强对流天气发生时的雷暴监测预警。通过对南京各测站2016~2017年夏季积累的大气电场资料的观测分析,发现南京本地在发生单体雷暴过程时,地面电场主要有以下几个变化特征:
1) 正向增大型:表明此时电场已经发生了瞬间倒转,一般发生在先导闪电发生后,如图2(a)江宁观测站一次雷暴过境时的电场变化,此时雷暴云底部和地面之间的电势差已经达到了击穿场强。
2) 负向增大型:图2(b)所示为江心洲站一次雷暴过程时地面电场负向增大型变化,当对流云团逐渐向测站靠近时,伴随雷暴云底部聚集的负电荷逐步累积,大气电场强度随即负向增加,如图中“T”所示的区域部分。
3) 正负交替型:当雷暴云已经进入电场仪预警范围内并发生对地放电时,主要表现为图2(c)所示的昆仑路测站的正负交替型电场特征。
此外多单体雷暴过境时的地面电场变化主要是由多个单体雷暴在经过不同的生消发展阶段并且对地放电时电场极性的频繁快速反转进而呈现出的综合效果 [13] [14] [15]。如图2(d)所示为六合观测站探测到的一次多单体雷暴发生时的电场变化特征曲线图。
(a) (b) (c) (d)
Figure 2. (a) A thunderstorm at Jiangning station is increasing the magnitude of the electric field; (b) negative increasing electric field during a thunderstorm at Jiangxinzhou station; (c) an alternating electric field during a thunderstorm at Kunlun Road Station (d) Luhe Station Multicellular thunderstorm process electric field
图2. (a) 江宁站一次雷暴过程正向增大型电场;(b) 江心洲站一次雷暴过境时负向增大型电场;(c) 昆仑路站一次雷暴过程正负交替型电场;(d) 六合站多单体雷暴过程电场
4. 雷暴过境时地面电场变化特征分析
4.1. 雷暴过程电场变化的阶段特征分析
一般而言,大气电场仪的有效探测范围为以站点为圆心,半径为15 km以内雷暴云电荷量变化情况,故而可以将距测站15 km范围内的区域叫做预警区域(Row),在对2016~2017年6~8月南京地区140个雷暴天气电场样本详细研究分析的基础上,根据发生闪电时的位置距电场测站的距离得出雷暴过程中的大气电场变化大致有如下阶段:
第一阶段:无雷暴云影响或者闪电发生在距电场仪较远处的阶段,此时电场的快变抖动主要围绕在零点附近,相应的电场幅值变化也不太明显;
第二阶段:雷暴云逐渐向测站靠近的阶段,闪电在15 km范围内发生,伴随云底聚集的电荷累积电场幅值渐渐增大,当先导闪电发生时,电场迅疾发生极性翻转;而后伴随雷暴云更进一步的临近,电场曲线表现为频繁上下跳动;
第三阶段:雷暴云继续发展移近到电场仪5 Km范围内,电场曲线变化抖动剧烈;
第四阶段:雷暴云离境阶段,这个阶段的地面电场经过一段时间震荡后,电场幅值也渐渐恢复到晴天时的电场值,雷暴过程进入消亡阶段。
以2016年7月15日南京江宁测站的一次雷暴过程电场资料为例,设闪电定位仪探测到某电场仪站点15 km内首次闪击的时刻为t1,探测到该站点10 km内首次闪击的时刻为t2,探测到该站点5 km内首次闪击的时刻为t3。则上述各个阶段的电场曲线变化特征见图3(a)。其电场幅值的变化情况也与雷暴云的逐渐临近密切相关,见图3(b)和图3(c)。由图3(a)可以看出15:20分之前雷暴云逐渐移近江宁观测站,随后电场反向增大,主要是受雷暴云所带正电荷所致;当15:10分左右雷暴云抵达测站附近区域时,此时大气发生了明显的击穿放电现象(首次闪击发生在15:12分),主要表现为图3(a)中较为明显的针状脉冲凸起,伴随第一次充放电的完成,从图中可以看出后续也发生了多次闪击,主要为负地闪,电场也主要呈现出明显的负向增大特征。在此过程中与之相对应的雷达回波的强度也是明显加强的,图3(b)所示的为此次雷暴过程临近测站时的回波变化发展情况,伴随闪电的临近,回波面积逐步增大,回波强度加强。之后随着雷暴云的过境与离开,电场曲线波动也会逐渐回到0轴附近,在15:48分以后(见图3(a)是雷暴云的消散阶段,从图3(c)上也能看出此阶段的雷达回波也是逐渐减弱并消散的。
(a) (b) (c)
Figure 3. Electric field variation and radar echo during a thunderstorm at Jiangning station on 15 July 2016; (a) the electric field curve and the distance from the flash point to stage T of a thunderstorm at Jiangning Station on 15 July 2016 are descending waveforms; (b) radar echo at 15:25 minutes within 10 km of the second stage of the station; (c) radar echo at 16:02 hours during phase IV thunderstorm cloud dispersal phase
图3. 2016年7月15日江宁站一次雷暴过程电场变化及雷达回波图;(a) 2016年7月15日江宁站雷暴过程电场曲线和闪击点距离T阶段为下降型波形;(b) 测站第二阶段10 km内15:25分雷达回波图;(c) 第四阶段雷暴云消散阶段16:02分雷达回波图
根据以电场仪测站为中心半径15 km的区域叫做预警区域(Row)定义的雷暴过程的四个阶段(见图3(a)中t1-t4阶段),当发生闪击时距江宁电场测站的各阶段的具体闪电情况分析如下:
第一阶段(t1-t2时段):闪电定位仪探测到的电场仪站点15 km内首次闪击发生在15:12 (北京时,下同),距离江宁站13.1 km,这一阶段共探测到42次闪击,距测站最近距离为9.6 km。
第二阶段(t2-t3时段):10~5 km内首次闪击发生在15:21,此时进入雷云临近阶段。
第三阶段(t3-t4时段):15:30分闪电定位仪探测距测站4.4 km处闪击,此时进入雷云到达阶段,一直持续到15:48,该时刻探测到5 km范围内最后一次闪击。
第四阶段(t4时刻之后数10 min):雷暴云消散阶段,电场曲线经过30 min左右震荡后渐趋于平稳。
4.2. 大气电场仪预警阶段特征分析
4.2.1. 电场幅值变化特征
以2016~2017年6~8月南京周边地区各站点的电场观测资料为样本(其中雷暴样本共计140个),结合所有夏季雷暴过程中闪电定位仪定位到的以南京各测站为中心,半径分别为15 km、10 km、5 km的首次闪击数据,对南京八个测站的电场仪报警阈值以及阈值的分钟变化率、预警时间、阈值持续时间、T阶段时间进行了统计分析,在此设电场强度绝对值到达预警阈值开始预警时刻为T,则对应到达一级预警阈值的时刻为T1,到达二级预警阈值的时刻为T2,到达三级预警阈值的时刻为T3;电场仪在T1、T2、T3时刻达到的预警阈值分别为E1、E2、E3,与之相对应的到达报警时刻1分钟内平均阈值分别为E1/T1,E2/T2,E3/T3。
设实际地闪发生的时刻为t,则距大气电场仪15 km、10 km、5 km半径各区域内发生第一次地闪的时刻分别为t1、t2、t3; 以上述标准参数的设定为基础,观察分析得到南京各站大气电场仪预警参数分析结果见表2:
Table 2. Distribution of early warning thresholds of electric field meters in Nanjing area
表2. 南京地区各测站电场仪预警阈值分布
从表2中可以看出,总体上南京各站在15 km、10 km、5 km范围内达到的平均阈值在1.33 kv/m、2.23 kv/m、3.06 kv/m左右;与之相对应的1分钟平均阈值分别为1.25 kv/m、1.93 kv/m、2.51 kv/m左右。其中,奥体站、高淳站、六合站达到的报警阈值相对较小。
大气电场的幅值大小可以作为雷电预警研究中的指标,相关研究人员 [16] 也根据各地发生雷暴天气时大气电场仪预警时所达到的幅值来定义当地雷电预警时的电场阈值。然而现实中以某一固定电场阈值统一当作一个地区所有站点的报警门限值势必会影响电场仪在使用过程中的预警效果。因此可以引入电场放大倍数的概念,即雷暴发生时,电场幅值变化的平均值与非雷暴天气时大气电场幅值的比值(一般默认为晴天大气电场幅值),可以在一定程度上弱化各测站因局地性差异带来的影响,进而有效反映各个站点在雷暴过程中电场幅值变化的特征,其定义见公式(1) [17]:
(1)
其中:ei为某站点雷暴天气中电场仪预警时电场幅值变化的平均值;Ei为该站晴天大气电场幅值的平均值。计算结果表明,140个雷暴样本中,电场放大倍数m的最小值为0.6,最大值为21.3,分别出现在高淳站和浦口站的雷暴过程。根据相关理论对晴天的定义,即不发生降水,且当天云量少于3/8,平均风速小于4 m/s,统计分析得到近三年南京各站点晴天大气电场幅值的平均值分布(见图4)。从中也可以看出南京各站晴天大气电场分布在0.106~0.388 kv/m之间。
Figure 4. The distribution of the average atmospheric electric field volt value in KV/M during the summer of 2016~2017 around Nanjing
图4. 2016~2017年夏季南京周边地区晴天平均大气电场伏值分布(单位:kv/m)
4.2.2. 预警时间特征
由于大气电场仪的有效探测范围为15 km以内雷暴云电荷量变化情况,故而可以将距测站15 km范围内的区域叫做预警区域(Row),可以将距测站15 km范围内的预警时间定义为T1,距测站10 km范围内的预警时间定义为T2,距测站5 km范围内的预警时间定义为T3,依据2016~2017年的雷暴过程样本计算分析可得南京各测站大气电场仪预警时间分布见表3。
Table 3. Distribution of early warning time of electric field meter in Nanjing area
表3. 南京地区各测站电场仪预警时间分布
由图5可以看出南京各测站15 km范围内的预警时间发生在23~44分钟之间,10 km范围内的预警时间在19~52分钟,5公里范围内的预警时间在13~53分钟,阈值持续时间大致在10~49分钟左右,T阶段时间在15~47分钟。
Figure 5. Distribution of early warning time of stations in Nanjing
图5. 南京各测站预警时间分布
5. 各站大气电场仪预警效果分析
对2016~2017年6~8月南京观测区多个地面大气电场测站(南京本站、昆仑路站、奥体站、江心洲站等)和江苏省气象局闪电定位系统以及新一代天气雷达等资料的综合分析来进行预警研究。同时为了提高大气电场仪的探测效率,重点对以电场仪为中心,半径为15 km的范围(预警范围(ROW))进行分析。 设电场强度绝对值到达预警阈值开始预警时刻为T,则对应到达一级预警阈值的时刻为T1;到达二级预警阈值的时刻为T2。设实际地闪发生的时刻为t,则距大气电场仪15 km、10 km、5 km半径各区域内发生首次地闪的时刻分别为t1、t2、t3。T − t为预警时间,若T − t > 0,可以认为预警成功,以X表示;若T − t < 0,则预警未成功,可以视为漏报,以Y事件表示;若T时刻后相关区域内并未发生闪电则可以视为空报,以Z事件表示。
通常用预报评分方法 [18] 中的探测概率、漏报率、虚警率以及临界成功指数等参数来评判预报效果,其具体计算方法如下:
(2)
(3)
(4)
(5)
在以上公式为理论基础的前提下,结合闪电定位仪以及大气电场仪数据对南京地区2016~2017年各站夏季雷暴过程的雷电预警效果进行了深入的分析研究,进而得到南京八个测站大气电场仪的预警探测效率(POD)、漏报率(PO)、空报率(FAR)、临界成功指数(CSI)的结果见图6(a)~(d)。
(a) (b) (c) (d)
Figure 6. Distribution of early warning detection efficiency (POD), false alarm rate (PO), false alarm rate (FAR) and critical success index (CSI) of atmospheric electric field instruments at various stations in Nanjing; (a) The early-warning detection efficiency of Atmospheric Electric Field Instrument in Nanjing Stations (POD); (b) The rate of missing report of Atmospheric Electric Field Instrument in Nanjing Stations (PO); (c) Air-to-air reporting rate of atmospheric electric field instruments at various stations in Nanjing (FAR); (d) Critical Success Index of atmospheric electric field instruments at various stations in Nanjing (CSI)
图6. 南京各测站大气电场仪的预警探测效率(POD)、漏报率(PO)、空报率(FAR)、临界成功指数(CSI)分布情况;(a) 南京各测站大气电场仪的预警探测效率(POD);(b) 南京各测站大气电场仪的漏报率(PO);(c) 南京各测站大气电场仪的空报率(FAR);(d) 南京各测站大气电场仪的临界成功指数(CSI)
从中可以看出:1) 南京本站、高淳站大气电场仪预警准确率总体上在15~10 km范围效果较好,其余各站基本上预警准确率随着预警半径范围的缩小,预警率也是增加的;
2) 高淳站与南京本站的大气电场仪在使用过程中发现其漏报率较高,分析其原因与这些站点的数据缺失有一定的关系;
3) 各站临界成功指数总体上在15~10 km范围内的探测效率较高。
6. 结论与讨论
1) 本文在结合闪电定位、雷达与大气电场仪数据的基础上对南京地区2016~2017年各站夏季雷暴过程的大气电场特征进行了深入的分析研究,得出电场变化与闪电发生时刻之间的关系;
2) 以各站夏季雷暴过程的电场观测资料为样本(其中雷暴样本共计140个),结合闪电定位仪定位距各测站不同范围内的首次闪击数据,对南京本地八个测站的电场仪报警阈值以及阈值的分钟变化率、预警时间、阈值持续时间、T阶段时间进行了统计分析,得出南京本地电场仪预警阶段的各项预警参数值;
3) 在此基础上对南京各测站大气电场仪的雷电预警效果进行了计算分析,发现南京本站、高淳站大气电场仪预警准确率总体上在10~15 km范围效果较好,其余各站基本上预警准确率随着预警半径范围的缩小,预警准确率增加;各站临界成功指数总体上在15 km范围内的探测效率较好;
4) 本文以2016~2017年夏季雷暴过程电场为样本资料,分析得出南京本地化预警参数值并探讨了南京本地各站的电场仪预警效果,由于目前积累的后续雷暴过程资料尚不丰富,所以未能从预报方程的角度去检验各站预警参数值并予以修正;接下来的工作中,将试图融入雷达预警参数因子等建立南京本地化的雷电预报方程,以便从更加精细化的角度来检验并修正本文得出的各项预警参数值。
基金项目
华东区域气象科技协同创新基金合作项目(QYHZ201604),江苏省气象局面上项目“多源观测资料在雷暴云特征分析和预报技术中的应用研究”(KM201909)。