三門灣近期圍填海工程對海洋環(huán)境綜合影響分析

陳培雄 張鶴 周鑫 楊士瑛

摘要：指出了三門灣位于浙江沿海中部，為半封閉式海灣，三面群山環(huán)抱，一面臨海，海灣呈西北一東南走向，是浙江省三大半封蔽型海灣之一。2000年以來三門灣內(nèi)圍填海速度更明顯加快，其對海洋自然環(huán)境的累積影響也逐漸增大。近年來針對三門灣圍填海的影響分析多以海床演變、水動力或者是水質(zhì)環(huán)境中單一因子影響分析為主，尚未有將這些因素進(jìn)行綜合評價。采用歷年現(xiàn)場采集調(diào)查資料成果對比的方法，首次從海底岸灘演變、海洋水動力環(huán)境、海洋水質(zhì)、沉積物環(huán)境幾個方面綜合分析了三門灣近期大型圍填海工程的海洋環(huán)境綜合效應(yīng)。同時利用三門灣實(shí)際沖淤狀況驗(yàn)證了海床沖淤數(shù)值計(jì)算模型，并對海床沖淤進(jìn)行了預(yù)測。通過三門灣海床整體淤積概況、水體的流速減小程度、沉積物重金屬含量變化以及水域富營養(yǎng)化歷年和季節(jié)性變化的綜合特征分析認(rèn)為：三門灣水域富營養(yǎng)化逐年增加的成因與陸域污染源的相關(guān)性較小，主要是因?yàn)槿T灣圍填海規(guī)模的持續(xù)發(fā)展導(dǎo)致灣內(nèi)水動力環(huán)境減弱而使得水體環(huán)境的自我稀釋能力下降引起的。

關(guān)鍵詞：三門灣;填海造地工程;海洋環(huán)境;流速;富營養(yǎng)化

中圖分類號：X703

文獻(xiàn)標(biāo)識碼：A

文章編號：1674-9944（2018）4-0104-06

1 引言

海灣是海陸交界處的重要地理單元，也是人類從事海洋經(jīng)濟(jì)活動及發(fā)展旅游業(yè)的重要基地。三門灣作為浙江三大典型海灣之一，是一個西北一東南方向的半封閉式海灣，灣內(nèi)灘涂面積為295 kniz，占三門灣海域面積（775 kmz）的38%[1]。三門灣海域建國以后陸續(xù)開始高灘圍涂和港汊圍堵，至今三門灣已圍涂面積（包括堵港中的水面面積）33.2萬畝（221.6 km2），占三門灣海域總面積的28.6%，總潮灘面積的75.1%。圍涂60多年來，平均每10年占據(jù)4.7%的三門灣海域面積，12.3%的三門灣潮灘總面積。其中，2003～2013年圍涂面積約10萬畝，占已圍涂總面積的30%左右，是圍涂面積增加較快的時期，尤其是灣的頂部。2003～2013年實(shí)施和規(guī)劃了一批大規(guī)模、大面積的填海造地工程：已竣工完成的主要有下洋涂工程、蛇蟠涂工程、晏站涂工程等，填海造地面積約1000 hmz以上;正在實(shí)施的填海造地工程有洋市涂區(qū)域農(nóng)業(yè)圍墾用海規(guī)劃填海工程（406.5413 hm²）和牛山火電廠填海造地工程（108.3644hm²）[2]。另外，三門灣近期及中期圍填海規(guī)劃已經(jīng)出臺，預(yù)計(jì)到2020年大規(guī)模填海造地項(xiàng)目達(dá)13.85km²[3]。因此，需要對圍填海工程對環(huán)境的影響進(jìn)行分析和評價。

2 國內(nèi)外研究現(xiàn)狀

國內(nèi)外學(xué)者對三門灣圍填海影響進(jìn)行了研究。為了進(jìn)一步分析三門灣大型填海造地工程對海灣水動力和海床的影響，謝亞力等依據(jù)1949～2003年實(shí)測資料分析了三門灣內(nèi)灘涂沖淤特性，概述了三門灣1949～2003年的填海造地情況及其對灣內(nèi)海床演變的影響[13]。楊輝等[4]胡研究發(fā)現(xiàn)三門灣貓頭水道自然狀態(tài)下因強(qiáng)潮作用維持沖淤平衡，圍填海影響下出現(xiàn)快速淤積，與其納潮量銳減、落潮流削弱有關(guān);夏小明，謝欽春[5]研究了人類活動對三門灣地貌變化尤其是潮灘和汊道發(fā)育的影響;許雪峰等[6]以蛇蟠涂和下洋涂圍墾工程為例，討論了前后兩個圍墾工程對灣內(nèi)潮流流速的影響。LinLei等認(rèn)為海灣內(nèi)的圍填工程通常會導(dǎo)致灣內(nèi)潮流減弱，主要動力機(jī)制是由于圍填?？s小了海灣的空間尺度，減小了駐波型潮波反射距離所致[7]。Yanagi T研究表明從1968年到1983年，由于海灣面積的縮小，日本的東京灣灣口流速減弱了約20%[8]。這類現(xiàn)象在日本的Ariake灣[9]、中國福建的多個海灣[10]、浙江的象山港[11]、遼寧的錦州灣[12]等地都有發(fā)現(xiàn)。這些研究多以模型模擬的方式定性進(jìn)行三門灣圍填海的影響分析，對于定量的研究內(nèi)容較少，尤其是將海床演變、水動力和生態(tài)環(huán)境進(jìn)行綜合效益評價的還沒有。本文采用現(xiàn)場采集調(diào)查資料成果以及歷史調(diào)查研究成果對比的方法，試圖從海底地形地貌、海洋水質(zhì)環(huán)境和海洋沉積物環(huán)境等幾個方面綜合分析三門灣近期大型圍填海工程所產(chǎn)生的環(huán)境影響，以期為綜合評估海灣圍填海工程的水動力和生態(tài)環(huán)境效應(yīng)提供參考，也為海灣的環(huán)境資源管理和可持續(xù)利用提供科學(xué)的依據(jù)。

3 材料與方法

水文測驗(yàn)為多船同步，即在冬、夏兩季大、中、小潮期間對13條水文泥沙測驗(yàn)垂線實(shí)施同步周日連續(xù)觀測。①流速、流向、懸沙含沙量：按六點(diǎn)法（即表層、0.2H、0.4H、0.6H、0.8H、底層）。②水溫：大、中、小潮27 h連續(xù)觀測，取樣時間間隔2h，按六點(diǎn)法（與海流觀測同）進(jìn)行。③鹽度、懸沙顆粒分析：在大、小潮第一次漲急、漲憩、落急、落憩時取樣，觀測層次采用三點(diǎn)法。④潮位觀測：從水文泥沙觀測前6 d開始，連續(xù)觀測1個月（冬、夏兩季各1個月），采樣間隔為10 min。⑤底質(zhì)顆粒分析：在大、小潮落憩時采樣。⑥氣象（風(fēng)速、風(fēng)向、氣溫、天氣現(xiàn)象）和海況：大、中、小潮連續(xù)觀測，間隔3h（即02、05、08、11、14、17、20和23時）觀測一次。水文測驗(yàn)及同步潮位觀測站位布置見圖1。

生態(tài)環(huán)境調(diào)查：海洋環(huán)境調(diào)查過程中的樣品采集、貯存、運(yùn)輸、預(yù)處理及分析測定過程均按《海洋調(diào)查規(guī)范》GB12763 - 2007和《海洋監(jiān)測規(guī)范》GB17378 - 2007中的要求進(jìn)行。水質(zhì)要素（pH、DO、NH4+-N、重金屬等）按照《海洋監(jiān)測規(guī)范》（GB17378-2007）、《海洋調(diào)查規(guī)范》（GB/T 12763 - 2007）、《水質(zhì)分析大全》（科學(xué)技術(shù)文獻(xiàn)出版社重慶分社，1989）規(guī)定的標(biāo)準(zhǔn)進(jìn)行測試分析。沉積物要素（底質(zhì)類型、有機(jī)質(zhì)含量、重金屬等）按照《海洋監(jiān)測規(guī)范》（ GB 17378 - 2007）進(jìn)行測試分析。海洋生物要素調(diào)查方法如表1。生物體質(zhì)量分析測試按照國標(biāo)《海洋監(jiān)測規(guī)范》GB17378.6 - 2007。

4 水動力影響

4.1 海床演變與沖淤影響分析

4.1.1 數(shù)值模型

潮流是泥沙運(yùn)動的載體，岸灘和水下地形對流場的分布和變化有直接影響，而潮流是造成泥沙輸運(yùn)，導(dǎo)致海床沖淤變化，從而對地貌形態(tài)進(jìn)行再塑造的最根本原因。因此，研究工程引起的海床沖淤變化，關(guān)鍵就是水動力場和含沙量場的模擬，兩者的分布和變化是海域泥沙運(yùn)動的主導(dǎo)因素。

分析采用荷蘭Delft水利研究所開發(fā)的Delft- 3D的Delft - Flow模塊，該模型既能模擬三維水動力，也能模擬二維流場，本文采用二維模式。通過平面二維數(shù)學(xué)模型計(jì)算分析三門灣內(nèi)的圍涂工程建成前后三門灣潮流場變化。

計(jì)算的水準(zhǔn)基面為1985國家高程基準(zhǔn)面。采用直接坐標(biāo)系右手系，三門灣內(nèi)水下地形采用2013年1月的實(shí)測地形資料，灣外區(qū)域則由2007年海圖拼接而成。

數(shù)學(xué)模型計(jì)算范圍應(yīng)包含研究的對象并且工程對邊界的影響足夠小，同時邊界條件容易取得。考慮到上述因素，計(jì)算范圍為整個三門灣、海游溪及亭旁溪，外邊界雀兒岙至韭山列島一線，見圖2。計(jì)算域的總面積約3268 km²。

模型采用正交曲線網(wǎng)格，為保證數(shù)模計(jì)算的精度，在水流和地形變化梯度比較大的區(qū)域適度加密網(wǎng)格，在圍墾區(qū)域附近對網(wǎng)格進(jìn)一步細(xì)化，最小網(wǎng)格5m，總網(wǎng)格數(shù)為480×522，網(wǎng)格布置見圖3。

4.1.2 海床沖淤影響結(jié)果和討論

海床沖淤模型采用2003～2013年的三門灣實(shí)際沖淤狀況進(jìn)行驗(yàn)證?？紤]到三門灣2003～2013年10年間規(guī)模較大的圍涂有蛇蟠圍涂、晏站圍涂和下洋涂圍涂，根據(jù)圍涂合攏時間，采用2006年實(shí)測地形，計(jì)算出蛇蟠和晏站圍涂工程在2006～2013的7年間產(chǎn)生的沖淤影響;采用2006年實(shí)測地形，計(jì)算出下洋涂圍涂工程在2006～2013的期間產(chǎn)生的沖淤影響，然后將兩者線性疊加，得到三門灣海域2006～2013年7年間的沖淤變化，并與這7年間的實(shí)際沖淤情況進(jìn)行對比，見圖4和圖5。

可見，實(shí)際沖淤與計(jì)算出的沖淤分布基本一致：貓頭水道、滿山水道、蛇蟠水道等都發(fā)生了不同程度的淤積，且貓頭水道內(nèi)的深潭淤積量較大，三門核電海域的一期碼頭前沿、小深潭和大深潭都發(fā)生了較大程度的淤積，模型基本復(fù)演了驗(yàn)證區(qū)域的地形變化，沖淤特性和變化量級相當(dāng)，表明海床沖淤模型能夠較好的模擬出三門灣實(shí)際的沖淤情況。

選取3條斷面和取排水口附近的2個敏感點(diǎn)進(jìn)行驗(yàn)證，斷面和敏感點(diǎn)布置見圖6，計(jì)算值與實(shí)測值的比較見表2?？梢?，各斷面、敏感點(diǎn)計(jì)算與實(shí)測的沖淤趨勢基本一致，沖淤幅度也比較接近，表明海床沖淤模型能夠較好的模擬出三門灣實(shí)際的沖淤情況。

經(jīng)驗(yàn)證，泥沙模型中主要參數(shù)取值為：沉降幾率α為0.2，沉降速度0.0005 m/s，泥沙于密度γ=725 kg/m3，區(qū)域含沙量取本次驗(yàn)證計(jì)算冬夏季平均值。

2003～2006年沖淤分布見圖7。由圖可見三門灣近期海床以淤積為主，淤積主要發(fā)生在下洋涂南側(cè)、蛇蟠水道、貓頭水道中部、滿山水道中部、貓頭山嘴北側(cè)水域淤積幅度最高可達(dá)2.0 m以上。沖刷主要發(fā)生在局部區(qū)域，如石浦水道西側(cè)、珠門港西側(cè)及貓頭山嘴南側(cè)局部區(qū)域，沖刷幅度約小于1.0 m。2006年4月至2013年沖淤分布見圖4。由圖可見，三門灣海床2006～2013年以淤積為主，下洋涂南側(cè)、蛇蟠水道東北側(cè)、貓頭水道、滿山水道、貓頭山嘴北側(cè)沿岸，淤積幅度平均在1.0m左右，最高可達(dá)到2.0 m以上。沖刷主要發(fā)生在滿山水道局部、貓頭水道區(qū)域及瀝洋港西側(cè)，沖刷幅度基本小于1.0 m?？傮w而言，三門灣歷史上處于緩慢淤積過程。2003～2013年，由于三門灣頂部區(qū)域大范圍、大規(guī)模的圍涂造地，減少了海域的納潮量，導(dǎo)致三門灣海床以淤積為主，淤積主要發(fā)生在下洋涂南側(cè)、蛇蟠水道、貓頭水道、滿山水道、貓頭山嘴北側(cè)沿岸，淤積幅度在2.0m以上。長期海域納潮量的減小以及淤積的加重會造成海水交換能力減弱，海水自凈能力驟減，水質(zhì)、底質(zhì)環(huán)境惡化，使海灣環(huán)境功能減弱。沖刷僅為局部地區(qū)，沖刷幅度總體小于1.0 m。

4.2 潮流流速變化

三門灣潮流流速分析采用2013年夏、冬兩季在三門灣布設(shè)的13個潮流周日連續(xù)觀測站進(jìn)行分析。從圖9和圖10可以看出圍涂工程實(shí)施后所引起的漲落潮流平均速變化主要分布在貓頭水道、滿山水道、蛇蟠水道、力洋港和青山港內(nèi)。其中貓頭水道的平均流速減小約0.04～0.07 m/s，滿山水道的平均流速減小約0.02～0.04 m/s，蛇蟠水道的平均流速減小約0.04～0.14m/s，力洋港平均流速減小約0.08～0.12 m/s，青山港平均流速減小約0.06～0.11 m/s。根據(jù)實(shí)測潮流資料分析表明：三門灣海域潮差大、潮流急。冬季實(shí)測最大漲、落潮流速分別為1.28 m/s和1.39 m/s;夏季實(shí)測最大漲、落潮流速分別為1.59 m/s和1.77 m/s;三門灣海域潮流場的基本特征為：流速的季節(jié)變化為夏季略大于冬季;落潮流速普遍要大于漲潮流速，最大潮流流速皆發(fā)生在落潮流中。將2003年春、秋調(diào)查數(shù)據(jù)和2013年冬、夏季調(diào)查數(shù)據(jù)進(jìn)行對比可知：三門灣的潮流場近年來具有相似的變化規(guī)律，流速季節(jié)變化為秋季略大于春季;夏季略大于冬季。2003年春季整個測區(qū)實(shí)測最大流速為1.78 m/s，秋季最大流速為2.02 m/s，均出現(xiàn)在三門灣口，最大流速均為落潮流。圍填海工程并沒有改變?nèi)T灣潮流場的變化規(guī)律，但由數(shù)據(jù)對比可以看出，由于2013年時一些大型填海在造地工程已經(jīng)實(shí)施，使得灣內(nèi)納潮量已逐年減小，因而導(dǎo)致流速亦逐年減小。灣內(nèi)納潮量和流速的減小將會減少灣內(nèi)海水交換、降低海水自凈能力，增大了灣內(nèi)水質(zhì)和底質(zhì)污染的可能性。

5 水環(huán)境質(zhì)量

5.1 海水水質(zhì)環(huán)境

根據(jù)2013年春、夏、秋、冬四季在三門灣區(qū)域采集的水質(zhì)監(jiān)測結(jié)果分析可知：2013年懸浮物變化范圍為（12～1269）mg/L，平均值為242 mg/L，總體上懸浮物含量有所下降，垂直分布上同樣為表層低于底層。pH值、總堿度變化范圍分別為7.84～8.15和（2.05～2.88）mmol/L，平均值分別為8.03和2.30 mmol/L，pH值總體上沒有變化，而總堿度則略有上升。DO、COD和BOD含量分別為（5.04～10. 14） mg/L、（0.22～3.84）mg/L和（0.11～1.88）mg/L，平均值分別為7.52 mg/L，1.26 mg/L和0.77 mg/L，DO、COD與BODs含量有所下降。表層海水石油類濃度變化范圍在0.008 mg/L～0.046 mg/L之間，含量略有上升。硝酸鹽態(tài)氮、亞硝酸鹽態(tài)氮、銨鹽態(tài)氮和活性磷酸鹽磷的濃度變化范圍分別為（ 0.258～0.853） mg/L，（ 0.000～0.060） mg/L、（0.003～0.036） mg/L和（0.018～0.052） mg/L，四季平均值分別為0. 533mg/L，0.011 mg/L、0.011 mg/L和0.037mg/L。

將2013年數(shù)據(jù)與2005年夏、秋、冬和2016年春季四季調(diào)查結(jié)果列表進(jìn)行對比（表3），結(jié)果可見：2005年三門灣水體懸浮物變化范圍較大，變化范圍為（18～4534） mg/L，平均值為419mg/L，總體表現(xiàn)為表層含量低于底層。pH值、總堿度變化范圍分別為（7.84～8.2和（1.72～2.45） mmol/L，平均值分別為8.04和2.12mmol/L。DO、COD和BOD5含量分別為（5.02～10.82） mg/L、（0.29～3.78）mg/L和（0.14～2.91） mg/L，平均值分別為7.96 mg/L，1.45 mg/L和0.77 mg/L。DO的垂直分布在夏季表現(xiàn)為表層高、底層底的特征，其它三個季節(jié)垂直變化則比較均勻。表層海水石油類濃度變化范圍在<6.5 μg/L（低于檢測限）～29.4 μg/L之間。硝酸鹽態(tài)氮、亞硝酸鹽態(tài)氮、銨鹽態(tài)氮和活性磷酸鹽磷的濃度變化范圍分別為（0.289～0.678） mg/L，（0.000～0.055）mg/L、（0.007～0.245） mg/L和（0.006～0.075） mg/L，四季平均值分別為0.494 mg/L，0.010mg/L、0.034 mg/L和0.029 mg/L。

由表3分析可見，2005年和2013年三門灣水質(zhì)環(huán)境調(diào)查數(shù)據(jù)無顯著差異。三門灣大型的填海造地工程對于水質(zhì)環(huán)境沒有形成明顯的影響。

5.2 水體富營養(yǎng)水平評價

水體營養(yǎng)鹽含量多少對該區(qū)域環(huán)境的影響可以用富營養(yǎng)化指數(shù)來評價，本文根據(jù)營養(yǎng)狀態(tài)指數(shù)（E）值大小來表示該海域富營養(yǎng)化水平，公式如下[14]：

E=CODX DIN×DIP×106/4500

（1）

其中COD為化學(xué)需氧量，DIN為溶解態(tài)無機(jī)氮，DIP為無機(jī)磷。

當(dāng)E≥1時，表明水體已達(dá)富營養(yǎng)化。

圖11展示了此次海域調(diào)查站位分布圖。表4顯示了本次調(diào)查海域四季的海水E值及其與歷史數(shù)據(jù)的比較，由評價結(jié)果可以看出，本次調(diào)查該海域四季水體均呈現(xiàn)明顯的富營養(yǎng)化狀態(tài)，E值在2.85～15.56之間，底層水體富營養(yǎng)化程度大于表層水。在季節(jié)變化上，富營養(yǎng)化程度秋季最為嚴(yán)重，冬季次之，秋、冬季E值明顯高于春、夏兩季。與歷史資料對比可以看出，三門灣海域水體富營養(yǎng)化程度呈現(xiàn)加重趨勢，本次調(diào)查四季E值與2005～2006年度相比，僅夏季E值略低，其它三季E值也明顯升高。值得注意的是，本次調(diào)查秋、冬季E值較前兩次調(diào)查增幅明顯，尤其秋季，底層水體E值平均為15.56，是2005～2006年底層水體E值的2.39倍，遠(yuǎn)遠(yuǎn)高于歷史同期數(shù)據(jù)，反映出圍填海工程后三門灣海域不僅富營養(yǎng)化程度逐漸加重，而且季節(jié)變化差異顯著，秋季富營養(yǎng)化程度為一年中最大，冬季次之，春、夏季稍弱。這主要是因?yàn)槿T灣海域圍填海工程實(shí)施后，灣內(nèi)納潮量較低，灣內(nèi)海水流速變慢，海水交換較少，降低了海水自凈能力，降低了海域水環(huán)境容量，使得污染物富集，污染物量增大。從對應(yīng)的數(shù)據(jù)也可以看出，水體中無機(jī)氮和磷酸鹽含量均是秋季最高，而COD含量也是秋季和冬季明顯較高。

5.3 沉積物環(huán)境質(zhì)量

三門灣沉積物環(huán)境調(diào)查于2013年夏冬兩季進(jìn)行。調(diào)查結(jié)果顯示三門灣海域夏季沉積物在灣內(nèi)以粉砂和砂質(zhì)粉砂為主，灣口基本為粉砂。冬季除了灣口部分站位為泥質(zhì)外，整個三門灣的沉積物均為粉砂。夏、冬兩季沉積物中值粒徑均值為7.83 μm，夏季以灣內(nèi)區(qū)最高，冬季則主要在中部出現(xiàn)高值。含水率平均為48.57%，中部海域的沉積物含水率在夏季較高，而冬季則是灣內(nèi)較高。沉積物除銅、總鉻外，其它重金屬參數(shù)與石油類均符合海洋沉積物質(zhì)量I類標(biāo)準(zhǔn)。在夏季的沉積物調(diào)查中，夏季共有4個站位銅含量略超I類，冬季則僅有1個站位略超I類?？傘t含量調(diào)查在夏季有3個站位略超I類標(biāo)準(zhǔn)，冬季則有8個站位超I類。銅、總鉻含量最大超標(biāo)（I類）指數(shù)分別為1.14（夏季）與1.14（冬季），但均符合沉積物II類標(biāo)準(zhǔn)。需要適當(dāng)展開，分析原因。

將2013年的調(diào)查數(shù)據(jù)與2005年春秋兩季調(diào)查的歷史數(shù)據(jù)進(jìn)行比對可知：2005年三門灣沉積物中重金屬銅變化范圍（20～31.5）×10^-6，均值23.3×10^-6;鉛變化范圍（23～36）×10^-6，均值31.7×10^-66;鋅變化范圍（96～119.5）×10^-6，均值107.8×10^-6;鎘的變化范圍（0.1～0.28）×10_6，均值為0.14×10^-6，汞變化范圍（0.05～0.1）×10^-6，均值0.08×10 ^-66;砷變化范圍（10.8～14.2）×10^-6，均值11.8×10^-6。而2013年調(diào)查結(jié)果顯示，銅的變化范圍（17～40）×10^-6，均值為27×10^-6;鉛的變化范圍（24～37）×10^-6，均值為32×10^-6;鋅的變化范圍為（18～103）×10^-6，均值71×10^-6;鎘的變化范圍（0.11～0.21）×10^-6，均值0.17×10^-6，汞的變化范圍（0.03～0.06）×10^-6，均值為0.05×10^-6;砷的變化范圍為（10.0～15.9）×10^-6，均值為13.6×10^-66。

分析表明：2005年至2013年期間，沉積物重金屬含量無顯著差異。三門灣大型的填海造地工程對于沉積物重金屬的含量沒有形成明顯的影響。三門灣沿岸附近居民以養(yǎng)殖、捕撈為主要收入來源，飼料的大量投放加上魚類糞便及居民生活廢物一起分解所產(chǎn)生的氮磷營養(yǎng)物質(zhì)，導(dǎo)致海區(qū)營養(yǎng)鹽含量偏高;三門灣沿岸大型工業(yè)不多，因此該海區(qū)水質(zhì)、沉積物基本未受重金屬污染（表5）。