一种自适应气候变化的水库多变量设计洪水推求方法技术

本发明专利技术公开了一种自适应气候变化的水库多变量设计洪水推求方法，本发明专利技术建立GCM与VIC的耦合模型，以预测未来气候变化情景下的水库径流过程，通过Copula函数构建了水库设计洪水峰量最可能组合法计算模型，并将未来气候变化情景下的水库预测径流资料作为设计洪水峰量最可能组合法计算模型的输入，计算水库在不同重现期水平下的联合设计值。本发明专利技术具有较强的统计基础，能够客观地反映实测样本的特征，并充分考虑各个特征量之间的内在相关性；考虑气候变化条件下的设计洪水设计值，更具实用性。

A design method of multi variable design flood based on adaptive climate change

The design flood calculation method of multi variable reservoir, the invention discloses an adaptive climate change, the invention establishes a coupling model of GCM and VIC, to predict the future climate change scenarios of reservoir runoff process, through the Copula function constructed reservoir design flood peak most likely group legal model, and the future climate change scenarios in reservoir the prediction of runoff data as design flood peak most likely model combination method of input, calculation of reservoir in different return period level of joint design. The invention has strong statistical basis, can objectively reflect the characteristics of the measured sample, and considering the correlation between the amount of each feature; design considerations for design flood under climate change conditions, more practical.

【技术实现步骤摘要】

本专利技术涉及水库防洪安全设计
，具体地说涉及一种自适应气候变化的水库多变量设计洪水推求方法。
技术介绍
设计洪水是防洪规划的基本依据，现阶段我国工程实践中主要采用单变量同频率法来推求设计洪水，该法假设洪峰和不同时段的洪量完全相关，即假定洪峰流量和各时段洪量都符合于同一设计频率，然后通过各水文变量的边缘分布来推求设计洪水设计值。这一假设未能充分考虑各个特征量之间的内在相关性。近年来，随着Copula函数在水文频率分析领域的深入发展，国内外学者提出了采用Copula函数构造洪峰和洪量的联合分布，从而推求设计洪水设计值的方法。例如：肖义、郭生练、刘攀等发表的基于Copula函数的设计洪水过程线方法[J].武汉大学学报(工学版)，2007,40(4):13-17；基于Copula函数构造了洪峰与7日洪量的两变量联合分布，并基于联合分布推求了隔河岩水库在两变量重现期下的设计洪水过程线。李天元、郭生练、闫宝伟等发表的基于多变量联合分布推求设计洪水过程线的新方法[J].水力发电学报,2013,32(3):10-14,38；应用Copula函数构造了洪峰与时段洪量之间的多变量联合分布，推求了三变量重现期下的设计洪水过程线，并将该方法应用在三峡水库的设计洪水研究中。然而，他们推求的峰量组合是联合重现期下等值线，仍然采用了同频率假设来确定洪峰和洪量的设计值。为了研究洪水峰量的内在相关性，公开号为CN104615907A的专利技术专利申请提供了一种基于多变量条件最可能组合推求设计洪水过程线的方法，该法以条件概率密度函数最大为原则，求解得到洪峰、洪量的条件最可能组合的估计值，以单变量条件重现期作为设计洪水的重现期，从而推求不同设计频率下的设计洪水过程线。在工程设计中，人们通过对实际发生洪水的时空特性规律分析，通常关心最可能发生的洪水模式，现阶段对水库最可能洪水设计值的研究较少。由于人类活动影响和下垫面变化，气候变化对水库的设计洪水计算产生较大的影响，而现有方法只注重当前入库径流规律，没有文献研究设计洪水对未来气候变化的自适应问题。
技术实现思路
本专利技术的目的是克服现有技术存在的不足，基于全球气候模式(GCM)与可变下渗能力水文模型(VIC)获取未来气候变化情景下的预测径流信息，并考虑峰量最可能组合方法，提供一种自适应气候变化的水库多变量设计洪水推求方法。本专利技术一种自适应气候变化的水库多变量设计洪水推求方法，包括以下步骤：步骤1，建立全球气候模式与可变下渗能力水文模型的耦合模型，以预测未来气候变化情景下的水库径流过程；步骤2，基于Copula函数建立水库设计洪水峰量最可能组合法计算模型；步骤3，将未来气候变化情景下的水库预测径流资料作为设计洪水峰量最可能组合法计算模型的输入，计算水库的设计洪水值。步骤2进一步包括以下子步骤：(2.1)构造边缘分布和联合分布函数；令Q表示洪峰流量；Wi(i＝1,2,...,n)表示不同时段的洪量，n表示水库推求设计洪水时需要的时段洪量数；采用P-III分布函数构建洪峰、洪量的边缘分布，分别记为FQ(q)、；密度函数分别记为fQ(q)、；采用线性矩法估计P-III分布函数的参数，通过Copula函数构造峰量的联合分布；(2.2)基于或重现期(OR重现期)推求水库设计洪水最可能值；采用或重现期(OR重现期)作为水库防洪标准的度量指标；通过Q和Wi(i＝1,2,…,n)的联合概率密度函数值f(x,y1,y2,…,yn-1,yn)来度量不同峰量组合发生的相对可能性大小；水库设计洪水峰量最可能组合法计算模型即是求解水库的洪峰和不同时段洪量在满足防洪标准条件下，联合概率密度函数值最大的设计值组合。所述步骤(2.1)中，以ArchimedeanCopula函数族中的Gumbel-HougaardCopula函数作为联合分布函数，构造洪峰与洪量的联合分布。所述步骤(2.2)中，采用拉格朗日乘数法求解最可能组合问题，给定联合重现期T，构造拉格朗日函数，分别对q,w1,…,wn和λ求偏导，并令其为0，可求得拉格朗日函数的所有极值点；采用多变量同频率组合计算得到的峰量设计结果作为初始解，采用调和平均数牛顿法进行迭代求解，得到水库设计洪水峰量最可能组合。所述步骤3具体为：将水库未来气候变化情景下的预测径流资料进行整理分析，通过年最大值取样方法，得到预测系列的年最大洪峰和洪量值，然后构造年最大洪峰和洪量的边缘分布函数FQ(q)、，代入步骤2中，从而可以求得考虑未来气候的水库多变量设计洪水设计值。与现有技术相比，本专利技术的有益效果在于：1、科学合理、贴近工程实际；本专利技术利用Copula函数建立洪峰与不同时段洪量的联合分布，基于联合概率密度最大原则，推导得到最可能组合法的洪峰和洪量设计值，具有较强的统计基础，能够客观地反映实测样本的特征，并充分考虑各个特征量之间的内在相关性。2、可为气候变化条件下水库推求设计洪水提供重要且可操作性强的参考依据：本专利技术建立GCM与VIC的耦合模型，预测未来气候变化情景下的水库径流过程，并将预测径流资料输入设计洪水计算模型，以适应未来气候变化，更具实用性。附图说明图1为本专利技术方法的具体流程图；图2为GCM与VIC模型的耦合原理图；图3为联合重现期水平下峰量等值曲面的示意图。具体实施方式下面结合附图对本专利技术进行进一步说明。本专利技术基于全球气候模式(GCM)与可变下渗能力水文模型(VIC)获取未来气候变化情景下的预测径流信息，并考虑峰量最可能组合方法，提供了一种自适应气候变化的水库多变量设计洪水推求方法，其具体流程详见图1。下面通过实施例，并结合附图，对本专利技术的技术方案做进一步具体说明：步骤1，建立GCM与VIC的耦合模型，以预测未来气候变化情景下的水库径流过程。相对于集总式水文模型，分布式水文模型可以将土壤、植被的空间分布、土地利用的区域特征以及降水、蒸发和径流的空间分布较好的反映出来，并将其建立在网格点上，更适合与GCM耦合。本实施例中，应用统计降尺度方法，应用于全球气候模式GCM，首先实现对未来降水变化的预测分析，然后通过降尺度方法处理GCM的输出(降水和气温)作为流域水文模型VIC的输入，模拟预测未来气候情景下汉江流域的径流过程。基于NCEP观测资料和统计降尺度方法的VIC模型与GCM耦合机本文档来自技高网...

【技术保护点】
一种自适应气候变化的水库多变量设计洪水推求方法，其特征在于包括以下步骤：步骤(1)，建立全球气候模式与可变下渗能力水文模型的耦合模型，以预测未来气候变化情景下的水库径流过程；步骤(2)，基于Copula函数建立水库设计洪水峰量最可能组合法计算模型；步骤(3)，将未来气候变化情景下的水库预测径流资料作为设计洪水峰量最可能组合法计算模型的输入，计算水库的设计洪水值。

【技术特征摘要】
1.一种自适应气候变化的水库多变量设计洪水推求方法，其特征在于
包括以下步骤：
步骤(1)，建立全球气候模式与可变下渗能力水文模型的耦合模型，
以预测未来气候变化情景下的水库径流过程；
步骤(2)，基于Copula函数建立水库设计洪水峰量最可能组合法计算
模型；
步骤(3)，将未来气候变化情景下的水库预测径流资料作为设计洪水
峰量最可能组合法计算模型的输入，计算水库的设计洪水值。
2.如权利要求1所述的方法，其特征在于步骤(2)进一步包括以下
子步骤：
(2.1)构造边缘分布和联合分布函数；令Q表示洪峰流量；Wi(i＝1,2,...,n)
表示不同时段的洪量，n表示水库推求设计洪水时需要的时段洪量数；采用
P-III分布函数构建洪峰、洪量的边缘分布，分别记为FQ(q)、密度
函数分别记为fQ(q)、采用线性矩法估计P-III分布函数的参数，通
过Copula函数构造峰量的联合分布；
(2.2)基于或重现期推求水库设计洪水最可能值；采用或重现期作
为水库防洪标准的度量指标；通过Q和Wi(i＝1,2,…,n)的联合概率密
度函数值f(x,y1,y2,…,yn-1,yn)来度量不同峰量组合发生的相...

【专利技术属性】
技术研发人员：郭生练，尹家波，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42