張亞敬，劉海源，劉鳴洋，陳漢寶

(交通運輸部天津水運工程科學(xué)研究所，港口水工建筑技術(shù)國家工程研究中心，工程泥沙交通行業(yè)重點實驗室，天津300456)

隨著港口工程大型化以及深水岸線資源的缺乏，離岸深水碼頭建設(shè)越來越多。重力式碼頭要求地基基礎(chǔ)好，同時深水海域浪大、流急導(dǎo)致施工困難，結(jié)構(gòu)斷面龐大且造價高，因此墩柱式透空結(jié)構(gòu)或樁基結(jié)構(gòu)碼頭成為相對較好的結(jié)構(gòu)形式。由于墩柱式透空結(jié)構(gòu)碼頭的沉箱墩并不是孤立存在，波浪傳播至墩群中時存在墩間反射、繞射，使得群墩之間的波高分布呈現(xiàn)出不同于單墩周圍的波高分布規(guī)律，從而影響碼頭面高程設(shè)計以及沉箱墩受力計算。我國現(xiàn)行規(guī)范對透空式碼頭面高程的設(shè)計并未考慮群墩效應(yīng)，定義群墩中波浪力衰減程度的群墩系數(shù)也僅適用于圓柱形墩，并不能反映矩形墩群中的波浪力變化規(guī)律。

關(guān)于群墩波浪力分布規(guī)律，已有許多研究者開展了相關(guān)研究。季新然[1]通過物理模型試驗分析多向不規(guī)則波與群墩作用時波浪在大尺度圓形墩柱上的正向力、橫向力和爬高，同時建立多向不規(guī)則波浪與大尺度墩柱及群墩結(jié)構(gòu)作用的計算方法，給出不同布置形式的群墩系數(shù)，發(fā)現(xiàn)了不規(guī)則波作用下的波浪俘獲現(xiàn)象；許佳[2]建立數(shù)值計算方法，對多向不規(guī)則波與大尺度群墩結(jié)構(gòu)的作用進(jìn)行模擬計算，分析不同布置形式群墩受橫向力和正向力的變化規(guī)律；任效忠等[3]通過物理模型試驗研究不規(guī)則波對準(zhǔn)橢圓沉箱群墩結(jié)構(gòu)的作用力，分析入射波向、周期、相對間距等因素與沉箱所受波浪力之間的關(guān)系，之后又進(jìn)一步開展模型試驗[4]研究規(guī)則波對準(zhǔn)橢圓沉箱群墩結(jié)構(gòu)的波浪力，發(fā)現(xiàn)與正向波浪力相比，群墩效應(yīng)對橫向波浪力的影響更大，對于大型結(jié)構(gòu)準(zhǔn)橢圓沉箱，波浪非線性的影響不容忽視；王永學(xué)[5]提出求解任意截面群墩結(jié)構(gòu)波浪荷載的數(shù)值計算方法，對2個方形等截面群墩的系數(shù)進(jìn)行研究，探討群墩截面幾何形狀不同與群墩間相互影響的關(guān)系；周枝榮等[6]通過物理模型試驗，研究矩形重力墩及墩間面板在斜向波浪作用下的波浪力大小及其群墩折減情況，提出矩形重力墩結(jié)構(gòu)上最大總水平波浪力簡化計算公式?？梢钥闯?，已有研究絕大多數(shù)揭示的是圓柱形群墩受力規(guī)律，針對矩形群墩的受力規(guī)律研究則遠(yuǎn)不及對圓柱形墩的研究豐富。因此，需要開展模型試驗研究，觀察和揭示矩形群墩式碼頭的墩間波高分布規(guī)律以及波浪力在墩間的衰減規(guī)律。

本文以茂名某墩式散貨碼頭工程為例，通過開展三維物理模型試驗，研究波浪作用下矩形沉箱墩群周圍的波高分布及沿程衰減規(guī)律，研究代表性沉箱墩所受水平、側(cè)向和浮托總力以及傾覆力矩，研究波浪力沿碼頭軸線方向的衰減規(guī)律。

1 工程概況

本工程位于茂名港博賀新港區(qū)大宗干散貨作業(yè)區(qū)中突堤東側(cè)岸線，中突堤西側(cè)岸線為已建的粵電項目配套碼頭工程。從南向北連續(xù)建設(shè)1個20萬噸級(M7)、1個10萬噸級(M8)和1個5萬噸級散貨泊位(M9)，水工結(jié)構(gòu)分別按照30萬、20萬和10萬噸級預(yù)留。擬建工程碼頭結(jié)構(gòu)總長度960 m，其中20萬噸級泊位長度447 m，10萬噸級泊位長度276 m，5萬噸級泊位長度237 m，棧橋結(jié)構(gòu)長度401.7 m，均采用透空式墩臺結(jié)構(gòu)，墩臺之間采用大跨度預(yù)應(yīng)力梁連接。距離碼頭前沿線102.4 m處為已建的粵電防波堤。工程位置如圖1所示，工程區(qū)域平面布置見圖2。

圖1 工程位置

圖2 工程平面布置(單位：m)

圖3 沉箱墩布置及波高傳感器分布

圖4 代表性沉箱墩(M7泊位)結(jié)構(gòu)斷面(單位：m)

表1 率定點波要素

2 試驗方法與試驗工況

2.1 模型制作

模型試驗在交通運輸部天津水運工程科學(xué)研究院綜合試驗廳中完成。模型按重力相似準(zhǔn)則設(shè)計，采用正態(tài)、定床模型。試驗場地為55 m×42.5 m×1 m(長×寬×高)的矩形水池(試驗水池邊界見圖2)，采用幾何比尺λL為55。模型中粵電防波堤、東西防波堤和本工程沉箱墩的護(hù)底以及沉箱墩與粵電防波堤之間的護(hù)坡結(jié)構(gòu)中的各種塊石按重力比尺挑選，粒徑級配符合《防波堤與護(hù)岸設(shè)計規(guī)范》[8]以及《波浪模型試驗規(guī)程》[9]的要求。本工程港池內(nèi)的水底地形包括航道、碼頭前沿以及天然地形，按照試驗規(guī)程采用水準(zhǔn)儀進(jìn)行高程控制，偏差控制在±1 mm以內(nèi)；工程區(qū)域所有結(jié)構(gòu)的平面尺寸及位置偏差控制在±1 cm以內(nèi)。矩形沉箱墩、預(yù)應(yīng)力箱梁和碼頭面板等根據(jù)設(shè)計圖紙采用木板組裝制作，并根據(jù)試驗規(guī)程和設(shè)計水位在沉箱墩的側(cè)立面和底板預(yù)留安裝壓強傳感器的孔位。

2.2 試驗設(shè)備及方法

造波采用可移動搖板式造波機及其控制系統(tǒng)。模型造波時，由計算機根據(jù)輸入的造波參數(shù)計算出目標(biāo)波浪的板前波浪信號，并按一定算法將其轉(zhuǎn)換成相當(dāng)于造波板運動速度和位置的數(shù)據(jù)，由伺服驅(qū)動器輸出脈沖信號控制伺服電機的轉(zhuǎn)速和轉(zhuǎn)動的角度，從而造出期望的波浪。本模型試驗研究采用不規(guī)則波進(jìn)行，不規(guī)則波根據(jù)《港口與航道水文規(guī)范》[10]采用JONSWAP譜進(jìn)行頻譜模擬，其解析式為：

Sη(f)=

(1)

(1.094-0.019 15lnγ)

(2)

Tp=Ts/[1-0.132(γ+0.2)-0.559]或

(3)

(4)

本試驗水位采用表1的設(shè)計水位。根據(jù)試驗要求，各個試驗水位依據(jù)給定的波浪要素進(jìn)行率定，率定誤差控制在0%～5%范圍內(nèi)，產(chǎn)生嚴(yán)格控制的正向偏差以保證實際造波波要素偏差的一致性。每組波要素的波列都保持波個數(shù)在100以上，將最后得到的造波參數(shù)存儲在計算機中。正式試驗前，均應(yīng)用小波持續(xù)作用一段時間。模型試驗采用TK-2008型動態(tài)水位測量系統(tǒng)(儀器精度為1.0 mm)對波高進(jìn)行采集分析。波高傳感器布置在每個沉箱的前沿、后沿以及之間(圖3)，以記錄波浪在沉箱之間的分布情況和沿程衰減規(guī)律。

各個矩形沉箱墩所受前后向、側(cè)向和浮托總力采用兩種測力方式取得，其一為利用2010型總力儀直接測量總力；其二為利用布置在沉箱側(cè)立面和底板的若干壓強傳感器，測得壓強值與其代表面積進(jìn)行積分得到各個方向的總力?？偭x測力原理為：在兩個平整的剛性矩形面板之間的4個角上分別布置4個三維測力天平，其中一個矩形面板與所測結(jié)構(gòu)剛性固接，另一個矩形面板與剛性足夠大的固定支架固接，所測結(jié)構(gòu)形心在平面上的投影與傳感器矩形面板的形心重合；結(jié)構(gòu)所受的力傳遞到與其固接的總力儀時，每個測力天平可測得其所受的3個分力Fxn、Fyn和Fzn；根據(jù)六分力原理，由所有測力天平上的共計12個分力可計算得出結(jié)構(gòu)整體所受的6個分力，即水平總力Fx、側(cè)向總力Fy、浮托總力Fz、x軸力矩Mx、y軸力矩My和z軸力矩Mz。各個力的關(guān)系見式(5)，力的分布見圖5。

圖5 總力儀原理

(5)

式中：Lx和Ly分別為每個力傳感器之間在x和y方向上的距離。

對于由壓強值積分得到總力的方式，依據(jù)《波浪模型試驗規(guī)程》和試驗技術(shù)要求，在沉箱結(jié)構(gòu)的底面和側(cè)面布置相應(yīng)數(shù)量的壓強傳感器，其中沉箱上部靜水面附近布置的傳感器密度比下部大，綜合考慮不同設(shè)計水位，并保證能夠準(zhǔn)確捕捉到波峰波谷變化導(dǎo)致的沉箱壁面所受波壓強的變化。波壓力數(shù)據(jù)通過SG2008型數(shù)據(jù)采集系統(tǒng)進(jìn)行采集分析，試驗中連續(xù)采集至少100個波作用的波壓力過程，模型采樣的時間間隔均為0.01 s。試驗時在靜水條件下，對所有測點標(biāo)零，在靜水面以下的測點以此時的靜水壓強作為對應(yīng)測點的零點，在靜水面以上的測點以此時的大氣壓強作為零點。試驗采集到的壓強值為測點實際壓強與標(biāo)零時測點對應(yīng)壓強的差值。將沉箱置于x、y、z坐標(biāo)軸所構(gòu)成的三維直角坐標(biāo)系內(nèi)，然后將各測點所代表的壓強和面積在該坐標(biāo)系x、y、z3個方向分別進(jìn)行投影，再由各測點測得壓強過程線，利用積分得到在x、y、z3個方向所受到的波浪力，最后統(tǒng)計其所受總的水平力Fx、側(cè)向力Fy和浮托力Fz。對于力的方向規(guī)定為：Fx與x軸正方向相同為正，反之為負(fù)，F(xiàn)y和Fz的規(guī)定與Fx相同。將沉箱沿前后和兩側(cè)方向的傾覆力矩的旋轉(zhuǎn)中心定義為沉箱底板兩個方向的中心線，則結(jié)合每個測點受力所對應(yīng)的力臂，可以得到沉箱在前后以及兩側(cè)的總傾覆力矩，即My和Mx。壓強傳感器的布置方式如圖6所示，圖中箭頭所指方向分別定義了x、y、z3個方向的總力以及沿前后和兩側(cè)傾覆力矩的正方向。同時，將總力儀得到的力和力矩與圖6所示坐標(biāo)系相統(tǒng)一。

表2 試驗工況

圖6 壓強傳感器布置方式以及本文定義坐標(biāo)系

2.3 試驗工況

本文針對該工程實例的不同泊位的沉箱間波高分布以及代表性沉箱墩所受的波浪力變化規(guī)律進(jìn)行研究。觀察在極端高水位、設(shè)計高水位和設(shè)計低水位，重現(xiàn)期為50和10 a的波浪條件下，沉箱間波浪傳播現(xiàn)象和波高沿程衰減規(guī)律，揭示各泊位代表性矩形沉箱墩所受波浪力的沿程變化規(guī)律。對于沉箱間波高分布，本模型試驗的沉箱墩上均不放置預(yù)應(yīng)力箱梁，以排除預(yù)應(yīng)力箱梁與水體接觸對沉箱間波高產(chǎn)生的干擾。對于波浪力，沿碼頭軸線方向選取M7泊位的系纜墩沉箱、第1個單沉箱、第1個雙沉箱南側(cè)沉箱和第1個雙沉箱北側(cè)沉箱，選取M8泊位第1個雙沉箱南側(cè)沉箱，選取M9泊位第1個單沉箱，以進(jìn)行墩群中波浪力沿程變化規(guī)律的研究。將上述代表性沉箱按順序進(jìn)行編號1#～6#，其位置如圖3中虛線框和標(biāo)記所示，試驗工況和組次見表2。

3 試驗結(jié)果與分析

3.1 沉箱間波高分布

試驗中對港區(qū)的波浪傳播情況進(jìn)行觀察，發(fā)現(xiàn)港池和航道存在挖深的情況下，加上粵電防波堤的影響，波浪經(jīng)過率定點后傳至本工程區(qū)域時，其傳播方向基本與碼頭軸線方向一致。通過波高傳感器所采集的數(shù)據(jù)發(fā)現(xiàn)，不同水位不同重現(xiàn)期波浪條件下，該碼頭沉箱墩區(qū)域波高沿程分布呈現(xiàn)出相似的規(guī)律，即：受沉箱側(cè)面的阻水和反射的影響，沿碼頭軸線方向，位于M7泊位南端5個沉箱墩范圍內(nèi)波高衰減幅度最大。以最不利的極端高水位重現(xiàn)期50 a的波浪作用為例，M7泊位1#墩南側(cè)測點H1%=7.38 m，隨后沿碼頭軸線方向逐漸衰減，至M7泊位第6個沉箱時H1%=3.07 m，波高衰減58.4%；至M8泊位的5#墩時H1%=1.87 m，波高衰減74.7%；至M9泊位所有沉箱H1%均在1.40 m左右，波高衰減81.0%。極端高水位重現(xiàn)期50 a波浪條件下代表性沉箱四周H1%分布沿程變化如圖7所示，其中圖7b)為歸一化的波高值，其定義為該沉箱的前沿、后沿、沉箱間前側(cè)或沉箱間后側(cè)位置的波高與M7泊位1#墩四周相應(yīng)位置波高的比值。對波高沿程衰減幅度最大的M7泊位南端5個沉箱區(qū)域的波高分布進(jìn)行觀察可知：受相鄰沉箱側(cè)壁的入射波與反射波疊加的影響，波高分布呈現(xiàn)兩沉箱間波高>沉箱前沿波高>沉箱后沿波高的規(guī)律；而隨著波高由南向北沿程衰減，其余位置各沉箱四周波高差別也逐漸減小。

圖7 極端高水位重現(xiàn)期50 a波浪條件下代表性沉箱H1%波高沿程變化

3.2 沉箱受力

根據(jù)表2所列試驗工況，對選取的不同泊位的沉箱進(jìn)行受力測量，結(jié)果見表3、4。其中表3為首先受到波浪作用的M7泊位系纜墩(1#墩)的受力結(jié)果，可以看出，由總力儀直接測量得到的沉箱在各個方向所受波浪力和力矩，與由壓強傳感器所測壓強值積分所得結(jié)果較為接近，且呈現(xiàn)出的隨設(shè)計水位和波浪重現(xiàn)期的變化規(guī)律是一致的，即沉箱所受的水平力Fx、側(cè)向力Fy、浮托力Fz以及側(cè)向傾覆力矩Mx和前后傾覆力矩My均與設(shè)計水位和波高呈正相關(guān)。以上結(jié)果也驗證了本文采用的模型試驗方法以及試驗過程。

表3 M7泊位1#墩受力結(jié)果

表4為所有選取的代表性沉箱(1#～6#墩)在極端高水位和重現(xiàn)期50 a的波浪條件下，所受波浪力的對比。極端高水位重現(xiàn)期50 a波浪條件下代表性沉箱受力沿程變化如圖8所示，其中縱坐標(biāo)為歸一化的各沉箱受力值，其定義為該沉箱在某個方向的總力(Fx、Fy或Fz)或傾覆力矩(Mx或My)與M7泊位1#墩所受相應(yīng)的力或力矩的比值?？梢钥闯觯煌料渌懿ɡ肆Φ难爻趟p幅度，呈現(xiàn)出與沉箱周圍波高沿程衰減情況一致的規(guī)律，即：作用在M7泊位3#和4#墩的波浪力相比M7泊位1#墩所受波浪力已經(jīng)發(fā)生了很大程度的衰減；從M7泊位3#和4#墩開始沿碼頭軸線方向向北，直到M8和M9泊位，沉箱所受波浪力不再有明顯衰減。

表4 極端高水位下重現(xiàn)期50 a波浪的不同泊位代表性沉箱墩受力結(jié)果

圖8 極端高水位重現(xiàn)期50 a波浪條件下代表性沉箱受力沿程變化

3.3 波浪力衰減規(guī)律

本文的矩形沉箱墩形成了串列式的墩群，波浪在矩形沉箱墩組成的墩群中的衰減導(dǎo)致的沉箱所受波浪力的衰減規(guī)律，尚未有試驗進(jìn)行研究。本文針對M7泊位1#～3#墩的串列結(jié)構(gòu)進(jìn)行群墩系數(shù)計算。由于波浪入射方向與碼頭軸線走向非常接近，故本文研究y方向的波浪側(cè)向總力Fy的變化規(guī)律。根據(jù)《港口與航道水文規(guī)范》中群墩系數(shù)的定義和計算方法，對以上選取的3個沉箱墩結(jié)構(gòu)所受y方向的波浪總力Fy的沿程衰減情況進(jìn)行計算，并與模型試驗所得結(jié)果進(jìn)行對比?！陡劭谂c航道水文規(guī)范》定義的參數(shù)k(由譜峰周期Tp對應(yīng)的波長Lp定義的波數(shù)，k=2π/Lp，得0.03)和D(圓柱體直徑，本文定義為矩形沉箱上垂直于碼頭軸線方向的尺度，即矩形沉箱沿碼頭前后方向的長度，取27.55 m)，對于本文矩形沉箱的kD=0.83。則根據(jù)規(guī)范計算所得群墩系數(shù)為：M7泊位2#墩(距離1#墩45 m)，群墩系數(shù)為0.99；M7泊位3#墩(距離1#墩90 m)，群墩系數(shù)為0.96。然而，本文模型試驗得到的對應(yīng)沉箱在y方向所受側(cè)向力Fy的歸一化結(jié)果分別為0.75和0.49，與規(guī)范計算所得結(jié)果相差較大。這是因為矩形沉箱墩存在較大面積的直立壁面，一方面阻水面積較大，另一方面波浪在沉箱間多次反射，致使波能衰減的幅度相比在圓柱形墩群中明顯更大。因此，針對矩形沉箱組成結(jié)構(gòu)的群墩系數(shù)計算，《港口與航道水文規(guī)范》的圓柱形墩的群墩系數(shù)計算方法不再適用。在后續(xù)工作中應(yīng)開展系統(tǒng)性的試驗研究，以得到適用于矩形沉箱的群墩系數(shù)計算方法。

4 結(jié)論

1)受沉箱側(cè)面的阻水和反射的影響，沿碼頭軸線方向，位于M7泊位南端5個沉箱墩范圍內(nèi)波高衰減幅度最大；受相鄰沉箱側(cè)壁的入射波與反射波疊加的影響，波高分布呈現(xiàn)兩沉箱間波高>沉箱前沿波高>沉箱后沿波高的規(guī)律，隨著波高由南向北沿程衰減，其余位置各沉箱四周波高差別逐漸減小。

2)不同沉箱所受波浪力的沿程衰減幅度，呈現(xiàn)出與沉箱周圍波高沿程衰減情況一致的規(guī)律。作用在M7泊位3#和4#墩的波浪力相比M7泊位1#墩所受波浪力已經(jīng)發(fā)生了很大程度的衰減；從M7泊位3#和4#墩開始沿碼頭軸線方向向北，直到M8和M9泊位，沉箱所受波浪力不再有明顯衰減。

3)矩形沉箱墩阻水面積較大，且波浪在沉箱間多次反射，致使波能衰減的幅度相比在圓柱形墩群中明顯更大。針對矩形沉箱組成結(jié)構(gòu)的群墩系數(shù)計算，規(guī)范中的圓柱形墩的群墩系數(shù)計算方法不再適用。