土壩滲流及穩定分析在水庫除險加固工程中的應用

溫 州，楊鴻剛，趙玉新

(黑龍江省水利水電勘測設計研究院，哈爾濱150080)

1 土壩的滲流及滲透穩定計算方法

目前，對于土壩的滲流計算方法分為有限元法和公式法，根據《小型水利水電工程碾壓式土石壩設計導則》(SL189—96)，針對小型水利工程可采用簡單實用的公式法進行計算。對于透水地基滲流計算采用公式如下:

式中:q 為總單位寬度滲流量，m3/s·m;qD為不透水堤基上求得相同排水的單位寬度滲流量，m3/s·m;k 為壩身滲透系數，m/s;k0為壩基滲透系數，m/s;H1為上游水位，m;H2為下游水位，m;h0為下游逸出點高程，m;m1為上游坡坡率;m2為下游坡坡率;T 為透水地基的深度。

浸潤線方程為:

通過計算，可求得特征水深h0。

不透水地基均質土堤坡面滲流比降計算:

滲出點:

堤坡與不透水面交點:

壩坡的抗沖刷能力計算:

壩坡抗沖刷破壞的臨界比降計算:

式中:γ′為土的浮容重;γw為水的容重;γ 為土的濕容重;φ 為土的水下內摩擦角;β 為壩坡的坡角;c 為土的凝聚力;K 為安全系數。

2 土壩的邊坡穩定計算方法

目前廣泛被采用的邊坡穩定計算方法有瑞典條分法、簡化的Bishop 法及JanBu 法等。針對小型土壩工程，可采用考慮條塊間作用力的簡化Bishop 法進行壩坡穩定分析。計算公式如下:

式中:W 為土條重量;Q、V 為分別為水平和垂直地震慣性力(向上為負，向下為正);μ 為作用于土條底面的孔隙壓力;α 為條塊重力線與通過此條塊底面中點的半徑之間的夾角;b 為土條寬度;c′、φ′為土條底面的有效應力抗剪強度指標;Mc 為水平地震慣性力對圓心的力矩;R 為圓弧半徑。

3 工程實例

文章針對某均質土壩的除險加固工程，采用上述滲流及邊坡穩定分析方法對土壩進行計算。該土壩壩頂寬度約6m，最大壩高約12m，壩頂長度290m，上游邊坡1:1.8 ～1:3.5，下游邊坡1:2 ～1:3。壩體上下游無護破，壩體斷面極不規整，壩軸線曲折、現狀壩坡不規則，上下游壩坡受水流沖蝕破壞嚴重、壩后無排水設施。其滲流計算參數見表1:

表1 滲流計算參數表 kN·m－3

滲流計算下游出逸點高度h0按透水地基考慮，據式(1)、(2)求出;下游壩坡滲透比降按不透水地基考慮，據式(3)、(4)求出，再由式(5)、(6)求得下游坡的允許滲透比降。滲流計算結果見表2。

表2 壩后坡滲流計算成果表

通過計算可以看出，壩后坡滲流出逸點較低，壩后坡出逸比降大于允許比降，故本次設計在下游坡坡腳處加設一層無紡布作為反濾層，防止在坡腳處出現流土等現象。土壩壩坡穩定計算應用理正軟件程序，采用簡化的Bishop 法，在滲流分析的基礎上將孔隙水壓力場引入到邊坡穩定計算中，由試驗得出有效抗剪強度指標c′、φ，采用有效應力法進行穩定滲流期的邊坡穩定分析。其物理力學指標見表3:

表3 計算采用的土質物理力學指標表

采用的計算代表斷面為0+250，現狀斷面上游坡比 為1∶1.5 ～1∶1.2.25，下 游 坡 比 為1∶1.67 ～1∶2.6，分別對現狀各種工況進行穩定計算，計算成果見表4。根據《小型水利水電工程碾壓式土石壩設計導則》(SL189—96)，采用本次設計的上游坡為1:2.5、下游坡為1∶2.0時，計算結果均滿足抗滑穩定最小安全系數正常運用條件下K ≥1.25;而現狀壩坡均不滿足規范要求，因此，本次的設計壩坡是可行的。當上游壩坡陡于1∶2.5時，放緩至1∶2.5;當下游壩坡陡于1∶2.0時，放緩至1∶2.0。

表4 壩坡穩定計算成果表

4 結 論

文章采用滲流計算中的公式法對實際工程中的土壩進行計算分析，得出下游坡腳處滲透比降＞允許滲透比降，需要在坡腳處進行反濾設計。采用簡化Bishop 法并考慮孔隙水壓力的影響，應用有效應力法對壩坡在各種工況下進行穩定分析，得出現有壩坡不滿足抗滑穩定要求，而本次設計采用的壩坡符合規范對抗滑穩定最小安全系數的要求。為工程的實施提供了理論依據。

［1］陳剛.水庫土壩滲流穩定分析及除險加固措施［J］.黑龍江水利科技，2014，42(02):67－69.