钱 勇,施俊杰,王延慧,刘成武,张小军,张建涛
(1.新疆气象灾害防御技术中心,新疆 乌鲁木齐830002;2.新疆气象台,新疆 乌鲁木齐830002;3.国网新疆电力有限公司电力科学研究院,新疆 乌鲁木齐830011;4.中国气象局乌鲁木齐沙漠气象研究所,新疆 乌鲁木齐830002)
闪电是指发生在不同极性电荷之间的长距离放电现象,所伴随的瞬时高温和强大的电磁辐射具有强大破坏力,是我国危害程度仅次于暴雨洪涝、气象地质灾害的第三大气象灾害[1]。随着科学技术的发展,对闪电的监测技术也不断地丰富,所得闪电数据资料日趋丰富。根据传感器所处位置不同,闪电观测资料获取主要分为两种方式,基于闪电辐射电磁场传播理论的地面闪电探测技术和基于光学成像原理的卫星闪电探测技术,两者各有优势,互为补充[2-4]。
我国地面闪电观测系统目前主要有中国气象局建设的国家级闪电定位系统(Active Directory Topology Diagrammer,ADTD)、美国的全球闪电定位网 (World-Wide Lightning Location Network,WWLLN)、电力部门闪电定位系统(Lightning Location System,LLS)、各省气象局自行建设的甚低频/甚高频三维闪电定位系统、大气电场监测系统等。根据卫星所处轨道高度的不同,卫星闪电探测分为极轨卫星探测和静止卫星探测,前者以热带测雨卫星TRMM(Tropical Rainfall Measuring Mission)携带的LIS(Lightning Imaging Sensor)为代表[5],后者以中国风云四号卫星(FY-4A)LMI(Lightning Mapping Imager)和美国的GOES-R GLM(Geostationary Lightning Mapper)为代表[6-7]。极轨卫星闪电成像仪受观测平台制约,对任一地点的观测时间十分有限,只能记录该地点全部闪电事件中的小部分,提供区域平均意义上的闪电分布信息。静止卫星闪电成像仪可连续、实时监测闪电活动,可实现对强对流天气系统的追踪和预警,是卫星闪电探测新的发展方向,但目前全球都处于研究阶段[8]。随着闪电探测技术不断发展与完善,积累了丰富的数据资料,国内外学者根据不同类型的数据资料,开展了雷电预警的研究、区域雷电活动特征的分析等。王娟等[9]统计分析了中国国家闪电监测网数据资料,给出了地闪的活动规律。Shvets、郄秀书、朱润鹏等[10-12]利用LIS资料分析了全球闪电分布情况。柴瑞等[13]分析了地面大气电场数据在雷电预警中的应用。曾庆锋等[14]对大气电场仪观测数据进行了质量控制和统一标定,并结合天气雷达、闪电定位仪和探空等多源数据资料,对地面大气电场的时空演变特征和电荷结构进行了分析,提出了利用单一来源观测资料及综合利用多源资料进行雷电临近预警,并评估了预警效果。
随着地面和卫星闪电监测技术的发展,积累了大量不同角度的数据资料,为闪电活动特征研究提供了基础。但目前的研究多是基于单一方式获取的闪电观测数据,孤立地分析某一区域雷电活动特征,其结果具有一定的局限性。不同地区由于闪电监测站点布站方式、站点密度、卫星观测角度的不同,闪电资料的可利用性也不同,因此分析某一区域闪电活动特征的前提是要弄清楚该区域各种闪电数据资料的可利用性,进而根据各种数据的特点,融合多源数据更加准确地进行分析。
新疆目前可利用的星地闪电观测资料有ADTD、LLS、WWLLN、LMIG、大气电场。新疆ADTD、LLS、大气电场仪分布如图1。本文主要通过上述多源闪电观测资料,根据2019年7月20日新疆大部分区域出现局部强对流天气过程,分析这5种资料在新疆地区的应用情况,为后续融合多源数据更加准确地分析新疆地区闪电特征、闪电资料并在雷电预警中的应用提供基础。
图1 新疆ADTD、LLS、大气电场仪站点分布
ADTD闪电监测系统由中国科学院空间科学与应用研究中心研制,利用到达时差法的原理对地闪进行定位,有效探测半径是150 km,可以给出闪电发生的时间、电流强度、极性、经纬度、陡度等相关参数[15]。新疆目前已建设49个监测站点,但由于地域辽阔仍然不能实现有效的覆盖,除北疆部分区域外,新疆大部分区域探测效率不理想。
新疆电力部门闪电定位系统LLS于2012年7月建设完成,目前共建设54个闪电定位探测站,单站探测范围约为100 km,定位误差≤1 km。系统定位算法采用时差法、定向定位法和定向时差联合法对地闪定位[16]。数据包括地闪发生时间、经纬度、电流、回击数、定位站数等要素。王延慧等[17]根据近5 a的数据资料对比分析了新疆ADTD和LLS两套闪电探测系统,得出了ADTD系统探测的地闪次数偏少、地闪密度偏小,而LLS系统监测的地闪强度偏小、地闪强度高值区分布较为分散。
全球闪电定位系统WWLLN是由美国研究和开发的一套旨在全球范围实时连续监测闪电活动的定位系统,该系统利用TOGA(Time of Group Arrival)进行闪电定位,WWLLN以监测强闪电为主,数据中大部分为强地闪,少部分为强云闪,能给出闪电发生的时间、经纬度。由于测站监测的是闪电辐射甚低频(VLF)电磁信号,该波段的电磁信号能以低衰减速度稳定传播,以至于测站间的分布距离可以达到上千千米。目前亚洲地区有6个测站,分别位于中国、日本、新加波境内。我国的3个测站分别位于南京、北京、兰州。WWLLN自投入商业化运营以来,其全球测站数目逐年增加,这也使得WWLLN能够更加有效地对全球范围内闪电活动进行连续观测[18]。新疆目前很少有人采用全球闪电定位系统的数据资料。
风云四号静止卫星(FY-4A)携带的闪电成像仪(Lightning Mapping Imager,LMI)星下点104.7°E,CCD面阵大小为400×600,空间分辨率为7.8 km。当闪电成像仪CCD面阵单个像元探测闪电光辐射亮度高于背景阈值时,定义为一次闪电“事件”,即像元观测区域云中闪电光辐射透过云层的发光现象,对应的CCD面阵像元中心位置即为这次闪电“事件”的位置[19-20]。同一帧CCD图像上的多个相邻微小像元探测到的闪电“事件”通过聚类算法组成一个“组”,对应于地闪的一次回击或云闪的一次K变化,满足一定阈值条件的多个“组”聚类为一次真实的“闪电”过程。数据资料来源于国家卫星中心网站提供LMI的L2级“事件”和“组”产品,“闪电”产品暂未公开,因此本文采用LMI的L2级“组”产品(LMIG)进行分析。
大气电场仪是用来测量大气电场及其变化的设备,有效探测半径约为15 km。它是利用导体在电场中产生感应电荷的原理来测量电场,雷暴地面电场是云内所有电荷共同作用的结果,当雷暴云移动时,可引起电场的缓慢低频变化,这种变化可实现对云内的电荷区进行评估,当雷电发生时,可引起电场的快速高频变化,根据电场变化极性,可大致判断雷电发生的区域。目前,国内部分省市已开展了大气电场仪与闪电定位仪结合的雷电监测预警技术的研究。秦微等[21]利用大气电场资料和闪电资料,根据大气电场的波形特征,提取与闪电相关性较高的因子,利用多元回归技术建立一套预报方程,根据该预报方程得出了最佳预警参数。谢屹然等[22]根据雷暴云近地面电场变化特征,探讨了大气电场仪临近预警方法,首次提出极性变化预警法,其命中率和提前预警时间都明显优于阈值预警法。目前乌鲁木齐市有12个观测站点,克拉玛依市有22站点已投入运行,电场分辨率<5 V/m,响应时间为1 s。
2019年7月20日全疆大部分区域出现了强对流天气过程,本文选取了这一天的闪电资料进行分析。图2分别是ADTD地闪监测产品、LMIG全闪监测产品、WWLLN全闪监测产品、LLS地闪监测产品、大气电场监测产品(地基闪电定位数据没有进行归闪处理)。由图2中可知,LLS的探测效率最高,其次是ADTD,LMIG和WWLLN虽属于全闪监测系统,但在北疆大部分区域探测的总闪频数要低于ADTD和LLS探测的地闪频数,说明LMIG和WWLLN对地闪的探测效率较低。而在南疆部分区域LMIG探测的总闪频数高于ADTD、低于LLS,相比于LLS这主要是由于南疆部分地区ADTD监测系统不完善,探测效率较低,这也充分说明地面闪电监测系统对站点的布设较为依赖。此外,LMIG监测的高密度区与ADTD、WWLLN、LLS监测结果存在一定的差异,需要作进一步分析。图2e是一次天气过程的地面大气电场变化过程,在起始阶段和结尾阶段波形快变抖动频率较中间阶段慢,在结尾阶段表现的较为明显。
图2 多源闪电数据监测产品
全天ADTD监测闪电频次3768次,LLS监测18327次,WWLLN监测786次,LMIG监测1206次。由于各观测系统监测闪电频数处于不同的量级,不便于直观的对比各个时段闪电频次的分布。因此,本文将闪电频次归一化处理,图3是归一化后ADTD、LLS、WWLLN、LMIG各时段的闪电频次分布。ADTD、LLS、WWLLN监测闪电发生时段主要集中在12—23时,而LMIG主要集中在16—23时。ADTD、LLS、WWLLN各时段闪电频次变化基本一致,16时是闪电高发时段,而LMIG监测闪电高发时段是22时。
图3 归一化后ADTD、LLS、WWLLN、LMIG各时段闪电频次分布
2019年7月20日乌鲁木齐市发生一次局部强对流天气过程。图4是19—23时(北京时间,下同)乌鲁木齐雷达站组合反射率部分图(数据经过质量控制)。19时对流云团由石河子经呼图壁、昌吉向乌鲁木齐移动,21时在乌鲁木齐市天山区附近对流云团发展至旺盛阶段,22时对流云团逐渐减弱。此次天气过程,属于局部小范围天气过程,降雨量不大,但雷暴强度是近几年罕见的,雷达监测到了此次天气过程基本完整的发展路径。
图4 2019年7月20日19—23时乌鲁木齐雷达站组合反射率
图5分别是ADTD、LLS、WWLLN、LMIG监测的北京时间18—00时以乌鲁木齐雷达站为中心周围150 km的闪电分布情况。沿着此次强对流天气移动路径上ADTD、LLS、WWLLN三者在时间和空间上基本对应一致,仅是探测效率有所差别(路径以外局部的闪电数据暂不考虑,由于雷达被地物遮挡或雷达监测的此次天气过程不全,无法参考雷达从时间和空间上进行比对)。而LMIG监测的闪电主要发生时段是22—00时,高密度区位于达坂城区域,且在乌鲁木齐天山区附近没有监测到闪电,监测结果与前三者在时间和空间上都存在差异。通过查询22—00时雷达组合发射率和基本反射率情况,发现在该区域没有对流天气过程,因此初步认为这可能是由于LMIG数据出现了偏差。
图5 北京时间18—00时以乌鲁木齐雷达站为中心周围150 km的闪电分布
考虑到雷达可能会受到地物的遮挡,进一步分析了风云四号卫星云顶亮温资料。卫星携带的多通道扫描成像辐射计AGRI可获取地球表面和云的多光谱、高精度定量观测数据和图像,采用双扫描镜结合三反射光学系统、线列阵探测器,可以全天空、实时监测地面的强对流天气过程,通过红外探测技术,能够感知对流云顶的亮温。AGRI分为两种扫描模式:全圆盘扫描(扫描时间15 min),中国区域扫描(扫描时间为4 min)。两种扫描方式交替进行。云顶亮温可用来表征对流活动,亮温越低,表明云顶越高,云层越厚,对流越旺盛[23],闪电通常是发生在强对流区域。图6是22:00—00:00乌鲁木齐地区云顶亮温部分图,与LMIG监测结果相对应区域的云顶气温分布在270~290 K(虚线圆圈区域),对应的摄氏温度为-4~17℃,根据每上升1 km温度下降6℃,可以判断该时段内此区域没有强对流云,风云四号卫星监测的云顶越高更可能是地表温度。进一步论证了LMIG数据存在一定偏差,这可能是由于新疆处于风云四号卫星监测的边缘,距离星下点较远,像元分辨率较大,导致出现偏差。因此,使用LMI数据前需要利用相关资料进行数据检验。
图6 22—00时乌鲁木齐地区部分云顶亮温
大气电场仪可以直观地看出监测区域内电场强度的变化情况,能够记录闪电发生前后雷暴云中闪电活动情况。图7是天山区气象站内的大气电场仪监测的此次天气过程地面电场强度和闪电的分布情况(以大气电场站点为中心,ADTD、LLS、WWLLN分别监测的周围40 km闪电分布)。当大气电场仪周围15 km有闪电发生时,电场强度值会出现快变抖动现象,随着闪电发生位置越来越近快变抖动会越频繁,而随着闪电发生位置越来越远时,电场值的快变抖动会变慢,且逐渐转为缓慢变化,这种现象主要是云层或者云地之间频繁的放电现象,导致云层电荷量的变化。另外,在雷暴云接近大气电场仪且伴随放电现象时,电场值的变化应由缓慢到频繁的快变抖动,而图7中ADTD、LLS初始阶段都出现了闪电成竖排列,且电场值变化出现不明显的现象,这主要是由于大气电场测站为了减小数据量,在没有雷暴云靠近时采样率为1 min 1个数据,当雷暴云靠近且伴随着电场值达到一定强度时触发采样率为1 min 1 s一个数据,这导致了初始阶段采样间隔过大,部分电场数据丢失,而初始阶段是目前利用大气电场仪开展雷电监测预警的关键阶段,因此这是大气电场仪目前需要改进的地方。
图7 天山区站大气电场强度变化及周围40 km闪电分布
单独利用大气电场仪开展雷电监测临近预警(旅游景区、机场、油库等)存在的主要问题:易受外界电磁场干扰,触发阈值难以判断,导致出现虚警太多。因此,目前国内外大多采用闪电定位和大气电场两种设备相结合的方法。大气电场仪可以对其上空一定半径内的云层带电状况进行监测,闪电定位仪能够监测闪电的发生时刻及位置,因而可以从宏观上看到闪电的分布及走向,但是新疆并不完全适用这种方法,因西北地区不是闪电的高发区,部分局部小范围对流天气过程中地闪频次只有几次,甚至一次,由于监测地闪频次太少,无法根据闪击点与大气电场之间存在的对应关系,找到较为理想的预报因子。因此,在新疆结合大气电场仪开展雷电临近预警需要借助其他手段或方法,如将星地多源数据融合再与大气电场仪结合。
(1)新疆星地多源闪电观测资料中LLS对地闪探测效率最高,LMIG和WWLLN对地闪探测效率较低。卫星全天空、实时监测,可以弥补南疆地面观测站覆盖不到的区域,两者可互为补充。
(2)强对流天气移动路径上ADTD、LLS、WWLLN三者定位结果在时间和空间上对应基本一致,仅是探测效率上有所差别。
(3)新疆部分LMIG数据存在时间和空间上的偏差,使用前需要利用相关资料进行数据检验。
(4)闪电发生位置距离大气电场仪较近时,电场值会出现频繁的快变抖动现象,而随着闪电发生位置越来越远,电场值的快变抖动频率会变慢,且逐渐转为缓慢变化。新疆大气电场仪起始阶段采样间隔过大导致部分数据丢失,不利于预警阈值的设定,需进一步改进。大气电场仪可以作为雷电临近预警的一种辅助手段,下一步将对存在偏差的LMIG数据进行分析,开展多源数据的融合,更加准确地分析新疆地区闪电活动特征,并结合大气电场仪建立雷电监测预警方法。