1. 课程导入与找矿预测问题界定
找矿预测的本质,是在一个面积广大、露头有限、勘查成本高昂的研究区内,回答一个空间决策问题:哪些位置最有可能赋存目标矿种,值得优先投入野外验证与工程揭露。地质、物探、化探、遥感等手段提供了大量与成矿有关的证据图层,但这些证据的度量单位彼此不可比——断裂缓冲距离以米计,化探异常以元素含量计,地层含矿性以类别代号计,遥感蚀变以指数值计。若直接把原始图层叠加,数值大的图层会淹没数值小的图层,叠加结果没有地质意义。
栅格重分类(Reclassify)正是解决这一问题的关键环节。它把不同量纲、不同连续性的证据图层,统一映射到同一套离散等级体系(例如 1~5 级或 0~9 的成矿有利度分值),使"距离断裂 300 m 以内"与"Cu 元素含量处于高背景区"可以在同一尺度上比较与加权。可以说,重分类是连接"原始地质数据"与"可计算预测模型"的转换器,也是整个找矿预测流程中主观地质认识注入最集中、对结果影响最大的步骤之一。
从地质数据到预测靶区的完整工作流
本课程采用的实战工作流可概括为五个阶段,各阶段的输入、处理与输出如下表所示(示例参数用于教学演示,实际项目需依据研究区地质特征调整):
| 阶段 | 主要输入 | 核心处理 | 输出 |
|---|---|---|---|
| 数据准备 | 地质图、断裂线、DEM、化探点位、遥感影像 | 坐标系统一、裁剪、栅格化、插值、缓冲区 | 统一投影与分辨率的证据栅格集 |
| 证据构建 | 证据栅格集 | 距离分析、密度分析、异常分级、蚀变指数计算 | 连续型证据图层 |
| 重分类 | 连续型证据图层 | 按地质意义设定分级阈值并赋有利度分值 | 离散等级证据图层 |
| 加权合成 | 离散等级证据图层 | 权重确定(专家打分/信息量/证据权)、栅格计算器叠加 | 成矿有利度连续面 |
| 靶区圈定 | 成矿有利度面 | 阈值分割、连通性聚合、面积与等级统计 | 分级预测靶区与优先排序 |
需要强调的是,重分类并非孤立操作。它上游依赖证据构建阶段产出的图层质量,下游决定加权合成的可解释性。若重分类阈值设定缺乏地质依据,后续无论采用多精细的权重算法,输出都只是"数学上合理、地质上无意义"的图案。
ArcGIS 平台上的工具落点
在 ArcGIS Pro 中,上述工作流主要落在 Spatial Analyst 与 Conversion Tools 两个工具箱内。与本课程直接相关的工具包括:欧氏距离(Euclidean Distance)用于生成断裂缓冲距离栅格;IDW 或克里金插值用于把化探点位转为连续异常面;重分类(Reclassify)与查找表重分类(Reclass by Table)用于等级映射;栅格计算器(Raster Calculator)用于加权叠加;重分类为整数与栅格转面用于靶区边界提取。ArcMap 用户可在 ArcToolbox 中找到同名工具,参数界面略有差异但逻辑一致。
实战数据集围绕一个假想的斑岩铜矿研究区组织,包含:1:5 万地质图(地层与岩体面要素)、断裂线要素、30 m 分辨率 DEM、1:2.5 万土壤化探点位(Cu、Mo、Au 三元素)、Landsat 8 影像衍生的铁染与羟基蚀变指数栅格。所有图层预先统一为 CGCS2000 高斯-克吕格投影,栅格分辨率统一重采样至 30 m,确保后续重分类与叠加时像元严格对齐——这是很多初学者容易忽略却直接导致结果错位的前提条件。
本节要界定的核心问题
围绕"重分类如何影响找矿预测结果",本课程将逐层展开三个问题:其一,分级依据从何而来——阈值不能凭软件默认的分位数机械切分,而应结合区域成矿模式、已知矿点分布与证据图层统计特征共同确定;其二,不同证据如何加权——专家打分法直观但主观性强,信息量法与证据权法可借助已知矿点做统计检验,二者适用条件不同;其三,结果如何验证与表达——预测靶区需要给出等级、面积、包含已知矿点的比例等可核查指标,而非仅一张彩色图。