，，

(浙江工業(yè)大學(xué) 建筑工程學(xué)院，浙江 杭州 310023)

隨著沿海地區(qū)經(jīng)濟(jì)建設(shè)的高速發(fā)展、海洋資源的開(kāi)發(fā)與利用，跨海大橋、深水港、海底電纜、石油輸送管道和防波堤等離岸工程的建設(shè)越來(lái)越多。然而，這些結(jié)構(gòu)物地基在受到上部荷載的同時(shí)，還受到波浪荷載的作用，導(dǎo)致海床土體中的孔隙水壓力上升，有效應(yīng)力降低，從而造成海床地基變得不穩(wěn)定甚至發(fā)生液化。海床一旦發(fā)生液化將失去穩(wěn)定性，對(duì)海上建筑物和海底設(shè)備造成極大的破壞。由于海底環(huán)境的復(fù)雜性，這些破壞一般來(lái)說(shuō)是不可逆的。對(duì)波浪荷載作用下海床的動(dòng)力響應(yīng)和液化，國(guó)內(nèi)外學(xué)者開(kāi)展了大量的研究工作。Henkel[1]分析了波浪作用下海床的動(dòng)力響應(yīng)，Zen[2]在此基礎(chǔ)上進(jìn)一步研究了波浪作用下海床的動(dòng)力響應(yīng)和液化判別相關(guān)理論和實(shí)驗(yàn)方法。Jeng[3]將Zen的液化判別準(zhǔn)則擴(kuò)展至三維情況。Nataraja[4]分析了海床地基失效案例，對(duì)地震引起的液化分析方法進(jìn)行了修正，分析了波浪引起的海床液化與波浪周期、水深和海床強(qiáng)度等的關(guān)系。目前，我國(guó)有關(guān)土壤液化問(wèn)題的研究主要集中在地震作用下砂土的液化，其判別方式主要采用各類(lèi)規(guī)范中有關(guān)液化判別的定義，需要結(jié)合現(xiàn)場(chǎng)實(shí)測(cè)資料才能對(duì)地基的液化情況進(jìn)行判別。

為進(jìn)一步分析波浪荷載作用下海床地基的液化及其影響因素，筆者通過(guò)求解線性波理論和Biot波動(dòng)方程，結(jié)合一定的初始條件和邊界條件，得到海床土體在前進(jìn)波作用下的位移、有效應(yīng)力和孔壓的解析表達(dá)式。利用Jeng提出的海床液化判別準(zhǔn)則，計(jì)算了不同波浪周期、水深、海床土體的物理力學(xué)性質(zhì)參數(shù)下海床的液化深度，分析了這些參數(shù)對(duì)液化的影響，總結(jié)其變化規(guī)律，提出海床抗液化的工程措施。

1 波浪荷載作用下海床的動(dòng)力響應(yīng)

1.1 泥水交界面處的波壓力

當(dāng)波浪在海床上傳播時(shí)，會(huì)在泥水交界面處作用一隨時(shí)間變化的波浪壓力。該壓力會(huì)在海床中產(chǎn)生超靜孔隙水壓力，并以一定的阻尼比和相位滯后的形式傳播到海底。超靜孔隙水壓力的增加和有效應(yīng)力的減小會(huì)導(dǎo)致海床剪切強(qiáng)度下降。當(dāng)超靜孔隙水壓力大于有效應(yīng)力的時(shí)候，可能會(huì)導(dǎo)致海床的失穩(wěn)或者發(fā)生液化。

圖1給出了一簡(jiǎn)諧前進(jìn)波在海床中沿著x方向傳播的示意圖。根據(jù)線性波理論，海床表面(z=0)的波壓力pb為

(1)

根據(jù)線性波理論，重力波的色散方程為

ω2=gktanh(kd)

(2)

式中：H為波高；d為水深；rw為水的重度；g為重力加速度；k=2π/L為波數(shù)，L為波長(zhǎng)；ω=2π/T為頻率，T為波浪周期。

隨著海床表面波壓力pb的變化，海床內(nèi)部的有效應(yīng)力、孔隙水壓力和土骨架的抗剪強(qiáng)度隨之呈現(xiàn)動(dòng)態(tài)變化。

圖1 波浪在海床上傳播示意圖Fig.1 Skeleton of wave propagation on a porous seabed

1.2 海床地基的動(dòng)力響應(yīng)

1.2.1 基本假定

在分析波浪荷載作用下海床地基的動(dòng)力響應(yīng)和液化時(shí)，做如下假設(shè)：

1) 海床是水平、完全飽和、各向同性、有限或無(wú)限厚度的。

2) 海床地基可看作是固相介質(zhì)和均勻分布的牛頓黏性可壓縮流體顆粒的液相介質(zhì)組成。其中孔隙是連通的，土骨架和孔隙流體都是可壓縮的。

3) 土骨架看作是線性、彈性和小變形的，且服從胡克定律。

4) 假設(shè)水波為周期性的簡(jiǎn)諧平面波，沿x方向進(jìn)行傳播，如圖1所示。在分析時(shí)，考慮土-水相互作用的慣性效應(yīng)和波傳播過(guò)程中的能量損失。

5) 海床表面的孔隙水壓力等于海床泥-水分界線處波浪產(chǎn)生的波壓力。

1.2.2 邊界條件

在海床表面，即z=0時(shí)，垂直有效應(yīng)力和剪應(yīng)力非常小，可視作零；海床表面的孔隙水壓力等于泥水交界面的波壓力。因此有

(3)

對(duì)于有限厚度海床，其海床底部視為剛性不透水層，即當(dāng)z=-h時(shí)有

(4)

對(duì)于無(wú)限厚度海床，隨著深度的無(wú)限增加，孔壓與位移逐漸消散，即當(dāng)z→-∞時(shí)有

ux=uz=pf→0

(5)

1.2.3 海床土體動(dòng)力響應(yīng)的解析表達(dá)式

Biot[5]在1956年建立了經(jīng)典的Biot波動(dòng)理論。該理論考慮了慣性效應(yīng)和固體與流體的壓縮性，可用于分析波浪荷載作用下海床的動(dòng)力響應(yīng)?；谇笆黾僭O(shè)和邊界條件式(3～5)，經(jīng)過(guò)繁瑣的代數(shù)運(yùn)算，可得到海床土體應(yīng)力和孔壓的解析表達(dá)式為

(6)

(7)

(8)

(9)

式中：μ，G分別為土體的泊松比和剪切模量；aj，bj，cj，λj分別為與波浪和海床物理力學(xué)性質(zhì)參數(shù)有關(guān)的系數(shù)，其具體表達(dá)式詳見(jiàn)參考文獻(xiàn)[6]。

2 波浪荷載作用下海床的液化

2.1 海床的液化及其判別

在沿海工程中，海床液化機(jī)理主要有兩種：殘余液化和瞬時(shí)液化。其中，殘余液化主要是由于波浪和地震作用使海床土體中產(chǎn)生了累積孔隙水壓力而引起的。根據(jù)有效應(yīng)力原理，累積孔隙水壓力的增加，則有效應(yīng)力相應(yīng)地減小。當(dāng)土骨架間的有效應(yīng)力減小至零的時(shí)候，土壤就出現(xiàn)了液化現(xiàn)象。對(duì)于瞬時(shí)液化，Zen等通過(guò)模型試驗(yàn)和現(xiàn)場(chǎng)觀測(cè)驗(yàn)證了瞬時(shí)液化主要是由于波浪作用對(duì)海床產(chǎn)生的壓力和海底的動(dòng)態(tài)孔隙水壓力之間產(chǎn)生的相位滯后，海底的瞬時(shí)液化在波谷以下的區(qū)域會(huì)呈周期性的出現(xiàn)和消失[2]。

目前，國(guó)內(nèi)外學(xué)者采用土性指標(biāo)來(lái)評(píng)估動(dòng)態(tài)荷載作用下土體液化，這些指標(biāo)包括黏粒含量、液限值、塑性指數(shù)、含水量和側(cè)限應(yīng)力等。但是，這些方法都是基于工程經(jīng)驗(yàn)，說(shuō)明土體有液化的可能，不能說(shuō)明土體一定會(huì)發(fā)生液化。同時(shí)，這些經(jīng)驗(yàn)方法很難用于數(shù)值計(jì)算。對(duì)于土體瞬時(shí)液化的計(jì)算，可采用以下幾種判別方法：

1) Okusa[7]認(rèn)為如果垂直有效應(yīng)力大于土體重力，海床將會(huì)發(fā)生液化，即

(10)

2) Zen等[2]提出了前進(jìn)水波作用下海床液化的判別標(biāo)準(zhǔn)，即

-(γs-γw)z≤pf-pb

(11)

3) Tsai[8]認(rèn)為有效正應(yīng)力等于0時(shí)，土壤將會(huì)發(fā)生液化，即

(12)

Jeng[3]將Zen等[2]的液化判別準(zhǔn)則擴(kuò)展至三維情況，即

(13)

4) 王棟等[9]基于Jeng的理論，在考慮殘余孔隙水壓力的情況下提出了液化判別準(zhǔn)則，即

(14)

5) Lin等[10]用液化參數(shù)LF(孔隙水壓力梯度/土浮重度)來(lái)判別液化的發(fā)生，即

(15)

當(dāng)LF≥1時(shí)，土壤發(fā)生液化。

6) Gao等[11]采用耦合剪切試驗(yàn)和擬合曲線的方法，提出了臨界剪切應(yīng)力比來(lái)判斷海床是否發(fā)生液化，即

(16)

7) Ye[12]考慮土壤的黏聚力和內(nèi)摩擦角，根據(jù)摩爾-庫(kù)侖定律提出了三維情況的液化判別方法，即

(17)

2.2 海床液化判別方法的選擇

3 算例分析與討論

3.1 計(jì)算參數(shù)

從前述理論分析可知：影響波浪荷載作用下海床動(dòng)力響應(yīng)和液化的因素主要包括波浪參數(shù)和海床的物理力學(xué)性質(zhì)指標(biāo)。表1給出了計(jì)算中所用的波浪和海床物理力學(xué)性質(zhì)參數(shù)。

表1 波浪和海床物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of wave and seabed

3.2 波浪對(duì)海床液化深度的影響

根據(jù)線性波理論，影響海床表面波壓力的大小包括波高、水深、波數(shù)以及周期。根據(jù)色散方程式(2)可知：波數(shù)、水深和周期有一定的關(guān)系。下面采用式(13)對(duì)不同水深d和周期T下海床瞬時(shí)液化深度進(jìn)行計(jì)算，其結(jié)果如圖2所示。

圖2 不同水深和波浪周期下的海床液化深度變化情況Fig.2 Variation of liquefaction depth for different water depth and wave period

圖2給出了不同水深d和波浪周期T下海床的液化深度分布情況。從圖2(a)可看出：海床液化深度先是隨著水深的增大而增大；當(dāng)水深大于5.0 m后，則隨著水深的增大而減小。在其他條件都相同的情況下，與水深為2.0 m的最大液化深度值相比，水深分別為5.0，8.0，10.0，15.0 m下的液化深度分別是其值的124%，120%，111%，80%。因此，存在一個(gè)臨界水深。當(dāng)超過(guò)臨界水深后，海床土體的最大液化深度則隨著水深的增大而減小。從圖2(b)可看出：海床最大液化深度隨周期的增大而增大。在一定的水波和海床土體力學(xué)參數(shù)下，周期分別為6.0，7.0，8.0，10.0 s下的液化深度分別是周期為4.0 s時(shí)液化深度的4.31，6.06，7.75，10.94倍。這是由于在其他條件都相同的情況下，周期越大，海床泥水分界線處的波壓力越大。因此，周期T越長(zhǎng)，海床土體的液化深度越大。

3.3 海床土性對(duì)液化深度的影響

海床土性指標(biāo)影響海床的液化深度。為說(shuō)明海床土性指標(biāo)對(duì)液化深度的影響，分別考慮不同滲透率K、泊松比υ、剪切模量G和流體壓縮性kf下的海床液化深度變化情況，其結(jié)果如圖3所示。

圖3(a)是不同滲透率K下海床液化深度的變化情況。隨著K的減小，海床的液化深度逐漸減小。當(dāng)K=10-8m2時(shí)，海床液化深度zL達(dá)到最大。在其他條件相同的情況下，其液化深度約為K=10-12m2時(shí)的1.2倍，這在一定程度上解釋了具有較高滲透性的砂土比低滲透性粉土和黏性土更容易液化。圖3(b)給出了不同泊松比下海床液化深度的變化情況。隨著泊松比的增大，海床液化深度逐漸減小。泊松比的大小反應(yīng)水平有效應(yīng)力與豎向有效應(yīng)力比值的大小。這說(shuō)明水平有效應(yīng)力對(duì)液化深度有一定影響，在計(jì)算中，需要考慮水平應(yīng)力的影響。由圖3(c)可知:隨著剪切模量G的增加，海床的液化深度也逐漸減小，對(duì)比G=5.0×106Pa的結(jié)果，可以發(fā)現(xiàn)最大液化深度zL在G=1.0×107，5.0×107，1.0×108，5.0×108Pa時(shí)分別減小39.54%，34.58%，29.10%和7.34%。圖3(d)是不同流體壓縮性kf下海床液化深度的變化情況。隨著kf的增大，海床液化深度逐漸增大。水中若含有一定的空氣后，流體的壓縮性將急劇增加。當(dāng)kf>2.0×108Pa后，海床液化深度增加很小，說(shuō)明此時(shí)海床液化深度受流體壓縮性的影響很小。

圖3 不同土性指標(biāo)對(duì)海床液化深度的影響Fig.3 Variation of liquefaction depth for different seabed parameters

3.4 海床液化的主要防治措施

從上述數(shù)值算例可知：海床液化深度受水波的周期、水深以及海床土體的剪切模量、滲透系數(shù)、泊松比以及土體中流體(主要指水和空氣)壓縮性的影響。因此，對(duì)于實(shí)際工程，針對(duì)這些影響因素提出如下海床液化防治方法：

1) 根據(jù)海床的臨界水深，對(duì)海上工程所在的一定區(qū)域內(nèi)進(jìn)行海水深度實(shí)時(shí)監(jiān)測(cè)，當(dāng)海水深度小于臨界水深時(shí)，需進(jìn)行海床基底穩(wěn)定性監(jiān)測(cè)。

2) 當(dāng)海床最大液化深度小于不液化的顆粒層厚度時(shí)，海床將不會(huì)發(fā)生液化[13]。較高滲透性的砂土比低滲透性粉土/黏性土更容易引發(fā)液化。因此，可采用滲透性低的黏性土層替換易液化的砂質(zhì)層。

3) 土體剪切模量越大，海床液化深度越小。因此，可通過(guò)對(duì)海床受力層的加固，增加土體的抗剪強(qiáng)度來(lái)防止或減弱海床的液化。

4) 波浪加速度導(dǎo)致海床的液化深度變得窄而淺，海床土體破壞范圍變小[14]。因此，可進(jìn)一步研究造波裝置，反向利用其功能適當(dāng)削減波浪對(duì)海床的泥-水交界面處的壓力，從而消弱或防止海床的液化。

4 結(jié) 論

基于線性水波理論和Biot波動(dòng)理論，獲得了波浪荷載作用下海床土體中孔隙水壓力和有效應(yīng)力的解析表達(dá)式，研究了波浪荷載作用下海床的液化和其主要影響因素，分析了水波的周期、水深和海床土體的強(qiáng)度、滲透率、泊松比、流體的壓縮性等水波和海床土體物理力學(xué)參數(shù)對(duì)海床瞬時(shí)液化的影響，得出不同參數(shù)下海床土體液化的變化規(guī)律，提出了海床液化的防治措施。