-
GMS
使用说明
第一章
概述
该文档包含了:
GMS
及各组成模块的使用指南。
第二章
用
T
INS
建立地面模型
2.1
导入顶点文件
(
1
)
p>
点击
图标,然后选择“
FILE
”菜单下的“
OPEN
”菜单项,打开
“
”文件;
(
2
)
然后选择
“
Build TIN
”
下的
“
Trian
gulate
”
菜单项;
建立三角网格
,
并自动利用
Delauney
进行优
化,避免病态矩阵;
(
3
)
p>
选择
“
Display
”
菜单中的
“
Display O
ptions
”
命令;
选中
“
Contours
”
和
“
TIN Boundary
”
复选框,绘制地面等值线图;
(
4
)
p>
选择
“
Display
”
菜单中的
“
Display O
ptions
”
命令;
取消
“
Contours
”
和
“
TIN Boundary
”
复选框,
选中
“
Tr
iangles
”
复选框,
按
“
OK
”
;
然后选择
“
View
”<
/p>
菜单下的
“
Oblique
Vie
w
”命令,再选择“
Display
”菜单中的“<
/p>
Shade
”命令绘制带阴影的曲面图;
编辑顶点:
(
1
)选择“
Display
”菜单中
的“
Display Options
”命令;选中“
Contours
”和“
Vertices
p>
”复选框,
并单击“
Contours
p>
”右边的的“
Options
”按钮;更改
等值线间隔;
(
2
< br>)拖动顶点:从“
Tool
Palette
”上选择“
Select
Vertices
”工具
,选择“<
/p>
View
”菜单下的
“
< br>Plan
View
“命令,然
后用鼠标选中顶点拖动。选择“
View
”菜单下的“
Oblique Vie
w
< br>”命令,也可以
进行拖动。
(
3
)用键盘进行精确移动:选中顶点,在工具条中的
X
、
Y
、
Z
坐标文本框中输入精确数值。
锁定顶点:
按下
SHIFT
键,
选择要锁定的顶点,
然后选择
“
Modify TIN
”
菜单下的
“
Lock/Unlock
Vertices
”
命令。锁定了的顶点为另一种颜色。
添加顶点
:在“
Plan
View
”模式下选择“
Create
Vertex
”
(下面一组工具中的
)
,然后
用鼠标单击来添加顶点。
<
/p>
删除顶点:
择“
Edit
”菜单下的“
Delete
”命令。然后利用“
Create Triangle
”工具建立新的
三角形。
在
“
Modify TI
N
”
菜单下的
“
Vertex Options
”
命令,
在弹出的对话框中选择
“
Retriangulate
after deleting
”
,可
以在删除顶点后自动连接新的三角形。
平滑
TIN
文件:
选择“
Edit
”菜单下的“
Delete All
”命令删除所有
TIN
,然后打开比较稀疏的新的
TIN
文件,然后
选择
“
Build TIN
”
菜单下的
“
TIN ->
Scatter Points
”
命令,
在弹出的消息框中选择
NO
;
选择<
/p>
“
Modify
TIN
”
菜单下的
“
Uniform
ly Subdivide TIN
”
命令,
< br>并选择合适的系数;
然后单击
“
2D Scatter Point
”
,
选择“
I
nterpolation
”菜单下的“
to Active
TIN
”命令,输入新的文件名,选中“
Map elevat
ions
”
复选框,即可。可以通过选择“
Edit
”菜单下的“
Delete All
”命令删除所有
TIN
。
1
GMS
使用说明
在
p>
GMS
中可以同时建立几个
TIN
,通过点击选择
TIN
图标
可以选种
TIN
;双击
TIN
p>
,则在其上
出现一个
A
,说明该
TIN
激活;可以通过“
Display
”菜单下的“
Hide
”或“
Show
”命令隐藏或显
示
p>
TIN
。
第三章
立体地层模型
TIN
模块
Borehole
模块
Solid <
/p>
模块
Map
模块
建立立体模型需要三步:
1
、
读入钻孔数据:
可以根据参考手册说
明的文件格式输入钻孔数据;
或者通过
“
Boreholes
”
菜单下的
“<
/p>
Borehole
Editor
”命令。
读入数据文件:单击“
”进入
Borehole
p>
模块,选择“
File
”菜单下的“
Open
”命令,选择
“
tutorialsolids
”下的“
< br>”文件。
改变
Z
轴刻度:选择“
View
”菜单下的“
Z
Magnification
”命令,输入比例系数。
显示钻孔名称:
< br>选择
“
Display
”
菜单中的
“
Display Options
”
命令或按钮
复选框。
2
、创建扩展多边形:
(
1
p>
)选择“
Plan
View
”
图标,选择“
Zoom
”工具
;
p>
选中
“
Hole
Names
”
(
2
)选择“
Display
p>
”菜单下的“
Drawing Grid Options
”命令,将“
Spacing
”设置为
50.0
,选
择
“
p>
Snap to grid
”
和
“
Display grid lines
”
p>
复选框,
将
“
Gr
id line spacing increment
”
设置
为
1
,
“
OK
”
。
(
3
)
p>
定义边界:
扩展多边形用
Map
模块中的特征对象工具输入。
首先,
确定一个包围
研究区的
单个弧段:
进入
Map
模块
,选择
“
p>
Create Arc
”
工具,
然后用鼠标按顺时针或逆时针方向依
次单击边界
角点,绘制研究区的边界,在最后双击起始点,封闭区域。
(
4
)创建
多边形:选择“
Feature
Objects
”菜单下的“
Build
Polygons
”命令,
(
5
)选择
“
Display
”菜单下的“
Dra
wing Grid Options
”命令,取消“
Snap
to grid
”和“
Display
grid
lines
”复选框,将“
Grid line
spacing increment
”设置为
1
,“
OK
”。
(
6
p>
)选择“
View
”菜单下的“
Oblique
Vie
w
”命令。
3
、创建地表
TIN
3.1
选择接触点:
单击“
”进入
Borehole
模块,选择“
Select
Contacts
”工具
2
,在按下
SHIFT
的同时单击
GMS
使用说明
每个孔的顶部。
3.2
创建
TIN
选择“
Boreholes
”菜单下的“
Co
ntacts -> TIN
”命令,然后选择“
Auto-
extrapolate only
”选
项,
< br>OK
。
3.3
隐藏
TIN
转到
TIN
模式
模块
。
,选择
Select TINs
工具
,并单击来选择
TIN
,选择
Hide
命令
,返回
< br>Borehole
3.4
创建绿色层
TIN
单击“
”进入
Borehole
模块,选择“
Select Contacts
”工具
,单击<
/p>
3G
孔的绿
/
蓝
界面。
3.4.1
自动选择接触面
:选择“
Boreholes
”菜单下的“
Auto Select
”命令
3.4.2
创建
TIN
p>
选择
“
Boreholes
”
菜单下的
“
Contact
s -> TIN
”
命令,
然后选择<
/p>
“
Auto-extrapolate with
trimming
to selected boreholes
”选项,
OK
。;输入名字“
to
pgreen
”。
3.4.2
藏
TIN
转到
TIN
模式
模块
。
,选择
Select TINs
工具<
/p>
,并单击来选择
TIN
,选择
Hide
命令
,返回
Bo
rehole
3.5
创建红色层
TIN
单击“
”进入
Borehole
模块,选择“
Select Contacts
”工具
,单击<
/p>
4G
孔的红
/
蓝
界面;选
择
“
Boreholes
p>
”
菜单下的
“
Au
to Select
”
命令;
关闭
p>
“
Match material above
”
选项;
选择
“
< br>Boreholes
”
菜单下的
“
Contacts -> TIN
”
命令,
然后选择
“
Auto-
extrapolate only
”
选项,
< br>OK
;
输入名字
“
topred
”
。
转到
p>
TIN
模式
模块
。
,选择
Select TINs
p>
工具
,并单击来选择
TIN
,选择
Hide
命令
,返回<
/p>
Borehole
3.6
创建蓝色层
TIN
单击“
”进入
Borehole
模块,选择“
Select Contacts
”工具
,单击<
/p>
8G
孔的红
/
蓝
界面(上
红下蓝)
;
选择
“
Boreholes
”
菜
单下的
“
Auto Select
”<
/p>
命令;
选择
“
B
oreholes
”
菜单下的
“
Contacts
->
TIN
”命令,然后选择“
Auto-extrapolate
with trimming to selected boreholes
”选项,
OK
;输入
名字“
topblue2
”。转到
TIN
模式
,返回
Borehole
模块
p>
3.7
创建底层
TIN
。
,选择
Select TINs
工具<
/p>
,并单击来选择
TIN
,选择
Hide
命令
3
GMS
使用说明
单击“
”
进入
Borehole
模块,选择“
S
elect Contacts
”工具
,按下
< br>SHIFT
键,依次选择
8G
孔
的蓝
/
红下,
6G
的底面,
7G
的蓝
/
红,
5G
的蓝
/
红下,
2G
的底面;
选择
“
Boreholes
”
菜单下的
“
Contacts
->
TIN
”命令,然后选择“
Auto-extrapolate
with trimming to selected boreholes
”选项,
OK
;输入
名字“
botblue2
”。转到
TIN
模式
,返回
Borehole
模块
p>
。
,选择
Select TINs
工具<
/p>
,并单击来选择
TIN
,选择
Hide
命令
3.12
创建红色立体
3.12.1
创建立体
选择名字为
的
TIN
,选择“
Build TIN
”菜单下的“
Extrude
TIN -> Solid
”命令,输入
“
red
”,在
Elevation
中输入
-45.0
。
3.12.2
遮蔽立体
选择
Shade
命令
,创
建阴影;选择
Refresh
命令
,则
返回线框图。
3.13
创建兰色层
选择名称为
的
TIN
,(可以用“
Edit
”菜单下的“
Select from List
”命令来选择),按
下
SHIFT
键,
选择名称为
“
botblue2<
/p>
”
的
TIN
,<
/p>
选择
“
Build TIN
”
菜单下的
“
Fill
Between TINs -> Solid
”
命令,输入“
blue2
”,
OK
< br>。
3.14
去掉兰色层
此时,兰色层和红色层占
据相同的空间,要创建准确的地层模型,必须在红色层中的兰色层
位置做一个孔,这可用
在红色层中减去兰色层来实现。
进入“
Solid
”模式
,选择
Selec
t Solids
工具
,选择名称为“
red
”的立体,按下
SHIFT
键,
选择名为“
blue2
”
的立体,
p>
选择“
Solids
”菜单下的“
Set
Operations
”命令,在“
Operand #1
p>
”
框架中选则
Delete
选项,选择减号
命令
,则返回线框图。
,输入名称。选择
Shade
命令
,创建阴影;选择
Refresh
3.15
创建绿色立体
转到<
/p>
TIN
模式
,选择
Select TINs
工具
,选择名为“
< br>topgreen
”的
TIN
,
按下
SHIFT
键,选
择名为“
topred
”的
TINS
< br>,选择“
Build
TIN
”菜单下的“
Fill Between TINs
-> Solid
”命令,输入
“
gr
een
”,
OK
。
3.16
创建顶部兰色立体
p>
转到
TIN
模式
,
选择
Select TINs
工具
,选
择名为“
topblue1
”的
TIN
,按下
SHIFT
键,选
择名为
“
topgreen
”
的
TINS
及名为
< br>“
topred
”
的
TIN
;
选择
“
Build TIN
”
菜单下的
“
Fill Between TINs
-> So
lid
”命令,输入“
blue1
”,
OK
。
3.17
显示立体图
选择
Shade
命令
4
GMS
使用说明
3.18
切面
3.18.1
创建切面
进入
So
lid
模式
3.18.2
隐藏立体
选择
Select Solids
工具
View
”模式
。
,
从
“
Edit
”
菜单中选择
“
Select All
”
,
选择
Hide
命令
,
< br>进入
“
Oblique
,选择
Create Cross
Section
工具
,选择
Plan
View
,
3.18.3
删除所有钻孔
进入
< br>Boreholes
模式
3.18.4
删除多边形
进入
Map
模式
,从“
Edit
”菜单中选择“
DeleteAll
”。
p>
,从“
Edit
”菜单中选择“
DeleteAll
”。
3.18.5
显示立体剖面图
选择
Shade
命令
< br>
3.19
层边界
p>
绿和兰色层边界在垂向上有突变,而不是渐变。边界的明显尖灭可以通过在创建立体之前用<
/p>
TIN
的编辑工具改变
TIN
顶点坐标来实现。
3.20
删除立体和
TIN
从“<
/p>
File
”菜单中选择“
New
”命令,选择“
No
”。
第
4
章
2
维地质统计
在
GMS
中,可以用
2D
Scatter Point
模块进行
2
维插值。
4.3
输入离散点集
文件格式参见
GMS
File
Formats
文件中的
Tabular Scatter
Point File
格式。
进入
2D Scatter
Point
模块
文件。
4.4
改变显示模式
进入
Display
Options
,选中
Data
colors
,选择
Scatter point symbo
ls
左侧的按钮,选择一种
,
从
“
File
”
菜单下
选择
“
Import
”
命令,
选择
“
2
”
图形,
OK
,
OK
。颜色表示了污染物的浓度。
从“
Data
”菜单下选择“
Color Ramp
Options
”命令,选中
Show Color
Legend
选项,来显示图
例。
4.5
创建限定网格
从“
Scatter
Points
”菜单下选择“
Bounding Grid
p>
”命令,输入
X
(
60
)、
Y
(
40
)方向上的数
目,
OK
,
OK
。
4.6
选择插值方式
4.7
线性插值
5
GMS
使用说明
选择“
I
nterpolation
”菜单下的“
Interp.
Options command
”命令,选择“
Linear
”选项,
OK
。
选择“
Interpolation
”菜单下的“
to 2D Grid
”命令,输入数据名“<
/p>
Linear
”,
OK
< br>。
4.8
绘制等值线
进入“
2D Grid
”模式
掉“
Nodes
”,
O
K
。
4.9
映射高程
进入“
Oblique View
”模
式
为“
linear
”的数据集。
p>
4.10
阴影化网格
选择
Shade
命令
4.11
深谷插值
进入
2D
Scatter Point
模块
,选择“
Interpolation
”菜单下的“
Interp.
Options command
”命
。选择“
Data
”菜单下的“
Ma
p Elevations
”命令,选择名称
,进入
Display Options
,选择“
Con
tours
”和
“
< br>Fringes
”,关
令,选择“
Clough-Tocher
”选项,
OK
< br>。
选择
“
Interpolation
”
菜单下的
< br>“
to 2D Grid
”
命令
,
输入数据名
“
ct
< br>”
,
选中
“
Map Elevations
”
(可以将
Z
方向的数值直自动赋上),
OK
。
选择
Shade
命令
4.12
简单
p>
IDW
插值
选择<
/p>
“
Interpolation
”
菜单下的
“
Interp. Options
command
”
命令,
选择
“
Inverse distance weighted
”
选项,
选择
“
Inverse distance weighted
”
选项右侧的
“
Options
”
的按钮,
在此对话框顶部的
“
Nodal
function
”
框架中选择
“
Constant
< br>”
,
在
“
Computation of interpolation weights
”<
/p>
框架中选择
“
Use all
points
”,
OK
,
CLOSE
。
选择“
Interpolation
”菜单下的“
to 2D Grid
”命令,输入数据名“
i
dw_const
”,
OK
。
选择
Shade
命令
4.14
带梯度面的
IDW
插值
选择
“<
/p>
Interpolation
”
菜单下的
“
Interp. Options command
”
命令,
选择
“
Inverse distance weighted
”
选项,
选择
“
Inver
se distance weighted
”
选项右侧的
p>
“
Options
”
的按钮,
在此对话框顶部的
“
Nod
al
function
”框架中选择“
Quadratic
”,
OK
,
p>
CLOSE
。
选
择“
Interpolation
”菜单下的“
to 2D Grid
”命令,输入数据名“
idw_
quad
”,
OK
。
选择
Shade
命令
4.15
截切
p>
由于插值的问题,出现复数,这是不合实际的。可以利用截切的方法进行修正。
选择“
Interpolation
”菜单下的“
Interp. Options
command
”命令,选择“
Truncate value
s
”选
项,选择“
Truncate
to specified range
”选项,输入最小和最大值,
< br>OK
。
选择“
Interpolation
”菜单下的“
to 2D
Grid
”命令,输入数据名“
idw_quad_
trunc
”,
OK
。
6
GMS
使用说明
选择
Sh
ade
命令
4.16
普通邻近插值
利用邻近点的拓扑关系进行插值,当离散点聚集时效果好。
p>
选择“
Interpolation
”菜单
下的“
Interp. Options
command
”命令,选择“
Natural neighb
or
”
选项,选则右侧的“
Optio
ns
”的按钮,在此对话框顶部的“
Bounding
pseudo-points
”框架中选择
“
< br>Quadratic
”,在“
Bounding
pseudo-
points
”框架中取消“
Extrapolate
beyond convex hull
”选项
(保证区外的浓
度为
0
),
OK
,
CLOSE
。
< br>选择“
Interpolation
”菜单下的“
to 2D Grid
”命令,输入数据名“
nn
”,
OK
。
选择
Shade
命令
4.17 Kriging
插值
p>
选择“
Interpolation
”菜单
下的“
Interp. Options command
”命
令,选择“
Kriging
”选项,
选
则其右侧的“
Options
”的按钮。
4.17.1
创建实验方差图
选择
“
Edit Variogram
s
”
按钮,
在
“
Variogram Editor
”
对话框右上部的
“
Experimental
variogram
editor
”区的“
< br>New
”按钮,
OK
。
4.17.2
创建模拟方差图
创建实验方差图后,则需要创建模拟方差图。
在“
Variogram
Editor
”对话框右下部的“
Nested
structure specification
”区的“
N
ew
”按钮,
在对话框左边名为“
Mo
del functions
”区中选择“
Gaussian<
/p>
”选项,在名为“
Model parameters
”
区中,为“
Contribution
”输入
1965.0
,为“
Range
”输入
63.0
,
p>
OK
,
CLOSE
,
CLOSE
。
选择“
Interpolation
”菜单下的“
to 2D Grid
”命令,输入数据名“
k
rig
”,
OK
。
选择
Shade
命令
4.18
转换数据集
进入“
2D
Grid
”模式
,在“
Data Se
t
”组合框中选择“
linear
”项
(等值线的网格变了,但是
高程值仍然是前一数据集中的数据,可以通过“
Map
Elevations
”命令修改)。选择“
Plan
View
”。
用颜色填充等值线:
进入
Display Options
,关掉“
Cells
”,选择“
Gri
d boundary
”,选择“
Contours
”右侧的
“
Options
”按钮,在对话框的左下角的“
Contour
method
”区中选择“
Color fill
between
contours
”
选项,
OK
。
4.19
利用数据计算器
有时利用数据计算器(
Data
Calculator
)比较插值所得的两组数据。
从“
Data
”菜单中选择“
< br>Data Calculator
”命令,选择数据集,输入表达式:
<
/p>
在
“
Expression
”
字段中输入
“
abs(h
-i)
”
,
在
“
Result
”
字段中输入
“
Difference
”
,
选择
“
Compute
”
按钮。
7
GMS
使用说明
选择“
Data Set
Info
”按钮,查看新的数据集的
统计信息。
OK
,
DONE
。
4.20
删除所有数据
从“
< br>File
”菜单中选择“
New
”命令
,
选择“
NO
< br>”。
第
5
章
3
维地质统计
5.2
所需模块和界面
3D Scatter Point
3D Grid
5.3
输入离散点集
文件格式参见
GMS
File
Formats
文件中的
Tabular Scatter
Point File
格式。
进入
3D Scatter Point
模块
文件,选择
Oblique View
< br>模式
,
从
“
File
”
菜单下选择
“
Import
”
命令,
选择
“
3
”
。
5.4
显示数据颜色
选择<
/p>
“
Display
”
菜单下的
“
Display Options
”
命令,
选中
“
Data colors
”
选项,
OK
。
选择
“
Data
”
菜单下的“
Col
or Ramp Options
”命令,选中“
Show
color legend
”选项,
OK
。
5.5
Z
轴放大
从“
Vi
ew
”菜单中选择“
Z
Magnification
”命令,输入放大倍数。
5.6
创建限界网格
从“
Scatter
Points
”菜单下选择“
Bounding Grid
p>
”命令,输入
X
(
30
)、
Y
(
10
)、
Z
(
10
)
方向上的数目,并确认选中“
M
esh Centered
”。
OK
。
5.7
简单
IDW
插值
选择<
/p>
“
Interpolation
”
菜单下的
“
Interp. Options
command
”
命令,
选择
“
Inverse distance weighted
”
选项,
选择
“
Inverse distance weighted
”
选项右侧的
“
Options
”
的按钮,
在此对话框顶部的
“
Nodal
function
”
框架中选择
“
Constant
< br>”
,
在
“
Computation of interpolation weights
”<
/p>
框架中选择
“
Use subset
of points
”
,
在
“
Computati
on of interpolation weights
”
区选择
“
Subset
”
按钮,
选择
“
Use
nearest ___ points
”
,
< br>并为点数输入
64
。
OK
,
CLOSE
。
<
/p>
选择“
Interpolation
”菜
单下的“
to 3D Grid
”命令,输入数据名“
idw_const
”,
OK
。
5.8
显示等值面
进入
“
3D Grid
”
模式
,
选择
“
Display Options
”
模式
,
关掉
“
Cells
”
选项,
选中
p>
“
Grid
s
hell
”、“
Fringes
”和“
Isosurfaces
”选项;选择“
Isosurfaces
”右边的“
Options
”按钮,在第一
个编辑区中输入
3000.0<
/p>
,选中值
1
和值
2
之间的“
Cap
”选项。
OK
,
OK
。
选择“
Display
”菜单下的“
Shading
Options
”命令,选择“
Hidden
surface
”选项,选择“
Use
light
source
”“
Smooth
features
”和“
Overlay edges
”,
OK
。
选择
Shade
命令
5.9
删除内部边
p>
从
“
Data
”<
/p>
菜单中选择
“
Iso-Surface
Options
”
命令,
在
对话框底部选中
“
Interior
edge removal
”
。
p>
选择
Shade
命令
8
GMS
使用说明
5.10
边界指定范围
改变显示选项以更清楚的显示颜色差异。
选择“
Data
”菜单下的“
Fr
inge Options
”命令,选择“
Color
fringe specified range
”选项,
在“
Minimum fringe value
”中输入
3000
,在“
Maximum fringe
value
”中输入
9000
。
选择
Shade
命
令
5.11
用
Z
轴选项
进入
3D Scatter
Point
模块
,
选择“
Interpolation
”菜单下的“
Interp. Options
command
”
命令,将“
z scale
”的值改为
8
。
p>
选择“
Interpolation
”菜单
下的“
to 3D Grid
”命令,输入数据名“
idw_const2
”,
OK
。
选择
Shade
命令
5.12
带梯度面的
IDW
插值
选择
“<
/p>
Interpolation
”
菜单下的
“
Interp. Options command
”
命令,
选择
“
Inverse distance weighted
”
选项,
选择
“
Inver
se distance weighted
”
选项右侧的
p>
“
Options
”
的按钮,
在此对话框顶部的
“
Nod
al
function
”框架中选择“
Gradient plane
”,
OK
,
CLOSE
。
选择“
Interpolation
”菜单下的“
p>
to 3D Grid
”命令,输入数据名“
idw_gp
”,
OK
。
选择
Shade
命令<
/p>
5.13
带
2
次函数的
IDW
插值
选择
“
In
terpolation
”
菜单下的
“
Interp. Options command
”
命令,
选择
“
Inve
rse distance weighted
”
选项,
选择
“
Inverse distance
weighted
”
选项右侧的
“
p>
Options
”
的按钮,
在此对话框顶部的
“
Nodal
function
”
框架中选择
“
Quadratic
”
,
在
“
Computation of nodal
function coefficients
”
区中选择
p>
“
Use
all points
”,
OK
,
CLOSE
。
选择“
I
nterpolation
”菜单下的“
to 3D Grid
”命令,输入数据名“
idw_quad
”,
OK
。
选择
Shade
命令
5.14
其他插值方式
5.15
显示污染羽的截面
进入“
3D
Grid
”模式
,选择
Side
View
命令
,选择
Create
Cross Section
工具
,
切一截面,然后进入
Oblique
View
模式
,选择
Display
Options
,关掉“
Iso-surfaces
”,
OK
。
选择
“
Data
”
菜单下
的
“
Cross Section Options
”
命令,
选择
“
Interior edge removal
”
和
“
Fringes
”
,
OK
。选择“
Data<
/p>
”菜单下的“
Fringe
Options
”命令,取消“
Color fringe
specified range
”,
OK
。
选择
Shade
命令
4.16
截切
Truncation
由于插值
的问题,出现复数,这是不合实际的。可以利用截切的方法进行修正。
9
GMS
使用说明
进入
3D
Scatter Point
,选择“
Interpolation
”菜单下的“
Interp.
Options command
”命令,
选择“
Truncate
values
”选项,选择“
Truncate to
min/max of active data set
”选项,
OK
。
选择“
< br>Interpolation
”菜单下的“
to 3D
Grid
”命令,输入数据名“
idw_quad_trunc
”,
OK
。
选择
Shade
命令
5.17
建立移动切面影像循环
用
Film
Loop
建立污染羽的动画显示。
5.17.1
显示模式
进入“
3D
Grid
”模式
,选择
Select
Cross Sections
工具
,选择切面,选择“
Edit
”菜
单下的
“
Delete
”
命令,
OK
。
选择
“
Data
”
菜单下的
“
p>
Color Ramp Options
”
命令,
关掉
“
Show color
legend
”选项,
OK
。选择“
Data
”菜单下的“
< br>Fringe Options
”命令,选择“
Color fringe
specified
range
”选项,
为“
Minimum fringe value
”输入
1000
,为“
Ma
ximum fringe value
”输入
15000.0
。
OK
。
5.17.2
设置影像循环
p>
选择“
Data
”菜单下的“
Film Loop
”命令,选择“
Setup
p>
”按钮,选择“
Animate cutting plane
p>
”
选项,选中“
Z
”选项,
OK
。
5.17.3
回放影像循环
p>
选择“
Play
”
按钮,选择
Stop
停止动画。
5.18
建立移动等值线影像循环
在
“
Film Loop
”
对话框中选择
“
Setu
p
”
按钮,
关掉
“
Animate cutting plane
”
选项,
选中
“
Anim
ate
iso-
surface
”选项,为“
Begin value
”输入
1000
,为“
End value
”输入
15000
。选择“
Cap above
”
选项,
选择“
Display values
”选项。
OK
。
选择“
Play
”
按钮,选择
Stop
停止动画。
选择“
Film Loop
”对话框中
的“
Done
”按钮退出。
第
6
章
MODFLOW
—网格化
p>
在
GMS
中,有两种方法可以构建
MODFLOW
模拟:网格化和概念模型法。前者直接利用
3D
网格,在单个单元上使用源/汇项和其他模型参数;后者利用
GMS
中
Map
模块中的
GIS
工具建立
概念模型,模型中的数
据拷贝到网格中。
实例:面积,
75000 ft ×
75000 ft
;参数见图。
6.4
创建网格
进入
“
3D Grid
”
模式
,
选择
“
Grid
”
菜单中的
“
Create Grid
”
命令,
在名为
“
X-dime
nsion
”
区域中的“
Length
”中输入
75000
,在“
Number cells
”中输入
15
;在名为“
Y-dimension
”区域中<
/p>
的
“
Length
”
中输入
75000
,
在
“
Number cells
”
中输入
15
;
< br>在名为
“
Z-dimension
”
区域中的
“
Length
”
中输入
15000
,
在“
Number cells
”中输入
3
;
OK
。
6.5 MODFLOW
模拟的初始化
p>
从“
MODFLOW
”菜单中选择“
New Simulation
”命令。
10
GMS
使用说明
6.6
基本模块
从“
MO
DFLOW
”菜单中选择“
Basic
Package
”命令。
6.6.1
标题
首先输入标题。在头文件的第一行
中输入:
–
Grid Approac
h
;在
第二行中输入你的名字和日期。
6.6.2
程序包
单击“
Packages
”按钮,在对话框中选择“<
/p>
Drain
”、“
Well
”、“
Recharge
”模块,在“
Solver
”
区中选择“
Strongly Implicit Procedure
”,
OK
。
6.6.3
单位
单击“
U
nits
”按钮,为长度单位输入“
ft
”,时间单位输入“
d
”,
OK
p>
。
6.6.4
IBOUND
数组
单击“
IBOUND
”按钮。
< br>
IBOUND>0
:有效单元
IBOUND=0
:无效单元
IBOUND<0
:定水头单元
6.6.5
初始水头
单击“
Starting
Heads
”按钮。
6.6.6
退出对话框
CLOSE
。
6.7
直接将
IBOUND
值赋给单元
利用“
Cell Attributes
”命令可以直接将
IBOUND
值赋给单元,首先选中单元(
最左侧两层的
单元)
6.7.1
显示最左列
选择
View the J Axis
6.7.2
选择单元
选择“
Select
Cells
”
工具,拖拉鼠标选择。
,
6.7.3
改变
IBOUND
值
p>
从
“
MODFLOW
”
菜单中选择
“
Cell Attr
ibutes
”
命令,
为
“
IBOUND
”
输入
p>
-1
;
OK
。
p>
在
“
Graphics
11
GMS
使用说明
Window
”外单击,取消单元选择。
选择“
View the K
Axis
”命令
,
6.7.4
检查值
p>
从“
MODFLOW
”菜单中选择“
Basic Package
”命令。单击“
IBOUND
”按钮。
6.8
BCF
模块
从“
MODFLOW
”菜单中选择“
BCF
Package
”命令。
6.8.1
层的类型
6.8.2
层的参数
6.8.3
顶层
选择“
Top
Elevation
”按钮,选择“
Constant ->
Layer
”按钮,输入
200
,
p>
OK
,
OK
。
p>
然后以同种方法输入:
Bottom
Elevation
(
-150
ft
)
、
Horizontal
Hydraulic Conductivity
(
50
ft/d
)
、
Vertical
Hydraulic Conductivity
(
5
ft/d
)。
6.8.4
中间层
然后以同种方法输入:
Bottom
Elevation
(
-400
ft
)
、
Horizontal
Hydraulic Conductivity
(
3
ft/d
)
、
Vertical
Hydraulic Conductivity
(
0.5
ft/d
)。
6.8.5
底层
然后以同种方法输入:
Bottom
Elevation
(
-700
ft
)
、
Horizontal
Hydraulic Conductivity
(
7
ft/d
)
、
Vertical
Hydraulic Conductivity
(
1.5
ft/d
)。
Close
。
6.9
Recharge
模块
从
“
MODFLOW
”
p>
菜单中的
“
Source/Sink
”
子菜单中选择
“
R
echarge Package
”
命令,
选择
“
Constant
->
Array
”按钮,输入值
0.003
,
OK
。
OK
。
6.10
Drain
模块
6.10.1
选择单元
选择“
Select
Cells
”
工具,安下
SHIFT
键,用鼠标选择。
6.10.2
指定
Drain
p>
从“
MODFLOW
”菜单中的“
Point Sources/Sinks
”命令,在“
Drain
”区选择“
Add
”选项,
选择在“
conductance
”中输入
80000
,
OK
。
6.10.3
指定
Drain
高程
< br>
从
“
MODFLOW
”
菜单中的
“
Sour
ce/Sink
”
子菜单中选择
“
p>
Drain Package
”
命令,<
/p>
在
“
elev.
”
栏中输入高程值,
OK
。
6.11
井模块
6.11.1
顶层井
选择
“
Select
Cells
”
工具,
安下
SHIFT
键,
用鼠标选
择。
从
“
MODFLOW
”
菜单中的
“
Point
Sources/Sinks
”命令,在“
Well
”区选择“
Add
”选项
,选择在“
Flow
”中输入
-432
000
,
OK
。
6.11.2
中间层井
单击左下侧的“
”进入下一层。
12
GMS
使用说明
选择
“
Select
Cells
”
工具,
安下
SHIFT
键,
用鼠标选
择。
从
“
MODFLOW
”
菜单中的
“
Point
Sources/Sinks
”命令,在“
Well
”区选择“
Add
”选项
,选择在“
Flow
”中输入
-432
000
,
OK
。
6.11.3
底层井
选择
“
Select
Cells
”
工具,
安下
SHIFT
键,
用鼠标选
择。
从
“
MODFLOW
”
菜单中的
“
Point
Sources/Sinks
”命令,在“
Well
”区选择“
Add
”选项
,选择在“
Flow
”中输入
-432
000
,
OK
。
6.12
保存模拟
从“
MODFLOW
”菜单中的“
Save As
”命令,输入文件名。
6.13
运行
MODFLOW
p>
从“
MODFLOW
”菜单中的“
Run MODFLOW
”命令,
OK
。
6.14
查看结果
从“
MODFLOW
”菜单中的“
Read Solution<
/p>
”命令,
OK
。
6.14.1
改变层
单击左下侧的“
”进入下一层。
单击左下侧的“
”进入上一层。
6.14.2
颜色填充等值线
p>
从“
Data
”菜单中选择“
Contour
Options
”命令,在对话框左下角许选择“
Color
fill between
contours
”。
6.14.3
显示颜色图例
从
“
Data
”
p>
菜单中选择
“
Color Ramp Op
tions
”
命令,
在对话框左下角许
选择
“
Show color
legend
”
。
第
7
章
利用概念模型建立
MODFLOW
7.1
问题描述
7.2
启动
GMS
7.3
所需模块和界面
3D Grid
模块
2D Scatter
Point
模块
Map
模块
MODFLOW
界面
7.4
导入背景图
7.4.1
读入图像文件
p>
选择“
File
”菜单下的“
Open
”命令,打开文件“
tutorialmo
dfmap
”,单击
Open
按
p>
钮。
7.4.2 IMAGE
文件
VSTIFF
文件
装入
TIFF
文件则要进行注册,也就
是选择
IMAGE
上的三个点并且输入其相应的实际坐标。
p>
7.5
特征目标(对象)
点、弧段、结点(弧段两端的点)、顶点(弧段中间的点)、面
7.6
创建局部源/汇图层(
Cov
erage
)
进入<
/p>
Map
模块
,
选
择
“
Feature Objects
”
菜单下的
“
Coverages
p>
”
命令,
将
“
p>
Default elevation
”
改
为
700
,将默认名字改为“
Sour
ces/Sinks
”,将“
Coverage type
p>
”改为“
MODF/
MT3D
/MODP
”,
选择“
Op
tions
”按钮,检查图层类型是否是
MODF/
MT3D
local sources/sinks
,
OK
。
7.6.1
定义单位
13
GMS
使用说明
单击“
U
nits
”按钮,为长度单位输入“
ft
”,时间单位输入“
d
”,
OK
p>
。
7.6.2
定义边界
选择“
Create Arc
”工具<
/p>
,沿边界造弧,(退一点时按
Backspace
键,取消造弧时按
ESC
键),最后在始点双击。
p>
7.6.3
拷贝边界
首先为每一层创建水力传导度图层:
选择“
Feature Objects
”菜单下的“
Coverages
”命令,选择“
Copy
”命令,将新图层名改为
“
Layer1
”
,
选择
p>
“
Options
”
按钮,
将
“
Coverage ty
pe
”
改为
“
MODF/
MT3D
/MODP layer attribu
tes
”
,
OK
。当图层
Layer 1
coverage
高亮度显示时,再次按“
Copy
”按钮,将新图层名
改为“
Layer2
”,
选择“
Options
”按钮,将“
Coverag
e type
”改为“
MODF/
MT
3D
/MODP layer attributes
”,在中
间对话框中,将“
Assign from
layer
”和“
to laye
”均
改为
2
,
OK
。
然后创建补给多边形图层:
选择“
Feature Objects
”菜单下的“
Coverages
”命令,选择“
Copy
”命令,将新图层名改为
“
Recharge
”
,
选择
“
Options
”
< br>按钮,
将
“
Coverage
type
”
改为
“
MODF/
MT3D
/MODP areal attri
butes
”
,
OK
< br>。
最后:
< br>选择“
Sources/Sinks
”图层,单击
OK
,退出
Coverages
对话框。
7.6.4
定义指定水头弧段
首先选择
“
Select
Vertices
”
工具
<-> No
d
”
命令;
选择
“
Select Arcs
”
工具<
/p>
,
选择两个顶点,
选择
< br>“
Feature Objects
”
< br>菜单下的
“
Vertices
,
选择弧段,
选择
“
< br>Feature Objects
”
菜单下的
“
Attributes
”
,
选择“
Specified Head
”,在底部对话框的第二个文本框内输入
2
,单击
OK
。
然后指定水头
值(为弧段两端点指定水头值,假设其间的水头线性变化):
选择“
Select Points/Nodes
”工具
入其他
2
个交点的水
头值。
7.6.5
定义排水沟弧段
选择“
Create
Arc
”工具
键),最后在终点双击。
选择“
Select Arcs
”工具
,选择弧段,选择“
Feature Objects
”菜单下的“
Attributes
”,
,沿排水沟造弧,
(退一点时按
Backsp
ace
键,取消造弧时按
ESC
,双击
南部边界左侧的端点,输入值(
697.0
),
OK
。输
选择“
Drain<
/p>
”,在文本框内输入
6000.0
,单击
OK
。
然后
指定排水沟的高程值(为弧段两端点指定高程值,假设其间的水头线性变化):
选择“
Select Points/Nodes
”工具
,双击弧段的端点,输入值,
OK
。
7.6.6
创建多边形
先选择
弧段,
然后选择
“
Feature O
bjects
”
菜单下的
“
Build Polygons
”
命令。
然后选择
“
Feature
< br>Objects
”菜单下的“
Display
Options
”命令,选择
Polygons
(Fill)
选项,可以显示区。
14
GMS
使用说明
7.6.7
创建井
选择
“
Create Point
p>
”
工具
,
在井位处
单击,
要准确定位则可以输入坐标值;
选择
“
Feature
Objects
”菜单下的“
Attributes
”,选择“
Well
”,在文本框内输入开采量,并确定开采层位,单
< br>击
OK
。
网格细化:
选择
“
Select Points/
Nodes
”
工具
,
< br>选择所有井,
选择
“
Featu
re Objects
”
菜单下的
“<
/p>
Attributes
”
,
单击“
Refine
Point
”按钮,选择“
Refine grid in X
direction
”选项,为
Base cell size
输入
75
,为
Max
cell size
输入
50
0
;选择“
Refine grid in Y
direction
”选项,为
Base cell size
输入
75
,为
Max cell size
输入
500
;
OK
。
7.7
刻画补给区
7.7.1
选择
Recharge<
/p>
图层
从
Coverage
组合框中选择
Recharge
图层
< br>
7.7.2
创建垃圾填埋场边界
选择“
Create
Arc
”工具
,选择“
Select
Points/Nodes
”工具
,依次选择
< br>4
个端点。
,选择端点,输入
准确的
坐标值,选择“
Select
Vertices
”工具
7.7.3
创建多边形
先选择
弧段,
然后选择
“
Feature O
bjects
”
菜单下的
“
Build Polygons
”
命令。
然后选择
“
Feature
< br>Objects
”菜单下的“
Display
Options
”命令,选择
Polygons
(Fill)
选项,可以显示区。
7.7.4
指定补给值
选择
Select Polygons
工具
,双击补给区,输入补给强度;双击外围区,输入外围区的补
给强度。
7.8
定义水力传导系数
7.8.1
顶层
从
Coverage
组合框中选择
Layer
1
图层。后选择“
Feature
Objects
”菜单下的“
Build Polygons<
/p>
”
命令。双击区,输入水平和垂直渗透系数。
7.8.2
底层
从
p>
Coverage
组合框中选择
Layer
2
图层。后选择“
Feature
Objects
”菜单下的“
Build Polygons<
/p>
”
命令。双击区,输入水平和垂直渗透系数。
从
Coverage
组合框中选择“<
/p>
Sources/Sinks
”图层。
7.9
装入网格边框
选择“
Feature
Objects
”菜单下的“
Grid
Frame
”命令,选择“
New Frame
”按钮,然后用拖
拉或输入数值确定边界。
7.10
创建网格
选择“
Feature
Objects
”菜单下的“
Map -> 3D
Grid
”命令,输入数值。
7.11
指定有效单元
选择“
Feature
15
GMS
使用说明
Objects
< br>”菜单下的“
Activate Cells in
Coverage(s)
”命令。
7.12
初始化
MODFLOW
p>
的值
进入
3D
Grid
模式
7.13
转换概念模型
转到
Ma
p
模块
,选择“
Feature
Objects
”菜单下的“
Map -> MODFLOW<
/p>
”命令;确定选中
,选择“
MODFLO
W
”菜单下的“
New
Simulation
”命令。
“
All applicable coverages
”选项,
OK
。
7.14
层的高程插值
7.14.1
导入地面高程离散数据文件
转到
2D Scatter Point
模式
,
选择
“
File
”
菜单下的
“
Import
”
命令,
打开
“
tutorialmodfmap
2
”文件。
7.14.2
计算初始水头数据集
通过将
地面高程减少
10ft
作为初始水头值。
选择
“
Data
< br>”
菜单下的
“
Data Cal
culator
”
命令,
在
“
Expression
”
字段中输入计算公式
“
a-10
”<
/p>
,
在“
Result
”字段输入结果名“
starting head
”,单击
“
Compute
”按钮,单击“
Do
ne
”按钮退出。
7.14.3
高程和水头插值
选择<
/p>
“
Interpolation
”
菜单下的
“
to MODFLOW Laye
rs
”
命令,
对话框上部左侧为数据文
件名,
右侧为其相应的用途(
MODFLOW
< br>中的数据名)。
选择数据集名及其相应的用途名,
单击“
Map
”按钮,单击“
OK
p>
”进行插值。
7.14.4
导入层底高程离散数据文件
选择“
File
”菜单下的“
Import<
/p>
”命令,打开“
tutorialmodfmap
2
”文件。
7.14.5
层底高程插值
选择“
Interpolation
”菜单下的“
to
MODFLOW Layers
”命令,对应关系已经默认。
7.14.6
调整显示
选择
Select Scatter Point Sets<
/p>
工具
,
在离散点符号周围拖拉一个矩形框
,
选择
Hide
命令
< br>。
关掉网格外框(模拟边界)的显示:
转到
Map
模块
选项。
7.14.7
显示模型剖面
转到
3D Grid
模块
,选择靠近模型中心的一个单元,选择
View J Axis
命令
,要更好的
,选择“
D
isplay
”菜单下的“
Display
Options
”命令,关掉“
Grid frame
”
显示,可以调整垂向比例,选择“
View<
/p>
”菜单下的“
Z Magnification
< br>”命令,输入放大值,
OK
。
通
过由下方的两个水平方向的箭头可以看不同列的剖面。
7.14.8
修正高程值
p>
选择“
MODFLOW
”菜单下的“
Check
Simulation
”命令,单击“
Run
Check
”按钮,选择“
Fix
Layer Errors
”按钮。则出现错误信息,并提供了
修正方法。选择“
Truncate to
bedrock
”按钮,
单击“
Fix
Affected Layers
”按钮,
OK
,
Done
。
16
GMS
使用说明
同样通过
可以看水平切面。
7.15
检查模拟
在运行前应该运行
MODFLOW
Model Checker
多次,
选择“
MODFLOW
”菜单下的“<
/p>
Check
Simulation
”命令,单击“
Run
Check
”按钮,
7.16
保存工程
7.17
运行
MODFLOW
选择“
MODFLOW
”菜单下的
“
Run
MODFLOW
”命令。
7.18
显示水头等值线
选择“
MODFLOW
”菜单下的“
Read Solution
”命令。
改变等值线颜色:
选择“
Data
”菜单下的“
Contour
Options
”命令,单击“
Contour color<
/p>
”按钮,单击“
Color
”
窗口,选择颜色,
OK
。
p>
单击右下部的
按钮可以看第二层的水头等值线。
7.19
显示侧面上的水位线
在图中任选一点,单击
7.20
查看流量均衡
进入
< br>Map
模块
,选择
Select
Arcs
工具
,单击最右侧的排水沟弧段,则通过该排水沟的<
/p>
。
总流量显示在底部倒数第二条带上。
通过按下
SHIFT
键可以选择多条弧段,
显示单元的流量:
转到
3D
Grid
模块
7.21
添加注释
进入
Map
模块
,选择
Create T
ext
工具
,在添加注释的地方单击,输入文本,选择是
可以选择注释,并进行编辑修改。
可以选择
,选择一些单元,选择“
Data
”菜
单下的“
Flow
Budget
”命令。
否填充和字体颜色等。利用
Select Drawing O
bjects
工具
通过创建
Creat
e Line
线工具
,
绘制线,
然后利用利用
Select Drawing Objects
工具
该线,双击,在
Arrowheads<
/p>
区内选择适当的箭头。
第
8
章
MODPATH
8.1
问题描述
8.2
开始
8.3
所需模块
8.4
打开工程
选择“
F
ile
”菜单下的“
Open
”命令,
打开“
tutorialmodfmapsample
”文件。
8.5
初始化
MODPATH
模拟
进入
3D
Grid
模块
,选择“
MODPATH
”菜单下的“
New
Simulation
”命令。
17
GMS
使用说明
8.6
指定参数
转到
Ma
p
模块
工具
,
在
GMS
窗口中的
Coverages
组合框中选择
“
Layer 1
”
,
选择
Selec
t
Polygons
,双击该层多边形,选中“
Aquifer
Porosity
”选项,并输入值,
OK
。在
GMS
窗口中的
,
双击该层多边形,
选中
“
Aquifer
Coverages
组合框中选择
“
Layer 2
”
,<
/p>
选择
Select Polygons
工
具
Porosity
”选项,并输入值,
OK
。
8.6.1
为单元指定参数
选择“
Feature
Objects
”菜单下的“
Map ->
< br>MODPATH
”命令,
OK
。
8.7
指定起始位置
8.7.1
选择单元
进入
3D
Grid
模块
,选择
进入第二层,选择
Zoom
工具
,选择
< br>Select Cells
工具
,
选择单元。
8.7.2
创建起始位置
选择<
/p>
“
MODPATH
”
菜单下的
“
Generate Particles
”
命令,
选择
“
p>
Distribute starting points on cell
faces
”选项,输入点数,
OK
。
保存显示:
p>
选择
Frame
工具
;单击
返回上一层。
8.8
指定跟踪方向
选择“
MODPATH
”菜单下的“
Gene
ral Options
”命令,在对话框的顶部,选择“
Ba
ckward
”
选项,
OK
。
8.9
保存模拟
选择“
M
ODPATH
”菜单下的“
Save As
”命令,输入保存的文件名“
”,单击“
< br>OK
”
按钮。
8.10
运行
MODPATH
p>
选择“
MODPATH
”菜单下的“
Run
MODPATH
”命令,单击“
p>
Save
”按钮。
8.11
查看结果
选择“
MODPATH
”菜单下的“
< br>Read Solution
”命令。
在图中任选一点,单击
,可从侧面看。同样通过
可以看
水平切面。
8.12
从垃圾场追踪颗粒
8.12.1
改变跟踪方向
选择“
MODPATH
”菜单下的“
Gene
ral Options
”命令,在对话框的顶部,选择“
Fo
rward
”
选项,
OK
。
8.12.2
删除原来的颗粒
选择
Select Starting Locations<
/p>
工具
,
用鼠标拖拉选择示踪源,
按
Backspace
键或选择
“
Edit
”
菜单下的“
p>
Delete
”命令删除,
OK
,
OK
。
18
GMS
使用说明
8.12.3
指定新的起始位置
在
GMS
窗口中的
Coverages
< br>组合框中选择“
Recharge
”,选择
Select Cells
工具
,拖拉鼠标选择
p>
单元。
选择“
M
ODPATH
”菜单下的“
Generate
Particles
”命令,选择“
Distribute
starting points
on w.t. surface
< br>”选项,
OK
。
8.13
保存模拟
p>
选择“
MODPATH
”菜单下的“
Save As
”命令,输入保存的文件名“
”,单击“
OK
”
按钮。
8.14
运行
MODPATH
选择“
MODPATH
”菜单下的“
Run
MODPATH
”命令,单击“
Sav
e
”按钮。
8.15
查看结果
选择“
MODPATH
”菜单下的“
Read
Solution
”命令。
第
9
章
MT3DMS
进行
MT3DMS
p>
模拟之前必须先建立水流模型(根据前面所介绍的
MODFLOW<
/p>
使用说明)。
9.5
建立溶质运移模型
9.5.1
初始化模拟
选择“
MT3D
”菜单下的“
New
Simulation
”命令。
9.5.2
基本运移模块
p>
选择“
MT3D
”菜单下的“
Basic Transport
Package
”命令。
标题:输入注释说明及作者姓名、日期等
物种:
选择“
Define Species
”按钮,单击“
New
”按钮,给物种命名,
< br>OK
。
模块:
单击“
P
ackages
”按钮,选中“
Advection
Package
”选项,选中“
Dispersion Pac
kage
”选项,
选中“
Source
/Sink Mixing Package
”选项,
OK
p>
。
应力期:
单击
“
Stress Periods
”
按钮,
将
“
Length
”
字段改为
365
(模拟时间长度)
;
保证
Transport step
size
的值为
0
(
0
表示
MT3DMS
自动计算适宜的时间步长);
OK
。
输出控制:
单击“
Output
Control
”按钮,选中“
Print or save
at specified interval
”选项,指定间隔大小,
OK
。
浓度边界数组
ICBUND
:
初始浓度数组:
HTOP
和厚度数组:从
MODFLOW
模型中自动获取。
孔隙度数组:
9.5.3
对流模块
p>
选择
“
MT3D
”
菜单下的
“
Advection Pa
ckage
”
命令,
在
“
Solution schem
”
< br>下拉框中选择
“
Third
order TVD scheme (ULTIMATE)
”
,
OK
。
9.5.4
弥散模块
p>
选择“
MT3D
”菜单下的“
Dispersion
Package
”命令,单击“
Longitudinal D
ispersivity
”按钮,
选择“
Constant -> Grid
”选项,输入纵向弥散系数值,
< br>OK
;分别为“
Ratio of trans.
dispersivity to
19
GMS
使用说明
long.
dispersivity
”和字段“
Ratio of
vert. dispersivity to long. dispersivity
”输入横
/
纵和垂
/
纵弥散度
比例值。
9.5.5
源
/
汇项模块
p>
选择单元,选择“
MT3D
”菜单下的“<
/p>
Point Sources/Sinks
”命令,打开“
Well
”选项,输入
浓度,
OK
。
9.5.6
保存模拟
选择“
MT3D
”菜单下的“
Save
As
”命令。
9.5.7
运行
MT3DMS
选择“
MT3D
”菜单下的“
Run
MT3DMS
”命令。
9.5.8
读入结果
选择“
MT3D
”菜单下的“
Read
Solution
”命令。
9.5.9
改变等值线选项
选择“
Data
”菜单下的“
Contour
Options
”命令,打开“
Color fill
between contours
”选项,
OK
。
9.5.10
设置影象循环
选择“
Data
”菜单下的“
Film Loop
”命令,单击“
Setup
”按钮,打开“
Display
clock
”选项,
选择“
Use
constant interval
”选项,并将“
Time
interval
”的值改为所要的值,
OK
< br>。
第
10
章
MT3DMS
——概念模型方法
10.4
导入工程
p>
选择“
File
”菜单下的“
Open
”命令,打开“
tutorialmodf
mapsample
”文件。
10.5
初始化
MT3DMS
模拟
转到
3D
Grid
模式
,选择“
MT3D
”菜单下的“
New Simulation
”命令,选择“
MT3D
”菜
单下的“
Basic Transport
Package
”命令。
10.5.1
指定单位
单击“”按钮,选择适当的单位(
slug for
mass, ppm for
concentration
)。
10.5.2
定义物种
单击“
Define
Species
”按钮,输入物种名。
10.5.3
定义应力期
单击“
Stress
Periods
”按钮;输入适当的值。
10.5.5
设定输出控制
单击“
Output
Control
”按钮;选择“
Print at
specified time
”选项,单击“
Times
p>
”按钮,单
击“
Initialize V
alues
”按钮,输入适当的值;
OK
。
10.5.6
选择模块
单击“
P
ackages
”按钮;选择“
Advection
package
”、“
Dispersion
package
”和“
Source/Sink
Mixing package
”选项,
OK
。
10.6
指定含水层参数
10.6.1
按区赋参数
转到<
/p>
Map
模块
,
在
GMS
窗口中的
Coverages<
/p>
组合框中选择
“
Layer 1
”
,
选择
Select
Polygons
20
GMS
使用说明
工具
,<
/p>
双击该层多边形,
选中
“
Aquifer Porosity
”
选项,
并输入值,
选择
“
Lo
ngitudinal dispersivity
”
选项,<
/p>
并输入数值,
OK
。
在
GMS
窗口中的
Coverag
es
组合框中选择
“
Layer 2<
/p>
”
,
选择
Sel
ect Polygons
工具
,
双击
该层多边形,
选中
“
Aquifer
Porosity
”
选项,
并输入值,
选择
“
Longitudinal d
ispersivity
”
选项,并输入数值,
OK
。
10.7
指定补给浓度值
在
< br>GMS
窗口中的
Coverages
组合框中选择“
Recharge
”,选择
Select Polygons
工具
,双击该层多边形,在“
Recharge
”区的“
Contaminant
concen
tration
”文本框内输入浓度值。
OK
< br>。
10.8
转换概念模型
选择“
Feature
Objects
”菜单下的“
Map ->
MT3DMS
”命令,保证“
All applicable
coverages
”
选中,
OK
p>
。
10.9
层的厚度
MT3DMS
需要
HTOP
数组来指定最上部含水层的顶板高程;
并为各层输入厚度值。由于此处
采用“
true
layer
”技术,系统自动赋值。
10.10
对流模块
转到
3D Grid
模式
,选择“
MT3D
”菜单下的
“
Advection Package
”
< br>命令,
采用“
Method of
Characteristics (MOC)
”求解方案,<
/p>
OK
。
10.11
弥散模块
选择“
MT3D
”菜单下的“
Dispersion Package
”命令。纵向弥散系数自动赋值,其他三个参
数则需要人工赋值:
Ratio
of transverse disp. to long. disp
.
:
0.2
Ratio of
vert. disp. to long. disp
.
:
0.2
Effective
molecular diff. Coefficient
:
0
dispersivity ratios
:
0.2
进入第
2
层,输入同样的值。
10.12
源
/
汇项模块
此处已经由概念模型自动赋值。
10.13
保存模拟
选择“
MT3D
”菜单下的“
Save As
”命令,输入文件名。
10.14
运行
< br>选择“
MT3D
”菜单下的“
R
un MT3DMS
”命令。
10.15
查看结果
选择“
MT3D
”菜单下的“
Read Solution
”命令,
OK
< br>。在窗口顶部的
TS
组合框内选择最后一
步。
用颜色填充等值线:
选
择
“
Data
”
菜单下的
“
Contour Options
”
命令,
选择
“
Contour specified
range
”选项。
10.16
设置影象循环
21
GMS
使用说明
选择“
D
ata
”菜单下的“
Film Loop
”命令,单击“
Setup
”按钮,打开“
< br>Display clock
”选项,
选择“
Use
constant interval
”选项,并将“
Time
interval
”的值改为所要的值,
OK
< br>。
10.17
模拟吸附和衰变
10.17.1
选中化学反应选项
Chemical Reactions
Package
选择
“<
/p>
MT3D
”
菜单下的
“
Basic Transport Package
”<
/p>
命令,
单击
“
P
ackages
”
按钮,
选择
“
Chemical
reaction pa
ckage
”选项,
OK
。
10.17.2
输入吸附和生物降解数据
选择“
MT3D
”菜单下的“
Chemica
l Reaction Package
”命令,在“
Sorp
tion
”区中选择
“
Linear
isotherm
”
选项,
在
“
Kinetic rate reaction
”
区中选择
“
First-
order irreversible kinetic reaction
”
选项,
分别为各层输入数据
(层的选择通过对
话框下部的框架顶部带有上下箭头的列表框)
,
OK
。
10.18
保存模拟
选择“
MT3D
”菜单下的“
Save
As
”命令。
10.19
运行
MT3DMS
选择“
MT3D
”菜单下的“
Run
MT3DMS
”命令。
10.20
查看结果
选择“
MT3D
”菜单下的“
Read Solution
”命令,
OK
。在窗口顶部的
TS
组合框内选择最后一
步。
10.21
绘制时间历时曲线
10.21.1
创建观测图层
p>
转到
Map
模块
,
选择“
Feature Objects
”菜单下的“
Coverages
”命令,单击“
New
p>
”按
钮,输入图层名,将图层的类型改为“
Observation
”,单击“
Options
”按钮,单击“
New
”按钮,
< br>输入名称,
OK
。
10.21.2
创建观测点
选择
Create Point
工具<
/p>
,在窗口中单击添加点,并通过编辑框输入准确的坐标,保持该点
处于被选中状态,
选择
“
Featur
e Object
”
菜单下的
“
Attributes
”
命令,
选中
“
Assign to specified
layer
”
选项,并选择
layer
2
,
OK
。
10.21.3
创建时间序列图
p>
然后,
创建时间序列图:
保持该点处于被选
中状态,
选择
“
Display
”
菜单下的
“
Obs
. Plot Options
”
命令,单击“
New
”按钮,在“
Plot
type
”列表框中选择“
Time series
”,
OK
。
p>
最后,选择“
Display
”菜单下的“
Show Plot Window
”命令,选择
Frame
命令
。
10.21.4
绘制多条曲线
p>
选择
“
Display
”
菜单下的
“
Obs. Plot
Options
”
命令,
单击
“
Options
”
按
钮,
选择
“
Use selected
data set
”选项,选择要显示的物种名称,选中“
p>
Display
”选项。
10.21.5
移动观测点
选择
Select Points/Nodes
工具
,拖动点。
第
11
章
SEAM3D
SEAM3D
是一个用于模拟涉及多酶化作用和多电子接受剂的复杂生物降解的反应溶质运移
22
GMS
使用说明
模型。它基于
MT3DMS
程序,并添加了生物降解和
NAPL
溶解模块。
11.1
问题描述
模拟
NA
PL
(包括:苯和甲苯(
benzene and
toluene
))的污染羽的分布。
11.2
启动
GMS
11.3
所需模块和界面
11.4
导入水流模型
转到
3D Grid
模式
,
选择
“
MODFLOW<
/p>
”
菜单下的
“
R
ead Simulation
”
命令,
打开
MODFLOW
模拟文件(
p>
),选择“
MODFLOW
”菜单下的“<
/p>
Read Solution
”命令查
看
水流模拟结果,
OK
。
11.5
初始化
SEAM3D
模拟
选择“
M
T3D
”菜单下的“
New
Simulation
”命令。
SE
AM3D
、
RT3D
和
MT3DMS
都用
MT3D
菜
单。
11.6
基本运移模块
选择<
/p>
“
MT3D
”
菜
单下的
“
Basic Transport Package<
/p>
”
命令,
在
“<
/p>
Model
”
区内选择
< br>“
SEAM3D
”
选项。
标题:输入注释说明及作者姓名、日期等
11.6.1
指定单位
p>
长度和时间单位在
MODFLOW
中已经定
义,在此只需定义质量和浓度单位:
单击“
< br>Units
”按钮,选择适当的单位。
11.6.2
设定应力期
单击
“
Stress Periods
”
按钮,
将
“
Length
”
字段改为
2000
(模拟时间长度)
;
保证
Transport step
size
的值为
0
(
0
表示
MT3DMS
自动计算适宜的时间步长);
p>
OK
。
11.6.3
选择模块
单击“
Packages
”按钮,选中如下模块:
Advection Package
Dispersion Package
Source/Sink Mixing Package
Chemical Reaction Package
Biodegradation Package
NAPL
Dissolution Package
OK
。
11.6.4
指定物种
选择
“
Define Species
”
按钮,
将
“
Nondegradable Tracers
”
改为
2
,
将
“
Hydrocarbon Substrates
”
改为
2
,在“
Microbial Processes
”区内打开“
SO4
reduction
”选项,在“
Products to t
rack
”区内
打开“
H2S
”选项,单击“
New
”按钮,给物种命名,<
/p>
OK
。
在
p>
“
Names
”
列
表框内选择
“
Tracer1
”
,
并改变其名字为
“
Tracer (Conservative)
”
;
选择
“
Tracer2
”
,
并改变其名字为
“
< br>Tracer (Non-Conservative)
”
;
选择
“
Substrate1
p>
”
,
并改变其名字为
“
Benzene
”
;
选择“
Substrate2
”,并改变其名字为“
Toluene
”;
OK
。
11.6.5
输出控制:
单击“
Output
Control
”按钮,选中“
Print at
specified times
”选项,单击“
Times<
/p>
”按钮,单
23
GMS
使用说明
击“
Initialize
Values
”按钮,为“
Initial time
step size
”输入适当的值(
100
< br>),为“
Maximum time
step siz
e
”输入适当的值(
100
),为“<
/p>
Maximum simulation time
”输入适当的
值(
2000
),
OK
,
OK
,
OK
。
11.6.6
输入孔隙度
单击“
Porosity
”按钮,单击“
Co
nstant -> Grid
”按钮,输入适当的值(
0.2
5
),
OK
。
11.6.7
初始浓度
必须为每一物种赋初始浓度(其中
O2
(
4.0
)
, SO4
(
9.0
)
,
和
H2S
(
0.01
< br>)
要赋背景值):
在列表框中
选择
“
O2
”
,
单击
“
Starting Conc
entration
”
按钮,
单击
p>
“
Constant -> Grid
”<
/p>
按钮,
输入适当的值,
OK
;依次输入
O2
(
4.0<
/p>
)
, SO4
(
9.0
)
,
和
H2S
(
0.01
)
。
11.7
对流模块
通常要输入适当的值。但对于该例子,默认值“
Third
order TVD scheme
(ULTIMATE)
”即可。
选择
“
MT3D
”
菜单下的
“
Advection Package
”
命令,
在
“
Solution schem
”
下拉框中选择
p>
“
Third
order TVD
scheme (ULTIMATE)
”,
OK
。
11.8
弥散模块
选择“
M
T3D
”菜单下的“
Dispersion
Package
”命令,单击“
Longitudinal D
ispersivity
”按钮,
选择“
Constant -> Grid
”选项,输入纵向弥散系数值,
< br>OK
;分别为“
Ratio of trans.
dispersivity to
long.
dispersivity
”和字段“
Ratio of
vert. dispersivity to long. dispersivity
”输入横
/
纵和垂
/
纵弥散度
比例值。
11.9
源
/
汇项模块
p>
通过工具
选择单元,
选择
< br>“
MT3D
”
菜单下的
“
Point Sources/Sinks
”
命令,
打开
“
Constant
head
”选项,依次在列表框内选择
p>
O2
(
4.0
)<
/p>
, SO4
(
9.0
)
,
和
H2S
< br>(
0.01
)
并输入相应的背景
浓度,
OK
。
11.10
化学反应模块
Chemical Reactions
Package
选择
“
MT3D
”
菜单下的
“
p>
Chemical Reaction Package
”
命令,
在
“
Sorpt
ion
”
区中选择
“
< br>Linear
isotherm
”选项,输入比重
p>
Bulk density
(
5e-7
p>
);从列表框中选择“
Tracer (Non-Conserva
tive)
”,
输入“
1st
sorption constant
”值(
5.0e-9
p>
),同样给“
Benzene
”和“
Toluene
”赋值(都为
5.0e-9<
/p>
),
OK
。
11.11 NAPL
溶解模块
p>
对污染物从污染羽中逐渐溶解到地下水中,
在
MT3DMS
中可通过定浓度单元、
注水井或补给
来模拟,但是此法并不理想;
NAPL
溶解模块则
比较理想。
11.11.1
选择单元
单击
Select Cells
工具<
/p>
,选择“
Grid
”菜单下的“
Find Cell
”命令,输入单元的(
I,
J, K
)
值,然后拖拉选择所要的单元(污染羽中的单元)。
11.11.2
指定浓度
选择“
MT3D
”菜单下的“
Point Sources/Si
nks
”命令,在对话框的底部选择“
NAPL
”选项,
为“
Initial concentrat
ion
”输入适当的值(
2000
),
为“
Dissolution rate
”输入适当的值(
p>
3.0
),
OK
。
11.11.3
输入
NAPL
数据
选择“
p>
MT3D
”菜单下的“
NAPL
Dissolution
Package
”命令,将“
Number of
tracers in NAPL
”
24
GMS
使用说明
的值改为适当的值(
2
),将“
Number of
hydrocarbons in NAPL
”的值改为适当的值(
2
),在
对话框的底部为“
Ine
rt fraction molecular weight
”输入适当的值(
150
),依次选择左侧列表框
内各项,分别
输入相应的“
Initial mass fraction
”
、
“
Solubility
”和“
p>
Molecular
weight
”值(
0.1,
4400, 110
),
OK
。
11.12
生物降解模块(
Biodegradation
package
)
需要输入:给水度(
yield coefficients<
/p>
),阻碍系数(
inhibition coefficient
s
),和控制苯和
甲苯降解的参数。
选择“
MT3D
”菜单下的“
Biodegradation
Package
”命令,
11.12.1
最小浓度
选择
“
El. Acc
.
”
表单,
为
“
Number of bio steps per transport
step
”
输入值
(
< br>2
)
,
为
“
Microcolony
minimum concen
tration
”输入值(
0.001
);从文本框内选择“
Benzene
”并输入其
Minimum
concentration
最小
浓度值(
0.01
),从文本框内选择“
Toluene
”并输入其
Minimum concent
ration
最小浓度值(
0.05
)
,
从文本框内选择
“
O2
”
并输入其
Minimu
m concentration
最小浓度值
(
0.1
)
,
从文本框内选择<
/p>
“
SO4
”
并输
入其
Minimum concentration
最小浓度值
(
0.2
),
11.12.2
电子接受体系数
选择“
Min. Conc
.
”表单,为“
SO4-O2
”的“
Inihibition coefficient
”输入值(
0.5
),选择
下部文本框内每一项并输入它们的“
Yield
coefficient
”值(
O2-Benzene:
0.5
;
O2-Toluene:
0.4
;
SO4-Benzene:
0.15
;
SO4-Toluene:
0.1
)
11.12.3
生成系数
generation
coefficient
输入
H2
S
的生成系数:选择“
Gen.
Coeffs.
”表单,为“
Inihibition coe
fficient
”输入值(
0.15
)
11.12.4
利用系数(
use
coefficients
)
输入<
/p>
电子接受体
的利用系数:选择“
Use
Coeffs.
”表单,选择左上部文本框内每一项并输入
它们
的“
Elec. acceptor use
coeff.
”值(
O2-Benzene:
3.0
;
O2-Toluene:
3.2
;
SO4-Benzene:
8.0
;
SO4-Toluene:
8.2
)
11.12.5
饱和常数(
saturation
constants
)
选择
“
.
”
表单,
选择第一个文本框内每一项并输入它们的
“
Hy
dro. half sat. const.
”
值(
O2-Benzene:
1.0
;
O2-Toluene:
2.0
;
SO4-Benzene:
1.0
;
SO4-Toluene: 2.0
< br>);选择第二个
文本框内每一项并输入它们的“
Elec
. acc. half sat.
const.
”值(
O2:
0.1
;
SO4:
0.5
)
11.12.6
速率
选择“
Rates
”表单,确保对话框的“
Death
Rate
”部分的“
Calculated by model
”选项被选
中;选择第二个文本框内每一项并输入它们的“
p>
Max. specific rate of substrate utilizati
on
”值
(
O2-Benzene:
0.5
;
O2-Toluene:
0.05
;
SO4-Benzene:
0.005
;
SO4-Toluene:
0.0
)
11.12.7
初始浓度
利用默认值(
0
)。
11.13
保存和运行
选择“
< br>MT3D
”菜单下的“
Save As
< br>”命令。选择“
MT3D
”菜单下的“
< br>Run MT3DMS
”命令
11.14
读入结果
选择“
MT3D
”菜单下的“
Read Solution
”命令。
11.15
设置等值线选项
选择“
Data
”菜单下的“
Contour
Options
”命令,打开“
Color fill
between contours
”选项,
25
GMS
使用说明
单击“
Color Ramp
Options
”按钮,打开“
Show color
legend
”选项,
OK
11.16
查看浓度等值线
p>
首先看最后时刻(
2000d
后)
conservative tracer
的解:从
GMS
窗口的“
Data
Set
”组合框内
选择“
Tracer
(Conservative)
”数据集,从“
TS
”中选择时刻,该图则显示了没有吸附没有化
学反应的浓度等值线。它可以
作为比较其他污染物的基准。
通过
“
Data
Set
”组合框显示其他物种。
non-conservative
tracer
表示有吸附,但无反应。
11.17
绘制时间历时曲线
首先:
转到
Map
模块
,选择“
Feature
Objects
”菜单下的“
Coverages
”命令,单击“
New
”按
钮,输入图层名,将图层的类型改为“
Observation
”,单击“
Options
”按钮,单击“
New
”按钮,
输入名称,
OK<
/p>
。
然后:创建观测点
选择
Create Point
工具<
/p>
,在窗口中
NAPL
的下游单击添加点。
然后,创建时间序列图:
保持该点处于被选中状态,选择“
Display
”
菜单下的“
Obs. Plot Options
”命令,
单击
“
New
”
按钮,在“
Plot
type
”列表框中选择“
Time series
”,
OK
。
p>
最后,选择“
Display
”菜单下的“
Show Plot Window
”命令,选择
Frame
命令
11.17.1
移动观测点
。
选择
Select
Points/Nodes
工具
11.17.2
绘制多条数据线
,拖动点。
选择
“
Display
”
菜单下的
p>
“
Obs. Plot Options
”
命令,
单击
“
Options
”
按钮,
选择
“
Use selected
data se
t
”选项,选择要显示的物种名称,选中“
Display
p>
”选项,
OK
;打开“
Curve legend
”选项,
OK
< br>。
11.18
其他显示
选择“
Data
”菜单下的“
Film
Loop
”命令,可以制作影象循环放映。
第
12
章
定义层的数据
12.1
开始
选择“
F
ile
”菜单下的“
New
”命令。<
/p>
12.2
所需模块
< br>/
界面
Map
模块
;
2D Scatter
Point
模块
;
3D
Grid
模块
。
12.3
利用“
True
Layer
”模式
对于多
层
MODFLOW
模型,
要根据每个层
的类型来赋参数。
GMS
可自动计算某些参数。
层的参数
可以通过两种方式输入:直接输入数组或“
t
rue layer
”模式。对于“
true layer
p>
”模式,用户需要输
入顶底板高程及水平和垂直渗透系数。
12.4
通过插值
26
GMS
使用说明
利用
GM
S
中“
2D Scatter Point
”模块的“
Interpolation
”菜单下的“
p>
to MODFLOW
”命令对二
维离散数
据进行插值,并通过“
Model
Checker
”检查错误。
12.5
实例
12.6
创建网格
转到
3D
Grid
模式
Origin
Length
# Cells
500
2500
50
,选择“
Grid
”菜单下的“
< br>Create Grid
”命令,输入如下值:
Y-Dimension
500
1500
30
Z-Dimension
0
600
3
X-Dimension
选择
View I
Axis
工具
12.7
例
1
——完整层
查看剖面:选择
View K
Axis
工具
返回平面。
如图
12.1
所示,该模型包含
< br>3
层,且
3
层都扩展到整个区。
12.7.1
导入离散点数据
离散点数据文件可以通过
Microsoft Excel
p>
输入,
保存成以
TAB
分隔的文本文件
(扩展名为
)
p>
,
如图
12.2
。
其中
ID
列可以不要,上一层的底等于下一层的顶。
然后导入离散点数据文件:
选择
p>
“
File
”
菜单
下的
“
Import
”
,
打开
“
2
< br>”
文件。
12.7.2
对高程值进行插值
在插值前,需要初
始化
MODFLOW
数据:选择“
MO
DFLOW
”菜单下的“
New
Simulation
”命
令。
转到
2D Scatter Point
模块
,
选择
“
Interpolation
”
菜单下的
“
to MODFLOW Layers
”
命令,
将对话框左上部列表内数据集与右侧的
MODF
LOW
数组对应,单击“
Map
”按钮
进行映射,结果在下
部列表中,若要修改映射,可单击“
Unm
ap
”按钮,取消映射,然后重新映射。
27
GMS
使用说明
12.7.3
查看结果
转到
3D
Grid
模式
,选择某一位置的单元,然后选择
View I Axis
工具
查看所选单元处
的剖面:选择
View
K Axis
工具
12.8
例
2
——嵌入层
如图
12.3
返回平面。
12.8.1
导入离散点数据
离散点数据文件可以通过
Microsoft Excel
p>
输入,
保存成以
TAB
分隔的文本文件
(扩展名为
)
p>
,
28
GMS
使用说明
如图
12
.2
。其中
ID
列可以不要,上一层的
底等于下一层的顶,但是右边的数据需要调整:使第
2
层的底<
/p>
高于第
2
层的顶(插值后再修改)。
p>
然后导入离散点数据文件:
导入新的
数据前,
先删除原有的数据。
转到
2D
Scatter
Point
模块
文件。
12.8.2
对高程值进行插值
<
/p>
选择“
Interpolation
”菜
单下的“
to MODFLOW Layers
”命令,将对话
框左上部列表内数据集
与右侧的
MODFLOW
数组对应,单击“
Map
”按钮进行映射,结果在下部
列表中,若要修改映射,可
单击“
Unmap
< br>”按钮,取消映射,然后重新映射。
12.8.3
修正层的数据
转到
3D Grid
模式
,
选择
“
MODFLOW<
/p>
”
菜单下的
“
C
heck Simulation
”
命令,
单击
“
Run Check
”
p>
,
选择
“
Edit
”
菜单下的
“
Delete All
”
;
选择
“
File
”
菜单下
的
“
Import
”
< br>,
打开
“
2
”
按钮,左上侧列表框内显示错误信息,单击“
Fix
Layer
Errors
”按钮;在右侧列表框内选择“
Layer
p>
2
”,选择“
Average
”选项,单击“
Fix Selected Layer
”按钮,
OK
。
12.9
例
3
——含水层出露
如图
12.4
12.9.1
导入离散点数据
离散点数据文件可以通过
Microsoft Excel
p>
输入,
保存成以
TAB
分隔的文本文件
(扩展名为
)
p>
,
如图
12.2
。
其中
ID
列可以不要,上一层的底等于下一层的顶,但是右边的
数据需要调整:使第
1
层的底
高于第<
/p>
2
层的顶(插值后再修改)。
p>
然后导入离散点数据文件:
导入新的数据前,
先删除原有的数据。
转到
2D Scatter
Point
模块
文件。
12.9.2
对高程值进行插值
<
/p>
选择“
Interpolation
”菜
单下的“
to MODFLOW Layers
”命令,将对话
框左上部列表内数据集
与右侧的
MODFLOW
数组对应,单击“
Map
”按钮进行映射,结果在下部
列表中,若要修改映射,可
单击“
Unmap
< br>”按钮,取消映射,然后重新映射。
12.9.3
修正层的数据
29
,
选择
“
Edit
”
菜单下的
“
Delete All
”
;
选择
“
File
”
菜单下的
“
Import
”
,
打开
“
2
”
GMS
使用说明
p>
转到
3D Grid
模式
,
选择
“
MODFLOW<
/p>
”
菜单下的
“
C
heck Simulation
”
命令,
单击
“
Run Check
”
p>
按钮,左上侧列表框内显示错误信息,单击“
Fix Layer
Errors
”按钮;在右侧列表框内选择“
Layer
1
”,选择“
Preserve
top
”选项,单击“
Fix Selected Layer
”按钮,
OK
。
12.9.4
显示结果
选择
View I
Axis
工具
12.10
例
4
——基岩切割
如图
12.5
查看剖面:选择
View K
Axis
工具
返回平面。
12.10.1
导入离散点数据
离散点数据文件可以通过
Microsoft Excel
p>
输入,
保存成以
TAB
分隔的文本文件
(扩展名为
)
p>
,
如图
12.2
。
其中
ID
列可以不要,上一层的底等于下一层的顶。
然后导入离散点数据文件:
导入新的
数据前,
先删除原有的数据。
转到
2D
Scatter
Point
模块
文件。
12.10.2
对高程值进行插值
选择“
Interpolation
”
菜单下的“
to MODFLOW Layers
”命令,将对
话框左上部列表内数据集
与右侧的
MODFLOW
数组对应,单击“
Map
”按钮进行映射,结果在下
部列表中,若要修改映射,可
30
,
选择
“
Edit
”
菜单下的
“
Delete Al
l
”
;
选择
“
File
”
菜单下的
< br>“
Import
”
,
打开
“
2
”
GMS
使用说明
单击“
Unmap
< br>”按钮,取消映射,然后重新映射。
12.10.3
修正层的数据
转到
3D Grid
模式
,选择“
MODFLOW
”菜单下的“
Check
Simulation
”命令,单击“
Run
Check
”按钮,左上侧列表框内显示错误信息,单击“
Fix Layer Errors
”按钮;在右侧列表框内选择
< br>“
Layer
3
”,选择“
Truncate to
bedrock
”选项,单击“
Fix Selected L
ayer
”按钮,
OK
。
12.10.4
显示结果
选择
View I
Axis
工具
查看剖面:选择
View K
Axis
工具
返回平面。
,选择
Shade
工具
选择“
V
iew
”菜单下的“
General
Mode
”命令,选择
Oblique View
工具
,选择
Rotate
工
具
。
第
13
章
非稳定流数据管理
对于非稳定流
p>
,
需要大量的数据,
GMS
提供了适当的工具对它们进行输入和管理。
13.1
问题描述
13.2
开始
13.3
所需模块
3D Grid
< br>模块,
Map
模块,
MODFL
OW
界面
13.4
读入工程
选择“
F
ile
”菜单下的“
Open
”命令,
打开“
tutorialtrans
”工程。
13.5
非稳定流数据
通过概念模型获得。
13.6
p>
在
Map
模式下输入非稳定流数据
13.6.1
指定非稳定补给速率
转到<
/p>
Map
模块
,
选
择
Select Polygon
工具
,
选择一个区
(区
2
< br>)
,
选择
“
Feature Objects
”
菜单下的“
Attributes
”命令,选择“
Recha
rge Rate
”下面的“
Transient
”选项,单击“
Transient
”
下面的区域启动
XY Series Editor
.
编辑器,单击“
New
”按钮创建新
的时间序列,输入新的名字
(
zone2
)
,并输入新的时间数据序列;也可以通过“
Import<
/p>
”按钮将已经准备好的文本文件导入。
13.6.2
导入开采量数据
开采井
的数据文件包括两个文件:第一个文件包括井名、滤管几何形状、井的
XY
坐标。第
二个文件包括开采时间安排,其文件各格式为:
导入开采井数据:从
GMS<
/p>
窗口的
coverage
组合框中选择“
Sources & Sinks
”图层,选择“
File
”菜
31
GMS
使用说明
单下的“
Import
”命令,选择文件(
)<
/p>
,
Open
,
Y
es
。双击井单元,查看数据。
13.6.3
指定重力给水度
p>
在
map
模块下选择
Hydraulic
Conductivity
图层,选择
Select
Polygon
工具<
/p>
,选择一个区
(区
2
和
1
)
,选择“
< br>Feature Objects
”菜单下的“
Attr
ibutes
”命令,选中“
Specific Yield<
/p>
”并输入
适当的值(
0.2
)
;选择一个区(区
3
和<
/p>
4
)
,选择“
F
eature Objects
”菜单下的“
Attribut
es
”命令,
选中“
Specific
Yield
”并输入适当的值(
0.15
)
。
13.7
初始化
MODFLOW
13.7.1
将
MODFLOW
设置为非稳定流
p>
转到
3D Grid
模块
,
选择
“
MODFLOW
”
菜单下的
“
BCF Package
”
命令,
p>
在
“
Model Type
”
中选择“
Transient
”
。
13.7.2
指定模拟参考时间
选择
“
MODFLOW
”
菜单下
的
“
BCF Package
”
命令,
单击
“
Str
ess Periods
”
按钮,
在<
/p>
“
Reference
Time
”区下选择“
Edit
”按钮,输入初始时间
年月日时分秒(
000
)
。注意其下的
时间
显示选项,保证“
Date/time
”
。
13.7.3
设置应力期
单击<
/p>
“
Initialize
”
按钮,
输入
No. of stress peri
ods
(
7
)
、
length
(
100
)
和
No. of time steps
(
1
)
;
然后在
stress periods
对话框中修改
应力期及时间步数:通过对话框中左侧的带上下箭头的列表框
改变应力期数,
在右侧的
No. of time steps
文
本框内输入步数;
依次输入各应力期的结束时间和期
间的步数。
13.8
转换概念模型
转到
< br>Map
模块
,选择“
Featu
re Objects
”菜单下的“
Map -> MODFL
OW
”命令,
OK
。
< br>
13.9
保存和运行
MODFLOW
转到
3D Grid
模块
,选择“
MODFLOW
”菜单下的“
Save As
”和“
Run
MODFLOW
”命
令。
13.10
读入结果并显示
p>
选择“
MODFLOW
”菜单下的“
Read Solution
”命令。
13.11
创建影象循环
选择“
Data
”菜单下的“
Film Lo
op
”命令,单击“
Setup
”按钮
。
13.12
非稳定观测数据
13.12.1
导入非稳定观测数据
选择“
F
ile
”菜单下的“
Open
”命令,
打开扩展名为“
*.map
”的文件
<
/p>
。其格式为:
32
GMS
使用说明
13.12.2
< br>导入文件:选择“
File
”菜单下的“
Import
”命令,打开文件
trans_
。出现
Observation
Coverage
Options
对话框,在“
Measurement Typ
e
”列表中选择“
trans_head
”
(保证该种测量模式
被用于残生校正目标)
。
13.12.3
创建非稳定流观测曲线
转到
Map
模块
,选择“
D
isplay
”菜单下的“
Obs. Plot Option
s
”命令,单击“
New
”按钮,
p>
通过下拉列表框将图形类型设置为
Error vs time s
tep
,
OK
。选择“
Display
”菜单下的“
Show Plot <
/p>
Window
”命令,选择“
Frame
Image
”命令
。
p>
Me
为平均误差;
mae
< br>为平均绝对误差;
rms
为误差平方的平均值的平方根。
此图显示了测量值与模拟值的误差。
时间序列图:
转到
< br>Map
模块
,选择两个观测井,选择“
< br>Display
”菜单下的“
Obs. Plot Op
tions
”命令,
单击“
New
p>
”按钮,通过下拉列表框将图形类型设置为
Time
series
,选中“
Curve legend
”选项,
OK
。选择“
Di
splay
”菜单下的“
Show Plot
Window
”命令,选择“
Frame
Image
”命令
。
第
14
章
由区域模型转为局域模型(
Regional to
Local Model Conversion
)
p>
对于许多模拟研究,确定合适的边界条件是非常困难的。可通过先粗后细的方法分两步进行<
/p>
模拟。
问题描述:
如图<
/p>
14.1
所示:大部分边界为隔水边界,河流的出入口为已知水头
边界;两层;局域模
型位于其中心位置。
由区域模型转为局域模型的基本目标是从区域模型中创建一个包含水头和层参数的二维离
散点集,建立局域模型,并将水头和层参数插值到局域模型。基本步骤为:
建立区域模型并进行求解,利用“
MODFLOW
Layers -> 2D Scatter Points
”
命令从区域模型
中创建二维离散数据集,为局域模型创建三维网格,并为其进行插值。<
/p>
33
-
-
-
-
-
-
-
-
-
上一篇:全新版大学英语阅读教程答案完全版
下一篇:ANSYS函数