地下水模型发展与应用汇编.pptx

地下水模型发展与应用 北京地下水环境调查监测高级研修班 邵景力 中国地质大学（北京） 水资源与环境学院 2011年11月2日 主要内容 地下水模型概述 地下水数值模型建模步骤 地下水模拟软件简介 地下水模型应用 经验和体会 地下水模型概述 地下水的功能和作用 资源：人类的主要供水水源 环境：重要的环境要素：地下水生态问题；海咸水入侵；地面沉降；地下水污染；….； 如何实现地下水的合理开发利用、保持地下水的环境功能？ 给出地下水影响因素与状态的定量关系—即建立地下水模型； 以地下水数值模型为工具，给出地下水控制、合理开发利用方案、……。 1 原型：研究的对象，客观存在的真实系统。 模型：是原型的替代物，对真实系统某些特性的仿真和模拟，是实际系统或过程的简化、抽象和类比表示。 Anderson：模型是近似表示野外情况的任意装置和手段。 作用：描述或表述实际系统的特性，再现实际系统的状态和过程。 类型： 物理实体模型 概念模型：文字、框图、示意图等 数学模型 1 地下水模型概述 时间状态上 稳态模型 非稳定模型 空间分布上 集中（中）式模型（随机模型）—黑箱模型，一般用统计方法求解 分布式模型（确定性模型）：一维、二维、三维 模拟的地下水状态 水流模型（饱和、非饱和） 溶质运移模型（辨析：水文地球化学模型） 热运移模型 地面沉降模型（地下水-介质应变模型） 求解方法上 统计方法（相关分析、回归分析、时间序列分析、地质统计学、频谱分析、） 解析法（泰斯公式、纽曼公式、…..） 数值法（有限差分法、有点单元法、边界元法、……） ……… 地下水模型概述：数学模型 1 连续介质假设：典型单元体 质量守恒定律 达西定律 1 地下水模型概述：控制方程 地下水模型概述：控制方程 地下水流连续性方程： 质量守恒+达西定律 非均质各向异性介质 非均质各向同性介质或非均质各向异性主方向与坐标一致 1 地下水模型及分类：控制方程 非饱和带水分运移方程 溶质运移方程的定解问题 质量守恒+达西定律+费克定律 1 地下热水运移方程 质量守恒+能量守恒+达西定律+热运移定律 强调一下地下水流模型的重要性 地下水模型概述：控制方程 1 如果能求解已经建立起的地下水连续性方程，那么，就能知道地下水的状态随着输入、边界条件和初始条件的变化情况，用地下水系统的参数和输入来表示地下水系统的输出。 建模方法 解析法 数值模拟 有限差分法 有限单元法 ………….. 地下水模型概述：数值模型 1 数值法是一种近似解法。 它以地下水运动的微分方程的定解问题为基础，将表示水位（或其他状态变量）随时间和空间连续变化的函数离散化，求得函数在有限节点（或结点）上的近似值。 用解析法可以求出任意时间和任意点的水位值，而数值法只能求出空间上有限个节点在有限时刻的水位近似值； 只要近似值能满足精度要求，就可用于解决实际水文地质问题。 用于水文地质计算的数值法主要有有限差分法和有限单元法两种。他们都是将描述地下水运动的定解问题用不同的方式离散化，使复杂的定解问题化成简单的代数方程组，再求出有限点上不同时刻的数值解。 1 地下水模型概述：数值模型 有限差分法基本思想： 用渗流区内选定的有限个离散点的集合来代替连续的渗流区，在这些离散点上用差商来近似代替导数，将描述求解问题的偏微分方程及其定解条件化为一组以有限个未知函数在离散点上的近似值为未知量的差分方程组，然后对差分方程组进行求解，得到所求解在离散点上的近似值。 发展历史 上世纪50年代石油流动领域 上世纪60年代中期拓宽用于解决地下水流问题； 优缺点： 优点：数学上直观易懂；有相应的高效解法；有成熟的商业软件。 不足：处理不规则边界、大角度倾斜含水体、各向异性介质等复杂条件难以处理；对于溶质运移、热运移等问题，不如有限元求解精度高。 1 地下水模型概述：数值模型 有限单元法基本思想： 建立在直接求解函数近似解的基础上。用有限个单元的集合来代替渗流区，选择简单的近似函数（常用多项式差值）表示单元内部的状态（水头、浓度、温度的未知函数）分布，运用Rayleigh-Ritz法、Galerkin法或均衡法等建立单元内未知变量的表达式，最后集合单元方程形成整个渗流区的代数方程组，并求解方程组得到未知变量（水头、浓度、温度等）在节点上的值。 发展历史 20世纪60年代后期引入地下水计算中。 优缺点 优点：程序的统一性（水流、溶质运移、…）；不规则边界处理灵活；单元大小比较随意（局部加密）；理论上计算精度较高。 缺点：局部质量不守恒，有时会影响计算精度。 1 地下水模型概述：数值模型 1 地下水模型概述：数值模型 1 地下水模型概述：数值模型 数值法解地下水模型流程图 1 地下水模型概述：数值模型 数值法优点 广