图1 研究区概况图
纸质出版日期:2024-06-25,
收稿日期:2024-01-14
扫 描 看 全 文
引用本文
阅读全文PDF
流域径流是变化环境的重要响应要素,在气候变化和人类活动的影响下,流域径流量已经发生了深刻变化,研究其演变及响应特征有助于提高人类对变化环境的应对能力。以北运河北京市域主要集水区为研究区,通过分析水文气象观测站近40年水文数据,分析北运河流域降水径流变化趋势,提取洪水场次;引入全球气候模式数据,利用降尺度模型进行降尺度处理,与水文模型结合对流域径流进行了预报分析。结果表明:流域历史降水趋势性变化不明显,径流量呈增大趋势;未来北运河流域(通县站)水资源的变化在气候变化影响下较为明显,随着温度升高、降水增加,流量增大,预测流量在21世纪中期达到新高,为385.1 m3/s;随着辐射强迫的持续增加,水文响应程度也相应增大,北运河流域大部分径流峰值出现在2050时期后期和2090时期。研究结果可为城市流域未来水资源利用、防洪减灾措施与双碳目标的实现路径的探索提供一定的参考。
Basin runoff is an important response factor to changing environment. Under the influence of climate change and human activities, watershed runoff has undergone profound changes. Studying its evolution and response characteristics is helpful to improve human’s ability to cope with changing environment. Taking the main catchment area of the Beijing municipal area of the North Canal as the study area, the variation trend of rainfall and runoff in the North Canal basin was analyzed by analyzing the hydrological data of hydrometeorological observation stations in the past 40 years, and the flood field was extracted. The global climate model data is introduced, the downscaling model is used for downscaling, and the runoff forecast is analyzed in combination with the hydrological model. The results show that the historical precipitation trend of the basin is not obvious, and the runoff is increasing. In the future, the change of water resources in the North Canal Basin (Tongxian Station) will be more obvious under the influence of climate change. With the increase of temperature and precipitation, the runoff will increase, and the discharge is predicted to reach a new high of 385.1m3/s in the middle of this century.With the continuous increase of radiative forcing, the degree of hydrological response increases correspondingly, and most of the peak runoff in the North Canal basin occurs in the late 2050 and 2090 periods. The research can provide some reference for the future water resource utilization, flood containment and disaster reduction measures and the realization path of dual carbon targets in urban watershed.
城市流域径流关系着城市生态环境、防灾减灾、水资源利用以及社会经济发展和稳定。百年间,全球气候发生了以变暖为主要特征的变化。探究径流对气候变化的响应,尤其是准确预估未来气候变化情景下径流变化对于水资源规划和管理以及水利高质量发展等都具有重要的科学与应用价值。研究气候变化对水文水资源的影响及响应应对措施是当今学界研究的重点问题之一。
关于水文气象要素的演变规律研究已相对成熟。1999年,叶陌生等[
国内外学者围绕气候变化下的径流预测开展了大量的工作,并形成了“未来气候情景-水文模拟-径流预测”的研究思路[
目前,利用全球气候模式数据和统计降尺度模型在城市流域开展径流预报的研究相对较少。本研究基于全球气候模式数据和更完整的原始水文数据,构建了适用于北京市域内北运河流域的HEC-HMS模型,预测了流域未来降雨径流变化情况,旨在丰富预测城市流域径流的方法和技术体系,为北运河流域应对气候变化、开展水资源战略规划和管理、研判治水方略[
本文的研究区为北运河北京市域主要集水区,北京市北运河流域分布在潮白河流域与永定河流域之间,流域面积为4 423 km2,是北京市内平原流域面积最大的流域(见
图1 研究区概况图
Fig. 1 Summary map of the study area
本研究涉及的数据包括输入模型的水文数据、气候模式数据等,数据来源及基本信息如下。
1)水文数据。收集北运河流域王家园水库、十三陵水库、马驹桥、榆林庄、温泉、高碑店、通县7个代表性水文监测站、雨量监测站点的不同历时时期的逐日、逐场次降雨径流摘录资料。主要包括:2010—2020年降雨径流摘录资料、1980—2021年逐年径流量统计资料、流域控制站控制面积等,将数据整理为小时尺度,整理为适合模型输入的格式。北运河上游为温榆河,至通县站(即通州北关闸)与通惠河相汇合后始称北运河,是北运河重要的监测站点之一,通县(北)闸上站的降水流量资料是研究北运河干流历史降水径流演变的重要资料来源(见
图2 1980—2021年通县观测站年平均降水量与径流量变化情况
Fig. 2 Changes of annual mean rainfall in Tongxian observation station from 1980 to 2021
序号 | 洪水场次 | 洪水历时/h | 洪峰流量/m3·s-1 | 洪峰模数m3/(s×km2) | 径流深/mm | 径流系数 | 峰型 |
---|---|---|---|---|---|---|---|
1 | 110623 | 30 | 143.40 | 0.09 | 3.46 | 0.11 | 单峰 |
2 | 110814 | 30 | 173.00 | 0.10 | 6.47 | 0.16 | 双峰 |
3 | 120721 | 73 | 449.52 | 0.27 | 27.99 | 0.16 | 单峰 |
4 | 140901 | 75 | 102.00 | 0.06 | 5.98 | 0.10 | 双峰 |
5 | 150721 | 82 | 106.00 | 0.06 | 9.46 | 0.18 | 双峰 |
6 | 150727 | 42 | 161.00 | 0.10 | 7.09 | 0.23 | 单峰 |
7 | 150904 | 44 | 241.00 | 0.14 | 12.01 | 0.17 | 单峰 |
8 | 160719 | 144 | 659.92 | 0.39 | 51.84 | 0.19 | 单峰 |
9 | 180716 | 54 | 440.00 | 0.26 | 27.23 | 0.38 | 单峰 |
10 | 180723 | 60 | 320.80 | 0.19 | 15.82 | 0.21 | 单峰 |
11 | 180807 | 36 | 320.00 | 0.19 | 13.94 | 0.28 | 单峰 |
12 | 190804 | 50 | 341.68 | 0.20 | 13.66 | 0.39 | 单峰 |
13 | 190909 | 47 | 342.00 | 0.20 | 18.22 | 0.31 | 单峰 |
14 | 200812 | 61 | 311.00 | 0.18 | 15.47 | 0.21 | 双峰 |
15 | 平均值 | 59 | 293.67 | 0.17 | 16.37 | 0.22 | - |
2)气候模式数据。本研究数据由CMIP6的气候模式BCC-CSM2-MR资料和北京市2000—2010年观测降雨资料构成,模式资料来源于国家气候中心。本研究选取基于SSP2(中度发展)的RCP 4.5路径下的人口密度和经济的未来预估数据(即SSP245典型情景)进行下一步研究,分析研究区未来不同时期气候变化特征。该数据的空间分辨率为1.125°×1.125°。本文通过同倍比法对比两者历史年份(2000—2010年)的降水数据,对未来气候情景数据进行了修正。
2.2.1 CMhyd降尺度模型
流域模型通常用于模拟未来气候条件对水文过程的影响[
偏差校正程序采用转换算法来调整气候模型输出。基本思想是确定观测和模拟的历史气候变量之间的偏差,以参数化偏差校正算法,该算法用于校正模拟的历史气候数据。偏差校正方法被假定为固定的,即用于当前气候条件的校正算法及其参数被假定为对未来条件也有效。因此,相同的校正算法被应用于未来的气候数据。
2.2.2 HEC-HMS模型
HEC-HMS(hydrologic engineering center-hydrologic modeling system)是一款用于模拟人工或者自然状态下流域降雨径流过程的分布式水文模型。HEC-HMS的主要功能是用于洪水预报计算,在模拟洪水过程中融合物理概念。该模型考虑了流域下垫面情况以及气候因素,可连续或场次模拟流域在各种情况下的降雨径流过程,有助于提高洪水模拟精度。作为一款适用于湿润、半湿润、干旱和半干旱地区的分布式水文模型,HEC-HMS模型具有强大的功能,不同模块中有不同的产汇流计算方法,多种计算方法结合使用,可以找到更适合研究流域的组合方法,使模拟结果更准确。基于研究区的自然地理条件,本研究选取以下4种方法组合。
1)SCS曲线数法。SCS曲线数法是一种经验计算流域降雨损失方法,由美国水土保持局基于水平衡原理率先提出,在短期降雨形成的径流过程中模拟效果更好[
(1)
(2)
累计净雨量公式转换为
(3)
由于土壤潜在最大截留S与集水区中间参数曲线数CN存在相互转换关系,S和CN转换公式如下。
CN是无量纲参数,其理论取值在0~100之间,取值大小与流域内土地利用、植被、土壤类型等下垫面情况密切相关,在降雨径流实际模拟过程中,CN可根据流域下垫面真实情况对照美国国家工程手册提供的CN值表查询,一般取值范围在40~98之间。
2)SCS单位线法。SCS单位线是一种参数化单位线模型[
(4)
峰值时间TP计算公式为
(5)
SCS单位线法只要输入集水区流域滞时tlag,与汇流时间tc的关系可用下式大致估算
(6)
3)退水曲线法。退水曲线法是一个基流指数衰减模型,源于对集水区任意时刻的基流量和初始基流可能存在某种联系的假设而建立,常用来解释集水区蓄水量的自然排水问题[
(7)
4)马斯京根法。马斯京根法是对水力学方程进行水文学简化,即将连续方程简化成水量平衡方程(河段入流减出流等于河段蓄水量变化)、将动量方程简化为槽蓄关系曲线(河段出流量与蓄水量成正比),由于使用起来较简单[
(8)
(9)
(10)
(11)
2.2.3 模拟结果评价指标
根据《水文情报预报规范》(GB/T 22482—2008,以下简称“规范”),选取洪水预报合格率(R)、洪峰流量相对误差(EQ)、峰现时差(ΔT)及纳什效率系数(Nash-Sutcliffe efficiency coefficient,NSE)作为模型模拟精度等级评定指标。评价标准为:EQ、ΔT绝对值越小,模拟结果越好;R、NSE的值越接近1,模型精度越高。各项评价指标的计算方法如下。
1)洪峰流量相对误差计算公式
(12)
2)峰现时差
(13)
3)纳什效率系数(NSE,式中以NNSE表示)计算公式
(14) 为实测流量的均值。
4)合格率计算公式
(15)
根据合格率及确定性系数可以评价模型模拟精度等级,规定如
本研究中模型涉及的参数有CN值、不透水面积比例、地表截留、洪峰滞时、初始流量、衰退系数以及河道汇流演算所需的槽蓄曲线坡度Κ和比重因子X。采用人工调参法和目标函数法结合的方式来进行参数率定,率定过程如下:CN值与洪峰滞时使用CNgrid和HEC-Geo HMS工具计算得出,其余参数首先通过人工试错的方式,结合典型流域特征以及其他学者在相似流域所作研究,在参数合理范围内逐一尝试不同的取值,不断重复该步骤直到模拟的洪水过程与实测结果较为吻合;在人工试错的基础上,使用HEC-HMS水文模型内置的单变量梯度搜索(Univarite-Gradient)自动优化方法,以均方根误差(RMS)为参数率定标准,以实测资料为目标值,通过对各参数的初设值进行多次重复叠代,最终得出调整优化后的模型参数。同样选取流域2010—2020年间的4场典型洪水,其中3场用于各子流域参数率定,1场用于模型验证,探索HEC-HMS模型在该流域的适用性,率定期的场次洪水模拟结果如
图3 率定期的场次洪水模拟
Fig. 3 Rate regular field flood simulation
流域参数率定结果见
峰现时间随着槽蓄曲线坡度Κ值的增大滞后,洪峰流量会随着比重因子X值的增大而减小。采用马斯京根法进行北运河流域(通县站)河道汇流演算,本研究最终率定的Κ值和比重因子X分别为2 h和0.2。
据流域验证期和率定期场次洪水模拟结果,将洪峰流量、峰现时间的实测值和模拟值进行对比,计算结果的相对误差及纳什效率系数见
洪水场次 | 洪峰流量/(m3·s-1) | 洪峰误差/% | 峰现时刻误差/h | Nash系数 | |
---|---|---|---|---|---|
实测 | 模拟 | ||||
20160719 | 659.9 | 753.4 | 14.19 | 4 | 0.67 |
20180723 | 320.8 | 380.0 | 18.45 | 2 | 0.64 |
20190804 | 341.7 | 376.7 | 10.26 | 1 | 0.80 |
20200812 | 311.0 | 387.6 | 24.6 | 2 | 0.68 |
由
基于CMhyd模型对气候模式数据进行降尺度和偏差校正处理,再根据校准后的HEC-HMS模型,将基于SSP2的RCP 4.5典型浓度发展路径下2015—2100年修正后的降水数据写入模型,其余参数不变,分析2030—2100年流域径流。
降尺度数据的时间范围是2010—2100年,本研究选取2000—2015年降雨数据为历史参照时段,对2030—2039年、2050—2059年和2090—2100年(分别代表本世纪近期、中期和末期)展开研究。径流模拟结果见
图4 HEC-HMS分布式水文模型径流模拟结果
Fig. 4 Results of runoff simulation by HEC-HMS distributed hydrological model
径流的变幅是气候变暖的一种响应结果,未来北运河流域(通县站)降水与径流的变化在气候变化影响下较为明显。结合SSP245未来气候情景下的径流预测结果,评估未来2030—2100年降水量,得到近期、中期与末期的日降水量最大值分别为131.38、187.71和173.39 mm。由于2030年左右温度升高相对不明显,降水量没有发生较大变化。伴随着设置情境中气温逐渐升高、辐射强迫增加,降水量也呈增大趋势。
分析未来气候情景下的径流变化可得:随着未来流域温度升高和降水增加,流量也在增大。由
径流演变趋势与决定水循环的气候条件和人类活动密切相关,如河川径流的流量和流速取决于各种环境因素,包括降水的数量和模式、气温、地形地貌和下渗能力。除了自然因素外,河流的水文状况还受到人类取用水资源、水库调蓄拦截、土地利用和覆被变化的影响。结合上述分析,综合对流域径流预报的模拟结果,针对北运河流域的防洪调度提出以下对策及建议。
1)关注北运河流域2038年、2055年、2095年关键年份洪水情况,及时通过温潮减河、运潮减河削减洪峰,做好防洪预案。根据流域径流预报结果,北运河在2038年夏季、2055年夏季、2095年左右会出现径流量高峰,其中2095年左右流域的径流量达到新高,发生5年一遇以下(含5年)洪水,需要通过温潮减河和运潮减河对北运河进行分洪。
2)充分考虑人类活动对径流的影响,做好防洪预案。北运河应重点保障城市副中心行政办公区排水排涝设施的安全运行,做到标准内洪水安全行洪。在保障行洪安全的前提下,保障通航水位,实现洪水资源化。遇中雨及以上预报或河道大流量泄水时,应停止一切在北运河河道的水上活动,确保人员安全上岸,并将游船固定在安全区域,以防冲入河道,影响下游水工建筑物安全。
3)控制碳排放,减少洪涝灾害发生。流域径流预报结果表明,未来气候变化背景下,北运河流域洪水的峰值流量都将可能增大,极端洪水的风险也相应增加,同时非汛期径流量也有一定的增加。通过综合碳排放气候情景,未来通过控制碳排放实现双碳目标,在一定程度上可降低发生洪涝灾害的风险。
叶佰生, 丁永建, 康尔泗, 等. 近40a来新疆地区冰雪径流对气候变暖的响应[J]. 中国科学(D辑:地球科学), 1999, 29(S1): 40-46. [百度学术]
TE M S, DING Y J, KANG E S, et al. Response of snow and ice runoff to climate warming in Xinjiang in recent 40 years[J]. Scientia Sinica(Terrae), 1999(S1): 40-46. [百度学术]
周婷, 于福亮, 李传哲, 等. 湄公河清盛站水文情势变化分析[J]. 水电能源科学, 2011, 29(11): 15-18,209. [百度学术]
ZHOU T, YU F L, LI C Z, et al. Analysis of hydrological regime changes of chiang saen station in Mekong River[J]. Water Resources and Power, 2011, 29(11): 15-18,209. [百度学术]
陈鑫, 刘艳丽, 刁艳芳, 等. 关河水库流域水文序列突变诊断与非一致性研究[J]. 中国农村水利水电, 2020(5): 28-32. [百度学术]
CHEN X, LIU Y L, DIAO Y F, et al. Research on the diagnosis of hydrologic series abrupt changes and inconsistency in Guanhe Reservoir Basin[J]. China Rural Water and Hydropower, 2020(5): 28-32. [百度学术]
CARLSON R F, MACCORMICK A J A, WATTS D G. Application of linear random models to four annual streamflow series[J]. Water Resources Research, 1970, 6(4): 1070-1078. [百度学术]
APOSTOLOPOULOS T K, GEORGAKAKOS K P. Parallel computation for streamflow prediction with distributed hydrologic models[J]. Journal of Hydrology, 1997, 197(1/2/3/4): 1-24. [百度学术]
张力, 王红瑞, 郭琲楠, 等. 基于时序分解与机器学习的非平稳径流序列集成模型与应用[J]. 水科学进展, 2023, 34(1): 42-52. [百度学术]
ZHANG L, WANG H R, GUO F N, et al. Integrated model and application of non-stationary runoff based on time series decomposition and machine learning[J]. Advances in Water Science, 2023, 34(1): 42-52. [百度学术]
SHAMSELDIN A Y, O’CONNOR K M, CONNOR . A real-time combination method for the outputs of different rainfall-runoff models[J]. Hydrological Sciences Journal, 1999, 44(6): 895-912. [百度学术]
HONG M, WANG D, WANG Y, et al. Mid-and long-term runoff predictions by an improved phase-space reconstruction model[J]. Environmental Research, 2016, 11(24): 1-14. [百度学术]
CHEN H, XU C Y, GUO S L. Comparison and evaluation of multiple GCMs, statistical downscaling and hydrological models in the study of climate change impacts on runoff[J]. Journal of Hydrology, 2012, 434-435: 36-45. [百度学术]
左其亭, 邱曦, 马军霞, 等. 黄河治水思想演变及现代治水方略[J]. 水资源与水工程学报, 2023, 34(3): 1-9. [百度学术]
ZUO Q T, QIU X, MA J X, et al. Evolution of flood control thought and modern flood control straregy in the Yellow River Basin[J]. Journal of Water Resources and Water Engineering, 2023, 34(3): 1-9. [百度学术]
王红瑞, 赵含, 李占玲, 等. 减碳与生物多样性若干问题评述[J]. 西北大学学报(自然科学版), 2023, 53(1): 17-24. [百度学术]
WANG H R, ZHAO H, LI Z L, et al. Review on some issues of carbon reduction and biodiversity[J]. Journal of Northwest University(Natural Science Edition), 2023, 53(1): 17-24. [百度学术]
孙才志, 段兴杰. 黄河流域水资源-能源-粮食系统生态可持续发展能力评价[J]. 人民黄河, 2023, 45(2): 85-90,96. [百度学术]
SUN C Z, DUAN X J. Ecological sustainable development capacity assessment of water-energy-food system in the Yellow River Basin[J]. Yellow River, 2023 45(2): 85-90,96. [百度学术]
袁冯, 张君枝, 王冀, 等. 气候变化背景下北京市短历时暴雨的强度及雨型变化特征[J]. 大气科学学报, 2020, 43(5): 802-809. [百度学术]
YUAN F, ZHANG J Z, WANG J, et al. Variation characteristics of intensity and pattern of short duration rainstorm in Beijing under the background of climate change[J]. Transactions of Atmospheric Sciences, 2020, 43(5): 802-809. [百度学术]
张力, 赵自阳, 王红瑞, 等. 气候变化下水文模拟不确定性若干问题讨论[J]. 水资源保护, 2023, 39(1): 109-118,149. [百度学术]
ZHANG L, ZHAO Z Y, WANG H R, et al. Discussion on several issues of uncertainty in hydrological simulation under climate change[J]. Water Resources Protection, 2023, 39(1): 109-118,149. [百度学术]
刘浩. 澜沧江流域历史径流变化归因和未来径流变化预估[D]. 上海: 华东师范大学, 2022. [百度学术]
程旭, 马细霞, 王武森, 等. HEC-HMS 模型参数区域化在河南省小流域适用性研究[J]. 水文, 2022, 42(1): 40-46,102. [百度学术]
CHENG X, MA X X, WANG W S, et al. Applicability Research of HEC-HMS model parameter regionalization in small basin of Henan Province[J]. Journal of China Hydrology, 2022, 42(1): 40-46,102. [百度学术]
RAD S, DAI J F, XU J X, et al. Lijiang flood characteristics and implication of karst storage through Muskingum flood routing via HEC-HMS, S. China[J]. Hydrology Research, 2022, 53(12): 1480-1493. [百度学术]
宋世凯, 王海龙, 郎子龙, 等. 基于HEC-HMS模型的太行山小流域山洪模拟研究:以邢台路罗川流域为例[J]. 灾害学, 2023, 38(1): 117-124. [百度学术]
SONG S K, WANG H L, LANG Z L, et al. Mountain flood simulation of small basin in Taihang Mountains using HEC-HMS model: A case study of Luluochuan River Basin[J]. Journal of Catastrophology, 2023, 38(1): 117-124. [百度学术]
WANG W C, TIAN W C, XU D M, et al.Muskingum models’ development and their parameter estimation: A state-of-the-art review[J]. Water Resources Management, 2023, 37(8): 3129-3150. [百度学术]
233
浏览量
255
下载量
0
CSCD
相关文章
相关作者
相关机构