林纪江，胡日军，王平，宋卓利，张文静，牛建伟，伊善堂，丁小雨

1.国家海洋局南海信息中心，广州 510310

2.自然资源部海洋环境探测技术与应用重点实验室，广州 510310

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

4.国家海洋局北海海洋工程勘察研究院，青岛 266061

5.中国科学院南海海洋研究所，广州 510301

6.广州海洋地质调查局，广州 510075

表层沉积物再悬浮是影响水体中悬浮泥沙浓度的一个重要因素，悬浮泥沙浓度的变化过程及运移趋势指示了不同的水动力环境，对海洋沉积、生态环境、岸线冲淤以及物质循环具有重要意义。前人对庙岛海峡周边海域开展了大量的研究工作，主要集中于大尺度的渤黄海悬浮体季节性分布、变化及运移和沉积动力、地质地貌等方面，认为渤海悬浮泥沙浓度分布的主要影响因素是水动力条件（海浪、海流等），悬浮泥沙输运主要动力形式包括平流输运、潮泵效应等[1-2]，底质类型对水体含沙量影响较大，细颗粒物质更易发生再悬浮[3-4]。庙岛群岛周边沉积物分布及沉积动力格局主要受潮流控制[5]，冬季是渤黄海海域悬浮体输运的主要季节，山东半岛沿岸在强动力下的再悬浮是水体悬浮体的主要来源[6-8]。对庙岛海峡周边小尺度海域悬浮泥沙输运机制与趋势的研究较少。研究区受到庙岛群岛和庙岛海峡的双重影响，其悬浮泥沙的分布和输运必然具有独特的性质，本文基于连续实测水文泥沙资料，利用不同数理公式和物质通量分解方法，分析探讨庙岛海峡周边悬浮泥沙分布、沉降及表层沉积物再悬浮特性，揭示悬浮泥沙输运的动力机制和运移特性，不仅可以丰富庙岛海峡海域沉积动力过程的相关研究内容，加深对黄渤海泥沙交换过程的理解，而且可以为周边海域生态修复、环境保护等提供科学参考。

1 区域概况

庙岛海峡位于渤海海峡的南部，北侧是庙岛群岛，南侧是山东半岛，是渤黄海水沙交换的主要通道。研究区海域属于海洋暖湿季风性气候，年平均气温13℃，累年平均降水量为586.7 mm，多年平均风速4.2 m/s，常风向为SSW向[9]，最大水深约37 m[10]。潮汐类型属于不规则半日潮，平均潮差1 m左右，最大潮差2.9 m[11]；潮流属于不规则半日潮，以往复流为主，落潮流向西，涨潮流向东，落潮平均历时略小于涨潮平均历时；波浪以风浪为主，涌浪次之，强浪向为NE向，最大波高4.1 m，常浪向为NNE向，次常浪向为ENE向[12]。庙岛海峡周边海底有登州浅滩和登州水道两个重要的地貌元素。登州浅滩位于蓬莱的西北侧，主要由4个沙洲组成，地形复杂，大致位于10 m等深线以外，登州水道位于登州浅滩的东北侧，是渤黄海水沙交换的主要通道，蓬莱西海岸近岸海底与登州浅滩之间有一条EN—WS向的海底洼地，比两侧区域深约1～2 m，坡度较缓，并无明显的边缘陡坎，近岸海底区域在2 m等深线附近一般发育有水下沙坝[9]。表层沉积物类型多样，以砂质粉砂分布最广，其次为粉砂，砾与砂质砾集中分布在登州水道内[5]。

2 资料与方法

2.1 资料来源与处理

2.1.1 海流

2011年7 月天津水运工程勘察设计院在庙岛海峡及其两侧海域布设9个站位（图1），进行大、小潮期的27 h海流连续观测，2015年10月中国海洋大学在庙岛群岛海域布设3个站位进行大潮期的27 h海流连续观测，采集流速、流向、水深等，测流间隔为1 h，单次测流历时2.5 min，2011年数据按“六点法”进行处理（即表层、0.2H、0.4H、0.6H、0.8H和底层（H为取样点的水深）），2015年数据按表层、中层、底层进行处理，对海流资料进行准调和分析，得出研究区的海流性质、余流的特征。

图1 研究区水深地形及调查站位Fig.1 Topography of the study area and survey stations

2.1.2 悬浮泥沙

悬浮泥沙采用CTD（日本）与测流同步进行测定，为保证含沙量数据的准确性，在测量过程中使用竖式采水器同步采取现场水样，经实验室处理（抽滤、烘干、称重）后，用其率定同步获取的浊度数据。2011年数据按“六点法”进行处理（即表层、0.2H、0.4H、0.6H、0.8H和底层（H为取样点的水深））。2015年数据按表层、中层、底层进行处理。悬浮泥沙粒度的测定采用Mastersizer2000型激光粒度仪，测量范围为0.02～2000 μm，重复试验的相对误差＜3%，粒级分辨率为0.1Φ，粒径间距为0.25 Φ。由于滤膜中粘附的悬浮泥沙样品不等量，仅对9个站位中粘附量大的部分滤膜进行测定，测定步骤同常规激光粒度仪粒径测定。

2.1.3 表层沉积物

2011年6—7 月天津水运工程勘察设计院在庙岛海峡周边海域通过蚌式取样器获取了107个表层沉积物样品，沉积物粒度测定采用NSY-Ⅲ型宽域粒度分析仪（河海大学）和SFY-B2000音波振动式半自动筛分粒度仪（中科院南京地理与湖泊研究所）。

2.2 研究方法

2.2.1 悬浮泥沙沉降速度

悬浮泥沙沉降速度计算公式参考Camenen[13]的半经验公式，通过A、B和m3个参数将泥沙颗粒的形状因子（csf）和圆度（P）引入泥沙沉降公式中，公式如下：

式（1）中，Ws－沉降速度；Re－雷诺数；v－流体动粘滞系数；d－颗粒直径；A，B，m－阻力系数方程中的系数；d*－无量纲的粒子直径。

2.2.2 泥沙起动流速

泥沙起动流速计算采用窦国仁1999年提出的泥沙起动流速公式[14]，各参数及含义见文献[14]。

2.2.3 再悬浮通量与沉降通量

再悬浮通量采用Partheniades提出的公式[15]:

沉降通量采用Krone提出的公式[16]:

式（5）中E是再悬浮常数，其范围一般为5×10-5～2×10-4[17]，再悬浮常数对再悬浮通量的计算影响较大，本文选取再悬浮常数的最小和最大值分别进行估算； τ0为底切应力； τcr为临界切应力；Cb为近底部悬浮泥沙浓度； ωs为悬浮泥沙的沉降速度。

2.2.4 悬浮泥沙单宽通量

悬浮泥沙单宽通量机制分解依据 Ingram[18]和Uncles 等[19]提出的用相对水深进行分解瞬时的物质输运量。潮周期Tt平均瞬时单宽悬沙输移通量T的计算式为：

式（7）中：T1－欧拉余流贡献项；T2－Stokes漂移输移量；T1+T2－平流输移；T3－潮汐与悬沙潮变化相关项；T4－悬沙与潮流变化相关项；理论上，，，ht各自潮周期平均值为零，但三者的相关产生了泥沙净输移，即潮汐捕集作用；T5－垂向流速变化和含沙量变化的相关项，为垂向净环流的贡献；T6，T7－时均量和潮汐振动切变引起的剪切扩散；T8－垂向潮振荡切变作用。

其中，垂线平均含沙量采用黄才安和梅小文[20]推荐的流量法计算：

式（8）中，垂线平均含沙量用S1表示，h代表垂线水深，s代表距底床y处的当地含沙量，u代表距底床y处的当地流量。

3 结果与讨论

3.1 海流与底质特征

研究区海域潮流性质以不规则半日潮流性质为主，实测潮流最大流速值为1.84 m/s，流向为284°，出现在大潮期A5站位的表层，潮流最小流速值为0.04 m/s，流向为164°，出现在小潮期A9站位的底层，垂线平均最大1.53 m/s，潮流流速值呈现自庙岛海峡至周边海域、自表层至底层递减的趋势。大、小潮期潮流流向基本一致（A8站位除外），A3站位潮流流向最为复杂，表现出旋转流的性质，A2、A5、B3站位受地形影响明显，潮流流向较单一，往复流特征显著（图2）。

图2 各站位垂向平均潮流矢量[22]Fig.2 The vector of vertical-averaged current of each station

各站各层余流流速相差较小，大、小潮期最大值分别为26.4、16.2 cm/s，均出现在A3站位表层，A2、A5站位自表层至底层余流逐渐增大，其余各站位余流大致自表层至底层逐渐减小；大、小潮期垂线平均余流方向A1、A2、A8站位不一致，其余站位基本一致（图3）。基于Mike-21数值模拟软件模拟的研究区海域余流场见图4[21]，由图4可见，庙岛海峡南侧余流由渤海指向黄海方向，庙岛海峡北侧余流由黄海指向渤海方向，并在中部产生小环流。

图3 各站位垂向平均余流矢量Fig.3 Vector of vertical-averaged residual current at each station

图4 研究区余流场数值模拟结果[21]Fig.4 The numerical simulation results of residual current field in the study area

根据Folk分类命名，研究区表层沉积分为砾（G）、砂质砾（sG）、砾质砂（gS）、含砾砂（(g)S）、砂质粉砂（sZ）、粉砂（Z）、砂质泥（sM）和泥（M）8种类型（图5），含砾样品中砂的含量较高（其中4个样品为砾（G），直径最大3 cm，最小0.5 cm），不含砾样品中粉砂含量较高（图6）。

图5 研究区沉积物类型分布Fig.5 Distribution of sediment types in the study area

图6 各类型沉积物粒级频率曲线N为样品数，不同颜色曲线表示不同站位。Fig.6 Grain size frequency curve of various types of sedimentsN is the number of samples, different color curves represent different stations.

3.2 悬浮泥沙分布及变化特征

3.2.1 时空分布特征及变化

实测最大悬浮泥沙浓度为155 mg/L，出现在A4站位的底层，研究区海域垂向平均悬浮泥沙浓度整体上较低，最大值为37.04 mg/L（大潮期A4站位）；庙岛群岛北侧海域表、中、底层悬浮泥沙浓度差别不大（图7），群岛中部表、中层差别不大，远低于底层悬浮泥沙浓度，其他海域自表层至底层递增明显；落潮段悬浮泥沙高于涨潮段（A6、B3站位除外），庙岛海峡、庙岛群岛海域悬浮泥沙浓度高于周边开阔海域，呈现自岸向海悬浮泥沙浓度递减趋势，这可能是海峡处、近岸处底部泥沙发生再悬浮所致。

图7 悬浮泥沙涨落潮段(左图)与表中底各层(右图)平均浓度分布[22]Fig.7 Distribution of suspended sediment concentration averaged during ebb and flow period (left) and in each layer (right) [22]

悬浮泥沙浓度随时间变化呈现出周期性波动，9～16 h尺度的周期变化在整个分析时段表现的非常稳定，具有全域性[22]，通过对比各站位悬沙浓度时间序列剖面（图8）发现，潮周期内庙岛海峡悬浮泥沙浓度峰值形态较平滑，高值持续时间较长，其他区域悬浮泥沙浓度峰值曲线形态较尖锐，高值持续时间较短。

3.2.2 潮周期内悬沙浓度变化特征

根据相关经验，潮流流速的变化直接影响水体悬浮泥沙浓度的变化，当流速增大时，海底表层沉积物受到的剪切应力相应变大，海底受到冲刷，细颗粒的泥沙物质进入水体，悬浮泥沙浓度增高；当流速减小到一定程度时，海底表层沉积物受到的剪切应力小于临界剪切应力，水体中的悬浮泥沙发生沉降，悬浮泥沙浓度也相应降低。对比潮流流速与悬沙浓度时间序列剖面（图8）发现，悬浮泥沙浓度潮周期变化与潮流流速周期变化相关性较好，潮流流速超过40 cm/s时，研究区海域发生明显的再悬浮现象，并且潮流流速越大，底部悬浮泥沙浓度越高，悬浮泥沙浓度峰值整体上滞后于潮流流速峰值0.5～2 h。潮周期内各站位的悬浮泥沙浓度随流速出现高低值的更迭，且存在分层现象。悬沙浓度高低值更迭表现为底层变化对高流速响应明显，表层变化对低流速响应明显，底层悬浮泥沙浓度大于表层。分层现象可能是研究区海域存在跃层的影响：当流速较高时，底质沉积物再悬浮，底层浓度增加，因跃层影响，表层接受自底部向上扩散的悬沙很少；当流速处于低值时，表层的悬浮泥沙会发生落淤，使表层悬浮泥沙浓度明显降低，因此表层变化对低流速反应明显。

3.3 泥沙沉降及再悬浮特性

3.3.1 悬浮泥沙沉降特征

研究区悬浮泥沙沉降速度计算结果表明，单颗粒沉降速度整体上不大，最大为0.030 6 mm/s，出现在研究区西部的A9站位；最小值为0.024 6 mm/s，出现在研究区中部的A6站位（表1）。有学者研究结果表明，泥沙絮凝沉降速度是单颗粒悬浮泥沙沉降速度的几十倍至上百倍[23-26]，在悬浮泥沙絮凝沉降中起主导作用的有悬浮泥沙粒径、盐度、水体流速、含沙量4个因子[27]，悬浮泥沙粒径＜16 μm絮凝沉降现象明显，悬浮泥沙粒径8 μm左右，絮凝沉降现象最为强烈[28]，水体中细颗粒物质发生絮凝沉降的必要条件是阳离子浓度为60～500 mmol/L[29]，研究区悬浮泥沙粒径为10.4～11.6 μm、盐度值30.87～31.21（其对应的水溶液阳离子浓度在前述区间内），均适宜悬浮泥沙的絮凝沉降。有学者[26,30]基于控制变量的室内试验证明，基于同等其他条件，在较低悬浮泥沙浓度（0～1 000 mg/L）情况下，亦可发生絮凝沉降，絮凝平均沉降速率与悬浮泥沙浓度正相关，但此时悬浮泥沙浓度对沉降速率影响较小。在水体流速小于30 cm/s的情况下，粒径为8～32 μm的悬浮泥沙可以絮凝沉降[31]，由图9可见，A4、A5站位流速小于30 cm/s的持续时间较短，其余站位表层流速小于30 cm/s的持续时间较长，表明研究区除庙岛海峡外其他海域水动力条件适合悬浮泥沙絮凝沉降，并以中、底层絮凝沉降为主。

图9 单日内各站位表中底层流速小于30 cm/s的持续时间（h）Fig.9 Duration of flow velocity less than 30 cm / s in a single day in the each layer of each station (h)

表1 A1—A9各站悬浮泥沙沉降速度Table 1 The settling velocity of suspended sediment at stations A1-A9

3.3.2 沉积物再悬浮特征

表层沉积物起动流速计算结果表明（表2），研究区内自庙岛海峡向渤海、黄海、庙岛群岛海域表层沉积物起动流速递减，研究区东侧表层沉积物起动流速自岸向海递减，西侧自岸向海先增大后减小，对比各站位流速时间序列剖面图（图8），A2、A4、B3站位最大底层潮流流速均达不到表层沉积物起动流速。根据同期获取的底质资料发现，A2站位表层沉积物主要为泥，颗粒粒级细，颗粒组分中黏土含量较高，颗粒间黏性大，需要的起动流速较大；A4站位表层沉积物主要为砾，颗粒粒级粗，自身的重力较大，需要的起动流速亦增加，B3站位位于庙岛群岛内部且水深较浅，潮流流速较小，动力不足。A6、A7、A8 3个站位底层最大流速超过起动流速的时间较久，指示其潮流流速容易使表层沉积物再悬浮。其余站位底层最大流速超过起动流速的时间较短，短暂的时段内不足以使表层沉积物再悬浮太久。

表2 各站位表层沉积物起动流速Table 2 The critical motion velocity of surface sediment at each station

图8 各站位潮流流速、含沙量时间序列剖面Fig.8 The current velocity and the variation of sediment concentration with time in each station

3.3.3 再悬浮通量与沉降通量

为了加深对研究区悬浮泥沙特性的理解，对研究区表层沉积物再悬浮通量与悬浮泥沙沉降通量进行了计算、对比分析，以求寻得二者对研究区悬浮泥沙浓度分布及运移的影响。为了更为直观地分析出二者单日内对水体中悬浮泥沙浓度的影响，将计算得到的再悬浮通量与沉降通量区间进行了算术平均计算，同时将单日24 h按照再悬浮与沉降特点划分为再悬浮时间与沉降时间，进而求出单日每平方米的再悬浮量与沉降量（表3）。研究区再悬浮通量整体较小，庙岛海峡最小，分区差异不显著，沉降通量明显小于再悬浮通量；单日表层沉积物再悬浮量最大为43.47 kg/m2，出现在A2站位，最小为2.30 kg/m2，出现在A5站位，悬浮泥沙沉降量较小，远小于再悬浮通量，表层沉积物再悬浮量、悬浮泥沙沉降量平面变化规律与再悬浮通量的平面变化具有较高的一致性。

表3 A1—A9各站位再悬浮通量与沉降通量Table 3 Resuspension and sedimentation fluxes of stations A1-A9

有学者认为，登州浅滩细颗粒的粉砂、黏土级物质被潮流带走同时输入较粗的细砂物质，但峡道地貌形成之后，海区地貌形态基本趋于稳定，海底侵蚀量已极为有限[4]。张伟[32]发现蓬莱西海岸沿岸泥沙呈离岸方向运移的趋势，并有可能成为登州浅滩泥沙来源的一部分，同时也有学者通过卫星遥感反演大面积海域，认为渤海中部和渤海海峡海域表层悬浮泥沙主要来源于表层沉积物再悬浮及其输运[33]。还有学者认为自20世纪80年代以来，西庄至栾家口一带黄土质海岸遭受严重侵蚀[12]，季节性河流径流将数量可观的以黄土为主的侵蚀物带入海中[34]，海岸侵蚀来沙、河流输沙、邻近海域搬沙及微弱的海底侵蚀供沙共同构成庙岛海峡周边海域的主要泥沙来源[5]。据以上计算结果及讨论分析，大致可以认为研究区除庙岛海峡以外的其他海域，表层沉积物再悬浮对水体中悬浮泥沙的影响较显著，尤其研究区东南侧的泥质分布区和西部登州浅滩的西南侧细颗粒物质区，在庙岛海峡的峡道效应及庙岛群岛的群岛效应双重影响下，加之山东半岛北部沿岸流的周期性推动，表层沉积物再悬浮与底层水体水沙交换的现象更加显著。

3.4 悬浮泥沙输移特征

3.4.1 悬浮泥沙输运机制

本文所提到的悬浮泥沙输运通量计算均是指在正常天气潮流作用下的计算结果（表4—6、图10—11），未考虑波浪、风暴潮的影响因素。结果显示，各站位中单宽悬沙净输运率（T）只有A9站位小潮期大于大潮期，其余站位均是大潮期大于小潮期，大、小潮期T最大值均出现在A4站位，T最小值均出现在A9站位。大、小潮期A1、A2、A9 3个站位T与（T1+T2）的方向差别较大，其余站位T与（T1+T2）的方向基本一致。潮周期平均输运率由T1—T8决定，但各分项贡献值差别较大，在庙岛群岛海域（T1+T2）平流输运项主导着悬浮泥沙的输运，其余各项均较小，在其他海域内（T1+T2）平流输运项在研究区悬浮泥沙输运过程中起了主导作用，T5垂向净环流贡献项对研究区悬浮泥沙输运的影响也较大，（T3+T4）和（T6+T7+T8）对研究区悬沙输运贡献值很小，本文讨论分析时（T3+T4）和（T6+T7+T8）贡献项忽略不计。大潮期的A3—A6、A8、B1、B2站位和小潮期的A3站位T2贡献值大于其他站位，表明这几个站位斯托克斯漂移效应明显，这与登州水道涨落潮流作用增强相关，实测资料亦显示，此处潮流流速大于研究区内的其他海区。A7站位T2贡献值较小，主要是因为登州浅滩的存在，地形影响使A7站位处涨落潮流作用减弱。登州浅滩和登州水道之间的泥沙交换、潮周期不对称输沙以及表层沉积物再悬浮导致了垂向流速变化和悬浮泥沙浓度变化的相关项（T5）的增大。

图10 大潮期各站位悬沙输运项与平均余流叠置图注：A9的T放大10倍，A4的T1+T2和T缩放为1/2，B1的T1+T2和T放大4倍，B3的T1+T2和T放大20倍，B2放大2倍。Fig.10 The superimposition of suspended sediment transport term and average residual current at each station during spring tideNote: T of A9 is magnified 10 times, T1+T2 and T of A4 are zoom 1/2, T1+T2 and T of B1 are magnified 4 times, T1+T2 and T of B3 are magnified 20 times, and B2 is magnified 2 times.

表4 A1—A9各站悬沙输运项及单宽悬沙净输运率（大潮）Table 4 Items of suspended sediment transport and the net transport rate of single-width suspended sediment at stations A1-A9 (spring tide)

表5 A1—A9各站悬沙输运项及单宽悬沙净输运率（小潮）Table 5 Items of suspended sediment transport and the net transport rate of single-width suspended sediment at stations A1-A9 (neap tide)

表6 B1—B3各站悬沙输运项及单宽悬沙净输运率（大潮）Table 6 Items of suspended sediment transport and the net transport rate of single-width suspended sediment at stations B1-B3 (spring tide)

对于平流输运项（T1+T2）与垂向净环流项（T5）在大、小潮期对各个站位的贡献又有所不同，总体来讲，平流输运项（T1+T2）在悬浮泥沙输运中占主导作用，除了A5站位平流输运项（T1+T2）与垂向净环流项（T5）悬浮泥沙输运方向差别较大外，其余各站悬浮泥沙输运方向基本一致。庙岛海峡（T1+T2）贡献项要大于东、西两侧海域，从实测资料的分析中也可以发现（图8），垂向悬浮泥沙浓度和潮流流速的梯度变化，庙岛海峡小于东、西两侧海域，庙岛海峡东、西两侧海域余流作用相对减弱。研究区悬浮泥沙输运整体以（T1+T2）贡献项为主，T5贡献项为辅，自庙岛海峡至东、西两侧海域T5贡献项的作用逐渐增加。研究区悬浮泥沙运移的主要动力来源于潮流的周期性变化，潮流对悬浮泥沙的搬运以往复搬运为主，落潮流在悬浮泥沙输运中发挥了更大的作用。

平面上研究区大、小潮期垂向平均悬浮泥沙净输运方向和垂向平均余流方向基本一致，悬浮泥沙净输运量大小变化趋势与垂向平均余流大小变化趋势亦基本一致。研究区内对悬浮泥沙净输运贡献最大的是余流，净输运率不大，主要由水道向两侧浅滩输沙。A3、A6站位净输沙方向指向长山尾附近的浅滩，A5、A7站位净输沙方向指向登州浅滩，A5、A6站位净输沙方向相反，峡道效应起了关键性作用。水深地形的复杂变化，岛屿对潮流场的制约作用，致使悬浮泥沙在庙岛海峡北侧浅滩与登州浅滩沉降堆积。

图11 小潮期各站位悬沙输运项与平均余流叠置图注：A7、A9各项放大2倍。Fig.11 The superimposition of suspended sediment transport term and average residual current at each station during neap tideNote: A7 and A9 are magnified 2 times.

3.4.2 悬浮泥沙输运趋势

庙岛群岛海域悬浮泥沙输移方向表现为自北向南，穿过庙岛群岛浅滩又分为两个方向，一是沿着庙岛群岛西侧岸线向北输运，另一方向则是沿着庙岛海峡北侧浅滩输运至黄海方向；庙岛海峡西侧则表现出远海向近岸输沙、近岸海域向庙岛海峡方向输沙的特点；庙岛海峡东侧悬浮泥沙输运相对简单，沿山东半岛北岸向东南方向输沙的同时，自近岸向远海输沙。庙岛海峡悬浮泥沙输运相对复杂，通过图10、11，并结合数值模拟的余流场（图4）和计算的表层沉积物通量（A4、A6、A7站位净输运方向为涨潮，即向东输沙，A5净输运方向为落潮，即向西输沙）的结果，发现在庙岛海峡的南侧由渤海向黄海方向输运物质，在庙岛海峡的北侧亦有明显的自黄海向庙岛海峡方向的物质输移，并且庙岛海峡内部也有小环流进行物质交换。北黄海悬浮泥沙通过庙岛海峡北侧进入渤海，一部分向登州浅滩方向输运，另一部分在长山岛涡旋的作用下通过庙岛海峡北侧进入黄海；渤海悬浮泥沙通过庙岛海峡南侧进入北黄海，并沿着山东半岛向东和东北向输运；庙岛群岛北侧悬浮泥沙沿着长山岛与大黑山岛之间狭道向登州水道输运，穿过峡道浅滩后一部分悬浮泥沙加入小环流进行物质交换（图12）。

图12 研究区悬浮泥沙输运方向示意图Fig.12 The transportation direction of suspended sediment in the study area

根据以上探讨，初步认为研究区悬浮泥沙运移趋势整体如下：① 研究区西南侧悬浮泥沙沿着蓬莱西海岸经庙岛海峡南侧进入黄海，到达赵格庄北侧海域分为两部分，一部分悬浮泥沙继续沿着山东半岛沿岸向前运移，另一部分悬浮泥沙垂直岸线向北黄海深处运移；② 来自黄海的少量悬浮泥沙，经庙岛海峡北侧进入，在南长山岛西南侧海域与来自渤海的水体碰撞产生小环流，部分悬浮泥沙进入庙岛群岛海域，部分悬浮泥沙继续跟随环流运移；③ 来自渤海的悬浮泥沙由于东来的黄海水体的干扰，部分南下运移与西部区西南侧悬浮泥沙一起经庙岛海峡南侧运往黄海；④ 庙岛群岛北侧悬浮泥沙为群岛内悬浮泥沙输运提供了来源，输运至庙岛群岛内部的悬浮泥沙，大部分会沉降淤积在庙岛群岛内，小部分随落潮流输运至群岛外侧；⑤ 庙岛海峡悬浮泥沙输运最为复杂，不只是将渤海悬浮物质输往黄海，亦有黄海悬浮物质经庙岛海峡向渤海运移，如本文的A4、A7站位净输运方向指向黄海，A3、A5站位净输运方向指向渤海。

4 结论

（1）研究区海域悬浮泥沙浓度潮周期变化与潮流流速周期变化相关性较好，底层悬沙浓度变化对高流速的响应比较明显，表层悬沙浓度变化对低流速响应比较明显，悬沙浓度存在垂向分层现象，底层大于表层，悬浮泥沙浓度峰值整体上滞后于潮流流速峰值0.5～2 h。

（2）研究区悬浮泥沙单颗粒沉降现象不明显，除庙岛海峡外其他海域较适合悬浮泥沙絮凝沉降，并以中、底层絮凝沉降为主，且表现出自表层至底层絮凝沉降作用逐渐加强趋势；表层沉积物再悬浮对近岸浅水区、庙岛群岛周边海域水体悬浮泥沙浓度的影响要强于研究区其他海域，庙岛群岛南、北两侧表层沉积物再悬浮供沙作用强于中部；表层沉积物再悬浮量、悬浮泥沙沉降量平面变化规律与再悬浮通量的平面变化具有较高的一致性。

（3）研究区悬浮泥沙输运整体以平流输运为主，垂向净环流为辅；潮流对悬浮泥沙的搬运以往复搬运为主，落潮流在悬浮泥沙输运中发挥了更大的作用；悬浮泥沙净输运方向与余流方向大致相同，净输运率不大，主要由水道向两侧浅滩输沙。

（4）研究区西部悬浮泥沙经庙岛海峡南侧输运至黄海，东部悬浮泥沙经庙岛海峡北侧输运至渤海，二者同时进行，部分悬浮物质经海峡内的小环流进行物质间的交换，庙岛群岛北侧悬浮泥沙为群岛内悬浮泥沙输运提供了来源。