1. 课程导入与ASTER数据认知

ASTER(Advanced Spaceborne Thermal Emission and Reflection Radiometer,先进星载热发射与反射辐射计)搭载于1999年发射的Terra卫星,是矿产遥感领域应用最广的多光谱数据源之一。它之所以在热液蚀变矿物识别中占据重要地位,根本原因在于其波段设置恰好覆盖了蚀变矿物在可见光—近红外(VNIR)、短波红外(SWIR)和热红外(TIR)三个区间的诊断性吸收与发射特征。理解这套波段"为什么这样切",是后续所有波段组合分析的前提。

一、ASTER的波段设置与矿物诊断意义

ASTER共14个波段,分为三个独立光学子系统:VNIR 3个波段(0.52–0.86 μm,空间分辨率15 m)、SWIR 6个波段(1.60–2.43 μm,30 m)、TIR 5个波段(8.125–11.65 μm,90 m)。其中SWIR的6个波段是热液蚀变识别的核心:B4(1.60–1.70 μm)、B5(2.145–2.185 μm)、B6(2.185–2.225 μm)、B7(2.235–2.275 μm)、B8(2.295–2.365 μm)、B9(2.360–2.430 μm)。这些窄波段分别对应Al-OH、Mg-OH、Fe-OH和碳酸盐的吸收位置差异,使得不同蚀变矿物在波段间呈现可区分的吸收深度组合。

具体而言,明矾石、高岭石等Al-OH矿物在B5/B6附近出现吸收;绿泥石、绿帘石等Mg-OH、Fe-OH矿物吸收位置偏向B7/B8;方解石等碳酸盐在B8附近有特征吸收。TIR波段则对石英、长石、碳酸盐等造岩矿物的发射率差异敏感,可辅助区分硅化带与围岩。VNIR的B1、B2、B3(含3B立体后视)主要用于铁氧化物识别(赤铁矿、褐铁矿在B1/B2吸收)以及影像配准与地形校正。

ASTER 14波段与热液蚀变矿物诊断特征对应 VNIR · 15 m B1 0.52–0.60 μm B2 0.63–0.69 μm B3 0.76–0.86 μm(含3B) → 铁氧化物、配准/地形 SWIR · 30 m(核心) B4 1.60–1.70 B5 2.145–2.185 B6 2.185–2.225 B7 2.235–2.275 B8 2.295–2.365 B9 2.360–2.430 → Al-OH / Mg-OH / Fe-OH / CO₃ TIR · 90 m B10–B14 8.125–11.65 μm 发射率光谱 → 石英、碳酸盐、硅化带 Al-OH 类 明矾石·高岭石·绢云母 吸收位:B5/B6 Mg-OH / Fe-OH 类 绿泥石·绿帘石·阳起石 吸收位:B7/B8 碳酸盐类 方解石·白云石 吸收位:B8 铁氧化物 赤铁矿·褐铁矿 吸收位:B1/B2 波段组合分析定位:波段比值/相对吸收深度 → 蚀变矿物填图 → 蚀变分带 → 找矿靶区圈定 输入:ASTER L1T/L2 反射率产品 输出:矿物丰度/分布图与勘查靶区

二、热液蚀变矿物识别的勘查意义

热液蚀变是含矿流体与围岩相互作用的产物,蚀变矿物组合与分带直接记录了流体的温度、酸碱度和成分演化。斑岩—浅成低温热液系统中,从中心向外常出现硅化—绢英岩化(绢云母+石英)—泥化(高岭石、明矾石)—青磐岩化(绿泥石+绿帘石+方解石)的分带格局,而矿体往往与特定蚀变带空间耦合。因此,识别蚀变矿物种类及其空间分带,本质上是在"读取"流体通道与沉淀环境,是缩小找矿范围最经济有效的手段之一。ASTER的SWIR波段能在区域尺度上快速圈定这类分带,成本远低于大面积地面填图。

三、波段组合分析在矿产勘查中的定位

需要明确:波段组合分析不是孤立的技术动作,而是"数据预处理—波段组合/指数计算—异常提取—地质解译—野外验证"链条中的核心分析环节。它的上游依赖辐射定标与大气校正得到的表观反射率,下游输出的是矿物异常图层,最终要与地质图、构造线、地球化学异常叠加才能形成靶区。常见组合包括比值法(如B5/B6增强Al-OH吸收、B7/B6区分Mg-OH与Al-OH)、相对吸收深度指数(如B6相对B5与B7连线的凹陷深度)以及主成分分析(如Crosta法对B1、B3、B4、B6做特征向量定向)。这些方法的选择依据,正是本节强调的波段—矿物对应关系。

学习路径建议:先掌握波段与矿物的物理对应(本节),再依次进入数据预处理、比值与指数计算、主成分分析、蚀变分带解译与靶区圈定,最后通过已知矿区案例完成从影像到靶区的完整闭环。每一步都应回到"这个波段组合在增强哪种矿物的哪个诊断特征"这一根本问题上。

四、一个可操作的入门示例

以某斑岩铜矿区ASTER L1T产品为例(示例参数,非实测):完成大气校正后,取B5、B6、B7三个波段计算绢云母指数 B5/B6 与绿泥石指数 B7/B6,再计算B6的相对吸收深度 D6 = 1 − B6 / (a·B5 + b·B7)(a、b为B5、B7到B6的线性插值权重)。将D6高值区与B5/B6高值区叠加,若二者空间重合且呈环带状围绕侵入体接触带分布,则指示绢英岩化带;其外侧B7/B6高值区对应青磐岩化带。此时应导出异常图层,与1:5万地质图套合,标记接触带附近的异常集中区作为野外查证优先区,并用便携式光谱仪对典型露头采样验证矿物种类。这一闭环把波段组合从"算出图"推进到"能判断、可验证"。