浮堤局部三维水流特性试验研究

黄 卫，钱红露，於思瀚, ，倪玉芳，段文刚

(1.长江科学院，湖北 武汉 430010；2.长江设计集团有限公司，湖北 武汉 430010；3.武汉大学 水利水电学院，湖北 武汉 430072)

浮堤常被用于水电站拦污导漂[1-3]、进水口不良流态破除[4]等水利工程中。若将浮堤推广至河道整治，相较于坝体，具有如下优点：（1）浮堤对河流泥沙输移的影响较小。天然河道中水流垂向流速“上大下小”，泥沙分布“上小下大”，在分汊河道治理中使用坝体（如丁坝）难以实现只增加水流分流比而不增加分沙比的目标，而浮堤调节水流能力在水面较强，在近底部分较弱，因此对泥沙输移的影响小。（2）浮堤对来水来沙条件的适应性较好。因为浮于水面，所以河床冲淤对浮堤的影响较小；而坝体通常为永久建筑物，不仅达成整治目标需要一定时间，而且在预期目标达成后，若来水来沙条件发生重大改变，坝体可能产生不利后果。（3）浮堤对原有河床及生态环境的影响较小。坝体的修建及其局部防冲措施都会显著改变河床基质[5]，从而改变水生生物栖息地环境。综上，在河道整治工程中，浮堤具有广阔的应用前景。然而，由于缺乏对浮堤作用机理的全面认识，浮堤在河道整治工程中的推广应用受到了严重制约。

目前，针对坝体作用下三维水流结构的研究较多[6-10]，而浮体作用下水流结构的研究多集中于表面流场，缺少对其局部三维水流结构的研究。童中山等[11]对水电站导漂建筑物研究现状进行了综述，指出了导漂浮排在设计和实际运行中存在的较多问题，尤其需对水流结构进行深入研究；段文刚等[4]通过物理模型试验，在研究旋涡形成原因的基础上，采用浮堤破除平面回流、不对称进流和侧向进流等不利流态引发的立轴旋涡；黄明海[12]为了防止钉螺吸附在漂浮物上通过大型引水工程扩散，采用浮排进行导漂从而避免漂浮物进入引水渠道，并应用数学和物理模型分别对浮堤作用下的表面水流结构和漂浮物运动轨迹进行研究，初步探明浮排的布置型式、来流流量和引水流量等因素对表面流场的影响。本文针对浮堤局部水流结构特征及其关键影响因素作用规律开展水槽试验研究，为浮堤在河道整治中的推广应用提供参考。

1 试验方法

试验在整体可变坡水槽中进行。水槽长14 m、宽0.3 m、高0.25 m，入口放置整流格栅以平顺水流，出口设置合页式尾门（图1）；侧壁及底板均由厚3.6 mm 的透明玻璃拼接而成；支撑系统与地面软连接，不同部件之间柔性连接，均可有效减小因供水泵工作引起的振动；侧壁粘贴方格纸，以准确读取浮堤上下游水位。采用高分辨率PIV 系统测量流速，该系统主要由图像采集系统和流场计算系统组成。根据试验布置，试验段流场采集窗口距水槽进口9 m，距水槽出口5 m，进、出口长度满足紊流充分发展和消除水流扰动的要求。

图1 水槽系统示意Fig.1 Experimental setup in the flume

试验研究不同几何尺寸（长度L，淹没深度d，厚度B）浮堤作用下的水流结构，试验工况见表1。试验中浮堤固定在左岸，不同工况浮堤布置见图2。试验水槽宽度W=0.3 m，底坡坡比为0.001。进口处水流为恒定流，流量为0.002 88 m3/s，水深h=0.05 m，计算断面平均流速为0.19 m/s。沿水槽宽度方向设置流场采集断面，沿水槽宽度方向布置的10 个测量断面如图2(a)所示，编号为1#～10#，断面到玻璃水槽左边（浮堤所在）壁的距离分别为2.5、5.0、7.5、10.0、12.5、15.0、17.5、20.0、22.5、27.5 cm。

表1 试验工况汇总Tab.1 Summary of the experiment cases

图2 测量断面及浮堤布置Fig.2 Experimental setup of the measured cross-sections and the floating dikes（unit: cm）

试验测量参数包括浮堤附近水位和局部纵剖面二维流速场，其中二维流场为主流方向和垂向上的瞬时流速u与w。根据二维流场测量数据计算垂向紊动强度与雷诺应力，以反映局部水流特性及浮堤对水流的影响。

（1）垂向紊动强度计算式：

式中：RSMV为垂向紊动强度（m/s）；w′为脉动流速（m/s）；RSM为均方差；n为统计的有效样本数。

（2）雷诺应力计算式：

式中：RS为 雷诺应力（Pa）；ρ为水流密度，取1 000 kg/m3。雷诺应力是由于流场中流速分布不均而产生的，雷诺应力越大流速分布越不均匀，即紊动越剧烈。

2 试验结果与分析

2.1 浮堤局部水流结构特征分析

通过对比无浮堤（工况1）和布置浮堤（工况2，基本工况）的水流结构，揭示浮堤对水流结构的影响。无浮堤（工况1）工况下，测量断面的流线几乎平行，近似于均匀流动，流速分布规律遵循一般水槽内流速分布规律。图3 为工况2 中浮堤周围的水流表面流态，由于布置了浮堤，表面水体过流断面束窄，表面水流流经浮堤后偏向右岸，并向下游延伸，这一现象的本质是流经浮堤的水体发生了边界层分离。

图3 工况2 浮堤局部水流表面流态Fig.3 Photograph of the surface flow feature in Case 2

图4 展示了工况2 各断面流速分布，可见流场发生了剧烈变化：（1）由于浮堤下方发生了边界层分离，紧邻浮堤下游出现了回流区，回流区沿水槽宽度方向长10 cm，与浮堤长度（10 cm）相符。（2）回流在断面2#和3#处最为显著；靠近水槽边壁和浮堤右侧（图4(a)和图4(d)）的回流较小，这是由于断面1#受到了水槽边壁边界层的影响，也表现在断面1#较断面2#浮堤下方流速更小；断面4#正好位于浮堤非近岸端，该断面受浮堤的影响产生回流，又受浮堤右侧水体流速较大的影响，使回流减弱。（3）回流内最大的负流速为0.05 m/s，明显小于主流断面平均流速0.18 m/s，这表明回流处水体流动缓慢。（4）由于浮堤对水流的阻挡作用，紧邻浮堤上游水体流速较小，约为0.1 m/s，靠近浮堤壁面处流速仅0.01 m/s。（5）浮堤和浮堤下游回流使得过水断面减小，浮堤下方和右侧水体流速均较大，其中距离左岸5.0～10.0 cm 之间断面下方和浮堤右侧（断面5#）水体流速最大达0.25 m/s。（6）浮堤及回流影响下，断面流速分布云图出现了分层现象（图4(e)），尽管图4(f～j)中流线近似平行，但是结合图3 可知，回流的影响具有三维特性，水流经过浮堤后向右岸偏移，左侧为回流，故右侧水体流速加快。（7）断面6#～8#之间位于浮堤上游底部水流流速较小，与上层水体存在差异，这表明受到了浮堤的影响，而断面9#、10#之间水流逐渐恢复为均匀流动，这说明工况2 中浮堤的影响范围为20.0～22.5 cm。

图4 工况2 典型纵剖面流速分布Fig.4 Distributions of velocity in Case 2

图5 展示了浮堤右端（W=10 cm）附近断面的垂向紊动强度及雷诺应力分布，可以观察到垂向紊动强度与雷诺应力的分布规律在定性上吻合，即垂向紊动强度较大的位置其相应的雷诺应力绝对值也较大。浮堤下游回流核心部分（图4）的紊动强度达到0.035 m/s，其雷诺应力绝对值达到1 Pa 以上（图5）。值得注意的是，尽管浮堤及浮堤下方边界层的存在抑制了垂向紊动，但是流向（X方向）上的紊动较为剧烈，因此浮堤正下方的垂向紊动强度虽然较小，但是雷诺应力较大且绝对值达到0.5 Pa。

图5 工况2 浮堤附近断面垂向紊动强度及雷诺应力分布Fig.5 Distributions of vertical turbulent intensity and Reynolds stress around the dike in Case 2

2.2 浮堤淹没深度影响

图6 为工况3 和4 浮堤附近断面流速分布。与图4 相对比可见，随着浮堤淹没深度比D增大：（1）浮堤下游回流的垂向尺度相应变大；（2）浮堤右侧（断面5#）水体流速明显增大，D为1/5（工况3）时最大流速达到0.22 m/s，D为2/5（工况2）时最大流速达到0.25 m/s，D为4/5（工况4）时最大流速则达到了0.30 m/s。这是由于浮堤淹没深度比D增大后，浮堤横向尺寸更大，使过水断面减小更为显著，因此断面流速增大。

图6 工况3 和4 浮堤附近断面流速分布Fig.6 Distributions of velocity around the dike in Cases 3 and 4

图7 为不同工况浮堤附近断面垂向紊动强度差异分布。图7 中的垂向紊动强度差异为正数，表示工况2 中相应位置的垂向紊动强度大于工况3 或4 中相同位置。对比图5 与7 可见：（1）工况3 和2 下，断面3#、4#中浮堤下游水体的垂向紊动强度较大，且工况3 表层水体的垂向紊动强度相较于工况2 更大，而中间部分水体的垂向紊动强度较小，这与工况3 回流更接近表层水体有关。工况4 浮堤下游水体的垂向紊动强度较工况2 和3 小；（2）工况3 和2 浮堤右侧（断面5#）水体的垂向紊动强度较小，而工况4 浮堤右侧部分区域垂向紊动强度较大，尤其是紧邻浮堤右侧的水体（图7(f)虚线框内），这是由于此处水体需绕过浮堤流向下游，在水平方向上发生了边界层分离，因此垂向紊动强度也较大。对比图7 和6 可知：（1）浮堤下游回流处垂向紊动强度相对其他区域均较大，当淹没深度比D达到4/5 时，虽然回流垂向尺度更大，但回流区域内流速绝对值较小，因而垂向紊动强度相对较小；（2）浮堤淹没深度比D增加使过水断面减小，故浮堤右侧水体流速较大，加之浮堤下方水体存在上升运动，且水流偏转角度较大，因此垂向紊动强度较大。

图7 工况2-3、工况2-4 浮堤附近断面垂向紊动强度差异分布Fig.7 Distributions of differences of vertical turbulent intensity around the dike between Cases 2-3 and Cases 2-4

图8 为工况3 和4 浮堤附近断面雷诺应力分布。从定性看，两种工况下的雷诺应力分布与工况2 规律相似，即回流和浮堤下方区域雷诺应力较大；分布规律与垂向紊动强度（图7）吻合。通过比较分析可知，浮堤淹没深度比D增大并未造成雷诺应力分布出现明显单调变化，淹没深度比为2/5 时，雷诺应力相较于淹没深度比1/5 时更大，但是淹没深度比为4/5 时的雷诺应力明显减小，故推测在淹没深度比为2/5 和4/5 之间出现极值。当淹没深度比D达到4/5 时，由于浮堤显著的阻碍作用，回流区内流动缓慢，同时过水断面显著减小，使浮堤右侧水体明显加速，故浮堤下游雷诺应力减小，其最大值也减小，而浮堤右侧水体雷诺应力明显增大。

图8 工况3 和4 浮堤附近断面雷诺应力分布Fig.8 Distributions of Reynolds stress around the dike in Cases 3 and 4

2.3 浮堤长度影响

图9 为工况5 和6 浮堤附近断面流速分布。与图4 对比可见：（1）工况2、5、6 中浮堤下游均产生回流，且沿水槽横向回流尺寸与浮堤长度大致相当，即浮堤长度比K（工况5 为1/6、工况2 为1/3、工况6 为1/2）越大，回流横向尺度越大；（2）工况2、5、6 中浮堤正下方、右侧水体流速均较大，这是由于过水断面减小使流速增加，且浮堤长度比K越大，过水断面减幅就越大，相应水体的流速增幅就越大，具体表现在工况5 浮堤正下方、右侧水体的最大流速达0.22 m/s，而工况2 和6 则达到0.25 m/s。

图9 工况5 和6 浮堤附近断面流速分布Fig.9 Distributions of velocity around the dike in Cases 5 and 6

图10 为工况5 和6 浮堤附近断面垂向紊动强度分布。由图10 可见：（1）浮堤下游回流所在区域垂向紊动强度较大；（2）对比图5，随着浮堤长度比K增大，工况5、2、6 浮堤下游回流的垂向紊动强度呈增大趋势，浮堤右侧水体的最大垂向紊动强度呈减小趋势。

图10 工况5 和6 浮堤附近断面垂向紊动强度分布Fig.10 Distributions of vertical turbulent intensity around the dike in Cases 5 and 6

图11 为工况5 和6 浮堤附近断面雷诺应力分布。由图11 可见：（1）浮堤下方和下游回流区域内水体雷诺应力较大，其余位置雷诺应力接近0；（2）对比图5，随着浮堤长度比K的增大，工况5、2、6 浮堤正下方雷诺应力最大绝对值比较接近，均大于0.5 Pa，无明显规律；而浮堤下游回流区雷诺应力达到1 Pa 的区域面积随浮堤长度比K的增加而增加。

图11 工况5 和6 浮堤附近断面雷诺应力分布Fig.11 Distributions of Reynold stress around the dike in Cases 5 and 6

2.4 浮堤厚度影响

图12 为工况7 和8 浮堤附近断面流速分布。对比图4 可见：随着浮堤厚度增加（工况7-工况2-工况8），浮堤下游的回流明显减弱，尤其是工况8 断面4#处已无回流，且浮堤下方水流平稳上升，流线顺畅延长至下游，这说明浮堤达到一定厚度后，浮堤对其下方水体具有较显著的稳流作用，对紊动有抑制作用；同时，浮堤下方水体流速减小，与浮堤下游回流减弱相符。图13 为工况7 和8 中断面6#、7#、8#处的流速分布，均可观察到浮堤右侧水体在15.0～17.5 cm 之间恢复均匀流状态，这表明浮堤对右侧水体的影响范围为15.0～17.5 cm。前文提到工况2 浮堤在水槽宽度方向的影响范围为20.0～22.5 cm，这说明浮堤的厚度与其在水槽宽度方向的影响范围之间不存在单调变化关系，可能存在一个临界厚度，在该厚度下浮堤沿水槽宽度方向的影响范围最大。同时，该临界厚度可能与水流条件、浮堤布置等有关，需进一步研究。

图12 工况7 和8 浮堤附近断面流速分布Fig.12 Distributions of velocity around the dike in Cases 7 and 8

图13 工况7 和8 浮堤右侧附近断面流速分布Fig.13 Distributions of velocity on the right side of the dike in Cases 7 and 8

图14 为工况2 与工况7、8 中浮堤附近断面垂向紊动强度差异分布，其中垂向紊动强度差异为正，表示工况2 相应位置的垂向紊动强度相较于工况7 或8 相同位置更大。对比图5 与14 可见：浮堤厚度减小时，浮堤下游回流区的垂向紊动强度呈增大趋势（图14(a)，14(b)）；浮堤厚度增加时，浮堤下游回流区的垂向紊动强度呈减弱趋势（图14(d)，14(e)）。这也体现了浮堤厚度越大，对水流紊动的抑制作用越显著。

图14 工况2-7、工况2-8 浮堤附近断面垂向紊动强度差异分布Fig.14 Distributions of differences of vertical turbulent intensity around the dike between Cases 2-7 and Cases 2-8

图15 为工况7 和8 浮堤附近断面雷诺应力分布。对比图5 可见：随着浮堤厚度增加，浮堤正下方和下游回流区水体的雷诺应力绝对值减小，这表明水流紊动剧烈程度降低，与图12 和14 展示的规律一致，即增加浮堤厚度会抑制水流紊动。

图15 工况7 和8 浮堤附近断面雷诺应力分布Fig.15 Distributions of Reynolds stress around the dike in Cases 7 and 8

3 结语

通过开展浮堤局部水流结构试验，探讨了浮堤及其尺寸对水流结构的影响，并得出如下结论：

（1）布置浮堤后，由于浮堤表面发生了边界层流动分离，表面水流流经浮堤后偏转向对岸，浮堤下游产生回流，且回流处的垂向紊动强度与雷诺应力绝对值较大；由于浮堤的挡水作用，过水断面面积减小，浮堤上游水流缓慢，而浮堤下方和右侧水流流速较大。

（2）浮堤淹没深度增大的影响：浮堤下游回流的垂向尺度相应增大，但淹没深度比D达到4/5 时垂向紊动强度明显减弱；同时，随着水槽过水断面减小，浮堤右侧水体流速明显增大，垂向紊动强度和雷诺应力绝对值也相应增大。

（3）浮堤长度增加的影响：浮堤下游回流的横向尺度相应增大，且垂向紊动强度和雷诺应力绝对值呈增大趋势；同时，随着水槽过水断面减小，浮堤右侧水体流速增大，垂向紊动强度减小。

（4）浮堤厚度增加的影响：浮堤下游回流明显减弱（当厚度为4 cm 时，断面4#处已无回流，且浮堤正下方水流平稳上升）；浮堤下方水体的流速呈减小趋势；浮堤在水槽宽度方向上的影响范围先增后减，说明存在临界厚度，在该厚度下浮堤在水槽宽度方向上的影响范围最大，有待进一步研究；同时，浮堤下游回流区水体的垂向紊动强度呈减弱趋势，浮堤正下方和下游回流区水体的雷诺应力绝对值减小。