朱佳寧, 王建新,2*, 馮成明

(1.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院, 烏魯木齊 830052;2.新疆農(nóng)業(yè)大學(xué)新疆水利工程與水旱災(zāi)害防治實驗室, 烏魯木齊 830052)

渠道岸坡相較于常見的邊坡,既有共同特點,也有其自身特點。其自身特點在于它的狀態(tài)與渠道內(nèi)水位的位置有較大關(guān)聯(lián)。國內(nèi)外已建的輸水工程因岸坡失穩(wěn)而導(dǎo)致災(zāi)害的常有發(fā)生:巴西圣保羅班加拉河流域岸坡在短期強降雨和長期累積降雨天氣下共發(fā)生65處滑坡失穩(wěn)[1];四川省東風(fēng)渠人工渠道[2]在強降雨天氣下,岸坡淺層土壤因抗剪強度降低,進而引起滑坡。河南省陸渾灌區(qū)干渠[3]、溪洛渡庫區(qū)金沙江右岸[4]和三峽木魚包[5]也曾因水位驟變導(dǎo)致了岸坡滑坡,給人民生命、財產(chǎn)安全帶來嚴重損害。因此對輸水明渠進行岸坡穩(wěn)定分析時,不僅需研究降雨對其產(chǎn)生的影響,而且需研究渠道水位狀態(tài)對其產(chǎn)生的影響。

目前,國內(nèi)外眾多學(xué)者針對渠道岸坡穩(wěn)定性已開展大量的研究。沈丹萍等[6]基于Geo-Studio軟件分析膨脹土渠道在不同降雨強度和降雨時長下的邊坡安全系數(shù)和位移的變化規(guī)律,并運用灰色關(guān)聯(lián)度對邊坡的影響因素進行敏感性分析;朱銳等[7]運用SEEP/W得出北疆輸水渠道膨脹土岸坡在各類工況下的暫態(tài)滲流場,并以此將暫態(tài)孔隙水壓力用于渠道岸坡的極限平衡分析中(SLOPE/W 模塊),進而得出在各類工況下的渠道岸坡安全系數(shù);Alonso等[8]采用控制變量法研究了在不同降雨強度和降雨歷時條件下對岸坡穩(wěn)定性的影響;Guo等[9]研究了渠道在水位升降作用下岸坡的安全系數(shù)變化規(guī)律;汪斌等[10]、李清華等[11]采用流固耦合的數(shù)值分析方法,對水位升降作用下岸坡滲流場及應(yīng)力場進行仿真模擬,并得出水位驟降和降雨是影響岸坡安全性的主要原因;文獻[12-13]研究了各類土體渠道在降雨、水位變化時期下的岸坡安全系數(shù)、變形特征和失穩(wěn)機理的變化規(guī)律,但以上研究分析主要集中在膨脹土、鹽漬土、黃土和壤土等渠道岸坡穩(wěn)定性的方面;王海林[14]、李紅娜[15]已對風(fēng)積土渠道有所研究,但主要集中在風(fēng)積土的壓實工藝和施工參數(shù)的選取,關(guān)于沙漠風(fēng)積土渠道的岸坡研究和應(yīng)用仍不成熟。鑒于風(fēng)積土粒徑細小、天然含水率低、結(jié)構(gòu)疏松、保水性差、無黏性等特殊的工程特點[16-17],因而對已建的新疆長距離輸水明渠風(fēng)積土岸坡穩(wěn)定性演變機理尚有深入研究探討的空間。

現(xiàn)以新疆沙漠輸水明渠風(fēng)積土岸坡為實例,選取具有代表性的試驗段風(fēng)積土岸坡作為研究對象,采用Midas GTS/NX軟件中的有限元強度折減法(strength reduction method,SRM)對渠道在未加固和加固兩種工況下的竣工期、降雨期、輸水期和水位驟降4個時期的風(fēng)積土岸坡進行穩(wěn)定性計算,以期在一定程度上為新疆沙漠風(fēng)積土輸水渠道工程建設(shè)提供參考。

1 有限元強度折減法原理

Midas GTS/NX的岸坡穩(wěn)定性分析主要采用了基于有限元法的強度折減法(SRM)[18-19]。

所謂有限元強度折減法,其求解過程與傳統(tǒng)極限平衡法原理本質(zhì)上是一致的,在進行數(shù)值計算過程中,具體是采用逐步折減抗剪強度指標,即原有材料黏聚力c、內(nèi)摩擦角φ分別考慮折減系數(shù)F,從而得到一組新的強度指標cF和φF,替換原來的抗剪強度指標c和φ,重復(fù)運算直至岸坡達到臨界破壞狀態(tài),此時的強度折減系數(shù)F即為該岸坡的安全系數(shù)FS,如式(1)～式(3)所示。

(1)

(2)

τfF=cF+σtanφF

(3)

式中:cF為強度折減后土體虛擬的黏聚力,kPa;φF為強度折減后土體虛擬的內(nèi)摩擦角,(°);τfF為折減后的抗剪強度,kPa。

2 邊坡二維數(shù)值實例分析

2.1 工程實例

新疆古爾班通古特沙漠輸水明渠,由于料源的匱乏,均由當(dāng)?shù)赝馏w填筑,全長為166.5 km,屬于大型長距離輸水工程[圖1(a)]。根據(jù)現(xiàn)場勘探,得出該地巖土體為風(fēng)積土,顏色為淡黃色,具有弱濕陷性,顆粒組成以細砂(0.075～0.25 mm)和中砂(0.25～0.50 mm)為主,含量分別占43.1%～57.8%和22.4%～30.3%,平均粒徑0.24 mm,幾乎無黏性。顆粒圓度不一,為次圓狀～棱角狀。礦物成分以石英顆粒為主,可視為單一成分的特殊性土[圖1(b)]。

圖1 新疆沙漠風(fēng)積土渠道實例

鑒于沙漠風(fēng)積土輸水渠道剖面為對稱結(jié)構(gòu),因而選取渠道剖面的一半進行數(shù)值模擬,并假定渠道岸坡的風(fēng)積土為均勻質(zhì)、各向同性和連續(xù)接觸的彈塑性材料。同時由于渠道建于沙壟、沙丘地帶,地勢起伏不平,一般高差在10～20 m,因此選取具有代表性的橫斷面如圖2所示。該渠道岸坡長19.85 m,高9.24 m,坡度為1∶2.5。風(fēng)積土輸水明渠位于沙漠地區(qū),地下水位埋深達40～150 m,建模時暫不考慮地下水位對其岸坡產(chǎn)生的影響。為查明沙漠風(fēng)積土渠道岸坡是否安全,需對渠道在竣工期、降雨期、輸水期、水位驟降4個時期的風(fēng)積土岸坡進行穩(wěn)定性分析。

圖2 渠道岸坡圖示

2.2 風(fēng)積土的物理力學(xué)指標

根據(jù)開展室內(nèi)土工試驗和查閱渠道現(xiàn)場地質(zhì)勘察報告[20-21],綜合分析得出風(fēng)積土的物理力學(xué)指標(表1)。

表1 風(fēng)積土的物理力學(xué)指標

2.3 岸坡二維模型的建立和網(wǎng)格劃分

Mohr-Coulomb屈服準則是巖土工程中最簡單、應(yīng)用最為普遍的屈服準則。因此在Midas GTS/NX軟件中進行沙漠明渠有限元模擬計算中岸坡風(fēng)積土采用Mohr-Coulomb本構(gòu)模型;同時為了提高網(wǎng)格與網(wǎng)格之間的耦合性,即在易滑帶附近尺寸按四周0.5 m→岸坡0.2 m進行加密過渡,并采用四邊形網(wǎng)格單元。渠道岸坡共劃分為1 642個節(jié)點,1 559個單元,根據(jù)表1的風(fēng)積土數(shù)據(jù),對模型施加對應(yīng)的材料參數(shù),各時期邊界條件詳見表2。

表2 各時期岸坡邊界條件

2.4 研究方案

(1)竣工期:對比在自然狀態(tài)下未加固與加固兩個工況下,風(fēng)積土明渠岸坡的安全系數(shù)、整體位移和等效塑性應(yīng)變區(qū)。

(2)降雨期:對比在西北沙漠地區(qū)強降雨天氣,即降雨強度為2.30×10-3mm/s,降雨時長6 h,總降雨量共50 mm時,未加固和加固兩個工況下風(fēng)積土明渠岸坡的安全系數(shù)、整體位移和等效塑性應(yīng)變區(qū)。

(3)輸水期:對比渠道位于初始水位時,未加固與加固兩個工況下風(fēng)積土明渠岸坡的安全系數(shù)、整體位移和等效塑性應(yīng)變區(qū)。

(4)水位驟降:對比水位在5 d內(nèi)從初始水位驟降到渠底后,未加固與加固兩個工況下風(fēng)積土明渠岸坡的安全系數(shù)、整體位移和等效塑性應(yīng)變區(qū)。

(5)鑒于在降雨和輸水過程中,風(fēng)積土岸坡含水率不斷變化,岸坡材質(zhì)c和φ不斷變化。基于強度折減法(SRM)對風(fēng)積土明渠岸坡進行穩(wěn)定性分析,以獲得岸坡安全系數(shù)和整體位移。

3 岸坡穩(wěn)定性

3.1 未加固時的岸坡數(shù)值模擬計算結(jié)果

在有限元分析時,對風(fēng)積土岸坡采用Mohr-Coulomb 彈塑性模型[22-23],并采取強度折減法進行穩(wěn)定性計算。

《建筑邊坡工程計算規(guī)范》(GB 50330—2013)[24]將渠道岸坡穩(wěn)定性歸納為穩(wěn)定、基本穩(wěn)定、欠穩(wěn)定和不穩(wěn)定4種狀態(tài),寒旱區(qū)長距離輸水明渠風(fēng)積土岸坡穩(wěn)定性岸坡穩(wěn)定安全系數(shù)FS劃分如表3所示。

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

根據(jù)Midas GTS/NX軟件已計算得出未加固工況下竣工期、降雨期、輸水期、水位驟降4個時期渠道岸坡整體位移和安全系數(shù)FS,其具體計算數(shù)值詳見表4。

表4 未加固岸坡分析結(jié)果

由表4得出,在未加固時,風(fēng)積土岸坡在竣工期由于坡體自重,造成整體位移約2.51 cm,安全系數(shù)為1.98;在強降雨工況下,由于雨水滲透岸坡土體內(nèi)部,使風(fēng)積土吸水飽和,容重增加,抗剪強度降低[25-26],使得岸坡整體位移增加至3.98 cm,安全系數(shù)下降至1.93,相較于竣工期安全系數(shù)降低了約2.6%;在輸水期,由于岸坡前緣水壓力增加,致使岸坡的安全系數(shù)上升至2.13,但相較于竣工期,由于坡體前緣水壓力增加,岸坡產(chǎn)生指向岸坡內(nèi)部的水平方向的變形,渠底產(chǎn)生豎直向下的變形[27],使得岸坡整體位移增加至3.38 cm;渠道水位驟降時期,由于岸坡前緣水壓力驟減,渠道發(fā)生卸荷回彈,渠底產(chǎn)生豎直向上的變形,同時岸坡產(chǎn)生指向坡外的水平方向的變形[28],此時岸坡整體位移約5.50 cm,安全系數(shù)為1.70,相較于輸水期安全系數(shù)下降了約20.18%?？⒐て?、降雨期和輸水期的最危險滑弧面大致沿岸坡分布;水位驟降的最危險圓弧有向渠底發(fā)展的趨勢,詳見圖3。

圖3 未加固工況下各時期岸坡整體位移圖

圖4為未加固工況下4個時期風(fēng)積土岸坡等效塑性應(yīng)變區(qū),竣工期、降雨期和輸水期時風(fēng)積土岸坡等效塑性應(yīng)變圖形大體相同,沿岸坡分布;水位驟降時期風(fēng)積土岸坡的等效塑性貫通區(qū)較其他3個時期有沿向渠底發(fā)展的趨勢,因此可大致確定邊坡潛在滑動面呈圓弧分布且已貫通。

圖4 未加固工況下各時期岸坡等效塑性應(yīng)變圖

經(jīng)上述有限元強度折減法數(shù)值分析,在未加固工況下竣工期、降雨期、輸水期、水位驟降4個時期岸坡均處于較穩(wěn)定狀態(tài)。鑒于渠道實際運行過程風(fēng)積土受水荷載作用下,產(chǎn)生壓縮、濕陷變形,其岸坡表層存在約20 cm厚松散砂基處于欠密室狀態(tài),常產(chǎn)生滲透流土破壞[29]。為進一步保證輸水渠道風(fēng)積土岸坡的安全性,方便管理,經(jīng)研究決定對風(fēng)積沙岸坡采取加固措施,以期為結(jié)構(gòu)提供穩(wěn)定基礎(chǔ),從而提高渠道的耐久性。

3.2 加固后的岸坡二維數(shù)值計算結(jié)果分析

為了保證新疆沙漠輸水明渠風(fēng)積土岸坡的穩(wěn)定和減少渠水滲漏損失,則要求渠道堤防加固材料達到一定的強度和防滲性能。墊層材料是渠道防滲結(jié)構(gòu)層穩(wěn)定的基礎(chǔ)。渠道設(shè)計前期,對風(fēng)積土岸坡及渠底采取砂礫石基層、土工格室、化學(xué)固化劑、沙袋和土工膜袋等多種墊層形式進行對比試驗[30],綜合考慮結(jié)構(gòu)穩(wěn)定性、環(huán)境友好性、施工難易性和經(jīng)濟合理性等多因素,最終選擇40 cm 厚的砂礫石作為墊層材料,并在其基礎(chǔ)上鋪設(shè)砂漿找平層、土工膜及C20預(yù)制混凝土襯砌板等措施來進行加固處理[圖1(c)、圖1(d)],典型加固橫斷面布置圖詳見圖5,具體材質(zhì)物理力學(xué)指標詳見表5。

表5 材料的物理力學(xué)指標

渠道岸坡加固工況與未加固工況所賦予的邊界條件相一致,根據(jù)Midas GTS/NX軟件中有限元強度折減法(SRM),研究渠道加固工況下竣工期、降雨期、輸水期、水位驟降4個時期的岸坡整體位移、等效塑性區(qū)和安全系數(shù)的變化,其加固工況下具體計算數(shù)值見表6。

表6 加固岸坡分析結(jié)果

由表6可知,加固工況下竣工期由于岸坡及渠底襯砌結(jié)構(gòu)的約束,整體位移量縮小至約1.33 cm,安全系數(shù)增加至3.16,相較于未加固工況下竣工期整體位移縮小了47.0%,安全系數(shù)增加約37.3%;降雨期由于岸坡和渠底已鋪設(shè)襯砌層結(jié)構(gòu),雨水大量從渠頂滲入,使渠頂表層土體吸水飽和,土體容重增大,抗剪強度降低,岸坡安全系數(shù)下降至3.13,相較于加固工況下竣工期安全系數(shù)縮小了0.9%;輸水期時,由于風(fēng)積土岸坡前緣水壓力增加,岸坡產(chǎn)生指向岸坡內(nèi)部的水平方向的變形,渠底產(chǎn)生豎直向下的變形,但由于襯砌結(jié)構(gòu)的約束,變形較小,坡體整體位移僅為1.21 cm,此時岸坡安全系數(shù)達到3.37,使得岸坡更加穩(wěn)定;水位驟降時期,由于坡體前緣水壓力驟減,渠道土體及襯砌產(chǎn)生卸荷回彈,此時岸坡整體位移增加至3.28 cm,安全系數(shù)為3.09,相較于加固工況下輸水期安全系數(shù)下降了8.3%,整體位移量增加了約63%。加固工況下竣工期、降雨期和輸水期的岸坡最危險滑弧面大致沿岸坡分布;水位驟降時期的岸坡最危險圓弧有向渠底發(fā)展的趨勢,詳見圖6。

圖6 加固工況下各時期岸坡整體位移圖

經(jīng)加固后的風(fēng)積土岸坡由于防滲襯砌結(jié)構(gòu)層的約束,等效塑性應(yīng)變區(qū)由原先圓弧狀變成大體沿襯砌層呈直線分布,原有圓弧狀塑性應(yīng)變區(qū)已不明顯,使滑移帶位移變形得到了較好的限制,避免其發(fā)生過大變形,使岸坡達到極限平衡失穩(wěn)狀態(tài),有效保證了岸坡的整體穩(wěn)定性,較大程度降低了發(fā)生滑坡災(zāi)害的可能,加固工況下岸坡等效塑性區(qū)詳見圖7。

圖7 加固工況下各時期岸坡等效塑性應(yīng)變圖

3.3 未加固與加固岸坡數(shù)值模擬結(jié)果對比

由圖8和圖9分析可知,未加固與加固兩個工況下,輸水期由于風(fēng)積土岸坡前緣水壓力增加,坡體前緣水壓陡增,有效提高渠道岸坡的安全系數(shù);對風(fēng)積土渠道岸坡產(chǎn)生較大位移的分別為降雨期和水位驟降兩個時期。經(jīng)過對渠道岸坡及渠底鋪設(shè)防滲襯砌結(jié)構(gòu)層后,相較于未加固時降雨期和水位驟降兩個時期,風(fēng)積土岸坡總體位移量分別減少了56.4%和40.4%,對渠道整體位移產(chǎn)生了有效約束。由圖8對比可知,未加固時,岸坡處于較為穩(wěn)定的狀態(tài),但影響風(fēng)積土岸坡穩(wěn)定性因素較多,經(jīng)綜合考慮對風(fēng)積土岸坡加固以后,岸坡的安全系數(shù)遠大于1.30,可保障新疆風(fēng)積土輸水明渠長期穩(wěn)定運行。

圖8 未加固與加固岸坡位移對比圖

圖9 未加固與加固岸坡安全系數(shù)對比圖

在未加固和加固兩種工況對比下,渠道水位驟降對岸坡安全系數(shù)影響最大,降雨次之,則要求當(dāng)?shù)毓ぷ魅藛T控制上下游輸水速率,避免水位驟升驟降,同時要定期監(jiān)測,修復(fù)渠道破壞之處,保障沙漠風(fēng)積土渠道的正常輸水。

4 結(jié)論

根據(jù)Midas GTS/NX軟件中的強度折減法,建立了新疆沙漠輸水渠道未加固與加固兩個工況下岸坡的模型,所得結(jié)論如下。

(1)通過分析渠道岸坡未加固和加固兩個工況下竣工期、降雨期、輸水期、水位驟降4個時期的岸坡穩(wěn)定性,得出兩個工況下4個時期岸坡均處于穩(wěn)定狀態(tài)。

(2)雖然未加固工況下渠道岸坡處于穩(wěn)定狀態(tài),但為預(yù)防風(fēng)積土岸坡長期受渠水沖刷和水位驟升驟降的影響,土體滲透性增強,岸坡隨時間的延長抗剪強度降低,從而會產(chǎn)生滑坡的可能。因此經(jīng)分析研究決定,對風(fēng)積沙岸坡進行加固處理,通過數(shù)值模擬,加固效果良好。

(3)未加固時,4個時期的岸坡等效塑性應(yīng)變圖呈圓弧狀,經(jīng)加固后,等效塑性應(yīng)變圖大體沿襯砌呈直線分布,說明已較好控制了滑移帶位移變形,使岸坡不再出現(xiàn)過大變形,從而發(fā)生極限平衡失穩(wěn)狀態(tài),保證了岸坡的整體穩(wěn)定,較大程度降低了滑坡災(zāi)害發(fā)生的可能性。

(4)通過對未加固和加固兩個工況4個時期的岸坡總體位移量和安全系數(shù)對比分析,得出水位驟降對岸坡影響最大,降雨次之,則要求當(dāng)?shù)毓ぷ魅藛T控制上下游輸水速率,避免水位驟升驟降,同時需要定期監(jiān)測,修復(fù)渠道因各項原因?qū)е缕茐闹?保障渠道的正常輸水。