劉天宇，楊建輝，楊 強(qiáng)，薛興橋，鄧 兵

(1.四川省冶金地質(zhì)勘查局水文工程大隊(duì)，四川郫縣 611730;2.中國(guó)地質(zhì)調(diào)查局水文地質(zhì)環(huán)境地質(zhì)調(diào)查中心，河北保定 071051)

根據(jù)防洪、抗旱、發(fā)電的需要，三峽庫(kù)區(qū)庫(kù)水位在145～175 m之間頻繁波動(dòng)，形成了變幅30 m的消落帶，在一定程度上改變了原有地質(zhì)環(huán)境的平衡狀態(tài)，加之庫(kù)區(qū)地理?xiàng)l件特殊、環(huán)境相對(duì)復(fù)雜、降雨量較大，加劇了滑坡等地質(zhì)災(zāi)害的發(fā)生，使三峽庫(kù)區(qū)成為地質(zhì)災(zāi)害重災(zāi)區(qū)［1－3］?；麦w水動(dòng)力學(xué)條件改變是導(dǎo)致水庫(kù)型滑坡發(fā)生的根本原因［4－6］，而滑坡體水動(dòng)力學(xué)條件改變與庫(kù)水位升降引起滑坡體滲流場(chǎng)變化密切相關(guān)，蓄水過程中引發(fā)的多起滑坡地質(zhì)災(zāi)害就證實(shí)了這一點(diǎn)。目前關(guān)于庫(kù)水位升降對(duì)滑坡穩(wěn)定性影響的研究相對(duì)較少，已有的成果多是采用定性描述及簡(jiǎn)單的極限平衡法進(jìn)行計(jì)算，沒有充分考慮地下水滲流場(chǎng)對(duì)滑坡穩(wěn)定的影響［7－9］，計(jì)算結(jié)果人為因素影響較大。因此，采集更加準(zhǔn)確的數(shù)據(jù)，建立相對(duì)細(xì)致的地質(zhì)模型，應(yīng)用長(zhǎng)時(shí)間的實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)作為模型校正依據(jù)，才能從整體上保證模擬結(jié)果的準(zhǔn)確性。

選擇三峽庫(kù)區(qū)李家坡滑坡作為研究對(duì)象，采用地學(xué)模擬軟件GeoStudio中的SEEP模塊和SLOPE模塊建立相應(yīng)模型，選用長(zhǎng)時(shí)間實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行校正，并用校正好的模型模擬庫(kù)水位以不同速率升降對(duì)滑坡穩(wěn)定性的影響。

1 李家坡滑坡概況

李家坡滑坡位于大寧河左岸，下部為龍水移民新村庫(kù)岸，滑坡體長(zhǎng)度約350 m，寬度約400 m，平均厚度約15 m，剪出口高程約150 m，體積約210萬(wàn)m3?；聟^(qū)基巖大面積出露，巴東組第一段(T12b)灰黃色粉—微晶泥質(zhì)灰?guī)r、泥灰?guī)r，內(nèi)陸湖相—瀉湖相碎屑巖和碳酸鹽巖沉積，總厚度50～100 m，強(qiáng)風(fēng)化層較厚;第四系殘坡積(Qel+dl)塊碎石、塊碎石黏土、巖屑碎石土夾碎塊石零星分布，上部結(jié)構(gòu)較松散—稍密，下部結(jié)構(gòu)中密，碎石、塊石成分主要為灰?guī)r、泥灰?guī)r、泥巖，一般呈次棱角狀及棱角狀，碎石大小為5～20 cm，塊石大小為0.3～1.5 m，與強(qiáng)風(fēng)化巖層呈過渡接觸?；禄瑒?dòng)帶為泥灰?guī)r強(qiáng)風(fēng)化和弱風(fēng)化接觸帶，與岸坡呈順層斜交，傾向345°，傾角35°～42°。岸坡坡角較陡，強(qiáng)風(fēng)化層較厚，在水庫(kù)蓄水上升過程中庫(kù)岸再造強(qiáng)烈。

2 庫(kù)水位升降對(duì)滑坡穩(wěn)定性影響模擬

2.1 計(jì)算流程

選用地學(xué)軟件GeoStudio中的SEEP、SLOPE模塊進(jìn)行滲流場(chǎng)模擬和滑坡穩(wěn)定性分析。根據(jù)勘探、試驗(yàn)、調(diào)查、監(jiān)測(cè)等資料，用SEEP模塊建立SEEP模型，用滑坡下部滑帶部位實(shí)時(shí)監(jiān)測(cè)的水分含量和孔隙水壓力數(shù)據(jù)校正該模型，監(jiān)測(cè)數(shù)據(jù)時(shí)長(zhǎng)為庫(kù)水位的1個(gè)升降周期，對(duì)不同工況進(jìn)行滲流場(chǎng)模擬;再將校正后的SEEP模型代入到SLOPE模塊中模擬庫(kù)水位以不同速度升降對(duì)滑坡穩(wěn)定的影響，分析不同工況下滑坡穩(wěn)定性的主要影響因素。SEEP模型基本參數(shù)見表1，SEEP模型校正模擬與實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)對(duì)比見圖1、2。

從圖1、2可以看出，模擬結(jié)果和實(shí)時(shí)監(jiān)測(cè)數(shù)據(jù)無論從總體趨勢(shì)還是各個(gè)時(shí)間點(diǎn)，總體吻合較好，說明建立的SEEP模型較接近真實(shí)狀況，模擬結(jié)果可靠，代入SLOPE模型后進(jìn)行穩(wěn)定性分析能較好地反映各個(gè)工況下滑坡穩(wěn)定性變化。

表1 SEEP模型基本參數(shù)

2.2 庫(kù)水位真實(shí)升降過程滑坡穩(wěn)定性計(jì)算

計(jì)算庫(kù)水位1個(gè)升降周期內(nèi)自然升降狀態(tài)下145、150、155、160、165、175 m 庫(kù)水位標(biāo)高對(duì)應(yīng)的滑坡穩(wěn)定性系數(shù)，計(jì)算結(jié)果見圖3，滑坡中下部?jī)蓚€(gè)深部位移監(jiān)測(cè)數(shù)據(jù)見圖4。

圖3 不同庫(kù)水位滑坡穩(wěn)定性系數(shù)

水庫(kù)型滑坡穩(wěn)定性主要受動(dòng)水壓力、靜水壓力、浮托力、巖土體物理參數(shù)(巖土體內(nèi)聚力c、內(nèi)摩擦角φ)改變共同影響。從圖3可以看出，在庫(kù)水位自然升降過程中滑坡穩(wěn)定性系數(shù)變化趨勢(shì)是先減小后增大，穩(wěn)定性的變化主要是看哪個(gè)影響因素占主導(dǎo)地位。庫(kù)水位真實(shí)升降過程中滑坡穩(wěn)定性系數(shù)均大于1.15，說明該滑坡體總體處于穩(wěn)定狀態(tài)。從圖4監(jiān)測(cè)數(shù)據(jù)可以看出，整個(gè)庫(kù)水位升降過程中滑坡下部滑帶部位(EL_LJP_1)最大位移為66 mm，整體位移量相對(duì)較小，滑坡處于勻速變形階段，滑坡基本穩(wěn)定，這與模型計(jì)算結(jié)果較為一致，說明所建模型較適用。

圖4 滑坡體下部滑帶部位移監(jiān)測(cè)數(shù)據(jù)

2.3 庫(kù)水位升降對(duì)滑坡穩(wěn)定性影響模擬

庫(kù)水位升降速率是水庫(kù)型滑坡地質(zhì)災(zāi)害發(fā)生的關(guān)鍵影響因素。模擬庫(kù)水位以不同速率升降對(duì)滑坡穩(wěn)定性的影響，找出不同庫(kù)水位、不同升降速率工況下滑坡穩(wěn)定性的主要影響因素，對(duì)于研究水庫(kù)型滑坡成災(zāi)機(jī)理和成災(zāi)模式等具有重要意義。穩(wěn)定性模擬工況見表2。

表2 穩(wěn)定性模擬工況

工況1條件下庫(kù)水位以不同速率上升時(shí)滑坡穩(wěn)定性系數(shù)變化見圖5，工況2條件下庫(kù)水位以不同速率下降時(shí)滑坡穩(wěn)定性系數(shù)變化見圖6，穩(wěn)定性變化規(guī)律見表3。

圖5 145 m標(biāo)高開始上升時(shí)穩(wěn)定性系數(shù)

圖6 175 m標(biāo)高開始下降時(shí)穩(wěn)定性系數(shù)

表3 145/175 m水位標(biāo)高下不同速率升降過程中穩(wěn)定性系數(shù)

由圖5知，在庫(kù)水位從145 m開始上升的過程中，上升速率大于1 m/d時(shí)，滑坡穩(wěn)定性系數(shù)變化趨勢(shì)是先快速增大到極值，再緩慢減小直至逐漸穩(wěn)定;上升速率小于1 m/d時(shí)，滑坡穩(wěn)定性系數(shù)的變化趨勢(shì)是先減小，再緩慢增大直至逐漸穩(wěn)定。庫(kù)水位上升速率大于1 m/d時(shí)，開始階段庫(kù)水的靜水壓力對(duì)滑坡穩(wěn)定性起主導(dǎo)作用，穩(wěn)定性系數(shù)持續(xù)增大，隨著庫(kù)水的入滲，在浮托力增大的同時(shí)，滑帶部位巖土體的物理力學(xué)參數(shù)(c、φ值)逐漸減小，滑坡穩(wěn)定性系數(shù)逐漸降低，并慢慢達(dá)到穩(wěn)定狀態(tài)。庫(kù)水位上升速率小于1 m/d時(shí)，開始階段浮托力及c、φ值改變對(duì)滑坡穩(wěn)定性的作用大于靜水壓力，滑坡穩(wěn)定性系數(shù)持續(xù)減小，隨著庫(kù)水位的上升，靜水壓力的作用逐漸增加，滑坡穩(wěn)定性系數(shù)也逐漸增大直至穩(wěn)定狀態(tài)。由表3知，上升速率越大，達(dá)到穩(wěn)定性極值的時(shí)間越短，至穩(wěn)定狀態(tài)的時(shí)間越長(zhǎng)，反之亦然，穩(wěn)定性系數(shù)最終變化值比初始狀態(tài)增大0.04左右。在庫(kù)水位上升過程中，滑坡整體處于穩(wěn)定狀態(tài)。

由圖6知，在庫(kù)水位從175 m下降的過程中，滑坡穩(wěn)定性變化趨勢(shì)均為先減小至極值，再緩慢增大至穩(wěn)定狀態(tài)。庫(kù)水位下降時(shí)，初始階段庫(kù)水的靜水壓力對(duì)滑坡穩(wěn)定性起主導(dǎo)作用，隨著靜水壓力逐漸減小、局部動(dòng)水壓力增大，穩(wěn)定性持續(xù)減小到極值;達(dá)到極值后，隨著庫(kù)水的疏干，滑帶部位巖土體的物理力學(xué)參數(shù)(c、φ值)逐漸增大，c、φ值增大對(duì)滑坡穩(wěn)定性的作用開始大于靜水壓力和動(dòng)水壓力，滑坡穩(wěn)定性系數(shù)逐漸增大直至達(dá)到穩(wěn)定狀態(tài)。由表3知，下降速度越大達(dá)到穩(wěn)定性極值(最小值)的時(shí)間就越短，其中下降速率大于1 m/d時(shí)，增大至穩(wěn)定狀態(tài)的時(shí)間相差不大，在45 d左右，穩(wěn)定性系數(shù)比初始狀態(tài)減小0.005左右。庫(kù)水位下降過程中，滑坡整體處于基本穩(wěn)定—穩(wěn)定狀態(tài)。

3 結(jié)語(yǔ)

李家坡滑坡屬典型的水庫(kù)型滑坡，庫(kù)水位升降對(duì)滑坡穩(wěn)定性的影響受動(dòng)水壓力、靜水壓力、浮托力、巖土體物理力學(xué)性質(zhì)(c、φ值)改變等多種因素影響，是個(gè)相對(duì)復(fù)雜的過程，穩(wěn)定性變化主要是看哪種影響因素占主導(dǎo)地位。

(1)李家坡滑坡在庫(kù)水位正常自然升降過程中滑坡穩(wěn)定性系數(shù)均表現(xiàn)為先減小后增大的趨勢(shì)，升降過程中滑坡穩(wěn)定性系數(shù)均大于1.15，滑坡總體處于穩(wěn)定狀態(tài)。

(2)庫(kù)水位以不同速率上升的過程中，上升速率越大，庫(kù)水位上升對(duì)滑坡穩(wěn)定性影響越大，上升到175 m前，庫(kù)水的靜水壓力對(duì)滑坡穩(wěn)定性起主導(dǎo)作用，穩(wěn)定性系數(shù)持續(xù)增大到極值，達(dá)到極值后浮托力逐漸增大，c、φ值逐漸減小，穩(wěn)定性系數(shù)也隨之逐漸減小直至達(dá)到相對(duì)穩(wěn)定狀態(tài)。

(3)庫(kù)水位以不同速率下降的過程中，下降速率越大，庫(kù)水位下降對(duì)滑坡穩(wěn)定性影響越大，初始階段庫(kù)水的靜水壓力、動(dòng)水壓力起主要作用，穩(wěn)定性減小到極值后，c、φ值增大，巖土體物理力學(xué)性質(zhì)變化開始起主要作用，穩(wěn)定性系數(shù)逐漸增大直至達(dá)到穩(wěn)定狀態(tài)。

(4)在庫(kù)水位升降過程中，李家坡滑坡體始終處于基本穩(wěn)定—穩(wěn)定狀態(tài);上升過程中上升速率越大達(dá)到最大穩(wěn)定性的時(shí)間越短，減小到穩(wěn)定狀態(tài)的時(shí)間越長(zhǎng)，上升速率大于1 m/d時(shí)最終穩(wěn)定性系數(shù)比初始狀態(tài)增大0.04左右;下降過程中下降速率越快到達(dá)最差穩(wěn)定性的時(shí)間越短，下降速率大于1 m/d時(shí)下降達(dá)到穩(wěn)定狀態(tài)的時(shí)間相差不大，在45 d左右，最終穩(wěn)定性系數(shù)比初始狀態(tài)減小0.05左右。

［1］曾剛.庫(kù)水升降作用下水庫(kù)庫(kù)岸滑坡穩(wěn)定性分析［J］.三峽大學(xué)學(xué)報(bào):自然科學(xué)版，2011，33(4):15 －18.

［2］肖學(xué)沛，李天斌，張志龍.三峽水庫(kù)區(qū)嚴(yán)家坡滑坡成因分析及穩(wěn)定性評(píng)價(jià)［J］.中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2005，16(1):39－43.

［3］連志鵬，譚建民，閆舉生，等.庫(kù)水位變化與降雨作用下庫(kù)岸斜坡穩(wěn)定性分析［J］.安全與環(huán)境工程，2011，18(2):14－22.

［4］薛星橋，張俊義，金裊豪，等.三峽庫(kù)區(qū)淌里滑坡變形特征及影響因素分析［J］.中國(guó)水土保持，2010(6):47－49.

［5］易慶林，易武，尚敏.三峽庫(kù)區(qū)某滑坡變形影響因素分析［J］.中國(guó)水土保持，2009(7):32 －34.

［6］周彬，朱曉強(qiáng)，楊達(dá)源.長(zhǎng)江三峽水庫(kù)庫(kù)岸消落帶地質(zhì)災(zāi)害防治研究［J］.中國(guó)水土保持，2007(11):43－45.

［7］吳爭(zhēng)光.庫(kù)水位變化對(duì)庫(kù)岸邊坡穩(wěn)定性影響研究［J］.災(zāi)害與防治工程，2009(1):1－6.

［8］朱冬林，任光明，聶德新，等.庫(kù)水位變化下對(duì)水庫(kù)滑坡穩(wěn)定性影響的預(yù)測(cè)［J］.水文地質(zhì)工程地質(zhì)，2002，29(3):6－9.

［9］劉才華，陳從新，馮夏庭.庫(kù)水位上升誘發(fā)邊坡失穩(wěn)機(jī)理研究［J］.巖土力學(xué)，2005，26(5):769 －773.