馬慧敏　何向東　張帥　劉吉永

摘 要：膨脹土（巖）工程特性特殊，易造成邊坡失穩(wěn)、滑塌，因此對膨脹土（巖）渠坡的穩(wěn)定性進行分析尤為重要。針對膨脹土（巖）渠段不同的破壞模式提出了相應的邊坡穩(wěn)定分析方法，并給出了邊坡穩(wěn)定分析時參數(shù)取值的經驗方法：對于有明顯裂隙面的裂隙強度控制下的渠坡滑動采用極限平衡法中的摩根斯坦-普拉斯法（Morgenstein-Price Method）進行計算（滑裂面為裂隙結構面），對于無明顯裂隙面的渠坡采用極限平衡法中的畢肖普法（Bishop）進行計算（搜索最危險滑動面）。以南水北調中線工程干渠某膨脹土（巖）渠段為例進行渠坡穩(wěn)定分析，指出：對于南水北調中線工程土體參數(shù)，弱膨脹土黏聚力取峰值強度，中、強膨脹土（巖）體強度采用分帶選取原則，在大氣影響帶（0～2.5 m）范圍內取殘余強度，在過渡帶（2.5～4.5 m）范圍內黏聚力按峰值強度的0.5～0.7倍取值、內摩擦角按峰值強度對應內摩撐角的0.7～0.8倍取值，4.5 m以下非影響帶黏聚力取峰值強度。

關鍵詞：膨脹土（巖）邊坡;穩(wěn)定分析;參數(shù)取值;南水北調中線工程

中圖分類號：TV91?? 文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.03.029

Analysis and Study on Slope Stability of Expanded Soil Channel Section of Large Water Pipeline

MA Huimin1， HE Xiangdong2， ZHANG Shuai1，LIU Jiyong1

（1.Henan Water & Power Engineering Consulting Co.， Ltd.， Zhengzhou 450016， China;

2.Henan Provincial Water Conservancy Research Institute， Zhengzhou 450008， China）

Abstract： The canal slope or bottom of the trunk canal of the South-to-North Water Diversion Project involves about one-third of the open channel section of the expansion rock（soil）， and the engineering geological conditions are complex. Due to its special engineering characteristics， expansive rock（soil） is prone to slope instability and collapse， and the problem of expansive rock（soil） encountered in waterconservancy projects（especially channel engineering） is more and more difficult to deal with. In addition， the amount of engineering， processing investment and processing are also more difficult. Therefore， the stability of the expansion rock（soil） slope is one of the key problems in the design of the expansion rock（soil） channel. In this paper， the corresponding slope stability analysis method was proposed for the different failure modes of expansive rock（soil） channel. The empirical method of the parameter value of slope stability analysis is given： For the channel slope sliding under the control of the fracture strength with obvious fissure surface， Pulasifa（Morgenstein-Price Method） in the limit equilibrium method is used to calculate（sliding surface is fracture structure surface） for drainage slopes with no obvious fissure surface， Bixiaopufa in the limit equilibrium method is used to calculate（search for the most dangerous sliding surface）; The strength of the soil body in the middle line of the South-to-North Water Diversion takes the peak strength of the weak expansion soil， and the principle of band selection is adopted for the strength of the medium and strong expansion（rock） soil body： The residual strength is taken within the range of the atmospheric influence zone（0 to 2.5 m）; In the transition zone（2.5-4.5 m）， the C value is 0.5 to 0.7 times the peak strength， and the φ value is 0.7 to 0.8 times; 4.5 m below non-effect band peak strength.

Key words： expansive soil slope; stable analysis; parameter value; Middle Route Project of South to North Water Transfer

南水北調中線工程從丹江口水庫引水至北京、天津，輸水距離長達1 400 km，年均調水約95億m3。干渠工程大多采用明渠輸水，渠坡或渠底涉及膨脹土（巖）的渠段累計約380 km，約占明渠段的1/3，工程地質條件復雜，膨脹土（巖）的特殊性主要反映在土體的脹縮性、裂隙性和超固結特性上，在水利工程中水的存在使其危害更加突出。膨脹土（巖）渠道設計需要解決的主要問題之一為膨脹土（巖）渠段的邊坡設計：綜合考慮膨脹土（巖）級別、土體結構與工程特性、環(huán)境地質條件、大氣影響深度等影響因素，建立合理的分析計算模型以及選取合理的土體力學參數(shù)進行膨脹土（巖）邊坡穩(wěn)定分析計算。

1 傳統(tǒng)膨脹土（巖）渠坡的穩(wěn)定分析理論和計算方法

計算膨脹土（巖）邊坡穩(wěn)定性的傳統(tǒng)方法有極限平衡法和有限元法。

極限平衡法是目前實際工程穩(wěn)定分析中最常采用的方法。坡比為1∶4～1∶6的膨脹土邊坡很多會發(fā)生淺層滑坡，采用通常的極限平衡法計算時往往顯示邊坡安全不會滑坡，這時不考慮膨脹土（巖）的強度變化、裂隙因素等而一味采用極限平衡法顯然不合適[1]。以往研究中，為反映水對膨脹土的影響，一般引入非飽和土強度隨土體含水量的增大而降低這一概念，通過改變土體強度進行膨脹土（巖）邊坡穩(wěn)定分析[2]。大量的工程實踐證明，不考慮膨脹土（巖）的基本特性、土體結構、環(huán)境地質條件以及邊坡的工程運行特點等因素，統(tǒng)一采用傳統(tǒng)的極限平衡法進行邊坡穩(wěn)定分析是不合適的。

有限元法是把土坡當成變形體，按照土的變形特性計算土坡的應力分布，然后引入圓弧滑動面的概念，驗算滑動土體的整體抗滑穩(wěn)定性。目前有限元法在膨脹土邊坡穩(wěn)定性分析方面的應用主要是通過非飽和土滲流場的有限元分析來定性解釋降雨期間或降雨后這類邊坡的失穩(wěn)機理[3-9]。

2 南水北調中線采用的膨脹土（巖）渠坡穩(wěn)定分析理論和計算方法

南水北調中線工程有近1/3的渠段涉及膨脹土（巖），膨脹土（巖）因其特殊的工程特性而易造成渠坡失穩(wěn)，對工程的安全運行影響很大，而且其處理難度、處理的工程量和投資較大。為此，國家“十一五”科技支撐計劃項目

“南水北調工程若干關鍵技術研究與應用”對膨脹土（巖）渠段的邊坡設計問題進行了專項研究，其研究成果在南水北調中線工程中進行了應用，主要研究成果如下。

（1）提出膨脹土（巖）渠坡的兩種主要破壞模式。從膨脹土（巖）工程的破壞現(xiàn)象和破壞特征中歸納出兩種主要的破壞模式，即裂隙強度控制下的渠坡滑動和膨脹作用下的渠坡滑動，前者破壞以沿裂隙面或結構面的滑動為主，多表現(xiàn)為重力式滑動，后者以膨脹變形為驅動力，多表現(xiàn)為淺表層膨脹變形后的牽引式滑動。

（2）針對不同的破壞模式提出渠坡穩(wěn)定分析理論和分析方法。對于裂隙強度控制下的渠坡滑動可以采用極限平衡理論;對于膨脹作用下的渠坡滑動，應采用考慮膨脹變形的有限元分析理論。

（3）提出穩(wěn)定分析參數(shù)測試和取值方法。對于裂隙強度控制下的渠坡穩(wěn)定，應采用裂隙面強度指標。

由于有限元法存在很多假定，本構關系不準確，因此在南水北調中線工程（沙河南至黃河南段以及黃河北至漳河南段）膨脹土（巖）邊坡穩(wěn)定分析時仍采用極限平衡法，對于有明顯裂隙面的裂隙強度控制下的渠坡滑動采用極限平衡法中的摩根斯坦-普拉斯法（Morgenstein-Price Method）進行計算（滑裂面為裂隙結構面），對于無明顯裂隙面的渠坡采用極限平衡法中的簡化畢肖普法（Bishop）進行計算（搜索最危險滑動面）。以下以簡化畢肖普法為例推導極限平衡法的計算公式，從圓弧滑動體內取出土條i進行分析。作用在條塊i上的力除了重力Wi外，滑動面上有切向力Ti和法向力Ni，條塊的側面分別有法向力Pi、Pi+1和切向力Hi、Hi+1（見圖1）。

圖1 畢肖普法條塊作用力分析

若條塊處于靜力平衡狀態(tài)，根據(jù)豎向力平衡條件（∑FZ=0），應有：

Nicos θi=Wi+ΔHi-Tisin θi（1）

根據(jù)滿足安全系數(shù)為Fs時的極限平衡條件，有：

Ti=1Fs（cili+Nitan φi）（2）

將式（2）代入式（1），整理得

式中：li為條塊i沿滑面底部的長度;ci、φi為條塊i土體黏聚力及內摩擦角。

考慮整個滑動土體的整體力矩平衡條件，各土條的作用力對圓心力矩之和為零。這時條間力Pi和Hi成對出現(xiàn)，大小相等，方向相反，相互抵消，對圓心不產生力矩?；瑒用嫔系恼龎毫i通過圓心，也不產生力矩。因此，只有重力Wi和滑動面上的切向力Ti對圓心產生力矩：

∑Widi=∑TiR（5）

進一步假定ΔHi=0，簡化得

Fs=∑1mθi（cibi+Witan φi）∑Wisin θi（6）

以上是簡化Bishop法的分析過程，其他方法的分析原理與之類似。Bishop法只考慮了力的平衡，而忽略了力矩的平衡條件，只適用于圓弧滑裂面。

Morgenstern等對極限平衡理論進行研究發(fā)展，建立了能同時滿足力矩和力的平衡來求解安全系數(shù)的方程，而且滑裂面可以是任意形狀。摩根斯坦-普拉斯法對相鄰條間作用力進行假定，給出了條間合力的作用位置，通過改變條間合力的作用方向以求得最佳解[10]。

3 南水北調中線膨脹土（巖）邊坡的穩(wěn)定分析參數(shù)和取值方法

渠坡穩(wěn)定計算涉及的主要指標是巖土體的抗剪強度。膨脹土（巖）強度指標的選取是個復雜的問題，目前還沒有統(tǒng)一的標準。膨脹土（巖）峰值強度較高，但是從失穩(wěn)的膨脹土（巖）邊坡反算出的抗剪強度卻遠遠低于峰值強度。這一現(xiàn)象早就引起了人們的注意，經過多年研究后，人們普遍認為膨脹土（巖）的脹縮性、裂隙性和超固結性是造成強度衰減的重要原因。目前，膨脹土（巖）的參數(shù)試驗值和設計實際取值脫節(jié)。一方面，由膨脹土（巖）原狀樣試驗得到的數(shù)值往往偏大，設計人員不敢直接采用;另一方面，對于某些已滑動的渠坡，即使采用膨脹土（巖）的殘余強度計算，安全系數(shù)仍然很大。實際工程中，在膨脹土（巖）地區(qū)即使邊坡坡比為1∶4～1∶6也會發(fā)生滑坡。如湖北鄖縣邊坡坡比約1∶5的地區(qū)長期存在緩慢的蠕變;南水北調中線工程南陽膨脹土試驗段中某區(qū)左岸開挖成1∶2邊坡后一年自然坡比已演化為1∶6，右岸1∶1.5邊坡已演變?yōu)?∶4左右。一些規(guī)范中提出采用長期強度或殘余強度進行渠坡設計，盡管這樣的設計偏于安全，但其取值沒有理論依據(jù)和工程實踐檢驗，且實際工程中這樣的取值仍可能發(fā)生滑坡。

針對以上膨脹土（巖）邊坡穩(wěn)定計算時抗剪強度的經驗取值法，“十一五”國家科技支撐項目“膨脹土（巖）渠坡破壞機理及分析方法研究”課題組對膨脹土（巖）渠坡穩(wěn)定分析參數(shù)測試及取值方法進行分析后提出，膨脹土（巖）的強度有土塊強度、裂隙面強度和土體強度之分。土塊強度是指土中不含裂隙面時的土體強度;裂隙面強度是指膨脹土（巖）中“非脹縮裂隙面”或地層結構面的強度。對于脹縮裂隙無裂隙面強度，當土體不存在裂隙時，土塊強度也是土體強度;對于非脹縮裂隙膨脹土的原狀土才談土塊強度和裂隙面強度。

南水北調中線工程膨脹土（巖）渠段設計中提出了裂隙強度和土體強度兩個概念，設計人員依據(jù)渠坡不同的破壞模式進行參數(shù)取值。對于存在明顯裂隙面的渠段采用裂隙面強度，其余采用土體強度。裂隙面強度采用由現(xiàn)場原狀樣測得的土工試驗數(shù)值并結合反算的邊坡穩(wěn)定抗剪強度指標進行取值。如：某渠段施工期在邊坡開挖卸荷后發(fā)生邊坡滑坡，此處地層結構為黏土巖，初步設計時黏土巖計算參數(shù)黏聚力為16 kPa，內摩擦角為17°，根據(jù)滑塌體滑裂面的位置對邊坡采用摩根斯坦-普拉斯法進行復核，計算工況為渠內無水、穩(wěn)定滲流，計算安全系數(shù)為1.854，滿足要求，但顯然不符合實際情況，為此進行滑裂面參數(shù)反演分析。邊坡穩(wěn)定安全系數(shù)采用1.0，反演出來的滑裂面力學參數(shù)為：黏聚力8.5 kPa，摩擦角10°。土體強度總體趨于中、強膨脹土（巖）取殘余強度或流變強度，弱膨脹土取峰值強度。土的應力—應變曲線上最大剪應力的強度為峰值強度，峰值后若變形繼續(xù)發(fā)展，偏應力將不斷降低，當變形很大時，其趨于穩(wěn)定值，稱為殘余強度，是反映土體滑動后滑動面的強度。針對南水北調中線工程典型土（巖）體，通過反復剪試驗得到的成果見表1[1]。

南水北調中線工程對中、強膨脹土（巖）的強度采用分帶法進行邊坡穩(wěn)定計算參數(shù)取值，根據(jù)以上試驗成果以及工程經驗，膨脹土（巖）抗剪強度的分帶選取原則：在大氣影響帶（0～2.5 m）范圍內，取殘余強度;在過渡帶（2.5～4.5 m）范圍內，黏聚力c按峰值強度的0.5～0.7倍取值，內摩擦角φ按峰值強度的0.7～0.8倍取值;4.5 m以下在過渡帶下限的非影響帶取峰值強度。

4 南水北調中線膨脹土（巖）渠坡典型案例分析

4.1 典型實例計算工況

選擇南水北調中線某段弱膨脹巖渠段典型斷面進行計算分析。該典型斷面屬軟巖丘陵區(qū)，地面高程91～134 m，區(qū)內沖溝發(fā)育。該段總干渠為深挖方段，挖深約25 m，設3級馬道，其中三級馬道高程為112.59 m、二級馬道高程為106.59 m、一級馬道高程為100.49 m，渠底高程為91.54 m，一級～三級馬道寬度分別為5、8、2 m?？偢汕O計一級～四級邊坡坡比分別為1∶2.5、1∶2.5、1∶2、1∶1.75，設計水位高于渠底7 m。選擇如下計算工況：①設計工況，計算內坡，渠道內設計水深，考慮渠坡內地下水滲流;②校核工況Ⅰ，計算內坡，施工期，渠內無水，考慮渠坡內地下水滲流;③校核工況Ⅱ，計算內坡，設計工況+地震。典型斷面土（巖）c、φ、濕密度、飽和密度見表2，計算得到設計工況、校核工況Ⅰ、校核工況Ⅱ邊坡設計安全系數(shù)分別為1.585、1.582、1.217。軟巖丘陵區(qū)典型斷面邊坡穩(wěn)定計算不同工況時最危險圓弧滑裂面見圖2。

不同工況時最危險圓弧滑裂面

該典型斷面為弱膨脹土（巖）渠段，抗剪強度指標均采用峰值強度，由計算結果可知，設計邊坡滿足規(guī)范要求，該段的處理措施為對膨脹土（巖）出露段進行換填處理，自通水以來該段邊坡穩(wěn)定，未出現(xiàn)影響正常通水的情況。

4.2 存在明顯裂隙面膨脹土（巖）渠坡穩(wěn)定分析

某渠段典型斷面為強膨脹巖渠段，地貌單元為崗地，地形起伏較大，挖深為15～28 m，邊坡主要為黃土狀重粉質壤土（Qalpl3）、重粉質壤土（Qalpl2）及黏土巖（N1L）。基面開挖后，對該段施工邊坡及探槽進行了地質編錄和取樣復核工作，地質情況顯示渠底板主要位于黏土巖（強膨脹性）中，存在深層裂隙結構面，且?guī)r層傾角傾向渠內。兩級計算邊坡坡比分別為1∶3.25和1∶3，一級、二級馬道寬分別為5、2 m。在考慮深層裂隙的情況下，當滑動面呈非圓弧形、不規(guī)則任意滑裂面時，邊坡宜采用摩根斯坦-普拉斯法（折線法）進行抗滑穩(wěn)定計算，裂隙面摩擦角為10°，黏聚力為8.5 kPa。計算工況如下：①設計工況，計算內坡，渠道內設計水深，考慮渠坡內地下水滲流;②校核工況，計算內坡，施工期，渠內無水，考慮渠坡內地下水滲流?？紤]深層裂隙情況時典型斷面土（巖）c、φ、濕密度、飽和密度見表3，計算得到設計工況、校核工況邊坡設計安全系數(shù)分別為1.83、0.99?？紤]深層裂隙情況時典型斷面邊坡穩(wěn)定計算成果見圖3。

該典型斷面為強膨脹巖渠段，抗剪強度指標浸潤線以上采用分帶法取峰值強度的折減，浸潤線以下取殘余強度。計算的校核工況中邊坡安全系數(shù)不滿足規(guī)范要求，因此在考慮裂隙面存在的情況下，對渠道邊坡進行加固很有必要。經分析，該段采用抗滑樁進行渠坡加固。根據(jù)滑裂面處抗剪強度指標，計算滑坡推力，接著根據(jù)地形、地質情況確定抗滑樁設計位置，然后擬定抗滑樁樁長、樁徑、樁間距等設計參數(shù)。該段采用抗滑樁進行深層處理后，對膨脹巖出露段表層進行換填處理，自通水以來該段邊坡穩(wěn)定性方面未出現(xiàn)影響正常通水的情況。

5 結 論

膨脹土（巖）的邊坡穩(wěn)定分析是膨脹土（巖）邊坡設計中一個重點，同時也是難點，目前對膨脹土（巖）渠段沒有公認的行之有效的邊坡穩(wěn)定分析方法。通過對南水北調中線工程部分膨脹土（巖）渠段采用的邊坡穩(wěn)定分析方法進行研究，提出了膨脹土（巖）渠段不同破壞模式下相應的邊坡穩(wěn)定分析方法，并給出了邊坡穩(wěn)定分析時參數(shù)取值的經驗方法。

參考文獻：

[1] 長江水利委員會長江科學院，長江水利委員會長江勘測規(guī)劃設計研究院，南水北調中線干線工程建設管理局，等.膨脹土（巖）渠坡破壞機理及分析方法研究：國家“十一五”科技支撐課題-膨脹土地段渠道破壞機理及處理技術研究-專題二研究報告[R].鄭州：河南省水利勘測設計研究有限公司，2010：2.

[2] 秦祿生，鄭健龍.膨脹土路塹邊坡雨季表層滑坍穩(wěn)定性分析[J].長沙交通學院學報，1999，15（1）：76-79.

[3] 朱伯芳.有限單元法原理與應用[M].2版.北京：中國水利水電出版社，1998：252.

[4] FREDLUND D G， MORGENSTERN N R， WIDGER R A. The Shear Strength of Unsaturated Soils[J]. Canadian Geotechnical Journal，1978，15（3）：313-312.

[5] 鄭穎人，趙尚毅，張魯渝.用有限元強度折減法進行邊坡穩(wěn)定分析[J].中國工程科學，2002，4（10）：57-61.

[6] 陳善雄，陳守義.考慮降雨的非飽和土邊坡穩(wěn)定性分析方法[J].巖土力學，2001，22（4）：447-450.

[7] 龔文惠，王元漢，鄭俊杰.浸水條件下膨脹土路塹邊坡的穩(wěn)定性分析[J].公路，2007（10）：96-100.

[8] 高潤德，彭良泉，王釗.雨水入滲作用下非飽和土邊坡的穩(wěn)定性分析[J].人民長江，2001，32（11）：25-27.

[9] 秦祿生，鄭健龍.膨脹土路基邊坡雨季失穩(wěn)破壞機理的應力應變分析[J].中國公路學報，2001，14（1）：25-30.

[10] 單儒嬌.極限平衡法和有限元法分析排土場邊坡穩(wěn)定性[D].沈陽：遼寧工程技術大學，2009：5.

【責任編輯 張華巖】

收稿日期：2019-06-28

作者簡介：馬慧敏（1981—），女，河南澠池人，高級工程師，主要從事水利水電工程研究工作

E-mail：mahuiminyi@126.com