用的三维地下水污染物运移数值模型.pdfVIP

维普资讯 第28卷 第1期 东 华 理 工 学 院 学 报 VoL28 N0．1 2005年3月 JOURNAL OF EAST CHINA INSTITI胍 OF TECHNOLOGY Mar．2o05 M D一通用的三维地下水污染物运移数值模型 朱玉水 ，段存俊2 (1．福建莆田市涵江区环境保护局，福建 莆田351111；2．南京勘察工程公司，江苏南京210000) 摘 要：定量研究污染物在地下水中的运移过程通常采用数值模拟方法。MT3D是一套基于有限差分方法的污染物运移 模拟软件，近年来在国外水文地质和水环境模拟等领域的研究中已经得到较为广泛的认可。MT3D比较全面地考虑了污 染物在地下水中的对流、弥散和化学反应等过程，可以灵活处理各种复杂的源汇项和边界条件，能够准确模拟承压、无压 和越流含水层中的污染物运移过程。MT3D具有模块化的程序结构、灵活的求解方法以及全面的模拟功能，非常适合实际 问题的研究，值得在国内推广使用。 关键词：地下水；数值模拟；运移模型；M33D 中图分类号：P641．2 文献标识码 ：B 文章编号 ：1000—2251(2005J0l一026—04 水作为一种极其重要的资源，在经济和社会发 用了对流一弥散方程来描述污染物在三维地下水 展中具有重要的战略价值。由于地表水资源分布 流中的运移，即 很不均衡而且容易受到污染，地下水资源的开发、 尺oC = 利用和保护研究被提到了日益重要的地位。近几 詈(篝)一蠢()+詈c一 十年来，由于工农业生产的发展、环境污染的预测 A(C+ 和防治、新能源的开发利用、核废料的贮存等方面 的需要，对地下水运动的研究已经越来越深人。除 式中，C是溶解于水中的污染物的浓度，ML一；是 了地下水资源的定量评价外，越来越多的研究涉及 阻滞因子，无量纲；t是时间，T；x是空间坐标，L；D 到地下水中污染物的运移以及污染含水层的治理 是水动力弥散系数张量，L2T～；vi是地下水渗透流 等问题。同时，由于计算方法的长足进步和高性能 速，LT～；吼是源 (正值)或汇(负值)的单位流量， 计算机的广泛采用，在解决实际生产问题时，越来 T～；c是源或汇的浓度，ML一；0是孔隙度，无量 越多地采用了数值模拟方法(朱学愚等，1990)。近 纲；A是一阶反应速率常数，T～； 是多孔介质的 年来，Zheng(1990)开发了一套用于模拟污染物在 比重，ML～；C是吸附在介质上 的污染物浓度， 地下水中运移的三维有限差分计算程序—— MM～。式(1)的等号右端从左至右依次为弥散项、 MT3D，逐渐引起国外水文地质工作者的重视，目前 对流项、源汇项和化学反应项。 已经成为污染物运移模拟领域事实上的标准。 求解式(1)中的对流一弥散方程需要首先确定 MT3D能灵活有效地处理各种源汇项和边界条件， 地下水渗透流速。根据Darcy定律，有 适用于承压、无压以及承压一无压含水层中的污染 物运移问题。MT3D中包括多种求解算法可供选 = 一 等差 (2) 择，因此特别适合实际问题的研究。 式中，K．j是渗透系数的主轴分量，LT～；h是地下水 位，L．地下水位可以通过求解三维地下水流动方 I MT3D的基本方程和主要功能 程来得到，即 MT3D是英文 Modular3-Dimensiona