河流动力学——第五章

sv
L dSv dy
,v
du L
dy
vsv
L2
dU dy
dSv dy
usv
x
dSv dx
vsv
y
dSv dy
y
L2
dU dy
最后可得二维水流中悬移质泥沙的扩散方程：
对流项
Sv t
U
Sv x
(usv x
)
(vsv y
)
ω
Sv y
vsv
y
dSv dy
usv
x
dSv dx
水流紊动引 起的扩散项
y]xt
g
S vt x
yt
U t S vt yt
[ g
Svt x
x
(
g
Svt x
)x]yt
[U t
S vt
(U t S vt x
)
x]yt
g
S vt y
xt
Vt Svt xt
Svtxt
x方向：
由于水体流动而进入和流出该隔离体的染色物 质的差值为 ：
Ut Svtyt
[Ut
S e v
( ya)/ y
S va
Rouse试验
为了检验悬移质含沙量的垂线分布公式是否合理， Rouse (1938)在圆筒中加一组等间距的格栅，使格栅在 园筒中作上下简谐振动，从而在较长一段垂向距离内得 到均匀的紊动流场。通过测量四种不同粒径的泥沙在不 同浓度的垂线含沙量分布。试验结果表明含沙量分布基 本符合式( 5-15 ) 的理论曲线。
Sv t
U
Sv x
(usv x
)
(vsv y
)
ω
Sv y
g
2 Sv x 2
2 Sv y 2
对流项
水流紊动引 起的扩散项
重力 沉降项
分子运动引 起的扩散项
由于紊流中流体微团随机运动的规模远大于分子 热运动的规模，即第四项与第二项相比要小得多，故一 般可以忽略分子扩散项的影响，则上式可化简为：
线为清水中的动量交换系数理论分布。
可见，在这一试验所研究的含沙水流中，浑水中的动
多沙河流中的泥沙输运大部分是以悬移运动的形式 进行的。
例如，长江泥沙以悬移质输沙为主，约占总输沙量的90% 。
在三峡水库蓄水卵前石，推移长质江宜昌站沙多质推年移平质 均的卵石悬推移移质 质
(沙D量>1约0m宜寸为昌滩m8站站)6年2万输t沙，量而约悬(2D87>为.6移9170万7m万质6tmt万)(Dt5，0=沙0.0质(D351推860=60m020移.2万万m1m质tt)m年)(D输50沙=0量(.D25则1054=.m.206达6.0亿m亿31到tm)t 年m5).输26
y
dSv dy
Sv
0
y
m
ku*(1
y)y h
h y κU* y h
dSV dy
SV
0
dSV SV
κU*
1 y
h
1
y
dy
悬浮指标
Z kU*
ln SV
ln
h
y
y Z
lnC
令 y＝a 为参考点，该点
泥沙浓度记为Sva ,则可得：
Z
Sv Sva
h y
y
h
a
a
悬浮指标
Z kU*
悬浮指标Z反映了重力作用与紊动扩散作用的比值，其 大小决定了泥沙在垂线上分布的均匀程度。Z值越小，悬移 质分布越均匀。
Einstein认为，Ｚ=5 可作为给定的 水流和泥沙条件下，泥沙是否进入悬 浮状态的临界判别值。
Einstein还提出，在有床面形态时， 上面推导过程中的剪切流速U* 都应代
之以沙粒阻力对应的剪切流速 U*。
Rouse公式的缺陷
Sv Sva
h y
y
h
a
a
Z
1、得出的只是相对含沙浓度，没有参考值的话，
紊流中沿水深不同高度处各层水体之间存在水团的紊动 交换，其结果是形成一个向上运动的泥沙通量qs1。
另一方面，由于泥沙比水重，往下沉降形成一个向下运 动的泥沙净通量qs2。
如果悬移质含沙量沿垂线出现稳定的时均泥沙浓度分布， 说明qs1与qs2达到了动平衡状态，即悬移质含沙量沿垂线分布 达到平衡状态。
研究悬移质含沙量沿垂线分布的有关理论：扩散理论和 重力理论。
Rouse试验表明，对某些特定的条件，可以取y为 常数。而对于天然河道，可以近似取为：
y 0.067u*h
进一步的研究表明，泥沙的扩散系数y 不是常数而 是空间位置的函数，但现有的理论还不能给出y 沿垂线
的分布规律。
最常用的方法是假定泥沙扩散系数y 与动量交换系 数 m 相等。剪切紊流中相邻流层所存在的剪切应力，
=0.4为Kárman常数
对上式微分可得流速梯度：
du dy
u* ky
对二维水流来说，剪切应力在垂线上呈线性分布，
0(1
y) h
床面处最大，水面处为0。
则：
y
m
du
0 (1
y) h
du
u*2 (1 u*
y) h
ku* (1
y )y
h
dy
dy
ky
Rouse公式（悬移质含沙量沿垂线分布公式）
式中负号表示溶解物质总是从浓度高的地方向低 的地方扩散。
在二维水流中，设Ut、Vt分别为纵垂向的瞬时流速, Svt为瞬时含沙浓度。如果在水中注入一团染色物质，由 于水流的扩散作用，这团染色物质随着水流向前运动的 同时，将不断向上下左右扩散。由于垂向上的时均流速 为零，所以垂向上的染色水体扩大与时均运动无关，而 是纯粹的扩散作用引起；纵向上的染色水体扩大与扩散 作用和时均剪切离散作用都有关。
亿t（，涪占陵长到江全河段部）输沙量(D的50=5918m.m2)%。 (D50=0.14mm)
(D50=0.034mm)
寸滩水文站资料统计表明，涪陵长江河段多年平均卵石推移
质(D50研=5究1m认m为)输：沙推量移为质2输8.沙97量万与t，悬沙移质质推输移沙质量(D相50比=0，.14在m平m原) 为 60地0万区t河，流而中悬仅移占质到(D总50=的0.输03沙4m率m的)输1%沙~量5为%4；.6丘亿陵t，地占区到河全流部为输沙 量5的%9~81.57%%；。山区河流为15%~30%。
dSv dy
Graf和Cellino利用声学颗粒流速仪技术进行了系列
水槽试验，得到了m和y值沿垂线分布的实测结果。试
验表明，在平衡输沙条件下(含沙量等于挟沙能力)，
存在床面形态时，沿垂线平均的>1，
❖没有床面形态时，沿垂线平均的<1。
右图为根据试验资料、按定义式
计算得到的泥沙扩散系数y和动量交 换系数m沿垂线分布(<1的情况)。图 中实点据为值y ，虚点据为m 值，实
数y与动量交换系数m 相等的缘故，不少研究者进行了
试验研究，结果表明两者并不相等，应作出如下修正：
y m
Z
kU*
计算比例系数的大小时，需要通过试验得出m和y值，
一般是通过测量流速和含沙浓度的脉动值、以及含沙浓度 的时均值Sv，再根据定义式进行计算如下：
m uv
dU dy
y vsv
§5.1 泥沙扩散方程
以悬移形式运动的泥沙颗粒具有较细的粒径，可 以随着水流的紊动在水体中随机运动。
在垂向上，泥沙颗粒的运动可以看成是两种运动的 叠加，即重力驱动下的沉降运动和水流紊动驱动下的 随机运动。
当颗粒的数量很大时，宏观上将形成泥沙垂向运 动的动态平衡，此时泥沙浓度在垂向上有一个稳定的 分布。
研究泥沙的扩散时，由于泥沙比水重，沙粒还将以 速度ω下沉。
在扩散物质到达的断面取一个长度为Δx，高度 为Δy，厚度为1的水体，分析在时间Δt内进入或者离 开这个水体上下左右四面染色物质的体积。
[εg
Svt y
y (εg
Svt y
)y]xt
[Vt Svt
(Vt Svt ) y
y]xt
[Sv
Svt y
河 流 动 力 学 概 论
黄河壶口瀑布
第五章 悬移质运动和水流挟沙力
5.1 泥沙扩散方程 5.2 悬移质含沙量的垂线分布
扩散理论、泥沙扩散系数的试验研究、重力理论
5.3 悬移质输沙率
Einstein公式、维利卡诺夫和张瑞谨公式、Bagnold 公式
5.4 水流挟沙力
理论公式、经验或半经验公式、天然河流挟沙力公式
沉降项
Sv t
U
Sv x
x
( x
Sv ) x
y
(
y
Sv ) y
Sv y
同样可推导出三维水流中悬移质泥沙的扩散方程：
Sv t
U
Sv x
x
( x
Sv ) x
y
(
y
Sv ) y
z
( z
Sv ) z
Sv y
§5.2 悬移质含沙量的垂线分布
泥沙的重力沉降使得含沙水流中沿垂线形成上清下浑的 浓度分布。
Sv t
U Sv x
(usv x
)
(vsv y
)
ω
Sv y
确定因水流紊动而引起的泥沙扩散输移率 usv 和 vsv ，一般有两种方法。
一是与 Fick 定律直接类比，即假定泥沙的扩散输移 率与泥沙的浓度梯度成正比：
usv
x
dSv dx
vsv
y
dSv dy
二是借用紊流模型中的混掺长度理论，设浓度和垂 向流速的脉动可写为：
在主流区可以看作主要是由于流体脉动导致的相邻流层
间动量的交换所引起的，可以仿照分子粘性应力的表达