曹永港，陳宜展，廖世智，歐陽(yáng)永忠，楊 威，張 燕，肖志建，馬 媛

(1.國(guó)家海洋局南海調(diào)查技術(shù)中心，廣東 廣州 510300；2.自然資源部海洋環(huán)境探測(cè)技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510300)

目前，世界各國(guó)面臨應(yīng)對(duì)氣候變化及環(huán)境持續(xù)惡化的挑戰(zhàn)，均將目光投向最清潔的石化能源—天然氣，紛紛看好其未來(lái)發(fā)展趨勢(shì)。中國(guó)作為世界上最大的能源消費(fèi)大國(guó)和天然氣消費(fèi)增長(zhǎng)最快的國(guó)家，積極參與未來(lái)全球性天然氣市場(chǎng)治理過(guò)程。直灣島LNG項(xiàng)目功能定位是為粵港澳及周邊區(qū)域增加一個(gè)新的可靠的供氣氣源，推動(dòng)清潔能源利用和節(jié)能減排，保障港澳地區(qū)天然氣的供應(yīng)。直灣島LNG接收站碼頭工程位于珠江口外伶仃海域擔(dān)桿列島西部直灣島北側(cè)，擬于東灣順著浪向布置1座裝卸碼頭，于西灣順著岸線依次布置1座裝船碼頭、1座罐箱(兼重件)碼頭、1座加注碼頭和1座工作船碼頭，見(jiàn)圖1。工程區(qū)實(shí)測(cè)資料稀少，工程后產(chǎn)生的水流條件及泥沙回淤情況是值得關(guān)注的問(wèn)題，碼頭泊位淤積問(wèn)題是碼頭設(shè)計(jì)和建設(shè)關(guān)注的重點(diǎn)。泥沙沖淤結(jié)合潮流數(shù)值模擬技術(shù)是目前被廣泛采用的科學(xué)技術(shù)手段之一，用于論證項(xiàng)目選址、工程方案合理性等。戴勇等[1]建立二維潮流數(shù)學(xué)模型模擬分析通用碼頭工程后的流場(chǎng)變化，利用劉家駒公式預(yù)測(cè)工程實(shí)施后泥沙回淤情況；張娜等[2]、張焯等[3]采用潮流泥沙數(shù)值模擬方法研究了LNG碼頭建設(shè)后的水動(dòng)力泥沙問(wèn)題，主要涉及工程實(shí)施后的水流及海床沖淤變化情況。另外，還有順岸式棧橋碼頭工程[4]、廈門雙魚島游艇碼頭工程[5]、舟山岱山島燕窩山碼頭工程[6]、杭州灣開(kāi)敞式長(zhǎng)棧橋碼頭工程[7]等論證研究。

圖1 工程平面布置

為滿足工程的平面布置設(shè)計(jì)需求，基于最新實(shí)測(cè)海流、懸沙、底質(zhì)及地形資料，開(kāi)展潮流數(shù)學(xué)模型試驗(yàn)，預(yù)測(cè)工程實(shí)施后的水流條件，利用經(jīng)驗(yàn)公式進(jìn)行泥沙回淤分析，為工程科學(xué)決策及建設(shè)實(shí)施提供可靠的技術(shù)依據(jù)。

1 水動(dòng)力及泥沙環(huán)境

2019年4月22日—23日在工程海域進(jìn)行16個(gè)站春季大潮水文調(diào)查，調(diào)查內(nèi)容包括海流、懸沙、表層沉積物等。同時(shí)收集附近大萬(wàn)山(21°56′N，113°43′E)、珠海(22°15′N，113°35′E)和赤灣(22°28′N，113°53′E)海洋站2019年4月的氣象與潮位資料。

1.1 潮汐和海流

觀測(cè)海區(qū)的潮汐主要表現(xiàn)為不規(guī)則半日潮特征。大萬(wàn)山站、珠海站和赤灣站最大潮差分別為1.92、2.16、2.45 m。圖2為2019年4月22—23日珠江口垂線平均海流矢量(v垂線=0.3v表+0.5v0.6H+0.2v底)。珠江口外附近海域的實(shí)測(cè)海流主要受珠江沖淡水、南海北部沿岸流和潮汐3種動(dòng)力過(guò)程共同作用。沖淡水自珠江上游河道而來(lái)，在科氏力作用下沿著珠江河口西側(cè)進(jìn)入南海北部；南海北部沿岸流主要影響西側(cè)站位中底層及桂山島東側(cè)站位，主要表現(xiàn)為東向流；潮汐作用為不規(guī)則半日潮特征，由口門外向上游河口方向逐漸增強(qiáng)，漲潮流為北向，落潮流為南向。

圖2 珠江口垂線平均海流矢量圖

1.2 含沙量特征

直灣島離珠江口河口較遠(yuǎn)，幾乎無(wú)內(nèi)河沙源。工程海域海水含沙量低，擬建碼頭附近C1站實(shí)測(cè)含沙量在0.000 1～0.013 5 kg/m3，大潮期平均含沙量為0.003 7 kg/m3，中值粒徑為0.019 mm。

1.3 表層沉積物特征

表層沉積物總體顆粒較細(xì)，中值粒徑范圍為0.009～0.012 mm，類型為細(xì)粉砂。沉積物中粉砂含量最高，平均值為77.30%；黏土含量次之，平均值為17.79%；砂含量第3，平均值為4.92%。

2 二維潮流數(shù)學(xué)模型及泥沙淤積

2.1 潮流數(shù)學(xué)模型

水流連續(xù)方程[8]：

(1)

x、y方向水平動(dòng)量方程分別為：

(2)

(3)

2.2 泥沙淤積計(jì)算

2.2.1正常天氣下回淤

直灣島港址所處海域開(kāi)闊，-20 m等深線距島最近距離約80 m，根據(jù)圖1可知，受直灣島掩護(hù)，西灣碼頭泊位和重件泊位等港池均屬于掩護(hù)式港池；東灣北側(cè)的接卸泊位港池則為開(kāi)敞式，碼頭和港池離開(kāi)陸岸有一定距離，碼頭和棧橋均為透空式高樁梁板結(jié)構(gòu)。這種類型的港口布置幾乎不改變近岸水流的流路，港池回淤計(jì)算公式[9]如下：

(4)

(5)

式中：P1和P2分別為掩護(hù)式和開(kāi)敞式港池泥沙年淤積強(qiáng)度(m/a)；K0為經(jīng)驗(yàn)系數(shù)，取值為0.14～0.17；K1為淤積系數(shù)，取0.35；ω0為黏性細(xì)顆粒泥沙絮凝沉速(m/s)，取0.000 6～0.000 8 m/s；v1、v2分別為工程實(shí)施前后平均流速(m/s)；S為當(dāng)?shù)貙?shí)測(cè)平均含沙量(kg/m3)；T為1 a回淤時(shí)間(s)；ρ0為泥沙干密度；A為港內(nèi)水下淺灘水域面積(m2)；A0為港內(nèi)總水域面積，包括港池和水下淺灘面積(m2)，本工程A與A0的面積取相同。

2.2.2風(fēng)浪天氣下回淤

首先進(jìn)行水體含沙量修正[10]，計(jì)算得到風(fēng)浪情況下的平均含沙量。風(fēng)浪作用下的泥沙可能一部分或大部分以底沙或浮泥形式進(jìn)入港區(qū)、碼頭泊位水域。工程海域表層沉積物平均中值粒徑為0.009～0.012 mm，泥沙輸移采用的浮泥輸沙率計(jì)算公式[11]如下：

(7)

(8)

(9)

式中：qsf為浮泥單寬輸沙率[(kg/(s·m)]；Mc為沖刷系數(shù)[(kg/(m2·s)]，由試驗(yàn)確定；ubmax為波浪的最大軌道速度(m/s)；uc為泥沙起動(dòng)流速；vm為潮流、波流(傳質(zhì)速度)和風(fēng)吹流的合成速度(m/s)；ub為波浪的平均軌道速度(m/s)；φ為波流傳質(zhì)速度方向與航道方向的夾角(°)。

泥沙參數(shù)參考全潮水文觀測(cè)的實(shí)測(cè)結(jié)果，大潮期懸沙中值粒徑平均值為0.019 mm，懸沙含量平均值為0.003 7 kg/m3，表層沉積物中值粒徑平均值為0.009～0.012 mm。

2.3 模型網(wǎng)格與相關(guān)參數(shù)

模型計(jì)算區(qū)域和網(wǎng)格見(jiàn)圖3，模型網(wǎng)格數(shù)為275 620個(gè)、節(jié)點(diǎn)數(shù)為139 540個(gè)。采用三角形非結(jié)構(gòu)化網(wǎng)格，碼頭工程區(qū)域進(jìn)行適當(dāng)加密。根據(jù)樁基實(shí)際位置情況，基于等效過(guò)水面積對(duì)樁基處網(wǎng)格進(jìn)行概化，模型最小網(wǎng)格為1 m。岸邊界采用可滑不可入條件，靠近陸地的水邊界根據(jù)實(shí)測(cè)潮位值得出，外海水邊界利用全球潮汐模型(TPXO9)通過(guò)15個(gè)分潮(M2、S2、N2、K2、K1、O1、P1、Q1、MM、MF、M4、MN4、MS4、2N2、S1)推算得到。糙率根據(jù)水深進(jìn)行插值計(jì)算獲取。

圖3 計(jì)算域及網(wǎng)格剖分

2.4 模型驗(yàn)證

具體驗(yàn)證過(guò)程見(jiàn)文獻(xiàn)[12]，選取2019年4月22日—23日工程海域春季大潮水文調(diào)查其中6個(gè)站進(jìn)行對(duì)比，由圖4、5可知，計(jì)算值與實(shí)測(cè)值吻合較好，滿足精度要求。

圖4 潮位驗(yàn)證

圖5 流速流向驗(yàn)證

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

3.1 潮流場(chǎng)分析

工程方案實(shí)施前后大潮期漲、落急流場(chǎng)見(jiàn)圖6，工程海域漲潮流偏W、落潮流偏E為往復(fù)流。漲潮流經(jīng)直灣島和二洲島之間的二門水道進(jìn)入直灣島北部海域后，向西側(cè)繞過(guò)直灣島東邊岬角折向直灣角工程海域，再向西與繞過(guò)直灣島西邊岬角的海流匯合繼續(xù)向北流動(dòng)。落潮流由西向東流經(jīng)直灣島北部海域，在直灣島東北方分成兩股，一股水流經(jīng)直灣島西邊岬角后繼續(xù)向西流動(dòng)；主流則在直灣島北部海域經(jīng)過(guò)港區(qū)后受二洲島阻擋，一部分流向二洲島北部海域，一部分經(jīng)二門水道流入南海。直灣角西灣海域受岬角掩護(hù)，潮流不大。東部碼頭位于直灣角東側(cè)，靠近直灣角岬角，位于漲落潮流通道，漲潮流向NW、落潮流向SE。工程實(shí)施后，東灣碼頭樁基對(duì)近岸沿岸流有一定的阻水作用，海流流過(guò)樁基后的海流流速有所減弱，樁基之間的流速增強(qiáng)。西灣北碼頭和西灣南碼頭位于西灣內(nèi)，流速較小，工程前后的流速都普遍小于0.10 m/s。西灣南部碼頭與岸線之間海域海水輻聚、流速增強(qiáng)，工程前流速約為0.01 m/s、工程后流速約為0.04 m/s。

圖6 工程方案實(shí)施前后大潮期漲、落急流場(chǎng)

工程建設(shè)后直灣島海域流場(chǎng)變化見(jiàn)圖7。碼頭工程建設(shè)后，管樁起到阻水作用，管樁附近海流流速減小，管樁之間和碼頭外圍水體輻聚、流速增大。東灣海域位于直灣島東側(cè)，靠近水道，海水流速相對(duì)較大，流速變化范圍也較大，主要表現(xiàn)為：外圍管樁繁密海域和近岸海域流速變小，碼頭中部海域海水幅聚、流速增大，流速變大和減小幅度都達(dá)到0.04 m/s。西灣北海域的流場(chǎng)變化表現(xiàn)為近岸管樁附近海域流速明顯減弱，流速減小最大值達(dá)到0.08 m/s。西灣南位于直灣島西側(cè)弱流區(qū)，碼頭建設(shè)后流速相對(duì)變化較小，表現(xiàn)為碼頭區(qū)流速明顯變小，碼頭外側(cè)海域流速小幅增大。從中可比較出，直灣島西側(cè)海域較東側(cè)海域流速小、屬弱流區(qū)，更適合碼頭建設(shè)。

圖7 工程前后平均流速差(工程后-工程前)等值線

3.2 泥沙回淤分析

考慮22個(gè)特征點(diǎn)[13](圖8)，代表了東灣碼頭、西灣北碼頭、西灣南碼頭、直灣島東側(cè)大回旋圓海域和直灣島西側(cè)小回旋圓海域的波浪以及泥沙回淤特征，其水深范圍為5～39 m。工程海域不利波向主要集中在NE與WSW向。5 a、10 a、50 a和100 a一遇有效波高范圍分別為2.8～4.6、3.0～4.9、3.5～5.6、3.8～6.1 m。年回淤計(jì)算考慮波浪各向年平均波高以及相應(yīng)的頻率分布。S向波浪對(duì)碼頭所在地的影響很小，1#、10#～12#點(diǎn)的不利波向?yàn)镹E向，2#～3#點(diǎn)的為ENE向，其余計(jì)算點(diǎn)的為WSW向。驟淤計(jì)算歷時(shí)取12 h，10 a一遇對(duì)應(yīng)各特征點(diǎn)的有效波高見(jiàn)表1。

圖8 特征計(jì)算點(diǎn)

表1 計(jì)算點(diǎn)最大有效波高

直灣島沿岸水深可滿足大型船舶?？?，碼頭工程建設(shè)不需要進(jìn)行清淤，1#～22#特征點(diǎn)不需要開(kāi)挖。工程后泥沙回淤主要考慮碼頭工程建設(shè)后泥沙的自然回淤和驟淤。正常天氣條件下，工程海域水體含沙量較低，約為0.003 7 kg/m3，計(jì)算特征點(diǎn)工程前后的流速變化很小，僅直灣島西南側(cè)的碼頭泊位會(huì)出現(xiàn)回淤，且年平均淤積不超過(guò)0.020 m，其余海域未出現(xiàn)回淤，年平均回淤厚度為0.016 m?？紤]50 a一遇風(fēng)浪后，工程海域水體含沙量有所增加，深水區(qū)域約為0.04～0.17 kg/m3，量級(jí)仍不大，回淤的程度較小。直灣島西南側(cè)的碼頭泊位水深小于10 m，受水流和波浪動(dòng)力的影響較東北側(cè)的開(kāi)敞水域要大，含沙量相對(duì)較高，可出現(xiàn)較大的回淤。計(jì)算大風(fēng)浪作用下的泥沙輸移，各計(jì)算特征點(diǎn)淤積最大的區(qū)域?yàn)?8#，50 a一遇風(fēng)浪的淤積厚度為0.32 m；淤積最小的區(qū)域?yàn)?#，50 a一遇的淤積厚度為0.02 m。見(jiàn)表2。

表2 工程后泥沙回淤預(yù)測(cè)結(jié)果

4 結(jié)語(yǔ)

1)工程海域漲潮流偏W、落潮流偏E，為往復(fù)流。直灣島西側(cè)海域較東側(cè)海域流速小，屬弱流區(qū)，更適合碼頭建設(shè)。工程后并未引起工程海域流場(chǎng)發(fā)生較大變化，僅在碼頭建設(shè)區(qū)域附近有較小程度變化，對(duì)周圍水動(dòng)力環(huán)境的影響較小。

2)正常天氣條件下僅直灣島西南側(cè)的碼頭泊位會(huì)出現(xiàn)回淤，年平均淤積不超過(guò)2 cm，其余海域不會(huì)出現(xiàn)回淤。極端天氣下工程海域短期回淤也較小，50 a一遇風(fēng)浪淤積厚度為0.32 m。

3)工程位置離河口較遠(yuǎn)，泥沙來(lái)源不豐富，水清沙少，泥沙回淤較小，工程后對(duì)周邊海域泥沙環(huán)境影響整體較小。