近矿蚀变形成的蚀变岩石与围岩在矿物种类、结构、颜色等方面都存在差异, 这些差异在遥感图像中多形成特定蚀变岩石的光谱异常, 这使得遥感蚀变异常信息提取成为可能。基于此, 王晓鹏等^[1]和Ramadan等^{[2, 3]}采用选择性主成分分析(TM1, 3, 4, 5波段)等方法提取铁染(Fe²⁺和Fe³⁺))蚀变异常信息, 采用TM5/7, 4/5, 3/1比值增强处理和TM1, 4, 5, 7主成分分析等方法进行羟基(OH^-)蚀变异常信息提取; 张玉君等^{[4, 5]}基于矿化蚀变波谱特征数学统计分析构建了羟基、铁染矿化蚀变“ 去干扰异常主分量门限化技术流程” , 并在我国矿产资源潜力评价中得到了大量推广应用。ASTER卫星成功发射后, 因其数据包含丰富的岩矿波谱信息、较低的价格而在地质领域广泛应用。Zhang等^[6], Galvã o等^[7], Gabr等^[8], Bedini 等^[9], Freek 等^[10]和Pour 等^{[11 ]}通过分析确定Fe²⁺/Fe³⁺, OH^-基团在ASTER数据中的诊断性光谱, 采用比值、主成分分析、光谱角填图等信息增强处理方法进行了泥化、碳酸盐化、硅化等矿化蚀变遥感异常提取。

由于存在植被、云雪、阴影等影响, 遥感矿化蚀变异常往往与复杂地质背景叠加在一起, 且多以微弱异常信息存在。近年来, 对遥感矿化蚀变异常进行的信息增强与提取已经取得了一些进展。Bertoldi等^[12]和赵志芳等^[13]采用ASTER数据探索了在地形高差大、阴影和云雪覆盖厚的背景下, 基于地形改正或干扰等级划分等提取遥感矿化蚀变弱异常及侵入岩体等信息的方法。分形和多重分形理论因其能较好地表征复杂地质背景下的地质异常, 近年来也逐渐被应用于遥感蚀变弱异常灰度信息分布特征研究及遥感矿化蚀变异常信息提取中(王倩等^[14], 梁珏琦等^[15], 刁海等^[16], 郑桂香等^[17], 赵迪裴等^[18]和 Shahriari等^[19])。

综上所述, 以往研究为遥感矿化蚀变异常信息的增强与提取奠定了较好基础, 但仍存在一些不足: ①对遥感矿化蚀变地质异常信息与背景干扰信息混合叠加这一实质现象重视不足, 导致蚀变信息提取中去除干扰信息时, 许多低弱蚀变异常信息也一并被屏蔽掉, 难以获取隐伏和微弱的遥感蚀变异常信息; ②分形和多重分形理论在遥感矿化蚀变异常提取中虽有探索性应用, 但只是重点关注了铁染等蚀变矿物集合体的蚀变光谱分析及蚀变信息增强与提取, 针对斑岩型铜矿具有的钾化硅化— 硅化绢云母化— 青磐岩化分带诊断性波谱特征分析及弱异常信息增强处理则尚显薄弱。基于此, 本文选取具有典型蚀变分带特征的普朗斑岩型铜矿进行剖析, 采用ASTER数据, 系统分析斑岩型铜矿遥感矿化蚀变分带— 特征矿物光谱特征; 基于多重分形模型, 探索斑岩型铜矿遥感矿化蚀变分带— 特征矿物弱异常增强与提取方法。

普朗斑岩型铜矿区位于滇西北云南省香格里拉县与四川省交界处, 青藏高原东南缘、横断山脉南段北端, 处于德格— 中甸陆块南部, 属义敦岛弧南部。主要出露地层为三叠系图姆沟组(T₃t); 火山岩主要为以安山岩为主的钙碱系列玄武岩— 安山岩组合, 发育有与印支期安山岩同源的中酸性浅成斑(玢)岩, 主要为三叠纪石英二长斑岩和三叠纪石英闪长玢岩; 矿化多集中于格咱岛弧带内, 沿图姆沟组分布的斑(玢)岩形成斑岩型铜矿床(图1)。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 研究区地质图及ASTER遥感图像Fig.1 Geological sketch map and ASTER image of study area

区内中酸性斑岩、富碱性斑岩分布广泛, 矿化蚀变分带特征明显, 表现为从外至内呈青磐岩化— 硅化绢云母化— 钾化硅化的特点; 且青磐岩化是近矿指示标志, 钾化硅化则是主矿体的重要指示标志。

遥感蚀变异常信息属深层次矿致弱异常, 异常分量特征具有非线性幂律关系^{[19, 20, 21]}, 即

A(> S₀)∝ S^-β , (1)

式中: S为能谱密度; A为大于能谱密度某一临界值(S₀)的面积; β 为滤波器的阈值, 不同的β 值对应不同的分形滤波器, 可通过双对数图体现, 所有的直线均服从幂律关系。因此, 遥感蚀变异常信息具有“ 广义自相似性” 与“ 局部奇异性” 特点, 即遥感异常分量的波谱信息在特征空间域具有局部不均一性和各向异性。成秋明等^{[22, 23, 24, 25, 26, 27, 28]}的研究表明, 将非线性理论和方法应用于深层次矿致弱异常信息提取, 可有效地增强矿致弱异常信息。据此, 可利用多重分形“ S-A” 法生成异常分量的双对数(log-log)图。log-log图中不同的直线段代表不同的分形关系, 2条直线的交点所对应的横坐标值(能谱密度值)被视为确定分形滤波器的阈值。借助这些阈值, 可以构造各种异常滤波器和背景滤波器, 通过傅立叶逆变换将其变换到空间域中, 以实现对复杂空间模式(如区域异常和局部异常、背景与弱异常)的分离。

为降低地表覆盖物及混合像元效应的影响, 本文以异常分量灰度值(DN值)及数量分别作为S-A法中的S(能谱密度)和A(累积面积), 求其对数, 制作log-log图; 以双对数图中不同直线段交点作为异常与背景分离的阈值, 获取复杂地质背景条件下不同蚀变分带的特征矿物异常信息。

本文选取ASTER数据的轨道号/行号为AST00167PRDAT0110, 获取时间为2008年8月16日。由于研究区地处山地高原, 受云、雪覆盖的影响, 遥感图像的获取受到制约, 本次研究使用的ASTER遥感数据是从存档数据中经反复挑选后得到的, 但仍有大量云雪覆盖, 影像质量较差。本文通过波段打包、辐射校正和几何纠正后, 经裁剪得到覆盖研究区的ASTER数据。在研究普朗铜矿蚀变分带— 特征蚀变矿物的波谱特征基础上, 结合ASTER数据波谱特征, 设计了斑岩型铜矿遥感矿化蚀变分带精细化信息提取方案^{[29, 30, 31]}。

4.1 石英及硅化蚀变带

分析石英的反射波谱曲线(图2(a)(b))可以看出, 石英在ASTER数据的B1和B2呈吸收谷, B4具有高反射, B12存在强反射, B11为强吸收, 故认为选取B2, B4, B11和B12进行选择性主成分分析, 可以提取绝大多数硅化蚀变信息。硅化异常分量反映在主成分分量PC2中, 其表征特征为ASTER数据B2和B11的贡献与B4和B12相反, 且B12具有高载荷。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 石英波谱曲线及石英异常分量log-log图Fig.2 Spectral curve and log-log figure of quartz abnormal component

利用S-A法生成硅化异常分量PC2的log-log图(图2(c))。选取阈值-2.337 242(即2条曲线的拟合直线的交点, 对应S=0.004 6), 设计滤波器, 可实现复杂地质背景与矿致硅化蚀变异常的分离。

4.2 钾长石及钾化蚀变带

分析钾长石的反射波谱曲线(图3(a)(b))可以看出, 钾长石在ASTER数据的B2具有高反射, B11存在反射峰, B6和B12为小吸收谷, 故选取B2, B6, B11和B12进行选择性主成分分析, 可以提取绝大多数钾化蚀变信息。钾化异常反映在主成分分量PC3中, 其表征特征为在ASTER数据的B6和B12的贡献与B2和B11相反, 且B11具有高载荷。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 钾长石波谱曲线及钾长石异常分量log-log图Fig.3 Spectral curve and log-log figure of orthoclase abnormal component

应用S-A法生成钾化异常分量PC3的log-log图(图3(c))。选取阈值-1.363 799(即2条曲线的拟合直线的交点, 对应S=0.043 3), 设计滤波器, 可进行地质背景与矿致钾化钾长石蚀变异常的分离。

4.3 绢云母与硅化绢云母化蚀变带

分析绢云母的反射波谱曲线(图4(a)(b))可以看出, 绢云母在ASTER数据的B6和B10具有强吸收, B7为强反射, B11为高反射, 故认为选取B6, B7, B10和B11进行选择性主成分分析, 可以提取绝大多数绢云母信息。绢云母异常反映在主成分分量PC2中, 其表征特征为在ASTER数据的B6和B10的贡献与B7和B11相反, 且B7具有高载荷。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 绢云母波谱曲线及绢云母异常分量log-log图Fig.4 Spectral curve and log-log figure of sericitization abnormal component

应用S-A法生成绢云母异常分量PC4的log-log图(图4(c))。选取阈值-1.086 065(即2条曲线的拟合直线的交点, 对应S=0.082 0), 设计滤波器, 进行地质背景与矿致硅化绢云母化蚀变异常的分离。

4.4 绿泥石及青磐岩化蚀变带

分析绿泥石的反射波谱曲线(图5(a)(b))可以看出, 绿泥石在ASTER数据的B5存在反射峰, 在B2和B8为吸收谷, 在B11为强吸收, 故认为选取B2, B5, B8和B11进行选择性主成分分析, 可以提取绝大多数绿泥石信息。绿泥石异常包括于主成分分量PC2中, 其表征特征为在ASTER数据的B5的贡献与B11相反, 且B5具有高载荷。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 绿泥石波谱曲线及绿泥石异常分量log-log图Fig.5 Spectral curve and log-log figure of chlorite abnormal component

应用S-A法生成绿泥石异常分量PC2的log-log图(图5(c))。选取阈值-2.691 826(即2条曲线的拟合直线的交点, 对应S= 0.002 0), 设计滤波器, 进行地质背景与矿致青磐岩化绿泥石蚀变异常的分离。

4.5 绿帘石及青磐岩化蚀变带

分析绿帘石的反射波谱曲线(图6(a)(b))可以看出, 绿帘石在ASTER数据的B5存在反射峰, 在B9为强反射, 在B8为吸收谷, 在B11为强吸收, 故认为选取B5, B8, B9和B11进行选择性主成分分析, 可以提取绝大多数绿帘石信息。绿帘石异常反映于主成分分量PC2, 其表征特征为在ASTERB5和B9的贡献与B11相反, 且B5具有高载荷。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 绿帘石波谱图及绿帘石异常分量log-log图Fig.6 Spectral curve and log-log figure of epidote abnormal component

应用S-A法生成绿帘石异常分量PC2的log-log图(图6(c))。选取阈值-1.427 003(即2条曲线的拟合直线的交点, 对应S= 0.037 4), 设计滤波器, 进行地质背景与矿致青磐岩化绿帘石蚀变异常的分离。

依据上述蚀变分带信息增强提取方案, 获得普朗矿区矿化蚀变分带遥感信息分布特征(图7)。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 遥感矿化蚀变分带信息Fig.7 Mineralized alteration zoning information extracted by remote sensing

钾化带大部分集中于中甸普朗铜矿东北侧及东侧, 部分位于地苏嘎矿矿区周边。钾化蚀变除在普朗铜矿主矿体部位发育较好外, 还在主矿体东侧有大面积出露; 硅化绢云母化在普朗铜矿南东侧独哥附近浓集; 青磐岩化带的绿泥石化与绿帘石化均主要发育于普朗铜矿西北侧; 硅化零星分布于矿区周围。

对提取的矿化蚀变分带异常信息进行了野外查证(查证了20个野外验证点, 详见图8)与室内岩矿鉴定分析(观测点岩矿鉴定照片见图9)。

图8

Fig.8

	Figure Option ViewDownloadNew Window
	图8 遥感矿化蚀变分带野外采样点分布图Fig.8 Map of field sampling point of remote sensing mineralized alteration zoning information

图9

Fig.9

	Figure Option ViewDownloadNew Window
	图9 2号(a)及1号(b)观测点岩矿鉴定照片Fig.9 Microscopic indentification photos for No.2(a) and No.1(b) field sampling points

野外查证和室内分析结果表明: 普朗铜矿主矿体东侧大面积的钾化硅化遥感异常带中见有钾长石和石英富集, 并见有细脉状、浸染状黄铁矿产出; 主矿体部位钾化硅化遥感异常显示较好, 尽管在野外观察为第四系覆盖, 但在其钻孔岩心中见有钾长石斑晶。

1)本文探索性地采用ASTER遥感数据, 在分析斑岩型铜矿钾化硅化— 硅化绢云母化— 青磐岩化等分带特征的基础上, 将多重分形非线性矿致异常增强处理方法引入遥感矿化蚀变分带信息提取中, 获取的弱异常信息经野外验证具较强的吻合性, 表明该方法可靠性强。

2)本文探索的遥感矿化蚀变分带信息增强提取部分成果经野外查证地表为第四系覆盖, 但在钻孔岩心中发现矿化蚀变异常, 印证了该方法可获取隐伏的地质异常信息, 能够实现弱异常与复杂地质背景的分离。

3)本文获取成果与以往研究成果在背景与异常多重分形特征分析方面存在差异, 估计是由于工作区地处山地高原地区, 使用的ASTER遥感图像云雪覆盖量较大、干扰因素较多所致。但是这个问题仍然值得今后深入探讨。