武雅潔， 高曉紅， 梁丙臣， 張 敬(. 中國海洋大學工程學院，山東 青島 26600； 2. 中國海洋大學山東省海洋工程重點實驗室，山東 青島 26600；. 中國石油大學(華東)石油工程學院，山東 青島 266580)

?

劉公島海域潮流泥沙數(shù)值模擬研究?

武雅潔1,2， 高曉紅1， 梁丙臣1,2， 張 敬3
(1. 中國海洋大學工程學院，山東 青島 266100； 2. 中國海洋大學山東省海洋工程重點實驗室，山東 青島 266100；3. 中國石油大學(華東)石油工程學院，山東 青島 266580)

本文以劉公島旅游碼頭擴建工程為例，探索平面二維數(shù)值模式MIKE21FM在海洋工程環(huán)境影響評價工作中的應用。在MIKE21FM的基礎上，建立了劉公島周圍海域潮流場數(shù)值模型和泥沙輸移數(shù)值模型，分析和預測劉公島旅游碼頭工程建設前后周邊海域的水動力條件和泥沙沖淤狀況。數(shù)值模擬的潮位、潮流流速和流向與實測資料吻合較好，表明MIKE21FM能精確有效地模擬劉公島及其周圍海域潮流運動過程和特征。對工程建設前后數(shù)值模擬結果比較分析，表明工程建成后對周邊海域的水動力特征和地形地貌環(huán)境影響甚微，為該工程的設計和規(guī)劃提供技術依據(jù)。模擬結果表明，MIKE21FM模式能夠模擬實際工程中的復雜流動，可以在海洋工程環(huán)境影響評價中推廣應用。

海洋工程；水動力模型；數(shù)值模擬；MIKE21FM；環(huán)境影響評價

利用數(shù)值模擬技術定量分析海洋工程對其周圍海域水動力和地形地貌的影響程度和范圍，以其有效性、靈活性和相對廉價性，受到越來越多的重視和廣泛的應用。近年來，水動力數(shù)值模擬普遍使用的方法有日本學者提出的ADI方法、漢堡大學研發(fā)的HAMSOM模型、POM模式、ECOM-si模型等，這些非商業(yè)模型有較強的開發(fā)潛力，在海洋環(huán)境評價領域用得較多。但是，他們的源程序和操作手冊都十分復雜、人機交互功能較弱，不便于普及推廣[1]。而丹麥水力研究所開發(fā)的MIKE系列軟件包括水動力、水質、泥沙、對流擴散和離子追蹤等模塊，在平面二維自由表面流數(shù)值模擬方面具有強大的功能，而且操作簡單[2]。

MIKE21FM是一個基于不規(guī)則網(wǎng)格的模型，以三角形網(wǎng)格為基礎，網(wǎng)格設計靈活且可隨意控制網(wǎng)格疏密。因此，可以在任意淺水區(qū)及海洋工程所在的重點區(qū)域局部加密網(wǎng)格，因此它在模擬岸線彎曲的情況中具有相當?shù)膬?yōu)勢。MIKE21FM模型考慮了風應力、斜壓作用以及陸地邊界徑流的影響，具有算法可靠、計算穩(wěn)定、前后處理功能強大等諸多優(yōu)點，已在全球70多個國家得到廣泛應用。近年來該軟件在國內已被成功應用于潮流場數(shù)值模擬[3-4]、洪水淹沒分析[5-6]、溫排水影響分析[7-8]、水質模擬與預測[9-10]等方面。

劉公島位于為威海灣口，將威海灣口分為南北兩口。劉公島旅游碼頭擴建工程泊位擬在現(xiàn)有劉公島旅游碼頭東側延長80m而形成。國內學者對威海灣附近海域的潮流特性和泥沙淤積作了一些研究。張偉[11]采用ECOMSED數(shù)值模型對威海灣南部楊家灣內的岸灘整治工程建設后的水動力環(huán)境及沖淤環(huán)境進行了數(shù)值模擬研究；蔡學石和王永學[12]建立了威海中心漁港的波流場、泥沙場數(shù)值計算模型，通過數(shù)值模擬研究探討了威海中心漁港修建前、后的波流場及海底沖淤情況的變化特征。但是由于劉公島特殊的地理位置，他們的模擬研究海域一般以劉公島為邊界，未涉及劉公島及其周圍海域的研究。而且更重要的是泥沙輸移僅考慮了波流共同作用，而沒有考慮風的作用。因此，本文以劉公島旅游碼頭擴建工程為例，在MIKE21FM基礎上建立潮流場和風浪流作用下的沙輸輸移數(shù)值模型，預測和分析擴建工程的建設對周邊海域的水動力環(huán)境及地形地貌環(huán)境產生的影響，探索平面二維數(shù)值模式MIKE21FM在海洋工程環(huán)境影響評價工作中的應用。

1 模型簡介

(1)模型控制方程[2]

質量守恒方程為

(1)

動量方程為

(2)

(3)

泥沙控制方程為

(4)

(2)定解條件

(5)

邊界條件：在潮灘區(qū)采用動邊界處理；水邊界采用預報潮位控制。固定邊界取法向流速為零，即：

(6)

2 水動力模擬預測和分析

2.1 計算域和網(wǎng)格設置

(1)計算域設置。本文所建立的海域數(shù)學模型計算域范圍見圖1，即為圖中A、B、C 3點以及岸線圍成的海域。模擬采用非結構三角形網(wǎng)格，用動邊界的方法對干、濕網(wǎng)格進行處理。為了清楚地反映工程實施對其附近海域水動力環(huán)境的影響，模擬中將工程區(qū)附近海域的計算網(wǎng)格進行局部加密，即D、E、F 3點與岸線圍成的海域，見圖2。為了更好的與現(xiàn)有岸線相吻合，對局部岸線進行了修正加密，整個模擬區(qū)域網(wǎng)格最小空間步長為24m。

(2)水深和岸界。

水深 選取中國人民解放軍海軍航海保證部制作的11 981(1∶4 000)、12 100(1∶250 000)海圖以及工程附近海域水深地形測量資料。

岸界 采用以上海圖中岸界、908山東省海岸線勘測資料以及工程附近海岸線勘測資料。

(3)模型水邊界輸。

開邊界 本次模擬的開邊界水位由黃渤海潮流模型提供，其開邊界潮位由下式輸入計算

(7)

其中：A0為平均海平面；fi、σi是第i個分潮(分別取M2、S2、O1和K1四個分潮)的交點因子和角速度；Hi和gi是調和常數(shù)，分別為分潮的振幅和遲角；(v0+u)i是分潮的幅角。

閉邊界 以大海域和工程周邊岸線作為閉邊界。

(4)計算時間步長。模型計算時間步長根據(jù)CFL條件進行動態(tài)調整，確保模型計算穩(wěn)定進行，最小時間步長0.8s。

(5)底床糙率和水平渦動黏滯系數(shù)。底床糙率通過曼寧系數(shù)進行控制，曼寧系數(shù)n取45～58m1/3/s。水平渦動黏滯系數(shù)采用考慮亞尺度網(wǎng)格效應的Smagorinsky(1963)公式計算水平渦黏系數(shù)。

2.2 海底沉積物采樣分析

為了解劉公島附近海域底質組分和泥沙粒徑情況，在劉公島附近海域進行了11個站位的海底沉積物采樣，站位布設見圖3，沉積物粒度的實驗室分析結果見表1。實驗室分析結果表明：(1)1～11號測站的泥沙粒徑組分主要是粉砂(0.016～0.004mm)和黏土(0.004～0.001mm)；(2)沒有測到粒徑大于2.0mm的礫石；(3)除靠近劉公島的5號測站和6號測點外，其它測站泥沙中的粘土組分以50%～70%居多，粉砂組分以25%～45%居多。

2.3 潮流數(shù)值模型及驗證

潮流場數(shù)值模型采用計算區(qū)域的潮位、潮流流速和流向的實測資料對模型進行驗證，其中潮位驗證采用2010年9月的實測資料，潮流驗證分別采用2010—2011年2次的實測資料，共設置了3個驗潮點。潮位、潮流驗證點的位置見圖4。

1#驗潮點(經(jīng)度122°08.000′E，緯度37°30.333′N)自2010年9月22日0時—2010年9月23日0時的大潮潮位驗證結果見圖5。可以看出，該海域的潮汐屬于不規(guī)則半日潮。一日潮位過程包括2個完整的漲潮、落潮過程，2次潮位過程的高低潮不等現(xiàn)象明顯，模擬結果與實測資料較吻合。

圖6～7分別為2#(經(jīng)度122°09′35.2″E，緯度37°31′7.5″N)和3#(經(jīng)度122°12′16.4″E，緯度37°28′20.1″N)潮流測點2010年12月21日10時—2010年12月22日10時的大潮流速和流向的驗證結果。由驗證結果可以看出，2#和3#站的大潮期間的流速和流向均與實測資料吻合較好。

圖1 大海域計算域及網(wǎng)格設置圖Fig.1 Computational domain and grid of large sea area

圖2 局部海域網(wǎng)格設置圖Fig.2 Computational grid of local sea area

監(jiān)測站位Stations礫石Gravel/mm砂Sand/mm粉砂Siltysand/mm黏土Clay/mm平均粒徑Averagegrainsize/mm中值粒徑Mediangrainsize/mm名稱及代號Name>2．02．0?0．0630．063?0．004<0．004DmD50106．07741．15152．7710．00410．0036粉砂質黏土203．06522．8974．0450．00180．0018粉砂質黏土300．00652．2747．7250．00370．0043黏土質粉砂400．77847．67851．5430．00340．0038粉砂質黏土5017．97663．76218．2620．01790．0275黏土質粉砂607．78472．39919．8180．01340．0175黏土質粉砂70030．50569．4950．00230．0025粉砂質黏土800．70131．8467．4590．00230．0023粉砂質黏土901．54729．57868．8750．00220．0023粉砂質黏土1001．36325．29473．3430．00170．0015粉砂質黏土1103．62833．90462．4680．00270．0027粉砂質黏土

圖3 劉公島附近海域沉積物采樣站位圖Fig.3 Sediment sampling stations near the Liugong island

圖4 潮位、潮流觀測點位置Fig.4 The locations of tidal level and tidal current

圖5 大潮1#站潮位驗證Fig.5 Validation of tidal level at 1# station

圖6 大潮2#站流速、流向驗證Fig.6 Validation of tidal current at 2# station

圖7 大潮3#站流速、流向驗證Fig.7 Validation of tidal current at 3# station

因此，本文在MIKE21FM模式基礎上建立的潮流場數(shù)值模型計算結果能較好地反映劉公島及其周圍海域潮流運動過程和特征，再現(xiàn)了該海區(qū)真實的流場情況。因此，本文所建立的數(shù)值模型和邊界設置合理，可以作為旅游碼頭擴建工程前后水動力和地形地貌沖淤環(huán)境影響預測和評價的基礎。

2.4 潮流數(shù)值模擬及分析

運用MIKE21FM對工程建設前后周圍海域潮流場變化情況進行數(shù)值模擬，得到漲、落急時流場圖如圖8～9?？梢钥闯?，由于工程量較小，工程建成后對威海灣和楊家灣的潮流場幾乎沒有影響。劉公島周圍海域的漲潮流整體還是由東往西流(如圖8)，一部分潮流在劉公島西南處轉向西北，一部分在楊家灣口轉向西南流向楊家灣內，流速普遍介于0.10～0.30m/s之間。工程附近漲潮流整體由東往西流，流速較小，小于0.15m/s。劉公島周圍海域的落潮流整體由西往東流(見圖9)，劉公島西南側海域潮流由西北往東南流，流速普遍介于0.10～0.3m/s之間。工程附近漲潮流整體由東往西流，流速比漲潮時較小，小于0.08m/s。

工程附近海域由于碼頭對潮流的阻礙作用，漲潮時，潮流自東向西流，遇到碼頭工程后，沿碼頭自北向南流，至碼頭南端向西偏轉，碼頭東、西側附近潮流與工程前相比有所減小(如圖8(a)和8(b))。落潮時，潮流與漲潮時方向相反，碼頭東、西側附近潮流與工程前相比也有所減小(如圖9(a)和9(b))，影響范圍在工程附近500m范圍內，對1km之外的海域的水動力幾乎沒有影響。

3 地形地貌沖淤數(shù)值模擬

利用沉積物取樣分析、海流觀測等方法，結合水深地形、工程地質、風速資料，運用MIKE21FM模型模擬和預測潮流、波浪(施加風)作用條件下工程建設對工程周圍海域地形地貌沖淤環(huán)境的影響。根據(jù)工程附近海域風資料的統(tǒng)計結果，將全年的大風引起的波浪與潮流、徑流共同作用于地形地貌沖淤模擬中。其中沉積物類型及粒度特征參數(shù)根據(jù)該區(qū)近期和歷史表層沉積物調查資料確定，限于篇幅，沉積物沉積和侵蝕計算公式參考文獻[2，13]。

工程建設前后工程海域的地形地貌沖淤數(shù)值模擬結果見圖10。由圖10(a)可以看出，威海灣和楊家灣整體呈微淤積狀態(tài)，淤積厚度小于4cm/a，劉公島東南側和西南側與陸地連線附近海域由于地形改變導致海域面積變窄，流速較大，處于微沖刷狀態(tài)，沖刷厚度小于6cm/a。工程附近海域東南角處于小范圍的微沖刷，程度小于4cm/a，其余附近海域處于微淤積狀態(tài)，淤積厚度小于4cm/a。比較工程建設前后工程海域的沖淤情況可以看出，工程建成后對威海灣和楊家灣整體沖淤影響很小，對工程附近200m范圍內的沖淤有一定的影響，工程建成后，擬建工程南端出現(xiàn)小范圍的微沖刷，程度小于2cm/a，工程東南側和南側淤積程度增大，最大淤積厚度小于4cm/a，比工程前增大2cm/a。隨著時間向后推移，預計沖淤幅度會越來越小，10a左右基本達到穩(wěn)定狀態(tài)，最大沖淤幅值在20cm左右。

4 結論

本文以劉公島旅游碼頭擴建工程對周邊海域的水動力環(huán)境及地形地貌環(huán)境產生的影響為主要研究目標，采用數(shù)值模擬的方法，分析工程建設后對劉公島周邊海域水動力和地形地貌環(huán)境的影響。本文首先建立了劉公島周邊海域淺水潮流數(shù)值模型，并利用歷史測量的潮位、潮流數(shù)據(jù)對模型結果進行校準和驗證，數(shù)值模型計算結果與實際觀測資料吻合較好。本文得出以下結論：

圖8 工程海域大潮期間漲急時流場圖Fig.8 Tidal current field at flood strength before and after the construction

圖9 工程海域大潮期間落急時流場圖Fig.9 Tidal current field at ebb strength before and after the construction

圖10 工程區(qū)域海域的地形地貌沖淤圖Fig.10 Tidal current field at ebb strength before and after the construction

(1)本文運用潮流數(shù)值模型對劉公島周邊海域的潮流場進行數(shù)值模擬，研究工程區(qū)域的潮流運動規(guī)律和工程建設后對周邊海域的水動力條件的影響程度。數(shù)值計算結果表明：由于工程量較小，工程建成后對威海灣和楊家灣的潮流場幾乎沒有影響，影響范圍在工程附近500m范圍內。

(2)在此基礎上，本文建立了劉公島周邊海域泥沙輸移數(shù)值模型，對風、浪、流作用下研究海域內懸沙漂移擴散過程進行數(shù)值模擬。耦合潮流模型和波浪模型進一步模擬預測工程建設引起的地形地貌環(huán)境的變化。數(shù)值計算結果表明：工程建成后對工程附近200m范圍內的沖淤有一定的影響，對整個威海灣和楊家灣地形地貌沖淤環(huán)境影響甚微。

綜上所述，由于工程規(guī)模較小，對水動力的影響1km范圍外非常小，因此對水動力和地形地貌的影響的范圍會集中在1km范圍內，對外海域影響不大。本文從水動力和地形地貌的角度分析了項目建設的可行性，為該工程的設計以及規(guī)劃提供了技術依據(jù)。本文算例的成功模擬表明MIKE21FM數(shù)值模型計算穩(wěn)定、精度高、計算結果可信，能有效地模擬實際工程中的復雜流動，滿足海洋工程應用要求，為海洋工程環(huán)境影響評價領域提供新的借鑒方法。

[1] 馮靜. MIKE21FM數(shù)值模型在海洋工程環(huán)境影響評價中的應用研究[D]. 青島: 中國海洋大學, 2011.

[2] DHI MIKE21 FM Hydrodynamic Module User Guide[M]. Denmark: DHI WATER & ENVIRONMENT, 2009.

[3] 陳雪峰, 王桂萱. MIKE21計算軟件及其在長興島海域改造工程上的應用[J]. 大連大學學報, 2007, 28(6): 94-98.

[4] 鄭志飛. Mike21模型在碼頭工程流場模擬中的應用研究[J]. 水利科技, 2010(2): 25-27.

[5] 劉紹青. 濟南市城區(qū)洪水淹沒模擬研究[D]. 濟南：山東大學, 2009.

[6] 修海峰, 吳聯(lián)志. 基于MIKE II的平原河網(wǎng)洪水演進水動力研究[J]. 南水北調與水利科技, 2012, 10(6)： 151-154.

[7] 王慶改, 戴文楠, 趙曉宏, 等. 基于Mike21FM的來賓電廠擴建工程溫排水數(shù)值模擬研究[J]. 環(huán)境科學研究, 2009, 22(3): 332-336.

[8] 嚴冰, 張娜, 趙洪波. 沿岸往復流海岸電廠平面二維溫排水數(shù)值模擬研究[J]. 水道港口, 2011, 32(4): 291-296.

[9] 王哲, 劉凌, 宋蘭蘭. Mike21在人工湖生態(tài)設計中的應用[J]. 水電能源科學, 2008, 26(5): 124-127.

[10] 李娜, 葉閔. 基于MIKE21的三峽庫區(qū)涪陵段排污口COD擴散特征模擬及對下游水質的影響[J]. 華北水利水電學院學報, 2011(1): 128-131.

[11] 張偉. 威海灣岸灘整治工程的水動力及沖淤環(huán)境研究[D]. 青島： 中國海洋大學， 2008.

[12] 蔡學石, 王永學. 波流共同作用下威海中心漁港泥沙沖淤變化數(shù)值模型研究[J]. 中國水運(下半月), 2011(12): 70-72.

[13] Bingchen Liang, Huajun Li, Dongyong Lee. Numerical Study of Three-Dimensional Suspended Sediment Transport in Waves and Currents[J]. Ocean Engineering, 2007, 34(11-12): 1569-1583.

責任編輯 陳呈超

Study on Tidal Current and Sediment of the Liugong Island Based on MIKE21FM

WU Ya-Jie1,2, GAO Xiao-Hong1, LIANG Bing-Chen1,2, ZHANG Jing3

(1. College of Engineering, Ocean University of China, Qingdao 266100, China; 2. The Key Laboratory of Ocean Engineering of Shandong Province, Ocean University of China, Qingdao 266100, China; 3. School of Petroleum Enginerring, China University of Petvoleum, Qingdao 266580,China)

Taking the construction of the tourism wharf of the Liugong island as an example, the application of MIKE21FM on environmental assessment of ocean engineering is investigated. In this work, numerical models of tidal current and sediment in the Liugong island′s adjacent waters are developed based on MIKE21FM in order to analyze and predict the impacts of the construction of the tourism wharf of the Liugong island on the hydrodynamic, geographic and geomorphic environment of its adjacent waters. The values of tidal current given by the simulation have a good agreement with observed data, which indicates that the numerical model provides a fair simulation of the coastal flow in the Liugong island and its adjacent waters. The simulation results before and after the construction shows that the construction have less influence on the hydrodynamic conditions and geomorphic environment. It reveals that the construction is feasible and it can provide some theoretical bases for planning and designing of the tourism wharf. The results of this work indicate that MIKE21FM can simulate the complex flows efficiently and effectively. It can be widely applied to forecast and estimate the oceanic hydrodynamics conditions as well as erosion and deposition features in the field of environmental assessment of ocean engineering.

ocean engineering； hydrodynamic model； numerical simulation； MIKE21FM； environmental impact assessment

國家自然科學基金項目(41302195)；中央高?；究蒲袠I(yè)務費專項項目(27R1102059A; 201213005)資助

2014-03-10；

2014-06-15

武雅潔(1980-)，女，講師。E-mail: yajiewu@ouc.edu.cn

P731.2

A

1672-5174(2015)08-107-06

10.16441/j.cnki.hdxb.20140042