王喜娜，黃華兵，班 亞，李 杰

(1.武漢大學(xué)測(cè)繪學(xué)院，湖北武漢430079;2.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州510275)

一、引 言

滑坡是指構(gòu)成斜坡的巖土體在重力作用下失穩(wěn)，沿著坡體內(nèi)部軟弱結(jié)構(gòu)面(帶)發(fā)生剪切而產(chǎn)生的整體性運(yùn)動(dòng)現(xiàn)象，俗稱垮山、地滑，是較為嚴(yán)重的一種山地災(zāi)害。滑坡災(zāi)害是我國(guó)山區(qū)最重要的自然災(zāi)害之一，往往給人類社會(huì)和自然環(huán)境帶來(lái)巨大的傷害和損失。因此加強(qiáng)滑坡災(zāi)害危險(xiǎn)性評(píng)估研究，以及進(jìn)行滑坡災(zāi)害的危險(xiǎn)性區(qū)劃制圖對(duì)于提高應(yīng)對(duì)滑坡災(zāi)害的防災(zāi)減災(zāi)能力具有重要的基礎(chǔ)意義［1］。

攀枝花地處四川省西南部，川滇交界處，界內(nèi)山區(qū)面積占總面積的88.38%，地貌類型復(fù)雜多樣，且是全國(guó)有名的礦區(qū)，因此加強(qiáng)地質(zhì)災(zāi)害防范，做好地質(zhì)災(zāi)害普查和危險(xiǎn)性區(qū)劃制圖是地區(qū)國(guó)土部門的重要基礎(chǔ)工作。

基于地理信息系統(tǒng)的滑坡危險(xiǎn)性區(qū)劃制圖是識(shí)別和管理潛在滑坡風(fēng)險(xiǎn)的有效工具。地理信息系統(tǒng)具有強(qiáng)大的地理信息存儲(chǔ)、處理、分析和展示功能。近20年，地理學(xué)家和工程技術(shù)人員運(yùn)用GIS識(shí)別滑坡，繪制滑坡災(zāi)害區(qū)劃圖，他們采用了多種方法，如綜合參數(shù)法(專家打分法、層次分析法)［1-3］、信息量法［4-5］、模糊綜合評(píng)判法［6］、神經(jīng)網(wǎng)絡(luò)法［7］、Logistical回歸模型［8］等方法。

筆者基于GIS技術(shù)和邏輯回歸方法，以四川省攀枝花市為研究區(qū)域，選取高程、坡度、坡向、土地利用類型、地層巖性、道路交通、河流水系7個(gè)影響因子作為評(píng)價(jià)指標(biāo)，建立邏輯回歸模型進(jìn)行預(yù)測(cè)，編制了滑坡危險(xiǎn)性區(qū)劃分布圖，圈定了滑坡災(zāi)害多發(fā)、易發(fā)區(qū)域。模型結(jié)果的預(yù)測(cè)情況顯示與實(shí)際滑坡災(zāi)害點(diǎn)在空間上的分布十分吻合。

二、研究區(qū)域概況

四川省攀枝花市位于 26°05'N—27°21'N，101°08'E—102°15'E，包括東區(qū)、西區(qū)、仁和區(qū)、米易縣、鹽邊縣5個(gè)區(qū)縣，研究區(qū)總面積7 440 km2，最高海拔4 195.5 m，最低海拔937 m。攀枝花市地處攀西裂谷中南段，屬浸蝕、剝蝕中山丘陵、山原峽谷地貌，山高谷深，盆地交錯(cuò)分布，地貌類型復(fù)雜多樣。該區(qū)域地質(zhì)災(zāi)害頻發(fā)，截至2010年，研究區(qū)內(nèi)共發(fā)現(xiàn)滑坡災(zāi)害點(diǎn)及滑坡災(zāi)害隱患點(diǎn)184處，具體分布如圖1所示。

攀枝花市的滑坡數(shù)量多，分布廣，規(guī)模大，而且發(fā)生的頻率和成災(zāi)的概率都大，給該區(qū)的工農(nóng)業(yè)生產(chǎn)﹑交通運(yùn)輸和人民生命財(cái)產(chǎn)造成了很大的危害。因此，加強(qiáng)地質(zhì)災(zāi)害防范是當(dāng)?shù)匕踩芾砉ぷ鞯闹刂兄亍?/p>

圖1 滑坡災(zāi)害分布圖

三、研究方法

1.影響因子的選擇

滑坡災(zāi)害的發(fā)育形成受內(nèi)外力綜合影響，本文根據(jù)攀枝花地區(qū)的地質(zhì)、地形調(diào)查資料，綜合考慮滑坡內(nèi)外動(dòng)力影響因素進(jìn)行滑坡影響因子的選擇。內(nèi)動(dòng)力因素包括地質(zhì)因素、地形因素的和水文因素，其中地質(zhì)因素用巖性來(lái)衡量，地形因素用坡度、坡向與高程來(lái)表達(dá)，水文因素用距水流距離來(lái)衡量。外動(dòng)力因素包括人工開挖、修坡建房等人類工程活動(dòng)，可利用距公路距離和土地利用類型來(lái)表達(dá)人類活動(dòng)對(duì)滑坡的影響。綜合上述考慮，選取地層巖性、坡度、坡向、高程、距水流距離、距公路距離和土地利用類型共7個(gè)滑坡影響因子進(jìn)行滑坡危險(xiǎn)性評(píng)價(jià)。

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

研究數(shù)據(jù)來(lái)源于中國(guó)山地災(zāi)害與環(huán)境數(shù)據(jù)庫(kù)，由中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所提供。原始數(shù)據(jù)包括攀枝花市1∶25萬(wàn)行政界線數(shù)據(jù)、攀枝花市1∶25萬(wàn)柵格地形圖數(shù)據(jù)(DEM)、攀枝花市1∶50萬(wàn)地質(zhì)類型數(shù)據(jù)、攀枝花市1∶10萬(wàn)土地利用數(shù)據(jù)(2010年)、攀枝花市1∶10萬(wàn)河流水系湖泊分布數(shù)據(jù)、攀枝花市遙感影像數(shù)據(jù)、攀枝花市道路交通數(shù)據(jù)、攀枝花市歷史滑坡災(zāi)害及災(zāi)害隱患點(diǎn)分布數(shù)據(jù)等。

3.數(shù)據(jù)預(yù)處理

首先將數(shù)據(jù)轉(zhuǎn)換到統(tǒng)一的地理坐標(biāo)系和投影坐標(biāo)系下;然后將土地利用類型和地質(zhì)巖性分布兩個(gè)矢量圖層?xùn)鸥窕?，生成的柵格圖像的分辨率與數(shù)字高程模型(DEM)的分辨率一致，即85 m×85 m;由河流線圖層和道路線圖層生成距離柵格圖，由DEM生成坡度柵格圖和坡向柵格圖;最后將所有柵格圖疊加，將發(fā)生滑坡的184個(gè)樣本點(diǎn)和未發(fā)生滑坡的148個(gè)樣本點(diǎn)共332個(gè)點(diǎn)作為采樣點(diǎn)，提取7個(gè)影響因子的屬性數(shù)據(jù)，編制 Logistic回歸數(shù)據(jù)表。圖2—圖8為各影響因子專題圖。

圖2 坡向?qū)ｎ}圖

圖3 土地利用專題圖

圖4 地層巖性專題圖

圖5 坡度專題圖

圖6 水系專題圖

圖7 道路交通專題圖

圖8 數(shù)字高程模型

4.危險(xiǎn)性評(píng)價(jià)模型

(1)理論依據(jù)

滑坡危險(xiǎn)性評(píng)價(jià)的理論依據(jù)是工程地質(zhì)類比法，即類似的工程地質(zhì)環(huán)境及其組合可能發(fā)生類似的滑坡作用?；聻?zāi)害危險(xiǎn)性評(píng)價(jià)模型經(jīng)歷了從定性模型、半定量模型到定量模型的發(fā)展。定性模型是基于專家的野外現(xiàn)場(chǎng)經(jīng)驗(yàn)判斷，直接從現(xiàn)場(chǎng)調(diào)查得出的易發(fā)性和危險(xiǎn)性結(jié)果;半定量模型雖然使用數(shù)值方法對(duì)地質(zhì)目標(biāo)的空間位置和屬性開展數(shù)值計(jì)算，并得到了量值，但評(píng)價(jià)因素權(quán)重確定仍依賴于專家經(jīng)驗(yàn);定量模型是根據(jù)過(guò)去區(qū)域滑坡災(zāi)害的分布情況，在數(shù)理統(tǒng)計(jì)理論的指導(dǎo)下，建立影響因素和地質(zhì)災(zāi)害之間的數(shù)學(xué)模型，在模型得到驗(yàn)證后，將其應(yīng)用到地質(zhì)環(huán)境相同或類似的區(qū)域，以評(píng)價(jià)地質(zhì)災(zāi)害的發(fā)生情況。定量模型包括二元與多元回歸模型、人工神經(jīng)網(wǎng)絡(luò)模型、分形理論模型、支持向量機(jī)等。邏輯回歸分析方法沒(méi)有關(guān)于分布類型、協(xié)方差陣等方面的嚴(yán)格假設(shè)，不要求數(shù)據(jù)的正態(tài)分布，其參數(shù)估計(jì)比多元判別分析更加穩(wěn)健，因此邏輯回歸方法在滑坡危險(xiǎn)性評(píng)價(jià)中有著更為廣泛的應(yīng)用。本文即選用邏輯回歸模型建立預(yù)測(cè)模型。

(2)邏輯回歸模型

邏輯回歸模型用于研究某一隨機(jī)事件發(fā)生的概率與某些因素之間的關(guān)系，如滑坡發(fā)生的概率與地層巖性、高程、河流水系分布等的關(guān)系。在邏輯回歸分析中，因變量Y是一個(gè)二分類變量，其取值Y=1和Y=0分別表示滑坡的發(fā)生和不發(fā)生。假定有p-1個(gè)因素 X1，X2，…，Xp-1影響 Y 的取值，π =P{ Y =1}為研究對(duì)象，以π(x)表示滑坡發(fā)生的概率，建立 π 與 X1，X2，…，Xp-1的 Logistic 模型

式中，f(X1，X2，…，XP-1)取為 X1，X2，…，Xp-1的線性函數(shù)，即

此時(shí)

即最終的邏輯回歸模型為

式中，β0代表回歸常數(shù);βk代表第k個(gè)影響因素的回歸系數(shù)(k=1，2，…，p-1)。

在高程、坡度、坡向、土地利用類型、地層巖性、距河流距離、距公路距離7個(gè)影響因子中，DEM、坡度圖、距河流距離圖和距道路距離圖的柵格值都屬于連續(xù)變量，在預(yù)測(cè)時(shí)需要對(duì)連續(xù)變量進(jìn)行合理的狀態(tài)劃分，這樣對(duì)優(yōu)化預(yù)測(cè)模型，提高預(yù)測(cè)效果有一定的作用。劃分區(qū)間的原則是區(qū)間之間的數(shù)據(jù)要有盡可能大的差異性，但區(qū)間數(shù)不宜太多。根據(jù)快速聚類方法劃分的影響因子分類分級(jí)，見表1。

本文選取了184個(gè)滑坡災(zāi)害點(diǎn)(含隱患點(diǎn))和148個(gè)未發(fā)生滑坡的地點(diǎn)作為樣本，共計(jì)332個(gè)樣本點(diǎn)，使用ArcGIS中Spatial Analyst工具箱Extraction工具中的Sample功能，提取樣本點(diǎn)上的各個(gè)影響因素屬性，作為試驗(yàn)數(shù)據(jù)，帶入SPSS中進(jìn)行二元邏輯回歸計(jì)算。

表1 影響因子分類分級(jí)表

5.計(jì)算結(jié)果及分析

由SPSS的邏輯回歸過(guò)程，首先得到似然比統(tǒng)計(jì)量K2的觀測(cè)值K2=68.753，自由度為7，檢驗(yàn)的p值為p=p{χ2(7)＞68.753}≤0.000 1。由此可見，x1，x2，…，x7對(duì)概率 P{ Y =1}的綜合影響是高度顯著的，建立的邏輯回歸模型是有意義的，即投入的自變量中，至少有一個(gè)自變量能有效預(yù)測(cè)因變量反映的事件發(fā)生的概率。各個(gè)影響因素的參數(shù)估計(jì)結(jié)果見表2。

表2 參數(shù)估計(jì)值及統(tǒng)計(jì)檢驗(yàn)

由于 β2、β3、β4、β5、β6的 Wald 檢驗(yàn)的 p 值較大，因此 β2、β3、β4、β5、β6對(duì) P Y=1{ }的影響不顯著，為了建立一個(gè)更緊湊有效的Logistic回歸方程，最終只保留x1、x7兩個(gè)自變量，得到如下的Logistic回歸方程

式中，x1代表高程因素;x7代表坡向因素。方程表明，攀枝花地區(qū)滑坡的發(fā)生受高程因素和坡向因素影響較大，這兩個(gè)因素對(duì)滑坡發(fā)生有較好的解釋效果。

由式(5)計(jì)算各個(gè)柵格單元的邏輯回歸預(yù)測(cè)值，并根據(jù)概率大小將回歸預(yù)測(cè)值劃分為5個(gè)等級(jí)，見表3。

表3 危險(xiǎn)性等級(jí)劃分

編制的攀枝花市滑坡災(zāi)害危險(xiǎn)性區(qū)劃圖如圖9所示，歷史滑坡點(diǎn)的頻數(shù)分布如圖10所示。由圖9可以看出，攀枝花境內(nèi)低、中、高和極高危險(xiǎn)性分級(jí)區(qū)域面積占總面積的比例分別為15.2%、23.03%、33.36%和20.98%，說(shuō)明攀枝花市約1/3地區(qū)屬于中、低等地質(zhì)災(zāi)害危險(xiǎn)區(qū)，1/2以上地區(qū)屬于地質(zhì)災(zāi)害高和極高危險(xiǎn)區(qū)。

圖9 攀枝花地區(qū)滑坡危險(xiǎn)性區(qū)劃圖

圖10 歷史滑坡點(diǎn)的頻數(shù)分布

根據(jù)圖9和圖10的區(qū)劃結(jié)果統(tǒng)計(jì)，危險(xiǎn)性級(jí)別越高的區(qū)域，滑坡面積占分級(jí)區(qū)域面積的比例就越高，即滑坡災(zāi)害密度越大。高危險(xiǎn)區(qū)和極高危險(xiǎn)區(qū)集中了大量的滑坡災(zāi)害，而低危險(xiǎn)區(qū)和中等危險(xiǎn)區(qū)的滑坡分布相對(duì)稀疏。這說(shuō)明本文劃分的危險(xiǎn)性級(jí)別區(qū)域符合歷史滑坡分布現(xiàn)狀，編制的攀枝花市滑坡危險(xiǎn)性區(qū)劃圖合理、有效。

四、結(jié)束語(yǔ)

本文基于GIS技術(shù)，利用攀枝花市基礎(chǔ)地理數(shù)據(jù)構(gòu)建了滑坡災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo)體系和評(píng)價(jià)模型，定量、客觀地完成了攀枝花市滑坡災(zāi)害危險(xiǎn)性評(píng)估和區(qū)劃，能夠較好地反映滑坡災(zāi)害危險(xiǎn)性的空間分布特征。

本研究得出的危險(xiǎn)性沒(méi)有時(shí)間因素，如果在此評(píng)價(jià)基礎(chǔ)上疊加對(duì)滑坡危險(xiǎn)性產(chǎn)生影響的時(shí)間因素，如地震因素或降雨因素，并以概率的形式與本文得到的危險(xiǎn)性相乘，則可得到更為完整的滑坡危險(xiǎn)性評(píng)價(jià)結(jié)果，也將會(huì)進(jìn)一步提高評(píng)價(jià)精度。

［1］趙寶云，胡艷香，趙忠營(yíng)，等.降雨誘發(fā)公路滑坡災(zāi)害危險(xiǎn)性評(píng)價(jià)動(dòng)態(tài)權(quán)重系統(tǒng)［J］.地球與環(huán)境，2011，39(3):376-381.

［2］肖擁軍，高華喜.庫(kù)區(qū)移民新城滑坡災(zāi)害危險(xiǎn)性區(qū)劃研究［J］.地球與環(huán)境，2011，39(3):388-392.

［3］童祥.基于ANP的鐵路隧道塌方風(fēng)險(xiǎn)評(píng)價(jià)研究［J］.隧道地下工程，2011(7):74-77.

［4］孟令超，盧曉倉(cāng)，史晨曉，等.基于信息量模型的達(dá)曲庫(kù)區(qū)滑坡危險(xiǎn)性分析［J］.災(zāi)害學(xué)，2009，24(4):31-34.

［5］朱良峰，吳信才，殷坤龍，等.基于信息量模型的中國(guó)滑坡災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究［J］.地球科學(xué)與環(huán)境學(xué)報(bào)，2004.26(3):52-56.

［6］黃俊寶.基于GIS的地質(zhì)災(zāi)害危險(xiǎn)性模糊數(shù)學(xué)評(píng)判［J］.福建地質(zhì)，2009，28(4):346-351.

［7］劉藝梁，殷坤龍，劉斌.邏輯回歸和人工神經(jīng)網(wǎng)絡(luò)模型在滑坡災(zāi)害空間預(yù)測(cè)中的應(yīng)用［J］.水文地質(zhì)工程地質(zhì)，2010，37(5):92-96.

［8］王衛(wèi)東，陳燕平，鐘晟.應(yīng)用CF和Logistic回歸模型編制滑坡危險(xiǎn)性區(qū)劃圖［J］.中南大學(xué)學(xué)報(bào):自然科學(xué)版，2009，40(4):1127-1132.

［9］石菊松，張永雙，董誠(chéng)，等.基于GIS技術(shù)的巴東新城區(qū)滑坡災(zāi)害危險(xiǎn)性區(qū)劃［J］.地球?qū)W報(bào)，2005，26(3):275-282.

［10］殷坤龍，張桂榮，陳麗霞，等.滑坡災(zāi)害風(fēng)險(xiǎn)分析［M］.北京:科學(xué)出版社，2010:30，86-88.

［11］牛全福.基于GIS的地質(zhì)災(zāi)害風(fēng)險(xiǎn)評(píng)估方法研究［D］.蘭州:蘭州大學(xué)，2011.

［12］刑秋菊，趙純勇，高克昌，等.基于GIS的滑坡危險(xiǎn)性邏輯回歸評(píng)價(jià)研究［J］.地理與地理信息科學(xué)，2004，20(3):49-51.

［13］夏瑞庚，朱元友，李新平.基于神經(jīng)網(wǎng)絡(luò)的巖質(zhì)邊坡穩(wěn)定性評(píng)估系統(tǒng)研究［J］.自然災(zāi)害學(xué)報(bào)，1996(1):98-104.

［14］范金城，梅長(zhǎng)林.數(shù)據(jù)分析［M］.北京:科學(xué)出版社，2002.

［15］劉斌，殷坤龍，劉藝梁.基于柵格數(shù)據(jù)的巴東新城區(qū)滑坡空間預(yù)測(cè)［J］.地質(zhì)科技情報(bào)，2010，29(3):133-138.