陳彥平
(山西省水文水資源勘測局,山西 太原 030001)
由于環境污染的日益加劇,工業污染物的排放量逐漸提高,導致我國氣候條件出現了溫度升高、降雨減少等現象[1- 2]。據有關研究表明,1980—2016年,全球地表平均溫度升高了0.85℃,21世紀末全球平均氣溫將上升1.5℃左右[3],這將對人類居住環境造成很大的影響。由于氣候變化,造成全球水文過程發生變化,冰川融化、海平面上升均會改變生物的生存環境。因此,研究氣候變化對水文過程的影響已成為世界水文水資源方向的研究重點[4]。
針對氣候變化對水文過程的影響,國內已有了部分研究。林嫻[5]研究了氣候變化和人類活動對武江流域水文過程的影響,指出氣候變化對水文過程變化的貢獻率高達94%;謝天明[6]分析了在氣候變化條件下塔河流域的徑流變化特征,指出溫度的變化是影響當地徑流變化最關鍵的因素;陳亞寧等[7]對塔里木河三源流區的氣象因子與徑流量的相關分析,認為徑流不同年份間較為穩定主要是因為降水和冰雪融水的互補,氣溫升高是影響徑流變化的最關鍵因素。
本文基于SWAT模型模擬不同氣候條件下對水文過程的影響,深入分析在不同氣候變化的條件下的水文徑流的變化趨勢,得出的結論可為研究區域水資源優化配置提供理論參考,達到區域水資源可持續發展的目的。
已有研究表明,在未來我國平均氣溫將上升2.0℃左右,而降雨將增加10%~20%[8- 9]。選擇山西省太原市為研究區域,以當地水文站數據為基礎,結合研究區域所在位置,本文假設降水量的幾個變化標準為±5%、±10%、±15%和不變的情況,溫度的幾個變化標準為±1℃、±2℃和不變的情況,模擬不同變化條件下水文徑流過程的響應規律。
SWAT模型符合水分循環系統規律[10]:
(1)
式中,sωt—土壤最終含水率,%;sω0—土壤初始含水率,%;Rday—總降水量,mm;Qsurf—地表總徑流量,mm;Ea—總蒸發量,mm;ωseep—土壤測流總量,mm;Qgw—地下徑流總量,mm。
本文模型模擬結果評價指標體系采用相對誤差(RE)、模型效率系數(Ens)和決定系數(R2)綜合分析精度,具體公式如下:
(2)
(3)
(4)

在基于SWAT模型進行水文模擬時,需調整模型參數、初始條件、邊界條件和限制條件,以提高模型模擬的精度,通過反復調試,可得出SWAT模型各參數取值,見表1。其中各參數的含義見文獻[11]。
SWAT模型調參后的模擬結果分析如圖1所示。圖1顯示,SWAT模型的模擬精度較高,模擬值的變化趨勢與實測值的變化趨勢基本一致,均呈現了先增加后減少的周期性變化,同時模擬值與實測值的擬合方程斜率為0.997,與標準值“1”十分接近,表明模擬值表現出了很高的精度,模擬值的決定系數R2為0.971,相對誤差RE僅為9.78%,而模型效率系數Ens達到了0.983,且模擬值與實測值的一致性達到了極顯著水平(P<0.01),通過以上分析可知,模型模擬結果的誤差較低,且一致性較高,驗證了調參后的SWAT模型的可靠性。
基于SWAT模型模擬未來不同氣候變化條件下,豐水年、平水年和枯水年不同水平年下水文徑流變化的響應過程如圖2所示。圖2顯示,降水量和溫度變化對徑流的影響呈反相關關系,即當氣溫降低、降雨增加時,徑流呈現增加的趨勢,在不同水平年,當降水量不變,溫度升高2℃,徑流變化了2%左右,當溫度不變,增加15%的降水量時,徑流變化了10%以上,以上數據表明,降雨對徑流的影響要明顯高于氣溫;當降水量增加幅度為5%、10%、15%時,年徑流的變化率為14.5%、17.3%、18.7%,而當降水量降低幅度為5%、10%、15%時,年徑流的變化率為-10.1%、-10.7%、-11.4%,表明年徑流量隨降水量減小而變化的下降幅度,明顯低于流域年徑流量隨降水量增加而變化的上升幅度,同時,徑流變化在雙重因素綜合影響下的效果更顯著。綜上所述,溫度的變化與降雨的變化均會影響水文過程的變化,而降雨是影響徑流最關鍵的因素。

圖1 模型模擬結果分析

圖2 不同水平年不同氣候變化條件徑流量的影響
根據參數率定后的SWAT模型進行不同氣候變化條件下的徑流過程模擬,模擬的結果可知,當降水量增加時,徑流量呈現增加的趨勢,而當氣溫增加時,徑流量呈現降低的趨勢,同時降水量對徑流量的影響要高于氣溫,而且徑流變化在雙重因素綜合影響下的效果更顯著。本文將降水變化梯度分為了±5%、±10%、±15%和不變的情況,溫度的幾個變化標準為±1℃、±2℃和不變的情況,在今后的研究中,可增加梯度的劃分層次,才能更加精確得出氣候變化對徑流影響的定量分析。