基于加權(quán)證據(jù)權(quán)模型的麗江—小金河斷裂中南段滑坡易發(fā)性評價

田鵬 常玉巧 白仙富 周青云 孫濱

摘要：麗江—小金河斷裂全新世活動強烈、地震頻發(fā)，沿斷裂帶的滑坡地質(zhì)災(zāi)害極為發(fā)育。以斷裂帶中南段兩側(cè)10km為研究區(qū)，根據(jù)地質(zhì)地理環(huán)境和滑坡發(fā)育特征，選取高程、坡度、坡向、距活動斷裂距離、距河流水系距離、距道路距離、工程地質(zhì)巖組、降雨量、土地利用類型以及地震動峰值加速度10個影響因子為評價指標，運用加權(quán)證據(jù)權(quán)模型，開展麗江—小金河斷裂中南段滑坡易發(fā)性評價，基于自然斷點法將滑坡易發(fā)程度劃分為高易發(fā)、中等易發(fā)、低易發(fā)和非易發(fā)4個級別，評價結(jié)果AUC值為0.81。結(jié)果顯示：①研究區(qū)內(nèi)滑坡受坡度、斷裂、水系、巖性因素的影響程度更高；②高易發(fā)區(qū)和中等易發(fā)區(qū)主要沿斷裂帶和金沙江等主要河流水系兩側(cè)分布，在玉龍縣、松坪鄉(xiāng)、大東鄉(xiāng)等周邊區(qū)域較集中；③西川鄉(xiāng)處于高易發(fā)區(qū)，但目前滑坡災(zāi)害點較少，應(yīng)加強關(guān)注。

關(guān)鍵詞：麗江—小金河斷裂；滑坡；加權(quán)證據(jù)權(quán)模型；易發(fā)性評價

中圖分類號：P315.94文獻標識碼：A文章編號：1000-0666（2023）03-0332-11

doi：10.20015/j.cnki.ISSN1000-0666.2023.0044

0引言

滑坡是我國最顯著的地質(zhì)災(zāi)害類型之一，它會對人類生命財產(chǎn)安全造成威脅，對重大基礎(chǔ)設(shè)施造成破壞，嚴重制約社會經(jīng)濟發(fā)展。云南地處歐亞板塊和印度板塊的碰撞接觸帶，受喜馬拉雅造山運動的影響，地質(zhì)條件復(fù)雜，滑坡地質(zhì)災(zāi)害頻發(fā)，而滇西北地區(qū)又是云南滑坡災(zāi)害最嚴重地區(qū)之一?；碌男纬蓹C制較為復(fù)雜，除與坡度、降雨等因素密切相關(guān)外，活動斷裂及地震亦是其重要影響因素。近年來，諸多學者對活動斷裂的滑坡地質(zhì)災(zāi)害開展了相關(guān)研究。如黃小龍等（2021）分析了程海斷裂古地震滑坡分布特征與斷裂活動性的關(guān)系，認為斷裂的強烈活動引發(fā)區(qū)內(nèi)強震，強震進一步誘發(fā)古地震滑坡；閆怡秋等（2021）對巴塘斷裂滑坡易發(fā)性進行了評價；楊志華等（2021）獲得該斷裂50年超越概率10%的潛在地震滑坡危險性區(qū)劃；郭長寶等（2015）研究了鮮水河斷裂帶的活動習性、地質(zhì)災(zāi)害發(fā)育特征、典型地質(zhì)災(zāi)害實例和穩(wěn)定性等，認為鮮水河斷裂對區(qū)內(nèi)地質(zhì)災(zāi)害有明顯的控制作用。前人研究成果表明，活動斷裂的展布及活動性對滑坡災(zāi)害的形成與穩(wěn)定性影響顯著。目前，對滇西北地區(qū)的滑坡地質(zhì)災(zāi)害已有豐富的研究成果，但對麗江—小金河斷裂等活動斷裂的地震高發(fā)區(qū)域滑坡易發(fā)性評價的研究仍需加強。

麗江—小金河斷裂是滇西北地區(qū)內(nèi)一條重要的活動斷裂，也是川滇菱形塊體內(nèi)重要的次級邊界斷裂（鄧起東等，2002）。沿斷裂曾發(fā)生過多次中強地震：1951年劍川6.2級地震、1976年和1998年寧蒗東北6.7級和6.4級地震、2022年寧蒗5.5級地震。斷裂周邊地區(qū)強震頻發(fā)，最大為1996年麗江7.0級地震，這次地震誘發(fā)了至少420處中小型崩塌和30處大中型滑坡，主要集中在主震震中東南部（唐川等，1997）。因此，有必要對麗江—小金河斷裂帶及其周邊地區(qū)的滑坡易發(fā)性進行評價研究。

本文整理滑坡地質(zhì)災(zāi)害調(diào)查成果，建立云南境內(nèi)麗江—小金河斷裂兩側(cè)10km區(qū)域內(nèi)的滑坡數(shù)據(jù)庫，選取高程、坡度、坡向、距活動斷裂距離等10個因子為滑坡易發(fā)程度的主控因素（白仙富等，2022；吉日伍呷等，2022；凌曉等，2021；王麗麗，2016），采用加權(quán)證據(jù)權(quán)模型，對麗江—小金河斷裂帶中南段滑坡易發(fā)性進行研究，以期為區(qū)域滑坡災(zāi)害防治提供參考。

1區(qū)域地質(zhì)背景及災(zāi)害數(shù)據(jù)

1.1地質(zhì)構(gòu)造

麗江—小金河斷裂位于川滇菱形塊體內(nèi)部，NE向延伸將川滇菱形塊體分為川西北和滇中2個次級塊體，對川滇菱形塊體區(qū)域地殼運動有重要影響（向宏發(fā)等，2002）。斷裂全長360km，傾向NW，傾角60°～80°，北起四川石棉西南，向西南經(jīng)木里、響水河、寶地、寧蒗、麗江至劍川，云南境內(nèi)長約260km（安曉文等，2018）。斷裂大致可分為3段：南段為麗江盆地以南，包含2條主要分支斷裂，其間為小型地塹或半地塹小盆地（李安等，2016）；中段從麗江盆地至木里，斷錯地貌清晰；北段為木里以北，走向由NE轉(zhuǎn)為NNE，受強烈侵蝕和切割，缺少第四系沉積（李環(huán)宇等，2020）。麗江—小金河斷裂全新世以來活動強烈，其周邊地區(qū)地震頻發(fā)，一直受諸多學者關(guān)注。徐錫偉等（2003）通過位錯地貌和年齡測試獲得斷裂中段全新世左旋滑動速率為（3.8±0.7）mm/a，垂直滑動速率為（0.65±0.14）mm/a?；贕PS資料，Shen等（2005）獲得斷裂左旋走滑速率為3mm/a。李寧等（2018）計算斷層閉鎖程度和滑動虧損速率特征，認為斷裂南段、中段具有發(fā)生較大地震的潛在危險。劉曉霞和邵志剛（2020）分析認為斷裂西南段的背景滑動速率明顯高于東北段。李安等（2016）、丁銳等（2018）通過探槽開挖，已揭露斷裂帶南段、中段全新世發(fā)生過3次古地震事件。斷裂的活動性和地震活動強度，將會加劇或誘發(fā)區(qū)域內(nèi)滑坡地質(zhì)災(zāi)害的發(fā)生。

1.2地層巖性

研究區(qū)地層發(fā)育有古近系、三疊系、二疊系、石炭系、泥盆系、奧陶系、前寒武系。斷裂帶北東區(qū)域地層發(fā)育較全，厚度較大，以淺海相沉積為主；南西區(qū)域缺失較多，厚度逐漸變薄，以濱海相為主。古近系為寶相寺組礫巖、砂巖夾粉砂巖、泥巖；三疊系為松桂組粉砂巖、頁巖；二疊系為黑泥哨組灰色含煤砂巖、頁巖夾灰?guī)r以及陽新組深色、淺色灰?guī)r，含燧石條帶或結(jié)核；石炭系主要為尖山營組上部灰?guī)r和黃龍組淺灰色厚層灰?guī)r、白云巖；志留系為青山組灰?guī)r；泥盆系為破松沖組灰?guī)r、底部泥頁巖；奧陶系為向陽組頁巖夾石英細砂巖，頂部夾含礫砂巖、細砂巖以及寶塔組灰色中厚層網(wǎng)紋灰?guī)r；前寒武系燈影組為含燧石條帶、團塊白云巖夾灰?guī)r，主要成窄條帶狀沿斷裂分布。二疊紀、三疊紀時期玄武巖大量噴發(fā)，區(qū)內(nèi)大量發(fā)育峨眉山玄武巖組地層（李環(huán)宇等，2020）。根據(jù)巖性硬度回彈值、類別及風化程度等指標，判斷研究區(qū)內(nèi)巖石主要為軟巖類，有利于滑坡發(fā)育。

1.3地質(zhì)災(zāi)害數(shù)據(jù)

筆者收集整理1∶5萬地質(zhì)災(zāi)害詳查和年度地質(zhì)災(zāi)害隱患核查、排查等數(shù)據(jù)，輔以遙感影像進行對比和修正，并選擇部分滑坡點開展野外調(diào)查工作，獲得麗江—小金河斷裂兩側(cè)10km區(qū)域內(nèi)的滑坡地質(zhì)災(zāi)害共214處，其中大型滑坡1處、中型滑坡8處，小型滑坡205處。研究區(qū)滑坡類型主要為降雨滑坡和地震滑坡，分布于河流水系兩側(cè)、斷裂帶及其交會部位（圖1）。研究區(qū)內(nèi)的滑坡點集中分布于1996年麗江7.0級地震Ⅸ、Ⅷ、Ⅶ度烈度區(qū)，Ⅵ度烈度區(qū)滑坡數(shù)量明顯減少。研究區(qū)斷裂上的許多滑坡體現(xiàn)今仍具有活動性，在多種因素誘發(fā)下極易再次活動。

2加權(quán)證據(jù)權(quán)模型

滑坡易發(fā)性評價方法分為定性和定量兩種方法。常用的定性方法有專家打分法（石菊松等，2007）和層次分析法（李學平，2001）。定量方法是根據(jù)現(xiàn)有滑坡數(shù)據(jù)和其它相關(guān)信息，運用數(shù)理統(tǒng)計模型，分析滑坡與各個影響因子的相關(guān)性。主要方法模型包括運用邏輯回歸模型（Lee，Evangelista，2006）、證據(jù)權(quán)模型（范強等，2014）、信息量模型（楊志華等，2017）等。本文采用定性和定量相結(jié)合的方法，建立集成層次分析法和加權(quán)證據(jù)權(quán)模型來評價研究區(qū)滑坡易發(fā)性。

2.1證據(jù)權(quán)模型

證據(jù)權(quán)模型是基于Bayes統(tǒng)計分析模式，通過對與滑坡形成相關(guān)的影響因素圖層的加權(quán)疊加復(fù)合獲得滑坡易發(fā)性，其優(yōu)勢為方法簡便、建模過程易懂，已被廣泛應(yīng)用于地學研究（張生元等，2009）。本文采用基于滑坡面積概率法的證據(jù)權(quán)模型，計算公式如下：

式中：w+i為當前影響因子級別內(nèi)滑坡形成的概率；w-i為當前影響因子級別外滑坡發(fā)生的概率；P為滑坡概率；B為分級后該級別下的滑坡面積；B為分級后該級別下無滑坡覆蓋面積；L為研究區(qū)總滑坡面積；L為研究區(qū)無滑坡覆蓋總面積。

根據(jù)研究區(qū)內(nèi)的滑坡點數(shù)據(jù)，分別計算每個影響因子不同級別內(nèi)的滑坡危險性權(quán)重。根據(jù)式（1）和（2）可得到：

式中：w+i表示該影響因子級別內(nèi)發(fā)生滑坡的正相關(guān)系數(shù)；w-i表示該影響因子級別內(nèi)發(fā)生滑坡的負相關(guān)系數(shù)；Npix1表示該影響因子級別內(nèi)發(fā)生滑坡的柵格數(shù)；Npix2表示該影響因子級別外發(fā)生滑坡的柵格數(shù)；Npix3表示該影響影子級別內(nèi)未發(fā)生滑坡的柵格數(shù)；Npix4表示該影響影子級別外未發(fā)生滑坡的柵格數(shù)。通過式（3）和（4）計算影響因子權(quán)重值的對數(shù)。兩權(quán)重之差是權(quán)重對比度，代表該影響因子級別內(nèi)滑坡發(fā)生的權(quán)重。

2.2層次分析法

層次分析法（AnalyticHierarchyProcess，簡稱AHP）是一種用于解決多因子決定目標的決策分析方法，對于復(fù)雜的決策思維，通過建立矩陣，每次只對兩個因子進行評判，從而提高對因子重要程度的評價準確性，再通過排序計算、一致性檢驗，獲得每個因子的權(quán)重值。在決策過程中定性和定量分析相結(jié)合，避免了人為因素導(dǎo)致經(jīng)驗判斷與實際情況不符的情況出現(xiàn)，最終的決策結(jié)果更科學合理（蒲娉璠，2016；郭長寶等，2009；劉振軍，2001）。

2.3加權(quán)證據(jù)權(quán)模型

為充分考慮不同因子對滑坡災(zāi)害發(fā)生影響程度的差異，將層次分析法與證據(jù)權(quán)模型相結(jié)合進行賦權(quán)，其公式如下：

式中：LSI為發(fā)生滑坡的加權(quán)證據(jù)權(quán)值，即為滑坡易發(fā)性評價指數(shù)；ai為第i個因子的權(quán)重值；Wf代表該影響因子級別內(nèi)滑坡發(fā)生的權(quán)重（閆怡秋等，2021）。

3滑坡易發(fā)性評價

滑坡的發(fā)生受各種地質(zhì)地理環(huán)境因素影響，其形成是一個復(fù)雜的動力學過程，滑坡在不同地區(qū)呈現(xiàn)出不同的特征。本文研究區(qū)處于滇西北地區(qū)，地貌類型復(fù)雜、地形起伏大，地勢整體從北到南逐漸降低。區(qū)域內(nèi)活動斷裂縱橫交錯，地層巖性類別多樣，河流水系發(fā)育，湖泊廣泛分布，交通道路密集（楊志華等，2021）。筆者根據(jù)研究區(qū)地質(zhì)條件、地形地貌和滑坡發(fā)育分布特征，選取10個因子作為評價指標開展滑坡易發(fā)性評價研究。

3.1評價單元確定

網(wǎng)格單元在GIS技術(shù)中能快速賦值和運算，操作便捷（劉傳正等，2007；趙曉燕，2018）。選擇合適的評價單元對滑坡易發(fā)性評價結(jié)果至關(guān)重要。本文研究區(qū)內(nèi)滑坡的影響因子數(shù)據(jù)量巨大，將研究區(qū)按25m×25m柵格大小進行劃分，獲得研究區(qū)評價單元格網(wǎng)。

3.2評價指標因子

3.2.1高程

本文使用ArcGIS軟件對研究區(qū)12.5m分辨率的DEM柵格數(shù)據(jù)做空間分析，得到研究區(qū)高程為1225～4216m，將高程分為5級：小于2000m、2000～2500m、2500～3000m、3000～3500m、大于3500m。研究區(qū)內(nèi)滑坡主要發(fā)育在高程小于3000m區(qū)域內(nèi)，占滑坡總數(shù)高達98%，3000m以上，滑坡數(shù)量明顯變少（圖3a、表1）。

3.2.2坡度

滑坡主要發(fā)生在斜坡上，坡度的大小影響斜坡上巖土體的穩(wěn)定性，同樣對斜坡體地表水的徑流、松散物質(zhì)的運移、坡體應(yīng)力分布、土地利用方式以及斜坡體的有效臨空面等都有影響。對研究區(qū)DEM數(shù)據(jù)作坡度分析，坡度主要分布在0°～81°。將坡度分為5級：0°～10°、10°～20°、20°～30°、30°～40°以及大于40°。分析結(jié)果顯示，研究區(qū)內(nèi)滑坡主要分布在坡度為10°～30°區(qū)域內(nèi)，占滑坡總數(shù)的73%（圖3b、表1），坡度大于40°的區(qū)域滑坡災(zāi)害極少發(fā)育，這可能是因為在研究區(qū)內(nèi)坡度超過40°易發(fā)生崩塌。

3.2.3坡向

由于不同的坡向接收到的太陽輻射能量不同，影響斜坡上的植被覆蓋率、巖石風化程度、降雨量和蒸發(fā)量的大小等，從而影響滑坡的發(fā)育。因此，坡向是滑坡易發(fā)性評價中非常重要的因素。將坡向分為9級，分別為平坦（無坡向）、北、北東、東、東南、南、西南、西、北西。研究區(qū)內(nèi)34%的滑坡主要分布在東至東南向的斜坡上（圖3c、表1），東和東南向斜坡日照時間長、風化作用強，容易導(dǎo)致斜坡失穩(wěn)。

3.2.4距活動斷裂距離

斷裂附近巖石較為破碎，受氣候影響易形成風化殼，從而降低坡體的完整性；斷裂的構(gòu)造活動會引發(fā)地震，這些均直接或間接影響滑坡的形成或演化。麗江—小金河斷裂是影響研究區(qū)滑坡類型和滑坡分布的重要因素。筆者對研究區(qū)斷裂兩側(cè)1000m、2000m、3000m、4000m、5000m、6000m和大于6000m作緩沖分析。從圖3d和表1可看出，研究區(qū)內(nèi)的滑坡災(zāi)害點沿斷裂延伸方向分布，其中20%的滑坡分布在距活動斷裂1000m區(qū)域內(nèi)。在5000m區(qū)域內(nèi)，隨距離增加，滑坡災(zāi)害發(fā)育減少；而大于5000m的區(qū)域，滑坡災(zāi)害仍較多，說明研究區(qū)的滑坡不只受活動斷裂影響，還受到其它地質(zhì)地理因素的綜合影響。

3.2.5距河流水系距離

河流水系對斜坡坡腳的沖蝕作用降低了坡體的穩(wěn)定性，為滑坡的發(fā)生提供了水動力條件，是影響滑坡發(fā)育的重要因素。通常河流湖泊兩側(cè)，滑坡災(zāi)害頻發(fā)。筆者對研究區(qū)河流兩側(cè)1000m、2000m、3000m、4000m、5000m以及大于5000m作緩沖分析，發(fā)現(xiàn)59%的滑坡分布在距離河流1000m區(qū)域內(nèi)（圖3e）。從表1中可看出，隨著與河流水系距離增加，滑坡分布越來越少，表明距離河流水系越近，坡腳受侵蝕作用的影響越大，斜坡失穩(wěn)的概率越高。

3.2.6距道路距離

道路建設(shè)易形成不穩(wěn)定邊坡，是滑坡災(zāi)害發(fā)育的重要誘發(fā)因素之一，因此本文選取道路為評價因子，對研究區(qū)內(nèi)的道路作空間距離緩沖分析。按照距道路500m、1000m、2000m、4000m、8000m以及大于8000m，共分為6級。結(jié)果顯示（圖3f、表1），在500m區(qū)域內(nèi)分布的滑坡數(shù)量占56%，距道路越遠，滑坡數(shù)量越少，說明道路導(dǎo)致斜坡失穩(wěn)的概率較高。

3.2.7工程地質(zhì)巖組

滑坡發(fā)育的物質(zhì)基礎(chǔ)是巖土體，不同的巖性具有不同的物理化學性質(zhì)，對滑坡的穩(wěn)定性和地表抗侵蝕能力影響顯著。據(jù)研究區(qū)1：20萬地質(zhì)圖資料，區(qū)域內(nèi)巖性多樣，根據(jù)地層巖性的工程地質(zhì)特性，將研究區(qū)工程巖組劃分為5大類：軟巖組、較軟巖組、較堅硬巖組、堅硬巖組和極堅硬巖組。結(jié)果顯示，研究區(qū)63%的滑坡發(fā)生在軟巖組內(nèi)（圖3g、表1）。

3.2.8降雨量

降雨量是滑坡地質(zhì)災(zāi)害的誘發(fā)因素之一，大氣降水滲入巖土體中，軟化巖土體，降低摩擦力，坡體容易失穩(wěn)發(fā)生滑動。研究區(qū)年平均降雨量710～951.8mm，因此將降雨量數(shù)據(jù)劃分為5個等級：小于750mm、750～800mm、800～850mm、850～900mm以及大于900mm。結(jié)果顯示，研究區(qū)內(nèi)降雨量大于900mm的區(qū)域分布的滑坡占比為87%（圖3h、表1），年均降雨量最大的區(qū)域，滑坡災(zāi)害點占比最高?；聻?zāi)害的發(fā)育與降雨量不等分級所占面積比有關(guān)，研究區(qū)內(nèi)降雨量大于900mm的區(qū)域面積最大，其滑坡災(zāi)害分布最多。降雨因素還包含雨季降雨量、累積降雨量、日最大降雨量等。研究區(qū)每年6—9月為雨季，雨季月平均降雨量為111～133mm，從北到南逐漸減小。統(tǒng)計1948—2016年研究區(qū)日最大降雨量數(shù)據(jù)作插值分析，獲得日最大降雨量為67～78mm。月最大降雨量南部多、北部少，與年降雨量、雨季月均降雨量分布規(guī)律有所差異。綜上認為，研究區(qū)滑坡分布受不同降雨特征綜合影響。

3.2.9土地利用類型

本文土地利用類型數(shù)據(jù)采用歐洲航天局（ESA）發(fā)布的2020年10m分辨率的全球土地利用（WorldCover）數(shù)據(jù)。研究區(qū)主要包括林地、灌木、草地、耕地、建筑用地、裸地（稀疏植被區(qū)）、開闊水域及草本濕地8類，滑坡主要集中分布在林地、草地、耕地和裸地（稀疏植被區(qū)），占比97%（圖3i、表1）。

3.2.10地震動峰值加速度

地震動作用于斜坡，引起斜坡變形和破壞，是誘發(fā)滑坡非常重要的一個因素。據(jù)《中國地震動峰值加速度區(qū)劃圖》（GB18306—2015），研究區(qū)地震動峰值加速度（PGA）為0.15、0.2、0.3?；轮饕性赑GA為0.2的區(qū)域內(nèi)，占滑坡總量的86%（圖3j、表1）。

3.3權(quán)重值計算

根據(jù)層次分析法計算原理，對本文中選取的10個評價指標因子建立層次分析結(jié)構(gòu)（表2）。對各評價因子建立判斷矩陣，對其兩兩比較，判斷各因子對滑坡災(zāi)害點分布的重要程度，進行標度打分，計算出各因子的權(quán)重值，并進行一致性檢驗。最終計算出矩陣的最大特征值為10.368，一致性比例CR值為0.027，小于0.1，說明判斷矩陣滿足一致性檢驗，計算所得權(quán)重具有一致性。

運用滑坡因子權(quán)重值（表2）和證據(jù)權(quán)值（表1），按照式（5），對10個評價因子的證據(jù)權(quán)值作加權(quán)計算，并進行空間疊加分析，獲得研究區(qū)滑坡易發(fā)性評價指數(shù)（LSI）為-1.25～0.87，LSI值越大，滑坡易發(fā)程度越高。

4滑坡易發(fā)性結(jié)果分析

4.1滑坡易發(fā)性結(jié)果檢驗

本文采用評價模型接收者操作特征（ROC）曲線對評價結(jié)果的合理性和準確性進行檢驗。圖4中曲線下面積（AUC）表示易發(fā)性預(yù)測成功率，其值越接近1，代表準確率越高（Hoyer，Kuss，2020；Walker，2019）。結(jié)果表明，麗江—小金河斷裂中南段滑坡易發(fā)性評價結(jié)果的AUC值為0.81，表明評價結(jié)果準確度較高。

4.2滑坡易發(fā)性評價結(jié)果

根據(jù)滑坡易發(fā)性評價指數(shù)，并按照《地質(zhì)災(zāi)害風險調(diào)查評價技術(shù)要求》（FXPC/ZRZYB-01），使用自然間斷法將研究區(qū)滑坡易發(fā)程度劃分為4個等級：高易發(fā)區(qū)（0.34≤LSI≤0.87）、中等易發(fā)區(qū)（-0.11≤LSI≤0.34）、低易發(fā)區(qū)（-0.45≤LSI≤-0.11）和非易發(fā)區(qū)（-1.25≤LSI≤-0.45），分別占研究區(qū)總面積的11%、27%、35%和27%。研究區(qū)內(nèi)的滑坡高易發(fā)區(qū)和中等易發(fā)區(qū)主要沿麗江—小金河斷裂以及區(qū)域內(nèi)金沙江等主要河流兩側(cè)分布，具體集中于玉龍縣、松坪鄉(xiāng)、大東鄉(xiāng)等周邊區(qū)域。其中十二欄桿大型滑坡以及躍進、拉八課、阿金華等中型滑坡均位于滑坡高易發(fā)區(qū)，研究區(qū)滑坡易發(fā)性評價結(jié)果與實際滑坡災(zāi)害點吻合度較好。

4.3評價結(jié)果討論

本文滑坡易發(fā)性評價結(jié)果顯示，研究區(qū)內(nèi)滑坡分布受坡度、距活動斷裂距離、距河流水系距離、工程地質(zhì)巖組因素的影響程度更高，滑坡高易發(fā)區(qū)和中等易發(fā)區(qū)主要沿活動斷裂和河流水系分布，其區(qū)域地質(zhì)環(huán)境特征為：①活動斷裂發(fā)育。研究區(qū)內(nèi)麗江—小金河斷裂和附近其它活動斷裂所在區(qū)域均為滑坡高易發(fā)區(qū)或中等易發(fā)區(qū)。距離斷裂越近，滑坡易發(fā)性越高，這是因為斷裂活動在滇西北地區(qū)易形成峽谷地貌，斷裂使巖性破碎，為滑坡發(fā)育提供有利條件。②河流水系發(fā)育。研究區(qū)內(nèi)有金沙江、漾濞江、漾弓江和寧蒗河穿過，沿河流兩側(cè)多為滑坡易發(fā)性較高區(qū)域，如松坪鄉(xiāng)西部的滑坡均沿金沙江分布，寧蒗縣附近的滑坡沿寧蒗河發(fā)育。③巖性主要為軟巖組。研究區(qū)軟巖組巖性分布面積北部大于南部，滑坡高易發(fā)區(qū)多集中于北部。如西川鄉(xiāng)滑坡高易發(fā)區(qū)面積大于玉龍縣西南區(qū)域的滑坡易發(fā)面積，雖然玉龍縣西南區(qū)域距離斷裂更近，且有漾濞江穿過，由于軟巖組巖性分布較少，滑坡易發(fā)性更弱。

研究區(qū)內(nèi)滑坡誘發(fā)因素主要為降雨量和地震。降雨型滑坡的發(fā)育與日最大降雨量、累積降雨量、年均降雨量等降雨特征密切相關(guān)。研究區(qū)年平均降雨量豐富，河流的沖蝕作用明顯，影響地形地貌坡度，為滑坡發(fā)生提供了有利條件，年均降雨量呈現(xiàn)從北東至南西逐漸減少的特征。日降雨量影響群發(fā)型滑坡的發(fā)生，一場大的降雨可能會誘發(fā)當日或后幾日發(fā)生滑坡，研究區(qū)日最大降雨量為67～78mm，南、北區(qū)域差異不明顯。累積降雨量主要對大、中型滑坡和基巖類滑坡影響更大，研究區(qū)內(nèi)雨季降雨量占年降雨量的約60%，北東區(qū)域雨季月均降雨量大于西南區(qū)域。研究區(qū)降雨型滑坡受多種降雨特征的綜合影響，滑坡分布呈現(xiàn)出北部比南部更稠密的特征。

地震型滑坡是由于地震產(chǎn)生地震動引起斜坡失穩(wěn)從而誘發(fā)滑坡，滑坡特征為規(guī)模大、地貌明顯、滑體滑移距離較遠。地震型滑坡受震源深度、發(fā)震構(gòu)造性質(zhì)、地震區(qū)地質(zhì)地貌和地形條件等因素影響。一般來說，距發(fā)震斷裂越近，滑坡密度越大，滑坡規(guī)模也越大（張永雙等，2016）。研究區(qū)內(nèi)40%的滑坡分布在麗江—小金河斷裂兩側(cè)3km區(qū)域內(nèi)，主要沿活動斷裂走向分布。斷層破碎帶風化作用強烈，多呈碎石土狀，受斷裂左旋運動或滑坡上部載體壓力作用后，易塑性流動，為滑坡的形成、發(fā)育提供了初始條件。斷裂帶上的滑坡體受地震影響后，其應(yīng)變水平會增加，地震對滑坡體變形失穩(wěn)起促進作用，滑坡體變形量增大，極易發(fā)生同震滑坡和震后滑坡。研究區(qū)地震型滑坡多為密集成群分布，發(fā)育在斜坡中上部，破壞程度高、危害性更大，易中斷交通、堰塞河道。

基于以上研究結(jié)果，研究區(qū)滑坡災(zāi)害防治工作應(yīng)注意以下幾點：若發(fā)生重大降雨事件，應(yīng)重點做好研究區(qū)北部及距離河流水系較近的滑坡高易發(fā)區(qū)和中等易發(fā)區(qū)的滑坡監(jiān)測、評估等工作；若發(fā)生地震，應(yīng)及時關(guān)注活動斷裂附近的高易發(fā)區(qū)和中等易發(fā)區(qū)的滑坡地質(zhì)災(zāi)害發(fā)展情況。

5結(jié)論

本文采用加權(quán)證據(jù)權(quán)模型對云南境內(nèi)麗江—小金河斷裂中南段兩側(cè)10km區(qū)域內(nèi)的滑坡易發(fā)性評價進行研究，主要結(jié)論如下：

（1）研究區(qū)內(nèi)共收集滑坡點214處，這些滑坡點均受麗江—小金河斷裂影響顯著，其中20%的滑坡分布在距離斷裂帶兩側(cè)1km區(qū)域內(nèi)，40%的滑坡分布在距離斷裂帶兩側(cè)3km區(qū)域內(nèi)。

（2）選取10個影響因子，采用加權(quán)證據(jù)權(quán)模型作滑坡易發(fā)性評價?；赂咭装l(fā)區(qū)、中等易發(fā)區(qū)、低易發(fā)區(qū)和非易發(fā)區(qū)分別占研究區(qū)總面積的11%、27%、35%和27%。ROC曲線檢驗?zāi)Ｐ偷腁UC=0.81，評價結(jié)果較好。本文選取的評價因子和模型方法適用于麗江—小金河斷裂中南段滑坡易發(fā)性研究。

（3）研究區(qū)內(nèi)滑坡發(fā)育主要受坡度、距活動斷裂距離、距河流水系距離、工程地質(zhì)巖組因素的影響，滑坡高易發(fā)區(qū)和中等易發(fā)區(qū)主要沿斷裂帶和金沙江等主要河流水系兩側(cè)呈條帶狀分布，在玉龍縣、松坪鄉(xiāng)、大東鄉(xiāng)等周邊區(qū)域較為集中。本文模型的計算結(jié)果與實際滑坡災(zāi)害點吻合度較高，較好地揭示了研究區(qū)滑坡易發(fā)性特征。

（4）西川鄉(xiāng)周邊區(qū)域處于高易發(fā)區(qū)和中等易發(fā)區(qū)，但實際滑坡災(zāi)害點較少，若降雨增加或發(fā)生地震，應(yīng)加強監(jiān)測此區(qū)域的滑坡災(zāi)害發(fā)展情況。

參考文獻：

安曉文，常祖峰，陳守軍，等.2018.云南第四紀活動斷裂［M］.北京：地震出版社.

白仙富，戴雨芡，葉燎原，等.2022.基于GIS和專家知識的滇西南地區(qū)滑坡敏感性模糊邏輯推理方法［J］.地震研究，45（1）：118-131.

鄧起東，張培震，冉永康，等.2002.中國活動構(gòu)造基本特征［J］.中國科學：地球科學，32（12）：1020-1030.

丁銳，任俊杰，張世民，等.2018.麗江—小金河斷裂中段晚第四紀古地震歷史［J］.地震地質(zhì)，40（3）：622-640.

范強，巨能攀，向喜瓊，等.2014.證據(jù)權(quán)法在區(qū)域滑坡危險性評價中的應(yīng)用——以貴州省為例［J］.工程地質(zhì)學報，22（3）：474-481.

郭長寶，杜宇本，張永雙，等.2015.川西鮮水河斷裂帶地質(zhì)災(zāi)害發(fā)育特征與典型滑坡形成機理［J］.地質(zhì)通報，34（1）：121-134.

郭長寶，張永雙，蔣良文，等.2009.基于GIS的滇藏鐵路麗江—香格里拉段工程地質(zhì)條件分區(qū)研究［J］.現(xiàn)代地質(zhì)，23（3）：545-552.

黃小龍，吳中海，劉鋒，等.2021.滇西北程海斷裂帶主要古地震滑坡及其分布特征的構(gòu)造解釋［J］.地學前緣，28（2）：125-139.

吉日伍呷，田宏嶺，韓繼沖，等.2022.基于不同機器學習算法的地震滑坡易發(fā)性評價——以魯?shù)榈卣馂槔跩］.昆明理工大學學報（自然科學版），17（2）：47-56.

李安，張世民，丁銳，等.2016.麗江—小金河斷裂南段全新世古地震研究［C］//中國地震局地殼應(yīng)力研究所.地殼構(gòu)造與地殼應(yīng)力論文集·27.北京：地震出版社，1-9.

李環(huán)宇，丁銳，張世民.2020.麗江—小金河斷裂南、中段地貌特征及構(gòu)造指示意義［J］.震災(zāi)防御技術(shù)，15（2）：380-391.

李寧，朱良玉，劉雷.2018.麗江—小金河斷裂帶現(xiàn)今閉鎖程度與地震危險性分析［J］.地震研究，41（2）：244-250.

李學平.2001.用層次分析法求指標權(quán)重的標度方法的探討［J］.北京郵電大學學報（社會科學版），3（1）：25-27.

凌曉，劉甲美，王濤，等.2021.基于致災(zāi)因子對稱法分級的信息量模型在地震滑坡危險性評價中的應(yīng)用［J］.國土資源遙感，33（2）：172-181.

劉傳正，劉艷輝，溫銘生，等.2007.長江三峽庫區(qū)地質(zhì)災(zāi)害成因與評價研究［M］.北京：地質(zhì)出版社.

劉曉霞，邵志剛.2020.麗江—小金河斷裂帶現(xiàn)今斷層運動特征［J］.地球物理學報，63（3）：1117-1126.

劉振軍.2001.層次分析法在青藏高原中西部航磁異常分類及找礦預(yù)測中的應(yīng)用［J］.物探與化探，25（3）：161-168.

蒲娉璠.2016.重慶市滑坡災(zāi)害時空分布特征與易發(fā)性評價研究［D］.上海：華東師范大學.

石菊松，徐瑞春，石玲，等.2007.基于RS和GIS技術(shù)的清江隔河巖庫區(qū)滑坡易發(fā)性評價與制圖［J］.地學前緣，14（6）：119-128.

唐川，黃楚興，萬曄.1997.云南省麗江大地震及其誘發(fā)的崩塌滑坡災(zāi)害特征［J］.自然災(zāi)害學報，6（3）：76-84.

王麗麗.2016.降雨型滑坡地質(zhì)災(zāi)善易發(fā)性評價中的特征處理方法［D］.杭州：浙江大學.

向宏發(fā)，徐錫偉，虢順民，等.2002.麗江—小金河斷裂第四紀以來的左旋逆推運動及其構(gòu)造地質(zhì)意義—陸內(nèi)活動地塊橫向構(gòu)造的屏蔽作用［J］.地震地質(zhì)，24（2）：190-198.

徐錫偉，聞學澤，鄭榮章，等.2003.川滇地區(qū)活動塊體最新構(gòu)造變動樣式及其動力來源［J］.中國科學：地球科學，33（S1）：151-162.

閆怡秋，楊志華，張緒教，等.2021.基于加權(quán)證據(jù)權(quán)模型的青藏高原東部巴塘斷裂帶滑坡易發(fā)性評價［J］.現(xiàn)代地質(zhì)，35（1）：26-37.

楊志華，郭長寶，吳瑞安，等.2021.青藏高原巴塘斷裂帶地震滑坡危險性預(yù)測研究［J］.水文地質(zhì)工程地質(zhì)，48（5）：91-101.

楊志華，郭長寶，姚鑫，等.2017.考慮地震后效應(yīng)的青藏高原東緣地質(zhì)災(zāi)害易發(fā)性評價［J］.中國地質(zhì)災(zāi)害與防治學報，28（4）：103-112.

張生元，成秋明，張素萍，等.2009.加權(quán)證據(jù)權(quán)模型和逐步證據(jù)權(quán)模型及其在個舊錫銅礦產(chǎn)資源預(yù)測中的應(yīng)用［J］.地球科學—中國地質(zhì)大學學報，34（2）：281-286.

張永雙，姚鑫，郭長寶，等.2016.地震擾動區(qū)重大滑坡泥石流災(zāi)害防治理論與實踐［M］.北京：科學出版社.

趙曉燕.2018.滇東北地區(qū)滑坡地質(zhì)災(zāi)害危險性評價研究［D］.昆明：云南大學.

HoyerA，KussO.2020.MetaanalysisoffullROCcurveswithflexibleparametricdistributionsofdiagnostictestvalues［J］.ResearchSynthesisMethods，11（2）：301-313.

LeeS，EvangelistaDG.2006.Earthquakeinducedlandslidesusceptibilitymappingusinganartificialneuralnetwork［J］.NaturalHazardsandEarthSystemSciences，6（5）：687-695.

ShenZK，LüJ，WangM，etal.2005.ContemporarycrustaldeformationaroundthesoutheastborderlandoftheTibetanPlateau［J］.JournalofGeophysicalResearch：SolidEarth，110（B11409）.

WalkerSP.2019.TheROCcurveredefinedoptimizingsensitivity（andspecificity）tothelivedrealityofcancer［J］.TheNewEnglandJournalofMedicine，380（17）：1594-1595.

FXPC/ZRZYB-01，地質(zhì)災(zāi)害風險調(diào)查評價技術(shù)要求［S］.

GB18306—2015，中國地震動參數(shù)區(qū)劃圖［S］.