CN107341134A - 一种对数值预报格点温度预报数据精细化处理的方法 - Google Patents

一种对数值预报格点温度预报数据精细化处理的方法 Download PDF

Info

Publication number
CN107341134A
CN107341134A CN201710500756.6A CN201710500756A CN107341134A CN 107341134 A CN107341134 A CN 107341134A CN 201710500756 A CN201710500756 A CN 201710500756A CN 107341134 A CN107341134 A CN 107341134A
Authority
CN
China
Prior art keywords
mrow
msub
value
temperature
day
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710500756.6A
Other languages
English (en)
Other versions
CN107341134B (zh
Inventor
郭铭博
高鸿飞
许方璐
杨林菲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Luoyang Meteorological Bureau
Original Assignee
Luoyang Meteorological Bureau
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Luoyang Meteorological Bureau filed Critical Luoyang Meteorological Bureau
Priority to CN201710500756.6A priority Critical patent/CN107341134B/zh
Publication of CN107341134A publication Critical patent/CN107341134A/zh
Application granted granted Critical
Publication of CN107341134B publication Critical patent/CN107341134B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Human Resources & Organizations (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Health & Medical Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Primary Health Care (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提供了一种对数值预报格点温度预报数据精细化处理的方法，具体如下：进行系统初始化后，获取实验数据，读取观测点及周围四个格点的经纬度坐标及观测点实况温度数据、周围四个格点温度预报值，进行误差分析和数据筛选，然后计算24节气每个节气的误差算术平均值和连续三天的同一预报时次误差的滑动平均值，并对所得数值进行反正切函数优化处理，再求得以上两个误差对预报结果影响的权重系数，然后按设定顺序依次计算得到各观测站点的误差滑动平均值、节气误差平均值及权重系数的集合，使用以上结果对各观测站点所代表县域内的数值预报格点温度数据进行修订计算，最后用修订数据对指定坐标点进行精细化插值计算和日最高温、日最低温的计算。

Description

一种对数值预报格点温度预报数据精细化处理的方法
技术领域
本发明涉及气象预报领域，尤其涉及一种对数值预报格点温度预报数据精细化处理的方法。
背景技术
目前，气象部门用于业务系统的预报方法主要有DMO方法、PP法、人工神经网络方法、MOS方法、相似预报方法、动力方法及卡尔曼滤波方法等。
传统的DMO方法就是通过插值把格点上的模式要素预报结果分析到具体的站点，得到站点上的要素预报，对于不是模式直接输出的要素，采用经验公式计算得到，其缺点在于对模式误差没有订正能力，预报精度完全依赖于模式，相对于形势场预报模式对要素预报的精度往往不是很高；
完全预报方法(PP法)是用历史资料中与预报对象同时间的实际气象参量作预报因子，建立统计方程，之后在假定数值预报的结果是“完全正确” (perfect)的前提下，用数值预报输出值代入到上述统计方程中，就可得到与预报时刻相对应的预报值；其缺点除含有统计关系造成的误差外，主要是无法考虑数值模式的预报误差，使最终的气温预报结果受到一定影响；
人工神经网络法是非线性方法的一种，它实际上是一个可自动实现2组变量间非线性映射关系的数据处理系统，其本质是优化计算中的梯度下降法，利用误差对于权、阈值的一阶导数信息，应用误差反传原理不断调整网络权值、阈值，使网络输出值与期望值之间的误差平方和达到最小或小于设定精度。
在当前人们对温度预报要求精度更高的情况下，传统的各种方法都已经不能适应目前的高精细化温度预报需求。
发明内容
为解决上述技术问题，本发明提供了一种对数值预报格点温度预报数据精细化处理的方法。
一种对数值预报格点温度预报数据精细化处理的方法，其中包括如下步骤：
1)、进行系统初始化；之后，进入步骤2)；
2)、获取实验数据，所述实验数据包括：1、地域范围内的EC细网格2米温度数值预报数据T；2、地域范围内各观测站的逐小时实况观测气温数据t；3、地域范围内各观测站的逐日最高气温数据tmax和最低气温数据tmin；将以上所述实验数据放入指定的系统数据存储区域；之后，进入步骤3)；
3)、观测站点的坐标经纬度值设为(I,J),对应EC细网格2米温度数值预报格点数据，该观测站点(I,J)的周围四个格点的经纬度值分别设为(I-1、J-1)、 (I+1、J-1)、(I+1、J+1)、(I-1、J+1)，所述观测站点逐小时实况气温数据的集合为 t1(I,J)，其周围四个格点位置的EC细网格2米温度数值预报数值的集合分别是 T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)，读取上述数据，对数据进行误差分析，并进行数据样本筛选，具体步骤如下：
31)、使用EC细网格2米温度数值预报数据T,对其预报时效内的每个预报数值集合按时间先后顺序，对当前观测站点周围四个格点位置的EC细网格2 米温度数值预报数值的集合T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)使用双线性二次插值法进行插值计算，从而得到和该观测站对应位置的每个对应时次的预报值集合T1(I,J),具体采用先纬向、后经向插值：
先在纬向I-1和I+1上进行线性一元一次插值，公式如下：
再在径向J上进行线性一元一次插值，公式如下：
式中：数据T1(I-1、J-1)、T1(I-1、J+1)、T1(I+1、J-1)、T1(I+1、J+1)分别表示当前观测站点(I,J)的周围最近相邻四个格点位置的EC细网格2米温度数值预报数值； T1(I-1、J)、T1(I+1、J)分别是纬度I-1、纬度I+1上的一次线性插值结果；T1(I、J)是进行上述插值运算后所得出的该观测站位置温度的预报值,之后进入步骤32)；
32)、利用步骤31)中计算得到的数据T1(I、J)，和其对应的实况数值t1(I,J)，计算两者误差ΔT1(I、J)＝t1(I、J)-T1(I、J)；之后，进入步骤33)；
33)、对所述观测站点所有预报时次的误差ΔT1(I、J)值进行筛选，只选取误差绝对值小于等于4的样本，对于|ΔT1(I、J)|大于4的按4取值，使 -4≤|ΔT1(I、J)|≤4；之后，进入步骤4)；
4)、利用步骤3)筛选后得到的误差ΔT1(I、J)，分别计算24节气每个节气的误差ΔT1(I、J)的算数平均值和连续三天的同一预报时次的误差ΔT1(I、J)的滑动平均值，并对计算所得数值使用反正切函数进行优化处理，具体计算方法为：同时进入步骤41)、42)：
41)、计算24节气每个节气的误差ΔT1(I、J)的算数平均值，具体方法如下：
其中n为同一个节气内的同一预报时次样本总数；
使用反正切函数对进行优化处理：
ST代表节气；
42)、计算连续三天的同一预报时次的ΔT(I、J)i滑动平均值，i是变量，表示一个具体日期的某一个预报时次，具体方法如下：
该式表示使用反正切函数对进行优化处理，
式中MA代表滑动平均；
5)、使用步骤3)得到的与所述观测站对应位置的每个对应时次的预报值集合T1(I,J)和步骤4)中得到的使用反正切函数优化后的24节气每个节气的误差ΔT1(I、J)的算数平均值使用反正切函数优化后的连续三天同一预报时次的误差ΔT1(I、J)的滑动平均值建立方程：
由此得到：
其中数据T1(I、J)为观测站点(I、J)的预报值，数据t1(I、J)为对应时间的该观测站点的气温实况值，根据步骤41)和步骤42)计算所得的数据和可以求得观测站点(I、J)的权重系数x1值，之后进入步骤6)；
6)、对于地域范围内其他观测站点按设定顺序重复步骤3)、步骤4)、步骤 5)，即可得到各观测站点的权重系数、滑动平均值、节气误差平均值的集合 之后进入步骤7)；
7)、根据步骤6)，计算得到各个观测站点节气的误差平均值滑动平均值和权重系数x值三个参数；各个观测站点均代表了各个县域内的气象条件，因此设定各个观测站所在县域内所有数值预报数据的格点均参照该观测站点的节气的误差平均值滑动平均值和权重系数x值三个参数；逐个读取区域内所有EC细网格格点位置最新时次的2米温度数值预报数据的值 TNEW，使用该格点所在县域观测站点的当天所在节气的误差平均值滑动平均值和权重系数x值三个参数，再进行计算：
TR为格点数值的修订值，将结果存入系统指定的存储区域，即完成对格点资料的订正；之后进入步骤8)；
8)、根据步骤7)得到的订正后的格点数值TR，然后对地域范围内指定坐标点(E,F)再次使用双线性二次插值法进行插值计算，从而得出该坐标点的预报值T(E、F)，同时根据该坐标点的预报值T(E、F)，做出地域范围内指定坐标点的日最高温度和最低温度的预报，具体步骤如下：
81)、在区域范围内指定坐标点(E,F)，其周边四个EC细网格格点的坐标分别为：(E-1,F-1)、(E-1,F+1)、(E+1,F+1)、(E+1,F-1)，使用EC细网格2米温度数值预报数据T，得到该指定坐标点(E,F)周边四个EC细网格格点的最新预报时次2 米温度预报数值的集合分别是T(E-1、F-1)、T(E-1、F+1)、T(E+1、F+1)、T(E+1、F-1)，使用双线性二次插值法进行插值计算，具体采用先纬向、后经向插值：
先在纬向E-1和E+1上进行线性一元一次插值，公式如下：
再在径向F上进行线性一元一次插值，公式如下：
式中：T(E-1、F)、T(E+1、F)分别是纬度E-1、纬度E+1上的一次线性插值结果；T(E、F)即是进行上述插值运算后所得出的指定坐标点(E,F)的温度预报值的集合；之后，同时进入步骤82)、83)；
82)、实现坐标点(E,F)最高温度预报的具体步骤如下：
821)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最高气温数据tmax建立最高温度公式：
tmax＝td14+α(td14-td08)
得到
其中d代表当天，08、14代表北京时的时次，td08、td14即为当天08时和 14时的实况观测温度值，tmax为当日最高温度，α为逐日最高温度的回归系数；
之后，进入步骤822)；
822)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，依此计算出历年来同一个节气里逐日最高温度回归系数α的算术平均值；之后，进入步骤823)；
823)、由步骤822)得到该观测站站点的逐日最高温度回归系数α的算术平均值各个观测站站点分别代表各个对应县域的气象条件，指定坐标点(E,F) 使用其对应的逐日最高温度回归系数α的算术平均值值，读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即分别为：当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时、14 时的预报值，T(E、F)H即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最高温度数值；之后，进入步骤824)；
824)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的14时、17时、当前日期前一天的20时三个时次数值：T(E、F)d14、T(E、F)d17、T(E、F)(d-1)20，与步骤823)所得到的当日的预计算最高温度数值T(E、F)H对比，选取最大值：
Tmax＝max{T(E、F)d14，T(E、F)d17，T(E、F)(d-1)20，T(E、F)H}；
Tmax即为当前日期坐标点(E、F)的最高温度,将结果存入系统指定的存储区域；
83)、实现坐标点(E,F)最低温度预报的具体步骤如下：
831)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最低气温数据tmin建立最低温度公式：
tmin＝td08+β(td14-td08)
得到
其中d代表当前日期，08、14代表北京时的时次，td08、td14即为当前日期08时和14时的实况观测温度值，β为逐日最低温度的回归系数；之后进入步骤832)；
832)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，计算出历年来同一个节气里逐日最低温度回归系数β的算术平均值之后进入步骤833)；
833)、由步骤832)得到该观测站站点的历年来同一个节气里逐日最低温度回归系数β的算术平均值各个观测站点都代表各个县域的气象条件，指定坐标点(E,F)也使用算术平均值读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即为当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时和14时的预报值，T(E、F)L即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最低温度数值；进入步骤834)；
834)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的02时、05时、20时三个时次数值：T(E、F)d02、 T(E、F)d05、T(E、F)d20，与步骤833)得到的当日的预计算最低温度值T(E、F)L对比，选取最小值：
Tmin＝min{T(E、F)d02，T(E、F)d05，T(E、F)d20，T(E、F)L}；Tmin即为当前日期坐标点(E、F)的最低温度,将结果存入系统指定的存储区域。
本发明所提供的一种对数值预报格点温度预报数据精细化处理的方法，是先通过插值把格点上的模式要素预报结果分析到具体的站点，得到站点上的要素预报，然后用历史资料中与预报对象同时间的实际气象要素值和前面得到的要素预报值做误差对比分析，建立统计方程；在预报误差影响因子的选取上，本发明方法还引入了滑动平均和节气的影响因子，对滑动平均影响因子和节气的影响因子两组变量通过非线性映射关系的数据处理，寻求滑动平均值因子和节气因子最优权重系数来进行建模，从而得到最佳的误差修订值；经试验证明，本发明方法对温度预报数据精细化的优化作用非常有效，具体优点如下：
(1)本发明方法可实现对数值预报产品中的格点温度要素预报产品进行细致订正，本发明方法还适用于当前气象部门使用的各类数值预报要素预报产品，并且数值预报模式升级后，本发明方法仍然适用；
(2)本发明方法对数值预报要素预报产品的修订，引入了滑动平均因子和不同节气的气候差异因子，用两者加权平均作为回归系数，这样既保证了结果相对于数值预报要素预报产品本身的准确性又体现了当地气候规律本身的客观性；
(3)本发明方法具有智能化自主学习功能，可以随着资料的进一步累积不断寻求最优参数自动优化预报结果，不需要人工干预自我升级。
附图说明
图1为本发明提供的一种对数值预报格点温度预报数据精细化处理的方法的流程图。
具体实施方式
本发明提供了一种对数值预报格点温度预报数据精细化处理的方法，如图 1所示流程图，具体包括如下步骤：
1)、进行系统初始化；之后，进入步骤2)；
2)、获取实验数据，所述实验数据包括：1、地域范围内的EC细网格2米温度数值预报数据T；2、地域范围内各观测站的逐小时实况观测气温数据t；3、地域范围内各观测站的逐日最高气温数据tmax和最低气温数据tmin；将以上所述实验数据放入指定的系统数据存储区域；之后，进入步骤3)；
3)、观测站点的坐标经纬度值设为(I,J),对应EC细网格2米温度数值预报格点数据，该观测站点(I,J)的周围四个格点的经纬度值分别设为(I-1、J-1)、 (I+1、J-1)、(I+1、J+1)、(I-1、J+1)，所述观测站点逐小时实况气温数据的集合为 t1(I,J)，其周围四个格点位置的EC细网格2米温度数值预报数值的集合分别是 T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)，读取上述数据，对数据进行误差分析，并进行数据样本筛选，具体步骤如下：
31)、使用EC细网格2米温度数值预报数据T,对其预报时效内的每个预报数值集合按时间先后顺序，对当前观测站点周围四个格点位置的EC细网格2 米温度数值预报数值的集合T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)使用双线性二次插值法进行插值计算，从而得到和该观测站对应位置的每个对应时次的预报值集合T1(I,J),具体采用先纬向、后经向插值：
先在纬向I-1和I+1上进行线性一元一次插值，公式如下：
再在径向J上进行线性一元一次插值，公式如下：
式中：数据T1(I-1、J-1)、T1(I-1、J+1)、T1(I+1、J-1)、T1(I+1、J+1)分别表示当前观测站点(I,J)的周围最近相邻四个格点位置的EC细网格2米温度数值预报数值； T1(I-1、J)、T1(I+1、J)分别是纬度I-1、纬度I+1上的一次线性插值结果；T1(I、J)是进行上述插值运算后所得出的该观测站位置温度的预报值,之后进入步骤32)；
32)、利用步骤31)中计算得到的数据T1(I、J)，和其对应的实况数值t1(I,J)，计算两者误差ΔT1(I、J)＝t1(I、J)-T1(I、J)；之后，进入步骤33)；
33)、对所述观测站点所有预报时次的误差ΔT1(I、J)值进行筛选，只选取误差绝对值小于等于4的样本，对于|ΔT1(I、J)|大于4的按4取值，使 -4≤|ΔT1(I、J)|≤4；之后，进入步骤4)；
4)、利用步骤3)筛选后得到的误差ΔT1(I、J)，分别计算24节气每个节气的误差ΔT1(I、J)的算数平均值和连续三天的同一预报时次的误差ΔT1(I、J)的滑动平均值，并对计算所得数值使用反正切函数进行优化处理，具体计算方法为：同时进入步骤41)、42)：
41)、计算24节气每个节气的误差ΔT1(I、J)的算数平均值，具体方法如下：
其中n为同一个节气内的同一预报时次样本总数；
使用反正切函数对进行优化处理：
ST代表节气；
42)、计算连续三天的同一预报时次的ΔT(I、J)i滑动平均值，i是变量，表示一个具体日期的某一个预报时次，具体方法如下：
该式表示使用反正切函数对进行优化处理，
式中MA代表滑动平均；
5)、使用步骤3)得到的与所述观测站对应位置的每个对应时次的预报值集合T1(I,J)和步骤4)中得到的使用反正切函数优化后的24节气每个节气的误差ΔT1(I、J)的算数平均值使用反正切函数优化后的连续三天同一预报时次的误差ΔT1(I、J)的滑动平均值建立方程：
由此得到：
其中数据T1(I、J)为观测站点(I、J)的预报值，数据t1(I、J)为对应时间的该观测站点的气温实况值，根据步骤41)和步骤42)计算所得的数据和可以求得观测站点(I、J)的权重系数x1值，之后进入步骤6)；
6)、对于地域范围内其他观测站点按设定顺序重复步骤3)、步骤4)、步骤 5)，即可得到各观测站点的权重系数、滑动平均值、节气误差平均值的集合 之后进入步骤7)；
7)、根据步骤6)，计算得到各个观测站点节气的误差平均值滑动平均值和权重系数x值三个参数；各个观测站点均代表了各个县域内的气象条件，因此设定各个观测站所在县域内所有数值预报数据的格点均参照该观测站点的节气的误差平均值滑动平均值和权重系数x值三个参数；逐个读取区域内所有EC细网格格点位置最新时次的2米温度数值预报数据的值 TNEW，使用该格点所在县域观测站点的当天所在节气的误差平均值滑动平均值和权重系数x值三个参数，再进行计算：
TR为格点数值的修订值，将结果存入系统指定的存储区域，即完成对格点资料的订正；之后进入步骤8)；
8)、根据步骤7)得到的订正后的格点数值TR，然后对地域范围内指定坐标点(E,F)再次使用双线性二次插值法进行插值计算，从而得出该坐标点的预报值T(E、F)，同时根据该坐标点的预报值T(E、F)，做出地域范围内指定坐标点的日最高温度和最低温度的预报，具体步骤如下：
81)、在区域范围内指定坐标点(E,F)，其周边四个EC细网格格点的坐标分别为：(E-1,F-1)、(E-1,F+1)、(E+1,F+1)、(E+1,F-1)，使用EC细网格2米温度数值预报数据T，得到该指定坐标点(E,F)周边四个EC细网格格点的最新预报时次2 米温度预报数值的集合分别是T(E-1、F-1)、T(E-1、F+1)、T(E+1、F+1)、T(E+1、F-1)，使用双线性二次插值法进行插值计算，具体采用先纬向、后经向插值：
先在纬向E-1和E+1上进行线性一元一次插值，公式如下：
再在径向F上进行线性一元一次插值，公式如下：
式中：T(E-1、F)、T(E+1、F)分别是纬度E-1、纬度E+1上的一次线性插值结果；T(E、F)即是进行上述插值运算后所得出的指定坐标点(E,F)的温度预报值的集合；之后，同时进入步骤82)、83)；
82)、实现坐标点(E,F)最高温度预报的具体步骤如下：
821)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最高气温数据tmax建立最高温度公式：
tmax＝td14+α(td14-td08)
得到
其中d代表当天，08、14代表北京时的时次，td08、td14即为当天08时和14时的实况观测温度值，tmax为当日最高温度，α为逐日最高温度的回归系数；
之后，进入步骤822)；
822)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，依此计算出历年来同一个节气里逐日最高温度回归系数α的算术平均值；之后，进入步骤823)；
这里“历年来同一个节气里逐日最高温度回归系数α的算术平均值”是指：节气是一个固定的时间点，两个节气之间约15天，用一个节气指代这15 天，每天一个回归系数α，一个节气15个回归系数α；历年来同一个节气指连续几年每一年的同一个节气，比如从2012年到2016年连续六年里每年的春分代表的15天，就是历年来春分这个节气，这五年来的春分这个节气的逐日回归系数的算数平均值是2012年的15个α、2013年的15个α、2014年的15个α、 2015年的15个α，2016年的15个α，总共75个α相加再除以75(5年乘以 15天)就是2012年到2016年来同一个节气里逐日回归系数的算数平均值；
823)、由步骤822)得到该观测站站点的逐日最高温度回归系数α的算术平均值各个观测站站点分别代表各个对应县域的气象条件，指定坐标点(E,F) 使用其对应的逐日最高温度回归系数α的算术平均值值，读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即分别为：当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时、14 时的预报值，T(E、F)H即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最高温度数值；之后，进入步骤824)；
824)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的14时、17时、当前日期前一天的20时三个时次数值：T(E、F)d14、T(E、F)d17、T(E、F)(d-1)20，与步骤823)所得到的当日的预计算最高温度数值T(E、F)H对比，选取最大值：
Tmax＝max{T(E、F)d14，T(E、F)d17，T(E、F)(d-1)20，T(E、F)H}；
Tmax即为当前日期坐标点(E、F)的最高温度,将结果存入系统指定的存储区域；
83)、实现坐标点(E,F)最低温度预报的具体步骤如下：
831)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最低气温数据tmin建立最低温度公式：
tmin＝td08+β(td14-td08)
得到
其中d代表当前日期，08、14代表北京时的时次，td08、td14即为当前日期08时和14时的实况观测温度值，β为逐日最低温度的回归系数；之后进入步骤832)；
832)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，计算出历年来同一个节气里逐日最低温度回归系数β的算术平均值之后进入步骤833)；
833)、由步骤832)得到该观测站站点的历年来同一个节气里逐日最低温度回归系数β的算术平均值各个观测站点都代表各个县域的气象条件，指定坐标点(E,F)也使用算术平均值读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即为当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时和14时的预报值，T(E、F)L即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最低温度数值；进入步骤834)；
834)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的02时、05时、20时三个时次数值：T(E、F)d02、 T(E、F)d05、T(E、F)d20，与步骤833)得到的当日的预计算最低温度值T(E、F)L对比，选取最小值：
Tmin＝min{T(E、F)d02，T(E、F)d05，T(E、F)d20，T(E、F)L}；Tmin即为当前日期坐标点(E、F)的最低温度,将结果存入系统指定的存储区域。
本发明所提供的一种对数值预报格点温度预报数据精细化处理的方法，是先通过插值把格点上的模式要素预报结果分析到具体的站点，得到站点上的要素预报，然后用历史资料中与预报对象同时间的实际气象要素值和前面得到的要素预报值做误差对比分析，建立统计方程；在预报误差影响因子的选取上，本发明方法还引入了滑动平均和节气的影响因子，对滑动平均影响因子和节气的影响因子两组变量通过非线性映射关系的数据处理，寻求滑动平均值因子和节气因子最优权重系数来进行建模，从而得到最佳的误差修订值。
经试验证明，本发明方法对温度预报数据精细化的优化作用非常有效，相比现有技术中的各种处理方法，本发明方法具有如下优点：
(1)本发明方法可实现对数值预报产品中的格点温度要素预报产品进行细致订正，本发明方法适用于当前气象部门使用的各类数值预报要素预报产品，并且数值预报模式升级后，本发明方法仍然适用；
(2)本发明方法对数值预报要素预报产品的修订，引入了滑动平均因子和不同节气的气候差异因子，用两者加权平均作为回归系数，这样既保证了结果相对于数值预报要素预报产品本身的准确性又体现了当地气候规律本身的客观性；
(3)本发明方法具有智能化自主学习功能，可以随着资料的进一步累积不断寻求最优参数自动优化预报结果，不需要人工干预自我升级。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内,因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (1)

1.一种对数值预报格点温度预报数据精细化处理的方法，其特征在于,包括如下步骤：
1)、进行系统初始化；之后，进入步骤2)；
2)、获取实验数据，所述实验数据包括：1、地域范围内的EC细网格2米温度数值预报数据T；2、地域范围内各观测站的逐小时实况观测气温数据t；3、地域范围内各观测站的逐日最高气温数据tmax和最低气温数据tmin；将以上所述实验数据放入指定的系统数据存储区域；之后，进入步骤3)；
3)、观测站点的坐标经纬度值设为(I,J),对应EC细网格2米温度数值预报格点数据，该观测站点(I,J)的周围四个格点的经纬度值分别设为(I-1、J-1)、(I+1、J-1)、(I+1、J+1)、(I-1、J+1)，所述观测站点逐小时实况气温数据的集合为t1(I,J)，其周围四个格点位置的EC细网格2米温度数值预报数值的集合分别是T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)，读取上述数据，对数据进行误差分析，并进行数据样本筛选，具体步骤如下：
31)、使用EC细网格2米温度数值预报数据T,对其预报时效内的每个预报数值集合按时间先后顺序，对当前观测站点周围四个格点位置的EC细网格2米温度数值预报数值的集合T1(I-1、J-1)、T1(I+1、J-1)、T1(I+1、J+1)、T1(I-1、J+1)使用双线性二次插值法进行插值计算，从而得到和该观测站对应位置的每个对应时次的预报值集合T1(I,J),具体采用先纬向、后经向插值：
先在纬向I-1和I+1上进行线性一元一次插值，公式如下：
<mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>J</mi> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>J</mi> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>J</mi> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>J</mi> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>J</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
再在径向J上进行线性一元一次插值，公式如下：
<mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>I</mi> <mo>-</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>I</mi> <mo>-</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>I</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>I</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> </mrow>
式中：数据T1(I-1、J-1)、T1(I-1、J+1)、T1(I+1、J-1)、T1(I+1、J+1)分别表示当前观测站点(I,J)的周围最近相邻四个格点位置的EC细网格2米温度数值预报数值；T1(I-1、J)、T1(I+1、J)分别是纬度I-1、纬度I+1上的一次线性插值结果；T1(I、J)是进行上述插值运算后所得出的该观测站位置温度的预报值,之后进入步骤32)；
32)、利用步骤31)中计算得到的数据T1(I、J)，和其对应的实况数值t1(I,J)，计算两者误差ΔT1(I、J)＝t1(I、J)-T1(I、J)；之后，进入步骤33)；
33)、对所述观测站点所有预报时次的误差ΔT1(I、J)值进行筛选，只选取误差绝对值小于等于4的样本，对于|ΔT1(I、J)|大于4的按4取值，使-4≤|ΔT1(I、J)|≤4；之后，进入步骤4)；
4)、利用步骤3)筛选后得到的误差ΔT1(I、J)，分别计算24节气每个节气的误差ΔT1(I、J)的算数平均值和连续三天的同一预报时次的误差ΔT1(I、J)的滑动平均值，并对计算所得数值使用反正切函数进行优化处理，具体计算方法为：同时进入步骤41)、42)：
41)、计算24节气每个节气的误差ΔT1(I、J)的算数平均值，具体方法如下：
<mrow> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>n</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mn>1</mn> <mi>n</mi> </munderover> <msub> <mi>&amp;Delta;T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> </mrow>
其中n为同一个节气内的同一预报时次样本总数；
使用反正切函数对进行优化处理：
ST代表节气；
42)、计算连续三天的同一预报时次的ΔT(I、J)i滑动平均值，i是变量，表示一个具体日期的某一个预报时次，具体方法如下：
该式表示使用反正切函数对进行优化处理，
式中MA代表滑动平均；
5)、使用步骤3)得到的与所述观测站对应位置的每个对应时次的预报值集合T1(I,J)和步骤4)中得到的使用反正切函数优化后的24节气每个节气的误差ΔT1(I、J)的算数平均值使用反正切函数优化后的连续三天同一预报时次的误差ΔT1(I、J)的滑动平均值建立方程：
由此得到：
<mrow> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mo>-</mo> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mn>1</mn> </msub> <msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mrow> <mi>S</mi> <mi>T</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mn>1</mn> </msub> <msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mn>1</mn> </msub> <msub> <mrow> <mo>(</mo> <mi>I</mi> <mo>,</mo> <mi>J</mi> <mo>)</mo> </mrow> <mrow> <mi>S</mi> <mi>T</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> </mrow>
其中数据T1(I、J)为观测站点(I、J)的预报值，数据t1(I、J)为对应时间的该观测站点的气温实况值，根据步骤41)和步骤42)计算所得的数据和可以求得观测站点(I、J)的权重系数x1值，之后进入步骤6)；
6)、对于地域范围内其他观测站点按设定顺序重复步骤3)、步骤4)、步骤5)，即可得到各观测站点的权重系数、滑动平均值、节气误差平均值的集合 ……；之后进入步骤7)；
7)、根据步骤6)，计算得到各个观测站点节气的误差平均值滑动平均值和权重系数x值三个参数；各个观测站点均代表了各个县域内的气象条件，因此设定各个观测站所在县域内所有数值预报数据的格点均参照该观测站点的节气的误差平均值滑动平均值和权重系数x值三个参数；逐个读取区域内所有EC细网格格点位置最新时次的2米温度数值预报数据的值TNEW，使用该格点所在县域观测站点的当天所在节气的误差平均值滑动平均值和权重系数x值三个参数，再进行计算：
<mrow> <msub> <mi>T</mi> <mi>R</mi> </msub> <mo>=</mo> <msub> <mi>T</mi> <mrow> <mi>N</mi> <mi>E</mi> <mi>W</mi> </mrow> </msub> <mo>+</mo> <mi>x</mi> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>M</mi> <mi>A</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mover> <mrow> <msub> <mi>&amp;Delta;T</mi> <mrow> <mi>S</mi> <mi>T</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow>
TR为格点数值的修订值，将结果存入系统指定的存储区域，即完成对格点资料的订正；之后进入步骤8)；
8)、根据步骤7)得到的订正后的格点数值TR，然后对地域范围内指定坐标点(E,F)再次使用双线性二次插值法进行插值计算，从而得出该坐标点的预报值T(E、F)，同时根据该坐标点的预报值T(E、F)，做出地域范围内指定坐标点的日最高温度和最低温度的预报，具体步骤如下：
81)、在区域范围内指定坐标点(E,F)，其周边四个EC细网格格点的坐标分别为：(E-1,F-1)、(E-1,F+1)、(E+1,F+1)、(E+1,F-1)，使用EC细网格2米温度数值预报数据T，得到该指定坐标点(E,F)周边四个EC细网格格点的最新预报时次2米温度预报数值的集合分别是T(E-1、F-1)、T(E-1、F+1)、T(E+1、F+1)、T(E+1、F-1)，使用双线性二次插值法进行插值计算，具体采用先纬向、后经向插值：
先在纬向E-1和E+1上进行线性一元一次插值，公式如下：
<mrow> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>F</mi> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>F</mi> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>F</mi> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>F</mi> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>F</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
再在径向F上进行线性一元一次插值，公式如下：
<mrow> <mi>T</mi> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>E</mi> <mo>-</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>E</mi> <mo>-</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>E</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mi>T</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> </mrow>
式中：T(E-1、F)、T(E+1、F)分别是纬度E-1、纬度E+1上的一次线性插值结果；T(E、F)即是进行上述插值运算后所得出的指定坐标点(E,F)的温度预报值的集合；之后，同时进入步骤82)、83)；
82)、实现坐标点(E,F)最高温度预报的具体步骤如下：
821)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最高气温数据tmax建立最高温度公式：
tmax＝td14+α(td14-td08)
得到
其中d代表当天，08、14代表北京时的时次，td08、td14即为当天08时和14时的实况观测温度值，tmax为当日最高温度，α为逐日最高温度的回归系数；
之后，进入步骤822)；
822)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，依此计算出历年来同一个节气里逐日最高温度回归系数α的算术平均值；之后，进入步骤823)；
823)、由步骤822)得到该观测站站点的逐日最高温度回归系数α的算术平均值各个观测站站点分别代表各个对应县域的气象条件，指定坐标点(E,F)使用其对应的逐日最高温度回归系数α的算术平均值值，读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
<mrow> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mi>H</mi> </msub> <mo>=</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>14</mn> </mrow> </msub> <mo>+</mo> <mover> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>S</mi> <mi>T</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mrow> <mi>E</mi> <mo>,</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>14</mn> </mrow> </msub> <mo>-</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mrow> <mi>E</mi> <mo>,</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>08</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> 3
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即分别为：当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时、14时的预报值，T(E、F)H即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最高温度数值；之后，进入步骤824)；
824)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的14时、17时、当前日期前一天的20时三个时次数值：T(E、F)d14、T(E、F)d17、T(E、F)(d-1)20，与步骤823)所得到的当日的预计算最高温度数值T(E、F)H对比，选取最大值：
Tmax＝max{T(E、F)d14，T(E、F)d17，T(E、F)(d-1)20，T(E、F)H}；
Tmax即为当前日期坐标点(E、F)的最高温度,将结果存入系统指定的存储区域；
83)、实现坐标点(E,F)最低温度预报的具体步骤如下：
831)、使用坐标点(E,F)所在县域内观测站的逐小时实况观测气温数据t和逐日最低气温数据tmin建立最低温度公式：
tmin＝td08+β(td14-td08)
得到
其中d代表当前日期，08、14代表北京时的时次，td08、td14即为当前日期08时和14时的实况观测温度值，β为逐日最低温度的回归系数；之后进入步骤832)；
832)、按照节气划分一年为24个时间段，读取对应节气时段的历史实况数据，计算出历年来同一个节气里逐日最低温度回归系数β的算术平均值之后进入步骤833)；
833)、由步骤832)得到该观测站站点的历年来同一个节气里逐日最低温度回归系数β的算术平均值各个观测站点都代表各个县域的气象条件，指定坐标点(E,F)也使用算术平均值读取步骤81)得到指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，代入公式：
<mrow> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mi>L</mi> </msub> <mo>=</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mi>E</mi> <mo>,</mo> <mi>F</mi> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>08</mn> </mrow> </msub> <mo>+</mo> <mover> <msub> <mi>&amp;beta;</mi> <mrow> <mi>S</mi> <mi>T</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mrow> <mi>E</mi> <mo>,</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>14</mn> </mrow> </msub> <mo>-</mo> <mi>T</mi> <msub> <mrow> <mo>(</mo> <mrow> <mi>E</mi> <mo>,</mo> <mi>F</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mi>d</mi> <mn>08</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
其中d代表当前日期，08、14代表北京时的时次，T(E、F)d08、T(E、F)d14即为当前日期指定坐标点(E,F)位置EC细网格2米温度资料在08时和14时的预报值，T(E、F)L即为计算出来的当前日期指定坐标点(E,F)位置当日的预计算最低温度数值；进入步骤834)；
834)、读取步骤81)得到的指定坐标点(E,F)的最新时次的EC细网格资料订正数值T(E、F)，选取当前日期的02时、05时、20时三个时次数值：T(E、F)d02、T(E、F)d05、T(E、F)d20，与步骤833)得到的当日的预计算最低温度值T(E、F)L对比，选取最小值：
Tmin＝min{T(E、F)d02，T(E、F)d05，T(E、F)d20，T(E、F)L}；Tmin即为当前日期坐标点(E、F)的最低温度,将结果存入系统指定的存储区域。
CN201710500756.6A 2017-06-27 2017-06-27 一种对数值预报格点温度预报数据精细化处理的方法 Expired - Fee Related CN107341134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710500756.6A CN107341134B (zh) 2017-06-27 2017-06-27 一种对数值预报格点温度预报数据精细化处理的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710500756.6A CN107341134B (zh) 2017-06-27 2017-06-27 一种对数值预报格点温度预报数据精细化处理的方法

Publications (2)

Publication Number Publication Date
CN107341134A true CN107341134A (zh) 2017-11-10
CN107341134B CN107341134B (zh) 2020-05-19

Family

ID=60221046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710500756.6A Expired - Fee Related CN107341134B (zh) 2017-06-27 2017-06-27 一种对数值预报格点温度预报数据精细化处理的方法

Country Status (1)

Country Link
CN (1) CN107341134B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108920485A (zh) * 2018-04-28 2018-11-30 深圳市雅码科技有限公司 一种格点气象服务产品精准制作分发系统及方法
CN109085664A (zh) * 2018-07-04 2018-12-25 山东省气象科学研究所 一种温度精细化预报偏差滑动订正方法
CN109934378A (zh) * 2018-12-30 2019-06-25 孙力勇 一种功率预测方法及装置
CN110068878A (zh) * 2019-04-22 2019-07-30 山东省气象科学研究所 一种气温智能网格最优集成预报方法
CN110837902A (zh) * 2018-08-15 2020-02-25 中国电力科学研究院有限公司 一种松弛逼近资料同化方法及系统
WO2020103677A1 (zh) * 2018-11-21 2020-05-28 国网青海省电力公司 数值天气预报的气象要素数据处理方法及装置
CN111352174A (zh) * 2020-03-20 2020-06-30 山东省气象科学研究所 一种数值天气预报及格点客观预报产品选优方法
CN112488382A (zh) * 2020-11-27 2021-03-12 清华大学 基于深度学习的enso预报方法
CN112748483A (zh) * 2020-12-24 2021-05-04 北京思湃德信息技术有限公司 一种基于深度学习的气温预报偏差订正方法及装置
CN113158578A (zh) * 2021-05-06 2021-07-23 北京邮电大学 基于机器学习的海洋低空波导预测方法
CN113239318A (zh) * 2021-05-17 2021-08-10 中国气象局乌鲁木齐沙漠气象研究所 一种区域数值预报模式的土壤湿度初值订正方法
CN114201561A (zh) * 2021-12-06 2022-03-18 宋英杰 一种节气气候时段的校正方法
CN116595806A (zh) * 2023-07-14 2023-08-15 江西师范大学 一种自适应温度数据补全方法
CN116805030A (zh) * 2023-05-15 2023-09-26 宁波市气象台 一种海浪浪高数值预报订正方法
CN117056661A (zh) * 2023-09-08 2023-11-14 华风气象传媒集团有限责任公司 一种气候三伏的确定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521403A (zh) * 2011-12-26 2012-06-27 南京成风大气信息技术有限公司 精细化气象信息服务系统
CN104298877A (zh) * 2014-10-13 2015-01-21 水利部交通运输部国家能源局南京水利科学研究院 一种可降低不确定性的气候变化情景修订方法
CN104849776A (zh) * 2014-12-08 2015-08-19 国家电网公司 一种结合动态修正的电网高低温精细化预警方法
CN105160425A (zh) * 2015-08-17 2015-12-16 国家电网公司 基于气象和地形因子的输变电设备覆冰灾害监测和预报方法
CN106202920A (zh) * 2016-07-08 2016-12-07 中国石油大学(华东) 一种单站海面气压的数值预报释用方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521403A (zh) * 2011-12-26 2012-06-27 南京成风大气信息技术有限公司 精细化气象信息服务系统
CN104298877A (zh) * 2014-10-13 2015-01-21 水利部交通运输部国家能源局南京水利科学研究院 一种可降低不确定性的气候变化情景修订方法
CN104849776A (zh) * 2014-12-08 2015-08-19 国家电网公司 一种结合动态修正的电网高低温精细化预警方法
CN105160425A (zh) * 2015-08-17 2015-12-16 国家电网公司 基于气象和地形因子的输变电设备覆冰灾害监测和预报方法
CN106202920A (zh) * 2016-07-08 2016-12-07 中国石油大学(华东) 一种单站海面气压的数值预报释用方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
邱学兴等: "乡镇精细化最高最低气温预报方法研究", 《气象与环境学报》 *
黄治勇等: "湖北省乡镇温度预报方法初探", 《气象》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108920485A (zh) * 2018-04-28 2018-11-30 深圳市雅码科技有限公司 一种格点气象服务产品精准制作分发系统及方法
CN108920485B (zh) * 2018-04-28 2022-03-25 深圳市雅码科技有限公司 一种格点气象服务产品精准制作分发系统及方法
CN109085664A (zh) * 2018-07-04 2018-12-25 山东省气象科学研究所 一种温度精细化预报偏差滑动订正方法
CN109085664B (zh) * 2018-07-04 2020-10-13 山东省气象科学研究所 一种温度精细化预报偏差滑动订正方法
CN110837902A (zh) * 2018-08-15 2020-02-25 中国电力科学研究院有限公司 一种松弛逼近资料同化方法及系统
CN110837902B (zh) * 2018-08-15 2024-07-26 中国电力科学研究院有限公司 一种松弛逼近资料同化方法及系统
WO2020103677A1 (zh) * 2018-11-21 2020-05-28 国网青海省电力公司 数值天气预报的气象要素数据处理方法及装置
CN109934378A (zh) * 2018-12-30 2019-06-25 孙力勇 一种功率预测方法及装置
CN110068878B (zh) * 2019-04-22 2021-04-09 山东省气象科学研究所 一种气温智能网格最优集成预报方法
CN110068878A (zh) * 2019-04-22 2019-07-30 山东省气象科学研究所 一种气温智能网格最优集成预报方法
CN111352174A (zh) * 2020-03-20 2020-06-30 山东省气象科学研究所 一种数值天气预报及格点客观预报产品选优方法
CN112488382A (zh) * 2020-11-27 2021-03-12 清华大学 基于深度学习的enso预报方法
CN112748483B (zh) * 2020-12-24 2022-08-19 北京思湃德信息技术有限公司 一种基于深度学习的气温预报偏差订正方法及装置
CN112748483A (zh) * 2020-12-24 2021-05-04 北京思湃德信息技术有限公司 一种基于深度学习的气温预报偏差订正方法及装置
CN113158578A (zh) * 2021-05-06 2021-07-23 北京邮电大学 基于机器学习的海洋低空波导预测方法
CN113239318A (zh) * 2021-05-17 2021-08-10 中国气象局乌鲁木齐沙漠气象研究所 一种区域数值预报模式的土壤湿度初值订正方法
CN113239318B (zh) * 2021-05-17 2022-12-20 中国气象局乌鲁木齐沙漠气象研究所 一种区域数值预报模式的土壤湿度初值订正方法
CN114201561A (zh) * 2021-12-06 2022-03-18 宋英杰 一种节气气候时段的校正方法
CN116805030A (zh) * 2023-05-15 2023-09-26 宁波市气象台 一种海浪浪高数值预报订正方法
CN116805030B (zh) * 2023-05-15 2024-02-13 宁波市气象台 一种海浪浪高数值预报订正方法
CN116595806A (zh) * 2023-07-14 2023-08-15 江西师范大学 一种自适应温度数据补全方法
CN116595806B (zh) * 2023-07-14 2023-10-10 江西师范大学 一种自适应温度数据补全方法
CN117056661A (zh) * 2023-09-08 2023-11-14 华风气象传媒集团有限责任公司 一种气候三伏的确定方法
CN117056661B (zh) * 2023-09-08 2024-06-04 华风气象传媒集团有限责任公司 一种气候三伏的确定方法

Also Published As

Publication number Publication date
CN107341134B (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN107341134A (zh) 一种对数值预报格点温度预报数据精细化处理的方法
CN111222698B (zh) 面向物联网的基于长短时记忆网络的积水水位预测方法
CN115507822B (zh) 水文循环变异驱动下的洪水风险预测方法
CN101480143B (zh) 一种预测灌区作物单产量的方法
Stathopoulos et al. Wind power prediction based on numerical and statistical models
CN110442937B (zh) 一种融合卫星遥感和机器学习技术的流域水文模拟方法
CN111027175B (zh) 基于耦合模型集成模拟的洪水对社会经济影响的评估方法
Lee et al. Generation of typical weather data using the ISO Test Reference Year (TRY) method for major cities of South Korea
CN110263866A (zh) 一种基于深度学习的电力用户负荷区间预测方法
Lima et al. Modeling and forecasting of Brazilian reservoir inflows via dynamic linear models
CN110555561A (zh) 一种中长期径流集合预报方法
Stauffer et al. Ensemble postprocessing of daily precipitation sums over complex terrain using censored high-resolution standardized anomalies
CN109242265B (zh) 基于误差平方和最小的城市需水量组合预测方法
KR101515003B1 (ko) 일사량 예측방법
CN105139093A (zh) 基于Boosting算法和支持向量机的洪水预报方法
CN103117546A (zh) 一种超短期风电功率滑动预测方法
CN112116149B (zh) 一种考虑预报不确定性关联演化特征的多站中长期径流滚动概率预测方法
CN104143043B (zh) 一种多功能气候数据获取方法
CN104849776A (zh) 一种结合动态修正的电网高低温精细化预警方法
CN106919645A (zh) 复杂地貌大景区的景点气象要素智能精细预测方法
CN115238947A (zh) 气候变化下旱涝急转事件的社会经济暴露度预估方法
CN115526413A (zh) 一种基于全连接神经网络日最高气温的预报方法
Suksamosorn et al. Post-processing of NWP forecasts using Kalman filtering with operational constraints for day-ahead solar power forecasting in Thailand
CN117172037B (zh) 一种分布式水文预报方法、装置、计算机设备及介质
Jose et al. Weather dependency of electricity demand: A case study in warm humid tropical climate

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200519

CF01 Termination of patent right due to non-payment of annual fee