基于跌坎压力自相似分布的明渠水深−流量关系推求

李乃稳 黄滟淳 陈小攀 李龙国 刘超

李乃稳, 黄滟淳, 陈小攀, 等. 基于跌坎压力自相似分布的明渠水深−流量关系推求 [J]. 工程科学与技术, 2024, 56(1): 218-227. doi: 10.15961/j.jsuese.202201219
引用本文: 李乃稳, 黄滟淳, 陈小攀, 等. 基于跌坎压力自相似分布的明渠水深−流量关系推求 [J]. 工程科学与技术, 2024, 56(1): 218-227. doi: 10.15961/j.jsuese.202201219
LI Naiwen, HUANG Yanchun, CHEN Xiaopan, et al. Quasi-theoretical End-depth–Discharge Relationship Based on the Pressure Distribution at the End of Free Overfall for Rough Rectangle Channels [J]. Advanced Engineering Sciences, 2024, 56(1): 218-227. doi: 10.15961/j.jsuese.202201219
Citation: LI Naiwen, HUANG Yanchun, CHEN Xiaopan, et al. Quasi-theoretical End-depth–Discharge Relationship Based on the Pressure Distribution at the End of Free Overfall for Rough Rectangle Channels [J]. Advanced Engineering Sciences, 2024, 56(1): 218-227. doi: 10.15961/j.jsuese.202201219

基于跌坎压力自相似分布的明渠水深−流量关系推求

基金项目: 国家自然科学基金面上项目(52079084;52379073);四川省重大科技专项(2018SZDZX0027);四川省科技支撑计划项目(2022YFS0472;2023YFS0377);四川大学–泸州校地合作项目(2020CDLZ–6)
详细信息
    • 收稿日期:  2022-11-08
    • 网络出版时间:  2023-01-31 12:00:00
  • 作者简介:

    李乃稳(1975—),博士,副教授. 研究方向:水力学;水环境与水生态. E-mail:linaiwen@scu.edu.cn

    通信作者:

    刘超, E-mail:liuchao@scu.edu.cn

  • 中图分类号: TV131.4

Quasi-theoretical End-depth–Discharge Relationship Based on the Pressure Distribution at the End of Free Overfall for Rough Rectangle Channels

  • 摘要: 明渠流量计量是高效灌溉、合理分配水资源和进行用水有效管理的重要手段。跌水是常见的明渠水流上下游衔接的水利设施。当跌落水流不受跌坎下游水流影响时形成自由跌流,可用于明渠流量的精确计量。受自由跌流影响,跌坎断面水流由缓流向急流转变,压力分布显著偏离静水压力分布规律。而获得跌坎压力分布及其准确描述方程则正是明渠自由跌流流量公式推求的关键,其准确与否决定了利用跌流进行流量计量的精度。本文首先通过模型试验和数值计算对矩形明渠自由跌流的断面压力分布特征进行研究,而后利用理论推导建立了明渠跌流的水深−流量关系公式。研究发现,明渠跌水跌坎断面压力呈典型自相似分布规律,其无量纲压力值分布倾向于一条曲线,压力最大值为0.238he (he为跌坎水深),出现在竖向距跌坎0.209he水深处,并且压力自相似分布规律与明渠流量、渠道底坡和边壁糙率无关。在此基础上,结合幂函数和尾迹函数提出了描述跌坎断面压力分布的函数方程,并给出了最大压力和最大压力位置的计算公式,然后根据动量方程和流量连续方程推求了明渠自由跌流的流量计算公式。最后,对本文和其他研究者提出的流量计算公式采用试验数据进行了对比分析,结果表明:本文提出的流量公式在明渠流量为 5~100 L/s,底坡为–0.0112~0.0534、糙率为0.093 0~0.019 3的研究范围内,具有更高的精度和适用性,计算误差在±5%之内,可适用于该流量、底坡和边壁糙率范围内矩形明渠自由跌流的流量精确计算。本文研究结果对于利用自由跌流进行明渠流量计量提供了理论与技术支撑,具有重要的工程应用价值。

     

    Abstract: Flow measurement is an important means for irrigation, efficient allocation of water resources and water diversion management. In the open channel, the overfall is usually applied as a facility connecting the upstream and downstream. When the overfall is not affected by the flow of downstream, the overfall is classified as a free one, which can be used as a discharge measurement by a single measurement of depth at the end of the channel. At this condition, the flow pattern changes from a subcritical flow into a supercritical one, and the pressure distribution in the end section deviates greatly from the hydrostatic pressure distribution with a skewed shape, which is the key to derive the discharge formula by energy equation and momentum equation. In this study, the characteristics of sectional pressure distribution of free overfall in a rectangular open channel were studied by model test and numerical calculation, and then the water depth-flow relationship formula of open channel with free overfall was established. It was found that the pressure of the end section was conformed to the self-similar distribution law, as the distribution of dimensionless pressure values tended to be the same curve, in which the maximum value was 0.238he (here he was water depth of end section) and it occurred at 0.209he distancing from the end fall. The above findings are independent of the flow discharge, bottom slope and side wall roughness. Based on these findings, a simple formula combining a power function and wake function was proposed to describe the pressure distribution at the end section. Then using the momentum equation and continuity equation, a quasi-theoretical end-depth–discharge relationship was proposed. With the known end depth, bed slope and Manning’s coefficient, the discharge could be direct solved by the proposed equations. The proposed equation was valid over the entire practical ranges of Q=5~100 L/s, S=–0.011 2~0.053 4, n=0.093 0~0.019 3, and they showed excellent agreement with the available experimental data in this paper and from others, with the precision of ±5%. The proposed formula for discharge is a useful computational tool for estimation of discharge in rough bed rectangular open channels. The results of this paper provide theoretical and technical supports for the open channel discharge measurement using free over flow, which is important to the application of free overfall as a discharge measure device.

     

  • 渠道量水是灌溉、水资源高效配置及引供水用水管理的重要手段。目前常用的渠道量水设施有如下几种[1]:利用水工建筑物量水,如渠槽、堰闸、渡槽等;利用流速仪量水和特设量水设备量水,如各种量水堰(矩形堰、三角堰等)或量水槽(巴歇尔量水槽)等。这些量水设施存在要么影响渠道过流能力、流量系数确定难度大,要么通过泥沙、漂浮物能力差,或者量测精度低等问题。因此,寻求一种结构简单、测量精确的渠道量水设施则成为目前节水社会建设背景下明渠输水和量水技术研究及发展的重点目标。

    在明渠中,跌流为常用的衔接上下游水利设施。当跌落水流不受跌坎下游水流影响时,定义此跌流为自由跌流,此时存在急流到缓流的流态转变,跌坎断面压力分布与静水压力分布规律之间存在巨大偏差[23]。但此时,跌坎水深与明渠流量之间存在直接对应关系,这种对应关系使得自由跌流可以用作明渠流量的计量。早在1936年,Rouse[4]就认识到了这种自由跌流的水力特征并提出了坎末水深比的概念,即跌坎水深与明渠临界水深之比,通过该值可计算渠道流量。之后,研究者为寻求准确适用的跌流流量计算方法,从理论分析、试验研究两个方面进行了大量研究 [5],提出了理论或经验或半经验流量计算公式[67]。基于自由跌流和薄壁堰流二者之间的相似性,一些研究者根据能量方程或动量方程 [89],提出了底坡为平坡和缓坡的近似理论计算方程。不同研究者提出的流量方程不同之处在于其在推求方程时,针对跌坎压力分布不再符合静水压力分布规律的问题,对能量方程或动量方程中跌坎断面竖向压力分布项采用不同的处理方式[10],如基于跌坎断面零压力分布[1112],跌坎断面水流流线曲率的Boussinesq近似规律 [1315],跌坎压力为上游平均流断面压力分布的折减系数[11,16],跌坎断面压力二次抛物线分布 [3, 1718]等。

    显然,跌坎断面压力分布是理论或半理论推求其流量方程的关键,假定与实际压力分布之间偏差大小势必会造成流量计算值与实测值存在较大偏差,进而影响利用自由跌流对明渠流量进行计量。

    本文利用试验和数值计算相结合的方法,首先利用动量方程,推求了明渠自由跌坎水深和流量之间的理论关系。之后,针对自由跌流水深−流量理论关系推求中的跌坎压力分布这个关键问题,在前期明渠自由跌流水力特性研究的基础上[19],根据试验和数值计算,提出了其压力分布自相似规律方程。在此基础上,提出了明渠跌流流量的理论计算方程,并利用本文和其他研究者实测数据与其他研究者提出的流量计算方程进行了比较,结果表明,本文方程计算更精确,并且适用性更广,能够为利用渠道跌坎进行流量计算提供理论依据和方法。

    水流在重力作用下在跌坎处跌落,形成自由跌落水流(图1(a)),其流态类似于薄壁堰流,水面线为典型的降水曲线,存在从缓流到急流、均匀流到非均匀流的流态转变,在跌坎临近出现临界水流。在跌坎处,水舌上方和下方均为大气压,即压力P=0,跌坎断面压力明显偏离明渠均匀流的静水压力分布规律,压力明显小于静水压力值 [13]

    图  1  矩形明渠自由跌坎过流示意图
    Fig.  1  Schematic view of a typical free overfall
    下载: 全尺寸图片

    图1所示,O−O断面水深为h0,跌坎e−e断面水深为heh为水深,u为流速,P为断面压力,py为断面水深y处的压强。以O−O与e−e断面及沿渠底与水流自由表面所围成的空间作为控制体,将该段明渠水流看作1维水流,则沿水流方向的动量方程如下:

    $$ {P_{\text{o}}}\cos \;\phi - {P_{\text{e}}}\cos \;\phi - [{F_{\mathrm{b}}}]_{{{{{x}} = - {{L}}}}}^{{{{{x}} = 0}}} + W\sin \;\phi = \rho Q({\beta _{\text{e}}}{v_{\text{e}}} - {\beta _{\text{o}}}{v_{\text{o}}}) $$ (1)

    式中:Po为O−O断面上总动水压力,Pa;Pe为e−e断面上总动水压力,Pa;$ [{F_{\mathrm{b}}}]_{{{{{x}} = - {{L}}}}}^{{{{{x}} = 0}}} $x=–Lx=0之间渠道对水体的摩擦阻力,N;W为水体自重,N;φ为渠道底坡与水平面夹角,对明渠而言,底坡$ S = \tan \;\phi $,则$ \cos \;\phi = 1/\sqrt {1 + {S^2}} $$ \sin \;\phi = S/\sqrt {1 + {S^2}} $;ρ为水体密度,kg/m3Q为流量,m3/s;vo为O−O断面平均流速,m/s;ve为e−e断面平均流速,m/s;β为动量修正系数,βo = βe = 1。

    断面O−O水压力分布为静水压力分布,表达式为:

    $$ {P_0} = \frac{1}{2}\rho gbh_0^2 $$ (2)

    式中:g为重力加速度,取9.81 m/s2b为矩形渠道断面宽度,m。

    断面O−O(x = –L)到跌水断面e−e(x = 0)之间的总摩擦阻力为:

    $$ \left[ {{F_{{{\mathrm{b}}}}}} \right]_{{{{{x}} = - {{L}}}}}^{{{{{x}} = 0}}} = b\int_{{{ - {{L}}}}}^{\text{0}} \tau {\text{d}}x $$ (3)

    式中:τ为控制体水体所受边界切应力,$ \tau = \dfrac{f}{8}\rho {v^2} $,其中,v为控制体水流流速,计算中取断面O−O断面和跌坎断面e−e两个过水断面的平均流速,即$ v = \dfrac{{{v_{\text{o}}} + {v_{\text{e}}}}}{2} $f为摩擦因数,$f = {\left( {2.344 - 2\lg \left( {{k_{\text{s}}}/h} \right)} \right)^{ - 2}}$(Kraijenhoff等[20]),h为控制体的水深, m,计算中取两个计算控制断面的平均水深,即$ h = \dfrac{{{h_{\text{o}}} + {h_{\text{e}}}}}{2} $ks为尼古拉兹等效砂粗糙度。则由式(3)可得:

    $$ \left[ {{F_{\text{b}}}} \right]_{{{{{x}} = - {{L}}}}}^{{{{{x}} = 0}}} = \frac{{bL\rho {{\left( {{v_0} + {v_{\text{e}}}} \right)}^2}}}{{32{{\left[ {2.344 - 2\lg \left( {\dfrac{{2{n^6}}}{{{{0.038\;2}^6}\left( {{h_{\text{0}}} + {h_{\text{e}}}} \right)}}} \right)} \right]}^2}}} $$ (4)

    式中,n为渠道边壁粗糙度系数,$n = 0.038\;2k_{\text{s}}^{{1/ 6}}$(Guo等 [21])。

    断面O−O和跌坎断面e−e之间的水体自重为:

    $$ W=\rho gb{\displaystyle \int _{-L}^{0}h\text{d}x}=\frac{1}{2}\rho gb({h}_{0}+{h}_{\text{e}})L $$ (5)

    断面O−O和跌坎断面e−e之间的距离L,可由式(6)计算[19]

    $$ L = 3.0{\left( {\frac{{\sqrt S }}{n}} \right)^{ - 0.80}} $$ (6)

    如明渠底坡SScSc为临界底坡,则在跌坎上游存在均匀流。在明渠均匀流中,有:

    $$ Q = AC\sqrt {RJ} $$ (7)

    式中:A为过水断面面积,$ A = b{h_{\text{o}}} $C为谢齐系数,$ C = \dfrac{1}{n}{R^{1/6}} $R为水力半径,$ R = \dfrac{{b{h_{\text{o}}}}}{{b + 2{h_{\text{o}}}}} $J为水力坡降,J = S;因此,

    $$ Q = \frac{{b{h_{\text{o}}}}}{n}{\left( {\frac{{b{h_{\text{o}}}}}{{b + 2{h_{\text{o}}}}}} \right)^{2/3}}{S^{1/2}} $$ (8)

    由以上公式推导而知,对于式(1),关键是获得跌坎断面e−e的压力分布Pe

    图1中,跌坎断面e−e的压力分布不符合静水压力分布。假定跌坎处压力为二次抛物线分布[13],对称轴为y = 0.5he,压力最大值为khe,其中k为压力系数,其值与渠道断面形状有关。

    结合边界条件,如在点y = 0 和ye = he处,有py = 0,则跌坎断面处的压力分布方程为:

    $$ {p_y} = - \frac{{4k}}{{{h_{\text{e}}}}}{{{y}}^2} + 4ky $$ (9)

    式中,y为距离跌坎边缘的距离,py 为点y处的压力。

    从式(9)可知,y = 0.5he时,即在0.5倍跌坎水深处,py达到最大值,并且有压力最大值pmax = khe。根据文献[13]研究可知,对于矩形明渠而言,k = 0.281。图1(a)中点化线即为由式(9)对应的压力分布。事实上,文献[2]试验和文献[2223]数值计算表明,跌坎压力分布为严重扭曲形态,而非二次抛物线形态,其跌坎压力最大值点也不是在0.5倍跌坎水深处。

    采用模型试验和数值计算相结合,本文对不同底坡和边壁糙率的矩形明渠自由跌流跌坎断面压力分布进行了详细研究,其试验和数值模型如图2所示,试验模型由钢化玻璃制作,并采用粒径分别为0.5~1.0 mm,1.0~2.0 mm,2.0~3.0 mm和3.0~5.0 mm的沙粒黏附于渠道边壁而获得不同的渠道曼宁糙率n值,经率定[24],其糙率系数对应为n = 0.009 3,0.012 5,0.014 8,0.017 1,0.019 3;试验中渠道底坡S = –0.011 2~0.053 4,Q=5~100 L/s,共计152个工况。

    图  2  试验装置
    Fig.  2  Experimental installation
    下载: 全尺寸图片

    数值模型采用VOF结合RANS kε紊流模型[2526],基于紊流平衡假设,可以得到近壁处的紊动能和耗散率值 [2223]。模型中尼古拉斯等效砂高度ks由经验公式$ n = 0.038\;2k_{\text{s}}^{1/6} $给出 [21]

    在跌坎附近不同断面的压力分布如图36所示。图36中,黑实线为式(9)计算得出的跌坎断面压力分布,x为跌坎断面e−e坐标0点时的横向坐标值(见图1(a))。图4(a)图6(a)中,散点数据为模型试验数据。

    图  3  不同断面压力分布(S = –0.011 2)
    Fig.  3  Pressure distribution at section of S = –0.011 2
    下载: 全尺寸图片
    图  4  不同断面压力分布 (S = 0)
    Fig.  4  Pressure distribution at section S = 0
    下载: 全尺寸图片
    图  5  不同断面压力分布 (S = 0.011 2)
    Fig.  5  Pressure distribution at section S = 0.011 2
    下载: 全尺寸图片
    图  6  不同断面压力分布 (S = 0.024)
    Fig.  6  Pressure distribution at section S = 0.024
    下载: 全尺寸图片

    图36可知,底坡S和流量Q不同时,从跌坎断面至上游各断面,对应断面压力分布形状几乎一致。随着水流临近跌坎,受自由跌流影响,压力逐渐偏离静水压强分布。而在均匀流中,各断面压力分布符合静水压强分布规律。最终,在跌坎断面x = 0 m处,压力分布完全不同于静水压力分布而呈现扭曲形态。从水面至最大压力点,压力从0逐渐增加至最大值,而后急剧减小至0,其分布形态也不同于式(9)所描述的二次抛物线形状。在大约0.20he处,压力达到最大值,而根据式(9),压力最大值应该出现在0.5he处。

    对底坡S = –0.0112和S = 0的明渠,由式(9)计算的压力最大值与数值计算值吻合良好,并且总压力值也与模拟计算值一致。但对于陡坡渠道(如S = 0.0112和0.0240),随着渠道底坡的增加,由式(9)计算的最大压力值与数值计算的压力最大值之间的偏差迅速增加。例如,在渠道底坡S=0和流量Q = 80.00 L/s 时,式(9)计算的压力最大值为393.76 Pa,而数值计算的压力最大值为372.71 Pa,二者偏差仅为5.3%。但当渠道底坡S = 0.011 2 和流量Q = 80.00 L/s 时,式(9)计算的压力最大值为343.43 Pa,而数值计算的压力最大值则为265.61 Pa,二者偏差为22.7%。当渠道底坡S = 0.0240和流量Q = 80 L/s 时,式(9)计算的压力最大值为273.76 Pa,数值计算的压力最大值则为182.75 Pa,二者偏差达到33.2%。

    因跌坎断面压力分布存在相似性,本文引入无量纲数$\hat p = {p_{\text{y}}}/{p_{{\text{max}}}}$$\theta = y/{h_{\text{e}}}$,对跌坎断面压力分布进行无量纲化。将$ \theta $作为X轴,$ \hat p $作为Y轴,则$ \hat p $$ \theta $之间关系曲线如图7所示。

    图  7  $ \hat {\boldsymbol{p}} $$ {\boldsymbol{\theta}} $断面分布
    Fig.  7  $ \hat {\boldsymbol{p}} $${\boldsymbol{ \theta}} $ distribution on the end section
    下载: 全尺寸图片

    图7可知,底坡S和流量Q都不同的矩形渠道,跌坎断面压力分布呈现非常良好的自相似特征,无量纲压力$ \hat p $倾向于为一条实线分布。在图7中存在3个特殊点,即点A、B、和M,其分别对应$ \hat{p} $(θ=1)=0, $ \hat{p} $(θ=0)=0和$\hat p$(θ=θm)=1,此处有${\theta _{\text{m}}} = {y_{\text{m}}}/{h_{\text{e}}}$。也就是说,在点M处(y = ym),压力py达到最大值,此时$ \hat{p} $(θm=1)。

    引入由幂函数$ {\theta ^\xi } $和尾迹函数$ {\left( {1 - \theta } \right)^\gamma } $组成的混合函数来表征图7中的无量纲压力分布,并假定存在如下函数形式:

    $$ \hat{p}(\theta )=\sigma {\theta }^{\xi }{\left(1-\theta \right)}^{\gamma } $$ (10)

    式中,$ \sigma $为系数,$ \xi $$ \gamma $为指数系数。这些系数的值可以由图7中的数据确定。

    对点M,$ \hat{p} $(θ=θm)=1, 由式(10)可得:

    $$ \sigma = \theta _{\mathrm{m}}^{ - \xi }{(1 - {\theta _{\mathrm{m}}})^{ - \gamma }} $$ (11)

    对式(10)求解$ \hat p $关于$ \theta $的微分,可得:

    $$ \frac{{{\kern 1pt} {\text{d}}\hat p}}{{{\text{d}}\theta }} = \sigma \xi {\theta ^{{{\xi - 1}}}}{\left( {1 - \theta } \right)^{\gamma}} - \sigma {\theta ^{\xi}}\gamma {\left( {1 - \theta } \right)^{{\gamma-}1}} $$ (12)

    在点M (${\theta _{\text{m}}} = {y_{\text{m}}}/{h_{\text{e}}}$)出现压力最大值,即有$ \hat{p}(\theta = {\theta }_{\text{m}})=1 $,式(10)存在极值。因此,有$\dfrac{{{\text{d}}\hat p}}{{{\text{d}}\theta }}\left( {\theta = {\theta _{\text{m}}}} \right) = 0$。应用此条件,由式(12)可得:

    $$ {\theta _{\mathrm{m}}} = \frac{\xi }{{\xi + \gamma }} $$ (13)

    将式(13)代入式(11),可得:

    $$ \sigma = \frac{{{{(\xi + \gamma )}^{\xi + \gamma }}}}{{{\xi ^\xi }{\gamma ^\gamma }}} $$ (14)

    从式(14)可知,系数$ \sigma $的值与指数$ \xi $$ \gamma $有关,可由图7中的$ \hat p $数据确定。当指数$ \xi = 0.36 $$ \gamma = 1.36 $时,式(10)计算值与图7中数据吻合最好。因此,由式(14)可得$ \sigma = 2.417 $

    将此3个系数值代入式(10),可得:

    $$ \hat p = 2.417{\theta ^{0.36}}{\left( {1 - \theta } \right)^{1.36}} $$ (15)

    即:

    $$ \frac{{{p_{\text{y}}}}}{{{p_{\max }}}} = 2.417{\left(\frac{y}{{{h_{\text{e}}}}}\right)^{0.36}}{\left( {1 - \frac{y}{{{h_{\text{e}}}}}} \right)^{1.36}} $$ (16)

    利用式(15),可得图7中方程与数据的标准偏差、标准误差和相关系数分别为0.325、0.044和0.991。最大压力点为${\theta _{\text{m}}} = 0.209$,此值与图7数据得值非常一致。由图7数据可知,当${\theta _{\text{m}}} = 0.202$时,压力达到最大值。

    由式(15)计算值与文献[2]、[22]的实验值对比如图8所示。

    图  8  式(15)与文献试验值对比
    Fig.  8  Comparison between results of Eq.(15) and previous studies
    下载: 全尺寸图片

    图8可知,方程式(15)计算值与文献[2]和[22]的试验值吻合良好,说明本文提出由幂函数和尾迹函数组成的混合函数可以用于准确描述跌坎处压力竖向分布规律。图8中计算值与试验值出现一定量偏差,是由于文献[2]的试验水槽长度较短、或者较低流量时水位测针对水流的干扰,或者陡坡渠道小流量时水面的剧烈波动而造成。

    试验测得跌坎断面e−e最大压力值pmax,其与跌坎水深he之间对应关系如图9所示,同时将文献[2]的试验值也一并列入。

    图  9  pmax(ρg)–1和跌坎水深he间关系
    Fig.  9  Relationship between pmax(ρg)–1 and end depth he
    下载: 全尺寸图片

    图9可知,跌坎断面最大压力值pmax与跌坎水深he间存在良好的线性关系,经线性回归,可得:

    $$ \frac{{{p_{\max }}}}{{\rho g}} = 0.238{h_{\text{e}}} $$ (17)

    式(17)的相关系数为R2 = 0.937,说明pmaxhe存在良好的线性关系。同时,由图9可知,由式(17)所得的最大压力计算值与文献[2]的试验值也吻合一致。

    将式(17)代入式(16),则可得跌坎断面压力py的计算公式:

    $$ {p_{\text{y}}} = 0.575\rho g{h_{\text{e}}}{\left(\frac{y}{{{h_{\text{e}}}}}\right)^{0.36}}{\left( {1 - \frac{y}{{{h_{\text{e}}}}}} \right)^{1.36}} $$ (18)

    跌坎断面e−e的总压力为:

    $$ {P_{\text{e}}} = b\int_0^{{h_{\text{e}}}} {{p_{{y}}}} {\mathrm{d}}y $$ (19)

    则:

    $$ {P_{\text{e}}} = b\int_0^{{h_{\text{e}}}} {0.575\rho g{h_{\text{e}}}{{\left(\frac{y}{{{h_{\text{e}}}}}\right)}^{0.36}}{{\left( {1 - \frac{y}{{{h_{\text{e}}}}}} \right)}^{1.36}}{\text{d}}y} $$ (20)

    对式(20)进行积分变换并求解得:

    $$ {P_{\text{e}}} = 0.145b\rho gh_{\text{e}}^2 $$ (21)

    将式(2)(4)~(6)和(21)代入式(1),并积分整理得:

    $$ \begin{split} & \frac{{{Q^2}}}{{gb}} \left[ {b\left( \frac{1}{{{h_{\text{e}}}}} - \frac{1}{{{h_0}}} \right) + \frac{{S{{\left(\dfrac{{\sqrt S }}{n}\right)}^{ - 0.80}}{{\left(\dfrac{1}{{{h_0}}} + \dfrac{1}{{{h_{\text{e}}}}}\right)}^2}}}{{10.63gb{{\left[ 2.344 - 2\lg \left( \dfrac{{2{n^6}}}{{{{0.038\;2}^6}({h_0} + {h_{\text{e}}})}} \right) \right]}^2}}}} \right] =\\ &\qquad \frac{{0.50h_0^2 - 0.145h_{\text{e}}^2}}{{\sqrt {1 + {S^2}} }} + \frac{{1.5{{\left(\dfrac{{\sqrt S }}{n}\right)}^{ - 0.80}}S({h_0} + {h_{\text{e}}})}}{{\sqrt {1 + {S^2}} }} \end{split} $$ (22)

    式(22)为Qh0he的隐函数,根据测量的he值并结合式(8)可迭代求取h0,进而利用式(8)计算不同底坡S和糙率n的明渠流量。

    式(22)在推求过程中利用了渠道均匀流式(7),均匀流只有在渠道底坡为陡坡,即渠道底坡SSc时,跌坎上游才会出现均匀流流态。而对于渠道底坡S <Sc的缓坡或平坡甚至逆坡渠道,则跌坎上游不会出现均匀流流态,因而式(22)不能适用于渠道底坡为缓坡、平坡或逆坡渠道的流量计算。

    对于底坡S为缓坡、平坡或逆坡的矩形明渠,水流自跌坎自由跌落时,其流态存在缓流到急流的流态转变,因而跌坎上游出现临界流,逆坡渠道水流流态等同视为平坡。此时,O−O断面为临界流断面,断面水深为临界水深hc,则ho=hc。对于临界流,有:

    $$ Q = b{g^{0.5}}{h^{\frac{3}{2}}_{\text{c}}} $$ (23)

    此时,可以将自由跌流近似为堰高为0的薄壁堰流[9,13],因此可以应用推导薄壁堰流流量计算公式的方法推求其流量计算公式,则有:

    $$ {\text{d}}Q = b\sqrt {2g} \sqrt {H - y} {\text{d}}y $$ (24)

    式中,H为O−O断面的总水头,$H = {h_{\text{c}}} + \dfrac{{v_{\text{c}}^2}}{{2g}}$,m。在跌坎处,水流受重力影响而流速增加,水深减小,水深he=CchcCc表示反映跌坎处的水流垂直收缩系数。因此,总流量可表示为:

    $$ Q = b{C_{\text{c}}}\int_0^{{h_{\text{c}}}} {\sqrt {2g} \sqrt {H - y} } {\text{d}}y $$ (25)

    对式(25)积分并整理,得:

    $$ Q = b{C_{\text{c}}}\frac{{2\sqrt {2g} }}{3}\left( {{H^{\frac{3}{2}}} - {{\left( {H - {h_{\text{c}}}} \right)}^{\frac{3}{2}}}} \right) $$ (26)

    对临界流O−O断面,有:

    $$ F{r_{\text{c}}} = \frac{{{v_{\text{c}}}}}{{\sqrt {g{h_{\text{c}}}} }} = 1 $$ (27)

    则:

    $$ H = {h_{\text{c}}} + \frac{{g{h_{\text{c}}}}}{{2g}} = 1.5{h_{\text{c}}} $$ (28)

    故由式(26)和(28)可得:

    $$ Q = 1.484b\frac{{2\sqrt {2g} }}{3}\frac{{{h_{\text{e}}}}}{{{h_{\text{c}}}}}h_{\text{c}}^{\frac{3}{2}} $$ (29)

    由式(23)和式(29)可得:

    $$ {C_{\text{c}}} = \frac{{{h_{\text{e}}}}}{{{h_{\text{c}}}}} = 0.715 $$ (30)

    式(30)所得he/hc值与文献[4]的试验值(0.715)、刘清朝[27]、杨永森等[28]的数值计算值(0.720,0.716)基本一致。由式(30)可知,对于缓坡和平坡明渠,自由跌流处水深he与上游临界流断面O−O水深hc之比值he/hc,即水流垂直收缩系数Cc值为常数。因此,由式(30)和式(23)可得:

    $$ Q = 5.179bh_{\text{e}}^{\frac{3}{2}} $$ (31)

    文献[9]假定跌坎处0压力分布,运用能量方程得到了矩形明渠流量计算公式:

    $$ Q = F{r_{\text{o}}}{\left( {\frac{{{{(2 + Fr_{\text{o}}^2)}^{\frac{3}{2}}} - Fr_{\text{o}}^3}}{{3F{r_{\text{o}}}}}} \right)^{\frac{3}{2}}}b{g^{\frac{1}{2}}}h_{\text{e}}^{\frac{3}{2}} $$ (32)

    式中,$ F{r_{\text{o}}} $为O−O断面(图1)水流弗劳德数。

    对于平坡或缓坡渠道,O−O断面$ F{r_{\text{o}}} = F{r_{\text{c}}} = 1 $,则由式(32)可得:

    $$ Q = 5.181bh_{\text{e}}^{\frac{3}{2}} $$ (33)

    与本文式(31)基本一致。

    Dey [29]假定渠底流线曲率为0,利用Boussinesq近似来寻求跌坎处的压力分布规律,进而运用动量方程,得到粗糙矩形明渠流量计算公式为:

    $$ Q = \hat Qb{g^{0.5}}k_{\text{s}}^{{\text{1}}{\text{.5}}} $$ (34)

    式中,$ \hat Q = 1/\tilde k_{\text{s}}^{1.5} $,其中,$ {\tilde k_{\text{s}}} = {\hat k_{{\text{se}}}}{\tilde h_{\text{e}}} $$ {\hat k_{{\text{se}}}} = {k_{\text{s}}}/{h_{\text{e}}} $$ {\tilde h_{\text{e}}} = {\hat h_{\text{e}}}/ {\hat S^{ \frac{3}{10}}} $$ {\hat h_{\text{e}}} = {h_{\text{e}}}/{h_{\text{o}}} $$ \hat S = S/{S_{\text{c}}} $

    Tiğrek等[30]在大量试验的基础上,采用量纲分析的方法提出如下流量计算经验公式:

    $$\left\{\begin{aligned} &Q = 5.549bh_{\text{e}}^{\frac{3}{2}}\text{,}S<S_{\mathrm{c}};\\ &Q = \frac{{b{g^{0.5}}{{({h_{\text{e}}})}^{1.5}}}}{{{{\left(0.773 - 0.018\dfrac{{\sqrt S }}{n}\right)}^{1.5}}}}\text{,}S>S_{{\mathrm{c}}} \end{aligned}\right.$$ (35)

    从上述明渠跌水流量计算公式推求可以看出,对于陡坡渠道(S>Sc),流量除了和宽度b和跌坎水深he有关外,还与渠道底坡S和边壁糙率系数n均有关。而对于缓坡或平坡及逆坡渠道(S<Sc),流量仅和宽度b和跌坎水深he有关,而与渠道底坡S和边壁糙率系数n无关。

    根据本文试验获得的he,由式(8)和式(22)及式(31)(32)(34)(35)计算得到的流量计算值Qc和实测值QT之间的比较,结果如图10所示。

    图  10  不同公式计算得到的流量Q计算值与实测值比较
    Fig.  10  Comparison of calculated and tested flow by different equations proposed by researchers
    下载: 全尺寸图片

    图10可知,基于压力分布和动量方程提出的流量公式计算所得的计算值与实测值基本一致,对缓平坡(S<Sc)渠道,相对误差在2%以内,而陡坡渠道(S>Sc),其相对误差91.5%的值控制也在5%之内,达到较高的计算精度。而式(32)(34)(35)则仅在缓平坡渠道中计算值较准确(5%内),而在陡坡渠道中,流量计算值与实测值偏差较大,最大计算误差甚至超过47.2%。

    由上述理论和试验数据对比分析可见,对于跌坎断面压力分布如何处理是推求矩形跌流流量公式的关键,其不同的压力分布假定对流量公式计算精度具有重要影响。本文提出的跌坎断面压力分布规律方程对渠道糙率n和底坡S具有适用性,即在自由跌流流态下,跌坎断面压力分布规律具有自相似性,可由幂函数和尾迹函数构成的混合函数表述,不受渠道糙率n和底坡S变化值的影响。并且,本文试验中,渠道边壁糙率n值在0.009 3~0.019 3之间,底坡S在–0.011 2 ~ 0.053 4之间,公式具有较宽的应用范围。而文献[29]假定跌坎断面压力分布Boussinesq近似、文献[9]假定跌坎断面0压力分布,都与实际压力分布偏差较大,因而造成计算误差也较大,甚至达到60%。文献[30]中,n = 0.009 1和n = 0.014 7,n值较少,故其经验公式计算值与本文试验值也偏差较大,误差在30%以内。

    为进一步证实本文明渠流量计算式(22)和式(31)的适用性,文献 [31]和[24]的试验数据用来比较其测量值与式(22)及(31)计算值之间的差异性,结果如图11所示。

    图  11  式(22)和(31)流量计算值文献实测值对比
    Fig.  11  Comparison of the calculated flow by Eqs.22 & 31 and tested from references
    下载: 全尺寸图片

    图11可知,本文公式计算值与文献[24]和[31]实测值吻合良好。本文公式的计算值与试验值87.5%的误差在5%之内,只有极少数值的计算误差超过10%,表明本文提出的流量计算公式具有较高的计算精度。

    本文通过理论分析、模型试验并结合数值计算的方法,对矩形明渠跌流的断面压力分布及流量−水深关系进行了深入研究,明渠流量范围为5~100 L/s, 底坡范围为–0.0112~0.0534,边壁糙率系数为0.0093~0.0193,包括逆坡,平坡、缓坡和陡坡的粗糙明渠。主要结论如下:

    1)跌坎断面受自由跌流影响,压力分布显著偏离静水压强分布规律,但具有高度的自相似性,而这种自相似性与明渠边壁糙率和底坡无关,并提出了跌坎压力分布方程。

    2)利用动量方程推求了矩形渠道自由跌流的流量计算公式。本文公式相较于文献 [929-30]流量计量公式具有更高的精度和更宽应用范围,并且与文献[2431]试验数据一致。

    3)本文公式不仅适用于逆坡,平坡、缓坡和陡坡明渠流量计量,而且还适用于光滑和粗糙明渠流量计量明渠自由跌流的流量计量,针对平坡或缓坡渠道(无论光滑壁面还是粗糙壁面),其流量计量更加准确(误差在2%以内),具有重要的工程应用价值。

    4)因在推求过程中考虑了断面压力、阻力和重力,式(22)为流量Q的隐函数,形式较为复杂,同时式(22)不适用于平坡或缓坡渠道,因而应与式(31)结合,才可解决底坡和糙率变化的明渠自由跌流进行流量计量的问题,这可能限制公式的应用,因此,有必要在大量试验和理论分析的基础上,寻求更具有适用性的统一简化形式。

    5)当跌坎下游水深增加到一定值,其对跌落水流形成壅托,造成临界水流消失,此时水流为非自由跌落水流,跌坎下游水深会影响跌坎水深和流速,进而影响流量计量,此种情况下的明渠跌坎流量变化及其计算需进一步研究。

  • 图  1   矩形明渠自由跌坎过流示意图

    Fig.  1   Schematic view of a typical free overfall

    下载: 全尺寸图片

    图  2   试验装置

    Fig.  2   Experimental installation

    下载: 全尺寸图片

    图  3   不同断面压力分布(S = –0.011 2)

    Fig.  3   Pressure distribution at section of S = –0.011 2

    下载: 全尺寸图片

    图  4   不同断面压力分布 (S = 0)

    Fig.  4   Pressure distribution at section S = 0

    下载: 全尺寸图片

    图  5   不同断面压力分布 (S = 0.011 2)

    Fig.  5   Pressure distribution at section S = 0.011 2

    下载: 全尺寸图片

    图  6   不同断面压力分布 (S = 0.024)

    Fig.  6   Pressure distribution at section S = 0.024

    下载: 全尺寸图片

    图  7   $ \hat {\boldsymbol{p}} $$ {\boldsymbol{\theta}} $断面分布

    Fig.  7   $ \hat {\boldsymbol{p}} $${\boldsymbol{ \theta}} $ distribution on the end section

    下载: 全尺寸图片

    图  8   式(15)与文献试验值对比

    Fig.  8   Comparison between results of Eq.(15) and previous studies

    下载: 全尺寸图片

    图  9   pmax(ρg)–1和跌坎水深he间关系

    Fig.  9   Relationship between pmax(ρg)–1 and end depth he

    下载: 全尺寸图片

    图  10   不同公式计算得到的流量Q计算值与实测值比较

    Fig.  10   Comparison of calculated and tested flow by different equations proposed by researchers

    下载: 全尺寸图片

    图  11   式(22)和(31)流量计算值文献实测值对比

    Fig.  11   Comparison of the calculated flow by Eqs.22 & 31 and tested from references

    下载: 全尺寸图片
  • [1] 杨春蕾,蔡守华,王滇红,等.灌区量水技术发展历程及研究进展[J].灌溉排水学报,2017,36(增刊2):34-38.

    Yang Chunlei,Cai Shouhua,Wang Dianhong,et al.Review on research of discharge measuring technique in irrigation district[J].Journal of Irrigation and Drainage,2017,36(Supp2):34-38
    [2] Rajaratnam N,MuralidharD.Characteristics of the rectangular free overfall[J].Journal of Hydraulic Research,1968,6(3):233–258. doi: 10.1080/00221686809500236
    [3] Rajaratnam N,MuralidharD.The trapezoidal free overfall[J].Journal of Hydraulic Research,1970,8(4):419–447. doi: 10.1080/00221687009500324
    [4] Rouse H.Discharge characteristics of the free overfall[J].Civil and Environmental Engineering,1936,6(4):257–260.
    [5] Dey S.Free overall in open channels:State-of-the-art review[J].Flow Measurement and Instrumentation,2002,13(5/6):247–264. doi: 10.1016/S0955-5986(02)00055-9
    [6] Anderson M V.Non-uniform flow in front of a free overfall[J].Acta Polytechnic Scandinavia,1967,42:1–24.
    [7] Hager W H.Hydraulics of plane free overfall[J].Journal of Hydraulic Engineering,1983,109(12):1683–1697. doi: 10.1061/(asce)0733-9429(1983)109:12(1683)
    [8] Anastasiadou-Partheniou L,Hatzigiannakis E.General end-depth–discharge relationship at free overfall in trapezoidal channel[J].Journal of Irrigation and Drainage Engineering,1995,121(2):143–151. doi: 10.1061/(asce)0733-9437(1995)121:2(143)
    [9] Ferro V.Theoretical end-depth-discharge relationship for free overfall[J].Journal of Irrigation and Drainage Engineering,1999,125(1):40–44. doi: 10.1061/(asce)0733-9437(1999)125:1(40)
    [10] Dey S,Kumar B R.Hydraulics of free overfall in Δ-shaped channels[J].Sadhana,2002,27(3):353–363. doi: 10.1007/BF02703656
    [11] Ahmad Z.Quasi-theoretical end-depth-discharge relationship for rectangular channels[J].Journal of Irrigation and Drainage Engineering,2003,129(2):138–141. doi: 10.1061/(asce)0733-9437(2003)129:2(138)
    [12] Nabavi S V,Beirami M K,Chamani M R,et al.Free overfalls in flat-based circular and U-shaped channels[J].Flow Measurement and Instrumentation,2011,22(1):17–24. doi: 10.1016/j.flowmeasinst.2010.11.001
    [13] Ahmad Z,Azamathulla H M.Quasi-theoretical end-depth–discharge relationship for trapezoidal channels[J].Journal of Hydrology,2012,456/457:151–155. doi: 10.1016/j.jhydrol.2012.06.025
    [14] Rashwan I M H,Idress M I.Brink as a device for measurement discharge for partially filled circular channel[J].Ain Shams Engineering Journal,2013,4(2):133–141. doi: 10.1016/j.asej.2012.08.003
    [15] Vatankhah A R.Direct solution for discharge in generalized trapezoidal free overfall[J].Flow Measurement and Instrumentation,2013,29:61–64. doi: 10.1016/j.flowmeasinst.2012.09.006
    [16] Vatankhah A R.Power-law free overfall in subcritical flow regime[J].Ain Shams Engineering Journal,2015,6(2):399–402. doi: 10.1016/j.asej.2014.10.013
    [17] Noori B M A,Ibrahim S S.Effect of bed and side slopes on flow measurements in trapezoidal free overfall channels[J].Arabian Journal for Science and Engineering,2016,41(10):4187–4194. doi: 10.1007/s13369-016-2129-z
    [18] Dai Shubing,Jin Sheng.Hydraulics of free overfall in steeply sloping rough rectangular channel:A general computational approach[J].Flow Measurement and Instrumentation,2019,69:101625. doi: 10.1016/j.flowmeasinst.2019.101625
    [19] 李乃稳,何孟,刘超,等.陡坡粗糙矩形明渠中自由跌流的水力特性[J].工程科学与技术,2019,51(6):83–89. doi: 10.15961/j.jsuese.201900647

    Li Naiwen,He Meng,Liu Chao,et al.Hydraulic characteristics of free overfall in a steep and rough rectangular channel[J].Advanced Engineering Science,2019,51(6):83–89 doi: 10.15961/j.jsuese.201900647
    [20] Kraijenhoff D A,Dommerholt A.Brink depth method in rectangular channel[J].Journal of the Irrigation and Drainage Division,1977,103(2):171–177. doi: 10.1061/jrcea4.0001140
    [21] Guo Yakun,Zhang Lixiang,Shen Yongming,et al.Modeling study of free overfall in a rectangular channel with strip roughness[J].Journal of Hydraulic Engineering,2008,134(5):664–667. doi: 10.1061/(asce)0733-9429(2008)134:5(664)
    [22] Ramamurthy A S,Junying Qu,Diep Vo.Volume of fluid model for an open channel flow problem[J].Canadian Journal of Civil Engineering,2005,32(5):996–1001. doi: 10.1139/l05-038
    [23] Ramamurthy A S,Junying Qu,Diep Vo.VOF model for simulation of a free overfall in trapezoidal channels[J].Journal of Irrigation and Drainage Engineering,2006,132(4):425–428. doi: 10.1061/(asce)0733-9437(2006)132:4(425)
    [24] Firat C E.Effect of roughness on flow measurements in sloping rectangular free overfall [D].Ankara:Middle East Technical University,2004.
    [25] Yakhot V,Orszag S A,Thangam S,et al.Development of turbulence models for shear flows by a double expansion technique[J].Physics of Fluids A:Fluid Dynamics,1992,4(7):1510–1520. doi: 10.1063/1.858424
    [26] Ferziger J H,Peric M.Computational method for CFD[M].3rd ed.New York:Springer,2002.
    [27] 刘清朝.自由跌流的水力特性[J].水利学报,1990,21(11):8–17. doi: 10.13243/j.cnki.slxb.1990.11.002

    Liu Qingchao.Hydraulic characteristics of free overfall[J].Journal of Hydraulic Engineering,1990,21(11):8–17 doi: 10.13243/j.cnki.slxb.1990.11.002
    [28] 杨永森,杨永全.自由跌流的水力特性[J].水科学进展,2000,11(2):148–152. doi: 10.3321/j.issn:1001-6791.2000.02.007

    Yang Yongsen,YangYongquan.Hydraulic character of free overfall[J].Advances in Water Science,2000,11(2):148–152 doi: 10.3321/j.issn:1001-6791.2000.02.007
    [29] Dey S.End depth in steeply sloping rough rectangular channels[J].Sadhana,2000,25(1):1–10. doi: 10.1007/BF02703802
    [30] Tiğrek Ş,Firat C E,Ger A M.Use of brink depth in discharge measurement[J].Journal of Irrigation and Drainage Engineering,2008,134(1):89–95. doi: 10.1061/(asce)0733-9437(2008)134:1(89)
    [31] Davis A C,Ellett B G S,Jacob R P.Flow measurement in sloping channels with rectangular free overfall[J].Journal of Hydraulic Engineering,1998,124(7):760–763. doi: 10.1061/(asce)0733-9429(1998)124:7(760)
图(11)

本文结构

    /

    返回文章
    返回