曹羽哲, 曹運江

(湖南科技大學, 湖南 湘潭 411201)

降雨型滑坡的發(fā)生頻率較高，分布范圍非常廣，斜坡破壞形式多種多樣，嚴重制約災害多發(fā)地區(qū)的社會經(jīng)濟發(fā)展，威脅災害區(qū)人民生命財產(chǎn)安全。因此，對區(qū)域開展降雨型滑坡災害穩(wěn)定性評價和滑坡預測預報，對國土規(guī)劃、防災減災、災害管理與決策等均具有非常重要的研究價值和現(xiàn)實意義[1]。

滑坡穩(wěn)定性分析是提高滑坡災害預警預報精度的先決條件。目前評價滑坡穩(wěn)定性所使用最多的方法是極限平衡法，如簡化Bishop法和瑞典條分法等，這些定值分析方法雖然能有效得出邊坡穩(wěn)定系數(shù)，但忽略了實際存在的不確定性因素及計算參數(shù)的隨機分布對滑坡產(chǎn)生的影響。大量的研究成果表明，滑坡災害是一個非常復雜的非線性系統(tǒng)，在不同區(qū)域、不同構造背景下，各影響因素對滑坡發(fā)生的促進或抑制作用及程度各不相同。因此，需要分析對滑坡的發(fā)生具有顯著影響作用的工況，建立合理、符合實際的地質(zhì)災害分析模型，以獲得較好的預測結果。由于在實際工程中，邊坡穩(wěn)定性計算參數(shù)分布較為離散，可以視為隨機變量，故客觀上的穩(wěn)定性系數(shù)亦為隨機變量，因而采用概率分析方法進行邊坡穩(wěn)定性評價更為合理。早在20世紀50年代，F(xiàn)reudental就提出將概率分析方法引入工程結構可靠性分析中[2]；通過不斷地學習探索，20世紀80年代，祝玉學將概率分析方法應用到土坡穩(wěn)定性分析中并出版了中國第一部相關著作《邊坡可靠性分析》[3]，從此，概率分析也逐漸成為一種熱門的穩(wěn)定性評價手段；對于滑坡發(fā)生概率，眾多學者采用了多種方法及概率模型進行了分析；王志恒等采用確定性概率系數(shù)模型來對高程、坡度、巖土類型及斷裂構造等降雨型滑坡孕災因子進行分析[4]。喬建平等提出了將滑坡穩(wěn)定性分析與滑坡發(fā)育率相結合的KP模型計算滑坡破壞概率的方法[5]。蒲書豪等提出將敏感性分析與概率分析法相結合的分析手段，以敏感性分析結果的針對性來彌補概率分析法參數(shù)選取上的不足[6]，以期得到更具有針對性、高效的計算結果，更好地為滑坡防治工程提供參考。杜時貴等采用破壞概率分析法，結合不平衡推力傳遞法，對巖石邊坡的穩(wěn)定性進行了評價[7]。

為更好地研究分降雨型滑坡發(fā)生破壞的概率，本文采用破壞概率分析法，考慮到穩(wěn)定性分析傳統(tǒng)方法的合理性和局限性，結合穩(wěn)定性分析中的剩余推力法，以韶關地區(qū)典型的降雨型滑坡為研究對象，對其天然工況和暴雨工況進行穩(wěn)定性分析與評價，以期更好地確定滑坡破壞概率，為滑坡預防治理工程提供參考。

1 研究區(qū)概況

韶關市地處廣東省北部，丘陵山地多，地勢起伏大，地形切割強烈，地質(zhì)構造較為復雜，地質(zhì)環(huán)境脆弱，是廣東省地質(zhì)災害多發(fā)地區(qū)之一。滑坡區(qū)位于韶關市甘棠工業(yè)園區(qū)西側(cè)的山坡上，邊坡主要為人類工程活動切坡形成的人工邊坡，由于人類活動加劇導致附近原有地質(zhì)環(huán)境發(fā)生改變，邊坡中一部分坡體相對于另一部份坡體產(chǎn)生相對位移導致其喪失原有穩(wěn)定性，從而形成滑坡。一般坡度為20°～40°。區(qū)域上有褶皺、斷裂構造發(fā)育。結合野外地質(zhì)災害調(diào)查和鉆孔揭露資料，研究區(qū)山體淺表層的殘積土層較厚，地表局部巖石裸露。研究區(qū)出露地層為由上至下為第四系(Q4 el+dl)殘坡積層，下石炭統(tǒng)梓門橋組(C1z)及測水組(C1c2)，巖層產(chǎn)狀65°∠45°?；轮骰瑒臃较?30°，長約91 m，寬約238 m，前后緣高差約30 m，滑坡后緣陡坎約1.1 m。該滑坡抗震設防烈度為Ⅵ度，地震動峰值加速度值為0.05 g，場區(qū)地震動反映譜特征周期為0.35 s，區(qū)域地殼穩(wěn)定?；麦w全貌見圖1。

圖1 韶關市工業(yè)園西側(cè)滑坡全貌

根據(jù)鉆探取心情況以及滑坡的已有變形破壞特征和滑坡微地貌形態(tài)綜合分析，滑坡體物質(zhì)組成為殘坡積土、全風化泥巖，其平均厚度約為5 m，體積約6.8萬m3，屬淺層大型土質(zhì)滑坡，滑動面應位于全風化泥巖層中，滑動面平均深度4～6 m。降雨為主要誘發(fā)滑坡破壞的因素，滑坡觸發(fā)原因總體歸納為由于其處于地形地貌環(huán)境的特定性結合相應的地層巖土體環(huán)境下，力學性質(zhì)較差的坡體再經(jīng)降雨入滲及人類活動影響其受力，導致其穩(wěn)定性進一步變差。由于滑坡下方為韶關市甘棠工業(yè)園廠區(qū)，若滑坡發(fā)生蠕滑破壞，不僅會對工業(yè)園廠區(qū)造成直接破壞，且也會影響到擬建工廠的建設，威脅人數(shù)1 037人，造成的潛在經(jīng)濟破壞或損失大于3 000萬元，采取及時合理的工程防治措施刻不容緩。

在對韶關市工業(yè)園西側(cè)滑坡進行治理前，需要對滑坡穩(wěn)定性進行分析評價，本文通過建立降雨型滑坡穩(wěn)定性系數(shù)計算模型，并采用破壞概率分析法來對滑坡在兩種不同工況下的破壞概率進行計算，從而更準確的評價滑坡穩(wěn)定性。

2 模型及方法

2.1 降雨型滑坡穩(wěn)定性系數(shù)計算模型

對于降雨型滑坡，首先要建立合適的降雨型滑坡穩(wěn)定性系數(shù)計算模型，以便更精準的分析與評價降雨型滑坡穩(wěn)定性。

本文沒有選取邊坡工程常用的簡化Bishop法及瑞典條分法，而采用剩余推力法進行穩(wěn)定性計算。對于滑面為任意形狀的滑坡穩(wěn)定性計算，采用剩余推力傳遞法較為合適。其考慮了滑體自重、坡面荷載、動水壓力、靜水壓力、滑動面處的浮托力、暴雨、地震和不同條塊滑面段抗剪強度參數(shù)差異對滑坡穩(wěn)定性的作用和影響，見圖2。

圖2 剩余推力法力學計算模型

剩余推力法計算穩(wěn)定性系數(shù)完整的計算公式如下：

(1)

式中:Ri為第i塊條塊的抗滑力;Ti為第i塊條塊的下滑力;ψj為第i塊的傳遞系數(shù);Rn為最后一個計算條塊的抗滑力;Tn為最后一個計算條塊的下滑力。

由于降雨型滑坡需要著重考慮暴雨工況下坡體受力情況，而在降雨過程中地表水及地下水的活動使得坡體內(nèi)部的水壓力難以精準確定。由此構建降雨型滑坡分析模型，由圖3所示。

從圖3可以看出，當對暴雨條件下的滑梯前緣條塊進行受力計算時，可將ABCDEFGHIA視為一個前緣受靜水壓力阻力作用的整體滑坡，該條塊所受重力應為DEFG內(nèi)所有水、巖重量總和，仿此考慮滑坡塹中受水作用的條塊受力；當水位下降時，孔隙水壓力變?yōu)镵IHGJ，按此水壓線計算滑坡穩(wěn)定性。即無論何種地表水、地下水活動情況，坡體內(nèi)部地下水滲流的作用均可視為作用于各個條塊邊界上的靜水壓力。這種分析模型就可以獲得任意形狀滑動面在復雜荷載作用下的滑坡推力[8]，并且滑動面條分后各條塊的受力均可快速確定，使得降雨型滑坡在暴雨工況下的穩(wěn)定性計算分析變得更為簡潔合理。

圖3 降雨型滑坡分析模型

2.2 滑坡破壞概率計算方法

在對滑坡Fs穩(wěn)定性系數(shù)進行計算時，往往采用的是傳統(tǒng)計算方法，計算滑坡穩(wěn)定性系數(shù)的公式中需要考慮各類因素的具體確定數(shù)值，比如研究區(qū)巖土體力學性質(zhì)參數(shù)C，φ值、各層巖土體的標準重度、條塊劃分的尺寸大小、傾斜角度以及所受各種力的確定值，由此方法計算得出的穩(wěn)定性系數(shù)Fs也為一個具體確定數(shù)值。但在實際工程條件下，雖然條塊劃分尺寸及傾斜角度實為可測定值，但各層巖土體物理力學參數(shù)分布較為離散，并且測試過程中存在誤差，各種受力的大小也會隨著外界條件的改變而改變，導致無法取得其準確數(shù)值并導致計算產(chǎn)生極大的誤差?，F(xiàn)將滑坡穩(wěn)定性系數(shù)帶入概率框架進行分析，穩(wěn)定性系數(shù)計算參數(shù)不為一個定值，而應視為符合某種規(guī)律分布的數(shù)值集合。當確定滑坡穩(wěn)定性系數(shù)的主要影響計算參數(shù)的概率分布，利用Monte-Carlo法，通過產(chǎn)生隨機數(shù)值帶入計算參數(shù)中，使其符合概率分布，將符合概率分布的計算參數(shù)帶入穩(wěn)定性系數(shù)公式計算得出符合概率分布的滑坡穩(wěn)定性系數(shù)[9]。對于此流程不斷重復若干次，計算得出符合概率分布的穩(wěn)定性系數(shù)足夠多時，通過大數(shù)定律可以將此時滑坡穩(wěn)定性系數(shù)屬于破壞的頻率視為破壞概率：

(2)

式中：N為滑坡穩(wěn)定性系數(shù)總體數(shù)量；M為滑坡穩(wěn)定性系數(shù)小于1的數(shù)量。

通過公式(2)，對計算得出的滑坡穩(wěn)定性系數(shù)總體數(shù)量N進行統(tǒng)計分析，便可以驗證滑坡穩(wěn)定性系數(shù)分布規(guī)律是否與主要計算參數(shù)分布規(guī)律一致。由于各個隨機變量對滑坡穩(wěn)定的影響均不相同，因此可以把影響程度較小的量作為定值以簡化計算。在不確定因素中以巖土體黏聚力C和巖土體內(nèi)摩擦角φ對滑坡穩(wěn)定性影響最大[10]，故選取這兩個參數(shù)為隨機變量，其他參數(shù)則作為定值來處理。由此可以建立狀態(tài)方程FS=f(C,φ)。

當確定狀態(tài)方程FS=f(C,φ)之后，下一步需要對巖土體黏聚力C和巖土體內(nèi)摩擦角φ的概率分布進行分析。首先根據(jù)各類工程經(jīng)驗及研究總結廣泛應用的絕大多數(shù)概率分布類型，主要有Bernoulli分布、Poisson分布、極值分布、中心極限分布及其他分布。結合各概率分布類型的特點、物理意義以及主要工程應用，發(fā)現(xiàn)中心極限分布與滑坡穩(wěn)定性主要影響因子分布是較為契合的，并且通過許多學者的研究也表明是合理可行的。對于中心極限分布，首先需要了解中心極限定理，其描述為：如果一個隨機變量由多種(數(shù)量較大)的，并且相互獨立的隨機因素綜合影響所構成，且其中每個隨機因素在總體影響中所起到的作用均十分微小，則此隨機變量往往近似服從正態(tài)分布[11]。其原理如公式(3)所示：

Y=X1+X2+…+Xn

(3)

式中:Y為隨機變量；X1,X2,Xn為隨機變量Y的相互獨立影響因素。由此可以看出,當n→∞時,隨機變量Y服從正態(tài)分布,由此可以看出,若隨機變量的影響因素互相獨立的影響隨機變量的分布,并為相互簡單疊加,即可以視為正態(tài)分布。

若各影響因素之間存在相互影響，不為簡單相加方式，可以將其認為共同聯(lián)合影響，類似于相乘關系，原理公式如式(4)所示：

Y=X1X2…Xn

(4)

對于公式(4)，兩邊同時取自然對數(shù)，可得：

lnY=lnX1+lnX2+…+lnXn

(5)

可以看出當n→∞時，lnY服從正態(tài)分布，故隨機變量Y符合對數(shù)正態(tài)分布。當隨機變量影響因素不僅僅互相獨立，而是共同作用來規(guī)范隨機變量的分布可以將其認為規(guī)律分布符合對數(shù)正態(tài)分布。

在確定好概率分布模型之后，需要確定巖土體黏聚力C和巖土體內(nèi)摩擦角φ的分布是否為服從正態(tài)分布。當確定滑坡巖土體參數(shù)黏聚力C和巖土體內(nèi)摩擦角φ的分布符合正態(tài)分布或?qū)?shù)正態(tài)分布后，隨機數(shù)值的抽樣可以通過以下公式進行，來保證數(shù)值的范圍準確及離散隨機[12]:

(6)

X=σR+μ

(7)

式中:R1，R2為[0,1]區(qū)間的均勻分布隨機數(shù);R為服從標準正態(tài)分布N(0,1)的隨機變量;X為服從均值為μ、方差為σ的正態(tài)分布隨機變量。

綜上所述，破壞概率分析法的基本計算步驟如下：首先根據(jù)現(xiàn)場勘察和巖土體試驗得出多組巖土體黏聚力C和巖土體內(nèi)摩擦角φ數(shù)值，并對其分布進行判別。當確定了C和φ的分布函數(shù)基本滿足正態(tài)分布之后，可利用Monte-Carlo法，采用公式(6)，(7)進行抽樣與變換，產(chǎn)生一系列滿足正態(tài)分布的C，φ隨機變量數(shù)值。隨后將產(chǎn)生的一系列C，φ隨機變量數(shù)值代入公式(1)中，根據(jù)具體工況對滑坡進行條塊劃分。確定其他基礎數(shù)據(jù)后，結合剩余推力法進行計算可得出一系列相應的穩(wěn)定性系數(shù)Fs的分布函數(shù)，將得到的穩(wěn)定性系數(shù)分布函數(shù)帶入到公式(2)中進行計算，得出滑坡破壞概率曲線，統(tǒng)計滑坡穩(wěn)定性系數(shù)小于1時的分布密度，求出其在所計算滑坡穩(wěn)定性系數(shù)中的占比，便可以求出滑坡破壞概率。計算方法流程見圖4。

圖4 降雨型滑坡破壞概率計算流程

采用破壞分析概率模型求出滑坡穩(wěn)定系數(shù)和破壞概率的數(shù)值后，可以參考按破壞概率大小滑坡穩(wěn)定性等級劃分表[13](表1)對滑坡穩(wěn)定性進行評判。

表1 按破壞概率大小邊坡穩(wěn)定性等級劃分

3 結果與分析

3.1 研究區(qū)滑坡天然及暴雨工況破壞概率計算

將采用上述所構建模型及破壞概率分析法對研究區(qū)滑坡進行破壞概率計算及穩(wěn)定性分析評價。根據(jù)室內(nèi)試驗，野外鑒別及現(xiàn)場原位測試，結合滑坡分布位置范圍，測得多組研究區(qū)滑坡在天然狀態(tài)和暴雨飽水狀態(tài)下的主要巖土體參數(shù)C，φ值，并根據(jù)統(tǒng)計獲得滑坡的同一滑帶土層主要巖土體參數(shù)C，φ在天然與暴雨飽水兩種不同工況下的分布，見表2。

表2 不同工況狀態(tài)下主要巖土體參數(shù)表

通過表2可以看出，巖土體參數(shù)C，φ分布基本服從正態(tài)分布，可以采用公式(6—7)進行隨機抽樣。本次分析累計產(chǎn)生500個隨機數(shù)，得到對應的500組巖土體參數(shù)C，φ值。

在確定巖土體參數(shù)分布之后，還需確定研究區(qū)滑坡主滑面，并繪制主滑面剖面圖；然后在剖面圖中將滑動帶劃分成若干相對規(guī)則的條塊，并確定各滑塊的面積s、單位滑面長度L、滑塊與水平面夾角θ以及地下水位線。本次分析將滑動帶分割成10塊相對規(guī)則的塊體，見圖5。在計算之前，條塊基本幾何參數(shù)可以直接通過AutoCAD在計算剖面上標注測量得出，將其作為定值進行運算。

圖5 滑坡滑動面條分結果

在確定好以上參數(shù)后，將天然狀態(tài)和暴雨飽水狀態(tài)下的500組巖土體參數(shù)C，φ值代入公式(1)中進行計算，得到天然狀態(tài)和暴雨飽水狀態(tài)兩種不同工況下的滑坡穩(wěn)定系數(shù)，根據(jù)所得到的穩(wěn)定系數(shù)可統(tǒng)計分析出滑坡穩(wěn)定系數(shù)在天然工況和暴雨工況下的分布圖，結合500組滑坡穩(wěn)定性系數(shù)在各個范圍的分布數(shù)量，可以得出累計分布概率密度曲線，經(jīng)過繪制便可得出滑坡破壞概率曲線圖，通過分析曲線圖，選取滑坡穩(wěn)定性小于1.0的分布頻率，此頻率即為當前工況條件下滑坡的破壞概率。

經(jīng)過計算統(tǒng)計得出天然工況下研究區(qū)穩(wěn)定性系數(shù)及破壞概率曲線，見圖6?？梢钥闯?，在天然工況下計算得出的滑坡穩(wěn)定性系數(shù)也基本滿足正態(tài)分布，經(jīng)統(tǒng)計其均值為1.106，可以認定天然工況下研究區(qū)滑坡穩(wěn)定性系數(shù)為1.106，屬于基本穩(wěn)定狀態(tài)；通過破壞概率曲線得出研究區(qū)滑坡在天然工況下的破壞概率為37.1%，參考按破壞概率大小滑坡穩(wěn)定性等級劃分表認為研究區(qū)滑坡天然狀態(tài)下的穩(wěn)定情況屬中等危險偏低危險。

圖6 天然工況下研究區(qū)滑坡穩(wěn)定性系數(shù)及破壞概率

當研究區(qū)滑坡經(jīng)歷持續(xù)降雨達到飽和狀態(tài)之后，其穩(wěn)定性系數(shù)及破壞概率曲線經(jīng)計算統(tǒng)計見圖7。通過對研究區(qū)滑坡穩(wěn)定性系數(shù)計算統(tǒng)計，發(fā)現(xiàn)暴雨飽和工況下滑坡穩(wěn)定性系數(shù)降低，暴雨工況下穩(wěn)定性系數(shù)降至1.026，滑坡整體處于欠穩(wěn)定狀態(tài)；破壞概率曲線表明，伴隨著穩(wěn)定性系數(shù)的降低，暴雨飽和工況下研究區(qū)滑坡破壞概率增加至55.3%，穩(wěn)定性情況屬中等危險偏高危險。

圖7 暴雨工況下研究區(qū)滑坡穩(wěn)定性系數(shù)及破壞概率

3.2 考慮巖土體變異系數(shù)及偏度系數(shù)下的滑坡破壞概率計算

對于室內(nèi)試驗測定的現(xiàn)場巖土體C，φ值集合，如果不考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)，便可以按上述方法進行正態(tài)分布抽樣，但如果考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)，需要確定其是按正態(tài)分布進行抽樣還是按對數(shù)正態(tài)分布進行抽樣。根據(jù)廖景高等[14]人的研究成果來看，正態(tài)分布及對數(shù)正態(tài)分布特性對比，完全且僅受偏度系數(shù)的控制影響，當隨機變量的變異系數(shù)小于30%并且偏度系數(shù)小于0.025時，其特征選擇服從正態(tài)分布；否則選擇對數(shù)正態(tài)分布。故此次也對天然黏聚力及暴雨飽水黏聚力采取對數(shù)正態(tài)分布抽樣，并對兩種工況進行對比計算不同分布下的滑坡破壞概率。

當考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)后，研究區(qū)滑坡天然狀態(tài)破壞概率為31.3%(圖8A)，暴雨工況下破壞概率為51.6%(圖8B)，與不考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)計算結果相比兩者分別相差15%及7%，并且不考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)計算的破壞概率比考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)破壞概率高。分析其中的緣由，認為是：當考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)，對巖土體參數(shù)數(shù)據(jù)進行對數(shù)正態(tài)分布處理之后，其分布較原來更加集中，分布范圍減小，導致計算的破壞概率較?。坏豢紤]巖土體參數(shù)變異系數(shù)，巖土體參數(shù)按照正態(tài)分布取值時，由于巖土體參數(shù)分布較為離散，其分布范圍較寬，計算的破壞概率會較高。推測兩者計算結果差距范圍在5%～20%。

圖8 考慮變異系數(shù)及偏度系數(shù)下研究區(qū)滑坡破壞概率曲線

3.3 傳統(tǒng)計算方法對比及降雨特征量確定

當采用破壞概率分析法得出滑坡穩(wěn)定性系數(shù)后，再采用傳統(tǒng)的瑞典條分法對研究區(qū)滑坡的穩(wěn)定性系數(shù)進行計算，計算結果及對比見表3?？梢钥闯觯茐母怕史治龇ńY合剩余推力法與瑞典條分法所計算的滑坡穩(wěn)定系數(shù)相差較小，均可以得出該滑坡在天然工況下基本穩(wěn)定，在暴雨工況下欠穩(wěn)定，計算結果與現(xiàn)場調(diào)查研究情況基本一致，由于調(diào)查研究前研究區(qū)存在持續(xù)降雨，滑坡整體受降雨影響飽和度較大，室內(nèi)試驗所測得滑坡巖土體物理力學參數(shù)基本與飽水狀態(tài)一致。坡表受滑動變形影響較松散，處于蠕動變形階段坡腳前緣水溝局部地段受滑坡變形影響造成損毀，并有加劇的趨勢。所以，研究區(qū)滑坡潛在的安全隱患非常突出。

表3 兩種不同計算方法下的研究區(qū)滑坡穩(wěn)定系數(shù)

針對研究區(qū)滑坡的提前預警預報，應該綜合考慮多種降雨情況，比如前期巖土體含水率、降雨過程中的歷史及強度變化及未來降雨特征量[15]。但是存在以下難點：(1) 關于巖土體含水率，僅能測定部分范圍，大范圍定量定點測定巖土體含水率現(xiàn)在暫時無法實現(xiàn)；(2) 降雨過程不僅一直在變化，還需要考慮到前期累計降水量及降雨時長有關，并不是每次降雨對地質(zhì)災害的影響均一致，時間越靠前、越短，影響越?。?3) 未來降雨特征量也需要考慮到累積雨量和降雨強度的變化。因此，考慮到概率定量分析切實可行性，結合韶關市降水歷史資料，在統(tǒng)計分析之后，主要選取韶關市多年平均月降水量[16]作為降雨外因影響因素進行分析。由圖9可知，韶關市4—9月的月均降雨量較大，對降雨型滑坡的影響加劇，可能誘發(fā)滑坡進一步發(fā)育，產(chǎn)生蠕動破壞。在雨季來臨之前，對研究區(qū)滑坡進行合理、有效的防治是十分有必要的。

圖9 韶關多年平均月降雨量柱狀示意圖

4 結 論

(1) 韶關市工業(yè)園西側(cè)滑坡在天然工況下基本穩(wěn)定，滑坡穩(wěn)定性系數(shù)為1.106，破壞概率范圍31.3%至37.1%，處于中等偏低危險等級；在暴雨工況下欠穩(wěn)定，滑坡穩(wěn)定性系數(shù)為1.026，破壞概率范圍為51.6%至55.3%，處于中等偏高危險等級。

(2) 破壞概率分析法結合剩余推力法所建立的滑坡穩(wěn)定性計算模型，與傳統(tǒng)計算方法相比，同樣能夠得到滑坡在不同工況下的穩(wěn)定情況，并且在計算中更為簡便合理。

(3) 與考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)相比，不考慮巖土體參數(shù)變異系數(shù)及偏度系數(shù)計算的滑坡破壞概率相對更高，兩者計算結果差距范圍在5%～20%左右。實際情況下建議根據(jù)巖土體參數(shù)情況，若變異系數(shù)及偏度系數(shù)較大，分布較離散，可以采用兩種方法同步進行計算破壞概率，確定出滑坡的破壞概率范圍，確保分析的合理性及有效性。

得出結論之后，針對此次所做出的結果和方法，提出一些展望：對于暴雨工況，本次只根據(jù)其飽水狀態(tài)的巖土體參數(shù)進行了隨機取樣，但是降雨對滑坡的影響是多樣的，還需要正確建立巖土體降雨入滲模型，得出滑坡入滲過程，并且已合理的方式確定滑坡破壞臨界雨量。在今后需要采用多種方法耦合進行分析，才能更好的分析出滑坡在不同降雨條件下的破壞概率。