摘" 要：為分析土石壩在滲流及變形長(zhǎng)期作用下壩基接觸帶性能及變化規(guī)律，基于水力耦合的分析方法，以某水庫(kù)工程不同時(shí)期得到的壩基和壩基接觸帶滲透系數(shù)，建立土石壩壩基接觸帶滲透性能與力學(xué)參數(shù)演化模型，把握壩體及壩基材料物理力學(xué)特性、滲透特性的演化規(guī)律，對(duì)水庫(kù)長(zhǎng)期處于正常蓄水位情況下的滲流場(chǎng)及應(yīng)力場(chǎng)進(jìn)行分析。研究表明，由于水流的長(zhǎng)期滲流作用，導(dǎo)致壩基接觸帶中的土壤流失，滲透系數(shù)增大，進(jìn)而導(dǎo)致通過(guò)壩坡面、兩岸山體的滲漏量增大?？紤]水力耦合對(duì)壩體變形穩(wěn)定的影響，研究發(fā)現(xiàn)，在水庫(kù)運(yùn)行期間，壩基接觸帶的強(qiáng)度是不斷降低的，接觸帶的屈服區(qū)增加，變形量不斷增大，安全系數(shù)隨運(yùn)行時(shí)間的增加不斷降低。

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

中圖分類(lèi)號(hào)：TU46" " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：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

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

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

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

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

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

根據(jù)A.F.Gangi于1978年提出的數(shù)學(xué)模型

式中：K為滲透系數(shù)；σ為有效應(yīng)力；C0為系數(shù)，一般取2～4；k0為表層或松散狀態(tài)下材料的滲透系數(shù)，可取材料的室內(nèi)試驗(yàn)值；Ek為塊石的體積模量。對(duì)于碾壓堆石結(jié)構(gòu)，不同部位的σ都是變化的，而應(yīng)力的變化將同時(shí)導(dǎo)致Ek和K的變化。

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

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

在小變形假設(shè)條件下，根據(jù)接觸帶顆粒的質(zhì)量守恒方程，可將堆石體的孔隙率表示為體積應(yīng)變的負(fù)指數(shù)函數(shù)

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

式中：β為修正系數(shù)。當(dāng)忽略接觸帶顆粒的變形時(shí)，β=1。

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

1.2" 壩基接觸帶水力耦合計(jì)算過(guò)程

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

2" 研究區(qū)域與資料分析

2.1" 工程概況

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

2.2" 壩基接觸帶性能演化

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

根據(jù)該水庫(kù)以往的現(xiàn)場(chǎng)地質(zhì)勘測(cè)工作，以有限的資料進(jìn)行相應(yīng)的分析歸納，近似采用一次函數(shù)對(duì)材料的滲透系數(shù)演化規(guī)律進(jìn)行擬合，壩基接觸帶的滲透系數(shù)隨時(shí)間的變化規(guī)律曲線及相應(yīng)擬合曲線如圖2所示。

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

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

式中：k為滲透系數(shù)；t為年份。

2.2.2" 壩基接觸帶力學(xué)參數(shù)演化

整理該水庫(kù)大壩歷年地質(zhì)勘測(cè)資料后，發(fā)現(xiàn)壩基接觸帶材料的黏聚力與摩擦角均隨時(shí)間增加而降低，但因地質(zhì)勘測(cè)資料較少，難以看出規(guī)律。但借鑒均質(zhì)土壩黏聚力與摩擦角隨時(shí)間變化規(guī)律并將其歸納后，可以認(rèn)為壩基接觸帶材料力學(xué)性能將隨時(shí)間推演至某一收斂值。

由于壩基接觸帶土體材料長(zhǎng)期處于浸潤(rùn)線以下，水的流動(dòng)帶走部分土體顆粒，導(dǎo)致土體內(nèi)部顆粒重排，壩基接觸帶的土體結(jié)構(gòu)會(huì)發(fā)生一定程度變化。

根據(jù)歷史監(jiān)測(cè)資料，反算相關(guān)折減系數(shù)，改變折減系數(shù)逐次計(jì)算，直至當(dāng)折減系數(shù)到達(dá)某一值時(shí)，當(dāng)年測(cè)點(diǎn)沉降計(jì)算值將極其貼近于監(jiān)測(cè)值。將反演得到的折減系數(shù)隨時(shí)間變化進(jìn)行函數(shù)擬合，擬合曲線如圖3所示。

擬合表達(dá)式

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

2.3" 滲流及穩(wěn)定性分析

2.3.1" 有限元模型

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

2.3.2" 計(jì)算參數(shù)及計(jì)算條件

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

該水庫(kù)模型邊界條件設(shè)置為：①上游壩面前庫(kù)水淹沒(méi)區(qū)、上游水位以下的上游壩面節(jié)點(diǎn)及上游水位下兩岸山體區(qū)域取上游定水頭邊界；②下游土石壩壩體、排水棱體及兩岸山體處等在下游水位以下的表面節(jié)點(diǎn)取下游定水頭邊界；③模型頂部未被水浸潤(rùn)區(qū)域設(shè)置為溢出邊界；④模型底部邊界取隔水邊界。

本文中水庫(kù)滲流場(chǎng)計(jì)算均取其正常運(yùn)行工況時(shí)水位，上游運(yùn)行水位為157.8 m，下游無(wú)水。根據(jù)滲流場(chǎng)分析結(jié)果，針對(duì)正常運(yùn)行工況下開(kāi)展壩區(qū)典型工況滲流場(chǎng)特性分析。為了更直觀地展示壩區(qū)滲流場(chǎng)的規(guī)律，選取了壩區(qū)中心剖面進(jìn)行滲流計(jì)算成果的展示。計(jì)算參數(shù)見(jiàn)表1。

2.3.3" 滲流計(jì)算成果分析

1）原始滲流場(chǎng)規(guī)律。

采用表1所示滲透系數(shù)對(duì)水庫(kù)壩體及壩基三維整體滲流場(chǎng)進(jìn)行計(jì)算，計(jì)算結(jié)果如圖5所示。

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

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

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

采用公式（4）對(duì)壩基接觸帶進(jìn)行演化分析，相關(guān)計(jì)算結(jié)果見(jiàn)表2。

從表中可以看出，隨著運(yùn)行時(shí)間的增加，通過(guò)壩坡、山體左右岸的滲漏量都有一定程度的增加，這是因?yàn)樵谒畮?kù)運(yùn)行的過(guò)程中，由于滲流的作用，壩基接觸帶內(nèi)的土壤細(xì)顆粒隨水流緩慢流失，導(dǎo)致滲透系數(shù)增大，進(jìn)而引起滲漏量的增大。圖6是水庫(kù)運(yùn)行不同時(shí)間階段典型剖面的滲透坡降矢量圖，從圖中可以看出，隨著水庫(kù)運(yùn)行時(shí)長(zhǎng)的增加，最大滲透坡降值不斷增加，漲幅約1.2%。

2.4" 變形穩(wěn)定計(jì)算成果分析

2.4.1" 計(jì)算參數(shù)及計(jì)算條件

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

2.4.2" 應(yīng)力計(jì)算成果分析

圖7為水庫(kù)長(zhǎng)期處于正常蓄水位條件下的不同階段下的塑性區(qū)分布圖，隨著運(yùn)行時(shí)間的增加，壩基接觸帶上的塑性屈服區(qū)由中心向上下游不斷增加的，這是由于壩基接觸帶的細(xì)部顆粒在滲流的作用下不斷流失，導(dǎo)致壩基接觸帶的強(qiáng)度不斷下降。

圖8為水庫(kù)長(zhǎng)期處于正常蓄水位條件下的不同階段下的最大剪切應(yīng)變?cè)隽繄D。從水庫(kù)開(kāi)始運(yùn)行至水庫(kù)運(yùn)行15年，通過(guò)FLAC3D自帶強(qiáng)度折減法計(jì)算得出水庫(kù)壩體的抗滑穩(wěn)定安全系數(shù)由1.98降至1.793，降低幅度約為9.4%。

2.5" 小結(jié)

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

1）水庫(kù)長(zhǎng)期處于正常蓄水位情況，由于水流的長(zhǎng)期滲流作用，導(dǎo)致壩基接觸帶中的土壤流失，滲透系數(shù)增大，進(jìn)而導(dǎo)致通過(guò)壩坡面、兩岸山體的滲漏量增大。

2）考慮水力耦合對(duì)壩體變形穩(wěn)定的影響，研究發(fā)現(xiàn)，在水庫(kù)運(yùn)行期間，壩基接觸帶的強(qiáng)度是不斷降低的，接觸帶的屈服區(qū)增加，變形量不斷增大，最大變形量為6.966 mm，心墻以下接觸帶出現(xiàn)少許應(yīng)力集中范圍，最大壓應(yīng)力為1.344 MPa。

3）計(jì)算所得壩體抗滑穩(wěn)定安全系數(shù)隨運(yùn)行時(shí)間的增加不斷降低，由1.98降至1.793，降低幅度約為9.4%。

3" 結(jié)論

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

在實(shí)際工程中，通常采用的垂直防滲措施有混凝土防滲墻、高壓噴射灌漿、垂直鋪塑、深層攪拌樁防滲墻、充填灌漿和帷幕灌漿等。加固處理設(shè)計(jì)時(shí)應(yīng)考慮不同的壩基接觸帶滲漏嚴(yán)重程度、接觸帶的破壞形式、當(dāng)?shù)夭牧虾褪┕l件等因素，采取一種或多種處理方法相結(jié)合的防滲處理措施。

參考文獻(xiàn)：

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

[2] 田東方，劉德富，王世梅，等.土質(zhì)邊坡非飽和滲流場(chǎng)與應(yīng)力場(chǎng)耦合數(shù)值分析[J].巖土力學(xué)，2009，30（3）：810-814.

[3] 劉祎，蔡國(guó)慶，李艦，等.非飽和土熱-水-力全耦合本構(gòu)模型及其驗(yàn)證[J].巖土工程學(xué)報(bào)，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] 陳益峰，周創(chuàng)兵，毛新瑩，等.水布埡地下廠房圍巖滲控效應(yīng)數(shù)值模擬與評(píng)價(jià)[J].巖石力學(xué)與工程學(xué)報(bào)，2010，29（2）：308-318.

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