天津海域圍填海工程對渤海灣水交換的影響研究

靳玉丹，張秋豐，李希彬，王魯寧，葉風(fēng)娟

（1.上海海洋大學(xué)，上海201306；2.國家海洋局天津海洋環(huán)境監(jiān)測中心站，天津300457）

天津海域圍填海工程對渤海灣水交換的影響研究

靳玉丹1，2，張秋豐1，2，李希彬2，王魯寧2，葉風(fēng)娟2

（1.上海海洋大學(xué)，上海201306；2.國家海洋局天津海洋環(huán)境監(jiān)測中心站，天津300457）

利用三維海洋數(shù)值模型FVCOM，進(jìn)行渤海灣三維水動力和水交換數(shù)值模擬，經(jīng)實(shí)測潮汐和潮流資料驗(yàn)證，模型模擬結(jié)果較好。然后采用該模型對渤海灣內(nèi)的水體水交換能力進(jìn)行定量研究。研究結(jié)果表明，在天津海域進(jìn)行圍填海工程之前渤海灣水體的半交換周期為300 d左右，圍填海之后，水體半交換周期延長25 d，渤海灣西部水體的水交換率下降可達(dá)10%，半交換周期延長92 d。尤其是天津沿海南部海域的水交換能力下降嚴(yán)重，圍填海之后其水體半交換周期延長可達(dá)200 d。渤海灣北部也有部分海域水交換周期延長達(dá)200 d。建議在進(jìn)行圍填海工程建設(shè)時(shí)，應(yīng)將工程對水體交換能力的影響納入考慮，避免因圍填海工程因素造成的惡劣環(huán)境影響。

渤海灣；圍填海工程；水交換影響

渤海灣是渤海西部的一個(gè)海灣，位于河北唐山、天津、河北滄州和山東黃河口之間，有海河、永定新河、獨(dú)流減河等河流注入，是京津的海上門戶，華北的海運(yùn)樞紐。渤海灣水深較淺，灘涂面積大，生態(tài)環(huán)境較為脆弱。特別是近年來經(jīng)濟(jì)發(fā)展迅速，圍填海工程大規(guī)模進(jìn)行，給渤海灣海域的水環(huán)境造成了很大的壓力，尤其是對該海域水交換的影響，已經(jīng)引起諸多學(xué)者的廣泛關(guān)注（李靜等，2015）。

在海灣的水交換過程研究方面已經(jīng)有很多學(xué)者涉足，并且進(jìn)行了大量的研究，如有些學(xué)者基于平均存留時(shí)間的概念研究了膠州灣的水交換能力，計(jì)算結(jié)果正確的反映了膠州灣的水交換情況（劉哲，2004；曹振東，2011）。還有學(xué)者建立了欽州灣海域保守物質(zhì)的對流-擴(kuò)散模型，基于半交換時(shí)間概念對該海域不同時(shí)期的水交換能力進(jìn)行了數(shù)值模擬，并進(jìn)行分析（張坤，2014）。另外有利用delft3D軟件中的How模塊，以保守物質(zhì)作為示蹤劑，建立對流-擴(kuò)散型的海灣水交換數(shù)值模型來模擬象山港水交換過程，研究象山港整體與外海之間水體交換特性（彭輝，2013）。對于渤海灣的水交換過程，前人也有研究，如采用二維水動力模型對渤海灣的水交換情況進(jìn)行了模擬，認(rèn)為當(dāng)時(shí)情況條件下對渤海灣水交換影響最大的因素是風(fēng)（何磊，2004）。還有人提出了關(guān)聯(lián)矩陣的方法，并利用該方法對渤海灣水交換進(jìn)行初步分區(qū)模擬，通過計(jì)算得到了渤海灣的半交換周期（孫健，2007）。應(yīng)用三維HAMSOM模型對渤海灣進(jìn)行模擬，得到了渤海灣的半交換時(shí)間（魏皓等，2002）。有關(guān)圍填海工程對渤海水交換的影響也有一些學(xué)者進(jìn)行研究，如利用水動力學(xué)模型和拉格朗日粒子跟蹤法，針對圍填海對渤海灣水交換的影響進(jìn)行研究，并分別從整體和局部分析得到：進(jìn)行大規(guī)模圍填海之后，渤海灣半交換時(shí)間延長，其主要原因很有可能就是因?yàn)閲詈９こ淌购乘w通道受阻（袁德奎等，2015）?；贛IKE3水動力模型耦合粒子追蹤模塊分析了近十年來渤海灣海岸工程對渤海灣水交換的影響。結(jié)果表明，隨著近十年來渤海灣沿岸圍填海工程的進(jìn)行，渤海灣水域平均水存留時(shí)間明顯增加，水交換率下降（王勇智等，2015）。

前人對水交換的研究大都使用拉格朗日水質(zhì)點(diǎn)追蹤方法，采用這種方法主要是研究水質(zhì)點(diǎn)本身的物質(zhì)交換，著眼于每個(gè)水質(zhì)點(diǎn)的位置隨時(shí)間變化的規(guī)律，重點(diǎn)描述每一個(gè)流體質(zhì)點(diǎn)自始至終的運(yùn)動過程；另外采用這種方法雖然很自然，很直觀，但實(shí)現(xiàn)起來卻非常困難，需要對成千上萬的流體質(zhì)點(diǎn)進(jìn)行跟蹤，效率比較低。而本文采用歐拉方法進(jìn)行水質(zhì)點(diǎn)的追蹤來研究污染物濃度的擴(kuò)散規(guī)律，這種方法著眼于空間點(diǎn)，主要研究空間中每一個(gè)點(diǎn)上流體運(yùn)動隨時(shí)間的變化情況。相比于拉格朗日方法，能夠使得海水充分稀釋和擴(kuò)散，可以合理的反映海水的污染物擴(kuò)散情況，更加的科學(xué)合理。本文利用三維水動力彌散水質(zhì)模型FVCOM，針對渤海灣沿岸尤其是天津近海近年來的大規(guī)模圍填?，F(xiàn)狀，對該海域保守物質(zhì)濃度輸運(yùn)及其交換能力的變化進(jìn)行數(shù)值模擬，并對其定量的研究，為渤海灣水環(huán)境的情況提供更加科學(xué)的認(rèn)知，同時(shí)為海洋工程的建設(shè)創(chuàng)造依據(jù)。

1 模型

1.1 數(shù)值模型簡介

本文采用由陳長勝研究組提出的三維海洋數(shù)值模型FVCOM（Finite Volume Coast and Ocean Model），該模型具有以下幾個(gè)特點(diǎn)：在垂直方向使用滓坐標(biāo)，可以更好的擬合渤海灣研究區(qū)域的水深情況；在水平方向采用三角形網(wǎng)格，這樣可以加密需要重點(diǎn)研究的區(qū)域，同時(shí)能夠準(zhǔn)確的擬合岸線復(fù)雜區(qū)域；另外該模型包含干/濕網(wǎng)格模塊，能夠?qū)┩康难蜎]和露出過程進(jìn)行準(zhǔn)確模擬，針對渤海灣灘涂面積大的特點(diǎn)可以保證流體守恒（CHEN et al，2006）。因此，利用FVCOM對渤海灣海域進(jìn)行數(shù)值模擬是十分合適的（李希彬等，2010）。

本文基于FVCOM所建立的數(shù)值模型由水動力模型及水交換模型這兩部分構(gòu)成。其中水動力模型提供驅(qū)動水交換模型所需的三維水動力場和湍流擴(kuò)散，水交換模型用于污染物擴(kuò)散的數(shù)值模擬（李希彬等，2012）。

1.2 模擬時(shí)期選擇

為展示天津海域圍填海工程對水交換能力的影響，本文選取2次圍填海工程規(guī)模變化較大的典型歷史時(shí)期，分別為2015年（大規(guī)模圍填海后）及2005年（大規(guī)模圍填海前），在兩個(gè)不同歷史時(shí)期期間，由于圍填海工程的大量進(jìn)行，天津海域面積減少約280 km2，渤海灣海域面積已減少近0.5%，如圖1所示，圖中灰色區(qū)域?yàn)?005-2015年期間通過圍填海工程增加的陸地區(qū)域。利用歷史數(shù)據(jù)及相應(yīng)歷史時(shí)期的海圖、高分辨率衛(wèi)星遙感圖像及水深資料等，分別建立相應(yīng)時(shí)期岸線和地形特征的三維水動力數(shù)值模型，使其能夠反應(yīng)對應(yīng)時(shí)期天津近海的水動力環(huán)境及水交換狀況。

1.3 模型計(jì)算區(qū)域及網(wǎng)格設(shè)置

本文研究的渤海灣水交換指的是該水域的水體通過對流擴(kuò)散等物理過程與周圍水域的水體進(jìn)行相互混合，因此所選的模擬計(jì)算區(qū)域是整個(gè)渤海（37°-41°N，117.5°-122.8°E），如圖2所示，以便于研究渤海灣與外界的水交換。

圖1 2005-2015年天津沿海岸線變化

圖2 模式網(wǎng)格圖

利用最新的衛(wèi)星遙感影像，將其進(jìn)行矢量化，并結(jié)合海圖資料得到岸線數(shù)據(jù)，將通過航保部獲得的水深資料進(jìn)行插值得到想要的水深數(shù)據(jù)。模型采用冷啟動，水位、流速初始條件均用零。

為了著重研究天津近岸圍填海工程對海域水交換的影響，在模型網(wǎng)格中，對渤海灣天津近海區(qū)域進(jìn)行加密，渤海中部區(qū)域網(wǎng)格相對較疏。模擬區(qū)域面積為450×500km2，區(qū)域內(nèi)最小網(wǎng)格步長為200m，最大網(wǎng)格步長為6 000 m。模型內(nèi)模時(shí)間步長20 s，外模時(shí)間步長2 s，垂向6個(gè)滓層。

為了減小開邊界對渤海灣內(nèi)研究區(qū)域的影響，將開邊界設(shè)置在渤海海峽以外，開邊界條件通過開邊界上56個(gè)點(diǎn)的調(diào)和常數(shù)進(jìn)行計(jì)算得到開邊界驅(qū)動水位。

2 模型結(jié)果驗(yàn)證

模型結(jié)果的驗(yàn)證，我們采用2015年10月18日-19日國家海洋局天津海洋環(huán)境監(jiān)測中心站在天津近海進(jìn)行的3個(gè)站位的海流監(jiān)測結(jié)果和塘沽海洋環(huán)境監(jiān)測站的驗(yàn)潮資料，將實(shí)測資料與模型計(jì)算結(jié)果相比較，繪制流向、流速曲線如下圖3所示，潮位驗(yàn)證結(jié)果如下圖4所示，圖中虛線代表觀測數(shù)據(jù)，實(shí)線為數(shù)值模擬結(jié)果。

從模擬值與實(shí)測值的對比可以看出，模型模擬結(jié)果較好，流速、流向模擬結(jié)果誤差均在允許范圍內(nèi)，潮汐模擬中振幅誤差較小，而且相位對比基本一致。因此，本文所采用的模型科學(xué)合理，可以用來模擬渤海灣的水交換情況，其模擬結(jié)果能夠正確的反映渤海灣的水動力特征。

3 水交換模擬結(jié)果及其對比

在計(jì)算污染物擴(kuò)散時(shí)，采用歐拉方法從理論上計(jì)算保守污染物擴(kuò)散，在模擬的初始時(shí)刻將渤海灣海域內(nèi)（118.95°E以西海域）保守污染物的濃度全部設(shè)置為1，灣外濃度全部設(shè)置為0（圖5），在與前文模型水動力條件相同的情況下進(jìn)行污染物擴(kuò)散的數(shù)值模擬，設(shè)置模型開始時(shí)刻進(jìn)行污染物的釋放。分別針對2005年和2015年兩段不同歷史時(shí)期的岸線特征，進(jìn)行800 d的污染物擴(kuò)散模擬，并進(jìn)行分析。

先對2015年情況下模擬計(jì)算的污染物擴(kuò)散結(jié)果進(jìn)行分析，由下圖6和圖7可知，在潮流系統(tǒng)作用下，污染物由渤海灣向外擴(kuò)散，擴(kuò)散的主要方向?yàn)槿R州灣及渤海中部，模擬800 d后，萊州灣底部部分海域的污染物濃度也達(dá)到了0.3以上。而由于渤海灣地區(qū)水深較淺、區(qū)域半封閉，潮流特征較弱，所以污染物擴(kuò)散能力較弱，模擬800 d后，仍有部分海域的污染物殘留濃度達(dá)0.6以上。渤海灣內(nèi)污染物濃度隨模擬時(shí)間而逐漸降低，約300多天時(shí)，污染物濃度降為初始濃度的一半。在模擬的后階段，其水交換率呈現(xiàn)下降趨勢。模擬區(qū)域內(nèi)污染物濃度隨時(shí)間遞減的過程中有著半月左右的周期振蕩，與潮汐大小潮的周期相同（圖8）。在大潮時(shí)，納潮量達(dá)到最大，對應(yīng)的水體體積會增加，污染物濃度會相應(yīng)減少，而小潮時(shí)納潮量減少，相應(yīng)污染物濃度出現(xiàn)增加。

圖3 潮流驗(yàn)證圖

圖4 潮位結(jié)果驗(yàn)證

圖5 污染物初始場設(shè)置

圖6 模擬擴(kuò)散800 d后污染物濃度分布（渤海2015）

圖7 模擬擴(kuò)散800 d后污染物濃度分布（渤海灣2015）

圖8 污染物濃度隨時(shí)間變化圖

為展現(xiàn)圍填海工程對污染物擴(kuò)散的影響，本文采用2015年?duì)顩r模擬的污染物擴(kuò)散結(jié)果減去采用2005年?duì)顩r模擬的結(jié)果，得到污染物擴(kuò)散能力之間的差值，如下圖9所示。模擬800 d后，相比2005年，2015年?duì)顩r下，整個(gè)渤海灣海域的污染物濃度均呈現(xiàn)增加的趨勢，即擴(kuò)散慢，水交換率降低。其中，渤海灣北部海域，尤其是天津南部的南港工業(yè)區(qū)和唐山東部附近海域的水交換率下降的更加顯著，此兩處不同年份的污染物殘留濃度變化量達(dá)到0.12以上，僅靠近渤海灣口處，污染物濃度殘留略呈現(xiàn)減小趨勢。

圖9 擴(kuò)散800 d后污染物濃度差值分布

為定量展現(xiàn)圍填海工程對污染物擴(kuò)散的影響，本文列表進(jìn)行分析，通過下表1發(fā)現(xiàn)，相比2005年，在進(jìn)行水交換模擬200、400、800 d后，采用2015年?duì)顩r下的渤海灣水交換率分別下降了4.22%、5.76%和6.99%。水交換率的下降，隨時(shí)間的延長有逐漸增大的趨勢，下降的水交換率將延長污染物擴(kuò)散時(shí)間，使渤海灣原本就較弱的擴(kuò)散能力更加降低，惡劣的海洋環(huán)境更加惡化。

表1 水交換率變化對比

根據(jù)渤海灣的地理特征、水動力場及污染物擴(kuò)散能力特征，我們將渤海灣劃分為3個(gè)區(qū)域，分別是位于渤海灣西北角的區(qū)域一，渤海灣中部區(qū)域二和渤海灣南部區(qū)域三（圖10）。

圖10 渤海灣區(qū)域劃分

其中，區(qū)域一位于渤海灣西北角，處于天津北部和河北唐山西部之間，該區(qū)域?yàn)┩棵娣e大，水深淺，水動力較弱，擴(kuò)散能力較弱；區(qū)域二位于渤海灣中部，主要為天津中部和南部以東海域，此區(qū)域由于水深較深，且從灣底至灣口正對渤海中部，受外部潮波影響直接，所以水動力條件較好，擴(kuò)散能力較強(qiáng)；區(qū)域三位于渤海灣南部，包括河北滄州、山東濱州和東營沿海海域，渤海灣底部東營老黃河口附近，由于M2分潮無潮點(diǎn)的存在，此處水動力條件較弱，污染物擴(kuò)散較慢，形成一條阻止污染物向外傳輸?shù)淖韪魩В撎幬鞑坑捎诖俗韪魩У淖韪糇饔?，只能通過阻隔帶的北部通道進(jìn)行水交換，但由于該海域水動力條件較好，整體水交換能力較強(qiáng)。

表2是兩個(gè)不同歷史時(shí)期下，渤海灣內(nèi)3個(gè)區(qū)域在不同模擬時(shí)間下的水交換率對比結(jié)果。在模擬的初始階段，區(qū)域二由于直面渤海中部海域，水交換通道暢通，因此水交換率較高；區(qū)域一由于水淺、灘涂面積大、水動力條件差，水交換率很低，模擬200 d時(shí)，水交換率仍在1%以下，尤其是圍填海工程完成后，在模擬200 d時(shí)，水交換率極低，僅達(dá)到0.27%，比圍填海前下降0.57%；區(qū)域三在模擬的初始階段，由于阻隔帶的存在，水交換率較低，在模擬200 d時(shí)，水交換率在10%以下。模擬800 d時(shí)，在持續(xù)的水動力條件作用下，3個(gè)海域的水交換率均有大幅提升，其中區(qū)域一的水交換率達(dá)到50%左右，區(qū)域二的水交換率達(dá)到70%左右，而區(qū)域三，尤其是區(qū)域三的中部和西部海域，在較好的SW-NE向潮流條件作用下，水交換率大幅提升，達(dá)到70%以上。值得注意的是，依據(jù)2015年條件模擬的3個(gè)區(qū)域的水交換率，均低于依據(jù)2005年條件模擬的水交換率，尤其是區(qū)域二，在模擬中段，水交換率甚至下降達(dá)11%，這主要是由于天津沿海的大部分圍填海工程均集中在此區(qū)域，包括天津港、臨港經(jīng)濟(jì)區(qū)、南港工業(yè)區(qū)等，大規(guī)模的圍填海工程占用了海域，使得此海域的納潮面積減少，納潮量降低，水交換率下降。區(qū)域一內(nèi)圍填海工程相對區(qū)域二較少，僅北疆電廠和濱海旅游區(qū)，其水交換率下降也較區(qū)域二少，模擬800 d后，水交換率下降6%。區(qū)域三中，圍填海工程少，其水交換率的變化也較少，尤其是模擬的中前期，其水交換率的變化在1%內(nèi)。

在對渤海灣的水交換率進(jìn)行計(jì)算和對比后，我們對渤海灣內(nèi)污染物濃度降為初始濃度一半的時(shí)間，即渤海灣海水的半交換周期進(jìn)行了計(jì)算，并繪制半交換周期圖，如圖11所示。渤海灣水體半交換周期的等值線分布與污染物殘留等值線分布大體類似，整個(gè)渤海灣水體的半交換周期表現(xiàn)從灣口向?yàn)硟?nèi)逐漸增大，仍然出現(xiàn)北部和東南部周期長、中部和西南部周期短的趨勢。同樣利用2015年?duì)顩r模擬的半交換周期結(jié)果減去采用2005年?duì)顩r模擬的結(jié)果，得到半交換周期的差值，如圖12所示，可以發(fā)現(xiàn)渤海灣北部和西部的中間位置，在大規(guī)模的圍填海工程后，水體的半交換周期是延長的，部分海域的半交換周期延長200 d以上；渤海灣南部海域的半交換周期是縮短的，其中部分海域的縮短時(shí)間可達(dá)20 d以上。

表2 渤海灣內(nèi)不同區(qū)域水交換率計(jì)算結(jié)果

圖11 渤海灣水體半交換周期分布圖

圖12 渤海灣水體半交換周期差值分布圖

計(jì)算之后，本文同樣定量的分析了渤海灣的半交換周期，并對兩個(gè)不同歷史時(shí)期不同區(qū)域的半交換周期進(jìn)行對比（見表3）。在2015年圍填海工程大規(guī)模建設(shè)的狀態(tài)下，整個(gè)渤海灣水體的半交換周期為398 d，而在2005年的狀態(tài)下，僅為366 d，半交換周期延長32 d。區(qū)域一水體交換能力最弱，其半交換周期達(dá)近800 d，由于圍填海工程的建設(shè)，半交換周期延長64 d，此區(qū)域沿岸主要為天津市北部和唐山市西部，可能也考慮到該區(qū)域的水交換能力較弱，天津市北部沿海主要規(guī)劃為旅游區(qū)和漁港區(qū)，并沒有規(guī)劃石油化工等重污染區(qū)，建議在此區(qū)域進(jìn)行生產(chǎn)活動時(shí)，要充分考慮海水的交換能力，嚴(yán)格控制污染物排放的時(shí)間和地點(diǎn)，同時(shí)在進(jìn)行相關(guān)規(guī)劃時(shí)，應(yīng)充分考慮到水交換率低、污染物難以擴(kuò)散的現(xiàn)實(shí)，進(jìn)行合理規(guī)劃，以防造成災(zāi)難性后果。區(qū)域二水體交換能力較強(qiáng)，整體半交換周期較短，圍填海工程后，其半交換周期由320 d延長至412 d。此區(qū)域的圍填海項(xiàng)目均為港口和石化工程，雖然此海域水交換能力相對較強(qiáng)，但仍需注意圍填海工程對水交換能力的負(fù)面影響，且應(yīng)密切關(guān)注，以防?；沸孤?，對生態(tài)環(huán)境造成災(zāi)難性影響和后果。區(qū)域三水體的交換能力強(qiáng)，水體半交換周期相對較短，圍填海工程對此海域的影響也較小，2015年?duì)顟B(tài)下，水體半交換周期為360 d，2005年?duì)顟B(tài)下，水體半交換周期為354 d，變化較小。

表3 渤海灣內(nèi)不同區(qū)域半交換周期計(jì)算結(jié)果

4 結(jié)論

本文以海洋數(shù)值模型為基礎(chǔ)，以圍填海工程對渤海灣水體的交換能力為關(guān)注點(diǎn)，通過三維水動力模型和歐拉擴(kuò)散方法，對渤海灣的污染物輸運(yùn)和水體半交換時(shí)間進(jìn)行了理論數(shù)值模擬，定量研究了渤海灣內(nèi)水體的水交換能力，主要得到以下結(jié)論：

（1）渤海灣海域水體的交換率總體較低，水體半交換周期達(dá)到300 d以上，在圍填海工程后，水體半交換周期延長32 d。

（2）圍填海工程后，渤海灣西部，尤其是天津沿海南部海域的水交換能力下降嚴(yán)重，其水體半交換周期延長可達(dá)200 d。渤海灣北部也有部分海域水交換周期延長達(dá)200 d。

（3）在圍填海工程最多的區(qū)域二，圍填海后，其水體的水交換率下降甚至可達(dá)10%，半交換周期延長92 d。而圍填海工程最少的區(qū)域三，其水體交換率和半交換周期的變化較少。

綜上所述，建議在日后進(jìn)行圍填海建設(shè)時(shí)，應(yīng)將水交換能力納入考慮，事先考慮圍填海工程對水體的交換能力造成的影響，優(yōu)先選擇水體交換能力較強(qiáng)的海域。且在排放污染物時(shí)，充分考慮水體的交換能力和排放時(shí)間，避免因人為因素造成難以擴(kuò)散的環(huán)境污染。

致謝：十分感謝國家海洋局天津海洋環(huán)境監(jiān)測中心站的專家老師為文章修改提出的寶貴意見。

Chen Ch Sh,Cowles G,Beardsley R C,2006 An unstructured grid,finitevolume coastal ocean model：FVCOM user mannual:Second edition.Oceanogvaphy,19(1):78-89.

Hainbucher D,Wei H,Pohlmann T,2004.Variability of the Bohai Sea circulation based on model circulation,Journal of Marine Systems,44(3-4):153-174.

Li X B,Sun X Y,Yu H M,et al,2010.Numerical Study of the Donghai Dam Impact on the Hydrodynamic Environment of Zhanjiang Bay.Marine Science Bulletin,12(01):16-29.

曹振東，2011.膠州灣水交換數(shù)值研究.中國海洋大學(xué).

何磊，2004，海灣水交換數(shù)值模擬方法研究.天津大學(xué).

李靜，宋軍，牟林，等，2015.冬季大風(fēng)影響下的渤黃海水交換特征.海洋通報(bào)，34（06）：647-656.

李希彬，孫曉燕，牛福新，等，2012.半封閉海灣的水交換數(shù)值模擬研究.海洋通報(bào)，31（03）：248-254.

劉哲，2004.膠州灣水體交換與營養(yǎng)鹽收支過程數(shù)值模型研究.中國海洋大學(xué).

彭輝，2013.象山港水交換數(shù)值研究.浙江大學(xué).

孫健，2007.海灣、近岸海域水交換研究的關(guān)聯(lián)矩陣方法及應(yīng)用.天津大學(xué).

王勇智，吳頔，石洪華，等，2015.近十年來渤海灣圍填海工程對渤海灣水交換的影響.海洋與湖沼，（3）：471-480.

魏皓，田恬，周峰，等，2002.渤海水交換的數(shù)值研究：水質(zhì)模型對半交換時(shí)間的模擬.青島海洋大學(xué)學(xué)報(bào)，32（4）：519-525.

袁德奎，李廣，王道生，等，2015.圍填海工程對渤海灣水交換能力影響的數(shù)值模擬.天津大學(xué)學(xué)報(bào)（自然科學(xué)與工程技術(shù)版），48（7）：605-613.

張坤，2014.欽州灣茅尾海納潮量和水交換能力的數(shù)值模擬研究.中國海洋大學(xué).

Effects of Tianjin reclamation projects on the Bohai Bay water exchange

JIN Yu-dan1,2,ZHANG Qiu-feng1,2,LI Xi-bin2,WANG Lu-ning2,YE Feng-juan2

（1.Shanghai Ocean University,Shanghai 201306,China;2.Tianjin Marine Environmental Monitoring Central Station,SOA,Tianjin 300457,China)

Based on the FVCOM（the Finite-Volume Coastal Ocean Model),numerical model of Bohai bay water exchange was established,and the observed data of tide and currents was used to validate the simulated result.The simulated results of the Bohai bay water exchange showed that half period of water exchange in Bohai Bay extends from 300 days to 325 days,when Tianjin reclamation projects are established.In western Bohai Bay,the water exchange rate even reduces by 10%,and half period of water exchange extends for 92 days.In southern Tianjin coastal area especially,water exchange capacity decreased significantly,and half period of water exchange extended for 200 days with the effect of reclamation projects,as well as in some areas of the northern Bohai Bay.Therefore,it was suggested that the effect of projects on water exchange and environment of coastal sea should be considered in the development of reclamation projects.

the Bohai Bay;reclamation project;effect of water exchange

P731.26

A

1001原6932（圓園17）05原園578原07

10.11840/j.issn.1001-6392.2017.05.014

2016-06-20；

2016-09-27

2012年天津市科技興海項(xiàng)目（KJXH2012-26）。

靳玉丹（1991原)，碩士研究生，主要從事物理海洋研究，電子郵箱：jinyudan2010@163.com。

（本文編輯：袁澤軼）