地下水和地震共同作用下的邊坡穩(wěn)定性分析

黃 帥 汪洪祥 王 榮 蔡德鉤 閆宏業(yè)

(1.中國地震局地殼應(yīng)力研究所，北京100085;2.中國港灣工程有限責(zé)任公司，北京100027;3.中國鐵道科學(xué)研究院，北京100081)

當(dāng)邊坡內(nèi)含有較高的地下水位時，其產(chǎn)生較大的孔隙水壓力，導(dǎo)致邊坡有效應(yīng)力降低，對邊坡的穩(wěn)定性產(chǎn)生不利影響。尤其是我國西部地區(qū)年降水量表現(xiàn)出普遍增加的趨勢，在降雨集中時期，會導(dǎo)致邊坡地下水位的升高，對邊坡的穩(wěn)定性極為不利。除了西部地區(qū)覆蓋大量砂質(zhì)土外，在已經(jīng)修建的高鐵沿線也存有大量的砂質(zhì)土，特別是濕地地區(qū)，例如京滬高鐵沿線經(jīng)過黃河濕地，存有大量的砂質(zhì)土邊坡。由于砂質(zhì)土邊坡滲透性比一般的黏性土大，在滲透力作用下極易引發(fā)邊坡失穩(wěn)破壞，嚴(yán)重影響著邊坡的穩(wěn)定性?；诖?，研究不同地下水位下邊坡穩(wěn)定性對鐵路的安全運營和建設(shè)具有重要的意義。

國內(nèi)外就不同的地下水位對邊坡地震動力響應(yīng)和穩(wěn)定性的影響進行了相關(guān)研究。國外的G. W.Jia［1］基于粉質(zhì)砂土邊坡的原型試驗，研究了地下水位上升和下降對邊坡破壞模式的影響規(guī)律，得出隨著坡外水位的下降，坡內(nèi)孔隙水壓力表現(xiàn)出明顯的滯后現(xiàn)象，產(chǎn)生朝向坡外較大的外推力，嚴(yán)重影響邊坡的穩(wěn)定性;Griffith［2］和Lane［3］基于強度折減法研究了水位變化對邊坡安全系數(shù)的影響規(guī)律，得出安全系數(shù)隨著水位增加表現(xiàn)出先減小后增大的趨勢;Dawson［4］則利用有限差分法和強度折減法對飽水土質(zhì)邊坡進行了計算。國內(nèi)的黃帥等［5］基于時程分析法研究了地下水位變化對邊坡地震穩(wěn)定性的影響;戴繪［6］研究了降雨過程與雨停后浸潤線變化規(guī)律;周靖人［7］研究了不同降雨強度下，地下水位埋深的變化對厚沖積層邊坡的安全系數(shù)的影響規(guī)律。然而這些研究大多針對地下水位單一荷載的影響，而較少考慮地震荷載的影響，而邊坡的防護設(shè)計中地震荷載是必須考慮的因素。因此，本研究基于邊坡的有限元模型考察地下水位變化和地震共同作用下邊坡的穩(wěn)定性。

1 邊坡地震穩(wěn)定性的數(shù)值仿真分析

1.1 有限元模型

邊坡的高度為12 m，坡度為35°，建立的仿真模型如圖1 所示，為保證計算的精度，本模型劃分的網(wǎng)格最大尺寸小于地震波最短波長的1/10 ～1/8。邊坡所在地區(qū)基本烈度為7 度，設(shè)計基本加速度為0.1g。模型采用莫爾－庫侖本構(gòu)模型。左側(cè)和右側(cè)邊界采用水平位移約束，底部邊界采用豎向位移約束，其他邊界為自由邊界。為了分析在地下水位變化對邊坡不同位置的位移、應(yīng)力、應(yīng)變和孔隙水壓力的影響規(guī)律，定義邊坡模型底部為零水位線，邊坡水位線最高位置為地下水位高度。

圖1 邊坡模型Fig.1 Model of the slope

基于室內(nèi)直剪試驗測的邊坡巖土(砂土)的物理參數(shù)如表1 所示。

表1 邊坡材料參數(shù)Table 1 Material parameters of slope

1.2 地震波的選取

邊坡所在地區(qū)基本烈度為7 度，場地為II 類場。參考日本《道路橋示方書－同解說v. 耐震設(shè)計篇》(2002)，選用遠(yuǎn)場地震波(Type I:T1－II－1 和T1－II－3)和近場地震波(Type II:T21－II－1 和T2－II－3)以及EIcentro 地震波進行分析。如圖2 所示。

圖2 地震波時程曲線Fig.2 Time history of the earthquake wave

2 地下水位對邊坡穩(wěn)定性影響分析

2.1 對邊坡抗剪性能的影響分析

對邊坡在不同地下水位下的剪應(yīng)力和剪應(yīng)變進行了分析，如圖3 ～圖5 所示(選取20 m 水深為例)。

圖3 EICentro 地震下邊坡剪切應(yīng)力應(yīng)變Fig.3 The shear stress and strain of slope under EICentro

圖4 T1－II－1 地震下邊坡剪切應(yīng)力應(yīng)變Fig.4 The shear stress and strain of slope under T1－II－1

圖5 T2－II－1 地震下邊坡剪切應(yīng)力應(yīng)變Fig.5 The shear stress and strain of slope under T2－II－1

由圖3 ～圖5 可知，在地震作用下，坡腳的剪應(yīng)力－剪應(yīng)變曲線呈現(xiàn)出滯回特性。隨著地震的持續(xù)作用，由于邊坡坡腳發(fā)生了明顯的塑性變形，使得其滯回曲線的平衡位置不斷地改變。隨著地下水位的升高，邊坡的剪應(yīng)力－剪應(yīng)變滯回圈表現(xiàn)出增大的趨勢，即邊坡的剪應(yīng)力－剪應(yīng)變滯回圈向遠(yuǎn)離坐標(biāo)中心方向延伸，說明地下水的存在對坡腳的抗剪能力有較大影響。主要原因是地下水的存在使得土層含水率增大，坡腳位置為地下水滲出位置，土體在滲透力和孔隙水壓力的共同影響下有效應(yīng)力降低，地震發(fā)生時邊坡極易發(fā)生剪切破壞。在不同地下水位下，遠(yuǎn)場地震對邊坡坡腳抗剪能力的影響要大于近場地震。

2.2 對邊坡主應(yīng)力的影響分析

提取了沿邊坡高度的最大主應(yīng)力，如圖6 ～圖8所示(選取22 m 水深為例)。

圖6 EICentro 地震作用下邊坡的第一有效主應(yīng)力Fig.6 The first principal effective stress of slope under EICentro

圖7 T1－II－1 地震作用下邊坡的第一有效主應(yīng)力Fig.7 The first principal effective stress of slope under T1－II－1

圖8 T2－II－1 地震作用下邊坡的第一有效主應(yīng)力Fig.8 The first principal effective stress of slope under T2－II－1

由圖6 ～圖8 可知，隨著地下水位的升高，在不同類型地震作用下邊坡的第一有效主應(yīng)力整體呈現(xiàn)增大的趨勢，EICentro 地震作用下水位高度24 m 時的第一有效主應(yīng)力在監(jiān)測點A、C、E、G、H、I、J 位置(如圖1(b))比0 m 水位(圖件從略)時分別增大了8.7%、3.27%、1.75%、1.50%、1.53%、1.50% 和0.4%。遠(yuǎn)場地震T1－II－1 作用下，監(jiān)測點A、C、E、G、H、I、J 位置的第一有效主應(yīng)力比0 m 水位時分別增大了10.53%、5.69%、3.19%、25.04%、10.48%、10.94%和9.88%;T1－II－3 作用下邊坡不同監(jiān)測點的第一有效主應(yīng)力在監(jiān)測點A、C、E、G、H、I、J 位置比0 m 水位時分別增大了4.0%、2.12%、0.8%、30.59%、11.57%、15.49%和21.24%。近場地震T2－II－1 作用下，監(jiān)測點A、C、E、G、H、I、J 位置的第一有效主應(yīng)力比0 m 水位時分別增大了0.21%、0.7%、11.46%、38.10%、19.96%、17.68% 和18.31%;T2－II－3 作用下監(jiān)測點A、C、E、G、H、I、J 位置的第一有效主應(yīng)力比0 m 水位時分別增大了9.40%、1.11%、8.84%、37.79%、24.42%、14.57%和7.46%。由此可知，隨著地下水位的增加，邊坡的第一有效應(yīng)力減小，在近場地震作用下其減小幅度要大于遠(yuǎn)場地震和EICentro 地震作用下的減小幅度;而在不同地下水位下，遠(yuǎn)場地震作用的邊坡第一有效主應(yīng)力要大于近場地震作用，說明隨著地下水位的增加，近場地震作用對邊坡穩(wěn)定性影響更大。

3 地下水位對邊坡動力安全系數(shù)的影響分析

計算了邊坡在T1－II－1、T1－II－3、T2－II－1 和T2－II－3 作用下無地下水及水深H=20 m 時的動力安全系數(shù)，并與《GB50111—2006 鐵路工程抗震設(shè)計規(guī)范》中的擬靜力法計算的邊坡安全系數(shù)進行了對比分析，如圖9 所示。

圖9 近遠(yuǎn)場地震作用下邊坡的安全系數(shù)Fig.9 The safety factor of the slope under near and far field earthquake

由圖9 可以看出，地震作用下邊坡的安全系數(shù)隨著地震的持續(xù)作用上下波動，而規(guī)范中的擬靜力法計算的邊坡安全系數(shù)為一定值。在邊坡無水時，擬靜力法計算的安全系數(shù)值位于動力安全系數(shù)波動值的中部偏下;有地下水影響時，隨著地下水位的升高，擬靜力法計算的安全系數(shù)值開始移向動力安全系數(shù)波動值的中部偏上，且動力安全系數(shù)的最小值比擬靜力法的靜力安全系數(shù)小。H=20 m 時，動力安全系數(shù)的最小值比擬靜力法的安全系數(shù)在T1－II－1 地震作用時減小了21.30%，T2 －II －1 地震作用時減小了18.13%。因此，在有地下水影響時，基于擬靜力法計算的邊坡安全系數(shù)進行邊坡地震穩(wěn)定性評價時有一定的風(fēng)險。地下水位H=20 m 時，邊坡的動力安全系數(shù)比無水時明顯降低，在T1－II－1 地震作用時最大減小了27.60%，T2－II－1 地震作用時減小了26.98%。故地下水的存在對邊坡穩(wěn)定性具有較大的影響。

4 結(jié) 論

(1)近遠(yuǎn)場地震作用下，邊坡坡腳的剪應(yīng)力剪應(yīng)變表現(xiàn)出滯回性能，且隨著地震的持續(xù)作用，由于邊坡發(fā)生了明顯塑性變形，其滯回平衡位置不斷改變。

(2)隨著水深的增加，地下水對邊坡坡腳抗剪能力的影響逐漸增大，且遠(yuǎn)場地震作用時，邊坡坡腳的剪應(yīng)力－剪應(yīng)變滯回延伸較大，說明遠(yuǎn)場地震對邊坡坡腳的抗剪能力的影響要大于近場地震。

(3)地下水位越高，對邊坡安全系數(shù)的影響越顯著。在邊坡無地下水時，擬靜力法計算的安全系數(shù)值位于動力安全系數(shù)波動值的中部偏下;而隨著地下水位的升高，擬靜力法計算的安全系數(shù)曲線開始移向動力安全系數(shù)波動值的中部偏上。因此，當(dāng)?shù)叵滤簧邥r，擬靜力法計算的邊坡安全系數(shù)進行邊坡穩(wěn)定性評估將具有一定的潛在風(fēng)險性。

［1］ Jia G W，Tony L T，Zhan Y M，et al. Performance of a large-scale slope model subjected to rising and lowering water levels［J］.Engineering Geology，2009，106(2):92-103.

［2］ Griffiths D V，Lane P A. Slopes stability analysis by finite elements［J］.Geotechnique，1999，49(3):388-403.

［3］ Lane P A，Griffiths D V. Assessment of stability of slopes under drawdown conditions［J］. Journal of Geotechnical and Geoenviromental Engineering，2000，126(5):443-450.

［4］ Dawson E M，Koth W H，Drescher A. Slope stability analysis by strength reduction［J］.Geotechnique，1999，49(6):835-841.

［5］ Huang Shuai，Yan Hongye，Cai Degou，et al.Dynamic response analysis of soil slope under strong earthquake［J］.Advance Material Research，2014，920:921-924.

［6］ 戴 繪，秦衛(wèi)星，張岳安.強降雨作用下山谷型尾礦壩浸潤線演化規(guī)律［J］.金屬礦山，2013(7):149-152.

Dai Hui，Qin Weixing，Zhang Yuean.Evolution law of phreatic line in valley-type tailing s dam under heavy rainfall［J］. Metal Mine，2013(7):149-152.

［7］ 周靖人，楊天鴻，于慶磊，等. 降雨對厚沖積層邊坡穩(wěn)定性影響［J］.金屬礦山，2014(9):11-14.

Zhou Jingren，Yang Tianhong，Yu Qinglei，et al.Impact of rainfall on stability of thick alluvium slope［J］.Metal Mine，2014(9):11-14.