输沙强度计算公式的统一.docxVIP

输沙强度计算公式的统一 1 输沙强度的预测 自法国杜沃斯（1879）提出的第一个过程沙输沙率公式以来，已得出100多年的历史。据不完全统计,目前各种输沙率公式已超过 100 个,但由于问题的复杂性,至今还没有得到一个能精确预报输沙率的公认的一般计算公式,根据不同的公式计算出的结果与实测资料相比往往有较大的出入。所以说输沙率问题的研究还远远未达到令人满意的程度,因而还有待进一步研究。 建立输沙率公式的基础是选择合适的水流强度指标。现有输沙率公式中的水流强度指标各式各样,归纳起来主要有三类,即切应力、流速与功率,而每一类又有若干不同的表示形式。不同的水流强度指标,究竟何者更适合预测输沙率,目前还存在很大争议。如对于推移质泥沙运动,Thompson 则认为切应力不如流速更能预报推移质输沙率;Carson 则认为利用沙粒切应力预报推移质输沙率比流速、功率等指标更佳。又如对于悬移质泥沙运动,Samaga 等认为总切应力是决定悬移质输沙率的指标;Swamee-Ojha 则提出沙粒切应力比总切应力更适于悬移质输沙率计算。Yang 则多次指出,无论是推移质还是悬移质泥沙运动,单位水流功率预测输沙强度要比平均流速、切应力、水流功率等指标优越。因此要研究输沙率,必须先研究水流强度指标。本文的目的有三,一是对于各种预测输沙强度的水流强度指标进行归纳,在此基础上总结出一个统一的水流强度指标形式;二是基于统一的水流强度指标,建立起一个一般形式的输沙强度计算公式;三是利用水槽资料对此一般形式的输沙强度计算公式进行率定,并与其它公式比较。 2 水流强度指数的分析 2.1 水资源分布及水流强度 现有输沙率公式中的水流强度指标各式各样,归纳起来主要有三类,即切应力、流速与功率,而每一类又有若干不同的表示形式。切应力类指标有多种具体的形式,早期的输沙公式大多是利用床面总切应力作为水流强度指标。总切应力的定义为: τ=γhJ(1)τ=γhJ(1) 式中,τ为总切应力,h为水深 (若有侧壁影响时,应取床面水力半径,下同),J为能坡,γ为水的容重。后来越来越多的学者认识到只有沙粒切应力才对泥沙运动起作用。若沙粒切应力用 Manning-Strickler 公式表述,则可分别得到基于水力半径分割法和能坡分割法的沙粒切应力: τ′1=γh′J=[√ρV7.68(h/d)1/6]1.5(γhJ)0.25(2)τ′1=γh′J=[ρ√V7.68(h/d)1/6]1.5(γhJ)0.25(2) τ′2=γhJ′=[√ρV7.68(h/d)1/6]2(3)τ′2=γhJ′=[ρ√V7.68(h/d)1/6]2(3) 式中,τ′1、τ′2分别为基于水力半径分割法的沙粒切应力及基于能坡分割法的沙粒切应力,h′为与沙粒阻力相对应的水深,J′为与沙粒阻力对应的能坡,ρ为水的密度,V为平均流速,d为泥沙粒径。 流速类水流强度也有多种形式,如近底流速、水面流速等。在实际问题中,用得更多的是平均流速V这一水流强度指标。前苏联及一些中国学者在研究输沙率问题时,利用了基于平均流速的 Velikanov 参数作为预测输沙浓度的水流强度指标。Velikanov 参数的定义为 R=V3ghω(4)R=V3ghω(4) 式中ω为泥沙沉速,g 为重力加速度。实际上,佛汝德数也是一个流速类水流强度指标: Fr=V(gh)12(5) 功率类水流强度指标:设长度为L的均匀流河段,其单位时间内所消耗的水流总功率为W=γQJL,基于此总水流功率,可导出两个功率类水流强度指标。一是 Bagnold 提出的单位面积水体水流功率(简称水流功率): w=WbL=γQJLbL=γhJV=τV(6) 式中,Q为流量;b为河宽。二是 Yang 建立的单位重量水体水流功率(简称单位水流功率): Ρ=WγAL=γQJLγAL=VJ(7) 2.2 水流强度指标t 由上述对各现有的水流强度指标的分析,我们发现不管是哪一个指标,最终由三项组成,一是水深,二是平均流速,三是能坡,因此我们可以构造出如下统一的水流强度形式: Τ=a1ha2Va3Ja4(8) 式中,T为统一的水流强度指标;a1为系数,该系数仅与泥沙粒径有关,a2、a3、a4、为指数。 若式 (8) 中的a1、a2、a3、a4取不同值时,可得上述各水流强度指标,如表 1 所示。这些水流强度指标系数与指数的变化范围还较大,主要是由于这些系数、指数均为各研究者根据某种理论或假定得出,并未得到实测资料的直接验证与支持,从而基于这些水流强度指标的输沙公式难免会产生一些与实际不完全符合的计算结果。 3 指数函数关系 输沙公式的本质在于建立输沙率gt与水流强度指标T之间的关系。一般均采用下列指数形式的函数关系: gt=AΤB(9) 式中A为系数;B为指数。由于水流强度指标与输沙率的量纲不尽相