CN111665575A - 一种基于统计动力的中长期降雨分级耦合预报方法及系统 - Google Patents

一种基于统计动力的中长期降雨分级耦合预报方法及系统 Download PDF

Info

Publication number
CN111665575A
CN111665575A CN202010463013.8A CN202010463013A CN111665575A CN 111665575 A CN111665575 A CN 111665575A CN 202010463013 A CN202010463013 A CN 202010463013A CN 111665575 A CN111665575 A CN 111665575A
Authority
CN
China
Prior art keywords
rainfall
forecasting
model
medium
forecast
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
CN202010463013.8A
Other languages
English (en)
Other versions
CN111665575B (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.)
Huazhong University of Science and Technology
Bureau of Hydrology Changjiang Water Resources Commission
Original Assignee
Huazhong University of Science and Technology
Bureau of Hydrology Changjiang Water Resources Commission
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 Huazhong University of Science and Technology, Bureau of Hydrology Changjiang Water Resources Commission filed Critical Huazhong University of Science and Technology
Priority to CN202010463013.8A priority Critical patent/CN111665575B/zh
Publication of CN111665575A publication Critical patent/CN111665575A/zh
Application granted granted Critical
Publication of CN111665575B publication Critical patent/CN111665575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Atmospheric Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Biomedical Technology (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Hydrology & Water Resources (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于统计动力的中长期降雨分级耦合预报方法及系统，属于降雨预报领域，包括：模型训练阶段：搜集目标区域和全球尺度的历史水文气象数据，并建立涵盖统计学模型和动力学模型的中长期降雨预报模型库，对于其中的每一个降雨预报模型，利用该模型对不同历史时间段的降雨量进行预报，并将预报结果划分为k个降雨量等级；基于神经网络建立各降雨量等级对应的耦合预报模型，并结合同期实测降雨量对其进行训练；降雨预报阶段：利用各降雨预报模型分别预报目标时间段的降雨量，确定相应的降雨量等级，并选取对应的耦合预报模型耦合各预报结果，得到最终预报结果。本发明能够提高中长期降雨预报的精度，并延长降雨预报的预见期。

Description

一种基于统计动力的中长期降雨分级耦合预报方法及系统
技术领域
本发明属于降雨预报领域，更具体地，涉及一种基于统计动力的中长期降雨分级耦合预报方法及系统。
背景技术
中长期水文预报是指根据前期气象水文要素，基于成因分析与数理统计的方法，对未来较长时间内的水文要素进行科学的预测。中长期降雨预报可以为中长期水文预报提供降雨数据，是中长期水文预报的基础。流域降雨过程是一种强烈依赖于大气环流、地表蒸发、全球气候因子变化等动力特性的非线性时变过程，涉及的变量众多，并受自然条件胁迫和人类活动干扰，因此流域降雨时空变异规律异常复杂，如何提高中长期降雨预报的精度，是中长期水文预报中一个亟待解决的难题。
现有的中长期降雨预报方法，主要利用特定的模型进行降雨预报，这些方法在一定程度上能够实现中长期降雨预报，但是降雨预报精度通常受制于输入的数据质量、选用的预报模型结构和优选的模型参数，具体表现在以下方面：①输入方面：以大气环流、海表温度、海冰、积雪等气象因子为降雨模型主要输入信息，降雨预报预见期受观测降雨时效性的限制；②模型方面：对于研究区域，通常采用一种降雨预报模型进行降雨过程模拟，忽略了不同降雨预报模型的流域适应性，预报结果容易出现较大误差，导致决策风险；③参数率定方面：优选模型参数的目标函数单一，难以全面反映流域内不同降雨过程的峰值、雨量以及演化过程差异；④预报精度方面：忽略了模型预报误差对预报模型的指导意义，未充分利用模型的互补性特征，进行预报结果的校正。
此外，大流域降雨预报，在空间维度存在地形复杂、空间差异性大的特点，在时间维度存在不确定性强的特点，且受上述一种或多种问题的协同制约，现行降雨预报难以满足流域水文预报、防灾减灾及水资源科学管理对降雨预报精度及预见期的需求。
发明内容
针对现有技术的缺陷和改进需求，本发明提供了一种基于统计动力的中长期降雨分级耦合预报方法及系统，其目的在于，提高中长期降雨预报的精度，并延长降雨预报的预见期。
为实现上述目的，按照本发明的一个方面，提供了一种基于统计动力的中长期降雨分级耦合预报方法，包括：
模型训练阶段：
搜集目标区域和全球尺度的历史水文气象数据，并建立涵盖统计学模型和动力学模型的中长期降雨预报模型库；中长期降雨预报模型库中每个降雨预报模型都用于预报预见期内任意一个时间段内的降雨量；
对于中长期降雨预报模型库中的每一个统计学模型，利用该模型根据所搜集的历史水文气象数据对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；对于中长期降雨预报模型库中的每一个动力学模型，利用该模型对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；其中，k≥1；
针对每一个降雨量等级，建立基于神经网络的耦合预报模型，用于耦合各降雨预报模型对于同一时间段的预报结果，作为相同时间段的最终预报结果；结合历史同期实测的降雨量，对各耦合预报模型模型进行训练，在训练完成后将各耦合预报模型作为各降雨量等级对应的中长期降雨耦合预报模型；
降雨预报阶段：
搜集目标区域和全球尺度的历史水文气象数据后，利用中长期降雨预报模型库中的各降雨预报模型分别预报目标时间段的降雨量，根据预报结果确定相应的降雨量等级，并选取对应的中长期降雨耦合预报模型耦合各预报结果，从而得到目标时间段内降雨量的最终预报结果。
本发明通过建立涵盖统计学模型和动力学模型的中长期降雨预报模型库，并进一步建立能够耦合多种降雨预报模型的预报结果的中长期降雨耦合预报模型，由此能够充分考虑不同降雨预报模型的流域适应性，利用模型的互补性特征，进行预报结果的校正，从而有效提高中长期降雨预报的预报精度；本发明通过聚类将降雨量划分为不同的降雨量等级，针对不同的降雨量等级，分别构建基于神经网络的中长期降雨耦合预报模型，由此能够充分考虑不同降雨量情形下各降雨预报模型的适应性特征，进一步提高中长期降雨预报的精度。基于高精度的中长期降雨预报结果，本发明可利用已获得的降雨预报结果对未来下一个时间段的降雨量进行预报，由此能够实现对未来多个时间段降雨量的预报，有效延长降雨预报的预见期。总的来说，本发明通过对降雨量进行等级划分，并针对各降雨量等级分别建立能够耦合多种降雨预报模型的预报结果的中长期降雨耦合预报模型，能够有效提高中长期降雨预报的精度，并有效延长降雨预报的预见期。
进一步地，对每一个降雨预报模型的预报结果都采用多种聚类算法进行聚类，并综合所有降雨预报模型对应的聚类结果确定每一个降雨预报模型对应的聚类结果，以使得不同降雨预报模型对同一时间段的降雨预报结果被划分到相同的降雨量等级。
本发明针对每个降雨预报模型都采用多种聚类算法对预报结果进行聚类，并综合所有预报模型对应的聚类结果进行筛选，由此能够使聚类效果稳定、分离度强且簇内紧凑度高。
进一步地，在任意时间段t的降雨量的最终预报结果为：
Wt *＝fi(Wt,1,Wt,2,…Wt,m,…Wt,M)；
其中，Wt,m为中长期降雨预报模型库中第m个降雨预报模型对时间段t的降雨预报结果，1≤m≤M，M为中长期降雨预报模型库中降雨预报模型的总数；i为各降雨预报模型对时间段t的降雨量预报结果所属的降雨量等级，fi()为第i个降雨量等级对应的中长期降耦合预报模型。
本发明利用神经网络模型耦合各降雨预报模型对于同一时间段的降雨预报结果，能够利用模型的互补性特征，进行预报结果的校正，提高预报精度。
进一步地，中长期降雨预报模型库中的降雨预报模型为：多元线性回归模型、随机森林模型、RegCM4模型和NCEP CFSv2；
其中，多元线性回归模型以多种大气-海洋-陆面过程物理因子为输入；随机森林以全球尺度下的多种大气-海洋-陆面遥相关因子为输入。
本发明所建立的降雨预报模型具体为多元线性回归模型、随机森林模型、RegCM4模型和NCEP CFSv2，其中，多元线性回归模型和随机森林模型为统计学模型，RegCM4模型和NCEP CFSv2为动力学模型；由于动力学模型反映了大气演化过程，因此，本发明所建立的中长期降雨耦合预报模型具有物理机制；由于多元线性回归模型和随机森林以目标区域内和全球尺度下的多种大气-海洋-陆面过程物理因子和遥相关因子为输入，本发明进行能够同时利用流域内以及全球尺度下的水文气象数据进行降雨预报，从而提高预报精度。
进一步地，多元线性回归模型为：
W＝αA+βO+γL+ε；
其中，W为多元线性回归模型的预报结果，A为大气因子集，O为海洋因子集，L为陆面因子集，α为大气因子的回归系数，β为海洋因子的回归系数，γ为陆面因子的回归系数，ε为残差项。
进一步地，本发明所提供的基于统计动力的中长期降雨分级耦合预报方法，在搜集目标区域的历史水文气象数据后，还包括：
计算所搜集的历史水文气象数据中各大气-海洋-陆面过程物理因子与目标区域内历史同期实测的降雨量之间的相关系数，并筛选出其中相关系数最大的topK个物理因子作为预报因子；
利用降雨预报模型进行降雨预报时，仅利用筛选出的预报因子进行预报；
其中，topK为正整数。
本发明仅筛选出历史水文气象数据中与历史实测降雨量之间的相关性较高的因子作为预报因子，用于后续的降雨预报，由此能够剔除无用的数据，防止过拟合，从而保证预报精度，还可以降低降雨预报的计算复杂度。
进一步地，相关系数的计算公式为：
Figure BDA0002511579310000051
其中，X表示因子序列，Y表示实测降雨量序列，Cov(X,Y)表示序列X和序列Y的协方差，Var(X)表示序列X的方差，Var(Y)表示序列Y的方差。
按照本发明的另一个方面，提供了一种基于统计动力的中长期降雨分级耦合预报系统，包括：计算机可读存储介质和处理器；
计算机可读存储介质用于存储可执行指令；
处理器用于读取计算机可读存储介质中存储的可执行指令，执行本发明提供的基于统计动力的中长期降雨分级耦合预报方法。
总体而言，通过本发明所构思的以上技术方案，能够取得以下有益效果：
(1)本发明通过对降雨量进行等级划分，并针对各降雨量等级分别建立能够耦合多种降雨预报模型的预报结果的中长期降雨耦合预报模型，能够有效提高中长期降雨预报的精度，并有效延长降雨预报的预见期。
(2)本发明针对每个降雨预报模型都采用多种聚类算法对预报结果进行聚类，并综合所有预报模型对应的聚类结果进行筛选，由此能够使聚类效果稳定、分离度强且簇内紧凑度高。
附图说明
图1为本发明实施例提供的基于统计动力的中长期降雨分级耦合预报方法示意图；
图2为本发明实施例提供的多元线性回归模型的降雨预报结果与同期实测降雨的对比示意图；
图3为本发明实施例提供的随机森林模型的降雨预报结果与同期实测降雨的对比示意图；
图4为本发明实施例提供的RegCM4模型的降雨预报结果与同期实测降雨的对比示意图；
图5为本发明实施例提供的NCEP CFSv2的降雨预报结果与同期实测降雨的对比示意图；
图6为本发明实施例提供的多种预报模型的预报结果及同期实测降雨的对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明中，本发明及附图中的术语“第一”、“第二”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。
为了提高中长期降雨预报的精度，并延长降雨预报的预见期，在本发明的一个实施例中，提供了一种基于统计动力的中长期降雨分级耦合预报方法，如图1所示，包括模型训练阶段和降雨预报阶段；
模型训练阶段包括：
搜集目标区域和全球尺度的历史水文气象数据，包括目标区域内的历史观测降雨、蒸发等水文气象数据，搜集大气环流指数、海温指数、太阳黑子、冷空气、台风数、北极海冰、青藏高原积雪等气候因子数据；建立涵盖统计学模型和动力学模型的中长期降雨预报模型库；中长期降雨预报模型库中每个降雨预报模型都用于预报预见期内任意一个时间段内的降雨量；
对于中长期降雨预报模型库中的每一个统计学模型，利用该模型根据所搜集的历史水文气象数据对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；对于中长期降雨预报模型库中的每一个动力学模型，利用该模型对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；其中，k≥1；
针对每一个降雨量等级，建立基于神经网络的耦合预报模型，用于耦合各降雨预报模型对于同一时间段的预报结果，作为相同时间段的最终预报结果；结合历史同期实测的降雨量，对各耦合预报模型模型进行训练，在训练完成后将各耦合预报模型作为各降雨量等级对应的中长期降雨耦合预报模型；
降雨预报阶段包括：
搜集目标区域和全球尺度的历史水文气象数据后，利用中长期降雨预报模型库中的各降雨预报模型分别预报目标时间段的降雨量，根据预报结果确定相应的降雨量等级，并选取对应的中长期降雨耦合预报模型耦合各预报结果，从而得到目标时间段内降雨量的最终预报结果。
在本实施例中，所建立的中长期降雨预报模型库中具体包括四个降雨预报模型，分别是多元线性回归模型、随机森林模型、RegCM4模型和NCEP CFSv2；
其中，多元线性回归模型以大气环流指数、海温指数、太阳黑子、冷空气、台风数、北极海冰、青藏高原积雪等多种大气-海洋-陆面过程物理因子为输入；随机森林以全球尺度下的亚洲纬向环流指数、南方涛动指数、亚洲区极涡强度指数等多种大气-海洋-陆面遥相关因子为输入；
在本实施例所建立的中长期降雨预报模型库中，多元线性回归模型和随机森林模型为统计学模型，RegCM4模型和NCEP CFSv2为动力学模型；在此所建立的四个降雨预报模型保证了中长期降雨预报模型库能够同时涵盖统计学模型和动力学模型，并且，这四个模型都已成功应用于中长期降雨预报，本实施例建立这四个降雨预报模型作为基础的降雨预报模型，能够保证模型的可靠性，保证了后续的耦合预报的可行性；
为了保证所建立的每一个降雨预报模型都能够可靠地应用于中长期降雨预报中，本实施例还利用历史水文气象数据分别对各降雨预报模型的预报精度进行了验证，元线性回归模型、随机森林模型、RegCM4模型和NCEP CFSv2的降雨预报结果与同期实测降雨的对比结果分别如图2～图5所示；多元回归和随机森林为统计模型，应用前要先利用历史数据拟合模型的参数，确定参数所采用的数据时间范围为率定期(1956-2006)，确定参数后要检查参数是否具有普遍性和合理性，就采用除率定期以外的数据进行验证(2007-2016)，如果计算精度与率定期接近，则证明所确定的参数合理，这段数据的时间范围为检验期；根据图2～图5所示的结果可知，本实施例所建立的四个降雨预报模型均具有较高的预报精度；
应当说明的是，本实施例中所建立的中长期降雨预报模型库，仅为本发明的一种可选的实施方式，不应理解为对本发明的唯一限定，只要能保证中长期降雨预报模型库同时涵盖统计学模型和动力学模型，且其中的各降雨预报模型均能可靠地应用于中长期降雨预报，本发明也可以建立不同的中长期降雨预报模型库。
本实施例中，在搜集目标区域的历史水文气象数据后，还包括：
计算所搜集的历史水文气象数据中各大气-海洋-陆面过程物理因子与目标区域内历史同期实测的降雨量之间的相关系数，并筛选出其中相关系数最大的topK个物理因子作为预报因子；其中，topK为正整数；
相关系数的计算公式为：
Figure BDA0002511579310000091
其中，X表示因子序列，Y表示实测降雨量序列，Cov(X,Y)表示序列X和序列Y的协方差，Var(X)表示序列X的方差，Var(Y)表示序列Y的方差；
Figure BDA0002511579310000092
和
Figure BDA0002511579310000093
分别为序列X和序列Y的均值；N表示所搜集的历史水文气象数据中，各种因子和历史降雨的序列长度；
利用降雨预报模型进行降雨预报时，仅利用筛选出的预报因子进行预报；
本实施例所建立的多元线性回归模型为：
W＝αA+βO+γL+ε；
其中，W＝[w1,w2,…wn]T为多元线性回归模型的预报结果，A＝[a1,a2,…an]T为大气因子集，O＝[o1,o2,…on]T为海洋因子集，L＝[l1,l2,…ln]T为陆面因子集，α＝[α1,α2,…,αn]T为大气因子的回归系数，β＝[β1,β2,…,βn]T为海洋因子的回归系数，γ＝[γ1,γ2,…,γn]T为陆面因子的回归系数，ε＝[ε1,ε2,…,εn]T为残差项，预报结果W、大气因子集A、海洋因子集O、陆面因子集L、回归系数α、β和γ以及残差项ε均为矩阵形式；n为多元线性回归方程的阶数，用于表示采用多少个历史值对未来降雨进行预测；
本实施例通过筛选出历史水文气象数据中与历史实测降雨量之间的相关性较高的因子作为预报因子，用于后续的降雨预报，能够剔除无用的数据，防止过拟合，从而保证预报精度，还可以降低降雨预报的计算复杂度。
本实施例中，对每一个降雨预报模型的预报结果都采用K-means、BRICH、分层聚类等多种聚类算法进行聚类，并综合所有降雨预报模型对应的聚类结果确定每一个降雨预报模型对应的聚类结果，以使得不同降雨预报模型对同一时间段的降雨预报结果被划分到相同的降雨量等级；
本实施例针对每个降雨预报模型都采用多种聚类算法对预报结果进行聚类，并综合所有预报模型对应的聚类结果进行筛选，由此能够使聚类效果稳定、分离度强且簇内紧凑度高；本实施例中，具体将降雨预报结果划分为三个等级；
对于每一个降雨预报模型，确定其聚类结果的同时，可提取各类预报降雨的最大值和最小值作为对应的降雨量等级的上、下约束边界，以便于后续确定任意一个预报结果所属的降雨量等级；对于任意第m(m＝1，2，3，4)个降雨预报模型，以Bup,i和Bdown,i分别表示第i(i＝1，2，3)个降雨量等级的上、下边界约束，则
Bup,i＝maxWi,j，
Bdown,i＝minWi,j，
其中，Wi,j为第i个降雨量等级中第j个降雨量值。
本实施例以不同模型参数适应性和多模型结构互补性为切入点，基于神经网络与多目标参数优化原理，耦合不同模型雨量等级适应性特点，在降雨预报中实现了多模型耦合互校，有效提高了流域中长期降雨的高精度预报；在本实施例中针对所划分的三个降雨量等级，分别基于神经网络模型建立了用于耦合各降雨预报模型对于同一时间段的预报结果的耦合预报模型；在模型训练过程中，将历史降雨预报结果及同期历史实测降雨划分为训练期和检验期，以不同预报降雨量级下不同模型模拟结果作为输入，以历史实测降雨作为模型输出，对耦合模型进行训练；进一步，采用检验期的预报降雨作为输入，对训练好的神经网络耦合模型进行检验，将耦合模型预报结果与实测降雨进行对比，计算耦合模型的预报结果精度；整个训练过程结束后，得到的耦合模型即为各降雨量等级对应的中长期降雨耦合预报模型，其表达式为：
Figure BDA0002511579310000111
其中，Wt,多元回归、Wt,随机森林、Wt,NCEP和Wt,regcm4分别表示利用多元线性回归模型、随机森林模型、NCEP CFSv2和RegCM4模型对时间段t的降雨预报结果，Wt,*表示这四个降雨预报结果；f1()、f2()和f3()分别表示第1、23个降雨量等级对应的中长期降耦合预报模型；
在本实施例中，可基于任意一种类型的神经网络建立耦合预报模型，例如，可以是BP神经网络，RBF神经网络，线性神经网络等常见的可用于数据分类的神经网络，也可以是其他的类型，在此将不作一一列举。
以汉江流域石泉站点以上区域的中长期降雨预报为例，采用上述方法实施例所提供的基于统计动力的中长期降雨分级耦合预报方法，进行降雨预报，降雨量等级划分结果如表1所示。
表1预报雨量精度等级划分结果
Figure BDA0002511579310000121
搜集上述四种降雨预报模型的同期历史降雨预报结果(2007年1月～2010年12月)，将上述历史同期降雨预报结果及同期历史实测降雨划分为训练期(2007年1月～2009年12月)和检验期(2009年1月～2010年12月)；
以不同预报降雨量级下不同模型模拟结果作为输入，以同期历史实测降雨作为模型输出，对基于神经网络的耦合预报模型进行训练；
采用检验期的预报降雨作为输入，对训练好的神经网络耦合预报模型进行检验，将耦合预报模型的预报结果与实测降雨进行对比，计算模型耦合结果的精度；各降雨预报模型的预报结果及耦合预报模型的预报结果如图6所示，相应的精度量化评价结果如表2所示。
表2各模型预报精度
Figure BDA0002511579310000122
由图6和表2可知，耦合预报模型预报结果的准确度要基本都高于任一单一降雨预报模型，不存在两项精度指标同时优于耦合预报模型的单一的降雨预报模型；算术平均法由于忽略了每种模型在不同条件下的适用性，导致其预报结果的精度相较于单一模型有所下降，加权平均法虽然能够在一定程度上改善降雨预报结果，但其精度水平仍低于耦合预报模型。由此可知，上述实施例所提供的基于统计动力的中长期降雨耦合分级预报方法能够有效提高中长期降雨预报的精度。
在本发明的另一个实施例中，提供了一种基于统计动力的中长期降雨分级耦合预报系统，包括：计算机可读存储介质和处理器；
计算机可读存储介质用于存储可执行指令；
处理器用于读取计算机可读存储介质中存储的可执行指令，执行上述方法实施例提供的基于统计动力的中长期降雨分级耦合预报方法。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于统计动力的中长期降雨分级耦合预报方法，其特征在于，包括：
模型训练阶段：
搜集目标区域和全球尺度的历史水文气象数据，并建立涵盖统计学模型和动力学模型的中长期降雨预报模型库；所述中长期降雨预报模型库中每个降雨预报模型都用于预报预见期内任意一个时间段内的降雨量；
对于所述中长期降雨预报模型库中的每一个统计学模型，利用该模型根据所搜集的历史水文气象数据对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；对于所述中长期降雨预报模型库中的每一个动力学模型，利用该模型对不同历史时间段的降雨量进行预报，并通过聚类将预报结果划分为k个降雨量等级；其中，k≥1；
针对每一个降雨量等级，建立基于神经网络的耦合预报模型，用于耦合各降雨预报模型对于同一时间段的预报结果，作为相同时间段的最终预报结果；结合历史同期实测的降雨量，对各耦合预报模型模型进行训练，在训练完成后将各耦合预报模型作为各降雨量等级对应的中长期降雨耦合预报模型；
降雨预报阶段：
搜集所述目标区域和全球尺度的历史水文气象数据后，利用所述中长期降雨预报模型库中的各降雨预报模型分别预报目标时间段的降雨量，根据预报结果确定相应的降雨量等级，并选取对应的中长期降雨耦合预报模型耦合各预报结果，从而得到所述目标时间段内降雨量的最终预报结果。
2.如权利要求1所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，对每一个降雨预报模型的预报结果都采用多种聚类算法进行聚类，并综合所有降雨预报模型对应的聚类结果确定每一个降雨预报模型对应的聚类结果，以使得不同降雨预报模型对同一时间段的降雨预报结果被划分到相同的降雨量等级。
对不同降雨预报模型的预报结果进行聚类时，保证对于同一时间段的降雨量预报结果被划分到相同的降雨量等级。
3.如权利要求1所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，在任意时间段t的降雨量的最终预报结果为：
Wt *＝fi(Wt,1,Wt,2,…Wt,m,…Wt,M)；
其中，Wt,m为所述中长期降雨预报模型库中第m个降雨预报模型对时间段t的降雨预报结果，1≤m≤M，M为所述中长期降雨预报模型库中降雨预报模型的总数；i为各降雨预报模型对时间段t的降雨量预报结果所属的降雨量等级，fi()为第i个降雨量等级对应的中长期降耦合预报模型。
4.如权利要求1-3任一项所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，所述中长期降雨预报模型库中的降雨预报模型为：多元线性回归模型、随机森林模型、RegCM4模型和NCEP CFSv2；
其中，多元线性回归模型以多种大气-海洋-陆面过程物理因子为输入；随机森林以全球尺度下的多种大气-海洋-陆面遥相关因子为输入。
5.如权利要求4所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，所述多元线性回归模型为：
W＝αA+βO+γL+ε；
其中，W为多元线性回归模型的预报结果，A为大气因子集，O为海洋因子集，L为陆面因子集，α为大气因子的回归系数，β为海洋因子的回归系数，γ为陆面因子的回归系数，ε为残差项。
6.如权利要求1-3所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，在搜集目标区域的历史水文气象数据后，还包括：
计算所搜集的历史水文气象数据中各大气-海洋-陆面过程物理因子与所述目标区域内历史同期实测的降雨量之间的相关系数，并筛选出其中相关系数最大的topK个物理因子作为预报因子；
利用降雨预报模型进行降雨预报时，仅利用筛选出的预报因子进行预报；
其中，topK为正整数。
7.如权利要求6所述的基于统计动力的中长期降雨分级耦合预报方法，其特征在于，相关系数的计算公式为：
Figure FDA0002511579300000031
其中，X表示因子序列，Y表示实测降雨量序列，Cov(X,Y)表示序列X和序列Y的协方差，Var(X)表示序列X的方差，Var(Y)表示序列Y的方差。
8.一种基于统计动力的中长期降雨分级耦合预报系统，其特征在于，包括：计算机可读存储介质和处理器；
所述计算机可读存储介质用于存储可执行指令；
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令，执行权利要求1-7任一项所述的基于统计动力的中长期降雨分级耦合预报方法。
CN202010463013.8A 2020-05-27 2020-05-27 一种基于统计动力的中长期降雨分级耦合预报方法及系统 Active CN111665575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010463013.8A CN111665575B (zh) 2020-05-27 2020-05-27 一种基于统计动力的中长期降雨分级耦合预报方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010463013.8A CN111665575B (zh) 2020-05-27 2020-05-27 一种基于统计动力的中长期降雨分级耦合预报方法及系统

Publications (2)

Publication Number Publication Date
CN111665575A true CN111665575A (zh) 2020-09-15
CN111665575B CN111665575B (zh) 2021-01-05

Family

ID=72384874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010463013.8A Active CN111665575B (zh) 2020-05-27 2020-05-27 一种基于统计动力的中长期降雨分级耦合预报方法及系统

Country Status (1)

Country Link
CN (1) CN111665575B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112232554A (zh) * 2020-09-28 2021-01-15 中国科学院空天信息创新研究院 一种基于bp神经网络的局域短临降水预报模型的构建方法
CN112800634A (zh) * 2021-04-07 2021-05-14 水利部交通运输部国家能源局南京水利科学研究院 耦合干湿状态辨识与多源信息融合的降水估计方法及系统
CN113435646A (zh) * 2021-06-25 2021-09-24 中国水利水电科学研究院 基于年降雨长期预报方法的矿区生态需水量计算方法
CN114648181A (zh) * 2022-05-24 2022-06-21 国能大渡河大数据服务有限公司 一种基于机器学习的降水预报修正方法及系统
CN115356789A (zh) * 2022-10-08 2022-11-18 南京气象科技创新研究院 一种梅雨期短时强降水分级预警方法
WO2023060466A1 (zh) * 2021-10-13 2023-04-20 中山大学 一种基于确定性系数的降水预报能力与遥相关作用关联方法
CN117111181A (zh) * 2023-09-05 2023-11-24 浙江省气象台 一种短时强降水概率预报方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08106448A (ja) * 1994-10-04 1996-04-23 Nippon Telegr & Teleph Corp <Ntt> 気象予測装置
KR101563244B1 (ko) * 2014-02-06 2015-10-26 부산대학교 산학협력단 장마 강수의 물리 통계 예측 시스템 및 방법
CN113379109B (zh) * 2017-11-21 2022-04-19 中国水利水电科学研究院 一种基于预测模型自适应的径流预报方法
CN109492825A (zh) * 2018-11-26 2019-03-19 中国水利水电科学研究院 基于互信息和主成分分析筛选因子的中长期预报方法
CN109711617B (zh) * 2018-12-25 2022-08-19 河海大学 一种基于blstm深度学习的中长期径流预测方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112232554A (zh) * 2020-09-28 2021-01-15 中国科学院空天信息创新研究院 一种基于bp神经网络的局域短临降水预报模型的构建方法
CN112800634A (zh) * 2021-04-07 2021-05-14 水利部交通运输部国家能源局南京水利科学研究院 耦合干湿状态辨识与多源信息融合的降水估计方法及系统
CN112800634B (zh) * 2021-04-07 2021-06-25 水利部交通运输部国家能源局南京水利科学研究院 耦合干湿状态辨识与多源信息融合的降水估计方法及系统
CN113435646A (zh) * 2021-06-25 2021-09-24 中国水利水电科学研究院 基于年降雨长期预报方法的矿区生态需水量计算方法
WO2023060466A1 (zh) * 2021-10-13 2023-04-20 中山大学 一种基于确定性系数的降水预报能力与遥相关作用关联方法
CN114648181A (zh) * 2022-05-24 2022-06-21 国能大渡河大数据服务有限公司 一种基于机器学习的降水预报修正方法及系统
CN115356789A (zh) * 2022-10-08 2022-11-18 南京气象科技创新研究院 一种梅雨期短时强降水分级预警方法
CN117111181A (zh) * 2023-09-05 2023-11-24 浙江省气象台 一种短时强降水概率预报方法及系统
CN117111181B (zh) * 2023-09-05 2024-04-09 浙江省气象台 一种短时强降水概率预报方法及系统

Also Published As

Publication number Publication date
CN111665575B (zh) 2021-01-05

Similar Documents

Publication Publication Date Title
CN111665575B (zh) 一种基于统计动力的中长期降雨分级耦合预报方法及系统
Roy et al. Predicting land cover change in a Mediterranean catchment at different time scales
CN113537600B (zh) 一种全过程耦合机器学习的中长期降水预报建模方法
CN104764868B (zh) 一种基于地理加权回归的土壤有机碳预测方法
CN114254561A (zh) 一种内涝预测方法、系统及存储介质
Körner et al. Introducing Gradient Boosting as a universal gap filling tool for meteorological time series
CN111767517B (zh) 一种应用于洪水预测的BiGRU多步预测方法、系统及存储介质
CN105974495B (zh) 利用分类拟合法预判目标区域未来平均云量的方法
CN111260111A (zh) 基于气象大数据的径流预报改进方法
Olsson et al. Statistical atmospheric downscaling of short-term extreme rainfall by neural networks
CN113361742A (zh) 一种基于水文模拟的区域综合干旱识别方法
CN114254802B (zh) 气候变化驱动下植被覆盖时空变化的预测方法
CN116805439A (zh) 一种基于人工智能和大气环流机制的干旱预报方法及系统
CN111445085A (zh) 一种考虑中大型水库工程蓄水影响的中长期径流预报方法
Liu et al. Identify optimal predictors of statistical downscaling of summer daily precipitation in China from three-dimensional large-scale variables
CN114372631A (zh) 一种基于小样本学习和lstm的缺资料地区径流预测方法
CN117526274A (zh) 极端气候下新能源功率预测方法、电子设备和存储介质
CN116663404A (zh) 一种耦合人工智能和贝叶斯理论的洪水预报方法及系统
CN117114176A (zh) 基于数据分析和机器学习的土地利用变化预测方法及系统
Lu et al. Uncertainty quantification of machine learning models to improve streamflow prediction under changing climate and environmental conditions
Chao et al. A spatio-temporal neural network learning system for city-scale carbon storage capacity estimating
CN115860165A (zh) 一种考虑初损的神经网络流域降雨径流预报方法及系统
CN112380778A (zh) 一种基于海温的气象干旱预报方法
Neykov et al. Linking atmospheric circulation to daily precipitation patterns over the territory of Bulgaria
Zhang et al. Calibration and uncertainty analysis of a hydrological model based on cuckoo search and the M-GLUE method

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