基于SWAT模型的內蒙古錫林河流域降水-徑流特征及不同水文年徑流模擬研究

段超宇, 張 生, 李錦榮, 張成福, 綦俊諭, 吳 用

(1.內蒙古農業大學 水利與土木建筑工程學院, 呼和浩特 010018; 2.水利部牧區水利科學研究所,呼和浩特 010020; 3.新布朗斯維克大學 林業與環境管理學院, 加拿大 弗雷德里頓 E3B 5A3)

基于SWAT模型的內蒙古錫林河流域降水-徑流特征及不同水文年徑流模擬研究

段超宇1, 張 生1, 李錦榮2, 張成福1, 綦俊諭3, 吳 用1

(1.內蒙古農業大學 水利與土木建筑工程學院, 呼和浩特 010018; 2.水利部牧區水利科學研究所,呼和浩特 010020; 3.新布朗斯維克大學 林業與環境管理學院, 加拿大 弗雷德里頓 E3B 5A3)

以我國北方寒旱區草原中典型流域——錫林河流域為研究對象，結合我國北方寒旱區草原型河流流域的氣候特點、下墊面條件及水文過程的特殊性，通過分析錫林河流域近30 a的降水—徑流變化特征，利用錫林河水文站的年徑流頻率曲線，采用一定保證率條件下的年徑流量作為劃分流域豐平枯水年的標準，對錫林河流域歷史水文年進行劃分。基于此基礎，以ARCSWAT2012為操作平臺，應用地理統計分析方法，對模型建立過程中所需的各數據庫進行基礎數據整備與參數化，分析影響模型徑流模擬精度的主要參數，進行模型的參數率定與結果精度評價，分別對錫林河流域各豐、平、枯水年進行徑流模擬，并對其結果進行了檢驗。結果表明：SWAT模型對細化流域降水—徑流分配的徑流模擬具有不同的精度，即平水年>偏豐水年>豐水年>偏枯水年>枯水年，體現出SWAT模型在中國北方寒旱區豐水年和平水年具有較好的可操作性。

徑流模擬; 水文年; SWAT; 錫林河流域

錫林郭勒草原位于內蒙古自治區中部，是我國溫帶典型草原的核心分布區和重要的草地畜牧業生產基地，同時也是我國北方草原牧區和農牧交錯帶的重要組成部分。它就像是我國北方一道綠色的生態屏障，能有效地阻止草原腹地的土壤侵蝕、沙化及來自中亞和我國西部的沙塵侵害，對于維持整個華北地區，特別是京、津地區的生態環境安全起著非常重要的作用[1]。近年來，錫林河流域在以水資源和礦產資源開發利用為標志的高強度人類社會經濟活動下，流域自然生態過程發生著明顯的變化，以天然植被為主體的生態系統，因人為改變水資源的時空格局受到嚴重影響[2]。流域內水資源分配的不均勻，導致中下游流域極度干旱，目前已出現斷流情況，上游產流區對中下游地區的社會經濟和環境變化有著舉足輕重的作用。因此，作為我國北方寒旱區草原資源科學研究中最典型的流域——錫林河流域，如何運用科學的方法進行水資源分布調查、儲量分析、優化管理及水資源的預測和保護已成為目前錫林河流域地區社會經濟環境可持續發展的大問題。

錫林河流域的水資源主要來源于降水，降水的多少對流域水資源的數量及時空分布有著重要的影響，同時也決定了流域內水資源的豐枯[3]。目前多數文獻對流域內的徑流模擬僅限于單一年份[4]或連續年內[5-7]的模擬，對細化水資源分配的不同水文年模擬提及較少。因此，本研究以國內外應用較為廣泛的SWAT(Soil and Water Assessment Tool)模型[8-10]為操作平臺，通過分析內蒙古錫林河域多年降水—徑流特征，對流域內不同的水文年份進行徑流模擬研究，從而探討SWAT模型在錫林河流域這樣特殊的自然環境條件下的適用性及不同水文年份內的徑流模擬精度，為錫林河流域水資源科學管理和可持續利用提供理論依據。

1 材料與方法

1.1 流域概況

錫林河屬內蒙古內陸河流，全長175 km，流域面積105 452 km2，地理坐標介于43°26′—44°39′N，115°32′—117°12′E之間。它發源于赤峰市克什克騰旗寶爾圖西南山頂，海拔高度1 334 m，自東向西分別流經赤峰市的克什克騰旗與錫林郭勒盟的阿巴嘎旗，在貝爾克牧場轉向西北經過錫林浩特市，最終注入查干諾爾沼澤地自然消失[11]。流域自東南向西北依次為多級玄武巖臺地、低山塔拉、山前洪積平原[12]。錫林河流域是我國北方典型的草原區域，流域內有90%的面積覆蓋著草原，植物生長期約150 d[13]。

1.2 降水-徑流特征及水文年劃分

1.2.1 降水—徑流特征 錫林河流域地處中國北方寒旱區，夏季溫暖降雨集中，冬季寒冷積雪豐富，春季融雪產生較強的徑流與夏季降水會形成明顯的雙峰型地表徑流。就流域內降水而言，錫林河流域在氣候類型上為內陸半干旱氣候向東南濕潤、半濕潤季風氣候的過渡帶[14]。據統計[15]，錫林河流域1971—2000年的年均降水量為292.6 mm，年潛在蒸發量為1 700 mm，干旱指數為2。5—9月是錫林河流域主要的降雨期，其余月份為積雪期。圖1為錫林河流域1971—2000年的多年降水量的年內分配圖。

圖1 錫林河流域降水分配

通過對錫林河水文站1971—2000年實測徑流資料統計可知，錫林河多年平均徑流總量1 958.5萬m3，年均徑流深5.08 mm。由年均徑流深僅占年降水量的1.8%可知，錫林河流域每年約有98.2%的降水用于蒸發。從圖2可以看出，錫林河年內徑流分配極不均勻，表現為明顯的雙峰型地表徑流特征。每年12月到次年2月，由于河水封凍，徑流量幾乎為零。4月份后，伴隨著氣溫的回暖，地表大量冰雪在太陽輻射的作用下開始融化，但此時由于土壤層中的溫度依舊低于地表溫度，土壤層處于凍結狀態，融雪水難以向土壤層入滲，導致大量融雪水匯入河道，形成春汛。經過計算，錫林河流域4月份融雪形成的徑流占全年總徑流的32.5%。5—9月是錫林河流域主要的汛期，地表徑流的產生全部以降雨為主。其中7月和8月集中了全年56%的降水。從圖2還可以看出，錫林河流域在7月和8月產生的徑流量相對4月份偏低，這是因為7月和8月是錫林河流域植被生長最旺盛的時期，此時的降水一部分用于植被蒸散發消耗，另一部分用于產生徑流。10月中旬后，錫林河流域進入到積雪期。

圖2 錫林河流域徑流分配

1.2.2 豐平枯水年劃分 根據GB/T50095—98水文基本術語和符號標準中對河川徑流豐平枯水年的劃分標準，本文通過利用錫林河水文站年徑流頻率曲線，采用一定保證率[16]的年徑流量作為劃分流域豐平枯水年的標準。年徑流量豐平枯水年劃分標準如表1所示。

表1 錫林河流域年徑流量水文年劃分標準

通過對錫林河水文站1971—2000年的徑流量進行水文年劃分可知，近30年內錫林河流域共出現11個豐水年(3個豐水年和8個偏豐水年)、3個平水年和16個枯水年(4個枯水年和12個偏枯水年)，分別占36.7%，10%和53.3%。同時，從時間尺度上還可以看出，流域徑流呈現出以每10年為周期的變化規律，即1971—1980年偏多，1981—1990年偏少，1991—2000年偏多。表2為錫林河流域水文年劃分結果。

表2 錫林河流域年徑流量水文年劃分結果

1.3 模型構建與數據準備

1.3.1 模型結構 SWAT(Soil and Water Assessment Tool)是由美國農業部農業研究中心(USDA-ARS)開發的流域尺度模型[17-19]。針對區域性水文評價而言，SWAT模型目前在我國各流域的水文過程中已經得到了非常廣泛的應用，但是相對于寒旱區降水時空分布不均、下墊面條件復雜、水文氣象監測站分布少等特點，SWAT模型在我國北方寒旱區的模擬仍舊存在很大的困難。本文以中國北方寒旱區典型流域——錫林河流域為研究對象，通過對其徑流進行模擬，為寒旱區水資源的利用與保護提供理論依據。

SWAT模型中對水文循環的模擬計算基于水量平衡方程，見公式(1)：

(1)

式中：SWt——土壤最終含水量(mm)；SW0——土壤前期含水量(mm)；t——時間步長(day)；Rday——第i天降雨量(mm)；Qsurf——第i天的地表徑流(mm)；Ei——第i天的蒸發量(mm)；Wseep——第i天存在于土壤剖面底層的滲透量和側流量(mm)；Qgw——第i天地下水出流量(mm)。

1.3.2 數據準備 構建SWAT模型所需要的主要基礎數據包括：數字高程模型(DEM)、土地利用類型(Landuse)、土壤空間分布(Soil)及水文氣象實測數據等。本研究在劃分研究區范圍時，考慮到近年來錫林河流域內水資源分配的不均勻性，位于錫林浩特市水庫以北的錫林河流域下游已出現極度干旱、斷流等現象，且流域內僅有錫林河水文站一個監測站，通過對錫林河流域實地考察，結合流域水文站位置及水文站以北的流域水系走向，以錫林河水文站為中心，東西各取1 000 km進行勘探人為確定流域的北面邊界。

DEM數據來源于中國科學院國際科學數據服務平臺(http:∥datamirror.csdb.cn/index.jsp)SRTM(Shuttle Radar Topography Mission)90 m分辨率原始高程數據(30 m分辨率的DEM拼接后圖像上有壞點，圖像信息發生丟失，不能繼續劃分流域)。通過應用ArcGIS 10.0對所獲得圖像進行拼接、投影變化等預處理最終生成錫林河流域DEM，見附圖6。

SWAT模型中建議土地利用類型不超過10種，如果數據類型超過10種，則需要對其進行重分類，目的是減少模型最后生成的HRU的數量。土地利用數據是內蒙古師范大學內蒙古自治區遙感與地理信息系統重點實驗室在中國科學數據服務平臺上下載TM影像并對其進行提取、合并、矢量轉化，同時運用ArcToolbox中Reclassify工具對2000年錫林河流域各土地利用類型重新合并，分類結果見附圖7，土壤類型數據來源于寒區旱區科學數據中心1∶100萬土壤數據庫，如附圖8所示。

本研究區地處牧區，地廣人稀，由內蒙古氣象站點分布可知，研究區內沒有氣象站點。針對這一現狀和特點，本文選擇與研究區相鄰的氣象站作為研究區氣象數據的輸入依據。研究區周圍的氣象站點有錫林浩特站(54102)、阿巴嘎旗站(53192)、西烏珠穆沁旗站(54012)及林西站(54115)，除錫林浩特站外其余三站都距研究區較遠。通過對各站氣象數據進行差異分析后可知，各站多年平均降水量為：林西站(389.6 mm)>西烏珠穆沁旗站(353.01 mm)>錫林浩特站(292.62 mm)>阿巴嘎旗站(256.33 mm)，四站的平均降雨量為322.89 mm。在查閱大量文獻后和咨詢當地水文站工作人員與牧民后，得出錫林河流域上游年降水量一般在300 mm左右，因此模型在進行氣象數據的輸入時，將四站氣象數據全部導入到天氣發生器，綜合考慮錫林浩特站更貼近研究區的實際降水情況，故選擇錫林浩特站(43°57′N，116°07′E)1971—2000年的日值數據作為模型模擬的依據。氣象數據來源于國家氣象信息中心氣象資料室中國地面氣候資料日值數據集，包括日平均氣壓、最高氣壓、最低氣壓、平均氣溫、最高氣溫、最低氣溫、平均相對濕度、最小相對濕度、平均風速、最大風速及風向、極大風速及風向、日照時數、降水量。水文數據采用錫林河水庫上游的錫林河水文站1971—2000年實測逐日徑流資料。

2 結果與分析

根據流域內1971—2000年的水文年份劃分結果，考慮到模型的率定期必須同時包括豐平枯不同的水文年份，將1971—1987年作為模型的率定期，1988—2000年作為模型的驗證期，其中1971—1973年為模型的預熱期。參數率定過程中，對流域所劃分的9個子流域在徑流的校核上遵循先上游后下游；先調整水量平衡再調整過程；先調整地表徑流再調整土壤水、蒸發和地下徑流的原則[20]。通過實地調查流域水文過程進行手動調參，最終確定影響不同水文過程的參數分別是：① 地表水分配參數CN2；② 植被蒸發參數IPET、CANMX、ESCO、CNOP；③ 土壤水參數SOL_AWC；④ 地下水參數ALPHA_BF、GWQMN、GW_REVAP、REVAPMN、RCHRG；⑤ 融雪參數SFTMP、SMTMP、SMFMX、SMFMN、TIMP、SNOCOVMX 。

2.1 豐水年徑流模擬

(1) 偏豐水年。由表2可知，率定期內的偏豐水年份共有4個，分別為1974年、1978年、1979年和1986年；驗證期內的偏豐水年份有3個，分別為1990年、1992年和1996年。模型對偏豐水年的徑流模擬結果顯示，率定期內的相關系數R2為0.62，納什系數Ens為0.59，驗證期內的R2為0.77，Ens為0.72，模擬結果較滿意。圖3為率定期和驗證期的偏豐水年內SWAT模型徑流擬合結果。

圖3 偏豐水年內SWAT模型徑流擬合結果

(2) 豐水年徑流模擬。同樣，對流域的豐水年份進行模擬。從模擬結果可以看出，豐水年的徑流模擬精度在率定期時(R2=0.80，Ens=0.78)較優于偏豐水期，而在驗證期時的模擬精度(R2=0.60，Ens=0.52)卻略低于偏豐水期。造成這樣的模擬結果可能是由于1998年全國范圍內降水量普遍偏大，模型對水量突增的特殊年份沒有較好的適應性。圖4為率定期和驗證期的豐水年內SWAT模型徑流擬合結果。

圖4 豐水年內SWAT模型徑流擬合結果

2.2平水年徑流模擬

由于錫林河流域在1971—2000年內平水年出現的頻率較少，僅有1972年、1991年和1997年這3個年份，且1972年又處于模型的預熱期內，因此模型在對平水年的徑流模擬評定中只對平水年做模型的驗證。驗證結果顯示，平水年的R2為0.797 9，Ens為0.65。率定期和驗證期的平水年內SWAT模型徑流擬合結果見圖5。

2.3 枯水年徑流模擬

通過對錫林河流域偏枯水年和枯水年進行模擬發現，SWAT模型不論在偏枯水年還是枯水年對徑流的模擬效果均不理想。枯水年夏季降雨偏少，蒸發量因植被的生長依舊很高，實測徑流遠遠低于模擬值。圖6和圖7分別是偏枯水年和枯水年內SWAT模型徑流擬合結果。

通過以上對錫林河流域豐平枯不同的水文年份進行徑流模擬的結果可以看出，模型在模擬流域豐水年和平水年時具有比較理想的模擬效果，在枯水年時模擬的精度不高。SWAT模型在不同水文年的模擬精度依次為：平水年>偏豐水年>豐水年>枯水年>

偏枯水年。

圖5 平水年內SWAT模型徑流擬合結果

圖6 偏枯水年內SWAT模型徑流擬合結果

圖7 枯水年內SWAT模型徑流擬合結果

3 結 論

基于SWAT模型對內蒙古錫林河流域上游不同水文年展開的徑流模擬研究表明，SWAT模型在細化流域降水—徑流模擬時具有不同的精確程度，即平水年>偏豐水年>豐水年>偏枯水年>枯水年。雖然SWAT模型從整體上看在該流域具有滿意的模擬效果，但相比南方半濕潤、濕潤區，SWAT模型在中國北方寒旱區的模擬效果上仍舊不盡如人意。造成這樣的原因可能有以下兩個方面：(1) 寒旱區降水—徑流分配不均勻，模擬時對蒸發、融雪等物理過程把握不清晰；(2) 該區域地理面積大，水文氣象監測站點分布少，數據稀缺。數據的準確性與詳細性會直接影響到模型對水文過程的模擬，對寒旱區而言今后應當加強水文氣象的監測。此外，今后還應當對SWAT模型的積雪融雪模塊加以改進，使SWAT模型在模擬我國北方寒旱區徑流時具有更高的模擬精度。

[1] 洪立國.錫林河流域草原利用現狀與保護對策研究[D].呼和浩特:內蒙古大學,2009.

[2] 王忠華.錫林郭勒盟草原退化沙化原因分析及其保護利用對策[J].內蒙古科技與經濟,2002,2(2):49-51.

[3] 范廣洲,呂世華,程國棟.華北地區夏季降水資源特征分析及其對氣候變化的響應[J].高原氣象,2002,21(1):45-50.

[4] 袁軍營,蘇保林,李卉,等.基于SWAT模型的柴河水庫流域徑流模擬研究[J].北京師范大學學報:自然科學版,2010,46(3):361-365.

[5] 李慧,雷曉云,靳晟.基于SWAT模型的山區冰雪融水河流的日徑流模擬研究[J].灌溉排水學報,2010,29(3):105-108.

[6] 于磊,朱新軍.基于SWAT 的中尺度流域土地利用變化對水文響應模擬研究[J].水土保持研究,2007,14(4):53-56.

[7] 龍銀平,張耀南,趙國輝,等.SWAT模型水文過程模擬的數據不確定性分析：以青海湖布哈河流域為例[J].冰川凍土,2012,34(3):660-667.

[8] Nie W, Yuan Y, Kepner W, et al. Assessing impacts of Landuse and Landcover changes on hydrology for the upper San Pedro watershed[J]. Journal of Hydrology,2011,407(1):105-114.

[9] 張蕾,盧文喜,安永磊,等.SWAT模型在國內外非點源污染研究中的應用進展[J].生態環境學報,2009,18(6):2387-2392.

[10] 康杰偉,李碩.SWAT模型運行結構與組織研究[J].地球信息科學,2007,9(5):76-82.

[11] 李楊,賈德彬,朱仲元,等.融雪的水量平衡模型在內蒙古錫林河流域水文模擬中的應用[J].水資源與水工程學報,2003,22(1):79-86.

[12] 仝川,雍偉義,烏云娜,等.錫林河流域近15年植被空間結構變化遙感分析[J].內蒙古大學學報,2001,32(9):562-566.

[13] 陳佐忠.草原生態系統研究(第3集)[M].北京:科學出版社,1988.

[14] 季勁鈞,劉青,李銀鵬.半干旱地區地表水平衡的特征和模擬[J].地理學報,2004,59(6):964-971.

[15] 包蘇雅.錫林河流域景觀格局動態變化研究[D].呼和浩特:內蒙古師范大學,2012.

[16] 溫川,殷戰陽.格爾木河格爾木站年徑流特征及豐平枯水年劃分[J].科技信息,2011(13):447-450.

[17] 胥彥,秦耀民,李懷恩,等.SWAT模型在陜西黑河流域非點源污染模擬中的應用[J].水土保持通報,2009,29(4):114-219.

[18] 郭曉軍,崔鵬,朱興華.典型泥石流流域蔣家溝的降雨—徑流模擬[J].水土保持通報,2011,31(1):176-179.

[19] 李麗嬌,薛麗娟,張奇.基于SWAT的西苕溪流域降雨—徑流關系及水量平衡分析[J].水土保持通報,2008,28(5):81-85.

[20] 朱新軍,王中根,李建新,等.SWAT模型在漳衛河流域應用研究[J].地理科學進展, 2006,25(5):105-111.

AnalysisonRainfall-RunoffCharacteristicsandSimulationoftheDifferentHydrologicYearRunoffofXilinRiverBasininInnerMongoliaBasedonSWATModel

DUAN Chao-yu1, ZHANG Sheng1, LI Jin-rong2, ZHANG Cheng-fu1, QI Jun-yu3, WU Yong1

(1.CollegeofWaterResourcesandCivilEngineering,InnerMongoliaAgriculturalUniversity,Hohhot010018,China; 2.InstituteofWaterResourcesforPastoralArea,InnerMongolia,Hohhot010020,China; 3.FacultyofForestry&EnvironmentalManagement,UniversityofNewBrunswick,FrederictonE3B5A3,Canada)

Xilin River is a typical grassland river located in Inner Mongolia, which is a cold-arid-region in the north of China. Its basin was investigated and some necessary experiments were carried out by means of hydrological processes. The hydrological nonlinear system theory, geological statistics, remote sensing and geographic information system were used to determine the parameters that are required by SWAT. In conclusion, using SWAT to simulate the hydrological processes in the cold-arid-regions with a characteristic of snowmelt-runoff could represent different accuracies. The best fitting result between the observed and simulated data was obtained for the normal hydrologic years, the worst for the dry years, between which the accuracies were partial wet years > wet years > partial dry year.

runoff simulate; hydrologic years; SWAT model; Xilin river basin

2013-12-24

：2014-01-08

國家自然科學基金資助項目(511069007,51169011,51169017,51269016,5133902);國際合作項目(2011DFA90710)

段超宇(1988—),女,山西省大同市,碩士研究生,主要從事水環境修復與保護研究。E-mail:duanchaoyu123@126.com

張生(1964—),男,內蒙包頭市人,教授,博士,主要從事水環境修復與保護研究。E-mail:shengzhang@imau.edu.cn

P333.1

：A

：1005-3409(2014)05-0292-06