中国综合性科技类核心期刊(北大核心)

中国科学引文数据库来源期刊(CSCD)

美国《化学文摘》(CA)收录

美国《数学评论》(MR)收录

俄罗斯《文摘杂志》收录

留言板

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

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

崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟

鲁佩仪 朱建荣 钱伟伟 袁琳

鲁佩仪, 朱建荣, 钱伟伟, 袁琳. 崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
引用本文: 鲁佩仪, 朱建荣, 钱伟伟, 袁琳. 崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
LU Peiyi, ZHU Jianrong, QIAN Weiwei, YUAN Lin. Numerical simulation of erosion and deposition at the water intake channel of the outer seawall sluice in the ecological restoration project area of Chongming Dongtan Bird Habitat[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
Citation: LU Peiyi, ZHU Jianrong, QIAN Weiwei, YUAN Lin. Numerical simulation of erosion and deposition at the water intake channel of the outer seawall sluice in the ecological restoration project area of Chongming Dongtan Bird Habitat[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013

崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟

doi: 10.3969/j.issn.1000-5641.201941013
基金项目: 上海市科委重点项目(17DZ1201902)
详细信息
    通讯作者:

    朱建荣, 男, 教授, 博士生导师, 研究方向为河口海洋学. E-mail: jrzhu@sklec.ecnu.edu.cn

  • 中图分类号: P731.23

Numerical simulation of erosion and deposition at the water intake channel of the outer seawall sluice in the ecological restoration project area of Chongming Dongtan Bird Habitat

  • 摘要: 在河口水动力模式的基础上, 耦合泥沙模块和底部冲淤方程, 建立了崇明东滩鸟类保护区海堤外侧引水渠冲淤三维数值模式. 模式上游闸门水通量由堰流公式计算, 内侧水位由随塘河水动力模式计算, 模式下游水位、盐度、泥沙浓度由长江河口大区域数值模式计算. 利用2018年7月28—30日优化工程区水闸开闸冲淤前后实测引水渠底部高程进行验证, 模式计算的高程与实测值吻合良好, 表明模式能较好地模拟引水渠冲淤变化. 冲淤计算结果表明, 闸下引水渠在自然状态下, 涵闸关闭60 d后基本淤平, 闸口淤积量最大, 淤积程度自闸口向外海逐渐减小. 引水渠自然淤积2个月后, 随塘河蓄水至3.0 m开闸放水, 引水渠平均冲淤量为96 mm; 随塘河蓄水至3.8 m开闸放水, 引水渠平均冲淤量为133 mm; 纳潮冲淤一天, 引水渠平均冲淤量为625 mm. 纳潮冲淤是最为有效的冲淤方式. 研究成果可为崇明东滩国家级自然保护区生态修复工程区的保护与管理、引水渠冲淤方案提供科技指导.
  • 图  1  长江口崇明东滩国家级自然保护区优化工程区C区位置(左)和布局(右)

    注: 白线范围之内为生态修复工程区

    Fig.  1  Location (left) and layout (right) of Chongming Dongtan National Nature Reserve

    图  2  闸外引水渠

    Fig.  2  Water intake channel outside the sluice

    图  3  闸下引水渠模式网格

    Fig.  3  Model grid of the water intake channel

    图  4  随塘河模式网格

    Fig.  4  Model grid of Suitang River

    图  5  长江河口模型计算区域和网格(左)及放大的崇明东滩区域网格(右)

    Fig.  5  Model domain and grid of Changjiang Estuary (left), enlarged view of the grid near Chongming Eastern Shoal (right)

    图  6  7月引水渠东边界水位(a)、盐度(b)和悬浮泥沙浓度(c)随时间变化

    Fig.  6  Temporal variation in water elevation (a), salinity (b), and suspended sediment concentration (c) outside the sluice in July

    图  7  引水渠开闸冲淤模式验证

    注: 蓝线为冲淤前测量高程, 红点为冲淤后测量高程, 黑线为模型计算冲淤后高程; 横坐标为与闸门间的距离; 基面为吴淞基面, 下同

    Fig.  7  Model validation for scouring and deposit at the water intake channel

    图  8  闸下引水渠自然状态下2个月后淤积厚度

    Fig.  8  Deposit thickness at the water intake channel after 2 months under natural conditions

    图  9  随塘河蓄水至3.0 m放水至2.0 m工况底床冲刷厚度分布(左)和底高程分布(右)

    Fig.  9  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.0 m of water and discharging down to 2.0 m

    图  10  随塘河蓄水至3.8 m放水至2.0 m工况底床冲刷厚度分布(左)和底高程分布(右)

    Fig.  10  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.8 m of water and discharging down to 2.0 m

    图  11  随塘河蓄水至3.0 m纳潮冲淤工况底床冲刷厚度分布(左)和底高程分布(右)

    Fig.  11  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.0 m of water and with the sluice opened to receive tidal water

  • [1] 杨元平, 唐子文, 史英标. 曹娥江大闸闸下淤积现状分析及减淤对策初探 [J]. 浙江水利科技, 2013(2): 35-39. DOI:  10.3969/j.issn.1008-701X.2013.02.008.
    [2] 金元欢, 沈焕庭. 我国建闸河口冲淤特性 [J]. 泥沙研究, 1991(4): 59-68.
    [3] 李大山, 任汝述. 波浪作用下异重流运动特性研究 [J]. 海洋工程, 1996(4): 54-59.
    [4] 窦国仁. 窦国仁论文集 [M]. 北京: 中国水利水电出版社, 2003.
    [5] 伍冬领, 林炳尧, 余大进, 等. 永宁江枢纽引航道和挡潮闸闸港冲淤研究 [J]. 泥沙研究, 2003(2): 69-72. DOI:  10.3321/j.issn:0468-155X.2003.02.012.
    [6] 曲红岭, 窦希萍, 赵晓冬, 等. 不同类型建闸河口闸下淤积模拟计算 [J]. 海洋工程, 2011, 29(2): 59-67. DOI:  10.3969/j.issn.1005-9865.2011.02.009.
    [7] 成建, 刘德荣. 浅析闸下港道淤积及治理问题 [J]. 山东水利, 2011(7): 24-26. DOI:  10.3969/j.issn.1009-6159.2011.07.011.
    [8] 徐雪松, 窦希萍, 陈星, 等. 建闸河口闸下淤积问题研究综述 [J]. 水运工程, 2012(1): 116-121. DOI:  10.3969/j.issn.1002-4972.2012.01.024.
    [9] 王文杰. 海泊河挡潮闸泥沙模型试验研究 [D]. 山东 青岛: 中国海洋大学, 2012.
    [10] 朱明成. 江苏淤涨型海岸闸下淤积问题初探 [J]. 江苏水利, 2014(6): 11-12.
    [11] 钱建平, 严琴琴, 夏栩. 小洋口闸下港槽淤积分析 [J]. 中国水运, 2014, 14(11): 307-309.
    [12] 陈其美. 新洋港闸下港道淤积分析及对策 [J]. 水利建设与管理, 2009, 29(9): 57-59. DOI:  10.3969/j.issn.1005-4774.2009.09.024.
    [13] 孙振华, 虞快. 崇明东滩鸟类自然保护区的建立及其意义 [J]. 上海建设科技, 1998(4): 24-26.
    [14] 孙鸿雁, 唐芳林, 周红斌, 等. 上海崇明东滩鸟类国家级自然保护区总体规划的启示 [J]. 林业建设, 2013(1): 22-26.
    [15] BLUMBERG A F, MELLOR G L. A description of a three-dimensional coastal ocean circulation model [G]// HEAPS N S. Three-Dimensional Coastal Ocean Models. Washington, DC: American Geophysical Union, 1987: 1-16.
    [16] QIU C, ZHU J R. Influence of seasonal runoff regulation by the Three Gorges Reservoir on saltwater intrusion in the Changjiang River Estuary [J]. Continental Shelf Research, 2013, 71: 16-26. DOI:  10.1016/j.csr.2013.09.024.
    [17] LYU H H, ZHU J R. Impact of the bottom drag coefficient on saltwater intrusion in the extremely shallow estuary [J]. Journal of Hydrology, 2018, 557: 838-850. DOI:  10.1016/j.jhydrol.2018.01.010.
    [18] LI L, ZHU J R, WU H. Impacts of wind stress on saltwater intrusion in the Yangtz Estuary [J]. Science China Earth Sciences, 2012, 55(7): 1179-1192.
    [19] LI L, ZHU J R, WU H, et al. A numerical study on the water diversion ratio of the Changjiang Estuary during the dry season [J]. Chinese Journal of Oceanology and Limnology, 2010, 28(3): 700-712. DOI:  10.1007/s00343-010-9114-2.
    [20] 朱建荣, 吴辉, 顾玉亮. 长江河口北支倒灌盐通量数值分析 [J]. 海洋学研究, 2011, 29(3): 1-7. DOI:  10.3969/j.issn.1001-909X.2011.03.002.
    [21] QIU C, ZHU J R. Assessing the influence of sea level rise on salt transport processes and estuarine circulation in the Changjiang River Estuary [J]. Journal of Coastal Research, 2015, 31(3): 661-670.
    [22] 刘高峰. 长江口水沙运动及三维泥沙模型研究 [D]. 上海: 华东师范大学, 2011.
    [23] 曹祖德, 王运洪. 水动力泥沙数值模拟 [M]. 天津: 天津大学出版社, 1994: 241-249.
    [24] 罗志发. 长江口及其邻近海域泥沙输运及其动力机制 [D]. 上海: 华东师范大学, 2018.
    [25] LIANG B C, LI H J, LEE D Y. Bottom shear stress under wave-current interaction [J]. Journal of Hydrodynamics Ser B, 2008, 20(1): 88-95. DOI:  10.1016/S1001-6058(08)60032-3.
    [26] 吴宇帆. 长江河口细颗粒泥沙沉降速度研究 [D]. 上海: 华东师范大学, 2016.
    [27] 刘红, 何青, 吉晓强, 等. 波流共同作用下潮滩剖面沉积物和地貌分异规律——以长江口崇明东滩为例 [J]. 沉积学报, 2008(5): 833-843.
    [28] 戴苒, 朱建荣. 长江口崇明东滩风况统计分析 [J]. 华东师范大学学报(自然科学版), 2015(4): 17-25.
    [29] 庞启秀, 白玉川, 杨华, 等. 淤泥质浅滩泥沙临界起动切应力剖面确定 [J]. 水科学进展, 2012, 23(2): 249-255.
    [30] 赵振兴, 何建京. 水力学 [M]. 2版. 北京: 清华大学出版社, 2010: 176.
    [31] 戴苒. 长江河口河床泥沙临界起动应力空间分布 [D]. 上海: 华东师范大学, 2016.
  • [1] 钱丽洁, 杨岩, 孙真荣.  反向旋转圆偏振双色激光场驱动下氩原子非顺序双电离研究 . 华东师范大学学报(自然科学版), 2020, (2): 55-63. doi: 10.3969/j.issn.1000-5641.201922010
    [2] 陈洁, 吴健平, 郭珮珺, 姚申君.  “保基”原则下上海市养老机构服务可达性评价 . 华东师范大学学报(自然科学版), 2020, (3): 119-128. doi: 10.3969/j.issn.1000-5641.201841038
    [3] 苏爱平, 吕行行, 吴宇帆.  南水北调工程对长江口盐水入侵和淡水资源的影响 . 华东师范大学学报(自然科学版), 2020, (3): 32-42. doi: 10.3969/j.issn.1000-5641.201941031
    [4] 朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行.  长江河口北支建闸对减轻盐水入侵的数值模拟 . 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
  • 加载中
图(11)
计量
  • 文章访问数:  24
  • HTML全文浏览量:  32
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-25
  • 网络出版日期:  2020-05-29
  • 刊出日期:  2020-05-01

崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟

doi: 10.3969/j.issn.1000-5641.201941013
    基金项目:  上海市科委重点项目(17DZ1201902)
    通讯作者: 朱建荣, 男, 教授, 博士生导师, 研究方向为河口海洋学. E-mail: jrzhu@sklec.ecnu.edu.cn
  • 中图分类号: P731.23

摘要: 在河口水动力模式的基础上, 耦合泥沙模块和底部冲淤方程, 建立了崇明东滩鸟类保护区海堤外侧引水渠冲淤三维数值模式. 模式上游闸门水通量由堰流公式计算, 内侧水位由随塘河水动力模式计算, 模式下游水位、盐度、泥沙浓度由长江河口大区域数值模式计算. 利用2018年7月28—30日优化工程区水闸开闸冲淤前后实测引水渠底部高程进行验证, 模式计算的高程与实测值吻合良好, 表明模式能较好地模拟引水渠冲淤变化. 冲淤计算结果表明, 闸下引水渠在自然状态下, 涵闸关闭60 d后基本淤平, 闸口淤积量最大, 淤积程度自闸口向外海逐渐减小. 引水渠自然淤积2个月后, 随塘河蓄水至3.0 m开闸放水, 引水渠平均冲淤量为96 mm; 随塘河蓄水至3.8 m开闸放水, 引水渠平均冲淤量为133 mm; 纳潮冲淤一天, 引水渠平均冲淤量为625 mm. 纳潮冲淤是最为有效的冲淤方式. 研究成果可为崇明东滩国家级自然保护区生态修复工程区的保护与管理、引水渠冲淤方案提供科技指导.

English Abstract

鲁佩仪, 朱建荣, 钱伟伟, 袁琳. 崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
引用本文: 鲁佩仪, 朱建荣, 钱伟伟, 袁琳. 崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
LU Peiyi, ZHU Jianrong, QIAN Weiwei, YUAN Lin. Numerical simulation of erosion and deposition at the water intake channel of the outer seawall sluice in the ecological restoration project area of Chongming Dongtan Bird Habitat[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
Citation: LU Peiyi, ZHU Jianrong, QIAN Weiwei, YUAN Lin. Numerical simulation of erosion and deposition at the water intake channel of the outer seawall sluice in the ecological restoration project area of Chongming Dongtan Bird Habitat[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
    • 在我国沿海大部分中小入海河口建闸以挡潮御卤、排洪蓄淡, 以及防台抗暴, 但河口建闸面临的一个重要问题是闸下淤积. 我国建闸河口闸下淤积非常严重, 显著影响上游河道的防洪排涝, 闸下淤积预报和减淤措施研究一直是水利工程中急需解决的问题之一. 在以海相来沙为主的潮汐河口建闸, 闸下淤积问题相当普遍. 海河流域35个河口建闸后, 22座发生严重淤积. 江苏沿海修建的58座大中型挡潮闸中, 5座基本淤死. 浙江省甬江支流姚江建闸后闸下河段及甬江河床发生严重淤积[1]. 金元欢和沈焕庭指出, 建闸河口淤积的主要原因是水动力条件发生了改变, 淤泥质河口与砂质河口冲淤特性差异很大, 大部分闸门第一年淤积速度最大[2]. 在泥沙含量较高的水域, 由于密度差的作用, 浑水会流向清水形成异重流. 李大山、任汝述基于对射阳河建闸前后径流量、潮动力条件等水力要素的比对, 对该地区泥沙来源和泥沙输运量进行分析[3]. 窦国仁得出射阳河发生闸下淤积的条件和影响因素[4]. 伍冬领等对永宁江枢纽引航道和挡潮闸下淤积等问题进行了深入研究, 认为该处闸港泥沙粒径较细, 易发生絮凝[5]. 曲红岭等采用数学模型研究了短引河河道型、短引河滩槽型和长引河河道型三种类型建闸河口闸下淤积特征, 结果表明闸下淤积与闸下河道特征密切相关[6]. 河口建闸后, 由于潮波变形及径流量的匮乏, 浮泥随纳潮水体进入闸下河道形成泥沙淤积. 成建、刘德荣结合江苏沿海一些挡潮闸, 全面分析了闸下港道淤积泥沙来源和泥沙淤积的原因, 比较全面地介绍了一些可行的治理措施[7]. 徐雪松等总结影响闸下淤积形态的多个因素, 包括潮流作用、闸下引河长度、河口泥沙的动力特性、闸下港道类型以及建闸位置等[8]. 王文杰认为闸下河道发生持续性淤积是由水流挟沙能力和含沙量的不对称所导致[9]. 杨元平等研究得出曹娥江大闸闸下淤积现状及减淤对策, 认为尖山河段主槽靠近曹娥江口, 曹娥江大闸闸下淤积不明显; 一旦出现主槽偏离南岸, 闸下滩面较宽时, 闸下淤积问题就会出现[1]. 朱明成分析了闸下淤积的影响因素, 论证了闸址条件不同时闸下淤积情况, 其中闸下河道类型、闸下河道长度、径流条件、植被条件和闸下河道入海方向都起到一定作用[10]. 钱建平等分析了小洋口闸下港槽淤积成因, 提出了减淤的措施[11].

      以往研究提出的闸下减淤的对策有: (1)水力冲淤. 将势能转化为动能, 利用水流的携沙能力达到清淤的效果. 水位清淤适用于闸门上下河道比降大且水源较为充足的河流, 开闸放水的时机也要考虑到潮汐能量的变化. 纳潮冲淤指在涨潮期利用河道纳潮, 待有一定水头差时开闸冲淤, 一般在上游水源不足或潮汐能量强的地区可考虑该方法[12]. 或者, 大潮期间打开闸门, 增强纳潮量和落潮流速, 达到清淤效果. (2)机械清淤. 机船拖耙清淤, 是以机动船拖动耙具耙起淤泥, 利用推进器, 搅动水体, 使部分耙起的淤泥变成悬沙, 随水流流走. 水力挖塘机组, 抓斗式或链斗式挖泥船亦可用于引河清淤. 也可采用绞吸式、抓斗式或链斗式挖泥船清淤. (3)工程措施. 整治闸上河道和闸下引河, 扩大排泄量; 闸下裁弯、切滩, 增加泄流能力; 以清顶浑, 对于较短的港槽可以在涨潮前2~3 h, 在引河中预先开闸放水, 减少进入引河的浑水. 选好闸址, 应考虑冲淤水源、港口的稳定性、闸下引河的长短等条件; 围垦时保留冲港滩面水; 植物拦沙减淤, 阻止滩面产沙、进入引河; 修筑导堤, 束水攻沙, 切断由风浪掀起的海滩泥沙补给源, 减少泥沙进入引河.

      长江口挟沙量巨大, 崇明东滩鸟类保护区新建大堤外区域的淤涨明显, 海堤水闸外侧引水渠极易淤积. 为了提高冲淤效果, 降低开闸能耗, 保障调控优化工程区生态用水安全, 就需要了解引水渠自然淤积情况, 以及开闸放水冲淤的效果. 本文对此开展研究, 给出在水闸关闭下引水渠自然淤平所需时间、开闸冲淤和纳潮冲淤效果, 研究成果可为崇明东滩保护区大型湿地修复工程水闸运行提供关键技术支撑, 对崇明东滩鸟类保护区湿地资源保育具有重要的理论意义和应用价值.

    • 崇明东滩位于长江口崇明岛东端, 是东亚-澳大利亚候鸟迁飞路线的中段[13]. 这里植被分布广, 有大量的底栖动物, 适宜水鸟停歇、越冬, 是候鸟迁徙的中转站. 1998年上海市政府成立崇明东滩鸟类保护区, 2005年成为国家级自然保护区[14]. 上海市政府以遏制互花米草扩张, 改善崇明东滩生态系统, 优化鸟类栖息地, 稳定鸟类食物来源为目标, 在崇明东滩鸟类国家级自然保护区内开展了规模巨大的互花米草生态控制及鸟类栖息地优化工程(以下简称优化工程). 该优化工程北边在八滧闸门, 南边接崇明东滩1998大堤中部, 西边以崇明东滩1998大堤为界, 东边界为新建的海塘大堤. 优化工程范围达240 km2, 于2016年12月竣工. 本文研究区域为优化工程区C区 (见图1).

      图  1  长江口崇明东滩国家级自然保护区优化工程区C区位置(左)和布局(右)

      Figure 1.  Location (left) and layout (right) of Chongming Dongtan National Nature Reserve

      该优化工程包括内部生境区、随塘河、大堤和涵闸. 生态修复工程区的南部C区, 内部生境区分为C1N、C1S、C2W、C2E、C3、C4六个不同的生境区块(见图1). 随塘河围绕内部生境区域一圈, 是内部生境与外界进行水交换的重要水域. 由于蒸发、降水、植物蒸散和土壤下渗, 生境区的水深随时间变化. 当水深不在候鸟合适的水深范围内, 就需要保护区管理人员从随塘河引水或排水, 以维持不同内部生境区适宜水深. 在调控优化工程区新建海堤上建有水闸, 当内部生境区和随塘河水位低、不能满足生态用水时, 需要开闸引进闸外崇明东滩水体; 当降水量大, 如梅雨季节, 内部生境区和随塘河水位超过临界值时, 需要开闸放水. 在一年中的大部分时间, 由于降水的补充和生态耗水, 内部生境区和随塘河的水位一般能满足要求, 大部分时间水闸是关闭的. 长时间的水闸关闭, 以及崇明东滩高泥沙浓度环境, 极易导致闸外引水渠淤积. 引水渠淤积后就难以开闸取水, 需要开闸冲淤. 为了提高闸下冲淤效率, 优化工程区在闸口外建立了宽50 m、长300 m的引水渠(见图2). 引水渠两岸是由挡水建筑物建成的导堤, 在高潮位时也不会被淹没. 闸门开启期间, 当随塘河水位低于外海时, 水体由引水渠自外海引入; 当随塘河水位高于外海时, 水体自引水渠向外海流出. 水流受到引水渠导堤的约束, 动能相对集中, 作用在底床泥沙的时间增长, 底切应力增大, 可提高冲刷效率.

      图  2  闸外引水渠

      Figure 2.  Water intake channel outside the sluice

    • 水动力模式采用作者课题组长期改进和应用的三维数值模式ECOM-si (Semi-implicit Estuarine, Costal and Ocean Model). 该模式是基于海洋数值模式POM (Princeton Ocean Model) 的基础上研发的[15], 可设置包括潮流、径流、风、波浪等各种动力因子. 改进的ECOM-si模式在长江口水动力过程和盐水入侵等方面研究中得到广泛的应用, 取得了大量研究成果[16-21]. 为了节省篇幅, 有关模式介绍可参考已发表的文献[16-18].

    • 三维泥沙数值模式是在上述水动力模式基础上耦合泥沙模块建立起来的. 本文中泥沙模块运用是基于刘高峰改进后的泥沙模块, 该模式考虑了泥沙的沉降过程和底部临界启动应力, 具体模块介绍和参数设置参照文献[22], 这里简单介绍泥沙控制方程. 悬沙输运模块控制方程为

      $$\frac{{\partial DJ{S_c}}}{{\partial t}} + \frac{{\partial DJ\hat U{S_c}}}{{\partial \xi }} + \frac{{\partial DJ\hat V{S_c}}}{{\partial \eta }} + \frac{{\partial J(\omega - {w_s}){S_c}}}{{\partial \sigma }} = \frac{1}{D}\cdot\frac{\partial }{{\partial \sigma }}\left( {{K_v}\frac{{\partial J{S_c}}}{{\partial \sigma }}} \right) + DJ{F_{sc}}.$$ (1)

      式中, $\xi $, $\eta $为网络2个水平坐标, σ为网络垂向坐标, D为水深, J为网络面积, Sc为悬浮泥沙浓度, ω为水体垂向流速, ωs为泥沙沉降速度, Kv为垂向湍流扩散系数, Fsc为泥沙水平湍流扩散项, $\widehat U $$\widehat V $分别为与$\eta $$\xi $正交的法向速度.

      模型中泥沙通量的海表边界条件为

      $$ - {w_s}{S_c} - \left( {\frac{{{K_v}}}{D}\cdot\frac{{\partial {S_c}}}{{\partial \sigma }}} \right) = 0,\;\;\;\;{\rm{\sigma }} = 0. $$ (2)

      海底泥沙通量的边界条件为

      $$ - {w_s}{S_c} - \left( {\frac{{{K_v}}}{D}\cdot\frac{{\partial {S_c}}}{{\partial \sigma }}} \right) = F,\;\;\;\;\;{\rm{\sigma }} = - 1. $$ (3)

      其中, $F = {q_{ero}} - {q_{dep}}$, F为水底净泥沙通量, ${q_{ero}}$是水底泥沙的冲刷通量, ${q_{dep}}$表示水底泥沙的沉积通量, 计算公式为[23]

      $$ {q_{ero}} = \left\{ {\begin{array}{*{20}{l}} {0,\quad\quad\quad\quad\;\;\;\tau < {\tau _e}},\\ {M\left( {\frac{\tau }{{{\tau _e}}} - 1} \right),\tau \geqslant {\tau _e}}, \end{array}} \right. $$ (4)
      $$ {q_{dep}} = \left\{ {\begin{array}{*{20}{l}} {0,\quad\quad\quad\quad\quad\quad\;\;\;\tau > {\tau _d}},\\ {\alpha '{w_s}{S_c}\left( {1 - \frac{{\tau }}{{{\tau _d}}}} \right),\tau \leqslant {\tau _d}}. \end{array}} \right. $$ (5)

      (4)式和(5)式中, $\tau $为模式计算所得的底部切应力, ${\tau _e}$为临界起动应力, ${\tau _d}$为临界沉降应力. M为底床冲刷系数, 表示了水底冲刷的难易程度, 取值范围为1 × 10–5~4 × 10–3 kg/(m2·s) . $\alpha '$为泥沙发生沉降的概率, 根据经验一般取值为0.67~0.84[24].

      波浪与潮流共同作用下的底切应力参考Liang and Li[25]中给出的计算公式. 泥沙沉速使用的是吴宇帆[26]经过多次沉降实验建立的絮凝沉降公式. 悬沙粒径的设定决定了静水沉速的计算结果, 对于三维水动力泥沙数值模式的模拟至关重要. 长江口及邻近海域悬沙组分多样, 在平面上和垂向上悬沙粒径变化明显, 需依靠实测数据确定该海域悬沙粒径. 刘红等在对崇明东滩悬浮泥沙进行采集分析后指出, 崇明东滩护花米草和泥滩表层沉积物较细, 悬浮泥沙中值粒径为8~16 μm, 自岸向海逐渐变粗, 模式中泥沙中值粒径统一取15 μm[27]. 详细的长江口泥沙絮凝沉降公式可参见文献[26].

      计算底部泥沙通量的关键参数是临界起动应力, 临界起动应力受到多种因素影响, 难以使用具体公式计算, 通过实测含沙量进行调试率定, 设置经验值是目前最有效的计算方法. 本文根据临界起动应力, 利用经验公式来换算临界沉降应力, 即${\tau _d} = 4/9{\tau _e}$[23].

      河床底部高程的变化直接受到底部泥沙净通量的影响, 泥沙净通量为负时, 河床发生冲刷, 河床底高程减小; 泥沙净通量为正时, 发生淤积, 河床底高程增大. 引水渠内局地高程变化$\Delta h $的计算公式为

      $$ \Delta h = \left( {F \cdot t} \right)/\rho . $$ (6)

      其中, F为水底净泥沙通量, 是床底泥沙的侵蚀通量和沉积通量的差值, $F = {q_{ero}} - {q_{dep}}$; t为时间 (s); ρ为表层沉积物的密度 (kg/m3).

    • 新建大堤水闸外侧引水渠淤积和冲刷是本文研究的要点. 引水渠西边接涵闸闸口, 东边直通崇明东滩, 导堤两边为优化工程建设挡水的堤坝, 最高点为5.0 m, 外海潮水基本不会漫过导堤, 因此可以将其作为固边界设定(见图2). 使用Delft3D网格生成软件设计的闸下引水渠正交化网格精度为2 m(见图3).

      图  3  闸下引水渠模式网格

      Figure 3.  Model grid of the water intake channel

      水表面风应力由长江口崇明东滩实测月平均风速和风向给出[28]. 刘红在对崇明东滩表层沉积物进行采集分析后指出, 崇明东滩护花米草和泥滩表层沉积物粒径由岸向海逐渐增大, 靠近外海处的中值粒径为33~84 μm, 本文取80 μm. 本文依据淤泥质海滩, 泥沙临界起动应力取值为0.05 N/m2[29].

      闸下引水渠侧向开边界为西侧的闸门和东侧的崇明东滩. 水闸闭合时为固边界, 打开时为开边界. 因闸门内外水位不连续和网格不同, 过闸水通量由堰流公式计算. 引水渠东侧开边界由长江口大区域模式计算所得的水位、盐度和泥沙浓度提供.

    • 建立随塘河模式是为了计算闸内随时间变化的水位, 采用堰流公式计算可靠的过闸水通量. 随塘河的网格见图4, 网格分辨率为10 m. 关闸时, 随塘河的流动主要由风应力驱动; 开闸时, 受过闸流量和风应力的共同作用.

      图  4  随塘河模式网格

      Figure 4.  Model grid of Suitang River

      依据宽顶堰闸孔淹没出流的计算公式计算过闸水通量[30], 即

      $$ {\rm{Q}} = {\rm{\varepsilon '}}\varphi eB\sqrt {2g\left( {{H_0} - {H_t}} \right)} . $$ (7)

      其中, ε′ 是闸门的平面收缩系数, 查表取值为0.661; φ是闸孔流速系数, 取值范围为0.85~0.95, 这里取0.95; e为闸门开度, 是闸门底部距离坎顶的铅直距离; B为闸门总净宽, 该处涵闸为五孔闸, 每孔宽度为3 m, 总净宽为15 m; H0为上游水位; Ht为下游水位(当外海水位高于随塘河时, 外海水位为上游水位, 否则为下游水位).

      由闸孔出流流量计算公式可知, 闸门水通量的计算同时受到闸内水位与闸外水位的制约. 闸内水位由随塘河模式计算, 闸外水位由长江口大区域模式计算.

    • 长江口大区域数值模式为引水渠模式提供水位、盐度和泥沙浓度开边界条件. 模式的网格设置见图5. 模式中外海开边界条件、风场、泥沙沉降速度、底部泥沙临界起动应力的设置和给出, 以及波浪场的计算等, 可参见罗志发博士的学位论文[24]. 模式的计算区域和网格见图5, 包括整个长江河口、杭州湾和邻近海域.

      图  5  长江河口模型计算区域和网格(左)及放大的崇明东滩区域网格(右)

      Figure 5.  Model domain and grid of Changjiang Estuary (left), enlarged view of the grid near Chongming Eastern Shoal (right)

      随塘河闸外海域7月闸外水位、盐度模拟计算结果过程见图6. 由于崇明东滩地势较高, 高潮位时淹没, 低潮位时露出, 此时水位即为滩地高程, 随时间不变, 盐度不连续. 随塘河三维泥沙模式采用长江口大区域三维泥沙模式计算的悬浮泥沙浓度作为泥沙边界条件, 7月随塘河闸外泥沙浓度见图6.

    • 本文三维水动力和泥沙数值模式已做了大量验证, 结果表明此模式能很好地模拟水文、流速、流向和盐度的变化, 也能较好地模拟泥沙浓度的变化. 为了节省篇幅, 有关模式验证等可参考文献[16-19, 24, 26, 31].

      本文采用2018年7月引水渠开闸前后底高程观测资料验证模式. 崇明东滩候鸟保护区优化工程区的涵闸开启调度情况为: 2018年7月27日21时开启涵闸, 从外海引水, 7月28日1时关闭涵闸; 7月28日21时30分开启涵闸, 从外海引水, 7月29日2时关闭涵闸; 7月29日22时开启涵闸, 从外海引水, 7月30日1时20分关闭涵闸. 闸门开启时, 过闸水流不断冲刷闸下引水渠. 采用三维红外激光扫描仪于7月27日开闸前测量闸下引水渠断面底高程, 于7月30日关闸后测量闸下引水渠断面底高程. 通过反复调整模型中底部泥沙临界启动应力, 确定最佳值. 实测高程与模型计算数据的比较见图7. 由实测数据可得引水渠内平均冲淤122 mm, 模型计算值为123 mm, RMSE(Root Mean Square Error)、CC(相关系数 Correlation Coefficent)、SS值 (Skill Score) 分别为0.062 5 m、0.74、0.851 1, 验证结果良好.

      图  6  7月引水渠东边界水位(a)、盐度(b)和悬浮泥沙浓度(c)随时间变化

      Figure 6.  Temporal variation in water elevation (a), salinity (b), and suspended sediment concentration (c) outside the sluice in July

      图  7  引水渠开闸冲淤模式验证

      Figure 7.  Model validation for scouring and deposit at the water intake channel

    • 在自然状态下, 引水渠西侧闸门关闭, 只有东侧与崇明东滩水域相连, 随涨落潮交换水体、盐分和泥沙. 模式设置计算时段为1月1日至3月1日, 引水渠东侧开边界由水位、盐度和泥沙浓度给出, 资料来自大区域模式计算的结果. 引水渠初始设置为无泥沙淤积, 底高程为1.0 m.

      模式计算的闸下淤积厚度的分布见图8. 涵闸关闭两个月后, 引水渠靠近闸门处基本淤平, 引水渠内淤积严重. 涵闸关闭时, 引水渠内受到外海潮涨潮落的影响, 高潮位时长江口携带大量悬浮泥沙的盐水进入闸下引水渠, 进入引水渠后水体动能逐渐减小, 水流携沙能力逐渐减小, 到达闸口时水流流速最小, 底切应力减小, 泥沙发生沉降, 难以再次起悬, 此时最易发生淤积. 闸下引水渠靠近外海侧受到外海潮动力的影响, 水动力泥沙易悬浮, 较难沉降, 此处淤积程度不高.

      图  8  闸下引水渠自然状态下2个月后淤积厚度

      Figure 8.  Deposit thickness at the water intake channel after 2 months under natural conditions

      由关闸后引水渠中的底高程变化可以看出, 引水渠闸口淤积量最大, 靠近外海淤积量较小, 引水渠内淤积程度自闸口向外海逐渐减小, 同时靠近西北侧一岸淤积量较大. 闸前水流流速很小, 几乎为零, 在异重流的作用下, 闸门附近的底部泥沙向闸下输运, 闸下底部切应力很小, 泥沙无法被落潮的水流带走, 基本全部沉降下来, 因此闸前淤积程度最高. 在科里奥利惯性力(科氏力)的作用下, 涨潮时, 进入引水渠内的水体向北偏移, 北部水体挟沙能力较大; 落潮时, 离开引水渠的水流向南偏移, 南部水体挟沙能力较大, 泥沙冲刷量较大. 引水渠南北两岸挟沙能力的不平衡导致了引水渠内淤积厚度不同.

    • 闸门淤积厚度较低时, 冲刷效果不明显, 频繁开启闸门进行冲淤耗时耗力, 不利于优化工程区的管理; 闸门淤积程度高时, 长时间不开启闸门则很有可能淤死. 本节分别选择自然淤积两个月后的淤积程度进行冲淤效果模拟, 分为3种工况进行: (1) 随塘河蓄水至3.0 m, 开闸放水至2.0 m; (2) 随塘河蓄水至3.8 m, 开闸放水至2.0 m; (3) 随塘河蓄水至3.0 m, 选择大潮时期开闸一天. 本节分别模拟和分析自然淤积两个月情况下3种不同工况的冲淤效果.

      模式初始底高程为淤积两个月后的引水渠内底高程, 冲淤前引水渠底高程图见图8. 随塘河蓄水至3.0 m, 放水至2.0 m, 开闸放水持续时间长达20 h后的冲淤效果见图9, 引水渠内平均总冲淤厚度为96 mm, 闸前冲刷效果最好. 这是因为靠近闸口处流态复杂, 流线急剧弯曲, 湍流动能大, 底床冲刷效果明显. 由于淤积时间长、引水渠内淤积程度高, 在相同的水力条件下, 水深较小, 底部切应力较大, 冲淤效果明显.

      随塘河蓄水至3.8 m, 放水至2.0 m, 开闸放水时长达22 h后的冲淤效果见图10. 引水渠闸前冲刷效果较好, 渠内平均总冲淤厚度为133 mm.

      随塘河蓄水至3.0 m, 外海高水位进水, 开闸一天后的纳潮冲淤效果见图11. 纳潮冲淤时, 引水渠内受到外海潮动力的作用, 引水渠内整体冲淤效果较好, 渠内平均总冲淤厚度为625 mm, 闸前冲淤效果比较好. 在冲淤的模式验证中, 3 d的平均冲淤幅度为122 mm, 远比这里的平均冲淤量值小. 原因是: (1)模式验证时涵闸只在夜间涨潮时开启, 3 d的开闸时间总和约为12 h, 而这里开闸时长为24 h, 少了一半. (2)验证时段前期, 台风安比2018年7月22日12时30分在上海崇明登陆, 大量降雨使得保护区水闸打开排涝, 冲刷了引水渠. 模式验证时引水渠初始淤积高程较低, 导致模式验证时总的冲淤幅度较小. (3)这里纳潮冲淤选择的外边界水位是7月最高潮位的一天, 比验证计算中的外边界水位更高.

      图  9  随塘河蓄水至3.0 m放水至2.0 m工况底床冲刷厚度分布(左)和底高程分布(右)

      Figure 9.  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.0 m of water and discharging down to 2.0 m

      图  10  随塘河蓄水至3.8 m放水至2.0 m工况底床冲刷厚度分布(左)和底高程分布(右)

      Figure 10.  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.8 m of water and discharging down to 2.0 m

      图  11  随塘河蓄水至3.0 m纳潮冲淤工况底床冲刷厚度分布(左)和底高程分布(右)

      Figure 11.  Distribution of bed scouring thickness (left) and bed height (right) at the water intake channel in the case of accumulating up to 3.0 m of water and with the sluice opened to receive tidal water

    • 本文在河口水动力模式的基础上, 耦合泥沙模块和底部冲淤方程, 建立了崇明东滩鸟类保护区海堤外侧引水渠冲淤三维数值模式. 模式上游闸门水通量由堰流公式计算, 内侧水位由随塘河水动力模式计算, 模式下游水位、盐度、泥沙浓度由长江河口大区域数值模式计算. 利用2018年7月28—30日优化工程区水闸开闸冲淤前后实测引水渠底部高程, 模式计算的高程与实测值吻合良好, 表明模式能较好地模拟引水渠冲淤变化.

      模式计算结果表明, 引水渠在自然状态下关闭60 d后闸口处基本淤平, 无法使用. 涵闸闸前淤积厚度约达2.0 m, 引水渠靠近外海处淤积厚度约为1.0 m, 需要开闸冲淤. 淤积厚度较高时, 冲淤效果明显. 仅在随塘河水力冲淤、无外海潮动力作用时, 闸前冲淤效果明显, 靠近外海处冲淤厚度几乎为0.

      引水渠自然淤积两个月后, 随塘河蓄水至3.0 m、放水至2.0 m, 引水渠内平均冲淤厚度为96 mm; 随塘河蓄水至3.8 m、放水至2.0 m, 引水渠内平均冲淤厚度为133 mm. 引水渠内淤积程度高, 开闸冲淤时, 流量相同的情况下, 水深较浅, 流速较大, 冲淤量也因此增加. 纳潮冲淤后, 引水渠内整体冲淤效果好, 平均冲淤厚度为625 mm.

      分析比较不同工况下引水渠内的冲淤厚度, 淤积程度越高, 冲淤效果越明显. 利用随塘河蓄水冲淤时, 闸前冲淤效果好, 总体冲淤厚度较小, 水位下降快, 水量损失较大, 不利于维护优化工程区闸内生态系统的稳定. 在高潮位纳潮冲淤时, 水动力强, 闸下引水渠整体冲淤效果明显, 引水渠内淤积的泥沙基本清除. 此时随塘河内水量损失较少, 有利于优化工程区内鸟类栖息地的建设与修复. 本文研究成果可为崇明东滩国家级自然保护区生态修复工程区的保护与管理、引水渠冲淤方案提供科技指导.

参考文献 (31)

目录

    /

    返回文章
    返回