山東威海濱海新城北部海域潮流數值模擬研究

王玉蓮 張坤 王振興 高興 王嫚

摘要： 山東省威海市濱海新城作為威海市的新行政中心，是威海市著力打造的宜居、宜業、宜游、宜學的現代化卓越新城。本文采用MIKE21FM平面二維數值模型研究威海市濱海新城北部海域的潮流場運動，并以分層流速流向實測資料對模型進行驗證。結果表明，研究區潮流場流速和流向數值模擬結果與實測資料變化基本一致。

關鍵詞： 潮流；MIKE21FM；數值模擬；濱海新城；山東威海

中圖分類號： U652.3 ????文獻標識碼： A ???doi：10.12128/j.issn.1672 6979.2023.03.014

引文格式： 王玉蓮，張坤，王振興，等.山東威海濱海新城北部海域潮流數值模擬研究［J］.山東國土資源，2023，39（3）：96 101. WANG Yulian， ZHANG Kun， WANG Zhenxing， et al. Numerical Simulation of Tidal in Northern Sea Area ?in Binhaixincheng in Eastern Weihai City［J］.Shandong Land and Resources，2023，39（3）：96 101.

0 引言

威海市地處山東半島最東端，北、南、東三面瀕臨黃海，北與遼東半島相對，西與煙臺市接壤，東與朝鮮半島相望［1 2］。該區地處中緯度，屬于北溫帶季風型大陸性氣候，四季變化和季風進退都較明顯，與同緯度的內陸地區相比，具有雨水豐富、年溫適中、氣候溫和的特點［3］。威海市東部濱海新城作為威海市的新行政中心，是威海市著力打造的宜居、宜業、宜游、宜學的現代化卓越新城。

研究區位于威海市濱海新城北部海域，潮汐類型主要以不正規半日潮為主，僅在成山角東北海域出現一小范圍的全日潮區域。海流主流向均為偏NW—SE向，其中漲潮流向以NW向為主，落潮流向以SE向為主。本文通過2021年10月23日—10月24日潮流實測數據進行研究區的潮流數值模擬研究。

1 潮流模型簡介

采用MIKE21FM平面二維數值模型［4］研究威海市東部濱海新城北部海域潮流場運動，MIKE21FM模型是由丹麥水力研究所開發［5］，用于對河流、湖泊、河口、海灣、海岸及海洋的水流、波浪、泥沙等的模擬［6］，在國內外都得到了廣泛應用。模型采用非結構三角網格剖分計算域［7］；采用標準Galerkin有限元法進行水平空間離散，在時間上采用顯式迎風差分格式離散動量方程與輸運方程［8］。

1.1 模型控制方程

質量守恒方程： ???ζ ?t + ? ?x （hu）+ ? ?y （hv）=0

動量方程：

?u ?t +u ?u ?x +ν ?u ?y - ? ?x （εx ?u ?x ）- ? ?y （εy ?u ?y ）-fv+ gu u2+ν2 ?CZ2H =-g ?ζ ?x

?ν ?t +u ?ν ?x +ν ?ν ?y - ? ?x （εx ?ν ?x ）- ? ?y （εy ?ν ?y ）+fu+ gν u2+ν2 ?CZ2H =-g ?ζ ?x

式中： ζ—水位；h—靜水深；H—總水深，H=h+ζ；u—x向垂向平均流速；ν—y方向垂向平均流速；g—重力加速度；f—科氏力參數（f=2ω sin φ，φ為計算海域所處地理緯度）；CZ—謝才系數，CZ= 1 n H 1 6 ，n—曼寧系數；εx、εy—x、y方向水平渦動粘滯系數。

1.2 定解條件

初始條件： ??ζ（x，y，t）｜t=t0=ζ（x，y，t0）=0 u（x，y，t）｜t=t0=ν（x，y，t）｜t=t0=0

邊界條件：開邊界采用《渤、黃、東海海洋圖集——水文分冊》中4個主要分潮（K1、O1、M2、S2）的潮汐調和常數進行預報［9］。

潮位預報方程為： ?η（t）=∑ N i=1 Ai（t）αi cos （ωit+fi+Fi（t））

式中： η（t）—水位，ωi—第i個分潮的潮汐角頻率，αi和fi分別為第i個分潮的振幅和遲角，Ai（t）和Fi（t）則分別為第i個分潮的振幅和遲角修正因子［10］。

2 計算域和網格設置

2.1 計算域設置

本項目所建立的海域數學模型計算域范圍見圖1［11］。模擬采用三角網格，用動邊界的方法對干、濕網格進行處理［12］。整個模擬區域內由10672個節點和20411個三角單元組成，最小空間步長約為20m。計算海域數值模擬網格分布見圖2。

2.2 水深和岸界

水深選取前人制作的海圖資料［13］。岸界采用以上海圖中岸界和附近海域最新遙感影像資料。

2.3 計算時間步長和底床糙率

模型計算時間步長根據CFL條件進行動態調整，確保模型計算穩定進行［14］，最小時間步長1s。底床糙率通過曼寧系數進行控制［15］，曼寧系數n取32～55m1/3/s。

2.4 水平渦動粘滯系數

采用考慮亞尺度網格效應的Smagorinsky（1963）公式計算水平渦粘系數［16］，表達式如下： ?A=cs2l2 2SijSij

式中： cs—常數；l—特征混合長度［17］；Sij為網格變形率，由Sij= 1 2 （ ?ui ?xj + ?uj ?xi ）（i，j=1，2）計算得到［18］。

3 模型驗證

為評估潮流模型的準確性，在計算域內選取2個潮位觀測站位和6個潮流觀測站位的觀測資料對潮流模型進行驗證。

其中，潮位觀測資料來源于國家海洋科學數據中心公布的威海和成山角2個站位的實測潮汐數據；潮流觀測資料為實測，觀測時間為2021年10月23日12：00～24日13：00（農歷九月十八—十九，大潮期）。驗證點位置見圖3、表1。

根據潮流觀測資料，各站位的平均流速為0.275～0.455m/s。從各站實測海流資料中，提取了大潮期間各站各層及各站垂向平均的漲、落潮流向平均流速、流向和漲、落潮流的最大流速、流向（表2），1號、2號、3號和6號站位實測海流表現為較強的往復性流動（圖4、圖5）［19］，海流主流向均為偏NW—SE向，其中偏NW向為漲潮流向，偏SE向為落潮流向，4號和5號受近岸地形影響，隨漲落潮的往復性流動不是很明顯。

6個站位水體在整體上呈垂向分布，表底層水體垂向差異較小，垂向混合較為均勻。6號站的垂向平均流速最大，1號站的垂向平均流速最低。漲潮流以NW向為主，最大漲潮流速位于5號站，可達1.518m/s；落潮流以SE向為主，最大落潮流速也于5號站，最高達1.710m/s。

3.1 潮位驗證

圖3為威海站和成山角站2個站位的潮位驗證曲線。可以看出，模擬得到的潮汐漲潮歷時、落潮歷時與實測結果基本一致，且高潮時刻和低潮時刻與實測結果的相位差不超過1小時；潮位模擬結果和實測結果的最大絕對誤差在20cm以內。潮位模擬結果和實測結果總體吻合較好，數值模擬結果可信。

3.2 潮流驗證

圖4（a）—圖4（f）分別是威海北部海域海流觀測點1號—6號站位流速、流向模擬值與實測值比較圖。可以看出，除5號站位以外，與實測結果相比，各站位潮流流速的最大誤差基本小于0.2m/s，漲潮過程、落潮過程、轉流時刻的相位誤差均小于1h，流速和流向數值模擬結果與實測資料相比基本一致［20］，吻合較好。5號站大部分時刻的實測值較凌亂，導致該站位模擬結果吻合度一般。

以上潮位和潮流驗證結果表明，潮流場模擬結果與實測結果的潮位高度、漲落潮趨勢、漲落潮流速均較為一致，模擬結果可靠，能夠較好地反映研究海域的潮流狀況［21 22］。

4 潮流計算結果分析

圖5為研究海域大潮期高潮、落急、低潮、漲急時刻潮流場分布圖。從圖中可以看出，計算海域整體基本為NW—SE方向的往復流。

漲潮時，潮流自成山頭以東海域流入，沿SE—NW方向流經研究海域，并向NW方向流出，沿岸海域基本為自E向W的順岸流。成山頭附近最大流速可達0.9～1.1m/s，研究區域內離岸海域最大流速范圍基本在0.3～0.8m/s，沿岸海域流速較小，基本小于0.2m/s。

落潮時潮流方向與漲潮時相反，沿NW—SE方向流經研究海域，并向SE方向往成山頭以東流出，沿岸海域基本為自W向E的順岸流。成山頭附近最大流速可達0.9～1.2m/s，研究區域內離岸海域最大流速范圍基本在0.36～0.85m/s，沿岸海域流速較小，基本小于0.2m/s。

高潮時和低潮時潮流處于轉向時刻，流速很小，普遍小于0.15m/s。

潮流場模擬結果受岸線地形和水深條件、觀測資料的準確性和細致程度等因素影響，同時受資料限制，本研究未考慮海表面風場對表層潮流的影響。在后續研究中可通過獲取更為細致的水深地形條件、更為準確的觀測資料，并在此基礎上進行模型的優化，可進一步提升模擬結果的準確性和可靠性。

5 結論

采用MIKE21FM平面二維數值模型對威海市東部濱海新城北部海域建立了數值模型，對水動力環境進行了細致研究，得出以下結論：

（1）威海濱海新城北部海域潮汐類型主要以不正規半日潮為主，僅在成山角東北海域出現一小范圍的全日潮區域。

（2）研究區海流主流向均為偏NW—SE向，其中漲潮流向以NW向為主，落潮流向以SE向為主，且離岸海域漲潮期間和落潮期間的流速大小基本一致。

（3）研究區海域表底層水體垂向差異較小，水體垂向混合強烈。

參考文獻：

［1］ ?王武修，竇佩軍，王裕平.青島港外航道潮流場數值預報研究［J］.中外船舶科技，2021（1）：24 29.

［2］ ?袁星芳.威海市地質環境質量評價［J］.山東國土資源，2021，37（3）：35 40.

［3］ ?王玉蓮，王振興，鐘振楠.威海市地下水防污性能評價［J］.山東國土資源，2014，30（2）：47 49，53.

［4］ ?王韜翔，徐家棟.污水尾水排放對地表水環境影響預測分析：以臺州某污水廠二期工程為例［J］.綠色科技，2020（12）：6.

［5］ ?曾凡.南寧青秀山蠟燭灣水系水動力與水質效應模擬研究［D］.長沙：中南林業科技大學，2015：1 20.

［6］ ?王靜.基于EFDC的胖頭泡蓄滯洪區洪水演進模擬研究［D］.大連：大連理工大學，2017：1 15.

［7］ ?張盼.萊州灣西南部現代沉積環境研究［D］.青島：中國海洋大學，2014：1 20.

［8］ ?馬強.我國綠色船廠構建研究［D］.大連：大連海事大學，2017：1 15.

［9］ ?王如云，李慧娟，蔣風芝.一個簡化的潮汐預報準調和分析方法［J］.河海大學學報（自然科學版），2009，37（5）：611 614.

［10］ ?何帥.茅尾海水質數值模擬及環境容量研究［D］.青島：中國海洋大學，2015：1 15.

［11］ ?呂立功.漳州LNG碼頭工程懸浮物影響預測研究［J］.船海工程，2014，43（5）：151 154.

［12］ ?蔡磊，王燕燕，朱宇新.天津港某液體化工碼頭船舶污染影響模擬研究［J］.港工技術，2021，58（2）：76 80.

［13］ ?張先鋒.靖海灣沉積動力環境研究［D］.青島：中國海洋大學，2013：1 20.

［14］ ?張立奎，吳建政，李巍然，等.萊州灣東北部人工島群建設對水動力環境影響的數值研究.第二十三屆全國水動力學研討會暨第十屆全國水動力學學術會議文集［C］.北京：海洋出版社，2011：618 624.

［15］ ?岳英潔，賀志鵬，冷星，等.萊州灣芙蓉島西側人工魚礁建設對周邊海域潮流場和水交換的影響研究［J］.海洋漁業，2022，44（1）：9 17.

［16］ ?安永寧，楊鯤，王瑩，等.MIKE21模型在海洋工程研究中的應用［J］.海岸工程，2013（3）：10.

［17］ ?陳維，毛晨浩，顧一凡，等.長江口潮波對海平面上升的響應研究［J］.中國航海，2022，45（1）：112 119.

［18］ ?朱永威.浙江海島地區灘涂圍墾資源環境效應研究與分析［D］.鄭州：華北水利水電大學，2019.

［19］ ?胡甘霖.勝利埕島油田埕北256井區地形及地質分析［C］.北京：中國石油學會，2013：1 20.

［20］ ?吳仁豪，蔡樹群，王盛安，等.大亞灣海域潮流和余流的三維數值模擬［J］.熱帶海洋學報，2007，26（3）：18 23.

［21］ ?趙洪波，王廣聚.田灣核電站海域潮流場數值模擬［J］.水道港口，2007，28（1）：5 9.

［22］ ?羅慶，劉麗紅，王雨蒙.MIKE21水動力學模型應用研究進展［J］.環境保護前沿，2020，10（4）：510 515.

Numerical Simulation of Tidal in Northern Sea Area ???in Binhaixincheng in Eastern Weihai City

WANG Yulian1， ZHANG Kun2， WANG Zhenxing3， GAO Guang4， WANG Man1

（1.No.6 Geological Brigade of Shandong Provincal Bureau of Geology and Mineral Resources， Shandong Weihai 264209， China; 2.Qingdao Zhonghaichangyang Environmental Technology Limited Corporation，Shandong Qingdao 266100，China; 3.Lu'nan Geo engineering Exploration Institute， Shandong Ji'ning 272100， China; 4.Weihai Marine and Fishery Monitoring and Disaster Reduction Center， Shandong Weihai 264209， China）

Abstract： ?As the new administrative center of Weihai city， Binhaixincheng in eastern coastal area is a modern and excellent new city which is strived to build. It is suitable for living， working， traveling and learning. In this paper， the two dimensional MIKE21FM plane numerical model has been used to study the flow field movement in the northern sea area of Binhaixing in eastern Weihai city， and the stratified flow direction measured data have been used to verify the model. It is showed that the numerical simulation results of flow velocity and flow direction in the study area are basically consistent with the measured data.

Key words： ?Tidal; MIKE21FM; numerical simulation; Binhaixincheng; eastern Weihai city