修冬紅

(河北省水利水電勘測設(shè)計(jì)研究院集團(tuán)有限公司，石家莊 050000)

1 概 述

堤壩是防洪系統(tǒng)的基本組成部分[1]。截至2018年，我國全國已建成5級(jí)及以上江河堤壩31.2×104km，累計(jì)達(dá)標(biāo)堤壩21.8×104km，達(dá)標(biāo)率為69.8%。其中，1級(jí)、2級(jí)達(dá)標(biāo)堤壩長度為3.4×104km，達(dá)標(biāo)率為80.5%。堤壩保護(hù)人口6.3億人，保護(hù)耕地4.1×104hm2。

堤壩變形是由于水的液壓、水通過堤壩及其底土的過濾作用以及堤壩材料的水化和干燥引起的[2-3]，導(dǎo)致抗剪強(qiáng)度降低，某些部位產(chǎn)生不均勻變形[4-5]、裂縫、下滑等，堤壩變形分為滑坡、崩岸等幾種形式。在洪水期間，在水位上升的影響下，這一過程可能還會(huì)加劇[6]。在水壓變化的條件下，變形過程的發(fā)展往往會(huì)導(dǎo)致堤壩的破壞。

本文針對(duì)水位高達(dá)4 m時(shí)，水壓增加和減少的情況下堤壩的數(shù)值模擬結(jié)果；并分析洪水條件下，堤壩中孔隙壓力的變化規(guī)律。同時(shí)，將數(shù)值計(jì)算的結(jié)果與堤壩監(jiān)測IT系統(tǒng)[7-8]下建造的試驗(yàn)性蓄水設(shè)施內(nèi)的現(xiàn)場測量結(jié)果進(jìn)行了比較。

2 研究區(qū)域的特點(diǎn)

研究區(qū)域由208 m×58 m×4.5 m(長×寬×高)的蓄水設(shè)施形成兩段相互連接的堤壩，見圖1，該結(jié)構(gòu)體積約為28.5 m3。

圖1 試驗(yàn)蓄水設(shè)施位置及平面圖

堤岸坡度可變，水側(cè)為1∶2和1∶2.5，陸側(cè)為1∶2。堤壩東側(cè)為粉質(zhì)黏土，西側(cè)為粉砂，過濾系數(shù)分別為1.83×10-5和5.24×10-5m/s。底土由粉質(zhì)黏土和黏質(zhì)砂、砂和礫石組成，與底土接觸的堤壩用不透水的膨潤土隔離。地下水條件通過鉆孔、動(dòng)態(tài)和靜態(tài)探測確定。

在現(xiàn)場測量中，對(duì)水位升降對(duì)堤壩行為的影響進(jìn)行模擬，并通過孔隙壓力傳感器測量水壓的變化。

3 方法論

3.1 計(jì)算方案

計(jì)算過程的示意圖見圖2。在計(jì)算中，假設(shè)水位在96 h內(nèi)持續(xù)上升至4 m，并在120 h內(nèi)下降。借助FLAC 2D軟件[9-10]，在平面應(yīng)變狀態(tài)下進(jìn)行數(shù)值模擬。選取I - I′截面和位于堤壩下部的4個(gè)監(jiān)測點(diǎn)，分析水位的增減對(duì)孔隙壓力變化的影響(圖1)。A點(diǎn)和B點(diǎn)位于蓄水設(shè)施東側(cè)的堤壩上，C點(diǎn)和D點(diǎn)位于蓄水設(shè)施西側(cè)的堤壩上。

3.2 物理模型

物理模型沿橫斷面I - I′制備，橫斷面貫穿試驗(yàn)蓄水設(shè)施中部，見圖3。假設(shè)垂直框架高11 m，水平框架長90 m。表1為6個(gè)巖土層的材料特性，第I層對(duì)應(yīng)于試驗(yàn)蓄水設(shè)施東側(cè)或西側(cè)的堤壩材料，其余5層位于試驗(yàn)蓄水設(shè)施的底土中。

圖2 數(shù)值模擬示意圖

圖3 沿橫截面I - I′的底土蓄水物理模型

表1 堤壩和底土層材料常數(shù)

3.3 計(jì)算模型

在計(jì)算中，采用莫爾-庫侖強(qiáng)度準(zhǔn)則，建立彈塑性力學(xué)模型。在堤壩水側(cè)施加靜水壓力，堤壩的陸側(cè)允許水流。采用0.25 m的規(guī)則網(wǎng)格對(duì)模型進(jìn)行離散，其中堤壩為0.5 m，第I-第III層為0.5 m，第IV和第V層為1 m。分析水位上升的影響：①48 h內(nèi)上升至2 m；②96 h內(nèi)上升至4 m；③120 h內(nèi)恢復(fù)至初始狀態(tài)。

4 結(jié)果與分析

孔隙壓力、剪應(yīng)力和水平位移等值線變化形式的數(shù)值模擬結(jié)果見圖4-圖6。水位上升至2 m，導(dǎo)致堤壩底部的孔隙壓力增加至20 kPa，見圖4(a)。在假定的洪水持續(xù)時(shí)間內(nèi)，試驗(yàn)蓄水設(shè)施陸側(cè)的剪應(yīng)力達(dá)到-15～20 kPa，見圖4(b)。堤壩開始變形，水平位移達(dá)到7 mm，見圖4(c)。在這種情況下，堤壩中的最大水流速度為0.33 m/d。

在這些條件下，水位上升至4 m，孔隙壓力增加至40 kPa，堤壩底部的變形過程增加，見圖5(a)。在假定的洪水持續(xù)時(shí)間內(nèi)，試驗(yàn)蓄水設(shè)施陸側(cè)的剪應(yīng)力達(dá)到-30～40 kPa，見圖5(b)。

注：黑色箭頭表示最大位移區(qū)(7 mm)。圖4 水位上升至2 m高度時(shí)堤壩內(nèi)的孔隙壓力場、剪應(yīng)力和水平位移

注：黑色箭頭表示最大位移區(qū)(14 mm)。圖5 水位上升至4 m高度時(shí)堤壩內(nèi)的孔隙壓力場、剪應(yīng)力和水平位移

在堤壩中，孔隙壓力值顯著降低至幾個(gè)kPa，這可能是由于水流速度增加至3.5 m/d。在堤壩的陸側(cè)，水平位移達(dá)到14 mm見圖5(c)。在堤壩的水側(cè)部分，觀察到剪切應(yīng)力增大，達(dá)到40 kPa，見圖5(b)。

由于水位降至1 m，孔隙壓力達(dá)到20 kPa，堤壩底部的變形過程減少，見圖6(a)。在假定的洪水持續(xù)時(shí)間內(nèi)，試驗(yàn)蓄水設(shè)施陸側(cè)的剪應(yīng)力為-15～15 kPa，見圖6(b)。在這些條件下，達(dá)到的最大位移為3 mm，見圖6(c)。堤壩內(nèi)的水流速度降至約0.13 m/d，見圖6(b)。

注：黑色箭頭表示最大位移區(qū)(3 mm)。圖6 水位下降至1 m高度時(shí)堤壩內(nèi)的孔隙壓力場、剪應(yīng)力和水平位移

一般來說，數(shù)值模型中的孔隙壓力大致隨水位變化而增大或減小，見圖7。通過比較數(shù)值模型和現(xiàn)場試驗(yàn)設(shè)施中A、B、C和D 4個(gè)監(jiān)測點(diǎn)的孔隙壓力變化，可以觀察到這些變化的特征是相似的，計(jì)算得到的孔隙壓力最大值與實(shí)測值基本一致。與數(shù)值模型中的堤壩響應(yīng)相比，試驗(yàn)堤壩僅在后期響應(yīng)才會(huì)出現(xiàn)差異(差異僅出現(xiàn)在試驗(yàn)堤壩的后期響應(yīng)與數(shù)值模型堤壩的響應(yīng)之間)(圖7)。據(jù)推測，這些差異與水流通過試驗(yàn)堤壩時(shí)的過濾侵蝕效應(yīng)有關(guān)。在水位下降的情況下，各監(jiān)測點(diǎn)的計(jì)算孔隙壓力值與實(shí)測孔隙壓力值相似。

在現(xiàn)場測量中，隨著水位由2 m上升至4 m，堤壩6個(gè)巖土層上的孔隙壓力均顯著提高。當(dāng)水位下降至1 m時(shí)，孔隙壓力顯著降低。需要強(qiáng)調(diào)的是，水位升降時(shí)間短，限制了堤壩變形的發(fā)展。在現(xiàn)場試驗(yàn)條件下，粉砂質(zhì)非對(duì)稱大堤更易變形，孔隙壓力值較高，約為5～7 kPa，見圖(8)。

圖7 在數(shù)值模型和現(xiàn)場試驗(yàn)中4個(gè)監(jiān)測點(diǎn)(A-D)在水位升高4 m時(shí)的孔隙壓力變化圖

圖8 I-VI剖面孔隙壓力變化在水位上升至2、4 m和下降至1 m的現(xiàn)場測量

5 結(jié) 論

通過水壓對(duì)堤壩影響的數(shù)值模擬結(jié)果，結(jié)論如下：

1) 在水位升高的情況下，數(shù)值模型中孔隙壓力隨時(shí)間的變化及其最大值與現(xiàn)場試驗(yàn)堤壩相似。兩個(gè)堤壩的行為差異僅出現(xiàn)在試驗(yàn)堤壩的后期響應(yīng)與數(shù)值模型堤壩的響應(yīng)之間。據(jù)推測，這些差異與水流通過實(shí)驗(yàn)堤壩時(shí)的過濾侵蝕效應(yīng)有關(guān)。

2) 在水位下降的情況下，各監(jiān)測點(diǎn)的計(jì)算和實(shí)測孔隙壓力隨時(shí)間的變化相似。

3) 在數(shù)值模型中，當(dāng)水位上升至4 m高度時(shí)，最大位移發(fā)生在堤壩陸側(cè)的底部。在試驗(yàn)堤壩中，觀測到該邊坡表面的出水情況，計(jì)算的剪切應(yīng)力為30 ～40 kPa，水平位移增加至14 mm。

4) 研究結(jié)果表明，水位升高引起的水壓對(duì)孔隙壓力變化有顯著影響，并可能導(dǎo)致堤壩破壞。