(中交上航局航道建設(shè)有限公司，浙江 寧波 315200)

伴隨著國家“一帶一路”、長江經(jīng)濟(jì)帶等戰(zhàn)略的提出和各項(xiàng)舉措的落地，浙江海洋經(jīng)濟(jì)建設(shè)頻傳利好，特別是寧波梅山港，具有得天獨(dú)厚的優(yōu)良水深和庇護(hù)條件，是不可多得的自然深水良港。寧波—舟山港梅山港區(qū)工程在建設(shè)過程中遇到新建圍堤工程預(yù)留龍口方案的比選難題，最后通過引進(jìn)水交換理論比選出了最佳方案，保障了工程的順利實(shí)施。

寧波—舟山港梅山港區(qū)工程位于寧波梅山港區(qū)1～5號(hào)集裝箱碼頭東側(cè)，毗鄰5號(hào)泊位，北側(cè)緊鄰鹽田大道。該工程主要包括圍堤、吹填砂及軟基處理等項(xiàng)目，其中軟基處理工程中，陸域縱深約為670～780m，陸域?qū)挾燃s為2150m，總面積約166.41萬m2，將通過陸域形成與地基加固后建設(shè)成為集裝箱堆場。工程地理位置如圖1所示。

圖1 工程地理位置圖

1 圍堤工程施工工藝

該工程新建圍堤施工分為袋裝砂基礎(chǔ)、水上、陸拋石、圍堤護(hù)面、防浪墻及堤頂?shù)缆返?，?個(gè)階段實(shí)施：

1.1 一階段圍堤基礎(chǔ)

在一階段圍堤最底部2層通長砂被鋪設(shè)完成后，插排水板進(jìn)行軟基處理，再進(jìn)行第3層通長砂被的施工。在通長砂被施工完畢后立即進(jìn)行圍堤底部軟體排，袋裝碎石墊層的施工；在軟體排、袋裝碎石墊層施工一段時(shí)間初步具備護(hù)底和棱體塊石拋填作業(yè)條件后立即進(jìn)行此2項(xiàng)工作的開展。

1.2 一階段圍堤堤心及護(hù)面

一階段圍堤堤心及護(hù)面采用堤心石拋填施工,分4級(jí)進(jìn)行。第1級(jí)堤心石拋填完成后，進(jìn)行對(duì)應(yīng)部位后坡倒濾層、土工膜的和平臺(tái)護(hù)面層的施工作業(yè)，在第1級(jí)平臺(tái)護(hù)面施工的同時(shí)安排拋石船進(jìn)行第2級(jí)堤心石的拋填施工，同時(shí)安排施工人員進(jìn)行第2級(jí)后坡倒濾層、土工布?jí)浩碌仁┕ぷ鳂I(yè)，第3級(jí)、第4級(jí)與第2級(jí)施工工藝基本相同。

1.3 一階段圍堤上部結(jié)構(gòu)

當(dāng)一階段圍堤結(jié)構(gòu)全部完成并沉降穩(wěn)定后進(jìn)行防浪墻及堤頂?shù)缆返氖┕ぷ鳂I(yè)，防浪墻計(jì)劃采取分段跳倉澆筑的方式進(jìn)行作業(yè)。堤頂?shù)缆吩诜览藟就瓿珊筮M(jìn)行施工作業(yè)。

1.4 二階段圍堤部分

二階段圍堤在一階段圍堤初步具備吹砂條件后同步展開，二階段圍堤施工順序及流程基本與一階段圍堤一致。

在新建圍堤實(shí)際施工過程中，涉及選擇預(yù)留龍口位置的問題，如何保證預(yù)留的龍口水流平緩不影響圍堤安全、使吹填砂充分沉淀、砂的流失率最少和取得最佳經(jīng)濟(jì)效益是一大難題。

2 FVCOM模型

2.1 模型簡介

FVCOM (Finite-Volume Coastal Ocean Model)是三維有限體積近海模型。有限差分、有限元模式使用的是微分形式，而 FVCOM 是對(duì)控制方程進(jìn)行積分求解。在非結(jié)構(gòu)化網(wǎng)格中(與有限元方法相同)，用通量進(jìn)行計(jì)算(與有限差分法相同)，可以從數(shù)學(xué)上求解這些積分方程。本文通過應(yīng)用FVCOM模型，模擬梅山港區(qū)周圍海域的工程不同預(yù)留龍口方案下的潮流場變化。

2.2 模型建立

圖2 地形圖

根據(jù)研究工程范圍大小，選定梅山島附近海域作為研究區(qū)域，地形圖如圖2所示。由于研究范圍較小，選用最小間距為50m的三角形網(wǎng)格，垂向分為5層，本文采用M2分潮驅(qū)動(dòng)，通過獲得開邊界調(diào)和常數(shù)，從而計(jì)算得到水位進(jìn)行驅(qū)動(dòng)。公式如下：

式中ξ0——水位,m;

hi和gi——第i分潮的調(diào)和常數(shù)；

wi——各分潮的角速度;

t——時(shí)間,s;

fi——分潮的交點(diǎn)因子;

v0i——分潮的天文初位相;

ui——分潮的交點(diǎn)頂正角;

N0——分潮數(shù)。

2.3 模型驗(yàn)證

驗(yàn)證資料采用2016年5月15日4∶00至5月16日4∶00實(shí)測的水位、流速、流向資料，監(jiān)測點(diǎn)如圖3所示。圖4～圖6給出了模擬結(jié)果與驗(yàn)證資料的對(duì)比圖，其中實(shí)線為模擬值，虛線為實(shí)測值，實(shí)測值為實(shí)測潮位與流速的垂向平均值，其中T1站潮位誤差平均值為3.8cm，T2站潮位誤差平均值為6.2cm，流速誤差較大，T1、T2站分別為25%、34%。分析發(fā)現(xiàn)模擬結(jié)果與實(shí)測結(jié)果基本吻合，驗(yàn)證結(jié)果趨勢基本一致，但部分點(diǎn)出現(xiàn)偏差較大，一方面是模型模擬的不足，僅用M2分潮作為驅(qū)動(dòng)條件進(jìn)行模擬；僅選擇驗(yàn)證點(diǎn)附近的三角網(wǎng)格，具體地理位置可能存在偏差；另一方面，實(shí)測資料是用人工測量與記錄，存在一定誤差?？傮w而言模型可用。

圖3 驗(yàn)證點(diǎn)地理位置

圖4 T1、T2站位水位驗(yàn)證

圖5 T1、T2站位流速驗(yàn)證

圖6 T1、T2站位流向驗(yàn)證

3 預(yù)留龍口多方案比選

為更好地指導(dǎo)施工，將水動(dòng)力環(huán)境對(duì)圍堤建設(shè)的影響降至最低，本文對(duì)預(yù)留龍口位置進(jìn)行方案比選。本工程新建圍堤共2813m，為保證數(shù)值模擬的合理性以及可操作性，預(yù)留龍口有5處位置可供選擇，如圖7所示進(jìn)行布置(按照新建圍堤長度滿足500m處布設(shè))。但是考慮到梅山港區(qū)新建圍堤以及內(nèi)部陸域形成施工是從西南側(cè)向東北側(cè)推進(jìn)，若選用龍口3、龍口4、龍口5位置進(jìn)行預(yù)留，大大增加了新建圍堤陸拋石、后坡以及護(hù)面施工等工序的施工難度，且與吹砂施工沖突，直接造成吹填砂原材料的流失。因此預(yù)留龍口位置最終可選的只有在新建圍堤K1+1800處(龍口1)，以及K2+500處(龍口2)。本文首先根據(jù)龍口1和龍口2預(yù)留龍口方案，模擬出2個(gè)龍口方案在漲急時(shí)刻以及落急時(shí)刻的最大流速，掌握2個(gè)龍口方案對(duì)新建圍堤沖刷情況，用流速作為方案比選依據(jù)，為工程決策提供理論依據(jù)。

圖7 龍口預(yù)留位置

3.1 預(yù)留龍口1流場分布

圖8 龍口1處漲急時(shí)刻流場

圖9 龍口1處落急時(shí)刻流場

圖8、圖9分別繪制了新建海堤工程預(yù)留龍口1時(shí)的漲急時(shí)刻與落急時(shí)刻流場。漲潮時(shí)海流從外海的南側(cè)傳播而來，向北輸運(yùn)，海流經(jīng)由佛度島時(shí)一分為二，然后在六橫島北側(cè)匯合流至外海，落潮流與漲潮流方向大致相反。在工程區(qū)域：漲潮時(shí)，潮流在新建海堤海側(cè)順著新建圍堤由西南向東北側(cè)運(yùn)輸，由龍口1處流入工程區(qū)域內(nèi)，潮流在新建海堤及陸側(cè)之間形成逆時(shí)針環(huán)流，在龍口處潮流流速最大，漲急時(shí)刻龍口處最大流速36.67cm/s；落潮時(shí)，潮流方向與漲潮潮流方向大致相反，在龍口處潮流流速最大，落急時(shí)刻龍口處最大流速25.41cm/s。

3.2 預(yù)留龍口2流場分布

圖10、圖11分別繪制了新建海堤工程預(yù)留龍口2時(shí)的漲急時(shí)刻與落急時(shí)刻流場。漲潮時(shí)海流經(jīng)龍口2處由海側(cè)流入工程區(qū)域內(nèi)，潮流在新建海堤及陸側(cè)之間形成逆時(shí)針環(huán)流，在龍口處潮流流速最大，漲急時(shí)刻龍口處最大流速29.27cm/s；落潮時(shí)，潮流方向與漲潮潮流方向大致相反，在龍口處潮流流速最大，落急時(shí)刻龍口處最大流速27.38cm/s。

圖10 龍口2處漲急時(shí)刻流場

圖11 龍口2處落急時(shí)刻流場

4 水交換

4.1 水交換概念

在該工程中，新建圍堤陸側(cè)需進(jìn)行吹填砂，可引入水交換的概念，通過模擬不同龍口方案下水交換的快慢，來反映工程區(qū)域內(nèi)吹填砂流失率的大小，即水交換越快，吹填砂的流失率越大。

在模型計(jì)算出流場的情況下，開啟拉格朗日質(zhì)點(diǎn)追蹤模塊，從某一初始刻開始，在整個(gè)圍堤工程區(qū)域內(nèi)均勻布設(shè)一定量的質(zhì)點(diǎn)，然后統(tǒng)計(jì)各時(shí)刻通過龍口流出的質(zhì)點(diǎn)數(shù)X(t),用X(t)與初時(shí)刻釋放的質(zhì)點(diǎn)相除，計(jì)算得出整個(gè)工程區(qū)域的水交換率，計(jì)U(X,t)=dX/dt，水交換率越大說明水交換越快。

4.2 不同預(yù)留龍口方案下的水交換率

在新建圍堤與陸側(cè)之間的海域內(nèi)均勻布設(shè)1000個(gè)質(zhì)點(diǎn)，在模型運(yùn)行10d穩(wěn)定后一次性釋放，模擬龍口方案1與龍口方案2下的質(zhì)點(diǎn)移動(dòng)趨勢，如圖12、13 所示。不同預(yù)留龍口方案下的水交換率如圖14所示。

圖12 龍口方案1在1h、5h、10h時(shí)質(zhì)點(diǎn)分布

圖13 龍口方案2在1h、5h、10h時(shí)質(zhì)點(diǎn)分布

圖14 龍口方案與龍口方案2水交換率對(duì)比

對(duì)比兩種龍口方案下的質(zhì)點(diǎn)輸運(yùn)趨勢發(fā)現(xiàn)，起始質(zhì)點(diǎn)均經(jīng)由龍口處向東南側(cè)輸運(yùn)，后受到佛渡島阻擋，大部分質(zhì)點(diǎn)有向西側(cè)運(yùn)輸?shù)内厔?。且預(yù)留龍口方案2下的水交換率明顯低于預(yù)留龍口方案1。

5 總 結(jié)

該工程預(yù)留龍口雖然有5處，但由于陸域形成吹砂施工是從西南側(cè)向東北側(cè)推進(jìn)，為延長吹砂的流程，使吹填砂充分沉淀，提高沉降效果，減少吹填流失量，龍口設(shè)置必須遠(yuǎn)離出砂口，即龍口位置距西南側(cè)越遠(yuǎn)越好；為了避免新建圍堤陸拋石、后坡以及護(hù)面施工等工序的施工難度，龍口3、龍口4、龍口5方案不可行，龍口1和龍口2非常合適，在2者中選擇1個(gè)即可。

為此，利用FVCOM模型對(duì)梅山島海域進(jìn)行模擬，通過模擬計(jì)算得到水動(dòng)力場數(shù)據(jù)，并與實(shí)際水文資料驗(yàn)證基本吻合，說明模型能夠較好的模擬出梅山島海域的水動(dòng)力。

從模擬實(shí)驗(yàn)結(jié)果可以看出：

a.預(yù)留龍口2處漲急與落急流速明顯小于預(yù)留龍口1處，因此海流對(duì)龍口處沖刷的影響更小，更有利于圍堤的安全穩(wěn)定性與成本控制。

b.預(yù)留龍口方案2下的水交換率明顯低于預(yù)留龍口方案1，因此預(yù)留龍口2時(shí)水交換更慢，工程內(nèi)部的吹填砂流失速率會(huì)更小，有利于降本增效。

綜上所述，預(yù)留龍口2(新建圍堤K2+500處)方案比預(yù)留龍口1(新建圍堤K1+1800處)方案更適合該工程。