波浪作用下可滲沙質(zhì)海床模型相似率研究*

孫天霆， 王登婷， 李巖汀， 劉清君， 黃 哲

(南京水利科學(xué)研究院， 港口航道泥沙工程交通行業(yè)重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京210024)

近年來， 石油、 天然氣和金屬礦物等海洋資源的開發(fā)以及海洋空間利用日益引起國內(nèi)外的高度重視， 海洋的開發(fā)利用已成為人類21 世紀(jì)支柱性技術(shù)產(chǎn)業(yè)之一， 從而為海洋工程科學(xué)與技術(shù)的發(fā)展提供了前所未有的機(jī)遇和挑戰(zhàn)。 然而海洋環(huán)境極端惡劣， 在海洋工程設(shè)施建設(shè)與使用中， 海床滲流除造成波浪衰減外， 還會(huì)在海床上建筑物(防波堤和平臺(tái)、墩柱等)基床上產(chǎn)生滲透壓力和滲流力， 影響建筑物的穩(wěn)定性， 嚴(yán)重的還會(huì)引起海床沙土液化， 造成建筑物傾斜和位移， 存在巨大的生命和財(cái)產(chǎn)損失風(fēng)險(xiǎn)。 因此， 研究波浪與可滲海床的相互作用具有重要的工程意義和學(xué)術(shù)價(jià)值。

國外早在20 世紀(jì)40 年代就對波浪周期荷載作用下海床土體動(dòng)力反應(yīng)等一系列海洋土力學(xué)問題展開了研究， 但我國直到20 世紀(jì)70 年代隨著北海油田的開發(fā)才真正引起了人們的重視， 開始對波浪作用下海床的動(dòng)力響應(yīng)進(jìn)行研究。 從70 年代至今， 為了驗(yàn)證各種理論(主要是彈性解)的正確性和適用范圍， 國內(nèi)外學(xué)者在理論分析的同時(shí)也進(jìn)行了大量的模型試驗(yàn)研究。 這些試驗(yàn)大都證明了在原來設(shè)想的波浪作用下海床中存在孔壓幅值隨深度的衰減和孔壓、 有效應(yīng)力的相位遲滯現(xiàn)象， 如Sleath[1]通過波浪與海床相互作用的水槽試驗(yàn)研究， 發(fā)現(xiàn)了海床孔隙水壓力與波形的不同步現(xiàn)象； Tsui 等[2]通過試驗(yàn)進(jìn)一步確認(rèn)了這種相位滯后現(xiàn)象； 高學(xué)平等[3]基于直立堤前海床沖刷形態(tài)相似， 給出了模型沙的比尺取值范圍， 并推導(dǎo)出堤前相對粗沙型沖刷坑的比尺關(guān)系； 王立忠等[4]采用波浪水槽模型試驗(yàn)的方法， 研究波浪荷載作用下沙質(zhì)和粉質(zhì)海床的孔壓響應(yīng)問題， 發(fā)現(xiàn)對于沙質(zhì)海床， 其內(nèi)部超靜孔壓不會(huì)出現(xiàn)累積，而粉質(zhì)海床孔壓累積現(xiàn)象非常明顯， 但土工布可以顯著降低粉質(zhì)海床的超靜孔隙水壓力累積， 起到防止海床液化的作用； 鐘佳玉等[5]采用波流水槽模型試驗(yàn)的方法， 研究了規(guī)則波和不規(guī)則波作用下沙質(zhì)海床的孔隙水壓力響應(yīng)問題， 主要考慮不同深度、 波高及周期對孔隙水壓力的影響；Zhang 等[6]進(jìn)行了規(guī)則波作用下， 波浪與沙質(zhì)海床相互作用的物理模型試驗(yàn)， 探討最大孔隙水壓力隨海床深度、 波周期的變化規(guī)律， 以及孔隙水壓的相位滯后現(xiàn)象和由于海床底部摩阻作用的波高衰減現(xiàn)象。 Zhang 等[7]進(jìn)行了規(guī)則波作用下， 波浪與均勻混合海床相互作用的物理模型試驗(yàn)， 發(fā)現(xiàn)含泥量對海床液化有重要影響， 在高含泥量條件下， 由于海床滲透性降低， 易發(fā)生液化。

但上述波浪物理模型試驗(yàn)只關(guān)注于模型試驗(yàn)本身， 提出的相關(guān)計(jì)算公式及影響因子變化規(guī)律也僅適用于模型， 試驗(yàn)結(jié)果并不能完全應(yīng)用于實(shí)際工程中。 這主要是由于波浪與建筑物相互作用的模型試驗(yàn)系統(tǒng)采用重力相似準(zhǔn)則， 即Froude 數(shù)相似， 并不能簡單適用于波浪與海床相互作用的水-沙系統(tǒng)。 若將海床同樣按重力相似進(jìn)行模擬， 則模型沙尺寸將過小， 如沙質(zhì)海床在實(shí)驗(yàn)室中將被過度縮小為比淤泥還細(xì)的形態(tài)， 二者力學(xué)特性有本質(zhì)不同。

近年來， 土工離心模型試驗(yàn)也逐漸應(yīng)用于波浪與海床相互作用的模擬研究中[8-10]， 其本質(zhì)是將按照相應(yīng)比尺縮小的土工模型置于高速旋轉(zhuǎn)的離心機(jī)中， 讓模型承受大于重力加速度的離心加速度的作用， 以補(bǔ)償由于模型縮尺帶來的土工構(gòu)筑物自重的損失， 使模型與原型土體的應(yīng)力狀態(tài)保持一致[11]， 這也是離心模型試驗(yàn)對土體模擬的最大優(yōu)勢所在。 但離心模型試驗(yàn)也具有一些局限性:離心模型試驗(yàn)中， 顆粒尺寸比尺一般為1∶1， 與其他幾何尺寸的比尺不同， 不滿足相似關(guān)系， 導(dǎo)致模型力學(xué)特性與原型存在一定差異， 即顆粒尺寸效應(yīng)； 擴(kuò)散、 滲流和蠕變時(shí)間的比尺也無法統(tǒng)一，須根據(jù)具體問題進(jìn)行近似取舍。 同時(shí)， 由于離心機(jī)模型箱尺度的局限， 離心設(shè)備無法很好地解決波浪的消能問題， 不能實(shí)現(xiàn)波浪推進(jìn)過程的物理時(shí)程模擬[12]， 其模型施工難度、 試驗(yàn)成本等也遠(yuǎn)高于常規(guī)波浪模型試驗(yàn)。

因此， 在采用重力相似準(zhǔn)則的波浪模型試驗(yàn)系統(tǒng)的基礎(chǔ)上， 提出一種波浪作用下實(shí)驗(yàn)室內(nèi)模擬海床的相似準(zhǔn)則和模擬方法， 對于相關(guān)波浪模型試驗(yàn)研究的開展具有重要意義。

1 模型設(shè)計(jì)

1.1 試驗(yàn)設(shè)備及量測儀器

波浪與海床相互作用的物理模型試驗(yàn)在南京水利科學(xué)研究院河流海岸研究所波浪水槽中進(jìn)行，見圖1。 該水槽長40 m、 寬0.8 m、 深1.0 m。 水槽的一端配有推板式不規(guī)則波造波機(jī)， 可產(chǎn)生規(guī)則波和不同譜型的不規(guī)則波， 水槽兩端均配有消浪緩坡用于吸收波浪。

圖1 試驗(yàn)波浪水槽(單位: m)

1.2 模型沙選取及試驗(yàn)波浪條件

模型沙的選取主要依據(jù)中值粒徑d50來區(qū)分。基于Jensen 等[13]、 Engelund[14]、 Le Mehaute[15]以及Burcharth 等[16]對水力梯度I 的推導(dǎo)， 給出原型與模型的多孔介質(zhì)粒徑比K=dp為原型的多孔介質(zhì)材料的粒徑，dm為模型的多孔介質(zhì)材料的粒徑)與模型幾何比尺λ 的關(guān)系:

式中: α0、 β0為常數(shù)； n 為多孔介質(zhì)的孔隙率；vp為原型多孔介質(zhì)中流體流動(dòng)速度； ν 為水的運(yùn)動(dòng)黏滯系數(shù)。

式(1)表明: 對于完全層流的情況， 多孔介質(zhì)中流體流動(dòng)速度vp較小， ξp→0、 K→λ1∕4； 對于完全紊流的情況， 多孔介質(zhì)中流體流動(dòng)速度vp較大， ξp→∞、 K→λ； 對于完全層流和完全紊流之間的過渡段， 多孔介質(zhì)粒徑比K 介于λ1∕4～λ。

假設(shè)波浪在海床中的傳播為完全層流的情況，則根據(jù)式(1)， 原型與模型的海床粒徑比與幾何比尺的關(guān)系為K==λ海床=λ1∕4。 根據(jù)實(shí)驗(yàn)室條件，選取模型幾何比尺λ=3.94， 則海床比尺λ海床=λ1∕4=1.41， 時(shí)間比尺λt= λ=1.99。

據(jù)此， 選取中值粒徑d50為0.31 mm(原型)和0.22 mm(模型)兩種級(jí)配均勻的模型沙， 并對試驗(yàn)原型值和模型值波浪條件進(jìn)行選取。 試驗(yàn)采用規(guī)則波， 最小入射波高H 為2.0 cm， 最小入射波周期T 為0.8 s， 均符合《波浪模型試驗(yàn)規(guī)程》[17]對原始入射波的規(guī)定， 避免了水的黏滯力和表面張力對試驗(yàn)測量精度的影響。 不同入射波高H、 周期T 及水深h 的試驗(yàn)組合見表1。

表1 試驗(yàn)波要素

1.3 模型及傳感器布置

試驗(yàn)?zāi)Ｐ图皞鞲衅鞑贾靡妶D2。 試驗(yàn)中海床試驗(yàn)段長度為2 m， 為消除試驗(yàn)段前后平臺(tái)對波浪的影響， 在中點(diǎn)處設(shè)置1 個(gè)波高傳感器， 用于采集淺水變形后穩(wěn)定的試驗(yàn)波浪要素； 波高傳感器正下方海床內(nèi)部沿高度方向設(shè)置1 排(4 個(gè))孔隙水壓力傳感器， 用于測量海床孔隙水壓力沿高度方向的分布情況。 試驗(yàn)的不同海床厚度hz以及各孔隙水壓力傳感器距泥面高度h1～h4見表2。

圖2 模型及傳感器布置(單位: m)

表2 海床厚度及各孔隙水壓力傳感器位置分布

2 試驗(yàn)結(jié)果及分析

根據(jù)模型幾何比尺λ=3.94， 選取海床比尺λ海床=λ1∕4=1.41 設(shè)計(jì)的海床孔隙水壓力的試驗(yàn)結(jié)果見圖3。 可以看出， 海床孔隙水壓力沿高度方向基本呈衰減的趨勢， 高度h1處傳感器由于距海床面較淺， 此處孔隙水壓力原型值(dp=0.31 mm)與模型換算值(dm=0.22 mm)受海床影響程度相對高度h2～h4處較小， 因此該處孔隙水壓力值與海床面波浪壓力理論值較接近， 吻合程度相對較好；高度h2和h3處傳感器在海床內(nèi)部， 受不同組次波高、 波長等因素影響， 海床面處波浪動(dòng)水壓力不同， 導(dǎo)致海床內(nèi)部孔隙水壓力值也有較大差異，數(shù)據(jù)點(diǎn)基本呈條狀帶分布； 高度h4處傳感器距海床面最深， 孔隙水壓力受海床多孔介質(zhì)的阻力作用影響充分， 導(dǎo)致不同組次該處孔隙水壓力值較接近， 數(shù)據(jù)點(diǎn)較密集。 但總體而言， 試驗(yàn)結(jié)果孔隙水壓力原型值(dp= 0.31 mm)與模型換算值(dm=0.22 mm)吻合程度一般， 除高度h1處外，高度h2～h4處模型換算值均比原型值偏大， 高度h1處孔隙水壓力模型換算值總體偏大約2.1%， 高度h2處孔隙水壓力模型換算值總體偏大約27.9%， 高度h3處孔隙水壓力模型換算值總體偏大約33.6%， 高度h4處孔隙水壓力模型換算值總體偏大約41.0%。

圖3 按模型幾何比尺與海床比尺關(guān)系為λ海床=λ1∕4的分層試驗(yàn)結(jié)果

因此， 按模型幾何比尺與海床比尺關(guān)系λ海床=λ1∕4設(shè)計(jì)的波浪與海床相互作用試驗(yàn)存在不合理之處。 這主要由于根據(jù)理論推導(dǎo)， λ海床=λ1∕4僅適用于完全層流的情況， 如地下水通過致密多孔介質(zhì)的滲流運(yùn)動(dòng)， 與達(dá)西定律適用的情況基本一致。但波浪作用下海床內(nèi)部的孔隙流動(dòng)雖然由層流運(yùn)動(dòng)占主導(dǎo)， 但仍介于完全層流和完全紊流之間的過渡區(qū)， 是一種“強(qiáng)層流、 弱紊流” 狀態(tài)， 其海床比尺λ海床范圍為λ1∕4～λ。

為此， 基于原型孔隙水壓力與模型值相似的條件， 提出一種實(shí)驗(yàn)室內(nèi)模擬沙質(zhì)海床的相似準(zhǔn)則和模擬方法， 仍認(rèn)為波浪作用下海床內(nèi)部的孔隙流動(dòng)由層流運(yùn)動(dòng)占主導(dǎo)， 但同時(shí)具備一定紊流的特征， 假設(shè)原型與模型的海床粒徑比K 與幾何比尺λ 的關(guān)系為K=λ海床=λ1∕3。 根據(jù)實(shí)驗(yàn)室條件， 選取模型幾何比尺λ = 2.80， 則海床比尺λ海床=λ1∕3=1.41， 時(shí)間比尺λt= λ =1.67， 選取中值粒徑d50為0.31 和0.22 mm 兩種級(jí)配均勻的模型沙。 按該模擬方法修改模型比尺關(guān)系后， 試驗(yàn)不同入射波高H、 周期T 及水深h、 海床厚度hz以及各孔隙水壓力傳感器距泥面高度h1～h4的組合見表3。

表3 比尺關(guān)系修改后試驗(yàn)波要素、 海床厚度及各孔隙水壓力傳感器位置分布

根據(jù)模型幾何比尺λ=2.80， 選取海床比尺λ海床=λ1∕3=1.41 設(shè)計(jì)的海床孔隙水壓力的試驗(yàn)結(jié)果見圖4。 可以看出， 海床孔隙水壓力沿高度方向基本呈衰減的趨勢， 高度h1處傳感器距海床面較淺， 受海床多孔介質(zhì)的阻力作用影響較小， 孔隙水壓力原型值(dp=0.31 mm)與模型換算值(dm=0.22 mm)與海床面波浪壓力理論值較接近， 吻合程度最好； 高度h2和h3處傳感器在海床內(nèi)部， 不同組次孔隙水壓力值有較明顯差異， 數(shù)據(jù)點(diǎn)呈條狀帶分布； 高度h4處傳感器距海床面最深， 孔隙水壓力受海床多孔介質(zhì)的阻力作用影響充分， 導(dǎo)致不同組次該處孔隙水壓力值較接近， 數(shù)據(jù)點(diǎn)較密集。 總體而言， 試驗(yàn)結(jié)果孔隙水壓力原型值(dp=0.31 mm)與模型換算值(dm= 0.22 mm)吻合程度較好， 相關(guān)系數(shù)R=0.96， 高度h1～h4處原型值與模型換算值的平均誤差為1.41%(圖5)。

圖4 按模型幾何比尺與海床比尺關(guān)系為λ海床=λ1∕3的分層試驗(yàn)結(jié)果

圖5 按模型幾何比尺與海床比尺關(guān)系為λ海床=λ1∕3的試驗(yàn)結(jié)果

因此， 按模型幾何比尺與海床比尺關(guān)系λ海床=λ1∕3設(shè)計(jì)的波浪與海床相互作用試驗(yàn)可較好地模擬波浪與海床相互作用時(shí)海床內(nèi)部孔隙流動(dòng)的“強(qiáng)層流、 弱紊流” 狀態(tài)。 該模擬方法可較好地用于可滲沙質(zhì)海床的實(shí)驗(yàn)室模擬。

3 結(jié)論

1)若假設(shè)波浪在海床中的傳播為完全層流的情況， 按原型與模型的海床粒徑比與幾何比尺關(guān)系K=λ海床=λ1∕4選取海床比尺， 海床孔隙水壓力原型值與模型值吻合程度較差， 模型換算值比原型值平均大34.1%。 這是由于波浪作用下海床內(nèi)部的孔隙流動(dòng)雖然由層流運(yùn)動(dòng)占主導(dǎo)， 但仍介于完全層流和完全紊流之間的過渡區(qū)， 是一種“強(qiáng)層流、 弱紊流” 狀態(tài)， 其海床比尺λ海床范圍為λ1∕4～λ。

2)基于原型孔隙水壓力與模型值相似的條件，提出一種新的實(shí)驗(yàn)室內(nèi)模擬沙質(zhì)海床的相似準(zhǔn)則和模擬方法， 仍認(rèn)為波浪作用下海床內(nèi)部的孔隙流動(dòng)由層流運(yùn)動(dòng)占主導(dǎo)， 但同時(shí)具備一定紊流的特征， 按模型幾何比尺與海床比尺關(guān)系λ海床=λ1∕3選取海床比尺可較好地模擬波浪與海床相互作用時(shí)海床內(nèi)部孔隙流動(dòng)的“強(qiáng)層流、 弱紊流” 狀態(tài)?？紫端畨毫υ椭蹬c模型值吻合程度較好， 相關(guān)系數(shù)R=0.96。

3)本文提出的實(shí)驗(yàn)室內(nèi)模擬海床的相似準(zhǔn)則和模擬方法可較好地用于波浪作用下海床孔隙水壓力原型值與模型值之間的換算， 解決水-沙系統(tǒng)的比尺效應(yīng)問題， 對于相關(guān)波浪模型試驗(yàn)研究的開展具有重要意義， 并可應(yīng)用于實(shí)際工程。