国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

渤黃海潮汐潮流精細(xì)化數(shù)值模擬及可視化預(yù)報(bào)系統(tǒng)

2011-12-23 08:44:22馮興如尹寶樹(shù)楊德周
海洋預(yù)報(bào) 2011年4期
關(guān)鍵詞:黃海潮汐調(diào)和

馮興如 ,尹寶樹(shù) ,楊德周

(1.中國(guó)科學(xué)院 海洋研究所,山東青島266071;2.中國(guó)科學(xué)院研究生院,北京100049;3.中國(guó)科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室,山東青島 266071)

渤黃海潮汐潮流精細(xì)化數(shù)值模擬及可視化預(yù)報(bào)系統(tǒng)

馮興如1,2,3,尹寶樹(shù)1,3,楊德周1,3

(1.中國(guó)科學(xué)院 海洋研究所,山東青島266071;2.中國(guó)科學(xué)院研究生院,北京100049;3.中國(guó)科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室,山東青島 266071)

利用有限元海洋模式ADCIRC(Advanced Circulation Model),建立了高分辨率的渤黃海二維潮汐潮流數(shù)值模型,該模型以M2,S2,K1,O1等8個(gè)分潮的水位作為驅(qū)動(dòng),模擬出了該8個(gè)分潮的潮汐潮流調(diào)和常數(shù);利用該調(diào)和常數(shù)預(yù)報(bào)的潮位和二維平均潮流與實(shí)測(cè)資料相比,符合較好;利用模擬得到的潮汐潮流調(diào)和常數(shù),以Matlab為平臺(tái),建立了能獨(dú)立運(yùn)行的可視化潮汐潮流預(yù)報(bào)系統(tǒng),該系統(tǒng)可以預(yù)報(bào)渤黃海任意時(shí)刻、任意地點(diǎn)的潮位和平均潮流,并且可以查詢渤黃海任意地點(diǎn)的部分潮汐特征值,該系統(tǒng)運(yùn)行快捷,使用方便,可為海上作業(yè)提供環(huán)境預(yù)報(bào)服務(wù)。

ADCIRC模式;潮汐潮流預(yù)報(bào)系統(tǒng);渤黃海

1 引言

潮汐潮流對(duì)人類在近岸的活動(dòng)有著很大的影響,因此潮汐潮流的預(yù)報(bào)是一個(gè)海洋研究中重要而久遠(yuǎn)的話題。傳統(tǒng)的方法一般采用潮汐潮流圖表,這種方法只能預(yù)報(bào)有限站位,使用不方便,且不適合外海使用,因此開(kāi)展區(qū)域海洋的微機(jī)化潮汐潮流預(yù)報(bào)系統(tǒng)就顯得非常有必要。方國(guó)洪等[1]在2005—2006年期間通過(guò)同化高度計(jì)和沿岸站的資料,得到了分辨率為5'×5'的潮汐潮流預(yù)報(bào)產(chǎn)品,該預(yù)報(bào)產(chǎn)品覆蓋了我國(guó)近海各海區(qū),可預(yù)報(bào)任意時(shí)刻和地點(diǎn)的潮汐潮流;中國(guó)海洋大學(xué)開(kāi)發(fā)的Chinatide[2]潮汐預(yù)報(bào)軟件,可以對(duì)中國(guó)近海各海區(qū)和日本海的一部分進(jìn)行潮汐預(yù)報(bào),其分辨率也是5'×5'。但是近岸地形和岸線復(fù)雜,5'×5'的分辨率有時(shí)不能滿足要求,因此在近岸開(kāi)展更高分辨率的潮汐潮流數(shù)值模擬及預(yù)報(bào)系統(tǒng)的開(kāi)發(fā)就非常有意義。

非結(jié)構(gòu)網(wǎng)格的海洋模式能夠?qū)Π毒€復(fù)雜,地形變化劇烈的海域進(jìn)行高分辨率的模擬,因此其越來(lái)越多地被用來(lái)進(jìn)行近岸和河口的數(shù)值研究。例如,王培濤等[3]利用基于三角形網(wǎng)格的ELCIRC(3D Eulerian-Lagrangian Circulation)海洋模型,對(duì)黃渤海的溫帶風(fēng)暴潮進(jìn)行了高分辨率的數(shù)值模擬;馮興如等[4]利用有限體積海洋模式FVCOM(Finite Volume Coast and Ocean Model)對(duì)龍口海域的潮汐潮流進(jìn)行了精細(xì)的數(shù)值研究。本文首先利用非結(jié)構(gòu)網(wǎng)格的有限元海洋模式ADCIRC(Advanced Circulation Model)[5]模擬得到渤黃海近岸高分辨率的潮汐潮流調(diào)和常數(shù),然后基于這些調(diào)和常數(shù),以matlab為平臺(tái)開(kāi)發(fā)了一套可視化的潮汐潮流預(yù)報(bào)系統(tǒng)。

2 ADCIRC模式簡(jiǎn)介及設(shè)置

ADCIRC是美國(guó)北卡羅來(lái)納州大學(xué)開(kāi)發(fā)的有限元海洋模式,該模式可以采用笛卡爾坐標(biāo),也可以采用球坐標(biāo),有二維和三維兩種運(yùn)行方式。該模式采用三角形網(wǎng)格,易于擬合邊界,可以對(duì)岸線復(fù)雜的海域進(jìn)行高分辨率的數(shù)值模擬。二維ADCIRC模型采用沿水深積分的時(shí)均連續(xù)方程和運(yùn)動(dòng)方程,其表達(dá)式為:

其中,t為時(shí)間,x和y分別為x軸和y軸方向,ζ為從平均海平面起算的水位高度,U和V分別表示x和y方向垂向平均流速,H為總水深,f為科氏力系數(shù),Ps為表面大氣壓力,ρ0為水密度,g為重力加速度,(η+γ)表示牛頓潮勢(shì)和固體潮作用,τsx和τsy表示表面風(fēng)應(yīng)力和波浪輻射應(yīng)力的x和y向分量,τbx和τby、Dx和 Dy、Bx和 By分別代表底部切應(yīng)力、擴(kuò)散項(xiàng)以及斜壓梯度的x和y向分量。

本文采用的是二維球坐標(biāo)系,模式的分辨率由渤黃海海中央的5'逐漸過(guò)渡到近岸的1'(見(jiàn)圖1),一共61357個(gè)節(jié)點(diǎn),117265個(gè)三角形網(wǎng)格。邊界上采用M2,S2,K1,O1,N2,S2,P1,Q1共8個(gè)分潮的調(diào)和常數(shù)預(yù)報(bào)的水位作為驅(qū)動(dòng),邊界上的調(diào)和常數(shù)來(lái)自NAO99潮汐模型[6]。 模式時(shí)間步長(zhǎng)為10 s。

3 模式結(jié)果與驗(yàn)證

模式一共運(yùn)行50天,采用后40天的數(shù)據(jù)進(jìn)行調(diào)和分析。得到了和邊界驅(qū)動(dòng)對(duì)應(yīng)的8個(gè)分潮的潮汐潮流調(diào)和常數(shù)。由于篇幅限制,本文只給出其中M2和K1分潮的同潮圖(見(jiàn)圖2),從圖2可以看出,該模式成功模擬出了半日分潮M2在渤黃海的4個(gè)無(wú)潮點(diǎn)和全日分潮K1在渤黃海的2個(gè)無(wú)潮點(diǎn),模式模擬的M2和K1的同潮圖,同《渤海黃海東海海洋圖集》[7](水文)里相應(yīng)的同潮圖基本一致。

圖1 模式計(jì)算網(wǎng)格

圖2 M2,K1分潮同潮圖,虛線為等振幅線,實(shí)線為等遲角線(格林威治遲角)

圖3 觀測(cè)點(diǎn)位置示意圖

圖4 水位模擬結(jié)果與實(shí)測(cè)結(jié)果對(duì)比圖

圖5 流速模擬結(jié)果與實(shí)測(cè)結(jié)果對(duì)比圖

為了進(jìn)一步驗(yàn)證模式結(jié)果,我們搜集了渤黃海四個(gè)觀測(cè)點(diǎn)的一些實(shí)測(cè)的水位和流速資料,將這四個(gè)觀測(cè)點(diǎn)分別標(biāo)記為A,B,C,D,其位置見(jiàn)圖3。這四個(gè)觀測(cè)點(diǎn)的資料情況是:A,D兩點(diǎn)僅有水位資料,B點(diǎn)僅有流速資料,C點(diǎn)既有流速又有水位資料。利用模擬的調(diào)和常數(shù)預(yù)報(bào)的潮位和平均潮流與實(shí)測(cè)資料的對(duì)比情況分別見(jiàn)圖4和圖5。

四個(gè)觀測(cè)點(diǎn)的對(duì)比結(jié)果顯示,調(diào)和常數(shù)預(yù)報(bào)的水位和流速與實(shí)測(cè)值的位相基本重合,水位絕對(duì)誤差的平均值為15 cm,而利用NAO99資料預(yù)報(bào)的水位與上述實(shí)測(cè)水位的絕對(duì)誤差平均值為17.8 cm,因此本文得到的調(diào)和常數(shù)在中國(guó)的近岸能更好一些;流速東分量絕對(duì)誤差的平均值8.95 cm/s,流速北分量絕對(duì)誤差的平均值為9.22 cm/s。我們認(rèn)為,這些誤差都在可以接受的范圍內(nèi)。模擬結(jié)果與實(shí)測(cè)結(jié)果的誤差可能來(lái)源于開(kāi)邊界的的驅(qū)動(dòng)、地形精度和模式參數(shù)的設(shè)置等。

4 可視化潮汐潮流預(yù)報(bào)系統(tǒng)的實(shí)現(xiàn)

為了快速方便地將模擬得到的潮汐潮流調(diào)和常數(shù)用來(lái)預(yù)報(bào)潮位和潮流,我們以Matlab為平臺(tái),開(kāi)發(fā)了能獨(dú)立運(yùn)行,可視化的潮汐潮流預(yù)報(bào)系統(tǒng),其界面見(jiàn)圖6。該系統(tǒng)分為輸入?yún)^(qū)、查詢功能區(qū)和預(yù)報(bào)功能區(qū)。

4.1 系統(tǒng)工作原理

圖6 渤黃海潮汐潮流預(yù)報(bào)系統(tǒng)界面

預(yù)報(bào)功能區(qū)工作原理為:系統(tǒng)根據(jù)用戶在輸入?yún)^(qū)輸入的經(jīng)緯度和時(shí)間信息,自動(dòng)尋找和輸入地點(diǎn)最近的網(wǎng)格點(diǎn),得到該地點(diǎn)的潮汐潮流調(diào)和常數(shù),然后自動(dòng)調(diào)用t_tide[9]函數(shù),預(yù)報(bào)所要求時(shí)刻的潮汐和平均潮流值,并顯示在界面的相應(yīng)位置。

4.2 系統(tǒng)使用方法

輸入?yún)^(qū)的使用:在輸入?yún)^(qū),用戶可以手動(dòng)輸入要預(yù)報(bào)的時(shí)間,經(jīng)緯度。如果想自動(dòng)輸入當(dāng)前時(shí)間,可以點(diǎn)“現(xiàn)在時(shí)間”按鈕。

查詢功能的使用:點(diǎn)擊“查詢”按鈕,系統(tǒng)便會(huì)根據(jù)用戶輸入的經(jīng)緯度查詢?cè)撐恢玫某毕愋?、潮流類型、平均大的潮差或者平均大潮差。?dāng)潮汐類型為“正規(guī)半日潮”時(shí)只顯示“平均大潮差”的結(jié)果,“平均大的潮差”顯示為“無(wú)”,當(dāng)潮汐類型為“正規(guī)半日潮”以外的類型時(shí),只顯示“平均大的潮差”的結(jié)果,“平均大潮差”顯示為“無(wú)”。

預(yù)報(bào)功能的使用:點(diǎn)擊“預(yù)報(bào)”按鈕,系統(tǒng)會(huì)根據(jù)用戶輸入的經(jīng)緯度和時(shí)間來(lái)預(yù)報(bào)該地點(diǎn)輸入時(shí)刻的水位和潮流、流向,并把結(jié)果顯示在相應(yīng)的文本框里,其中流速方向的規(guī)定為:流向北為0°,東為90°,南為180°,西為270°。同時(shí)在界面的左下方,會(huì)顯示輸入時(shí)刻之后12小時(shí)的潮位變化曲線,在該曲線圖里,會(huì)自動(dòng)標(biāo)出輸入時(shí)刻之后的12小時(shí)里,出現(xiàn)最大和最小潮位的時(shí)刻以及對(duì)應(yīng)的水位值。

5 總結(jié)與未來(lái)工作展望

本文開(kāi)發(fā)的潮汐潮流預(yù)報(bào)系統(tǒng)具有近岸高分辨率,運(yùn)行快捷方便等特點(diǎn),可以為海上作業(yè)提供海洋環(huán)境預(yù)報(bào)和查詢服務(wù),未來(lái)的工作主要集中在:(1)搜集更多的實(shí)測(cè)數(shù)據(jù),不斷改善模式模擬的結(jié)果;(2) 重點(diǎn)區(qū)域再加密,在一些重要的港口和城市進(jìn)行分辨率更高的模擬;(3) 將模式從二維發(fā)展到三維,使系統(tǒng)可以預(yù)報(bào)任意時(shí)刻和深度的潮流;(4) 外拓模式的開(kāi)邊界,以保證邊界上調(diào)和常數(shù)的準(zhǔn)確性。

致謝:感謝美國(guó)北卡羅來(lái)納州大學(xué)提供ADCIRC模式源代碼。

[1]方國(guó)洪,魏澤勛,王永剛.我國(guó)潮汐潮流區(qū)域預(yù)報(bào)的發(fā)展[J].地球科學(xué)進(jìn)展,2008,23(4):331-336.

[2]李孟國(guó),鄭敬云.中國(guó)海域潮汐預(yù)報(bào)軟件Chinatide的應(yīng)用[J].水道港口,2007,28(1):65-68.

[3]王培濤,董劍希,趙聯(lián)大等.黃渤海精細(xì)化溫帶風(fēng)暴潮數(shù)值預(yù)報(bào)模式研究及應(yīng)用[J].海洋預(yù)報(bào),2010,27(4):1-8.

[4] 馮興如,楊德周,尹寶樹(shù).FVCOM在龍口海域潮汐潮流模擬中的應(yīng)用研究[J].海洋科學(xué),2010,34(6):94-99.

[5]Luettich R A,Hu J S,Westerrink J J.The development of the direct stress solution technique for three-dimensional hydrodynamic models using finite elements[J].International Journal for Numerical Methods in Fluids,1994,(19):295-319.

[6]Matsumoto K,Takanezawa T,Ooe M.Ocean tide Models Developed by Assimilating TOPEX/POSEIDON Altimeter Data into Hydrodynamical Model:A Global Model and a Regional Model around Janpan[J].Journal of Oceanography,2000,56:567-581.

[7]海洋圖集編委會(huì).渤海黃海東海海洋圖集(水文)[M].北京:海洋出版社,1993.429-432.

[8]方國(guó)洪,鄭文振,陳宗鏞等.潮汐和潮流的分析和預(yù)報(bào)[M].北京:海洋出版社,1986.216-233.

[9]Pawlowicz R,Beardsley B,Lentz S.Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE[J].Computers and Geosciences,2002,28:929-937.

Refined simulation and establishment of a tide and tidal current forecasting system in the Bohai Sea and the Yellow Sea

FENG Xin-gru1,2,3,YIN Bao-shu1,3,YANG De-zhou1,3

(1.Institute of Oceanology,the Chinese Academy of Sciences,Qingdao 266071 China;2.Graduate School,the Chinese Academy of Sciences,Beijing 100049 China;3.Key Laboratory of Ocean Circulation and Wave,the Chinese Academy of Sciences,Qingdao 266071 China)

A refined two-dimensional tide and tidal current model is built based on the finite element ocean model ADCIRC(Advanced Circulation Model).The model is forced by 8 tide constituents of M2,S2,K1,O1,etc,and the corresponding 8 tide harmonic constituents are calculated using the model result.The tide and tidal current forecasted by the calculated harmonic constituents agrees well with the observation.A visual tide and tidal current forecasting system which can be operated independently is established based on Matlab.The system can forecast the tide and depth-integrated tidal current in the Bohai Sea and the Yellow Sea.Some of the tide characteristics also can be shown through the system.The tide and tidal current forecasting system can be operated fast and conveniently and provide forecasting service for the offshore operation.

ADCIRC model;tide and tidal current forecasting system;Bohai Sea and Yellow Sea

P731

A

1003-0239(2011)04-0065-05

2010-11-01

中科院創(chuàng)新項(xiàng)目(KZCX2-YW-BR-215-3);中科院創(chuàng)新項(xiàng)目群(KZCX2-YWQ07-01);國(guó)家“908”補(bǔ)充調(diào)查項(xiàng)目(908-01-BC12);“908”“風(fēng)暴潮災(zāi)害對(duì)沿海地區(qū)社會(huì)經(jīng)濟(jì)發(fā)展綜合評(píng)價(jià)”(908-02-03-02)

馮興如(1984-),男,博士研究生,主要從事浪流相互作用研究,E-mail:fengxingru07@qdio.ac.cn

猜你喜歡
黃海潮汐調(diào)和
潮汐與戰(zhàn)爭(zhēng)(上)
五味調(diào)和醋當(dāng)先
東方濕地 黃海明珠
黃海簡(jiǎn)介
絕美海灘
從“調(diào)結(jié)”到“調(diào)和”:打造“人和”調(diào)解品牌
調(diào)和映照的雙Lipschitz性質(zhì)
潮汐式灌溉控制系統(tǒng)的設(shè)計(jì)及應(yīng)用
電子制作(2017年9期)2017-04-17 03:00:56
三角恒等變換去哪兒了
干法紙的潮汐
生活用紙(2016年6期)2017-01-19 07:36:25
通榆县| 阜南县| 丰台区| 砚山县| 南川市| 金乡县| 白山市| 马边| 许昌市| 基隆市| 都昌县| 通海县| 南宫市| 长兴县| 桂林市| 大姚县| 鹤庆县| 高台县| 天长市| 祁门县| 西丰县| 浑源县| 伊宁市| 小金县| 高州市| 霍邱县| 碌曲县| 永新县| 象山县| 海口市| 天长市| 巴塘县| 宁津县| 昭苏县| 永丰县| 石泉县| 宜兴市| 溧阳市| 从江县| 简阳市| 巴林右旗|