南中國海臺風(fēng)浪數(shù)值模擬研究
——以臺風(fēng)“珍珠”為例

趙紅軍,宋志堯,徐福敏,程 晨

(1.河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點實驗室,江蘇 南京 210098;2.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室,江蘇 南京 210097)

南中國海臺風(fēng)浪數(shù)值模擬研究
——以臺風(fēng)“珍珠”為例

趙紅軍1,宋志堯2,徐福敏1,程 晨1

(1.河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點實驗室,江蘇 南京 210098;2.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室,江蘇 南京 210097)

以QSCAT/NCEP混合風(fēng)資料和Myers經(jīng)驗?zāi)Ｐ惋L(fēng)場構(gòu)造臺風(fēng)風(fēng)場,并以之作為驅(qū)動風(fēng)場,建立一個基于第三代海浪模式SWAN的兩重嵌套臺風(fēng)浪數(shù)值模擬模型。以0601號臺風(fēng)珍珠為例,對南中國海至廣東的臺風(fēng)浪進(jìn)行數(shù)值模擬研究。將數(shù)值模擬結(jié)果與臺風(fēng)期間Jason-1衛(wèi)星高度計觀測資料和近岸浮標(biāo)實測資料(波高、波向和波周期)作了較為詳細(xì)地比較,并分析臺風(fēng)浪要素的時空分布。結(jié)果顯示臺風(fēng)浪要素的數(shù)值模擬值與實測值吻合良好,表明SWAN模型能夠較好地再現(xiàn)大洋和近岸臺風(fēng)浪的時間發(fā)展過程和空間分布特征。

QSCAT/NCEP混合風(fēng);SWAN模型;臺風(fēng)浪;Jason-1衛(wèi)星高度計

Abstract:With the driven wind incorporating the Myers typhoon model and the QSCAT/NCEP blended wind,a two-nested domain typhoon wavemodel is developed based on the third generationwavemodel SWAN.Asa case study,typhoon Chanchu(0601)is selected to simulate waves from South China Sea to Guangdong province.Analysisof the numerical typhoonwave distribution in time and space ismade and comparisonsof themodel datawith thoseof observation by the Jason-1 satellite altimeter and the nearshore-locatedwave recorder show good results.It is concluded that the SWANmodel could reasonably reproduce the featuresof typhoonwaves inocean and coastalwaters spacially and temporally.

Key words:QSCAT/NCEP blended wind;SWAN model;typhoon waves;Jason-1 satellite altimeter

南海戰(zhàn)略位置重要,油氣資源豐富,地處熱帶和亞熱帶,面積寬廣,地形復(fù)雜,是西太平洋和印度洋之間的主要通道之一。南海是我國近海臺風(fēng)浪出現(xiàn)頻率最大和最為嚴(yán)重的海區(qū)之一,據(jù)文獻(xiàn)[1]的統(tǒng)計,1966～1990年發(fā)生于該海域有效波高大于6m的臺風(fēng)浪共189次,平均每年出現(xiàn)7.6次。臺風(fēng)伴隨著巨浪會造成重大的海難事故和海洋工程的破壞,所以對重點海域的臺風(fēng)浪進(jìn)行正確預(yù)報,深入研究臺風(fēng)浪的生成和發(fā)展

海浪數(shù)值模擬始于20世紀(jì)50年代Gelci[2]基于二維波譜能量傳播方程建立數(shù)值模式。隨著對各種物理過程描述的不斷深入和參數(shù)化形式的不同,模式經(jīng)歷了第一代到第三代的演變,其中SWAN模型[3]就是第三代數(shù)值海浪模式,它全面地考慮了多種物理過程(包括風(fēng)生浪、白帽損耗、底摩擦損耗、深度變淺引起的波浪破碎、波—流相互作用、非線性波—波相互作用等),已經(jīng)被較為廣泛地應(yīng)用于海浪的數(shù)值模擬。文獻(xiàn)[4-7]分別使用SWAN模式對臺灣島和長江口等附近海域的臺風(fēng)浪進(jìn)行了數(shù)值模擬研究,均取得了較好的模擬效果。

臺風(fēng)風(fēng)場的準(zhǔn)確性對臺風(fēng)浪數(shù)值模擬效果具有顯著地影響。目前采用較多的風(fēng)場或是對稱型經(jīng)驗臺風(fēng)模型,如藤田公式[8]、Myers公式[9]等,或是在經(jīng)驗臺風(fēng)風(fēng)場模型中加入再分析風(fēng)場資料,或是大中尺度氣象數(shù)值模式的輸出風(fēng)場。文獻(xiàn)[10-11]分別對Holand經(jīng)驗?zāi)Ｐ惋L(fēng)場和NCEP再分析風(fēng)場相加構(gòu)造的臺風(fēng)風(fēng)場下臺風(fēng)浪模擬效果進(jìn)行了研究,他們的研究成果均表明再分析風(fēng)場可有效的彌補(bǔ)經(jīng)驗?zāi)Ｐ惋L(fēng)場的背景風(fēng)精度不足的缺陷。通過QSCAT/NCEP混合風(fēng)和Myers經(jīng)驗?zāi)Ｐ惋L(fēng)相加的方法構(gòu)造了臺風(fēng)風(fēng)場,以之作為第三代數(shù)值海浪模式SWAN的驅(qū)動風(fēng)場,波浪模型采用自嵌套的方法,對0601號強(qiáng)臺風(fēng)“珍珠”經(jīng)過中國南海和登陸廣東時的臺風(fēng)浪進(jìn)行了數(shù)值模擬研究。

1 海浪模式和臺風(fēng)風(fēng)場

1.1 海浪模式

第三代海浪譜模型SWAN以二維動譜密度表示隨機(jī)波,動譜密度N(σ,θ)與能譜密度E(σ,θ)之間關(guān)系為N(σ,θ)=E(σ,θ)/σ,在球坐標(biāo)系下 ,動譜平衡方程可表示為:

式中:λ是經(jīng)度,φ是緯度。方程左邊第一項為動譜密度隨時間的變化率,第二和第三項表示動譜密度在地理坐標(biāo)空間x和y方向的傳播,第四項表示動譜密度在相對頻率σ空間的傳播,第五項表示動譜密度在譜分布方向θ空間的傳播;方程右邊的S代表以譜密度表示的源匯項,包括風(fēng)能輸入、白浪、破碎、海底摩擦、非線性波-波相互作用等物理過程。

SWAN模式采用全隱式BSBT的有限差分格式[3],該格式具有良好的數(shù)值穩(wěn)定性,但數(shù)值耗散較大,一般僅適用于較小區(qū)域波浪場的數(shù)值計算。Rogers等[12]在總結(jié)和比選多種數(shù)值方法的基礎(chǔ)上,采用了S&L和SORDUP兩種數(shù)值格式發(fā)展了SWAN模型,減少了模式在較大空間尺度計算時的數(shù)值耗散,使之適用于各種空間尺度波浪數(shù)值模擬[14]。

1.2 臺風(fēng)風(fēng)場

模型風(fēng)場采用Myers對稱型臺風(fēng)風(fēng)場,其氣壓場的分布形式[9]為

式中:p0是臺風(fēng)中心氣壓;p∞是臺風(fēng)外圍氣壓,一般取為1 013.3 hpa;r是計算點至臺風(fēng)中心的距離;r0是最大風(fēng)速半徑,可表示為[15]

式中:φ是緯度,Vf是移行風(fēng)風(fēng)速。

通過梯度風(fēng)和氣壓場的關(guān)系,可得到梯度風(fēng)場如下

式中:f是科氏力參數(shù)(f=2ωsinφ,ω是地球自轉(zhuǎn)角速度),ρa(bǔ)是空氣密度。

移行風(fēng)場采用宮崎正衛(wèi)公式[16],其分布形式為

式中:Vx和Vy分別是臺風(fēng)中心移動速度的正東分量和正北分量。

將式(4)和(5)相加,得到臺風(fēng)的模型風(fēng)場為

式中:c1和c2是訂正系數(shù),β是梯度風(fēng)與海面風(fēng)的夾角,θ是計算點和臺風(fēng)中心的連線與x軸的夾角。

在臺風(fēng)中心附近,經(jīng)驗?zāi)Ｐ涂奢^好地反映臺風(fēng)大風(fēng)區(qū)的風(fēng)場特征,但一般僅限于幾百公里的范圍;在臺風(fēng)外圍,由于同時受到臺風(fēng)和其它天氣系統(tǒng)的影響,風(fēng)場一般與經(jīng)驗?zāi)Ｐ筒顒e較大。文獻(xiàn)[10-11]分別使用NCEP再分析資料描述了臺風(fēng)外圍的風(fēng)場特征,他們的研究均表明NCEP再分析資料可有效彌補(bǔ)經(jīng)驗?zāi)Ｐ惋L(fēng)場在臺風(fēng)外圍的不足,則采用NCAR提供的QSCAT/NCEP混合風(fēng)資料作為臺風(fēng)背景風(fēng)場。QSCAT/NCEP混合風(fēng)是由QSCAT衛(wèi)星散射計觀測資料和NCEP再分析資料在時間和空間上混合得到,其時間分辨率為6 h,與NCEP再分析資料相同;空間分辨率為0.5°×0.5°,高于NCEP再分析風(fēng)場1.875°×1.905°。QSCAT/NCEP混合風(fēng)場和臺風(fēng)經(jīng)驗?zāi)Ｐ偷暮铣煞椒ㄈ缦?/p>

式中:VM是經(jīng)驗?zāi)Ｐ惋L(fēng)場;VQ是背景風(fēng)場;e是權(quán)重系數(shù),保證了兩個風(fēng)場的平滑連接。根據(jù)文獻(xiàn)[15],e=c4/(1+c4),c=r/(10×R)。

2 實例和試驗設(shè)計

以2006年第1號強(qiáng)臺風(fēng)“珍珠”為例,研究其在中國南海和登陸廣東時的臺風(fēng)浪。該次臺風(fēng)于2006年5月13日02時進(jìn)入我國南海時增強(qiáng)為臺風(fēng)并沿14°N向西移動,到15日02時中心氣壓為950 hPa,近中心最大風(fēng)速為45m/s,發(fā)展成強(qiáng)臺風(fēng),隨即以90°的右折穩(wěn)定地向正北方向移動直趨廣東,于18日02時在廣東饒平登陸,登陸時中心氣壓960 hPa,近中心最大速為35 m/s,表1所列為該臺風(fēng)的要素,圖1顯示了該臺風(fēng)過程的移動路徑。由于該臺風(fēng)在中國南海強(qiáng)度大(12級以上風(fēng)力維持約116小時)、范圍廣(7級大風(fēng)半徑范圍為400～540 km),所以粵東沿岸深水海區(qū)從14日08時開始就受到3 m以上大浪的影響,在17日08時臺風(fēng)中心穿過20°N時即被納入到8m以上狂浪的影響范圍。

表1 0601號臺風(fēng)“珍珠”的要素表Tab.1 Time series of parameters for typhoon No.0601

采用SWAN40.72版本對該次臺風(fēng)浪過程進(jìn)行數(shù)值模擬研究。計算區(qū)域分I區(qū)和II區(qū)大小兩個區(qū)域,如圖 1 所示 ,其中 I區(qū)的計算范圍為 105°E～125°E、10°N～26°N;II區(qū)的計算范圍為 112.5°E～116.5°E、20°N～23°N,通過I區(qū)的運行結(jié)果提供嵌套區(qū)域II區(qū)的波譜邊界條件。模型采用球坐標(biāo)系下的非靜態(tài)模式,時間和空間上的離散方法為S&L格式,譜空間的離散為中心差分格式。對于 I區(qū),模型的空間分辨率為4′×4′,時間步長為10 min,計算時段從5月13日02時至5月18日02時,共計121 h;對于II區(qū),模型的空間分辨率為1′×1′,時間步長為3min,計算時段從5月15日02時至5月18日02時,共計73 h。在頻率和方向的二維譜空間分辨率上,I區(qū)和II區(qū)相同,頻率的計算從0.01～0.5,以對數(shù)分布劃分為25個;方向的分段為60個,分辨率為6°。風(fēng)能的輸入考慮線性增長和指數(shù)增長兩部分,其中線性增長采用Caraler和Malanotte-Rizzoli的表達(dá)式[17],指數(shù)成長采用Komen等的研究成果[18];海底摩擦造成的能量損耗采用Collins公式[19],其它參數(shù)的設(shè)置均取模型默認(rèn)值。

3 數(shù)值結(jié)果分析

圖2是2006年5月13日02時至5月18日02時數(shù)值模擬的I區(qū)最大有效波高分布圖。該圖表明受臺風(fēng)系統(tǒng)的影響,臺風(fēng)浪大浪區(qū)發(fā)生在臺風(fēng)中心大風(fēng)區(qū)附近,臺風(fēng)移動路徑右側(cè)的有效波高大于左側(cè),且隨著臺風(fēng)的西進(jìn)和北上,臺風(fēng)浪大浪區(qū)也逐漸西行和北移。臺風(fēng)移動過程中數(shù)值模擬的最大有效波高為15.7m,發(fā)生在中沙群島東北洋面附近。為驗證模式對大洋臺風(fēng)浪的模擬效果,圖3將模擬值與Jason-1衛(wèi)星高度計觀測值進(jìn)行了對比分析。在模式運行時段內(nèi)Jason-1衛(wèi)星處于第160個運行周期,有三條軌道靠近臺風(fēng)的移動路徑,分別是第T114、T127和T153軌道,它們的生成時間分別為15日05:37、15日17:48和16日18:09,軌道位置和軌道生成時刻對應(yīng)的臺風(fēng)中心位置(C114、C127和C153)如圖1所示。在剔除衛(wèi)星數(shù)據(jù)中的無效數(shù)據(jù)和粗差數(shù)據(jù)后,將衛(wèi)星觀測波高與模式計算結(jié)果進(jìn)行比較。圖3(a)～(c)分別顯示了沿軌道T114、T127和T153的有效波高分布。由圖可知模式較好地模擬了這三條軌道上的波高變化趨勢,但大部分觀測點的模擬值都稍小于觀測值,尤其是在距離臺風(fēng)中心較遠(yuǎn)的地方,這可能是由于本文構(gòu)造的臺風(fēng)系統(tǒng)偏弱和背景風(fēng)場偏小的原因所造成的。

圖1 數(shù)值模擬的計算范圍及水深等值線Fig.1 Sketch of the computational domain and water depthdistribution

圖2 數(shù)值模擬I區(qū)最大有效波高等值線Fig.2 Themaximum significantwave height contours in the simulated area I

圖4(a)～(d)是數(shù)值模擬的嵌套區(qū)域II在不同時刻的波浪場和風(fēng)場分布圖。由圖4(c)可知在大洋上臺風(fēng)移動路徑的右側(cè),強(qiáng)風(fēng)的持續(xù)作用使得波浪的傳播方向和風(fēng)向較為一致;而在臺風(fēng)移動路徑的左側(cè),由于風(fēng)浪和涌浪的共同作用,波浪的傳播方向偏離風(fēng)向,并且越向臺風(fēng)外圍波向與風(fēng)向之間的偏角就越大。圖4(a)～(b)表明在臺風(fēng)中心進(jìn)入嵌套區(qū)域 II之前,近岸水域持續(xù)受到SE向波浪的影響,波向與風(fēng)向差別較大,這是由于大洋涌浪傳至近岸后波向主要受控于地形變化所致。圖4(d)顯示在臺風(fēng)登陸后,臺風(fēng)大風(fēng)區(qū)附近的近岸水域產(chǎn)生了與風(fēng)向相一致的離岸風(fēng)浪,但是在臺風(fēng)大風(fēng)區(qū)外圍,地形對波浪傳播的影響使得波向偏離風(fēng)向,有時二者甚至相反。值得注意的是在珠江口水域,因為受到口門附近眾多島嶼掩護(hù)的影響,無論是臺風(fēng)登陸前、登陸時和登陸后數(shù)值模擬的波高都很小,造成這一小波高的原因可能是相對珠江口水域嵌套區(qū)域II具有較大的空間分辨率,亦或是SWAN模式相解耦的方法對繞射的計算只是定性地描述[20],所以筆者建議對珠江口水域臺風(fēng)浪的數(shù)值模擬宜嵌套一個基于緩坡方程的細(xì)網(wǎng)格數(shù)值模擬模型,以更好地考慮地形變化和波浪的繞射作用。

為檢驗?zāi)Ｊ綄杜_風(fēng)浪的數(shù)值模擬效果,圖4將模擬的有效波高、平均波周期和波向與近岸浮標(biāo)的實測值進(jìn)行了比較分析。測波浮標(biāo)位于廣東平海灣 ,地理坐標(biāo)為22°33′56″N,114°43′48″E,如圖4中A點所示 ,該浮標(biāo)記錄了臺風(fēng)“珍珠”期間測波點附近的波浪要素變化。圖5分別顯示了浮標(biāo)所在位置處數(shù)值模擬的有效波高、平均波周期和波向與浮標(biāo)實測值的比較,由圖可知波向的模擬值和浮標(biāo)觀測值吻合很好,有效波高的模擬值略偏小,平均波周期的模擬值略偏大,這可能是由于浮標(biāo)所在位置部分方向受到遮擋,數(shù)值高頻波不能傳至浮標(biāo)所在位置所致。

圖3 數(shù)值模擬的有效波高與Jason-1衛(wèi)星的軌道觀測波高比較Fig.3 Comparisonsof the simulated significantwave heightswith data on Jason-1 tracks

圖4 不同時刻數(shù)值模擬的嵌套II區(qū)波浪場和風(fēng)場分布Fig.4 Time seriesof the numericalwave height vector and wind vector distribution in the simulated area II

圖5 數(shù)值模擬的有效波高、平均波周期和波向與浮標(biāo)實測值的比較Fig.5 Comparisonsof the numerical significantwave heights,themean wave periods and thewave directionswith those observed by thewave recorder

4 結(jié) 語

首先通過QSCAT/NCEP混合風(fēng)和Myers經(jīng)驗?zāi)Ｐ惋L(fēng)相結(jié)合構(gòu)造了臺風(fēng)風(fēng)場,然后基于SWAN模型采用自嵌套的方法對0601號強(qiáng)臺風(fēng)“珍珠”經(jīng)過中國南海和登陸廣東時所引起的臺風(fēng)浪進(jìn)行了數(shù)值模擬研究。為了驗證模式對大洋臺風(fēng)浪和近岸臺風(fēng)浪的模擬效果,將模式計算值分別與Jason-1衛(wèi)星高度計觀測資料和近岸浮標(biāo)實測資料(波高、波向和波周期)作了較為詳細(xì)的比較分析,結(jié)果顯示臺風(fēng)浪要素的數(shù)值模擬值與實測值吻合良好,表明SWAN模型可較好地再現(xiàn)大洋和近岸臺風(fēng)浪的時間發(fā)展過程和空間分布特征。

由于珠江口附近島嶼眾多,地形復(fù)雜,波浪繞射效應(yīng)明顯,所以筆者建議在珠江口水域嵌套一個基于緩坡方程的細(xì)網(wǎng)格波浪數(shù)值模擬模型,由SWAN為其提供邊界條件,以更合理地考慮地形變化和波浪繞射對臺風(fēng)浪傳播的影響;此外,以QSCAT/NCEP混合風(fēng)描述的臺風(fēng)外圍風(fēng)場可能偏弱,因此進(jìn)一步采用氣象數(shù)值模式研究臺風(fēng)的生成及發(fā)展過程,較為準(zhǔn)確地刻畫臺風(fēng)的風(fēng)場結(jié)構(gòu),并與SWAN模型耦合研究臺風(fēng)浪的生成和發(fā)展機(jī)制,亦是本文的后續(xù)工作。

志謝:QSCAT/NCEP混合風(fēng)場資料和Jason-1衛(wèi)星高度計資料分別來源于National Center for Atmospheric Research(NCAR)Data Support Section(DSS)和NASA Jet Propulsion Laboratory(JPL)Physical Oceanography Distributed Active Archive Center(PODAAC)。

[1] 許富祥,余宙文.中國近海及其鄰近海域災(zāi)害性海浪監(jiān)測和預(yù)報[J].海洋預(yù)報,1998,15(3):63-68.

[2] Gelci R,CazaléH,Vasal J.Prévision de la houle.La méthode des densités spectroangulaires[J].Bull.Inform.ComitéCentral O-ceanogr.,D′Etude C?tes,1957,9:416-435.

[3] Booij N,Ris R C,Holthuijsen H.A third-generation wavemodel for coastal regions,Part I,Model description and validation[J].J.Geophys.Res.,1999,104(C4):7649-7666.

[4] Qu SW,Liau JM,Hsu TW.Simulating typhoonwaves by SWAN wavemodel in coastalwatersof Taiwan[J].Ocean Eng.,2002,29:947-971.

[5] 陳 希,沙文鈺,閔錦忠.臺灣島鄰近海域臺風(fēng)浪的數(shù)值模擬研究[J].海洋預(yù)報,2002,19(4):1-10.

[6] Xu FM,Will P,Zhang JL,et al.Simulation of typhoon-driven waves in the Yangtze estuary withmultiple-nested models[J].China O-cean Eng.,2005,19(4):613-624.

[7] 胡克林,丁平興,朱首賢,等.長江口附近海域臺風(fēng)浪的數(shù)值模擬—以鹿沙臺風(fēng)和森克拉臺風(fēng)為例[J].海洋學(xué)報,2004,26(5):23-33.

[8] Fajita T.Pressure distribution in typhoon[J].Geophys.Mag.,1952,23:437.

[9] Myers V A.Characteristicsof United Stateshurricanespertinent to levee design for Lake Okechobeem[R].FL Hydromet Report32.Government Printing Office,1954.

[10] 張志旭,齊義泉,施 平,等.衛(wèi)星高度計波高資料的同化試驗分析[J].海洋學(xué)報,2003,25(5):21-28.

[11] 蔣小平,鐘 中,張金善,等.臺風(fēng)浪模擬預(yù)報中的風(fēng)場比較研究[J].海洋通報,2007,26(2):11-19.

[12] RogersW E,Kaihatu JM,Petit H A H,et al.Diffusion reduction in an arbitrary scale third generation wind wavemodel[J].Ocean Eng.,2002,29:1357-1390.

[13] RogersW E,Hwang PA,Wang DW.Investigation of wave growth and decay in the SWAN model:Three regional-scale application[J].J.Phys.Oceanogr.,2003,33:366-389.

[14] Graham H E,Nunn D E.Meteorological conditions pertinent to standard project hurricane[R].Atlantic and Gulf Coasts of United States,National Hurricane Research Project,Report No.3,USWeather Service,1959.

[15] CarrⅢL E,Elsberr YL.Modelsof tropical cyclonewind distribution and beta-effectpropagation for application to the tropical cyclone track forecasting[J].Mon.Wea.Rev.,1997,125:3190-3209.

[16] 宮崎正衛(wèi).海洋物理Ⅲ[M].臺中:東海大學(xué)出版社,1977:311-312.

[17] Cavaleri L,Malanotte-Rizzoli P.Wind wave prediction in shallow water:theory and application[J].J.Geophys.Res.,1981,86(C11):10961-10973.

[18] Komen GJ,Hasselmann S,Hasselmann K.On the exustence of a fully developed wind-sea spectrum[J].J.Phys.Oceanogr.,1984,14:1271-1285.

[19] Collins J I.Prediction of shallow water spectra[J].J.Geophys.Res.,1972,77(C15):2693-2707.

[20] Holthuijsen L H,Herman A,Booij N.Phase-decoupled refraction-diffraction for spectralwavemodels[J].Coastal Eng.,2004,49:291-305.

Numerical simulation of typhoon waves in the South China Sea——a case study of typhoon Chanchu

ZHAO Hong-jun1,SONG Zhi-yao2,XU Fu-min1,CHENGChen1
(1.Key Laboratory of Coastal Disastersand Defence,Ministry of Education,HohaiUniversity,Nanjing 210098,China;2.Key Laboratory of Virtual Geographic Environment,Ministry of Education,Nanjing Normal University,Nanjing 210097,China)

P731.22

A

2009-10-23

海岸災(zāi)害及防護(hù)教育部重點實驗室開放研究基金資助項目;水利部公益性行業(yè)科研專項資助項目(200901032);國家自然科學(xué)基金資助項目(50779015);廣東省交通廳科研專項資助項目(2007-31)

趙紅軍(1980-),男,天津薊縣人,講師,主要從事河口海岸動力學(xué)和水波動力學(xué)研究。E-mail:hjzhao998@yahoo.com機(jī)制,具有重要的現(xiàn)實意義。

1005-9865(2010)03-0128-07