非饱和-饱和区域中重金属污染物运移的数值模拟.pdfVIP

非饱和-饱和区域中重金属污染物运移的数值模拟1 梁冰，姜利国 辽宁工程技术大学，辽宁阜新 (123000) E-mail：lbwqx@163.com 摘 要：为了定量整体描述非饱和-饱和区域内水流及溶质迁移过程，在基于非饱和流动理论 及溶质运移理论基础上，建立了非饱和-饱和区域内水流方程与重金属污染物运移方程耦合 的数学模型，采用V-G 模型来描述土壤水分特征函数，用Lagrange-Quadratic 方法进行求解 并给出模拟结果；分析结果表明，在非饱和-饱和区域内的压力水头、饱和度以及污染物阻 滞因子等参数均随含水率的变化而改变，这是区别于饱和流的主要特点。 关键词：饱和-非饱和流，重金属污染物，V-G 模型 1. 引言 随着工农业生产的迅速发展和社会人口的不断增长，土壤污染问题的日趋严重，各种施 加在地表的污染物随入渗的水经过非饱和带进入地下水中，因此，非饱和带中污染物迁移问 题己引起人们的高度重视。由于土壤一旦被重金属所污染，不仅影响土质，而且关系到植物、 动物甚至人类的健康，而且很难彻底消除。由于重金属污染物具有多源性、隐蔽性、一定程 度上的长距离传输性和污染后果的严重性，因此土壤中重金属污染问题受到各国的广泛关 注。据有关文献记载，最早开始针对土壤介质中重金属污染的研究始于上世纪 80 年代。当 时美、英等国通过大量的室内及野外土柱试验，在确定了非饱和带垂向一维弥散系数和衰减 系数后，逐步开始研究重金属离子在非饱和带的迁移转化规律，考虑土壤液相和固相浓度的 分配系数，并借助于 Henry ，Freundlich 和 Langmuir 的等温吸附模式来表示液相和固相浓度 吸附和解吸问题。在数学模型求解方面也在不断发展，由非饱和带的简单解析解发展到在非 饱和-饱和带区域中考虑复杂因素的数值解，求解的初始条件和边界条件也在不断改进，使 之更加接近于污染物迁移的实际情况。国内学者对重金属污染物在非饱和土壤中迁移转化也 进行了大量的研究工作，如朱济成（1977）、吴岳（1981）、聂永丰（1987）等针对氮、磷、 重金属和有机农药污染方面进行了试验和数值模拟研究[1] 。但是由于土壤及地下水污染具有 隐蔽性以及溶质运移具有很多不确定性，而且人们对于溶质在地下水中的弥散机理和各种化 学反应过程都了解不够，已有的研究成果对于污染物在非饱和-饱和土壤中迁移转化的整体 模拟还不多，大都示对这一系统的分散、零碎的研究，并没有将这一过程联系起来综合研究 [2] 。为此本文在已有的理论基础上，力求用数值方法模拟出重金属污染物在非饱和-饱和带 土壤中迁移转化规律。 2. 重金属污染物迁移转化特征 重金属在进入地下环境系统后所发生的迁移转化过程是一个复杂的物理、化学和生物过 程。因此，在研究非饱和-饱和区域内重金属污染物的迁移转化规律时，必须考虑各个过程 及其影响因素，其迁移转化的特征主要分以下几种： （1）物理迁移 重金属污染物随水流进入非饱和-饱和区域内以后，会随水流的运动而发生扩散、对流 1本课题得到教育部博士点基金项目“煤矿开采对阜新水环境污染的机理及控制研究”（20050147006 ）的 资助。 - 1 - 等作用，从而被机械搬运。 （2 ）物理化学迁移和化学迁移 地下环境系统中的重金属污染物与土壤无机胶体结合，发生非专性吸附或专性吸附；或 被土壤中有机胶体络合或整合，或者由有机胶体表面吸附；另外，重金属化合物的溶解和沉 淀作用，是土壤环境中重金属元素化学迁移的主要形式。 （3 ）生物迁移 土壤环境中重金属的生物迁移，主要是指植物通过根系从土壤重吸收某些化学形态的重 金属，并在植物体内积累起来。另外，土壤微生物的吸收以及土壤动物啃噬重金属含量较高 [3] 的表土，也是重金属发生生物迁移的一种途径 。 3. 数学模型的建立 地下水流方程 在分析非饱和-饱和流时以压力水头