英曉明，丁平興

（華東師范大學 河口海岸學國家重點實驗室，上海 200062）

洋山港海域水體和懸沙輸運機制研究

英曉明，丁平興

（華東師范大學 河口海岸學國家重點實驗室，上海 200062）

對洋山港海域水體和懸沙長期輸運速度分離，并進行機制分解，分析各種物理過程對水體和懸沙長期輸運速度所起的作用。結(jié)果表明，水體和懸沙輸運存在明顯差異，直接描述懸沙凈輸運更準確合理；歐拉余流在水體長期輸運中占主導地位，平流輸運和潮泵輸運是洋山港海域懸沙輸運的主要動力因素，由于潮動力混合較強，垂向輸運作用很小。

水沙長期輸運速度；機制；洋山港海域

Abstract：The long-term transport velocity of water and suspended sediment are separated, mechanism decomposition is made on them, and effects of all kinds of physical process by observed data is analyzed in Yangshan Harbor waters in this paper.The results show that, the obvious difference exists between water and suspended sediment residual transport, and it is more accurate and reasonable to calculate residual sediment transport directly.Euler transport plays predominant role in water long-term transport, and horizontal advection and tidal pumping take primary role in suspended sediment transport in the study area, while the contribution of vertical transport is small due to strong mixed tide.

Keywords：Long-term transport velocity of water and suspended sediment; mechanism; Yangshan Harbor Waters

1 引 言

在河口海岸地區(qū)，海流作用下的懸沙輸運引起地形沖淤變化，進而影響岸灘演變、航道疏浚及水下岸坡的穩(wěn)定性，懸沙的長期輸運成為河口海岸研究的核心問題之一[1,2]。潮流及其余流對河口海岸地區(qū)懸沙的輸運起著非常重要的作用,由于潮流漲落潮不對稱輸沙，因而余流在某種程度上可以指示懸沙的凈輸運。然而，余流并不能完全代表懸沙凈輸運。實際上，河口存在潮周期平均物質(zhì)輸運與潮周期平均水體輸運方向不同的現(xiàn)象[3]。在懸沙空間分布較均勻的地方，懸沙和水體的長期（一個或若干個潮周期）輸運特征基本一致；在懸沙空間分布差異大的區(qū)域，懸沙和水體長期輸運存在明顯差異[4]。通量分析是研究物質(zhì)長期輸運較好的方法[5,6]。懸沙長期輸運速度比水體長期輸運速度描述懸沙輸運更準確合理[4]。朱首賢等[3]從理論上分析了河口地區(qū)存在物質(zhì)和水體長期輸運分離的自然條件。通過ADCP走航觀測研究，洋山港地區(qū)斷面存在周期內(nèi)凈潮通量與凈懸沙通量相反的現(xiàn)象[7]，說明存在水體和懸沙輸運分離現(xiàn)象。

李玉中等[8]通過1996年實測資料，證實了洋山港海域下層水體存在余流分離和匯聚現(xiàn)象。然而，關于洋山港海域懸沙長期輸運的研究尚未見報道。洋山港海域懸沙長期輸運有何特征？水體長期輸運（余流）能否代表懸沙長期輸運？本文將運用通量分析方法揭示洋山港海域水體和懸沙長期輸運的特征，并對機理進行研究。

2 研究區(qū)域概況

洋山港海域位于杭州灣崎嶇列島內(nèi)，具有多島嶼和多汊道輸水輸沙的特征。峽道由以大、小洋山為主的南北島鏈組成，呈現(xiàn)為向西擴展的喇叭形狀（圖1）。港區(qū)最大水深約87 m（黃海高程）；開闊處淺灘的水深為10 ～ 14m，具有建造深水港的優(yōu)良條件。自2002年以來，洋山深水港區(qū)相繼實施了堵汊工程、陸域填充工程和港池疏浚等一系列重大工程。港區(qū)一期碼頭和二期碼頭分別于2005年和2006年年底開始運營。目前，小洋山北島鏈僅剩顆珠山-蔣公柱汊道。洋山港海域?qū)儆诜钦?guī)半日潮，以M2分潮為主，該水域在東海潮波傳播、長江口沖淡水擴散和島鏈峽道地形的綜合制約下, 具有波浪弱、潮流強和懸沙濃度高的特點。該海區(qū)平均潮差2.79 m，年平均波高為0.4 m，冬夏半年平均波高相差不大[9]。落潮平均流速為0.90 ～ 1.20m/s，漲潮均流速為0.80 ～1.20 m/s，最大可能潮流流速在3.00 m/s以上，潮流為往復流[10]。洋山港海域泥沙主要為長江直接擴散泥沙和潮流攜帶的海域泥沙，年平均懸沙濃度約為1.3 kg/m3，冬高夏低。

圖1 洋山港工程及測量站位Fig.1 Sketch of Yangshan Harbor projects and survey stations

3 資料與研究方法

本文資料主要來源于交通部天津水運工程科學研究所2006年5月的定點水文泥沙觀測資料。測量位置見圖1。觀測資料包括大、中和小潮各25 h的數(shù)據(jù)，垂向為六點法，對應水深為表層、0.2 H、0.4H、0.6 H、0.8 H和底層，時間間隔為一小時。為綜合考慮懸沙和水體的長期輸運效應，對大、中和小潮的資料進行全潮分析，計算懸沙和水體長期輸運速度及其相關項。

為使水體長期輸運和懸沙長期輸運能夠進行比較，本文采用朱首賢[3,6]定義的水體和懸沙長期輸運速度方法來分析洋山港海域二、三維水體和懸沙長期輸運特征及進行機理分析。

瞬時流速、含沙量和水深由式(1)—(4)[11,12]表達:

4 結(jié)果和討論

4.1 水體和懸沙輸運速度特征及比較

截止2006年5月，港區(qū)3個堵汊工程均已經(jīng)竣工，測量站位分布較廣，具有一定的代表性，因此選擇此次資料進行分析。二維懸沙和水體長期輸運速度如圖2(a) 所示，垂向給出表層、中層和底層懸沙和水體長期輸運速度，分別如圖2(b-d) 所示。

港區(qū)北側(cè) N1測站二維懸沙和水體長期輸運方向基本一致，前者是后者量值的2.1倍；中層和底層懸沙和水體長期輸運速度方向基本一致，但表層懸沙和水體長期輸運速度方向相差39°，且表層量值上前者是后者的3.5倍，表層和底層懸沙輸運方向相差114°，表層和底層水體輸運方向相差134°。N2測站二維懸沙和水體長期輸運方向相差69°，底層懸沙和水體長期輸運速度方向基本一致，主要是表層和中層兩者速度方向差異顯著，表層和底層懸沙長期輸運速度方向相差 119°，表層和底層水體輸運速度方向相差95°。N3測站二維懸沙和水體輸運速度方向相差146°，量值上后者是前者的1.7倍，垂向上懸沙和水體長期輸運速度方向均差異顯著，但垂向上懸沙輸運速度方向相差不大，量值上表層是底層的5.7倍。

港區(qū)西側(cè)W1測站二維懸沙和水體長期輸運速度大小和方向相差不大，垂向上亦如此，表層懸沙輸運速度方向和底層相差63°。W2測站二維懸沙和水體長期輸運方向基本一致，前者是后者量值的2.2倍；垂向上懸沙和水體長期輸運方向相差不大，表層和底層懸沙輸運速度方向相差96°，量值上前者是后者的5.4倍，表層和底層水體輸運速度方向相差138°。W3測站二維懸沙和水體長期輸運速度相差不大，垂向上底層懸沙和水體長期輸運速度方向相差101°，表層懸沙輸運速度方向和底層相差69°，量值上前者是后者的 3倍，表層和底層水體輸運速度方向相差135°。

港區(qū)南側(cè)S1測站二維懸沙和水體長期輸運速度差異不大，垂向上亦如此，垂向上表層和底層懸沙輸運方向相差不大，量值上前者是后者的1.8倍。

港區(qū)東側(cè)P4測站二維懸沙和水體長期輸運方向相差很小，量值上前者是后者的2.5倍，垂向上表層和中層懸沙和水體長期輸運方向相差很小，但底層懸沙和水體長期輸運方向相差 159°，量值上后者是前者的1.5倍，表層懸沙輸運速度量值是底層的7.3倍，表層水體和中下層水體輸運方向相反。E1測站二維懸沙和水體長期輸運速度相差不大，垂向上亦如此，表層懸沙輸運速度量值是底層的 1.9倍。E2測站二維懸沙和水體長期輸運速度相差不大，表層和底層懸沙輸運方向相差73°，表層和底層水體輸運方向相差39°。

港區(qū)顆珠山汊道 K1測站二維懸沙和水體輸運方向相差55°，量值上后者是前者的1.6倍，垂向上底層懸沙和水體輸運速度基本一致，但表層和中層懸沙和水體輸運速度方向差異顯著，方向差異與兩者垂向平均差異相差不大，底層懸沙輸運速度量值是表層的1.9倍。K2測站二維懸沙和水體輸運方向相差不大，垂向上亦如此，表層和底層懸沙輸運速度方向相差64°，量值上相差不大。

港區(qū)雙連山汊道SC1測站二維懸沙和水體輸運方向基本一致，垂向上亦如此，底層懸沙輸運速度量值是水體的2.6倍。SC2測站二維懸沙和水體輸運速度相差很小，垂向上亦如此，表層懸沙輸運速度量值是底層的3.6倍。港區(qū)大山塘汊道SC3測站二維懸沙和水體輸運方向相差54°，量值相近，垂向上兩者方向差異與垂向平均差異值相差不大。

港區(qū)內(nèi)部 K4測站二維懸沙和水體輸運方向相差59°，量值上前者是后者的1.8倍，垂向上兩者差異變化不大。P1測站二維懸沙和水體輸運速度相差不大，垂向上兩者差異變化不大，僅表層兩者方向相差稍大為 40°。P2測站二維懸沙和水體輸運速度量值相近，方向相差25°，垂向上兩者差異由表至底減小，表層兩者方向相差為39°，底層兩者方向基本一致。

以上分析說明，懸沙和水體長期輸運速度大小和方向存在明顯差異，在描述洋山港海域懸沙長期輸運時，懸沙長期輸運速度比水體長期輸運速度更準確合理。

在洋山港北側(cè)N1和N2測站、西側(cè)W2和W3測站和東側(cè)P4測站，表層和底層水體輸運方向相反或差異較大。而懸沙輸運不僅在N1、N2、W2、W3測站，而且在W1、K1、K2、P2和E2測站，表層和底層輸運方向相反或差異較大，但P4測站表層和底層懸沙輸運方向確相差很小。在東側(cè)E1和E2測站水體輸運方向相反。

在中、南部P1、P2和W1測站形成懸沙輸運方向輻聚區(qū)，此區(qū)域位于2006年4月到2007年4月形成的港區(qū)淤積帶（圖3），港區(qū)南側(cè)落潮流和落潮沙占優(yōu)，而中軸漲潮流和漲潮沙占優(yōu)，兩股流向切變區(qū)流速必然減少，加之懸沙匯聚是造成港區(qū)淤積的原因之一。

圖2 懸沙和水體長期輸運速度Fig.2 Water and suspended sediment long-term transport velocity

圖3 2006年-2007年地形沖淤分布（據(jù)三航院資料繪制）Fig.3 Deposition/erosion distribution from 2006 to 2007

表1 二維水體長期輸運速度及其分解項（量值：cm/s 方向：°）Tab.1 2-D long-term transport velocity of water and its components (magnitude: cm/s direction: °)

表2 二維懸沙長期輸運速度及其分解項（量值：cm/s 方向：°）Tab.2 2-D long-term suspended sediment transport velocity and its components (magnitude: cm/s direction: °)

4.2 二維懸沙和水體長期輸運速度機制分析

二維水體長期輸運速度分量列在表1中。由表1可以分析各物理過程對水體輸運速度的貢獻。在大部分測站（W1、K4、K1、K2、P2、E1、E2、P4、N1、N3、SC1、SC2和 S1），歐拉余流在水體長期輸運中占主導地位；在 W3測站歐拉余流量值是斯托克斯漂移速度1.3倍，兩者方向相反；在W2和SC3測站斯托克斯漂移速度大于歐拉余流，兩者方向相反；在P1測站歐拉余流和斯托克斯漂移速度量值相近，方向相差92°。洋山港海域相當多測站余流超過0.1 m/s，該區(qū)域?qū)儆谟嗔鬏^強區(qū)域。由于島嶼附近地形集聚效應，因此W1、E1和E2等測站余流值較大。

二維懸沙長期輸運速度分量列在表2中。從表2中可以看出，T1、T2、T5占主要地位，T4和T6僅在部分測站作用相對重要，T3的作用很小，T7和T8幾乎不起作用。在大部分測站（與水體輸運一樣），T1量值最大，說明洋山海域歐拉輸運是主導因素。在W3測站歐拉輸運和斯托克斯漂移輸運方向相反，量值前者是后者的1.3倍。在W2測站歐拉輸運和斯托克斯漂移輸運方向相反，量值后者是前者的1.4倍。在P1測站歐拉輸運和斯托克斯漂移輸運方向相差92°，兩者量值相差不大。在SC3測站歐拉輸運和斯托克斯漂移輸運方向相反，兩者量值相差不大。

T5在有些測站（W1、K4、K2、SC2）與T2量值相近，在很多測站（K1、P1、P2、E1、E2、P4、N3、N2、N1、SC1和S1）T5比T2量值都要大的多，說明洋山海域潮泵輸運作用地位僅次于歐拉輸運。T4僅在W3、W1、K4、P1和SC3測站相對重要，T6在K4測站相對重要。

為直觀了解懸沙在不同區(qū)域的輸運機制，將平流輸運項（T1+T2）、潮泵輸運項（T3+T4+T5）和垂向輸運項（T6+T7+T8）與懸沙長期輸運速度（Uc）繪制在一起，如圖4所示。平流輸運項是洋山港海域泥沙輸運的主導因素，潮泵輸運所起作用居于第二位，垂向剪切擴散和切變作用不大。由于洋山港海域是強潮混合區(qū)，垂向混合作用大，因此垂向作用不大。

5 結(jié) 語

本文對洋山港海域懸沙和水體長期輸運分離計算，并對懸沙和物質(zhì)輸運速度進行機制分解。研究結(jié)果表明，在洋山港海域懸沙輸運速度比水體輸運速度描述懸沙輸運更準確合理。

圖4 懸沙長期輸運速度及三個合成分量Fig.4 Long-term transport velocity of suspended sediment and its components

機制分解結(jié)果表明，歐拉余流在水體長期輸運中占主導地位，平流輸運和潮泵作用是洋山海域懸沙輸運的主要動力因素，由于潮動力混合較強，垂向輸運作用很小。懸沙和水體的分層輸運分析表明，洋山海域懸沙和水體輸運垂向存在差異，例如，在N2測站表層水沙和底層水沙輸運方向均相反，在P4測站表層水體和底層水體輸運方向相反，而懸沙輸運方向相同。港區(qū)中南部懸沙凈輸運指向大山塘北側(cè)，可能是該地區(qū)淤積的原因之一。由于測站偏少，利用實測資料得到的懸沙長期輸運尚不能明確說明港區(qū)內(nèi)部地形沖淤演變機理，通過補充觀測或者利用泥沙數(shù)值模擬研究港區(qū)懸沙長期輸運變化，將會更好的分析地形沖淤演變機理。

[1]Rose P C, Thorne P D.Measurements of suspended sediment transport parameters in a tidal estuary [J].Cont Shelf Resh, 2001, (21): 1 551-1 575.

[2]陳建勇, 戴志軍, 陳吉余, 等.杭州灣北岸弧形岸段懸沙凈輸移分析—以龍泉—南竹港岸段為例 [J].泥沙研究, 2009, 2: 53-59.

[3]朱首賢, 丁平興, 沙文鈺, 等.河口物質(zhì)和水體長期輸運分離的理論分析和觀測驗證Ⅰ物質(zhì)和水體長期輸運分離的理論分析 [J].海洋學報, 2008, 30(6): 24-29.

[4]張文靜, 朱首賢, 沙文鈺.二維懸沙長期輸運速度的定義和機理分析 [J].海洋工程, 2003, 21(3): 56-61.

[5]Dyer K R.Coastal and Estuarine Sediment Dynamics [M].A Wiley Interscience Publication John Wiley and Sons, 1985:253-254.

[6]朱首賢.流、浪模式和物質(zhì)長期輸運分離研究 [D].華東師范大學博士學位論文, 2005.

[7]虞志英, 張志林.上海市國際航運中心洋山深水港區(qū)西港區(qū)海床沖淤趨勢和增深研究 [R].華東師范大學河口海岸學國家重點實驗室,長江水利委員會長江口水文水資源勘測局, 2009, 3.

[8]李玉中，陳沈良.洋山港海域余流分離和會聚現(xiàn)象研究 [J].水利學報, 2003, 5: 24-34.

[9]俞航, 陳沈良, 谷國傳.崎嶇列島海區(qū)水沙特征及近期沖淤演變[J].海岸工程, 2008, 27(1): 10-20.

[10]劉蔡胤.上海洋山港港區(qū)海域潮流泥沙數(shù)值模擬及水下地形變化特征分析 [D].大連理工大學碩士學位論文, 2008.

[11]Uncles R J, Elliott R C A, Weston S A.Observed fluxes of water, salt and suspended sediment in a partly mixed estuary [J].Estuar Coast Shelf Sci, 1985, 20(2): 147-167.

[12]沈健, 沈煥庭, 潘定安, 等.長江河口最大渾濁帶水沙輸運機制分析 [J].地理學報, 1995, 50(5): 411-420.

Research on transport mechanism of water and suspended sediment in the Yangshan Harbor waters

YING Xiao-ming, DING Ping-xing

(State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai 200062, China)

P731.2

A

1001-6932(2011)02-0135-07

2010-02-23；收修改稿日期：2010-10-18

國家海洋局公益性項目 ( 200905001-9 )。

英曉明 ( 1979－ )，男，博士，主要從事河口海岸動力學研究。電子郵箱：xmying2005@yahoo.com.cn。