陳文銳

（廣東中灝勘察設(shè)計(jì)咨詢有限公司,廣東 肇慶 526060）

1 引言

氣候變化增加了極端天氣的頻率,長期的氣候變化和海洋溫度是影響降雨模式或強(qiáng)度的關(guān)鍵因素[1]。特別是在山區(qū),高強(qiáng)度降雨經(jīng)常產(chǎn)生土壤侵蝕和相關(guān)的碎屑滑坡,持續(xù)時(shí)間較長的降雨使地表徑流滲入斜坡,從而提高地下水水位,破壞邊坡穩(wěn)定性,引發(fā)滑坡[2,3]。山體滑坡已成為主要我國工程建設(shè)中的主要危害之一。我國大多數(shù)破壞性滑坡均是由長時(shí)間或強(qiáng)降雨誘發(fā)的,大多數(shù)由降雨引起的滑坡均為淺層滑坡,滑面深度約數(shù)米深,且滑動(dòng)迅速。多數(shù)降雨引起的滑坡在沿著陡坡向下移動(dòng)時(shí)轉(zhuǎn)變?yōu)槟嗍?尤其是在西南山區(qū)溝谷地帶,滑坡往往演化為泥石流,造成嚴(yán)重的次生地質(zhì)災(zāi)害[4]。因此,研究滑坡的形成機(jī)理,對(duì)于滑坡災(zāi)害的防治具有重要意義。劉威等[5]以金沙江沃達(dá)滑坡為例,通過離散元數(shù)值模擬軟件分析其滑動(dòng)過程,研究結(jié)果表明,滑坡運(yùn)動(dòng)速度最高可達(dá)38 m/s,滑坡破壞范圍約2.8 km2；賀匯文等[6]采用有限差分?jǐn)?shù)值模擬軟件FLAC對(duì)某路塹邊坡進(jìn)行穩(wěn)定性分析,研究結(jié)果表明,該邊坡失穩(wěn)后,后部坡體仍處于不穩(wěn)定狀態(tài),將繼續(xù)下滑；唐建新等[7]通過UDEC離散元數(shù)值模擬軟件分析某順層巖質(zhì)斜坡失穩(wěn)特征,研究結(jié)果表明,該巖質(zhì)邊坡失穩(wěn)是由于坡體中的軟弱結(jié)構(gòu)面受擾動(dòng),從而發(fā)生滑動(dòng)；汪磊等[8]以浙江麗水市某滑坡為例,利用Geo Studio軟件對(duì)不同的降雨工況下坡體水位進(jìn)行模擬,研究結(jié)果表明,強(qiáng)降雨作用會(huì)顯著降低邊坡的穩(wěn)定性；謝振華[9]基于FLAC3D數(shù)值模擬軟件,以貴州省某礦山排土場為研究對(duì)象,對(duì)降雨強(qiáng)度與邊坡安全系數(shù)的關(guān)系進(jìn)行分析,研究結(jié)果表明,降雨是影響排土場斜坡穩(wěn)定的主要因素；尚敏[10]等利用Geo Studio巖土數(shù)值模擬軟件模擬某滑坡在降雨作用下的滲流場,研究結(jié)果表明,該滑坡失穩(wěn)的主要原因是由于滑坡內(nèi)土體在降雨作用下達(dá)到飽和,導(dǎo)致坡體內(nèi)孔隙水壓力的急劇增加,從而發(fā)生失穩(wěn)。

然而,以上研究均為分析降雨作用對(duì)滑坡的影響機(jī)制,缺少在不同工況下滑坡失穩(wěn)破壞范圍的對(duì)比分析,本文以某邊坡為例,分析其在降雨和天然工況下的變形破壞特征及滑動(dòng)范圍。

2 工程概況

研究區(qū)所在地具有亞熱帶溫帶過度型氣候特征,全年溫暖時(shí)間較長,多年平均氣溫17.1℃ ,一般4月～10月平均氣溫在16℃以上,11月至翌年3月平均氣溫 6℃左右。區(qū)內(nèi)降雨量豐富,年平均降雨量1619.0 mm。研究區(qū)所在地巖土層由第四系全新統(tǒng)素填土,粉質(zhì)粘土及粉砂巖組成。素填土厚為0.60 m～5.20 m,平均約2.44 m。粉質(zhì)黏土層厚為1.6 m～2.4 m,平均厚約1.9 m。粉砂巖厚2.5 m～4.4 m,平均厚度為3.5 m。根據(jù)現(xiàn)場原位測試及室內(nèi)土工試驗(yàn)結(jié)果,各層巖土相關(guān)力學(xué)參數(shù)見表1。

表1 各土層物理力學(xué)參數(shù)

3 降雨條件下滑坡變形破壞特征數(shù)值模擬

針對(duì)降雨條件下的非飽和土的滲流和入滲分析,有限元數(shù)值模擬計(jì)算方法是常用的分析方法之一,ANASYS中的滲流耦合模塊可以較好的用于降雨作用下邊坡穩(wěn)定性進(jìn)行分析。然而,由于實(shí)際滑坡地質(zhì)條件的復(fù)雜性,以及巖土體參數(shù)的不確定性,導(dǎo)致數(shù)值模擬得到的結(jié)果有時(shí)難以解釋。降雨作用下邊坡失穩(wěn)數(shù)值分析涉及許多與土體性質(zhì)、坡度幾何形狀、地下水條件、初始條件和降雨特征相關(guān)的參數(shù)。根據(jù)現(xiàn)場勘測,該滑坡地質(zhì)剖面圖見圖1。

圖1 滑坡剖面圖

由圖1可以看出,滑坡滑動(dòng)前的原始地形較為平緩,平均坡度為5°～15°。原始地形在縱向上有三個(gè)陡峭的臺(tái)階一級(jí)臺(tái)階深20 m～40 m,前緣坡度18°～33°；二級(jí)臺(tái)階發(fā)育高陡坡,高差35 m～49 m,陡坡前緣坡度達(dá)到25°～30°,坡中部沿坡向發(fā)育深溝,深約10 m～20 m。三級(jí)臺(tái)階高差35 m～45 m,坡度20°～30°。這三個(gè)陡坡構(gòu)成了滑坡良好的自由條件,為滑坡的啟動(dòng)提供了有利條件。

如圖2所示,利用Ansys建立了1∶1的滑坡滑動(dòng)前二維地質(zhì)模型,將其網(wǎng)格化并導(dǎo)入Flac3D中,模擬滑坡自然和暴雨條件下的變形破壞特征,揭示滑坡的變形破壞機(jī)制。斷面計(jì)算模型全長1000 m,高300 m,共1130 個(gè)計(jì)算單元。

圖2 滑坡二維數(shù)值模型

剪切破壞面上的強(qiáng)度是逐漸降低的,即滑坡潛在剪切帶不同部分土體的抗剪強(qiáng)度并非同時(shí)發(fā)展到最大值,而是在受集中應(yīng)力最大部位的土體首先發(fā)生剪切或者壓縮破壞,進(jìn)入塑形狀態(tài),隨著降雨入滲,土體抗剪強(qiáng)度降低,導(dǎo)致塑性區(qū)范圍不斷擴(kuò)大,直至整個(gè)滑面完全進(jìn)入塑性區(qū)。利用快速拉格朗日三維數(shù)值計(jì)算計(jì)算出的剪應(yīng)變?cè)隽?根據(jù)數(shù)值模擬得到的最大剪應(yīng)變范圍,可以確定邊坡內(nèi)部的最容易失穩(wěn)的部位,即邊坡的最易潛在滑動(dòng)面,并計(jì)算其安全系數(shù),從而可以判斷邊坡是否穩(wěn)定。根據(jù)數(shù)值模擬結(jié)果得到天然和降雨工況下,滑體內(nèi)剪應(yīng)力分別見圖3、圖4。

圖3 天然狀態(tài)下邊坡剪應(yīng)力云圖

圖4 降雨條件下邊坡剪應(yīng)力云圖

由圖3可知,天然狀態(tài)下,該邊坡在200 m～400 m范圍剪應(yīng)力較大,剪應(yīng)力值在5.5×102MPa～8.5×102MPa之間,可以看出潛在滑面上的剪應(yīng)力已超過土體的最大抗剪強(qiáng)度,發(fā)生失穩(wěn)滑動(dòng)；在500 m～600 m范圍內(nèi)存在另一個(gè)潛在滑動(dòng)區(qū),該滑動(dòng)區(qū)內(nèi)剪應(yīng)力值在1.5×102MPa～2.5×102MPa,剪應(yīng)力小于滑帶土抗剪強(qiáng)度,仍可保持穩(wěn)定狀態(tài)。

由圖4可知,降雨?duì)顟B(tài)下,該邊坡在100 m～400 m范圍剪應(yīng)力較大,剪應(yīng)力值在6.5×102MPa～8.7×102MPa之間,由于降雨入滲導(dǎo)致滑帶土抗剪強(qiáng)度降低,滑面上的剪應(yīng)力超過土體的最大抗剪強(qiáng)度,從而形成滑坡；在400 m～600 m范圍的潛在滑動(dòng)區(qū)內(nèi),剪應(yīng)力值在1.8×102MPa～3.0×102MPa之間,同樣受降雨入滲的影響,導(dǎo)致土體抗剪強(qiáng)度下降,形成局部滑坡。

目前,多數(shù)研究主要采用控制變量法,即針對(duì)單一影響因素進(jìn)行分析,并確定控制參數(shù)。一般而言,坡度陡、高差大的斜坡,更有可能受降雨影響而發(fā)生失穩(wěn)下滑[11]。根據(jù)數(shù)值模擬結(jié)果,可以得出天然狀態(tài)和降雨?duì)顟B(tài)下,坡體內(nèi)的位移變化情況,見圖5。

由圖5可知,天然工況和降雨工況下,位移均隨著時(shí)步增加,且變化趨勢一致。在3000 時(shí)步內(nèi),天然工況下位移呈線性增加,最大位移值為3.12 m。在4000 時(shí)步時(shí)發(fā)生微小下降,隨后位移繼續(xù)增大,在6000 時(shí)步位移值為3.52 m,最后保持穩(wěn)定,說明坡體變形保持穩(wěn)定狀態(tài)。降雨工況下,位移在3000 時(shí)步增加到2.95 m,在6000 時(shí)步位移值為3.52 m,隨后繼續(xù)增加,說明此時(shí)邊坡已經(jīng)發(fā)生失穩(wěn)下滑。

圖5 不同工況下坡面監(jiān)測點(diǎn)位移值

4 結(jié)論

本文利用有限差分軟件FLAC,對(duì)某滑坡在天然和降雨工況下的變形破壞特征進(jìn)行模擬,主要結(jié)論如下：

（1）天然狀態(tài)下,該邊坡在200 m～400 m范圍剪應(yīng)力較大,剪應(yīng)力值在5.5×102MPa～8.5×102MPa之間,最大滑動(dòng)距離為3.52 m。

（2）降雨?duì)顟B(tài)下,由于降雨入滲導(dǎo)致滑帶土抗剪強(qiáng)度降低,剪應(yīng)力值在6.5×102MPa～8.7×102MPa之間,導(dǎo)致該邊坡在100 m～400 m范圍內(nèi),400 m～600 m范圍內(nèi)均發(fā)生了滑坡。

（3）對(duì)比兩種工況下的數(shù)值模擬結(jié)果可知,降雨會(huì)導(dǎo)致滑坡發(fā)生的概率增加,且破壞范圍變大,因此在工程設(shè)計(jì)中應(yīng)加以重視。