在气候变化的情况下,研究未来流域各水均衡要素及地下水的补给变化情况,对未来流域的地下水资源的开发利用非常重要。本文利用SWAT模型对西安黑河流域未来气候变化情景下,水均衡要素及地下水补给情况进行模拟。结果表明:1) SWAT模型对西安黑河流域水资源模拟适用性良好;2) 在未来的50年,流域的地表径流量与地下径流量逐年减小,对降水的响应程度都较高;蒸散发量和降水量逐渐增加,受温度和降水等多种气候因素的影响;3) 未来50年的地下水补给量整体上有所下降,降水量在短期内是影响地下水补给量的主要因素;而在长时间系列内,温度等其他因素也会较强的影响地下水补给量。 Studying the impact of future climate change on water equilibrium factors and groundwater re-charge is very important for the exploitation and utilization of groundwater resources in the basin. Based on SWAT (Soil and Water Assessment Tool) model, the variation of water equilibrium factors and groundwater recharge in Xi’an Heihe River Basin under the future climate change scenarios is simulated. The result indicated that: 1) The SWAT model is of higher adaptability and accuracy in water resources simulation of the Heihe River Basin. 2) In the next 50 years of the basin, the amounts of surface and groundwater runoffs decrease year by year, and they have a higher degree of response to precipitation; evapotranspiration and precipitation gradually increase, and they are affected by various climatic factors such as temperature and precipitation. 3) The amount of groundwater recharge shows a decline over the next 50 years. Precipitation is the main factor affecting groundwater recharge in the short run; however, other factors such as temperature also strongly affect groundwater recharge in a long time series.
王磊,霍艾迪*,王菊翠,吉玉洁,耿继祥,曹哲,李聪怡
长安大学环境科学与工程学院,陕西 西安
收稿日期:2018年5月4日;录用日期:2018年5月22日;发布日期:2018年5月29日
在气候变化的情况下,研究未来流域各水均衡要素及地下水的补给变化情况,对未来流域的地下水资源的开发利用非常重要。本文利用SWAT模型对西安黑河流域未来气候变化情景下,水均衡要素及地下水补给情况进行模拟。结果表明:1) SWAT模型对西安黑河流域水资源模拟适用性良好;2) 在未来的50年,流域的地表径流量与地下径流量逐年减小,对降水的响应程度都较高;蒸散发量和降水量逐渐增加,受温度和降水等多种气候因素的影响;3) 未来50年的地下水补给量整体上有所下降,降水量在短期内是影响地下水补给量的主要因素;而在长时间系列内,温度等其他因素也会较强的影响地下水补给量。
关键词 :SWAT模型,西安黑河流域,地下水资源,气候变化
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
水是生命之源,是地球上生物赖以生存的根本要素,目前我国地表水资源已经不能满足地球生物的生存发展,地下水资源成为了水资源开发的又一重点。地下水资源变化受气候变化和人类活动两方面的影响。而20世纪80年代以来,许多关于全球气候变化的研究,发现在过去的近100年间,全球平均气温上升了0.3˚C~0.6˚C [
自上世纪70年代,国内外学者开始更多的关注气候变化对地下水影响。国外21世纪初开始利用水量平衡模型研究不同时空尺度的地下水补给对气候变化的响应,模型中包含了地表径流、地下水水位等水均衡要素,但并未考虑植被、土壤、坡度等下垫面条件的影响 [
目前为止,已经有学者开始建立数值模型来研究不同气候情景下地下水资源的变化,尹玉龙 [
陕西黑河是渭河右岸的支流,发源于秦岭最高峰太白山,全河长125.8公里,由西南流向东北,至周至县马召镇附近的家庄出峪后由东北的石马村入渭河,流域位于秦岭以北渭河以南的陕西省周至县境内,地理位置介于N33˚42'~N34˚13'、E107˚43'~E108˚24',总流域面积2258 km2 (见图1),黑河流域地势南高北低,自西南向东北倾斜,上游为秦岭山区,中游为秦岭山前洪积裙,下游为渭河一级阶地。
黑河流域属暖温带半干旱半湿润大陆性季风气候,具有春暖干燥、夏季燥热多雨、秋季湿润、冬季寒冷少雨雪的气候特点。多年平均气温13.2˚C,多年平均降水量700 mm,但降雨量分布不均匀,降雨主要集中在6~10月,降雨强度大,雨量集中,总降雨量占全年的70%以上,其他月份则干旱少雨,降雨量较少。
黑河的河川径流多年平均径流量8.17亿m3,主要来源于大气降水,其次有源头的冰川融水补给,占约总径流量的4% [
图1. 黑河流域气象站分布图、河网水系图及子流域划分图
酸盐岩为主,岩性多变、厚度大、范围广,是地下水的主要开采层区。中游地区地表水和地下水转换频繁:河流从山上流出,进入山前冲积扇,水资源一部分用于于农、林业的灌溉以及人畜饮用和工业用水,其余则沿河床下泄,补给地下水。被引灌的河水,相当一部分下渗补给了地下水,在细土平原一带出露成为泉水,或者再向前回归河流,或者再被引灌,连同打井抽取的地下水,再进入补给地表水 [
本研究模型建立所需的基础数据资料包括:流域的DEM、河网水系、土地利用、土壤类型、水文站以及2005~2013年坝下、小王涧、双庙、康家娅子、白羊滩、南天门、钓鱼台、金井、陈河乡、东老君岭、老水磨、沙梁子、板房子、中滩坪、麦场等15个气象站点(见图1)的日降水、最高最低气温、风速、相对湿度等气象资料(见表1)。此外,为了分析未来气候变化对地下水资源的影响,选择了在未来大气温室气体中等排放情景A1B模式下,根据IPCC气候数据中心提供的数据,采用PCA (Principal Component Analysis, PCA)与逐步线性回归相结合的统计降尺度方法模拟黑河流域各气象站未来2013~2062年的降水、最高气温、最低气温等气象数据。
SWAT模型的产汇流模拟包括大气降水、地表径流、壤中流、蒸散发、地下水、河网汇流等过程。SWAT模型首先根据DEM数据、土地利用类型、土壤性质和气象数据将研究区划分为若干个水文响应单元,然后在每个水文单元上分别建立水文物理概念模型,在每个单元上进行坡面产流、汇流计算,最后通过河网将水文响应连接起来进行整个研究区的计算 [
模型中水量平衡表达式:
S W t = S W 0 + ∑ i = 1 t ( R d a y − Q s u r f − Q g w − E a − W s e e p ) (1)
式中: S W t 为土壤最终含水量(mm); S W 0 为土壤前期含水量(mm);t为时间步长(d); R d a y 为第i天的降水量(mm); Q s u r f 为第i天的地表径流(mm); Q g w 为第i天的地下水含量(mm); E a 为第i天的蒸发量(mm); W s e e p 为第i天通过土壤剖面底部进入地下含水层的水分通量(mm)。
模型对地下水补给量的模拟采用下列表达式来计算:
W r c h r g , i = [ 1 − e ( − 1 / δ g w ) ] ⋅ W s e e p + e ( − 1 / δ g w ) ⋅ W r c h r g , i − 1 (2)
式中: W r c h r g , i , W r c h r g , i − 1 为第i天和i-1天的含水层补给量(mm); δ g w 为补给滞后时间(d); W s e e p 为第 天通过土壤剖面底部进入地下含水层的水分通量(mm)。
数据类型 | 数据属性 | 精度 | 数据来源 |
---|---|---|---|
DEM | 高程,坡度,坡长 | 网格:30 m | 国际科学数据服务平台 |
土地利用 | 土地利用分类 | 1:100万 | 西部数据共享中心 |
土壤类型 | 土壤分类及其理化数据 | 1:100万 | 西部数据共享中心 |
气象数据 | 日降水,最高、最低气温, 辐射强度,湿度,风速 | 15个(2005~2013年) | 陕西省周至县15个气象观测站 |
水文数据 | 径流(月) | 1个(2005~2013年) | 陕西省周至县陈河水文站 |
表1. 西安黑河流域SWAT模型输入数据
SWAT模型的构建需要经过以下几个步骤:1) 利用DEM数据生成流域河网:设置集水面积阈值为2000 km2;2) 根据生成的流域河网上的节点来划分子流域:划分为42个子流域(见图1),最小子流域面积541 km2,最大子流域面积12,583 km2;3) 根据土地利用数据和土壤类型数据确定水文响应单元(HRU):设置土地利用的阈值和土壤类型的阈值均为5%,划分为107个水文响应单元(HRU);4) 气象数据输入:利用流域15个站点2005~2013年的日降水,最高、最低气温,辐射强度,湿度,风速最为输入数据;5) 写入SWAT输入数据表。
SWAT-CPU (SWAT Calibration Uncertainty Procedures)可以对模型输入的参数进行敏感性分析,本文采用其中的SUFI-2 (Sequential Uncertainty Fitting, ver.2)法对模型参数进行率定和验证,该算法已经在西安黑河流域进行了应用,运行速度快,率定精度高,能够提高参数自动率定效率,使得SWAT模型自动率定的整体效果大幅提高 [
采用Nash-Suttcliffe模型效率系数Ens和决定系数R2对模型的模拟结果进行适用性评价 [
R 2 = [ ∑ i = 1 n ( Q m , i − Q m ¯ ) ( Q s , i − Q s ¯ ) ] 2 ∑ i = 1 n ( Q m , i − Q m ¯ ) 2 ∑ i = 1 n ( Q s , i − Q s ¯ ) 2 (3)
E n s = 1 − ∑ i = 1 n ( Q m , i − Q s , i ) 2 ∑ i = 1 n ( Q m , i − Q m ¯ ) 2 (4)
式中: Q m , i 为实测径流值; Q m ¯ 为平均实测径流值; Q s , i 为模拟径流值;
当评价指标达到一定的标准才能说明模型在该流域有较好的适用性 [
本研究利用陈河乡水文站的月径流数据的对流域模型参数进行敏感性分析与率定,根据敏感性评判因子t (t的绝对值越大表明参数对模型的模拟结果影响越明显),最终确定了敏感性较高的6个参数:scs径流曲线系数(CN2.mgt)、土壤有效含水量(SOL_AWG.sol)、基流衰退常数(ATPHA_RF.gw)、地下水滞后系数(GW_DELAY.gw)、地下水再蒸发系数(GW_REVAP.gw)、浅层地下水回流阀值(GWQMN.gw),说明水文条件和土壤的物理性质对模型的结果影响较大。
通过SUFI-2算法对参数进行敏感性分析和率定,得到较好的模型参数,把参数带回SWAT模型重新模拟计算,得到率定期和验证期的Nash-Suttcliffe模型效率系数Ens和相关系数R2结果,见表2,率定期月径流的相关系数R2为0.87,模型效率系数Ens为0.86;验证期月径流的相关系数R2为0.82,模型效率系数Ens为0.84,其相关系数R2和模型效率系数Ens都高于0.8。
对率定期和验证期的实测径流值和模拟径流值进行相关分析研究,相关分析曲线如图2所示,可以发现模拟值和实测值拟合程度较高,模型可信,由此证明该模型可以较好的对黑河流域进行径流模拟。
SWAT模型可以改变输入数据,将根据未来大气温室气体中等排放情景A1B模式下模拟的未来50年的气象数据输入,驱动SWAT模型,模拟未来50年各子流域的径流变化,进而研究气候变化对黑河流域水资源的影响。
黑河流域的水循环过程对水量水质状况产生重要影响,研究其水循环过程有着重要作用。本文利用SWAT模型获取了2013~2062年黑河流域年降水量、地表径流、地下径流、蒸散发变化情况(见图3)。结果表明,黑河流域的年均降水量、地表径流量、地下径流、蒸散发量分别为786.10 mm,78.78 mm,10.63 mm,628 mm,故降水是流域水资源的主要输入项,蒸散发是流域水资源的主要输出项。地表径流和地下径流在未来50年整体变化趋势都表现为下降,但幅度较小,分别约为0.1 mm/10a,0.6 mm/a,地下径流下降趋势为地表径流的六倍;降水量和蒸散发量50年间都呈现上升的趋势,分别为1.3 mm/a,0.6 mm/a。这表明流域的产汇流能力随着时间的推移逐渐减弱,而由于温度和降水的影响,流域的蒸散发增强。
将未来50年划分成三个时段,分别为2013~2025年、2016~2050和2051~2062年,第一时段和第三时段降水量相对较小,第二时段降水量相对较大,并且变化幅度较大;地表径流和地下径流出现相似的变化情况,这表明地表径流和地下径流对降水的响应程度较高,流域的水资源主要来自于大气降水。蒸
站名 | 时间 | 决定系数 | 效率系数 |
---|---|---|---|
率定期 | 20050101~20081231 | 0.87 | 0.86 |
验证期 | 20090101~20130430 | 0.82 | 0.84 |
表2. 陈河乡水文站点径流模拟结果
图2. 率定期和验证期模拟和实测径流相关分析图
图3. 黑河流域2013~2062年平均水均衡要素变化曲线
散发量在50年变化情况较复杂多变,这是因为蒸散发量是温度和降水的共同响应的结果,未来极端气候频发导致了蒸散发的复杂多变。
利用SWAT模拟的2013~2062年大气温室气体中等排放情景A1B模式气候变化情景下流域地下水补给的变化(见图4)。在未来50年的时期内,地下水的平均补给强度在15 mm左右,地下水的补给量整体上有所下降,但是不明显,下降速度为0.8 mm/a。地下水补给量和降水零变化情况基本一致,降水量增加时,地下水补给量也增加,并且可以看出当降水量小于某个值时,地下水的补给就会变为零,所以可以说降水量对地下水补给具有很重要的影响。但是地下水补给量和降水量在整体趋势上表现不一致,说明降水量在短期内(15年左右)是影响地下水补给量的主要因素,而在长时间系列内(50年以上),由于温度等其他因素的影响导致地下水补给量与降水量表现出不一致的趋势。
在2029年之前地下水补给变化幅度较小,呈现逐渐上升的趋势,这主要是因为这一时期蒸散发量逐渐减小;2029~2035年间地下水补给量为50年最少的时期,有些年份补给量甚至为0,这主要是由于这个时期为枯水期降水量较小且蒸散发较大导致的;2036~2045年这十年间,地下水补给量都呈现较高的值,这是由丰水期降水量较高引起;2046年地下水补给量开始突然的降低,随后几年补给量也呈现较低的值,这与降水量较小有关系。所以说地下水补给量在各个短时期内主要受降水量控制,部分的受到温度等其他因素的影响。
本文利用SWAT模型,并对模型参数进行率定和验证,通过设定A1B气候情景模式模拟黑河流域未来50年的水均衡要素及地下水补给量的变化情况,初步探讨气候变化对地下水补给量的影响。研究结果表明:
图4. 黑河流域2013~2062年地下水补给量变化曲线
1) SWAT模型对西安黑河陈河水文站2005~2013年月径流过程模拟精度是可以接受的,表明SWAT模型能够用于西安黑河流域水资源变化情况模拟。
2) 在未来的50年,黑河流域的年均降水量、地表径流量、地下径流量、蒸散发量分别为786.10 mm,78.78 mm,10.63 mm,628 mm,地表径流量与地下径流量逐年减小,蒸散发量和降水量逐渐增加;地表径流和地下径流对降水的响应程度都较高,蒸散发量是温度和降水等多种气候因素的共同影响的结果。
3) 西安黑河流域未来50年的地下水补给量整体上有所下降,地下水的年平均补给强度在15 mm左右。降水量是影响地下水补给量的主要因素,而在长时间系列内,温度等其他因素也会较强的影响地下水补给量。
值得注意的是,本研究中的未来气候变化是由气候模式模拟得到的,存在一定的不确定性;本文中用到的SWAT模型是通过2005~2013年实测径流数据进行校准和验证的,但必定仍然存在一定的误差,两者叠加必然会使结果受到一定的影响,但这种不确定性是所有水文模拟共同客观存在的。但是这个结果仍然在黑河流域未来时期地下水补给量对气候变化的响应研究方面有一定意义,可以为决策者管理黑河流域水资源状况提供一定的理论依据。
长安大学2016年国家级大学生创新创业训练计划(201610710077)。
王 磊,霍艾迪,王菊翠,吉玉洁,耿继祥,曹 哲,李聪怡. 气候变化下基于SWAT模型的地下水补给变化研究 Impact of Climate Change on Change of Groundwater Recharge Based on SWAT Model[J]. 自然科学, 2018, 06(03): 180-188. https://doi.org/10.12677/OJNS.2018.63028