基于GIS的彭州市崩塌滑坡地質災害地形因素分析

范 娟，楊武年,2，南聰強，翁中銀

（1.成都理工大學國土資源信息技術國土資源部重點實驗室遙感與GIS研究所，四川成都610059；2.成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059；3.成都理工大學地球科學學院，四川成都610059）

基于GIS的彭州市崩塌滑坡地質災害地形因素分析

范 娟1，楊武年1,2，南聰強1，翁中銀3

（1.成都理工大學國土資源信息技術國土資源部重點實驗室遙感與GIS研究所，四川成都610059；2.成都理工大學地質災害防治與地質環境保護國家重點實驗室，四川成都610059；3.成都理工大學地球科學學院，四川成都610059）

基于ArcGIS平臺，利用DEM數據資料，選擇位于四川西北部的彭州市作為研究區域，提取了區內高程、地形起伏等地形因子；統計了5.12震后區內39個滑坡點，58個崩塌點，建立了地形地貌與崩塌滑坡地質災害之間的關系：這類地區對應的地貌類型主要是海拔高程較高的山坡地帶。實驗證明，利用GIS技術，結合數字地面模型，進行崩塌、滑坡等地質災害的地形因子相關性分析，結果可靠，對防災減災具有重要的借鑒意義。

ArcGIS；滑坡災害；崩塌災害；高程；地形起伏度

2008年5月12日14時28分，四川省發生里氏8.0級強烈地震，地震造成了巨大的生命財產損失。隨后沿破裂帶發生了大量的余震，形成地震帶。位于余震區兩端的汶川和青川2縣4.0級以上余震次數達147次，占4.0級以上余震總次數的49.49%，6.0級以上余震2縣各有3次。地震同時引發了一系列的崩塌、滑坡、泥石流、地面塌陷、堰塞湖等次生地質災害。隨著空間信息技術的發展,GIS、RS、GPS等技術的應用日趨廣泛。鑒于3S技術應用到地震災害分析有著強大的優勢，能及時準確地獲取相關信息，所以本文選取地震災區彭州市為研究區，利用3S技術對研究區地震前后影像數據做對比，利用DEM提取各地形因子，對彭州市地震后崩塌、滑坡災害點進行了分析，從而得出災害點的區域分布特點，為研究區的災后重建規劃提供了科學依據。

1 研究區概況

彭州市位于四川省會成都西北部，距成都市區25km，人口78萬，面積1420km2，是成都市人口第三，面積第二的小康市。彭州是古蜀國建都立業的核心地區，自秦漢以來，建縣設郡達2000多年，有著悠久的歷史文化，素有“天府金盆”、“蜀漢名區”之美譽。

彭州縣級市，是成都市19個區（市）縣之一。周邊有都江堰縣級市、郫縣、新都區及德陽市的什邡縣級市、廣漢縣級市，阿壩藏族羌族自治州的汶川縣與茂縣。它是成都市管轄的省屬縣級市，設28個鎮，市政府駐天彭鎮。區內4000多m的海拔落差，山丘壩俱全的地形和溫暖濕潤和氣候，決定了彭州豐富的礦藏、林業、藥材、水利和旅游資源。

2 數據與方法

2.1 基于ArcGIS的數據分析流程

本文使用的高程數據為DEM數據，滑坡崩塌數據來源于全國地質災害實地調查資料和中國滑坡崩塌數據庫系統。本研究所用的平臺為ArcGIS9.3。在提取高程、起伏度之前應該對數據進行預處理，包括數據格式轉換、投影變換、崩塌滑坡點的提取及其與DEM高程值的匹配，數據處理流程如圖1所示。

圖1 數據處理流程圖

2.2 運用Excel的統計分析

將彭州市的地圖數據在ArcGIS進行坡度、坡向分析[1-4]。根據提取出的災害點的高程值、起伏度在Excel中進行統計分析，并得出災害點的高程及起伏度的分布特點和規律。

3 滑坡崩塌點的高程分布特征

3.1 滑坡點高程統計分析

通過提取每個滑坡災害點的高程值，并根據具體情況以100 m、200 m、300 m為區間統計滑坡點個數，得到地形高程、滑坡點頻率統計的頻率直方圖（見圖2）。整個研究區的海拔高度在800～3100m之間。通過表1可以看出滑坡災害點主要集中在1 400～1 800 m之間，其間的災害點大約占了研究區總數的75%。

表1 滑坡點高程統計分析表

圖2 滑坡頻率統計直方圖

3.2 崩塌點高程統計分析

通過提取每個崩塌災害點的高程值，并根據具體情況以200m為區間統計崩塌點個數，得到地形高程、崩塌點統計的頻率直方圖（見圖3）。整個研究區的海拔高度在800～2500m之間。通過表2可以看出崩塌災害點主要集中在2個區間段，分別為1 400～1 600 m與2 000～2 200 m。這2個區間段的崩塌災害點大約分別占了研究區總數的24.1%與16.2%。

4 滑坡崩塌與地形起伏度對應關系

滑坡地形起伏度的計算關鍵在于統計單元的選擇，也就是定義所指的某一確定面積的值。隨著統計單元半徑的增大，地形起伏度的值也隨之增大，但增到一定程度即趨于穩定。最佳的統計單元才能真實反映地形起伏[5]。

表2 崩塌點高程統計分析表

圖3 崩塌累計頻率直方圖

對于本文的研究地區，不僅需要計算最適合本研究區的地形起伏度統計單元，而且需要找到能夠反映地形起伏度和滑坡發育關系的最佳統計單元，即這個最佳統計單元不僅反映研究區地形地貌，而且反映滑坡點對應的真實地形起伏度。本文采用10×10鄰域分析窗口，得到研究區的地形起伏分布。利用GRID模塊中SETMASK命令提取出研究區內每一個滑坡災害點所對應的地形起伏度值，并做出地形起伏度值/滑坡個數的統計曲線[4，5]。

地形起伏度的計算，可先求出一定范圍內海拔高度的最大值和最小值,然后求其差值。地形起伏度和滑坡災害的發生存在很大的相關性[6-8]。

4.1 災害點地形起伏度分析

通過提取每個崩塌災害點的起伏度，并根據具體情況以50 m為區間統計崩塌點起伏度，得到崩塌點統計的頻率直方圖（見圖4）。整個研究區的起伏度在100～600m之間。通過圖5可以看出崩塌點的起伏度主要集中在200～300 m區間段，其崩塌災害點大約占了研究區總數的51.3%。崩塌出現在這個高程帶與研究區本身的地形有直接關系。

圖4 崩塌點地形起伏分布頻率直方圖

圖5 災害點地形起伏分布統計圖

4.2 滑坡災害點分布圖

通過提取每個滑坡災害點的起伏度，并根據具體情況以100 m為區間統計崩塌點起伏度，得到滑坡點統計的頻率直方圖（見圖6）。整個研究區的起伏度主要集中在2個區間段。整個區域的起伏度在100～600m之間。通過圖7可以看出滑坡點的起伏度主要集中在150～200m和250～300m區間段。這2個區間段的崩塌災害點大約分別占了研究區總數的28.3%和33.3%。這個結果說明，在區域范圍內地形起伏度和滑坡發育存在很好的相關性，這個結果對滑坡災害空間預測具有重要的應用價值。

圖6 滑坡點地形起伏度頻率統計直方圖

圖7 滑坡災害點分布圖

5 結語

盡管滑坡崩塌的發生是邊坡失穩的結果，與各種外部因素及地質條件有關。但地形地貌參數統計分析表明，滑坡崩塌地質災害在空間分布上與高程、坡度、坡向等地形地貌參數也有著密切的關系[9-11]。在研究區，這種相關性主要體現在：滑坡主要集中發育的高程區間為 1 400～1 800 m，地形起伏度范圍為 150～200 m和250～300m的地區；而崩塌集中發生的高程范圍為1400～1600m與2 000～2 200 m，地形起伏度范圍為200～300 m的地區。由此可以得出崩塌和滑坡常發生的高程區間及地形起伏度范圍。所以我們在建筑物選址的時候應盡量避免將其建在上訴地形區間。只有避開這些區間才能盡量減少地質災害給人們的生活、財產等方面帶來的損失。地質災害防治應貫徹“以防為主，防治結合”的原則，建立和健全地質災害群測群防體系和預警應急反應系統是預防地質災害重要保證。

GIS技術為滑坡崩塌地質災害的地形地貌特征分析提供了實用、有效的方法，對區域地質災害評價具有良好的實用價值。

需要指出的是，滑坡數據屬于歷史調查數據，鑒于目前區域滑坡災害調查的局限性，有很多地區屬于有災無害，可能有許多滑坡點沒有記錄，這會影響到統計結果的精度。

[1] Helena M itasova，Russell S.Harmon，Katherine J.Weaver，et al. Scientific Visualization of L and scapes and L and forms[J].Geomorphology,2012,137(1):122-137

[2] Jiménez-Perálvarez J D,Irigaray C,Hamdouni R El,et al.Building Models for Automatic L and slide-susceptibility Analysis,Mapping and Validation in ArcGIS[J].Natural Hazards,2009,50(3): 571-590

[3] 王慶安,張翔,毛竹.汶川地震生態環境破壞嚴重區地形高程和起伏度分布特點[J].科技導報,2009,27（16）:53-59

[4] 郭芳芳,楊農,張岳橋,等.基于GIS的滑坡地質災害地貌因素分析[J].地質力學學報，2008,14（1）：87-96

[5] 王秋良,薛宏交,楊強,等.北川地區地震災害的地形因子分析[J].大地測量與地球動力學，2008,28（6）:68-72

[6] 郭芳芳,楊農,孟暉,等.地形起伏度和坡度分析在區域滑坡災害評價中的應用[J].中國地質，2008,35(1):131-143

[7] 汪寶山.地質災害特征及影響[J].四川氣象，2003（4）:11-12

[8] 金曉媚,劉金韜.四川省萬縣市滑坡群災害災情評估[J].工程地質學報，1999,7(1):25-29

[9] 周中生,張建華,胡善銓,等.重慶市近年地質災害的特點及防治對策[J].地下空間，2000,20(1):27-31

[10]張建國,魏平新.廣東省主要地質災害發育特點與防治對策[J].中國地質災害與防治學報，2003,14(4):44-48

[11]任玉環,劉亞嵐,魏成階,等.汶川地震道路震害高分辨率遙感信息提取方法探討[J].遙感技術與應用，2009,24(1):52-58

Analysis of Terrain Factors for Landslide Hazards and Collapse Hazards in Pengzhou Based on GIS

by FAN Juan

In this paper,we used GIS technologies and DEM data to analyze Pengzhou geological hazards based on ArcGIS platform.The study area was located in Pengzhou,northwest of Sichuan Province. Aided by ArcGIS,we figured out the terrain factors including elevation and Relief Amplitude,and then analyzed the relationship between landslides,collapse and the terrain factors.The main conclusions are as followed.Such region corresponds to the type of landscape is the hillside area of higher elevation.The results show that the use of GIS technology,combined with digital terrain model,correlation analysis of the digital terrain model of theterrain,landslides and other geologicaldisasters factor,test analysis results are reliable,and have a certain reference to the disaster prevention and mitigation.

ArcGIS,landslides hazards,collapse hazards,elevation, topographic relief amplitude

2012-03-30

項目來源：國家自然科學基金資助項目（41071265）;地質災害防治與地質環境保護國家重點實驗室自主研究基金資助項目(SKLGP2009Z004);2010年度高等學校博士學科點專項科研基金資助項目 (2010512211 0006）。

P208

B

1672-4623(2012)06-0019-03

范娟，碩士，研究方向為地圖學與地理信息系統。