余昭里 穆秀麗
(江西省九江市水利局 九江 332000 江西省九江市水利電力規劃設計院 九江 332000)
SWAT模型國內徑流模擬研究進展
余昭里 穆秀麗
(江西省九江市水利局 九江 332000 江西省九江市水利電力規劃設計院 九江 332000)
隨著SWAT(Soil and Water Assessment Tool)模型的發展與更新,其功能不斷完善,國內應用該模型進行水文模擬的研究也日益增多。在簡述SWAT模型徑流模擬基本原理的基礎上,通過探究國內徑流模擬研究的相關文獻,得出近幾年研究方向,主要包括徑流模擬的適用性研究、模擬結果影響因素研究、模型參數敏感性分析與率定研究等三個方面。最終總結了該模型在徑流模擬中存在的問題和未來的研究趨勢。
SWAT模型 徑流模擬 影響因素 研究進展
1.1產流模擬的原理
SWAT模擬的流域水文過程分為水循環的陸面部分(產流和坡面匯流部分)和水循環的水面部分(河道匯流部分)。前者控制每個子流域內主河道的水、泥沙、營養物質和化學物質等的輸入量;后者決定水、泥沙等物質從河網向流域出口的輸移運動。
模型開始運算時,首先根據地形將流域劃分為若干子流域,然后根據土壤、土地利用及地形特征進一步將子流域細分為若干水文響應單元(HRU)。HRU是SWAT模型運算的最基本單位,通過計算各個HRU內的徑流量,模型提供徑流曲線數法(SCS法)Green-Ampt模型用于計算地表徑流,最后采用變動存儲系數模型或馬斯京根法演算得到流域總徑流量。這些計算基于降水、地表徑流、蒸散、壤中流、滲透、地下水回流和河道運移損失等水平衡過程。
1.2產流模擬所需輸入數據
模擬過程需要輸入的數據總體分為圖數據和表數據兩大類,其中圖數據包括DEM數據、土地利用/土地覆蓋圖和土壤類型圖。SWAT通常采用將整個流域離散的方式,首先利用DEM數據將研究流域細分為若干子流域,然后根據不同的土壤類型和土地利用圖將子流域進一步劃分為水文響應單元(Hydrological Response Unit HRU)。表數據主要包括氣象數據(包括最高最低氣溫、日降雨量、相對濕度、太陽輻射、風速等)和水文數據表格(日流量、月流量等),在產流模擬過程中可不用考慮施肥時間和數量等信息。
徑流模擬在國內已取得了一定的研究成果,總結近幾年相關文獻可知,國內關于SWAT模型徑流模型方面的研究主要分為三類:(1)不同流域的徑流模擬的適用性研究;(2)不同影響因素對SWAT模型徑流模擬結果的影響研究,主要包括如下幾個方面:不同DEM分辨率、不同土地利用/覆被變化、不同子流域劃分方法和氣候變化的水文響應研究;(3)模型參數的率定與驗證研究。
2.1不同流域SWAT模型徑流模擬的適用性研究
SWAT的徑流模擬已應用于不同國家、不同流域,而徑流模擬結果的精度和可用性也不盡相同,但總體趨勢表明,SWAT模型在徑流模擬中適用性普遍滿足要求;國內研究進展主要集中于中小型流域,而模型的開發者開發時主要是基于大、中型流域。目前,國內王中根等(2008)在大流域海河流域構建SWAT模型,為SWAT在大流域應用提供范例;而中型流域研究較多,何磊磊等以我國淮河上游各支流組成的流域為研究對象,發現在半濕潤地區SWAT模型能很好地模擬出徑流的年內、年際變化;尹竹峰等構建了渾河流域分布式水文模型,模擬結果較好,表明SWAT模型在整個渾河流域內的徑流模擬中具有良好的適用性。已有作者對環江和漢江流域進行了SWAT模擬徑流的研究,結果表明無論是驗證期還是率定期R2和ENS的值均高于評價標準(ENS>0.5,R2>0.6),評價精度較高。可見我國的大中型流域SWAT模型徑流模擬的適用性較強。而對于我國一些偏小的流域,比如三峽地區張家沖流域、昆明松華壩流域等的徑流模擬研究結果表明,在較小流域的徑流模擬中,SWAT模型也發揮著良好的水文評價效果。
2.2不同影響因素對SWAT模型徑流模擬結果的影響研究
不同影響因素對徑流模擬產生的影響結論不一,影響程度也有所不同。在國內,關于影響因素的分析研究主要有如下幾個方面:
(1)不同DEM分辨率和不同子流域劃分下的水文響應研究:王艷君等以秦淮河流域為研究對象,結果表明,產流量隨著子流域劃分數的增加幾乎沒有變化,最大相對誤差不超過5%;DEM分辨率的變化對徑流的模擬影響也不大,尤其在DEM格網單元100m時,各種格網單元模擬結果的相對偏差最大約為1%。吳軍等認為不同分辨率DEM流域地形分析計算機自動提取得到的最長河道相差較小,但河道總長、坡度等
相差較大,進而影響了分布式水文模型徑流模擬的效率。邱臨靜等對黃土丘陵溝壑區杏子河流域進行了研究,發現當DEM柵格分辨率在20~150m之間時,SWAT能有效地模擬年河川徑流、地表徑流、基流及產沙量,各水文要素模擬結果的R2和Ens分別在0.93和0.51以上,而當柵格分辨率大于150m時,各水文要素的模擬效果存在差異;而子流域劃分對流域產流模擬影響較小。李曼曼等為了明確洱海流域最佳的子流域劃分方案,得到的結論為隨著子流域的增加,年均徑流深增加但敏感性較小,相對偏差僅0.89%。同時胡連伍等也得到類似的結論,表明不同的子流域劃分層次在閾值范圍內對徑流過程的影響較小。綜合上述研究不難看出,DEM分辨率和不同子流域劃分方法對徑流模擬的結果影響不大,所以模擬結果對該因素的敏感性不高。
(2)不同土地利用/覆被變化下的水文響應研究:陳軍鋒等應用SWAT模型模擬了長江上游梭磨河流域不同土地覆被情景下的多年降水徑流關系,得知隨著土地覆被狀況由無植被到全是有林地覆被,徑流深減小,蒸發量增加,枯季徑流深減小幅度明顯小于雨季的減小幅度,而且雨季初期徑流深減小的幅度大于雨季后期。歐春平等發現土地利用/覆被變化對SWAT模型水循環模擬結果的影響很大。李昱等以遜畢拉河雙河屯測站以上流域為研究區,定量分析土地利用變化對徑流量的影響,發現增加耕地面積會增大總徑流、直接徑流,減小基流;減小林地面積會增大直接徑流,減小基流;森林涵養水源能力大于耕地。15°以上坡度的耕地退耕還林還草可以增大基流,減小直接徑流,有利于水土保持。羅巧等研究證明增加林地和草地面積將減少徑流而耕地和建設用地的增加導致徑流的增加。李佳等對長江源土地利用/覆被情景變化對徑流影響研究也表現出類似的結果。王學對白馬河流域設定了四種不同的土地利用/覆被情景計算流域四種主要的土地利用類型對徑流量的貢獻系數,得到的有林地、灌木林地和疏林地、居民點及建設用地的貢獻系數均為正值,說明三者均有促進產流的作用,而有林地對流域產流的貢獻最大;耕地的貢獻系數為負,說明耕地具有抑制徑流產生的作用。總體來說,隨著近幾年國內關于土地利用/覆被變化對徑流影響研究的發展,更加確定該因素對徑流模擬的影響顯著,只是在不同的流域影響方式稍有不同,所以這一敏感因素對不同地區土地利用管理起著不可忽視的作用,研究意義重大。
2.3模型參數敏感性分析與率定驗證
參數的敏感性分析是SWAT模型至關重要的一步,直接影響著模型模擬的精度,隨著不同版本的更新,該模塊也在不斷地更新與改進。對此,國內外已有不少研究者進行了研究。尤其是國內對該方面的研究在近幾年增加的比較明顯。
相關文獻顯示國內對參數敏感性分析最常用的方法是LH-OAT敏感性分析法。LH-OAT方法是由Morris于1991年提出的,結合了LH抽樣法和OAT敏感度分析,同時兼備這兩種方法的優點。尹雄銳等發現部分參數敏感性存在空間變異性,分析主要原因在于氣候和下墊面的空間異質性導致了流域上下游產流模式存在差異;所以即使是同一流域的不同子流域,或干流的不同區域都有可能得到不一樣的參數敏感性分析結果。白薇等證明了在進行參數敏感性分析、參數率定時,計算耗費的時間隨著子流域個數的增加而增加;若只對部分子流域進行模擬可以有效地加快率定的速率,而且結果顯示精度較高。還有研究表明,相對于LH-OAT法,SWAT-CUP敏感性分析方法的速度較快,但其精度稍弱。
總體來說,參數的敏感性分析和率定的精度高低直接影響徑流模擬的精度,它可以減少模型模擬結果的不確定性。但無論運用何種方法進行分析,結果都會存在一定的差異。特別是在不同下墊面等外因的作用下,差異都較為明顯,因此,該方面有待進一步的研究。
盡管SWAT模型在不斷的更新與改進,但仍存在一些實質性問題有待解決與提高。其突出不足之處包括以下幾個方面:
(1)徑流模擬的基礎數據缺乏,難以收集整理。SWAT模型模擬徑流所需數據涉及氣象、水文、土壤分布、土地利用等多個方面,為了更好地滿足模擬精度要求,所涉及的數量也日益增多,建立模型的數據庫工作量也日益增大。正因如此,這些數據往往不全,由于數據庫構建不完善會導致模擬結果不理想,不同地區的資料全面程度也各有不同,許多地區資料匱乏甚至無資料,導致徑流模擬研究無法進行,這一現象在國內尤為突出。并且由于該模型是由美國研發的,很多基礎數據需要進行轉化為美國制才能使用,其轉化過程增加了模擬的難度,同時對精度也有一定的影響。加強對無資料地區進行徑流模擬研究,已成為一個重要研究趨勢。
(2)徑流模擬的不確定性。物理系統接受實際的數據,并對這些實際的輸入作出反應,然而這種反應本身就包含不確定性。由于氣象數據,包括降雨、溫度、蒸發等監測數據,還有土壤的物理屬性等數據都存在一定的時空差異性,這些輸入數據的隨機性和系統誤差是導致模型不確定性的主要因素之一。同時,通常用來與模擬結果作比較的徑流量和地下水位等監測數據,也存在很大的隨機和系統誤差。最后參數值的不確定性大,存在非最優參數值所導致的誤差;這些不確定性因素都是模擬結果精度誤差的主要來源。如何進一步提高模型精確度仍然是研究的熱點。
(3)考慮徑流模擬影響因素比較局限。主要只包括氣候變化和土地利用/覆被變化(LUCC)等方面,結合考慮水利工程比如水庫以及城市化對水文的影響還很少。大部分研究都是基于自然下墊面狀況進行模擬,而在SWAT模型中加入城市模塊研究較少,在國內幾乎為空白,鑒于這些影響因素對水文循環產生著重要影響,該方面的研究將是一個明顯的趨勢■