主河道岸坡角和支流汇入量对泥沙冲淤的影响

付梁其; 李勋章; 何晓静; 胡腾腾; 潘磊; 刘飞诗; Liang-qi FU; Xun-zhang LI; Xiao-jing HE; Teng-teng HU; Lei PAN; Fei-shi LIU

中国农村水利水电 ›› 2022 ›› (1) : 28-31,38.

水环境与水生态

主河道岸坡角和支流汇入量对泥沙冲淤的影响

作者信息 +

The Influence of the Bank Slope Angle and Tributary Inflow on Sediment Erosion and Deposition

Author information +

稿件信息 +

摘要

为了分析常州城市河道淤积的机制，采用了量纲分析和水槽实验的方法，研究支流河道汇入主流河道过程中，主河道岸坡角、支流与主流的流量比、密度弗劳德数对汇流区冲刷的影响。研究将支流水槽与主流水槽的汇流角定为90°，采用不同的流量比和密度弗劳德数对四种不同岸坡角（45°、60°、75°、90°）进行对比实验。结果表明：支流河道中的流量与主流河道的流量比值越小时，岸坡角对冲刷深度的影响越大；对于任何给定的密度弗劳德数，越小的岸坡角可导致越小冲刷深度。依据拟合得到的两个经验关系式，岸坡角对泥沙淤积高度的影响大于对冲刷深度的影响；密度弗劳德数对无量纲冲刷深度影响最大，岸坡角对无量纲淤积高度的影响最大。

Abstract

In order to analyze the mechanism of the siltation of urban river in Changzhou， this paper uses the method of dimensional analysis and water channel experiment to study the influence of bank slope angle， flow ratio and density Froude number on the scour of the confluence area during the process of tributary river flowing into the mainstream river. In this study， the confluence angle of tributary flume and mainstream flume is set at 90°， and four different bank slope angles （45°， 60°， 75°， 90°） are compared with different flow ratios and density Froude numbers. The results show that the smaller the ratio of flow in tributary to the flow of the main stream， the greater the influence of the bank slope angle on scour depth is. For any given density Froude number， the smaller bank slope angle can lead to the smaller scour depth. According to the two empirical relations obtained by fitting， the influence of the bank slope angle on the height of sediment deposition is more than that on the scour depth， and the density Froude number has the greatest influence on dimensionless scour depth， and the bank slope angle has the greatest influence on the dimensionless sedimentation height.

关键词

支流 / 岸坡 / 密度弗劳德数 / 冲刷深度 / 淤积高度

Key words

tributary / bank slope / density Froude number / scour depth / siltation height

基金

水利部水科学与水工程重点实验室开放研究基金项目(YK914012)

水利部公益性行业科研专项经费项目(201401083-2)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

付梁其 , 李勋章 , 何晓静 , 胡腾腾 , 潘磊 , 刘飞诗. 主河道岸坡角和支流汇入量对泥沙冲淤的影响[J].中国农村水利水电, 2022(1): 28-31,38

Liang-qi FU , Xun-zhang LI , Xiao-jing HE , Teng-teng HU , Lei PAN , Fei-shi LIU. The Influence of the Bank Slope Angle and Tributary Inflow on Sediment Erosion and Deposition[J].China Rural Water and Hydropower, 2022(1): 28-31,38

0 引言

水流对河道的冲刷会产生泥沙淤积，从而造成水工建筑物运行受阻，成为了城市河道修筑、疏浚过程中必须解决的技术问题。在支流与主河道汇合的汇流处，循环涡流的形成和流速的增加会导致汇流后的河道冲刷和泥沙沉积^［1］。汇流处的流型很复杂，受许多因素的影响，如：支流与主流的比率（以下简称流量比）、主流与支流河道之间的动量传递速率等^［2，3］。水流汇合处可分为六个水动力区域：①合流上游拐角处的停滞区；②合流点处的流动偏转区；③合流下游拐角附近的分离区；④主河道中出现最大流速区（分离区导致的水流收缩）；⑤汇合后河道最大流速区下游的水流恢复区；⑥两个合流水流之间形成的剪切层。据前人研究^［4，5］，汇流处的河床剪应力增强是由于汇合点处主流和支流的组合，以及强涡产生的螺旋运动和动能。在文献［4，5］中，关于汇流区的水流流态和泥沙输移的研究已有不少报道，研究结果表明，将汇流角从15°增加到90°可使得底流速增加1.3倍，较大的汇流角会显著增加最大冲刷深度。另外，当采用较大的汇流角和流量比时，泥沙淤积的高度增加。

河流生态保护和泥沙淤积治理已上升为重大的国家战略^［6］，揭示河流泥沙淤积的机制是各地生态经济区发展的基础。对于泥沙淤积研究，李泽刚^［7］分析了黄河入海泥沙扩散、淤积的形成与演变等规律，提出了河口段河道整治和拦门沙治理等稳定黄河河口流路的对策。王恺忱^［8］分析了黄河下游及河口河段冲淤特性，进行了黄河河口延伸改道对下游河道的影响研究，用实测资料证实了河口尾闾入海流路的绝对长度是河口河段和下游河道冲淤幅度与发展趋势的制约因素。胡春宏等^［9，10］分析了水沙运动与演变基本规律，提出了黄河口治理的方向与措施。余欣等^［11］通过试验研究了尾闾河道出汊机制，提出了出汊阈值、触发条件以及河口流路稳定综合判别指标。王崇浩等^［12-14］利用一、二维连接数学模型研究了河口潮流与泥沙输移过程，对河口演变对下游河道反馈的影响进行了分析。上述泥沙淤积的研究成果有效地指导了各类河道的疏浚和整治^{［15，16］}。由此，本文以常州市金坛区城市河道“畅流活水”工程为背景，研究在支流河道汇入主流河道过程中，主河道岸坡角、支流流量比和弗劳德数对汇流区冲刷的影响。

1 研究方法

1.1 理论分析

Ghobadian和Shafai-Bajestan对考虑不同汇流角（60°、75°和90°）的冲刷模式进行了试验研究^［17］，并提出了汇流后最大冲刷深度（

d_{s e}

）和泥沙淤积高度（

H_{s e}

）的近似计算方法，如下所示：

\frac{d_{s e}}{y_{3}} = 0.029 Q_{r}^{0.73} F r_{g 3}^{2.03} α^{0.35}

（1）

\frac{H_{s e}}{y_{3}} = 0.27 Q_{r}^{0.638} - 0.559 F r_{g 3}^{0.321} + 0.21 l n (α)

（2）

式中：

d_{s e}

为最大冲刷深度；

Q_{r}

为支流与主流的流量比；

y_{3}

为汇流后河道的水深；

H_{s e}

为汇流后的泥沙淤积高度；

α

为汇流角；

F r_{g 3}

为汇流后河道的密度弗劳德数。

F r_{g 3}

的表达式如下：

F r_{g 3} = \frac{V_{3}}{\sqrt[]{g (G_{s} - 1) d_{50}}}

（3）

式中：

V_{3}

为汇合后河道的流速；

G_{s}

为泥沙颗粒的比重（本研究中取2.63）；

d_{50}

为泥沙粒径中值。

基于以上理论，本文开展了量纲分析，并尽可能地全面考虑可能对汇流处冲刷有影响的不同变量。量纲关系式如下所示：

\begin{array}{l} f (d_{s e}, H_{s e}, d_{50}, ρ_{s}, Q_{b}, B_{1}, B_{2}, B_{3}, Q_{1}, Q_{2}, Q_{3}, \\ y_{1}, y_{2}, y_{3}, ρ, g, σ, μ, α, θ, S_{0}) = 0 \end{array}

（4）

式中：

ρ_{s}

为泥沙颗粒密度；

Q_{b}

为推移质；

B

为渠道的底宽；

Q

为渠道的流量；

ρ

为水体密度；

g

为重力加速度；

σ

为表面张力；

μ

为动力黏度；

θ

为主渠道岸坡角；

S_{0}

为汇流处主渠道纵向河床坡度。下标1、2和3分别指上游、支流和汇合后河道。

利用BuckinghamΠ定理对汇流处冲刷的量纲分析结果如下式所示：

\frac{d_{s e}}{y_{3}}, \frac{H_{s e}}{y_{3}} = f (\frac{Q_{2}}{Q_{3}}, \frac{Q_{b}}{Q_{3}}, \frac{B_{2}}{B_{3}}, α, θ, S_{0}, R_{e 3}, W_{e 3}, F r_{g 3})

（5）

式中：Re为雷诺数；We为韦伯数。

由于实验中黏度和表面张力的影响可以忽略不计，因此本研究中没有考虑雷诺数和韦伯数。此外，汇流处上游不允许推移质运移，因此

Q_{b}

可以删除。根据Gurram^［18］研究结果，在亚临界流条件下，河床坡度对汇流处的流型影响可以忽略不计。由于所有实验都是在亚临界条件下进行的，因此本研究不考虑

S_{0}

。实验采用恒定的汇合角

α

，即90度汇合，则式（5）可简化为式（6）：

\frac{d_{s e}}{y_{3}}, \frac{H_{s e}}{y_{3}} = f (\frac{Q_{2}}{Q_{3}}, θ, F r_{g 3})

（6）

1.2 实验设置

在理论基础上，本文开展了水槽试验研究，共进行了64次试验。图1中，主水槽长6 m，宽50 cm，高50 cm；支流水槽长5 m，宽24 cm，高50 cm；支流水槽与主水槽垂直连接，即汇合角为90度。两个水槽的纵向河床坡度设为零，木片制成的底板安装在水槽底部上方14 cm处，底板之间的空隙由沙子填充，然后在水槽中加载沙床。为了研究岸坡对汇流处冲刷的影响，本研究在主水槽中采用了45°、60°、75°和90°的不同倾角斜墙。图2展示了实验室中泥沙的冲刷与淤积。主水槽和支水槽都安装了供水机泵与闸阀，每台机泵的供水流量都经过试验，能够灵活控制流量，故而主流和支流的流量可以准确分配。

图1 实验装置示意图

Fig.1 Schematic diagram of experimental device

Full size|PPT slide

图2 泥沙的冲刷与淤积（以45°岸坡为例）

Fig.2 Erosion and deposition of sediment （45°bank slope）

Full size|PPT slide

实验采用的总排水量均为25 L/s，具体流量由位于主水槽下游端的校准矩形堰来测量。本研究中汇合处连接的3个通道的流量值见表1。为了调节流量，首先在主渠道调节流量，同时关闭支流渠道出口闸门。然后，位于连接点的支流河道闸门逐渐打开，以建立支流流量。在主水槽下游端使用大量叠梁来调整尾水深度，调整范围为0.152~0.235 m，从而将密度弗劳德数调整在1.28~2.58的理想范围内。此外，还使用了一个精度为 ±0.1 mm的超声波测量仪来测量水面高程。

表1 汇流连接处的三个通道流量值

Tab.1 Three channel flow values at confluence connection

Q ₁/（m³·s^-1）	Q ₂/（m³·s^-1）	Q ₃/（m³·s^-1）	Q ₂/Q ₃
0.020 15	0.004 85	0.025 0	0.194
0.018 1	0.006 9	0.025 0	0.276
0.015 0	0.010 0	0.025 0	0.400
0.011 2	0.013 8	0.025 0	0.552

在实验中，均匀沙被用作沉积物材料，d ₅₀=1.28 mm，标准偏差为1.32 mm。根据式（6）所提出的3个无量纲参数，本研究在实验中采用了4个岸坡角和4个流量比，以评估流量比、汇合后水流条件和河岸坡度对汇合处冲刷的影响。此外，在实验室条件下，通过控制水流深度，对每一流量比采用四个弗劳德数和密度弗劳德数。每次试验结束时，在冲刷坑内建立平衡条件。在初步试验中，θ=45°时获得了最长的平衡时间。因此，在θ=45°的条件下，在72 h内进行了一系列实验，以确保达到平衡条件。对于最大和最小流量比，在不到6 h的时间内获得了几乎90%的最大冲刷深度。然后，在不干扰泥沙床的情况下小心地将水槽中的水排干；随后，使用分级杆测量冲刷坑的尺寸和泥沙淤积的高度。在图2中，在45°岸坡角的实验室水槽中显示了泥沙冲刷和淤积的形成。最后，利用SPSS统计软件进行多元回归分析，得到最大冲刷深度与泥沙淤积高度在上述参数下的无量纲经验关系。

2 结果与讨论

在每个实验中，汇流处观察到一对漩涡并向相反方向移动，从而引发河床物质的冲刷。同时，由于汇合处的水流收缩产生的高速导致了河床物质的输送，并在合流后河道中形成泥沙淤积。最后，随着水流断面的恢复（即增大），流速降低。

2.1 最大冲刷深度

图3为无量纲最大冲刷深度（

d_{s e} / y_{3}

）随密度弗劳德数的变化情况。如图3所示，虽然在所有流量比下，

d_{s e} / y_{3}

随Fr _g3增加而增大，但在最大流量比下，

d_{s e} / y_{3}

的增加率最高，见图3（d）。这一结果与Ghobadian和Shafai-bajestan^［17］提出的矩形横截面主河道的实验结果一致。对比图3（a）~（d）可以看出，对于给定的密度弗劳德数（Fr _g3）和岸坡角（θ），流量比的增加导致冲刷坑更深，即产生更大的

d_{s e} / y_{3}

。该结果说明了流量比的增加在汇流处分离区扩张中的作用，更大的流量比可导致更高的流速和河床剪应力，从而导致更深的冲刷坑。另一方面，最小的岸坡角（θ=45°）对应的

d_{s e} / y_{3}

最小，这可归因于岸坡角越小时，支流进入主河道的水流越缓，使得水流与河床物质的接触较小，冲刷最小。最后，岸坡角对冲刷深度的影响随Q_r 的增大而显著减小。换言之，当在汇流处建立适当的流量比时，河流的岸坡角才会影响冲刷深度。因此，本文认为河岸坡度对汇流处冲刷深度的影响在实践中不能忽略。通过多元回归方法（R ²=0.83）得出公式（7）以近似计算无量纲最大冲刷深度，公式（7）如下：

\frac{d_{s e}}{y_{3}} = 0.121 Q_{r}^{1.104} F r_{g 3}^{2.208} {(s i n θ)}^{0.451}

（7）

图3 无量纲最大冲刷深度随密度弗劳德数的变化

Fig.3 Variation of dimensionless maximum scour depth with density Froude number

Full size|PPT slide

虽然岸坡角（θ）对冲刷深度有影响，但通过比较式（7）右侧3个变量的指数，可以得出θ对最大冲刷深度的影响小于密度弗劳德数（Fr _g3）和流量比（Q_r ）的影响。此外，图3表明，对于较低的密度弗劳德数，冲刷深度被略微低估，而随着密度弗劳德数和岸坡角的增大，误差增大。

本研究认为，Ghobadian和Shafai-Bajesta的经验公式忽略了不同岸坡的影响，仅考虑了恒定岸坡（θ=90°）下不同汇流角的作用。而本研究的无量纲最大冲刷深度的散点图能够比较和验证前人研究的成果和假设，同时，本研究所提计算公式具有考虑河岸坡度的优点。

2.2 泥沙淤积最大高度

在每个实验的过程中，水槽内形成了一个泥沙淤积坝并逐渐向下游移动。当支流槽与主流槽的流量比增大时，水流变得平缓，而流量比和密度弗劳德数的减小则导致了较尖锐的淤积高度。图4考虑了不同的流量比和岸坡角，并根据密度弗劳德数绘制了汇流处淤积的最大无量纲高度（

H_{s e} / y_{3}

）变化情况。如图所示，当流量比增大到Q_r =0.552时，θ=75°的实验结果与θ=90°的实验结果接近。此外，随着θ的增加，当Q_r =0.194、0.276、0.400和0.552时，

H_{s e} / y_{3}

的平均增长率分别为41.7%、29.4%、27.8%和23.4%。通过多元回归方法（R ²=0.80）得出公式（8）以近似计算无量纲最大淤积高度，公式（8）如下：

\frac{H_{s e}}{y_{3}} = 0.296 Q_{r}^{0.359} F r_{g 3}^{- 1.046} {(s i n θ)}^{1.782}

（8）

图4 不同岸坡角和流量比下最大无量纲淤积高度与密度弗劳德数的关系

Fig.4 Relationship between maximum dimensionless deposition height and density Froude number under different bank slope angle and discharge ratio

Full size|PPT slide

根据式（8）右侧的指数可以得出，岸坡角对淤积最大无量纲高度的主要影响是隐含的，这一点与岸坡角对无量纲最大冲刷深度的影响不同。值得注意的是，较大的流量比（Q_r ）和岸坡角（θ）会导致较高的泥沙淤积高度，而较大的密度弗劳德数（Fr _g3）会降低河道泥沙淤积的最大高度。图4表明，密度弗劳德数越高，误差越小。此外，岸坡角的增加导致了淤积高度估计误差的增大。应用相同的假设，公式（8）得出的结果与Ghobadian和Shafai-Bajestan的结果非常接近。随着Fr _g3的增加，泥沙淤积的最大高度一般出现在下游方向离合流口较远的地方。根据不同的Q_r 和θ值，本文可得对应4个Fr _g3值的4个最大淤积高度的位置。结果显示，在岸坡倾斜角度较小（θ=45°和60°）的情况下，在倾斜的河岸上近处产生了泥沙淤积。同时，当使用小流量比（Q_r =0.194和0.276）时，以及θ=45°和60°时，离汇流处较远的位置形成了泥沙淤积。相反，在θ=75°和90°的情况下，当实验室水槽中使用较小的流量比时，泥沙淤积出现在了汇合点，但流量比的逐渐增大会使泥沙淤积位置再次向下游移动。Shakibainia等^［19］指出，靠近汇合通道内壁的泥沙淤积的形成是由于流动分离区的压力降低所致。在当前的研究中，观察结果表明，当θ减小时，在距离汇合后河道起点较远的地方会产生泥沙淤积。另外，对于较缓的岸坡角，汇流下游似乎有一个较小的分离区，这是导致在距汇流后河道起点较远的位置形成泥沙淤积的原因。为检验本研究的有效性，本文使用Ghobadian和Shafai-bajestan的实验数据来验证本文的公式（7）和公式（8），对比结果见图5和图6。如图所示，两者拟合效果较好，误差较小。

图5 Ghobadian & Shafai-bajestan实验数据验证公式（7）

Fig.5 Validation of formula （7） by Ghobadian & Shafai-bajestan

Full size|PPT slide

图6 Ghobadian & Shafai-bajestan实验数据验证公式（8）

Fig.6 Validation of formula （8） by Ghobadian & Shafai-bajestan

Full size|PPT slide

3 结论

本文研究了主河道岸坡角（θ）、流量比（Q_r ）和密度弗劳德数（Fr _g3）3个变量对汇流处冲刷和泥沙淤积的影响，得出以下结论：①在汇流处确定小流量比时，岸坡角对冲刷深度的影响最大；②对于给定的汇流处支流量，汇流后Fr _g3的增加使得最大冲刷深度增加，并降低了淤积高度；③岸坡角对淤积高度的影响大于对冲刷深度的影响；④本研究中，θ=45°和Q_r =0.194的试验效果最好，使最大冲刷深度减少了46%。 □

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	杨春瑞，邓金运，齐永铭，等. 中小洪水调度对三峡水库泥沙淤积的长期影响［J］. 水电能源科学，2020，38（6）：34-37. 本文引用 [1]

2	李霞，侯极. 入汇角60°时入汇口泥沙淤积试验研究［J］. 水电能源科学，2011，29（8）：75-76，212. 本文引用 [1]

3	胡成，杨校礼，杨林，等. 双弧形导沙坎相对高度对取水口冲沙效果的影响［J］. 水电能源科学，2021，39（3）：70-73. 本文引用 [1]

4	NAZARI-GIGLOU A， JABBARI-SAHEBARI A， SHAKIBAEINIA A， BORGHEI S M. An experimental study of sediment transport in channel confluences［J］. International Journal of Sediment Research， 2016，31（1）：87-96. 本文引用 [2]

5	YUAN S， TANG H， XIAO Y， et al. Water flow and sediment transport at open-channel confluences： an experimental study［J］. Journal of Hydraulic Research， 2017，56（3）：333-350. 本文引用 [2]

6	卜慧，邵骏，欧阳硕，等. 改进SCS模型在老挝南乌河流域中的应用［J］. 人民长江，2018，49（22）：88-92．本文引用 [1]

7	李泽刚. 黄河近代河口演变基本规律与稳定流路治理［M］. 郑州：黄河水利出版社，2006：164-191. 本文引用 [1]

8	王恺忱. 黄河河口的演变与治理［M］. 郑州：黄河水利出版社，2010：195-259. 本文引用 [1]

9	胡春宏，曹文洪. 黄河口水沙变异与调控I：黄河口水沙运动与演变基本规律［J］. 泥沙研究，2003，28（5）：1-8. 本文引用 [1]

10	胡春宏，曹文洪. 黄河口水沙变异与调控Ⅱ：黄河口治理方向与措施［J］. 泥沙研究，2003，28（5）：9-14. 本文引用 [1]

11	余欣，吉祖稳，王开荣，等. 黄河河口演变与流路稳定关键技术研究［J］. 人民黄河，2020，42（9）：66-70. 本文引用 [1]

12	王崇浩，曹文洪，张世奇. 黄河口潮流与泥沙输移过程的数值研究［J］. 水利学报，2008，39（10）：1 256-1 263. 本文引用 [1]

13	王崇浩，张世奇，曹文洪. 黄河口演变对下游河道反馈影响的计算研究［J］. 水利水电技术，2007，38（9）：16-20.

14	王崇浩，安催花，曹文洪，等. 黄河口清水沟流路冲淤发展与使用年限预测［J］. 水利学报，2008，39（6）：709-713. 本文引用 [1]

15	钟磊，周迎，刘畅，等. 基于GIS的多年土壤侵蚀对安基山水库淤积的影响［J］. 中国农村水利水电，2016（10）：6-10. 本文引用 [1]

16	谭超，黄本胜，邱静，等. 波流联合作用下广东阳西电厂海域二维泥沙数值模拟研究［J］. 中国农村水利水电，2019（2）：24-29. 本文引用 [1]

17	GHOBADIAN R， SHAFAI-BAJESTAN M. Investigation of sediment patterns at river confluence［J］. Journal of Applied Sciences， 2007，7（10）：1 372-1 380. 本文引用 [2]

18	GURRAM S K， KARKI K S， HAGER W H. Subcritical junction flow［J］. Journal of Hydraulic Engineering， 1997，123（5）：447-455. 本文引用 [1]

19	SHAKIBAINIA A， TABATABAI M R， ZARRATI A R. Three-dimensional numerical study of flow structure in channel confluence［J］. Canadian Journal of Civil Engineering， 2010，37（5）：772-781. 本文引用 [1]

PDF(1175 KB)

访问

引用

详细情况

段落导航

收稿日期	出版日期
2021-04-23	2022-01-15
发布日期
2022-02-13

选择文件类型/文献管理软件名称

选择包含的内容