沙漠高含沙洪水沙坝堵河的实验模拟
1.
2.
Experimental simulations on river jamming by hyper-concentrated flows from desert tributaries
1.
2.
收稿日期: 2025-11-03 修回日期: 2026-02-09
| 基金资助: |
|
Received: 2025-11-03 Revised: 2026-02-09
作者简介 About authors
张晓培(1999—),女,河南汝州人,硕士研究生,研究方向为生态环境学E-mail:
关键词:
Keywords:
本文引用格式
张晓培, 拓万全.
Zhang Xiaopei, Ta Wanquan.
0 引言
沙漠季节性河道常常爆发超高浓度的高含沙洪水。该高含沙洪水具有高含沙(1 000~1 600 kg·m-3)、粗粒径、强侵蚀、陡涨快落、历时短等特性[1-4]。高含沙洪水入河,常引起河道淤堵、分流、迁徙、决口改道、阻塞成湖等灾害[1,5-17],一直是沙漠河流治理的重点和难点。典型的沙漠河道是鄂尔多斯高原的十大孔兑。十大孔兑高含沙洪水发生频繁,入黄泥沙通量大,对黄河“几字弯”河道水沙关系和河道演变产生严重影响,已发育“新悬河”近300 km。“新悬河”河段洪凌灾害频发,多次出现决口灾害事件:1989年“7·12”十大孔兑高含沙洪流堵塞黄河决口事件,2003年9月洪峰流量1 300 m3·s-1导致大河弯决口事件,2008年3月凌汛流量1 600 m3·s-1造成杭锦旗决口事件等。近年来,随着黄河径流量快速减小,水沙关系恶化,径流输沙能力降低,主河槽严重萎缩,进一步抬高“新悬河”河床[18-20]。特别是,受全球气候快速变化和强人类活动作用和影响,频发的大风沙尘暴和极端暴雨会加大孔兑风沙水沙入黄通量,进一步加剧“新悬河”决口风险。
但是,沙漠高含沙洪水入河历时短、破坏力强,难以在实地进行观测。因此,高含沙洪水沙坝堵河机理一直不明。一些研究在室内水槽应用粉煤灰材料实验模拟高含沙洪水入河沙坝形成过程[18]。由于粉煤灰密度是泥沙密度的1/4,在水流中运移规律与真实泥沙过程不同,其模拟结果不能给出沙坝堵河确定性结论。因此,为正确认识高含沙洪水沙坝堵河机理,本文应用自然粗泥沙和大尺度弯道进行原型过程模拟实验,揭示不同入汇条件下沙坝的形成机制与淤堵积类型,确定沙坝堵河的关键阈值,为沙坝堵河机理认识和精准防控提供关键数据。
1 实验方法
沙坝堵河模拟实验在罕太川银肯沙自建的弯道水槽中完成,主要包括蓄水池、稳流水槽、弯道试验段、高含沙流入汇坡道和水沙排除道五部分(图1)。蓄水池长20 m、宽10 m、高5 m,设计有出水闸门控制弯道试验段的径流量。在闸门与弯道试验段之间设计稳流水槽,长10 m、宽1 m、深1.2 m,由混凝土制成。稳流水槽与弯道试验段之间安装超声波水位测量计,自动记录水位变化。弯道试验段宽3 m、中心线夹角120°。高含沙流发生坡道是由自然沙漠沙丘沙经过平整形成倾角12°的沙坡坡道。应用流量0.167 m3·s-1抽水泵(管径
图1
图1
沙坝堵河模拟实验结构图
Fig.1
Illustration for experimental simulation on river jamming by hyper-concentrated flows
实验设计分为5个控制径流(径流深0.2、0.4、0.6、0.8、1.0 m)和4个高含沙流激发流量(2泵、3泵、4泵和5泵)进行不同组合,进行20种不同情景沙坝堵河实验模拟。由于高含沙流携沙量大,每次实验控制时间2 min。每次实验后,用挖掘机重新铺设高含沙流发生沙漠坡道,并用RIEGL VZ-400扫描记录弯道试验段沙坝分布。之后,清除试验段沙坝,开始下一次模拟实验。实验时高含沙流通过人工采样测量其泥沙含量。应用固定在长1.0 m的不锈钢杆的塑料瓶(体积为480 mL)采集高含沙流发生时悬移质泥沙样品。采样间隔30 s,每次4个取样。对样品进行烘干,称重,计算高含沙流的泥沙含量。通过水位测定,计算稳流水槽过水断面的水力半径和平均流速:
式中:R是湿周长;h是监测水位;A是过水面积;V是流速;S是弯道比降;n是曼宁粗糙系数0.024。
2 实验结果
2.1 高含沙流入汇口相干涡流调控附岸沙坝形成发育
图2
图2
沙漠坡道模拟高含沙流泥沙浓度
Fig.2
Suspended sediment loads of hyper-concentrated flows on a slope covered with aeolian sands
图3
图3
高含沙流入汇导致附岸沙坝形成发育
Fig.3
Formation of attached-bank sandbars induced by hyper-concentrated flows
图4
图4
不同情景模拟条件下弯道附岸沙坝变化过程
Fig.4
Change of attached-bank sandbars response to different flow conditions and hyper-concentrated flows
2.2 高含沙流促使附岸沙坝堵河临界
在30次不同情景下高含沙流入汇弯道水流中,有17次发生沙坝完全堵塞河道现象,另外13次形成的沙坝尚未完全堵塞弯道(表1),其中17次是由ACW-B沙坝增大造成的“下堵河”情景,另外3次是由CW-B沙坝造成的“上堵河”情景。由于每次模拟的高含沙流泥沙浓度变化不显著,分析高含沙流与弯道水流径流比QHCF/Q0与河道淤堵情景关系指出,当QHCF/Q0>18.0时,发生“上堵河”情景。“下堵河”情景相对计较复杂,实验结果表明,当0.5<QHCF /Q0<6.0时,“下堵河”情景容易发生。实验中,没有观测到“上堵河”和 “下堵河”现象同时出现情景。
表1 30次不同高含沙流入河沙坝堵河情景实验模拟
Table 1
| WP | WH/cm | HCFC/(kg·m-3) | QHCF/Q0 | SB/m3 | HCFS/m3 | SB/HCFS/% | UJ | DJ |
|---|---|---|---|---|---|---|---|---|
| 2 | 10 | 961.34 | 12.12 | 7.6789 | 34.0708 | 22 | 否 | 否 |
| 20 | 1 059.54 | 3.89 | 8.4992 | 23.8024 | 36 | 否 | 否 | |
| 30 | 872.12 | 2.02 | 11.9705 | 19.7517 | 61 | 否 | 是 | |
| 40 | 934.75 | 1.27 | 8.2514 | 15.1008 | 54 | 否 | 是 | |
| 50 | 1 012.70 | 0.89 | 9.7061 | 29.1571 | 33 | 否 | 是 | |
| 60 | 930.08 | 0.67 | 9.8417 | 28.2447 | 35 | 否 | 是 | |
| 70 | 906.41 | 0.53 | 7.9007 | 37.5569 | 21 | 否 | 否 | |
| 80 | 927.46 | 0.43 | 9.5297 | 25.5026 | 37 | 否 | 否 | |
| 3 | 10 | 962.79 | 18.18 | 7.6169 | 15.0006 | 51 | 是 | 否 |
| 20 | 885.39 | 5.83 | 13.1337 | 28.5397 | 46 | 否 | 是 | |
| 30 | 901.88 | 3.02 | 12.8224 | 49.7279 | 26 | 否 | 是 | |
| 40 | 781.22 | 1.91 | 12.7980 | 38.6005 | 33 | 否 | 是 | |
| 50 | 829.34 | 1.34 | 12.8495 | 33.8284 | 38 | 否 | 是 | |
| 60 | 843.77 | 1.01 | 6.2532 | 19.7201 | 32 | 否 | 是 | |
| 70 | 807.03 | 0.79 | 7.7346 | 35.3857 | 22 | 否 | 是 | |
| 80 | 810.07 | 0.65 | 9.7617 | 31.3569 | 31 | 否 | 是 | |
| 4 | 10 | 763.71 | 24.24 | 6.4118 | 40.8566 | 16 | 是 | 否 |
| 20 | 821.44 | 7.78 | 8.6510 | 49.4304 | 17 | 否 | 否 | |
| 30 | 843.32 | 4.03 | 15.8663 | 62.7564 | 25 | 否 | 是 | |
| 40 | 778.79 | 2.55 | 11.2110 | 32.4595 | 34 | 否 | 是 | |
| 50 | 836.94 | 1.79 | 5.9749 | 49.3597 | 12 | 否 | 是 | |
| 60 | 997.72 | 1.34 | 3.2203 | 51.4487 | 6 | 否 | 否 | |
| 70 | 779.18 | 1.06 | 6.6479 | 40.6257 | 16 | 否 | 否 | |
| 80 | 822.22 | 0.86 | 17.6220 | 56.4646 | 31 | 否 | 是 | |
| 5 | 10 | 893.72 | 30.30 | 5.7588 | 38.9918 | 15 | 是 | 否 |
| 20 | 934.18 | 9.73 | 5.7647 | 61.3761 | 9 | 否 | 否 | |
| 40 | 876.55 | 3.18 | 6.3266 | 63.4180 | 9 | 否 | 否 | |
| 60 | 916.16 | 1.68 | 3.5362 | 15.5460 | 23 | 否 | 否 | |
| 70 | 845.15 | 1.32 | 6.9077 | 32.9272 | 21 | 否 | 否 | |
| 80 | 855.49 | 1.08 | 12.2517 | 44.6597 | 27 | 否 | 是 |
3 讨论
本实验模拟发现高含沙流入汇口相干涡流导致附岸沙坝堵河的新机制。该机制与应用粉煤灰模拟的异重流堵河机制认识不同[21]。粉煤灰模拟的主要缺点在于材料密度(700 kg·m-3)太小、颗粒粒径(30~40
本实验模拟指出沙坝堵河与汇流比密切相关。当高含沙流泥沙含量保持基本稳定时,汇流比是高含沙洪水沙坝堵河的主要决定因素。“下堵河”型是高含沙流沙坝堵河的主要形式,实验结果表明,当0.5<QHCF /Q0<6.0时,“下堵河”情景容易发生。“上堵河”情景一般发生较少,需要大流量高含沙洪水入汇,一般需要QHCF /Q0>18.0时,才能发生,这在自然河流中发生的可能性很小。因此,沙坝堵河主要发生的是“下堵河”,而不是之前认为的“上堵河”。
本实验模拟揭示出高含沙洪水入汇与低含沙水流入汇在交汇区存在显著的动力差异(图5)。低含沙水流入汇常常形成相对低流速带,高流速带与下游交汇角岸壁之间形成逆时针环流分离区。模拟结果显示,高含沙流汇入与弯道水流相互作用,在交汇线附近的两侧河岸形成相干涡流对,相干涡流对引起其控制区域水流流速快速下降,大量泥沙快速淤积形成逆时针回流附岸沙坝ACW-B和顺时针回流附岸沙坝CW-B。尽管ACW-B沙坝与低含沙流引起的逆时针分离区形成的沙坝类似,但是,其成因机制却完全不同。因此,未来在沙坝堵河机理研究和精准防控方面,需要考虑相干涡流导致附岸沙坝堵河的新机制、新效应。
图5
图5
低含沙流与高含沙流入河交汇区动力机制差异
Fig.5
Difference in sandbar formation induced by hyper-concentrated flow and low-concentrated flow
4 结论
高含沙流入汇口存在相干涡流对,调控交汇线两侧河岸泥沙快速沉积形成顺时针回流附岸沙坝CW-B和逆时针回流附岸沙坝ACW-B,是造成沙坝堵河的新动力机制。
沙坝堵河主要发生的是“下堵河”,而不是之前认为的“上堵河”。当0.5<QHCF /Q0<6.0时,“下堵河”情景容易发生。当QHCF /Q0>18.0时,“上堵河”现象才能发生。
高含沙流输入弯道中的泥沙约有28%在入汇口沉积形成沙坝,是沙坝堵河的主要原因。
参考文献
Aeolian process-induced hyper-concentrated flow in a desert watershed
[J].
毛不拉孔兑暴雨径流产输沙特征研究
[J].
高流速含沙浑水体流变特性试验研究
[J].
黄河下游高含沙洪水演进及河床冲淤过程的数值模拟
[J].
“揭河底”冲刷现象的分析
[J].
黄河上游沙漠宽谷段高含沙支流对干流的淤堵影响
[J].
高含沙洪水“揭河底”的判别指标及其条件
[J].
黄河高含沙洪水“揭河底”机理探讨
[J].
黄河上游高含沙水流交汇区河床形态及冲淤变化规律
[J].
黄河上游干支流交汇区沙坝淤堵形成条件
[J].
高含沙水流特性和滩槽冲淤特征浅述
[J].
黄河“揭河底”判别指标理论研究及验证
[J].
Lacustrine sediments at Narabeb in the Central Namib Desert,Namibia
[J].
Waterholes and their significance in the Anastomosing Channel System of Cooper Creek,Australia
[J].
Mapping abandoned river channels in mali through directional filtering of thematic mapper data
[J].
黄河泥沙百年演变特征与近期波动变化成因解析
[J].
An experimental study of fluvial processes at asymmetrical riverconfluences with hyperconcentrated tributary flow
[J].
/
| 〈 |
|
〉 |

甘公网安备 62010202000688号