周子宜

(廣東河海工程咨詢有限公司，廣東 茂名 525000)

1 工程概況

雞姆塘水庫位于高州市鎮(zhèn)江鎮(zhèn)，是一個以農(nóng)業(yè)灌溉為主，兼顧防洪功能的小(2)型水庫工程，總庫容為13.08萬m3。

雞姆塘水庫現(xiàn)狀主要建筑物包括大壩、溢洪道、輸水涵等，其中，大壩為均質(zhì)土壩，溢洪道為正槽開敞式溢洪道，輸水涵管為壩下涵管形式。主要建筑物級別均為5級，次要建筑物級別為5級。溢洪道位于大壩壩肩右側(cè)，為開敞式溢洪道，無閘控制，過流寬度為2.20 m，兩側(cè)擋墻為磚砌結(jié)構(gòu)，高1.30 m。輸水涵管位于大壩左側(cè)，采用斜臥管式進(jìn)水口，壩下涵管為Φ800鋼筋混凝土涵管，直徑為0.80 m，長度為26 m。雞姆塘水庫大壩經(jīng)除險加固后，大壩總長115 m，最大壩高為8 m，壩頂寬度約5 m。

雞姆塘水庫大壩由素填土組成。壩基主要由粉質(zhì)黏土層和全風(fēng)化變質(zhì)砂巖層組成。粉質(zhì)黏土層厚度2.9～3.4 m，層底高程為30.37～30.48 m，基礎(chǔ)承載力特征值為130 kPa，天然含水量為29.2%～31.3%，平均值為30.3%；孔隙比0.89～0.95；塑性指數(shù)為13.40～13.87，平均值為13.64。全風(fēng)化變質(zhì)砂巖層為磚紅色、黃褐色，全風(fēng)化砂土狀，巖體風(fēng)化劇烈，巖芯多呈散體狀，標(biāo)貫擊數(shù)為34擊/30 cm。該層揭露厚度2.2～2.9 m，層底高程為24.38～24.67 m。

本文對經(jīng)過除險加固后的雞姆塘水庫大壩進(jìn)行滲流與壩坡穩(wěn)定計算，具體而言，選取雞姆塘水庫大壩典型斷面，采用二維滲流有限元法進(jìn)行滲流穩(wěn)定計算分析，并采用剛體極限平衡法中的簡化畢肖普法對壩坡進(jìn)行穩(wěn)定計算分析。雞姆塘水庫大壩典型斷面如圖1所示。

2 大壩壩頂高程復(fù)核

雞姆塘水庫大壩等級為5級，根據(jù)《小型水利水電工程碾壓式土石壩設(shè)計規(guī)范》(SL 189—2013)[1]，壩頂在水庫靜水位以上的超高計算如式(1)：

Y=R+A

(1)

式中：Y為壩頂超高，m；R為波浪沿著壩坡的最大爬高，m；A為安全加高，m，對于 4級、5級壩，設(shè)計工況取0.5 m，校核工況取0.3 m。

2.1 計算參數(shù)的選取

(1)風(fēng)速W。庫區(qū)無實(shí)測氣象資料，采用高州市多年平均年最大風(fēng)速為14.6 m/s。根據(jù)規(guī)范SL 189—2013規(guī)定：正常運(yùn)用條件下，采用多年平均最大風(fēng)速的1.5倍，即21.9 m/s。非常運(yùn)用條件下，采用多年平均最大風(fēng)速14.6 m/s。

(2)吹程D為80 m。

(3)水域平均深度Hm。設(shè)計水位條件下，水庫靜水位40.05 m，水域平均深度6.84 m；校核水位條件下，水庫靜水位40.25 m，水域平均深度7.03 m。

(4)波浪爬高R。根據(jù)規(guī)范SL 189—2013，對于4級、5級壩，應(yīng)先計算平均波浪爬高值Rm，再根據(jù)平均波高與壩前迎水面水深的比值(hm/H)，查得Rp/Rm比值，從而求得累計頻率爬高值Rp。雞姆塘水庫按正向來波單坡考慮，迎水坡坡度m為2.25，平均波浪爬高Rm計算如式(2)：

(2)

式中：KΔ為斜坡糙率滲透性系數(shù)，取0.90；KW為經(jīng)驗(yàn)系數(shù)；Lm為平均波長，m。其中，平均波高h(yuǎn)m、平均波長Lm采用莆田試驗(yàn)站公式計算[2]。

2.2 計算結(jié)果

壩頂高程等于水庫靜水位與超高之和，并按不同運(yùn)行條件計算，取計算結(jié)果的較大值，計算結(jié)果如表1 所示。水庫大壩工程所在地區(qū)地震基本烈度為Ⅶ度，地震安全加高值取0.6 m。

表1 大壩壩頂高程計算結(jié)果

由表1可知，對水庫大壩三種工況進(jìn)行壩頂高程復(fù)核，水庫設(shè)計工況時，水庫壩頂高程為41.01 m，校核工況時，壩頂高程為40.79 m，地震工況時，壩頂高程為41.03 m。雞姆塘水庫大壩經(jīng)除險加固后的壩頂高程為41.10 m，壩頂高程復(fù)核滿足要求。

3 大壩滲流穩(wěn)定計算

3.1 滲流計算工況

按規(guī)范(SL 274—2020)[3]規(guī)定，對以下水位組合的工況進(jìn)行計算：

工況一：上游正常蓄水位39.55 m時形成的穩(wěn)定滲流情況。

工況二：上游設(shè)計洪水位40.05 m時形成的穩(wěn)定滲流情況。

工況三：上游校核洪水位40.24 m時形成的穩(wěn)定滲流情況。

工況四：正常蓄水位39.55 m降至死水位35.30 m時形成的非穩(wěn)定滲流情況。

3.2 地質(zhì)參數(shù)

大壩典型斷面各區(qū)巖土層的滲透系數(shù)根據(jù)相關(guān)勘測資料進(jìn)行取值，大壩滲流穩(wěn)定計算所選用的地質(zhì)參數(shù)見表2。

表2 各土層主要物理力學(xué)指標(biāo)建議值

3.3 滲流計算方法

采用二維滲流有限元法進(jìn)行滲流計算分析，計算程序?yàn)槔碚龓r土中的有限元滲流計算軟件，計算方法為有限元法，主要依據(jù)是非飽和土理論[4]、達(dá)西定律[5]等。

2.3.1 二維滲流有限元法原理

對于穩(wěn)定滲流，符合達(dá)西定律的非均勻各向異性二維滲流場，水頭勢函數(shù)滿足微分方程，如式(3)[6]：

(3)

式中：φ=φ(x，y)為待求水頭勢函數(shù)；x，y為平面坐標(biāo)；kx，ky分別為x，y方向的滲透系數(shù)。水頭φ還需要滿足一定的邊界條件，一般包括以下幾種邊界條件：

(1)上游邊界水頭已知φ0的情況如式(4)：

φ=φ0

(4)

(2)逸出邊界水頭和位置高程z已知的情況如式(5)：

φ=z

(5)

(3)某邊界上滲流流量q已知的情況如式(6)：

(6)

式中：lx和ly為邊界表面向外法線在下x，y方向的余弦。

將滲流場用有限元離散，假定單元滲流場的水頭勢函數(shù)φ為多項(xiàng)式，由微分方程及邊界條件確定問題的變分形式，可得到穩(wěn)定滲流有限元計算如式(7)所示：

[H]{φ}={F}

(7)

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

求解以上方程組便可以得到節(jié)點(diǎn)水頭，在此基礎(chǔ)上求得單元的水力坡降、流速等物理量。求解滲流場的關(guān)鍵步驟是確定浸潤線的位置。本文計算軟件采用的是理正巖土中的有限元滲流計算軟件，其用改進(jìn)平方根法直接求解以上滲流計算方程組，從而可以確定多種邊界條件下的大壩浸潤線位置。

3.4 大壩滲流穩(wěn)定計算成果

根據(jù)擬定的邊界條件對雞姆塘水庫大壩進(jìn)行滲流分析計算，計算結(jié)果見表3。

表3 雞姆塘水庫大壩滲流計算成果

《水閘設(shè)計規(guī)范》(SL 265—2016)[7]規(guī)定：當(dāng)滲流出口處設(shè)有濾層時，允許滲透比降按土的允許滲透比降可增大30%。由滲流計算結(jié)果可知：各個工況下，雞姆塘水庫大壩的滲流出逸位置均在下游排水棱體內(nèi)，出逸比降均小于設(shè)有反濾層時的允許水力比降，因此，雞姆塘水庫大壩的滲流是安全的。

4 大壩壩坡穩(wěn)定計算

4.1 壩坡穩(wěn)定計算工況

雞姆塘水庫大壩壩坡穩(wěn)定計算斷面與滲流穩(wěn)定計算的斷面一致。根據(jù)規(guī)范SL 189—2013[7]，結(jié)合該水庫運(yùn)行中實(shí)際可能出現(xiàn)的不利情況，擬定計算工況如表4所示：

表4 雞姆塘水庫大壩壩坡穩(wěn)定計算工況

4.2 壩坡穩(wěn)定計算方法

計算采用剛體極限平衡法中的簡化畢肖普法，計算軟件采用理正巖土系列軟件中的“邊坡穩(wěn)定分析軟件”。

瑞典條分法[8]和簡化畢肖普法[9]是在巖土工程領(lǐng)域計算邊坡穩(wěn)定系數(shù)應(yīng)用最為廣泛的兩種方法，都屬于極限平衡法的一種。其中，簡化畢肖普法是在邊坡土體中取一圓弧滑動體，當(dāng)土體內(nèi)部所受剪應(yīng)力大于其最大抗剪應(yīng)力時，土體就會沿著此圓弧滑動面發(fā)生破壞。

4.3 大壩壩坡穩(wěn)定計算成果

各種工況下，大壩的邊坡穩(wěn)定最小安全系數(shù)計算結(jié)果見表5。

表5 雞姆塘水庫大壩壩坡穩(wěn)定計算結(jié)果

由表5可知：雞姆塘水庫大壩壩坡穩(wěn)定最小安全系數(shù)出現(xiàn)在考慮地震作用，正常蓄水位時形成穩(wěn)定滲流的情況下，最小邊坡穩(wěn)定安全系數(shù)為1.24，小于規(guī)范規(guī)定的允許值1.25，因此，雞姆塘水庫大壩壩坡在正常蓄水位、設(shè)計洪水位、校核洪水位以及考慮地震作用下的正常蓄水位形成穩(wěn)定滲流的情況下，其上下游壩坡均能保持穩(wěn)定[10]。

5 結(jié) 語

對經(jīng)除險加固后的雞姆塘水庫大壩進(jìn)行滲流和壩坡穩(wěn)定計算分析，評估其能否安全穩(wěn)定運(yùn)行，經(jīng)分析，雞姆塘水庫大壩滿足滲流和壩坡穩(wěn)定要求，安全穩(wěn)定。