飞机场污染物运移模拟.ppt

研究区域的背景和所要解决的问题 明确研究区域的地质条件、含水层特性、补给排泄条件、边界条件 模型所需要的特性参数（如渗透率，抽水量等） 水流模型：模型的离散化（时间，空间）；源汇项（补排条件）处理；边界条件处理 溶质运移模型：溶质迁移和反应；源汇项（如污染物注入）处理；边界条件处理 模型运行：根据需求选择不同算法 （WHS,GMG,SAMG) 输出显示：根据研究的问题绘制相关曲线（如浓度穿透曲线），图幅（如等水头线图，三维流场图） 问题概念 模型建立 计算输出 稳定流与非稳定流 应力期与时间步长 应力期循环，时间循环，迭代循环 边界条件 稳定流与非稳定流 稳定流（steady flow）：指水在渗流场内运动过程中各个运动要素（水位，流速，流向等）不随时间改变的水流运动。 非稳定流（Unsteady flow）：任一点的运动要素(如流速、压强、密度等)随流程、时间等不断发生变化的水流。 应力期与时间步长及对应循环 应力期（Stress Periods）：根据外应力（Stress），如抽水量，蒸发量等确定的模拟区间。在每一个应力期内，所有的外应力保持不变。 时间步长（Time Steps): 在每个应力期内划分的若干时间段。在进行非稳定流模拟时，通过对有限差分方程组的迭代求解来计算每一个时间段的水头场。 每个模拟包括三个循环： 应力期循环 时间段循环 迭代求解循环 Modflow运行流程 边界条件 为了获得物理问题（各种微分方程）的唯一解，必须对计算域边界设定各种参数值. 注意：给定水头不是定水头 Modflow里面对于边界条件的处理 根据计算单元的性质，将模拟单元划分为：定水头计算单元；无效计算单元（Inactive）；变水头计算单元（Active Cells） 无效计算单元相当于该单元的渗透系数为0，不允许地下水通过，也称为无渗流或不透水单元（no-flow cells）。可当做水文地质边界和模拟区域内的弱透水层（不透水层）。 定水头计算单元和无效计算单元是用户预先设定的，在计算过程中保持不变。其余的为变水头计算单元，其水流场随时间和空间发生变化。 其他边界条件（如定流量边界）则等同于源汇项一起进行处理。 可设置水文地质边界和不透水边界 污染源：飞机加油站 不连续含水层 供水水井 研究区域：Waterloo 河外的飞机场附近 地貌特征：飞机加油站；市内的供水井区域以及一个不连续的隔水层。 背景资料：在过去的十年里，飞机燃料已经周期性地泄露到加油站区，而且经过自然渗透导致一部分污染物进入到上层含水层。 研究问题：渗漏的燃料污染物对市内供水井的潜在影响 550m3/d 400m3/d 潜水含水层（沙土、砾石） 弱透水层（粘土，泥沙） 承压含水层（沙土、砾石） 1e-10m/sec 2e-4m/sec 2e-4m/sec 定水头边界 定水头边界 概化实际的水文地质（污染运移）问题，使用Visual Modflow创建一个地下水流模型（溶质运移模型） 对模拟区域内进行初始条件、边界条件、源汇项、迭代计算方法进行设置。 运行查看结果并进行可视化输出显示。 案例学习要求： 1）全屏打印关键步骤结果，复制到word文档。 关键步骤：几何模型、渗流场、浓度图 2）图片要小，彩色打印，2页纸。 3）12月底交，占10%。 4）12月有几次答疑。 Visual Modflow 4.1简介 1 模型的思路和步骤 3 地下水模拟相关概念 4 运行环境及界面操作 2 算例的概念模型及学习目标 5 算例的具体操作（演示） 6 Modflow是美国地质调查局（USGS)于80年代开发的一套专门用于孔隙介质中三维有限差分地下水流数值模型。 不足：基于DOS的操作系统无交互界面，数据准备复杂，要依靠第三方软件进行结果的后处理，且软件之间无交互性。 MODFLOW MODPATH MT3D PEST 标准的可视化操作界面 强大的模型输入输出支持 三维后处理显示 4.2 PP 79 PP 435 PP 76 管理工程文件 运行已建立模型 查看模型结果 创建模型 数值引擎 导入使用USGS标准MODFLOW创建的模型 从这里找到保存的工程 从这里导入标准MODFLOW模型 屏幕截图 需要先打开模型文件 主工具栏 编辑窗口 操作的对象 所选单元 文件保存，显示偏好等 当前视图，编辑网格 同时在这里可以选择编辑任意参数 从这里进入编辑网格 从这里进入编辑井位信息 同时在这里可以选择编辑任意参数 同时在这里可以选择编辑任意参数 从这里进入编辑所列参数 同时在这里可以选择编辑任意参数 从这里进入编辑边界条件 同时在这里可以选择编辑任意参数 从这里进入编辑示踪粒子与水均衡 编辑模型时的一些辅助性工具 3D浏览器 加载底图 大小、移动 垂向放