水库三角洲河道输沙之研究(俞维升 李鸿源)

热度118票 浏览38次 【共0条评论】【我要评论 时间:1999年7月01日 13:57

水库三角洲河道输沙之研究水利论文{M)S? wX.g&SR#}

俞维升水利论文J?T3}2zv*s'X

李鸿源

G;Rc Px9[,u0

(农业工程研究中心 台湾)

%s)c*q(B4z/d:`Q1C2~a0

(台湾省水利处)

gxH+J5x0

摘 要水库进口附近三角洲河道之水流系近似均匀流,并非真正之均匀流。三角洲河道动床之自我调整作用,因水库水位之回水影响,其输沙特性有异于明渠均匀流。传统上,将明渠均匀流平衡条件下之输沙公式应用于估算三角洲河道之输沙能力,似有其不确定性。本研究进行一系列低水流区之输沙试验,并纳入文献[10]Hotchkiss 之试验数据一并分析。结果显示,在水库之回水影响下,河床自我调整作用使其输沙能力大于假设明渠均匀流之情况,其中Engelund Hansen公式低估低水流区之输沙能力达10倍;Parker公式亦低估,尤其是无因次颗粒剪应力小于0.08时差误更大。依据试验成果,本研究并建议一经验公式,以供估算水库三角洲河道之输沙能力。水利论文+f1F2ai2k(^:TPj

关键词三角洲 泥沙运动 边壁校正 水库淤积水利论文$\+LU^5{iC-J

1 前言水利论文O6J j3^3y.MKtw ?,U

    泥沙随着水流进入水库后,由于坝体的壅水,向下游水深渐增、流速渐减,水流挟沙能力因而降低,造成粗粒泥沙呈现超载状态并发生淤积,逐渐形成三角洲。依据三角洲床面的坡度变化,自上游往下游可将三角洲区分为尾部段、顶坡段及前坡段,顶坡段的渠道称为三角洲河道,能够到达三角洲前坡段下游的泥沙一般皆为细颗粒,水库三角洲纵向剖面示意如图1。

Z5\E5@6N0

990302t1.gif (2188 bytes)

'Ju+\u ~&Z!zPv0

图1 水库淤积示意图[11]水利论文{7a_Xu^ D9D'FeD

Schematic diagram of reservoir sedimentation

o'Lt6_)m\4X \$ZK1W0

990302t2.gif (1718 bytes)

MyD7|1{@b*I0

图2 河川来沙超载造成底床升高之示意图[6]水利论文+w8?w1\ h

Schematic diagram of raising river bed due to overloading sediment水利论文/wFP1eg

    归纳以往模拟水库顶坡段泥沙冲淤的数值模式,计算输沙量的方法有二种,一种是针对河床推移质采用平衡条件下的输沙公式,针对悬移质以泥沙扩散方程式求解[1~3],其中泥沙扩散方程式可反映悬移质的不平衡输沙行为。另一种是采用平衡条件下的总输沙量公式[4,5]水利论文/q%N'J2Zi gE7|

    明渠水流输沙公式的数据来源,不外乎现场实测及水槽试验量测,通常要求量测数据时的水流条件是均匀流,量测段必须达到水面平行底床,水深是正常水深,在上述均匀流条件下,运动泥沙与床面泥沙的交换达到平衡。严谨而言,现有输沙公式大多只适用于均匀流且泥沙冲淤平衡的情况。

3ua#s!vt2] X0

    Hotchkiss and Parker指出在明渠中淤积造成底床升高的情况有两种[6],一种是上游来沙量大于该断面输送能力所造成,此情况下,该断面自我调整增加河床坡度以增加输沙能力,底床升高过程示意如图2,另一种是下游回水效应使该断面水流输沙能力降低,此情况下,该断面为增加输沙能力的自我调整作用受回水影响,在三角洲河道即为此种情况,第二种情况之底床升高如图1所示,明显与第一种情况不同。第一种情况,在底床升高的过程水流仍为均匀流,仍可引用平衡输沙公式;第二种情况,水流非均匀流且受回水影响,其输沙特性为本文研究重点.Hotchkiss and Parker指出三角洲河道仅是近似均匀流[6](near uniform flow),并非真正的均匀流,本文认为以往将平衡条件的输沙公式应用于计算三角洲河道输沙量,存在不确定性。水利论文'E+@5G*z(] s

    Sugio and Okabe指出文献[7]的试验之三角洲是由河床推移质淤积所形成[8]。Graf指出文献[2]的数值成果显示三角洲主要是河床推移质淤积造成,悬移质淤积主要发生在前坡段下游[9]。Graf整理前人研究,认为三角洲主要是由河床推移质淤积所形成[9]。Hotchkiss以长12.2m,宽15cm,高38cm的水槽做顶坡段输沙试验[10],福氏数Fr=0.33~0.64,在这种福氏数情况下,通常不会发生高水流区(upper flow regime)的床形,劳斯数(Rouse number)R0=Ws/κU*,其中Ws是颗粒沉降速度,κ是冯卡门常数(von Karman constant),U*是水流剪力速度,当R0>5时可不考虑悬移质[11],Hotchkiss的试验R0=1.09~1.43存在悬移质,雷诺剪力参数(Reynolds shear number)Re*皆大于35,Re*=U*d50/ν,其中ν是水流的运动粘滞系数,通常Re*>11.6时,沙粒突出层流次层(laminar sublayer),不会发生沙纹(ripple),底床无因次剪应力θb=0.384~1.461,Celik and Rodi提出当Re*≥0.6时θb只要大于0.25即有泥沙再悬浮发生[12],所以Hotchkiss的试验应会有悬移质,此与前述R0<5存在悬移质相符。由上述试验条件可知Hotchkiss的试验之床形属低水流区(lower flow regime)中输沙量较大的沙垄(dune)情况。至于低水流区中悬移质很少以河床推移质为主的低输沙量情况,则是本文水槽试验重点,本文并将文献[10]Hotchkiss的沙垄数据一并纳入,以分析三角洲河道的输沙特性,试验成果并与平衡输沙公式比较,进而分析前人以平衡输沙公式计算三角洲河道输沙量的误差。

lpr*UvH0

2 试验设备与方法

'D!ptT[:I o0

    本研究在预备试验时,曾进行如图1般底床有一固定坡度的水槽试验,检验结果发现尾部段的溯源淤积将使三角洲河道的来沙量无法维持一固定值。所以本研究设计如图3的水槽,该水槽长3m,其中1m是明渠段、2m是水库段,宽15cm,底床坡度接近水平。试验水槽在明渠段与水库段交界处有一落差10cm的阶梯,试验时间内控制明渠段无淤沙,如此可确保试验时间内三角洲河道有固定的输沙量。为使明渠段水流稳定进入水库段,控制顶坡段淤沙高程接近明渠段底床高程。

?yLZzm h^~0
    本研究试验泥沙有三种,平均粒径(mean diameter)dm分别是1.341mm、0.759mm及0.398mm,各种泥沙之几标准偏差(geometric standard deviation)σg皆小于1.4,可视为均匀粒径,各种泥沙的粒径特性如表1所示。泥沙比重为2.6635。本文共进行13组试验,皆为量流且定量加沙的试验,试验条件如表2所示,流量有三种分7别是3142l/s、4169l/s、5217l/s,含沙体积浓度在0.00006~0.0007间,加沙量在0.83g/s~5.83g/s间。试验步骤依序为1.升高尾水板蓄水;2.水面线稳定后,自上游施放固定的流量;3.水面线稳定后,以固定的速率加沙;4.调整尾水板高度使明渠段无淤沙,并使顶坡段淤沙高程接近明渠段底床高程;5.继续加沙,待顶坡段长度大于60cm后,量水面高程;6.停止供水及加沙,徐徐将尾水板放倒,以淤沙完全不移动为原则;7.量底床高程,每个断面侧向量取8个位置的高程。

990302t3.gif (2493 bytes)

IO}:N(y U!P/V%CZ0

图3 试验水槽布置水利论文/i.NI S(JV%[

Layout of experiment flume水利论文VI'M`K\+n0y5C

  表1 试验泥沙的粒径特性

4S`)s9L3`c"\8ECU0

Particle properties of experimental material(单位mm)

+_)],dKb F ~+[0

项目dmd50d84.1d15.9σg

第一种泥沙1.33411.3301.6101.0501.238
第二种泥沙0.7590.7390.8780.6321.179
第三种泥沙0.3970.3770.5000.3061.278

表2 试验条件与成果

;B5FyFlA0

Experimental conditions and results

J"O+]cwh(j&{ W0

Run水利论文Q]_)~]!V;|R F

(1)水利论文.t[3m @O%J#Xh4_

Q
(Cm3/s)
(2)
Gs
(g/s)
(3)

99030232.gif (835 bytes)

j.G}?H~m0
(by volume)
(4)
d50
(cm)
(5)
H
(cm)
(6)

S0水利论文_.]%bXo;oV

(7)水利论文q e{u@N/g


C13141.91.670.0002000.03775.2840.0069
C23141.93.330.0003980.03774.7470.0062
C33141.90.830.0000990.03776.1690
C43141.91.110.0001330.03775.8470.0015
C55217.01.670.0001200.13307.7950.0043
C64168.53.330.0003000.13306.1980.0064
C74168.51.670.0001501.13305.9480.0045
C85217.00.830.0000600.13306.6790.0077
C94168.53.330.0003000.07395.9880.0080
C104168.51.670.0001500.07396.7480.0223
C115217.05.830.0004200.13305.8260.0154
C125217.05.830.0004200.07396.3710.0058
C133141.95.830.0006970.03774.4390.0079

3 试验结果

Z(o;HK,dH4T-G2IA0

3.1 量测数据

EMs%z"@x/[0

    典型的淤积剖面如图4所示,图4是分别在不同时刻停水量测的淤积剖面,可看出淤积剖面很稳定,在三角洲顶点及顶坡段与明渠交界处(即图4中A与B位置)底床因受局部流场影响而较不稳定,但在位置A与B之间,底床坡度与水深接近不变,可量得顶坡段代表性之水深及底床坡度。水利论文 gv3U c/k ]c

990302t4.gif (1817 bytes)水利论文Z r7mc,m)e @5eO T

图4 本研究典型的淤积剖面变化历程水利论文0i;c[C)I

Typical sedimentation profile varied in time

p*xq m4_0

    本文每个断面侧向量8个位置的底床高程,虽然宽深比约为3有边壁效应存在,但由所量测数据可知侧向上底床变化不大。另外,由柏努力方程式(Bernoulli equation)可知底床坡度与水深的斜率同级,所以量测值中水深明显较底床坡度有较高的可信度,亦即底床坡度很容易因量测技术及数据处理上的些微差异而产生较大的变化,本研究在初步分析中发现所量测的底床坡度不具规律性,所以本文量测之底床坡度仅作参考,未当作分析数据。实验结果水深H在4.747~6.748cm间,顶坡段底床坡度S0在0~1.54%间,示如表2。

*q&ts9R,d K/w%W0

3.2 计算影响泥沙运动的有效剪应力

QcUM&CpO0

    由于水槽边壁是丙烯酸树脂(acrylic resin)制,底床是泥沙,不同材质使得剪应力在湿周上并非均匀分布,且水槽宽15cm,宽深比约为3,所以边壁与底床不同材质对泥沙输移的影响不可忽略,应进行边壁校正,从总应力τ中分离出与输沙有关的底床剪应力τb。此外,底床剪应力中克服床形起伏的形状阻力所消耗能量与河床质输沙量无关。仅克服泥沙颗粒表面摩擦阻力所消耗能量与河床质输沙量有关,所以需将底床剪应力中与河床质输沙量有关的粒径剪应力τg分离出来。本文边壁校正采用1942年Johnson提出的方法[11],分离粒径剪应力采用1952年Einstein and Barbarossa提出的方法[11],本文并提出校正剪应力的步骤,校正剪应力的流程详述如下:

#H,|3sD$z"a+ct0

    水槽边壁是丙烯酸树脂制的光滑边壁,流过光滑面的水流摩擦因子,可由慕迪图(Moody diagram)中光滑紊流之摩擦因子与雷诺数关系曲线求得,Granger提出下式足以表现该关系曲线[13]

'cl:SOj T2C@n0

990302e1.gif (1337 bytes)

'p"q9?G\{H0

(1)

F8qs"a.Y P0
式中 w为与边壁有关,Rew为雷诺数,Cw=fw/8,fw为达西-韦斯巴摩擦因子,Cw为修正的摩擦因子,剪应力τw=ρCwU2,其中ρ为水的密度,U为水流平均速度。

0BCIRRow0    (1)式中Rew=4rwU/ν,rw是水力半径,ν是水的运动粘滞系数,因为rw为未知数所以无法直接求得Cw值。Johnson针对边壁校正作以下假设[11]

o2f$nu'g7uJGX:C0

A=Ab+Aw

$H3e wBj)J:fy0

(2)水利论文uj1N'~)a*i)O%K}$b

U=Ub=Uw

E9d/n#sawk0

(3)

*t"m2M*m*P0

Sf=Sfb=Sfw水利论文#gcV{] p-| D4~[T

(4)水利论文 L iks4y'r Ps

上面诸式中b表示与渠床有关,无下标的参数是综合渠床与边壁整体效应的参数,A为水下横断面积,Sf为能  量坡度。

ei|Zt r7J\qQ0

Re/Cf=4rU/vCf=4U2*U/gSf/vCf=4*U2Cf/gSf*U/vCf=4U3/vgSf

"R_2E+E9R0

(5)水利论文2e5}ax;p,K

同理

1Q7rz0Z9fS(eR;@s|0

Reb/Cb=4U3b/vgSfb水利论文k`W5A7U

(6)

$_3}h;T:J9Nwl,n}7^0

ReW/ CW=4U3W/ vgSfW水利论文K,Irs5Qan

(7)

/Y {*leP~s jF0
由(2)~(4)式及(5)~(7)式可知下式成立

Re/Cf=Reb/Cb=Rew/Cw水利论文X QIVsM

(8)

{J!?e _ c)me0
由(8)式可知,(1)式中Rew可由下式求得

Rew=ReCwCf

{.~k%o7B*D3We0

(9)

eK7Dh2p0sy0
(9)式中Cf至目前步骤仍为未知数。依据基本的定义

U2*=grSf

yM"s5@6M*j9F+h0

(10)

#J0sY0m-yy'Ke0

U2*=CfU2

(W0}v3s j:wGUc0

(11)

$E4ho4VL1YX0
由(10)式及(11)式得

grSf=CfU2水利论文_2I8kE#M+C

(12)

'}*U nflG0m0
又r=A/P,P是湿周,所以

A=CfU2P/gSf水利论文b2I \z y

(13)

(_P"C1r-N!B'L0@4A0
由(2)、(3)、(4)及(13)式可得

PCf=PbCb+PwCw

L1ve h/cS0

(14)水利论文1O7g S0ua5]HSc-u

(14)式中Cb至目前步骤仍为未知数。Engelund and Hansen提出渠床剪应力与颗粒剪应力关系图[14],该图中的关系曲线可以下式表示(15)

θg=θb水利论文w3W)F)V7bw@KG

θb≤0.0615 床形为平床与沙纹

v pY Fp:J0

θg=0.06+0.4θ2b

4l5\?N"Yq Bv0

0.0615≤θb≤0.541 床形为沙垄水利论文%waj(\Y HW dp

θg=θb水利论文/`Q:EJ\9xx

θb>0.541 床形为平床水利论文q8hc[hU/O5m4`jG

上式中,θ=U2*/(s-1)gd50,其中s是泥沙比重,d50是中值粒径,θ称为希尔无因次剪力参数,下标g表示与颗粒有关的参数。(15)式提出了Cg与Cb的关系,针对Cg本文采用Einstein Barbarossa方法计算[11]

990302e16.gif (1270 bytes)水利论文$V,P(NGT |

(16)

o}rX8T0
上式 ks是渠床糙度高,本文令ks=d50,X是水流未达完全粗糙紊流的修正因子,X与ks/δ有关,δ是层流次层厚度,δ=11.6ν/u*g,由(11)式及(10)式知990302e17.gif (1157 bytes)。White,Milli and Crabbe分析文献[15]建议的X与ks/δ关系[16],提出(17)

X=1.90+0.7383ln(ks/δ)水利论文[Nq6t&A8K+D

0.265≤ks/δ≤0.5

%|\bc!dvC%O2m0

X=1.615- 0.407{|ln(ks/δ)|}1.6水利论文7V#q|1pK_J[ Et @

0.5≤ks/δ≤2.35水利论文Z;y4u1Bfk

X=1+0.926{1-0.434ln(ks/δ)}2.43水利论文(zj5?E7l` t-u

2.35≤ks/δ≤10

W2j$p eCn g:^ a e0
X=1.0

10<ks水利论文-?:Va0{.J0P&E;E

当ks/δ<0.265为平滑水流,Einstein建议以下式求Cg[15]

990302e18.gif (1371 bytes)水利论文J?5d9Bl e:W,F'zOa

(18)

"oD `!nL H J0
(16)~(18)式中rg可由(12)式求得,即

rg=CgU2/gSf水利论文P[#S%d @4v9ad7l

(19)水利论文[1^#B"h4U vT'k1k

由本文的试验量测无法得知Sf值,但Sf可由(12)式及(14)式联立解,即

Sf=U2/gr*PbCb+PwCw/P水利论文.Ix&v!p0Y;I W&E

(20)

-QB @xPT0
水槽宽度为b,水深为h,所以(20)式可改写如下

Sf=U2/gr×(bCb+2hCw)/(b+2h)水利论文,`Y,CdNDq]

(21)

5rB8E K8@;nT+O0B0
前述(1)、(9)、(14)、(15)、(16)、(19)及(21)式中共有7个未知数:Cw、Cb、Cg、Cf、Sf、rg、Rew,以迭代法可求解,本文建议之求解步骤如下
1)假设Cg=0.0025~0.0065

iY*yE}Vnc02)计算δ值,δ=11.6v/u*g=11.6v/990302e23.gif (917 bytes)水利论文#V1m G-Y0Tj

;n @%\i.`5z8]03)计算ks/δ值,取ks=d50

0|4[/_IX%^T0

/v?f-QLo}04)由(17)式、(18)式计算X值水利论文1]dy9k2Ws.`4{$K

水利论文U7eX V#}g0i

5)计算Cf

5Z$hRRo5O GPv0
    由(19)式及(12)式得

rg/Cg=r/Cf水利论文-[5s d C~~O i

(22)水利论文ke4P C,A3@*ct:Sh

    将(22)式代入(16)式得

Cf=12.3rCgX/ks*e*-1/2.5990302e23.gif (917 bytes)水利论文:MKX Kc"p0qHL j

(23)水利论文|5H;k5LWK

   将假设之值与前述计算所得X值代入(23)式,可得Cf

6)计算Cw水利论文vfbP y\1ua

将(9)式代入(1)式得

Cf=ReCw/6.9*e*-1/2.211990302e24.gif (915 bytes)

?t;d0EVo0

(24)水利论文p:}W{,{^[ e{3]

在Cw=0.0025~0.0035范围内任假设一Cw值,由(24)式可得C′f值,此C′f值与步骤5之Cf值比较,若误差大于0.00005,则假设一新的Cw值,原则上C′f值>Cf则降低Cw假设值,C′f<Cf则增加Cw假设值。以RUN C1为例,C′f/Cf与Cw关系如图5所示。
7)计算Cb值,将Cw、Cf值代入(14)式,得Cb
8)计算θb

θb=U2*b/(s-1)gd50=CbU2(s-1)gd50水利论文0kbBGhx5]

(25)水利论文&o2T+E/iU*\

9)计算θg值,将θb代入(15)式,得θg
10)计算Cg

θg=CgU2/(s-1)gd50

(]|`(x*?7]h3H%h0

(26)水利论文)O7Km{5_&`1t(kRY

由(25)式及(26)式知

Cg=θg/θb*Cb水利论文 ^6]4[ RL]yL

(27)水利论文0T;z^bO9g@n/z9Z

 11)检核步骤10之C′g与步骤1之Cg值是否误差在0.00005之内,若误差大于0.00005则重新假设Cg值并重复步骤2~步骤11。原则上C′g>Cg则增加步骤1之Cg假设值,C′g<Cg则降低Cg假设值。以RUN C1为例,Cg~C′g关系如图6所示。

990302t5.gif (2033 bytes)

9Rc3G? A8w0

图5 Cw与C′f/Cf关系(Run C1)水利论文O [&@ w[

Relationship between Cw and C′f/Cf(Run C1)水利论文9@ I5hL[

990302t6.gif (2157 bytes)

O&Zc6C Z/f5["t0

图6 Cg与C′g的关系(Run C1)

9m"|"^ ~6wx0

Relationship between Cgand C′g(Run C1)

r+J-e wi8@0

3.3 床形水利论文1MEu'}yqsum+wZ)v

    剪应力校正后的各项水力参数列于表3,输沙参数列于表4,对照表4θg与(15)式,可知本文中六组属平床与沙纹,其余七组属沙垄,θg最大值为0.134,而文献[10]Hotchkiss的θg值在0.108与0.271间,也是沙垄情况,但其输沙量较大。所以本文与文献[10]Hotchkiss的试验涵盖低水流区的平床、沙纹与沙垄,但未包括沙垄中输沙量很大的情况。

[1Z b_ rYH7md*]0

表3 各组试验之水力参数

&D:k{$\"F(Ic,T0

Hydraulic parameters for all experimental cases水利论文/h*U6@#Z]


Run

dk GG p yn#|+r&J0

(1)

ZL*W7^ dkL.? ZY0
r
(cm)

(2)

4Qg{.Cy w+u{k{2_0

U*

A7rp%[-c+R }rl_0

(cm/s)水利论文A5l/N'Dmc6g0o

(3)

@C _ IY;~ q6t0

rg水利论文p%T)@I N

(cm)水利论文 by(q+Ef+fqGS&EE-Z

(4)水利论文T$S2LtA3yh

U*g水利论文;W'{!?:Q_N/E\ V"rk

(cm/s)水利论文}.A1{(QBFYd

(5)水利论文8zK,a2Bz

Sf

Rx*T j ^ XY7h0

(6)

gO AV)D F}&lu&s%K0

Cf水利论文X N_+{;`]7nO o

(7)

;bm)q;dk5o0

Cb水利论文'i}0XvOM4R(Jc

(8)水利论文 v3_E j]`7s_

Cg

&lq oh"}7f/D?0

(9)

^0^U3hnB~z-M0

Cw

q"Q!jI$A F I0

(10)

^ @z1i;w0

ks

COMGB*G4]!l0

(11)

F-jAa[#vO-tf'm0

C13.1003.5901.4142.4240.004240.008230.011650.003780.003240.83
C22.9074.1341.2722.7340.006000.008760.008800.003810.003260.94
C33.3851.7973.5281.8350.000970.002790.002930.002930.002670.63
C43.2861.9023.4201.9410.001120.002820.002950.002950.002660.66
C53.8222.5804.7442.8740.001780.003370.004110.004110.002563.23
C63.3945.9411.8174.3470.010610.011240.017810.005970.003165.02
C73.3172.7804.0873.0850.002380.003510.004300.004300.002643.57
C83.5332.4374.3472.7030.001720.003410.004170.004170.002653.12
C93.3302.5663.7802.7340.002020.003040.003070.003500.002571.79
C103.5522.2304.0602.3840.001430.002890.003320.003320.002581.56
C113.2796.7191.6314.7380.014050.012710.019930.006260.003225.48
C123.4455.2901.7063.7230.008290.009400.014750.004590.003062.44
C132.7894.4791.2072.9470.007340.008990.012440.003880.003261.01

表4 各组试验之输沙参数水利论文a]!b#z@:Av&}

Parameters of sediment transport for all experimental cases

2zyN5H }0

Run水利论文Jn_ B8lG

(1)

5BZ8_p8H'^0

Qs(g/s)

1R V4C;T.]fGk*`0

(3)

2i f,H:w CU1J{0

q*s水利论文 `I2]4t s c},e

(4)水利论文4H:M'p)X3_&I

Re*

oL A;@;y Y!LV0

(5)水利论文$kNv#t G\ y

θb水利论文[2W Gb-m?

(6)

,Y([p2A TG8\0

θg水利论文|%yZaDF+qO

(7)

+I`3jT n!D v-[4X0

C11.670.1319.620.2850.093
C23.330.26110.860.3730.116
C30.830.0657.280.0530.053
C41.110.0877.700.0590.059
C51.670.02138.540.0380.038
C63.330.04258.290.2590.087
C71.670.02141.370.0430.043
C80.830.01036.250.0330.033
C93.330.09920.750.0610.061
C101.670.05018.090.0460.046
C115.830.07463.540.3280.103
C125.830.17328.250.3580.111
C135.830.45711.700.4310.134

4 分析

"B3e]2\ vU2v:~0

4.1 边壁摩擦因子的特性

+T9bLekvU!h%]F0

分析边壁摩擦因子Cw值与表3及表4中各参数的关系,发现存在临界的希尔无因次剪力θgc水利论文&bm:V0J)hB

θg<θgc时Cw=0.00262±0.00005

MVu3D*m:_0

(28)

s v&~/a b#AlC _0

θg>θgc时Cw=0.00320±0.00008水利论文6L(};ab0xff?

(29)

+k k0Y5BMAZ#j!b|;Y0

Cw与θg之关系如图7所示。本研究之数据显示θgc=0.061~0.087,上式提供了以试误法计算Cw值的一个依据。水利论文5p C$QH"w"DJpC/p

4.2 三角洲河道之输沙能力

8] C9biw_u0

    本研究采用以剪应力为基础的输沙公式来探讨三角洲河道之输沙能力,从剪应力观点建立的输沙公式可明显区分为两类,一类不考虑起动临界剪应力,如Engelund Hansen公式

p'b s6j.Cg4B0

Cgq*s=aθbg

k8Z4F]5?*Q-i$S[3h:BCTw0

(30)

:y+g z j:I3iP0
上式中a=0.05,b=2.5,q*s=qs/990302e31.gif (1048 bytes),qs单位是L2/T。水利论文i;a6X@wzy

    另一类考虑起动临界剪应力,如Meyer Peter and Muller公式水利论文3Vv}1lgd m[$wj

q*s=8(θg-0.047)1.5

Xf9w0A4G!d$[0

(31)

i&g\:d)O? };D1@Mo0

990302t7.gif (1740 bytes)

p:r P bN3p#N A0

图7 Cw与θg关系水利论文5};m0sJ!b\3z

Relationship between Cw and θg水利论文Y?2a%}.ID#h

在希尔图上需Re*>500无因次起动剪应力才不变,Meyer-Peter and Muller公式采用无因次起动剪应力为0.047,所以Meyer-Peter and Muller公式适用Re*>500的情况,本研究Re*=7.2~63.5

TAG: 水库 三角洲 李鸿源 俞维升 河道
顶:12 踩:19
【已经有87人表态】
17票
极差
12票
很差
7票
较差
10票
稍差
13票
稍好
11票
较好
8票
很好
9票
极好
下一篇:三峡工程引航道冲沙流量的合理规模(周建军)
上一篇:渭河下游近期水沙特性及冲淤规律(张翠萍 张原锋 高际萍)
查看全部回复【已有0位网友发表了看法】