长江河口盐水入侵的原因及对策.docxVIP

长江河口盐水入侵的原因及对策 该河口是河流和海洋之间混合的过渡地带，含有丰富的物质资源。作为中国第一大河的长江在我国最大城市上海市入海, 不仅为临近沿海城市港口、航运、渔业等的发展提供了便利, 同时还为上海市的居民生活及工农业生产提供了丰富的淡水资源。但是由于长江口其特殊的地理位置及“三级分汊, 四口入海”的格局, 河口地区时常有盐水入侵现象发生, 枯季大潮期间尤为严重。特别是近年来随着长江流域重大水利工程三峡及南水北调的实施及长江口北槽深水航道治理工程的完工, 必然导致长江上游来水的季节分布及长江口北槽河势发生变化, 这都会对长江口盐度的时空分布及盐水入侵格局产生一定的影响。 我国对长江口盐水入侵的研究始于20世纪60~70年代, 经过几十年的发展, 众多学者利用实测资料分析及数值模拟等方法对长江口盐水入侵的时空分布、北支盐水倒灌及其通量、盐水入侵对重大工程及水文气象条件变化的响应等方面进行了深入的研究。 文章依据长江口枯季的水文盐度实测资料, 利用丹麦DHI公司开发的MIKE21Flow Model FM无结构网格水动力 (HD) 模块及其附属的温盐 (TS) 输移模块, 建立“长江口—杭州湾”二维水流盐度数学模型, 在验证良好的基础上, 对枯季长江口流场及盐度的时空变化特征进行详细分析, 并对北支盐水倒灌过程及倒灌盐通量进行了初步计算。 1 模型构建及网格单元 根据以往的研究, 长江口枯季大潮间盐水入侵现象最为严重, 最远可上溯至江阴附近, 同时, 长江口与杭州湾之间的水体和盐度交换对精确模拟长江口盐水入侵至关重要。因此, 模型范围为上游至江阴验潮站, 北至江苏吕四港附近, 南至浙江象山县, 东部外海开边界距离北槽深水航道约120 km的大范围海域。模型模拟时间为2002年2月28日~3月9日共10 d, 验证资料为长江口水域同步全潮水流盐度实测资料, 模型计算均以1985国家高程基准为基准面。模型验证点见图1。模型共包含45 712个网格单元, 23 977个网格节点。网格空间步长外海开边界处最大为7 050 m, 北支青龙港及北槽深水航道处最小分别约为300 m、450 m;时间输出步长为30 s, 计算时间步长由CFL数控制。具体的模型网格见图2。 1.1 控制方程式 动量方程 盐度输运方程 1.2 模型建立及边界条件 (1) MIKE21模型中空间离散采用有限体积法;时间离散可以采用一阶显式欧拉法或者二阶龙格—库塔方法, 由于模型范围较大, 为减少计算时间, 模型的时间离散采用一阶显式欧拉法求解。 (2) 底部糙率根据模型区域水深及床面形态, 经反复调试后取为0.01~0.02之间。 (3) 模型水平涡粘系数A计算中cs取为0.28;盐度水平扩散系数AH主要是由水流的紊动引起的, 模型中取AH=A。 (4) 模型的外海开边界由全球潮波模型预报潮位驱动, 上游边界采用江阴验潮站实测潮位控制;固壁边界法向流速为0, 对于潮滩较大区域, 模型中引入干湿网格判别技术;对于盐度边界条件, 南边界自西向东采用15‰~30‰线性插值, 东边界自南向北为30‰~35‰线性插值, 北边界自西向东为25‰~35‰线性插值, 上游开边界取盐度梯度为0。 (5) 模型的初始条件, 水动力模型采用冷启动方式, 初始水位、流速及流向皆为0;对于盐度输运模型, 经过9个月的反复计算后得到的盐度场作为初始盐度。 1.3 长江口管理特征 由于实测点位较多, 潮位的验证选取上游的徐六泾及北槽顶端的横沙站进行验证, 因为要计算北支倒灌盐通量, 北支涌潮现象的精确模拟对计算结果影响很大, 选取青龙港站进行验证;潮流在徐六泾附近、南北支各选一点进行大、小潮验证;盐度仅选择代表性站位进行大潮的验证。 由图3~图6可知, 潮位、潮流的计算结果与实测资料较比较吻合, 能够反映长江口涨潮历时小于落潮历时, 并且由口门至口内涨潮历时缩短、落潮历时延长的基本特征, 其中北支青龙港站的潮位模拟结果涨潮历时明显小于落潮历时, 潮位模拟结果与实测值也比较吻合, 说明该模型对北支涌潮现象的模拟比较准确。模型的盐度模拟结果与实测值基本吻合, 能够反映长江口盐度的变化趋势, 由Y4、Y5、Y8、Y9、Z5的盐度模拟结果对比以及Z1、Z6、Z9的盐度, 模型较好地模拟了长江口南支盐水入侵形势及北支的高盐水状况, 结果可以用其来分析长江口枯季期间的盐水入侵现象。 由图7可知, 模型很好地模拟了长江口流态。对整个区域来说, 杭州湾流速最大, 外海水域开阔, 流速较小。长江口水流北支由于大潮时涌潮存在, 涨潮流速特别大, 落潮流相对较小。对于南支水域, 总体来说落潮流速大于涨潮流速, 且越往上游这种现象越明显。同时, 由图7可以看出, 北槽与北港会通过长兴岛与横沙岛之间的横沙通道进行