基于JASON-2測高交叉點的東海有效波高SWAN模擬分析

謝友鴿，郭金運*，朱金山，劉 新，李國偉

（1.山東科技大學(xué)測繪科學(xué)與工程學(xué)院，山東 青島 266590；2.威海市海洋環(huán)境監(jiān)測中心，山東 威海 264209）

基于JASON-2測高交叉點的東海有效波高SWAN模擬分析

謝友鴿1，郭金運1*，朱金山1，劉 新1，李國偉2

（1.山東科技大學(xué)測繪科學(xué)與工程學(xué)院，山東 青島 266590；2.威海市海洋環(huán)境監(jiān)測中心，山東 威海 264209）

以JASON-2衛(wèi)星在上升和下降軌跡中所形成的交叉點的東海局部海域作為研究區(qū)域，以3個交叉點作為研究的實驗點，SWAN模式基于1'-0.5'（間隔為0.1'）的三角形水深網(wǎng)格的分辨率，對東海局部區(qū)域模擬2011年1-8月份的波浪數(shù)值，并將模擬的結(jié)果與JASON-2衛(wèi)星的有效波高進行對比，結(jié)果表明，采用局部加密的三角形網(wǎng)格能夠提高波浪數(shù)值模擬的精度，在采用有限三角形網(wǎng)格分辨率情況下，模擬精度隨著網(wǎng)格分辨率的規(guī)律性增加呈現(xiàn)不規(guī)則遞增現(xiàn)象，同時，設(shè)備運算效率會出現(xiàn)指數(shù)似的降低。通過對運算結(jié)果的統(tǒng)計可知，采用0.6'的三角形水深網(wǎng)格的分辨率獲得模擬結(jié)果中精度幅度提高最大，同時加密的三角形網(wǎng)格主要對水深梯度較大的東南區(qū)域作用效果明顯。

swan模式；JASON-2衛(wèi)星；有效波高；三角形網(wǎng)格

海洋是高度復(fù)雜的非線性動力系統(tǒng)，具有混沌和不穩(wěn)定的特性[1]。海洋中的海浪不僅是海洋動力學(xué)研究的重點，而且在海氣相互作用等研究領(lǐng)域占有重要位置。海浪的數(shù)值模擬結(jié)合海浪實測一直是研究海浪的主要手段和工具，通過對模擬結(jié)果分析和海浪預(yù)報，可以對真實海浪觀測的不足作補充[2]。

SWAN模式是第三代海浪數(shù)值模擬的主要工具，目前具有穩(wěn)定性高、適應(yīng)性好，模擬精度高等優(yōu)點[3]。從SWAN模式建立之初至今，科學(xué)工作者對SWAN模式進行了大量的研究與改進工作，使得SWAN模式逐漸走向成熟。20世紀90年代末，Ris等[4]模擬了北海南部的德國和荷蘭的Haringvliet、Norderneyer Seegat和Friesche Zeegat水域的波浪場，并對SWAN模式進行了分析評定，結(jié)果表明SWAN模式模擬結(jié)果穩(wěn)定可靠。此后，新西蘭的Richard M.Gorman和 Cameron G.Neilson[5]利用SWAN模式模擬了潮間帶河口淺水區(qū)域的風浪的生成和傳播過程,模擬的結(jié)果與實測的結(jié)果非常接近。采用基于能量平衡方程的SWAN模式，不僅得到了國際上的廣泛認可，在中國也得到了大量的應(yīng)用，大量的模擬實例表明SWAN模式在中國海域的海浪數(shù)值模擬是可靠穩(wěn)定的。一些學(xué)者利用SWAN模式對遼東灣海浪進行了模擬，并與浮標實測數(shù)據(jù)對比分析，表明SWAN模式能夠較好地模擬遼東灣近岸水域海浪，同時也分析了在南風過程中，此區(qū)域內(nèi)底摩擦耗散主要作用的范圍[6]。SWAN模式可以較好地模擬理想情況下和實際風浪情況下中國海域內(nèi)的風浪生成和傳播過程，該模式可以應(yīng)用于中國近岸海域地區(qū)的風浪預(yù)報[7-8]。同時，SWAN模式在開發(fā)海浪波能方面也得到了廣泛的應(yīng)用?；谝燥L場驅(qū)動的第三代海浪數(shù)值模式SWAN模擬多年南海波浪能的分布，并對模擬的波浪能資源進行綜合分析研究，尋找中國南海海域蘊藏著的豐富適宜開發(fā)的波能資源[9-10]。

SWAN模式自發(fā)布以來經(jīng)過多次改進以獲得理想的波浪參數(shù)估計，目前采用多嵌套技術(shù)、提高三角形水深網(wǎng)格分辨率和改進物理機理參數(shù)項等方法以達到提高模擬精度的目的?；诮Y(jié)構(gòu)網(wǎng)格模式WaveWatch3和非結(jié)構(gòu)網(wǎng)格海浪數(shù)值模式SWAN的多重嵌套網(wǎng)格，許多研究者建立了一套適用于海域海浪數(shù)值預(yù)報系統(tǒng)，在采用此嵌套網(wǎng)格的情況下，預(yù)報精度均能達到海浪精細化業(yè)務(wù)化預(yù)報的要求[11-12]。SWAN模式中Phillips線性增長項、白冠耗散項和風能輸入項等都與波浪成長有很大關(guān)系，一些學(xué)者通過對這些影響波浪推算的因素進行了分析探討并在此基礎(chǔ)上做出了相應(yīng)的改進，結(jié)果顯示改進的SWAN模式模擬的海浪數(shù)值更接近真實值[1,13]，這使得國際上通用的SWAN波浪模式在中國海域有更好的適用性。SWAN模式基于相同的三角形網(wǎng)格和矩形網(wǎng)格分辨率的模擬中，通過對比分析結(jié)果表明采用同樣分辨率的三角形網(wǎng)格模擬的效果較好，在此基礎(chǔ)上局部加密的三角形網(wǎng)格使得模擬的精度更好[14-15]。

加密的不規(guī)則三角形網(wǎng)格不僅能較好地刻畫復(fù)雜地形，而且SWAN模式采用此網(wǎng)格格式模擬的有效波高能夠獲得更高的精度和準確度[14]。不規(guī)則三角形網(wǎng)格可以改進網(wǎng)格復(fù)雜不規(guī)則海域地形的空間匹配,特別對水深梯度較大的海域地形,但基于不同分辨率的三角形網(wǎng)格運算耗時不等，網(wǎng)格分辨率的增加使得SWAN模式運行效率降低[15]。考慮到計算機耗時時間、不規(guī)則三角形網(wǎng)格分辨率對模擬精度的影響，選取合適的不規(guī)則三角形網(wǎng)格分辨率是至關(guān)重要的，為了分析不規(guī)則三角形網(wǎng)格分辨率、水深梯度、設(shè)備的運算效率與模擬精度的關(guān)系，本文主要選取了2011年1-8月份的中國東海陸架淺水域中JASON-2衛(wèi)星軌跡交叉點作為模擬的實驗點。由于獲取東?，F(xiàn)場浮標等觀測資料的缺乏，為了說明SWAN模式模擬的有效性，本文將模擬結(jié)果與2011年1-8月份的JASON-2衛(wèi)星觀測的有效波高數(shù)值大小進行對比分析，JASON-2衛(wèi)星在空間上采用臨近兩點線性插值獲得與SWAN模式相適應(yīng)的實驗位置。

1 SWAN控制方程

SWAN模式以二維動譜密度表示隨機波，因在流場中，動譜密度守恒，而能譜密度不守恒，動譜密度N(σ,θ)為能譜密度E(σ,θ)與相對頻率σ之比，在球坐標系下，動譜密度方程[16]

其中笛卡爾坐標系轉(zhuǎn)換為球坐標系的過程中，需要滿NR2cosφdσdθdφdλ=N dσdθd x d y

式中：R表示地球半徑；λ表示經(jīng)度；φ表示緯度。

式中：S為以譜密度表示的源項；Sin代表風輸入項；Sds代表由于底摩擦，白浪，破碎等引起的能量耗損；Snl代表波與波之間相互作用引起的波能量變化。

2 數(shù)值模擬

2.1 資料

主要選取東海部分海域作為研究區(qū)域，區(qū)域范圍26°N～32°N,120°E～126°E,其中水深和海岸線數(shù)據(jù)來源于美國國家地球物理資料中心NGDC (http://www.ngdc.noaa.gov/mgg/global/global.html)提供的ETOPO1,ETOPO1中所有的水深數(shù)據(jù)都是基于WGS-84坐標系，這與SWAN模型采用WGS-84坐標系進行模擬相一致，其中水深數(shù)據(jù)是以網(wǎng)格形式排列，空間分辨率為1'x1'。如圖1所示，東海海域的水深分布由東南向西北方向呈現(xiàn)階梯狀，與該區(qū)域內(nèi)的海岸線分布方向較匹配，并且海域水深由東南方向向陸地邊界逐漸變淺。該區(qū)域水深深度在0～200m左右。

海域風場數(shù)據(jù)是由NASA(美國國家航天局)提供的CCMP(Cross-Calibrated,Multi-Platform)衛(wèi)星遙感海面風場數(shù)據(jù)。風速區(qū)域范圍26°N～32°N,120°E～126°E，風場時間分辨率為6 h，空間分辨率為0.25° ×0.25°，輸入的時間長度為2011年1月1號0時-2011年9月3號18時，給出的是距海面10m處風場沿徑向向緯向的速度分量。

2.2 實驗點位置

JASON-2衛(wèi)星測量的有效波高具有較高的測量精度，與美國國家浮標數(shù)據(jù)中心（NDBC）實測數(shù)據(jù)對比分析的精度（均方根誤差）達到了0.28 m[23]，JASON-2衛(wèi)星是TOPEX/Poseidon和JASON-1的后繼雷達測高衛(wèi)星,主要用于測量有效波高和海平面，有效波高的測量精度達到2～3 cm[17]。JASON-2衛(wèi)星不僅可以結(jié)合TOPEX/Poseidon、JASON-1衛(wèi)星共同探測中國海平面時空變化，而且還可以利用JASON-2星載GPS數(shù)據(jù)采用簡化動力學(xué)方法進行JASON-2衛(wèi)星精密定軌[18-19]。圖1給出了東海水深與JASON-2軌跡分布，該衛(wèi)星經(jīng)過研究區(qū)域的pass共有5條，其中選取JASON-2測高衛(wèi)星在上升和下降地面軌跡的交點，如圖1所示交叉點A，B，C，作為研究SWAN模式模擬時間域上有效波高的實驗點。交叉點A是pass51和pass62的交點，交叉點B是pass51和pass240的交點，交叉點C是pass127和pass62的交點，交叉點坐標分別為A(123.315 9°E, 29.725 4°N)，B(121.898 2°E,26.887 2°N),C(124.733 3°E, 26.889 1°N)。

圖1 東海水深與JASON-2軌跡分布

2.3 模擬過程

SWAN模式在給定風場、水深情況下能夠獲得比較理想的波浪要素[13]，模擬的數(shù)值可與浮標實測數(shù)據(jù)或衛(wèi)星測高數(shù)據(jù)進行對比分析。其模擬過程如圖2所示。

圖2 流程圖

在該模擬過程中，模式參數(shù)的設(shè)定和水深網(wǎng)格點數(shù)據(jù)的加密是重要的步驟。其中就這2個步驟給予了描述。

模式計算中，最小截斷頻的率為0.01 Hz，最大截斷頻率為1.4 Hz,波數(shù)空間為15°的分辨率。在swan模式中，空間傳播的波使用上風格式即空間backward格式（BSBT格式），模式中共開啟2種類型的耗散機制，即波浪破碎（使用了大量實驗數(shù)據(jù)分析的破碎參數(shù)的平均值0.79）和底摩擦耗散（采用了經(jīng)驗性的JONSWAP模式中的0.038m2s-3底摩阻系數(shù)），其他參數(shù)項使用了SWAN模式手冊中的默認值，比如，三幅波-波相互作用項選取了由Madsen and Sorenen提供的相互作用系數(shù)J[22]。為了使該時刻的輸出值不受初始場的影響，提前一定時間開始計算，開始時間是2010年12月29號0時0分，結(jié)束時間為2011年9月3號18時，由于JASON-2衛(wèi)星采集數(shù)據(jù)的時間大約是238 h一個周期，人為的將SWAN模式的計算時間和數(shù)據(jù)輸出時間間隔設(shè)為2 h，數(shù)據(jù)輸出的開始時間與衛(wèi)星采集數(shù)據(jù)的時間相匹配，這就使得SWAN模式每個數(shù)據(jù)輸出的時刻更接近于JASON-2衛(wèi)星采集數(shù)據(jù)的時間，降低了由于時間的延遲引起的誤差。

在原有水深分辨率1'×1'的規(guī)則網(wǎng)格的基礎(chǔ)上進行局部三角形不規(guī)則網(wǎng)格加密，分別以交叉點的位置為中心，以邊長為6 km的正方形區(qū)域為加密的海域，進行逐步的加密，如圖1所示正方形包含的交叉點A,B和C區(qū)域就是所需加密的海域。分別選取3個加密的海域，創(chuàng)建3個區(qū)域邊界特征弧段，設(shè)置邊界特征弧段的分辨率，并使用線性內(nèi)插的方式進行插值，這就使得3個區(qū)域的水深數(shù)據(jù)在經(jīng)度和緯度方向上都得到了與邊界特征弧段相同的分辨率。這3個區(qū)域空間水深數(shù)據(jù)加密的不規(guī)則三角形網(wǎng)格的分辨率取值是1'-0.5'，取值間隔為0.1’。交叉點A所在的區(qū)域內(nèi)水深深度在40～60m左右，等值線相對來說比較稀疏，水深梯度較小。B區(qū)域內(nèi)水深深度是60～100 m左右，C區(qū)域內(nèi)水深深度在100～140m左右，C區(qū)域和B區(qū)域內(nèi)水深等值線分布相對較密，這2個區(qū)域內(nèi)的水深梯度比A區(qū)海域要大。

在整個計算區(qū)域中，一邊為海岸線，其它3邊為開邊界，3個驗證的海域距離計算區(qū)域邊界較遠，在該加密的區(qū)域中計算出的波浪數(shù)值不受開邊界的影響。

3 結(jié)果分析與討論

為了研究SWAN模式基于不同網(wǎng)格水深數(shù)據(jù)的分辨率對運算耗時的影響，本文主要選取了6個（1'-0.5'）不同的網(wǎng)格分辨率進行實驗，在加密的三角形水深網(wǎng)格分辨率從低到高的情況下，相同的計算機配置 (Pentium(R)Dual-Core CPU E6600 @3.06GHz 3.07GHz，4.00GB，Intel(R)G41 Express Chipset(MicrosoftCorporation-WDDM 1.1)，SWAN模式模擬結(jié)果所耗時時間如表1所示，表1中為了研究耗時時間的變化趨勢，定義了耗時時間的上升幅度。

其中上升幅度定義如下：

式中：t為當前時刻；ft為當前分辨率下計算機運算耗時時間的上升幅度；Δt為當前分辨率下的運算耗時時間；Δt-1為基于上一個三角形網(wǎng)格分辨率設(shè)備耗時時間。

如表1所示，三角形網(wǎng)格分辨率從1'-0.5'規(guī)則的遞增，計算機耗時時間就會大幅度延長，上升幅度呈現(xiàn)遞增的趨勢。出現(xiàn)這樣的現(xiàn)象是因為加密的三角形網(wǎng)格水深數(shù)據(jù)的節(jié)點增多，例如，0.5'分辨率的三角形網(wǎng)格水深數(shù)據(jù)節(jié)點數(shù)為61 835，1'分辨率的三角形網(wǎng)格水深數(shù)據(jù)節(jié)點數(shù)為35 675，前者大概是后者的1.7倍，數(shù)據(jù)節(jié)點的增加會使得單機運算降低。

為了分析計算機耗時時間的發(fā)展趨勢，主要選取表1中不同的三角形網(wǎng)格分辨率數(shù)據(jù)和耗時時間數(shù)值作為對象，繪制時間序列圖分析結(jié)果如圖3所示，同一設(shè)備下的運算耗時時間是隨著不規(guī)則三角形網(wǎng)格分辨率的增加呈現(xiàn)指數(shù)似形式的增長。

表1 三角形網(wǎng)格分辨率與運算耗時時間的關(guān)系

圖3 同一設(shè)備下耗時時間變化的趨勢圖

經(jīng)分析知，Jason-2衛(wèi)星測高數(shù)據(jù)在2011年1-8月份的匹配數(shù)據(jù)是25個?；诓煌切尉W(wǎng)格分辨率下模擬的有效波高在時間上的變化，通過編程提取與JASON-2衛(wèi)星數(shù)據(jù)在時間和空間上相一致的模擬值，將3個區(qū)域中不同三角形網(wǎng)格分辨率下的模擬值與衛(wèi)星實測值進行對比分析，對比結(jié)果如圖4所示，圖4中橫軸表示JASON-2衛(wèi)星的實測值，縱軸表示SWAN模擬的有效波高，單位都是m，其中圖4中的第一行從左至右的6個小圖，依次代表的是a區(qū)域內(nèi)的從低到高的三角形網(wǎng)格分辨率下的模擬值與實測值的散點圖，第二行代表的是b區(qū)域內(nèi)的，第三行指的是c區(qū)域內(nèi)的。

圖4 不同分辨率的三角形網(wǎng)格在3個區(qū)域中的模擬結(jié)果與Jason-2衛(wèi)星實測值的散點分布圖

由圖4可知，大多數(shù)點基本上都在對角線附近，這表明基于不同三角形網(wǎng)格分辨率的SWAN模式在3個區(qū)域中模擬的有效波高與衛(wèi)星測高值均具有較好的相關(guān)性。

由于模式的參數(shù)化和近似計算使得SWAN模式不可能完全準確地反演波浪的生成，模擬結(jié)果和衛(wèi)星測高在數(shù)值上仍存在著顯著的偏差。在同一個區(qū)域，隨著三角形網(wǎng)格分辨率的遞增，大多數(shù)點趨向于對角線的位置略有靠近但不明顯。在不同的區(qū)域中同樣分辨率的三角形網(wǎng)格模擬的波高數(shù)值存在著顯著的不同，相對于a區(qū)域的散點分布，水深梯度較大的b區(qū)域和c區(qū)域的更為密集且靠近于對角線。由圖4還可以看出，b區(qū)域和c區(qū)域的散點分布較相接近，在2 m以內(nèi)的有效波高，模擬值大部分大于衛(wèi)星實測值，也就是說散點基本上都在對角線以上，說明模擬值較實測值有高估趨勢。當有效波高大于2m，大部分的散點在對角線以下分布較多，表明在b區(qū)域和c區(qū)域處模擬值較實測值有低估趨勢。

為了定量地分析和評定不同三角形網(wǎng)格分辨率下的SWAN模式模擬的計算結(jié)果，采用了均方根誤差這個統(tǒng)計指標作為本文的誤差分析標準，主要采用2011年1-8月份SWAN模式模擬的有效波高數(shù)據(jù)，選取與JASON-2衛(wèi)星測高在時間域上相適應(yīng)的波高數(shù)值，取模擬值和真實值之間的殘差值獲取該模擬值的均方根誤差，并把此均方根誤差作為本文中所要分析的模擬精度。水深梯度間接的反應(yīng)了海底地形的坡度[20],文中主要選取不同的三角形水深網(wǎng)格分辨率代表海底地形刻畫的程度。

表2 不同分辨率中模擬精度的統(tǒng)計分析

由表2可見，SWAN模式模擬的波高精度存在明顯的區(qū)域差異。采取同樣的三角形水深網(wǎng)格分辨率，c區(qū)域模擬精度最高，其次是a區(qū)域，b區(qū)域最低。在c區(qū)域中，模擬精度提高幅度的平均值是0.56%，一般是a區(qū)域精度上升的1.5倍以上，同等于b區(qū)域內(nèi)的精度上升幅度。結(jié)合圖1和表2可知，SWAN模式模擬的精度存在著由水深梯度小的a區(qū)域向水深梯度大的c區(qū)域遞增的規(guī)律?；谙嗤切嗡罹W(wǎng)格分辨率，除去區(qū)域不同，模擬精度不同之外，在三角形水深網(wǎng)格分辨率不同的情況下，區(qū)域相同，SWAN模擬的精度也會存在明顯的差異現(xiàn)象。在文中3個加密的區(qū)域，隨著水深網(wǎng)格分辨率規(guī)則的遞增，SWAN模式模擬的波高精度并不是規(guī)則的沿著網(wǎng)格分辨率的提高而變化，a區(qū)域和b區(qū)域在精度上升幅度的趨勢上具有相同的變化規(guī)律，減小增加再減小，但c區(qū)域內(nèi)精度上升幅度存在著顯著的不同，呈現(xiàn)增加減小的現(xiàn)象。同一地區(qū)，雖然在精度上升幅度上出現(xiàn)顯著差異，但總體而言，在同一水深和水深梯度下，SWAN模式模擬的精度隨著三角形水深網(wǎng)格分辨率的增加而提高，即在有限水深情況下，能夠更細致地刻畫海底地形的加密的三角形網(wǎng)格，與水深梯度有較大的關(guān)系，而與水深深度沒有明顯的關(guān)系，更一步地說明除去風、水流、水深等因素對海浪傳播有很大影響外[21]，海底地形（水深梯度）刻畫的程度同樣可以影響模擬的精度。因為加密的三角形網(wǎng)格是對于海底水深梯度的描述，不再涉及到水深深度的大小，當三角形網(wǎng)格的分辨率增加時模擬的有效波高的精度增加，模擬結(jié)果的分析可知，在此文中，水深深度不是影響模擬結(jié)果的主要因素，海底的水深梯度是關(guān)鍵性的影響。

3個交叉點中精度最大提高幅度的位置基本上都是在基于三角形水深網(wǎng)格分辨率為 0.6'的SWAN模式下模擬的結(jié)果，隨之精度上升幅度減小，即在3個區(qū)域內(nèi)精度基本上沒有明顯的提高，出現(xiàn)這種現(xiàn)象可能原因是海域地形內(nèi)0.6'三角形網(wǎng)格的分辨率就足以刻畫該地形區(qū)域內(nèi)的水深梯度，由于相鄰網(wǎng)格中水位差而影響SWAN模式模擬的精度因素就會降低。整體來講，相對于未加密的三角形網(wǎng)格分辨率1'來說，基于0.6'三角形網(wǎng)格分辨率取得的精度上升的幅度中，C區(qū)域和B區(qū)域同等于最大，A區(qū)域最小，分別為2.8%和1.8%。出現(xiàn)這樣的結(jié)果主要與這3個區(qū)域內(nèi)的水深梯度有著密切的關(guān)系，由于A區(qū)域內(nèi)水深等值線較稀疏，水深梯度較小，采用局部加密的三角形網(wǎng)格獲得的有效波高的精度不會太明顯，而B區(qū)域和C區(qū)域水深等值線較密集，SWAN模式得到的結(jié)果逐漸接近真實值。因此來說，水深梯度較大的區(qū)域，采用局部加密三角形網(wǎng)格模擬的有效波高趨向于真實值才越明顯。

結(jié)合表2和圖4，再次說明，同一海域地形，采用局部加密的三角形網(wǎng)格能夠使SWAN模式模擬的有效波高趨近于真實值，空間域上，基于不同三角形網(wǎng)格分辨率的SWAN模式，水深梯度增加，模擬精度變化的偏差就越大且向著精度較好的方向發(fā)展。

4 結(jié)論

基于SWAN模式，采用不同水深三角形網(wǎng)格分辨率，對東海有效波高進行了模擬，利用JASON-2測高交叉點的實測數(shù)據(jù)進行對比分析，得出如下結(jié)論。

（1）采用局部加密的三角形網(wǎng)格能夠使得SWAN模式模擬的精度提高，隨著三角形網(wǎng)格分辨率的提高，模擬的有效波高效果就會越好，但同時，計算機運行耗時時間就會越長，近似呈現(xiàn)指數(shù)似增加，即運算效率就會降低。

（2）在本研究區(qū)域內(nèi)，通過對采用1'～0.5'加密的三角形網(wǎng)格所模擬的有效波高的統(tǒng)計量可知，采用的0.6'的三角形網(wǎng)格的分辨率獲得模擬結(jié)果中精度幅度提高是最大的，說明0.6'分辨率的三角形網(wǎng)格足以對該區(qū)域內(nèi)的水深梯度進行精細化。

（3）通過對交叉點A，B和C區(qū)域內(nèi)的水深和水深梯度，以及三角形網(wǎng)格的分辨率和所獲得的結(jié)果進行分析可知，局部加密的三角形網(wǎng)格對于水深梯度較大的海域具有很好的適用性。

基于水深三角形網(wǎng)格加密的SWAN模式對復(fù)雜海域地形下的波浪模擬能力較好，對于平坦地形，作用不明顯，當水深梯度較大時，模型結(jié)果的精度就會越高就越接近真實值，但當三角形網(wǎng)格的分辨率已經(jīng)能夠很準確地描述該地形的基本特征時，若再提高分辨率，模型的精度會處于穩(wěn)定的狀態(tài)，不再提高，反而會使模型運行的時間加長，運算效率降低。

[1]王輝,劉娜,李本霞,等.海洋可預(yù)報性和集合性預(yù)報研究綜述[J].地球科學(xué)進展,2014,29(11):1212-1225.

[2]齊義泉,朱伯承,施平,等.WWATCH模式模擬南海海浪場的結(jié)果分析[J].海洋學(xué)報,2003,25(4):1-9.

[3]Team S.SWAN Scientific and Technical Documentation,WAN Cycle IIIVersion 40.81[R].The Netherlands:Delft University of Technology,2010.

[4]Ris R C,Holthuijsen L H,Booij N.A Third-Generation Wave Model for Coastal Region 2.Verification[J].Journal of Geophysical Research,1999,104(C4):7667-7681.

[5]Gorman RM,Neilson CG.Modeling ShallowWaterWaveGeneration and Transformation in an IntertidalEstuary[J].CoastalEngineering, 1999,36:197-217.

[6]王道龍,華鋒,江志輝.SWAN近岸海浪模式在遼東灣的應(yīng)用[J].海洋科學(xué)進展,2010,28(3):285-291.

[7]梅嬋娟,趙棟梁,史劍.兩種海浪模式對中國黃海海域浪高模擬能力的比較[J].海洋預(yù)報,2008,25(2):92-98.

[8]王殿志,張慶河,時鐘,等.渤海灣風浪場的數(shù)值模擬[J].海洋通報,2004,23(5):10-17.

[9]鄭崇偉,鄭宇艷,陳洪春.基于SWAN模式的近10年南海北部波浪能源資源研究[J].亞熱帶資源與環(huán)境學(xué)報,2011,6(2):54-59

[10]宗芳伊,吳克儉.基于近20年的SWAN模式海浪模擬結(jié)果的南海波浪能分布、變化研究[J].海洋湖沼通報,2014(3):1-12

[11]馮芒,張文靜,李巖,等.臺灣海峽及近岸區(qū)域精細化海浪數(shù)值預(yù)報系統(tǒng)[J].海洋預(yù)報,2013,30(2):42-48.

[12]徐麗麗,肖文軍,石少華,等.WaveWatchIII和SWAN模型嵌套技術(shù)在業(yè)務(wù)化海浪預(yù)報系統(tǒng)中的應(yīng)用及檢驗 [J].海洋通報, 2015，34(3):283-294.

[13]楊德周,尹寶樹,徐艷青,等.SWAN淺水波浪模式在渤海的應(yīng)用研究-Phillips線性增長比例系數(shù)的改進 [J].水科學(xué)進展, 2005,16(5):710-714.

[14]宗芳伊,魏來,吳克儉.基于SWAN模式的2種網(wǎng)格對嵊泗海域海浪模擬的比較研究[J].海洋湖沼通報,2014,(4):10-18.

[15]史劍,王璞,鐘中,等.兩種網(wǎng)格下的swan模式對黃渤海海域模擬比較[J].海洋預(yù)報,2011,28(4):48-57.

[16]梁書秀,孫昭晨,尹洪強,等.基于SWAN模式的南海臺風浪推算的影響因素分析[J].海洋科學(xué)進展,2015,33(1):19-30.

[17]郭金運,常曉濤,黃金維,等.基于衛(wèi)星測高交叉點的海洋表面地轉(zhuǎn)流速度[J].地球物理學(xué)報,2010,53(11):2582-2589.

[18]秦建,郭金運,孔巧麗,等.Jason-2衛(wèi)星星載GPS數(shù)據(jù)cm級精密定軌[J].武漢大學(xué)學(xué)報:信息科學(xué)版,2014,39(2):137-141.

[19]郭金運,王建波,胡志博,等.由TOPEX/Poseidon和Jason-1/2探測的1993-2012中國海海平面時空變化[J].地球物理學(xué)報, 2015,58(9):3103-3120.

[20]歐陽明達,孫中苗,翟振和.基于重力地質(zhì)法的南中國海海底地形反演[J].地球物理學(xué)報,2014,57(9):2756-2765.

[21]馮芒,沙文鈺,朱首賢.近岸海浪幾種數(shù)值計算模型的比較[J].海洋預(yù)報,2003,20(1):52-59.

[22]The SWAN Team.SWAN TechnicalDocumentation[R].The Netherlands:DelftUniversity of Technology,2013.

[23]葉小敏,林明森,宋慶君.基于現(xiàn)場觀測數(shù)據(jù)的衛(wèi)星雷達高度計海面風速和有效波高真實性檢驗方法研究[J].遙感技術(shù)與應(yīng)用,2014,29(1):26-32.

Simulative Analysis on the SignificantWave Height over the East China Sea by SWAN Modelwith Jason-2 Satellite Altimetric Crossover Points

XIE You-ge1,GUO Jin-yun1,ZHU Jin-shan1,LIU Xin1,LIGuo-wei2
1.College of Geodesy and Geomatics,Shandong University of Science and Technology,Qingdao 266590,Shandong Province,China; 2.WeihaiMarine EnvironmentalMonitoring Center,Weihai264209,Shandong Province,China

To research the effects of different triangular grid resolutions employed for the wave field variations over the East China Sea,this paper selects three crossover points of the Jason-2 satellite along the ascending and descending tracks and its surrounding local seas as the testing areas from January to August,2011.Through comparing significantwave height determined by the Jason-2 satellite,it is showed that the simulating precision gradually and disorderly increases,with the computing efficiency exponentially loweringwith improving triangular grid resolution from 1'to 0.5'at the interval of 0.1'.At the same time,statistical data show that the largest rising amplitude of precision appears for the triangular grid resolution of 0.6'.In addition,the refined triangularmesh produces its effects mainly in the southeastern area of the East China Sea where the submarine topography gradient is very large.

SWAN model;Jason-2 satellite;significantwave height;triangular grid

P731；P228

A

1003-2029（2017）01-0024-07

10.3969/j.issn.1003-2029.2017.01.005

2016-01-20

國家自然科學(xué)基金資助項目(41374009);山東省自然科學(xué)基金資助項目(ZR2013DM009);國家科技基礎(chǔ)性工作專項資助項目(2015FY310200)

謝友鴿(1990-),女,主要從事海洋大地測量研究。E-mail:1126396365@qq.com

郭金運，教授、博導(dǎo)，主要從事海洋大地測量、空間大地測量和物理大地測量等研究。E-mail:jinyunguo1@126.com