和田河干流河道径流损失与植被生态需水关系研究
2020-12-04龙爱华邓晓雅徐海量张广朋
龙爱华,张 沛,李 江,邓晓雅,王 杰,徐海量,张 继,张广朋
(1.中国水利水电科学研究院,北京 100038;2.新疆水利水电规划设计管理局,新疆 乌鲁木齐 830000;3.石河子大学 水利建筑工程学院,新疆 石河子 832003;4.中国科学院 新疆生态与地理研究所,新疆 乌鲁木齐 830001)
1 研究背景
河道径流损失量(简称河损)是径流行进过程中由于水面蒸发、河道渗漏和洪水漫溢等因素影响下两断面之间径流量的差值,是河流水量的重要组成部分。在干旱区由于地势平坦,河道上下游比降小,河道较宽、水流较缓,河损占比较湿润地区大。河道渗漏是干旱区径流性水资源补给地下水与土壤水的主要方式,也是维持天然植被正常生长繁殖的重要生态水文过程[1-2]。因此,定量分离河损的组成,对研究干旱区植被的生态需水量具有重要意义。
和田河发源于昆仑山脉,由玉龙喀什河和喀拉喀什河汇流而成,是塔里木河的三源流之一,年均约1/4的径流下泄注入塔里木河,形成了塔里木河流域唯一贯通塔克拉玛干沙漠的绿色走廊。因此,和田河既承担着向塔里木河输水的任务,又承担着维系天然植被与生态走廊的重任,其健康持续意义重大。玉龙喀什河与喀拉喀什河汇形成和田河干流,自阔什拉什汇合口至肖塔断面约319 km,为沙漠无人区。该河段河道宽阔平坦,根据遥感测量河道宽度500~3500 m之间,平均河宽约为1750 m。由于河段缺乏整治工程,河道泥沙淤积严重,河势多变,过流能力不足,同时加上两源流上控制性水利枢纽尚不完备,防洪能力严重不足,在降水与径流量均呈增加趋势时[3-4],洪灾发生频率较高[5],因此和田河中下游河道整治一度被列为塔里木河综合治理的重要工程任务。然而,部分学者提出和田河流域尤其是干流地区分布着濒危的灰胡杨对区域环境改善和维持具有重要生态意义,灰胡杨的生长主要依靠和田河渗漏补给的地下水,一旦对和田河河道进行加固和整治,河道渗漏量减少,灰胡杨的生态需水的满足程度会受到严重影响,并且认为近年来河道两岸植被的退化是因为植被生态需水量得不到有效保障而造成的。上述疑问实质是需科学合理地回答“和田河的河道渗漏水量是多少?能否满足两岸植被生态需水?”这两个基本问题。为此,本研究以和田河干流为研究靶区,分析河段来水量与河道损耗水量之间的关系,定量分离河道损耗水量,探究河道渗漏水量与天然植被生态需水之间的定量响应过程,以期为流域河道管理和最严格的水资源管理制度提供技术支撑,为新一轮的塔里木河流域综合治理提供决策参考。
2 资料与方法
2.1 研究区概况和田河是昆仑山北坡的第二大河流,东邻克里雅河流域,西与叶尔羌河流域接壤,两河出山口多年平均(1962—2012)径流量约44.08亿m3[6]。和田河流域属于典型的暖温带大陆性气候,年降水量为5.4~89.6 mm,年蒸发潜力为2159~3137 mm[7]。流域内有5个水文站,分别为玉龙喀什河上的同古孜洛克(出山口)与艾格利牙、喀拉喀什河上的乌鲁瓦提(出山口)与吐直鲁克、以及和田河干流入塔里木河的肖塔。整个流域95%以上的人类活动和社会水资源消耗集中于玉龙喀什河和喀拉喀什河段,和田河干流无种植用地分布及人类取水活动(图1)。
2.2 数据来源本研究使用的水文数据来源于塔里木河流域管理局与和田河流域管理局,地面气象站观测数据(日值)来源于中国气象网(http://www.cma.gov.cn/),遥感数据源于 Landsat 8 OLI_TIRS 卫星,下载于地理空间数据云网站(http://www.gscloud.cn/),空间分辨率30 m。
2.3 研究方法
2.3.1 遥感图像处理 本研究选择2015年8月份Landsat 8 OLI_TIRS卫星遥感数据,根据研究区的范围,利用ENVI 5.3遥感图像处理软件,采用5、7、1波段对获得的遥感影像进行假彩色合成用于凸显植被和水体(植被表现为红色,水体为蓝色或黑色,沙漠为浅棕黄色,盐碱地则呈现由青白到青灰色调等)。随后对图像进行波段融合、图像镶嵌、图像裁剪等处理,再根据国家土地资源遥感宏观调查的分类体系,并结合统计资料、地形图及Google Earth同期地面资料,建立解译标志,应用图像处理软件,采用人机交互的监督分类方法进行解译,按照国土资源部发布的《土地利用现状分类标准》(GB/T 21010-2017)最后生成土地覆被数据。
图1 研究区概况及和田河干流假彩色遥感图
2.3.2 河道水面蒸发计算 河道水面蒸发是河道水量损耗的主要组成之一。和田河干流河道水面蒸发量采用面积定额法计算:
式中:WR为水面蒸发量,108m3;EФ为直径20 cm水面蒸发皿观测的水面蒸发量,mm,取柯坪、阿拉尔国家气象站蒸发数据的平均值;k为水面折算系数,借鉴塔里木河干流的相关研究成果,取阿拉尔站与博斯腾湖的平均值0.62;AR为河段水面面积,m2,采用遥感解译的方法获得。
2.3.3 植被生态需水计算 潜水蒸发法适合于干旱区植被生存主要依赖于地下水的情况。和田河流域平原区降水稀少,两岸多为中旱生的非地带性植被,主要依靠地下潜水维持生命。因此,本研究在干流植被景观类型遥感解译的基础上,选择潜水蒸发法来估算天然植被需水量。潜水蒸发法即某一植被类型在某一潜水位的面积乘以该潜水位下的潜水蒸发量。计算公式为[8]:
式中:W为植被需水量,108m3;Ai为植被类型i的面积,hm2;Wgi为植被类型i在某一地下水埋深时的潜水蒸发量,mm。
植被的面积(Ai)通过遥感解译获得,以2015年8月份的植被分布为基准;潜水蒸发量(Wgi)是潜水蒸发法计算植被生态需水量的关键,以阿维利扬诺夫公式计算较为常见,计算公式如下[9-10]:
式中:a、b为经验系数,其中a为0.52,b为2.51[11];hi为植被类型i的地下水埋深,m;hmax为潜水蒸发极限埋深,m;EΦ为20 cm蒸发皿蒸发量,mm。
2.3.4 河损计算 本研究采用水量平衡方程计算河道损失水量。
式中:WL为河道损失水量;WU和WD分别为研究区河道上、下断面水量;WP为河段降水量;WO为区间河段人类活动引水量;WI为河段回归水量。
3 结果与分析
3.1 河道水量耗散规律分析基于水文站监测数据,绘制月尺度(2007—2016多年平均)和田河干流径流过程线,并分析干流年际来水量与肖塔下泄水量的相关关系,结果如图2所示。其中,和田河干流来水量为艾格利牙与吐直鲁克水文站来水量之和。
图2 和田河干流来水径流过程
由图2(a)可以看出,和田河干流来水量主要集中在7—10月份,而当年11月—翌年6月河道无水经过,河床裸露。除2016年以外,当干流来水量较大时(当月超过8×108m3),肖塔下泄水量的峰值与上断面有显著的响应关系。从年际尺度来看,和田河干流来水量与肖塔下泄水量呈显著正相关关系(图2(b)),根据线性拟合曲线的结果计算可知,当干流来水量小于5.23×108m3时,肖塔断面将不会有水下泄到塔里木河。
以年为尺度计量和田河干流径流量以及河损如表1所示。
从表1可知,和田河干流2007—2016年水量损失介于4.29×108~9.45×108m3之间,单位河长河损量介于134.48×104~296.24×104m3/km之间,其中,2009年河段水量损耗最小,2010年河段河损水量最大。以2011年塔里木河流域近期综合治理基本完成为界限,对比分析2007—2011年和2012—2016年的不同时段径流量和河损量,可见近期综合治理之后干流来水量有明显增加,河损量和单位河长河损较前一阶段下降了20.32%,肖塔断面径流量增加了123.96%;河道输水效率(肖塔断面下泄水量与干流来水量之比)在近期综合治理之前为49%,之后提高到73%。
3.2 河道水量耗散定量分离本研究将和田河干流河段河道损耗分为蒸发和渗漏两部分,水面蒸发采用式(1)计算方法而得,河道渗漏水量为河损量与蒸发量的差值,具体情况见表2。
根据表2可知,2007—2016年和田河干流蒸发水量介于2.51×108~4.68×108m3之间,均值为3.80×108m3,占河损总量的58.88%;渗漏水量介于1.18×108~5.04×108m3,均值为2.66×108m3,占河损总量的41.20%。在塔里木河流域近期综合治理完成后,河道渗漏水量有所下降,而水面蒸发水量增加了21%,导致蒸发水量占河损总量的72.64%。
表1 和田河干流河段逐年核算情况
表2 和田河干流河损定量分离
将和田河干流来水量与水面蒸发和河道渗漏进行曲线拟合,如图3(a)所示,和田河干流河道水面蒸发量和来水量呈显著正相关,原因是来水量越大,河道水面越宽,导致蒸发量越大;图3(b)和图3(c)表明河道渗漏水量与当年来水量无显著相关关系,但与前一年来水量显著相关:随着前一年来水量的增加,河道渗漏水量呈现“先减后增”的变化趋势。其主要原因可能是河道垂向入渗速度随河床初始含水率的增大而减小,后期随着河道来水量增大,河道湿周增加,侧向渗入量随着增大。由此可知,维持较为稳定的年际河道来水量可有效调控河道渗漏损耗。
3.3 和田河干流天然植被生态需水根据遥感解译结果可知,和田河干流林地和高、中覆盖草地主要分布在河道两岸3 km以内,低覆盖草地分布在河道两岸7 km之内,干流天然植被空间分布情况见图4。
根据公式计算和田河干流天然植被生态需水量,选取柯坪(51720)、阿拉尔(51730)两个国家级气象站的2000—2016年20 cm蒸发皿蒸发量数据的多年平均值(1425 mm),作为研究区的潜在蒸发数据;植被影响系数借鉴参考《中国塔里木河水资源与生态问题研究》[12]中的相关成果,见表3。
根据遥感解译结果及地下水埋深情况,分析不同植被类型的面积及其潜水埋深情况,并根据计算公式,得到和田河干流生态需水量,见表4所示。
图3 来水量与蒸发量和渗漏量的关系
图4 和田河干流不同植被类型分布情况
表3 不同潜水埋深下植被影响系数
表4 和田河干流不同植被类型面积与生态需水量
根据表4可知,和田河干流生态需水量为1.37×108m3,其中高覆盖度草地和中覆盖度草地需水量最多,分别占总需水量的39%和32%,疏林地生态需水量最少,不足0.01×108m3。由此可知,和田河干流河道多年平均(2007—2016)渗漏水量约为2.65×108m3,近乎干流植被生态需水量2倍,因此依靠河道渗漏水量完全可满足和田河干流天然植被的生态需水量。
4 讨论
4.1 和田河干流生态需水量根据已有研究成果,塔里木河干流有林地、灌木林地、疏林地、高覆盖度草地、中覆盖度草地和低覆盖度草地蒸发定额分别为 1432.4、497.3、392.1、2343.8、867.3、392.1 m3/hm2[13],基于此计算可得和田河干流天然植被生态需水量为1.28×108m3。
对比分析两种方法计算的和田河干流天然植被生态需水总量结果(见表5)可知,二者总量差异较小,但不同植被类型计算结果差异较大,其中仅有灌木林地生态需水量差异大于50%,高、中、低覆盖度草地生态需水量差异分别为45.33%、38.69%和43.74%。本研究取两种方法计算的不同植被类型需水量的最大值之和及最小值之和作为干流天然植被生态需水量的范围,可知和田河干流天然植被生态需水量介于0.99×108~1.66×108m3之间。
表5 不同计算方法下和田河干流不同植被生态需水量及其比较 (单位:×108m3)
另外,有学者研究得出和田河流域天然植被生态需水量为5.52×108m3,其中玉龙喀什河、喀拉喀什河与和田河干流3条河段植被分布面积相当,且两源流植被质量优于和田河干流[8],因此可以推断,干流植被生态需水量小于和田河流域天然植被生态需水总量的1/3(1.84×108m3),也小于本研究中多年平均河道渗漏水量2.65×108m3。
根据干旱区天然植被的生态需水过程可知,4—6月份是胡杨返青期,7—9月份为胡杨落种期。胡杨落种期内实现1~2次漫溢,可以促进植被群落的更新繁育,因此,保护与修复天然植被不仅仅要总量满足,更应该在满足其生长过程需水,河道整治对径流进行年内调控可进一步满足胡杨返青期的植被需水。
4.2 河道损失与植被生态需水之间的关系据以往研究成果可知,胡杨和柽柳等荒漠河岸林植被在不同水分胁迫条件下对土壤水与地下水的利用策略虽然有所不同[14],但二者所占比重远远超过了90%,凝结水利用所占比例不足10%,且主要是浅根系草本植物[15],因此,在干旱与半干旱地区,地下水和土壤水是植被满足其生理需求的主要水分来源。地下水与土壤水的补给来源则主要是河流的侧向渗漏,并且土壤水与地下水的主要排泄方式是垂向的潜水蒸发和植被蒸腾,侧向的排泄量则相对很小。鉴于长时期来看流域地下水储量恒定不变,因此可认为由地表水转化的地下水与土壤水为潜水蒸发与植被蒸腾所消耗[2]。
在整个塔里木河流域,由于大量泥沙在河道频繁改道的过程中沉积,形成了巨厚的第四系松散沉积层,其岩性主要为粉砂、粉土及黏性土的互层,地质条件较为复杂[16],河流地表水与地下水相互转换较为频繁,难以利用地表水-地下水转换模型模拟出河道渗漏水量,因此本研究采用河损分离的方法推算河道渗漏水量,即河道损失总量减去河道水面与潜水蒸发量。基于此,本研究所得出的河道渗漏水量即地表水转化为地下水量约为河道损失量的41.2%,这与王希义等[17]在塔里木河下游采用同位素计算的地表水与地下水的转化系数43.36%相近。
前文中提及部分学者担心进行河道整治后,河道宽度缩减、河岸硬化,虽然输水能力将得到提升,水面无效蒸发量大大减小,但渗漏补给天然植被生态水量也会相应减少。根据地下水运动达西定律,河道渗漏水量的大小与河道面积、渗透系数和地表水与地下水的水力坡度成正比,为满足干流河道天然植被生态保护与修复的需水量,可以通过增加地表水与地下水的水力坡度、延长过水时间等措施增加河道补给天然植被生态需水量。另外,结合玉龙喀什河与喀拉喀什河山区控制性水利枢纽,通过开展精细化生态调度,更有利于河道两岸天然植被的生态保护与修复。
4.3 和田河干流河道损失和田河作为塔里木河干流三源流之一,向干流阿拉尔断面下泄水量的多年平均值为11.20×108m3(2007—2016),占同期阿拉尔断面径流量的25%。秦胜英等[3]核算了和田河干流的河道蒸发和渗漏量,认为干流水面多年平均蒸发量为0.64×108m3,河道渗漏量0.065×108m3,漫溢水量为6.56×108m3,这一结果与本研究有较大差异,主要是前者认为干流存在较大的漫溢水量,并且蒸发量与渗漏水量较本研究偏小。
已有大量研究表明,河道漫溢对天然植被的生长和繁育具有重要作用[18],如果河道发生大规模漫溢,必然会出现大面积新增天然植被。然而和田河流域天然植被面积近20年来不仅没有增加,反而呈现略有减少趋势[19]。因此本研究认为文献[3]中和田河干流所谓河水漫溢并非是真正的漫溢,而是由于河道较宽,水面较大,河道中间形成了诸多沙洲,遥感影像上将这种现象认为是漫溢。
以往有研究在计算河道水面蒸发过程中,计算了河道水面宽度与来水量的关系,从而计算了不同来水情况下河道水面蒸发量。然而和田河干流由于河道较宽,多数情况下,河道水量处于漫流状态,水面宽度小于河道宽度。虽然水面蒸发面积较小,但河床潜水位与河道水面水位持平,因此河床潜水蒸发与水面蒸发相差无几,可以当做水面蒸发来计算。文献[3]中水面蒸发量较小主要原因是河床潜水蒸发部分未计入蒸发项所致。
5 结论与建议
本研究基于水文站监测数据、遥感影响资料和蒸发模型,定量分析了和田河干流河道损失组成与过程,并对天然植被的生态需水量进行了核算,结果表明:(1)多年平均来水条件下,和田河干流河损量约为6.46×108m3(2007—2016),河损率为44.48%,单位河长河损量为202.60×104m3/km。从河损的组成来看,水面蒸发量为3.5×108m3,占河损量的58.88%,渗漏水量为2.65×108m3,占41.12%;(2)和田河干流天然植被的生态需水量1.37×108m3,河道渗漏水量可以满足植被生态水量需求。进行河道整治,不仅可以减少水面蒸发,提高和田河输水能力,更有利于开展生态水量调度,保护和恢复干流天然植被。
当前和田河源流之一喀拉喀什河上游已建有乌鲁瓦提水库,总库容3.47亿m3,兴利库容2.24亿m3;玉龙喀什河上游山区玉龙喀什水利枢纽工程设计库容5.36亿m3,建成后对和田河径流调节能力达到20%,基本满足流域社会经济需水与天然植被生态需水的调度要求。建议河道治理建设中设置生态闸口用于调度,适时适量开展天然植被漫溢灌溉。