区域饱和-非饱和多孔介质溶质运移简化模型.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 区域饱和非饱和多孔介质的溶质运移简化模型 林琳 杨金忠 史良胜 周发超 武汉大学水资源与水电工程科学国家重点实验室 湖北武汉 中国科学院南京土壤研究所 江苏南京 摘要 在迭代求解区域饱和 非饱和拟三维渗流模型的基础上 建立区域溶质运移简化模型 该模型将溶质运动在 非饱和带简化为一维垂直运动 饱和带简化为二维剖面运动 在非饱和带与饱和带分别形成对流 弥散方程 并通 过维持近饱和带溶质质量守恒组成统一的系统 采用有限单元法数值处理技巧 在同一个矩阵方程中求解饱和与 非饱和含水层的浓度变化 该模型无需迭代计算 大大减少工作量 为求解区域溶质运移问题提供了新方法 算 例验证表明 区域溶质运移简化模型技术可行 结果可信 关键词 区域溶质运动 简化模型 饱和 非饱和 多孔介质 有限元法 中图分类号 文献标识码 研究背景 多孔介质中溶质运移理论及模型的应用研究对于环境岩土工程 改良盐渍土和防止土壤次生盐渍 化等农业生产活动 精确地计算和预测土壤中养分和污染物运移规律 水资源开发利用以及控制地下水 污染等都有着重大的现实意义 应用数学模型模拟可溶性污染物在地下含水层中迁移时 需要将 水流模型和水质模型耦合起来 求出污染物浓度的时空变化规律 以此来预测地下水污染的瞬时动态与 扩展范围 为制定合理有效的地下水污染防治措施 选择最佳治理方案提供科学依据 二维饱和非饱和水流和溶质运移软件 可以处理各种边界条件 应用非常广泛 是一个综合性地下水模拟软件 其中的 模块用来模拟饱和 非饱和流环境下水流和溶质运移的 三维有限元耦合模型 是基于三维 有限元的地下水模拟可视化软件包 能够解决完 全瞬时 半瞬时 稳态的饱和非饱和地下水流动与溶质运移模拟 完全三维饱和非饱和地下水流和 溶质运移数学模型可以得出完备的模拟结果 对于区域问题 模型缺点是需要剖分大量的网格 花费巨 大的时间 很难满足现实中人们对计算效率的要求 分别建立地下水 土壤水溶质运移子模型 然后将 二者耦合 实际是准三维模型构造 是解决区域水流和溶质运移问题的理想方案 采用有限 差方法建立准三维饱和非饱和水流模型 在剖面二维情况下进行了算例验证 计算速度已较完全二维 模型有所提高 溶质模型没有进行简化 仍然采用完全二维形式 从模拟产汇流过程出发建立 模型 对土壤水模型和地下水模型在变动潜水面上的耦合算法进行了深入的讨论 但是没有提 到溶质模型在整个土壤剖面的耦合问题 总之 区域饱和非饱和水流和溶质运移简化数学模型的构建 以及耦合算法的研究必须得到重视并提出新的思路 考虑区域地下水流运动特征 区域饱和非饱和拟三维地下水渗流模型认为非饱和带渗流以垂向运 动形式为主 建立一维水流方程 饱和带则主要是横向运动 沿含水层厚度平均得到水平二维方程 将非 收稿日期 基金项目 国家自然科学研究基金 作者简介 林琳 女 河南新乡人 博士生 主要从事地下水资源与水环境研究 饱和含水层的水量变化作为饱和含水层的垂向补排量 使二者成为互相联系的子系统 本文在耦合 求解上述模型的基础上 以溶质运移理论为基础 建立区域溶质运移简化模型 在非饱和含水层忽略水 平水流运动 从而忽略横向溶质运移 建立垂向一维对流弥散方程 而在饱和含水层仍然建立剖面二维 对流弥散方程 饱和及非饱和含水层之间通过加入相应的溶质通量项及溶质质量变化项达到维持质量 和能量守恒之目的 从而将二者统一在一个迦辽金有限单元方程里进行数值求解 避免相应拟三维水流 模型需迭代求解的复杂过程 节省计算时间 减少计算工作量 提高计算效率 区域溶质运移数学模型 基本模型 根据 定律和质量守恒原理 导出剖面二维对流弥散基本微