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吸收加以区分,避免“碳酸盐—绿泥石”混淆。
二、蚀变分带与典型矿床的对应关系
蚀变矿物不是孤立出现的,它们按温度、酸碱度和流体成分呈规律分带。以斑岩铜矿为例,自岩体中心向外依次为钾化带(黑云母、钾长石)、绢英岩化带(绢云母、石英、黄铁矿)、泥化带(高岭石、明矾石、叶蜡石)和青磐岩化带(绿泥石、绿帘石、方解石)。遥感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 μm | B2/B1、B4/B3 | 比值、主成分分析 | 氧化带、黄铁矿风化帽 |
| Al–OH(绢云母/高岭石) | 2.17~2.22 μm | B5、B6 | (B5+B7)/B6 | 绢英岩化、泥化带 |
| Mg–OH(绿泥石) | 2.30~2.35 μm | B8、B9 | (B7+B9)/B8 | 青磐岩化带 |
| 碳酸盐 | 2.31~2.35 μm | B8、B9 | (B6+B9)/B8 | 大理岩化、接触带 |
把这些对应关系记牢,后续无论是做波段选择、指数构造还是异常分级,都能做到“知其所以然”,而不是机械套用软件默认参数。