潮汐和潮流影響下蘇北輻射沙洲海域波浪模擬分析

湯志華，鄭曉琴，曹翔宇，顧云碧，肖文軍，張蓓

(1.河海大學海岸災害及防護教育部重點實驗室，江蘇南京210098; 2.國家海洋局東海預報中心，上海200081; 3.河海大學企業管理學院，江蘇常州213022 )

?

潮汐和潮流影響下蘇北輻射沙洲海域波浪模擬分析

湯志華1，鄭曉琴2，曹翔宇3，顧云碧2，肖文軍2，張蓓2

(1.河海大學海岸災害及防護教育部重點實驗室，江蘇南京210098; 2.國家海洋局東海預報中心，上海200081; 3.河海大學企業管理學院，江蘇常州213022 )

摘要：基于第三代海浪模型SWAN，采用自嵌套的方法提供譜邊界條件，對影響蘇北輻射沙洲海域的一次冷空氣過程和一次臺風過程作用下的波浪進行了模擬。考慮到沙洲海域強潮水動力環境，分析了潮位和潮流的變化對該海域波浪的影響。結果表明，沙洲處波高和波周期受潮位影響顯著，受潮流影響弱，具有潮周期起伏的特點，而波向受潮位潮流影響不顯著；考慮高潮位后，以弶港為界，南北輻射沙洲波高顯著增加的區域與波浪傳播方向有關：波浪由北向南傳播，相差不大，波浪由東向西傳播，北部明顯大于南部。

關鍵詞：輻射沙洲；SWAN模型；波浪模擬；潮位；潮流

1　引言

輻射沙洲位于射陽河口至長江口北岸近岸淺水區，南北延伸200 km，東西橫跨90 km，由70多條沙脊和分布其間的潮流通道組成的。各條沙脊高低不等，形態各異，沙脊之間有深槽相隔，深槽坡陡水深各沙脊大小不等，寬窄不一，地形極其復雜。輻射沙洲為江蘇省發展提供了潛在的后備土地資源，適合發展大農業、旅游業、港口運輸業、風電場產業和圍墾造陸等[1]。波浪作為輻射沙洲區域的主要海岸動力之一，在塑造獨特的輻射沙脊地形方面起著至關重要的作用[2]。在江蘇沿海大開發的背景下，研究該海域的波浪特征對海洋資源的開發利用以及沿岸的經濟建設具有重要的意義。

在海浪模擬中，對臺風浪的模擬尤為重視，臺風浪的破壞性眾所周知，不少學者進行了研究[3-5]，但冷空氣風浪往往被輕視。研究表明冷空氣帶來的大浪往往具有嚴重的破壞性，且臺風浪的頻率遠不如冷空氣出現的頻率高[6]。冷空氣浪數值模擬相對于臺風浪數值模擬較少。周兆黎等[7]曾利用WW3模式對影響南海的三次冷空氣過程進行過模擬，發現在冷空氣活動末期以及近岸地區，模擬效果略偏差一些，但誤差仍在可接受的范圍內。

輻射沙洲海域的波浪模擬也有不少研究。邱桔斐[8]通過第三代淺水波浪數值模型SWAN模擬江蘇沿海風浪場，并對其特征進行分析研究；陳波等[9]采用SWAN模型對江蘇如東附近海域風浪場進行了數學模擬，并分析了SWAN模型在該海域風浪場模擬方面的適用性；楊耀中等[10]利用波浪折射繞射模型模擬了蘇北輻射沙洲區的波浪傳播變形過程；李杰等[11]利用風浪譜模型對蘇北輻射沙洲海域波浪的傳播過程進行數值模擬。這些研究對于了解蘇北輻射沙洲區的波浪特征是很重要的，但依然有需要進一步研究的地方，主要為：蘇北輻射沙洲區屬于強潮海域，潮差大，潮流強，潮波屬正規半日潮型，受太平洋前進潮波系統和南黃海旋轉潮波系統的控制，形成以弶港為中心的輻射狀潮流場，弶港附近岸邊出現特大潮差[2]。輻射沙洲外部海域主要表現為旋轉流特征，越往沙洲中心區域，往復性增強。各潮汐通道內漲落潮流沿著水道方向作定向往復運動，越靠近輻射沙洲頂部，往復性越明顯，在強天氣過程影響下，考慮潮汐潮流的影響下的波浪場模擬的研究較少。實際上，在有較強潮流存在的海岸河口水域,流場對波浪場的影響是不容忽視的[12]。肖文軍等[13]考慮潮汐和流影響進行了長江口波浪場數值計算，王彪等[14]進行了長興島海區波流相互作用數值模擬研究，這些均表明考慮潮影響的波浪模擬效果更好；另一方面，蘇北輻射沙洲區地形復雜，已有的研究缺乏精細的水深地形資料，難以準確模擬刻畫該海域的波浪特征。

本文采用高精度實測地形數據，基于物理過程考慮較為全面的第三代海浪數值模式SWAN，以WRF風場為驅動場，采用自嵌套方法，對襲擊蘇北輻射沙洲海域的一次冷空氣過程和臺風過程引起的海浪進行了數值模擬研究?？紤]到沙洲海域強潮水動力環境，分析了潮位和潮流對該海域波浪的影響。

2　模型設置

海浪模型采用近年來廣泛使用的第三代近岸海洋數值模型SWAN。在有水流影響時，波能密度不守恒但波作用量守恒，SWAN模型采用動譜平衡方程[15]作為描述海浪的控制方程，并且采用無結構網格[16]，能很好地擬合復雜的岸線地形，在潮流、地形、風場等復雜環境影響下，仍能準確地模擬出波浪場。在波浪向近岸傳播的過程中，近岸波浪場復雜多變，影響因素復雜，如折射繞射、淺水變形、波浪破碎和波流相互作用等現象，SWAN模型在緩坡方程模型的基礎上，采用相位平均的方法加入繞射作用，使得模式在淺水海域的計算精度得以提高。有關SWAN模型進一步詳細的介紹，可參考文獻[15]。

本文SWAN模型采用大小區域嵌套：大區域如圖1a，計算范圍為115°—134°E，16°—41°N，包括29193個計算節點，15423個三角形；小區域網格如圖1b，包括了整個江蘇沿海海域、部分山東沿海海域以及長江口、杭州灣附近海域，外?？臻g分辨率最大，在蘇北輻射沙洲海域進行局部加密處理，其網格分辨率約為300—700 m，最高空間分辨率為100 m左右，三角形個數為114166，網格點數為57892個。模型采用球坐標系下的非定常模式，大區域為小區域提供波浪譜邊界條件，頻率從0.041—1 Hz，以對數分布劃分為60個；方向的分段為72個，分辨率為5°。模式初始條件選擇由計算時刻的風場和JONSWAP經驗公式[17]給定的模式，時間步長為10 min?？紤]底摩阻、三相波非線性相互作用、四相波非線性相互作用、波浪破碎和波浪繞射等物理過程的影響，參數的設置為默認值。風場模型采用中尺度大氣模式WRF，由東海預報中心提供，時間分辨率為1 h。

天文潮模型選用FVCOM模型。FVCOM是無結構網格的、有限體積的、三維原始方程的海洋模式。模型方程包括動量方程、連續方程、溫鹽守恒方程以及狀態方程，通過采用耦合了Mellor 和Yamada的2.5階湍封閉模型來對方程進行封閉。水平方向上三角網格，而在垂向上采用的是σ坐標，為了節省計算機機時，其采用內外模態交替計算的方法。關于模型的詳細介紹和設置詳見文獻[18]，本文的潮流模型結果是基于文獻[18]的研究成果，由于文獻[18]已經對其做過詳細的驗證，本文不再重新驗證。

研究海域地形復雜，對該海域波浪的精確模擬離不開高精度的水深地形資料。本文模型所用水深來自“908”專項以及國家海洋公益性項目《蘇北淺灘“怪潮”災害監測預警關鍵技術研究及示范應用》進行的高精度測量1:5000和1:10000的水深數據，長江口臨近海域采用2009年實測，外海海域采用ETOPO全球1′×1′水深數據。模型局部水深如圖1c所示，從圖中可看出，水深地形資料對該海域的地形刻畫得比較精細。

3　海浪模擬分析

3.1冷空氣過程模擬分析

該海域淺水浮標（洋口港，地理位置如圖1c）從2012年夏季開始進行觀測，本文選擇兩年來較強的一次冷空氣過程進行模擬。冷空氣模擬時間為2012年12月27日8時—2012年12月31日0時，農歷時間為11月15日—11月19日，處于天文大潮期間。

為了分別考慮潮位和潮流對波浪的不同影響程度，將波浪模擬分為4種：不考慮潮位潮流影響、單獨考慮潮位影響、單獨考慮潮流影響以及綜合考慮潮位潮流影響，將4種模擬與實測進行對比并分析。

圖1　計算網格及研究區域

圖2　洋口港浮標冷空氣過程不同要素圖

圖2為洋口港的對比圖，洋口港浮標處水深為10 m，周圍地形復雜，除了東邊水深較深外，其余周邊都是沙脊，水深較淺，波浪傳播影響因素復雜。為了減少初始狀態的影響，只分析12月28 日20:00以后的結果。有效波高的模擬值與實測值變化趨勢基本保持一致。此次波浪過程模擬基本反映了冷空氣影響前期、中期的海浪變化，但在冷空氣活動末期，模擬效果較差，這與周兆黎等[17]的研究結論是一致的，并且波浪衰減較慢。

圖2a顯示該點潮位對波高的演化過程影響很大，而潮流的影響較小。在潮位影響下，波高具有明顯的潮周期性起伏，高潮位下波浪增大，低潮位下波浪減小。由于該點周圍水深較淺，波浪的發展嚴重地受到水深的限制，水深的增加減少了底摩擦耗散，可以產生更大的波浪。在實測波高峰值處，該時刻潮位較高，考慮潮位后，模擬值與實測值吻合得更好，與不考慮潮位比較，波高可增加達0.3 m。在波高峰值處，潮位和潮流對波高的影響程度不同，雖然高潮位（見圖2d）使波高增大，但潮流的作用使波高減少，即潮位和潮流對波高的影響是相反的，互相抵消的；但由于潮流較小（見圖2d），其影響更弱，相對來說潮位的影響更大，潮位和潮流對波高的效應不能完全抵消，潮位和潮流總的影響使波高增大。冷空氣消亡階段，潮位和潮流對波高的影響更為明顯，考慮潮位和潮流的模擬波高與實測波高的變化趨勢非常一致，低潮位時，波高減少最多達0.3 m。圖2b顯示考慮潮位潮流后周期的變化規律與波高的變化規律是一致的，也是潮位占主導作用，也具有潮周期性起伏，考慮潮位潮流后周期的模擬更為合理，模擬趨勢與實測趨勢更為一致,最大差異可達0.9 s。從圖2c可以看出，考慮潮位潮流后波向（此處波向為笛卡爾坐標系中的去向，東為零度，逆時針為正，北為90°）變化不大，最大差異約為10°，其中潮流的作用比潮位大，但沒有波高和波周期的類似潮周期變化規律。

圖3　冷空氣過程波要素空間分布圖（單位：m）

從圖2a中可看出，12月30日2時左右波浪波高最大，該時刻潮位也較高，潮位潮流影響最大，圖3a顯示該時刻波浪由北向南傳播，此時不同影響下波高空間差異分布如圖3b、圖3c和圖3d所示。從圖中可以看出，潮流和潮位對波浪的影響在不同區域存在差異性。各潮汐通道存在弱的負作用區域，波高稍微減少，一般減少0.1 m，最大減少達0.2 m，而由于潮位增高波高略有增加，潮汐通道雖然水深較深，但波浪傳播還受到周圍沙洲的影響，潮位高時，沙洲處波高增加，潮汐通道波高也會增加。除了北部一些潮汐通道特別是西洋通道外，潮位及潮流對波高的總的效應是增加的，最大增加不超過0.3 m。淺灘沙脊處，潮流的影響很小，波高差異不超過0.1 m，潮位的影響是顯著的，高潮位時波高可增大0.4—0.8 m。而在沙洲外部海域，潮位影響不大，波高增大不超過0.1 m，相對而言潮流的影響稍大，一般波高增加不超過0.3 m。

3.2臺風過程模擬分析

臺風案例選取1210號臺風“達維”，1210號臺風在江蘇連云港登陸，臺風路徑如圖4，移動路徑為西北向，最強級別為臺風級別。模擬時間為2012 年7月31日8時—2012年8月4日8時，農歷時間為6月13日—6月17日，處于天文潮中大潮期。

此次臺風過程中，8月2日中午離洋口港浮標最近。從圖5a中可以看出，臺風來臨前，波高慢慢增大，臺風靠近，波高減小，臺風遠離后，波高又增加，較好地模擬了此次臺風過程中波高的變化趨勢，但由于此次臺風中心風速低估了，嚴重地低估了最靠近臺風期間的波浪。由于達維為登陸型臺風，其登陸后強度慢慢減小，對研究海域的影響時間長，從圖5a和圖5b中可以更加明顯地看出波高和周期受潮的影響，具有典型的潮周期性，相對于潮流來說，仍然是潮位的影響更大，這與之前冷空氣得出的結論是一致的。高潮位時，波高最大增加0. 3 m，周期最大增加1.3 s。圖5c顯示潮位和潮流對波向的影響仍然是不大，最大差異約為20°，還是潮流占主導作用。

圖4　1210號臺風“達維”路徑圖

圖5　洋口港浮標臺風過程不同要素圖

如圖5a所示，8月3日0時潮位潮流影響大，該時刻潮位高、潮流?。ㄒ妶D5d），臺風臨近登陸，圖6a顯示波浪由東向西傳播。為了節省篇幅這里僅討論潮位影響下的波浪場。從圖6b可以看出，以弶港為界，北部沙洲海域海浪受潮位影響顯著的區域明顯大于南部沙洲海域，而在冷空氣情況下南北差異并不明顯（見圖3c），這是輻射沙洲特殊的地形導致的。沙洲具有明顯的南北不對稱性，北部沙洲規模大，分布密集且較為連續，沙脊向東和東南方向延伸，延伸距離長，而南部沙洲規模小，分布較分散，沙脊向東和東南方向延伸，延伸距離較短[1]。此次臺風期間，波浪從外海傳播（見圖6a），由于地形原因導致北部波浪耗散多，而南部波浪耗散少，故考慮高潮位，北部沙洲海域波高顯著增加的區域比南部沙洲海域大。但在冷空氣期間，風向為偏北風，波浪由北向南傳播（見圖3a），南部沙洲海域波浪受到北部沙洲海域地形阻擋，波浪發展受到限制，考慮高潮位后，兩者差別不明顯。

圖6　臺風過程波要素空間分布圖

4　結論

本文采用高精度實測地形數據，基于第三代近岸海浪模型SWAN，以WRF風場為驅動場，采用自嵌套的方法提供譜邊界條件，對蘇北輻射沙洲海域在一次冷空氣過程和臺風過程作用下的波浪進行了模擬?？紤]到蘇北輻射沙洲海域潮動力強，潮流模型選用已經驗證過的適用該海域的FVCOM模型，分析了潮位和潮流對該海域波浪的影響。本文所得出主要結論如下：

（1）在輻射沙洲這樣的強潮海域，為了得到更為準確合理的波浪模擬結果，要考慮潮位和潮流的影響；

（2）沙洲處波高和波周期受潮位影響顯著，受潮流影響弱，具有潮周期起伏的特點，而波向受潮位潮流影響不顯著，其中潮流的作用比潮位大；

（3）對波高而言，潮汐通道處，潮流和潮位均有影響，西洋通道潮流影響大，而在沙洲外海，潮位影響很弱，相對來說潮流影響更大；

（4）考慮高潮位后，以弶港為界，南北輻射沙洲波高顯著增加的區域與波浪傳播方向有關：當波浪由北向南傳播時，如在冷空氣影響下，南北差異不大；當波浪由東向西傳播時，如在臺風影響下，北部明顯大于南部。

致謝：感謝東海預報中心提供的計算資源和觀測數據。

參考文獻：

[1]李孟國.輻射沙洲研究開發的進展[J].水道港口, 2011, 32(4): 229-243.

[2]張東生,張君倫,張長寬,等.潮流塑造-風暴破壞-潮流恢復-試釋黃海海底輻射沙脊群形成演變的動力機制[J].中國科學(D 輯), 1998, 28(5): 394-402.

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

[4] Ou S-H, Liau J-M, Hsu T-W, et al. Simulating typhoon waves by SWAN wave model in coastal waters of Taiwan [J]. Ocean Engineering, 2002, 29(8): 947-971.

[5]徐福敏,黃云峰,宋志堯.東中國海至長江口海域臺風浪特性的數值模擬研究[J].水動力學研究與進展, 2010, 23(6): 604-607.

[6] Zheng C W, Zhuang H, Li X, et al. Wind energy and wave energy resources assessment in the East China Sea and South China Sea [J]. Science China Technological Sciences, 2012, 55(1): 163-173.

[7]周兆黎,楊顯宇. WWATCHⅢ模式模擬南海冷空氣過程海浪場的有效性檢驗[J].熱帶海洋學報, 2008, 27(2): 1-6.

[8]楊耀中,馮衛兵.南黃海輻射沙脊群波浪場數值模擬[J].河海大學學報(自然科學版), 2010, 38(4): 457-461.

[9]邱桔斐.江蘇沿海風、浪特征研究[D].南京:河海大學, 2005.

[10]陳波,蔣昌波,陳漢寶.江蘇如東附近海域風浪場的數值模擬[J].水運工程, 2007(1): 20-24.

[11]李杰,袁建忠.風浪譜模型在蘇北輻射沙洲海域的應用研究[J].水利經濟, 2012, 30(4): 53-56.

[12]徐福敏,張長寬,茅麗華,等.一種淺水波浪數值模型的應用研究[J].水動力學研究與進展, 2000, 15(4): 429-434.

[13]肖文軍,丁平興,胡克林.潮汐和流影響下長江口波浪場數值計算[J].海洋工程, 2008, 26(4): 45-52.

[14]王彪,沈永明,王亮.長興島海區波流相互作用數值模擬研究[J].海洋工程, 2012, 30(3): 87-96.

[15] Booij N, Ris R C, Hojthuijsen L H. A third-generation wave model for coastal regions: 1. Model description and validation[J]. Journal of Geophysical Research, 1999, 104(C4): 7649-7666.

[16] Zijlema M. Computation of wind-wave spectra in coastal waters with SWAN on unstructured grids[J]. Coastal Engineering, 2010, 57(3): 267-277.

[17] Hasselmann K, Barnett T P, Bouws E, et al. Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP) [J]. Deutches Hydrographisches Zeitschrift, 1973, 8(12): 1-95.

[18]張蓓.蘇北輻射沙洲海域潮波及風暴增水特征研究[D].南京:河海大學, 2012.

Numerical simulation and analysis of waves with tide and currents in the radial sandbanks

TANG Zhi-hua1, ZHENG Xiao-qin2, CAO Xiang-yu3, GU Yun-bi2, XIAO Wen-jun2, ZHANG Bei2
(1. Key Laboratory of Coastal Disaster and Defence Ministry of Education, Hohai University, Nanjing 210098 China; 2. Forecast Centers of East China Sea, SOA, Shanghai 200081 China; 3. School of Business Administration, Hohai University, Changzhou 213022 China)

Abstract：Based on the third generation numerical wave model SWAN, using self-nested method to provide the spectral-type boundary, the process of wave caused by a cold air and typhoon in the radial sandbanks is simulated. Considering the strong hydrodynamic environment, the characteristics of wave influenced by changes of tide and current is analyzed in this area. The results show that, in the sandbanks, wave height and period are affected significantly by tide, and weakly by current, which have a tidal cycle. But wave direction is affected weakly both by tide and current. Considering the high tide level, demarcated by Jiang Port, the area of the south and north radial sandbanks, where wave height increases significantly, is related to wave direction. The difference of wave height is little when wave spreads from north to south. While wave spreads from east to west, the wave height in north area is significantly higher than that in south area.

Key words：radial sandbanks ; SWAN model; wave simulation; tide; current

通訊作者：鄭曉琴（1982-），女，工程師，碩士研究生，主要從事海洋環境預報研究。E-mail：christian1001@163.com

作者簡介：湯志華（1990-），男，碩士研究生，主要從事河口海岸動力學研究。E-mail：tangzhhu@163.com

基金項目：國家海洋局海洋公益性行業科研專項（201405022-3）；上海市科委項目（14DZ1205204）

收稿日期：2014-09-11

DOI:10.11737/j.issn.1003-0239.2015.02.004

中圖分類號：P731.22

文獻標識碼：A

文章編號：1003-0239（2015）02-0024-07