第一作者简介: 刘文兰(1984-),女,在读硕士研究生,主要从事遥感技术在矿产资源勘查与评价方面的研究。E-mail: lty102627@163.com。
遥感地质异常信息提取是遥感地质学的重要研究内容之一。以老挝为例,对研究区ETM+遥感图像采用几何纠正、图像镶嵌、彩色合成、空间变换、图像增强等技术手段,提取地质构造信息; 利用掩模技术去除干扰信息,采用主成分分析和阈值分割技术进行铁染和羟基蚀变异常信息提取; 结合区域地质、地球物理、地球化学等多元地学信息综合分析,开展研究区遥感找矿预测,圈定了5处遥感找矿远景区。研究成果为加快老挝地质找矿工作提供了基础数据和决策依据。
The extraction of remote sensing geological anomaly information is one of the important research methods in GARS. In this paper,taking Laos as an example and using such means as geometric rectification,image mosaicking,color composite,space transformation and image enhancement processing to the ETM+ image in the study area,the authors extracted the geological structural information. The interference information was removed by mask technique,and the iron staining and hydroxyl alteration information was extracted by using principal component analysis and threshold segmentation techniques. Five perspective areas for mineral exploration were delineated on the basis of a synthetic analysis of the multi-geoscience information. The results can provide a reference and effective means for accelerating the geological prospecting work in Laos.
随着矿产资源勘查工作的不断深入, 易于发现的露天矿、浅部矿和已知类型矿床日趋减少, 找矿难度越来越大[1]。遥感技术以其宏观、经济、快速等技术特点, 尤其是在提取地质构造信息和遥感蚀变异常信息方面具有得天独厚的技术优势, 被越来越多地应用到矿产资源勘查和资源评价, 尤其是在境外地质找矿中发挥着越来越重要的作用。
线性构造在遥感图像上往往具有清晰的线性形迹, 常常表现为直线状、弧形、波状色线等特征, 具有一定规模时则构成色调异常带。线性构造提取常用的方法有滤波、主成分分析、多重主成分分析等[2]。20世纪70年代末, 以美国Hunt等[3]为代表的研究人员先后发表了一系列关于矿物、岩石的可见光— 近红外光谱特性等方面的研究成果, 较系统地研究了岩石矿物的光谱特性, 分析了产生这些光谱特性的原因, 并讨论了不同岩类(包括近矿围岩蚀变矿物)光谱特性的变化规律及影响因素。利用这些近矿围岩蚀变矿物特定的光谱特征, 通过遥感数字图像处理, 可有效地提取遥感蚀变异常信息。国内外研究人员已先后提出了多种提取遥感蚀变异常信息的方法, 如比值法、主成分分析法、光谱角法、“ 掩模— 主成分分析— 异常分级(MPH)” 法[4, 5, 6]、“ 去干扰— 主成分分析— 异常筛选” 法[7, 8, 9]、微量信息处理方法、彩色空间变换法及克罗斯塔法等[10, 11, 12], 取得了明显的遥感蚀变信息提取和找矿效果。
本文以老挝为例, 分别从遥感构造信息提取、遥感蚀变异常信息处理以及多元地学信息综合分析等方面, 在研究区开展了成矿地质条件综合分析和找矿预测研究。研究结果证实, 本文方法在老挝及同类地区遥感找矿中是可行的。
老挝地跨昌都— 思茅— 南帮和印支成矿省, 其大地构造位置处于印度板块和太平洋板块之间、欧亚板块与印度板块缝合线的东侧, 在地质历史上, 构造运动与岩浆活动十分强烈。由于老挝地处冈瓦纳和欧亚两大陆相互离散拼合区[13], 受该构造环境影响, 断裂性质极其复杂, 以张剪性为主, 且多呈北东(NE)、北西(NW)向展布, 形成了一系列相互交错的线性、环形构造, 在遥感图像上表现尤为突出。岩浆侵入作用分布较为广泛, 岩体形成时代主要有元古宙、华力西期、印支期、燕山期和新生代, 其中以印支期岩体最为发育(图1)。强烈的构造运动与岩浆活动为矿产资源的形成与分布提供了必要条件和物质基础。尽管老挝的地质工作程度较低, 已开展的地质填图、区域地球化学填图和航空地球物理调查等工作有限, 但现有的工作成果已经反映出老挝具有巨大的矿产资源开发潜力。
在漫长的地质演化历史过程中, 老挝的区域成矿作用无论在矿产种类还是在矿床类型上均呈现出多样性, 主要为热液型、矽卡岩型、沉积型、砂矿型和风化淋滤型, 其次为岩浆型、伟晶岩型和火山岩型, 并且在矿种上具有一定特色。目前, 老挝已经发现许多金属矿床、矿点和矿化点, 钾盐、锡、铁和蓝宝石是老挝的特色矿产, 石油和天然气被大国石油公司看成是最有前景的能源矿产。
根据现有的工作成果可知[14], 老挝的区域成矿时期比较集中(主要集中在华力西晚期、印支期、燕山晚期、晚第三纪和第四纪)。铜、铅锌、钼、锡、金、铁等内生多金属矿化主要与中生代岩浆活动有关, 其中铁、金、铜、铅锌等矿化主要与花岗闪长岩有关, 钼、锑、锡矿化一般与石英二长岩及花岗岩有关。此外, 与碳酸盐岩和砂岩接触的新生代流斑岩也引起金、银、铜、铅锌等弱矿化。上古生界和中生界地层中有同生层状矿床, 如赋存于二叠系地层中的巨大的石灰岩矿床、赋存于侏罗系砂岩中的弱铜矿化等。晚石炭世、二叠纪和三叠纪为成煤期。新生代发育有厚层蒸发岩, 形成石膏(硬石膏)和钾盐类矿床。第四纪冲积物和残坡积物中, 往往有金、锡石和宝石等富集并达到工业品位, 这些矿床一般都产在距原生矿床或矿源岩不远的范围内, 是目前矿产开发的主要对象。此外, 还有一些矿化富集源于表生风化作用, 如新生代玄武岩风化形成的铝土矿、蚀变花岗岩风化形成的高岭土、含锡花岗岩风化形成的含锡红石土、超基性岩风化形成的菱铁矿等[15]。
本次遥感地质解译工作选取的遥感数据源是具有中等空间分辨率和光谱分辨率的美国陆地卫星ETM+图像, 成像时间为1999— 2001年, 共计18景。图像清晰, 对比度好, 纹理结构特征明显, 利于地质构造解译及蚀变异常信息提取。
采用的DEM数据是基于ASTER立体相对生成的, 可作为构造信息提取的有效补充。
利用遥感图像处理软件(ENVI, ERDAS), 对研究区的ETM+图像进行大气校正和几何纠正, 启用拼接功能对18景ETM+图像和DEM数据分别进行镶嵌与拼接, 通过矢量边界对拼接后的数据进行裁剪。在后续研究中均以裁剪后的区域为研究对象, 开展空间变换、图像增强等处理, 尽可能地突出感兴趣信息(如线性构造信息、隐伏构造信息、弱异常信息等), 掌握图像的有效特征以获取高质量的影像特征和相关信息, 以便于地质人员进行分析和解释。
遥感图像波段的选取对提取地质构造信息非常重要, 不同的波段组合在不同地区呈现出不同的效果。根据OIF(optimum index factor)指数[16]最佳波段组合原理, 结合研究区的矿化蚀变特征, 本次研究最终选用ETM7(R), ETM4(G), ETM2(B)波段组合并与全色波段进行数据融合的假彩色合成图像作为遥感地质解译的基础图像(图2); 利用ASTER数据自动生成全区的DEM数据(图3), 并在此基础上制作研究区山体阴影图像, 从立体的角度观测和推断地质构造及断裂的空间展布, 使构造信息的解译更加全面和直观, 可作为常规遥感地质解译的重要补充。
阴影、水体、植被及白泥地等地物是遥感蚀变异常信息提取的干扰信息, 在进行蚀变信息提取之前要对其进行去除, 以免蚀变信息失真; 根据不同地物的不同波谱特征信息而采用不同的去除方法, 并尽可能地减少对遥感蚀变异常提取产生影响[17]。
老挝地处热带, 常年有植被覆盖, 势必对遥感蚀变异常信息的提取造成一定的干扰。本次研究主要通过提取植被覆盖区植物的地球化学异常来进行蚀变异常信息的提取, 即利用TM1— 5, 7波段做掩模主成分分析, 根据第5, 6主成分的贡献决定对异常结果的选用, 并结合ETM3/ETM4比值进行相关异常信息提取; 根据植被的波谱特性, 对上述结果进行对比分析, 最大限度地去除植被对异常提取结果的干扰。而由于阴影与水体有相似的波谱特征, 因此可以将二者一并去除, 采用比值ETM6/ETM1去除水体和阴影。同时根据白泥地的波谱特征, 采用ETM3高端切割去除白泥地[18]。
构造条件是矿床形成的重要因素之一[19], 根据遥感影像特征和地质成矿理论, 提取的构造信息主要有断裂构造和环形构造。其中断裂构造对矿产的控制作用十分明显, 矿床常常富集在与主干构造斜交或平行的断层或节理中, 并且以线性体或环形体等结构信息反映在遥感影像上。这些成矿有利部位既有利于成矿元素的迁移与富集, 又提供了成矿空间, 有利于形成规模较大、品位较高的矿床。各种不同类型金属矿床的分布与线性构造的关系, 大体上有以下4种形式: ①表示构造单元边界的深、大断裂(带)的巨大型、大型主干线性构造(带), 控制着成矿带或矿化集中带的展布; ②交切于上述线性构造(带)的稍次级区域性大断裂(带)的线性构造(带)及其切割的地质块体, 往往是成矿亚带、成矿区的出现部位; ③在前述线性构造叠合地带上的次级线性构造带, 往往控制着矿田、矿区的分布; ④规模较小、但密集程度较高的次级线性构造, 通常是矿床、矿体的产出部位。
本次遥感构造解译以ETM+图像解译为主, 结合利用DEM数据制作的山体阴影图进行综合解译(图4), 主要的方法有两种: ①直接解译, 根据岩石或者地层发生位移, 构造线发生错位, 地层重复或者缺失, 陡崖的直线状分布, 岩体、岩脉、火山口呈线状分布等地质解译标志进行构造解译; ②间接解译, 根据山脊、湖泊、沼泽错位, 对头沟、对头河的出现, 直线状沟谷、洼地, 不同地貌区沿直线相接, 水系特征和地表水体异常等[20]地貌解译标志进行构造解译。
从图4可以看出, 研究区的构造主要分布在中部和北部, 有成群成带分布特征, 整体以NE向、NW向为主, 北部局部地区近东西(EW)向断裂较为发育。其中波里坎晒— 他曲断裂、沙拉湾— 塞公断裂等NW向断裂带横贯全区, 这些断裂控制着区域的沉积建造、岩浆活动和变质作用, 也控制着与其有关的矿产, 是重要的控矿构造。其强烈的构造运动和热液活动可使原来相对分散的铜、铅锌、铁、金元素得到进一步聚集, 并在有利的构造空间部位富集成矿。从线性构造密集度可以看出, 本次研究共圈定了6处线性构造密集区, 这些区域往往是成矿的有利地段, 既有利于成矿元素的迁移、富集, 又提供了成矿空间, 是今后开展矿产资源勘查的重点地区。
根据前人研究成果可知[21], 研究区的区域成矿作用时间比较集中, 主要集中在华力西晚期、印支期、燕山晚期、早第三纪和第四纪, 其中铁、有色金属和稀有金属等矿产主要集中在燕山晚期(晚侏罗世至早白垩世); 以热液型与矽卡岩型矿产为主, 主要的蚀变矿物有黄钾铁矾、绿泥石、绿帘石、高岭土等, 而这些蚀变矿物主要与铁染异常和羟基异常有关。据蚀变矿物的波谱特征可知, 黄钾铁矾、绿泥石等蚀变矿物分别在ETM+1, 3, 4, 5波段和ETM+1, 4, 5, 7波段具有特殊的光谱特征, 利用不同波段对不同地物类型的强反射与强吸收特性, 分别对这4个波段进行主成分分析[22]; 通过对各主分量进行分析和定量化评价, 得出异常主要集中在第4主成分中; 运用阈值分割技术, 得到研究区铁染蚀变和羟基蚀变的异常分布图[23](图5)。在完成遥感异常筛选的同时, 通过对区域成矿条件的综合分析与研究, 对区域遥感异常进行分类, 对可能属于同一类型、空间位置关系密切或成矿地质背景相近的遥感异常, 用闭合的人工包络线进行归类, 将其定义为具有一定找矿意义的遥感异常。本次研究共圈定了遥感找矿异常区27处。
![]() | 图5 老挝铁染蚀变与羟基蚀变异常分布图Fig.5 Maps of distribution of iron staining and hydroxyl alteration anomalies of Laos |
从图5可以看出, 研究区内的遥感异常以铁染异常为主, 其分布与区内的构造格架大体一致, 多分布在断裂构造和环形构造较为发育的地区; 大多数位于冲沟内植被覆盖区的羟基异常主要为植被干扰因素未完全消除所致; 而铁染异常和羟基异常套合度较好地段往往是成矿有利地区, 对在这些地区寻找浅层低温热液型、矽卡岩型以及火山岩型矿床具有一定的指示作用, 可为遥感找矿远景区的圈定提供基础资料和参考依据。同时考虑到老挝地区常年有植被覆盖, 且第四系松散堆积物较为发育, 尽管最大限度地对其进行了筛选和剔除, 但对矿致异常和岩性异常的提取仍存在着一定的影响和干扰, 须结合地、物、化等综合信息进行遥感找矿预测工作。
本次研究通过遥感构造解译和异常信息提取获得了大量的遥感信息, 结合区域地质、矿产等多元地学信息, 建立了以遥感信息为主的多元信息遥感找矿预测技术方法流程(图6)。
按照上述遥感找矿预测技术流程, 以遥感构造信息和蚀变异常信息为主, 圈定出线性构造密集区和遥感异常区; 并结合所收集的地质、矿产等资料, 对老挝进行了以遥感信息为主的多元信息综合分析和找矿预测。考虑到老挝目前的区域成矿研究程度较低, 还不可能划分出比较细致的成矿区或成矿带, 因此本次研究主要以遥感信息为主, 结合己知资料, 共圈定出乌多姆塞铜、金、铅锌多金属找矿区, 南俄河水库西铁、铜、锡、铅锌银多金属找矿区, 川圹铁、锰及贵金属找矿区, 波里坎晒东铜、铅锌找矿区和占巴塞铜、铅锌矿找矿区等5处找矿远景区(图7)。其中乌多姆塞、川圹和占巴塞找矿区在前人资料[21]中有一些反映, 而南俄河水库西和波里坎晒东找矿区是根据本次研究新提出的找矿区, 研究成果为加快老挝地质找矿工作提供了基础数据和决策依据。
1)利用ETM+遥感图像的ETM7(R), ETM4(G), ETM2(B)波段组合进行假彩色合成, 并与全色波段进行数据融合, 获取的图像更有利于目视解译和遥感构造信息的提取; 与此同时, 利用DEM生成山体阴影图像能从立体的角度观测构造信息和断裂的空间展布, 使得构造信息的解译更加全面, 可作为常规遥感地质解译的重要补充。
2)遥感蚀变异常信息提取不受人为因素影响, 可作为矿产资源潜力评价的基础性参量; 结合遥感地质解译结果、区域典型矿床、蚀变带的吻合程度等进行综合分析, 可为开展区域矿产资源勘查和资源评价提供基础数据和参考依据。
3)由于矿床的形成过程是极其复杂的, 因此利用遥感图像进行地质构造解译和遥感蚀变异常信息提取、开展矿产资源勘查与评价还只是一种战略手段; 加之研究区内植被和第四系覆盖严重、交通不便等因素, 导致预测结果可能存在一定的偏差。但是, 随着遥感数据空间分辨率的不断提高, 以及高光谱、高分雷达数据的不断发展, 遥感技术在地质找矿中必将发挥更大的作用。
志谢: 中国地质大学(北京)王功文教授提供了研究所需的老挝1∶100万比例尺的地质图件,在此表示诚挚的谢意。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|