水库三角洲河道输沙之研究(俞维升 李鸿源)
摘 要水库进口附近三角洲河道之水流系近似均匀流,并非真正之均匀流。三角洲河道动床之自我调整作用,因水库水位之回水影响,其输沙特性有异于明渠均匀流。传统上,将明渠均匀流平衡条件下之输沙公式应用于估算三角洲河道之输沙能力,似有其不确定性。本研究进行一系列低水流区之输沙试验,并纳入文献[10]Hotchkiss 之试验数据一并分析。结果显示,在水库之回水影响下,河床自我调整作用使其输沙能力大于假设明渠均匀流之情况,其中Engelund Hansen公式低估低水流区之输沙能力达10倍;Parker公式亦低估,尤其是无因次颗粒剪应力小于0.08时差误更大。依据试验成果,本研究并建议一经验公式,以供估算水库三角洲河道之输沙能力。
关键词三角洲 泥沙运动 边壁校正 水库淤积
1 前言
泥沙随着水流进入水库后,由于坝体的壅水,向下游水深渐增、流速渐减,水流挟沙能力因而降低,造成粗粒泥沙呈现超载状态并发生淤积,逐渐形成三角洲。依据三角洲床面的坡度变化,自上游往下游可将三角洲区分为尾部段、顶坡段及前坡段,顶坡段的渠道称为三角洲河道,能够到达三角洲前坡段下游的泥沙一般皆为细颗粒,水库三角洲纵向剖面示意如图1。
|
|
归纳以往模拟水库顶坡段泥沙冲淤的数值模式,计算输沙量的方法有二种,一种是针对河床推移质采用平衡条件下的输沙公式,针对悬移质以泥沙扩散方程式求解[1~3],其中泥沙扩散方程式可反映悬移质的不平衡输沙行为。另一种是采用平衡条件下的总输沙量公式[4,5]。
明渠水流输沙公式的数据来源,不外乎现场实测及水槽试验量测,通常要求量测数据时的水流条件是均匀流,量测段必须达到水面平行底床,水深是正常水深,在上述均匀流条件下,运动泥沙与床面泥沙的交换达到平衡。严谨而言,现有输沙公式大多只适用于均匀流且泥沙冲淤平衡的情况。
Hotchkiss and Parker指出在明渠中淤积造成底床升高的情况有两种[6],一种是上游来沙量大于该断面输送能力所造成,此情况下,该断面自我调整增加河床坡度以增加输沙能力,底床升高过程示意如图2,另一种是下游回水效应使该断面水流输沙能力降低,此情况下,该断面为增加输沙能力的自我调整作用受回水影响,在三角洲河道即为此种情况,第二种情况之底床升高如图1所示,明显与第一种情况不同。第一种情况,在底床升高的过程水流仍为均匀流,仍可引用平衡输沙公式;第二种情况,水流非均匀流且受回水影响,其输沙特性为本文研究重点.Hotchkiss and Parker指出三角洲河道仅是近似均匀流[6](near uniform flow),并非真正的均匀流,本文认为以往将平衡条件的输沙公式应用于计算三角洲河道输沙量,存在不确定性。
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的沙垄数据一并纳入,以分析三角洲河道的输沙特性,试验成果并与平衡输沙公式比较,进而分析前人以平衡输沙公式计算三角洲河道输沙量的误差。
2 试验设备与方法
本研究在预备试验时,曾进行如图1般底床有一固定坡度的水槽试验,检验结果发现尾部段的溯源淤积将使三角洲河道的来沙量无法维持一固定值。所以本研究设计如图3的水槽,该水槽长3m,其中1m是明渠段、2m是水库段,宽15cm,底床坡度接近水平。试验水槽在明渠段与水库段交界处有一落差10cm的阶梯,试验时间内控制明渠段无淤沙,如此可确保试验时间内三角洲河道有固定的输沙量。为使明渠段水流稳定进入水库段,控制顶坡段淤沙高程接近明渠段底床高程。
本研究试验泥沙有三种,平均粒径(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个位置的高程。 |
|
表1 试验泥沙的粒径特性 |
Particle properties of experimental material(单位mm) |
项目 | dm | d50 | d84.1 | d15.9 | σg |
第一种泥沙 | 1.3341 | 1.330 | 1.610 | 1.050 | 1.238 |
第二种泥沙 | 0.759 | 0.739 | 0.878 | 0.632 | 1.179 |
第三种泥沙 | 0.397 | 0.377 | 0.500 | 0.306 | 1.278 |
表2 试验条件与成果 |
Experimental conditions and results |
|
|
|
|
|
|
| |||||||||||||||||||||
C1 | 3141.9 | 1.67 | 0.000200 | 0.0377 | 5.284 | 0.0069 | |||||||||||||||||||||
C2 | 3141.9 | 3.33 | 0.000398 | 0.0377 | 4.747 | 0.0062 | |||||||||||||||||||||
C3 | 3141.9 | 0.83 | 0.000099 | 0.0377 | 6.169 | 0 | |||||||||||||||||||||
C4 | 3141.9 | 1.11 | 0.000133 | 0.0377 | 5.847 | 0.0015 | |||||||||||||||||||||
C5 | 5217.0 | 1.67 | 0.000120 | 0.1330 | 7.795 | 0.0043 | |||||||||||||||||||||
C6 | 4168.5 | 3.33 | 0.000300 | 0.1330 | 6.198 | 0.0064 | |||||||||||||||||||||
C7 | 4168.5 | 1.67 | 0.000150 | 1.1330 | 5.948 | 0.0045 | |||||||||||||||||||||
C8 | 5217.0 | 0.83 | 0.000060 | 0.1330 | 6.679 | 0.0077 | |||||||||||||||||||||
C9 | 4168.5 | 3.33 | 0.000300 | 0.0739 | 5.988 | 0.0080 | |||||||||||||||||||||
C10 | 4168.5 | 1.67 | 0.000150 | 0.0739 | 6.748 | 0.0223 | |||||||||||||||||||||
C11 | 5217.0 | 5.83 | 0.000420 | 0.1330 | 5.826 | 0.0154 | |||||||||||||||||||||
C12 | 5217.0 | 5.83 | 0.000420 | 0.0739 | 6.371 | 0.0058 | |||||||||||||||||||||
C13 | 3141.9 | 5.83 | 0.000697 | 0.0377 | 4.439 | 0.0079 | |||||||||||||||||||||
3 试验结果
3.1 量测数据 典型的淤积剖面如图4所示,图4是分别在不同时刻停水量测的淤积剖面,可看出淤积剖面很稳定,在三角洲顶点及顶坡段与明渠交界处(即图4中A与B位置)底床因受局部流场影响而较不稳定,但在位置A与B之间,底床坡度与水深接近不变,可量得顶坡段代表性之水深及底床坡度。 |
|
本文每个断面侧向量8个位置的底床高程,虽然宽深比约为3有边壁效应存在,但由所量测数据可知侧向上底床变化不大。另外,由柏努力方程式(Bernoulli equation)可知底床坡度与水深的斜率同级,所以量测值中水深明显较底床坡度有较高的可信度,亦即底床坡度很容易因量测技术及数据处理上的些微差异而产生较大的变化,本研究在初步分析中发现所量测的底床坡度不具规律性,所以本文量测之底床坡度仅作参考,未当作分析数据。实验结果水深H在4.747~6.748cm间,顶坡段底床坡度S0在0~1.54%间,示如表2。
3.2 计算影响泥沙运动的有效剪应力
由于水槽边壁是丙烯酸树脂(acrylic resin)制,底床是泥沙,不同材质使得剪应力在湿周上并非均匀分布,且水槽宽15cm,宽深比约为3,所以边壁与底床不同材质对泥沙输移的影响不可忽略,应进行边壁校正,从总应力τ中分离出与输沙有关的底床剪应力τb。此外,底床剪应力中克服床形起伏的形状阻力所消耗能量与河床质输沙量无关。仅克服泥沙颗粒表面摩擦阻力所消耗能量与河床质输沙量有关,所以需将底床剪应力中与河床质输沙量有关的粒径剪应力τg分离出来。本文边壁校正采用1942年Johnson提出的方法[11],分离粒径剪应力采用1952年Einstein and Barbarossa提出的方法[11],本文并提出校正剪应力的步骤,校正剪应力的流程详述如下:
水槽边壁是丙烯酸树脂制的光滑边壁,流过光滑面的水流摩擦因子,可由慕迪图(Moody diagram)中光滑紊流之摩擦因子与雷诺数关系曲线求得,Granger提出下式足以表现该关系曲线[13]
(1) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
式中 w为与边壁有关,Rew为雷诺数,Cw=fw/8,fw为达西-韦斯巴摩擦因子,Cw为修正的摩擦因子,剪应力τw=ρCwU2,其中ρ为水的密度,U为水流平均速度。 (1)式中Rew=4rwU/ν,rw是水力半径,ν是水的运动粘滞系数,因为rw为未知数所以无法直接求得Cw值。Johnson针对边壁校正作以下假设[11]
|
|
|
3.3 床形
剪应力校正后的各项水力参数列于表3,输沙参数列于表4,对照表4θg与(15)式,可知本文中六组属平床与沙纹,其余七组属沙垄,θg最大值为0.134,而文献[10]Hotchkiss的θg值在0.108与0.271间,也是沙垄情况,但其输沙量较大。所以本文与文献[10]Hotchkiss的试验涵盖低水流区的平床、沙纹与沙垄,但未包括沙垄中输沙量很大的情况。
表3 各组试验之水力参数 |
Hydraulic parameters for all experimental cases |
|
|
|
|
|
|
|
|
|
|
| |||||||||||||||||||||||||||||||||
C1 | 3.100 | 3.590 | 1.414 | 2.424 | 0.00424 | 0.00823 | 0.01165 | 0.00378 | 0.00324 | 0.83 | |||||||||||||||||||||||||||||||||
C2 | 2.907 | 4.134 | 1.272 | 2.734 | 0.00600 | 0.00876 | 0.00880 | 0.00381 | 0.00326 | 0.94 | |||||||||||||||||||||||||||||||||
C3 | 3.385 | 1.797 | 3.528 | 1.835 | 0.00097 | 0.00279 | 0.00293 | 0.00293 | 0.00267 | 0.63 | |||||||||||||||||||||||||||||||||
C4 | 3.286 | 1.902 | 3.420 | 1.941 | 0.00112 | 0.00282 | 0.00295 | 0.00295 | 0.00266 | 0.66 | |||||||||||||||||||||||||||||||||
C5 | 3.822 | 2.580 | 4.744 | 2.874 | 0.00178 | 0.00337 | 0.00411 | 0.00411 | 0.00256 | 3.23 | |||||||||||||||||||||||||||||||||
C6 | 3.394 | 5.941 | 1.817 | 4.347 | 0.01061 | 0.01124 | 0.01781 | 0.00597 | 0.00316 | 5.02 | |||||||||||||||||||||||||||||||||
C7 | 3.317 | 2.780 | 4.087 | 3.085 | 0.00238 | 0.00351 | 0.00430 | 0.00430 | 0.00264 | 3.57 | |||||||||||||||||||||||||||||||||
C8 | 3.533 | 2.437 | 4.347 | 2.703 | 0.00172 | 0.00341 | 0.00417 | 0.00417 | 0.00265 | 3.12 | |||||||||||||||||||||||||||||||||
C9 | 3.330 | 2.566 | 3.780 | 2.734 | 0.00202 | 0.00304 | 0.00307 | 0.00350 | 0.00257 | 1.79 | |||||||||||||||||||||||||||||||||
C10 | 3.552 | 2.230 | 4.060 | 2.384 | 0.00143 | 0.00289 | 0.00332 | 0.00332 | 0.00258 | 1.56 | |||||||||||||||||||||||||||||||||
C11 | 3.279 | 6.719 | 1.631 | 4.738 | 0.01405 | 0.01271 | 0.01993 | 0.00626 | 0.00322 | 5.48 | |||||||||||||||||||||||||||||||||
C12 | 3.445 | 5.290 | 1.706 | 3.723 | 0.00829 | 0.00940 | 0.01475 | 0.00459 | 0.00306 | 2.44 | |||||||||||||||||||||||||||||||||
C13 | 2.789 | 4.479 | 1.207 | 2.947 | 0.00734 | 0.00899 | 0.01244 | 0.00388 | 0.00326 | 1.01 | |||||||||||||||||||||||||||||||||
表4 各组试验之输沙参数 |
Parameters of sediment transport for all experimental cases |
|
|
|
|
|
| ||||||||||||
C1 | 1.67 | 0.131 | 9.62 | 0.285 | 0.093 | ||||||||||||
C2 | 3.33 | 0.261 | 10.86 | 0.373 | 0.116 | ||||||||||||
C3 | 0.83 | 0.065 | 7.28 | 0.053 | 0.053 | ||||||||||||
C4 | 1.11 | 0.087 | 7.70 | 0.059 | 0.059 | ||||||||||||
C5 | 1.67 | 0.021 | 38.54 | 0.038 | 0.038 | ||||||||||||
C6 | 3.33 | 0.042 | 58.29 | 0.259 | 0.087 | ||||||||||||
C7 | 1.67 | 0.021 | 41.37 | 0.043 | 0.043 | ||||||||||||
C8 | 0.83 | 0.010 | 36.25 | 0.033 | 0.033 | ||||||||||||
C9 | 3.33 | 0.099 | 20.75 | 0.061 | 0.061 | ||||||||||||
C10 | 1.67 | 0.050 | 18.09 | 0.046 | 0.046 | ||||||||||||
C11 | 5.83 | 0.074 | 63.54 | 0.328 | 0.103 | ||||||||||||
C12 | 5.83 | 0.173 | 28.25 | 0.358 | 0.111 | ||||||||||||
C13 | 5.83 | 0.457 | 11.70 | 0.431 | 0.134 | ||||||||||||
4 分析
4.1 边壁摩擦因子的特性
分析边壁摩擦因子Cw值与表3及表4中各参数的关系,发现存在临界的希尔无因次剪力θgc
θg<θgc时Cw=0.00262±0.00005 | (28) |
θg>θgc时Cw=0.00320±0.00008 | (29) |
Cw与θg之关系如图7所示。本研究之数据显示θgc=0.061~0.087,上式提供了以试误法计算Cw值的一个依据。
4.2 三角洲河道之输沙能力
本研究采用以剪应力为基础的输沙公式来探讨三角洲河道之输沙能力,从剪应力观点建立的输沙公式可明显区分为两类,一类不考虑起动临界剪应力,如Engelund Hansen公式
Cgq*s=aθbg | (30) |
上式中a=0.05,b=2.5,q*s=qs/![]() 另一类考虑起动临界剪应力,如Meyer Peter and Muller公式
|
|
在希尔图上需Re*>500无因次起动剪应力才不变,Meyer-Peter and Muller公式采用无因次起动剪应力为0.047,所以Meyer-Peter and Muller公式适用Re*>500的情况,本研究Re*=7.2~63.5