4. 钻孔测斜数据复核与轨迹重建

在矿产资源估算中，钻孔见矿数据通常以“孔口坐标 + 测点深度 + 测斜数据”的形式进入三维模型。若直接假定钻孔为铅直孔，在深孔或陡倾斜地层中会带来数米到数十米的平面偏移，进而影响矿体边界、品位插值和资源量分类。本节围绕测斜方法、偏斜计算、轨迹插值、轨迹与剖面交切关系检查及异常处理展开，目标是让读者能够独立完成一条钻孔的测斜复核与轨迹重建，并判断其是否可用于剖面复核。

一、测斜方法与数据形式

常用测斜手段包括单点/多点陀螺测斜、磁通门电子测斜仪以及随钻测量（MWD）。无论采用哪种仪器，最终记录的核心参数都是测点深度、倾角（Dip，相对水平面，向下为正）和方位角（Azimuth，自正北顺时针量取）。现场记录通常还包含磁偏角修正、仪器编号、测量日期与质量标记。复核时首先要确认三件事：方位角是否已归算到真北；倾角符号约定是否统一（向下为正还是向上为正）；测点深度是井深（沿孔轴）还是垂深。这三项约定一旦混用，后续轨迹计算会系统性偏移。

一个常见的教学示例：某孔孔口坐标 X=500000 m、Y=4000000 m、Z=1200 m，设计为斜孔，开孔方位 135°、倾角 -60°（即向下 60°）。实测测斜数据如下表（示例数据，用于说明计算方法，非实测成果）：

二、偏斜计算与轨迹插值

轨迹重建的本质是把离散测点之间的孔段近似为一段圆弧或直线，逐段累加得到各测点的三维坐标。工程上最常用的是最小曲率法：把相邻两测点之间的孔段视为一段圆弧，其曲率由两点倾角与方位角之差确定。设上测点倾角 θ₁、方位角 φ₁，下测点 θ₂、φ₂，孔段长度 ΔL，则全角变化量（狗腿角）β 满足：

cos β = cos(θ₂ − θ₁) − sin θ₁ · sin θ₂ · [1 − cos(φ₂ − φ₁)]

当 β 很小时用直孔段近似，避免除零。孔段在三个方向上的增量按比例因子 RF = (2/β)·tan(β/2) 修正：

ΔZ = (ΔL/2)(cos θ₁ + cos θ₂)·RF
ΔN = (ΔL/2)(sin θ₁cos φ₁ + sin θ₂cos φ₂)·RF
ΔE = (ΔL/2)(sin θ₁sin φ₁ + sin θ₂sin φ₂)·RF

以第 0→1 段为例：ΔL=50 m，θ₁=60°、θ₂=59.2°（取向下为正），φ₁=135°、φ₂=136.5°。代入得 β≈1.53°，RF≈1.00006。计算得 ΔZ≈−25.0 m、ΔN≈−17.6 m、ΔE≈+17.7 m。逐段累加后，200 m 处测点相对孔口的偏移约为 N 方向 −70 m、E 方向 +72 m，垂深约 −98 m。若按铅直孔处理，200 m 处将错误地落在孔口正下方 200 m 处，平面偏差超过 100 m——这正是必须做测斜复核的原因。

插值密度应与矿体厚度和剖面间距匹配。一般要求轨迹折点在剖面附近有测点控制，测点间距不宜大于 50 m；在见矿段附近应加密到 10～25 m，避免用长孔段圆弧跨越矿体边界造成“假交切”。

三、轨迹与剖面交切关系检查

剖面复核的核心是判断钻孔轨迹与剖面面的交点位置是否落在见矿区间内。设剖面走向方位为 α，剖面线过点 P₀，则剖面法向量 n=(sin α, −cos α, 0)（水平面内）。对轨迹上每个插值点 P，计算有符号距离 d = n·(P−P₀)。当相邻两插值点的 d 异号时，说明轨迹穿过剖面，用线性插值求出交点深度与坐标，再与采样表对照：交点深度应落在见矿段起止深度之间，且交点平面坐标与剖面图上标注的见矿位置一致。若交点落在矿段之外，或同一矿段被轨迹多次穿越，说明测斜数据或剖面方位存在问题，需要回溯。

上图左侧为平面投影，对比铅直假设与实测轨迹的偏移；右侧为沿剖面线展开的垂直剖面，标出见矿段与轨迹交切点。复核时应确认交切点落在红色见矿段内，且交点深度与采样表一致。

四、异常识别与处理

测斜复核中常见的异常有四类。第一类是狗腿度突变：相邻测点全角变化超过 3°/30 m 时，需检查是否为仪器故障或记录错误，必要时补测。第二类是方位角跳变：在磁性地层中磁通门仪器易受干扰，表现为方位角突然偏转数十度，应结合陀螺复测或区域磁偏角修正判断。第三类是测点缺失：深部无测斜数据时，不得简单外推，应按最后测点趋势延伸并标注“推测段”，在资源量分类中降级处理。第四类是轨迹与剖面交切异常：交点落在矿段外或同一矿段多次穿越，通常提示剖面方位错误或测斜数据未归算到真北。

五、复核成果与判断依据

一条钻孔完成复核后，应输出：归算到统一坐标系的测点表、逐段轨迹增量与累计坐标、狗腿度统计、剖面交切点深度与坐标、异常标记及处理说明。判断轨迹是否可用于剖面复核的依据是：测点间距满足要求、狗腿度在允许范围内、交切点与见矿段一致、平面偏移量已量化并记录。只有满足这些条件，钻孔见矿数据才能进入后续的剖面圈连与资源量估算环节。