首页 / 专利库 / 地球科学 / 水文学 / 一种串并联耦合的多模型水文预报方法

一种串并联耦合的多模型文预报方法

阅读:833发布:2020-05-31

专利汇可以提供一种串并联耦合的多模型文预报方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及 水 文学中的水文预报领域,公开了一种串并联耦合的多模型水文预报方法。充分发挥了实时校正与组合预报两种传统方法的优势,建立了一种串并联耦合的水文预报模型,按耦合方式的不同,共分为三种形式,即先串后并法、先并后串法以及一体化耦合方法。先串后并方法是先对水文模型的计算结果进行 串联 校正,而后再对多个模型的校正结果进行并联校正;先并后串方法,是先对多个水文模型的计算结果进行并联校正,再对组合预报的结果进行串联校正;一体化耦合方法是基于预报误差最小原则,建立串并联耦合校正结果的目标函数,同时率定串联校正和并联校正的参数。本发明能最大限度地减少预报误差,提高预报 精度 ,可广泛应用于水文实时预报中。,下面是一种串并联耦合的多模型文预报方法专利的具体信息内容。

1.一种串并联耦合的多模型文预报方法,其特征在于,包括如下步骤:
步骤1:选择预报断面,获取该断面历史时期实测降雨和流量资料,将资料分为率定期和检验期两部分;采用率定期的数据,率定水文模型的参数,采用检验期的数据,检验水文模型的效果;建立多个水文预报模型,分别预报该断面的径流量;
步骤2:建立串并联耦合的水文预报模型,按耦合方式的不同,共分为三种形式,即先串后并法、先并后串法以及一体化耦合模型;
步骤3:选取平均绝对误差、平均相对误差、确定性系数和均方根误差等常用精度评定指标,对三种形式的串并联耦合模型的结果进行精度评定。
2.如权利要求1所述的方法,其特征在于,所述步骤2中的先串后并法具体为:
①首先,对单个水文预报模型的结果进行校正,完成串联预报的过程;
基于实测和步骤1中各水文模型的预报流量,得到各预报模型的预报误差序列;依据预报误差序列的自相关特性,建立误差自回归校正模型;具体如下:
令实测序列为Qt=(Q1,Q2,...,Qn),n为序列的长度,预测序列为
预测序列与实测序列的差值为误差序列,记为et=(e1,e2,...,en);误差自回归模型的数学表达式为:
et=θ1et-1+θ2et-2+…+θqet-q+ξ
式中,θi(i=1,2,…,q)为自回归模型的参数;采用上式,建立q阶自回归模型AR(q);
q阶自回归模型AR(q)的AIC(q)值为:
选择在检验期AIC值最小的q阶自回归模型参与以下计算:通过建立的自回归模型预测t时刻的预报误差et,采用得到的预报误差值对预报结果进行校正;针对每个水文预报模型,利用上述方法,得到各个预报模型在该断面预报流量的串联预报结果;
②其次,综合每个单一模型串联预报的结果,完成并联校正过程;
基于m个水文模型串联预报序列 建立并联校正模型,如下所示:
式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1;由此得到了先串后并耦合模型的计算结果。
3.如权利要求1所述的方法,其特征在于,所述步骤2中的先并后串法具体为:
①基于m个水文模型的预报序列Q1t,Q2t,...,Qmt,建立并联校正模型,如下所示:
Ft=ω1Q1t+ω2Q2t+…+ωmQmt
式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1;采用最小二乘法计算并联校正中各模型的权重;
②依据并联校正的结果,得到并联校正的误差序列;利用自回归模型原理,建立误差自回归串联校正模型AR(q),基于AIC准则确定模型的阶数,对多个预报模型的并联校正结果进行串联校正;由此得到了先并后串耦合模型的计算结果。
4.如权利要求1所述的方法,其特征在于,所述步骤2中的一体化耦合方法具体为:
利用断面实测径流量和各个水文预报模型模拟预报的径流量,建立一体化耦合模型;
记水文模型个数为m,实测序列为Qt,各个水文模型的预报结果为Q1t,Q2t,...,Qmt;对各个模型采用误差自回归模型预测t时刻的误差序列et,表达式如下:
式中,(α11,α21),(α12,α22),…,(α1m,α2m)为回归系数;采用上述预报误差值对预报结果进行校正,得到各模型串联校正后的流量序列
在各模型串联预报结果的基础上,利用最小二乘法,进行多个模型的并联校正,并联校正的结果为Ft:
式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1;
采用整体优化的方法同时求得到串联校正与并联校正的参数;令并联校正结果的误差为Ft-Qt,其误差的期望记为E(Ft-Qt)2;为了使误差的期望值E(Ft-Qt)2最小,得到如下目标函数:
式中,其参数包含回归系数(α11,α21),(α12,α22),…,(α1m,α2m)以及权重系数ω1,ω2,…,ωm两部分,共3m个待定参数;
引入粒子群智能优化算法同时求解这些参数;得到这些参数后,先采用串联校正方法,对各水文模型的预报流量进行串联校正,得到串联校正结果;而后,利用并联校正方法,对各水文模型实时校正后的结果进行并联校正;由此得到了一体化耦合方法的计算结果。
5.如权利要求1-4任一项所述的方法,其特征在于,所述步骤3中的精度评定指标包括以下一种或多种:平均绝对误差、平均相对误差、确定性系数以及均方根误差。
6.如权利要求5所述的方法,其特征在于,所述平均绝对误差MAE如下式所示:
式中,Qi为实测值, 为预报值,n为预报序列长度。
7.如权利要求5所述的方法,其特征在于,所述平均相对误差MRE如下式所示:
式中,Qi为实测值, 为预报值,n为预报序列长度。
8.如权利要求5所述的方法,其特征在于,所述确定性系数DC如下式所示:
式中, 为实测值在预报序列长度内的平均值,Qi为实测值, 为预报值,n为预报序列长度。
9.如权利要求5所述的方法,其特征在于,所述均方根误差RMSE如下式所示:
式中,Qi为实测值, 为预报值,n为预报序列长度。
10.如权利要求1或2所述的方法,其特征在于,所述步骤1中的水文预报模型包括如下的一种或多种:新安江模型、水箱模型以及API模型。

说明书全文

一种串并联耦合的多模型文预报方法

技术领域

[0001] 本发明属于水文学中的水文预报领域,涉及一种串并联耦合的多模型水文预报方法。该方法能够为水利管理部实现高精度的水文预报提供理论依据与技术支撑

背景技术

[0002] 水文预报通过对未来水文情势(如洪峰流量)做出科学预测,特别是对灾害性水文现象做出准确预报,从而实现防洪减灾以及水资源的合理开发利用。提高水文预报的精度是水文预报工作的重要内容,对防洪减灾、保护人民生命财产安全、充分发挥水利工程效益改善生态环境等起着至关重要的作用。
[0003] 现行的提高水文预报精度的方法主要有实时校正和组合预报两种方法。实时校正方法利用了预报误差序列自身的相关性特征,预测未来时刻的误差值,从而实现对水文预报结果的实时校正(如专利CN201010106038.9—洪水预报系统中实时校正模型的优选方法)。组合预报方法是在当前应用广泛的各种水文预报模型的基础上,利用模型之间的互补性,融入了模型加权法的思想,选取对应用流域较为适用的几种水文模型建立组合预报方案,进行水文预报(如专利CN200910234628.7—一种不同机制水文模型组合的水文预报方法)。

发明内容

[0004] 现行方法未将实时校正与组合预报两种方法紧密结合起来,建立更为系统和科学的水文预报校正模型。本发明提出了实时校正与组合预报耦合的水文预报方法,即一种串并联耦合的多模型水文预报方法,以便充分发挥实时校正与组合预报的优势,从而最大限度地提高水文预报的精度和水平。
[0005] 本发明在现有实时校正与组合预报模型的基础上,创造性地结合了二者的优势,研发了一种串并联耦合的多模型水文预报方法,依据耦合方式的不同,共提出三种有效形式,即先串后并法、先并后串法,以及一体化耦合方法。本发明最大限度地减少了预报误差、提高了预报精度,从而实现了高精度的水文实时预报。
[0006] 为解决上述技术问题,本发明采用如下技术方案:
[0007] 步骤1:选择预报断面,获取该断面历史时期实测降雨和流量资料,将资料分为率定期和检验期两部分;采用率定期的数据,率定水文模型的参数,采用检验期的数据,检验水文模型的效果。建立多个水文预报模型,如新安江模型、水箱模型以及API(Antecedent Precipitation Index)模型等,分别预报该断面的径流量。
[0008] 步骤2:建立串并联耦合的水文预报模型。
[0009] 本发明研发了一种串并联耦合的多模型水文预报方法,按耦合方式的不同,共分为三种形式,即先串后并法、先并后串法以及一体化耦合模型。
[0010] (1)先串后并法
[0011] ①首先,对单个水文预报模型的结果进行校正,完成串联预报的过程。
[0012] 基于实测和步骤1中各水文模型的预报流量,得到各预报模型的预报误差序列。依据预报误差序列的自相关特性,建立误差自回归校正模型。具体如下:
[0013] 令 实 测 序 列 为 Qt=(Q 1,Q2,...,Qn),n为 序 列的 长 度,预 测 序 列 为预测序列与实测序列的差值为误差序列,记为et=(e1,e2,...,en)。误差自回归模型的数学表达式为:
[0014] et=θ1et-1+θ2et-2+…+θqet-q+ξ
[0015] 式中,θi(i=1,2,…,q)为自回归模型的参数。采用上式,建立q阶自回归模型AR(q)。为了使所选取的自回归模型对数据的拟合效果最好,通过AIC准则来确定模型的阶数。q阶自回归模型AR(q)的AIC(q)值为:
[0016]
[0017] 选择在检验期AIC(Akaike information criterion)值最小的q阶自回归模型参与以下计算。通过建立的自回归模型预测t时刻的预报误差et,采用得到的预报误差值对预报结果进行校正。针对每个水文预报模型,利用上述方法,得到各个预报模型在该断面预报流量的串联预报结果。
[0018] ②其次,综合每个单一模型串联预报的结果,完成并联校正过程。
[0019] 基于m个水文模型串联预报序列 建立并联校正模型,如下所示:
[0020]
[0021] 式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1。采用最小二乘法计算并联预报的权重。由此得到了先串后并耦合模型的计算结果。
[0022] (2)先并后串法
[0023] ①基于m个水文模型的预报序列Q1t,Q2t,...,Qmt,建立并联校正模型,如下所示:
[0024] Ft=ω1Q1t+ω2Q2t+…+ωmQmt
[0025] 式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1。采用最小二乘法计算并联校正中各模型的权重。
[0026] ②依据并联校正的结果,得到并联校正的误差序列。利用上述自回归模型原理,建立误差自回归串联校正模型AR(q),基于AIC准则确定模型的阶数,对多个预报模型的并联校正结果进行串联校正。由此得到了先并后串耦合模型的计算结果。
[0027] (3)一体化耦合方法
[0028] 利用断面实测径流量和各个水文预报模型模拟预报的径流量,建立一体化耦合模型。记水文模型个数为m,实测序列为Qt,各个水文模型的预报结果为Q1t,Q2t,...,Qmt。对各个模型采用误差自回归模型预测t时刻的误差序列et,表达式如下:
[0029]
[0030] 式中,(α11,α21),(α12,α22),…,(α1m,α2m)为回归系数。采用上述预报误差值对预报结果进行校正,得到各模型串联校正后的流量序列
[0031] 在各模型串联预报结果的基础上,利用最小二乘法,进行多个模型的并联校正,并联校正的结果为Ft:
[0032]
[0033] 式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足条件ω1+ω2+…+ωm=1。
[0034] 采用整体优化的方法(一体化方法)同时求得到串联校正与并联校正的参数。令并联校正结果的误差为Ft-Qt,其误差的期望记为E(Ft-Qt)2。为了使误差的期望值E(Ft-Qt)2最小,得到如下目标函数:
[0035]
[0036] 式中,其参数包含回归系数(α11,α21),(α12,α22),…,(α1m,α2m)以及权重系数ω1,ω2,…,ωm两部分,共3m个待定参数。
[0037] 引入粒子群智能优化算法(Particle Swarm Optimizer,PSO)同时求解这些参数。得到这些参数后,先采用串联校正方法,对各水文模型的预报流量进行串联校正,得到串联校正结果;而后,利用并联校正方法,对各水文模型实时校正后的结果进行并联校正。由此得到了一体化耦合方法的计算结果。
[0038] 步骤3:选取平均绝对误差、平均相对误差、确定性系数和均方根误差等常用精度评定指标,对三种形式的串并联耦合模型的结果进行精度评定。
[0039] 与现有技术相比,本发明具有以下优点和效果:
[0040] (1)已有技术只采用实时校正技术或组合预报方法提高水文预报的精度,而未将二者结合起来考虑。本发明在综合实时校正与组合预报优势的基础上,提出了一种串并联耦合的多模型水文预报方法,校正结果要明显的优于单纯实时校正或组合预报的结果,最大限度地减少了预报误差,提高了水文预报的精度和水平。
[0041] (2)已有技术采用最小二乘法来计算串联校正与并联校正的参数,这种参数估计方法应用到耦合模型时,缺乏对校正效果的整体把握。本发明利用粒子群优化算法,同时率定串并联耦合校正的参数,即一体化耦合模型。与先串后并和先并后串法相比,一体化耦合模型的预报精度最高。附图说明
[0042] 图1为本发明串并联耦合的多模型水文预报方法流程图

具体实施方式

[0043] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0044] 本发明提出了一种串并联耦合的多模型水文预报方法,依据耦合方式的不同,模型共有三种形式即先串后并法、先并后串法,以及一体化耦合模型。借助预报误差序列的自相关性以及水文模型之间的互补性,减少了预报误差,提高了预报精度,从而实现了高精度的水文实时预报,为科学防控洪水提供了理论依据和技术支持。如图1所示,为本发明串并联耦合的多模型水文预报方法具体流程,所述方法包括以下步骤:
[0045] 步骤1:选择预报断面,获取该断面历史时期实测降雨和流量资料,将资料分为率定期和检验期两部分;采用率定期的数据,率定模型的参数,采用检验期的数据,检验模型的效果。建立多个水文预报模型,如新安江模型、水箱模型以及API模型等,分别预报该断面的径流量。
[0046] 步骤2:建立串并联耦合的水文预报模型。
[0047] 本发明提出了一种串并联耦合的水文预报方法,按耦合方式不同,共分为三种形式,即先串后并法、先并后串法,以及一体化耦合模型。
[0048] (1)先串后并法
[0049] ①基于实测流量和预报流量,得到预报误差序列。依据预报误差序列,建立误差自回归串联校正模型。
[0050] 具体建模方法如下:
[0051] 令 实 测 序 列 为 Qt=(Q 1,Q2,...,Qn),n为 序 列的 长 度,预 测 序 列 为预测序列与实测序列的差值为误差序列,记为et=(e1,e2,...,en)。误差自回归模型的数学表达式为:
[0052] et=θ1et-1+θ2et-2+…+θqet-q+ξ
[0053] 式中,q为误差自回归模型的阶数,θ1,θ2,...,θq为回归系数,ξ是均值为零、方差为某值的白噪声信号
[0054] 采用上式,建立q阶自回归模型AR(q),q=2,3,…,30。基于率定期的数据,利用最小二乘法率定模型的参数;基于检验期的数据,检验模型的效果。
[0055] 通过AIC准则来选取模型的阶数。q阶自回归模型AR(q)的AIC(q)值为:
[0056]
[0057] 选用检验期AIC值最小的自回归模型作为最优模型,参与校正。
[0058] 最后,借助所建自回归模型预测t时刻的预报误差et,并对预报结果进行校正。针对每个水文预报模型,利用上述方法,得到各个预报模型在该断面预报流量的串联校正结果。
[0059] ②根据断面实测流量以及①中各水文模型的串联校正结果,建立并联校正模型。设参与并联校正的水文模型个数为m,各水文模型校正后的预报结果分别为可表示为:
[0060]
[0061] 并联校正的结果通过下式来计算:
[0062]
[0063] 其中,ω1,ω2,...,ωm分别为各个水文模型的耦合权重,且满足条件ω1+ω2+…+ωm=1。
[0064] 采用最小二乘法来确定并联校正的权重,具体如下。
[0065] 计算各个模型的误差序列,得到各个模型的误差序列如下所示:
[0066]
[0067] 并联校正结果的误差序列记为(Ft-Qt),用符号E来表示期望值,误差序列平方的2
期望值记为E(Ft-Qt),其计算式如下:
[0068]
[0069] 为了使并联校正结果误差平方的期望值最小,求解该目标下各个模型权重值的问题转化为求解如下线性规划问题:
[0070]
[0071] 引入拉格朗日乘子λ,构建目标函数:
[0072] L(ω1,ω2,...,ωm,λ)=E(ω1e1t+ω2e2t+…+ωmemt)2+λ(ω1+ω2+…+ωm-1)[0073] 目标函数分别对ω1,ω2,...,ωm,λ求偏导,并令其偏导等于0,得到如下方程组:
[0074]
[0075] 其中,
[0076] 通过求解方程组,即可得到各个水文模型的耦合权重ω1,ω2,...,ωm的值。
[0077] (2)先并后串方法
[0078] ①基于m个水文模型的预报序列Q1t,Q2t,...,Qmt,建立并联校正模型,如下所示:
[0079] Ft=ω1Q1t+ω2Q2t+…+ωmQmt
[0080] 式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足ω1+ω2+…+ωm=1。利用上述最小二乘法原理求解耦合权重ω1,ω2,...,ωm,由此得到了多个预报模型的并联校正结果。
[0081] ②根据断面实测径流量和多个预报模型在该断面预报径流量的并联预报结果,建立误差自回归实时串联校正模型。利用上述自回归模型原理,对多个预报模型的并联校正结果进行串联校正。由此得到了先并后串法的计算结果。
[0082] (3)一体化耦合方法
[0083] ①利用断面实测径流量和各个水文预报模型模拟预报的径流量,建立一体化耦合模型。记水文预报模型个数为m,实测序列为Qt,各个水文模型的预报结果为Q1t,Q2t,...,Qmt。对各个模型采用误差自回归模型进行校正,记各个模型的误差校正表达式如下:
[0084]
[0085] 在获取各个模型的串联校正结果的基础上,利用最小二乘法,进行多个模型的并联校正。并联校正的结果如下所示:
[0086]
[0087] 式中,ω1,ω2,...,ωm分别为各个水文预报模型的耦合权重,且满足ω1+ω2+…+ωm=1。
[0088] 组合预报结果的误差记为Ft-Qt,其期望的平方记为E(Ft-Qt)2,为了使E(Ft-Qt)2最小,得到如下目标函数:
[0089]
[0090] 式中,其参数包含回归系数(α11,α21),(α12,α22),…,(α1m,α2m)以及权重系数ω1,ω2,…,ωm两部分,共3m个待定参数。
[0091] 引入粒子群智能优化算法(Particle Swarm Optimizer,PSO)同时求解这些参数。得到这些参数后,先采用串联校正方法,对各水文模型的预报流量进行串联校正,得到串联校正结果;而后,利用并联校正方法,对各水文模型串联校正后的结果进行并联校正。由此得到了一体化耦合方法的计算结果。
[0092] 步骤3:选取平均绝对误差、平均相对误差、确定性系数和均方根误差等常用的预报精度评定指标,对原模型以及串并联耦合校正的结果进行精度评定。各个精度评定指标的定义如下:
[0093] ①平均绝对误差(MAE)
[0094] 平均绝对误差如下式所示:
[0095]
[0096] 式中,Qi为实测值, 为预报值,n为预报序列长度。
[0097] ②平均相对误差(MRE)
[0098] 平均相对误差如下式所示:
[0099]
[0100] ③确定性系数(DC)
[0101] 确定性系数如下式所示:
[0102]
[0103] 式中, 为实测值在预报序列长度内的平均值,DC反映了洪水预报过程与实测过程之间的吻合程度,DC值越高,吻合程度越高。
[0104] ④均方根误差(RMSE)
[0105] 均方根误差如下式所示:
[0106]
[0107] 采用新安江、水箱、神经网络以及灰色模型预报三峡水库入库流量序列,其精度评定结果如表1所示;利用本发明所提串并联耦合模型对预报结果进行校正,其精度评定结果如表2所示。
[0108] 表1单一水文模型预报精度评定结果
[0109]
[0110] 表2耦合水文预报模型精度评定结果
[0111]
[0112] 本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈