• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊

留言板

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

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

非恒定流驱动下床面形态变化特征

王乐 AlanJ S Cuthbertson 张尚弘 王永强

王乐, AlanJ S Cuthbertson, 张尚弘, 王永强. 非恒定流驱动下床面形态变化特征[J]. 水科学进展, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
引用本文: 王乐, AlanJ S Cuthbertson, 张尚弘, 王永强. 非恒定流驱动下床面形态变化特征[J]. 水科学进展, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
WANG Le, Alan J S Cuthbertson, ZHANG Shanghong, WANG Yongqiang. Study on the characteristics of bed forms developed under unsteady flows[J]. Advances in Water Science, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
Citation: WANG Le, Alan J S Cuthbertson, ZHANG Shanghong, WANG Yongqiang. Study on the characteristics of bed forms developed under unsteady flows[J]. Advances in Water Science, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010

非恒定流驱动下床面形态变化特征

doi: 10.14042/j.cnki.32.1309.2021.05.010
基金项目: 

国家自然科学基金资助项目 52009041

国家自然科学基金资助项目 52061135104

详细信息
  • 中图分类号: TV122

Study on the characteristics of bed forms developed under unsteady flows

Funds: 

the National Natural Science Foundation of China 52009041

the National Natural Science Foundation of China 52061135104

  • 摘要: 为深入认识非恒定流驱动下床面形态的演变规律及尺度变化趋势,采用水泵控制系统生成了一系列历时与洪峰流量不同的洪水过程来冲刷由粗沙组成的实验动床(中值粒径d50=1.95 mm),研究中运用非恒定流的形态、非恒定性及水流做功参数来量化洪水过程;统计分析河床高程数据得到床面形态的波长、波高及陡度。实验结果表明:在不同尺度的非恒定流洪水作用下,实验动床表面分别形成了沙垄、交错边滩以及介于两者之间的过渡型床面形态;洪水的非恒定性对床面形态尺度的影响最大,其次是水流做功参数,而非恒定流的形态影响较为微弱;在此基础上构建了用于描述非恒定流驱动条件与床面形态尺度响应的定量关系。
  • 图  1  研究床面形态的水槽实验及主要观测仪器

    Figure  1.  Experimental flume and main equipment employed for studying bed forms developed over a mobile bed

    图  2  实验序列U1洪水作用下的床面形态与高程变化

    Figure  2.  Detailed bed elevation maps measured for post-hydrograph bed surface in representative tests within experimental series U1

    图  3  非恒定流形态、非恒定性及水流做功参数对床面形态波长(平均值+标准偏差)的影响

    Figure  3.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form wavelength (averaged value±standard deviation)

    图  4  非恒定流形态、非恒定性及水流做功参数对床面形态波高(平均值+标准偏差)的影响

    Figure  4.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form height (averaged value±standard deviation)

    图  5  非恒定流形态、非恒定性及水流做功参数对床面形态陡度(平均值+标准偏差)的影响

    Figure  5.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form steepness (averaged value±standard deviation)

    图  6  非恒定流综合参数ΓHGWk0.4与床面形态陡度的关系

    Figure  6.  Relationship between hydrograph combined parameter ΓHGWk0.4 with bed-form steepness ψ

    表  1  实验中的非恒定流条件

    Table  1.   Summary of unsteady flow hydrographs in the present experiment

    序列 组次 Qp/(l·s-1) Hp/cm ΔT/s S ΓHG/10-4 Wk η
    U1 U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
    S1b 58.0 12.01 14 400 0.002 1.389 231.51 0.40
    S1c 58.0 12.01 14 400 0.002 1.389 231.51 2.50
    U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
    V1b 50.0 10.98 11 400 0.002 1.467 143.64 1.00
    V1c 40.0 9.60 7 800 0.002 1.577 68.47 1.00
    V1d 35.0 8.86 6 000 0.002 1.657 41.32 1.00
    V1e 28.0 7.75 3 600 0.002 1.776 15.26 1.00
    U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
    U1b 50.0 10.98 18 000 0.002 0.929 226.12 1.00
    U1c 45.0 10.30 21 600 0.002 0.674 233.75 1.00
    U1d 41.0 9.75 25 200 0.002 0.507 233.61 1.00
    U1e 34.0 8.71 36 000 0.002 0.263 231.57 1.00
    U2 U2a/V2a 58.0 12.01 7 200 0.002 2.778 117.15 1.00
    V2b 40.0 9.60 3 900 0.002 3.154 34.24 1.00
    V2c 28.0 7.75 1 800 0.002 3.553 7.63 1.00
    U2a/V2a 58.0 12.01 7 200 0.002 2.778 117.15 1.00
    U2b 38.0 9.30 14 400 0.002 0.789 117.15 1.00
    U2c 28.0 7.80 27 000 0.002 0.237 117.49 1.00
    U3 U3a/V3a 58.0 12.00 3 600 0.002 5.555 58.57 1.00
    V3b 40.0 9.60 1 950 0.002 6.308 17.12 1.00
    V3c 28.0 7.80 900 0.002 7.105 3.81 1.00
    U3a/V3a 58.0 12.00 3 600 0.002 5.555 58.57 1.00
    U3b 38.0 9.31 7 200 0.002 1.579 58.57 1.00
    U3c 28.0 7.75 13 500 0.002 0.474 58.74 1.00
    注: Qp为洪峰流量;Hp为洪峰流量对应的水深;S为坡度。
    下载: 导出CSV
  • [1] PARSONS D R, BEST J L, ORFEO O, et al. Morphology and flow fields of three-dimensional dunes, Rio Paraná, Argentina: results from simultaneous multibeam echo sounding and acoustic Doppler current profiling[J]. Journal of Geophysical Research: Earth Surface, 2005, 110(F4): F04S03. doi:  10.1029/2004JF000231
    [2] SIMON A, RINALDI M. Disturbance, stream incision, and channel evolution: the roles of excess transport capacity and boundary materials in controlling channel response[J]. Geomorphology, 2006, 79(3/4): 361-383. http://www.sciencedirect.com/science/article/pii/S0169555X06002583
    [3] AHANGER M A, ASAWA G L, LONE M A. Experimental study of sediment transport hysteresis[J]. Journal of Hydraulic Research, 2008, 46(5): 628-635. doi:  10.3826/jhr.2008.3185
    [4] GRAF W H, QU Z. Flood hydrographs in open channels[J]. Proceedings of the Institution of Civil Engineers-Water Management, 2004, 157(1): 45-52. doi:  10.1680/wama.2004.157.1.45
    [5] 周俊伟, 包为民, 孙逸群, 等. 河道一维非恒定流微分结构摩阻公式及其试验论证[J]. 水科学进展, 2017, 28(4): 571-578. doi:  10.14042/j.cnki.32.1309.2017.04.011

    ZHOU J W, BAO W M, SUN Y Q, et al. Analysis of 1-D unsteady flow friction in open channel flows and experimental demonstration[J]. Advances in Water Science, 2017, 28(4): 571-578. (in Chinese) doi:  10.14042/j.cnki.32.1309.2017.04.011
    [6] 刘春晶, 曲兆松, 李丹勋, 等. 明渠非恒定流推移质输沙试验研究[J]. 水力发电学报, 2006, 25(2): 31-37. doi:  10.3969/j.issn.1003-1243.2006.02.007

    LIU C J, QU Z S, LI D X, et al. Experimental study on bed-load transport in unsteady open channel flows[J]. Journal of Hydroelectric Engineering, 2006, 25(2): 31-37. (in Chinese) doi:  10.3969/j.issn.1003-1243.2006.02.007
    [7] 程小兵, 王艳华, 乐培九. 非恒定流清水冲刷输沙规律初步研究[J]. 泥沙研究, 2011(3): 10-16. https://www.cnki.com.cn/Article/CJFDTOTAL-NSYJ201103003.htm

    CHENG X B, WANG Y H, YUE P J. Preliminary study of bed load sediment transport of unsteady flow[J]. Journal of Sediment Research, 2011(3): 10-16. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NSYJ201103003.htm
    [8] 马爱兴, 陆永军, 陆彦, 等. 非恒定流作用下砾石推移质输移特性试验研究[J]. 水利学报, 2013, 44(7): 800-809. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201307009.htm

    MA A X, LU Y J, LU Y, et al. Experimental study on gravel bed-load transport in unsteady flow[J]. Journal of Hydraulic Engineering, 2013, 44(7): 800-809. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201307009.htm
    [9] 吴国茂, 胡江, 杨胜发. 明渠非恒定流下输沙对比研究[J]. 水运工程, 2015(5): 117-121. doi:  10.3969/j.issn.1002-4972.2015.05.021

    WU G M, HU J, YANG S F. Comparative study on sediment transport under unsteady flow of open channel[J]. Port & Waterway Engineering, 2015(5): 117-121. (in Chinese) doi:  10.3969/j.issn.1002-4972.2015.05.021
    [10] LEE K T, LIU Y L, CHENG K H. Experimental investigation of bedload transport processes under unsteady flow conditions[J]. Hydrological Processes, 2004, 18(13): 2439-2454. doi:  10.1002/hyp.1473
    [11] BOMBAR G, ELÇI Ş, TAYFUR G, et al. Experimental and numerical investigation of bed-load transport under unsteady flows[J]. Journal of Hydraulic Engineering, 2011, 137(10): 1276-1282. doi:  10.1061/(ASCE)HY.1943-7900.0000412
    [12] 段自豪, 陈杰, 蒋昌波, 等. 非恒定流作用下的推移质泥沙输移实验研究[J]. 中国科学: 技术科学, 2019, 49(11): 1372-1382. https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201911012.htm

    DUAN Z H, CHEN J, JIANG C B, et al. Experimental study on bedload sediment transport under unsteady flow[J]. Scientia Sinica Technologica, 2019, 49(11): 1372-1382. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201911012.htm
    [13] WANG L, CUTHBERTSON A J S, ZHANG S H, et al. Graded bed load transport in sediment supply limited channels under unsteady flow hydrographs[J]. Journal of Hydrology, 2021, 595: 126015. doi:  10.1016/j.jhydrol.2021.126015
    [14] MARTIN R L, JEROLMACK D J. Origin of hysteresis in bed form response to unsteady flows[J]. Water Resources Research, 2013, 49(3): 1314-1333. doi:  10.1002/wrcr.20093
    [15] 孙东坡, 刘明潇, 张晓雷, 等. 冲积性河流河床冲淤调整对洪水泥沙过程的响应: 以黄河游荡型河段为例[J]. 水科学进展, 2014, 25(5): 668-676. http://skxjz.nhri.cn/article/id/2452

    SUN D P, LIU M X, ZHANG X L, et al. Alluvial river adjusting response to the flood-sediment process: for the wandering reach of Yellow River as an example[J]. Advances in Water Science, 2014, 25(5): 668-676. (in Chinese) http://skxjz.nhri.cn/article/id/2452
    [16] WATERS K A, CURRAN J C. Linking bed morphology changes of two sediment mixtures to sediment transport predictions in unsteady flows[J]. Water Resources Research, 2015, 51(4): 2724-2741. doi:  10.1002/2014WR016083
    [17] REDOLFI M, WELBER M, CARLIN M, et al. Morphometric properties of alternate bars and water discharge: a laboratory investigation[J]. Earth Surface Dynamics, 2020, 8(3): 789-808. doi:  10.5194/esurf-8-789-2020
    [18] 赵维阳, 杨云平, 张华庆, 等. 三峡大坝下游近坝段沙质河床形态调整及洲滩联动演变关系[J]. 水科学进展, 2020, 31(6): 862-874. doi:  10.14042/j.cnki.32.1309.2020.06.006

    ZHAO W Y, YANG Y P, ZHANG H Q, et al. Adjustment patterns and causes of the morphology of sandy riverbed downstream of the Three Gorges Dam[J]. Advances in Water Science, 2020, 31(6): 862-874. (in Chinese) doi:  10.14042/j.cnki.32.1309.2020.06.006
    [19] GUNSOLUS E H, BINNS A D. Effect of morphologic and hydraulic factors on hysteresis of sediment transport rates in alluvial streams[J]. River Research and Applications, 2018, 34(2): 183-192. doi:  10.1002/rra.3184
    [20] 马爱兴, 陆彦, 陆永军, 等. 明渠非恒定流流速分布及推移质运动研究进展[J]. 水科学进展, 2012, 23(1): 134-142. http://skxjz.nhri.cn/cn/article/doi/CNKI:32.1309.P.20120104.2014.013

    MA A X, LU Y, LU Y J, et al. Advances in research on velocity distribution and bed-load transport in unsteady open-channel flow[J]. Advances in Water Science, 2012, 23(1): 134-142. (in Chinese) http://skxjz.nhri.cn/cn/article/doi/CNKI:32.1309.P.20120104.2014.013
    [21] YEN C L, LEE K T. Bed topography and sediment sorting in channel bend with unsteady flow[J]. Journal of Hydraulic Engineering, 1995, 121(8): 591-599. doi:  10.1061/(ASCE)0733-9429(1995)121:8(591)
    [22] WANG L. Bed-load sediment transport and bed evolution in steady and unsteady flows[D]. Edinburgh: Heriot-Watt University, 2016.
  • [1] 肖诗云, 郭晓阳, 李娜.  洪水水流压力及其影响因素模型试验研究 . 水科学进展, 2021, 32(3): 418-426. doi: 10.14042/j.cnki.32.1309.2021.03.010
    [2] 赵汗青, 唐洪武, 闫静, 戴会超, 刘志武.  淹没植物明渠床面冲淤及其对水流运动的影响 . 水科学进展, 2021, 32(2): 250-258. doi: 10.14042/j.cnki.32.1309.2021.02.010
    [3] 刘易庄, 杨树清, 蒋昌波, 隆院男, 邓斌.  洞庭湖萎缩对湖内洪水影响 . 水科学进展, 2019, 30(4): 485-495. doi: 10.14042/j.cnki.32.1309.2019.04.004
    [4] 曹斌, 邹燚, 夏建新.  不同粒度颗粒在非恒定管流中跟随性试验 . 水科学进展, 2017, 28(3): 356-363. doi: 10.14042/j.cnki.32.1309.2017.03.005
    [5] 周俊伟, 包为民, 孙逸群, 张乾.  河道一维非恒定流微分结构摩阻公式及其试验论证 . 水科学进展, 2017, 28(4): 571-578. doi: 10.14042/j.cnki.32.1309.2017.04.011
    [6] 程永舟, 潘昀, 蒋昌波, 李青峰.  破碎波作用下沙质海床床面形态变化试验 . 水科学进展, 2014, 25(2): 253-259.
    [7] 张原锋, 申冠卿, Verbanck M A.  黄河下游床面形态判别方法探讨 . 水科学进展, 2012, 23(1): 46-52. doi: CNKI:32.1309.P.20111125.1635.006
    [8] 赖锡军, 傅国圣, 孙波.  非恒定水流计算的最优控制问题及其变分求解 . 水科学进展, 2008, 19(4): 537-545.
    [9] 关铁生, 王家祁.  20世纪暴雨和洪水极值的年代际变化 . 水科学进展, 2006, 17(1): 89-94.
    [10] 范杰, 王长德, 管光华, 崔巍.  渠道非恒定流水力学响应研究 . 水科学进展, 2006, 17(1): 55-60.
    [11] 冯小香, 张小峰, 崔占峰.  垂向二维非恒定流及悬浮物分布模型研究 . 水科学进展, 2006, 17(4): 518-524.
    [12] 马宗伟, 许有鹏, 李嘉峻.  河流形态的分维及与洪水关系的探讨——以长江中下游为例 . 水科学进展, 2005, 16(4): 530-534.
    [13] 陈前海, 方红卫, 王光谦.  三门峡库区一维非恒定非均匀泥沙输移数学模型 . 水科学进展, 2004, 15(2): 160-164.
    [14] 唐勇, 胡和平, 田富强, 卢麾.  鳌江流域洪水预报调度决策支持系统 . 水科学进展, 2003, 14(3): 295-299.
    [15] 许炯心.  黄河下游洪水的泥沙输移特征 . 水科学进展, 2002, 13(5): 562-568.
    [16] 张红武, 黄远东, 赵连军, 江恩惠.  黄河下游非恒定输沙数学模型——Ⅱ模型验证 . 水科学进展, 2002, 13(3): 271-277.
    [17] 魏根群, 陈璧宏, 宿晓辉.  反调节水库非恒定流数值模拟 . 水科学进展, 2001, 12(4): 491-498.
    [18] 程永光, 索丽生.  计算明渠非恒定流的一个Lattice BGK模型 . 水科学进展, 2000, 11(4): 362-367.
    [19] 魏一鸣, 杨存键, 金菊良.  洪水灾害分析与评估的综合集成方法 . 水科学进展, 1999, 10(1): 25-30.
    [20] 王家祁, 沈国昌, 耿雷华, 顾文燕, 姚惠明.  中国热带气旋暴雨洪水的分布和水文特性 . 水科学进展, 1995, 6(2): 121-126.
  • 加载中
图(6) / 表 (1)
计量
  • 文章访问数:  353
  • HTML全文浏览量:  44
  • PDF下载量:  63
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-29
  • 网络出版日期:  2021-07-19
  • 刊出日期:  2021-09-30

非恒定流驱动下床面形态变化特征

doi: 10.14042/j.cnki.32.1309.2021.05.010
    基金项目:

    国家自然科学基金资助项目 52009041

    国家自然科学基金资助项目 52061135104

  • 中图分类号: TV122

摘要: 为深入认识非恒定流驱动下床面形态的演变规律及尺度变化趋势,采用水泵控制系统生成了一系列历时与洪峰流量不同的洪水过程来冲刷由粗沙组成的实验动床(中值粒径d50=1.95 mm),研究中运用非恒定流的形态、非恒定性及水流做功参数来量化洪水过程;统计分析河床高程数据得到床面形态的波长、波高及陡度。实验结果表明:在不同尺度的非恒定流洪水作用下,实验动床表面分别形成了沙垄、交错边滩以及介于两者之间的过渡型床面形态;洪水的非恒定性对床面形态尺度的影响最大,其次是水流做功参数,而非恒定流的形态影响较为微弱;在此基础上构建了用于描述非恒定流驱动条件与床面形态尺度响应的定量关系。

English Abstract

王乐, AlanJ S Cuthbertson, 张尚弘, 王永强. 非恒定流驱动下床面形态变化特征[J]. 水科学进展, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
引用本文: 王乐, AlanJ S Cuthbertson, 张尚弘, 王永强. 非恒定流驱动下床面形态变化特征[J]. 水科学进展, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
WANG Le, Alan J S Cuthbertson, ZHANG Shanghong, WANG Yongqiang. Study on the characteristics of bed forms developed under unsteady flows[J]. Advances in Water Science, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
Citation: WANG Le, Alan J S Cuthbertson, ZHANG Shanghong, WANG Yongqiang. Study on the characteristics of bed forms developed under unsteady flows[J]. Advances in Water Science, 2021, 32(5): 751-758. doi: 10.14042/j.cnki.32.1309.2021.05.010
  • 冲积河流中存在着形态各异、结构复杂的床面形态, 如沙垄、交错边滩等, 当床面形态出现时, 会产生形态阻力, 构成对近床面水流结构的扰动, 影响泥沙输移的强度及特征[1-3]。实际上, 河流系统中的泥沙运动及河床演变主要发生在以洪水为代表的强非恒定流过程中。目前, 围绕明渠非恒定流开展的研究主要集中在3个方面。第一, 明渠非恒定流水动力学研究, 如Graf和Qu[4]以洪水中的水动力要素为关注对象, 通过理论及实验方法分析水流参数极值出现的先后顺序;周俊伟等[5]运用量纲分析得到了以微分结构表达的摩阻改进公式。第二, 非恒定流作用下泥沙运动特征为主的研究, 如刘春晶等[6]的实验表明非恒定流的推移质输沙强度与水流强度变化存在不同步的现象;程小兵等[7]、马爱兴等[8]由实验观测认为非恒定流不存在传统意义上的饱和输沙与平衡输沙, 输沙率始终存在滞后现象;吴国茂等[9]分析了Lee等[10]和Bombar等[11]的输沙数据后提出一个综合水沙运动的非恒定性系数来进行拟合;段自豪等[12]分析实验结果得出非恒定流作用下的泥沙颗粒更易起动的结论;Wang等[13]在实验中证实了非恒定流非均匀沙床中不同粒径组间最大输沙强度的滞后性。第三, 非恒定流河床演变的研究, 如Martin和Jerolmack[14]模拟了流量增加或减小时的河床响应, 基于床面形态陡度不变且输沙无绳套特性的矛盾假设建立了床面形态调整时间的预测模型;孙东坡等[15]研究表明黄河下游主槽冲刷总量随含沙量增大而减小, 滩地淤积量随含沙量增大而增大;Waters和Curran[16]实验观测到沙波在小洪水事件中形成;Redolfi等[17]用理论模型描述了洪水波中交错边滩随时间的变化特征;赵维阳等[18]分析三峡大坝下游的河床形态时发现, 随着来沙量减少, 河床呈冲刷态势, 洲滩面积也在同步减少。显然, 现阶段以非恒定流输沙成果居多,而非恒定流河床演变的研究则只关注某一特定类型的床面形态(如交错边滩)或以河道断面形态为主, 真正涉及多尺度洪水过程与多类型床面形态的关联性研究并不多见[19-20]

    鉴于此, 本文拟通过水槽实验来生成一系列具有典型代表性的单波非恒定流过程, 探究实验动床在不同尺度非恒定流驱动下的床面形态变化特征, 定量分析非恒定流洪水尺度与床面形态尺度之间的响应关系, 研究成果将为揭示洪水作用下的河床形态变化趋势及河流演变规律奠定基础。

    • 非恒定性参数(ΓHG)是描述非恒定流的基本参数, 是基流与洪峰流间的水深差(ΔH)与洪水总历时(ΔT)的函数, 其表达式为

      $$ {\mathit{\Gamma }_{{\rm{HG}}}} = \frac{1}{{{u_{*{\rm{b}}}}}}\frac{{\mathit{\Delta }H}}{{\mathit{\Delta }T}} $$ (1)

      式中: u*b为基流剪切流速, m/s;ΔTTRTF, ΔTR、ΔTF分别为洪水上升段与下降段历时。另外, 洪水过程的总水量也会影响到泥沙运动及河床形变。因此, Yen和Lee[21]提出了非恒定流的水流做功概念, 其表达式为

      $$ W_{\mathrm{k}}=\frac{u_{* \mathrm{~b}}^{2} V_{\mathrm{ol}}}{g H_{\mathrm{b}}^{3} B} $$ (2)

      式中: Vol为洪水过程中所包含的总水量, m3, 不包括基流量;Hb为基流水深, m;B为水槽宽度, m。此外, Wang[22]用形态参数来判别对称性与非对称性三角形洪水过程, 其被定义为洪水上升段和下降段历时的比值, 即

      $$ \eta = \frac{{\mathit{\Delta }{T_{\rm{R}}}}}{{\mathit{\Delta }{T_{\rm{F}}}}} $$ (3)

      η=1.0时为对称性洪水过程;当η ≠ 1.0时则为非对称性洪水过程。

    • 研究非恒定流对床面形态的影响, 需要统计分析河床高程数据来确定不同类型的床面形态多维尺度。一般而言, 定量描述床面形态大小的参数主要包括波长(λ)、波高(D)以及陡度(ψ)[1]。波长用于表征相邻2个波峰所代表的最大高程点(或2个相邻波谷所代表的最低高程点)之间的水平距离;对应的波高定义为上述2个波峰所指代的空间范围内高程最高点(波峰)与最低点(波谷)之间的垂直间距;陡度为波高与波长的比值, 即

      $$ \psi=\frac{D}{\lambda} $$ (4)

      显然, 波长、波高是有量纲参量, 陡度为量纲一参数, 三者被广泛用于描述沙垄等水沙界面普遍存在的床面形态。

    • 实验在英国赫瑞瓦特大学水力学实验室的大型水槽内进行。水槽长22 m, 宽0.75 m, 高0.50 m, 断面为矩形, 底板为不锈钢, 边壁由透明玻璃制成。如图 1所示, 该水槽为水流循环的封闭系统, 坡度可手动调节, 水泵设置在水槽入口底端, 可输出最大流量为100 l/s, 水槽下游出口处布设有可调节的百叶窗尾门。实验中使用的测量仪器包括测试段上下游的电子水位计、超声波流量计、调坡仪及变频器等主要设备。其中, 流量由超声波流量计测得, 其测量精度可以达到±0.01 l/s。实验用沙分为A、B、C、D 4种, 其中, d为沙粒直径。

      图  1  研究床面形态的水槽实验及主要观测仪器

      Figure 1.  Experimental flume and main equipment employed for studying bed forms developed over a mobile bed

    • 实验水流为典型非恒定流洪水。非恒定流洪水系列分为3个序列, 即表 1所给出的U1、U2、U3序列, 主要考虑不同的形态参数、非恒定性及水流做功对床面形态变化的影响, 以实验序列U1为例(见表 1), 在非恒定流实验组S1a—S1c中, 其他参数保持不变, η由0.4增大到2.5;同理, 实验组V1a—V1e中只有Wk变化, 而在实验组U1a—U1e中ΓHG变化。实验用沙参照了Lee等[10]的实验设计, 泥沙粒径范围为1.0~2.8 mm, 中值粒径d50=1.95 mm。

      表 1  实验中的非恒定流条件

      Table 1.  Summary of unsteady flow hydrographs in the present experiment

      序列 组次 Qp/(l·s-1) Hp/cm ΔT/s S ΓHG/10-4 Wk η
      U1 U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
      S1b 58.0 12.01 14 400 0.002 1.389 231.51 0.40
      S1c 58.0 12.01 14 400 0.002 1.389 231.51 2.50
      U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
      V1b 50.0 10.98 11 400 0.002 1.467 143.64 1.00
      V1c 40.0 9.60 7 800 0.002 1.577 68.47 1.00
      V1d 35.0 8.86 6 000 0.002 1.657 41.32 1.00
      V1e 28.0 7.75 3 600 0.002 1.776 15.26 1.00
      U1a/V1a/S1a 58.0 12.01 14 400 0.002 1.389 234.29 1.00
      U1b 50.0 10.98 18 000 0.002 0.929 226.12 1.00
      U1c 45.0 10.30 21 600 0.002 0.674 233.75 1.00
      U1d 41.0 9.75 25 200 0.002 0.507 233.61 1.00
      U1e 34.0 8.71 36 000 0.002 0.263 231.57 1.00
      U2 U2a/V2a 58.0 12.01 7 200 0.002 2.778 117.15 1.00
      V2b 40.0 9.60 3 900 0.002 3.154 34.24 1.00
      V2c 28.0 7.75 1 800 0.002 3.553 7.63 1.00
      U2a/V2a 58.0 12.01 7 200 0.002 2.778 117.15 1.00
      U2b 38.0 9.30 14 400 0.002 0.789 117.15 1.00
      U2c 28.0 7.80 27 000 0.002 0.237 117.49 1.00
      U3 U3a/V3a 58.0 12.00 3 600 0.002 5.555 58.57 1.00
      V3b 40.0 9.60 1 950 0.002 6.308 17.12 1.00
      V3c 28.0 7.80 900 0.002 7.105 3.81 1.00
      U3a/V3a 58.0 12.00 3 600 0.002 5.555 58.57 1.00
      U3b 38.0 9.31 7 200 0.002 1.579 58.57 1.00
      U3c 28.0 7.75 13 500 0.002 0.474 58.74 1.00
      注: Qp为洪峰流量;Hp为洪峰流量对应的水深;S为坡度。
    • 实验前, 水槽内铺设厚度均匀的平整沙床(其厚度为11 cm);实验初期, 为使水流能够浸润测试沙床而不造成扰动, 将水槽下游的尾门关闭, 调控水泵生成最小流量(Q=3.0 l/s), 待床面被水流浸没后, 打开尾门, 缓慢增大水泵输出流量, 直至达到基流水平(Q≈17.0 l/s), 基流运行约30 min后, 使用非恒定流控制系统生成预设洪水过程;非恒定流测试结束后, 继续运行基流, 并缓慢关闭尾门, 以产生回水效应, 这样既可保护床面形态,又可增加水深, 有利于借助声学多普勒水流仪(Acoustic Doppler Velocimeter, ADV)来测量河床高程, 精度可达到±0.1 mm。

    • 形态参数变化的非恒定流作用下形成的床面形态如图 2所示。显然, 浅滩或深槽犬牙交替式分布在动床两侧, 形成了交错边滩类型为主的床面形态(S1a, S1b)。在3种形态变化的非恒定流洪水作用后, 河床高程整体降低, 并没有出现局部净淤积(河床高程变化量Δz>0)的情况;冲刷最大深度(≈40 mm)出现在非对称型洪水S1b作用的沙床下游段, 与流量梯度差异对最大输沙率的影响有关。具体而言, 上升段历时较短的洪水过程涉及的正向流量梯度较大, 附加应力增大导致局部河床冲刷更为强烈。

      图  2  实验序列U1洪水作用下的床面形态与高程变化

      Figure 2.  Detailed bed elevation maps measured for post-hydrograph bed surface in representative tests within experimental series U1

      非恒定性变化的洪水过程(U1a, U1c, U1d)对床面形态的影响见图 2。随着洪水非恒定性降低, 即洪峰流量减小且洪水历时延长, 水流对河床的冲刷能力降低, 其塑造大尺度床面形态的能力逐渐减弱。相应地, 床面形态类型也发生了变化, 由较为规则的交错边滩逐渐变为非规则的边滩/沙垄混合的床面形态。图 2也呈现了水流做功变化的洪水过程(V1a, V1c, V1d)对床面形态的影响, 随着Wk的减小, 河床冲刷深度相应减小;同样, 对应的床面形态类型也发生了类似变化, 由交错边滩逐渐变为较规则的沙垄(如V1d)。

    • 实验中各类床面形态的波长、波高、陡度与非恒定流基本参数间的关系详见图 3。显然, 床面形态的波长并没有随洪水形态变化而表现出明显的趋势性变化, 波长维持在1.8~2.8 m范围内;然而, 实验观测得到的床面形态波长与ΓHGWk表现出一定关联性, 具体来说, 随着水流的ΓHG增大, 波长逐渐减小;随着Wk的增大, 波长缓慢增大。相比于形态参数与水流做功, 洪水的非恒定性对床面形态波长的影响更为明显。

      图  3  非恒定流形态、非恒定性及水流做功参数对床面形态波长(平均值+标准偏差)的影响

      Figure 3.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form wavelength (averaged value±standard deviation)

      图 4所点绘的床面形态波高与非恒定流参数关系表明, 随着洪水非恒定性及水流做功的增大, 波高也在逐渐增大(图 4(b), 图 4(c)),其中, 波高与洪水非恒定性具有显著正相关性。洪水形态与床面形态的波高之间无明显的响应关系。

      图  4  非恒定流形态、非恒定性及水流做功参数对床面形态波高(平均值+标准偏差)的影响

      Figure 4.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form height (averaged value±standard deviation)

      基于平均波长与波高可得到对应的床面形态的陡度, 其与非恒定流参数的关系由图 5呈现。显然, 床面形态的陡度随着水流非恒定性的增大而增大, 两者有显著相关性;陡度随着水流做功的增大而缓慢增大, 两者也具有相关性。而在洪水形态变化时, 床面形态的陡度并没有呈现趋势性变化, 两者之间无明显关联, 说明洪水形态对床面形态陡度的影响较微弱。

      图  5  非恒定流形态、非恒定性及水流做功参数对床面形态陡度(平均值+标准偏差)的影响

      Figure 5.  Impact of hydrograph shape, unsteadiness, and total water work on bed-form steepness (averaged value±standard deviation)

      值得注意的是, 不同尺度洪水作用下, 床面形态尺度与类型之间存在着良好的对应关系, 当波长与波高变化较小时, 床面形态也基本相似, 如不同形态洪水影响下形成的交错边滩(图 3(a), 图 4(a));当波长与波高变化显著时, 床面形态类型会发生相应改变, 如洪水非恒定性及水流做功变化而驱动形成的交错边滩、沙垄等床面形态, 见图 3(b)图 3(c)图 4(b)图 4(c)

    • 综上所述, 非恒定流基本参数(η, ΓHG, Wk)中, ΓHGWk对床面形态的尺度(波长、波高、陡度)有影响, 而洪水形态的影响微弱。对于床面形态尺度而言, 陡度是一个床面形态的综合描述参数, 包含了波长与波高, 且陡度是无因次量。因此, 有必要重点分析包含有非恒定性与水流做功的非恒定流综合判数与床面形态尺度的综合描述参数(即陡度)之间的关联性, 建立表征两者响应关系的定量表达式。

      由前文分析知, 非恒定性对床面形态尺度影响最大, 水流做功的影响次之。基于此, 可将2个参数的相对影响用一个综合判数来表示, 其形式为

      $$ \chi = {\mathit{\Gamma }_{{\rm{HG}}}}{W_{\rm{k}}}^\beta $$ (5)

      式中: β是待定参数, 其数值小于1.0。通过对数据的回归分析发现, 当β=0.4时, χ 与床面形态陡度呈现最佳拟合关系(图 6), 其相互关系为线性关系式, 即

      图  6  非恒定流综合参数ΓHGWk0.4与床面形态陡度的关系

      Figure 6.  Relationship between hydrograph combined parameter ΓHGWk0.4 with bed-form steepness ψ

      $$ \psi=3.90 \chi+0.0017 $$ (6)

      由拟合关系(6)知, 床面形态陡度作为反映床面形态的综合无因次量, 其与回归分析得到的非恒定流综合判数紧密关联, 拟合线性关系式时, 相关系数(R2)达到0.82。当然, 式(6)是基于本研究有限实验数据所得, 是否适用于更复杂的野外河流形态演变情形, 还需要更多实测数据进一步验证。

    • 本文运用非恒定性、水流做功、形态参数来量化多尺度非恒定流洪水过程, 研究非恒定流驱动下床面形态演变特征, 通过波长、波高及陡度来定量表征床面形态尺度。得到以下主要结论:

      (1) 在系列洪水作用下实验动床表面分别形成了交错边滩、沙垄及沙垄/边滩相混合的过渡型床面形态, 交错边滩常在非恒定性参数与水流做功均较大的洪水作用下形成, 而沙垄形成所对应的洪水条件正好相反。

      (2) 非恒定性参数对床面形态的陡度影响最大, 其次是水流做功参数, 而洪水形态的影响较为微弱;通过回归分析得到了非恒定流的综合判数, 其与床面形态陡度具有良好的关联性, 可建立表征非恒定流对床面形态陡度影响的定量描述关系。

参考文献 (22)

目录

    /

    返回文章
    返回