1. 课程导入与靶区评价逻辑

矿产勘查靶区加权叠加评价,本质上是在一个统一的空间参照下,把多种来源、多种尺度、多种可信度的地质证据图层,按照勘查目标和成矿模型赋予不同权重,逐像元综合成一个可比较的"找矿有利度"分值,并据此圈定优先部署的勘查靶区。它要回答的不是"哪里有一张异常图",而是"在有限预算下,先去哪里、为什么先去那里"。因此,产出形式通常是一张连续的有利度栅格图、若干按阈值或自然断点划分的靶区分级面,以及配套的靶区属性表(面积、主要证据组合、平均有利度、建议工作手段)。决策意义在于把专家经验显式化、可复现化:权重和证据选择一旦写入工程,就能被审阅、被修改、被同行复核,而不是停留在某位地质师的手绘圈定上。

从地质认识到加权综合,完整工作流可以概括为五个环节。第一是成矿模型与找矿标志的建立:明确目标矿种与矿床类型,例如斑岩型铜矿或构造蚀变岩型金矿,据此列出控矿要素——赋矿地层与岩性、控矿断裂及其缓冲距离、岩浆岩体接触带、蚀变异常、化探元素组合异常、重磁与激电异常等。第二是证据图层的准备与规范化:把地质图、断层线、岩体边界、化探点位、物探栅格统一到同一投影坐标系,统一像元大小与范围,处理缺失值与边界。第三是证据图层转化为证据栅格:线状要素做距离缓冲或核密度,点状异常做插值或密度面,分类图层做重分类打分,使所有图层都变成同一量纲下的"有利度"面。第四是权重确定与加权叠加:用专家打分、层次分析或证据权法给出权重,逐像元求加权和。第五是靶区圈定、分级与复核:按分位数或已知矿点验证确定阈值,提取连片面,剔除面积过小或落在禁采区的图斑,形成分级靶区并回填属性。

从地质认识到加权综合:证据流与QGIS承担角色 ① 成矿模型 矿种·类型·控矿要素 找矿标志清单 ② 证据图层 统一投影/像元/范围 地质·断裂·化探·物探 ③ 证据栅格化 缓冲/密度/插值/重分类 统一为有利度 0–1 ④ 加权叠加 AHP/专家打分/证据权 逐像元加权求和 ⑤ 靶区圈定 阈值分级·连片 属性回填·复核 证据图层示例(同一范围、同一像元) 赋矿地层 分类→重分类 控矿断裂 距离缓冲 岩体接触带 距离缓冲 化探组合异常 插值/密度面 激电/磁异常 栅格重分类 遥感蚀变 主成分/比值 已知矿点 验证用 加权综合(逐像元) S = Σ ( wᵢ × rᵢ ) ,Σwᵢ = 1 wᵢ:证据权重;rᵢ:该像元在第 i 层的有利度(0–1) 权重示例(假设值) 断裂 0.30 · 化探 0.25 · 地层 0.20 岩体 0.15 · 物探 0.10 QGIS 承担的角色 数据管理:GeoPackage 统一存放矢量/栅格,投影与像元对齐,图层样式与元数据 空间分析:距离缓冲、核密度、栅格重分类、栅格计算器加权求和、分区统计 制图输出:分级配色、图例、靶区编号标注 成果:有利度图 + 靶区面 + 属性表

QGIS 在这条链路中承担三类角色。其一是数据管理:建议把全部证据图层装入一个 GeoPackage,矢量与栅格同库管理,统一设置为项目坐标系(例如 CGCS2000 3 度带或矿区本地投影),并明确像元大小(如 25 m 或 50 m)与统一范围,避免后续叠加时出现错位或重采样失真。其二是空间分析:用"距离缓冲"把断裂线转为距断裂带距离面,用"核密度"或"插值"把化探点位转为连续异常面,用"重分类"把分类图层映射为 0–1 有利度,最后用栅格计算器按权重逐层相乘再相加。其三是制图与成果表达:对有利度栅格做分级配色,叠加靶区边界与编号,输出图件与属性表供立项论证使用。

需要强调的是,加权叠加不是"把图层加起来"这么简单,它的科学性取决于三个判断。第一,证据是否与成矿模型一致:如果目标为构造蚀变岩型金矿,那么断裂缓冲距离的合理区间应由已知矿点到断裂的统计距离确定,而不是随意取 500 m。第二,权重是否可解释:无论用层次分析法构造判断矩阵,还是用证据权法计算正负权重,都应记录权重来源与理由,便于复核。第三,结果是否可验证:应保留一部分已知矿点不参与权重确定,作为检验集,观察高有利度区是否覆盖这些矿点,以此判断模型是否过拟合或证据选择是否偏离。

教学建议:第一次练习可先用 4–5 个证据图层、固定权重跑通全流程,重点体会"统一像元—栅格化—加权求和—阈值分级"四步的衔接,再逐步引入权重计算与验证集检验。

下面给出一个可直接套用的靶区属性表字段设计,用于把叠加结果落成可汇报的成果表。字段只保留靶区评价真正需要的信息,示例数值为假设值,实际应以工程计算结果为准。

字段含义示例值
target_id靶区编号A-01
area_km2靶区面积(km²)3.6
mean_score平均有利度(0–1)0.78
max_score最高有利度0.91
key_evidence主要证据组合断裂缓冲+Cu异常+蚀变
grade分级(Ⅰ/Ⅱ/Ⅲ)Ⅰ
suggest_work建议工作手段1:1万地质填图+土壤剖面