4. 矿体地质模型与域划分

在资源量估算中,矿体地质模型与域划分(Domaining)是连接“地质认识”与“数值估算”的桥梁。域划分的本质,是把连续、渐变、非均质的矿化体,按某种可解释的地质或品位规则切分为若干内部相对均一、边界可追踪的单元,使每个单元内的样品统计特征、变异函数结构和估算参数尽量一致。域划错了,后续无论采用普通克里格、指示克里格还是距离反比,都会把系统偏差带入块体模型。

一、矿化域圈定的三条常用路径

实际工作中,矿化域圈定通常有三条路径,可单独使用也可组合使用。

1. 地质边界圈定。以岩性接触面、蚀变带界线、构造破碎带、氧化—原生分带界面等地质标志为依据,在剖面和平面图上手工或半自动连线圈定。其优点是地质意义明确、可解释性强,适合矿化明显受岩性/构造控制的矿床(如矽卡岩型、脉状矿床)。缺点是边界位置依赖地质人员判断,主观性较大,且对渐变型矿化不敏感。

2. 品位边界圈定。以边界品位(cut-off grade)或某一指示阈值为准,用样品品位在剖面上圈定等值线或包络线。常见做法是先按工业指标确定边界品位(例如 Cu 0.2%),再结合最小可采厚度、夹石剔除厚度等指标圈定矿体。优点是客观、可复现;缺点是品位边界受样品分布和块段尺寸影响大,容易产生“锯齿状”边界,且高品位样品会显著外推边界。

3. 统计/地质统计学方法。利用品位累积概率曲线拐点、指示克里格、序贯指示模拟或聚类分析(如 k-means 对多元素样品分组)辅助确定域边界。这类方法适合矿化连续渐变、地质界线不清的情形,但需要足够样品量,且结果需地质解释校验。

工程实践中,推荐“地质定框架、品位定细节、统计做校验”的组合策略:先用地质疑定大域(如氧化带/原生带、不同岩性单元),再在域内用品位或统计方法细分高/低品位亚域,最后用变异函数和交叉验证检验域内均一性。

二、地质边界与品位边界的对比

两种边界并非对立,而是回答不同问题。地质边界回答“矿化为什么在这里”,品位边界回答“矿化到什么程度算矿”。下表从几个关键维度对比。

对比维度地质边界品位边界
圈定依据岩性、蚀变、构造、分带等地质标志边界品位、指示阈值、统计拐点
客观性/可复现性依赖地质判断,主观性较强规则明确,可复现性好
对渐变矿化的适应性较弱,易“一刀切”较强,可反映品位过渡
对样品稀疏区的稳健性较好,可借助地质外推较差,边界易被个别高值样品拉动
与估算方法的匹配适合地质统计学建模、地质体约束适合指示克里格、品位阈值建模
主要风险边界主观、域内可能仍非均质边界锯齿、高值外推、域内样品不足

需要强调的是,边界品位本身是一个经济—技术参数,会随金属价格、采选成本、回收率变化而调整。因此品位边界圈定的矿体是“动态”的,而地质边界相对稳定。在方案对比中,应明确所采用的边界品位取值及其依据,避免不同方案因边界品位不一致而失去可比性。

三、域划分对估算方案的影响

域划分直接决定样品如何分组、变异函数如何拟合、搜索椭球如何设置、块体如何赋值。其影响主要体现在四个方面。

(1)样品分组与统计参数。域划分后,每个域内样品独立统计。若域划得过粗,域内可能混合了高品位与低品位样品,品位分布出现双峰,变异函数呈现“块金效应大、变程短”的假象;若域划得过细,单域样品数不足,变异函数拟合不稳定,甚至无法拟合。经验上,单域有效样品数宜不少于 30~50 个,否则应合并相邻域或借用区域化变量结构。

(2)变异函数结构。不同域的矿化连续性往往不同。例如某斑岩铜矿示例(教学假设):氧化带域变程约 60 m,原生带域变程约 120 m。若不分域统一拟合,会得到一个“折中”的变程,导致氧化带块体估值过于平滑、原生带估值过于保守。分域拟合后,各域采用各自变程与各向异性方向,估值更贴合地质实际。

(3)搜索椭球与边界处理。域边界是硬边界还是软边界,直接影响搜索策略。硬边界下,搜索椭球不跨越域界,避免相邻域样品“污染”;软边界下允许跨域但降权。对于品位突变明显的接触带,宜采用硬边界;对于渐变过渡带,可采用软边界或设置过渡域。

(4)资源量分类与不确定性。域边界附近往往是估值不确定性最高的区域。域划分越细、边界越复杂,边界附近的块体越难获得足够样品支持, measured/indicated/inferred 的转换关系也随之变化。方案对比时,应把“域划分方案”作为一个独立变量,固定其他参数(块体尺寸、搜索参数、边界品位),观察资源量总量与分类比例的变化。

域划分 → 估算参数 → 资源量结果 的影响链 域划分方案 地质边界 / 品位边界 硬边界 / 软边界 域数量与粒度 域内估算参数 样品分组与统计 均值 / 方差 / 分布形态 变异函数拟合 块金 / 变程 / 各向异性 搜索椭球设置 跨域策略 / 样品数下限 块体赋值与分类 估值 / 分类 / 边界处理 资源量结果 总量 / 平均品位 分类比例 局部估值偏差 交叉验证 / 边界校验反馈 示例对比(教学假设,非实测):同一矿床两种域划分方案 方案A:仅按地质边界分2域 域内品位双峰,变程折中,边界附近估值偏平滑 方案B:地质域+品位亚域共4域 各域变程独立,分类比例更合理,边界需样品支撑

四、方案对比时的操作要点

若要把域划分作为方案对比变量,建议按以下顺序推进,确保对比有效:

  1. 固定其他变量。块体尺寸、边界品位、搜索参数、估算方法(如均用普通克里格)保持一致,只改变域划分方案。
  2. 建立域划分记录表。记录每个域的圈定依据、边界类型(硬/软)、样品数、品位均值与方差、变异函数参数(块金、变程、各向异性方位)。
  3. 做域内均一性检验。对每个域做品位直方图与累积概率曲线,若出现明显双峰或长尾,说明域内仍非均质,需细分或调整边界。
  4. 交叉验证。采用留一法或 k 折交叉验证,比较不同域方案下估值误差的均值(偏差)与均方根误差(精度),并重点检查域边界附近块体的误差。
  5. 结果对比与决策。对比资源量总量、平均品位、各分类级别比例,结合地质合理性判断哪一方案更可信,而非单纯追求误差最小。
常见误区:为追求“边界好看”而过度细分域,导致单域样品不足、变异函数无法拟合;或为省事只划一个大域,掩盖了品位分布的双峰结构。两种极端都会使资源量分类失去意义。域划分应服务于“域内均一、域间可区分”这一核心目标。

综上,矿化域圈定不是一次性制图工作,而是贯穿建模、估值、验证的迭代过程。地质边界提供可解释的框架,品位边界提供可量化的约束,二者结合并辅以统计校验,才能得到既符合地质认识、又满足估算精度要求的域模型。在方案对比中,应把域划分视为与估算方法同等重要的变量,用统一的对比框架和验证指标来评估其影响。