第一讲：地下水流数值模拟软件介绍课件.ppt

VM 地下水流数值模拟软件(Visual MODFLOW) 中国地质大学(武汉）环境学院 2011.10 课程安排 1、熟悉VM软件的计算流程与数据准备与输入 1.1上课 1.2上机 2、数据计算与调参 2.1上课 2.2上机 3、模型预测、出图写报告 3.1上课 3.2上机 什么是数值模拟 数值模拟也叫计算机模拟。它以电子计算机为手段，通过数值计算和图像显示的方法，达到对工程问题和物理问题乃至自然界各类问题研究的目的。 地下水数值模拟的目的和意义 地下水数值模拟的基本目的是预测地下水未来动态，为水质和水量的评价提供理论依据。 任务：模拟地下水的流向及地下水水头与时间的关系。 数值模拟：机遇与挑战 提 纲 地下水数值模型建模步骤 地下水流数值模拟资料需求 MODFLOW软件 实例 1、地下水数值模型建模步骤 模拟步骤 建立概念模型 建立数学模型 数值方法及软件（编程） 参数 拟合模拟：模型校正与检验 参数敏感性分析 预测模拟 一、概念模型（模型概化） 根据详细的地形地貌、地质、水文地质、构造地质、水文地球化学、岩石矿物、水文、气象、工农业利用情况等 模拟的区域: 含水层类型: 潜水（无压）、承压、混合、多层 维数:一维、二维、三维 水流状态：稳定流/非稳定流、饱和流/非饱和流 介质状况: 均质和非均质/各向同性和各向异性 孔隙/裂隙/双重介质 流体的密度差 边界条件和初始条件 必要时需进行一系列的室内试验与野外试验, 以获取有关参数, 如渗透系数、弥散系数、分配系数、反应速率常数等。 二、数学模型 三维地下水流动问题控制方程 三、数值方法及软件（或编程） 绝大部分数学模型是无法用解析法求解的，数值化就是将数学模型转化为可解的数值模型。 软件 有限差分法: MOFLOW系列 GMS 中MODFLOW Visual MODFLOW Processing MODFLOW 有限单元法: FEFLOW 积分有限差分法: TOUGH2，TOUGH REACT 四、模型参数 含水层参数：渗透系数，弹性释水系数（重力给水度），孔隙度等 源汇项： 大气降水入渗系数（分区、数值） 蒸发排泄系数 地表水体（河流、湖泊、水库等）水位、底面高程、底面岩性特征（厚度、渗透系数等） 渠系灌溉入渗系数 人工开采（点状、面状） 边界条件 初始条件 参数的不确定性 五、模拟：模型校正（参数识别） 将模拟结果与实测结果比较,进行参数调整, 使模拟结果在给定的误差范围内与实测结果吻合。 调参过程是一个复杂而辛苦的工作, 所调整的参数必须符合模拟区的具体情况。 人机交互与自动调参相结合。尽管自动调参程序(如PEST ) , 也不能代替人的工作。 五、模拟：模型检验　 模型验证是在模型校正的基础上, 进一步调整参数, 使模拟结果与第二次实测结果吻合, 以进一步提高模型的置信度。 六、灵敏度分析 校正后的模型受参数值的时空分布、边界条件、水流状态等不确定度的影响。灵敏度分析就是为了确定不确定度对校正模型的影响程度。 七、预测 用校正的参数值进行预测, 预测时需估算未来的水流状态。 后续检查与模型的重建（完善） 后续检查在模拟研究结束数年后进行。收集新的野外数据以确定预测结果是否正确。如果模拟结果精确, 则该模型对该模拟区来说是有效的。由于场址的唯一性, 故模型只对该模拟区有效。 后续检查应在预测结束足够长的时间后进行, 以便有足够的时间发生明显的变化。 模型的再设计 一般来说, 后续检查会发现系统性能的变化, 从而导致概念模型和模型参数的修改。一般来说, 所有模拟研究都应该进行到第五步, 即校正灵敏度分析。 2.地下水流数值模拟资料需求 空间展布相关资料 地表高程、分层数据 边界位置 含水层参数：渗透系数，弹性释水系数（重力给水度），孔隙度等 源汇项： 大气降水入渗系数（分区、数值） 蒸发排泄系数 地表水体（河流、湖泊、水库等）水位、底面高程、底面岩性特征（厚度、渗透系数等） 渠系灌溉入渗系数 人工开采（点状、面状）：开采井位置、井结构、开采量动态 边界条件：边界类型、水头或流量 初始条件：统测水位 水位动态观测资料：观测孔位置、结构、水位时间变化 3、MODFLOW软件 MODFLOW的由来 MODFLOW是英文名称Modular Three-Dimensional Finite-difference Ground-water flow model (三维有限差分地下水流模型)的简称。由美国地质调查局于80年代开发出的一套专门用于孔隙介质中地下水流动数值模拟的软件。自问世以来, MODFLOW已经在