技术领域
[0001] 本
发明涉及数值天气预报技术领域，更具体地，涉及一种基于动力松弛逼近的闪电定位资料同化方法。
背景技术
[0002] 闪电定位仪是一种能够提供
对流云中闪电发生
位置、时间、强度信息的气象探测设备，具有探测
精度高、
覆盖范围广以及不受地形影响的优点。闪电活动与对流性天气系统的动力过程和微物理过程密切相关，因此闪电观测信息可用于数值预报模式资料同化，用于改善初始场的动力条件和
水汽条件分布。
[0003] 目前，闪电资料同化方法分为两类，第一类方法是同化，通过变分同化或松弛逼近，由闪电观测信息反演的降水率、水汽或
冰相粒子，改善初始场水物质的分布；第二类方法是电观测间接影响垂直运动，该方法将闪电观测和数值预报模式对流参数化方案相结合，利用闪电观测控制对流参数化方案的触发。然而，由于高
分辨率数值模式不再需要对流参数化方案，因此第二类方案已不再可用；而第一类方法主要针对水物质场的初始化，缺乏与动力场的联系，容易出现水物质与动力场不匹配的问题，进而导致水物质过快降落。此外，考虑到雷达回波的同化应用也针对水物质场的初始化，使得闪电定位和雷达回波两类资料应用的同质化，导致不能满足充分挖掘多源资料的综合初始化应用效果。
发明内容
[0004] 本发明为克服上述
现有技术所述的闪电资料同化中容易出现水物质与动力场不匹配的
缺陷，提供一种基于动力松弛逼近的闪电定位资料同化方法。
[0005] 为解决上述技术问题，本发明的技术方案如下：
[0006] 一种基于动力松弛逼近的闪电定位资料同化方法，包括以下步骤：
[0007] S1：采集闪电观测资料并对其进行预处理，得到其闪电
频率；
[0008] S2：通过反演公式对所述闪电频率进行反演，得到其最大垂直速度；
[0009] S3：根据组合反射率值分类，对归一化系数进行统计；
[0010] S4：根据所述最大垂直速度和所述统计的归一化系数，计算获得一维垂直速度估算廓线；
[0011] S5：在数值预报模式积分过程中，采用松弛逼近法将所述垂直速度估算廓线引入数值模式初始场，并结合雷达回波反演水物质场初始化，得到与动力场匹配的初始场。
[0012] 本技术方案中，首先基于闪电频率和垂直速度的动力统计关系，利用闪电定位观测反演获得最大垂直速度，然后利用模式输出数据统计的归一化系数将其扩展为垂直速度廓线，最后在数值预报模式积分过程中，通过松弛逼近方法引入所述垂直速度估算廓线信息，用于改善初始场的垂直动力场分布，并结合雷达回波反演水物质场初始化，得到与动力场匹配的初始场，可更好地挖掘多源观测资料的总和初始化应用效果。本技术方案基于闪电和对流云垂直运动的物理关系，通过动力松弛逼近实现了对闪电定位观测的同化，同时可结合雷达回波反演云水物初始化应用，能够发挥更好的多源资料综合应用效果。
[0013] 优选地，S1步骤中，对所述闪电观测资料进行预处理的步骤包括：按照一定空间分辨率和
时间窗，统计落在网格内的闪电频次，得到所述闪电频率。
[0014] 优选地，S2步骤中，所述最大垂直速度Wmax的计算公式如下：
[0015] Wmax＝a(Fr)b
[0016] 其中，Fr表示闪电频率，a和b为随着下垫面类型变化的经验系数。
[0017] 优选地，当所述S2步骤中所采集的闪电观测资料场景为陆地时，所述经验系数a和b取值为14.66和0.22；当所述S1步骤所采集的闪电观测资料场景为海洋时，所述经验系数a和b取值为8.69和0.22。
[0018] 优选地，S3步骤中，统计归一化系数的计算公式如下：
[0019]
[0020] 其中，σL(zk)表示L类垂直第zk层的归一化系数，Wm(i,j)为格点坐标为(i,j)的气柱内的最大垂直速度；W(i,j,zk)为垂直第zk层的垂直速度；M表示统计样本总数，即M＝I×J×N；I和J表示水平格点数，N为预报输出次数。
[0021] 优选地，S4步骤中，一维垂直速度估算廓线West的计算公式如下：
[0022] West(zk)＝σL(zk)·Wmax。
[0023] 优选地，S5步骤中，所述松弛逼近公式如下：
[0024]
[0025] 其中， 分别表示积分时步为t+1和t时的垂直速度值；Wb表示初始时刻垂直速度值，且当t＝0时 nt为积分步数。
[0026] 与现有技术相比，本发明技术方案的有益效果是：基于闪电和对流云垂直运动的物理关系，通过动力松弛逼近引入由闪电观测反演的垂直速度信息，解决了水物质与动力场不匹配的问题，改善了初始场垂直运动的合理性缺陷，提高了降水的模拟或预报准确率；具有计算量小、应用范围广的特点，可嵌入应用于各种数值天气预报模式。
附图说明
[0027] 图1为本
实施例的一种基于动力松弛逼近的闪电定位资料同化方法的
流程图。
[0028] 图2为闪电定位资料同化前后预报的1小时降水与观测比较图。
[0029] 图3为闪电定位资料同化前后预报的
站点降水时间序列与观测对比图。
具体实施方式
[0030] 附图仅用于示例性说明，不能理解为对本
专利的限制；
[0031] 为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；
[0032] 对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
[0033] 下面结合附图和实施例对本发明的技术方案做进一步的说明。
[0034] 如图1所示，为本实施例的一种基于动力松弛逼近的闪电定位资料同化方法的流程图。
[0035] 本实施例提出一种基于动力松弛逼近的闪电定位资料同化方法，包括以下步骤：
[0036] S1：采集闪电观测资料并对其进行预处理，得到其闪电频率。
[0037] 本步骤中，对所采集的闪电观测资料间隔12分钟统计落在0.03°×0.03°网格内的闪电频次，得到其闪电频率。
[0038] S2：通过反演公式对所述闪电频率进行反演，得到其最大垂直速度。
[0039] 本步骤中，最大垂直速度Wmax的计算公式如下：
[0040] Wmax＝a(Fr)b
[0041] 其中，Fr表示闪电频率，a和b为随着下垫面类型变化的经验系数，且当所采集的闪电观测资料场景为陆地时，经验系数a和b取值为14.66和0.22；所采集的闪电观测资料场景为海洋时，经验系数a和b取值为8.69和0.22。
[0042] S3：根据组合反射率值分类，对归一化系数进行统计。
[0043] 本步骤中，先根据组合反射率值CR对各层垂直速度进行分类，其中，将18≤CR<30划分为第一类L1；将30≤CR<40为第二类L2；将CR≥40为第三类L3。
[0044] 本步骤中对归一化系数进行统计的计算公式如下：
[0045]
[0046] 其中，σL(zk)表示L类垂直第zk层的归一化系数，Wm(i,j)为格点坐标为(i,j)的气柱内的最大垂直速度；W(i,j,zk)为垂直第zk层的垂直速度；M表示统计样本总数，即M＝I×J×N；I和J表示水平格点数，N为预报输出次数。
[0047] S4：根据所述最大垂直速度和所述统计的归一化系数，计算获得一维垂直速度估算廓线。
[0048] 本步骤中，一维垂直速度估算廓线West的计算公式如下：
[0049] West(zk)＝σL(zk)·Wmax。
[0050] S5：在数值预报模式积分过程中，采用松弛逼近法将所述垂直速度估算廓线引入数值模式初始场，并结合雷达回波反演水物质场初始化，得到与动力场匹配的初始场。
[0051] 本步骤中，所采用的松弛逼近公式如下：
[0052]
[0053] 其中， 分别表示积分时步为t+1和t时的垂直速度值；Wb表示初始时刻垂直速度值，且当t＝0时 nt为积分步数，本实施例中的积分步数取值为12。
[0054] 在具体实施过程中，利用数值天气预报模式GRAPES meso，选用2018年5月7日广东省的一次飑线为例进行发明方法的实施演示。
[0055] 本实施例所述的闪电定位资料来源于粤港澳闪电定位系统，该定位系统拥有17个探测子站，可实时提供云对地闪电回击的经纬度、时间、雷
电流幅度等，其中发生闪电的经纬度、时间、次数被用于本次实施例。对上述闪电的经纬度、时间、次数等闪电观测资料进行处理，统计间隔12分钟统计落在0.03°×0.03°网格内的闪电频次，得到其闪电频率Fr。然后对上述闪电频率Fr通过反演公式对所述闪电频率进行反演，得到其最大垂直速度Wmax，其中，本实施例中经验系数a和b的取值为14.66和0.22。
[0056] 本实施例中，利用ECMWF全球高分辨率预报场生成区域数值天气模式GRAPES meso的初始条件和边界条件。其中，GRAPES meso模式网格选取两层网格，单向嵌套，第一层网格水平分辨率为0.03°，格点数为913×513，第二层网格水平分辨率为0.01°，格点数为1201×901。本实施例选择2017年5月7日00时-5月8日00时(UTC)的时间范围内进行一次飑线的预报，其中选择2017年5月7日08时0.01°分辨率的预报场作为本实施案例的初始场。在数值预报模式积分过程中，在上述初始场积分的前12时步中，采用松弛逼近法将所述垂直速度估算廓线引入数值模式初始场，并对飑线进行数值预报。
[0057] 如图2、3所示，分别为本实施例的闪电定位资料同化前后预报的1小时降水与观测比较图，以及闪电定位资料同化前后预报的站点降水时间序列与观测对比图。由图可知，应用本实施例所提出的一种基于动力松弛逼近的闪电定位资料同化方法，提高了0-3小时的强降水预报精度，其中邻域空间检测评分提高了约0.13，即相对于不同化试验提高约46％。本实施例对站点降水预报也有明显的改善效果，其中10时的小时降水量绝对误差由同化前的23mm减小为7mm。
[0058] 本实施例基于闪电和对流云垂直运动的关系，通过动力松弛逼近法实现了对闪电定位观测的同化，能够有效增强垂直运动，解决了初始场中水物质与动力场不匹配的问题，改善了初始场垂直运动的合理性缺陷，提高了降水的模拟或预报的准确率。此外，本实施例可结合雷达回波反演云水物初始化应用，能够发挥更好的多源资料综合应用效果。
[0059] 相同或相似的标号对应相同或相似的部件；
[0060] 附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制；
[0061] 显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的
基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何
修改、等同替换和改进等，均应包含在本发明
权利要求的保护范围之内。
