陈卫东
摘要:本文主要以中国南岭地区为例,利用卫星遥感数据为信息源,通过断裂构造信息提取和矿化蚀变信息提取进行成矿预测。
关键词:遥感找矿信息信息提取断裂构造信息
0引言
花岗岩发育地区地形条件复杂,通行不便,植被较发育,地面地质工作困难,而且地质工作程度相对较低,而采用遥感技术(遥感找矿法)可以有效地弥补上述不足。花岗岩地区与锡矿有关的矿化蚀变、与脉状锡矿体关系密切的断裂构造、花岗斑岩脉等,都可以通过遥感找矿信息提取技术有效地提取出来,从而大量增加了地质信息。笔者以中国南岭成矿带花山地区为例,探讨了花岗岩地区利用遥感找矿信息提取技术预测锡矿找矿远景区的方法。
1遥感找矿信息与地物波谱分析
1.1遥感找矿信息提取的地质基础对山地地区典型锡矿田(区)地质特征、成控矿特点分析可知:锡矿床以构造蚀变岩型、接触带矽卡岩型和残留顶盖矽卡岩型最具找矿前景,其次是石英斑岩脉型。区内蚀变以矽卡岩化、云英岩化、电气石化、硫化物化、钾化、硅化为主,次为萤石化、绢云母化、高岭土化等,其中与锡矿化关系密切的是矽卡岩化、云英岩化、电气石化、黄铁矿化。因此,本区各类型锡矿床的遥感找矿信息提取的地质基础为:
1.1.1接触带矽卡岩型和残留顶盖矽卡岩型锡(铅锌)矿床。主要提取矽卡岩化蚀变信息,在地表多数矽卡岩常风化成为铁帽,其波谱特征与花岗岩和灰岩等其它沉积岩有较明显的区别。
1.1.2构造蚀变岩型锡矿床。主要提取云英岩化、电气石化、硫化物化、钾化、硅化等蚀变信息,结合进行断裂构造信息提取。锡矿体往往明显地受断裂控制。
1.1.3石英斑岩型锡矿床。由于石英斑岩脉在地貌上一般呈比较明显的线(带)状正地形,而且,石英斑岩本身也受一定的断裂构造控制。所以,对此类矿床重点进行线状影像信息提取。这是遥感信息的优势所在。
1.2典型地物波谱特征分析为了研究花岗岩地区锡多金属矿化及有关围岩、蚀变岩(矿)石的光谱特性,分别在多个山区等已知锡多金属矿区采集了有代表性的岩石、矿石、矿化蚀变标本,利用IRIS型波谱测试仪进行了反射波谱测试。绘制典型岩、矿石样本波谱曲线图。由波谱曲线圈看出,产于外接触带的矽卡岩、铁帽、灰岩具有大致相同的波谱特征,处于低反射率区,曲线平滑:沉积的中细粒砂岩处于中反射率区,在TM1和TM7有小幅吸收谷:分布于岩体内接触带以及岩体内部的矽卡岩、云英岩化花岗岩处于高反射率区,且在TM1和TM7有较深的吸收谷。
分析表明,本区主要找矿标志铁帽、矽卡岩化、云英岩化等蚀变与地层、岩浆岩等围岩之间存在明显的反射率差异,也具有较明显的波谱特性,具备遥感找矿信息提取的物质基础和理论条件。
2断裂构造信息提取与解译
2.1断裂构造信息提取技术方法断裂构造信息提取的基本方法包括空间卷积变换、主成分分析、假彩色合成、HIS变换、人机交互式解译等技术方法。
2.2断裂构造解译结果通过线状影像增强处理和人机交互式解译,新增断裂1355条。特别是花山岩体内部断裂极其发育,近EW向、近SN向(NNE向)、NE向、NE向都有发育,既有规模较大的NE向、NE向和近SN向断裂,也有近EW向、近SN向和NE向的断裂密集带。断裂的形成时代应该与岩体的形成时代相同或较晚,因此,一般认为,岩体内断裂主要形成于燕山晚期,这与南岭地区锡多金属主要成矿期是一致的。
3遥感异常信息提取
3.1遥感异常信息提取方法遥感异常信息提取的重点是:矽卡岩化、铁帽(褐铁矿化)、云英岩化等。研究认为采用矿化蚀变信息增强处理(比值方案选取)、矿化蚀变信息提取(主成分析)、光谱角分类等技术方法得到的结果较为理想。
3.1.1比值增强处理根据以上地物波谱特征分析,考虑到消除照度和阴影等因素对信息提取的影响,采用比值方法可以取得较好的效果。又由于矿物岩石蚀变一般在TM5有较高的反射率,而一些干扰信息如植被在TM4有较小的反射率,水体或含水分较高的粘土等在TM1有较高的反射率。所以,在设计比值时一般将TM5作分子,其它TM波段作分母。优选4个增强弱矿化蚀变信息的比值(TM5/TM1、TM5/TM4.TM5/TM7、TM7/M1),这些比值均能不同程度地增强某种或几种蚀变信息,压抑某些干扰信息。如TM5/TM1增强富铁类岩石,对铁帽反应敏感,有助于提取矽卡岩化、云英岩化、磁铁矿化等信息:TM5/TM4增强富铁和富二氧化硅岩类(两者在0.85μm处有吸收带),有助于提取褐铁矿化、硅化、黄铁矿化信息;TM5/TM7增强与粘土化有关的矿化蚀变信息。有助于提取云英岩化、磁铁矿化、电气石化、矽卡岩化蚀变信息;TM7/TM1增强氧化铁帽信息,有助于提取铁帽。
3.1.2主成分分析选择上述4个比值经过主成分分析,可以得到4个主分量(PC1、PC2、PC3、PC4)。选择集中了主要信息量的PC1、PC2、PC3合成彩色图像,比原始图像的KL变换结果更好地显示弱信息。
3.1.3光谱角分类法采用光谱角分类法是对岩石进行波谱形态识别的主要方法之一,它将光谱数据视为多维空间的矢量,利用解析方法计算像元光谱与光谱数据库光谱或像元训练光谱之间矢量的夹角,根据夹角的大小来确定光谱间的相似程度,以达到识别地物的目的。由于光谱角对增益不敏感,因此,能够减弱因地形和照度引起的增益变化。
以上述合成彩色图像作为分类依据,训练各类异常的光谱样区,取其在各对应波段的光谱反射率的平均值建立3维空间的矢量样本模型,计算各像元在3个波段的灰度值数据矢量与样本模型之间的夹角,根据夹角的大小确定光谱间的相似程度,夹角越小,相似程度越高:反之亦然。矢量夹角(θ)通过矢量间相似系数(即矢量夹角的余弦cosθ)计算获得。通过分析给定像元与样本模型之间光谱相似程度量值,或设定光谱角度阈值,试验认为本区取值0.05可较好地区分矿化蚀变遥感异常与非矿化蚀变遥感异常。
分类后处理,包括过滤(SIEVE)处理,最小阀值为2:膨胀(CLUMP)处理,窗口尺寸采用3×3。光谱角分类能快速客观全面地反映遥感异常的光谱分布,但不能识别“异物同谱”和“同物异谱”现象,有可能漏画或多画遥感异常;异常图斑细碎。需进行一定的人工修改。
3.2遥感异常信息提取结果与评价分类后的遥感异常图,按照其分布的位置与花岗岩体的关系,分为3种类型:即分布于岩体内的遥感异常,有21处;分布于接触带的遥感异常,有7处;分布于岩体外围的遥感异常,有11处。按照遥感异常所代表的矿化蚀变类型,可以分为矽卡岩化遥感异常、云英岩化等蚀变遥感异常、粘土化遥感异常以及铁帽等4类。
各类型异常约39处。经过野外调查以及相关地质找矿单位的已知矿产资料查证,已知矿化地段15处。将已见矿的遥感异常的处数除以遥感异常的总处数,再乘以100%,初步计算出了遥感异常分布区的见矿率评价指标为:到本项目野外工作结束时为止,遥感异常的总见矿率为38.4%;分布于花岗岩体内的遥感异常见矿率为47.5%:分布于花岗岩体接触带的遥感异常的见矿率为43.0%;分布于岩体外围的遥感异常的见矿率为18.2%。从岩体内部到远离岩体,遥感异常的见矿率由高到低的分布规律可以看出锡矿最有利的找矿地段在岩体内部,其次在岩体的内外接触带。从目前遥感异常的总见矿率仅为38.4%可知,本区锡矿找矿的前景还很大。
4遥感成矿远景预测
4.1预测方法及预测图编制遥感找矿预测的基本方法是以遥感异常为主,辅以多源信息(多专题信息)叠加综合分析法。主要专题图层及叠加顺序为:地理信息图层、地层岩浆岩解译图层、地质构造解译图层、化探异常图层、遥感异常图层、已知矿点图层。最后,在上述叠合图上综合分析,圈定遥感找矿远景区,根据遥感异常的类型和预测依据划分远景区的类型和级别。
4.2预测结果及结论本区共划分了23个遥感找矿远景区,其中Ⅰ级遥感异常远景区10个;Ⅱ级遥感异常远景区4个;Ⅲ级遥感异常远景区9个。
5结束语
综上所述,采用遥感找矿异常信息提取与断裂构造信息提取、多源信息综合分析的遥感找矿方法在花岗岩地区寻找锡矿可以发挥较好的作用。