半封閉海灣共振周期計(jì)算與研究

楊萬(wàn)康, 伊小飛, 陳忠彪

?

半封閉海灣共振周期計(jì)算與研究

楊萬(wàn)康1, 伊小飛2, 陳忠彪3

(1. 國(guó)家海洋局第二海洋研究所工程海洋學(xué)重點(diǎn)實(shí)驗(yàn)室, 浙江杭州 310012; 2. 廣東海洋大學(xué)海洋與氣象學(xué)院, 廣東湛江 524088; 3. 南京信息工程大學(xué)海洋科學(xué)學(xué)院, 江蘇南京 210044)

三門灣海域水波共振現(xiàn)象較為顯著, 容易引發(fā)海洋災(zāi)害。針對(duì)三門灣這一典型半封閉海灣, 分別從理論推導(dǎo)和數(shù)值模擬兩方面研究了三門灣的共振周期。研究發(fā)現(xiàn): 理論推導(dǎo)結(jié)果可以用于對(duì)海灣整體共振周期的粗略估計(jì), 數(shù)值模擬結(jié)果更加精確, 也能充分反映海灣內(nèi)不同水域的共振特征。數(shù)值模擬結(jié)果表明, 三門灣內(nèi)各水域普遍存在3~4個(gè)共振周期模態(tài), 灣內(nèi)各水道第一模態(tài)和第三模態(tài)共振周期數(shù)值解與理論值較為接近, 石浦水道與外海相聯(lián)通, 共振周期不顯著。三門灣內(nèi)各水域共振周期第一模態(tài)對(duì)應(yīng)的振幅最大, 然后依次遞減, 但灣頂?shù)那嗌礁邸⑸唧此赖诙B(tài)共振周期及對(duì)應(yīng)的振幅值與第一模態(tài)相差較小, 基本呈現(xiàn)雙峰結(jié)構(gòu)。共振振幅由灣外向三門灣頂部逐漸增加, 尤其是灣頂處振幅增益比較顯著。通過(guò)本研究可以為三門灣的防災(zāi)減災(zāi)提供科學(xué)依據(jù)。

半封閉海灣; 共振; 理論推導(dǎo); 數(shù)值模擬

共振是海洋中普遍存在的現(xiàn)象。尤其是對(duì)于沿岸的海灣地區(qū), 共振現(xiàn)象會(huì)造成沿岸水位的暴增, 造成嚴(yán)重的海洋災(zāi)害, 因此是海洋動(dòng)力學(xué)研究的熱點(diǎn)之一。

海灣的共振周期取決于海灣的水深和形狀, 關(guān)于海灣共振周期的研究, 19世紀(jì)20年代就有了著名的梅立恩公式, 這一公式的不足之處是要求海灣形狀和深度具有很規(guī)則的變化或簡(jiǎn)單地看作常數(shù), 這與實(shí)際情況有所出入。Proudman[1]認(rèn)為芬迪灣較強(qiáng)的半日潮與潮波共振有關(guān), 即1/4波長(zhǎng)共振理論; 方國(guó)洪等[2]通過(guò)研究發(fā)現(xiàn)北部灣較強(qiáng)的全日潮是由于共振造成的, 通過(guò)理論推導(dǎo)得到海灣共振周期約為28.9 h, 與全日潮的周期較為接近, 因此全日潮振幅的增益較為顯著。Cummins等[3]認(rèn)為哈德遜海峽的共振情況與半波長(zhǎng)共振較為接近, 考慮了昂加瓦灣和科氏力的影響, 最終估算共振周期約為9.5~10.5 h。王鐘桾等[4]利用二維數(shù)值模式, 計(jì)算了芝罘灣的固有周期, 與實(shí)測(cè)資料吻合較好。Webb[5]提出了一種邊值方法用于計(jì)算海灣的共振周期, 通過(guò)在開邊界處添加不同周期的波動(dòng), 然后選取灣內(nèi)振幅極大值對(duì)應(yīng)的周期為共振周期。Grean等[6]利用數(shù)值模型計(jì)算了胡安德富卡海峽和喬治亞海峽的共振周期約為16 h, 認(rèn)為Helmholtz共振模型更符合該海區(qū)情況。Wilson[7]認(rèn)為, 沿岸假潮現(xiàn)象在許多情況下是一種強(qiáng)迫共振或者近似共振現(xiàn)象, 當(dāng)長(zhǎng)周期波進(jìn)入灣口時(shí), 是比較明顯的強(qiáng)迫振動(dòng), 當(dāng)外強(qiáng)迫消失后, 就轉(zhuǎn)變?yōu)樽杂烧鹗? 而且是以海灣的固有周期繼續(xù)振動(dòng)。Wong[8]在長(zhǎng)島海灣的研究中發(fā)現(xiàn), 由于共振效應(yīng), 半日分潮的振幅增加了4倍; Clarke等[9]通過(guò)簡(jiǎn)化拉普拉斯方程, 推導(dǎo)了共振發(fā)生的前提條件; Bertin等[10]采用大陸架共振理論模型和數(shù)值實(shí)驗(yàn)兩種方法研究發(fā)現(xiàn)數(shù)值模型計(jì)算結(jié)果能夠較好的匹配實(shí)際現(xiàn)象。一般而言, 對(duì)于實(shí)際海域共振周期的計(jì)算, 理論方法主要基于1/4波長(zhǎng)和 Helmholtz模型共振; 數(shù)值計(jì)算方面主要有簡(jiǎn)正模態(tài)方法和邊值方法[11]。

三門灣位于浙江象山縣南田島和三門縣牛頭山之間, 是一個(gè)西北-東南方向的半封閉海灣, 從灣口到灣頂長(zhǎng)約40 km, 灣內(nèi)多數(shù)地區(qū)寬10 km, 除東南方向通過(guò)石浦水道與外海連通外, 其余則為陸地和海島包圍, 具體地理位置如圖1所示。屬于較為典型的半封閉海灣, 潮汐周期受到外海的支配, 當(dāng)潮波進(jìn)入海灣后, 受到灣頂?shù)姆瓷? 從而形成駐波, 產(chǎn)生共振。此前關(guān)于三門灣共振特征的研究較少, 本文將分別采用理論推導(dǎo)和數(shù)值模擬方法來(lái)研究分析三門灣內(nèi)的共振周期, 通過(guò)建立高分辨的數(shù)值模型, 對(duì)三門灣內(nèi)各水道的共振周期特征進(jìn)行精細(xì)化研究。

圖1 三門灣地理位置及水深(單位: m)示意圖

1 理論推導(dǎo)

據(jù)文獻(xiàn)[12], 早在1828年, Merian在研究湖泊共振時(shí)提出了1/4波長(zhǎng)共振理論, 針對(duì)只存在一個(gè)封閉波節(jié)的封閉水域, 共振周期計(jì)算公式為:

其中,、、、分別表示海灣長(zhǎng)度、任意整數(shù)、重力加速度、水深。

經(jīng)典的梅里恩公式要求研究區(qū)域?yàn)橥耆忾]的海灣、湖泊, 由于三門灣為半封閉海灣, 因此需要對(duì)梅里恩公式進(jìn)行一定的改進(jìn)才能適用于三門灣海域。

根據(jù)上述兩個(gè)方程可以得到

這是典型的波動(dòng)方程, 解的形式一般為

利用關(guān)系式

則代入公式(6)可以得到:

所以共振發(fā)生的條件是

假設(shè)三門灣長(zhǎng)度為38 km, 平均水深取為7 m, 其共振周期理論結(jié)果如表1所示, 由理論推導(dǎo)的結(jié)果可知, 三門灣的共振周期理論解有無(wú)數(shù)個(gè), 周期最大的為5 h, 其余分別為1.67, 1 h, ……, 由于理論計(jì)算假設(shè)水深為固定值, 三門灣實(shí)際形狀也簡(jiǎn)化為矩形海灣, 因此結(jié)果可能有一定的誤差, 后面將利用更為準(zhǔn)確的數(shù)值模擬方法推算三門灣的共振周期。

2 數(shù)值模擬方法

理論方法通過(guò)概化海灣形狀、水深地形等條件, 可以用于許多海域共振周期的粗略估計(jì)。但是由于推導(dǎo)過(guò)程中采用了大量的近似及假設(shè)條件, 與實(shí)際情況存在一定的誤差, 尤其是理論推導(dǎo)方法不能反映海灣內(nèi)各局部地區(qū)的共振周期特征, 因此根據(jù)海灣內(nèi)的實(shí)際情況, 本文采用數(shù)值模擬方法對(duì)三門灣內(nèi)的共振周期特征進(jìn)行研究分析。

表1 三門灣共振周期理論解

2.1 模型設(shè)置與驗(yàn)證

本文數(shù)值模型采用了MIKE21數(shù)學(xué)模型, 該模型可以很好地模擬河流、海灣及海洋的水流、波浪、泥沙及水質(zhì)環(huán)境等, 在近岸港灣地區(qū)得到廣泛應(yīng)用[13-14]。本模型計(jì)算區(qū)域如圖2所示, 網(wǎng)格采用非結(jié)構(gòu)三角形網(wǎng)格, 共有100 079個(gè)網(wǎng)格單元, 海灣內(nèi)分辨率為20 m左右, 水深地形采用實(shí)際工程測(cè)量水深與海圖資料, 計(jì)算步長(zhǎng)采用動(dòng)態(tài)步長(zhǎng), 范圍為0.01~1 s。模型驗(yàn)證采用三門灣內(nèi)的實(shí)測(cè)潮位站數(shù)據(jù), T1、T2、T3潮位站位置如圖3所示。模型驗(yàn)證結(jié)果如圖4所示, 可以看到模擬振幅和實(shí)測(cè)振幅符合較好, 相位誤差較小, 表明該模型能準(zhǔn)確反映海灣內(nèi)的水位變化特征。

圖2 研究區(qū)域網(wǎng)格劃分示意圖

2.2 數(shù)值模擬結(jié)果與分析

各代表計(jì)算點(diǎn)隨開邊界入射波周期變化的振幅曲線如圖5所示, 由圖可知, 各計(jì)算點(diǎn)振幅曲線存在多個(gè)極大值, 各曲線振幅極大值對(duì)應(yīng)于共振周期的不同模態(tài), 其中共振周期數(shù)值最大的為第一模態(tài), 隨著模態(tài)的遞增, 共振周期的數(shù)值依次遞減。不同模態(tài)的共振周期統(tǒng)計(jì)結(jié)果如表2所示。由數(shù)值模擬結(jié)果可知:

圖3 海灣內(nèi)代表性計(jì)算點(diǎn)及實(shí)測(cè)潮位站分布圖

圖4 模擬潮位與實(shí)測(cè)潮位對(duì)比

圖5 各計(jì)算點(diǎn)振幅隨周期的變化

表2 各水道共振周期數(shù)值模擬結(jié)果(單位: h)

1) 三門灣內(nèi)各水域的實(shí)際共振周期模態(tài)要遠(yuǎn)少于理論推導(dǎo)得到的模態(tài)個(gè)數(shù), 至多存在4~5個(gè)模態(tài), 這是由于數(shù)值模擬考慮了實(shí)際的岸線阻擋及底摩擦消耗等因素, 尤其是理論推導(dǎo)中的高頻共振在實(shí)際條件下不會(huì)發(fā)生。

2) 除了白礁水道和石浦水道外, 其他區(qū)域的共振周期第一模態(tài)與理論值相差較小, 青山港和蛇蟠水道與理論值一致, 為5.0 h, 健跳水道和三門灣中心位置都為5.5 h。石浦水道由于與外海相聯(lián)通, 并非一個(gè)完整的半封閉水道, 因此共振周期第一模態(tài)不是非常顯著, 大約為6.5 h。白礁水道位于三門灣內(nèi)部, 但是其走向與三門灣走向夾角較大, 幾乎呈現(xiàn)垂直狀態(tài), 相當(dāng)于是兩個(gè)半封閉海灣的嵌套, 因此共振周期要大于三門灣的理論值。

3) 數(shù)值模擬結(jié)果的第二模態(tài)共振周期與理論值相差較大, 灣頂處的青山港, 蛇蟠水道第二模態(tài)周期值為3.0~3.17 h, 而理論推導(dǎo)無(wú)此共振周期結(jié)果。石浦水道是因?yàn)榕c外界水流的相互作用存在一個(gè)2.5 h的周期。數(shù)值模擬結(jié)果顯示各水道的第三模態(tài)周期為1.67 h, 與理論值計(jì)算結(jié)果基本一致。白礁水道、健跳水道和三門灣中心第4模態(tài)周期為0.67 h, 青山港和蛇蟠水道第5模態(tài)共振周期為0.33 h。

4) 從各個(gè)計(jì)算點(diǎn)的振幅來(lái)看, 第一模態(tài)振幅最大, 然后依次遞減, 但位于灣頂?shù)那嗌礁? 蛇蟠水道不僅共振周期有一定的特殊性, 其第二模態(tài)的共振振幅與第一模態(tài)振幅相差較小, 基本呈現(xiàn)雙高峰結(jié)構(gòu), 其他水道則無(wú)此特征, 隨著模態(tài)的增加, 振幅迅速減小。

5) 由于第一模態(tài)為三門灣內(nèi)所有水域普遍具有的共振特征, 而且振幅最為顯著, 因此本文計(jì)算了在第一模態(tài)共振周期下, 整個(gè)海灣的振幅增益如圖6所示。由圖可知, 共振振幅由灣外向三門灣內(nèi)部逐漸增加, 假定開邊界處振幅為1.0 m, 在灣頂處振幅達(dá)到了3.5~4.0 m, 石浦水道共振振幅為1.5~2.0 m, 健跳水道和白礁水道振幅為2.5~3.0 m。通過(guò)本研究可知如果三門灣外海產(chǎn)生的擾動(dòng)周期正好與海灣的共振周期一致, 即使外海的波動(dòng)振幅較小, 由于較強(qiáng)的振幅增益, 也會(huì)對(duì)三門灣內(nèi)沿岸造成嚴(yán)重的災(zāi)害。

圖6 第一模態(tài)共振周期對(duì)應(yīng)的振幅增益(單位: m)

3 結(jié)語(yǔ)

本文分別采用理論推導(dǎo)和數(shù)值方法研究了三門灣的共振特征, 通過(guò)建立高分辨率的數(shù)值模型, 能夠?qū)θT灣內(nèi)各水道的共振周期進(jìn)行精細(xì)化研究, 數(shù)值模擬結(jié)果比理論解更加準(zhǔn)確。研究發(fā)現(xiàn)三門灣內(nèi)各水域普遍存在4~5個(gè)共振周期模態(tài), 各水道第一模態(tài)和第三模態(tài)共振周期與理論值比較接近, 石浦水道由于與外海相聯(lián)通, 并非一個(gè)完整的半封閉水道, 因此共振周期不是非常顯著。三門灣內(nèi)各水域共振周期第一模態(tài)對(duì)應(yīng)的振幅最大, 然后依次遞減, 但灣頂?shù)那嗌礁? 蛇蟠水道, 其第二模態(tài)共振周期及對(duì)應(yīng)的振幅與第一模態(tài)相差較小, 基本呈現(xiàn)雙峰結(jié)構(gòu)。通過(guò)計(jì)算在第一模態(tài)共振周期下整個(gè)海灣的振幅增益可知, 共振振幅由灣外向三門灣內(nèi)部逐漸增加, 尤其是灣頂處振幅會(huì)有顯著的增加, 因此如果外海的擾動(dòng), 正好與海灣的共振周期一致, 由于較強(qiáng)的振幅增益, 可能會(huì)在三門灣內(nèi)造成嚴(yán)重的災(zāi)害, 因此下一步需要結(jié)合實(shí)際觀測(cè)資料進(jìn)行更加細(xì)致深入的研究。

[1] Proudman J. Oscillations of tide and surge in an estuary of finite length[J]. Journal of Fluid Mechanics, 1957, 2(04): 371-382.

[2] Cui X, Fang G, Teng F, et al. Estimating peak response frequencies in a tidal band in the seas adjacent to China with a numerical model[J]. Acta Oceanology Sinica, 2015, 34(1): 29-37.

[3] Cummins P F, Karsten R H, Arbic B K. The semi-diu-rnal tide in Hudson strait as a resonant channel oscillation[J]. Atmosphere-ocean, 2010, 48(3): 163-176.

[4] 王鐘桾, 龍寶森. 半封閉海灣固有振動(dòng)周期的數(shù)值計(jì)算[J]. 海岸工程, 1998, 4: 1-4. Wang Zhongjun, Long Baosen. Computations of natural periods of a half-closed bay using a numerical model[J]. Coastal Engineering, 1998, 4: 1-4.

[5] Webb D J. On the shelf resonances of the Gulf of Carpentaria and the Arafura Sea[J]. Ocean Science, 2012, 8(5): 733.

[6] Crean P B, Murty T S, Stronach J A. Mathematical modelling of tides and estuarine circulation: The coastal seas of southern British Columbia and Washington state[M]. Berlin: Springer-Verlag, 1988: 471-472.

[7] Wilson B S. The Encyclopedia of oceanography[M]. New York: Academic Press, 1966: 804-817.

[8] Wong K C, Moses Hall J E. On the relative importance of the remote and local wind effects to the sub tidal variability in a coastal plain estuary, Journal of Geophysical Research: Oceans, 1998, 103(C9): 18393- 18404.

[9] Clarke A J, Battisti D S. The effect of continental shelves on tides[J]. Deep Sea Research Part A. Oceano-graphic Research Papers, 1981, 28(7): 665-682.

[10] Bertin X, Bruneau N, Breilh J F, et al. Importance of wave age and resonance in storm surges: The case Xynthia, Bay of Biscay[J]. Ocean Modelling, 2012, 42: 16-30.

[11] 崔欣梅, 方國(guó)洪, 滕飛, 等. 芬迪灣、緬因?yàn)澈褪趥愃篂吵毕舱竦臄?shù)值研究[J]. 海洋科學(xué)進(jìn)展, 2015, 33(1): 11-18. Cui Xinmei, Fang Guohong, Teng Fei, et al.Numerical study on tidal resonances in the bay of Fundy, the Gulf of Maine and the Gulf of Saint Lawrence[J]. Advances in Marine Science, 2015, 33(1): 11-18.

[12] Proudman J. Dynamical Oceanography[M]. London: Methuen, 1953: 85-89.

[13] 劉瀟, 馮秀麗, 劉杰. 港口工程影響下萊州灣西南側(cè)海域水動(dòng)力演化特征[J]. 海洋科學(xué), 2016, 40(3): 138-145. Liu Xiao, FengXiuli, Liu Jie. Hydrodynamic evolution characteristics of southwest Laizhou Bay under the effect of port construction[J]. Marine Sciences, 2016, 40(3): 138-145.

[14] 周廣鎮(zhèn), 馮秀麗, 劉杰, 等. 萊州灣東岸近岸海域規(guī)劃圍填海后沖淤演變預(yù)測(cè)[J]. 海洋科學(xué), 2014, 38(1): 15-19. Zhou Guangzhen, Feng Xiuli, Liu Jie, et al.Prediction of erosion evolution and deposition in the east coast of the Laizhou Bay after the implementation of the coastal planning[J]. Marine Sciences, 2014, 38(1): 15-19.

Calculation and study of resonant period in semi-closed bay

YANG Wan-kang1, YI Xiao-fei2, CHEN Zhong-biao3

(1. Key Laboratory of Engineering Oceanography, the Second Institute of Oceanography, State Oceanic Administration, Hangzhou 310012, China; 2. College of Ocean and Meteorology, Guangdong Ocean University, Zhanjiang 524088, China; 3. College of Marine Science, Nanjing University of Information Science and Technology, Nanjing 210044, China)

The wave resonance phenomenon is significant in Sanmen Bay and leads to many marine disasters. As a typical semi-closed bay, we investigated the resonance period of Sanmen Bay by theoretical formulation and numerical simulation. We found the theoretical formula to be of use in estimating the resonance period of the bay but the numerical simulation results were more accurate and reasonable, because the differences in the resonance characteristics of the bay waters can be described in greater detail. The numerical results show there to be 4–5 resonant period modes in Sanmen Bay, and the numerical simulation solutions of the first and third modal resonance periods are close to the theoretical values. The Shipu waterway is connected with the offshore region, so its resonance period is not significant. The first mode amplitude of the resonance period is the largest and then it decreases with increasing mode, but there is only a slight difference between the first and second modes with respect to the period and amplitude values of the Qingshan port and Shepan waterway, which show a bimodal structure. The resonant amplitude gain gradually increases from outside to the inside of the bay, and is especially significant at the top of bay. Through this study, we have established a scientific basis for Sanmen Bay disaster prevention.

semi-closed bay; resonance; theoretical derivation; numerical simulation

(本文編輯: 李曉燕)

P731.23

A

1000-3096(2017)07-0071-07

10.11759/hykx20170311001

2017-03-11;

2017-05-16

廣東海洋大學(xué)近海海洋變化與災(zāi)害重點(diǎn)實(shí)驗(yàn)室開放基金(GLOD1405)

[Guangdong Key Laboratory for Coastal Ocean Variation and Disaster Prediction open Foundation, No.1405]

楊萬(wàn)康(1987-), 男, 山東青州人, 漢族, 工程師, 主要從事海洋預(yù)報(bào)和數(shù)值模擬工作, 電話: 0571-89715123, Email: yangwankang@126.com

Mar. 11, 2017