深厚覆蓋層上高土石壩防滲墻裂縫滲流分析模型及應用研究

嚴 俊，溫彥鋒，璩愛玉，唐瑜蓮

(1.中國水利水電科學研究院流域水循環模擬與調控國家重點實驗室，北京100048；2.環境保護部環境規劃院，北京100048；3.浙江省水利水電咨詢中心，浙江杭州310020)

深厚覆蓋層上高土石壩防滲墻裂縫滲流分析模型及應用研究

嚴 俊1，溫彥鋒1，璩愛玉2，唐瑜蓮3

(1.中國水利水電科學研究院流域水循環模擬與調控國家重點實驗室，北京100048；2.環境保護部環境規劃院，北京100048；3.浙江省水利水電咨詢中心，浙江杭州310020)

深厚覆蓋層上修建混凝土防滲墻在其施工和運行期由于多種原因而容易發生裂縫，可能影響覆蓋層地基的滲透穩定性，對工程安全造成難以估量的影響。為有效地求解深厚覆蓋層上高土石壩防滲墻裂縫滲流問題，在改進節點虛流量法的基礎上引入無厚度的縫單元來模擬防滲墻裂縫滲流行為。以瀑布溝深厚覆蓋層工程滲控為例，計算了混凝土防滲墻裂縫未充填和充填下的三維滲流場，得到不同裂縫條件下滲流場的水頭分布、浸潤線以及滲漏量，系統分析了深厚覆蓋層地基中混凝土防滲墻垂直裂縫條件下的穩定滲流場規律和特點。結果表明，該方法能對滲流逸出點和浸潤線能準確定位，并能很好地模擬混凝土防滲墻裂縫對壩區滲流場的影響，可以為該類問題的分析提供參考。

深厚覆蓋層；防滲墻；裂縫；滲流分析

1 研究背景

深厚覆蓋層是近年來我國在開發西南山區水電工程常常遇到的一種復雜的地質條件，其厚度多達幾十米，有的甚至達到幾百米[1]。目前國內對于該類地基的滲流控制多采用垂直防滲方式，而混凝土防滲墻則是其中首選的、最可靠的工程措施[2]。但是，混凝土防滲墻在施工和運行期常常會發生開裂：如在防滲墻施工過程中，由于地基條件復雜[3- 4]、粘土泥漿質量差[5]、施工工藝[6]等問題常會導致墻體產生劈裂裂縫、下部開叉現象；在大壩建成后，防滲墻可能會因不均勻沉降和變形產生拉裂縫[7- 9]；水庫蓄水后，墻體受類似于彎壓組合荷載作用、在壩體及基巖約束下將產生撓曲，可能會出現縱向裂縫；地震條件下，混凝土防滲墻也可能因動荷載的作用產生規模不等的裂縫[10-12]。混凝土防滲墻的這些裂縫，尤其是貫穿性的橫縫、縱縫，隨著裂縫開度的增加，將逐漸削弱防滲墻局部的截滲性能、加劇裂縫處的滲流，對墻體的下游造成沖刷，影響覆蓋層地基的滲透穩定性[13-14]，對工程安全造成難以估量的影響[15]。此外，混凝土防滲墻開裂后，在上覆土體壓力以及地下水滲透力的作用下，細顆粒可能會進入裂縫中，對防滲墻裂縫起到一定的愈合作用，改善裂縫局部的地下水流態[16]，因此，迫切需要研究深厚覆蓋層地基內混凝土防滲墻開裂后的分析模型，以掌握其對地基滲流的影響規律。

目前，裂縫滲流計算模型多采用空間薄層單元的等效連續體模型方法進行模擬。該方法可以利用等效連續體模型較成熟的理論基礎和豐富的運用經驗，但其存在的缺點是，裂縫單元的尺寸很小，與周邊單元尺寸存在較大的差異，容易影響求解的精度；而對裂縫周邊單元進行加密以適應裂縫尺寸又容易增加求解的規模。為解決這些問題，并有效地求解深厚覆蓋層上高土石壩防滲墻裂縫滲流問題，本文在改進節點虛流量法的基礎上引入無厚度的縫單元來模擬防滲墻裂縫，從而達到既能求解高土石壩整體滲流問題，又能突出裂縫對整個滲流場的影響，并結合工程示例來驗證該方法的可靠性。

2 無壓滲流場求解的改進節點虛流量法

文獻[17]中提出了結點虛流量法，為求不變網格的有自由面滲流問題開辟了新思路：該方法以自由面為界將整個計算域Ω分為實域Ω1(自由面以下)和虛域Ω2(自由面以上)。由于自由面和逸出點的位置都是未知的，需要通過多步迭代來求解，其相應的有限元求解的支配方程如下：

[K]{h}={Q}-{Q2}+{ΔQ}(令{ΔQ}=[K2]{h})

(1)

式中，[K]、[K2]分別為引入邊界條件后的計算域全域及虛域分別貢獻的整體滲透矩陣；{h}為未知節點的水頭列陣；{Q}、{Q2}分別為已知水頭節點、內部源匯項和流量邊界對計算域的全域、虛域域貢獻的流量列陣。

3 無厚度的防滲墻裂縫滲流分析模型

文獻[18]中用無厚度裂縫單元模型來模擬混凝土面板裂縫，取得了很好的效果，本文引入該無厚度裂縫單元來模擬混凝土防滲墻的裂縫：縫面水力隙寬遠小于裂縫另外兩向的尺寸，法向透水能力遠大于混凝土本體的透水能力，縫面方向水頭損失非常小，因此地下水在裂縫中的運動可以簡化為準二維的滲流。此時裂縫單元的傳導矩陣計算公式為

(r，s=1，2，…，n)

(2)

式中，sf為裂縫單元域；r和s均為單元節點；Nr和Ns為裂縫縫面單元插值函數；n為縫面單元的節點數。

在得到每個裂縫單元的傳導矩陣后，與其他部位的單元一樣，裂縫單元的滲流特性主要通過水頭連續條件和流量平衡條件組裝進整體滲透矩陣[K]，其中，水頭連續性條件為裂縫單元任一點的水頭與裂縫壁上的水頭相同；流量平衡條件是滿足改進節點虛流量法的流量平衡條件。

4 工程應用分析

4.1 工程背景

瀑布溝水電站壩基覆蓋層最大厚度75.36 m，采用2道混凝土防滲墻全封閉防滲，墻厚1.2 m，中心間距14 m。上游墻670.00 m高程以下最大深度76.85 m，下游墻670.00 m高程以下最大深度75.55 m。上游側防滲墻670.00 m高程以上頂部直接插入大壩心墻，插入深度10 m；下游側防滲墻頂部設3.5 m×4 m的灌漿兼觀測廊道。

為減少不均勻沉降，防止壩體開裂，在心墻與兩岸基巖接觸面上鋪設3 m厚高塑性粘土，在防滲墻頂、廊道周圍和心墻底部也鋪設高塑性粘土。

4.2 有限元分析模型及條件

4.2.1 有限元分析模型

圖1 瀑布溝壩區滲流場三維網格計算模型

圖2 計算模型中壩體典型剖面網格情況

為弄清混凝土防滲墻裂縫對地基乃至壩區的整體滲流場的影響程度，在分析中對大壩中的主要水工結構物也進行了較精細的模擬，其中包括壩體結構(心墻、反濾層等)、壩殼、防滲墻、防滲帷幕、廠房(包括廠房上游面的排水孔等)、各灌漿及排水廊道等，對地基深厚覆蓋層也進行了較細致的模擬。圖1為三維有限元分析模型，圖2為網格中壩體典型斷面圖。模型生成后主要由六面體8結點等參元和局部區五面體6結點過渡性等參元組成，以前者為主，共有67 754個單元和74 592個結點。

表1 滲流計算參數

4.2.2 有限元分析邊界及參數

分析域四周截取邊界條件分別假定為：上游截取邊界、下游截取邊界以及底邊界均視為隔水邊界面；左岸截取邊界和右岸截取邊界則均取為第一類邊界條件，即已知水頭邊界條件；對于地表邊界，壩軸線上游側，低于河或庫水位的地方為已知水頭邊界(上游正常蓄水位為850.00 m)，高于河或庫水位的地方為滲流可能逸出面；在壩軸線下游側，同樣低于下游水位的地方為已知水頭邊界條件(正常尾水位670.00 m)，高于下游水位的地方均為滲流可能逸出面。計算中所采取的各種材料的滲透系數取值見表1。

4.2.3 計算分析方案

為弄清混凝土防滲墻裂縫對地基乃至壩區的整體滲流場的影響程度，考慮如下混凝土防滲墻裂縫滲流特性計算方案：

(1)方案1，水庫正常運行方案。上游庫水位850.00 m，下游尾水位670.00 m；覆蓋層中防滲墻、巖體中防滲帷幕、壩內心墻、反濾層、過渡層以及排水孔幕等按設計布置情況進行模擬，防滲墻按設計要求布置成兩排。

(2)方案2～4。考慮防滲墻垂直開裂度的影響，裂縫為貫穿縫的，縫寬分別為1、2、5 mm。缺陷防滲墻段不考慮充填，滲透系數取為砂礫石覆蓋層的10倍，即0.8 cm/s，其余各種計算條件同方案1。

(3)方案5～7。考慮防滲墻垂直開裂度的影響，裂縫為貫穿縫的，縫寬分別為1、2、5 mm。缺陷防滲墻段考慮砂礫石充填，滲透系數與砂礫石相似，為8.0×10-2cm/s，其余各種計算條件同方案1。

4.2.4 裂縫模擬情況

主要考慮墻體垂直向開裂的情況出現，計算分析這些垂直向開裂的壩基。整個壩基滲流場宏觀特性的影響：開裂狀態為上、下游防滲墻同時在對應高程位置處開裂，裂縫的隙寬分別為1、2、5 mm，且分別考慮裂縫未充填和充填，其中未充填時，裂縫部位的滲透系數取為砂礫石覆蓋層的10倍，即0.8 cm/s；裂縫充填情況主要考慮被覆蓋層中的砂礫石料充填，裂縫部位的滲透系數與砂礫石相似，滲透系數取為0.08 cm/s。計算分析中考慮最不利的情況，即考慮裂縫為貫穿性裂縫。

4.3 計算結果分析

4.3.1 水庫正常運行時的滲流場分布

圖3為正常運行期方案(方案1)的壩址滲流場滲流自由面等高程線的分布，圖4為典型剖面及局部滲流水頭的等值線分布。從圖3、4中可以看出，整個滲流場的水頭分布規律合理，水頭等值線形態、走向和密集程度都較準確地反映了相應區域防滲或排水滲控措施的特點、滲流特性和邊界條件，計算域內的主要防滲和排水措施都得到了細致的模擬，滲控效果也及時得到了正確反映：①左岸。在巖體及防滲帷幕的共同作用下，帷幕后的地下水水位已經跌至720.00 m左右，廠房前端的排水廊道及排水孔幕的消減水頭的作用明顯，在防滲帷幕和排水幕的聯合作用下，廠區上游來水得到了有效控制，自由面下降明顯，廠房上游端地下水位已跌至630.00 m左右。②壩中部位。等水頭線走勢有明顯的規律性，主要集中在滲透系數相對較小的防滲心墻和壩基上部的弱透水層內，自由面在防滲心墻內中急劇下降。壩體心墻表現出良好的抗滲能力，壩體中的水頭等值線基本都集中分布在相對弱透水能力的防滲心墻內，心墻防滲作用明顯，有效地降低和控制了壩體自由面的位置，下游逸出點的高程基本保持在671.00 m。壩基中，由于布置了兩道防滲墻，等水頭線明顯集中在兩道防滲墻周圍，防滲墻上下游地下水位下降明顯，防滲墻消減了的水頭損失約170 m，占上下游水位差180 m的94.44%，墻體承擔了相對集中的滲透水頭差，可見防滲墻及其下端布置的防滲帷幕消減水頭的作用明顯，對于深厚覆蓋層地基中地下水滲流有很好的控制作用。③右岸。由于防滲帷幕的滲透性與右岸巖體的滲透性接近，可以看作一均質的防滲體，右岸地下水逐漸下降，水頭等值線呈明顯的規律性分布。

圖3 方案1下的滲流場自由面等水頭線分布(單位：m)

圖4 方案1下典型剖面在正常運行條件下滲流場等水頭線分布(單位：m)

圖5 方案2～4下典型剖面在正常運行條件下滲流場等水頭線分布(單位：m)

圖6 各方案下典型部位滲漏量計算結果

4.3.2 垂直裂縫(未充填)條件下水庫滲流場分布

在方案2～4中考慮防滲墻在壩中典型斷面處出現未被充填的垂直裂縫，裂縫隙寬逐漸增加，分別為1、2 mm和5 mm，圖5分別為3種方案下的典型斷面滲流水頭分布，圖6為各方案下典型位置的滲漏量計算結果。從圖6、7可以看出，庫區其他部位處的等水頭線分布與方案1的等水頭線分布相比，在山體中的等水頭線分布基本一致，壩基中等水頭線仍然主要集中在防滲墻上，壩體中等水頭線分布出現小的變化。但是，隨著防滲墻裂縫隙寬的增加，引起壩體心墻下游側出逸點逐漸抬高，此處的滲透梯度增加，通過防滲墻和防身帷幕的流量稍有增加，壩基中等水頭線分布出現變化：方案2、3下，裂縫隙寬分別為1 mm和2 mm，且裂縫未被充填，此時壩基巖體中的等水頭線分布也與方案1甚為一致，略有變化的部位是在覆蓋層防滲墻附近，兩道防滲墻消減的水頭減小，同時下游防滲墻上等水頭線減少一條，減少的等水頭線同樣出現在墻體下游側的覆蓋層中，致使此處覆蓋層中滲透梯度增加，而且隨著隙寬的增加，出現在下游覆蓋層中等水頭線逐漸向下游移動，下游逸出部位的滲透梯度逐漸增大，只是未出現超過其允許滲透梯度的現象；在方案4下，裂縫隙寬增加至5 mm，此時防滲墻雖能發揮一定的防滲作用，但是防滲墻的防滲效果削弱，下游防滲墻上等水頭線同樣減少一條，墻體上的滲透梯度減小，同時減少的等水頭線出現在下游逸出部位的覆蓋層中，這些部位中出現了超過允許滲透梯度的現象，可見裂隙隙寬達到5 mm時容易對下游逸出部位造成沖刷，對樞紐的長期滲透穩定不利。

圖7 方案5～7下典型剖面在正常運行條件下滲流場等水頭線分布(單位：m)

4.3.3 垂直裂縫(充填)條件下水庫滲流場分布

在方案5～7中，考慮防滲墻上的垂直裂縫被砂礫石充填，圖7分別為3種方案下的典型斷面滲流水頭分布。從圖7可以看出，當防滲墻出現這些被充填的垂直裂縫后，對山體、壩體中的等水頭線分布均沒有什么影響。在充填后的裂縫隙寬為1 mm和2 mm時，隨著裂縫寬度的增加，壩基巖體、防滲墻附近的等水頭線分布和滲透梯度均沒有出現明顯的變化；但是當裂縫隙寬擴大到5 mm時，防滲墻周圍的等水頭線分布情況相對有了改變，等水頭線仍主要集中在防滲墻上，與方案1相比，防滲墻消減水頭減至150 m，占上下游水位差的83.33%，上下游防滲墻仍為壩基覆蓋層中的主要防滲體；同時，下游防滲墻上的等水頭線減少了一條，下游防滲墻滲透梯度減小，下游覆蓋層中出現等水頭線，但覆蓋層滲透梯度變化不大，仍處于安全狀態。

5 結 論

混凝土防滲墻是深厚覆蓋層地基的首選的、最可靠的滲控工程措施，其在施工和運行中容易產生各種裂縫，可能會對庫區周邊的地下水分布造成影響，威脅地基覆蓋層的滲透穩定行。本文采用無厚度的防滲墻裂縫分析模型，結合瀑布溝典型深厚覆蓋層工程，對混凝土防滲墻裂縫充填與未充填條件下的滲流狀態進行了分析：

(1)在混凝土防滲墻出現較小的垂直裂縫(如小于5 mm)時，其影響范圍較小，對整個工程區的滲流狀態影響不大，但是在裂縫處，防滲墻的截滲效果減弱，雖仍能較好地發揮滲控作用，但覆蓋層地基中的水頭分布發生變化，局部滲透比降增大，而且地基中滲漏量也會隨著裂縫的增大而增加；在裂縫被覆蓋層土石料充填時，其裂縫處滲流特性將受充填料的影響，隨著裂縫的增大，覆蓋層局部滲透比降增大，地基中滲漏量也逐漸增大，但較未充填時略有緩和。

(2)無厚度的防滲墻裂縫分析模型能夠很好地模擬地基中混凝土防滲墻裂縫下的滲流特性，而且能夠與工程區其他部位的求解模型較好地結合，計算方法簡便可行，是一種很好的防滲墻裂縫滲流模擬的方法。

[1] 許強, 陳偉, 張倬元, 等. 對我國西南地區河谷深厚覆蓋層成因機理的新認識[J]. 地球科學進展, 2008, 23(5)： 50- 59．

[2] 王清友, 孫萬功, 熊欣, 等. 塑性混凝土防滲墻[M]. 北京： 中國水利水電出版社, 2008: 5- 6．

[3] 于利均. 防滲墻施工中裂縫產生的原因和對策研究[J]. 人民黃河, 2006, 28(12)： 77- 78．

[4] 宋瑞華, 曲苓. 克孜爾水庫右壩肩防滲墻裂縫成因分析及處理措施[J]. 水利水電工程設計, 2009, 28(3)： 3- 5．

[5] 馬希正. 臨海市牛頭山水庫混凝土防滲墻裂縫檢測及處理[J]. 浙江水利科技, 2006(2)： 44- 45．

[6] 蔡長治, 胡春環. 同馬大堤防滲墻施工產生劈裂裂縫的原因[J]. 華北水利水電學院院報, 2004, 25(4)： 24- 26．

[7] 祁偉強, 彭云楓, 袁玉琳, 等. 壩基混凝土防滲墻應力變形三維有限元分析[J]. 水電能源科學, 2012(8)： 63- 66．

[8] 賈華, 何順賓, 伍小玉, 等. 長河壩心墻堆石壩地基防滲墻應力變形分析[J]. 水資源與水工程學報, 2008, 19(3)： 72- 75．

[9] 李存法. 高壓噴射灌漿防滲墻三維有限元應力應變分析[J]. 水利水電技術, 1993(8)： 11- 14．

[10] 高江林, 陳云翔. 基于滲流與應力耦合的防滲墻與壩體相互作用的數值模擬[J]. 水利水運工程學報, 2013(2)： 58- 63．

[11] 李濤, 張曉麗, 李永鑫. 透水壩基薄防滲墻三維有限元模擬分析[J]. 人民黃河, 2005, 27(6)： 53- 55．

[12] 李少明, 秦忠國. 林溪水庫主壩防滲墻應力變形的有限元分析[J]. 人民長江, 2012, 43(19)： 73- 75．

[13] 盛金昌, 趙堅, 速寶玉. 混凝土防滲墻開裂對壩基滲透穩定性的影響[J]. 水利水電科技進展, 2006, 26(1)： 23- 26．

[14] 許國安, 揚凱宏, 武桂生. 小浪底水庫壩基防滲墻裂縫對滲流控制的影響[J]. 人民黃河, 1993(9)： 28- 32．

[15] 李郴娟, 廖華勝, 李連俠, 等. 閘壩防滲墻開裂的數值計算與分析[J]. 四川水力發電, 2008, 27(6)： 89- 91．

[16] 許國安, 武桂生. 混凝土防滲墻出現裂縫的滲流分析研究[J]. 水力發電, 1994(11): 31- 35．

[17] YAN Jun, LIU Sihong. The improved node virtual method for the seepage problem with free surface[C]∥Proceedings of International Conference on Electric Technology and Civil Engineering, 2011: 2760- 2765．

[18] 陳守開, 嚴俊, 李健銘. 面板堆石壩垂直縫破壞下三維滲流場有限元模擬方法[J]. 巖土力學, 2011, 32(11)： 3473- 3478．

ResearchonCrackSeepageAnalysisModelofHighEarth-rockDamonDeepOverburdenandItsApplication

YAN Jun1, WEN Yanfeng1, QU Aiyu2, TANG Yulian3

(1. State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing 100048, China; 2. Chinese Academy for Environmental Planning, Beijing 100048, China; 3. Zhejiang Consulting Center of Water Conservancy and Hydropower Technology, Hangzhou 310020, Zhejiang, China)

The concrete anti-seepage wall built on deep overburden may appear lots of cracks during construction and operating period for various reasons, which would affect the seepage stability of deep overburden and the safety of project. In order to solve the seepage problem of high earth-rock dam on deep overburden with cracks in concrete anti-seepage wall, a crack model with no thickness is introduced based on improved node virtual method. The deep overburden project in Pubugou Hydropower Station is chosen as an example and related three-dimensional seepage field under cracks in concrete wall both un-filled and filled by overburden soil are calculated. The water-head distribution, free surface and seepage under different cracks are obtained. The law and characteristic of steady seepage field of deep overburden with vertical cracks in concrete anti-seepage wall are systematically analyzed. The calculation results show that this crack model can accurately locate seepage escape point and free surface, and can well simulate the influence of cracks in anti-seepage wall to whole seepage field．

deep overburden; concrete anti-seepage wall; crack; seepage analysis

TU43

A

0559- 9342(2017)09- 0039- 06

2017- 03- 09

國家自然科學基金項目(51409278)；國家973計劃課題(2014CB047004)

嚴俊(1984—)，男，湖北襄陽人，博士，主要從事工程滲控、多場耦合效應研究．

(責任編輯焦雪梅)