水位升降條件下馱英水庫上游圍堰滲流及穩定性分析

勞乃勝，王金龍，周欣華

（1.崇左市左江治旱工程管理中心，廣西 崇左 532200；2.長江水利委員會長江科學院 水利部巖土力學與工程重點實驗室，武漢 430010；3.國家大壩安全工程技術研究中心，武漢 430010）

0 引言

圍堰作為河道中修建的臨時擋水建筑物，主要用來維護永久水工建筑物的干地施工，其自身安全直接關系到主體建筑物的施工安全、工期及造價等。因此，圍堰自身穩定性對于整個工程極其重要［1］。

滲流問題與土石圍堰的安全運行有著緊密的聯系，很多土石圍堰的失事都是水體的滲透造成的。壩體出現滲漏問題后會出現不利的滲流場分布，從而導致土石圍堰發生滲透變形，更嚴重的會導致圍堰堰體邊坡滑移或垮塌[2]。此外，水位驟升或驟降都將使土石圍堰堰體內滲流場在短時間內發生較大的變化，進而影響堰體尤其是上游迎水坡的穩定性。長期蓄水后水位驟降，土石圍堰上游迎水坡由于孔隙水壓力來不及消散形成向邊坡滲流，當水庫內的水位降低后土壩內的浸潤面依然處于較高位置，不但造成土壩內土料的容重隨著水位變化而改變，還產生非穩定滲流現象和滲透壓力，使得土體顆粒間的有效應力快速變小，降低壩體土料的抗剪強度，從而影響到壩體的安全穩定性[3]。庫水位驟升對上游迎水坡的穩定也是不利的。因此，有必要對庫水位驟升驟降工況下土石圍堰進行滲流及穩定性分析，為圍堰安全評價提供參考[4]。

1 工程概況

馱英水庫主要建筑物有瀝青混凝土心墻堆石攔河壩、右岸溢洪道、右岸河道電站發電引水隧洞、右岸泄洪隧洞、左岸灌溉發電引水系統、渠首電站廠房、壩后河道電站廠房及場內永久交通道路等。

攔河壩采用瀝青砼心墻堆石壩壩型，最大壩高為72.2 m，壩頂總長為225 m，壩頂寬8 m，壩頂高程為233.2 m，壩頂防浪墻頂高程234.4 m，壩基面最低高程為161 m。壩體材料分區：從上游往下游依次分為上游干砌石護坡、主堆石IA區、主堆石IVA區、次堆石IIA 區、過渡層IIIA 區；瀝青混凝土心墻；過渡層IIIB 區、主堆石IB 區、次堆石表面保護層堆石IVB區、次堆石IIB區及下游干砌石護坡等。

攔河壩上游圍堰填筑施工斷面見圖1。堰體上游采用土工膜防滲，堰體填料采用強風化料，上游設計坡比為1∶2.5，下游設計綜合坡比為1∶1.8。在高程178～183 m 范圍內填料的粘粒含量較高，透水性相對較弱。

圖1 上游圍堰計算斷面圖

依據導流洞泄洪全年10年一遇庫水位過程，庫前水位降落條件為：0～第55 小時為水位上升過程，水位由176.0 m 升高至205.83 m，第55～第100 小時為水位降落過程，水位由205.83 m降落至182.97 m。本文針對施工期圍堰擋水度汛工況開展滲流與穩定性分析。

2 基本理論

2.1 滲流基本理論

根據達西定律和質量守恒連續性原理，以壓力水頭為未知量，二維飽和-非飽和滲流的基本方程為：

式中：kr為相對滲透率；kij為飽和滲透張量；h為壓力水頭；x2為正向向上的鉛直坐標；C為比容水度；β為系數，非飽和區β=0，飽和區β=1；Ss為單位貯存量；t為時間；S為源匯項。

定解條件包括初始條件和邊界條件。

2.2 穩定性分析理論

采用簡化Bishop法對壩坡進行穩定性分析，由圓弧滑面確定的滑體被分成n 個垂直條塊。第i 條塊寬度bi，底面傾角αi，條塊體重力Wi，水平條間作用力Ei和Ei+1，垂直條間作用力（條間剪力）Xi和Xi+1，條底法向作用力Ni，條底剪力Ti，水壓力的合力Ui，滑面內擦角φi，黏聚力ci。滑面安全系數為Fs。A.W.Bishop 假設整個滑面抗剪強度按同一比例調用（亦稱Bishop安全系數）。

穩定安全系數的計算公式：

（1）檢驗結果比較。47例肺部孤立性肺結節患者均接受CT增強掃描與常規病理檢查，比較2種檢查方式的檢驗結果。

依據非飽和土滲流理論和極限平衡理論，本文運用有限元分析軟件GeoStudio 對庫水位升降條件下圍堰滲流場及穩定性進行模擬分析。

3 滲流及穩定性分析

3.1 模型建立

圍堰汛期非穩定滲流主要是分析上游水位上升和降落過程中，堰體內水位上升和消落過程，分析的重點區域是堰體，由于基礎開挖至巖土分界線，因此堰基部分的滲流場分布對于降落期堰體的水頭分布影響較小，計算模型只考慮堰體部分。

堰體從上游至下游分為4個填筑分區：護坡、反濾墊層、過渡料（含防滲土工膜）、堆石料。堆石料依據填料檢測和試驗成果，又分為下部、中部、上部3個分區。

3.2 參數條件

堆石料建基面至178 m 高程為下部，飽和滲透系數為1.0×10-5cm/s；178～183 m 高程為中部，飽和滲透系數為1.15×10-6cm/s；183 m至堰頂為上部，飽和滲透系數為1.6×10-3cm/s。考慮防滲土工膜失效的不利工況，土工膜所在的過渡料滲透系數取值為1.0×10-5cm/s。圍堰滲流及穩定性計算中材料介質的基本物理力學參數見表1。

填料名稱護坡墊層過渡料堆石料上部（183 m至堰頂）堆石料中部（175～183 m高程）堆石料下部（175 m以下）干密度/（g/cm3）2.06 2.06 2.06抗剪強度指標c/（kPa）15.3 15.3 15.3 φ/（°）26.7 26.7 26.7滲透系數/（cm/s）1.0×10-3 1.0×10-3 1.0×10-5 2.11 26.8 27.3 1.6×10-3 1.90 15.3 26.7 1.15×10-6 2.11 26.8 27.3 1.0×10-5

3.3 計算工況

本次計算的運行工況主要為施工渡汛期，考慮圍堰擋10年一遇洪水，同時計算洪水過后上游水位的退水工況，計算不同退水時間的非穩定滲流浸潤線時上、下游堰體壩坡的抗滑安全系數，為了了解退水過程上游邊坡局部抗滑穩定性，增加退水過程定折線滑弧局部穩定性計算工況。具體計算工況見表2。

方案編號F1 F2 F3條件及說明上游水位176 m，下游水位164，0 m，穩定滲流；結果作為非穩定滲流的初始條件。上游水位205.83 m，下游水位164，0 m，穩定滲流；結果可作為水位升降過程非穩定滲流的對照，并可作為堰體穩定計算的不利水位條件。上游水位對應導流洞泄洪全年10年一遇庫水位過程，上游水位先升后降，歷時100 h。

3.4 滲流計算分析

穩定滲流場水頭等值線分布，非穩定滲流堰體浸潤線升降過程見圖2。

方案F1為穩定滲流，上游水位為176.0 m，下游水位為164.9 m，堰體自由面較平緩，下游坡面出逸點高程為165.59 m，堰腳處出逸比降約為0.55。

圖2 滲流計算滲流場成果圖

方案F2 為穩定滲流，上游水位為205.83 m，下游水位為164.9 m，此工況堰體水位最高，可作為堰體穩定分析不利的水位條件，下游坡面出逸點高程為186.58 m，出逸比降為0.415，位于堆石料上部。

3.5 穩定性分析

根據上述計算斷面、計算參數、水位等計算條件和計算方法，對圍堰上、下游邊坡施工度汛期的運行工況進行了穩定性計算分析，各工況的抗滑安全系數見表3。最小抗滑安全系數對應的滑弧位置見圖3、圖4。

在不考慮圍堰土工布防滲條件下的堰體邊坡穩定性計算成果表明：

（1）圍堰汛期正常擋水形成穩定滲流狀態時（水位205.83 m），上游邊坡的抗滑安全系數為2.65，滿足規范要求，而下游邊坡的抗滑安全系數只有0.97，不滿足規范要求，主要是由于上游坡未考慮土工膜防滲，坡體浸潤線位置與出逸點均較高，屬于極端工況。

（2）漲水過程中，考慮較短時間漲至汛期正常擋水水位205.83 m 時，隨著水位上升，上游邊坡的抗滑安全系數增大，下游邊坡的抗滑安全系數基本維持不變，為1.44，上、下游邊坡的抗滑安全系數均滿足規范要求，主要是由于漲水時間較短，坡體未形成穩定滲流狀態，浸潤線位置相對較低。

（3）退水期考慮圍堰上游短時間內漲至汛期正常水位205.83 m后即時退水，退水期坡體滲流場處于非穩定滲流狀態，坡體浸潤線與漲水過程的浸潤線位置基本一致，隨著水位的下降，上游邊坡的抗滑安全系數減小，而下游坡的抗滑安全系數基本維持不變，整個退水過程圍堰上、下游邊坡的抗滑安全系數，均滿足規范要求。圍堰上游坡綜合坡比較小，定折線局部滑弧的抗滑安全系數均大于對應整體滑弧安全系數，說明計算退水過程中出現淺層脫坡的可能性不大。但仍需注意在有土工膜的實際條件下，由于土工膜阻止坡內水的出溢，而在退水過程可能發生的局部破壞。

工況編號F2備注穩定滲流期堰前水位及浸潤線205.83 m水位穩定滲流期上升20小時水位182.81 m上升40小時水位183.87 m上升55小時水位205.83 m退水20小時水位190.13 m穩定標準1.25抗滑安全系數漲水期1.25退水期F3 1.44滑面特性整體滑動整體滑動整體滑動整體滑動整體滑動局部滑動整體滑動局部滑動整體滑動局部滑動退水40小時水位183.87 m 1.25下游坡0.97 1.44 1.44 1.44 1.44非穩定滲流期退水45小時水位183.0 m上游坡2.65 1.93 2.10 4.35 2.25 3.19 1.95 2.74 1.91 2.74 1.44

圖3 方案F2汛期高水位穩定滲流期典型滑弧圖

圖4 方案F3汛期庫水升降過程圍堰上游邊坡典型滑弧圖

4 結論

（1）汛期高水位穩定滲流工況堰體水位最高，堆石料出逸比降為0.415，上游邊坡的抗滑安全系數為2.65，滿足規范要求；而下游邊坡的抗滑安全系數只有0.97，不滿足規范要求，該安全系數值反映上游坡考慮土工膜防滲失效且長期未有效退水的狀態，坡體浸潤線位置與出逸點均較高，屬于極端工況。汛期在該條件下建議考慮對圍堰下游坡進行反濾保護和壓坡處理。

（2）施工期度汛水位先升后降，非穩定滲流計算結果表明：由于上游側過渡區滲透性較低，庫水位上升時間較短，坡體未形成穩定滲流狀態，浸潤線位置相對較低。上游水位變動，只對過渡料上游側區域影響較大，對下游側的堰體影響很弱，歷時100 h 先升后降的非穩定滲流過程中，堰體水位基本與初始時刻相同。

（3）漲水過程中，隨著水位上升，上游邊坡的抗滑安全系數增大，下游邊坡的抗滑安全系數基本維持不變。退水過程中，隨著水位的下降，上游邊坡的抗滑安全系數減小，而下游坡的抗滑安全系數基本維持不變，整個退水過程圍堰上、下游邊坡的抗滑安全系數，均滿足規范要求。

（4）圍堰上游坡綜合坡比較小，定折線局部滑弧的抗滑安全系數均大于對應整體滑弧安全系數，說明計算退水過程中出現淺層脫坡的可能性不大。但仍需注意在有土工膜的實際條件下，由于土工膜阻止坡內水的出溢，而在退水過程可能發生的局部破壞。