地下水模拟软件MODFLOW的特点及其与GIS的整合.pdf

地下水模拟软件MODFLOW的特点 及其与GIS的整合 王仕琴邵景力崔亚莉刘明柱 (中国地质大学水资源与环境学院，北京100083) 摘要：本文从地下水流数值模拟的软件出发，介绍了MODFLOW三个版本的发展过程，分析了最新版本 和源汇项子程序包，其计算处理功能从水文地质过程的机理上也更加符合实际水文地质条件。文中重点阐述了根据 程序软件输入输出数据的格式进行软件前后处理的工作流程。MODFLOW程序软件运行时其输入文件具有固定的格 式，通过分析各子程序包输入输出文件数据的物理意义，可以为程序软件与GIS的整合找到良好的数据接口，从而 为我国地下水资源评价信息化集成提供可靠依据。由此指出MODFLOW程序软件在水资源评价和管理的信息化进程 中具有重要作用。 关键词：地下水MODFLOW子程序包数据结构GIS 1引省 现代地下水研究中，地下水数值模拟已成为重要的工具。随着计算机、数值模拟和GIS技术的 发展，地下水数值模型发展和应用日臻完善，同时为了用户使用的方便，基于不同方法的数值模拟 软件被相继开发并得到广泛的使用。目前常用的地下水模型软件有基于有限单元法的FEMFLOW…，有 基于限差分的GMS“l、，Visual MODFLOW。5’，PMWIN等软件。其中，后者大都是在地下水数值模拟程序 MODFLOW的基础上开发的，是目前在国际上最流行的地下水流模拟程序，被广泛应用于水文地质、 环境地质、生态规划以及水资源利用等各个行业部门中。本文将重点介绍MODFLOW一2000版本的发展 过程、不同模块程序设计特点以及它在现代水资源信息集成中发挥的重要作用。 2 MODFLOW水流数值模拟软件的发展过程 Three—DimensionalFinite—differenceGround～waterFlowModel MODFLOW是英文Modular (模 块化三维有限差分地下水流动模型)的简称，它是一套用有限差分的方法模拟孔隙介质中地下水流动 以用来解决实际中大量的地下水流问题，程序更有实用性。第三个版本为MODFLOW-96，一方面对以前 的用于解地下水流方程的软件功能进行了改进和增加，另一方面，美国地调局还开发了附属于MODFLOW 环境是在VisualStudio的编译环境下进入Fortran编辑器，其中加入了C语言的某些功能。同时， 融入了对水文地质条件机理的新理解，其输入数据的处理中增加了利用参数定义数据一项，方便了用 户，且更加有利于以后的观测过程、参数估值过程以及灵敏度分析过程的数据改动。 822 3 MODFLOW-2000程序设计结构特点 同以前的版本一样，最新版本程序设计依然具有程序设计结构的模块化、离散方法的简单化和 求解方法的多样化的特点哺1。但是在新版本中，除了用于解地下水流方程的地下水流(GWF)过程和 Fi 一个定义全局变量的过程，即Global过程，其输入数据文件包括：文件名文件(即Namele，定 义了整个模拟中激活的各个子程序包程序)、离散文件(即DiscretizationFile，包含在时间上和 空间上对模型的离散化以及描述模型结构的信息等)、乘数数组文件(MultiplierFile)和分区数 组文件(ZoneFile)。其中后两个文件是对输入矩阵形式数据的一种改进，可以利用文件中的“参 数”来定义矩阵数组中的数，这里的参数不同于模型的渗透系数、释水系数等参数，仅相当于一个 乘数因子，对于模型的每一个子程序包输入数据都可以使用参数来定义数组数据。例如：若输入某 个含水层的渗透系数数组，而这个数组是由参数定义的，这样我们可以对在分区数组文件中对渗透 系数数组分区，每一个分区里的值由乘数数组中相应的乘数因子值与该数组本身的值相乘得到。这 项功能主要利于后续的模型校正过程，可以通过改变乘数因子或分区值来改变模型的参数。而不需 要改变模型参数数组的每一个值。针对模型的校正，还增添了三个过程，分别是：观测过程(OBS)一 ——解计算值与观测值之间的方差，权重差等；灵敏度分析(SEN)