許靈靜, 楊　波??, 江文勝

(中國海洋大學 1. 海洋與大氣學院； 2.海洋環(huán)境與生態(tài)教育部重點實驗室，山東 青島 266100)

?

臺風對杭州灣水交換影響的定量研究?

許靈靜1, 楊波1??, 江文勝2

(中國海洋大學 1. 海洋與大氣學院； 2.海洋環(huán)境與生態(tài)教育部重點實驗室，山東 青島 266100)

摘要：基于ADCIRC(Advanced Circulation Model for Shelves, Coasts, and Estuaries)模式建立了杭州灣臺風風暴潮流二維數(shù)值模型，結(jié)合粒子追蹤模型，分別研究了天文潮和臺風過程作用下杭州灣的水交換情況，定量分析了臺風對杭州灣水交換的影響;統(tǒng)計7場典型路徑的臺風，結(jié)果表明：臺風影響比較劇烈時，與天文潮的單獨作用相比，灣內(nèi)、外參與交換的水體范圍分別可增加約1.6倍和7倍，灣內(nèi)流出的水體影響范圍可以增加31%；灣內(nèi)水交換速度可提高1.7倍。研究還發(fā)現(xiàn)不同路徑的臺風對杭州灣內(nèi)外水體交換的影響有所不同。

關(guān)鍵詞：杭州灣；水交換；臺風；數(shù)值模擬

引用格式：許靈靜, 楊波, 江文勝. 臺風對杭州灣水交換影響的定量研究[J]. 中國海洋大學學報(自然科學版)， 2016， 46(6): 29-41.

XU Ling-Jing, YANG Bo, JIANG Wen-Sheng. Quantitative study on water exchange affected by typhoons in the Hangzhou Bay[J]. Periodical of Ocean University of China, 2016, 46(6): 29-41.

杭州灣位于浙江省北部，北臨上海，西起澉浦，東接中國東海，并與長江口毗鄰。隨著城市經(jīng)濟的發(fā)展，大量的污染物質(zhì)如生活和工業(yè)廢水、疏浚物等被排入灣中，灣內(nèi)的污染物日益增多，富營養(yǎng)化日趨嚴重[1-2]，給灣內(nèi)環(huán)境造成很大壓力，杭州灣已成為中國近岸污染最嚴重的區(qū)域之一[3]。灣內(nèi)污染物的輸移擴散和物理自凈能力主要通過水交換的強弱來衡量[4]，因此研究杭州灣水交換對于保護杭州灣，進行灣內(nèi)環(huán)境規(guī)劃具有重要的指導意義。

目前，通過箱式模型、水質(zhì)模型、質(zhì)點追蹤模型等，學者們已對國內(nèi)外海灣的水交換進行了大量研究[5-10]。對杭州灣的水交換，國內(nèi)也已有一些研究：如張淑珍等[11]選取杭州灣內(nèi)4個擬定的排污口進行研究，發(fā)現(xiàn)杭州灣水交換最活躍的地方為灣口中部區(qū)域；劉新成等[12]計算了杭州灣和長江口的歐拉余流，定量地研究了長江口與杭州灣的水交換范圍，指出長江水向杭州灣的擴散主要集中在近岸南匯嘴至8m等深線以內(nèi)的范圍；Kong等[13]利用二維數(shù)值模型，分析了影響長江口和杭州灣水交換和沉積輸運的各種機制如潮位的變化、潮流動能、余流、粒子追蹤以及潮灘等，研究表明長江口與杭州灣之間存在一個周期性的水交換，凈的水交換主要發(fā)生在0～6m等深線范圍內(nèi)，凈水通量的最大值出現(xiàn)在2m等深線附近，水體的輸運路徑和范圍與沉積物的不同；姚炎明等[14]通過淡水組分法計算了包括杭州灣在內(nèi)的整個錢塘江河口的水體沖洗時間，認為沖洗時間與徑流量之間存在指數(shù)函數(shù)的關(guān)系；Wang等[15]研究發(fā)現(xiàn),長江與杭州灣的水交換受大小潮的調(diào)節(jié)，并隨季節(jié)發(fā)生變化，冬天西北風或北風的作用促進了長江水向杭州灣的輸運，水齡比夏天明顯減小。從前人的研究中可以看出，其研究多側(cè)重于長江口與杭州灣的水交換，對杭州灣內(nèi)的水體情況研究較少。

水交換的眾多研究表明：潮流是影響水交換的基本動力因素，風是影響水交換的重要動力要素，特別是臺風過程，對水交換的影響非常顯著。陳春華等[16]曾經(jīng)對8410和9111臺風在瓊州海峽和?？跒骋鸬乃粨Q進行研究，結(jié)果表明在一場臺風過程中風暴流輸運的水量可以使灣內(nèi)的污染物濃度降低15%左右。而前人對杭州灣水交換動力機制的探討主要集中于天文潮、季風和徑流，對于極端天氣情況下水交換的研究很少。

杭州灣是一個受臺風影響比較頻繁的地區(qū)，對中國氣象局熱帶氣旋資料中心(tcdata.typhoon.gov.cn)[17]提供的最佳路徑數(shù)據(jù)集進行統(tǒng)計，發(fā)現(xiàn)自建國以來平均每年約有1.7場臺風影響到杭州灣地區(qū)，因此研究臺風對杭州灣水交換的影響具有重要的科學意義。

本文通過數(shù)值模擬的方法探討不同臺風情境下杭州灣的水交換情況，定量地給出了臺風影響的程度，為改善杭州灣海水環(huán)境提供參考。

1模型介紹

1.1 ADCIRC模型

ADCIRC (Advanced Circulation Model for Shelves, Coasts, and Estuaries)由美國北卡羅來納州大學的Luettich[18]等人主持開發(fā)，基于非結(jié)構(gòu)三角網(wǎng)格，能很好地擬合復雜的岸線，采用空間有限元和時間有限差分相結(jié)合的方法既提高了區(qū)域的精度又提升了計算的效率。目前該模式已被廣泛應用于潮汐和風暴潮的預報中[19-22]。

雖然自然界中的水體運動都是三維問題，但杭州灣水深較淺，水平尺度遠遠大于垂向尺度，同時流速在垂向的大小和變化都遠小于水平方向，因此為簡化模型、減少計算量，本文采用不考慮斜壓作用的二維深度積分模型來模擬杭州灣的流場。笛卡爾坐標系下基本控制方程為：

(1)

(2)

(3)

1.2臺風模型

要模擬臺風引起的風暴流場，需將臺風過程中的氣壓場和風場作為模式的強迫場。臺風的位置、中心氣壓等基本信息可以從中國氣象局熱帶氣旋最佳路徑數(shù)據(jù)集[17]中得到，最大風速WR采用ATKINSON-HOLLIDY[24]提出的經(jīng)驗計算公式，最大風速半徑R采用MEF經(jīng)驗公式[25]，分別為：

WR=3.029×(P∞-P0)0.644，

(4)

(5)

其中：P∞為無窮遠處的大氣壓(取1013.25hPa)；P0為臺風的中心氣壓(單位hPa)；φ為臺風中心的緯度(單位為度)；V為臺風中心的移動速度；M為起算半徑，取45km。

臺風風場和氣壓場分布采用Jelesnianski模型[26]：

當0≤r≤R時

(6)

當r>R時

(7)

1.3 水交換模型

本文采用拉格朗日粒子追蹤模型進行水體運動的模擬，以Takeoka[27]提出的剩余函數(shù)概念表示水交換量。

2模型設(shè)置及驗證

2.1 ADCIRC模型設(shè)置

考慮到杭州灣的潮波能量主要由東中國海傳入[28]，將中國東海和黃海的部分區(qū)域(119.2°E～128.7°E，25.7°N～36.8°N)作為模型計算區(qū)域(見圖1)。模擬區(qū)域采用變分辨率網(wǎng)格，在開邊界網(wǎng)格步長設(shè)為10km，岸線處500m，島嶼周邊進行加密處理，杭州灣口處的大小洋山分辨率精細到150m。模型節(jié)點數(shù)132 089個，三角形單元259 262個，網(wǎng)格較好地擬合了計算區(qū)域內(nèi)的復雜岸線(見圖1)。模式中岸線資料結(jié)合衛(wèi)星遙感圖片和海圖數(shù)字化得到；水深數(shù)據(jù)來源于Choi(Laboratory for Coastal and Ocean Dynamics Studies Sungkyunkwan Univ.)提供的1’×1’地形數(shù)據(jù)與海軍航保部的3張海圖插值疊加。其中3張海圖分辨率為1∶250 000、1∶80 000、1∶35 000，分別對應中國東海舟山群島及附近、南匯嘴至火山列島、上海洋山港港區(qū)。

模型由天文潮、風場、氣壓場和徑流驅(qū)動。4個分潮(M2，S2，K1，O1)的振幅和遲角從OTPS(Oregon State University Tidal Prediction Software)中提取得到[29]；底應力采用混合非線性底摩擦公式，最小底摩擦系數(shù)取0.001；風應力是風速的函數(shù)。大通水文站測得長江洪峰期間的多年平均流量為56 200m3·s-1，百年一遇的洪峰流量為92 600m3·s-1[30]，由于模擬的時間為夏季豐水期，模型區(qū)域中長江的西邊界位于大通站下游約400km處的江陰，綜合考慮多方面因素并進行試驗對比前人的結(jié)果，選取73 000m3·s-1作為模式的長江徑流量輸入；錢塘江屬于山溪型河流，其流量的變率比較大，徑流量比長江的小很多，并且前人的研究表明錢塘江對杭州灣頂水交換的影響不明顯[31]，因此本模型不予考慮。

2.2 粒子追蹤模型的設(shè)置

為了研究杭州灣水交換的情況，將121.8300°E作為灣內(nèi)外的界限，在區(qū)域28.452 5°N～32.497 5°N，120.002 5°E～123.497 5°E內(nèi)(見圖2中的黑色虛框)均勻布放粒子，其間距為0.005(°)×0.005(°)，然后利用模型得到的流場對這些粒子進行追蹤。

(A：吳淞、B：高橋、C：中浚、D：蘆潮、E：綠華山、F：灘滸(圓點)；C1～C4為測流站的位置(菱形)；黑色虛線框為追蹤粒子布放區(qū)域。A：Wusong, B：Gaoqiao, C：Zhongjun, D：Luchao, E：Lvhuashan, F：Tanhu, (dots). C1～C4 are current stations (diamonds). Dashed box stands for the area where particles are released.)

圖2驗潮站及粒子布放區(qū)域

Fig.2Tidal stations and released area of particles

2.3 模型驗證

本文選取位于長江口和杭州灣附近的6個驗潮站及4個測流站進行天文潮潮位和流速的驗證，站位點見圖2，其中圓點(A～F)為潮位站，菱形(C1～C4)為潮流測站。

2.3.1 水位驗證圖3為2002-09-01—30日模擬的潮位與潮汐表中潮位數(shù)據(jù)的對比，可以看出模擬的潮位與潮汐表中的比較符合，變化趨勢較為一致。對潮位站的潮汐表水位數(shù)據(jù)和模擬得到的水位數(shù)據(jù)調(diào)和分析，得到調(diào)和常數(shù)的對比結(jié)果見表1。從表中可以看出4個分潮振幅的絕對誤差大多在10cm以下，最大不超過15cm；遲角的絕對誤差控制在20°以內(nèi)。

計算臺風情景時，風場和氣壓場是由Jelesnianski模型得到的，模式中其余設(shè)置均與只計算天文潮時一致，以此對2000-08-29—31日0012號派比安臺風過程進行了模擬(見圖4)。圖4左側(cè)一組圖為高橋、灘滸2個站位的風暴潮總水位與實測數(shù)據(jù)的對比，可以看出臺風過程中模式對水位的模擬也比較準確。圖4右側(cè)組圖給出了高橋和灘滸站實測增水與模擬增水的對比，模型較好地再現(xiàn)了水位變化過程。

2.3.2 流場驗證圖5，6分別為2009年6、7月大潮和小潮期間流速流向的驗證結(jié)果，實測數(shù)據(jù)源于姚炎明等[14]同期的觀測，可以看出流的轉(zhuǎn)向時刻與實測比較接近，模擬的大潮流速和流向與實測結(jié)果基本吻合，小潮期間模擬的與實測的流向比較吻合，流速略偏小。

對2000年8月28日漲潮和落潮潮流場進行模擬(見圖7)。從漲憩和落憩的潮流矢分布可以看出，東海潮波傳入杭州灣以后，潮流基本成往復狀態(tài)，漲潮流向西略偏北，基本與岸線平行，進入灣內(nèi)后受喇叭地形的約束逐漸轉(zhuǎn)向西南；落潮流向東偏南，流向與漲潮流相反。以上研究表明該模式能較好地反映天文潮以及臺風情況下所研究區(qū)域的動力情況。

2.3.3 拉格朗日余流場從2000年8月8日的高潮起釋放粒子，追蹤75h，得到杭州灣及其附近的拉格朗日余流場(見圖8)。從圖中可以看出，在南匯嘴東部，長江口余流分為向東流向外海的和貼著南匯嘴東灘流入杭州灣的兩支；流入杭州灣的順時針渦旋有利于長江水沿杭州灣北岸輸入灣內(nèi)；杭州灣北岸余流在金山附近南下后呈逆時針流向灣外，這是控制灣內(nèi)余流的主渦旋；在灣頂南岸存在順時針的渦旋，可將灣頂?shù)乃驗持胁枯斔筒㈦S逆時針主渦流向灣外；此外在舟山群島各島嶼附近有小的渦旋。這些結(jié)構(gòu)與前人[32-35]的研究基本一致。

3模擬結(jié)果與分析

為研究杭州灣內(nèi)與灣外的水體交換情況，本節(jié)將從灣內(nèi)外水體交換的來源、剩余函數(shù)和出流影響區(qū)域3個方面對杭州灣的水交換進行討論。

3.1 天文潮作用下的水交換

選取長江口(121.973 5°E，31.057 5°N)、杭州灣中部(121.872 5°E，30.272 5°N)以及杭州灣南部(121.8325°E，30.1075°N)3個粒子分別在高潮時和低潮時釋放，追蹤162h后的路徑見圖9。從圖中可以看出，由于杭州灣的潮流運動形式為顯著的往復流，且各分潮的旋轉(zhuǎn)率很小[28]，故粒子路徑成往復狀態(tài)；并且拉格朗日粒子追蹤的結(jié)果顯著依賴于初始釋放時刻，初始釋放時刻不同，粒子的運動路徑也不同。為了研究杭州灣內(nèi)外水體交換的來源情況，選擇在高潮時釋放粒子并追蹤72h。

3.1.1 杭州灣內(nèi)外水體交換的來源在追蹤過程中，只要灣外的粒子在某時刻進入到了灣內(nèi)，就將其初始位置標記在圖10(a)中，稱為灣外參與交換的水體來源；同樣只要灣內(nèi)的粒子到過灣外，其初始位置就標記在圖10(b)中，稱為灣內(nèi)參與交換的水體來源。從圖10可以看出在漲落潮流的作用下，杭州灣外可以參與交換的粒子初始位置集中在長江口南岸以及灣內(nèi)外分界線的附近，灣內(nèi)可以進行交換的粒子初始位置全集中在分界線附近。前人的研究也表明，長江口與杭州灣的水交換主要發(fā)生在近岸[12]，而杭州灣灣口中部區(qū)域是水交換最活躍的地方[11]，這與本文的結(jié)果一致。

3.1.2 杭州灣內(nèi)的剩余函數(shù)值為了進一步了解杭州灣內(nèi)外的水交換情況，計算了杭州灣的剩余函數(shù)時間序列。自2002年9月1日起模擬60d，得到灣內(nèi)水體剩余函數(shù)的時間序列(見圖11)。結(jié)果顯示，隨著時間的推移，剩余函數(shù)呈明顯的震蕩下降趨勢，2個月后，灣內(nèi)剩余函數(shù)降到48%。經(jīng)過計算，半交換時間約為46d；而當剩余函數(shù)為71%，即交換量達到29%時，約需要25.5d。

本文著重研究臺風過程中杭州灣水交換的情形，因此在1個臺風風暴潮過程內(nèi)研究初始釋放時間對剩余函數(shù)值的影響是有必要的。以2002年9月份為例，選取不同的大小潮高低潮組合做了8組敏感性試驗(見表2)，每個試驗計算72h。試驗表明無論是大潮期間還是小潮期間，高潮作為初始釋放時刻時灣內(nèi)粒子的剩余量都比低潮作為初始釋放時刻時要少：高潮時釋放72h后剩余函數(shù)值約為82%～87%，而從低潮釋放時剩余函數(shù)值約為93%～95%。所以選取不同潮時釋放粒子對水交換結(jié)果存在著較大影響，為統(tǒng)一表現(xiàn)灣內(nèi)外的水交換情況，本文的研究選擇高潮時作為初始時刻來釋放粒子。

Note：①Tidal moment；②Higher high tide；③Lower high tide；④Higher low tide；⑤Lower low tide；⑥Spring tide；⑦Neap tide

3.2 臺風影響下的水交換結(jié)果

自1949年到2014年，共有100多場臺風影響杭州灣地區(qū)，從中選取7條典型路徑的臺風進行研究(見圖12)，分別歸為南北向、東西向和轉(zhuǎn)向型的臺風。南北向路徑中，0807、0004、1109號臺風中心分別經(jīng)過灣西側(cè)，灣中部和灣東側(cè)；東西向路徑中0216、0008、1210臺風中心分別在灣南側(cè)、灣中部和灣北側(cè)；轉(zhuǎn)向型臺風選擇了1個個例即0515號臺風。為了研究臺風過程中的水交換，將初始時刻定義在臺風中心進入模型計算區(qū)域后第一個高潮時，此時開始釋放粒子。

3.2.1臺風影響下灣內(nèi)外交換的水體來源以0008

號臺風(見圖13左)為例，研究臺風影響下杭州灣內(nèi)外進行交換的水體來源。在臺風作用下，可以進入灣內(nèi)的粒子數(shù)明顯增多，范圍也擴展到崇明島附近(見圖13中)；能夠從灣內(nèi)流出的粒子量也顯著增加，范圍擴展到了杭州灣的南部區(qū)域(見圖13右)，這說明臺風的確大大地促進了杭州灣的水交換。

統(tǒng)計7場臺風對杭州灣水體交換來源情況的影響(見表3)，發(fā)現(xiàn)除0807號臺風以外，灣外參與交換的水體量都是增加的，增加量達43%～704%。灣內(nèi)參與交換的水體量也都是增加的，增加量達11%～161%。0807號臺風的情況將在后面的3.3節(jié)中進行詳細分析。

3.2.2 臺風影響下杭州灣內(nèi)剩余函數(shù)的變化首先以0008號臺風為例，研究天文潮和臺風情況下灣內(nèi)剩余函數(shù)的變化(見圖14)。從圖中可以看出，3天后灣內(nèi)水體的剩余量分別約89.2%和70.8%，即在天文潮和有臺風影響的情況下交換量分別為10.8%和29.2%，交換量提高了約1.7倍。在3.1節(jié)中提到，4個分潮作用下交換量達到29%約需要25.5d，與此處的3d相比，交換時間縮短了約7.5倍。研究其他幾場臺風，發(fā)現(xiàn)與天文潮時相比要交換相同的量所用時間可縮短約0.6～7.5倍。統(tǒng)計7場臺風過程中的剩余函數(shù)(見表4)，對比天文潮和臺風情況下的結(jié)果發(fā)現(xiàn)，在天文潮作用下，72h后水交換量約為10%～15%，而臺風到來時，交換量則可達11%～30%，交換量增加約0.2～1.7倍。

Note：①Tides；②Typhoons；③Increased multiples

3.2.3 灣內(nèi)水體對灣外的影響范圍為了進一步了解臺風對杭州灣水交換的影響，以0008臺風為例，模擬了杭州灣水體出流后的去向(見圖15(a))及其對灣外的影響面積(見圖15(b))。從圖15(a)可以看出，臺風情況下，灣內(nèi)水體影響范圍不僅僅局限在灣口中部區(qū)域，還擴展到南部舟山群島附近。圖15(b)中實線為每個時刻灣內(nèi)水體在灣外的影響面積；虛線為整個過程中灣內(nèi)水體能夠影響到灣外的所有區(qū)域的面積，稱為掃海面積。天文潮作用下各個時刻的影響范圍出現(xiàn)周期性震蕩，與漲落潮流有關(guān)。臺風作用下灣內(nèi)水體對灣外的影響范圍迅速增加。對比兩幅圖可知，臺風過程比天文潮過程影響面積約增加31%。

(圖中實線表示每個時刻對灣外的影響面積;虛線表示總影響面積。The solid lines stand for the affected acreage at every moment; The dotted lines stand for the total affected acreage.)

圖15天文潮(左)和臺風(右)情況下灣內(nèi)水體對灣外的影響范圍(a)和面積(b)

Fig.15The area affected by the water outflow from Hangzhou Bay under the tidal condition(left) and typhoon-tidal condition (right)

3.3 不同路徑的臺風對杭州灣水體來源與出流的影響

鄭東升等[36]研究風場對渤海灣水交換的影響時發(fā)現(xiàn)：風向決定風場對水交換是否有影響，風速大小決定影響的程度。不同路徑的臺風在杭州灣附近的風向不相同，因此不同路徑的臺風對水交換的影響也不相同。7場臺風過程除0004和0807號在計算過程中約為10～11級風，其余5場約為12～13級風。本文在這里只考慮風向?qū)λ粨Q的影響，風速的影響在以后的研究中討論。0008臺風過程中灣內(nèi)水體的來源與出流情況見圖13，其余6場臺風則如圖16所示。

在3.2節(jié)臺風對灣內(nèi)水體來源的影響研究中發(fā)現(xiàn)，在0807臺風的作用下，可以進入灣內(nèi)的水體量減少了約6%，與其路徑平行的0004號臺風增量也只有43%。分析兩場臺風發(fā)現(xiàn)，當臺風在臺灣島附近登陸北上或者不登陸繼續(xù)北上時，臺風路徑基本與杭州灣的走向垂直，臺風中心位于杭州灣左側(cè)，臺風移動過程中，風場由東南向順時針轉(zhuǎn)為西南向，抑制長江水北下，促進水從北部及灣口中部區(qū)域流出。從圖16(a)和(d)可以看出，0004和0807號臺風對長江水進入杭州灣有顯著的抑制作用，而促進灣內(nèi)水體從杭州灣北岸流出。

1109臺風也屬于南北向路徑的臺風，與另外2條南北向路徑臺風不同的是：其臺風中心位于杭州灣的東側(cè)，臺風過程中影響杭州灣的風場由東北向逆時針轉(zhuǎn)向為西向促進長江水進入杭州灣；在灣南部的風場為東南向，促進灣內(nèi)的水從南部流出(見圖16(e))。

3個東西向的臺風中，0008和0216臺風的中心位于灣的西側(cè)，在臺風移動的過程中，杭州灣附近的臺風由北向順時針轉(zhuǎn)為東北向，最終轉(zhuǎn)為東南向；而1210臺風的中心位于杭州灣右側(cè)，在杭州灣附近產(chǎn)生的風場由北向逆時針轉(zhuǎn)為西南。3場臺風產(chǎn)生的風暴余流場都與1109相似，灣北部東南向，灣南部西南向，因而灣內(nèi)水體的來源與流出區(qū)域也比較相似(見圖13，圖16(b),16(f))。

轉(zhuǎn)向型臺風0515在移動過程中，杭州灣附近的風場由東北順時針轉(zhuǎn)為西南，產(chǎn)生的風暴余流場比較混亂，但灣北部有向南的分量，促進長江水進入杭州灣，而對于能夠從灣內(nèi)出來的水體范圍無明顯規(guī)律(見圖16(c))。

綜上，雖然臺風在移動的過程中，風場的方向不盡相同，但是只有南北向從臺灣島附近登陸或者不登陸直接北上的臺風會促進杭州灣內(nèi)水體從杭州灣北部流出，抑制長江水進入杭州灣；其余路徑的臺風都促進長江水進入杭州灣，除轉(zhuǎn)向型路徑外，其余路徑還促進灣內(nèi)的水從南部流出，只是由于臺風強度的不同，影響的大小也不一樣。

4結(jié)論

在杭州灣及其臨近海域二維潮汐潮流模式的基礎(chǔ)上，結(jié)合拉格朗日粒子追蹤方法，模擬了不同臺風過程中杭州灣水體交換的情況，得出結(jié)論如下：

(1)大多臺風能夠促進灣外水進入灣內(nèi)，灣外參與交換的水體增加量因臺風強度和路徑的不同最少是原來的0.43倍，最大可以達到7.04倍；臺風也會促進灣內(nèi)水體的出流，灣內(nèi)參與交換的水體增加量至少是原來的0.11倍，最大可達到1.61倍。

(2)當一場臺風影響杭州灣時，72h的水交換量可以提高0.2～1.7倍。與天文潮作用相比，要交換相同的量，臺風情景時所需要的時間可縮短約0.6～7.5倍。臺風作用下杭州灣內(nèi)水體對灣外的影響區(qū)域發(fā)生改變，影響面積也顯著增大。

(3)不同路徑的臺風對杭州灣參與交換的水體來源及灣內(nèi)水體的出流影響區(qū)域不同，南北向從臺灣島附近登陸或者不登陸直接北上的臺風促進杭州灣內(nèi)的水從灣北部流出，抑制長江水進入杭州灣；轉(zhuǎn)向型路徑促進長江水進入杭州灣，對灣內(nèi)水體出流的區(qū)域影響不明顯；其余路徑臺風都促進長江徑流進入杭州灣并促進灣內(nèi)的水從北部流出。

參考文獻：

[1]高生泉, 陳建芳, 金海燕, 等. 杭州灣及鄰近水域營養(yǎng)鹽的時空分布與富營養(yǎng)化特征[J]. 海洋學研究, 2011(3): 36-47.

Gao S, Chen J, Jin H, et al. Characteristics of nutrients and eutrophicationin the Hangzhou Bay and its adjacent waters[J]. Journal of Marine Sciences, 2011(3): 36-47.

[2]張健, 施青松, 鄔翱宇, 等. 杭州灣豐水期主要污染因子的分布變化及成因[J]. 東海海洋， 2002(4): 35-41.

Zhang J, Shi Q, Wu A, et al. Distribution characteristic analysis of main pollution factor in rainy season in the Hangzhou Bay [J]. Donghai Marine Science, 2002(4): 35-41.

[3]國家海洋局. 2014年中國海洋環(huán)境狀況公報[R]. 北京: 國家海洋局, 2015.

State Oceanic Administration People's Republic of China. Bulletin of Marine Environmental Status of China for the year of 2014[R]. Beijing: State Oceanic Administration People's Republic of China, 2015.

[4]姚炎明, 彭輝, 杜雅杰, 等. 象山港分區(qū)水交換數(shù)值研究[J]. 海洋學報(中文版), 2014(1): 126-130.

Yao Y, Peng H, Du Y, et al. Numerical study on water exchange in the Xiangshangang Bay[J]. Acta Oceanologica Sinica(in Chinese), 2014(1): 126-130.

[5]高抒, 謝欽春. 狹長形海灣與外海水體交換的一個物理模型[J]. 海洋通報， 1991(3): 1-9.

Gao S, Xie Q. Physical model for water exchange process in Xiangshan Bay[J]. Marine Science Bulletin, 1991(3): 1-9.

[6]Kitheka J U. Coastal tidally-driven circulation and the role of water exchange in the linkage between tropical coastal ecosystems[J]. Estuarine, Coastal and Shelf Science, 1997, 45(2): 177-187.

[7]魏皓, 田恬, 周鋒, 等. 渤海水交換的數(shù)值研究-水質(zhì)模型對半交換時間的模擬[J]. 青島海洋大學學報, 2002(4): 519-525.

Wei H, Tian T, Zhou F, et al. Numerical study on the water exchange of the Bohai Sea: simulation of the half -life time by dispersion model[J]. Journal of Ocean University of Qingdao, 2002(4): 519-525.

[8]Luff R P T. Calculation of water exchange times in the ICES-boxes with an Eulerian dispersion model using a half-life time approach[J]. Deutsche Hydrografische Zeitschrift, 1995, 47(4): 287-299.

[9]Signell R P, Butman B. Modeling tidal exchange and dispersion in Boston Harbor[J]. Journal of Geophysical Research: Oceans, 1992, 97(C10): 15591-15606.

[10]Bilgili A, Proehl J A, Lynch D R, et al. Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study[J]. Estuarine, Coastal and Shelf Science, 2005, 65(4): 607-624.

[11]張淑珍, 俞光耀, 王化桐. 海灣物理自凈能力分析及水質(zhì)預測實例[J]. 山東海洋學院學報, 1988(S1): 53-72.

Zhang S, Yu G, Wang H. Analysis of physical self-purification capacity and prediction of water quality in the bay[J]. Journal of Shandong College of Oceanology, 1988(S1): 53-72.

[12]劉新成, 盧永金, 潘麗紅, 等. 長江口和杭州灣潮流數(shù)值模擬及水體交換的定量研究[J]. 水動力學研究與進展, 2006(2): 171-180.

Liu X, Lu Y, Pan L, et al. Tidal current numerical simulating and water exchange research in Yangtze Estuary and Hangzhou Bay[J]. Journal of Hydrodynamics, 2006(2): 171-180.

[13]Kong J, Song Z, Xia Y, et al. Characteristics of water and sediment exchange Between Yangtze Estuary and Hangzhou Bay[J]. China Ocean Engineering, 2007, 21(2): 255-266.

[14]姚炎明, 朱斌, 李佳. 錢塘江河口水體沖洗時間的計算[J]. 浙江大學學報(理學版), 2012(6): 711-716.

Yao Y, Zhu B, Li J. The computation of flushing time in Qiantang Estuary[J]. Journal of Zhejiang University(Science Edition), 2012(6): 711-716.

[15]Wang Y, Shen J, He Q, et al. Seasonal variations of transport time of freshwater exchanges between Changjiang Estuary and its adjacent regions[J]. Estuarine, Coastal and Shelf Science, 2015, 157: 109-119.

[16]陳春華, 侍茂崇, 孫士才, 等. 瓊州海峽和?？跒撑_風引起的水交換研究[J]. 海洋環(huán)境科學, 1997(3): 9-11.

Chen C, Shi M, Sun S, et al. Study on the water exchange by typhoon in Qiongzhou Strait and Haikou Bay[J]. Marine Environmental Science, 1997(3): 9-11.

[17]Ying M, Zhang W, Yu H, et al. An overview of the China meteorological administration tropical cyclone database[J]. Journal of Atmospheric and Oceanic Technology, 2014, 31(2): 287-301.

[18]Luettich R A, Westerink J J. ADCIRC User Manual: A (Parallel) Advanced Circulation Model for Oceanic, Coastal and Estuarine Waters[R]. Morehead City: University of North Corolina at Chapel Hill Institute of Marine Sciences, 2006.

[19]郭洪琳. 甌江口風暴潮漫灘數(shù)值預報研究[D]. 北京： 國家海洋環(huán)境預報研究中心, 2011.

Guo H. A numerical prediction study of storm surge and flood plain on the delta of Oujiang[D]. Beijing: National Marine Environmental Forecasting Center, 2011.

[20]馮興如. 重點區(qū)臺風風暴潮數(shù)值預報與災害評估研究[D]. 青島: 中國科學院研究生院(海洋研究所), 2012.

Feng X. Study on the hurricane storm surge disaster assessment and numerical prediction technique[D]. Qingdao: Institute of Oceanology, Chinese Academy Sciences. 2012.

[21]趙長進. 長江口及其鄰近海區(qū)無結(jié)構(gòu)網(wǎng)格風暴潮模式的建立與應用[D]. 上海: 華東師范大學, 2014.

Zhao C. Development and analysis of the unstructured numerical model for storm surge in the Changjiang Estuary[D]. Shanghai: East China Normal University. 2014.

[22]Feng J, Jiang W, Bian C. Numerieal prediction of storm surge in the Qingdao area under the impact of climate change[J]. Journal of Ocean University of China, 2014, 13(4): 539-551.

[23]Luettich R A, Westerink J J, Scheflher N W. An Advanced Three-Dimensional Circulation Model for Shelves, Coasts, and Estuaries. Report 1. Theory and Methodology of ADCIRC-2DDI and ADCIRC-3DL[R]. Washington, DC: Coastal Engineering Research Center Vicksburg Ms, 1992.

[24]Atkinson G D, Holliday C R. Tropical cyclone minimum sea level pressure/maximum sustained wind relationship for the western north Pacific[J]. Monthly Weather Review, 1977, 105(4): 421-427.

[25]潘麗紅. 臺風條件下上海地區(qū)典型海堤防御能力評價研究[D]. 上海: 華東師范大學, 2011.

Pan L. Assessment of defense capability of sea dikes in Shanghai for typhoon[D]. Shanghai: East China Normal University. 2011.

[26]Jelesnianski C P. A numerical calculation of storm tides induced by a tropical storm impinging on a continental shelf[J]. Monthly Weather Review, 1965(93): 343-358.

[27]Takeoka H. Fundamental concepts of exchange and transport time scales in a coastal sea[J]. Continental Shelf Research, 1984, 3(3): 311-326.

[28]谷國傳, 李身鐸, 胡方西. 杭州灣北部近岸水域水文泥沙特性[J]. 東海海洋, 1987, 5(4): 15-24.

Gu G, Li S, Hu F. The hydrologic characteristics and silts transportation in nearshore waters of northern Hangzhou Bay[J]. Donghai Marine Science, 1987, 5(4): 15-24.

[29]Maritin P J, Smith S R, Posey P G, et al. Use of the Oregon State University Tidal Inversion Software (OTIS) to Generate Improved Tidal Prediction in the East-Asian Seas[R].Arlington: Naval Research Laboratory Oceanography Division Stemis Space Center, 2009.

[30]黃華. 長江口及杭州灣風暴潮三維數(shù)值模擬[D]. 上海: 華東師范大學, 2006.

Huang H. 3-D numerical simulation of storm surges in the Changjiang Estuary and Hangzhou Bay[D]. Shanghai: East China Normal University, 2006.

[31]郜會彩, 李若華, 朱軍政. 強潮河口水體交換能力敏感因素分析[J]. 浙江水利科技, 2013(5): 5-7.

Gao H, Li R, Zhu J. Sensitive factor analysis of water exchange capacity in strong tide estuary[J]. Zhejiang Hydrotechnics, 2013(5): 5-7.

[32]李身鐸, 孫衛(wèi)陽. 杭州灣潮致余流數(shù)值研究[J]. 海洋與湖沼, 1995(3): 254-261.

Li S, Sun W. Numerical modeling residual currents in Hangzhou Bay[J]. Oceanologia et Limnologia Sinica, 1995(3): 254-261.

[33]朱首賢, 丁平興, 史峰巖, 等. 杭州灣、長江口余流及其物質(zhì)輸運作用的模擬研究Ⅱ冬季余流及其對物質(zhì)的運輸作用[J]. 海洋學報(中文版), 2000(6): 1-12.

Zhu S, Ding P, Shi F, et al. Numerical study on residual current and its effect on mass transport in the Hangzhou Bay and the Changjiang EstuaryⅡ. The residual current and its effect on mass transport in winter [J]. Acta Oceanologica Sinica (in Chinese), 2000(6): 1-12.

[34]孔俊. 長江口、杭州灣水沙交換特性初步研究[D]. 南京： 河海大學, 2005.

Kong J. Preliminary study on the flow and sediment flux exchange characteristics between Changjiang Estuary and Hangzhou Bay[D]. Nanjing: Hohai University, 2005.

[35]鄒濤, 高會旺, 孫文心, 等. 長江口、杭州灣及其鄰近海區(qū)Lagrange環(huán)流的數(shù)值模擬研究Ⅰ——正壓環(huán)流[J]. 中國海洋大學學報(自然科學版), 2009(1): 153-159.

Zou T, Gao H, Sun W, et al. Numerical simulation of Lagrange residual current in the Changjiang Estuary and their adjacent sea Ⅰ: Barotropic circulation[J]. Periodical of Ocean University of China, 2009(1): 153-159.

[36]鄭東升, 孫健, 陶建華. 一種基于能量形式的風場數(shù)據(jù)統(tǒng)計方法及其在渤海灣水交換方面的應用[J]. 海洋環(huán)境科學, 2013(2): 280-283.

Zheng D, Sun J, Tao J. A statistical method of wind field based on the form of energy and its application in Bohai Bay water exchange[J]. Marine Environmental Science, 2013(2): 280-283.

責任編輯龐旻

Quantitative Study on Water Exchange Affected by Typhoons in the Hangzhou Bay

XU Ling-Jing1, YANG Bo1, JIANG Wen-Sheng2

(Ocean University of China, 1.College of Oceanic and Atmospheric Sciences; 2.The Key Laboratory of Marine Environment and Ecology, Ministry of Education, Qingdao 266100, China)

Abstract：A 2D numerical model of storm surge is established based on ADCIRC model in the Hangzhou Bay. The Particle Tracking Method is used to study the water exchange according to the model results under the tidal condition and typhoon-tidal condition respectively. The influences of typhoons on water exchange in Hangzhou Bay are quantitatively analyzed. Results show that when the bay is affected by strong typhoons, the outflow water area inside the bay can increase 1.6 times compared to conditions that only forced by tide, while the inflow area outside the bay can increase 7 times. And the affected area of outflow water can increase 31%, the water exchange rate can increase 1.7 times compared with those under tidal conditions. It also reveals that typhoons with different paths have different influences on the water exchange of Hangzhou Bay.

Key words:Hangzhou Bay; water exchange; typhoon; numerical simulation

基金項目：? 海洋公益性行業(yè)科研專項經(jīng)費項目(201205017-3)；國家自然科學基金委員會-山東省人民政府聯(lián)合資助海洋科學研究中心項目(U1406401)資助

收稿日期：2015-08-05；

修訂日期：2015-11-12

作者簡介：許靈靜(1991-)，女，碩士生。E-mail：lingjingxu2016@126.com ??通訊作者：E-mail: yang.bo@ouc.edu.cn

中圖法分類號：P731.27

文獻標志碼：A

文章編號：1672-5174(2016)06-029-13

DOI：10.16441/j.cnki.hdxb.20150280

Supported by Public Science and Technology Research Funds Projects of Ocean (201205017-3)； NSFC-Shandong Joint Fund for Marine Science Research Centers (U1406401)