深覆蓋層上的某面板堆石壩抗震穩定分析

邢建營，呂小龍，杜全勝

(1.黃河勘測規劃設計研究院有限公司，河南 鄭州 450003；2.水利部水利水電規劃設計總院，北京 100120)

我國已建水庫大壩總數約9.8萬座，其中90%以上為土石壩。在水能資源豐富的西南和西北地區，地質條件較理想的壩址基本已完成，部分建在深厚覆蓋層上的大壩在采取適當措施后，建設和運行條件良好[1- 4]。相關文獻[5- 7]對建在深覆蓋層上的大壩進行了有益的探索和實踐。隨著國家清潔能源戰略的持續推進，一大批擬建的土石壩將修筑于深厚覆蓋層上[8- 9]。

準確獲取覆蓋層上堆石壩的動力響應，進行抗震穩定分析是保證工程安全的重要手段之一[10- 13]。本文以某面板堆石壩為例，系統研究了壩體加速度、動位移、安全系數等動力反應，總結了高堆石壩地震響應規律。

1 工程概況

該工程位于黃河一級支流的最后一段峽谷出口處，是黃河下游防洪體系的重要組成部分。擋水建筑物采用混凝土面板堆石壩，壩址處岸坡陡峻，河谷呈“U”型。壩頂高程288.5m，最大壩高122.5m，壩頂長度481.0m，壩頂寬10.0m，上、下游壩坡分別為1∶1.5和1∶1.6。趾板置于覆蓋層上，布置在面板的周邊，與防滲面板通過設有止水的周邊縫連接，形成壩基以上的防滲體，河床趾板上游壩基采用混凝土防滲墻截滲。趾板與防滲墻之間采用連接板連接，連接板長度4.0m，厚度0.9m。防滲墻兩端和底部嵌入基巖0.5m，為防止產生過大的應力集中，嵌入處設置石渣柔性支座。典型剖面如圖1所示。

大壩基礎座落在含漂石及泥的砂卵石覆蓋層上，覆蓋層平均厚度約30m，最大厚度41.87m。根據河床鉆孔資料，在壩軸線附近及下游內存在14個砂層透鏡體，其中6個透鏡體分布在地面以下8m以內，4個分布在地面20m以下。覆蓋層中發現4層較連續的黏性土夾層，最厚達12m，順河延伸350～800m，對壩基穩定、變形起控制作用。

根據地質勘查資料，黏性土夾層累計厚度5～20m，占覆蓋層總厚度的1/6～1/2，壓縮系數為0.1～0.2MPa-1，屬中低壓縮性土；砂卵礫石層的壓縮系數為0.01～0.068MPa-1，屬低壓縮～不可壓縮；根據標貫擊數及相對密度，砂層透鏡體相當于中密～密室。

2 動力計算

2.1 有限元模型

根據壩體分區、施工及加載過程，考慮防滲墻的連接型式，剖分大壩的三維有限元網格。模型以8節點六面體等參單元為主，單元總數為9625，節點總數為11302。

2.2 動力本構模型及參數

動力計算采用等效線性粘彈性模型，即假定壩體堆石料和壩基卵石層為粘彈性體，采用等效剪切模量G和等效阻尼比λ這兩個參數反映筑壩料的非線性和滯后性，并表示為剪切模量與阻尼比與動剪應變的關系[14- 16]。其中，最大動剪模量Gmax根據靜力狀態確定，可表示為：

圖1 大壩典型剖面

(1)

表1 筑壩料最大動剪模量參數

有限元動力分析時，動模量G和阻尼比λ通過對室內動力三軸試驗得到的G/Gmax及λ與動剪應變γ的關系插值獲得，動三軸試驗結果見表2。

2.3 地震波輸入

根據地震局提供的壩址地震資料，基巖輸入地震動采用超越概率100年2%的峰值強度為201gal，地震動持續時間取為24s。人工合成地震波時程如圖2所示。地震波輸入方向：x為沿順河向水平輸入；y為豎向輸入，根據水工建筑物抗震設計規范，將峰值折減2/3；z方向為壩軸向輸入。

圖2 地震動時程曲線

3 抗震穩定性分析

通過有限元計算，可求出各單元的應力情況，根據摩爾-庫倫破壞準則，將局部安全系數小于1的區域連接起來，即可判斷出最危險復合滑動面的位置。在滑動面上，用總抗滑力和總滑動力的比值確定最小安全系數Fs，min。通過分析整個地震持續時間內，安全系數與時間的關系，評價壩體的抗震穩定性。

表2 筑壩料在不同動剪應變下的和

計算出壩體內各單元的靜應力和動應力后，局部安全系數可按下式計算：

(2)

式中：σ1、σ3—任意時刻的最大、最小主應力。

4 計算結果分析

4.1 加速度分析

地震期間，壩體典型剖面D0+170的加速度等值線如圖3所示。壩體地震反應呈現典型的“鞭鞘效應，即隨著壩高增加，加速度逐漸增大，在壩頂附近達到最大。其中，順河向加速度最大值為7.5m/s2，放大系數為3.75。豎向加速度最大值為7m/s2，放大系數為3.5。

圖3 壩體加速度反應等值線

4.2 動位移反應

圖4 壩體動位移等值線(m)

典型剖面動位移等值線見圖4。其中，順河向最大位移為10cm，豎向最大動位移為4.5cm，最大值均發生在下游壩頂附近。

4.3 抗震穩定分析

當t=10s，t=18s和t=22s時，壩體安全系數等值線如圖5所示，壩體內單元安全系數均大于1。從分布規律上看，由于受水荷載作用，上游側壩體的安全系數高于下游壩體。

圖5 某瞬時壩體安全系數分布

為了分析大壩的抗震安全性，在壩基、1/4和3/4高程的壩軸線處選取監測點，給出整個地震時程內的安全系數，如圖6所示。地震期間，壩基和1/4壩高處，單元安全系數均大于1。在3/4壩高處，存在短時間內安全系數小于1的情況，然而安全系數小于1的持續時間很短，且地震慣性力為往復荷載，單元內應力在超過抗剪強度后，慣性力轉向，安全系數立即增大。因此，認為壩體抗震穩定能夠滿足要求。

圖6 壩體內某單元安全系數時程

5 結論

(1)隨著壩高增加，加速度逐漸增大，在壩頂附近達到最大。其中，順河向加速度最大值為7.5，放大系數為3.75。豎向加速度最大值為7，放大系數為3.5。

(2)壩體地震反應呈現典型的鞭鞘效應。在壩頂附近，加速度、動位移等地震響應都達到最大。因此，壩頂是堆石壩抗震安全的薄弱部位，應在壩頂附近采取緩坡、加筋、漿砌石護坡等加固方案。

(3)總體上，壩體內單元安全系數大于1。從分布規律上看，由于受水荷載作用，上游側壩體的安全系數高于下游壩體。

(4)根據壩體單元局部安全系數可知，在3/4壩高處，短時間內安全系數小于1。然而持續時間很短，且地震慣性力為往復荷載，單元內應力在超過抗剪強度后，慣性力轉向，安全系數立即增大。因此，大壩整體抗震穩定能夠滿足安全要求。