楊 奎

(浙江省水利水電勘測設(shè)計院 惠州分院，廣東 惠州 516000)

1 概 述

在我國水電開發(fā)的大背景下，許多大型水電站相繼投入建設(shè)。水電站的建設(shè)將給當?shù)亓饔虻沫h(huán)境帶來一定的影響，水電站蓄水后，邊坡兩岸巖體受到水壓力及水力侵蝕的作用，巖體變得疏松，影響水電站邊坡的整體穩(wěn)定性，蓄水過程帶來的庫岸邊坡變形破壞等穩(wěn)定性問題逐漸成為制約水電開發(fā)的工程地質(zhì)問題之一。大型水利工程蓄水后將形成大量的涉水邊坡，導致水電站邊坡的穩(wěn)定性易受到區(qū)域環(huán)境條件的影響，如氣候、降水等。水電站邊坡的穩(wěn)定性是自然因素與人為活動共同作用的產(chǎn)物，具有復雜的系統(tǒng)性，存在較多的不確定性因素。

邊坡安全系數(shù)是巖質(zhì)邊坡穩(wěn)定性評價應用最廣泛的指標。實際工程中，通常采用抗滑力與滑動力的比值來定義安全系數(shù)。雖然形式簡單，但帶有較多的經(jīng)驗成分，并且它是建立在確定性的概念之上。巖質(zhì)邊坡巖體的穩(wěn)定性受巖的巖性、強度、構(gòu)造、地下水位的高低、結(jié)構(gòu)面的產(chǎn)狀和位置、爆破震動等因素影響，導致穩(wěn)定性計算中含有若干不確定性。為克服現(xiàn)有穩(wěn)定性計算的缺點，建立在不確定性概念之上的概率分析法被引入到巖質(zhì)邊坡的穩(wěn)定性評價中，結(jié)合公式計算法與概率分析法的巖質(zhì)邊坡穩(wěn)定性分析方法可靠性更高。

Rosenblueth法又稱概率矩點估計法(以下簡稱R法)，是由Rosenblueth于1975年提出的一種近似的求多個隨機變量矩估計的方法，其基本原理較之改進一次二階矩陣及蒙特卡洛法要相對簡單且應用方便。Rosenblueth法的基本思想：對于狀態(tài)函數(shù)中的各隨機變量，利用各變量的均值與方差表示的某些特定的點值，通過狀態(tài)函數(shù)的響應，求得狀態(tài)函數(shù)的各階矩。通過Rosenblueth法求得各個評價指標，由于Rosenblueth法原理簡單、實用性強，對于巖質(zhì)高陡邊坡的穩(wěn)定性評價可靠性較高。水電站邊坡工程具有較強的復雜性及不確定性，利用可靠度理論，應用蒙特卡洛法和改進Rosenblueth法進行水電站邊坡可靠性分析。

2 可靠度理論

對于水電站巖質(zhì)邊坡穩(wěn)定性分析，根據(jù)水電站邊坡的整體結(jié)構(gòu)特點、失穩(wěn)破壞機理、巖體力學狀態(tài)以及環(huán)境作用的影響，建立安全系數(shù)狀態(tài)函數(shù)，綜合考慮黏聚力、內(nèi)摩擦角、巖石塊體結(jié)構(gòu)等隨機變量，對于本文中的各隨機變量可認為服從正態(tài)分布或者對數(shù)正態(tài)分布，建立狀態(tài)函數(shù)表達式。公式如下：

Z=F(X1，X2，X3，…，Xn)

(1)

式中:X1、X2、X3、…、Xn為水電站巖質(zhì)邊坡的巖體的黏聚力、內(nèi)摩擦角、密度、水電站泄洪蓄水期水位等邊坡力學性質(zhì)等隨機變量，文中的各隨機變量可認為服從正態(tài)分布或者對數(shù)正態(tài)分布。

通常影響水電站邊坡穩(wěn)定性的主要依賴于巖體的抗剪強度，即內(nèi)摩擦角和黏聚力。邊坡穩(wěn)定性方法主要包括極限平衡法、Morgenstern-Price法、不平衡推力法等。文中選取Morgenstern-Price法中邊坡安全系數(shù)計算方程為狀態(tài)函數(shù)：

作用于條塊上的力包括：巖塊自重、滑動面的分布水壓力和表面荷載如地震荷載，處于蓄水狀態(tài)的庫水對部分條塊產(chǎn)生水壓力,單獨條塊主要受條塊的自重Gi、外部荷載Wi、條塊間法向力Ei及Ei-1、與底面的垂直距離分別是zi及zi-1、水壓力Ui、抗剪力Fp、Fs為邊坡穩(wěn)定安全系數(shù)，假定法向作用力與切向作用力之間的關(guān)系為Y=μg(x)X,μ為力平衡系數(shù)及g(x)為巖體條塊間法向作用力與切向作用力的關(guān)系函數(shù)。一般為常值函數(shù)或者半正弦函數(shù)，則巖體條塊間的剪切力為μgi-1Ei-1、μgiEi。Rosenblueth法是一種近似求隨機函數(shù)數(shù)學期望的方法，對于未知分布函數(shù)的隨機變量，只需要在變量分布區(qū)間上取對稱點，一般可用整體均值加減標準差，即式(2)確定狀態(tài)函數(shù)取值點：

X1=μx1+σx1

X2=μx2-σx1

(2)

利用概率分布的知識，對于m個相互獨立的隨機變量，狀態(tài)函數(shù)取值點有2m個，根據(jù)所有的取值點與狀態(tài)函數(shù)對應關(guān)系，求得2m個狀態(tài)函數(shù)，也就求得2m個安全系數(shù)。對于n個隨機變量的一般情況，設(shè)功能函數(shù)可由式(2)表示，其中Xn為對功能函數(shù)產(chǎn)生影響的各隨機變量。由Rosenblueth方法可以求得功能函數(shù)的均值與方差：

(3)

特殊的，對于本文中所確定的隨機變量，通過對邊坡力學性質(zhì)的分析，可以認為隨機變量是不相關(guān)，對于n個相互獨立的隨機變量，每一取值點出現(xiàn)的概率相等。

(4)

根據(jù)中心距與原點矩的估計，可以導出安全系數(shù)概率分布的4階表達式，由此可估計出概率分布的空間形態(tài)和位置。如果狀態(tài)變量存在相關(guān)關(guān)系，則其相互組合的出現(xiàn)概率不一，出現(xiàn)概率大小應取決于狀態(tài)變量之間的相關(guān)系數(shù)，則：

式中：i=1，2,3,...,n；當xi=xi1時，ei=1，當xi=xi2時，ei=-1；φij為變量xi,xj之間的相關(guān)系數(shù)，則狀態(tài)函數(shù)可以寫為：

3 水電站邊坡工程算例

水電站邊坡位于廣東省境內(nèi)，引水口穿越陡峭山體段，水電站上下游水位界限相距320 m，正常蓄水高程為150 m。從河岸邊坡上游往下游170～230 m，為花崗巖露出段，風化侵蝕較為嚴重，岸坡體較為完整無明顯沖溝及裂縫，坡度約40～48°，為高陡巖質(zhì)邊坡；左岸主要巖石為粉質(zhì)砂巖，巖體較為疏松，坡頂上部有小型裂縫，裂縫寬度約為3～4.2 m，坡度為46～52°，是水電站易發(fā)生滑坡的最大危險區(qū)域。將Morgenstern-Price法中的抗剪強度參數(shù)以及巖體重度作為狀態(tài)函數(shù)的隨機變量。水電站邊坡安全系數(shù)采用Morgenstern-Price法和不平衡推力法共同求解，將滑動面所受的力，按滑動面法向與切向分解，根據(jù)力的平衡原理以及力矩平衡方程，得到水電站邊坡安全系數(shù)如下，計算結(jié)果見表1。

(5)

根據(jù)不平衡推力法計算巖體條塊間法向力：

Ei=Ei-1ζi+FsGh-cili-Gvtanφi

(6)

式中：Ei為第i個條塊推力，kN；傳遞系數(shù)ζi=cos(θi-1-θi)-tanφisin(θi-1-θi)；Gh為滑體重量沿滑動面方向分量，kN；Gv為滑體重量垂直滑動面方向分量，kN。

表1 改進Rosenblueth法計算結(jié)果Tab.1 Improved Rosenblueth calculation results

續(xù)表1

對于n個相互獨立的隨機變量，每一取值點出現(xiàn)的概率相等，由于兩個隨機變量相互獨立:

求得狀態(tài)函數(shù)的均值:

狀態(tài)函數(shù)的方差為：

狀態(tài)函數(shù)服從正態(tài)分布，則可靠指標為

邊坡發(fā)生破壞概率為

Pf=1-Φ(β)=1-0.735=26.5%

采用蒙特卡洛法進行對比計算(運行迭代次數(shù)為3 000次)，對比結(jié)果見表2。

表2 水電站邊坡不同計算方法結(jié)果Tab.2 Different calculation methods of the slope of the hydropower station results

根據(jù)表2可得，改進Rosenblueth法與蒙特卡洛法所求得的狀態(tài)函數(shù)均值、方差、可靠指標最大相差約6%，可以驗證改進Rosenblueth法的對水電站邊坡分析的有效性。水電站邊坡發(fā)生滑坡的概率為26.5%。

4 結(jié) 論

在巖質(zhì)邊坡穩(wěn)定性可靠度分析的方法中，改進一次二階矩法、蒙特卡羅法由于在計算過程中需要對狀態(tài)函數(shù)求偏導，但狀態(tài)函數(shù)往往具有復雜的表達形式，蒙特卡洛還需要多次的迭代，往往出現(xiàn)不收斂的情況，大大降低了兩種方法的適用性。通過Rosenblueth方法與蒙特卡洛法求解對比，改進Rosenblueth方法是一種原理簡單且計算結(jié)果較為精確的方法，對于復雜不宜求偏導的狀態(tài)函數(shù)方程顯示較強的適用性。水電站高陡巖質(zhì)邊坡穩(wěn)定性分析是一項復雜的系統(tǒng)工程，水電站邊坡的穩(wěn)定性具有重要的經(jīng)濟效益、社會效益及環(huán)境效益，需要對其進行深入的分析與了解。通過上述分析可知，改進Rosenblueth方法可以作為高陡巖質(zhì)邊坡穩(wěn)定性分析的一種判據(jù)。