加寬堤防滲流與變形特征數(shù)值分析

郝 偉，樊恒輝，黨進(jìn)謙

(西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，陜西 楊凌 712100)

0 引 言

為發(fā)展現(xiàn)代化水利，保障防洪安全，近年來(lái)我國(guó)在大江大河段建設(shè)了多條高標(biāo)準(zhǔn)堤防道路，出現(xiàn)了許多對(duì)原有河堤進(jìn)行拓寬改造的情況[1]。普通道路是在路基的單側(cè)或雙側(cè)拓寬，但由于堤防的單側(cè)臨水特殊性，堤防拓寬只能在背水側(cè)進(jìn)行。根據(jù) “前堵后排”的設(shè)計(jì)原則，堤防加寬填料多采用滲透系數(shù)大、排水性好的砂性土，原有堤防以黏性土為填料較為常見(jiàn)。由于舊堤和新填堤身填料的不同，在季節(jié)性河水位的影響下，加寬后堤身滲流場(chǎng)將會(huì)出現(xiàn)變化，新、舊堤身之間水分的相互補(bǔ)給改變了路堤內(nèi)孔壓場(chǎng)，這將影響到堤防的穩(wěn)定性。

目前研究多集中在公路拓寬路基沉降變形上[2-5]。對(duì)于堤防的滲流和變形研究多以普通未加寬堤防為主，且多考慮堤防邊坡穩(wěn)定性和滲流穩(wěn)定性[6-9]。上述研究模型試驗(yàn)未能考慮到加寬堤防新老填料的不同特性，沒(méi)有對(duì)比舊堤與拓寬后堤防的滲流變化特征，未能耦合滲流-應(yīng)力特征考慮在水位影響下的新填堤防變形特性。本文運(yùn)用飽和-非飽和滲流有限元法，利用GeoStudio中的Seep滲流分析軟件和Sigma巖土應(yīng)力變形分析軟件對(duì)典型沿河加寬堤段滲流和變形進(jìn)行模擬和耦合計(jì)算。

1 分析理論和數(shù)值模型

1.1 分析理論

(1)滲流控制方程。Seep中的二維滲流的控制方程為：

(1)

(2)應(yīng)力計(jì)算模型。Sigma中在給定的一個(gè)時(shí)間增量上的有限元方程為：

式中：[B]為應(yīng)變-位移矩陣；[C]為本構(gòu)矩陣；{a}為節(jié)點(diǎn)x和y方向位移增量的列矢量；〈N〉為插值函數(shù)行矢量；A為沿單元邊界的面積；V為單元體積；b為體積力密度；p為面荷載增量；{Fn}為集中節(jié)點(diǎn)力增量。

1.2 數(shù)值計(jì)算模型

原堤防頂寬8 m，新填筑堤防寬12 m，堤高5 m，堤坡坡率1∶3，地基厚10 m，地基較路基左右各延伸10 m，如圖1所示。

圖1 加寬堤防數(shù)值模型Fig.1 Numerical model of widened levee

1.3 計(jì)算參數(shù)及邊界條件

文獻(xiàn)[10]通過(guò)測(cè)試不同土樣，認(rèn)為路基土所處自然條件下，土-水特性曲線(xiàn)滿(mǎn)足全范圍預(yù)估路基濕度的需要，F(xiàn)redlund & xing模型更符合實(shí)際情況，因此本試驗(yàn)采用Fredlund & xing作為計(jì)算模型：

式中：θw為體積含水量；θs為飽和體積含水量；hs為基質(zhì)吸力；a，n，m為擬合參數(shù)，本文取C(h)=1。

假設(shè)各個(gè)分區(qū)材料為均質(zhì)，根據(jù)常見(jiàn)的填料類(lèi)型，參考文獻(xiàn)[10]和相關(guān)經(jīng)驗(yàn)選用對(duì)應(yīng)土質(zhì)的模型參數(shù)值，見(jiàn)表1。

2 滲流計(jì)算結(jié)果與分析

2.1 含水率變化特征

經(jīng)過(guò)數(shù)值計(jì)算，常水位下原堤防和加寬堤防濕度場(chǎng)云圖見(jiàn)

表1 各分區(qū)材料滲流計(jì)算參數(shù)Tab.1 Soil hydraulic parameters of material for different zones

圖2。圖2表明，水從臨水側(cè)邊坡滲入堤身內(nèi)部，濕度等勢(shì)線(xiàn)由臨水側(cè)向背水側(cè)呈左高右低的斜線(xiàn)分布，這是由于河流水從臨水側(cè)邊坡入滲，隨著距離的增加，浸潤(rùn)線(xiàn)逐漸下降，含水率逐漸降低。舊堤身被飽和濕度線(xiàn)分為上、下兩部分，下部為濕度飽和區(qū)，上部為濕度漸變區(qū)，漸變區(qū)內(nèi)濕度由下而上逐漸減小，濕度梯度也由下而上逐漸減小。加寬堤身含水率等勢(shì)線(xiàn)在交界面處不連續(xù)，新堤身未達(dá)到飽和狀態(tài)，但濕度分布與舊堤身相似，整體濕度小于舊堤身。

繪制距離堤頂以下水平30 cm層位、80 cm層位和150 cm層位濕度變化圖，見(jiàn)圖3。從圖3可以看出，從臨水側(cè)至背水側(cè)含水率曲線(xiàn)逐漸降低，土體含水量逐漸減少，這種趨勢(shì)無(wú)論在舊堤身還是新堤身都一致，舊堤整體濕度明顯大于新堤。另外在新、舊堤身交界面上濕度急劇下降，水分從舊堤身入滲至相對(duì)而言含水率較小、滲透系數(shù)較大的新堤身中，因此含水率出現(xiàn)急劇變化。交界面處30 cm層位、80 cm層位和150 cm層位含水率分別下降了15.2%、14.8%和13.7%，表明交界面處高度越高、含水率越低的層位，含水率的下降越多，但這種含水率變化的差異不大。不論是加寬堤防含水率等勢(shì)線(xiàn)不連續(xù)，還是交界面處含水率的急劇變化，都是由于新、舊堤身的填料差異而引起的。

圖2 常水位作用下濕度場(chǎng)云圖Fig.2 Maps of humidity field at ordinary water level

圖3 堤身各層位濕度分布Fig.3 Humidity distribution at different widened levee horizons

2.2 總水頭變化特征

圖4為常水位下原堤防和加寬堤防總水頭等勢(shì)線(xiàn)云圖，比較發(fā)現(xiàn)，總水頭等勢(shì)線(xiàn)在靠近交界面處發(fā)生轉(zhuǎn)折，表明靠近交界面總水頭減小。究其原因，新填加寬堤身土質(zhì)滲透系數(shù)大于舊路基，在新舊路基結(jié)合面處，相對(duì)含水率低的新堤將含水率高的原堤防進(jìn)行吸濕作用，從而導(dǎo)致交界面處原堤防含水率降低，水頭減小，這與含水率變化特征分析一致。新填堤身由于滲透系數(shù)大，距河流水位較遠(yuǎn)，整體濕度小，水頭等勢(shì)線(xiàn)梯度小于原堤身。

圖5為原堤防堤頂外邊緣(1斷面)、加寬堤防中線(xiàn)(2斷面)、加寬路基堤頂外邊緣(3斷面)在堤身高度范圍內(nèi)總水頭變化曲線(xiàn)，和在相同位置原堤防總水頭變化線(xiàn)對(duì)比，顯示堤防加寬后靠近交界面處的水頭有急劇下降的過(guò)程，而原堤防在上述斷面總水頭基本保持一致。交界面處總水頭減少率隨著斷面而降低。交界面處1斷面總水頭下降了6.8%，2斷面下降了4.6%，3斷面下降了1.3%，表明總水頭的降低率隨著橫向距離的增加和高度的降低而減少，這也與上述含水率在交界面較高處下降顯著的特性一致。根據(jù)含水率變化特征，堤身范圍內(nèi)濕度自上而下降低，交界面處含水率差值自上而下降低，較小的差值使得新堤吸濕程度下降，水分遷移效果不顯著，從而導(dǎo)致總水頭自上而下降低的程度減小。

圖4 常水位作用下總水頭云圖Fig.4 Maps of total head at ordinary water level

2.3 滲透坡降變化特征

圖6為常水位下原堤身和加寬堤身滲透坡降云圖，可以看出加寬堤身在三處位置出現(xiàn)較大的滲透坡降：臨水側(cè)邊坡水位高度處、交界面靠近頂面處、新填堤身坡腳處，而舊堤僅在臨水側(cè)邊坡水位高度處出現(xiàn)最大滲透坡降。根據(jù)總水頭變化云圖，這三處總水頭等勢(shì)線(xiàn)分布最為密集，導(dǎo)致單位距離內(nèi)的水頭差增大，從而使得滲透坡降出現(xiàn)較大值。對(duì)于交界面靠近頂面處和加寬路基坡腳處出現(xiàn)較大的滲透坡降，根據(jù)含水率變化特征和總水頭變化特征分析可知，該位置自原堤身向新填堤身方向含水率和總水頭均急劇下降，從而引起滲透坡降出現(xiàn)較大值。

圖5 堤防加寬前后交界面處總水頭變化Fig.5 The variation of total head at interface before and after levee widened

圖6 常水位作用下滲透坡降云圖Fig.6 Maps of seepage slope at ordinary water level

重點(diǎn)分析交界面處的滲透坡降變化情況，圖7為堤防加寬前后交界面處滲透坡降變化率。整個(gè)交界面處滲透坡降均有較大變化，三個(gè)斷面處出現(xiàn)的最大滲透坡降值較原堤防的滲透坡降分別增大了923.1%、589.8%和150.1%，同樣是自上而下變化率有所減少，與含水率和總水頭變化率一致。

為進(jìn)一步了解加寬堤防交界面和坡腳處產(chǎn)生較大滲透坡降的原因，將新堤材料的滲透系數(shù)分別修改為與舊堤相同(5×10-6m/s)和舊堤的1/10(5×10-7m/s)，在其他條件不變的情況下重新計(jì)算，計(jì)算結(jié)果見(jiàn)圖8。比較圖8與圖6(b)可以看出，隨著新、舊堤滲透系數(shù)比值的減小，交界面和坡腳處的滲透坡降明顯減小，說(shuō)明該位置的滲透坡降與新、舊堤材料的滲透系數(shù)比值相關(guān)，比值越大，滲透坡降越大，反之亦然。

圖7 堤防加寬前后交界面處滲透坡降變化Fig.7 The variation of seepage slope at interface before and after levee widened

圖8 新堤材料不同滲透系數(shù)時(shí)滲透坡降云圖Fig.8 Maps of seepage slope at different permeability coefficient of new levee materials

3 變形計(jì)算結(jié)果與分析

認(rèn)為舊堤防已經(jīng)固結(jié)穩(wěn)定，在未加寬前本身不再出現(xiàn)變形，數(shù)值計(jì)算時(shí)首先對(duì)舊堤防進(jìn)行原始初始應(yīng)力場(chǎng)分析，然后添加新堤防進(jìn)行應(yīng)力-應(yīng)變和滲流耦合分析。

3.1 豎向位移變化特征

計(jì)算不同水位下加寬堤防在滲流和應(yīng)力耦合作用下的變形特性，見(jiàn)圖9。從圖9可以看出，加寬堤防的最大沉降位置在新填堤身的外側(cè)堤肩處。不同洪水位下最大豎向位移量在0.041 7～0.043 4 m(見(jiàn)圖10)，約為堤高的0.8%左右。從數(shù)值可見(jiàn)盡管水位對(duì)最大沉降量的影響不大，不同水位下的最大沉降量差值僅為0.001 67 m，但可以發(fā)現(xiàn)洪水位在2 m時(shí)沉降量最小，當(dāng)水位下降或者升高時(shí)沉降量都會(huì)有所增加。

圖9 不同水位下加寬堤防豎向位移分布云圖Fig.9 Maps of vertical displacement of widened levee at different water levels

圖10 不同水位下加寬堤防最大沉降量Fig.10 Maximum settlement of widened levee at different water levels

由于不同水位下加寬堤防豎向位移分布云圖相似，僅以3 m水位為例，繪制加寬堤防整個(gè)基底的豎向變形曲線(xiàn)，見(jiàn)圖11。堤防加寬后整個(gè)基底的變形如馬鞍形狀，基底的最大變形出現(xiàn)在近新、舊堤防的相接處，方向向下，大小約為0.013 m；原堤基底中部和新堤基底外邊緣處變形方向向上，大小約為0.004 m。這是由于在新填筑堤防自重作用下堤基被壓縮，同時(shí)孔隙水壓力場(chǎng)有所變化，導(dǎo)致新填堤防堤基兩側(cè)應(yīng)力較大，從而出現(xiàn)向上隆起現(xiàn)象。

圖11 加寬堤防基底豎向變形量Fig.11 Vertical deformation of widened levee base

3.2 水平位移變化特征

圖12為不同洪水位下加寬堤防水平位移分布云圖?？梢钥闯龀霈F(xiàn)舊堤向左側(cè)、新堤向右側(cè)的位移運(yùn)動(dòng)趨勢(shì)。舊堤的最大水平位移大約出現(xiàn)在加寬后堤防中線(xiàn)的下部堤身處，最大位移量為0.006～0.009 m；向右的最大水平位移量出現(xiàn)在新填堤防堤底以下的堤基內(nèi)，數(shù)值為0.007～0.008 m?？梢?jiàn)水平位移量值很小，不同水位下平均最大水平位移量?jī)H為最大豎向位移量的19%，水位高低對(duì)最大水平位移量的影響很小，以3斷面為例，繪制水平位移分布圖，如圖13。

圖12 不同水位下加寬堤防水平位移分布云圖Fig.12 Maps of horizontal displacement of widened levee at different water levels

圖13 不同水位下斷面3水平位移隨深度變化Fig.13 The variation of horizontal displacement vs. depth of the third section at different water levels

從圖13可以發(fā)現(xiàn)，水位高低對(duì)新堤和堤基上部的水平位移幾乎沒(méi)有影響，而對(duì)下側(cè)堤基的有一定影響。盡管數(shù)值不大，高水位時(shí)堤基下部有持續(xù)向右的位移趨勢(shì)，當(dāng)水位降低時(shí)，這種向右的位移趨勢(shì)逐漸又變向左側(cè)，且水位越低，這種趨勢(shì)越明顯。

4 結(jié) 語(yǔ)

(1)在常水位作用下，由于新舊填料透水性函數(shù)參數(shù)不同，加寬堤防的滲流特征發(fā)生明顯變化。非飽和狀態(tài)下在相同基質(zhì)吸力時(shí)，新堤滲透系數(shù)較大，導(dǎo)致交界面處新堤吸濕舊堤水分，降低了交界面處的含水率和總水頭，并產(chǎn)生較大滲透坡降。

(2)新舊堤防交界面處的含水率變化率、總水頭變化率和滲透坡降變化率隨著路基高度的增加而增加，尤其是交界面處和加寬堤防坡腳處的滲透坡降的變化率十分 顯著，這與新、舊填料滲透系數(shù)比值相關(guān)。

(3)加寬堤防的最大沉降位置出現(xiàn)在新堤的外側(cè)堤肩處，沉降量與洪水位高度有關(guān)系，在水位升高或降低的過(guò)程中最大沉降量存在一個(gè)最小值，新堤填筑后整個(gè)基底的變形如馬鞍形狀；堤防加寬后會(huì)出現(xiàn)舊堤堤身和新堤堤基向相反水平方向位移的趨勢(shì)，水位高度對(duì)位移量的影響較小。

(4)針對(duì)上述沿河加寬堤防的滲流特征，若交界面處的計(jì)算滲透坡降大于允許值時(shí)，可采取適當(dāng)措施，如采用與舊堤滲透系數(shù)相近的材料等，防止可能在交界面處發(fā)生滲透破壞。在加寬堤防前在適當(dāng)?shù)奈恢眠M(jìn)行地基加固、在新堤填筑過(guò)程中設(shè)置土工織物等措施保證加寬堤防的沉降穩(wěn)定性。

□