简介：SURPAC是地质勘探、矿产资源评价与环境监测领域广泛应用的建模工具，这份Word文档系统梳理了SURPAC建模的核心知识点，适合初次接触SURPAC建模的工程师、地质技术人员及相关专业学生。内容涵盖项目数据在微机中的存贮组织、原始数据组织、剖面解译与中段解译流程，并重点讲解文件组织规范、命名方法及线串号分配规则；在剖面解译部分，详细介绍了勘探线、地层、矿体、小矿体、剖面坐标线、钻孔剖面及用于打印的剖面图等类型文件的处理要点，帮助读者理解建模过程中的数据分类与ID命名逻辑。资源共1个docx文件，压缩包大小约17KB，轻量便携，便于快速查阅。目前已有487人学习下载，适合作为日常备课、项目实操前的速查手册。
1. 用 SURPAC 软件建模：从钻孔数据到三维矿体的完整链路
“用 SURPAC 软件建模”这句话，外行听着像点一下按钮就能出模型，干过的人知道，它是一条从原始钻孔数据一路走到资源量估算的链路：建库、校验、组合样、圈定矿体、连接实体、填充块体并估值。最花时间的往往不是画线，而是数据整理和参数选择；最容易让结果翻车的也不是软件崩溃，而是坐标单位、缺样标记、搜索椭球体这些看不见的细节。这篇东西按一套主流工作流把每一步拆开，给步骤、给参数，也给踩坑记录。地质、采矿、资源估算方向的从业者都可以照着走一遍，从零把一个矿区的三维模型搭起来，并且知道什么时候该怀疑计算出来的数字。
2. 建模第一步：把散乱钻孔整理成可用的地质数据库
2.1 地质数据库为什么必须建在 SURPAC 内部
SURPAC 建模的地基不是 Excel 里的数据表，而是软件内部的地质数据库。虽然 SURPAC 允许从电子表格和 CSV 导入数据，但它能识别的不是一张大宽表，而是由孔口表、测斜表、岩性表、化验表四张互相通过 hole_id 关联的关系表组成的数据结构。
其中孔口表存钻孔的坐标和终孔深度，测斜表记录孔深方向变化，岩性表存分段岩性代号，化验表存分段品位分析结果。四张表各司其职，后面画剖面时按深度自动对应，做组合样时按钻孔钻孔分别处理，这些操作都必须依赖这种结构。我见过不少团队把所有信息塞进一张超宽表格的做法，数据是全的，但 SURPAC 无法把测斜和化验段按深度对齐，后面要么手动修正，要么重新建库，非常被动。
数据库表结构按常见做法可以梳理成这样：
| 表名 | 关键字段 | 作用 |
|---|---|---|
| collar | hole_id, x, y, z, total_depth | 孔口坐标与深度 |
| survey | hole_id, depth, azimuth, dip | 孔身轨迹 |
| lithology | hole_id, from_m, to_m, rock_code | 岩性分层 |
| assay | hole_id, from_m, to_m, grade | 分段品位值 |
如果说坐标是空间基准，hole_id 就是数据血缘。只要某个钻孔在四张表里的孔号命名不一致，这个孔就会被当成独立的新孔，后续一切操作都会把它带偏。建库时建议先做一张孔号清单，四张表都从这个清单派生，比事后慢慢排查快得多。
2.2 从 CSV 导入后的第一轮体检：用 SQL 查重复孔和异常品位
数据从各方汇集过来，经常是不同人整理的版本，字段名不完全一致，单位也可能混用。我的做法是先把各分表的 Excel 另存为 CSV，再在 SURPAC 的数据库导入向导里做字段映射，字段名一律改成英文小写，避免之后写查询语句时被中文字段绊住。
导入完成后不要急着画剖面，必须先做一轮完整性检查，直接在 SURPAC 的数据库查询窗口里执行标准 SQL 是最快的。下面这段脚本是我每次建库都会跑一遍的体检：
-- 检查孔口表是否有重复孔号 SELECT hole_id, COUNT(*) AS cnt FROM collar GROUP BY hole_id HAVING COUNT(*) > 1; -- 检查化验表是否有负品位 SELECT hole_id, from_m, to_m, grade FROM assay WHERE grade < 0 ORDER BY hole_id, from_m;第一段查询把孔口表里重复的 hole_id 列出来，重复孔会造成后续实体模型连接时线串错乱。第二段查负的品位值，这个值在数据库里通常代表缺样标记，正常化验结果不可能是负数。排查时要区分“真缺样”和“零品位”，两者对后期估算的影响完全不同。
2.3 样品组合：一场必须做的按长度加权平均
原始化验样长度是不一致的，有的样长 0.5 米，有的样长达 3 米，如果直接拿原始样段参与统计，长样段会对品位分布统计产生过度影响。SURPAC 的做法是先按一个统一组合长度把孔深切段，再在每个组合段内按长度加权平均，这个操作业内叫样品组合，可以把它理解为“把不等长的样品归一化成等长的品位序列”。
常见做法是选 1 米或 2 米作为组合长度。薄矿脉用 1 米，厚大矿体用 2 米。在用 SURPAC 的组合功能正式跑之前，我一般会先用一段标准 SQL 做一次预组合，检查数据合理性，也顺便确认缺样是否被正确排除：
-- 按 2 米组合长度做加权平均品位 SELECT hole_id, FLOOR(from_m / 2.0) AS combo_id, SUM(grade * (to_m - from_m)) / NULLIF(SUM(to_m - from_m), 0) AS avg_grade, SUM(to_m - from_m) AS combo_len FROM assay WHERE grade IS NOT NULL AND grade >= 0 AND to_m > from_m GROUP BY hole_id, combo_id ORDER BY hole_id, combo_id;这段 SQL 的逻辑说明：FLOOR 把起始深度除以 2.0 并向下取整，得到组合编号，比如从 1.2 米到 3.1 米的样段会归到编号 0 和编号 1 两个组合里；avg_grade 用每个样段的长度作为权重做加权平均，等价于“品位乘以长度再除以总长度”；NULLIF 处理除零，防止空洞段造成计算结果异常；WHERE 里剔除了空值、负值和反向段。跑完这一步，如果发现组合后的长度远超 2 米，就要回头查原始样段是否有跨孔连接的问题。
组合样做完还要处理特高品位，常见做法是按品位分布的 98 百分位进行截断，或者按均值加 3 倍标准差做封顶。这一步直接影响后续估算结果，却经常被当作“可以后面再说”而跳过，实际上错过这个封顶点，一块高品位矿段就可能把整片块体模型的品位抬高一个等级。
3. 实体模型：在勘探线剖面上圈矿体并把剖面连成三维
3.1 线串与实体：SURPAC 建模的两层数据结构
数据库准备好之后，建模进入核心阶段：圈定矿体边界并生成三维实体模型。SURPAC 里要先把矿体边界数字化成线串，再把多条线串连成一个封闭的实体模型。线串是一系列有顺序的三维坐标点连接而成，可以闭合也可以不闭合；实体模型则是由若干三角形面构成的封闭外壳，才是真正能算体积的对象。
这两层结构经常被混为一谈。有人把矿体的每层边界画成了线串，就想直接算体积，结果是软件只能计算三角形网格面积，体积永远是 0。必须先确认线串已经封闭，再通过三角网连接生成实体。90% 的建模时间其实都花在线串的编辑调整上，实体生成反而是最后一个操作，确认了这层关系，后面才不会在软件里乱找按钮。
3.2 在勘探线上数字化矿体：从显示设置到闭合线串
线串数字化要在剖面视图里完成。打开剖面显示功能，填入当前勘探线的方位角和宽度，接着把钻孔加载进来，再按品位阈值给化验段上色。颜色分级的间隔建议先按区间下限设，比如铜矿按 0.2%、0.4%、0.8% 分四档，这样可以直观看到矿化边界在哪里。
随后新建一条线串文件，沿着矿体边界逐点点击，转弯处加密节点，直线段可以适当放宽，最终必须把首尾点连起来形成闭合线串。保存这条线串时用统一命名规则，比如 BODY-1-W1 表示 1 号矿体西侧第一条剖面边界。完成一个剖面后切回三维视图，把线串和钻孔轨迹叠加显示，确认它和矿化段的上下盘位置吻合。
这个过程有两点必须注意：第一，矿体边界要对应地质规律，不是把品位高于阈值的样段全部框进去，要考虑岩性、构造和矿化连续性的综合因素；第二，同一矿体在相邻剖面上的边界形态要尽量相似，突变太大时后面连接三角网必然出问题。这一步属于典型的“软件好教、经验难学”，初学者往往会圈出一个边缘极其复杂的边界，导致后续实体模型难以收敛。
3.3 相邻剖面连三角网：节点顺序、封闭性检查与体积计算
所有剖面线串都画好后，进入实体模型生成阶段。在 SURPAC 里选中同一矿体相邻的两条剖面线串，执行连接三角网操作，软件会自动在两条线串的对应节点之间生成三角面片。生成效果取决于线串的基本形态，节点数差不多、走向一致、闭合方向也一致时，三角网会很规整。
最容易出的问题是线串方向不一致。一条剖面顺时针闭合，另一条逆时针闭合，连接出来的三角面就会交叉扭曲，表面看着如乱麻。解决方法是选线串时先做方向检查，把所有边界线串统一成同一方向的闭合方向。生成实体后必须做封闭性检查，SURPAC 会列出模型的开放边数量，正常应该是 0。如果存在开放边，体积计算就会漏掉一部分区域，检查工具会把所有出问题的边高亮出来，逐个修复后重新检查。
体积计算放在实体模型封闭之后做才有效。我一直保留的习惯是：算完体积后，把实体模型的多层剖面切出来和原始线串叠在一起看一遍，确认没有飘出钻孔范围，再继续做块体模型。实体阶段出现偏差，后面块体模型再精细都会被整体带偏。
4. 块体模型与品位估算：把实体变成可计的储量块
4.1 块体尺寸：主块、次分块与矿体产状的匹配
实体模型是连续的外壳，块体模型则是往这个外壳里填充一个个立方体单元，每个单元都要被估算出一个品位值。块体尺寸和矿体产状、工程间距必须匹配，这是整个建模流程中最容易被“先随便设一个”的参数。
块体尺寸过大会把边界切割得坑坑洼洼，尺寸过小又会让计算量激增。按行业常见做法，主块尺寸一般取工程间距的三分之一到二分之一，比如钻孔间距 50 米，主块可以取 10 米或 15 米。薄矿脉要单独处理，主块太厚会把围岩混进矿石里，建议沿矿体厚度方向加密。
| 矿体类型 | 主块尺寸 | 次分块尺寸 | 适用条件 |
|---|---|---|---|
| 薄矿脉 | 2 m × 2 m × 2 m | 0.5 m × 0.5 m × 0.5 m | 坑道密集取样 |
| 中等厚度矿体 | 5 m × 5 m × 5 m | 1 m × 1 m × 1 m | 探槽加钻孔结合 |
| 厚大矿体 | 10 m × 10 m × 10 m | 2 m × 2 m × 2 m | 钻探网度 50 m |
次分块的作用是让块体模型贴近实体边界。主块填不满边界时，软件会把边界处的块继续细分，直到能拟合实体表面。次分块尺寸一般取主块的二分之一到四分之一。这一层设置直接决定边界附近的资源量归属，忽视它会造成边界块品位归属混乱。
4.2 变异函数参数：块金值、基台值、变程如何指导搜索半径
如果选用克里格法做品位估算，变异函数分析就是绕不开的环节。SURPAC 的变异函数模块先按实验变异函数把数据点算出来，再拟合一条理论曲线，常用的拟合模型是球状模型。球状模型有三个核心参数：块金值、基台值、变程。
| 参数 | 含义 | 常见判断 |
|---|---|---|
| 块金值 C0 | 小尺度上的随机变异 | 过高说明数据质量差或矿化随机性强 |
| 基台值 C0+C1 | 空间相关性的总方差上限 | 用于归一化处理 |
| 变程 a | 相关性消失的距离 | 决定后续搜索半径的基准 |
三个参数里，块金值的解释最需要谨慎。块金值太高，往往不是矿体本身的问题，而是化验数据里有大量缺样或错样造成的假随机性。变程则直接决定后续品位估算时的搜索半径，一般取变程的 1 到 1.5 倍，超过这个范围搜索进来的样品基本没有参考意义。沿矿体走向、倾向、垂直三个方向要分别计算变异函数，因为矿体在这三个方向上的连续性差异很大。
很多人把变异函数当成黑匣子，跑完一条曲线就直接用，其实判断这个结果是否合理有一个很朴素的标准：沿矿体走向的变程应该明显大于垂直方向，如果走倾向和走向的变程几乎一样，说明数据分布存在系统性问题，先回头检查钻孔位置是否过于集中。
4.3 选 IDW 还是克里格：参数表与常见取值
距离反比加权和克里格是 SURPAC 里两种最常用的品位估算方法。IDW 按距离的 n 次方加权，简单直接，不依赖变异函数模型，适合数据稀疏或勘探前期快速出结果；克里格利用变异函数做无偏最优估计，同时输出估值方差，适合正式资源量估算。
| 参数 | IDW 常用值 | 克里格常用值 |
|---|---|---|
| 搜索半径 | 变程的 1～1.5 倍 | 变程的 1～1.5 倍 |
| 最小样品数 | 4 | 4～6 |
| 最大样品数 | 16 | 12 |
| 每个扇区最大样品数 | 4 | 2～3 |
| 扇区数 | 4～8 | 4～8 |
搜索椭球体的方向和矿体产状必须一致。矿体走向 45 度、倾伏角 20 度，搜索椭球体的主轴就不能仍然保持正东正北。这个参数经常被按默认设置直接跑，结果高品位块体的位置明显偏移到钻孔连线之外，后续剖面切出来对不上，返工成本很高。扇区数的作用是防止某一方向的样品过度集中，把一个方向的样品数压得太高会让估值偏向局部高值，扇区数增加会降低这种偏差。
5. 建模避坑：五个让模型翻车的细节
5.1 孔口坐标一个错，模型整体错位
现象：实体模型和实际钻孔轨迹在三维视图里差出几十米，矿体跑到地形外面。
原因：最常见的是单位混用。孔口表一个字段写的米，另一个字段用的是英尺，SURPAC 按统一单位导入后，某个方向的坐标被放大了对应倍数。这类问题不容易靠肉眼看出来，因为整体几何形状是完整的。
解决：导入前用 SQL 查一下最大最小坐标范围，对比已知矿区范围的合理尺度，发现异常立刻停。建立一张孔口坐标抽查表，随机选几个钻孔和测量成果比对单孔坐标，校验通过后再进入建模。
5.2 缺样值被当成零品位，矿体品位被稀释
现象：组合后品位偏低，块体模型里大片区域出现低于边界品位的数值，但钻孔编录明明显示那是矿石。
原因：缺样段在原始表格里被填成了 0，导入后无法从“零品位”和“真缺样”之间区分出来，组合脚本把缺失段当成无矿段参加加权平均。
解决：导入前把所有缺样统一改成 -9999，并在后续 SQL 和 SURPAC 组合参数里设置“忽略负值”。组合样算完后统计每个组合段的有效样长比例，如果大量组合段的有效长度远小于理论长度，就说明缺样处理不到位。
5.3 组合长度与矿体厚度不匹配，品位被平滑过头
现象：薄矿脉的块体模型品位分布变得很均匀，峰值明显被压低，边界也模糊。
原因：组合长度取得比矿体厚度还大，一个组合段同时跨了矿石和围岩，加权平均自然会把品位拉平。
解决：先统计矿体厚度分布，再定组合长度。矿体平均厚度 1.2 米时，组合长度取 0.5 米或 1 米，而不是默认的 2 米。同时查看组合后品位的直方图，如果分布曲线变得过于规则，就要缩短组合长度重新算。
5.4 实体封闭性检查过不了，体积算不准
现象：实体模型生成成功，但计算体积时结果明显偏大或偏小，和实际开挖量对不上。
原因：实体模型存在开放边，软件在计算体积时对未封闭的部分做了不合理填充，或者把某些边界重复计算。
解决：每次生成实体后运行封闭性检查，把所有开放边高亮显示，逐条修复。修复原则是先补漏点，再调整三角形连接方式，避免用一个大三角面强行覆盖复杂边界。封闭性检查结果必须为 0 时才允许进入体积计算环节。
5.5 搜索椭球体主轴没跟矿体走向，高品位位置偏移
现象：品位估算完成后，块体模型的高品位区域整体偏离钻孔矿化段的位置，剖面上看高品位块跑到钻孔下方十几米。
原因：搜索椭球体的方位角、倾伏角仍然保持系统默认，没有按矿体的实际产状做旋转。主轴方向错了，搜索范围沿错误方向延伸，估算出的品位分布自然跟着偏移。
解决：先做一盘矿体产状统计分析，把走向、倾向、倾角整理出来，把搜索椭球体的主轴方向设置成矿体走向，倾伏向设置成矿体的整体倾伏向。跑完估算后再切一个沿走向的剖面验证，叠加钻孔品位看分布是否对齐。
6. 进阶验证：用剖切面与钻孔对比给模型“体检”
块体模型算完，先别急着出表，做一次剖切面对比再提交结果。这个习惯帮我发现了不少参数配置问题，属于投入最小、收益最高的验证手段。
切片验证的具体操作：选中已完成的块体模型，设定一个剖面方向，比如按矿体走向切成垂直剖面，在剖面上叠加钻孔轨迹和化验品位分级显示，然后把块体模型的估算品位按同样的色阶显示。两者对比时，主要看三个点：矿化富集中心是否重合，高品位区是否在边界处被截断，以及边界品位的过渡是否连续。任何一处有明显错位，都要倒回去查搜索椭球体方向和扇区参数。
剖面方向不要只做一条。沿走向、倾向各切两刀，再看一个水平切片，检查矿体在平面上的品位展布。没有钻孔控制的外推区，要核对实体模型是否有不合理膨胀。这一类检查本质上是在和野外一手数据做对拍，比任何统计报表都直观。
做完切片验证，再分别用 IDW 和克里格各估一遍，对比两套结果的吨位量和平均品位差异。差异小于 3% 说明模型稳定；差异太大则说明数据支撑不足或搜索参数不合理，应当以地质规律为依据确定取信哪一套结果。这个差距也是向协作方解释不确定性的一个定量依据。
我做模型的习惯是每一步都留中间文件，剖面线串、实体模型、块体模型分版本存档，绝不覆盖。参数调整全部记录在案，这样排错时能快速定位是哪个环节引入的问题。资源量估算不是一次成型的事，多留一个版本就多一条后悔药。希望这套流程对你有帮助，祝建模顺利。
本文还有配套的精品资源，点击获取
