王佳妮, 王云琦?, 李耀明, 魏 商, 李 成, 王玉杰, 戚海梅

(1.北京林業(yè)大學(xué)水土保持學(xué)院,100083,北京;2.北京林業(yè)大學(xué)重慶縉云山三峽庫區(qū)森林生態(tài)系統(tǒng)定位觀測(cè)研究站,100083,北京;3.北京市水生態(tài)保護(hù)與水土保持中心,100086,北京;4.自然資源部重慶典型礦區(qū)生態(tài)修復(fù)野外科學(xué)觀測(cè)研究站 重慶地質(zhì)礦產(chǎn)研究院,401120,重慶)

重慶市由于其地形地貌的復(fù)雜性、極端氣候的頻發(fā)性以及三峽工程的建設(shè),致使其滑坡災(zāi)害頻發(fā),具有“點(diǎn)多、面廣”的分布特點(diǎn)[1]。2015年6月24日重慶巫山突發(fā)大面積滑坡造成6人傷亡,2001年5月1日重慶武隆縣發(fā)生山體滑坡造成79人死亡。因此,進(jìn)行重慶市的滑坡災(zāi)害易發(fā)性評(píng)價(jià)具備一定的必要性,為政府決策、提前預(yù)測(cè)、防災(zāi)減災(zāi)提供科學(xué)依據(jù)。

近年來由于遙感(remote sensing, RS)和地理信息系統(tǒng)(geographic information system, GIS)技術(shù)的迅速發(fā)展,對(duì)地質(zhì)災(zāi)害的評(píng)價(jià)由定性分析逐漸推向定量分析,大幅度提高評(píng)價(jià)和預(yù)測(cè)的精準(zhǔn)度。目前常用的定量分析模型包括:頻率比模型、信息量模型、證據(jù)權(quán)模型、加權(quán)頻率比模型、集成學(xué)習(xí)模型、隨機(jī)森林模型、人工神經(jīng)網(wǎng)絡(luò)模型、邏輯回歸模型等。毛華銳等[2]運(yùn)用頻率比法繪制渝東北三峽庫區(qū)滑坡敏感性圖;田春陽等[3]采用信息量理論開展遼寧省西豐縣的地質(zhì)災(zāi)害評(píng)價(jià)工作;郭子正等[4-5]分別采用證據(jù)權(quán)法、加權(quán)頻率比模型對(duì)三峽庫區(qū)萬州區(qū)進(jìn)行滑坡易發(fā)性等級(jí)預(yù)測(cè);周超等[6]基于集成學(xué)習(xí)和徑向基神經(jīng)網(wǎng)絡(luò)耦合模型(RBNN-Adaboost)評(píng)價(jià)三峽庫區(qū)龍駒壩滑坡敏感性;吳潤澤等[7]采用隨機(jī)森林模型對(duì)三峽庫區(qū)湖北段進(jìn)行易發(fā)性評(píng)價(jià);Yesilnacar等[8]對(duì)比分析了邏輯回歸模型和神經(jīng)網(wǎng)絡(luò)模型在土耳其Hendek地區(qū)的易發(fā)性評(píng)價(jià)結(jié)果;Sun等[9]采用貝葉斯優(yōu)化算法優(yōu)化隨機(jī)森林模型中的參數(shù),提高了三峽庫區(qū)奉節(jié)縣的滑坡評(píng)價(jià)精度;Sahana等[10]對(duì)比了邏輯回歸、模糊邏輯、頻率比3種模型在印度Rudraprayag地區(qū)的滑坡評(píng)價(jià)結(jié)果,結(jié)果表明模糊邏輯模型更適合評(píng)價(jià)Rudraprayag地區(qū)。上述模型都具各自的優(yōu)點(diǎn)和缺點(diǎn),不同地區(qū)適用的模型不同。由于信息量模型概念明確、客觀性強(qiáng)、適用范圍廣[11],使其在地質(zhì)災(zāi)害評(píng)價(jià)中被廣泛應(yīng)用,且評(píng)價(jià)結(jié)果得到認(rèn)可。

筆者以重慶市為研究對(duì)象,選取地形、地質(zhì)、水文環(huán)境、人類活動(dòng)4大類指標(biāo),采用相關(guān)性檢驗(yàn)剔除相關(guān)系數(shù)較大的指標(biāo)后,建立研究區(qū)滑坡災(zāi)害易發(fā)性的評(píng)價(jià)指標(biāo)體系;基于信息量模型確定各評(píng)價(jià)指標(biāo)對(duì)歷史滑坡災(zāi)害提供的信息量值以及確定主要影響因子;再基于GIS平臺(tái),將各評(píng)價(jià)指標(biāo)的信息量值進(jìn)行疊加得到研究區(qū)的總信息量值,進(jìn)行易發(fā)性分區(qū);最后統(tǒng)計(jì)歷史滑坡災(zāi)害點(diǎn)在各易發(fā)性分區(qū)的分布情況,計(jì)算信息量模型的預(yù)測(cè)精度。

1 研究區(qū)概況

重慶市位于我國內(nèi)陸西南部、長江上游地區(qū),下轄巫山、巫溪、奉節(jié) 、開州、云陽等38個(gè)區(qū)縣。地跨E 105°11′～110°11′、N 28°10′～32°13′,總面積約8萬2 422 km2。研究區(qū)屬于亞熱帶濕潤季風(fēng)氣候區(qū),平均相對(duì)濕度在70%～80%,年均降雨量1 000～1 350 mm,集中分布于汛期5—9月,多災(zāi)害性暴雨天氣。區(qū)內(nèi)地勢(shì)東高西低,最高海拔2 786 m,最低海拔52 m,地質(zhì)構(gòu)造復(fù)雜多樣,地形以丘陵、山地為主。由于其地形、地質(zhì)、水文環(huán)境等因子的特殊性,造成其滑坡災(zāi)害頻發(fā)。據(jù)統(tǒng)計(jì),2000—2020年期間,區(qū)內(nèi)共發(fā)生滑坡8 756次,其中主要分布在奉節(jié)(1 288次)、巫山(783次)、云陽(602次)、萬州(490次)、忠縣(467次)、江津(444次)、涪陵(431次)、開州(417次)8個(gè)區(qū)縣(圖1)。

2 評(píng)價(jià)指標(biāo)與信息量模型

2.1 評(píng)價(jià)指標(biāo)選取

根據(jù)前人研究[12],地質(zhì)災(zāi)害的影響指標(biāo)分為地形、地質(zhì)、水文環(huán)境、人類活動(dòng)因子4大類。筆者選取地形因子包括高程、坡度、坡向;地質(zhì)因子包括地層巖性、距構(gòu)造距離;水文環(huán)境因子包括距水系距離、降雨量、土壤類型;人類活動(dòng)因子包括土地利用類型、距道路距離、歸一化植被指數(shù)(normalized difference vegetation index, NDVI)。

根據(jù)野外實(shí)地調(diào)查與評(píng)價(jià)精度的要求,收集各指標(biāo)因子數(shù)據(jù),數(shù)據(jù)源見表1。采用ArcGIS 10.2中的柵格表面工具從高程數(shù)據(jù)中提取坡度、坡向;基于地質(zhì)圖、構(gòu)造綱要圖中提取巖性、構(gòu)造信息;采用ArcGIS 10.2中的多環(huán)緩沖工具對(duì)構(gòu)造、水系、道路等數(shù)據(jù)進(jìn)行處理;通過ArcGIS 10.2中的克里金插值法對(duì)多年原始降雨量數(shù)據(jù)插值得到重慶市平均月降雨柵格數(shù)據(jù);采用ArcGIS 10.2中地圖代數(shù)的柵格計(jì)算器能求得重慶市2000—2020年的NDVI多年均值。在基于遙感影像數(shù)據(jù)確定全區(qū)滑坡分布的情況下,根據(jù)野外實(shí)地調(diào)查、滑坡詳勘資料等對(duì)滑坡歷史清單信息進(jìn)行更新補(bǔ)充。

2.2 信息量模型

滑坡災(zāi)害的發(fā)生受多種因子耦合作用影響,區(qū)域位置不同,各因子的影響程度不同。信息量法是一種統(tǒng)計(jì)預(yù)測(cè)方法,通過統(tǒng)計(jì)各因子在歷史滑坡災(zāi)害中所提供的信息量,對(duì)區(qū)域內(nèi)滑坡易發(fā)性進(jìn)行分區(qū)具有依據(jù)性。信息量值越大,指標(biāo)因子對(duì)滑坡的影響程度越大,在進(jìn)行分區(qū)時(shí)滑坡災(zāi)害發(fā)生的可能性越大。

式中:F(i,L)為指標(biāo)因子(或范圍)i對(duì)區(qū)域內(nèi)滑坡災(zāi)害發(fā)生所提供的信息量值;Ni為指標(biāo)因子(或范圍)i內(nèi)已發(fā)生的滑坡次數(shù);N為研究區(qū)已發(fā)生的滑坡總次數(shù);Mi為指標(biāo)因子(或范圍)i的面積,km2;M為研究區(qū)總面積,km2。

完成各指標(biāo)因子(或范圍)i的信息量計(jì)算后,通過ArcGIS 10.2中的柵格疊加功能將所有指標(biāo)因子(或范圍)i的信息量值進(jìn)行疊加得到研究區(qū)的總信息量值

式中n為指標(biāo)因子數(shù)量。

3 指標(biāo)因子分析

3.1 地形因子

高程、坡度屬于連續(xù)型變量,依據(jù)自然斷點(diǎn)法將高程、坡度指標(biāo)分為6個(gè)等級(jí),滑坡災(zāi)害及各指標(biāo)的等級(jí)空間分布見圖2,對(duì)應(yīng)信息量值見表2。高程是影響滑坡發(fā)育的重要因素,具體體現(xiàn)在:直接影響坡體內(nèi)部的受力面分布以及受力是否均勻,以及間接影響溫差、植被分布、人類活動(dòng)等指標(biāo)。重慶市高程分布呈現(xiàn)東高西低,90.50%的滑坡發(fā)生在[3 m, 957 m)區(qū)間內(nèi)(圖2a)。當(dāng)高程位于[957 m,1 294 m)、[1 294 m,1 739 m)、[1 739 m,2 786 m]區(qū)間時(shí),滑坡數(shù)量分別為704、117、8,滑坡數(shù)量極劇減少,原因?yàn)楫?dāng)高程大于957 m時(shí),降水和人類活動(dòng)較少,導(dǎo)致滑坡不易發(fā)生。坡度會(huì)直接影響斜坡的穩(wěn)定性,對(duì)滑坡災(zāi)害的發(fā)生起著決定作用[13]。由圖2b可知,重慶市坡度分布為東北高、西北低,81.06%的滑坡發(fā)生在[7°, 28°)之間,當(dāng)坡度≥37°時(shí),幾乎沒有滑坡發(fā)生(1.18%)。坡向是指臨空面的朝向,不同坡向直接影響光熱的再分配,影響水分狀況和植被分布,從而影響斜坡穩(wěn)定性。坡向?qū)儆陔x散型變量,本文等級(jí)劃分中,坡向分為平地、北、東北、東、東南、南、西南、西、西北9個(gè)等級(jí)?；略谄降厣戏植驾^少(滑坡次數(shù)為3),南、西南、東南的信息值較大,分別為0.100、0.092、0.036,表明滑坡更易發(fā)生在南面[14]。

圖2 重慶市地形因子Fig.2 Topographic factors of Chongqing city

3.2 地質(zhì)因子

地層巖性屬于離散型變量,直接將類型劃分為12個(gè)等級(jí)。重慶市主要以砂巖、石灰?guī)r、泥巖、頁巖為主,其次為泥灰?guī)r、礫巖、花崗巖,局部分布有千枚巖、板巖、玄武巖(圖3a)。不同類型的巖石物理力學(xué)特性存在較大差異[15],例如泥巖、頁巖、板巖、千枚巖的巖性差,所構(gòu)成的斜坡易風(fēng)化產(chǎn)生泥化礦物,致使坡體結(jié)構(gòu)面抗剪強(qiáng)度降低,破壞坡體穩(wěn)定性。重慶市的砂巖、泥巖主要分布在西北部,石灰?guī)r、頁巖主要分布在東北部和南部。其中,91.53%的滑坡發(fā)生在砂巖、泥巖、頁巖、石灰?guī)r4個(gè)巖層上,但板巖的信息量值最大,表明重慶市板巖分布少,但只要板巖存在的區(qū)域發(fā)生滑坡的概率最大。距構(gòu)造距離屬于連續(xù)型變量,基于ArcGIS的多環(huán)緩沖功能將距離分為7個(gè)等級(jí),地質(zhì)構(gòu)造通過影響坡體內(nèi)部軟弱結(jié)構(gòu)面的發(fā)展與演化間接影響斜坡穩(wěn)定。據(jù)圖3b和表2可知,48.20%的滑坡發(fā)生在距構(gòu)造2 500 m的范圍之外,但其信息量值最小。距構(gòu)造距離[500 m,1 000 m)區(qū)間的信息量值最大,表明在距構(gòu)造[500 m, 1 000 m)范圍內(nèi)最易發(fā)生滑坡。

圖3 重慶市地質(zhì)因子Fig.3 Geological factors of Chongqing city

3.3 水文環(huán)境因子

距水系距離屬于連續(xù)型變量,基于ArcGIS的多環(huán)緩沖功能將距水系距離分為7個(gè)等級(jí)。水流的下切作用以及入滲作用會(huì)破壞邊坡內(nèi)部應(yīng)力分布,據(jù)圖4a和表2可知重慶市58.92%的滑坡發(fā)生距水系2 500 m范圍之外,其信息量值最小。距水系距離[200 m, 500 m)的信息量值最大,表明重慶市在距水系[200 m, 500 m)距離內(nèi)最易發(fā)生滑坡。降雨量屬于連續(xù)型變量,依據(jù)自然斷點(diǎn)法將降雨量分為6個(gè)等級(jí);降雨是滑坡災(zāi)害產(chǎn)生的主要誘發(fā)因素[16],其影響主要體現(xiàn)在4方面:1)降雨對(duì)地表的擊濺作用破壞坡面表層結(jié)構(gòu);2)降雨入滲會(huì)在坡面淺表層形成暫態(tài)飽和區(qū),降低土體有效應(yīng)力,加速坡體變形;3)雨水的入滲提高坡體含水量,增加坡體自重;4)未入滲的雨水形成地表徑流沖刷坡面。據(jù)圖4b和表2可知,重慶市的月均降雨量主要分布于[134 mm, 157 mm]之間,東南、西北降雨量較大;其中,6 526次(74.53%)滑坡發(fā)生在降雨量[134 mm, 148mm)的區(qū)域,[134 mm, 139mm)區(qū)間的信息量值最大,表明在月均降雨量位于[134 mm, 139mm)區(qū)間的地域,最易發(fā)生滑坡災(zāi)害。

圖4 重慶市水文環(huán)境因子Fig.4 Hydro-environmental factors of Chongqing city

土壤類型屬于離散型變量,直接將類型劃分為8個(gè)等級(jí);不同類型的土壤存在不同的物理力學(xué)特性。據(jù)圖4c和表2可知,重慶市主要以黃壤和紫色土為主,分別占47.94%和39.36%,黃壤土主要分布在西北部,紫色土主要分布在東北部和南部。89.47%的滑坡發(fā)生在黃壤和紫色土,黃壤的信息量值最大,表明重慶市滑坡災(zāi)害主要發(fā)生在黃壤地區(qū)。

3.4 人類活動(dòng)因子

土地利用類型屬于離散型變量,直接將類型劃分為6個(gè)等級(jí)。土地利用會(huì)影響水土流失、降雨入滲、地表結(jié)構(gòu)特征,不合理的土地利用加速滑坡發(fā)育。據(jù)圖5a和表2可知,重慶市的土地利用主要以耕地和林地為主,分別占45.48%和40.76%;耕地主要分布在西北部和西南部,林地主要分布在東北和東南部。其中,84.43%的滑坡發(fā)生在耕地和林地上,兩者的信息量值分別為0.159和-0.270,耕地更易誘發(fā)滑坡發(fā)生。距道路距離屬于連續(xù)型變量,基于ArcGIS的多環(huán)緩沖功能將距道路距離分為7個(gè)等級(jí)。道路建設(shè)破壞了原有的自然地質(zhì)條件,破壞坡體自身平衡狀態(tài)。據(jù)圖5b和表2可知,50.54%的滑坡發(fā)生在距道路2 500 m范圍之外,但其信息量值最小;[0, 200 m)的信息量值最大,表明在距道路200 m范圍之內(nèi)最易發(fā)生滑坡。NDVI屬于連續(xù)型變量,依據(jù)自然斷點(diǎn)法將NDVI分為6個(gè)等級(jí),植被的護(hù)坡作用通過力學(xué)機(jī)制、水文機(jī)制發(fā)揮作用。重慶市的NDVI分布在[8%, 74%]之間,主要呈現(xiàn)東北高、西北低的走勢(shì)(圖5c)。NDVI 為[8%, 38%)時(shí),信息量值最大,表明植被覆蓋較少的地方更易發(fā)生滑坡災(zāi)害。

圖5 重慶市人類活動(dòng)因子Fig.5 Human activity factors of Chongqing city

4 重慶市滑坡易發(fā)性評(píng)價(jià)

4.1 指標(biāo)因子相關(guān)性

依據(jù)工程地質(zhì)類比法的思想,各指標(biāo)因子之間存在一定的相關(guān)性。若選取所有指標(biāo)進(jìn)行評(píng)價(jià)分析,各指標(biāo)之間的影響權(quán)重可能會(huì)疊加,影響評(píng)價(jià)結(jié)果的準(zhǔn)確性。因此,采用ArcGIS中的波段集統(tǒng)計(jì)工具對(duì)11個(gè)指標(biāo)進(jìn)行相關(guān)性檢驗(yàn),得到各指標(biāo)之間的相關(guān)性系數(shù)分布在-0.310～0.648之間,以弱相關(guān)和無相關(guān)為主(圖6)。其中,共有3組指標(biāo)因素的相關(guān)性較強(qiáng):高程與坡度(0.440)、高程與NDVI(0.648)、坡度與NDVI(0.510)。由于高程分別和坡度、NDVI指標(biāo)相關(guān),且相關(guān)系數(shù)較大,因此易發(fā)性評(píng)價(jià)過程中剔除高程指標(biāo)。

4.2 基于信息量模型的指標(biāo)因子分析

通過直接分類、自然斷點(diǎn)、多環(huán)緩沖3種方法確定各指標(biāo)因子的等級(jí)劃分后,結(jié)合相關(guān)性分析結(jié)果,選取坡度、坡向、地層巖性、距構(gòu)造距離、距水系距離、降雨量、土壤類型、土地利用類型、距道路距離、NDVI 10個(gè)指標(biāo)作為信息量模型的評(píng)價(jià)指標(biāo)體系。將各指標(biāo)因子的信息量值進(jìn)行統(tǒng)計(jì)并進(jìn)行排序(表2),排序前5位的分別為:距水系距離[200 m,500 m)(信息量值0.809)、NDVI[8%,38%)(信息量值0.563)、坡度[13°,20°)(信息量值0.500)、距水系距離[0,200 m)(信息量值0.429)、距水系距離[500 m,1 000 m)(信息量值0.428)。表明重慶市水系、NDVI、坡度是影響滑坡發(fā)育的主要影響因子。

4.3 重慶市滑坡災(zāi)害易發(fā)性評(píng)價(jià)

將評(píng)價(jià)指標(biāo)體系中的10個(gè)指標(biāo)因子的信息量通過ArcGIS的柵格疊加分析得到重慶市的總信息量值分布于-5.516～3.584之間。通過自然斷點(diǎn)法將總信息量劃分為5類:不易發(fā)區(qū)[-5.516,-2.233)、低易發(fā)區(qū)[-2.233,-1.087)、中易發(fā)區(qū)[-1.087,-0.239)、高易發(fā)區(qū)[-0.239,0.573)、極高易發(fā)區(qū)[0.573,3.584]。

繪制滑坡災(zāi)害易發(fā)性分區(qū)圖,研究區(qū)主要以高易發(fā)區(qū)和中易發(fā)區(qū)為主(圖7)。重慶市極高易發(fā)區(qū)面積為1萬4 230.14 km2,占研究區(qū)面積的17.26%,已發(fā)生滑坡次數(shù)3 627次,占滑坡總數(shù)41.42%。極高易發(fā)區(qū)在重慶市各個(gè)區(qū)縣均有分布,主要集中分布于渝中、南岸、大渡口、巴南、江北、云陽。根據(jù)圖7可以看出,極高易發(fā)區(qū)的空間分布與水系分布高度吻合。高易發(fā)區(qū)面積為2萬6 222.76 km2,占研究區(qū)面積的31.82%,已發(fā)生滑坡次數(shù)3 220次,占滑坡總數(shù)36.77%。高易發(fā)區(qū)在重慶市各個(gè)區(qū)縣均有分布,主要集中分布于酉陽、忠縣、九龍坡、彭水、渝北區(qū)域。中易發(fā)區(qū)面積為2萬3 260.00 km2,占研究區(qū)面積的28.22%,已發(fā)生滑坡次數(shù)1 351次,占滑坡總數(shù)15.43%。低易發(fā)區(qū)面積為1萬5 293.23 km2,占研究區(qū)面積的18.55%,已發(fā)生滑坡次數(shù)501次,占滑坡總數(shù)5.72%。不易發(fā)區(qū)面積為3 415.88 km2,占研究區(qū)面積的4.14%,已發(fā)生滑坡次數(shù)57次,占滑坡總數(shù)0.65%。其中,渝中、南岸、巴南、大渡口、江北、云陽分別有96.97%、86.46%、80.94%、78.53%、77.46%和71.56%的區(qū)域處于極高易發(fā)區(qū)和高易發(fā)區(qū),應(yīng)加強(qiáng)滑坡災(zāi)害防治措施。

圖7 重慶市滑坡災(zāi)害易發(fā)性分區(qū)Fig.7 Landslide susceptibility mapping of Chongqing city

5 討論

基于信息量模型確定重慶市滑坡災(zāi)害的3大主要影響指標(biāo)為水系、NDVI和坡度。水系對(duì)邊坡影響主要體現(xiàn)在3方面[17]:1)河流對(duì)岸坡坡腳的沖刷、掏蝕作用,為滑坡的發(fā)生提供了良好的臨空條件;2)水位的周期性變化,使岸坡土體長期處于干濕交替狀態(tài),破壞岸坡土體強(qiáng)度,使岸坡變形特征產(chǎn)生不可逆的變化;3)水的入滲作用會(huì)增加邊坡含水量,對(duì)巖土體起到軟化作用。當(dāng)距水系距離特別近時(shí),水同時(shí)也會(huì)對(duì)邊坡土體產(chǎn)生靜水壓力,促使邊坡短時(shí)間內(nèi)不會(huì)發(fā)生失穩(wěn),這可能也是距水系[200 m,500 m)信息量值比[0,200 m)大的原因,是否正確還需進(jìn)一步細(xì)化研究。植被對(duì)邊坡穩(wěn)定存在正負(fù)兩方面作用[18],正作用包括:截留雨水、減弱地表沖刷、根系錨固、調(diào)節(jié)水文;負(fù)作用主要包括:增加坡體自重、增加土壤滲透系數(shù)、向坡體轉(zhuǎn)移風(fēng)荷載[19]。根據(jù)表2,NDVI 5個(gè)等級(jí)對(duì)應(yīng)的信息量值分別為0.563、-0.143、-0.126、0.319、-0.059和-0.552,表明當(dāng)NDVI 為[8%,38%)和[54%,59%)時(shí),植被的負(fù)作用占主要地位,促進(jìn)滑坡的發(fā)生;當(dāng)NDVI處于[38%,54%)和[63%,74%]時(shí),植被的正作用占據(jù)主要地位。坡度決定坡體內(nèi)部應(yīng)力分布狀態(tài)、地表徑流狀態(tài)、地下水補(bǔ)給情況,以及決定由重力所產(chǎn)生的下滑力大小和滑坡規(guī)模,坡度越大斜坡穩(wěn)定性越差[13]。根據(jù)表2,坡度6個(gè)等級(jí)對(duì)應(yīng)的信息量值分別為-0.879、0.082、0.500、0.269、-0.300和-1.164,信息量值隨坡度的升高逐漸增加,在[13°,20°)區(qū)間達(dá)到最大值,隨后隨坡度增加逐漸減小。后面信息量值逐漸減小的原因可能是,當(dāng)坡度≥20°時(shí),人類工程設(shè)施(修建公路、房屋修建、水利工程)減少、雨水入滲量減少,滑坡災(zāi)害不易發(fā)生。

此外,基于遙感和GIS平臺(tái),采用信息量模型確定各指標(biāo)因子為滑坡災(zāi)害發(fā)生提供的信息量值,進(jìn)行重慶市滑坡災(zāi)害易發(fā)性分區(qū),是一項(xiàng)有利于防治重慶市地質(zhì)災(zāi)害的工作。在易發(fā)性評(píng)價(jià)結(jié)果中,易發(fā)性越高的地區(qū)滑坡密度越大,則代表評(píng)價(jià)結(jié)果越合理。在本次評(píng)價(jià)結(jié)果中,存在78.20%滑坡分布于極高易發(fā)區(qū)和高易發(fā)區(qū),且滑坡密度隨易發(fā)性分區(qū)等級(jí)的增加而增加,表明本次研究重慶市滑坡易發(fā)性分區(qū)的結(jié)果有效且合理(表3)。采用歷史滑坡點(diǎn)落在極高易發(fā)區(qū)、高易發(fā)區(qū)的比例表征信息量模型預(yù)測(cè)重慶市滑坡災(zāi)害的精度,本研究評(píng)價(jià)精度為78.20%。

表3 重慶市滑坡易發(fā)性分區(qū)結(jié)果Tab.3 Results of landslide susceptibility zoning

6 結(jié)論

1)以重慶市為研究對(duì)象,選取高程、坡度、坡向、地層巖性、距構(gòu)造距離、距水系距離、降雨量、土壤類型、土地利用類型、距道路距離、NDVI共11個(gè)指標(biāo)作為研究滑坡的影響因子。通過相關(guān)性分析剔除高程因子后,構(gòu)建重慶市滑坡災(zāi)害易發(fā)性評(píng)價(jià)指標(biāo)體系?；?000—2020年的滑坡災(zāi)害數(shù)據(jù),采用信息量模型繪制重慶市滑坡災(zāi)害易發(fā)性分區(qū)圖,采用自然斷點(diǎn)法將研究區(qū)劃分為極高易發(fā)區(qū)、高易發(fā)區(qū)、中易發(fā)區(qū)、低易發(fā)區(qū)、不易發(fā)區(qū)。研究區(qū)內(nèi)主要以高易發(fā)區(qū)、中易發(fā)區(qū)為主,信息量模型評(píng)價(jià)精度達(dá)78.20%,能夠較好反映重慶市的滑坡災(zāi)害情況。

2)根據(jù)信息量模型確定各指標(biāo)因子的信息量值,影響重慶滑坡災(zāi)害的主要指標(biāo)為水系、NDVI和坡度。其中,在距水系距離[200 m,500 m)范圍、NDVI[8%,38%)、坡度[13°,20°)的區(qū)域,發(fā)生滑坡災(zāi)害的概率最大。重慶市的渝中、南岸、巴南、大渡口、江北、云陽地區(qū)均有超過70%的區(qū)域處于極高易發(fā)區(qū)和高易發(fā)區(qū),應(yīng)加強(qiáng)滑坡災(zāi)害的防治工作。