一种欧亚大陆冷锋自动识别订正方法

本发明涉及天气系统自动识别技术领域，具体涉及一种欧亚大陆冷锋自动识别订正方法。

背景技术：

锋面是冷暖气团的交界面，其附近往往伴随着重要的天气系统和天气现象。我国冷锋活动频繁，常引起灾害性天气，特别是近年来欧亚大陆极端冷事件频发，锋面活动得到了广泛的关注。锋面分析在天气预报的业务中是一项重要的工作。业务中的其他分析项目目前基本已经实现自动化分析，但锋面分析仍以人工分析为主。人工分析要花费预报员宝贵的业务时间，且具有主观性。因此，实现锋面自动识别是气象业务现代化的迫切需求，对气象科研工作也具有十分重要的意义。

由于锋线是一条不规则且没有固定二维结构的线条，所以其识别难度要高于有闭合等值线的天气系统。虽有不少工作提出了一些锋面自动识别算法，但这些算法的主要关注点在大洋风暴轴区域，并不能直接应用于欧亚大陆上的锋面识别。因此，需要提出一种适用于欧亚大陆冷锋的自动识别订正方法，帮助实现业务预报中的锋面分析自动化。

技术实现要素：

为了解决上述问题，本发明提供一种欧亚大陆冷锋自动识别订正方法，依次通过确定高空冷锋锋区、初始锋面复选点以及订正锋面复选点的设计思路，能够实现锋面的自动识别，在一定程度上消除了人工分析锋面的主观性，并对天气预报业务工作中的锋面分析自动化做出积极贡献。

为了实现以上目的，本发明采取的一种技术方案是：

一种欧亚大陆冷锋自动识别订正方法，包括如下步骤：s10下载的850hpa风场、温度场资料，用grads将资料转换成二进制格式，并按年、月储存获得气象参数；s20使用所述气象参数定位高空冷锋锋区，根据所述冷锋锋区确定锋区的暖边界，获得初始锋面复选点；s30将所述初始锋面复选点拟合平滑，绘制锋面，得到自动识别的初始冷锋数据集，并将自动识别的所述初始冷锋数据集的初始冷锋经纬度输出为文本格式；s40使用matlab中max()函数筛选出所述初始冷锋经纬度所在格点中位于最东、最西、最南、最北的格点，分别向东、西、南、北方向扩展5°，获得冷锋订正范围；s50在所述冷锋订正范围，以西北风向的格点为起始点，用起始点的风向减去其同纬度东侧格点的风向，向东计算同一纬度上西北风的风向转变度数；s60将所述风向转变度数大于0的格点定义为具有西北风逆转特征的格点，筛选其中最东侧、最南侧的点，记为订正锋面复选点；s70将所述订正锋面复选点进行拟合平滑，得到订正后的地面冷锋数据集。

进一步地，所述步骤s20包括：s21使用分辨率为0.25°×0.25°的era-5再分析资料，在经纬度范围为10°～70°n、0～160°e的区域，计算热锋面参数τ为在850hpa的温度值，选择|tfp|≤2×10^-11k/m²区域作为锋区；s22计算温度平流其中u、v为850hpa的纬向水平风、经向水平风，为u、v的复合风，x为纬向，y为经向，t为850hpa的温度进行五点平滑100次后的温度值，在所述锋区内选择的区域为冷锋锋区；s23所述冷锋锋区的格点为第一初选点，将所述第一初选点在一个大格点中的总个数占大格点中所有小格点的总数的比例设为第一密度，将所述第一密度大于0.05的大格点设为第二初选点，获得完整锋区，将所述完整锋区中所述第二初选点位于东侧和南侧的格点保留，确定锋区的暖边界，将暖边界第二初选点中位于东侧和南侧的第一初选点保留，获得初始锋面复选点。

进一步地，所述步骤s23以满足|tfp|≤2×10^-11k/m²、的分辨率为0.25°×0.25°格点作为冷锋锋区的第一初选点，2.5°×2.5°所在的范围为一个大格点，计算一个2.5°×2.5°的大格范围内所述第一初选点的个数获得所述第一密度。

进一步地，所述步骤s40输入地面风场资料，通过matlab中atan()函数计算所述冷锋订正范围内各格点上的风向，筛选出西北风向的格点，以风向为西北风的格点为起始点，用起始点的风向减去其同纬度东侧格点的风向，向东计算同一纬度上西北风的风向转变度数，计算公式如下：

δai＝ai-ai+1

其中，δai为西北风的风向转变度数，ai为向为西北风的格点的风向，ai+1为西北风的格点同纬度东侧格点的风向。

本发明的上述技术方案相比现有技术具有以下优点：

本发明的一种欧亚大陆冷锋自动识别订正方法，依次通过确定高空冷锋锋区、初始锋面复选点以及订正锋面复选点的设计思路，能够实现锋面的自动识别，在一定程度上消除了人工分析锋面的主观性，并对天气预报业务工作中的锋面分析自动化做出积极贡献。

附图说明

下面结合附图，通过对本发明的具体实施方式详细描述，将使本发明的技术方案及其有益效果显而易见。

图1所示为本发明一实施例的一种欧亚大陆冷锋自动识别订正方法的流程示意图；

图2所示为本发明一实施例的步骤s20的流程示意图；

图3所示为本发明一实施例的第一初选点在分辨率为0.25°×0.25°小格点中的位置图；

图4所示为本发明一实施例的第二初选点在分辨率为2.5°×2.5°大格点中的位置图；

图5所示为本发明一实施例的暖边界位置图；

图6所示为本发明一实施例的初始锋面复选点的位置图；

图7所示为本发明一实施例的步骤s40确定订正范围的流程示意图(以东边界为例)；

图8所示为本发明一实施例的获得的初始冷锋与地面风场关系的实例图；

图9所示为本发明一实施例的西北风逆转度数与初始冷锋、订正后冷锋实例图；

图10所示为本发明一实施例的自动识别初始锋、订正过后的地面锋、人工分析锋的对比图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本实施例提供了一种欧亚大陆冷锋自动识别订正方法，如图1所示，包括如下步骤：s10下载的850hpa风场、温度场资料，用grads将资料转换成二进制格式，并按年、月储存获得气象参数。s20使用所述气象参数定位高空冷锋锋区，根据所述冷锋锋区确定锋区的暖边界，获得初始锋面复选点。s30将所述初始锋面复选点拟合平滑，绘制锋面，得到自动识别的初始冷锋数据集，并将自动识别的所述初始冷锋数据集的初始冷锋经纬度输出为文本格式。s40使用matlab中max()函数筛选出所述初始冷锋经纬度所在格点中位于最东、最西、最南、最北的格点，分别向东、西、南、北方向扩展5°，获得冷锋订正范围。s50在所述冷锋订正范围，以西北风向的格点为起始点，用起始点的风向减去其同纬度东侧格点的风向，向东计算同一纬度上西北风的风向转变度数。s60将所述风向转变度数大于0的格点定义为具有西北风逆转特征的格点，筛选其中最东侧、最南侧的点，记为订正锋面复选点。s70将所述订正锋面复选点进行拟合平滑，得到订正后的地面冷锋数据集。

所述步骤s10中，所述850hpa温度场、850hpa风场、地面风场资料来自于欧洲中期天气预报中心(https://cds.climate.copernicus.eu/#！/search？text＝era5&type＝dataset)的era-5再分析资料。资料时间分辨率为12小时，取一日两次，分别为世界时00时、12时；空间分辨率为0.25°×0.25°。

如图2所示，所述步骤s20包括：s21使用分辨率为0.25°×0.25°的era-5再分析资料，在经纬度范围为10°～70°n、0～160°e的区域，计算热锋面参数τ为在850hpa的温度值，选择|tfp|≤2×10^-11k/m²区域作为锋区。s22计算温度平流其中u、v为850hpa的纬向水平风、经向水平风，为u、v的复合风，x为纬向，y为经向，t为850hpa的温度进行五点平滑100次后的温度值，在所述锋区内选择的区域为冷锋锋区。如图3～6所示，s23所述冷锋锋区的格点为第一初选点，将所述第一初选点在一个大格点中的总个数占大格点中所有小格点的总数的比例设为第一密度，将所述第一密度大于0.05的大格点设为第二初选点，获得完整锋区，将所述完整锋区中所述第二初选点位于东侧和南侧的格点保留，确定锋区的暖边界，将暖边界第二初选点中位于东侧和南侧的第一初选点保留，获得初始锋面复选点。所述第一初选点在分辨率为0.25°×0.25°格点中构成的锋区构成的锋区常常不成片、有中断，不能绘制冷锋；故将研究范围10°～70°n、0～160°e的区域划分为0.25°×0.25°的大格点，获得完整、成片的冷锋锋区。所述步骤s23以满足|tfp|≤2×10^-11k/m²、的分辨率为0.25°×0.25°格点作为冷锋锋区的第一初选点，2.5°×2.5°所在的范围为一个大格点，计算一个2.5°×2.5°的大格范围内所述第一初选点的个数获得所述第一密度。

所述步骤s30中需要将初始锋面复选点进行二阶多项式拟合平滑，绘制锋面，删除长度低于800km的冷锋，得到自动识别的初始冷锋数据集。此时得到的初始冷锋数据集是根据850hpa的气象要素得到的，但实际业务预报中的冷锋位于地面，因此需要结合地面气象要素对初始冷锋进行订正。所述s40，如图7，以东边界为例给出了确定冷锋订正范围的流程图。

所述步骤s50输入地面风场资料，通过matlab中atan()函数计算所述冷锋订正范围内各格点上的风向。如图8所示，填色部分为地面风向，锋线由初始锋面复选点构成，图中可以明显看到在此锋线东边存在风向逆转，即西北风向东，故此锋线应向东订正。s60筛选出西北风向的格点，以风向为西北风的格点为起始点，用起始点的风向减去其同纬度东侧格点的风向，向东计算同一纬度上西北风的风向转变度数，计算公式如下：

δai＝ai-ai+1

其中，δai为西北风的风向转变度数，ai为向为西北风的格点的风向，ai+1为西北风的格点同纬度东侧格点的风向。如果δai>0,认为该格点具有西北风逆转特征。筛选具有西北风逆转特征的格点中最东侧、最南侧的点，将其作为订正复选点；当同一纬度上出现两个满足条件的格点时，取东侧的点。当初始格点东侧、南侧没有风向逆转时保留原有格点位置，如图9所示。图9中填色部分为风向转变度数大于0的格点。两条锋线分别为初始冷锋和订正后地面冷锋。

s70将所述订正锋面复选点进行拟合平滑，如图10中黑色实线，得到订正后的地面冷锋，图10中深灰色实线为人工分析得到的冷锋，浅灰色实线为初始冷锋。

以上所述仅为本发明的示例性实施例，并非因此限制本发明专利保护范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。