一个砂岩型铀矿区的三维地质建模，第一版交上去就被退了回来：剖面图上一条砂体横穿隔水层，没有断开。查了两天，软件没问题，插值参数也没问题，问题出在开工前——三个批次的钻孔，用了两套坐标系。
一、返工的账，都记在建模前那两三天
项目规模不算大。3 km² 的建模范围，286 个钻孔，9 条勘探线剖面，成图比例尺 1:2000。第一版模型提交后，剖面上三处砂体与钻孔柱状对不上，最深的一处差了 11 m，甲方一眼就看出来，退回意见只有一句：砂体不该穿过隔水层。
代价很实在。数据要重新导表、重新投影、重跑剖面矢量化，模型要重算变差函数、重调搜索椭球，前后搭进去两周；按四个人算，等于八个工日的有效工时被抹掉。而这批钻孔的数据，本来在开工第一天下午就能查完。
把工作量摊开看，三维地质建模这件事里，建实体和插值只占一小半，大头在数据。钻孔、剖面、平面图、样品化验、测井解释，五类数据来自五个部门、跨四个年份，格式从 Excel 到 CAD 再到纸质扫描件，全都有。【作者测算：按近几年经手的 20 余个矿区项目统计，因数据前置问题产生的返工时间约占总工期的 30%–40%，非官方数据】
判断门槛其实很低：清点数据时，只要发现来源超过 3 个年份，或者坐标系来源记录缺失，就先别建实体，把体检做完。这套体检也挡不住一种情况——原始钻孔没有测斜数据、孔斜全靠经验修正，这种孔越往深部偏得越多，体检阶段看不出来，只能标记为低置信度，建模时降权处理。
二、坐标系统与孔深闭合：最便宜的两项
坐标系统的反例最好认。把全部孔位投到平面图上，如果东南角有一片孔整体往西挪了百米量级、横竖都对不上勘探线，基本就是北京 54 与西安 80 混用了。另一种是投影带号写错，表现是孔位整体沿经线方向被拉开，越往东偏差越大。
做法是先在表里加一列，逐批记录坐标系来源和投影带；再按批次分组批量转换；转完抽 5 个孔，用实测 GPS 点校核。平面偏差超过 5 m 就回头查参数，不要靠整体平移去凑。这一招在早期测量记录不全的批次上会失效：只有整片坐标、没有分孔原始观测时，只能整片平移，残余误差得靠邻孔对比慢慢吸收。
这一列来源记录属于钻孔数据整理的第一道工序，缺了它，后面所有校核都无从谈起。甲方问建模需要什么数据，第一项答案就是它：每一批孔的坐标系、投影带、高程基准，都要落成可查的一列。
孔深闭合的错也一眼能看出来：柱状图上孔口标高 1186 m、孔深 420 m，终孔标高却标着 768 m，中间那 2 m 就是凭空消失的。批量查更快，直接算「孔口标高 − 孔深 − 终孔标高」，不为零的行全部筛出来。单孔闭合差超过 0.5 m，或全矿区超过 1% 的孔存在闭合差，就回退原始班报表重新录入，别用程序平差糊过去——平差能让数字对得上，却把真实的孔深错误一并抹平。
三、倾角方位、样品重复、地层代号
倾角方位错的表现很直观：剖面生成出来，钻孔全都笔直竖着。砂岩型铀矿的矿体常沿层间氧化带展布，孔斜哪怕只差 3.5°，到 400 m 深度上就是二十几米的水平位移，砂体边界直接跑偏。测斜记录里有、库表里没录，出来的结果一样歪。
做法是把测斜记录逐孔对上，看终孔坐标与设计孔位差了多少；与测斜记录差超过 1°，或者终孔偏移超过 5 m，就按测斜曲线重建孔轴，不要拿设计孔斜当实测。补测斜按孔计价，一根孔几百元，比重画一条剖面便宜得多。
样品重复的表现是同一深度段出现两条分析结果，品位一个 0.031%、一个 0.068%，谁对谁错没人说得清。成因多半是补采样时只加了新行、没删旧行。按采样编号和深度段做双键去重，重复段保留化验日期较晚、送样单位一致的那条；重复段长度超过 1 m 就整批重导，别手工删几行，删掉的往往是最该留的那条。
地层代号统一最琐碎。同一层砂岩，三个年份的报告里可能分别写成 K₁、下白垩统砂岩、T3，建模软件会把它们当成三套互不相干的地层，实体一建就是一摞薄饼。代号种类数超过原始报告层级数的 1.2 倍，就先建一张映射表，一列旧名、一列新名、一列依据来源，定稿签字之后再导库。这张表挡得住反复改口，挡不住新钻孔带来的新叫法——每进一批新数据，还得回头补一次。
四、剖面套合与拓扑错误
剖面数据矢量化是返工率最高的一环。反例是这样：把剖面的平面位置叠到孔位图上，剖面线没有穿过孔口符号的中心，偏出去三十来米，剖面上画的地层界面和钻孔柱状自然对不上，越往下越离谱。
做法是先定剖面线的平面坐标，再把孔口投影到线上，投影距离写进属性表；扫描剖面按控制点配准、按图面比例尺换算深度；收尾再抽查一遍孔位套合。套合只能管住线位与孔位的关系，剖面线本身画错了，这一步查不出来。
拓扑错误在三维里表现为悬边、自相交和空洞，二维阶段看不出来，一建实体就露馅。常见成因是矢量化时线头没接上，或者同一条界面画了两遍。七项检查里，只有这一项必须等三维地质建模的实体出来才看得见，前置检查只能查一半。
判定看抽样结果：抽 10 条剖面，若有 2 条以上孔位套合偏差超过图面 1 mm（1:2000 下约 2 m），或者单层拓扑错误面超过 3 处，整批重做矢量化，不要在建模软件里逐点修补。修补的成本比重做高，补过的面后期还没法追溯。七类检查的门槛与处置动作汇总成一张表，可以直接当台账用：
表1 七类数据体检的判定阈值与处置动作
| 检查项 | 一眼可见的错 | 判定阈值 | 处置动作 |
|---|---|---|---|
| 坐标系统 | 孔位整片偏移、与勘探线错位 | 平面偏差 > 5 m | 按批次回查投影参数后重投影 |
| 孔深闭合 | 孔口标高 − 孔深 ≠ 终孔标高 | 单孔 > 0.5 m | 回退原始班报表重录 |
| 倾角方位 | 剖面上钻孔近乎垂直 | 与测斜记录差 > 1° | 补录测斜或标记降权 |
| 样品重复 | 同一深度段两条分析结果 | 重复段长度 > 1 m | 按采样编号与深度双键去重 |
| 地层代号 | 同一层出现多种叫法 | 种类 > 层级数 1.2 倍 | 建映射表签字后导库 |
| 剖面套合 | 剖面线未穿孔口符号 | 图面偏差 > 1 mm | 整批重做矢量化 |
| 拓扑错误 | 悬边、自相交、空洞 | 单层错误面 > 3 处 | 修面合格后再建实体 |
五、关键判断清单
- 数据来源跨 3 个年份以上、坐标系来源无记录 → 先做体检，不建实体。
- 全矿区闭合差超标孔数超过 1% → 停建，回退原始资料重录，不做程序平差。
- 剖面套合偏差 2 条以上超限，或单层拓扑错误超过 3 处 → 重做矢量化，不进建模软件修补。
- 地层代号种类超过层级数 1.2 倍 → 先定映射表并签字，再谈插值与变差函数。
- 无测斜数据的钻孔占比超过两成 → 模型整体降为概略级，交付时写明置信度和影响范围。
六、边界与风险
体检这套流程管不了深部。钻孔控制深度以下的界面靠地质认识外推，数据再干净也换不来精度。甲方要求 300 m 以深的地层界面误差控制在 10 m 以内时，靠现有钻孔密度做不到，该谈的是加密钻孔或物探约束，不是继续调插值参数。
另一类风险在数据权属和时效。三年前的钻孔、五年前的剖面和今年的新增孔拼在一起，中间几年补勘的孔可能已经改变了原有地层划分方案，体检查得出代号不统一，查不出哪个方案才是现行——这得找矿区地测科书面确认。
还有一种走不通的情况：原始剖面只有纸质扫描件，控制点模糊到配准残差超过 3 mm。这种剖面硬做剖面数据矢量化，误差会直接传进三维地质建模的模型里，比不建更糟。代价是要么放弃该剖面、把模型范围缩小，要么补做实测剖面。
七、落地建议
- 第一步，开工前抽 1 天，把七类检查做成一张 Excel 台账，286 个孔逐行过，红灯项当天退回数据提供方，并记下退回时间。
- 第二步，坐标系、地层代号两张映射表定稿并签字，作为后续所有三维地质建模流程的唯一数据依据；建模需要什么数据的争议，到这里一次性谈完。
- 第三步，先拿 10 条剖面里的 2 条做试点，套合与拓扑过关后再批量矢量化，钻孔数据整理的脚本挂进同一套流程，避免整批返工。
- 第四步，把体检台账作为交付附件之一，甲方复核时有据可查，下一轮更新模型不用从零再查。
八、结语
九成的三维地质建模返工，都落在开工前那两三天。七类检查不需要额外软件，一张台账、一遍投影、几次抽查就能过完。真正的门槛不是技术，是别人催着出图的时候，肯不肯先把数据按住不动。体检做完，再谈三维地质建模流程里的插值与变差函数才有意义。
