基于SWAT 模型的湘江流域上游徑流量模擬研究

易 勇，曾 鳴

（湖南省長沙生態環境監測中心，湖南 長沙 410001）

湘江是湖南省境內最大的河流，被稱為湖南的母親河。近些年全球氣候變暖的背景下，湘江流域出現洪澇的極端天氣災害頻繁。本文只選取氣候因子來研究氣候變化對湘江流域水資源的分布影響。研究湘江流域上游的徑流量分布特征，對合理規劃和管理湘江流域水資源具有重要的現實意義[1～2]。此外，利用SWAT 水文模型模擬湘江流域上游徑流量的研究案例比較少。因此，本文通過收集湘江流域上游2008～2017 年土地利用和水文氣象數據，運用SWAT 模型對該流域徑流量進行模擬預測。

1 流域概況

湘江是洞庭湖水系中流域面積最廣的河流，地處東經110°30′～114°01′，北緯24°31′～29°01′之間，東部為幕阜山脈-羅霄山脈，南部與兩廣相鄰，西部為衡山山脈，北接洞庭湖入長江。湘江發源于湖南永州藍山縣，主要流經湖南省的永州、衡陽、株洲和湘潭等城市，河流在岳陽市匯合后進入洞庭湖[2～4]。湘江從源頭到入湖口全長為856 km，流域面積大約為94 660 km2。本文選擇湘江流域上游區域進行研究（見圖1）。

圖1 研究區域地理位置圖

2 數據處理及流域劃分

2.1 數據預處理

1）DEM 數據。DEM 從地理空間數據云網站下載，數據精度為90 m。

2）土地利用數據。本文利用的2015 年土地利用數據是從資源與環境數據云平臺下載，再經過重分類處理為城市用地和農業用地等6 種土地利用類型。

3）土壤數據。土壤數據從寒區旱區科學數據中心獲得，通過數據裁剪分類處理，土壤類型共有27 種類型。

4）氣象水文數據。氣象數據是從中國氣象科學數據網下載的，再經過天氣發生器處理，包括溫度和降雨數據；水文數據由湖南省水利廳提供，數據尺度為2008～2017 年日徑流量（歸陽站）。

2.2 流域空間劃分

1）流域水系提取。DEM 數據導入SWAT 模型中，根據高程、坡向和坡度識別并計算水流方向，最終生成各流域水系。

2）子流域的劃分。在SWAT 模型中的Stream Definition 模塊中設置面積閾值，模型自動識別出35 個子流域。

3）水文響應單元劃分。在子流域劃分的基礎上進行子流域參數的計算劃分出164 個水文響應單元。

3 結果與分析

3.1 模型模擬效果分析

用SWAT-CUP 模型進行徑流量的率定，先對2008～2012 年的徑流量數據進行率定，再將參數復制到模型驗證模塊，對2013～2017 年的徑流量再模擬60次進行驗證。結果顯示率定期和驗證期的評價系數值都大于0.5，說明用SWAT 模型模擬能較好反應該流域的水文信息（見圖2）。

圖2 研究區域實測與模擬徑流量對比及相關關系圖

3.2 模擬結果分析

將模擬生成的2050 年氣象數據輸入SWAT 模型進行計算模擬，得出2015～2050 年湘江流域各子流域的徑流量。湘江流域上游地表徑流量平均值將由2015 年的1 078 m3/s 上升至2050 年的1 296 m3/s，每年以6.23 m3/s 的速率增長（見圖3）。

圖3 2015 年與2050 年各子流域徑流量趨勢對比

4 結論與討論

通過搭建SWAT 模型框架識別出了35 個子流域和164 個水文單元，結合2008～2012 年實測徑流數據進行參數率定，并通過2013～2017 年實測數據進行驗證。模擬結果表明，SWAT 模型對湘江上游徑流的模擬效果較好，為湘江流域水文資源合理開發與利用提供了合理的科學依據。

模擬結果表明，湘江流域上游2015～2050 年間地表徑流量上升，地表徑流量平均值由2015 年的1 078 m3/s上升至2050 年間的1 296 m3/s，年增長速率為6.23 m3/s。