摘 要：混凝土面板和接縫止水是面板堆石壩壩體主要防滲結(jié)構(gòu)，面板開裂或止水失效會(huì)顯著影響裂縫或失效部位局部范圍滲流場分布，影響大壩正常服役性能。為探究面板局部開裂和止水結(jié)構(gòu)失效引起的滲流安全問題，設(shè)計(jì)不同面板開裂和止水失效位置，建立鉛山抽蓄電站混凝土面板壩三維滲流有限元分析模型，研究其滲流位勢分布、滲透坡降和滲透流量演化規(guī)律，探明面板開裂和止水失效的影響，并開展臨時(shí)度汛條件下墊層料滲透系數(shù)敏感特性分析。結(jié)果表明：隨著面板止水結(jié)構(gòu)失效、開裂位置降低，通過面板滲透流量增大。臨時(shí)度汛條件下，隨墊層滲透系數(shù)增大，墊層最大平均滲透坡降減小。墊層滲透系數(shù)從1. 0×10-5 m/s增大至5. 0×10-5 m/s時(shí)，墊層最大平均滲透坡降減小了1. 94。建議墊層滲透系數(shù)宜大于5. 0×10-5 m/s，可提高臨時(shí)度汛條件下墊層滲透穩(wěn)定安全。研究成果可為類似抽蓄電站面板堆石壩防滲體系設(shè)計(jì)和施工提供理論依據(jù)。

關(guān)鍵詞：混凝土面板壩；面板失效；臨時(shí)度汛；滲流分析；鉛山抽蓄電站

中圖分類號：TV223. 4 文獻(xiàn)標(biāo)識碼：A 文章編號：1001-9235（2025）03-0039-07

Three-Dimensional Seepage Characteristics of Yanshan Concrete Slab Dam under Slab"Failure and Temporary Flood Conditions

YUE Xiao1， CHEN Dewen1， HUANG Xinyi2， GAN Lei2*， SHEN Zhenzhong2

Abstract： Concrete slabs and water stop joints are the main seepage control structures of the slab rockfill dams. Slab cracking or water stop failure will significantly affect the distribution of seepage fields in the local area of the cracks or the failure site， which will threaten the normal service performance of the dam. To investigate the seepage safety problems caused by local slab cracking and failure of water stop structures， different positions of slab cracking and water stop failure were designed， and a three-dimensional seepage finite element analysis model of the concrete slab dam of Yanshan pumping and storage power station was established. The evolution law of seepage potential distribution， seepage gradient， and seepage flow rate was studied. The effects of slab cracking and water stop failure were explored， and the sensitivity of the cushion material's permeability coefficient under temporary flood conditions was analyzed. The results show that the seepage flow rate through the slab increases as the water stop structure fails， and the cracking position decreases. Under temporary flood conditions， the maximum average seepage gradient of the cushion decreases with the increase in permeability coefficient. When the permeability coefficient of the cushion increases from 1. 0 × 10?5 m/s to 5. 0 × 10?5 m/s， the maximum average seepage gradient of the cushion decreases by 1. 94. It is suggested that the permeability coefficient of the"cushion should be greater than 5. 0 × 10?5 m/s， which can improve the stability and safety of cushion seepage under temporary flood conditions. The research results can provide a theoretical basis for the design and construction of a slab rockfill dam anti-seepage system of similar pumping and storage power stations.

Keywords： concrete slab dam; slab failure; temporary flood; seepage analysis; Yanshan pumping and storage power station

混凝土面板堆石壩是抽水蓄能電站常采用的優(yōu)勢壩型，具有施工便捷、適應(yīng)性良好、造價(jià)低等特點(diǎn)。壩區(qū)滲流安全是保障大壩安全運(yùn)行和電站效益發(fā)揮的前提。混凝土面板和止水結(jié)構(gòu)是面板堆石壩的重要防滲結(jié)構(gòu)，一旦面板局部開裂或者垂直縫止水失效，會(huì)威脅大壩整體安全和效益［1］。

混凝土面板不僅可防止?jié)B漏，還可傳遞水壓力。面板堆石壩面板開裂主要是由壩體不均勻變形、溫度應(yīng)力、干縮應(yīng)力、施工原因等原因引起［2-5］。姜媛媛等［6］發(fā)現(xiàn)當(dāng)面板局部發(fā)生水平或豎向裂縫時(shí)，裂縫附近墊層過渡料區(qū)中會(huì)出現(xiàn)集中水頭損失和較大滲透坡降。林起明等［7］探究了大壩在各種裂縫情況下滲流特性，認(rèn)為面板任意位置開裂，其壩體填料滲透坡降會(huì)顯著提高，易產(chǎn)生滲透破壞，且滲流量急劇增加。張曉強(qiáng)等［8］以班多面板堆石壩為例，分析比較了面板接縫止水失效位置、失效長度和失效寬度對壩體壩基及止水失效局部區(qū)域滲流場特性的影響，指出壩體內(nèi)浸潤面、滲透流量隨接縫失效尺寸的增大而提高。溫鵬業(yè)［9］認(rèn)為相比于裂縫止水結(jié)構(gòu)失效引起的壩體滲流場變化，面板滲透系數(shù)對壩體滲流場影響更大一些。熊璐等［10］認(rèn)為面板接縫失效后，墊層區(qū)能擋住70%左右的水頭，成為防滲的第二道防線。

面板堆石壩施工周期較長，往往跨多個(gè)汛期，且面板施工順序靠后，在水位高、變化快的汛期常需采用壩體臨時(shí)斷面擋水度汛和墊層直接擋水度汛，以保障汛期面板堆石壩滲透穩(wěn)定安全［11］。因此墊層設(shè)計(jì)至關(guān)重要，不僅對面板起柔性支承作用且起臨時(shí)擋水作用。墊層直接擋水度汛時(shí)，壩體水頭差、坡降、流速隨著墊層區(qū)、過渡區(qū)和主堆石滲透系數(shù)增加而減?。?2］。張鳳財(cái)?shù)龋?3］計(jì)算猴子巖堆石壩防滲面板出現(xiàn)裂縫時(shí)的滲流場，指出墊層滲透系數(shù)對滲流場和滲流量的影響較大，但未考慮墊層滲透系數(shù)。喬蓓等［14］

提出單縫失效時(shí)，垂直縫位置、縫寬、縫長及墊層滲透系數(shù)與各滲流控制要素呈正相關(guān)。Gan等［15］提出了一種新的滲流反演技術(shù)來預(yù)測裂隙或斷層發(fā)育的巖體的滲透系數(shù)。

綜上所述，墊層作為第二道防滲生命線或臨時(shí)度汛時(shí)起直接擋水作用已有所認(rèn)識，但缺乏臨時(shí)度汛條件下墊層滲透系數(shù)變化對壩體滲流場分布影響的定量研究，本文以鉛山抽水蓄能電站混凝土面板壩為研究對象，建立其三維滲流有限元分析模型，開展面板止水結(jié)構(gòu)失效、局部開裂和臨時(shí)度汛條件下墊層滲透系數(shù)敏感性分析，研究面板止水結(jié)構(gòu)失效、局部開裂及不同失效位置引起的滲流影響，確定臨時(shí)度汛條件下墊層料合理的滲透參數(shù)，以期為類似工程提供借鑒。

1 基本理論

本文采用三維有限元滲流分析程序［16-18］，飽和滲流有限元控制微分方程如下：

?s?h?h

式中：h為總水頭；kisj為飽和滲透系數(shù)張量；Ss為彈性貯水率，飽和土體的Ss為一個(gè)常數(shù)；Q為源匯項(xiàng)。

飽和穩(wěn)定-非穩(wěn)定滲流微分方程的定解條件如下［19-20］：

a））初始條件

h（xi，0） =h（xi，t0）， i= 1，2，3（2） b））邊界條件

-kij

（5）式中：ni為邊界面外法線方向余弦；t0為初始時(shí)刻；h1為已知水頭；qn為已知流量；Γ1為已知水頭邊界；Γ2為已知流量邊界；Γ3為飽和逸出面邊界。滲流邊界見圖1。

2 實(shí)例分析

2. 1 計(jì)算模型

鉛山混凝土面板堆石壩位于江西省上饒市鉛山縣，壩頂高程1 048. 00 m，最大壩高為93. 0 m，壩頂長463. 0 m，壩頂寬8 m。防浪墻頂高程1 049. 20 m，壩體上、下游壩坡坡比為1. 0∶1. 4。上水庫正常蓄水位1 044. 00 m，相應(yīng)死水位為1 015. 00 m。采用“控制斷面超單元有限元網(wǎng)格自動(dòng)剖分法”，截取從混凝土面板壩上游坡腳以上703 m至下游坡腳以下314 m，左、右岸邊界分別截取至左壩肩以左459 m、右壩肩以右384 m；頂高程按實(shí)際地形考慮，底高程截至800 m，至新鮮巖體。模型截取范圍見圖2。

模型采用笛卡爾坐標(biāo)系，具體如下：原點(diǎn)取大地坐標(biāo)（x，y）=（567 793. 049，3 095 415. 779）；取X軸為順河流方向，垂直于壩軸線，上游指向下游為正；Y軸為壩軸線方向，右岸指向左岸為正；Z軸為垂直方向，向上為正，與高程一致。

根據(jù)建筑物布置、巖體分層、斷層構(gòu)造以及計(jì)算要求等信息，建立三維有限元計(jì)算模型。取控制斷面15個(gè)。首先形成三維超單元網(wǎng)格，加密細(xì)分后形成三維有限元網(wǎng)格，庫盆開挖及工程竣工后生成的有限元網(wǎng)格結(jié)點(diǎn)總數(shù)為36 136個(gè)，單元總數(shù)為34 632個(gè)。庫盆開挖及工程竣工后計(jì)算模型三維有限元網(wǎng)格見圖3。對河床段部分模型進(jìn)行加密，建立精細(xì)的三維滲流有限元模型。加密細(xì)分部分的有限元網(wǎng)格見圖4，節(jié)點(diǎn)數(shù)為5 940，單元總數(shù)5 640。墊層區(qū)域局部網(wǎng)格進(jìn)行了加密處理，加密后的有限元網(wǎng)格見圖5。

建立模型時(shí)，除特殊結(jié)構(gòu)外，各主要建筑物（或結(jié)構(gòu)）均按實(shí)際尺寸模擬。防滲帷幕僅考慮最深的一排，滲透系數(shù)按實(shí)際考慮。對上水庫壩址區(qū)巖體按不同呂榮值（Lu）進(jìn)行分層模擬，分別為微透水層（＜1 Lu），微弱透水層（1～3 Lu），弱透水層下段（3～5 Lu），弱透水層上段（＞5 Lu）；由于計(jì)算區(qū)域內(nèi)無對計(jì)算影響較大斷層，無需對斷層進(jìn)行模擬。面板垂直止水結(jié)構(gòu)失效或開裂裂縫采用流量等效原則進(jìn)行模擬［21］。

2. 2 計(jì)算參數(shù)和工況

2. 2. 1 反演參數(shù)及邊界條件

結(jié)合工程實(shí)際地質(zhì)情況以及反演分析成果，確定上水庫各巖層滲透系數(shù)見表1，表2為壩體各材料分區(qū)滲透參數(shù)［22-23］。

模型計(jì)算邊界主要有水頭邊界、不透水邊界和出滲邊界。具體如下：已知水頭邊界包括水庫蓄水位以下的給定水頭邊界，以及給定地下水位的地下分水嶺、截取邊界；出滲邊界包括水庫蓄水位以上的左、右岸和其他庫盆山坡的迎水面，所有與大氣接觸的邊界；不透水邊界包括模型底面以及模型四周（上下游及兩側(cè)）截取邊界。

2. 2. 2 計(jì)算工況

為研究面板止水結(jié)構(gòu)不同失效位置、面板不同開裂位置和墊層不同滲透系數(shù)對壩體滲流特性的影響，設(shè)計(jì)了正常運(yùn)行期不同高程位置處止水部分失效、面板部分開裂和臨時(shí)度汛條件下墊層滲透系數(shù)敏感性分析3組試驗(yàn)方案。D-1方案模擬的是面板垂直止縫水部分失效工況；D-2方案模擬的是面板部分開裂工況；D-3模擬的是墊層滲透系數(shù)敏感性分析工況，表3為計(jì)算工況［9，24］。

2. 3 計(jì)算結(jié)果

2. 3. 1 垂直縫止水失效和面板開裂工況

各工況主要防滲體和重要部位巖體的最大滲透坡降見表4。當(dāng)止水失效和面板開裂時(shí)，面板滲透坡降變化在0. 4%以內(nèi)，基巖的滲透坡降穩(wěn)定在0. 170左右，因而對壩內(nèi)堆石體等各分區(qū)填筑體基本沒有影響，壩體的滲透穩(wěn)定性滿足要求。

各工況通過混凝土面板滲流量見表5。由表5可知，當(dāng)止水失效時(shí)，滲流量變化可增加12. 23%；面板部分開裂時(shí)，滲流量變化可高至40. 89%。與防滲體系完好相比，止水失效、面板部分開裂時(shí)，滲流量明顯增大。因此，由于裂縫的存在，通過混凝土面板滲流量明顯增大，與防滲面板完好時(shí)相比，該變化流量主要由止水失效和面板開裂部位產(chǎn)生。當(dāng)止水失效，且失效位置高程從990 m到1 030 m時(shí)，滲流量變化從12. 23%降低至4. 18%。當(dāng)面板開裂，且失效位置高程從990 m至1 030 m時(shí)，滲流量變化從40. 89%顯著降低至15. 13%。因而，面板止水結(jié)構(gòu)、面板裂縫不同失效位置對滲透流量變化較為敏感，相同失效縫長和失效縫寬時(shí)，面板止水結(jié)構(gòu)、面板縫失效部位作用的水頭越大，通過面板的滲透流量也越大。隨著面板止水結(jié)構(gòu)失效、面板開裂位置降低，通過面板的滲流量增大。

面板垂直縫止水結(jié)構(gòu)失效對滲透坡降的影響主要體現(xiàn)在失效裂縫局部范圍內(nèi)（圖6）。當(dāng)止水失效位于990 m高程時(shí)，失效影響局部滲流場順坡向范圍為9. 970 m，面板法向方向上形成局部飽和區(qū)，飽和區(qū)深度約為1. 232 m；而當(dāng)止水失效處于高程1 030 m時(shí)，失效影響局部滲流場順坡向范圍為6. 783 m，飽和區(qū)深度約為0. 554 m。可見在面板止水結(jié)構(gòu)失效區(qū)下部墊層內(nèi)，隨著面板止水結(jié)構(gòu)失效位置高程降低，飽和區(qū)范圍越大。當(dāng)面板開裂處位于高程990 m時(shí)，裂縫影響局部滲流場順坡向范圍為4. 320 m，面板法向方向上形成局部飽和區(qū)，飽和區(qū)深度約為0. 497 m；而當(dāng)面板開裂處處于高程1 030 m時(shí)，裂縫影響局部滲流場順坡向范圍為1. 435 m，飽和區(qū)深度約為0. 438 m。飽和區(qū)在面板開裂失效區(qū)下部墊層內(nèi)，且隨面板開裂失效位置高程降低，墊層內(nèi)飽和區(qū)范圍越大。

2. 3. 2 臨時(shí)度汛工況

臨時(shí)度汛條件下，墊層不同滲透系數(shù)工況下Y = 637 m剖面位勢分布見圖7。壩體浸潤面較高，絕大部分壩體均處于飽和狀態(tài)，總體呈現(xiàn)向下游降低的趨勢。浸潤面在墊層上下游側(cè)形成突降，但對距兩岸較遠(yuǎn)位置處的地下水滲流場影響不大。因此，墊層滲透系數(shù)對在墊層上、下游側(cè)浸潤面位勢分布較敏感，而對較遠(yuǎn)處的地下水滲流場不敏感。

表6中給出了墊層下游浸潤面最高位置。由表6可知，可知隨著墊層滲透系數(shù)的增大，墊層后浸潤面最高點(diǎn)高程從984. 33 m升高至987. 70 m，且削減水頭百分率從7. 43%降低至3. 62%。因此，隨著墊層滲透系數(shù)的增加，上水庫庫區(qū)和兩岸壩肩以及周邊較小范圍內(nèi)的地下水位升高，墊層削減水頭百分率下降。由于面板未施工，墊層直接面臨水頭作用，墊層滲透系數(shù)越大，削減水頭率越低，而墊層后的過渡層和堆石區(qū)滲透系數(shù)較大，庫水會(huì)快速通過壩體流向下游。

墊層和堆石體的最大平均滲透坡降見表7。與D-3-3相比，墊層滲透系數(shù)從1. 0×10-5 m/s增至5. 0×10-5 m/s，墊層平均最大滲透坡降減小了1. 94，墊層滲透系數(shù)對墊層平均最大滲透坡降較為敏感。并且隨著墊層滲透系數(shù)的增大，墊層的平均最大滲透坡降降低。此外，堆石區(qū)最大平均滲透坡降較小，小于其允許滲透坡降，滿足滲透穩(wěn)定要求。

墊層作為大壩第二道防滲體系，其滲透系數(shù)選取至關(guān)重要，建議墊層滲透系數(shù)取大于5. 0×10-5 m/s，可一定程度上提高臨時(shí)度汛條件下墊層滲透穩(wěn)定安全。

3 結(jié)論

a））在止水失效和面板開裂前后，面板、基巖滲透坡降基本分別穩(wěn)定在177. 50和0. 17左右，當(dāng)止水失效、面板開裂時(shí)，滲流量分別增加12. 23%和40. 89%。因此，由于裂縫的存在，通過混凝土面板滲流量明顯增大。隨著面板止水結(jié)構(gòu)失效、面板開裂位置降低，通過面板的滲流量增大。

b））臨時(shí)度汛工況下，隨著墊層滲透系數(shù)的增大，墊層后浸潤面最高點(diǎn)高程抬升了3. 37 m，且削減水頭百分率從7. 43%降低至3. 62%。墊層滲透系數(shù)越大，削減水頭率越低。

c））臨時(shí)度汛工況下，墊層滲透系數(shù)對在墊層上、下游側(cè)浸潤面位勢分布、墊層平均最大滲透坡降較敏感，而對較遠(yuǎn)處的地下水滲流場不敏感。隨著墊層滲透系數(shù)增大，墊層平均滲透坡降減小。墊層滲透系數(shù)從1. 0×10-5 m/s增大至5. 0×10-5 m/s時(shí)，墊層平均最大滲透坡降減小了1. 94。建議墊層滲透系數(shù)宜大于5. 0×10-5 m/s，可提高臨時(shí)度汛工況下墊層滲透穩(wěn)定安全。

參考文獻(xiàn)：

［1］李銳青. 混凝土面板堆石壩開裂統(tǒng)計(jì)及堆石壩滲流分析研究［D］. 西安：西安理工大學(xué)，2022.

［2］張國新，厲易生. 堆石壩面板收縮性貫穿裂縫的理論分析及防裂措施［J］. 水力發(fā)電學(xué)報(bào)，2005，24（3）：30-33.

［3］ARICI Y. Investigation of the cracking of CFRD face plates［J］. Computers and Geotechnics， 2011， 38（7）： 905-916.

［4］WANG Z J， LIU S H， VALLEJO L， et al. Numerical analysis of the causes of face slab cracks in Gongboxia rockfill dam［J］. Engineering Geology， 2014， 181： 224-232.

［5］JIA J S. A technical review of hydro-project development in China［J］. Engineering， 2016， 2（3）： 302-312.

［6］姜媛媛，唐瑜，辛俊生. 多諾水電站面板堆石壩滲流特性研究［J］. 水電站設(shè)計(jì)， 2014，30（3）：32-35.

［7］林起明，龔軍. 面板開裂對面板堆石壩壩體滲流特性的影響分析［J］. 吉林水利，2012（10）：1-4，7.

［8］張曉強(qiáng)，盧廷浩，周愛兆. 面板壩垂直縫失效滲流場有限元模擬［J］. 水利科技與經(jīng)濟(jì)，2006，12（12）：801-802，809.

［9］溫鵬業(yè). 超高面板堆石壩止水失效分析［D］. 西安：西安理工大學(xué)，2023.

［10］熊璐，黨發(fā)寧，張海鋒，等. 某混凝土面板堆石壩止水失效后的滲流場分析［J］. 水電能源科學(xué)，2015，33（6）：80-83，75.

［11］馬國杰，朱晟，王匯明. 某抽水蓄能電站上庫堆石壩防滲面板允許裂縫寬度分析［J］. 水電能源科學(xué)，2023，41（2）：86-89.

［12］黨嬌特 .特高面板堆石壩施工期度汛滲透穩(wěn)定研究［D］. 西安：西安理工大學(xué)，2023.

［13］張鳳財(cái)，姜媛媛，竇向賢. 猴子巖水電站面板堆石壩三維滲流計(jì)算分析［J］. 水力發(fā)電，2018，44（11）：80-83.

［14］喬蓓，江德軍，柳玉蘭. 面板堆石壩垂直縫失效后集中滲流特性研究［J］. 水電能源科學(xué)，2023，41（4）：123-126.

［15］GAN L， CHEN G Y， SHEN Z Z. A new approach to permeability inversion of fractured rock masses and its engineering application ［J］. Water， 2020， 12（3）： 734-750.

［16］ZHANG W B， SHEN Z Z， CHEN G Y， et al. Optimization design and assessment of the effect of seepage control at reservoir sites under karst conditions： A case study in Anhui Province， China［J］. Hydrogeology Journal， 2021， 29（5）： 1831-1855.

［17］楊金孟，趙蘭浩，沈振中， 等. 某抽水蓄能電站上水庫防滲帷幕深度優(yōu)選研究［J］. 水資源與水工程學(xué)報(bào)，2021，32（2）：209-216.

［18］沈振中，江流，沈長松. 復(fù)合土工膜缺陷滲漏試驗(yàn)的飽和-非飽和滲流有限元模擬［J］. 水利學(xué)報(bào)，2009，40（9）：1091-1095.

［19］朱岳明，龔道勇. 三維飽和-非飽和滲流場求解及其逸出面邊界條件處理［J］. 水科學(xué)進(jìn)展，2003，14（1）：67-71.

［20］姜媛媛. 飽和—非飽和滲流影響下非連續(xù)性巖體邊坡穩(wěn)定分析方法研究［D］. 南京：河海大學(xué)，2005.

［21］岑威鈞，王蒙，楊志祥.（ 復(fù)合）土工膜防滲土石壩飽和-非飽和滲流特性［J］. 水利水電科技進(jìn)展，2012，32（3）： 6-9.

［22］祝林彥，楊金孟，趙邵峰. 面板堆石壩主要分區(qū)滲透系數(shù)反演及敏感性分析［J］. 水電能源科學(xué)，2024，42（2）：99-103.

［23］王相烜，劉得潭，何傳凱，等. 長河壩水電站特高心墻堆石壩雙防滲墻滲流控制特性反演分析［J］. 水利水電科技進(jìn)展，2023，43（5）：88-93.

［24］呂興棟，李家正. 面板堆石壩混凝土面板裂縫現(xiàn)狀、成因與防裂技術(shù)進(jìn)展［J］. 長江科學(xué)院院報(bào)，2021，38（11）：127-134，141.