蘇興寬

（吉林省第二松花江防汛機(jī)動(dòng)搶險(xiǎn)隊(duì)，吉林 吉林 132101）

我國(guó)是世界上水庫(kù)數(shù)量最多的國(guó)家，同時(shí)我國(guó)水庫(kù)的病險(xiǎn)問(wèn)題也日益凸顯［1］。我國(guó)病險(xiǎn)水庫(kù)不僅數(shù)量龐大，而且分布廣泛［2］。在過(guò)去的“十五”期間，我國(guó)加大了對(duì)大中型病險(xiǎn)水庫(kù)除險(xiǎn)加固的投資力度，大型水庫(kù)的病險(xiǎn)率已由1999年底的42%下降到14%，中型水庫(kù)病險(xiǎn)率由41%下降到25%，取得了較好效果［3］。而根據(jù)統(tǒng)計(jì)，我國(guó)中小型病險(xiǎn)水庫(kù)數(shù)量達(dá)5.74萬(wàn)座，占總數(shù)量的95.3%。因此，小型病險(xiǎn)水庫(kù)除險(xiǎn)加固是當(dāng)前我國(guó)病險(xiǎn)水庫(kù)除險(xiǎn)加固工作的重點(diǎn)，而準(zhǔn)確定量地評(píng)估加固后壩體的滲透穩(wěn)定更是水庫(kù)除險(xiǎn)加固的工作重點(diǎn)［4-6］。

本文以腰石水庫(kù)為例，首先分析了大壩目前存在的問(wèn)題，并提出相應(yīng)的整治建議，研究了該病險(xiǎn)水庫(kù)均質(zhì)壩在加固后的滲流穩(wěn)定情況，重點(diǎn)分析不同工況下壩體滲流場(chǎng)、滲透比降、滲流量、壩坡安全系數(shù)和最不利滑裂面，從而復(fù)核大壩整治設(shè)計(jì)方案的合理性，并為后期壩體的運(yùn)行管理提供借鑒。

1 工程概述

腰石水庫(kù)位于吉林省磐石市煙筒山鎮(zhèn)腰石屯南0.5km，飲馬河支流大力河上游，控制流域面積8.3km2，是一座以防洪、灌溉為主，結(jié)合養(yǎng)魚(yú)等綜合利用的小型水庫(kù)。該水庫(kù)是以灌溉為主，兼有防洪等綜合利用的?。?）型水利工程，大壩等主要建筑物為5級(jí)建筑物，次要建筑物與臨時(shí)建筑物均為5級(jí)。該工程壩址以上集雨面積8.3km2，總庫(kù)容33.64萬(wàn)m3，有效庫(kù)容25.65萬(wàn)m3，設(shè)計(jì)控制灌溉面積26.67hm2。樞紐工程由大壩、溢洪道、取水設(shè)施3部分組成。大壩為均質(zhì)土壩，壩頂高程384.50m，最大壩高12m，壩頂寬2.5m，壩頂長(zhǎng)100m。大壩內(nèi)坡為六棱塊襯砌，雜草叢生；大壩壩頂未硬化，下游未襯護(hù)、坡面不規(guī)整。

2 大壩問(wèn)題分析及整治建議

2.1 現(xiàn)狀調(diào)查

通過(guò)對(duì)水庫(kù)大壩進(jìn)行現(xiàn)場(chǎng)安全檢查，并對(duì)大壩工程質(zhì)量、運(yùn)行管理、結(jié)構(gòu)安全、滲流安全進(jìn)行分析評(píng)價(jià)，對(duì)大壩防洪標(biāo)準(zhǔn)、抗震安全進(jìn)行復(fù)核。發(fā)現(xiàn)水庫(kù)大壩主要存在以下問(wèn)題：①大壩上下游壩坡雜草叢生，下游坡面不平整，無(wú)周邊排水溝；②大壩下游壩坡穩(wěn)定安全系數(shù)不滿足SL189—2013《小型水利水電工程碾壓式土石壩設(shè)計(jì)規(guī)范》相關(guān)要求。

2.2 整治建議

針對(duì)大壩下游壩坡不穩(wěn)的情況，可采用兩種方案進(jìn)行處理。

（1）方案1：加寬硬化壩頂，將壩軸線下移，培厚下游壩體。

（2）方案2：對(duì)下游壩體進(jìn)行削坡處理，放緩下游壩坡。

方案2將進(jìn)一步減少壩頂寬度，對(duì)壩頂結(jié)構(gòu)及交通要求不利。故本工程最終采用方案1為大壩加固設(shè)計(jì)推薦方案。

本次考慮對(duì)下游壩坡培厚至1∶1.9，并采用框格梁草皮護(hù)坡進(jìn)行加固處理；加寬壩頂路面至4.5m，重建防浪墻，拆除原下游排水棱體，并新建排水棱體，增設(shè)排水溝和下游坡腳混凝土擋墻整治。

3 研究方法

3.1 非飽和滲流有限元理論

根據(jù)達(dá)西滲透定律，巖土體非穩(wěn)定滲流場(chǎng)的控制性微分方程可表述為：

式中 ?x， ?y分別為x，y軸方向；kx，ky分別為x，y方向的滲透系數(shù)（m/s）；H為水頭（m）；Q為邊界滲漏量（m3/s）；Mw為水土特征曲線的斜率；γw為水重度，取1000kg/m3。進(jìn)一步采用（Glerkin）加權(quán)余量法，則可推導(dǎo)出有限元滲流方程［7-9］，如式（2）：

式中 τ為單元厚度（m）；A為單元面積（m2）；λ為Mwγw；B為水力梯度矩陣；C為單元導(dǎo)水率矩陣；H為水頭向量；N為插值函數(shù)向量；q為單元邊界單位滲流量（m3/s）；L為單元邊長(zhǎng)（m）。

3.2 塊體極限平衡法

壩體的穩(wěn)定性計(jì)算主要考慮水位以上坡體的滲流壓力，而水位以下部分，將坡面以上水體視為一種有重度無(wú)強(qiáng)度的特殊材料，根據(jù)該思路本次壩坡穩(wěn)定計(jì)算采用Bishop簡(jiǎn)化條分法［10］。

結(jié)合大壩工程地質(zhì)條件，采用GeoStudio軟件計(jì)算分析庫(kù)水位在不同水位時(shí)壩體內(nèi)瞬時(shí)浸潤(rùn)線，并將計(jì)算結(jié)果導(dǎo)入邊坡穩(wěn)定性計(jì)算模塊中，運(yùn)用塊體極限平衡法進(jìn)行穩(wěn)定性計(jì)算，即可獲得大壩壩坡穩(wěn)定安全系數(shù)與水位變動(dòng)速率條件、水位高度與時(shí)間的相互關(guān)系［11］。結(jié)合現(xiàn)場(chǎng)取樣及室內(nèi)土體試驗(yàn)報(bào)告，壩體二維滲流計(jì)算主要物理力學(xué)參數(shù)如表1。

表1 壩體材料主要物理力學(xué)參數(shù)

3.3 計(jì)算工況

3.3.1 滲流分析計(jì)算工況

根據(jù)SL189—2013《小型水利水電工程碾壓式土石壩設(shè)計(jì)規(guī)范》規(guī)定，滲流計(jì)算應(yīng)考慮水庫(kù)運(yùn)行中的各種不利組合，并結(jié)合水庫(kù)的水位特征和運(yùn)行情況，對(duì)各種水位情況下做了穩(wěn)定滲流計(jì)算和非穩(wěn)定滲流計(jì)算。其滲流計(jì)算工況如下：

（1）上游正常蓄水位與相應(yīng)下游水位。（2）上游設(shè)計(jì)洪水位與相應(yīng)下游水位。（3）上游校核洪水位與相應(yīng)下游水位。（4）上游校核洪水位降至正常水位。

（5）上游正常蓄水位降至死水位。

由于大壩下游為養(yǎng)殖魚(yú)塘，下游水位高程均采用373.10m。大壩正常蓄水位382.00m，設(shè)計(jì)洪水位382.98m，校核水位383.44m，死水位375.66m。

3.3.2 壩坡穩(wěn)定分析計(jì)算工況

本工程區(qū)地震烈度為6°，不考慮地震。計(jì)算工況包括正常運(yùn)用條件及非常運(yùn)用條件。各計(jì)算工況如下：

3.3.2.1 正常運(yùn)用條件

（1）工況1：由正常蓄水位形成穩(wěn)定滲流期的上、下游壩坡。

（2）工況2：設(shè)計(jì)洪水位形成穩(wěn)定滲流期的上、下游壩坡。

（3）工況3：校核水位形成穩(wěn)定滲流期的上、下游壩坡。

3.3.2.2 非正常運(yùn)用條件

（1）工況4：由校核水位驟降至正常蓄水位形成非穩(wěn)定滲流期的上游壩坡。

（2）工況5：由正常蓄水位驟降至死水位形成非穩(wěn)定滲流期的上游壩坡。

根據(jù)實(shí)際臥管的最大放水速率，選取庫(kù)水位驟降速率為1.5m/d。

4 大壩滲流穩(wěn)定分析

4.1 滲流

圖1為整治后的腰石水庫(kù)大壩在5種運(yùn)行工況下的浸潤(rùn)線分布。分析可知：在正常蓄水位、設(shè)計(jì)洪水位和校核水位3種工況下，壩體形成穩(wěn)定滲流，浸潤(rùn)線均通過(guò)排水棱體迅速降低，并排至下游魚(yú)塘。在2種驟降工況下，由校核水位驟降至正常蓄水位形成非穩(wěn)定滲流時(shí)，各時(shí)間段內(nèi)浸潤(rùn)線均在上游位置最高，說(shuō)明此時(shí)壩體內(nèi)“退水”較快，與水位下降基本保持同步，而在正常蓄水位驟降至死水位形成非穩(wěn)定滲流時(shí)，在驟降過(guò)程中壩體內(nèi)的浸潤(rùn)線明顯高于上游水位，說(shuō)明此時(shí)壩體內(nèi)“退水”慢，壩體內(nèi)水相對(duì)于庫(kù)水位下降更慢，上游部位的壩體土處于非飽和狀態(tài)，孔隙水壓力較大，根據(jù)有效應(yīng)力原理可知，這將對(duì)上游壩坡穩(wěn)定產(chǎn)生不利影響。

圖1 水庫(kù)大壩各工況壩體浸潤(rùn)線分布

表2為整治后的腰石水庫(kù)大壩不同運(yùn)行工況下的壩體滲漏量及滲透比降，分析可知：上游水位越高，壩體單寬滲流量越大（校核洪水位時(shí)為5.32×10-7m3/s），滲透比降也越大，在校核洪水位時(shí)滲透比降最大為0.51，但均小于試驗(yàn)允許值0.62，故認(rèn)為整治后的腰石大壩不會(huì)發(fā)生滲透破壞。

表2 滲流量及滲透比降計(jì)算成果

4.2 壩坡穩(wěn)定分析

根據(jù)上述對(duì)不同工況下大壩體的滲流計(jì)算結(jié)果，進(jìn)一步計(jì)算了對(duì)應(yīng)工況的上、下游壩坡穩(wěn)定安全系數(shù)及不利滑裂面，如圖2。分析可知，5種工況下大壩的上、下游最不利滑裂面基本相同，其中上游壩坡滑裂面起于下游坡面頂部，從上游坡腳剪出，而下游壩坡滑裂面起于上游坡面上部，從下游坡腳剪出。

圖2 水庫(kù)大壩各工況壩體抗滑穩(wěn)定計(jì)算結(jié)果

從表3可知，在壩體形成穩(wěn)定滲流時(shí)，隨著庫(kù)水位的升高，上游壩坡安全系數(shù)逐漸增大，而下游壩坡有逐漸減小的趨勢(shì)。當(dāng)壩體在非穩(wěn)定滲流工況（庫(kù)水位驟降）時(shí)，上游壩坡安全系數(shù)出現(xiàn)驟降，其中校核水位驟降至正常蓄水位（工況4）安全系數(shù)為1.592，而正常蓄水位驟降至死水位時(shí) （工況5）安全系數(shù)僅1.123，說(shuō)明此時(shí)上游壩坡穩(wěn)定性最差。但各工況下大壩上、下游壩坡抗滑穩(wěn)定最小安全系數(shù)均大于規(guī)范最低要求值，說(shuō)明整治后的大壩壩坡穩(wěn)定滿足規(guī)范要求。

表3 大壩壩坡抗滑穩(wěn)定安全系數(shù)統(tǒng)計(jì)

5 結(jié)語(yǔ)

（1）在正常蓄水位驟降至死水位形成非穩(wěn)定滲流時(shí)，壩體內(nèi)的浸潤(rùn)線明顯高于上游水位，上游部位的壩體土處于非飽和狀態(tài)，對(duì)上游壩坡穩(wěn)定不利。

（2）各工況下大壩的滲透比降均小于允許值0.62，整治后的大壩不會(huì)發(fā)生滲透破壞。

（3）當(dāng)壩體形成穩(wěn)定滲流時(shí)，庫(kù)水位升高使上游壩坡安全系數(shù)明顯增大，同時(shí)下游壩坡安全系數(shù)逐漸減小。

（4）不同水位下降速率時(shí)壩體的最不利滑裂面基本相同，庫(kù)水位下降速率和上游水位高程對(duì)大壩上、下游壩坡滑裂面的分布沒(méi)有影響。

（5）正常蓄水位驟降至死水位時(shí)上游壩坡安全系數(shù)僅為1.123，此時(shí)壩坡穩(wěn)定性最差。各工況下大壩壩坡抗滑穩(wěn)定最小安全系數(shù)均大于規(guī)范最低要求值，說(shuō)明整治后的大壩壩坡穩(wěn)定滿足規(guī)范要求。

［1］楊杰，鄭成成，江德軍，等.病險(xiǎn)水庫(kù)理論分析研究進(jìn)展［J］.水科學(xué)進(jìn)展， 2014， 25（1）：148-154.

［2］溫鵬，童學(xué)衛(wèi)，趙春，等.全國(guó)小型病險(xiǎn)水庫(kù)除險(xiǎn)加固進(jìn)度管理系統(tǒng)研究［J］.中國(guó)水利水電科學(xué)研究院學(xué)報(bào)，2016，14（5）：334-339.

［3］盛金保，劉嘉炘，張士辰，等.病險(xiǎn)水庫(kù)除險(xiǎn)加固項(xiàng)目潰壩機(jī)理調(diào)查分析［J］.巖土工程學(xué)報(bào)，2008， 30（11）：1620-1625.

［4］時(shí)鐵城，阮建飛，張曉.庫(kù)水位驟降情況下土石壩壩坡穩(wěn)定分析［J］.人民黃河， 2014（2）：93-94.

［5］羅春，羅園，王鋒.庫(kù)水位驟降條件下均質(zhì)土石壩滲流穩(wěn)定分析［J］.江西水利科技， 2016，42（3）：202-207.

［6］尹吉娜，柴瑞，楊杰，等.庫(kù)水位驟降下除險(xiǎn)加固后土壩安全性態(tài)響應(yīng)及評(píng)價(jià) ［J］.水電能源科學(xué)，2015 （7）：76-79.

［7］周躍峰，譚國(guó)煥，甄偉文，等.非飽和滲流分析在FLAC3D中的實(shí)現(xiàn)和應(yīng)用［J］.長(zhǎng)江科學(xué)院院報(bào)， 2013， 30（2）：57-61.

［8］Huang M， Jia C Q.Strength reduction FEM in stability analysis of soil slopes subjected to transient unsaturated seepage［J］.Computers&Geotechnics， 2007， 36（1）：93-101.

［9］Huang M S， Jia C Q.Stability analysis of soil slopes subjected to unsaturated transient seepage ［J］.Chinese Journal of Geotechnical Engineering， 2006， 28 （2）：202-206.

［10］Chen X，Wu Y， Yu Y， et al.A two-grid search scheme for large-scale 3-D finite element analyses of slope stability［J］.Computers&Geotechnics， 2014， 62（62）：203-215.

［11］李海洋.基于Auto-bank軟件的閘基滲流分析［D］.青島：中國(guó)海洋大學(xué)，2015.