曹 洪 洋，王 禹，滿 兵

(1.石家莊經(jīng)濟(jì)學(xué)院勘查技術(shù)與工程學(xué)院，河北 石家莊 050031；2.河北省地礦局水文工程地質(zhì)勘查院，河北 石家莊 050021)

?

基于GIS的區(qū)域群發(fā)性降雨型滑坡時(shí)空預(yù)報(bào)研究

曹 洪 洋1，王 禹1，滿 兵2

(1.石家莊經(jīng)濟(jì)學(xué)院勘查技術(shù)與工程學(xué)院，河北 石家莊 050031；2.河北省地礦局水文工程地質(zhì)勘查院，河北 石家莊 050021)

以滑坡災(zāi)害突出的雅安市雨城區(qū)為例，綜合考慮降雨強(qiáng)度、前期降雨量及下墊面(地形、巖性、植被覆蓋等)構(gòu)建了基于GIS分析獲取的易發(fā)指數(shù)+BP型神經(jīng)網(wǎng)絡(luò)時(shí)空預(yù)報(bào)模型。首先通過試驗(yàn)確定了模型的網(wǎng)絡(luò)參數(shù)和網(wǎng)絡(luò)結(jié)構(gòu)，然后通過危險(xiǎn)性區(qū)劃圖獲取降雨型滑坡易發(fā)指數(shù)，并利用GIS的空間插值功能和雨量站數(shù)據(jù)獲取相應(yīng)降雨型滑坡的雨量數(shù)據(jù)，將量化后下墊面的易發(fā)指數(shù)和降雨數(shù)據(jù)作為神經(jīng)元輸入層數(shù)據(jù)。將模型應(yīng)用于研究區(qū)，其中46個(gè)降雨型滑坡數(shù)據(jù)作為訓(xùn)練樣本，10個(gè)降雨型滑坡數(shù)據(jù)作為檢驗(yàn)樣本，預(yù)測精度達(dá)到90%，顯示該模型對(duì)于降雨型滑坡的時(shí)空預(yù)報(bào)精度較高。

降雨型滑坡；時(shí)空預(yù)報(bào)；地理信息系統(tǒng)(GIS)；BP神經(jīng)網(wǎng)絡(luò)模型

各種自然因素引發(fā)的滑坡中，降雨型滑坡發(fā)生的頻率最高，分布的地域最廣，造成的災(zāi)害最嚴(yán)重，也是滑坡預(yù)報(bào)研究的主要對(duì)象[1,2]。降雨引起的滑坡具有群發(fā)性和區(qū)域性的特點(diǎn)，一次強(qiáng)降雨可以引發(fā)幾十、數(shù)百甚至上千處滑坡。建立預(yù)報(bào)預(yù)警系統(tǒng)是防范群發(fā)性降雨型滑坡災(zāi)害的一個(gè)有效途徑，這也是目前地質(zhì)災(zāi)害的一個(gè)研究熱點(diǎn)[3]。

近幾年來龍門山斷裂活動(dòng)頻繁，造成破壞性較大的5·12汶川地震和4·20雅安地震，由其引發(fā)的次生地質(zhì)災(zāi)害也較為嚴(yán)重，大部分地質(zhì)災(zāi)害的發(fā)生與降雨有很大關(guān)系，對(duì)龍門山斷裂帶影響區(qū)域地質(zhì)災(zāi)害規(guī)律的研究也成為研究熱點(diǎn)[4,5]。本文研究區(qū)域雅安市雨城區(qū)即位于斷裂帶影響的區(qū)域。

目前對(duì)于降雨型滑坡預(yù)測預(yù)報(bào)研究主要集中在對(duì)降雨閾值的研究，而忽略了下墊面(地質(zhì)要素)和外部因素共同作用的影響。針對(duì)此問題，本文擬采用基于GIS分析技術(shù)的易發(fā)指數(shù)+BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型，充分考慮上述兩種因子，研究滑坡靜態(tài)因子、降雨與滑坡之間的響應(yīng)關(guān)系。擬將區(qū)域滑坡危險(xiǎn)性易發(fā)指數(shù)和降雨參數(shù)作為降雨型滑坡預(yù)報(bào)的輸入因子，結(jié)合研究區(qū)域歷史滑坡資料對(duì)模型的訓(xùn)練及預(yù)測，進(jìn)而達(dá)到對(duì)降雨型滑坡時(shí)空預(yù)測的目的。

1 國內(nèi)外降雨與降雨型滑坡關(guān)系研究概況

降雨型滑坡的發(fā)生既與降雨的強(qiáng)度和歷時(shí)分布有關(guān)，也與地層巖性、地質(zhì)構(gòu)造、地形地貌等環(huán)境因素有關(guān)，隨著位置不同而變化，難以事先對(duì)所有地點(diǎn)進(jìn)行治理或監(jiān)測。因此對(duì)這類滑坡災(zāi)害的防范需要通過及時(shí)預(yù)報(bào)提供警示。

降雨引發(fā)的滑坡是雨水與作為天氣系統(tǒng)下墊面的地質(zhì)、地形、土壤等因素共同作用的結(jié)果，對(duì)同一地區(qū)，如何根據(jù)給定的降雨條件，綜合利用地質(zhì)、地形、土壤等信息細(xì)化對(duì)滑坡預(yù)報(bào)的空間分辨率，既是減災(zāi)實(shí)踐的需要，也是一個(gè)重要的科學(xué)問題。地質(zhì)、地形、土壤等條件在空間上是非均勻變化的，只有在地質(zhì)、地形、土壤等因素滿足一定條件的地方，降雨才會(huì)引起滑坡。因此降雨引發(fā)滑坡的預(yù)報(bào)是一個(gè)非線性問題，它涉及復(fù)雜的系統(tǒng)工程。要實(shí)現(xiàn)對(duì)降雨引發(fā)滑坡的精細(xì)化預(yù)報(bào)，除需要提高降雨預(yù)報(bào)的時(shí)間和空間分辨率，更需要考慮下墊面因素的作用，并引入和發(fā)展新的理論、方法[3]。

國內(nèi)外現(xiàn)有的滑坡預(yù)報(bào)系統(tǒng)大多是以降雨強(qiáng)度-歷時(shí)閾值方法為基礎(chǔ)建立的，一般很少考慮地層巖性、斷裂構(gòu)造、地形地貌、植被土壤、土地利用等環(huán)境因素的影響，因此不能實(shí)時(shí)、動(dòng)態(tài)地表示出潛在滑坡的位置，難以提供有效的警示。

為解決這個(gè)問題，有學(xué)者采用事先編制的滑坡敏感度圖來指示潛在滑坡的空間分布?；旅舾卸葓D在我國稱滑坡易發(fā)程度區(qū)劃圖，對(duì)潛在滑坡的空間概率提供的是一種靜態(tài)描述，難以在滑坡敏感度圖與基于降雨閾值的警報(bào)系統(tǒng)之間實(shí)現(xiàn)實(shí)時(shí)聯(lián)動(dòng)，仍不能實(shí)時(shí)、動(dòng)態(tài)地指示潛在滑坡的位置。

國內(nèi)外學(xué)者對(duì)降雨與滑坡的關(guān)系研究主要從兩方面開展：1)研究降雨強(qiáng)度引發(fā)滑坡的閾值，以某一時(shí)段或某幾個(gè)時(shí)段的組合雨量值作為滑坡發(fā)生的一個(gè)界限值。杜榕桓等[6]研究發(fā)現(xiàn)：淺層堆積滑坡多在暴雨后10～12 h出現(xiàn)，16%的滑坡在暴雨后24 h內(nèi)出現(xiàn)，大量滑坡(包括大型和巨型滑坡)則在暴雨后28～30 h出現(xiàn)，另外77%滑坡在降雨后48 h內(nèi)產(chǎn)生，85%滑坡在降雨后72 h內(nèi)產(chǎn)生。李曉[7]對(duì)重慶一帶的地質(zhì)、地貌特點(diǎn)、降雨的侵蝕強(qiáng)度等進(jìn)行研究，分別以日降雨量25 mm、50 mm、100 mm、150 mm、200 mm和250 mm為閾值研究了上述因素之間的關(guān)系，得到了當(dāng)?shù)匕l(fā)生地表侵蝕或觸發(fā)地質(zhì)災(zāi)害有隨日降雨強(qiáng)度變化的規(guī)律。2)研究降雨時(shí)段與滑坡發(fā)生的關(guān)系，即滑坡與前期降水的關(guān)系。Lumb早在1975 年的研究中就發(fā)現(xiàn)，滑坡與前期降水有明顯的聯(lián)系[8]；但Brand 等的研究則認(rèn)為，對(duì)于大型滑坡，前期降水的多少并不是最根本的誘發(fā)因素[9]。晏同珍等[10]對(duì)災(zāi)害性滑坡的降雨臨界值進(jìn)行了討論。有的學(xué)者建議以一次降雨過程總量作為發(fā)生大量滑坡的依據(jù)，香港地區(qū)滑坡研究資料則提出以小時(shí)降雨量作為發(fā)生災(zāi)害性滑坡的臨界值。李軍等[11]的研究結(jié)論為：不同類型的滑坡，前期降水與滑坡體積的關(guān)系不盡相同，體積較大的滑坡相對(duì)于體積較小的滑坡與前期降水間有較明顯的關(guān)系。

綜合國內(nèi)外學(xué)者關(guān)于降雨與滑坡關(guān)系的研究成果，雖然在某些方面還有爭執(zhí)，但在以下方面達(dá)成了共識(shí)：1)環(huán)境因子是影響暴雨型滑坡災(zāi)害的一個(gè)重要因素；2)降雨強(qiáng)度、前期降水與滑坡災(zāi)害存在著某種關(guān)系；3)不同地區(qū)滑坡災(zāi)害存在著不同的降雨閾值，即降雨閾值與研究的范圍有關(guān)，并不存在一個(gè)適用于各個(gè)地區(qū)的降雨閾值。上述學(xué)者重點(diǎn)研究了降雨與滑坡發(fā)生之間的關(guān)系，而忽略了內(nèi)在因素和觸發(fā)因素(降雨)共同作用這一原則。降雨引發(fā)的地質(zhì)災(zāi)害是降雨作用于地面環(huán)境后，地質(zhì)-氣象兩系統(tǒng)中諸因素間的組合效應(yīng)超過某一限值導(dǎo)致的結(jié)果，所以要從這兩方面因素入手進(jìn)行研究。

2 基于BP型神經(jīng)網(wǎng)絡(luò)的降雨型滑坡預(yù)測模型的構(gòu)建

降雨型滑坡災(zāi)害空間預(yù)測的理論基礎(chǔ)是工程地質(zhì)類比法，即認(rèn)為類似的工程地質(zhì)環(huán)境可以發(fā)生類似的滑坡[12]。BP型人工神經(jīng)網(wǎng)絡(luò)是一個(gè)超大規(guī)模非線性連續(xù)時(shí)間自適應(yīng)信息處理系統(tǒng)，能夠通過學(xué)習(xí)樣本完成預(yù)測功能，符合滑坡災(zāi)害風(fēng)險(xiǎn)預(yù)測的基本原理和思路，因此筆者試圖將此模型引入降雨型滑坡預(yù)測中。

2.1 神經(jīng)網(wǎng)絡(luò)模型概述

目前得到廣泛應(yīng)用的人工神經(jīng)網(wǎng)絡(luò)(Artificial Neural Network，ANN)又稱連接機(jī)制模型(Connectionism Model)，是由大量與自然神經(jīng)細(xì)胞類似的人工神經(jīng)元廣泛互連而成的網(wǎng)絡(luò)。網(wǎng)絡(luò)的信息處理由神經(jīng)元間的相互作用實(shí)現(xiàn)，知識(shí)與信息的存儲(chǔ)表現(xiàn)為網(wǎng)絡(luò)元件互連間分布式的物理聯(lián)系，網(wǎng)絡(luò)的學(xué)習(xí)和識(shí)別決定于各神經(jīng)元連接權(quán)系的動(dòng)態(tài)演化過程[13]。

BP(Back Propagation)網(wǎng)絡(luò)是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一，訓(xùn)練機(jī)制為誤差逆?zhèn)鞑ニ惴ǎ捎米钏傧陆捣ǖ膶W(xué)習(xí)規(guī)則。通過反向傳播不斷調(diào)整其權(quán)值和閾值，達(dá)到網(wǎng)絡(luò)的誤差平方和最小。其基本原理可以從兩方面理解：1)人工神經(jīng)元。人工神經(jīng)元是人工神經(jīng)網(wǎng)絡(luò)的基本組成單位，它能完成生物神經(jīng)元最基本的3種處理過程：評(píng)價(jià)輸入信號(hào)，決定每個(gè)輸入信號(hào)的強(qiáng)度；計(jì)算所有輸入信號(hào)的權(quán)重之和，并與處理單元的閾值進(jìn)行比較：若權(quán)重和大于閾值，則人工神經(jīng)元被激發(fā)產(chǎn)生輸出信號(hào)，否則沒有輸出。2)人工神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。典型的人工神經(jīng)網(wǎng)絡(luò)模型由輸入層、中間層、輸出層組成，每一層包含許多神經(jīng)元。網(wǎng)絡(luò)通過輸入層接受輸入信號(hào)并經(jīng)過它傳送給中間層，輸入層并不處理數(shù)據(jù)。中間層也稱隱含層，其處理單元與外界沒有直接聯(lián)系，因此在人工神經(jīng)網(wǎng)絡(luò)中類似于黑箱，中間層可以為0～12層，為簡化通常設(shè)計(jì)中間層為1層或2層。網(wǎng)絡(luò)通過輸出層將結(jié)果輸出。

2.2 研究技術(shù)路線

本模型具體研究思路為：采用BP型神經(jīng)網(wǎng)絡(luò)作為數(shù)學(xué)模型，綜合考慮降雨因素和地質(zhì)因素對(duì)滑坡災(zāi)害實(shí)時(shí)預(yù)測預(yù)報(bào)等級(jí)研究。在同一地區(qū)不同的滑坡危險(xiǎn)性等級(jí)區(qū)段對(duì)同一次降雨過程的反應(yīng)是不同的，在空間預(yù)測的基礎(chǔ)上通過疊加降雨數(shù)據(jù)，從而確定在不同降雨量下的危險(xiǎn)性預(yù)測預(yù)報(bào)等級(jí)，即網(wǎng)絡(luò)輸入層參數(shù)為兩類：降雨參數(shù)和易發(fā)指數(shù)。

2.3 預(yù)測模型的構(gòu)建

基于上述技術(shù)路線，本模型首先設(shè)計(jì)了BP網(wǎng)絡(luò)配置。目前BP神經(jīng)網(wǎng)絡(luò)研究中對(duì)隱含層數(shù)目及學(xué)習(xí)步長等參數(shù)的選取還沒有理論指導(dǎo)。本研究采用從訓(xùn)練樣本中選取子樣本的方法，進(jìn)行多種網(wǎng)絡(luò)結(jié)構(gòu)的配置，從中選取較好的一種網(wǎng)絡(luò)配置。表1是網(wǎng)絡(luò)配置的參數(shù)，表2是在相同參數(shù)情況下不同隱含層個(gè)數(shù)的測試結(jié)果。

表1 網(wǎng)絡(luò)配置參數(shù)

Table 1 Neural network configuration parameters

參數(shù)激勵(lì)函數(shù)終止條件訓(xùn)練參數(shù)設(shè)置對(duì)數(shù)形式迭代次數(shù)50000動(dòng)量項(xiàng)0．2容錯(cuò)性(%)95權(quán)重調(diào)節(jié)速度0．3

表2 網(wǎng)絡(luò)結(jié)構(gòu)配置測試結(jié)果

Table 2 Testing result of network structure

實(shí)驗(yàn)編次輸入層隱含層輸出層誤差17711．921527811．871437911．8669471011．8248571111．8526

學(xué)習(xí)樣本中包含兩類數(shù)據(jù)：一是降雨數(shù)據(jù)，通過分析前人的研究，確定由臨近滑坡一次持續(xù)降雨量的強(qiáng)度和前期降雨兩個(gè)參數(shù)來研究降雨與滑坡的關(guān)系。前期降雨的時(shí)間段考慮為前1天、前5天、前10天、前15天和前1個(gè)月，包含6個(gè)降雨參數(shù)的數(shù)據(jù)分別為X1～X6。二是區(qū)域地質(zhì)因子的易發(fā)指數(shù)X7。輸出層個(gè)數(shù)為1個(gè)，對(duì)于輸出數(shù)值可按滑坡劃分為危險(xiǎn)和不危險(xiǎn)2個(gè)等級(jí)，數(shù)值化分別為1、0。

最后選用網(wǎng)絡(luò)結(jié)構(gòu)為輸入層神經(jīng)元個(gè)數(shù)7個(gè)、隱含層神經(jīng)元個(gè)數(shù)10個(gè)、輸出層神經(jīng)元個(gè)數(shù)1個(gè)。圖1為基于BP型神經(jīng)網(wǎng)絡(luò)滑坡災(zāi)害預(yù)測的網(wǎng)絡(luò)配置圖。

圖1 滑坡災(zāi)害預(yù)測模型網(wǎng)絡(luò)配置

Fig.1 Network arrangement diagram of landslide hazard prediction model

3 研究實(shí)例

3.1 研究區(qū)域概況

本文選取雅安市雨城區(qū)作為研究區(qū)，并以2003年7-9月發(fā)生在此地區(qū)56處降雨滑坡為研究對(duì)象。繼2008年5月12日四川汶川M8.0級(jí)大地震后，2013年4月20日08時(shí)02分46秒在四川蘆山發(fā)生M7.0級(jí)地震，兩次地震均發(fā)生在龍門山地震帶上[4]。與蘆山縣相鄰的雨城區(qū)同樣位于地震破壞帶內(nèi)，研究區(qū)內(nèi)松散堆積物廣布，加之氣候類型屬于亞熱帶濕潤季風(fēng)氣候，雨量充沛，降雨集中，因此本區(qū)域是降雨型滑坡的集中地。

3.2 基于GIS的降雨數(shù)據(jù)處理

為了研究降雨與滑坡之間的關(guān)系，在研究區(qū)域內(nèi)布置了20臺(tái)雨量計(jì)(圖2)。由于并不是所有的滑坡點(diǎn)都設(shè)有雨量計(jì)，為了獲得每個(gè)滑坡點(diǎn)的降雨量數(shù)據(jù)，利用GIS技術(shù)做如下處理：首先通過GIS表面插值功能，對(duì)某個(gè)時(shí)段20個(gè)雨量站點(diǎn)降雨數(shù)據(jù)進(jìn)行插值處理，形成一個(gè)連續(xù)表面，生成等值線圖。通過疊加滑坡點(diǎn)圖層和等值線圖，便可得到相應(yīng)時(shí)間段內(nèi)每個(gè)滑坡點(diǎn)的降雨量數(shù)值。分別統(tǒng)計(jì)滑坡前持續(xù)降雨量和前1天、前5天、前10天、前15天、前1個(gè)月降雨數(shù)值。

圖2 雨量計(jì)平面位置

Fig.2 The plane position plan of the hyetometers

3.3 基于GIS的易發(fā)指數(shù)數(shù)據(jù)獲取

本模型除了考慮降雨數(shù)據(jù)外還考慮了影響降雨型滑坡發(fā)生的下墊面因子，筆者在前期研究中通過GIS空間分析功能及多因子疊加方法獲取了易發(fā)指數(shù)[14,15]，易發(fā)指數(shù)中充分考慮了斜坡類型、地層巖組、高程、年降雨量、坡度、植被覆蓋、坡向、道路、水系、線性構(gòu)造等滑坡靜態(tài)影響因素。通過疊加滑坡與易發(fā)指數(shù)圖層，得到每個(gè)滑坡對(duì)應(yīng)的易發(fā)指數(shù)，此指數(shù)經(jīng)過處理作為模型的輸入值。

3.4 模型訓(xùn)練及預(yù)測

數(shù)據(jù)庫中共有232個(gè)歷史滑坡數(shù)據(jù)[14]，從中選擇與降雨有密切關(guān)系的56個(gè)滑坡作為模型的學(xué)習(xí)訓(xùn)練數(shù)據(jù)。在網(wǎng)絡(luò)模型運(yùn)算過程中，選取其中46個(gè)滑坡作為學(xué)習(xí)樣本用來訓(xùn)練網(wǎng)絡(luò)模型，當(dāng)?shù)螖?shù)大于50 000次時(shí)結(jié)束訓(xùn)練，最后控制總誤差為0.3386。將其余的10個(gè)滑坡點(diǎn)作為檢驗(yàn)樣本輸入訓(xùn)練好的模型中，如果精度達(dá)到要求，則說明網(wǎng)絡(luò)模型的建立是成功的。表3是檢驗(yàn)樣本的輸入值及學(xué)習(xí)結(jié)果，表4是檢驗(yàn)樣本的預(yù)測精度。

表3 檢驗(yàn)樣本輸入值

Table 3 Input values of the test samples

滑坡序號(hào)X1X2X3X4X5X6X7196152．3501．2190．0181．3173．367．50．3197152．3501．2190．0181．3173．367．50．3198152．3501．2190．0181．3173．367．50．2200181．3525．5212．4204．3195．882．20．2207203．9639．9260．3249．4242．2129．10．2210196．9627．2253．2242．3235．0124．80．2211210．8644．4267．3256．7249．8133．80．3212191．4639．1248．5237．1229．2120．70．39327．1282．2213．1167．197．615．30．2995．7121．261．559．422．25．70．2

表4 檢驗(yàn)樣本預(yù)測值

Table 4 Predicted values of the test samples

滑坡序號(hào)1961971982002072102112129399理想輸出值0000011000實(shí)際輸出值00000．9960．9800．9900．14300

3.5 結(jié)果分析

分析預(yù)測結(jié)果可以看到，207號(hào)滑坡預(yù)測結(jié)果與實(shí)際情況不符，211號(hào)預(yù)測結(jié)果較好。對(duì)比207號(hào)與211號(hào)滑坡的輸入數(shù)據(jù)，二者差別較小，說明其地質(zhì)條件和觸發(fā)因素比較相似，其輸出結(jié)果應(yīng)該也比較接近，考察模型的輸出值可以看到輸出值都接近于1，因此出現(xiàn)誤差的原因不是模型本身的問題。繼續(xù)考察207號(hào)李家壩滑坡，之所以在8·25降雨中沒有發(fā)生大面積的變形移動(dòng)，是因?yàn)榇嘶略谧冃纬跗谟幸恍┖喴椎闹ёo(hù)措施。模型中并沒有考慮人為的影響因素，導(dǎo)致輸出值誤差較大，需要在今后的工作中進(jìn)一步完善模型。

210號(hào)、211號(hào)和212號(hào)滑坡實(shí)際輸出值與理想輸出值有一定的差別，在實(shí)際處理中，應(yīng)將實(shí)際輸出值在0和1附近一定區(qū)間的數(shù)值都分別歸類為0和1，即分別定性為危險(xiǎn)程度小和大?？傮w看，在10個(gè)預(yù)測樣本中9個(gè)值符合要求，預(yù)測精度達(dá)到90%。

4 結(jié)論與建議

采用BP神經(jīng)網(wǎng)絡(luò)方法和GIS分析技術(shù)建立的群發(fā)性降雨型滑坡災(zāi)害預(yù)測模型，充分考慮了降雨和下墊面(地形、巖性、植被覆蓋等)兩方面因素，實(shí)現(xiàn)了滑坡區(qū)劃和降雨之間的聯(lián)動(dòng)，為降雨型滑坡的預(yù)報(bào)提供了一種新的理論方法。

本模型的預(yù)測精度會(huì)受訓(xùn)練樣本分布和降雨數(shù)據(jù)插值的影響。隨著降雨滑坡數(shù)據(jù)增多，可以加大網(wǎng)絡(luò)的訓(xùn)練樣本以提高預(yù)測精度。本模型是基于歷史滑坡的學(xué)習(xí)構(gòu)建的，因此在推廣應(yīng)用時(shí)，首先應(yīng)進(jìn)行本地區(qū)滑坡樣本的學(xué)習(xí)。

本研究所用的滑坡數(shù)據(jù)和降雨數(shù)據(jù)由中國地質(zhì)環(huán)境監(jiān)測院災(zāi)害室提供，此致謝忱！

[1] 陳偉,許強(qiáng).地質(zhì)災(zāi)害降雨預(yù)警基準(zhǔn)研究[J].地球與環(huán)境,2011,39(3):393-398.

[2] 趙衡,宋二祥.誘發(fā)區(qū)域性滑坡的降雨閾值[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2011,41(5):1481-1487.

[3] 李長江,麻士華,朱興盛,等.降雨型滑坡預(yù)報(bào)的理論、方法及應(yīng)用[M].北京:地質(zhì)出版社，2008.

[4] 蘭恒星,周成虎,高星,等.四川雅安蘆山地震災(zāi)區(qū)次生地質(zhì)災(zāi)害評(píng)估及對(duì)策建議[J].地理科學(xué)進(jìn)展,2013,32(4)：499-504.

[5] 齊信,唐川,陳州豐,等.汶川地震強(qiáng)震區(qū)地震誘發(fā)滑坡與后期降雨誘發(fā)滑坡控制因子耦合分析[J].工程地質(zhì)學(xué)報(bào),2012,20(4):522-531.

[6] 杜榕桓.長江三峽工程庫區(qū)滑坡與泥石流研究[M].成都：四川科學(xué)技術(shù)出版社,1991.

[7] 李曉.重慶地區(qū)的強(qiáng)降雨過程與地質(zhì)災(zāi)害的相關(guān)分析[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，1995,6(3):39-42.

[8] LUMB P.Slope failure in Hong Kong[J].Quarterly Journal of Engineering Geology,1975,8:31-65.

[9] BRAND E W,PREMCHITT J,PHILLIPSON H B.Relationship between rainfall and landslide in Hong Kong[A].Proceedings of the Fourth International Symposium on Landslide[C].Toronto，1984.377-384.

[10] 晏同珍,楊順安.崩塌滑坡災(zāi)害群發(fā)性初探[J].災(zāi)害學(xué),1987,2(3):23-24.

[11] 李軍,周成虎.香港地區(qū)滑坡體積與前期降水關(guān)系分析[J].自然災(zāi)害學(xué)報(bào),2002,11(2):37-45.

[12] BP模型在區(qū)域滑坡災(zāi)害風(fēng)險(xiǎn)預(yù)測中的應(yīng)用[J].巖土力學(xué),2005,26(9):1409-1413.

[13] 王偉.人工神經(jīng)網(wǎng)絡(luò)原理——入門與應(yīng)用[M].北京：北京航空航天大學(xué)出版社，1995.

[14] 曹洪洋,邊亞東.一種新的區(qū)域滑坡影響因子敏感性分析方法研究與應(yīng)用[J].湖南科技大學(xué)學(xué)報(bào),2009,24(2):49-52.

[15] 曹洪洋,陳繼華,高彥麗.基于GIS的斷層構(gòu)造對(duì)滑坡災(zāi)害發(fā)育控制作用分析——以雅安市雨城區(qū)為例[J].地球與環(huán)境,2012,40(4):595-598.

Spatio-Temporal Prediction of the Regional Group-Occurring Rainfall-Induced Landslides Based on GIS

CAO Hong-yang1,WANG Yu1,MAN Bing2

(1.CollegeofProspectingTechniques,ShijiazhuangUniversityofEconomics,Shijiazhuang050031; 2.HydrogeologyandEngineeringGeologySurveyInstitute，GeologyandMineralExplorationBureauofHebeiProvince,Shijiazhuang050021,China)

The new method of spatio-temporal prediction model is proposed based on the rainfall-induced historical landslides of Yucheng County,Yaan City,which influence factors of rain intensity,precipitation and underlying surface are fully considered.The prediction model of BP neural network and landslide easy-happening index based on GIS are constructed.Firstly,on the basis of the experiments，the network parameters and the network structure are confirmed.Secondly,easy-happening values of landslides underlying surface are obtained from the hazard zonation data.Every landslide rainfall data is performed by the spatial interpolation of GIS.The quantified underlying surface and precipitation are used as the input layer data of the model.Finally,the new method of landslide prediction is applied to the researched region.Forty-six landslides data of that are used as the training samples and ten landslides are used in the test samples.It shows 90 percent are correct.It is proved the new method of landslide prediction can predict the occurrence time of landslides with higher precision.

rainfall-induced landslides;spatio-temporal prediction;Geographic Information System(GIS);BP neural network model

2013-12-27；

2014-07-21

國家自然科學(xué)基金項(xiàng)目(41301015)；國土資源公益性行業(yè)科研專項(xiàng)項(xiàng)目(201011003)

曹洪洋(1975-)，男，博士，副教授，主要從事環(huán)境地質(zhì)和巖土工程教學(xué)和科研工作。E-mail:hongyangc@126.com

10.3969/j.issn.1672-0504.2015.01.022

P642.22

A

1672-0504(2015)01-0106-04