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

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

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

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

俄罗斯《文摘杂志》收录

留言板

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

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

长江河口北支建闸对减轻盐水入侵的数值模拟

朱建荣 鲁佩仪 唐川敏 陈晴 吕行行

朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行. 长江河口北支建闸对减轻盐水入侵的数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
引用本文: 朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行. 长江河口北支建闸对减轻盐水入侵的数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
ZHU Jianrong, LU Peiyi, TANG Chuanmin, CHEN Qing, LÜ Hanghang. Numerical simulation of saltwater intrusion mitigation by building a sluice in the North Branch of the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
Citation: ZHU Jianrong, LU Peiyi, TANG Chuanmin, CHEN Qing, LÜ Hanghang. Numerical simulation of saltwater intrusion mitigation by building a sluice in the North Branch of the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017

长江河口北支建闸对减轻盐水入侵的数值模拟

doi: 10.3969/j.issn.1000-5641.201941017
基金项目: 国家自然科学基金(41676083); 上海市科委重点项目(17DZ1201902); 上海教委高峰学科“岛屿大气与生态”
详细信息
    通讯作者:

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

  • 中图分类号: P751

Numerical simulation of saltwater intrusion mitigation by building a sluice in the North Branch of the Changjiang Estuary

  • 摘要: 枯季长江河口盐水入侵的最大特色是北支盐水倒灌, 它是南支东风西沙、太仓和陈行水库盐水的唯一来源, 也是青草沙水库盐水的主要来源. 考虑潮汐和气候态1月和2月的径流量与风况, 采用已严格验证过的长江河口盐水入侵三维数值模式, 模拟和分析北支上段建闸前后盐水入侵的变化. 模拟结果表明: 在北支上段建闸后, 整个南支全为淡水, 北支盐水倒灌南支的现象消失, 北支上段盐度明显下降; 在东风西沙、太仓和陈行水库取水口盐度接近0; 在青草沙水库取水口盐度大幅下降, 几乎所有时间盐度都低于0.45, 全为淡水. 数值试验中, 闸门的运行方式采用两种方案: 全天落潮流期间开闸、夜里涨潮流期间关闸、白天涨潮流期间开闸, 以及全天落潮流期间开闸、夜里和白天涨潮流期间关闸. 两者的试验结果中南支盐度变化几乎一致, 原因在于前者的运行方式已经使得北支上段盐水入侵大幅减弱, 出现盐度接近0.45的淡水区域; 即使白天涨潮流期间开闸, 其间进入南支的也是淡水, 并且增加了南支向海的总余流. 从数值模拟的结果和闸门运行的成本考虑, 推荐前者的北支建闸运行方案. 北支建闸极大地提高了上海东风西沙、陈行和青草沙水库取水时间, 同样极大地提高了江苏太仓水库的取水时间, 保障了两地的供水安全.
  • 图  1  长江口形势图

    注: 图中标注了水库和北支上段闸门位置

    Fig.  1  Map of the Changjiang Estuary

    图  2  模式计算区域及网格(a), 南北支分汊口(b)和南支下段局部放大网格(c)

    Fig.  2  Model domain and grids (a), enlarged views of the model grid at the bifurcation between the North Branch and South Branch (b), and in the lower reaches of the South Branch (c)

    图  3  大潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

    注: 图(a)和(b)为数值试验1, 图(c)和(d)为数值试验2, 图(e)和(f)为数值试验3; 红色等值线为盐度0.45(饮用水盐度标准)

    Fig.  3  Distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during spring tide

    图  4  南北支上段放大的大潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

    注: 图(a)和(b)为数值试验1, 图(c)和(d)为数值试验2, 图(e)和(f)为数值试验3

    Fig.  4  Enlarged views of the distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during spring tide at the bifurcation between the North Branch and South Branch

    图  5  小潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

    注: 图(a)和(b)为数值试验1, 图(c)和(d)为数值试验2, 图(e)和(f)为数值试验3

    Fig.  5  Distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during neap tide

    图  6  南北支上段放大的小潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

    注: 图(a)和(b)为数值试验1, 图(c)和(d)为数值试验2, 图(e)和(f)为数值试验3

    Fig.  6  Enlarged views of the dustribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during neap tide at the bifurcation between the North Branch and South Branch

    图  7  青草沙水库取水口水位随时间的变化(a); 东风西沙水库(b)、太仓水库(c)、陈行水库(d)和青草沙水库(e)表层盐度随时间的变化

    注: 黑线为数值试验1, 红线为数值试验2, 绿线为数值试验3, 红色虚线为盐度0.45(饮用水盐度标准)

    Fig.  7  Temporal variation of water elevation at the water intake of the Qingcaosha Reservoir (a); temporal variation of surface salinity at the water intake of the Dongfengxisha Reservoir (b). Taicang Reservoir (c), Chenhang Reservoir (d), and Qingcaosha Reservoir (e)

  • [1] 陆忠民, 关许为, 吴彩娥. 上海长江水源地规划与建设 [J]. 河海大学学报(自然科学版), 2010, 38(s2): 33-35.
    [2] 沈焕庭, 茅志昌, 朱建荣. 长江河口盐水入侵 [M]. 北京: 海洋出版社, 2003: 15-74.
    [3] 宋志尧, 茅丽华. 长江口盐水入侵研究 [J]. 水资源保护, 2002(3): 27-30. DOI:  10.3969/j.issn.1004-6933.2002.03.009.
    [4] WU H, ZHU J, CHEN B, et al. Quantitative relationship of runoff and tide to saltwater spilling over from the North Branch in the Changjiang Estuary: A numerical study [J]. Estuarine Coastal and Shelf Science, 2006, 69(1/2): 125-132.
    [5] 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.
    [6] 罗小峰, 陈志昌. 径流和潮汐对长江口盐水入侵影响数值模拟研究 [J]. 海岸工程, 2005, 24(3): 1-6. DOI:  10.3969/j.issn.1002-3682.2005.03.001.
    [7] 李路. 长江河口盐水入侵时空变化特征和机理 [D]. 上海: 华东师范大学, 2011.
    [8] LI L, ZHU J R, WU H. Impacts of wind stress on saltwater intrusion in the Yangtze Estuary [J]. Science China Earth Sciences, 2012, 55(7): 1178-1192. DOI:  10.1007/s11430-011-4311-1.
    [9] LI L, ZHU J, WU H, et al. Lateral saltwater intrusion in the North Channel of the Changjiang Estuary [J]. Estuaries and Coasts, 2014, 37(1): 36-55. DOI:  10.1007/s12237-013-9669-1.
    [10] ZHU J, DING P, ZHANG L, et al. Influence of the deep waterway project on the Changjiang Estuary [M]// WOLANSKI E. The environment in Asia Pacific harbours. Dordrecht, South Holland: Springer, 2006: 79-92.
    [11] 杨桂山, 朱季文. 全球海平面上升对长江口盐水入侵的影响研究 [J]. 中国科学, 1993(1): 69-76.
    [12] QIU C, ZHU J. 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.
    [13] 韩乃斌. 长江口南支河段氯度变化分析 [J]. 水利水运工程学报, 1983(1): 77-84.
    [14] 茅志昌, 沈焕庭, 姚运达. 长江口南支南岸水域盐水入侵来源分析 [J]. 海洋通报, 1993(3): 17-25.
    [15] 徐建益, 袁建忠. 长江口南支河段盐水入侵规律的研究 [J]. 水文, 1994(5): 1-6.
    [16] 茅志昌, 沈焕庭, 肖成猷. 长江口北支盐水倒灌南支对青草沙水源地的影响 [J]. 海洋与湖沼, 2001, 32(1): 58-66. DOI:  10.3321/j.issn:0029-814X.2001.01.010.
    [17] 顾玉亮, 吴守培, 乐勤. 北支盐水入侵对长江口水源地影响研究 [J]. 人民长江, 2003, 34(4): 1-3. DOI:  10.3969/j.issn.1001-4179.2003.04.001.
    [18] 王国峰, 乐勤. 长江口北支盐水入侵对陈行水库取水口的影响 [J]. 交通与港航, 2003, 17(4): 21-22. DOI:  10.3969/j.issn.1001-599X.2003.04.008.
    [19] 朱建荣, 吴辉, 顾玉亮. 长江河口北支倒灌盐通量数值分析 [J]. 海洋学研究, 2011(3): 1-7. DOI:  10.3969/j.issn.1001-909X.2011.03.002.
    [20] LYU H, ZHU J. 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.
    [21] 朱建荣, 鲍道阳. 近60年来长江河口河势变化及其对水动力和盐水入侵的影响I: 河势变化 [J]. 海洋学报, 2016, 38(12): 11-22.
    [22] 鲍道阳, 朱建荣. 近60年来长江河口河势变化及其对水动力和盐水入侵的影响II: 水动力 [J]. 海洋学报, 2017, 39(2): 1-15.
    [23] 鲍道阳, 朱建荣. 近60年来长江河口河势变化及其对水动力和盐水入侵的影响III: 盐水入侵 [J]. 海洋学报, 2017, 39(4): 1-14.
    [24] 吴辉, 朱建荣. 长江河口北支倒灌盐水输送机制分析 [J]. 海洋学报(中文版), 2007, 29(1): 17-25.
    [25] 朱建荣, 吴辉, 李路, 等. 极端干旱水文年(2006)中长江河口的盐水入侵 [J]. 华东师范大学学报(自然科学版), 2010(4): 1-6.
    [26] 王绍祥, 朱建荣. 不同潮型和风况下青草沙水库取水口盐水入侵来源 [J]. 华东师范大学学报(自然科学版), 2015(4): 65-76.
    [27] 朱建荣, 顾玉亮, 吴辉. 长江河口青草沙水库最长连续不宜取水天数 [J]. 海洋与湖沼, 2013, 44(5): 1138-1145.
    [28] 朱建荣, 吴辉. 长江河口东风西沙水库最长连续不宜取水天数数值模拟 [J]. 华东师范大学学报(自然科学版), 2013, 28(5): 1-8.
    [29] BLUMERG A F, MELLOR G L. A description of a three-dimensional coastal ocean circulation model [J]. Three-dimensional coastal ocean models, 1987(4): 1-16.
    [30] BLUMBERG A F. A primer for ECOM-si [R]. Technical report of HydroQual, 1994: 66.
    [31] CHEN C, ZHU J, ZHENG L, et al. A non-orthogonal primitive equation coastal ocean circulation model: application to Lake Superior [J]. Journal of Great Lakes Research, 2004, 30: 41-54. DOI:  10.1016/S0380-1330(4)70376-7.
    [32] WU H, ZHU J. Advection scheme with 3rd high-order spatial interpolation at the middle temporal level and its application to saltwater intrusion in the Changjiang Estuary [J]. Ocean Modelling, 2010, 33(1/2): 33-51.
    [33] 海洋图集编委会. 渤海、黄海、东海海洋图集(水文) [G]. 北京: 海洋出版社, 1992: 13-168.
  • [1] 张梦霞, 郑艳玲, 尹国宇, 董宏坡, 韩平, 高娟, 刘程, 常永凯, 刘敏, 侯立军.  纳米银对河口潮滩硝酸盐异化还原成铵过程的影响 . 华东师范大学学报(自然科学版), 2020, (3): 68-77. doi: 10.3969/j.issn.1000-5641.201941005
    [2] 张红丽, 尹国宇, 郑艳玲, 高娟, 高灯州, 常永凯, 刘程.  沉积物再悬浮对长江口潮滩上覆水体脱氮过程的影响 . 华东师范大学学报(自然科学版), 2020, (3): 78-87. doi: 10.3969/j.issn.1000-5641.201941007
    [3] 陆佳玉, 葛建忠, 丁平兴.  潮控型分汊河口分流过程探讨 . 华东师范大学学报(自然科学版), 2020, (3): 1-12. doi: 10.3969/j.issn.1000-5641.201941015
    [4] 鲁佩仪, 朱建荣, 钱伟伟, 袁琳.  崇明东滩鸟类栖息地优化工程区海堤水闸外侧引水渠冲淤数值模拟 . 华东师范大学学报(自然科学版), 2020, (3): 43-54. doi: 10.3969/j.issn.1000-5641.201941013
    [5] 唐川敏, 朱建荣.  长江河口水位上升对流场和盐水入侵的影响 . 华东师范大学学报(自然科学版), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
    [6] 苏爱平, 吕行行, 吴宇帆.  南水北调工程对长江口盐水入侵和淡水资源的影响 . 华东师范大学学报(自然科学版), 2020, (3): 32-42. doi: 10.3969/j.issn.1000-5641.201941031
  • 加载中
图(7)
计量
  • 文章访问数:  23
  • HTML全文浏览量:  41
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-14
  • 网络出版日期:  2020-05-29
  • 刊出日期:  2020-05-01

长江河口北支建闸对减轻盐水入侵的数值模拟

doi: 10.3969/j.issn.1000-5641.201941017
    基金项目:  国家自然科学基金(41676083); 上海市科委重点项目(17DZ1201902); 上海教委高峰学科“岛屿大气与生态”
    通讯作者: 朱建荣, 男, 教授, 研究方向为河口海洋学. E-mail: jrzhu@sklec.ecnu.edu.cn
  • 中图分类号: P751

摘要: 枯季长江河口盐水入侵的最大特色是北支盐水倒灌, 它是南支东风西沙、太仓和陈行水库盐水的唯一来源, 也是青草沙水库盐水的主要来源. 考虑潮汐和气候态1月和2月的径流量与风况, 采用已严格验证过的长江河口盐水入侵三维数值模式, 模拟和分析北支上段建闸前后盐水入侵的变化. 模拟结果表明: 在北支上段建闸后, 整个南支全为淡水, 北支盐水倒灌南支的现象消失, 北支上段盐度明显下降; 在东风西沙、太仓和陈行水库取水口盐度接近0; 在青草沙水库取水口盐度大幅下降, 几乎所有时间盐度都低于0.45, 全为淡水. 数值试验中, 闸门的运行方式采用两种方案: 全天落潮流期间开闸、夜里涨潮流期间关闸、白天涨潮流期间开闸, 以及全天落潮流期间开闸、夜里和白天涨潮流期间关闸. 两者的试验结果中南支盐度变化几乎一致, 原因在于前者的运行方式已经使得北支上段盐水入侵大幅减弱, 出现盐度接近0.45的淡水区域; 即使白天涨潮流期间开闸, 其间进入南支的也是淡水, 并且增加了南支向海的总余流. 从数值模拟的结果和闸门运行的成本考虑, 推荐前者的北支建闸运行方案. 北支建闸极大地提高了上海东风西沙、陈行和青草沙水库取水时间, 同样极大地提高了江苏太仓水库的取水时间, 保障了两地的供水安全.

English Abstract

朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行. 长江河口北支建闸对减轻盐水入侵的数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
引用本文: 朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行. 长江河口北支建闸对减轻盐水入侵的数值模拟[J]. 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
ZHU Jianrong, LU Peiyi, TANG Chuanmin, CHEN Qing, LÜ Hanghang. Numerical simulation of saltwater intrusion mitigation by building a sluice in the North Branch of the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
Citation: ZHU Jianrong, LU Peiyi, TANG Chuanmin, CHEN Qing, LÜ Hanghang. Numerical simulation of saltwater intrusion mitigation by building a sluice in the North Branch of the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
    • 长江河口为一特大型分汊河口, 径流量巨大, 淡水资源丰富. 河口旁边有我国经济发达地区上海和江苏, 人口密集, 随着社会经济的快速发展, 对用水数量和质量提出了更高的要求. 目前, 在长江口已建成了陈行水库、青草沙水库、太仓水库和东风西沙水库(见图1). 陈行水库位于南支南岸, 建成于1996年, 有效库容为8.60 × 106 m3, 日供水量2.06 × 106 m3, 上海西北地区约有200万人因此而受益[1]. 该水库目前供水占上海全部用水的20%. 青草沙水库于2010年建成, 有效库容为4.35 × 108 m3, 日供水量7.19 × 106 m3, 受益人口约有1 300万人. 青草沙水库面积大约为西湖面积的10倍, 是世界上最大的河口水库. 苏州太仓水库建成于2012年, 有效库容为15.0 × 106 m3, 日供水量6.0 × 105 m3. 东风西沙水库建成于2014年, 有效库容为8.9 × 106 m3, 日供水量2.1 × 105 m3, 崇明岛约有70万人受益. 目前, 上海大约80%的淡水取自长江口水源地. 在长江口建设水源地面临的最大问题是枯季盐水入侵. 盐水入侵季节, 若取水口盐度高于0.45 psu(饮用水盐度标准), 水库就不宜从长江取水.

      在河口, 河流淡水和海洋咸水交汇, 盐水入侵是河口的一个普遍现象. 在长江河口, 盐水入侵主要受径流量和潮汐的控制[1-7], 但也受风[8]、地形[9]、流域和河口工程[5, 10], 以及海平面上升[11-12]的影响. 发生在长江口盐水入侵最突出的现象是枯季大潮期间北支盐水倒灌进入南支[2, 4, 13-20]. 从20世纪50年代到2000年, 北支的自然演变和潮滩的人为圈围, 导致北支上段变为几乎垂直于南支, 而北支下段变为喇叭口形状[21]. 北支河势的演变阻碍了径流进入北支, 尤其在枯季, 导致北支潮差大于南支[22]. 北支巨大的潮差产生了显著的水平环流, 即在低径流量和大潮期间北支盐水净向陆流动, 进入南支后向海流动. 这个余输运就是形成北支盐水倒灌进入南支的动力成因[2, 4, 8, 23]. 由于落潮期间北支上段大片滩涂露出水面, 大潮涨潮期间从北支进入南支的盐水只有小部分退回到北支. 倒灌进入南支的盐水在径流作用下向下游输运, 在随后的中潮和小潮期间到达南支的中下段, 影响太仓、陈行和青草沙水库[24-26]. 东风西沙、太仓和陈行水库取水口的盐水全部来自北支盐水倒灌[2, 13-18]. 青草沙水库取水口的盐度则来自北支倒灌和北港向陆的盐水入侵[2, 12, 16-18, 26].

      水库的持续供水时间取决于有效库容和日供水量, 而有效库容的设计取决于最长连续不宜取水时间. 以1978—1979年特枯水文年为计算年(径流量保证率为97%), 青草沙水库和东风西沙水库的最长不宜取水时间分别为68 d和26 d[27-28]. 在水库的运行中, 较难预测径流量未来变化, 也就难于预测盐水入侵未来变化. 目前, 长江河口盐水入侵预报时效在10 d之内, 预报结果可靠, 但对长时间尺度预报来说较难正确预报径流量和风况, 盐水入侵预报可靠性就较差. 这就对水库的取水和淡水资源的安全构成挑战, 需要提出减缓盐水入侵、保障水库取水安全的对策.

      图  1  长江口形势图

      Figure 1.  Map of the Changjiang Estuary

      本文对北支建闸方案开展数值模拟, 对比分析建闸前后盐水入侵变化和南支水源地取水口盐度变化, 为长江口淡水资源的充分利用和保障供水安全提供科技指导.

    • 本文数值模式采用基于POM模式[29]发展起来的ECOM-si[30], 之后从正交曲线坐标发展到非正交曲线坐标[31]. Wu和Zhu研发了一个物质平流格式的HSIMT-TVD[32], 这个格式以三阶精度有效消除数值频散和大幅降低数值耗散. 模式已成功地应用于长江河口水动力和盐水入侵的研究中[5, 8, 20].

      模式范围包括整个长江河口、杭州湾和邻近海区, 纬向从117.5°E~125°E, 径向从27.5°N~33.7°N(见图2a). 总的网格数在水平方向上加密到337 × 225, 在垂向方向上10个均匀的σ层, 以较好地分辨河口地形. 最小的网格距在南北支分汊口约为100 m, 以更好地模拟北支盐水倒灌(见图2b); 分辨率在南支下段约为200 m(见图2c). 模式网格的正交性和平滑性良好, 拟合了分汊口河道岸线和深水航道工程的形状. 在口门的分辨率为300~500 m, 随后逐渐减低到外海开边界处的10 km. 地形资料采用2017年整个长江口实测资料. 模式积分时间步长设置为40 s.

      外海开边界条件来自NaoTide数据集(http://www.miz.nao.ac.jp/), 由16个分潮驱动: M2, S2, N2, K2, K1, O1, P1, Q1, MU2, NU2, T2, L2, 2N2, J1, M1和OO1. 大通实测月均径流量作为模式河流边界条件. 1月和2月为全年径流量最低期间, 也是盐水入侵最为严重月份. 1950年以来大通1月和2月平均径流量分别为11 500 m3/s和11 800 m3/s. 风数据采用NECP资料, 分辨率0.5° × 0.5°, 10年半月平均. 长江口外盐度初始条件来自黄海东海海洋图集(水文)[33], 口门内来自多次观测资料. 模式从1月1日开始计算, 至2月28日结束计算, 从2月下旬模拟结果中给出大潮和小潮期间涨憩和落憩时刻的盐度分布. 涨憩和落憩参考点为北支建闸处.

      图  2  模式计算区域及网格(a), 南北支分汊口(b)和南支下段局部放大网格(c)

      Figure 2.  Model domain and grids (a), enlarged views of the model grid at the bifurcation between the North Branch and South Branch (b), and in the lower reaches of the South Branch (c)

      上述数值模式已经在长江河口验证了许多次, 结果表明, 此模式能成功模拟长江河口水动力过程和盐水入侵[5, 8, 20]. 因论文篇幅限制, 模式验证过程的详细描述可参见上述文献.

      本文设置3个数值试验. 数值试验1: 不考虑北支建闸. 数值试验2: 考虑北支建闸(闸门位置见图1). 长江河口为典型的半日潮, 一天有两次涨潮和两次落潮, 北支盐水入侵是由涨潮期间带来的, 考虑夜里涨潮流期间关闸, 白天涨潮流期间开闸, 夜里和白天落潮流期间开闸. 数值试验3: 考虑北支建闸, 与数值试验2唯一不同的是夜里和白天涨潮流期间均关闸. 本文数值模式为原代码程序, 能方便地设置闸门开启和关闭的条件. 与一些商业软件无法修改源代码程序相比, 具有明显的优点.

    • 由大潮期间涨憩和落憩时刻垂向平均盐度分布可见(见图3), 在潮汐和枯季气候态径流量和风况作用下, 在北支上段建闸前(数值试验1)北支被高盐水所占据, 大潮期间南支上段出现高盐水团, 这是由北支盐水倒灌产生的, 南支下段出现盐度低于0.45的淡水区域. 小潮期间, 已经倒灌的盐水团和下游的淡水区域在径流作用下向下游移动. 在北港、北槽和南槽拦门沙区域, 出现强烈的盐度锋面. 在这3个分汊河道, 盐水入侵南槽最强、北槽次之, 北港最弱. 北港北侧出现明显的向东北输运的冲淡水, 这是由潮泵输运作用造成的[9]. 在数值试验2中, 尽管在白天涨潮流期间开闸, 整个南支全为淡水, 北支盐水倒灌南支的现象消失. 北支上段盐度明显下降, 小潮落憩时刻长兴岛北侧北港全为淡水. 在数值试验3中, 白天和夜里涨潮流期间均关闸, 北支盐水倒灌现象彻底消失, 4个水库的取水口均为淡水, 北支盐度下降更明显. 为更清楚地显示北支上段建闸前后盐度分布和变化情况, 对北支和南支上段作放大展示(见图4). 建闸前北支上段为高盐水占据. 在数值试验2中, 涨潮流期间白天开闸、夜里关闸, 南支上段没有盐度高于0.45的区域, 北支上段盐度下降, 落憩时刻更显著. 在数值试验3, 涨潮流期间白天和夜里均关闸, 北支上段涨憩和落憩时刻均为淡水, 落憩时刻淡水区域更大.

      图  3  大潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

      Figure 3.  Distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during spring tide

      在小潮期间(见图5), 在北支建闸前, 南支上段和南支南侧为盐度小于0.45的淡水区域, 大潮期间北支倒灌进入南支的盐水在径流作用下向下游移动, 影响陈行和青草沙水库. 在北港、北槽和南槽拦门沙区域, 等盐度线向海移动, 相比于大潮期间盐水入侵减弱. 北港北侧向东北输运的冲淡水消失, 原因在于小潮期间潮泵运输作用减弱, 北风作用相对占优势[9]. 在数值试验2中, 南支全域均为淡水, 北支上段盐水入侵减弱. 在数值试验3中, 与数值试验2比较, 南支盐度变化不明显, 但北支上段和中段盐度下降更为显著. 从北支和南支上段局部放大盐度分布来看(见图6), 上述盐度变化更为清晰. 建闸前北支上段被高盐水占据; 在数值试验2中, 北支上段为淡水; 在数值试验3中, 淡水区域已移至北支中段. 由此可见, 建闸后, 不仅完全消除了南支来自北支的盐水入侵, 在北支上段和中段也出现了持续的淡水.

      图  4  南北支上段放大的大潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

      Figure 4.  Enlarged views of the distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during spring tide at the bifurcation between the North Branch and South Branch

      从青草沙水库取水口水位和4个水源地取水口盐度随时间变化可见(见图7), 长江口潮汐为半日潮, 具有日不等现象, 半月的大小潮变化明显. 在东风西沙水库取水口, 建闸前最高盐度达到1.5, 在一个大小潮周期内约三分之一的时段盐度低于0.45. 建闸后, 在数值试验2中和数值试验3中盐度均接近0. 在太仓水库取水口, 建闸前大约有一半时间盐度小于0.45, 建闸后盐度接近0. 在陈行水库取水口, 情况与太仓水库一致. 在青草沙水库取水口, 建闸前最高盐度约为1.2, 在一个大小潮周期内大约五分之二时段能取到淡水, 建闸后几乎全部时间均能取到淡水. 值得注意的是, 数值试验3的盐度反而比数值试验2的盐度高, 原因在于数值试验3白天和夜里涨潮流期间均关闭水闸, 彻底消除了北支水体倒灌进入南支, 使得南支向海的总余流减小, 北港海外向陆的盐水入侵增强.

      图  5  小潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

      Figure 5.  Distribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during neap tide

      从上述模拟的盐度平面分布和水库取水口盐度变化过程可以看到, 数值试验2和数值试验3的结果几乎一致. 原因在于在数值试验2中, 尽管白天涨潮流期间开闸, 但因夜里涨潮流期间关闸、白天和夜里落潮流期间开闸, 使得北支上段盐水入侵大幅减弱, 出现盐度接近0.45的淡水. 这样, 即使白天涨潮流期间开闸, 其间进入南支的也是淡水, 与来自徐六泾上游的河水无异. 从数值模拟的结果和闸门运行的成本考虑, 本文推荐数值试验2中的北支建闸运行方案.

      图  6  南北支上段放大的小潮涨憩(左侧)和落憩(右侧)时刻垂向平均盐度分布

      Figure 6.  Enlarged views of the dustribution of vertically averaged salinity at the flood slack (left panel) and ebb slack (right panel) during neap tide at the bifurcation between the North Branch and South Branch

    • 长江河口盐水入侵的最大特色是北支盐水倒灌, 它是南支东风西沙、太仓和陈行水库盐水入侵的唯一来源, 也是青草沙水库盐水入侵的主要来源. 减轻北支盐水倒灌是延长水源地取水口取水时间和保障供水安全的有效方法. 考虑潮汐和气候态1月和2月的径流量和风况, 本文采用已严格验证的长江河口盐水入侵三维数值模式, 模拟和对比北支上段建闸前后盐水入侵的变化. 建闸方案分为2个, 全天落潮流期间打开闸门, 一个是涨潮流期间夜里关闸、白天开闸(数值试验2), 另一个是涨潮流期间夜里和白天均关闸(数值试验3). 模拟结果再现了建闸前(数值试验1)长江河口北支被高盐水占据、北支盐水倒灌、拦门沙盐度锋面和空间差异等特征. 在数值试验2中, 整个南支全为淡水, 北支盐水倒灌南支的现象消失, 北支上段盐度明显下降. 在东风西沙、太仓和陈行水库取水口盐度接近0. 在青草沙水库取水口, 盐度大幅下降, 几乎所有时间盐度都低于0.45. 在数值试验3中, 北支盐水倒灌现象彻底消失, 北支盐度下降更明显, 北支上段和中段为淡水区域. 小潮期间北支上段和中段淡水范围比大潮期间大. 建闸后, 不仅完全消除了南支来自北支的盐水入侵, 在北支上段和中段也出现了持续的淡水. 在东风西沙、太仓和陈行水库取水口, 盐度接近0, 与数值试验2结果一致. 在青草沙水库取水口, 几乎所有时间都能取淡水, 但盐度反而比数值试验2的盐度高. 原因在于数值试验3的水闸运行方式彻底消除了北支水体倒灌进入南支, 使得南支向海的总余流减小, 北港海外向陆的盐水入侵增强.

      图  7  青草沙水库取水口水位随时间的变化(a); 东风西沙水库(b)、太仓水库(c)、陈行水库(d)和青草沙水库(e)表层盐度随时间的变化

      Figure 7.  Temporal variation of water elevation at the water intake of the Qingcaosha Reservoir (a); temporal variation of surface salinity at the water intake of the Dongfengxisha Reservoir (b). Taicang Reservoir (c), Chenhang Reservoir (d), and Qingcaosha Reservoir (e)

      数值试验2和数值试验3的结果几乎一致. 原因在于数值试验2的运行方式已经使得北支上段盐水入侵大幅减弱, 出现盐度接近0.45的淡水. 即使白天涨潮流期间开闸, 其间进入南支的也是淡水. 从数值模拟的结果和闸门运行的成本考虑, 本文推荐数值试验2中的北支建闸运行方案.

      北支建闸后可极大地提高上海的陈行水库、青草沙水库和东风西沙水库取水时间, 保障供水安全, 同样可极大地提高江苏太仓水库的取水时间. 北支建闸工程对沪苏两地均是有益的.

参考文献 (33)

目录

    /

    返回文章
    返回