萊州灣東岸近岸海域規(guī)劃圍填海后沖淤演變預測

周廣鎮(zhèn), 馮秀麗, 劉 杰, 劉 瀟, 徐 芳

(中國海洋大學海洋地球科學學院, 山東 青島 266100)

萊州灣東岸近岸海域規(guī)劃圍填海后沖淤演變預測

周廣鎮(zhèn), 馮秀麗, 劉 杰, 劉 瀟, 徐 芳

(中國海洋大學海洋地球科學學院, 山東 青島 266100)

運用Mike21模型中的水動力模塊對研究區(qū)的潮流場進行了模擬, 分析研究區(qū)的沖淤趨勢; 利用經驗公式, 對波浪和潮流進行耦合, 計算出極端條件下的最大沖刷深度, 對比規(guī)劃用海實施前后的最大沖刷深度, 以對該區(qū)域的海岸防護、工程選址等提供科學依據。結果表明: 受到填海工程的影響, 龍口灣內有輕微淤積, 界河至石虎咀近岸海域的波浪和潮流受工程影響而變小, 沖刷強度變弱。

萊州灣東岸; 水動力演變; 最大沖刷深度

山東半島藍色經濟區(qū)規(guī)劃的 9個集中集約用海核心區(qū), 萊州灣有4個, 僅位于萊州灣東岸的龍口灣海洋裝備制造業(yè)集聚區(qū)一期工程, 填海面積接近 50 km2, 如此大規(guī)模的圍填海, 必定會導致萊州灣海域潮流、波浪、泥沙沖淤狀況的改變, 進而影響萊州灣東岸附近海域的沖淤變化, 以及岸線的變化。近50 a來許多學者對該區(qū)的侵蝕災害和泥沙運移做過大量研究[1-3], 研究萊州灣東岸近岸海域的沖淤變化, 對指導該區(qū)海岸防護、工程選址和沿岸地區(qū)的經濟發(fā)展規(guī)劃具有重要意義。

1 自然概況

1.1 研究區(qū)地質地貌

萊州灣東岸位于山東半島西北部, 北界為屺坶島及其連島沙壩, 東為龍口港, 向西南經招遠延伸至萊州的石虎咀(研究位置如圖1所示)。屬于穩(wěn)定的平原港灣海岸, 其平面形態(tài)呈對數螺旋形。本段海岸沿岸堤低,海灘較寬, 坡度緩, 物質細。構造上位于郯廬大斷裂東側, 為上升區(qū), 屬于膠東隆起及膠北抬凸的西北翼, 沿岸多低緩丘陵, 沙嘴、沙壩、潟 湖十分發(fā)育, 是典型的對數螺旋型砂質海岸。研究區(qū)地勢由東南向西北傾斜,入海河流主要有界河、朱橋河和王河等[4]。

1.2 工程概況

龍口灣臨港高端制造業(yè)聚集區(qū)一期工程龍口部分, 南至龍口界河以北300 m, 北至龍口恒河入?？谀蟼? 西至約-8.5 m等深線, 東至海岸線, 集中集約填海約37 km2。共有6個人工島和1個突堤式人工島填海, 人工島群和陸地間水道寬度約 500 m, 島間 NE-SW 向縱向主水道寬度也為 500 m, 北部NW-SE向橫向水道寬度300 m, 南部NW-SE向橫向水道寬度200 m。規(guī)劃用海區(qū)共形成人工岸線54.8 km,占用自然岸線3.9 km。招遠部分填海約9 km2。

圖1 填海工程位置及研究區(qū)水深地形圖Fig. 1 Location of the reclamation and water depth of the study area

2 數值模擬簡介

2.1 模型介紹及參數設置

MIKE21是丹麥水動力學研究所開發(fā)的一種平面二維數值模型, 采用擬合陸邊界、網格設計靈活且可隨意控制網格疏密的非結構三角網格剖分計算域。因其具有友好的界面、可靠的算法、強大的前后處理功能等優(yōu)點, 已經在世界70多個國家得到應用, 并得到世界公認。MIKE21FM 采用標準 Galerkin有限元法進行水平空間離散, 在時間上采用顯式迎風差分格式離散動量方程與輸運方程[5]。

采用三角形非結構網格, 開邊界處由M2, S2, O1, K1四個主要分潮提供水位開邊界條件。四分潮調和常數由渤、黃、東海海洋圖集各分潮同潮圖內插得到。為了能清楚了解本工程附近海域的潮流狀況, 將本工程附近海域進行局部加密, 加密區(qū)域計算范圍為37.447°～37.686°N, 120.012°～120.317°E。整個模擬區(qū)域(填海之前)共有7 024個節(jié)點和12 304個三角單元組成。

2.2 潮位驗證

國家海洋局北海環(huán)境監(jiān)測中心于2007年12月25日19: 00～2007年12月26日20: 00(大潮期間)在工程周邊海域分別進行了3個站位26 h單周日潮流同步連續(xù)觀測。對于半日潮海區(qū), 利用大潮期間的潮位觀測資料進行調和分析計算, 可以確保得到較為正確的結果。因此, 本次模擬選用大潮期間的潮流觀測數據作為驗證資料。驗證結果表明, 對應觀測點上潮位和潮流模擬結果與實測潮流資料吻合較好, 能夠較好地反映工程周邊海域潮流狀況。

3 數值模擬結果分析

圖 2是研究區(qū)海域大潮期間落潮中間時潮流場現狀, 潮流整體由SW往NE流, 屺姆島最大流速約為75 cm/s; 石虎咀最大流速約25 cm/s, 招遠岸段附近海域潮流流速整體介于18～35 cm/s之間。

通過數值模擬結果(圖 3)可以看出, 工程前后,研究區(qū)西北部流速大小基本不變; 受到工程影響,龍口灣流速減小 4～8 cm/s, 填海工程招遠部分至石虎咀的近岸區(qū)域流速減小4～12 cm/s, 工程形成的半圓形區(qū)域, 流速減小 24 cm/s; 填海工程西側由于挑流作用流速變大 0～8 cm/s, 填海工程招遠部分與岸之間的水道以及兩工程之間的水道, 流速增大幅度在4 cm/s左右。

4 極限沖刷深度預測計算

4.1 計算原理

海岸地區(qū)大規(guī)模的泥沙輸移大多是波浪和水流共同作用下完成的。根據力學平衡原理, 將波浪和潮流在海底產生的切應力耦合, 與沉積物顆粒臨界起動剪切力相比較, 如果切應力大于臨界起動剪切力,研究區(qū)海床將繼續(xù)沖刷, 然后根據沖刷后的水深重新計算波浪要素, 計算浪流切應力, 再次比較, 直到切應力等于臨界起動剪切力, 就停止計算, 保留下此點的水深, 即為極限沖刷深度。

波流共同作用下的床面切應力:

圖2 工程前后研究區(qū)海域潮流場對比Fig. 2 Comparison of tide current field of the study area

圖3 填海前后流速變化圖Fig.3 Flow velocity changes after the coastal planning

曹祖德[6]等利用波浪槽試驗結果, 推導出不同流態(tài)、泥沙粒徑時的臨界起動力τb:

式中:ν為水體動粘滯系數;ρs為泥沙顆粒密度;ρ為水的密度;d為沉積物中值粒徑。

本文利用招遠岸段沉積物調查資料、908渤海灣沉積物調查資料、海灣志第3冊龍口港的資料, 得到研究區(qū)的中值粒徑分布。利用surfer8.0, 對數據進行網格化, 每隔120 m取一個點, 得到研究區(qū)的中值粒徑, 利用上述公式, 計算出萊州灣東岸附近海域沉積物臨界啟動力。

4.2 極限沖刷深度預測結果

4.2.1 N向風

規(guī)劃用海實施后, 研究區(qū)無沖刷區(qū)域變大, 屺坶島前端沖刷強度不大, 龍口灣到填海工程之間的區(qū)域沖刷作用不明顯。波浪和潮流在傳播的過程中受到填海工程的阻擋, 發(fā)生繞射和反射, 工程前端沖刷強度增大。春雨漁港西側1 000 m左右至石虎咀的近岸海域, 最大沖刷深度0～4 m, 界河口西側也有輕微的沖刷。

4.2.2 NE向風

規(guī)劃用海實施后, 研究區(qū)無沖刷區(qū)域變大, 屺坶島和龍口灣最大沖刷深度無較大變化, 但是界河到石虎咀近岸海域最大沖刷深度變小, 填海工程西側最大沖刷深度變大。

兩種工況下的計算結果見圖4和表1。通過經驗公式計算最大沖刷深度發(fā)現, 受到填海工程的影響,龍口灣流速和有效波高變小, 不利于泥沙的啟動和運移, 最大沖刷深度變化不大; 填海工程前端由于挑流和破浪作用, 潮流和波浪作用較強, 有利于泥沙的啟動和運移, 最大沖刷深度增大; 界河至石虎咀近岸海域的波浪和潮流受工程影響而變小, 最大沖刷深度較填海前變小。所得結論與數值模擬得到的沖淤趨勢一致。

5 結論

表1 工程前后八級N和NE向風最大沖刷深度對比表Tab.1 The comparison between the maximum scour depth before and after the engineering N and NE direction 8 degree wind

本文運用Mike21模擬研究區(qū)的潮流特征, 分析規(guī)劃填海前后研究區(qū)海洋環(huán)境的變化趨勢; 利用經驗公式計算出極端條件下的最大沖刷深度, 分析規(guī)劃填海實施后研究區(qū)的沖淤演變規(guī)律。研究成果主要有: (1) 工程前后, 研究區(qū)西北部流速大小基本不變; 受到工程影響, 龍口灣流速減小4～8 cm/s, 填海工程招遠部分至石虎咀的近岸區(qū)域流速減小 4～12 cm/s,工程形成的半圓形區(qū)域, 流速減小 24 cm/s; 填海工程西側由于挑流作用流速變大0～8 cm/s, 填海工程招遠部分與岸之間的水道以及兩工程之間的水道,流速增大 4 cm/s。(2) 通過經驗公式計算最大沖刷深度發(fā)現, 受到填海工程的影響, 龍口灣最大沖刷深度變化不大; 填海工程前端最大沖刷深度增大;界河至石虎咀近岸海域最大沖刷深度較填海前變小,與數值模擬得到的流場變化趨勢基本一致。

圖4 八級N, NE向風作用下的最大沖刷深度Fig. 4 The maximum scouring depth after and before the engineering

[1] 莊振業(yè), 陳衛(wèi)民, 許衛(wèi)東, 等. 山東半島若干平直砂岸近期強烈蝕退及后果[J]. 青島海洋大學學報, 1989, 19 (1): 90-98.

[2] 董衛(wèi)衛(wèi). 萊州灣東岸沖淤演變分析與防護[D]. 青島:中國海洋大學, 2008.

[3] 杜國云, 孫祝友. 渤海灣萊州灣東岸海岸帶人為地質災害研究[J]. 地質災害與環(huán)境保護, 2005, 16(03): 225-230.

[4] 中國海灣志編纂編委會. 中國海灣志: 第三分冊-山

東半島北部和東部海灣[M]. 北京: 海洋出版社, 1991: 1-79.

[5] Janssen P, Lionello P, Zambresky L. On the interaction of wind and waves [J]. Philosophical transactions of the royal society of London. Series A, Mathematical and Physical Sciences, 1989, 329(1604): 289-301.

[6] 曹祖德, 孔令雙, 焦桂英. 波、流共同作用下的泥沙起動[J]. 海洋學報, 2003, 3: 113-119.

[7] Silvester R. Coastal Engineering 2: Sedimentation, Estuaries, Tides, Effluents, and Modelling [M]. Netherlands: Elsevier Scientific Publishing Company, 1974: 72-77.

[8] Hsu J R C, Silvester R, Xia Y M. Static equilibrium bays: new relationships [J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1989, 115(3): 285-298.

(本文編輯: 劉珊珊 李曉燕)

Prediction of erosion evolution and deposition in the east coast of the Laizhou Bay after the implemention of the coastal planning

ZHOU Guang-zhen, FENG Xiu-li, LIU Jie, LIU Xiao, XU Fang
(College of marine geosciences, Ocean University of China, Qingdao 266100, China)

May, 17, 2012

east coast of the Laizhou Bay; hydrodynamic evolution; the limit depth of scouring

In this paper, the hydrodynamic model of Mike21 was used to simulate the flow current, deposition characteristics of the east coast of the Laizhou Bay, to predict its erosion and deposition evolution, Coupling the wave and current, to calculate the limit depth of scouring using empirical formula. By comparing the limit depth of scouring before and after planning reclamation, we found that the Longkou Bay would have a slight deposition. The wave and current of coastal water between the Jiehe River and Shihu area would become smaller and erosion intensity would become weaker, under the influence of engineering. All these analyses have shown that the coastal water between the Jiehe River and Shihu area would change from erosion to deposition when the planning reclamation is implemented.

P751

A

1000-3096(2014)01-0015-05

10.11759/hykx20120517002

2012-05-17;

2013-11-02

海洋公益科研專項(201005009)

周廣鎮(zhèn)(1986-), 男, 山東棗莊人, 碩士在讀, 主要從事海洋地質環(huán)境與工程研究, E-mail: anchor.ocean@163.com; 馮秀麗(1962-), 通信作者, 電話: 0532-66782057, E-mail: fengxiuli@ouc.edu.cn