楊宇軒，霍志濤，方仕達(dá)，柯 超，楊超越，汪 洋

(1.中國(guó)地質(zhì)大學(xué)(武漢)工程學(xué)院，湖北 武漢 430074；2.中南地質(zhì)科技創(chuàng)新中心，湖北 武漢 430205；3.中國(guó)地質(zhì)大學(xué)(武漢)地質(zhì)調(diào)查研究院，湖北 武漢 430074)

滑帶土的力學(xué)強(qiáng)度參數(shù)對(duì)滑坡穩(wěn)定性的影響很大。葛云峰等分別采用正交數(shù)值模擬試驗(yàn)和動(dòng)態(tài)反演分析，研究了滑坡滑動(dòng)面力學(xué)參數(shù)對(duì)滑坡變形特征的影響；Tang等對(duì)三峽庫(kù)區(qū)滑坡的活動(dòng)規(guī)律和滑動(dòng)面滲透性的相關(guān)性進(jìn)行了研究；何晨等對(duì)三峽庫(kù)區(qū)黃土坡滑坡滑帶土在不同pH值條件下的溶出特征和釋放規(guī)律進(jìn)行了研究；何坤等基于現(xiàn)場(chǎng)調(diào)查、無(wú)人機(jī)航測(cè)、數(shù)值模擬分析等對(duì)四川省玄武巖古滑坡復(fù)活機(jī)制進(jìn)行了分析，結(jié)果表明由于滑帶長(zhǎng)期受地下水影響，滑帶土體強(qiáng)度逐漸削弱，抗剪強(qiáng)度降低，從而誘發(fā)了古滑坡復(fù)活。

三峽庫(kù)區(qū)水位的年復(fù)一年的波動(dòng)會(huì)導(dǎo)致部分古滑坡的復(fù)活或引起滑坡間歇性活動(dòng)。對(duì)于間歇活動(dòng)的復(fù)活型古滑坡，一般用滑帶土的殘余強(qiáng)度來(lái)分析滑坡的穩(wěn)定性。環(huán)剪儀因具有剪切面積固定、允許試樣沿一個(gè)方向連續(xù)剪切的特點(diǎn)，常常用于測(cè)試滑坡滑帶土的殘余強(qiáng)度。許多學(xué)者利用環(huán)剪儀研究了滑坡滑帶土的含水率對(duì)其殘余強(qiáng)度的影響，如Wei等研究了不同風(fēng)化層花崗巖滑坡滑帶土殘余強(qiáng)度隨含水率的變化特征，結(jié)果表明隨著滑帶土含水率的降低，黏土有顯著的軟化趨勢(shì)，但砂土有硬化現(xiàn)象，且在低含水率狀態(tài)下，含水率變化對(duì)滑帶土內(nèi)聚力和內(nèi)摩擦強(qiáng)度有顯著的影響；古鵬翔等通過(guò)室內(nèi)環(huán)剪試驗(yàn)研究了邊坡滑帶土在不同含水率下的蠕變特性，建立了考慮含水率變化的滑帶土蠕變本構(gòu)模型；吳迪等對(duì)殘積土進(jìn)行了環(huán)剪試驗(yàn)，得到了滑帶土含水率與其殘余強(qiáng)度指標(biāo)之間的關(guān)系，并證明大變形下殘積土具有浸水軟化的特性；Zhou等、吳岱諠等、翟棟梁等和Lian等分別對(duì)根-土復(fù)合材料、金坪子滑坡滑帶土、典型黃土和灞橋滑坡滑帶土進(jìn)行了不同含水率下的環(huán)剪試驗(yàn)，結(jié)果都表明滑帶土的殘余強(qiáng)度受其含水率的影響較大，但滑帶土的殘余強(qiáng)度隨含水率的變化規(guī)律不同。受庫(kù)水位波動(dòng)的影響，三峽庫(kù)區(qū)滑坡不同位置處滑帶土的含水率不同，滑坡不同位置處滑帶土的殘余強(qiáng)度也有差異，因此開展不同含水率下滑帶土殘余強(qiáng)度的研究對(duì)于分析三峽庫(kù)區(qū)滑坡穩(wěn)定性具有重要的意義。

本文基于環(huán)剪試驗(yàn)，開展了不同含水率和不同法向應(yīng)力作用下四方碑滑坡滑帶土的殘余強(qiáng)度特性分析，并結(jié)合滑坡滲流模擬的地下水水位對(duì)滑帶土殘余強(qiáng)度進(jìn)行分區(qū)取值，進(jìn)而開展了滑坡穩(wěn)定性計(jì)算。該方法考慮了滑帶土力學(xué)指標(biāo)在空間上的差異性，對(duì)于復(fù)活型滑坡滑帶土強(qiáng)度參數(shù)的選取以及提高滑坡穩(wěn)定性評(píng)價(jià)的精度具有重要的意義。

1 四方碑滑坡基本概況和變形特征

圖1 四方碑滑坡工程地質(zhì)平面圖Fig.1 Geological map of Sifangbei landslide

圖2 四方碑滑坡工程地質(zhì)剖面圖Fig.2 Engineering geological profile of Sifangbei Landslide

根據(jù)地表位移監(jiān)測(cè)點(diǎn)監(jiān)測(cè)結(jié)果，2019年1月至4月滑坡累計(jì)位移曲線表現(xiàn)為近水平狀，呈現(xiàn)基本穩(wěn)定狀態(tài)；2019年4月至9月滑坡累計(jì)位移曲線隨時(shí)間呈階躍型上升，9月至11月滑坡累計(jì)位移曲線又變?yōu)榛痉€(wěn)定狀態(tài)，見圖3。四方碑滑坡滑帶土在周期性庫(kù)水位變動(dòng)的作用下，含水率發(fā)生變化，進(jìn)而引起滑帶土抗剪強(qiáng)度發(fā)生變化，滑坡常呈現(xiàn)間歇性活動(dòng)特點(diǎn)，據(jù)此將滑坡活動(dòng)分為活動(dòng)期和間歇期?；麻g歇性活動(dòng)演化不僅受庫(kù)水位變化引起的孔隙水壓力變化的影響，而且也與滑帶土抗剪強(qiáng)度的變化關(guān)系密切。

圖3 四方碑滑坡累計(jì)位移曲線Fig.3 Cumulative displacement curves of Sifangbei landslide

2 四方碑滑坡滑帶土殘余強(qiáng)度特性分析

2. 1 滑帶土基本性質(zhì)

滑帶土試樣取自重慶市萬(wàn)州區(qū)四方碑滑坡鉆孔，測(cè)得滑坡后緣和中部的滑帶土質(zhì)量含水率分別為16.0%和20.5%，滑坡前緣的滑帶土由于位于地下水水位以下，故取其飽和含水率為24.0%，3種滑帶土試樣含水率對(duì)應(yīng)的飽和度分別為0.67、0.85和1.00。按照《土工試驗(yàn)方法標(biāo)準(zhǔn)》(GB/T 50123—2019)要求，對(duì)滑帶土試樣進(jìn)行了基本物理力學(xué)性質(zhì)指標(biāo)的測(cè)試，其測(cè)試結(jié)果見表1。

表1 四方碑滑坡滑帶土基本物理力學(xué)性質(zhì)指標(biāo)

2. 2 滑帶土環(huán)剪試驗(yàn)方案

環(huán)剪試驗(yàn)采用德國(guó)Wille Geotechnic公司生產(chǎn)的ARS型全自動(dòng)閉合回路控制環(huán)剪儀，剪切盒內(nèi)環(huán)直徑為100 mm，外環(huán)直徑為150 mm，有效高度為30 mm，其剪切面積為98.13 cm。本次試驗(yàn)利用過(guò)2 mm篩的重塑土樣進(jìn)行環(huán)剪試驗(yàn)，控制重塑土樣干密度為1.77 g/cm，即與原狀土樣干密度相同；配置含水率為16.0%、20.5%、24.0%的滑帶土試樣，分別對(duì)應(yīng)滑坡后緣、中部和前緣滑帶土的含水率。四方碑滑坡滑體厚度區(qū)間較大，厚度為6～38 m，根據(jù)滑帶土上覆自重應(yīng)力綜合確定試驗(yàn)的法向應(yīng)力分別為200 kPa、325 kPa、450 kPa、575 kPa和700 kPa，并根據(jù)滑坡變形特征確定試驗(yàn)的剪切速率為0.06 mm/min。試驗(yàn)采用單級(jí)剪切程序，試驗(yàn)前先讓滑帶土試樣在各法向應(yīng)力下固結(jié)穩(wěn)定，然后采取固定的剪切速率對(duì)滑帶土試樣進(jìn)行排水剪切，直至滑帶土試樣達(dá)到殘余狀態(tài)為止。

2. 3 滑帶土殘余強(qiáng)度特性

不同含水率和不同法向應(yīng)力作用下四方碑滑坡滑帶土的剪應(yīng)力-剪切位移曲線，見圖4。

由圖4可見，四方碑滑坡滑帶土的殘余強(qiáng)度隨著法向應(yīng)力的增加而增大，并在相同法向應(yīng)力作用下，隨著滑帶土試樣含水率的增加，其殘余強(qiáng)度均減小。

不同含水率下四方碑滑坡滑帶土的殘余強(qiáng)度參數(shù)，見表2。

表2 不同含水率下四方碑滑坡滑帶土的殘余強(qiáng)度參數(shù)

圖4 不同含水率和不同法向應(yīng)力作用下四方碑滑坡滑帶土的剪應(yīng)力-剪切位移曲線Fig.4 Shear stress-shear displacement curves of sliding zone soil of Sifangbei landslide under different moisture content and different normal stress

由表2可知，與含水率為20.5%的滑帶土相比，含水率為16.0%的滑帶土的殘余黏聚力減小了33.51%，殘余內(nèi)摩擦角減小了38.33%；與含水率為24.0%的滑帶土相比，含水率為16.0%的滑帶土的殘余黏聚力減小了37.86%，殘余內(nèi)摩擦角減小了53.92%，說(shuō)明隨著含水率的增加，滑帶土的殘余黏聚力和殘余內(nèi)摩擦角都減小，且滑帶土的殘余內(nèi)摩擦角衰減程度更大。

通過(guò)數(shù)值擬合發(fā)現(xiàn)：四方碑滑坡滑帶土的殘余強(qiáng)度與法向應(yīng)力呈現(xiàn)較好的線性關(guān)系，符合Mohr-Coulomb準(zhǔn)則(見圖5)；該滑坡滑帶土的殘余內(nèi)摩擦角與含水率之間呈現(xiàn)良好的負(fù)指數(shù)關(guān)系(見圖6)，這與王煒等得出的結(jié)論類似。

圖5 四方碑滑坡滑帶土殘余強(qiáng)度與法向應(yīng)力的擬合 關(guān)系曲線Fig.5 Fitting curves of residual strength and normal stress of sliding zone soil of Sifangbei landslide

圖6 四方碑滑坡滑帶土殘余內(nèi)摩擦角與含水率的 擬合關(guān)系曲線Fig.6 Fitting curve between residual internal friction angle and moisture content of sliding zone soil of Sifangbei landslide

3 滑帶土殘余強(qiáng)度空間差異性和滑坡穩(wěn)定性分析

3. 1 滑坡滲流場(chǎng)分析

本文選取四方碑滑坡1-1′剖面為研究對(duì)象，采用GeoStudio軟件對(duì)該滑坡滲流特征和穩(wěn)定性進(jìn)行了數(shù)值模擬計(jì)算與分析，建立了滑坡滲流特征和穩(wěn)定性計(jì)算模型，見圖7。該滑坡滑動(dòng)面為不透水邊界，滑坡表面以水庫(kù)水位為分界線，其上為零流量邊界，其下為變水頭邊界，其水頭隨著庫(kù)水位的變化而變化，依據(jù)現(xiàn)有水文孔監(jiān)測(cè)數(shù)據(jù)設(shè)定初始地下水水位線。由單環(huán)滲水試驗(yàn)測(cè)得的滑帶土滲透系數(shù)為0.3 m/d，選取2017年10月28日到2018年10月28日的庫(kù)水位變化作為滲流模擬的分析工況。

圖7 四方碑滑坡滲流特征和穩(wěn)定性計(jì)算模型圖Fig.7 Seepage characteristics and stability calculation model of Sifangbei landslide

3.2 滑帶土殘余強(qiáng)度分區(qū)取值和滑坡穩(wěn)定性分析

由前述室內(nèi)環(huán)剪試驗(yàn)可知，四方碑滑坡滑帶土殘余強(qiáng)度受含水率的影響較大。受庫(kù)水位波動(dòng)的影響，庫(kù)區(qū)滑坡不同位置處滑帶土的含水率不同，基于滑坡滲流模擬的地下水水位對(duì)滑帶土進(jìn)行分區(qū)，計(jì)算滑坡穩(wěn)定性時(shí)不同區(qū)域的滑帶土取不同的計(jì)算參數(shù)。地下水水位線以上的滑坡后緣滑帶土含水率較小，計(jì)算滑坡穩(wěn)定性時(shí)滑帶土抗剪強(qiáng)度取參數(shù)1(見表2)；庫(kù)水位以下的滑坡前緣滑帶土按飽和狀態(tài)考慮，計(jì)算滑坡穩(wěn)定性時(shí)滑帶土抗剪強(qiáng)度取參數(shù)3(見表2)；處于地下水水位和庫(kù)水位之間的滑坡中部滑帶土，含水率介于滑坡前緣滑帶土和滑坡后緣滑帶土之間，計(jì)算滑坡穩(wěn)定性時(shí)滑帶土抗剪強(qiáng)度取參數(shù)2(見表2)。參考四方碑滑坡勘察報(bào)告進(jìn)行其他材料參數(shù)取值，滑體黏聚力取25 kPa，滑體內(nèi)摩擦角取13°，滑體重度取21 kN/m，滑帶重度取20.3 kN/m。

在滑坡滲流模擬的基礎(chǔ)上，采用極限平衡法中的摩根斯頓-普賴斯法(簡(jiǎn)稱MP法)計(jì)算175 m和145 m庫(kù)水位以下四方碑滑坡的穩(wěn)定性，水位線采用滑坡滲流模擬得出的地下水水位線，滑帶土不同區(qū)段抗剪強(qiáng)度參數(shù)取不同值。通過(guò)計(jì)算得到的175 m和145 m庫(kù)水位以下滑坡的穩(wěn)定性系數(shù)分別為1.14和1.17，根據(jù)滑坡穩(wěn)定性狀態(tài)劃分標(biāo)準(zhǔn)(見表3)可知，該滑坡處于基本穩(wěn)定和穩(wěn)定狀態(tài)，與通過(guò)野外調(diào)查得到的滑坡穩(wěn)定性現(xiàn)狀相符。若不對(duì)滑帶土進(jìn)行分區(qū)，計(jì)算滑坡穩(wěn)定性時(shí)全部采用滑帶土的飽和殘余強(qiáng)度參數(shù)，計(jì)算得到的175 m和145 m庫(kù)水位以下滑坡的穩(wěn)定性系數(shù)分別為1.04和1.08，根據(jù)表3判斷該滑坡處于欠穩(wěn)定和基本穩(wěn)定狀態(tài)，這低估了滑坡的穩(wěn)定性現(xiàn)狀。

表3 滑坡穩(wěn)定性狀態(tài)劃分標(biāo)準(zhǔn)

4 結(jié) 論

本文以四方碑滑坡滑帶土為例，通過(guò)環(huán)剪試驗(yàn)分析了不同含水率和不同法向應(yīng)力作用下滑帶土的殘余強(qiáng)度特性，并基于滑坡滲流模擬結(jié)果對(duì)滑帶土殘余強(qiáng)度的空間差異性進(jìn)行分區(qū)取值來(lái)分析該滑坡的穩(wěn)定性，主要得到如下結(jié)論：

(1) 不同含水率滑帶土試樣，其殘余強(qiáng)度與法向應(yīng)力呈現(xiàn)較好的線性關(guān)系，符合Mohr-Coulomb準(zhǔn)則；法向應(yīng)力相同時(shí)，隨著滑帶土含水率的增加，滑帶土試樣的殘余強(qiáng)度減小。隨著滑帶土含水率的增加，滑帶土試樣的殘余黏聚力和殘余內(nèi)摩擦角都減小，且滑帶土殘余內(nèi)摩擦角的衰減程度更大。通過(guò)數(shù)值擬合發(fā)現(xiàn)，滑帶土殘余內(nèi)摩擦角與含水率之間呈現(xiàn)良好的負(fù)指數(shù)關(guān)系。

(2) 基于滑坡滲流模擬的地下水水位對(duì)四方碑滑坡滑帶土進(jìn)行分區(qū)，計(jì)算滑坡穩(wěn)定性時(shí)不同區(qū)域的滑帶土取不同的計(jì)算參數(shù)，結(jié)果表明：175 m和145 m庫(kù)水位以下該滑坡處于基本穩(wěn)定和穩(wěn)定狀態(tài)，滑坡穩(wěn)定性計(jì)算結(jié)果與通過(guò)野外調(diào)查得到的滑坡穩(wěn)定性現(xiàn)狀相符，表明本文方法具有可靠性。