沈秋池，蘇 燕

（1．福建船政交通職業(yè)學(xué)院道路工程系，福州350007；2．福州大學(xué)土木工程學(xué)院，福州350116）

庫(kù)水和降雨耦合作用下公路滑坡研究

沈秋池1，蘇 燕2

（1．福建船政交通職業(yè)學(xué)院道路工程系，福州350007；2．福州大學(xué)土木工程學(xué)院，福州350116）

以福建山區(qū)庫(kù)岸公路滑坡為對(duì)象，采用室內(nèi)模型試驗(yàn)和有限元數(shù)值模擬方法，研究庫(kù)岸公路降雨滑坡形成機(jī)理，探討庫(kù)水位下降與降雨聯(lián)合作用下導(dǎo)致臨水邊坡發(fā)生失穩(wěn)的模式和原因。研究表明庫(kù)水位變動(dòng)與降雨的聯(lián)合作用將改變沿坡體內(nèi)的滲流場(chǎng)，降低滑坡土體的力學(xué)強(qiáng)度，使得滑坡前緣由于浮力減重和浸泡軟化作用導(dǎo)致抗滑力減小，滑坡中后部下滑力增大，產(chǎn)生滑坡現(xiàn)象。

水位降落；降雨；滑坡；模型試驗(yàn)；有限元

邊坡失穩(wěn)是目前發(fā)生頻率較高的地質(zhì)災(zāi)害之一，水庫(kù)庫(kù)區(qū)時(shí)常有庫(kù)水位漲落和降雨發(fā)生往往會(huì)導(dǎo)致庫(kù)岸公路滑坡失穩(wěn)，從而造成的人員傷亡和經(jīng)濟(jì)損失巨大［1］。如何科學(xué)地確定在降雨和庫(kù)水漲落作用下的邊坡穩(wěn)定性問(wèn)題，是工程中亟需解決的問(wèn)題。

目前，國(guó)內(nèi)外學(xué)者在庫(kù)岸公路降雨邊坡的滲流理論、數(shù)值模擬分析、模型試驗(yàn)、危險(xiǎn)性評(píng)價(jià)等方面取得了一定的研究成果，但還有許多需要深入研究和改進(jìn)的地方［2－4］。本文擬通過(guò)模型試驗(yàn)和數(shù)值分析結(jié)合的方法，模擬臨水邊坡在水位降落和降雨聯(lián)合作用下的失穩(wěn)情況，通過(guò)數(shù)碼攝像、動(dòng)態(tài)應(yīng)變采集儀來(lái)記錄水位下降和降雨聯(lián)合作用下邊坡產(chǎn)生裂縫的形成發(fā)展過(guò)程、滑動(dòng)面的形態(tài)、滑坡失穩(wěn)過(guò)程，揭示臨水邊坡在水位下降和降雨聯(lián)合作用下的失穩(wěn)方式和原因。

1 試驗(yàn)裝置及模型邊坡

福建山區(qū)滑坡體多為基巖表面的淺層覆蓋土體，本試驗(yàn)據(jù)此將滑坡地質(zhì)結(jié)構(gòu)設(shè)置為基巖型邊坡，在模型箱底部鋪設(shè)混凝土以模擬基巖。試驗(yàn)?zāi)Ｐ桶ɑ履Ｐ拖?、人工降雨系統(tǒng)、動(dòng)態(tài)應(yīng)變采集儀。模型箱尺寸為285×100×130cm3，坡度為30°，模型箱底座和骨架采用鋼板和型鋼，為便于觀察，四周采用有機(jī)玻璃（圖1）。模型箱底部設(shè)有排水閥門(mén)，通過(guò)控制閥門(mén)開(kāi)度，可以滿足庫(kù)水位下降速率的要求。人工降雨系統(tǒng)采用農(nóng)用噴頭雨滴發(fā)生器，距離試驗(yàn)坡面4m，能夠達(dá)到模擬天然降雨終極速度的要求［3］。降雨噴頭組采用直徑20mm的PVC管進(jìn)行串聯(lián)，由蓄水池通過(guò)水泵供水，并利用閥門(mén)和分流管道來(lái)調(diào)節(jié)降雨強(qiáng)度。動(dòng)態(tài)應(yīng)變采集儀包括孔隙水壓力計(jì)、土壤水分計(jì)，分別用于監(jiān)測(cè)庫(kù)水下降和降雨聯(lián)合作用下邊坡內(nèi)部孔隙水壓力、含水量的變化，其布置如圖2所示。

圖1 滑坡模型試驗(yàn)裝置

圖2 孔隙水壓力計(jì)和土壤水分計(jì)布置圖

邊坡模型試驗(yàn)采用福建山區(qū)殘積亞黏土，天然重度為18．8kN／m3，根據(jù)三軸試驗(yàn)測(cè)得其黏聚力為11．2kPa，內(nèi)摩擦角為19．6°，由變水頭滲透試驗(yàn)測(cè)得其滲透系數(shù)k＝6．7×10－6cm／s。邊坡成型采用擊實(shí)法，分層鋪土，每層厚度為5cm。每分層鋪好夯實(shí)后，用環(huán)刀取樣，測(cè)其干密度和含水量，保證每層土密實(shí)度一致。當(dāng)夯實(shí)至所需埋設(shè)土壤水分計(jì)和孔隙水壓力計(jì)的分層時(shí)，按照布置圖緩慢插入儀器，并整理固定好儀器線路，避免滑坡過(guò)程中儀器損毀。模型邊坡填筑完成后，靜置72h，使邊坡自然沉降成型。

2 試驗(yàn)過(guò)程及結(jié)果分析

2．1 水位和降雨變化過(guò)程

開(kāi)始試驗(yàn)前，緩慢勻速地向邊坡水庫(kù)內(nèi)注水，使庫(kù)水位逐漸上升至34cm。庫(kù)水上漲使坡面含水量升高，非飽和土基質(zhì)吸力降低；同時(shí)，涉水前緣土體孔隙水壓力增大，有效應(yīng)力降低，導(dǎo)致邊坡涉水前緣出現(xiàn)局部崩塌現(xiàn)象。坡體中部土體在庫(kù)水的浸泡作用下軟化，導(dǎo)致黏聚力和內(nèi)摩擦角減小，土體抗剪強(qiáng)度降低而產(chǎn)生局部變形。當(dāng)模型邊坡達(dá)到飽和穩(wěn)定后，開(kāi)啟排水閘門(mén)，使坡外庫(kù)水位按0．08cm／min勻速下降，歷時(shí)425min。隨著排水孔水頭的降低，水庫(kù)流量也不斷減少。在庫(kù)水位下降過(guò)程中降雨，分為2次，每次降雨歷時(shí)7min，分別從第47分鐘至54分鐘、147分鐘至154分鐘，降雨強(qiáng)度均為110mm／h。

2．2 滑坡漸進(jìn)破壞

當(dāng)試驗(yàn)進(jìn)行到第40分鐘時(shí)，邊坡距水面約50cm處出現(xiàn)平行于坡肩的張拉裂縫，裂縫長(zhǎng)度約50cm；距水位10cm處出現(xiàn)與水平方向呈45°的剪切裂縫，如圖3所示。隨著庫(kù)水位持續(xù)下降和降雨的聯(lián)合作用，裂縫的深度和寬度不斷增大。第一場(chǎng)人工降雨結(jié)束后，滑坡尚未啟動(dòng)。在庫(kù)水下降和水分蒸發(fā)的作用下，試驗(yàn)進(jìn)行到第70分鐘時(shí)滑坡啟動(dòng)，剪切裂縫和張拉裂縫下方土體向下運(yùn)動(dòng)，運(yùn)動(dòng)范圍為1～3cm，滑坡體在坡腳處堆積，支撐住上部土體，阻擋滑坡的進(jìn)一步發(fā)展，坡體暫時(shí)達(dá)到穩(wěn)定。當(dāng)試驗(yàn)進(jìn)行到第152分鐘，即第2場(chǎng)雨第2分鐘時(shí)，由于孔隙水壓力驟然增大導(dǎo)致土體抗剪強(qiáng)度急劇降低，坡體再次滑動(dòng)，土體沿張拉裂縫產(chǎn)生向下塊體滑動(dòng)。此時(shí)滑動(dòng)土體上方的土體失去支撐，產(chǎn)生滑坡臨空面，在自身重力和水體滲透力的作用下，在斜坡上部再次出現(xiàn)張拉裂縫，促使土體再次發(fā)生塊體滑動(dòng)，如此重復(fù)破壞產(chǎn)生漸進(jìn)向后滑動(dòng)破壞。邊坡破壞的最終形態(tài)如圖4所示。

圖3 坡面上的張拉裂縫和剪切裂縫

2．3 含水量和孔隙水壓力變化

圖4是邊坡中部CH2土壤水分計(jì)隨時(shí)間變化的過(guò)程圖，庫(kù)水剛開(kāi)始下降時(shí)，含水量基本不變，分析原因是坡內(nèi)水位降落速率滯后于庫(kù)水位的下降和蒸發(fā)作用。經(jīng)歷第一場(chǎng)降雨（7min），土體含水量迅速增加，向飽和狀態(tài)靠攏，降雨結(jié)束后，由于水向下滲入和蒸發(fā)作用，含水量緩慢下降。當(dāng)遭遇第二場(chǎng)降雨（7min）時(shí)，在第152分鐘，土體沿張拉裂縫產(chǎn)生向下塊體滑動(dòng)，體積含水量急劇下降，這時(shí)由于邊坡失穩(wěn)后土壤水分計(jì)露出土體（如圖5所示），無(wú)法繼續(xù)測(cè)量。

圖4 坡外水位CH2含水量隨時(shí)間變化過(guò)程

圖5 滑坡最終形態(tài)

圖6是邊坡中部CH5土壤孔隙水壓力隨時(shí)間變化過(guò)程圖，在庫(kù)水下降初期，孔隙水壓力隨著庫(kù)水的下降而減小，在第一場(chǎng)降雨3min后，孔隙水壓力由0．3kPa突增到0．8kPa，而邊坡出現(xiàn)裂縫破壞，隨著降雨和庫(kù)水的作用，裂縫越來(lái)越大，這驗(yàn)證了孔隙水壓力的增大會(huì)使土體有效應(yīng)力減小，進(jìn)而坡體產(chǎn)生破壞。在第二場(chǎng)降雨中，孔隙水壓力從－0．3kPa突增到0．7kPa，此時(shí)坡體從裂縫處開(kāi)始發(fā)生漸進(jìn)向后的滑動(dòng)破壞。試驗(yàn)過(guò)程中，孔隙水壓力實(shí)測(cè)結(jié)果和試驗(yàn)現(xiàn)象都顯示出土體表面逸出點(diǎn)有明顯的坡內(nèi)指向坡外的滲流，這是導(dǎo)致發(fā)生滑坡的重要原因。并且水位降落過(guò)程中，坡內(nèi)水位降落速度滯后于坡前水位降落速度，這會(huì)引起坡體內(nèi)水向外排除，這種滲流力和庫(kù)水降落差造成的壓力卸載對(duì)庫(kù)水邊坡的穩(wěn)定性非常不利。

圖6 坡外水位CH5孔隙水壓力隨時(shí)間變化過(guò)程

3 數(shù)值模擬

3．1 有限元模型的建立

基于物理模型試驗(yàn)的尺寸，建立了Seep／w滲流分析的有限元模型，數(shù)值模型尺寸和有限元?jiǎng)澐值木W(wǎng)格如圖7所示。有限元單元采用結(jié)構(gòu)化四邊形單元，單元尺寸為0．02m×0．02m，共1 675個(gè)節(jié)點(diǎn)，1 576個(gè)單元組成。

圖7 數(shù)值模型及網(wǎng)格劃分示意圖

在數(shù)值模擬中，參數(shù)選取的正確性直接決定了分析結(jié)果的合理與否，而這些參數(shù)主要是通過(guò)巖土體的物理力學(xué)實(shí)驗(yàn)、滲透實(shí)驗(yàn)、工程地質(zhì)類比以及經(jīng)驗(yàn)值等綜合確定的。本數(shù)值模擬的滲流和穩(wěn)定性計(jì)算過(guò)程中涉及到的滑坡巖土體的參數(shù)主要根據(jù)模型試驗(yàn)確定：天然重度取18．8kN／m3，黏聚力C＝11．2kPa，內(nèi)摩擦角φ＝19．6°，滲透系數(shù)K＝6．7× 10－6cm／s，飽和含水量為68%。

在滲流分析中，非飽和材料參數(shù)主要有土水特征曲線和滲透系數(shù)函數(shù)。采用上述模型參數(shù)，由Fredlund＆Xing模型擬合得到非飽和土滲透系數(shù)與基質(zhì)吸力的關(guān)系，如圖8所示。根據(jù)土體的滲透系數(shù)與體積含水量有關(guān)特性，并假定土體的殘余體積含水量為飽和含水量的10%，通過(guò)Van Genuchten模型對(duì)整個(gè)土體體積含水量方程進(jìn)行積分［5］，得到土的體積含水量與基質(zhì)吸力的關(guān)系，如圖9所示。根據(jù)圖8～9可知，土體的體積含水量越小，基質(zhì)吸力越大，致使?jié)B透系數(shù)越低，這是土體內(nèi)部的空氣氣泡阻礙水體在土體內(nèi)部流動(dòng)的結(jié)果。

圖8 滲透系數(shù)和基質(zhì)吸力的關(guān)系

圖9 含水量和基質(zhì)吸力的關(guān)系

在Seep／w模塊中，在一個(gè)節(jié)點(diǎn)上，只能給定水頭或流量邊界條件。Seep／w對(duì)流量邊界的假定是：若邊界流量小于土體滲透系數(shù)，采用流量邊界條件；若邊界流量大于土體滲透系數(shù)，坡面有徑流產(chǎn)生，則采用水頭邊界條件。本次數(shù)值分析中，根據(jù)降雨入滲強(qiáng)度為110mm／h和表層單元網(wǎng)格長(zhǎng)度，換算成邊界流量為6．11×10－5cm／s，大于土體滲透系數(shù)，因此采用水頭邊界條件。在庫(kù)水下降和降雨過(guò)程中，滲流表面和降雨表面的水頭邊界條件是動(dòng)態(tài)變化的，采用水位變化函數(shù)來(lái)表示，函數(shù)設(shè)置與模型試驗(yàn)一致，另外底部基巖設(shè)為不透水邊界。

3．2 數(shù)值模擬結(jié)果分析

圖10給出了在庫(kù)水下降速率為0．08cm／min，降雨強(qiáng)度110mm／h時(shí)的浸潤(rùn)線變化圖。由圖10可知，滑坡體浸潤(rùn)線的彎向由滑坡體內(nèi)向滑坡體外彎曲，說(shuō)明隨著庫(kù)水位的下降，岸坡地下水位也在逐漸下降，滑坡體內(nèi)地下水位的下降存在滯后現(xiàn)象。且隨著降雨時(shí)間的延續(xù)，滑坡土體內(nèi)非飽和區(qū)域變小，飽和區(qū)域變大，土壤含水率變大，因而浸潤(rùn)線不斷上升。隨著降雨的進(jìn)行，坡體淺層的基質(zhì)吸力逐漸減小，并隨著降雨的停止逐漸趨于穩(wěn)定值。基質(zhì)吸力的喪失會(huì)導(dǎo)致滑坡的非飽和抗剪強(qiáng)度的降低，進(jìn)而影響邊坡的穩(wěn)定性，發(fā)生局部滑坡。在降雨進(jìn)行的過(guò)程中和結(jié)束后的一段時(shí)間內(nèi)，雨水都在緩慢地入滲到滑坡體內(nèi)，淺層坡體基質(zhì)吸力逐漸喪失，并局部飽和，以至于出現(xiàn)滑坡現(xiàn)象。

圖10 浸潤(rùn)線變化圖

通過(guò)對(duì)庫(kù)水下降孔隙水壓力等值線云圖（圖11）分析發(fā)現(xiàn)：在庫(kù)水下降過(guò)程中疊加降雨的作用，邊坡前緣出現(xiàn)較大的孔隙水壓力，導(dǎo)致基質(zhì)吸力降低，引起坡體前緣變形增大。同時(shí)在淺層形成暫態(tài)飽和區(qū)和局部積水，增加坡體自重，加之邊坡上存在的張拉裂隙，降雨入滲至裂隙，使滑帶土含水量增加甚至局部飽和，導(dǎo)致抗剪強(qiáng)度降低，發(fā)生滑坡破壞。

將上述模擬結(jié)果導(dǎo)入slope／w模塊進(jìn)行邊坡穩(wěn)定性計(jì)算，得到：隨著水位降落與降雨的持續(xù)作用，邊坡的安全系數(shù)不斷降低，且在降雨初期變化速率較快，分析原因可能由于在庫(kù)水下降前期，坡前水壓力的卸載，致使坡內(nèi)和坡外形成的孔壓差對(duì)邊坡的不利影響大于坡內(nèi)摩擦力增加的有利影響，所以安全系數(shù)呈現(xiàn)驟降現(xiàn)象。

圖11 孔隙水壓力分布圖

4 結(jié)語(yǔ)

本文通過(guò)物理模型試驗(yàn)和有限元模擬了水位降落與降雨聯(lián)合作用下對(duì)臨水邊坡穩(wěn)定性的影響，得到的結(jié)論如下：

1）坡體在庫(kù)水和降雨的聯(lián)合作用下，由于坡內(nèi)的水位降落值滯后于庫(kù)水位下降值，坡體出現(xiàn)45°的剪切裂縫和平行于坡肩的張拉裂縫。隨著庫(kù)水位的下降和降雨作用，裂縫的深度和寬度越來(lái)越大，同時(shí)坡體前緣產(chǎn)生局部破壞，使其上方土體產(chǎn)生滑坡臨空面，進(jìn)而發(fā)生進(jìn)一步失穩(wěn)破壞現(xiàn)象。

2）由于降雨滲透使土體含水量增大，導(dǎo)致土體的容重增加和孔隙水壓力增大，基質(zhì)吸力降低，抗剪強(qiáng)度降低。滑坡前緣主要表現(xiàn)為由于浮力減重和浸泡軟化作用使抗滑力減小，滑坡中后部下滑力增大使邊坡安全系數(shù)降低，邊坡出現(xiàn)失穩(wěn)現(xiàn)象。

［1］卿篤干．五強(qiáng)溪庫(kù)岸滑坡穩(wěn)定性特征與信息系統(tǒng)研究［D］．長(zhǎng)沙：中南大學(xué)，2011．

［2］Borga M，Dalla Fontana G，Da Ros D，et al．Shallow landslide hazard assessment using aphysically based model and digital elevation data［J］．Environmental Geology，1998，2：81－88．

［3］羅先啟，劉德富，吳劍，等．雨水及庫(kù)水作用下滑坡穩(wěn)定模型試驗(yàn)研究［J］．巖石力學(xué)與工程學(xué)報(bào)，2005，24（14）：2476－2483．

［4］賈官偉，詹良通，陳云敏．水位驟降對(duì)邊坡穩(wěn)定性影響的模型試驗(yàn)研究［J］．巖石力學(xué)與工程學(xué)報(bào)，2009，28（9）：1798－1803．

［5］Genuchten MTh Van．A closed－form equation for predicting for the hydraulic conductivity of unsaturated soils［J］．Soil Science Society of America Journal，1980，44：892－898．

The Research on Highway Landslide Influenced by Reservoir Water and Rainfall

SHEN Qiu－chi，etc．
（Department of Road Enieerning，F(xiàn)ujian Chuanzheng Communication College，F(xiàn)uzhou350007，China）

Taking mountainous rainfall－induced landslide in Fujian province as the study subject，combining indoor artificial rainfall landslide experiment with finite element numerical simulation method，the mechanism of rainfall－induced landslide of highway bank slope has been studied in this article．The water slope failure mode and causes of stability have been discussed with the effect of reservoir water level fluctuation and the changes of rainfall intensity on highway bank slope．The research shows that the seepage field along the slope body has been changed，and the mechanical strength of soil has been reduced．The anti－slide force in front of slope will decrease with the sliding force increasing caused by buoyancy weight loss and soaking softening，and the sliding force in the rear of the slope will increase．This will lead to landslide．

reservoir water level；rainfall；landslide；model test；finite element method

TU46

A

1009－8984（2016）02－0010－04

10．3969／j．issn．1009－8984．2016．02．003

2016－04－08

福建省自然科學(xué)基金項(xiàng)目（2014J01186）

沈秋池（1973－），男（漢），福建詔安，副教授主要研究滑坡災(zāi)害。