水位變動速度對堆積體穩(wěn)定性的影響

水位變動速度對堆積體穩(wěn)定性的影響

在國庫水位上升或下降的過程中，水庫兩側(cè)斜坡上（不同河流體積的差異）地下水滲流場及其物理力學(xué)性質(zhì)的變化是一個動態(tài)適應(yīng)過程。在這個過程中，如何預(yù)測地下水滲流場和邊坡巖石物理力學(xué)性質(zhì)的動態(tài)變化，以及對邊坡穩(wěn)定性的影響一直是工程和學(xué)術(shù)界面臨的一個重要問題。針對庫水位升降對庫區(qū)斜坡穩(wěn)定性的影響,近年來國內(nèi)外眾多學(xué)者展開了一些探索和研究,取得了一定的成果。如張文杰、柳群義、馬崇武等從非飽和土理論角度探討了水位升降對斜坡滲流場及其穩(wěn)定性的影響;李邵軍等對水位升降條件下斜坡失穩(wěn)模式展開了離心試驗研究;曹玲等以千將坪滑坡為例,研究了水位升降過程中干-濕循環(huán)條件下滑坡滑帶土強(qiáng)度特性及其變形特性;徐文杰等從流-固耦合的角度采用數(shù)值計算方法對庫水位升降條件下斜坡穩(wěn)定性狀況進(jìn)行了分些和探討;何曉英等對周期性浸泡條件下松散堆積體的微觀物性特征進(jìn)行了試驗研究。綜合目前研究成果,國內(nèi)外學(xué)者認(rèn)為庫水位升降對斜坡穩(wěn)定性的影響主要表現(xiàn)在以下幾個方面:1)水位升降對斜坡滲流場的影響;2)水位升降對斜坡物理力學(xué)性質(zhì)的影響;3)水位升降對斜坡失穩(wěn)破壞模式的影響。既然庫水位升降是一個動態(tài)過程,其對斜坡滲流場、物理力學(xué)性質(zhì)等方面的影響是一個動態(tài)過程,那么不同的水位升降速度對斜坡滲流場、物理力學(xué)性質(zhì)的影響程度也應(yīng)該存在一定差別,目前關(guān)于水位升降速度對斜坡穩(wěn)定性影響的研究成果鮮見報道。中國許多水庫都需要承擔(dān)發(fā)電、防洪以及調(diào)解徑流量的綜合功效,因此如何合理地控制庫水位升降速度,在保證發(fā)揮其最大社會經(jīng)濟(jì)效益的同時,避免庫區(qū)斜坡出現(xiàn)影響沿岸生產(chǎn)、生活活動的大規(guī)模失穩(wěn)破壞事件成為一項重要課題。作者以瀾滄江某電站庫區(qū)巨型堆積體為例,探討庫水位升降速度對其地下水滲流場及其穩(wěn)定性的影響,以期拋磚引玉,為相關(guān)研究提供參考。瀾滄江某電站庫區(qū)發(fā)育一巨型冰水堆積體,堆積體距大壩500m的常年回水區(qū)(或變動回水區(qū)),沿河方向?qū)捈s800m,垂直河流方向長約460m,總方量約6.6×106m3(如圖1所示)。電站設(shè)計死水位為2150m,正常蓄水位為2250m,堆積體前緣大部位于庫水位變動范圍內(nèi)。由于堆積體所在河段在水位變動過程中沿庫縱向水流速度較小,因此忽略沿水庫縱向的水流速度對堆積體穩(wěn)定性的影響。堆積體物質(zhì)構(gòu)成主要是形成于第四系上中更新統(tǒng)的冰水堆積物,下伏基巖為二疊系下統(tǒng)吉東龍組砂板巖。1堆積體物理力學(xué)性質(zhì)應(yīng)是由水位動一般在庫水位上升或下降過程中,對庫區(qū)上述堆積體穩(wěn)定性所構(gòu)成的影響主要表現(xiàn)在兩個方面,一是庫水位變動所導(dǎo)致的堆積體滲流場的調(diào)整和變化,二是水位變動過程中,伴隨庫水浸入、地下水浸潤線抬升或地下水排出、地下水浸潤線下降引起的堆積體物理力學(xué)性質(zhì)的變化。基于上述認(rèn)識,探討庫水位變動速度對堆積體穩(wěn)定性影響的思路可分述為以下兩個部分。1.1基于gler東南角加權(quán)余量法的二維滲流方程庫水位變動期間,地下水滲流場的調(diào)整是一個非常復(fù)雜的過程,國內(nèi)外眾多學(xué)者在這方面進(jìn)行了廣泛研究,然而在確定地下水浸潤線時大多是建立在半經(jīng)驗、半理論基礎(chǔ)上,因此不同的人和單位提出的計算方法和公式往往缺乏普適性,并且計算得出的結(jié)果大多是庫水位變動到某一水位穩(wěn)定下來后的地下水穩(wěn)定滲流浸潤線。有限元方法的引入為解決這一問題提供了一條可行途徑,尤其是在求解庫水位變動期間,地下水瞬時滲流場時具有獨(dú)到的優(yōu)勢。Fredlund等提出對于各向同性巖土體非穩(wěn)定滲流場的控制性微分方程可描述為:??x(k?H?x)+??y(k?H?y)+Q=mwγw?H?t(1)??x(k?Η?x)+??y(k?Η?y)+Q=mwγw?Η?t(1)式中,k為滲透系數(shù);H為水頭,隨時間而變化;Q為邊界滲流量;mw為水土特征曲線的斜率;γw為水容重。根據(jù)上述微分方程,采用Glerkin加權(quán)余量法可推導(dǎo)出適合有限元計算的2維滲流方程:τ∫A([B]T[C][B])dA{H}+τ∫A(λ〈B〉T〈N〉)dA[H],t=qτ∫L〈N〉TdL。(2)式中,[B]為水力梯度矩陣,[C]為單元導(dǎo)水率矩陣,{H}為水頭向量,〈N〉為插值函數(shù)向量,q為單元邊界單位滲流量,τ為單元厚度,t為時間,λ值為mwγw,A為單元面積,L為單元邊長。根據(jù)上述理論,采用SEEP/W有限元計算軟件,計算出不同水位變動速度條件下,不同庫水位、不同水位變動時刻的堆積體內(nèi)瞬時地下水浸潤線以及堆積體內(nèi)地下水滲流場的變化特征。1.2土體浸漬線含水層穩(wěn)定性計算方法水位變動期間堆積體穩(wěn)定性計算模型可概括為如圖2所示,穩(wěn)定性計算可分為兩部分,對于庫水位以上坡體考慮滲流壓力,而對于庫水位以下坡體,將坡面以上水體看成一種有重度無強(qiáng)度的特殊材料,根據(jù)該思路,按簡化Bishop條分法,堆積體穩(wěn)定性系數(shù)可按下式計算:ks={∑i=1n1mai[cilicosαi+(W1i+W2i′+Disinβi)tanφi]+∑j=n+1m1mai[cjljcosαj+(W1j+W2j?Uj)tanφj]}/{∑i=1n[(W1i+W2i′)sinαi+Dicos(αi?βi)]+ks={∑i=1n1mai[cilicosαi+(W1i+W2i′+Disinβi)tanφi]+∑j=n+1m1mai[cjljcosαj+(W1j+W2j-Uj)tanφj]}/{∑i=1n[(W1i+W2i′)sinαi+Dicos(αi-βi)]+∑j=n+1m(W1j+W2j)sinαj}∑j=n+1m(W1j+W2j)sinαj},ma=cosα+tanφkssinα。ma=cosα+tanφkssinα。式中,c為土體內(nèi)聚力;φ為內(nèi)摩擦角;W1i為條塊浸潤線以上土體重量;W2i′為條塊浸潤線以下土體浮重量;W1j為條塊上覆庫水重量;W2j為條塊飽和重量;α為滑面傾角;β為地下水浸潤線傾角;l為條塊滑面長度;n為庫水位以上坡體條塊數(shù)量;m坡體條塊總數(shù);D為滲流壓力,D=γwhwlcosαsinβ。按照該計算方法,將根據(jù)1.1節(jié)方法得到的地下水瞬時浸潤線導(dǎo)入穩(wěn)定性計算模型中,進(jìn)行穩(wěn)定性計算,即可獲得在不同水位變動速度條件下,不同庫水位、不同水位變動時刻堆積體的穩(wěn)定性系數(shù)。通過對地下水瞬時滲流場和穩(wěn)定性計算成果的分析,得出水位變動速度對堆積體穩(wěn)定性影響的基本認(rèn)識。水位變動過程中,隨著堆積體內(nèi)地下水的浸入和排出,其物理力學(xué)性質(zhì)的變化將是一個復(fù)雜的動態(tài)過程,目前關(guān)于這方面的研究成果仍不夠成熟,因此沒有討論水位變動期間,堆積體在飽和與非飽和變化過程中巖土體強(qiáng)度的變化及其對堆積體穩(wěn)定性的影響,在分析過程中,為方便起見,假設(shè)地下水瞬時浸潤線以上巖土體的物理力學(xué)參數(shù)取天然值,瞬時浸潤線和庫水位以下的巖土體取有效值。2不同水位變動速度下根據(jù)相關(guān)試驗資料,計算過程中巖土體的滲流參數(shù)見表1。計算過程中通過改變庫水位的上升或下降速度,從而計算得到不同水位變動速度條件下,不同庫水位、不同水位變動時刻的堆積體內(nèi)瞬時地下水浸潤線,計算過程中設(shè)定庫水位變動速度分別為:0.5、1、2、3和5m/d。2.1水位上升過程中的水力梯度圖3~5是堆積體在不同水位上升速度條件下地下水瞬時浸潤線的變化、調(diào)整過程和規(guī)律,以及水位變動過程中堆積體瞬時浸潤線平均水力梯度與時間、水位之間的關(guān)系。根據(jù)計算結(jié)果,可以看出水位上升過程中水位上升速度對堆積體滲流場的影響主要表現(xiàn)在以下幾個方面:1)水位上升過程中,受地下水滲流影響,在達(dá)到正常蓄水位之前,堆積體內(nèi)瞬時浸潤線的抬升明顯滯后,瞬時浸潤線呈明顯上凹型,并且隨著水位上升速度的增大這種特點(diǎn)愈發(fā)強(qiáng)烈。2)從圖4、5可以看出(圖中水力梯度負(fù)值表示地下水流動方向指向堆積體坡內(nèi)),隨著庫水位上升速度的增大,在達(dá)到正常蓄水位2250m以前,堆積體瞬時浸潤線的反向滲透平均水力梯度隨時間增大和庫水位的升高呈逐漸增大趨勢,在達(dá)到正常蓄水位滯后,平均水力梯度隨著時間的增大逐漸降低。另外圖6還表明,2210m水位似乎是一個臨界水位,在達(dá)到此水位之前,水力梯度隨水位增長的速度較平緩,但超過該水位后,水力梯度隨水位進(jìn)入陡增階段,并且水位上述速度越大,趨勢越明顯。處于2210m水位時,坡體浸水深度約位于坡高的1/3。2.2地下水在死水位之前的水力梯度變化從圖6~8所示水位下降過程中,隨著水位下降速度的變化,堆積地下水滲流特點(diǎn)可發(fā)現(xiàn)以下規(guī)律:1)水位下降過程中,受地下水滲流的相對滯后性影響,在達(dá)到設(shè)計死水位前,地下水瞬時浸潤線斜率(即水力梯度)一般呈逐漸增大特點(diǎn),一旦達(dá)到設(shè)計死水位,隨著時間增長逐漸減小,直至形成穩(wěn)定滲流場。2)從圖8所示的水位下降速度與瞬時浸潤線平均水力梯度曲線來看,并不是水位下降速度越大,坡體瞬時浸潤線水力梯度越大,而是存在一個臨界水位下降速度值,使得在不同水位條件下,瞬時浸潤線的水力梯度獲得最大值,從該堆積體來看,該水位下降速度臨界值為2m/d,超過該臨界值后,隨著水位下降速度的增大,平均水力梯度反而減小。3地下水水土流失對斜坡穩(wěn)定性的影響水位變動期間,影響堆積體穩(wěn)定性的因素多且復(fù)雜,如水位上升期間有庫水位以下坡體孔隙水壓力增大、水位以上地下水反向滲流、地下水位線抬升導(dǎo)致的巖土體強(qiáng)度降低等等;水位下降期間,有滲透壓力增大、地下水排出巖土體強(qiáng)度恢復(fù)、庫水位以下坡體孔隙水壓力降低等等。在這些因素中,有些對斜坡穩(wěn)定有利,而有些則是不利的,水位變動期間斜坡穩(wěn)定性狀況就是這些因素綜合作用的結(jié)果。根據(jù)對堆積體坡體結(jié)構(gòu)特征的分析,結(jié)合穩(wěn)定性計算結(jié)果表明,水位變動期間堆積體可能的失穩(wěn)模式有兩種,一種是沿基覆界限整體滑移,一種是堆積體前緣出現(xiàn)局部失穩(wěn),據(jù)相關(guān)試驗成果,堆積體物理力學(xué)參數(shù)如表2所示。3.1水位上升過程中穩(wěn)定性系數(shù)的變化從圖9、10可見,水位上升期間,當(dāng)水位上升速度較小時,堆積體局部和整體穩(wěn)定性系數(shù)都表現(xiàn)出先降低后升高,至最高水位達(dá)到最大值,而后又隨時間逐漸減小,而當(dāng)水位上升速度較大時,穩(wěn)定性系數(shù)一般表現(xiàn)為先升高,至最高水位達(dá)到最大值,而后隨時間減小,對于該堆積體而言,該臨界水位上升速度約3m/d。堆積體穩(wěn)定性系數(shù)一般在最高水位形成穩(wěn)定性滲流場后達(dá)到最小值,但當(dāng)水位上升速度非常小時(<0.5m/d)堆積體前緣局部穩(wěn)定性系數(shù)在蓄水至2210m左右取得最小值。圖11、12表明(計算結(jié)果為水位上升至不同水位時穩(wěn)定性系數(shù),不包括形成穩(wěn)定性滲流期間穩(wěn)定性系數(shù)的變化),在水位以較低速度(<3m/d)上升過程中,堆積體穩(wěn)定性系數(shù)取得最小值的水位隨著水位上升速度減小逐漸增大。上述分析結(jié)果說明,水位上升速度較大對堆積體穩(wěn)定性有利,過小的水位上述速度(<0.5m/d)更易導(dǎo)致堆積體前緣失穩(wěn)。水位上升期間,當(dāng)速度較大時,應(yīng)注意蓄水初期的穩(wěn)定性狀況,而速度較小時,堆積體在較高水位的穩(wěn)定性狀況往往較差。3.2堆積體穩(wěn)定性與水位的關(guān)系從水位下降至形成穩(wěn)定性滲流場期間堆積體穩(wěn)定性系數(shù)的變化來看(圖13、14),隨著水位下降,堆積體穩(wěn)定性系數(shù)都表現(xiàn)出降低趨勢,且隨著下降速度的增大,穩(wěn)定性系數(shù)下降越強(qiáng)烈,達(dá)到最低值后,隨著時間增長,堆積體穩(wěn)定性系數(shù)逐漸回升。圖15、16所示的堆積體穩(wěn)定性與水位關(guān)系曲線表明(計算結(jié)果為水位下降至不同水位時穩(wěn)定性系數(shù),不包括形成穩(wěn)定性滲流期間穩(wěn)定性系數(shù)的變化),水位下降過程中,普遍存在一中間水位取得堆積體穩(wěn)定性系數(shù)的最小值,超過該水位后,隨著水位進(jìn)一步下降,堆積體穩(wěn)定性系數(shù)逐漸回升,僅當(dāng)水位下降速度很小時(<0.5m/d),上述中間危險水位才會消失。并且隨著水位下降速度增大,堆積體穩(wěn)定性系數(shù)取得最小值的水位逐漸降低。總體而言,水位下降速度增大對堆積體穩(wěn)定性不利,水位下降速度由0.5m/d增大至5m/d,堆積體穩(wěn)定性系數(shù)可降低15%~20%,堆積體穩(wěn)定性系數(shù)最小值往往出現(xiàn)在水位下降至堆積體下部1/3~1/5坡高處。4水位下降過程中堆積體穩(wěn)定性最小值采用有限元數(shù)值計算和剛體極限平衡分析相結(jié)合的方法,對瀾滄江某堆積體在水位變動期間的水位變動速度對穩(wěn)定性的影響進(jìn)行了初步探討,結(jié)果表明較大的水位上升速度對其穩(wěn)定性有利,當(dāng)速度較小時,水位上升過程中堆積存在危險水位,若水位上升速度過小(<0.5m/d),堆積體穩(wěn)定性最小值并不是出現(xiàn)在最高水位處,而是出現(xiàn)在中間過程的2210m水位。較小的水位下