坎式阶梯溢流堰的冲刷特性研究.pdf

第10 期2023年10 月文章编号：16 7 3-9 0 0 0（2 0 2 3)10-0 0 0 7-0 4陕西水利Shaanxi WaterResourcesNo.10October,2023坎式阶梯溢流堰的冲刷特性研究赵海，高海静，裘丽恩，钟绍柱，毛欣？（1.杭州市富阳区林业水利局，浙江杭州3 1140 0；2.浙江水利水电学院，浙江杭州3 10 0 18）【摘要坎式阶梯溢流堰，即在传统阶梯溢流堰的台阶上增加尾坎结构，达到减小流速、减少冲刷的目的。为研究阶梯水平面上的尾坎位置与冲刷特性之间的关系，运用Flow-3D软件对坎式阶梯溢流堰的冲刷特性进行数值模拟，得到尾坎与台阶水平面末端出口之间的距离为0.5 m和10 m条件下，坎式阶梯溢流堰的水流流态、最大冲坑深度和位置、冲坑长度等。结果表明尾坎布置越靠近阶梯溢流堰的台阶末端，消能效果越好，冲刷深度更小。【关键词坎式阶梯；溢流堰；冲刷特性；数值模拟【中图分类号TV132+.2【文献标识码BStudy on Scour Characteristics of Sill-Type Stepped Overflow WeirZhao Hai,Gao Haijing,Qiu Lien,Zhong Shaozhu,Mao Xin?(1.Forestry and Water Resources Bureau of Fuyang District,Hangzhou 311400,Zhejiang;2.Zhejiang University of WaterResources and Electric Power,Hangzhou 310018,Zhejiang)Abstract:The sill-type stepped overflow weir is to add the tail sill structure on the steps of the traditional steppedoverflow weir,so as to reduce the flow velocity and reduce the erosion.In order to study the relationship between theposition of the tail sill on the horizontal plane of the staircase and the erosion characteristics,Flow-3D software was used tonumerically simulate the erosion characteristics of the cascade overflow weir,and the flow pattern,maximum crater depth andlocation,and crater length of the cascade overflow weir were obtained under the condition that the distance d between the tailsill and the end outlet of the horizontal plane of the staircase was 0.5 m and 10 m.The results show that the closer the tail sill isto the end of the step of the stepped overflow weir,the better the energy dissipation effect is and the less the erosion depth is.Key words:Sill-type step;overflow weir;scour characteristics;numerical simulation蚀破坏的直接原因。对不同泄槽角度、消能池底高度、尾1引言水深度和流速条件下的冲刷过程进行了研究。确定了冲刷深每个时期的水利建设都是特定自然和社会背景的产物。近年来，伴随快速发展的水利行业，阶梯式溢流堰在各种水利工程中得到了普遍的运用。坎式阶梯式溢流堰是近几年成长起来的新颖的消能工型式，这种新型阶梯式溢流堰是在传统阶梯式溢流堰的末端设置一个尾坎而形成的。早在5 0 0 0多年前古埃及用溢流堰的结构设计过早期坝，随着科技的发展，逐渐有学者开始采用数值模拟的方法研究台阶溢流堰上复杂的水流水力特性。阶梯式溢流堰可以采用粗糙坡度台阶和设置尾坎来消耗部分能量。然而，在大流量、高水头、高流速条件下，阶梯式溢流堰的外部凸起形状是造成建筑物空度依赖于台阶几何，随台阶高度的增加而缩小，而随流量以及溢流堰泄槽角度的增加而放大。彭勇等3-4 研究发现可以通过设置前置掺气坎提高阶梯溢流堰消能率，从而解决了大单宽流量下阶梯溢流堰运行难题。后又有王业红、肖兴斌结合模型试验研究表明5，泄洪道阶梯式后，结合加气减蚀的办法，在梯级平台坝面加气可减少冲蚀。吴守荣等学者对平滑溢流堰、阶梯式溢流堰和前加气坎式阶梯式溢流堰的水力特性对比着来研究。认为采纳前加气，坎式阶梯式溢流堰不仅可以保证阶梯式溢流堰的高能量耗散比率，还可以维护保证台阶表面免受空化的破坏，解决了溢流堰形状优化及下收稿日期 2 0 2 2-11-11【基金项目】2 0 2 1年浙江省重点水利科技项目（RB2119）；2 0 2 1年浙江省一般水利科技项目（RC2143）；国家级大学生创新创业训练计划项目（2 0 2 2 1148 10 3 4）【作者简介赵海（197 8-），男，浙江杭州人，高级工程师，主要从事工程规划、技术审查、工程管理等工作。【通讯作者毛欣（2 0 0 1），男，浙江义乌人，本科在读，主攻方向：水力学。7第10 期2023年10 月游抗冲刷问题。尽管前人已有大量文献研究了水流单宽流量、阶梯高度、尾水水深以及坡度等因素对阶梯溢流堰冲刷特性的影响8,阶梯溢流堰被普遍应用于各种水利工程实践中，但目前对坎式阶梯溢流堰冲刷特性的相干论述非常少。基于此，本文通过CAD软件建立尾坎与台阶9 水平面出口之间的距离d分别为0m、5 m、10 m 的坎式阶梯溢流坝模型导人Flow-3D软件进行数值模拟计算。对结果进行分析，研究阶梯水平面上的尾坎位置对冲刷特性的影响，建立它们之间的关系，以便相关结论应用到坎阶梯溢流堰的建造防护中，拓展水利工程界对坎式阶梯溢流堰消能的认识。2试验装置和方法2.1几何模型本文是利用Flow-3D制作几何模型。选用AutoCAD进行建模，利用Flow-3D软件来模拟流动特性。实验模型由上游库区、宽顶堰、带尾坎阶梯段、出水渠组成0,底坡26.6。其中，宽顶堰宽W=52m，长Lerest=101m，弯角是一个r=6m的倒角；设计的是10 阶模型坎式阶梯溢流堰阶梯，每个台阶尺寸相同：长l=20m,高h=10m；台阶水平面上尾坎高度w=5m，厚度lw=2m，尾坎距台阶出口末端距离为d。实验中尾坎位置有三种,距台阶出口距离d分别为0 m,5m,10m。Ler图1坎式阶梯溢流堰模型示意图2.2数学模型本次研究采取了最为通用的RNGk-8模型，准确度高且应用领域广，适合用来本次对坎式阶梯溢流堰上的模拟计算。RNGk-8模型的控制方程如下。连续方程为：pui十=0atxi动量方程为：opu;Oxixjk方程为：a(pk)a(pu.k)十at8方程为：a(p8)a(pu;)十at式中：t为时间；ui、u,为速度分量；xi、x j 为坐标分量；G为由平均速度梯度引起的动能产生项；、别为密度和分子黏性系数；p为修正的压强；为紊流黏性系数；模型常数 Cis=1.42、C2 e=1.6 8、O,=0.7 17 9、0。=0.7 17 9。对于自由面的捕捉，本文采用VOF方法：.8.陕西水利Shaanxi WaterResources式中：F为体积函数，其最小值、最大值分别为0 和1,分别表示给定单元中全部为气和全部为水。采用结构网格，共2 90 110 个。上游人口的边界条件采用压力入口条件，设置六种水位，分别是110 m、115 m、12 0 m、130m、13 5 m、140 m。下游出口设置为自由出流条件，两侧边界和底部边界设置为不透水的边墙。本文中采用的三种泥沙粒径分别为1mm,6.9mm,10.5mm,其密度为3 6 5 0 kg/m，设置三种泥沙的比例为3：3：4,推移质系数为8,携沙系数为0.018。图2计算模型总体网格划分图3结果与分析3.1水流流态以下选取尾坎距台阶出口距离d为5 m的模型为例，模拟堰上水头H为2 5 m时坎式阶梯溢流堰的整体流动情况，分析水体流动的过程，时长为10 s。初始设置12 5 m深的水体，无进口流速，观察水体流动情况。图3 为d=0时坎式阶梯溢流堰在不同时段（每隔45 s）的流态图。pressure and vectors0.000E+00140.0Z-80.0-40.0FLOW-3D t=0.0 y=2.650E+01 ix=3 to 4 kz=3 to 7pressure and vectors-1.457E+055.409E+05140.0OuQuj十(u+uOx;Ox;从+x;xu+xNo.10October,2023aFaF+=0atx;98.0236.0X(a)t=0 sZxOx;ak+Gi-p88+Ciepx;0.000E+00374.0512.0max=3.82E+011.228E+06-80.01-40.0FL0W-3D t=45.003265 y=2.650E+01 ix=3 to 4 kz=3 to 7-1.618E+05140.0Z80.0-40.0FL0W-3D t=90.012100 y=2.650E+01 ix=3 to 4 kz=3 to 7650.098.0236.0X(b)t=45 spressure and vectorsmax=3.43E+015.329E+0598.0236.0X(c)t=90 s374.0374.0512.01.228E+06512.0650.0650.0第10 期2023年10 月-1.561E+05140.0Z-80.0-40.0FL0W-3D(=150.00984 y=2.650E+01 ix=3 to 4kz=3 to 7图3 坎式阶梯溢流堰不同时段的流态图由图4可知，在t=45s时，水流已铺满在溢流堰表面，并且水流无飞溅情况，初步接触到泥沙区域，沙面开始出现凹陷；当水流运行至90 s时，泥沙块的凹陷程度明显加大且水流趋于平稳；在第15 0 s时，水流冲刷最终形成一个相对稳定的冲刷坑，且在台阶面的水流也处于稳定流动状态。3.2消能率用上下游断面总能量之差与上游断面总能量的比值作为阶梯溢流堰的消能率，在这里选用的计算公式为：AEn=100%=EE式中：E，为上游过水断面的总水头；E,为下游过水断面的总水头。0.90.80.70.60.510图4不同尾坎位置在各堰上水头下的消能率由图4可以得出，坎式阶梯溢流堰的消能率随堰上水头H的增大而降低，堰上水头H增大，沿程水深增加，导致阶梯结构相对于水深的影响逐渐减小，因此，堰上水头H越大消能率越低。堰上水头H为2 5 m和3 0 m时，d分别取0 m、5 m、10m的消能率差别较小，但无论在何种工况下,d=10m时的消能率都比d=0时要低，消能效果更差，进而推测冲刷力度更大。3.3最大冲坑深度和位置水流会在泥沙表面产生大于泥沙临界起动力的床面切应力，致使泥沙起动，经过一段时间的冲刷，溢流堰下游形成冲刷坑。图5 为尾坎距台阶出口末端距离d为0 m、5 m、10 m 时，坎式阶梯溢流堰在堰上水头H为2 5 m时水流稳定后，冲刷时长为15 0 s时x-z方向的冲坑地形图。堰上水头H一定时，随着尾坎距台阶出口末端距离d的增加，冲刷形成的冲坑深度D将会变大，冲坑位置X也会逐渐后移。进一步分析不同堰上水头H条件下，不同尾坎位置的坎式阶梯溢流堰的最大冲坑深度D和位置X,见图6 和图7。由图6 可以看出，当尾坎位置相同时，随着堰上水头H的增大，沿程水头增加，导致阶梯结构相对于水深的影响逐渐减小，消能效果减弱，最大冲坑深度D也逐步增大；当堰陕西水利Shaanxi Water Resourcespressure and vectorsmax=3.26E+015.357E+051.228E+0698.0236.0X(d)t=150 sEi-E2100%d=0 md=5 md=10 m1520H/mNo.10October,2023上水头H达到15 m后，水头保持一定时，尾坎距台阶末端距离d=0时，冲坑的深度最小；d=0m时，冲坑的深度D最小；d=5m略大于d=10m的冲坑深度接近但远高于d=0时各堰上水头的冲坑深度。由图7 可以得出，尾坎位置相同时，随着堰上水头H的374.0512.02530650.035增大，冲坑x坐标逐渐增大而后移；研究尾坎位置对冲刷的影响，列向看图表，可以发现在堰上水头H达到15 m之后，流态与冲刷坑的形成稳定后，尾坎距台阶水平出口距离越大，冲坑的x位置坐标越大，即冲刷破坏力度越大。terrainelevation contours-42-6140.0-80.0-40.0-53(1)140.0Z-80.0-40.0-53140.0Z-80.0-40.0图5 不同距离时坎式阶梯溢流堰xz方向的冲坑地形图605040302011010图6堰上水头与尾坎位置对冲坑深度的影响12010080604010图7 堰上水头与尾坎位置对冲坑位置的影响2998.0236.0X(a)d=0 mterrain elevation contours-1498.0236.0X(b)d=5 mterrain elevation contours-1598.0(c)d=10 m1520H/m1520H/m65374.0512.02462374.0512.02462236.0374.0X2525100650.0100650.0100512.0650.0d=0 m-d=5 md=10 m3035d=0 md=5 md=10 m30359第10 期2023年10 月在这六种工况下，冲坑深度D的变化在水头较低（堰上水头H在10 m15m）的情况下，三个不同的尾坎位置由于水头过低，冲至泥沙的效果不明显，三个位置形成的最大冲深度和冲坑位置X差别不大。随着水头的加大，开始呈现一定的变化规律,d=5m和d=10m时的最大冲坑深度较为接近，同时变化幅度较d=O时也较大。而三种模型的冲坑位置X变化幅度相当。3.4冲坑长度图8 为不同尾坎位置条件下，冲坑长度随堰上水头H的变化。当尾坎位置相同时，随着堰上水头H的增大，冲坑长度也逐步增大；当堰上水头H一定时，尾坎距台阶末端距离d越大，冲坑的长度就越大。尾坎距水平出口距离d越大，冲坑长度的变化幅度就越大，在图中表示为冲坑长度增长越多，即形成的冲坑越深，冲刷破坏力度随尾坎位置d的增大而增大。因此，在实际的水利建设工程中，坎式阶梯溢流堰的尾坎布置在靠近尾坎出口位置，以减小水流的冲刷。从地形地貌图中可以看到冲坑形态的横剖面图，冲坑均呈曲面下凹态势，可以得出：随着上游堰上水头H的增加，最大冲坑深度D先随之增大，当达到一定深度后，基本不再随水头的增大而发生改变。尾坎位置的变化对冲坑长度和深度的影响较大，对冲坑位置X变化影响的幅度较小,d=O时的冲坑最小。所以在日常水利建设中，将尾坎设置在台阶末端的阶梯溢流堰最为广泛。14012010080604010图8 堰上水头与尾坎位置对冲坑长度的影响陕西水利Shaanxi WaterResources4结论（1）坎式阶梯溢流堰冲坑位置随堰上水头和尾坎位置距离的增大而后移。当尾坎位置相同时，随着堰上水头的增大，最大冲坑深度也逐步增大；尾坎在台阶水平出口时冲坑深度最小，冲刷破坏能力最小。（2）当尾坎位置相同时，随着堰上水头的增大，冲坑长度也逐步增大；当堰上水头一定时，尾坎距台阶末端距离d越大，冲坑的长度就越大。1徐荣峰.水利工程建设与生态环境可持续发展探讨J.城市建设理论研究，2 0 14,0 0 0（0 3 6):4447-4448.2A.I.瓦西列夫斯基.世界各国水电工程发展的经验J.水利水电快报，2 0 0 4,0 2 5（0 0 5）：2 9-3 3.3彭勇，张建民，许唯临，等.前置掺气坎式阶梯溢流堰掺气水深及消能率的计算J.水科学进展，2 0 0 9,2 0（0 1）：6 3-6 8.4陈剑刚，张建民，许唯临.前置掺气坎式阶梯溢流堰体型特点及工程应用试验研究J.四川大学学报（工程科学版）,2 0 10,42(06):6-11.5王业红，肖兴斌.阶梯溢流堰大单宽流量的水流特性及消能效果.长江职工大学学报，1999，（0 1）：45-49.6吴守荣，张建民，许唯临，等.前置掺气坎式阶梯溢流堰体型布置优化试验研究J.四川大学学报（工程科学版），2 0 0 8，（0 3）:3 7-42.7刘洋.高坝溢流堰掺气减蚀试验研究J.水利科技与经济，2 0 17,(1).8】严维，栗帅，赵德亮.阶梯消能工的研究与应用.价值工程，2013.9贾洪涛.尾坎位置对坎式阶梯溢流堰水力特性影响研究 .人民d=-0 md=5 md-10m1520H/mNo.10October,2023参考文献长江,2 0 19,5 0(0 0 7):18 9-19 4.10贾洪涛。台阶形式对阶梯溢流堰水力学特性的影响研究。中253035国农村水利水电，2 0 19,No.444（10）：15 2-15 6+16 1.11刘成林，陈宇豪.基于Flow-3d的水平射流冲刷泥沙数值模拟J.人民长江,2 0 16.（上接第6 页）2008,609(1):241-251.3结论2 Castro J J,Santana-Ortega S.A general theory on fish aggregation to（1）根据过鱼设施的类别和特点，针对水利、水电工程floating objects:An alternative to the meeting point hypothesisJ.高水头高坝，集运鱼系统较其他过鱼设施更优。Reviews in Fish Biology&Fisheries,2001.(2结合某水利枢纽工程，提出一种新型集运鱼系统的3饶振祺.彭水水电站工程环境保护措施实施情况及取得的效果设计，详细介绍了集鱼系统、提升转运系统、分选系统、运J.办公室业务，2 0 12（S2)：3.输系统和放流系统的设计方案。4 鲁玉军，王军，张楚.马岭水利枢纽工程集鱼系统设计与施工技术.云南水力发电，2 0 19,3 5（4)：4.参考文献5 Chandler J A.Feasibility of Reintroduction of Anadromous Fish Above1 Williams J G.Mitigating the effects of high-head dams on theor Within the Hells Canyon Complex.Columbia River,USA:experience from the trenchesJ.Hydrobiologia,:10.