西藏某水庫工程土石壩穩定性分析

姚 莉 李冰妮 白興鋒

(西藏自治區水利電力規劃勘測設計研究院，西藏 拉薩 850000)

?

西藏某水庫工程土石壩穩定性分析

姚 莉 李冰妮 白興鋒

(西藏自治區水利電力規劃勘測設計研究院，西藏 拉薩 850000)

以西藏某水庫工程碾壓式瀝青混凝土心墻壩為研究對象，利用STAB軟件，分析了該工程土石壩穩定性與壩體材料參數的敏感性，結果表明，庫水位驟降情況下壩坡穩定系數較低，其他工況基本滿足規范要求；在對壩體材料力學參數進行90%折減后，壩坡穩定系數相應降低，有多重工況基本不滿足規范要求，建議進一步合理設計壩坡坡角，并嚴格控制施工質量。

巖土工程，水庫，壩坡，穩定性

西藏某水庫位于西藏山南地區乃東縣索珠鄉境內雅魯藏布江左岸一級支流旺曲上游丁拉溝，壩址處地理坐標為東經91°52′36″，北緯29°30′25″，距山南地區澤當鎮41.0 km，距拉薩市186 km。攔河壩為碾壓式瀝青混凝土心墻堆石壩，壩頂高程4 179.35 m，防浪墻頂高程4 180.55 m，壩頂寬8 m，最大壩高67.6 m，壩頂長度365.0 m。滲控措施由瀝青混凝土心墻、防滲墻及灌漿帷幕組成。土石壩作為一種堆石材料壩，大量專家學者進行了穩定性方法及穩定性研究[1-4]。張明等針對在水利工程中土石壩邊坡穩定分析和設計相對落后的現狀，以概率統計和結構可靠度理論為基礎，對土石壩邊坡穩定可靠度分析和設計問題，包括土石壩抗剪強度參數粘聚力和內摩擦角的統計方法和統計特性，土石壩邊坡穩定可靠度分析方法，設計表達式中標準值和分項系數的取值問題等，進行了分析研究[5]。李南生等人將非線性統一強度理論應用于土石壩邊坡穩定分析中，考慮壩體各應力分量包括中間主應力對土質材料強度的影響，以區別于傳統土石壩穩定分析中僅考慮剪切和拉伸極限強度的極限平衡法。提出了一個土石壩邊坡穩定優化分析模型，并采用最優化遺傳算法對土石壩邊坡進行斷裂面優化搜索[6]。孫麗麗等人應用有效應力原理，以初始應力方式將滲流場導入到位移應力分析中，基于強度對壩體邊坡整體進行彈塑性分析[7]。

文中利用中國水利水電科學研究院所編制的土質邊坡穩定分析程序STAB對西藏某水利水庫工程土石壩壩坡穩定性進行了分析研究，并對壩體材料參數敏感性進行了分析，研究結果可為該工程壩坡設計提供參考，為類似工程穩定性計算提供一定的借鑒。

1 計算原理及方法

假定滑動面為一圓弧，其半徑為R。任取一土條i，其上作用力有土條自重Wi；作用于土條底面的抗滑剪力Ti、有效法向反力Ni′及孔隙水壓力uili，假定這些力的作用點都在土條底面中點。除此以外，在土條兩側還分別作用有法向力Ei和Ei+1及切向力Xi和Xi+1，Xi+1-Xi=ΔXi。Ei及Xi分別為法向及切向條間作用力，Wi為土條自重，Qi為水平向地震慣性力，Ni，Ti分別為土條底部的總法向力(包括有效法向力及孔隙水壓力)和切向力。根據每一條豎向力的平衡條件有:

Wi+Xi-Xi+1-Tisinαi-Nicosαi=0

(1)

或

Nicosαi=Wi+Xi-Xi+1-Tisinαi

(2)

按照安全系數的定義及摩爾—庫侖準則，Ti可寫成式(3)：

(3)

將式(3)代入式(2)得土條底部總法向力為：

(4)

(5)

在極限平衡時，各土條對圓心的力矩之和應當為0，此時土條間作用力將相互抵消，因此得：

∑Wixi-∑TiR+∑Qiei=0

(6)

將式(3),式(4)代入式(6)，且xi=Rsinαi，最后得到安全系數的公式為：

(7)

其中，Xi及Xi+1是未知的，為使問題得解，畢肖普又假定各土條之間力均略去不計，也就是假定ΔX=0，這樣式(7)可簡化成：

(8)

因為在mai內有Fs因子，所以在求Fs時要進行迭代。在計算時，一般可以先假定Fs=1，求出mai(或假定mai=1)，再求Fs，再用此Fs求出新的mai及Fs，如此反復迭代，直至假定的Fs與算出的Fs非常接近為止。

2 計算斷面及壩體材料參數

2.1 計算斷面

計算斷面為該土石壩最大斷面位置，即樁號0+000 m段。

2.2 壩體材料參數

根據地勘基本資料，壩體各區堆石料計算參數見表1。

表1 各區計算參數

3 計算結果分析

3.1 計算工況

結合規范要求，共計算了4種不同的工況，具體如表2所示。

表2 計算工況

3.2 計算結果

對該水庫工程土石壩采用了常用的簡化畢肖普方法進行了淺層及深層滑坡穩定計算，此外，還考慮了90%的強度折減情況。各種不同工況下，大壩淺層壩坡穩定性計算結果如表3所示，其中工況1和工況3滑裂面如圖1,圖2所示。

表3 各工況上、下游壩坡淺層的抗滑穩定安全系數

工況上游安全系數下游安全系數正常折減正常折減11.4391.2971.2721.14621.3821.244——30.6610.595——41.1461.0321.1291.018

從表3可以看出各工況下上、下游壩坡的抗滑穩定安全系數基本上都大于允許安全系數，但是在水庫驟降條件下，壩坡穩定安全系數較低，但是需要說明的是，計算壩坡淺表層的穩定基本上沒有特別大的意義，當考慮下游護坡等措施后，淺表層的壩坡穩定問題基本上不存在問題。另外，當考慮大壩材料參數的90%折減情況下，壩坡淺層滑動穩定安全系數有所降低，在折減后，除穩定滲流期工況以及施工完建工況上游壩坡滿足規范要求外，其他情況下壩坡穩定均不滿足規范要求。

表4 各工況上、下游壩坡深層的抗滑穩定安全系數

各種不同工況下，大壩深層壩坡穩定性計算結果如表4所示，其中工況1和工況3滑裂面如圖3，圖4所示。

從表4可以看出各工況下上、下游壩坡的抗滑穩定安全系數都大于允許安全系數，但是在水庫驟降條件下，壩坡穩定安全系數較低。另外，當考慮大壩材料參數的90%折減情況下，壩坡淺層滑動穩定安全系數有所降低，在折減后，壩坡穩定均滿足規范要求。

4 結語

利用中國水利水電科學研究院所編制的土質邊坡穩定分析程序STAB對西藏某水利水庫工程土石壩壩坡穩定性進行了分析研究，結果表明，庫水位驟降情況下壩坡穩定系數較低，其他工況基本滿足規范要求。在對壩體材料力學參數進行90%折減后，壩坡穩定系數相應降低，有多重工況基本不滿足規范要求，建議對壩坡坡角進一步合理設計，并嚴格控制施工質量。

[1] 陳祖煜.土質邊坡穩定分析——原理·方法·程序[M].北京:中國水利水電出版社,2003.

[2] 鄭穎人,趙尚毅,張魯渝.用有限元強度折減法進行邊坡穩定分析[J].中國工程科學,2002,4(10):58-61.

[3] 吳海真,顧沖時.有限元強度折減法在土石壩邊坡穩定分析中的應用[J].水電能源科學,2006,24(4):55-63.

[4] 趙劍明.高土石壩地震反應的三維非線性有效應力分析[D].北京:中國水利水電科學研究院,1999.

[5] 張 明,劉金勇,麥家煊.土石壩邊坡穩定可靠度分析與設計[J].水力發電學報,2006,25(2):103-107.

[6] 李南生,唐 博,談風婕,等.基于統一強度理論的土石壩邊坡穩定分析遺傳算法[J].巖土力學,2013,34(1):243-249.

[7] 孫麗麗,李宗利.基于強度折減法的土石壩邊坡穩定分析與探討[J].水資源與水工程學報,2009,20(5):113-116.

Analysis of earth-rock dam slope stability of a reservoir project in Tibet

Yao Li Li Bingni Bai Xingfeng

(TibetHydropowerPlanningandDesignInstitute,Lhasa850000,China)

Using the software of STAB, the earth-rock dam slope of a reservoir project was studied. The dam material parameter sensitivity analysis was executive also. The results show that the stability of slope is lower when the water sharply drops in reservoir, when other conditions basically satisfy the code requirements. When the dam mechanical parameters of the material reduction of 90%, slope stability coefficient corresponding decrease, there are a variety of conditions does not meet the code specifications. Recommended slope angle of slope design is reasonable, and strictly control the quality of construction.

geotechnical engineering, reservoir， dam slope， stability

1009-6825(2016)22-0071-02

2016-05-22

姚 莉(1970- )，女，工程師； 李冰妮(1983- )，女，助理工程師； 白興鋒(1978- )，男，工程師

TV641

A