福宁湾及附近海域悬沙粒度与影响因素

毕云天，胡日军,2，陈娟娟，李毅，伊兆晗，陈晓磁，朱龙海,2，尹砚军，刘波

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

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

3. 华能霞浦核电有限公司，宁德 352100

粒度是悬浮泥沙研究的重要参数，它与泥沙搬运、沉积等过程关系密切[1-2]。悬沙粒度特征一直是海岸、河口地区沉积动力研究的重要内容之一，研究者们主要关注悬沙粒径及组分特征[3-4]、悬沙粒度与悬沙锋[5]、悬沙粒度与浓度的关系[6]、悬沙粒度与再悬浮[7-9]、悬沙絮凝体粒径及絮凝沉速等[10-13]，认为悬沙粒度变化与泥沙来源、底质再悬浮及水动力环境关系密切，含沙量和絮凝作用也会对悬沙粒度产生影响。研究悬沙粒度的变化对于分析泥沙再悬浮和沉降作用[14]、追溯沉积物的来源及其形成过程[15-16]，理解生物地球化学循环具有重要意义[17-18]，同时悬沙粒径还是影响絮凝沉降的重要因子[19]。因此，了解悬沙粒度及其变化是进一步研究泥沙再悬浮、沉降和输运规律的基础。

东海属宽广陆架海，地势平坦，岸线蜿蜒曲折，沿岸发育有长江及众多中小型河流[]，为东海提供了丰富的陆源沉积物质供应[20]。东海陆架沉积物在复杂动力因素（沿岸流、台湾暖流和潮流等）的作用下不断分异，形成了东海内陆架泥质沉积区[21]，主要分布在60 m 等深线以浅的浙闽沿岸海域[22-23]。受地形地貌、水动力条件、物质来源等因素的影响，在浙闽沿岸地区发育了一系列淤泥质海湾，这些淤泥质海湾由于其独特的沉积动力环境吸引了众多学者的关注[24-28]。

福宁湾属浙闽沿岸典型的淤泥质海湾，沉积物颗粒较细[29]。受长江及浙闽沿岸众多中小型河流的物源影响，泥沙来源复杂；岸线曲折，岬角突出，岛屿众多，近岸地形多变；近年来海湾开发活动日益加剧，人类活动影响频繁，沉积动力环境复杂。此前对于福宁湾海域沉积动力环境的研究主要集中在潮流场与冲淤特征、表层沉积物粒度组成与分布、悬浮泥沙输运特征等方面[30-32]，对于悬沙粒度的研究涉及的内容较少。本文基于2019 年12 月（冬季）和2020 年6 月（夏季）在闽北福宁湾及附近海域获得的水文泥沙实测数据，从不同时间尺度分析悬沙粒度特征，探讨影响悬沙粒度季节差异的主要因素以及悬沙粒度与浓度的关系，以加深对浙闽淤泥质海湾沉积动力环境的理解。

1 研究区概况

研究区（26°36'～26°57'N、120°05'～120°31'E）位于福建东北部福宁湾附近海域（图1），湾内水深不足12 m，湾外水深最大不超过30 m[32]。沉积物类型主要为粉砂和黏土质粉砂，中值粒径大于6.5 Ф。潮汐类型为规则半日潮，以M2分潮为主，平均潮差大于4 m[33]。湾内潮流往复性明显，涨潮流向多为NWW 向，落潮流向多为SEE 向，由湾内向湾外潮流的旋转性逐渐加强[31]。研究区附近无大型河流注入，仅发育有罗汉溪、长溪三河等源短流小的山溪性小河。研究区季风特征明显，冬季盛行强而稳定的偏北风，夏季则以偏南风为主，受季风的影响，冬季波高明显大于夏季[34]。低温低盐的闽浙沿岸流主要分布在长江口以南的浙闽沿岸，在冬季偏北风期间为一支较强的南向流；夏季偏南风期间，闽浙沿岸流流速很小，甚至消失[35]。高温高盐的台湾暖流在闽浙沿岸流的东侧海域终年向北流动，位置稳定，年际变化不大[36]。

图1 研究区附近流系（a）及调查站位分布（b）附近流系改自文献[37]，MZCC：闽浙沿岸流，TWC：台湾暖流。Fig.1 Distributions of current systems and stations in the study areaNearby current systems are modified from the reference [37], MZCC: Min-Zhe Coastal Current, TWC: Taiwan Warm Current.

2 数据与方法

2.1 数据来源

浙江省水利河口研究院在闽北福宁湾及附近海域布设了8 个同步观测站（图1），分别在冬、夏季大小潮共进行了4 次连续26 h 的水文泥沙观测（表1），4 次观测站位均相同。观测内容包括水深、海流（流速、流向）、含沙量、温度、盐度、悬沙和底质粒度、海面基本气象参数（气温、海况、风速风向）。

表1 各测次观测时间统计Table 1 Time statistics of four observations

水深和海流的观测使用美国RDI 公司的剖面流速仪ADCP，观测间隔为30 min，盲区0.5 m，层宽根据水深情况设置成0.25～1.0 m 不等，水深直接采用经过吃水改正的ADCP 水深。按照6 点法（表层、0.2H、0.4H、0.6H、0.8H、底层）提取流速、流向数据。

与测流同步使用XCL 横式取样器在上述层位采集含沙量水样，每小时（整点）取样一次，每次采水量不少于1 L；水温和盐度的观测使用SST CTD75型温盐深仪，观测时间、观测层位均与含沙量观测同步；另外，在白天潮的4 个潮流特征时段（涨急、涨憩、落急、落憩）各取一个悬沙粒度分析样品，取样层位为0.6H层（中下层），单个水样体积为10 L。采用抓斗式取样器在每个调查站位各取一个底质粒度分析样品。

2.2 分析方法

2.2.1 含沙量及粒度测试

含沙量的室内分析采用烘干称重法，具体步骤如下：① 烘干烘杯，先将烘杯洗净，放入温度为100～110℃烘箱中烘2 h，后移入干燥器内冷却至室温，再用1/10 000 电子天平称重，得到烘杯重。②水样经过48 h 沉淀后，抽取部分样品上层清液（一般为涨憩时）进行氯化物分析，根据氯化物浓度决定是否对样品进行洗盐处理。凡氯化物浓度超过250 mg/L 的，均作洗盐处理。③ 用少量清水将浓缩水样全部冲入烧杯中，加热至无流动水时，移入烘箱，在温度为100～110℃烘箱中烘干，烘干所需时间由试验确定。④ 烘干后的沙样，及时移入干燥器中冷却至室温，用1/10 000 电子天平称重后减去烘杯重，得到每个样品的含沙量。

悬沙及底质粒度测试采用英国Malvern 公司的Mastersizer 3000 激光粒度仪。测量范围为0.02～2 000 μm，中值粒径测量误差不超过1%，实际测量中校正系数、平均粒径、分选系数的内检误差均满足要求。粒度参数采用McManus 公式[38]计算，各个粒度参数评价采用McManus 矩值法粒度参数等级表。沉积物命名采用Shepard[39]分类方法。

2.2.2 水体层化指标

梯度Richardson 数[40]（Ri）可以用来判断水体的层化与湍流混合状态，如式（1）:

式中，z为水深（m），ρ为水体密度（kg/m3），u为流速北分量（m/s），v为流速东分量（m/s），g 为重力加速度（取9.8 m/s2）。一般认为当Ri＞0.25 时，水体处于稳定的层化状态，湍流被抑制；当Ri＜0.25 时，层化效应随着Ri 减小而减小，流速剪切更容易诱发湍流混合。

2.2.3 悬沙输运通量计算

根据泥沙输运通量公式[41]，分别计算冬、夏两季8 个站位26 h（两个潮周期）单宽面积上的潮平均纵向净输沙率Ts（g/s），如式（2）、（3）：

式中，k=1,…,6，为各层位编号，Uk（m/s）为各水层内流速u的深度平均值，Ck（g/m3）为各水层悬沙浓度的深度平均值，Δzk（m）为各层位高度，T为潮周期（h），Tr(k)（g）为潮周期内穿过每一个水层单位宽度面积的净输沙量。

3 结果

3.1 水文泥沙特征

冬季观测期间，研究区平均风速一般不超过6 m/s，常风向为NE 向，以二级海况为主，波高小于0.5 m，属于较平静天气；夏季观测期间，研究区平均风速一般不超过4 m/s，常风向为SW 向，以一级海况为主，波高小于0.1 m，属于平静天气。与冬季相比，夏季观测期间天气状况较好。

3.1.1 潮流

观测期间潮流流速整体不大（表2），冬季研究区平均流速，大潮期为0.32 m/s，小潮期为0.28 m/s；夏季研究区平均流速，大潮期为0.33 m/s，小潮期为0.22 m/s，流速的季节变化不大，大小潮差异明显。潮流在湾内较弱，在湾外较强，冬、夏季最大流速均出现在湾外7#站，分别为0.74 m/s 和0.92 m/s。

表2 研究区流速特征值Table 2 Characteristic value of current velocity in the study area m/s

3.1.2 悬沙浓度

研究区冬季悬沙浓度为0.020～0.897 kg/m3，夏季为0.060 kg/m3～0.178 kg/m3；各站位垂线平均悬沙浓度，冬季为0.060～0.178 kg/m3，夏季为0.009～0.032 kg/m3。与冬季同时期相比，夏季悬沙浓度大幅减少（图2），据统计，各站位冬季大潮悬沙浓度是夏季大潮的5.6～11.6 倍，冬季小潮悬沙浓度是夏季小潮的3.9～13.2 倍。大小潮变化上，大潮悬沙浓度大于小潮悬沙浓度。

图2 不同潮周期、不同季节悬沙浓度对比Fig.2 Comparison of suspended sediment concentration (SSC) in different tides and seasons

3.1.3 水体层化

本文利用取对数后的梯度Richardson 数（lg(Ri/0.25)）来表征水体层化状态。当lg(Ri/0.25)＜0，表明层化较弱，水体处于湍流混合状态；当lg(Ri/0.25)＞0，表明层化较强，水体湍流混合受到抑制。本文在冬、夏季大潮期各取典型站位（2#、4#、7#）进行分析，并给出了lg(Ri/0.25)与悬沙浓度随时间变化的剖面图（图3）。图3 显示，研究区水体层化具有明显的季节差异，冬季水体垂向混合较强，lg(Ri/0.25)＞0 的区域仅出现在个别观测时间段的某一层位；夏季水体层化明显，lg(Ri/0.25)＞0 的区域由水体表层向下延伸至近底层，并在潮周期内几乎呈连续分布。

图3 大潮期2#、4#、7#站lg(Ri/0.25)与悬沙浓度时间序列Fig.3 lg(Ri/0.25) and SSC time series of 2#, 4# and 7# stations in spring tide

夏季水体受层化结构的影响，泥沙在垂向上的扩散受到抑制。如图3 所示，夏季各站位强层化区（lg(Ri/0.25)＞0）基本与悬沙浓度低值区相对应，0.01 kg/m3等值线随着层化边界的变化上下移动，层化区域内并未观察到悬沙浓度的潮周期变化，底层泥沙基本不能穿过水体层化边界进入中、上层。相反，冬季水体基本不存在层化结构，悬沙在垂向上的扩散较为充分，整个剖面上悬沙浓度的潮周期变化显著。

3.2 悬沙粒度特征

3.2.1 悬沙粒度的时空分布

受现场观测条件的限制，本文所分析的悬沙粒度数据，均取自水体0.6H层（中下层）。粒度分析结果显示，正常天气下，研究区悬沙粒度较细，中值粒径普遍大于6.5 Ф（表3）。粒度组成以粉砂为主，黏土次之，两者含量之和达90%，基本不含砂或砂含量很小，分选较差—差，峰态呈很宽平分布，偏态以近对称为主（表4）。

表3 研究区悬沙中值粒径变化范围Table 3 Variation range of suspended sediment median grain size in the study area Φ

表4 悬沙粒度参数Table 4 Suspended sediment grain size parameters

（1）平面分布

悬沙粒径平面分布见表3，冬季各站位中值粒径差异较小，平均中值粒径为7.36～7.53 Ф。湾内和湾口中值粒径较为接近，平均中值粒径分别为7.48 Ф和7.45 Ф，湾外平均中值粒径为7.36 Ф，与湾口和湾内相比略有增大。夏季各站位中值粒径差异较大，平均中值粒径为6.45～7.31Ф，中值粒径最大值出现在湾外8#站，最小值出现在湾口6#站。

（2）季节变化

悬沙粒度季节变化明显（表4），表现为冬季细、夏季粗的季节分布特征。冬季中值粒径为7.1～7.7 Φ，均值为7.45 Φ。4～8 Φ（4～63 μm）的粉砂含量最丰富，为57.2%～71.0%，并且含有一定量8～12 Φ（＜4 μm）的黏土，为26.8%～42.6%，砂含量最少，不超过10%。夏季，中值粒径为5.7～7.5 Φ，均值为6.97 Φ。粉砂、黏土、砂含量分别为53.1%～79.2%、12.6%～37.4%、0.8%～20.8%。与夏季相比，冬季悬沙中的黏土含量明显增加，从25.9%增加到34.6%，粉砂和砂含量分别从68.7%和5.5%减少到64.1%和1.3%。

与夏季相比，冬季悬沙中增加的细粒级部分主要为粒径＞7 Ф（＜8 μm）的极细粉砂和黏土，粒径为6～7 Ф（8～16 μm）的细粉砂含量变化不大，粒径＜6 Ф（＞16 μm）的粒级部分显著减少。从8#站冬、夏季大潮期悬沙粒径分布来看（图4），冬季悬沙粒径分布整体向细粒方向移动，悬沙分选变好，峰态变窄。

图4 大潮期8#站悬沙粒径分布Fig.4 Suspended sediment grain size distribution of 8# station in spring tide

（3）大小潮变化

受潮动力变化的影响，大潮悬沙中值粒径明显大于小潮（表4）。冬季大、小潮期间中值粒径均值分别为7.37 Φ 和7.53 Φ，大潮约为小潮的1.1 倍；夏季大、小潮期间中值粒径均值分别为6.77 Φ 和7.17 Φ，大潮约为小潮的1.3 倍。冬季大潮期间由于粗颗粒物质的加入，砂、粉砂含量增加，黏土含量相应减少，悬沙粒度变粗，与小潮相比悬沙分选变差，峰态变宽。

夏季则表现出不同的规律，大潮期间粉砂组分增加尤为明显，含量平均增加约为8%，黏土和砂组分均有不同程度的减少，悬沙粒径分布与小潮相比更为集中，悬沙分选变好，峰态变窄。粉砂组分具有起动流速小、沉降速度大的特征，在潮流和波浪作用下，泥沙运动活跃，易于悬浮和沉降[42]。夏季大潮期间流速增加明显，水体挟沙力增强，“易悬易沉”的粉砂组分大量进入水体，对含沙量的贡献显著；对于粒径较粗的砂组分，不仅起动需要较高的流速，而且在搬运过程中，由于流速减小也易于沉降，所以很难长时间保持悬浮状态，多数是处于间歇性的悬浮状态；黏土组分的黏聚力强，难以起动，因此悬沙中黏土组分的含量也较少[5]。

（4）潮周期变化

由于本文在潮周期尺度上只在白天潮的4 个特征时段（涨急、涨憩、落急、落憩）获得了悬沙粒度数据，悬沙粒径在时间上精度较低，如果对每个站位单独分析，则很难观察到明显的变化规律，因此，本文将所有站位作为一个整体分析。

考虑到潮周期内悬沙粒径变化可能与水动力有关，本文对中值粒径与流速进行相关性分析（图5），结果表明，冬季大、小潮期间，随着流速的增大，中值粒径只表现出微弱增大的趋势，夏季两者的趋势性关系更为复杂。无论冬季还是夏季，中值粒径与流速的对应关系较为散乱，两者并没有表现出显著的相关性。悬沙粒径变化与泥沙的再悬浮和沉降过程有关[43]，流速虽然是重要的影响因素，但在多数情况下它并不是主导因素[44]。本次悬沙粒度取样层位位于0.6H层，泥沙从底床向上扩散需要一定的时间，当流速从极大值开始减小，但还没有低于再悬浮的临界流速，泥沙也不会立刻发生沉降，这可能会造成悬沙粒径的变化滞后于流速变化，同时，这一过程还会受到水体混合状态的影响，较为复杂的动力过程可能是导致悬沙粒径与流速没有明显相关关系的重要原因。

图5 悬沙中值粒径与流速的相关性Fig.5 Relationship between suspended sediment median grain size and current velocity

3.2.2 粒级-标准偏差曲线

从全部样品各粒级含量的粒级-标准偏差曲线可以清楚地看出粒度组成的变化，其理论研究基于历史上沉积环境变化引起颗粒物的粒度差异，这种差异在受沉积动力和物质来源显著影响的现代沉积环境中理应存在[45]。标准偏差大的粒级组分，其在各样品中的含量变化大，可以作为指示海洋环境变化的敏感粒级[46]。

本文对冬、夏季悬沙样品的粒度数据分析后得出粒级-标准偏差曲线（图6），可以看出，冬季粒级-标准偏差曲线呈现明显的双峰分布，对环境最为敏感的两个高峰值粒级分布范围分别为2.8～3.4 μm和21.1～25.1 μm，其中2.8～3.4 μm 粒级成分为粗黏土，21.1～25.1 μm 粒级成分为中粉砂。低峰值粒级分布范围为10.8～13.2 μm，对环境的敏感度最差，粒级成分为细粉砂。这种分布特征与前人在浙闽沿岸泥质沉积区柱状样以及表层沉积物中获得的敏感粒级分布十分接近[47-48]。与冬季不同的是，夏季每一粒级的标准偏差都要大于冬季，粒级-标准偏差曲线虽然变化趋势与冬季基本一致，但双峰分布特征并不明显，在细敏感粒级处表现为一个不明显的次峰，这意味着控制此粒级的海洋环境要素可能发生了变化。

图6 悬沙的粒级-标准偏差曲线Fig.6 The standard deviation curve of suspended sediment grain size

3.3 泥沙输运通量

为探究研究区泥沙净输运趋势，本文利用同步得到的悬沙浓度和流速剖面计算了两个半日潮周期（26 h）内各站位单宽面积上的潮平均纵向净输沙率。结果表明，冬季大、小潮期间，由湾内向湾外净输沙率不断增大（图7a、图7b），湾内1#、2#站净输沙率较小，量值均在35 g/s 以下，方向总体指向近岸；湾口3#、4#、5#、6#站净输沙率与湾内相比有明显增大，量值为42.9～137.9 g/s，大致沿岸由东北指向西南；湾外7#、8#站净输沙率最大，量值均超过100 g/s，方向与湾口站位基本一致。这一特征显示，冬季研究区悬沙总体输运趋势为由东北向西南，这与冬季闽浙沿岸流的方向基本一致。薛碧颖等[21]在闽北近岸海域冬季悬浮体输运通量的研究中也得到了类似的结果，认为冬季闽浙沿岸流受东北风的驱动，沿岸向西南方向流动，使得近岸悬浮体的输运方向总体指向西南。

图7 各站位观测26 h 悬沙净单宽输运通量Fig.7 Sediment fluxes at 8 stations based on 26-hour observations

夏季大、小潮期间，净输沙率与冬季相比显著减小（图7c、7d），各站位量值均在10 g/s 以下，由湾内向湾外净输沙率逐渐增大，方向总体由东南指向西北。这与薛碧颖等[21]在夏季的研究结果存在一定的出入，这可能与夏季闽浙沿岸流强度减弱有关，所计算的净输沙率在近岸主要受潮流作用的影响。

4 讨论

4.1 悬沙粒度季节变化的影响因素

研究区悬沙粒度具有冬季细、夏季粗的季节差异，冬季由于悬沙中细颗粒物质（主要指极细粉砂和黏土组分）的相对增加，粒度细于夏季。底质再悬浮和泥沙来源是造成悬沙粒度季节差异的重要因素[49]，本文着重对这两个影响因素进行分析，探讨其对悬沙粒度季节变化的影响。

4.1.1 底质再悬浮

据本文统计，冬季研究区平均流速约为0.30 m/s，夏季约为0.28 m/s，潮流动力的季节变化很小。研究区水深较浅，冬季在偏北风的作用下波浪成长迅速，据相关研究统计，福宁湾湾口长表岛附近海域冬季H1/10波高平均值约为1.10 m，夏季约为0.87 m[50]，冬季波浪动力强于夏季。

研究区沉积物类型以黏土质粉砂为主（表5），粒度较细，中值粒径均值约为7 Ф（16 μm），粒径＜30 μm 的黏性泥沙，累计含量占比为84.6%～94.0%，并含有相当比例的黏土组分（＜4 μm），含量普遍超过25%，属于淤泥质泥沙的范畴[42,51]。不同粒径的泥沙在潮流和波浪作用下的起动遵循Shields 曲线[52-53]，即当粒径大于100 μm 时，起动流速随粒径增大而增大；当粒径小于100 μm 时，由于细颗粒间的黏聚力作用，起动流速随粒径减小而增大。由于黏土物质的黏聚力强，其再悬浮需要较强的水动力条件。由本文3.1.3 节可知，夏季水体普遍存在层化结构，同时在较弱的波浪作用下沉积物再悬浮微弱，研究区平均悬沙浓度低于0.020 kg/m3。冬季水体湍流混合和波浪作用与夏季相比均明显增强，动力增强并克服黏聚力，再悬浮了更多的细颗粒物质，悬沙粒度整体变得更细[54]。

表5 研究区底质粒度特征Table 5 Characteristics of sediment grain size in the study area

此外，波浪作用深度还与波高呈正相关[55]，冬季更强的波浪可以作用于更深的海床。研究区表层沉积物中值粒径随着水深的增加而减小（表5），当波浪增强时，其作用深度增大，海底深处粒度更细的沉积物就会悬浮起来，这是冬季悬沙粒度细于夏季的另一个重要原因。

4.1.2 泥沙来源

研究表明，浙闽沿岸诸河、长江以及海岸基岩都是福宁湾附近海域沉积物的潜在物源[56-57]。福建地区多发育侵蚀型海岸[58-60]，研究区范围内较粗的砂质粉砂和粉砂等沉积物，主要来源于在较强波浪作用下海岸基岩风化产生的粗颗粒物质[50]。福宁湾沿岸还有罗汉溪、长溪三河等源短流少的山溪性小河，以及东冲半岛几条季节性小溪流，流量变化大，每遇洪水季节，有大量泥沙冲入湾内，这些泥沙一般颗粒较粗，以细砂、中砂为多，搬运距离不远，主要在湾顶局部区域堆积。沉积物中的细颗粒物质主要为外来输入物质，系长江入海泥沙与浙闽沿岸诸河的泥沙经搬运后在此沉积[61]。粒级-标准偏差曲线显示，悬沙中对环境最为敏感的两个粒级组分分布范围分别为2.8～3.4 μm 和21.1～25.1 μm，这两种主要敏感粒级的出现，应是受上述粗细不同的物质影响，其中粗敏感粒级组分可经落潮流由岸边带入湾内，后经再悬浮进入水体；细敏感粒级组分为黏土，来源于异地河流输入物质（长江与浙闽沿岸诸河沉积物），东海悬浮体“夏储冬输”的季节性输运格局决定了到达本海域异地河流来源物质的输入强度与闽浙沿岸流的强度有关，具有明显的季节性变化[47,62]。

闽浙沿岸流是东海内陆架泥沙输运的重要动力[63-64]，起源于长江口一带，主要分布在长江口以南浙闽沿岸50 m 以浅的海域[65]。冬季，长江口与杭州湾附近海域受潮流、波浪作用显著，源于洪季就近沉积下来的细颗粒泥沙发生强烈的再悬浮，悬沙浓度急剧升高，部分区域表层悬沙浓度可达3 kg/m3[66-67]；大量细颗粒泥沙被闽浙沿岸流夹带南下[68-70]，同时，沿岸流在南下过程中还能携带浙闽沿岸诸河沉积物，最远能到达泉州湾附近[71]。这些泥沙在搬运过程中由于重力分异作用逐渐细化，曾有研究统计，在杭州湾附近，悬沙中值粒径为9.9 μm，当被搬运到福宁湾北侧的沙埕港附近，悬沙中值粒径已减小至3.7 μm[72]，达黏土级范畴，与粒级-标准偏差曲线上的细敏感粒级基本对应，呈现出明显的峰值。夏季，闽浙沿岸流势弱，台湾暖流势强，研究区南部闽江入海的细颗粒物质也会随着终年北上的台湾暖流到达研究区附近，但其量值大大减少[56]，在粒级-标准偏差曲线上只呈现出一个不明显的峰值。泥沙净输运通量的计算结果也显示（图7），冬季研究区附近海域悬沙输运总体趋势与闽浙沿岸流方向基本一致，悬沙由东北向西南，由湾外向湾内输运，而在夏季悬沙输运方向转为由东南向西北，输运通量与冬季相比显著减小。

从冬、夏季细粒级组分的级配（＜63 μm）差异来看（图8），冬季总细粒级悬沙中平均有65.7%的＜8 μm 粒级的颗粒，夏季＜8 μm 粒级的颗粒只占总细粒级悬沙量的53.1%，同时冬季含沙量也远高于夏季，因此＜8 μm 粒级的细颗粒较夏季增加非常可观。其中，除本地再悬浮的物质外，闽浙沿岸流向研究区输送的黏土级颗粒（＜4 μm）也占重要组成部分，进一步加大了悬沙粒度“冬细夏粗”的季节差异。

图8 研究区悬沙细粒级部分（＜63 μm）级配Fig.8 Percentage of fine suspended sediment (＜63 μm) in the study area

4.2 悬沙粒度与浓度的关系

4.2.1 中值粒径与悬沙浓度

就单个站位而言，悬沙粒径和悬沙浓度易受天气状况和波浪等因素的影响，两者趋势性关系不明显，而将多个站位整体分析，两者能表现出较好的相关关系[73]，因此本文将研究区作为一个系统整体进行分析。相关性分析结果显示（图9），夏季大潮期，中值粒径与悬沙浓度在0.99 的置信水平上呈显著的线性相关，悬沙粒径随悬沙浓度的增加而减小，小潮期这种线性关系依旧存在。相反，冬季中值粒径与悬沙浓度的相关性较差，随着悬沙浓度的增加，悬沙粒径只呈现出微弱增大的趋势，这种季节差异可能与悬沙来源和混合作用有关[6]。如前所述，冬季悬沙部分来自较强的底质再悬浮，部分来自闽浙沿岸流向研究区的输送，悬沙粒度组成较为复杂；同时悬沙在垂向上的混合较为剧烈（图3），悬沙浓度在潮周期内的相对变化较大，增加了悬沙粒度组成变化的复杂性，导致中值粒径与悬沙浓度没有显著的相关关系。夏季，由于闽浙沿岸流强度大大减弱，其向研究区的泥沙输送基本中断，水体中的悬沙主要来自本地沉积物的弱再悬浮作用，悬沙粒度组成相对单一；同时夏季悬沙在垂向上的混合较弱（图3），悬沙浓度在潮周期内的相对变化较小，悬沙粒度组成变化相对简单，中值粒径与悬沙浓度存在显著的线性关系。

图9 悬沙中值粒径与浓度的相关性Fig.9 Relationship between suspended sediment median grain size and SSC

值得注意的是，夏季无论大小潮，悬沙粒径均随悬沙浓度的减小而增大（图9b）。通常认为，当水动力减小，悬沙以单颗粒的形式发生重力分选沉降时，由于较粗的颗粒首先沉降，悬沙粒径将逐渐减小，以响应动力条件的改变[74]。相反的是，夏季悬沙中值粒径随着悬沙浓度的减小而增大。研究表明，在以黏性细颗粒（＜30 μm）为主的悬沙中会发生动水絮凝沉降，且表现为粒径越细，絮凝越强[75-76]。研究区悬沙中约有90%的颗粒集中在30 μm 以下，当挟沙水流小于某一临界流速时，悬沙发生动水絮凝沉降，与重力分选沉降（以单颗粒的形式由粗到细依次沉降）不同，悬沙中的细颗粒组分相互结合在一起，大小混杂同时落淤[77]。由于沉降的主体物质是细颗粒，悬沙粒径会随着这些细颗粒的沉降而增大[78]。

4.2.2 组分与悬沙浓度

悬沙浓度发生变化，悬沙组分必然会发生改变，不同组分对悬沙浓度变化的响应程度也会有所差异。研究区悬沙以粉砂为主，含量高达70%，对悬沙浓度的贡献度最大，但相关性分析结果表明（图10），在夏季大潮期，粉砂组分含量与悬沙浓度的相关性较差，粉砂并非是响应悬沙浓度变化的敏感组分；相反，砂和黏土在悬沙中的含量虽然较少，但与悬沙浓度的相关性较好，砂与悬沙浓度在0.99 的置信水平上呈显著负相关，其含量随着悬沙浓度的减小而增加，而黏土与悬沙浓度在0.99 的置信水平上呈显著正相关，其含量随着悬沙浓度的减小而减少，这种现象在小潮期依旧存在。黏土组分颗粒较细，絮凝较强，比其他组分更容易结合形成絮凝体而沉降[76]。当沉降发生时，由于黏土组分作为细颗粒絮凝沉降的重要部分，随着沉降的进行其含量趋于减少。相反，砂组分随着沉降的进行含量却趋于增加，推测砂组分可能并没有参与细颗粒絮凝沉降过程，其含量变化是由悬沙中其他组分含量变化引起的。

图10 夏季悬沙各组分与浓度的相关性Fig.10 Relationship between suspended sediment populations and SSC in summer

5 结论

（1）正常天气下，研究区悬沙粒度较细，中值粒径普遍大于6.5 Φ。悬沙粒度具有冬细夏粗的季节分布特征，冬季＞7 Ф（＜8 μm）粒级部分与夏季相比增加明显；受潮动力影响，悬沙粒度大小潮变化明显；潮周期内悬沙粒径变化规律不明显。悬沙粒级-标准偏差曲线显示两种主要的敏感粒级分布范围为2.8～3.4 μm 和21.1～25.1 μm。

（2）底质再悬浮和泥沙来源是造成悬沙粒度季节差异的主要原因。冬季较强的波浪和水体垂向混合导致大量细粒沉积物的再悬浮，以及闽浙沿岸流携带的细颗粒泥沙的输入共同造成了研究区悬沙“冬细夏粗”的季节差异。

（3）冬季悬沙粒径与浓度不存在相关性，夏季两者存在显著的线性关系，悬沙粒径随悬沙浓度的增加而减小，这种季节差异可能与悬沙来源和混合作用有关。夏季悬沙粒径的粗化可能与细颗粒发生絮凝沉降有关。不同组分对悬沙浓度变化的响应程度不同，夏季，黏土组分与悬沙浓度呈显著正相关，砂组分与悬沙浓度呈显著负相关，粉砂组分对悬沙浓度变化的响应最差。