nav emailalert searchbtn searchbox tablepage yinyongbenwen piczone journalimg journalInfo journalinfonormal searchdiv searchzone qikanlogo popupnotification paper paperNew
2025, 12, v.56 67-86
基于Copula函数和蒙特卡洛法的黄河上游多站径流相关分析与随机模拟
基金项目(Foundation): 国家重点研发计划(2023YFC3206700); 国家自然科学基金项目(52479032,U2243232)
邮箱(Email): lixiang@iwhr.com;
DOI: 10.13928/j.cnki.wrahe.2025.12.006
摘要:

【目的】径流具有随机性、非平稳性、时间连续性和空间异质性。考虑径流时空二维相关特征,建立了一种多站径流序列重构模型,能够模拟径流时空多变化场景。【方法】优选了合适的边缘分布和Copula函数,建立了相同站点不同时段、不同站点相同时段的径流相关关系,进而采用蒙特卡洛法中的逆变换抽样随机模拟多站径流序列。以黄河上游为例,开展了唐乃亥和兰州两大控制性水文站月尺度径流相关分析与随机模拟。【结果】结果表明,随着模拟次数增加,模拟径流序列与实测径流序列的一致性不断提高;当模拟次数为500次时,两站月径流平均相对误差最大为2.07%和3.28%;当模拟次数达到10 000次时,两站各月径流平均相对误差分别低于0.34%和0.37%。【结论】综上,基于Copula函数和蒙特卡洛法的径流序列重构模型能够推广应用于更多站点,为延长少资料区水文序列,建立多时空径流遭遇场景下具有水力与电力联系的水库群调度规则,以及流域水资源配置方案等提供可靠的数据支持。

Abstract:

[Objective]Runoff exhibits characteristics of randomness, non-stationarity, temporal continuity, and spatial heterogeneity. This paper considers the spatiotemporal two-dimensional correlation characteristics of runoff and establishes a multi-station runoff sequence reconstruction model that can simulate various spatiotemporal scenarios of runoff.[Methods]This study optimally selected appropriate marginal distributions and Copula functions to establish the runoff correlations at the same station over different time periods, as well as across different stations at the same time period. Subsequently, the inverse transform sampling method from Monte Carlo simulations was used to randomly simulate multi-station runoff sequences. The method was applied to the Upper Yellow River, where correlation analysis and stochastic simulation of monthly runoff were conducted for the two control hydrological stations at Tangnaihai and Lanzhou.[Results]The results indicated that as the number of simulations increased, the consistency between the simulated and observed runoff series improved. When the number of simulations was 500, the maximum mean relative error(MRE) of monthly runoff at Tangnaihai and Lanzhou stations reached 2.07% and 3.28%, respectively. When the number of simulations reached 10,000, the MRE for all monthly runoff at both stations was lower than 0.34% and 0.37%, respectively.[Conclusion]The runoff series reconstruction model based on the Copula function and Monte Carlo method is applicable to more stations. It provides reliable data support for extending hydrological sequences in areas with sparse data, establishing operational rules for reservoir systems with hydraulic and electrical connections, and for formulating basin water resource allocation schemes under various spatiotemporal runoff scenarios.

参考文献

[1] 左其亭.黄河流域生态保护和高质量发展研究框架[J].人民黄河,2019,41(11):1-6.ZUO Q T.Research framework for ecological protection and high-quality development in the Yellow River Basin[J].Yellow River,2019,41(11):1-6.

[2] 莫淑红,甄晓彤,陈丽丽,等.基于abcd和Budyko模型的佳芦河干湿季径流变化归因分析[J].中国水利水电科学研究院学报(中英文),2023,21(6):501-510.MO Shuhong,ZHEN Xiaotong,CHEN Lili,et al.Attribution analysis of dry and wet season runoff changes in the Jialu River based on abcd model and Budyko model[J].Journal of China Institute of Water Resources and Hydropower Research,2023,21(6):501-510.

[3] WANG Z,LI M W,ZHANG X,et al.Prediction of long-term future runoff under multi-source data assessment in a typical basin of the Yangtze River[J].Journal of Hydrology:Regional Studies,2024,56:102053.

[4] NAZEER A,MASKEY S,SKAUGEN T,et al.Simulating the hydrological regime of the snow fed and glaciarised Gilgit Basin in the Upper Indus using global precipitation products and a data parsimonious precipitation-runoff model[J].Science of The Total Environment,2022,802:149872.

[5] 李国英.为以中国式现代化全面推进强国建设、民族复兴伟业提供有力的水安全保障:在2024年全国水利工作会议上的讲话[J].水利发展研究,2024,24(1):1-10.LI G Y.Improved water security for China’s efforts to build itself into a stronger country and rejuvenate the Chinese nation on all fronts by pursuing Chinese modernization:Speech at the 2024 National Water Conservancy Work Conference[J].Water Resources Development Research,2024,24(1):1-10.

[6] 李国英.深入贯彻落实党的二十大精神扎实推动新阶段水利高质量发展:在2023年全国水利工作会议上的讲话[J].水利发展研究,2023,23(1):1-11.LI G Y.Thoroughly implement the spirit of the 20th National Congress of the Communist Party of China and solidly promote the high-quality development of water conservancy in the new stage:Speech at the National Water Conservancy Work Conference in 2023[J].Water Resources Development Research,2023,23(1):1-11.

[7] 吴丰昌.我国水体污染控制与治理成效、科技支撑与展望[J].水利发展研究,2023,23(12):1-8.WU F C.Effectiveness,scientific and technological support,and prospects for water pollution control and management in China[J].Water Resources Development Research,2023,23(12):1-8.

[8] 刘扬,王立虎,杨礼波,等.基于QR-RBL的短期径流预测研究[J].水利水电技术(中英文),2023,54(1):87-94.LIU Y,WANG L H,YANG L B,et al.Study on short-term runoff forecast based on QR-RBL[J].Water Resources and Hydropower Engineering,2023,54(1):87-94.

[9] 吴昊昊,倪晋,曾兰婷.NPDM及改进模型在淮河流域的月径流模拟研究[J].水利水电技术(中英文),2024,55(3):113-126.WU H H,NI J,ZENG L T.Simulation of monthly runoff in Huaihe River Basin based on NPDM and improved model[J].Water Resources and Hydropower Engineering,2024,55(3):113-126.

[10] JIA B J,ZHOU J Z,TANG Z Y,et al.Effective stochastic streamflow simulation method based on Gaussian mixture model[J].Journal of Hydrology,2022,605:127366.

[11] BRUNNER M,GILLELAND E.Stochastic simulation of streamflow and spatial extremes:A continuous,wavelet-based approach[J].Hydrology and Earth System Sciences,2020,24(8):3967-3982.

[12] THALLI MANI S,KOLLURU V,AMAI M,et al.Enhanced streamflow simulations using nudging based optimization coupled with data-driven and hydrological models[J].Journal of Hydrology:Regional Studies,2022,43:101190.

[13] 王权森,李安强,卢程伟.基于条件重采样的多站洪水过程随机模拟方法研究[J].水电能源科学,2023,41(7):98-101.WANG Q S,LI A Q,LU C W.Research on stochastic simulation method of multi-station flood series based on conditional resampling[J].Water Resources and Power,2023,41(7):98-101.

[14] 陈大春,何英,曹伟.基于非参数核密度估计模型的乌鲁木齐河月径流随机模拟[J].水文,2014,34(2):66-70.CHEN D C,HE Y,CAO W.Urumqi River monthly runoff stochastic simulation based on non-parameter kernel density estimation model[J].Journal of China Hydrology,2014,34(2):66-70.

[15] 康玲,郭金垒,周丽伟,等.基于Copula函数的多站洪水过程随机模拟研究[J].水文,2024,44(5):32-39.KANG L,GUO J L,ZHOU L W,et al.Stochastic simulation study of the flood process based on multi-dimensional copula function[J].Journal of China Hydrology,2024,44(5):32-39.

[16] 王慧颖,刘舒,李敏,等.基于Copula函数的区域洪水频率分析[J].水利水电技术(中英文),2023,54(8):30-42.WANG H Y,LIU S,LI M,et al.Regional flood frequency analysis based on copula function[J].Water Resources and Hydropower Engineering,2023,54(8):30-42.

[17] 姚礼双,彭杨,廉博勋,等.三峡区间年最大降雨与上游洪水遭遇风险研究[J].水力发电学报,2023,42(12):48-60.YAO L S,PENG Y,LIAN B X,et al.Study on coincidence risks for annual maximum rainfall in Three Gorges region and flood events in Upper Yangtze River[J].Journal of Hydroelectric Engineering,2023,42(12):48-60.

[18] 洪柳天骄,陈伏龙,王统霞,等.新疆气象干旱过程识别及多维变量联合分布特征分析[J].水利水电技术(中英文),2023,54(8):1-15.HONG L,CHEN F L,WANG T X,et al.Identification of drought process and analysis of multi-dimensional variable joint distribution in Xinjiang[J].Water Resources and Hydropower Engineering,2023,54(8):1-15.

[19] NELSEN R B.An Introduction to Copulas[M].New York:Springer New York,2006.

[20] WANG P Y,ZHANG J,WANG M L,et al.Stochastic simulation of daily runoff in the middle reaches of the Yangtze River based on SVM-Copula model[J].Systems Science & Control Engineering,2019,7(1):452-459.

[21] 杨成铭,吴辉明.基于 Copula 函数的深圳市短历时降雨和潮位联合分布概率及风险分析[J].水利发展研究,2025,25(11):204-211.YANG Chengming,WU Huiming.Joint distribution possibility of short duration rainfall and tide level in Shenzhen based on Copula function and related risk analysis[J].Water Resources Development Research,2025,25(11):204-211.

[22] KHEDUN C P,MISHRA A K,SINGH V P,et al.A copula-based precipitation forecasting model:Investigating the interdecadal modulation of ENSO’s impacts on monthly precipitation[J].Water Resources Research,2014,50(1):580-600.

[23] LATIF S,SIMONOVIC S P.Trivariate probabilistic assessments of the compound flooding events using the 3-D fully nested Archimedean (FNA) copula in the semiparametric distribution setting[J].Water Resources Management,2023,37(4):1641-1693.

[24] 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.

[25] HUANG K,FAN Y R.Parameter uncertainty and sensitivity evaluation of copula-based multivariate hydroclimatic risk assessment[J].Journal of Environmental Informatics,2021,38(2):131-144.

[26] 王颖,于忱,王红瑞,等.基于条件混合三维Copula函数的多支流干流年最大流量模型研究[J].应用基础与工程科学学报,2021,29(1):64-77.WANG Y,YU C,WANG H R,et al.Study of flood coincidence probability for multi-tributary based on mixed conditional 3D-copula function[J].Journal of Basic Science and Engineering,2021,29(1):64-77.

[27] 卢韦伟,陈璐,周建中,等.基于多元分布函数的区域洪水频率分析[J].水文,2015,35(5):6-10.LU W W,CHEN L,ZHOU J Z,et al.Regional flood frequency analysis based on multivariate distribution function[J].Journal of China Hydrology,2015,35(5):6-10.

[28] 潘国勇,张钰荃,郝曼秋,等.基于Copula函数的武澄锡虞区雨潮组合风险分析[J].水电能源科学,2021,39(8):22-25.PAN G Y,ZHANG Y Q,HAO M Q,et al.Risk analysis of rain tide combination in wuchengxiyu district based on copula function[J].Water Resources and Power,2021,39(8):22-25.

[29] 林乐曼,周倩倩.基于Copula函数的并联水库群洪水地区组成方法[J].水电能源科学,2021,39(8):89-93.LIN L M,ZHOU Q Q.Research on flood region composition method of parallel reservoirs based on copula function[J].Water Resources and Power,2021,39(8):89-93.

[30] SADEGH M,RAGNO E,AGHAKOUCHAK A.Multivariate Copula Analysis Toolbox (MvCAT):Describing dependence and underlying uncertainty using a Bayesian framework[J].Water Resources Research,2017,53(6):5166-5183.

[31] CHEN J,BRISSETTE F P,LUCAS-PICHER P,et al.Impacts of weighting climate models for hydro-meteorological climate change studies[J].Journal of Hydrology,2017,549:534-546.

[32] BRUNNER M I,SEIBERT J,FAVRE A C.Bivariate return periods and their importance for flood peak and volume estimation[J].WIREs Water,2016,3(6):819-833.

[33] WEN Y L,ZHOU L W,KANG L,et al.Drought risk analysis based on multivariate copula function in Henan Province,China[J].Geomatics,Natural Hazards and Risk,2023,14(1):2223344.

[34] 王芸,赵鹏祥.黄河流域极端气候事件的时空变异特征研究[J].西北林学院学报,2021,36(3):190-196.WANG Y,ZHAO P X.Temporal and spatial variation characteristics of extreme climate events in the Yellow River Basin[J].Journal of Northwest Forestry University,2021,36(3):190-196.

[35] 陶长琪,杨海文.空间计量模型选择及其模拟分析[J].统计研究,2014,31(8):88-96.TAO C Q,YANG H W.Spatial econometric model selection and its simulation analysis[J].Statistical Research,2014,31(8):88-96.

[36] SALVADORI G,DURANTE F,DE MICHELE C,et al.A multivariate copula-based framework for dealing with hazard scenarios and failure probabilities[J].Water Resources Research,2016,52(5):3701-3721.

[37] 郭生练,闫宝伟,肖义,等.Copula函数在多变量水文分析计算中的应用及研究进展[J].水文,2008,28(3):1-7.GUO S L,YAN B W,XIAO Y,et al.Multivariate hydrological analysis and estimation[J].Journal of China Hydrology,2008,28(3):1-7.

[38] 王丽霞.概率论与随机过程理论、历史及应用[M].北京:清华大学出版社,2012.WANG L X.Probability Theory and Stochastic Processes:Theory,History and Applications[M].Beijing:Tsinghua University Press,2012.

基本信息:

DOI:10.13928/j.cnki.wrahe.2025.12.006

中图分类号:TV121

引用信息:

[1]米子甲,李想,沈延青,等.基于Copula函数和蒙特卡洛法的黄河上游多站径流相关分析与随机模拟[J].水利水电技术(中英文),2025,56(12):67-86.DOI:10.13928/j.cnki.wrahe.2025.12.006.

基金信息:

国家重点研发计划(2023YFC3206700); 国家自然科学基金项目(52479032,U2243232)

引用

GB/T 7714-2015 格式引文
MLA格式引文
APA格式引文
检 索 高级检索