蒋孜伟 ,许光祥
(1. 重庆交通大学 河海学院,重庆 400074;2. 重庆交通大学 水利水运工程教育部重点实验室,重庆 400074)
对未建大坝的天然河道进行流场模拟时,为了获得较好的来流边界,通常将进口选在较为顺直的河段,并且假定进口流速或单宽流量均匀分布。但当大坝建成后,计算域进口位置只能选择大坝断面,但发电、控泄发电、敞泄停机等不同工况下的电站泄流(即计算域进口来流边界)大不相同。特别是在开展引航道口门区通航水流条件研究时,由于引航道紧靠坝下,如仍采用无坝情况的均匀来流边界,必然会与实际情况造成较大误差。电站无论以何种方式泄流,距坝下一定距离必然会出现流速大小及分布恢复到无坝全断面泄流工况的断面,这个断面与大坝之间的距离为恢复距离,并将此断面定义为恢复断面。分析不同泄流方式下的水流恢复距离,对分析电站分散泄流对引航道口门区通航水流条件的影响具有重要意义。
引航道口门区水流呈现出复杂的三维特性[1-2]。水流恢复距离的大小受泄水闸的开启方式[3-9]、泄洪闸的位置[10]、电站的调节方式[11]、船闸的冲泄水方式[12]等影响,目前针对水流恢复距离等相关研究成果尚不多见。本文针对景洪电站发电、控泄发电、敞泄停机等工况下的水流沿程分布进行数值模拟研究,以探究水流恢复距离,进而对电站的调度运行提供参考。
景洪电站位于澜沧江干流下游,云南省西双版纳州首府景洪市北郊约5 km处,介于思茅港—中缅243号界碑5级航段上,是澜沧江中下游梯级开发规划中的第6级电站。该工程是以发电为主,兼顾航运,并具有防洪、旅游等综合利用效益的大型水利枢纽工程。水库正常蓄水位602.0 m,死水位591.0 m。枢纽主体工程包括:左右岸混凝土重力坝段、泄洪闸段、冲沙底孔、坝后式电站厂房、右岸垂直升船机等建筑物,船闸级别为4级,电站与升船机采用异岸布置于河床右侧。下游引航道紧邻河道右侧凸岸,长约601 m。最高通航流量为7 100 m3/s,最低通航流量为504 m3/s。当日均入库流量在405~504 m3/s时,电站下泄流量不小于405 m3/s;当日均入库流量小于405 m3/s时,按天然入库流量控制下泄。景洪电站坝址处于河道窄弯处,坝轴线位于两反向河道间的过渡段(图1),平面上呈烟斗形,下游引航道紧邻河道右侧凸岸,其左侧河床较高,岸坡较陡,水深较浅,枯水水花强烈,加之其上游受枢纽发电、泄洪等影响,引航道口门区的水流条件将更加复杂,因此有必要对其口门区的水流条件开展研究。根据《内河通航标准》[13]并结合本工程的船闸等级,船闸引航道口门区的水流表面最大流速限制如下:① 纵向流速不大于2.0 m/s;②横向流速不大于0.3 m/s;③ 回流流速不大于0.4 m/s。
图 1 景洪电站枢纽总平面布置(单位:m)Fig. 1 General layout drawing of Jinghong Hydropower Station (unit: m)
对于电站泄流边界的给定,采用较为简单的方法,进口单元节点切向流速Vτ均设为0,法向流速Vn则根据节点所在区域具体确定,如:
发电工况:只有电站厂房段才有水流,则设厂房段的节点Vn=V1,其他段节点的Vn=0。即图2(a) 中节点1~5 的 Vn=V1,6~13 的 Vn=0。
控泄发电工况:电站厂房和部分泄洪段有水流,则设厂房段的节点Vn=V1,开启的泄洪闸段节点Vn=V2,其他段节点的 Vn=0。即图 2(b)中节点 1~5 的 Vn=V1,6~9 的 Vn=V2,10~13 的 Vn=0。
敞泄停机工况:电站厂房、升船机和非溢流坝段无水流,泄流主要集中在溢流坝段,则泄洪闸段全部节点 Vn=V3,其他段节点的 Vn=0。即图 2(c)中节点 1~5 的 Vn=0,6~9 的 Vn=V3,10~13 的 Vn=0。
选用6个流量和13种边界共78组工况组合(表1)进行了流场模拟,用以分析恢复距离。
图 2 电站泄流边界处理示意Fig. 2 Diagram of discharge boundary processing of hydropower station
表 1 电站坝下非全断面泄流恢复距离计算工况Tab. 1 Calculation conditions of non-full section discharge recovery distance downstream of hydropower station dam
设部分泄流边界工况下某断面某节点的流速为u1i,无坝全断面泄流工况下同断面同节点的流速为u0i,则二者的均方差与全断面泄流工况下同断面平均流速V0的比值为:
均方差越小,数据越聚集,反之,数据越离散。当某断面流速的均方差5%时,则可认为部分泄流工况与全断面泄流工况的水流条件一致,该断面称为恢复断面,其距大坝的长度即为恢复距离L。
图3是Q=3 000 m3/s时典型工况的计算结果。从图3可知,泄流边界所占的比例越小,水流的恢复距离越长。根据5%的判断条件,可得各工况的恢复距离(表2)。
从表2和图4可以看出,恢复距离具有如下规律:
(1)两端为泄流边界的恢复距离均大于中间为泄流边界的恢复距离,整体比较而言,左端为泄流边界的恢复距离最长,这主要是由于受双侧或单侧扩散以及边界条件的影响。
(2)随着流量增大,动量增强,恢复距离越长,但其变化规律性不强,这主要是由于受天然河流的河床、河岸以及人工建筑边界的影响。
(3)泄流宽度占的比例越大,恢复距离越短,但其变化规律性不强,这与复杂的地下边界可能存在关系。
图 3 部分工况沿程变化情况(Q=3 000 m3/s)Fig. 3 changes along the river under partial conditions (Q=3 000 m3/s)
表 2 各工况的恢复距离Tab. 2 Recovery distance of each working condition
图 4 恢复距离随流量的变化Fig. 4 Recovery distance varies with flow
由表 2 可知,Q=504,1 500,3 000,5 000,7 100 和 8 967 m3/s时的恢复距离分别为 423.6,473.8,540.4,602.2,675.5和739.3 m。由表3可以看出,Q=504,1 500和3 000 m3/s时,=1.0和=0.2工况的最大纵流、最大横流及最大回流是基本一致的,此时的水流恢复距离小于引航道长度,所以无论上游泄流工况如何改变,引航道口门区的通航水流条件都不会受其影响。Q=5 000,7 100和8 967 m3/s时,=1.0和=0.2的最大纵流、最大横流及最大回流变化甚为明显,此时的水流恢复距离大于引航道长度,所以上游泄流工况的变化会引起引航道口门区的通航水流条件明显改变。
表 3 =1.0与=0.2的水流条件对比Tab. 3 Comparative analysis table of flow conditions as=1.0 and =0.2
表 3 =1.0与=0.2的水流条件对比Tab. 3 Comparative analysis table of flow conditions as=1.0 and =0.2
?
(1)本文基于平面二维水流数学模型计算了景洪电站不同泄流工况下的水流恢复距离,结果表明,两侧泄流的恢复距离均大于中段泄流的恢复距离,泄流宽度占的比例越大,恢复距离越短,泄流流量越大,恢复距离越长。通过对全断面泄流工况和极限泄流工况在引航道口门区水流条件的对比可知,当水流恢复距离大于引航道长度时,相对应的工况会造成引航道口门区水流条件的明显变化,此时需考虑电站泄流方式对下游通航的影响。
(2)计算工况相较于景洪电站的调度计划而言显得“极限”,这种“极限”体现在与电站的调度运行方式相比泄流边界更为集中,通航水流条件更为复杂,所以电站实际调度运行方式的水流恢复距离比计算工况得出的水流恢复距离要短,对通航的影响也应更微弱。同时本文计算得到的水流恢复距离可为景洪电站的调度运行是否对下游引航道口门区的通航水流条件造成影响提供参考。