楊 娜

(河北省子牙河河務(wù)中心，河北 衡水 053000)

0 前 言

土石混合體壩基在水壩工程中常是設(shè)計首選結(jié)構(gòu)形式，但是有關(guān)土石混合體壩基的滲透破壞特點，相關(guān)研究成果相對較少。胡亞元等才用有限元方法研究了管涌作用對堤壩滲透系數(shù)和滲透速度的影響規(guī)律，主要基于改進后的SPV方程結(jié)合有限元中的Newton-Raphson求解法和網(wǎng)格單元剛度矩陣有效及時調(diào)整法對二維條件下的土石壩水流滲透問題進行了研究，結(jié)果表明管涌現(xiàn)象顯著增大了土石壩的滲透系數(shù)和滲透速度[1-2]。李新華等深入分析影響堤壩滲流的各關(guān)鍵參數(shù)后，結(jié)合一具體工程實例中監(jiān)測的滲流數(shù)據(jù)，應(yīng)用Matla神經(jīng)網(wǎng)絡(luò)方法建立了能夠有效預(yù)測滲流參數(shù)變化的模型，為滲流參數(shù)在堤壩中的變化預(yù)測提供了借鑒[3-5]。陸忠民等對上海市一水庫建設(shè)階段存在多種潛在不安全擾動因素影響條件下的實際困難的解決進行了研究，總結(jié)了一套針對潮汐河口位置的堤壩建設(shè)方案，為類似復(fù)雜地質(zhì)條件下的堤壩建設(shè)提供了參考[6-7]。黃海均等對于常見強弱互層的堤壩基礎(chǔ)結(jié)構(gòu)形式的滲透破壞特點進行了研究，采用PFC3D軟件結(jié)合“反演法”，首先對壩基中各土層的顆粒細(xì)觀參數(shù)進行了正確標(biāo)定，得出壩基滲透變形破壞發(fā)展演化過程中各土層細(xì)顆粒的運動特征，揭示了強弱互層堤壩的滲透破壞微觀變形發(fā)展機制[8-9]。武亞軍等研究并成功制配出用于與黏性土滲透變形內(nèi)部變化過程的透明質(zhì)黏土材料，并進行了多種配比條件下的透明黏性土滲透試驗，得出黏土的滲透系數(shù)隨著黏質(zhì)成分的含量提高而變小，而黏土的壓縮性會隨著黏質(zhì)成分的含量提高而變大[10-11]。

以下對位于山西省一水庫堤壩進行數(shù)值計算分析，對比當(dāng)下游邊界設(shè)置與不設(shè)置滲溝時，堤壩內(nèi)的水力梯度、流線形式和孔壓變化特征，采用能夠有效進行飽和-非飽和分析的Geo-studio軟件建立實體堤壩模型，討論滲漏邊界設(shè)置與否相應(yīng)的各水力關(guān)鍵參數(shù)的變化趨勢，總結(jié)了滲漏邊界的出現(xiàn)與否對各關(guān)鍵水力參數(shù)的影響特征，為類似地質(zhì)條件下的堤壩設(shè)計和工程維護提供了一定參考。

1 工程地質(zhì)

該水庫堤壩位于遼寧省內(nèi)，是為該地區(qū)汛期階段調(diào)整各小區(qū)域水量分布而建設(shè)，壩基主要由土石混合體材料組成，土石混合體材料相關(guān)物理力學(xué)參數(shù)如表1所示，壩基剖面模型如圖1所示。預(yù)計汛期儲水高度為11m，堤壩基礎(chǔ)長52m，高12m，頂寬4m。

2 Geo-studio滲流分析

Geo-studio軟件在滲流分析上有著獨特的優(yōu)勢，其軟件自帶有相關(guān)的水力特征曲線估計算法，含水量與滲透系數(shù)估計算法及相應(yīng)的孔壓與含水量估計算法，輸入相應(yīng)的參數(shù)，軟件即可自行計算得出，內(nèi)置有多種前人研究成果中的各種水力相關(guān)函數(shù)模型，為各個工程領(lǐng)域涉及到流體滲流問題的分析提供了極大便利。

圖1 壩基剖面圖

表1 土石混合體物理力學(xué)參數(shù)

3 堤壩模型構(gòu)建

在Geo-studio軟件中建立如圖2所示堤壩模型，通過設(shè)置材料特性，并賦給相應(yīng)區(qū)域的材料，指定水位面位置，按照0.5m的單元網(wǎng)格密度劃分模型，模型底部右端45-52m處設(shè)置或不設(shè)置滲漏邊界，進行計算。

圖2 水庫堤壩模型

4 模擬結(jié)果

下文主要通過分析設(shè)置與不設(shè)置滲漏邊界條件時，壩基內(nèi)水力梯度變化特征、水流流線特征及零水壓頭的分布位置的特點，得出設(shè)置與不設(shè)置滲漏邊界條件時，以上這些參數(shù)的變化規(guī)律。

4.1 不設(shè)置滲漏邊界時

圖3所示是不設(shè)置滲漏邊界時，邊坡內(nèi)水力梯度等值線在壩基內(nèi)的分布位置，和幾個主要流線的形式，以及零水頭壓力線在壩基中的分布位置，可見，從壩基上游開始發(fā)生滲透，水力梯度在最初時，坡降較小，越向下游，坡降越大，水流在通過壩基進行滲透時，靠近上游處，以垂直加水平的混合運動方式為主，向下游滲流時，迅速過渡到以水平運動為主，零壓力水頭線以一定斜角趨勢幾乎成直線形式指向下游。

圖3 不設(shè)置滲漏邊界時的水力參數(shù)圖

圖4 不設(shè)置滲漏邊界時滲透率和基質(zhì)吸力關(guān)系

圖4所示為在不設(shè)置滲漏邊界時的土石混合體材料的滲透率(水力傳導(dǎo)度)和基質(zhì)吸力關(guān)系曲線。由圖中所示可知，當(dāng)滲透率隨著土石混合體材料中的基質(zhì)吸力的不斷增大而減小，當(dāng)達(dá)到1kPa的臨界基質(zhì)吸力值時，滲透率出現(xiàn)迅速陡降。

4.2 設(shè)置滲漏邊界時

圖5所示是設(shè)置滲漏邊界時，邊坡內(nèi)水力梯度等值線在壩基內(nèi)的分布位置，和幾個主要流線的形式，以及零水頭壓力線在壩基中的分布位置，可見，從壩基上游開始發(fā)生滲透，水力梯度在最初時，坡降較小，越向下游，坡降越大，與不設(shè)置滲漏邊界情況相比，整體坡降較大。水流在通過壩基進行滲透時，靠近上游處，以垂直加水平的混合運動方式為主，向下游滲流時，迅速過度到以水平運動為主，這一點與不設(shè)置滲漏邊界情況相似，零壓力水頭線不再穿過下游斜坡，先以一定斜角趨勢也幾乎成直線形式指向下游，后在滲漏邊界處迅速回轉(zhuǎn)向下至滲漏邊界，說明滲漏邊界的設(shè)置能夠有效預(yù)防壩基發(fā)生整體的失穩(wěn)滲漏破壞。

圖5 設(shè)置滲漏邊界時的水力參數(shù)圖

圖6所示為在設(shè)置滲漏邊界時的土石混合體材料的滲透率(水力傳導(dǎo)度)和基質(zhì)吸力關(guān)系曲線。

由圖6可知，當(dāng)滲透率隨著土石混合體材料中的基質(zhì)吸力的不斷增大而減小，當(dāng)達(dá)到1kPa的臨界基質(zhì)吸力值時，滲透率出現(xiàn)迅速陡降，這與不設(shè)置滲漏邊界時的情況相同，可見滲漏邊界的設(shè)置對滲透率和基質(zhì)吸力的關(guān)系沒有影響。

5 結(jié) 論

1)水力梯度在靠近上游處，坡降較小，越向下游，坡降越大，水流在通過壩基進行滲透時，越靠近上游處，以垂直加水平方向向下游運動為主，而后迅速過渡到以水平運動為主。

2)與不設(shè)置滲漏邊界情況相比，設(shè)置滲漏邊界時，零壓力水頭線不再穿過下游斜坡，說明滲漏邊界的設(shè)置能夠有效預(yù)防壩基發(fā)生整體的失穩(wěn)滲漏破壞。

3)建議在相似水位和壩基結(jié)構(gòu)模式下進行堤壩的設(shè)計時，采用下游適當(dāng)設(shè)置滲漏邊界的措施。