FVCOM在長江口水動力數值模擬中的應用

李春良，倪曉雯，梅國永

（1.山東省交通規劃設計院，山東 濟南 250031；2.山東省建筑科學研究院，山東 濟南 250031）

?

FVCOM在長江口水動力數值模擬中的應用

李春良1，倪曉雯2，梅國永2

（1.山東省交通規劃設計院，山東 濟南 250031；2.山東省建筑科學研究院，山東 濟南 250031）

本文以“大江河口濕地演變退化的評價體系”項目為背景，應用FVCOM潮流及形態動力學模型建立長江口三維潮流數值計算模型，建立了包括長江口、杭州灣及鄰近海域大范圍的三維潮流數值模型，基于Linux平臺下的并行計算使得變尺度大范圍河口地區的模擬效率得到了很大的提高。運用實測潮位、流速、流向對模型的相似性進行驗證，計算驗證結果與實測值比較吻合，模擬流場能夠比較好的反映長江口地區往復流場和口外區域順時針旋轉流特征，可以用于長江口潮流的進一步研究。

長江河口區；潮流界；三維；數值模擬

1 長江口概況

1.1 概況

長江口是長江在東海入海口的一段水域，屬于較為典型的潮汐河口，潮區界位于安徽省銅陵市和蕪湖市之間，距離長江河口約640 km，潮流界在江蘇省江陰市以下，長度約240 km。按照河口地區潮流潮汐特征，通常把上自安徽省銅陵市大通鎮，下至水下三角洲地區前緣（約東經123°）的河段稱為河口區［1］。在徐六經以下，由于科氏力的作用，落潮流偏向南、漲潮流偏向北的現象較為明顯，長江在崇明島西側被分為南支和北支，在吳淞口以下南支又被長興島和橫沙島分為兩支，即南港和北港，南港又因為九段沙的分割分為南槽和北槽，河槽自西往東呈現較有規律的分叉。最終南、北支，南、北港，南、北槽呈三級分汊、四口入海的格局［2］（圖1）。

圖1　長江口河勢

1.2 潮汐與潮流

長江河口地區屬于不正規淺海半日潮，潮汐現象主要受外海的潮流潮波影響，潮汐日不等現象較為顯著［2］。潮波進入長江河口地區后，受到岸灘、河底河床抬高和上游徑流等因素的影響，潮波在上游的變形要大于下游，上游潮位要高于下游，上游潮差要小于下游，越往上游漲潮歷時逐步縮短，落潮歷時相應逐步延長［3］。

長江河口以外海域的潮流多為順時針旋轉流，潮位和流速在時間節點上基本能夠保持一致。潮流在河口進入各個河槽以后，因受到河底和岸灘約束多呈往復流形態，水位和流速隨時間發生變化，不再同步且相位發生偏差。長江口的南支和北支相對而言，南支無論從漲潮量還是落潮量上都大于北支，南支的落潮流占主導地位，但是北支從 1958年以后，徑流量分配比逐步減少，漲潮量大于落潮量成為趨勢，演變為漲潮流占主導地位［2］。

2 FVCOM海洋數值模型

FVCOM（Finite Volume Coast and Ocean Model）是美國麻州大學陳長勝教授研究課題組編程開發的，是基于無結構三角形網格架構、有限體積、三維預報原始方程近岸海洋數值模型。其計算原理是通過在水平方向上建立非結構化三角形網格，運用有限體積法解控制方程，通過積分的運算方式計算三角形控制體單元中的通量［4］。

2.1 FVCOM主要特征

1）垂向混合系數

FVCOM的垂向混合系數是基于二階湍流閉合模型來明確的，其目的是為了盡量擺脫人為干擾因素的影響。由湍流假設理論推導形成閉合湍流模型后移植到FVCOM模型中來［4］。

該模式在垂直方向上采用了s坐標系，對于處理地形的顯著變化具有良好的適應性，這對于研究受潮汐影響的河口海岸水動力學非常有必要。

3）差分格式和模態分離

模型采用顯式水平時間差分和隱式垂向時間差分，從而在水平方向和垂直方向都有比較高的分辨率。其算法采用的是時間分裂算法，其中，外模應用二維方程求解，選取的時間步長較短；內模應用三維方程來求解，所選取的時間步長較長，采用內外模式能夠節省大量的計算時間［4］。

4）干濕判別和并行計算

干濕處理法目前已經廣泛應用于河口潮間帶的水動力模型數值模擬中，通過移動邊界在不同時間步長時的干濕來判斷，干是不過水狀態，濕是過水狀態。模型運用MPI平行計算模塊，在方便調試和運行的同時極大的節約了運算時間。

2.2 原始方程組

模式控制方程組由不可壓連續、動量、溫度、密度和鹽度方程［4～6］等組成聯合方程組：

2.3非結構網格設計

本模型在計算區域內采用的非結構化三角形網格，計算范圍內由三角形單元組成，分別有節點、中心和三邊組成，并且每個單元都不重合。把變量和等設在三角單元節點上，u、v設在三角單元中心，這樣的設定將高度、速度、溫度和研通量的計算準確性大大提高了。通過選取所有相鄰三角單元中心所圍成切面的凈通量來計算節點上的變量值，通過選取進出三角單元三邊的凈通量來計算中心點的變量值［7］。

圖2　非結構網格中參量位置示意

2.4 Linux與MPI

計算程序以 Fortran77為藍本編譯，在河海大學聯想機群網絡系統上運行計算。機群操作系統采用局域網連接和星型拓撲結構［8］，集群所有節點上的操作系統均采用RedHat Linux 9.0。

圖3　MPI的執行過程

MPI是可以被 Fortran77/C/Fortran90/C++調用的函數庫系統，其具有可移植性、高效性、靈活性和易用性的特點，目前MPICH是使用最廣泛的版本，它既可以在純 Linux環境下運行，也可以在Windows環境下運行，并且使用比較方便。本文的計算就是基于 Linux平臺，運用 MPICH進行的FVCOM并行計算。

3 數值模型的計算與驗證

3.1 模擬范圍及網格

長江口徑流和潮汐等水動力作用比較強，模型的模擬計算范圍應注意潮流界在河口地區的變動范圍和潮汐作用對河口地區的綜合影響，減小選取的邊界條件對研究區域的干擾。模型選定范圍自西往東從江陰開始，東到外海50 m等深線，南北邊界介于呂四港和舟山島南側之間，包括長江口河口區域、杭州灣及相鄰部分外海區域［9］，鑒于模型計算范圍較大，為提高運算速度和效率，對部分復雜岸邊界進行了簡化，從大范圍的角度分析長江口的三維潮流動力作用。模型的計算范圍示意以及模型范圍的地形見圖4、圖5。

圖4　模型計算范圍示意

圖5　模型計算范圍地形高程填充

圖6　計算網格

圖7　深水航道網格加密示意

模型采用 SMS軟件生成全三角形網格，共有52 085個單元，27 099個節點，將水深自動插值到網格節點上，網格尺度250～10 000 m。計算網格見圖 6。模型在上邊界、崇明島、長興島、深水航道和舟山群島附近都進行了加密，加密區域的網格尺度為250～500 m，深水航道段加密網格見圖7。由于FVCOM模式采用的是干濕判別法，所以為了節省計算時間，我們將計算范圍中的邊界和島嶼進行了整合，利于模型計算更加穩定。

3.2 地形資料及模型參數

模型的地形采用 2005年實測水下地形圖，計算基準面采用 1985年國家高程基面。采用冷啟動方式進行初始啟動，設定初始為0的水位和流速，上游邊界采用實地測量的流量數值進行控制輸入，外海邊界采用東中國海潮波數值模型提供的潮位控制輸入，在時間序列上連續模擬潮流場。結合網格尺度，我們估算時間步長為2.0 s。

FVCOM為三維大范圍海洋數值模式，底部摩阻由 brough_ud.F子代碼控制，在計算過程中，可以根據自己需要修改代碼分塊設定糙率。

3.3 模型驗證及結果分析

模型采用實測水文站統計的數據進行對比驗證，通過實地測量的流速流向值、潮位值和垂直方向的平均流速值與模型預測數據進行對比分析。

1）潮位驗證分析

對比各潮位觀測站驗證曲線，除位于北支中段的青龍港站外，絕大多數觀測點模型推算潮位值與實地觀測值比較接近，誤差低于10 cm，能有效控制在誤差范圍內。共選擇 11個潮位對比站點，整個模擬區域站點布置見圖8。

圖8　潮位驗證站

經過對驗證曲線進行分析可以得出：堡鎮、橫沙、連興港和高橋站模型推算出的潮位最高值比實地測量的略小，潮位最低值比實地測量的略大，推算的潮差比實地測量數值略小。經分析，認為主要因為南北支地形和自然條件比較復雜，個別河段河底高程起伏較大，而且FVCOM模式采用的插值地形與實測地形也必然存在一定差異，所以造成個別站點的部分數據出現誤差是正常的。此外，還有一個比較重要的原因是在模型的底摩阻系數選擇上，我們為了考慮大范圍的潮波傳播過程，在大河段采用統一底摩阻系數，這無疑也影響了模擬的精度。由于外海邊界是由東中國海潮波數值模型提供，其提供邊界與實際的真實數據必然存在的一些誤差，以上幾個原因導致個別站點小潮計算值偏高（徐六經站、南門站和石洞口站）。除此之外，其余觀測點的驗證結果比較理想。我們可以得出，FVCOM數值模型推算的結果比較可靠，與選取模擬范圍的實際潮波傳播過程相吻合［10］。

2）流速驗證分析

我們選取9個流速驗證觀測點（位置見圖9），對觀測點的表、底層以及垂直方向平均流速模擬值進行提取，與實地觀測值進行對比分析，并將流向圖與實際進行對照。

圖9　流速驗證點示意

對比分析流速驗證圖，FVCOM數值模型模擬結果中提取的站點的表層、底層及垂直方向平均流速流向與實地測量結果驗證較好，對流速和流向的變化能夠比較準確的反映。河口地區因受漲潮落潮和徑流相互作用表現為往復流形態，模擬得出的落、漲潮流速，落、漲潮歷時均與實際情況相符。北支相對于南支來說，因為北支地形地貌條件比較復雜，水量較小，模擬的精度也必然受到影響，個別流速流向與實際出現偏差［10］。

3）流場驗證分析

最后根據模型視頻輸出模式形成漲落潮流場（見圖10），通過觀察可以得出，長江河口地區南、北支，南、北港，南、北槽漲落潮的水流流向和流態變化均能清楚的呈現出來。模擬范圍河段內潮流的往復流形態和口門外順時針旋轉流形態也能夠較好的呈現。

圖10　長江口落、漲急流場

4 結 語

本文對FVCOM海洋數值模式基本理論以及模型在Linux操作系統環境下的構建和應用進行了簡要介紹，并對利用MPI并行在校園網絡平臺上進行計算的原理進行了較為詳細的闡述。

利用FVCOM數值模型對長江河口大范圍的三維潮流進行模擬，并根據提取的節點表、底以及垂向平均流速與實地測量的數據進行對比分析，模擬的結果吻合度較高，往復流形態和長江口門地區順時針旋轉流的特征明顯，說明該模型在潮流水動力特性方面驗證合理，可以用來模擬長江口潮流水動力狀況。

［1]禔李來， 李誼純. 長江口整治工程對鹽水入侵影響研究［J］. 海洋工程， 2005.

［2]丁遠靜. 長江口擬建工程對枯季水動力環境的影響研究［D］. 河海大學， 2007.

［3]鄭金海， 諸裕良. 長江河口鹽淡水混合的數值模擬計算［J］. 海洋通報， 2001， 20（4）: 1-10.

［4]朱學明. 中國近海潮汐潮流的數值模擬與研究［D］. 中國海洋大學， 2009.

［5]Haosheng Huang， Changsheng Chen， Jackson O Blanton，et al. A numerical study of tidal asymmetry in Okatee Creek， South Carolina［J］. Estuarine， Coastal and Shelf Science， 2008.

［6]Changsheng Chen， Robert C Beardsley， Geoffrey Cowles. An unstructured grid， finite-volume coastal ocean model FVCOM user manual［J］.

［7]汪守東. 基于 Lagrange追蹤的海上溢油預報模型研究［D］.大連理工大學， 2008.

［8]陳林. 基于 Linux機群的大型結構并行有限元方法研究［D］. 河海大學， 2006.

［9]謝銳. EFDC模型在長江口及相鄰海域三維水流模擬中的開發應用［J］. 水動力學研究與進展: A輯， 2010.

［10]王振華. 海平面上升對長江口水動力及污染物運動軌跡影響的數值模擬研究［D］. 大連理工大學， 2011.

Application of FVCOM in Numerical Simulation of Hydrodynamic Force at Yangtze River Estuary

Li Chunliang1，Ni Xiaowen2，Mei Guoyong2
（1.Shandong Provincial Communications Planning and Design Institute， Jinan Shandong 250031， China；2.Shandong Academy of Building Research， Jinan Shandong 250031， China）

Based on the evaluation system for the evolution and degradation of great river estuarine wetlands， 3D numerical model is built to calculate the tidal current at Yangtze River Estuary by using FVCOM tidal flow and form dynamics model. By now， 3D tidal current models have been established for Yangtze River estuary， Hangzhou Bay and the surrounding sea area within a large range. The parallel computation on Linux platform can improve the variable-metric simulation of river estuarine area within a large range. The above models are tested from the aspects of the measured tidal level， velocity and direction of current. The calculation results agree very well with the measured values. The simulated flow field shows the reversing flow filed at Yangtze River delta and the clockwise rotation characteristics outside the river estuary well， which applies to the further study of the tidal current at Yangtze River estuary.

Yangtze River estuary area； tidal current limit； three-dimensional； numerical simulation

TV148+.5

A

1004-9592（2016）03-0001-05

10.16403/j.cnki.ggjs20160301

2016-01-14

李春良（1983-），男，工程師，主要從事港口及航道工程設計和海岸動力數值模擬工作。