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

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

垂线平均含沙量两种表述方法之比较水利论文+_3sN+Nri#S

黄才安 梅小文

HLQ T!yc f0

(扬州大学水利学院)水利论文;~ | zV A L#\2u'M[e+~

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

"g^H,EG `%E%i.}O0t0

关键词悬移质 悬浮指标 流量平均含沙量 面积平均含沙量

Y0@2u#nR0

1 问题的提出

w\J'y9BO hP2m0

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

!t+zC+f_ v,}0

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

S1=990113e1.gif (989 bytes)/990113e2.gif (972 bytes)水利论文7W.q0k2S$d

(1)水利论文llMt/~#e

S1=1/h·990113e3.gif (970 bytes)

5TI{1},O_|0

(2)水利论文X)j'gu/r5^*duL\

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

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

I$dt OI\T zV7_{0

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

+\B R#x|0

dS/dx=-[αωs/q(S-S*)]

/uh.BD(T@,A+n0

(3)水利论文tc L{ p

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

yi;y ^ d|0

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

e,f0\^3T E0

2 流量平均与面积平均含沙量的差异水利论文S7Y$|8DZY\ e0s|

  垂线流速分布采用著名的普朗特对数流速分布公式水利论文'e{i%[*j2X |\pD

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

U XVD'e"I us0

(4)水利论文Mw0N|!I(gkG9@

或写成水利论文J6wr&|g ^FW

u/U=[1+1/kC0(1+lnζ)]

l'B,mj;s6] X V6nt O0

(5)

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

  垂线含沙量分布采用著名的罗斯公式水利论文!W^w3X:a)gB

S/sa={ζa(1-ζa)·[(1-ζ)/ζ]}Z水利论文8_iE"]:xKk(Q

(6)

*VH'La,~N_0

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

H z/r`0z6F,h,h)K0

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

*s@A3?g+hoz0

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

q]4Y_#aJ0

S1=∫10usdy/∫10udy=sa[ζa/(1-ζa)]Z10{(1+1/kC0)·[(1-ζ)/ζ]Z+1/kC0·[(1-ζ)/ζ]Zlnζ}dζ水利论文Z4c?v&W

(7)

@q IJ M+GNQ0

若令

;f]-Re M z)F Q NT0

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

da(`9z Fil:w0

(8)

Bq;]v6g;z#\0

10[(1-ζ)/ζ]Zlnζdζ=J2水利论文&U-yBj|3Qr#D

(9)水利论文_/e7wk&N+e$R

则式(7)变为

G,k/U1Up D0

S1=saa/(1-ζa)Z]{(1+1/kC0)J1+1/kC0J2}水利论文YpN!Y"Vy:}$h

(10)

4W0G+W:nYI|4\ v0
  又利用式(6)可得面积平均含沙量

S2=∫10sdζ=[ζa/(1-ζa)]Z·sa10[(1-ζ)/ζ]Z

_n TW'R%J&Xq8]y(G0

(11)

)qp(Wm#q5dXtn0
同样将式(8)代入上式可得

S2=saa(1-ζa)]ZJ1

7~ KA0kLxQ0H }@5K`0

(12)

+u?!l%{\7I l7}%m O N0

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

{7v&oY1jo0

β=S1/S2=1+1/kC0(1+J2/J1)

:iG,{-Pm[1s`0

(13)

:H'x*RBp*dKZ(P!E#Y E0

由于J1、J2仅与Z有关,因而β也仅为Z与kC0的函数,即水利论文A;Ut$}dc%nj@

β=f(kC0,Z)水利论文 YN9aH,a2j3l

(14)水利论文)h-{ NqLNM

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

dPKv3hloF)L0

β=S1/S2<1水利论文:pPAyGThG

(15)水利论文XI Eg8Xg'z~z

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

!p.{6Q&VL&`0

  根据文献[3]中提供的J1、J2计算曲线,可计算β与Z、kC0之间的关系,如图1所示。水利论文M6@|5h;H

  从图1中可以发现β具有如下变化规律:水利论文 F.d$]0CcAM

  (1)随着Z值的增大,S1与S2相差愈大。悬浮指标Z是反映垂线含沙量分布均匀程度的物理指标,Z愈大,含沙量垂线分布愈不均匀;Z愈小,含沙量垂线分布愈均匀。所以说含沙量分布愈均匀,β愈小,反之,β愈大。水利论文o`;@5`#A'j i^ ]

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

990113t1.gif (2485 bytes)水利论文V)a OGX1Y0[

图1 β与Z、kC0的关系图

TA$t QAk L0

Relationship among β and Z、kC0

J.j%\p%]6a5W~'P0

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

W.@.P#k r1D0

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

3 S1与S2的应用探讨

g8n/a0N/n-B8G:J0

3.1 在悬移质输沙率计算中的应用水利论文6p9_s'u!V \YF&i

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

,w?'M'MU0

gs=∫h0usdy

(qIx6e+u.|$J'A:v-E0

(16)

)M!rkBF3TN1}0

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

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

v zc}&t%W'Z3v6i0

若以垂线平均含沙量表示,则为水利论文v0m7z"RE&xQ

990113t2.gif (2135 bytes)

7fSpG~/g(T3]*xw0

图2 式(17)验证结果

3W~9N,n/F o0

Verification of equation(17)水利论文,b-Z#oKe

gs=qS1=h990113e5.gif (869 bytes)S1990113e5.gif (869 bytes)hS2水利论文f-? f/@7f Dq+y

(17)

6b)w(H"nFw'B0

3.2 在悬移质扩散方程中的应用

7j7My7S-~%m#Kv0

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

990113e6.gif (855 bytes)(qS)/990113e6.gif (855 bytes)x+γ′990113e6.gif (855 bytes)Z0/990113e6.gif (855 bytes)t=0水利论文Xkv0G @e8h J

(18)水利论文|JD8ZD|g4v(~k

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

4 结语水利论文#T'DtM|3Jvr^

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

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

参考文献水利论文2}GY M_

1 张瑞瑾。河流泥沙动力学。北京:水利电力出版社,1989:185~203.

CL3l6\,| w g0

2 窦国仁。潮汐水流中的悬沙运动及冲淤计算。水利学报,1963(4):13~23.

/u[ g7S![[0

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

}Nf{}H{g I0

4 华东水利学院等四校合编。河流动力学。北京:人民交通出版社,1980:211~228.水利论文o6n0J.} A S*L-I0s

 水利论文M'ZcE7TyH/B

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

广告投放

广告投放