，，

(陜西水環(huán)境工程勘測設計研究院，陜西 西安 710018)

本文根據(jù)陜西省延川縣袁家溝水庫工程大壩滲流計算，通過AUTO BANK軟件計算大壩的滲流穩(wěn)定，水庫壩基坐落在強透水的強風化巖石層上，這就需要下游壩坡進行強有力的排水措施排出壩面及壩內的滲水，本文通過棱體排水和棱體排水與褥墊排水綜合型式，布設兩種方案，結果表明采用綜合排水型式能夠有效排出壩內水，有效降低浸潤線，構建模型后計算過程簡單，根據(jù)計算結果確定壩坡排水型式，為今后類似工程提供相關設計參考。

1 工程概況

袁家溝水庫位于延川縣城西南約15 km處的袁家溝溝口，距文安驛鎮(zhèn)1.0 km。袁家溝屬清澗河二級支流，文安驛川右岸一級支流，流域面積24.5 km2，河道長12.3 km，河道平均比降13.6‰，袁家溝水庫壩址位于袁家溝河入文安驛川河口以上約0.3 km處，壩址以上控制流域面積24.3 km2。該工程任務以城鎮(zhèn)供水為主。水庫樞紐工程由大壩、溢洪道及輸水放空建筑物等三部分組成?？値烊?82萬 m3，死庫容220萬 m3，起始興利庫容125萬 m3，50年使用年限淤積庫容320萬 m3，終極興利庫容25萬 m3，滯洪庫容137萬 m3。本工程屬Ⅳ等小(1)型工程，相應主要建筑物：大壩、溢洪道、輸水放空洞等為4級，次要建筑物為5級，臨時建筑物為5級。水庫樞紐工程防洪標準按30年一遇設計，300年一遇校核，相應洪峰流量分別為322 m3/s和517 m3/s。

水庫樞紐工程由大壩、溢洪道、輸水放空建筑物三部分組成。其中大壩為碾壓式均質土壩，壩頂高程為918.5 m，正常蓄水位912.90 m，設計洪水位915.30 m，校核洪水位917.10 m，死水位908.30 m，壩頂寬度6 m，最大壩高35.5 m，壩頂總長為189 m。壩體上游邊坡為1：3；下游壩坡坡比均為1：2.75。按壩體結構要求，該均質壩分為壩體、護坡、排水棱體三個區(qū)。其中壩體采用黃土狀壤土進行填筑，排水棱體和褥墊采用石料和混凝土粗細骨料進行加工，上游坡面采用15 cm厚C25鋼筋混凝土護坡，下設20 cm級配碎石墊層、20 cm細砂墊層，上游護坡頂部高程與壩頂齊平為918.50 m，底部交至現(xiàn)狀地面高程，下游坡面采用C20混凝土拱格內植草皮護坡。為方便日后管理，在上游壩坡設一馬道，高程為909 m，寬2 m；下游坡比1:2.75，設一級戧臺，高程為900.50 m，戧臺寬2 m，最大壩底寬度為200.7 m。

排水體設計：按照《小型水利水電工程碾壓式土石壩設計規(guī)范》(SL189-2013)6.9規(guī)定：壩體排水設施型式應綜合下列因素確定：(1)壩型及壩體、壩基材料的性質。(2)壩基的工程地質和水文地質條件。(3)下游水位。(4)排水設施的材料及施工條件。(5)壩址區(qū)的氣候條件。壩體排水形式一般采用棱體排水、貼坡排水、褥墊排水及綜合排水四種。經(jīng)綜合考慮因下游水位較高，貼坡排水體雖能防止發(fā)生滲透破壞但是不能有效的降低浸潤線。本次壩體排水采用棱體排水、棱體排水和褥墊排水綜合排水兩種型式進行比較后確定。

排水棱體頂部高程由下游最大水深和壩體浸潤線距壩面距離大于最大凍土深度兩種工況確定。本次設計因溢洪道下游校核水位為884.7 m，則下游水頭為1.7 m。壩體浸潤線出逸點高程為883.6 m，本地區(qū)最大凍土深度為0.98 m。結合現(xiàn)狀地形經(jīng)綜合確定壩體排水棱體頂高程為887.0 m。棱體內、外坡比均為1:1.5，頂部寬度2.0 m。

設計時對大壩進行滲流穩(wěn)定計算，選定棱體排水、棱體排水和褥墊排水綜合排水兩種型式計算大壩的滲流穩(wěn)定性，最終通過計算確定本次采用排水棱體與褥墊排水綜合排水形式。褥墊深入壩體的長度根據(jù)規(guī)范、施工難易程度和地質情況等綜合考慮確定褥墊深入壩體內70 m，褥墊采用砂礫料填筑，反濾層級配及厚度同排水棱體，褥墊厚度根據(jù)類似工程經(jīng)驗確定為50 cm厚。

2 模型建立

2.1 計算軟件介紹

計算采用河海大學滲流穩(wěn)定分析軟件AUTO BANK ，該軟件主要根據(jù)《碾壓土石壩設計規(guī)范》(SL 274-2001)、《土質邊坡穩(wěn)定分析-原理、方法、程序》(陳祖毓，中國水利水電出版社)等為編制依據(jù)，采用極限平衡法進行滲流穩(wěn)定計算，通過軟件進行建模，定義模型的材料參數(shù)，包括容重、滲透系數(shù)、水上和水下土體的物理力學指標等，再給定模型的水頭邊界條件進行滲流計算，通過定義模型的初始邊界條件，最終輸出結果。

2.2 計算原理

大壩滲流計算符合達西定律的非均各向異性二維滲流場，水頭勢函數(shù)滿足微分方程：

(1)

式中：φ=φ(x,y)為待求水頭勢函數(shù)；x，y為平面坐標；Kx，Ky為x，y軸方向的滲透系數(shù)。

水頭φ還必須滿足一定的邊界條件，經(jīng)常出現(xiàn)以下幾種邊界條件：

(1) 在上游邊界上水頭已知φ=φn

(2)

(2) 在逸出邊界水頭和位置高程相等φ=z

(3)

(3) 在某邊界上滲流量q已知

(4)

其中l(wèi)x，ly為邊界表面向外法線在x，y方向的余弦。

將滲流場用有限元離散，假定單元滲流場的水頭函數(shù)勢φ為多項式，由微分方程及邊界條件確定問題的變分形式，可導得出線性方程組：

[H]{φ}={F}

(5)

式中：[H]為滲透矩陣；{φ}為滲流場水頭；{F}為節(jié)點滲流量。

求解以上方程組可以得到節(jié)點水頭，據(jù)此求得單元的水力坡降，流速等物理量。求解滲流場的關鍵是確定浸潤線位置，Autobank采用節(jié)點流量平衡法通過迭代計算自動確定浸潤線位置和滲流量。

2.3 計算工況

根據(jù)《碾壓式土石壩設計規(guī)范》SL274-2001及本工程實際情況，確定大壩計算工況為：

(1)正常蓄水位穩(wěn)定滲流；

(2)正常蓄水位912.90 m降落至放水管底部高程896.00 m降落期非穩(wěn)定滲流。

2.4 計算材料參數(shù)

根據(jù)地質報告成果，確定均質土壩壩體、壩基各分區(qū)滲流計算采用參數(shù)指標見表1。

表1 滲流計算參數(shù)表

2.5 計算斷面選擇

滲流計算時選用最大典型斷面進行滲流計算(圖1)。

圖1 大壩標準橫斷面圖

3 計算結果及分析

根據(jù)大壩典型斷面，計算大壩的滲流特性(見圖2)。

圖2 有限元網(wǎng)格剖分圖

3.1 棱體排水下大壩滲流分析

工況一：正常蓄水位穩(wěn)定滲流(見圖3)

圖3 工況一水頭等值線和水力坡降圖

根據(jù)正常蓄水位下水頭等值線圖和水力坡降等值線圖可以看出，水力坡降在上游最大，下游排水棱體處最小，最大值為0.57，最小值為0.008 25。

工況二：正常蓄水位912.90 m降落至放水管底部高程896.00 m降落期非穩(wěn)定滲流(見圖4)

圖4 工況二水頭等值線和水力坡降圖

根據(jù)水位降落期水頭等值線圖和水力坡降等值線圖可以看出，水力坡降在上游水位變動區(qū)最大，下游排水棱體處最小，最大值為0.382，最小值為0.054 7。

3.2 棱體排水和褥墊排水下大壩滲流分析

工況一：正常蓄水位穩(wěn)定滲流(見圖5)

圖5 工況一水頭等值線和水力坡降圖

根據(jù)正常蓄水位下水頭等值線圖和水力坡降等值線圖可以看出，水力坡降在上游最大，下游排水棱體處最小，最大值為0.772，最小值為0.089 7。

工況二：正常蓄水位912.90 m降落至放水管底部高程896.00 m降落期非穩(wěn)定滲流(見圖6)

圖6 工況二水頭等值線和水力坡降圖

根據(jù)水位降落期水頭等值線圖和水力坡降等值線圖可以看出，水力坡降在上游水位變動區(qū)最大，下游排水棱體處最小，最大值為0.545，最小值為0.061 4。

3.3 計算結果分析

兩種排水型式下計算正常蓄水位工況和水位降落期工況下大壩的浸潤線和水力坡降，可以看出，僅設排水棱體時，大壩在下游壩坡浸潤線下降比較慢，而排水棱體和褥墊排水時，壩體浸潤線在褥墊排水處呈明顯下降趨勢，壩內水能夠快速的從褥墊排出壩外，增強了大壩的滲透穩(wěn)定性。同時從等值線圖可以看出，增設褥墊排水后，等值線值在該處明顯變小，說明褥墊增設后水力坡降變小，有利于壩內水的排出。最終確定袁家溝水庫壩坡排水型式采用棱體排水和褥墊排水相結合的綜合排水型式。

4 結語

通過進行建模計算，對袁家溝水庫工程確定合理的壩坡排水型式，最終進行比較后確定采用棱體排水和褥墊排水的綜合排水型式，在實際工程中，有多種綜合排水型式，若遇下游水位較高時，可采用棱體排水和貼坡排水相結合的排水型式；若遇壩基透水性較強可選用本次計算分析的棱體排水和褥墊排水相結合的排水型式。因此，排水型式的選擇，必須結合壩基排水的需要和形式、針對不同的地基和外界因素，選擇合適的壩坡排水型式。