樊長剛,盛焰正 ,陳雙慶

(1.中國市政工程西南設(shè)計研究總院有限公司，四川 成都 610081;2.中交第二公路勘察設(shè)計研究院有限公司,湖北 武漢 430000；3.中南大學(xué) 土木工程學(xué)院,湖南 長沙 410018)

1 概述

紅砂巖通常指侏羅紀到新近紀的紅色、紅褐色陸相沉積巖系，約占全國陸地總面積的8.61%[1]，主要呈粒狀碎屑結(jié)構(gòu)和泥狀膠結(jié)結(jié)構(gòu)2種形式，在湖南省俗稱“紅砂巖”。湖南省湘西地區(qū)是紅砂巖分布較為集中的地區(qū)，由于紅砂巖往往分布在亞熱帶季風(fēng)氣候區(qū) ，熱量豐富，春夏降水比較集中，雨量充沛，紅砂巖容易在水、氧氣、風(fēng)等因素下風(fēng)化，極大影響了紅砂巖的抗剪強度與其邊坡的穩(wěn)定性[2]。學(xué)者們已經(jīng)開展通過一系列試驗研究了紅砂巖抗剪強度隨環(huán)境因素的變化規(guī)律，例如陳文杰[3]考慮了紅砂巖崩解性，獲取了抗剪強度參數(shù)隨含水率的變化規(guī)律。徐志華[4]等開展了干濕循環(huán)作用下紅砂巖強度劣化特性試驗研究，定量分析了干濕循環(huán)作用對紅砂巖強度力學(xué)參數(shù)的影響，并提出了基于黏聚力的損傷變量。劉新榮[5]等指出當經(jīng)歷20次干濕循環(huán)后，紅砂巖的c和φ分別降低58%和34%。上述試驗結(jié)果為考慮環(huán)境影響的邊坡穩(wěn)定性分析和支護設(shè)計提供了試驗資料，但紅砂巖化學(xué)成分和結(jié)構(gòu)的復(fù)雜性導(dǎo)致紅砂巖的性質(zhì)具有很強的空間差異性，即使在同一地區(qū)，其抗剪強度范圍區(qū)間也分布很廣，這對需要明確巖土材料抗剪強度參數(shù)的傳統(tǒng)邊坡穩(wěn)定性分析帶來很大的挑戰(zhàn)[6]。實際工程中，勘察設(shè)計人員往往在試驗數(shù)據(jù)甚至缺乏試驗數(shù)據(jù)的基礎(chǔ)上，根據(jù)經(jīng)驗確定某一定值來進行穩(wěn)定性計算，結(jié)果帶有較大的主觀性。

可靠性理論常用于解決某一系統(tǒng)內(nèi)的不確定性問題，通過計算邊坡的可靠性指標和破壞概率，能夠反映邊坡的各種不確定影響因素，避免了安全系數(shù)使用過程中的絕對化。鑒于此，針對湖南省某公路的風(fēng)化紅砂巖樣品進行了抗剪強度參數(shù)統(tǒng)計分析，建立了正態(tài)分布模型，并開展了利用蒙特卡洛法和羅森布魯斯法的邊坡可靠性分析，研究了支護前后可靠性指標的變化規(guī)律，并將可靠性方法分析結(jié)果與傳統(tǒng)定值方法進行了對比，旨在為今后類似工程的可靠性分析提供參考。

2 邊坡可靠性分析理論

2.1 可靠性狀態(tài)方程與指標

常規(guī)的邊坡可靠性理論可以參照文獻[7]。類似結(jié)構(gòu)可靠性，采用傳統(tǒng)的邊坡安全系數(shù)進行狀態(tài)判斷，即安全系數(shù)表示功能函數(shù)：

(1)

與上式對應(yīng)的失效概率為：

(2)

因此，安全系數(shù)功能函數(shù)對應(yīng)的可靠性指標β為：

(3)

式中：μz為功能函數(shù)Z的均值；σz表示Z的標準差。邊坡可靠性指標與破壞概率的關(guān)系如表1所示。

表1 邊坡可靠性指標與失效概率的關(guān)系Table 1 Relationship between slope reliability index and failure probability性能等級可靠性指標失效概率性能等級可靠性指標失效概率高5.02.871×10-7差2.00.022 75好4.02.169×10-5糟糕1.50.066 81中等以上3.00.001 35危險1.00.158 66中等以下2.50.006 21

2.2 可靠性分析方法

關(guān)于邊坡可靠性研究已取得了許多成果，可靠性指標的計算方法包括蒙特卡洛法[8]、一次二階矩法[9]、Rosenblueth法[10]等。

利用蒙特卡洛法(MCS)方法計算邊坡失效概率的估計值可表示為[11]：

(4)

式中：xi(i=1，2，…，N)為模擬的隨機樣本。

蒙特卡洛法[12]又稱為統(tǒng)計實驗法或隨機模擬法，它是一種數(shù)值計算方法，其中蒙特卡洛法需要預(yù)知各基本變量的分布形式和參數(shù)特征，采用蒙特卡洛法計算須隨機生成N組隨機數(shù)，而每組隨機數(shù)的生成均是以隨機變量的分布形式和參數(shù)特征為基礎(chǔ)的，它只是一種數(shù)值計算方法，沒有固定的求解公式，只能求得邊坡在給定條件和參數(shù)下的解，但這恰恰回避了邊坡可靠度分析中數(shù)學(xué)模型的建立困難。

羅森布魯斯法又稱統(tǒng)計矩法，相對蒙特卡洛法，可以大大減少可靠性分析中的計算量，具體理論詳見文獻[13]和[14]。對于邊坡問題，假設(shè)模糊參數(shù)為c和φ值，則功能函數(shù)Z=g(c，φ)。對于每個參數(shù)取2個計算點，則有：

Z++=g[(μc+σc)，(μφ+σφ)]

Z+-=g[(μc+σc)，(μφ-σφ)]

Z-+=g[(μc-σc)，(μφ+σφ)]

Z--=g[(μc-σc)，(μφ-σφ)]

(5)

對于普通的可靠性分析，羅森布魯斯法的一階矩M1即均值μz，二階矩M2即方差：

(6)

式中：下標j=1、2、3、4分別代表j=++、+-、-+、--。

P++=P--=0.25(1+ρXY)，P+-=P-+=0.25(1-ρXY)，ρXY為相關(guān)系數(shù)。

當引入模糊理論時，利用模糊判決中的普通加權(quán)平均法，則m個λ水平下考慮模糊隨機性的狀態(tài)函數(shù)的一階原點矩(均值)為：

(7)

二階矩(方差)為：

(8)

3 實例分析

3.1 工程概況

湖南省湘西某高速公路屬低山丘陵地貌，地形起伏較大，山體自然邊坡約為10°～40°。邊坡大部分呈上下2層的“二元結(jié)構(gòu)”，上部為泥質(zhì)粉砂巖，巖質(zhì)軟，厚度15～20 m，結(jié)構(gòu)大部分已經(jīng)破壞，風(fēng)化裂隙很發(fā)育，巖體很破碎，根據(jù)《公路工程地質(zhì)勘察規(guī)范》(JTG C20-2011)[15]，屬于強風(fēng)化巖；下部為粉砂巖，結(jié)構(gòu)部分破壞，風(fēng)化裂隙發(fā)育，巖體被切割成巖塊，屬于中風(fēng)化巖。二者皆為該高速公路項目的控制性地層，可劃于紅砂巖的范疇，下文簡稱“強風(fēng)化紅砂巖”和“中風(fēng)化紅砂巖”,見圖1。

圖1 代表性巖樣

地勘階段在邊坡沿線周邊2 km范圍內(nèi)進行了大量鉆孔取樣，該范圍內(nèi)無大型地質(zhì)構(gòu)造，地下水位較深，巖樣含水率主要受天然降雨入滲所控制，而受地下水影響可能性不大。參數(shù)獲取與統(tǒng)計按照《公路土工試驗規(guī)程》(JTG 3430—2020)[16]、 《公路工程結(jié)構(gòu)可靠性設(shè)計統(tǒng)一標準》(JTG 2120—2020)[17]開展。圖2為不同位置紅砂巖抗剪強度參數(shù)的變化情況，可以看出，紅砂巖的抗剪強度具有明顯的不確定性，這是由礦物成分、含水率等因素的不均勻所造成的，一般而言，石英含量越高，抗剪強度參數(shù)越高，而含水率越高，抗剪強度參數(shù)則越低。盡管各個樣品的抗剪強度參數(shù)有較大差異，但整體上均值比較穩(wěn)定，沒有發(fā)生漂移，具有正態(tài)分布的特征，風(fēng)化程度越高，不確定性更明顯。

對上述結(jié)果進行整理，得到抗剪強度參數(shù)正態(tài)分布曲線如圖3所示。從圖3可以看出，巖石抗剪強度參數(shù)c、φ值均近似服從正態(tài)分布，強風(fēng)化紅砂巖參數(shù)的不確定性更大。

(a) 黏聚力

(b) 內(nèi)摩擦角

(a) c的參數(shù)分布

(b) φ的參數(shù)分布

用于可靠性分析的案例邊坡為四級坡，高度約29 m。邊坡內(nèi)匯水面較小，地表水不發(fā)育，地下水貧乏，水文地質(zhì)條件簡單，對邊坡穩(wěn)定性影響小。一級、二坡比為1∶0.75，三、四級坡比為1∶1。一至三級坡采用錨桿防護，錨桿長9 m；四級坡采用植草防護。在支護前，坡體發(fā)生了2次小規(guī)模失穩(wěn)破壞。圖4為該邊坡支護前后的照片。

(a) 支護前

(b) 支護后

3.2 蒙特卡洛法計算結(jié)果分析

在蒙特卡洛法基礎(chǔ)上，利用簡化Bishop條分搜索滑面并計算功能函數(shù)值(即安全系數(shù))，共進行了800次抽樣計算。天然未支護狀態(tài)和采用錨桿支護狀態(tài)的計算結(jié)果如圖5所示。

從圖5可以看出，支護前后最危險滑面均從坡頂開始發(fā)育，在第一級坡的坡面中上部出露，是否支護對滑面的形態(tài)影響不大。圖6為2種工況下安全系數(shù)的分布概率直方圖，可知邊坡的安全系數(shù)基本上符合正態(tài)分布。對圖中安全系數(shù)小于1的計算次數(shù)進行統(tǒng)計，邊坡在未支護情況下安全系數(shù)小于

(a) 未支護

(b) 已支護

1的次數(shù)為180次，即邊坡失效破壞概率為22.6%，說明該邊坡在不采取任何支護措施的情況下，發(fā)生破壞的概率比較大，這與強風(fēng)化紅砂巖抗剪強度較低有關(guān)；而邊坡采取錨桿支護措施后，安全系數(shù)小于1的次數(shù)為38次，即邊坡的失效破壞概率為4.8% ，說明錨桿加固后邊坡穩(wěn)定性顯著提高。

(a) 未支護

(b) 已支護

參照羅文強[18]等的報道，從失效概率角度出發(fā)，可以將邊坡穩(wěn)定性分為5個等級，見表2。

表2 基于失效概率的邊坡穩(wěn)定性等級Table 2 Slope stability grade based on failure probability穩(wěn)定性描述失效概率/%穩(wěn)定等級不穩(wěn)定>90 1穩(wěn)定性差60～902欠穩(wěn)定20～603基本穩(wěn)定5～204穩(wěn)定<5 5

根據(jù)之前的分析，支護前邊坡失效概率Pf=22. 6%，處于欠穩(wěn)定狀態(tài)，這與現(xiàn)場出現(xiàn)失穩(wěn)的情況相符，進行支護后，邊坡失效概率Pf=4.8%，表明支護后邊坡處于穩(wěn)定狀態(tài)，體現(xiàn)了進行錨桿支護的工程意義。

3.3 羅森布魯斯法計算結(jié)果分析

考慮到計算量，僅采用羅森布魯斯法計算支護后邊坡的可靠性。表3為采用羅森布魯斯法時的計算參數(shù)。下標1代表強風(fēng)化紅砂巖，下標2代表中風(fēng)化紅砂巖，例如c1-代表強風(fēng)化紅砂巖c值在小于統(tǒng)計期望值一側(cè)模糊估計值。

表3 計算參數(shù)Table 3 Parameters for Calculation隸屬度不確定性參數(shù)c1-/kPac1+/kPaφ1-/(°)φ1+/(°)c2-/kPac2+/kPaφ2-/(°)φ2+/(°)0.1011.230.88.2225.58385224.2635.740.2011.930.18.8424.9638.551.524.6735.330.3012.629.49.4624.34395125.0834.920.4013.328.710.0823.7239.550.525.4934.510.50142810.723.1405025.934.10.6014.727.311.3222.4840.549.526.3133.690.7015.426.611.9421.86414926.7233.280.8016.125.912.5621.2441.548.527.1332.870.9016.825.213.1820.62424827.5432.46

不同隸屬度下安全系數(shù)和可靠性指標的變化規(guī)律如圖7所示?？梢钥闯?，可靠性指標與安全系數(shù)同時隨著參數(shù)隸屬度的增加而提高。根據(jù)式(3)、式(7)和式(8)，邊坡的可靠性指標β=1.70，失效概率為Pf=6.3%。由于羅森布魯斯法考慮了巖土參數(shù)的模糊隨機性，因此計算出的失效概率是大于蒙特卡洛法的。

圖7 不同隸屬度下的安全系數(shù)均值與可靠性指標Figure 7 Safety factor mean and reliability index under different memberships

3.4 與傳統(tǒng)方法對比

為了進一步對比可靠性分析法與傳統(tǒng)定值安全系數(shù)法的區(qū)別，選擇了3名工作經(jīng)驗在5 a以上的公路勘察設(shè)計人員，其根據(jù)工作經(jīng)驗分別給出了強風(fēng)化紅砂巖和中風(fēng)化紅砂巖抗剪強度參數(shù)經(jīng)驗值，并利用簡化Bishop法進行經(jīng)驗值參數(shù)下的穩(wěn)定性計算，結(jié)果見表4。

表4 計算參數(shù)Table 4 Calculation Parameters 工程人員編號強風(fēng)化紅砂巖c/kPa中風(fēng)化紅砂巖c/kPa強風(fēng)化紅砂巖φ/(°)中風(fēng)化紅砂巖φ/(°)安全系數(shù)(支護前)安全系數(shù)(支護后)1號126015431.281.512號155018401.231.453號105020351.031.18

根據(jù)《公路路基設(shè)計規(guī)范》JTG D30—2015[19]，對于高速公路正常工況，安全系數(shù)為1.20～1.30視為穩(wěn)定。受到礦物成分、含水率等因素影響，現(xiàn)場紅砂巖抗剪強度具有明顯的不確定性，采用傳統(tǒng)確定性安全系數(shù)法時，由于個人主觀經(jīng)驗原因，1號和2號工程人員估計的計算參數(shù)使得未支護前邊坡處于穩(wěn)定狀態(tài)，與現(xiàn)場情況有所偏差，導(dǎo)致邊坡設(shè)計支護手段不足而出現(xiàn)失穩(wěn)破壞危險。而3號工程人員的計算結(jié)果偏保守，即支護之后安全系數(shù)仍然未達穩(wěn)定要求，可能造成支護措施的過度加強而引起新的工程成本開支。 同時可以看出，采用傳統(tǒng)安全系數(shù)法計算時，支護前邊坡的安全系數(shù)都大于1，但實際邊坡出現(xiàn)了失穩(wěn)破壞的現(xiàn)象，這從傳統(tǒng)安全系數(shù)法理論上是難以解釋的。而可靠性理論認為任何邊坡的失穩(wěn)破壞都是概率問題，只是支護前失效概率高，支護后失效概率低，從而為參數(shù)不確定條件下的邊坡穩(wěn)定性分析提供了更科學(xué)的理論支撐。

4 結(jié)論

本文以湘西某紅砂巖區(qū)高速公路邊坡建設(shè)為例，在考慮巖土參數(shù)空間不確定性的前提下，利用蒙特卡洛法、羅森布魯斯法進行了邊坡可靠性分析，得到如下結(jié)論：

a.根據(jù)勘察統(tǒng)計數(shù)據(jù)，紅砂巖的抗剪強度具有明顯的不確定性，這是由礦物成分、含水率等因素的不均勻所造成的，一般而言，石英含量越高，抗剪強度參數(shù)越高，而含水率越高，抗剪強度參數(shù)則越低。

b.根據(jù)蒙特卡洛法計算結(jié)果，支護前紅砂巖邊坡失效概率Pf=22.6%，處于欠穩(wěn)定狀態(tài)，這與現(xiàn)場出現(xiàn)失穩(wěn)破壞情況相符；進行支護后，邊坡失效概率Pf=4.6%，表明支護后邊坡處于穩(wěn)定狀態(tài)，失穩(wěn)概率很小。

c.根據(jù)羅森布魯斯法計算結(jié)果，支護后紅砂巖邊坡的可靠性指標β=1.70，失效概率為Pf=6.3%，相對蒙特卡洛法，采用羅森布魯斯法計算得出失效概率更大，偏于保守。

d.可靠性分析理論為參數(shù)不確定條件下的邊坡穩(wěn)定性分析提供了更科學(xué)的理論支撐，可避免工程人員因為個人經(jīng)驗對巖土參數(shù)做出過高或過低的估計，從而導(dǎo)致設(shè)計支護強度不足。尤其在大數(shù)據(jù)時代，可以利用超高次數(shù)計算來獲取高精度的可靠度計算結(jié)果，具有十分廣闊的應(yīng)用前景。