米熱尼沙·麥麥提

(喀什地區(qū)水利水電勘測設(shè)計院，新疆 喀什 844000)

1 概 述

中國是世界上地震頻發(fā)的國家之一，地震災(zāi)害給人民的財產(chǎn)安全和生命安全造成巨大的威脅[1]。由地震引發(fā)的滑坡災(zāi)害就是其中之一，而降雨又加速了這種災(zāi)害的發(fā)生概率。例如1999年的臺灣集集大地震[2]后的強(qiáng)降雨誘發(fā)了新的滑坡等地質(zhì)災(zāi)害，降雨誘發(fā)的滑坡面積是地震直接誘發(fā)的滑坡面積的3倍；汶川大地震也誘發(fā)了1.5萬處滑坡[3]，說明震后的降雨將會誘發(fā)更大范圍的滑坡。國內(nèi)外許多專家學(xué)者都對地震、降雨或者降雨耦合地震工況下的邊坡穩(wěn)定性進(jìn)行大量的研究。Premchit[4]在1991年通過研究發(fā)現(xiàn)，在香港24 h降雨量超過200 mm時一定會有滑坡發(fā)生。Keefer[5]1984年通過對40個歷史地震進(jìn)行研究分析得出, 地震震級不大于4級時不會有滑坡發(fā)生；當(dāng)震級達(dá)到9.2級時, 將會在500 000 km2范圍內(nèi)發(fā)生滑坡。徐文杰[6]等選用ABAQUS軟件, 從三維空間出發(fā)對地震誘發(fā)肖家橋滑坡的失穩(wěn)進(jìn)行了相關(guān)的模擬分析, 克服了二維分析方法帶來的誤差, 且與現(xiàn)實較符。孫軍杰[7]認(rèn)為某種特定條件下 (降雨充沛地區(qū)發(fā)生特大地震時) , 地震和降雨的持續(xù)影響效應(yīng)可能會對邊坡體穩(wěn)定性產(chǎn)生耦合削弱作用, 由此造成的與滑坡有關(guān)的潛在災(zāi)害不容忽視。

為了分析研究降雨及地震耦合作用對土石壩邊坡穩(wěn)定性的影響, 本文以某水利樞紐中的黏土心墻壩為工程背景, 采用Geostudio軟件模擬分析土石壩壩坡在不同庫水位下地震與降雨耦合情況下的穩(wěn)定性, 為今后類似工程研究提供科學(xué)依據(jù)。

2 理論背景

2.1 非飽和理論

庫水位驟降時，上下游堆石區(qū)呈現(xiàn)出典型的飽和-非飽和滲流特性，特別是上游堆石區(qū)內(nèi)的孔壓力出現(xiàn)較大的波動，飽和-非飽和滲流[8]的張量形式如下：

(1)

2.2 非線性材料模型

本文進(jìn)行動力模擬時，采用Geostudio中Quake/w中的非線性材料模型。本文主要研究對象為上下游堆石區(qū)，QUAKE/w提供的粗粒土的剪切模量表達(dá)式如下：

(2)

(3)

2.3 壩坡永久變形

本文采用1965年Newmark[9]提出的地震下邊坡的永久變形分析方法，即同時計算多個潛在滑動面并且計算出最大變形邊坡滑動面，在獲得最大變形的滑動面后Geoslope計算出安全系數(shù)Fs為1時的屈服加速度ay以及該潛在滑動面隨時間變化的平均加速度a，對(a-ay)關(guān)于時間二次積分即可得到潛在滑動面的Newmark邊坡永久位移大小。

3 計算模擬

3.1 工程概況

某水利樞紐工程中的攔河壩為黏土心墻土石壩，從上游到下游依次為上游圍堰、上游堆石區(qū)、過渡層、黏土心墻、下游堆石區(qū)以及下游堆石棱體。壩頂高程為156.20 m，正常蓄水位為152.0 m，死水位為125 m，典型剖面見圖1。

圖1 心墻壩典型剖面圖

3.2 計算模型

模型網(wǎng)格圖見圖2。地基厚度為50 m，長305 m，地面高程為0 m，則死水位為15 m，正常蓄水位為42 m。N1、N2、N3上游堆石區(qū)內(nèi)不同位置的監(jiān)測點。網(wǎng)格單元類型主要為三角形和四邊形。為了更準(zhǔn)確反映上下游堆石體內(nèi)不同監(jiān)測點的孔隙水壓力變化，取上下游堆石體和心墻部分的網(wǎng)格尺寸為1.5 m，其余部分取3.5 m，全局共劃分3 854個單元，3 896個節(jié)點。

3.3 土體參數(shù)

黏土心墻土石壩各區(qū)域的物理力學(xué)參數(shù)見表1。土水特征曲線采用Fredlund&Xing模型， 描述土水特征曲線 (SWCC曲線) 是衡量土體內(nèi)部滲透系數(shù)(體積含水率)與基質(zhì)吸力之間的關(guān)系的重要曲線，見圖3。

圖2 模型網(wǎng)格圖

材料干密度/kg·(m3)-1最大剪切模量 /GPa泊松比黏聚力/kPa內(nèi)摩擦角/(°)滲透系數(shù)/cm·s-1阻尼比心墻1 61050.3086263×10-60.05堆石區(qū)1 50050.2880345.79×10-30.03過渡層1 55050.3348267.08×10-30.03壩基土1 60050.2553203.53×10-30.03上游圍堰1 58050.3537.47285.23×10-30.03排水棱體1 68050.2550342.11×10-10.03基巖1 67050.3535241.00×10-90.06

圖3 土水特征函數(shù)曲線

根據(jù)可研性研究報告，庫區(qū)的地震動峰值加速度為0.10 g，相應(yīng)的地震烈度為Ⅶ度。本次數(shù)值模擬采用國際上常用的EI-Centro水平方向的地震加速度時程曲線，為縮短計算時間截取地震振幅較大的前10 s作為輸入波，在對地震波進(jìn)行濾波以及基線校正處理的基礎(chǔ)上，將其峰值調(diào)為0.10 g(g=10 m/s2)，圖4為校正后的加速度時程曲線。

圖4 校正后水平向加速度時程曲線

4 計算結(jié)果分析

4.1 不同驟降速率上下游堆石區(qū)孔隙水壓力分析

本文研究的是某黏土心墻土石壩上下游堆石區(qū)在不同庫水位工況下發(fā)生降雨時的滲透特性，以及在降雨過程中遇到地震作用下的上下游壩坡穩(wěn)定性。降雨強(qiáng)度為0.1m/d，降雨持時為2天，發(fā)生在0～2天；地震最大輸入加速度為0.1 g，持續(xù)時間為10 s，發(fā)生時刻為第一天降雨結(jié)束時刻；數(shù)值模擬時長為10天。N1、N2、N3、N4、N5和N6分別為上下游堆石區(qū)內(nèi)的監(jiān)測點。從圖5可以看出，在死水位工況下，上游堆石區(qū)內(nèi)的監(jiān)測點孔壓總體上呈現(xiàn)先增大后減小的變化規(guī)律。由于監(jiān)測點離上游庫水位的距離不同，各個監(jiān)測點的孔壓在降雨過程中的變化也存在差異。監(jiān)測點離庫水位越遠(yuǎn)，初始孔壓力越小，降雨后孔壓變化越劇烈。例如，N1、N2和N3的初始孔壓力分別為-265.6、-176.7和-18.4 kPa；在降雨發(fā)生的一天內(nèi)，最大孔壓力上升到-13.6、-12.2和-9.7 kPa,分別增大95%、93%和34%。這是因為降雨強(qiáng)度比較大，上游堆石區(qū)的滲透系數(shù)也比較大，大部分降雨入滲到土體內(nèi)，導(dǎo)致土體的孔壓力迅速升高。監(jiān)測點發(fā)生最大孔壓力時刻并非停雨時刻，而是在降雨期間。停雨后，各監(jiān)測點的孔壓力緩慢減小，直至計算結(jié)束時刻，分別為-41.8、-41.3和-17.1 kPa。對于下游堆石區(qū)的監(jiān)測點來說，監(jiān)測點離壩體浸潤線的距離越遠(yuǎn)，初始孔壓力越小，降雨后孔壓力變化越劇烈，達(dá)到孔壓力最大值的時間越長。例如，N4、N5和N6在降雨期間孔壓力達(dá)到最大值的時間分別為0.25、1.75和2天。停雨后，監(jiān)測點的孔壓力緩慢減小，計算結(jié)束時刻孔壓力相差不大。對比上下游同一高度的監(jiān)測點孔壓力可以看出，下游監(jiān)測點的初始孔壓比上游相同高度的監(jiān)測點小，而且降雨期間孔壓力的變化更加劇烈。

圖5 死水位下監(jiān)測點孔壓力圖

由圖6可以看出，在正常蓄水位工況下，上游堆石區(qū)內(nèi)監(jiān)測點的孔壓總體上保持不變。這是因為當(dāng)庫水位為正常蓄水位時，上游監(jiān)測點都位于飽和區(qū)，降雨對其幾乎沒有影響。對于下游各監(jiān)測點來說，其變化規(guī)律與死水位工況下的相似。對比不同庫水位下下游同一監(jiān)測點可以看出，死水位下監(jiān)測點的孔壓力小于正常蓄水位，但是到達(dá)孔壓最大值的時間相對較短，計算結(jié)束時刻孔壓力更大。

圖6 正常蓄水位下監(jiān)測點孔壓力圖

4.2 不同庫水位工況下降雨偶遇地震上下游堆石區(qū)抗震穩(wěn)定分析

地震可能發(fā)生于降雨過程中的任意時刻，所以限于篇幅，本文假定降雨1天結(jié)束時刻發(fā)生地震。從圖7(a)可以看出，總體上正常蓄水位工況下上游壩坡安全系數(shù)比死水位工況下大，而且最小安全系數(shù)均大于1。由此可見，庫水位越高，上游壩坡在地震作用下的穩(wěn)定性越好。從圖7(b)可以看出，正常蓄水位和死水位工況下，上游壩坡的永久變形分別為0和0.025 m。

圖7 上游動力響應(yīng)圖

從圖8(a)可以看出，對于下游壩坡而言，死水位工況下的壩坡安全系數(shù)比高水位下大，高水位和低水位下下游壩坡永久變形分別為0.241和0.039 m。由此可見，庫水位越高，下游壩坡在地震作用下的穩(wěn)定性越差。對比上下游壩坡安全穩(wěn)定系數(shù)曲線可得，庫水壓力可以加強(qiáng)上游壩坡的穩(wěn)定性，但是水位越高，下游壩坡的穩(wěn)定性越差。這是由于庫水位越高，作用在坡面的水壓力越大，增加了上游壩坡的抗滑力。但是庫水位越高，下游壩坡的地下水位線越高，下游堆石區(qū)的飽和區(qū)越大，土體的抗剪強(qiáng)度減小，下游壩坡的抗滑力減小。綜上可知，低水位下黏土心墻土石壩的上游邊坡為抗震穩(wěn)定薄弱區(qū)，高水位下黏土心墻土石壩的下游邊坡為抗震穩(wěn)定薄弱區(qū)。

圖8 下游動力響應(yīng)圖

5 結(jié) 論

1) 在不同水位工況下發(fā)生降雨時，上下游堆石區(qū)內(nèi)的監(jiān)測點孔壓總體上呈現(xiàn)先增大后減小的變化規(guī)律，上游堆石區(qū)內(nèi)的監(jiān)測點離庫水位越遠(yuǎn)，孔壓力變化越劇烈；下游堆石區(qū)內(nèi)的監(jiān)測點離地下水位越遠(yuǎn)，孔壓力變化越劇烈。

2) 在降雨期間遇到地震時，高水位工況下，上游壩坡的安全系數(shù)總體上比低水位工況大，而下游壩坡恰恰相反。