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

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

垂线平均含沙量两种表述方法之比较水利论文9J T{W{i

黄才安 梅小文

b g7Y [9s.bP5s6G0

(扬州大学水利学院)水利论文,N_(b r3I0DCt"RS

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

关键词悬移质 悬浮指标 流量平均含沙量 面积平均含沙量水利论文Thr ].icJ;^#W

1 问题的提出水利论文:X,^5~QCBK"U(G

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

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

zu1X*I {mV0

S1=990113e1.gif (989 bytes)/990113e2.gif (972 bytes)水利论文FW ulZ+d-? g&Dc

(1)水利论文 ?#xo%Pe E.]m

S1=1/h·990113e3.gif (970 bytes)水利论文v6Dzc D'o:wh4q

(2)水利论文!ly.]Ir d ~D!f

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

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

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

O.k`b;hc@2h;N0

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

o3A!q Wb?0

(3)

(U?:c M)i,MD0

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

4S rv@ t h*d _b$t0

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

2 流量平均与面积平均含沙量的差异水利论文^C FxVp0l"qP,y

  垂线流速分布采用著名的普朗特对数流速分布公式

z Y|n'o+Q9ut0

u/U*=(1/k)·lnζ+C水利论文S0G*T6t3u-Ov fiR#D

(4)

El { BLW R/AM0

或写成水利论文'u_Nrh Q4V j-NJ-u

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

x.dbN8|(`{;C:m(vM0

(5)

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

  垂线含沙量分布采用著名的罗斯公式水利论文 pWEL ^,a

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

Y\.c&R"L[2@3\0

(6)水利论文(OtSL;^o

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

0uKn s.EG1A0

式(5)及(6)所表示的流速及含沙量垂线分布规律,尽管还存在一些不足之处,但由于两式的计算精度较高,且有一定的理论基础,因而得到了广泛的应用。水利论文"vZD|u!o#@H2[(\$W

利用式(5)及式(6),可得流量平均含沙量水利论文m!S2O AIcQZ:C%\

S1=∫10usdy/∫10udy=sa[ζa/(1-ζa)]Z10{(1+1/kC0)·[(1-ζ)/ζ]Z+1/kC0·[(1-ζ)/ζ]Zlnζ}dζ水利论文IvB*AR ^

(7)水利论文 n])C8h8Ze\1m7Q;mg?

若令水利论文(\ d rz o7B

10[(1-ζ)/ζ]Zdζ=J1水利论文7Th"k(rB2[E[1k

(8)

\'By"@ TH0

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

'Gq7v5P[$y0

(9)

yMS(Y*\&[O&O0

则式(7)变为水利论文UZ5R&@[?

S1=saa/(1-ζa)Z]{(1+1/kC0)J1+1/kC0J2}水利论文'J*{A$W j?7vM

(10)水利论文)\7h1U.U @*~*v,I

  又利用式(6)可得面积平均含沙量

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

kP~'U X%])k5ew4c0

(11)

x)f_rR/Z(C0
同样将式(8)代入上式可得

S2=saa(1-ζa)]ZJ1

U;ge Zam W"Uzx0

(12)水利论文5zi~%AQ0rW-|

  流量平均含沙量与面积平均含沙量之比为水利论文k5z2? L"w0r&p

β=S1/S2=1+1/kC0(1+J2/J1)水利论文 puAg QJ t

(13)

)A@Z}2C~0

由于J1、J2仅与Z有关,因而β也仅为Z与kC0的函数,即

#{J.g6R(N p$m0

β=f(kC0,Z)水利论文/mU;s9U7Izn

(14)水利论文E3k*z(J$E!S^

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

!dlqWI8^;Ytbi L0

β=S1/S2<1水利论文w1@ TAP2x8|2O|

(15)

X Ci-NL!QOn0

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

  根据文献[3]中提供的J1、J2计算曲线,可计算β与Z、kC0之间的关系,如图1所示。

V k_|q.N:gdpb0

  从图1中可以发现β具有如下变化规律:

{y7\R6s0

  (1)随着Z值的增大,S1与S2相差愈大。悬浮指标Z是反映垂线含沙量分布均匀程度的物理指标,Z愈大,含沙量垂线分布愈不均匀;Z愈小,含沙量垂线分布愈均匀。所以说含沙量分布愈均匀,β愈小,反之,β愈大。水利论文@5l1\+^ ?4T \

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

990113t1.gif (2485 bytes)水利论文6Xo&UAZY!O

图1 β与Z、kC0的关系图

C(V+[HQmh E"zH0

Relationship among β and Z、kC0水利论文g qWj8{(`dF

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

*i k~ kh |}O~0

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

3 S1与S2的应用探讨

N%{u ~qX+S0

3.1 在悬移质输沙率计算中的应用水利论文5zp8| Y4I9G4S

  悬移质输沙率的计算,常使用下式水利论文|#AR ? ~&[.q'IlG }8q

gs=∫h0usdy水利论文0[0s"m1f$D UZ

(16)

P8x"d*e6]0n;X1[%\0

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

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

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

X_8\ w3[VP LvG0

990113t2.gif (2135 bytes)

5bIHz.@P"@0

图2 式(17)验证结果水利论文FvMl/n \)g8A

Verification of equation(17)水利论文 bv4b!U_$x%~ Qn

gs=qS1=h990113e5.gif (869 bytes)S1990113e5.gif (869 bytes)hS2水利论文j$?w$pIZ5pbuy*U

(17)水利论文5c'UJ8`)\hJR

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

*u1q/H-P4]s,TQ0

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

t5scyc&ll2H0

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

Q1H4L svD0

(18)

Z$n4U`M]0k%gH0

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

4 结语

|{ ~0f1Ur.B M2_2p0

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

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

参考文献水利论文 W5ZH:By%`'r v |c

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

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

3 H.A.爱因斯坦(钱宁译).明渠水流的挟沙能力。北京:水利出版社,1956:20~22.水利论文~;Hz?2W f*u:c

4 华东水利学院等四校合编。河流动力学。北京:人民交通出版社,1980:211~228.

-~JW7ov:Z5I^0

 水利论文-IO H N @*{'{iT-g P%i

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

广告投放

广告投放