王啟亮，孟朝霞

(1.山西水利職業(yè)技術(shù)學(xué)院，山西運(yùn)城 044004;2.運(yùn)城學(xué)院，山西運(yùn)城 044004)

地震滑坡風(fēng)險(xiǎn)分析研究

王啟亮1，孟朝霞2

(1.山西水利職業(yè)技術(shù)學(xué)院，山西運(yùn)城 044004;2.運(yùn)城學(xué)院，山西運(yùn)城 044004)

為定量評價(jià)地震滑坡風(fēng)險(xiǎn)管理，對地震滑坡風(fēng)險(xiǎn)分析各因子取值方法進(jìn)行了量化分析。針對承災(zāi)體在坡下及距離災(zāi)害體有一定距離的兩種情況分別討論。重點(diǎn)對一次地震過程中滑坡潛在滑移距離及擴(kuò)展范圍進(jìn)行了計(jì)算，從而對風(fēng)險(xiǎn)分析中滑坡到達(dá)承災(zāi)體概率這一因子進(jìn)行了量化。研究成果可為地震滑坡風(fēng)險(xiǎn)管理的量化分析提供理論依據(jù)，為滑坡風(fēng)險(xiǎn)管理由定性評價(jià)向定量評價(jià)發(fā)展提供參考。

風(fēng)險(xiǎn)管理;承災(zāi)體;地震滑坡;定量風(fēng)險(xiǎn)分析

0 引言

國內(nèi)外研究學(xué)者對滑坡風(fēng)險(xiǎn)進(jìn)行了大量的定性分析及少量量化分析的探討［1－4］，但對于地震滑坡的風(fēng)險(xiǎn)分析研究還較少。尤其是2008年“5·12”汶川地震后誘發(fā)的大量滑坡等地震地質(zhì)災(zāi)害，不僅對人民生命財(cái)產(chǎn)造成重大損失，而且對震前地質(zhì)災(zāi)害的預(yù)防尤其是震前滑坡地質(zhì)災(zāi)害風(fēng)險(xiǎn)評估提出了挑戰(zhàn)。由于目前中國尚無較成熟的滑坡風(fēng)險(xiǎn)管理指南性文件出現(xiàn)，以指導(dǎo)滑坡尤其是地震滑坡風(fēng)險(xiǎn)管理。本文對地震滑坡風(fēng)險(xiǎn)的量化分析進(jìn)行了探討，并以Newmark假設(shè)為基礎(chǔ)，重點(diǎn)對滑坡到達(dá)承災(zāi)體的概率取值方法進(jìn)行了探討。不僅能為滑坡災(zāi)害的防災(zāi)減災(zāi)提供參考，且能為地震滑坡的風(fēng)險(xiǎn)管理提供理論支持。

1 風(fēng)險(xiǎn)分析方法

滑坡災(zāi)害風(fēng)險(xiǎn)管理包括危險(xiǎn)性分析、危害性分析、風(fēng)險(xiǎn)分析、風(fēng)險(xiǎn)評價(jià)和風(fēng)險(xiǎn)管理等5部分［5］。其中風(fēng)險(xiǎn)分析是風(fēng)險(xiǎn)管理中最重要的組成部分。它分為定性分析和定量分析兩種分析方法，涉及到風(fēng)險(xiǎn)源的識別、災(zāi)害發(fā)生概率估計(jì)、承災(zāi)體的確定及其易損性估計(jì)、風(fēng)險(xiǎn)計(jì)算等過程。

2005年5 ～6月加拿大溫哥華“滑坡災(zāi)害風(fēng)險(xiǎn)管理國際會議”確定了風(fēng)險(xiǎn)計(jì)算的公式如下:

式中:PLOL——人口傷亡年概率;

PL——滑坡發(fā)生的概率;

PT:L——滑坡到達(dá)承災(zāi)體的概率;

PS:T——承災(zāi)體的時(shí)空概率;

VD:T、Vprop:S——承災(zāi)體對于滑坡事件的易損性;

Rprop——財(cái)產(chǎn)的年度風(fēng)險(xiǎn);

E——承災(zāi)體的價(jià)值。

當(dāng)涉及到多個(gè)滑坡，且滑坡災(zāi)害之間是相互獨(dú)立的個(gè)體時(shí)，其計(jì)算公式可表述為:

式中:n——滑坡災(zāi)害的數(shù)量。

2 風(fēng)險(xiǎn)分析因子的量化分析

式(1)和式(2)中各因子的量化取值方法如下:

(1)災(zāi)害概率(PL)根據(jù)n年內(nèi)發(fā)生災(zāi)害的概率來取值。n取值依據(jù)為政府及業(yè)主對特定承災(zāi)體在滑坡災(zāi)害威脅下的時(shí)間要求。對單個(gè)災(zāi)害體來說，發(fā)生概率為:PL=1/n;

(2)時(shí)空概率(PS:T)是一個(gè)條件概率。是在災(zāi)害發(fā)生時(shí)承災(zāi)體在災(zāi)害影響區(qū)內(nèi)的概率，其值介于0和1之間。對于承災(zāi)體為建筑物而言，若其在滑坡體上或在其滑動路徑上，則其時(shí)空概率值為1;若對承災(zāi)體為人員狀況而言，時(shí)空概率指的是一段時(shí)間內(nèi)(通常以一年作為計(jì)量單位)人員處于建筑物中的時(shí)間。

(3)易損性(VD:T和Vprop:s)是在災(zāi)害影響區(qū)內(nèi)，一個(gè)或多個(gè)對象受到災(zāi)害影響或破壞的程度。它也是個(gè)條件概率，前提是滑坡災(zāi)害發(fā)生，且承災(zāi)體位于滑坡體上或滑坡路徑上。對于承災(zāi)體為財(cái)產(chǎn)而言，通常用0－1來描述承災(zāi)體受損的程度;對承災(zāi)體為人而言，通常描述為人處于滑坡體上或其滑動路徑上的概率，亦用0－1來表示。易損性評價(jià)分為面評價(jià)和點(diǎn)評價(jià)。面評價(jià)是指采用區(qū)域承災(zāi)體綜合指數(shù)來反映該地區(qū)總易損性。對于地震滑坡而言，同樣可用人口密度、單位面積工農(nóng)業(yè)產(chǎn)值、交通干線密度以及土地資源豐度來評價(jià)易損性因子。汶川地震形成的以縣域或鄉(xiāng)鎮(zhèn)為單位的承災(zāi)體易損性評價(jià)即可采用面評價(jià)的方法，評價(jià)方法可參考文獻(xiàn)2。點(diǎn)評價(jià)是針對單一災(zāi)害體對特定承災(zāi)體造成的易損性的評價(jià)，評價(jià)時(shí)可根據(jù)滑坡災(zāi)害體的體積，及它們撞擊房屋時(shí)的可能速度。根據(jù)專家打分的形式量化災(zāi)害體撞擊房屋后房屋內(nèi)人員傷亡的概率［3］。

(4)承災(zāi)體價(jià)值(E)可根據(jù)實(shí)體的市場價(jià)格予以估算。特殊承災(zāi)體價(jià)值計(jì)算方式可參考文獻(xiàn)4。

(5)災(zāi)害到達(dá)承災(zāi)體的概率PT:L取值取決于災(zāi)害體與承災(zāi)體各自的位置及災(zāi)害體可能的運(yùn)動路徑。當(dāng)滑坡災(zāi)害發(fā)生破壞建筑物且人員未能幸免于難時(shí)，其取值標(biāo)準(zhǔn)如表1［5］。

當(dāng)滑坡距離承災(zāi)體(人員及建筑物)較遠(yuǎn)時(shí)，在地震荷載作用下滑坡到達(dá)承災(zāi)體概率的計(jì)算主要基于Newmark假設(shè)，其計(jì)算方式如下。

表1 坡底建筑物受到破壞的概率Table 1 Destructive probability of buildings with a slope behind

假設(shè):①滑坡塊體為剛性摩擦塊體;②一次地震坡體水平波動震蕩形成永久位移的時(shí)間占地震總持時(shí)的1/3;③臨界平衡狀態(tài)為穩(wěn)定系數(shù)Fs=1;④不考慮動孔隙水壓力作用;⑤臨界加速度值不隨地震作用而改變。

地震過程中，動荷作用會增加臨空面方向的負(fù)荷，下滑力隨時(shí)間增加的量(ΔN)會消除靜態(tài)狀態(tài)下下滑力N與最大抗滑力Rmax之間的差值，并致使滑坡塊體達(dá)到短暫的臨界平衡狀態(tài)。此時(shí)ΔN可表示為質(zhì)量與坡體表面地震加速度的乘積，即ΔN=mac。為簡化計(jì)算，可用穩(wěn)定系數(shù)來表述ac，即:

當(dāng)滑面為直線型時(shí)，N=mgsinα

式中:g——重力加速度;

α——滑面傾角;

m——塊體質(zhì)量。

故可得

式中:Fs——穩(wěn)定系數(shù)。

當(dāng)滑塊在地震動荷載作用下地表運(yùn)動時(shí)，若塊體相對坡面而言無位移，則在某一時(shí)刻t時(shí)原來在坡面A點(diǎn)的滑塊應(yīng)到達(dá)B點(diǎn)，由于慣性作用，滑塊將產(chǎn)生一個(gè)相對坡面的位移U(t)，此時(shí)滑塊位于C點(diǎn)(圖1)。在計(jì)算動荷下滑坡體位移時(shí)，需要確定的一個(gè)重要參數(shù)就是地面加速度時(shí)間函數(shù)a(t)(圖2)。

圖1 滑塊滑動示意圖Fig 1Sketch of slip mass

圖2 地震滑坡位移計(jì)算示意圖Fig 2Sketch for displacement calculating

由圖2可見，滑坡在動荷作用下的運(yùn)動狀態(tài)可分為未滑動(0－tA)、開始滑動(此時(shí)t=tA)、滑動階段(tA－tC)及停止滑動(t=tC)幾個(gè)時(shí)段，當(dāng)a(t)＜aC時(shí)，滑塊處于穩(wěn)定階段;a(t)＞aC時(shí)，滑塊具有剩余加速度Δa，此時(shí)Δa=a(t)－aC，隨地震波的傳播過程，Δa逐漸增大至Δamax，隨后減至aC以下(圖2中B點(diǎn))。此時(shí)滑塊速度可用塊體相對加速度Δa通過對時(shí)間積分的形式計(jì)算得到，即B點(diǎn)塊體速度為:

其位移可相應(yīng)求得［6］:

為簡化計(jì)算，假設(shè)地面震動加速度a= amaxsinωt，其中amax為地面最大加速度，ω為圓頻率，此時(shí)的臨界加速度ac=g(Fs－1)sinα，根據(jù)Newmark假設(shè)，可推導(dǎo)得到一次地震由于坡體水平波動振蕩產(chǎn)生總位移計(jì)算公式如下［6］:

由上式可以看出，地震持時(shí)越長，地震周期越長，其位移量越大。承災(zāi)體距離滑坡越近，災(zāi)害到達(dá)承災(zāi)體的概率越大。其易損性值越趨向于1。當(dāng)滑坡距離承災(zāi)體(人員及建筑物)較遠(yuǎn)時(shí)，在地震荷載作用下滑坡到達(dá)承災(zāi)體的概率，可依據(jù)位移計(jì)算公式予以計(jì)算后，進(jìn)行量化分析。

3 結(jié)論

充分的滑坡定量風(fēng)險(xiǎn)分析，目前還存在一些難點(diǎn)。如對滑坡災(zāi)害特征的把握具有很大的主觀性;滑坡災(zāi)害發(fā)生概率分析具有不確定性;承災(zāi)體時(shí)空概率的確定具有隨機(jī)性等因素，是導(dǎo)致目前滑坡風(fēng)險(xiǎn)定性分析較成熟，而定量分析較難推廣的主要原因。近幾年地震滑坡的大量出現(xiàn)已經(jīng)造成了重大的人員及財(cái)產(chǎn)損失。本文通過對地震滑坡風(fēng)險(xiǎn)分析各因子進(jìn)行量化分析研究，希望能促進(jìn)滑坡風(fēng)險(xiǎn)分析由定性研究向定量研究拓展。同時(shí)，在地震多發(fā)區(qū)于地震前進(jìn)行滑坡災(zāi)害風(fēng)險(xiǎn)分析，結(jié)合地質(zhì)災(zāi)害防治工程便易做到防患于未然，大大降低地震災(zāi)害由于滑坡造成的人員傷亡及財(cái)產(chǎn)損失。

［1］Fell·R，Ho·K·K·S，Lacasse·S，et al..A Framework for Landslide Risk Assessment and Management［C］// Hunger·O，F(xiàn)ell·R，Couture·R，et al..Proceedings of theInternationalConferenceonLandslideRisk Management.Vancouver，Canada，31May－3June，2005.Taylor and Francis，London:3－25.

［2］張梁.地質(zhì)災(zāi)害風(fēng)險(xiǎn)評價(jià)理論與方法［J］.中國地質(zhì)礦產(chǎn)經(jīng)濟(jì)，1996，(4):43－44.

［3］王樹豐，張茂省，唐亞明，等.延安寶塔山景區(qū)滑坡地質(zhì)災(zāi)害風(fēng)險(xiǎn)評估［J］.工程地質(zhì)學(xué)報(bào)，2009，17(5):628－631.

［4］王樹豐，張茂省，陳志新，等.基于滑坡風(fēng)險(xiǎn)管理的寶塔山景區(qū)景點(diǎn)價(jià)值核算［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2009，20(2):28－30.

［5］張茂省，唐亞明.地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查方法與實(shí)踐［J］.地質(zhì)通報(bào)，2008，27(8):1206－1207.

［6］毛彥龍，胡廣韜，趙法鎖，等.地震動觸發(fā)滑坡滑動的機(jī)理［J］.西安工程學(xué)院學(xué)報(bào)，1998，20(4):45－48.

Abstract:For quantitatively assessing risk management of a slope under earthquake.Impact，quantitative analysis was studied by choosing factors of slope risk management.For two scenarios，objects just under the slop and objects have a distance to the slop.Sepecially，sliding distance and its incidence during one earthquake were calculated.Then the factors of probability for the slope reaching the objects was quantitatively assessed.Results of the discussion can provide not only theoretical foundation for quantitative analysis for risk management of a slope under earthquake impact，but reference for the process from a qualitative evaluation to a quantitative assessment for risk management. Key words:risk management;objects at risk;slope under earthquake;QRA

Risk analysis for slopes induced under by earthquake impact

WANG Qi－liang1，MENG Zhao－xia2
(1.Shanxi Water Conservation Professional Technology Institute，Yuncheng044004，China; 2.Yuncheng University，Yuncheng044004)

1003－8035(2010)03－0014－03

P642

A

2010－05－11;

2010－06－16

國土資源大調(diào)查項(xiàng)目(1212010740907)

王啟亮(1964—)，男，山西臨猗人，副教授，研究方向?yàn)榄h(huán)境地質(zhì)。

E－mail:wql.976@163.com