1. 课程导入与靶区评价逻辑
矿产勘查靶区加权叠加评价,本质上是在一个统一的空间参照下,把多种来源、多种尺度、多种可信度的地质证据图层,按照勘查目标和成矿模型赋予不同权重,逐像元综合成一个可比较的"找矿有利度"分值,并据此圈定优先部署的勘查靶区。它要回答的不是"哪里有一张异常图",而是"在有限预算下,先去哪里、为什么先去那里"。因此,产出形式通常是一张连续的有利度栅格图、若干按阈值或自然断点划分的靶区分级面,以及配套的靶区属性表(面积、主要证据组合、平均有利度、建议工作手段)。决策意义在于把专家经验显式化、可复现化:权重和证据选择一旦写入工程,就能被审阅、被修改、被同行复核,而不是停留在某位地质师的手绘圈定上。
从地质认识到加权综合,完整工作流可以概括为五个环节。第一是成矿模型与找矿标志的建立:明确目标矿种与矿床类型,例如斑岩型铜矿或构造蚀变岩型金矿,据此列出控矿要素——赋矿地层与岩性、控矿断裂及其缓冲距离、岩浆岩体接触带、蚀变异常、化探元素组合异常、重磁与激电异常等。第二是证据图层的准备与规范化:把地质图、断层线、岩体边界、化探点位、物探栅格统一到同一投影坐标系,统一像元大小与范围,处理缺失值与边界。第三是证据图层转化为证据栅格:线状要素做距离缓冲或核密度,点状异常做插值或密度面,分类图层做重分类打分,使所有图层都变成同一量纲下的"有利度"面。第四是权重确定与加权叠加:用专家打分、层次分析或证据权法给出权重,逐像元求加权和。第五是靶区圈定、分级与复核:按分位数或已知矿点验证确定阈值,提取连片面,剔除面积过小或落在禁采区的图斑,形成分级靶区并回填属性。
QGIS 在这条链路中承担三类角色。其一是数据管理:建议把全部证据图层装入一个 GeoPackage,矢量与栅格同库管理,统一设置为项目坐标系(例如 CGCS2000 3 度带或矿区本地投影),并明确像元大小(如 25 m 或 50 m)与统一范围,避免后续叠加时出现错位或重采样失真。其二是空间分析:用"距离缓冲"把断裂线转为距断裂带距离面,用"核密度"或"插值"把化探点位转为连续异常面,用"重分类"把分类图层映射为 0–1 有利度,最后用栅格计算器按权重逐层相乘再相加。其三是制图与成果表达:对有利度栅格做分级配色,叠加靶区边界与编号,输出图件与属性表供立项论证使用。
需要强调的是,加权叠加不是"把图层加起来"这么简单,它的科学性取决于三个判断。第一,证据是否与成矿模型一致:如果目标为构造蚀变岩型金矿,那么断裂缓冲距离的合理区间应由已知矿点到断裂的统计距离确定,而不是随意取 500 m。第二,权重是否可解释:无论用层次分析法构造判断矩阵,还是用证据权法计算正负权重,都应记录权重来源与理由,便于复核。第三,结果是否可验证:应保留一部分已知矿点不参与权重确定,作为检验集,观察高有利度区是否覆盖这些矿点,以此判断模型是否过拟合或证据选择是否偏离。
下面给出一个可直接套用的靶区属性表字段设计,用于把叠加结果落成可汇报的成果表。字段只保留靶区评价真正需要的信息,示例数值为假设值,实际应以工程计算结果为准。
| 字段 | 含义 | 示例值 |
|---|---|---|
| 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万地质填图+土壤剖面 |