本次介绍的研究题为《Threshold of surface erosion of cohesive sediments》（《粘性沉积物表面侵蚀的阈值》），由Dake Chen、Jinhai Zheng、Chi Zhang*、Dawei Guan、Yuan Li及Huiming Huang共同完成。研究者来自南京水利科学研究院、河海大学港口海岸与近海工程学院以及水文水资源与水利工程科学国家重点实验室。该研究作为一篇原创性研究（Original Research）文章，于2022年9月15日发表在学术期刊《Frontiers in Marine Science》上。

一、 研究背景 该研究隶属于海岸与海洋工程及沉积动力学领域，核心科学问题是粘性沉积物在水动力作用下的侵蚀起始条件。粘性沉积物主要由粘土和粉砂颗粒组成，普遍存在于河口、海岸及湖泊等水生系统中。它们不仅是底栖生物的重要栖息地，也是有机碳储存和生物地球化学循环的关键场所，其输运过程深刻影响着水质、水生生态系统及地貌演变。

侵蚀是沉积物动力学的控制性过程之一，而侵蚀阈值则是描述侵蚀起始临界水动力条件的核心参数，对于精确模拟粘性沉积物侵蚀过程、解决相关工程与环境问题至关重要。然而，与无粘性泥沙（如砂）相对成熟的希尔兹曲线（Shields diagram）理论相比，粘性沉积物的侵蚀阈值预测仍面临巨大挑战。主要难点在于其颗粒间存在复杂的电化学引力（内聚力），并且这些细颗粒在水中通常会凝聚成结构松散的絮团（aggregates）或集合体，形成具有分形（fractal）特征的结构。现有研究多采用经验公式，将侵蚀阈值与干容重、含水量、塑性指数等物性参数关联，但这类公式通用性差，且物理机制不清晰。另一种基于力学平衡的解析方法虽能揭示物理本质，但如何量化颗粒间的内聚力，尤其是如何将沉积物床面由分形絮团构成这一关键结构特征纳入理论模型，尚未得到很好解决。因此，本研究旨在建立一个能够考虑分形絮团结构的粘性沉积物表面侵蚀阈值理论公式，并量化絮团有效自重对侵蚀阈值的贡献率，以弥补当前研究的不足。

二、 研究详细流程 本研究采用理论推导与实验数据验证相结合的方法，具体流程可分为以下几个步骤：

第一步：理论模型构建与公式推导。 这是研究的核心创新部分。研究团队首先从粘性沉积物床面的微观结构出发，引入分形理论来描述絮团。分形维数（fractal dimension, *f*）是描述絮团空间填充能力的关键参数，其值介于1（松散线状）到3（致密实体）之间。利用分形关系，可以推导出絮团大小（*d_a*）和密度（*ρ_a*）与原始颗粒直径（*d_p*）、固体体积分数（*ϕ_s*，反映沉积物压实程度）以及分形维数（*f*）之间的数学关系。

接着，研究者量化了颗粒间的内聚力。他们摒弃了传统的基于经验幂函数的做法，而是选择了物理机制更清晰的范德华力（van der Waals force）作为内聚力的基本来源。通过一系列几何和统计假设，他们将平均范德华力与沉积物的固体体积分数联系起来，从而建立了一个依赖于沉积物压实度的内聚力表达式。

然后，研究对位于水-床界面的一个代表性絮团进行受力与力矩平衡分析。如图1所示，絮团受到水流拖曳力（*F_D*）和上举力（*F_L*）的作用，同时受到其水下有效自重（*G_a*）以及周围颗粒通过范德华力施加的合力（*F_c*，抵抗侵蚀）的作用。通过分析絮团绕支点转动的临界力矩平衡条件，并代入前述的絮团分形参数和内聚力表达式，最终推导出一个预测表面侵蚀临界剪切应力（*τ_cr*）的理论公式（即文中的公式19）。该公式的关键形式为： q_cr = q_cr0(d_a) * [ (ρ_s - ρ_w)g d_p * (第一项：代表絮团有效自重的贡献) + (C / (ρ_s - ρ_w)g d_p^2) * (第二项：代表内聚力的贡献) ]* 其中，q_cr 为临界希尔兹参数，q_cr0(d_a)* 是对应絮团尺寸的无粘性泥沙临界希尔兹参数，C 是一个综合反映沉积物内聚强度的系数（与哈马克常数、结合水层厚度等因素相关）。该公式清晰地表明，粘性沉积物的侵蚀阈值来源于两部分：絮团（被视为一个整体）的水下有效自重提供的阻力，以及颗粒间内聚力提供的阻力。

第二步：模型应用与参数确定——两种方法。 为了验证和应用推导出的公式，研究者从已发表的文献中系统收集了四类粘性沉积物的侵蚀阈值实验数据，包括三组高岭土（kaolinite）、两组石英粉（quartz）、九组湖塘淤泥（lake and pond mud）以及九组海岸淤泥（coastal mud），共计23个数据集。每个数据集都包含了同一类沉积物在不同固体体积分数（或容重）下对应的临界剪切应力测量值。

研究采用两种方法将理论公式与实验数据拟合，以确定公式中的关键参数——分形维数 f 和内聚系数 *C*。 * 方法一（常量 *f*）： 将分形维数 f 视为一个不随沉积物状态变化的常数。通过非线性回归，为每个数据集找到一组最优的 (f, C) 值，使得公式计算值与实验值吻合最佳。结果显示，虽然使用常数值 f 也能获得较好的拟合效果（图2-4），但这种方法获得的 f 值物理意义不明确，且不足以充分约束模型。 * 方法二（变量 *f*）： 基于分形几何和前人研究，假设分形维数 f 是絮团尺寸与原始颗粒尺寸之比的函数（公式24），进而推导出 f 与固体体积分数 ϕ_s 和原始颗粒直径 d_p 相关的具体表达式（公式25）。其中引入了一个系数 *b*。在此框架下，再次通过非线性回归为每个数据集确定最优的 (b, C) 值。结果显示，这种方法不仅拟合效果良好，而且获得的 f 值范围（2.01~2.80）与文献中报道的粘性沉积物床面分形维数范围一致，更具物理合理性。

第三步：数据分析与规律提炼。 基于方法二获得的结果，研究者进行了深入分析： 1. 分形维数 (f) 的规律： 他们发现系数 b 的值随原始颗粒直径 (d_p) 增大而线性减小（图5，公式26）。结合公式25和26，可以计算出 f 随 ϕ_s 增大而增大（压实度越高，结构越致密），随 d_p 增大而减小（图6）。这一发现量化了分形维数的影响因素。 2. 内聚系数 (C) 的规律： C 值具有明显的“场地或沉积物特异性”。分析显示，高岭土和海岸淤泥的 C 值量级相近（约10^-5 J/m²），石英粉的 C 值低1-2个量级（约10^-6 J/m²），而富含有机质（11~19%）的湖塘淤泥的 C 值则高出1-2个量级（约10^-3 J/m²）。这反映了不同矿物成分和有机质含量对内聚强度的显著影响。基于最小化预测误差的原则，研究为每类沉积物提供了 C 的参考值（guideline value），便于实际应用。 3. 絮团特性与贡献率分析： 利用确定的分形维数公式，可以估算床面絮团的平均直径（主要在几十到几百微米量级）以及单个絮团包含的原始颗粒数。更重要的是，研究首次量化了絮团有效自重对侵蚀阈值的贡献率。计算表明（图9）：对于原始颗粒细小、固体体积分数高的沉积物（如 d_p = 0.004 mm, ϕ_s = 0.25-0.35），絮团自重的贡献率很低（仅0.6~2.1%），侵蚀阻力主要来自内聚力；相反，对于颗粒较粗、压实度较低的沉积物（如 d_p = 0.02 mm, ϕ_s = 0.10-0.15），絮团自重的贡献率可达33~53%，不可忽略。研究据此提出了一个临界贡献率（30%）及相应的判别条件（公式28）。 4. 公式简化： 对于自重贡献率低于30%的情况（即内聚力主导），理论公式可简化为一个仅含一个经验系数 A 的更简洁形式（公式30），该公式同时适用于粘性细颗粒和非粘性粗颗粒，展示了良好的普适性（图11）。

三、 主要研究结果 1. 成功建立了粘性沉积物表面侵蚀阈值的理论预测公式。 该公式（公式19）首次系统地将沉积物床面的分形絮团结构、基于范德华力的内聚力量化以及絮团受力力矩平衡分析相结合，从机理上阐明了侵蚀阈值的组成。 2. 揭示了分形维数的变化规律。 提出了分形维数 f 与原始颗粒直径 d_p 和固体体积分数 ϕ_s 的函数关系（公式25, 26），指出 f 随 ϕ_s 增加而增加，随 d_p 增加而减小。 3. 量化了内聚强度系数 C 并提供了参考值。 明确了 C 的场地特异性，并基于大量数据为高岭土、石英、含有机质淤泥和无明显有机质淤泥分别给出了具有实用价值的 C 参考值。 4. 首次精确量化了絮团有效自重对侵蚀阈值的贡献率。 纠正了以往认为粘性沉积物侵蚀阻力主要或完全来自内聚力的片面认识，明确指出在特定条件（颗粒较粗、床面较松散）下，絮团自重的贡献十分显著。并给出了判断自重贡献是否可忽略的定量标准（贡献率30%）。 5. 完成了广泛的数据验证。 理论公式及其衍生的参数化模型成功应用于来自不同环境、具有不同物性的23组实验数据，预测结果与实测值吻合良好（图2-4, 7, 10, 11），证明了模型的有效性和鲁棒性。 6. 提出了适用于不同情况的简化公式。 在自重贡献可忽略的条件下，推导出了形式简洁的单参数公式（公式30），便于工程应用，且能统一描述从粘性到非粘性泥沙的侵蚀阈值趋势。

四、 研究结论与价值 本研究通过创新的理论建模和系统的数据验证，成功发展了一套用于预测粘性沉积物（分形絮团构成）表面侵蚀阈值的完整理论框架。该框架不仅提供了具有明确物理意义的预测公式，还揭示了分形维数的内在变化规律、内聚强度的材料依赖性，并首次定量厘清了絮团自重与内聚力在抵抗侵蚀中的相对重要性。

科学价值： 研究将分形理论、胶体化学（范德华力）与沉积物侵蚀力学紧密结合，深化了对粘性沉积物微观结构-宏观力学性质之间关联的理解，推动了粘性沉积动力学从经验描述向机理预测的发展。提出的理论模型为后续研究床面絮团的侵蚀、输运和再沉降全过程提供了重要的理论基础。

应用价值： 研究所提供的公式及参数参考值，可直接用于改进海岸、河口、水库等环境中粘性沉积物侵蚀和输运的数值模型，提高对地貌演变（如岸滩侵蚀、航道淤积）、水环境质量（如污染物再悬浮）以及水生生态系统影响的预测精度。简化公式（公式30）尤其便于工程师在涉及粘性泥沙问题的项目中快速估算侵蚀阈值。

五、 研究亮点 1. 理论创新性强： 首次在粘性沉积物侵蚀阈值理论中综合考虑了分形絮团结构和范德华力内聚力，推导过程严谨，物理图像清晰。 2. 量化分析深入： 不仅建立了预测模型，更深入量化了模型关键参数（f, *C*）的变化规律及其物理含义，特别是首次精确量化了絮团自重的贡献率，解决了该领域的一个长期疑问。 3. 验证数据广泛： 使用了涵盖人工材料（高岭土、石英）和天然沉积物（湖泥、海泥）的多元数据集进行验证，增强了结论的可靠性和普适性。 4. 成果实用性好： 研究最终提供了包含完整参数化方案的理论模型和便于使用的简化公式，并给出了关键系数的参考值，架起了理论与工程应用之间的桥梁。

六、 其他有价值内容 研究还讨论了“结合水层重叠产生的附加压力”作为另一种“假内聚力”的观点，但因在大多数实验水深较浅的工况下其影响较小，且与电化学内聚力的相对大小尚未量化，故未纳入本次模型。这为未来的研究留下了一个可探索的方向。此外，研究承认有机质等产生的粘附力（adhesion）未予考虑，富含有机质的湖泥表现出异常高的 C 值也印证了这一点，提示在有机质丰富的环境中，需要进一步扩展模型以包含粘附力机制。