李怡飛, 劉延國(guó), 梁麗萍, 李景吉

(1.成都理工大學(xué) 地球科學(xué)學(xué)院, 成都 610059; 2.西南科技大學(xué), 四川 綿陽(yáng) 621010;3.地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(成都理工大學(xué)), 成都 610059; 4.成都理工大學(xué) 生態(tài)環(huán)境學(xué)院, 成都 610059)

青藏高原崛起促成了其東南緣的高山峽谷地貌的形成。雅江縣位于青藏高原東南緣，該區(qū)域地質(zhì)活動(dòng)強(qiáng)烈，構(gòu)造復(fù)雜，并具有典型的高原深谷地貌特征[1]。縣內(nèi)峽谷兩側(cè)山勢(shì)陡峭，再加上降水、植被覆蓋度低以及強(qiáng)烈的人類工程活動(dòng)等其他因素影響，此區(qū)域是泥石流、滑坡等地質(zhì)災(zāi)害高發(fā)區(qū)。

地質(zhì)災(zāi)害對(duì)人類生產(chǎn)活動(dòng)、區(qū)域經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境都產(chǎn)生了極為重要的影響[2-3]，對(duì)地質(zhì)災(zāi)害進(jìn)行研究能夠?yàn)槿祟愐?guī)避風(fēng)險(xiǎn)、生態(tài)環(huán)境保護(hù)提供可靠的依據(jù)。目前地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)方法眾多，張曉東、李莉等[4-5]采用層次分析法，對(duì)鹽池縣、平武縣進(jìn)行研究，結(jié)果較為可靠，但層次分析法還是有一定的主觀性，需要對(duì)研究區(qū)域成災(zāi)機(jī)理等問(wèn)題進(jìn)行深入研究才能克服不足。楊樂(lè)[6]、呂遠(yuǎn)強(qiáng)等[7]分別采用灰色關(guān)聯(lián)分析法、人工神經(jīng)網(wǎng)絡(luò)法進(jìn)行地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)研究。吳柏清等[8]采用信息量方法，對(duì)九龍縣進(jìn)行研究，評(píng)價(jià)結(jié)果與實(shí)際調(diào)查成果相符合；張波等[9-14]也使用信息量法進(jìn)行地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)并得到了很好的結(jié)果，可見(jiàn)此方法易于實(shí)現(xiàn)并且結(jié)果客觀可靠。

本文以雅江縣為研究區(qū)域，運(yùn)用GIS與信息量模型相結(jié)合的方法進(jìn)行地質(zhì)災(zāi)害(泥石流、滑坡、崩塌)危險(xiǎn)性評(píng)價(jià)，一定程度上能夠揭示青藏高原東南緣區(qū)域的地質(zhì)災(zāi)害危險(xiǎn)性區(qū)劃特征，并且可以幫助當(dāng)?shù)卣玫剡M(jìn)行地質(zhì)災(zāi)害防治規(guī)劃以及區(qū)域規(guī)劃建設(shè)，減少泥石流、滑坡、崩塌等地質(zhì)災(zāi)害造成的人員和經(jīng)濟(jì)損失，為區(qū)域經(jīng)濟(jì)發(fā)展提供保障。

1 研究區(qū)概況

雅江縣(29°03′—30°30′N，100°19′—101°26′E)處于青藏高原東南緣、四川省甘孜藏族自治州南部，全縣轄區(qū)面積7 681.5 km2。該縣地處川西北丘狀高原山區(qū)和橫斷山脈中段，地勢(shì)北高南低?？h域內(nèi)西南部是極高山地貌(面積占比5%)，中部為河谷地貌(面積占比2.5%)，西北和東北部為山原地貌(面積占比92.5%)，其中大部分地區(qū)海拔3 000 m以上，有35座山峰海拔5 000 m以上。雅江縣最高點(diǎn)在東緣努其巴山峰，海拔5 252 m，最低點(diǎn)在牙衣河鄉(xiāng)布林永河口，海拔2 266 m。研究區(qū)內(nèi)雅礱江自北向南貫穿，河床海拔不足2 600 m，斜坡坡頂高程可達(dá)3 500 m，相對(duì)高差900多米，坡度多在30°以上，岸坡陡峭[15]。

2 數(shù)據(jù)來(lái)源與分析方法

2.1 數(shù)據(jù)來(lái)源

本研究所需使用的數(shù)據(jù)主要包括工程巖組、土壤類型、土地利用類型、數(shù)字高程模型(DEM)、歸一化植被指數(shù)(NDVI)、降水量等多個(gè)數(shù)據(jù)集。其中，工程巖組數(shù)據(jù)來(lái)源于中國(guó)地質(zhì)調(diào)查局1∶20萬(wàn)地質(zhì)圖；土壤類型數(shù)據(jù)使用南京土壤所1∶100萬(wàn)中國(guó)土壤數(shù)據(jù)庫(kù)(http:∥vdb3.soil.csdb.cn/)；土地利用類型基于2018年土地變更調(diào)查數(shù)據(jù)庫(kù)并結(jié)合《土地利用現(xiàn)狀分類》(GB/T21010-2017)生成；DEM以及NDVI數(shù)據(jù)是從地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/)下載；1989—2018年降水量數(shù)據(jù)下載自國(guó)家氣象信息中心(http:∥data.cma.cn/)。

2.2 地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)方法

2.2.1 評(píng)價(jià)因子選取 地質(zhì)災(zāi)害危險(xiǎn)性受多種因素影響和制約，并且不同因素對(duì)其影響和制約程度也不同[16-18]。合理地選擇評(píng)價(jià)因子才能更加準(zhǔn)確地做出地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)，每個(gè)研究區(qū)具體情況不盡相同，都具有其自身特點(diǎn)，因此評(píng)價(jià)因子也需要謹(jǐn)慎考量，本研究區(qū)評(píng)價(jià)因子的選擇不僅參考了大量已有研究成果，還結(jié)合了雅江縣當(dāng)?shù)乇阑鞯刭|(zhì)災(zāi)害的實(shí)際情況和具體特點(diǎn)。最終選擇工程巖組、土壤類型、土地利用類型、坡度、坡向、NDVI、距河流距離、距道路距離和降水量等9個(gè)因子進(jìn)行雅江縣地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)，各評(píng)價(jià)因子分級(jí)情況見(jiàn)圖1?；谲浖嗀rcGIS 10.6，由DEM數(shù)據(jù)生成雅江縣坡度和坡向數(shù)據(jù)；距河流距離和距道路距離因子分別根據(jù)雅江縣河流和道路分布圖利用緩沖區(qū)功能分析計(jì)算提?。唤邓繑?shù)據(jù)根據(jù)1989—2018年雅江縣周邊6個(gè)氣象站點(diǎn)多年平均降水量，運(yùn)用Kriging方法進(jìn)行插值得到雅江縣內(nèi)降水空間分布圖。

2.2.2 信息量法 信息量法源于信息理論，是一種對(duì)數(shù)化隨機(jī)分布概率以解決度量問(wèn)題的定量統(tǒng)計(jì)預(yù)測(cè)方法，早已被廣泛應(yīng)用于滑坡、斜坡穩(wěn)定性的空間預(yù)測(cè)等環(huán)境地質(zhì)研究中。信息量法的應(yīng)用思路可簡(jiǎn)單理解為通過(guò)計(jì)算某種地質(zhì)因素提供的研究對(duì)象的信息量值，根據(jù)其大小來(lái)評(píng)價(jià)其與地質(zhì)災(zāi)害發(fā)生的關(guān)系[8-14]。信息量法在地質(zhì)災(zāi)害評(píng)價(jià)中的使用，使得各評(píng)價(jià)因子與地質(zhì)災(zāi)害的關(guān)系緊密程度被更直觀、數(shù)字化地表達(dá)。

信息量是條件概率運(yùn)算，但在實(shí)際應(yīng)用時(shí)各評(píng)價(jià)因子對(duì)地質(zhì)災(zāi)害發(fā)生的貢獻(xiàn)率是用樣本頻率來(lái)計(jì)算：

(1)

式中:S為研究區(qū)評(píng)價(jià)單元總數(shù);N為研究區(qū)有地質(zhì)災(zāi)害分布的單元總數(shù);Si為研究區(qū)內(nèi)含有評(píng)價(jià)因素Xi的單元數(shù);Ni為分布在因素Xi內(nèi)特定類別內(nèi)的地質(zhì)災(zāi)害單元數(shù)。

單元的危險(xiǎn)性等級(jí)由總信息量值I來(lái)確定，各評(píng)價(jià)因子單元的綜合信息量值越大，表示其貢獻(xiàn)率越大，則越有利于地質(zhì)災(zāi)害的發(fā)生，意味著該區(qū)域危險(xiǎn)性等級(jí)越高，反之，貢獻(xiàn)率低不利于地質(zhì)災(zāi)害發(fā)生，其危險(xiǎn)性也就越低。需要計(jì)算出所有單元格綜合信息量值，根據(jù)其值大小進(jìn)行類別劃分，得到等級(jí)不同的危險(xiǎn)區(qū)。

圖1 地質(zhì)災(zāi)害評(píng)價(jià)因子分級(jí)

(2)

式中:I為評(píng)價(jià)單元總的信息量值;n為參評(píng)因子數(shù)。

根據(jù)前面的分析結(jié)果，選取土壤類型、坡度、坡向、工程巖組、降雨量、距道路距離、距河流距離、歸一化植被指數(shù)(NDVI)和土地利用類型9個(gè)指標(biāo)因子開(kāi)展地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)。

首先，利用ArcGIS軟件平臺(tái)，將各評(píng)價(jià)因子的基礎(chǔ)數(shù)據(jù)轉(zhuǎn)化為30 m×30 m的柵格數(shù)據(jù)，再分析各評(píng)價(jià)因子與地質(zhì)災(zāi)害點(diǎn)的空間位置關(guān)系并根據(jù)公式(1)計(jì)算每個(gè)評(píng)價(jià)因子下的每個(gè)等級(jí)信息量值，最后利用公式(2)以及疊加9個(gè)評(píng)價(jià)因子含有信息量值的圖層，依據(jù)其綜合信息量值大小重分類，得到賦有綜合信息量值的研究區(qū)地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)圖。

3 評(píng)價(jià)結(jié)果與驗(yàn)證

3.1 地質(zhì)災(zāi)害分布特征

雅江縣地質(zhì)災(zāi)害點(diǎn)基于谷歌衛(wèi)星影像(2018年)通過(guò)目視解譯獲取，并經(jīng)實(shí)地調(diào)查驗(yàn)證，地質(zhì)災(zāi)害主要類型包括滑坡、泥石流和崩塌，共計(jì)224處災(zāi)害點(diǎn)。利用GIS的空間分析功能統(tǒng)計(jì)上述評(píng)價(jià)因子與地質(zhì)災(zāi)害數(shù)量和密度的關(guān)系，見(jiàn)圖2。從工程巖組來(lái)看，地質(zhì)災(zāi)害密度最大的是硬質(zhì)巖體，約為0.03個(gè)/km2。土壤類型中，地質(zhì)災(zāi)害主要發(fā)生在半淋溶土類型，有地質(zhì)災(zāi)害點(diǎn)130個(gè)，災(zāi)害密度約為0.19個(gè)/km2。土地利用類型中，地質(zhì)災(zāi)害在水域及水利設(shè)施用地發(fā)育較多，地質(zhì)災(zāi)害密度約為0.82個(gè)/km2。坡度大于20°更易發(fā)生地質(zhì)災(zāi)害，200個(gè)地質(zhì)災(zāi)害點(diǎn)都處于20°以上坡面，其中>35°時(shí)地質(zhì)災(zāi)害密度最大，約為0.05個(gè)/km2。南坡地質(zhì)災(zāi)害更加發(fā)育，災(zāi)害密度約為0.04個(gè)/km2。距河流距離小于400 m地質(zhì)災(zāi)害最集中，200 m內(nèi)的地質(zhì)災(zāi)害密度約為0.27個(gè)/km2，200～400 m的地質(zhì)災(zāi)害密度約為0.20個(gè)/km2。距道路距離小于200 m時(shí)地質(zhì)災(zāi)害更集中，100 m內(nèi)與100～200 m地質(zhì)災(zāi)害密度均大約為0.09個(gè)/km2。NDVI在0.6～0.8地質(zhì)災(zāi)害密度最大，約為0.04個(gè)/km2，有地質(zhì)災(zāi)害點(diǎn)135處。年均降水量在68～70 mm有地質(zhì)災(zāi)害99處，地質(zhì)災(zāi)害密度約為0.05個(gè)/km2。

圖2 評(píng)價(jià)因子與地質(zhì)災(zāi)害關(guān)系

3.2 危險(xiǎn)性評(píng)價(jià)結(jié)果與驗(yàn)證

3.2.1 危險(xiǎn)性評(píng)價(jià) 根據(jù)上文信息量公式計(jì)算得出各分級(jí)指標(biāo)的信息量值，見(jiàn)表1。在ArcGIS軟件中，依據(jù)總信息量值采用自然斷點(diǎn)分級(jí)法對(duì)雅江縣進(jìn)行地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)，分為極低、低、中、高、極高危險(xiǎn)區(qū)5個(gè)分區(qū)[8]，為了能直觀地看出地質(zhì)災(zāi)害點(diǎn)分別落在何種級(jí)別危險(xiǎn)性區(qū)域，特將地質(zhì)災(zāi)害點(diǎn)與危險(xiǎn)性分區(qū)圖疊加成圖(圖3)。

表1 各因素信息量

圖3 地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)

利用ArcGIS的空間統(tǒng)計(jì)功能統(tǒng)計(jì)5個(gè)等級(jí)的危險(xiǎn)性分區(qū)面積，分區(qū)面積與雅江縣面積之比為：極低危險(xiǎn)區(qū)20.67%，低危險(xiǎn)區(qū)41.50%，主要是高山區(qū)，人口密度低，受人類工程活動(dòng)影響?。恢形ｋU(xiǎn)區(qū)23.41%，主要是在道路附近這種人類工程活動(dòng)比較強(qiáng)烈的地方；高危險(xiǎn)區(qū)9.89%，極高危險(xiǎn)區(qū)4.53%，基本位于河流兩側(cè)，主要發(fā)生泥石流災(zāi)害。

經(jīng)信息量模型分區(qū)結(jié)果和統(tǒng)計(jì)情況(表2)表明，研究區(qū)共170處地質(zhì)災(zāi)害位于高危險(xiǎn)區(qū)以及極高危險(xiǎn)區(qū)，占總地質(zhì)災(zāi)害的75.89%。其中高危險(xiǎn)區(qū)地質(zhì)災(zāi)害與總災(zāi)害比為18.30%；極高危險(xiǎn)區(qū)地質(zhì)災(zāi)害與總災(zāi)害比為57.59%。中等及以下等級(jí)危險(xiǎn)區(qū)地質(zhì)災(zāi)害發(fā)育較少。因此，最終結(jié)果可以看出地質(zhì)災(zāi)害點(diǎn)與危險(xiǎn)性分區(qū)吻合度較高，利用GIS和信息量法得出的危險(xiǎn)性分區(qū)結(jié)果比較理想。

表2 危險(xiǎn)性分區(qū)結(jié)果統(tǒng)計(jì)

3.2.2 評(píng)價(jià)精度驗(yàn)證 受訓(xùn)者操作特征曲線(receiver operating characteristics，ROC)是評(píng)價(jià)模型精準(zhǔn)度的常用方法[19-20]。它以預(yù)測(cè)結(jié)果的每個(gè)值作為可能的判斷閾值，以假陽(yáng)性率(特異度)為橫坐標(biāo)，即代表研究區(qū)中非地質(zhì)災(zāi)害實(shí)例經(jīng)試驗(yàn)診斷為陰性(不真實(shí)地質(zhì)災(zāi)害)的概率；以真陽(yáng)性率(靈敏度)為縱坐標(biāo)，即代表真實(shí)地質(zhì)災(zāi)害經(jīng)試驗(yàn)診斷為陽(yáng)性(發(fā)生地質(zhì)災(zāi)害)的概率。ROC曲線線下面積(Area Under Curve)取值范圍為[0.5,1]，中間一條斜線為參考線，其線下面積0.5，ROC曲線越靠近左上角，其線下面積越大，表示模型精度越高。

研究區(qū)內(nèi)共有地質(zhì)災(zāi)害點(diǎn)224處，為了更客觀、公正地進(jìn)行模型精度驗(yàn)證，運(yùn)用ArcGIS在研究區(qū)內(nèi)隨機(jī)生成224個(gè)測(cè)試樣本點(diǎn)，將這些測(cè)試樣本點(diǎn)以及實(shí)際地質(zhì)災(zāi)害點(diǎn)的概率數(shù)據(jù)導(dǎo)入SPSS軟件中并繪制出ROC曲線，即圖4，可以看出ROC曲線很接近左上角，AUC值為0.909，非常接近1?？梢?jiàn)在本研究區(qū)構(gòu)建的信息量模型精度較高，所得地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)結(jié)果可靠。

圖4 地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)結(jié)果ROC圖

4 討論與結(jié)論

4.1 討 論

從工程巖組來(lái)看，地質(zhì)災(zāi)害更易發(fā)生在硬質(zhì)巖體中，新鮮巖石力學(xué)強(qiáng)度相對(duì)較高，但風(fēng)化后強(qiáng)度顯著降低，常發(fā)生滑坡、崩塌等工程地質(zhì)問(wèn)題。土壤類型中，地質(zhì)災(zāi)害主要發(fā)生在半淋溶土類型，半淋溶土發(fā)生層分化明顯，表層具有不同厚度的腐殖質(zhì)層，其下游粘化土層，透水性相對(duì)較差，不利于泄洪。土地利用類型中，水域及水利設(shè)施用地地質(zhì)災(zāi)害更多，一是因?yàn)樗蚨酁榈偷貏?shì)、溝谷深壑，二是河流沖刷容易造成巖(土)體滑動(dòng)。坡度大于20°地帶更有利于地質(zhì)災(zāi)害的發(fā)育，較大的地面或溝谷坡度為其發(fā)生提供了基礎(chǔ)條件。地災(zāi)更集中在南坡，來(lái)自南方的暖濕氣流在南坡造成較為豐富的降水。距河流距離小于400 m、距道路距離小于200 m時(shí)地質(zhì)災(zāi)害更集中，河流沖刷或道路的使用使得巖(土)體更易松動(dòng)。NDVI在0.6～0.8，年均降水量在68～70 mm更易發(fā)生地質(zhì)災(zāi)害，雖然植被覆蓋度較高，但植被更喜在陽(yáng)坡(南坡)生長(zhǎng)，更多的降水也使得土壤更加松軟、巖體更易滑動(dòng)。

地質(zhì)災(zāi)害評(píng)價(jià)有定性評(píng)價(jià)與定量評(píng)價(jià)，定性評(píng)價(jià)受主觀影響，而定量評(píng)價(jià)是根據(jù)大量數(shù)據(jù)所計(jì)算得出的結(jié)果，更加客觀，因此，信息量法在地質(zhì)災(zāi)害評(píng)價(jià)中受到廣泛應(yīng)用[21-24]。本研究還存在許多不足，評(píng)價(jià)因子和模型選取的不同可能也會(huì)造成最終結(jié)果的差別，沒(méi)有運(yùn)用多個(gè)模型相互補(bǔ)充和驗(yàn)證，但是最終評(píng)價(jià)結(jié)果與青藏高原高山峽谷地貌區(qū)其他縣特征比較吻合，有一定的區(qū)域代表性。深受深切河谷地貌和人類工程活動(dòng)影響，地質(zhì)災(zāi)害主要分布在斷裂發(fā)育、巖體較破碎的河流兩側(cè)以及道路旁邊。

目前，雅江縣還沒(méi)有一個(gè)完善的地質(zhì)災(zāi)害評(píng)價(jià)系統(tǒng)，本文的危險(xiǎn)性分區(qū)為其地質(zhì)災(zāi)害評(píng)價(jià)系統(tǒng)添磚加瓦，也為以后其他人員的相關(guān)研究提供參考。青藏高原東南緣地區(qū)地質(zhì)結(jié)構(gòu)復(fù)雜，水系發(fā)達(dá)，同時(shí)還受其他多種因素影響，地質(zhì)災(zāi)害頻發(fā)，但是關(guān)于地質(zhì)災(zāi)害的研究還不是很成熟，本文也為青藏高原東南緣地區(qū)的地質(zhì)災(zāi)害危險(xiǎn)性評(píng)價(jià)提供一個(gè)參考依據(jù)。

4.2 結(jié) 論

(1) 根據(jù)各評(píng)價(jià)因子信息量值來(lái)看，研究區(qū)致災(zāi)評(píng)價(jià)因子中土地利用類型、距河流距離和距道路距離所占比重較大，是地質(zhì)災(zāi)害的主控因子。從各評(píng)價(jià)要素提供的信息量值來(lái)看，硬質(zhì)工程巖組，土壤類型為半淋溶土，土地利用類型為水域及水利設(shè)施用地，坡度>20°，坡向?yàn)槟?，NDVI值在0.6～0.8，距離河流400 m以內(nèi)，距離道路200 m以內(nèi)，降水量在68～70 mm，地質(zhì)災(zāi)害危險(xiǎn)性越高。

(2) 地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)結(jié)果表明，研究區(qū)14.42%被劃為高危險(xiǎn)區(qū)以及極高危險(xiǎn)區(qū)，研究區(qū)內(nèi)224處地質(zhì)災(zāi)害點(diǎn)，其中170處落在高危險(xiǎn)區(qū)和極高危險(xiǎn)區(qū)中，占總災(zāi)害點(diǎn)的75.89%。中危險(xiǎn)區(qū)面積占23.41%，有35處地質(zhì)災(zāi)害點(diǎn)，占總災(zāi)害點(diǎn)數(shù)15.63%。低危險(xiǎn)區(qū)與極低危險(xiǎn)區(qū)面積占研究區(qū)面積62.17%，共有災(zāi)害點(diǎn)19個(gè)，占總災(zāi)害點(diǎn)數(shù)的8.48%。AUC值為0.909，表明利用信息量法進(jìn)行研究區(qū)危險(xiǎn)性區(qū)劃較為精確。