感潮河段盐水入侵的河口水质模型 Numerical simulation of water quality in tidal river network under the action of unsteady saline water intrusion.pdfVIP

水 利 学 报 2009年5月 SHUIUXUEBAO 第40卷第5期 文章编号：0559．9350(2009)05．0520-09 感潮河段盐水入侵的河口水质模型 茅泽育，赵雪峰，江春波 (清华大学水利水Fg-r程系，北京100084) 摘要：盐潮入侵将严重影响内河水质。本文考虑强混合型河口潮汐的特点，采用同步不同格式对河口盐潮入侵问 题的流场及浓度场进行数值计算。采用欧拉一拉格朗日混合解法将对流一扩散方程中的扩散项和对流项分裂求 解，应用高精度时间序列插值的特征差分隐格式求解对流项，扩散项采用传统的有限差分法。在汉口处，将描述河 网汉口水质变化的三维对流扩散方程在汉口控制体体积积分，变量间的插值分布采用Panankar提出的幂函数，得 到新的汉口离散方程。计算与实测结果比较表明，本文提出的方法有效可行。 关键词：河口水质；海水入侵；盐度分布；感潮河网；数值计算 中图分类号：X143 文献标识码：A 1研究背景 河口感潮河段和一般内河河网最显著的区别是受潮汐影响，流量变化剧烈，并随海潮带来大量氯离 子(CL。)，导致其水质发生变化。盐水入侵大体分两种形式，对于无潮海河口，海水沿底向河121延伸，形 成稳定的盐水楔异重流，其位置视径流量大小而上下移动；对于有潮海河口，盐淡水之间以混合型为主， 盐水沿河口入侵后与淡水不同程度混合，不存在稳定的盐水楔，情况比盐水楔异重流复杂。特别当河口 潮差较大，潮流作用较强而径流作用相对较弱时，咸淡水之间存在强烈混合，因此在水平方向仍有密度 梯度，而垂直方向的密度梯度甚小。一些研究结果表明¨。31，对于强混合型的潮汐河口狭窄河网的海水 入侵，可简化为纵向一维模型。 正确确定河口盐度纵向分布较为复杂。盐度的纵向分布受制于上游径流和海洋潮汐，径流有洪枯 之分，潮汛有大小之别，组合情况甚多。盐水人侵范围、盐淡水混合程度均随潮流强弱而变。潮流速随 时间和地点而变，含盐度的分布也随之而变。流场计算的正确与否对确定盐度的纵向分布影响很大。 以往，感潮河网的流场计算大多采用水文学上的连续方程与简化的动量方程进行计算。显然，简化方法 没有完整方程组精确，其解缺乏普遍性。而且多数简化方法不能考虑邻近支流人流或分流的影响，另一 方面，由于问题的复杂性，人们常应用一定简化研究盐度纵向分布的解析方法，通过解析解揭示盐淡水 混合的内涵，增加对盐淡水混合认识。虽然解析法较简便，有一定优越性，但考虑的往往是简单的准恒 定情况，并假定河口为矩形断面。一般的河口都难以满足这些条件，对于淡水流量变化和日潮不等河口 的盐水人侵(盐淡水混合)，迄今尚无完善的解析解。近年来，人们开始研究河口不恒定海水入侵的数值 计算，如文献[1]采用有限差分法求解雷诺形式的二维水流N—S方程，研究盐水楔的入侵发展。 数值求解盐度纵向分布关键在于求解盐度守恒方程。盐度守恒方程为对流一扩散型方程，数值求 解该类型方程的主要困难是对流占优情况。对于扩散占优情况，通常的数值方法和格式，例如二阶中心 差分格式、线行插值的Galerkin有限元等，均能给出较好结果。而对于对流占优情况，通常的数值方法 收稿日期：2008-05-04 基金项目：国家973项目(2006CB403304) 作者简介：茅泽育(1962一)，男，浙江人，教授，主要从事水动力学研究。E-mail：mamyu@岫咄Ⅱ．edu．∞ ·———520·-—— 万方数据 则会给出在快速变化区的振荡型解，甚至出现负浓度等。为了消除这种不应有的振荡，人们发展了各种 不同的方法和格式，如迎风有限差分格式、迎风有限元格式、“反扩散”办法，以及采用高阶空间一时间有 限差分或有限元格式。但这些办法不是过分复杂就是效果欠佳，且常导致较大的数值扩散和耗散。 在汉口处理方面，目前各种水质数值模型一般均假定汉口各断面浓度相等。当一条浓度较大的河 流与相对洁净的河流交汇时，这个假定将给数值计算带来较大误差。河网的特殊性就在于汉口，随着汇 流方式不同，水流及含有物汇合的方向也不同，不再象单向河道那样存在一个主流方向，因此用一维对 流扩散方程描述含有物输移是不够的，而必须采用三维对流扩散方程。 2数学模