异侧竖缝式鱼道水力特性的数值模拟_李洋.pdf

第 卷第期 年月水电能源科学 ：异侧竖缝式鱼道水力特性的数值模拟李洋，韩雷，田振华，狄高健，李树航，叶昆河（黑龙江省水利科学研究院，黑龙江 哈尔滨 ；黑龙江大学水利电力学院，黑龙江 哈尔滨 ）摘要：竖缝式鱼道因其能够适应不同水位变化、消能效果明显、水流流态稳定等特点而逐渐引起水利工程界的重视。通过数值模拟的方法探究了池室长宽比对异侧式竖缝鱼道水力特性的影响，结果表明，池室长宽比的变化对单位体积耗散率影响较小；主流区沿程衰减率先随池室长宽比的增大而增大，后保持不变。池室长宽比越大，池室内最大流速越靠近竖缝处；不同池室长宽比，池室主流流速与池室最大流速（主流流速沿程变化最大值）比值在一定范围内。池室长宽比在 之间，可在鱼道主流区内获得较好的水流流态、较大的回流区及较好的消能效果，有利于研究区内鱼类的洄游。关键词：池室长宽比；异侧式；竖缝鱼道；水力特性中图分类号：文献标志码：文章编号：（）收稿日期：，修回日期：基金项目：黑龙江省自然科学基金联合引导项目（）作者简介：李洋（），男，工程师，研究方向为水力学及河流动力学，：1引言在水利工程建设过程中，由于人为截断水流，迫使具有洄游习性的鱼类无法洄游繁衍，导致鱼类减少甚至灭绝，因而建立一些过鱼工程帮助鱼类洄游繁衍十分必要。过鱼设施包括鱼闸、升鱼机、集运渔船和鱼道，目前应用最广的过鱼工程是鱼道，按照形式可分为槽式鱼道、池室鱼道和特殊形态的鱼道，池室鱼道主要包括溢流堰式、竖缝式、底孔式，其中竖缝式鱼道能够适用于不同的水位变化，具有流速、紊动能变化较小的特点，适用于大多数鱼类洄游，。目前，关于竖缝式鱼道的研究成果较丰富，但多为同侧竖缝式鱼道，而关于异侧竖缝式鱼道在整体结构设计和细节优化上的研究较少，关于实际鱼道工程运用的研究更少。现有关于异侧竖缝式鱼道的研究，在进行模型试验时，并未考虑具体区域鱼类的极限流速，从而针对具体鱼类提出适宜的参数范围。鉴此，本文充分考虑具体区域的鱼类极限流速，通过数值模拟的方法研究池室长宽比对异侧竖缝式鱼道水流结构的影响，定量分析主流流速沿程轨迹线、主流流速沿程衰减规律、竖缝附近最大流速、回流区分布规律，给出适合研究区域内的四大家鱼上溯的池室长宽比，并与已有研究结果进行优劣比较，以期为鱼道设计提供参考依据。2鱼道结构国内某水利工程设计鱼道总长为 ，坡比:，为便于计算分析，以其中 长的鱼道为例进行分析研究。考虑到工程所在区域的 鱼 类 以 四 大 家 鱼 为 主，感 应 流 速 为 ，极限流速为 ，因而设计流速采用 。由文献 可知，当竖缝宽度与池室宽度之比为 之间时，鱼道池室内可获得较为理想的水流流态，竖缝宽度 ，池室宽度 ，。长隔板长 ，横向导板长 ，竖缝相对位置（竖缝中心与较近一侧池室边壁距离与池室宽度的比值）为 。隔板和导板厚均为 。以池室长度为研究变量，其变化范围在 之间，每增加 为新的池室长度，共选取 种不同池室长度进行研究，对应的池室长宽比范围为 ，具体尺寸见图。3数学模型建立3.1网格剖分及边界条件采用 边 长 为 的 规 则 六 面 体 划 分 网图鱼道尺寸布置（单位：）格，其中节点 个，网格 个。具体剖分情况见图。图鱼道网格剖分 进口 采 用 流 量 为 边 界 条 件，流 量 为 ；出口采用水深为边界条件且设置为压力出口，水深为 ；边墙及底部设置为无滑移固壁边界。追踪自由液体表面采用 方法，通过计算水和气的体积分数来表征物体的形态，其控制方程为：（）（）式中，为水的体积分数；为气体的体积分数；为时间；为方向流速；为方向的坐标。3.2基本控制方程鱼道采用专业软件 标准模型并耦合 模块对水流自由表面进行捕捉，；其中三维水流模型的控制方程如下。连续方程为：（）动量方程为：（）（）紊流模型方程为：（）（）（）（）（）（）（）（）紊动能为：（）（）紊动黏性系数为：（）式中，为水的密度；、均为流速分量；、均为坐标分量；为压力；、分别为分子黏性系数和紊动黏性系数；为重力加速度在方向的分量；为紊动能；为紊动能耗散率；为平均速度梯度引起的紊动能产生项；、均为模型参数，分别取为 、。迭代计算的各项残差值均为，时间步长设置 。3.3模型验证为验证数学模型参数设置的合理性和网格剖分的正确性，以有机玻璃为材料，制作 组比尺为:的隔板物理模型。其糙率为 （对应原型糙率为 ），与实际混凝土工程糙率略有不同。考虑到鱼道中竖缝附近流速最大，水头损失较大，池室边壁影响较小，可忽略糙率不同带来的影响。以第、个隔板之间的池室左下角为原点，选取沿着水流方向距离原点 处为测量断面，在断面上等距离选取 个点（测量流速位置对应原型水面以下 处）测量流速（流速通过无线旋浆流速仪测定），具体测点布置见图，对比物理模型实际测量结果与数值模拟结果可知，选 取 断 面 个 测 量 点 流 速 拟 合 误 差 在 范围内，说明数值模拟参数设置正确，网格剖分合理。在此基础上，对 种长宽比的鱼道进行数值模拟计算，为简洁起见，本文只水电能源科学 年第 卷第期李洋等：异侧竖缝式鱼道水力特性的数值模拟图测点分布（单位：）选取具有代表性的模拟结果进行分析说明。4池室水力特性分析4.1池室流态分析竖缝式鱼道池室内存在两种流态，即由于竖缝射流而出现的位于池室中间位置的主流区流态和主流区两侧的回流区流态，两侧回流区流速方向与主流区流速方向相反，主流区是鱼类在上溯过程中的主要通道，因而主流弯曲程度对鱼类的洄游影响较大。根据计算结果（图），当 时（为池室长，为池室宽），相邻池室主流区流态衔接处弯曲程度较大，曲率半径较小，扭曲明显，在隔板竖缝处的主流方向变化剧烈，单个池室主流区流态呈一斜线型，相邻池室流区流态呈现型。主要是因为池室长宽比较小时，竖缝射流并未在池室中充分消散便迅速进入下一竖缝。当 时，单个池室主流区流态呈半抛物线型，相邻池室呈现 型，池室长宽比（）越大，型越明显，这与文献 的结论一致；图不同池室长宽比下流速分布图 当 时，主流区出现贴壁流现象，鱼类在上溯过程中会受伤，因而在设计时池室长宽比宜小于 。从相邻池室主流区衔接处弯曲程度及避免出现贴壁流现象，在设计过程中，池室长宽比宜在 之间。4.2池室最大流速分布与流速沿程衰减规律池室最大流速是决定过鱼效果的关键参数之一，通常出现在竖缝附近（图），其不能超过鱼类的极限流速，具体位置和大小取决于池室结构设计。池室的最大流速也影响鱼道消能效果，消能效果通常以单位水体消能率作为衡量指标：（）（）式中，为水密度，取为 ；为竖缝宽度，取为 ；?为各池室最大流速平均值；为池室宽度，取为 ；为池室长度，根据池室长宽比确定；为隔板厚度，取为 ；为流速系数，取为 。图为不同池室长宽比最大流速与单位水体消能率变化规律。由图（）可知，池室最大流速随着池室长宽比的增大而增大，当 时，流速增长较快，池室长宽比每增大 ，最大流速增加 ；当 时，流速增长较慢，池室长宽比每增大 ，最大流速仅增加 。由图（）可知，单位水体消能率随着池室长宽比的增大先增后减，在池室长宽比 时，单位水体消能率最大。主要是因为当 时，虽然池室长度逐渐增大，但池室最大流速增加更快；当 时，池室长度逐渐增大，但主流区水体与池室边壁接触，出现了贴壁流现象，消能效果较好，因而池室最大流速增加较慢，单位水体消能率逐渐减小。图不同池室长宽比最大流速变化和单位水体消能率 分析不同池室长宽比下的主流区流速沿程衰减率（），具体变化规律见表，其中 为池室主流流速沿程变化的波谷值，为池室主流流速沿程变化的波峰值。由表可知，当 时，主流区流速沿程衰减率随池室长宽比的增大而增大；当 时，衰减率基本保持不变。主要是因为池室长宽比较小时，竖缝射流未与池室水体充分扩散便迅速进入下一竖缝，使得流速降低较小；池室长宽比较大时，竖缝射流与池室水体较充分混合，使得流速在池室中央明显降低。表主流区最大流速沿程衰减率 池室长宽比（）池室长宽比（）注：池室长宽比 后，沿程衰减率不再变化，因而后面部分未列出。本次设计流速为 ，因而池室内最大流速不宜大于 ，即池室长宽比宜小于 ；考虑到主流流速沿程衰减率，池室长宽比不宜小于 ；随着池室长宽比的变化，单位水体消能率先增大再减小，但均小于 ，可忽略其对池室水力特性的影响。4.3主流流速沿程变化与主流轨迹以上游隔板与右侧池室边界交点为原点，在不同池室长宽比下，垂直水流方向上每隔 选取一个横断面，每个断面选取 个等距离的点进行流速分析，对每个断面的最大流速进行连接，即得到主流流速沿程变化曲线，见图（）；对每个断面的最大流速的位置进行连接，即得到主流轨迹线，见图（）。由图（）可知：不同池室长宽比，主流流速沿程变化曲线略有不同，当 时，主流流速先增大，后逐渐降低，即将进入下一竖缝时，流速快速增加，即出现一次波峰和一次波谷；当 时，主流流速先增大，后逐渐降低，距离选取原点 处，主流流速出现第一次波谷，随后主流流速又开始增加，并再次减小且出现第二次波谷，在进入下一竖缝前迅速增加，即出现两次波峰和两次波谷。与第一次出现的波谷相比，第二次出现的波谷流更小，这与董志勇等在 时，前半池主流流速逐渐减小，后半图主流流速沿程变化曲线和主流轨迹线 池主流流速逐渐增加不同，主要是因为本文鱼道设计时增设了横向隔板，主流区又出现了较大范围贴壁流现象。主流刚贴壁时，流速便出现第一次波谷，贴壁中间位置流速出现波峰，之后水流到达下一竖缝前，逐渐降低，使得池室主流流速整体发生改变。随着池室长宽比的增大，池室主流流速波峰处流速先明显增加后趋于不变，波谷流速先缓慢增加后趋于不变（保持在 左右，见图（）中右侧斜线部分），这与同侧竖缝式鱼道得到的逐步减小趋势相反。池室最大流速并未出现在竖缝处，而是偏下游一些，这与文献 结论一致，因而池室最大流速要比竖缝处流速略大，且池室长宽比越大，池室最大流速（池室主流流速沿程最大值）位置越靠近竖缝处，主要是池室长宽比越大，进入下一竖缝前水流越平稳，进入竖缝时流速迅速增加，因而位置靠近竖缝处。池 室 主 流 流 速 与 池 室 最 大 流 速 比 值 均 在 之间，因此测量出竖缝附近池室最大流速，便知道主流流速沿程分布情况，进而确定主流流速是否符合鱼类上溯。由图（）可知：沿着水流方向距离原点 范围内，不同池室长宽比，池室最大流速的位置基本不变，超出该范围，池室长宽比越大，最大流速的位置越偏离水流入口边壁（即轴）。主要是因为距离竖缝较近时，竖缝位置射流效果相同，但当距离竖缝越远，竖缝射流作用越弱，随着池室长宽比的增加，池室后半程越远离竖缝，因而主流不会迅速进入下一竖缝，而是在池室内充水电能源科学 年第 卷第期李洋等：异侧竖缝式鱼道水力特性的数值模拟分扩散。当 时，主流轨迹出现平行于池室边壁的情况，随着池室长宽比的增加，主流轨迹随之延长，流速轨迹越向水流入口对面池室边壁靠拢。主流轨迹平行池室边壁部分的长度随着池室长宽比的增加而增加，从 增至 时，平行池室边壁部分的长度由 增加到 。分析可知，当 ，主流区出现贴壁流现象，平行池室边壁的主流轨迹越长，鱼类在上溯过程中受伤的概率越大。4.4回流区水力特性分析回流区的面积分布、最大流速分布影响着鱼类的短暂休息和上溯效果，因此有必要分析回流区的水力特性。表（表中，左、右分别为左、右侧回流区面积；左、右分别为左、右侧回流区最大流速）为不同池室长宽比两侧回流区分布规律和最大流速分布规律。由图、表可知：池室内部存在两个回流区，即靠近上游隔板回流区（简称左侧回流区）和靠近下游隔板回流区（简称右侧回流区）。随着的增大，左侧回流区面积（左）逐渐减小，右侧回流区面积（右）逐渐增大，当 时，两侧回流区大小相当。对于不同池室长宽比，其左右两侧回流区的最大流速都在逐渐增大，左侧回流区最大流速（左）为 ，右侧回流区最大流速（右）为 。表不同池室长宽比的回流区参数 池室长宽比（）左左（）右右（）左右 一般两侧回流区面积相差不大，鱼类可在两侧回流区内休息，考虑到回流区流速不能超过鱼类的感应流速，否则鱼类难以持续上溯，池室长宽比宜小于 。5结论 不同池室长宽比下，单位水体消能率均较小；池室长宽比越大，池室内最大流速越靠近竖缝处；主流流速沿程与池室最大流速（主流流速沿程最大值）的比值均在 之间。当 时，竖缝两侧相邻池室主流区流态衔接处弯曲程度较小，且不会出现主流贴壁现象；当 时，主流流速沿程衰减率最大，且池室最大流速小于设计流速。综合考虑，当 时，最适合研究区内鱼类上溯洄游，这与同侧竖缝式鱼道池室长宽比宜在 范围内可获得良好的流态有一定差异。参考文献：王琲，杨文俊，陈辉生态水工建筑物 鱼道的建设及研究进展人民长江，（）：董志勇，冯玉平，异侧竖缝式鱼道水力特性及放鱼试验研究水力发电学报，（）：董志勇，江录兵，毛斌，等异侧布置竖缝式鱼道紊流结构试验研究水生态学杂志，（）：张国强，孙双科竖缝宽度对竖缝式鱼道水流结构的影响水力发电学报，（）：，：张羽，张凌峰，杜庆君，等异侧竖缝宽度对鱼道水流结构的影响研究水力发电，（）：赵彬如，戴会超，戎贵文，等竖缝位置对竖缝式鱼道水 力 特 性 的 影 响 水 利 水 电 科 技 进 展，（）：，李秀萍，周济人，傅朝康，等竖缝宽度对双侧竖缝式鱼道水力特性的影响研究中国农村水利水电，（）：徐体兵，孙双科竖缝式鱼道水流结构的数值模拟水利学报，（）：曹庆磊，杨文俊，陈辉异侧竖缝式鱼道水力特性试验研究河海大学学报（自然科学版），（）：（下转第 页）由表可知，各模型下预测值与实际值之间误差均相对较小，表明机器学习算法应用于水库排沙预测有效，在一定程度上可实现综合考虑不同影响因素的水库排沙预测。其中，算法模型的平均绝对误差 为 ，均方根误差 为 均最小，同时结合其决定系数最高，表明针对万家寨水库排沙预测过程，在现有数据条件下，算法模型预测精度优于其他三种算法模型。4结论 利用万家寨水沙系列数据，基于不同机器学习算法，建立综合考虑入库流量、入库含沙量、出库流量、坝前水位、坝前水位差、累计淤积量影响的水库排沙预测模型。不同机器学习算法模型得到的预测出库含沙量与实际出库含沙量的分布情况基本一致，预测出库含沙量与实际出库含沙量之间相关性良好，各模型决定系数均在 以上。基于 算法所建立模型决定系数最高为 ，且平均绝对误差 为 ，均方根误差 为 均最小，表明相对于其他模型，算法模型在水库排沙预测方面具有更高的准确性和精度。参考文献：王梦雅，刘丽冰，熊桂龙，等面向袋式除尘器的大数据挖掘 优化算法研究 电子测量与仪器学报，（）：周鑫，谢晖，付山，等基于 算法的中心带孔圆板拉深翻孔变形方式的研究 锻压技 术，（）：刘泉声，王栋，朱元广，等支持向量回归算法在地应力场反演中的应用岩土力学，（增刊）：吴芳，李映雪，张缘园，等基于机器学习算法的冬小麦不同生育时期生物量高光谱估算 麦类作物学报，（）：，（，；，）：，：；（上接第 页），（，；，）：，；，；（），：；水电能源科学 年