4. 地质剖面数据整理

剖面绘制的质量，八成取决于数据整理是否扎实。露天采矿的剖面图不是“画”出来的，而是由钻孔数据库经界线提取、品位插值后“算”出来的。本节围绕一条完整的数据链展开：钻孔数据库建立 → 地层与矿体界线提取 → 品位数据插值 → 地质剖面图生成，重点讲清每一步的输入、关键处理与结果判读。

一、钻孔数据库建立：把原始记录变成可计算的结构

钻孔数据通常来自地质编录、测斜、化验与岩性描述四类原始资料。建立数据库的核心是把它们统一到“孔口坐标 + 测斜轨迹 + 分段属性”的结构上。工程上常用三张表组织：孔口表（collar）记录孔号、X/Y/Z 坐标、孔深、开孔方位与倾角；测斜表（survey）记录各测点深度及方位角、倾角；分段表（interval）记录从某深度到某深度的岩性、地层代号、品位等属性。

关键处理是孔迹计算。以测斜数据为基础，按中角度法或最小曲率法把“沿孔深”的采样位置换算为三维坐标。例如某孔孔口坐标 (X=500.00, Y=1200.00, Z=180.00)，开孔倾角 -60°、方位角 90°，第一测段长 30 m，则该测段终点近似为：ΔX = 30·cos(-60°)·sin(90°) = 15.00 m，ΔY = 30·cos(-60°)·cos(90°) = 0，ΔZ = 30·sin(-60°) ≈ -25.98 m，即终点约 (515.00, 1200.00, 154.02)。以上为示例参数，实际应以测斜表逐段递推。

建库时还要做三项校验：孔号唯一且与化验单一致；分段区间首尾相接、不重叠、不超孔深；品位低于检出限时统一记为“未检出”而非 0，避免插值时被误当作真实零值。校验通过后，数据库才能作为后续界线提取与插值的唯一数据源。

二、地层与矿体界线提取：从分段属性到空间界面

界线提取的本质是找出相邻钻孔之间“岩性/地层代号发生变化”的深度位置，并把这些位置在剖面上连成界面。做法分两步：先在单孔内识别界线点——某分段表中 0～12 m 为第四系、12～35 m 为强风化岩、35 m 以下为矿体，则 12 m 与 35 m 即为两个界线深度；再按孔迹计算把界线深度换算为三维点坐标。

随后在剖面方向上，把同一界面的界线点按空间顺序连接。连接不是简单直线，而要遵循地质规律：界面在剖面上应连续、不穿越自身、不出现无依据的突变；当相邻孔间距较大（如 100 m 以上）时，宜结合产状趋势做合理外推，并在图上标注“推断界线”。矿体界线还需区分地质界线与可采矿体边界：前者由岩性控制，后者还要叠加工业指标（如边界品位）约束，二者在剖面上可能不重合，必须分别提取、分别成图。

三、品位数据插值：把离散样品变成连续品位场

钻孔样品是沿孔轴离散分布的，剖面图需要的是沿剖面线的连续品位。常用做法是沿剖面线插值：先确定剖面线方程，把邻近钻孔的样品按垂直剖面方向投影到剖面线上，形成“距离—品位”样本序列，再选用插值方法。

工程上常用反距离加权（IDW）：待估点品位 = Σ(wᵢ·gᵢ)/Σwᵢ，其中 wᵢ = 1/dᵢ²，dᵢ 为待估点到第 i 个样品的距离。示例：某待估点距三个样品分别为 20 m、40 m、60 m，品位 1.2%、0.8%、0.5%，则 w 分别为 0.0025、0.000625、0.000278，加权品位 ≈ (0.0030+0.0005+0.000139)/(0.003403) ≈ 1.07%。该结果属于示例计算，实际应结合样品长度加权与搜索半径限制。

插值必须设置搜索半径与样品数量上限，避免远距离样品“拉平”局部高品位；对矿体内部与围岩应分别插值，不能跨界线混算。插值完成后应检查品位等值线是否出现不合理的“牛眼”或锯齿，必要时调整幂次或改用普通克里金。

四、地质剖面图生成：把数据链落到图面

剖面图生成是把上述成果按统一坐标系落到图面：先定剖面线方位与里程，建立“水平距离—高程”剖面坐标系；再依次绘制地形线、地层界线、矿体界线、品位等值线或品位色块，最后叠加钻孔投影、样品位置与图例。生成后要做三项检查：剖面线与钻孔位置是否吻合；界线在相邻剖面间是否可对比；品位高值区是否与样品集中区一致。只有数据链闭合，剖面图才具备工程量统计与开采设计的使用价值。