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

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

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

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

俄罗斯《文摘杂志》收录

留言板

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

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

长江河口水位上升对流场和盐水入侵的影响

唐川敏 朱建荣

唐川敏, 朱建荣. 长江河口水位上升对流场和盐水入侵的影响[J]. 华东师范大学学报(自然科学版), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
引用本文: 唐川敏, 朱建荣. 长江河口水位上升对流场和盐水入侵的影响[J]. 华东师范大学学报(自然科学版), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
TANG Chuanmin, ZHU Jianrong. Influence of water level rise on currents and saltwater intrusion in the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
Citation: TANG Chuanmin, ZHU Jianrong. Influence of water level rise on currents and saltwater intrusion in the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001

长江河口水位上升对流场和盐水入侵的影响

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

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

  • 中图分类号: P731.2

Influence of water level rise on currents and saltwater intrusion in the Changjiang Estuary

  • 摘要: 基于长江河口水动力和盐水入侵三维数值模式ECOM-si, 通过数值模拟分析长江河口水位上升对径流、潮流和风生流的影响, 以及在多种动力因子综合作用下对流场和盐水入侵的影响. 数值实验结果表明, 在长江河口水位上升30 cm 的情况下, 各河道横截面面积增大, 向海的径流流速减小; 潮流随水深增深略微增大; 枯季北风作用产生的北港进、南港出的水平风生环流加强, 在北支向陆的风生流有所加强. 水位上升后, 北支盐水入侵增强; 南支中段盐度变化不明显; 北港、北槽盐水入侵随水位增加变化最为显著, 小潮期间盐度增大值大于1, 大潮期间增幅有所降低, 北港北汊受水深增加盐水入侵变化最为强烈; 南槽口门处滩地由于水位增加, 非线性效应减弱, 盐度不同程度的降低. 水位上升后南支水源地三个水库取水口盐度均有所上升, 减少了可取水时间, 不利于供水安全.
  • 图  1  长江河口形势图

    注: 黑点为水库取水口位置, 三角形为模式输出点

    Fig.  1  Map of the Changjiang Estuary

    图  2  模型计算区域和网格(a), 及放大的南支分汊口区域网格(b)和南北槽分汊口区域网格(c)

    Fig.  2  Model domain and grid (a), enlarged view of grid near the bifurcation of the South and North Branch (b) and near the bifurcation of the South and North passage (c)

    图  3  只有径流作用情况下水位上升前垂向平均流场分布(a)和水位上升后垂向平均流场变化(b)

    Fig.  3  Distribution of vertically averaged currents before the water level rise (a) and the difference in vertically averaged currents between the results before and after the water level rise (b) under only the force of river discharge

    图  4  只有潮汐作用情况下在模式输出点S1(上)、S2(中) 和S3(下)流速和流向随时间变化

    注: 红线为水位上升前, 黑线为水位上升后

    Fig.  4  Temporal variations in current speed and direction at site S1 (upper panel), S2 (middle panel), and S3 (lower panel)

    图  5  只有风作用的情况下, 水位上升前流场垂向平均分布图(a)、水位上升后流场变化(b)

    Fig.  5  Distribution of vertically averaged current before the water level rise (a) and the difference in vertically averaged current between the results before and after the water level rise (b) under only the force of wind

    图  6  水位上升后大潮期间涨潮垂向平均分布变化图(a)、落潮垂向平均分布变化图(b)

    Fig.  6  Variation in the vertically averaged current after the elevation rise in flood (a) and ebb (b) during spring tides

    图  7  水位上升前大潮期间盐度垂向平均分布图(a)及水位上升后盐度变化(b)

    Fig.  7  Distribution of the vertically averaged salinity before the water level rise (a) and the difference in vertically averaged salinity between the results before and after the water level rise (b) during spring tide

    图  8  水位上升后小潮期间涨潮垂向平均分布变化图(a)、落潮垂向平均分布变化图(b)

    Fig.  8  Variation in the vertically averaged current after the elevation rise in flood (a) and ebb (b) during neap tides

    图  9  水位上升前小潮期间盐度垂向平均分布图(a)及水位上升后盐度变化(b)

    Fig.  9  Distribution of the vertically averaged salinity before the water level rise (a) and the difference in vertically averaged salinity between the results before and after the water level rise (b) during neap tide

    图  10  东风西沙水库、陈行水库、青草沙水库取水口表层(左)和底层(右)盐度随时间的变化

    注: 红色实线为水位上升前盐度, 黑色实线为水位上升后盐度; 灰色虚线为0.45饮用水标准

    Fig.  10  Variation in salinity on the surface (right) and bottom (left) at the water intake of Dongfengxisha Reservoir, Chenhang Reservoir, and Qingcaosha Reservoir

  • [1] QIU C, ZHU J R, GU Y L. Impact of seasonal tide variation on saltwater on saltwater intrusion in the ChangJiang River Estuary [J]. Chinese Journal of Oceanology and Limnology, 2012, 30(2): 342-351. DOI:  10.1007/s00343-012-1115-x.
    [2] 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.
    [3] LÜ 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.
    [4] 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.
    [5] LI L, ZHU J R, WU H, et al. Lateral saltwater intrusion in the North Channel of the Changjiang Estuary [J]. Estuaries & Coasts, 2014, 37(1): 36-55.
    [6] ZHU J R, WU H, LI L. Hydrodynamics of the Changjiang Estuary and Adjacent Seas [M]//Ecological Continuum from the Changjiang (Yangtze River) Watersheds to the East China Sea Continental Margin. [S.l.]: Springer International Publishing, 2015: 19-45.
    [7] 杨桂山, 朱季文. 全球海平面上升对长江口盐水入侵的影响研究 [J]. 中国科学(B辑), 1993, 23(1): 69-76.
    [8] HANSEN N D, KOTERA A, SAKAMOTO T, et al. Sensitivity of salinity intrusion to sea level rise and river flow change in Vietnamese Mekong Delta impacts on availability of irrigation water for rice cropping [J]. Journal of Agricultural Meteorology, 2008, 64(3): 167-176. DOI:  10.2480/agrmet.64.3.4.
    [9] MACCREADY P, GEYER W R. Advances in estuarine physics [J]. Annual Review of Marine Science, 2010(2): 35-58.
    [10] BLUMBERG A F, MELLOR G L. A description of a three-dimensional coastal ocean circulation model [M]//Three-Dimensional Coastal Ocean Models. Washington, DC: American Geophysical Union (AGU), 2013: 1-16.
    [11] 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.
    [12] 朱建荣, 吴辉, 顾玉亮. 长江河口北支倒灌盐通量数值分析 [J]. 海洋学研究, 2011, 29(3): 1-7. DOI:  10.3969/j.issn.1001-909X.2011.03.002.
    [13] 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.
    [14] 李永平, 秦曾灏, 段义宏. 上海地区海平面上升趋势的预测和研究 [J]. 地理学报, 1998, 53(5): 393-403. DOI:  10.3321/j.issn:0375-5444.1998.05.002.
    [15] WU H, ZHU J R. 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: 33-51. DOI:  10.1016/j.ocemod.2009.12.001.
    [16] 李路. 长江河口盐水入侵时空变化特征和机理 [D]. 上海: 华东师范大学, 2011.
  • [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, (2): 55-63. doi: 10.3969/j.issn.1000-5641.201922010
    [6] 苏爱平, 吕行行, 吴宇帆.  南水北调工程对长江口盐水入侵和淡水资源的影响 . 华东师范大学学报(自然科学版), 2020, (3): 32-42. doi: 10.3969/j.issn.1000-5641.201941031
    [7] 朱建荣, 鲁佩仪, 唐川敏, 陈晴, 吕行行.  长江河口北支建闸对减轻盐水入侵的数值模拟 . 华东师范大学学报(自然科学版), 2020, (3): 13-22. doi: 10.3969/j.issn.1000-5641.201941017
  • 加载中
图(10)
计量
  • 文章访问数:  76
  • HTML全文浏览量:  58
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-08
  • 网络出版日期:  2020-05-29
  • 刊出日期:  2020-05-01

长江河口水位上升对流场和盐水入侵的影响

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

摘要: 基于长江河口水动力和盐水入侵三维数值模式ECOM-si, 通过数值模拟分析长江河口水位上升对径流、潮流和风生流的影响, 以及在多种动力因子综合作用下对流场和盐水入侵的影响. 数值实验结果表明, 在长江河口水位上升30 cm 的情况下, 各河道横截面面积增大, 向海的径流流速减小; 潮流随水深增深略微增大; 枯季北风作用产生的北港进、南港出的水平风生环流加强, 在北支向陆的风生流有所加强. 水位上升后, 北支盐水入侵增强; 南支中段盐度变化不明显; 北港、北槽盐水入侵随水位增加变化最为显著, 小潮期间盐度增大值大于1, 大潮期间增幅有所降低, 北港北汊受水深增加盐水入侵变化最为强烈; 南槽口门处滩地由于水位增加, 非线性效应减弱, 盐度不同程度的降低. 水位上升后南支水源地三个水库取水口盐度均有所上升, 减少了可取水时间, 不利于供水安全.

English Abstract

唐川敏, 朱建荣. 长江河口水位上升对流场和盐水入侵的影响[J]. 华东师范大学学报(自然科学版), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
引用本文: 唐川敏, 朱建荣. 长江河口水位上升对流场和盐水入侵的影响[J]. 华东师范大学学报(自然科学版), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
TANG Chuanmin, ZHU Jianrong. Influence of water level rise on currents and saltwater intrusion in the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
Citation: TANG Chuanmin, ZHU Jianrong. Influence of water level rise on currents and saltwater intrusion in the Changjiang Estuary[J]. Journal of East China Normal University (Natural Sciences), 2020, (3): 23-31. doi: 10.3969/j.issn.1000-5641.201941001
    • 长江河口地区是中国人口密集度最高、经济最发达的地区之一. 长江是我国第一大河, 全长约6 300 km. 长江河口是中国最大的河口, 地形复杂(见图1). 河口是上游下泄淡水与外海高盐水交汇的区域, 深受盐水入侵的影响. 盐水入侵与自然因子、生态环境、淡水资源分布和城市供水等密切相关. 研究长江河口水位上升对长江河口水动力过程和盐度入侵具有重要的理论意义和应用价值.

      图  1  长江河口形势图

      Figure 1.  Map of the Changjiang Estuary

      潮汐、径流量和风应力是影响长江河口盐水入侵的主要因素. 周期性潮流变化和涨落潮造成的水位变化在口内盐度分布中扮演着重要角色, 在北支盐水倒灌现象中尤为显著[1]. 长江径流量具有显著的季节性变化, 在低径流的情况下长江口盐水入侵加剧, 所有汊道的等盐度线向陆收缩, 对长江口水源地取水安全造成不利影响[2-3]. 长江口外风场的季节性变化显著, 冬季盛行偏北风, 夏季盛行偏南风. 冬季北风的大小和持续时间会显著影响长江河口盐水入侵[4-5].

      近50年来, 我国海平面呈波动上升的趋势, 根据国家海洋局2013年以来发布的《中国海平面公报》, 1980—2012年长江河口沿海相对海平面上升幅度超过150 mm, 即平均每年海平面上升4.5 mm左右. 径流量变化会造成长江河口水位的变化, 洪季高径流量期间水位比枯季低径流量期间水位高约50 cm[6]. 冬季强偏北风产生沿岸增水波, 导致长江口水位上升. 长江河口水位上升会导致咸淡水交换加强[7], 盐度梯度和向陆斜压压强梯度力增大, 使得外海高盐水向口内运动, 盐水入侵增强[8-9]. 长期的水位上升将加深河口地区的水深、淹没潮滩、影响河口水动力过程、加剧外海盐水入侵, 对长江河口淡水资源有重要影响. 研究长江河口水位上升对流场和盐水入侵的影响可以为长江口水源地的建设和水库取水调度提供理论依据, 保障供水安全.

      本文应用长江河口三维水动力数值模式, 从多个动力因子着手, 逐一对比分析河口区域水位上升对流场和盐水入侵程度的影响.

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

      本文研究长江河口水位上升对流场和盐水入侵的影响, 模式计算范围包括整个长江河口、杭州湾和邻近海域(见图2a). 河流上游边界设置在大通, 动量方程的边界条件采用大通水文站实测径流量. 模式采用非正交坐标曲线网格, 较好地拟合长江河口岸线, 并对局部地区如南北支分汊口、深水航道等进行了加密处理 (见图2b2c). 外海开边界由潮位和余水位驱动, 潮位考虑由16个主要分潮M2、S2、N2、K2、K1、O1、P1、Q1、MU2、NU2、T2、L2、2N2、J1、M1及OO1的调和常数合成得到, 余水位由黄海东海大区域模式计算结果提供. 初始温盐场在长江口外通过《渤海黄海东海海洋图集(水文)》数字化数据插值得到, 口内由多年枯季实测资料插值获得.

      图  2  模型计算区域和网格(a), 及放大的南支分汊口区域网格(b)和南北槽分汊口区域网格(c)

      Figure 2.  Model domain and grid (a), enlarged view of grid near the bifurcation of the South and North Branch (b) and near the bifurcation of the South and North passage (c)

    • 长江河口盐水入侵主要受制于径流量、潮汐和风况. 长江径流量具有明显的季节性变化, 5—10月为长江洪季, 11月到次年4月为枯季. 一般一年中1、2月份径流量最小, 自1950年以来大通水文站1月和2月的平均值分别为11 200 m3/s和11 800 m3/s. 枯季代表性月份为1月和2月, 这期间长江口盐水入侵最为严重[12]. 为了数值实验更具有普遍性, 选取1、2月份大通水文站月平均径流量11 550 m3/s, 即1—2月长江径流量气候态平均值作为上游开边界条件. 同样, 长江河口风场也取1—2月气候态平均值, 为北风5 m/s. 数值模拟时段为1月1日至2月28日.

      对于长江河口水位上升值, 主要考虑海平面变化. 根据李永平等[14]的研究成果及国家海洋局发布的中国海平面上升公告等, 得出至2050年长江口及邻近海域相对水平面上升值约为30 cm. 本文综合上述因素设置余水位(水位上升)为30 cm, 即通过模式全域增加水深30 cm来体现.

      为分析长江河口水位上升对流场和盐水入侵的影响, 设置4组数值实验. 先仅考虑径流、潮汐、风单个动力因子, 分别讨论水位上升30 cm前后对流场的影响. 最后综合考虑多个动力因子模拟和分析水位上升30 cm前后对流场、盐水入侵和水库取水口盐度的影响.

    • 在径流量为11 550 m3/s和水位未上升前, 长江河口径流流速最大可达18 cm/s, 整体呈现出由陆向海河势转宽流速逐渐减小的变化趋势(见图3a). 极少长江水进入北支, 径流流速约为2 cm/s. 南支径流受科氏力作用右偏, 沿南岸向海, 至青草沙分为南、北港, 北港为径流的主要入海口. 北港径流受口门浅滩影响, 一部分沿北港北汊流出, 流速快速减小; 另一部分继续向东并向右偏转. 南港径流被九段沙分流后, 沿深水航道向海保持约10 cm/s的速度, 而南槽部分由于滩地众多, 流速不断减小.

      图  3  只有径流作用情况下水位上升前垂向平均流场分布(a)和水位上升后垂向平均流场变化(b)

      Figure 3.  Distribution of vertically averaged currents before the water level rise (a) and the difference in vertically averaged currents between the results before and after the water level rise (b) under only the force of river discharge

      长江河口水位上升30 cm后, 因为水位上升、河道横断面增大, 在径流量不变的情况下, 径流流速减小, 变化量小于0.2 cm/s, 且流速大的区域变化更显著(见图3b). 而北支则因为过水断面的增大, 流入的径流量增加, 导致流速增大, 最大可增加1.5 cm/s左右.

    • 潮流受水位上升影响变化微小, 平面图难以区别, 故本文选取模式结果3个输出点S1、S2、S3(位置见图1), 对潮流流速和流向随时间变化过程作比较分析. S1位于北支上段, 大、小潮流速差异大, 水位上升前, 大潮期间流速最大约为1.3 m/s, 小潮期间最大仅为0.5 m/s左右(见图4). 水位上升后, 流速有较明显的增加, 且大潮期间增幅大于小潮期间, 最大增幅约为0.2 m/s. 点S2在青草沙水库取水口附近, S3位于北港拦门沙, 最大流速均在2 m/s左右, 水位上升后流速有微小增加.

      图  4  只有潮汐作用情况下在模式输出点S1(上)、S2(中) 和S3(下)流速和流向随时间变化

      Figure 4.  Temporal variations in current speed and direction at site S1 (upper panel), S2 (middle panel), and S3 (lower panel)

    • 在枯季平均北风作用下, 苏北沿岸和长江口外南向风生流流速较大, 北支口门附近受浅滩影响, 向陆流减弱, 少量水体从北支向口内输移, 最终从南北支分汊口进入南支(见图5a); 另一部分进入北港, 北港北汊流速最大约0.25 m/s, 水体上溯至青草沙上端, 向南转向进入南港, 在南北港间形成风生水平环流. 南港至深水航道段流速维持在0.1 m/s, 向海流出, 汇入南下的沿岸流中, 而南槽流速在九段沙附近迅速减小.

      图  5  只有风作用的情况下, 水位上升前流场垂向平均分布图(a)、水位上升后流场变化(b)

      Figure 5.  Distribution of vertically averaged current before the water level rise (a) and the difference in vertically averaged current between the results before and after the water level rise (b) under only the force of wind

      水位上升30 cm后, 北风作用导致北支倒灌加强, 从北支进入南支的水体增多(见图5b). 南北港水平环流增强, 北港口门附近流速增大0.17 m/s, 南支下段流速增大约0.2 m/s. 进而南港下泄的流速增大, 最大增值为0.2 m/s左右. 深水航道由于水深较大, 河道阻碍较少, 为南港区域主要下泄通道, 向海流速显著增强, 约0.15 m/s; 南槽向海流减弱, 水位上升后, 沿岸流增强, 在科氏力作用下右偏进入南槽, 流速变化值可达0.25 m/s.

    • 选取大小潮周期内各3个完整潮周期的数据进行统计分析, 涨、落潮分别取6个涨潮、6个落潮的平均值. 本文取青草沙水库取水口作为流速涨、落急和涨、落憩的参考点.

      在多个动力因子综合作用下, 大潮期间, 水位上升之后涨潮流、落潮流均增强, 涨潮变化程度大于落潮(见图6). 北支涨潮流整体增加约0.02 m/s, 部分水体进入南支. 南支中段以下流速增幅明显, 越靠近口门增幅越大, 口门附近流速增大0.1 m/s左右, 其中深水航道涨潮流变化最大, 约为0.16 m/s. 落潮流时, 北支落潮流由陆向海表现为先增大后减小再增大, 北支中上段流速增大约0.05 m/s. 南支落潮流变化趋势与北支相似, 北港在北汊附近流速增幅较大, 为0.07 m/s左右, 深水航道北岸外侧落潮流流速变化最大, 约0.1 m/s. 深水航道流速均增大约0.02 m/s, 而南槽流速变化最大值出现在南槽上段, 约为0.09 m/s.

      图  6  水位上升后大潮期间涨潮垂向平均分布变化图(a)、落潮垂向平均分布变化图(b)

      Figure 6.  Variation in the vertically averaged current after the elevation rise in flood (a) and ebb (b) during spring tides

      大潮期间, 整个北支都被高盐水占据(见图7a), 北支上端盐度变化剧烈, 分汊口处盐度锐减, 南支上段盐度值大于0.45, 再现北支盐水倒灌现象. 北港是径流的主要出海口, 出口门的冲淡水向北延伸至北支口门附近. 大潮期间南港盐水入侵略强于北港, 但从九段沙开始下端盐度有所降低. 水位上升后, 盐水入侵增强主要体现在北支和南、北港(见图7b). 北支盐度增加约0.2, 北支上段增加0.4左右. 北支倒灌增强, 南支上段盐度略微增加. 南、北港上段盐度基本不变, 中下段外海高盐水入侵明显增强, 北港盐度增大值及增大范围均大于南港. 北港中下段为盐度锋面, 盐度梯度较大, 尤其是拦门沙区域盐度值增加约0.5. 南、北槽整体盐度增加, 增大约0.3. 而横沙东滩后端和崇明东滩随水位上升盐度明显降低, 崇明东滩盐度降低约1.5, 横沙东滩降低约0.4. 吴辉[15]指出由于潮滩水深较浅, 在涨落潮期间会形成非线性效应, 横沙东滩后端和崇明东滩区域形成局地环流, 促进高盐水汇集, 大潮期间由于潮汐强度较大, 效应显著. 水位上升后, 滩地附近水深加深, 相应的非线性效应减弱. 这就解释了本文模拟的水位上升后潮滩盐度下降的现象.

      图  7  水位上升前大潮期间盐度垂向平均分布图(a)及水位上升后盐度变化(b)

      Figure 7.  Distribution of the vertically averaged salinity before the water level rise (a) and the difference in vertically averaged salinity between the results before and after the water level rise (b) during spring tide

      长江河口水位上升后小潮期间涨、落潮流平面分布(见图8)与大潮期间相似, 但流速变化更小. 涨潮流依旧在深水航道区域变化最大, 深水航道上段流速增大约0.04 m/s, 下段最大值可达0.2 m/s左右. 北港北汊流速增大值最大约为0.05 m/s, 而拦门沙外涨潮流略微减小. 南槽九段沙以外流速大面积增大, 靠近深水航道南岸流速增大0.1 m/s左右. 落潮时, 南、北支上段落潮流显著增强. 北支增大0.02 m/s, 而南支上端增大约0.05 m/s. 选取青草沙水库取水口为参考点时口外部分区域已经开始涨潮. 横沙东滩以外靠近深水航道北岸外侧落潮流变化最大, 落潮流流速增大0.3 m/s左右.

      图  8  水位上升后小潮期间涨潮垂向平均分布变化图(a)、落潮垂向平均分布变化图(b)

      Figure 8.  Variation in the vertically averaged current after the elevation rise in flood (a) and ebb (b) during neap tides

      小潮期间, 北支盐度值依然增加且增加量比大潮期间略有上升, 北支上段盐度增加值均可达1以上(见图9). 北港在小潮期间受水位上升影响更为剧烈, 北港中下段连接北支口门, 整个区域盐度增加, 北港北汊附近盐度变化最大, 约1.5左右. 李路[16]对北港北汊受盐水入侵的影响进行了数值模拟, 认为该汊道加深后盐水入侵加强. 本文模拟结果与这一结论一致. 河口水位上升后, 北港北汊水深增加, 从而增加沿该汊道进入北港的水、盐通量. 南港及南、北槽上段盐度增大0.3左右. 深水航道下段盐度值增加0.1左右; 南槽九段沙下段及横沙东滩往外盐度值降低, 最小区域减小1左右. 原因同样是浅滩水位上升导致水深增加潮汐非线性效应减弱产生的.

      图  9  水位上升前小潮期间盐度垂向平均分布图(a)及水位上升后盐度变化(b)

      Figure 9.  Distribution of the vertically averaged salinity before the water level rise (a) and the difference in vertically averaged salinity between the results before and after the water level rise (b) during neap tide

      东风西沙水库、陈行水库和青草沙水库3个水库取水口分布在南支不同区域, 其表、底层盐度随时间的变化见图10, 水位上升后盐度值均有增加. 可以看出水位上升前后盐度值及变化趋势基本相同, 底层变化明显比表层剧烈. 东风西沙水库和陈行水库位于南支上段, 盐度只受北支倒灌下移的盐水团影响. 东风西沙水库由于更靠近南北支分汊口, 北支倒灌发生时, 盐度最大值高于陈行水库取水口. 表层盐度小至0.05大到1.75, 有明显的半日周期性, 盐度落憩时最大、涨憩时最小, 水位上升导致的盐度变化最小; 底层盐度变化与表层相似, 但盐度增大值最大约0.8. 陈行水库取水口盐度比东风西沙水库低, 但因为落潮流受科氏力作用右偏沿南岸向海, 因而水位上升后盐度增长大于东风西沙水库取水口, 盐度值在大潮后中潮达到最大, 最大值约1, 此时盐度差值也达到最大, 约为0.3. 又因为其地理位置离南北支分汊处更远具有延迟性, 盐度最大值比涨憩略晚几小时. 青草沙水库位于南北港分汊口, 小潮后中潮盐度值最大, 表层盐度最大值在3左右, 底层为4.5左右. 水位上升后盐度增大, 不利于青草沙水库取水口取水.

      图  10  东风西沙水库、陈行水库、青草沙水库取水口表层(左)和底层(右)盐度随时间的变化

      Figure 10.  Variation in salinity on the surface (right) and bottom (left) at the water intake of Dongfengxisha Reservoir, Chenhang Reservoir, and Qingcaosha Reservoir

    • 本文采用三维数值模式ECOM-si, 数值模拟长江河口水位上升30 cm前后流场和盐度的变化, 从动力因子上分析水位上升对长江河口的影响. 该模式包括潮流、径流、风应力等多种动力因子. 选取多年1月、2月平均径流作为上游开边界条件和冬季平均5 m/s北风计算海表面风应力. 分别从流场和盐度的分布及其变化定量分析了潮流、径流、风应力对长江河口水位上升的响应, 分析了水位上升对东风西沙水库、陈行水库和青草沙水库取水口取水的影响.

      长江河口水位上升30 cm, 各河道垂向横截面面积增大, 径流向海流速整体减小. 相反, 潮流随水深增深增强. 风生流加强, 使得北支水体、盐分倒灌加剧和南北港水平风生环流加剧. 长江河口大、小潮期间, 涨落潮整体流速随水深增加而增大, 其中北港的流场分布受水深变化影响最大.

      长江河口水位上升, 进入北支的径流量增加, 北支上段盐度有不同程度的降低, 但潮流和风生流受水位上升影响更大, 总体上北支受外海盐水入侵增强. 南支中段盐度主要来源于随北支倒灌向下的高盐水团, 水深增加的情况下, 盐度无显著变化. 北港中下段在水位上升之后, 大、小潮盐度变化相似, 盐度增加最为明显, 且北港北汊由于水深增加盐水入侵尤为显著. 南北槽中上段在大潮期间盐度略微增加, 口门处滩地由于水位增加, 非线性效应减弱, 盐度下降. 水位上升后南支水源地3个水库取水口盐度均有所上升, 减少了可取水时间, 不利于水库取水.

参考文献 (16)

目录

    /

    返回文章
    返回