程 艳，吴礼国，徐 红，李 维

(四川省交通勘察设计研究院有限公司，四川 成都 610017)

枢纽在运行过程中，闸坝的调度运行方式是根据上游来流量及电站发电需求进行调节，一般运行方式为小流量时调度闸坝控制下泄流量，船闸通航，流量增大至分界流量后闸坝敞泄，船舶停航。闸坝的调度运行方式不同，枢纽上下游的水流条件会有所不同，而船闸引航道口门区水流特性是影响船舶航行安全和通航效率的重要因素[1]。

船闸通航主要集中在小流量情况下，此时枢纽通过开启不同闸孔数及闸孔高度，控制下泄流量，维持坝前正常蓄水位，使电站正常发电。而泄水闸有多种开启方式，包括集中、分散、均匀、区段、局部(均匀)或分段间隔等[2]。泄水闸不同的开启方式，使河道水能量再分配，改变水流流速和流向，造成回流范围和强度的差异，影响船闸引航道及口门区的水流和航行条件，尤其是下引航道。本文以数学模型为依托，研究风洞子航运枢纽工程泄水闸调度运行方式对船闸下引航道口门区的通航条件影响。

1 模型建立与验证

1.1 模型基本原理

河道水流一般可视为不可压缩牛顿流体，其水流运动规律可用Navier-Stokes方程组描述，对该方程进行时均化得到其Reynolds形式，包括质量守恒方程和动量守恒方程[3]。从Reynolds方程出发，考虑到工程所在河段水平尺度远大于垂直尺度，水力参数在垂直方向的变化明显小于水平方向的变化，可忽略垂向加速度，设定压强服从静水分布；不计垂直方向的流动时间和空间的微分，将方程沿水深积分，得到平面二维的河道水流运动控制方程[4]。

水流连续方程：

(1)

水流动量方程：

(2)

(3)

式中：Zw为水面高度；Zb为床面高度；t为时间；H为水深；q1、q2为x、y方向的单宽流量；qm为单位面积的旁侧入流流量(为正)或出流流量(为负)，式(1)假定水流密度在整个模拟区域内保持不变[5]；β为各向同性的动量修正系数；ρ为水流密度；pa为水面大气压力；Ω为科里奥利力参数；τbx、τby分别为x、y方向的河床剪切应力；τsx、τsy分别为x、y方向的表面风剪切应力；τxx、τxy、τyx、τyy为由紊流引起的剪切应力。

1.2 模型建立与验证

本文选用风洞子航运枢纽工程坝轴线以下1.8 km范围内的实测地形，采用三角形网格对计算区域进行离散，既可以较好地描述河道边界和河道中建筑物边界，又能达到较高精度，网格最大长度约20 m，模型中对船闸引航道、电站、闸坝区等部位进行网格加密，网格最小长度约2 m。研究区域网格划分见图1。

图1 研究区网格划分

为保证模型参数选取的合理性，在模型建立后对其进行验证，验证资料为实测流量1 150 m3/s时沿程瞬时水位，验证结果见图2。可以看出，模型水位与原型最大水位偏差为0.06 m，差值较小，建立的模型可做进一步研究。

图2 模型与原型沿程水位

2 方案拟定

2.1 风洞子航运工程常规调度方式

渠江风洞子航运枢纽工程泄洪冲沙闸布置于主河槽，为开敞式宽顶平底闸。共布置20孔泄水冲沙闸，从左岸至右岸为1#～20#，单孔净宽15 m，堰顶高程230.0 m，闸顶高程260.6 m。本枢纽为径流式日调节运行，枢纽运行调度方案主要考虑电站本身的日调节作用和枢纽冲排泥沙的要求。枢纽拟定的调度运行方式为：当入库流量小于发电引用流量843 m3/s时，水库水位在242.6 m(死水位)～243.0 m(正常蓄水位)之间消落发电；当入库流量大于843 m3/s且小于6 500 m3/s时，电站引用流量以外多余流量通过冲沙闸下泄，水库水位维持正常蓄水位243.0 m发电；当入库流量大于6 500 m3/s时，水库闸门逐步开启泄洪冲沙，电站停机，枢纽平面布置见图3。

2.2 研究方案

图3 风洞子航运枢纽工程平面布置(单位：m)

根据枢纽调度运行方式，选取流量小于6 500 m3/s、泄水闸控泄时的4级流量进行研究，此时电站满发，电站下泄流量843 m3/s。为全面分析泄水闸开启方式对船闸下游引航道口门区及连接段水流条件的影响，从闸坝开启位置和闸坝开启孔数两方面进行研究，各方案及边界条件见表1。

表1 风洞子航运枢纽工程闸坝调度方案

3 结果与分析

在船闸下引航道口门区及连接段布置17个断面，共148个测点(图4)，用以观测各工况条件下测点范围内的最大纵向、最大横向及回流流速，探讨泄水闸不同开启方式对下引航道口门区及连接段的通航水流条件影响。根据《船闸总体设计规范》[6]规定，Ⅰ～Ⅳ级船闸口门区的最大流速限值为：平行于航线的纵向流速vy≤2.0 m/s，垂直于航线的横向流速vx≤0.30 m/s，回流流速v0≤0.4 m/s。

3.1 闸门开启位置对船闸下游水流条件的影响

图4 研究区测点布置

流量3 000 m3/s、闸门开启不同位置时下游河道的流速流态分布见图5。可以看出，闸门集中开启位置不同，下引航道口门区及连接段的流速流态有所不同。当集中开启左岸1#～4#闸门时，下游流场如图5a)所示，下泄水流集中在左岸，水流出闸孔后与电站尾水汇合顺流而下，右侧有大范围的回流。水流在导航堤头开始向右偏转，口门区及连接段大部分处于回流区，流速较小，河道断面流速分布不均匀，主流偏右岸。当集中开启中间9#～12#闸门时，下游流场如图5b)所示，下泄水流集中在中部，两侧各形成一个小范围的回流区，口门区有小范围的回流，流速较小。当集中开启右岸17#～20#闸门时，下游流场如图5c)所示，开启闸孔远离电站，下泄水流集中在右岸，在河道中间形成一个小范围的回流，两股水流在导航堤处汇合后顺流而下，主流仍在右岸，口门区及连接段大部分处于回流区，流速较小。

图5 Q=3 000 m3/s时风洞子航运枢纽下游流速流态分布

开启不同位置闸门时船闸下引航道口门区及连接段测点流速对比见图6。各级流量下集中开启左岸闸孔的最大流速与开启中间闸孔相差不大，开启右岸闸孔的较小，最大流速值都没有超过规范限值，见图6a)；集中开启左岸闸孔的横向最大流速小于集中开启右岸闸孔，开启中间闸孔的最大，其中当流量等于6 000 m3/s，开启中间和右岸8孔闸门时，横向最大流速超过规范规定的限值，最大值为0.328 m/s，见图6b)；各级流量下的回流流速相差不大，开启左岸和中间闸孔时基本一致，开启右岸时稍大，但都在规范限值以内，见图6c)。

图6 不同闸门开启位置测点流速对比

综合分析，当流量小于6 500 m3/s，风洞子闸坝集中开启左岸闸门时，横向流速整体较小，没有超过规范限值，下引航道口门区及连接段的流速流态更有利于船舶安全航行。

3.2 闸门开启数量对船闸下游水流条件的影响

流量4 500 m3/s、左岸闸门开启不同数量时下游河道的流速流态分布见图7。可以看出，下游流速流态基本都满足船舶安全航行的要求。但闸门开启数量不同，下游河道及口门区的回流范围有所不同，流速流态也有所差别。下泄水流集中在左岸，水流出闸孔后与电站尾水汇合顺流而下。泄水闸下游右侧河槽是回流区，主流在导流墩处穿入航道，从船闸口门区以下逐渐向右过渡，口门区及连接段的航道左侧处于回流区，连接段最大回流流速为0.293 m/s。

图7 Q=4 500 m3/s风洞子航运枢纽下游流速流态分布

开启左岸不同数量闸门时船闸下引航道口门区及连接段测点流速对比见图8。可以看出，各级流量下开启闸门数量越多，测点范围内的最大流速和横向最大流速越大，回流流速越小。当流量为6 000 m3/s，开启8孔闸门时最大流速为1.329 m/s，横向最大流速为0.298 m/s，接近规范限值；开启9孔闸门时最大流速为1.391 m/s，横向最大流速为0.316 m/s，超过规范限值要求。造成这种现象的主要原因为水流出闸孔后与电站尾水汇合顺流而下，开启闸门越多泄水闸单宽流量越小，闸孔下游水位越平稳，与电站出水水流的相互抵冲作用越弱，此时电站出流对船闸引航道口门区及引航道的水流条件影响越明显，从而造成开启闸门越多，测点范围内水流流速和横向流速较大。

图8 不同闸门开启数量测点流速对比

综合分析，当流量小于6 500 m3/s，风洞子闸坝在满足闸孔过流条件时，开启闸门孔数越少，下游引航道口门区及连接段的最大流速和横向流速越小，没有超过规范限值，水流条件更有利于船舶安全航行。

4 结论

1)风洞子航运枢纽工程闸门调度运行方式不同，下游引航道口门区及连接段的水流条件不同。相同流量级下，集中开启左侧闸门时，下游引航道口门区及连接段的水流条件最优。相比于集中开启中间闸孔和右侧闸孔，横向流速更小，更有利于船舶安全航行。

2)相同流量级下，在满足闸孔过流条件时，开启左侧闸门孔数目越少，下游引航道口门区及连接段的最大流速和横向流速越小，水流条件更有利于船舶安全航行。

3)在枢纽设计过程中，可通过建立数学模型，不断试验闸孔调度方式对下游引航道口门区及连接的水流条件影响，为初步拟定闸坝调度运行方式、保证船舶顺利安全通航提供依据。