1、452景洪水电站水力驱动式升船机输水系统数值仿真分析计算摘 要：本文在景洪水力驱动式升船机输水系统供水管路布置型式的基础上，根据水流、平衡筒及承船厢运行原理，建立输水系统数值仿真模型，仿真输水管道及竖井充泄水流量过程，计算分析了输水系统的水力特性及承船厢的运行特性。本文还计算比较不同的平衡筒与竖井间隙情况下，承船厢出入水时竖井水位响应过程及承船厢运行特性，从而为为该种型式升船机的设计和应用提供参考依据。关键词：水力驱动式升船机 输水系统 仿真 1 前言景洪水电站为澜沧江干流中下游河段规划梯级中第六级，是兼顾发电、航运的综合性水利枢纽工程。景洪水电站航运过坝建筑物的规模按级航道、300t 级船型
2、设计，远期考虑 500t 船只通航过坝，其通航运行条件是为：上游最高通航水位 602.00m、最低通航水位 591.00m；下游最高通航水位544.90m，下游最低通航水位 535.58m。年通航 310 天，每天工作 21 小时，年货运量 1358000t（双向）。通过多方经济技术比较，通航过坝方式采用水力驱动升船机。水力驱动式升船机是一种新型的通航过坝建筑物，在安全性及机构的简化方面跟其他型式升船机相比具有较大的优势，但目前在国内外尚无这种型式的升船机投入实际运行，无实际经验可循，由于“水力驱动”是该升船机与其他型式升船机的主要区别，因此建立输水系统的水力学非恒定流数学模型，对升船机输水过
3、程进行数值仿真就显得尤为必要，通过对不同方案输水系统的仿真计算分析和比较，可以初步研究输水系统水力特性及升船机的运行特性，为升船机输水系统的合理布置提供科学依据。由于仿真模型具有参数可调的特点，通过该仿真系统还可研究不同竖井与平衡筒之间的间隙大小对升船机运行的影响，为该型式升船机的关键参数的确定提供参考数据，并与物理模型试验相互验证和补充。2 输水系统数值仿真模型2.1 模型的建立2.1.1 输水系统能量方程 竖 井， , 4， Q， A， 3， 2， 1， Q， A， 0， v，（ 下 游 水 库 ）上 游 水 库图 2.1 景洪升船机输水系统概化图升船机输水系统概化如图 2.1，先考虑竖井
4、充水（承船厢下降）的情况，根据能量守恒 Bernouli方程可得：453242421200 .)( sv gAQgAQhH（1）. 410 dtLdtLtgA其中：H上游水库水位； h竖井水位；A 各管道计算断面面积； As竖井横截面积；L各管道段惯性长度；Q 各管道段的流量； 各管道段的阻力系数； 阀门阻力系v数。由于系统供水管路按等惯性布置，在阻力系数相同情况下，分流口可按平均对称分流处理，方程（1）可整理为：24240216208 )(.)( sv gAQgAgAhH（2）dtLL401304.为统一和方便起见，可令竖井输水流量： 4Q阻力项系数： 240216208.)(),( gAg
5、Afv惯性项系数： 401304.),( LLf方程（2）可写为：（3）dtQAfgAfhHs),(2),(2同时考虑竖井充、泄水的情况，方程（4）可整理为：（4）tLffs),(1),(2其中正号表示充水，负号表示泄水。2.1.2 水流、平衡筒及承船厢系统整体运动方程水力驱动式升船机的原理就是通过向竖井内充泄水，增大或减小平衡筒所受的浮力来驱动平衡筒运动，进而带动承船厢系统下降或上升（见图 2.2） ，如果不考虑钢丝绳的重量及弹性作用，整个系统的运动遵循牛顿第二定律，分别对平衡筒和承船厢作运动受力分析可得：（5）cppffcpw MgKXhMgAdtX4)1()1(2)(42 454其中：
6、竖井水位； 平衡筒底面位置高程； 单个平衡重质量； hXpMcM该钢丝绳承担的船箱质量。 钢丝绳受力； 作用在平衡筒上浮力； 水体质量密gFfFw度。 系统摩擦系数，按系统阻力为 100T 估算，并以承船fK厢及平衡筒的重量作为参考力， 取 0.0121。fK令 cpwMgA4cppff g)1()1(2方程（7）可统一写为：（6） )(2Xhdt2.1.3 水流连续方程由图 2.2 可知，竖井水流需满足连续性条件，则竖井内输水流量可用下式表示：（7） dthAdtQpss )()(其中： 竖井横截面积； 平衡筒底面积； sAQ竖井输水流量。方程（7）可整理为： （8） dthAdtXQpsp
7、)(2.1.4 输水系统数值仿真模型联合输水系统水流能量方程（4） 、整体运动方程（6）及连续方程（8） ，可建立输水管道、平衡筒及承船厢系统的整体仿真模型，通过仿真计算，可以得到输水主管道及竖井的供水流量响应过程，同时还可以模拟水力驱动下的平衡筒及承船厢系统的运行过程。（9）)( ),(21),(2 2XhdtdthAtQdtQALfgfhHpsps承 船 厢平 衡 筒 竖井FfFgMcMphX图 2.2 水流、平衡筒及承船厢系统受力分析图4552.2 阻力系数景洪升船机输水管道系统按等惯性布置，通过各级分流口向 16 个竖井供水，按分流口可划分为5 个输水管道段，由各管道段的沿程阻力系数、
8、局部阻力系数和惯性长度，就可以确定输水系统的阻力项系数及惯性项系数。3 升船机输水系统运行仿真计算3.1 输水系统布置及运行参数景洪升船机输水系统设计方案和修改方案输水系统布置及整体运行参数见表 3.1。两方案的主管路及第一级分管路直径均为 2.2m。设计方案第二至第四级分管路直径从 1.6m 逐级减小至 0.8m，竖井横截面为 6.3m* 6.3m（边长 *边长）的正方形，四周导圆半径为 1.5m，横截面积为 37.76m2，平衡筒与竖井间隙为 0.2m；修改方案第二至第四级分管路直径均为 1.6m，竖井横截面采用直径为 6.9m 的圆形，横截面积为 37.39m2，平衡筒与竖井间隙为 0.
9、1m。上游最高通航水位：602.00m，上游最低通航水位：591.00m；下游最高通航水位： 544.90m，下游最低通航水位： 535.58m。表 3.1 升船机输水系统布置及整体运行参数名称 参数 设计方案 修改方案各级管道直径(m) 2.2、 2.2、1.6 、1.2、0.8 2.2、2.2、1.6、1.6、1.6阀门启闭方式 线性匀速启闭 线性匀速启闭输水管路系统 阀门启闭时间（s ） 60 60横截面形状 6.3 6.3-R1.5（边长边长，四角导圆半径为 R） 6.5 圆形横截面积（m 2） 37.76 33.18竖井数量 16 16横截面形状 5.95.9-R1.3 （边长边长，
10、四角导圆半径为 R） 6.3 圆形横截面积（m 2） 33.36 31.17重量（t） 344 344平衡筒数量 16 16横截面积（m 2） 69.112.0=829.2（长宽）净重（t） 800正常水深（m） 2.5m正常水重（t） 1800误载水深（m） 2.5 0.2m承船厢最大水重（t） 19443.2 两方案输水系统仿真计算结果根据输水系统布置及运行参数，再联合求解输水系统数值仿真模型，在上游最高通航水位602.00m 及下游最低通航水位 535.80m 的运行条件下，模拟了升船机设计及修改方案的充、泄水过程，图 3.1（a） 、 （b）为设计方案竖井充、泄水水力特性曲线及承船厢运
11、行速度曲线；图 3.2（a） 、 （b）为修改方案竖井充、泄水水力特性曲线及承船厢运行速度曲线，其仿真计算统计结果见表 3.2。456表 3.2 两种方案输水系统充、泄水仿真计算结果设计方案 修改方案方案工况特征参数 充水 泄水 充水 泄水主管道最大输水流量(m 3/s) 37.70 35.49 40.24 37.20输水管道最大流速（m/s） 9.9 9.3 10.6 9.8平衡筒最大淹没水深(m) 8.59 8.09 9.22 8.66正常运行中 0.062 0.059 0.076 0.072竖井水面升降速度(m/s) 承船厢入(出)水 0.248 0.242 0.384 0.359正常运
12、行中 0.125 0.119 0.152 0.140承船厢最大运行速度(m/s) 承船厢入(出)水 0.047 0.067 0.069 0.094阀门开启 0.009 0.005 0.008 0.007阀门关闭 -0.002 -0.003 -0.002 -0.003承船厢最大运行加速度(m/s2) 承船厢入(出)水 -0.060 0.065 -0.046 0.073总运行时间(s) 687 781 561 649充泄水体积（m 3） 20693.6 18049.2从表 3.2 可以看出：（1）设计方案输水系统主管道最大充、泄水流量分别为 37.70 m3/s 和 35.49 m3/s，管道中最
13、大流速分别为 9.9m/s 和 9.3m/s，相应竖井水位上升、下降速度分别为 0.062 m/s 和 0.059 m/s，承船厢最大运行速度分别可达 0.125 m/s 和 0.119 m/s；输水系统充、泄水运行总时间分别为 687s 和781s，一次充泄水体积约为 20693m3；阀门启闭过程中承船厢运行加速度最大为 0.009 m/s2，承船厢入（出）水时其运行加速度随之骤然减小或增大，其峰值分别为-0.06 m/s 2 和 0.065 m/s2 。图 3.1(a) 竖井及主管道充水水力特性曲线（设计方案）00.050.10.150.20.250 100 200 300 400 500
14、 600 700 800t(s)竖井水面上升速度(m/s)05101520253035400 100 200 300 400 500 600 700 800t(s)Q0(m3/s)5205305405505605705805906006100 100 200 300 400 500 600 700 800t(s)高程(m) 竖 井 水 位平 衡 重 底 部 高 程承 船 厢 位 置 高 程-0.0200.020.040.060.080.10.120.140 100 200 300 400 500 600 700 800t(s)承船厢下降速度(m/s)457图 3.1(b) 竖井及主管道泄水水力
15、特性曲线（设计方案）图 3.2 (a) 竖井及主管道充水水力特性曲线（修改方案）5205305405505605705805906006100 100 200 300 400 500 600 700 800 900t(s)高程(m) 竖 井 水 位平 衡 重 底 部 高 程承 船 厢 位 置 高 程00.050.10.150.20.250.30 100 200 300 400 500 600 700 800 900t(s)竖井水面下降速度(m/s)05101520253035400 100 200 300 400 500 600 700 800 900t(s)Q0(m3/s)-0.0200.0
16、20.040.060.080.10.120.140.160 100 200 300 400 500 600 700 800 900t(s)承船厢上升速度(m/s)5205305405505605705805906006100 100 200 300 400 500 600t(s)高程(m) 竖 井 水 位平 衡 重 底 部 高 程承 船 厢 位 置 高 程-0.0500.050.10.150.20.250.30.350.40.450 100 200 300 400 500 600t(s)竖井水面上升速度(m/s)0510152025303540450 100 200 300 400 500 6
17、00t(s)Q0(m3/s)-0.0200.020.040.060.080.10.120.140.160 100 200 300 400 500 600t(s)承船厢下降速度(m/s)458图 3.2(b) 竖井及主管道泄水水力特性曲线（修改方案）（2）修改方案主管道最大充、泄水流量分别为 40.24m3/s 和 37.20 m3/s，管道中最大流速分别为 10.6m/s 和 9.8m/s，相应竖井水位上升、下降速度分别为 0.076 m/s 和 0.072 m/s，承船厢最大运行速度分别可达 0.152 m/s 和 0.140 m/s；输水运行总时间分别为 561s 和 649s，一次充泄水
18、体积为18049.2m3；阀门启闭过程中承船厢瞬时加速度最大为 0.008 m/s2；承船厢下降入水时瞬时加速度峰值为负值-0.046 m/s 2，承船厢上升出水 时其峰值为正值 0.073 m/s2。由于锥形阀门阻力系数较大，正常输水运行中两个方案竖井水位最大升降速度均小于 0.1m/s，承船厢最大运行速度小于 0.2 m/s 。与设计方案相比，由于增大了第三、四级分水管路的管径，修改方案输水能力较设计方案略有提高，系统输水总运行时间有所缩短。 两方案正常运行过程中（空气中）承船厢加速度均小于 0.01 m/s2，满足运行要求，承船厢入（出）水时其最大瞬时加速度则超过 0.01 m/s2 ，
19、下一节将对竖井与平衡筒的间隙面积比对承船厢入（出）运行的影响进行初步的计算分析。4 输水系统不同间隙面积比计算比较从两种方案（设计方案和修改方案）输水系统仿真计算分析结果来看，正常运行工况下承船厢速度曲线平滑下降，但在承船厢入（出）水时其速度骤然下降或上升，这是由于承船厢受力状态突然发生变化（浮力突然增加或消失）而引起其加速度的突变，这就要求竖井水位必须随之迅速上升或下降，从而改变平衡筒所受浮力来调节整个系统达到动态平衡，所以在承船厢入（出）水时竖井与平衡筒间隙内的淹没水深能否随之快速变化成为系统是否能在短时间内稳定下来的关键。为研究竖井与平衡筒不同间隙大小情况下水体对承船厢入（出）水时系统突
20、然 “失稳”的响应能力，在修改方案的基础上(竖井直径 =6.5m)，取间隙面积（ ）与竖井横截面面积（ ）的fDAsA比值 作为控制参数，将间隙大小 L 从 0.01m 逐步增至 0.8m，间隙面积比 A / As 相应从sA/0.006 增至 0.432，对修改方案输水系统充、泄水过程进行了仿真计算，表 4.1 给出了承船厢入（出）水时系统水力及运行特征值统计结果，图 4.1图 4.2 给出了竖井泄水时输水系统水力特性曲线及承船厢运行特性曲线。5205305405505605705805906006100 100 200 300 400 500 600 700t(s)高程(m) 竖 井 水
21、位平 衡 重 底 部 高 程承 船 厢 位 置 高 程00.050.10.150.20.250.30.350.40 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s)05101520253035400 100 200 300 400 500 600 700t(s)Q0(m3/s)-0.0200.020.040.060.080.10.120.140.160.180 100 200 300 400 500 600 700t(s)承船厢上升速度(m/s)459表 4.1 不同间隙面积比承船厢入（出）水系统水力及运行特征值充水 泄水平衡筒直径Dp(m)间隙大小L(
22、m)间 隙面积比A/As竖井水面上升速度(m/s)承船厢运行加速度(m/s2)系统振荡时间(s)竖井水面下降速度(m/s)承船厢运行加速度(m/s2)系统振荡时间(s)6.48 0.01 0.006 0.436 -0.02 1.0 0.438 0.022 3.46.46 0.02 0.012 0.428 -0.025 3.8 0.427 0.039 4.16.4 0.05 0.031 0.441 -0.047 6.0 0.4 0.052 7.46.3 0.1 0.061 0.384 -0.068 8.8 0.36 0.073 12.26.1 0.2 0.119 0.325 -0.081 12.
23、2 0.301 0.08 22.05.7 0.4 0.231 0.226 -0.092 18.0 0.228 0.081 44.85.3 0.6 0.335 0.175 -0.102 20.0 0.18 0.077 71.04.9 0.8 0.432 0.14 -0.104 24.8 0.15 0.07 106.0从计算统计结果可以看出：（1）竖井充水承船厢下降入水时，承船厢速度突然减小，竖井与平衡筒的间隙面积比越大，承船厢的负加速度越大，间隙面积比从 0.006 增大为 0.432 时，承船厢瞬时加速度从-0.02m/s 2增大至-0.104 m/s2；相应竖井水面上升速度也瞬时增大，但间隙
24、面积比从 0.006 增大为 0.432 时，竖井水面瞬时上升速度却从 0.436m/s 减小至 0.14m/s，这反过来表明间隙面积比越小，竖井水位越能迅速上升增大平衡筒浮力，从而抵消承船厢入水突然所受浮力的影响。（2）竖井泄水承船厢上升出水时，承船厢运行速度突然减小, 间隙面积比从 0.006 增大为 0.231时，承船厢瞬时加速度从 0.022m/s2增大至 0.081 m/s2；间隙面积比继续增大至 0.432，承船厢瞬时加速度略有减小，处于 0.070.081m/s 2之间；竖井水面瞬时下降速度则随间隙面积比的增大而逐步减小，间隙面积比从 0.006 增大为 0.432 时，竖井水面
25、瞬时下降速度从 0.438m/s 减小至 0.150m/s，同样表明间隙面积比越小，竖井间隙内水流对承船厢系统状态突变的响应能力越强。从图 4.1图 4.2 还可以得知，在承船厢出水后，随着间隙面积比的增大，竖井水位下降速度曲线、承船厢上升速度及加速度曲线逐步出现上下振荡的调整过程，且间隙面积比越大，其初始振幅越大，持续的时间也越长。竖井泄水承船厢出水时，随着间隙面积比从 0.006 增大为 0.432，系统从开始振荡至最终稳定的时间从 1.0s 增大为 106s，特别是间隙面积比大于 0.119 后，其波动时间急剧增大，振荡特征也愈加明显，对承船厢系统的运行稳定的影响也越大。5 结论与建议通
26、过对景洪升船机设计及修改两种方案输水系统进行数值仿真分析和比较，同时以竖井与平衡筒间隙面积与竖井面积的比值作为控制参数，在修改方案的基础上，计算比较了不同间隙面积比对升船机运行特性的影响，得出主要结论及建议如下：（1）景洪水力驱动式升船机输水系统能保证承船厢的稳定运行，竖井及设置动滑轮的平衡筒的布置型式可适应上下游通航水位变幅较大的条件。（2）升船机输水系统修改方案的主管道最大充、泄水流量分别为 40.24m3/s 和 37.20 m3/s，输水能力略高于设计方案；承船厢最大运行速度相应分别可提高至 0.152 m/s 和 0.140 m/s，升船机总运行时间则分别缩短为 561s 和 649
27、s；一次充、泄水最大耗水量 18049.2m3，较设计方案有所减小。两种方案升船机正常运行过程中承船厢最大加速度均小于 0.01 m/s2，但承船厢入（出）水时其最大瞬时加速度则超过 0.01 m/s2。（4）对升船机修改方案不同间隙面积比情况下的仿真计算表明：间隙面积比越大，承船厢出、入水时速度瞬时增大或减小得越大，即瞬时加速度峰值就越大；间隙面积比越小，承船厢出、入水时竖460井水面下降或上升速度就越大，竖井间隙内水流对承船厢系统状态突变的响应能力越强。同时，在承船厢设计布置型式上，工作阀门线性开启，间隙面积比大于 0.119 后，竖井泄水承船厢出水后其速度有较明显的振荡现象，竖井水位下降
28、速度也随之发生上下波动。因此，为减小承船厢入（出）水时瞬时加速度峰值，在设计及施工允许的前提下，竖井与平衡筒的间隙应选取较小值为好。（5）由于承船厢底部设计型式为平底，承船厢出、入水时其所受浮力骤然减小或增大是造成系统运行振荡的原因，故可考虑将承船厢底部设置为圆锥形等型式，从而逐步改变承船厢出、入水时其所受浮力，减缓其速度的突然变化率，同时还可考虑阀门变速开启方式（如先慢开至船厢出水后再快开），从而控制减小船厢出水后的瞬时加速度，下一步可在这两方面作进一步的研究（6）由于仿真模型的阻力系数等是根据升船机输水系统布置型式选取的经验值，尚需要物理模型试验来确定和校核，最终结果也应以试验结果为准。0
29、0.050.10.150.20.250.30.350.40.450.50 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s) A/As=0.00600.050.10.150.20.250.30.350.40.450 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s) A/As=0.01200.050.10.150.20.250.30.350.40.450 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s) A/As=0.03100.050.10.150.20.250.30.350.
30、40 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s) A/As=0.06100.050.10.150.20.250.30.350 100 200 300 400 500 600 700t(s)竖井水面下降速度(m/s) A/As=0.11900.050.10.150.20.250 200 400 600 800t(s)竖井水面下降速度(m/s) A/As=0.23100.020.040.060.080.10.120.140.160.180.20 200 400 600 800t(s)竖井水面下降速度(m/s) A/As=0.33500.020.040.
31、060.080.10.120.140.160.180 200 400 600 800t(s)竖井水面下降速度(m/s) A/As=0.432图 4.1 不同间隙面积比竖井水面下降速度曲线（泄水）461参考文献1须清华、张瑞凯等通航建筑物应用基础研究中国水利水电出版社，1999 年 4 月2刘亚辉、罗家麟船闸廊道复合管非稳定流的数值计水利水运科学研究，3（1） ，19973董曾南、余常昭著.水力学. 清华大学水力学教研组编高等教育出版社， 1995（第 4 版）4肖田元，张燕云，陈加栋. 系统仿真导论.北京：清华大学出版社，2000。-0.01-0.00500.0050.010.0150.020
32、.0250 100 200 300 400 500 600 700t(s)承船厢加速度(m/s2) A/As=0.006-0.02-0.0100.010.020.030.040.050 100 200 300 400 500 600 700t(s)承船厢加速度(m/s2) A/As=0.012-0.03-0.02-0.0100.010.020.030.040.050.060 100 200 300 400 500 600 700t(s)承船厢加速度(m/s2) A/As=0.031-0.06-0.04-0.0200.020.040.060.080 100 200 300 400 500 600
33、 700t(s)承船厢加速度(m/s2) A/As=0.061-0.06-0.04-0.0200.020.040.060.080.10 100 200 300 400 500 600 700t(s)承船厢加速度(m/s2) A/As=0.119-0.08-0.06-0.04-0.0200.020.040.060.080.10 200 400 600 800t(s)承船厢加速度(m/s2) A/As=0.231-0.08-0.06-0.04-0.0200.020.040.060.080.10 200 400 600 800t(s)承船厢加速度(m/s2) A/As=0.335-0.08-0.06-0.04-0.0200.020.040.060.080 200 400 600 800t(s)承船厢加速度(m/s2) A/As=0.432图 4.2 不同间隙面积比承船厢上升加速度曲线（泄水）
