黄河流域黏性泥沙群体沉速计算分析

DOI:10.12170/20221013002
丰青，肖千璐，郑艳爽，等. 黄河流域黏性泥沙群体沉速计算分析[J]. 水利水运工程学报，2024（1）：8-14. （FENG Qing, XIAO Qianlu, ZHENG Yanshuang, et al. Examination of cohesive sediment settling velocity calculations in the Yellow River Basin[J].Hydro-Science and Engineering, 2024（1）: 8-14. （in Chinese ））
黄河流域黏性泥沙群体沉速计算分析
丰 青1, 2，肖千璐2，郑艳爽2，张晓华2
(1. 北京大学 环境科学与工程学院，北京 100084； 2. 黄河水利委员会黄河水利科学研究院，河南 郑州 450003)
摘要: 针对黄河流域由黏性矿物组成的黏性泥沙群体沉降特征展开研究。

基于泥沙颗粒受力分析和能量守恒
定律，以处于稳态的泥沙颗粒为参照物，构建了泥沙颗粒水沙界面能量耗散物理图形；建立了考虑挟沙水流运动
黏性系数和紊动黏性系数的泥沙群体沉速显式表达式，并采用黄河流域典型黏性泥沙红色砒砂岩颗粒开展了泥
沙群体沉速试验研究。

结果表明：红色砒砂岩组成的含沙水体中泥沙群体沉速整体随含沙量增大而减小并逐渐
趋于定值；当垂向平均含沙量为0.34~0.45 kg/m 3时，泥沙群体沉速表现为随含沙量的增加而增大。

采用试验数
据对建立的公式进行验证，公式能较合理地表达水流黏性和尾流分离对泥沙颗粒沉降阻力的影响。

研究结果可
为黄河流域黏性泥沙淤积和输移规律研究提供参考。

关　键　词：黏性泥沙； 沉降速度； 能量耗散； 黏性系数； 黄河流域
中图分类号：TV142.1 文献标志码：A 文章编号：1009-640X(2024)01-0008-07
泥沙颗粒沉降特性是泥沙运动力学研究的关键基础问题之一。

与单颗粒泥沙沉降不同，泥沙颗粒群体沉降特性受到泥沙级配组成、水流含沙量及颗粒相互作用等的影响[1]。

黄河流域跨越多种地貌单元和气候带，黄河泥沙具有来源广、沙量多的特征。

根据黄河泥沙不同来源区颗粒级配分析结果可知，泥沙颗粒级配由小于0.005 mm 到1.000 mm 连续存在，河道中不同颗粒组成泥沙的水沙动力特征差异较大[2-3]。

王兆印等[4]对黄河干支流泥沙进行野外采样分析，认为黄河流域泥沙中矿物成分及含量是影响水体中悬浮颗粒流变特性的主导因素，黏性矿物含量对高含沙水流流变特性及沉速的影响比粒径分布的影响更大；另一方面，泥沙在河道中发生沿程分选，挟沙水流中悬浮泥沙沉降表现出显著的非均匀泥沙沉降特征，且不同含沙量及不同矿物组成条件下泥沙沉降规律可能发生变化[5]。

这表明仅考虑泥沙粒径变化不能准确反映黄河流域泥沙的动力学效应，还应充分考虑泥沙矿物组成尤其是黏性沙矿物组成及含沙量的影响。

水体絮凝结构的形成与水体内泥沙颗粒pH 值和盐度直接相关[6]，考虑到黄河流域干支流非感潮河段的水体和土壤基本呈碱性，在含沙量不是特别高的条件下基本无法形成絮凝态结构[7]。

因此，在进行黄河流域泥沙颗粒群体沉速研究过程中，需主要考虑挟沙水流黏性特征及含沙量的影响。

S v 国内外学者采用试验研究和理论推导等方法围绕挟沙水流中含沙量对泥沙颗粒群体沉速影响开展了深入研究，提出了不同类型的泥沙颗粒群体沉速计算公式[8-9]。

沙玉清[10]以体积含沙量=0.3为临界值，考虑不同流区阻力系数和雷诺数关系变化，建立了泥沙颗粒群体沉速分段表达式；窦国仁[11]引入泥沙颗粒雷诺数对泥沙颗粒尾流产生的形状阻力进行分析，推导了含有泥沙颗粒沉速的尾流形状阻力表达式，采用该式求解沉速需要进行较复杂的迭代试算；郜国明等[12-13]采用深度取样锥形筒流态化法开展了非均匀泥沙群收稿日期：2022-10-13
基金项目：国家自然科学基金资助项目（51909100，52009047）；中央级公益性科研院所基本科研业务费专项项目（HKY-
JBYW-2020-03，HKY- JBYW-2019-07）
作者简介： 丰　青（1986—），男，河南郑州人，高级工程师，博士，主要从事河流动力学及泥沙环境学研究。

E-mail ：***************通信作者： 肖千璐（E-mail ：****************）
第 1 期
水　利　水　运　工　程　学　报No. 12024 年 2 月HYDRO-SCIENCE AND ENGINEERING Feb. 2024
体沉降速度试验研究，建立了包含悬沙浓度、中值粒径及非均匀性系数的泥沙群体沉速表达式。

上述研究多适用于非黏性泥沙，而黄河流域泥沙来源区中的泥沙多含黏性矿物成分，特别是黄河中游砒砂岩区的红色砒砂岩富含蒙脱石，具有较强的水溶性和膨胀性，泥沙颗粒群体沉降规律与非黏性沙明显不同。

通过野外采样分析可知，该类泥沙中值粒径小于0.05 mm ，但在实验室静态水体沉降过程中并未形成明显的絮凝结构。

在天然暴雨产流条件下，红色砒砂岩溶于水后会明显提升水体黏滞性，造成泥沙沉降特性及含沙水流特别是高含沙水流输送特性发生显著变化，低含沙水流条件下的泥沙群体沉速计算方法不再适用，需进一步考虑挟沙水流黏性和紊动特征变化对泥沙颗粒沉降过程中能量耗散的影响。

本文从泥沙颗粒沉降过程的受力分析出发，考虑泥沙颗粒沉降引起的运动黏性切应力和紊动黏性切应力的影响，引入水流综合黏性系数，建立反映水流黏性及颗粒阻力特征影响的泥沙沉降速率表达式，并采用黄河流域典型黏性泥沙红色砒砂岩颗粒开展泥沙群体沉速试验对构建公式进行验证分析，可为黄河流域黏性泥沙淤积和输移规律研究提供参考。

1 泥沙颗粒沉速理论分析
针对挟沙水流中含沙量对水流黏性及泥沙沉降模式的影响，张红武等[14]基于水沙两相连续性假定提出在固定单位水柱内，泥沙颗粒沉降将引起同体积水体向上运动，根据质量守恒定律，有：
S v ω=(1−S v )(V up −ω)(1)
ωV up 式中：S v 为体积含沙量；为泥沙群体沉速；为泥沙颗粒运动速度。

根据式（1）可得：
ω=V up (1−S v )(2)
V up S v =0ω=V up =ω0ω0当为定值，假设时，有，其中为清水中的泥沙颗粒沉速。

考虑含沙水体黏性变化，张红武等[14]提出泥沙群体沉速计算公式如下：
ω=(ν0/νm )n 3V up (1−S v )(3)
ν0νm n 3n 3=3.18式中：为清水运动黏性系数；为浑水运动黏性系数；为系数。

张红武等[14]采用实测数据对上式进行率定计算，推荐采用。

V up 通过欧拉和拉格朗日不同参照物对比分析和泥沙颗粒的受力分析，泥沙颗粒在下沉过程中，其下部水体和侧面水体不断补入泥沙颗粒原有位置。

试验结果表明，泥沙颗粒下落过程中其尾部可能出现负压区域但未出现真空区域。

考虑到泥沙颗粒具有双电层结构，下落过程中补充进入原位置处的水体密度仍为自由水密度，并非束缚水内的高密度状态。

以水体内固定位置处为参照物，根据动量守恒定理，也能获得式（1）形式的表达，但经过泥沙颗粒运动状态及受力特征分析可知，速度并非为定值，即使在清水条件下，泥沙颗粒粒径差异和泥沙颗粒沉降过程引起的尾流强度变化导致置换水体运动速度有所不同。

若以泥沙颗粒及其置换水体为参照物，当泥沙颗粒下沉过程中进入稳态，即泥沙颗粒水下有效重力等于阻力时，根据能量守恒定律，单位长度内泥沙颗粒损耗的势能除转化为泥沙颗粒的动能外，泥沙颗粒下沉速度与置换水体之间的速度差被泥沙颗粒周围的黏性切应力摩擦消耗。

若以水面或河床床面为参照物，泥沙颗粒沉降过程的能量耗散可采用下式表达：12(π6D 3)ρ(ω0−V up )2=−πD 2ρνΣ∂V up ∂n (4)
D ρνΣνm 式中：为混合沙代表粒径，在进行泥沙群体沉速计算中一般取中值粒径；为水体密度；为水流复合黏性系数，为运动黏性系数和紊动黏性系数的耦合值，在泥沙颗粒沉降过程中未出现边界层分离诱发的紊流阻力情况下可取运动黏滞系数；当泥沙颗粒出现边界层分离时，清水条件下单颗粒泥沙可通过窦国仁采用
第 1 期丰 青，等：黄河流域黏性泥沙群体沉速计算分析9
ρνΣ∂V up /∂n ∂V up /∂n 的分离角试验结果确定水流复合黏滞系数的阈值[15]；而浑水水流泥沙颗粒之间由于具有较强的相互作用，复合黏滞系数不易通过单颗粒沉降分析获取，大多采用室内试验率定的方式确定其数值。

实际上，为作用于泥沙颗粒表层的黏性切应力，流速梯度方向经过泥沙颗粒的重心。

根据泥沙颗粒侧向摩擦的物理图形，可将进一步表示为：
∂V up ∂n =∂V up ∂D ∂D ∂n =D ∂V up ∂D (5)
将式（5）代入式（4）并积分，有：
V up =ω0−12νΣ/(D +c 2)(6)
c 2式中：为积分常数。

泥沙颗粒在水体沉降过程中，考虑最小的泥沙颗粒不能小于泥沙颗粒的薄膜水厚度，即极限情况下的厚度值为薄膜水厚度，此时有：
V up =ω0−12νΣ/(D +δ)(7)
δ式中：为泥沙颗粒薄膜水厚度。

联立式（3）与式（7），则可获得考虑泥沙颗粒沉降引起的黏性切应力变化的沉降速度表达式：
ω=(ν0νm )n 3(ω0−12νΣD +δ)(1−S v )(8)
νm δ=4×10−7m νm 式（8）中含沙水流黏性系数可采用费祥俊[16]提出的黏性系数公式进行计算，泥沙颗粒薄膜水厚度采用韩其为[17]推荐的数值进行计算。

浑水运动黏滞系数的计算式为：νm =ν0[1−k s S v S vm ],S vm =0.92−0.2lg ∑P i D i ,k s =1+2(S v S vm
)0.3(1−S v S vm )4(9)S vm P i S vm ∑P i /D i 式中：为极限体积比含沙量；和分别为某一粒径级的平均直径及其相应的质量百分比。

假设各粒
径级泥沙密度差异不大，则质量百分比约等于体积百分比，式中
可理解为混合沙单位质量泥沙比表
面积的一种表达形式。

2 泥沙颗粒沉速试验
试验装置采用定制的搅拌式沉降桶，沉降桶高2.0 m 、内径1.0 m ，搅拌桶底部安装三叶倾斜剪切桨进行含沙水体的掺混搅拌，外部由变齿变频电机提供搅拌动力。

搅拌桶上部设置加筋环，以增强搅拌桶的侧向稳定性，上部固定十字测架，在测架上焊接挂架并将测量设备固定于挂架底端（图1），测量设备可沿径向和垂向移动。

试验过程中，根据各测点空间位置操作挂架沿滑道完成指定距离的移动，并进行对应时间的含沙量测量。

丝杠升降控制器
测架与玻璃之间
的固定接口，
可采用鱼嘴夹
或其他固定扣
升降控制器距离玻
璃约 10 cm 加粗丝杠
直径: 与玻璃桶
直径相当
钢制十字固定架 图 1 搅拌筒上置十字测架和仪器挂架
Fig. 1 The configuration of cross measuring frame and instrument pylons
10水 利 水 运 工 程 学 报2024 年 2 月
泥沙颗粒群体沉速试验采用Campbell Scientific 公司生产的OBS 3+浊度仪和CR1000型数采控制器测量含沙量，通过控制挂架的伸缩长度完成系列点位的含沙量测量。

正式试验前开展了不同含沙量条件下OBS 3+浊度仪标定试验，结果如图2所示。

600300900 1 20010
02040
30OBS 输出浊度/NTU
(a) 非黏性泥沙含沙量/(k g /m 3)200100300100204030OBS 输出浊度/NTU (b) 黏性泥沙含沙量/(k g /m 3)
图 2 非黏性泥沙（花园口沙样）和黏性泥沙（红色砒砂岩）含沙水流浊度值与含沙量标定Fig. 2 Calibration diagram for OBS 3+ probe in non-cohesive sediments (sampled at Huayuankou) and cohesive sediments (red Pisha sandstone) within sediment-laden flow
试验中黏性泥沙采用黄土高原鄂尔多斯南部乌兰木伦河上中游滩岸处的红色砒砂岩，并同时取黄河下游花园口河段汛后滩面的非黏性泥沙作为对比研究。

测定黏性泥沙中值粒径为0.044 mm ，沙样涵盖各粒径的泥沙颗粒，级配曲线连续性明显且较为平缓，具有显著的非均匀特征；采用X-射线荧光光谱仪测定黏性泥沙矿物成分，结果表明蒙脱石质量占比接近50%，其他成分为石英、斜长石、钾长石等非黏性矿物。

考虑到河道中水流含沙量特征及浊度仪测量含沙量的有效量程，试验中沉降桶含沙水体各组次垂向平均含沙量分别设置为2.0、6.0、12.0、16.0 kg/m 3。

试验时先将清水注入搅拌桶直至液面距搅拌桶底高度1.2 m ，再加入提前称重的试验用沙，启动剪切桨片充分搅拌，确保桶内水沙混合均匀。

随后关闭搅拌机，在桶内水体旋转基本停止后开始测量桶内含沙量垂向分布，测点垂向间距0.1 m ，测量时间间隔3~5 min 。

每组次试验开始前先进行搅拌桶清理，以降低不同类型泥沙混合引起的试验误差。

根据各试验组次不同时刻含沙量垂向分布测量结果，对泥沙群体沉速变化规律进行分析。

为反映不同含沙量条件下沉降桶内泥沙沉降规律，以平均含沙量为2.0 kg/m 3和12.0 kg/m 3组次为例，分析含沙量变化对不同时刻沉降桶内含沙量垂向分布及沉降特性的影响。

由图3可见，在测量开始阶段，垂向含沙量差异较大，且含沙量较低时垂向分布更不均匀；随着时间的推移，含沙量变化速率逐渐减缓，低含沙量条件下浓度减缓速率大于高含沙条件，泥沙颗粒群体沉速逐渐趋于稳定值。

根据试验测得的非黏性泥沙和黏性泥沙不同时刻含沙量垂向分布数据，采用麦克劳林公式计算不同沉降时间内不同水深处泥沙群体瞬时沉速，断面平均沉速取含沙量达到初始含沙量50%时的沉降时间内的沉速平均值。

以黏性泥沙平均含沙量为16.0 kg/m 3的组次为例，计算不同水深处沉速和断面平均沉速，结果见表1。

若以最底层泥沙群体沉速作为该试验组次的泥沙群体沉速，计算表明含沙量为16.0 kg/m 3条件下泥沙群体沉速为0.014 6 cm/s 。

由试验结果可知，红色砒砂岩组成的含沙水体中泥沙群体沉速整体随含沙量增加而减小并逐渐趋于定值，而当平均含沙量为0.34~0.45 kg/m 3时，泥沙群体沉速出现随含沙量增加而增大的现象。

分析试验过程中泥沙颗粒垂向分选变化过程，当含沙量较低时，水流黏性仍可保持与清水相似特征，泥沙颗粒的分选作用使得沉速较大的泥沙群集中于含沙量较大的近底处，造成含沙量较大的区域出现泥沙颗粒群体沉速较大的反常现象。

这说明含沙量很低的条件下，泥沙颗粒非均匀特征伴随的各粒径沉速垂向分选未体现泥沙群体特征，仅表现出各粒径颗粒混合沉降的综合统计结果。

随着含沙量的增加，泥沙颗粒对水流黏性影响及颗粒间的相互作用成为制约泥沙颗粒群体沉速的主要因素。

第 1 期丰 青，等：黄河流域黏性泥沙群体沉速计算分析11
3 泥沙沉速计算式验证
利用沉降试验获得的泥沙群体颗粒沉速数据对上文建立的公式进行验证。

各组次泥沙颗粒群体沉速试验值与计算值对比见表2。

由表2可见，式（8）考虑了水流运动黏性切应力和紊动切应力影响，能够更完整地反映泥沙颗粒群体沉降过程中的能量损耗程度，可较好地呈现黏性泥沙群体沉降规律。

为区分泥沙颗粒沉降过程中非
均匀性引起的垂向分选对验证分析的影响，以每组
黏性泥沙沉降试验测点处浓度降至初始浓度50%的时间和含沙量变化计算泥沙群体沉速，并以此作为各组次不同含沙量对应的颗粒群体沉速数据。

由于泥沙粒径垂向分选的影响，验证计算过程中需进行不同组次对应稳态泥沙粒径试算。

先给定各泥沙粒径引起的水流复合黏滞系数初值，再利用垂向平均沉速确定各试验组次的水流复合黏滞系数。

试验水体垂向平均含沙量分别为2.0、6.0、12.0和16.0 kg/m 3时，相应的水流复合黏滞系数为7.15、6.10、4.57和4.35 mm 2/s 。

可见，随着含沙量的增加，水流复合黏滞系数呈减小趋势。

相同剪切速率条件下泥沙颗粒附近水沙界面处能量耗散强度逐渐减小，泥沙颗粒沉降过程中与向上运动的同体积水体之间相对位移有所减小。

表 1 垂向平均含沙量16.0 kg/m 3组次黏性泥沙群体沉速Tab. 1 Calculated settling velocity of cohesive sediment group for sediment concentration 16.0 kg/m 3
距底高度/m
积分计算时间段内固定位置处泥沙群体沉速/(cm/s)(0, 320] s (320, 640] s (640, 960] s 1.0
0.044 30.031 60.024 50.9
0.040 90.030 00.022 40.8
0.037 90.027 10.020 00.7
0.037 00.026 30.019 40.6
0.036 50.024 80.018 80.5
0.032 20.024 90.019 10.4
0.030 00.025 60.018 70.30.028 50.025 50.018 9表 2 泥沙颗粒群体沉速试验值与计算值对比Tab. 2 Comparison between the experimental and calculated settling velocity of sediment particles
各组次垂向平均含沙量/(kg/m 3）测点含沙量/(kg/m 3)沉速试验值/(cm/s)沉速计算值/(cm/s)
相对误差/
%各组次垂向平均含沙量/(kg/m 3）测点含沙量/(kg/m 3)沉速试验值/(cm/s)沉速计算值/(cm/s)相对误差/%2.0
1.5150.010 70.010 9 1.78
12.010.2690.017 00.016 0−5.722.0400.009 50.009 0 −5.5011.0230.013 20.015 215.192.3350.007 60.007 9 4.35
12.1500.013 10.014 0 7.106.0
4.9500.016 40.016 5 0.53
16.09.9160.038 70.037 9−2.145.6470.015 30.015 2 −0.82
12.6690.035 10.035 7 1.756.0560.014 40.014 4 0.2016.6990.031 10.032 5 4.476.5990.013 60.013 5 −1.00
21.5350.030 10.028 6−4.996.8400.012 90.013 0 1.13
22.2430.029 10.028 0−3.6712.0 4.2900.022 40.023 7 5.82
23.7170.026 10.026 9 2.916.8290.020 40.020 2 −1.07
24.6460.025 50.026 1 2.45
9.5230.019 90.016 9−15.22 1.00.5 1.5
2.00.4
0.200.6
1.2
1.0
0.8含沙量/(kg/m 3)
(a)
平均含沙量
2.0 kg/m
3距
底高度/m t = 320 s t = 640 s t = 960 s t = 1 280 s t = 1 600 s 10515200.40.200.61.21.00.8含沙量/(kg/m 3
)(b)
平均含沙量
12.0 kg/m
3距底高度/m t = 200 s t = 400 s t = 600 s t = 800 s
图 3 不同含沙量条件下不同时刻沉降桶含沙量垂向分布
Fig. 3 Variations in sediment concentration profiles at different time intervals with varying average sediment concentrations 12水 利 水 运 工 程 学 报
2024 年 2 月
4 结 语
针对黄河流域含有黏性矿物成分的黏性泥沙沉降特征进行分析，以泥沙颗粒为参照物进行了受力分析和能量损耗分析，建立了考虑挟沙水流运动黏性系数和紊动黏性系数的泥沙群体沉速显式表达式，并采用黄河流域典型黏性泥沙红色砒砂岩颗粒开展了泥沙群体沉速试验研究，得出主要结论如下：
（1）黏性泥沙进入水体不仅改变了水流的黏滞性，对泥沙颗粒水沙界面处的紊动剪切强度也将产生显著影响，黏性泥沙群体沉速降低是水流黏性、密度和能量耗散强度变化的综合表现。

（2）天然条件下黄河流域典型黏性矿物组成的黏性泥沙未出现显著的絮凝团结构，红色砒砂岩组成的含沙水体中泥沙群体沉速整体随含沙量增加而减小并逐渐趋于定值；当平均含沙量为0.34~0.45 kg/m 3时，出现泥沙群体沉速随含沙量增加而增大的现象。

采用试验数据对建立的公式进行了验证，公式能较合理地表达水流黏性和尾流分离对泥沙颗粒沉降阻力的影响。

（3）黄河流域黏性泥沙颗粒级配跨度较大，受试验条件限制泥沙沉降过程中可能出现垂向泥沙颗粒分选，对试验结果产生一定影响。

天然河道中，经过一定距离的水沙界面交换，河道较强的沿程分选作用下悬沙级配在一定范围内较为稳定，此时泥沙群体沉降速度计算可采用不考虑悬沙颗粒垂向分选影响的结果。

参　考　文　献：
钱宁, 万兆惠. 泥沙运动力学[M ]. 北京: 科学出版社, 1983: 60-72. （QIAN Ning, WAN Zhaohui. Mechanics of sediment
transport [M ]. Beijing: Science Press, 1983: 60-72. （in Chinese ））[ 1 ]路新飞. 含沙量和泥沙粒径对水流特性的影响[D ]. 郑州: 华北水利水电大学, 2018: 53-58. （LU Xinfei. Influence of sediment concentration and diameter on flow characteristics [D ]. Zhengzhou: North China University of Water Resources and Electric Power, 2018: 53-58. （in Chinese ））[ 2 ]
张晓华. 黄河宁蒙河段“新悬河”治理对黑山峡水库再造床能力的需求研究[R ]. 郑州: 黄河水利科学研究院, 2022: 35-48. （ZHANG Xiaohua. Study on the capacity of Heishanxia Reservoir to rebuild the bed by harnessing the “new suspended river ” in the Ningmeng reach of the Yellow River [R ]. Zhengzhou: Yellow River Institute of Hydraulic Research, 2022: 35-48. （in Chinses ））[ 3 ]
王兆印, 王文龙, 田世民. 黄河流域泥沙矿物成分与分布规律[J ]. 泥沙研究，2007(5)：1-8. （WANG Zhaoyin, WANG Wenlong, TIAN Shimin. Mineral composition and distribution of the sediment in the Yellow River Basin [J ]. Journal of Sediment Research, 2007(5): 1-8. （in Chinese ））[ 4 ]
田世民, 王兆印, 李志威, 等. 黄土高原土壤特性及对河道泥沙特性的影响[J ]. 泥沙研究，2016(5)：74-80. （TIAN Shimin,WANG Zhaoyin, LI Zhiwei, et al. Effects of soil property of Loess Plateau on sediment characteristics in the Yellow River [J ].Journal of Sediment Research, 2016(5): 74-80. （in Chinese ））[ 5 ]
杨铁笙, 熊祥忠, 詹秀玲, 等. 粘性细颗粒泥沙絮凝研究概述[J ]. 水利水运工程学报，2003(2)：65-77. （YANG Tiesheng,XIONG Xiangzhong, ZHAN Xiuling, et al. On flocculaton of cohesive fine sediment [J ]. Hydro-Science and Engineering,2003(2): 65-77. （in Chinese ））[ 6 ]
金德泽, 吴慧文, 师福东. 泥沙絮凝问题分析[J ]. 东北水利水电，1998，16(11)：25-27, 37. （JIN Deze, WU Huiwen, SHI Fudong. Analysis of sediment flocculation [J ]. Water Resourse & Hydropower of Northeast China, 1998, 16(11): 25-27, 37. （in Chinese ））[ 7 ]
VAN RIJN L C. Principles of sediment transport in rivers, estuaries and coastal seas [M ]. Den Haag: Aqua Publications, 1993.[ 8 ]
SOULSBY R. Dynamics of marine sands: a manual for practical applications [M ]. Oxford: Thomas Telford Press, 1997: 23-41.[ 9 ]
沙玉清. 泥沙运动学引论[M ]. 西安: 陕西科学技术出版社, 1996: 78-92. （SHA Yuqing. Introduction to the dynamics of sediment transport [M ]. Xi’an: Shaanxi Science and Technology Press, 1996: 78-92. （in Chinese ））[10]
窦国仁. 窦国仁论文集[M ]. 北京: 中国水利水电科学出版社, 2003: 243-259. （DOU Guoren. Dou guoren’s collected papers [M ]. Beijing: China Water & Power Press, 2003: 243-259. （in Chinese ））[11]
郜国明, 郭选英, 曾芹, 等. 非均匀颗粒群体沉速试验方法的分析[J ]. 泥沙研究，2000(3)：68-71. （GAO Guoming, GUO
Xuanying, ZENG Qin, et al. Analysis on experimental method of group settling velocity of non-uniform grains [J ]. Journal of Sediment Research, 2000(3): 68-71. （in Chinese ））
[12]第 1 期丰 青，等：黄河流域黏性泥沙群体沉速计算分析13
14水　利　水　运　工　程　学　报2024 年 2 月郜国明. 非均匀沙在紊流中沉降的试验研究[D]. 北京: 清华大学, 1997: 21-30. （GAO Guoming. Experimental study on [13]
settling velocity of non-uniform sediment in turbulent flow[D]. Beijing: Tsinghua University, 1997: 21-30. （in Chinese））张红武, 张清. 黄河水流挟沙力的计算公式[J]. 人民黄河，1992(11)：7-9, 61. （ZHANG Hongwu, ZHANG Qing. Formula of [14]
sediment carrying capacity of the Yellow River[J]. Yellow River, 1992(11): 7-9, 61. （in Chinese））
窦国仁. 河口海岸全沙模型相似理论[J]. 水利水运工程学报，2001(1)：1-12. （DOU Guoren. Similarity theory of total [15]
sediment transport modeling for estuarine and coastal regions[J]. Hydro-Science and Engineering, 2001(1): 1-12. （in Chinese））
费祥俊. 黄河中下游含沙水流粘度的计算模型[J]. 泥沙研究，1991(2)：1-13. （FEI Xiangjun. A model for calculating [16]
viscosity of sediment carrying flow in the middle and lower Yellow River[J]. Journal of Sediment Research, 1991(2): 1-13. （in Chinese））
韩其为. 水库淤积[M]. 北京: 科学出版社, 2003: 121-124. （HAN Qiwei. Reservoir sedimentation[M]. Beijing: Science Press, [17]
2003: 121-124. （in Chinese））
Examination of cohesive sediment settling velocity calculations in
the Yellow River Basin
FENG Qing1, 2, XIAO Qianlu2, ZHENG Yanshuang2, ZHANG Xiaohua2
(1. College of Environmental Science and Engineering,Peking University,Beijing 100084,China; 2. Yellow River Institute of Hydraulic Research, Yellow River Conservancy Commission, Zhengzhou 450003, China)
Abstract: Examining the settling velocity characteristics of cohesive sediment, comprising cohesive minerals within the Yellow River Basin, involves constructing a physical graph depicting the energy dissipation at the water-sediment interface of sediment particles. This construction incorporates force analysis and the energy conservation of sediment particles. The explicit expression for sediment group settling velocity is derived, utilizing the kinematic viscosity coefficient and turbulent viscosity coefficient of sediment-laden flow. This formulation effectively captures the influence of flow viscosity and wake separation on sediment particle settling resistance. Validation of the theoretical framework is achieved through settling velocity tests conducted on representative cohesive sediment, specifically red arsenic sandstone particles, within the Yellow River. The results indicate a monotonic decrease in the group settling velocity of red arsenic sandstone in sediment-laden flow with increasing sediment concentration, ultimately reaching a constant value. Conversely, for vertically averaged sediment concentrations ranging from 0.34 to 0.45 kg/m3, the group settling velocity consistently increases with sediment concentration. The theoretical analysis aligns with experimental observations, and the formulated expression is substantiated by collected data. This comprehensive study serves as a valuable reference for understanding the dynamics of deposition and transport of cohesive suspended sediments in the Yellow River Basin.
Key words: cohesive sediments; settling velocity; energy dissipation; viscosity coefficient; the Yellow River Basin。