梁麗霞，楊富元，張忠揚，馬章懷，王一博

（1.甘肅省引洮工程水資源利用中心，蘭州 730000；2.蘭州大學資源環(huán)境學院，蘭州 730000）

為改善水資源時空分布的差異性，實現(xiàn)水資源的科學高效利用，人工引水是解決水資源短缺和時空調(diào)控的主要措施之一［1］。我國水源的分布有明顯的空間差異性和時間上的不平衡性［2］。干旱地區(qū)水資源短缺和時空差異給社會經(jīng)濟發(fā)展、生態(tài)環(huán)境保護和區(qū)域高質量發(fā)展帶來了嚴重的挑戰(zhàn)［3］。水利工程是水資源有效調(diào)控和管理措施之一，可提高水資源利用率，實現(xiàn)水資源科學均衡的有效分配，使有限的水資源得到合理化利用，促進社會經(jīng)濟和生態(tài)環(huán)境可持續(xù)發(fā)展［4］。

洮河是黃河上游的第一大一級支流，其干流河道長673 km，流域面積25 527 km2，自產(chǎn)水資源總量47.2 億m3。引洮工程是解決甘肅省中部地區(qū)干旱缺水主要工程［5］，引洮一期工程自2015年建成通水，給受水區(qū)提供了一定的水資源，改善了受水區(qū)缺水狀況，促進了受水區(qū)的發(fā)展和人民生活水平的提高［5］。引洮一期供水工程是大型跨流域調(diào)水工程，該工程對受水區(qū)的水資源時空變化和區(qū)域水文循環(huán)產(chǎn)生一定影響。在國際上，大型水利工程也是調(diào)配水資源的主要工程措施［6］。美國、前蘇聯(lián)、印度等國家都是采用人為修建引水工程的措施改善干旱地區(qū)缺水的問題［7］。但是，這種工程擾動自然狀態(tài)的水資源時空分布，對取水水源及受水區(qū)的水資源水量、水質及生態(tài)環(huán)境產(chǎn)生影響，出現(xiàn)土地荒漠化和地表、地下水質惡化等水環(huán)境問題［8，9］。

本研究以引洮一期工程受水區(qū)及周邊區(qū)域水資源的變化為研究對象，通過地統(tǒng)計學方法和相關模型對比分析，揭示引洮一期工程對受水區(qū)水資源的變化影響，為進一步深入研究引洮總體工程實施的影響提供科學依據(jù)和技術支撐。

1 研究資料和研究方法

1.1 區(qū)域概況和資料

本文以引洮一期工程受水區(qū)為研究區(qū)，西起臨洮，東至會寧，南起隴西，北至榆中（34.84°N～36.45°N，103.47°E～105.52°E），受水區(qū)范圍見圖1。受水區(qū)地勢整體由西北向東南呈遞減態(tài)勢，海拔最高處在興隆山區(qū)域3 936 m，平均海拔2 088 m。年降水量300～450 mm，年蒸發(fā)量1 400～1 700 mm，干旱指數(shù)3.0～4.0，氣溫日變化和季節(jié)性變化較大，夏季日照時間長，太陽輻射充足，地面蒸發(fā)量要遠大于降水量，導致地表極度缺水，水土流失嚴重［10］，水資源嚴重缺乏。受水區(qū)屬于山丘、溝壑和盆地相互交錯地區(qū)。植被類型和群落結構單一，主要是旱生灌木、草原及稀疏林地。人類活動主要以典型的傳統(tǒng)旱作農(nóng)業(yè)為主。

圖1 研究區(qū)位置圖Fig.1 Location map of the study area

1.2 數(shù)據(jù)與方法

1.2.1 WRF模擬潛在蒸散發(fā)和降水的方法介紹

選用WRF V4.0模式［11］在引洮一期工程受水區(qū)進行氣象水文的模擬。模式采用雙層嵌套，粗細網(wǎng)格水平方向格距分別為15 km 和5 km，格點數(shù)89×89、75×75，網(wǎng)格中心35.70°N，104.95°E。模式垂直方向設有27 層地形追隨垂直坐標，積分步長為90 s。

根據(jù)WRF 模型輸出結果對受水區(qū)潛在蒸散發(fā)和降水進行計算。受水區(qū)降水數(shù)據(jù)可由WRF 模式直接輸出，并通過Arc-Map 10.2 進行裁剪繪制。潛在蒸散發(fā)ET0由模式輸出的氣溫、風速、輻射等數(shù)據(jù)，進一步根據(jù)糧農(nóng)組織（FAO）的指導方針［12］，使用Penman-Monteith 方程［公式（1）］進行計算，將計算結果通過ArcMap 10.2 進行裁剪繪制。該方法在全球范圍內(nèi)進行了測試，可以在任何氣候和土地利用中使用，無需局部校正［13-15］。具體公式如下：

式中：ET0為潛在蒸散，mm/d；Rn為凈輻射，MJ/（m2·d）；G 為土壤熱通量，MJ/（m2·d）；T 為日平均氣溫，℃；u2為2 m 高度的風速，m/s；es為飽和水汽壓，kPa；ea為實際的水蒸氣壓力，kPa；Δ 為水汽壓曲線斜率，kPa/℃；γ為濕度計常數(shù)，kPa/℃。

1.2.2 模型輸入數(shù)據(jù)

本文采用FNL 全球數(shù)據(jù)（https：//rda.ucar.edu/）作為WRF 模型的輸入數(shù)據(jù)進行模擬。FNL 數(shù)據(jù)是美國國家環(huán)境預報中心（NCEP）發(fā)布的全球再分析數(shù)據(jù)，數(shù)據(jù)精度為1°×1°，間隔時間為6 h。該產(chǎn)品來自全球數(shù)據(jù)同化系統(tǒng)（GDAS），該系統(tǒng)持續(xù)從全球電信系統(tǒng)（GTS）和其他來源收集觀測數(shù)據(jù)，用于多種分析。參數(shù)包括地表壓力、海平面壓力、位勢高度、溫度、海面溫度、土壤溫度、u、v 風速等。研究選擇引洮一期工程通水后的2015、2017、2020年分別進行模擬。

1.2.3 徑流量以及地下水位研究方法

運用統(tǒng)計分析法，對引洮工程供水前及供水后徑流量和地下水位的實測資料進行對比，然后進行圖表的繪制，通過對圖表的分析得出引洮工程對研究區(qū)的影響。

1.2.4 徑流量以及地下水位數(shù)據(jù)

河道徑流量數(shù)據(jù)來源甘肅省水文資料、供水工程實際監(jiān)測資料和實際現(xiàn)場調(diào)查的資料，地下水位的資料來源于甘肅省水文局引洮一期工程的定點實際觀測資料。

2 結 果

2.1 蒸散發(fā)的變化

蒸散發(fā)是陸地水文循環(huán)的關鍵過程，受到氣候、降水及下墊面條件等因素的影響。通過對受水區(qū)蒸散發(fā)的時空變化進行分析（圖2），發(fā)現(xiàn)蒸散發(fā)量在山區(qū)較低，在平緩的盆地和平原區(qū)較高，在夏季較高而在冬季較低。2017年的年平均蒸散發(fā)量比2015年降低了12.5%，其中，春、夏和秋三季的蒸散發(fā)量較2015年分別降低了21.6%、12.5%和11.1%，而冬季蒸散發(fā)量則升高了6.7%。在年內(nèi)夏季蒸發(fā)強烈，冬季蒸發(fā)微弱，2015年夏季的蒸散發(fā)量比冬季高6.1 倍，而2017年夏季的蒸散發(fā)量比冬季高5.2 倍。結果表明，受水區(qū)的蒸散發(fā)量在供水后較供水前有一定幅度的降低，這是由于引洮一期工程對受水區(qū)地下水和地表水的補充，導致受水區(qū)的土壤溫度有所下降，從而降低了土壤蒸發(fā)量。受水區(qū)2020年的蒸散發(fā)量與2015年的相比，山區(qū)蒸散發(fā)量有所增加，尤其在6-7月時段的變化最為明顯。說明受水區(qū)的植被生態(tài)環(huán)境改善，植物蒸騰量增強，另一方面由于引洮一期工程供水使土壤的含水量增加而引起蒸散發(fā)量的增加。

圖2 引洮工程覆蓋區(qū)潛在蒸散發(fā)變化分布圖Fig.2 Distribution map of the potential evapotranspiration changes in the coverage area of Yintao Project

2.2 降水變化

降水是受水區(qū)水資源的主要來源。運用WRF 模型對引洮工程運行前的2015年及運行后的2017年和2020年受水區(qū)降水的時空變化進行分析（圖3），結果顯示，2017年降水量比2015年增加了75.36%，春季增加了10.03%，而秋季則有所降低，在冬季基本沒有變化。2015年的夏季降水量比冬季多了10.8倍，最大值出現(xiàn)在9月，為77.9 mm，而2017年的夏季降水量比冬季多17.2 倍，最大值則出現(xiàn)在8月，為204.96 mm。2020年的降水量較2015年整體減少；在空間上，受水區(qū)的降水主要集中在南部地區(qū)，其他地區(qū)降水相對較少。受水區(qū)的降水在空間和時間上具有變異性，說明降水的變化主要受氣候的影響，而受引洮工程的影響相對較弱。

圖3 受水區(qū)近3年的月降水變化分布特征圖Fig.3 Distribution characteristics of monthly precipitation change in water receiving area in recent 3 years

2.3 徑流的變化

河道徑流量是地表水資源量的重要指標，主要來源于夏季的降水。受水區(qū)內(nèi)2018-2020年第一季度（1-3月份）的徑流量變化如表1所示。

表1 2018-2020年研究區(qū)域內(nèi)1-3月徑流量變化表Tab.1 Table of river discharge changes in the study area from January to March，2018-2020

在1-3月份，除靖遠站2020 河道徑流量較2019年減少（21萬m3），其他各站1-3月河道徑流量均逐年增加，其中堯甸站2020年1-3月徑流量較2018年、2019年分別增加323.9 萬、195萬m3；大羊營站較2018年、2019年分別增加597 萬、294 萬m3；郭城驛站則是分別增加381萬、174 萬m3。另外，臨洮站和西河站2018年1-3月均無河道徑流量，而2019年則分別增加了207萬和116 萬m3的河道徑流量，這表明自2019年開始運行后，該區(qū)段內(nèi)河道的流量逐漸增加，實現(xiàn)從無到有的情況，且兩站徑流量2021年較2019年同期分別增加174萬和221 萬m3，由此可見引洮一期工程的供水，使各河道的徑流量有了明顯的上升。

通過對2018年和2021年4-6月徑流量的分析表明（表2），2021年洮河流域徑流量較2018年增加了196.5 萬m3；臨洮（三）站增加了290.4 萬m3。祖厲河流域的西河（二）站2018年的徑流量為0，2021年較2018年增加了232.1 萬m3，徑流量實現(xiàn)從無到有的狀況。另外，祖厲河流域大羊營（二）站、郭城驛站和靖遠站2021年河道的徑流量較2018年分別增加了154.6萬、262.2萬和486 萬m3。由此可見，由于引洮一期工程的投入運行，區(qū)域內(nèi)各河道的徑流均呈現(xiàn)出大幅度的增加趨勢，水資源明顯增加。

表2 2018和2021年研究區(qū)域內(nèi)河道徑流量（4-6月）變化表 萬m3Tab.2 Table of river discharge（April-June）changes in the study area in 2018 and 2021

2.4 地下水的影響

引洮供水一期工程運行對受水區(qū)的水資源影響不僅僅是對地表水資源量的影響，同時對受水區(qū)的地下水也產(chǎn)生了強烈影響。內(nèi)官-香泉盆地是引洮一期供水工程的主要灌區(qū)，通過分析該盆地供水前后供水井的水位變化可以初步研究引洮一期工程供水對地下水資源的影響。根據(jù)甘肅省定西水文水資源勘測局的監(jiān)測資料，對黃家坡、清溪、先鋒3 個研究樣地的觀測井水位2018-2021年4-6月份的平均水位對比分析結果顯示（表3），內(nèi)官盆地的黃家坡、清溪、先鋒觀測井地下水位在2018-2021年時段內(nèi)整體呈現(xiàn)上升的變化趨勢，其中黃家坡觀測井2021年地下水位分別比2018年、2019年和2020年同期的地下水位上升了2.76、2.50和1.71 m；先鋒站2021年的觀測井水位比2018年、2019年和2020年同期的地下水位分別上升了6.54、1.08 和0.89 m；清溪站2021年監(jiān)測井水位比2018年、2019年和2020年同期的地下水位分別上升了1.93、1.14 和2.65 m。段耀峰［16］對香泉盆地、臥龍川盆地、內(nèi)官盆地3 個盆地在2001-2010年年均水位變化的研究結果顯示，3 個盆地的地下水位下降速率分別為0.66～0.99、0.12～0.92、0.21～1.01 m/a，說明研究區(qū)在引洮一期工程運行前地下水位呈現(xiàn)快速的下降趨勢，地下水的排泄量嚴重大于補給量，呈現(xiàn)負均衡狀態(tài)。而由上述分析可知引洮一期供水工程運行以后，受水區(qū)的地下水得到明顯補給，地下水位回升明顯，這表明引洮一期工程增強了對受水區(qū)地下水的補給，區(qū)域水資源矛盾得到有效緩解。

表3 內(nèi)官盆地各監(jiān)測井歷年同期地下水月平均水位分析表Tab.3 Analysis table of monthly average groundwater-level of monitoring Wells in Neiguan Basin

3 討 論

在全球變化背景下，干旱半干旱地區(qū)水資源時空差異越來越突出，引水工程對受水區(qū)水資源的影響主要表現(xiàn)在降水、蒸發(fā)、地表和地下水資源的變化上［17］。引洮一期工程的通水運行大大減輕了受水區(qū)的用水壓力，也促進了當?shù)氐纳鐣?jīng)濟發(fā)展和生態(tài)環(huán)境的逐步改善。本研究中降水和蒸散發(fā)的變化具有明顯的不確定性，2015-2020年的變化中蒸散發(fā)量表現(xiàn)出先減小后增加的變化，而降水表現(xiàn)出先增加后減少的變化過程，而且在季節(jié)上的變化又有所不同（圖2 和圖3），這說明區(qū)域氣候和水文循環(huán)過程主要受區(qū)域氣候和大氣環(huán)流的作用，受引洮工程的作用較為微弱。

自2015年引洮一期工程供水以來，洮河和祖厲河徑流量都明顯增加，部分斷流的河段恢復流水，受水區(qū)的天然降水非常稀少，不足以補給河道徑流，兩個河道徑流增加都是由于引洮一期工程的引水改變了河道徑流量。證明引水工程是調(diào)控天然水資源時空不均衡的有效方法，也是實現(xiàn)生態(tài)環(huán)境恢復治理的主要技術手段。受水區(qū)的降水補給地下水非常弱，降水主要是對淺層包氣帶水的補給，深層地下水補給非常少。通過對內(nèi)官-香泉盆地地下水位變化的調(diào)查分析發(fā)現(xiàn)，在引洮一期工程供水前，內(nèi)官-香泉盆地的地下水均衡均呈負均衡狀態(tài)，地下水排泄量均大于補給量。2019年內(nèi)官-香泉盆地地下水資源量均大于地下水資源開發(fā)利用量，說明引洮一期工程運行以來，受水區(qū)地下水資源量在逐漸增加，使受水區(qū)的地下水資源量總體上得到了補給。因此，可得出引洮一期工程運行明顯增加了受水區(qū)地下水的補給量，地下水資源得到恢復和增加，將對受水區(qū)的干旱生態(tài)環(huán)境治理恢復產(chǎn)生積極影響。

4 結 論

基于以上的研究可以得到以下結論：

（1）引洮一期工程受水區(qū)蒸散發(fā)2017年相比2015年為減少，這說明了供水對地表水和地下水的補給使土壤溫度降低從而減少了蒸發(fā)量，而2020年相比2015年為增加，這說明了受水區(qū)的植被生態(tài)環(huán)境改善，植物蒸騰量增強以及引洮一期工程供水使土壤的含水量增加而致使蒸散發(fā)量增強。

（2）在引洮一期供水后，受水地區(qū)的河道徑流明顯增加，典型洮河和祖厲河流域的徑流在全時段都有增加，其中洮河流域徑流量年均增加量為161.95～190.5萬m3；祖厲河流域徑流量年均增加量為80.5～248.5 萬m3。引洮一期工程受水區(qū)的地表徑流有著明顯的改善，實現(xiàn)了區(qū)段內(nèi)河道的水資源恢復過程和生態(tài)環(huán)境逐漸恢復。

（3）引洮一期供水受水區(qū)的地下補給非常明顯，全流域的地下水超采現(xiàn)象得到了遏制，其中黃家坡觀測井地下水位增幅為0.26～1.71 m/a；清溪站觀測井地下水位增幅為0.79～2.65 m/a；先鋒站觀測井地下水位增幅為0.79～5.46 m/a。