黄河口潮滩泥沙絮凝研究

于上，谢卫明*，何青，王宪业，赵中豪，郭磊城，徐凡

( 1. 华东师范大学 河口海岸学国家重点实验室，上海 200041)

1 引言

絮凝是泥沙的重要特性，小颗粒泥沙会团聚成为大的絮凝体，进而改变沉降速度等特性[1–6]。在河口海岸区域，超过90%的泥沙都是以絮凝体的形式运输[7]。研究泥沙絮凝特性及影响因素，对深化认识河口海岸泥沙输运规律有重要意义。已有诸多研究对天然水体中的絮凝特性与影响因素进行了探讨：陈锦山等[8]、程江等[9–10]和唐建华[11]利用LISST 观测了长江口细颗粒泥沙的絮凝特性，发现长江口絮团平均粒径为86.3 μm，受水动力影响明显；Xia 等[12]、邓智瑞等[13]、林建良等[14]对珠江口絮凝体特性也进行了研究，发现珠江口絮团粒径略大于长江口，絮团粒径受盐度影响较小但受动力影响较大；长江中下游淡水环境中也发现有絮凝现象的存在，淡水絮团的粒径略低于河口絮团[15–16]。综合现场观测与室内实验发现，影响絮凝的因素包括水体紊动剪切强度、悬沙浓度、盐度、泥沙粒径、泥沙级配以及生物作用等[8,13,17–20]。

过往对絮凝的现场观测研究主要集中在河流与河口区域，对潮滩区域絮凝特性的研究相对较少。但是另一方面，潮滩上存在明显的絮凝现象，且絮凝对潮滩泥沙运移与冲淤变化有明显影响[21–23]。因此，研究潮滩区域泥沙絮凝对进一步揭示河口海岸区域泥沙运动规律有重要意义。近年来，已经有一些研究者在美国、荷兰的潮滩上进行了絮凝观测研究，但是这些研究普遍集中于中潮（潮差在2～4 m 之间）或强潮河口（潮差大于4 m），鲜见在弱潮河口潮滩进行的泥沙絮凝观测[21,23–24]。

基于以上背景，本文选取黄河口南部潮滩开展絮凝研究。黄河是著名的高含沙量河流，2019 年黄河口利津站实测输沙量为2.71×108t，大量泥沙在河口输移沉降。黄河口属于弱潮河口，潮差在0.7～1.7 m 之间[25]。因此，黄河口区域适于开展泥沙絮凝研究。通过现场观测获取黄河口絮团粒径、悬沙浓度、波浪、潮流动力等数据，可以研究弱潮河口波流共同作用下潮滩泥沙的絮凝规律与控制因素，丰富对我国河口海岸地区水沙运动规律的认识。

2 数据和方法

2.1 研究区域概况

研究区域选定在黄河口南部、莱州湾西岸的潮滩上，观测点位于垦东大堤外侧，属于中国山东省东营市（图1a）[26]。该区域潮滩坡度较小，退潮时有大片潮间带露出水面。研究区域的波浪主要受季风控制，全海区的风浪以波浪为主[25]。潮汐类型属于不正规混合半日潮，平均高潮位为1.6 m，平均低潮位为0.6 m（理论深度基准面）。该区域沉积物来源既包括黄河携带的陆相泥沙，也包括涨潮流携带的海相泥沙[27–28]。对采集观测架附近表层沉积物于实验室内进行粒度测量，测得中值粒径为64.66 μm。

图1 研究区域及仪器布置示意图Fig. 1 Sketch of study area and instrument layout

2.2 观测方法与过程

实施现场观测的时间为2019 年11 月18–23 日。观测点位于黄河口南部、莱州湾西岸潮滩上，距垦东大堤垂线距离约2 400 m，经纬度坐标为37°34′40″N，118°57′01″E。观测仪器被固定于观测架上，安置仪器前将观测架锤入泥滩约1 m 以尽可能减小振动，避免对仪器观测产生影响（图1b）。

观测架上仪器布置情况如下：ADV，测定观测期间的水体紊动数据及滩面高程变化；OBS，测定浊度与水深数据；LISST-200（Laser In Situ Scattering and Transmissometry），采集现场悬浮体粒径及体积浓度数 据；RBRsolo³ D|wave16（示 意 图 中 简 记 为RBR），采集波浪数据。以上仪器数量均为1 台。OBS、RBRsolo³ D|wave16 及LISST-200 探 头 距 床 面 高 程 均为5 cm；ADV 探头距床面高程为15 cm，该仪器的水动力采样点被设置为探头下方10 cm，即测量距床面高程为5 cm 处的水体紊动数据。4 台仪器时间同步，采样间隔均为5 min。各仪器探头位置接近，可以认为反映同一点的水沙特征。

安置及回收仪器时，均同步采集了观测站附近的表层沉积物，用于分析床沙粒径组成及对OBS 进行标定。研究时间内风速、风向收集自山东东营垦利气象站（37°34.8′N ，118°33′E）实测数据。

2.3 数据处理与计算方法

2.3.1 悬沙浓度

观测期间的悬沙浓度（SSC）根据OBS 测定的浊度（NTU）数据反演后获得。在实验室内使用现场获取的表层沉积物配置不同浓度的悬沙浊液并测定对应的浊度，获得浊度与含沙量间的相关关系[29]。标定结果显示悬沙浓度与浊度间关系可以用关系式SSC=0.006 9 NTU–0.240 8 表 示， 相 关 系 数 的 平 方R2=0.89（数据点数量N=13）表示。

2.3.2 表层沉积物粒度

安置及回收仪器时各采集观测站附近的表层沉积物样品2 份，在实验室中每份样品测验2 次，共测验8 次。测验时随机从样品袋中取出5 g 左右泥沙，加入分散剂六偏磷酸钠之后超声波震荡30 min，使絮团充分破碎之后，使用Coulter LS-100Q 激光粒度仪测定[30–31]。计算8 次测验所得的中值粒径平均值，作为观测站处表层沉积物中值粒径，计算结果为64.66 μm。

2.3.3体积浓度及中值粒径

本研究使用美国Sequoia 公司的现场激光粒度仪LISST-200来观测天然水体中的絮团。射入水中的激光束会被悬浮颗粒所散射，颗粒物粒径不同，对激光的散射角度也不同。大颗粒以小角度散射、小颗粒以大角度散射[1]。LISST-200 将向前散射的角度分为36 级，分别对应1～500 μm 之间的36 个粒级。运用Mie 散射理论，可以反推得到36 个粒级的悬浮体的体积浓度，并计算得到悬浮体中值粒径[1,9]。综合前人研究发现，潮滩环境中约80%的泥沙以絮团的形式输运[1,9]。综合郭超和何青[16]利用LISST-100 在荷兰Scheldt 河口Kapellebank 潮滩、邓智瑞等[13]在珠江口磨刀门海域、程江等[10]在长江口海域、Guo 和He[15]在长江流域的观测实践方法，以上在不同水体环境中开展的絮凝研究将LISST 观测到的悬浮体视为絮团。因此在本研究中，也将LISST-200 的观测数据整体上视为絮团进行处理。

2.3.4有效密度

有效密度是反映絮团结构的重要参数，且对沉降速度有明显影响[4,17]。本研究采用Fettweis 方法计算有效密度[3]，公式为

式中，ρw与ρs分别为水与组成絮团的分散态泥沙颗粒密度（单位：kg/m3），分别取为1 030 kg/m3与2 650 kg/m3；SSC 为悬沙浓度（单位：kg/m3）；TVC 为颗粒总体积浓度，即36 个粒级组的絮团体积浓度之和（单位：μL/L）。

2.3.5沉速

沉速可以认为是絮团重力W与下沉过程中受到的阻力F达到平衡时等速沉降的值，采用Stokes 公式计算[5]：

式中，Dm为絮团平均粒径（单位：μm）；μ为水的动力黏滞系数，根据水温选取相应的值（单位：Pa·s）；g为重力加速度，取为9.8 m/s2。

2.3.6紊动剪切率

在泥沙絮凝研究中，通常使用紊动剪切率G（s-1）来描述水体的紊动强度，其计算方法为[17,32]

式中，u*为摩阻流速（单位：m/s-1）；z和H分别为距离床面高程与水深（单位：m）；ν为水体运动黏度，根据水温选取相应的值（单位：m2/s）；κ为卡门常数，取为0.4。水流与波浪共同引起的床面切应力τcw（单位：N/m2）为[21,33]

式中，τc为水流引起的切应力，使用改进的湍流动能法（Turbulent Kinetic Energy，TKE）计算[34]：

τw为风浪引起的底部剪切力，计算方法为[35–36]

式中Aδ为波质点轨道流速（单位：m/s）；Hs为有效波高；k为波数；Ts为波浪有效周期（单位：s）；L为波长（单位：m）；fw为波浪摩擦系数；Rew与r分别为波浪雷诺数（无量纲）与相对糙率（无量纲），分别由下式求得[35–36]：

式中，ks为床面粗糙高度，取2.5d，d为潮滩表层沉积物中值粒径。根据实验室测验结果，观测期间天气较为平静，潮滩表层沉积物粒度没有明显变化。d取观测前、后的表层沉积物中值粒径平均值，为64.66 μm。根据式（4）至式（13）计算得到床面剪切力τcw之后，可由式（14）求得摩阻流速，并代入式（3）得到紊动剪切率[5]。

2.3.7临界侵蚀切应力及临界沉降剪切应力

根据前人在黄河口潮滩的研究，李华国等[37]提出的模型能够较好地预测潮滩表面沉积物的临界剪切应力：

式中，c为黏结力系数，取2.9×10-5kg/m；γd、γdc分别为表层沉积物的干容重及稳定干容重（单位：kg/m3）；其余符号含义同前[37]。γd的值通过“环刀法”采样之后在实验室中测量，为1 490 kg/m3；稳定干容重γdc指经过长期的一般压密后的沉积物干容重，根据张耀哲和王敬昌[38]提出的经验公式计算为

式中，γs为泥沙颗粒干容重（单位：kg/m3），取2 650 kg/m3；其余符号含义同前。计算得γdc的值为1 373.85 kg/m3，计算结果接近前人的实测值，结果可信[38–39]。根据式（15）、式（16），计算得观测架附近表层沉积物临界侵蚀切应力τe为0.086 N/m2。

实验室分析结果显示，研究区表层沉积物成分为黏土（8.13%）、粉砂（46.95%）和砂（44.92%）。根据Lumborg[40]实验结果，研究区表层沉积物悬沙临界沉降剪切应力τd取0.05 N/m2。也有研究表明，τd约为临界侵蚀切应力τe的0.5 倍[41]。两种方法得到的值接近，证明结果可信。

2.3.8相关性检验

本文选用适用性较强的Spearman 秩相关系数ρ来检测水文变量之间的相关关系，这里的“秩”指的是顺序或者排序。所以，计算时先对由n个样本组成的时间序列x1，x2，…，xn，和由n个样本组成的时间序列y1，y2，…，yn进行排序（同时为升序或者降序），得到排序后的集合X和Y，Xi、Yi分别是xi、yi在X、Y中的排行，i=1，2，3…，n，则ρ的计算公式为

计算所得的相关系数ρ值，一般取小数点后两位来表示，取值范围在–1.00～1.00 之间，负数表示负相关，正数表示正相关。ρ值对应的两变量之间的相关程度为：

ρ= ±0～±0.30 时，不相关；

ρ= ±0.30～±0.50 时，微相关；

ρ= ±0.50～±0.80 时，显著相关；

ρ= ±0.80～±1.00 时，高度相关。

3 水沙动力特性与变化过程

3.1 潮流变化过程

研究区域潮型为不正规混合半日潮，只有高潮情况下才会淹没观测架，观测时段内观测架一共被淹没4 次，涨落潮历时大致相等。11 月19–22 日，分别记作T1–T4，最大水深分别为0.19 m、0.62 m、0.53 m与0.72 m（图2a）。T2–T4 潮位接近，且均显著大于T1。观测期间的潮流呈现明显的旋转流特性。涨潮阶段由约N100°顺时针旋转至约N50°，落潮阶段由约N50°顺时针旋转至约N200°（图2b）。

研究区域流速存在两种不同的随时间变化模式。T1、T2 流速呈现“单峰型”，潮周期内只出现一次流速极值，涨潮平均流速约为落潮平均流速的4 倍，最大流速出现在涨潮中期，分别为0.52 m/s 与0.46 m/s。T3、T4 流速呈现“双峰型”，涨潮、落潮各出现一次流速极值，涨落潮平均流速接近，涨落潮中期均存在明显的流速峰值。T3 涨潮最大流速为0.49 m/s，落潮最大流速为0.50 m/s；T4 涨潮最大流速为0.60 m/s，落潮最大流速为0.53 m/s（图2c）。

3.2 风浪变化过程

T1 期间波浪动力较弱，有效波高Hs低于0.05 m。T2–T4 期间波浪动力较强，有效波高Hs最大值分别为0.24 m、0.21 m 及0.27 m。T2–T4 期间，涨落潮初期流速、流向波动较大，这是因为波浪动力较强，在水深较浅的情况下，波浪对流速流向的影响较显著（图2d）。T1 期间有效波浪周期Ts在1.24～3.21 s 之间，平均为2.32 s。T2–T4 期间有效波浪周期Ts在1.16～2.11 s 之间，平均为1.68 s，高于T1 期间。

T1 期间平均风速为1.33 m/s，最低为0.5 m/s，最高为2.1 m/s，风向由向岸（328°）转变为离岸（73°）。T2–T4 期间平均风速为2.67 m/s，最高为4.60 m/s，最低为0.10 m/s，风向为离岸风，平均风向为153.94°（图2e）。

图2 观测期间水动力、波浪与风速随时间变化序列Fig. 2 Time series of hydrodynamic force, wave and wind during the observation period

3.3 床面剪切力及紊动剪切率变化过程

T1 期间，波流共同引起的床面剪切力（τcw）平均为0.09 N/m2，τcw最大达到0.34 N/m2，中潮期间τcw平均为0.23 N/m2，最大可以达到0.99 N/m2。T1 期间潮流引起的床面切应力τc平均为0.05 N/m2，最大达到0.21 N/m2；T2–T4 期间τc平均为0.14 N/m2，最大达到0.68 N/m2。观测期的4 d 中，总床面剪切力τcw与潮流引起的床面剪切力τc的变化趋势接近，峰值都出现在落潮中期（图3a）。

黄河口为弱潮河口，研究区属于开敞性海岸，波浪没有遮蔽，波浪引起的床面剪切力τw较高[42]。T1 期间τw在0.90～1.19 N/m2之间，平均为0.98 N/m2；T2–T4 期间τw在0.99～1.30 N/m2之间，平均为1.11 N/m2，T2–T4 期间τw高于T1。波浪引起的床面剪切力变化图线呈“U”形，涨潮初期与落潮后期τw较高，水深较大时的τw反而较小。说明水深较浅时，波浪对床面切应力的影响会明显增加（图3a）。

紊动剪切率G的变化规律与床面切应力大致相同，在落潮中期出现峰值。T1 期间紊动剪切率在1.88～9.63 s-1之间，平均为3.33 s-1；T2–T4 期间紊动剪切率在2.51～28.48 s-1之间，平均为9.00 s-1，T2–T4期间水体紊动剪切率大于T1（图3b）。

3.4 悬沙浓度变化过程

研究区域水体悬沙浓度较低，T1 期间水体平均悬沙浓度为0.08 kg/m3，最大可以达到0.26 kg/m3。T2–T4 期间，落潮中后期水体悬沙浓度会明显增加，出现一个明显的悬沙浓度峰。T3–T4 期间水体平均悬沙浓度为0.10 kg/m3，最大可以达到0.40 kg/m3。悬沙浓度峰出现的时间与床面切应力峰值出现的时间接近，且出现悬沙浓度峰时，床面切应力显著高于临界剪切应力τe（0.086 N/m2），说明这一悬沙浓度峰的出现源于泥沙再悬浮（图3c）。

图3 观测期间床面切应力与悬沙浓度（SSC）随时间变化序列Fig. 3 Time series of bed shear stress, bottom turbulent shear rate and suspended sediment concentration (SSC) during the observation period

T2 期间水体悬沙浓度显著高于T3 和T4，水体平均悬沙浓度为0.38 kg/m3，最大可以达到0.62 kg/m3。T2 期间，悬沙浓度随涨潮不断增加，在最高潮位前后保持峰值，后随水深降低而不断减小，在落潮中后期出现第2 个悬沙浓度峰。T2 期间两个悬沙浓度峰值出现的时间与床面切应力峰值出现的时间不一致，悬沙浓度变化过程与床面切应力变化过程不同步。但是T2 期间悬沙浓度增加的阶段，与流速峰值出现的时间基本吻合。推测T2 期间的悬沙浓度峰值并不是来源于泥沙再悬浮。可能是由于涨潮水流从附近携带大量泥沙，并在高水位期间沉降造成的[43]（图3c）。

观测中发现，研究区潮滩存在大量的“泥丸”型结构，即大型的泥沙聚集体（图4）。这种聚集体呈球状，可以随潮汐运动在滩面迁移滚动，但又容易破坏，一旦分散会导致所在局部悬沙浓度大幅上升。T2 期间出现悬沙浓度峰值，可能是发生了“泥丸”的迁移与分散。

图4 研究区域的“泥丸”型泥沙聚合体Fig. 4 The “pill-like” sediment aggregates in the study area

4 絮凝特性与变化

4.1 絮团粒径

T1 期间絮团粒径先增大后减小，在最高潮位阶段达到最大，为264.44 μm。T2–T4 期间絮团粒径呈现缓慢递减趋势，中潮超过60%的时间内，絮团中值粒径都低于100 μm。但是在T2–T4 每个落潮的中后期絮团中值粒径会出现一个峰值，最大可达到203.83 μm。小潮期间絮团平均粒径为123.90 μm，大于中潮期间的91.81 μm（图5a）。

4.2 絮团体积浓度

T1 期间絮团总体积浓度平均为221.05 μL/L，最大为810.22 μL/L；T2–T4 期间絮团总体积浓度平均为390.77 μL/L，最大为1 100.29 μL/L。在T1 最高潮位及T2、T3 的落潮中期，絮团总体积浓度均出现一个峰值（图5b）。

图5 观测期间絮团粒径、体积浓度、有效密度与沉速随时间变化序列Fig. 5 Time series of floc size, volume concentration, effective density and settling velocity during the observation period

4.3 絮团有效密度与沉速

T1 期间黄河口潮滩絮团有效密度在56.55～452.44 kg/m3之间，平均为257.03 kg/m3。T2–T4 期间絮团有效密度逐日递减，T2 期间平均有效密度为427.77 kg/m3，T3 期 间 絮 团 平 均 密 度 为319.38 kg/m3，T4 期间平均有效密度为166.32 kg/m3（图5c）。T1 期间，絮团沉速在0.99～1.86 mm/s 之间，平均为1.41 mm/s。T2–T4 期间，絮团沉速在0.17～4.07 mm/s之间，平均为0.94 mm/s。总体来看，中潮期T2–T4 期间絮团的沉速大于小潮期T1（图5d）。

5 分析与讨论

5.1 絮团粒径变化的影响因素

从泥沙运动学的角度，天然水体中的絮凝主要受到悬浮泥沙粒径、水动力环境、悬沙浓度的控制[17]。本研究期间受交通与天气条件限制，未采集悬沙样品，无法分析悬浮泥沙粒径对絮凝的影响，但可分析水动力与悬沙浓度与絮凝的相关性。采用相关性分析方法，分别计算各个潮周期内水体紊动剪切、流速及悬沙浓度与絮团粒径之间的Spearman 系数，计算结果及对应的相关性如表1。絮团粒径与流速的关系主要为不相关，与悬沙浓度及底部紊动剪切率则存在着较明显的相关关系，说明黄河口潮滩絮团粒径主要受到底部紊动剪切率与悬沙浓度的影响，与流速关系不大。

表1 各个潮周期絮团粒径与紊动剪切率、悬沙浓度及流速的Spearman 系数Table 1 Spearman correlation coefficient of floc particle size and bottom turbulent shear rate, suspended sediment concentration and flow velocity

5.2 水体紊动剪切促进絮凝的上限

研究发现，不同潮周期中，黄河口潮滩絮团粒径与紊动剪切率的关系不尽相同。T1 期间，絮团粒径与紊动剪切率变化趋势基本一致，粒径随紊动剪切率增大而增大，紊动剪切率减小时粒径也随之减小。但是T2–T4 各个潮周期内，水体紊动剪切率与絮团粒径呈现负相关趋势，紊动强度较大时对应粒径较小。

在影响絮凝的外部因素中，水体紊动是最重要的因素之一，水体紊动强度的增大可以增加颗粒的碰撞频率，促进颗粒聚合；但当紊动强度增大到一定程度时，反而会增加絮团破碎频率抑制絮凝[2,44]。T1 期间为小潮，水动力较弱，底部紊动剪切率平均仅为3.33 s–1。T2–T4 期间底部紊动剪切率较高，平均为9.00 s–1，最高可达到28.48 s–1。在T1 期间紊动剪切较弱的情况下，水体紊动增加对絮凝起促进作用；T2–T4 期间紊动强度较高，水体紊动增加抑制絮凝。

为了进一步探究黄河口潮滩水体紊动影响絮凝的上限，首先对实测数据进行等距分箱以体现规律性，分箱间距为4[22–23]。发现絮团粒径随紊动强度增加先增后减。选择不同的临界紊动剪切率值（Gl），将观测结果分为“絮团粒径随紊动强度增加而增大”与“絮团粒径随紊动强度增加而减小”阶段，并对两个阶段分别使用幂函数拟合（图6）。

图6 絮团粒径与底部紊动剪切率关系Fig. 6 Variation of floc size with bottom turbulent shear rate

结果发现，当Gl=3.76 s-1时，取得的拟合结果最为理想（表2）。T1 期间，水体紊动剪切率低于3.76 s-1，T2–T4 期间则高于这一值，统计分析得到的结果与实测数据能够相互验证。因此可以认为，黄河口当地水体紊动促进絮凝的上限为Gl=3.76 s-1，水体紊动强度低于这一数值时水体紊动促进絮凝，高于这一数值时水体紊动破坏絮凝。

表2 不同Gl 对应拟合结果对比Table 2 Fitting results of different Gl

5.3 悬沙浓度对絮凝影响

悬沙浓度对絮凝的影响尚未有确切定论，一些实验发现，悬沙浓度增加增大了颗粒碰撞概率，对絮凝起促进作用[2,20,45]；但也有研究认为，悬沙浓度增加导致絮团碰撞更频繁，造成絮团破坏[17,20,45–46]。参考郭超[17]的研究方法，对黄河口潮滩絮团粒径与悬沙浓度之间的关系整体采用幂函数拟合，发现二者呈负相关关系，相关系数为R²=0.778 8。说明在黄河口潮滩，悬沙浓度增加对絮凝起抑制作用（图7）。

图7 絮团粒径与悬沙浓度（SSC）关系Fig. 7 Variation of floc size with suspended sediment concentration

在5.1 节的讨论中，发现T2–T4 期间悬沙浓度与絮团粒径呈负相关，但是T1 期间悬沙浓度与絮团粒径没有明显相关性。这是因为絮凝同时受到悬沙浓度与水动力共同控制，T1 期间水体紊动强度先增后减，悬沙含量则缓慢降低。絮团粒径的变化趋势与水体紊动强度类似，也是先增后减。这说明相比于悬沙含量变化造成的影响，水体紊动剪切对絮凝的影响更大。悬沙含量降低对絮凝的促进作用，不足以抵偿水体紊动强度变化对絮凝的影响，水体紊动剪切才是絮团粒径变化的主导因素。因此在T1 期间，絮团粒径与水体紊动剪切强度显著相关，但是与悬沙浓度几乎没有相关性。

5.4 絮团有效密度与沉速变化规律

黄河口潮滩絮团粒径与有效密度呈负相关关系，粒径更大的絮团对应的有效密度更小。过往研究中，通常使用 Δp=·Dam的关系来描述絮团粒径与有效密度的关系，a值在–0.46（Sternberg 等[6]在美国加利福尼亚海域测得）到–1.58（Xia 等[12]在中国珠江口测得）之间[5,10,19]。郭超[17]、Guo 等[21]在荷兰Kapellebank 潮滩常态天气下测得a=–0.53，风暴天气下测得a=–1.28。

分别对黄河口潮滩小潮（T1）与中潮（T2–T4）期间的絮凝特性进行研究，计算得小潮期间a=–1.52，中潮期间a=–0.79，属于上述文献中所列的取值范围（图8）。无论在何种动力条件下，有效密度都与絮团粒径成反比，粒径越大的絮团对应的有效密度越低。在水体紊动剪切促进絮凝的情况下，泥沙聚合体在由小变大的过程中，先由易絮凝的细颗粒泥沙发生絮凝，而后捕获粒度相对更大的泥沙，形成的絮凝体越高阶，其孔隙率越大，有效密度越低[10]。在水体紊动剪切抑制絮凝的情况下，结构疏松、有效密度低的大絮团更容易受到紊动影响破碎，形成粒径更小但是有效密度更高的小泥沙聚合体[10,17,19]。

图8 絮团粒径与有效密度关系Fig. 8 Variation of floc size with effective density

由于絮凝受到泥沙特性、当地生物化学特性等多方面影响，a值与水动力强弱的关系目前尚无定论。Guo 等[21]在荷兰Kapellebank 潮滩的研究中，水动力较强的风暴天气下计算得的a值低于水动力相对弱的常态天气；但是本研究与其他观测结果中，也发现了相反的结论[9,16]。本研究在黄河口观测的结果显示，小潮期间a值低于中潮；程江等[9]在长江口徐六泾的观测结果也类似，小潮期间的a值（–0.646）低于大潮期间（–0.368）。

根据式（2），絮团的沉速与絮团粒径及有效密度有关。过往研究发现，沉速与粒径之间的关系与“有效密度-粒径”关系类似，可以用 ωs=c·Ddm的形式来表示。c取值在0.290～1.632 之间；d取值在0.000 022 8～0.008 20 之 间[5–6,10,12]。黄 河 口 潮 滩 絮 团沉速与粒径之间相关性较好，沉速与粒径间关系可以用ωs=0.0028·D1m.2679表示，属于过往文献中所列取值范围（图9）。但是对有效密度与沉速关系进行分析时发现，二者之间不存在明显相关性。说明黄河口潮滩絮团沉速主要受到絮团粒径的影响，与有效密度的关系不明显。

图9 絮团粒径与沉速关系Fig. 9 Variation of floc size with settling velocity

6 结论

本文通过研究黄河口南部潮滩絮团特征的变化，得到如下几点认识：

（1）黄河口潮滩絮团粒径在25.42～264.44 μm 之间，平均为95.20 μm。小潮潮周期T1 期间絮团总体积浓度平均为221.05 μL/L，平均有效密度为257.03 kg/m3，平均沉速为1.41 mm/s；中潮潮周期T2–T4期间絮团总体积浓度平均为390.77μL/L，平均有效密度为304.49 kg/m3，平均沉速为0.94 mm/s。

（2）水体紊动较低时对黄河口潮滩絮凝起促进作用，较高时则抑制絮凝。随紊动强度增加，黄河口潮滩絮团粒径先增后减。分析得底部紊动剪切率促进絮凝的上限值为Gl=3.76 s-1。

（3）悬沙浓度增加对黄河口潮滩泥沙絮凝起抑制作用，但是悬沙浓度变化对絮凝的影响低于水体紊动剪切变化对絮凝的影响。

（4）黄河口潮滩絮团有效密度与粒径呈现Δp=b·Dam的关系。T1 潮周期a值为–1.52，T2–T4 潮周期a值为–0.79。

（5）黄河口潮滩絮团沉速主要受到泥沙粒径影响。絮团沉速与粒径之间相关性较好，但与有效密度之间不存在明显相关关系。

致谢：本文野外观测工作得到水利部黄河水利委员会黄河水利科学研究院黄河河口研究所、山东黄河三角洲国家级自然保护区等单位的大力支持，在此表示诚挚的感谢！