• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Copula函数在水文水资源中的研究进展与述评

刘章君 郭生练 许新发 许世超 成静清

刘章君, 郭生练, 许新发, 许世超, 成静清. Copula函数在水文水资源中的研究进展与述评[J]. 水科学进展, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
引用本文: 刘章君, 郭生练, 许新发, 许世超, 成静清. Copula函数在水文水资源中的研究进展与述评[J]. 水科学进展, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
LIU Zhangjun, GUO Shenglian, XU Xinfa, XU Shichao, CHENG Jingqing. Application of Copula functions in hydrology and water resources: a state-of-the-art review[J]. Advances in Water Science, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
Citation: LIU Zhangjun, GUO Shenglian, XU Xinfa, XU Shichao, CHENG Jingqing. Application of Copula functions in hydrology and water resources: a state-of-the-art review[J]. Advances in Water Science, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015

Copula函数在水文水资源中的研究进展与述评

doi: 10.14042/j.cnki.32.1309.2021.01.015
基金项目: 

国家自然科学基金资助项目 51909112

国家自然科学基金资助项目 51879192

详细信息
    作者简介:

    刘章君(1991-), 男, 江西吉安人, 博士, 主要从事水文分析与计算研究。E-mail:liuzhangjun@whu.edu.cn

    通讯作者:

    郭生练, E-mail:slguo@whu.edu.cn

  • 中图分类号: P333;TV12

Application of Copula functions in hydrology and water resources: a state-of-the-art review

Funds: 

the National Natural Science Foundation of China 51909112

the National Natural Science Foundation of China 51879192

  • 摘要: Copula函数是一种灵活构造多变量联合分布和处理多变量问题的有效工具,在水文水资源领域应用广泛;随着研究的广度和深度不断推进,显示出其良好的适用性和不可替代的优越性。重点综述了近10年来Copula函数在水文事件多变量频率分析、水文事件遭遇组合分析、水文随机模拟、水文模型与预报以及其他问题中的最新应用研究进展与不足。从应用条件、参数估计、不对称性、拟合优选、尾部特性、多变量重现期选取、成果抽样误差等方面阐述了Copula函数应用需要注意的问题。展望了未来Copula函数在水文水资源领域应用中的研究重点和发展方向,为今后的研究提供借鉴和参考。
  • 表  1  Copula函数在水文水资源应用领域的主要研究问题及相应变量

    Table  1.   Main research problems and relevant variables of Copula function applied in the field of hydrology and water resources

    研究方面 研究问题 研究变量
    水文事件多变量频率分析 降雨事件 降雨强度、降雨总量和降雨历时等
    洪水事件 洪峰、洪量和洪水历时等;不同分期内的洪水;发生时间和量级
    干旱事件 干旱烈度、干旱强度、干旱历时、干旱间隔时间等
    水文事件遭遇组合分析 洪水遭遇 干支流不同站点的洪水
    设计洪水地区组成 水库断面和区间洪水
    复合洪涝事件 暴雨和外江水位;暴雨、河洪和潮位
    丰枯遭遇 不同站点/地区的降水或径流
    水文随机模拟 降雨随机模拟 降雨总量、降雨历时、降雨强度、降雨间隔时间等
    洪水随机模拟 洪峰、洪量和洪水过程线形状
    径流随机模拟 相邻时刻径流;不同站点径流
    水文模型与预报 水文预报不确定性评估 预报流量和实际流量
    水文模型输入、参数、输出相关性模拟 水文模型输入、参数、输出
    其他问题 水量水质组合 水量和水质指标
    降雨—径流关系变异诊断 降雨和径流
    水文气象站网优化设计 不同站点日降雨量
    水库随机优化调度 相邻时段的水库入库径流
    河流水沙关系 径流与泥沙
    下载: 导出CSV
  • [1] 闫宝伟, 潘增, 薛野, 等.论水文计算中的相关性分析方法[J].水利学报, 2017, 48(9):1039-1046. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201709005.htm

    YAN B W, PAN Z, XUE Y, et al. Modeling dependence and correlation in hydrological calculation[J]. Journal of Hydraulic Engineering, 2017, 48(9):1039-1046. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201709005.htm
    [2] NELSEN R B. An introduction to copulas[M]. 2nd ed. Berlin:Springer, 2006.
    [3] 郭生练, 闫宝伟, 肖义, 等. Copula函数在多变量水文分析计算中的应用及研究进展[J].水文, 2008, 28(3):1-7. https://www.cnki.com.cn/Article/CJFDTOTAL-SWZZ200803002.htm

    GUO S L, YAN B W, XIAO Y, et al. Application of Copula function in multivariate hydrological analysis and estimation[J]. Journal of China Hydrology, 2008, 28(3):1-7. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWZZ200803002.htm
    [4] de MICHELE C, SALVADORI G. A Generalized Pareto intensity-duration model of storm rainfall exploiting 2-Copulas[J]. Journal of Geophysical Research:Atmospheres, 2003, 108(D2):4067. doi:  10.1029/2002JD002534
    [5] 熊立华, 郭生练.两变量极值分布在洪水频率分析中的应用研究[J].长江科学院院报, 2004, 21(2):35-37. doi:  10.3969/j.issn.1001-5485.2004.02.011

    XIONG L H, GUO S L. Application study of a bivariate extremal distribution in flood frequency analysis[J]. Journal of Yangtze River Scientific Research Institute, 2004, 21(2):35-37. (in Chinese) doi:  10.3969/j.issn.1001-5485.2004.02.011
    [6] SCHWEIZER B. Introduction to copulas[J]. Journal of Hydrologic Engineering, 2007, 12(4):346. doi:  10.1061/(ASCE)1084-0699(2007)12:4(346)
    [7] ZHANG L, SINGH V P. Bivariate rainfall frequency distributions using Archimedean copulas[J]. Journal of Hydrology, 2007, 332(1/2):93-109.
    [8] 许月萍, 童杨斌, 富强, 等.几种Copulas模拟不同历时降雨量的影响分析[J].浙江大学学报(工学版), 2009, 43(6):1107-1111. doi:  10.3785/j.issn.1008-973X.2009.06.024

    XU Y P, TONG Y B, FU Q, et al. Impact analysis for rainfall depthsimulation of different durations through several Copulas[J]. Journal of Zhejiang University(Engineering Science), 2009, 43(6):1107-1111. (in Chinese) doi:  10.3785/j.issn.1008-973X.2009.06.024
    [9] VANDENBERGHE S, VERHOEST N E C, de BAETS B. Fitting bivariate copulas to the dependence structure between storm characteristics:a detailed analysis based on 105 year 10 min rainfall[J]. Water Resources Research, 2010, 46(1):489-496.
    [10] BALISTROCCHI M, BACCHI B. Modelling the statistical dependence of rainfall event variables through copula functions[J]. Hydrology and Earth System Sciences, 2011, 15(6):1959-1977. doi:  10.5194/hess-15-1959-2011
    [11] JUN C, QIN X S, GAN T Y, et al. Bivariate frequency analysis of rainfall intensity and duration for urban stormwater infrastructure design[J]. Journal of Hydrology, 2017, 553:374-383. doi:  10.1016/j.jhydrol.2017.08.004
    [12] LI H S, WANG D, SINGH V P, et al. Non-stationary frequency analysis of annual extreme rainfall volume and intensity using Archimedean Copulas:a case study in Eastern China[J]. Journal of Hydrology, 2019, 571:114-131. doi:  10.1016/j.jhydrol.2019.01.054
    [13] ZHANG L, SINGH V P. Bivariate flood frequency analysis using the copula method[J]. Journal of Hydrologic Engineering, 2006, 11(2):150-164. doi:  10.1061/(ASCE)1084-0699(2006)11:2(150)
    [14] SALVADORI G, de MICHELE C, DURANTE F. On the return period and design in a multivariate framework[J]. Hydrology and Earth System Sciences, 2011, 15(11):3293-3305. doi:  10.5194/hess-15-3293-2011
    [15] GRÄLER B, van den BERG M J, VANDENBERGHE S, et al. Multivariate return periods in hydrology:a critical and practical review focusing on synthetic design hydrograph estimation[J]. Hydrology and Earth System Sciences, 2013, 17(4):1281-1296. doi:  10.5194/hess-17-1281-2013
    [16] VOLPI E, FIORI A. Hydraulic structures subject to bivariate hydrological loads:return period, design, and risk assessment[J]. Water Resources Research, 2014, 50(2):885-897. doi:  10.1002/2013WR014214
    [17] 李天元, 郭生练, 刘章君, 等.基于峰量联合分布推求设计洪水[J].水利学报, 2014, 45(3):269-276. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201403003.htm

    LI T Y, GUO S L, LIU Z J, et al. Design flood estimation based on bivariate joint distribution of flood peak and volume[J]. Journal of Hydraulic Engineering, 2014, 45(3):269-276. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201403003.htm
    [18] 陈子燊, 黄强, 刘曾美.基于非对称Archimedean Copula的三变量洪水风险评估[J].水科学进展, 2016, 27(5):763-771. doi:  10.14042/j.cnki.32.1309.2016.05.014

    CHEN Z S, HUANG Q, LIU Z M. Risk assessment of trivariate flood based on asymmetric Archimedean Copulas[J]. Advances in Water Science, 2016, 27(5):763-771. (in Chinese) doi:  10.14042/j.cnki.32.1309.2016.05.014
    [19] BALISTROCCHI M, ORLANDINI S, RANZI R, et al. Copula-based modeling of flood control reservoirs[J]. Water Resources Research, 2017, 53(11):9883-9900. doi:  10.1002/2017WR021345
    [20] 刘章君, 郭生练, 许新发, 等.两变量洪水结构荷载重现期与联合设计值研究[J].水利学报, 2018, 49(8):956-965. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201808007.htm

    LIU Z J, GUO S L, XU X F, et al. Bivariate structure load return period and joint flood quantile estimation[J]. Journal of Hydraulic Engineering, 2018, 49(8):956-965. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201808007.htm
    [21] XIONG L H, JIANG C, XU C Y, et al. A framework of change-point detection for multivariate hydrological series[J]. Water Resources Research, 2015, 51(10):8198-8217. doi:  10.1002/2015WR017677
    [22] YIN J B, GUO S L, HE S K, et al. A copula-based analysis of projected climate changes to bivariate flood quantiles[J]. Journal of Hydrology, 2018, 566:23-42. doi:  10.1016/j.jhydrol.2018.08.053
    [23] QI W, LIU J G. A non-stationary cost-benefit based bivariate extreme flood estimation approach[J]. Journal of Hydrology, 2018, 557:589-599. doi:  10.1016/j.jhydrol.2017.12.045
    [24] JIANG C, XIONG L H, YAN L, et al. Multivariate hydrologic design methods under nonstationary conditions and application to engineering practice[J]. Hydrology and Earth System Sciences, 2019, 23(3):1683-1704. doi:  10.5194/hess-23-1683-2019
    [25] MICHAILIDI E M, BACCHI B. Dealing with uncertainty in the probability of overtopping of a flood mitigation dam[J]. Hydrology and Earth System Sciences, 2017, 21(5):2497-2507. doi:  10.5194/hess-21-2497-2017
    [26] 尹家波, 郭生练, 吴旭树, 等.两变量设计洪水估计的不确定性及其对水库防洪安全的影响[J].水利学报, 2018, 49(6):715-724. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201806008.htm

    YIN J B, GUO S L, WU X S, et al. Uncertainty of bivariate design flood estimation and its impact on reservoir flood prevention[J]. Journal of Hydraulic Engineering, 2018, 49(6):715-724. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201806008.htm
    [27] LIU Y R, LI Y P, MA Y, et al. Development of a Bayesian-copula-based frequency analysis method for hydrological risk assessment:the Naryn River in Central Asia[J]. Journal of Hydrology, 2020, 580:124349. doi:  10.1016/j.jhydrol.2019.124349
    [28] 肖义, 郭生练, 刘攀, 等.分期设计洪水频率与防洪标准关系研究[J].水科学进展, 2008, 19(1):54-60. doi:  10.3321/j.issn:1001-6791.2008.01.009

    XIAO Y, GUO S L, LIU P, et al. Seasonal flood frequency analysis and flood prevention standard[J]. Advances in Water Science, 2008, 19(1):54-60. (in Chinese) doi:  10.3321/j.issn:1001-6791.2008.01.009
    [29] CHEN L, GUO S L, YAN B W, et al. A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence[J]. Hydrological Sciences Journal, 2010, 55(8):1264-1280. doi:  10.1080/02626667.2010.520564
    [30] 覃光华, 曹泠然, 王文圣, 等.以前期水量为条件的管运洪水推求[J].工程科学与技术, 2019, 51(5):17-24. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201905003.htm

    QIN G H, CAO L R, WANG W S, et al. Calculation of manage-flood based on the previous flood volume[J]. Advanced Engineering Sciences, 2019, 51(5):17-24. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201905003.htm
    [31] 陆桂华, 闫桂霞, 吴志勇, 等.基于Copula函数的区域干旱分析方法[J].水科学进展, 2010, 21(2):188-193. http://skxjz.nhri.cn/article/id/279

    LU G H, YAN G X, WU Z Y, et al. Regional drought analysis approach based on Copula function[J]. Advances in Water Science, 2010, 21(2):188-193. (in Chinese) http://skxjz.nhri.cn/article/id/279
    [32] ZHANG Q, XIAO M Z, SINGH V P, et al. Copula-based risk evaluation of hydrological droughts in the East River basin, China[J]. Stochastic Environmental Research and Risk Assessment, 2013, 27(6):1397-1406. doi:  10.1007/s00477-012-0675-9
    [33] AYANTOBO O O, LI Y, SONG S B, et al. Probabilistic modelling of drought events in China via 2-dimensional joint Copula[J] Journal of Hydrology, 2018, 559:373-391. doi:  10.1016/j.jhydrol.2018.02.022
    [34] SERINALDI F, BONACCORSO B, CANCELLIERE A, et al. Probabilistic characterization of drought properties through copulas[J]. Physics and Chemistry of the Earth:Part A/B/C, 2009, 34(10/11/12):596-605.
    [35] CHEN L, SINGH V P, GUO S L, et al. Drought analysis using Copulas[J]. Journal of Hydrologic Engineering, 2013, 18(7):797-808. doi:  10.1061/(ASCE)HE.1943-5584.0000697
    [36] CHANG J X, LI Y Y, WANG Y M, et al. Copula-based drought risk assessment combined with an integrated index in the Wei River basin, China[J]. Journal of Hydrology, 2016, 540:824-834. doi:  10.1016/j.jhydrol.2016.06.064
    [37] KANG L, JIANG S W. Bivariate frequency analysis of hydrological drought using a nonstationary standardized streamflow index in the Yangtze River[J]. Journal of Hydrologic Engineering, 2019, 24(2):05018031. doi:  10.1061/(ASCE)HE.1943-5584.0001749
    [38] GUO Y, HUANG S Z, HUANG Q, et al. Copulas-based bivariate socioeconomic drought dynamic risk assessment in a changing environment[J]. Journal of Hydrology, 2019, 575:1052-1064. doi:  10.1016/j.jhydrol.2019.06.010
    [39] GU L, CHEN J, YIN J B, et al. Projected increases in magnitude and socioeconomic exposure of global droughts in 1.5 and 2℃ warmer climates[J]. Hydrology and Earth System Sciences, 2020, 24(1):451-472. doi:  10.5194/hess-24-451-2020
    [40] 闫宝伟, 郭生练, 陈璐, 等.长江和清江洪水遭遇风险分析[J].水利学报, 2010, 41(5):553-559. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201005008.htm

    YAN B W, GUO S L, CHEN L, et al. Flood encountering risk analysis for the Yangtze River and Qingjiang River[J]. Journal of Hydraulic Engineering, 2010, 41(5):553-559. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201005008.htm
    [41] 陈璐, 郭生练, 张洪刚, 等.长江上游干支流洪水遭遇分析[J].水科学进展, 2011, 22(3):323-330. http://skxjz.nhri.cn/article/id/1735

    CHEN L, GUO S L, ZHANG H G, et al. Flood coincidence probability analysis for the upstream Yangtze River and its tributaries[J]. Advances in Water Science, 2011, 22(3):323-330. (in Chinese) http://skxjz.nhri.cn/article/id/1735
    [42] HUANG K D, CHEN L, ZHOU J Z, et al. Flood hydrograph coincidence analysis for mainstream and its tributaries[J]. Journal of Hydrology, 2018, 565:341-353. doi:  10.1016/j.jhydrol.2018.08.007
    [43] PENG Y, SHI Y L, YAN H X, et al. Coincidence risk analysis of floods using multivariate Copulas:case study of Jinsha River and Min River, China[J]. Journal of Hydrologic Engineering, 2019, 24(2):05018030. doi:  10.1061/(ASCE)HE.1943-5584.0001744
    [44] FENG Y, SHI P, QU S M, et al. Nonstationary flood coincidence risk analysis using time-varying copula functions[J]. Scientific Reports, 2020, 10(1):3395. doi:  10.1038/s41598-020-60264-3
    [45] BENDER J, WAHL T, MVLLER A, et al. A multivariate design framework for river confluences[J]. Hydrological Sciences Journal, 2016, 61(3):471-482. doi:  10.1080/02626667.2015.1052816
    [46] 甘富万, 黄宇明, 张华国, 等.河流支流入汇口处水利工程防洪设计水位研究:以珠江流域西江桂平航运枢纽为例[J].湖泊科学, 2020, 32(1):198-206. https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX202001019.htm

    GAN F W, HUANG Y M, ZHANG H G, et al. Design water level for flood control of water conservancy projects at the inlet and confluence of branches in river basins:a case study from Guiping shipping hub, Xijiang River, Pearl River basin[J]. Journal of Lake Sciences, 2020, 32(1):198-206. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX202001019.htm
    [47] 闫宝伟, 郭生练, 郭靖, 等.基于Copula函数的设计洪水地区组成研究[J].水力发电学报, 2010, 29(6):60-65. https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201006011.htm

    YAN B W, GUO S L, GUO J, et al. Regional design flood composition based on Copula function[J]. Journal of Hydroelectric Engineering, 2010, 29(6):60-65. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201006011.htm
    [48] 刘章君, 郭生练, 李天元, 等.梯级水库设计洪水最可能地区组成法计算通式[J].水科学进展, 2014, 25(4):575-584. http://skxjz.nhri.cn/article/id/2441

    LIU Z J, GUO S L, LI T Y, et al. General formula derivation of most likely regional composition method for design flood estimation of cascade reservoirs system[J]. Advances in Water Science, 2014, 25(4):575-584. (in Chinese) http://skxjz.nhri.cn/article/id/2441
    [49] 宋松柏, 王小军.基于Copula函数的水文随机变量和概率分布计算[J].水利学报, 2018, 49(6):687-693. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201806005.htm

    SONG S B, WANG X J. Probability distribution calculation of the sum of hydrological random variables based on Copula function approach[J]. Journal of Hydraulic Engineering, 2018, 49(6):687-693. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201806005.htm
    [50] 郭生练, 熊丰, 尹家波, 等.水库运用期的设计洪水理论和方法[J].水资源研究, 2018, 7(4):327-338. https://www.cnki.com.cn/Article/CJFDTOTAL-SZYY201804001.htm

    GUO S L, XIONG F, YIN J B, et al. Theory and method of design flood in reservoir operation period[J]. Journal of Water Resources Research, 2018, 7(4):327-338. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SZYY201804001.htm
    [51] XIONG F, GUO S L, LIU P, et al. A general framework of design flood estimation for cascade reservoirs in operation period[J]. Journal of Hydrology, 2019, 577:124003. doi:  10.1016/j.jhydrol.2019.124003
    [52] 熊丰, 郭生练, 陈柯兵, 等.金沙江下游梯级水库运行期设计洪水及汛控水位[J].水科学进展, 2019, 30(3):401-410. doi:  10.14042/j.cnki.32.1309.2019.03.010

    XIONG F, GUO S L, CHEN K B, et al. Design flood and control water levels for cascade reservoirs during operation periodin the downstream Jinsha River[J]. Advances in Water Science, 2019, 30(3):401-410. (in Chinese) doi:  10.14042/j.cnki.32.1309.2019.03.010
    [53] 郭生练, 熊丰, 王俊, 等.三峡水库运行期设计洪水及汛控水位初探[J].水利学报, 2019, 50(11):1311-1317. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201911003.htm

    GUO S L, XIONG F, WANG J, et al. Preliminary exploration of design flood and control water level of Three Gorges Reservoir in operation period[J]. Journal of Hydraulic Engineering, 2019, 50(11):1311-1317. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201911003.htm
    [54] ZSCHEISCHLER J, WESTRA S, van den HURK B J J M, et al. Future climate risk from compound events[J]. Nature Climate Change, 2018, 8(6):469-477. doi:  10.1038/s41558-018-0156-3
    [55] 刘曾美, 陈子燊.区间暴雨和外江洪水位遭遇组合的风险[J].水科学进展, 2009, 20(5):619-625. doi:  10.3321/j.issn:1001-6791.2009.05.003

    LIU Z M, CHEN Z S. Risk study of the bivariate encounter of interzone rainstorm and flood level of the outer river[J]. Advances in Water Science, 2009, 20(5):619-625. (in Chinese) doi:  10.3321/j.issn:1001-6791.2009.05.003
    [56] 涂新军, 杜奕良, 陈晓宏, 等.滨海城市雨潮遭遇联合分布模拟与设计[J].水科学进展, 2017, 28(1):49-58. doi:  10.14042/j.cnki.32.1309.2017.01.006

    TU X J, DU Y L, CHEN X H, et al. Modeling and design on joint distribution of precipitation and tide in the coastal city[J]. Advances in Water Science, 2017, 28(1):49-58. (in Chinese) doi:  10.14042/j.cnki.32.1309.2017.01.006
    [57] BEVACQUA E, MARAUN D, HOBÆK HAFF I, et al. Multivariate statistical modelling of compound events via pair-copula constructions:analysis of floods in Ravenna (Italy)[J]. Hydrology and Earth System Sciences, 2017, 21(6):2701-2723. doi:  10.5194/hess-21-2701-2017
    [58] ZELLOU B, RAHALI H. Assessment of the joint impact of extreme rainfall and storm surge on the risk of flooding in a coastal area[J]. Journal of Hydrology, 2019, 569:647-665. doi:  10.1016/j.jhydrol.2018.12.028
    [59] 康玲, 何小聪.南水北调中线降水丰枯遭遇风险分析[J].水科学进展, 2011, 22(1):44-50. http://skxjz.nhri.cn/article/id/1779

    KANG L, HE X C. Risk analysis of synchronous-asynchronous encounter probability of rich-poor precipitation in the Middle Route of South-to-North Water[J]. Advances in Water Science, 2011, 22(1):44-50. (in Chinese) http://skxjz.nhri.cn/article/id/1779
    [60] DU H H, WANG Y M, LIU K L, et al. Exceedance probability of precipitation for the Shuhe to Futuan Water Transfer Project in China[J]. Environmental Earth Sciences, 2019, 78(7):240. doi:  10.1007/s12665-019-8207-2
    [61] 吴海鸥, 涂新军, 杜奕良, 等.基于Copula函数的鄱阳湖水系径流丰枯遭遇多维分析[J].湖泊科学, 2019, 31(3):801-813. https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX201903018.htm

    WU H O, TU X J, DU Y L, et al. Multi-dimensional analysis of wetness-dryness encountering of streamflow based on the Copula function in Lake Poyang basin[J]. Journal of Lake Sciences, 2019, 31(3):801-813. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX201903018.htm
    [62] JIANG C, XIONG L H, XU C Y, et al. Bivariate frequency analysis of nonstationary low-flow series based on the time-varying copula[J]. Hydrological Processes, 2015, 29(6):1521-1534. doi:  10.1002/hyp.10288
    [63] AHMADI F, RADMANEH F, SHARIFI M R, et al. Bivariate frequency analysis of low flow using copula functions (case study:Dez River basin, Iran)[J]. Environmental Earth Sciences, 2018, 77(18):643. doi:  10.1007/s12665-018-7819-2
    [64] GYASI-AGYEI Y. Copula-based daily rainfall disaggregation model[J]. Water Resources Research, 2011, 47(7):W07535.
    [65] GYASI-AGYEI Y, MELCHING C S. Modelling the dependence and internal structure of storm events for continuous rainfall simulation[J]. Journal of Hydrology, 2012, 464/465:249-261. doi:  10.1016/j.jhydrol.2012.07.014
    [66] CALLAU PODUJE A C, HABERLANDT U. Short time step continuous rainfall modeling and simulation of extreme events[J]. Journal of Hydrology, 2017, 552:182-197. doi:  10.1016/j.jhydrol.2017.06.036
    [67] GAO C, XU Y P, ZHU Q, et al. Stochastic generation of daily rainfall events:a single-site rainfall model with Copula-based joint simulation of rainfall characteristics and classification and simulation of rainfall patterns[J]. Journal of Hydrology, 2018, 564:41-58. doi:  10.1016/j.jhydrol.2018.06.073
    [68] BÁRDOSSY A, PEGRAM G G S. Copula based multisite model for daily precipitation simulation[J]. Hydrology and Earth System Sciences, 2009, 13(12):2299-2314. doi:  10.5194/hess-13-2299-2009
    [69] CALLAU PODUJE A C, BELLI A, HABERLANDT U. Dam risk assessment based on univariate versus bivariate statistical approaches:a case study for Argentina[J]. Hydrological Sciences Journal, 2014, 59(12):2216-2232. doi:  10.1080/02626667.2013.871014
    [70] de MICHELE C, SALVADORI G, CANOSSI M, et al. Bivariate statistical approach to check adequacy of dam spillway[J]. Journal of Hydrologic Engineering, 2005, 10(1):50-57. doi:  10.1061/(ASCE)1084-0699(2005)10:1(50)
    [71] LIU Z J, XU X F, CHENG J Q, et al. Hydrological risk analysis of dam overtopping using bivariate statistical approach:a case study from Geheyan Reservoir, China[J]. Stochastic Environmental Research and Risk Assessment, 2018, 32(9):2515-2525. doi:  10.1007/s00477-018-1550-0
    [72] LI F, ZHENG Q. Probabilistic modelling of flood events using the entropy Copula[J]. Advances in Water Resources, 2016, 97:233-240. doi:  10.1016/j.advwatres.2016.09.016
    [73] 闫宝伟, 郭生练, 刘攀, 等.基于Copula函数的径流随机模拟[J].四川大学学报(工程科学版), 2010, 42(1):5-9. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201001003.htm

    YAN B W, GUO S L, LIU P, et al. Streamflow simulation based on Copula function[J]. Journal of Sichuan University (Engineering Science Edition), 2010, 42(1):5-9. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201001003.htm
    [74] CHEN L, SINGH V P, GUO S L, et al. Copula-based method for multisite monthly and daily streamflow simulation[J]. Journal of Hydrology, 2015, 528:369-384. doi:  10.1016/j.jhydrol.2015.05.018
    [75] CHEN L, QIU H Y, ZHANG J H, et al. Copula-based method for stochastic daily streamflow simulation considering lag-2 autocorrelation[J]. Journal of Hydrology, 2019, 578:123938. doi:  10.1016/j.jhydrol.2019.123938
    [76] MADADGAR S, MORADKHANI H, GAREN D. Towards improved post-processing of hydrologic forecast ensembles[J]. Hydrological Processes, 2014, 28(1):104-122. doi:  10.1002/hyp.9562
    [77] 刘章君, 郭生练, 李天元, 等.贝叶斯概率洪水预报模型及其比较应用研究[J].水利学报, 2014, 45(9):1019-1028. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201409002.htm

    LIU Z J, GUO S L, LI T Y, et al. Comparative study of Bayesian probabilistic flood forecasting models[J]. Journal of Hydraulic Engineering, 2014, 45(9):1019-1028. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201409002.htm
    [78] MADADGAR S, MORADKHANI H. Improved Bayesian multimodeling:integration of copulas and Bayesian model averaging[J]. Water Resources Research, 2014, 50(12):9586-9603. doi:  10.1002/2014WR015965
    [79] LIU Z J, GUO S L, XIONG L H, et al. Hydrological uncertainty processor based on a copula function[J]. Hydrological Sciences Journal, 2018, 63(1):74-86. doi:  10.1080/02626667.2017.1410278
    [80] XIONG L H, YU K X, GOTTSCHALK L. Estimation of the distribution of annual runoff from climatic variables using copulas[J]. Water Resources Research, 2014, 50(9):7134-7152. doi:  10.1002/2013WR015159
    [81] FAN Y R, HUANG G H, BAETZ B W, et al. Development of a copula-based particle filter (CopPF) approach for hydrologic data assimilation under consideration of parameter interdependence[J]. Water Resources Research, 2017, 53(6):4850-4875. doi:  10.1002/2016WR020144
    [82] WANG W G, FU J Y. Global assessment of predictability of water availability:a bivariate probabilistic Budyko analysis[J]. Journal of Hydrology, 2018, 557:643-650. doi:  10.1016/j.jhydrol.2017.12.068
    [83] YU K X, ZHANG X, LI P, et al. Probability prediction of peak break-up water level through vine copulas[J]. Hydrological Processes, 2019, 33(6):962-977. doi:  10.1002/hyp.13377
    [84] WANG W Z, DONG Z C, LALL U, et al. Monthly streamflow simulation for the headwater catchment of the Yellow River basin with a hybrid statistical-dynamical model[J]. Water Resources Research, 2019, 55(9):7606-7621. doi:  10.1029/2019WR025103
    [85] 郭爱军, 黄强, 王义民, 等.基于Archimedean Copula函数的流域降雨-径流关系变异分析[J].水力发电学报, 2015, 34(6):7-13. https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201506002.htm

    GUO A J, HUANG Q, WANG Y M, et al. Detection of variations in precipitation-runoff relationship based on Archimedean Copula[J]. Journal of Hydroelectric Engineering, 2015, 34(6):7-13. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SFXB201506002.htm
    [86] XU P C, WANG D, SINGH V P, et al. A two-phase copula entropy-based multiobjective optimization approach to hydrometeorological gauge network design[J]. Journal of Hydrology, 2017, 555:228-241. doi:  10.1016/j.jhydrol.2017.09.046
    [87] LEI X H, TAN Q F, WANG X, et al. Stochastic optimal operation of reservoirs based on copula functions[J]. Journal of Hydrology, 2018, 557:265-275. doi:  10.1016/j.jhydrol.2017.12.038
  • [1] 郭生练, 熊立华, 熊丰, 尹家波.  梯级水库运行期设计洪水理论和方法 . 水科学进展, 2020, 31(5): 734-745. doi: 10.14042/j.cnki.32.1309.2020.05.010
    [2] 王红川, 杨氾, 潘军宁.  考虑极端天气时港珠澳大桥人工岛设计波浪要素——II.设计波要素校核方法 . 水科学进展, 2019, 30(6): 902-907. doi: 10.14042/j.cnki.32.1309.2019.06.014
    [3] 周平, 周玉良, 金菊良, 蒋尚明, 吴成国.  水文双变量重现期分析及在干旱中应用 . 水科学进展, 2019, 30(3): 382-391. doi: 10.14042/j.cnki.32.1309.2019.03.008
    [4] 徐翔宇, 许凯, 杨大文, 郦建强.  多变量干旱事件识别与频率计算方法 . 水科学进展, 2019, 30(3): 373-381. doi: 10.14042/j.cnki.32.1309.2019.03.007
    [5] 张宇亮, 蒋尚明, 金菊良, 吴志勇, 周玉良.  基于区域农业用水量的干旱重现期计算方法 . 水科学进展, 2017, 28(5): 691-701. doi: 10.14042/j.cnki.32.1309.2017.05.006
    [6] 涂新军, 杜奕良, 陈晓宏, 柴苑苑, 卿颖.  滨海城市雨潮遭遇联合分布模拟与设计 . 水科学进展, 2017, 28(1): 49-58. doi: 10.14042/j.cnki.32.1309.2017.01.006
    [7] 陈子燊, 黄强, 刘曾美.  变化环境下西江北江枯水流量联合分布分析 . 水科学进展, 2015, 26(1): 20-26. doi: 10.14042/j.cnki.32.1309.2015.01.003
    [8] 刘章君, 郭生练, 李天元, 徐长江.  梯级水库设计洪水最可能地区组成法计算通式 . 水科学进展, 2014, 25(4): 575-584.
    [9] 刘文琨, 裴源生, 赵勇, 肖伟华.  区域气象干旱评估分析模式 . 水科学进展, 2014, 25(3): 318-326.
    [10] 陈璐, 郭生练, 周建中, 熊立华, 闫宝伟.  长江上游多站日流量随机模拟方法 . 水科学进展, 2013, 24(4): 504-512.
    [11] 谢 华, 罗 强, 黄介生.  基于三维copula函数的多水文区丰枯遭遇分析 . 水科学进展, 2012, 23(2): 186-193. doi: CNKI: 32.1309.P.20120224.2003.015
    [12] 梁忠民, 郭彦, 胡义明, 汤春义.  基于copula函数的三峡水库预泄对鄱阳湖防洪影响分析 . 水科学进展, 2012, 23(4): 485-492. doi: CNKI:32.1309.P.20120614.2158.008
    [13] 康玲, 何小聪.  南水北调中线降水丰枯遭遇风险分析 . 水科学进展, 2011, 22(1): 44-50.
    [14] 陈璐, 郭生练, 张洪刚, 闫宝伟, 刘心愿.  长江上游干支流洪水遭遇分析 . 水科学进展, 2011, 22(3): 323-330.
    [15] 陆桂华, 闫桂霞, 吴志勇, 康燕霞.  基于copula函数的区域干旱分析方法 . 水科学进展, 2010, 21(2): 188-193.
    [16] 赖锡军.  水动力学模型与集合卡尔曼滤波耦合的实时校正多变量分析方法 . 水科学进展, 2009, 20(2): 241-248.
    [17] 刘曾美, 陈子燊.  区间暴雨和外江洪水位遭遇组合的风险 . 水科学进展, 2009, 20(5): 619-625.
    [18] 肖义, 郭生练, 刘攀, 熊立华, 方彬.  分期设计洪水频率与防洪标准关系研究 . 水科学进展, 2008, 19(1): 54-60.
    [19] 方彬, 郭生练, 肖义, 刘攀, 武见.  年最大洪水两变量联合分布研究 . 水科学进展, 2008, 19(4): 505-511.
    [20] 谢华, 黄介生.  两变量水文频率分布模型研究述评 . 水科学进展, 2008, 19(3): 443-452.
  • 加载中
计量
  • 文章访问数:  25
  • HTML全文浏览量:  2
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-06-02
  • 网络出版日期:  2020-11-16
  • 刊出日期:  2021-01-30

Copula函数在水文水资源中的研究进展与述评

doi: 10.14042/j.cnki.32.1309.2021.01.015
    基金项目:

    国家自然科学基金资助项目 51909112

    国家自然科学基金资助项目 51879192

    作者简介:

    刘章君(1991-), 男, 江西吉安人, 博士, 主要从事水文分析与计算研究。E-mail:liuzhangjun@whu.edu.cn

    通讯作者: 郭生练, E-mail:slguo@whu.edu.cn
  • 中图分类号: P333;TV12

摘要: Copula函数是一种灵活构造多变量联合分布和处理多变量问题的有效工具,在水文水资源领域应用广泛;随着研究的广度和深度不断推进,显示出其良好的适用性和不可替代的优越性。重点综述了近10年来Copula函数在水文事件多变量频率分析、水文事件遭遇组合分析、水文随机模拟、水文模型与预报以及其他问题中的最新应用研究进展与不足。从应用条件、参数估计、不对称性、拟合优选、尾部特性、多变量重现期选取、成果抽样误差等方面阐述了Copula函数应用需要注意的问题。展望了未来Copula函数在水文水资源领域应用中的研究重点和发展方向,为今后的研究提供借鉴和参考。

English Abstract

刘章君, 郭生练, 许新发, 许世超, 成静清. Copula函数在水文水资源中的研究进展与述评[J]. 水科学进展, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
引用本文: 刘章君, 郭生练, 许新发, 许世超, 成静清. Copula函数在水文水资源中的研究进展与述评[J]. 水科学进展, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
LIU Zhangjun, GUO Shenglian, XU Xinfa, XU Shichao, CHENG Jingqing. Application of Copula functions in hydrology and water resources: a state-of-the-art review[J]. Advances in Water Science, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
Citation: LIU Zhangjun, GUO Shenglian, XU Xinfa, XU Shichao, CHENG Jingqing. Application of Copula functions in hydrology and water resources: a state-of-the-art review[J]. Advances in Water Science, 2021, 32(1): 148-159. doi: 10.14042/j.cnki.32.1309.2021.01.015
  • 相关性普遍存在于水文事件内部属性间, 例如暴雨的历时和强度、洪水的洪峰和洪量、干旱的历时和烈度等。此外, 还存在由于外部原因(如空间联系、因果关系等)引起的相关性, 如洪水的遭遇和地区组成、降雨和径流的关系等[1]。深入了解和研究水文现象中的这些相关性规律, 对于提高水文分析计算成果的精度和可靠性大有裨益。

    早期的多变量水文分析计算研究主要采用边缘分布相同的多变量分布模型(如多变量正态分布、多变量对数正态分布和多变量Gumbel等)和Meta-Gaussian分布模型[2]。这些传统模型通常基于变量之间的线性相关关系而建立, 对于非线性、非对称的随机变量难以很好地描述;另外部分模型假定变量服从相同边缘分布, 同时对相关性大小也存在一定的限制, 影响了其适用性和成果准确性。实际上, 水文变量的相关性非常复杂, 包括线性相关和非线性相关;水文变量边缘分布既可能服从正态分布, 也可能服从偏态分布。因此, 如何构建非线性、非正态条件下的多变量联合概率分布无疑具有很大的挑战性。Copula函数是多变量联合分布构建理论与方法的重大突破, 可以将边缘分布和相关性结构分开来研究, 且对边缘分布类型没有任何限制, 形式灵活多样, 可以描述变量间非线性、非对称的相关关系[3]。1959年Sklar在统计学中提出Sklar定理, 直到20世纪90年代末才开始广泛应用于金融、财经、保险、精算和风险分析等领域。2003年, 意大利水文学者de Michele[4]将Copula函数引入水文水资源领域, 首次建立了边缘分布为广义Pareto分布的降雨强度和降雨历时的联合概率分布。2004年, 熊立华和郭生练[5]采用Gumbel-Hougaard Copula函数构建了长江流域某站点洪峰和洪量的联合分布。近十几年来, Copula函数引起了国内外水文水资源界的高度关注和兴趣, 研究的广度和深度不断推进。大量的研究和应用实践表明, Copula函数作为构造多变量联合分布的一种有效工具, 完整地保留了变量间相关性信息, 具有很强的灵活性和良好的适用性[6]

    本文重点综述近10年来Copula函数在水文事件多变量频率分析、水文事件遭遇组合分析、水文随机模拟、水文模型与预报以及其他问题中的最新研究进展, 对未来发展方向进行展望, 以期为继续开展Copula函数在水文水资源中的应用研究提供借鉴和参考。

    • Copula函数应用基本流程包括: ①根据实际研究需要选定随机变量;②确定各随机变量的边缘分布;③计算相关系数及显著性检验;④确定各备选Copula函数的参数;⑤ Copula函数的拟合检验与优选;⑥推求联合概率分布和条件概率分布;⑦结合实际问题开展相关后续研究。各步骤相关具体方法可以参见文献[2-3]。表 1列出了Copula函数在水文水资源应用领域的主要研究问题及变量。

      表 1  Copula函数在水文水资源应用领域的主要研究问题及相应变量

      Table 1.  Main research problems and relevant variables of Copula function applied in the field of hydrology and water resources

      研究方面 研究问题 研究变量
      水文事件多变量频率分析 降雨事件 降雨强度、降雨总量和降雨历时等
      洪水事件 洪峰、洪量和洪水历时等;不同分期内的洪水;发生时间和量级
      干旱事件 干旱烈度、干旱强度、干旱历时、干旱间隔时间等
      水文事件遭遇组合分析 洪水遭遇 干支流不同站点的洪水
      设计洪水地区组成 水库断面和区间洪水
      复合洪涝事件 暴雨和外江水位;暴雨、河洪和潮位
      丰枯遭遇 不同站点/地区的降水或径流
      水文随机模拟 降雨随机模拟 降雨总量、降雨历时、降雨强度、降雨间隔时间等
      洪水随机模拟 洪峰、洪量和洪水过程线形状
      径流随机模拟 相邻时刻径流;不同站点径流
      水文模型与预报 水文预报不确定性评估 预报流量和实际流量
      水文模型输入、参数、输出相关性模拟 水文模型输入、参数、输出
      其他问题 水量水质组合 水量和水质指标
      降雨—径流关系变异诊断 降雨和径流
      水文气象站网优化设计 不同站点日降雨量
      水库随机优化调度 相邻时段的水库入库径流
      河流水沙关系 径流与泥沙
    • 降雨、洪水、干旱等很多水文事件往往具有多个特征属性, 考虑这些属性的相关性并进行多变量频率分析可以更加科学合理、全面地评估水文事件的整体特征。水文事件多变量频率分析研究的是发生在同一个站点/区域的同一个水文事件的多个特征属性之间的相关性。

    • 降雨事件通常包括降雨强度、降雨总量和降雨历时等特征属性, 分析这些特征属性变量的相关性, 并建立联合概率分布, 有助于更加全面地认识降雨事件的变化规律。Zhang和Singh[7]采用4种Archimedean Copula函数分别对降雨强度和降雨总量、降雨历时和降雨总量及降雨强度和降雨历时进行两变量联合分析;许月萍等[8]比较分析了4种Copula函数在模拟不同历时降雨量的联合分布的差异及其对降雨频率分析的影响; Vandenberghe等[9]基于比利时105 a的10 min降雨数据, 采用二维Copula函数拟合了暴雨特征量之间的相关性结构;Balistrocchi和Bacchi[10]采用Copula函数分析了暴雨总量和暴雨历时之间的统计相关性;Jun等[11]通过Copula函数建立暴雨强度和暴雨历时的两变量联合分布来计算场次暴雨事件的发生概率;Li等[12]分别采用时变GEV分布和Copula函数来模拟年最大降雨总量和降雨强度的边缘分布和联合分布, 进行非一致性条件下两变量降雨频率分析。

      工程实践中, 中小流域或城市排水暴雨洪水计算时常遇到流量资料不足情况, 采用Copula函数构建暴雨强度、暴雨历时的联合分布, 分析确定当地雨量-频率-历时关系或暴雨强度公式, 是提高用暴雨资料推求设计洪水精度的一条有效途径。然而, 目前降雨事件多变量频率分析的系列抽样方法多采用年最大取样法, 可能会遗漏一些数值较大的暴雨, 造成小重现期部分明显偏小, 需要加强年多次取样法在降雨事件多变量频率分析的应用研究, 以扩大暴雨信息量, 增强样本代表性。

    • 洪水事件属性包括洪峰、洪量和洪水历时等, 考虑这些属性的相关性并进行多变量频率分析可以更好地掌握洪水事件特征, 并应用于设计洪水估算和和风险分析。这方面的研究最初主要是分析计算洪峰、洪量和洪水历时等特征属性的联合和条件概率分布、联合和条件重现期。Zhang和Singh[13]基于Copula函数建立洪峰、洪量和历时的二维联合概率分布, 并对条件重现期和联合重现期进行了探讨。在此基础上, 如何在多变量框架下进行洪水联合设计值估算成为学者们关心的热点问题, 核心是多变量重现期的定义和联合设计值的选取[14]。多变量重现期以“或”(“OR”)、“且”(“AND”)和Kendall重现期应用最为广泛[15-18], 最近兴起的结构荷载重现期也逐步引起了学者们的注意, 联合设计值选取的常用方法是最可能组合法[19-20]。针对变化环境导致的洪水非平稳性特征, 也有学者开展了非一致性多变量洪水频率分析, 包含洪峰、洪量边缘分布的非一致性, 以及洪峰、洪量相关结构的非一致性[21-24]。同时, 多变量洪水分析计算的不确定性评估也是一个前沿方向, 不仅涉及边缘分布的不确定性, 还要考虑相关性结构(Copula函数)的不确定性。两者又均涉及到样本抽样、线型(类型)选取和参数估计不确定性, 评估综合不确定性非常复杂[25-27]。此外, 为了在水库控制运用中充分考虑洪水的季节性变化规律, 提高防洪兴利综合效益, 基于Copula函数的分期设计洪水和水库管运洪水也是学者们关注的焦点[28-30]

      目前洪水事件的多变量频率分析仍然面临一些挑战。首先, 多变量重现期的概念比单变量情形复杂, 尚无公认和通用的定义方法, 基本原则是要根据具体研究问题的物理背景而定。其次, 对于同一防洪标准对应无穷多种洪峰、洪量组合, 如何科学合理地选择设计值已成为另一个关键问题, 当前应用最广泛的是最可能组合法。如何在多变量框架下进行洪水联合设计值估算及风险评估仍然是一个有争议的问题, 需要进一步研究完善。最后, 由于洪水通常关心的是上尾部稀遇频率的设计值, 多变量洪水分析模型比单变量模型具有更大的抽样不确定性, 导致模型成果不确定性很大。

    • 描述干旱事件的特征属性主要有干旱烈度、烈度峰值、干旱强度、干旱历时和干旱间隔时间等, 对多个干旱特征属性进行多变量频率分析可以更加科学地评估干旱风险。Copula函数早期主要被应用在对干旱特征值的二元分析上, 其中以干旱历时和干旱烈度的联合概率分布应用最广[31-33]。为了更全面地反映干旱事件的多维特征, 多元Copula函数也被应用于干旱问题的多特征分析中[34-36]。近年来, 针对气候变化和人类活动导致的干旱非一致性问题, 有学者开展了变化环境下的多变量干旱频率分析研究, 通常采用时变Copula函数、滑动窗口或全球气候模型(GCMs)气候情景分析等方法[37-39]

      尽管目前采用的干旱指标很多, 但这些指标大多数只是从某一方面(降水亏缺或地表径流减少等)来评估旱情严重程度, 未能全面地评估干旱风险, 难以反映旱情在地区间和季节上的差异。因此, 通过融合多种旱情指标, 构建能反映干旱发生机理和过程的综合指标仍然是今后研究的重点之一。此外, 受气候变化和人类活动的影响, 亟待开展基于非一致水文序列资料的水文干旱事件重现期和设计值估算及不确定性评估方法研究, 以适应变化环境带来的影响。

    • 水文事件遭遇组合分析主要针对具有物理联系和相关性的不同时空水文事件变量之间, 包括洪水遭遇中干支流不同站点的洪水遭遇、设计洪水地区组成中水库断面和区间洪水的遭遇组合、复合洪涝事件中多个致灾因子遭遇、不同站点或地区的丰枯遭遇等。

    • 洪水遭遇是指干流与支流或支流与支流的洪峰在相差较短的时间内到达同一河段的现象。闫宝伟等[40]和陈璐等[41]分别采用二维和四维Copula函数研究了最大洪水发生时间及发生量级遭遇的风险特征;Huang等[42]提出了一种干支流洪水过程遭遇概率分析方法;Peng等[43]研究了考虑历史洪水信息的洪水遭遇概率;Feng等[44]采用时变Copula函数探讨了气候变化影响下的非一致性洪水遭遇风险。此外, 河流交汇处的洪水遭遇组合和工程水文设计问题也受到学界关注[45-46]

      目前在洪水遭遇概率计算方面做了很多有益的探索, 但大多数研究仅分析了年最大洪水遭遇情况, 没有考虑其他有可能危及保护目标安全的较大洪水遭遇情况, 导致计算的遭遇概率与实际情况不完全相符。另外, 如何将洪水遭遇分析成果应用于干支流水库群防洪补偿调度、错峰调度, 以发挥其实际价值, 也是值得继续探讨的课题。

    • 设计洪水地区组成主要研究当防洪控制断面发生设计洪水时, 其上游水库断面、区间流域等各分区来水的分配及遭遇组合的形式。闫宝伟等[47]提出的基于Copula函数的最可能地区组成法统计意义明确、方案数唯一, 应用前景良好[48], 也有学者提出了基于Copula函数的改进离散求和法克服了传统离散求和法需要进行变量独立性转换的问题[49]。郭生练等[50]在此基础上提出了水库运行期设计洪水计算理论和方法, 并在金沙江下游梯级水库、三峡水库运行期设计洪水及汛控水位推求中得到了应用[51-53]

      洪水的空间分布情况对水库群防洪系统设计与调度运用具有重要影响, 设计洪水地区组成分析是不可回避的关键问题。在水库群防洪联合调度中采用运行期设计洪水及汛控水位作为依据, 并与实时洪水预报结合起来, 是一个具有潜力和实用价值的方向。然而, 目前实际中应用最广泛的方法是将设计断面设计洪量分配给上游各分区以研究水库调洪影响, 这类方法假定水库调蓄后设计断面的洪水与天然情况同频率, 没有充分考虑洪水地区组成的随机特性和水库调洪影响, 此外还隐含假设了设计断面洪水的洪峰、洪量同频率, 其合理性需要进一步探究。

    • 洪涝灾害有时是由多个致灾因子遭遇组成的复合水文事件[54]。例如, 中国的许多沿江平原地区和三角洲河网地区大多地势低平, 暴雨涝水能否顺利排除不仅取决于设计暴雨标准, 还受到外江水位的极大制约, 设计暴雨和外江水位的遭遇问题本质上就是一个两变量联合分布问题。此外, 沿海地区城市的外江通常为感潮河段, 其洪水位既受上游河道洪水的作用, 也受下游河口潮位的影响, 任一时刻的洪水位都是上游河道洪水和下游河口潮位共同作用的结果。因此, 可采用Copula函数构造区间暴雨与外江水位、暴雨与潮位、上游洪水与河口潮位等多个致灾因子的联合概率分布进行组合风险分析[55-58]

      多致灾因子影响下的复合洪涝事件分析最终目的是要研究多个荷载作用下的荷载效应变量的概率分布问题, 而当前的大多数研究局限于多个荷载本身的联合概率。例如, 感潮河段洪水位是上游河道洪水和下游河口潮位共同作用的结果, 堤防工程结构破坏机制是河段洪水位超过了某一指定高程(如堤防设计洪水位)。因此, 在今后复合洪涝事件研究中应试图建立荷载效应变量与荷载变量的功能函数, 基于随机变量函数的概率分布理论推算出荷载效应变量的概率分布, 并据此进行工程设计和风险评估。

    • 丰枯遭遇实质是气象水文变量的多区域组合问题, 通常涉及到调水工程水源区和受水区的降水或径流组合、多条河流的径流组合或枯水流量组合等形式。分析水源区和受水区的丰枯遭遇对于跨流域水资源调配、多水源联合调度等都具有重要意义, 康玲和何小聪[59]、Du等[60]分别采用Copula函数研究了南水北调中线工程和山东省沭河至付疃河调水工程的丰枯遭遇问题。此外, 很多学者在多条河流的径流组合、不同站点的枯水流量组合等方面作了有益的探索[61-62], Ahmadi等[63]还进行了基于时变Copula函数的两变量非一致性枯水频率联合分析。

      目前丰枯遭遇研究采用时间尺度一般为全年、汛期或非汛期, 对于降雨径流的变化描述不够精细, 后续研究应尽可能细化到月尺度以增强实用性。此外, 当前的丰枯遭遇研究通常采用天然径流量资料, 然而由于人类活动(水库调蓄、土地利用等)影响了天然径流的时空分布特征, 原有的天然丰枯遭遇规律可能已经发生变化, 需要加强考虑人类活动影响的丰枯遭遇研究, 更加真实地反映遭遇情况。

    • 单站随机模拟重点是描述水文事件多个特征属性间的相关性结构或水文时间序列的自相关特性, 核心是如何构造有效的联合分布或条件分布。多站随机模拟则要在此基础上考虑站点间的空间相关性, 同时模拟时间、空间2个相关结构。Copula函数是刻画变量相关性的有效工具, 近年来也应用于水文随机模拟研究中。从随机模拟的对象来看, 主要可以分为降雨随机模拟、洪水随机模拟和径流随机模拟。

    • 降雨事件随机模拟首先基于Copula函数模拟降雨特征属性, 如降雨总量、降雨历时、降雨强度、降雨间隔时间的相关性结构, 然后模拟降雨的时程分配模式, 最后对两者进行耦合。大多数研究主要针对单站点降雨事件随机模拟, 时间尺度上可以实现日、小时和分钟[64-67。Bárdossy和Pegram[68]利用多维Copula函数同时模拟各站点间降雨量和降水发生的空间相关性, 开发了多站点降雨随机模拟模型。

      近些年来, 由于人类活动的强烈干扰, 径流形成的条件发生显著的变化, 破坏了水文资料系列的一致性, 过去的实测资料无法反映未来的变化规律, 在人类活动影响后新的径流形成条件下通过暴雨资料推求设计洪水是一种可行途径。因此, 很有必要将建立的能考虑暴雨特征量间内在相关性的随机降水模型与流域水文模型进行有机耦合, 得到长系列的流量过程, 从中统计得到洪水频率曲线, 推求变化环境下的设计洪水。目前这方面的研究刚刚起步, 尚未进入实用阶段。

    • 洪水事件随机模拟一般分为2个模块, 即基于Copula函数的洪峰、洪量的联合随机模拟和洪水过程线形状的模拟[69]。de Michele等[70]采用二维Copula函数随机模拟水库洪峰、洪量进而生成大量入库洪水过程线, Liu等[71]在此基础上建立了大坝洪水漫顶风险率计算模型。也有学者直接采用三维Copula函数构建了洪量、洪峰和洪水历时的联合概率分布随机模拟洪水过程线[72]

      目前基于Copula函数的洪水随机模拟研究主要集中在单站洪水过程模拟, 对于多站洪水同步模拟研究较少, 而多站洪水随机模拟是水库群系统地区组成分析和运行期设计洪水及汛控水位计算的有力工具, 可以较好地克服水库调蓄后设计断面的洪水与天然情况同频率假定的问题, 建议加强这方面的研究。

    • 径流随机模拟的核心是描述径流时间序列的自相关结构, 因此构造径流随机模拟模型的重点就转移到如何构造有效的联合分布或条件分布。闫宝伟等[73]建立了基于Copula函数的一阶非平稳时间序列季节性模型;Chen等[74]在此基础上进一步提出了可以考虑二阶自相关的日流量随机模拟方法;Chen等[75]将基于Copula函数的径流随机模拟推广到多站点的月或日流量序列。

      径流随机模拟技术的正确使用, 将有助于得到比传统方法更为可靠的结果, 但是目前基于Copula函数的径流随机模拟研究大都局限于径流样本随机生成理论方法的研究, 应积极推动随机模拟生成径流系列在水资源工程规划设计和管理运用实际中的应用。另外, 由于在径流随机模拟中所使用的模型缺乏物理机制, 存在较大的模型不确定性, 亟待加强模型与参数的不确定性研究工作。

    • 水文模型与预报中大量存在变量相关性问题, 例如评估水文预报结果的不确定性时需要模拟预报和实际流量之间的相关性, 水文模型的输入、参数和输出之间也存在相关性[76]。应用Copula函数解决水文模型预报相关问题主要包括水文预报不确定性评估和水文模型输入、参数和输出之间相关性模拟等。

    • 采用Copula函数构建预报流量和实际流量的联合概率分布, 然后基于给定预报流量时实际流量的条件概率分布分析评估水文预报不确定性[76], 或用条件概率分布推导描述贝叶斯先验分布和似然函数, 再基于得到的后验分布对水文预报不确定性进行评价[77]。Madadgar和Moradkhani[78]提出了耦合Copula函数的贝叶斯模式平均(BMA)方法;刘章君等[77, 79]提出了基于Copula函数的单变量和多变量水文不确定性处理器(HUP)。

      近几年, 基于Copula函数的水文预报不确定性评估方法取得了快速发展, 应用非常灵活, 可以实现非线性和非正态条件下的概率水文预报。但仍存在有待进一步深入研究的问题:一方面, 基于Copula函数的概率水文预报应逐步实现从单站点、单变量向多站点、多变量转变, 以期考虑水文变量间的时空相关性并降低预报不确定性;另一方面, 要加强概率水文预报结果在洪水风险预警、水库调度风险决策中的应用, 实现预报与风险决策过程的有机结合, 真正体现概率水文预报的价值和效益。

    • 通过利用Copula函数建立水文模型各输入、各参数或输入和参数之间的联合分布, 捕捉其相关性结构, 得到的结果比假定它们相互独立更加合理。Xiong等[80]基于Copula函数和概率分布推导法得到了年径流的概率分布。Fan等[81]提出了基于Copula函数的粒子滤波方法用于水文资料同化。Wang和Fu[82]基于Copula函数建立了二元概率Budyko方法框架, 并应用于全球水资源可利用量的概率评估。此外, 也有学者以洪水或径流作为因变量, 采用Copula函数直接建立因变量与自变量(影响因子)的联合分布, 基于条件概率分布进行洪水或径流模拟预测[83-84]

      不同水文模型各输入、各参数或输入和参数之间的相关性各具特点, 建模时采用能准确捕捉相关性特征的Copula函数极为重要, 如何有效地识别最优Copula函数仍有待进一步研究。采用Copula函数建立洪水(或径流)及其影响因子的条件概率分布涉及的主要影响因子通常有多个, 便于构建高维Copula函数的藤Copula模型得到普遍应用, 然而目前藤Copula函数的结构和参数不确定性分析考虑较少, 有必要加强这方面的研究。

    • 除了上述几个方面外, 国内外学者还在降雨—径流关系变异诊断[85]、水文气象站网优化设计[86]、水库随机优化调度[87]等方面进行了探索性研究, 拓展了Copula函数在水文水资源领域的应用范围。

    • 虽然Copula函数是描述变量间相关性结构和多变量水文分析的一种有效数学工具, 但要有一个客观的认识, 既要认识到其优越性, 也不能不加分析判断盲目使用, 在实际应用中应注意以下几个方面的问题。

      (1)  Copula函数的应用条件。Copula函数不是一种万能工具, 并不是所有关于相关性和多元分布的问题都可以套用。应用Copula函数的前提条件是确保所分析的变量之间确实存有物理上的相关性, 否则会得到错误的结论。

      (2)  Copula函数的参数估计。以水文领域中最常采用的阿基米德Copula函数族为例, 目前二维情形可以采用Kendall秩相关系数法;对于三维及以上的Copula函数, 大多采用极大似然法。依据边缘分布估计方法的不同, 极大似然方法又可分为全参数、两阶段和半参数极大似然法, 其中半参数法不受边际分布的影响, 稳健性较好, 推荐使用。也有学者提出贝叶斯方法和非参数方法进行参数估计, 但目前实际中应用较少。

      (3)  对称型和非对称型Copula函数。对称型Copula函数要求变量两两之间的相关性结构相同或非常相似, 要求变量之间的相关系数非常接近, 但多变量水文分析实际中两两变量的相关性结构往往并不相同。对于具有不对称相关性的高维随机变量, 单参数难以真实反映其复杂的不对称相关结构。非对称型Copula函数具有更加灵活的参数和结构形式, 更适用于拟合高维的水文随机变量, 在水文多变量分析计算实际应用中更加合适。

      (4)  Copula函数的选择。Copula函数本身不具有水文物理基础, 不能从物理意义上确定各随机变量的理论联合分布模型。实际应用中采用不同Copula函数, 计算结果有较大差异, 因此选取合适的Copula函数尤其关键。实际工作中一般选取与已有实测资料配合最优的Copula函数, 常用的拟合优度准则有均方根误差准则(RMSE值最小)、Akaike信息准则(AIC值最小)以及Cramér-von Mises检验法(P值大于0.05通过拟合检验, 统计量Sn最小作为最优Copula函数)。

      (5)  Copula函数的尾部特性。不同水文事件表现出的尾部相关特性差异较大, 不考虑尾部相关性可能会高估或低估风险。例如, 通常认为Gumbel-Hougaard Copula函数能够拟合洪水的上尾相关性, Clayton Copula函数可以拟合枯水的下尾相关性, 若没有考虑尾部相关性而误用了其他的Copula函数, 对于洪水和枯水都会造成高估或低估灾害事件的风险。因此, 拟合优度只是选择Copula函数的必要条件, 而不是充分条件, 应根据水文事件特点选择可以准确描述尾部相关特性的Copula函数。

      (6)  多变量重现期的选取。单变量情况下, 重现期和设计值存在唯一的对应关系, 计算方法清晰简单。在多变量情况下, 由于涉及到各变量的多种组合, 目前对多变量重现期的定义也相应有很多种, 最常用的有“或”(“OR”)、“且”(“AND”)、Kendall和生存Kendall重现期以及最近兴起的结构荷载重现期。这些多变量重现期表征的物理意义和危险区域的划分方式不同, 实际应用中尚没有公认、通用的选择方法, 基本原则是要根据具体研究问题的物理背景而定, 选择能够正确表征水文失事机理的多变量重现期, 作为水文事件联合设计和风险评估的依据。

      (7)  Copula函数模型成果的抽样误差。相同样本容量下, 基于Copula函数的多变量水文分析计算模型比单变量模型具有更大的抽样不确定性, 相同模型成果许可误差下需要的样本容量更大(如单变量模型需要几十年资料, 相应多变量模型则可能需要几百年)。因此, 虽然Copula函数在理论上是构建多变量水文分析计算模型的一种有效工具, 但在应用层面上, 仍然难以摆脱工程实践中样本容量较小的制约, 应该重视模型成果不确定性分析。

    • 自Copula函数理论与方法引入水文水资源领域以来, 许多学者对其进行研究和探讨, 在水文事件多变量频率分析、水文事件遭遇组合分析、水文随机模拟、水文模型与预报及其他问题中的应用研究迅速发展。经过十几年在水文水资源领域中的应用, 已证实Copula函数是一种灵活构造多变量联合分布和处理多变量问题的有效工具, 并显示出良好的适用性和不可替代的优越性。鉴于水文现象的高度复杂性, 仍存在不少需要继续深入探讨的问题, 建议进一步开展以下几个方面的研究工作。

      (1)  继续完善Copula函数的理论和方法。首先, 高维Copula函数的构建及其参数估计方法仍是Copula理论研究中的难点之一, 藤Copula可以将高维Copula分解为若干二元Copula, 在应用中降低计算量, 使得高维问题得到简化, 是今后研究的重点。其次, 目前应用最广泛的阿基米德和椭圆类Copula函数都有其自身的局限性, 后续研究可以拓展新的Copula函数簇, 以适应复杂多变的水文现象。最后, 与其他水文分析计算方法相结合也是未来的一种发展趋势, 如Copula函数与熵理论、贝叶斯理论、马尔可夫链等结合。

      (2)  深入开展非一致性条件下的Copula函数建模方法研究。由于气候变化及人类活动(水库调蓄、土地利用和植被变化)的影响, 应用Copula函数建模时单个水文变量的边缘分布和水文变量之间的相关结构都可能发生变化。时变Copula函数模型的求解难度增大, 仍需在实用化方面进一步研究。因此, 建议深入开展非一致性条件下的Copula函数建模方法研究, 以期得到更加科学实用的结论。

      (3)  建立Copula函数多变量模型的不确定性分析框架。基于Copula函数建立的多变量模型往往比单变量模型具有更大的不确定性, 不仅涉及边缘分布的不确定性, 还要考虑相关性结构(Copula函数)的不确定性。两者又均涉及到样本抽样、线型(类型)选取和参数估计不确定性, 评估综合不确定性非常复杂, 目前相关研究仍处于起步阶段, 亟待建立一套多变量模型不确定性分析框架。

      (4)  推进Copula函数在水文水资源领域更深、更广层次的应用研究。深化Copula函数在水文事件多变量频率分析、水文事件遭遇组合分析、水文随机模拟、水文模型与预报等现有领域的应用, 改进和完善当前研究存在的问题和不足。此外, 进一步拓展Copula函数的应用范围, 例如水文系列插补延长、山洪灾害动态临界雨量计算、水位—流量关系、相应水位(流量)预报等。

参考文献 (87)

目录

    /

    返回文章
    返回