羅路廣 裴向軍 黃潤(rùn)秋 裴 鉆 朱 凌

(地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué))， 成都 610059， 中國(guó))

0 引 言

2017年九寨溝“8·8”地震是四川省繼汶川MS8.0級(jí)地震、蘆山MS7.0地震以來(lái)的又一次強(qiáng)震(Fan et al.，2018)，誘發(fā)了大量的崩滑地質(zhì)災(zāi)害，邊坡巖土體結(jié)構(gòu)受到不同程度的震裂損傷，加之區(qū)內(nèi)降雨量大、人類(lèi)活動(dòng)頻繁等特點(diǎn)使失穩(wěn)概率急劇增高。震后九寨溝景區(qū)經(jīng)過(guò)短暫的開(kāi)放后，受強(qiáng)降雨誘發(fā)崩塌、滑坡、泥石流等地質(zhì)災(zāi)害的威脅，再次臨時(shí)閉園。因此，基于GIS平臺(tái)的滑坡易發(fā)性評(píng)價(jià)可以為危險(xiǎn)范圍預(yù)測(cè)及重點(diǎn)防治區(qū)域區(qū)劃等方面提供理論依據(jù)與技術(shù)支撐。

滑坡易發(fā)性評(píng)價(jià)為滑坡危險(xiǎn)性和風(fēng)險(xiǎn)評(píng)價(jià)的基礎(chǔ)，它是指特定區(qū)域內(nèi)多個(gè)影響因子組合導(dǎo)致滑坡發(fā)生的可能性。國(guó)內(nèi)外學(xué)者利用某種單一數(shù)學(xué)方法及模型進(jìn)行地質(zhì)災(zāi)害易發(fā)性分析和評(píng)價(jià)(胡瑞林等， 2013)，如信息量模型、確定性系數(shù)法(CF)、支持向量機(jī)法(SVM)、證據(jù)權(quán)法及邏輯回歸法(Logistic)、層次分析法(AHP)、頻率比模型(FR)等，均取得了較好的應(yīng)用效果(許沖等， 2009； Pradhan et al.，2010； 張艷玲等， 2012； 劉麗娜等， 2014； 韓振華等， 2015； 許英姿等， 2016； 戴嵐欣等， 2017)，然而仍存在易發(fā)性影響因子的量化、各因子之間及同一因子不同特征值之間對(duì)地質(zhì)災(zāi)害易發(fā)性的影響值確定的人為因素干擾等缺陷。因此，兩種方法甚至多方法、多模型相互耦合是近年來(lái)新的探究方向，如基于地理加權(quán)回歸模型的云南省地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)(饒品增等， 2017)、采用聚類(lèi)分析和最大似然法的汶川極重災(zāi)區(qū)震后地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)(胡凱衡等， 2012)、基于CF和SVM建立的CF-SVM耦合易發(fā)性評(píng)估模型(李遠(yuǎn)遠(yuǎn)等， 2018)、加權(quán)信息量模型(范林峰等， 2012)、信息量與Logistic回歸耦合模型(樊芷吟等， 2018)等。沈玲玲等(2016)應(yīng)用模糊邏輯法、信息量模型及Shannon熵改進(jìn)的信息量模型，對(duì)岷縣MS6.6級(jí)地震誘發(fā)滑坡進(jìn)行了易發(fā)性評(píng)價(jià)； 黃發(fā)明等(2018)為保證SVM等機(jī)器學(xué)習(xí)模型所選的非滑坡柵格單元是真正的“非滑坡”，提出了基于自組織映射(SOM)的聚類(lèi)分析和SVM耦合的滑坡易發(fā)性評(píng)價(jià)模型，結(jié)果顯示SOM-SVM模型具有比單獨(dú)SVM模型更高的成功率和預(yù)測(cè)率； 郭子正等(2019)基于邏輯回歸-模糊層次分析方法(LR-FAHP)的加權(quán)頻率比模型(WFR)，對(duì)三峽庫(kù)區(qū)萬(wàn)州區(qū)內(nèi)滑坡進(jìn)行易發(fā)性等級(jí)預(yù)測(cè)，與單一的LR、FAHP和FR模型相比，WFR模型能將滑坡易發(fā)性評(píng)價(jià)精度提升4%～9%。

諸多研究結(jié)果表明，相比單一模型，多模型耦合在評(píng)價(jià)精度、合理性及成功率方面具有明顯的優(yōu)越性。傳統(tǒng)確定性系數(shù)模型可客觀地反映評(píng)價(jià)因子內(nèi)部不同級(jí)別對(duì)滑坡易發(fā)性程度的貢獻(xiàn)值，但是忽略了各因子之間關(guān)系的復(fù)雜性及其對(duì)滑坡易發(fā)性影響的差異性，不能很好地確定因子間的相對(duì)權(quán)重，而邏輯回歸方法根據(jù)影響因子與歷史災(zāi)害點(diǎn)的關(guān)系基于統(tǒng)計(jì)結(jié)果可較好地確定滑坡易發(fā)性影響因子之間的相互權(quán)重大小。因此，本文基于遙感影像解譯和野外調(diào)查成果，分析地質(zhì)災(zāi)害空間發(fā)育分布規(guī)律及其影響因素，選取地震峰值加速度(PGA)、距斷層距離、高程、坡度、坡向、坡面曲率、地層巖性、公路及水系等9個(gè)因子，將確定性系數(shù)模型與邏輯回歸模型相結(jié)合開(kāi)展九寨溝景區(qū)滑坡易發(fā)性評(píng)價(jià)，為次生災(zāi)害防災(zāi)減災(zāi)提供參考。

1 滑坡易發(fā)性評(píng)價(jià)方法

1.1 確定性系數(shù)模型

確定性系數(shù)模型是由Shortliffe et al.(1975)提出的一個(gè)概率函數(shù)，由于評(píng)價(jià)流程相對(duì)簡(jiǎn)單，精度較高，在地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)中應(yīng)用較為廣泛，其前提是假設(shè)已發(fā)生的地質(zhì)災(zāi)害和未來(lái)將發(fā)生的災(zāi)害處于相同的地質(zhì)環(huán)境條件下。計(jì)算公式為：

(1)

式中:Pa為滑坡在因子分類(lèi)數(shù)據(jù)a中發(fā)生的條件概率，即a條件下滑坡面積與數(shù)據(jù)分類(lèi)a的面積比值比表示；Ps為整個(gè)研究區(qū)的滑坡發(fā)生概率，即滑坡面積與研究區(qū)總面積的比值(本文Ps=0.0059)。

CF的取值范圍為[-1，1]， 0～1表示滑坡發(fā)生的確定性高，該地質(zhì)環(huán)境條件下滑坡易發(fā)生，CF值越接近于1，表示該單元對(duì)滑坡易發(fā)性越敏感； 反之，- 1～0表示滑坡發(fā)生的確定性低，該地質(zhì)環(huán)境條件下滑坡不易發(fā)生； “0”表示該地質(zhì)環(huán)境下不能確定是否會(huì)發(fā)生滑坡。

1.2 邏輯回歸模型

邏輯回歸模型是二項(xiàng)分類(lèi)因變量常用的統(tǒng)計(jì)分析模型，它描述的是因變量滑坡是否發(fā)生(0代表不發(fā)生， 1代表發(fā)生)和多個(gè)致災(zāi)因子(x1，x2，x3，…xn)之間的關(guān)系。該模型中自變量可以是連續(xù)的也可以是離散的，不需要滿足正態(tài)的頻率分布。邏輯回歸函數(shù)表達(dá)式為：

P=1/[1+e-(a+β1x1+…+βixi)]

(2)

式中：P為滑坡發(fā)生概率，范圍為0～1；α為邏輯回歸計(jì)算出的一個(gè)常數(shù)項(xiàng)；β為邏輯回歸計(jì)算而得出的回歸系數(shù)；i為評(píng)價(jià)因子種類(lèi)數(shù)目。

將式(2)兩邊取自然對(duì)數(shù)，In(P/(1-P))作為因變量，將影響因子xi(i=1，2，…，n)作為自變量，得：

In(P/(1-P))=α+β1x1+β2x2…+

βnxn=α+βx

(3)

1.3 確定性系數(shù)與邏輯回歸耦合模型

鑒于確定性系數(shù)模型未能考慮各評(píng)價(jià)因子對(duì)滑坡易發(fā)性影響的差異性問(wèn)題，為提高評(píng)價(jià)精度，本文將確定性系數(shù)法與邏輯回歸模型相耦合，利用邏輯回歸法根據(jù)影響因子與歷史災(zāi)害點(diǎn)的關(guān)系統(tǒng)計(jì)得到各評(píng)價(jià)因子的權(quán)重，采用確定性系數(shù)法計(jì)算各因子不同特征變量對(duì)應(yīng)的CF值，并引入權(quán)重值將所有影響因子的CF值進(jìn)行加權(quán)求和，滑坡易發(fā)性指數(shù)I的表達(dá)式為：

(4)

式中：I為滑坡加權(quán)易發(fā)性指數(shù)；ωi為第i個(gè)評(píng)價(jià)因子的權(quán)重；Ii為第i個(gè)評(píng)價(jià)因子不同級(jí)別對(duì)應(yīng)的CF值。

2 研究區(qū)概況

九寨溝國(guó)家級(jí)自然保護(hù)區(qū)(圖1)地處青藏高原東緣的深切割高山峽谷地帶，是以高山湖泊群、瀑布群和鈣華灘流為主體的國(guó)家重點(diǎn)風(fēng)景名勝區(qū)及世界自然文化遺產(chǎn)中心，屬于高原寒溫帶-亞寒帶季風(fēng)氣候，年平均大氣降雨量為704.3mm，降水多集中在5～9月。九寨溝溝域面積為655.49 km2，為嘉陵江支流白水河的主要支流，由則查洼溝(長(zhǎng)約32 km，流域面積221.10 km2)、日則溝(長(zhǎng)約28.55 km，流域面積245.88 km2)及樹(shù)正溝(長(zhǎng)約13.4 km)構(gòu)成，3條支溝在平面上呈“Y”字型展布，形成一個(gè)完整的樹(shù)枝狀水系。地貌以高寒山地與峽谷為主，地勢(shì)南高北低，最高海拔4828 m，最低海拔2000 m，平均海拔3000 m以上，平均坡度超過(guò)30°。出露地層主要有：泥盆系、石炭系、二疊系、三疊系灰?guī)r、白云巖及第四系堆積物。

圖1 研究區(qū)位置與滑坡分布圖

研究區(qū)附近主要有岷江斷裂、塔藏?cái)嗔押突⒀罃嗔训谋蔽鞫?圖1所示“8·8”地震推測(cè)發(fā)震斷層； 李渝生等， 2017)，美國(guó)地質(zhì)調(diào)查局(USGS)公布“8·8”地震峰值加速度(PGA)為0.08～0.26 g。區(qū)內(nèi)構(gòu)造作用復(fù)雜，斷層發(fā)育(圖1)，在構(gòu)造應(yīng)力作用下，巖體節(jié)理裂隙發(fā)育，為滑坡的產(chǎn)生提供了條件。2017年九寨溝“8·8”地震誘發(fā)了大量滑坡，根據(jù)震前、震后高清遙感影像對(duì)比及現(xiàn)場(chǎng)復(fù)核，目視解譯滑坡1047處，滑坡總面積為3.88 km2。

3 九寨溝景區(qū)滑坡空間發(fā)育分布規(guī)律影響因子分析

滑坡的發(fā)生是斜坡自身內(nèi)部基礎(chǔ)地質(zhì)條件與外界環(huán)境因素共同作用的結(jié)果。前者主要包括地層巖性、地形地貌及地質(zhì)構(gòu)造等，后者為水文地質(zhì)環(huán)境、人類(lèi)工程活動(dòng)等觸發(fā)因素。基于前人的研究成果及研究區(qū)地質(zhì)條件，選取地質(zhì)構(gòu)造因子(地震峰值加速度、距斷層距離)、地形地貌因子(高程、坡度、坡向、坡面曲率)、地層巖性、其他因子(公路及水系)等分析其對(duì)滑坡發(fā)育分布的影響(圖2)?；旅婷芏仁侵富旅娣e與各級(jí)別因子總面積的比值，滑坡數(shù)密度是指單位面積內(nèi)的滑坡數(shù)量。

圖2 滑坡分布與影響因子關(guān)系

3.1 地質(zhì)構(gòu)造因子

3.1.1 地震峰值加速度

地震峰值加速度(PGA)是反映地震時(shí)地表震動(dòng)強(qiáng)度的重要參數(shù)，震中附近值較高，隨距震中距離增加而逐漸衰減。圖2a顯示滑坡主要分布在0.24 g和0.26 g 地震峰值加速度區(qū)域，占滑坡總面積的94.11%，隨著PGA由0.2 g增加至0.26 g，滑坡面密度及數(shù)密度迅速增加。

3.1.2 距斷層距離

區(qū)域斷層構(gòu)造對(duì)于地質(zhì)體中節(jié)理裂隙的發(fā)育程度起著控制性作用，將巖體切割成塊狀或碎裂狀，破壞了斜坡的穩(wěn)定性。滑坡分布與距斷層距離分布統(tǒng)計(jì)結(jié)果如圖2b所示，滑坡主要分布在距斷層2 km范圍內(nèi)，隨著距斷層距離的增加，滑坡數(shù)密度衰減迅速，表明巖體越破碎，滑坡越容易發(fā)生。

3.2 地形地貌因子

地形地貌條件是影響滑坡發(fā)育的重要因素，因此，在前人研究成果基礎(chǔ)上，選擇高程、坡度、坡向、坡面曲率等因子(由10 m×10 m的DEM數(shù)據(jù)提取獲得)作為分析滑坡發(fā)育分布的影響因素。

3.2.1 高 程

高程是坡體應(yīng)力值大小的重要影響因素，應(yīng)力會(huì)隨著坡高的增加而增加，影響著滑坡的勢(shì)能。圖2c為滑坡分布與海拔的關(guān)系，由圖2c可知滑坡主要分布在2600～3200 m范圍，該高程范圍的面積僅占研究區(qū)面積的19.67%，卻包含了69.21%的滑坡面積，其他高程范圍滑坡面積及數(shù)量顯著減少。其中高程2600～2800 m范圍面密度最高達(dá)到2.786%，數(shù)密度為6.97個(gè)·km-2，高程2800～3000 m范圍面密度次之， 3800 m以上基本沒(méi)有滑坡發(fā)生。

3.2.2 坡 度

坡度是滑坡易發(fā)性評(píng)價(jià)的重要因子，平緩的山坡由于低剪應(yīng)力使其發(fā)生滑坡概率變小，坡度越高，應(yīng)力也隨之增加，發(fā)生災(zāi)害的概率也將變大。圖2d可知，滑坡主要分布于30°～60°范圍，尤其在40°～50°達(dá)到峰值，隨著坡度增大，災(zāi)害面密度呈現(xiàn)迅速增加，當(dāng)坡度> 60°時(shí)，災(zāi)害面密度高達(dá)4.043%，數(shù)密度為8.79個(gè)·km-2，符合滑坡易發(fā)生在高坡度范圍的規(guī)律。

3.2.3 坡 向

地震滑坡坡向上的分布情況與發(fā)震斷層位置存在相關(guān)性，許強(qiáng)等(2010)將其受到地震波傳播和斷裂錯(cuò)動(dòng)的影響作用歸結(jié)為“背坡面效應(yīng)”與“斷層錯(cuò)動(dòng)方向效應(yīng)”。由圖2e可知，滑坡主要發(fā)生在北東至南坡向范圍內(nèi)，滑坡面積占滑坡總面積的73.68%。推測(cè)發(fā)震斷層(圖1)走向?yàn)楸蔽鳌蠔|，北東向至南范圍與斷層走向垂直或大角度相交且多處于背坡面，與汶川地震滑坡的分布呈現(xiàn)出相似的規(guī)律性。

3.2.4 曲 率

坡面曲率是對(duì)坡表一點(diǎn)扭曲變化程度的定量度量指標(biāo)，正值表示斜坡為凸坡，負(fù)值表示凹坡，曲率為0或接近于0表示坡面越平坦。由斜坡曲率統(tǒng)計(jì)關(guān)系圖(圖2f)可見(jiàn)，坡面曲率越大，滑坡面積及其面密度、數(shù)密度越大，雖表現(xiàn)得不顯著，但仍表明凹凸不平的斜坡較相對(duì)平坦的斜坡更有利于滑坡發(fā)生。

3.3 地質(zhì)因子

地層是滑坡發(fā)生的物質(zhì)基礎(chǔ)，不同性質(zhì)的巖石因堅(jiān)硬程度和巖體結(jié)構(gòu)的差異，直接制約滑坡發(fā)育類(lèi)型與規(guī)模，對(duì)斜坡的變形破壞起著重要的作用。基于1︰5萬(wàn)地質(zhì)圖，根據(jù)地質(zhì)年代劃分為泥盆系(D)、石炭系(C)、石炭二疊系(Cp)、二疊系(P)、三疊系(T)及第四系(Q)。統(tǒng)計(jì)結(jié)果見(jiàn)圖2g，石炭系巖組在所有地層中災(zāi)害最為發(fā)育，滑坡面積占總滑坡面積的68.6%，滑坡數(shù)密度為2.21個(gè)·km-2，第四系與石炭二疊系滑坡面密度次之，其次為二疊系、泥盆系及三疊系。

3.4 其他因子

3.4.1 公 路

景區(qū)內(nèi)觀光線路進(jìn)行大規(guī)模的切坡、開(kāi)挖等工程改變了斜坡的應(yīng)力狀態(tài)，加劇了滑坡的發(fā)生。在GIS平臺(tái)中將研究區(qū)每隔300 m建立多環(huán)緩沖區(qū)，統(tǒng)計(jì)結(jié)果見(jiàn)圖2h，滑坡面密度與數(shù)密度均隨著與公路距離的增加迅速衰減，表明人類(lèi)活動(dòng)如公路建設(shè)對(duì)滑坡的發(fā)育具有一定的影響。

3.4.2 水 系

河流水系對(duì)兩岸存在不同程度的沖刷和浸潤(rùn)現(xiàn)象，導(dǎo)致岸坡形成高陡的臨空面，從而誘發(fā)滑坡。由圖2i可知，滑坡主要分布在距水系600 m范圍內(nèi)，滑坡面密度隨距水系距離增加由0～300 m范圍處的1.19%降低至大于2100 m范圍處的0.04%，規(guī)律性較為明顯。

4 GIS支持下的滑坡易發(fā)性評(píng)價(jià)

4.1 評(píng)價(jià)因子選取與分級(jí)

根據(jù)前文中對(duì)各個(gè)因子對(duì)研究區(qū)滑坡發(fā)生的影響作用分析，共確定4類(lèi)9個(gè)指標(biāo)因子，構(gòu)成本次滑坡易發(fā)性評(píng)價(jià)的指標(biāo)體系，據(jù)前人經(jīng)驗(yàn)與統(tǒng)計(jì)結(jié)果在GIS環(huán)境下對(duì)各個(gè)因子進(jìn)行分級(jí)，結(jié)果見(jiàn)表1和圖3。

圖3 滑坡易發(fā)性評(píng)價(jià)因子分級(jí)圖

表1 因子分級(jí)和CF、WCF值

4.2 評(píng)價(jià)因子權(quán)重的確定

將是否發(fā)生滑坡作為因變量(0表示滑坡不發(fā)生， 1表示滑坡發(fā)生)，將解譯滑坡導(dǎo)入SPSS軟件進(jìn)行二元邏輯回歸分析，因子權(quán)重由大到小依次為坡度(0.2316)、坡向(0.2245)、高程(0.1872)、地層(0.1263)、距斷層距離(0.0742)、曲率(0.0414)、地震峰值加速度(0.0395)、距公路距離(0.0384)及距水系距離(0.0365)。

4.3 評(píng)價(jià)結(jié)果

根據(jù)式(1)計(jì)算各因子分類(lèi)級(jí)別的CF值，計(jì)算結(jié)果見(jiàn)表1，將各因子的CF值進(jìn)行疊加，得出基于確定性系數(shù)模型的滑坡易發(fā)性分區(qū)圖(圖4)。二元邏輯回歸模型得到的滑坡易發(fā)性分區(qū)結(jié)果見(jiàn)圖5。將CF值與權(quán)重值W代入式(4)，得到加權(quán)易發(fā)性指數(shù)WCF(表1)，在ArcGIS中將計(jì)算結(jié)果利用自然間斷法把滑坡易發(fā)性劃分為極高易發(fā)區(qū)、高易發(fā)區(qū)、中度易發(fā)區(qū)和低易發(fā)區(qū)(圖6)，分區(qū)面積占比見(jiàn)表2。其中低易發(fā)區(qū)面積為397.99 km2，占總面積的60.72%； 中度易發(fā)區(qū)面積為158.51 km2，占總面積的24.18%。高易發(fā)區(qū)和極高易發(fā)區(qū)面積分別為64.83 km2和34.17 km2，占比為9.89%和5.21%。

圖4 確定性系數(shù)模型評(píng)價(jià)結(jié)果

圖5 Logistic回歸模型評(píng)價(jià)結(jié)果

圖6 耦合模型評(píng)價(jià)結(jié)果

表2 不同評(píng)價(jià)模型下的易發(fā)性分區(qū)統(tǒng)計(jì)結(jié)果

將確定性系數(shù)模型(圖4)、邏輯回歸模型(圖5)和耦合模型(圖6)獲得的易發(fā)性評(píng)價(jià)結(jié)果進(jìn)行對(duì)比發(fā)現(xiàn)，由3種方法所計(jì)算出的滑坡易發(fā)性評(píng)價(jià)結(jié)果與滑坡點(diǎn)分布規(guī)律可知分區(qū)結(jié)果主要有以下特點(diǎn)：

(1)滑坡高-極高易發(fā)區(qū)主要沿溝谷分布，日則溝比樹(shù)正溝、則查洼溝滑坡更為發(fā)育，尤其集中在熊貓海附近，說(shuō)明地震峰值加速度對(duì)滑坡的分布有明顯的控制作用。

(2)基于耦合模型相比確定性系數(shù)與邏輯回歸模型獲得的易發(fā)性分區(qū)圖中，更少的滑坡分布在低易發(fā)區(qū)，更多的滑坡分布在高-極高易發(fā)區(qū)，具有更好的實(shí)用性。

4.4 模型精度評(píng)價(jià)及檢驗(yàn)

滑坡易發(fā)性評(píng)價(jià)結(jié)果是否準(zhǔn)確直接關(guān)系到模型的可靠性，因此有必要對(duì)評(píng)價(jià)結(jié)果進(jìn)行精度檢驗(yàn)。受試者工作特征曲線(receiver operating charac-teristic curve，簡(jiǎn)稱(chēng)ROC曲線)是滑坡易發(fā)性評(píng)價(jià)模型性能檢驗(yàn)的常用方法，以未發(fā)生滑坡的單元被正確預(yù)測(cè)的比例(假陽(yáng)性率)為橫坐標(biāo)，以發(fā)生滑坡的單元被正確預(yù)測(cè)的比例(真陽(yáng)性率)

為縱坐標(biāo)繪制曲線，曲線越靠近左上角，說(shuō)明模型分類(lèi)的準(zhǔn)確率越高。如圖7所示，AUC值為ROC曲線以下至橫坐標(biāo)的面積，是評(píng)價(jià)模型準(zhǔn)確度的指標(biāo)，取值范圍為0.5～1。由表2可以看出， 3種模型下滑坡分別有95.35%、93.37%及96.65%在中及以上易發(fā)區(qū)，仍有極少數(shù)滑坡分布在低易發(fā)區(qū)，這是由于滑坡受眾多因素影響，本文模型僅考慮了最主要影響因素。AUC值分別達(dá)到了0.832、0.813和0.847，說(shuō)明耦合模型比單一模型評(píng)價(jià)及預(yù)測(cè)能力更好。

圖7 不同評(píng)價(jià)模型下的ROC曲線

5 結(jié) 論

(1)本文以九寨溝國(guó)家級(jí)自然保護(hù)區(qū)為研究對(duì)象，通過(guò)遙感解譯和野外調(diào)查，發(fā)現(xiàn)九寨溝景區(qū)滑坡受斷層控制作用明顯，并沿溝谷及公路、水系相對(duì)發(fā)育?；贕IS平臺(tái)，采用確定性系數(shù)法、邏輯回歸及其耦合模型進(jìn)行了滑坡易發(fā)性區(qū)劃研究。

(2)基于前人研究成果及災(zāi)害發(fā)育分布規(guī)律，選取地質(zhì)構(gòu)造、地形地貌、地層巖組、其他因子等4類(lèi)共9個(gè)因子采用3種模型建立滑坡易發(fā)性評(píng)價(jià)模型，其中坡度、坡向、高程、地層巖性4個(gè)因子對(duì)景區(qū)滑坡易發(fā)性影響相對(duì)較大。

(3)根據(jù)耦合模型下滑坡易發(fā)程度分區(qū)結(jié)果，將研究區(qū)劃分為低易發(fā)區(qū)(60.72%)、中度易發(fā)區(qū)(24.18%)、高易發(fā)區(qū)(9.89%)和極高易發(fā)區(qū)(5.21%)，采用耦合模型得到的易發(fā)區(qū)劃結(jié)果與解譯及野外調(diào)查成果存在較好的一致性。結(jié)果表明采用耦合模型AUC值達(dá)到了0.847，比確定性系數(shù)模型、邏輯回歸模型等單一模型評(píng)價(jià)結(jié)果更加合理、精度更高。

(4)由于地震滑坡具有長(zhǎng)期效應(yīng)，根據(jù)易發(fā)性評(píng)價(jià)結(jié)果，熊貓海及老虎海周邊為滑坡高-極高易發(fā)區(qū)，地震產(chǎn)生的大量松散固體碎屑物質(zhì)在強(qiáng)降雨條件下極易轉(zhuǎn)化為新的滑坡或泥石流形成災(zāi)害鏈，為保障游客安全，建議針對(duì)重點(diǎn)防治區(qū)域進(jìn)行有效的工程治理、監(jiān)測(cè)及風(fēng)險(xiǎn)管控，降低成災(zāi)風(fēng)險(xiǎn)后開(kāi)放此景點(diǎn)。