劉明，席小慧，雷利元，張笑，畢遠(yuǎn)溥

(遼寧省海洋水產(chǎn)科學(xué)研究院，遼寧 大連116023)

近年來(lái)，遼寧海岸帶圍填海工程項(xiàng)目不斷開(kāi)發(fā)，在帶來(lái)經(jīng)濟(jì)效益的同時(shí)，海灣面積大量減少，截灣取直、灣口束狹改變了海灣的潮流系統(tǒng)，海底淤積嚴(yán)重，生態(tài)功能退化，環(huán)境惡化、資源衰退，造成海灣資源嚴(yán)重破壞［1］。研究海岸帶開(kāi)發(fā)建設(shè)過(guò)程中對(duì)海域環(huán)境的影響，對(duì)如何合理開(kāi)發(fā)錦州灣海岸帶資源、有效保護(hù)海洋生態(tài)環(huán)境至關(guān)重要。

目前，有學(xué)者對(duì)錦州灣海域開(kāi)展了不同領(lǐng)域的專題研究，如對(duì)重金屬特性［2-3］、沉積環(huán)境變化趨勢(shì)［4］、沉積物生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［5-6］、懸浮物分布特性［7］、近岸海域富營(yíng)養(yǎng)化趨勢(shì)［8］、浮游生物分布規(guī)律［9］以及重金屬對(duì)海域生物的影響［10］等方面進(jìn)行了研究。在水動(dòng)力方面，對(duì)水文特征與泥沙狀況［11-12］、波浪和海流特征［13］進(jìn)行了初探;在環(huán)境容量方面，僅對(duì)其他海灣如膠州灣［14-15］等進(jìn)行了研究。但目前尚未見(jiàn)從水動(dòng)力方面研究圍填海工程對(duì)海灣水環(huán)境影響的報(bào)道。為此，本研究中作者選取2000、2005、2010年錦州灣圍填海工程引起的海灣面積及海岸線的變化信息，通過(guò)數(shù)值模擬方法計(jì)算近年來(lái)海灣納潮量、海灣水交換指標(biāo)的變化情況，分析海岸帶圍填海工程對(duì)錦州灣海域水動(dòng)力的影響，并結(jié)合水質(zhì)監(jiān)測(cè)數(shù)據(jù)，對(duì)海灣水環(huán)境變化趨勢(shì)進(jìn)行分析，以期為海灣資源能夠得到有效治理修復(fù)和保護(hù)以及海岸帶資源得到有效的開(kāi)發(fā)利用提供科學(xué)依據(jù)。

1 數(shù)值模型的建立與研究方法

1.1 數(shù)值模型的建立

1.1.1 潮流場(chǎng)數(shù)學(xué)模型 考慮到研究區(qū)域的范圍，且水深相比于平面尺度較小，因此，采用平面二維潮流數(shù)學(xué)模型進(jìn)行研究。其基本控制方程如下:

連續(xù)方程

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

其中:u、v 為x、y 方向的平均流速分量;H 為水深，H=h+ζ;f 為科氏力系數(shù)，f=2ωsin φ，ω 為地轉(zhuǎn)角速度，φ 為緯度;ve為效黏性系數(shù)，ve=vt+v，vt為紊動(dòng)黏性系數(shù);τbx、τby分別為底部切應(yīng)力在x、y方向的分量;τsx、τsy分別為表層切應(yīng)力在x、y 方向的分量。本研究計(jì)算中暫不考慮表面風(fēng)應(yīng)力的影響。

1.1.2 納潮量模型 納潮量是海灣潮流特征變化的總體反映，其大小決定海灣與外海的交換強(qiáng)度，反映了海灣的自凈能力。對(duì)納潮量的研究對(duì)于海灣整體規(guī)劃、合理開(kāi)發(fā)利用和資源修復(fù)及保護(hù)具有重要意義。海灣接納潮水的體積就是該海灣的納潮量［16-17］。通過(guò)計(jì)算錦州灣各個(gè)網(wǎng)格的面積以及網(wǎng)格點(diǎn)的高低潮時(shí)水位差，得到灣內(nèi)可容納海水的體積差，從而計(jì)算錦州灣的納潮量。計(jì)算公式［18］為

其中:P 為平均潮差條件下的納潮量;h 為平均潮差;S 為平均水域面積(即平均高潮位與平均低潮位水域面積之均值)。

1.1.3 水交換模型 水交換數(shù)學(xué)模型的基本控制方程［19］為

其中:C 為深度平均的示蹤染料濃度;u、v 分別為x、y 方向的流速(m/s);Dx、Dy 分別為x、y方向的擴(kuò)散系數(shù)，可根據(jù)Elder 公式計(jì)算;S、E 分別代表源和匯，即單位時(shí)間內(nèi)污染物的排放量和轉(zhuǎn)移量。

1.1.4 模型范圍及邊界條件 數(shù)值計(jì)算范圍為錦州灣海域(圖1)。模擬范圍在25 km，網(wǎng)格設(shè)置為200 ×150 個(gè);水深由海圖讀取并訂正到平均海平面，時(shí)間步長(zhǎng)為30 s，Manning 系數(shù)取0.026，開(kāi)邊界條件由M2、S2、K1和O1 四個(gè)分潮調(diào)和常數(shù)控制，初始水位和流速全部為零。

本計(jì)算中采用交錯(cuò)式矩形網(wǎng)格對(duì)方程進(jìn)行離散，在離散過(guò)程中采用隱式方向交替差分格式(ADI 法)［20］。根據(jù)ADI 的差分運(yùn)算規(guī)則，得到動(dòng)力學(xué)方程的差分公式。在平流項(xiàng)的處理中采用迎風(fēng)格式離散。

1.2 研究方法

選取并對(duì)比2000、2005、2010年圍填海工程引起的海灣面積及海岸線的變化信息(圖1)，通過(guò)數(shù)值模擬方法分別計(jì)算不同年份相鄰大小潮時(shí)的潮流場(chǎng)?？紤]到潮流場(chǎng)受不同分潮潮波疊加影響，不同時(shí)間的潮汐海流不同，本研究中設(shè)計(jì)方案A和方案B 兩種模擬方案。

1)方案A(現(xiàn)狀模擬):采用同年岸線邊界條件和潮汐海流邊界分別模擬錦州灣2000、2005、2010年海域潮流場(chǎng)現(xiàn)狀。

2)方案B(方案模擬):假設(shè)2000—2005年岸線未發(fā)生變化，即采用2000年岸線和2005年潮汐海流邊界，以模擬說(shuō)明2005年錦州灣岸線變化對(duì)海域潮流場(chǎng)的影響;假設(shè)2005—2010年岸線未發(fā)生變化，即采用2005年岸線和2010年潮汐海流邊界，以模擬說(shuō)明2010年錦州灣岸線變化對(duì)海域潮流場(chǎng)的影響。

2 結(jié)果

2.1 錦州灣岸線和海域面積變化

從表1可見(jiàn):2000—2005年錦州灣岸線長(zhǎng)度、灣口寬度和海域面積均有所減少;2005—2010年錦州灣岸線長(zhǎng)度有增加，而灣口寬度和海域面積均有所減少。而2000—2010年自然岸線長(zhǎng)度由19.16 km 減少到16.99 km，圍海養(yǎng)殖面積增加1.48 km2，城鎮(zhèn)建設(shè)面積增加2.02 km2，臨海工業(yè)區(qū)面積增加3.06 km2，港口面積增加15.8 km2［21］。

表1 錦州灣岸線和面積變化Tab.1 The coastline and area information in Jinzhou Bay

2.2 錦州灣海域潮流場(chǎng)模擬

通過(guò)對(duì)2000—2010年錦州灣海域潮流場(chǎng)的模擬，可以發(fā)現(xiàn)海灣流場(chǎng)整體趨勢(shì)并未產(chǎn)生明顯變化。圖2 為錦州灣海流實(shí)測(cè)數(shù)據(jù)與計(jì)算值的驗(yàn)證圖，說(shuō)明具有較好的吻合性。

圖3 為2010年模擬一個(gè)潮周期漲急時(shí)刻和落急時(shí)刻的流場(chǎng)圖，表明了錦州灣海域的潮流分布特征及其演變規(guī)律。

2.3 海灣納潮量的變化趨勢(shì)

從表2可見(jiàn):2000年錦州灣納潮量現(xiàn)狀模擬平均納潮量為2.48 ×108m3;2005年平均納潮量為2.12 ×108m3;2010年平均納潮量為1.47 ×108m3?？梢钥闯?，2000—2005年平均納潮量減少了14.52%，2005—2010年平 均 納 潮 量 減 少 了30.66%，錦州灣納潮量呈現(xiàn)逐年遞減的趨勢(shì)。

2.4 海灣水交換率的變化趨勢(shì)

根據(jù)水交換率公式計(jì)算錦州灣海域水交換模擬結(jié)果如表3所示，2000年錦州灣大潮現(xiàn)狀模擬水交換率平均為8.02%;2005年交換率平均為7.52%;2010年水交換率平均為7.15%。可以看出，2000—2005年平均水交換率減少了6.58%，2005—2010年平均水交換率減少了4.92%，錦州灣水交換率呈現(xiàn)逐年遞減的趨勢(shì)。

圖1 錦州灣2000、2005、2010年海岸線圖Fig.1 The coastline of Jinzhou Bay in 2000，2005 and 2010

圖2 錦州灣海流數(shù)據(jù)驗(yàn)證圖Fig.2 The current verification diagram in Jinzhou Bay

表2 2000—2010年錦州灣的納潮量Tab.2 Storage capacity of tidal water in Jinzhou Bay from 2000 to 2010 108m3

圖3 錦州灣海域漲急時(shí)刻和落急時(shí)刻的潮流場(chǎng)Fig.3 The critical moment of tide current in Jinzhou Bay

表3 2000—2010年錦州灣的水交換率Tab.3 The rate of water exchange in Jinzhou Bay from 2000 to 2010 %

3 討論

3.1 圍填海工程對(duì)海灣的整體影響

錦州灣海流表現(xiàn)為較強(qiáng)的往復(fù)性流動(dòng)，海流主流向?yàn)槠玈 -N 向，其中偏S 向?yàn)闈q潮流向，偏N向?yàn)槁涑绷飨?。本研究中通過(guò)模擬對(duì)照站位的潮流變化數(shù)據(jù)，說(shuō)明圍填海工程造成的岸線變化對(duì)錦州灣潮流場(chǎng)的影響。由模擬結(jié)果(表4)可知，2000—2005年圍填海工程造成錦州灣各站位流速減小，流向呈順時(shí)針變化;2005—2010年，由于受灣口北側(cè)位置的圍填海工程影響，造成北側(cè)灣口處(P1)流速變大，流向呈逆時(shí)針變化，而其他站位(P2和P3)流速均減小，流向呈順時(shí)針變化。

表4 錦州灣模擬方案計(jì)算結(jié)果Tab.4 The estimated result of the tide in Jinzhou Bay

3.2 圍填海工程對(duì)海灣納潮量的影響

根據(jù)模擬方案對(duì)比表明，2000—2005年岸線變化導(dǎo)致錦州灣納潮量平均減少了11.92%，2005—2010年岸線變化導(dǎo)致錦州灣納潮量平均減少了29.74%。近年來(lái)，納潮量減小的主要原因是由于海灣圍填海工程直接導(dǎo)致了海灣面積的減小，而納潮量的減少意味著減少了海灣海水的容量，即降低了錦州灣的環(huán)境容量，減少了容納水質(zhì)污染物的能力，將對(duì)錦州灣內(nèi)的水質(zhì)環(huán)境造成嚴(yán)重影響。

3.3 圍填海工程對(duì)海灣水交換率的影響

統(tǒng)計(jì)結(jié)果表明，2000—2005年岸線變化導(dǎo)致錦州灣水交換率平均減少了17.51%，2005—2010年岸線變化導(dǎo)致錦州灣水交換率平均減少了9.28%。近年來(lái)海灣水交換率減小的原因:一是由于灣口的寬度逐漸縮短，2000—2005年灣口寬度由10.1 km 減小到9.6 km，2005—2010年灣口寬度由9.6 km 減小到9.4 km;二是由于錦州灣的平均水交換率較低，圍填海工程造成海灣面積減小，海水受灣頂?shù)摹绊斖小弊饔迷斐伤粨Q活躍區(qū)整體向?yàn)晨谄啤Ｋ粨Q率的減少意味著海灣海水與灣外海水交換能力減弱，即降低了海灣的自凈能力，將對(duì)錦州灣內(nèi)的水質(zhì)環(huán)境造成嚴(yán)重影響。

3.4 圍填海工程對(duì)海灣水質(zhì)的影響趨勢(shì)

由表5可以看出，研究水域的化學(xué)耗氧量(COD)的質(zhì)量濃度為1.13 ～2.10 mg/L，2006年和2010年，COD 含量均超過(guò)二類水質(zhì)標(biāo)準(zhǔn);無(wú)機(jī)氮的質(zhì)量濃度為0.16 ～1.07 mg/L，超二類水質(zhì)標(biāo)準(zhǔn)主要在2006、2008、2010年;磷酸鹽的質(zhì)量濃度為0.0030 ～0.0389 mg/L，超三類水質(zhì)標(biāo)準(zhǔn)主要在2006年，其他年份未超標(biāo);油類質(zhì)量濃度為0.05 ～0.39 mg/L，2004年超過(guò)三類水質(zhì)標(biāo)準(zhǔn)，其他年份均在三類水質(zhì)標(biāo)準(zhǔn)范圍內(nèi)。表明錦州灣受到COD、無(wú)機(jī)氮、磷酸鹽和石油類不同程度的污染。

圍填海工程帶來(lái)的經(jīng)濟(jì)效益和新的發(fā)展空間將吸引更多的人們投資并開(kāi)展生產(chǎn)活動(dòng)，增大污染物的排放量，將造成水質(zhì)進(jìn)一步污染。根據(jù)對(duì)錦州灣近幾年海域水動(dòng)力的模擬，納潮量和水交換率都在逐年減小，灣內(nèi)環(huán)境容量和自凈能力均逐年降低，從水質(zhì)監(jiān)測(cè)數(shù)據(jù)趨勢(shì)分析看，各種污染物濃度逐年加強(qiáng)，證實(shí)了本研究的模擬結(jié)果。

表5 錦州灣水質(zhì)監(jiān)測(cè)數(shù)據(jù)Tab.5 The data of water quality in Jinzhou Bay

4 建議

為充分、合理地利用海灣資源，實(shí)現(xiàn)社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展，作者提出以下幾點(diǎn)建議:

1)海灣開(kāi)發(fā)利用應(yīng)充分考慮到圍填海工程對(duì)海域水動(dòng)力、海灣納潮量、環(huán)境容量指標(biāo)的影響;禁止在灣口岬角附近建立海岸工程。

2)限期整治和關(guān)閉污染嚴(yán)重的入海排污口;嚴(yán)格控制陸源污染物的入?？偭?，如有必要逐步采取并實(shí)施排污權(quán)交易方式，實(shí)施對(duì)海灣污染物的總量控制。

3)開(kāi)展海灣綜合治理與環(huán)境建設(shè)，保育自然岸線，修復(fù)已受損的海岸帶和海灣生態(tài)系統(tǒng)，充分利用生態(tài)系統(tǒng)的調(diào)控能力和采取必要的修復(fù)措施并舉，逐步恢復(fù)海灣典型的自然生態(tài)系統(tǒng)。

4)針對(duì)具有重要意義的海岸帶工程項(xiàng)目進(jìn)行合理開(kāi)發(fā)，并使其嚴(yán)格符合海洋功能區(qū)劃，嚴(yán)格評(píng)價(jià)其影響，建議在環(huán)境影響評(píng)價(jià)和海域使用論證中除了要考慮工程對(duì)周邊海域的影響，還應(yīng)考慮工程對(duì)海灣容量的影響。

5)結(jié)合錦州灣的實(shí)際情況，建立海灣保護(hù)和管理?xiàng)l例，實(shí)現(xiàn)錦州灣海域的可持續(xù)發(fā)展。

［1］劉偉，劉百橋.我國(guó)圍填?，F(xiàn)狀、問(wèn)題及調(diào)控對(duì)策［J］.廣州環(huán)境科學(xué)，2008，23(2):26 -30.

［2］Wang J，Liu R H，Yu P，et al.Study on the pollution characteristics of heavy metals in seawater of Jinzhou Bay［J］.Proscenia Environmental Sciences，2012，13:1507 -1516.

［3］Wang S F，Jia Y F，Wang S Y，et al.Fractionation of heavy metals in shallow marine sediments from Jinzhou Bay［J］.Journal of Environmental Sciences，2010，22(1):23 -31.

［4］李風(fēng)昌，鮑永恩.淺析錦州灣沉積環(huán)境變化趨勢(shì)［J］.海洋環(huán)境科學(xué)，1994，13(3):27 -36.

［5］王淑瑩，賈永鋒，王少鋒.錦州灣及附近河口沉積物中砷含量、分布及形態(tài)［J］.生態(tài)學(xué)雜志，2009，28(5):895 -900.

［6］范文宏，張博，陳靜生，等.錦州灣沉積物中重金屬污染的潛在生物毒性風(fēng)險(xiǎn)評(píng)價(jià)［J］.環(huán)境科學(xué)學(xué)報(bào)，2006，26(6):1000-1005.

［7］苗豐民，張斐然，孫樹(shù)德.錦州灣及其河口區(qū)的懸浮顆粒物質(zhì)［J］.海洋學(xué)報(bào)，1990，12(5):610 -619.

［8］趙軍，穆云俠.錦州灣近岸海域富營(yíng)養(yǎng)化趨勢(shì)的探討［J］.遼寧城鄉(xiāng)環(huán)境科技，2000，20(3):3 -5.

［9］宋倫，周遵春，王年斌，等.錦州灣夏季浮游動(dòng)物的群落特征［J］.水產(chǎn)科學(xué)，2006，25(8):408 -412.

［10］宛立.錦州灣重金屬污染狀況及其對(duì)菲律賓蛤仔的遺傳損傷［D］.大連:大連海事大學(xué)，2008.

［11］符文俠，劉國(guó)賢.錦州灣水文特征與泥沙狀況分析［J］.海洋通報(bào)，1994，13(3):42 -52.

［12］符文俠，鮑永恩，何寶林.錦州灣沉積物特征與物質(zhì)來(lái)源的探討［J］.黃渤海海洋，1994，12(2):47 -56.

［13］劉恒魁.錦州灣附近水域波浪和海流特征初探［J］.海洋湖沼通報(bào)，1995(1):1 -6.

［14］Liu Z，Wei H，Liu G S，et al.Simulation of water exchange in Jiaozhou Bay by average residence time approach［J］.Coastal and Shelf Science，2004，61(1):25 -35.

［15］Han H Y，Li K Q，Wang X L，et al.Environmental capacity of nitrogen and phosphorus pollutions in Jiaozhou Bay:Modeling and assessing［J］.Marine Pollution Bulletin，2011，63(5):262 -266.

［16］葉海桃，王義剛，曹兵.三沙灣納潮量及灣內(nèi)外的水交換［J］.河海大學(xué)學(xué)報(bào)，2007，35(1):96 -98.

［17］楊世倫，陳啟明，朱駿，等.半封閉海灣潮間帶部分圍墾后納潮量計(jì)算的商榷——以膠州灣為例［J］.海洋科學(xué)，2003，27(8):43-47.

［18］季小梅，張永戰(zhàn)，朱大奎.樂(lè)清灣近期海岸演變研究［J］.海洋通報(bào)，2006，25(1):44 -53.

［19］沈林杰，陳道信，黃惠明.溫州圍墾工程對(duì)河口水交換能力的影響［J］.海洋學(xué)研究，2009，27(4):72 -76.

［20］潘海，方國(guó)洪.海洋流動(dòng)動(dòng)力學(xué)的一種交替方向隱式二維數(shù)值模擬［J］.海洋學(xué)報(bào)，1995，17(5):21 -31.

［21］王偉偉，殷學(xué)博，吳英超，等.海岸帶開(kāi)發(fā)活動(dòng)對(duì)錦州灣環(huán)境影響分析［J］.海洋科學(xué)，2010，34(9):94 -96.