山东烟台夹河中, 下游地区海水入侵三维水质数值模拟研究.pdfVIP

第 卷第 期 1 ! 地学前缘（中国地质大学，北京） 9:M #1 N: # ! 年 月 ! % （ ， ） )=5 # ! @=5AB CD734D3 E5:4A735; +B74= F47G35;7A6 :H I3:;D734D3; J37K74L 山东烟台夹河中、下游地区海水入侵 三维水质数值模拟研究 ! ! 成建梅 ，陈崇希 ，吉孟瑞 ，孙桂明 （ 中国地质大学 环境地质研究所，湖北 武汉 ； 山东省第三地质队 地质工程勘察院，山东烟台 ） !# $%’$ # ($ 摘 要：建立了三维变密度对流 弥散水质数学模型来研究山东省烟台夹河中、下游地区咸淡水界面 的运移规律。以四面体为基本离散单元，推导出三维海水入侵变密度水质模型求解的数值方法，其中 水流方程求解时运用了迦辽金有限单元法。溶质运移方程求解时运用了欧拉 拉格朗日混合方法，将 对流项与弥散项分离，用传统迦辽金有限元方法求解弥散项；采用自适应)*+,))*+ 法求解对流项， 以消除人工过量和数值弥散。根据地下水的潮汐效应观测信息，确定了含水系统的海底延伸边界；利 用该地区地下水水头及水质长观资料识别了模型的水文地质参数，探讨了夹河地区海水入侵的原因： 认为夹河下游地区滨海地带地下水过量开采是造成烟台地区海水入侵的主要原因。此外，海水随潮 定期地倒灌进入夹河，通过局部岩性天窗侵入淡水含水层加剧了沿夹河河床两侧地下水的咸化。同 时还预测了几种情况下地下水的水质演化趋势，为防止和减轻夹河地区海水入侵提供合理、科学的依 据。 关键词： 海水入侵；三维变密度模型；欧拉 拉格朗日混合法；模型预测 中图分类号： 文献标识码： 文章编号： （ ） -(( . !/ %! ! ! !’0 ( ［］ 目前，在世界范围内许多沿海地区遭受海水入 单层二维模型。26=:54 提出了与密度相依赖的 侵的危害。在中国，至少有!! 个城市和地区发生了 地下水流方程和运移方程，建立了滨海多层含水层 ［］ 不同程度的海水入侵。滨海地区海水入侵常引发水 中海水入侵过渡带三维有限元模型。薜禹群% 建立 质恶化和土地盐渍化等环境问题，大大地约束沿海 了三维特征有限元模型，研究了山东龙口一滨海含 ［］ 地区经济的发展。对已退化的环境进行治理，则需 水层中的海水入侵。吴吉春$ 建立了反映水 岩阳 要花费更多的代价。通常，研究海水入侵问题时运 离子交换作用的海水入侵数学模型。 用两种模型：突变界面模型和过渡带模型。早期，许 多学者将突变界面模型用于海水入侵平面及剖面二 ! 数学模型 维模拟研究。事实上，淡水与海水是可混溶的，它们 之间存在一个咸淡水过渡带，且多数情况下过渡带 过渡带海水入侵的数学模型中必须有两个偏微