介绍以及PMWIN的使

1、地下水数值模拟技术的重要性

地下水数值模拟技术是当代水文地质学研究和应用的最重要手段之一。这是因为地下水数值模拟技术具有特殊的优越性：1.1地下水数值模拟技术可以通过系统地分析和整理已有的水文地质、物探、钻探、调查、监测、同位素、遥感等资料信息，建立地下水流动和污染物运移的数值模型，从而在总体上刻画地下水及污染物在空间与时间上的分布与运移规律；1.2地下水数值模拟技术具有超前预报预测地下水水位变化、污染物运移、包气带水分变化等功能，从而可为工程方案设计（如污染治理方案、供水方案等）提供科学的依据；1.3

地下水数值模拟技术与计算机可视化技术相结合，可以及时的以动画与数字化生动的表达地下水运动、地下水水位变化、污染物扩散范围等，得到对地下水资源及环境污染的实时管理与监控，从而进行风险评价分析与调控2常见地下水模拟软件简介2.1MODFLOW、综合已有的MODFLOW、MODPATH、MT3D、RT3D和WinPEST等地下水模型而开发,可进行三维水流模拟、溶质运移模拟，但无法进行热量的运移模拟。MODFLOW可以模拟水井、河流、溪流、排泄、蒸散和补给对非均质地质条件和复杂边界条件的水流系统的影响，Modflow的可视化版本:visualmodflow和processingmodflow2.2Feflow由德国Wasy水资源规划系统研究所研制开发的地下水模型软件包。它采用迦辽金有限单元法进行复杂二维和三维稳定、非稳定流和污染物运移模拟，能够解决下列地下水模拟问题：完全瞬时、半瞬时、稳态地下水流动与溶质运移；随时间变化的实体属性和约束边界条件；饱和与不饱和流动；包含潜水面的承压与不承压含水层；带有非线性吸附作用、衰变、对流、弥散的化学质量运移；考虑贮存、对流、热散失、热运移的流体和固体热量运移；密度变化的流动(海水入侵等)。2.3GMS、GMS是地下水模拟系统(Groundwatermodelingsystem)的简称,是美国一所大学综合已有的MODFLOW、MODPATH、MT3D、PEST、等地下水模型而开发的可视化三维地下水模拟软件包。可使用有限差分或有限元模型对地下水流进行模拟、可进行溶质运移模拟。建立三维地层实体,任意切割剖面,产生逼真的图像、进行钻孔数据管理、二维(三维)地质统计。但也无法模拟热量运移。Modflow是GMS的核心模块，因此GMS软件也存在着Modflow的缺陷，对混合井流以温度流的模拟有待完善HST3D:一个三维热及溶液运移模型(3DHeat&SoluteTrasportModel)。可以模拟三维空间地下水流及有关的热、溶液运移,进行地质废物处置、填埋物浸出、盐水入侵、淡水回灌与开采、放射性废物处理、水中地热系统和能量储藏等问题的分析。但HST3D迄今为止还是DOS版本，三维可视化方面有待提高。两大地下水模拟软件比较Modflow以其求解方法的简单实用、适应范围的广泛及可视化功能的强大正成为最有影响的地下水模拟平台环境,其使用范围还将会不断扩大。然而，ModFlow由于采用矩形网格进行剖分，因而对于处理复杂地质体中的地下水三维渗流场模拟方面存在着不足，没有有限元三角剖分灵活多变。实践也证明,它往往并不适合非饱和流动、密度变化的流动(海水入侵)、热对流等问题。

Feflow使用有限单元法，用户可以快速构建模型来精确地进行复杂三维地质体的地下水流及运移分析，在这方面其功能要强于ModFlow。并且由于Feflow是为满足专门从事复杂地下水模拟工程的专家对技术的要求而设计的,因此,在这些VisualModflow及其它主流软件陷于困境的问题中,它更可以发挥其独特的作用。Feflow为德国软件，与平常使用的美国软件有一定差异，上手可能稍慢些。

MODFLOW的应用Modflow在全世界地下水模拟应用中占了40%以上,仅visualModflow的全世界用户就达到了1万3千多家.Feflow主要在加拿大和欧洲的一些国家应用1MODFLOW的发展MODFLOW:Themodularfinite–differencegroundwaterflowmodel

由美国地质调查局(ＵＳＧＳ)开发的用来模拟地下水流动和污染物迁移等特性的计算机程序，Modflow使用有限差分方法。其主要局限是仅在DOS模式下运行。

主要的版本：83—88—96—2000---20052MODFLOW的特点程序结构的模块化离散方法的简单化引进了应力期概念求解方法的多样化

2.1程序结构的模块化MODFLOW采用模块化结构（ModularStructure），即由一个主程序（MainProgram）和若干个相对独立的子程序包（Package）组成，各个子程序包又由不同的模块（Module）组成，供主程序随时调用。例如河流子程序包用来模拟河流与含水层之间水力联系；井流子程序包用来模拟抽水井和注水井对含水层的影响。MDFLOW的这种模块化结构使得其程序易于理解、操作、修改和添加MODFLOW问世以来，不断有新的子程序包被开发出来，例如:

用来模拟抽水引起地面沉降的子程序包(Leake和Prudic,1998)，用来模拟水平流动障碍(Horizontalflow-barrier)的子程序包(Hsieh和Freckleton,1993)等。新子程序的加入，使MODFLOW的应用范围不断扩大。优点2.1.1主程序

对于任意单元(I,j,k),MODFLOW的主程序还可以写成方程2:根据达西定律可以写出向邻单元间的渗流量:同理,可以写出其它5个面的渗流量:方便大家应用,引入水力传导系数:公式3-8可以写成10-15:将公式10-15代入上式:进行时间离散,向后差分得:公式23可以改写成:到此,对任意单元(I,j,k)的计算方程推导完毕,方程24中包含?个未知数.对于建立的整个模型,总共n个单元,可以写出?个类似的方程,共含有?个未知数.那么我们就可以利用线性代数的知识进行求解.2.1.2子程序MODFLOW采用有限差分法对地下水流进行数值模拟。差分法易于程序的普及和数据文件的规范。其主要缺点是当对某些单元网格加密时，会增加许多额外不必要的计算单元，延长程序的运行时间，随着计算机速度的迅速提高，计算机受网格数量的限制越来越小，差分法的优势越来越大，MODFLOW解决地下水流运动问题已经将含水层剖分到多达500×500×80个网格单元。2.2离散方法的简单化2.2.1三维空间的离散2.2.1单元中心法2.3引进了应力期(STRESSPERIOD)概念它将整个模拟时间分为若干个应力期，每个应力期又可再分为若干个时间段(TimeStep)。在同一应力期，各时间段既可以按等步长，也可以按一个规定的几何序列逐渐增长。而在每个应力期内，所有的外部源汇项的强度应保持不变这样就简化、规范了数据文件的输入，而且使得物理概念更为明确。最多可分成1000个应力期强隐式法(SIP)strongimplicitprocedure逐次超松弛迭代法(SOR)successiveoverrelaxation预调共轭梯度法(PCG)Pre-conditionedconjugateGradient2.4求解方法的多样化大量实际工作表明，只要恰当使用，MODFLOW也可以用来解决裂隙介质中的地下水流动问题。不仅如此，经过合理的概化，MODFLOW还可以用来解决空气在土壤中的流动（Guo,1995）；将MODFLOW与溶质运移模拟的软件结合起来，还可以用来模拟诸如海水入侵等地下水密度发生变化的问题(Guo和Benett,1997)。三維地下水模式PMWIN3D-GroundwaterModelingwithPMWINPMWIN的简介它是MODFLOW的可视化工具几乎包含了MODFLOW的所有功能和扩充模块,例如:MODFLOW2000,PMPATH,ZoneBudget,MT3D/RT3D,WinPEST

PMWIN的特点(1)面向用户的完全可视化菜单设计,符合一般操作习惯,容易上手;菜单功能区简单划分为输入、运行和输出3大部分,使用方便(2)软件可非常容易地在计算机上确定模型区域和完全自动剖分,可在原有剖分基础上任意扩大剖分范围或加密网格而不破坏原有的输入数据.(3)可视化的为边界条件赋值,智能化判别补排边界的性质,无须事先确定;允许模型设置根据地下水埋深变化调整潜水蒸发量及泉水流量.PMWIN的特点(4)有丰富的边界类型供选择.(5)对任意划定范围,能进行分区水量均衡计算,解决了水资源评价与规划中水资源量分区计算难题.(6)扩展模块能解决诸多的特殊水文地质问题.如地面沉降、湖泊萎缩、承压水顶托补给、泉水溢出减少等.适于复杂条件的应用,是我国目前开发的同类软件所不具备的.1:地下水运动模拟_____Modflow20002:粒子运动轨迹和传播时间模拟_____PMpath3:污染物质在地下水中的运动和转化模拟_____MT3D:目前国内外应用最为广泛的溶质运移数值模拟软件,zheng1990

4:水文地质参数的估计和优化_____Pest澳大利亚WatermarkComputing公司开发的PEST新版本，是功能强大的模型独立的参数估计程序，广泛应用于地下与地表水文地质学、地球物理、化学、结构与地质技术工程以及其它许多领域的模型校正和数据插值。5:区域的水平衡计算_____ZoneBudget,PMWIN的功能介绍1模拟地下水区域流场及制定地下水资源规划和管理方案;

2模拟矿区露天开采或地下开采对区域地下水的影响及其最优对策方案;

3模拟由于近海岸地下水开采或者矿区抽排地下水引起的海水或深部盐水入侵问题;

4模拟饱和带地下水流及其温度分布问题;

PMWIN的用途介绍5模拟污染物在地下水中迁移过程及其时间空间分布规律（分析和评价工业污染物及城市废物堆放对地下水资源和生态环境的影响，研究最优治理方案和对策）;

6模拟地下水过量开采引起的地面沉降等;PMWIN的边界条件处理在输入参数的时候:定水头边界:赋值为-1无流边界:赋值为0其余的边界:赋值为1层数以及边界条件栅格数目和大小各层的相对高程建立模型的空间结构输入水文地质参数应力期水平和垂向导水系数给水度和有效空隙率确定子模型运行提取结果进行分析河流、井、水库、蒸发、补给等水平衡或者三维显示PMWIN的使用过程详图PMWIN在西安广运潭项目的初步应用要求：评价有无防渗措施对周围地下水的影响，预测防渗层出现身楼以后对周围地下水的影响1）无防渗措施的模拟结果2）有防渗措施的模拟结果3）防渗措施出现问题的模拟结果实例讲解1:Step1:空间结构的划分Step2:模型参数输入对话框Step3:时间的划分Step4:子程序包的调入对话框4.1蒸发子程序包的调入4.2水库或者湖泊子程序的调入

3-314.３补给子程序的调入4.４河道子程序的调入运行结果３年以后井的影响范围１００天以后井的影响范围３年以后污染物的影响范围监测井中的浓度＿时间关系曲线实例演示2:空间划分1年后的地下水位分布填充以后的显示空间流场分布剖面的流场分布确定出不同时间长度以后污染的范围:1年3年