地下水模拟技术及应用培训

1、地下水数值模拟技术与应用地下水数值模拟技术与应用第四讲第四讲 地下水流模型地下水流模型MODFLOW简介简介中科院计算所培训中科院计算所培训主要内容3.1 主要功能和主要模块3.2目前国内外主要地下水模拟软件简介3.3 Visual Modflow简介（简单操作）3.4 GMS简介3.5地下水数值模拟发展趋势3.6一点体会和经验4.1 主要功能和模块MODFLOW（Modular three-dimensional finite-difference ground-water flow model）是由美国地质调查局上世纪80年代开发的基于有限差分法的孔隙介质中三维地下水流模拟程序，现已推出M

2、ODFLOW88、96、2000和2005四个版本，已成为功能完善、扩展性强、应用最为广泛的地下水流模拟程序。主要功能模拟稳定、非稳定地下水流非均质含水层，简单的各项异性含水层二维、三维地下水流模拟承压、潜水、承压-无压含水层地下水流模拟适用条件饱和空隙介质地下水流，地下水流动满足达西定律。4主要特点：采用FORTRAN语言编程，可下载源程序，可根据需要对程序改编；采用模块化结构，MODFLOW程序可分为一个主程序和若干个高度独立的子程序（模块，modules），若干相关的子程序整合形成具有特定功能的子程序包（pakages）。使程序易于理解和修改，便于二次开发和增加新的模块和子程序包，对其功

3、能进行扩展；离散单元的简单化。采用矩形不等距网格离散，便于用户对模拟区剖分和准备输入数据，输出的计算结果也比较规范化。Upscale技术；在时间离散上，引入应力期的概念，便于模拟期内时间段的划分和时间步长的设定；求解方法多样化；输出格式的标准化和多样化；资源丰富，可在有关网站上下载源程序、参考手册、操作手册等；有较多的软件产品，广大的用户。子程序包名称子程序包名称英文英文缩写缩写功能功能基本子程序包基本子程序包BASBAS指定边界条件、时间段长度、初始条件以及输出指定边界条件、时间段长度、初始条件以及输出计算单元间渗流子程序包计算单元间渗流子程序包BCFBCF计算有限差分方程组各项的值，即网格

4、间流量和储计算有限差分方程组各项的值，即网格间流量和储存变化量存变化量外应外应力子力子程序程序包包水井子程序包水井子程序包WELWEL定义点井的滤水管位置和井流量定义点井的滤水管位置和井流量补给子程序包补给子程序包RCHRCH定义含水层的面状补排强度定义含水层的面状补排强度河流子程序包河流子程序包RIVRIV将河流处理为三类边界将河流处理为三类边界排水沟子程序包排水沟子程序包DRNDRN如果地下水位高于排水沟底板，则定义为一类边界，如果地下水位高于排水沟底板，则定义为一类边界，如果低于底板，则无补排量。如果低于底板，则无补排量。蒸发蒸腾子程序包蒸发蒸腾子程序包EVTEVT按阿维杨诺夫公式处理地

5、下水的蒸发蒸腾量按阿维杨诺夫公式处理地下水的蒸发蒸腾量通用水头子程序包通用水头子程序包GHBGHB将边界定义为三类边界条件将边界定义为三类边界条件求解求解子程子程序包序包SIPSIP子程序包子程序包SIPSIP强隐式法求解线性方程组强隐式法求解线性方程组SSORSSOR子程序包子程序包SSORSSOR超松弛迭代法求解线性方程组超松弛迭代法求解线性方程组PCGPCG子程序包子程序包PCGPCG预调共轭梯度法求解线性方程组预调共轭梯度法求解线性方程组MODFLOW主要子程序包功能表主要子程序包功能表名称名称功能功能备注备注ETSETS蒸散发子程序包，允许用户定于蒸散发蒸散发子程序包，允许用户定于蒸

6、散发-地下水位关系。地下水位关系。已整合到已整合到MODFLOW-2000和和MODFLOW-2005中。中。STRSTR河流河流-含水层交换量子程序包，河流量计算考虑了沿河流路径与地下水交换。含水层交换量子程序包，河流量计算考虑了沿河流路径与地下水交换。同上同上SFR2SFR2河流入渗子程序包，可模拟河流非饱和入渗补给含水层。河流入渗子程序包，可模拟河流非饱和入渗补给含水层。同上同上HUF2HUF2计算水文地质单元渗流子程序包，可将实际介质的水文地质属性（渗透性）赋到计算水文地质单元渗流子程序包，可将实际介质的水文地质属性（渗透性）赋到个网格上。个网格上。同上同上IBS1/2IBS1/2基于

7、基于MODFLOW开发的有粘土夹层的含水层形变子程序包，可计算由于地下水开发的有粘土夹层的含水层形变子程序包，可计算由于地下水储存量变化引起的含水层形变量。储存量变化引起的含水层形变量。同上同上SUBSUB在在IBS1/2基础上改进的地面沉降和含水层压缩子程序包，可单独计算出粘性土索基础上改进的地面沉降和含水层压缩子程序包，可单独计算出粘性土索性压缩和含水层弹性压缩量。性压缩和含水层弹性压缩量。同上同上UZFUZF模拟非饱和带模拟非饱和带-饱和带垂向水流子程序包。饱和带垂向水流子程序包。源程序已整合到源程序已整合到MODFLOW-2005中。中。LGRLGR局部网格加密（局部网格加密（LGR）

8、功能，允许用户在网格更大更粗糙的母模型中模拟更高）功能，允许用户在网格更大更粗糙的母模型中模拟更高精度的局部网格（称为子模型）中的水流。精度的局部网格（称为子模型）中的水流。源程序已整合到源程序已整合到MODFLOW-2005中。中。MODPATHMODPATH地下水中质点运移路径的三维质点示踪模型，可进行正向示踪和反向示踪。地下水中质点运移路径的三维质点示踪模型，可进行正向示踪和反向示踪。/softwar/software/MODPATH/e/MODPATH/MT3DMSMT3DMS基于基于MODFLOW开发的考

9、虑对流、弥散、源汇项、化学作用的多组分地下水溶开发的考虑对流、弥散、源汇项、化学作用的多组分地下水溶质运移模型质运移模型/mt3d//mt3d/RT3DRT3D在在MT3D基础上模拟地下水中多组分反应，适合于模拟自然衰减和生物恢复。基础上模拟地下水中多组分反应，适合于模拟自然衰减和生物恢复。/rt3/rt3d.htmd.htmSEAM3DSEAM3D在在MT3D模型基础上开发的碳氢化合物降解模型，可模拟多达模型基

10、础上开发的碳氢化合物降解模型，可模拟多达27种物质的运移和种物质的运移和相互作用。相互作用。MOC3DMOC3D基于特征值法的三维溶质运移模型，可考虑对流、弥散、混合、线性吸附和放射基于特征值法的三维溶质运移模型，可考虑对流、弥散、混合、线性吸附和放射性衰减等作用。性衰减等作用。MODFLOWMODFLOW- -SURFACTSURFACT饱和饱和-非饱和水流和污染质运移模型非饱和水流和污染质运移模型http:/www.scientificsoftwarehttp:/ 国内外流行的地下水模型软件名称主要功能特点开发单位、下载或购买网址GMS(Groundwater Modeling Syste

11、m，windows)整合了MODFLOW、MODPATH、MT3D、FEMWATER、RT3D、SEEP2D、SEAM3D、UTCHEM、PEST、UCODE、NUFT等模型和程序包，等可进行水流、溶质运移、反应运移模拟；建立三维地层实体，进行钻孔数据管理、二维(三维)地质统计，与ARCGIS有良好的接口。使用界面友好，前、后处理功能及三维可视化效果优良，目前已成为国际上最受欢迎的地下水模拟软件。适用于孔隙介质三维地下水模拟，与GIS接口好，是目前国内最常用的软件地下水流和溶质运移模拟软件，但软件较复杂，价格偏贵。由美国Brigham Young University的环境模型研究实验室开发，

12、网址：/gms.htm。分销商为：http:/www.ems- Modflow（windows）综合已有的MOD-FLOW、MODPATH、MT3D、RT3D和WinPEST 等地下水模型而开发的可视化地下水模拟软件，可进行三维水流模拟、溶质运移模拟和反应运移模拟。合理的菜单结构、友好的界面和功能强大的可视化特征和极好的软件支撑使之成为许多地下水模拟专业人员选择的对象。适用于孔隙介质三维地下水模拟，简单易学，价格相对便宜，是目前国内最流行的地下水流和溶质运移模拟软件之一。加拿大Waterloo水文地质公司开发研制http:/www.Waterloo

13、www. S 有多种理论模式进行选择。运用达西、泊松以及manning-strickler理论的离散单元分析。与GIS接口好，剖分灵活，可进行复杂的（各向异性）地下水流、溶质和热运移模拟。专业性强，价格偏贵。是目前国内常用的地下水模型软件之一。Feflow是由德国Wasy水资源规划系统研究所研制开发，网址：http:/ www.wasy.de/english/produkte/feflow/ index.html#。Groundwater Vistas（windows）基于MODFLOW开发的在功能上类似于GMS的地下水模型软件，功能强大，用户界面优化，可视化程度高，可用蒙特卡洛方法结合mod

14、flow、MT3D进行参数敏感性分析、地下水风险评价等。结合么特卡罗方法的随机地下水数值模拟美国Environmental Simulations, Inc.公司开发：http:/www.ground - .为汉字界面和帮助，有利于国内专业人员学习和使用。中国水利水电研究院开发IGW（Interactive Ground Water， windows）交互式二维、三维地下水流和溶质运移数值模拟软件，可设定确定性、蒙特卡洛和条件蒙特卡洛随机条件，运用网格逐步加密技术（upscale）进行局部加密。网格局部加密可在区域模型中嵌套局部精细刻画的地下水模型美国密西根州立大学土木工程和环境工程系开发：h

15、ttp:/www.egr. / igw/，可下载可执行文件、使用手册等。PMWIN基于MODFLOW开发的地下水流和溶质运移模型，软件中整合了IBS模块，可进行含水介质的弹性和塑性变形量计算。简单易学，可进行抽水引起的地面沉降计算。由台湾人Wen-Hsing Chiang开发：http:/ GeoTrans公司开发：http:/www.geotran- (3D Heat&Solute Transport Model，DOS)基于有限差分法的三维热及溶液运移模型。可以模拟三维空间地下水流及有关的热、溶液运移，进行地质废物处置、填埋物浸出、盐水入侵、淡水回灌与开采、放射性废物

16、处理、水中地热系统和能量储藏等问题的分析。目前国内常用于浅层地下水系统热运移模拟。美国地质调查局开发：/projects/GW_Solute/hst/，可下载可执行程序和使用说明。TOUGH2（DOS）通用的中渗流、多组分溶质运移和热运移数值模拟程序，适用二维和三维、饱和及非饱和多孔介质和裂隙介质问题，主要用于地热工程、核废料处理、环境评价与治理等方面。基于Fortran77开发的，适于在任何平台上运行，目前国内多用于地下热水运移模拟。由美国加州大学伯克利分校劳伦斯实验室开发，/TOUGH2/index.

17、htmlHYDRUS-2D/3D（windows）二维、三维饱和和非饱和带水分、溶质和能量运移模型，可进行剖面二维模拟，用户界面优化，输出的可视化表达较好。特别适合有植被条件下的土壤水-地下水流和溶质云以模拟。美国加州大学河滨分校、盐湖研究所国际地下水模型中心和捷克PC-Progress公司联合开发：http:/ www.pc- International Ltd开发：http:/www.geo-slope. com/products/，可从以上网址了解有关信息、下载学生版。MIKE SHE（windows）综合的分布式流域模型，可完整地模拟陆地水循环中地表、非饱和带和地下水的水量和溶质运移过

18、程。与AcrGIS和MODFLOW有良好的接口，界面友好，操作简单，可视化功能强。可应用于水环境影响评价、地下水地表水联合使用、地下水和地表水管理等方面。特别适合地下水-地表水联系密切条件下水资源系统的模拟。丹麦DHI国际咨询和研究机构开发：http:/ Visual Modflow介绍加拿大Waterloo大学水文地质公司开发研制，http:/www.Waterloo Modflow演示3.4 GMS简介IntroductionGMS Groundwater Modeling System is the most sophisticated and comprehensive groundw

19、ater modeling software available today. GMS is a comprehensive package that provides tools for every phase of a groundwater simulation including site characterization, model development, post-processing, calibration, and visualization.Developed by Brigham Young university4功能: 多个软件的集成MODFLOW - 2D, 3D

20、 Groundwater flow modelMODPATH - a particle tracking model based on MODFLOWMT3DMS - 3D Mass Transport modelRT3D - alternate chemical reaction based on MT3DMSSEAM3D - biodegradation model based on MT3DMSFEMWATER 3D FE Flow and MT model. Un & saturatedNUFT - 3D multi-phase flow and MT model.UTCHEM

21、 - multi-phase flow and MT model. (DNAPL) SEEP2D - 2D finite element groundwater model. PEST - Parameter ESTimation utilityUCODE - parameter estimation utility4Development- GMS简介4Development- GMS简介The backward tracking option in MODPATH can be used to delineate capture zones for wells for a given ti

22、me period4Development- GMS简介MT3DMS is a modular three-dimensional transport model for the simulation of advection, dispersion, and chemical reactions of dissolved constituents in groundwater systems. MT3DMS uses a modular structure similar to the structure utilized by MODFLOW. MT3DMS is used in conj

23、unction with MODFLOW in a two-step flow and transport simulation. Heads and cell by cell flux terms are computed by MODFLOW during the flow simulation and are written to a specially formatted output file. This file is then read by MT3DMS and utilized as the flow field for the transport portion of th

24、e simulation.4Development- GMS简介4Development- GMS简介RT3D is a modified version of MT3DMS that utilizes alternate chemical reaction packages. Numerous pre-defined reactions are available including instantaneous aerobic degradation, BTEX degradation with multiple electron acceptors, sequential anaerobi

25、c degradation of PCE/TCE, and combined aerobic/anaerobic degradation of PCE/TCE. For special cases, an option is provided for creating user-defined reactions. 4Development- GMS简介SEAM3D is a reactive transport model used to simulate complex biodegradation problems involving multiple substrates and mu

26、ltiple electron acceptors. It is based on the MT3DMS code. In addition to the regular MT3DMS packages, SEAM3D includes a Biodegradation package and NAPL Dissolution package.4Development- GMS简介4Development- GMS简介FEMWATER is a three-dimensional finite element ground water model. It can be used to simu

27、late flow and transport in both the saturated and unsaturated zone. Furthermore, the flow and transport simulation can be coupled to simulate density dependent problems such as salinity intrusion.4Development- GMS简介4Development- GMS简介NUFT is a 3D multi-phase non-isothermal flow and transport model.

28、It is ideally suited for vadose zone problems. It can be used to simulate a wide variety of remediation strategies including stream injection, vapor extraction, and air sparging.4Development- GMS简介UTCHEM is a multi-phase flow and transport model. UTCHEM is ideally suited for the simulation of surfac

29、tant-enhanced aquifer remediation (SEAR), and DNAPL transport. UTCHEM is a mature model that has been widely used. The UTCHEM solver is efficient and robust. Future versions of the GMS/UTCHEM interface will include the simulation of biodegradation. 动画演示4Development- GMS简介SEEP2D is a 2D finite elemen

30、t groundwater model. SEEP2D is designed to be used on profile models (XZ models) such as cross-sections of earthen dams or levees. The output from SEEP2D can be used to plot complete flow nets.4Development- GMS简介4Development- GMS简介4Development- GMS简介4Development- GMS简介4Development- GMS简介优点适用面广，几乎可用于

31、建立所有类型的地下水模型前处理程序好：源汇项、各类边界均有单独模块处理有大量的参考手册和源程序代码，便于根据自己的需要修改技术支持好，软件更新快对硬件要求不高价格合理缺点功能强大，但操作过于复杂。计算结果的输出、可视化表达需进一步完善。无法处理一些特殊条件，如：蒸发、排水沟、多和方向的各向异性、.4Development- GMS简介3.5 地下水模型的发展趋势 地下水模拟技术展望逐级加密技术（upscale）高精度并行算法最优参数估计（不确定性理论和方法、人工智能算法）替代模型理论与方法（孙讷正等）饱和带非饱和带耦合模型地下水地表水耦合模型流域模型（大气水、地表水、包气带水、地下水）多项流模

32、拟（核废料、CO2储存研究）地下水模拟与水资源管理系统整合研究地下水模拟软件发展趋势界面化(交互式图形界面)图形界面、下拉菜单、对话框通用化一维、二维、三维水流、水质运移、热运移、地面沉降饱和带、非饱和带地下水-地表水耦合模型前后处理的智能化自动剖分离散点插值水文地质结构识别错误的自动检测均衡计算图形、表格表示前后处理的智能化自动剖分离散点插值水文地质结构识别错误的自动检测均衡计算图形、表格表示基于GIS技术（概念模型）边界、参数、源汇项、几何参数12345678IIIIII面源（汇）分区点井分布钻孔分布（各层标高、水位标高）参数分区底图边界及地表水体Conceptual models overlay with meshes to transform the properties of conceptual models3.6 一点体会和经验地下水流数值模型地下水溶质运移数值模型地下水热运移数值模型53.6.1 地下水流数值模型要求的资料（水文地质概念模型）水文地质结构；边界条件和边界值；水文地质参数；地下水流场，地下水长观资料；地下水的补排项；准三维流（多层）与三维流区域模型与局部地段（水源地模型）的差异及关