陸罕芳， 顧峰峰， 陳學恩??， 戚定滿

(1.中國海洋大學海洋與大氣學院，山東 青島 266100； 2.上海河口海岸科學研究中心，上海 201201)

舟山群島位于長江口東南側(cè)，北面連接長江口海域，西面與杭州灣相連，東面連接東中國海。群島島嶼眾多，岸線地勢復雜多變。舟山群島海域(見圖1)既是杭州灣水體與東海水體交換的主要通道，也是杭州灣向東海輸運泥沙的重要通道，還是東海沿岸流攜帶泥沙向南輸運的過渡地帶[1]。在舟山群島海域，開闊海域沉積類型以粉砂為主，島嶼附近則是以細砂為主，黏土主要來源于從長江口直接進入杭州灣的泥沙，外海輸沙的影響很小[1-2]。表層沉積物總體上呈東粗(砂)西細(泥)、粒徑由北向南逐漸變細的特點[3-4]。舟山群島海域懸沙含量東南方向低、西北方向高，自杭州灣向外海迅速降低[1，5]。落潮平均懸沙含量一般高于漲潮平均懸沙含量，大潮間懸沙平均含量大于小潮間，冬季懸沙平均含量大于夏季。潮流是控制泥沙輸運的主要動力因子[1，6-7]。懸沙輸運具有往復搬運、凈輸沙量小的特點[1]。舟山群島的存在對海流變化、泥沙分布和泥沙輸運有顯著影響，具有明顯的“群島效應”[5]。由杭州灣向外海，泥沙起動流速越來越大，在群島內(nèi)，泥沙基本沿水道走向輸運[1]。

一般天氣條件下，舟山群島海域波浪主要為風浪，呈現(xiàn)東高西低特征：在群島以東，平均波高1 m左右，在群島以西，由于群島阻礙作用，波浪向西傳播時迅速衰減[8]。舟山群島海域夏秋臺風(熱帶風暴)影響頻繁。臺風天氣條件下，舟山群島對臺風浪的阻擋效果亦顯著。臺風期間舟山東部海域波高大、周期長且涌浪影響明顯，西部海域波高較小且以風浪影響為主[9]。已有研究表明，波浪能通過影響流速、底部切應力和紊動擴散系數(shù)來影響泥沙的起動、懸沙的輸運和擴散[10-13]。故而泥沙沖淤研究中，特別是基于數(shù)值模型模擬研究時考慮波浪的影響是十分有必要的。尹超等[14]采用ECOMSED模式進行了不同波浪動力條件影響下山東半島東側(cè)黑泥灣海域沖淤演化的數(shù)值模擬，結(jié)果表明，隨機波浪要素的引入對沖淤演變模擬結(jié)果有較大的影響，引入波浪要素后的模擬結(jié)果與海圖更為接近。舟山群島海域泥沙輸運數(shù)值模擬研究方面，二維泥沙模型主要有Delft 2D[15]、SWEM2D[16]以及潮流模型如ADCIRC耦合懸沙輸運基本方程[17-18]等；三維泥沙模型主要有ECOMSED[19]、Delft 3D[20]FVCOM[21-23]等。但上述數(shù)值模擬研究或只針對舟山群島海域中某一特定小區(qū)域，或缺少對舟山群島整體海域泥沙沖淤特征詳細的分析，且均未分析波浪作用對泥沙輸運的影響。

(b2：嵊山驗潮站，b1：佘山驗潮站； a2：浮標Q206，a1：浮標Q207； 1～9：懸沙濃度測站；c*：流速測站；d：春曉油管水深測點)b2: Shengshan tide station, b1:Sheshan tide station; a2:buoy Q206,a1:buoy Q207; 1～9:sediment suspended concentration(SSC); c*: velocity; d:Chunxiao pipeline water depth).)

綜合前人的研究可知，現(xiàn)有的舟山群島海域泥沙輸運的研究多為利用一次或多次有限站點的調(diào)查采樣數(shù)據(jù)來分析舟山群島海域整體的懸沙和沉積物的分布特征、泥沙輸運方向和通量以及動力機制，或者是針對某個具體水道、具體工程范圍的泥沙輸運數(shù)值模擬研究。但整個舟山群島海域一般天氣條件下的泥沙輸運特征尚不清晰，波浪對舟山群島海域泥沙沖淤影響以及不同的波浪條件下產(chǎn)生的影響差異有待進一步深入研究。

因而，本文在長江口與杭州灣及其鄰近海域建立了一個高分辨的三維水沙鹽模型SWEM，并以舟山群島海域為重點研究海域；研究分別模擬了一般天氣條件下，該海域在潮汐、徑流、風和波浪影響下的地形沖淤與懸沙含量分布特征，并對一般天氣條件下和臺風天氣條件下波浪對該海域泥沙沖淤過程的影響及其差異進行了和分析比較。

1 研究方法與數(shù)據(jù)簡介

SWEM是上海河口海岸科學研究中心自主研發(fā)的三維潮流泥沙數(shù)值模型[24]。該模型基于無結(jié)構(gòu)三角網(wǎng)格和有限體積離散，具有很好的復雜邊界適應能力和守恒性；對于流場計算采用了歐拉—拉格朗日追蹤法，計算效率較高，可以滿足較大范圍、復雜邊界區(qū)域的研究需要。本文模式計算網(wǎng)格范圍(見圖2)120.31°E～124.13°E，28.46°N～33°N，包括長江口、杭州灣、舟山群島及其鄰近海域，其中長江口區(qū)域采用了1∶2000的實測地形資料，其余采用了該區(qū)域2010年的海圖地形資料。

圖2 SWEM模式水平網(wǎng)格Fig.2 SWEM model horizontal mesh

模式采用冷啟動，對2010年全年進行了逐月模擬。開邊界采用16個主要分潮驅(qū)動；模式考慮長江、錢塘江徑流量；長江上游采用大通站2010年的日平均流量，錢塘江采用多年平均的年平均流量；采用WOA09氣候態(tài)逐月特征溫度均一場；泥沙邊界條件取開邊界零通量，泥沙初始場取為長江口口門外小范圍(121.4°E～122.3°E，31.2°N～31.8°N)海域0.3 g/mL，其他海域為零初始場；泥沙沉降速度與臨界起動應力隨著海域泥沙粒徑由西向東逐漸增大，沉降速度：0.25～0.5 mm/s，臨界起動應力：0.4～0.8 Pa；沖刷率為0.000 2 kg/m3,沉降概率為0.2；考慮波流共同作用下的底部切應力。

CCMP(Cross Calibrated Multi-Platform)風場數(shù)據(jù)集與浮標QF206(見圖1:a1)風場數(shù)據(jù)(見圖3:a,b)對比結(jié)果顯示，該數(shù)據(jù)集在本研究海域可信度較高。故而一般天氣條件下的風場數(shù)據(jù)采用2010年CCMP(Cross Calibrated Multi-Platform)數(shù)據(jù)集(0.25°×0.25°，6 h/次)插值到計算網(wǎng)格上。選取2014年經(jīng)過舟山群島海域的“鳳凰”臺風(見圖4)，并基于臺風特征因素，通過常用的Jelesnianski-I[25]風場模型重建得到臺風風場。相應天氣下的波浪場由SWAN海浪模式給出。SWAN模型波浪驗證結(jié)果見圖3c。模式實驗設(shè)置如下(見表1)：一般天氣條件下考慮波浪(實驗1-1)，一般天氣條件下不考慮波浪(實驗1-2)，臺風天氣條件下考慮波浪(實驗2-1)，臺風天氣條件下不考慮波浪(實驗2-2)。

圖3 CCMP風場(a、b)和SWAN波浪驗證(c) Fig.3 CCMP data compare with the buoy data (a，b)and SWAN significant wave high verification(c)

此外，由于缺少相應時段的實測資料，根據(jù)可用實測資料的時段，本文選取2006年10月、2010年1月、2010年5月、2012年7月分別模擬，進行泥沙、水位和流速驗證。模式驅(qū)動條件均為該時段的潮汐、CCMP風場(2012年采用NECP后報風場)、對應風場下的波浪場、日平均長江徑流與錢塘江多年的年平均徑流。

圖4 臺風路徑示意圖Fig.4 Typhoon path

表1 實驗設(shè)置Table1 Model experiments setup

2 模式結(jié)果驗證

模式結(jié)果驗證分為水動力驗證和泥沙(懸沙含量、地形沖淤)驗證兩部分。所用實測數(shù)據(jù)測點見圖1。潮位實測數(shù)據(jù)采用佘山、嵊山測站(見圖1:b1,b2)2010年5月的逐時水位數(shù)據(jù)(見圖5)。流速數(shù)據(jù)采用2010年1月7日—9日浮標QF207(見圖1:a2)的垂向分層逐時流速數(shù)據(jù)(見圖6)與2012年7月20日14時至21日16時測點SW1(見圖1:c)的垂向平均逐時流速數(shù)據(jù)(見圖7)。

泥沙含量驗證(見圖8)采用2006年10月青島海洋地質(zhì)研究所和中國海洋大學在該海域進行的9個站位(見圖1:1～9)、25 h周日平均數(shù)據(jù)[1]。地形沖淤變化采用春曉油管海域2007年和2013年兩次相隔6年的實測地形數(shù)據(jù)(見圖1:d)之差作為實測數(shù)據(jù)。因此，本文將2010年全年地形沖淤變化的模擬結(jié)果乘以6倍，與觀測作對比(見圖9)。

(S1: 佘山站；S2: 嵊山站。 S1: Sheshan station；S2: Shengshan station.)

圖6 模式垂向分層流速結(jié)果驗證Fig.6 Model vertical velocity verification

圖7 模式平均流速結(jié)果驗證Fig.7 Model mean velocityverification

圖8 模式懸沙含量結(jié)果驗證Fig.8 Model SSC verification

圖9 模式地形沖淤變化結(jié)果驗證Fig.9 Model micro-topography change verification

從圖5～9的模式驗證結(jié)果來看，模式結(jié)果和實測數(shù)據(jù)吻合較好，流速與實測數(shù)據(jù)基本接近，模式懸沙含量模擬結(jié)果與實測結(jié)果變化趨勢一致，均呈現(xiàn)杭州灣1～3測點懸沙含量高，舟山群島海域4～9測點懸沙含量較低的趨勢，模式模擬的地形沖淤變化與實測地形沖淤變化趨勢一致。綜上表明，模式結(jié)果能較好的反映研究海域的水動力、懸沙分布以及地形沖淤特征，可以用來進行后續(xù)研究分析，其中計算的誤差主要和該區(qū)域較為復雜的地形及風浪條件有關(guān)。

3 結(jié)果分析討論

3.1 一般天氣條件下地形沖淤和懸沙含量時空分布特征

為研究一般天氣條件下地形沖淤變化和懸沙含量的分布特征，選取2、5、7、8、10與12月作為典型月份。取模式逐月運行結(jié)果后15 d的懸沙含量進行時間平均，作為當月潮周期平均的懸沙含量分布，與當月地形沖淤變化模擬結(jié)果繪制圖10。

舟山群島海域島嶼間東西方向水道基本表現(xiàn)為沖刷地形，如螺頭水道、舟山—岱山水道等。沖刷量表現(xiàn)為2月份開始增強，在7月達到最大值10 cm，秋季開始減弱，至12月或1月最弱。島嶼間南北方向的水域如岱山島—大長途間水域和朱家尖東側(cè)等海域則表現(xiàn)為較強的淤積特征。淤積量基本表現(xiàn)為秋季開始增強，在10或12月達到最大，最大淤積量分別為7和14 cm，春季開始減弱，至5月達到最小值。海域懸沙含量由西北向東南迅速降低,懸沙含量在群島西北側(cè)海域一般為1 g/L左右，在群島東南側(cè)海域則降低為0.1 g/L，體現(xiàn)了群島“慮沙”的特點。冬季懸沙平均含量大于夏季。整體上看3月份懸沙含量最大，此時1 g/L 等值線深入群島內(nèi)部，8月份懸沙含量最小，此時1 g/L等值線遠離舟山群島海域。懸沙含量分布特征與胡日軍[1]、黃慧明等[5，26]、壽瑋瑋[27]和周鴻權(quán)等[7]的研究結(jié)果吻合。懸沙含量的分布與潮汐余流場分布吻合，表明了潮汐是控制懸沙輸運的主要因素。來自長江口的懸沙隨著潮波運動，進入舟山群島海域，結(jié)合局地掀沙，在潮波與地形作用下或沿著幾個主要的水道向外海擴散，或在流速較小的區(qū)域沉積。

3.2 一般天氣條件下波浪對泥沙沖淤過程的影響

波浪模擬結(jié)果表明波浪自舟山群島東側(cè)傳入，受到群島阻礙作用迅速減弱。波浪分布自西向東呈帶狀分布，等值線密集。海域秋冬季節(jié)的波浪強于春夏季節(jié)，其中5月最弱，10和12月最強(具體的SWAN波浪模擬結(jié)果本文未給出)。

((填充圖為地形沖淤變化/m；等值線為懸沙含量。藍色0.1:0.2:0.5 g/L；綠色1:0.5:1.5 g/L；紅色2 g/L。The fill color, /m and SSC distribution (the contour). Blue 0.1:0.2:0.5 g/L； Green 1:0.5:1.5 g/L； Red 2 g/L) .)

由考慮波浪作用(實驗1-1)的模式結(jié)果與不考慮波浪作用(實驗1-2)的模式結(jié)果，得到波浪對舟山群島海域泥沙沖淤過程的影響的分布變化(見圖11)。從時間分布上看，舟山群島海域秋冬季受波浪影響較大，春夏季受影響較小。波浪對泥沙沖淤過程影響較顯著的區(qū)域為朱家尖東側(cè)海域。波浪對該海域泥沙沖淤過程的影響表現(xiàn)為促進泥沙淤積和增加懸沙含量。在該海域，波浪影響顯著的月份為1—3月和10—12月，其中10月的影響最大，波浪影響增加的淤積量為10.7 mm，約占總體淤積量的10%。在其他月份的影響很小，特別在4、5月波浪的影響幾乎可以忽略。結(jié)合Soulsby底應力的計算公式[28]與黃慧明等[26]給出的2007年11月朱家尖東北側(cè)海域水質(zhì)點的運動軌跡圖，推測該海域受波浪影響形成促進淤積、增加懸沙含量特征的原因可能是：相對于不考慮波浪作用，考慮波浪作用情況下整體海域的底應力增大，使得海域水體懸沙含量有所增加。隨著水體輸運到朱家尖東側(cè)海域的懸沙也隨之增多，從而朱家尖東側(cè)海域懸沙含量增加。同時在該海域，水質(zhì)點停留時間長，故而泥沙沉降增多。

3.3 臺風天氣條件下波浪對泥沙沖淤過程的影響

為研究在波浪條件變化劇烈的情況下，波浪對舟山群島海域泥沙沖淤過程的影響是否會發(fā)生較大變化，進行了臺風天氣條件下舟山群島海域泥沙沖淤過程的模擬。臺風自2014年9月22日15時進入舟山群島海域，到2014年9月23日21時離開舟山群島海域，一共歷時30 h。在臺風過程中，有效波高高值區(qū)集中在舟山群島外側(cè)(見圖12)，并隨著臺風中心向北移動。最大有效波高達到6 m。

由實驗2-1與實驗2-2得出2014年“鳳凰”臺風風浪對舟山群島海域泥沙沖淤過程的影響，同時選取2010年9月份波浪的影響圖，繪制圖13，對比討論臺風天氣條件下和一般天氣條件下波浪對舟山群島海域泥沙沖淤過程的影響。

((填充圖為地形沖淤變化/m；等值線為懸沙含量；藍色0.01:0.02:0.05 g/L；綠色0.1:0.05:0.15 g/L；紅色0.2 g/L).The fill color, /m and SSC distribution (the contour). Blue 0.01:0.02:0.05 g/L； Green 0.1:0.05:0.15 g/L； Red 0.2 g/L.).)

圖12 臺風天氣條件下有效波高(填充圖,/m)和風場分布特征(等值線, m/s)Fig.12 Significant wave high (the fill color, /m) and wind speed distribution (the contour，m/s) under typhoon weather condition

圖13 臺風天氣條件下波浪引起的地形沖淤變化與泥沙含量分布變化(a)和一般天氣條件下9月波浪引起的地形沖淤變化與懸沙含量分布變化(b) (填充圖為地形沖淤變化/m；等值線為懸沙含量g/L)Fig.13 Micro-topography change(the fill color,/m)and SSC distribution(the contour,g/L) induced by wave influence under typhoon weather condition(a) and in September under general weather condition(b)

由圖13可知，在臺風天氣條件下，舟山群島大部分海域受波浪影響，地形沖刷增強。波浪對泥沙沖淤過程的影響在朱家尖東側(cè)海域仍然較為顯著。在該海域，波浪對影響表現(xiàn)為促進泥沙的淤積和增加懸沙含量。在臺風天氣條件下，波浪引起的朱家尖東側(cè)海域的最大淤積量為4.2 mm。在一般天氣條件下，波浪9月在朱家尖東側(cè)海域的影響不顯著，最大淤積量僅為0.8 mm。臺風天氣條件下波浪在朱家尖東側(cè)海域引起的淤積量為一般天氣條件下的6倍，且引起的淤積范圍也有所增大。對比一般天氣情況，臺風天氣的懸沙含量有所增加?？傮w而言，在臺風天氣條件下，波浪對除朱家尖東側(cè)海域外的其他舟山海域的泥沙沖淤過程的影響為促進沖刷和增加含沙量，在朱家尖東側(cè)海域則依舊表現(xiàn)為促進淤積和增加含沙量。

4 結(jié)論

本文基于SWEM三維水沙鹽模式，考慮潮汐、徑流、風場和波浪影響，在長江口與杭州灣及其鄰近海域建立高分辨率的三維水沙鹽模型?；谀Ｊ浇Y(jié)果，本文研究了舟山群島海域在一般天氣條件下地形沖淤與懸沙含量的分布特征。分析了一般天氣條件下和臺風天氣條件下波浪對舟山群島海域泥沙沖淤過程的影響。得出主要結(jié)論如下：

(1)一般天氣條件下，舟山群島海域內(nèi)部東西方向水道為沖刷特征，舟山群島東側(cè)海域以及舟山群島內(nèi)部南北方向水道為淤積特征。以沖刷為主的海域，沖刷強度春夏逐漸增強，至6月左右沖刷強度最大，然后開始減弱，至10月達到最小值。以淤積特征為主的海域，淤積強度秋季開始增強，10或12月淤積強度最大，春季開始減弱，在5月達到最小值。

(2)一般天氣條件下，舟山群島海域懸沙含量分布由西北向東南迅速降低。懸沙含量3月份最大，8月份最小，冬季懸沙平均含量大于夏季。懸沙含量的分布表明，潮汐是控制懸沙輸運的主要因素。

(3)一般天氣條件下，波浪對舟山群島海域泥沙沖淤過程影響較為顯著的海域為朱家尖東側(cè)海域。在該海域波浪的影響為促進淤積和增加懸沙含量。波浪影響的強度秋冬季增強，夏季最小。波浪影響的最大淤積量約占總淤積量的10%。臺風天氣條件下，波浪對除朱家尖東側(cè)海域外的其他舟山群島海域的泥沙沖淤過程的影響為促進沖刷和增加含沙量，在朱家尖東側(cè)海域則依舊表現(xiàn)為促進淤積和增加含沙量。

本文模擬了一般天氣下舟山群島海域的地形沖淤與懸沙含量分布特征，并探究了一般天氣與臺風天氣下波浪對該海域泥沙沖淤過程的影響。但海域受波浪的影響與具體的波浪因素如波高、波向、周期等之間的關(guān)系還有待進一步研究。

致謝：國家超級計算濟南中心提供了“神威藍光”計算機系統(tǒng)的支持。