天然和降雨工況下某邊坡失穩(wěn)機制數(shù)值模擬研究

陳文銳

（廣東中灝勘察設計咨詢有限公司,廣東 肇慶 526060）

1 引言

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

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

2 工程概況

研究區(qū)所在地具有亞熱帶溫帶過度型氣候特征,全年溫暖時間較長,多年平均氣溫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)土工試驗結(jié)果,各層巖土相關(guān)力學參數(shù)見表1。

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

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

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

圖1 滑坡剖面圖

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

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

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

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

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

圖4 降雨條件下邊坡剪應力云圖

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

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

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

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

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

4 結(jié)論

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

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

（2）降雨狀態(tài)下,由于降雨入滲導致滑帶土抗剪強度降低,剪應力值在6.5×102MPa～8.7×102MPa之間,導致該邊坡在100 m～400 m范圍內(nèi),400 m～600 m范圍內(nèi)均發(fā)生了滑坡。

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