宋 昊

（上海市崇明區(qū)農(nóng)村水利管理所,上海 202150）

1 引言

土質(zhì)防滲斜心墻壩是我國20 世紀(jì)中期發(fā)展起來的一種土石壩壩型,隨后在我國小浪底等土石壩中廣泛采用這種壩型[1]。相比傳統(tǒng)的心墻壩,土質(zhì)防滲斜心墻壩具有壩坡穩(wěn)定性強,上游防滲效果好等優(yōu)點[2]。為評價土石壩的滲透性能,一般要進行穩(wěn)定-非穩(wěn)定滲流分析,進而對土石壩的滲透穩(wěn)定性進行評估。通常相關(guān)學(xué)者的評價指標(biāo)包括滲流量與出逸點滲透坡降。然而,傳統(tǒng)土石壩滲流計算分析,一般主要針對均質(zhì)土壩、心墻壩、混凝土面板堆石壩的滲流計算分析較多,對黏土斜心墻壩穩(wěn)定—非穩(wěn)定滲流研究報道較少[3]。因此,本文以某壤土斜心墻壩為工程背景,開展?jié)B透穩(wěn)定計算,基于計算指標(biāo)來評價滲透穩(wěn)定及滲漏損失是否滿足安全要求。

2 工程概況

某水庫為一灌溉、供水和防洪的綜合利用水利樞紐工程,水庫設(shè)計總庫容為528 萬m3。擋水建筑物由非過水土壩和過水土壩兩壩段組成,3 級建筑物。設(shè)計標(biāo)準(zhǔn)為百年一遇洪水設(shè)計、千年一遇洪水校核。土壩全長163.5 m,其中過水土壩壩段位于河床部位,長72.0 m,溢流前沿凈寬50.0 m。左右兩岸為非過水土壩壩段,累計長度為91.5 m。最大壩高36.8 m。土壩防滲墻采用壤土斜心墻,上游側(cè)邊坡為1∶3,下游側(cè)邊坡1∶2.55。防滲斜墻頂部最小有效厚度為1.5 m,底部有效厚度為4.5 m。土壩壩殼由砂礫料填筑,壤土斜心墻上游側(cè)砂礫料壩殼厚1.2 m,上游壩坡與壤土斜心墻上游側(cè)邊坡相同,為1∶3。下游側(cè)砂礫料壩殼邊坡非過水土壩壩段為1∶2,過水土壩壩段為1∶3。

3 計算分析方法

3.1 計算原理

（1）滲流計算原理

符合達西定律的二維非均質(zhì)各向異性土體非穩(wěn)定滲流的基本方程為：

穩(wěn)定滲流的基本方程為：

式中：h為水頭函數(shù)；x、y為計算坐標(biāo)；t為時間；kx、ky為以x、y軸為主軸方向的滲透系數(shù)；為給水度。

式（1）～式（2）即為二向穩(wěn)定和穩(wěn)定滲流計算的基本微分方程和定解條件,聯(lián)立求解即可確定土體的穩(wěn)定滲流場和非穩(wěn)定滲流場。

（2）滲流壓力計算原理

根據(jù)滲流計算的流網(wǎng)可以得到水頭損失的沿程變化,繪制出沿地下輪廓的滲流壓力分布圖。

（3）滲透坡降計算原理

根據(jù)《水利水電工程地質(zhì)勘察規(guī)范》流土型的滲透變形的判別可采用以下方法：

式中：Jcr為土的臨界水力比降；Gs為土粒比重；n為土的孔隙率。

設(shè)計時采用的容許水力比降等于臨界水力比降除以安全系數(shù)：安全系數(shù)一般取1.5～2.0；當(dāng)滲流穩(wěn)定對水工建筑物的危害較大時,取2.0 的安全系數(shù)；對特別重要的工程也可用2.5 的安全系數(shù)。

3.2 計算模型

本論文中采用有限元法對大壩進行滲流計算,運用Geo-Stvdio軟件中seep/w模塊,建立大壩現(xiàn)狀最大橫斷面的二維有限元模型。根據(jù)大壩現(xiàn)狀最大橫斷面（樁號0+100）建立有限元模型,模型范圍為：大壩向上游與下游各取40 m,建基面以下取40 m。采用四邊形單元為主、輔以三角形單元來進行有限元網(wǎng)格剖分,共包括5321 個單元,5409 個節(jié)點,得到的有限元網(wǎng)格見圖1。

圖1 大壩滲流計算有限元網(wǎng)格劃分

3.3 計算參數(shù)

壩殼砂礫料、壩基砂礫石、壩體防滲斜墻、防滲心墻、壩體反濾層、壩基片麻巖、回填碎石渣及堆石排水體等滲流計算材料參數(shù)見表1。

表1 壩體和壩基滲流計算參數(shù)表

3.4 計算工況

參考《碾壓式土石壩設(shè)計規(guī)范》（SL 274-2001）規(guī)定,滲流計算應(yīng)考慮運行中出現(xiàn)的不利工況條件,按穩(wěn)定滲流考慮,并考慮以下滲流計算工況：

（1）正常蓄水位工況：上游正常蓄水位288.80 m,下游對應(yīng)水位取248.00 m；

（2）設(shè)計洪水位工況：上游設(shè)計洪水位290.30 m,下游對應(yīng)水位取248.00 m；

（3）校核洪水位工況：上游校核洪水位291.60 m,下游對應(yīng)水位取248.00 m；

（4）水位驟降工況Ⅰ：上游水位按水位降落速度0.5 m/d由正常蓄水位288.80 m驟降至死水位265.00 m,下游對應(yīng)水位取248.00 m；

（5）水位驟降工況Ⅱ：上游水位按水位降落速度1.0 m/d由正常蓄水位288.80 m驟降至死水位265.00 m,下游對應(yīng)水位取248.00 m。

4 計算結(jié)果分析

本次正常蓄水位工況、設(shè)計洪水位工況及校核洪水位工況下滲流計算的滲透坡降、滲流量及壩體出逸點計算成果見表2,計算的浸潤線位置及等勢線圖見圖2～圖4。由表2 及圖2～圖4 可以看出,在上下游水位差的作用下,滲流區(qū)域主要集中在斜心墻區(qū)域,沒有出現(xiàn)明顯的集中滲流。由于防滲墻的作用,壩體內(nèi)部的浸潤線大致沿著斜墻部位快速下降,壩后浸潤線很低,正常蓄水位、設(shè)計洪水位和校核洪水位工況的壩體出逸點高程變化正常,隨著上游水位的升高略有增大。同時,計算結(jié)果表明,隨著上游水位增大,正常蓄水位工況、設(shè)計洪水位工況及校核洪水位工況下壩體滲流量分別為29.91 m3/（d·m）、30.05 m3/（d·m）、30.95 m3/（d·m）,表現(xiàn)為逐漸增大的趨勢。此外,壩體出逸段的滲透坡降也隨著上下游的水位差逐漸增大,壩體出逸面是壩基片麻巖,故壩體出逸點滲透穩(wěn)定。壩體及壩基與堆石排水體接觸面的滲透坡降小于允許值,故壩體及壩基滲透穩(wěn)定。

圖2 正常蓄水位工況浸潤線位置及等勢線圖

圖3 設(shè)計洪水位工況浸潤線位置及等勢線圖

圖4 校核洪水位工況浸潤線位置及等勢線圖

表2 各工況滲流分析結(jié)果

本次水位驟降工況Ⅰ、Ⅱ工況下滲流計算的滲透坡降、滲流量及壩體出逸點計算成果見表2,取計算時間段任意時刻的浸潤線位置及等勢線圖見圖5～圖6。由表2及圖5～圖6可以看出：降落水位即為浸潤線與上游壩坡交點的高程。不同時段上游壩體內(nèi)的浸潤線基本與水位降落速度同步,基本均呈現(xiàn)為“Z”字形斜線；相比較而言,水位驟降工況Ⅰ的上游壩體內(nèi)浸潤線高于水位驟降工況Ⅱ的浸潤線,且隨著降落水位的逐步降低,兩者浸潤線的高差相應(yīng)增大。當(dāng)降至死水位時,驟降工況Ⅰ和工況Ⅱ的浸潤線相同。根據(jù)以上分析可以得出,隨著上游水位降落速度由0.5 m/d增加至1.0 m/d,壩體滲流量減小36.8%,滲透坡降基本一致,均滿足滲透穩(wěn)定要求。

圖5 水位驟降工況Ⅰ等勢線及浸潤線位置圖

圖6 水位驟降工況Ⅱ等勢線及浸潤線位置圖

5 結(jié)論

（1）通過分析不同工況下的水頭等值線、浸潤線、滲流量等結(jié)果,發(fā)現(xiàn)正常蓄水位、設(shè)計洪水位及校核洪水位工況下隨著上下游水頭差的增大,滲流量與出逸點滲透坡降均有所增大。水位驟降工況下,隨著上游水位降落速度增加,壩體滲流量增大,滲透坡降變化不大。

（2）由于斜墻壩壩體內(nèi)部設(shè)置了壤土斜墻,有效地降低了浸潤線,各工況下壩體浸潤線和出逸點的位置正常,壩體滿足滲透穩(wěn)定要求。