GEOLOGICALBULLETINCHINAVol.34,No.4Apr.,2015收稿日期：2014-09-09；修订日期：2015-03-19资助项目：中国地质调查局项目（编号：41、12）作者简介：独文惠(1987-)，女，硕士，助理工程师，从事遥感在地学方面的应用研究。E-mail:通讯作者：**文(1972-)，男，博士，副教授，从事矿产普查与勘探、矿产资源评价的科研和教学工作。E-mail:基于ETM数据和GIS技术的缅甸铜金矿成矿预测独文惠，王功文，陈永清，郭娜娜，郝应龙，冯DUWenhui,WANGGongwen,CHENYongqing,GUONa'na,HAOYinglong,FENGYuan中国地质大学（北京）地质过程与矿产资源国家重点实验室，北京100083StateKeyLaboratoryGeologicalProcessesMineralResources/ChinaUniversityGeosciences(Beijing),Beijing100083,China摘要：在利用遥感数据和地理信息系统（GIS）技术的矿产资源勘查与评价中，矿化蚀变信息提取和构造信息提取起着至关重要的作用。

以缅甸铜金矿为例，利用ETM遥感图像,采用比值法、阈值分割和主成分分析（PCA）提取羟基蚀变信息；利用数字高程模型和遥感数据提取与铜金矿化有关的线环构造，结合地质资料和上述信息，利用证据权和分形方法进行整合，确定该研究区的成矿远景区。结果发现，已知铜金矿均分布于成矿靶区中，综合信息可为下一步的矿产资源勘探工作提供决策依据。关键词：主成分分析(PCA)；蚀变信息；线环构造；GIS；铜金矿化；缅甸中图分类号：P618.41；P618.51文献标志码：A文章编号：1671-2552（2015）04-0786-06DuInterpretationanalysismineralpotentialtargetsbasedGIStechnology:casestudygoldmineralizationMyanmar.GeologicalBulletinChina,2015,34(4):786-791Abstract:structureinformationusingremotesensingdatageographicalInfor⁃mationSystem(GIS)technologyplaysimportantrolemineralresourcesprospectingassessment.Choosinggoldmineauthorsadoptbandratio,thresholdsegmentationprincipalcomponentanalysis(PCA)hydroxylalterationinformationusingETMremotesensingimages.Digitalelevationmodel(DEM)(30mspatialreso⁃lution)extractlinearcircularfaultsgoldmineralization.Combin⁃inggeologicaldataaboveinformation,authorsuseweightfractalmethodfavorableore-formingzonesresearchresultsshowknowncoppergolddepositslocatedwithinhighgradepotentialtargets,integratedinformationcanmineralresourcesdecision-makingfurtherexploration.Keywords:principalcomponentanalysis(PCA);alterationinformation;linearcircularfaults;GIS;coppergoldmineraliza⁃tion;Myanmar自20世纪80年代以来，GIS技术逐渐发展起来，为矿产资源勘查多源致矿信息的综合利用提供了平台。

采用GlS技术进行矿产资源预测评价华纳公司，就是通过各种空间分析方法对各种地质空间数据进行综合，了解地质实体之间的相互关系，掌握各种地质现象与成矿规律的耦合关系，最终对矿产资源进行有效的评价。因此，将多源地学信息和GIS结合，是进行矿产预测的必然趋势和主要手段[1-2]为了贯彻实施“走出去”战略，中国地质调查局成都地质调查中心自2002年以来华纳娱乐,华纳国际，持续开展东南亚中南半岛地区（缅甸、泰国、老挝、柬埔寨、越南）的地质矿产研究。总体上缅甸金矿储量，缅甸矿产资源相当丰富华纳公司华纳公司，但是地质调查与研究的工作研究程度较低，更有不少区域(特别是山区)还属于地质工作的空白第34数据和GIS技术的缅甸铜金矿成矿预测787区。遥感技术越来越多地应用到矿产资源勘查和资源评价，尤其是在提取地质构造信息和遥感蚀变异常信息方面具有优势，在境外地质找矿中发挥着更大的作用。在境外找矿中，采用遥感手段，可以快速有效地开展大区域尺度矿产资源勘查与成矿预测。本次研究以缅甸铜金矿为例，利用ETM遥感影像，采用比值法、阈值分割和主成分分析（PCA）提取蚀变信息，进行图像增强处理后提取与铜金矿化有关的线环构造。基于GIS平台，结合地质概况，利用证据权和分形方法进行整合，确定成矿远景区。

研究结果发现，已知铜金矿床与成矿预测区有较高的相符程度。地质概况缅甸是中南半岛成矿域巨型成矿单元的一部分，成矿地质环境好，矿产资源丰富，现已发现铜、金、铅锌等多个矿床，其中铜金矿床主要形成于中生代、中晚新生代造山带和岩浆弧环境。缅甸境内地层发育齐全，元古宇—新生界均有分布。在东部的中生代盆地中划分出3个大的沉积阶段，依次为石炭纪—二叠纪的白云岩和灰岩沉积，三叠纪—侏罗纪的灰岩、泥灰岩、粉砂岩和页岩沉积及晚侏罗世—白垩纪的页岩、粉砂岩和砂岩（红层）沉积。在西部和中部的古近纪—新近纪盆地中，分布着4个主要的沉积序列，依次为新统—上新统。掸邦边界断层将西部的印—缅，古近纪—新近纪地槽与东部的印支—掸邦台地隔缅甸位于喜马拉雅山南缘东部，古近纪—新近纪时，印度板块向印支地块的俯冲作用使缅甸地区形成东、西分带和边缘高、中间低的区域构造格局遥感数据源及预处理本次地质解译工作的数据源主要采用USGS网站的ETM数据，其空间分辨率在可见光、近红外、中红外波段为30m，热红外波段为60m，全色波段为15m。研究区共有44景影像覆盖，ETM影像图像质量清晰，云量低于10%，植被覆盖低于45%，符合本次研究的需要。

ETM数据和DEM数据主要用于提取与铜金矿化有关的蚀变信息和构造信息。 蚀变信息和构造信息是本次研究的重要组成部分， 为了有效地提取影像的波谱信息和找矿有利信息， 需要对影像进行预处理。ETM 经过辐射校正、几何校正后，进行镶嵌。 遥感构造信息提取线环构造往往与导矿与成矿相关联，线环构造 的解译为地质解译的重要组成部分。DEM镶嵌、裁 剪后，提取山体阴影图。对遥感影像进行滤波、拉 伸等增强处理，结合山体阴影图提取线环构造，为 后期成矿预测提供基础数据（图2）。 遥感蚀变信息提取4.1 干扰信息的去除 ETM 数据主要用于提取与铜金矿化有关的蚀变信息，为了有效地提取影像的波谱信息和找矿有 利信息，对ETM 影像进行辐射校正、几何校正等预处理后缅甸金矿储量，运用比值法消除水体（TM7/TM1）、盐碱地 （TM4/TM3）、植被（TM5/TM4）等影响，用高端或 低端切割的方法对去除云、阴影区、雪等信息的干 扰。干扰信息的掩膜处理，对伪蚀变信息的去除、 蚀变信息的准确提取有重要意义。 4.2 图像增强处理 在蚀变信息提取中，主成分分析（PCA）是应用 较为广泛、较为经典的方法。基于对ETM 波谱信息的分析，分别对TM1、TM4、TM5、TM7 波段和 TM1、TM3、TM4、TM5 波段进行主成分分析，提取 羟基和铁染蚀变信息。

利用TM1、TM4、TM5和TM7四个波段进行主 成分分析（PCA），对代表羟基化物主分量的判断准 则是：构成该主分量的本征向量，其TM5 系数应与 TM7及TM4的系数符号相反，TM1一般与TM5系 数符号相同。用 TM1、TM3、TM4 和TM5 四个波 段进行主成分分析，对代表铁染物主分量的判断准 则是：构成该主分量的本征向量，其TM3 系数应与 TM1 及TM4 的系数符号相反，TM3 一般与TM5 数符号相同。主要采用均值+N倍标准差的方法进行分级， 对于羟基异常，N一般取2~3；铁染异常，N一般取 1.5~2.5 。对铁染、羟基主分量进行阈值分割，得到分级蚀变信息图像，一般分为3级。 788 GEOLOGICALBULLETIN CHINA2015 4.3蚀变结果分析 在将分景的蚀变成果图镶嵌、裁减后缅甸金矿储量，得到缅 甸蚀变成果图（图3）。分析缅甸羟基与铁染的蚀变 成果图分布情况可知，中部地区蚀变信息成块聚集 分布，与典型铜金矿的位置基本相符，南部地区有 的成带状分布，其他各个地区都有零星分布，这可 能与局部蚀变有关。 成矿远景区的划分与分析5.1 证据权法和分形相结合的缅甸找矿远景区圈定 证据权法是由加拿大数学地质学家Agter⁃ bergn [8-9] 提出的，是一种利用矿产形成后的后验概率 来圈定研究区有利成矿部位的数学模型。

最初是 基于二值图像的。它采用一种统计分析模式，以 贝叶斯条件概率为基础，通过对一些与矿产形成 相关的地学信息的叠加复合分析来进行矿产远景 区的预测 [10] 假设研究区被划分成T个像元单位，其中有D个矿点，则随机选取一个像元单位是矿点的概率 将每一种地学信息作为成矿预测的一个证据因子，其对成矿预测的贡献由权重值来确定，定义 分别为证据因子存在区和不存在区的权重值，对于原始数据缺失区域，其权重值为0；B 为这个因子存在的网格单元总数，为不存在的网格 单元数。用C表示证据层与矿床（点）证据层的相 关程度，C定义为： 证据权重法要求各证据因子之间相对于矿点分布满足条件独立缅甸金矿储量，对于n个证据因子，若它们关于 矿点条件独立，则研究区任一网格单元k