一种基于天气雷达的短临智能外推方法技术

本发明专利技术提供了一种基于天气雷达的短临智能外推方法，所述方法包括对风暴单体历史个例的收集、识别和追踪，将获取的n+1项风暴单体轨迹序列作为范例存储于数据库中，当进行天气预测时，将获取的n项实时风暴单体序列与数据库中的所有范例进行相似度计算，找出相似度最高的范例，将该范例的第n+1项作为下个时次的风暴单体，以此计算出所有单体轨迹下个时次的状态，构成所预测的下个时次的天气状况，同时将新获取的实时风暴单体序列作为范例，存储于数据库中，以进行自学习。所述方法综合机器学习，构建起基于范例推理的短临外推专家系统，并结合客观预报的方法对天气变化进行相似性预测，提高了灾害性短时天气的预警水平。

【技术实现步骤摘要】
一种基于天气雷达的短临智能外推方法
本专利技术涉及气象学天气预测
，特别涉及一种基于天气雷达的短临智能外推方法。
技术介绍
对流天气的临近预报是指对未来0～2小时的对流天气系统及其所伴随的灾害性天气的发生、发展、演变和消亡的预报，强对流天气是影响我国的主要灾害性天气之一，具有空间尺度小、生命史短、突发性强、发展演变迅速、破坏力大等特点，因此，针对强对流天气的临近预报一直是气象预报工作的重点和难点。目前传统的天气预测方法有两种，即数值天气预报(NWP)法和基于光流矢量的天气预测方法。数值天气预报法是一种先对大气进行物理建模，再通过模拟物理模型进行预测的方法，该方法的优势体现在对于长时间范围内的天气预测比较准确，但是对0～2小时内的天气状况无法预测，因此不宜使用NWP方法进行短临预报。另一种在实际中得以应用的天气预测方法是基于光流矢量的方法，该方法是先通过两个雷达回波图像来估算光流矢量，该光流矢量可以理解为云的运动方向，然后根据最后一张雷达图进行合理的外推，得到最终的预测结果，该方法在0～2小时内的短临预测结果准确度更高。以上基于雷达资料的短临外推预报方法的研究思路主要是围绕着单个天气变化过程中的对流单体演变过程，并未考虑不同天气变化过程中的不同对流单体间的类比，且缺乏自学习能力。
技术实现思路
专利技术目的：针对上述现有技术存在的问题和不足，本专利技术提供了一种基于天气雷达的短临智能外推方法，综合机器学习方法，构建起基于范例推理的短临外推专家系统，并结合客观预报的方法对天气变化进行相似性预测，提高了灾害性短时天气的预警水平。技术方案：本专利技术提供了一种基于天气雷达的短临智能外推方法，包括以下步骤：步骤a：收集某个地区天气个例雷达基数据；步骤b：提取每一时刻雷达基数据中的所有风暴单体，并提取、追踪其特征，获得若干条长度不等风暴单体轨迹序列；步骤c：将获得的若干条长度不等风暴单体轨迹序列，切割成长度为n+1的风暴单体轨迹范例存储于数据库中；步骤d：当进行天气预测时，读取当前n个时次的雷达基数据，利用SCIT识别、追踪，获得若干条长度为n的风暴单体轨迹序列，将该序列分别与数据库中所有历史风暴单体的前n项的轨迹范例进行相似度计算，找出相似度最高的范例，将该范例的第n+1项作为下个时次的风暴单体。按此方法计算出所有单体轨迹下个时次的状态，构成所预测的下个时次的天气状况，写为雷达基数据；步骤e：按照步骤d进行迭代预测，直至推算出未来2小时全部的天气状况；步骤f：将新获取的风暴单体序列切割成长度为n+1的范例，存储于数据库中，进行自学习，丰富范例库。所述步骤b中提取风暴单体的方法为风暴单体识别追踪算法，提取的风暴单体特征包括当前季节、风暴体序列号、单体质心坐标、所在的仰角、质量权重体积、最大反射率、移动方向、移动速度。所述步骤d中，在进行相似度计算时，所筛选的历史风暴单体的形成位置位于地理相似判断范围内，所述地理相似判断范围为以待预测风暴单体的形成位置为圆心，半径不大于预设纬距的范围。针对风暴单体的移动方向，设定与该移动方向偏差预设角度的范围，作为移向相似判断范围，针对风暴单体的移动速度，设定预设大小的上下波动范围，作为移速相似判断范围。有益效果：与现有技术相比，本专利技术的优点在于：1、本专利技术提供了一种基于天气雷达的短临智能外推方法，不仅考虑天气系统本身的演变信息，还将同一地区历史相似天气系统的演变信息考虑在内，提高了强对流短临天气预报的预警能力和准确性。2、本专利技术所采用的预测方法只需要一台配置较高的小型服务器即可，相对于传统数值预报需要大型计算机进行辅助运行而言，大大缩短了预测所需要的时间，其次在预测过程中不需要进行大量的参数输入，大大提高了气象预测的工作效率。附图说明图1为基于范例推理的雷达基数据外推预报方法流程图；图2为风暴单体识别流程图；图3为质心法追踪流程图；图4为二维风暴分量的识别示意图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。如图1所示，基于机器学习的短临外推方法的预报准确性依赖从范例库中提取的历史对流单体演变轨迹范例，范例质量好则预报准确性高。因此能否从提取库中选楚最优的范例取决于不同对流单体演变序列的相似性度量算法设计，考虑到对流单体具有很强的季节性特征，拟采用基于非平稳季节ARIMA模型来进行拟合，更科学地实现不同序列片段间的相似性度量。各模块的方法如下：(1)天气个例雷达基数据的收集、质量控制；(2)基于雷达基数据的对流单体识别与跟踪；(3)基于对流单体演变的时间序列进行特征选择；(4)构建基于时序片段的对流单体演变范例库；(5)设计适合对流单体演变序列的相似性度量算法；(6)基于类比方法进行对流单体演变的短临外推；(7)业务检验与对流单体演变范例库的自学习机制。本专利技术设计的一种天气变化预测方法在实际应用过程当中，基于历史气象雷达基数据，针对风暴单体的移动路径进行预测，获得未来0～2小时预测的天气变化过程，具体包括如下步骤：步骤a：收集某个地区天气个例雷达基数据。步骤b：如图2所示，从一块地区的历史气象雷达基数据中识别风暴单体，具体步骤如下：步骤b1：首先在体积扫描的每个径向上搜索强度大于一定阀值(默认为30dBz)的连续点，合并成有一定权重长度的一维风暴段；步骤b2：如图4所示，在PPI(像素密度)层中按方位距离、长度重叠的相关性将风暴段合并成有一定权重面积(默认为10平方千米)的二维风暴分量；步骤b3：最后按空间垂直相关将风暴分量合并成具有一定权重体积的三维风暴体，识别完成后，获得风暴单体；步骤b4：提取所识别风暴单体的特征，包括季节、质心位置、质量权重体积、最大反射率等，具体详见下表1所示：表1风暴单体特征值步骤c：如图3所示，采取质心追踪法追踪风暴单体的轨迹，所谓质心跟踪法，主要是根据前一时刻风暴的运动矢量和当前风暴与前一时刻风暴可能存在的位置关系对风暴单体的轨迹进行追踪，同时根据位置关系结果估计两个风暴单体的合并、分裂、新生和消亡。该算法的基本思想是：存在相似特征(包括体积、强度等)的风暴单体优先配对，近距离的风暴单体优先配对，但两者之间的最大速度不得超过限定值。风暴跟踪根据连续时间内的多个体扫识别出的风暴单体及其特征，通过对最后两个体扫匹配风暴体实现风暴跟踪。风暴单体的轨迹追踪的具体步骤如下：步骤c1：用前一时刻t0每一个风暴体的运动矢量估计其在当前t1时刻可能的质量权重中心位置，如果该风暴体在t0时刻扫描体中是第一次测得，则用缺省的运动矢量，即：其中Xexp、Yexp为缺省的运动矢量；vpast为上个时次的速度标量；θ为速度的方位角；Dt＝t1-t0，表示两个风暴单体扫描间的时间间隔。关键问题是如何给出t0时刻新生风暴体的初始速度，特别是当t0时刻是程序处理的第一个体扫描资料，无历史路径时，就必须较为准确的给出每一个风暴体的初始速度。通过分别计算两时刻风暴体的群体质量权重中心，来获取新生风暴体初始速度缺省值，方法如下：t0时刻所有风暴体的群体质量权重中心为：t1时刻所有风暴体的群体质量权重中心为：则t0时刻任一风暴的速度缺省值为：其中，Vz0i，Vz1i分别为t0和t1时刻任一风暴的反射率因子权重体积，Xz0i，Yz0i，Xz1i，Yz1i，分别为t0和t1时刻任本文档来自技高网...

【技术保护点】
1.一种基于天气雷达的短临智能外推方法，包括以下步骤：步骤a：收集天气个例雷达基数据；步骤b：提取每一时刻雷达基数据中的所有风暴单体，并提取、追踪其特征，获得若干条长度不等风暴单体轨迹序列；步骤c：将获得的若干条长度不等风暴单体轨迹序列，切割成长度为n+1的风暴单体轨迹范例存储于数据库中；步骤d：在进行天气预测时，读取当前n个时次的雷达基数据，利用SCIT识别、追踪，获得若干条长度为n的风暴单体轨迹序列，将该序列分别与数据库中所有历史风暴单体的前n项的轨迹范例进行相似度计算，找出相似度最高的范例，将该范例的第n+1项作为下个时次的风暴单体，按此方法计算出所有单体轨迹下个时次的状态，构成所预测的下个时次的天气状况，写为雷达基数据；步骤e：按照步骤d进行迭代预测，直至推算出未来2小时全部的天气状况；步骤f：将新获取的风暴单体序列切割成长度为n+1的范例，存储于数据库中，进行自学习，丰富范例库。

【技术特征摘要】
1.一种基于天气雷达的短临智能外推方法，包括以下步骤：步骤a：收集天气个例雷达基数据；步骤b：提取每一时刻雷达基数据中的所有风暴单体，并提取、追踪其特征，获得若干条长度不等风暴单体轨迹序列；步骤c：将获得的若干条长度不等风暴单体轨迹序列，切割成长度为n+1的风暴单体轨迹范例存储于数据库中；步骤d：在进行天气预测时，读取当前n个时次的雷达基数据，利用SCIT识别、追踪，获得若干条长度为n的风暴单体轨迹序列，将该序列分别与数据库中所有历史风暴单体的前n项的轨迹范例进行相似度计算，找出相似度最高的范例，将该范例的第n+1项作为下个时次的风暴单体，按此方法计算出所有单体轨迹下个时次的状态，构成所预测的下个时次的天气状况，写为雷达基数据；步骤e：按照步骤d进行迭代预测，直至推算出未来2小时全部的天气状况；步骤f：将新获取...

【专利技术属性】
技术研发人员：耿焕同，戴中斌，周山胜，韩伟民，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32