陸凡，沈良朵*，高郁，王晉寶

(1.浙江海洋大學(xué) 船舶與海運(yùn)學(xué)院，浙江 舟山 316021)

1 引言

在經(jīng)濟(jì)全球化快速發(fā)展的大背景下，世界各國(guó)對(duì)海洋資源的開(kāi)發(fā)利用越來(lái)越重視。特別地，我國(guó)大部分沿海城市在港口碼頭、物流集散中心和綠色石化基地的建設(shè)投入逐步加大[1]。然而，部分以圍填海方案來(lái)獲取工業(yè)用地的建設(shè)方式由于缺少科學(xué)的指導(dǎo)依據(jù)，對(duì)周邊海域的生態(tài)環(huán)境、水動(dòng)力條件和泥沙沖淤造成了不良影響[2]?；诖?，目前通過(guò)采用數(shù)值模擬分析法以水工建筑物對(duì)周邊海洋環(huán)境的影響范圍為目標(biāo)進(jìn)行科學(xué)預(yù)測(cè)，能夠合理有效地解決如何綠色開(kāi)發(fā)海洋資源的問(wèn)題。

圍填海工程前能夠科學(xué)地模擬出建設(shè)海域相關(guān)的潮流運(yùn)動(dòng)形態(tài)差異和泥沙沖淤的顯著地區(qū)，對(duì)工程實(shí)施有著積極的推動(dòng)作用。因此，一些國(guó)內(nèi)外專(zhuān)家學(xué)者借助數(shù)值模擬進(jìn)行了相關(guān)研究。任一晗等[3]基于FVCOM 三維水動(dòng)力數(shù)值模型，選取1984 年、2010 年、2019 年3 個(gè)代表年份，探討圍墾工程影響下舟山群島海域潮流結(jié)構(gòu)與潮能分布的時(shí)空變化狀況；Safavi等[4]采用MIKE3 平臺(tái)構(gòu)建了烏魯米耶湖的三維數(shù)值模型，研究了受湖中堤道分隔影響湖中南北部水中鹽度的變化情況；竇明等[5]利用FVCOM 三維海洋數(shù)值模型，模擬嵊泗大洋山圍墾工程建設(shè)填海前后的流場(chǎng)及泥沙運(yùn)動(dòng)，對(duì)該工程附近海域潮流及各層泥沙分布變化進(jìn)行了探討；劉嘉星和劉長(zhǎng)根[6]利用驗(yàn)證后的三維潮流數(shù)學(xué)模型ROMS 探討了天津港附近海域的潮流特征、潮流的分層特點(diǎn)和圍海造地引起的潮流場(chǎng)變化情況；劉金鵬等[7]以龍口市人工島為例，利用驗(yàn)證合理的潮流場(chǎng)和波浪場(chǎng)數(shù)學(xué)模型，模擬分析了在極限風(fēng)況條件下人工島群周?chē)Ｓ蛩h(huán)境的變化特征及人工島岸線(xiàn)的波高特征值的變化規(guī)律。Zamani 和Koch[8]以伊朗西南部Maroon 水庫(kù)為研究對(duì)象，通過(guò)與實(shí)測(cè)數(shù)據(jù)進(jìn)行比較的基礎(chǔ)上，對(duì)這兩種三維水動(dòng)力模型的適用性進(jìn)行了對(duì)比研究；Soltani等[9]通過(guò)Mike21，2011 Version 模型，模擬了海岸帶的波浪和水流，以確定奈班德海灣地區(qū)的泥沙輸移速率和方向；Elshemy 等[10]使用MIKE21 建模系統(tǒng)開(kāi)發(fā)了曼扎拉湖的水動(dòng)力模型和水質(zhì)模型，研究了未來(lái)氣候變化對(duì)湖泊水動(dòng)力和水質(zhì)特征的影響。

雖然前人對(duì)圍墾工程的數(shù)值模擬研究已經(jīng)取得不錯(cuò)的結(jié)果，但大多數(shù)都是以FVCOM 三維模型或者M(jìn)IKE 二維水動(dòng)力模型對(duì)相關(guān)海域的水動(dòng)力和溫鹽變化進(jìn)行了細(xì)致研究，而采用MIKE3 和MIKE21 這兩種水動(dòng)力模型對(duì)國(guó)內(nèi)近海岸島嶼圍墾工程的對(duì)比性研究較少?；诖?，本文以舟山綠色石化基地三期圍填海工程為研究對(duì)象，基于MIKE3 平臺(tái)構(gòu)建了工程前后三維水動(dòng)力數(shù)學(xué)模型，并考慮到MIKE21 中的二維模型也可以分析水平方向上的潮流運(yùn)動(dòng)情況，故對(duì)兩種模型進(jìn)行了潮位驗(yàn)證結(jié)果精度的比較分析。但結(jié)合到舟山群島海域水道縱橫、地形結(jié)構(gòu)復(fù)雜的特性，水平方向上的二維研究無(wú)法呈現(xiàn)出三維模型展現(xiàn)出的流場(chǎng)空間結(jié)構(gòu)變化，因此在三維模型下分析了采用“裁彎取直”和水道封堵方案下對(duì)潮流場(chǎng)特征的影響，并從表層、中間層和底層分別對(duì)工程前后的流速流向影響進(jìn)行三維垂向剖析研究。

2 工程概況

本工程位于舟山群島和岱山以西的大、小魚(yú)山島海域，北瀕岱衢洋、南接灰鱉洋，周邊島嶼岸線(xiàn)曲折，港灣眾多。水下地形南北深、東西淺，海底高程基本在-65～-1 m，水下地形起伏明顯[11]。三期圍填海工程地處于121°55′11″～121°59′37″E 之間，其用海面積約為18.41 km2，成陸面積為17.09 km2。工程涉及海域潮流運(yùn)動(dòng)以往復(fù)流為主，而且基本為平行岸線(xiàn)方向的往復(fù)流。漲潮流主要為來(lái)自東海的前進(jìn)波經(jīng)由舟山群島之間的狹道進(jìn)入杭州灣，其中岱山海域的漲潮流主流較為平順，整體上呈現(xiàn)由東南流向西北方向。落潮流則為杭州灣的落潮水流經(jīng)舟山群島之間的狹道流入外海，落潮流基本與漲潮流一致，岱山周邊落潮流整體較為平順。海區(qū)內(nèi)泥沙主要來(lái)源于長(zhǎng)江口以及杭州灣外泄泥沙隨浙東沿岸流的輸移擴(kuò)散，含沙量總體較大。

3 三維水動(dòng)力模型

3.1 控制方程

三維模型使用的是分層網(wǎng)格，即在水平域中使用非結(jié)構(gòu)網(wǎng)格，而在垂直域中使用結(jié)構(gòu)化網(wǎng)格。垂直網(wǎng)格是分為sigma 坐標(biāo)或組合的sigma/z 級(jí)坐標(biāo)。本文采用的是考慮了自由曲面的sigma 坐標(biāo)。使用sigma坐標(biāo)最重要的優(yōu)勢(shì)是它們能夠準(zhǔn)確地表示水深測(cè)量，并在海床附近提供一定的分辨率[12]。采用單元中心的有限體積方法進(jìn)行了原始方程的空間離散化，方程式如下：

式中，t為時(shí)間；x、y、z為笛卡爾坐標(biāo)；η為表面高度；d為靜止水深；h=η+d為總水深；u、v、w為x、y、z方向的速度分量；f=2Ωsinφ，sinφ為科里力參數(shù)（Ω為角轉(zhuǎn)速，φ為地理緯度）；g為重力加速度；ρ為水密度；sxx、sxy、syx、syy為輻射應(yīng)力張量的分量；νt為垂直湍流（或渦流）黏度；pa為大氣壓力；ρ0為水的參考密度；S為由于點(diǎn)源引起的排放量大?。籾s和νs為水被排放到環(huán)境水中的速度；Fu和Fv分別為x和y方向上的水平應(yīng)力分項(xiàng)。

3.2 模型建立

3.2.1 計(jì)算區(qū)域及網(wǎng)格

考慮到工程區(qū)域附近漲落潮的流路經(jīng)杭州灣和部分東海海域，因此所選取區(qū)域范圍遠(yuǎn)大于三期工程計(jì)算區(qū)域，北界至32°N 左右，南界至28.5°N，東到125°E 附近，潮流界頂點(diǎn)江陰作為長(zhǎng)江上游邊界，澉浦潮位站附近作為杭州灣上游邊界，計(jì)算區(qū)域位于120°～126°E 之間。三維模型采用sigma 網(wǎng)格，垂向均勻分10 層，工程附近海域網(wǎng)格進(jìn)行加密，網(wǎng)格最小尺度為18 m，工程外海區(qū)域，網(wǎng)格距為8 000 m，不同尺度網(wǎng)格之間通過(guò)設(shè)置實(shí)現(xiàn)平滑過(guò)渡，共計(jì)52 007 個(gè)節(jié)點(diǎn)，99 195 個(gè)單元。對(duì)于二維模型，共計(jì)41 041 個(gè)節(jié)點(diǎn)，78 473 個(gè)單元。模型計(jì)算區(qū)域及水深如圖1 所示。

圖1 模型區(qū)域和水深Fig.1 Computational domain and water depth

3.2.2 模型參數(shù)設(shè)置

模型計(jì)算時(shí)間步長(zhǎng)設(shè)置根據(jù)克朗數(shù)CFL 條件進(jìn)行動(dòng)態(tài)調(diào)整，在確保模型計(jì)算穩(wěn)定進(jìn)行下，平均時(shí)間步長(zhǎng)為2 s，最短時(shí)間步長(zhǎng)為0.01 s。干濕水深判別設(shè)置為干水深0.005 m，淹沒(méi)水深0.05 m，濕水深0.1 m。這有利于保護(hù)模型穩(wěn)定運(yùn)行。水平渦黏系數(shù)經(jīng)率定取0.28，垂直渦黏系數(shù)使用對(duì)數(shù)公式確定。通過(guò)粗糙高度來(lái)表征底床糙率，率定后取0.003 m。邊界采用水位控制，即用潮位預(yù)報(bào)的方法得到開(kāi)邊界條件。外海開(kāi)邊界潮位根據(jù)由16 個(gè)主要分潮（M2、S2、N2、K2、K1、O1、P1、Q1、MU2、NU2、T2、L2、2N2、J1、M1和OO1）的調(diào)和常數(shù)，以開(kāi)邊界潮位預(yù)報(bào)公式計(jì)算得到：

式中，A0為平均海面；Fi和 (ν0+u)i為天文要素；Hi和 gi為調(diào)和常數(shù)。在閉邊界取流速的法向?qū)?shù)為0，在潮灘區(qū)采用漫灘邊界處理。

4 模型驗(yàn)證

由于本工程周邊缺少實(shí)測(cè)數(shù)據(jù)，故本次模型驗(yàn)證采用了距魚(yú)山島較近的小洋山海域，由中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司在水文測(cè)驗(yàn)期間的潮流實(shí)測(cè)數(shù)據(jù)。其中，分大中小潮對(duì)臨時(shí)測(cè)站T1 和T6 進(jìn)行潮位驗(yàn)證（大潮測(cè)流時(shí)間為2015 年5 月4 日至5 日；中潮測(cè)流時(shí)間為2015 年5 月7 日至8 日；小潮測(cè)流時(shí)間為2015 年5 月12 日至13 日）；同樣地，對(duì)定點(diǎn)測(cè)站N1、N14 以及走航測(cè)站W(wǎng)P05 流速流向數(shù)據(jù)進(jìn)行潮流分層驗(yàn)證，各測(cè)站位置見(jiàn)圖2。

圖2 驗(yàn)證點(diǎn)位置Fig.2 Verification point location

采用MIKE3 的水動(dòng)力模型計(jì)算數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證，潮位驗(yàn)證結(jié)果如圖3 所示。圖中顯示，從整體上看T1 測(cè)站和T6 測(cè)站潮位擬合較好，其中小潮階段擬合效果最好，相對(duì)誤差結(jié)果在±10%左右，個(gè)別高潮位處相對(duì)誤差稍大?？傮w來(lái)看，三維模型模擬的潮位驗(yàn)證是符合要求的。

圖3 潮位驗(yàn)證Fig.3 Tide level verification

對(duì)潮流測(cè)站N1、N14 和WP05 分為表層（第10層）、中層（第5 層）和底層（第1 層）進(jìn)行流速流向驗(yàn)證。在分層驗(yàn)證下，可以清楚看到表層小潮階段的流速擬合結(jié)果是相對(duì)較大的，其他各層擬合結(jié)果都是極好的。表層模擬結(jié)果如圖4 所示，各測(cè)站僅在高潮位由于流量加大的影響導(dǎo)致流速擬合結(jié)果有0.5 m/s左右的誤差，其他時(shí)間段擬合相對(duì)誤差控制在8%左右，表層流向模擬結(jié)果令人滿(mǎn)意。中間層模擬結(jié)果如圖5 所示，各測(cè)站除高潮位處流速存在0.15～0.3 m/s 的誤差，其余各時(shí)間點(diǎn)相對(duì)誤差在10%以?xún)?nèi)，流向擬合結(jié)果較滿(mǎn)意。底層模擬結(jié)果如圖6 所示，大潮和中潮階段的模擬流速僅在N1 測(cè)站的高潮位時(shí)有0.35 m/s 的誤差，其余擬合較好，小潮階段除在個(gè)別時(shí)間點(diǎn)有0.45 m/s 的誤差外，其余各點(diǎn)存在0.15～0.25 m/s 的誤差范圍；底層流向模擬結(jié)果較好，流向相對(duì)誤差在7%以?xún)?nèi)，極個(gè)別點(diǎn)處由于測(cè)量過(guò)程存在失誤而出現(xiàn)的誤差，不予考慮。

圖4 表層流速流向驗(yàn)證Fig.4 Verification of flow velocity and direction in the surface layer

圖5 中間層流速流向驗(yàn)證Fig.5 Verification of flow velocity and direction in the middle layer

圖6 底層流速流向驗(yàn)證Fig.6 Verification of flow velocity and direction in the bottom layer

為對(duì)模型驗(yàn)證結(jié)果有更加科學(xué)性的認(rèn)識(shí)，此處通過(guò)Willmott[13]所提出的skill 方程進(jìn)行評(píng)估，表達(dá)式如下：

式中，Pi是模擬值；Oi是 實(shí)測(cè)值；是實(shí)測(cè)平均值。通過(guò)d值判斷模型效率，d＞0.65 時(shí)，模型效率極好；0.5＜d≤ 0.65，表示模型效率非常好；0.2＜d≤ 0.5，表示模型效率好；d≤ 0.2，表示模型效率差。

結(jié)果表明，N1 和WP05 測(cè)站各層的大中小潮期間流速d值在0.80～0.97 之間，流向d值區(qū)間為0.73～0.99，其模型效率都是極好的。限于篇幅，僅具體列出了N14 測(cè)點(diǎn)各層的模型效率值，具體見(jiàn)表1。

表1 模型效率系數(shù)Table 1 The index of agreement

5 三維結(jié)果分析

5.1 MIKE3 與MIKE21 對(duì)比分析

對(duì)MIKE3 與MIKE21 兩種擬合效果進(jìn)行對(duì)比分析，分別將這兩種水動(dòng)力模型的計(jì)算數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證，潮位對(duì)比驗(yàn)證結(jié)果如圖7 所示，由圖可知：兩種模型的潮位計(jì)算值與實(shí)測(cè)值的吻合度較高，特別是在小潮期間，其最大誤差僅有0.15 m 左右。大潮和中潮期間，三維數(shù)值結(jié)果與實(shí)測(cè)值在個(gè)別高潮處有0.5 m 左右的誤差，總體上看，在將二維模型中的初始條件同樣地設(shè)置于三維模型的情況下，二維數(shù)值結(jié)果更接近于實(shí)測(cè)值。

圖7 潮位對(duì)比Fig.7 Tide level comparison

考慮到理論方面上MIKE3 中間層的潮流模擬結(jié)果與MIKE21 的平均流速流向的結(jié)果相似，故在此處進(jìn)行了流速流向的對(duì)比驗(yàn)證。中間層模擬結(jié)果對(duì)比如圖8 所示，圖中各測(cè)站流向驗(yàn)證結(jié)果良好。本文選擇的水深0.6H的實(shí)測(cè)數(shù)據(jù)與垂向平均流速基本保持一致，在大潮和中潮階段兩種模型擬合結(jié)果較好，最大相對(duì)誤差在20%左右，相較而言，小潮階段的高潮位處誤差稍大，相對(duì)誤差在25%左右，其余各時(shí)間點(diǎn)兩者模擬結(jié)果均良好。通過(guò)這4 組數(shù)據(jù)的對(duì)比發(fā)現(xiàn)：二維模型的大部分計(jì)算值更加接近于實(shí)測(cè)值，尤其是在漲急時(shí)刻，原因可能是二維模型是取值于垂向水深平均的流速進(jìn)行模擬，而本文三維模型未考慮風(fēng)場(chǎng)應(yīng)力的影響，表層流速計(jì)算誤差會(huì)對(duì)三維流速計(jì)算結(jié)果有影響?？傮w上看，三維模型中間層的模擬結(jié)果與二維計(jì)算結(jié)果是最接近的，且中間層擬合效果最為出色，模擬精度基本滿(mǎn)足《海岸與河口潮流泥沙數(shù)值模擬技規(guī)程》的要求，模型可用于工程后的預(yù)測(cè)等各項(xiàng)工作[14-15]。

5.2 工程前后潮流場(chǎng)分析

從計(jì)算的潮流場(chǎng)分析得出，漲急時(shí)刻外海潮波沿東南向西北方向前進(jìn)。在大魚(yú)山的東南側(cè)受島嶼阻擋，潮流分成兩股，分別沿著大魚(yú)山島東南側(cè)向南北流入杭州灣。大魚(yú)山島區(qū)域漲潮時(shí)刻受峽道效應(yīng)，南岸流速大于北岸，落潮時(shí)刻則北岸較大，工程南北兩端的磯頭岸線(xiàn)對(duì)漲落潮流起到了挑流的作用，而且兩側(cè)都為深槽，因此漲落潮流速相對(duì)較大，但整體還較為平整。其中，漲潮時(shí)平均流向?yàn)?30°～300°，落潮時(shí)平均流向?yàn)?0°～110°，漲落潮平均流向約為180°，基本符合往復(fù)流性質(zhì)。

考慮到圍墾工程所處的大小魚(yú)山島周邊范圍內(nèi)島嶼錯(cuò)落，水道縱橫，潮流運(yùn)動(dòng)情況復(fù)雜，故對(duì)工程前后計(jì)算域內(nèi)大潮漲落時(shí)刻的表層流場(chǎng)變化進(jìn)行具體的對(duì)比分析，同時(shí)給出了工程后流速的等值線(xiàn)圖以更直觀(guān)看出流速大小變化，如圖9 所示。整體上看，此圍填海工程前后對(duì)工程周邊海域的潮流運(yùn)動(dòng)沒(méi)有較大影響，是明顯的往復(fù)流形式。在大潮漲急階段，由于受?chē)鷫üこ獭敖貜澣≈薄狈椒ǖ挠绊?，處于工程?hào)|北側(cè)的岸線(xiàn)外伸，東北側(cè)岸線(xiàn)存在一個(gè)向內(nèi)的45°轉(zhuǎn)角，受岸線(xiàn)角度和洋流反射作用，較于工程前平行于岸線(xiàn)的潮流運(yùn)動(dòng)方向，此處產(chǎn)生一個(gè)直徑約為700 m 的逆時(shí)針旋渦；魚(yú)山島的西北側(cè)和東南側(cè)受?chē)钣绊?，潮流運(yùn)動(dòng)遇到外擴(kuò)海岸線(xiàn)，此處潮流運(yùn)動(dòng)方向變化呈現(xiàn)出30°轉(zhuǎn)角；大范圍流場(chǎng)前后變化不大。大潮落急階段，工程西北角處的一處小島嶼圍填時(shí)與主島進(jìn)行了連接處理，將此處原先經(jīng)狹道效應(yīng)的高速流與經(jīng)岸線(xiàn)折射過(guò)來(lái)的洋流形成的旋轉(zhuǎn)流影響削弱，故旋轉(zhuǎn)流轉(zhuǎn)變成了與岸線(xiàn)彎口一致的運(yùn)動(dòng)形式；魚(yú)山島東北側(cè)、南側(cè)以及西側(cè)落潮流受外延海岸線(xiàn)的挑流作用，工程后流場(chǎng)運(yùn)動(dòng)均出現(xiàn)30°～45°的轉(zhuǎn)角變化，這表明處于迎流面的岸線(xiàn)變化對(duì)流場(chǎng)的改變起主導(dǎo)作用。

圖9 表層潮流場(chǎng)對(duì)比（上：漲急 下：落急）Fig.9 The comparison of tidal current field on the surface layer (up:rising;down:falling)

5.3 垂直流場(chǎng)分析

為更立體的表現(xiàn)出工程區(qū)域的流速流向情況，本節(jié)同時(shí)選取了N1、N14 和WP05 這3 個(gè)典型的測(cè)站數(shù)據(jù)進(jìn)行垂直結(jié)構(gòu)的分析。由圖10 分析得：N1 測(cè)站的流速在垂直方向上的實(shí)測(cè)值和漲落急階段的計(jì)算值各自都較接近，并在中間層附近的水深0.6H處僅有4.6%左右的相對(duì)誤差，計(jì)算與實(shí)測(cè)的最大誤差出現(xiàn)在表層附近，相對(duì)誤差在-25%左右，漲落急各分層流速較為一致，整體上都呈現(xiàn)出流速隨水深增加而減小，且幅度較大，可能是受到底部海床的摩阻力較大的影響；對(duì)于各分層流向趨于一致，且與實(shí)測(cè)數(shù)據(jù)擬合較好。N14 測(cè)站處漲落急時(shí)刻的流速計(jì)算值與實(shí)測(cè)值最大誤差出現(xiàn)在表層，其相對(duì)誤差為-28.01%，同時(shí)此測(cè)站的漲急流速計(jì)算值大于落急時(shí)刻計(jì)算值，這可能是由于此處網(wǎng)格布置較疏的影響，加上缺少風(fēng)場(chǎng)實(shí)測(cè)資料，未考慮風(fēng)場(chǎng)的原因所導(dǎo)致。WP05 測(cè)站可能處在洋山港與殼子山之間的狹窄水道里，由寬闊海域進(jìn)入該水道，過(guò)水?dāng)嗝嫱蝗粶p小，故同樣是在表層處的流速與實(shí)測(cè)數(shù)據(jù)有-30%的相對(duì)誤差，而0.6H以下水深的計(jì)算流速與實(shí)測(cè)值的誤差范圍為13～22 cm/s，且漲急時(shí)刻流速小于落急流速，同樣地各分層的流向是近乎一致的。通過(guò)對(duì)垂直流場(chǎng)的分析，可以清楚地得到流速與水深呈負(fù)相關(guān)，流向受水深影響近乎不計(jì)。而針對(duì)表層擬合較大問(wèn)題，初步考慮是受到了風(fēng)場(chǎng)和波浪對(duì)潮流的驅(qū)動(dòng)作用。王世澎等[16]在二維情況下波浪對(duì)潮流場(chǎng)的數(shù)值分析中指出，對(duì)于流速過(guò)程，波浪在沒(méi)有破碎的情況下對(duì)流速過(guò)程的影響比對(duì)潮位過(guò)程的影響大?？紤]到圍填海工程后，周?chē)Ｓ蚴芫植刻袅髯饔茫又ɡ溯椛鋺?yīng)力影響，波浪會(huì)在近岸破碎產(chǎn)生沿岸流疊加到原有潮流場(chǎng)上，對(duì)表層流速影響較大。針對(duì)此問(wèn)題，將會(huì)在以后的研究中通過(guò)波浪耦合進(jìn)行率定工作。

圖10 測(cè)站垂向流速流向分布Fig.10 Flow velocity and direction in vertical direction

6 小結(jié)

本文以舟山綠色石化三期工程為背景，采用MIKE3建立寧波舟山及其周邊附近海域三維水動(dòng)力模型。通過(guò)與實(shí)測(cè)數(shù)據(jù)進(jìn)行潮流驗(yàn)證，對(duì)比了二維模型與三維模型驗(yàn)證結(jié)果的區(qū)別，并主要對(duì)圍墾工程前后的水動(dòng)力變化情況進(jìn)行了詳細(xì)分析。最終得到如下幾點(diǎn)結(jié)論：

（1）以MIKE21 FM 的模型參數(shù)設(shè)置于MIKE3 中，二維潮位驗(yàn)證效果更好一點(diǎn)，但兩種模型的驗(yàn)證都是符合要求的；兩種模型流速擬合的最小相對(duì)誤差都控制在10%左右，通過(guò)4 組數(shù)據(jù)的對(duì)比，更直觀(guān)地驗(yàn)證了二維數(shù)模數(shù)據(jù)是取值于垂向平均流速的理論結(jié)果，而三維中間層的模擬結(jié)果與二維計(jì)算結(jié)果則是非常接近。

（2）圍墾工程前后并未大面積改變大魚(yú)山島附近的潮流運(yùn)動(dòng)，僅在工程?hào)|北側(cè)、東南側(cè)和西北角由于圍墾工程將小島嶼之間的水道連接起來(lái)，并在“截彎取直”的圍墾方案影響下，使得局部區(qū)域潮流流態(tài)發(fā)生改變。

（3）以MIKE3 對(duì)工程附近垂向分布上的流速流向進(jìn)行了分層結(jié)構(gòu)分析，整體上流速都呈現(xiàn)出隨水深的增加而減小，可能由于未考慮風(fēng)場(chǎng)和波浪場(chǎng)的影響，導(dǎo)致部分測(cè)站表層計(jì)算數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)誤差稍大，對(duì)于流向各層十分一致。