平原水庫土石壩飽和砂土地基地震液化分析

胡南雄 任旭華 張繼勛

(河海大學(xué) 水利水電學(xué)院， 南京 210098)

我國是一個(gè)地震發(fā)生頻率較高的國家，地震活動(dòng)的分布范圍十分廣泛且強(qiáng)度較大，基本烈度在Ⅶ度以上區(qū)域的比例在60%以上[1]，大量的水庫工程都建造在地震區(qū)域．其中平原水庫工程一般位于大江、大河下游沖積平原地區(qū)，這類水庫大壩的地基條件往往較差，多數(shù)為軟土和透水地基，由此經(jīng)常產(chǎn)生壩基滲漏、壩基地震液化等普遍性問題．由地震引發(fā)的飽和砂土液化往往會(huì)導(dǎo)致土石壩壩基沉陷，進(jìn)而引起壩體裂縫、滑坡等嚴(yán)重的區(qū)域性地質(zhì)災(zāi)害[2]．對(duì)于地震引發(fā)的砂土液化問題，李颯[3]等以某典型平原水庫為背景，在室內(nèi)及現(xiàn)場(chǎng)試驗(yàn)的基礎(chǔ)上對(duì)可液化土層進(jìn)行判別并考慮滲流對(duì)可液化層土體的影響．嚴(yán)祖文[4]等對(duì)水庫大壩地基地震液化特性進(jìn)行研究并用動(dòng)力有限元軟件分析土體的液化過程．陳育民[5]等運(yùn)用二維有限元程序?qū)δ车袒吧喜康谭赖卣鹨夯冃芜M(jìn)行分析，探討了堤防的殘余變形、有效應(yīng)力等動(dòng)力響應(yīng)．但是對(duì)于不同庫水位下水庫大壩壩基可液化土層在地震作用下的變化過程目前鮮有研究．本文運(yùn)用FLAC3D軟件對(duì)某中型水庫大壩壩基可液化土層進(jìn)行數(shù)值模擬，分析了不同庫水位作用下，壩基可液化土層的孔隙水壓力變化過程及壩體典型節(jié)點(diǎn)的動(dòng)力響應(yīng)，研究結(jié)果可為類似工程提供技術(shù)支持．

1 研究方法

1.1 動(dòng)孔壓模型

FLAC3D可以模擬砂土在動(dòng)力作用下的孔壓積累直至土體液化的過程，其內(nèi)置的Finn模型可以描述孔壓積累效應(yīng)．具體方程如下：

(1)

(2)

式中，γ為土體的剪應(yīng)變，εvd為累積體積應(yīng)變，C1、C2、C3和C4為模型常量．

1.2 液化判別標(biāo)準(zhǔn)

FLAC3D模擬地震作用下土體液化時(shí)常采用超孔壓比來描述，公式如下[7]

(3)

2 計(jì)算模擬

2.1 工程概況

山東省某中型平原水庫庫區(qū)地貌上屬微傾斜低平原區(qū)的沖擊-海積平原亞區(qū)，地勢(shì)低平，易受海潮影響，由于海水浸漬，多濕洼地，土壤鹽漬嚴(yán)重．水庫圍壩為復(fù)合土工膜防滲體斜墻砂壤土均質(zhì)壩，壩頂高程為15 m，地面高程為3.5 m，水庫設(shè)計(jì)最高蓄水位為12.5 m，死水位為3.9 m．

圍壩典型斷面如圖1所示．壩址區(qū)的地震動(dòng)峰值加速度為0.15g，相應(yīng)的地震基本烈度為Ⅶ度．在工程正常運(yùn)用時(shí)1、2和3層土壤處于飽和狀態(tài)，粒徑小于 0.005 mm 的顆粒pc含量質(zhì)量百分率均小于相應(yīng)于地震設(shè)防烈度Ⅶ度時(shí)的界限值 16%，在地震基本烈度Ⅶ度條件下，存在振動(dòng)液化的可能性．利用勘察場(chǎng)地內(nèi)鉆孔的標(biāo)準(zhǔn)貫入試驗(yàn)對(duì)可液化土層復(fù)判，綜合評(píng)定1、2和3土層為可液化土層，其余各層為非液化土層．

圖1 地質(zhì)剖面圖

2.2 計(jì)算模型

本文取圍壩典型斷面3+050處50m長(zhǎng)壩段進(jìn)行建模和計(jì)算．模型上下游面距離壩體中心都為100 m，地基厚度為40 m．模型原點(diǎn)設(shè)在壩體底部中心，X軸順河向，指向下游為正；Y軸沿壩軸線，指向左岸為正；Z軸為豎直向，向上為正．壩體高度為11.5 m，地面高程為0 m，正常蓄水位為9 m，死水位為0.4 m．上游壩坡面設(shè)有復(fù)合土工膜防滲體，上游壩腳位置還設(shè)有防滲墻．壩體及地基各土層分布如圖2所示．

圖2 計(jì)算網(wǎng)格圖

壩體采用分層填筑，其中考慮1層細(xì)砂土、2層粉砂土和3層粉砂土為可能液化土層．同時(shí)為了便于觀測(cè)，在2層粉砂土內(nèi)及壩體各處設(shè)置了數(shù)個(gè)監(jiān)測(cè)點(diǎn)，用于監(jiān)測(cè)地震作用下該點(diǎn)所在單元的動(dòng)力響應(yīng)情況．

2.3 材料參數(shù)

根據(jù)現(xiàn)場(chǎng)勘測(cè)資料并結(jié)合室內(nèi)試驗(yàn)測(cè)得各土層及壩體的物理力學(xué)參數(shù)見表1.水的密度為1 000 kg/m3，體積模量為200 MPa，抗拉強(qiáng)度為0.本文在模擬可液化土層在地震作用下的液化過程時(shí)采用FLAC3D內(nèi)置的Finn本構(gòu)模型，該模型的模型常數(shù)近似取為相對(duì)密度為45%的結(jié)晶二氧化硅常數(shù)：C1=0.8，C2=0.79，C3=0.45，C4=0.73．

2.4 靜力計(jì)算

本文先對(duì)地基土層進(jìn)行初始地應(yīng)力場(chǎng)計(jì)算，而后進(jìn)行壩體分層填筑靜力計(jì)算，最后對(duì)模型設(shè)置初始孔壓場(chǎng)，得到準(zhǔn)確的初始應(yīng)力場(chǎng)和初始孔壓場(chǎng)．初始孔壓云圖及初始地應(yīng)力云圖如圖3所示．由理論計(jì)算可得，在正常蓄水位狀況下模型上游底部的孔壓力為490 kPa，壩基初始應(yīng)力為666.6 kPa；由FLAC3D計(jì)算的值分別為482.5 kPa和670.4 kPa．因此可以認(rèn)為靜力計(jì)算是正確的．

表1 土層物理力學(xué)參數(shù)

圖3 地基初始應(yīng)力及初始孔壓力云圖(單位：Pa)

2.5 動(dòng)力輸入

為吸收地震波在邊界上的反射，本文對(duì)模型施加自由場(chǎng)邊界，施加后模型四周將生成二維和一維網(wǎng)格．同時(shí)，在模型底部邊界設(shè)定為靜態(tài)邊界條件，消除地震波在模型底部的反射．本文采用局部阻尼進(jìn)行動(dòng)力計(jì)算，土體的阻尼比按經(jīng)驗(yàn)方法取為5%，故本工程的局部阻尼取為0.157(局部阻尼為0.05π)．

庫區(qū)的地震動(dòng)峰值加速度為0.15g，相應(yīng)的地震烈度為Ⅶ度．采用國際上常用的EI-Centro南北方向(X向)，截取地震振幅較大的30 s作為此次地震輸入波，并對(duì)地震波進(jìn)行濾波和基線校正處理后將其峰值加速度調(diào)為0.15g，調(diào)整后的加速度如圖4所示．

圖4 校正后的EI-Centro水平向加速度

3 計(jì)算結(jié)果及分析

3.1 可液化土層超孔壓比分析

圖5～8為地震結(jié)束時(shí)，庫區(qū)處于正常蓄水位和死水位時(shí)模型的超孔壓比云圖．嚴(yán)格意義上來說，只有當(dāng)超孔壓比EPWPR=1時(shí)，表示可液化土體完全液化．但是由文獻(xiàn)[8，9]及實(shí)際工程經(jīng)驗(yàn)可知，當(dāng)超孔壓比未達(dá)到1時(shí)，飽和土就可能發(fā)生部分液化，這同樣會(huì)給建筑物帶來安全隱患．所以，本文將超孔壓比EPWPR=0.6作為土體液化和不液化的臨界值：超孔壓比大于0.6視為液化，低于0.6為不液化．

圖5 正常蓄水位和死水位下模型超孔壓比

由圖5(a)可知，在正常蓄水位情況下，地震結(jié)束時(shí)刻上游距壩體較遠(yuǎn)的可液化土層發(fā)生局部液化，而上游壩腳附近以及壩體下方的可液化土層基本不發(fā)生液化，下游壩腳以及下游距離壩體較遠(yuǎn)的可液化土層也發(fā)生了局部液化．這是因?yàn)楫?dāng)庫水位為正常蓄水位時(shí)，施加在迎水面復(fù)合土工膜防滲體上的水荷載較大，使得上游壩基可液化土層的初始圍壓較大，該處可液化土體不易液化；壩體自身的重力作用使壩體下方可液化土層的孔壓難以達(dá)到上覆壓重，因此在地震作用下不易發(fā)生液化．當(dāng)上游水位較高時(shí)，在滲流的作用下，下游地下水位抬高，壩體的浸潤(rùn)線升高，使得下游壩腳以及下游距離壩體較遠(yuǎn)的可液化土層的初始有效圍壓減小，土體易發(fā)生液化．

如圖5(b) 所示在死水位情況下，地震結(jié)束時(shí)刻上游距壩體較遠(yuǎn)的可液化土層發(fā)生較大面積的液化，發(fā)生液化的區(qū)域距離地面更近且液化面積較正常蓄水位大；上游壩腳附近也發(fā)生了較小面積的液化，壩體下方及下游壩腳附近的可液化土層基本不發(fā)生液化，下游遠(yuǎn)離壩體的可液化土層發(fā)生小面積液化．這是由于在死水位時(shí)，作用在上游迎水面的水荷載基本消失了，上游壩基可液化土體的初始圍壓大大減小，且該區(qū)域土體受到來自壩體初始剪應(yīng)力向上分力作用，有效上覆壓力進(jìn)一步減小，在地震作用下就容易發(fā)生液化．低水位情況下，下游壩腳附近可液化土層的初始有效圍壓較大，土體不易發(fā)生液化．

圖6為典型單元E1、E2和E3在不同水位下的超孔壓比時(shí)程曲線(具體監(jiān)測(cè)單元見圖2)．E1單元的超孔壓比在地震開始發(fā)生的5s內(nèi)迅速增加，之后超孔壓比保持在相對(duì)穩(wěn)定的狀態(tài)且在正常蓄水位和死水位下的超孔壓比都超過了0.6，該處壩基可液化土體發(fā)生了液化，這個(gè)結(jié)果與超孔壓比云圖顯示的結(jié)果一致．E2單元位于壩體正下方，在正常蓄水位和死水位情況下超孔壓比都遠(yuǎn)小于0.6，可知該處土體不會(huì)發(fā)生液化．同理，E3單元在正常蓄水位狀況下超孔壓比大于0.6，在死水位情況下小于0.6，與超孔壓比云圖結(jié)果一致．

圖6 各典型單元超孔壓比

3.2 各節(jié)點(diǎn)加速度地震響應(yīng)分析

壩體各監(jiān)測(cè)節(jié)點(diǎn)加速度時(shí)程曲線如圖7所示(具體監(jiān)測(cè)節(jié)點(diǎn)見圖2)．從圖7(a) 可以看出，在地震發(fā)生的前5 s內(nèi)，正常蓄水位和死水位狀況下節(jié)點(diǎn)N1的加速度分別達(dá)到了各自最大值2.382 m/s2和2.409 m/s2，是地震峰值加速度的1.588倍和1.606倍．隨著地震的持續(xù)作用，死水位情況下壩體上游壩腳附近可液化土層發(fā)生局部液化，地面節(jié)點(diǎn)的運(yùn)動(dòng)受到液化區(qū)域的影響，液化砂土吸收了地震波的高頻成分，減小了地面節(jié)點(diǎn)的加速度[10]，因此在地震波5～30 s區(qū)間節(jié)點(diǎn)N1的加速度明顯比輸入波小．從圖7(b) 可見，正常蓄水位和死水位狀況下壩頂?shù)湫凸?jié)點(diǎn)N2的最大加速度比上下游壩腳處的節(jié)點(diǎn)都要大，分別達(dá)到了2.913 4 m/s2和2.824 2 m/s2，是地震峰值加速度的1.942倍和1.883倍，說明壩頂位置存在鞭梢效應(yīng)．從圖7(c) 可以看出對(duì)于下游壩腳處的節(jié)點(diǎn)N3，在地震波7～30 s區(qū)間正常蓄水位情況下節(jié)點(diǎn)N3的加速度較輸入波小，這同樣是因?yàn)樵撎幍目梢夯翆影l(fā)生了液化，液化區(qū)域?qū)Φ孛孢\(yùn)動(dòng)產(chǎn)生了影響．

圖7 各典型節(jié)點(diǎn)加速度地震響應(yīng)

3.3 殘余變形分析

圖8～9是地震作用結(jié)束時(shí)刻庫區(qū)處于不同水位時(shí)豎向殘余變形和水平殘余變形云圖．從圖8(a)可以看出，當(dāng)庫水位為死水位時(shí)，地震結(jié)束時(shí)刻上游和下游壩腳處都產(chǎn)生了向上的豎向殘余變形，即壩體周圍土體出現(xiàn)了不同程度的隆起現(xiàn)象．壩體的豎向殘余變形比壩基土層大，從壩頂?shù)綁蔚棕Q向變形逐漸減少，其中壩頂?shù)呢Q向殘余變形最大，約為18 cm．從圖8(b)可以看出，當(dāng)庫水位為正常蓄水位時(shí)，地震結(jié)束時(shí)刻下游壩腳也產(chǎn)生了向上的豎向殘余變形．模型的最大豎向變形同樣發(fā)生在壩頂，變形趨勢(shì)與死水位情況較為相似．

圖8 正常蓄水位和死水位下模型豎向殘余變形

如圖9(a)所示，在庫區(qū)水位處于死水位時(shí)，壩體的水平變形呈兩極分化趨勢(shì)：上游壩腳的水平位移指向上游，最大變形約為22 cm；下游壩腳的水平位移指向下游，最大變形約為24 cm．如圖9(b)所示，在正常蓄水位情況下壩體的整體水平位移指向下游，水平位移從上游到下游逐漸增大，最大水平變形發(fā)生在下游壩腳處，達(dá)到了26 cm．這可能是因?yàn)閹靺^(qū)處于高水位時(shí)，作用在壩體迎水面復(fù)合土工膜防滲體上的水荷載較大，使壩體向下游移動(dòng)．

圖9 正常蓄水位和死水位下模型水平殘余變形

綜上所述，不管庫區(qū)處于正常蓄水位還是死水位，下游壩腳處都出現(xiàn)了隆起現(xiàn)象，壩頂中心出現(xiàn)了較大的沉陷位移，壩體的水平位移較豎直位移大，因此需要對(duì)上下游壩腳處等其他薄弱部位進(jìn)行加固以減小其變形．

4 結(jié) 論

本文采用有限差分軟件FLAC3D中的Finn本構(gòu)模型，以某中型平原水庫土石壩為背景，討論了不同庫水位下壩基地震液化問題以及土石壩變形，得出如下結(jié)論：

1)FLAC3D可以較好地模擬地震作用下可液化土層孔隙水壓力的變化過程，其中超孔壓比可以作為判別土體液化的依據(jù)，并且具有很強(qiáng)的直觀性．

2)壩基可液化土層同一深度的不同單元在地震作用下孔壓的發(fā)展水平各不相同．庫水位較高時(shí)，迎水面上的有效荷載大，一定程度上抑制了上游可液化土層的液化，但是由于下游地下水位抬高，減小了下游壩腳處可液化土體的上覆有效應(yīng)力，使該處的土體發(fā)生液化．庫水位較低時(shí)，迎水面有效荷載大大減小，上游可液化土層易發(fā)生液化；下游壩腳處可液化土層初始有效應(yīng)力較大，所以該處土體不易發(fā)生液化．

3)壩基液化區(qū)域會(huì)吸收地震波的高頻成分，一定程度上減小了地表節(jié)點(diǎn)的加速度．壩頂?shù)募铀俣确逯底畲螅霈F(xiàn)了鞭梢效應(yīng)．

4)正常蓄水位下壩體的整體水平位移指向下游，死水位情況下壩體的水平位移呈現(xiàn)兩極化；不同水位下的壩頂豎向殘余變形都比較大，上下游壩腳出現(xiàn)隆起現(xiàn)象．因此工程上為了保證大壩的安全必須對(duì)可液化土層進(jìn)行處理，并在地基薄弱部位和壩體周圍增加抗液化措施，保證工程在地震作用下的安全性．

參考文獻(xiàn)：

[1] 李紅軍．高土石壩地震變形分析與抗震安全評(píng)價(jià)[D]．大連：大連理工大學(xué)，2008．

[2] 鄧益兵，周 健，劉文白，等．水庫圍堤地基液化有效應(yīng)力動(dòng)力分析[J]．巖土力學(xué)，2010，32(s2)：292-296．

[3] 李 颯，劉松良，要明倫，等．地震作用下滲流對(duì)平原水庫土石壩穩(wěn)定的影響[J]．巖石力學(xué)與工程學(xué)報(bào)，2008，27(4)：858-864．

[4] 嚴(yán)祖文，李敬梅．水庫大壩地基地震液化特性及動(dòng)力有限元分析[J]．水力發(fā)電學(xué)報(bào)，2006，25(6)：45-49．

[5] 陳育民，仉文崗，劉漢龍．堤防地震液化數(shù)值模擬及動(dòng)力反應(yīng)分析[J]．解放軍理工大學(xué)(自然科學(xué)版)，2011，12(3)：239-244．

[6] 錢家歡，殷宗澤．土工原理與計(jì)算[M]．北京：中國水利水電出版社，1996．

[7] 陳育民．FLAC/FLAC3D基礎(chǔ)與工程實(shí)例[M]．北京：中國水利水電出版社，2008．

[8] 顧寶和，石兆吉．地震液化效應(yīng)的綜合評(píng)價(jià)[J]．工程地質(zhì)學(xué)報(bào)，1995，3(3)：1-10．

[9] 韓 石，貢金鑫，張艷青．地震作用下重力式碼頭地基液化及變形[J]．水利水運(yùn)工程學(xué)報(bào)，2013，35(4)：45-54．

[10] 李培振，任紅梅，呂西林，等．液化地基自由場(chǎng)振動(dòng)臺(tái)模型試驗(yàn)研究[J]．地震工程與工程振動(dòng)，2008，28(2)：171-178．