黄智敏, 付 波, 陈卓英
(广东省水利水电科学研究院 广东省水动力学应用研究重点实验室, 广东 广州 510635)
低水头水闸运行特点为泄洪流量大、水位差较小,其泄洪消能多采用底流消能方式。在给定的水闸泄流水力参数条件下,根据工程具体条件,可先计算和确定水闸下游消力池的池末端尾坎顶高程、池深等,然后再计算和确定消力池长度。水闸下游消力池末端尾坎顶高程、池深等计算和确定方法可参考文献[1-4] ,但目前国内有关水闸下游消力池合理池长的研究成果仍较少[5-10]。
水闸下游消力池池长通常可按文献[1] 建议的公式计算,但由于水闸下游消力池的水流条件较复杂,水闸出流条件(自由出流、淹没出流)对闸下游流态和水跃长度影响较大,若不正确区分水闸下游的流态,则计算得出的消力池池长明显不合理,会增加工程投资和造成浪费。
本文结合广东省鉴江高岭拦河水闸重建工程下游消力池水力模型试验研究成果[11],在水闸闸门全开泄流(即闸门全开敞泄,下同)的下游出流流态判别的基础上,提出了低水头水闸下游消力池合理长度的计算和选取方法。
常规水闸下游消力池长度计算公式为公式(1)~(4)[1,12],相应的消力池体型和水力参数见图1。
图1 水闸下游消力池体型和水力参数示意图
消力池长度Lsj
Lsj=Ls+βLj
(1)
水跃长度Lj
Lj=6.9(h2-h1)
(2)
其中:
(3)
(4)
式中:Ls为消力池斜坡段水平投影长度,m;β为水跃长度校正系数,可取0.7~0.8;h1为消力池进口收缩断面水深,m;q为泄流单宽流量,m3/(s·m);E0为以消力池底板为基准计算的总水头,m;φ为流速系数;g为重力加速度,m/s2;h2为跃后水深,m;α为水流动量校正系数,可采用1.0~1.05;b1为消力池首端宽度,m;b2为消力池末端宽度,m。
由公式(1)~(4)分析,在已知水闸泄流单宽流量q、闸上游水位Zu、闸下游水位Z与流量Q关系等水力参数条件下,就可以计算出消力池底板总水头E0、消力池的水跃参数(如h1、h2及Lj) 等,从而可进一步计算出消力池长度Lsj。
由公式(2)分析,正常水跃的跃长Lj随消力池跃后水深h2与跃前水深h1的差值(h2-h1)增大而增加,而h1、h2和(h2-h1) 等随泄流单宽流量q增加而增大,因此,消力池水平段长度βLj随q增加而增大,其通常可由水闸最大泄洪流量Qm来确定。对于平原地区河流或闸上、下游水位差较小的水闸,在水闸闸门全开泄流、最大泄洪流量Qm运行时,其上、下游水位差较小(计算闸孔总净宽时,平原区水闸的过闸水位差可采用0.1~0.3 m[1]),出闸水流多呈波状的明渠流流态,闸下游消力池无水跃产生。此时,若不正确判别水闸出流流态,笼统按Qm计算的消力池长度会明显偏大,从而造成工程量增加和投资浪费。
通常,水闸泄流调控运行有两种方式:一是在水闸上游为正常蓄水位条件下,水闸闸门局部开启、控泄运行,此运行方式水闸上、下游水位差相对较大;二是当上游洪水来流量较大时,为了确保水闸工程的安全运行,水闸闸门需全部开启敞泄,此运行方式的闸上、下游水位差相对较小。这两种运行方式的水闸下游出流流态是有差异的。
图2是高岭拦河水闸重建工程水力模型试验研究的下游消力池推荐方案剖面图,拦河闸共设16孔闸,单孔闸净宽15 m,闸墩厚为2.5 m和3.0 m两种,闸下游消力池池底高程-1.00 m、水平段长度22.5 m。拦河闸下游消力池水力模型试验在60 cm宽的玻璃水槽中进行,水力模型按重力相似律设计为正态,几何比尺Lr=46.67。为了分析下游消力池长度的合理性,进行了水闸闸门两种运行开启工况的试验比较:(1)水闸闸门全开、泄放最大洪水流量(校核洪水流量,P=0.5%)运行;(2)水闸上游为正常蓄水位(8.20 m)、控泄运行的最大闸门开度em运行。
图2 高岭水闸下游消力池布置剖面图(单位:m)
在最大洪水流量(P=0.5%,闸孔泄流单宽流量q=23.99 m3/(s·m))、不同闸下游河道水位Z的水力模型试验水力参数见表1,闸下游运行的流态和流速分布见图3~5。
(1) 当闸下游河道水位Z=5.00~6.00 m时,出闸水流为自由出流,以消力池底为基准的下游水深(ht+d)
表1 高岭水闸及消力池运行水力参数 m
注:(1)ht为下游河床水深; (2)hs为以闸室堰顶为基准的下游河床水深(见图1)。
图3 水闸下游消力池流态和流速分布示意图(Z=5.00、6.00m)
(2)当闸下游河道水位上升为Z=7.00~8.00 m时,其下游水深(ht+d)>h2、且hs/H<0.8 (H为闸室堰顶的上游水深,hs为以闸室堰顶为基准的下游河床水深),此时,出闸水流仍为自由出流,消力池内逐渐形成稳定的水跃,池内水跃长度随闸下游相对水深hs/H增加而减小,即池内跃尾断面随hs/H增加往上游移动(见表1和图4)。
(3) 当Z=8.50 m(hs/H=0.8)时,出闸水流为自由出流和淹没出流的临界状态 (或称淹没出流[12]),池内跃尾断面位于约桩号0+053 m断面。随着闸下游水位Z继续增加,池内水跃长度减短,跃尾断面向消力池上游移动。当闸下游河道水位上升至Z=9.00~10.00 m (hs/H=0.88~0.98,淹没出流),消力池内无水跃、呈波状的明渠流流态;受水闸闸室平底堰的影响,下游消力池上游段底部出现局部回流区,主流在消力池水流的上部(见表1和图5)。
因此,高岭水闸在下泄校核洪水流量(P=0.5%,Z=12.16 m) 运行条件下,闸上、下游水位差△Z<0.1 m,出闸水流无水跃产生,闸下游呈波状的明渠流态[11](见图6)。
图4 水闸下游消力池流态和流速分布示意图(Z=7.00、8.00m)
图5 水闸下游消力池流态和流速分布示意图(Z=8.50、9.00、10.00m)
图6 水闸下游消力池流态和流速分布示意图(Z=12.16m)
在水闸上游正常蓄水位8.20 m、闸门最大开度em=0.75 m(Z=3.52 m)的泄流运行工况下,其下游消力池的流态和流速分布见图7。测试的消力池水跃跃尾断面桩号约为0+058 m,比水闸泄放校核洪水频率流量(P=0.5%)、闸下河道水位Z=8.50 m(即hs/H=0.8,自由出流和淹没出流的临界状态)相应的水跃长度要长。
图7 水闸下游消力池流态和流速分布示意图
采用相关的水闸下游消力池计算公式[1,12],对高岭水闸下游消力池长度进行计算,计算结果见表2,对计算结果分析如下:
(1) 在闸上游正常蓄水位8.20 m、闸门控泄最大开度em=0.75 m(q=4.93 m3/(s·m)、Z=3.52 m)条件下,计算的下游水跃长度Lj=20.63 m,其消力池水平段长度可取0.8Lj=16.5 m。
(2) 在校核洪水频率流量(P=0.5%、闸门全开泄洪、Zu=12.22 m、Z=12.16 m)条件下,计算的下游水跃长度Lj=43.06 m,其消力池水平段长度可取0.8Lj=34.45 m。
由表2可见,校核洪水频率流量(P=0.5%)泄流的消力池下游水深(ht+d)明显大于其计算的跃后水深h2,消力池内产生高淹没水跃;且由于闸下游相对水深hs/H>0.8,水闸为淹没出流。因此,可判断水闸下游无法产生正常的水跃、闸下游为波状的明渠流(见图6)。
(3) 经综合分析各种泄流工况的水闸下游流态和池长等,并取一定的安全系数,高岭水闸重建工程下游消力池水平段长度选取为22.5 m ,为闸上游正常蓄水位8.20 m、闸门控泄最大开度em=0.75 m计算的水跃长度Lj的1.09倍(见图2)。
表2 水闸下游消力池长度计算表
注:h1、Fr1为消力池水流收缩断面的水深和佛氏数。
在水闸闸门全开、泄洪流量一定的条件下,水闸下游流态和水跃长度与其下游河道水位Z有关。在水闸下游消能设计Z~Q关系确定的条件下,水闸下游消力池水平段长度βLj应考虑以下的因素进行选取:
(1)水闸闸门全开、宣泄最大洪水流量Qm运行时,若hs/H<0.8(水闸出流为自由出流),应采用Qm,按文献[1] 等计算和选取消力池长度。
(2)水闸闸门全开、宣泄最大洪水流量Qm运行时,若hs/H≥0.8(水闸出流为淹没出流),闸下游水跃长度较短或呈无水跃状态,此时,不应采用Qm计算和选取消力池长度,而应采用水闸上游为正常蓄水位、闸门控泄最大开度em的泄量计算和选取消力池长度。
在此运行条件计算消力池长度时,为了确保工程的安全运行,闸门控泄最大开度em对应的闸下游河道水位Z应选取低一级闸门开度泄流稳定之后的下游水位[13](如高岭水闸重建工程在正常蓄水位8.20 m条件下,闸门控泄的每一级闸门开度间隔设置为△e=0.25 m,其闸门控泄最大开度em=0.75 m运行对应的Z应选取闸门开度e=0.5 m泄流稳定之后的下游河道水位);同时,为了确保工程的安全运行,建议消力池水平段长度βLj的水跃长度校正系数β取1.0。
在实际计算中,为了方便正常蓄水位、闸门控泄条件下消力池长度的计算和选取,通常可根据闸门控泄的各级开度e分别计算消力池的长度(计算中,各级闸门开度e运行的Z为低一级闸门开度泄量对应的下游河道水位;同时,水跃长度校正系数β取1.0),选取计算的最大值为消力池长度。
(3) 按照文献[1],对于大型工程的下游消力池布置和体型尺寸,应通过水力模型试验论证之后确定。
本文根据水闸闸门全开、宣泄最大洪水流量等运行条件的水力模型试验研究,提出了低水头水闸下游消力池长度的计算方法:
(1) 水闸闸门全开、宣泄最大洪水流量Qm运行时,若hs/H<0.8(水闸出流为自由出流),应采用Qm,按水闸设计规范[1]等计算和选取消力池的长度。
(2) 水闸闸门全开、宣泄最大洪水流量Qm运行时,若hs/H≥0.8(水闸出流为淹没出流),闸下游水跃长度较短或呈无水跃状态,此时,应采用水闸上游正常蓄水位、闸门控泄最大开度em泄量的运行工况,计算和选取消力池长度。
本文研究成果已应用于工程实践,在确保工程安全运行的前提下,可大大节省工程投资,工程效益较为显著,成果可供工程设计和运行参考。