垂线平均含沙量两种表述方法之比较(黄才安 梅小文)

热度144票 浏览85次 【共0条评论】【我要评论 时间:1999年1月01日 14:31

垂线平均含沙量两种表述方法之比较水利论文N u9D _6D%t._1I

黄才安 梅小文水利论文EQ7JR6?&j"[

(扬州大学水利学院)

_ s lMaBi2A%]&Q0

摘 要表达沿垂线平均含沙量的方法通常有流量法与面积法两种,但在实际运用中常发生混淆。本文讨论了两种垂线平均含沙量之间的关系,指出面积平均含沙量要大于流量平均含沙量。含沙量及流速分布越不均匀,两者差别就越大,当悬浮指标Z≤0.3时两者相差不到10%。最后还对两种垂线平均含沙量的应用作了初步探讨。水利论文&Z L F&V/i2\

关键词悬移质 悬浮指标 流量平均含沙量 面积平均含沙量水利论文4~~8u pe s ^I

1 问题的提出水利论文g:pvzW

  悬移质垂线平均含沙量是泥沙运动力学中的一个重要内容,是研究悬移质输沙率与水流挟沙力的基础,同时也是研究河床变形问题的依据,在生产实践中具有重要意义。水利论文mnFX`7C2Y M

  为了研究、计算上的方便,常需考虑悬移质含沙量的垂线平均值,由于含沙量及流速分布在垂线上是不均匀的,因而存在多种垂线平均含沙量的表达方式。垂线平均含沙量的表达通常有所谓的流量法与面积法两种,它们的定义分别为[1]水利论文{"J| X*W]

S1=990113e1.gif (989 bytes)/990113e2.gif (972 bytes)

Z,lrM Sl._0t-oa0

(1)

;h?O5VW0

S1=1/h·990113e3.gif (970 bytes)水利论文x6JpS'D/yNS

(2)水利论文~gz,o R{G

式中 S1、S2分别为按流量法与面积法计算的垂线平均含沙量(以下分别简称为流量平均含沙量与面积平均含沙量),u、s为距河底y处的当地流量与含沙量,h为垂线水深。

^'t"zGji;J F5y0

  两种垂线平均含沙量均具有一定的物理意义。前者S1表示的是流量加权平均值,其物理意义是单位时间内通过单宽断面的悬移质输沙率与单宽流量之比,即分部输沙率平均;而后者S2表示的是面积加权平均值,其物理意义是单位时间内通过单宽断面的沙量与单宽面积之比,即分部沙量平均。在实际应用中,两种垂线平均含沙量也各具特色,流量法由于同时考虑了流速和含沙量这两个因素,且流速为向量因素,这样求得的垂线平均含沙量S1,常用于与悬移质输沙率有关的计算中;而面积法仅考虑含沙量这一非向量因素,这样的垂线平均含沙量S2一般多用于理论研究中及一维或二维垂线平均悬移质扩散计算中。

q4i0E8m.E^0

  但在一些关于悬移质泥沙运动的研究文献中,经常会出现将S1与S2混淆,甚至是相互代用的现象。例如在计算悬移质输沙率gs时,常用单宽流量乘以面积平均含沙量S2来推求,实际上此处应以流量平均含沙量S1代替S2。又如在一维悬移质扩散方程中[2]水利论文(ojP)__O f.Ej&b.fe

dS/dx=-[αωs/q(S-S*)]水利论文6y MO@u@9idC-L

(3)水利论文aS/z f-NRWv3R"OIm

式中α为泥沙恢复饱和系数,ωs为泥沙沉速,q为单宽流量。式中的垂线平均含沙量S及挟沙能力S*应取S1还是S2,常出现混乱。同样在垂线平均的平面二维悬移质扩散方程中,也常有类似的混乱之处。

Mr JniXY b6\ gm;]0

  无论是在数学形式上,还是在物理意义上,流量平均含沙量S1与面积平均含沙量S2是两个截然不同的概念,一般不可相互代用。本文试图对这两种垂线平均含沙量之间的差异作一分析,以求能澄清某些概念上的问题。水利论文6It0__z8o(Y

2 流量平均与面积平均含沙量的差异

gl zBV7]H0

  垂线流速分布采用著名的普朗特对数流速分布公式水利论文?(v.K%h#`

u/U*=(1/k)·lnζ+C

Fq gBH+p0

(4)水利论文(ZHz$u!O

或写成

9yu`NVHW0

u/U=[1+1/kC0(1+lnζ)]水利论文g}k_hW4D1G

(5)水利论文 p{jpT

式中 k为卡门常数,u是相对水深ζ(y/h)处的点流速,是垂线平均流速,U*为摩阻流速,C0为无因次谢才系数(C0=C/990113e4.gif (901 bytes),C=/U*)。

  垂线含沙量分布采用著名的罗斯公式

7Wt"X$Iy/M yY0

S/sa={ζa(1-ζa)·[(1-ζ)/ζ]}Z

t\y'O3|!xwI?,O0

(6)水利论文5t(r6m2Rp F7I

式中 s为相对水深ζ处的点含沙量,sa为ζ=ζa处参考含沙量,Z=ωskU*为悬浮指标。

_9j&N0H3R}+N0

式(5)及(6)所表示的流速及含沙量垂线分布规律,尽管还存在一些不足之处,但由于两式的计算精度较高,且有一定的理论基础,因而得到了广泛的应用。

XE v D { l9J*~5Z^0

利用式(5)及式(6),可得流量平均含沙量

(c5\6~-b%Zb:Y2z0

S1=∫10usdy/∫10udy=sa[ζa/(1-ζa)]Z10{(1+1/kC0)·[(1-ζ)/ζ]Z+1/kC0·[(1-ζ)/ζ]Zlnζ}dζ水利论文x;MzG V{3f A

(7)

{MX lo XH0

若令

aur8wR2Y F[.I0

10[(1-ζ)/ζ]Zdζ=J1

H4s lg dSH5s{[0

(8)

(w1tQ"Y^,e0

10[(1-ζ)/ζ]Zlnζdζ=J2

2T3x3TnJ*|q{0

(9)

byz `%O Sx0

则式(7)变为水利论文Va"u-^!^xt

S1=saa/(1-ζa)Z]{(1+1/kC0)J1+1/kC0J2}

O.u6R.~3i8D3N0

(10)

YzWo&}T"I DP0
  又利用式(6)可得面积平均含沙量

S2=∫10sdζ=[ζa/(1-ζa)]Z·sa10[(1-ζ)/ζ]Z水利论文.N7U@9B _,FP

(11)

`-X2c2U N{:{0
同样将式(8)代入上式可得

S2=saa(1-ζa)]ZJ1

;SE\%Bf-T0

(12)

N!lMf P&m?0

  流量平均含沙量与面积平均含沙量之比为

b&J _G5y(F,^7hO^0

β=S1/S2=1+1/kC0(1+J2/J1)水利论文B+f}R x"k#gB})y)@y

(13)

7s#b_f!m(l0

由于J1、J2仅与Z有关,因而β也仅为Z与kC0的函数,即水利论文9g#d#L;t-eo*v0V

β=f(kC0,Z)水利论文$l%Jb;Z*f E}d:E

(14)水利论文#h Kj ~Zd:p5Q`.e

  从数学上可以证明,函数J1与J2具有如下性质:①J1>0,J2<0;②|J1|<|J2|。所以有

z `;^"{w"P6}s%WE~0

β=S1/S2<1水利论文D3b7ci#C7KU OO

(15)

"|+I9I-t+~b{"eF0

式(5)说明面积平均含沙量永远大于流量平均含沙量。此差异完全是由流速分布的垂线不均匀性引起的。

8V&?%uj"N4u;T+~Bx7T9x0

  根据文献[3]中提供的J1、J2计算曲线,可计算β与Z、kC0之间的关系,如图1所示。水利论文!U N(`/u@d ]fY

  从图1中可以发现β具有如下变化规律:水利论文3~N`L wX#oe#X

  (1)随着Z值的增大,S1与S2相差愈大。悬浮指标Z是反映垂线含沙量分布均匀程度的物理指标,Z愈大,含沙量垂线分布愈不均匀;Z愈小,含沙量垂线分布愈均匀。所以说含沙量分布愈均匀,β愈小,反之,β愈大。水利论文:y9H&Z0~9KW

  (2)随着kC0值的增大,S1与S2相差愈小。kC0实际上是反应垂线流速分布均匀程度的一个物理指标,这里称之为流速分布指标,kC0愈大,流速分布愈均匀;kC0愈小,流速分布愈不均匀。因此,流速分布愈均匀,β愈小,反之,β愈大。

zM4ez+} d;s0

990113t1.gif (2485 bytes)水利论文WW(m'P9sIP)q

图1 β与Z、kC0的关系图水利论文UN"Fm ZH {#[

Relationship among β and Z、kC0

OcEe` L(I:W-CF0

  (3)一般情况下流速分布指标kC0=5~10,而悬浮指标Z=0.01~5,因而从图中可知,悬浮指标Z对β的影响要远大于流速分布指标kC0对β的影响,即含沙量的垂线分布均匀程度对β的影响要大于流速分布均匀程度对β的影响。水利论文eN0_4hs2|VQ'TZ ^ gY

  (4)通常情况下,当Z≤0.3时,S1与S2两者相差不到10%,此时可近似地以S1代替S2或以S2代替S1,对计算结果影响不大;当Z>0.3时,S1与S2两者相差较大,不可相互代用。

H Qy'q&j9cU|C0

3 S1与S2的应用探讨水利论文)Pq#s*|i,];N,i

3.1 在悬移质输沙率计算中的应用水利论文Y_b-m"L1b

  悬移质输沙率的计算,常使用下式

IA E+]"l K8Ac9sd#a V.D0

gs=∫h0usdy水利论文!d;r9J:X0VmIZ9X x9i

(16)水利论文!fH-v[%Yg&o

  利用式(17),对长江下游计134条垂线进行了验证计算,结果如图2所示,其中gs计算值系用β990113e5.gif (869 bytes)hS2求得,gs实测值系水文测验值(根据式(16)推求)。

%s,Vx U'In{0

  由图2可见,用式(17)计算悬移质输沙率有一定的精度。也就是说利用S1与S2均可计算悬移质输沙率,两者计算结果并无不同,仅是计算方法有所区别。同时利用S2计算悬移质输沙率使计算来得更为简单,从而为悬移质输沙率的计算提供了新的途径。水利论文.qyH_ up ^

若以垂线平均含沙量表示,则为

E%`oTE!RJ0

990113t2.gif (2135 bytes)

%RY'lME#iec3^'^0

图2 式(17)验证结果

!M(f;q|,{(a0

Verification of equation(17)水利论文|O-I%UK Ky

gs=qS1=h990113e5.gif (869 bytes)S1990113e5.gif (869 bytes)hS2

G T;Rr?#V0

(17)

R4sQ u?.Z,RU?B0

3.2 在悬移质扩散方程中的应用水利论文`8a3W"wT(J:j5?

  在一维及平面二维悬移质不平衡输沙计算中,都要涉及到垂线平均含沙量及挟沙力的概念。这里仅以一维悬移质扩散方程式(3)为例说明其中的垂线平均含沙量应取S1还是S2。式(3)中的垂线平均含沙量及挟沙力,使用面积平均含沙量S2或流量平均含沙量S1似无本质区别,仅要求等式两边均采用S1或均采用S2,因为两者之间仅相差一个β值,等式两边可约去。但在使用中应注意两点,其一是挟沙力S*一般均是利用实测水文资料加以率定的,而水文资料刊布的垂线平均含沙量多为流量平均含沙量S1;其二是由悬移质扩散方程求出的垂线平均含沙量沿程分布规律,大多应用于诸如下式[4]的悬移质河床变形方程中水利论文c} X5t;a

990113e6.gif (855 bytes)(qS)/990113e6.gif (855 bytes)x+γ′990113e6.gif (855 bytes)Z0/990113e6.gif (855 bytes)t=0

6n*z~/I&H.Y0

(18)

-E7\/v p-K%w0H7ex0

式中 γ′为泥沙干容重,Z0为河底高程。鉴于这两点考虑,作者推荐在式(3)中均使用流量平均含沙量S1,若要使用面积平均含沙量S2,则需作相应换算。水利论文 ^&m@0JZr

4 结语

RnH*Z/B A.Uq0

  目前泥沙研究中,通常采用的流量平均含沙量S1与面积平均含沙量S2,无论是数学表达形式,物理意义,还是具体计算结果都不相同,切勿误用。在对两者差异分析的基础上,本文认为随着悬浮指标Z的增大,两者的差别越来越大,且S2永远大于S1,当Z≤0.3时,两者相差不足10%,此时两者可近似替代,而当Z>0.3时,两者不可相互替代。

QJj"r(o"HC&Y0

  在实际应用中,特别是一些与悬移质输沙率有关的问题及与水流挟沙力有关的不平衡输沙问题中,要特别注意垂线平均含沙量两种表述方式的选用。一般来说,选用任何一种表述方式均可,当选用面积平均含沙量时,要进行适当转换,而当选用流量平均含沙量时,只会使问题得到简化,并使计算更为方便。因此,作者推荐采用流量平均含沙量。

-]&N)I ]5B\0

参考文献水利论文KNx;Tq:XYv5G)av

1 张瑞瑾。河流泥沙动力学。北京:水利电力出版社,1989:185~203.水利论文)iOZ+X!t'v F

2 窦国仁。潮汐水流中的悬沙运动及冲淤计算。水利学报,1963(4):13~23.水利论文Zzc`2MA\6[#E

3 H.A.爱因斯坦(钱宁译).明渠水流的挟沙能力。北京:水利出版社,1956:20~22.水利论文 lviSg0Y'I8v-Y

4 华东水利学院等四校合编。河流动力学。北京:人民交通出版社,1980:211~228.水利论文F;`%l@4|&B3xpzp

 

(D/E6b;p/a1ejw0
TAG: 含沙量 黄才安 梅小文
顶:21 踩:28
【已经有95人表态】
26票
极差
9票
很差
11票
较差
12票
稍差
9票
稍好
6票
较好
12票
很好
10票
极好
下一篇:流域年均含沙量的PP回归预测(李祚泳 邓新民 侯宇光)
上一篇:潮汐河口悬移质动床实物模型的理论与实践(熊绍隆 胡玉棠)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放