第一章 快速入门 qD(dAU
本章主要内容为: p v*f]Yzx
¨ 安装Visual MODFLOW 的硬件要求 ]2mfby
¨ Visual MODFLOW 的安装 WHR6/H
¨ Visual MODFLOW的基本运行过程 m>^#:JK
1.1硬件要求 UmP\;
为了运行Visual MODFLOW,你至少需要如下的系统配置: r4pR[G._
l CPU:486DX或奔腾(推荐奔腾) j'LO'&sQ(
l RAM:32MB(推荐64 MB) M^y5 Dep
l 软驱:用于软件安装,高密(1.44MB),3.5’’ `6Yk-5
l 硬盘:至少有40 MB的可用空间 T:VFyby\w
l 鼠标:Microsoft或其它兼容鼠标 Y-P?t+l
l 操作系统:Windows 95/Windows NT 4.0(Service Pack 3)操作系统 kXZG<?
如果在你的系统配置下软件出现了问题,请先确信你的软件是否严格按照安装说明来安装的(见安装Visual MODFLOW)。如果问题仍没有解决,请参阅设置或疑难解答,同时与你的硬件商联系。最后请阅读本章最后部分了解一下如何与Waterloo Hydrogeologic公司进行联系。 jY$3
1.2安装Visual MODFLOW 10..<v7
Visual MODFLOW必须安装在你的硬盘或网络盘上运行。请阅读本章开始有关硬件要求部分,在安装以前确保你的系统满足这一要求。用以下可执行文件安装Visual MODFLOW: B"Ttr+
Setup32.exe 在Windows 95/NT下安装Visual MODFLOW for Windows e8hwXz
这里讲述的安装过程假定Visual MODFLOW是从A盘(源盘)安装到C盘(目标盘)。 ,Y27uey{wa
缺省安装目录为\vmodnt ,但是也可以修改它。 cG4$)q;q
注意:安装盘必须去掉写保护 C<q@C!A
1.2.1 Windows 95/Windows NT 4.0下的安装过程 Z:<6Ck
按如下步骤在Windows 95/Windows NT 4.0下安装Visual MODFLOW : 0t0m?rVW
[1] 插入1号软盘。 e`t-:~'
[2] 进入Windows。 HAGpM\Qa
[3] 单击开始按钮,然后单击弹出菜单的运行项并在对话框中键入A:\Setup.EXE。Visual MODFLOW将开始安装。 tQ(gB_
[4] 单击OK按钮继续。 `j(-y`fo
[5] 键入安装路径或单击OK按钮选择缺省目录开始复制Visual MODFLOW文件。 lm(k[]@
[6] 安装完成以后,单击OK按钮完成安装。 )uH#+IU
现在在你的屏幕上就可以看到一个名为“WHI Software”的程序组,窗口内会有一个VmodWin2.70的图标。你现在就可以使用Visual MODFLOW了。要从你的系统中删除Visual MODFLOW,只要双击Uninstall VmodWin2.70图标即可。它会从你的计算机中删除所有的程序文件,但不会删除工程文件。 Qr]`flQ8
1.3启动Visual MODFLOW ),U>AiF]
安装完毕后,双击“WHI Software”程序组中的Visual MODFLOW快捷图标或依次打开开始/程序/WHI Software/VmodWin2.70,并单击它就可以启动Visual MODFLOW了。你第一次启动Visual MODFLOW时将出现一Visual MODFLOW注册窗口,输入你的公司名和软件盒上提供的系列号。以后你的公司名称将出现在你的Visual MODFLOW打印输出中,系列号将出现在帮助/关于菜单中。 %8! }" Xa
单击Next按钮以识别你所输入的资料,一旦识别成功,你就不必再次输入了。 l#!6
tw+e?
1.4菜单简介 oz?pE[[tm
Visual MODFLOW的集成环境由四个主要窗口组成:主窗口、输入窗口、运行窗口、输出窗口,再加上系统设置和帮助窗口。主窗口包括如下菜单选项: u*<knZ~ty
File Input Run Output Setup Help 8Rd*`]@[pk
File 选择一种文件功能或退出Visual MODFLOW。 2!g7F`/B
Input 进入输入窗口修改当前的Visual MODFLOW数据文件。 >WYradLUi
Run 进入运行窗口修改工程或批处理模式的运行时间参数和运行数值模型。 r(
8!SVX
Output 进入输出窗口数值模拟结果做后处理。 ]"-c?%L
Setup 选择你想用的数值引擎。 c>yqq'
Help 关于Visual MODFLOW的主要帮助窗口和一般信息。 LVxR*O
你的Microsoft 兼容鼠标键功能如下: +jGSD@32>
左键: 这是一个常规“单击”键。在输入框上部按下它并拖动鼠标,字符将变为高亮显示,可为后面输入的字符覆盖。 y[I)hSD=
右键: 该键根据情形的不同有不同的功能。例如,在设置属性和边界条件时,它用来闭合多边形或结束一条线。在网格设置时,它用来精确地定位网格的行和列。 >Ef{e6
当要输入数据和文本时,数字和字母键才有效;所有其它键在Visual MODFLOW中都是无效的。 |E9iG
启动Visual MODFLOW以后,你将进入主窗口,要继续进行你的工作,你应该单击File菜单,它有六个子菜单: (z.eXo P@>
1、 单击Open,打开一个已经存在的文件;或 rCgoU
xW`
2、 单击New,键入一新文件名创建一个新的数据文件; J(=io_\bO
3、 单击Save as,以不同的文件名保存当前模型; :6y;U
4、 单击Import,装入一个已经存在的MODFLOW模型或通过下拉菜单底部的文件路径来选择一个已经存在的文件; rnS&^
5、 单击Export,以AutoCad (dxf)或增强的Windows元文件格式(emf)输出文件。 b=Oec%Adx
你可以双击文件名或单击文件名然后按下回车键来选择一个文件。 RP~ hi%A
为了获得更多的信息请查阅第三章与创立一个新的数据文件,修改一个已经存在的数据文件和装入一个MODFLOW模型有关的部分,或按照第二章一步一步地进行。 F4
:#okt
1.5窗口概览 h][$1b&B
打开一个文件并选择输入、运行或输出,将出现类似于下图的窗口。 Zp
<^|=D
上部菜单按钮: 提供选项,这将根据情况不同而不同 qfl #ki`,
边部选项按钮: 包括视图选项和特别针对当前窗口或模型的功能选项。视图选项如下: 1_6oM/?'
View Column 按列方向上显示一剖面图 Tt\G y
View Row 按行方向上显示一剖面图 ?S~j2 J]
View layer 从剖面图切换到平面图 !ui:0_
Goto 显示某一特定的行、列或层 zJX _EO
Previous 显示前一行、列或层 }q D0-
Next 显示后一行、列或层 [7@9wa1v!
立方体导航块: 用十字丝表示的模拟区一个简单的3D示图,用于确定空间位置。 Vs@H>97,G
坐标区: 显示当前光标所在的实际世界坐标,并显示当前光标所在的单元。 U6yZKK
功能按钮: 输入、运行和输出窗口的通用功能。 o-r00H|
F1 特定窗口或当前区的帮助。 'fr~1pmx#3
F2 为未来的3-D选项保留。 |++\"g
F3 以原文件名保存文件。 K *xca(6
F4 弹出一个文件选择窗口,你可以选择输入一.DXF格式的图形文件,该图形可以在图层窗口上打开或者关闭。 B"=w9w]
F5 按鼠标左键让你满屏放大一窗口。按右键将弹出一个窗口允许你满屏放大一特定区。 V"DilV$v
F6 使屏幕重新显示整个模拟范围。 fsz:A"0H
F7 在模型上单击鼠标左键,拖动该点到新的位置后再单击左键。 8a,pDE
F8 在按行列显示时指定垂向放大倍数。 { bD:OF
F9 图层窗口允许你打开或关闭各种图形。 Auk#pO#
F10 返回主窗口。 vG'I|OWg
状态栏: 每个按钮的功能和使用描述 Z[?zaQ$
CAD(计算机辅助设计)环境是Visual MODFLOW的最大特征之一。它让用户轻松快速地设置复杂问题,同时允许用户监视输入,因此避免了在数据输入过程中经常产生的许多错误。 mo
CAD环境必须要有一个Microsoft(或兼容)鼠标。在所有的CAD环境中以下的键都是有效的: -wrVhCd~g]
功能键 如上所述。 mq6TwM
回车键 选择当前高亮度显示区。 9HBRWh6
Esc键 取消当前窗口。 s?~lMm' !
光标键 有限的功能(如选择菜单项)。 =@y
?Np^A
Tab键 单元之间移动。 %8$ldNhV
1.5.1图层功能 |cIv&\ x
Visual MODFLOW的显示是以一系列层的方式来组织的。各种特怔,如网格线、观测井和导水系数分区的描画是独立的,并且可以独立地打开或关闭,使你的图件或屏幕具有更好的可读性。这种不同的层称为图层(Overlay),它可以让你对数据进行交叉检验。 cPbAR'
带有你的计算结果的图层只能在输出(Output)屏幕中使用。除此之外,输入(Input)和输出(Output)的图层是一致的。 $BUm,
当你在主菜单、输入、输出和运行窗口之间切换时,这些图层都有一些特殊的缺省设置。在主菜单窗口中仅显示模型边界的底图.DXF和.BMP图层。当你首次进入输入窗口时自动加上网格(Grid)图层,但当以输入模式进行操作时,缺省的图层就和前后操作有关了。在输出窗口中,缺省的图层也和前后操作有关,但是当你操作时以前的图层依然是可见的,除非你明确地将其关闭。 SKeX~uLz
按F9键或用鼠标单击屏幕底部的F9按钮,将显示图层对话框。 -0f,qNF
各种各样的图层都可以让你对数据进行交叉检验,或者从你的图件中增删一些特性来使用你的图件更具可读性。输入的图件也可以从图层窗口中加以删除。 "#G`F
你可以用设置(Settings)按钮更改某些图层的显示特怔。分区类的特怔如导水系数和补给,既可用填充区块也可用外廓线的形式显示。DXF图件的上色以及计算所得的特怔如等值线和流速矢量的选项,都可通过设置按钮来设置。 (,OF<<OH
图层在表列里出现的顺序将决定图层在屏幕上的描绘顺序。描绘顺序与表列里的顺序相反(也就是说表列中最前的特怔是屏幕上最后才描绘的特怔)。将Default按钮切换为User Defined后,图层表列可以用图层控制对话框右下方的四个按钮来作调整。 TAp8x
1.6 Visual MODFLOW中的打印功能 =u
3YRqz
Visual MODFLOW可以让你打印当前模型的资料。从主窗口、输入、运行或输出窗口中选择File,再选择Print。你会看到如下对话框: Ze"m;T
Project Information(项目信息)指图件标题栏的文字信息。当Full Image(全图像)被选中时,就不会打印标题栏。 +\)a p
对一个连贯的图形窗口,在选择放大(Zoom In)后可用鼠标右键为你的模型指定视图的区域。该特怔取得顶底的网格坐标后再对左右边界的网格坐标进行计算,形成一个矩形视域。它还可以让用户转回到上一次缩放窗口或改变当前的缩放窗口。 j`:D BO&)\
程序将会自动地选择缺省的Windows打印机为当前的打印设备。要改变当前的打印机,单击Select,你会看到如下对话框: l _g JC.
通过选择打印机名称(Name)输入框旁的箭头,你就可以将当前打印机切换到任意一台目前已安装的Windows打印机上。选择Properties会有一个标准的Windows打印选项对话框,可让你调整纸张、字体和设备的设置。 z[+Sb;
1.6.1打印到DXF文件 =:;K nS
如果你要将结果输出到DXF文件中以便于在诸如AutoCAD之类的CAD程序中进行后处理,你可以将你的结果“打印”到DXF文件中去。 Df5!z \dx
要将结果输出到DXF文件中,先从主菜单上选择File再选择Print to DXF file。下面的对话框将提示你输入DXF的文件名。 k<}3_
1.7设置 HvLx
我们用Setup选项为不同的Visual MODFLOW软件包选择不同的数值引擎。要进入设置,点击顶部菜单的Setup,就会出现一个下拉式菜单。 _:7:ixN[Ie
1.7.1数值引擎 O
8XHaVLg3
Visual MODFLOW允许用户为MODFLOW和MT3D选择不同的数值引擎(FORTRAN可执行程序)。选择Setup后,从下拉式菜单中选择Numeric Engines,就会出现数值引擎(Numeric Engines)对话框。 /_qW?LKG/
MODFLOW的数值引擎只有一个即USGS MODFLOW 96 from WHI。 J^y?nE(j
MODPATH和分区均衡都只有一个数值引擎:分别是USGS MODPATH from Waterloo Hydrogeologic和 Zone Budget from Waterloo Hydrogeologic。 ]8/g[Ii
运移的数值引擎有四个可供选用: \qz! v
l EPA MT3D,v.1.1(公共版权); V_/.]zQA
l DOD MT3D,v.1.1(公共版权); dCpDA a3
l MT3D, v.1.86 (Papadopulos & Associate Inc.); W^YaC
(I
l MT3D96 (Papadopulos & Associate Inc.)(只有在单独购买了MT3D96后才可用)。 ]\m>N]P]
如果MODFLOW运行后MT3D的数值引擎发生了变化,就要重新运行MODFLOW,更新projectname.flo 文件,使之与MT3D的数值引擎保持一致。 yS1i$[JV
1.7.2 Windows95/NT 交换文件 X.+|o@G
Windows95/NT管理系统的虚拟内存。要改变这些设置请询问你的系统管理员。 / }Pj^^6A<
1.8疑难解答 y;<F|zIm
1.8.1常见问题 7j8lhrM}^
问:为何Visual MODFLOW有时反应迟缓? 6uIgyO*;k
答:在数据量很大时,如果Visual MODFLOW将RAM用完,它就会占用你的硬盘作为临时空间。对许多操作而言,将不用的数据移入移出临时空间增加了运行时间。Visual MODFLOW拥有的RAM越多,对临时空间的需求就会减少。 aM+Am,n`@
另外,你可能需要为计算机购买一些RAM了。 >!Dp'6
问:OUT OF MEMORY- CANNOT CONTINUE 这条错误信息是什么意思? PUmgcMt
答: Visual MODFLOW已经用掉了运行你的数据所需的内存。增加这个值就可以让Visual MODFLOW占用一些硬盘作为不常用的数据的临时存储空间。请注意,在运行前,你的硬盘上至少有这么多KB的自由空间。 bs0[ a 1/
问:为何在和模拟区比较时我的DXF图是HUGE(巨型)? >R}G
答:Visual MODFLOW要求你的DXF图形单位是十进制。你要输入的DXF文件可能是用英尺和英寸这种工程或建筑尺度。为了准备Visual MODFLOW 的DXF文件,你应该先用AUTOCAD或其他恰当的CAD软件编辑一下你的DXF文件。首先,你应该把单位转变成DECIMAL形式。这将把你的图形转变成十进制英寸形式。然而,这也将把你的图形文件放大成12倍。为了弥补这一点,你需要用SCALE命令把整个图以0.0833333(1/12)系数缩小,缩放完成之后,你需要用DXFOUT或用类似的命令把整个图形再保存一次。文件现在就可以以正确的尺寸输入到Visual MODFLOW中去,Visual MODFLOW的未来版本将可以识别这种情况,并自动的缩放图形。 Jq?Fi'2F%
问:我想重编译USGS MODFLOW。我的Visual MODFLOW还能用吗? #Q'#/\5
答:Visual MODFLOW使用Watcom F77编译器来编译MODFLOW源文件。另外,MODFLOW源文件还可以用Lahey F77编译器编译,它可以和Visual MODFLOW兼容。应该知道,用Fortran 编译器编译的非格式化二进制文件是不必兼容的。如果使用其它编译器,在对.HDS和.DDN进行后处理时可能会出现问题。如果你想使用其它编译器或在快速工作站运行MODFLOW,请向我们询问有关如何读取你特有的二进制文件的Visual MODFLOW用户定制的修改问题。 Mj`g84
问:Visual MODFLOW不能创立必要的XXXX Kbytes的交换文件,为什么? 72 ZoN<c
答:MODFLOW的数据文件可能非常大,经常需要比你的PC机所能提供的RAM存储空间还要大。然而,在MODFLOW处理你的问题时, Visual MODFLOW通过利用你的硬盘空间来暂时存储不需要的数据的方法来产生更多的“内存”。但是,根据你的程序大小,当你使用虚拟内存时,如果不能得到所需的总内存,MODFLOW可能仍然不能处理你的数据文件。 uHq;z{ 2GI
交换文件的大小等于虚拟内存减去可用内存,如果它大于可用的硬盘空间,那么Visual MODFLOW就会给出这条消息。 -2'1KAk-W
在这一点上你有三种选择: 055C1RV%
1. 尽力按正常去运行模型。有时,某个问题可能并不需要虚拟内存参数规定的内存的大小。 8SA"
bH:
2. 按指定空出硬盘空间。 6aLRnH"Ud
3. 在 vmod.vmc文件中把虚拟内存参数改小。然而请注意,如果你的问题所需内存大于规定的值时将产生一个错误消息:DOS4Gpro Error 1307¾not enough memory when attemping to run MODFLOW。 +Fb+dU
1.9如何获得技术支持 % {-r'Yi%
如果你有任何问题,请随时和Waterloo Hydrogeologic公司联系。我们对已注册使用Visual MODFLOW的用户提供免费技术服务。在简介那一章中你可以找到我们的地址。 _L+j6N.h1
一般情况下,我们会在一个业务天中回答所有技术服务问题,但是可能要两个业务天来作答复。 0n}v"61q
在信中除了简单描述你的问题外请写上你的:
ne:
'aq
l 名字 f,x;t-o+R
l 公司和办公室 }U '
l 电话和传真号码 YSrjg|k*
l email地址(如果你有的话) _F8THYg (
l 产品序列号 [bw1!X3
l 版本注册日期 aWPf3Q
你的版本注册日期可以在主菜单窗口中的HELP菜单的About Visual MODFLOW子菜单中找到。 $[Q;{Q
如果你的问题是“我要如何去做...”之类的话,附上这些资料就足够了。但是,如果你是报告一种故障(特别是导致系统崩溃的故障)或有关该软件的技术问题,请另外提供如下信息: "
sC]z}
l 计算机的硬件配置。 ?)T@qn+
l 操作系统(Win 95 /WinNT)。 yi:}UlO
l 计算机的RAM大小。 Fv*Et-8tN5
l 计算机的剩余硬盘空间大小。 1R;@v3
l 详细描述问题是在什么时候和如何发生的。 &u~Pp=kv
如果我们不能重现你所遇到的问题或不能解决此问题,我们将要求你给我们寄一份你的工程文件的副本。 "j&'R#$&d
sF>O=F-7
第二章 使用教程 S*gm[ZLQ
本章是Visual MODFLOW的使用教程,它用一个具有完整的输入、输出文件的模型例子(Airport)来让你熟悉后处理能力和Visual MODFLOW的功能。该问题的数值模拟结果(MODFLOW、MODPATH和MT3D)已经完成,你可以用它来评价一下该例模拟结果的可视化输出能力。 N\];{pe>
该章将向你显示设计并运行模型以及如何将结果可视化显示的一些必要的步骤。本教程是以按一步一步的方式组织的,你可以选择感兴趣的章节而不必做完全部练习。 XhA tf@n
2.1实例模型介绍 ltlo$`PR
本场地位于Waterloo城外的一个机场附近。场地表层地质情况是:上下各有一层砂砾含水层,中间有一粘土和粉砂弱透水层将它们分开。相关的场地地物有一个飞机加油场、一个生活供水井场和一块弱透水层不连续区(含水层天窗)。如下图所示。 ,,!P-kK$
生活供水井场有两眼井。东井抽水量稳定在550m3/d,西井为400 m3/d。在过去十年中,在加油场中不断地有飞机燃料的泄漏,其自然下渗已在上层含水层中产生了一个污染羽。本教程将向你显示如何一步步地为这个场地建立一个地下水流动和污染物运移的模型。这个模型将显示燃料污染物对生活供水井的潜在影响。 >bd@2au9!
在讨论时,我们将该场地在平面上的视图指定为上北下南左西右东。如图所示,地下水在一个由上部潜水含水层、中间弱透水层和下部承压含水层所组成的三层系统中自北向南(平面图上是从上到下)流动。上下含水层的水力传导系数为2e-4m/s,弱透水层的水力传导系数为1e-10m/s。 iCLH
2.2有关的术语和记号 D(OJr5Gg
在本教程中,将用到以下的术语和记号: BeN]D
TYPE 输入给定的词或值 J?EDz,
SELECT 在所指的地方按下鼠标左键 >JAWcT)d
« 按下〈Tab〉键 o2'Wu:Y"
¿ 按下〈Enter〉键 c&I
+ 在所指的地方单击鼠标左键 no W]E}nN
++ 在所指的地方双击鼠标左键 \m%c"'[
2.3启动Visual MODFLOW t0GJ$])
在Windows环境下,你只需简单地双击WHI Software程序组里的Vmod图标,即可启动Visual MODFLOW。 |j81?4<)v
现在你就可以看到Visual MODFLOW的启动屏幕。要进入主菜单: -Lu&bVt<>
+ OK m"'}{3$%
2.4打开文件 !\$V?*p7
在本教程中,你可以从任意地方开始,在章节之间跳来跳去,或者离开,或者从你停止的地方开始。要创建一个新文件,请跳至模块I。要打开一个已有的模型,在顶部菜单条上: ]Vmo>
+ File TQmrL
会有一个下拉菜单提示你从一个标准文件选项的表列中进行选择。 T2/:C7zL
+ Open - U\'Emu4
将会出现一个文件选择窗口。 #]6{>n1*+w
++ Tutorial文件夹 Uqb]&2
如果你要开始学习第2节: ~\AF\n%
+ Grid1.vmf d#]hqy
如果你要开始学习第3-7节: awwSgy
+ Grid2.vmf 7F6B
如果你要开始学习第8-10节: Mn/@?K?y
+ Airport.vmf hl7 z1h
如果你要开始学习第11-14节: K,I
+ Output.vmf g&85L$
如果你要继续本教程: `CH,QT7e
+ filename.vmf(你在第1节创建的文件) ~bCA8
+ OK打开文件 o{! :N> (
几秒钟后,主菜单就会出现,同时显示出模型网格和场地地图,如下图所示。 XL^05
2.5模块I:模型输入 A+gS'DZ9C
2.5.1第1节:设置数值引擎 6o7t eX
在进入输入屏幕以前,很重要的一项工件是选择运行模型所需的数值引擎。从顶部菜单条上: Ei):\,Nv
+ Setup &e@)yVLL
+ Numeric Engines B~N3k
新模型的缺省数值引擎如下所示: !&>`
重要说明: nvT@'y+
Visual MODFLOW包括免费的MT3D的公共版本(EPA的1.1版和DOD的1.5版)。但是1.10版没有1.86版稳定,它们都没有MT3D96那么快而灵活。本教程的飞机场模型练习是用MT3D v.1.10和MT3D96做的。对那些需要精确模拟污染物运移的专业咨询公司来说,所有的模拟都应该用MT3D96来做。 4k-Ak6s
要换一种数值引擎的话,单击你想要的改变的那种数值引擎的选择框旁边的向下箭头键,然后在下拉对话框中选择。对于本教程中的模拟练习,你可以使用DOD1.5或MT3D96,如果你已分别购买了的话。 L/r_MtN
+ Transport Numeric Engineer xE
w\'tH
+ MT3D by Papadoplulos_Associates Inc. [#q]B=JB
+ OK eu9*3'@A
2.5.2第2节:生成模型网格 p/?o^_s
第一个模块将教你如何利用Visual MODFLOW模拟环境通过必要的步骤生成一个新的模型数据文件。 {\(MMTQ
+ File nAEyL+6U
+ New [GI~ &
此时将出现一个创建新文件的窗口。 -]HZ?@
++ 在Tutorial文件夹中 (X,i,qK/
在文件名输入框中输入一个新的数据文件名(比如,键入‘Tutorial’),然后选择Save,就生成了一个新的数据文件。 jt on \9
Visual MODFLOW会自动在文件名后加入vmf扩展名。 ~cCMLK em
然后,将出现一个一个选择单位的窗口。在下图中用鼠标为你的特定系统指定度量单位。 P+}~6}wJE
+ 3Create model using base Map 5EcVW|(
下一步你必须说明.DXF底图文件的目录和文件名。 !43!JfD
+ Browse v={{$=/t
在你安装Visual MODFLOW的目录(缺省为:c:\VMODNT\tutorial)下找到并选中下面这个文件: &]#L'D!"
+ sitemap.dxf ^ls@Gr7`P
+ Open 3@Mh* \;\b
下一步你应该输入模型的行数、列数以及层数。 0H}tb}4
Enter Number of Columns (j)(列数) 40Ö <[cpaZT,
Enter Number of Rows (i)(行数) 40Ö P,{Q k~iu
Enter Number of Layers (k)(层数) 3Ö (Z(S?`')
Enter Zmin(最小的高程值) 0Ö 'UIFP#GtFO
Enter Zmax(最大的高程值) 18Ö 31a,i2Q4
在上述窗口的Units(单位)部分,用鼠标从它们各自的下拉列表选择如下选项: *uhQP47B
Length(长度) meters M}x%'=Pox
Time(时间) days Giv,%3'
Conductivity(传导系数) m/second GXVx/)H
Mass(质量) kg #f2k*8"eAF
Concentration(浓度) mg/L ,Gt!nm_
Pumping Rate(抽水速率) m3/day bTc>-e,
Recharge(补给速率) mm/year FN-/~Su~J
+ Create (确认这些设置) q>.7VN[
vE
之后,将出现一个模型区域选择(Select Model Region)窗口, 你可以在这里定义模拟区的范围。Visual MODFLOW会从sitemap.dxf(场地底图)中读取坐标的最大值和最小值,并默认它是处于模拟区域的中心。现在的Visual MODFLOW允许用户转动dxf图像、使用局部坐标系、设置DXF图的范围。如果用位图作为底图,则可以对图形进行地理参照处理和旋转。在屏幕中的原来的数字上输入以下数值: -[L\:'Gp5
Enter Start Point X 0Ö q()o|V
Enter Start Point Y 0Ö 'F/oR/4,
Enter Angle 0Ö _xaum
Enter Grid Co-ordinates X1 0Ö {Ya$Q#l
Enter Grid Co-ordinates Y1 0 Ö K=[7<b,:3
Enter Grid Co-ordinates X2 2000 Ö cERIj0~
Enter Grid Co-ordinates Y2 2000Ö G<$:[ +w
+ OK 接受网格尺寸 })/P[^
在模拟区内将会自动生成一个大小均一的4054053的有限差分网格,同时在屏幕上将会出现一个如下图所示的场地底图。 5qtk#FB
从这儿你可以进入输入部分的任何一节,看看Visual MODFLOW是如何处理各种输入参数组的。 x|rc[e%k
前面已经提到,如果是想看看结果的显示,是没有必要将本教程所写的所有步骤全部运行一遍的。本教程已提供了每个模块所需的所有文件。 G&g;ROgY
2.5.3第3节:细化模型网格 gu(:'5cX
本节介绍了对感兴趣的地段如供水井周围地区、含水层天窗区与加油场之间地段进行模型网格细化的步骤。由于单元的大小代表了井的尺寸,因而用更细的网格来模拟井的影响将会更接近实际情况。另外,如果在井周围出现了水位下降,用小的网格剖分可使水面更平滑。 j#5a&Z
如果你是从别的小节转到此地,你可以跳过框中的步骤。 XIl<rN@-
+ File bkV<ZUW|;
+ Open 5E'/8xp bB
此时针对出现一个文件选择窗口。 xSmG,}3mF
++ Grid_1.Vmf (Tutorial目录下) |IcW7(
打开文件后,就会出现Visual MOFLOW的主菜单。 nHdQe
首先我们要细化供水井周围的网格。从主菜单的顶部菜单条上: T7,]^
1
+ Input 1&=)Bxg4
就进入了输入模块中,其中模型的组成部分(Building Blocks)在顶部以菜单方式给出。当你首次进入输入模块时,Visual MODFLOW会装入网格(Grid)输入窗口。屏幕左边的最上面的四个按钮(View Column、View Row、View Layer、Goto、Next和Previous)在每个屏幕上都会出现,使你在任何时候都可以将模型的显示方式从平面变为剖面。屏幕左边其它的按钮有各种进行模型网格修改的功能。要细化X方向的网格, @s LN
+ Edit Columns U6M~N0)Yr
然后将鼠标移至网格的任一处。注意在模型网格上有一个高亮度的竖线会随鼠标移动。这条竖线可以用来在模拟区的任何位置加上一列。在本练习中,我们通过按鼠标右键引出一个添加垂直网线(Add Vertical Line)窗口,以指定的间隔增加网格线。如图所示。 OX7=g$S 1
+ ¤Every spaced gridlines from: \)~d,M}kK
单击其后的输入框,输入以下数值: P51M?3&=l
from: 500Ö r N$0qo
to: 1600Ö ,>a!CnK=
at intervals of: 25Ö , _e[P
+ OK接受这些值 fvC,P#z'|
+ Close离开 Pao^>rj
现在我们要细化从加油场至供水井Y方向的网格。 J\@6YU[A
+ Edit Rows eP-|3$
然后将鼠标移至网格的任一处,按鼠标右键,将会有一个弹出窗口要求你输入网格的信息。 Njc@5*rJ&
+ ¤Every spaced gridlines from: z{rV|vQ
单击其后的输入框,输入以下数值: Dp([r
from: 400Ö DhyR
to: 1900Ö EK';\}
at intervals of: 25Ö +8Lbz^#
+ OK接受这些值 &u>dKf)5
+ Close离开 -l= 4{^pK
细化后的网格如图所示。 B)JMughq_
下一步将向你展示如何在剖面上观察模型,以及向模型添加层位。要在剖面上观察模型,在左边菜单上 FH,]'
+ View Column YRu@;
`
将鼠标移至网格的任一处,当你在屏幕上移动光标时,光标对应位置的列就会以高亮度显示。如要观察某一列,在该列上单击左键。Visual MODFLOW将模型的显示方式从平面变为剖面视图。此时模型还没有进行垂向放大,剖面图看起来就象一条粗线,根本不能区分出三个层来。为了能恰当地显示出这三个层,你必须对剖面图进行垂向放大。在窗口底部选择: LiG$M{ 0
+ F8-Vertical Exaggeration {18hzhs
就会出现一个垂向放大窗口,供你输入垂向放大系数。 e(0OZ_ w
Type: 40 Y9)j1~
+ OK i0=U6S:#
如下图所示,这模型的三个层就会在屏幕上显示出来。从图上可以看出在整个剖面上每一层的厚度都是一样的。 8o~
NJ 6
在Visual MODFLOW中,你可以从Surfer的.grd文件或以空格分隔的x, y, z的 ASCII文件中输入变化的层位标高。本例中,我们将输入一个ASCII格式的x, y, z文件来生成一个有坡度的地表及层位,同时后者的厚度也是变化的。 PQ&*(G
+ Import Surface BHY8G06
此时将出现如下窗口: wMW<lT=;
+ Choose Filename ) D@j6r
显示一个所有ASCII文件的选择窗口。你将输入的地表高程文件名为airpt_gs.asc。要选择此地表高程文件, lFzVd
N
++ airpt_gs.asc ~|+ ~/
应确认以下缺省设置。 sW]_Ky.]
+ OK QLd*f[n
这是一个有一定坡度的地表,北部边界高程为20米,南部边界高程为17米。Visual MODFLOW从ASCII文件中读取X,Y,Z数据,并利用五个最近采样数据点利用距离平方倒数插值方法求出第1层的每个单元顶板的块中心高程。 =
lo.LFV
现在要输入层1底板的高程。 'ITq\1z
+ Import Surface 6
VEB2F
+ Choose Filename sA2-3V<t8
++ airpt_b1.asc 第1层的底板高程数据文件。 rwP#Yj[BK+
+ ¤ Import bottom elevation of layer(输入第...层的底板高程):(默认为1) X)b$CG
+ Minimum Layer Thickness 0.75 ?
:H+j6+f
+ OK 接受其它的默认参数值 , PlH|
几秒钟以后,第一层底板就会显现出来,它具有一个变化的表面。下一步要输入第二层底板高程。从左侧菜单上选择: `_\KN_-%Vu
+ Import Surface gb_Y]U
+ Choose Filename(选择输入的文件) b>-DX
++ airpt_b2.asc 第2层的底板高程数据文件。 j?K$w`
+ ¤ Import bottom elevation of layer:向框中输入2。 QG5WsuT
+ Minimum Layer Thickness 0.75 "hk#pQ
+ OK 接受其它的默认参数值 `2.c=,S{
第二层的底板就会显现出来,它也具有一个变化的表面。接着输入第三层底板高程。从左侧菜单上选择: _uIS[%4g
+ Import Surface $.`o
+ Choose Filename(选择输入的文件) CRb*sfKDL
++ airpt_b3.asc 第3层的底板高程数据文件。 =5|7S&{
+ ¤ Import bottom elevation of layer:向框中输入3。 @DysM~I
+ Minimum Layer Thickness 0.75 c1f"z1Z
+ OK 接受其它的默认参数值 {-Y;!
第三层的底板就会显现出来,它具有一个变化的表面。 PT/Nz+
为对模型网格进行更细的垂向离散,你可以对每一层进行细分。从左边的菜单中选择: J K/{IkF
+ Edit Layers ]\pi!oa
把鼠标移到模型横剖面中。在横剖面中将会有一条水平方向的线随着鼠标移动。要细分某个层位,只要让这条线移到第一层大致的垂向中央并单击鼠标左键。在此处就会加入一条分层线。 :&BPKqKp
为了细分中部承压层,把这条水平线移到弱透水层中单击鼠标右键。将出现一个添加水平层位(Add Horizontal Layer)窗口,如下所示。 rW?WdEg
+ ¤ Split current layer into ‘2’ evenly spaced layers(将当前层分为2个等大小的层) ~q`!928Gu
+ OK U6.aoqb%
对模型的底层重复这一过程。 f?5>V
当你完成这些操作后: )?rq8VO
+ Close 退出 <h#*wy:o2
现在,模型的横剖面上应该有六个层位,类似下图: FV^jCseZ
用左边菜单的Next和Previous按钮在模型列剖面中移动。 \|HtE(uCM1
接着选择View Row,在剖面的一高亮垂直列中用鼠标选择一行。使光亮区在模型剖面的行上移动。 )9eIo&Nl
要返回到模型区的平面视图状态,选择View Layer按钮,然后使模型顶层变成高亮显示,然后单击鼠标左键。就会显示飞机场的平面图。 :kgh~mx5LF
2.5.4第4节:添加井位 CS[]T9|_
这一节的目的是一步一步地指导你怎样在模型中加入抽水井。 Q|_F
P:
如果你已经阅读完了其它小节,你可以跳过该框中的步骤。 :c*"Dx'D
+ File im>(^{{r&
+ Open zhn?;Fi
将出现一个文件对话窗口。 :60vbO
++ grid_2.vmf(细化网格后的模型) "Z@P&jl
此时文件就会打开,出现Visual MODFLOW主菜单窗口。要进入输入模块: * 6}M.`.-
+ Input R$&;
从顶部菜单中选择 B!RfPk1B<*
+ Well -`L`kL<
再从下拉菜单中选择: 3%POTAw%
+ Pumping Wells ktyplo#F
系统将要求你保存你的数据 <pX?x3-'
+ Yes BE?]P?r?
现在开始加入井位,将供水井周围区域(模拟区的右下部)放大: v-8>@s jy8
+ F5-Zoom In K6v~!iiK$
把光标移到供水井的左上方,单击鼠标左键。然后,拖出一个方框围住供水井,再次单击左键,闭合放大窗口。 lSX1|,B7:]
+ Add Well Q>1BOH1by
把光标移到西边那眼井上,单击鼠标左键。就会出现一个井位编辑窗口,你可在此输入井位资料。 45yP {+/-Q
在标有Well Name的输入框中单击鼠标: ;$D,w
Type: Supply Well 1 MzvhE0ab
+ Add Screen C*Q7@+&
把光标移到窗口左边的井孔中,在高程为5米处单击一下鼠标。然后把鼠标移到高程为0.3米处再单击一次。就会在井中高程大约为5-0.3米处加入一滤水段。另外你也可以在方框中直接输入数值来精确确定滤水段的位置。 Cv#aBH'N
因为Visual MODFLOW计算的是每个单元的水头值,而不是井的水头值,因此Visual MODFLOW允许你把井当作观测点。然后你可以输入抽水井的半径,Visual MODFLOW会去计算校正后的水头值。这在与抽水井的实测降深进行校验时非常有用。 s#$t!F??9
下面我们要输入井位抽水的时间安排,在标有Stop的列下的方框内单击鼠标左键,并输入以下信息: 0+\725DJ
Stop[day]: 7300Ö A]bb*a1
Rate[m3/d]: -400Ö Z`-$b~0
+ OK 1<!P:@(
如果某些所需的数据没有输入的话,Visual MODFLOW会提醒你马上完成。 +/B
下一步是设置第二眼抽水井的参数,我们要使用从一眼井向另一眼井复制属性的快捷方法来完成。 cK$yr)7
+ Copy Well A{_CU-,
把光标移到西边的那眼井上,单击鼠标。然后把光标移到东边井位上,再次单击一次,就把井复制完成。下一步我们要编辑复制过来的井的资料。 alJ0gc2?
+ Edit Well ,hzRqFg2
单击井位标记,选择新的那眼井。就会出现一个井位编辑窗口。 e:G~P
u`
在标有Well Name的输入框中单击鼠标: 24wDnDyh
Type: Supply Well 2 <#0i*PM_
在标有Rate[m3/d]的方框中单击鼠标: vQ<
~-E
Type: -550 UBd+,]"f
+ OK 接受这些井参数 S*l/
Sa@
+ F6-Zoom Out +yL; ?+s>=
2.5.5第5节:设置模型的属性 vP{i+s18B
本节将向你显示如何一步步地设置一个由水力传导系数相差很大的层位所组成的模型的属性。 2IYzc3Z{9
如果你已经完成了其它小节,你可以跳过框中的步骤。 9a_P 9s3w
+ File 3,p!Fun:r
+ Open ZA=J`->k
出现一个文件选择窗口。 ,i_+Z
|Ls
++ Tutorial文件夹 {\>4)TA
++ grid_2.vmf &N.pW=%,N
打开文件,出现Visual MODFLOW主菜单。 O>)eir7
+ Input KR.;X3S}
+ Properties 0H_Ai=G
+ Conductivity SS/vw%
+ Yes (保存你的文件) t|@5,J
如果你正在创建一个新的模型,在此将出现一个默认K值(Default K Property) 窗口,让你输入水力传导系数、储水系数、给水度以及空隙度的初值。模型会把这些值赋给模型中的每个单元。 o@A`AA9
我们可以使用Assign Single、Assign Polygon或Assign Windows来输入具有不同的水力传导系数的区域或层位。 WsbVO|C
用给定的格式输入如下值。 ?;tPqOs&
Hydraulic Conductivity in x (m/sec) Kx: 2e-4Ö xa`xHh{0
(X方向水力传导系数Kx) -'c
qepC{T
Hydraulic Conductivity in y (m/sec) Ky: 2e-4(输入Kx后自动出现)Ö RxP~%oADw
(Y方向水力传导系数Ky) her>L3G-E
Hydraulic Conductivity in z (m/sec) Kz: 2e-4Ö bqn(5)% {
(Z方向水力传导系数Kz) 45 biy(qa
Specific Storage (Ss): 1e-4Ö C).\ J !
(单位储水系数) HD<$0M|
Specific Yield (Sy): 0.2Ö ]v<8l4p;
(单位给水度) }8^qb5+!3
Effective porosity (Eff.Por): 0.15Ö &l cfX\y
(有效孔隙度) "e-z2G@z
Total Porosity (Total Porosity): 0.15 .5 E)dU
(总孔隙度) ltOS()[X
+ OK 接受这些数值。 *o1US
当你输入Kx的值时,Ky就自动给出了,因为Visual MODFLOW假设模型是水平各向同性的。 L\mF[Kd#+T
现在输入第三和第四弱透水层的水力传导系数值。 p7\LLJ y
+ Goto r QzdHA
显示转到层弹出窗口 ,OlS>>,
Type: 3 [w}- )&c
+ OK J>R$K
Assign Windows让你可以在一个矩形窗口内设置一个不同的水力传导系数。现在,让我们给第三、四层弱透水层输入一个较小的水力传导系数。 ET^?>YsA
+ Assign Windows POk5+^
把鼠标移到网格的西北角,在单元中间单击。然后把鼠标移到东南角,点击单元中部。这就形成了一个覆盖整个层面的窗口。 ?[#w*Am7
出现一个属性设置(Assign Property)对话框。 cPcH
8Vd
+ New (整个网格将变成蓝色) ,LZA\XC
输入下列弱透水层的水力传导系数: l]wfL;u
把鼠标移到Kx(m/s)框中单击并输入如下值: bF9.k
Kx(m/s) 1e-8Ö q=^;lWs4
Ky(m/s) 自动设置为1e-8Ö (?A
c`H
Kz(m/s) 1e-9Ö Wm1dFf.>
+ OK 接受这些值 e
oE)Mq
从底部菜单栏中选择F6-Zoom Out更新屏幕显示。 wHbmK
现在把第三层的水力传导系数复制到第四层。在左边菜单上选择: &wDZ@{h
+ Copy Layer bB0/FiY7o
出现一个复制窗口。在窗口中: >}wFePl
+ þ Copy All Properties(复制所有的值) wMkHx3XD
+ Layer 4(它就会以高亮度显示) I3
6@x`f
+ OK 把第三层的K值复制到第四层 b
B#QIXY/L
尽管在含水层天窗地段代表弱透水层的模型层已收缩到非常薄(0.2m),但是这对于重新设置K值来表现弱透水层的不连续性是一个很好的练习。 QcWg
在这个特例中,含水层天窗区在底图上就已有体现。然而,在许多情况下是没有这种资料的,你必需依靠某些参数比如层厚来确定这个不连续的地段。要显示第三层的厚度可在下部菜单条上选择F9-Overlay按钮。这时就会出现一个图层控制(Overlay Control)窗口(如下图),其中有一个表列,上面包含以字母顺序排列的所有在设置模型输入参数可能出现的图层。 ,$!F,c
向下浏览表列直到找厚度图层(Layer Thickness Overlay),然后双击它,在它的左边就会出现一个“*”,表示它现在是一个活动图层。 :Xy51p`.;]
+ OK 将显示层厚的等值线图 4wv0~T$;x
现在将等值线收缩到0.5m的最小厚度的地段放大。 ;8=Bee4
+ F5-Zoom In Tf=1p1!3
将鼠标移动到含水层天窗地段的左上部,单击左键。在该区上拉出一个框,再单击一次。屏幕上就会出现含水层天窗地段放大后的图形。 S?Cd,WxT
+ Assign Single(对单个单元进行设置) XRV~yBIS
将会出现一个属性设置框,其缺省值是最后输入的K值(Property #2)。按下New键下的向下箭头。Property #1 就会是K的当前值。此时千万别选择OK。将鼠标移到由0.5m等值线所圈定的表示含水层天窗的地段。按下鼠标,绕着该区拖动直到区内单元格成为白色显示(如图)。 <:AA R2=
如果你已对区外的单元重新作了设置,那么,松开左键,按下右键就可以将单元恢复为它们原先的值的颜色。完成单元格的“描画”后,选择Assign Single窗口中的OK键。 ?o)?N8U
现在你要将该值复制到第四层。 nKd'5f1
+ Copy Layer(在左侧菜单条上) fhn$~8[_A
就会出现一个复制窗口,其缺省设置是Copy only property #1(仅复制第一种属性) '-(Z.e~e
+ Layer 4 激亮该层 7
2i&-`&4
+ OK 把层3中的K的Property #1值复制到层4上 [.Fm-$M-
要检查K值是否已经复制到层4,在左边菜单上选择Next进入层4。 ){mqo%{SO
+ F6-Zoom Out 回到整个模拟区的全屏显示状态。 tVf):}<h
要在剖面上观察模型,从左边菜单选择View Column键,移到含水层天窗地段。通过左键选择要观察的列。就会出现一个如图所示的剖面图。 B4HMs$>
用左侧菜单上的Next和Previous按钮可在整个模型中移动。通过屏幕左下侧的立方体导航器可以观察剖面的相对位置。 9e.v[K~
选择View Layer,然后点中模型剖面中的第1层,回到平面显示方式。 zA?AX1%Wa
然后Visual MODFLOW将会计算模型顶层的层厚度等值线。 jHob{3
要去掉层厚等值线,可选择底部菜单的F9-Overlay。浏览图层表列,双击Layer Thickness Overlay,将它变为不活动即可。 ?*UWg[
+ OK 'h;qI&
2.5.6第6节:设置模型边界条件 -+i7T^@|
本节介绍了设置各种模型边界条件的步骤。 mb`h
如果你已经完成了其它小节,你可以跳过框中的步骤。 16N|
+ File cxQ %tL+S&
+ Open hVyeHbx
出现一个文件选择窗口。 Rt,po
++ 在Tutorial文件夹中 w_"-rGV
++ grid_2.vmf \PZ;y=]p}
打开文件,出现Visual MODFLOW主菜单。 pm9sI4S
+ Input +BM[@?"hrh
要设置补给条件,你必须是在处于模型顶层的平面视图情况下。检查一下左下角的立方体导航器,看看你目前所处的层位,如果不在第1层,就用Next,Previous或Goto进入到第1层。 \'zloBU
我们首先要设置的边界条件是含水层的补给通量。 ;y7+ Q
+ Boundaries 'ahz@+lO
+ Recharge 5{ !"}
如果你是正在建立新模型,此处会出现一个缺省补给(Default Recharge)窗口,提示你输入补给量的初始值。Visual MODFLOW自动地将这个补给量值赋给模型的整个顶层。在框中输入: #zgO_H
Stop Time[days]: 7300Ö -ert42fN
Recharge[mm/yr]: 100 PB*G#2W
+ OK接受赋值 7k%!D"6_R
+ F5-Zoom-In v.sjWF
将光标移到加油场的左上角,单击鼠标左键。沿加油场拉出一个框,然后再点一次闭合放大窗口。 9Y:.v@:}0
现在我们要给加油场设置一个较高的补给值: w87$p821
+ Assign Window X%og}Cfi
将光标移到加油场的某个角上。单击鼠标左键,沿加油场拖出一个窗口,然后再单击鼠标左键。就会出现Assign Recharge窗口。 cb|hIn\>7
+ New t@ri`?0w
设置了补给条件的窗口将变为Property #2属性的蓝色,给Property #2输入以下值: UCj<FN `
Stop Time [days]: 7300Ö dEI!r1~n
Recharge [mm/yr]: 250 M<=e~';H
+ OK接受赋值 hAds15 %C
+ F6-Zoom-Out将模拟区恢复为满屏显示。 V@y&n1?6
下一步是给承压和非承压含水层的南北边界设置定水头边界条件。 7$Wbf4
从顶部菜单上选择: ;*+H&
+ Boundaries OaT]2o
+ Constant Head Z;tWV%F5
+ Yes(保存输入的补给量数据) JSylQ201
首先我们为上部潜水含水层沿模拟区北部的边界设置的边界条件。 0k_3]Li=(
+ Assign Line YUTh*`1k<
将鼠标移至网格的西北角。用鼠标左键点击该单元中心。然后用鼠标右键点击该单元中心。此时单元格组成的水平线将以高亮度显示,同时出现一个如图所示的定水头边界输入窗口。 )p#L "r^)
输入以下数值: b&Laxki
Code #: 1 ePR9r}
+ Stop Time A42!%>PB
Stop Time: 7300Ö u|\?6fz
Start Point: 19Ö kaoiSL<[6
End Point: 19 p/l">d]+
+ OK接受输入 ?|Z~mE
此时,粉色线将会变成暗红色,表示常水头边界值已经设置完毕。从左边菜单上选择: '6g-]rE[
+ Copy Layer GjN6Af~}
会出现一个复制窗口,在Copy only code #(仅复制第...种属性)缺省设置旁边的框中值为1。 nW K7*
+ Layer 2(该层会高亮度显示) RFSwX*!
+ OK 将定水头边界1号属性复制到第2层 a3A3mBw
我们现在要输入下部承压含水层沿模拟区北部的边界的边界条件。 ~>_UTI
+ Goto(在边上的菜单条上) u"wWekB
此时会有一个Goto Layer的弹出窗口。在Layer you wish to go to(你想去的层位)框中,1是缺省值,并高亮度显示。 =w:)AWZ
Type: 5 48 0M|^
+ OK进入第5层 A@:U|)+4
+ Assign Line /2*BdE[yG
将鼠标移至网格的西北角。用鼠标左键单击该单元中心。然后用鼠标右键单击网格东北角的单元的中心。该线会以高亮度显示,同时出现一个定水头边界输入弹出窗口。输入以下数值: i6P$>8jBQ-
Code #: 2 mQSn*;9\T3
+ Stop Time pp{p4Z
Stop Time: 7300Ö M} ri>o
Start Point: 18Ö K3WaBcm
End Point: 18 G^.tAO5:f
+ OK接受输入 <R~;|&o,$
+ Copy Layer
}J-e:FUF#
会出现一个复制窗口,选择Copy only code #并在旁边框中输入“2”。 D#`>p
+ Layer 6(该层会高亮度显示) D
dCcsYm,
+ OK 将定水头边界2号属性复制到第6层 Z( "-7_
下一步我们要设置下部承压含水层沿模拟区南部的边界的边界条件。 =&Tuh}
+ Assign Line \MdieO*
将鼠标移至网格的西南角。用鼠标左键单击该单元中心。然后用鼠标右键单击东南角单元的中心。该线会以高亮度显示,同时出现一个定水头边界输入弹出窗口。输入以下数值: LCzeE7x
Code #: 3 |QR9#Iv
+ Stop Time 6hcK%0z
Stop Time: 7300Ö ,!PV0(F(
Start Point: 16.5Ö -II03 S1
End Point: 16.5 S,avvY.U\
+ OK接受输入 lOe|]pQ.,
+ Copy Layer 4M6[5RAW{
会出现一个复制窗口,选择Copy only code #并在旁边框中输入“3”。 vCK+v
r!
+ Layer 6(该层会高亮度显示) AGO"),
+ OK 将定水头边界3号属性复制到第6层 ^ZRZ0:rZ
设置完定水头值后 zKaj<Og
+ View Column T5I#7LN#
将鼠标移至模拟区,选择通过含水层天窗地段的列。将会出现一个与下图相似的剖面图。 0V^I.S/q
同时它还可以显示表示水力传导系数属性的颜色。选择F9-Overlay,将显示图层窗口, 此窗口显示了一系列能显示和关闭的图层。使用鼠标双击Conductivity Overlay。将在Conductivity Overlay旁边显示一星号,表明它被激活了。 B(} 'yY@%u
+ OK 显示水力传导系数层 X8R:9q_
现在从左边菜单栏选择View Layer并单击模型横剖面的顶层,返回显示模拟区的平面图。 Zu>-y#Bw
下一步沿着模拟区南部边界的顶层输入河流边界条件。 Z#^|h0
+ Boundaries JY;u<xl
+ River 切换到河流输入窗口 c^puz2
+ Assign Line -XK;B--c
按照图层布置,(从网格的西南部开始)沿着河流单击鼠标左键使河流数字化。当到达东南边界时单击鼠标右键。 BgDWl{pm
此时会出现一个弹出窗口将让你输入河流的资料。 cq3Z}Cp
键入以下数据: (ov&iNx
Code # 4 TiG?r$6v%
þ Assign to appropriate Layer M`@AS L:u
+ Stop Time 5~im.XfiVx
Stop Time 7300Ö j+q)
Start Point River Stage 16.0Ö oe3=QE
Start Point River Bottom 15.5Ö l%?4L/J)#
Conductance 1000 Ö J~oxqw}
End Point River Stage 15.5 Ö BSJS4+,E
End Point River Bottom 15.0Ö Aho*E9VW
Conductance 1000 ]%FP*YU4O
+ OK 接受这些值 /eXiWa sQ
+ F6-Zoom Out 刷新屏幕 #T{)y
定义河流后,将显示一条用网格单元表示河流边界条件的蓝线,如下图所示。 3c(mZ
MA}}w&
注意:技术角度上来说,这些传导系数值是错误的,因为在MODFLOW 中传导系数的定义如下: G]EI!-y
传导系数= 河床的水力传导系数 ´河流面积 SXO.|"M
河床厚度 ;_|4c7
细化网格时产生不同表面积的单元,相应地应输入不同的传导系数值。为了节约时间,在该例中不给单个单元赋予不同的传导系数值。 HD2C^V2@M
2.5.7第7节:设置质点 7sWe32
本节将一步一步地指导你输入向前追踪的质点以确定污染物的途径。 v<<ATs%w
如果你已经完成了其它小节,你可以跳过框中的步骤。 }; M@JMu,
+ File P>_9>k@;Q
+ Open ZJF+./vN
出现一个文件选择窗口。 F#RN m5
++ 在Tutorial文件夹中 '|), ?
++ grid_2.vmf mpVD;)?JmM
打开文件,出现Visual MODFLOW主菜单。 R?/xH=u>
+ Input N`3^:EJL8
在顶部菜单栏上选择: \&ZEIAe
+ Particles Doq}UWp
此时程序会要求你保存数据: #n15_cd
+ Yes B.;@i;7L
我们在此要做的第一步是放大加油场地段。 XzqB=iX
+ F5-Zoom In 在底部菜单栏 "USzk7=&.
在加油场的左上角单击鼠标左键,然后在该区上拖出一窗口,再单击鼠标左键闭合放大窗口。 KsK]y,^Z
+ Add Line jvD_{r
把光标移到加油场的左边并单击鼠标左键。拖出一根线到加油场的右边,然后再单击一次。将显示一添加质点窗口。缺省的质点数是10。把质点数变成5。 sDTw</@
+ OK 设置好加油场中的五个质点的线 Gm9
通过加油站区的绿色质点线表示是向前追踪的质点。现在返回到模拟区的满屏显示状态。 s,!vBSn8
+ F6-Zoom Out .3V L
2.5.8第8节:设置模型的运移参数 RP[`\
如果你有MT3D1.1, 1.5或1.86,请进入第8A节。 1^V.L+0s]
如果你购买的是MT3D96,请进入第8B节。 }|N88PN
2.5.8.1第8A节:MT3D1.1, 1.5或1.86 DHuvHK0#
本小节将一步一步地向你介绍如何输入运行MT3D所需的污染物运移参数。 "L1LL
iS
如果你已经完成了其它小节,你可以跳过框中的步骤。 05\0g9
+ File K :+q9;g
+ Open .3XiL=^~Qp
出现一个文件选择窗口。 5[r}'08b
++ 在Tutorial文件夹中 f#m@eb
++ grid_2.vmf vWrTB
打开文件,出现Visual MODFLOW主菜单。 }Qe(6'l_
+ Input ZWzr8oY)
在顶部菜单栏上选择: W/ WP }QM
+ MT3D VAUd^6Xdwx
此时程序会要求你保存数据。 1>Vq<z
+ Yes L;L_$hu)
如果你以前没有进入过MT3D窗口,程序就会要求你输入质量运移参数所用的质量和浓度的单位。 Z(tO]tQE
然后程序会要求你输入整个模拟区的背景浓度。 U:eX^LE7
Type: 0 NZADHO@0
+ OK 设置好背景浓度 ycki0&n3
下一步我们要为模型的每一层输入一个污染物延迟系数。 P?c V d2Y
+ Chem. reaction (顶部菜单栏上) /-m)
+ Options/Defaults (边部菜单栏上) YGs'[On8
此时屏幕上会将显示一个化学反应(Chemical Reaction)窗口,如图所示。 "lC>_A
在本例子中我们只使用吸附等温线选项。 2Q@Jp`#,4
+ þ Sorption isotherm J kAd3ls
我们要使放射性衰变或生物降解选项(Radioactive decay or biodegradation Option)无效。如果它前面是þ,请在方框中单击鼠标去掉ü。注意现在窗口右下角的两个小窗口变灰了,表明它们已变为无效。 c=bK_Z_
现在我们要选择吸附等温线的类型。 cXf/
+ ¤ Linear isotherm c&wiTvRV
注意因为选择了线性吸附等温线,窗口右上角的Second Sorption Constant(吸附方程第二常数项)窗口已经无效了。 u2$.EM/iae
下一步,你必须输入每层的土壤容重值。 #.fJ
M:"tG
在窗口左下部的标有Bulk density of porous medium (Kg/m3) (孔隙介质容重)的窗口中的第一个方框中双击鼠标。方框左边的数对应于层号。每层的土壤容重值都有设为2650 Kg/m3。 i O? f&u
现在我们将输入用于计算污染物延迟系数的Kd值。 )+8r$ i
用鼠标单击First Sorption Constant(吸附方程第一常数项)窗口中第一个方框。 L]"$dF
Type: 5.7e-5Ö Vkd_&z7
所有六层都输入该值。 5 $$Cav
+ OK 接受这些化学反应参数 k $8Zg*)
2.5.8.2第8B节:MT3D96 0kls/^ 0,
本小节将向你介绍如何设置运行MT3D96所必需的污染物运移参数。 JdUdl_Dz
如果你已经完成了其它小节,你可以跳过框中的步骤。 ,O-_Pv
+ File nR@mm
j
+ Open +F-Y^):
出现一个文件选择窗口。 rUyGTe(@h
++ 在Tutorial文件夹中 9mv0} I
++ grid_2.vmf B"#pvJN
打开文件,出现Visual MODFLOW主菜单。 a_{'I6a*,
+ Input <gH-`3J6
从底部工具条上选择: 1=sL[I 7<
+ F10-Main Menu F0zaA
此时程序会要求你保存数据。 )D'#>!Y
+ Yes qTwl\dcncC
+ Setup sKR%YK
"A
+ Numeric Engines KHC Fz
新模型的缺省数值引擎如下图所示: On!+7is'
重要提示: 0:4>rYBC
Visual MODFLOW 包括免费的MT3D1.10(EPA)和MT3D1.5(DoD)。然而,MT3D1.10没有MT3D1.86那么稳定,且这些版本都没有MT3D96的速度快或有弹性。本教程中的机场练习是用MT3D1.10和MT3D96都做过。对那些要求进行精确的污染物运移模拟的所有专业咨询业务都应使用MT3D96。 kRb JK
要改变某个数值引擎,单击该引擎旁的下拉箭头键,在下拉框中选择你想要的数值引擎即可。在我们这次模拟中,我们将使用MT3D96(如果你没有购买MT3D96,你可以使用具有MMOC技术的MT3D1.10或MT3D1.5)。 0c!^=(
+ Transport Numeric Engine <&O*'
<6C
+ MT3D96 by Papadopulos _Associates Inc. 4oryTckS
+ OK 接受数值引擎 L>E{~yh
+ Input设置 b2[U3)|oO
+ MT3D P_8!Gp
然后你可以输入整个模拟区的背景(缺省)浓度(如果你事先没有进入过MT3D菜单)。 s~B)xYmyB'
Type: 0 WJa7
+ OK 接受背景浓度 _Ye.29
下一步你应输入每层的污染物延迟系数。 %Ny1H/@Q1+
+ Chem. reaction `nEqw/I
+ Options/Defaults c~OPH
0,
将出现化学反应窗口,如图所示。 AS~!YR
在本例中我们只使用吸附等温线选项。 <
]+Mdy
+ þ Sorption isotherm &sBD0R(a
在方框中单击鼠标去掉ü,使放射性衰退或生物降解选项无效。 v.TgB)
现在你必须选择吸附等温线的类型。 /dvronG
+ ¤ Linear isotherm 4tUt"N
+ OK S*J\YcqSC
MT3D96允许你按单元输入化学反应参数。缺省时,MT3D96给模型的每一层都赋予模型开始建立所定义的值。 'o0o.&/=
+ Database &OXx\}>MW
你必须为每一种化学反应属性输入土壤容重。 !^0vi3I
在标有Bulk density (Kg/m3) 的窗口的第一个方框中双击鼠标。方框左边的数对应于层号。把每层的土壤容重都设为2650 Kg/m3 。 8)pL0bg
现在我们将输入用于计算污染物延迟系数的Kd值。 #^9k&t#!6
在第一项系数(1st coeff.)窗口下的第一个方框中单击鼠标。 ;XjXv'
Type: 5.7e-5 5P\A++22Y
给所有的层都赋予这些值。 Cd]A1<6s
+ OK 接受这些化学反应参数。 P,WQN[(+
2.5.9第9节:标志污染源 Df *<3G
加油场的污染源将以污染物下渗源的补给浓度(Recharge Concentration)的方式给出。 :(yut
从顶部菜单条上选择: C|w<mryx
+ Boundaries vZ$E
[EG}
+ Recharge Conc. 83'+q((<
检查一下左下角的立方体导航器,确信你已经处于第一层上。如果不在第1层,就用Next,Previous或Goto进入到第1层。 3ZyvX]@_
+ F5-Zoom In YPY,gR
同时在加油场上拉出一个放大窗口。 .;ofRx<
然后,给进入模型的加油场的天然补给量设置一个浓度: mQ`2c:Rn&7
+ Assign Polygon NW3qs`$-(
沿加油场周长移动鼠标,并同时用鼠标左键单击该区的角点,然后单击鼠标右键,该区就会变成阴影图像。此时就会出现一个补给通量浓度设置(Assign Recharge Flux Concentration)窗口,输入以下数值: =Bm|9A1
Code #: 1Ö -Q?c'e
Stop Time[days]: 7300Ö '!?t+L%gO
Recharge Conc.[mg/l]: 5000Ö Gx}`_[-
+ OK接受这些补给浓度值 HyKA+7}
+ F6-Zoom Out回到模拟区全屏显示状态 Yh@2m9
MT3D最后一个需要输入的参数是弥散资料。从上部菜单条上选择: {.KD#W
$5
+ Dispersion +(afO~9
如果你以前没有输入过弥散资料,程序就会让你输入纵向弥散度和横向弥散度比值。 cTRtMk%^
Long. Dispersion[m](纵向弥散度): 0.1Ö dw5"}-D
Horiz. to Long. Ratio(水平横向/纵向): 0.1Ö c$BH`" <*
Vert. to Long. Ratio(垂直横向/纵向): 0.00001Ö Ij =NcP
+ OK接受弥散资料 wpi$-i`
在这个例子中,模型中每层的弥散系数都是均匀的。 vZ57
S13
最后,我们要向模型中加入三个观测井以监测整个模拟区内不同位置的浓度。 JhTr{8{
+ Calibrate @V5i
此时程序会提示你保存数据。 jrz.n4Y`
+ Yes W(4$.uZ)
+ Add Obs. }gt)cOaY
将鼠标移至模拟区的任意一点,单击左键。将会有一个编辑观测点(Edit Observation Point)窗口让你输入以下信息: cjO%X
Observation Name: OW1Ö JXq!v:w6
X Location: 760Ö xJU]py~o
Y Location: 1667Ö Vky]In=
+ Time Observed Conc.[mg/l] -vY5h%7kf
500Ö 500¿ IhZn
1000Ö 1600¿ 9y5JV3
2000Ö 2700¿ " Y^9g/
7300Ö 3300 }BL7P-km
+ OK接受这些数值 zb"4_L@m2
+ Add Obs. yqL" YD
再将鼠标移至模拟区的另外任意一点,单击左键。将会出现一个Edit Observation Point窗口让你输入以下信息: T{m) = (q
Observation Name: OW2Ö +r;t]
X Location: 760Ö 8Lx1XbwK
Y Location: 1350Ö 5$Kj#9g-#
+ Time Observed Conc.mg/l -OU{99$aS
7300Ö 800¿ _-543B}
+ OK接受这些数值 'kY/=*=Q
+ Add Obs. ~tRGw^<9
再将鼠标移至模拟区的另外任意一点,单击左键。将会出现一个Edit Observation Point窗口让你输入以下信息: *l"CIG'
Observation Name: OW3Ö ^E8qI8s
X Location: 760Ö U] GD6q
Y Location: 900Ö (EosLn
h0
+ Time Observed Conc.mg/l oaPWeM+
7300Ö 100 La'6k
+ OK接受这些数值 |mHf7gCX
要进入下一模块,请在选底部菜单条上选择F10-Main Menu,再选择Yes保存数据。