枯季三峡库区泄流对缓解长江口盐水入侵数值模拟研究.docVIP

枯季三峡库区泄流对缓解长江口盐水入侵数值模拟研究 　　摘要：为研究枯季三峡库区泄流后对长江河口咸潮入侵的影响，利用三维数值模式ECOM-si分别模拟了枯季1、2月份径流量为大通站实测多年平均状况下和三峡库区受上游水库群调蓄后保持天然入库量的情况下泄流后长江河口的盐水入侵特征，模拟结果表明：①三峡库区泄流对缓解长江口咸潮入侵有一定的作用，表现在径流量增加，咸潮入侵变弱；②上游放水后在大潮涨憩期间对北支盐水入侵影响不大，但使得南支北港和南港部分地区在大潮涨憩期间出现了淡水；③从等盐度线垂直分布方向上看，泄流后大小潮涨憩时刻的等盐度线较之前相比都向海洋方向移动；④从取水时间上来看，三个水库在泄流后不宜取水天数和取水时间都发生了变化，不宜取水时间分别缩短了1～2 d，且青草沙水库盐度的日变化幅度也有所缩小。 　　关键词：咸潮入侵；数值模型；ECOM模式；长江河口 　　中图分类号：TU991 　　文献标识码：A 文章编号2017 　　1 引言 　　河口咸潮入侵是河口普遍现象，受潮汐、径流、风应力、混合、地形和口外环流等的综合作用，动力过程复杂，又与河口环流、泥沙絮凝和生态环境等密切相关[1]。长江河口为特大型多级分汊河口，受多种动力因子的作用，咸潮入侵还受北支倒灌的影响。开展河口咸潮入侵的研究，对丰富河口海岸学理论具有重要科学意义。 　　长江河口地区经济发达，人口密集，对淡水资源的需求量大。2010年后，青草沙水库建成，使得长江河口水源地的青草沙水库、陈行水库和东风西沙水库供水达到上海市总水量的80%以上[2]。从长江河口取水会遇到河口咸潮入侵问题，开展长江河口咸潮入侵的研究?λ?源地水库运行和保障上海用水安全具有重要的应用意义。笔者利用数值模拟的手段，对三峡泄流压咸做了定量分析，对应对长江口枯季咸潮入侵有一定的借鉴意义。 　　2 ECOM数值模式概述 　　本研究采用的三维数值模式ECOM-si（semi-implicit estuarine，coastal，and ocean model）是基于POM模式（Princeton Ocean Model；Blumberg and Mellor，1987）发展起来的，已得到了广泛的应用[3，4]。模式可设置包括风、波浪、潮流??径流等各种动力因子。 　　ECOM-si模式中海洋原始控制方程组主要包括：动量方程、连续方程、温度方程、盐度和密度方程。在不可压缩流体、Boussinesq和静力近似条件下，给出非正交坐标系下河口海岸海洋控制方程组。引入水平非正交曲线和垂向σ坐标系[5]，ξ=ξ（x，y），η=η（x，y），σ=z-ξH+ξ，其中x，y和z分别为向东、向北和向上的笛卡尔坐标轴，垂向坐标σ从海底-1（z=-h）变至海表面0（z=ξ），ξ为海表面波动，H为总水深。海洋控制方程组如下[6]。 　　3 数值模式的设置 　　模式使用水平曲线非正交网格（图1），范围包括整个长江河口、杭州湾，上游开边界设在枯季潮区界大通水文站，这样可直接采用大通水文站的实测径流量资料来给出模式的径流边界条件。模式网格较好地拟合了长江河口的岸线，并主要对南北支分汊口以及北槽深水航道工程区域的网格进行了局部加密。水深采用2010年实测资料。 　　研究先设计一个控制数值试验，模式计算时段从2011年12月1日开始计算，至2012年 2 月底结束，输出2月计算结果作分析和比较。12月至翌年1月和2月径流量取大通1950 年以来月平均值，分别为13600、11300 和11650 m3/s。风场取NCEP多年月平均风场。外海开边界考虑16个主要分潮（ M2，S2，N2，K2，K1，O1，P1，Q1，U2，V2，T2，L2，2N2，J1，M1，OO1） ，由各分潮调和常数合成得到，资料从全球潮汐数值模式 NAOTIDE 中计算结果得到。相关研究表明：在枯季，由于受上游电站消落补水影响，三峡入库流量较未受上游电站调蓄影响偏多约2000 m3 /s[7]。因此，设计对照试验将1、2两月大通水文站月平均径流量的基础上加2000 m3 /s来表征三峡库区泄流后大通站的径流量，从而来研究枯季三峡库区放水对长江口咸潮入侵的影响。由于篇幅原因的限制，模型的验证结果请参考文献[8]。 　　4 结果分析 　　为便于分析咸潮入侵沿程变化，给出北支、南支-北港、南支-南港纵向断面的分布（图2）。 　　4.1 径流量为多年月平均的情况 　　在枯季径流量为1、2月多年平均的情况下，大潮涨憩时刻长江口表层盐度分布见（图3a），北支盐水倒灌进入南支，使得南支盐度从上段到下段呈现出“高―低―高”的变化趋势。陈行水库和青草沙水库盐度低于国家饮用水盐度标准0.45（GB5749-85），南