长江科学院院报 ›› 2021, Vol. 38 ›› Issue (11): 18-24.DOI: 10.11988/ckyyb.20200796

• 水资源 • 上一篇    下一篇

基于Copula函数的洪水峰量联合设计方法研究

苗正伟1, 丁志宏2, 路梅1, 杨学智1   

  1. 1.河北水利电力学院 水利工程学院,河北 沧州 061001;
    2.中水北方勘测设计研究有限责任公司,天津 300222
  • 收稿日期:2020-08-07 修回日期:2020-12-03 出版日期:2021-11-01 发布日期:2021-11-08
  • 作者简介:苗正伟(1981 -),男,山东聊城人,副教授,硕士,主要从事水文水资源方面的研究。E-mail:newstart2017@sina.com
  • 基金资助:
    河北省自然科学基金项目(E2020412219);河北省水利科技计划项目(2017-62);河北省教育厅青年基金项目(QN2019088)

Joint Design of Flood Peak and Volume Based on Copula Function

MIAO Zheng-wei1, DING Zhi-hong2, LU Mei1, YANG Xue-zhi1   

  1. 1. Department of Hydraulic Engineering, Hebei University of Water Resources and Electric Engineering, Cangzhou 061001, China;
    2. China Water Resources Beifang Investigation, Design and Research Co., Ltd., Tianjin 300222, China
  • Received:2020-08-07 Revised:2020-12-03 Online:2021-11-01 Published:2021-11-08

摘要: 为明晰洪水峰量联合设计的特点,以岗南水库洪水为例,基于Gumbel Copula函数,分析了AND、OR、Kendall、生存Kendall 4种重现期的优缺点,采用极大似然法、同频率法、条件最可能组合法3种方法计算了联合设计值。结果表明: ①AND和OR重现期在危险域和安全域的识别上存在局限性;相对而言,Kendall重现期更合理,但其安全域是无界的,这与实际不符;生存Kendall重现期则界定了有界的安全域,使得重现期的概念在逻辑上更科学合理。② 3种设计值计算方法的差别不大,但从简单实用角度出发,推荐采用同频率法计算设计值。③不同重现期标准的设计值差别比较明显,基于OR重现期计算的设计值总是最大的,生存Kendall、Kendall重现期设计值次之,AND重现期设计值最小。④推荐采用生存Kendall重现期进行两变量洪水设计,因其有比较严谨的理论基础,且设计结果兼顾了安全性与经济性。⑤两变量联合设计值与单变量设计值的差异受变量间相关性的影响较大,且变量相关性越弱,差异越大。研究显示,基于生存Kendall重现期、采用同频率法计算设计值是目前较为科学合理的洪水峰量联合设计途径。

关键词: 洪水, 洪量, Kendall函数, Copula函数, 重现期

Abstract: With the flood peak and flood volume of Gangnan reservoir as a case study, the “OR” return periods, “AND” return periods, Kendall return periods (hereinafter referred to as KRP) and survival Kendall return periods (SKRP) of bivariate joint distribution were computed by using the optimally fitted Gumbel copula. The joint design quantiles in a bivariate environment was calculated based on the maximum likelihood method, conditional most likely combination, and bivariate equivalent frequency combination. The main conclusions of this study are summarized as follows: (1) The traditional “AND” and “OR” approach in multivariate return periods definition is limited in identifying safe and dangerous regions, while the KRP based approach would be more rational, but the safe region of KRP is unbounded, which is inconsistent with reality. As specific features, the proposed approach of SKRP yields a bounded safe region (a natural request in applications), where all the variables of interest are finite and limited. Consequently, the SKRP is logically the most reasonable. (2) The design values calculated using these three methods see no considerable difference, but from the point of view of simplicity and practicality, the bivariate equivalent frequency combination is recommended. (3) The design values based on the return period definitions differ obviously, among which OR return period yields the maximum design value, followed by SKRP, KRP, and AND in sequence. (4) SKRP is recommended for bivariate flood design because of its rigorous theoretical basis, and its design results give consideration to both safety and economy. (5) The difference between bivariate joint design value and univariate design value is greatly affected by the correlation between variables. The weaker the correlation between variables, the greater the difference. In summary, KRP based bivariate equivalent frequency combination is currently the more scientific and rational approach to the joint design of flood peak and volume.

Key words: flood, flood volume, Kendall function, Copula, return period

中图分类号: