宋崗碼頭岸坡穩(wěn)定性分析

劉震 李開(kāi)強(qiáng) 賈靈慈 姚科

(華北水利水電大學(xué)，河南 鄭州 450046)

1 工程概況

1.1 水文地質(zhì)條件

工程區(qū)位于丹江口水庫(kù)東側(cè),岸坡大致為東西走向。高程154m以下地形較為平緩,坡度約為5°～15°；高程154～170m段地形較陡,坡度約為35°～55°；高程168～171m以上地形平坦。工程區(qū)域在大地構(gòu)造上位于秦嶺處于褶皺系東南方向,工程區(qū)總體處于相對(duì)穩(wěn)定的地塊上,地震不易發(fā)。地處北大巴山加里東冒地槽褶皺帶和南秦嶺印支冒地槽褶皺帶處,東部處在南陽(yáng)至襄樊地段。場(chǎng)地區(qū)巖層產(chǎn)狀平緩,傾角處在0°～6°,裂隙以垂直層面的短小裂隙為主。場(chǎng)地區(qū)地表水體主要為丹江水庫(kù),場(chǎng)區(qū)內(nèi)最低排泄基準(zhǔn)面為庫(kù)面。場(chǎng)地區(qū)地下水類型主要為松散砂巖孔隙水。根據(jù)勘察結(jié)果,宋崗碼頭坍岸段岸坡主要為軟質(zhì)巖,部分為巖土混合岸坡。巖土混合岸坡長(zhǎng)170m,約占岸坡總長(zhǎng)的16%,分布在岸坡西段,即海事局所在岸坡段。岸坡上部為厚度1～4.7m粘土；下部為泥巖、粉砂質(zhì)泥巖。巖質(zhì)岸坡長(zhǎng)909m,約占岸坡總長(zhǎng)的84%,主要分布在岸坡中段和東段。本次初步設(shè)計(jì)治理范圍僅包括海事局至老泵站段,治理長(zhǎng)度約為400m,根據(jù)岸坡物質(zhì)組成的差異性及岸坡地形特征,自西向東,由宋崗碼頭至宋崗泵站將岸坡分為2段,其地勢(shì)大致特征如圖1所示。

圖1 宋崗碼頭庫(kù)岸圖

岸坡主要為巖土混合岸坡,縱長(zhǎng)170m。高程154m以上,岸坡坡度35°～55°。坡頂覆蓋約有1～3m第四系中更新統(tǒng)沖積粘土,具有弱膨脹性；岸坡基巖產(chǎn)狀大約為水平,屬新近系,上部基巖為泥巖、泥質(zhì)粉砂巖,具有中等膨脹性,中下部基巖為厚層疏松砂巖,疏松砂巖結(jié)構(gòu)不穩(wěn)定,容易產(chǎn)生塌陷。高程154m以下,岸坡平緩其坡度約為10°～15°,覆蓋0.5～3m第四系上更新統(tǒng)沖積粘土,岸坡基巖以新近系粉砂質(zhì)泥巖為主,具有中等膨脹性,產(chǎn)狀近水平。

由于目前岸坡已經(jīng)進(jìn)行了部分治理,故下文模擬將主要以改性后自然狀態(tài)進(jìn)行。其剖面圖如圖2所示。

圖2 岸坡剖面圖

1.2 岸坡膨脹土特性

作為一種高塑性粘性土,膨脹土常在天然狀態(tài)下處于非飽和狀態(tài),對(duì)岸坡的穩(wěn)定性有極大影響。膨脹土一般由比較親水的礦物組成,具有很強(qiáng)的脹縮性。

膨脹土主要成分包括蒙脫石、伊利石和高嶺石等親水類礦物,因此吸水膨脹和失水收縮2個(gè)變形特性就是其基本特性之一。膨脹土是一種擁有多裂隙性和強(qiáng)脹縮性的非飽和土,在天然狀態(tài)下一般較為堅(jiān)硬。膨脹土在許多國(guó)家都有分布,而我國(guó)更是分布眾多。

自1938年的美國(guó)開(kāi)墾局因一個(gè)基礎(chǔ)工程首次提出膨脹土問(wèn)題[1],人們開(kāi)始逐漸關(guān)注起膨脹土問(wèn)題。由于膨脹土遇水反復(fù)變形的特性,深受膨脹土影響的工程包括公路、鐵路、岸坡和房屋地基等。由于膨脹土的強(qiáng)脹縮性和多裂隙性所帶來(lái)危害的長(zhǎng)期性、潛在性和難以處理,導(dǎo)致的經(jīng)濟(jì)損失是難以估計(jì)的。有資料顯示,美國(guó)每年需要花費(fèi)高達(dá)20億美元來(lái)彌補(bǔ)膨脹土帶來(lái)的損失,全世界每年由于膨脹土問(wèn)題導(dǎo)致的經(jīng)濟(jì)損失則高達(dá)150億美元[2]。因此考慮膨脹土對(duì)于水庫(kù)岸坡影響成為不可忽視的工程考量。

我國(guó)膨脹土在長(zhǎng)江流域和黃河流域分布比較集中,而在其分布廣闊地區(qū)又以水庫(kù)岸坡工程較多,由于膨脹土的顯著脹縮特性,故要根據(jù)膨脹土的力學(xué)及物理學(xué)性質(zhì)注意膨脹土對(duì)于工程的不良影響。在我國(guó),膨脹土的分布面積廣泛,主要分布在湖北、河南、廣西、陜西、貴州等地[3],其中主要集中于中西區(qū)域。膨脹土具有很高的膨脹潛勢(shì),這與其含水量的大小及變化有關(guān)。如果其含水量保持不變,則不會(huì)有體積變化。在工程施工中,建造在含水量保持不變的粘土上的構(gòu)造物不會(huì)遭受由膨脹而引起的破壞。當(dāng)黏土的含水量發(fā)生變化,就會(huì)立即產(chǎn)生垂直和水平2個(gè)方向的體積膨脹。所以對(duì)膨脹土岸坡的穩(wěn)定性分析對(duì)于膨脹土相關(guān)工程的應(yīng)用是有一定意義的。

岸坡膨脹土由于其特殊的性質(zhì)極易造成滑坡,而滑坡作為一種危害程度高、規(guī)模大、次數(shù)頻繁且有規(guī)律可循的自然災(zāi)害。這種現(xiàn)象指的是研究區(qū)域的巖土體,在多種地質(zhì)條件的共同作用下,導(dǎo)致原有的穩(wěn)定狀態(tài)向不穩(wěn)定狀態(tài)進(jìn)行轉(zhuǎn)化,在坡體自重或者坡體上方荷載的作用下,巖土體沿著軟弱面向下發(fā)生滑動(dòng)導(dǎo)致失穩(wěn)現(xiàn)象的發(fā)生。這種不良地質(zhì)災(zāi)害的頻繁發(fā)生,致使道路中斷、交通受阻、河道因山體滑落發(fā)生堵塞、多數(shù)村莊住宅被無(wú)情的摧毀,人們因此流離失所失去生命。所以在相關(guān)工程中應(yīng)極力避免滑坡的不良影響,保證人民的生命財(cái)產(chǎn)安全。

1.3 Geostudio極限平衡分析

對(duì)于庫(kù)岸邊坡失穩(wěn)的研究主要有數(shù)值仿真擬和模型試驗(yàn)2種方法,數(shù)值模擬以低成本、單快速等特性而被廣泛運(yùn)用,但需精確選擇計(jì)算邊界條件及土體物理力學(xué)特性等參數(shù)；模型試驗(yàn)觀性強(qiáng)、可靠性高,且可與數(shù)值仿真相互驗(yàn)證。膨脹土岸坡穩(wěn)定分析目前在理論上可分為2大類,常規(guī)的飽和土理論,根據(jù)水分對(duì)膨脹土邊坡失穩(wěn)的影響程度,人為地改變膨脹土強(qiáng)度來(lái)進(jìn)行分析；運(yùn)用非飽和土理論,利用膨脹土強(qiáng)度在吸水后的衰減特性來(lái)分析降雨對(duì)邊坡穩(wěn)定性的影響。在分析計(jì)算方法上,研究岸坡穩(wěn)定性的傳統(tǒng)方法主要有極限平衡法、極限分析法、有限元法等[4]。

極限平衡法是工程分析領(lǐng)域常用的分析方法之一,其主要理論分析依據(jù)是莫爾-庫(kù)倫強(qiáng)度理論,主要思想是建立每個(gè)條塊力平衡方程和力矩平衡方程,從而求解邊坡的穩(wěn)定性系數(shù)[5]。假定巖土體為剛體,不產(chǎn)生變形但傳遞力通過(guò)靜力(或力矩)平衡分析,獲取滑動(dòng)面的反力,進(jìn)而進(jìn)行穩(wěn)定性分析。同時(shí),極限平衡法包含了很多種方法,在進(jìn)行穩(wěn)定性計(jì)算分析時(shí),需要考慮工程的具體類型來(lái)選取適合的計(jì)算方法。為了更好地反映實(shí)際情況,本文選用了幾種方法進(jìn)行對(duì)比分析,從而能夠得出比較適宜的極限平衡計(jì)算方法。20個(gè)世紀(jì)70年代,Fredlund教授等通過(guò)大量研究應(yīng)用,開(kāi)發(fā)出了GeoStudio數(shù)值模擬軟件。軟件的應(yīng)用主要面向巖土、采礦、水利、地質(zhì)等領(lǐng)域。經(jīng)過(guò)多年的開(kāi)發(fā),GeoStudio軟件已經(jīng)是解決巖土工程方面問(wèn)題的專業(yè)軟件。目前的版本主要包含8大模塊,SLOPE/W、SEEP/W、SIGMA/W、QUAKE/W、TEMP/W、CTRAN/W、AIR/W、VADOSE/W。這 些模塊可以模擬幾乎所有的巖土問(wèn)題,如常見(jiàn)的邊坡的穩(wěn)定性分析問(wèn)題,地下水滲流問(wèn)題等。其中SLOPE/W模塊是一種基于極限平衡理論用于對(duì)邊坡穩(wěn)定性的計(jì)算而開(kāi)發(fā)的一種模擬分析軟件；SEEP/W模塊是基于有限元方法用于解決土體滲流方面問(wèn)題的軟件。GeoStudio在模型建立方面優(yōu)勢(shì)突出,只需要建立一個(gè)模型就可以供所有模塊使用；高效的建模功能,支持快速建立模型,可以復(fù)制模型,而且還可以調(diào)用上級(jí)目錄的分析結(jié)果,用來(lái)實(shí)現(xiàn)一個(gè)模型多種工況下的分析；后處理功能也非常的豐富與強(qiáng)大,分析結(jié)果便于輸出和整理。

2 不同工況下岸坡穩(wěn)定性模擬

2.1 自重狀態(tài)下岸坡穩(wěn)定性分析

本節(jié)采用Geostudio軟件中SLOPE/W模塊對(duì)滑坡穩(wěn)定性進(jìn)行分析,該模塊計(jì)算原理基于極限平衡法,包括瑞典條分法、Bishop法、Janbu法、Morgenstern-Price法等。首先按照物探數(shù)據(jù)建立模型,后采用模擬計(jì)算對(duì)自然狀態(tài)下該岸坡進(jìn)行穩(wěn)定性的分析,該模擬建立為高程130m以上的宋崗碼頭岸坡模型,從上至下土層主要分為膨脹土層、泥巖層、泥質(zhì)粉砂巖層,后將各個(gè)土層參數(shù)輸入,確定大致滑移面的滑動(dòng)范圍,計(jì)算得出自然狀態(tài)下的岸坡穩(wěn)定性系數(shù)。以此系數(shù)分析自然狀態(tài)下宋崗碼頭岸坡是否處于穩(wěn)定狀態(tài),是否會(huì)產(chǎn)生滑坡。使用Geostudio的計(jì)算穩(wěn)定性系數(shù)功能進(jìn)行穩(wěn)定性計(jì)算。計(jì)算參數(shù)選取見(jiàn)表1。

表1 巖土物理參數(shù)

對(duì)滑坡進(jìn)行穩(wěn)定性計(jì)算,強(qiáng)度準(zhǔn)則選用Mohr-Coulomb準(zhǔn)則,分別使用Morgenstern-Price法、Bishop法、Janbu法3種方法進(jìn)行計(jì)算,得到自重狀態(tài)下計(jì)算結(jié)果見(jiàn)圖3～5和表2。

表2 滑坡穩(wěn)定性系數(shù)對(duì)比表

圖3 Bishop法穩(wěn)定性系數(shù)

在自然條件下Bishop計(jì)算岸坡穩(wěn)定性系數(shù)還是處于較高的水平,整體岸坡處于穩(wěn)定狀態(tài),見(jiàn)圖3?？赡艹霈F(xiàn)的滑移面處于膨脹土層的斷層,符合組成斜坡的巖、土體只有被各種構(gòu)造面切割分離成不連續(xù)狀態(tài)時(shí),才有可能向下滑動(dòng)的條件。

由自然條件下Janbu法計(jì)算穩(wěn)定性系數(shù)可知,此岸坡仍處于穩(wěn)定狀態(tài),見(jiàn)圖4?；泼嫒匀皇俏挥谂蛎浲翆拥臄鄬犹?且滑移面較小。

圖4 Janbu法穩(wěn)定性系數(shù)

在自然條件下Morgenstern-Price法計(jì)算穩(wěn)定性系數(shù)與Bishop法計(jì)算結(jié)果十分相近,岸坡都處于穩(wěn)定狀態(tài),見(jiàn)圖5。

圖5 Morgenstern-Price法穩(wěn)定性系數(shù)

軟件采用系統(tǒng)自動(dòng)搜索最危險(xiǎn)滑面的方法確定滑面的位置,根據(jù)以上計(jì)算結(jié)果可以看出,計(jì)算出的穩(wěn)定性系數(shù)大小由于使用的方法不同有所差異。計(jì)算結(jié)果按數(shù)值大小排列為Bishop法＞Morgenstern-Price法＞Janbu法,使用3種不同計(jì)算方式計(jì)算出的滑動(dòng)面形狀基本都為一較小的圓弧形,位于滑坡膨脹土層斷層處。對(duì)比表3數(shù)據(jù),可以得出結(jié)論此岸坡在自然狀態(tài)下處于穩(wěn)定狀態(tài)。

表3 滑坡穩(wěn)定性狀態(tài)劃分

2.2 降雨工況下岸坡穩(wěn)定性分析

在降雨工況下,雨水滲入土體,使土體的飽和度增加,劣化了巖土體的抗剪強(qiáng)度,使得滑坡體自重增加,加大了滑坡體自重應(yīng)力沿滑移面的滑動(dòng)力,水浸入土體,對(duì)滑動(dòng)面起到了潤(rùn)滑的作用,進(jìn)一步增加了滑坡失穩(wěn)的可能性。模型兩側(cè)地下水位以上按零流量邊界處理,滑坡其他區(qū)域設(shè)為不透水的邊界條件。在滑坡表面設(shè)置降雨邊界,SEEP/W模塊會(huì)根據(jù)降雨強(qiáng)度和滲透系數(shù)的關(guān)系進(jìn)行計(jì)算。在降雨強(qiáng)度大于滲透系數(shù)時(shí),按水頭邊界處理；當(dāng)降雨強(qiáng)度小于滲透系數(shù)時(shí),按流量邊界處理。

本節(jié)主要研究了滑坡在不同降雨強(qiáng)度條件下穩(wěn)定性系數(shù)變化的情況。根據(jù)當(dāng)?shù)貧庀髼l件,結(jié)合降雨強(qiáng)度等級(jí)劃分,設(shè)定了以下幾種降雨工況:小雨(10mm·24h-1)、中雨(25mm·24h-1)、大 雨(50mm·24h-1)、暴 雨(100mm·24h-1)、大暴雨(200mm·24h-1),見(jiàn)表4。

表4 降雨強(qiáng)度等級(jí)劃分

圖6 10mm降雨穩(wěn)定性系數(shù)

由降雨分析可知,在10mm降雨條件下,岸坡穩(wěn)定性系數(shù)有所下降,證明降雨對(duì)于岸坡穩(wěn)定性有負(fù)面影響,故以此為例可進(jìn)行不同降雨條件下岸坡安全性系數(shù)模擬。

根據(jù)穩(wěn)定性計(jì)算可以知道,岸坡在降雨影響下安全性系數(shù)會(huì)降低,在100mm降雨下處于基本穩(wěn)定狀態(tài),而在200mm大暴雨下處于不穩(wěn)定狀態(tài),所以可知宋崗碼頭岸坡治理需要注意大暴雨情況。

表5 降雨工況穩(wěn)定性計(jì)算成果表

3 結(jié)論

根據(jù)現(xiàn)場(chǎng)對(duì)岸坡的工程地質(zhì)勘察,建立了典型剖面的數(shù)值計(jì)算模型,利用Geostudio軟件分別使用Morgenstern-Price法、Bishop法、Janbu法3種方法進(jìn)行在自重狀態(tài)、不同降雨強(qiáng)度工況下的穩(wěn)定性分析模擬計(jì)算,得到了穩(wěn)定性系數(shù)。獲得的主要結(jié)論如下。

天然狀態(tài)下宋崗碼頭岸坡穩(wěn)定性系數(shù)為1.998,滑坡處于穩(wěn)定狀態(tài),故在對(duì)于岸坡治理是有效的。

降雨對(duì)滑坡的穩(wěn)定性具有一定的削弱作用,暴雨會(huì)使滑坡出現(xiàn)失穩(wěn)情況。隨著降雨量的增加,滑坡的穩(wěn)定性系數(shù)也在不斷的降低。具體影響作用表現(xiàn)為:在降雨過(guò)程中,隨著雨水的下滲,坡體內(nèi)體積含水量增加,使得坡體自重增加,滑動(dòng)面的粘聚力和內(nèi)摩擦角都會(huì)降低,最終共同影響了滑坡的穩(wěn)定性。

目前宋崗碼頭岸坡處于穩(wěn)定狀態(tài),但在分析時(shí)可看出對(duì)于穩(wěn)定性影響較大的是降雨影響。所以應(yīng)該加強(qiáng)降雨情況下對(duì)于岸坡穩(wěn)定性影響的控制。