2. 遥感找矿地质基础

遥感找矿的核心逻辑,是把“肉眼看不见的矿物成分变化”转化为“传感器可量测的光谱差异”。热液活动在围岩中形成的蚀变带,往往会生成一批在可见光—短波红外(VNIR–SWIR,约0.4~2.5 μm)区间具有诊断性吸收特征的矿物。这些矿物在特定波长处吸收光子,使反射率曲线出现“凹陷”,而凹陷的位置、深度和宽度,正是识别矿物种类与相对丰度的依据。理解这套“矿物—波长—吸收特征”的对应关系,是后续选择波段、构造指数、设定阈值的前提。

一、三类关键蚀变矿物的波谱特征

在斑岩—浅成低温热液—矽卡岩等典型矿床的蚀变分带中,铁染矿物、羟基矿物和碳酸盐矿物构成了遥感蚀变提取的三大主线,它们分别对应不同的诊断吸收位置。

铁染矿物包括褐铁矿、赤铁矿、针铁矿、黄钾铁矾等。其电子跃迁过程在可见光区产生强吸收:Fe³⁺的晶体场吸收通常落在0.85~0.95 μm附近,并伴随0.48~0.52 μm的蓝紫光吸收,导致富铁岩石在真彩色影像上呈褐红、橙黄色调。因此铁染异常常用波段比值(如红/蓝、或近红外/红)增强,表现为高比值亮异常。

羟基矿物包括高岭石、伊利石、绢云母、明矾石、绿泥石、蒙脱石等。Al–OH的合频与倍频吸收集中在2.17~2.22 μm,Mg–OH则偏向2.30~2.35 μm。高岭石双吸收在2.16 μm与2.21 μm,明矾石在2.17 μm附近,绢云母(白云母族)在2.20 μm附近,绿泥石除2.25 μm外还在2.33~2.35 μm有特征。这类矿物是绢英岩化、泥化、青磐岩化的标志,也是SWIR提取的主力对象。

碳酸盐矿物包括方解石、白云石、菱铁矿等。CO₃²⁻的合频吸收位于2.31~2.35 μm,方解石约2.34 μm、白云石约2.32 μm,与绿泥石的Mg–OH位置接近,因此实际提取中需结合2.20 μm附近有无Al–OH吸收加以区分,避免“碳酸盐—绿泥石”混淆。

蚀变矿物诊断吸收位置与矿床蚀变分带对应关系 λ/μm 0.5 0.9 1.4 1.9 2.2 2.35 2.5 铁染矿物 Fe³⁺:0.48~0.52 / 0.85~0.95 μm 羟基矿物 Al–OH:2.17~2.22 μm Mg–OH:2.30~2.35 μm 碳酸盐矿物 CO₃²⁻:2.31~2.35 μm 斑岩型矿床蚀变分带(由内向外) 钾化带 黑云母/钾长石 绢英岩化带 绢云母/石英 泥化带 高岭石/明矾石 青磐岩化带 绿泥石/绿帘石/方解石 围岩 未蚀变 判读要点: ① 2.20 μm 附近 Al–OH 吸收 → 绢云母/高岭石/明矾石,指示绢英岩化与泥化; ② 2.32~2.35 μm 吸收需结合 2.20 μm 有无吸收区分碳酸盐与绿泥石; ③ 铁染异常(0.9 μm 附近)常为氧化带/黄铁矿风化帽标志,与羟基异常叠合处找矿意义更强。

二、蚀变分带与典型矿床的对应关系

蚀变矿物不是孤立出现的,它们按温度、酸碱度和流体成分呈规律分带。以斑岩铜矿为例,自岩体中心向外依次为钾化带(黑云母、钾长石)、绢英岩化带(绢云母、石英、黄铁矿)、泥化带(高岭石、明矾石、叶蜡石)和青磐岩化带(绿泥石、绿帘石、方解石)。遥感SWIR数据能同时捕捉2.20 μm的绢云母/高岭石异常与2.33 μm的绿泥石/碳酸盐异常,从而在影像上勾出“内Al–OH、外Mg–OH/CO₃”的环状分带,这正是斑岩系统最典型的找矿遥感标志。

浅成低温热液金矿则以明矾石—高岭石—伊利石组合的高级泥化带为特征,明矾石在2.17 μm的强吸收可作为酸性蚀变中心的定位依据;矽卡岩型矿床外围常见方解石、白云石大理岩化,2.34 μm碳酸盐异常可辅助圈定接触带。可见,同一套SWIR波段数据,通过不同吸收位置的组合判读,可以服务于不同矿种、不同尺度的找矿目标。

三、蚀变异常作为找矿标志的地质意义

蚀变异常之所以能作为找矿标志,根本原因在于热液蚀变是成矿流体活动的“化石记录”:流体带来成矿物质的同时,也改造了围岩的矿物组合,蚀变晕的范围通常比矿体大数倍至数十倍,因而成为有效的“放大靶区”。遥感提取的蚀变异常,本质上是蚀变矿物相对富集区的空间表达,其地质意义体现在三个层面:一是指示流体通道与热源位置,蚀变强度高值区往往靠近岩体或断裂交汇处;二是反映分带结构,环状或带状异常组合可推断剥蚀深度与矿化中心;三是约束矿化类型,Al–OH与Fe³⁺异常叠合、伴生高级泥化,指向浅成低温或斑岩型金铜矿。

需要强调:遥感蚀变异常是“矿物学异常”而非“元素异常”,它回答的是“哪里发生了热液蚀变”,而不是“哪里一定有矿”。异常必须结合地质构造、物探、化探与野外查证综合解释,才能转化为找矿靶区。

四、从波谱特征到提取方法的衔接

明确了诊断吸收位置,后续提取方法的选择就有了依据。对于铁染,可选用波段比值或主成分分析(如对TM/OLI的2、4、5波段做特征向量判读);对于羟基与碳酸盐,则依赖SWIR波段(ASTER的5、6、7、8、9波段,或Sentinel-2的11、12波段)构造相对吸收深度指数,如用(B6+B8)/B7增强2.20 μm吸收、用(B7+B9)/B8增强2.33 μm吸收。阈值设定时,应结合影像直方图与已知矿点统计,采用均值加n倍标准差或分形方法分级,避免单一固定阈值造成漏提或误提。

下表给出本节可直接参考的矿物—波段—指数对应关系示例(波段编号以ASTER SWIR为例,参数为教学示例,实际应用需按数据定标结果调整):

蚀变类型诊断吸收ASTER波段常用增强方式地质指示
铁染0.85~0.95 μmB2/B1、B4/B3比值、主成分分析氧化带、黄铁矿风化帽
Al–OH(绢云母/高岭石)2.17~2.22 μmB5、B6(B5+B7)/B6绢英岩化、泥化带
Mg–OH(绿泥石)2.30~2.35 μmB8、B9(B7+B9)/B8青磐岩化带
碳酸盐2.31~2.35 μmB8、B9(B6+B9)/B8大理岩化、接触带

把这些对应关系记牢,后续无论是做波段选择、指数构造还是异常分级,都能做到“知其所以然”,而不是机械套用软件默认参数。