蔣買(mǎi)勇，劉賢鵬，楊曉雅，李金友

（1.湖南水利水電職業(yè)技術(shù)學(xué)院，湖南長(zhǎng)沙 410131；2.武漢大學(xué)水利水電學(xué)院，湖北武漢 430072；3.長(zhǎng)沙縣水利局，湖南長(zhǎng)沙 410100）

0 引言

滲流問(wèn)題對(duì)于土石壩安全而言是一個(gè)至關(guān)重要的影響因素，近年來(lái)在處理土石壩滲流問(wèn)題上，越來(lái)越多的研究者選擇采用數(shù)值分析的方法，數(shù)值分析中的有限元法得以廣泛應(yīng)用。傳統(tǒng)的二維滲流計(jì)算方法無(wú)法從宏觀實(shí)際把握大壩、壩基及兩岸繞壩的滲流規(guī)律，預(yù)測(cè)水庫(kù)大壩除險(xiǎn)加固后的真實(shí)滲流狀況，三維滲流計(jì)算克服了二維計(jì)算的弊端，能適用于復(fù)雜邊界、多種介質(zhì)的滲流分析，且具有較高的精度［1，2］。

運(yùn)用有限元法解決土石壩的滲流問(wèn)題，實(shí)際是利用有限元軟件的計(jì)算功能來(lái)求解滲流微分方程，模擬出土石壩的滲流場(chǎng)以及得到滲流自由面等，從而計(jì)算出壩體的滲透坡降和總的滲漏量［3-5］。

隨著監(jiān)測(cè)技術(shù)的發(fā)展和應(yīng)用，近年來(lái)很多新建和經(jīng)過(guò)除險(xiǎn)加固的大壩都埋設(shè)有監(jiān)測(cè)儀器，在工程運(yùn)行管理中，對(duì)監(jiān)測(cè)資料的整編和分析也越來(lái)越重視。許多研究者以大壩實(shí)測(cè)資料為基礎(chǔ)，結(jié)合土石壩滲流理論和有限元計(jì)算結(jié)果，進(jìn)行大壩滲流性態(tài)分析，也常被應(yīng)用于驗(yàn)證有限元計(jì)算模型的合理性［6-8］。

本文通過(guò)三維滲流專(zhuān)業(yè)軟件以某水庫(kù)主壩為原型建模，采用有限元法對(duì)大壩的滲流場(chǎng)進(jìn)行三維模擬，結(jié)合大壩滲流監(jiān)測(cè)資料對(duì)有限元計(jì)算結(jié)果進(jìn)行全面的反饋分析。

1 工程概況

1.1 工程歷次病險(xiǎn)癥狀及處置情況

某水庫(kù)大壩為均質(zhì)土壩，壩頂高程102.40 m，寬7.5 m，壩長(zhǎng)159.0 m，最大壩高30.6 m。大壩上游在76.00 m 高程以下為拋石護(hù)腳，76.00 m 高程以上至84.00 m，設(shè)干砌石護(hù)坡，高程84.00 m 至102.40 m 之間坡面設(shè)砼預(yù)制塊護(hù)坡；下游壩腳設(shè)有排水棱體，大壩剖面圖見(jiàn)圖1。1980 年6 月庫(kù)水位為85.5 m時(shí)，主壩右岸條形山下游山坡84.0 m 高程處出現(xiàn)滲漏點(diǎn)，后設(shè)置三角形量水堰及進(jìn)行地勘探查。1984 年2 月完成地勘工作，并觀測(cè)庫(kù)水位96 m時(shí)，滲漏量為7～12 L/s。1984 年10 月完成防滲帷幕灌漿。經(jīng)處理后，條形山滲漏有所減小，但又有逐年加大趨勢(shì)。1985 年觀測(cè)最大滲漏量為2.5 L/s，1989 年5.6 L/s，1995 年10.0 L/s，1998年16.0 L/s。大壩歷經(jīng)多次維修加固，受資金和技術(shù)限制，除險(xiǎn)加固不徹底，險(xiǎn)情依然存在，大壩一直處于帶病運(yùn)行狀態(tài)。2000 年安全鑒定，確定為Ⅲ類(lèi)壩；2006 年10 月至2008 年4月，完成了水庫(kù)大壩除險(xiǎn)加固工程，主壩及條形山進(jìn)行了防滲加固處理，主壩防滲采用高噴灌漿和帷幕灌漿進(jìn)行加固處理，條形山防滲采用帷幕灌漿進(jìn)行處理。灌漿設(shè)單排孔，孔距2 m，灌漿終止線深5 Lu 線以下5 m，與大壩防滲帷幕連成整體。整個(gè)帷幕基本完全封閉，除險(xiǎn)加固后實(shí)測(cè)滲漏量為0.32 L/s。

圖1 某水庫(kù)大壩剖面圖（單位：m）Fig.1 Section of a reservoir dam

水庫(kù)兩岸山體雄厚，地表高程較高，不存在低于水庫(kù)正常高水位的鄰谷，水庫(kù)區(qū)庫(kù)岸巖性主要為第四系殘坡積堆積（Qedl）棕紅色～紅褐色含碎石粉質(zhì)黏土，碎石含量大于30%，粒徑一般3～5 cm。黏性土呈可塑～硬塑狀，透水性差，庫(kù)岸總體穩(wěn)定，局部存在的一些松散堆積物的輕微蠕變，對(duì)水庫(kù)影響甚微，土體呈可塑～硬塑狀，具一定抗沖蝕能力。庫(kù)區(qū)沿岸雖溪流發(fā)育，但水庫(kù)周邊及各溪流內(nèi)植被覆蓋較好，水土流失現(xiàn)象不嚴(yán)重，雨季造成水土流失現(xiàn)象較輕微，水庫(kù)淤積輕微。

1.2 大壩滲流監(jiān)測(cè)概況

根據(jù)壩基地質(zhì)條件、除險(xiǎn)加固情況，大壩滲流監(jiān)測(cè)設(shè)計(jì)3個(gè)橫斷面，即斷面0+035、0+072、0+115。0+072 為原河床最大斷面處，0+035、0+115均為除險(xiǎn)加固前滲漏較為嚴(yán)重的部位。

0+035 斷面共設(shè)3 孔3 個(gè)測(cè)點(diǎn)，軸線后設(shè)1 孔1 個(gè)測(cè)點(diǎn)；第2孔位于92.3 m 高程壩坡上，測(cè)點(diǎn)位置設(shè)于壩體與壩基接觸面處；第3孔位于81.8 m高程壩坡，排水棱體前設(shè)1測(cè)點(diǎn)。

0+072斷面共設(shè)3孔6個(gè)測(cè)點(diǎn)，灌漿軸線后設(shè)1孔2個(gè)測(cè)點(diǎn)，壩體內(nèi)設(shè)1 測(cè)點(diǎn)，壩基帷幕后設(shè)1 測(cè)點(diǎn)；第2 孔位于92.3 m 高程壩坡上，壩體與壩基各設(shè)1 個(gè)測(cè)點(diǎn)；第3 孔位于81.8 m 高程壩坡上，位于排水棱體前，壩體與壩基各設(shè)1個(gè)測(cè)點(diǎn)。

0+115 斷面共設(shè)3 孔3 個(gè)測(cè)點(diǎn)，軸線后設(shè)1 孔1 個(gè)測(cè)點(diǎn)；第2孔位于92.3 m 高程壩坡上，測(cè)點(diǎn)位置設(shè)于壩體與壩基接觸面處。第3孔位于81.8 m高程壩坡，排水棱體前設(shè)1測(cè)點(diǎn)。

各滲流監(jiān)測(cè)點(diǎn)布置信息詳見(jiàn)表1。

表1 滲流監(jiān)測(cè)點(diǎn)布置信息表Tab.1 Seepage monitoring point layout information table

2 有限元計(jì)算

2.1 計(jì)算模型

基于上述監(jiān)測(cè)資料分析，為評(píng)估除險(xiǎn)加固項(xiàng)目對(duì)壩體進(jìn)行高噴灌漿和壩下灌漿帷幕的防滲效果，重點(diǎn)研究除險(xiǎn)加固前滲漏較為嚴(yán)重的部位，建立了壩體壩基三維滲流模型，供對(duì)壩體的除險(xiǎn)加固效果評(píng)價(jià)。三維整體有限元模型以正東方向?yàn)閄軸正方向，正北方向?yàn)閅軸正方向，Z軸豎直向上。模型對(duì)均質(zhì)壩壩體、壩基及壩體內(nèi)高噴灌漿、壩下帷幕、左右岸防滲帷幕、及各主要地層滲透單元和河床覆蓋層等進(jìn)行了較為全面的模擬，計(jì)算模型范圍如下：

（1）沿高程方向（Z向）：底部計(jì)算邊界為沿建基面向下延伸約3倍壩高；

（2）東西向（X向）：計(jì)算范圍自壩軸線中心點(diǎn)往東西向各取375 m（均大于10倍壩高），即模型東西向總長(zhǎng)750 m；

（3）南北向：計(jì)算范圍自壩軸中心點(diǎn)往南北方向各取375 m（均大于10倍壩高），即模型南北向總長(zhǎng)750 m。

整個(gè)有限元模型均采用8 結(jié)點(diǎn)等參單元進(jìn)行網(wǎng)格剖分，可以較好地保證滲流分析的精度及成果的可靠性，模型共剖分399 739 個(gè)單元，75 302 個(gè)結(jié)點(diǎn)，模型采用四面體單元［4］。模型及有限元網(wǎng)格劃分如圖2、3所示。

圖2 三維滲流整體有限元網(wǎng)格圖Fig.2 3D integral finite element mesh diagram of seepage flow

圖3 典型壩段壩體結(jié)構(gòu)有限元網(wǎng)格及材料分區(qū)示意圖Fig.3 Diagram of finite element grid and material partition of typical dam section

2.2 計(jì)算參數(shù)

（1）計(jì)算參數(shù)選取。結(jié)合工程地質(zhì)報(bào)告和相關(guān)工程經(jīng)驗(yàn)等資料，本次滲流計(jì)算采用的綜合參數(shù)如表2。

表2 某水庫(kù)大壩三維滲流有限元計(jì)算參數(shù)選用情況Tab.2 Selection of 3D seepage finite element calculation parameters for a reservoir dam

（2）局部滲透穩(wěn)定分析。根據(jù)《水庫(kù)大壩安全評(píng)價(jià)導(dǎo)則》SL258-2017，大壩滲流安全評(píng)價(jià)可采用現(xiàn)場(chǎng)檢查法、監(jiān)測(cè)資料分析法、計(jì)算分析法和經(jīng)驗(yàn)類(lèi)比法，宜綜合使用。本文主要使用計(jì)算分析法并結(jié)合現(xiàn)場(chǎng)檢查法與監(jiān)測(cè)資料分析法進(jìn)行綜合分析。通過(guò)理論方法或數(shù)值模型計(jì)算大壩滲流量、水頭、滲流壓力、滲透坡降等水力要素及其分布，繪制流網(wǎng)圖，評(píng)判防滲體的防滲效果，以及關(guān)鍵部位滲透坡降是否小于允許滲透坡降，滲流出逸點(diǎn)高程是否在貼坡反濾保護(hù)范圍內(nèi)作為評(píng)判依據(jù)。

心墻土料為黏性土，其滲透變形形式主要為流土與接觸流失。本工程心墻位于壩體內(nèi)部，內(nèi)部滲流方向自上而下，下游設(shè)反濾層保護(hù)，可能發(fā)生滲透變形形式為接觸流失，心墻土料允許滲透坡降參照規(guī)范及經(jīng)驗(yàn)取4～6。

2.3 計(jì)算邊界條件

上游水庫(kù)淹沒(méi)區(qū)內(nèi)的表面結(jié)點(diǎn)及下游水庫(kù)淹沒(méi)區(qū)內(nèi)的表面結(jié)點(diǎn)均取為定水頭邊界，具體水位按不同蓄水位確定。模型上、下游邊界取隔水邊界，模型底部邊界也取隔水邊界。

2.4 計(jì)算工況

工況一：水庫(kù)上游為正常蓄水位98.70 m，相應(yīng)下游水位為71.80 m。

工況二：水庫(kù)上游為設(shè)計(jì)洪水位100.34 m，相應(yīng)下游水位71.80 m。

2.5 計(jì)算結(jié)果及分析

結(jié)合大壩特征和現(xiàn)有滲流觀測(cè)設(shè)施布置情況以及大壩安全監(jiān)測(cè)所反映出的滲透破壞薄弱部位信息，擬對(duì)0+035斷面、0+072 斷面、0+115 斷面這三處典型剖面進(jìn)行滲流計(jì)算分析，計(jì)算結(jié)果如圖4～12所示。

圖4 0+035斷面滲透坡降圖Fig.4 Permeability gradient diagram of 0+035 section

（1）0+035斷面。根據(jù)剖面滲透坡降圖可知，坡降較大的部位主要集中在壩體高噴灌漿所形成的防滲墻以及壩下帷幕，正常蓄水位情況下坡降值集中在0.33～0.75 之間，在設(shè)計(jì)洪水位情況下坡降值集中在0.44～0.80 之間，其他部位未見(jiàn)較大坡降。

圖5 0+035斷面等水頭線圖Fig.5 Equal head diagram of 0+035 section

圖6 0+035斷面壓力水頭及自由面分布圖Fig.6 Section pressure head and free surface distribution diagram of 0+035 section

根據(jù)剖面等水頭線圖可知，等水頭線在高噴灌漿所形成的防滲墻以及壩下帷幕處比較密集，下游相對(duì)密集，上游較為稀疏。

根據(jù)剖面壓力水頭及自由面分布圖可知，工況二的滲流自由面相較工況一有所抬升，且滲流自由面在高噴灌漿所形成的防滲墻以及壩下帷幕處急劇降低，在壩體黏土填筑區(qū)域從上游至下游均勻下降，并在排水棱體區(qū)域溢出。以上滲流狀況反映了該斷面在采取除險(xiǎn)加固后防滲墻以及壩下帷幕防滲性能得到了較大改善。

（2）0+072斷面。根據(jù)剖面滲透坡降圖可知，坡降較大的部位主要集中在近建基面附近的壩下帷幕，工況一的坡降值集中在0.4～1.0 之間，工況二的坡降值集中在0.53～1.2 之間，其他部位未見(jiàn)較大坡降。主要原因考慮是工程初建時(shí)期對(duì)于基礎(chǔ)的清理深度不夠，加之建設(shè)期采用土墻和散碎巖體回填壩基且采用農(nóng)用器具進(jìn)行夯實(shí)，材料和工藝都存在技術(shù)層面的缺陷，致使壩基回填區(qū)和弱風(fēng)化巖層之間存在滲漏薄弱層。灌漿帷幕厚度小于壩內(nèi)高噴防滲墻厚度，從而在灌漿帷幕穿過(guò)上述滲漏薄弱層區(qū)域出現(xiàn)較大滲透坡降。此處出現(xiàn)的較大滲透坡降僅分布在帷幕及其極近的下游側(cè)，反應(yīng)帷幕防滲體系運(yùn)行狀態(tài)良好。

根據(jù)剖面等水頭線圖可知，等水頭線在高噴灌漿所形成的防滲墻以及壩下帷幕處比較密集，其他部位較為稀疏。

圖7 0+072斷面滲透坡降圖Fig.7 Permeability gradient diagram of 0+072 section

圖8 0+072斷面等水頭線圖Fig.8 Equal head diagram of 0+072 section

圖9 0+072斷面壓力水頭及自由面分布圖Fig.9 Section pressure head and free surface distribution diagram of 0+072 section

根據(jù)剖面壓力水頭及自由面分布圖可知，工況二的滲流自由面相較工況一有所抬升，且不同工況下的滲流自由面均在高噴灌漿所形成的防滲墻以及壩下帷幕處急劇降低，在壩體黏土填筑區(qū)域從上游至下游均勻下降，并在排水棱體區(qū)域溢出，工況二中排水棱體后方設(shè)置的塘內(nèi)滲漏量較工況一有所增加。此斷面反映出防滲體系運(yùn)行狀態(tài)良好，滲透規(guī)律合理。

（3）0+115斷面。根據(jù)剖面滲透坡降圖可知，工況一的滲透坡降值集中在0.2～1.2之間，工況二的滲透坡降值集中在0.2～1.3 之間，其他部位未見(jiàn)較大滲透坡降，滲透坡降最值出現(xiàn)在灌漿帷幕上部。

根據(jù)剖面等水頭線圖可知，等水頭線在高噴灌漿所形成的防滲墻以及壩下帷幕處比較密集，其他部位較為稀疏，分布合理。

圖10 0+115斷面滲透坡降圖Fig.10 Permeability gradient diagram of 0+115 section

圖11 0+115斷面等水頭線圖Fig.11 Equal head diagram of 0+115 section

圖12 0+115斷面壓力水頭及自由面分布圖Fig.12 Section pressure head and free surface distribution diagram of 0+115 section

根據(jù)剖面壓力水頭及自由面分布圖可知，工況二的滲流自由面相較工況一有所抬升，壩后塘內(nèi)滲漏量有所增加。滲流自由面在高噴灌漿所形成的防滲墻以及壩下帷幕處下降明顯，在壩體黏土填筑區(qū)域下降均勻，自由面從防滲墻至下游呈現(xiàn)開(kāi)口向下的拋物線形狀，傾斜方向與下游坡面傾向逐漸貼合，兩面之間距離也逐漸靠近，溢出點(diǎn)位于排水棱體上部。該斷面下游側(cè)尤其是排水棱體上部區(qū)域存在潛在滲流溢出點(diǎn)，滲流溢出點(diǎn)位置偏上的原因一方面是除險(xiǎn)加固的規(guī)模和質(zhì)量，另一方面是加固材料老化和劣化。這與實(shí)際觀測(cè)到右岸下游該斷面附近出現(xiàn)的滲漏現(xiàn)象以及安全監(jiān)測(cè)機(jī)房漏水等事件相佐證。

3 基于監(jiān)測(cè)資料的滲流對(duì)比分析

0+035 斷面、0+072 斷面、0+115 斷面滲壓計(jì)監(jiān)測(cè)值過(guò)程線分別見(jiàn)圖13～15。根據(jù)監(jiān)測(cè)值，上游水位在90 m 高程附近波動(dòng)，滲壓計(jì)測(cè)值平穩(wěn)，滲壓計(jì)反映出各處滲透壓力與上游水位大致呈正相關(guān)。

圖13 0+035斷面滲壓計(jì)監(jiān)測(cè)值過(guò)程線圖Fig.13 Process diagram of the monitoring value of the osmometer of 0+035 section

參照工況一等水頭線圖，0+072 斷面埋設(shè)高程在69.5 m 的SZ2-1-2、SZ2-2-2 和SZ2-3-2 三支滲壓計(jì)在建基面以下，經(jīng)過(guò)帷幕阻滲作用，滲透水頭下降顯著，SZ2-1-2 布設(shè)軸距為3 m，處于帷幕后坡降較大區(qū)域，滲透壓力大，測(cè)值在77 m水頭，沿著滲徑水頭線逐漸降低，尤其是SZ2-3-2（排水棱體上游側(cè)附近）水頭值為70 m，表明防滲系統(tǒng)運(yùn)行狀況良好。埋設(shè)高程在73.3 m 的SZ2-2-1（軸距27.1 m）測(cè)值在78 m 高程附近，同高程軸距為56.2 m 布設(shè)的SZ2-3-1 測(cè)值降至74 m。各埋設(shè)部位滲壓計(jì)測(cè)值符合工況一等水頭線圖相應(yīng)點(diǎn)位分布規(guī)律，表明有限元計(jì)算結(jié)果能反映出較為真實(shí)的土石壩滲流場(chǎng)情況。

4 結(jié)語(yǔ)

通過(guò)對(duì)某水庫(kù)主壩在正常蓄水位和設(shè)計(jì)洪水位兩種工況下的滲流性態(tài)進(jìn)行有限元計(jì)算，以及與滲流監(jiān)測(cè)資料進(jìn)行對(duì)比分析可知。

（1）采用三維有限元法對(duì)大壩進(jìn)行滲流計(jì)算，其結(jié)果符合實(shí)測(cè)資料的分布規(guī)律，能為大壩滲流問(wèn)題分析提供可靠的依據(jù)。

（2）除險(xiǎn)加固后的高噴灌漿防滲墻和帷幕灌漿防滲體系能正常運(yùn)行，且土石壩整體的滲流規(guī)律合理。

（3）在水庫(kù)運(yùn)行中，應(yīng)注意大壩右岸下游測(cè)尤其是棱體排水上部區(qū)域的潛在滲流溢出點(diǎn)，加強(qiáng)對(duì)該區(qū)域的安全監(jiān)測(cè)，避免影響大壩安全。

圖14 0+072斷面滲壓計(jì)監(jiān)測(cè)值過(guò)程線圖Fig.14 Process diagram of the monitoring value of the osmometer of 0+072 section

圖15 0+115斷面滲壓計(jì)監(jiān)測(cè)值過(guò)程線圖Fig.15 Process diagram of the monitoring value of the osmometer of 0+115 section