龙口湾水动力特征及其对人工岛群建设的响应

2022-02-18费成鹏胡日军雒敏义姜胜辉张晓东朱龙海刘波

海洋地质与第四纪地质 2022年1期

费成鹏，胡日军,2，雒敏义，姜胜辉,2，张晓东,2，朱龙海,2，刘波

1. 中国海洋大学海洋地球科学学院，青岛 266100

2. 中国海洋大学海底科学与探测技术教育部重点实验室，青岛 266100

3. 青岛市即墨区自然资源局，青岛 266200

海湾具有独特的自然环境和优越的地理位置，对人类生存及社会发展具有极其重要的意义[1]，然而海湾是一个相对封闭的海洋环境，其水交换能力较差，生态环境较为脆弱，对人工构筑物建设等人类活动的干扰较为敏感[2]。人工构筑物的建设可以有效增加自然岸线、提高海洋资源利用效率，但人工构筑物的建设会改变原有海域的地形条件，直接或间接地影响海洋生态环境[3]，导致所在海域潮流、潮位、波浪等水动力条件发生变化[4-5]，甚至经常引起更大程度的极端事件以及海平面上升现象[6]。因此，为了保护海湾环境实现可持续发展，研究人工岛建设等人类活动对海湾环境的影响尤为重要[1,7]。

围填海工程建设等人类活动对水动力环境的影响，一直以来吸引着国内外学者的广泛关注。近岸潮汐的变化与围填海工程的建设密切相关[8]，Byun 等[9]认为在潮汐占主导地位、壅塞的河口或港湾环境中建造堤坝或海堤等围填海工程，可能导致潮汐振幅显著增加和水流速度下降。围填海工程的建设不会导致潮流特征大范围改变，仅在工程区域附近有较大变化[10]，流向受工程影响会发生偏转，工程防波堤等会引起挑流作用的区域流速增大，其余区域流速普遍减小[11-12]。法国兰斯河口世界第二大潮汐电站[13]以及加拿大芬迪湾潮汐动力潟湖工程[14]对附近潮流也产生了类似的影响。海岸的形态对于波浪能的耗散具有关键作用[15]，围填海工程造成的海岸形态的改变往往会中断与波浪条件相关联的沿岸泥沙输运，同时会破坏海岸抵抗风暴的能力[16]，但设计良好的防波堤会有效减小入港波浪的有效波高，从而达到保护港口的目的[17-18]。围填海造成的海域面积减小会直接导致纳潮量减少[19-20]，例如马尔代夫国际机场扩建工程[21]、雅加达湾巨型海堤[22]以及胶州湾围填海工程[23]等。水交换也会受到围填海工程的影响[24]，岸线大规模的改变会导致水交换能力减弱，通过修复岸线对水交换能力具有促进作用[25]，岸线变迁引起的水动力环境的改变是影响水交换率的主要因素[26]。此外，还有学者在水动力环境变化的基础上探讨了冲淤环境以及地貌演变对海岸工程建设的响应[27-29]，同时对于围填海工程造成的水质污染[30-31]、海洋生物受损[32-33]等方面的研究也有很多学者涉猎。然而，前人的研究多集中于河口或者较为开阔的海湾或海域中围填海工程的建设，对于小海湾中大规模离岸人工岛群建设对海洋环境影响的相关研究却鲜有报道，因此，本文针对龙口湾水动力特征及其对人工岛群建设的响应进行较为全面的研究。

龙口湾面积约84.13 km2，龙口人工岛面积约35.23 km2，约占龙口湾面积的41%，因此，在龙口湾这类小海湾中建设大规模的围填海工程，势必对其海洋生态环境造成明显影响。刘星池[34]等通过建立龙口湾水沙数值模型来预测人工岛建设对海洋环境的影响，人工岛建设会导致龙口湾潮流流速流向均发生不同程度的变化，冲淤特征表现为人工岛北侧和西南侧海域以冲淤为主，西侧海域以侵蚀为主[35]，人工岛外悬浮泥沙浓度和悬沙通量大于岛内水道[36]，同时表层沉积物组合特征也产生了重要变化[37]。但是，目前对于人工岛群建设对龙口湾水环境的影响研究主要集中在潮流场变化，缺乏其对水动力环境的系统研究；同时，龙口人工岛的实际建设方案与设计方案存在明显差异，导致前人对于该区域的研究存在误差，不能完全反映人工岛群建设对龙口湾水动力特征的影响。因此，本文利用龙口湾以及附近海域的水深地形、潮汐潮流等实测资料，利用Mike21 数学模型，分析了人工岛建设前后龙口湾的水动力特征，并在此基础上探讨了人工岛群建设对龙口湾潮流、波浪、潮位、纳潮量以及水交换的影响。

1 研究区概况

龙口湾为莱州湾的一个附属海湾，湾廓呈半圆形，是典型的连岛坝成因的次生海湾[38]。20 世纪90 年代以来龙口湾海岸开发活动不断增多，极大地改变了龙口湾的海岸形态和地形地貌[39]，其中以龙口人工岛规模最大。除航道外，湾内水深不足10 m，湾外水深一般为10～20 m（图1）。屺坶岛以南的连岛海岸以基岩海岸或人工岸线为主，人工岛群以南主要为砂质海岸。表层沉积物主要以粉砂和砂为主，湾外沉积物粒径较粗[40]。潮汐性质为不规则半日潮，潮流性质以不规则半日潮流为主，潮流运动形式主要为往复流[41]。波浪以风浪为主，湾内常浪向为SW 向，强浪向为WSW 向，湾外常浪和强浪向均为NE 向[42]。

图1 研究区地理位置Fig.1 Location of the study area

2 研究方法

2.1 潮流场数值模拟

本文潮流场模拟采用Mike21 中的Flow Model FM HD 模块[43]进行，模拟采用非结构三角网格剖分计算域，计算域范围如图2 所示，坐标范围为36°59′15.743″～40°59′21.417″N、117°32′22.881″～122°39′30.992″E，覆盖整个渤海海域及部分北黄海海域[44]。人工岛群建设前后除工程区域外其他位置的网格一致，增强了前后对比的准确性和可靠性。

图2 研究区计算网格图Fig.2 Gridding of the study area

对位于潮滩区的干、湿网格采用动边界的方法进行处理。为能清楚了解研究区所在海域潮流场特征，对该海域网格作局部加密处理。水深资料采用中国人民解放军海军航海保证部制作的海图1∶100 万（10011 号）、1∶15 万（11370 号、11570 号、11710 号、11840 号、11910 号、11932 号）水深资料及研究区附近海域最新实测水深资料。

2.2 波浪场数值模拟

波浪场模拟采用Mike21 中的SW 浅水波浪数值模块，该模型广泛适用于大范围或者局部区域的波浪预报和分析以及不同历史条件下近海、海岸和港口结构物设计过程中的波浪情况预报[45]。

波浪数值模拟的计算域及网格、岸界和水深资料与本文潮流场数值模拟设置一致，分别模拟了研究区N 向和SW 向六级风（12 m/s）作用24 小时下的波浪场状况，以此来探讨人工岛群建设对龙口湾波浪场的影响。

2.3 纳潮量计算

纳潮量是一个水域可以接纳潮水的体积，海湾的纳潮量不仅是衡量海湾开发价值的一个水文指标，而且也是反映湾内外水交换的一个重要参数[46]，本文对纳潮量的计算采用叶海桃[47]等对三沙湾纳潮量的算法，其公式如下：

式中：P为纳潮量，单位：m3；

ΔH为潮差，单位：m；

A0为最低潮位下水域面积，单位：m2；

ΔHi′为潮滩上第i个网格高潮位时的水深，单位：m；

Ai为第i个网格上最低潮位时的水域面积，单位：m2。

由于龙口人工岛最西侧填海区域超出了龙口湾界线（37°32′15.982″～37°40′12.993″N、120°13′26.819″～120°13′43.580″E），无法采用原龙口湾的界线来计算纳潮量，因此本文选取从屺坶岛至石虎咀断面和海岸线封闭的区域为本文计算区域（图1），以下简称“计算区域”。

2.4 水交换计算

水交换是海洋环境科学研究的一个基本命题，水交换率的计算是研究海湾自净能力的重要指标和手段[48]。本文对水交换率的计算采用保守污染物浓度扩散的方法[49]，在Mike21 水动力数值模型中，初始时刻将计算区域内示踪剂浓度设为1，外海域设为0，其他水动力条件保持不变，得到水交换率等于（1-浓度值）×100%。水交换计算区域与纳潮量计算区域一致，计算时长30 d。

3 结果

3.1 模型验证

本文收集了大连、旅顺、鲅鱼圈、曹妃甸、大口河、潍坊港、北隍城、八角、烟台港、龙口港等10 个潮位站历史观测资料并采用傅氏分析方法[50]进行调和分析，选用M2、S2、K1、O1四个分潮的调和常数预报出大潮期的潮位变化，同时结合中国海洋大学2017 年5 月10—11 日在工程附近进行的2 个站位的潮位观测资料进行潮位验证；采用2017 年5 月10— 11 日（大潮）6 个站位27 小时海流同步连续观测资料进行流向流速验证，验证结果表明，模拟结果和实测值吻合较好。限于篇幅限制，本文仅列出研究区周边的潮位验证曲线（C4 站位和龙口港）和潮流验证曲线（C4、C5 和C6 站位），见图3 和图4。

图3 潮位验证曲线Fig.3 Tide verification curve

图4 流速流向验证曲线Fig.4 Current velocity and direction validation curve

3.2 潮流场特征

3.2.1 人工岛建设前潮流场

潮流场数值模拟结果表明，在涨急时刻，湾外潮流流向整体上为西南向，受地形影响较小，流向变化较小，湾内潮流方向主要为西南向和南向，在屺坶岛西南侧潮流转为东南向，进入湾内则变为东向，向南逐渐转为西南向；流速在屺坶岛西侧最大，可达0.68 m/s，龙口湾内流速较小，整体小于0.2 m/s，且越靠近岸线流速越小。落急时刻的整体规律与涨急时刻相反，潮流带为明显的往复流形式，湾外潮流流向整体上为东北向，湾内潮流流向在龙口湾内受地形影响发生逆时针旋转，转为西北向至屺坶岛流向湾外；流速在屺坶岛西北侧最大，可达0.8 m/s，龙口湾内流速分布与涨急时刻类似（图5）。

图5 人工岛建设前涨急和落急时刻潮流场分布图（大潮期）Fig.5 Current fields during flood and ebb of tide before construction of artificial island (spring tide)

3.2.2 人工岛建设后潮流场

人工岛建设后，龙口湾潮流场特征发生了明显的改变。涨急时刻，湾外潮流流向整体上为西南向，流向变化较小，湾内潮流流向较为复杂，整体上围绕人工岛外侧流出湾内，小部分进入岛内水道，在人工岛防波堤内存在一个漩涡；湾外流速较大，为0.2～0.8 m/s，流速由北向南逐渐减小，湾内流速整体小于0.1 m/s。落急时刻潮流流向整体上与涨急时刻相反，湾外潮流流向整体上为东北向，湾内潮流在人工岛南部逆时针转为西北向，在人工岛防波堤处一部分向北流出湾内，一部分转为东方向进入龙口港围绕岸线流向湾外；湾外流速为0.2～0.8 m/s，湾内流速整体上小于0.1 m/s（图6）。

图6 人工岛建设后涨急和落急时刻潮流场分布图（大潮期）Fig.6 Current fields during flood and ebb of tide after construction of artificial island (spring tide)

3.3 波浪场特征

3.3.1 SW 向6 级风作用下波浪场特征

人工岛建设前，SW 向6 级风（风速12 m/s）作用12 小时情况下，湾外有效波高较大，大部分区域有效波高大于1.2 m，受SW 向风和地形的影响，波高等值线形态向南凸出。湾内有效波高整体小于1.2 m，近岸人工构筑物的掩蔽作用[51]较为明显，在龙口港内以及屺坶岛南侧部分半封闭区域形成波影区，有效波高小于0.4 m；受湾内水深影响，波高等值线变化率较快，且逐渐趋于与岸线平行（图7）。

图7 SW 向6 级风下人工岛建设前后波浪场分布图Fig.7 Wave fields before and after construction of artificial island under southwest wind

人工岛建设后，SW 向6 级风作用12 小时情况下，湾外的波浪特征并未发生明显的变化，有效波高整体上大于1.2 m，波高等值线形态向南凸出。湾内有效波高整体减小，屺坶岛南侧局部区域有效波高能达到1.3 m 左右，其他区域有效波高小于1.1 m。人工岛的建设使得龙口湾的地形变得更为复杂，除龙口港内和屺坶岛附近形成波影区外，人工岛水道内也形成波影区，波影区有效波高小于0.4 m。

3.3.2 N 向6 级风作用下波浪场特征

人工岛建设前，N 向6 级风（12 m/s）作用12 小时情况下，湾外有效波高整体较大，大部分区域有效波高大于1.65 m，受N 向风和地形的影响，波高等值线形态向南凸出。湾内有效波高整体小于1.2 m，波高等值线形态向东凸出，受湾内水深影响，波高等值线变化率较快，且逐渐趋于与岸线平行（图8）。

图8 N 向6 级风下人工岛建设前后波浪场分布图Fig.8 Wave fields before and after construction of artificial island under north wind

人工岛建设后，N 向6 级风作用12 小时情况下，湾外的波浪特征并未发生明显的变化，有效波高整体上大于1.65 m。湾内有效波高整体减小，在人工岛防波堤处波高骤减，形成波影区，湾内有效波高整体小于1.35 m，人工岛建设区域有效波高形态变化较大，人工岛北侧区域有效波高形态变化较小。

3.4 纳潮量

本文分别计算了人工岛建设前后大潮期和小潮期的纳潮量，结果如表1 所示，人工岛建设前大潮期纳潮量1.3620×108m3，小潮期纳潮量9.1227×107m3，平均纳潮量1.1371×108m3；人工岛建设后大潮期纳潮量1.1749×108m3，小潮期纳潮量7.8660×107m3，平均纳潮量9.8075×107m3。相比工程建设前后，大潮期纳潮量减小了13.74%，小潮期纳潮量减小了13.78%，平均减小了13.75%，可见人工岛的建设对计算区域纳潮量的影响很大。

表1 人工岛建设前后纳潮量Table 1 Tidal prism before and after construction of artificial island

3.5 水交换

根据计算区域30 d 的水交换率计算结果（图9），人工岛建设前，30 d 平均水交换率为62.58%，在计算边界附近水交换能力最强，水交换率为80%～90%，水交换率由边界处向湾内依次递减，越靠近岸线减小速率越快，在龙口港附近最弱，水交换率小于20%。人工岛建设后，30 d 平均水交换率59.82%，在计算边界附近水交换率在80%左右，在龙口港附近及人工岛内部水域水交换率迅速减小，龙口湾内水交换率小于20%，局部区域小于5%，人工岛内部水道水交换率为0～80%，南部和西部区域可达80%，中间水道局部区域水交换率小于5%。

图9 人工岛建设前后30 d 水交换率分布图Fig.9 30 days water exchange rate before and after construction of artificial island

4 讨论

4.1 人工岛的建设对潮流场的影响

（1）潮流流向变化

人工构筑物的建设会引起潮流流向的偏转[52]，本文中将人工岛建设前后的流场进行了叠加并将变化明显的区域进行了局部放大（图10、图11）。可见人工岛建设后，在涨急时刻人工岛西北侧（A区域）流向顺时针偏转，人工岛北侧区域（B 区域）由原来的开放海域变成半封闭海域，潮流流向在此变化较为复杂，整体上顺时针偏转，人工岛西侧（C 区域）由于人工岛防波堤的挑流作用形成了一个逆时针旋转的漩涡，人工岛南侧（D 区域）潮流流向逆时针偏转（图10）。在落急时刻人工岛西北侧（A 区域）流向顺时针偏转，人工岛北侧（B 区域）流向整体上逆时针偏转，人工岛西侧（C 区域）和人工岛南侧（D 区域）流向逆时针偏转（图11）。

图10 人工岛建设前后涨急时刻流向对比图（大潮期）Fig.10 Current direction during flood of tide before and after construction of artificial island (spring tide)

图11 人工岛建设前后落急时刻流向对比图（大潮期）Fig.11 Current direction during ebb of tide before and after construction of artificial island (spring tide)

（2）潮流流速变化

本文分别选取了工程建设前后大潮期涨急时刻和落急时刻的潮流场进行了流速对比（图12）。通过对比分析，龙口人工岛建设前后，研究区的潮流场特征发生了较为明显的变化，尤其是在人工岛附近区域变化更为明显，具体表现在：人工岛西部靠近防波堤的区域由于挑流作用[53]流速增大，最大流速变化超过0.4 m/s，变化率可达60%，涨急时刻流速变化相比落急时刻变化较大，距离防波堤越远，流速变化越小。李池鸿[54]等基于Mike21 模型对新建的洋浦神头港区南港区码头工程进行了研究，结果表明工程后防波堤口门附近涨、落潮都会因挑流的影响局部流速会变大。人工岛的建设阻挡了部分区域原有潮流的流动，致使潮流流速减小，在人工岛北部，流速整体减小，涨急时刻减小范围大多为0.08～0.16 m/s（变化率30%～70%），极少数区域超过0.16 m/s；落急时刻流速变化较大，流速减小超过0.16 m/s（变化率约70%）的区域大幅增加；在人工岛内部水道以及人工岛南部区域，流速减小，涨落急时刻变化较为相似，流速减小范围整体上为0.08～0.16 m/s（变化率30%～70%），在界河口的一小部分区域流速减小范围超过0.16 m/s。潮流流速的变化可能是引起人工岛北侧、西南侧和内部水道产生淤积以及人工岛西侧发生冲刷[35]的主要原因。

图12 人工岛建设前后涨急和落急时刻流速对比图（大潮期）Fig.12 Current velocity during flood and ebb of tide before and after construction of artificial island (spring tide)

（3）余流场变化

为了分析人工岛建设对龙口湾余流场的影响，本文选取了潮流数模结果中一个完整大小潮周期（15 天）的流速流向计算了人工岛建设前后研究区的欧拉余流（图13），并进行了差值对比（图14）。人工岛建设前，龙口湾湾外除屺坶岛附近区域外余流流向整体上为北向，流速普遍在0.05 m/s 左右，在屺坶岛西北侧和西南侧分别形成一个顺时针漩涡和逆时针漩涡，流速为0.05～0.15 m/s，龙口湾西侧存在一个顺时针漩涡，其流速在0.05 m/s 左右；湾内余流流向较为复杂，余流从龙口湾西北侧进入后分为两支，一支整体呈逆时针运移在龙口湾北侧流出，另一支呈顺时针运移在龙口湾南侧流出。

图13 人工岛建设前后余流特征Fig.13 Residual current characteristics before and after the construction of artificial island

图14 人工岛建设前后余流流速对比图Fig.14 Comparison of residual flow velocity before and after the construction of artificial island

人工岛建设后，湾外余流场未发生明显变化，湾内余流场变化较为明显，具体表现在人工岛西侧形成一个逆时针旋转的新漩涡，且流速增大，增大范围为0.03～0.4 m/s，人工岛北部余流整体上分为两支，一支逆时针运移在屺坶岛附近流出龙口湾，流速较人工岛建设前减小，另一支顺时针运移在人工岛防波堤处流出，流速略有增大，人工岛西南侧区域余流整体上向西南方向运移，流速较人工岛建设前减小。

4.2 人工岛建设对波浪场的影响

海域波浪场的分布特征与水深、地形和风速风向等要素密切相关，在岸线走向、海底地形、风速风向的影响下，波浪在传播过程中会发生一定的折减、绕射等衰减现象[55-56]。本文为了更好地分析人工岛的建设对波浪场的影响，将人工岛建设前后的有效波高做了差值对比（图15），在SW 向6 级风作用下，湾外有效波高未发生明显变化；湾内变化主要集中在人工岛北部及内部水道区域，其中人工岛内部区域有效波高变化最为明显，人工岛对该区域的掩蔽作用较强，导致有效波高在此明显减小，有效波高减小范围为0.3～1.2 m（变化率24%～96%）；人工岛北部海域人工岛的掩蔽作用较小，该区域有效波高减小，减小范围为0.15～0.6 m（变换率12%～48%）。在N 向6 级风作用下，湾外有效波高未发生明显变化，湾内变化主要集中在人工岛水道及人工岛西南侧区域，人工岛的掩蔽作用导致该区域的有效波高相对人工岛建设前整体减小，在人工岛防波堤处有效波高减小可达1.2 m，变化率可达96%。

图15 人工岛建设前后波浪场对比图Fig.15 Comparison of wave field before and after the construction of artificial island

波浪对比结果表明，龙口人工岛的建设导致工程附近区域的波浪有效波高整体减小，其影响范围主要集中在人工岛周边。顾杰[57]等研究结果也表明工程实施前后的波高变化仅集中在工程区域附近，工程的建设使该区域波高显著减小。

4.3 人工岛的建设对潮位的影响

围填海的建设会导致潮位产生一定的变化[58]，为此，我们在人工岛北部、西部、南部以及人工岛内部共选取了12 个代表点（图1）分别计算了大潮期一个潮周期内的最大潮差变化，计算结果如表2所示。通过对比分析发现，在人工岛北部潮差减小，最大潮差变化在-0.017 m 左右，这主要是由于人工岛的建设使得该区域的潮汐动力减弱，致使潮差变小[59]；在人工岛西侧潮差变化较小，呈现出离人工岛越远变化越小的趋势；人工岛南部潮差增大，距人工岛1000 m 的位置潮差增大约0.014 m，远离人工岛变化逐渐变小；人工岛内部水道由于壅水作用[60]潮差变化最为明显，最大潮差变化可达0.047 m。胶州湾跨海大桥的建设对胶州湾潮差的影响普遍小于0.01 m[61]，相比龙口人工岛建设引起的潮差变化较小，其主要原因是龙口人工岛占据的过水面积比例更大，对水流特征的影响更明显。

表2 人工岛建设前后代表点潮位变化Table 2 Tide changes before and after construction ofartificial island (spring tide)

4.4 人工岛的建设对纳潮量的影响

人工岛群的建设直接占用了计算区域的海域面积，改变了原有水动力环境，龙口人工岛建设后，本文计算区域的海域面积减少20.68%，大潮期纳潮量减少13.74%，小潮期纳潮量减少13.78%，平均纳潮量减少13.75%。本文整理了前人对于不同海湾纳潮量的研究成果（表3），通过对比发现，莱州湾和罗源湾围填海造成的海域面积减少对纳潮量的影响相对较小[62-63]（纳潮量变化率小于海域面积变化率），锦州湾、芝罘湾、湛江湾以及象山港海域面积减小引起的纳潮量变化较大[64-67]（纳潮量变化率大于海域面积变化率），可见，本文计算区域纳潮量的变化与莱州湾和罗源湾较为相似。总体而言，人工岛群的建设是造成计算区域纳潮量减少的主要原因，这与诸多海湾受围填海影响导致纳潮量减小的结果是相符的。

表3 不同海湾纳潮量变化对比Table 3 Variation of tide prism in different bays

4.5 人工岛的建设对水交换的影响

人工岛群建设后，本文研究区域的地形岸线变得更为复杂，对原有的水交换能力产生了一定的影响，尤其是在人工岛群附近及其内部水道，水交换率有了明显的改变，我们将人工岛建设后的30 d 水交换率做了差值对比（图16），在人工岛建设后，计算区域30 d 平均水交换率减小2.76%，具体表现在人工岛北部靠近屺坶岛的区域和人工岛南部界河附近区域水交换率明显变大，在人工岛内部水道及北部附近区域水交换率明显减小。本文收集了其他海湾水交换率的变化情况[62,66-68]（表4），通过对比发现，水交换率的变化与海域面积的变化并无明显相关关系，尤其是湛江湾和锦州湾，在减少更多面积的情况下，平均水交换率变化幅度反而减小。因此，围填海造成的海域面积减小对水交换率的直接影响不大，其引起的水动力环境的改变是影响水交换率变化的主要原因[66]。