考慮空間變異特性的邊坡可靠度分析簡化方法

董少亮　李　輝　朱營營

(青建集團(tuán)股份公司，山東青島　266071)

?

考慮空間變異特性的邊坡可靠度分析簡化方法

董少亮李輝朱營營

(青建集團(tuán)股份公司，山東青島266071)

摘要:介紹了一種衡量邊坡穩(wěn)定性的簡化處理方法，采用蒙特卡洛法與極限平衡法，對某均質(zhì)邊坡進(jìn)行了可靠度分析，并將簡化方法的結(jié)果與已有結(jié)果進(jìn)行了對比，結(jié)果表明:失效概率誤差控制在20%以內(nèi)，可以嘗試采用簡化方法對邊坡進(jìn)行可靠度分析。關(guān)鍵詞:邊坡，可靠度，土體參數(shù)，空間變異性

0　引言

在過去的幾十年里邊坡穩(wěn)定分析引起了相當(dāng)多的關(guān)注。以往的邊坡穩(wěn)定確定性分析中，一般采用譬如Spencer法、Bishop法、Morgenstern-Price法、SARMA法、不平衡推力法等計算邊坡的整體穩(wěn)定性，該方法通常將同一個土層的土體參數(shù)看作是確定值，進(jìn)而確定唯一的安全系數(shù)，所以確定性分析得到的安全系數(shù)并不能客觀地反映邊坡的真實安全程度［1-3］。隨后人們又提出了可靠度分析方法，該方法以“可靠度指標(biāo)”或者邊坡“失效概率”來衡量邊坡穩(wěn)定性。一般采用矩方法、蒙特卡羅模擬等。其中蒙特卡洛法是一種依據(jù)統(tǒng)計抽樣理論，利用計算機研究隨機變量的數(shù)值計算方法，該法受問題條件限制的影響較小，應(yīng)用越來越廣泛，在目前是一種相對精確的方法［4-10］。然而，常規(guī)的可靠度分析中，通常將土層的參數(shù)值看作是符合一定分布規(guī)律的隨機變量，土層內(nèi)各點的參數(shù)是一致的，忽略了土體具有空間變異性這一重要特性，也即:沒有考慮同一土層不同空間位置上土體參數(shù)的相關(guān)性和差異性。

為了分析土體參數(shù)的空間變異特性對邊坡破壞概率的影響，本文提出一種簡化的處理方法。該法根據(jù)土層的厚度與土體材料的相關(guān)長度的比值，將土層細(xì)化為一定數(shù)量的獨立土層。本文擬利用蒙特卡洛法和極限平衡理論對某均質(zhì)邊坡進(jìn)行系統(tǒng)可靠度分析，與已有結(jié)果進(jìn)行對比分析，希望得到有益結(jié)論，對邊坡可靠度分析的研究提供一定的參考。

1　考慮土體參數(shù)空間變異性的可靠度分析

1．1蒙特卡洛模擬

本文擬采用蒙特卡洛法［11］進(jìn)行邊坡可靠度分析，其中粘聚力和摩擦角c，φ是影響邊坡穩(wěn)定性的兩個關(guān)鍵性指標(biāo)。功能函數(shù)F為隨機變量c，φ的函數(shù)［12］。即:

其中，F(xiàn)s為安全系數(shù)，可以用Bishop法、Spencer法等計算得到。將隨機生成的N個樣本，依次代入功能函數(shù)F，可得到N個功能函數(shù)值，統(tǒng)計功能函數(shù)F＜0的個數(shù)，記為M。那么邊坡的失效概率大約為:

當(dāng)N取值足夠大時，得到的頻率即可認(rèn)為是邊坡的失效概率，但在實際應(yīng)用中，N通常取一大值，具體取值與要分析的失效概率成反比。

1．2土體參數(shù)空間變異性及簡化方法

土體相關(guān)系數(shù)公式如下［13］:

其中，Rij為相關(guān)系數(shù);zi-zj為土體中的兩個點的距離;λ為相關(guān)長度。對于不同的土體有不同的相關(guān)長度，這個長度值就相當(dāng)于土體中兩點相關(guān)性的一個臨界值，當(dāng)土體中兩點之間的距離小于這個臨界值時，這兩點之間的物理屬性相關(guān)性較強;當(dāng)土體中兩點之間的距離大于這個臨界值時，這兩點之間的物理屬性相關(guān)性較弱。

由式(3)可知，如果兩點間的距離超過了相關(guān)長度，則兩點間的相關(guān)系數(shù)就會急劇下降，近似地作為不相關(guān)處理;如果兩點間的距離小于相關(guān)長度，則兩點間的相關(guān)系數(shù)增大，兩點高度相關(guān)。因此，我們根據(jù)土層的厚度與土體材料的相關(guān)長度的比值，將土層細(xì)化為一定數(shù)量得到獨立土層，進(jìn)而進(jìn)行可靠度分析。

2　算例分析

該算例采用拉丁超立方抽樣方法，隨機抽取1 000個樣本，利用Bishop法和Spencer法計算安全系數(shù)Fs。對每次抽樣，均進(jìn)行最小安全系數(shù)搜索，然后將安全系數(shù)Fs代入功能函數(shù)進(jìn)行統(tǒng)計計算。其中Bishop法和Spencer法都選擇正態(tài)分布、圓弧滑動面，

考慮文獻(xiàn)［13］中的算例，該邊坡為均質(zhì)粘性邊坡，邊坡剖面圖如圖1所示，土體的密度為2 000 kg/m3，粘聚力c的均值為40 kPa，標(biāo)準(zhǔn)差為10 kPa，摩擦角為0°。

圖1　算例邊坡剖面圖

表1　計算結(jié)果匯總表

依次假定相關(guān)長度為1 m，2 m，2．2 m，2．5 m，2．9 m，3．3 m，4 m，5 m，6．7 m，10 m，20 m，相應(yīng)的獨立土層數(shù)目為20，10，9，8，7，6，5，4，3，2，1。用SLIDE軟件模擬的結(jié)果如表1所示。

將本文用蒙特卡洛模擬在SLIDE軟件中實現(xiàn)的簡化方法結(jié)果與文獻(xiàn)［13］中結(jié)果對比，如圖2所示。

圖2　失效概率折線圖

本文簡化方法結(jié)果與文獻(xiàn)［13］中結(jié)果的相對誤差如圖3所示。

圖3　失效概率誤差圖

由圖3可以發(fā)現(xiàn)，本文簡化方法結(jié)果與文獻(xiàn)［13］隨機場模擬的結(jié)果的相對誤差控制在20%以內(nèi)。

圖4給出了不同相關(guān)長度下，系統(tǒng)可靠度分析需要考慮的滑動面的個數(shù)。

圖4　滑動面?zhèn)€數(shù)

由圖4可見，當(dāng)相關(guān)長度較小時，需要考慮的危險滑動面?zhèn)€數(shù)較多，隨著相關(guān)長度的增加，需要考慮的危險滑動面?zhèn)€數(shù)減少，這說明，當(dāng)土體材料相關(guān)性較強時(相關(guān)長度較小時)，系統(tǒng)可靠度分析需要考慮較多的危險滑動面，然而隨著土體材料相關(guān)性的減弱(即相關(guān)長度變大)，系統(tǒng)可靠度分析所需的危險滑動面較少，這與文獻(xiàn)［8］中的結(jié)論一致。

3　結(jié)語

本文采用蒙特卡洛模擬方法，利用SLIDE軟件，采用簡化方法對某粘性土質(zhì)邊坡進(jìn)行了考慮土體參數(shù)空間變異性的可靠度分析。結(jié)果表明:簡化方法與文獻(xiàn)［13］中算例相比，失效概率誤差控制在20%以內(nèi)，隨著相關(guān)長度的增加，邊坡系統(tǒng)可靠度分析中需要考慮的危險滑動面逐漸減少。簡化方法可以嘗試用于考慮土體參數(shù)空間變異性的邊坡系統(tǒng)可靠度分析。

參考文獻(xiàn):

［1］熊贊民，高全臣，王春來．SLIDE在深基坑支護(hù)可靠性分析中的應(yīng)用［J］．礦業(yè)研究與開發(fā)，2008，28(2):39-41．

［2］任芳．基于蒙特卡洛法的邊坡可靠度分析［J］．福建建科科技，2010(6):13-14．

［3］祝玉學(xué)．邊坡可靠性分析［M］．北京:冶金工業(yè)出版社，1993．

［4］李亮，褚雪松，于廣明．蒙特卡羅法在邊坡失效概率計算中的應(yīng)用［J］．青島理工大學(xué)學(xué)報，2012，33(6):6-9．

［5］呂惠，李亮，路世豹，等．土坡系統(tǒng)可靠度分析［J］．山西建筑，2013，39(22):47-48．

［6］褚雪松，李亮，呂慧．考慮土體材料二維空間變異特性的邊坡可靠度分析［J］．水利水電科技進(jìn)展，2015，35(2): 32-36．

［7］李亮，褚雪松，鄭榕明．Rosenblueth法在邊坡可靠度分析中的應(yīng)用［J］．水利水電科技進(jìn)展(CSCD)，2012，32(3):53-56．

［8］褚雪松，李亮．偽蒙特卡羅法及其在邊坡可靠度分析中的應(yīng)用［J］．土木建筑與環(huán)境工程，2013，35(6):33-39．

［9］褚雪松，王旭春，張勇強，等．基于代表性滑動面的邊坡系統(tǒng)可靠度分析［J］．煤炭學(xué)報，2014，39(6):1077-1083．

［10］李亮，褚雪松，袁長豐．快速蒙特卡羅法及其在土坡可靠度分析中的應(yīng)用［J］．煤炭學(xué)報，2013，38(9):1576-1582．

［11］陳祖煜．土質(zhì)邊坡穩(wěn)定分析——原理、方法、程序［M］．北京:中國水利水電出版社，2010．

［12］楊繼紅，劉漢東，秦四清，等．考慮土性參數(shù)空間變異性的邊坡可靠度分析［J］．工程地質(zhì)學(xué)報，2007，15(2):205-207．

［13］Wang Yu，Cao Zijun，Au Siu-Kui．Practical reliability analysis of slope stability by advanced Monto Carlo simulations in a spreadsheet［J］．Canadian Geotechnical Journal，2011，48(1): 162-172．

Simplified approach to reliability analysis of a soil slope stability taking into consideration of spatial variability

Dong ShaoliangLi HuiZhu Yingying
(Qingjian Group Corporation Ltd，Qingdao 266071，China)

Abstract:Introduces the simplified treatment for slope stability，Monte-Carlo method and limit equilibrium method is used to perform the reliability analysis of homogeneous soil slope．The obtained results are compared with those from the existing results．The comparison shows that the simplified method can yield the result within 20%error，so the simplified method can be used for slope reliability analysis．

Key words:slope，reliability，soil parameters，spatial variability

作者簡介:董少亮(1984-)，男，工程師

收稿日期:2015-10-27

文章編號:1009-6825(2016)01-0079-02

中圖分類號:TU413．62

文獻(xiàn)標(biāo)識碼:A