基于SWAT模型分析嘉陵江流域藍、綠水資源量的時空變化特征

楊 昊， 倪福全， 鄧 玉， 沈敬雄

(四川農業大學 水利水電學院， 四川 雅安 625014)

1 研究背景

在全球氣候變暖和人類活動加劇的影響下，水資源短缺和分配不平衡問題日益突出，已經成為制約我國經濟社會發展的重要因素[1-2]，而傳統的水資源評價已經不能滿足解決當前水安全問題的要求，為了更好地理解水資源在陸地生態環境中的作用，瑞典科學家Falkenmark[3]在評價半干旱地區水資源對農業生產過程的影響中首次提出藍水、綠水的概念。藍水、綠水概念為水資源管理的研究及決策方式提供了新思路、新方法，引起了國內外學者的廣泛關注，逐步成為水文與水資源領域的重要研究方向[4-8]。

藍水關系到人類活動用水需求和流域可用水量，綠水關系到流域農業和生態系統需求水量[9-10]，探究流域的藍、綠水時空分布特征，對流域的水資源規劃管理具有重要的價值和意義。SWAT(soil and water assessment tool)模型由于其效率高、精度高、模擬周期較長和易于操作且能直接輸出藍、綠水分量，在國內外常被用于大尺度流域水資源的研究。例如：聶超等[11]利用SWAT模型定量分析了岷沱江流域氣候和土地利用對徑流的影響；谷新晨等[12]運用CMADS(China Meteoroligical Assimilation Driving Datasets for the SWAT model)與SWAT模型對瑪納斯河流域的水文過程進行了模擬；Nasiri等[13]使用SWAT模型對伊朗薩馬爾坎流域進行了水量分析；吳小宏等[14]結合SWAT模型對涇河長系列水沙變化規律進行了研究；Liu等[15]對老哈河流域土地覆蓋和土地利用對藍、綠水的影響進行了研究；馮暢等[16]總結了國內外綠水資源的研究進展和發展趨勢以及局限性因素，分析了氣候變化和人類活動對綠水資源的影響；趙安周等[17]應用SWAT模型討論了渭河流域氣候變化和人類活動對藍、綠水資源的影響；Schuol等[18]利用SWAT模型對非洲可利用藍、綠水的時空分布進行了評價；呂樂婷等[19]利用SWAT模型結合SUFI-2算法評估了細河流域藍、綠水的變化；李文婷等[20]基于SWAT模型對黃河源區藍、綠水的時空分布特征進行了研究。

嘉陵江流域位于成渝雙城經濟圈，是國家實施的長江經濟帶和一帶一路倡議的重要組成部分，是實現我國新發展格局的重要組成區域。目前，對嘉陵江流域的藍、綠水研究尚屬空白，因此本研究建立嘉陵江流域SWAT模型，對嘉陵江流域水資源量進行模擬研究，進一步對嘉陵江流域藍水資源量和綠水資源量進行定量評估，探討研究區藍水資源量和綠水資源量的時空分布特征，以期為嘉陵江流域藍水、綠水的管理提供理論依據，促進成渝雙城經濟圈綠色健康發展。

2 資料來源與研究方法

2.1 研究區概況

嘉陵江發源于秦嶺北麓的陜西省鳳縣代王山，地處北緯29°20′～34°30′、東經102°45′～108°55′，主要支流有涪江、西漢水、白龍江、八渡河和渠江等，全長1 345 km，干流流域面積為3.92×104km2，總流域面積為16×104km2[21]，是長江流域中流域面積最大的支流，也是長江流域中流量第二大的支流。截止到2018年底，嘉陵江流域人口約為4 720萬人，GDP總量約為16 807×108元。流域主要農作物為冬小麥、玉米、水稻、土豆和油菜等[22-23]，土地利用類型主要有林地、草地和耕地，土壤類型主要為黃壤土、紫色土和水稻土等。嘉陵江流域概況見圖1。

圖1 嘉陵江流域概況

2.2 數據來源

本文選擇嘉陵江流域范圍內、外共20個氣象站點，各氣象站點的數據源于國家氣象局整編的1981-2020年的逐日氣象數據，其中包括降水、最高氣溫、最低氣溫、相對濕度和風速等。水文數據采用的是北碚水文站1987-2018年的逐月徑流數據。此外還采用了數字高程(digital elevation model, DEM)、土地利用和土壤數據。采用的數據類型及其來源詳見表1。

2.3 SWAT模型

SWAT模型是美國農業部開發的流域尺度的分布式水文模型[24]，已廣泛應用于流域管理和水資源決策等領域。該模型可以直接輸出藍水和綠水流量的各分項，因而被認為是估算藍水和綠水流量的有效方法[25]。

表1 研究采用的數據類型及其來源

SWAT模型的適用性一般采用擬合優度確定系數R2、納什效率系數NSE及偏差百分比PBIAS作為評價標準進行評判，R2和NSE反映了模擬值與實測值的擬合度，一般認為當R2>0.6、NSE>0.5時模型的擬合精度合格，R2越接近于1則模擬效果越好[26]；PBIAS是用來衡量模擬值與觀測值的平均趨勢，PBIAS越接近于0則模擬性能越好，其最優值為0，若為負值表明模擬值被高估，若為正值表明模擬值被低估[27]。

2.4 藍、綠水量計算方法

SWAT模型對流域的水資源評估有較高的準確度[9]，由模型的輸出結果可以直接得出研究流域的產水量、綠水流量、綠水儲量、深層含水層補給量等各項分量。本研究采用面積加權法計算嘉陵江各子流域與全流域藍水資源量與綠水資源量[28]。綠水系數GWC(green water coefficient)表示綠水資源量在水資源量中的占比，能反映出研究流域的水資源分配情況。藍、綠水資源量與綠水系數計算公式如下：

BW=WYLD+DA_RCHG

(1)

GW=ET+SW

(2)

(3)

式中：BW、GW分別為藍水、綠水資源量，mm；WYLD為產水量，mm；DA_RCHG為深層含水層補給量，mm；ET為綠水流量，mm；SW為綠水儲量，mm；GWC為綠水系數。

2.5 Mann-Kendall趨勢檢驗和突變分析

Mann-Kendall(M-K)檢驗是一種非參數的突變檢驗，廣泛應用于水文和氣候長時間序列的檢驗與分析，其檢驗結果受少數異常值影響小，且能夠明確突變的區域和突變開始的時間[29]。本研究通過Mann-Kendall檢驗分析嘉陵江流域藍、綠水的顯著性及變化趨勢。

3 結果與分析

3.1 SWAT模型的率定與驗證

選取流域出口北碚水文站的徑流量實測值對SWAT模型進行校準和驗證，以1987-2000年作為模型模擬參數的率定期，以2001-2018年為模型模擬的驗證期，率定期和驗證期擬合結果分別見圖2、3。由圖2、3可知，在模型率定期和驗證期，北碚水文站月徑流量的模擬值與實測值較為吻合，均呈現出周期性的波動過程，并與年內降水量變化趨勢一致。率定期和驗證期月徑流量模擬值與實測值擬合的決定系數R2分別為0.95、0.93，NSE分別為0.94和0.92，率定期PBIAS為4.9，驗證期PBIAS為-2.7，說明SWAT模型在率定期低估了某些洪水期的峰值流量，而在驗證期高估了某些洪水期的峰值流量。R2、NSE與PBIAS的值均表明SWAT模型在本研究區的適用性較好。

圖2 1987-2000年模型率定期北碚站月徑流量模擬值與實測值擬合結果

圖3 2001-2018年模型驗證期北碚站月徑流量模擬值與實測值擬合結果

3.2 嘉陵江流域藍、綠水資源量時間變化

SWAT模型完成率定和驗證后，根據模型的輸出結果對嘉陵江流域的水資源量進行進一步的統計計算，結果見圖4。由圖4可見，1981-2020年間流域降水量(P)在662～1 198 mm變化，藍水資源量在166～644 mm變化，綠水資源量在497～597 mm變化。研究流域多年平均降水量為931 mm，其中藍水資源量為386 mm，占水資源總量的41.5%；綠水資源量為545 mm，占水資源總量的58.5%。在整個研究期內，流域的降水量呈一定的周期性波動，藍水資源量與降水量的波動規律相似；而綠水資源量也呈一定的波動變化，但波動的幅度較小。線性趨勢結果顯示，1981-2020年嘉陵江流域降水量與綠水資源量有所增加，藍水資源量有所減少，但均未通過M-K置信度為95%的顯著性檢驗，表明增加和減少的趨勢均不明顯。

圖4 1981-2020年研究區年降水量、藍水資源量和綠水資源量的變化過程

3.3 嘉陵江流域藍、綠水資源量空間分布

將嘉陵江流域劃分為21個子流域，以各子流域為統計單位，計算并繪制出整個流域1981-2020年共40年的平均降水量及平均藍、綠水資源量的空間分布圖，如圖5所示(圖中數字1～21為各子流域的編號)。

由圖5(a)可見，嘉陵江流域上游支流西漢水所在的1號子流域、白水江所在的3號子流域以及白龍江所在的4號子流域的降水量最少，其中西漢水子流域的降水量最少，為480 mm，流域中下游西部地區的涪江支流所在的7、13、14、15、16號子流域降水量并不豐富，嘉陵江干流中下游所在的17、18號子流域降水量為平均水平，流域中下游東部區域的渠江支流所在的9、10、11、12號子流域降水量最大，其中10號子流域降水量最大，為1 288 mm。由此可見，嘉陵江流域降水量空間分布極不均勻，這與周珍[30]研究得出的嘉陵江流域降水量時空分布特征基本一致。

分析圖5(b)可知，1981-2020年研究區平均藍水資源量為386 mm，其分布極不均勻。嘉陵江流域上游支流西漢水所在的1號子流域、白水江所在的3號子流域以及白龍江所在的4號子流域的藍水資源量最為貧瘠，渠江上游所在的10、11、12號子流域藍水資源量最豐富，西漢水所在的1號子流域藍水資源量僅有34 mm，而渠江上游10號子流域的藍水資源量為724 mm，前者不到后者的5%。白水江、白龍江、西漢水、涪江上游所在的1、3、4、8號子流域的面積之和占研究區總面積的26.9%，但藍水資源量僅占整個研究區的3.1%，渠江上游所在的10、11、12號子流域的面積之和占研究區總面積的14.11%，藍水資源量卻占整個研究區的23.14%，其余子流域藍水資源量均在400～500 mm之間。

分析圖5(c)可知，1981-2020年研究區平均綠水資源量為545 mm，其中白水江所在的3號子流域綠水資源量最少，為434 mm，渠江上游所在的10號子流域綠水資源量最豐富，為612 mm，兩區域差值為178 mm?？傮w而言，研究區綠水資源量分布不均勻，上游地區和下游地區西部的綠水資源量較為貧瘠，中下游東部的渠江上游所在的9、10、11、12號子流域綠水資源量較為豐富。

3.4 典型水文年嘉陵江流域藍、綠水資源量變化特征

3.4.1 典型水文年份的確定 本文采用《水文情報預報規范》(GB/T 22482—2008)中的年徑流距平百分率(Pa)作為劃分研究流域徑流豐、平、枯水年的指標[31]，具體劃分標準見表2。

1981-2020年嘉陵江流域共計出現豐水年6次，其中1983年雨量最多；平水年10次，2000年最為典型；枯水年12次，1997年最為干旱；偏豐水年與偏枯水年分別出現8次和3次?；谏鲜鰯祿?，本文選取1983年作為典型特豐水年，2000年作為典型平水年，1997年作為典型特枯水年，研究流域在各典型水文年的藍水資源量和綠水資源量變化情況。

3.4.2 典型水文年藍、綠水資源量的變化 研究流域綠水資源量在特豐水年為554.60 mm、平水年為545.53 mm、特枯水年為496.66 mm，各典型水文年的綠水資源量變化相對不大；藍水資源量在特豐水年為643.70 mm，平水年為456.71 mm，特枯水年為165.61 mm，特豐水年藍水資源量是特枯水年的3.89倍；典型水文年間的深層含水層補給量和產水量相對變化最大，而綠水儲量和綠水流量相對變化較小。在特枯水年綠水系數最大約為75%，遠大于特豐水年綠水系數的46%，表明枯水年的綠水資源量對維持生態平衡尤為重要。嘉陵江流域典型藍、綠水資源量的統計結果見表3。

表2 研究區豐、平、枯水年劃分標準

3.4.3 典型水文年藍、綠水資源量的空間分布 基于SWAT模型輸出的數據，結合ArcGIS繪制研究區各典型水文年降水量、藍水資源量和綠水資源量的空間分布，分別見圖6～8。

由圖6可見， 研究區在特豐水年的降水量分布極不均勻，在特枯水年的降水量分布較為均勻。各典型水文年降水量在涪江上游所在的8號子流域變化最小，特豐水年與特枯水年的降水量在渠江上游所在的10、9號子流域變化最大，差值分別為1 210、1 134 mm，差異顯著；西漢水所在的1號子流域、白水江所在的3號子流域、涪江中下游所在的15、16號子流域和嘉陵江流域出口所在的20、21子流域在特豐水年與特枯水年的降水變化量均為200～300 mm；嘉陵江干流上游所在的2、5、6號子流域、涪江上游所在的7號子流域、涪江中游所在的13、14號子流域和渠江中上游所在的11、12號子流域在特豐水年與特枯水年的降水變化量均為500～600 mm。渠江中下游所在的18號子流域特豐水年降水量比平水年降水量少97 mm，其他20個子流域特豐水年降水量均大于平水年降水量；渠江上游所在的10號子流域在特豐水年與平水年的降水變化量最大，為573 mm。各子流域平水年與特枯水年降水量變化范圍為75～807 mm，其中，嘉陵江干流上游所在的2號子流域變化量最小，為75 mm，渠江上游所在的9號子流域變化量最大，為807 mm；除此之外，渠江上游所在的10號子流域、嘉陵江干流中下游所在的17號子流域的降水變化量也較大，差值分別為637、682 mm。

圖7 研究區典型水文年藍水資源量空間分布

圖8 研究區典型水文年綠水資源量空間分布

由圖7可見，特豐水年研究區藍水資源量分布極不均勻，總體呈現出上游少、中游西部地區較平均、中下游東部地區最豐富的特點；特豐水年白水江所在的3號子流域、白龍江所在的4號子流域、西漢水所在的1號子流域的藍水資源量均小于100 mm，嘉陵江流域中下游所在的渠江支流上游9、10、11、12號子流域藍水資源量遠大于1 000 mm。特枯水年研究區的藍水資源量較均勻，研究區上游各子流域均小于100 mm，嘉陵江流域中下游地區除渠江上游所在的10號子流域外其余地區藍水資源量均在101～400 mm之間。平水年與特枯水年相比，除中游東部地區藍水資源量明顯增加外，其他區域藍水資源量與特枯水年空間分布基本一致。特豐水年渠江中上游所在的9、10、11、12號子流域藍水資源量占全流域的39%，而白水江、白龍江和西漢水所在的1、3、4號子流域藍水資源量僅占0.79%；同時，特豐水年與特枯水年相比，渠江上游所在的9、10號子流域藍水資源變化量超過了1 000 mm，而白水江、白龍江、西漢水和涪江上游所在的1、3、4、8號子流域藍水資源量變化不明顯；平水年與特枯水年相比，渠江中上游9、10、11、12子流域、嘉陵江干流中下游所在的17號子流域以及渠江中下游所在的18號子流域藍水資源量變化幅度均超過400 mm，其他子流域變化均小于200 mm。

由圖8可見，特豐水年除西漢水所在的1號子流域和白水江所在的3號子流域之外，其他區域綠水資源量分布較為均勻；特枯水年西漢水、白水江、白龍江所在的1、3、4號子流域綠水資源量較少，渠江中上游所在的11、12號子流域較豐富，其他區域分布較為均勻；平水年全流域綠水資源量分布較為均勻。整體來看，綠水資源量在各典型水文年的空間分布變化不大，特豐水年與平水年相比，僅1、3、4、5、19號子流域綠水資源量變化幅度超過100 mm，7、10～14號子流域在特豐水年的綠水資源量小于平水年；平水年與特枯水年相比，1、3、4、8號子流域綠水資源量變化幅度超過100 mm，7、12～14號子流域在特枯水年的綠水資源量大于平水年； 4號子流域在各典型水文年的綠水資源量變化幅度最大，15號子流域變化幅度最小。

綜上所述，嘉陵江流域降水量、藍水資源量、綠水資源量在豐、平、枯水年時空分布趨勢特征基本一致，且藍水資源量和綠水資源量的空間分布與降水量的空間分布也基本一致。流域在各典型水文年的降水量和藍水資源量的時空分布不均勻，流域東部地區在豐、平、枯水年的降水量和藍水資源量均大于其他地區。嘉陵江流域的豐水年降水量是枯水年降水量的1.81倍，特豐水年藍水資源量是特枯水年藍水資源量的3.89倍，綠水資源量在豐、平、枯水年變化幅度不大。在特枯水年，整個流域出現藍水資源量短缺的現象，除流域下游東部地區外，其他大部分地區的藍水資源量均不足200 mm；在特豐水年，藍水資源量最豐富，除流域上游部分地區外，大部分地區的藍水資源量均在400～870 mm之間。除流域上游部分地區外，大部分地區的綠水資源量在豐、平、枯水年均為460～650 mm。

4 討 論

研究流域多年平均降水量為931 mm，是流域水資源的主要組成部分，其中藍水資源量為386 mm，占水資源總量的41.5%；綠水資源量為545 mm，占水資源總量的58.5%。研究流域各典型水文年的藍水資源量和綠水資源量分布與降水量分布基本一致。綠水系數在不同典型水文年差異明顯，特枯水年為74.99%，特豐水年為46.28%，說明干旱年份蒸散發消耗的水資源量大于特豐水年。特豐水年降水量是特枯水年降水量的1.81倍，特豐水年藍水資源量是特枯水年藍水資源量的3.89倍，而不同典型水文年的平均綠水資源量基本無變化趨勢，表明特豐水年的降水量大部分都轉為藍水資源量。

嘉陵江流域大部分區域屬亞熱帶濕潤季風氣候，夏季高溫，雨熱同期，冬、夏季節氣溫和降水量差異較大。流域中下游盆地區域海拔較低，上游山區海拔較高，整體地形起伏大，上、下游高程差異導致嘉陵江流域上下游溫差較大，有明顯的梯度特征。同時，秦嶺山脈影響了部分來往嘉陵江流域的水汽輸送，位于秦嶺山脈的大巴山腳下的巴中、達州、萬源地區成為降水的集中地，流域的氣溫、氣流和降水量受到高程的影響，這與楊軍等[32]、王瓊[33]、范利杰等[34]、陳桂亞等[35]對長江流域氣候與降水量的研究結果基本一致。近年來，隨著全球變暖和溫室效應的不斷加劇，嘉陵江流域氣溫呈現出明顯的上升趨勢[36-37]，流域中上游山區多為林地和草地，下游平原多為耕地，氣溫升高導致了中上游用于植物蒸騰的綠水增加，各類型土地覆被對藍、綠水利用率的不同也導致了流域內藍水資源量和綠水資源量的空間分布不均。嘉陵江流域水力資源較豐富，目前已建、在建及規劃建設的梯級電站共有24座[38]，除區域內降水量的影響外，這些水利樞紐蓄水截流導致相應工程區的水資源量逐漸增加，越靠近水利樞紐的子流域的藍水和綠水資源量越大，這也造成了藍水資源量和綠水資源量分布的不均勻。

隨著國家“一帶一路”倡議、“長江經濟帶”戰略規劃的不斷推進以及國家新發展格局的構建，區域經濟的增長更是整個國民經濟發展的重要內容。嘉陵江流域經濟是川渝經濟帶的重要組成部分，但區域內工業發展相對滯后，總體經濟實力薄弱。要保持流域內經濟的可持續發展，須保障水資源安全，川渝地區應在水資源戰略規劃上協同好各區域間的水資源問題，通過分析嘉陵江流域水資源現狀，在人類活動日益加劇的情況下，厘清變化環境下不同區域藍水資源量和綠水資源量的分布狀況，以及流域水利樞紐的修建對藍水資源量和綠水資源量的影響，從而對嘉陵江流域的水資源進行科學地規劃與管理。

5 結 論

(1)本研究建立了嘉陵江流域SWAT模型，并對模型進行了率定與驗證，率定期R2為0.95，ENS為0.94，PBIAS為4.9；驗證期R2為0.93，ENS為0.92，PBIAS為-2.7，模型模擬效果較好。

(2)嘉陵江流域多年平均降水量為931 mm，其中藍水資源量為386 mm，占水資源總量的41.5%；綠水資源量為545 mm，占水資源總量的58.5%。1981-2020年流域內降水量和綠水資源量呈增加趨勢，藍水資源量呈減少趨勢，但變化趨勢均不顯著。典型水文年之間的藍水資源量變化顯著，而綠水資源量變化趨勢不顯著，特枯水年、平水年、特豐水年的藍水資源量分別為165.61、456.71和643.70 mm，綠水系數在典型水文年的變化差異較大，特豐水年為46.28%，平水年為54.43%，特枯水年為74.99%。

(3)嘉陵江流域藍水資源量和綠水資源量呈現出與降水量基本相同的分布規律，呈現從上游到下游先減小后增大、西少東多的空間格局。