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

(南京水利科學研究院 水文水資源與水利工程科學國家重點實驗室, 江蘇 南京 210029)

海河流域主要河口閘下通道風暴潮增水研究

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

(南京水利科學研究院 水文水資源與水利工程科學國家重點實驗室, 江蘇 南京 210029)

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

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

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

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

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

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

1 風暴潮數(shù)學模型建立與驗證

1.1 控制方程

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

連續(xù)方程：

運動方程：

式中：x、y為笛卡爾坐標系坐標；λ、φ為球坐標系坐標，分別為經(jīng)度和緯度；R為地球半徑，m；t為時間變量，s；η為相對于參考基面的水位，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為水平紊動黏性系數(shù)，由Smagorinsky提出的亞格子法進行計算[11]；τsx、τsy分別為表面風應(yīng)力在x、y方向上的分量，N/m2；τbx、τby分別為底部切應(yīng)力在x、y方向上的分量，N/m2。

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

1.2 模型建立

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

1.3 參數(shù)選取

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

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

圖2 9216臺風路徑與典型時刻風場和氣壓場Fig. 2 The wind and air pressure of Typhoon 9216 at typical moment and its path

1.4 模型驗證

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

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

2 渤海灣造陸工程前后風暴潮最高潮位分析

2.1 工程前

渤海灣造陸工程實施前，風暴潮最高潮位受淺水變形作用從渤海灣口向內(nèi)逐漸抬升，渤海灣灣頂沿線風暴潮最高潮位分布規(guī)律表現(xiàn)為由北向南逐漸增大，在南港工業(yè)區(qū)至濱州港(即套爾河口)一帶風暴潮最高潮位最高，如圖4所示?！叭涌凇遍l下區(qū)域風暴潮最高潮位沿程分布規(guī)律是從外海至閘下逐漸抬升，如圖5所示。“三河口”中獨流減河口閘下潮位最高、永定新河口次之、海河口較永定新河口略低，分別為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)計算網(wǎng)格 Fig. 6 The numeric grid in the three estuaries channel under tide gate after the reclamation in Bohai Bay

2.2 工程后

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

臺風比較項目最高潮位平均增加值/m海河口永定新河口獨流減河口9216臺風工程后較工程前(自然地形工況)0.120.110.16工程后較工程前(自然地形與泄洪工況)0.140.110.20工程后較工程前(開挖地形工況)0.060.060.10工程后較工程前(開挖地形與泄洪工況)0.070.090.11

2.2.1 閘下通道自然地形工況

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

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

2.2.2 閘下通道開挖地形工況

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

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

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

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

3 結(jié) 語

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

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

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

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

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

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

[2] 吳少華, 王喜年, 宋珊, 等. 天津沿海風暴潮災(zāi)害概述及統(tǒng)計分析[J]. 海洋預(yù)報, 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ū)一次氣旋型風暴潮的成因分析[J]. 安徽農(nóng)業(yè)科學, 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] 李大鳴, 徐亞男, 白玲, 等. 渤海灣溫帶風暴潮數(shù)值預(yù)報模型[J]. 天津大學學報, 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)勝. 風暴潮耦合數(shù)值模式在渤海海域中的應(yīng)用[J]. 水運工程, 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] 李艷蕓, 李紹武. 風暴潮預(yù)報模式在渤海海域中的應(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] 夏波, 吳相忠, 張慶河. 渤海灣西南海域潮汐風暴潮數(shù)值模擬[J]. 長沙理工大學學報：自然科學版, 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] 吳少華, 王喜年, 戴明瑞, 等. 渤海風暴潮概況及溫帶風暴潮數(shù)值模擬[J]. 海洋學報, 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暴潮影響的研究[J]. 水動力學研究與進展, 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)勝, 陳晨, 李鑫, 等. 海岸工程對渤海灣風暴潮高潮位影響分析[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)星遙感海面風場數(shù)據(jù)的渤海風浪模擬研究[J]. 海洋通報, 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風場模型在中國華南沿海地區(qū)的適用性研究[J]. 工程力學, 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è)科研專項(201301067)；南京水利科學研究院中央級公益性科研院所基本科研業(yè)務(wù)費專項(Y215018)

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

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