1、第 卷第 期 年 月河 海 大 学 学 报（自 然 科 学 版）（）：基金项目：中央级公益性科研院所基本科研业务费专项（）；引江济淮工程科技创新课题（）；江西省交通运输厅科技项目（，）作者简介：缪婧娴（），女，硕士研究生，主要从事水力学及河流动力学研究。：通信作者：王晓刚（），男，正高级工程师，博士，主要从事水力学及河流动力学研究。：引用本文：缪婧娴，王晓刚，辛沛，等进鱼口水深不足时竖缝式鱼道水力特性及进出口允许最大水深差河海大学学报（自然科学版），（）：，（），（）：进鱼口水深不足时竖缝式鱼道水力特性及进出口允许最大水深差缪婧娴，王晓刚，辛 沛，唐柱闩（河海大学水利水电学院，江苏 南京；南京
2、水利科学研究院水工水力学研究所，江苏 南京；安徽省引江济淮集团有限公司，安徽 合肥）摘要：为研究竖缝式鱼道进鱼口水深不足时鱼道进出口允许最大水深差，构建了鱼道整体数学模型，对鱼道流速、流态、水面线等进行了数值模拟分析。结果表明：当进鱼口水深不足时，上游段水深与出鱼口水深一致，流速稳定；中间段水深开始缓慢下降，流速逐渐变大；下游段水深急剧下降，流速迅速变大，在进鱼口处出现跌水，流速出现最大值。当鱼道坡度相同时，出鱼口水深与进出口允许最大水深差呈线性关系，出鱼口水深越大，进出口允许最大水深差越大；当出鱼口水深相同时，鱼道坡度与进出口允许最大水深差也呈线性关系，鱼道坡度越陡，进出口允许最大水深差越小
3、。关键词：竖缝式鱼道；水深差；水位变幅；水面线；水力特性；整体模型中图分类号：；文献标志码：文章编号：（），（，；，；，）：，：；闸坝等工程措施破坏了河流连通性，严重威胁河流生态系统健康。鱼道作为一种保护鱼类资源的水利工程，具有为鱼类提供迁徙通道和连接河流上下游栖息地的重要作用，通过鱼道维持河流的连通性和生物河 海 大 学 学 报（自 然 科 学 版）第 卷的多样性能够有效防止珍稀鱼类的灭绝。由于水文情势改变或者水电站调度运行，鱼道上下游水位处于持续变化中。通常设计要求鱼道进鱼口水深大于出鱼口水深，但实际工程中进鱼口水深小于出鱼口水深的情况时有发生。大量工程实践表明，鱼道进鱼口水深不足时，鱼道
4、上游段水面平稳，但下游段水位会持续下降，会在进鱼口段发生明显跌落，甚至出现鱼类难以逾越的跌水。竖缝式鱼道是应用最为广泛的鱼道布置形式之一，已在池室尺寸布置、鱼道坡度和池室水力特性方面取得了相对成熟的研究成果。尺寸布置方面，等较早针对 种不同布置形式的竖缝式鱼道进行了水力特性研究，结果表明流量和池室内水深呈线性关系，并且水流以射流形式通过竖缝，当鱼道池室长为 倍竖缝宽度、宽为 倍竖缝宽度时鱼道的消能效果最好；等利用二维 模型模拟分析了偏转角度、竖缝尺寸、隔板尺寸等影响因素对池室内水流结构的影响。鱼道坡度方面，等对、和 种坡度的竖缝式鱼道进行了水流结构分析，结果表明，当鱼道坡度为 时池室内水流流态
5、具有明显的二维特性，为 和 时池室内水流流态表现为三维特性；曹庆磊等和刘志雄等也通过物理模型试验和数值模拟分析发现，在鱼道坡度小于 时竖缝式鱼道内垂向流速变化不大，具有明显的二维特性；郭维东等研究发现，竖缝式鱼道常规池室底板坡度影响主流的消能速率，底板坡度越大消能速率越大。池室水力特性方面，等利用超声波多普勒测速仪检测了水流结构和量化速度分布，结果发现流场可划分为以最大速度为特征的主流区和由低速与水平涡流确定的回流区；等利用二维数学模型模拟了竖缝式鱼道池室内水流流态，采用结构性与非结构性相结合的四边形网格并基于浅水方程的二元解对竖缝式鱼道的流场结构进行了分析。鱼道上下游的运行水位会直接影响过鱼
6、季节的适宜过鱼条件，如水电枢纽机组调度极易造成鱼道上下游水位变幅过大的问题，使得鱼道上下游水位组合不合理，不利于鱼类上溯。等基于 模型对不同鱼道入口深度的流场分布进行了定量分析，确定了适合鱼类通行的鱼道进鱼口水深条件。等利用分离涡（）模型模拟了进鱼口流态，发现低水位时产生的漩涡对鱼类上溯有不利影响。段鸿锋等对某水电站鱼道进口布置进行了优化，解决了低水位时没有较好的诱鱼水力条件的问题。魏萍等提出了一种优化的堰孔组合形式以解决进鱼口水位不衔接问题。金志军等从吸引流和竞争流角度出发，提出了一种鱼道进口概化模型。竖缝式鱼道消能效果较充分，能适应较大水位变幅，一般用于能适应较复杂流态的大、中型鱼类。此外
7、，设计多个进鱼口也可以较好地适应大水位变幅情况，如丰满水电站尾水最大变幅为.，根据下游水位变幅，丰满水电站过鱼设施布置了不同高程的 个进鱼口，其中进口最低水深为.；犍为航电枢纽鱼道进口运行水位为.，为适应水位变幅较大的特点，设置了 个进鱼口。目前国内外关于鱼道沿程水面线变化的研究较少，且数值模拟研究主要集中于模拟 或 个标准池室的水流结构，缺乏对鱼道进出口允许水深差的研究，不能适应实际工程中大水位变幅的需要，且鲜有关于进鱼口水深不足时鱼道运行水位组合的研究。为此，本文依托鱼道水位变幅达 的某工程开展进鱼口水深不足时竖缝式鱼道水力特性及进出口允许最大水深差研究。数学模型及验证某鱼道布置在电站厂房
8、的左侧边坡上，全长，设有 个进鱼口和 个出鱼口。主要建筑物包括进鱼口、普通池室、休息池和出鱼口等。主要过鱼对象是异齿裂腹鱼、双须叶须鱼和拉萨裂腹鱼等。本文以该鱼道池室细部结构为原型构建数学模型。控制方程鱼道坡度小于 的竖缝式鱼道水流结构具有明显的二维特性，因此本文基于 软件建立二维鱼道水动力模型，其中水动力控制方程由连续性方程和垂向积分的、方向动量方程组成：（）|（）（）|（）第 期缪婧娴，等 进鱼口水深不足时竖缝式鱼道水力特性及进出口允许最大水深差|（）（）|（）式中：为水深，；为水位，；为时间，；、为平面坐标，；为谢才系数，；为重力加速度，；为风摩擦系数；、为风速及其在、方向上的分量，；为
9、科氏力，；为大气压强，；为水的密度，；、为有效剪切力分量，；、分别为、方向的单宽流量，（），其中?，?，?、?分别为、方向上沿水深的平均流速，。根据水动力控制方程，利用有限体积法求解即可得到每一时刻在（，）处的水位、水深以及流速?、?。模型域与网格剖分模型长度应包含鱼道中水深不相同的所有池室长度，模型模拟结果才能够准确反映原型鱼道池室水力特性。为模拟全部水深影响范围，模型总长，包含 个普通池室和 个休息池。以休息池为间隔，每 个池室组成一个池室段。池室自进鱼口向出鱼口以 编号，净长，净宽，竖缝宽度为。休息池净长为普通池室的 倍，净宽与普通池室一致。为了更直观地揭示水面线变化特征，暂不考虑转弯段
10、的设置。鱼道进口门槽处桩号，鱼道整体和池室布置见图。图 鱼道布置示意图 边界条件可以是水位控制条件或者流量控制条件，为探究鱼道不同上下游水位组合对进出口允许最大水深差的影响，模型边界出鱼口采用上游水位、进鱼口采用下游水位控制。鱼道边壁及底部采用混图 非结构化网格示意图 凝土浇筑，故糙率系数为定值.。模型计算时间步长为.，总步数为 步，计算时间总长为 ，干湿判定条件选用软件默认值，不考虑地球自转产生的科氏力的作用，模型为冷启动。模型采用非结构化网格（图），在保证计算精度的前提下，对模型进行了网格剖分，采用 网格，整个模型节点数为，网格数为。模型验证 利用所依托工程物理模型试验资料对进鱼口水深不足
11、时数学模型进行验证，验证工况如表 所示。由于流量和水位是一一对应的关系，物理模型采用流量水位边界。鱼道物理模型底坡为，包含 个池室。表 数学模型验证工况 工况出鱼口水位 进鱼口水位 出鱼口水深 进鱼口水深 鱼道下泄流量（）工况说明鱼道运行最小水深鱼道运行中间水深鱼道运行上游最大水深河 海 大 学 学 报（自 然 科 学 版）第 卷图 鱼道沿程水深试验值与模拟值对比 水深验证鱼道内水位高低是影响水流结构的关键因素。沿程水深验证结果（图）显示，水面线沿程下降，数值模拟结果和试验值基本一致。工况 上下游水深差为.，模拟值和试验值的误差最小，最大误差小于.；工况 最大误差小于.；工况 最大误差小于.。
12、模拟误差较大的区域均集中在进鱼口段，考虑到该段水位发生大幅下降，水面有明显波动，因而测量时会产生一定的误差。其余位置处误差均小于。流态流速验证典型池室平面流态见图，可见池室中存在两种不同的流态，一种连通着两级竖缝的主流区，一种是流速较小的形成封闭流线的回流区。主流区明确，鱼类可以感知上溯方向；有一定面积的回流区，保证鱼类有休息区域。主流在池室中整体呈 形，分布在池室左侧 的范围之间。回流区流速普遍低于 ，主流左侧回流区面积较小，为逆时针流动；右侧的回流区面积较大，为顺时针流动。数值模拟的回流区范围略大于物理模型试验得到的回流区范围，可能是因为包含了几乎静水的回流区域，而物理模型试验存在一定的人
13、为误差。图 鱼道池室内流态对比 表 竖缝流速试验值与模拟值对比 工况池室编号水深 竖缝流速（）试验值模拟值工况 工况 工况 不同工况和池室水深时对应的竖缝流速对比见表，可见试验值和模拟值吻合较好。以工况 为例，水流从 号池室到 号池室的过程中，池室水深从.下降至.，物理模型试验结果表明池室竖缝流速从.增大至.，数值模拟结果显示池室竖缝流速从.增大至.，模拟误差为.。水面线沿程变化和竖缝流速验证结果表明模型参数设置合理，模拟结果可以很好地反映实际水流结构。模拟结果与分析 鱼道沿程水力特性鱼道上下游水位变化会引起池室水深的壅高或者降低。竖缝流速是判断鱼类能否成功上溯的关键，为揭示水面线和竖缝流速沿
14、程变化规律，模拟了出鱼口水深均为、进鱼口水深分别为（工况）、（工况）和（工况）种工况鱼道池室的水流结构。水面线工况、水面线和水深沿程变化如图 所示。由图（）可知，进鱼口水深不足时，水面线沿程变化可分为 段：上游段水面线和鱼道底坡基本平行；中间段水深开始缓慢下降；下游段水深迅速减小，越靠近进鱼口相邻池室水深差越大，水面坡降陡增，形成跌水。以工况 为例，上游段为距离进口门槽 以外区域，此段池室水深与出鱼口水深一致，维持在.不变；中间段为距离进口门槽 范围区域，此段池室水深从.缓慢降至.，在长达 的鱼道范围内水深差仅为.；下游段为距离进口门槽 以内区域，此段池室水深从.降至.，在 的鱼道范围内水深差
15、达.，水面陡降，因此需第 期缪婧娴，等 进鱼口水深不足时竖缝式鱼道水力特性及进出口允许最大水深差重点关注下降段是否会形成流速屏障，阻碍鱼类上溯。此外，由图（）可以看出，上游段水深不变，无明显波动。工况 因进鱼口水深不足，水深影响的池室长度范围（水深低于出鱼口水深）约为 ，工况 和工况 水深影响的池室长度范围约分别为 和。可见进出口水深差越大，水深影响的池室长度范围越大，跌水越显著。图 鱼道水面线和水深沿程变化 流速工况、各池室竖缝流速沿程分布见图，可见 种工况鱼道内流速顺着水流方向总体上呈现增大趋势。竖缝流速沿程变化曲线可分为 段：上游段很稳定，流速在小范围内波动，与进鱼口水深无关，种工况下竖
16、缝平均流速均为 ，最大流速达到 ；中游段流速开始逐渐变大；下游段变化较大，因水位大幅度下降，竖缝流速越来越大，在进鱼口处流速出现最大值。因此，进鱼口水深不足时，出鱼口水深是鱼类上溯的关键。当出鱼口水深一定时，进出口水深差存在一允许最大值，超过该值，鱼道进鱼口将出现流速屏障甚至跌水，阻碍鱼类上溯。竖缝断面弗劳德数 沿程变化如图 所示，进鱼口水深不足时，弗劳德数顺着水流方向沿程增大；水深稳定时，弗劳德数稳定在；水深下降时弗劳德数增大，且水深下降得越快，弗劳德数增幅越大。但弗劳德数全程始终小于，因此竖缝处水流为缓流时有利于鱼类上溯。图 竖缝流速沿程变化 图 竖缝断面弗劳德数沿程变化 鱼道进出口允许最大水深差在鱼道坡度为 不变的情况下，通过数学模型试算出鱼口水深一定时，不同进鱼口水深组合，并逐渐减小进鱼口水深，从而得到进鱼口允许最小水深，以此来探究进出口允许最大水深差。四大家鱼的爆发游泳速度为 ，且我国竖缝式鱼道设计流速一般不大于 ，相对于竖缝宽度，鱼类体宽较小，故以竖缝剖面 区域流速小于 为判断依据，模拟得出进鱼口允许最小水深。鱼道坡度为 时不同上下游水位组合对应的进出口允许最大水深差试算结