一种基于SVR和VAR的地面风矢量临近预报方法

技术领域

本发明涉及机器学习领域，尤其涉及一种基于SVR和VAR的地面风矢量临近预报方法。

背景技术

风矢量（包括风速和风向）是地面气象站的重要观测要素，地面风矢量的变化能够显著影响社会生产和人民生活。例如在交通运输、尤其是航空领域，地面风的增强或骤变不仅影响飞机起降的频率，而且严重威胁生命财产安全。又如在能源领域，电力部门需要依据地面风的大小来调整风力发电计划，降低风机发电功率异常变化给输电网带来的潜在风险。因此，迫切需要准确率高、时间分辨率高的地面风矢量预报，为政府和民众采取防灾避灾措施提供重要支撑。

风矢量预报按照预报时效的长短，主要依靠数值模式或统计技术来实现。在短期（未来1~3天）预报方面，风矢量预报首要由数值模式的地面风预报场提供。而在临近预报方面，由于数值模式的局限性，比如计算耗时多、运行初始阶段模式因处于调整过程而预报误差较大等问题，基本通过统计技术来制作风矢量的预报。一方面，当前对地面风速的临近预报大多采用传统统计技术中的时间序列分析方法，例如自回归（AR）模型、滑动平均（MA）模型、自回归滑动平均（ARMA）模型等，主要考虑风在时间上的自相关性，即当前和过去一段时间内风速特征与未来一段时间内风速变化的关联性。在数学上，既由于风速的取值不能为负，又由于风速不满足高斯分布，不满足以上传统统计模型的基本假设，会给风速预报带来较大的误差；在物理上，以上常用方法中仅考虑了前期风速特征对后期风速变化的影响，而忽略了前期其他气象变量（如气压、温度等）与后期风速变化的联系，降低了风速预报的精度。另一方面，与地面风速预报相比，目前尚缺乏对地面风向的客观预报方法。为了解决以上问题，首先需要将风速预报扩展为风速和风向的预报，即将预报的地面风由标量还原为矢量，其次需要考虑地面其他气象要素对风矢量的作用。

在计量经济学领域广泛采用一种向量自回归（Vector Autoregressive model，简称VAR）模型，能够实现前期所有变量对后期所有变量的滞后回归。在VAR模型中，前、后期变量之间的关系是完全线性的，即前期不同变量对后期某一变量的影响是独立的、线性可叠加的。相比AR模型，可以将VAR模型引入地面风矢量的预报方法中，考虑地面多种气象要素对风矢量变化的调控作用。然而，地面气象要素间的关系以为非线性相关为主，传统的VAR模型无法捕捉。针对这一缺陷，需要将VAR模型中显式的线性方案改进为隐式的非线性方案，在模型估计前不预先设定变量间的关系方程，而是通过计算由模型自身给出。近年来，随着机器学习领域研究的不断深入和应用的不断扩展，机器学习内的相关技术已经能够考虑和描述变量间的非线性关系，例如SVR技术。SVR是将支持向量机应用在回归问题上，通过将原始变量映射到高维空间，实现将原始空间内的非线性关系用映射空间内的线性关系来表达，从而能够定量描述原始变量间的非线性相关性。此外，SVR因为不预先假设变量的统计分布，所以适用于非高斯分布的气象要素变量。

发明内容

为解决上述问题,本发明提出一种基于SVR和VAR的地面风矢量临近预报方法，本发明针对气象要素的特征将SVR技术进行优化，并引入和改进VAR模型，实现对地面风矢量的预报。

技术方案：为实现本发明的目的，本发明所采用的技术方案是：一种基于SVR和VAR的地面风矢量临近预报方法，其特征在于：包括以下步骤：

步骤1：地面多种气象要素的收集；

步骤2：地面风矢量的分解和质量控制；

所述步骤2中的质量控制公式为：

其中

步骤3：基于单个气象要素SVR模型的显著自变量提取；

所述步骤3中SVR模型估计的关系式为：

其中

其中

步骤4：基于所有显著自变量SVR和VAR的非线性模型估计；

所述步骤4中基于SVR技术的VAR模型估计的关系式为：

其中

步骤5：非线性模型的递推预报。

作为本发明进一步改进，所述步骤5中递推预报公式为：

其中

基于SVR和VAR的地面风矢量临近预报方法，具体步骤和特征如下：步骤1：地面多种气象要素的收集；步骤2：地面风矢量的分解和质量控制；步骤3：基于单个气象要素SVR的显著自变量提取；步骤4：基于所有显著自变量SVR和VAR的非线性模型估计；步骤5：非线性模型的递推预报。本发明提出一种基于SVR和VAR的地面风矢量临近预报方法，本发明针对气象要素的特征将SVR技术进行优化，并引入和改进VAR模型，实现对地面风矢量的预报。

附图说明

图1是方法架构和执行流程图。

具体实施方式

下面结合附图和实施例对本发明的技术方案作进一步的说明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。

本发明提出一种基于SVR和VAR的地面风矢量临近预报方法，本发明针对气象要素的特征将SVR技术进行优化，并引入和改进VAR模型，实现对地面风矢量的预报。

下面结合附图与具体实施方式对本发明做进一步详细描述。本发明工作流程图如图1所示，具体为：

首先，收集站点的历史观测数据，至少包括过去一个季度（90天）的地面多种气象要素观测值，观测间隔为10~15分钟。由于地面风的常规观测仅提供风速

分解后利用公式（1）对

接着，通过收集到的过去时段内的所有样本，针对每一种地面气象要素，分别利用公式（2）建立其与滞后的

然后，利用公式（4）分别建立滞后的

最后，利用公式（5）和公式（6）进行预报递推，其中的非线性模型参数由上一步估计得到。按照时间顺序进行递推，即先预报公式等号左边

以上所述，仅是本发明的较佳实施例而已，并非是对本发明作任何其他形式的限制，而依据本发明的技术实质所作的任何修改或等同变化，仍属于本发明所要求保护的范围。