降雨誘發(fā)某氣田非飽和滑坡機制分析

張鵬程，牟曉東

(中石化石油工程設計有限公司,山東 東營 257000)

0 引言

近幾年來，中國油氣行業(yè)發(fā)展迅速，地面設施建設日新月異。隨著油氣地面建設的不斷發(fā)展，建設環(huán)境越發(fā)惡劣，人類活動與自然環(huán)境矛盾日益凸顯，地質(zhì)災害已呈多發(fā)態(tài)勢[1]?；率怯蜌廨斔凸艿拦こ讨凶顬槌Ｒ姷牡刭|(zhì)災害，其中降雨對滑坡穩(wěn)定性影響受到了國內(nèi)外專家的格外重視[2]。Guzzetti[3-4]通過處理海量滑坡數(shù)據(jù)，得到誘發(fā)滑坡的降雨界限值；Floris[5]建立修正后的降雨量臨界滑坡模型；Crozier[6]提出前期有效降雨量的計算式。國內(nèi)學者在研究降雨對滑坡穩(wěn)定性影響方面也取得了許多成果，如吳樹仁[7]等通過對三峽水庫近幾年的滑坡災害資料進行統(tǒng)計分析，通過對滑坡狀態(tài)預警及滑坡時、空間預報判據(jù)完成三峽水庫滑坡預警研究；吳剛等[8]通過對都江堰五里坡滑坡進行分階段研究，發(fā)現(xiàn)了降雨入滲后地下水對斜坡的作用過程。自2020年雨季以來，四川地區(qū)連續(xù)降雨，導致油氣田地面設施建設區(qū)域發(fā)生嚴重地質(zhì)災害，其中滑坡、泥石流等最為嚴重，給國家資源和人民的生命財產(chǎn)造成嚴重威脅[9-10]。該氣田某閥室也嚴重受此次降雨影響，無法維持正常運營和周邊人民的生活秩序。

1 實例概況

1.1 滑坡概況及強度參數(shù)

研究區(qū)原為斜坡地貌，地形坡度較緩，由于修建下部閥室，對坡體進行開挖切坡，形成了高2.5 m～3.0 m的陡坎，致使坡腳處應力集中分布，同時坡體出露地層為崩坡積形成的粉質(zhì)黏土，其滲透性一般?；聫钠矫嫘螒B(tài)上看，呈不規(guī)則的“圈椅”狀，主滑方向為221°。前緣位于閥室后墻處，相對最大高差為36.3 m?；瑒雍缶壱研纬闪藢捈s30 cm的下錯裂縫；中下部為一斜坡地貌，坡度為28°～33°。該滑坡體后緣存在下錯裂縫，并存在持續(xù)發(fā)展的趨勢。一旦坡體進一步滑塌失穩(wěn)，將直接危害斜坡前緣平臺處該閥室的安全及正常運營，可能導致閥室輸氣管道的破壞進而造成天然氣泄漏，威脅閥室的工作人員及周邊居民的生命財產(chǎn)安全。根據(jù)GB/T 32864—2016滑坡防治工程勘查規(guī)范表B.2續(xù)，本滑坡為小型淺層牽引式土質(zhì)滑坡。具體土層分布及物理力學設計參數(shù)見表1。

表1 巖、土體物理力學性質(zhì)設計參數(shù)

1.2 現(xiàn)場監(jiān)測數(shù)據(jù)

場地內(nèi)共布置監(jiān)測點5個，其中滑坡后緣1個，滑床2個，前緣2個，監(jiān)測點具體見圖1。

現(xiàn)場監(jiān)測數(shù)據(jù)見圖2,圖3。

根據(jù)土壤含水率現(xiàn)場監(jiān)測數(shù)據(jù)可知：2020年7月1日～7月19日，該地區(qū)一直處于間歇式降雨狀態(tài)，滑坡體在間歇性降雨狀態(tài)下于7月2日開始發(fā)生蠕變，直至位移量達到60 mm，位移速度保持穩(wěn)定，持續(xù)變形。在間歇性降雨影響下，該滑坡于7月16日蠕變階段結束，發(fā)生劇烈滑動，位移達到102 mm，并在短短幾個小時迅速增大至201 mm，發(fā)生滑坡。其滑動情況屬于應力小于坡體極限強度，而大于長期強度，在發(fā)生破壞之前經(jīng)過一段較長時間的蠕動過程[11]。

2 非飽和滑坡模型建立

2.1 土-水特征曲線

由于該滑坡在雨季前長時間處于穩(wěn)定狀態(tài)，為完全還原滑坡發(fā)生環(huán)境，本次計算認為降雨初始階段滑坡處于非飽和狀態(tài)。進行非飽和滑坡分析之前必須獲得該地區(qū)土體的土-水特征曲線，以確定該土體含水率與基質(zhì)吸力的關系。試驗獲取土-水特征曲線周期長，實際工程中難以應用；使用GeoStudio中SEEP/W模塊內(nèi)置樣本函數(shù)可快速得到土-水特征曲線[12]?？紤]到該滑坡對氣田閥室的嚴重影響，本文參照四川地區(qū)某粉質(zhì)黏土滑坡研究文獻確定該地區(qū)土-水特征曲線，采用SEEP/W內(nèi)置的樣本函數(shù)法進行擬合。根據(jù)勘察資料，土樣飽和體積含水率為22.3%，得到土-水特征曲線見圖4。

2.2 滲透系數(shù)

非飽和土中，由于顆粒間孔隙由孔隙水和空氣充填，此時，空氣充填的孔隙就成了非傳導通道，水流只能在孔隙水充填的孔隙中流動，使得非飽和土滲透系數(shù)小于飽和土[13]。雖然滲透系數(shù)隨著含水率的變化，但它與含水率有著Van Genuchten滲透系數(shù)公式的對應關系[14]。土-水特征曲線與式(1)相結合可取得非飽和土滲透系數(shù)：

(1)

其中，kw為非飽和滲透系數(shù)；ks為飽和滲透系數(shù)；a，m，n均為擬合參數(shù)，其中n=(1-m)-1；hw為壓力水頭。

2.3 模型建立

本文采用GeoStudio中的SLOPE/W，SEEP/W和SIGMA/W模塊進行降雨模擬計算。其中，SLOPE/W模塊根據(jù)極限平衡法進行穩(wěn)定性分析；SIGMA/W模塊基于有限元進行應力應變分析；SEEP/W則可以進行孔隙水壓力、總水頭、壓力水頭重分布計算。三者可以互相耦合，實現(xiàn)降雨對滑坡穩(wěn)定系數(shù)、應力、變形分布計算。模型建立情況如圖5～圖7所示，因現(xiàn)場跡象明顯，建立模型通過指定后緣拉張裂縫、前緣剪切裂縫位置指定滑動面位置，其下臥地層設為穩(wěn)定地層，邊界為固定約束，即模型左側邊界、下邊界側X,Y方向均不產(chǎn)生變形；根據(jù)現(xiàn)場降雨量監(jiān)測數(shù)據(jù)確定模型降雨強度。

3 滑坡形成機制分析

本文采用GeoStudio中的SLOPE/W,SEEP/W和SIGMA/W模塊進行數(shù)值模擬。其中SLOPE/W模塊是基于極限平衡法進行滑坡穩(wěn)定系數(shù)計算；SEEP/W模塊是采用有限元方法進行水土耦合分析；SIGMA/W模塊是采用有限元方法進行土體、剛性構件應力應變分析[15]。同時3個模塊可以互相結合，實現(xiàn)降雨對滑坡穩(wěn)定性影響分析。

3.1 孔隙壓力重分布

采用SEEP/W進行自然狀態(tài)及降雨條件下孔隙水壓力分布情況分析。根據(jù)現(xiàn)場情況確定0.3 m均勻降雨方式。土-水特征曲線通過SEEP/W中非飽和材料輸入，滲透系數(shù)則通過式(1)進行計算?？紫秹毫Ψ植家妶D8,圖9。

由圖9可知，降雨環(huán)境下，該滑坡地表滲透土層較?。?.3 m降雨量的條件下，地表土體滲透厚度僅為0.7 m。因此，對于該滑坡，降雨主要加劇坡腳應力集中分布，削弱滑坡穩(wěn)定性。由此可得，該滑坡治理過程中，坡腳排水效果至關重要，如果坡腳排水效果不好，則極易引起坡腳剪切破壞。

3.2 滑坡形變量

SEEP/W模擬結束后，對雨后的滑坡變形進行分析計算(見圖10，圖11)。首先將降雨模擬結果導入SIGMA/W模塊中，采用耦合應力/孔隙水壓力進行滑坡變形分析。

由圖可知自然狀態(tài)下滑坡剪應力集中分布于坡腳，符合實際情況。此時滑坡形變量最大為32.53 mm，處于基本穩(wěn)定狀態(tài)。通過選擇3.1最后階段的孔隙水壓力分布進行降雨滑坡變形計算。結果表明：降雨對滑坡穩(wěn)定性影響極大，0.5 d共0.3 m降雨量下滑坡最大位移量達到192.33 mm，直接導致滑坡災害發(fā)生。分析其原因主要有兩點：1)坡面積水增大了坡體重量，是導致滑坡發(fā)生的主要原因之一；2)降雨增大了坡體內(nèi)部水頭，使得土體內(nèi)部發(fā)生滲流(圖11中黑色箭頭表示滲流方向)，使得下部低水頭坡體變形增大，導致牽引式滑坡形成。由此可得，該滑坡治理設計中，坡體截排水效果直接影響滑坡穩(wěn)定性，坡面一旦積水，既施加了外部荷載，也增大了滑坡內(nèi)部水頭，引發(fā)土體內(nèi)部滲流，導致滑坡發(fā)生。

3.3 滑坡穩(wěn)定系數(shù)

將降雨模擬結果導入SLOPE/W模塊中進行滑坡穩(wěn)定系數(shù)計算(見圖12，圖13)。

根據(jù)SLOPE/W計算結果，該滑坡自然狀態(tài)下穩(wěn)定系數(shù)為1.001，處于基本穩(wěn)定狀態(tài)；降雨0.5 d，降雨量0.3 m工況下，穩(wěn)定系數(shù)為0.795，處于不穩(wěn)定狀態(tài)。這一模擬結果與結論不完全符合，其原因是：對于濕陷性黃土、粉土而言，體積含水率在一定范圍內(nèi)，基質(zhì)吸力對土體的強度有增大作用，該作用的大小隨含水率呈某一函數(shù)關系，且存在一個含水率的臨界值，含水率小于臨界值時，隨著含水率增大，對強度的增強作用逐漸降低，當超過臨界值時，含水率增大會減小滑坡穩(wěn)定系數(shù)，但是在該滑坡中，基質(zhì)吸力對其黏聚力影響較小，降雨直接導致土體強度降低，穩(wěn)定系數(shù)減小。由此可得，在該滑坡治理設計中，應充分考慮坡體支護強度，盡量避免考慮基質(zhì)吸力對滑坡穩(wěn)定性影響。

4 結論

本文基于四川某滑坡監(jiān)測數(shù)據(jù)，采用GeoStudio對該滑坡產(chǎn)生機制進行分析，通過工勘資料確定滑坡地層分布情況及各地層物理力學參數(shù)，并參考已發(fā)文獻確定該處滑坡非飽和參數(shù)，綜合考慮現(xiàn)場降雨情況及土體非飽和特性，使得模型最大程度貼近實際情況，確保本次滑坡機制分析實際有效。得到結論如下：

1)降雨對該滑坡穩(wěn)定性影響極大，具體表現(xiàn)為：降雨難以滲入土體深部，只能提高淺層土體含水率，并加劇坡腳應力集中分布，誘發(fā)滑坡發(fā)生；坡面積水對滑坡變形有較大影響，其作用機理為：a.坡面積水施加荷載加重滑坡變形；b.通過增大土體內(nèi)部水頭使得土體內(nèi)部形成滲流，引發(fā)牽引式滑坡發(fā)生。因此，該滑坡治理設計中，坡面截排水措施作用效果至關重要，直接影響治理效果。

2)基質(zhì)吸力對該滑坡強度影響較小，降雨對土體黏聚力影響較小。因此，在治理設計過程中，對于該滑坡應盡量保守設計，避免考慮基質(zhì)吸力對滑坡穩(wěn)定性影響。