2. 国家电网公司高压电器设备现场实验技术重点实验室,武汉 430077;
3. 武汉区域气候中心,武汉 430074
2. State Grid Key Laboratory of On-site Test Technology on High Voltage Power Apparatus, Wuhan 40077;
3. Wuhan Regional Climate Center, Wuhan 430074
南水北调工程是我国继三峡工程之后又一超大型水利工程[1]。该工程旨在解决我国北方地区的缺水问题,分为东、中、西三条线路,其中,中线是从位于长江北面支流汉江的丹江口水库引水,沿着宽100 m、深8 m、总长1 240 km的水渠向沿岸的华北平原供水,最后将水送到北京和天津两市,该线已于2014年12月12日正式送水。丹江口水库位于湖北省丹江口市与河南省南阳市淅川县之间,水源来自汉江及其支流丹江。作为南水北调中线工程的水源区,丹江口水库的入库流量直接影响工程的可调水量[1],水库所在地区的水资源状况[2]、径流预报[3]、水文气候特征[4-5]、人类活动影响[6]等,一直备受关注。近年来,受全球气候变暖影响,丹江口水库水源区降水年际变率增大,入库径流呈减少趋势[6],尤其是2012年以来水源区持续5 a来水偏少,2016年入库径流较常年偏少4成左右,大大增加了调水难度。因此,积极开展水库水源区降水径流研究,对开展水源区入库径流预报进而更好地服务南水北调水资源调度具有重要意义。
国内外对流域径流变化的模拟,很多学者主要通过水文模型来实现,尤其是分布式水文模型,因其既能对流域内水文过程的机理进行详细描述,又可有效模拟水文过程[7],且能获取流域大量空间信息[8]。另外,该模型可根据未来气候变化情景改变气象条件输入,实现对气候变化引起的水量平衡变化的模拟[9];还可通过改变下垫面输入,实现人类活动对径流影响贡献率的定量分析[10]。以往也有学者[1, 11]利用分布式水文模型SWAT开展过丹江口水库径流模拟,但此类模型需要建模者清晰了解研究区域水文内部过程,熟悉有关地理信息系统和遥感技术,且需要进行大量资料输入,其数据处理较为繁琐,这给模型的建立带来较大困难。同时,由于模型中对产、汇流过程及部分经验参数需要验证,而对大多数流域搜集或验证所有参数是不可能的,目前国内外虽有一些分布式水文模型基于流域信息进行参数估计,但都给模型带来了一定的不确定性。
除水文模型模拟外,传统的基于观测数据的统计分析也是常用方法,其中气候弹性模型就是经典方法之一。气候弹性概念最早由Schaake[12]提出,多用于评价气候变化对径流的影响[13]。该模型输入数据即气候要素和径流观测值,最初仅考虑降水对径流的弹性,之后Fu等[14]和Ma等[15]将年平均气温引入到气候弹性模型中,以确定径流对降水、气温的弹性系数,建立考虑降水和气温影响的气候弹性模型。该模型因其原理和方法简单、数据收集便捷且模拟精度高而被广泛用于人类活动和气候变化对流域径流影响的解析研究。如:陈莹等[16]借鉴气候弹性模型研究了闽江流域径流演变与气候变化的关系,并指出该流域径流对降水弹性系数呈下降趋势,与区域人类活动的影响加剧存在重要联系;杨娜等[17]利用气候弹性模型分析了丹江口水库入库径流减少的原因,认为气候变化是导致径流减少的主要原因。然而,利用气候弹性模型开展丹江口水库水源区径流模拟评估及未来预测,迄今尚属空白。为此,本文基于丹江口水库水源区降水、气温及水库入库径流特征,建立气候弹性模型,模拟丹江口水库入库径流,并预估其未来演变,以期为南水北调中线调水工程的风险分析与合理调度管理提供参考依据。
1 资料与方法 1.1 资料说明本文使用的资料包括: (1) 1961—2015年丹江口水库月平均入库径流资料,该资料取自丹江口水库管理处水情信息中心。(2)同期丹江口水库水源区41个国家气象站逐日降水量、气温观测资料,其观测站点分布见图 1。(3) 2012年1月1日—2061年12月31日全球气候模式BCC-CSM1-1数据(分辨率1°×1°),该数据来源于国家气候中心;其数据要素包括平均气温和降水,排放情景包括RCP2.6、RCP4.5、RCP8.5三种典型排放浓度路径,其中RCP8.5为CO2排放参考范围90百分位数的高端路径,其辐射强迫高于SRES中高排放(A2)情景和化石燃料密集型(A1F1)情景,RCP2.6和RCP4.5都是中间稳定路径,且RCP4.5的优先性大于RCP2.6。
依据Schaake[12]提出的方法,定义径流对气候要素变化的弹性系数(ε)为
$ \varepsilon = \frac{{{\rm{d}}Q/Q}}{{{\rm{d}}X/X}} $ | (1) |
式(1)中,Q为径流量;X为气候要素,如降水或气温。
将式(1)改写为离散型,即
$ \varepsilon = \frac{{\Delta {Q_i}/\bar Q}}{{\Delta {X_i}/\bar X}} $ | (2) |
故
$ \Delta {Q_i}/\bar Q = \varepsilon \Delta {X_i}/\bar X $ | (3) |
其中,Q、X分别为径流、气候要素的平均值,均采用算术平均法求取;ΔQi、ΔXi分别为径流和气候要素相对各自平均值的变化量。
因此,ε可看作是ΔXi/X和ΔQi/Q之间的线性回归系数,可采用最小二乘法进行估计,即
$ \varepsilon = \frac{{\bar X}}{{\bar Q}}\frac{{\sum {\left( {{X_i}/\bar X} \right)\left( {{Q_i}/\bar Q} \right)} }}{{\sum {{{\left( {{X_i}/\bar X} \right)}^2}} }} = {\rho _{Q - X}}{C_Q}/{C_X} $ | (4) |
式(4)中,ρQ-X是径流与降水或气温的相关系数,CQ、CX分别是径流和气候要素(降水、气温)的变差系数,计算公式如下
$ {\rho _{Q - X}} = \frac{{\sum {\left( {{Q_i}/\bar Q} \right)\left( {{X_i}/\bar X} \right)} }}{{\sqrt {\sum {{{\left( {{X_i}/\bar X} \right)}^2}} .\sum {{{\left( {{Q_i}/\bar Q} \right)}^2}} } }} $ | (5) |
$ {C_Q} = \frac{{\sqrt {\frac{{\sum {{{\left( {{Q_i}/\bar Q} \right)}^2}} }}{n}} }}{{\bar Q}} $ | (6) |
$ {C_X} = \frac{{\sqrt {\frac{{\sum {{{\left( {{X_i}/\bar X} \right)}^2}} }}{n}} }}{{\bar X}} $ | (7) |
如果同时考虑降水和气温的影响,参考文献[14]、[15],可将式(3)改写为
$ \Delta {Q_i}/\bar Q = {\varepsilon _P} \cdot \Delta {P_i}/\bar P + {\varepsilon _t} \cdot \Delta {T_i}/\bar T $ | (8) |
式(8)中,εp为径流对降水的弹性系数,εt为径流对气温的弹性系数;ΔPi、ΔTi分别为i时刻降水和气温的变化量;ΔQi为i时刻径流的变化量;p、T分别为降水和气温的平均值。故i时刻的径流量为
$ {Q_i} = \bar Q + \Delta {Q_i} $ | (9) |
通过式(4)确定径流对降水、气温的弹性系数;然后,再依次将其代入式(8)和式(9)中,即可得到径流模拟公式。
2 丹江口水库入库径流变化特征分析统计结果表明,丹江口水库水源区1961—2015年年均降水量为863.7 mm,整体以-9.1 mm·(10 a)-1的速率减少,未通过0.05显著性水平检验,其减少趋势不显著;年均气温为14.2 ℃,近50多年来以0.13 ℃· (10 a)-1的速率显著增加,增温主要发生在1990年代以后(图略)。
从1961—2015年丹江口水库入库径流年序列演变看(图 2a),该水库年均入库径流为1 131.5 m3·s-1,1960—1980年代,入库径流较为丰沛,1980年代最多为1 395.3 m3·s-1,自1990年代以来,入库径流维持在来水量较枯阶段,平均入库径流减为870.0 m3·s-1;21世纪以来,入库径流又趋于增加,但从2012年开始连续5 a水量偏枯。年均入库径流最小值(447.8 m3·s-1)出现在1999年,最大值(2 476.3 m3·s-1)出现在1964年,入库径流最小的前10 a中有8 a出现在1985年以后。利用M-K检验方法对1961—2015年丹江口水库水源区入库径流序列进行突变检验(图 2b),结果发现,年入库径流在1978、1983、1985年均存在突变点。结合径流本身演变特征和前人的相关研究结果[17-18],可以确定1985年为1961—2015年期间丹江口水库入库径流的突变年。根据这一结论,将分析时段分为1961—1985年和1986—2015年两段,分别以1961—1985年为基准期、1986—2015年为变化期,分阶段建立气候变化弹性模型。
基于1985年为1961—2015年期间丹江口水库入库径流突变年,统计丹江口水库水源区基准期(1961—1985年)和变化期(1986—2015年)入库径流、降水、气温的平均值和变异系数,同时计算降水量和气温的弹性系数以及入库径流与降水量、气温的相关系数,其结果见表 1。
从表 1中看到,丹江口水库水源区基准期平均降水量为873.8 mm,平均气温为14.0 ℃,水库年入库径流为1 288.7 m3·s-1,此三要素中入库径流的变异系数最大,平均气温的最小,降水量的居中,表明入库径流波动更为剧烈;入库径流与降水呈正相关(相关系数为0.95)、与气温呈负相关(相关系数为-0.61);基准期径流对降水和气温的弹性系数分别为1.04、-1.02,表明当丹江口水库水源区历年降水增加或减少1%时,其径流平均增加或减少1.04%,但气温增加1%时,其径流减少1.02%。
将降水、气温的变异系数与逐年观测值以及降水、气温、径流的平均值代入式(8)和式(9),可建立丹江口水库水源区基准期入库径流气候弹性模拟公式,即
$ \begin{array}{l} {Q_i} = 1288.7 + 1288.7 \cdot \left( {0.92 \cdot \frac{{\left( {{P_i} - 873.8} \right)}}{{873.8}}} \right) + \\ \;\;\;\;\;\;\left. {0.60 \cdot \frac{{\left( {{T_i} - 14} \right)}}{{14}}} \right) \end{array} $ | (10) |
针对模拟结果,采用Nash-Suttcliffe系数(Ens)和决定系数(R2)进行验证评价,其定义分别为
$ {R^2} = \frac{{{{\left[ {\sum\limits_{i = 1}^n {\left( {{Q_{\rm{o}}} - {Q_{{\rm{avg}}}}} \right)\left( {{Q_i} - {Q_{{\rm{avg}},i}}} \right)} } \right]}^2}}}{{\sum\limits_{i = 1}^n {{{\left( {{Q_{\rm{o}}} - {Q_{{\rm{avg}}}}} \right)}^2}} \sum\limits_{i = 1}^n {{{\left( {{Q_i} - {Q_{{\rm{avg}},i}}} \right)}^2}} }} $ | (11) |
$ {E_{{\rm{ns}}}} = 1 - \frac{{\sum\limits_{i = }^n {{{\left( {{Q_{\rm{o}}} - {Q_p}} \right)}^2}} }}{{\sum\limits_{i = 1}^n {{{\left( {{Q_{\rm{o}}} - {Q_{{\rm{avg}}}}} \right)}^2}} }} $ | (12) |
其中:Qo为实测值;Qi为模拟值;Qavg为实测平均值;Qavg, i为模拟平均值;n为实测数据个数。R2可用于实测值与模拟值之间的数据吻合程度评价,R2=1表示非常吻合,R2小于1时,其值越小,反映出数据吻合程度越低。Ens若为负值,说明模型模拟平均值比直接使用实测平均值的可信度更低。通常取Ens > 0.5作为径流模拟效率的评价标准,但Ens值需要在0.7以上才认为比较准确[19, 20]。当决定系数(R2)越高、Ens越大时,表明模拟效果越好。
基准期模拟入库径流和实测入库径流间的Ens=0.74,R2=0.92,达到模拟精度检验标准,依据该结果绘制1961—1985年逐年入库径流模拟值变化序列(图 3),其模拟值和观测值吻合度较好,多数年份二者较为接近,仅部分入库径流较高年份(如1964年、1983年)模拟值较观测值每秒偏少500 m3以上。总体上,所建立的气候弹性模型能较好地模拟出丹江口水库入库径流的演变趋势,但模拟值与观测值相比,其年际间变化相对平稳。
从表 1中看到,丹江口水库水源区变化期平均降水量为835.6 mm,平均气温为14.3 ℃,水库年径流量为1 000.5 m3·s-1,三要素的变异系数分别为0.87、0.56、0.96,较之基准期均有所降低,入库径流与降水、气温的相关性也有所降低,径流对降水的弹性系数略有上升,对气温的弹性系数绝对值大幅下降,比较而言,降水对径流的影响较气温更大。
据此,建立1986—2015年期间入库径流气候弹性模拟公式,即
$ \begin{array}{l} {Q_i} = 1000.5 + 1000.5 \cdot \left( {0.86 \cdot \frac{{\left( {{P_i} - 835.6} \right)}}{{835.6}}} \right) + \\ \;\;\;\;\;\;0.56 \cdot \frac{{\left( {{T_i} - 14.3} \right)}}{{14.3}} \end{array} $ | (13) |
变化期模拟入库径流和实际观测入库径流间Ens= 0.71,R2=0.91。依据绘制出的模拟与观测径流曲线看出(图 4),1980年代末至1990年代初期,入库径流变化相对较平稳,模拟值和观测值较接近;进入21世纪后入库径流年际间波动较大,尤其是近5 a入库径流年际间变化较大,模拟值和观测值差距较大,如2010、2011年模拟值较观测值每秒偏少600 m3左右。
根据1961—2015年丹江口水库入库径流观测值和模拟值距平序列可看出(图 3),入库径流模拟值序列呈现-65.9 m3·s-1·(10 a)-1的减少趋势,与观测序列-64.6 m3·s-1·(10 a)-1变化趋势接近,因为本身模拟值是基于观测值的统计分析,观测序列和模拟序列平均值一致,两序列变化趋势接近也体现出了建立气候弹性模型开展径流模拟的优势。除趋势和实测值较为一致外,对于年际间入库径流正、负距平转变模拟效果也比较好。1985年以前入库径流多为正距平,尤其1964、1983年,实际入库径流较年均值每秒高出1 000 m3以上,模拟值也显示这两年是序列中入库径流最多的年份。1960年代初期及1980年代初期入库径流偏多,尤其是1979—1985年连续6 a入库径流居于高位;1985年以后入库径流多为负距平,1999、1997年为负距平大值年,1990—1995年、2006—2009年、2002—2015年均为连续4 a及以上偏少年,模拟序列均可很好地呈现出入库径流的这种多-少演变特征。由此可知,在该研究区域,利用气候变化弹性模型虽然只考虑了气候要素,但仍能较好地模拟出丹江口水库入库径流的变化特征,借此思路,可初步开展该水库未来入库径流模拟。
利用未来RCP2.6、RCP4.5、RCP 8.5三种情景下丹江口水库上游降水、气温预估数据,建立2006—2100年该流域平均降水、气温逐年变化序列。首先,将该序列和历史1961—2015年重叠期即2006—2015年观测降水、气温数据进行对比,三种情景下模拟降水较观测值偏多,气温较观测值偏低,但整体模拟效果较好(见中国气象局2016年气候变化专项《丹江口水源区和受水区旱涝配置特征及调水适应对策研究》工作报告)。然后,计算降水(气温)观测值与模拟值的平均比值(差值),利用比值(差值)订正其余年份的降水(气温)模拟值。再将订正结果输入丹江口水库变化期入库水量模拟公式(10)得到模拟流量,利用重叠期实际观测流量对模拟流量进行差值订正,就可得到最终模拟流量序列。
图 6给出未来RCP2.6、RCP4.5、RCP 8.5三种情景下2016—2100年丹江口水库上游径流模拟序列。从中可知,其年均降水量分别为1 046.0 mm、1 072.0 mm和1 055.1mm,年均气温分别为15.4 ℃、15.8 ℃和16.8 ℃,年均径流量分别为734.3 m3·s-1、1 046.5 m3·s-1和963.7 m3·s-1。三种情景下入库径流变化趋势不一,其中,RCP4.5情景下径流量以0.8 m3·s-1· (10 a)-1的速率呈微弱的增加趋势,RCP2.6情景下以2.91 m3·s-1·(10 a)-1的速率呈增加趋势,RCP8.5情景下则以-14.2 m3·s-1· (10 a)-1的速率呈现减少趋势。
表 2给出未来RCP2.6、RCP4.5、RCP8.5三种情景下各年代丹江口水库水源区降水、气温及径流变化。从中可见,降水、气温较1986—2015年均有所增加,RCP4.5情景下降水增幅最大,RCP8.5情景下气温增幅最明显,尤其是21世纪后期增幅较前期更显著;受降水、气温影响,RCP2.6情景下丹江口入库径流21世纪各个年代均偏少270 m3·s-1以上,后期波动起伏较大,偏少幅度最小出现在2071—2080年,偏少幅度最大出现在2081—2090年;RCP4.5情景下各年代入库径流一致偏多,相比之下,2021—2040年、2071—2090年水量较丰沛,其偏多幅度保持在50 m3·s-1以上;RCP8.5情景下入库径流总体以偏少为主,2040年以前变化较复杂,2040年以后偏少幅度逐年代递增。RCP8.5情景下径流变异系数为0.062,为三种情景下最大,加之各年代径流变幅波动较大,表明该情景相比其他两种情景径流年际间波动最剧烈,不利对于水资源进行管理调度。
本文针对丹江口水库水源区水文气象特征开展研究,并建立气候弹性模型,开展径流模拟及未来径流预测。主要得出以下结论:
(1) 丹江口水库水源区近50多年来降水呈微弱减少趋势,气温以0.13 ℃·(10 a)-1的速率呈显著增加,增温主要发生在1990年代以后,入库径流以-64.3 m3·s-1·(10 a)-1的速率呈减少趋势,1960—1980年代为入库径流偏多期,1990年代以来为偏少期,21世纪以来入库径流又趋于增加。
(2) 水库入库径流与降水呈正相关,与气温呈负相关,降水对入库径流影响更明显;径流对降水和气温的弹性系数在1961—1985年期间分别为1.04、-1.02,在1986—2015年期间分别为1.05、-0.75,基于气候弹性模型建立公式,开展的径流模拟能很好地呈现水库入库径流演变特征。
(3) 未来RCP2.6、RCP4.5情景下径流量分别呈0.8 m3·s-1·(10 a)-1、2.91 m3·s-1· (10 a)-1的增加趋势,RCP8.5情景下则呈现-14.2 m3·s-1· (10 a)-1的减少趋势,且RCP8.5情景下年际间波动最剧烈,不利于水资源调度。
本研究借鉴气候弹性模型,建立丹江口水库入库径流模拟公式,实现了基于观测数据的径流序列模拟,为丹江口水库入库径流预测提供了新方法,可充分发挥气象部门资料优势,更好开展径流模拟预测,为区域水资源调度提供参考依据,后期可考虑研究区域蒸散发、辐射、风速、相对湿度等气象因子,完善模型,更好地进行径流模拟预测。
[1] |
陈杨. 基于ArcGISDesktop丹江口入库径流分析[D]. 武汉: 华中科技大学, 2010 http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D139941
|
[2] |
郭海晋, 金蓉玲. 丹江口水库上游水资源现状及变化趋势[J]. 自然资源, 1997, 34(1): 28-33. |
[3] |
陆华友, 邵东国, 郭元裕. 丹江口水库径流长期预报研究[J]. 武汉水利电力大学学报, 1996, 29(6): 6-10. |
[4] |
刘明宽, 徐元顺, 朱理国, 等. 丹江口水库流域致洪过程的水文气候特征分析[J]. 暴雨灾害, 2001, 20(4): 15-16. |
[5] |
李英, 李致家, 张利升, 等. 丹江口水库北调水量的不确定性分析[J]. 武汉大学学报(工学版), 2008, 41(3): 34-37. |
[6] |
李凌程, 张利平, 夏军, 等. 气候波动和人类活动对南水北调中线工程典型流域径流影响的定量评估[J]. 气候变化研究进展, 2014, 10(2): 118-126. |
[7] |
徐翔宇. 气候变化下典型流域的水文响应研究[D]. 北京: 清华大学, 2012 http://cdmd.cnki.com.cn/Article/CDMD-10003-1013016560.htm
|
[8] |
杨大文, 李翀, 倪广恒, 等. 分布式水文模型在黄河流域的应用[J]. 地理学报, 2004, 59(1): 143-154. DOI:10.11821/xb200401018 |
[9] |
王苗, 刘敏, 夏智宏. 基于SWAT模型未来气候变化对洪湖流域水资源影响研究[J]. 气象与环境学报, 2016, 32(4): 39-47. |
[10] |
夏智宏, 刘敏, 王苗, 等. 1990s以来气候变化和人类活动对洪湖流域径流影响的定量辨识[J]. 湖泊科学, 2014, 26(4): 515-521. DOI:10.18307/2014.0404 |
[11] |
祝诗学, 梁忠民, 戴昌军, 等. 丹江口水库流域月尺度降雨与径流预报研究[J]. 南水北调与水利科技, 2016, 14(1): 96-101. |
[12] |
Schaake J C. From climate to flow, in Climate Chang and U S Water Resources chap 8, edited by Waggoner P E[M]. New York: John Wiley, 1990, 177-206.
|
[13] |
Dooge J C I. Sensitivity of Runoff to Climate Change: A Hortonian Approach[J]. Bulletin of the American Meteorological Society, 1992, 73(12): 2003-2013. |
[14] |
Fu G, Charles S P, Chiew F H S. A two-parameter climate elasticity of streamflow index to assess climate change effects on annual streamflow[J]. Water Resources Research, 2007, 43(11): 2578-2584. |
[15] |
Ma H, Yang D W, Tan S K, et al. Impact of climate variability and human activity on streamflow decrease in the Miyun Reservoir catchment[J]. Journal of Hydrology, 2010, 389(3): 317-324. |
[16] |
陈莹, 陈兴伟. 基于弹性分析的闽江流域径流演变与气候变化关系[J]. 福建师范大学学报(自然科学版), 2011, 27(6): 101-105. |
[17] |
杨娜, 赵巧华, 闫桂霞, 等. 气候变化和人类活动对丹江口入库径流的影响及评估[J]. 长江流域资源与环境, 2016, 25(7): 1129-1134. |
[18] |
王元超, 王旭, 雷晓辉, 等. 丹江口水库入库径流特征及其演变规律[J]. 南水北调与水利科技, 2015, 13(1): 15-19. |
[19] |
盛春淑, 罗定贵. 基于AVSWAT丰乐河流域水文预测[J]. 中国农学通报, 2006, 22(9): 493-496. |
[20] |
郭晓军, 崔鹏, 朱兴华. 典型泥石流流域蒋家沟的降雨-径流模拟[J]. 水土保持通报, 2011, 31(1): 175-179. |