使用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按钮查看（顶层）。