(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111572156.3 (22)申请日 2021.12.21 (71)申请人 武汉大学 地址 430072 湖北省武汉市武昌区珞珈山 武汉大学 申请人 中国长江三峡集团有限公司 (72)发明人 巩钰　刘攀　曹光荣　张玮　 刘旖旎　夏倩　 (74)专利代理 机构 武汉科皓知识产权代理事务 所(特殊普通 合伙) 42222 代理人 肖明洲 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 50/06(2012.01) G06Q 50/26(2012.01) (54)发明名称 一种水库两阶段实时防洪风险计算方法及 系统 (57)摘要 本发明公开了一种水库两阶段实时防洪风 险计算方法及系统， 根据预报信息， 推导预见期 内水库多个时段库容的联合分布， 对 预见期内防 洪风险进行计算； 通过历史资料进行常规调度， 拟合预见期外风险与预见期末库容的函数关系， 对预见期外防洪风险进行评估； 通过连接函数 (Copula函数)联合预见期内外风险， 计算总防洪 风险。 本发明采用解析方法， 对水库防洪风险进 行了准确评估， 避免了传统的模拟风险评估 方法 对样本容量的依赖， 可广泛应用于水库防洪风险 的评估， 为科 学决策提供依据。 权利要求书2页 说明书5页 附图3页 CN 114331072 A 2022.04.12 CN 114331072 A 1.一种水库两阶段实时防洪风险计算方法， 其特 征在于， 包括以下步骤： 步骤1： 根据预报信息， 推导预见期内水库多个时段库容的联合分布， 对预见期内防洪 风险进行计算； 当入库径流预报服从正态分布且可用一阶自回归模型描述 时， 若第一时段入流与自回 归模型各时段残差独立， 推求预 见期内库容 服从多元正态分布； 预 见期内风险计算式如下： 其中， RS为预见期内风险， Vc为临界库容， g()为预见期内库容联合概率密度函数， Vi为 时段i末库容， H为预 见期长度； 步骤2： 通过历史资料进行常规调度， 拟合预见期外风险与预见期末库容的函数关系， 对预见期外防洪风险进行评估； 通过常规防洪调度得预见期末库容VH和预见期外风险RL′的关系， 预见期外平均风险计 算式如下： 其中， c和s是方程的参数， 通过适线法求得； σVH是预见期末库容的方差， 表示预见期 末库容的均值， P(VH)表示预见期末库容 为VH的概率； 步骤3： 通过 连接函数 联合预见期内外风险， 计算总防洪风险； 其中， θ 是Gumbel ‑Hougaard Copula函数的参数。 2.一种水库两阶段实时防洪风险计算系统， 其特 征在于， 包括以下模块： 模块1， 用于根据预报信息， 推导预见期内水库多个时段库容的联合分布， 对预见期内 防洪风险进行计算； 当入库径流预报服从正态分布且可用一阶自回归模型描述 时， 若第一时段入流与自回 归模型各时段残差独立， 推求预 见期内库容 服从多元正态分布； 预 见期内风险计算式如下： 其中， RS为预见期内风险， Vc为临界库容， g()为预见期内库容联合概率密度函数， Vi为 时段i末库容， H为预 见期长度； 模块2， 用于通过历史资料进行常规调度， 拟合预见期外风险与预见期末库容的函数关 系， 对预见期外防洪风险进行评估；权　利　要　求　书 1/2 页 2 CN 114331072 A 2通过常规防洪调度得预见期末库容VH和预见期外风险RL′的关系， 预见期外平均风险计 算式如下： 其中， c和s是方程的参数， 通过适线法求得； σVH是预见期末库容的方差， 表示预见期 末库容的均值， P(VH)表示预见期末库容 为VH的概率； 模块3， 用于通过 连接函数 联合预见期内外风险， 计算总防洪风险； 其中， θ 是Gumbel ‑Hougaard Copula函数的参数。权　利　要　求　书 2/2 页 3 CN 114331072 A 3