蔡德鉤，黃 帥，閆宏業(yè)，陳 鋒，張千里，姚建平

(1.高速鐵路軌道技術國家重點實驗室，北京 100081;2.中國鐵道科學研究院 鐵道建筑研究所，北京 100081; 3.中國地震局地殼應力研究所，北京 100085)

地下水位上升對邊坡穩(wěn)定性影響的擬靜力分析

蔡德鉤1，2，黃 帥3，閆宏業(yè)1，2，陳 鋒1，2，張千里1，2，姚建平1，2

(1.高速鐵路軌道技術國家重點實驗室，北京 100081;2.中國鐵道科學研究院 鐵道建筑研究所，北京 100081; 3.中國地震局地殼應力研究所，北京 100085)

利用有限元模型分析了不同影響因素(坡角、黏聚力、內(nèi)摩擦角、地下水位高度、設計峰值加速度)對邊坡穩(wěn)定性影響的顯著性，并研究了地下水位變化對邊坡穩(wěn)定性和變形的影響規(guī)律。結(jié)果表明:內(nèi)摩擦角對邊坡穩(wěn)定性影響最大，其次為設計峰值加速度、坡角、黏聚力、地下水位高度。地下水位的變化可使邊坡安全系數(shù)降低15%，因此地下水的影響不容忽略。當水位達到較高位置(地下水位高度達22 m)時，水平最大位移比無水時增加58.33%。隨著地下水位的升高邊坡的最大水平位移開始表現(xiàn)出邊坡后緣下挫并逐步向前緣擴展的趨勢，使得邊坡的變形最大值由邊坡的最高處開始向坡腳處轉(zhuǎn)移。當?shù)叵滤簧咭欢ǜ叨葧r，邊坡的最大水平位移會出現(xiàn)一個突變，位移最大區(qū)域向邊坡內(nèi)部擴展，說明地下水的影響使得邊坡滑動面向邊坡內(nèi)部延伸，其滑動規(guī)模增大，這可作為邊坡滑坡預測的重要指標。

砂質(zhì)邊坡 地下水位 擬靜力 安全系數(shù) 位移

由于地表水的滲入、施工生活用水隨意排放、雨季施工等原因引起土質(zhì)高邊坡地下水位上升，可能導致水位上升前較穩(wěn)定的坡體產(chǎn)生滑塌。根據(jù)對大量邊坡失穩(wěn)原因的分析，幾乎所有的崩塌、滑坡都與水有關，特別是地下水的存在及變化往往是形成崩塌、滑坡等邊坡失穩(wěn)的主要條件［1-3］。如瑞典、挪威和加拿大的天然滑坡，大都是由于覆蓋層中地下水位的增加引起孔隙水壓力增大而誘發(fā)的。因此，在實際工程中應特別重視工程的水環(huán)境，控制地下水位的超標抬升。

目前邊坡的穩(wěn)定性分析方法主要包括擬靜力法、Newmark滑塊位移法、時程分析法和有限元分析法。擬靜力法簡便、易于操作，在工程實踐中被廣泛應用。此外，擬靜力法一般結(jié)合邊坡極限平衡法和強度折減法進行邊坡的穩(wěn)定性分析。國內(nèi)外對地下水位變化對邊坡穩(wěn)定性的影響進行了大量的研究。Griffith、Lane等［4-5］基于自己開發(fā)的有限元軟件，利用強度折減法分析了水位變化對邊坡安全系數(shù)的影響。賈官偉等［6］研究了水位驟降引起臨水邊坡滑坡的原因及失穩(wěn)模式。趙煉恒等［7］利用安全系數(shù)指標對地下水位變化影響下的均質(zhì)邊坡進行了穩(wěn)定性分析。黃帥等［8-12］基于時程分析法研究了地下水位變化對邊坡穩(wěn)定性的影響。然而大部分研究主要集中在地下水位變化對邊坡安全系數(shù)的影響，而較少就地下水位變化對邊坡的敏感性、破壞模式和位移的影響進行研究?；诖?，本文采用擬靜力法進行邊坡的穩(wěn)定性分析。研究地下水位變化對邊坡安全系數(shù)、變形和破壞模式的影響，明確水位變化對邊坡穩(wěn)定性和變形的影響規(guī)律。

1 砂質(zhì)邊坡有限元模型的建立

1.1 工程概況

本文以某鐵路CKl62+075—CKl63+390段為研究對象，該段鐵路處于鐵路終端，低山區(qū)路基深路塹挖方集中地段，不同形式的路塹主要集中在此段。選其典型的路塹斷面為研究對象，坡長為17 m，設計高度約為12 m，坡度35°，為砂質(zhì)土。鐵路沿線以黃綿土、紅土、風砂土為主。水位一般在2.0 m以上，水位隨季節(jié)的變化而變化。

1.2 邊坡的有限元模型

邊坡的高為12 m，邊坡坡度為35°，建立邊坡的仿真模型如圖1所示，圖中H為地下水位高度。為保證計算的精度，有限元分析模型網(wǎng)格的最大尺寸小于輸入地震波最短波長的1/10～1/8。模型采用莫爾—庫侖本構(gòu)模型。有限元模型底部和左右邊界采用黏彈性動力人工邊界，其它邊界為自由邊界。

為了分析地下水位變化對邊坡不同位置位移的影響規(guī)律，分別對邊坡不同位置設置了監(jiān)測線，如圖2所示。采用篩分試驗法對砂土顆粒進行級配分析，測出邊坡所用砂土的顆粒分布特征如圖3所示。通過室內(nèi)試驗測定的砂土的主要物理參數(shù)如表1所示。

圖1 邊坡計算模型

圖2 邊坡監(jiān)測線

圖3 砂土的篩分試驗及顆粒分布曲線

表1 試驗砂質(zhì)土的物理參數(shù)

2 邊坡穩(wěn)定性影響因素的敏感性分析

對邊坡而言，不同參數(shù)對其穩(wěn)定性影響的顯著性不同，故有必要確定各參數(shù)對邊坡穩(wěn)定性影響的敏感性大小。在實際問題中考慮1個因素或2個因素對計算結(jié)果的顯著性分析，可以選用一元或二元方差分析。而本文邊坡的穩(wěn)定性需考慮多個因素對其穩(wěn)定性的影響，可采用正交試驗的方法進行分析。研究不同參數(shù)對邊坡穩(wěn)定性的影響規(guī)律，具體計算方案如表2所示。

表2 各因素取值和水平

為了明確各參數(shù)對邊坡穩(wěn)定性的影響規(guī)律，采用正交試驗進行數(shù)值分析，研究坡角(A)、黏聚力(B)、內(nèi)摩擦角(C)、地下水位高度(D)、設計峰值加速度(E)等因素對邊坡穩(wěn)定性的影響規(guī)律，即以邊坡穩(wěn)定安全系數(shù)為指標進行多因素單指標計算分析。進行16次正交試驗，即以L16(45)正交試驗表(表3)進行數(shù)值分析。通過數(shù)值分析求出各因素不同組合下邊坡的安全系數(shù)，并對其進行極差分析，如圖4所示。

由圖4(a)可知，邊坡的安全系數(shù)隨著黏聚力、內(nèi)摩擦角的增大均表現(xiàn)出增大的趨勢，其中黏聚力和內(nèi)摩擦角對安全系數(shù)的影響幅度最大。而隨著邊坡坡角、地下水位高度和設計峰值加速度的增加邊坡的安全系數(shù)減小，且地下水位由水平1增加到水平4時，安全系數(shù)降低了15%，說明地下水位的影響不可忽略。由圖4(b)可知，邊坡穩(wěn)定安全系數(shù)計算值最大時的方案為A2B3C4D1E2，即在此水平下計算的邊坡安全系數(shù)最大，此時邊坡最穩(wěn)定。邊坡穩(wěn)定安全系數(shù)計算值最小時的方案為A2B2C1D4E3，即在此方案下邊坡的穩(wěn)定性最差。

圖4(b)可知，極差從大到小的順序依次為RC，RE，RA，RB，RD?？芍诳紤]的邊坡各參數(shù)中內(nèi)摩擦角對邊坡穩(wěn)定性影響最大，其次為設計峰值加速度、坡角、黏聚力、地下水位高度。盡管相對其它4個因素，地下水對邊坡安全系數(shù)影響較小，但是其它4個因素，都是邊坡穩(wěn)定性影響的敏感因素，且在進行計算時并沒有考慮地震作用下的動孔隙水壓力影響，且從圖4(a)也可以看出，地下水位的變化對邊坡安全系數(shù)的影響比較大，因此地下水的影響不容忽略。

表3 正交試驗數(shù)值分析

圖4 邊坡穩(wěn)定性影響因素均值和極差

3 地下水位上升對邊坡穩(wěn)定性的影響分析

3.1 地下水位上升對安全系數(shù)的影響分析

定義地下水位變化對邊坡穩(wěn)定性的影響系數(shù)ξ為

由表4可知，隨著地下水位的升高邊坡的安全系數(shù)呈現(xiàn)減小的趨勢。當?shù)叵滤惠^低時，地下水的升高對邊坡安全系數(shù)的影響不大。主要是由于在低水位時期，邊坡內(nèi)部的地下水位尚未與邊坡內(nèi)最小安全系數(shù)的滑動面相交，因而對此邊坡而言整體安全系數(shù)仍未改變。當?shù)叵滤怀掷m(xù)上升接近最小安全系數(shù)滑動面時，此邊坡的安全系數(shù)開始有下降的現(xiàn)象。隨著地下水位繼續(xù)升高，邊坡的安全系數(shù)逐漸減小，直到水位達到20 m水深時地下水位變化對邊坡的影響系數(shù)明顯增大。由此可知，在水位達到某一臨界水位時，邊坡的安全系數(shù)開始急劇下降。因此在實際工程中，應重視地下水對工程穩(wěn)定性的削弱作用。

采用強度折減法計算邊坡的安全系數(shù)。為了研究地下水位升高對邊坡破壞模式的影響規(guī)律通過強度折減法，分別提取了不同水位下邊坡發(fā)生滑動時的剪應變云圖，如圖5所示。

表4 水位變化對邊坡安全系數(shù)的影響

圖5 不同水位下邊坡發(fā)生滑動時的應變值

由圖5可知，隨著水位的升高，邊坡發(fā)生剪切破壞的最大剪應變逐漸增加，水位14 m時，邊坡發(fā)生滑動時的最大剪應變達到76%，且最大剪應變出現(xiàn)在坡腳位置。此外，當?shù)叵滤惠^高時，坡腳處應力集中明顯，最大剪應變多分布在坡腳附近，邊坡坡面屈服區(qū)域接近貫通至坡腳，最危險滑動面沿著邊坡塑性區(qū)發(fā)生滑移。此外，隨著地下水位上升，最大剪應變分布區(qū)域和塑性區(qū)域向邊坡深處擴展。

圖6 不同水位下邊坡坡面的位移變化規(guī)律

圖7 不同水位下邊坡坡內(nèi)的位移變化規(guī)律

3.2 地下水位上升對邊坡位移的影響分析

分別提取了邊坡在不同水位下不同位置的水平位移和豎向位移，如圖6至圖11所示。

圖8 不同水位下邊坡坡底的位移變化規(guī)律

圖9 不同水位下邊坡坡中的位移變化規(guī)律

圖10 不同水位下邊坡坡頂?shù)奈灰谱兓?guī)律

圖11 不同水位下邊坡坡后的位移變化規(guī)律

由圖6可知，隨著水位的升高，沿著邊坡高度方向坡面的水平位移表現(xiàn)出先增加后減小的趨勢，豎向位移表現(xiàn)出增加趨勢，且在邊坡頂部達到最大值。而水平位移在邊坡某一高度處出現(xiàn)最大值，且隨著水位的升高邊坡的最大位移出現(xiàn)位置相比無水時偏低，說明地下水的升高使得邊坡的變形最大值由邊坡的最高處開始向坡趾處轉(zhuǎn)移。因此，在地下水位較高時坡趾是薄弱位置。當?shù)叵滤贿_到20 m時沿邊坡坡面的最大水平位移增加幅度明顯提高，說明當?shù)叵滤烈欢ǜ叨葧r，邊坡的最大水平位移會出現(xiàn)一個突變，此時隨著水位的繼續(xù)增加邊坡將會出現(xiàn)滑坡。

由圖7可知，隨著水位的升高，邊坡最大水平位移出現(xiàn)在地面以上(＞12 m)，即邊坡坡面臨空時位移較大。豎向位移同樣是表現(xiàn)出增加的趨勢，在邊坡頂部達到最大值。

由圖8可以看出，邊坡坡底的水平位移隨著X方向坐標的增加，水平位移同樣表現(xiàn)出先增加后減小的趨勢。在X方向坐標增加到某一位置時出現(xiàn)一個突變，說明此位置接近邊坡滑移面位置。且隨著水位的升高，位移最大值出現(xiàn)位置偏向邊坡內(nèi)部，說明地下水的影響使得邊坡滑動面向邊坡內(nèi)部延伸，其滑動規(guī)模增大。也可從邊坡的位移云圖得到進一步的驗證，此可作為邊坡滑坡預測的一個重要指標。而豎向位移由坡內(nèi)向坡外增加趨勢不明顯，且在邊坡坡面位置增加到最大值，增加幅度更為明顯。

由圖9可知，隨著Y方向坐標的增加，邊坡的水平位移表現(xiàn)出增加的趨勢，且隨著水位的增加邊坡的位移也呈現(xiàn)增加的趨勢。而豎向位移呈現(xiàn)出先增大后減小的趨勢。

由圖10可知，邊坡坡頂?shù)乃轿灰坪拓Q向位移隨著X方向坐標增大，水平位移呈現(xiàn)增加的趨勢，在坡面達到最大值;而豎向位移表現(xiàn)出減小的趨勢。此可由邊坡的浸潤線位置得到解釋，由于靠近邊坡內(nèi)部地下水位較高，其孔隙水壓力大，使得邊坡有效應力降低，從而產(chǎn)生了較大的沉降。

由圖11可知，隨著X方向坐標增大，水平位移呈現(xiàn)減小的趨勢，而豎向位移在無水時表現(xiàn)出增加的趨勢，當存在地下水時其豎向位移在邊坡坡趾前方某一位置出現(xiàn)最大值，尤其是水位的增加其突變更為明顯。之后隨著X方向坐標的繼續(xù)增大，呈現(xiàn)減小的趨勢。因此，在坡趾前修建鐵路時，尤其是在地下水位較高時應該引起注意，防止路基沉降過大對高速鐵路運營帶來災難性的事故。

無水時邊坡最大變形位于邊坡坡面中上部(邊坡的最大水平位移主要出現(xiàn)在坡面，且最大水平位移可達到3.0 cm)。隨著水位的升高邊坡的水平位移表現(xiàn)出增加的趨勢，當水位達到22 m時，邊坡的最大水平位移達到了4.75 cm，比無水時的最大水平位移增加了58.33%。由此可知，地下水對邊坡變形的影響顯著，地下水位的升高應引起注意。從邊坡的最大位移出現(xiàn)位置可以看出，隨著地下水位的升高，邊坡的最大水平位移開始出現(xiàn)在邊坡中上部，且表現(xiàn)出逐步向邊坡內(nèi)部擴展的趨勢。說明隨著水位的升高，滑動區(qū)幾乎整個處于滲流面之下。在地下水滲透力作用下，坡趾部位的土體有效應力下降使得坡趾位置的位移開始增大，直至邊坡發(fā)生滑坡破壞。

4 結(jié)論

本文基于擬靜力法研究了地下水位變化對邊坡穩(wěn)定性和變形機制的影響規(guī)律，得到的主要結(jié)論為:

1)基于正交試驗方法進行了邊坡穩(wěn)定性的主要影響因素(黏聚力、內(nèi)摩擦角、坡角、地震峰值加速度和地下水位高度)的敏感性分析，其中內(nèi)摩擦角對邊坡穩(wěn)定性影響最大，其次為設計峰值加速度、坡角、黏聚力、地下水位高度。盡管相對其它4個因素地下水對邊坡安全系數(shù)影響較小，但是其它4個因素，都是邊坡穩(wěn)定性影響的敏感因素，且在進行計算時并沒有考慮地震作用下的動孔隙水壓力影響。通過影響因素的顯著性分析，地下水位的變化可使邊坡安全系數(shù)降低15%，因此地下水的影響不容忽略。

2)隨著水位的升高邊坡的水平位移表現(xiàn)出增加的趨勢，當水位達到22 m時，最大水平位移比無水時增加了58.33%。由此可知，地下水對邊坡變形的影響顯著，水位的升高應引起注意。從邊坡的最大位移出現(xiàn)位置可以看出，隨著地下水位的升高，邊坡的最大水平位移開始表現(xiàn)出邊坡后緣下挫并逐步向前緣擴展的趨勢。說明隨著水位升高至出現(xiàn)滲流面的情形時，滑動區(qū)幾乎整個處于滲流面之下。在地下水滲透力作用下，坡趾部位的土體有效應力下降使得坡趾位置的位移開始增大，直至邊坡發(fā)生滑坡破壞。

3)隨著水位的升高，邊坡最大位移出現(xiàn)位置相比無水時偏低，說明地下水的升高使得邊坡的變形最大值由邊坡的最高處開始向坡趾處轉(zhuǎn)移。因此，在地下水位較高時坡趾是薄弱位置。當?shù)叵滤簧烈欢ǜ叨葧r，邊坡的最大水平位移會出現(xiàn)一個突變，此時隨著水位的繼續(xù)增加邊坡將會出現(xiàn)滑坡。且隨著水位的升高，位移最大值出現(xiàn)位置偏向邊坡內(nèi)部，說明地下水的影響使得滑動面向邊坡內(nèi)部延伸，其滑動規(guī)模增大。

［1］HUANG M S，JIA C Q.Strength Reduction FEM in Stability Analysis of Soil Slopes Subjected to Transient Unsaturated Seepage［J］.Computers and Geotechnics，2009，36(1/2):93-101.

［2］TONG Fuguo，TIAN Bin，LIU Defu.A Coupling Analysis ofSlope Runoff and Infiltration under Rainfall［J］.Rock and Soil Mechanics，2008，29(4):1035-1040.

［3］RAHARDJO H，LEE T T，LEONG E C，et al.Responses of a Residual Soil Slope to Rainfall［J］.Canadian Geotechnical Journal，2005，42(3):340-351.

［4］GRIFFITHS D V，LANE P A.Slopes Stability Analysis by Finite Elements［J］.Geotechnique，1999，49(3):387-403.

［5］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.

［6］賈官偉，詹良通，陳云敏.水位驟降對邊坡穩(wěn)定性影響的模型試驗研究［J］.巖石力學與工程學報，2009，28(9):1798-1803.

［7］趙煉恒，羅強，李亮，等.地下水位變化對邊坡穩(wěn)定性影響的上限分析［J］.公路交通科技，2010，27(7):1-4.

［8］黃帥，宋波，牛立超，等.地震作用下動孔隙水壓力對邊坡永久位移影響的簡便計算方法［J］.建筑結(jié)構(gòu)學報，2014，35(3:):215-221.

［9］宋波，黃帥，蔡德鉤，等.地下水位變化對砂土邊坡地震動力響應的影響研究［J］.巖石力學與工程學報，2014，33(增1):2698-2706.

［10］HUANG Shuai，YAN Hongye，CAI Degou，et al.Dynamic Response Analysis of Soil Slope under Strong Earthquake［J］. Advance Material Research，2014，919(4):921-924.

［11］黃帥，宋波，蔡德鉤，等.近遠場地震下高陡邊坡的動力響應及永久位移分析［J］.巖土工程學報，2013，35(增2): 768-773.

［12］宋波，黃帥，蔡德鉤，等.地震和地下水耦合作用下砂土邊坡穩(wěn)定性研究［J］.巖土工程學報，2013，35(增2):862-868.

Quasi-static analysis on influence of rising of groundwater level on side-slope stability

CAI Degou1，2，HUANG Shuai3，YAN Hongye1，2，CHEN Feng1，2，ZHANG Qianli1，2，YAO Jianping1，2
(1.State Key Laboratory for Track Technology of High-speed Railway，Beijing 100081，China; 2.Railway Engineering Research Institute，China Academy of Railway Sciences，Beijing 100081，China; 3.Institute of Crustal Dynamics，China Earthquake Administration，Beijing 100085，China)

With the help of finite element model，different factors(slope angle，adhesion force，angle of internal friction，underground water level，and designed peak acceleration)and the possible impact they might generate to slope stability were studied，with special focus being given to relation among underground water level，slope stability and deformation.The results show that angle of internal friction has proven to be the most influential factor in this study，which is followed by designed peak acceleration，slope angle，adhesion force and underground water level.It needs to be noted that the last factor cannot be overlooked，as its change may lead to a 15%drop in slope's safety factor.If the water level reaches a relatively higher level(H=22 m for instance)，the maximum horizontal displacement is increased by 58.33%，compared with the case where H stands at 0.As the water level climbs up，the maximum horizontal displacement starts to take place at the forward edge，rather than the trailing edge of the slope; in other words，the maximum value，which is first observed at the peak of the slope，now moves to the toe.A sudden change is noticed in the slope's maximum horizontal displacement，as the underground water level reaches certain point.The region where the maximum values first are found now expand to the interior part of the slope，which indicates the expansion of the sliding scale.Therefore，it should be considered as an important index in the analysis of slope stability.

Sand slope;Underground water level;Quasi-static force;Safety factor;Displacement

U213.1+3

A

10.3969/j.issn.1003-1995.2015.01.13

1003-1995(2015)01-0056-07

(責任審編 趙其文)

2014-09-15;

2014-12-03

國家高技術研究發(fā)展計劃專項經(jīng)費資助(2011AA11A102);鐵道科學技術研究發(fā)展中心科研項目(J2014G006);高速度鐵路軌道技術國家重點實驗室開放課題基金資助(2012SKL01)

蔡德鉤(1978—)，男，浙江文成人，副研究員，博士。