1. 引言
遥感找矿以现代成矿理论为指导，以遥感物理为基础，通过遥感图像处理、解译和遥感信息地面成矿模式研究，同时配合野外地质调查及验证和室内样品分析，以保证遥感找矿的有效性。遥感技术在地质找矿中主要体现在利用影像的空间信息和光谱信息提取与成矿有关的地质体信息，如控矿构造、控矿地层等和蚀变异常信息 [1] 。相较于传统的物化探找矿手段，遥感技术具有宏观、信息量丰富、立体感强、地形地貌特征明显、便于定位等优势，但是随着地质勘查工作的开展，出露地表的矿产明显较少，勘查目标从地表矿、浅部矿、易识别矿转向寻找隐伏矿、难识别矿，找矿勘探的难度日益增大。因此将遥感技术与传统地质方法、现代信息技术相结合是遥感地质应用的发展方向 [2] 。
早期的地质遥感探测主要是依据影像中形状、大小、色调和纹理等特征标志识别、勾绘各种地质构造行迹的形态、产状、分布规律、组合关系及其成因联系等，以分析区域内控矿规律。随着遥感技术的快速发展，解译矿床围岩蚀变信息成为更有利的探测方式。不同蚀变矿物在可见近红外–短波红外光谱区间具有不同的波谱特征。遥感数据多个波段信息的获取，可以鉴别矿物具有的诊断性波谱特征，从而提取矿化蚀变信息。国内外学者基于蚀变岩中所含离子的诊断性波谱特征，利用各种遥感数据开展了多种蚀变遥感异常提取方法研究。比较成熟的矿化蚀变信息提取方法有波段代数运算法、主成分分析、光谱角法及基于图像光谱的矿物端元信息提取等 [3] [4] [5] 。
本文以位于招平断裂带构造中段的大尹格庄、曹家洼、姜家窑、夏甸等金矿富集区域为研究对象，基于Landsat-8影像通过数据预处理、主成分变换解译了研究区矿化蚀变信息，并分析其与招平断裂带构造的空间耦合性，最终圈定研究区找矿远景区，为后续矿产开采提供技术支持。
2. 研究区概况
招远金矿区位于山东半岛西北部。区内黄金资源蕴藏丰富，储量全国第一，分布有特大型金矿、大型金矿、中型金矿和小型金矿多座。招平断裂带经历了早期的韧性剪切和晚期的脆性剪切，并且至少经历了左旋–右旋–左旋的3次活动，力学性质则呈现出扭性–扭张性–扭性–压扭性的变化，是招远金矿区的主要控矿构造，控制了矿田、矿床、矿点的空间分布 [6] [7] [8] 。断裂活动不仅为金矿的形成提供储矿空间，而且为金矿热液的形成提供活化热能和大量流体。在断裂走向和倾向的变化处，热液富集，是金矿体出现的有利部位。
区内金矿化主要有石英脉性金矿化和蚀变岩性金矿化两种，在时空域上的分布规律性十分明显。其空间分布规律主要体现在矿田内的重要工业矿床或主要工业矿脉(体)严格受招平断裂主构造控制，并具有对应产出、等距分布、分段富集和矿体侧伏等特征 [9] 。在时间域上，区内的石英脉性金矿化和蚀变岩性金矿化大致分为黄铁矿–石英阶段，石英–黄铁矿阶段，块状黄铁矿阶段，多金属硫化物阶段和石英–碳酸盐–黄铁矿五个阶段。所以在招远金矿集中区内，招平断裂带附近的石英脉型和蚀变岩性可以作为找矿标志。
本文研究区主要包括招平断裂带主断裂南柳家沟–夏甸段中段及其周围区域。研究区范围内断裂走向NE，倾向SE，倾角40˚左右。研究区构造地质情况如图1所示。
1. 破碎带；2. 第四系；3. 闪长岩；4. 中粗粒二长花岗岩；5. 伟晶不等粒花岗岩；6. 弱片麻状中细粒二长花岗岩；7. 片麻状细粒含黑云花岗闪长岩；8. 片麻状中细粒含角闪黑云英云闪长岩；9. 弱片麻状细粒含石榴二长花岗岩；10. 断层；11. 黑云角闪变粒岩：12. 大理岩；13. 斜长角闪岩；14. 矿床(点)；15. 乡镇；16. 村庄
Figure 1. Geologic significance in the study area [10]
图1. 研究区构造地质图 [10]
3. 遥感数据预处理
Landsat-8卫星于2013年2月11日成功发射。卫星上携带两个传感器，分别是OLI陆地成像仪(Operational Land Imager)和TIRS热红外传感器(Thermal Infrared Sensor)在空间分辨率和光谱特性等方面，Landsat-8数据与Landsat 1-7基本保持一致，具体波段设置见表1。与Landsat-7上的ETM传感器相比，OLI陆地成像仪做了部分调整：Band 5的波段范围调整为0.845 µm~0.885 µm，排除了0.825 µm处水汽吸收的影响；Band 8全色波段范围较窄，从而可以更好区分植被和非植被区域。新增的两个波段为Band 1蓝色波段(0.433 µm~0.453 µm)和Band 9短波红外波段(1.360 µm~1.390 µm)，前者主要应用于海岸带观测，后者应用于云检测。Landsat-8上携带的TIRS热红外传感器主要用于收集地球两个热区地带的热量流失，可以了解观测地带的水分消耗 [11] 。
Table 1. Band parameters of Landsat-8 data
表1. Landsat-8数据波段参数
本文采用的Landsat-8影像2014年5月26日获取。由于原始遥感影像受到传感器自身和传输过程中各种因素的影响，应用前需要经过辐射定标、几何校正、大气校正、干扰信息去除等预处理 [12] [13] 。研究采用ENVI平台中相关模块对研究区Landsat-8影像进行预处理操作。
4. 遥感矿化蚀变信息提取
4.1. 矿化蚀变信息波谱特征分析
由于研究区所在的大型断裂带周围蚀变分带明显，并呈不对称分布，具有很好的构造蚀变–矿化分带，从主裂面向下盘，构造强度逐渐减弱，引起蚀变和矿化同步减弱，形成了绢英岩–绢英岩化花岗岩–钾化花岗岩蚀变分带和黄铁绢英岩化碎裂岩型–黄铁绢英岩化花岗质碎裂岩型(网脉型)–黄铁绢英岩化花岗岩型矿化分带。结合前期相关研究，本研究将研究区泥化蚀变作为主要遥感蚀变找矿标志，铁化蚀变作为辅助性信息 [14] 。
泥化蚀变矿物的波谱特征受羟基基团的影响，在2.2 µm处表现为较强的羟基基团的吸收特征，在2.0~2.5 µm处呈现比硅酸盐、碳酸盐低的反射率。图2(a)中为USGS波谱库中几种典型泥化蚀变矿物的波谱曲线。铁化蚀变矿物主要包含二价铁离子及三价铁离子，从几种典型含铁离子矿物的波谱曲线中可以看出，含三价铁离子类矿物在0.5 µm附近呈现一个拐点，形成清晰的吸收肩，表现出红色不透明物质的特征，其诊断性吸收谱带位于0.8 µm~0.9 µm处。含二价铁离子类矿物在1 µm处具有吸收特征，且在1 µm~1.5 µm处有宽吸收带(图2(b))。
Figure 2. Spectral curves of typical minerals (from USGS). (a) Spectral curves of hydroxyl minerals; (b) Spectral curves of iron minerals
图2. 典型矿物波谱曲线(取自USGS光谱库)。(a) 含羟基类矿物波谱曲线；(b) 含铁离子类矿物波谱曲线
结合研究区泥化蚀变和铁化蚀变指示性矿物，绢云母化、绿泥石化、黄铁矿化中含有羟基基团和三价铁离子。对照Landsat-8相应波段，含有羟基基团的矿物在B7波段呈吸收特征，在B6波段有较高的反射率；含三价铁离子的矿物在B2、B5波段呈吸收特征，在B4、B6和B7波段有较高的反射率。
4.2. 矿化蚀变信息提取
本文以Landsat8影像为数据源，采用主成分分析法对研究区内的泥化及铁化蚀变信息进行提取。主成分分析法在数学上又称KL变换，是多元统计中最常用的一种降维方法，可以减少信息间的冗余，在遥感矿化蚀变信息提取中有较为成熟的应用 [3] [15] 。
依据图2中不同矿物波谱曲线特征，分别选取B2，B5，B6，B7和B2，B4，B5，B6进行主成份运算以提取泥化和铁化蚀变信息。含羟基类矿物所在的分量特征应该是B2，B6波段贡献与B5，B7波段相反，含铁离子类类矿物所在的分量特征为B2，B5波段贡献与B4，B6波段相反。
分析蚀变矿物主成分分析特征向量表(表2)，泥化蚀变矿物主成分变换表中第4主分量(PC4)符合含羟基类矿物波谱特征，然而在B5波段处应表现为负值，因此将PC4取反后作为泥化蚀变矿物信息；铁化蚀变类矿物主成分变换表中第2主分量(PC2)符合含铁离子类矿物波谱特征，且符号符合规律，将PC2作为铁化蚀变矿物信息。为进一步突出蚀变信息分布、去除孤立噪点，对蚀变类矿物所在主分量影像进行5 × 5高斯低通滤波。最后对影像以其平均值加标准差为阈值，将大于阈值的区域作为矿化蚀变区进行拉伸，得到羟基和铁化蚀变矿物信息分布图，底图采用Landsat-8全色波段(图3)。
分析蚀变矿物强度分布图，研究区羟基蚀变信息强烈区主要沿招平主断裂带条状分布，其中大尹格庄金矿、夏甸主井矿和陡崖曹家金矿周边羟基蚀变信息最为丰富，其他地区也有团状聚集。铁化蚀变信息强烈区主要分布在范家庄金矿和小后夼金矿之间，在研究东北部也有片状分布。总体上而言，泥化蚀变与铁化蚀变的空间套合性不明显，前者与已有金矿分布吻合度较高，后者与已有金矿分布吻合度较低，与文献 [14] 研究结果相一致。
Figure 3. The alteration information in the study area. (a) Argillic alteration; (b) Iron alteration
图3. 研究区泥化、铁化蚀变信息解译。(a) 泥化蚀变信息；(b) 铁化蚀变信息
Table 2. Eigenvector of principal component analysis
表2. 蚀变矿物主成分分析特征向量表
5. 找矿远景区预测
以泥化蚀变空间分布为主要遥感找矿标志，辅以铁化蚀变空间分布特征，圈定找矿远景区9处(图4)。其中I级找矿远景区4处，II、III级找矿远景区各2处。I-1、I-2、I-3及I-4找矿远景区分别位于大尹格庄金矿北侧、曹家洼金矿西南侧、姜家窑金矿西北侧及夏甸主井矿附东南侧，可能为已有金矿床矿体深部延伸，浅层地表泥化蚀变信息丰富，并有少量铁化蚀变信息，为优先推荐找矿远景区。II级2处找矿远景区分别位于陡崖曹家北部及青草沟金矿周边。相较于I级四处找矿远景区，II级找矿远景区内的泥化蚀变信息稍弱，但仍位于招平断裂带主断裂附近，开采潜力较大。III级2处找矿远景区位于大尹格庄金矿东部、范家庄金矿和小后夼金矿之间，泥化蚀变信息较少，主要为铁化蚀变，有一定的开采价值。
Figure 4. The prospect areas in the study area
图4. 研究区找矿远景区圈定
6. 结论
遥感找矿已成为矿产资源勘查和监测不可替代的高技术手段。本文利用Landsat-8影像多波段、较好空间分辨率和影像信息丰富的特征，利用主成分分析法选择不同的波段组合，提取了大尹格庄至夏甸金矿一线金矿富集区域的泥化、铁化蚀变信息。由蚀变信息分布图可知，研究区内泥化蚀变信息较为强烈，主要沿主断裂带附近带状分布，与研究区已有金矿分布基本一致；铁化蚀变信息与泥化蚀变信息套合不明显，主要分布在后仓、范家庄和小后夼附近。综合研究区蚀变信息分布共圈定9处找矿远景区。其中I级找矿远景区4处，主要为已有金矿深部延伸，为优先推荐开采区；II级找矿远景区也位于已有金矿周边地区，但其泥化蚀变信息相对较弱。III级2处找矿远景区与招平断裂段距离相对较远，但其铁化蚀变信息丰富，也是推荐开发区域。
然而需要注意的是，研究区农田、水体等地物类型复杂，为提取地表矿化蚀变信息造成了较大干扰。本文将植被、水体信息进行了掩膜去除，但不可避免的也会丢失相应区域的矿化蚀变信息。此外，Landsat-8数据波段信息有限，不能识别更细化的矿化蚀变信息。因此，后续研究可以通过采用其他多源遥感数据如ASTER、MODIS、Hyperion等，提取蚀变矿物、植被异常及地热异常等与热液蚀变相关的多元信息，提高遥感矿物探测准确度。总体而言，遥感技术为地质勘查提供了新的途径，但其获取信息的准确性、可靠性受到多种因素的影响，目前仍主要作为辅助手段之一应用于地质领域。本文基于Landsat-8数据获取的相关信息，后续还需要结合物探、化探等多种技术手段，进一步圈定找矿远景区，为研究区进一步矿产开发提供技术支持。