基于證據權法與確定系數法的滑坡危險區劃對比分析
——以虹口鄉為例

劉光輝，劉漢湖*，姚 昆，夏 潔

(成都理工大學 a.地球科學學院，b.國土資源部地學空間信息技術重點實驗室，成都 610059)

基于證據權法與確定系數法的滑坡危險區劃對比分析
——以虹口鄉為例

劉光輝a,b，劉漢湖a,b*，姚 昆a,b，夏 潔a,b

(成都理工大學 a.地球科學學院，b.國土資源部地學空間信息技術重點實驗室，成都 610059)

都江堰虹口鄉為國家級旅游景區，汶川地震后，該區滑坡頻發，極大威脅該地區游客及居民人身財產安全，開展該地區滑坡危險區預測具有重要意義。這里以地質圖、地形圖、DEM高程數據、QuickBird影像為數據源，選取地層巖性、構造、高程、坡度、坡向、水系6個評價因子，基于GIS平臺，用證據權法、確定系數法對研究區滑坡危險性分析方法進行對比分析，得出最優的滑坡危險區劃研究方法。分析表明：確定系數法的準確率為0.7732，而證據權法的準確率為0.7216。經過野外的調查、研究，確定系數法預測的危險性區域與實際情況更加地吻合，有著較高的精度，因此確定系數法更加適合虹口鄉的滑坡危險區劃研究。

證據權法； 確定系數法； GIS； 滑坡

0 引言

滑坡是一種十分常見的地質災害，合理的滑坡預測方法能增加預測的準確性，降低滑坡造成的危害，減少生命財產的損失[1]。現在主要的滑坡預測方法分為三類：①知識驅動方法；②數據驅動方法；③確定性方法[2]。這里選擇的證據權法與確定系數法都是數據驅動方法。證據權法[2-6]最早用于醫療診斷，后由加拿大數學地質學家Agterberg引入地質領域。確定系數[7-9]最早是由E.H.Shortliffe和B.G.Buchanan提出，后由D.Heckerman進行改進用在災害評價中。兩種方法均為數學的概率學模型，在GIS柵格分析功能[10-11]的支持下，計算出每個評價因子的權重值。證據權法是需要得到每個評價因子的正、負權重值，通過正、負權重的疊加來獲取評價因子的最終權重值。確定系數法是利用一個概率函數直接計算出每個評價因子的權重值，且這個權重值的值域是[-1,1]。前人在滑坡危險性區劃研究中并沒有對確定系數法與證據權法進行預測準確率的對比分析，這里以虹口鄉為研究區，分析出證據權法與確定系數法在虹口鄉這種地形復雜，地質條件特殊，滑坡災害頻發的地區，該方法有更高的滑坡預測精度。證據權法與確定系數法，更加地適合虹口鄉的滑坡危險區劃研究。

1 數據

虹口鄉位于四川省都江堰市東北的白沙河中游，東接彭州市，南連向峨鄉、金鳳鄉、白沙鄉，西南靠龍池鎮，西北和北鄰阿壩州汶川縣(圖1)。研究區內水系單一，水網密集，在地形上中間沿河流起伏較小，向兩側坡度變陡。在地質構造上，受龍門山斷裂帶的影響，虹口鄉溝谷發育，坡陡坡長，厚層砂巖及厚層礫巖場形成后陡坎或陡崖，斜坡穩定性較差，規模不大的崩塌、滑坡比較發育，處于斷裂構造部位有大型的滑坡或崩塌。虹口鄉在地形，地質，地貌等條件上都容易誘發滑坡，特殊的地質條件導致虹口鄉成為地質災害的多發區，其地質災害主要以崩塌和滑坡為主。

圖1 虹口鄉滑坡區域位置圖Fig.1 Location map of Hongkou landslide area

本次滑坡的危險性區劃研究數據包括：1∶50 000 地質圖、1∶50 000 地形圖、ASRGTM 30 m空間分辨率數字高程模型、Quickbird影像。

2 研究方法

2.1 證據權法

證據權法是一種綜合各種證據來支持一種假設的定量方法，該方法以貝葉斯統計模型為基礎，通過對已經發生的滑坡與影響滑坡形成的影響因子( 坡度、坡向、地層巖性等) 進行空間關聯分析，求取各個評價因子對滑坡發生的貢獻率( 即權重) ，最后將各個評價因子賦予相應權重并疊加，得到地質災害危險性指數，該指數高低代表該區發生地質災害可能性的大小。其計算公式如式(1)～式(3)所示。

(1)

(2)

Wf=W+-W-

(3)

2.2 確定系數法

確定系數法是一種結合確定性系數(CF)來分析影響滑坡發生的各個因素敏感性的方法。確定系數法是一種利用概率函數的方法，其計算公式見式(4)。

(4)

式中：PPa為影響因子分級中事件發生的條件概率，在這里可表示為影響因子各級范圍內滑坡發生的柵格數與所有柵格數的比值；PPs為整個研究區內事件發生的先驗概率，表示為研究區內滑坡發生的柵格數與研究區總柵格數的比值。

2.3 評價因子

基于前人研究及歷史滑坡數據的綜合分析[12-15]，主要選取斷裂帶、地層巖性、水系、坡度、坡向、高程6個評價因子。并將各個評價因子分成不同的等級，見表1。

3 結果分析

3.1 權重計算

為便于權重的計算，將研究區的評價單元分成5 m×5 m的柵格，再根據證據權法和確定系數法分別計算出每個評價因子的權重值，計算結果見表2。

3.2 危險性指數計算

因為證據權法和確定系數法都是數據驅動方法，利用概率方法建立模型進行運算，每個評價因子之間必須滿足條件獨立性的假設。計算之前首先驗證各個評價因子兩兩之間的獨立性。這里選用的是卡方(X2)檢測方法對各個評價因子進行獨立性檢測。檢測結果表明評價因子之間均滿足獨立性條件。利用GIS的空間分析功能將各個評價因子進行加權分析，得到每個評價單元的危險性指數(LHI),其計算公式如下：

表1 評價因子及屬性分級

表2 滑坡危險性權重計算表

LHI=Wf巖性+Wf水系+Wf坡度+Wf高程+

Wf坡向+Wf斷裂帶

(5)

LHI=CFf巖性+CFf水系+CFf坡度+

CFf高程+CFf坡向+CFf斷裂帶

(6)

運用證據權法將6個評價因子按照公式(5)進行計算，得到滑坡的危險性指數的范圍是-4.958～7.1。而運用確定系數法將6個評價因子按照公式(6)進行計算，得到滑坡的危險性指數范圍為-3.391～2.605。

滑坡危險性評價的成功率檢驗采用成功率曲線法，這是一種定量的滑坡檢測方法(圖2)。在圖2中AUC曲線下面積越大，說明滑坡評價單元危險性指數越高，即模型預測效果越好。由圖2可以計算出，證據權法的曲線下面積為72.16，確定性系數法的曲線下面積為77.32。由此可知，虹口鄉的滑坡危險區劃研究中，確定系數法預測的結果要比證據權法預測的結果成功率高。

圖2 成功率曲線圖Fig.2 Success rate curve of landslides hazards

3.3 危險性分區

根據證據權法與確定系數法兩種方法計算得到的綜合危險性指數，將研究區滑坡危險性指數分為極高、高、中、低4類，并統計每類滑坡的面積百分比(表3)。

從表3可以看出：證據權法分析結果表明極高危險區面積百分比為21.7%，滑坡面積的76.8%分布在高危險區與中危險區；確定系數法分析結果表明高危險區所占的比例為31.3%，滑坡面積的80%分布在高危險區與中危險區。在研究區的滑坡危險性區劃研究中確定系數法的精確度高于證據權法，下面給出確定系數法的滑坡危險區劃的分區圖(圖2)。研究區總共解譯出滑坡117處，極高危險區有34處滑坡，占滑坡總個數的29%，高危險區有53處滑坡，占總滑坡個數的45%。

表3 證據權法與確定系數法危險區劃對比圖

圖3 研究區確定系數法危險性分區圖Fig.3 Risk zoning map of the study area determin-ation coefficient method

4 結論

滑坡預測準確率直接影響生命財產安全，針對不同區域的不同地形、地質、地貌的特點，滑坡的預測應采用適合的評價模型提高預測準確率。這里選取證據權法、確定系數法2種研究模型，對研究區域滑坡危險性區劃進行研究分析。這2種方法都是數據驅動方法，在各個評價因子選擇及權重賦值的過程中有效地減小人為影響。研究結果表明：

1)滑坡的分布與水系、坡度、巖性有著明顯的關系。滑坡在距水系小于400 m時，水系對滑坡的影響較大，占總滑坡數的65.2 %，特別是在100 m～300 m之間影響最大，占總滑坡數的34 %。坡度對滑坡的影響在30°～40°之間時影響最大，其滑坡數量占總滑坡的65.7%。從圖3可以看出，在圖3的左面滑坡分布較多，該處巖層風化厚度較大，坡度主要為30°～40°，且大多斜坡為順層斜坡，其滑坡數占總滑坡的77.9%。

2)在研究區的分區結果中，高危險區、中危險區所占面積較大，證據權法、確定系數法分析結果面積百分比分別達到76.8%和80%，證據權法的高危險區所占比例為21.7%，確定系數法的高危險區所占比例為31.3%，在高危險區確定系數法的預測準確率明顯高于證據權法。

3)由成功率曲線法得到確定系數法的準確率為0.773 2，而證據權法的準確率為0.721 6，從成功率曲線可以看出，在危險系數從10%～40%之間，確定系數法的曲線明顯的高于證據權法，預測出的滑坡數確定系數法比證據權法多。確定系數法的預測結果成功率相較于證據權法更高，從整體分析確定系數法更加適用于研究區的滑坡危險區劃研究。

[1]PARISE M,JIB SON R W.A seismic landslide susceptibility rating of geologic units based on analysis of characteristics of landslides triggered by the17 January,1994 Northridge,California earthquake[J].Engineering Geology,2000,58(3/4):251-270.

[2]范強,巨能攀,向喜瓊,等.證據權法在區域滑坡危險性評價中的應用——以貴州為例[J].工程地質學報,2014,03(17):124-129.FAN Q,JU N P,XIANG X Q,et al.Landslides susceptibility assessment using weights of evidence in the Tongzi watershed,China[J].Journal of Catastrophology,2015,30(17):124-129.(In Chinese)

[3]王志旺,李瑞有,王湘桂.證據權法在滑坡危險度區劃研究中的應用[J].巖土工程學報，2007,29(8):1268-1273.WANG Z W,LI D Y,WANG X G.Zonation of landslide hazards based on weights of evidence model[J].Chinese Journal of Geotechnical Engineering,2007,29(8):1268-1273.(In Chinese)

[4]張艷玲,南征兵,周平根.利用證據權法實現滑坡易發性區劃[J].水文地質工程地質,2012,39(2):121-125.ZHANG Y L,NAN Z B,ZHOU P G.Division of landslide susceptibility based on weights of evidence model[J].Hydrogeology & Engineering Geology,2012,39(2):121-125.(In Chinese)

[5]許沖,戴福初,徐錫偉.基于GIS平臺與證據權的地震滑坡易發性評價[J].地球科學：中國地質大學學報,2011，36(6):1155-1164.XU C,DAI F C,XU X W.Earthquake triggered landslide susceptibility evaluation based on GIS platform and weight of evidence modeling[J].Earth Science-Journal of China University of Geosciences,2011,36(6):1155-1164.(In Chinese)

[6]許沖,徐錫偉,于貴華.基于證據權方法的玉樹地震滑坡危險性評價[J].地震地質,2013,35(1):151-164.XU C,XU X W,YU G H.The Yushu earthquake triggered landslide hazard evaluation based on weight of ecidence method[J].Seismology and Geology,2013,35(1):151-164.(In Chinese)

[7]許沖,戴福初,姚鑫.基于GIS與確定性系數分析法的汶川地震滑坡易發性評價[J].工程地質學報，2010,18(1):15-26.XU C,DAI F C,YAO X.Study on Wenchuan earthquake-induced landslide susceptibility evaluation based on GIS platform and certainty factor analysis method[J].Journal of Engineering Geology,2010,18(1):15-26.(In Chinese)

[8]劉麗娜,許沖,陳劍.GIS支持下基于CF方法的2013年廬山地震滑坡因子敏感性分析[J].工程地質學報，2014,22(06):1176-1186.LIU L N,XU C,CHEN J.Landslide factor sensitivity analyses for landslides triggered by 2013 lushan earthquake using GIS platform and certainty factor method[J].Journal of Engineering Geology,2014,22(6):1176-1186.(In Chinese)

[9]許沖,戴福初,姚鑫.基于GIS的汶川地震滑坡災害影響因子確定性系數分析[J].巖石力學與工程學報,2010,29(SI):1-9.XU C,DAI F C,YAO X.Gis based certainty factor analysis of landslide triggering factors in WenChuan earthquake[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(SI):1-9.(In Chinese)

[10]溫守欽,李仁鋒,任群智.GIS技術在地質災害區劃中的應用[J].中國地質,2005,32(3):512-517.WEN S Q,LI R F,REN Q Z.Application of GIS technology in the regionalization of geological hazards[J].Geology in China,2005,32(3):512-517.(In Chinese)

[11]戴福初,李軍.地理信息系統在滑坡災害研究中的應用[J].地質科技情報,2000,19(1):91-96.DAI F C,LI J.Applications of geographical information systems in landslide studies[J].Geological Science and Technology Information,2000,19(1):91-96.(In Chinese)

[12]陳曉利,冉虹流,祁生文.1976年龍陵地震誘發滑坡的影響因子敏感性分析[J].北京大學學報,2009,45(1):104-110.CHEN X L,RAN H L,QI S W.Triggering factors susceptibility of earthquake-indueced landslides in 1976 Longling Earthquake[J].Acta Scientiarum Naturalium Universitatis Pekinensis,2009,45(1):104-110.(In Chinese)

[13]郭芳芳,楊農,孟暉,等.地形起伏和坡度分析在區域滑坡災害評價中的應用[J].中國地質,2008,35(1):131-143.GUO F Fang,YANG N,MENG H,et al.Application of the relief amplitude and slope analysis to regional landslide hazard assessments[J].Geology in China,2008,35(1):131-143.(In Chinese)

[14]陳曉利,鄧儉良,冉虹流.汶川地震滑坡崩塌的空間分布特征[J].地震地質,2011,33(1):191-202.CHEN X L,DENG J L,RAN H L.Analysis of landslides triggered by wenchuan earthquake[J].Seismology and Geology,2011,33(1):191-202.(In Chinese)

[15]蘭恒星,伍法權,周成虎.基于GIS的云南小江流域滑坡因子敏感性分析[J].巖石力學與工程學報,2002,21(10):1500-1506.LAN H X,WU F Q,ZHOU C H.Analysis on susceptibility of GIS based landslide triggering factors in Yunnan Xiao jiang watershed[J].Chinese Joumal of Rock Mechanices and Engineering,2002,21(10):1500-1506.(In Chinese)

Based on the weights of evidence method and certainty factor method of landslide danger division comparative analysis——A case study in Hongkou,chine

LIU Guang-huia,b,LIU Han-hua,b*,YAO Kuna,b,XIA Jiea,b

(Chengdu University of Techndogy,a.College of Eong sciences;b.key Laboratory of Geological spatial Information Technology of Ministry of Land and Resources,Chengdu 610059,China)

Hongkou township in Dujiangyan is a national scenic spot.After Wenchuan earthquake,landslides frequently occur in this area,which greatly threatens personal and property safety of tourists and residents there.Thus,it becomes significantly important to develop forecast to the landslide dangerous area of this region.This article based on GIS platform,takes geologic map,topographic map,DEM elevation data,QuickBird image as data sources,selects 6 evaluation factors of stratigraphic lithology,structure,elevation,slope,slope direction and water system,applies to weight of evidence method,determination coefficient method to make comparative analysis to the analytical methods of landslide dangers in research area and find out the best research method to landslide risk area.The analysis shows that,accuracy rate of determination coefficient method is 0.7732,while weight of evidence method has 0.7216 accuracy rate.After field research,the risk area predicated by method of determination coefficient can match practical situation more with high accuracy.Thus determination coefficient method is more fitness for the research of landslide risk area in Hongkou Township.

weights of evidence; certainty factor; GIS; landslides

2015-10-28 改回日期：2015-12-03

國家自然科學基金(41102225)

劉光輝(1988- )，男，碩士,主要研究方向為遙感地質，E-mail：lghyaogan@126.com。

*通信作者：劉漢湖(1978-)，男，博士，副教授，主要研究方向為遙感地質，E-mail：Liuganhu@cdut.cn。

1001-1749(2016)06-0848-06

TU 431