论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 27513阅读
  • 238回复

[实例分析]Flac3d 实例分析教程(刘 波 韩彦辉(美国)编著) [复制链接]

上一主题 下一主题
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 24楼 发表于: 2009-02-05
        
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 25楼 发表于: 2009-02-05
看看,谢谢先!!!
离线wenhanwan

发帖
80
土币
2671
威望
40
原创币
0
只看该作者 26楼 发表于: 2009-02-07
好东西一定要顶上去。
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 27楼 发表于: 2009-02-10
        
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 28楼 发表于: 2009-02-10
非常感谢楼主
离线gengdj123

发帖
889
土币
4326
威望
3161
原创币
0
只看该作者 29楼 发表于: 2009-02-11
第一章    快速入门
本章主要内容为:
¨    安装Visual MODFLOW 的硬件要求
¨    Visual MODFLOW 的安装
¨    Visual MODFLOW的基本运行过程
1.1硬件要求
为了运行Visual MODFLOW,你至少需要如下的系统配置:
l    CPU:486DX或奔腾(推荐奔腾)
l    RAM:32MB(推荐64 MB)
l    软驱:用于软件安装,高密(1.44MB),3.5’’
l    硬盘:至少有40 MB的可用空间
l    鼠标:Microsoft或其它兼容鼠标
l    操作系统:Windows 95/Windows NT 4.0(Service Pack 3)操作系统
如果在你的系统配置下软件出现了问题,请先确信你的软件是否严格按照安装说明来安装的(见安装Visual MODFLOW)。如果问题仍没有解决,请参阅设置或疑难解答,同时与你的硬件商联系。最后请阅读本章最后部分了解一下如何与Waterloo Hydrogeologic公司进行联系。
1.2安装Visual MODFLOW
Visual MODFLOW必须安装在你的硬盘或网络盘上运行。请阅读本章开始有关硬件要求部分,在安装以前确保你的系统满足这一要求。用以下可执行文件安装Visual MODFLOW:
Setup32.exe  在Windows 95/NT下安装Visual MODFLOW for Windows
这里讲述的安装过程假定Visual MODFLOW是从A盘(源盘)安装到C盘(目标盘)。
缺省安装目录为\vmodnt ,但是也可以修改它。
注意:安装盘必须去掉写保护
1.2.1 Windows 95/Windows NT 4.0下的安装过程
按如下步骤在Windows 95/Windows NT 4.0下安装Visual MODFLOW :
[1]    插入1号软盘。
[2]    进入Windows。
[3]    单击开始按钮,然后单击弹出菜单的运行项并在对话框中键入A:\Setup.EXE。Visual MODFLOW将开始安装。
[4]    单击OK按钮继续。
[5]    键入安装路径或单击OK按钮选择缺省目录开始复制Visual MODFLOW文件。
[6]    安装完成以后,单击OK按钮完成安装。
现在在你的屏幕上就可以看到一个名为“WHI Software”的程序组,窗口内会有一个VmodWin2.70的图标。你现在就可以使用Visual MODFLOW了。要从你的系统中删除Visual MODFLOW,只要双击Uninstall VmodWin2.70图标即可。它会从你的计算机中删除所有的程序文件,但不会删除工程文件。
1.3启动Visual MODFLOW
安装完毕后,双击“WHI Software”程序组中的Visual MODFLOW快捷图标或依次打开开始/程序/WHI Software/VmodWin2.70,并单击它就可以启动Visual MODFLOW了。你第一次启动Visual MODFLOW时将出现一Visual MODFLOW注册窗口,输入你的公司名和软件盒上提供的系列号。以后你的公司名称将出现在你的Visual MODFLOW打印输出中,系列号将出现在帮助/关于菜单中。
单击Next按钮以识别你所输入的资料,一旦识别成功,你就不必再次输入了。
1.4菜单简介
Visual MODFLOW的集成环境由四个主要窗口组成:主窗口、输入窗口、运行窗口、输出窗口,再加上系统设置和帮助窗口。主窗口包括如下菜单选项:
File            Input        Run        Output        Setup    Help
File        选择一种文件功能或退出Visual MODFLOW。
Input    进入输入窗口修改当前的Visual MODFLOW数据文件。
Run    进入运行窗口修改工程或批处理模式的运行时间参数和运行数值模型。
Output    进入输出窗口数值模拟结果做后处理。
Setup    选择你想用的数值引擎。
Help    关于Visual MODFLOW的主要帮助窗口和一般信息。
你的Microsoft 兼容鼠标键功能如下:
左键:    这是一个常规“单击”键。在输入框上部按下它并拖动鼠标,字符将变为高亮显示,可为后面输入的字符覆盖。
右键:    该键根据情形的不同有不同的功能。例如,在设置属性和边界条件时,它用来闭合多边形或结束一条线。在网格设置时,它用来精确地定位网格的行和列。
当要输入数据和文本时,数字和字母键才有效;所有其它键在Visual MODFLOW中都是无效的。
启动Visual MODFLOW以后,你将进入主窗口,要继续进行你的工作,你应该单击File菜单,它有六个子菜单:
1、    单击Open,打开一个已经存在的文件;或
2、    单击New,键入一新文件名创建一个新的数据文件;
3、    单击Save as,以不同的文件名保存当前模型;
4、    单击Import,装入一个已经存在的MODFLOW模型或通过下拉菜单底部的文件路径来选择一个已经存在的文件;
5、    单击Export,以AutoCad (dxf)或增强的Windows元文件格式(emf)输出文件。
你可以双击文件名或单击文件名然后按下回车键来选择一个文件。
为了获得更多的信息请查阅第三章与创立一个新的数据文件,修改一个已经存在的数据文件和装入一个MODFLOW模型有关的部分,或按照第二章一步一步地进行。
1.5窗口概览
打开一个文件并选择输入、运行或输出,将出现类似于下图的窗口。
上部菜单按钮:    提供选项,这将根据情况不同而不同
边部选项按钮:    包括视图选项和特别针对当前窗口或模型的功能选项。视图选项如下:
View Column    按列方向上显示一剖面图
View Row    按行方向上显示一剖面图
View layer    从剖面图切换到平面图
Goto    显示某一特定的行、列或层
Previous    显示前一行、列或层  
Next    显示后一行、列或层
立方体导航块:    用十字丝表示的模拟区一个简单的3D示图,用于确定空间位置。
坐标区:    显示当前光标所在的实际世界坐标,并显示当前光标所在的单元。
功能按钮:    输入、运行和输出窗口的通用功能。
F1    特定窗口或当前区的帮助。
F2    为未来的3-D选项保留。
F3    以原文件名保存文件。
F4    弹出一个文件选择窗口,你可以选择输入一.DXF格式的图形文件,该图形可以在图层窗口上打开或者关闭。
F5    按鼠标左键让你满屏放大一窗口。按右键将弹出一个窗口允许你满屏放大一特定区。
F6    使屏幕重新显示整个模拟范围。
F7    在模型上单击鼠标左键,拖动该点到新的位置后再单击左键。
F8    在按行列显示时指定垂向放大倍数。
F9    图层窗口允许你打开或关闭各种图形。
F10    返回主窗口。
状态栏:    每个按钮的功能和使用描述
CAD(计算机辅助设计)环境是Visual MODFLOW的最大特征之一。它让用户轻松快速地设置复杂问题,同时允许用户监视输入,因此避免了在数据输入过程中经常产生的许多错误。
CAD环境必须要有一个Microsoft(或兼容)鼠标。在所有的CAD环境中以下的键都是有效的:
功能键    如上所述。
回车键    选择当前高亮度显示区。
Esc键    取消当前窗口。
光标键    有限的功能(如选择菜单项)。
Tab键    单元之间移动。
1.5.1图层功能
Visual MODFLOW的显示是以一系列层的方式来组织的。各种特怔,如网格线、观测井和导水系数分区的描画是独立的,并且可以独立地打开或关闭,使你的图件或屏幕具有更好的可读性。这种不同的层称为图层(Overlay),它可以让你对数据进行交叉检验。
带有你的计算结果的图层只能在输出(Output)屏幕中使用。除此之外,输入(Input)和输出(Output)的图层是一致的。
当你在主菜单、输入、输出和运行窗口之间切换时,这些图层都有一些特殊的缺省设置。在主菜单窗口中仅显示模型边界的底图.DXF和.BMP图层。当你首次进入输入窗口时自动加上网格(Grid)图层,但当以输入模式进行操作时,缺省的图层就和前后操作有关了。在输出窗口中,缺省的图层也和前后操作有关,但是当你操作时以前的图层依然是可见的,除非你明确地将其关闭。
按F9键或用鼠标单击屏幕底部的F9按钮,将显示图层对话框。
各种各样的图层都可以让你对数据进行交叉检验,或者从你的图件中增删一些特性来使用你的图件更具可读性。输入的图件也可以从图层窗口中加以删除。
你可以用设置(Settings)按钮更改某些图层的显示特怔。分区类的特怔如导水系数和补给,既可用填充区块也可用外廓线的形式显示。DXF图件的上色以及计算所得的特怔如等值线和流速矢量的选项,都可通过设置按钮来设置。
图层在表列里出现的顺序将决定图层在屏幕上的描绘顺序。描绘顺序与表列里的顺序相反(也就是说表列中最前的特怔是屏幕上最后才描绘的特怔)。将Default按钮切换为User Defined后,图层表列可以用图层控制对话框右下方的四个按钮来作调整。
1.6 Visual MODFLOW中的打印功能
Visual MODFLOW可以让你打印当前模型的资料。从主窗口、输入、运行或输出窗口中选择File,再选择Print。你会看到如下对话框:
Project Information(项目信息)指图件标题栏的文字信息。当Full Image(全图像)被选中时,就不会打印标题栏。
对一个连贯的图形窗口,在选择放大(Zoom In)后可用鼠标右键为你的模型指定视图的区域。该特怔取得顶底的网格坐标后再对左右边界的网格坐标进行计算,形成一个矩形视域。它还可以让用户转回到上一次缩放窗口或改变当前的缩放窗口。
程序将会自动地选择缺省的Windows打印机为当前的打印设备。要改变当前的打印机,单击Select,你会看到如下对话框:
通过选择打印机名称(Name)输入框旁的箭头,你就可以将当前打印机切换到任意一台目前已安装的Windows打印机上。选择Properties会有一个标准的Windows打印选项对话框,可让你调整纸张、字体和设备的设置。
1.6.1打印到DXF文件
如果你要将结果输出到DXF文件中以便于在诸如AutoCAD之类的CAD程序中进行后处理,你可以将你的结果“打印”到DXF文件中去。
要将结果输出到DXF文件中,先从主菜单上选择File再选择Print to DXF file。下面的对话框将提示你输入DXF的文件名。
1.7设置
我们用Setup选项为不同的Visual MODFLOW软件包选择不同的数值引擎。要进入设置,点击顶部菜单的Setup,就会出现一个下拉式菜单。
1.7.1数值引擎
Visual MODFLOW允许用户为MODFLOW和MT3D选择不同的数值引擎(FORTRAN可执行程序)。选择Setup后,从下拉式菜单中选择Numeric Engines,就会出现数值引擎(Numeric Engines)对话框。
MODFLOW的数值引擎只有一个即USGS MODFLOW 96 from WHI。
MODPATH和分区均衡都只有一个数值引擎:分别是USGS MODPATH from Waterloo Hydrogeologic和 Zone Budget from Waterloo Hydrogeologic。
运移的数值引擎有四个可供选用:
l    EPA MT3D,v.1.1(公共版权);
l    DOD MT3D,v.1.1(公共版权);
l    MT3D, v.1.86 (Papadopulos & Associate Inc.);
l    MT3D96  (Papadopulos & Associate Inc.)(只有在单独购买了MT3D96后才可用)。
如果MODFLOW运行后MT3D的数值引擎发生了变化,就要重新运行MODFLOW,更新projectname.flo 文件,使之与MT3D的数值引擎保持一致。
1.7.2 Windows95/NT 交换文件
Windows95/NT管理系统的虚拟内存。要改变这些设置请询问你的系统管理员。
1.8疑难解答
1.8.1常见问题
问:为何Visual MODFLOW有时反应迟缓?
答:在数据量很大时,如果Visual MODFLOW将RAM用完,它就会占用你的硬盘作为临时空间。对许多操作而言,将不用的数据移入移出临时空间增加了运行时间。Visual MODFLOW拥有的RAM越多,对临时空间的需求就会减少。
另外,你可能需要为计算机购买一些RAM了。
问:OUT OF MEMORY- CANNOT CONTINUE 这条错误信息是什么意思?
答: Visual MODFLOW已经用掉了运行你的数据所需的内存。增加这个值就可以让Visual MODFLOW占用一些硬盘作为不常用的数据的临时存储空间。请注意,在运行前,你的硬盘上至少有这么多KB的自由空间。
问:为何在和模拟区比较时我的DXF图是HUGE(巨型)?
答:Visual MODFLOW要求你的DXF图形单位是十进制。你要输入的DXF文件可能是用英尺和英寸这种工程或建筑尺度。为了准备Visual MODFLOW 的DXF文件,你应该先用AUTOCAD或其他恰当的CAD软件编辑一下你的DXF文件。首先,你应该把单位转变成DECIMAL形式。这将把你的图形转变成十进制英寸形式。然而,这也将把你的图形文件放大成12倍。为了弥补这一点,你需要用SCALE命令把整个图以0.0833333(1/12)系数缩小,缩放完成之后,你需要用DXFOUT或用类似的命令把整个图形再保存一次。文件现在就可以以正确的尺寸输入到Visual MODFLOW中去,Visual MODFLOW的未来版本将可以识别这种情况,并自动的缩放图形。
问:我想重编译USGS MODFLOW。我的Visual MODFLOW还能用吗?
答:Visual MODFLOW使用Watcom F77编译器来编译MODFLOW源文件。另外,MODFLOW源文件还可以用Lahey F77编译器编译,它可以和Visual MODFLOW兼容。应该知道,用Fortran 编译器编译的非格式化二进制文件是不必兼容的。如果使用其它编译器,在对.HDS和.DDN进行后处理时可能会出现问题。如果你想使用其它编译器或在快速工作站运行MODFLOW,请向我们询问有关如何读取你特有的二进制文件的Visual MODFLOW用户定制的修改问题。
问:Visual MODFLOW不能创立必要的XXXX Kbytes的交换文件,为什么?
答:MODFLOW的数据文件可能非常大,经常需要比你的PC机所能提供的RAM存储空间还要大。然而,在MODFLOW处理你的问题时, Visual MODFLOW通过利用你的硬盘空间来暂时存储不需要的数据的方法来产生更多的“内存”。但是,根据你的程序大小,当你使用虚拟内存时,如果不能得到所需的总内存,MODFLOW可能仍然不能处理你的数据文件。
交换文件的大小等于虚拟内存减去可用内存,如果它大于可用的硬盘空间,那么Visual MODFLOW就会给出这条消息。
在这一点上你有三种选择:
1.    尽力按正常去运行模型。有时,某个问题可能并不需要虚拟内存参数规定的内存的大小。
2.    按指定空出硬盘空间。
3.    在 vmod.vmc文件中把虚拟内存参数改小。然而请注意,如果你的问题所需内存大于规定的值时将产生一个错误消息:DOS4Gpro Error 1307¾not enough memory when attemping to run MODFLOW。
1.9如何获得技术支持
如果你有任何问题,请随时和Waterloo Hydrogeologic公司联系。我们对已注册使用Visual MODFLOW的用户提供免费技术服务。在简介那一章中你可以找到我们的地址。
一般情况下,我们会在一个业务天中回答所有技术服务问题,但是可能要两个业务天来作答复。
在信中除了简单描述你的问题外请写上你的:
l    名字
l    公司和办公室
l    电话和传真号码
l    email地址(如果你有的话)
l    产品序列号
l    版本注册日期
你的版本注册日期可以在主菜单窗口中的HELP菜单的About Visual MODFLOW子菜单中找到。
如果你的问题是“我要如何去做...”之类的话,附上这些资料就足够了。但是,如果你是报告一种故障(特别是导致系统崩溃的故障)或有关该软件的技术问题,请另外提供如下信息:
l    计算机的硬件配置。
l    操作系统(Win 95 /WinNT)。
l    计算机的RAM大小。
l    计算机的剩余硬盘空间大小。
l    详细描述问题是在什么时候和如何发生的。
如果我们不能重现你所遇到的问题或不能解决此问题,我们将要求你给我们寄一份你的工程文件的副本。

第二章    使用教程
本章是Visual MODFLOW的使用教程,它用一个具有完整的输入、输出文件的模型例子(Airport)来让你熟悉后处理能力和Visual MODFLOW的功能。该问题的数值模拟结果(MODFLOW、MODPATH和MT3D)已经完成,你可以用它来评价一下该例模拟结果的可视化输出能力。
该章将向你显示设计并运行模型以及如何将结果可视化显示的一些必要的步骤。本教程是以按一步一步的方式组织的,你可以选择感兴趣的章节而不必做完全部练习。
2.1实例模型介绍
本场地位于Waterloo城外的一个机场附近。场地表层地质情况是:上下各有一层砂砾含水层,中间有一粘土和粉砂弱透水层将它们分开。相关的场地地物有一个飞机加油场、一个生活供水井场和一块弱透水层不连续区(含水层天窗)。如下图所示。
生活供水井场有两眼井。东井抽水量稳定在550m3/d,西井为400 m3/d。在过去十年中,在加油场中不断地有飞机燃料的泄漏,其自然下渗已在上层含水层中产生了一个污染羽。本教程将向你显示如何一步步地为这个场地建立一个地下水流动和污染物运移的模型。这个模型将显示燃料污染物对生活供水井的潜在影响。
在讨论时,我们将该场地在平面上的视图指定为上北下南左西右东。如图所示,地下水在一个由上部潜水含水层、中间弱透水层和下部承压含水层所组成的三层系统中自北向南(平面图上是从上到下)流动。上下含水层的水力传导系数为2e-4m/s,弱透水层的水力传导系数为1e-10m/s。
2.2有关的术语和记号
在本教程中,将用到以下的术语和记号:
TYPE    输入给定的词或值
SELECT    在所指的地方按下鼠标左键
«    按下〈Tab〉键
¿    按下〈Enter〉键
+    在所指的地方单击鼠标左键
++    在所指的地方双击鼠标左键
2.3启动Visual MODFLOW
在Windows环境下,你只需简单地双击WHI Software程序组里的Vmod图标,即可启动Visual MODFLOW。
现在你就可以看到Visual MODFLOW的启动屏幕。要进入主菜单:
+    OK
2.4打开文件
在本教程中,你可以从任意地方开始,在章节之间跳来跳去,或者离开,或者从你停止的地方开始。要创建一个新文件,请跳至模块I。要打开一个已有的模型,在顶部菜单条上:
+    File
会有一个下拉菜单提示你从一个标准文件选项的表列中进行选择。
+    Open
将会出现一个文件选择窗口。
++    Tutorial文件夹
如果你要开始学习第2节:
+    Grid1.vmf
如果你要开始学习第3-7节:
+    Grid2.vmf
如果你要开始学习第8-10节:
+    Airport.vmf
如果你要开始学习第11-14节:
+    Output.vmf
如果你要继续本教程:
+    filename.vmf(你在第1节创建的文件)
+    OK打开文件
几秒钟后,主菜单就会出现,同时显示出模型网格和场地地图,如下图所示。
2.5模块I:模型输入
2.5.1第1节:设置数值引擎
在进入输入屏幕以前,很重要的一项工件是选择运行模型所需的数值引擎。从顶部菜单条上:
+     Setup
+     Numeric Engines
新模型的缺省数值引擎如下所示:
重要说明:
Visual MODFLOW包括免费的MT3D的公共版本(EPA的1.1版和DOD的1.5版)。但是1.10版没有1.86版稳定,它们都没有MT3D96那么快而灵活。本教程的飞机场模型练习是用MT3D v.1.10和MT3D96做的。对那些需要精确模拟污染物运移的专业咨询公司来说,所有的模拟都应该用MT3D96来做。
要换一种数值引擎的话,单击你想要的改变的那种数值引擎的选择框旁边的向下箭头键,然后在下拉对话框中选择。对于本教程中的模拟练习,你可以使用DOD1.5或MT3D96,如果你已分别购买了的话。
+     Transport Numeric Engineer
+     MT3D by Papadoplulos_Associates Inc.
+     OK
2.5.2第2节:生成模型网格
第一个模块将教你如何利用Visual MODFLOW模拟环境通过必要的步骤生成一个新的模型数据文件。
+     File
+     New
此时将出现一个创建新文件的窗口。
++    在Tutorial文件夹中
在文件名输入框中输入一个新的数据文件名(比如,键入‘Tutorial’),然后选择Save,就生成了一个新的数据文件。
Visual MODFLOW会自动在文件名后加入vmf扩展名。
然后,将出现一个一个选择单位的窗口。在下图中用鼠标为你的特定系统指定度量单位。
+    3Create model using base Map
下一步你必须说明.DXF底图文件的目录和文件名。
+    Browse
在你安装Visual MODFLOW的目录(缺省为:c:\VMODNT\tutorial)下找到并选中下面这个文件:
+    sitemap.dxf
+     Open
下一步你应该输入模型的行数、列数以及层数。
Enter Number of Columns (j)(列数)    40Ö
Enter Number of Rows (i)(行数)    40Ö
Enter Number of Layers (k)(层数)    3Ö
Enter Zmin(最小的高程值)    0Ö
Enter Zmax(最大的高程值)    18Ö
在上述窗口的Units(单位)部分,用鼠标从它们各自的下拉列表选择如下选项:
Length(长度)    meters
Time(时间)    days
Conductivity(传导系数)    m/second
Mass(质量)    kg
Concentration(浓度)    mg/L
Pumping Rate(抽水速率)    m3/day
Recharge(补给速率)    mm/year
+    Create          (确认这些设置)
之后,将出现一个模型区域选择(Select Model Region)窗口, 你可以在这里定义模拟区的范围。Visual MODFLOW会从sitemap.dxf(场地底图)中读取坐标的最大值和最小值,并默认它是处于模拟区域的中心。现在的Visual MODFLOW允许用户转动dxf图像、使用局部坐标系、设置DXF图的范围。如果用位图作为底图,则可以对图形进行地理参照处理和旋转。在屏幕中的原来的数字上输入以下数值:
Enter Start Point X    0Ö
Enter Start Point Y    0Ö
Enter Angle    0Ö
Enter Grid Co-ordinates X1    0Ö
Enter Grid Co-ordinates Y1    0 Ö
Enter Grid Co-ordinates X2    2000 Ö
Enter Grid Co-ordinates Y2    2000Ö
+    OK 接受网格尺寸
在模拟区内将会自动生成一个大小均一的4054053的有限差分网格,同时在屏幕上将会出现一个如下图所示的场地底图。
从这儿你可以进入输入部分的任何一节,看看Visual MODFLOW是如何处理各种输入参数组的。
前面已经提到,如果是想看看结果的显示,是没有必要将本教程所写的所有步骤全部运行一遍的。本教程已提供了每个模块所需的所有文件。
2.5.3第3节:细化模型网格
本节介绍了对感兴趣的地段如供水井周围地区、含水层天窗区与加油场之间地段进行模型网格细化的步骤。由于单元的大小代表了井的尺寸,因而用更细的网格来模拟井的影响将会更接近实际情况。另外,如果在井周围出现了水位下降,用小的网格剖分可使水面更平滑。
如果你是从别的小节转到此地,你可以跳过框中的步骤。
+   File
+   Open
此时针对出现一个文件选择窗口。
++   Grid_1.Vmf  (Tutorial目录下)
打开文件后,就会出现Visual MOFLOW的主菜单。
首先我们要细化供水井周围的网格。从主菜单的顶部菜单条上:
+    Input
就进入了输入模块中,其中模型的组成部分(Building Blocks)在顶部以菜单方式给出。当你首次进入输入模块时,Visual MODFLOW会装入网格(Grid)输入窗口。屏幕左边的最上面的四个按钮(View Column、View Row、View Layer、Goto、Next和Previous)在每个屏幕上都会出现,使你在任何时候都可以将模型的显示方式从平面变为剖面。屏幕左边其它的按钮有各种进行模型网格修改的功能。要细化X方向的网格,
+    Edit Columns
然后将鼠标移至网格的任一处。注意在模型网格上有一个高亮度的竖线会随鼠标移动。这条竖线可以用来在模拟区的任何位置加上一列。在本练习中,我们通过按鼠标右键引出一个添加垂直网线(Add Vertical Line)窗口,以指定的间隔增加网格线。如图所示。
+    ¤Every spaced gridlines from:
单击其后的输入框,输入以下数值:
from:    500Ö
to:    1600Ö
at intervals of:    25Ö
+    OK接受这些值
+    Close离开
现在我们要细化从加油场至供水井Y方向的网格。
+    Edit Rows
然后将鼠标移至网格的任一处,按鼠标右键,将会有一个弹出窗口要求你输入网格的信息。
+    ¤Every spaced gridlines from:
单击其后的输入框,输入以下数值:
from:    400Ö
to:    1900Ö
at intervals of:    25Ö
+    OK接受这些值
+    Close离开
细化后的网格如图所示。
下一步将向你展示如何在剖面上观察模型,以及向模型添加层位。要在剖面上观察模型,在左边菜单上
+    View Column
将鼠标移至网格的任一处,当你在屏幕上移动光标时,光标对应位置的列就会以高亮度显示。如要观察某一列,在该列上单击左键。Visual MODFLOW将模型的显示方式从平面变为剖面视图。此时模型还没有进行垂向放大,剖面图看起来就象一条粗线,根本不能区分出三个层来。为了能恰当地显示出这三个层,你必须对剖面图进行垂向放大。在窗口底部选择:
+     F8-Vertical Exaggeration
就会出现一个垂向放大窗口,供你输入垂向放大系数。
Type:    40
+     OK
如下图所示,这模型的三个层就会在屏幕上显示出来。从图上可以看出在整个剖面上每一层的厚度都是一样的。
在Visual MODFLOW中,你可以从Surfer的.grd文件或以空格分隔的x, y, z的 ASCII文件中输入变化的层位标高。本例中,我们将输入一个ASCII格式的x, y, z文件来生成一个有坡度的地表及层位,同时后者的厚度也是变化的。
+    Import Surface
此时将出现如下窗口:
+     Choose Filename
显示一个所有ASCII文件的选择窗口。你将输入的地表高程文件名为airpt_gs.asc。要选择此地表高程文件,
++    airpt_gs.asc
应确认以下缺省设置。
+    OK
这是一个有一定坡度的地表,北部边界高程为20米,南部边界高程为17米。Visual MODFLOW从ASCII文件中读取X,Y,Z数据,并利用五个最近采样数据点利用距离平方倒数插值方法求出第1层的每个单元顶板的块中心高程。
现在要输入层1底板的高程。
+    Import Surface
+    Choose Filename
++    airpt_b1.asc  第1层的底板高程数据文件。
+    ¤ Import bottom elevation of layer(输入第...层的底板高程):(默认为1)
+    Minimum Layer Thickness   0.75
+    OK    接受其它的默认参数值
几秒钟以后,第一层底板就会显现出来,它具有一个变化的表面。下一步要输入第二层底板高程。从左侧菜单上选择:
+    Import Surface
+    Choose Filename(选择输入的文件)
++    airpt_b2.asc  第2层的底板高程数据文件。
+    ¤ Import bottom elevation of layer:向框中输入2。
+    Minimum Layer Thickness   0.75
+    OK    接受其它的默认参数值
第二层的底板就会显现出来,它也具有一个变化的表面。接着输入第三层底板高程。从左侧菜单上选择:
+    Import Surface
+    Choose Filename(选择输入的文件)
++    airpt_b3.asc  第3层的底板高程数据文件。
+    ¤ Import bottom elevation of layer:向框中输入3。
+    Minimum Layer Thickness   0.75
+    OK    接受其它的默认参数值
第三层的底板就会显现出来,它具有一个变化的表面。
为对模型网格进行更细的垂向离散,你可以对每一层进行细分。从左边的菜单中选择:
+    Edit Layers
把鼠标移到模型横剖面中。在横剖面中将会有一条水平方向的线随着鼠标移动。要细分某个层位,只要让这条线移到第一层大致的垂向中央并单击鼠标左键。在此处就会加入一条分层线。
为了细分中部承压层,把这条水平线移到弱透水层中单击鼠标右键。将出现一个添加水平层位(Add Horizontal Layer)窗口,如下所示。
+    ¤ Split current layer into ‘2’ evenly spaced layers(将当前层分为2个等大小的层)
+    OK
对模型的底层重复这一过程。
当你完成这些操作后:
+    Close  退出
现在,模型的横剖面上应该有六个层位,类似下图:
用左边菜单的Next和Previous按钮在模型列剖面中移动。
接着选择View Row,在剖面的一高亮垂直列中用鼠标选择一行。使光亮区在模型剖面的行上移动。
要返回到模型区的平面视图状态,选择View Layer按钮,然后使模型顶层变成高亮显示,然后单击鼠标左键。就会显示飞机场的平面图。
2.5.4第4节:添加井位
这一节的目的是一步一步地指导你怎样在模型中加入抽水井。
如果你已经阅读完了其它小节,你可以跳过该框中的步骤。
+           File
+           Open
将出现一个文件对话窗口。
++       grid_2.vmf(细化网格后的模型)
此时文件就会打开,出现Visual MODFLOW主菜单窗口。要进入输入模块:
+           Input
从顶部菜单中选择
+    Well
再从下拉菜单中选择:
+    Pumping Wells
系统将要求你保存你的数据
+    Yes
现在开始加入井位,将供水井周围区域(模拟区的右下部)放大:
+    F5-Zoom In
把光标移到供水井的左上方,单击鼠标左键。然后,拖出一个方框围住供水井,再次单击左键,闭合放大窗口。
+    Add Well
把光标移到西边那眼井上,单击鼠标左键。就会出现一个井位编辑窗口,你可在此输入井位资料。
在标有Well Name的输入框中单击鼠标:
Type:    Supply Well 1
+    Add Screen
把光标移到窗口左边的井孔中,在高程为5米处单击一下鼠标。然后把鼠标移到高程为0.3米处再单击一次。就会在井中高程大约为5-0.3米处加入一滤水段。另外你也可以在方框中直接输入数值来精确确定滤水段的位置。
因为Visual MODFLOW计算的是每个单元的水头值,而不是井的水头值,因此Visual MODFLOW允许你把井当作观测点。然后你可以输入抽水井的半径,Visual MODFLOW会去计算校正后的水头值。这在与抽水井的实测降深进行校验时非常有用。
下面我们要输入井位抽水的时间安排,在标有Stop的列下的方框内单击鼠标左键,并输入以下信息:
Stop[day]:    7300Ö
Rate[m3/d]:    -400Ö
+    OK  
如果某些所需的数据没有输入的话,Visual MODFLOW会提醒你马上完成。
下一步是设置第二眼抽水井的参数,我们要使用从一眼井向另一眼井复制属性的快捷方法来完成。
+    Copy Well
把光标移到西边的那眼井上,单击鼠标。然后把光标移到东边井位上,再次单击一次,就把井复制完成。下一步我们要编辑复制过来的井的资料。
+    Edit Well
单击井位标记,选择新的那眼井。就会出现一个井位编辑窗口。
在标有Well Name的输入框中单击鼠标:
Type:    Supply Well 2
在标有Rate[m3/d]的方框中单击鼠标:
Type:    -550
+    OK  接受这些井参数
+    F6-Zoom Out
2.5.5第5节:设置模型的属性
本节将向你显示如何一步步地设置一个由水力传导系数相差很大的层位所组成的模型的属性。
如果你已经完成了其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   Tutorial文件夹
++   grid_2.vmf
打开文件,出现Visual MODFLOW主菜单。
+    Input
+    Properties
+    Conductivity
+    Yes  (保存你的文件)
如果你正在创建一个新的模型,在此将出现一个默认K值(Default K Property) 窗口,让你输入水力传导系数、储水系数、给水度以及空隙度的初值。模型会把这些值赋给模型中的每个单元。
我们可以使用Assign Single、Assign Polygon或Assign Windows来输入具有不同的水力传导系数的区域或层位。
用给定的格式输入如下值。
Hydraulic Conductivity in x (m/sec) Kx:    2e-4Ö  
(X方向水力传导系数Kx)
Hydraulic Conductivity in y (m/sec) Ky:    2e-4(输入Kx后自动出现)Ö  
(Y方向水力传导系数Ky)
Hydraulic Conductivity in z (m/sec) Kz:    2e-4Ö  
(Z方向水力传导系数Kz)
Specific Storage (Ss):    1e-4Ö
(单位储水系数)
Specific Yield (Sy):    0.2Ö
(单位给水度)
Effective porosity (Eff.Por):    0.15Ö
(有效孔隙度)
Total Porosity (Total Porosity):    0.15
(总孔隙度)
+    OK 接受这些数值。
当你输入Kx的值时,Ky就自动给出了,因为Visual MODFLOW假设模型是水平各向同性的。
现在输入第三和第四弱透水层的水力传导系数值。
+    Goto
显示转到层弹出窗口
Type:    3
+    OK
Assign Windows让你可以在一个矩形窗口内设置一个不同的水力传导系数。现在,让我们给第三、四层弱透水层输入一个较小的水力传导系数。
+    Assign Windows
把鼠标移到网格的西北角,在单元中间单击。然后把鼠标移到东南角,点击单元中部。这就形成了一个覆盖整个层面的窗口。
出现一个属性设置(Assign Property)对话框。
+    New (整个网格将变成蓝色)
输入下列弱透水层的水力传导系数:
把鼠标移到Kx(m/s)框中单击并输入如下值:
Kx(m/s)    1e-8Ö
Ky(m/s)    自动设置为1e-8Ö
Kz(m/s)    1e-9Ö
+    OK  接受这些值
从底部菜单栏中选择F6-Zoom Out更新屏幕显示。
现在把第三层的水力传导系数复制到第四层。在左边菜单上选择:
+    Copy Layer
出现一个复制窗口。在窗口中:
+    þ Copy All Properties(复制所有的值)
+    Layer 4(它就会以高亮度显示)
+    OK 把第三层的K值复制到第四层
尽管在含水层天窗地段代表弱透水层的模型层已收缩到非常薄(0.2m),但是这对于重新设置K值来表现弱透水层的不连续性是一个很好的练习。
在这个特例中,含水层天窗区在底图上就已有体现。然而,在许多情况下是没有这种资料的,你必需依靠某些参数比如层厚来确定这个不连续的地段。要显示第三层的厚度可在下部菜单条上选择F9-Overlay按钮。这时就会出现一个图层控制(Overlay Control)窗口(如下图),其中有一个表列,上面包含以字母顺序排列的所有在设置模型输入参数可能出现的图层。
向下浏览表列直到找厚度图层(Layer Thickness Overlay),然后双击它,在它的左边就会出现一个“*”,表示它现在是一个活动图层。
+    OK 将显示层厚的等值线图
现在将等值线收缩到0.5m的最小厚度的地段放大。
+    F5-Zoom In
将鼠标移动到含水层天窗地段的左上部,单击左键。在该区上拉出一个框,再单击一次。屏幕上就会出现含水层天窗地段放大后的图形。
+    Assign Single(对单个单元进行设置)
将会出现一个属性设置框,其缺省值是最后输入的K值(Property #2)。按下New键下的向下箭头。Property #1 就会是K的当前值。此时千万别选择OK。将鼠标移到由0.5m等值线所圈定的表示含水层天窗的地段。按下鼠标,绕着该区拖动直到区内单元格成为白色显示(如图)。
如果你已对区外的单元重新作了设置,那么,松开左键,按下右键就可以将单元恢复为它们原先的值的颜色。完成单元格的“描画”后,选择Assign Single窗口中的OK键。
现在你要将该值复制到第四层。
+    Copy Layer(在左侧菜单条上)
就会出现一个复制窗口,其缺省设置是Copy only property #1(仅复制第一种属性)
+    Layer 4 激亮该层
+    OK 把层3中的K的Property #1值复制到层4上
要检查K值是否已经复制到层4,在左边菜单上选择Next进入层4。
+    F6-Zoom Out 回到整个模拟区的全屏显示状态。
要在剖面上观察模型,从左边菜单选择View Column键,移到含水层天窗地段。通过左键选择要观察的列。就会出现一个如图所示的剖面图。
用左侧菜单上的Next和Previous按钮可在整个模型中移动。通过屏幕左下侧的立方体导航器可以观察剖面的相对位置。
选择View Layer,然后点中模型剖面中的第1层,回到平面显示方式。
然后Visual MODFLOW将会计算模型顶层的层厚度等值线。
要去掉层厚等值线,可选择底部菜单的F9-Overlay。浏览图层表列,双击Layer Thickness Overlay,将它变为不活动即可。
+    OK
2.5.6第6节:设置模型边界条件
本节介绍了设置各种模型边界条件的步骤。
如果你已经完成了其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   在Tutorial文件夹中
++   grid_2.vmf
打开文件,出现Visual MODFLOW主菜单。
+    Input
要设置补给条件,你必须是在处于模型顶层的平面视图情况下。检查一下左下角的立方体导航器,看看你目前所处的层位,如果不在第1层,就用Next,Previous或Goto进入到第1层。
我们首先要设置的边界条件是含水层的补给通量。
+    Boundaries
+    Recharge
如果你是正在建立新模型,此处会出现一个缺省补给(Default Recharge)窗口,提示你输入补给量的初始值。Visual MODFLOW自动地将这个补给量值赋给模型的整个顶层。在框中输入:
Stop Time[days]:    7300Ö
Recharge[mm/yr]:    100
+    OK接受赋值
+    F5-Zoom-In
将光标移到加油场的左上角,单击鼠标左键。沿加油场拉出一个框,然后再点一次闭合放大窗口。
现在我们要给加油场设置一个较高的补给值:
+    Assign Window
将光标移到加油场的某个角上。单击鼠标左键,沿加油场拖出一个窗口,然后再单击鼠标左键。就会出现Assign Recharge窗口。
+    New
设置了补给条件的窗口将变为Property #2属性的蓝色,给Property #2输入以下值:
Stop Time [days]:    7300Ö
Recharge [mm/yr]:    250
+    OK接受赋值
+    F6-Zoom-Out将模拟区恢复为满屏显示。
下一步是给承压和非承压含水层的南北边界设置定水头边界条件。
从顶部菜单上选择:
+    Boundaries
+    Constant Head
+    Yes(保存输入的补给量数据)
首先我们为上部潜水含水层沿模拟区北部的边界设置的边界条件。
+    Assign Line
将鼠标移至网格的西北角。用鼠标左键点击该单元中心。然后用鼠标右键点击该单元中心。此时单元格组成的水平线将以高亮度显示,同时出现一个如图所示的定水头边界输入窗口。
输入以下数值:
Code #:    1
+    Stop Time
Stop Time:    7300Ö
Start Point:    19Ö
End Point:    19
+    OK接受输入
此时,粉色线将会变成暗红色,表示常水头边界值已经设置完毕。从左边菜单上选择:
+    Copy Layer
会出现一个复制窗口,在Copy only code #(仅复制第...种属性)缺省设置旁边的框中值为1。
+    Layer 2(该层会高亮度显示)
+    OK 将定水头边界1号属性复制到第2层
我们现在要输入下部承压含水层沿模拟区北部的边界的边界条件。
+    Goto(在边上的菜单条上)
此时会有一个Goto Layer的弹出窗口。在Layer you wish to go to(你想去的层位)框中,1是缺省值,并高亮度显示。
Type:    5
+    OK进入第5层
+    Assign Line
将鼠标移至网格的西北角。用鼠标左键单击该单元中心。然后用鼠标右键单击网格东北角的单元的中心。该线会以高亮度显示,同时出现一个定水头边界输入弹出窗口。输入以下数值:
Code #:    2
+    Stop Time
Stop Time:    7300Ö
Start Point:    18Ö
End Point:    18
+    OK接受输入
+    Copy Layer
会出现一个复制窗口,选择Copy only code #并在旁边框中输入“2”。
+    Layer 6(该层会高亮度显示)
+    OK 将定水头边界2号属性复制到第6层
下一步我们要设置下部承压含水层沿模拟区南部的边界的边界条件。
+    Assign Line
将鼠标移至网格的西南角。用鼠标左键单击该单元中心。然后用鼠标右键单击东南角单元的中心。该线会以高亮度显示,同时出现一个定水头边界输入弹出窗口。输入以下数值:
Code #:    3
+    Stop Time
Stop Time:    7300Ö
Start Point:    16.5Ö
End Point:    16.5
+    OK接受输入
+    Copy Layer
会出现一个复制窗口,选择Copy only code #并在旁边框中输入“3”。
+    Layer 6(该层会高亮度显示)
+    OK 将定水头边界3号属性复制到第6层
设置完定水头值后
+    View Column
将鼠标移至模拟区,选择通过含水层天窗地段的列。将会出现一个与下图相似的剖面图。
同时它还可以显示表示水力传导系数属性的颜色。选择F9-Overlay,将显示图层窗口,  此窗口显示了一系列能显示和关闭的图层。使用鼠标双击Conductivity Overlay。将在Conductivity Overlay旁边显示一星号,表明它被激活了。
+    OK 显示水力传导系数层
现在从左边菜单栏选择View Layer并单击模型横剖面的顶层,返回显示模拟区的平面图。
下一步沿着模拟区南部边界的顶层输入河流边界条件。
+    Boundaries
+    River  切换到河流输入窗口
+    Assign Line
按照图层布置,(从网格的西南部开始)沿着河流单击鼠标左键使河流数字化。当到达东南边界时单击鼠标右键。
此时会出现一个弹出窗口将让你输入河流的资料。
键入以下数据:
Code #     4
þ Assign to appropriate Layer
+    Stop Time
Stop Time    7300Ö
Start Point River Stage    16.0Ö
Start Point River Bottom    15.5Ö
Conductance    1000 Ö
End Point River Stage    15.5 Ö
End Point River Bottom    15.0Ö
Conductance    1000
+    OK  接受这些值
+    F6-Zoom Out  刷新屏幕
定义河流后,将显示一条用网格单元表示河流边界条件的蓝线,如下图所示。

注意:技术角度上来说,这些传导系数值是错误的,因为在MODFLOW 中传导系数的定义如下:
传导系数=    河床的水力传导系数    ´河流面积
    河床厚度    
细化网格时产生不同表面积的单元,相应地应输入不同的传导系数值。为了节约时间,在该例中不给单个单元赋予不同的传导系数值。
2.5.7第7节:设置质点
本节将一步一步地指导你输入向前追踪的质点以确定污染物的途径。  
如果你已经完成了其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   在Tutorial文件夹中
++   grid_2.vmf
打开文件,出现Visual MODFLOW主菜单。
+    Input
在顶部菜单栏上选择:
+    Particles
此时程序会要求你保存数据:
+    Yes
我们在此要做的第一步是放大加油场地段。
+    F5-Zoom In 在底部菜单栏
在加油场的左上角单击鼠标左键,然后在该区上拖出一窗口,再单击鼠标左键闭合放大窗口。
+    Add Line
把光标移到加油场的左边并单击鼠标左键。拖出一根线到加油场的右边,然后再单击一次。将显示一添加质点窗口。缺省的质点数是10。把质点数变成5。
+    OK  设置好加油场中的五个质点的线
通过加油站区的绿色质点线表示是向前追踪的质点。现在返回到模拟区的满屏显示状态。
+    F6-Zoom Out
2.5.8第8节:设置模型的运移参数
如果你有MT3D1.1, 1.5或1.86,请进入第8A节。
如果你购买的是MT3D96,请进入第8B节。
2.5.8.1第8A节:MT3D1.1, 1.5或1.86
本小节将一步一步地向你介绍如何输入运行MT3D所需的污染物运移参数。
如果你已经完成了其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   在Tutorial文件夹中
++   grid_2.vmf
打开文件,出现Visual MODFLOW主菜单。
+    Input
在顶部菜单栏上选择:
+    MT3D
此时程序会要求你保存数据。
+    Yes
如果你以前没有进入过MT3D窗口,程序就会要求你输入质量运移参数所用的质量和浓度的单位。
然后程序会要求你输入整个模拟区的背景浓度。
Type:    0
+    OK 设置好背景浓度
下一步我们要为模型的每一层输入一个污染物延迟系数。
+    Chem. reaction (顶部菜单栏上)
+    Options/Defaults (边部菜单栏上)
此时屏幕上会将显示一个化学反应(Chemical Reaction)窗口,如图所示。
在本例子中我们只使用吸附等温线选项。
+    þ Sorption isotherm
我们要使放射性衰变或生物降解选项(Radioactive decay or biodegradation Option)无效。如果它前面是þ,请在方框中单击鼠标去掉ü。注意现在窗口右下角的两个小窗口变灰了,表明它们已变为无效。
现在我们要选择吸附等温线的类型。
+    ¤ Linear isotherm
注意因为选择了线性吸附等温线,窗口右上角的Second Sorption Constant(吸附方程第二常数项)窗口已经无效了。
下一步,你必须输入每层的土壤容重值。
在窗口左下部的标有Bulk density of porous medium (Kg/m3) (孔隙介质容重)的窗口中的第一个方框中双击鼠标。方框左边的数对应于层号。每层的土壤容重值都有设为2650 Kg/m3。
现在我们将输入用于计算污染物延迟系数的Kd值。
用鼠标单击First Sorption Constant(吸附方程第一常数项)窗口中第一个方框。
Type:    5.7e-5Ö
所有六层都输入该值。
+    OK  接受这些化学反应参数
2.5.8.2第8B节:MT3D96
本小节将向你介绍如何设置运行MT3D96所必需的污染物运移参数。
如果你已经完成了其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   在Tutorial文件夹中
++   grid_2.vmf
打开文件,出现Visual MODFLOW主菜单。
+    Input
从底部工具条上选择:
+    F10-Main Menu
此时程序会要求你保存数据。
+    Yes
+    Setup
+    Numeric Engines
新模型的缺省数值引擎如下图所示:
重要提示:
Visual MODFLOW 包括免费的MT3D1.10(EPA)和MT3D1.5(DoD)。然而,MT3D1.10没有MT3D1.86那么稳定,且这些版本都没有MT3D96的速度快或有弹性。本教程中的机场练习是用MT3D1.10和MT3D96都做过。对那些要求进行精确的污染物运移模拟的所有专业咨询业务都应使用MT3D96。
要改变某个数值引擎,单击该引擎旁的下拉箭头键,在下拉框中选择你想要的数值引擎即可。在我们这次模拟中,我们将使用MT3D96(如果你没有购买MT3D96,你可以使用具有MMOC技术的MT3D1.10或MT3D1.5)。
+    Transport Numeric Engine
+    MT3D96 by Papadopulos _Associates Inc.
+    OK 接受数值引擎
+    Input设置
+    MT3D
然后你可以输入整个模拟区的背景(缺省)浓度(如果你事先没有进入过MT3D菜单)。
Type:    0
+    OK 接受背景浓度
下一步你应输入每层的污染物延迟系数。
+    Chem. reaction
+    Options/Defaults
将出现化学反应窗口,如图所示。
在本例中我们只使用吸附等温线选项。
+    þ Sorption isotherm
在方框中单击鼠标去掉ü,使放射性衰退或生物降解选项无效。
现在你必须选择吸附等温线的类型。
+    ¤ Linear isotherm
+    OK
MT3D96允许你按单元输入化学反应参数。缺省时,MT3D96给模型的每一层都赋予模型开始建立所定义的值。
+    Database
你必须为每一种化学反应属性输入土壤容重。
在标有Bulk density (Kg/m3) 的窗口的第一个方框中双击鼠标。方框左边的数对应于层号。把每层的土壤容重都设为2650 Kg/m3 。
现在我们将输入用于计算污染物延迟系数的Kd值。
在第一项系数(1st coeff.)窗口下的第一个方框中单击鼠标。
Type:    5.7e-5
给所有的层都赋予这些值。
+    OK  接受这些化学反应参数。
2.5.9第9节:标志污染源
加油场的污染源将以污染物下渗源的补给浓度(Recharge Concentration)的方式给出。
从顶部菜单条上选择:
+    Boundaries
+    Recharge Conc.
检查一下左下角的立方体导航器,确信你已经处于第一层上。如果不在第1层,就用Next,Previous或Goto进入到第1层。
+    F5-Zoom In
同时在加油场上拉出一个放大窗口。
然后,给进入模型的加油场的天然补给量设置一个浓度:
+    Assign Polygon
沿加油场周长移动鼠标,并同时用鼠标左键单击该区的角点,然后单击鼠标右键,该区就会变成阴影图像。此时就会出现一个补给通量浓度设置(Assign Recharge Flux Concentration)窗口,输入以下数值:
Code #:    1Ö
Stop Time[days]:    7300Ö
Recharge Conc.[mg/l]:    5000Ö
+    OK接受这些补给浓度值
+    F6-Zoom Out回到模拟区全屏显示状态
MT3D最后一个需要输入的参数是弥散资料。从上部菜单条上选择:
+    Dispersion
如果你以前没有输入过弥散资料,程序就会让你输入纵向弥散度和横向弥散度比值。
Long. Dispersion[m](纵向弥散度):    0.1Ö
Horiz. to Long. Ratio(水平横向/纵向):    0.1Ö
Vert. to Long. Ratio(垂直横向/纵向):    0.00001Ö
+    OK接受弥散资料
在这个例子中,模型中每层的弥散系数都是均匀的。
最后,我们要向模型中加入三个观测井以监测整个模拟区内不同位置的浓度。
+    Calibrate
此时程序会提示你保存数据。
+    Yes
+    Add Obs.
将鼠标移至模拟区的任意一点,单击左键。将会有一个编辑观测点(Edit Observation Point)窗口让你输入以下信息:
Observation Name:    OW1Ö
X Location:    760Ö
Y Location:    1667Ö
+    Time            Observed Conc.[mg/l]
500Ö                500¿
1000Ö                1600¿
2000Ö                2700¿
7300Ö                3300
+    OK接受这些数值
+    Add Obs.
再将鼠标移至模拟区的另外任意一点,单击左键。将会出现一个Edit Observation Point窗口让你输入以下信息:
Observation Name:    OW2Ö
X Location:    760Ö
Y Location:    1350Ö
+    Time            Observed Conc.mg/l
7300Ö                800¿
+    OK接受这些数值
+    Add Obs.
再将鼠标移至模拟区的另外任意一点,单击左键。将会出现一个Edit Observation Point窗口让你输入以下信息:
Observation Name:    OW3Ö
X Location:    760Ö
Y Location:    900Ö
+    Time            Observed Conc.mg/l
7300Ö                100
+    OK接受这些数值
要进入下一模块,请在选底部菜单条上选择F10-Main Menu,再选择Yes保存数据。
离线gengdj123

发帖
889
土币
4326
威望
3161
原创币
0
只看该作者 30楼 发表于: 2009-02-11
2.6模块II:运行Visual MODFLOW
下面我们将学习用于Visual MODFLOW的MODFLOW、MODPATH和MT3D的某些运行选项的选择。
2.6.1第10节:流动模拟运行选项
如果你已经完成了模块I中的所有小节或模块II中的其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   Tutorial文件夹
++   airport.vmf
打开文件,出现Visual MODFLOW主菜单。
从主菜单上选择:
+    Run
你就会进入运行选项屏幕,此时会有一个窗口提示你选择稳定(Steady-state)还是渐变(Transient)模拟。
+    ¤Steady state
+    OK
在Visual MODFLOW中,运行选项分为两个独立部分:流动模拟的运行选项和MT3D模拟的运行选项。当进入运行选项屏幕时,Visual MODFLOW在默认条件下将进入流动模拟的运行选项。从顶部菜单上选择
+    Basic
+    Initial Heads
新模型的默认初始条件是¤Constant by layer(按层给定定水头条件)。Visual MODFLOW将根据每层的标高和边界条件为各层确定一个初始估计值。对于简单的问题而言,这一般都会得到一个收敛的解。然而,一个好的估计值往往可以显著地减少找到收敛值所需的迭代次数。注意,用于输入初始水头估计值的文件可以有多种格式可供选择。因为我们要做的是个简单问题,所以用¤Constant by layer做初始水头估计就已经足够了。
+    OK
此时就会出现一个初始水头条件(Initial Heads)窗口,上面列出了模型中每一层所用的初始水头值。它们应该都是一样的。
+    OK
现在,检查一下Visual MODFLOW所提供的数值求解器的选项。
+    Solver
默认求解器是WHS Solver,这是由Waterloo Hydrogeologic Inc.开发并拥有的求解器。这是目前最快、最稳定的MODFLOW求解器。在本例中,用WHS Solver来求流动解。
+    OK接受
此时会出现一个WHS求解器参数设置(WHS Solver Parameters)窗口,其中列出了所有默认的求解器设置。
+    OK接受
还有其它一些流动模拟的运行选项,如果你有时间的话,可以看看这些本教程中未提到的选项。
2.6.2第11节:MT3D v.1.5和MT3D96模拟的运行选项
如果你已经完成了模块I中的所有小节或模块II中的其它小节,你可以跳过框中的步骤。
+    File
+    Open
出现一个文件选择窗口。
++   Tutorial文件夹
++   airport.vmf
打开文件,出现Visual MODFLOW主菜单。
+     Run
+     ¤Steady state
+     OK
从顶部菜单上选择:
+    MT3D
进入MT3D运行选项屏幕。再从顶部菜单上选择:
+    Advection
此时就会出现一个如下图所示的对流参数设置(Advection Parameters)窗口。
+    ¤Method of Characteristics(MOC)
特征值法(MOC)和改进的特征值法(MMOC)都能求得污染物运移的稳定解,但是MMOC有很大的数值弥散现象。在模块III:输出的可视化中我们将对这两种方法作一个比较。
+    OK接受其余的默认参数
此时将返回到MT3D运行选项屏幕。
+    Output/Time Steps
出现一个输出和时间步控制(Output and Time Steps Control)窗口,如图所示。
MT3D输出控制的默认设置是仅保存模拟结束时的一组浓度值。此处我们要作一些修改。首先,单击þSave simulation results at the end of simulation only(仅模拟结束时保存一次模拟结果)的前面的框,去掉ü,即不选该项,然后填入并核对下列各项:
Simulation Time:    7300 daysÖ
单击标明Max # of transport steps(最大运移步数):的框
type:    10000Ö
+    ¤ specified times
在¤ specified times(指定时间)标签下的框内输入下列数值:
1825¿
3650¿
7300
+    OK接受这些参数
以上就是MT3D的运行选项部分的教程。需要再次说明的是,还有许多本教程未提及的运行选项,如果你有时间的话,可以看看这些选项。
虽然本教程可以让模型进入运行,但可能由于时间的关系你还不想就在此时运行。下面将向你显示模型是如何开始运行的。
从上部菜单栏中选择Translate/Run,将出现运行设置(Run Setting)窗口,其中有以下一些可设置的项目:
选择恰当的用于模拟的数值引擎,Visual MODFLOW 将把程序输入部分提供的信息转换成带有合适扩展名(.BAS, .WEL, .BCF等)的标准MODFLOW、MODPATH和MT3D输入文件。然后Visual MODFLOW将自动地执行模型计算分析。模拟完成,得到一组收敛解后,它就会返回Visual MODFLOW的主菜单。注意:运行Visual MODFLOW的MT3D1.10时,运行时间可能很长。对于不同的计算机运行所花时间大致如下:
                              MT3D 1.5               MT3D 96
l    486                      10-12小时                3-4小时  
l    奔腾100                  5-6小时                 1.5-2.5小时
l    奔腾 150                  3-4小时                 1-2小时    
+    OK      开始运行
如果不想运行该模型而是继续我们的教程,就从底部菜单栏中选择F10-Main Menu按钮。
如果你选择了运行模型,Visual MODFLOW将启动Win32 MODFLOW程序组。在本手册第四章:运行选项中你可找到有关Win32 MODFLOW程序组的信息。
2.7模块III:输出的可视化
如果你已运行了模型,你可以将你的模拟结果可视化显示,如果你没有运行模型,你还可以将包含在Visual MODFLOW安装盘中的AIRPORT模型模拟结果可视化显示。
如果你已经完成了模块I和II并且运行了模型,你可以跳过框中的步骤。
从主菜单:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+     OK
+     Output
此时转到Visual MODFLOW的输出屏幕上。在默认情况下,输出屏幕将显示不带颜色填充的等势线图,如下图所示。
从顶部菜单上选择Contours就可以看看程序能绘制哪些等值线。其下拉菜单会列出以下一些可以选择的等值线。
选择Head Equipotentials回到等势线图上。
2.7.1第12节:等势线和等值线绘制选项
如果你已经完成了模块I和II并且运行了模型,你可以跳过框中的步骤。
从主菜单上:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
要选择等水头线的等值线绘制选项:
+    Options(在左边的菜单条上)
就会出现一个如图所示的等水头线绘制选项(Equipotential Overlay Contouring Options)窗口,但没有右边的颜色充填条。要激活颜色充填:
+    þ Colour shading
在标明间隔(Interval:)的文本框中,将值从0.5改为0.25。在标着标注(Labels)的框中,将小数位数改为2。
你也可以通过单击标着等值线生成精度/速度(Contouring resolution/speed)的按钮来改变等值线的绘制速度。每点一次,等值线的绘制速度增大2倍。在本次练习中,按钮上应该是High/Slow。
+    OK接受
等值线的位置应该与前一个图非常相似,但由于等值线的精度降低了,这些线会有一些轻微的摆动。上图没有加入颜色充填。
其它没有通过菜单激活的等值线选项可通过在模拟区内按下鼠标右键来激活。此时会有一个带有增添、删除、移动等值线和标注的选项的弹出窗口(见下图)。
选择Add Contour选项,将鼠标移至模拟区的任意一处并单击鼠标左键。在你鼠标点下的位置就会加入一条等值线。要再加入一条,只需再点一下鼠标即可。要再次弹出等值线选项窗口,只需简单地再按一下鼠标右键即可。这一次我们选择Move Label,然后将鼠标移至你想要移动的等值线标注的位置上。在标注上按下鼠标不放,然后将其沿等值线拖至想放的位置,放开鼠标,这样就设定了标注的新位置。标注都移动完后,再按一下鼠标右键。
下面,让我们看看场地的剖面情况。从左边菜单上选择View Column,然后将鼠标移至模拟区内。在模拟区中间位置按下鼠标左键选择一列。
要去掉用户设置上去的等值线,在剖面上单击鼠标右键弹出等值线绘制选项窗口。
+    Delete all custom
+    View Layer回到平面模式
在进入下一节前,必须关掉色彩充填功能。在左侧菜单条上
+    Options
+    Colour shading
+    OK
2.7.2第13节:流速矢量及其选项
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
要察看流速矢量选择:
+    Velocities
你将进入流速矢量输出选项屏幕,见下图。流速矢量将根据其默认的设置进行绘制,即矢量相对大小与流速的大小一致。
要让流速矢量与大小无关
+    Direction(左边菜单上)
这样所有的箭头的大小都将是一样的,它们只代表了流动方向。
如果想改变显示的矢量箭头数目
+    Options(左边菜单上)
并输入:
Vectors:    40
这就指定了一行上的流速矢量的数目。
þAutoscale
¤Variable Scale
+    OK接受
注意颜色:在默认情况下,红色表示朝外(也就是说,从层的角度看是向上的)的流速,蓝色表示向里(也就是说,从层的角度看是向下的),绿色表示与平面平行(也就是说,从层的角度看是水平的)。
再让我们看看场地的剖面情况。从左边菜单上选择View Column,然后将鼠标移至模拟区内。单击鼠标左键,选择任意一列。就会出现一个与下图类似的在剖面上同时显示等势线和流速矢量的图。
为显示水力传导系数图层的离散化情况,从底部菜单栏中选择F9-Overlay按钮,将出现一图层控制(Overlay Control)窗口,该窗口按字母顺序列出所有可以打开或关闭的图层。把鼠标置于Conductivity Overlay(水力传导系数图层)之上并双击鼠标左健。这样就会在Conductivity Overlay旁出现一个星号(*),意味着它被激活了。
+    OK 显示水力传导系数图层剖面。
从左边菜单选择View Layer,单击模型横剖面中的第一层,返回到模型平面视图。
要从屏幕中消除流速矢量,从底部菜单栏中选择F9-Overlay按钮,出现图层控制窗口后,把鼠标置于Velocity Overlay上并单击使之高亮度显示,再单击标有ON的按钮使之变为OFF。这样Velocity Overlay旁的星号将会消失,意味着它已无效了。
+    OK  显示没有流速矢量的窗口
2.7.3第14节:流线及其选项
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
+    Pathlines
默认条件下的屏幕将显示所有流线在模型当前层上的投影(如下图所示)。
要显示当前层的流线,从左边菜单栏中选择Segments按钮 。这样图上就会只显示模型当前活动层的流线。
用Next按钮可以顺序向下显示模型的每一层来查看质点流线的位置。
+    Projections   观察所有流线。
注意流线上标有箭头表明流向。这些箭头也作为时间标志表示一质点到达某一目的地的时间长度。要确定每一时间标志的时间间隔,从左边菜单栏中选择Options,此时将显示一流线选项(Pathlines Options)窗口,如图所示。
注意:þUse time marks选项是有效的。在标有¤Regular every ‘1000’ days方框中显示了标志时间间隔。
为了显示在10,000天内流线能走多远;
+    ¤ Time Related 并在方框中输入10,000。
+    OK 显示10,000天内与时间相关的流线
现在显示横剖面上的流线。
+    View Column
在屏幕上移动光标,使高亮度条带位于不连续弱含水层区附近,并单击它。
现在我们返回到平面视图显示。
+    View Layer
在屏幕上移动光标使高亮度条带位于第一层,并单击它。
要从图中删除流线,可从左边菜单栏选择F9-Overlay按钮,并向下滚动图层列表直到你看到Particle Overlay和Pathline Overlay为止。在它们上面双击鼠标去掉星号,使这两个图层失效。
+    OK
2.7.4第15节:MT3D浓度等值线
如果你已经完成了模块I和II并且运行了模型,或者来自模块III的某一节,那么你可以跳过框中的步骤。
从主菜单:
+    File
+    Open
出现一个包含一系列Visual MODFLOW输入文件的文件选择窗口。
++   Tutorial文件夹
++   output.vmf
+    OK
+    Output
在这一节中,我们将对特征值法和改进的特征值法进行比较。
+    Contours
+    MT3D Concentration
此时将会出现浓度等值线图,该图显示了第一个输出时间1825天时的浓度情况。因为用对数坐标及颜色充填来表现浓度等值线会更简单一些,因而你可以按如下步骤来自行设计等值线图:
+    Options
¨ Automatic contour levels
þ Use cudtom contour levels  
þ Colour shading
在方框中输入如下用户自定义的等值线间隔:
0.0
100.0
500.0
1000.0
2000.0
4000.0
+    OK
尽管这不是一个严格的对数坐标,但它还是能提供一种更好的浓度等值线的表示法。
下图是用特征值法(MOC)和改进的特征值法(MMOC)所计算出的在1825天时的浓度等值线。
上图是用改进的特征值法(MMOC)所计算出的在1825天后的解。
上图是用特征值法(MOC)所计算出的在1825天后的污染羽。
用MMOC来计算的污染物运移模拟的数值弥散性比MOC要高。尽管这两种模拟在数值上都是稳定的,但我们仍建议用户利用自己的专业知识来判断到底多大的数值弥散是允许的。
这些MT3D模拟指定了三个输出时间:1825天、3650天和7300天。从左边菜单栏选择Time可以查看输出时间的序列。
+    7300
+    OK   绘制7300天以后的浓度等值线图。
下面两图分别显示了用MMOC和MOC所做的第7300天的污染羽分布。
上图显示了用MMOC所做的第7300天的污染羽。
上图显示了用MOC所做的第7300天的污染羽。
在7300天以后,MMOC模拟的数值弥散量比第1825天要大得多。而MOC模拟则表明污染物向不连续弱含水层中的运移情况。
对单个或多个观测点都可以画出其浓度穿透曲线。从左侧菜单条上
+    Select来选择观测井
会有一个弹出窗口提示你选一种选择观测井的方法(Single、Box和Polygon)
+    Box
沿三个观测井用鼠标拉一个框后,这些井就会变为红色,表示它们已被选中。
从左边菜单上
+    Conc. Vs Time
就会显示一个包括每个井的浓度对时间关系图。观测到的浓度在图上以星号表示。下图是MOC模拟的浓度-时间图。穿透曲线的不规则是MOC技术的离散性的结果。对于高浓度,稳定状态的浓度可以通过时间上的平均来求得近似值。而在低浓度的情况下,由于噪声的影响,很难确定一个平均浓度。
作为对比,下图是MMOC模拟的浓度-时间图。其曲线要光滑一些是因为MMOC中增加了数值弥散的结果。下面的稳定浓度也显示出作为使用MMOC的结果,出现了数值弥散。
+    Close
整个教程到此为止就结束了。正如前面所说,这些说明的目的是给你提供一个Visual MODFLOW的主要特怔和功能的介绍。如果你有时间的话,我们鼓励你回过头去试试那些使用说明中还未涉及的特怔和功能。

第三章    输入
本章主要内容如下:
¨    建立一个新的Visual MODFLOW模型
¨    向Visual MODFLOW输入已有的MODFLOW模型
¨    编辑已有的Visual MODFLOW模型
如果你以前没有用过Visual MODFLOW,又想尽快地熟悉Visual MODFLOW的用户界面,我们建议你先看完第一章,再跟着做第二章的练习。
Visual MODFLOW的数据有两种方法生成。我们推荐使用Visual MODFLOW的图形环境,我们将在本节中描述有关内容。但是你也可以文本编辑器创建和修改数据文件。有时,在你对模型已经比较熟悉,并且只需对现有的文件作一点小小的修改时,用文本编辑器可能会更快一些。为方便起见,所有Visual MODFLOW数据文件的结构和形式都有将在第五章中加以描述。然而,如果你使用文本编辑器,Visual MODFLOW并不会去检查数据的一致性。Visual MODFLOW仅监测和检查在图形环境下输入的数据的一致性和合理性。因此,在用文本编辑器修改数据文件时,我们建议你要特别小心。
3.1新模型的建立
如果你正要创建一个新的Visual MODFLOW模型,请继续做下去。但是,如果你要用一个已经存在的模型,请跳过这一节到基本输入那一节。
在Visual MODFLOW 中,要创建一个新的模型的话就在File菜单中选择New即可。输入一个新的文件名。如果操作不成功的话将显示一个错误信息(比如已存在与此相同的文件名)。
然后将显示如下对话框让你去
l    设置网格的默认大小,
l    设置底图(.BMP或.DXF文件),以及
l    选择模型的运算单位。
在对话框中,你必需输入行、列和层数,以及模拟区最大、最小X,Y,Z坐标值。如果你已选择使用.DXF文件的坐标值,则该文件中的最大、最小X和Y坐标值将作为缺省的最大、最小X和Y坐标值。总层数和模型的最大、最小标高也应在此输入。Visual MODFLOW将在此基础上生成一个均匀的网格。而后,象在网格(Grid)一节所描述的那样可以细化网格,但是它的总体大小不变。
注意:运算单位一旦确定就不能再加以修改。
你必须为你的数据选择单位。单位不必和国际单位或英制单位一致;允许使用比如水力传导系数的单位用厘米/秒而抽水量用克/分之类的单位。在你输入运算单位后,选择Create。
3.1.1底图
Visual MODFLOW程序的所有模块都支持使用底图。这些底图既可以是DXF(AutoCAD v14)也可以是BMP格式。在选择模拟区域(Select Model Region)窗口中,对于DXF和BMP文件都会出现如下图标。
       放大:单击你想缩大区域的一角,再单击该区的另一角。
       缩小:缩至显出整个模型。
       新图:弹出一个选择模型底图(Select Map for Model)窗口,你可在此指定一个新的DXF或BMP文件作为底图。
       重设区域大小:允许用户在区域窗口的一角处单击鼠标并拖动它到所需的大                小。双击该区将在Rotate和Resize Region之间切换。
       旋转区域:允许用户在区域窗口的一角处单击鼠标并拖动它到所需的角度。双击该区将在Rotate和Resize Region之间切换。
       对齐区域:按零度角排列区域。
       扩充区域:把区域扩充到整个dxf图。
       把鼠标转化为一般的光标。
3.1.1.1 DXF格式的底图
使用DXF文件作为底图时,你既可以用鼠标也可以通过在文本框中键入值方法来重新设置区域大小。图象的坐标直接从DXF文件中读取,并显示在世界坐标系(World Co-ordinates)的文本框中。你也可以输入和世界坐标不同的局部坐标,它要在网格坐标系(Grid Co-ordinates)的文本框中输入。最后一点,允许用户改变模型的行列数并显示网格。
输入一般模型创建信息后,你将进入输入窗口。
3.1.1.2 BMP格式的底图
使用BMP文件作为底图时,因为文件中没有包含坐标信息,所以这种BMP文件必须具有地理参考系。将会出现如下窗口提示用户如何对图形作地理参考系统。
在你设置模拟区域以前,BMP图形应该有两个具有已知世界坐标的参考点。在图形上单击鼠标,然后在右边活动的文本框中输入坐标设置地理参考点。重复一次输入另一个地理参考点。一旦这两个坐标被接受,其余的文本框就会激活,你也就可以设置网格坐标、旋转角度、行列数,并显示网格。
在建立图形的地理参考点时,会用到以下图标:
      设置地理参考点
      删除地理参考点
使用BMP作为底图时,你既可以用鼠标也可以通过在文本框中键入值的方法来重新设置区域的大小。图象的坐标直接从BMP文件中读取,并显示在世界坐标系的文本框中。你也可以输入和世界坐标不同的局部坐标,它要在网格坐标系的文本框中输入。最后一点,允许用户改变模型的行列数并显示网格。
输入一般模型创建信息后,你将进入输入窗口。
3.2载入原有的模型
在主屏幕的顶部菜单中选择File、Import后,你就可以载入已经存在的MODFLOW或FLOWPATH模型。
3.2.1 MODFLOW
在输入MODFLOW数据时,Visual MODFLOW将尽力从MODFLOW文件的数据中重建你的数据集。然而,由于数据或Visual MODFLOW中的某些特定的设置,下面所说的顺序可能会有所不同。你必须核对你的数据以确保其中所有信息都是正确的。
下面列出了一些目前已知的在载入MODFLOW数据时会存在的限制。
l    层间越流(VCONT):Visual MODFLOW要求输入每一层的明确的层厚和水力传导系数。然后,Visual MODFLOW在转换时生成层间越流项。这将会给那些现有的具有不同越流项的隐含弱透水层的MODFLOW模型带来严重的问题,因为它们的弱透水层的厚度和水力传导系数已合并在越流项之中了。
l    垂向导水系数:在载入MODFLOW数据文件时,.BCF包中的VCONT参数被省去了。在Visual MODFLOW中垂向和水平向水力传导系数之比设为1:10。这个值可在数据文件载入后进行修改。
l    层底标高:如果Visual MODFLOW在数据文件碰到类型为0或2的层位时,就会提示用户输入层底标高。对于顶层,将会提示用户输入层顶标高。如果层底标高有变化将会出现一些问题。你可以把所有层位都设置为类型3,这样在重新生成MODFLOW 的数据文件时就会使数据的损失降到最小。
l    单位储水系数、单位给水度和空隙度:在输入时每层都只能设置一个单位储水系数、单位给水度和空隙度的值。但以后你可以逐单元地修改这些值。
l    洪水演算包: Visual MODFLOW不支持洪水演算包,但是如果存在是不会忽略的(见下文)。
l    弱透水层压缩性包: Visual MODFLOW不支持它,如果MODFLOW 数据文件存在该项也会被忽略。
用Visual MODFLOW运行洪水演算程序
    虽然Visual MODFLOW支持大多数常用的MODFLOW的程序包,但它目前还不支持用户对洪水演算程序的输入数据进行创建和可视化显示。但Visual MODFLOW可以利用用户所提供的适当的洪水演算的输入数据文件来生成合适的运行MODFLOW模拟所需的MODFLOW.IN和FILENAME.BAS文件。
    第四章介绍了如何在模型模拟中包含洪水演算程序的方法。在Visual MODFLOW对输入数据进行转换时,程序对工作目录进行检查,寻找是否存在洪水演算程序的输入数据文件(filename.STR,其中filename是模型文件的名字)。如果存在,Visual MODFLOW就会用这个数据文件来运行MODFLOW(假设它与模型数据文件的其它部分是兼容的)。这个.STR文件必须包含一个与.BAS文件中给出的模型网格及应力期相兼容的数据输入结构。只要模型网格或应力期有什么改变,都必须手工调整.STR数据文件以反映模型网格或应力期的变化。            
如果在弹出菜单中选择MODFLOW,将会出现一个选择文件的窗口,要求你输入MODFLOW模型的.BAS文件名。程序会自动识别.BAS文件中用于特定模拟的程序包,同时程序将读取如下带有标准扩展名的标准MODFLOW文件并输入数据:
                modelname.BAS
                modelname.BCF
                modelname.DRN
                modelname.EVT/.EVP
                modelname.GHB
                modelname.HFB
                modelname.RCH
                modelname.RIV
                modelname.STR
                modelname.SIP
                modelname.WEL
其中所有文件的modelname必须相同。
为MODFLOW模型选择了一个MODFLOW.BAS后,程序会要求你为这个新建的Visual MODFLOW模型选择一个名字。如果你选择了一个和MODFLOW模型一样的名字,在数据转换时将覆盖掉所有的MODFLOW源文件。
此时,输入工具包(Import Utility)将会要求你输入模型的单位。
这类似于创建一个新模型时的单位选择。缺省时间单位直接从原始的.BAS文件中读取,同时与此相应地读入并保存所有的数据。然而用户可以用上述窗口来修改Visual MODFLOW显示水力传导系数和补给的单位。选择好模型单位后,将出现如下窗口,问你模型的基点是否从.DXF文件中读取。
如果你回答是,将出现一个选择.DXF文件的窗口。选择了.DXF文件后,或你在上面的窗口中回答了否,将出现如下的弹出窗口,问你是确认还是输入模型基点。默认值是.DXF的原点或零。
如果你的数据中少了一些文件,输入工具包将会通知你同时中止输入。在试图输入同一个数据集以前,请找到或创建必需的文件并把它们放在.BAS文件的目录。如果你不想创建这些文件或不想把它们包含在输入中,修改.BAS文件使它不包括这些文件的设备号。
此时,输入工具包将要求你输入每层的顶底板标高。如果层的类型是1或3,就不要在此编辑这些数据。
输入层位标高后,程序会要你输入指定的和计算出来的水力传导系数的有效数字的位数。这个值仅在输入时发生作用,它通过减少水力传导系数计算的位数而减少了输入时间。
下一步将显示水力传导系数数据库,如下图所示。生成的水力传导系数是唯一的,Visual MODFLOW将这些值无重复列出。应仔细核对这些值,如果有必要,修改一下。
最后的窗口要用户输入单位储水度、单位给水度和空隙度的缺省值。这些变量的值没有输入,对所有的层位这些值在开始时都将设置为常量。
3.2.2 FLOWPATH 5.x
如果在弹出菜单中选择了FLOWPATH,将会出现一个文件选择窗口,要求你输入FLOWPATH 5.x模型的名字。程序将读入所有具有标准扩展名的标准FLOWPATH 5.x文件并输入数据。然而,因为数据文件或Visual MODFLOW中的某些特定设置,下列顺序可能会有所不同。你必须核对你的数据文件确保其中所有信息都是正确的。
指定了文件名后,出现的第一个窗口将要求你输入单位储水度和单位给水度的缺省值,因为在FLOWPATH 5.x中不使用这些值。
接着,输入工具包将显示由FLOWPATH 文件确定的水力传导系数数据。应仔细核对这些值,如果有必要,修改一下。
因为FLOWPATH 5.2模拟的是二维地下水流,没有垂向分量。显示的下一个窗口要求输入垂向水力传导系数因子。
最后一个窗口显示由FLOWPATH 5.2文件确定的单位储水度、单位给水度和空隙度数据。应仔细核对这些值,如果有必要,修改一下。
3.3基本输入
如果你要编辑一个已有的数据文件,启动Visual MODFLOW后在主菜单选择File项的Open。从文件选择窗口中选择要用的文件。
打开一个已有模型或如上节所述建立了一个新的模型后,从主菜单中选择Input。程序将会转到Flow输入窗口,如下图所示。
在这个窗口内你可以定义三维有限差分网格、井位、属性、边界条件、流线的起始点、模型校验的观测点、分析模型内流量的分区以及附加的注释。另外,MT3D输入屏幕(从Flow输入窗口中进入)可用来输入初始浓度、运移边界条件、模型校验的实测浓度以及弥散系数,如下所示。
创建一个流动模型需要一个完整的数据集,包括网格和边界形状、边界条件和含水层属性。校验资料、流线的设置、区域均衡计算资料和注释是可选的。另外,MT3D运移模拟需要指定初始条件、一个浓度边界条件和弥散系数分布情况。用于模拟校正的浓度实测值是可选的。
Flow输入窗口的顶部菜单栏包含如下输入选项:
File  Grid  Wells  Properties  Boundaries  Particles  Annotate  Zbud  MT3D  Tools  Help
File    保存文件,打印或返回主菜单
Grid    修改网格(默认屏幕)
Wells    插入、修改或删除抽水井和观测井
Properties    增加或修改水力和储水特性
Boundaries    输入或修改流动边界条件
Particles    插入或修改用于质点追踪的质点
Annotate    注释模型
Zbud    插入或修改区域均衡数据
MT3D    切换到MT3D输入窗口
Tools    调用单元监测器和预测工具
Help    获取一般帮助
MT3D输入窗口的顶部菜单栏包含如下输入按钮:
File  Initial Concentration  Boundaries  Calibrate  Dispersion  Chem.Reaction  Flow  Help
File    保存文件,打印或返回主菜单
Initial Concentration    输入或修改初始浓度分布
Boundaries    输入或修改TRANSPORT边界条件
Calibrate    输入或修改TRANSPORT校正数据
Dispersion    输入或修改弥散数据
Chem.Reaction    设置吸附和放射性衰变系数
Flow    切换到Flow输入窗口
Help    获取一般帮助
3.4文件菜单(File)
当你从输入菜单栏单击File时将出现如下菜单选项:
Save    保存数据集
Print    在当前输出设备上打印屏幕并选择打印设备
Print to DXF file    把屏幕打印到.DXF文件
Main Menu    返回主菜单
3.5网格菜单(Grid)
网格设计窗口允许用户定义和离散化模拟区域。用户可以增加网格线或删除单元或者将某些单元排除在计算之外。
当你在输入窗口中选择Grid选项,就会激活一个‘grid’边菜单栏,该菜单具有以下选项:
Edit Columns    打开一个弹出窗口选择列编辑选项,包括增添、删除、移动、输入或输出列。(仅在行或层显示方式下有效)。
Edit Rows    打开一个弹出窗口选择行编辑选项,包括增添、删除、移动、输入或输出行。(仅在列或层显示方式下有效)。
Edit Layer    打开一个弹出窗口选择层编辑选项,包括增添、删除、移动、输入或输出列。(仅在行或层显示方式下有效)。
Inactive Cells    打开一个下拉式菜单设置模型的无效区域,有以下选项:
Mark Poly. Active    把一个多边形区域设置为流动计算的有效单元。按鼠标右键封闭多边形。
Mark Inactive    把一个多边形区域设置为流动计算的无效单元。按鼠标右键封闭多边形。
Mark Single    按鼠标左键使单元无效,右键使单元有效。
Copy Single    按左键复制单个无效单元到其它层、行或列。
Copy Polygon    复制一个无效单元区域到其它层、行或列。
Import Surface    从ASCII文件或SURFER的‘.grd’文件输入一个变化的表面来定义地表或指定层位的底板。同时提示用户输入层的最小厚度和用于内插的观测点数。
X smoothing    在X方向上平滑网格奇点
Y smoothing    在Y方向上平滑网格奇点
3.5.1网格线编辑(行或列)
允许的最大网格是:X方向上500根网线(499列),Y方向上500根网线(499行),Z方向上60层。因此未知水头值的最大数量是15,000,000。如果这个数量还不能解决你的问题,请与我们联系。
如果你的计算机没有足够的内存,Visual MODFLOW将通过Windows的内存管理使用你的硬盘作虚拟内存。
要编辑网格请从边上的工具条中选择Edit Column或Edit Row。此时会弹出一个带有如下选项的菜单:
Add    以图形方式或通过弹出对话框在特定的位置上增加网线
Delete    删除一个以高亮度显示的网线
Move    移动一个以高亮度显示的网线
Refine by    用所选择的因数细化两网线间的网格(比如细化因数是二,将把每行或列等分成两行或两列)
Coarsen by    用所选的因数粗化两网线间的网格(比如粗化因数是二,将把两行或两列合并成一行或一列)
Import form file    从.txt文件载入网线坐标
Export to file    输出网线坐标到.txt文件
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 31楼 发表于: 2009-02-11
      
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 32楼 发表于: 2009-02-11
thanks a lot!!!
离线hengshan168

发帖
3617
土币
3261
威望
6709
原创币
0
只看该作者 33楼 发表于: 2009-02-11
thanks a lot!!!   
离线jjann

发帖
774
土币
1424
威望
1628
原创币
0
只看该作者 34楼 发表于: 2009-02-11
除了感谢,还是感谢!
离线pxb_ruby

发帖
586
土币
88
威望
606
原创币
0
只看该作者 35楼 发表于: 2009-02-13
非常感谢楼主的无私分享
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      浙公网安备 33010602003799号 浙ICP备14021682号-1

      工业和信息化部备案管理系统网站