基于GIS的地下水流数值模拟参数自动提取研究.pdfVIP

中国地理信息系统协会理论与方法学术会议(2005)论文集 基于GIS的地下水流数值模拟参数自动提取研究木 陈锁忠 徐网谷 张磊 (南京师范大学地理信息科学江苏省重点实验室，南京，210097) 摘要：本文根据地下水流数值模拟有关参数的空间分布特点，将其划分为点状、线状与面状空间分布参数 三种类型。利用GIS的空间分析功能，研究基于GIS的点状、线状与面状空间分布参数自动提取的技术路 线。通过点、线、面三者之间空间位置关系的拓扑分析，确定地下水开采井、水位动态监测井、线状与面 状水系、模拟区域边界、面状参数分布区与模型窄间离散网格之间的空间位置关系。经过矢量数据与栅格 数据互相转换，实现地F水流数值模拟有关参数的自动提取。结合地F水流有限差分数值模拟有关参数数 据文件的结构，将基于GIS提取的栅格结构数据文件自动转换成模型所需结构的数据文件，达到模型参数 数据文件的自动组织与提高地下水流数值模拟的效率日的。 关键词：地下水流；数值模拟；模型参数；自动提取；GIS 1引言 计算机技术的推广应用，给区域地下水运动数值模拟带来了发展契机。以20世纪60年代为 界线，60年代以前地下水流模拟主要采用解析解方法；60年代以来，地下水流的模拟方法逐步 由解析解转化为数值解，目前数值解法已成为地下水流模拟的主流【2¨31。地下水流数值模拟常采 用有限著与有限元两种数值解法，20世纪80年以前，地下水流数值模拟主要基于DOS操作系 统进行，模拟模型所需的计算参数需要水文地质技术人员手．J：准备，模型参数则以数据文件的形 式赋值。由于地下水流数值模拟涉及的参数较多，模拟区域的范围较广，手工获取计算参数与赋 值占用大量的模拟时间，导致地下水流数值模拟的周期较长，模拟结果的时效性较差。20世纪 80年代以来，为了改进DOS操作系统下地下水流模拟的不足，美国、加拿大、德国与荷兰等国 家的水文地质研究人员相继开发J，地下水流模拟系统，例如ModFlow、FeFlow、GMS等地下水 流模拟系统【4J【51。地下水流模拟系统的开发大大提高了模型的自动化程度，从模型的时空离散， 到计算参数的赋值，以及模拟结果表达等方法都有了较大的改进，提高了地下水流模拟的效率， 缩短J，模拟周期。但是，已开发的地下水流模拟系统参数赋值主要通过两种方法实现，一是对剖 分单元或节点统一赋值，二是逐个单元或节点赋值，仍然是借用鼠标或键盘的手工赋值方法，而 且对各个计算单元或节点赋值时，需要参照纸质的参数分布图件／4‘能赋值，在计算单元或节点与 计算参数分区的空间拓扑关系定位上存在问题，导致系统不能自动提取的有关计算参数值，形成 模型所需的数据文件。本文主要针对上述地下水流模拟系统中参数赋值存在的问题，研究地下水 流有限差分数值模拟中计算参数的快速自动提取方法，旨在改进目前地下水流有限筹分数值模拟 中参数赋值存在的问题，提高模拟效率。文中的地下水流可视化模拟系统基于北京超图公司开发 GIS平台自主开发，适用于孔隙地下水流系统的数值模拟。 的SuperMap 2计算参数空间分布类型 地下水流模拟模型中涉及到含水层基本条件、初始条件、边界条件、地下水开采和地下水位 动态等信息，它们主要包括：含水层的顶、底板高程，初始水位、边界类型、边界水位、边界流 基金项目：国家自然科学基金资助 作者简介：陈锁忠，1964．10，男，江苏丹阳，博士，教授。主要从事水文地质与科学计算可视化研究，发 学重点实验室，电话：025802 221 中国地理信息系统协会理论与方法学术会议(2005)论文集 量、河流与水库的水位、河床与湖瓮的渗透系数、河床与湖盆的高程，地下水开采最、动态监测 井水位，含水层的渗透系数、导水系数、承雎含水层的贮水系数、潜水含水层的给水度、降水入 渗系数、潜水极限蒸发深度、弱透水层的垂向渗透系数等。按各个参数的空间几何分布特征分类， 上述地下水流模拟计算参数可分为三类，即点状分布参数、线状分布参数和面状分布参数，各计 算参数的具体归属如表1所示。 表1 地下水流模拟参数空间几何分布特征分类一览表 参数类型 计算参数 点状分布参数 地下水开采量或注水量