周明浪，邵新民，羅美芳

(1.溫州市地質(zhì)環(huán)境監(jiān)測站，浙江溫州 325027;2.浙江省地質(zhì)環(huán)境監(jiān)測院，浙江杭州 310007)

0 引言

進入21世紀以來，滑坡地質(zhì)災害預警預報工作在國內(nèi)外發(fā)展很快。預警方法主要有面狀預警和點狀預警兩種，目前國內(nèi)外研究的現(xiàn)狀是以面狀預警為主，但是，預警單元面積大，尤其在中國，預警單元常是上萬平方千米。近年來，為提高精度，采用預警單元較小的“點狀”預警工作得到加強。點狀預警目的是提高預警準確度和水平，為政府在地質(zhì)災害防治管理方面提供技術(shù)支撐。

溫州市山區(qū)面積大，臺風暴雨頻繁，山體穩(wěn)定性差，工程建設強度大，突發(fā)性地質(zhì)災害頻發(fā)，是浙江省地質(zhì)災害最嚴重的市，災害點總數(shù)占全省四分之一以上。眾多隱患點使得溫州成為浙江省地質(zhì)災害防治任務最艱巨、防災形勢最為嚴峻的地區(qū)。針對嚴重的防災形勢，溫州市政府決定開展小范圍的地質(zhì)災害預警預報工作，即滑坡“點狀”地質(zhì)災害預警預報?；境霭l(fā)點是采用不同于省級層面的面狀預警方式，縮小預警范圍，提高預警精度，即:將每個行政村作為一個預警單元，這樣能改變傳統(tǒng)盲目的人員大轉(zhuǎn)移為有目的人員撤離，將極大的提高溫州市防預地質(zhì)災害的能力。

1 研究現(xiàn)狀

2003年，中國地質(zhì)環(huán)境監(jiān)測院將大陸劃分為7個大區(qū)、28個預警區(qū);2004年，進一步將大陸共分為74個預警區(qū)域［1］，分別建立了每個預警區(qū)的臨界降雨判據(jù)，開展全國地質(zhì)災害預警工作。這兩年，中國劃分預警單元面積達十幾或幾十平方千米，預警面積過大。2003年4月7日，中國國土資源部和中國氣象局簽訂《關(guān)于聯(lián)合開展地質(zhì)災害氣象預報預警工作協(xié)議》，同年6月1日，全國地質(zhì)災害氣象預報預警工作正式啟動。浙江省從2003年開始，每年汛期(4月15日～10月31日)由浙江省國土資源廳和浙江省氣象局共同開展突發(fā)性地質(zhì)災害氣象預(警)報工作。據(jù)此，溫州市地質(zhì)災害預警預報工作由此起步。

國家和省級地質(zhì)災害預警預報系統(tǒng)均是宏觀的面狀預警，通常預警面積過大(圖1、圖2)。溫州市市級地質(zhì)災害預警預報以行政村為預警單元，大大縮小了預警面積。這種預警預報可以近似稱為“點狀”地質(zhì)災害預警預報，并以滑坡最具代表性。2005年10月底到2006年3月“點狀預警”思想征求了國內(nèi)有關(guān)專家的意見，從降水預報、動態(tài)雨量數(shù)據(jù)獲取、技術(shù)與管理、信息發(fā)布等綜合能力方面來考慮，溫州市具備建設“點狀”系統(tǒng)建設與運行的能力。因此，溫州市國土資源局在2006年申請立項，市一級地質(zhì)災害預警預報系統(tǒng)建設由此展開。2006年3月，溫州市國土資源局提出建立溫州市地質(zhì)災害預警預報系統(tǒng)的需求;2006年7月至2007年9月，完成了災害數(shù)據(jù)、降雨量數(shù)據(jù)等資料采集和初步分析;2007年10月至2009年6月，因資金沒有落實，預警預報工作暫停開展;2009年1月，招投標成功，同年6月首付資金到位，預警預報系統(tǒng)建設工作開始啟動;2010年6月，初步完成系統(tǒng)建設;2010年7月2011年底，完成系統(tǒng)測試、試運行等工作;2012年至今，正式運行。

圖1 中國地質(zhì)災害預警預報圖Fig.1 The map of warning for geo-hazards in Chinal

圖2 浙江省地質(zhì)災害預警預報圖Fig.2 The map of warning for geo-hazards in Zhejiang province

2 預警方法

2.1 基本思想

基本理論思想是:以近似于“點狀”的行政村為預警單元，實地調(diào)查典型的行政村，并全面收集有關(guān)地形、地貌、基礎地質(zhì)、歷年地質(zhì)災害區(qū)劃調(diào)查等資料，對每個預警單元(行政村)分析其發(fā)生滑坡地質(zhì)災害潛在能力的大小(潛勢度)與災害發(fā)生的降雨閾值，從而建立預警系統(tǒng)，并應用于政府在地質(zhì)災害防災工作中。這種預警以降雨大小為變量，從而有針對性的開展預警工作。

2.2 預警方法的選擇

滑坡地質(zhì)災害預警預報通常是以建立一個系統(tǒng)為載體，而系統(tǒng)的模型選擇就是預警方法的關(guān)鍵步驟。建立預警模型的基本方法通常有三種:統(tǒng)計分析法、機理分析法和其它監(jiān)測方法［1］。統(tǒng)計分析主要是根據(jù)歷史降雨數(shù)據(jù)和地質(zhì)災害發(fā)生情況，進行統(tǒng)計與對比分析，得到地質(zhì)災害和降雨之間的定性、半定量或定量的關(guān)系。機理分析法是采用水文學方法，從降雨-滲流-災害發(fā)生的過程出發(fā)，分析滑坡災害的形成機理。其它監(jiān)測方法指專業(yè)儀器監(jiān)測等，本文不作介紹。

2.2.1 統(tǒng)計分析法

目前，國內(nèi)研究前期有效雨量的精確統(tǒng)計工作正在逐步開展［2］，也有很多學者采用統(tǒng)計學模型分析滑坡災害的時空分布與降雨過程的關(guān)系［3-5］，例如宋光齊等采用類比法預報地質(zhì)災害，該方法將區(qū)域進行單位網(wǎng)格劃分，然后計算每個單元格致災因素的概率值［6］。

2.2.2 機理分析法

機理分析法通常分為抽象模型法、數(shù)值模擬法、室內(nèi)模擬實驗法和現(xiàn)場模擬實驗法［1］。典型的抽象模型法是通過計算安全系數(shù)Fs來確定斜坡失穩(wěn)與否，從而建立降雨量和斜坡失穩(wěn)的關(guān)系。數(shù)值模擬法是運用地下水位、降雨強度等數(shù)值模擬分析斜坡的穩(wěn)定性。當前地質(zhì)災害預警均要分析非降雨因素，因此對地質(zhì)災害潛勢度等研究也日益深入。室內(nèi)模擬實驗法和現(xiàn)場模擬實驗法在日本研究較多，我國現(xiàn)階段一般不采用。

總之，不管是面狀預警還是點狀預警，目前國內(nèi)普遍采用的是綜合統(tǒng)計分析法和機理分析法。劉傳正根據(jù)每個預警區(qū)(中國大陸共74個)的潛勢度、前期雨量和預報雨量三者相結(jié)合來建立第二代預警系統(tǒng)［7］，本論文所論證的滑坡地質(zhì)災害預警方法也采用這種理論和思想。

2.3 建立預警系統(tǒng)的研究步驟

2.3.1 實地調(diào)查及資料分析

根據(jù)溫州近十年滑坡災害發(fā)育的特征，實地調(diào)查典型滑坡地質(zhì)災害的現(xiàn)場情況，從災害發(fā)生的規(guī)模、時間、原因以及誘發(fā)災害的各類環(huán)境因子出發(fā)，調(diào)查該災害點周邊地質(zhì)環(huán)境條件，核對所收集的各類資料，分析各種因素對滑坡災害發(fā)生的影響程度，為劃分預警單元奠定基礎(表1)。

表1 溫州市滑坡預警數(shù)據(jù)采集表Table 1 Data collection table of landslide early warning in Wenzhou city

2.3.2 建立預警單元地質(zhì)模型

誘發(fā)滑坡災害的因素較多，形成機理較為復雜。通常情況下，這些因素可以分為影響內(nèi)因和誘發(fā)外因。影響內(nèi)因就是滑坡本身存在的地質(zhì)環(huán)境因素，包括地形、地貌、巖性、構(gòu)造等。誘發(fā)外因通常是指降雨條件和人為工程建設，其中降雨條件不列入地質(zhì)環(huán)境因子中，作為單獨的變量來進行預警工作。本文所謂的滑坡環(huán)境因子體系除了影響內(nèi)因和誘發(fā)外因之外，還包括特定因子，這些特定因子是指和居民生命財產(chǎn)相關(guān)的各類因素。由于滑坡是一種不良地質(zhì)作用，預警所關(guān)注的滑坡的前提條件是滑坡必須對居民生命和財產(chǎn)產(chǎn)生威脅，如果發(fā)生或可能發(fā)生的滑坡對居民無影響的，則無需預警。因此，本文引入特定因子。

圖3 滑坡預警單元包含的地質(zhì)環(huán)境因子Fig.3 Geological environmental factors about landslide early warning unit

滑坡地質(zhì)環(huán)境因子原始數(shù)據(jù)體系共有16個因子(表1和圖3)，這16個地質(zhì)環(huán)境因子并不是相互獨立的，有些因子之間存在必然的聯(lián)系。我們通過均一性分析、相關(guān)性分析與地質(zhì)分析等方法，最終保留8個地質(zhì)環(huán)境因子(表2)。

表2 滑坡地質(zhì)環(huán)境因子體系Table 2 Geological environmental factors about landslide

2.3.3 建立預警統(tǒng)計模型

8個環(huán)境因子通過敏感性系數(shù)標準化計算后，根據(jù)數(shù)據(jù)類型和特征，選擇主成分分析(PCA)統(tǒng)計模型。通過主成分分析，共提取了7個主成分，前6個主成分的方差總體貢獻率已達到91.663%，它們的方差貢獻率分別為:29.798%，16.343%，12.753%，12.316%，11.473%及8.98%。第一主成分主要反映滑坡體所處建筑區(qū)環(huán)境，第二主成分主要反映滑坡體環(huán)境，第三、四主成分主要反映地質(zhì)環(huán)境，因此，各環(huán)境類別對滑坡的相對貢獻程度排序為:建筑區(qū)環(huán)境＞滑坡體環(huán)境＞地質(zhì)環(huán)境(狹義)。

為了定量化地表示每個滑坡點的潛勢度大小，需計算每個滑坡點綜合得分值。根據(jù)主成分分析的原理，前6個主成分可以用下式來表示:

其中，F(xiàn)是表示每個滑坡點的潛勢度分值，SC1—SC8指的是地貌位置、建筑區(qū)高差、建筑區(qū)最大坡度、滑坡體所在坡面的高差、滑坡體所在坡面的坡度、滑坡體所在坡面的形態(tài)、基巖巖性分類及是否斷層影響區(qū)內(nèi)等八個因子相對敏感系數(shù)標準化后的值。

將樣本數(shù)據(jù)代入后，可以得出每個滑坡點的潛勢度分值，分值的相對大小反映滑坡災害影響因子對滑坡發(fā)生的總的貢獻程度。分值越小，說明發(fā)生滑坡的潛在能力越小，反之，分值越大，說明發(fā)生滑坡點的潛在能力越大，越有利于滑坡災害的發(fā)生。將各滑坡樣本潛勢度分值從小到大排序后，得到一條連續(xù)的曲線(圖4)。

圖4 滑坡潛勢度(PP)分布曲線圖Fig.4 Distribution curve of landslide potentiality parameter

根據(jù)圖4，按等樣本劃分為5個等級，每個等級里包含已發(fā)生(極有可能發(fā)生)滑坡的百分含量，級別越高，百分含量越大。由此，建立5個等級地質(zhì)災害潛勢度分級(PP)，PP1表示發(fā)生地質(zhì)災害的可能性最小，PP5表示可能性最大。

2.3.4 確定降雨閾值(RT)

溫州市降雨類型可分為三類:鋒面雨、臺風雨和東風波暴雨。

(1)臺風雨和東風波暴雨降雨閾值

從降雨類型分析，97%滑坡災害的降雨誘發(fā)因素都來源于臺風雨和東風波雨。根據(jù)誘發(fā)滑坡歷史降雨特征，以及通過2009年“莫拉克”臺風雨誘發(fā)滑坡降雨量分析、模擬和測試，建立了6h累計降雨閾值(圖5、表3)。圖5中，淺色區(qū)域是降雨環(huán)境1區(qū)，表示引發(fā)滑坡的降雨閾值較低;灰色區(qū)域是降雨環(huán)境2區(qū)，表示引發(fā)滑坡的降雨閾值較高，深色是降雨環(huán)境3區(qū)，表示引發(fā)滑坡的降雨閾值最高。

表3 滑坡6h模式的降雨閾值Table 3 Rainfall threshold of landslide 6 hour type

圖5 滑坡6小時模式的降雨閾值分區(qū)圖Fig.5 Partition of rainfall threshold about landslide 6 hour type

(2)鋒面雨降雨閾值

鋒面雨降雨閾值不同于臺風暴雨或者東風波暴雨閾值，在鋒面雨的作用下所發(fā)生的滑坡通常是要考慮前期降雨的影響。從降雨預報的角度，是以24h預報作為1日，其閾值是以日為單位。通過2009年8月9日“莫拉克”臺風事件，同時考慮歷史滑坡實時有效雨量的關(guān)系，其閾值下限為215mm，并建立了5級滑坡有效降雨閾值(表4)。

表4 滑坡24h模式的有效降雨閾值Table 4 Valid rainfall threshold of landslide 24 hour type

2.3.5 建立預警判據(jù)矩陣

據(jù)預警原理，將潛勢度等級(PP)與誘發(fā)滑坡降雨閾值(RT)相組合，從而得到預警判據(jù)矩陣。判據(jù)矩陣交點的預警等級，表示實質(zhì)發(fā)生滑坡概率的大小。根據(jù)《溫州市突發(fā)性地質(zhì)災害應急預案及操作手冊》，突發(fā)性地質(zhì)災害應急響應等級分為5級，1級最輕，5級最嚴重，1～5級分別表示發(fā)生的可能性小、較小、較大、大、很大，分別用亮度表示。因此，所建立的預警判據(jù)也分為5級，和《預案》一致，以便更好的為政府提供技術(shù)支撐。24h滑坡降雨閾值(RT)共有5級，6h滑坡降雨閾值(RT)共有4級，而潛勢度等級(PP)均是5級，因此共組合形成24h預警判據(jù)等級圖和6h預警判據(jù)等級圖(圖6和圖7)。系統(tǒng)建立的5個等級中，1、2級不預警，3～5級預警。

圖6 24h模式的預警判據(jù)圖Fig.6 Warning criterion of 24 hour type

圖7 滑坡6h模式的降雨閾值分區(qū)圖Fig.7 Warning criterion of 6 hour type

3 滑坡災害預警方法的應用

3.1 建立滑坡預警預報系統(tǒng)

通過預警方法理論研究后，結(jié)合MapGIS軟件，建立溫州市滑坡地質(zhì)災害預警系統(tǒng)。溫州市滑坡地質(zhì)災害預警系統(tǒng)的菜單包括以下幾個內(nèi)容:實時雨量、預報雨量、滑坡預警預報、短信發(fā)布、預警預報結(jié)果輸出，常用工具、查詢分析、系統(tǒng)設置。

系統(tǒng)總體框架分為三部分:預警數(shù)據(jù)庫、預警分析、預警結(jié)果顯示及發(fā)布(圖8)。

預警數(shù)據(jù)庫中保存的數(shù)據(jù)包括:原始的雨量數(shù)據(jù)(分為預報雨量和實測雨量)、預警點信息(包括預警點的潛勢度等級以及其它信息)、預警分析的中間結(jié)果數(shù)據(jù)(包括有效雨量計算的結(jié)果等)、預警分析結(jié)果以及預警發(fā)布方案。

圖8 預警系統(tǒng)框架圖Fig.8 Frame diagram of warning system

預警分析功能包括:有效雨量計算、雨量分級、預警級別判定。有效降雨量計算的公式、雨量分級方案以及預警級別判定方案都可以在配置文件中由用戶進行修改。

根據(jù)所使用的雨量數(shù)據(jù)類型，分為預報預警和實時預警兩類;預報時段分為24h和6h兩種，因此預警結(jié)果可分為以下四種:

滑坡24h預報預警:基于前4d前期有效雨量和未來24h降水預報;

滑坡24h實時預警:基于前期120h天實時雨量的實時預警;

滑坡6h預報預警:基于未來6h降水預報;

滑坡6h實時預警:基于前期6h實時雨量。

24h模式涉及有效雨量的概念，不涉及降水分區(qū)，6h模式不涉及有效雨量的概念，只是6h累計雨量，但是涉及降水分區(qū)。

實時預警每小時可以制作，而預報預警(預報)只能在4個時間段制作，分別是:凌晨2時、早上8時、下午2時和晚上8時。

在汛期，每天必須開展的工作是預報預警工作，一般都在早上8點發(fā)布未來24h的地質(zhì)災害預報信息。而實時預警則是根據(jù)預報預警的結(jié)果或者是否收到市防汛辦的強降雨短信為依據(jù)，再考慮是否進行實時預警工作。一般來說，在臺風或者強降雨時段，均會開展實時預警工作。

3.2 在政府防災工作中的應用

3.2.1 預警操作流程

(1)遠程讀取實時雨量數(shù)據(jù)(溫州市水文站提供)。

(2)導入預報雨量(溫州市氣象臺提供);此步驟是預報預警必須經(jīng)歷的，如果做實時預警可以省略此步驟。

(3)生成預警結(jié)果。

(4)導出預警結(jié)果。如果出現(xiàn)3級或者3級以上預警單元，可導出相關(guān)人員的手機號碼和具體的預警單元(行政村點)，并將生成的結(jié)果文件通過MAPGIS軟件生成預警預報成果圖(圖9)。點如果沒有出現(xiàn)三級點，此步驟可以省略。

(5)向公眾發(fā)布預警信息。以2012年09月11日00時所做的滑坡6h實時預警為例，經(jīng)過上述步驟后，生成的預警結(jié)果有三個主要文件，這三個文件分別是“2012年09月11日00時滑坡06小時實時預警成果圖”、“2012年09月11日00時滑坡06小時實時預警行政村點一覽表”、“2012年09月11日00時滑坡06小時實時預警短信列表”。需要公開發(fā)布的是“2012年09月11日00時滑坡06小時實時預警成果圖”、“2012年09月11日00時滑坡06小時實時預警行政村點一覽表”這兩個文件，將這兩個文件制作Word樣版，掛在溫州市國土資源局外網(wǎng)，供公眾瀏覽，同時將預警短信發(fā)給相關(guān)人員，具體接收短信的人員是各預警行政村的村長、鄉(xiāng)鎮(zhèn)長、國土所所長、縣局和市局國土資源管理人員等。

圖9 2012年09月11日00時滑坡06小時實時預警成果圖Fig.9 Real-time warning results figure of landslide 6 hour type in 2012-9-1100:00:00

3.2.2 預警程序

溫州滑坡地質(zhì)災害預警方法以預警預報系統(tǒng)為平臺。這個系統(tǒng)平臺將地質(zhì)環(huán)境因子總影響因素(潛勢度)作為常量(常量做到系統(tǒng)中)，實際工作中以降雨量作為變量，以此來預報或者預警。但是，降雨在時空上是一個復雜的動態(tài)過程，預報預警就具有不確定性，因此需要實時預警來糾正預報預警的實際偏差。

滑坡災害預警在政府工作中的應用分為四個流程:啟動預警、發(fā)布預警、實時(監(jiān)測)預警和解除預警。從啟動預警到發(fā)布預警、解除預警，從預報預警到實時預警，從滑坡24h預警模式到6h預警模式，以及降雨類型和降雨階段的不同，均需要設置判斷條件;啟動預警、實時預警與發(fā)布預警在特殊情況下，可能直接進入實時預警并發(fā)布。

3.2.3 預警響應

當預警結(jié)果最高等級為3級時，由溫州市地質(zhì)環(huán)境監(jiān)測站主管負責人審核簽字同意后發(fā)布預警預報信息;當最高等級為4級時，由市地質(zhì)環(huán)境監(jiān)測站主管負責人審核簽字，送溫州市國土資源局地質(zhì)環(huán)境處負責人簽字(特殊情況可電話、口頭匯報)同意后發(fā)布預警預報信息;當最高等級為5級時，由溫州市國土資源局分管地質(zhì)災害的領(lǐng)導簽發(fā)后再發(fā)布。

當預警結(jié)果發(fā)布后，相關(guān)單位和人員要根據(jù)不同的預警等級采取不同的防御措施。各級地質(zhì)災害主管部門及居民點負責人對地質(zhì)災害預防工作應參照當天最新的滑坡災害預報預警或?qū)崟r預警結(jié)果，對于在等級預報為3級以上的居民點，參考預警等級采取應對措施(表5)，按照突發(fā)性地質(zhì)災害應急預案的安排，對居民點附近危險地段進行嚴密的監(jiān)測，做好預防工作。

表5 滑坡災害預報預警等級與響應措施Table 5 Grade of warning for geo-hazard and response measures

4 結(jié)論

(1)溫州市滑坡地質(zhì)災害的主要類型是淺層滑坡，其誘發(fā)因素是降雨，用地質(zhì)災害區(qū)域預警理論適合溫州市滑坡預警方法。

(2)環(huán)境因子分析結(jié)果:“公路”和“梯田”具有顯著的相似性，考慮溫州山區(qū)村村通公路，這兩個因子對于二分類數(shù)據(jù)而言，意義不大;“建筑區(qū)大小”因子的作用已隱含在“地形環(huán)境”之中;“長度”與“相對高差”相關(guān)性良好，去掉2個長度因子;“殘坡積巖性”、“殘坡積厚度”和“基巖巖性”無法給二分類數(shù)據(jù)中的“否”(未發(fā)生滑坡災害的單元)數(shù)值賦值，因此，共去掉8個原始環(huán)境因子，最終保留剩余的8個因子。

(3)利用主成分分析不僅可以起到降維作用，同時在主成分分析中提取的各個主成分都是各個因子的某種組合，可以更客觀地反映所有因子的共同影響。對8個環(huán)境因子進行分析后，各環(huán)境類別對滑坡的相對貢獻程度排序為:建筑區(qū)環(huán)境＞滑坡體環(huán)境＞地質(zhì)環(huán)境(狹義)。

(4)一次降雨過程，首次發(fā)布滑坡預警(實時)的警戒線為50個，但是管理者可以根據(jù)實際的需要來改變警戒線發(fā)布的臨界個數(shù)。

(5)溫州市滑坡地質(zhì)災害預警方法是滑坡預警工作的理論依據(jù)，但實際預警工作的反饋結(jié)果同樣有助于完善此理論，隨著對溫州滑坡災害研究和調(diào)查程度的不斷深入，溫州市滑坡預警系統(tǒng)會逐漸提高其預警的準確度和精確度。

［1］劉傳正，劉艷輝，溫銘生，等.中國地質(zhì)災害區(qū)域預警方法與應用［M］.北京:地質(zhì)出版社，2009:11-23.LIU Chuanzheng，LIU Yanhui，WEN Mingsheng，et al.Method and application of regional warning for geohazards in Chinal［M］.Beijing:Geological Publishing House，2009:11-23.

［2］張桂榮，殷坤龍，劉禮領(lǐng)，等.基于WEB的浙江省降雨型滑坡預警預報系統(tǒng)［J］.地球科學─中國地質(zhì)大學學報，2005，30(2):250-254.ZHANG Guirong，YIN Kunlong，LIU Liling， et al.Rainfall landslide early warning and forecast system based on WEB in Zhejiang province［J］.Earth Science-Journal of China University of Geosciences，2005，30(2):250-254.

［3］鄭明新，王全才.降雨滑坡的時空預測新方法［J］.中國地質(zhì)災害與防治學報，1997，8(3):45-49.ZHENG Mingxin，WANG Quancai，The new method of temporal prediction about rainfall landslide［J］.The Chinese Journal of Geological Hazard and Control，1997，8(3):45-49.

［4］姚學祥，徐晶，薛建軍.基于降水量的全國地質(zhì)災害潛勢預報模式［J］.中國地質(zhì)災害與防治學報，2005，16(4):97-102.YAO Xuexiang，XU Jing，XUE Jianjun.The potentiality prediction model of geological hazard based on precipitation in China［J］.The Chinese Journal of Geological Hazard and Control，2005，16(4):97-102.

［5］謝守益，張學年，許兵.長江三峽庫區(qū)典型滑坡降雨誘發(fā)的概率分析［J］.工程地質(zhì)學報，1995，3(2):60-69.XIE Shouyi，ZHANG Xuenian，XU Bing.Probabilistic analysis of the Three Gorges landslide induced by rainfall［J］.Journal of Engineering Geology，1995，3(2):60-69.

［6］宋光齊，李云貴，鐘沛林.地質(zhì)災害氣象預報預警方法探討——以四川省地質(zhì)災害氣象預報預警為例［J］.水文地質(zhì)工程地質(zhì)，2004:33-36.SONG Guangqi，LI Yungui，ZHONG Peilin.Discussion on geological disaster weather forecasting and early warning method—The weather forecast and early warning of geological disasters in Sichuan Province as an example［J］.Hydrogeology ＆ Engineering Geology， 2004:33-36.

［7］劉傳正，劉艷輝.地質(zhì)災害區(qū)域預警原理與顯式預警系統(tǒng)設計研究［J］.水文地質(zhì)工程地質(zhì)，2007，34(6):109-115.LIU Chuanzheng，LIU Yanhui.Theory ofregional warning for geo-hazards and research of explicit early warning system design［J］.Hydrogeology ＆ Engineering Geology，2007，34(6):109-115.