南驍聰，李 寧

（1.山東高速工程檢測(cè)有限公司，山東 濟(jì)南 250002；2.山東高速股份有限公司,山東 濟(jì)南 250101）

引言

我國(guó)約70%都為山區(qū)，降雨和地震等頻發(fā)，誘發(fā)了大量的滑坡災(zāi)害，平均年經(jīng)濟(jì)損失超過(guò) 20 億元，造成超過(guò)千人的生命損失。此外，善后處理以及事后整治等費(fèi)用也達(dá)數(shù)十億元，滑坡間接產(chǎn)生的其他經(jīng)濟(jì)損失（如交通中斷）更是無(wú)法估計(jì)。為了減輕滑坡造成的后果，決策者需要對(duì)滑坡災(zāi)害有更好的認(rèn)識(shí)，并對(duì)滑坡風(fēng)險(xiǎn)進(jìn)行管理。要管理風(fēng)險(xiǎn)，首先必須對(duì)風(fēng)險(xiǎn)進(jìn)行量化。滑坡災(zāi)害定量風(fēng)險(xiǎn)評(píng)估作為一種更加合理有效的風(fēng)險(xiǎn)量化手段，一直是國(guó)際上倡導(dǎo)和推廣的有效途徑之一。因此，開(kāi)展滑坡災(zāi)害定量風(fēng)險(xiǎn)評(píng)估對(duì)保護(hù)人民生命和財(cái)產(chǎn)安全以及維護(hù)社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展都有著重要的現(xiàn)實(shí)意義。

1 滑坡定量風(fēng)險(xiǎn)評(píng)估框架

滑坡風(fēng)險(xiǎn)評(píng)價(jià)經(jīng)歷了定性-半定量-定量的過(guò)程。其中，對(duì)于定性及半定量風(fēng)險(xiǎn)評(píng)價(jià)，主要是通過(guò)搜集歷史數(shù)據(jù)采用統(tǒng)計(jì)方法，并未充分考慮災(zāi)害成因機(jī)制及運(yùn)移特征，且易損性因子的賦值基本依賴于專家經(jīng)驗(yàn)法，對(duì)災(zāi)害風(fēng)險(xiǎn)的量化，較少考慮災(zāi)害自身的物理機(jī)制。風(fēng)險(xiǎn)定量評(píng)價(jià)不僅考慮影響災(zāi)害強(qiáng)度的物理力學(xué)參數(shù)的不確定性，同時(shí)也應(yīng)考慮影響受災(zāi)體易損性因素的不確定性，量化災(zāi)害風(fēng)險(xiǎn)分布。隨著數(shù)值模擬手段的廣泛運(yùn)用，并結(jié)合概率理論等數(shù)學(xué)方法，使得對(duì)滑坡災(zāi)害風(fēng)險(xiǎn)量化成為可能。風(fēng)險(xiǎn)評(píng)價(jià)不再僅是定性地劃分等級(jí)，而是直接得到某一災(zāi)害情境下的損失概率和損失值?；嘛L(fēng)險(xiǎn)評(píng)估需要解決的問(wèn)題:（1）滑坡發(fā)生的概率；（2）滑坡的運(yùn)動(dòng)過(guò)程；（3）財(cái)產(chǎn)和人的易損性；（4）財(cái)產(chǎn)和人的風(fēng)險(xiǎn)。

通常用概率方法定量描述滑坡事件的不確定性，單體滑坡定量風(fēng)險(xiǎn)評(píng)估數(shù)學(xué)模型：

式中：R—單體滑坡風(fēng)險(xiǎn)值；P1—滑坡在給定回歸期內(nèi)以某種規(guī)模失穩(wěn)的概率，P2—滑坡物質(zhì)以某種運(yùn)動(dòng)強(qiáng)度到達(dá)某位置的概率；P3—承災(zāi)體在空間上出現(xiàn)在該位置的概率；P4—承災(zāi)體損壞或人員傷亡的概率；C—承災(zāi)體的價(jià)值或人口數(shù)量。

2 單體滑坡定量風(fēng)險(xiǎn)評(píng)價(jià)方法

2.1 發(fā)生概率評(píng)價(jià)方法

對(duì)于單體滑坡發(fā)生概率的基本思路是在傳統(tǒng)的定值分析方法（Slope/W, Slide, Analysis和Flac2D/3D等數(shù)值模擬方法）的基礎(chǔ)上，將輸入?yún)?shù)定義為服從某種概率分布的隨機(jī)變量，得到穩(wěn)定性系數(shù)的概率分布，從而求解出穩(wěn)定系數(shù)小于1的概率，即滑坡發(fā)生的概率。在進(jìn)行滑坡失穩(wěn)概率分析時(shí), 通常需建立極限狀態(tài)函數(shù)（又稱功能函數(shù)）：

式中：X—隨機(jī)變量；FS（X）—利用數(shù)值模擬在X參數(shù)下計(jì)算出的穩(wěn)定性系數(shù)。

g（X）＜0的概率即滑坡發(fā)生概率。

一旦極限狀態(tài)函數(shù)確定，即可以通過(guò)蒙特卡羅模擬方法估計(jì)滑坡發(fā)生的概率，該方法是求解失效概率最簡(jiǎn)單、最有效的方法。（1）對(duì)影響滑坡穩(wěn)定性計(jì)算的隨機(jī)變量（如粘聚力、內(nèi)摩擦角和孔隙水壓力等）進(jìn)行隨機(jī)抽樣；（2）將隨機(jī)抽樣值逐個(gè)代入極限狀態(tài)函數(shù)中計(jì)算出功能函數(shù)值；（3）統(tǒng)計(jì)功能函數(shù)值小于零的個(gè)數(shù)，并根據(jù)總采樣次數(shù)確定滑坡的發(fā)生頻率，當(dāng)采樣次數(shù)足夠大時(shí)，得到的滑坡發(fā)生的頻率將逼近滑坡發(fā)生的概率，從而求解滑坡發(fā)生概率。然而，蒙特卡羅模擬方法的計(jì)算量非常大，當(dāng)遭遇小概率問(wèn)題時(shí)（如10-3），需要計(jì)算次數(shù)將達(dá)105次以上，占用大量計(jì)算時(shí)間。為了減少計(jì)算量，一些高效的可靠度方法已經(jīng)應(yīng)用于滑坡發(fā)生概率的計(jì)算，如一階可靠度方法、二階可靠度方法以及響應(yīng)面方法等。

2.2 滑坡運(yùn)動(dòng)過(guò)程評(píng)價(jià)方法

劃定滑坡危險(xiǎn)區(qū)范圍是滑坡風(fēng)險(xiǎn)評(píng)估的基礎(chǔ)。要求準(zhǔn)確預(yù)測(cè)滑坡的運(yùn)動(dòng)過(guò)程，獲得表征滑坡破壞潛力的參數(shù)，從而用于單體滑坡的定量風(fēng)險(xiǎn)評(píng)估。用于滑坡風(fēng)險(xiǎn)評(píng)估參數(shù)主要包括：運(yùn)動(dòng)距離、運(yùn)動(dòng)速度、堆積厚度以及沖擊力。隨著滑坡運(yùn)動(dòng)過(guò)程數(shù)值模擬的發(fā)展，使得滑坡運(yùn)動(dòng)過(guò)程的精細(xì)化模擬得以實(shí)現(xiàn)，并能直接輸出風(fēng)險(xiǎn)評(píng)估所需參數(shù)。目前較為常用的數(shù)值模擬方法有連續(xù)介質(zhì)力學(xué)方法和離散元方法。由于大多數(shù)滑坡、巖崩、碎屑流、泥石流都具有像流體一樣的性質(zhì)，Savage和Hutter提出了基于等效流體假設(shè)的流體力學(xué)方法，基于傳統(tǒng)的N-S方程，采用深度積分的方式，選擇合適的流體力學(xué)模型，實(shí)現(xiàn)連續(xù)介質(zhì)力學(xué)方法對(duì)滑坡動(dòng)態(tài)運(yùn)動(dòng)過(guò)程的模擬。目前該方法被廣泛運(yùn)用于滑坡的運(yùn)動(dòng)過(guò)程模擬中，并開(kāi)發(fā)各種不同的計(jì)算軟件，如DAN3D，MassMov2D，F(xiàn)LO-2D和RAMMS等。該方法具有計(jì)算參數(shù)少、計(jì)算效率高的特點(diǎn)。與連續(xù)介質(zhì)力學(xué)方法相比，離散元法能夠反映出滑坡運(yùn)動(dòng)過(guò)程中顆粒之間的分離與聚集效果，可以有效模擬顆粒大變形運(yùn)動(dòng)，是分析理解滑坡破壞機(jī)制的寶貴工具。該方法由美國(guó)學(xué)者Cundall和Strack于1971年首次提出，經(jīng)過(guò)近50年的發(fā)展，已開(kāi)發(fā)形成多款大型商業(yè)軟件，如PFC2D與PFC3D等。

2.3 易損性評(píng)價(jià)方法

在定量風(fēng)險(xiǎn)評(píng)估中，易損性通常以在一定災(zāi)害強(qiáng)度作用下影響區(qū)內(nèi)承災(zāi)體的損失程度或損失價(jià)值來(lái)衡量。它是一種條件概率，假定發(fā)生了具有一定規(guī)模的滑坡，承災(zāi)體在滑坡的路徑上。根據(jù)滑坡與工程邊坡聯(lián)合技術(shù)委員會(huì)對(duì)易損性的定義，易損性的值為0（沒(méi)有損失）和1（完全損失）之間的無(wú)量綱系數(shù)。易損性主要包括物質(zhì)易損性、經(jīng)濟(jì)易損性、環(huán)境易損性和社會(huì)易損性等4類。其中，物質(zhì)易損性主要考慮有形資產(chǎn)（如基礎(chǔ)設(shè)施和建筑物等），反映的是承災(zāi)體的自然特性；而后三者則涉及無(wú)形資產(chǎn)（如經(jīng)濟(jì)、人口以及資源環(huán)境），反映了承災(zāi)體的社會(huì)特性。易損性評(píng)估需要了解滑坡事件和承災(zāi)體之間的相互作用，即易損性曲線。目前，易損性研究大多根據(jù)滑坡災(zāi)害歷史資料，運(yùn)用統(tǒng)計(jì)模型構(gòu)建易損性曲線；也有一部分研究通過(guò)對(duì)建筑物進(jìn)行有限元模擬，構(gòu)建沖擊壓力與建筑物形變位移的關(guān)系，從而建立易損性曲線，這種更加基于物理的易損性曲線建立方法是未來(lái)研究的一種趨勢(shì)。

3 滑坡運(yùn)動(dòng)過(guò)程分析

在滑坡風(fēng)險(xiǎn)評(píng)估中，滑坡運(yùn)動(dòng)過(guò)程分析是非常重要的一部分。某黃土滑坡位于甘肅省內(nèi)，處于黃河的第四階地上。由于長(zhǎng)期的灌溉作用，灌溉水下滲，黃土強(qiáng)度軟化，發(fā)生滑坡?；麻L(zhǎng)217 m，寬176 m，滑坡平均厚度約20 m，滑坡體積大約44.4×104 m3，滑距770 m。數(shù)值模擬方法使用基于連續(xù)理論的深度積分方法，并通過(guò)二維剖面建立計(jì)算模型。基底阻力模型使用庫(kù)倫摩擦模型，包括黏聚力和摩擦角兩個(gè)參數(shù)?？紤]到黃土強(qiáng)度發(fā)生軟化，對(duì)于該滑坡，計(jì)算模型的黏聚力設(shè)置為10 kPa，摩擦角設(shè)置為8.5°。計(jì)算網(wǎng)格間距2 m，計(jì)算總時(shí)間40 s，時(shí)間步設(shè)置為0.001 s?；氯^(guò)程運(yùn)動(dòng)過(guò)程模擬結(jié)果見(jiàn)圖1?？梢园l(fā)現(xiàn)，數(shù)值模擬可以細(xì)致地重現(xiàn)滑坡的運(yùn)動(dòng)過(guò)程。根據(jù)最后的堆積結(jié)果（40 s時(shí)）可以發(fā)現(xiàn)，滑坡上部仍然有部分殘留在滑面上，與實(shí)際結(jié)果相符。另外，最終模擬運(yùn)動(dòng)距離為750 m，與實(shí)際運(yùn)動(dòng)距離770 m相差不大。

對(duì)于風(fēng)險(xiǎn)評(píng)估，不僅要關(guān)心滑坡是否會(huì)失穩(wěn)，也應(yīng)該關(guān)心滑坡失穩(wěn)后會(huì)造成什么結(jié)果，即關(guān)心滑坡失穩(wěn)后潛在的影響范圍?；逻\(yùn)動(dòng)過(guò)程分析非常關(guān)鍵，比如，重大交通工程遇到滑坡時(shí)，如果能夠?qū)υ摶逻\(yùn)動(dòng)距離進(jìn)行預(yù)測(cè)，了解潛在的影響范圍，可進(jìn)行有效避讓，而不是進(jìn)行大規(guī)模的工程防治，將大大節(jié)約成本。因此，未來(lái)的研究工作需要集中 在如何對(duì)滑坡運(yùn)動(dòng)過(guò)程進(jìn)行合理預(yù)測(cè)方面。

圖1 滑坡運(yùn)動(dòng)全過(guò)程數(shù)值模擬

4 發(fā)展趨勢(shì)

目前滑坡災(zāi)害的風(fēng)險(xiǎn)定量評(píng)價(jià)主要是通過(guò)搜集歷史數(shù)據(jù)，采用統(tǒng)計(jì)方法實(shí)現(xiàn)對(duì)危險(xiǎn)性、易損性取值量化。存在的主要問(wèn)題是對(duì)滑坡災(zāi)害的危險(xiǎn)性以及易損性評(píng)價(jià)，較少考慮災(zāi)害自身的物理機(jī)制。因此在未來(lái)的風(fēng)險(xiǎn)定量評(píng)價(jià)中，應(yīng)更多考慮災(zāi)害自身的物理機(jī)制量化滑坡災(zāi)害的風(fēng)險(xiǎn)。另外，將風(fēng)險(xiǎn)評(píng)估結(jié)果用于指導(dǎo)變形監(jiān)測(cè)、滑坡失穩(wěn)后潛在威脅范圍的劃定、滑坡災(zāi)害預(yù)警及應(yīng)急預(yù)案的制定等，以及如何將單體滑坡災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)與風(fēng)險(xiǎn)管理相結(jié)合，制定相關(guān)災(zāi)害管理辦法與流程。

5 結(jié)語(yǔ)

滑坡定量風(fēng)險(xiǎn)評(píng)估方法將不確定性與受到損失有機(jī)結(jié)合，為風(fēng)險(xiǎn)量化提供了一種合理有效的手段。評(píng)述了單體滑坡定量風(fēng)險(xiǎn)評(píng)價(jià)的關(guān)鍵組成部分（滑坡發(fā)生概率評(píng)價(jià)、滑坡運(yùn)動(dòng)過(guò)程評(píng)價(jià)以及易損性評(píng)價(jià)），為工程人員在實(shí)際工程中對(duì)滑坡風(fēng)險(xiǎn)管理提供依據(jù)。