強度參數統計特性對邊坡穩定可靠性的影響

桂勇　羅嗣海　鄧通發

摘要：由于對邊坡強度參數統計特性考慮不全面，使得邊坡可靠性指標計算值偏小，對應的失效概率偏大，常常超過10%，設計中難以采納，造成地基可靠度規范的制定和執行進展緩慢。以某高速公路全風化花崗巖土質高邊坡為工程背景，在分析土質參數統計特性的基礎上，采用極限平衡理論和蒙特卡羅模擬法，系統地分析了土質強度參數的均值不定性、變異性、相關性、區間特性和空間變化性等對邊坡穩定可靠性的影響。結果表明：土性參數的各種統計特性對邊坡穩定可靠性均具有不同程度的影響。可靠性指標計算值隨抽樣距離的減小而增加，隨c或φ均值的增大而增加，隨c和φ的變異系數的增大而減小，隨c和φ相關系數絕對值的增大而增加，考慮區間性后可靠性指標計算值也明顯變大。即變異性對邊坡穩定不利，而相關性、區間性和空間變化性對邊坡穩定有利。因此，準確而全面地考慮土質參數的統計特性，尤其是在參數變異性和相關性的基礎上加以考慮土性參數空間變化性和區間性會更加符合工程實際，且計算結果趨于安全，有利于地基可靠度規范的推廣運用和邊坡工程的安全評價。

關鍵詞：邊坡穩定可靠性；強度參數；統計特性；區間性；空間變化性

中圖分類號：TU 443 文獻標志碼：A 文章編號：1674-4764（2015）04-0067-10

Abstract： Without comprehensive consideration of the strength parameters statistical properties， the slope reliability index calculation value become smaller and the failure probability become larger， often more than 10%， which was difficult to adopt in design.Hence， the foundation reliability norms formulate and implement made slow progress. Take a completely weathered granite soil slope of an expressway as engineering background， and based on the analysis of the statistical characteristics of soil parameters， the limit equilibrium theory and monte carlo simulation method were used to systematically analyze the soil strength parameters mean uncertainty， variability， correlation and interval characteristics and spatial variability of impact on the reliability of the slope stability. The results showed that all kinds of statistical features of soil parameters had varying degrees impact on slope stability reliability. The reliability index calculation value increased with the decrease of the sampling distance， increased with the increase of the mean values of

c or φ， decreased with the increase of the variation coefficientc and φ， increased with the increase of the correlation coefficient absolute value of c and φ， and the reliability index calculation value became obvious bigger when interval characteristics were considered. The findings indicated that variability was adverse to the slope stability and the correlation， interval characteristics and spatial variability were good for the slope stability. Therefore， with accurate and comprehensive consideration of the soil parameters statistical features， especially with the consideration of space variability and interval characteristics on the basis of the consideration of parameters variability and correlation will be more in line with the engineering practice， and the calculation results tend to be more reliable. It is advantageous to the adopt reliability norm and safety evaluation of slope engineering.

Key words： slope stability reliability， strength parameters， statistical property， interval， spatial variability

邊坡系統具有顯著的模糊性和不確定性，影響邊坡穩定的因素也常常是不確定的，定值法分析沒有考慮這種特性，結果往往偏于危險，比如，安全系數大于設計規范值的邊坡也有可能失穩[1，6]，具有相同安全系數的邊坡，其可靠性卻可能不一樣[2]。可靠度分析是一種基于概率論理論的分析結構物安全程度的分析方法，能夠有效地考慮邊坡系統內實際存在的變異性和相關性，能夠說明“邊坡有多穩”。

20世紀70年代以來，在各類工程結構的設計中全面引入了概率極限狀態的設計原則和方法，但由于對邊坡土質強度參數統計特性的考慮不全面，使得可靠性指標計算值偏小，失效概率計算值偏大，往往超過10%，設計中難以采納，即使對已證明是安全穩定的邊坡工程也是如此[3]。這種理論與實際相背離的情況使得目前國內外對于地基可靠度規范的制定和執行進展緩慢[3]，因此，系統地分析土質強度參數統計特性對邊坡穩定可靠性的影響具有重要的理論意義和實踐價值。

不少學者對此進行了研究，取得了一些研究成果。陳鵬等[4]基于因素敏感性采用Rosenblueth法進行邊坡穩定可靠性分析，發現將敏感性較大的因素作為隨機變量時所得的破壞概率較大，且其變異性也會對結果產生重大影響。李宗坤等[5]探討了參數c、φ的相關性和變異性對土石壩邊坡穩定可靠性的影響，得出參數負相關性對壩坡穩定有利而參數變異性對壩坡穩定不利的結論。李萍等[6]基于點估計法進行了黃土邊坡可靠性分析，發現強度參數的變異性是影響邊坡失效概率大小的主要因素。蔣水華等[7]采用Sarma法和隨機響應面法進行邊坡穩定可靠性分析，發現潛在滑體的可靠指標隨著c、φ間負相關性增加而增大。譚曉慧等[8]采用有限元法和Mohr-Coulomb屈服準則分析了參數離散性和相關性等對邊坡穩定性的影響，發現c、φ對可靠性指標的相對影響比其他參數的影響要大得多，并且隨著c、φ間負相關系數的增加，其對可靠性指標的相對影響也相應增加。朱紅霞等[3]基于隨機場理論和方差折減函數進行了天津港口邊坡可靠度分析，發現φ的變異性對可靠度指標的影響較c更為敏感。鄭敏洲等[9]采用Excel規劃求解工具分析了花崗巖殘積土邊坡穩定可靠性，結果表明邊坡穩定可靠度不僅與土體參數的均值有關，而且與其相關性和變異性有關，且變量的相關性對可靠度影響較大。陳群等[10]分析了土石壩壩料強度參數的均值不定性、相關性和變異性對壩坡失穩概率的影響，表明失穩概率隨土質強度參數均值的增大而減小，隨變異系數的增大而增大，并且兩個強度參數的相關性對壩坡的失穩概率也有影響顯著。

不難看出，現有的研究著重于強度參數的變異性和相關性對穩定可靠性的影響上，而對強度參數的其他統計特性，如參數的區間特性和空間變化性則缺乏考慮。以某高速公路全風化花崗巖土質高邊坡為工程背景，在分析土質參數統計特性的基礎上，利用極限平衡理論和蒙特卡羅模擬法，系統地分析了土質強度參數的均值不定性、變異性、相關性、區間特性和空間變化性等統計特性對邊坡穩定可靠性的影響。

1 理論分析

1.1 蒙特卡羅法基本原理

可靠度計算方法主要有一次二階矩法、Rosenbluth點估計法和蒙特卡羅法等[11-14]。其中蒙特卡羅模擬法是以數理統計原理為基礎的一種抽樣計算理論，該法回避了可靠度分析中數學上的困難，且不受分析條件的限制，隨機變量可以具有任意分布的特點，只要模擬的次數足夠的多，計算精度就能相對的高，因此，蒙特卡羅模擬法在許多領域均具有廣泛的應用。蒙特卡羅法假定隨機變量的概率分布以及變量之間的相關結構都是已知的，生成足夠多組服從給定概率分布規律的隨機變量數據。然后對每組隨機變量分別運行數值模型，再用數理統計的方法評價這些隨機變量所對應的模擬結果即可獲得模擬結果的可靠度。其用于邊坡穩定可靠度分析的基本原理如下：

2 算例分析

2.1 工程概況

計算實例邊坡位于大慶至廣州高速公路江西里仁至龍南段K53+190左坡，該地區地質情況詳見文獻[15]。邊坡共分為5階，前4階坡高為8 m，第5階坡高7.5 m，邊坡總高39.5 m。前3階坡度為1∶1，后2階坡度為1∶1.25。采用有限元軟件Geostudio的Slope/W和Seep/W模塊進行分析，并考慮了地下水位的影響，以年均降雨量為0.3 mm/h的穩態滲流計算結果作為地下水分布的初始狀態。即先用Seep/W模塊進行穩態滲流分析，再將其結果導入Slope/W模塊進行邊坡穩定可靠度分析。有限元模型見圖2，初始壓力水頭分布見圖3，模型共計5 000個單元和5 034個節點。

2.2 全風化花崗巖強度參數統計特性

由于地基中不同地點的巖土在成土環境、礦物成分、顆粒級配組成和應力歷史等方面存在差異，造成地基土中點與點之間的土性會存在不同程度的差異，這種差異在空間上一般表現為漸變性。同時，在土質強度參數測量中，由于試驗設備和操作人員的不同，加上存在人為擾動和試驗誤差，土性參數測試條件與實際情況有出入等因素會導致測量參數的離散性。試驗表明，土性指標與其所處的空間位置密切相關，即土層中任意兩點的同一參數指標（如相鄰土體φ與φ、c與c之間）存在自相關性，土層中同一點的不同參數指標（如相同土體的c和φ之間）存在互相關性[16]，且其相關系數的大小在-0.72～0.35之間[17]，這種相關性會隨著地基中兩點間距離的增大而逐漸減弱，當超過一定的距離后便可認為參數的相關性為0，從而使得土性指標表現出空間變化性。

背景工程邊坡巖土體為全風化花崗巖，簡化等效為各向同性連續介質。根據現場鉆孔資料和室內實驗統計資料，邊坡土體的飽和滲透系數為2.1×10-4 cm/s，飽和含水率為37.5%，殘余含水率為11.9%，重度為19.2 kN/m3，黏聚力c和內摩擦角φ視為符合正態分布隨機變量，其統計特性見表1。

3 計算結果及分析

3.1 空間變化性的影響

參數相關性會隨著地基中兩點間距離的增大而逐漸減弱，當超過一定的距離后便可認為參數的相關性為0，即此時參數各自獨立，因此，當邊坡的潛在滑動面長度顯著且材料單一時，空間變化性就顯得尤為重要。圖4為背景工程邊坡的典型滑動面，在進行蒙特卡羅抽樣試驗時，有3種處理方法：1）對整個滑面的每一層土的參數只取樣一次；2）沿著滑面對每個土條的參數取樣；3）沿著滑面，每間隔一定抽樣距離取樣土的參數。方法1）實際上沒有考慮土質參數的空間變化性，方法2）取樣間距過小，人為地將空間變化性考慮過大，因此，對土只取樣一次和對每個土條都取樣一次都不太符合實際，間隔一定的距離取樣的方法3）比較合適，但方法1）、2）可以用來計算失效概率的上、下限范圍。

為了研究土性參數空間變化性對邊坡穩定可靠性的影響，不考慮表1參數取值的區間特性，取黏聚力c和內摩擦角φ為符合正態分布的隨機變量，計算不同抽樣距離下邊坡穩定可靠性指標，得到可靠性指標隨抽樣距離的變化曲線（圖5），從圖中可見，隨著抽樣距離的增加，可靠性指標總體趨勢是下降的，一開始急劇下降，隨后下降趨勢變緩，最后趨于一個定值。另外，可靠性指標的變化范圍較大，為區間[1.775，9.078]，由此可見，土性參數空間變化性對邊坡穩定可靠性具有較大的影響，考慮空間變化性有利于邊坡的穩定，正確認識和模擬空間變化性是正確采用可靠度計算邊坡穩定性的前提。

根據背景工程的鉆孔資料和室內實驗統計資料，本邊坡的抽樣距離取8 m（即下面的分析均按8 m計算），計算的邊坡穩定系數為1.061，失穩概率為0，可靠性指標為5.295，邊坡處于穩定狀態，但由于計算穩定系數小于1.2，故設計中采用了一定的工程防護措施和防排水措施。

3.2 區間性的影響

從表1可知，試驗得來的土質參數具有最小值和最大值，即在一個區間范圍內，邊坡可靠度研究應可慮這個特點才更加符合工程實際。以表1中黏聚力c為例，c的概率分布密度曲線在不考慮區間和考慮區間兩種情況下是有很大差異的，分別見圖6和圖7。為了研究土質參數區間性對邊坡穩定可靠性的影響，在后面的分析中，均分考慮區間性和不考慮區間性兩種情況進行了對比計算，計算結果及分析見下文。

3.3 均值的影響

保持黏聚力c和內摩擦角φ的變異系數不變，分別改變各自的均值μ進行邊坡穩定性概率分析。在表1參數基準均值的基礎上，分別減小1倍、2倍和3倍標準差，增大1倍、2倍標準差，共獲得6組不同的均值數據，并分考慮區間和不考慮區間兩種情況分別計算。計算工況和結果如表2所示，邊坡的失效概率和可靠性指標隨參數均值的變化曲線見圖8和圖9。

由圖可知，當c和φ的變異系數不變時，邊坡的失穩概率隨c或φ的均值的增大而減小，可靠性指標隨c或φ的均值的增大而增大。其中，φ的均值對邊坡穩定性指標的影響要比c的均值的影響大，即φ的邊坡穩定敏感性要比c大。從參數區間性來看，區間性考慮與否對邊坡的失穩概率的影響不大，但對可靠性指標的影響較大，考慮區間性后，失效概率和可靠性指標隨參數均值的變化曲線變得更陡，即區間性使得均值對邊坡穩定性的影響更加靈敏。

從表2中計算結果可知，c和φ的均值對概率分析所得的安全系數的均值影響較大，安全系數的均值隨c和φ的均值的增大而增大。

3.4 變異性的影響

保持參數c和φ的均值不變，分別改變各自的變異系數進行邊坡穩定可靠性概率分析，并分考慮區間和不考慮區間兩種情況分別計算，計算工況及結果見表3。邊坡的失穩概率和可靠性指標隨c和φ的變異系數的變化曲線見圖10和圖11。由圖可知，在c和φ的均值不變的情況下，當不考慮區間時，邊坡的失穩概率隨c和φ的變異系數的增大而增大，其中φ的變異系數對邊坡失穩概率的影響比c的變異系數的影響大，可靠性指標總體趨勢是隨c和φ的變異系數的增大而減小，可見，參數變異性對邊坡穩定不利；當考慮區間時，失穩概率計算值接近于0，可靠性指標規律性不強，但均明顯大于不考慮區間時的可靠性指標值，即考慮區間后對邊坡穩定有利。

3.5 相關性的影響

在參數c和φ均值和變異系數都不變的情況下，改變c和φ的相關系數進行邊坡穩定可靠性概率分析，并分考慮區間和不考慮區間兩種情況分別計算，計算工況及結果見表4，邊坡的失穩概率和可靠性指標隨c和φ的相關系數的變化曲線見圖12、13和圖14。

由于c和φ互相關的相關系數范圍在-0.72～0.35之間[17]，由圖可知，在各種變異系數情況下，當c和φ的相關系數為0時，邊坡的失穩概率最大而可靠性指標最小，隨著相關系數絕對值的增大，失穩概率逐漸變小而可靠性指標逐漸變大，且相關系數為負數時變化更加劇烈，可見，參數相關性對邊坡穩定有利。從參數區間性來看，考慮區間時失穩概率計算值接近為0，而可靠性指標的變化規律更不考慮區間時類似，但數值上明顯大于不考慮區間時的結果，進一步說明考慮區間后對邊坡穩定有利。另外，從不同變異系數對應曲線的對比分析可進一步說明變異系數越大，邊坡的失穩概率也越大。

從表4中計算結果可知，參數c和φ的相關系數對邊坡穩定可靠度分析所得的安全系數影響很小。

4 結 論

以某高速公路全風化花崗巖土質高邊坡為工程背景，在分析土質參數統計特性的基礎上，利用極限平衡理論和蒙特卡羅模擬法，系統地分析了土質強度參數的均值不定性、變異性、相關性、區間特性和空間變化性等統計特性對邊坡穩定可靠性的影響。

1）土性參數空間變化性對邊坡穩定可靠性具有較大的影響，隨著抽樣距離的增加，可靠性指標總體趨勢是下降的，一開始急劇下降，隨后下降趨勢變緩，最后趨于一個定值，正確認識和模擬空間變化性是正確采用可靠度計算邊坡穩定性的前提。然而，目前抽樣距離的計算尚無定論，建議在綜合分析邊坡潛在滑坡面的長度及所通過的土層的數目等因素確定。

2）強度參數c和φ的均值不定性、變異性和相關性對邊坡穩定可靠性的均具有一定的影響，其中，φ的影響更加顯著。可靠性指標隨c或φ均值的增大而增加；失穩概率隨c和φ的變異系數的增大而增大，可靠性指標隨c和φ的變異系數的增大而減小，參數變異性對邊坡穩定不利；當c和φ的相關系數為0時，邊坡的失穩概率最大而可靠性指標最小，隨著相關系數絕對值的增大，失穩概率逐漸變小而可靠性指標逐漸變大，且相關系數為負數時變化更劇烈，參數相關性對邊坡穩定有利。

3）考慮參數區間性后，各組曲線的變化趨勢大致不變，但是計算的失效概率明顯小于不考慮區間時的結果，而可靠性指標明顯大于不考慮區間時的結果，即考慮區間性的結果要趨于安全，對邊坡穩定有利。

4）安全系數受強度參數均值的影響較大，而對參數的空間變化性、區間性、變異性和相關性均不敏感，由此可見，作為定值法的安全系數不能反應邊坡穩定可靠性的情況。

5）準確而全面地考慮土質參數的統計特性，尤其是在參數變異性和相關性的基礎上加以考慮土性參數空間變化性和區間性會更加符合工程實際，計算結果趨于安全，將有利于地基可靠度規范的制定和推廣運用，有利于更科學合理地評價邊坡工程的安全性。

參考文獻：

[1]姚愛軍， 薛延河. 復雜邊坡穩定性評價方法與工程實踐[M]. 北京： 科學出版社， 2008.

[2]Park H， Westb T R. Development of a probabilistic approach for rock wedge failure [J]. Engineering Geology， 2001， 59（3-4）： 233-251.

[3]朱紅霞，閆澍旺.天津港地基土隨機場特性及可靠度分析[J].巖石力學與工程學報，2007（Sup2）：3959-3965.

Zhu H X， Yan S W. Random field model and reliability analysis of foundation soil in tianjin port[J]. Chinese Journal of Rock Mechanics and Engineering， 2007（Sup2）： 3959-3965.（in Chinese）

[4]陳鵬，徐博侯.基于因素敏感性的邊坡穩定可靠度分析[J].中國公路學報，2012（4）：42-48.

Chen P， Xu B H. Reliability analysis of slope stability based on factors sensitivity[J]. China Journal of Highway and Transport， 2012（4）： 42-48.（in Chinese）

[5]李宗坤，楊艷榮.土體強度參數統計特性對土石壩邊坡可靠性的影響分析[J].中國農村水利水電，2011（8）：146-148.

Li Z K， Yang Y R. The effect of strength parameters on reliability index of soil slope stability[J]. China Rural Water and Hydropower， 2011（8）： 146-148.（in Chinese）

[6]李萍，王秉綱，李同錄，等.陜西地區黃土路塹高邊坡可靠度研究[J].中國公路學報，2009（6）：18-25.

Li P， Wang B G， Li T L， et al. Study of reliability for loess cutting high slopes in Shaanxi province[J]. China Journal of Highway and Transport， 2009（6）： 18-25.（in Chinese）

[7]蔣水華，李典慶.基于隨機響應面法和Sarma法的邊坡可靠度分析[J].鐵道工程學報，2011（7）：21-27.

Jiang S H， Li D Q. Analysis of reliability of slope stability with stochastic response surface and sarma method method[J]. Journal of Railway Engineering Society， 2011（7）： 21-27.（in Chinese）

[8]譚曉慧，王建國，劉新榮，等.邊坡穩定的有限元可靠度計算及敏感性分析[J].巖石力學與工程學報，2007（1）：115-122.

Tan X H， Wang J G， Liu X R， et al. Finite element reliability computation and sensitivity analysis of slope stability[J]. Chinese Journal of Rock Mechanics and Engineering， 2007（1）： 115-122.（in Chinese）

[9]鄭敏洲，簡文彬，吳茂明.花崗巖殘積土邊坡穩定性可靠度分析[J].巖石力學與工程學報，2005（Sup2）：5337-5340.

Zheng M Z， Jian W B， Wu M M. Reliability analysis of stability of granite residual soil slope[J]. Chinese Journal of Rock Mechanics and Engineering， 2005（Sup2）： 5337-5340.（in Chinese）

[10]陳群，唐岷，朱分清.強度參數的不確定性對土石壩壩坡失穩概率的影響[J].巖土工程學報，2008（11）：1594-1599.

Chen Q， Tang M， Zhu F Q. Influence of uncertainty of strength parameters on instability probability of Embankment dam slopes[J]. Chinese Journal of Geotechnical Engineering， 2008（11）： 1594-1599.（in Chinese）

[11]趙明華， 曾昭宇， 劉曉明， 等. 考慮軸向橫向荷載共同作用的基樁可靠度[J]. 建筑科學與工程學報， 2005， 22（2）： 57-60.

Zhao M H， Zeng Z Y， Liu X M， et al. Reliability of pile foundation under vertical and lateral loads [J]. Journal of Architecture and Civil Engineering， 2005， 22（2）： 57-60.（in Chinese）

[12]李志純， 朱道立. 隨機動態交通網絡可靠度分析與評價[J]. 交通運輸工程學報， 2008， 8（1）： 106-112.

Li Z C， Zhu DL. Analysis and evaluation of reliability for stochastic and dynamic transport network [J]. Journal of Traffic and Transportation Engineering， 2008， 8（1）： 106-112.（in Chinese）

[13]闕云， 凌建明， 曾四平. 加筋土路堤內部穩定的可靠指標分析[J]. 交通運輸工程學報， 2006， 6（3）： 37-41.

Que Y， Ling J M， Zeng S P. Reliability index analysis of internal stability of reinforced embankment [J]. Journal of Traffic and Transportation Engineering， 2006， 6（3）： 37-41.（in Chinese）

[14]桂勇，鄧通發，羅嗣海.基于蒙特卡羅邊坡穩定二元體系的建立與應用[J].巖土力學，2014（7）：1979-1986.

GuiI Y， Deng T F， Luo S H. Establishment of slope stability dual index system based on Monte Carlo simulation and its application[J]. Rock and Soil Mechanics， 2014（7）： 1979-1986.（in Chinese）

[15]桂勇，鄧通發，羅嗣海，等.考慮力學水文雙重效應的生態邊坡長期穩定性分析[J].公路，2014（4）：209-216.

Gui Y， Deng T F， Luo S H， et al. Study of long-term stability of ecological slope with double effect of mechanics and hydrology[J]. Highway， 2014（4）： 209-216.（in Chinese）

[16]Krahn J. Stability modeling with SLOPE/W 2007 version， fourth edition： An engineering methodology[R]. Alberta： Geo-slope International Ltd.， 2007.

[17]李亮， 張丙強. c、φ相關性對邊坡整體穩定性的影響[J]. 鐵道科學與工程學報， 2004， 1（1）： 62-68.

Li L， Zhang B Q. Influence of correlations of c and φ on slope stability [J]. Journal of Railway Science and Engineering， 2004， 1（1）： 62-68.（in Chinese）

（編輯 胡 玲）