[发明专利]一种考虑水库调洪演算的多变量洪水重现期计算方法

权利要求书

1.一种考虑水库调洪演算的多变量洪水重现期计算方法，其特征在于包括如下步骤：

步骤1，收集水库历年最大洪水过程资料和洪峰、洪量系列；

步骤2，将步骤1中的历年最大洪水过程输入水库按照调度规则进行调洪演算，从输出的库水位过程线统计得到年最高水位系列；

步骤3，根据步骤1中的洪峰、洪量系列和步骤2中得到的年最高水位系列，利用遗传规划进行符号回归，推求水库调洪函数的显式表达式，并确定指定年最高水位值对应的安全区域；

步骤4，根据步骤1中的洪峰、洪量系列，选取边缘概率分布函数线型，并估计边缘概率分布函数的参数；

步骤5，采用Copula函数构造洪峰、洪量的联合概率分布函数，并估计Copula函数的参数，求解联合概率密度函数；

步骤6，对步骤5得到的联合概率密度函数在步骤3中确定的安全区域内进行积分运算，获取指定年最高水位的累积概率；

步骤7，根据步骤6中得到的累积概率计算指定年最高水位的单变量重现期，作为考虑水库调洪演算的多变量洪水重现期。

2.如权利要求1所述的方法，其特征在于：所述步骤2中，水库调洪演算采用四阶龙格-库塔算法进行数值求解。

3.如权利要求1所述的方法，其特征在于：所述步骤4中，将P-III型分布作为洪峰和洪量的边缘概率分布函数线型，采用线性矩法估计边缘概率分布函数的参数。

4.如权利要求1所述的方法，其特征在于：所述步骤5中，采用Gumbel-Hougaard Copula函数构造洪峰和洪量的联合概率分布函数，采用Kendall秩相关性系数法估计Gumbel-Hougaard Copula函数的参数。

5.如权利要求1所述的方法，其特征在于：所述步骤6中，利用蒙特卡洛模拟方法计算定积分的值。