留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

天然河道丁坝群局部冲刷三维数值模拟

戚蓝 曾庆达 吉顺文

引用本文:
Citation:

天然河道丁坝群局部冲刷三维数值模拟

    作者简介: 戚 蓝(1955—),女,山东威海人,教授,主要从事水工结构及流体分析研究。E-mail:lanqi0000@163.com.
    通讯作者: 曾庆达(E-mial:694542180@qq.com); 
  • 基金项目: 国家自然科学基金资助项目(51479178)
  • 中图分类号: TV135;TV863

Three-dimensional numerical simulation of local scour based on spur dikes in a natural river

    Corresponding author: ZENG Qingda, 694542180@qq.com;
  • 摘要: 为研究天然河道中丁坝长度对河床的冲刷影响,基于天然河道形状和实际丁坝群布置,提出适用于天然河道形状的丁坝群局部冲刷计算流体力学(CFD)三维数值模拟方法,建立了滦河马良子段丁坝群三维模型。采用RNG k-ε模型模拟计算,用VOF方法进行自由表面的追踪,并对丁坝局部冲刷坑形态及冲刷深度进行了数值模拟。对比实测数据,数值模拟结果与实测结果非常吻合。针对丁坝群中局部冲刷最严重的丁坝,对其长度进行调整,调整后冲坑深度明显改善。
  • 图 1  马良子段丁坝群布置

    Figure 1.  Layout of spur dike group along Maliangzi reach

    图 2  计算区域网格剖分

    Figure 2.  Meshes of computational domain

    图 3  坝头附近冲刷深度横纵断面对比

    Figure 3.  Comparison of transverse and longitudinal sections of the scouring depth near the spur dike head

    图 4  11#丁坝附近绕流流场及流速分布

    Figure 4.  Flow field and velocity distribution around No.11 spur dike

    图 5  丁坝群冲刷形态对比

    Figure 5.  Comparison of scour patterns around spur dikes

    图 6  测点布置

    Figure 6.  Layout of measuring points

    表 1  3种方案各测点数值对比

    Table 1.  Comparison of values of measuring points in three schemes

    测点 水面高度/m 流速/(m·s-1) 冲深/m 冲深减小值/m
    方案1 1 25.45 6.40 13.98
    2 26.20 6.00 7.00
    3 25.63 6.49 7.00
    4 25.30 6.50 10.70
    方案2 1 25.20 6.90 11.00 2.98
    2 26.00 5.60 6.50 0.50
    3 25.80 6.60 7.80 0.80
    4 25.90 6.50 9.50 1.20
    方案3 1 25.90 5.40 5.30 8.68
    2 26.00 3.50 5.00 2.00
    3 25.80 5.00 5.40 1.60
    4 25.80 4.00 4.70 6.00
    下载: 导出CSV
  • [1]

    TANG X L, DING X, CHEN Z C. Large eddy simulations of Three-dimensional flows around a spur dike[J]. Tsinghua Science & Technology, 2006, 11(1): 117-123.
    [2] 喻涛. 非恒定流条件下丁坝水力特性及冲刷机理研究[D]. 重庆: 重庆交通大学, 2013.

    YU Tao. Study on hydraulic characteristics and scour mechanism of spur dike in unsteady flow[D]. Chongqing: Chongqing Jiaotong University, 2013. (in Chinese)
    [3]

    BEY A, FARUQUE M A, BALACHANDAR R. Two-dimensional scour hole problem: role of fluid structures[J]. Journal of Hydraulic Engineering, 2007, 133(4): 414-430. doi: 10.1061/(ASCE)0733-9429(2007)133:4(414)
    [4] 吴学文, 詹义正. 非均匀河床上的丁坝局部冲刷问题[J]. 江西水利科技,2006,32(3):139-141. (WU Xuewen, ZHAN Yizheng. Discussion of the local erosion of the bankhead on the uneven river bed[J]. Jiangxi Hydraulic Science & Technology, 2006, 32(3): 139-141. (in Chinese) doi: 10.3969/j.issn.1004-4701.2006.03.005

    [5] 马永军. 丁坝回流尺度控制方法和减少回流淤积的实验研究[D]. 北京: 清华大学, 2003.

    MA Yongjun. Experimental study of control method of horizontal circulation and sedimentation behind a spur dike[D]. Beijing: Tsinghua University, 2003. (in Chinese)
    [6] 苏伟, 王平义, 喻涛, 等. 不同结构形式丁坝水毁过程分析[J]. 水运工程,2012(11):118-123. (SU Wei, WANG Pingyi, YU Tao, et al. Analysis of damage process of different types of spur dikes[J]. Port & Waterway Engineering, 2012(11): 118-123. (in Chinese) doi: 10.3969/j.issn.1002-4972.2012.11.023

    [7]

    TINGSANCHALI T, MAHESWARAN S. 2-D depth-averaged flow computation near Groyne[J]. Journal of Hydraulic Engineering, 1990, 116(1): 71-86. doi: 10.1061/(ASCE)0733-9429(1990)116:1(71)
    [8]

    MOLLS T, CHAUDHRY M H. Depth-averaged open-channel flow model[J]. Journal of Hydraulic Engineering, 1995, 121(6): 453-465. doi: 10.1061/(ASCE)0733-9429(1995)121:6(453)
    [9]

    JIA Y F, SCOTT S, XU Y C, et al. Three-dimensional numerical simulation and analysis of flows around a submerged rigid vegetation[J]. Advances in Water Resources, 2009, 32(4): 487-492. doi: 10.1016/j.advwatres.2008.11.014
    [10]

    DUAN J G, NANDA S K. Two-dimensional depth-averaged model simulation of suspended sediment concentration distribution in a groyne field[J]. Journal of Hydrology, 2006, 327(3-4): 426-437. doi: 10.1016/j.jhydrol.2005.11.055
    [11]

    NAGATA N, HOSODA T, NAKATO T, et al. Three-dimensional numerical model for flow and bed deformation around river hydraulic structures[J]. Journal of Hydraulic Engineering, 2005, 131(12): 1074-1087. doi: 10.1061/(ASCE)0733-9429(2005)131:12(1074)
    [12]

    OUILLON S, DARTUS D. Three-dimensional computation of flow around Groyne[J]. Journal of Hydraulic Engineering, 1997, 123(11): 962-970. doi: 10.1061/(ASCE)0733-9429(1997)123:11(962)
    [13] 杨兰, 李国栋, 李奇龙, 等. 丁坝群附近流场及局部冲刷的三维数值模拟[J]. 水动力学研究与进展(A辑),2016,31(3):372-378. (YANG Lan, LI Guodong, LI Qilong, et al. 3D numerical simulation of flow field and local scour around spur dikes[J]. Chinese Journal of Hydrodynamics (SerA), 2016, 31(3): 372-378. (in Chinese)

  • [1] 陈一梅张梦成 . 基于SVM的丁坝群束水攻沙功能预测. 水利水运工程学报, 2019, (3): 25-31. doi: 10.16198/j.cnki.1009-640X.2019.03.004
    [2] 程永舟吕行王文森黄筱云夏波 . 波流作用下淹没圆柱局部冲深影响因素分析. 水利水运工程学报, 2019, (6): 69-76. doi: 10.16198/j.cnki.1009-640X.2019.06.008
    [3] 潘存鸿汪求顺潘冬子 . 钱塘江涌潮对风场响应的三维数值研究. 水利水运工程学报, 2019, (6): 77-85. doi: 10.16198/j.cnki.1009-640X.2019.06.009
    [4] 陈一梅范丽婵李鑫 . 基于BP神经网络的丁坝坝头冲刷坑变化趋势预测. 水利水运工程学报, 2019, (6): 125-131. doi: 10.16198/j.cnki.1009-640X.2019.06.014
    [5] 刘明维曾丽琴陈刚沈立龙牟治忠王明镜 . 内河框架码头桩柱绕流水动力特性的二维数值模拟. 水利水运工程学报, 2020, (1): 74-83. doi: 10.12170/20181204005
    [6] 李健王平义许百强贺小含王梅力 . 黄河上游连续弯曲多滩险河道整治方案数值分析. 水利水运工程学报, 2020, (2): 1-6. doi: 10.12170/20190901001
    [7] 戚蓝肖厅厅张芝永黄君宝 . 涌潮水流CFD数值模拟. 水利水运工程学报, 2019, (3): 32-40. doi: 10.16198/j.cnki.1009-640X.2019.03.005
    [8] 陈述李雨晨许慧尚倩倩陈梓偲 . 长江下游河段桥墩压缩冲刷预测研究. 水利水运工程学报, 2019, (3): 16-24. doi: 10.16198/j.cnki.1009-640X.2019.03.003
    [9] 周舟曾诚周婕王玲玲丁少伟 . 等宽明渠交汇口流速分布特性数值模拟. 水利水运工程学报, 2020, (1): 32-39. doi: 10.12170/20190501005
    [10] 欧阳澍刘成林程永舟王文森王小明 . 透空潜坝附近水流特性的数值模拟. 水利水运工程学报, 2019, (1): 109-118. doi: 10.16198/j.cnki.1009-640X.2019.01.014
    [11] 刘建华杨中华李达李明 . 考虑底部过流的透水框架群阻力系数研究. 水利水运工程学报, 2020, (2): 1-6. doi: 10.12170/20190216001
  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  262
  • HTML全文浏览量:  162
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-20
  • 刊出日期:  2020-02-01

天然河道丁坝群局部冲刷三维数值模拟

    通讯作者: 曾庆达, 694542180@qq.com
    作者简介: 戚 蓝(1955—),女,山东威海人,教授,主要从事水工结构及流体分析研究。E-mail:lanqi0000@163.com
  • 1. 天津大学 水利工程仿真与安全国家重点实验室,天津 300072
  • 2. 浙江省水利河口研究院 工程安全研究中心,浙江 杭州 310020

摘要: 为研究天然河道中丁坝长度对河床的冲刷影响,基于天然河道形状和实际丁坝群布置,提出适用于天然河道形状的丁坝群局部冲刷计算流体力学(CFD)三维数值模拟方法,建立了滦河马良子段丁坝群三维模型。采用RNG k-ε模型模拟计算,用VOF方法进行自由表面的追踪,并对丁坝局部冲刷坑形态及冲刷深度进行了数值模拟。对比实测数据,数值模拟结果与实测结果非常吻合。针对丁坝群中局部冲刷最严重的丁坝,对其长度进行调整,调整后冲坑深度明显改善。

English Abstract

  • 20世纪50年代末,滦河开始兴建护岸丁坝工程。丁坝具有导流、护岸、防冲和稳定河势的作用,是保滩护岸工程中常见的水工建筑物[1]。但滦河上的丁坝群建设没有很好规划,丁坝长度、结构形式各异,而且大部分丁坝出现了不同程度的水毁。其中马良子段以险著称,若丁坝失效,滦河主流将直冲小埝,危及马良子村的安全。因此,模拟马良子丁坝群的局部冲刷,并基于丁坝冲刷规律对其优化十分必要。

    目前,丁坝的局部冲刷是众多专家学者研究的重点,丁坝建成后由于丁坝上游壅水,会在上游坝根出现一小的立轴回流;至坝头附近,出现下潜流,产生立轴漩涡并相互碰撞、破碎、向下游运移导致坝头局部冲刷剧烈。丁坝下游则出现较大的立轴回流及由此于下游坝根处引发次生反向小回流。这些部位会产生较大的切应力,导致河床发生局部冲刷。喻涛[2]测量了丁坝横向冲坑的发展过程。Bey等[3]认为丁坝冲坑的产生是由于丁坝的存在使得丁坝附近湍流对河床底部有冲蚀作用,并认为最大冲深和冲坑形状同样值得研究,因为冲坑形状一定程度上也会影响丁坝水毁失稳。吴学文等[4]认为丁坝局部冲刷主要是坝头的绕流作用,由此建立了坝头绕流挤压流动物理图示,考虑了非均匀沙的起动问题,得出了非均匀沙河床上丁坝局部冲深公式。事实上,天然河道中的丁坝三维水流情况十分复杂,许多参数难以计量,马永军[5]分析以往研究结果发现,各专家学者们都只是把回流长度、水流边界、冲深直接联系起来,并没有找出各参数之间的关系。苏伟等[6]对不同结构形式丁坝进行物理模型试验,发现丁坝的结构形式对水毁影响明显。计算机技术的发展,为复杂水流问题的求解提供了计算手段,Tingsanchahi等[7]利用二维平均水深模型并引入校正因子改善k-ε模型,以求解河床底部应力。Molls等[8]通过结合恒定的涡黏性湍流模型改进了不稳定二维平均水深模型。Jia等[9]利用二维平均水深模型计算丁坝附近水流流线和流速。Duan等[10]用改进二维模型对天然河道中泥沙输运过程进行仿真。结果表明,二维模型只能近似计算泥沙水平运动,无法较好地模拟泥沙在河道中的运输。越来越多的学者发现丁坝附近水流表现出很复杂的三维特性[11],因此,开始丁坝的三维数值模拟研究。Nagata等[11]利用三维水利模型计算了雷诺平均方程和斯托克斯方程,成功模拟了丁坝周围的泥沙输运过程。Ouillon等[12]利用三维k-ε模型研究了丁坝周围三维水流结构和自由液面形状,求解了坝后冲刷区水流的水力特性。杨兰等[13]通过使用Flow-3D软件,对上挑丁坝群附近流场及局部冲刷进行了三维数值模拟,并得出了可适用于计算丁坝群的湍流模型和推移质输沙率模型。但目前对丁坝群的三维数值模拟或物理模型模拟大部分针对的是直槽式水槽,而非实际工程中复杂的天然河道。本文首先基于不同丁坝结构形式的三维水槽模型模拟,对滦河上各种丁坝的冲刷规律进行探索,作为调整丁坝的依据。然后对天然河道的丁坝群进行三维数值模拟,对比分析冲刷数据与实测地形,其结果可为工程实际提供一定的参考。

    • 丁坝群流场三维流动控制方程为连续性方程和动量方程:

      $ \partial {u_i}/\partial {x_i} = 0 $

      (1)

      $ \partial {u_i}/\partial t + {u_j}\partial {u_i}/\partial {x_j} = - \frac{1}{\rho }\frac{{\partial P}}{{\partial {x_i}}} + \frac{\partial }{{\partial {x_j}}}\left( {2v{S_{ij}} - \overline {u' {_i}{u' _j}} } \right) $

      (2)

      $ {S_{ij}} = \frac{1}{2}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right) $

      (3)

      $ \overline {u' {_i}{u' _j}} = {v_{\rm{t}}}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right) - \frac{2}{3}k{\delta _{ij}} $

      (4)

      式中:uii方向的分速度;uii方向的脉动速度;P为压力;Sij为应变率张量;$\overline {u' {_i}{u' {_j}}}$为雷诺应力张量;ρ为流体密度;v为动力黏度;vt为湍流黏度;k为湍动能;δij为克罗内克符号(δij=1,i=jδij=0,ij)。

    • 标准k-ε模型用于强旋转流或带有弯曲壁面流动时容易出现失真。RNG k-ε模型通过修正湍动黏性系数和在ε方程中增加了反映主流的时均应变率[7],可以很好地模拟高应变率和流线弯曲程度较大的流动,而天然河道的复杂地形以及丁坝与水流之间的相互作用会带来剧烈的水流变形和破碎,且许多专家学者认为RNG k-ε模型对丁坝得到的结论更为可靠。

      RNG k-ε控制方程:

      $ \frac{{\partial k}}{{\partial t}}{\rm{ }} + {u_j}\frac{{\partial k}}{{\partial {x_i}}} = {\tau _{ij}}\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{\partial }{{\partial {x_j}}}\left[ {\frac{1}{\rho }\left( {\mu + \frac{{{\mu _{\rm{t}}}}}{{{\sigma _k}}}} \right)\frac{{\partial k}}{{\partial {x_j}}}} \right] - \varepsilon $

      (5)

      $ \frac{{\partial \varepsilon }}{{\partial t}} + {u_j}\frac{{\partial \varepsilon }}{{\partial {x_i}}} = {C_{1\varepsilon }}\frac{\varepsilon }{k}{\tau _{ij}}\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{\partial }{{\partial {x_j}}}\left[ {\frac{1}{\rho }\left( {\mu + \frac{{{\mu _{\rm{t}}}}}{{{\sigma _k}}}} \right)\frac{{\partial \varepsilon }}{{\partial {x_j}}}} \right] - {C_{2\varepsilon }}\rho \frac{{{\varepsilon ^2}}}{k} $

      (6)

      式中:k为湍动能,$k = \dfrac{1}{2}\overline {{u' _i}{u' _j}} $ε为湍动能耗散,$\varepsilon = v\overline {\left( {\dfrac{{\partial {u_t}' }}{{\partial {x_k}}}} \right)\left( {\dfrac{{\partial {u_t}' }}{{\partial {x_k}}}} \right)} $μt为涡黏性系数,μt=Cμρk2/ετij为雷诺应力,${\tau _{ij}} = \dfrac{{{\mu _{\rm{t}}}}}{\rho }\left( {\dfrac{{\partial {u_i}}}{{\partial {x_j}}} + \dfrac{{\partial {u_j}}}{{\partial {x_i}}}} \right) - \dfrac{2}{3}\rho k{\delta _{ij}}$;模型常数Cμ=0.085,C1ε=1.42,C2ε=1.68,σk=σε=0.717 9。

    • 目前,研究丁坝局部冲刷的规律主要有5个理论途径,分别是流速、切应力、能量平衡、统计法则及沙波运动,而湍流流场和泥沙输运对于丁坝群局部冲刷坑的发育影响在理论层面还不是十分成熟。基于前人对丁坝冲刷规律的试验以及泥沙运动理论分析,选用以希尔兹(Shields)数为基础的泥沙推移质输沙率模型作为此次丁坝群冲刷三维数值模拟的理论模型。

      与杨兰等[13]的泥沙输运模型相同,模型计算中不考虑悬移质对冲刷坑发育的影响,只分析河床底部推移质的冲刷作用,并认为推移质为均匀的泥沙颗粒。采用的推移质输沙率公式以切应力为主要参数,且切应力与输沙率的关系成正比。

      水流切应力:

      $ \tau = \rho U_*^2 $

      (7)

      泥沙起动的临界切应力:

      $ {\tau _{\rm{c}}} = \rho {U}_{*{\rm{c}}}^2$

      (8)

      希尔兹数:

      $\theta = \tau /(g\left( {{\rho _{\rm{s}}} - \rho } \right)d) $

      (9)

      临界希尔兹数:

      $ {\theta _{\rm{c}}} = {\tau _{\rm{c}}}/(g\left( {{\rho _{\rm{s}}} - \rho } \right)d) $

      (10)

      式中:U*为摩阻流速;U*c为临界摩阻流速;ρρs分别为水流和泥沙密度;d为泥沙平均粒径。

      假定单宽推移质输沙率计算式为:

      ${q_{\rm{b}}} = {\rho _{\rm{s}}}\left( {{\text{π}} /6} \right){d^3}{u_{\rm{b}}}p/{d^2} $

      (11)

      式中:qb为单宽推移质输沙率;ub为推移质的平均输运速度;p为泥沙起动概率。对于实际工程地质情况的土工试验,可以根据上式得到ubp,即可求解qb

      $ {F_{\rm{D}}} = {C_{\rm{D}}}\left( {{\text{π}}/4} \right){d^2}\rho {\left( {a{U_*} - {u_{\rm{b}}}} \right)^2}/2 $

      (12)

      动摩擦力为:

      $ {f_{\rm{D}}} = g\left( {{\rho _{\rm{s}}} - \rho } \right){\text{π}} {d^3}\beta /6 $

      (13)

      FD=fD并联立(7),得到:

      $\frac{{{u_{\rm{b}}}}}{{{U_*}}} = a\left( {1 - \sqrt {\frac{{4\beta }}{{3{a^2}{C_{\rm{D}}}\theta }}} } \right) $

      (14)

      $ {\text{令}}{\theta _0} = 4\beta /\left( {3{a^2}{C_{\rm{D}}}} \right),{\text{则}}{u_{\rm{b}}}/{U_*} = a\left( {1 - \sqrt {{\theta _0}/\theta } } \right) $

      (15)

      式中:aU*为推移质运动的水流速度,当离沙床较近时,a=6~10;CD为推力系数。当θ0=0时,ub=0,θ0相当于止动相对切应力,应小于临界相对切应力θc。由此,式(15)可写为:

      $ {u_{\rm{b}}}/{U_*} = a\left( {1 - 0.7\sqrt {{\theta _0}/\theta } } \right) $

      (16)

      取泥沙临界切应力和运动颗粒所受切应力之和为切应力,即:

      $ \tau = {\tau _{\rm{c}}} + ng\left( {{\rho _{\rm{s}}} - \rho } \right){\text{π}} {d^3}\beta /6 $

      (17)

      式中:n为泥沙颗粒总数,而单位面积河床泥沙颗粒总数1/d2p的关系为:p=n/d-2。联立式(9)和(10),得到床面泥沙颗粒的起动概率为

      $ p = 6\left( {\theta - {\theta _c}} \right)/({\text{π}} \beta) $

      (18)

      综上所述,影响单宽推移质输沙率的两个重要因素为ubp,将ubp的表达式代入(11),并取动摩擦系数β=0.08,常数a=9.3。可得如下推移质输沙公式:

      $ \left\{ \begin{array}{l} {u_{\rm{b}}} = 11.6{\rho _{\rm{s}}}d\left( {\theta - {\theta _{\rm{c}}}} \right)\left( {\sqrt \theta - 0.7{{\sqrt \theta }_{\rm{c}}}} \right) \times \sqrt {gd\left( {{\rho _{\rm{s}}} - \rho } \right)/\rho } {\rm{ }},\theta \geqslant {\theta _{\rm{c}}}\\ {q_{\rm{b}}} = 0,\theta < {\theta _{\rm{c}}} \end{array} \right. $

      (19)

      式(19)即为所采用的泥沙运动模型。可以发现,对于某些泥沙颗粒,推移质输沙率只与希尔兹数有关,只有当希尔兹数大于临界希尔兹数时,泥沙才会起动被冲刷,采用θc=0.05。

    • 实际工程丁坝段如图1所示,可见11#丁坝相比周围的丁坝明显较长,不够合理。这是由于丁坝有束窄河床的作用,河床被束窄后,河道断面面积减小;虽然丁坝前有一定程度的壅水,但断面面积仍有大比例收缩,导致相同流量下流速迅速增大,特别是河道底部的切应力会迅速增大,下切河床造成冲刷。同时,该河段自9#丁坝后河道宽度明显缩窄,流速更大,由于滦河段河床地质条件差,基本是粉细砂组成,更容易遭到冲刷。

      图  1  马良子段丁坝群布置

      Figure 1.  Layout of spur dike group along Maliangzi reach

      对比河底实测地形图发现,11#丁坝后有很深的冲刷坑,最深处超过13 m,因此选择11#丁坝作为主要的研究对象,研究不同丁坝长度对丁坝群的冲刷影响。

      计算模型长1 400 m,宽1 000 m,高30 m,其中11#丁坝上游400 m至下游300 m为动床计算区域,其余为定床。河床主槽糙率为0.02,平均水深为5 m,为非淹没丁坝群,计算丁坝群位于左岸,计算河道模型平面轮廓与天然河道一致,考虑到河道宽度远远大于水深,因此河道断面形状取矩形。原始11#丁坝长为34.5 m,动床所铺泥沙厚度为20 m,泥沙平均粒径d=0.018 mm,密度ρs=2 650 kg/m3。计算区域采用矩形网格划分为ABC区,A区和C区为定床,B区为动床。对动床部分加密网格,总网格数量约133万,网格分块划分见图2

      图  2  计算区域网格剖分

      Figure 2.  Meshes of computational domain

      计算时采用有限体积法离散控制方程,湍流模型为RNG k-ε模型。上游设为速度进口,根据上游实测流速和平均深度,给定模型流量为8 430 m3/s,平均水深5 m,出口设为自由出流,上表面采用VOF法捕捉液体自由面,固体边界为无滑移面,近壁面采用标准壁面函数进行处理。对比11#坝头附近的冲刷深度与实测地形图,实测与模拟冲刷地形拟合较好,如图3所示。可见,在复杂天然河道形状下,可以利用RNG k-ε模型和推移质输沙率模型模拟丁坝的局部冲刷。

      图  3  坝头附近冲刷深度横纵断面对比

      Figure 3.  Comparison of transverse and longitudinal sections of the scouring depth near the spur dike head

    • 滦河上马良子段丁坝结构形式简单,多为下挑式丁坝,挑角约70°~85°,变化幅度不大,且坝头形式单一,为直立式丁坝。通过丁坝设计资料以及地形勘察资料发现,马良子段11#丁坝设计明显过长,设计长度34.5 m,比相邻丁坝长约10 m,而且10#~12#丁坝在河流收缩段。测量数据也表明11#丁坝坝后是冲刷最严重的区域,最大冲深达13.98 m。因此丁坝长度不合理很可能是造成严重冲刷的原因。本文将原坝长方案定为方案1,将11#丁坝缩短5 m和10 m,对应方案2和3,对丁坝周围水力特性及冲刷特性进行研究。

    • 方案1~3计算域水流流场与流速如图4所示。可见,随着11#丁坝坝长的减小,水流随之平顺,丁坝后端坝田附近形成的漩涡减小。同时,在原有丁坝群间距不变的情况下,缩短11#丁坝依然可以保持丁坝群对水流的有效影响,即坝前水流流入坝后坝田区与坝田内水流相互作用,与岸边不直接接触,主河槽流速较两岸大,方案1中过11#坝后主河槽流速明显增加,方案2和3增速不明显。

      图  4  11#丁坝附近绕流流场及流速分布

      Figure 4.  Flow field and velocity distribution around No.11 spur dike

    • 丁坝群冲刷形态对比如图5所示,可以发现调整方案1和2中丁坝群最大冲深坑出现在11#丁坝坝头附近,并且在主流区形成一狭长冲刷带;而方案3丁坝群则没有明显冲刷坑,冲刷深度较为一致。这说明缩短丁坝群中过长的丁坝可以有效减轻其坝头的局部冲刷,丁坝群调整后河床冲刷明显减轻。

      图  5  丁坝群冲刷形态对比

      Figure 5.  Comparison of scour patterns around spur dikes

      根据实测冲刷地形在11#丁坝坝头附近设置测点,其中测点1是原坝后冲坑最深处,测点2和4分别是测点1上、下游100 m位置,测点3为测点1的右边5 m处,测点布置如图6所示。丁坝群冲刷数值对比如表1所示。可见,缩短过长的丁坝长度可以有效减小丁坝局部冲刷深度,从而降低丁坝的水毁风险。

      表 1  3种方案各测点数值对比

      Table 1.  Comparison of values of measuring points in three schemes

      测点 水面高度/m 流速/(m·s-1) 冲深/m 冲深减小值/m
      方案1 1 25.45 6.40 13.98
      2 26.20 6.00 7.00
      3 25.63 6.49 7.00
      4 25.30 6.50 10.70
      方案2 1 25.20 6.90 11.00 2.98
      2 26.00 5.60 6.50 0.50
      3 25.80 6.60 7.80 0.80
      4 25.90 6.50 9.50 1.20
      方案3 1 25.90 5.40 5.30 8.68
      2 26.00 3.50 5.00 2.00
      3 25.80 5.00 5.40 1.60
      4 25.80 4.00 4.70 6.00

      图  6  测点布置

      Figure 6.  Layout of measuring points

    • 基于天然河道形状,利用推移质输沙率模型和VOF自由表面处理法,对非淹没丁坝群的绕流流场、流速、局部冲刷进行三维数值模拟研究,并与实测资料进行对比,得到以下主要结论:丁坝群中某一丁坝长度较周围丁坝明显偏长易引起河道水流集中,坝后漩涡加剧,且丁坝长度愈长,阻水能力愈大,坝头冲刷越剧烈。最大冲刷发生在坝头处,且在冲坑内略偏坝轴线下游,上游部分坡度较陡,范围较小,下游部分坡度较缓,范围较大。说明丁坝群设计时应合理规划丁坝长度,丁坝群内各个丁坝长度不协调可能会形成新的水害。

      由于实际工程中丁坝的冲刷坑形成非常复杂,本次模拟的冲刷坑下游冲刷与实际还有一些偏差,对丁坝局部冲刷机理和实际水流过程有待进一步研究。

参考文献 (13)

目录

    /

    返回文章
    返回