GMS软件学习日志

qiang · 发表于 2009-12-10 20:23

GMS学习笔记－1
1.岩性好像不能超过5个。（据同学说的，还没实践）
2.改等值线间隔：在Display Option里，的Contours选项按钮，在右侧的Contour Interval
里选择Specified Interval，在正面设置等值线间隔数值。
3.显示属性：Ctrl+D.
4.钻孔数据导入文件样例如：horizon文件夹下的holes.txt文件。
5.钻孔之间的连线可以通过顶点的TIN数据来拟合，即在做地层的顶底板时可以直接用相应的地层高程TIN数据生成地层厚度的变化情况。在编辑切面连线时先打断直线，然后选择导入的高程数据后再生成切面即可。
6.缩放工具，按shift后为缩小，否则为放大。
GMS学习笔记－2
基于TIN的表层模型

通过XYZ坐标数据，生成不规则三角网格。高程或其它基于TIN的数据可以通过等值线显示出。
TIN用来建立裸体模型和3D有限元格网。
1.导入顶点数据
a.切换到TIN模型
b.打开文件tutfiles\tins\verts.gpr
2.网格化
点击Build Tin下的Triangulate
3.建立等值线
点击Display Option，打开Contours和Tin Boundary选项，关闭Triangle edges和Vertices，即可看到等值线。
也可以用光源来显示顶面：在显示属性里关掉Contours和Tin Boundary，打开Triangle faces。切换到Oblique View，点击Display下的Lighting Options，把Ambient Light改到0.7后，即可看到效果。
4.编辑Tin
点击显示选项按钮，打开Vertices和Contours，在Contours右边点击Option，设置间隔为20，点击OK，OK返回。
a.拖动顶点
①在工具栏点击Select Vertices图标
切换到Plan View视图（x,y视图，更改X,Y坐标）
点击一个顶点可拖动更改位置。不能拖也边界。
在Oblique View视图可以拖动顶点以更改其Z坐标轴的位置。
②也可以点击一个顶点后，直接在上面的坐标显示框里更改坐标具体数值，按回车或tab键确认。
b.锁定顶点
选择好顶点后，点击Modify TIN下的Lock/Unlock Vertices命令即可。
c.添加顶点
切换到Plan View视图，在工具栏里选择Create Vertex工具
然后在图形里直接在要添加的位置点击即可。
d.删除
选择顶点后在菜单栏里点删除即可。默认情况下删除顶点后与此顶点相关的网格也会被删除，但可以设置选项开关使删除一个顶点后自动重新建立网格。
点击Modify TIN下的Vertex Options，打开Retriangulate after deleting选项即可。
5.平滑TIN
当顶点较少时，生成的等值线可能折线很显明，可以通过加密网格的办法使生成的等值线光滑。
a.复制顶点
点击Build TIn下的Tin->2D Scatter Points，点击OK，在询问是否删除已存在的TIN时选择“否”。
b.细分TIN
点击Modify TIN下的Uniformly Subdivide TIN，在弹出的窗口，增加Factor数会加密当前网格。
c.插值高程
加密网格后等值线并没有变化，通过以下步骤使等值线平滑：
切换到2D Scatter Point模型，
选择Interpolation|Interpolate->Active TIN。在弹出的对话窗口里可以选择插值方法。点击OK后即可看到插值后的平滑等值线。
GMS学习笔记－4
使用t-prog
1.读取钻孔数据
直接导入数据即可（位于Tutfiles\t-progs\LH 3D.gpr）
2.建立3D Grid
切换到Map模块
点击Feature Objects | Grid Frame，选择New Frame建立框架。以确定位置范围及方向。
此时会看到一个网格框架，其大小和方向可以通过两种方式编辑：1）编辑网格对话框里的数值，2）用鼠标拖动更改。其中右下角的圆圈用于更改方向。
在工具栏点击
后，点击新建的框架，可看到框架的编辑节点，此时即可用鼠标拖动节点修改，也可双击框架弹出框架的对话框，在里面输入精确的数值。
修改完后，点击Feature Objects | Map ->3D Grid命令，分别对X,Y,Z的Number cells输入70，50和20，点击OK。显示生成的3D Grid。
3.初始化T-progs模拟
下一步为初始化T-PROGS模型，并对方位角，背景岩性以及模拟的岩性进行赋值。
切换到Borehole模块，点击T-progs | New Simulation命令。弹出对话框对岩性进行设置。
方位角是基于X方向的。如果经XY平面上是各向异性，方位角要按各向异性来设置；如果不是各向异性，方位角应于网格的X轴（行或者J方向）相同。默认情况下，方位角与网格的X轴一致。这个值为我们输入到网格框架的相反数（前面输入的是40，此时显示的为－40），这是因为从X轴来讲网格的方位角是逆时针方向的，但从Y轴讲是顺时针方向。
在对话框里，第一列显示分析时用到哪些岩性。默认情况下，与钻孔相关的所有岩性都被选中。第二列为岩性的背景值，默认情况下以占有最大比例的岩性作为背景值。
点击Next。进入垂直Markov Chains对话框。
第一剖分为垂直转移趋势，基于钻孔数据。走向和倾向可以从垂直数据里得出。
点击左上侧的Compute按钮，利用GAMEAS工具计算转移概率。弹出计算过程对话框并进行计算，计算完成后，点击Close关闭。
返回后，下边各项数据都用计算后的结果进行了修改。在右下角显示的数值为转移概率值。
右键点击右上角的图表可以弹出图形菜单。
点击Edit Transition rates选项，可以直接编辑转移概率和平均比例（在正面的两个表单里）。
点击Edit embedded transition probabilities选项，GMS分析钻孔数据并计算嵌入的转移可能性，并在右侧表单里显示结果。
一般情况下选择Fit curves to a discrete lag选项，程序自动计算最合适的。
在Lag #里输入17并按Tab键，这个数值能很好的适应测量的转移数据和Markov Chains。
点击Next后，进入走向和倾向上的转移概率设置窗口。可以使用默认值直接点击Next。
4.运行TSIM
点击T-Progs | Run TSIm命令，把模拟名称改为sim3d，在Number of realization里输入5.保留其它默认值，点击Ok后进行计算，此时会弹出一个计算过程窗口，计算完成后点击close关闭。此时会激活MCMOD和TSIM工具进行计算，从MCMOD和TSIM计算的结果会显示在过程窗口的底部。点击CLOSE关闭。
5.查看结果
从T-PROGS计算的结果可以看到。
切换到3D Grid模块。在右边的数据树里可以看到计算的系列。点开Material Sets文件夹，点开sim3d文件夹即可看到。
6.第二阶段，单层网格
a.建立单层网格
切换到map模块，点击Feature Objects | Map ->3D Grid，点击OK。分别对X,Y,Z的Number cells输入70，50和1，点击OK。
网格层高程可以利用2D Scatter Point模块的to MODFLOW Layers命令从散点生成。
b.运行TSIM
利用同样的转移概率数值，因此直接运行TSIM。切换到Borehole模块，点击T-PROGS | Run TSIM命令。
取名为sim2d，为Number of realizations输入5.其它默认，此时会提示只有一层，每个钻孔在计算中都只认为是一个单样。当GAMESA计算结果时，关闭。
7.第三阶段，生成多层HUF数据系列
a.关掉单元表面，在属性里设置（在3D Grid模块下）
b.建立网格
切换到Map模块，点击Feature Objects | Map -> 3D Grid。点击OK。在Number of cells里分别输入70，50，4。OK。
c.初始Modflow
切换到3D Grid，点击Modflow | New Simulation。点击Packages，在Flow Package剖分，选择HUF包，点击OK，OK。
d.插值地层高程
通过导入一系列散点来进行插值。
导入数据：Open，选择topo.txt，打开Heading row。在GMS data type里确认选择2D Scatter Points。在Type行里，确保高程列（第四列）上为Data Set类型。
插值高程：切换到2D Scatter Point模块，点击Interpolation | to MODFLOW Layers命令。在MODFLOW Layer Arrays列表，选择Top Elevations Layer 1顶，然后点击Map按钮。
查看结果：在数据树里取消topo，切换到Borehole模块，在数据树里隐藏所有钻孔，切换到3D Grid模块，在网格中央附近选择一个单元，点击View I Axis按钮或View J Axis按钮查看（顶层）。
GMS学习笔记－5
插值地层数据
通过散点数据的高程来插值地层。利用MODFLOW的Model Checker来完成检查与修复。总共分四种情况。
1.完整的地层
完整的地层时，不用检查，直接插值即可。
点击Open，打开tutfiles\layerdata\points.gpr。
点击3D Grid模块，点击Select Cell工具，选择网格中居中位置的单元，点击View I Axis按钮切换到Front View。
在进行插值前，需要先初始化Modflow。点击MODFLOW | New Simulation命令。点击OK。
切换到2D Scatter Point模块，在右侧的数据树里确保case1为活动数据。点击Interpolation | to MODFLOW Layers命令。点击OK后，可以看到插值的结果。
2.具有尖灭时
按上述方法，切换到2D Scatter Point模块，确保Case 2为活动数据，点击Interpolation | to MODFLOW Layers命令，点击OK。
修正错误数据：切换到3D Grid模块，点击MODFLOW | Check Simulation命令，点击Run Check命令，可以看到有警告显示。点击Fix Layer Errors命令进行修正。
在Fix Layer Errors对话框里，有几种不同的方式来进行修正。此时在layer2里列出了许多错误。左侧显示的是修正选项。此时选择Average选项，此选项的意思是，计算第个单元格的顶和底高程的平均值，并把高于平均的高程移到平均值以上附近，低于平均的高程移到平均值以下附近。
在右侧选中Layer 2，选中Average选项，点击Fix Selected Layer按钮，点击OK，Done退出。可以看到修正的结果。
3.具有露头时
按前面同样的方法进行插值后，运行Check Simulation后，在Fix Layer Errors对话框里选择Preserve top选项进行修正。操作步骤同前。
4.基岩切割时
按前述方法，激活case4。在Check Simulation的对话框里选择Truncate to Bedrock选项进行地层修正。