唐 磊，孫林云，孫 波，肖立敏

(南京水利科學(xué)研究院 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210029)

海河流域主要河口閘下通道風(fēng)暴潮增水研究

唐 磊，孫林云，孫 波，肖立敏

(南京水利科學(xué)研究院 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210029)

建立東海-黃海-渤海大范圍風(fēng)暴潮二維數(shù)學(xué)模型,對(duì)9216臺(tái)風(fēng)風(fēng)暴潮過(guò)程進(jìn)行模擬，分析渤海灣造陸工程前后風(fēng)暴潮增水對(duì)海河流域主要河口(海河口、永定新河口和獨(dú)流減河口)閘下區(qū)域最高潮位的影響，并探討有效降低河口閘下通道內(nèi)風(fēng)暴潮潮位的開(kāi)發(fā)方式。研究表明：渤海灣圍填海工程實(shí)施后，閘下通道內(nèi)維持自然地形時(shí)，通道內(nèi)最高潮位高于工程前；閘下通道實(shí)施地形開(kāi)挖后，通道內(nèi)最高潮位依然高于工程前，但較自然地形條件下抬高幅值有所減小。閘下通道結(jié)合開(kāi)發(fā)建設(shè)港池航道等工程實(shí)施地形開(kāi)挖，使得通道內(nèi)水深增加，可降低通道內(nèi)風(fēng)暴潮增水幅度。

渤海灣；造陸工程；海河流域；主要河口；閘下通道；風(fēng)暴潮增水

Abstract: A wide-range 2D storm surge numerical model covering the East China Sea, the Yellow Sea and the Bohai Sea was established. The process of the storm surge of Tyhoon 9216 was the simulated in the model. The impact of the reclamation project on the highest tidal level (HTL) in the downstream channel of the tidal gate at the major estuaries of the Haihe River Basin during the storm surge was analyzed. And the development method to effectively reduce the storm tidal level in the downstream channel of the tidal gate was discussed. Results show that when maintaining the original terrain in the downstream channel of the tidal gate, the HTL in the channel is higher than that before the reclamation project; when excavation of the terrain in the channel is implemented, the HTL is still higher than before, but the amplitude of the storm surge in the channel will be less than that of the original terrain in the downstream channel. So the excavation of terrain to increase water depth by the development project of the harbor basin and waterway can weaken the storm surge in the channel.

Keywords: Bohai Bay; reclamation project; Haihe River Basin; major estuaries; downstream channel of tidal gate; storm surge

渤海灣是渤海三大海灣之一，是一個(gè)向西凹入的弧形淺水海灣，灣口以河北唐山市樂(lè)亭縣大清河口至山東的舊黃河口一線為界，面積約1.75×104km2[1]。渤海灣水深一般小于20 m，灣內(nèi)地形自灣頂向渤海中央緩傾，灘面坡度平緩，沉積物以細(xì)顆粒黏性土為主，為淤泥質(zhì)海岸，局部區(qū)域沉積物為細(xì)粉沙。渤海灣是我國(guó)風(fēng)暴潮災(zāi)害最嚴(yán)重的地區(qū)之一，一年四季均有發(fā)生，夏季有臺(tái)風(fēng)風(fēng)暴潮災(zāi)害發(fā)生，春、秋、冬季均有災(zāi)害性溫帶風(fēng)暴潮發(fā)生[2]。據(jù)天津市水利志記載，自1890年至1985年間，天津沿海共發(fā)生有較大影響的風(fēng)暴潮20次，1985年以來(lái)共發(fā)生了約16次由臺(tái)風(fēng)或溫帶氣旋影響在渤海灣產(chǎn)生的較嚴(yán)重風(fēng)暴潮災(zāi)害，其中典型風(fēng)暴潮有7203、9216、9711與1210臺(tái)風(fēng)風(fēng)暴潮和2003年10月寒潮風(fēng)暴潮，9216臺(tái)風(fēng)風(fēng)暴潮對(duì)該地區(qū)影響最為顯著，給社會(huì)經(jīng)濟(jì)和人民生命財(cái)產(chǎn)帶來(lái)巨大損失。

針對(duì)渤海灣風(fēng)暴潮研究是工程界和學(xué)術(shù)界長(zhǎng)期以來(lái)關(guān)注的熱點(diǎn)話題。文獻(xiàn)[2]采用塘沽海洋站實(shí)測(cè)潮位資料對(duì)天津沿海風(fēng)暴潮災(zāi)害發(fā)生時(shí)間、天數(shù)等進(jìn)行統(tǒng)計(jì)分析，并采用耿貝爾(Gumbel)方法計(jì)算出了不同重現(xiàn)期的風(fēng)暴潮增水值和高潮位值，為防潮工程設(shè)計(jì)提供了重要參數(shù)依據(jù)。徐明娥等[3]對(duì)2008年8月渤海灣內(nèi)發(fā)生的一次風(fēng)暴潮成因進(jìn)行研究。大量學(xué)者[4-9]采用風(fēng)暴潮數(shù)學(xué)模型對(duì)作用于渤?；虿澈硟?nèi)風(fēng)暴潮過(guò)程進(jìn)行模擬和預(yù)報(bào)，此類研究并未考慮海岸線變化對(duì)風(fēng)暴潮的影響。近十年來(lái)，為適應(yīng)經(jīng)濟(jì)發(fā)展需要，渤海灣內(nèi)開(kāi)展了大量的圍填海造陸工程，如曹妃甸、濱海旅游區(qū)、天津港、臨港經(jīng)濟(jì)區(qū)、南港工業(yè)區(qū)、黃驊港等，海岸線發(fā)生了較大變化，章衛(wèi)勝等[10]建立了渤海風(fēng)暴潮二維數(shù)學(xué)模型，研究分析了圍填海工程實(shí)施對(duì)渤海灣內(nèi)及工程前沿風(fēng)暴潮高潮位的影響情況。海河流域主要入海河口有海河口、永定新河口和獨(dú)流減河口(簡(jiǎn)稱“三河口”)，隨著天津?yàn)I海新區(qū)圍填海工程的實(shí)施，“三河口”新的河口延伸至海圖水深3～9 m，與防潮閘之間形成約15～20 km不等的閘下通道，風(fēng)暴潮作用下河口閘下通道潮位變化對(duì)河口行洪安全的影響是防洪的關(guān)鍵問(wèn)題。

建立東海-黃海-渤海大范圍風(fēng)暴潮數(shù)學(xué)模型，對(duì)9216臺(tái)風(fēng)風(fēng)暴潮過(guò)程進(jìn)行模擬與驗(yàn)證，計(jì)算分析圍填海造陸工程對(duì)海河流域主要河口閘下通道內(nèi)風(fēng)暴潮增水的影響，探討有效降低河口(三河口)閘下通道內(nèi)風(fēng)暴潮增水的開(kāi)發(fā)方式。

1 風(fēng)暴潮數(shù)學(xué)模型建立與驗(yàn)證

1.1 控制方程

采用靜壓假定下的不可壓縮淺水流動(dòng)方程描述風(fēng)暴潮運(yùn)動(dòng)，即納維爾—斯托克斯(Navier-Stokes)方程。本研究主要針對(duì)平面尺度較大的海域潮流計(jì)算，故采用垂線平均后的二維水流基本方程，球坐標(biāo)系下的方程表達(dá)形式如下：

連續(xù)方程：

運(yùn)動(dòng)方程：

式中：x、y為笛卡爾坐標(biāo)系坐標(biāo)；λ、φ為球坐標(biāo)系坐標(biāo)，分別為經(jīng)度和緯度；R為地球半徑，m；t為時(shí)間變量，s；η為相對(duì)于參考基面的水位，m；h=h0+η為全水深，m；U、V分別為x、y方向上的垂線平均流速，m/s；f為科氏力系數(shù)(f=2ωsinφ，ω為地球自轉(zhuǎn)角速度)；ρ0為水體參考密度，kg/m3；g為重力加速度，m/s2；pa為大氣氣壓，hPa；A為水平紊動(dòng)黏性系數(shù)，由Smagorinsky提出的亞格子法進(jìn)行計(jì)算[11]；τsx、τsy分別為表面風(fēng)應(yīng)力在x、y方向上的分量，N/m2；τbx、τby分別為底部切應(yīng)力在x、y方向上的分量，N/m2。

其中，ρa(bǔ)為空氣密度，kg/m3；W10為海面上10 m處風(fēng)速，m/s；wx、wy分別為W10在x、y方向的分量；cd為風(fēng)拖曳力系數(shù)，其經(jīng)驗(yàn)公式形式[12]如下：

1.2 模型建立

本次建立的風(fēng)暴潮數(shù)學(xué)模型采用球坐標(biāo)系，模型范圍117.5°E～127°E、31.7°N～41.0°N，包含東海、黃海和渤海，開(kāi)邊界設(shè)置在長(zhǎng)江北支口至韓國(guó)濟(jì)州島和濟(jì)州島至韓國(guó)本土岸線。模型網(wǎng)格示意如圖1(a)，網(wǎng)格節(jié)點(diǎn)數(shù)為26 024，網(wǎng)格單元數(shù)為48 725，最大網(wǎng)格尺度為0.19°，位于模型開(kāi)邊界處，最小網(wǎng)格尺度為0.002°，位于海河流域主要河口閘下通道內(nèi)。大范圍水下地形來(lái)自美國(guó)國(guó)家海洋和大氣管理局(NOAA)網(wǎng)站的公開(kāi)數(shù)據(jù)，渤海灣內(nèi)水下地形采用收集到的海圖和局部實(shí)測(cè)地形資料。將水下地形數(shù)據(jù)插值到模型網(wǎng)格節(jié)點(diǎn)上，得到模型計(jì)算地形，如圖1(b)所示。

1.3 參數(shù)選取

風(fēng)場(chǎng)和氣壓場(chǎng)是風(fēng)暴潮數(shù)學(xué)模型中的關(guān)鍵輸入?yún)?shù)，采用的風(fēng)場(chǎng)CCMP(Cross-Calibrated, Multi-Platform)資料來(lái)自ESE(NASA Earth Science Enterprise)，它結(jié)合了ADEOS-II、QuikSCAT、TRMM TMI、AMSR-E、SSM/I幾種資料，利用變分方法得到。CCMP風(fēng)場(chǎng)具有較高的精度和時(shí)空分辨率，其空間分辨率為0.25°×0.25°，時(shí)間分辨率為6 h?？臻g范圍為78.375°S～78.375°N，0.125°E～359.875°E，時(shí)間范圍為1987年7月-2011年12月。張鵬等[13]采用渤海海區(qū)某觀測(cè)站的實(shí)測(cè)風(fēng)速風(fēng)向數(shù)據(jù)對(duì)CCMP風(fēng)場(chǎng)進(jìn)行檢驗(yàn)，結(jié)果表明，CCMP風(fēng)場(chǎng)和實(shí)測(cè)風(fēng)場(chǎng)風(fēng)速風(fēng)向相關(guān)性較好，能夠代表渤海海面的風(fēng)場(chǎng)狀況。本次選取20.125°N～43.125°N、115.125°E～135.125°E范圍的風(fēng)場(chǎng)數(shù)據(jù)。氣壓場(chǎng)根據(jù)臺(tái)風(fēng)路徑、中心氣壓等參數(shù)采用經(jīng)驗(yàn)公式計(jì)算得到，具體計(jì)算方法參考文獻(xiàn)[14-15]，典型時(shí)刻風(fēng)場(chǎng)如圖2所示。

圖1 風(fēng)暴潮數(shù)學(xué)模型Fig. 1 Sketch of storm surge numerical model

圖2 9216臺(tái)風(fēng)路徑與典型時(shí)刻風(fēng)場(chǎng)和氣壓場(chǎng)Fig. 2 The wind and air pressure of Typhoon 9216 at typical moment and its path

1.4 模型驗(yàn)證

首先，采用國(guó)家海洋信息中心發(fā)布的《潮汐表》中渤海灣內(nèi)各潮位站預(yù)報(bào)值，對(duì)天文潮過(guò)程進(jìn)行驗(yàn)證，為節(jié)省篇幅，圖3(a)僅給出塘沽海洋站天文潮驗(yàn)證結(jié)果；其次，采用收集到的9216臺(tái)風(fēng)風(fēng)暴潮期間塘沽海洋站潮位實(shí)測(cè)過(guò)程，對(duì)風(fēng)暴潮期間該站位的潮位和增減水過(guò)程進(jìn)行驗(yàn)證，分別如圖3(a)和圖3(b)所示。從驗(yàn)證結(jié)果來(lái)看，模型計(jì)算的潮位和增水幅度與實(shí)測(cè)結(jié)果均較為接近，最高潮位計(jì)算值(3.16 m)與實(shí)測(cè)值(3.20 m)誤差為-0.04 m，表明該風(fēng)暴潮數(shù)學(xué)模型對(duì)現(xiàn)場(chǎng)風(fēng)暴潮過(guò)程的模擬具有一定的精度，可運(yùn)用于后續(xù)相關(guān)研究。

圖4 渤海灣造陸工程前9216臺(tái)風(fēng)風(fēng)暴潮最高潮位分布Fig. 4 The distribution of highest tide level in Bohai Bay before the reclamation during Typhoon 9216 storm surge

2 渤海灣造陸工程前后風(fēng)暴潮最高潮位分析

2.1 工程前

渤海灣造陸工程實(shí)施前，風(fēng)暴潮最高潮位受淺水變形作用從渤海灣口向內(nèi)逐漸抬升，渤海灣灣頂沿線風(fēng)暴潮最高潮位分布規(guī)律表現(xiàn)為由北向南逐漸增大，在南港工業(yè)區(qū)至濱州港(即套爾河口)一帶風(fēng)暴潮最高潮位最高，如圖4所示。“三河口”閘下區(qū)域風(fēng)暴潮最高潮位沿程分布規(guī)律是從外海至閘下逐漸抬升，如圖5所示?！叭涌凇敝歇?dú)流減河口閘下潮位最高、永定新河口次之、海河口較永定新河口略低，分別為3.34 m、3.20 m和3.18 m。

圖5 閘下通道最高潮位沿程分布Fig. 5 The distribution of highest tide level in downstream channel of tidal gate

圖6 渤海灣造陸工程后三河口通道內(nèi)計(jì)算網(wǎng)格 Fig. 6 The numeric grid in the three estuaries channel under tide gate after the reclamation in Bohai Bay

2.2 工程后

渤海灣造陸工程實(shí)施后，“三河口”防潮閘閘下均形成較長(zhǎng)通道，沿岸陸地邊界較工程前發(fā)生變化，在前文建立的風(fēng)暴潮數(shù)學(xué)模型基礎(chǔ)上，對(duì)河口區(qū)域的計(jì)算網(wǎng)格進(jìn)行加密重構(gòu)，如圖6。計(jì)算分析河口閘下通道內(nèi)自然和開(kāi)挖地形兩種不同工況下，通道內(nèi)潮位沿程分布情況；并計(jì)算分析“三河口”防潮閘同時(shí)提閘泄洪條件下，通道內(nèi)潮位的變化情況。各計(jì)算工況閘下通道潮位沿程分布和渤海灣內(nèi)分布情況分別如圖5和圖7，閘下最高潮位值和通道最高潮位平均增加值分別如表1和表2所示。

表1 渤海灣造陸工程后三河口閘下最高潮位Tab. 1 The highest tide levels in three channels under tidal gate after the reclamation in Bohai Bay

表2渤海灣造陸工程后三河口通道最高潮位平均增加

Tab.2AverageincrementofhighesttidelevelsinthreeestuarieschannelsundertidalgatesafterthereclamationinBohaiBay

臺(tái)風(fēng)比較項(xiàng)目最高潮位平均增加值/m海河口永定新河口獨(dú)流減河口9216臺(tái)風(fēng)工程后較工程前(自然地形工況)0.120.110.16工程后較工程前(自然地形與泄洪工況)0.140.110.20工程后較工程前(開(kāi)挖地形工況)0.060.060.10工程后較工程前(開(kāi)挖地形與泄洪工況)0.070.090.11

2.2.1 閘下通道自然地形工況

“三河口”在防潮閘閘下均配套有清淤槽工程，造陸工程實(shí)施后，形成的閘下通道內(nèi)維持自然地形(即閘下清淤槽開(kāi)挖而通道內(nèi)不開(kāi)挖)，通道內(nèi)最高潮位因潮波淺水變形效應(yīng)增強(qiáng)而較工程前有所抬高。在9216臺(tái)風(fēng)風(fēng)暴潮作用下，海河口、永定新河口和獨(dú)流減河口閘下通道最高潮位較工程前分別平均抬高0.12 m、0.11 m和0.16 m。獨(dú)流減河口閘下潮位最高、永定新河口次之，海河口較永定新河口略低。

閘下通道內(nèi)維持自然地形，海河口、永定新河口和獨(dú)流減河口三河口同時(shí)提閘泄洪，流量分別800 m3/s、4 640 m3/s和3 600 m3/s。閘下通道內(nèi)最高潮位不僅較工程前有所抬高，下泄洪水在風(fēng)暴潮頂托作用下，使得該工況通道內(nèi)潮位較自然地形和不泄洪情況下進(jìn)一步抬高。

2.2.2 閘下通道開(kāi)挖地形工況

造陸工程實(shí)施后，結(jié)合閘下清淤和港池、航道與碼頭前沿開(kāi)挖，閘下通道內(nèi)實(shí)施地形開(kāi)挖后，通道內(nèi)最高潮位較工程前仍有所抬高，但較河口閘下通道自然地形情況下的抬高幅度有所減小，獨(dú)流減河口閘下潮位仍為三河口最高。在9216臺(tái)風(fēng)風(fēng)暴潮作用下，海河口、永定新河口和獨(dú)流減河口閘下最高潮位分別平均抬高0.06 m、0.06 m和0.10 m。

閘下通道內(nèi)開(kāi)挖地形，三河口同時(shí)泄洪，通道內(nèi)最高潮位較工程前有所抬高，與閘下通道自然地形泄洪工況比較，其潮位抬高幅值有所減小。

綜上所述，渤海灣造陸工程實(shí)施后，閘下通道內(nèi)維持自然地形時(shí)，最高潮位高于工程前；閘下通道結(jié)合開(kāi)發(fā)建設(shè)港區(qū)航道等工程實(shí)施地形開(kāi)挖后，最高潮位依然高于工程前，但較自然地形下抬高值有所減小?？梢?jiàn)，閘下通道結(jié)合開(kāi)發(fā)建設(shè)港池航道等工程實(shí)施地形開(kāi)挖，使得通道內(nèi)水深增加，可降低通道內(nèi)風(fēng)暴潮增水幅度。

圖7 渤海灣造陸工程后9216臺(tái)風(fēng)風(fēng)暴潮最高潮位分布Fig. 7 The distribution of highest tide level in Bohai Bay after the reclamation during Typhoon 9216 storm surge

3 結(jié) 語(yǔ)

建立了東海-黃海-渤海大范圍風(fēng)暴潮數(shù)學(xué)模型，采用了CCMP風(fēng)場(chǎng)和經(jīng)驗(yàn)公式計(jì)算所得氣壓場(chǎng)作為模型外部驅(qū)動(dòng)力的輸入?yún)?shù)，對(duì)9216臺(tái)風(fēng)風(fēng)暴潮過(guò)程進(jìn)行模擬，并采用潮位和增減水實(shí)測(cè)資料對(duì)模型進(jìn)行驗(yàn)證。在此基礎(chǔ)上，分析了工程前渤海灣內(nèi)最高潮位分布情況，并計(jì)算分析了工程后海河流域主要河口閘下通道內(nèi)風(fēng)暴潮潮位變化情況，探討了降低通道內(nèi)風(fēng)暴潮增水的開(kāi)發(fā)方式。得到如下主要結(jié)論：

1)渤海灣造陸工程前，風(fēng)暴潮最高潮位受淺水變形作用從渤海灣口向內(nèi)逐漸抬升，渤海灣灣頂沿線風(fēng)暴潮最高潮位分布規(guī)律表現(xiàn)為由北向南逐漸增大，在獨(dú)流減河口附近至濱州港(即套爾河口)一帶潮位最高?！叭涌凇遍l下區(qū)域風(fēng)暴潮最高潮位沿程分布規(guī)律是從外海至閘下逐漸抬升，其中獨(dú)流減河口閘下潮位最高、永定新河口次之、海河口較永定新河口略低。

2)渤海灣造陸工程實(shí)施后，三河口閘下均形成較長(zhǎng)的通道，若通道內(nèi)維持自然地形，其內(nèi)最高潮位因潮波淺水變形較工程前有所抬高。三河口閘下最高潮位表現(xiàn)為獨(dú)流減河口最高、永定新河口次之，海河口較永定新河口略低。在此條件下提閘泄洪，閘下通道內(nèi)潮位不僅較工程前有所抬高，并且較河口閘下通道自然地形和不泄洪情況下會(huì)進(jìn)一步抬高。

3)在三河口形成的通道內(nèi)進(jìn)行地形開(kāi)挖，通道內(nèi)最高潮位較工程前仍有所抬高，但較閘下通道自然地形情況下的抬高幅度有所減小。在此條件下提閘泄洪，通道內(nèi)潮位較不泄洪情況下有所抬高，但較自然地形條件下泄洪時(shí)的潮位有所降低。

渤海灣造陸工程實(shí)施后，“三河口”閘下較長(zhǎng)的通道內(nèi)風(fēng)暴潮增水有所增加，對(duì)河口行洪安全產(chǎn)生一定程度的影響；如閘下通道內(nèi)開(kāi)發(fā)建設(shè)港池航道等工程，實(shí)施地形開(kāi)挖增加通道內(nèi)水深，可降低風(fēng)暴潮增水幅度，進(jìn)而減緩其對(duì)河口行洪安全的影響。本研究成果可為其他類似河口灘涂開(kāi)發(fā)建設(shè)提供借鑒。

[1] 孫湘平. 中國(guó)近海區(qū)域海洋[M]. 北京: 海洋出版社, 2006:1-2. (SUN Xiangpin. China offshore area [M]. Beijing: Ocean Press, 2006:1-2. (in Chinese))

[2] 吳少華, 王喜年, 宋珊, 等. 天津沿海風(fēng)暴潮災(zāi)害概述及統(tǒng)計(jì)分析[J]. 海洋預(yù)報(bào), 2002, 19(1):29-35. (WU Shaohua, WANG Xinnian, SONG Shan, et al. The general status and statistics of storm surges in the coast area of Tianjin [J]. Marine Forecasts, 2002, 19(1):29-35. (in Chinese))

[3] 徐明娥, 趙玉潔, 劉建軍, 等. 濱海新區(qū)一次氣旋型風(fēng)暴潮的成因分析[J]. 安徽農(nóng)業(yè)科學(xué), 2013, 41(5):2141-2142, 2171. (XU Minge, ZHAO Yujie, JIU Jianjun, et al. Analysis of a storm surge in Binhai New area [J]. Journal of Anhui Agri. Sci., 2013, 41(5):2141-2142, 2171. (in Chinese))

[4] 李大鳴, 徐亞男, 白玲, 等. 渤海灣溫帶風(fēng)暴潮數(shù)值預(yù)報(bào)模型[J]. 天津大學(xué)學(xué)報(bào), 2011, 44(9):840-846. (LI Daming, XU Yanan, BAI Ling, et al. Numerical forecast model of variable zone storm surge for Bohai Bay [J]. Journal of Tianjin University, 2011, 44(9):840-846. (in Chinese))

[5] 李鑫, 張金善, 章衛(wèi)勝. 風(fēng)暴潮耦合數(shù)值模式在渤海海域中的應(yīng)用[J]. 水運(yùn)工程, 2009(10):25-31. (LI Xin, ZHANG Jinshan, ZHANG Weisheng. Application of storm surge numerical coupling model on Bohai Sea [J]. Port & Waterway Engineering, 2009(10):25-31. (in Chinese))

[6] 李艷蕓, 李紹武. 風(fēng)暴潮預(yù)報(bào)模式在渤海海域中的應(yīng)用研究[J]. 海洋技術(shù), 2006, 25(1):101-106. (LI Yanyun, LI Shaowu. Application research of a storm surge prediction model in Bohai Sea [J]. Ocean Technology, 2006, 25(1):101-106. (in Chinese))

[7] 夏波, 吳相忠, 張慶河. 渤海灣西南海域潮汐風(fēng)暴潮數(shù)值模擬[J]. 長(zhǎng)沙理工大學(xué)學(xué)報(bào)：自然科學(xué)版, 2008, 5(3):64-70. (XIA Bo, WU Xiangzhong, ZHANG Qinghe. Modeling of tide and storm surge in the southwest coast of the Bohai bay [J]. Journal of Changsha University of Science and Technology, Natural Science, 2008, 5(3):64-70. (in Chinese))

[8] 吳少華, 王喜年, 戴明瑞, 等. 渤海風(fēng)暴潮概況及溫帶風(fēng)暴潮數(shù)值模擬[J]. 海洋學(xué)報(bào), 2002, 24(3):28-34. (WU Shaohua, WANG Xinian, DAI Mingrui, et al. The general status of storm surges and the simulation of extra-tropical storm surges in the Bohai Sea [J]. ACTA Oceanologica Sinica, 2002, 24(3):28-34. (in Chinese))

[9] 李大鳴, 徐亞男, 宋雙霞, 等. 波浪輻射應(yīng)力在渤海灣海域?qū)︼L(fēng)暴潮影響的研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2010, 25(3):374-382. (LI Daming, XU Yanan, SONG Shuangxia, et al. The study on the effect of wave-radiation stress on storm surge in Bohai Bay [J]. Chinese Journal of Hydrodynamics, 2010, 25(3):374-382. (in Chinese))

[10] 章衛(wèi)勝, 陳晨, 李鑫, 等. 海岸工程對(duì)渤海灣風(fēng)暴潮高潮位影響分析[J]. 海洋工程, 2012, 30(2):72-78. (ZHANG Weisheng, CHEN Chen, LI Xin, et al. The impact of coastal engineering on high water level of storm surges in Bohai Bay [J]. The Ocean Engineering, 2012, 30(2):72-78. (in Chinese))

[11] SMAGORINSKY J. General circulation experiment with the primitive equations [J]. Monthly Weather Review, 1963, 91(3):99-164.

[12] WU Jin. Wind-stress coefficients over sea surface and near neutral conditions—A revisit [J]. Journal of Physical Oceanography, 1980, 10(5):727-740.

[13] 張鵬, 陳曉玲, 陸建忠, 等. 基于CCMP衛(wèi)星遙感海面風(fēng)場(chǎng)數(shù)據(jù)的渤海風(fēng)浪模擬研究[J]. 海洋通報(bào), 2011, 30(3):266-271. (ZHANG Peng, CHEN Xiaoling, LU Jianzhong, et al. Research on wave simulation of Bohai Sea based on the CCMP remotely sensed sea winds [J]. Marine Science Bulletin, 2011, 30(3):266-271. (in Chinese))

[14] THOMPSON E F, CARDONE V J. Practical modeling of hurricane surface wind fields [J]. Journal of Waterway, Coastal, and Ocean Engineering, 1996, 122(4):195-205.

[15] 肖玉鳳, 肖儀清, 段忠東, 等. CE風(fēng)場(chǎng)模型在中國(guó)華南沿海地區(qū)的適用性研究[J]. 工程力學(xué), 2010, 27(10):251-256. (XIAO Yufeng, XIAO Yiqing, DUAN Zhongdong, et al. The applicability of CE wind-field model in south China coastal region [J]. Engineering Mechanics, 2010, 27(10):251-256. (in Chinese))

Study on storm surge in downstream channels of tidal gates in the major estuaries of Haihe River Basin

TANG Lei, SUN Linyun, SUN Bo, XIAO Limin

(Nanjing Hydraulic Research Institute, State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Nanjing 210029, China)

1005-9865(2016)04-0055-07

P731.23

A

10.16483/j.issn.1005-9865.2016.04.008

2015-12-09

水利部公益性行業(yè)科研專項(xiàng)(201301067)；南京水利科學(xué)研究院中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(Y215018)

唐 磊(1985-)，男，江蘇靖江人，工程師，博士，主要從事港口航道工程及工程泥沙研究。E-mail：tlhppy@163.com

孫林云(1959-)，男，上海人，教授級(jí)高級(jí)工程師，主要從事河口海岸工程泥沙問(wèn)題及河口綜合利用規(guī)劃研究。 E-mail：lysun@nhri.cn