氣候變化背景下濟南市長清區干旱災害風險評價研究

摘要：綜合干旱類型，從致災因子的危險性、敏感性、易損性及防災減災能力4個指標出發，分析氣候變化背景下的山東省濟南市長清區干旱災害風險變化。結果表明，長清區干旱危險性分布特征總體表現為由西北向東南遞減的趨勢；干旱災害孕災體敏感性總體呈現出由西南向東北遞減的趨勢；干旱災害易損性分布基本以行政區為單位，高易損性分布在除平安街道辦事處外的西部地區；干旱災害的防災減災能力總體呈現東北高、西南低的分布特征；干旱綜合風險性呈現西南高、東北低的分布規律。

關鍵詞：長清區；區域干旱災害；風險評價研究

中圖分類號：P429；S423 文獻標識碼：A 文章編號：0439-8114（2016）20-5233-05

DOI：10.14088/j.cnki.issn0439-8114.2016.20.017

Abstract： According to drought types， using 4 indicators of hazard， sensitivity， vulnerability and disaster prevention and mitigation capacity， the change of drought disaster risk was analyzed in the context of climate change in changqing district， Ji′nan city， Shandong province. The results showed that， in changqing district， drought hazard overall decreased from the northwest to the southeast. Sensitivity decreased from southwest to northeast. Distribution of vulnerability was based on the administrative area as a unit， high vulnerability distributed in western region except Ping′an street. Disaster prevention and mitigation capacity of drought disaster was higher in the northeast and lower in southwest. Drought comprehensive risk was higher in southwest and lower in Northeast.

Key words：changqing district； regional drought disaster； risk assessment research

近幾十年來，世界各地的自然災害頻繁發生，大范圍的異常氣候事件時有出現，給人民的生命財產和社會經濟帶來嚴重損失，也給世界的糧食、水資源和能源安全造成了一系列危機。其中干旱、洪澇是氣象災害中最嚴重的，尤其是干旱，因其出現頻率高、持續時間長、波動范圍大等，對國民經濟特別是對農業生產有嚴重影響，歷來受到人們的廣泛關注。干旱及災害的危害不容忽視，在對歷史干旱分析的基礎上，了解和掌握干旱及災害的規律、進行災害風險評估以及災害預測預警來積極應對頻繁發生的干旱災害、減少旱災損失是一項十分重要且很有意義的研究工作，這也是以保證水資源安全來達到保證糧食安全和經濟社會快速發展需求的目的。本研究從致災因子的危險性、敏感性、易損性及防災減災能力4個指標出發，分析氣候變化背景下的山東省濟南市長清區干旱災害風險變化。

1 研究方法

1.1 災害風險評估原理

基于自然災害風險形成理論，氣象災害風險是由危險性（致災因子）、敏感性（孕災環境）、易損性（承災體）和抗災能力4部分共同形成[1]。每個因子又是由一系列子因子組成。其表達式為：

災害風險=f[危險性（可能導致災害的因素均可稱為致災因子）、敏感性（大氣環流和天氣系統、水文系統、土壤因素、地形地貌、植被）、易損性（致災因子作用的對象，是蒙受災害的實體）、防災減災能力（防御和減輕氣象災害的各種管理措施和對策，包括管理能力、減災投入、資源準備等）] （1）

1.2 因子標準化

由于所選因子的量綱不同，因此要將因子進行標準化[1，2]。本研究根據具體情況，采用極大值標準化和極差標準化方法。

極大值標準化：

極差標準化：由于研究中所選指標與因子成反比，且低溫存在負值，對其標準化時將極差標準化方法進行修正，將原公式分子中的最小值修訂為最大值，使得結果符合實際。

式中，Xij為第i個因子的第j項指標；X′ij為去量綱后的第i個因子的第j項指標；Xmin、Xmax為該指標的最小值和最大值。

1.3 加權綜合評價法

加權綜合評價法綜合考慮了各個因子對總體對象的影響程度，是把各個具體指標的優劣綜合起來，用一個數值化指標加以集中，表示整個評價對象的優劣，因此，這種方法特別適用于對技術、策略或方案進行綜合分析評價和優選，是目前最為常用的計算方法之一[3]。

式中，Cvj為評價因子的總值，Qvij為對于因子j的指標i（Qij≥0），Wci為指標i的權重值（0≤Wci≤1），通過層次分析法（AHP）計算得出，m是評價指標個數。

對于綜合風險指數，表達式為：

式中，Cv為綜合風險指數，Cvj為危險性、敏感性和易損性風險因子，Cvk為f防災減災能力。

1.4 層次分析法

層次分析法（Analytic Hierarchy Process，AHP）是對一些較為復雜、模糊的問題做出決策的簡易方法，它特別適用于那些難于完全定量分析的問題[3，4]。它是美國運籌學家、匹茲堡大學Saaty教授于20世紀70年代初提出的一種簡便、靈活而又實用的多準則決策方法。層次分析法是一種定性與定量相結合的決策分析方法。決策法通過將復雜問題分解為若干層次和若干因素，在各因素之間進行簡單的比較和計算，就可以得出不同方案重要性程度的權重，為最佳方案的選擇提供依據。

1.5 氣象資料插值方法

根據實際情況，分別采用線性回歸方法、趨勢面擬合方法及Kriging插值方法。其中線性回歸方法是對氣象站中所缺省數據進行的插值，Kriging插值是對氣象因子進行的GIS空間插值方法[5]。原理分別如下：

1.5.1 線性回歸方法 線性回歸模型描述兩個要素之間的線性相關關系，如兩個要素間存在顯著的相關關系，就可以建立二者之間的線性回歸方程，其表達式為：

1.5.2 趨勢面擬合方法 趨勢面分析是利用數學曲面模擬地理系統要素在空間上的分布及變化趨勢的一種數學方法，實質上是通過回歸分析原理，運用最小二乘法擬合一個二維非線性函數，模擬地理要素在空間上的分布規律，表達式為：

1.5.3 Kriging插值方法 Kriging就是根據一個區域內外若干信息樣品的某些特征數據值，對該區域作出一種線性無偏和最小估計方差的估計方法。從數學角度來說，是一種求最優線性無偏內插估計量的方法。該方法的適用范圍為區域化變量存在空間相關性，即如果變異函數和結構分析的結果表明區域化變量存在空間相關性，則可以利用Kriging方法進行內插或外推。其實質是利用區域化變量的原始數據和變異函數的結構特點，對未知樣點進行線性無偏、最優估計。克里格方法是通過對已知樣本點賦權重來求得未知樣點的值，表示為：

式中，Z（x0）為未知樣點的值，Z（x1）為未知樣點周圍的已知樣本點的值，？棕i為第i個已知樣本點對未知樣點的權重，n為已知樣本點的個數。與傳統插值法最大的不同是，在賦權重時，克里格方法不僅考慮距離，而且通過變異函數和結構分析，考慮了已知樣本點的空間分布及與未知樣點的空間方位關系。

2 干旱災害風險評價指標

長清區干旱災害風險區劃危險性、敏感性、易損性、防災減災能力4個指標[6-8]選取及權重值見圖1，其權重值通過層次分析法（AHP）確定[4]。

在干旱災害孕災體敏感性研究中，根據土壤質地和土地利用類型對干旱的影響進行評分，評分見表1和表2。

3 干旱災害危險性因子評價

干旱災害危險性等級劃分標準見表3。

由圖2可知，干旱危險性分布特征總體表現為由西北向東南遞減的趨勢，危險性整體水平不高，高、中危險性面積不大，分別占2.73%、16.38%。北部的邊緣帶狀區域基本屬于高、中危險性地區。高危險性主要在孝里鎮西南部；中危險性主要分布在孝里鎮東北部，歸德鎮、文昌街道辦事處和平安街道辦事處的西部；輕、低危險性主要分布在長清區的中部和東部，分別占32.34%、48.55%，主要包括萬德鎮、張夏鎮、五峰山街道辦事處、馬山鎮、雙泉鄉、崮云湖街道辦事處和歸德鎮、文昌街道辦事處、平安街道辦事處的東部。

4 干旱災害孕災環境敏感性評價

干旱孕災環境敏感性等級劃分標準見表4。

由圖3可知，干旱災害孕災體敏感性總體呈現出由西南向東北遞減的趨勢，分布空間性較強，沒有規則的分布界限。高孕災體敏感性主要分布在西南部邊緣地區，包括孝里鎮、雙泉鄉南部、歸德鎮西南角，占長清區總面積的18.74%；馬山鎮、歸德鎮大部分地區、雙泉鄉北部、萬德鎮西部、張夏鎮西南部以及文昌街道辦事處局部地區為中敏感性區域，占長清區總面積的34.02%；長清區東北部為輕、低敏感性區域，輕敏感性區域主要分布在萬德鎮中部、五峰山街道辦事處大部、張夏鎮東部、崮云湖街道辦事處西部和中部、平安街道辦事處東南部，占長清區總面積的30.85%；低敏感性區域占長清區總面積的16.40%，主要分布在萬德鎮東部、張夏鎮北部。崮云湖街道辦事處東南部和平安街道辦事處的局部地區。

5 干旱災害承災體易損性評價

干旱災害承災體易損性等級劃分標準見表5。

由圖4可知，干旱災害易損性分布基本以行政區為單位。高易損性分布在除平安街道辦事處外的西部地區，主要有文昌街道辦事處、歸德鎮、孝里鎮，占全區面積的31.29%；這些地區易損性高，其原因主要是由于這些地區的糧食和蔬菜的面積較大，人口密度較高。其他地區相比而言易損性較低。中易損性主要分布在平安街道辦事處、張夏鎮、萬德鎮，占全區面積的37.97%。中部地區的干旱災害易損性相對較低，馬山鎮為低易損性，占全區面積的7.11%。崮云湖街道辦事處、五峰山街道辦事處、雙泉鄉屬于輕易損性，占全區面積的23.64%。

6 干旱災害防災減災能力

干旱災害防災減災能力等級劃分標準[9]見表6。

由圖5可知，干旱災害的防災減災能力總體呈現東北高、西南低分布特征，其中東部的平安街道辦事處、文昌街道辦事處、歸德鎮、萬德鎮防災減災能力最高，為高防災減災能力區，占46.60%；中防災減災能力區分布在崮云湖街道辦事處、張夏鎮，占全區的19.70%，這些地區相比弱、低防災減災能力地區應對干旱災害的抵御、恢復能力強，潛在損失較小。弱防災減災能力區包括五峰山街道辦事處、馬山鎮、孝里鎮。雙泉鄉的防災減災能力最低。

7 干旱災害綜合風險評價

干旱災害綜合風險等級劃分標準見表7。

由圖5可知，干旱綜合風險性呈現西南高、東北低的分布規律。孝里鎮、雙泉鄉綜合風險性最高，占全區面積的20.30%。這些地區發生干旱災害的可能性最大，可能受到的潛在損失最大，恢復能力較差。中風險性分布在文昌街道辦事處、五峰山街道辦事處和歸德鎮的西南部局部地區，中風險性面積占23.17%；低風險性在全區分布面積占23.17%，主要分布在歸德鎮大部地區和張夏鎮；輕風險性區域占34.25%，分布面積最大，主要分布在萬德鎮、崮云湖街道辦事處、平安街道辦事處。

8 小結與討論

長清區干旱危險性分布特征總體表現為由西北向東南遞減的趨勢；干旱災害孕災體敏感性總體呈現出由西南向東北遞減的趨勢；干旱災害易損性分布基本以行政區為單位，高易損性分布在除平安街道辦事處外的西部地區；干旱災害的防災減災能力總體呈現東北高、西南低的分布特征；干旱綜合風險性呈現西南高、東北低的分布規律。

了解旱災對社會經濟和生態環境帶來的損失，能為今后制定區域發展規劃提供必要的支持，因而對旱災損失進行快速、精確評價的要求日益迫切。對生態環境脆弱帶的研究起步較晚，對外界干擾的抵抗能力較低、受到外力作用發生退化后恢復得較慢的地區，受人為和自然雙重作用的強烈影響，構成了獨特的孕災環境。當前干旱災害的風險評價方興未艾，研究中還存在一些亟待解決的問題。通常各種自然災害并非獨立發生，而是在某一種自然災害發生后誘發一種或幾種自然災害，而且彼此相互作用，形成復雜的災害鏈或災害群，這就是自然災害的群發性。因此，進行風險評估時，不僅要對單災種進行研究，更重要的是從系統的角度綜合研究區域災害，通過運用災害風險評價原理、因子標準化、加權綜合評價法、層次分析法、氣象資料插值系統的研究討論，提高干旱災害風險評價的可操作性。

參考文獻：

[1] 羅 培.基于GIS的重慶市干旱災害風險評估與區劃[J].中國農業氣象，2007，28（1）：100-104.

[2] 李 艷.河南省干旱承險脆弱性綜合評價研究[D].鄭州：鄭州大學，2011.

[3] 林 霞.遼寧省氣象災害風險評價[D].遼寧大連：遼寧師范大學，2009.

[4] 楊帥英，郝芳華，寧大同.干旱災害風險評估的研究進展[J].安全與環境學報，2004，4（2）：79-82.

[5] 蘇占勝，陳曉光，黃 峰，等.基于GIS的寧夏氣候要素推算及農業氣候資源分析[J].干旱地區農業研究，2008，26（4）：242-248.

[6] 唐 明.旱災風險分析的理論探討[J].中國防汛抗旱，2008（1）：38-40.

[7] 趙 靜，張繼權，嚴登華，等.基于格網GIS的豫北地區干旱災害風險區劃[J].災害學，2012，27（1）：55-58.

[8] 楊豐政，景元書.基于GIS的徐水縣干旱災害風險區域[J].氣象與減災研究，2012，35（1）：54-58.

[9] 張 會，張繼權，韓俊山.基于GIS技術的洪澇災害風險評估與區劃研究——以遼河中下游地區為例[J].自然災害學報，2005， 14（6）：141-146.

[10] GB/T20481-2006，氣象干旱等級[S].