霞浦海堤工程對福寧灣潮流泥沙的影響分析

呂 江 華

(河海大學海洋學院，南京 210098)

0 引 言

海洋水利工程的建設在一定程度上會對周圍海域的水文生態環境產生影響[1]，并給生態環境和原有港口航道造成不同程度的影響[2-5]，因此分析工程建設對工程海域的水文環境影響分析具有重要意義。海堤工程是海岸的一個重要防護工程，避免海浪等一些自然現象對沿海人民的工作和生活產生危害[6]，而海堤工程的建設會改變原本岸線形狀，改變原有海洋環境，因此分析海堤工程建設對海域的環境影響的數值模擬具有重要意義[7]。影響分析主要涉及潮流和泥沙，潮流與泥沙的數值模擬起步于20世紀60年代，大體經過從一維到三維的發展[8]。由于科技的進步，90年代三維模型得到廣泛應用，三維模型具有應用范圍廣、模擬精度高等優點，能更好地模擬研究區域的實際情況。在以往影響評估的相關研究中，采用的模型有MIKE、ROMS、FVCOM等[9]，其中FVCOM相比較其他幾種模型而言，采用水平三角形網格，在近岸復雜地形模擬中更加準確。

霞浦縣小皓村位于三沙鎮與霞浦縣之間，本項目位于東海海域西南方向福寧灣海域(圖1)，海堤全長431 m，頂高程為6.00 m，后方回填形成陸域2.452 2 萬m2，形成的陸域布置為新農村建設，灘涂綠地和陸域道路。本研究運用FVCOM海洋數值模式建立了福寧灣三維潮流數學模型，通過實測水位、流速、流向和含沙量對模型進行驗證，模擬計算結果與實測資料吻合較好，驗證后的模型可用于預測霞浦海堤工程對工程海域水文環境的影響。

圖1 項目位置Fig.1 The project location

1 模型的建立

1.1 FVCOM模式簡介

本文所應用的FVCOM(Finite Volume Coast and Ocean Model)數值模式是由美國馬薩諸大學與伍茲霍爾海洋研究所聯合開發的三維海洋模式[10]，該模型垂向采用σ坐標，能夠精確模擬不規則的海底地形水深變化，水平方向采用非結構的三角網格，采用有限體積的方法，通過積分的方式計算重疊水平三角形控制體的通量來解控制方程。這種有限體積方法很好地將有限元方法處理海灣岸邊界復雜曲折的優點和有限差分方法簡單的離散結構、高效的計算效率結合起來，對于近岸河口等具有復雜地形的區域可以更好地保證質量、動量、鹽度和熱量的守恒性。FVCOM分別與POM以及ECOM-si在理想化實驗和近岸河口區域模擬中的對比結果[11]說明，岸線擬合的程度以及數值方法的選擇在邊界較為復雜的近岸河口海域的潮流數值模擬中具有重要的作用,這也正是FVCOM的優勢所在[12]。

1.2 模型建立

本模型計算范圍為整個福寧灣海域，圖2為模型計算區域和水深布置圖。計算區域劃分為高質量、非結構的三角形網格，工程海域網格進行加密處理。計算時間步長為10 s，垂向 坐標分均勻6層，計算區域三角網格節點數為73 593 個，網格單元數為143 948 個，最小網格間距為5 m。圖3為工程前后工程海域網格和水深分布圖。溫鹽初始場采用常值，模型初始水位是NOAA發布的ETOPO1地形高程數據所提供，模式從靜止狀態啟動[13]。開邊界自西南起七尺門到東北三沙角連線，采用水位強迫作為開邊界，本模型的開邊界(水邊界)水位是通過TMD(Tide Model Driver)同化衛星高度計資料得到各主要分潮調和常數后再預報得到的，考慮了M2、S2、K2、N2、O1、K1、P1和Q1這8個主要分潮。

1.3 模型驗證

模型驗證數據采用2013年6月霞浦縣三沙小皓臨時站工程水文泥沙測驗分析報告[14]。該觀測進行3個大潮海流連續站(C1、C2、C3)各28 h的海流觀測，以及一個15 d的水位連續觀測站(T1)。觀測站位的詳細信息見圖4和表1。

圖2 模型計算區域及水深布置Fig.2 Calculation area and depth

圖3 工程前后工程海域網格和水深分布圖Fig.3 The grids and depth of the calculation area before and after the project

圖4 模型驗證觀測站位圖Fig.4 The observation stations

項目站號北緯東經潮位T126°55′32．91″120°8′40．71″C126°55′24．98″120°8′26．28″海流C226°54′4．51″120°8′22．72″C326°54′58．58″120°9′43．00″

T1站水位以及C1、C2、C3站海流流速和流向觀測值與計算值時間過程對比如圖5～圖7所示(以表層為例)。從水位過程對比曲線看出，水位計算值與觀測值吻合良好，潮漲歷時與落潮歷時基本一致，潮位最大誤差0.49 m，平均絕對誤差0.023 m，相對誤差為0.47%。從流速對比圖可以看出3個觀測站中最大平均絕對誤差為0.021 m/s，最大相對誤差為9.88%，均小于10%，漲落潮時最大流速值以及出現時間點計算值與觀測值吻合良好。從流向模型對比圖看出平均流向偏差為1°～8°，流向相對流速大小而言模擬結果符合校核，高低潮及轉流時間也基本一致(表2為各站流速、流向偏差表)。

圖5 潮位計算值和實測值驗證圖Fig.5 Verification of the calculated values and the observed of the tide

圖6 流速計算值和實測值驗證圖Fig.6 Verification of the calculated values and the observed of the flow velocity

圖7 流向計算值和實測值比較圖Fig.7 Verification of the calculated values and the observed of the flow direction

站點T1潮位/mC1流速/(m·s-1)流向/(°)C2流速/(m·s-1)流向/(°)C3流速/(m·s-1)流向/(°)平均絕對誤差0．0230．0141．2100．0196．0900．0217．070相對誤差/%0．479．880．608．603．127．443．43

通過潮位和潮流的對比驗證可以看出，模擬結果與實測過程吻合較好，滿足規范要求[15]，所建模型采用的物理參數和計算參數基本合理，計算方法可靠，能夠較為準確地刻畫霞浦海堤附近海域水動力情況，可以繼續利用此模型模擬工程區附近海域潮流場變化。

2 工程對海域潮流場的影響分析

2.1 工程前后流場變化分析

為分析霞浦海堤工程建設對周邊海域潮流場的影響，依據數值模擬結果分別給出福寧灣大范圍和網格加密區域中工程前后大、小潮漲落潮流場圖(以大潮為例見圖8～圖10)。

由圖8可見，工程實施前，該海域流速普遍較低，漲急時刻形成了比較明顯的沿岸流，最大流速在笊籬山南部海域達到0.84 m/s，灣內流速基本在0.28 m/s以下。同時由于該海域淺海灘涂面積較大，在漲落潮階段有大片潮灘露出。總體而言工程建設前漲急流速大于落急流速，沿岸流的流場特征明顯。

通過工程實施前后流場對比結果圖可見，工程實施后對工程海域灣內小海域流場有較大的影響，對外圍區域影響較小。海堤附近潮流流向發生調整，新流向與新岸線走向基本一致，對大區域流場的流向也無顯著影響。

2.2 工程前后流速對比點流速變化

為了定量分析工程實施對工程海域水動力的影響，在海堤外側布設了5個流速比較點。圖11為工程海域工程前后流速對比點分布圖。

圖8 計算區域漲潮流場圖和落潮流場圖Fig. 8 The flow field of the calculation area during high tide and ebb

圖9 工程前后工程附近海域漲潮流場對比Fig. 9 Comparison of the flow field in the project seas during high tide before and after the project

圖10 工程前后工程附近海域落潮流場對比Fig.10 Comparison of the flow field in the project seas during ebb before and after the project

圖11 流速比較點分布Fig.11 Comparison points of the flow velocity

根據模型計算的工程前后大小潮期間(以大潮為例)漲落潮過程中的各流速比較點上的平均流速、流向和最大流速、流向比較數據可知，工程的實施對工程海域的流速帶來局部影響。所選的5個點中，漲潮時3號和5號點工程后流速變化較大，最大增幅可達0.024 m/s，其他3個點的流速變化較小，變化的幅度不超過0.01 m/s；在落潮時只有2號點的流速變化較大，流速增幅為0.025 m/s，其他4個點的流速變化均比較小，流速變幅不超過0.01 m/s。此外，各點流向均與工程前發生改變，由于岸線改變，在海堤附近海域流向發生逆時針偏轉，新流向大致保持與岸線平行。

總體來看，霞浦海堤工程建設只對以灣內500 m海域的水流產生影響。主要變化區域集中在海堤工程沿岸及拐角處，漲潮時刻在海堤拐角處南側海域和工程西端流速變化相對較大，落潮時只在西側拐角處流速變化較大，工程海域流向與工程前發生改變，由于岸線改變，在海堤附近海域流向發生逆時針偏轉，新流向大致保持與岸線平行，而距離工程較遠的外圍海域流速流向變化均不明顯。

3 工程附近海域沖淤變化分析

福寧灣海域懸浮泥沙觀測包括大、小潮兩次測驗，布設3個水文泥沙觀測站，站位布設位置與海流觀測站位一致(圖4)，該海域海水含沙濃度較低，實測最高含沙量為0.139 6 kg/m3，大潮總平均含沙量為0.041 1 kg/m3水域表現為自表層向底層明顯遞增的變化特性。霞浦海堤工程建設后，在一定程度上改變了工程海域的水動力環境，使海域潮流場的動態平衡格局發生改變，使該海域出現沖淤變化。

工程區海域年回淤淤強公式采用半經驗半理論公式：

(1)

含沙量計算按劉家駒提出的風浪和潮流綜合作用的挾沙能力公式：

(2)

圖12 含沙量計算值和實測值驗證圖Fig.12 Verification of the calculated values and the observed of the Sediment concentration

選取C1、C2和C3三個測站的含沙量進行模式驗證分析(以大潮垂向平均為例)。通過將計算值與實測值相比較，整體趨勢基本一致，誤差較小(圖12)。總體看來，所用模型中物理參數以及計算相關參數的選取較合理，計算方法有效，能夠很好地模擬工程海域的泥沙運動。

霞浦海堤的建設改變了岸界，進而改變了工程海域的水動力環境，從而使得工程海域沖淤環境發生了變化，根據工程前后的水動力數值模擬結果，計算出工程附近海域年沖淤強度分布(圖13)，由圖13可知，工程的建設對工程附近海域泥沙沖淤環境的影響比較微弱，影響主要集中在新岸線拐角處海域，結合該海域流場變化分布可知，泥沙濃度與流場之間存在較好的對應關系，流速增大的位置呈沖刷狀態，主要發生在海堤南側海域，最大沖刷厚度為0.012 m/a，流速減小的位置呈淤積狀態，主要發生在海堤新岸線拐角處，但其年淤積厚度均小于0.004 m/a。海堤工程對周邊海域有一定影響作用，部分區域存在輕微沖刷以及輕微淤積，總體影響范圍較小。

圖13 工程建設后年淤積厚度(正值為淤積，負值為沖刷)Fig.13 Siltation intensity after the project (The positive value is siltation, the negative value is scour)

4 結 語

基于非結構網格、干-濕判別技術和有限體積法的FVCOM海洋數值模式，建立了福寧灣海域的三維潮流場數學模型，通過實測潮位、流速、流向和含沙量對模型進行了驗證，模擬值與實測資料吻合較好，用驗證后的模型預測霞浦海堤工程對該海域潮流場和沖淤環境的影響，得到下列主要結論。

(1)福寧灣海域屬于正規半日潮型的強潮海域，潮差較大。工程建設前各站位因地形的影響，漲落潮表現不一致，整體體現漲潮流速大于落潮流速，小潮期間因漲潮影響減落。港灣內觀測海域的潮流系潮溝和岸形制約，沿著岸線方向漲落潮。

(2)工程實施前，漲潮總體流向為外海向隴頭小皓村海岸線，落潮反之。工程實施后，海堤附近潮流流向與工程前對比發生調整，由于岸線改變，在海堤附近海域流向發生逆時針偏轉，新流向與新岸線走向基本一致。

(3)工程的實施對工程海域的流速帶來局部影響。所選的5個流速比較點中，漲潮時3號點和5號點即海堤東西兩側拐角處工程后流速變化較大，最大增速為0.024 m/s，其他3個點的流速變化較小，變化的幅度不超過0.01 m/s；在落潮時只有2號點的流速變化較大，流速變幅達到0.025 m/s，其他4個點的流速變化均比較小，流速變幅不超過0.01 m/s。工程的實施對外圍海域流場無顯著影響。

(4)工程建設后，流速增大的位置呈沖刷狀態，主要發生在海堤南側海域，最大沖刷厚度為0.012 m/a，流速減小的位置呈淤積狀態，主要發生在海堤新岸線拐角處，但其年淤積厚度均小于0.004 m/a。海堤工程對周邊海域有一定影響作用，部分區域存在輕微沖刷以及輕微淤積，總體影響范圍較小。

□

參考文獻：

[1] 具杏祥,蘇學靈. 水利工程建設對水生態環境系統影響分析[J]. 中國農村水利水電, 2008,(7):8-11.

[2] 曹宇峰,林春梅,余麒祥,等. 簡談圍填海工程對海洋生態環境的影響[J].海洋開發與管理,2015,32(6):85-88.

[3] 謝 挺,胡益峰,郭鵬軍. 舟山海域圍填海工程對海洋環境的影響及防治措施與對策[J]. 海洋環境科學,2009,28(s1):105-108.

[4] Zhang W, Liu R, Qian W, et al. Influence of large-scale coastal engineering on hydrodynamics and sediment transport[J]. Journal of Waterway & Harbor, 2014.

[5] Gui F U, Jiufa L I, Dai Z, et al. Study of tidal flow numerical simulation on the reclamation project of nanhuizui beach[J]. Transactions of Oceanology & Limnology, 2007.

[6] 劉 寧,孫東亞,黃世昌,等. 風暴潮災害防治及海堤工程建設[J]. 中國水利, 2008,(5):9-13.

[7] 蔡 磊, 王燕燕,劉西平,等. 天津某海堤工程對水動力環境影響數值模擬[J]. 港工技術, 2017,54(1):17-22.

[8] 遠 航,于定勇. 潮流與泥沙數值模擬回顧及進展[J]. 海洋科學進展, 2004,22(1):97-106.

[9] 鄭沛楠,宋 軍,張芳苒,等. 常用海洋數值模式簡介[J]. 海洋預報,2008,25(4):108-120.

[10] Chen C, Liu H, Beardsley R C. An Unstructured Grid, Finite-Volume, Three-Dimensional, Primitive Equations Ocean Model: Application to Coastal Ocean and Estuaries[J]. Journal of Atmospheric & Oceanic Technology, 2003,20(20):159-186.

[11] Chen C, Huang H, Beardsley R C, et al. A finite volume numerical approach for coastal ocean circulation studies: Comparisons with finite difference models[M]. Journal of Geophysical Research: Oceans (1978-2012), 2007:83-87.

[12] 宋德海,鮑獻文, 朱學明.基于FVCOM的欽州灣三維潮流數值模擬[J]. 熱帶海洋學報,2009,28(2):7-14.

[13] 吳修廣, 劉光生, 程文龍. 基于FVCOM的杭州灣三維泥沙數值模擬[J]. 水利水運工程學報, 2011(4):86-96.

[14] 福建中海檢測技術有限公司. 霞浦縣三沙小皓二級漁港工程水文泥沙測驗分析報告[Z]. 2013.

[15] JTS/T 231-2-2010，海岸與河口潮流泥沙模擬技術規程[S].