渐扩式消力池体型优化研究

何志亚; 向鹏鹏; 洪彰华; 刘文超; Zhi-ya HE; Peng-peng XIANG; Zhang-hua HONG; Wen-chao LIU

中国农村水利水电 ›› 2022 ›› (5) : 190-194.

水电建设

渐扩式消力池体型优化研究

作者信息 +

Research on the Optimization of Gradual Expanding Stilling Basin

Author information +

稿件信息 +

摘要

在俄垤水库的水工模型试验中发现其渐扩式消力池内发生远驱式水跃，流态较差，无法满足消能防冲要求。为解决此问题，采用数值分析软件FLOW-3D探究了渐扩式消力池体型参数变化对池内流态、淹没度、消能率等的影响变化规律，结果表明：对于发生远驱式水跃的渐扩式消力池，增加尾坎高度可以显著提高淹没度，选择合适的尾坎高度不仅能改善池内流态还可将水跃消能率控制在合理范围内；增大扩散角和池长，对淹没度和消能率影响较小且于改善流态无益。根据结果推荐了合理的体型方案并通过物理模型验证。

Abstract

In the hydraulic model test of Edie reservoir， it was found that the repelled downstream hydraulic jump occurred in the gradual expanding stilling basin， and the flow state was poor， which could not meet the requirements of energy dissipation and scouring prevention. In order to solve this problem， numerical analysis software FLOW-3D was used to explore the influence of the body type parameters of the gradual expanding stilling basin on the flow state， inundation degree， energy dissipation rate and so on in the stilling basin， The results show that： for the gradual expanding stilling basin with the repelled downstream hydraulic jump， increasing the height of the tail can significantly improve the inundation degree， the appropriate height of the tail sill can not only improve the flow state in the stilling basin， but also control the energy dissipation rate to a reasonable range； the increase of diffusion angle and pool length has little effect on the inundation degree and energy dissipation rate and is not beneficial to improving the flow state. Based on the results， a reasonable shape parameters was recommended and verified by physical model. The research results can provide a reference for design optimization in the same type of engineering.

关键词

消力池 / FLOW-3D / 数值模拟 / 模型试验

Key words

stilling basin / FLOW-3D / numerical simulation / model test

基金

国家自然科学基金项目(50809052)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

何志亚 , 向鹏鹏 , 洪彰华 , 刘文超. 渐扩式消力池体型优化研究[J].中国农村水利水电, 2022(5): 190-194

Zhi-ya HE , Peng-peng XIANG , Zhang-hua HONG , Wen-chao LIU. Research on the Optimization of Gradual Expanding Stilling Basin[J].China Rural Water and Hydropower, 2022(5): 190-194

0 引言

对于水工建筑物，水跃消能是最早也是最完善的消能方式，它通过急流向缓流过渡时发生的水跃产生的表面漩滚和强烈的紊动消除余能，具有消能效果好，无雾化等优点，适用于中低水头，在国内外的水利工程得到了广泛运用。

随着消能技术的蓬勃发展，消力池的形式也日益多样化。对于低弗劳德数水跃，一般采用多级消力池^［1］。谭高文^［2］等通过数值分析结合物理模型试验对二级消力池体型参数及其敏感性进行了分析。向家坝在设计中采用了跌坎消力池，实现了底流消能在高水头、大单宽流量的运用^［3］。王海军^［4］等采用水力学试验方法也证明跌坎型消力池能有效解决高水头、大单宽流量消力池水力学指标过高的问题。我国林禀南院士和龚振瀛^［5］提出了堰顶收缩新技术，并在此基础发展出宽尾墩和消力池联合消能技术。潘艳华^［6］等通过模型试验给出了宽尾墩体型参数的选择原则。王均星^［7］等给出了Fr<2.0时的综合式消力池的体型设计形式并提出了一种新型消能尾坎。对于一些特定条件的水工建筑物，水流需由较窄的泄槽过渡到较宽阔的河槽，通常采用扩散消力池进行衔接，但由于边墙的扩散，水跃末端尾水宽度加大，降低了跃后水深，无法形成典型水跃^［8］。扩散水跃属于空间三维问题，各研究单位所给扩散水跃计算公式应根据实际情况适用。张志昌^［9］等给出了渐扩综合式消力池深度、坎高和尾坎作用力的计算公式，可供参考，而渐扩综合式消力池体型参数变化对水力特性具体的影响还未有研究。

依托云南省勐龙河上游的俄垤水库溢洪道改建工程，主要研究了渐扩式消力池特征体型参数变化对消能效果的影响，分析了各参数组合方案下渐扩式消力池消能时水力特性参数（池内流态、沿程流速、消能率等）的变化规律，从而选择了合适的消力池体型来满足消能防冲要求和节省工程投资。

1 工程概况

俄垤水库改建后的溢洪道为岸边侧槽开敞式溢洪道，溢洪道轴线长375.98 m，宽4.00 m，堰顶高程与水库正常蓄水位 1 570.25 m齐平，尾水渠高程1 501.91 m，总落差达68.34 m。消能防冲水位1 572.04 m，设计水位1 572.32 m。

初设方案采用渐扩式矩形断面消力池进行底流消能，底宽由4 m扩至8 m。池深3.0 m，池长38.0 m。边墙高度6.9 m。

初设方案下消力池布置如图1所示。

图1 初设方案消力池布置示意图（单位：m）

Fig.1 Profile and layout of the initial design of stilling basin

Full size|PPT slide

2 初设方案模型试验成果

物理模型试验按照重力相似准则^［10］，选用几何比尺1∶30制作整体模型，相应的流速比尺

λ_{v} = λ_{L}^{1 / 2}

=5.48，流量比尺

λ_{Q} = λ_{L}^{5 / 2}

=4 929.50，糙率比尺

λ_{n} = λ_{L}^{1 / 6}

=1.76。

消能防冲水位下，水流在消力池内发生了远驱式水跃，整体流态较差。随着上游库水位上升至设计水位，下泄流量进一步增加，池内主流直冲消力池尾坎，高高跃起，溅起水花远超边墙高度。初设方案下的消力池体型需考虑优化，优化后的理想流态应该为池末水深适当壅高，使池内发生稍有淹没的水跃。

不同水位下池内流态如图2所示。

图2 初设方案不同水位下消力池流态

Fig.2 Flow pattern of the initial design stilling basin at different water levels

Full size|PPT slide

3 数值模拟优化

3.1 数值模型

为探究渐扩式综合消力池体型特征参数变化对水力特性的影响，选取渐扩式消力池扩散角、池长、尾坎高度的不同参数组合进行数值模拟研究。不同组合下的消力池体型参数如表1所示，其中组次1为初设方案。

表1 不同组合消力池体型参数

Tab.1 Body parameters of different combinations of the stilling basin

组次	扩散角/（°）	池长/m	尾坎高/m	影响因素
1	3.0	38	0
2	3.5	38	0	扩散角
3	4.0	38	0
4	3.0	39	0	池长
5	3.0	40	0
6	3.0	41	0
7	3.0	38	0.5	尾坎高度
8	3.0	38	1.0
9	3.0	38	1.5

数值模型Flow-3d软件中采用RNG k-ε紊流模型，自由液面追踪采用TruVOF法，采用有限体积法进行离散。由于溢洪道的纵横比过大，模拟效果不佳且计算时间长，故仅模拟消力池（包括进口上游25.0 m）及尾水渠段。根据上下游高差进行理论计算并结合物理模型实测数据对进入消力池的水流给定一个初始速度。本次模拟计算工况为设计水位工况，在消力池进口上游25.0 m处设置边界条件为体积流量，流量为67.1 m³/s；水流初始速度为23.3 m/s，顺泄槽向下游；下游设为自由出流边界，在尾水渠尾部设置流量监测点，用以校验模型流量。

3.2 网格划分及无关性验证

初设方案物理模型试验中得到溢洪道泄槽尾部水深在0.9 m左右，需较高的网格精度来模拟。模型整体共划分4个网格块，网格块1为消力池上游，网格块2为消力池部分，网格块3为消力池尾坎部分，网格块4为出水渠部分，如图3所示。

图3 网格划分示意图

Fig.3 Schematic diagram of mesh

Full size|PPT slide

针对初设方案体型（组次1），模型网格块进行解析时选择了4组网格精度进行模拟计算。根据计算结果，组合3采用的网格尺寸（网格块1为0.15 m、网格块2为0.22 m、网格块3、4为0.25 m）解析精度已能满足模拟要求。该组合下模拟结果与物理模型实测数据的对照如表2所示。

表2 数值模型与物理模型数据对照表

Tab.2 Comparison between numerical model and physical model data

参数	消力池进口		跃前断面			跃后断面			消力池出口
参数	水深/ m	流速/ （m·s^-1）	水深/ m	位置	流速/ （m·s^-1）	水深/ m	位置	流速/ （m·s^-1）	水深/ m	流速/ （m·s^-1）
数值模型	0.84	20.89	1.27	9.20	20.51	6.20	32.50	2.73	5.92	3.02
物理模型	0.80	22.75	1.43	9.53	21.35	6.19	34.33	4.32	6.19	3.15

注：跃前、跃后断面位置以消力池进口为0 m，下同。

通过初设方案体型的计算结果与物理模型实测数据对比得出：后续数值模拟可按照该网格尺寸组合进行模拟分析。根据计算时的稳定性时间步长变化曲线及监测面流量变化可知池内流态在t=85 s左右趋于稳定，为取得可靠结果选取t=120~160 s的数据进行分析。

3.3 数值模型成果及分析

3.3.1 流态分析

各模拟组次流速水面线图如图4所示。

图4 各组次方案流速水面线

Fig.4 Velocity & water surface line of each scheme

Full size|PPT slide

由图4可知，各模拟方案中随着扩散角增大，尾坎过水断面变宽，水跃发生位置略微后移，尾坎处水深明显减小，且角度过大时水流入池后不能贴壁下泄；随着池长增加，尾坎处水深无明显变化但跃前断面后移，未见池内流态明显改善，且池长增加导致工程量增加；随着尾坎高度增加，跃前断面不断前移，跃前水深增大，池末水深也有明显增加。由此可见尾坎高度对渐扩式消力池流态改善显著。根据模拟的池内流态分析［图4（h）与图4（i）对比］，该工程特征流量下最佳尾坎高度应该在第8、9组之间，可以在池内形成稍有淹没的流态，需进一步通过数值模拟验证。

3.3.2 底部流速分布

消力池内关键部位底部流速及水力参数如表3所示。

表3 各组次方案流速及特征水力参数

Tab.3 Flow velocity and characteristic hydraulic parameters of each scheme

组次	入池流速/ （m·s^-1）	收缩断面			跃前断面			出池流速/ （m·s^-1）
组次	入池流速/ （m·s^-1）	流速/（m·s^-1）	水深/m	Froude数	流速/（m·s^-1）	水深/m	Froude数	出池流速/ （m·s^-1）
第1组	20.89	21.95	0.66	8.63	20.51	1.27	5.81	3.02
第2组	21.39	22.09	0.64	8.82	21.70	1.23	6.25	2.86
第3组	20.95	21.58	0.63	8.68	20.32	1.12	6.13	2.62
第4组	21.39	22.07	0.58	9.26	21.18	1.29	5.96	3.38
第5组	20.95	22.10	0.57	9.35	18.91	1.37	5.16	3.97
第6组	20.96	21.35	0.55	9.20	17.79	1.25	5.08	4.04
第7组	21.38	22.07	0.74	8.20	20.56	1.29	5.78	6.60
第8组	21.38	22.01	0.87	7.54	20.68	1.32	5.75	5.73
第9组	17.04	22.34	0.92	7.44	21.37	1.30	5.99	6.46

由特征水力参数可知：随着扩散角增大，尾坎过水断面变宽，出池流速减小，收缩断面Froude数位于略小于9.0，水跃对尾水深度敏感；而池长的增加使得收缩断面Froude数均大于9.0，池内已非稳定水跃；尾坎高度增加后，收缩断面Froude数逐渐减小，水跃逐渐稳定。

3.3.3 消能率分析

通过模拟结果所得的流速水深计算消力池的消能率，采用消力池进口上游10 m断面和消力池尾坎断面的能量差来计算，示意图见图5。

图5 消力池消能率计算断面示意图

Fig.4 Schematic diagram of computation section of energy dissipation ratio of the stilling basin

Full size|PPT slide

水跃消能率根据公式（1）计算：

\{\begin{matrix} E_{1} = h_{1} + \frac{a_{1} v_{1}^{2}}{2 g} \\ E_{2} = h_{2} + \frac{a_{2} v_{2}^{2}}{2 g} \\ η = \frac{E_{1} - E_{2}}{E_{1}} \times 100 % \end{matrix}

（1）

式中：

h_{1}

、

h_{2}

分别为断面1-1、断面2-2相对于消力池底部的位能；

a_{1}

、

a_{2}

为动能修正系数，取1.0；

v_{1}

、

v_{2}

分别为断面1-1、断面2-2断面平均流速；

g

为重力加速度，取9.8 m²/s。

水跃淹没度

σ_{j}

采用式（2）计算：

σ_{j} = h_{T} / h_{c 1}^{''}

（2）

式中：

h_{T}

为池末水深；

h_{c 1}^{''}

为跃后断面水深。

各组次方案下水跃消能率及淹没度模拟结果如表4所示。

表4 各组次方案消力池水跃消能率及淹没度

Tab.4 Energy dissipation rate and submergence degree of water jump of each schemes

组次	1	2	3	4	5	6	7	8	9
水跃消能率/%	70.72	71.45	72.26	71.27	71.48	72.35	71.15	68.48	65.50
淹没度 $σ_{j}$	0.83	0.84	0.84	0.81	0.81	0.80	0.87	1.08	1.27

在此次数值模拟条件下，扩散角与消力池长度的增加都能使消能率增大，尾坎高度的变化对消能率影响显著；扩散角的增大和池长的增加对淹没度影响不大，但尾坎高度的增加能明显增加淹没度。

3.3.4 数模推荐方案

由上述数值模拟结果得到，尾坎高度对消力池内流态和消能率的影响明显。为了在消力池内水力特征参数满足规范要求的条件下，尽可能让池内最大水深增幅不至于过大，选择保持消力池原设计池长和扩散角度不变，仅在消力池末端增加尾坎。尾坎高度1.2 m时池内流态达到最佳，可作为推荐方案（如图6所示）。

图6 推荐方案流速水面线

Fig.6 Velocity & water surface line of recommended scheme

Full size|PPT slide

推荐消力池体型参数如下：池长38.00 m，池宽由进口4.00 m直线渐变增加至出口8.00 m，池末尾坎高度1.2 m，坎顶宽1.0 m，下游面坡度1∶1，除增添尾坎外消力池原有的体型参数保持不变。考虑到数值模拟对流态细节模拟的局限性，边墙高度暂定8.0 m。推荐方案剖面图如图7所示，平面布置图与图1（b）一致。

图7 推荐方案纵剖面（单位：m）

Fig.7 Longitudinal profile of recommended scheme

Full size|PPT slide

4 推荐方案物理模型试验

推荐方案物理模型试验中，在消能防冲水位及设计水位下，消力池内均发生淹没水跃，池内流态如图8所示，水跃特征值参数如表5所示。

图8 推荐方案模试验型消力池流态

Fig.8 Flow pattern of the recommended scheme in stilling basin

Full size|PPT slide

表5 消力池水跃特征参数

Tab.5 Characteristic parameters of hydraulic jump in stilling basin

工况	参数
工况	收缩断面流速/（m·s^-1）	收缩断面Froude数	跃后水深实测值/m	池末水深h_T /m	淹没系数 $σ_{j}$	水跃消能率/%
消能防冲水位	23.10	8.56	6.05	6.69	1.11	68.4
设计水位	24.76	8.95	6.51	7.56	1.16	58.5

从图8可以看出，推荐方案下消力池内水跃跃前位置均位于消力池进口。相较于初设方案，跃前断面水深明显增加，尾坎处未出现水流撞击尾坎而溅起水花的现象，池内流态较好。

由表5可知在消能防冲水位和设计水位下，水跃收缩断面弗劳德数略小于9.0，池内水跃基本稳定；池末水深也较初设方案略有增加。池内水跃淹没度

σ_{j}

位于1.10左右，消力池内能形成稍有淹没的水跃。水跃消能率较高，消能充分。池末边墙高度相较初设方案增幅不大，推荐方案体型能满足工程要求。

5 结语

通过FLOW-3D结合俄垤水库扩建工程实例探究了渐扩式消力池体型参数对池内水力特性的影响，通过数值模拟发现：对于发生远驱式水跃的渐扩式消力池，增加尾坎高度可以显著提高淹没度，选择合适的尾坎高度不仅能改善池内流态还可将水跃消能率控制在合理范围内；增大扩散角和池长，对淹没度和消能率影响较小且于改善流态无益。

本文的推荐方案经济方便，且技术上可行，可作为类似工程的借鉴。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	吴子荣，包中进. 多级消力池水力特性的研究及其应用［J］. 水利水电技术， 1998（12）：17-19. 本文引用 [1]

2	谭高文，韩昌海，余凯文. 低弗劳德数水流二级消力池体型参数试验研究［J/OL］. 水科学进展， 2018，29（6）：848-857.DOI：10.14042/j.cnki.32.1309.2018.06.010 . 本文引用 [1]

3	孙双科，柳海涛，夏庆福，等. 跌坎型底流消力池的水力特性与优化研究［J］. 水利学报， 2005（10）：1 188-1 193. 本文引用 [1]

4	王海军，赵伟，杨红宣，等. 跌坎型底流消能工水力特性的试验研究［J］. 水利水电技术， 2007（10）：39-41. 本文引用 [1]

5	林秉南，龚振瀛，王连祥. 突泄坝址过程线简化分析［J/OL］. 清华大学学报（自然科学版）， 1980（1）：17-31.DOI：10.16511/j.cnki.qhdxxb.1980. 01. 003 . 本文引用 [1]

6	李中枢，潘艳华，韩连超，等. 宽尾墩联合消能工体型选择及水力特性的研究［J］. 水科学进展， 2000（1）：82-87.DOI：10.14042/j.cnki.32.1309.2000.01.015 . 本文引用 [1]

7	刘洪滨，王亚洲，王均星，等. 低Froude数新型消能尾坎数值模拟研究［J］. 中国农村水利水电， 2020（8）：221-226，230. 本文引用 [1]

8	吴宝琴，张志恒.矩形扩散水跃水力计算新公式［J］.水利水电工程设计，2001（2）：42-44. 本文引用 [1]

9	张志昌，傅铭焕，贾斌，等. 渐扩综合式消力池尾坎高度和作用力的计算［J］. 水资源与水工程学报， 2014，25（6）：137-141. 本文引用 [1]

10	水工（常规）模型试验规程：SL 155-2012 ［S］. 水利部水利水电规划设计总院，2012. 本文引用 [1]

PDF(1505 KB)

访问

引用

详细情况

段落导航

摘要
Abstract
关键词
Key words
基金
引用本文
0 引言
1 工程概况
图1 初设方案消力池布置示意图（单位：m）
2 初设方案模型试验成果
图2 初设方案不同水位下消力池流态
3 数值模拟优化
3.1 数值模型
表1 不同组合消力池体型参数
3.2 网格划分及无关性验证
图3 网格划分示意图
表2 数值模型与物理模型数据对照表
3.3 数值模型成果及分析
3.3.1 流态分析
图4 各组次方案流速水面线
3.3.2 底部流速分布
表3 各组次方案流速及特征水力参数
3.3.3 消能率分析
图5 消力池消能率计算断面示意图
表4 各组次方案消力池水跃消能率及淹没度
3.3.4 数模推荐方案
图6 推荐方案流速水面线
图7 推荐方案纵剖面（单位：m）
4 推荐方案物理模型试验
图8 推荐方案模试验型消力池流态
表5 消力池水跃特征参数
5 结语
参考文献

收稿日期	出版日期
2022-03-14	2022-05-15
发布日期
2022-06-28

选择文件类型/文献管理软件名称

选择包含的内容