非均匀沙起动规律研究(李荣 李义天 王迎春)

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

非均匀沙起动规律研究

;D#DlKgFQ0

李 荣 李义天

P5W.yj3F6rT @3W0

王迎春

p&v"u(D5w {@M/k0

(武汉水利电力大学)水利论文1IX SC8\p

(吉林市水利水电勘测设计院)

!rG!i(j7q0

摘 要非均匀沙运动规律与均匀沙存有明显的不同。本文从非均匀沙运动特点出发,简要阐述了非均匀沙的运动机理,利用滚动平衡方程式建立了能反映非均匀沙粗化、细化过程和床沙受暴露、隐蔽作用的起动流速公式,利用实测资料率定计算公式中的系数。与其它起动流速公式相比,本文推导的起动流速公式较好地反映了非均匀沙的起动特性:较细的非均匀床沙较粒径均匀的同粒径床沙难于起动,而较粗的非均匀床沙较粒径均匀的同粒径床沙易于起动。

1mT-ysM R0

关键词非均匀沙 起动流速 相对暴露度

&Mr%@(xu-H"~T'A.s0

1 引言水利论文uB+[.XFs0D

  一般情况下,天然河流中床沙均属非均匀沙,对于山区宽级配卵石夹沙河床,非均匀性更强,拣选系数990105e01.gif (998 bytes)高达10左右。由于非均匀沙运动机理的复杂性,以往众多学者把主要精力集中在均匀沙运动机理和输沙特性的研究上,而暂时回避了泥沙的非均匀性,并取得了一些成果。但泥沙研究的最终归宿应是研究能反映天然泥沙特性的非均匀沙运动机理和输沙特性。在非均匀沙研究方面,国内外学者均进行了一些研究。不同学者从各自的途径,考虑不同影响因素,得到了不同的起动流速公式。主要途径有两类:一类是基于力学角度推导出来的起动流速公式,另一类是基于随机理论推导出来的公式。

1Z'zrlH-x-d!Eygl0

1.1 基于力学角度推导出来的起动流速公式

K6W d&QB9N0

 秦荣昱从力学角度来建立和分析非均匀沙的起动流速公式。秦在分析非均匀沙起动力学特性时引进了附加阻力R,并近似假定R与混合平均抗剪力τc成正比,即R=φτcαd2利用滚动平衡推导得起动流速公式[1]。秦荣昱非均匀沙起动流速公式用于均匀沙时与沙莫夫均匀沙起动流速公式一致,公式引进了附加力在一定程度上反映了床沙粗化、细化对非均匀沙起动的影响。但对宽级配床沙而言,秦荣昱公式在计算粗颗粒泥沙时起动流速较实测值大,而计算细颗粒时略偏小。分析其原因,主要是公式中对粗颗粒泥沙受暴露,细颗粒泥沙受隐蔽的特性未能体现完全。而颗粒的隐暴特性对非均匀沙的起动影响较大。水利论文5O$x:i!\A g

类似于秦荣昱的思想,张启卫引进了一个附加力R(d<dm时,表现为附加阻力,d>dm时,R表现为附加动力),R表示为R=K(γs-γ)d3dm/d·ln(dm/d),利用力的平衡关系推导得非均匀沙的起动流速公式[2]。张启卫公式在一定程度上考虑了粗细沙之间的隐暴作用,但未能考虑床沙粗化细化过程对起动的影响,因而所得公式不能全面反映床沙调整过程中非均匀沙的起动规律。陈媛儿、谢鉴衡在研究非均匀沙起动规律时限于讨论粒配范围较宽的非均匀沙中不同粒径级的临界起动条件,不考虑床沙级配的调整,得到非均匀沙的起动流速公式[3]。但陈、谢公式中,当d=dm或d<dm时,出现间断或负值,这不符合泥沙起动的规律,而且公式也只限于不考虑床沙组成变化,河床遭受到一定程度冲刷情况下的粗颗粒泥沙的起动。水利论文:t!rZkmDM(vC

  另外,耶格阿扎罗夫[4]和林太造[5]从起动拖曳力出发研究非均匀沙的起动。但由公式可知,耶、林在推导非均匀沙起动时仍沿用了均匀沙的起动公式,只是略加修正,这也不可能全面反映非均匀沙的起动规律。

vL}{ Fk0

1.2 基于随机理论推导的起动公式水利论文%zVl|})R

  泥沙起动既应满足一定的力学规律,同时也是一种随机过程。在分析泥沙运动规律时,由于泥沙问题的复杂性,欲对每个随机变量都进行分析目前尚有一定困难。孙志林主要从对非均匀沙起动起决定作用的水流脉动和颗粒间的相互位置这两个随机变量出发,运用统计理论建立了反映非均匀沙起动特点的起动流速公式[6]。然而,我们不难发现,对水流脉动和颗粒间位置的研究其本质仍是对颗粒受力的分析研究,另外由于近底水流结构的复杂性及起动标准的不统一,想用随机理论来较好地解决非均匀沙起动问题目前仍存在相当大的难度。因此目前对非均匀沙的研究基于力学理论的居多。水利论文mY LR%h_

 由上面分析可知,秦荣昱公式的特点是引进了附加阻力项,该项表征床沙组成变化对起动粒径d0的影响,这就概括了克拉米提出的床沙颗粒从停止不动image4.gif (852 bytes)开始起动image4.gif (852 bytes)全部起动的整个过程的起动机理和物理意义。起动公式中通过dm,m,d90等的变化基本反映了床沙组成粗化、细化作用对起动粒径的影响。与其它公式相比,秦荣昱公式较好地反映了非均匀沙起动的阶段性和不恒定的特性。但秦荣昱公式在与实测资料的吻合程度上表现为粗颗粒泥沙起动流速较实测值偏大,细粒泥沙偏小的问题,反映了秦公式对泥沙受隐暴作用缺乏全面的考虑。本文在秦荣昱研究基础上,引进相对暴露度,试图对非均匀沙起动规律作进一步的探讨和完善。水利论文0[XgK0`

2 非均匀沙运动机理

Y*w-C W#m^/q0

  根据已有非均匀沙的研究成果[7],非均匀沙与均匀沙在起动机理,运动特性等方面存在显著差异,这也是进一步研究非均匀沙的根本原因。非均匀沙与均匀沙相比有如下几个显著特点:

}2R4q$oUK0

2.1 起动的阶段性

)rQ7M+w ]Bke[0

  根据Kramer的观点[8]和各种试验资料,非均匀沙起动存在四个不同阶段,即床沙颗粒不动状态;床沙中细小颗粒开始运动而粗颗粒泥沙静止不发生输移的状态;床沙中中等颗粒泥沙分选起动以及床沙全部冲刷输移阶段。水利论文r,MdS g U"Y N+r6A'i

2.2 非均匀沙起动的不恒定性水利论文a6o~+C9o w,Xhj T^

  在Kramer定义的泥沙起动第三阶段中,随着细颗粒泥沙的冲刷外移,床沙中泥沙组成变粗,床沙粗化,可动泥沙颗粒逐渐减少,这样势必会降低河道的输沙能力(即发生粗化过程);从另一方面讲,如果河段上游来沙充分且组成较细,河床发生淤积,床沙组成逐渐变细,有利于提高河道的挟沙能力(即发生细化过程)。故粗化、细化过程是天然河道冲淤交替的反映,同时反过来也影响河道的冲淤特性。由此可见:非均匀沙的起动和输沙过程是床沙粗化、细化的过程,具有不恒定的特性。

@:bhv [ Nu0

2.3 非均匀沙的隐暴特性

e#Xef6b }dO0

  随着水流强度的增加,床沙在水流作用下发生分选。细颗粒泥沙被冲刷外移,填充在粗颗粒泥沙的副流区,受到粗颗粒泥沙的阻挡和隐蔽,除受阵发性向上漩涡的作用而可能逸出隐蔽区外,细颗粒泥沙处于被保护,静止不动的状态;与细颗粒相反,粗颗粒泥沙随着床沙分选而逐渐从床沙中“突现”出来,水流对粗颗粒泥沙的作用增强。水利论文9@Mr*e s8?$zG

3 起动公式的建立水利论文#p&TIq0E$Zq.y

  依据非均匀沙起动机理,同均匀沙相比,非均匀沙除受水下有效重力(ω),水流对泥沙的上举力(FL)和正面推力(FD)外,还应受到因床沙粗化、细化作用及隐暴作用所附加的力R(对粗颗粒来讲,R是促进粗沙运动的力,对细沙则是阻碍运动的力)。R不仅与床沙组成粗细,床沙的粗化、细化过程有关,更重要的应反映床沙受隐暴的特性。

w%s:Rs}tGLY0

  从刘兴年[9]定义的床沙暴露高度(如图1)出发,本文引入相对暴露度的概念,其定义为

gu\l5@LP+i:K0

η=(dm-d0)/dm水利论文x-K{ f9s,e&t

式中η为相对暴露度;dm为床沙平均粒径;d0为其中任一组泥沙粒径。由定义可知:均匀沙时η=0。粒径大于床面平均粒径的泥沙,其暴露度为正,受到暴露作用。相反,粒径小于平均粒径,则受到隐蔽作用。这种相对暴露度的定义,不仅反映了某一颗泥沙在床面上的位置,也反映了与该颗粒泥沙相邻颗粒在床面上的位置。相邻颗粒间的位置关系正好反映了非均匀沙大小颗粒间的相互作用关系,因而反映了非均匀沙颗粒在床面上的位置特性。

?]"M,Fb'WOM0

990105t1.gif (1683 bytes)水利论文 Jh'@,~G+RI(A

图1 床沙暴露高度示意图水利论文(EEn;F4j

Exposed degree of bed sediment

_*{+e~.P V0

  由于现有对作用力R的研究很少,尚无精确理论公式,通过分析,附加力R应是相对暴露度η,床沙组成P(用d90,d60,d10反映),床沙粗化细化因子(用紧密系数m表示)等的函数。写成函数关系式为

xWc"\d`,eN0

R=f(η,P,m,……)

h7l?Q \4H;M0

近似假定附加力R与混合沙的平均抗剪力τc成正比。按照[1]推导得τc=Kmm(γs-γ)dm,式中Km为无因子系数;τc为混合沙的平均抗剪力;m为紧密系数,与不均匀度有关,可根据水工手册确定。于是附加力R写为

@[b4C*VE4x4s0

R=φτcαd02η=φKms-γ)mdmαd02η

,V/q/D|s nch2G S0

(1)水利论文c6{G#R7t

式中φ为系数;α为面积系数。

.K/o[;a/axz;W0

  本文定义的附加力R与[1]定义的附加阻力有所不同。本文的附加力R可以是附加动力也可以为附加阻力。当d0>dm时表现为动力,反之表现为附加阻力。而[1]在任何情况下,R只可能是阻碍泥沙运动的力。文献[2]中附加力R在d>dm时表现为附加动力,d<dm时表现为附加阻力,在这个概念上与本文相同,张启卫在定义R时也着重考虑了粗细沙的隐暴特性,但与[2]相比,本文定义的R中多了一项反映床沙级配和床沙变化的项,即不均匀系数m。故本文定义的附加力R从泥沙受力机理上来讲更合理些。另外由于本文定义的附加力是在秦荣昱研究基础上引进相对暴露度的概念,因此本文的附加力R是反映床沙起动阶段性、不恒定性和隐暴特性的一个综合作用力。水利论文%RR,\G sl

  据[10],泥沙在水下所受有效重力、正面推力及上举力分别写为水利论文&w+\)tw3Cf![A\'iKq

ω=(γsαd03水利论文 zH#P$JO2S!l%O

(2)

3{ }${2?M9e0

FD=CDαd02γ(ub2/2g)  水利论文5wb'L%?A"iB

(3)水利论文g*[^8^;H4Q"O1_duB

FL=CLαd02γ(ub2/2g)

WDt/P-n0

(4)

gD8t!I2km0

式中α0、α1、α2分别为与力相应的面积系数,利用指数流速公式u=(n+1)/n·U(y/h)1/6取y=D90处的流速u=ub,将公式(3)、(4)中的ub换成断面平均流速U。

Q EaS7n*B[X0

FD=CDα1d02γ(D90/h)1/3(ub2/2g)

7P1iI-@p1Q b#q!UJ0

(5)

I}T"pi9N4L0

FL=CLα2d02γ(D90/h)1/3(ub2/2g)水利论文[(l'_3N"]

(6)水利论文z8r\'wX]h,?$l

据此建立滚动平衡式

u4q[xD1t R0

FDLD+FLL2=ωL3+RLL4

~;f'z h4~m0

(7)水利论文'N[qd m!H*@

式中 L1、L2、L3、L4分别为FD、FL、ω、R的力臂,并在此设L11d0、L22d0、L33d0、L44d0,将L1、L2、L3、L4和(1)、(2)、(5)、(6)式代入(7)式,整理得水利论文Te s)Yz

Image55.gif (1980 bytes)

pp8XqpDK0

(8)水利论文1g x&Q%W2Z(j

Ψ1=CDα1β1+CLα2β2

wu4m Y?0

Ψ2=2gα0β3·s-γ)/γ水利论文&k5JHQBlb

Ψ3=2gαβ4φKm·s-γ)/γ水利论文S8R+X-IL;[W*PqN

(8)式转化为水利论文%WErco-I`

Image59.gif (1537 bytes)水利论文yk sB&ie&j

(9)水利论文%[eeJ0u'q

令 B131,B221则上式写为

Image60.gif (1438 bytes)水利论文 \n*wN'j!GW)W

(10)

aEU9oZ0

式中 B1、B2为待定综合系数,须经实测资料确定。整理Lane和Carlson,冷魁的粗化试验资料共47组次经曲线拟合如图2,取B1=20.85,B2=15.26。水利论文*U @o4c0x

  式(10)为不均匀沙各个粒径组的起动流速成公式。式中m,dm,D90等反映了床沙组成,粗化细化及暴露度对起动流速的影响。当均匀沙时式(10)可写成均匀沙起动公式水利论文S%r;r'[b

U=3.91d1/3h1/6水利论文H4P1n W&I'D"G}

(11)

RI$Z W$b w5q0

990105t2.gif (2476 bytes)

{*o)a1G;zo9Q S0

图2 公式系数率定图

c$x CVK`$|R\G sP0

Determination of coefficient of formula

+Q'l`T(K2x+@0

  上式与沙莫夫均匀沙起动公式一致,系数略编小。当为粗颗粒泥沙时,非均匀沙起动条件小于均匀沙同粒径泥沙的起动条件,细颗粒泥沙起动条件大于均匀沙同粒径泥沙的起动条件。比较式(10)与式(11)可知:公式(10)多了B1mdmη项,该项正好反映了床沙组成,粗化细化作用及隐暴作用等因素,由此可知公式(10)较好地反映了非均匀沙起动是受床沙级配、粒化程度及隐暴度等综合作用的规律。

l2\s&G4CO)c&i4Q0

4 公式的检验与比较

}:C-~tf7w)b0

  为了验证公式(10)的合理性,本文利用孙志林在武汉水利电力大学做的非均匀沙起动资料,Gessler水槽试验资料,Little和Mayer资料以及若干组野外实测资料在粒配及粒化程度相同或基本相近的情况下对公式(10)进行了检验,并在相同或相近的情况下同时对秦荣昱非均匀沙起动流速公式,张启卫起动流速公式以及谢鉴衡公式进行计算、比较,如表1所示。水利论文Q7?Zp jlY2a j,m"m

表1 起动流速(单位cm/s)

:_%jg/av.X0

Threshol d velocity and threshold grain size水利论文S.N;s*cK[


来源

D0(mm)水利论文 eXa-fnj/m

0.4

6w:hdD0M7X*|0

1.0水利论文 nRE K(N"U3q G)o!|

2.0水利论文c6EX8b"IX

3.2水利论文9x~5XaY%j9i

4.0水利论文{9JA&w3V6q1k^c

4.8水利论文5h ]"JY$b `

6.0水利论文9Tx [aPS$^e+k


实测值水利论文'?4v l WM [

18.5-20.3

I'gY/ScyiM0

24-28.5

MuO$B A5Pd0

35.0水利论文w/t'pXB&UEK

39.2水利论文,]+y ol"c1f"XD

40.1-42.0

F }3P@#k;\ l0

46.6

M:D5r/i\ gOR0

50.8-53.0水利论文 R Vz X~v{C

秦公式水利论文']v+t:K Q6?

15.40

L"Q{"d7{0

21.5

6_#o&c {HvDH'i0

40.0

8s\$~b-S0

47.4水利论文HJa[O"Ak/e7A

48.5水利论文 xu:aN+`

52.2水利论文*xbJ&_-U2v%w2X

53.5水利论文'FJ%{1g,II+? ^'[Mw9U#S

张公式

rFjiru)^0

20.3水利论文 lexE:w

30.0水利论文-[a*fZ#g!l-C(b8B

38.2

1U5es*k7B\0~5Tjx%Z0

47.6

EA.J:PI?R0

50.0水利论文{ p q a+hI

54.1水利论文"h {Gb K2v*SA?

58.5水利论文,Xg1dM;}0Kyj7^4b

谢公式

\HLa!a:MIXH.T0

18.0水利论文6g_,WM a2G wQ

23.4水利论文l)f:W;V.x

29.0水利论文a1^&z7M&P] ?

32.0

8c-Fo1i T0O:{&d0

34.0

3G,B k6A7P(O3YqO/`0

35.6

{.Yz-O'Mq?0

38.0水利论文,q;S,l$NpY

本文公式水利论文5Ls g c+UL

19.8

mc\(z)wh^3t I0

26.3

S#K)oN8@/C a'I0

34.8

Em)Q7u_7~RS0

40.1

1{(I6MP\'QY? E!DH;Y0

40.5水利论文Y4x*A(w-N S

45.8

jk R.C$jR&\0

51.2水利论文.izch-h3h


  由表1不难看出:(1)秦荣昱公式粗粒径泥沙部分计算结果较实测值偏大,细沙部分又比实测值略偏小,表明秦荣昱公式对粗颗粒泥沙受暴露,细颗粒泥沙受隐蔽的特性反映的不够完全。(2)张启卫公式在细沙部分与实测值吻合较好,但当用于粗沙计算,或用于河床遭受一定冲刷,床沙级配调整后的起动流速计算时,计算值比实测值大,由前面的分析可知,张启卫附加力公式中主要考虑粗细沙的隐暴性,而对床沙级配及其变化对泥沙起动的影响缺乏具体体现,这正是导致偏差的原因。(3)陈媛儿、谢鉴衡起动流速公式与实测资料偏差较大,计算值比实测值小,这与模型的简化假设有关。(4)从表可以知道,本文公式(10)计算值与实测资料吻合较好,基本反映了非均匀沙起动的基本规律:较细的非均匀床沙较粒径均匀的床沙难于起动,而较粗的非均匀床沙较粒径均匀的床沙易于起动。因此本文公式可用于非均匀沙不同粒径起动流速的计算以及一定水力条件下最大起动粒径的确定。

2T9b+i3@J"K^0

5 结论

6S9m&PO2A0

  通过以上分析研究,本文获得如下几点认识:水利论文,s/})Rcb6W

  1.非均匀沙起动具有阶段性,非恒定性和粗细颗粒的隐暴等特性。

3Sh U-N5b.\#sm XU0

  2.非均匀沙中粗颗粒泥沙同均匀沙中同粒径沙相比较易起动,而细颗粒泥沙恰好相反。

N-ZYc7L`"b l~{0

  3.通过引入相对暴露度的概念,推导并建立了能反映非均匀沙起动的阶段性,非恒定性和粗细颗粒的隐暴等特点的起动流速公式。公式(10)与实测资料吻合较好。水利论文9[,{i5M5xm4t

  4.由于泥沙起动规律和近底水流结构的复杂性,加上现有实测资料较少,对宽级配(m>0.70)的非均匀沙,公式(10)的实用性有待进一步的检验。

U an P,GO;Y5q5{0

参考文献水利论文5o.? j$w\

1 秦荣昱。不均匀沙的起动规律。泥沙研究,1980(复刊号).水利论文+R(zSa8_(N"@

2 张启卫。非均匀沙起动流速。全国泥沙基本理论研究学术讨论会论文集,1992,11.水利论文7k4L._{4Xg,?.v

3 陈媛儿,谢鉴衡。非均匀沙起动规律初探。武汉水利电力学院院报,1988(3).水利论文 fXD]-s

4 Egiazaroff, i.v., Calculation of Nonuniform Sediment Concentration, J.H.D.ASCE,V.91,No. Hy4,July,1965.

a%t oU9v H)_V!S"z0

5 Hagashi, T.S., Qzaki and T.Ichibashi, Study on Bed Load Transport of Sediment Mixture, Proc. 24th. Jap.Conference on Hydraulics, 1980.水利论文J,VG3l b6Aq*o2AE/eW `

6 孙志林。非均匀沙输移的随机理论。武汉水利电力大学博士论文,1996年。

U p S"LF]e0

7 秦荣昱,王崇浩。河流推移质运动理论及应用。中国铁道出版社,1996.水利论文&N-p] Ry#^.}m Rj9W

8 Kramer, H., Sand Mixtures and Sand movement in Fluvial Models. Trans. ASCE, Vol.100,1935.水利论文 d;Ov^E*d(_,t

9 刘兴年等。卵石河道宽级配推移质输移特性研究。第二届全国泥沙基本理论研究学术讨论会论文集,中国建材工业出版社,1995年8月。水利论文!m$BW:G3KN&kkz'}#D O.hB

10 武汉水利电力学院。河流泥沙工程学(上).水利出版社,1982年。水利论文 `2`,VQ{3t f3q4M

 

$r6P*D/ruJh2I1w0
TAG: 李义天 王迎春 李荣 规律 起动
顶:28 踩:24
【已经有111人表态】
17票
极差
11票
很差
14票
较差
14票
稍差
17票
稍好
11票
较好
13票
很好
14票
极好
下一篇:低含沙水流运动特性综述(刘兆存 徐永年)
上一篇:水利立交地下涵洞输水输沙特性研究(徐金环)
查看全部回复【已有0位网友发表了看法】