GoCAD三维地质建模全流程：从数据准备到网格导出

上篇我写到了 GoCAD 的安装、基本界面认识和钻孔数据准备，很多朋友在评论区提到“数据已经导进去了，下一步该干嘛”，说实话这个阶段是最容易迷茫的。数据一堆堆躺在树上，工具面板密密麻麻，鼠标点来点去不知道从哪下手。这篇“下”我就接着上篇的进度，把从“有数据”到“出模型”的完整过程拆开揉碎讲清楚：断层怎么建、层面怎么生成、SGrid 网格怎么搭、属性怎么往网格里塞，以及最后怎么导出成果。

这篇文章的定位是给已经装好 GoCAD、也把基础数据整理得差不多的朋友看的。如果你还在安装和导入阶段卡着，先去翻上篇，把几个核心坑填平再回来。下面的内容全程基于 GoCAD 2015 及以上版本的菜单结构，不同版本界面表述略有差异，但核心逻辑完全一致。我会尽量按真实操作顺序走，关键是解释清楚每一步“为什么这么做”，因为你只有理解了意图，才能在遇到报错时知道怎么改。

1. 建模前的数据准备与格式转换

1.1 坐标系统统一与投影带核查

很多人建模建到一半发现层面扭曲得不能看，十有八九是坐标系统出了问题。这里我必须再说一次：GoCAD 不会主动帮你去判断坐标系统，它显示什么就是什么。如果一个钻孔用的是经纬度（单位是度），另一个钻孔用的是投影坐标（单位是米），放在同一个视图里，点位会跑到天上去。

我自己的习惯是在导入之前先用 Excel 或者 Python 把数据检查一遍。看三列：X 值范围是不是在一个合理的区间内，Y 值范围是不是正常，Z 值有没有出现明显的异常（比如负值代表高程基准面以下）。如果坐标数据里混了多个投影带，哪怕都是米制，模型也会出现错位，这时候唯一的办法就是回到原始数据重新处理。

具体来说，项目里常见的是地方独立坐标系或者国家标准坐标系。无论哪种，只要所有数据源统一到同一个坐标系下就行。检查方法也很简单：把钻孔的孔口坐标和地形点云画在一起，如果两者没有重叠，先别急着调模型，回头查数据源。

提示：整张表里出现个别点位飞到千万级数值，多半是经纬度被人为乘了比例系数或者加了带号前缀。宁可多花半小时清理，也不要拿脏数据进模型。

1.2 钻孔表结构设计与分层数据清洗

GoCAD 的 Well 结构虽然支持你直接右键新建，但更稳妥的做法是在外部把数据表做扎实。我的通用模板分成两张表：一张是孔位表（Well：Well ID, X, Y, Z, 深度, 孔深），另一张是分层表（Well Marker：Well ID, 层位名, 顶深, 底深, 岩性编码）。两张表通过 Well ID 关联，导入后用 Query 查询功能做联查。

这里有一个容易被忽略的细节：分层表中的顶深、底深，一定要用“孔深”而不是“高程”。孔深是从地面往下量的长度，高程是从海平面往上算的数值，两者相差一个孔口高程。GoCAD 在处理 Well Marker 时，默认你给的是深度值，它会用孔口高程换算成 Z 坐标。如果你把高程直接填进去，那个标志点就会被放到完全错误的位置。

数据清洗方面，最常见的问题是层位名不统一。“①层粉质粘土”“粉质粘土①”“ZC-1层”这种写法如果不做归一化，后面建层位的时候就等着一个个改吧。我的做法是准备一张字典表，把现场记录的别名统一映射成标准层位代码。这个工作看起来笨重，但能避免你建出来一个层位图层里出现几十个孤立小碎块。

1.3 点云与曲面数据的格式转换实操

地形数据、层面散点数据，经常以 txt、csv、Surfer GRD 或者 ESRI shapefile 的形式提交。直接拖进 GoCAD 不一定被识别，需要做一步转换。

我自己最常用的是把 txt/csv 转成 GoCAD 的 ASCII 格式，也就是 .gp 文件。格式非常简单，第一行写对象类型，后面跟坐标点列。手工整理容易错，推荐用一段小脚本完成。下面给一个参考，基于 Python 3（pandas）：

import pandas as pd # 读取原始点云或者散点文件 df = pd.read_csv("input_points.csv", encoding="gbk") # 只保留需要的三列 df = df[["X", "Y", "Z"]].dropna() # 写出 GoCAD ASCII Points 文件的头部 with open("output.gp", "w", encoding="ascii") as f: f.write("GOCAD PPoints 1\n") f.write("HEADER {\n") f.write("name:GroundPoints\n") f.write("}\n") f.write("TFACE\n") for _, row in df.iterrows(): f.write("VRTX 0 {} {} {}\n".format(row["X"], row["Y"], row["Z"])) f.write("END\n") print("转换完成，共 {} 个点".format(len(df)))

这个脚本生成的文件可以直接用 File > Import > GOCAD ASCII 导入。如果你的数据里有重复点，建议先用 df.drop_duplicates() 处理一下；如果点特别密（超过几十万个），可以先做抽稀，比如按坐标网格每 0.5 米留一个点，否则后面生成曲面时卡到你想哭。

Surfer GRD 或者 TIF 高程模型也可以导入，GoCAD 提供了专门的导入过滤器，导入后自动生成 Surface 对象。要是导入后表面朝向反了，记得用 Surface 菜单里的 Invert Orientation 翻转一下，这不影响几何形状，但会影响后续体积计算和剪切操作的正负号。

2. 核心建模流程：断层建模与层面构模

2.1 从断层数据到断层面对话框的完整操作

断层建模是地层建模里最“师傅领进门”的环节。如果你手里有断层解释数据，比如地震解释的断层多边形、野外测量的断层面产状，那么目标就是生成一张或多张连续、互不交错的断层面。

GoCAD 里面建断层有两条路：一条是直接从点数据拟合，称为 Surface from points；另一条是使用专门的 Fault 对象，从断点数据构建。

我先讲通用思路。在 Objects 目录下右键新建一个 Surface 对象，然后利用 Edit > Surface Edition > Create surface from scattered points，选中断层散点数据。这时候对话框里有个关键参数叫网格大小（Grid Size）。网格太大，断层面会像一块块平板拼的；网格太小，边缘会出现大量毛刺。我一般先根据点密度估算一个值，比如点间距的中位数乘以 1.5，跑一版看效果再做微调。

断层建模还有个特殊参数叫“走向延伸长度”。这是因为野外实测断层往往只露出地表一小段，但地下深处它是延续的。你不能让断层面在边缘处突然截止，那样后面切割地层时会把地层切出一个突兀的豁口。解决方法是把断层线沿走向方向外推一定距离，这个距离我习惯取断层长度的 20%-30%，具体看你要不要保留断层尖灭形态。

注意：断层建模最忌讳把多条产状不一致的断层揉进同一个 Surface 里。一个 Surface 只能代表一个断层面，要是数据里出现多个断层，分别建，后面再设置交切关系。

2.2 层位数据插值：带断层约束的曲面生成

层位建模的核心思路是用钻孔里的分层点、地表的地形点，拟合出每一个层位的顶底面曲面。这听上去像数学插值题，但地层建模讲究的是“构造合理”，不是“数值精确”。

GoCAD 建层面用的工具在 Surface 菜单下，叫 Surface from points with faults。为什么特别推荐带断层的版本？因为如果你先建了断层面，再用这个工具，它可以保证生成的层面在断层两侧是断开的，而不是硬生生穿过断层表面。这一点在逆断层或者高角度正断层地区极其关键。

实操时，你需要在对话框里指定三样东西：

• 参与插值的散点集，选择某个层位的所有分层点；
• 断层约束，勾选对应的断层面对象；
• 边界范围，选择一个闭合的 Polygon 或者 Rectangle 限定拟合区域。

参数区有两个值需要重点调：一个是平滑因子（Smooth factor），另一个是容差（Tolerance）。平滑因子我默认给 0.5，值越大曲面越圆润但越偏离原始点，值越小曲面越贴点但越毛糙。容差则控制拟合曲面上的点到原始点的最大允许偏差，这个我一般给 0.1 到 0.3 米，如果跑出来超差了就说明原始数据里有打架的点，回头查分层表。

生成好的层面对象名称建议带前缀，比如 “Top_Quaternary”、“Base_Quaternary”，别用默认的 surface-1 这种名字。项目模型做到后期会有几十个对象，命名混乱直接导致自己都找不到数据。

2.3 层面相交与边界裁剪的处理

多个层面生成之后，它们天然是互相穿越的。比如第四系底界和三叠系顶界，理论上在所有位置都应该是叠置关系，但因为插值算法不同，拟合出来两个面下部可能有交叉。

处理这种相交问题的标准动作是用一个“边界混合面”作为辅助。假设你把它理解成一块海绵蛋糕的顶部和底部：顶部曲面和底部曲面通过同一个侧边曲面围起来，就形成一个封闭层壳体。GoCAD 里这个操作叫 Closed Surfacic Model，也就是用若干个 Surface 围成一个封闭体。

我做区域地层体时，会先画一个矩形范围，然后用 Project/Clip 把每个层位面都裁剪到这个矩形范围内，再给各个层位面之间添加侧面封边面，最后组合成封闭体。这个环节很容易出“不闭合”的错，排查方法也很直观：用 View > Show Boundary 查看每个面的边界线，如果边界线之间没有完全重合，就没法闭合。

3. 网格化：从层面体到 SGrid

3.1 SGrid 的基本逻辑与行列深度设置

到了这一步，你已经有了一个由曲面构成的“壳子”，下面要做的是往壳子里塞网格，这就是 GoCAD 里的 SGrid 对象。你可以把它想象成一个魔方：魔方整体占据一块体积，但内部被切成了整齐的小方块。每个小方块就是一个单元（Cell），后续的岩性、孔隙度、渗透率等属性值就挂在单元上。

右键新建 SGrid 时，第一个对话框会让你指定区域范围和层位约束。区域范围需要你框选一个边界多边形（Polygon），这个范围不能超出你已经建好的层面范围，否则网格会出现空腔。层位约束就是你希望以哪些曲面作为网格顶底边界，通常从上到下按地层顺序依次选入。

行列深度设置有三个数字：I 方向网格数、J 方向网格数、K 方向网格数。前两个是水平方向的划分，第三个是垂直方向的层数。层数不是越密越好，我的建议是：水平方向上控制每个单元在 10 米到 25 米左右，垂直方向上每个地层至少分 3-5 层。之前做过一个小区块，面积约 1.5 平方公里，我设成 40×40×30，总共 4.8 万个单元，跑起来完全不卡。你要是做到百万级单元，每次旋转视角都能明显感到卡顿。

3.2 网格的约束控制与断层截断处理

SGrid 对话框里有个容易忽略但极其重要的选项，叫 Boundary Adjustment 或者类似的“贴合边界”功能。它会让边界上的网格单元自动收缩或者扩张，以贴合你的控制面（断层和层面）。如果你不做这一步，生成的网格是规整长方体的，断层位置会出现台阶状锯齿，后续属性建模和体积计算都会受到影响。

断层截断的意思，是让网格单元在断层两侧完全断开，两侧各自独立编号。操作上是在 SGrid 创建向导的 Fault 选项卡里勾选对应的断层面，然后指定“断开”模式。这样生成的网格，在做储量计算时，就能分别统计断层上盘和下盘体积。

这一步出现的典型报错是“Unable to compute valid grid due to constraint conflict”，本质上是断层面和层面之间的拓扑关系有冲突。排查方法：把断层面和被约束的层面都显示出来，用剖面视图（Clipping Plane）沿着断层走向切一刀，检查面与面之间是不是有交错。如果有，先回到层面建模阶段把相交关系理顺，再重新生成网格。

3.3 网格质量检查：正交性、负体积与单元退化

网格生成后不要急着填属性，先做一次网格质量体检。GoCAD 里可以查看单元体积分布，我习惯用颜色渲染单元体积，重点检查有没有体积为零或负值的单元。负体积单元本质上就是单元绕了：网格节点的顺序反转了，这类单元会导致后面的储量计算和数值模拟全面崩盘。

另一个指标是单元正交性。网格单元理想情况下是长方体，但贴断层和层面时，单元会被压扁或者拉伸。当单元的内角偏离直角太多时，数值模拟软件会报告收敛困难。我给自己定的标准是：平均正交性不低于 0.4（1 是完全正交），低于 0.3 的区域必须局部加密网格重新生成。

如果检查中发现某一层被严重压缩，比如厚度只剩 5 厘米，这通常是你生成层面的时候把两层曲面做到几乎贴在一起了。这时候别在网格里硬调，回到层面模型把曲面的 Z 值关系修好，再做网格。

4. 属性建模：把化验数据塞进每个网格

4.1 属性对象的创建与数据关联方式

SGrid 建好了，里面是空的，现在要往里填属性。属性建模的目的是给每一个网格单元赋予一个物理值，比如岩性编码、含水率、渗透率、品位等。这个属性信息来自稀疏的钻孔样品，因此本质上是“用已知点推断未知点”的数学过程。

GoCAD 的做法是右键 SGrid 对象，新建一个属性（Property）。属性创建后，你有两个输入来源：一个是钻孔的属性曲线数据（Well Log），另一个是离散的点集数据（比如化验品位分析点）。

我推荐把样品数据整理成一张带 XYZ 坐标的表格，然后用 Point Set 导入。这样不但可以用于属性插值，还可以先做出散点视图，直观看到高品位区或者低渗透区的空间分布趋势，再据此判断要不要调整插值参数。

4.2 插值算法选择：克里金、距离反比加权与最近邻

GoCAD 属性建模的插值算法主要在 Properties 菜单下。我用得最多的有三种：克里金（Kriging）、距离反比加权（IDW）、最近邻法（Nearest Neighbor）。

克里金其实是个家族，里面细分有简单克里金、普通克里金、协同克里金。对于地质属性建模，普通克里金是最常用的。它的最大优点是利用了变差函数（Variogram）来描述空间相关性，理论上给出的估值方差最小。缺点是上手门槛高，变差函数调不好，结果反而不如 IDW。

距离反比加权（IDW）就直白多了，离未知点越近的已知点权重越大，公式是权重的 p 次幂。p 值给 2 是最常见的，p 越大，远距离点的影响衰减越快，模型就越“局部”。

最近邻法，基本就是把每个单元赋予给它最近的那个已知点的值，结果看起来像马赛克。我只在给岩性这类离散变量做分区赋值时用它，连续型物性参数不用这个方法。

我自己的经验是：连续型参数（如孔隙度、渗透率、品位）先试普通克里金，如果变差函数拟合效果差，就退回 IDW p=2；离散型参数（岩性、风化程度分级）用最近邻法或者直接按距离场分类。下表是三个算法的适用场景对比：

4.3 变差函数的直观理解与调参经验

提到变差函数，很多初学朋友直接懵了，我试着用大白话拆解一下。变差函数做的事情很简单：统计所有已知点对之间，距离越近的点对，属性值是不是越接近。把你的已知点投影到某个方向上去，每隔一定距离数一数这些点对的方差，把方差画成随距离变化的曲线，就是实验变差函数。

这里有两个关键参数你必须懂：一个叫变程（Range），它表示在这个距离范围内，点与点之间才有相关性；超过这个距离，两边数值就完全是独立随机了。另一个叫基台（Sill），是变差函数在变程处趋于平稳的那个值，大致对应数据的整体方差。

调参时我一般这样做：先让 GoCAD 自动计算实验变差函数，看曲线的大致形态，再用球状模型（Spherical）拟合。如果你发现变程特别小，比如只有 10 米，但钻孔距离已经到 50 米了，那说明钻孔网度太稀，插值结果可信度很低，这时候最优解是加密钻孔或者采集更多样品，而不是硬调参数。

提示：变差函数是有方向性的。如果地层有明显的层理，沿层理方向的变程往往远大于垂直方向的变程。这时候要设置各向异性，不能拿同一个变程用到所有方向，否则插出来的属性会把层理抹平。

4.4 插值结果验证与交叉验证方法

属性插值完，第一件事不是高兴，而是验证。我有个习惯：把 10%-15% 的样品点抽出来不参与插值，插完之后再把这些点位的估值和实际化验值对比。算一下绝对误差和相对误差，误差过大就要回头调参数，或者怀疑数据里有录入错误。

GoCAD 里还有更省事的方法，就是用 Cross Validation 工具，它会自动轮流把每个已知点当作未知点去估，输出每个点的估值误差。误差大的点位会高亮显示，你顺着去看，常常能抓到异常值或者分层错误。

我经历过的典型案例：某个钻孔的样品位数在一个位置突然冒出一个高得离谱的品位值，之前散点图没仔细看，交叉验证时那个点误差爆表，一查原始记录发现是样品编号错行。这就是不做验证就会带进模型的系统性错误。

5. 模型可视化、出图与成果导出

5.1 显示样式、色标与三维场景调优

模型建完之后，关键是把三维视图调到一个能看的程度。这里所谓“能看”不是审美问题，而是要让不同层位、不同属性值在屏幕上一眼就能区分开。

在 GoCAD 里，每个对象都有独立的显示属性。配置步骤是：右键对象 > Properties > Visualization Parameters，调整颜色、透明度和光照。

我的常见做法是：断层显示为红色半透明，层面用带透明度的渐变色，属性网格用色标渲染。色标设置的技巧是不要用彩虹色阶，推荐用地质行业通用的蓝-绿-黄-红序列，冷色代表低值、暖色代表高值。你还要注意设定数据范围，不要让色标被几个异常值拉满，否则正常区域的颜色差异会非常小。这时直接在色标属性里锁定最小最大值，或者去掉超出数据的 5% 分位数。

三维场景调优方面，记得打开垂直比例（Vertical Exaggeration），尤其在平原地区，地形起伏只有十几米时，不做垂直放大根本看不出层次。我常设在 5-10 倍之间，视项目而定。

5.2 剖面切割与标准图件导出流程

剖面图是地质报告中必不可少的图件。GoCAD 里切割面用的是 Clipping Plane 功能。你会画一个剖面位置线（Line），然后以它为基准生成一个垂直剖面，系统会把所有对象沿着这个面切开显示。

要导出一张清晰可用的剖面图，我的流程是：先把视图调成剖面视角，然后关闭网格轴线和多余的对象，只保留要展示的层面、断层和属性网格。接着设置 View 的背景颜色为白色（默认灰色容易让打印效果脏兮兮的），打开图例（Legend）显示色标和比例尺。

最后用 File > Export > Image，选 PNG 格式，分辨率设 300 dpi，导出的图片直接能放进报告。不要用截图键截图，那样分辨率不够，放大就糊。剖面图在报告里通常需要标注孔号和孔口高程，如果 GoCAD 里嫌麻烦，我一般都在 CAD 或 PPT 里二次叠加标注，效率更高。

5.3 成果导出：SGrid 输出与通用格式转换

建模的最终目的是多端使用。模型要给数值模拟软件用，就得导成塔层（LGR）或者角点网格格式；要给绘图软件出图，就导成 DXF；要给别人只读浏览，就导成 PDF 或者图片。

GoCAD 里最常用的导出选项包括：

• SGrid 导出成 CMG、Eclipse 等数值模拟格式，导出前检查单元数与属性命名是否规范，避免一方软件里出现空格或者非法字符。
• 层面导出成 Surfer GRD 或者 ESRI TIN，用于二维绘图软件叠加分析。
• 点集和井数据导出成 CSV，回算给数据库系统归档。

有一个细节，导出到数值模拟软件之前，记得做一次属性离散化。原生的连续属性每层可能有几十个数值分段，数值模拟软件通常要求每个网格单元只能有一个代表值。GoCAD 里可以用 Property 菜单下的 Upscale 工具，按哪个主要体积占比的原则给每个单元赋一个代表值。这个环节偷懒的话，模拟结果肯定是要被返工重搞的。

注意：每次导出前都要重新检查一遍属性所在的空间范围。右键属性对象查看 Z 范围，如果出现比钻孔深度还深得多的值，说明插值算法在局部外推过度了，这类区域要么裁剪掉，要么把超范围值置为无数据（Void），不要带进下游。

6. 常见问题与排查技巧实录

6.1 高频报错信息速查表

自己在建模过程中踩坑无数，也帮同事排查过不少问题，整理了一张高频问题速查表。这里列的“报错信息”是表现层看到的提示，原因是实际逻辑，解决办法是我亲自验证过有效的。

6.2 性能卡顿与崩溃处理经验

大场景数据在 GoCAD 里拖动视角非常慢，这是所有桌面版三维软件的通病。我的优化心得大概可以归纳为三条：抽稀显示、分块建模、关闭全局渲染。

先说抽稀显示。很多卡顿源于你把百万级点集全部显示在屏幕上。可以在显示设置里把点的显示方式改成每隔 N 个点显示一下，或者只显示包裹体（Bounding Box），等到需要精细处理再切换回来。

再说分块建模。如果你负责的是几十平方公里的大区域模型，不要试图把全区的数据一次性做成一个 SGrid，那样电脑再强也会很吃力。按区块或者按层位分别建立子模型，最后用 Data Merge 工具拼起来。宁可拼接时多花功夫，也不要让主程序直接跑十几小时最后崩溃。

最后，操作 SGrid 这类大对象时，把动态透视（Picking）功能关掉。有时候你只是想旋转视角，GoCAD 却在实时计算鼠标下面那个单元格的属性，CPU 瞬间就满了。在 Settings 里把 Pick 功能设为手动，需要查询单元数据时再打开。

6.3 安装与授权常见坑

最后必须提一嘴安装和授权，因为后台私信被问爆了，基本集中在“双击图标没反应”“打开闪退”“显示无法连接授权服务”这三类。

“双击没反应”大部分是安装目录有中文路径或者权限不够，把安装目录换到纯英文路径，比如 C:\GoCad2015 下。“打开闪退”常见原因是显卡驱动太老，去显卡官网更新 OpenGL 驱动即可。“无法连接授权服务”则先检查本机 License Manager 服务有没有启动。在 Windows 服务管理器里找到对应的许可证服务，右键启动，再重启软件基本能解决。

提示：如果软件启动后一直卡在初始化界面，大概率是系统里残留了旧版本配置，删掉安装目录下的 cache 文件夹再重启一次，注意这个操作不会影响你存放在其他路径的模型文件。

写在最后的体会

这套流程前前后后我自己跑了不下五六个项目，从最初的钻孔表整理到最终的属性网格导出，每一步其实都有绕弯的地方。要说印象最深，一定不是某个高深的算法，而是“数据质量决定模型质量”这句话的分量。你花两周调出来的变差函数，可能还不如花两小时把钻孔分层表里的错层位改掉来得有效。所以我始终建议，无论你想把 GoCAD 玩得多花哨，先在数据清洗上多投入时间。建模软件只是一个翻译器，把现场钻探、测试、测绘的数据翻译成一个三维的地质认知框架，输入的数据错了，后面再怎么调参数也只是在把错误装饰得更漂亮。

如果你手头正在跑自己的项目，卡在某一步不知道什么意思，欢迎带着截图和数据描述来问。下一篇我预计会写一写模型精度评估和不确定性分析，毕竟三维地质模型的价值到底有多大，不能光凭“看着像”说了算。