陳海軍，姜紹材，王飛

（1.遼寧省海洋環(huán)境預(yù)報總站,遼寧沈陽110001；2.國家海洋局北海海洋環(huán)境監(jiān)測中心站，廣西北海536000；3.遼寧省海洋與漁業(yè)廳信息中心,遼寧沈陽 110001）

ADCIRC模式在渤海M2分潮模擬中的應(yīng)用研究

陳海軍1，姜紹材2，王飛3

（1.遼寧省海洋環(huán)境預(yù)報總站,遼寧沈陽110001；2.國家海洋局北海海洋環(huán)境監(jiān)測中心站，廣西北海536000；3.遼寧省海洋與漁業(yè)廳信息中心,遼寧沈陽 110001）

利用有限元方法的ADCIRC（Advanced Circulation Model）海洋模式，建立了渤海高分辨率的二維潮汐潮流模型，模式結(jié)果與實測資料吻合良好。模式成功模擬出了M2分潮在渤海的2個無潮點(diǎn)和3個圓流點(diǎn)，位置與前人的研究結(jié)果基本一致。M2分潮流在渤海中央為順時針旋轉(zhuǎn)的旋轉(zhuǎn)流，在遼東灣、渤海灣和萊州灣基本為往復(fù)流。M2分潮的潮汐和潮流的振幅都表現(xiàn)為遼東灣最大，萊州灣最小。M2分潮的潮致余流主要表現(xiàn)在近岸和島嶼附近，最大可達(dá)10 cm/s。

ADCIRC；渤海；M2分潮

1 引言

隨著計算機(jī)的發(fā)展，數(shù)值模擬在海洋科學(xué)的研究中發(fā)揮著越來越重要的作用。有結(jié)構(gòu)網(wǎng)格的數(shù)值模式在模擬岸線復(fù)雜的區(qū)域時會有很大限制，非結(jié)構(gòu)網(wǎng)格的海洋模式由于其能很好的擬合岸線，并且能夠在關(guān)心的區(qū)域靈活加密而越來越多的被用于近海和河口的數(shù)值模擬，因此，本文采用計算網(wǎng)格為三角形的ADCIRC[1]（Advanced Circulation Model）海洋模式來研究渤海M2分潮的傳播特征。

ADCIRC（Advanced Circulation Model）是美國北卡羅來納州大學(xué)開發(fā)的有限元海洋模式，該模式可以采用笛卡爾坐標(biāo)，也可以采用球坐標(biāo)，有二維和三維兩種運(yùn)行方式。該模式采用三角形網(wǎng)格，易于擬合邊界，可以對岸線復(fù)雜的海域進(jìn)行高分辨率的數(shù)值模擬。目前國內(nèi)利用ADCIRC模式進(jìn)行的研究工作還不是很多，比較有代表的是夏波[2]利用ADCIRC在渤海進(jìn)行了風(fēng)暴潮的數(shù)值模擬，曹永華等[3]利用該模式研究了黃驊港的二維流場。

專門針對渤海的潮汐潮流數(shù)值模擬，前人已經(jīng)做了很多工作，如黃祖珂[4]利用根據(jù)二維非線性潮波微分方程模擬了渤海幾個主要分潮的潮汐余流；呂咸青等[5]利用伴隨模式模擬了渤海的M2分潮。本文利用ADCIRC模式，對渤海的M2潮波進(jìn)行高分辨率的模擬，并且比較系統(tǒng)地研究了M2潮波在渤海的各種傳播特征，對于進(jìn)一步認(rèn)識M2潮波在渤海的傳播情況具有重要意義。

2 模式設(shè)置

本文的計算區(qū)域為渤海海域，該區(qū)域的水深場見圖1，由圖1可以看出，地形變化劇烈的地方都集中在近岸和島嶼附近，因此本文在形成計算網(wǎng)格時對近岸和島嶼附近進(jìn)行了加密，模式的分辨率由渤海中央的5'逐漸過渡到近岸的1'。最后采用的計算網(wǎng)格見圖2，本計算網(wǎng)格一共包含了三角形網(wǎng)格節(jié)點(diǎn)10100個，三角形單元18742個。本文采用的是二維球坐標(biāo)系，邊界上采用M2分潮的調(diào)和常數(shù)預(yù)報的水位來作為驅(qū)動，邊界上M2分潮的調(diào)和常數(shù)來自NAO99潮汐模型[6]，該模型將5年的T/P衛(wèi)星海面高度計資料和219個沿岸的潮汐觀測資料同化到水動力模型中，其模型區(qū)域分為全球和日本周邊兩部分，我們采用的是NAO99日本周邊區(qū)域的潮汐模型結(jié)果。全場底摩擦系數(shù)均取為0.0014，模式時間步長為60 s，一共運(yùn)行7天，采用后4天的數(shù)據(jù)進(jìn)行調(diào)和分析，得到了全場M2分潮的潮汐潮流調(diào)和常數(shù)。

圖1 渤海水深和觀測站點(diǎn)分布

圖2 模式計算網(wǎng)格

3 模式驗證

為了驗證模式的計算結(jié)果，我們搜集了渤海沿岸16個驗潮站（見圖1）的調(diào)和常數(shù)，與我們的模式結(jié)果進(jìn)行對比，對比結(jié)果見表1。

從表1可以看出，模式結(jié)果與實測結(jié)果吻合良好，可以用來進(jìn)一步分析M2分潮在渤海的傳播情況。誤差的主要來源可能有以下幾種情況：第一，岸線以及地形的精度不夠高，尤其是在近岸區(qū)域，可能會引起比較大的誤差；第二，邊界條件上采用的M2分潮的調(diào)和常數(shù)的誤差，也會影響計算結(jié)果；第三，全場采用一致的底摩擦系數(shù)，不是很合理，也會造成計算結(jié)果的誤差。以后的工作會從以上三方面入手，同時搜集更多海上觀測數(shù)據(jù)，與模型進(jìn)行對比，將結(jié)果模擬的更準(zhǔn)確。

表1 調(diào)和常數(shù)的計算值與實測值的比較

4 模式結(jié)果分析

4.1 M2分潮潮汐結(jié)果分析

根據(jù)計算結(jié)果，進(jìn)行調(diào)和分析，得到了渤海的M2分潮的同潮圖（見圖3），由圖3可以看出，模式成功模擬出了M2分潮在渤海的兩個無潮點(diǎn)，一個位于秦皇島外海，一個位于老黃河口附近，M2分潮波繞著這兩個無潮點(diǎn)逆時針旋轉(zhuǎn)。本文模擬出的無潮點(diǎn)位置與前人的計算結(jié)果對比見表2：

表2 渤海M2分潮無潮點(diǎn)位置

由表2可以看出，本文得出的M2分潮在渤海的無潮點(diǎn)位置與方國洪的結(jié)果最為接近。同時由圖3可以看出，M2分潮的等振幅線分別由渤海灣、遼東灣和萊州灣的灣口向灣頂逐漸增加。相比較而言，遼東灣的振幅最大，渤海灣次之，萊州灣最小。

圖3 M2分潮同潮圖

4.2 M2分潮潮流結(jié)果分析

在本節(jié)，文章將分別從M2分潮的潮流橢圓要素、圓流點(diǎn)及潮致余流等方面來討論M2分潮在渤海的潮流特征。

4.2.1 M2分潮在渤海的潮流橢圓要素

M2分潮在渤海的潮流橢圓及潮流的旋轉(zhuǎn)率分布分別見圖4和圖5。由圖5可知，M2分潮潮流的旋轉(zhuǎn)率在渤海中央最大，可達(dá)-0.8，同樣，在圖4中渤海中央的潮流橢圓已經(jīng)很接近一個圓，因此，M2分潮流在渤海中央為順時針旋轉(zhuǎn)的旋轉(zhuǎn)流。圖5中，遼東灣和渤海灣的灣口各有一條旋轉(zhuǎn)率為0的等值線，除了遼東灣和渤海灣的灣頂，整個遼東灣和渤海灣的旋轉(zhuǎn)率幾乎為0，同樣，圖4中遼東灣和渤海灣的潮流橢圓圖也幾乎為直線，因此在遼東灣和渤海灣，M2分潮流主要為往復(fù)流，方向大體為：渤海灣為東-西向，遼東灣為西南-東北向。萊州灣內(nèi)龍口市西側(cè)有一旋轉(zhuǎn)率極大值大約在-0.8附近，潮流方向按順時針旋轉(zhuǎn)，萊州灣的其他區(qū)域旋轉(zhuǎn)率比較小，大體在-0.2—0.2之間，基本為往復(fù)流。M2分潮的最大可能潮流等振幅線見圖6，最大可能潮流等振幅線的變化趨勢與潮汐的等振幅線變化趨勢相反，表現(xiàn)為由遼東灣、渤海灣和萊州灣的灣口向灣頂逐漸減小。相比較而言，遼東灣的M2分潮流最大，渤海灣次之，萊州灣最小。

圖4 潮流橢圓分布圖

圖5 M2分潮旋轉(zhuǎn)率分布圖 圖6 M2分潮最大可能潮流振幅分布

圖7 M2分潮最大潮流同潮時線圖 圖8 潮汐余流場分布

4.2.2 M2分潮在渤海的圓流點(diǎn)

圖7是M2分潮的最大潮流同潮時圖，由圖7可以看出，模式成功模擬出了M2分潮在渤海的三個圓流點(diǎn)。同潮流時線繞圓流點(diǎn)旋轉(zhuǎn)，旋轉(zhuǎn)方向為逆時針。本文模擬出的無潮點(diǎn)位置與前人的計算結(jié)果對比見表3。

由表3可以看出，在秦皇島附近的圓流點(diǎn)位置，本文計算結(jié)果同F(xiàn)ANG、趙保仁等和萬振文等的結(jié)果比較一致，在萊州灣口的圓流點(diǎn)位置，本文更接近FANG的結(jié)果。在山東半島北部，對于圓流點(diǎn)的個數(shù)，不同的學(xué)者得出的結(jié)論不一樣，本文可能是受計算區(qū)域的限制，只算出了一個圓流點(diǎn)，該位置與趙保仁等的計算結(jié)果中的一個比較接近。此外，通過對比圖7與圖5可以發(fā)現(xiàn)，圓流點(diǎn)的位置，都基本對應(yīng)著旋轉(zhuǎn)率的極大值。

4.2.3 M2分潮在渤海的潮致余流

潮致余流通常是由于潮流的非線性項受到側(cè)向岸線和底摩擦的作用，一部分周期性能量變成了非周期的能量。本文對模式輸出的一個M2周期的潮流數(shù)據(jù)進(jìn)行平均，剔除周期性信號，得到了渤海的M2分潮的潮致余流分布，如圖8所示。由圖8可以看出，M2分潮的余流在渤海普遍比較小，余流主要在近岸和島嶼附近比較明顯，最大值大約在10 cm/s左右。渤海中央有一個順時針旋轉(zhuǎn)的余流，但是量級比較小，大約在2 cm/s，這與趙保仁[4]的結(jié)果比較一致。

表3 渤海M2分潮圓流點(diǎn)位置

5 結(jié)論

ADCIRC的模擬和實測結(jié)果的對比表明，該模式可以很好地刻畫渤海的潮汐潮流現(xiàn)象。模式成功地模擬出了M2分潮在渤海的兩個無潮點(diǎn)和3個圓流點(diǎn)，位置與前人的結(jié)果基本一致。對模式結(jié)果分析發(fā)現(xiàn)：M2分潮的振幅在遼東灣最大，渤海灣次之，萊州灣最小，分布規(guī)律都是由灣口向灣頂逐漸增大；渤海的中央為順時針旋轉(zhuǎn)流，旋轉(zhuǎn)率可達(dá)-0.8，遼東灣、渤海灣和萊州灣內(nèi)主要為往復(fù)流，旋轉(zhuǎn)率在0附近；最大可能潮流等振幅線的變化趨勢與潮汐的等振幅線變化趨勢相反，表現(xiàn)為由遼東灣、渤海灣和萊州灣的灣口向灣頂逐漸減小，潮流的強(qiáng)弱同樣是遼東灣最強(qiáng)，萊州灣最弱。M2分潮的潮致余流主要在近岸比較明顯，最大可達(dá)10 cm/s。

[1]Luettich R,Westerrink J.ADCIRC usermanual:a(parallel)advanced circulation modelfor oceanic[R].Coastaland Estuarine Water,2006.

[2]夏波.風(fēng)暴潮過程中的波流耦合數(shù)值模式研究[D].天津大學(xué),2005.

[3]曹永華.黃驊港海域二維流場分析[J].中國港灣建設(shè),2005,3:1-4.

[4]黃祖珂.渤海的潮汐余流[J].海洋湖沼通報,1992,3:1-8.

[5]呂咸青,方國洪.渤海M2分潮的伴隨模式數(shù)值實驗[J].海洋學(xué)報,2002,24(1):17-24.

[6]Matsumoto K,Takanezawa T,Ooe M.Ocean tide models developed by assimilating TOPEX/POSEIDON altimeter data into hydrodynamical model:A global dodel and a regional model around Janpan[J].Journal of Oceanography,2000,56:567-581.

[7]沈育疆.東中國海潮汐數(shù)值計算[J].山東海洋學(xué)院學(xué)報,1980,10(3):28-35.

[8]Fang G.Tide and tidal current charts for the marginal seas adjacentto China[J].Chinese JournalofOceanology and Limnology,1986,4(1):1-16.

[9]陳國珍.渤、黃、東海海洋圖集（水文分冊）[M].北京:海洋出版社,1992.

[10]趙保仁,方國洪,曹德明.渤、黃、東海潮汐潮流的數(shù)值模擬[J].海洋學(xué)報,1994,16(5):1-10.

[11]萬振文,喬方利,袁業(yè)立.渤、黃、東海三維潮波運(yùn)動數(shù)值模擬[J].海洋與湖沼,1998,29(6):611-616.

[12] 王凱,方國洪,馮士筰.渤海、黃海、東海M2潮汐潮流的三維數(shù)值模擬[J].海洋學(xué)報,1999,21(4):1-13.

[13]林琿,閭國年,宋志堯等.東中國海潮波系統(tǒng)與海岸演變模擬研究[M].北京:科學(xué)出版社,2000.

Application ofADCIRC in M2 tide modeling in the Bohai Sea

CHEN Hai-jun1，JIANG Shao-cai2，WANG Fei3
(1.Marine Environmental Forecasting Central Station of Liaoning Province,Shenyang 110001,China;2.Beihai Environmental Monitoring Station of State Oceanic Administration,Beihai 536000,China;3.Ocean and Fishery Department of Liaoning Province Information Center,Shenyang 110001,China)

With the application of ADCIRC,a high resolution of two-dimensional tidal model in the Bohai Sea was built up,and the model results are conformable to the field observations.Two amphidromic points and three tidal current non-phase points for M2 tide were successfully modeled,whose positions were basically conformable to the previous studies.The style of the M2 tidal current is rotary current in the middle of the Bohai Sea,while rectilinear current in the Liaodong,Bohai and Laizhou Bay.The amplitude of the M2 tide and tidal current are largest in the Liaodong Bay,while smallest in the Laizhou Bay.The M2 tidal reduced residual current is only evident along the coastal region and at the area near the islands,which can reach as much as 10 cm/s.

ADCIRC;Bohai Sea;M2 tide

P731

A

1003-0239（2011）04-0070-06

2010-12-07

陳海軍（1982-），男，助理工程師，從事海洋環(huán)境預(yù)警報和海洋防災(zāi)減災(zāi)工作。E-mail:gh476007@163.com