南澳羊嶼村圍填海拆除工程對水環境的影響

徐婉明，鄧偉鑄，倪錦鋒，鄭兆勇

(1.國家海洋局南海規劃與環境研究院，廣東 廣州 510300；2.珠江水資源保護科學研究所，廣東 廣州 510611)

1 研究背景

進入21世紀，我國沿海地區土地資源稀缺和經濟快速發展之間的矛盾日益突出。為了解決這一矛盾，沿海地區實施了一系列的圍填海工程[1-9]。目前，圍填海工程已經成為人類利用海洋空間資源，向海洋拓展生存和生產空間的一種重要手段[10]。

近幾年，隨著“保護優先、生態用海、集約用?！钡瓤茖W用海理念的推廣，圍填海工程對海洋生態環境的影響得到了國內眾多學者的關注。龔旭東等[11]以東山灣為例，從潮流、納潮量和水交換3個方面探討了半封閉海灣圍填海對水動力環境的影響。戴霞等[12]利用天津大學水流泥沙數學模型軟件模擬分析了渤海灣灣頂海域在天津港圍填海工程和曹妃甸圍填海工程聯合作用下的潮流運動特征。田艷等[13]從海域空間、水動力和生態環境3個方面分析了2000-2005、2005-2010和2010-2013年萊州灣3個階段圍填海工程對海洋環境產生的累積影響。陳培雄等[14]采用歷年現場采集調查資料成果對比的方法，從海底岸灘演變以及海洋水動力環境、水質環境和沉積物環境幾個方面綜合分析三門灣近期大型圍填海工程的海洋環境綜合效應。王勇智等[15]通過MIKE 3模型分析羅源灣3個典型圍填海時期岸線與地形條件下的納潮量和水交換能力變化。崔磊等[16]從環境的變化趨勢以及浮游植物、浮游動物和魚類種類及多樣性的變化入手，研究2010-2014年淇澳島水域南部半潮礁圍填海工程對當地水域生態系統造成的影響。張威等[17]基于Lagrangian粒子追蹤技術模擬分析前灣填海工程對整個膠州灣污染物輸運的影響規律及范圍。鄧澤林[18]運用環境指標變量動態對比分析法、灰色關聯分析法和Pearson相關分析法研究沙埕港海灣環境質量與圍填海工程的關聯性。

為了減輕圍填海對海洋生態環境的負面影響，2018年國家發布了《國務院關于加強濱海濕地保護嚴格管控圍填海的通知》(國發〔2018〕24 號)，其中圍填海歷史遺留問題區域的生態保護修復方案編制工作受到沿海各省的高度重視。南澳羊嶼村圍填海項目屬于廣東省圍填海歷史遺留問題。本文采用二維潮流模型和物質輸移模型，通過對南澳羊嶼村圍填海拆除方案實施前后水動力和水體交換情況進行數值模擬，分析研究了南澳羊嶼村圍填海拆除工程對水動力和水體交換的影響，為南澳羊嶼村圍填海項目生態修復方案的制定提供科學依據和指導。

2 工程概況

南澳羊嶼村圍填海工程項目位于廣東省汕頭市南澳島北部后江灣海域。該項目從1992年開始投入建設，至今已施工完成。圍填海的區域形成較為封閉的一個環抱式港池，口門寬度僅有10 m(見圖1)。受圍填海區的阻水作用，港池內水動力減弱，水體交換能力降低，導致水質環境容易惡化。為了改善港池內水動力情況，提高水體交換能力，擬對南澳羊嶼村圍填海項目實施拆除工程，拆除圍填海長度約100 m。南澳羊嶼村圍填海工程項目概況見圖1。

3 模型建立

3.1 模型基本方程

在笛卡爾直角坐標系下，基于Bousinesq渦黏假定和靜水壓假設，沿垂向平均的二維潮流和物質輸移基本方程表述如下。

連續方程：

(1)

動量方程：

(2)

(3)

對流擴散方程：

(4)

式中：x、y為直角坐標系坐標,m；t為時間,s；u、v為流速，m/s；η為水位,m；h為水深,m；Ax、Ay為水流紊動粘性系數,m2/s；f為科氏參量；g為重力加速度,m/s2；C為謝才系數；c為保守物質濃度,mg/L；Dx、Dy為物質紊動擴散系數,m2/s。

3.2 計算范圍及模型設置

模型計算范圍東至宮口半島，西至新津河口，包含后江水道、柘林灣和大埕灣水域，計算域東西寬約48 km，南北長約45 km(見圖2)。模型采用不規則三角形網格，空間分辨率在5～1 000 m之間。水深地形采用潮州市、汕頭市和福建省海域的海軍海道測量局電子海圖。外海邊界采用從全球潮汐模型TPXO7.2提取的潮位，考慮了M2、S2、N2、K2、K1、O1、P1、Q1、M4等9個調和分潮。保守物質初始濃度在港池內取1 mg/L，其他水域取0 。模擬時段為連續30 d，時間步長設置為30 s。

3.3 模型驗證

潮流模型驗證采用南澳島北部海域多個臨時觀測站的潮位過程和潮流過程，其中潮位觀測站2個(C1和C2)，潮流觀測站4個(V1～V4)，各觀測點分布見圖2。驗證時間為2017年7月11日14:00-7月12日15:00。各觀測站潮位和潮流驗證結果分別見圖3、4。圖3、4表明，C1和C2站計算潮位過程與實測潮位過程吻合較好，V1～V4站模擬流速、流向過程與實測數據在數值和相位上較為一致?？傮w上，本文建立的潮流模型能正確地反映項目海域的潮流運動規律，可以應用于工程水域的水動力和水體交換研究。

4 結果與分析

4.1 水動力影響分析

4.1.1 流速流向變化 圍填海拆除方案實施后，口門變寬，附近水域的流速、流向將會發生一定程度的變化。以大潮為例，圖5給出了圍填海拆除前后漲、落急流場對比圖。

圖1南澳羊嶼村圍填海工程項目概況 圖2模型計算范圍及驗證點分布圖

圖3 C1和C2觀測站潮位驗證結果

圖4 V1～V4觀測站潮流驗證結果

從圖5中可以看出，圍填海拆除100 m后，漲潮時水流經過口門后直接向長堤內側水域擴散，無需繞過堤頭填海區；落潮時長堤內側水體直接流向口門，不形成回流。漲、落潮時口門附近水域流速變化較大，長堤堤頭附近水域以流速增大為主，短堤內外兩側水域以流速減小為主，流速增加最大值0.20 m/s，流速減小最大值0.25 m/s。港池末端水域及港池外東側海域漲落急流速流向基本不變。

4.1.2 納潮量變化 納潮量既是衡量海灣開發價值的重要指標，也是反映海灣生命力的重要參數[19]，其大小影響著灣內與外部海域的水體交換強度，決定著海灣的自凈能力，對海灣的水環境和生態系統產生顯著影響[20]。常用的納潮量計算方法是統計一個漲潮或者落潮周期內通過特定斷面的海水通量[21]。

為了分析研究圍填海拆除方案實施前后港池的納潮量變化情況，本文在口門處選取一個斷面(見圖1)進行海水通量統計，其納潮量變化情況如表1所示。由表1可以看出，實施圍填海拆除方案后，港池內水域面積擴大，口門斷面平均納潮量增大至76.8×104m3，納潮量增加30%。

4.2 水體交換影響分析

水體半交換時間指某海域保守物質質量濃度經過對流擴散稀釋到初始質量濃度一半所需的時間[22]，是表征半封閉海灣水體交換能力的重要指標。圖6給出了圍填海拆除方案實施前后港池內水體半交換時間平面分布。

表1 圍填海拆除方案實施前后口門斷面納潮量變化

圖5 圍填海拆除方案實施前后漲、落急流場對比

從圖6中可以看出，由口門往港池內部水體半交換時間逐漸增加；港池末端東南角水體半交換時間最長，西南角次之，口門最短?，F狀條件下，港池內全部水域完成水體半交換時間為12.8 d；實施圍填海拆除方案后，港池內全部水域完成水體半交換時間為11.1 d，較圍填海拆除前縮短1.7 d。港池水體交換時間統計結果顯示，圍填海拆除100 m后，保守物質平均濃度下降至一半所需時間較拆除前縮短，前5 d保守物質平均濃度震蕩幅度較拆除前大，15 d內保守物質平均濃度下降速率較拆除前快。

5 討 論

圍填海拆除后，原有圍填海區退陸還海，圍填海拆除前口門附近順著堤頭填海區曲折岸線流動的潮流只能改變流路直接經過口門斷面進出港池，口門附近水域流態變化較大，從而導致流速變化也較大。另一方面，圍填海被拆除導致港池納潮量增大，口門作為潮流進出港池的唯一通道，流速相應也會增大。通過港池內各水域流速流向對圍填海拆除方案實施的響應程度看，港池末端水域受影響不大，而口門附近水域受影響較大。這主要是因為口門是半封閉港池內部與外海進行水體交換的主要通道，且被拆除的圍填海正位于此處，因此口門水域對圍填海拆除的響應較為敏感和劇烈。

從納潮量變化的角度看，圍填海拆除100 m后，口門寬度增加，流速增大，使得流入流出口門斷面的海水通量提高30%。實施圍填海拆除方案后，港池內水域面積變大，這必然導致港池納潮量隨之增大。

由于后江灣潮汐屬于不規則半日潮，潮汐作用較強，潮流為灣內的主要水動力條件，且海灣內無大的河流注入，因此南澳羊嶼村圍填海水域的水體交換能力主要受潮汐作用的影響。圍填海拆除方案實施前、后港池保守物質平均濃度變化對比見圖7。

圖7 圍填海拆除方案實施前、后港池保守物質平均濃度變化對比

由圖7可看出，港池內保守物質濃度隨時間遞減過程中有周期性的震蕩，這與潮汐的周期性變化和作用密切相關。水動力模擬結果表明，圍填海拆除方案實施后，港池內不存在因拆除前堤頭填海區的阻擋作用形成的弱回流區，港內與外海的水體交換較拆除前順暢，口門附近潮流流速增大，港池容納海水的能力大幅提高，因此港池內全部水域完成水體半交換時間縮短，水體交換速率增大。

綜上所述，圍填海拆除工程實施后，可以導致口門附近水動力條件增強，港池環境容量增大，水體交換能力提高，從而使港內水質環境和生態問題得到改善。張瑋等[23]研究也表明，環抱式港池的水體交換能力主要取決于水流流速大小，增大水體交換通道過水斷面寬度，可增加進出潮量，增強港池的水動力條件，提高其水體交換能力。因此，南澳羊嶼村圍填海拆除方案可以達到一定的生態修復目的。

6 結 論

本文通過建立二維潮流模型和物質輸移模型，對南澳羊嶼村歷史遺留圍填海拆除工程開展水環境影響研究，主要結論如下：

(1)圍填海拆除100 m后，口門附近水域流速、流向變化較大，長堤堤頭周邊水域以流速增大為主，短堤內外兩側水域以流速減小為主；港池末端水域及港池外東側海域流速流向基本不變。

(2)現狀條件下，港池納潮量為59.1×104m3；圍填海拆除方案實施后，港池內水域面積增大，納潮量顯著增加，增幅達30%。

(3)港池水體半交換時間的空間分布由口門向內部逐漸增加，其中水體半交換時間最大值出現在港池末端東南角。部分歷史圍填海被拆除后，港池內全部水域完成水體半交換時間為11.1 d，較圍填海拆除前縮短1.7 d。

(4)實施圍填海拆除工程，可以增強口門附近水動力條件，增大港池環境容量，提高水體交換能力，改善港內水質環境和生態問題，為編制南澳羊嶼村歷史遺留圍填海項目生態保護修復方案提供科學的依據和客觀指導。

本文對南澳羊嶼村圍填海拆除工程產生的水環境影響進行了模擬分析，下一步可對實施圍填海拆除方案引起的泥沙沖淤變化以及波浪要素變化開展研究，以便更全面地評估南澳羊嶼村圍填海拆除工程對附近海洋環境的影響，編制更合理有效的生態保護修復方案。