基于SWRI指數的朝陽市水文干旱特征分析

冀連華

(遼寧省朝陽水文局，遼寧 朝陽 122000)

水文干旱是指河川徑流低于其正常值或含水層水位降低的現象，其主要特征是在特定面積、特定時段內可利用水量的短缺。傳統的水文干旱識別觀點認為，如果在一段時期內，流量持續低于某一特定的閥值，則發生了水文干旱，其中閥值的選擇可以依據流量的變化特征，或者根據水需求量來確定；目前，基于該觀點已經衍生了諸多的水文干旱判斷指標，如地表供水指數(SWSI)、標準干旱嚴重指數(SWDSI)、標準徑流指數(SRI)、Parlmer水文干旱強度指數(pHD)等。由以上定義與指標原理可知，傳統判別方法是以流量為單一依據，沒有考慮地下水的影響。在全球氣候變化與人類活動的影響下，地下水在區域供水水源中的比重不斷提高，傳統的水文干旱識別方法已經無法真實反映區域水文干旱情況。針對該問題，翟家齊等人構建了一個能夠反映流域地表水和地下水資源變化的干旱指標，即標準水資源指數(SWRI)，并成功應用該方法評估了海河北系的流域水文干旱特征及其演變規律。以行政區劃為水文干旱研究單元更有利于各地區水資源管理與干旱應對方案的制定，但目前SWRI指數多應用于流域尺度，本文以朝陽市為例，基于SWRI指數對全市各區劃的水文干旱進行識別與特征分析。

1 分析方法

1.1 研究區水資源概況

朝陽市位于遼寧省的西部，屬北溫帶大陸季風氣候區，為半干旱半濕潤易旱地區。全市多年平均降水量96.35億m3，地表水資源量14.71億m3，地下水資源量7.13億m3，水資源總量14.92億m3。以行政分區劃分：朝陽縣的水資源總量最大，占全市的23.8%，龍城區、北票市、建平縣、凌源市、喀左縣的水資源總量所占全市的比重依次為3.3%、21.1%、13.2%、23.3%、12.7%；雙塔區的水資源總量最小，僅占全市的2.6%。由于朝陽市的行政區劃在歷史上曾出現過變動，因此評估單元的劃分均以當前狀況為準。

1.2 水文模型

采用分布式水文模型輸出的評價單元地表和地下水資源量月過程為基本輸入值，構建水文干旱評價指標，即SWRI指數。本文主要采用SWAT(Soil and Water Assessment Tool)，該模型是美國農業部農業研究局在CREAMS和SWRRB模型的基礎上發展起來的，可依據DEM數據、土地利用數據、土壤數據等將研究區域劃分為多個只包含一種土地利用類型和土壤類型的水文響應單元(HRU)，可模擬流域水循環中幾乎所有的水文過程。趙安周等人應用SWAT模型結合Palmer干旱指數提出了干旱分析模型SWAT- PDI，對渭河流域干旱的時空演變規律和發生頻率進行了分析。

1.3 SWRI指數

SWRI指數的確定是首先假設評價單元地表和地下水資源總量服從Γ分布，并計算其概率分布密度函數和概率，然后對概率函數進行正態標準化處理，處理后函數、即式(5)的形態參數Z值即為所求的SWRI值。該方法的主要控制方程包括。

圖2 朝陽市各觀測點地下水位實測與模擬值

xij=Wij

(1)

(2)

(3)

F(xij=0)=k/N

(4)

(5)

式中，wij—第i個評估單元第j月的地表和地下水資源總量；k、N—xij=0的樣本數和總樣本數；β和γ—概率分布密度函數的尺度參數和形狀參數；f(xij)—隨機變量xij的Γ分布密度函數；F(xij

SWRI值與干旱等級的劃分標準可表示為：

(6)

2 分析結果

2.1 模型驗證

分布式水文模型具有一定的不確定性，因此分別采用徑流過程和地下水位對模型進行驗證。

徑流過程驗證采用朝陽站1956～2014年間的徑流量實測數據，結果見圖1。由圖1可知，SWAT模型可準確地模擬出朝陽地區徑流量的變化過程趨勢，數值量級也較為接近，實測的該時間段多年平均徑流量值為72480.73萬m3，模擬值為72512.65萬m3，兩者僅相差0.044%。在各年份中，最大誤差出現于1995年，誤差率為14.67%，即誤差均控制在15%以內，滿足精度要求。

圖1 朝陽站年徑流量實測與模擬值

地下水位驗證采用朝陽地區西合當、西高杖子、郭杖子、太平莊、西湖臺、和東營子6個觀測點的地下水位數據，比較結果見圖2。由圖2可知，SWAT模型同樣可準確地模擬出朝陽地區地下水位的變化過程趨勢，數值量級也同樣較為接近，各觀測點的平均誤差率分別為4.58%、4.27%、2.53%、4.10、5.45%、和5.01%；各觀測點的最大誤差率分別為9.79%、8.07%、5.00%、8.21%、9.94%、和9.87%，均低于10%，滿足精度要求。

2.2 水文干旱識別

基于SWAT模型輸出數據，利用MATLAB提供的函數包對公式組(1)～(5)進行求解，得到朝陽市雙塔區、龍城區、北票市、凌源市、朝陽縣、建平縣、和喀左縣在1956～2015年間的Z值，即SWRI指數值，再將結果按照年份分類匯總，結果見圖3。

圖3 朝陽市SWRI指數值

由圖3可知，朝陽市各地區的SWRI指數值具有基本一致的變化趨勢，約在1993年之前各地的SWRI指數值更為接近，之后則差異較大，主要因為近20年朝陽市社會經濟發展迅速，劇烈的人類活動影響了水資源的自然分布狀態。全市大部分SWRI值介于-0.5～-1.5之間，即干旱等級主要是輕旱和重旱，這與該地區為半干旱地區的特性相一致。總體來看，1981年朝陽市各地區的水文干旱SWRI值均較低，雙塔區、龍城區、北票市、凌源市、朝陽縣、建平縣、喀左縣的SWRI指數(等級)分別為-2.03(特旱)、-2.11(特旱)、-1.62(重旱)、-1.51(重旱)、-1.79(重旱)、-2.06(特旱)、-1.26(中旱)。據記載[8]，1981年朝陽地區發生春夏秋三季連旱，全區嚴重干旱18萬hm2，絕收面積達16.67萬hm2，與本文評估基本保持一致。綜上，采用SWAT模型、基于SWRI指數所評估的水文干旱結果與事實較為相符，即該方法可合理地評估出城市尺度各區劃的水文干旱情況。

2.3 水文干旱特征

基于表3數據對朝陽市的水文干旱特征進行分析，得到朝陽市水文干旱的特征值，結果列于表1。

根據表1結果可知，朝陽市各地區水文干旱最小SWRI指數的平均值為-2.05，其中最小的為-2.23、出現于北票市，干旱等級為“特旱”。水文干旱最大SWRI指數的平均值為3.47，其中最大的為3.95、出現于建平縣，無水文干旱事件。全市的SWRI指數平均值為-0.44，略高于-0.50，因此干旱等級為“無旱”。采用各區劃的平均SWRI指數值表征朝陽市各地區的干旱嚴重程度，按干旱程度由低到高排序分別為龍城區、雙塔區、建平縣、喀左縣、北票市、凌源市、朝陽縣。

表1 朝陽市水文干旱特征

變異系數可表示SWRI指數的時間變化過程離散程度，可由標準差與平均值的絕對值之比來表示，根據表1結果可知，朝陽全市的平均變異系數為2.71，其中龍城區的水文干旱程度離散度最高，為3.17，凌源市最低，為2.38。朝陽全市的SWRI指數傾向率平均值為-0.0255，即在1956～2015年間，全市的SWRI指數平均每年下降0.0255，其中下降最快的為朝陽縣、其傾向率為-0.0272，最慢的為凌源市、其傾向率為-0.0237。

由表1數據可知，朝陽市各地區的SWRI指數傾向率均為負值，說明朝陽市各地區的水文干旱程度均呈加重趨勢。在此情況下，應基于本文或其它有效方法積極加強朝陽市各地區的水文干旱特征分析，對旱情進行準確的預測預報，從而實現由被動抗旱向主動防旱抗旱的轉變。此外，也應繼續加強抗旱工程措施建設，水土保持生態環境建設和抗旱綜合管理體系建設。

3 結論

采用SWAT模型模擬了朝陽市各區劃的地表和地下水資源量變化過程，計算了各地區的SWRI指數，并以此分析了朝陽市的水文干旱特征，結果表明：SWAT模型可準確地模擬出朝陽地區徑流量和地下水位的變化過程趨勢，為SWRI指數的計算提供輸入數據；SWRI指數可準確評估城市尺度各區劃的水文干旱情況；朝陽市各地區按干旱程度由低到高排序分別為龍城區、雙塔區、建平縣、喀左縣、北票市、凌源市、朝陽縣；各地區的水文干旱程度均呈加重趨勢，因此應積極采取合理的措施來應對新的水文干旱形勢。

[1] 吐爾洪, 艾力. 塔河流域干旱災害特征與成因分析[J] . 水利規劃與設計, 2015(11): 18- 21.

[2] 周游. 不同干旱指數在遼寧西部地區的對比及適用性分析[J]. 水利規劃與設計, 2016(01): 43- 45.

[3] 翟家齊, 蔣桂芹, 裴源生, 等. 基于標準水資源指數(SWRI) 的流域水文干旱評估——以海河北系為例[J]. 水利學報, 2015, 46(06): 687- 698.

[4] 胡加琪. 朝陽市干旱成因及減災措施[J]. 東北水利水電, 2006, 24(03): 67- 68.

[5] 江凈超, 朱阿興, 秦承志, 等. 分布式水文模型軟件系統研究綜述[J]. 地理科學進展, 2014, 33(08): 1090- 1100.

[6] 師戰偉. 基于水文模型的大凌河流域綜合干旱指數研究[J]. 水利規劃與設計, 2016(04): 48- 51.

[7] 趙安周, 劉憲鋒, 朱秀芳. 基于SWAT 模型的渭河流域干旱時空分布[J]. 地理科學進展, 2015, 34(09): 1156- 1166.

[8] 解文娟, 楊曉光, 楊婕, 等. 氣候變化背景下東北三省大豆干旱時空特征[J]. 生態學報, 2014, 34(21): 6232- 6243.