地下水数值模拟 .pdfVIP

地下水数值模拟

一、地下水数值模拟简介地下水数值模拟是各种数值方法在水文地质计算中的应用，是目前定量研究地下水

水资源量的重要手段。20世纪60年代以来，随着计算机技术的迅速发展，数值方法作为一种求解近似解

的方法被广泛用于地下水水位预报和资源评价中。数值方法是采用离散化的方法来求解数学模型，从而得到

研究区域内有限个离散点上的未知函数值。离散化的方法是将研究区域划分成为若干个较小的子区域或称为

单元，即化整为零，这些单元的集合体代表的研究区域，即又积零为整。虽然所得解为数值解(即是数值的集

合，是数学模型的近似解)，但是只要将单元大小和时段长短划分得当，即对空间步长和时间步长取值合适，

计算所得的数值解便可较好的逼近实际情况而满足计算精度的要求。由于数值方法可以较好的反映复杂条件

下的地下水流状态，具有较高的仿真度，因此在理论和实际应用方面都得到了较快的发展。

数值法求解地下水流数学模型的基本步骤如下：

(1)将研究区域按照某种规则进行剖分或称离散化。剖分的原则和剖分后形成的子区域形状取决于所采

用的数值方法，从而将研究区域划分为若干个子区域单元。对于非稳定流问题，尚需将计算时间也进行离散

化，即将计算时间离散为若干个时段。

(2)将每个小单元作为地下水的小均衡域，并定义特征点上的各种物理量。(3)建立某一个时

段内结点之间制约各种物理量的关系式，关系式一般表达

为代数方程。

(4)利用初始条件和边界条件(即初边值问题)，建立在某一个划分时段内边界结点与内部结点的关

系式。

(5)求解上述(3)、(4)所构成的代数方程组，就可求得某一计算时刻，研究区域上各离散点的水

位H值，其集合{H}即是渗流区域上某一时刻地下水水位H的近似解，单元剖分的越小，{H}的仿真度就

越高。

(

(6)重复3)〜(5)，可计算下一时刻的水头{H}集合值

由于建立代数方程组的方法不同，也就产生了各种不同的离散化方法，即不同的数值方法。目

前，地下水流计算常用的数值方法有有限差分法（FDM）、有限单元法（FEM）、边界元法（BEM）、

有限分析法（FAM）等。

三、国内外研究进展

1、国外地下水数值模拟研究进展

目前，国外该领域的研究主要针对数值模拟法的薄弱环节，提出新的思维方法，采用新的数学工具，分

析不同尺度下的变化情况，合理地描述地下水系统中大量的不确定性和模糊因素。

（1）该领域科学家在地下水系统数值模拟的工作程序、步骤方面达成了一致，强调对水文地质条件

合理概化的重要性，并深入探讨尺度转换问题和量化不确定因素问题。

（2）国外开发了许多功能多样的地下水系统数值模拟软件，以其模块化、可视化、交互性、

求解方法多样化等特点得到广泛的使用，尤其是MODFLOV，据美国地质调查局统计，MODFLOW几

乎占地下水系统数值模拟软件总应用次数的一半，这些年其功能更是不断完善。地理信息系统（GIS）

与地下水模型的整合强化了数据的输入、传递、方案调整和空间分析等。遥感（RS）提供了判断地

质边界、地貌单元和估算地表蒸发等的工具。地下水系统数值模拟模型与相关领域模型的耦合更扩展了其发展

空间，可以解决更多的实际问题。

2、国内地下水数值模拟研究进展

我国自1973年以来在地下水的数值模拟方面发展很快，它的应用已遍及与地下水有关的各个领域和

各个产业部门。高校、科研院所与生产部门相结合，已运用数值模拟解决了很多国民经济建设中急需解决的各

类问题，其中包括水资源评价问题；地下水污染问题，水岩作用和生物降解作用的模拟；非饱和带水分和盐分

运移问题；海水入侵、高浓度咸水/卤水入侵问题；热量运移和含水层贮能问题；地下水管理与合理开发、井

渠合理布局和渠道渗漏问题；地下水-地表水联合评价

调度问题；地面沉降问题；参数的确定问题。所建立的模型囊括了包括识别模型、预报模型以及管理模型