楊子希，袁春光，張 弛*，張華慶，歐陽群安

(1.河海大學(xué) 港口海岸與近海工程學(xué)院，南京 210098；2.交通運輸部天津水運工程科學(xué)研究所 港口水工建筑技術(shù)國家工程實驗室 工程泥沙交通行業(yè)重點實驗室，天津 300456)

在海洋環(huán)境中，佇立在海床上的樁柱常受到波浪、水流等動力因素的作用，產(chǎn)生嚴(yán)重的沖刷，據(jù)相關(guān)資料統(tǒng)計，由于海床沖刷導(dǎo)致結(jié)構(gòu)失穩(wěn)破壞占工程事故的24.6%[1]。因此海床的沖刷深度與范圍是在工程設(shè)計中需要考慮的重要因素。

對于單樁的沖刷問題而言，前人已在恒定流或者波浪單獨作用的條件下展開了大量研究，這方面的文獻可以參考Melville[2-3]、Breusers[4]、Ettema[5]、Raudkivi[6-8]、Sumer[9-11]等。研究表明了恒定流作用時，控制沖刷深度ds的主要參數(shù)有相對水深h/D(D為樁柱直徑)，希爾茲參數(shù)θ，相對粒徑d/D，級配σ等；波浪作用時則為Kc數(shù)及相對直徑D/Lw(Lw為波長)。

值得注意的是，上述國內(nèi)外研究成果幾乎集中在沙質(zhì)海床，而對極細沙海床的研究成果開展較少。我國海岸工程中將粒徑處于0.063～0.125 mm的泥沙稱為極細沙，極細沙易起動，易沉降，且在一定動力條件下還可發(fā)生液化現(xiàn)象[20]。依據(jù)Sumer[11]的研究，極細沙底床的樁墩沖刷深度最大可達到沙質(zhì)底床的兩倍。上述結(jié)論是在波浪單獨作用下得出的，參考WangRK、Eadie、李林普、袁春光等針對沙質(zhì)底床的研究成果，在波、流聯(lián)合作用下，其局部沖刷深度也會更大[12-13,21-22]。

本文依托粵電湛江新寮海上風(fēng)電項目，采用正態(tài)系列模型與經(jīng)驗公式相結(jié)合的方法，旨在確定單樁結(jié)構(gòu)在極細沙海床上的沖刷深度及不同防護措施的防護效果，同時驗證經(jīng)典經(jīng)驗公式在極細沙海床下的適用性,為工程設(shè)計提供參考。

1 模型設(shè)計及試驗方法

對于沖刷物理模型來說，常采用動床模型或者系列模型。對于動床模型來說，如果能保證泥沙運動相似以及重力相似，則可以基本消除比尺效應(yīng)，且試驗周期短、成本低。然而在實際操作中，滿足泥沙運動相似十分困難——通常為了滿足泥沙的運動相似，要求試驗?zāi)Ｐ蜕沉捷^細重率較小，而過細的模型沙又會帶來絮凝及黏結(jié)力等一系列問題，反過來又使運動相似難以滿足，尤其是對于擬建工程區(qū)域的極細沙海床，在波流聯(lián)合作用下極其活躍，是常用的模型沙所難以模擬的。同時對于動床模型來說，還需要現(xiàn)場沖刷資料來對實驗結(jié)果進行驗證，由于該工程還在建設(shè)中，無法得到實際的沖刷資料，因此對于本試驗來說，基于比尺縮放模型沙的動床模型很難解決問題。

對于由泥沙運動不相似所帶來的結(jié)果偏差問題，一般可通過做多組幾何比尺不同的試驗并將試驗結(jié)果外延來消除偏差，也就是系列模型方法。系列模型的優(yōu)勢在于可不對泥沙進行縮尺而采用原型沙當(dāng)作試驗沙，并且也無需實測資料為試驗提供驗證數(shù)據(jù)。

1.1 系列模型原理[23-25]

沙玉清假定模型試驗所研究的物理量與影響這一物理量的其他變量之間存在如下指數(shù)關(guān)系

X=Khα1Yα2Zα3

(1)

式中：X為所研究的物理量，可以是沖淤深度，也可以是沖淤時間；h為深度；Y、Z為深度以外的其他水力、泥沙因素變量；α1、α2、α3分別為這些變量的指數(shù)；K為系數(shù)。這里自變量僅取3個，當(dāng)有較多自變量時，也可同樣納入。

對上述方程式寫成比尺關(guān)系式，其結(jié)果為

(2)

沙玉清認為系數(shù)比尺λK是由模型縮小引起的，假定它與深度比尺成指數(shù)關(guān)系

(3)

這樣，上述比尺關(guān)系式將轉(zhuǎn)化為

(4)

式中：α為全部水力、泥沙因素變量比尺均轉(zhuǎn)化為深度比尺后的指數(shù)和。對于原型而言，λh=1，λx=1，故應(yīng)有C=1，最后得

(5)

(6)

(7)

綜上分析可見，當(dāng)采用一定幾何比尺的模型進行試驗時，則λh為已知，Xm可以通過觀測求得，式中的未知量尚有α、Xy兩個。如果能做兩個模型就可建立兩個方程式，聯(lián)解后便可求得Xy，或?qū)⑺肵m及λh的數(shù)據(jù)點繪在雙對數(shù)紙上，連接1、2兩點，與λh=1或lgλh=0的縱軸相交，交點縱坐標(biāo)即為所求原型物理量的數(shù)值Xy。

當(dāng)模型沙采用原型沙或者模型沙與原型沙遵循同一起動規(guī)律時，則泥沙起動流速相似比尺一般式為

(8)

在滿足重力相似條件下流速比尺為

(9)

則得完全滿足相似理論的正態(tài)動床模型幾何比尺為

(10)

將有關(guān)參數(shù)代入上式，即可得出經(jīng)系列模型延伸到比尺λh，即可得出原型情況下的沖刷深度。即當(dāng)采用原型沙時，得出λh為1；當(dāng)采用模型沙時，λh不等于1。

1.2 模型設(shè)計

擬建海上風(fēng)電項目工程場址位于湛江市新寮島東北面海域，水深在6～9 m(海底面高程-7.0～-11.5 m)。項目規(guī)劃裝機總?cè)萘繛?03.5 MW，擬布置37臺5.5 MW單圓樁風(fēng)電機組(直徑為7.5 m)。本次研究的沖刷對象即為該風(fēng)電基礎(chǔ)。

在波流共同作用下近岸地區(qū)的沖刷模型應(yīng)滿足幾何相似，重力相似，水流運動相似，波浪相似和泥沙起動相似等，同時試驗沙與原型沙休止角相似，保證了沖刷坑的幾何形態(tài)相似。根據(jù)該項目物理模型試驗要求，考慮到實驗室條件以及造波和生流能力，模型的比尺設(shè)計為30，40，50。

為了滿足實際需要，試驗在1：30比尺下研究了碎石以及沙被兩種防護措施的效果。碎石防護采用3種方案，碎石直徑依次增大，其中方案一采用d50=2 mm的碎石進行覆蓋，防護層頂整體高出床面5 mm，呈橢圓形，長軸1 m，短軸0.8 m；方案二底層采用d50=2 mm碎石圓形覆蓋，半徑0.5 m，厚度17 mm，表層加蓋厚27 mm，半徑0.4 m，d50=5 mm的碎石；方案三底層采用d50=2.5 mm碎石圓形覆蓋，厚度17 mm，表層加蓋厚27 mm，d50=10 mm的碎石，覆蓋半徑均為0.5 mm。沙被方案整體呈八角形，由四塊沙被搭接而成，整體長、寬均為1.2 m，厚度17 mm，布置詳見2.2節(jié)。

圖1 泥沙級配曲線Fig.1 Size frequency curve of model sand

1.3 模型沙的選擇

1.4 試驗條件與過程

依據(jù)實測水文資料，該海域觀測期最大垂線平均流速為0.93 m/s。在風(fēng)電場設(shè)計過程中考慮在臺風(fēng)影響條件下50 a一遇可能的最大流速為2.14 m/s。風(fēng)電廠址處水深在6～9 m，此處波高衰減很快，并且臺風(fēng)影響期間波浪已發(fā)生破碎。50 a一遇波浪+極端高水位(5.47 m，理基)作用下H13%=5.28 m，Tp=13.1 s；50 a一遇波浪+極端低水位(-0.34 m，理基)作用下H13%=3.22 m，Tp=11.8 s；5 a一遇波浪+極端高水位作用下H13%=4.82 m，Tp=10.9 s。

試驗在交通運輸部天津水運工程科學(xué)研究院展開，試驗水槽長70 m，寬3 m，擬建模型試驗段長4 m，寬2 m。海床地形按照特征實測地形制作，特征試驗基礎(chǔ)布置在模型中部。動床鋪設(shè)深度0.5 m，按照系列試驗最大比尺，相當(dāng)于原型深度15.0 m。為了模擬對單樁基礎(chǔ)刷最不利的情況，試驗采用單向流，水流流向與波浪傳播方向一致。試驗工況如表1所示。

測量儀器采用Vectrino三維點式多普勒流速儀、超聲波水位儀、BG-1型波高傳感器及Z+F三維激光掃描儀，分別對特征點流速、流向、水位、波高地形等進行測量。同時在模型迎波面、背波面和兩側(cè)分別貼有輔助標(biāo)尺，在每次試驗前后分別讀取標(biāo)尺上的床面高度，檢驗三維激光掃描儀的測量精度。

表1 試驗工況Tab.1 Experimental condition

2-a 工況12-b 工況22-c 工況3圖2 不同工況下沖刷坑形態(tài)示意圖(單位:m)Fig.2 Schematic diagram of scour pit morphology under different conditions

2 試驗結(jié)果分析

2.1 無防護下沖刷分析

2.1.1 沖刷坑幾何形態(tài)

圖2為不同動力條件作用后圓柱周圍的沖淤情況。對于工況1、2而言，其波流動力條件較強，沖刷深度比較接近，沖刷坑形態(tài)均接近于水滴形；工況3減弱了波浪與水流強度，沖淤形態(tài)與工況1、2有明顯不同，沖刷坑集中在樁兩側(cè)，而樁前、后沖刷深度與遠處床面接近，無明顯沖刷坑，類似于動床沖刷。這是因為隨著波流強度減弱，馬蹄渦的強度減弱，此時造成沖刷的原因主要是圓柱兩側(cè)的束水沖刷。

圖3 系列模型延伸圖(單位：m)(極端高水位+50 a一遇流速、波浪)Fig.3 Extrapolation map of the scour depth in series models

三組工況下，沖刷均最先出現(xiàn)在圓樁迎波面兩側(cè)夾角與軸線約呈45°處，并且隨著沖刷深度的增加，沖刷坑逐漸向迎波面和背波面發(fā)展。由于極細沙易起動的特點，且實驗中水流強度較大，底床泥沙起動后，直接懸揚并由水流帶向下游，圓樁上下游并未觀察到泥沙落淤。

2.1.2 最大沖刷深度

通過比對5組工況發(fā)現(xiàn)，樁周沖刷深度初期發(fā)展較快，當(dāng)發(fā)展到一定階段后沖刷減慢，逐漸達到極限深度，其中樁前及兩側(cè)沖刷深度接近，背水側(cè)沖刷深度略小。根據(jù)1.1節(jié)方法繪制雙對數(shù)曲線(圖3)，得到“極端高水位+50 a一遇流速+50 a一遇波浪”下單樁基礎(chǔ)最大沖刷深度為15.09 m。

2.2 不同方案防護效果比較

圖4展示了不同方案的防護效果，三組碎石防護在試驗進行到2 h(模型值)后，均出現(xiàn)了不同程度的破壞情況，其中方案一碎石防護層破壞的最為嚴(yán)重。隨著拋石粒徑的增加，方案二、三的破壞程度逐漸減小，但對于方案三而言圓樁周圍沖刷坑也深入底床，表面d50=10 mm粗碎石層的邊緣發(fā)生了明顯的移動破壞，底層黃色細碎石部分裸露。

4-a 碎石防護方案一(碎石粒徑d50=2 mm,厚5 mm)

4-b 碎石防護方案二(表層碎石粒徑d50=5 mm,厚27 mm;底層碎石粒徑d50=2 mm,厚17 mm)

4-c 碎石防護方案三(表層碎石粒徑d50=10 mm,厚27 mm;底層碎石粒徑d50=2.5 mm,厚17 mm)

4-d 沙被方案圖4 防護措施及其比較(單位:m)Fig.4 Protective measures and their comparison

對于沙被防護方案而言，當(dāng)試驗進行了13 h(模型值)，依然起到了較好的防護效果，表現(xiàn)為：(1)單圓樁風(fēng)電基礎(chǔ)周圍底床受到砂被的保護作用未發(fā)現(xiàn)沖刷現(xiàn)象；(2)在波流作用下，沙被各處保持完好，沒有出現(xiàn)“掀起”或者“整體移動”現(xiàn)象；(3)由沙被引起的二次沖刷現(xiàn)象比較微弱，沙被邊緣的二次沖刷深度與床面在波流作用下形成的自然沖刷深度基本一致，沙被造成的二次沖刷與床面的自然沖刷相互融合，無法區(qū)分二次沖刷坑邊緣位置。

3 國內(nèi)外經(jīng)典經(jīng)驗公式計算及同類沖刷試驗結(jié)果對比

3.1 國內(nèi)外經(jīng)典經(jīng)驗公式計算

國內(nèi)外很多學(xué)者對波流共同作用下的大直徑單圓樁基礎(chǔ)沖刷進行了試驗研究，其經(jīng)驗公式總結(jié)如下

(1)Sumer公式[14]

(11)

式中：ds為樁基沖刷深度；dsc為相同水深、流速條件下，水流單獨作用時的平衡沖刷深度；D為直徑；Vc為水流單獨作用下的近底流速(可用距離底床0.5D處的流速表示)；Um為波浪單獨作用下近底水質(zhì)點運動流速最大值；Ucw表示水流導(dǎo)致的流速占波流共同作用形成流速的比例；Lw為波長，該公式僅適用于動床；D/Lw<0.2的細樁。

(2)韓海騫公式[26]

(12)

式中：h為全潮最大水深；B為全潮最大水深下的平均阻水寬度；Fr為弗勞德數(shù)；k1為基礎(chǔ)樁平面布置系數(shù)，條帶型k1=1，梅花型k1=0.862；k2為基礎(chǔ)樁垂直布置系數(shù)，直樁k2=1，斜樁k2=1.176，該式適用于潮流沖刷。

(3)Qi WenGang公式[16]

(13)

該式適用范圍：0.1

(4)陳海鷗公式[27]

(14)

式中：V為斷面平均流速，該式適用于D/Lw>0.2的大直徑樁。

(5)李林普公式[21]

(15)

α=0D/Lw≥0.5；α=-0.102D/Lw<0.5

該式適用于0.07

在30比尺下，5組公式的計算結(jié)果與試驗值的對比見表2。綜合比較上述計算結(jié)果，Sumer公式與Qi公式整體計算結(jié)果與試驗值比較接近。

表2 經(jīng)驗公式計算結(jié)果Tab.2 Calculation results of empirical formula

3.2 計算結(jié)果偏差分析

目前針對單向流作用下的圓樁局部沖刷研究成果較多，但對于近岸波流共同作用下圓樁局部沖刷公式較少，有些研究只考慮到了潮流的作用而忽略了波浪的作用，波浪條件相對于水動力條件較弱時，計算結(jié)果比較接近，遇到波浪作用不可忽視時，計算結(jié)果偏差較大，例如韓海騫公式。

在其余四家具有代表性的公式中，有公式以相對直徑、相對水深等作為主要水力指標(biāo)，也有公式使用Kc數(shù)、相對速率作為指標(biāo)。這些公式雖然考慮了波流共同作用，但它們多由試驗數(shù)據(jù)整理而來，具有一定的適用范圍。對于李林普、陳海鷗公式而言，其試驗以相對直徑D/L>0.15的大直徑單圓樁為背景，而本次試驗圓樁的相對直徑在0.05～0.06，屬于細直徑圓樁，由于大直徑圓樁與細直徑圓樁在沖刷機理上的不同(細直徑圓柱沖刷與樁前下降水流與樁周渦旋強度有關(guān)，大直徑圓柱沖刷與波浪反射有關(guān))，當(dāng)這些公式用來計算細直徑圓柱沖刷時，可能存在較大誤差。

Sumer與Qi的公式以相對速率及Kc數(shù)為主要水力指標(biāo)，計算結(jié)果與試驗值誤差較小，一方面可能是因為試驗設(shè)計的Kc數(shù)在0.4～26，涵蓋了絕大多數(shù)近岸工程中可能出現(xiàn)的情況。另一方面在Qi的試驗中，相對直徑在0.02～0.08，也屬于細直徑圓樁沖刷。Sumer的文章中雖然沒有詳細的列出試驗參數(shù)致使無法確定其相對直徑，但通過在袁春光[17]計算發(fā)現(xiàn)，Sumer公式與Qi公式在預(yù)測大直徑圓柱沖刷時均出現(xiàn)了結(jié)果為0的情況，說明Sumer公式與Qi公式都僅適用于計算細直徑圓樁；李林普與陳海鷗公式在預(yù)測細直徑圓柱沖刷時出現(xiàn)了負值，因此這兩家公式也僅適用于預(yù)測大直徑圓樁沖刷。

對于沙質(zhì)底床，實測資料表明沖刷深度隨粒徑的變細而有增大的趨勢，但無論是哪一家公式均未能正確反映這一現(xiàn)象，尤其是對于極細沙而言，其起動流速位于希爾茲曲線的底部，在波流作用下異常活躍。因此無論是相對直徑、相對水深還是Kc數(shù)、相對速率均無法反映泥沙粒徑和級配的影響。

4 結(jié)論

(1)波流動力條件較強時，極細沙海床在波流作用下的單樁沖刷深度類似于水滴型，隨著動力條件減弱，沖刷坑集中在圓樁兩側(cè)、樁前后沖刷深度接近于動床沖刷。且在三組工況下，沖刷均最先出現(xiàn)在圓樁側(cè)前方，并且隨著沖刷深度的增加，沖刷坑逐漸向迎波面和背波面發(fā)展。

(2)試驗最不利工況為“極端高水位+50 a一遇流速+50 a一遇波浪”(Kc=3.11，相對速率Ucw=0.53)，該條件下單樁基礎(chǔ)極限沖刷深度為15.09 m，達到了風(fēng)電基礎(chǔ)直徑的2倍。

(3)通過比對碎石及沙被兩種防護方案，在試驗設(shè)置的碎石粒徑范圍內(nèi)，碎石防護效果不及沙被防護效果好，圓樁受沙被的保護未發(fā)生沖刷破壞，沙被也未出現(xiàn)“掀起”或“整體移動”現(xiàn)象。

(4)通過對比以往國內(nèi)外學(xué)者提出的經(jīng)驗公式，認為應(yīng)綜合考慮經(jīng)驗公式的適用范圍與實際情況之間的差別，并分析其對沖刷結(jié)果可能帶來的影響趨勢。