一种提高GPT2w模型计算气象参数精度的方法与流程

一种提高gpt2w模型计算气象参数精度的方法
技术领域
1.本发明涉及气象应用研究领域，尤其涉及一种提高gpt2w模型计算气象参数精度的方法。


背景技术：

2.gpt是英文global pressure and temperature(全球气温气压)的简称，gpt系列模型可用于计算全球任意位置的气温、气压和水汽压等各种气象参数。应用最广且精度较高的全球气温气压模型主要为gpt2w模型。gpt2w模型通过输入站点处的经度、纬度、高程及年积日来计算距离站点最近的四个格网点上的气象参数并改正到站点高程处，最后利用双线性插值得到站点处的气象参数，其模型表达式如式(1)所示：
[0003][0004]
式(1)中，doy为年积日，x为格网点的气温、气压等气象参数，a0为对应参数的平均值，a1和b1为对应参数的年变化振幅，a2和b2为对应参数的半年变化振幅，以上参数可以从模型提供的格网参数文件中获得。
[0005]
gpt2w模型在计算每个年积日对应的气象参数时要对格网参数文件进行分类初始化存储，这导致每次进行气象参数的计算时将会消耗大量时间(计算7000个年积日对应的气象参数大约消耗3至4小时)，并且该模型仅考虑了气象参数的年周期振幅和半年周期振幅，这导致在不同的月份模型计算出的气象参数精度和稳定性有所差异并在个别月份计算出的气象参数与实际参数值具有较大差异。


技术实现要素：

[0006]
发明目的：针对上述现有技术的缺陷，本发明提供一种提高gpt2w模型计算得到的气温、气压参数精度和计算效率的方法。
[0007]
技术方案：本发明提高gpt2w模型计算气象参数精度的方法包括如下步骤：包括如下步骤：
[0008]
(1)对气象站点实测数据进行预处理，将实测数据中的粗差剔除，保存无误的实测数据；
[0009]
(2)将预处理后的实测数据加载到gpt2w模型中计算气象参数并保存，
[0010]
(3)根据步骤(2)中的气象参数计算各气象站点各月份的年均气温、气压误差值；
[0011]
(4)将步骤(3)中获取的各气象站点各月份的年均气温、气压误差值加入各气象站点对应的经度、纬度并保存，生成改正值文件；
[0012]
(5)对(4)中生成的改正值文件进行距离加权改正，得到新的格网参数文件并将其进行分类存储；
[0013]
(6)向gpt2w模型导入测站点的经度、纬度、高程及观测时间，按观测时间对应月份读取步骤(5)中生成的格网参数文件计算气温、气压参数。
[0014]
步骤(1)具体步骤如下：
[0015]
(11)筛选出实测数据中不含气温、气压值的观测时间点并将其剔除。
[0016]
(12)将经过步骤(1)1筛选得到的实测数据分别按照气温、气压值依次从大到小和从小到大排列，筛选出不符合实际的错误气温、气压值并将其剔除。
[0017]
步骤(2)具体步骤如下：
[0018]
(2)1、在gpt2w模型中导入其中一年某个省份中实测数据的经度、纬度、高程及观测时间并计算气温、气压参数；
[0019]
(2)2、将步骤(2)1中处理好的数据按原省份测站号进行保存，此时数据中包括实测气温、气压数据和模型计算出的气温、气压数据；
[0020]
(2)3、依次将各年份的数据重复(2)1和(2)2的操作，直至得到全部年份的数据。
[0021]
步骤(3)具体步骤如下：
[0022]
(31)将经步骤(2)中处理的数据每个观测时间点的实测气温、气压值与模型计算得到的气温、气压值相减得到每个观测时间的气温、气压误差值；
[0023]
(32)筛选出每年中各月份的气温、气压误差值，对各月份的气温、气压误差值取平均值，得到各月份的气温、气压月均误差；
[0024]
(33)将每年中相同月份的气温、气压月均误差值取平均值得到各月份的年均气温、气压误差值。
[0025]
步骤(4)中，改正值文件是通过将步骤(3)中得到的各气象站点各月份的年均气温、气压误差值加入各气象站点对应的经度、纬度并保存生成。
[0026]
步骤(5)中，距离加权改正的具体步骤如下：
[0027]
(51)导入步骤(4)中得到的各格网点各月份的气温、气压改正值文件，并依次读出每个格网点的经纬度；
[0028]
(52)根据每个格网点的经、纬度，利用经、纬度取整来找出距离原格网点最近的四个格网点，其具体步骤为：将每个格网点的经、纬度取整并加上0.5找到最近的格网点，并将原格网点的经、纬度与最近的格网点经、纬度相减得到经、纬度差值，当差值为正时表示为1，差值为负时表示为
‑
1，基于原格网点与最近的格网点的经纬度差值，找到其余三个距离原格网点最近的格网点；
[0029]
(53)找到距离原格网点最近的四个格网点对应的气温、气压平均参数a0，依次将各格网点各月份的气温、气压改正值依据原格网点到最近的四个格网点之间的距离进行加权分配加入到a0，得到全年各月份对应的新格网参数文件，加权分配公式如下所示：
[0030][0031][0032][0033]
[0034][0035][0036][0037][0038]
上式中，d1、d2、d3、d4分别为原格网点到最近的四个格网点的距离，lon0为原格网点经度，lat0为原格网点纬度，lat1、lat2、lat3、lat4分别为离原格网点最近的四个格网点的纬度，lon1、lon2、lon3、lon4分别为离原格网点最近的四个格网点的经度，delta为各格网点各月份的气温、气压改正值，delta1、delta2、delta3、delta4分别为离原格网点最近的四个格网点的距离加权改正值。
[0039]
步骤(5)中，分类存储是指将12个月份对应的新格网参数文件里的各格网点气温、气压参数对应的平均值a0、年变化振幅a1和b1、半年变化振幅a2和b2分类并加上格网点对应的经、纬度按月份进行保存。
[0040]
有益效果：与现有技术相比，本发明具有如下显著优点：显著提高计算出的气象参数精度和稳定性，省去计算所需站点的气象参数时对格网参数文件进行分类初始化存储的时间，提升模型的计算效率，气象参数的计算过程更加方便、快捷、计算结果精确，具有较高的科研和应用价值。
附图说明
[0041]
图1为本发明的方法流程图；
[0042]
图2为2012年改正前后气温误差分布对比图，图中实心点表示气温误差；
[0043]
图3为2013年改正前后气温误差分布对比图，图中实心点表示气温误差；
[0044]
图4为2014年改正前后气温误差分布对比图，图中实心点表示气温误差；
[0045]
图5为2015年改正前后气温误差分布对比图，图中实心点表示气温误差；
[0046]
图6为2016年改正前后气温误差分布对比图，图中实心点表示气温误差；
[0047]
图7为2012年改正前后气压误差分布对比图，图中实心点表示气压误差；
[0048]
图8为2013年改正前后气压误差分布对比图，图中实心点表示气压误差；
[0049]
图9为2014年改正前后气压误差分布对比图，图中实心点表示气压误差；
[0050]
图10为2015年改正前后气压误差分布对比图，图中实心点表示气压误差；
[0051]
图11为2016年改正前后气压误差分布对比图，图中实心点表示气压误差；
[0052]
图12为2016年改正前后气温月均误差对比图(安徽58321)；
[0053]
图13为2016年改正前后气压月均误差对比图(安徽58321)。
具体实施方式
[0054]
下面结合附图和实施例对本发明的技术方案作进一步说明。
[0055]
如图1所示，本发明的提高gpt2w模型计算气象参数精度的方法包括如下步骤：
[0056]
(1)对气象站点实测数据进行预处理，将实测数据中的粗差剔除，保存无误的实测数据；
[0057]
(2)将预处理后的实测数据加载到gpt2w模型中计算气象参数并保存，
[0058]
(3)根据步骤(2)中的气象参数计算各气象站点各月份的年均气温、气压误差值；
[0059]
(4)将步骤(3)中获取的各气象站点各月份的年均气温、气压误差值加入各气象站点对应的经度、纬度并保存，生成改正值文件；
[0060]
(5)对(4)中生成的改正值文件进行距离加权改正，得到新的格网参数文件并将其进行分类存储；
[0061]
(6)向gpt2w模型导入测站点的经度、纬度、高程及观测时间，按观测时间对应月份读取步骤(5)中生成的格网参数文件计算气温、气压参数。
[0062]
以2012
‑
2016年(共5年)中国陆地区域102个国家气象站的气温、气压数据为例：
[0063]
首先，对气象站点实测数据进行预处理，将实测数据中的粗差剔除，保存无误的实测数据；筛选出实测数据中不含气温、气压值的观测时间点并将其剔除，将筛选后剩余的实测数据分别对气温、气压值依次从大到小和从小到大排列，筛选出不符合实际的错误气温、气压值并将其剔除。
[0064]
2016年安徽省58321号气象站点部分实测数据如表1所示：
[0065]
表1 2016年安徽省58321号气象站点部分实测数据
[0066]
省份地市站名站号经度纬度海拔观测时间气温气压水汽压值安徽合肥合肥58321117.3
°
31.785
°
27m2016/12/31 23:002.0℃1024.9hpa6.2hpa安徽合肥合肥58321117.3
°
31.785
°
27m2016/12/31 22:002.8℃1025.1hpa6.4hpa安徽合肥合肥58321117.3
°
31.785
°
27m2016/12/31 21:003.0℃1025.4hpa6.1hpa安徽合肥合肥58321117.3
°
31.785
°
27m2016/12/31 20:004.1℃1025.1hpa6.1hpa安徽合肥合肥58321117.3
°
31.785
°
27m2016/12/31 19:004.8℃1024.9hpa5.8hpa
[0067]
将2012
‑
2016年的气象站点实测数据按年份分别进行保存，依次将每年不同省份的实测数据加载到gpt2w模型并计算气象参数且保存，其具体步骤如下：
[0068]
在gpt2w模型中导入其中一年某个省份中实测数据的经度、纬度、高程及观测时间并计算气温、气压参数；将计算得到的数据按原省份测站号进行保存，此时数据中包括实测气温、气压数据和模型计算出的气温、气压数据；依次将5年的实测数据重复第一步和第二步的操作。
[0069]
处理好的数据如表2所示，表2中为2016年安徽省58321号气象站点部分gpt2w模型计算的数据。
[0070]
表2 2016年安徽省58321号气象站点部分gpt2w模型计算完成数据
[0071][0072][0073]
计算各气象站点各月份的年均气温、气压误差值。
[0074]
对上一步骤中处理好的数据进行处理，将每个观测时间点的实测气温、气压值与模型计算得到的气温、气压值相减得到每个观测时间的气温、气压误差值；筛选出每年中各月份的气温、气压误差值，对各月份的气温、气压误差值取平均值得到各月份的气温、气压月均误差；将每年中相同月份的气温、气压月均误差值取平均值得到各月份的年均气温、气压误差值。
[0075]
获取气温、气压改正值文件：将各气象站点各月份的年均气温、气压误差值加入各气象站点对应的经度、纬度并保存，得到102个格网点各月份的气温、气压改正值文件。
[0076]
对原格网参数文件进行改正：
[0077]
导入上一步骤中得到的102个格网点各月份的气温、气压改正值文件，并依次读出每个格网点的经纬度。
[0078]
根据每个格网点的经、纬度，对经、纬度取整来找出距离格网点最近的四个格网点。
[0079]
以经度差值为
‑
1、纬度差值为举例：此时将最近格网点的经度减1、纬度保持不变即可得到第二个格网点的经、纬度；将最近格网点的经度减1、纬度加1即可得到第三个格网点的经、纬度；将最近格网点的经度保持不变、纬度加1即可得到第四个格网点的经、纬度。其余经、纬度差值以此类推即可得到距离原格网点最近的四个格网点。
[0080]
找到距离原格网点最近的四个格网点中对应的气温、气压平均参数a0，依次将各格网点各月份的气温、气压改正值依据原格网点到最近的四个格网点之间的距离进行加权分配加入到a0，进而得到各月份对应的新格网参数文件(共12个)。对气温、气压改正值进行
距离加权分配到a0中如式(2)至式(9)所示：
[0081][0082][0083][0084][0085][0086][0087][0088][0089]
式(2)至式(9)中，d1、d2、d3、d4分别为原格网点到最近的四个格网点的距离，lon0为原格网点经度，lat0为原格网点纬度，lat1、lat2、lat3、lat4分别为离原格网点最近的四个格网点的纬度，lon1、lon2、lon3、lon4分别为离原格网点最近的四个格网点的经度,delta为各格网点各月份的气温、气压改正值，delta1、delta2、delta3、delta4分别为离原格网点最近的四个格网点的距离加权改正值。
[0090]
将新格网参数文件分类存储：将12个月份对应的新格网参数文件里的各格网点气温、气压参数对应的平均值a0、年变化振幅a1和b1、半年变化振幅a2和b2分类并加上格网点对应的经、纬度按月份进行保存，能够使得gpt2w模型在计算所需站点的气象参数时省去了对格网参数文件进行分类初始化存储的时间。
[0091]
其中部分1月份对应的格网参数文件如表3所示。
[0092]
表3部分1月份对应格网参数文件
[0093]
latlonp:a0a1b1a2b2t:a0a1b1a2b2unduhs39.5
°
106.5
°
86385.53hpa542
‑
142
‑
106
‑
75276.01k
‑
15.1
‑
3.3
‑
1.10.1
‑
40.511397.0639.5
°
107.5
°
85453.53hpa466
‑
161
‑
111
‑
67275.01k
‑
14.9
‑
3.2
‑
1.10.1
‑
37.41486.9539.5
°
108.5
°
86435.00hpa524
‑
158
‑
105
‑
59282.30k
‑
14.9
‑
3.2
‑
1.10
‑
33.711357.9939.5
°
109.5
°
87597.00hpa573
‑
157
‑
101
‑
51282.60k
‑
14.9
‑
3.2
‑
10
‑
29.871247.53
[0094]
gpt2w模型导入测站点的经度、纬度、高程及观测时间，根据观测时间选择对应月份的格网参数文件，即可实现利用gpt2w模型高效计算精度更高的气温、气压参数。
[0095]
表4、6为2012
‑
2016年原gpt2w模型在中国陆地区域102个气象站点得到的气温、气压误差结果统计，表5、7为改进后的gpt2w模型在中国陆地区域102个气象站点得到的气温、气压误差结果统计。由表4至表7可知，改进后的gpt2w模型计算气温、气压的精度有大幅提升，2012
‑
2016年间每年的气温、气压误差均优于改进前，其中原gpt2w模型得到的5年平均气温、气压误差分别为
‑
0.45℃、2.05hpa，改进后的gptw模型得到的5年平均气温、气压误差分别为
‑
0.01℃、
‑
0.32hpa；改进后gpt2w模型得到的5年气温、气压最大正负误差均有所减小，说明改进后的gpt2w模型计算的气温、气压值更为稳定。
[0096]
表4 2012
‑
2016年中国陆地区域气温误差结果统计(gpt2w)
[0097][0098][0099]
表5 2012
‑
2016年中国陆地区域气温误差结果统计(gpt2w改正后)
[0100][0101]
表6 2012
‑
2016年中国陆地区域气压误差结果统计(gpt2w)
[0102][0103]
表7 2012
‑
2016年中国陆地区域气压误差结果统计(gpt2w改正后)
[0104][0105]
图2至图11为2012
‑
2016年gpt2w模型改进前后得到的气温、气压误差在中国陆地区域的分布图，由图可见在改进前大部分气象站点的气温年均误差值在
‑
1.5℃到1℃之间且存在地域性差异，其中，内蒙古、黑龙江、西藏东南部地区和新疆大部分地区产生的误差值较大，在改进后上述地区的气温误差值有显著减小且气温误差值浮动均在0℃上下；在改进前以东经105
°
经线作为分界线，小于东经105
°
经线的陆地区域各气象站点气压年均误差值较大，在改进后小于东经105
°
经线的陆地区域各气象站点气压年均误差值明显减小。
[0106]
图12至图13为gpt2w模型改进前后2016年安徽省各月份的气温、气压误差折线图，图中折线上的实习点分别表示气温月均误差和气压月均误差，由图可见每月改进后的gpt2w模型计算的气温、气压值精度有显著提升，且每月的气温、气压误差值浮动均在0刻度线附近。