黃華聰,賈 曉,2,路川藤

(1.上海河口海岸科學(xué)研究中心,上海 201201; 2.上海交通大學(xué)水動(dòng)力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200240;3.南京水利科學(xué)研究院河流海岸研究所,江蘇 南京 210029)

影響長(zhǎng)江口深水航道驟淤的非常態(tài)天氣過(guò)程Ⅰ:臺(tái)風(fēng)的路徑特征及數(shù)值驗(yàn)證

黃華聰1,賈 曉1,2,路川藤3

(1.上海河口海岸科學(xué)研究中心,上海 201201; 2.上海交通大學(xué)水動(dòng)力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200240;3.南京水利科學(xué)研究院河流海岸研究所,江蘇 南京 210029)

為歸納長(zhǎng)江口深水航道臺(tái)風(fēng)期驟淤的發(fā)生規(guī)律及特征,分析了發(fā)生驟淤時(shí)刻的氣象條件與對(duì)應(yīng)的波浪條件。研究發(fā)現(xiàn),牛皮礁站的波能與驟淤具有較好的相關(guān)性;從臺(tái)風(fēng)路徑上分析,長(zhǎng)江口東側(cè)過(guò)境臺(tái)風(fēng)對(duì)航道的驟淤影響顯著。結(jié)合歷史臺(tái)風(fēng)路徑,選取3個(gè)典型路徑的臺(tái)風(fēng),選擇藤田-高橋圓形經(jīng)驗(yàn)風(fēng)場(chǎng)和CFSR(climate forecast system reanalysis)風(fēng)場(chǎng)的混合風(fēng)場(chǎng)復(fù)演了臺(tái)風(fēng)場(chǎng),然后采用SWAN模型模擬了不同路徑臺(tái)風(fēng)期間的波況,最后以牛皮礁站的淺水波能流為判別參數(shù),分析不同路徑臺(tái)風(fēng)對(duì)長(zhǎng)江口深水航道驟淤的影響。研究表明長(zhǎng)江口東側(cè)過(guò)境的臺(tái)風(fēng)是較易產(chǎn)生較大波能并進(jìn)一步誘發(fā)驟淤的典型臺(tái)風(fēng)路徑,這一分析結(jié)果與2010年以來(lái)的驟淤實(shí)測(cè)臺(tái)風(fēng)路徑結(jié)果吻合。

長(zhǎng)江口深水航道;驟淤;臺(tái)風(fēng)路徑;藤田-高橋經(jīng)驗(yàn)風(fēng)場(chǎng);CFSR風(fēng)場(chǎng);SWAN模型;波能

驟淤是港口航道工程十分關(guān)心的水沙事件,其因淤積量大且淤積時(shí)間短,給疏浚造成很大的壓力,甚至?xí)绊懙酱暗陌踩ê?。長(zhǎng)江口深水航道位于長(zhǎng)江口攔門(mén)沙河段,其含沙量高,鹽度適于絮凝淤積,且不時(shí)遭受非常態(tài)天氣,諸如臺(tái)風(fēng)、寒潮的影響。在長(zhǎng)江口深水航道二期工程以來(lái),屢次觀測(cè)到驟淤現(xiàn)象[1-2]。

已有多位學(xué)者從不同方面闡述了長(zhǎng)江口深水航道的淤積機(jī)制,比如研究最大渾濁帶的絮凝因素[3]、長(zhǎng)江口整合的動(dòng)力場(chǎng)[4]及南導(dǎo)堤越堤水沙[5]等對(duì)淤積的潛在影響,但就長(zhǎng)江口的驟淤問(wèn)題,研究較少。目前驟淤計(jì)算方法主要源于黃驊港粉砂質(zhì)海岸航道的淤積計(jì)算[6-7],對(duì)于長(zhǎng)江口大型河口區(qū)域的航道淤積適用性不強(qiáng)。

筆者在現(xiàn)有的驟淤量統(tǒng)計(jì)基礎(chǔ)上,從氣象要素和波浪要素與驟淤量的表征相關(guān)性入手展開(kāi)分析,通過(guò)構(gòu)建風(fēng)場(chǎng)模型和SWAN模型,模擬特定臺(tái)風(fēng)過(guò)程下長(zhǎng)江口深水航道口門(mén)處的波能,研究發(fā)生驟淤時(shí)刻的氣象條件與對(duì)應(yīng)的波浪條件。

1 長(zhǎng)江口深水航道的驟淤概述

圖1 長(zhǎng)江口河勢(shì)及周邊水文測(cè)站位置Fig.1 River map of Yangtze Estuary and the location of surrounding hydrological stations

長(zhǎng)江口12.5 m深水航道橫貫長(zhǎng)江口攔門(mén)沙,周邊河勢(shì)及長(zhǎng)江口主要水文站點(diǎn)見(jiàn)圖1。航道自2010年3月貫通以來(lái),在發(fā)揮巨大經(jīng)濟(jì)效益和社會(huì)效益的同時(shí),回淤量時(shí)空分布高度集中的問(wèn)題也非常突出,每年需投入大量的維護(hù)疏浚費(fèi)用。長(zhǎng)江口12.5 m深水航道回淤包括常態(tài)回淤和由臺(tái)風(fēng)、寒潮大風(fēng)等惡劣自然因素引起的驟淤。臺(tái)風(fēng)、寒潮引起的航道驟淤由于作用時(shí)間短、淤積范圍集中,對(duì)航道水深的維護(hù)造成極大困難。

驟淤本質(zhì)是特定條件下產(chǎn)生的快速淤積過(guò)程。實(shí)際上巨型河口深水航道內(nèi)的驟淤,因?yàn)槠渌锢磉^(guò)程(典型如潮流)的作用以及疏浚作業(yè)的持續(xù)性,無(wú)法得到嚴(yán)格準(zhǔn)確的驟淤量。根據(jù)前后對(duì)比法得到的驟淤量,反映的是一段時(shí)間內(nèi)淤積的增加,可以準(zhǔn)確反映驟淤的發(fā)生與否。

據(jù)長(zhǎng)江口深水航道管理局統(tǒng)計(jì),2010—2014年航道回淤總量年均約為8.6×107m3,其中常態(tài)回淤年均約為7.6×107m3,驟淤量年均1.0×107m3左右。

2 波能與驟淤的關(guān)系

河流及海岸中的泥沙起懸、落淤能夠受到臺(tái)風(fēng)的影響,主要是以風(fēng)浪為中間媒介,將風(fēng)能傳遞到水體。研究表明,在波浪潮流作用下,泥沙極易揚(yáng)起,風(fēng)浪停息后,泥沙又極易沉積,造成港口航道驟淤[8-9]。有部分學(xué)者嘗試分析波浪動(dòng)力與長(zhǎng)江口深水航道淤積的相關(guān)性[10]。近年來(lái)波能流(或稱波功率,簡(jiǎn)稱為波能)作為波浪綜合作用指標(biāo)參數(shù),通常被用來(lái)直接統(tǒng)計(jì)其與淤積的表征相關(guān)性。劉猛等[11]采用深水波能公式計(jì)算了長(zhǎng)江口波能與淤積的相關(guān)性。本文以淺水波能為特征參數(shù),分析驟淤與波能的相關(guān)性。

波能是波浪在傳播方向上單位時(shí)間內(nèi)通過(guò)單位波峰寬度上的能量。以深水規(guī)則波為例,其等于單寬波峰線長(zhǎng)度的波浪作用在與波峰線平行的一個(gè)豎直面上各點(diǎn)的力與波浪的水質(zhì)點(diǎn)通過(guò)此豎直面上各相應(yīng)點(diǎn)的流速之積的總和在一個(gè)波周期內(nèi)的平均值。

[12],單個(gè)淺水規(guī)則波的波能:

(1)

式中:ρ——水體密度;g——重力加速度;H——波高;T——周期;k——波數(shù);h——水深。

式(1)表達(dá)的是時(shí)域中,單位時(shí)間內(nèi)特定波高周期下的波能,而在實(shí)際測(cè)量中不可能測(cè)量時(shí)域上每一波浪的尺度參數(shù)。以長(zhǎng)江口3個(gè)水文站為例(圖1),能夠表達(dá)波況的參數(shù)僅有最大波高Hmax、有效波高Hs、有效周期Ts、平均波高Hm及平均周期Tm等一系列頻域參數(shù)。

為了合理評(píng)估牛皮礁站水域的波能,將頻域上的波要素按照相應(yīng)的時(shí)間比例拓展在時(shí)域上。即將特征時(shí)間段(測(cè)站為1 h間隔的散點(diǎn)數(shù)據(jù))分成2段,其中1/3時(shí)間段的波能計(jì)算采用Hs和Ts,剩余時(shí)間的波能計(jì)算采用Hm和Tm,然后將各自得到的波能乘以相應(yīng)的時(shí)間,得到給定頻域波要素下特征時(shí)間段內(nèi)的總波能。以長(zhǎng)江口水文站為例,1 h的總波能為

Q=1 200Ps+2 400Pm

(2)

式中:Ps、Pm——依據(jù)有效波參數(shù)、平均波參數(shù)計(jì)算的波能。

牛皮礁站(圖1)位于長(zhǎng)江口深水航道口門(mén)內(nèi)約2 km處,能夠較好地代表輸入北槽的波能。以發(fā)生臺(tái)風(fēng)頻次較高的2012年為例,牛皮礁站實(shí)測(cè)數(shù)據(jù)計(jì)算的波能與發(fā)生驟淤的時(shí)間區(qū)間的對(duì)應(yīng)關(guān)系見(jiàn)圖2。由圖2可以看出,驟淤時(shí)間區(qū)間內(nèi),牛皮礁站的實(shí)測(cè)波能均較大。如果以牛皮礁站Hs>1.0 m為一次臺(tái)風(fēng)過(guò)程的統(tǒng)計(jì)時(shí)段,對(duì)該時(shí)段內(nèi)的波能進(jìn)行積分,就可以得到單次臺(tái)風(fēng)過(guò)程單位面積水體的總波能。

由歷年的驟淤量和波能相關(guān)關(guān)系(圖3)可以看出,驟淤和波能存在很好的相關(guān)性。該現(xiàn)象在以前諸多研究中得到過(guò)一定程度的分析,主要認(rèn)為是波能的增加引起海床孔隙水壓變化導(dǎo)致泥沙起懸增多等諸多因素導(dǎo)致。

圖2 牛皮礁站波能與驟淤時(shí)間區(qū)間的對(duì)應(yīng)關(guān)系Fig.2 Correspondenc between wave energy and time intervals of sudden siltation in Niupijiao Station

圖3 牛皮礁站波能與驟淤量的相關(guān)關(guān)系Fig.3 Correlation analysis between wave energy and sudden siltation in Niupijiao Station

類(lèi)型名稱氣壓/hPa十級(jí)風(fēng)圈半徑/km牛皮礁最大風(fēng)速/(m·s-1)最大風(fēng)速時(shí)相應(yīng)的風(fēng)向/(°)東側(cè)過(guò)境布拉萬(wàn)95515019.2350南側(cè)登陸?？?5018016.6101南側(cè)登陸西側(cè)過(guò)境麥德姆99015.6173

圖4 3種典型臺(tái)風(fēng)的路徑Fig.4 Paths of three typical typhoons

3 典型路徑臺(tái)風(fēng)及風(fēng)場(chǎng)驗(yàn)證

為了分析何種類(lèi)型的臺(tái)風(fēng)能夠使輸入北槽的波

能明顯增加,選取3種典型路徑的臺(tái)風(fēng),采用圓形梯度風(fēng)場(chǎng)耦合CFSR背景風(fēng)場(chǎng)的混合風(fēng)場(chǎng)數(shù)據(jù),驗(yàn)證了長(zhǎng)江口周邊的風(fēng)場(chǎng)。

3.13種典型路徑臺(tái)風(fēng)

根據(jù)臺(tái)風(fēng)的登陸位置和影響范圍,在2010年以后發(fā)生的臺(tái)風(fēng)內(nèi),選取影響長(zhǎng)江口水域的3種典型路徑的臺(tái)風(fēng)過(guò)程(表1),其路徑(參考浙江臺(tái)風(fēng)網(wǎng))見(jiàn)圖4。

3.2圓形風(fēng)場(chǎng)模式概述

進(jìn)行氣旋風(fēng)場(chǎng)模擬時(shí),首先需選擇合適的臺(tái)風(fēng)氣壓場(chǎng)模式。本文中采用國(guó)家海洋環(huán)境預(yù)報(bào)中心的臺(tái)風(fēng)暴潮模式[13],當(dāng)距臺(tái)風(fēng)中心距離r<2R0(R0為臺(tái)風(fēng)最大風(fēng)速半徑)時(shí)采用藤田模式,當(dāng)r≥2R0時(shí)采用高橋模式計(jì)算氣壓場(chǎng),移行風(fēng)速采用Veno Takeo公式。實(shí)踐證明,在合理選擇最大風(fēng)速半徑的前提下,這種方法基本能夠模擬海面上的氣旋風(fēng)場(chǎng)。

3.3CFSR風(fēng)場(chǎng)概述

CFSR(NCEP climate forecast system reanalysis)是由NCEP(national centers for environmental prediction)設(shè)立的全球高精度氣象后分析系統(tǒng),其可以給出風(fēng)浪計(jì)算所關(guān)心的高地理精度(0.3°×0.3°)的海面氣象預(yù)估場(chǎng)。自2011年以來(lái),CFSR 已經(jīng)升級(jí)至CFSV2(climate forecast system version 2),其可以認(rèn)為是CFSR的拓展版本,可以給出精度達(dá)0.2°×0.2°的氣象數(shù)據(jù)[14-15]。

CFSR氣象數(shù)據(jù)實(shí)際經(jīng)過(guò)諸多氣象衛(wèi)星數(shù)據(jù)的同化處理,大范圍的數(shù)據(jù)可信,并且其精度基本可以滿足圓形風(fēng)場(chǎng)外圍的風(fēng)場(chǎng)同化[16],在相關(guān)研究中已經(jīng)有較廣泛應(yīng)用(歷史文獻(xiàn)通常引用了其歷史資料,稱之為NCEP背景風(fēng)場(chǎng))。

3.4風(fēng)場(chǎng)驗(yàn)證

將圓形風(fēng)場(chǎng)和CFSR背景風(fēng)場(chǎng)進(jìn)行耦合,可以得到適用的混合風(fēng)場(chǎng)。以牛皮礁站和浦東機(jī)場(chǎng)(數(shù)據(jù)源自weather underground)為例,3個(gè)臺(tái)風(fēng)的合成風(fēng)場(chǎng)與實(shí)測(cè)風(fēng)場(chǎng)的數(shù)據(jù)對(duì)比見(jiàn)圖5和圖6。其中牛皮礁站點(diǎn)風(fēng)速經(jīng)過(guò)高度折算,自實(shí)測(cè)的25 m高程風(fēng)速折算至標(biāo)準(zhǔn)10 m高程,高度折減系數(shù)0.88[17]。由圖5和圖6看出, 2個(gè)站點(diǎn)的風(fēng)速、風(fēng)向過(guò)程模擬較好,可以以此為輸入條件計(jì)算長(zhǎng)江口水域的波浪場(chǎng)。

圖5 牛皮礁站風(fēng)場(chǎng)驗(yàn)證Fig.5 Validations of wind fields in NPJ Sites

圖6 浦東機(jī)場(chǎng)站風(fēng)場(chǎng)驗(yàn)證Fig.6 Validations of wind fields in PVG Sites

4 海浪數(shù)值模擬及波能統(tǒng)計(jì)

采用SWAN模型計(jì)算長(zhǎng)江口深水航道口門(mén)處的分向波能,并分析影響顯著的臺(tái)風(fēng)路徑。

4.1模型簡(jiǎn)介

SWAN模型是時(shí)下最新的第3代海浪數(shù)學(xué)模型[18-20],其在WAM模型基礎(chǔ)上改進(jìn)了算法并且加入了額外的淺水控制項(xiàng)。SWAN模型對(duì)不同水深的適應(yīng)性較好,且較為全面地考慮了包括能量輸入、損耗和非線性相互作用等在內(nèi)的一系列源匯項(xiàng);模型以不規(guī)則譜型的方向譜表示真實(shí)海浪的隨機(jī)性特征,模擬結(jié)果更接近真實(shí)海浪。SWAN模型采用動(dòng)譜平衡方程描述風(fēng)浪生成及其在近岸區(qū)的演化過(guò)程。

SWAN模型在對(duì)波浪成長(zhǎng)的譜型沒(méi)有任何先驗(yàn)限制的條件下,求解了動(dòng)譜能量平衡方程,且在深水域使用了相同的源項(xiàng),包括風(fēng)能輸入、海浪成長(zhǎng)、白帽耗散及四波相互作用;而在淺水域中,源項(xiàng)補(bǔ)充了包括底摩阻損耗、三波相互作用和水深引起的破碎等。

SWAN模型采用全隱式差分格式,且可以調(diào)節(jié)系數(shù)的取值決定譜空間的差分格式是偏于迎風(fēng)格式或中心格式,即決定在譜空間和方向空間的數(shù)據(jù)離散格式和收斂性強(qiáng)弱,在淺水海域的模擬中更穩(wěn)定且節(jié)省時(shí)間。

4.2海浪模型驗(yàn)證

以長(zhǎng)江口2個(gè)典型水文站的波浪過(guò)程為例,3個(gè)典型臺(tái)風(fēng)路徑下,牛皮礁站和長(zhǎng)江口站的波高和周期驗(yàn)證情況見(jiàn)圖7和圖8。

由圖7和圖8可見(jiàn),SWAN模型計(jì)算的臺(tái)風(fēng)浪過(guò)程與實(shí)測(cè)基本吻合,其計(jì)算結(jié)果能夠較好地代表長(zhǎng)江口水域在臺(tái)風(fēng)期間的波況。

圖7 牛皮礁站波況過(guò)程驗(yàn)證Fig.7 Validations of wave conditions in NPJ Sites

圖8 長(zhǎng)江口站波況過(guò)程驗(yàn)證Fig.8 Validations of wave conditions in CJK Sites

4.3波能時(shí)序分布

根據(jù)數(shù)值計(jì)算,3個(gè)臺(tái)風(fēng)情況下,長(zhǎng)江口站波能與牛皮礁站波能沿時(shí)間序列的分布見(jiàn)圖9。

圖9 牛皮礁站和長(zhǎng)江口站波能變化特征Fig.9 Characteristics of wave energy changes during three typical typhoons(NPJ and CJK Sites)

圖10 引起長(zhǎng)江口明顯驟淤的臺(tái)風(fēng)路徑(2010年后)Fig.10 Paths of typhoons inducing sudden siltation in Yangzi Estuary (after 2010)

由圖9可以看出,盡管“?？迸_(tái)風(fēng)與“布拉萬(wàn)”臺(tái)風(fēng)在長(zhǎng)江口氣象站的風(fēng)速相近,均為18 m/s左右,但其在長(zhǎng)江口形成的波能較小。上述現(xiàn)象的原因是南側(cè)登陸臺(tái)風(fēng)在嵊泗東部海域波高較大,但受東南側(cè)嵊泗列島的遮蔽,其產(chǎn)生的SE向浪在向長(zhǎng)江口水域傳播過(guò)程中波能衰減明顯。

根據(jù)上述研究,強(qiáng)度相近的臺(tái)風(fēng),路徑在東側(cè)過(guò)境的臺(tái)風(fēng)更容易在長(zhǎng)江口產(chǎn)生較大的波能并進(jìn)一步誘導(dǎo)驟淤的產(chǎn)生。根據(jù)長(zhǎng)江口航道管理局統(tǒng)計(jì)的數(shù)據(jù),自2010年以來(lái),引發(fā)驟淤的臺(tái)風(fēng)多數(shù)為東側(cè)過(guò)境臺(tái)風(fēng),部分由東側(cè)過(guò)境臺(tái)風(fēng)與其他路徑臺(tái)風(fēng)疊加影響,此與本文的研究結(jié)論相吻合,見(jiàn)圖10。

5 結(jié) 語(yǔ)

a. 北槽口內(nèi)牛皮礁站的波能與驟淤存在相關(guān)性。

b. 基于經(jīng)驗(yàn)風(fēng)場(chǎng)和CFSR風(fēng)場(chǎng)的耦合數(shù)據(jù),構(gòu)建了3種典型路徑的臺(tái)風(fēng)風(fēng)場(chǎng),并采用SWAN模型模擬了長(zhǎng)江口的波浪場(chǎng),研究結(jié)果表明該種方法能夠較好地復(fù)演長(zhǎng)江口水域臺(tái)風(fēng)期間的波況。

c. 研究表明長(zhǎng)江口東側(cè)海域過(guò)境臺(tái)風(fēng)更有可能造成北槽波能集中,為影響驟淤的典型天氣過(guò)程,數(shù)值計(jì)算結(jié)果與實(shí)測(cè)結(jié)果相吻合。

需要指出,本文針對(duì)的是汛期臺(tái)風(fēng)期間的驟淤與氣象過(guò)程的相關(guān)性分析。而在歷史上,寒潮大風(fēng)期間,長(zhǎng)江口深水航道也發(fā)生過(guò)偶然性驟淤現(xiàn)象,2010—2015年間共發(fā)生2次。寒潮大風(fēng)引起的驟淤和臺(tái)風(fēng)期間的驟淤?gòu)目臻g分布至誘發(fā)機(jī)制上均可能存在差異,本文沒(méi)有統(tǒng)計(jì)在內(nèi),今后將深入研究。

參考文獻(xiàn)：

[ 1 ] 趙德招,劉杰,吳華林. 近十年來(lái)臺(tái)風(fēng)誘發(fā)長(zhǎng)江口航道驟淤的初步分析[J]. 泥沙研究,2012,37(2):54-60. (ZHAO Dezhao, LIU Jie, WU Hualin. Preliminary analysis of typhoon-induced sudden sedimentation in navigation channel in Yangtze Estuary over last decade[J]. Journal of Sediment Research, 2012,37(2):54-60.(in Chinese))

[ 2 ] 楊世倫,丁平興,趙慶英. 開(kāi)敞大河口灘槽沖淤對(duì)臺(tái)風(fēng)的響應(yīng)及其動(dòng)力泥沙機(jī)制探討:以長(zhǎng)江口南匯邊灘—南槽—九段沙系統(tǒng)為例[J]. 海洋工程,2002, 20(3): 69-75.(YANG Shilun,DING Pingxing,ZHAO Qingying. Morphodynamic response of a large river mouth to typhoons[J]. The Ocean Engineering,2002, 20(3): 69-75.(in Chinese))

[ 3 ] 張志忠,阮文杰,蔣國(guó)俊. 長(zhǎng)江口動(dòng)水絮凝沉降與攔門(mén)沙淤積的關(guān)系[J]. 海洋與湖沼,1995,26(6): 632-638. (ZHANG Zhizhong,RUAN Wenjie,JIANG Guojun. The relationship between flocculating settlement in flowing water and deposition in the Changjiang River Estuary Mouth Bar Area[J]. Oceanologia et Limnologia Sinica,1995,26(6): 632-638.(in Chinese))

[ 4 ] 孔令雙,戚定滿,曹祖德,等. 長(zhǎng)江口航道淤積三層模式的建立與應(yīng)用[J]. 泥沙研究,2010,36(3): 1-8. (KONG Lingshuang,QI Dingman,CAO Zude,et al. Three-layer model for sedimentation in the deep navigation channel of the Yangtze Estuary[J]. Journal of Sediment Research,2010,36(3): 1-8.(in Chinese))

[ 5 ] 劉猛,吳華林,李為華,等. 長(zhǎng)江口北槽南導(dǎo)堤越堤輸沙運(yùn)動(dòng)的時(shí)空變化特征研究[J]. 海洋學(xué)報(bào)(中文版), 2013,35(2): 161-167. (LIU Meng,WU Hualin,LI Weihua, et al. The temporal and spatial variations of sediment flux overtopping the South Leading Jetty in the North Passage of the Changjiang Estuary[J]. Acta Oceanologica Sinica,2013,35(2): 161-167. (in Chinese))

[ 6 ] 侯志強(qiáng),楊華. 黃驊港外航道驟淤分析[J]. 水道港口,2004(4): 213-215. (HOU Zhiqiang,YANG Hua. Analysis of sudden siltation of outer channel of Huanghua Harbor[J]. Journal of Waterway and Harbor, 2004(4): 213-215.(in Chinese))

[ 7 ] 羅肇森. 大風(fēng)期黃驊港外航道的驟淤估算及防淤減淤措施探討[J]. 水運(yùn)工程,2004(10): 69-73. (LUO Zhaosen. Estimation of sudden sedimentation in outer channel of Huanghua Port during rough climate and measures for reduction of sedimentation[J]. Port & Waterway Engineering, 2004(10): 69-73.(in Chinese))

[ 8 ] 曹祖德. 日本熊本港的淤積分析[J]. 港口工程,1990(3): 28-36. (CAO Zude.Deposition analysis of Kumamoto Harbor in Japan[J].Harbour Engineering,1990(3): 28-36.(in Chinese))

[ 9 ] 羅肇森. 波、流共同作用下的近底泥沙輸移及航道驟淤預(yù)報(bào)[J]. 泥沙研究,2004,30(6): 1-9. (LUO Zhaosen. Sediment transport under the coexisting action of waves and currents and the prediction of sudden sedimentation in navigation channel[J]. Journal of Sediment Research,2004,30(6): 1-9.(in Chinese))

[10] 徐福敏,張長(zhǎng)寬. 臺(tái)風(fēng)浪對(duì)長(zhǎng)江口深水航道驟淤的影響研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展,2004,19(2): 137-143. (XU Fumin,ZHANG Changkuan. Study of the effect of storm waves on the rapid deposition of the Yangtze River Estuary channel[J]. Journal of Hydrodynamics, 2004,19(2): 137-143. (in Chinese))

[11] 劉猛,李為華. 長(zhǎng)江口北槽深水航道回淤量變化宏觀動(dòng)力原因分析[J]. 水運(yùn)工程,2013(3): 129-139. (LIU Meng, LI Weihua. Analysis of dynamic reasons for change of back silting amount in deepwater channel project at north passage of the Yangtze estuary[J]. Port & Waterway Engineering,2013(3): 129-139. (in Chinese))

[12] 白志剛,魏茂興. 基于時(shí)域模型的一種新的波浪能計(jì)算方法[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版),2014,47(6):491-497. (BAI Zhigang,WEI Maoxing. A new method for calculating wave power based on time domain model[J]. Journal of Tianjin University (Science and Technology), 2014,47(6):491-497.(in Chinese))

[13] 王喜年,尹慶江,張保明. 中國(guó)海臺(tái)風(fēng)風(fēng)暴潮預(yù)報(bào)模式的研究與應(yīng)用[J]. 水科學(xué)進(jìn)展, 1991,2(1): 1-10. (WANG Xinian, YIN Qingjiang,ZHANG Baoming.Research and application of a forecasting model of typhoon surges in China Seas[J]. Advances in Water Science,1991,2(1): 1-10.(in Chinese))

[14] ZHANG L, KUMAR A, WANG W. Influence of changes in observations on precipitation: a case study for the climate forecast system reanalysis (CFSR)[J]. Journal of Geophysical Research: Atmospheres, 2012, 117(D8):393-407.

[15] HODGES K I, LEE R W, BENGTSSON L. A comparison of extratropical cyclones in recent reanalyses ERA-Interim, NASA MERRA, NCEP CFSR, and JRA-25[J]. Journal of Climate,2011, 24(18): 4888-4906.

[16] 邱大洪. 工程水文學(xué)[M]. 北京: 人民交通出版社, 1999.

[17] 賈曉,潘軍寧,BARDUR Niclasen. SWAN模型風(fēng)能輸入項(xiàng)的改進(jìn)與驗(yàn)證[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010,38(5): 585-591. (JIA Xiao,PAN Junning,BARDUR Niclasen. Improvement and validation of wind energy input in SWAN model[J]. Journal of Hohai University (Natural Sciences), 2010,38(5): 585-591.(in Chinese))

[18] 李永博,吳克儉,畢凡,等. 基于SWAN模式的成山頭海域波浪能資源評(píng)估[J]. 海洋湖沼通報(bào), 2013(3): 1-9. (LI Yongbo,WU Kejian,BI Fan,et al. Wave energy resources assessment in the Chenshan Cape Sea during the last 20 years by using SWAN model[J]. Transaction of Oceanology and Limnology, 2013(3): 1-9.(in Chinese))

[19] 宗芳伊,魏來(lái),吳克儉. 基于SWAN模式的2種網(wǎng)格對(duì)嵊泗海域海浪模擬的比較研究[J]. 海洋湖沼通報(bào),2014(4): 10-18. (ZONG Fangyi,WEI Lai,WU Kejian. Research on wave simulation comparisons of Shensi Waters based on two kinds of grids using SWAN wave model[J]. Transaction of Oceanology and Limnology,2014(4): 10-18. (in Chinese))

[20] 胡克林,丁平興,朱首賢,等. 長(zhǎng)江口附近海域臺(tái)風(fēng)浪的數(shù)值模擬:以鹿沙臺(tái)風(fēng)和森拉克臺(tái)風(fēng)為例[J]. 海 洋 學(xué) 報(bào),2004,26(5):23-33.(HU Kelin,DING Pingxing,ZHU Shouxian,et al. Numerical simulation of typhoon around the waters of the Changjiang Estuary: a case study of Typhoon Rusa and Typhoon Sinlaku[J]. Acta Oceanologica Sinica,2004,26(5):23-33.(in Chinese))

TheextraordinaryweatherprocessinducingsuddensiltationindeepwaternavigationchannelofYangtzeEstuaryIthetyphoonpathcharacteristicsandnumericalverification

HUANGHuacong1,JIAXiao1,2,LUChuanteng3

(1.ShanghaiEstuarineandCoastalScienceResearchCenter，Shanghai201201，China；2.MOEKeyLaboratoryofHydrodynamics,SchoolofNAOCE,ShanghaiJiaoTongUniversity,Shanghai200240,China; 3.RiverandHarborEngineeringDepartment,NanjingHydraulicResearchInstitute,Nanjing210029,China)

To summarize the occurrence regularity and characteristics of sudden siltation during typhoon period in the deepwater navigation channel of Yangtze Estuary, the meteorological conditions and corresponding wave conditions during sudden siltation are analyzed. The results show that the wave energy correlates well with the siltation in NPJ station. The analysis of typhoon paths suggests that the sudden siltation in the navigation channel is significantly influenced by the taphoon passing through the east part of Yangtze Estuary. In this study, three typical typhoon paths are selected with reference to the historic typhoon paths in Yangtze Estuary, and the typhoon field was replicated by combing the empirical circular(Fujita and Takahashi)wind field and CFSR(climate forecast system reanalysis)wind field. SWAN model is then used to simulate the wave conditions under different typhoons paths. Finally, the effects of different typhoon paths on the sudden siltation are analyzed by taking the shallow water wave energy in NPJ station as the discriminant parameter. The study indicates that the typhoon passing through the east area of the Yangtze Estuary is found to be a typical typhoon path, which can readily induce relatively large wave energy and subsequent sudden siltation. And such analysis result agrees well with that of measured typhoon path causing sudden siltation since 2010.

deep water channel of Yangtze Estuary; sudden siltation; typhoon path; Fujita-Takahashi wind field; CFSR wind field; SWAN wave model; wave energy

10.3876/j.issn.1000-1980.2017.05.009

2016-07-26

國(guó)家自然科學(xué)基金青年基金(41206072);2016年港口航道泥沙交通行業(yè)重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(Yn216001)

黃華聰(1994—),男,江西贛州人,實(shí)習(xí)研究員,主要從事港口航道、海岸與近海工程研究。E-mail:huanghc1994@126.com

TV148;P732

A

1000-1980(2017)05-0432-07