范遠(yuǎn)芳，黃俊寶，王國民，陳紅旗

(1.福建省測繪院，福建福州 350003;2.福建省地質(zhì)環(huán)境監(jiān)測中心，福建福州 350002;3.中國地質(zhì)環(huán)境監(jiān)測院(國土資源部地質(zhì)災(zāi)害應(yīng)急技術(shù)指導(dǎo)中心)，北京 100081)

據(jù)全國地質(zhì)災(zāi)害通報(bào)，近十年來約53萬人成功避讓突發(fā)地質(zhì)災(zāi)害。這與2003年以來開始的全國地質(zhì)災(zāi)害氣象預(yù)警預(yù)報(bào)工作，是密不可分的。目前該項(xiàng)業(yè)務(wù)已覆蓋全國地質(zhì)災(zāi)害易發(fā)區(qū)，其預(yù)警預(yù)報(bào)的基礎(chǔ)是對降雨引發(fā)地質(zhì)災(zāi)害風(fēng)險(xiǎn)的識別評估?，F(xiàn)階段預(yù)警預(yù)報(bào)多采用基于歷史數(shù)據(jù)統(tǒng)計(jì)規(guī)律的雨量閾值模型，其核心是對目標(biāo)區(qū)域內(nèi)階段性靜態(tài)風(fēng)險(xiǎn)的一般分布情況的識別。由于滑坡、崩塌和泥石流具有隱蔽性和突發(fā)性，且由于降雨落區(qū)與強(qiáng)度的預(yù)測水平，災(zāi)害預(yù)報(bào)預(yù)警精度仍然較低［1］，有時(shí)會出現(xiàn)過度預(yù)警或者盲目響應(yīng)的現(xiàn)象［2］，進(jìn)而導(dǎo)致預(yù)警信號發(fā)出后，臨機(jī)決策面臨著響應(yīng)級別模糊、靶區(qū)不明和時(shí)機(jī)不好把握等難題。解決這一問題，迫切需要進(jìn)行動態(tài)風(fēng)險(xiǎn)識別評估，即在一次降雨過程中，要對災(zāi)害風(fēng)險(xiǎn)動態(tài)分析、實(shí)時(shí)評估，不斷優(yōu)化預(yù)警等級、區(qū)域和時(shí)機(jī)。

依據(jù)監(jiān)測信息如何及時(shí)識別災(zāi)害風(fēng)險(xiǎn)一直是研究熱點(diǎn)，這一點(diǎn)對預(yù)警與響應(yīng)一體化有現(xiàn)實(shí)意義。在2013年4月召開的第六屆全國地質(zhì)災(zāi)害防治學(xué)術(shù)會議上，有關(guān)動態(tài)風(fēng)險(xiǎn)評估等議題備受關(guān)注。在國外，相關(guān)動態(tài)風(fēng)險(xiǎn)評價(jià)已成為風(fēng)險(xiǎn)評價(jià)領(lǐng)域內(nèi)的發(fā)展方向，動態(tài)風(fēng)險(xiǎn)評價(jià)方法的研究已引起工程界的高度重視［3］，主要見于隧道施工、石油化工、核能、火災(zāi)、環(huán)境保護(hù)、航天工程、經(jīng)濟(jì)等領(lǐng)域［4-5］，并取得較好的應(yīng)用效果。已有的研究表明，風(fēng)險(xiǎn)評價(jià)的動態(tài)性主要體現(xiàn)在兩個方面:一是根據(jù)環(huán)境和自身狀態(tài)的改變，對原有風(fēng)險(xiǎn)評價(jià)模型做出調(diào)整;二是實(shí)現(xiàn)各種風(fēng)險(xiǎn)情景下評價(jià)信息的更新與反饋以及評價(jià)過程的自動化。前者是動態(tài)風(fēng)險(xiǎn)評價(jià)的內(nèi)涵，確保了模型真實(shí)地反映動態(tài)與復(fù)雜的系統(tǒng)。而后者確保了數(shù)據(jù)的動態(tài)性。評價(jià)方法主要有動態(tài)過程評價(jià)、動態(tài)概率評價(jià)、動態(tài)綜合評價(jià)和多智能體評價(jià)等方法［6］。在地質(zhì)災(zāi)害方面，靜態(tài)風(fēng)險(xiǎn)評價(jià)研究已有大量深入研究，成果應(yīng)用于風(fēng)險(xiǎn)治理規(guī)劃方面［7］，而動態(tài)風(fēng)險(xiǎn)評價(jià)方面尚缺乏研究。僅見于文獻(xiàn)［8］，通過分析庫水位升降及降雨作用引發(fā)地下水位的波動從而導(dǎo)致岸坡穩(wěn)定性的變化，按照單體滑坡動態(tài)穩(wěn)定性的理念，進(jìn)行了區(qū)域動態(tài)風(fēng)險(xiǎn)評價(jià)。由于其實(shí)質(zhì)為單體動力預(yù)警預(yù)報(bào)，在區(qū)域上推廣受到限制。本文擬以2012年臺風(fēng)“蘇拉”期間福建省德化縣潯中鎮(zhèn)地質(zhì)災(zāi)害預(yù)警預(yù)報(bào)為例，探索區(qū)域地質(zhì)災(zāi)害動態(tài)風(fēng)險(xiǎn)評價(jià)模型、方法與流程，以期拋磚引玉。

1 研究區(qū)地質(zhì)環(huán)境條件

德化縣地處福建省東南部，地處“閩中屋脊”戴云山脈中段，總體地勢由中部微向四周傾斜，呈層狀梯級下降，戴云山脈西坡陡峻，東坡較緩。地貌以低中山、中低山為主，山前丘陵交錯，其間偶有山間河谷盆地，村莊大都分布于河谷盆地。潯中鎮(zhèn)為德化縣人民政府所在地，面積57.20 km2，境內(nèi)巖漿活動頻繁，侵入巖較發(fā)育，以燕山期侵入巖為主，印支期、華力西期、加里東期零星分布。殘坡積土雙層土體廣泛分布各丘陵臺地、中低山山前地帶坡腳處，多為殘坡積各類粘性土，大多為雙層土體，一般厚2～4 m，局部地段厚度大于10 m，結(jié)構(gòu)松散，為本區(qū)地質(zhì)災(zāi)害的主要致災(zāi)體。中低山之間的盆地邊緣地帶則是地質(zhì)災(zāi)害分布密集區(qū)，以淺層小型滑坡崩塌為主(表1)。

表1 德化縣歷史崩滑統(tǒng)計(jì)表Table 1 The statistics of historical landslides of Dehua

區(qū)內(nèi)氣候?qū)賮啛釒降貧夂?，雨量充沛，潮濕多霧。多年平均總降雨量1867.2 mm。降雨量在時(shí)間上可劃分為三個降雨季節(jié)，其中7～9月為臺風(fēng)雷陣雨季，多年平均降雨量782.2 mm，約占全年總降雨量的39.9%，表現(xiàn)為臺風(fēng)強(qiáng)降雨，雨量大而集中，強(qiáng)度強(qiáng)，降雨的年際變化大，是地質(zhì)災(zāi)害易發(fā)期，局部區(qū)域具有多點(diǎn)同發(fā)、群發(fā)的特點(diǎn)。

2012年第9號強(qiáng)臺風(fēng)“蘇拉”于7月28日08時(shí)在菲律賓馬尼拉以東約630 km的西北太平洋洋面上生成，8月3日6:50在寧德福鼎秦嶼鎮(zhèn)登陸。臺風(fēng)強(qiáng)降雨影響德化縣潯中鎮(zhèn)，并加劇了地質(zhì)災(zāi)害風(fēng)險(xiǎn)水平。

2 地質(zhì)災(zāi)害風(fēng)險(xiǎn)動態(tài)過程評價(jià)

地質(zhì)災(zāi)害動態(tài)風(fēng)險(xiǎn)評價(jià)是依據(jù)災(zāi)害地質(zhì)條件、誘發(fā)因素和承災(zāi)體的實(shí)時(shí)狀態(tài)，對災(zāi)害風(fēng)險(xiǎn)的連續(xù)性識別，服務(wù)于應(yīng)急管理。

2.1 評價(jià)模型

在現(xiàn)有的動態(tài)風(fēng)險(xiǎn)評價(jià)方法中，動態(tài)概率風(fēng)險(xiǎn)評價(jià)方法使用的難度較大;動態(tài)綜合評價(jià)方法重視后果評價(jià)，對于強(qiáng)調(diào)預(yù)警預(yù)報(bào)的地質(zhì)災(zāi)害風(fēng)險(xiǎn)評價(jià)適用性不強(qiáng);基于多智能體的評價(jià)方法，需要構(gòu)建災(zāi)害各微觀要素之間相互作用，尚不具備理論基礎(chǔ)。相比而言，動態(tài)過程評價(jià)能實(shí)現(xiàn)各種風(fēng)險(xiǎn)情景下評價(jià)信息的更新與反饋以及評價(jià)過程的自動化，可操作性強(qiáng)，隨著評價(jià)次數(shù)的增加，數(shù)據(jù)的更新與修正也能使結(jié)果趨于完善和真實(shí)。為此，這里采用動態(tài)過程評價(jià)方法研究。

公式(1)所示的地質(zhì)災(zāi)害風(fēng)險(xiǎn)評估框架模型［9］。在一次降雨過程中，動態(tài)因子有實(shí)時(shí)降雨、災(zāi)情、險(xiǎn)情和承災(zāi)體狀態(tài)，依托群測群防監(jiān)測體系，均可以實(shí)現(xiàn)實(shí)時(shí)反饋更新。在靜態(tài)模型基礎(chǔ)上考慮上述動態(tài)因子，依據(jù)動態(tài)過程評價(jià)方法，建立動態(tài)風(fēng)險(xiǎn)評價(jià)模型，見公式(2)。在危險(xiǎn)性動態(tài)方面，誘發(fā)因素降雨分解為有效雨量［10］和實(shí)時(shí)雨強(qiáng);災(zāi)情險(xiǎn)情是指過程中已經(jīng)出現(xiàn)的災(zāi)害后果和災(zāi)害跡象，直接顯示危險(xiǎn)性狀態(tài)，根據(jù)相似原理由點(diǎn)狀監(jiān)測信息輻射到評價(jià)區(qū)域，對區(qū)域易發(fā)性進(jìn)行動態(tài)更新。在易損性方面，以承災(zāi)體狀態(tài)修正靜態(tài)易損值。通過因子動態(tài)分析，提出修正模型如下:

式中，Rt為總風(fēng)險(xiǎn)，H為危險(xiǎn)值，S為易發(fā)性，R為降雨，V為易損值，D為險(xiǎn)情，B為承災(zāi)體狀態(tài)，下標(biāo)d表示動態(tài)因子，沒有下標(biāo)表示靜態(tài)因子。

2.2 評價(jià)方法

按照上述更新的實(shí)時(shí)動態(tài)因子，采用地理信息系統(tǒng)進(jìn)行動態(tài)風(fēng)險(xiǎn)評價(jià)區(qū)劃。評價(jià)流程見圖1，包括5個主要評價(jià)步驟:

圖1 評價(jià)流程Fig．1 The flow of assessment

第一步是原型分析，進(jìn)行評價(jià)區(qū)描述。

第二步是動態(tài)因子監(jiān)測與信息采集更新;本文受數(shù)據(jù)與監(jiān)測時(shí)間間隔所限，采用1 h間隔進(jìn)行評價(jià)。若實(shí)現(xiàn)數(shù)據(jù)監(jiān)測、采集、傳輸、運(yùn)算的自動化，則可以實(shí)現(xiàn)風(fēng)險(xiǎn)評價(jià)過程的自動化。

第三步是依據(jù)災(zāi)情險(xiǎn)情監(jiān)測信息，依據(jù)地質(zhì)環(huán)境相似度計(jì)算，更新易發(fā)性分區(qū);一旦出現(xiàn)險(xiǎn)情或?yàn)?zāi)情，易發(fā)性區(qū)劃將在原來基礎(chǔ)上進(jìn)行更新。在后續(xù)過程中如果再次出現(xiàn)險(xiǎn)情或?yàn)?zāi)情，將以最新的易發(fā)性區(qū)劃為基礎(chǔ)進(jìn)行進(jìn)一步更新，如此反復(fù)，實(shí)現(xiàn)易發(fā)性分區(qū)的動態(tài)。隨著更新次數(shù)的增加，易發(fā)性分區(qū)將愈趨于準(zhǔn)確。

第四步是依據(jù)災(zāi)情險(xiǎn)情及避險(xiǎn)防范狀態(tài)，采用事件樹方法更新人員和道路兩項(xiàng)主要易損性指標(biāo)。事件樹可以從初始事件開始，對事件發(fā)展可能導(dǎo)致的各種后果進(jìn)行分析。風(fēng)險(xiǎn)演變過程中，降雨、險(xiǎn)情與災(zāi)情與人員狀態(tài)和道路狀態(tài)兩項(xiàng)因子為因果關(guān)系。降雨、險(xiǎn)情與災(zāi)情的變化將自動引發(fā)人員狀態(tài)與道路狀態(tài)的變化，將導(dǎo)致易損性區(qū)劃的動態(tài)變化。

第五步是基于模糊綜合評判的動態(tài)風(fēng)險(xiǎn)區(qū)劃。

依據(jù)研究區(qū)域地質(zhì)環(huán)境與地質(zhì)災(zāi)害特征，德化縣98.6%的滑坡為土質(zhì)小型滑坡。采用德化縣地質(zhì)災(zāi)害調(diào)查數(shù)據(jù)與各因子進(jìn)行統(tǒng)計(jì)分析，結(jié)果表明研究區(qū)滑坡主要受坡度、坡型、表土層厚度以及坡向控制。其中土體厚度3～4 m最易發(fā)，地質(zhì)災(zāi)害多發(fā)生于25°～35°凹凸坡坡度區(qū)，主要集中在坡向北向與西向(圖2)。故在危險(xiǎn)性因子方面，易發(fā)性通過坡度、坡向、表土層厚度、坡形加以反映，通過實(shí)施災(zāi)情險(xiǎn)情監(jiān)測信息進(jìn)行相似度修正;誘發(fā)因素降雨［11］，按照實(shí)際群測群防監(jiān)測周期及速報(bào)頻次，采用德化、三班、雷峰、蓋德和鳳洋五個雨量站臺風(fēng)“蘇拉”期間1 h時(shí)間間隔的監(jiān)測數(shù)據(jù)，更新有效降雨量與雨強(qiáng)。在易損性因子方面，從應(yīng)急管理的可操作性出發(fā)，重點(diǎn)考慮人與道路的承載狀態(tài)。評價(jià)指標(biāo)權(quán)重通過層次分析法確定［12］。確定評價(jià)指標(biāo)如表2所示。

(1)易發(fā)性動態(tài)相似度計(jì)算

當(dāng)某地出現(xiàn)變形跡象或發(fā)生崩滑災(zāi)害等險(xiǎn)情時(shí)，獲取險(xiǎn)情所在地區(qū)的地質(zhì)環(huán)境因素，包括表土層厚度、坡度、坡向、坡型四個作為參照因子。應(yīng)用ArcGIS對單個圖層各個單元與相應(yīng)的參照因子比對，獲取與參照因子相同的區(qū)域，同理，其它三個圖層采用相同的操作。將獲得的四個圖層進(jìn)行疊加分析，具有四個因子相同的單元賦值為1，其它單元賦值為0。將新圖層與易發(fā)區(qū)圖進(jìn)行疊加，將單元屬性值為1的單元作為高易發(fā)單元，進(jìn)行易發(fā)性區(qū)劃動態(tài)更新。

(2)承災(zāi)體狀態(tài)事件樹分析

事件樹分析(Event Tree Analysis，簡稱ETA)是一種邏輯演繹方法，由決策樹分析(簡稱DTA)發(fā)展而來。它由初始事件根據(jù)事件發(fā)展的時(shí)間順序進(jìn)行初始事件可能導(dǎo)致的各種后果進(jìn)行分析。具有可以定性與定量地進(jìn)行系統(tǒng)故障原因以及后果識別作用［13］。在降雨過程中，承災(zāi)體所處的狀態(tài)是隨著降雨、險(xiǎn)情、災(zāi)情狀況動態(tài)變化的。下面對人員狀態(tài)與道路狀態(tài)進(jìn)行事件樹分析，并進(jìn)行指標(biāo)更新。

圖2 滑坡與各因子統(tǒng)計(jì)分析Fig．2 The statistical analysis of landslides and each factors

表2 崩滑地質(zhì)災(zāi)害動態(tài)風(fēng)險(xiǎn)評價(jià)指標(biāo)與權(quán)重Table 2 The index and weight of landslide geological hazard dynamic risk evaluation

人員狀態(tài)變化事件樹。降雨為激發(fā)事件，在降雨激發(fā)下可能發(fā)生一系列子事件:降雨達(dá)到相應(yīng)等級臨界值，降雨臨界值通過歷史滑坡與降雨的統(tǒng)計(jì)數(shù)據(jù)得到;地表出現(xiàn)宏觀變形破壞、地下水異常等險(xiǎn)情;發(fā)生崩滑災(zāi)害，出現(xiàn)災(zāi)情等。選分支1說明，降雨達(dá)到臨界值2，根據(jù)應(yīng)急預(yù)案，人員應(yīng)當(dāng)轉(zhuǎn)移，而非隱患點(diǎn)人員將處于自然狀態(tài)。當(dāng)出現(xiàn)險(xiǎn)情或?yàn)?zāi)情時(shí)，人員轉(zhuǎn)移。否則繼續(xù)為自然狀態(tài)。下一時(shí)刻基于前一時(shí)刻的結(jié)果，仍以降雨為初始事件，已經(jīng)為轉(zhuǎn)移狀態(tài)的人員將不再進(jìn)行判斷，對自然狀態(tài)的人員繼續(xù)分析子事件，如此循環(huán)，進(jìn)行動態(tài)分析(圖3)。道路狀態(tài)變化事件樹分析。采用相同的方法進(jìn)行道路狀態(tài)分析(圖4)。

(3)模糊綜合評判

在處理難以用精確數(shù)學(xué)方法進(jìn)行描述的復(fù)雜系統(tǒng)問題時(shí)，表現(xiàn)出了獨(dú)特的優(yōu)越性，所以被廣泛應(yīng)用于各個學(xué)科領(lǐng)域。根據(jù)評價(jià)精度以及建立模糊集合所要滿足的對稱性原則，將評判等級標(biāo)準(zhǔn)分為“較低風(fēng)險(xiǎn)”、“中等風(fēng)險(xiǎn)”，“較高風(fēng)險(xiǎn)”和“高風(fēng)險(xiǎn)”四個等級，即評語集為:V={Ⅰ Ⅱ Ⅲ Ⅳ}={較低風(fēng)險(xiǎn)中等風(fēng)險(xiǎn)較高風(fēng)險(xiǎn) 高風(fēng)險(xiǎn)}，各風(fēng)險(xiǎn)等級定性描述如表3所示。從評價(jià)指標(biāo)體系與模糊評語集出發(fā)，將定性與定量化的指標(biāo)按照統(tǒng)一量綱的原則，確定簡易的評判標(biāo)準(zhǔn)。對于定性指標(biāo)，首先劃分出因子的風(fēng)險(xiǎn)級別，采用經(jīng)驗(yàn)賦值的方法確定隸屬度。對于降雨量，人口密度等定量指標(biāo)的隸屬函數(shù)，采用梯形分布函數(shù)來刻畫，即第i個指標(biāo)的風(fēng)險(xiǎn)等級I～I(xiàn)V隸屬度計(jì)算公式如下:

圖3 人員狀態(tài)變化事件樹Fig．3 The event tree of personnel status changed

圖4 道路狀態(tài)變化事件樹Fig．4 The event tree of road state changed

表3 動態(tài)風(fēng)險(xiǎn)定性描述Table 3 The qualitative description of dynamic Risk

其中，ui(xi)為各評價(jià)因子的隸屬函數(shù);xi為評價(jià)因子的實(shí)際取值;ai，bi，ci為第i個評價(jià)因子各等級標(biāo)準(zhǔn)值域的界限值。對于中間越優(yōu)的指標(biāo)，先采用x'=|xi-xmid|運(yùn)算求得x'后代入上面隸屬度函數(shù)進(jìn)行計(jì)算，值得注意的是，ai，bi，ci的值應(yīng)為x'的各等級中間值。評判各因子由于各有側(cè)重，需對各因子賦予權(quán)重，然后采用模糊綜合評判模型進(jìn)行三級模糊綜合評判。最后對評判結(jié)果作歸一化處理，根據(jù)最大隸屬原則確定樣本歸屬。為了將各因子對地質(zhì)災(zāi)害動態(tài)風(fēng)險(xiǎn)的影響程度都反映出來，需采用M(·，+)模糊算子進(jìn)行模糊變換。

根據(jù)以上方法，采用GIS進(jìn)行動態(tài)風(fēng)險(xiǎn)評價(jià)建模實(shí)現(xiàn)，包括網(wǎng)格劃分、柵格模型量化、有權(quán)疊加運(yùn)算。

2.3 評價(jià)結(jié)果分析

以圖5所示三個時(shí)刻的風(fēng)險(xiǎn)評價(jià)結(jié)果為例(表4)，8月2日15時(shí)降雨等級較低，風(fēng)險(xiǎn)普遍偏低，石鼓與后所兩處隱患點(diǎn)所處區(qū)域?yàn)橹械蕊L(fēng)險(xiǎn);3日10時(shí)降雨轉(zhuǎn)為暴雨，風(fēng)險(xiǎn)普遍增高，高風(fēng)險(xiǎn)區(qū)域廣泛分布在石山村、仙境村、石鼓村、龍翰村、鳳池、富東大部分地區(qū)、以及鳳洋南部、蒲坂村周邊地區(qū)。根據(jù)應(yīng)急預(yù)案，對高風(fēng)險(xiǎn)區(qū)域人員進(jìn)行了避險(xiǎn)防范部署;3日23時(shí)，高風(fēng)險(xiǎn)區(qū)域面積擴(kuò)大。雖然本時(shí)刻降雨已由暴雨轉(zhuǎn)為中到大雨，但地質(zhì)災(zāi)害的發(fā)生具有滯后效應(yīng)，龍翰村9組等發(fā)生多處滑坡。根據(jù)風(fēng)險(xiǎn)動態(tài)評價(jià)結(jié)果，優(yōu)化預(yù)案，并提高了應(yīng)急準(zhǔn)備的針對性。

圖5 三個時(shí)刻風(fēng)險(xiǎn)評價(jià)結(jié)果圖Fig．5 The results of the risk assessment map of the three moments

3 結(jié)論與認(rèn)識

表4 三個時(shí)刻風(fēng)險(xiǎn)評價(jià)結(jié)果對比表Table 4 The comparison table of three-time risk assessment results

在一次降雨過程中，對地質(zhì)災(zāi)害風(fēng)險(xiǎn)動態(tài)的識別是提高預(yù)警精度的關(guān)鍵因素之一。本文嘗試采用過程評價(jià)方法，進(jìn)行降雨型地質(zhì)災(zāi)害風(fēng)險(xiǎn)動態(tài)評價(jià)。通過易發(fā)性動態(tài)相似度更新有效雨量實(shí)時(shí)監(jiān)測、承災(zāi)狀態(tài)等過程分析，可建立動態(tài)風(fēng)險(xiǎn)評價(jià)修正模型。評價(jià)流程為:區(qū)域原型分析→監(jiān)測信息采集→易發(fā)性更新→易損性更新→動態(tài)風(fēng)險(xiǎn)區(qū)劃。需要指出的是，由于參評要素實(shí)時(shí)監(jiān)測采集條件尚有限，使該方法目前僅適用于防災(zāi)減災(zāi)基礎(chǔ)較好的重點(diǎn)區(qū)域。

［1］劉振國，謝必如.徐紹史在重慶三峽庫區(qū)專題調(diào)研時(shí)強(qiáng)調(diào)充分認(rèn)識地質(zhì)災(zāi)害三大特點(diǎn)依靠群測群防做好防治工作.［EB/OL］.http://www.mlr.gov.cn/xwdt/jrxw/201008/t20100802_729661.htm.2010-08-02/2013-12-2.

［2］劉維.以人為本防治地災(zāi)——《全國地質(zhì)災(zāi)害防治“十二五”規(guī)劃》解讀.［EB/OL］.http://www.mlr.gov. cn/xwdt/jrxw/201204/t20120412 _1083503.htm.2012-04-12/2013-12-2.

［3］Aldemir T，Siu N.Special issue on reliability and safety of dynamic process systems［J］.Reliability Engineering and Systems Safety，1996，52(3):181-347.

［4］繆春生.動態(tài)風(fēng)險(xiǎn)評價(jià)方法的研究進(jìn)展［A］.第一次全國工程風(fēng)險(xiǎn)分析技術(shù)學(xué)術(shù)會議論文集［C］.南京:2006:418-429.MIAO Chunsheng.TheOverview:DynamicRisk Assessment［A］.The first National Engineering Risk Analysis Technology Symposium［C］.Nanjing:2006:418-429.

［5］易立新，惠中玉.火災(zāi)風(fēng)險(xiǎn)評價(jià)的隨機(jī)過程模型［J］.中國安全科學(xué)學(xué)報(bào)，1997，7(5):6-9.YI Lixin，HUI Zhongyu.Stochastic model for fir risk assessment［J］.China Safety Science Journal，1997，7(5):6-9.

［6］范遠(yuǎn)芳，張小趁，褚宏亮，等.區(qū)域地質(zhì)災(zāi)害風(fēng)險(xiǎn)靜動態(tài)評估對比研究［A］.第九屆全國工程地質(zhì)大會論文集［C］.青島:科學(xué)出版社，2012:510-515.FAN Yuanfang，ZHANG Xiaochen，CHU Hongliang，et al. Comparative study of static and dynamic assessment of regional geo-hazard risk［A］.Proceedings ofthe 9th nationalconference on engineering geology［C］.Qingdao:Science Press，2012:510-515.

［7］董穎.地質(zhì)災(zāi)害風(fēng)險(xiǎn)評估理論與實(shí)踐［M］.北京:地質(zhì)出版社，2009.DONG Ying.The Chinese journal of geological hazard and control［M］. Beijing:GeologicalPublishing House，2009.

［8］馮文凱，石豫川，柴賀軍，等.三峽庫區(qū)渝巴公路馬道子滑坡岸段動態(tài)風(fēng)險(xiǎn)評價(jià)［J］.2007，22(1):54-59.FENG Wenkai，SHI Yuchuan，CHAI Hejun，et al.Dynamic risk evaluation about Madaozi landslide’s bank slope of Yu-ba highway in Three Gorges Area［J］.Journal of Catastrophology，2007，22(1):54-59.

［9］朱良峰，殷坤龍，張梁，等.地質(zhì)災(zāi)害風(fēng)險(xiǎn)分析與GIS技術(shù)應(yīng)用研究［J］.地理學(xué)與國土研究，2002，18(4):10-13.ZHU Liangfeng，YIN Kunlong，ZHANG Liang，et al.Risk analysis of geo-hazard and application of the GIS technique［J］.Geography and Territorial Research，2002，18(4):10-13.

［10］唐紅梅，唐云輝，高陽華.重慶地區(qū)降雨型滑坡相關(guān)性分析及預(yù)報(bào)模型［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2013，24(4):16-22.TANG Hongmei，TANG Yunhui，GAO Yanghua.Correlation analysis and prediction model for rainfallinduced landslide in Chongqing area［J］.The Chinese Journal of Geological Hazard and Control，2013，24(4):16-22.

［11］杜惠良，鈕學(xué)新，殷坤龍，等.浙江省滑波，泥石流多發(fā)區(qū)氣象預(yù)警研究［J］.高原氣象，2006，25(1):151-158.DU Huiliang，NIU Xuexin，YIN Kunlong，et al.Study of meteorological forecasting and warning for landslide and debris flow in Zhejiang province［J］.Plateau Meteorology，2006，25(1):151-158.

［12］唐亞明，程秀娟，薛強(qiáng)，等.基于層次分析法的黃土滑塌風(fēng)險(xiǎn)評價(jià)指標(biāo)權(quán)重分析［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2012，23(4):40-46.TANG Yaming，CHENG Xiujuan，XUE Qiang，et al.Weights analysis of Loess collapse risk assessing factors based on analytical hierarchy process［J］.The Chinese Journal of Geological Hazard and Control，2012，23(4):40-46.

［13］解家畢，孫東亞.事件樹法原理及其在堤壩風(fēng)險(xiǎn)分析中的應(yīng)用［J］.中國水利水電科學(xué)研究院學(xué)報(bào)，2006，4(2):133-137.XIE Jiabi，SUN Dongya.The constructing principle of event tree method and its application in risk analysis for dykes and dams［J］.Journal of China Institute of Water Resources and Hydropower Research，2006，4(2):133-137.