3. 矿体模型与块体模型构建

露天境界优化的本质，是在一个被离散化的三维空间里反复回答“哪些块该采、采到哪里停”的问题。因此，矿体模型与块体模型的质量，直接决定了优化结果的可靠性。本节围绕地质数据库整理、矿体圈定、块体模型建立、品位与岩性赋值、模型精度影响五个环节展开，重点讲清“数据如何变成可计算的块体”这一过程。

一、地质数据库整理：一切建模的起点

地质数据库通常由钻孔测量表（Collar）、测斜表（Survey）、化验表（Assay）和岩性表（Lithology）构成。整理的核心是建立统一的空间参照与编码体系：先明确坐标系与高程基准（例如采用矿区独立坐标系、高程为1985国家高程基准），再统一孔口坐标、测斜方位与倾角、样品起止深度和品位单位（如 Cu 用 %、Au 用 g/t）。

实际操作中要重点处理三类问题：一是测斜缺失，可用相邻测点线性插值补全，但需在数据库中标记为“插值”；二是样品区间重叠或超界，须按孔深排序后裁剪到有效孔段；三是品位低于检出限的记录，应按“检出限的一半”赋值而非直接置零，否则会系统性低估品位。整理完成后，应做一次孔深—样品深度一致性校验，确保每条化验记录都能落到唯一孔段上。

二、矿体圈定：从见矿段到三维实体

矿体圈定是把“一维孔段”扩展为“三维体”的过程。常用做法是：先按工业指标（如边界品位 Cu≥0.2%、最小可采厚度 2 m、夹石剔除厚度 2 m）在单孔上划分矿段，再在剖面和平面上按地质规律连接矿段，最后生成封闭的三维矿体实体（Wireframe）。

圈定中要区分“地质圈定”与“品位圈定”：前者依据岩性接触带、断层等地质界面，后者依据品位等值面。对于露天矿，建议以地质圈定为主、品位圈定为辅，避免把连续矿化体切得过碎。圈定完成后必须检查实体是否闭合、有无自相交，并统计矿体体积与平均品位，与地质储量报告对比，偏差一般应控制在 5%～10% 以内。

三、块体模型建立流程

块体模型是把矿体空间切分为规则长方体单元（Block）的离散模型。建立流程一般为：确定模型范围（覆盖矿体并留出 1～2 个块的外扩边界）→ 设定块尺寸→ 建立母块与子块（Sub-blocking）→ 用矿体实体约束块体（实体内的块标记为矿，实体外标记为废石）→ 逐块赋值品位、岩性、密度。

块尺寸的选择是核心权衡。以某铜矿示例：矿体走向长约 800 m、倾向延深 400 m，若采用 10 m×10 m×10 m 块，模型约 3.2 万块，边界较粗糙；若采用 5 m×5 m×5 m，块数增至约 25.6 万块，边界贴合度明显提高，但估值与优化计算量约增加 8 倍。工程上常用“母块 10 m + 子块 2.5 m”的混合方案，在矿岩接触带自动细分，兼顾精度与效率。

四、品位与岩性赋值方法

品位赋值常用普通克里格（OK）或距离反比加权（IDW）。以 OK 为例，需先做变异函数拟合：沿走向、倾向、厚度方向分别拟合球状模型，得到变程（如走向 120 m、倾向 80 m、厚度 25 m），据此确定搜索椭球的方向与大小。估值时限制每个块搜索到的样品数（如 4～12 个），并约束单孔样品占比不超过 50%，防止一条高品位孔“污染”大片块体。

岩性赋值相对简单，通常用“实体归属法”：块中心落在某岩性实体内部即赋该岩性；对跨界面块，可按体积占比或最近样品法处理。密度则按岩性查表赋值，例如斑岩 2.65 t/m³、围岩 2.70 t/m³、氧化带 2.30 t/m³。赋值后应做一次“块体平均品位 vs 组合样平均品位”的对比校验，偏差超过 10% 需检查搜索参数或变异函数。

五、模型精度对境界优化的影响

块体模型精度通过三条路径影响最终境界：其一，品位估值偏差改变块价值，高估会使境界外扩、低估会使境界收缩；其二，块尺寸过大会把矿岩界面“抹平”，导致境界在接触带附近出现系统性偏差；其三，岩性与密度赋值错误会直接影响剥采比与成本核算。

实践中可用“块尺寸敏感性分析”检验稳健性：分别用 10 m、5 m、2.5 m 三种块尺寸建模并优化，若最终境界体积差异小于 3%，说明模型精度足够；若差异超过 8%，则应缩小块尺寸或加密采样。此外，建议保留一份“未估值原始块体模型”，便于在优化结果异常时回溯是数据问题还是参数问题。