雅安市名山區地質災害危險性評價的構建

【摘 要】 地質災害的發生受控于多種影響因子，同時影響因子中各屬性類別對地質災害的發育分別產生促進或抑制作用，因此，構建合適的指標體系對滑坡危險性評價具有重要意義。本文選取雅安市名山區作為研究區，參照前人研究，選取7個影響因子，并對各影響因子及其各屬性類進行敏感性分析，建立評價指標體系，得到名山區危險性分區圖。

【關鍵詞】 地質災害 敏感性分析 危險性評價指標 名山區

本文選擇四川省雅安市名山區作為研究區，獲取名山區112個地質災害隱患點，并參照前人研究選取7個影響因子，根據唐川[1]在對德國波恩地區滑坡因子進行分析時提出的敏感性分析方法，對各影響因子及其各屬性進行敏感性分析，實現對名山區地質災害的危險性評價。

1 研究數據與研究方法

1.1 研究數據

本研究中地質災害隱患點數據來源于四川省雅安市人民政府。基于名山區30m分辨率的DEM數據生成名山區坡度圖、坡向圖及河流水系圖;利用地理空間數據云上下載的2017年7月Landsat 8影像（研究區內均無云）計算NDVI;所用到的名山區的地質年代圖來源于1：20萬的地質圖矢量化;所涉及的路網數據是通過Google Earth下載得到。

1.3 敏感性指數IL分析法

敏感性分析法是利用ArcGIS對屬性數據進行預處理，繼而對地質災害發育的影響因子進行篩選，得出關鍵影響因子，建立區域地質災害危險性評價指標體系。

式中：IL表示地質災害危險性指數;Pmn（L）表示影響因子m的n個屬性類別中存在地質災害隱患點的柵格個數;Pm（L）表示影響因子m中的所有地質災害隱患點的柵格個數;Pmn（G）表示影響因子m的n個屬性類別中不存在地質災害隱患點的柵格個數;Pm（G）表示影響因子m所占的總柵格個數。

2 評價指標體系的構建

2.1 影響因子敏感性分析

本文選取河流、NDVI、坡度、地質年代、坡向、人類工程活動以及高程7個影響因子[2]，并對這7個影響因子的各屬性類進行相對頻率組合定量計算，最終確定其滑坡敏感性指數IL。具體計算過程是：①在ArcGIS中將選取的7個影響因子的屬性用30 m×30 m柵格單元表示;②將名山區內112個地質災害隱患點圖層添加到ArcGIS中不同屬性圖層的柵格單元中，提取地質災害隱患點所處位置中各個影響因子的屬性數據，用柵格單元個數表示，得出影響因子m的n個類別中存在地質災害隱患點的柵格個數Pmn（L）、影響因子m中所有滑坡點的柵格個數 Pm（L）、影響因子m的n個類別中不存在地質災害隱患點的柵格個數Pmn （G）以及影響因子m 所占的總柵格個數 Pm（G） ; ③根據公式（1） 至公式（3） 計算得出各個影響因子的各屬性類別的敏感性; ④根據影響因子各屬性類別的敏感性可得出各影響因子的平均敏感性。

根據敏感性計算方法及Gokceoglu[3]等的研究成果，說明某個因子中某一屬性類別的IL值大于1，則地質災害隱患點與該屬性類別成正相關，其對地質災害隱患點的發育具有較明顯的促進作用，地質災害隱患點在該屬性類內發生的概率較大;若IL值小于1，則反之。當IL值接近于1時，說明該因子與地質災害隱患點相關性不大，接近于整個區域內的平均水平。根據各影響因子各屬性類別的敏感性得出各影響因子的平均敏感性，若某個影響因子的平均敏感性值小于1，則不作為該地區地質災害隱患點危險性評價的參評因子。

3.2 關鍵影響因子篩選

利用各個影響因子的各屬性類別的IL值，對于選的7個影響因子進行平均敏感性計算，得出：地質年代（2.600）>坡向（2.020）>人類工程活動（1.839）>高程（1.139）>河流（0.954）>NDVI（0.796）>坡度（0.619）。由于河流、NDVI和坡度的平均敏感度小于1，因此予以剔除。最終確定地質年代、坡向、人類工程活動以及高程4個因子作為名山區地質災害危險性評價的參評因子。

4 地質災害危險性評價

根據敏感性分析篩選出4個參評因子，并對其進行歸一化和數值化處理，每個因子按照4-5級危險度劃分，由此可作出相應因子的評價圖[4]。

在ArcGIS中利用Raster Calculate的疊加工具進行疊加，得到的數值范圍是0 ～8.5。在ArcGIS中將其重新分類，劃分為無危險區（0-1.759）、輕度危險區（1.759-2.246）、中度危險區（2.426-4.089）以及高危險區（4.089-8.414），最后得到名山區的地質災害危險性區劃圖。

5 結論

（1）本文利用地質災害的敏感性指數對影響名山區地質災害發育的影響因子進行了詳細分析，最終得出，在坡度介于15-35°，距離高速公路1 km以內，在水系200 m以內，NDVI大于0.3，高程在544-650 m、850-1470 m，坡向介于90-135°、225-350°，地質年代為三迭系、下三系、侏羅系以及白堊系，對研究區地質災害的發育起到促進作用。

（2）從得出的結論中可以看出，中、高危險區主要集中在名山區的東南方向和西北方向，占總面積的18.19%，其中高危險區僅占總面積的0.31%，故可得出名山區總體較為安全。

【參考文獻】

[1] 唐川，JorgGrunert.滑坡災害評價原理和方法研究[J].地理學報，1998（S1）：149-157.

[2] 齊識. 白龍江流域滑坡危險度評價技術研究[D].蘭州大學，2014.

[3] C. G?kceoglu，H. Aksoy. Landslide susceptibility mapping of the slopes in the residual soils of the Mengen region （Turkey） by deterministic stability analyses and image processing techniques[J]. Engineering Geology，1996，44（1）.

[4] 唐川，朱靜.基于GIS的山洪災害風險區劃[J].地理學報，2005（01）：87-94.

作者簡介：劉賢（1997-），女，漢族，福建省南平市人，學生，工學碩士，單位：成都理工大學地球科學學院測繪工程專業，研究方向：地圖制圖學與地理信息系統