考慮土石壩壩基接觸帶性能演化的滲流穩定分析

摘" 要：為分析土石壩在滲流及變形長期作用下壩基接觸帶性能及變化規律，基于水力耦合的分析方法，以某水庫工程不同時期得到的壩基和壩基接觸帶滲透系數，建立土石壩壩基接觸帶滲透性能與力學參數演化模型，把握壩體及壩基材料物理力學特性、滲透特性的演化規律，對水庫長期處于正常蓄水位情況下的滲流場及應力場進行分析。研究表明，由于水流的長期滲流作用，導致壩基接觸帶中的土壤流失，滲透系數增大，進而導致通過壩坡面、兩岸山體的滲漏量增大。考慮水力耦合對壩體變形穩定的影響，研究發現，在水庫運行期間，壩基接觸帶的強度是不斷降低的，接觸帶的屈服區增加，變形量不斷增大，安全系數隨運行時間的增加不斷降低。

關鍵詞：土石壩；壩基接觸帶；滲透性能；演化模型；水力耦合

中圖分類號：TU46" " " 文獻標志碼：A" " " " " 文章編號：2095-2945（2024）17-0076-06

Abstract： In order to analyze the performance and variation law of dam foundation contact zone of earth-rockfill dam under the long-term action of seepage and deformation， based on the analysis method of hydraulic coupling， based on the permeability coefficient of dam foundation and dam foundation contact zone obtained in different periods of a reservoir project， the evolution model of permeability and mechanical parameters of earth-rockfill dam foundation contact zone is established， and the evolution law of physical and mechanical properties and permeability characteristics of dam body and dam foundation materials are grasped. The seepage field and stress field of the reservoir under the condition of normal water level for a long time are analyzed. The study shows that due to the long-term seepage of water flow， the soil in the contact zone of the dam foundation is lost， and the permeability coefficient increases， which leads to the increase of leakage through the dam slope and the mountains on both sides. Considering the influence of hydraulic coupling on the deformation and stability of the dam， it is found that during the operation of the reservoir， the strength of the contact zone of the dam foundation decreases， the yield zone of the contact zone increases， the deformation increases， and the safety factor decreases with the increase of running time.

Keywords： Earth-rockfill dam; dam foundation contact zone; permeability; evolution model; hydraulic coupling

土石壩壩基接觸帶指自然基礎與壩體填土接觸部分，重點指土石壩壩體防滲心墻與壩基防滲帷幕結合部位，壩體防滲體與壩基防滲體由于組成成分、結構、力學特性等方面差異，兩者在長期運行過程中（如庫水位漲落）時常出現變形不協調、滲漏情況各異等現象，造成土石壩滲漏破壞[1]。為了保證水庫大壩的安全運行和人民生命財產安全，需要從宏觀及微觀角度對土石壩壩基接觸帶進行破壞機理研究。隨著計算技術的發展，相關的分析軟件已逐漸完善，如GEO-STUDIO、ABAQUS、ANSYS等，國內也有部分學者針對水力耦合進行了研究，田東方等[2]考慮了坡面徑流的影響，依據試驗結果提出了坡面徑流-非飽和土水力耦合分析方法，并編寫了對應的有限元程序，但未考慮到滲流場對非飽和土體強度參數的影響。劉祎等[3]以平均土骨架應力、修正吸力和溫度為應力變量，以比體積和飽和度為應變變量，建立了非飽和土的熱-水-力全耦合本構模型。本文通過建立有限元法計算滲流與有限差分法計算應力場的雙場間接耦合方法，可對特定計算參數進行公式修正，提高交替迭代過程中程序代碼的結果準確性，彌補了FLAC3D軟件對滲流場計算時，計算效率低、耗時長的缺點。

1" 土石壩壩基接觸帶水力耦合分析方法研究

1.1" 壩基接觸帶水力耦合分析方法

滲流與應力（或變形）的耦合現象是客觀存在的。壩基接觸帶受到上部填筑荷載的作用必然會被壓實，孔隙度減小，而塊體之間的擠壓、研磨，可能使塊石邊緣壓裂甚至壓碎。實驗研究表明，接觸帶的滲透特性及變形特性決定于其成因、物質組成、結構特征及厚度等。由于壩基接觸帶上、下2層土體材料性質差異較大，具備發生管涌與流土的條件，需進行必要的防滲處理。

由于壩基接觸帶結構松散、不連續，在水平和垂直方向有較大變化、成因類型復雜、物理力學性質呈較大不均勻性的特點，并且壩基接觸帶在上部壩料的填筑以及大壩的蓄水過程中，其應力狀態將發生變化，粗顆粒發生擠壓、研磨，甚至出現尖端破裂現象，從而導致接觸帶土體材料孔隙率和滲透特性發生變化。因此，接觸帶的滲透特性不僅與顆粒級配和初始碾壓狀態密切相關，還受到應力狀態的影響[4-5]。接觸帶的滲透特性在本質上受控于其內部連通的孔隙率或密實度，目前還沒有提出關于壩基接觸帶一般性的耦合機制。但是從壩基接觸帶的顆粒級配上來說，壩基接觸帶主要還是以粗顆粒（粒徑）為主，可以借鑒類似堆石體耦合機制來探討接觸帶的水力耦合機制。如陳益峰等[6]通過對水布埡面板堆石壩的研究，建立了高面板堆石壩非線性變形與非穩定滲流的水力耦合模型。

根據A.F.Gangi于1978年提出的數學模型

式中：K為滲透系數；σ為有效應力；C0為系數，一般取2～4；k0為表層或松散狀態下材料的滲透系數，可取材料的室內試驗值；Ek為塊石的體積模量。對于碾壓堆石結構，不同部位的σ都是變化的，而應力的變化將同時導致Ek和K的變化。

借助于多孔介質的Kozeny-Carman 方程，則接觸帶的滲透系數k與孔隙率n之間的關系可表示為

式中：k0為接觸帶初始的滲透系數；n0為初始孔隙率。

在小變形假設條件下，根據接觸帶顆粒的質量守恒方程，可將堆石體的孔隙率表示為體積應變的負指數函數

n=1-（1-n）exp（-βεv）， （3）

式中：β為修正系數。當忽略接觸帶顆粒的變形時，β=1。

但為了考慮顆粒變形及尖端破裂等因素對孔隙率的影響，可根據接觸帶的級配曲線，采用2009年周偉等建立的隨機散粒體不連續變形模型，通過數值模擬給出接觸帶在不同圍壓條件下的偏應力-軸向應變和孔隙率-體積應變等關系曲線，從而得到接觸帶與圍壓的最佳擬合公式。

1.2" 壩基接觸帶水力耦合計算過程

在雙場耦合分析方面，目前大致有2類方法：一是將兩場分開計算，然后通過兩場的交叉迭代達到耦合的目的，即間接耦合法；二是建立以應力場和滲流場為未知量的數學模型，通過直接求得解答達到完全耦合的目的，即直接耦合法[7]。本文采用迭代法即間接耦合法，滲流場和應力場耦合分析計算主要內容為：分別進行初始的滲流場與應力場，將滲透力與應力場疊加，確定新的滲流場與應力場，循環迭代直至滿足收斂條件，具體流程如圖1所示。

2" 研究區域與資料分析

2.1" 工程概況

本項目以湖南省某水庫為例，開展對土石壩壩基接觸帶的一系列研究。該水庫大壩為土壩，最大壩高為15 m，壩頂高程為159.5 m，水庫正常蓄水位為157.8 m，死水位為145.5 m，是一座小（1）型水利工程。

2.2" 壩基接觸帶性能演化

2.2.1" 壩基接觸帶滲透性能演化

根據該水庫以往的現場地質勘測工作，以有限的資料進行相應的分析歸納，近似采用一次函數對材料的滲透系數演化規律進行擬合，壩基接觸帶的滲透系數隨時間的變化規律曲線及相應擬合曲線如圖2所示。

壩基接觸帶滲透系數隨時間變化的擬合曲線方程為

k=1.179×10-5 t-0.023， （4）

式中：k為滲透系數；t為年份。

2.2.2" 壩基接觸帶力學參數演化

整理該水庫大壩歷年地質勘測資料后，發現壩基接觸帶材料的黏聚力與摩擦角均隨時間增加而降低，但因地質勘測資料較少，難以看出規律。但借鑒均質土壩黏聚力與摩擦角隨時間變化規律并將其歸納后，可以認為壩基接觸帶材料力學性能將隨時間推演至某一收斂值。

由于壩基接觸帶土體材料長期處于浸潤線以下，水的流動帶走部分土體顆粒，導致土體內部顆粒重排，壩基接觸帶的土體結構會發生一定程度變化。

根據歷史監測資料，反算相關折減系數，改變折減系數逐次計算，直至當折減系數到達某一值時，當年測點沉降計算值將極其貼近于監測值。將反演得到的折減系數隨時間變化進行函數擬合，擬合曲線如圖3所示。

擬合表達式

k=1.00-0.012 57（t-2 000），2 000?t 。 （5）

2.3" 滲流及穩定性分析

2.3.1" 有限元模型

通過收集該水庫主壩壩區的地質地貌、水文信息等，依據樞紐布置方案，建立包括左岸山體，右岸條形山及主壩部分的地形地貌創建三維可視化模型。其中，以主壩壩軸線為x軸，垂直壩軸線方向為y軸，高程作為z軸。有限元模型共有單元177 241個，節點97 592個。有限元模型范圍如下：①取大壩上游邊界距壩軸線垂直距離約87 m，大壩下游側邊界距壩軸線垂直距離約84 m；②模型最低高程取正常蓄水位157.8 m以下約48 m，模型底部高程約為110 m；③整個計算模型上、下游邊界相距約244 m，左邊界相距大壩中心線121 m，右邊界相距123 m。主要模型如圖4所示。

2.3.2" 計算參數及計算條件

本文選取水庫工程三維滲流分析模型以大壩和壩基滲控設計方案為主要參考依據，在三維有限元模型中，滲流計算模型范圍：XY 245.48 m×172.09 m。模型中水工建筑物包括主壩壩體、排水棱體、防滲墻及防滲帷幕。

該水庫模型邊界條件設置為：①上游壩面前庫水淹沒區、上游水位以下的上游壩面節點及上游水位下兩岸山體區域取上游定水頭邊界；②下游土石壩壩體、排水棱體及兩岸山體處等在下游水位以下的表面節點取下游定水頭邊界；③模型頂部未被水浸潤區域設置為溢出邊界；④模型底部邊界取隔水邊界。

本文中水庫滲流場計算均取其正常運行工況時水位，上游運行水位為157.8 m，下游無水。根據滲流場分析結果，針對正常運行工況下開展壩區典型工況滲流場特性分析。為了更直觀地展示壩區滲流場的規律，選取了壩區中心剖面進行滲流計算成果的展示。計算參數見表1。

2.3.3" 滲流計算成果分析

1）原始滲流場規律。

采用表1所示滲透系數對水庫壩體及壩基三維整體滲流場進行計算，計算結果如圖5所示。

從圖5中可以看出，在主壩的典型剖面段，自由面呈逐步降低趨勢，自由面在壩體自上游向下游有一次明顯的急劇下降段，由于壩體存在沖抓回填與灌漿帷幕，自由面在帷幕處急劇下降，說明帷幕防滲效果十分顯著，其他部分相對平緩，在下游排水棱體處水位下降有所加快，排水棱體有效地降低了浸潤線。下游壩坡在高程145.3 m左右出現滲流溢出點。從圖中還可以看出，自由面在壩體內部的分布特征完全滿足光滑連續和單調下降這2個基本幾何性質，因而滲流計算成果在理論上是正確的。

在水庫正常蓄水位作用下，最大滲透坡降出現在防滲帷幕處，達到6.48；在排水棱體處滲透坡降達到了0.72左右；大壩內部其他地方及地層內的滲透坡降均小于0.4。

2）壩基接觸帶滲透性能演化計算分析。

采用公式（4）對壩基接觸帶進行演化分析，相關計算結果見表2。

從表中可以看出，隨著運行時間的增加，通過壩坡、山體左右岸的滲漏量都有一定程度的增加，這是因為在水庫運行的過程中，由于滲流的作用，壩基接觸帶內的土壤細顆粒隨水流緩慢流失，導致滲透系數增大，進而引起滲漏量的增大。圖6是水庫運行不同時間階段典型剖面的滲透坡降矢量圖，從圖中可以看出，隨著水庫運行時長的增加，最大滲透坡降值不斷增加，漲幅約1.2%。

2.4" 變形穩定計算成果分析

2.4.1" 計算參數及計算條件

在本次計算過程中，選擇了FLAC3D中內置的Mohr-Coulomb模型。開展計算時，因模型四周均位于山體中，故在四周邊界施加法向約束，模型底部在Z=110 m施加固定約束。值得注意的是，本文應力變形計算以水庫建成蓄水后為初始狀態，計算的大壩變形均為蓄水后的累積變形。

2.4.2" 應力計算成果分析

圖7為水庫長期處于正常蓄水位條件下的不同階段下的塑性區分布圖，隨著運行時間的增加，壩基接觸帶上的塑性屈服區由中心向上下游不斷增加的，這是由于壩基接觸帶的細部顆粒在滲流的作用下不斷流失，導致壩基接觸帶的強度不斷下降。

圖8為水庫長期處于正常蓄水位條件下的不同階段下的最大剪切應變增量圖。從水庫開始運行至水庫運行15年，通過FLAC3D自帶強度折減法計算得出水庫壩體的抗滑穩定安全系數由1.98降至1.793，降低幅度約為9.4%。

2.5" 小結

通過對土石壩壩基接觸帶漸進破壞機理與安全控制措施進行研究，基于水力耦合分析方法，以某水庫工程為例，建立了壩基接觸帶性能演化模型。主要結論如下：

1）水庫長期處于正常蓄水位情況，由于水流的長期滲流作用，導致壩基接觸帶中的土壤流失，滲透系數增大，進而導致通過壩坡面、兩岸山體的滲漏量增大。

2）考慮水力耦合對壩體變形穩定的影響，研究發現，在水庫運行期間，壩基接觸帶的強度是不斷降低的，接觸帶的屈服區增加，變形量不斷增大，最大變形量為6.966 mm，心墻以下接觸帶出現少許應力集中范圍，最大壓應力為1.344 MPa。

3）計算所得壩體抗滑穩定安全系數隨運行時間的增加不斷降低，由1.98降至1.793，降低幅度約為9.4%。

3" 結論

目前針對土石壩壩基接觸帶安全控制措施較多，主要以各類防滲墻、灌漿防滲及土工合成材料防滲等垂直防滲為主。采取此類防滲措施，一般不需要放空水庫。反之，如果采取水平防滲措施，則必須放空水庫，才能徹底進行。而水庫長期蓄水后，總會有些淤積，給水平防滲處理帶來一定的困難；在保持壩基滲透穩定和截滲方面，水平防滲也不如垂直防滲徹底。

在實際工程中，通常采用的垂直防滲措施有混凝土防滲墻、高壓噴射灌漿、垂直鋪塑、深層攪拌樁防滲墻、充填灌漿和帷幕灌漿等。加固處理設計時應考慮不同的壩基接觸帶滲漏嚴重程度、接觸帶的破壞形式、當地材料和施工條件等因素，采取一種或多種處理方法相結合的防滲處理措施。

參考文獻：

[1] 劉妍.灌漿技術在壩基接觸帶處理中的應用[J].河南科技，2020（1）：99-101.

[2] 田東方，劉德富，王世梅，等.土質邊坡非飽和滲流場與應力場耦合數值分析[J].巖土力學，2009，30（3）：810-814.

[3] 劉祎，蔡國慶，李艦，等.非飽和土熱-水-力全耦合本構模型及其驗證[J].巖土工程學報，2021，43（3）：547-555.

[4] GALLIPOLI D， WHEELER S J， KARSTUNEN K. Modelling the variation of degree of saturation in a deformable unsaturated soil[J].Géotechnique，2003，53（1）：105-112.

[5] GALLAGE C， UCHIMURA T. Effects of dry density and grain size distribution on soil-water characteristic curves of sandy soils[J].Soils and Foundations，2010，50（1）：161-172.

[6] 陳益峰，周創兵，毛新瑩，等.水布埡地下廠房圍巖滲控效應數值模擬與評價[J].巖石力學與工程學報，2010，29（2）：308-318.

[7] 王凱，鄒斌，楊麗虹，等.幸福水庫大壩滲流安全分析及評價[J].人民珠江，2019，40（10）：110-115.