高光谱遥感地空综合预测方法在新疆卡拉塔格地区铜金矿床找矿中的应用

2017-01-17孙永彬王瑞军魏本赞汪冰董双发李存金李名松

河北遥感 2016年3期

关键词：褐铁矿红山矿床

孙永彬王瑞军魏本赞汪冰董双发李存金李名松

（核工业航测遥感中心，河北石家庄 050002）

技术应用

高光谱遥感地空综合预测方法在新疆卡拉塔格地区铜金矿床找矿中的应用

孙永彬1王瑞军1魏本赞1汪冰1董双发1李存金1李名松1

（核工业航测遥感中心，河北石家庄 050002）

新疆卡拉塔格地区位于东天山成矿带，地质矿产资源丰富，铜金及多金属矿床集中，找矿潜力较大。铜金矿床地表矿化蚀变现象明显，发育的褐铁矿化+黄钾铁矾化+绢云母化形成地物反射光谱特征异常，为高光谱遥感地空综合预测方法用于探测铜金矿化蚀变信息成为可能，并成为寻找铜金矿床的有效技术方法。本文在卡拉塔格地区开展了航空高光谱和地面高光谱在铜金多金属矿床找矿研究，分别采用HyMap航空成像光谱仪和FieldSpec Pro FR光谱仪，获取了研究区HyMap航空高光谱分辨率数据、地面准同步定标数据、地面高光谱数据，提取并筛选了基于HyMap航空高光谱影像的成矿有利蚀变信息，在此基础上对该区红山铜金矿床进行了航空-地面高光谱综合剖析，并结合矿区地质背景，建立了红山铜金矿床高光谱遥感地空综合找矿模型，经对比分析在该区圈定了两处找矿预测区。说明高光谱遥感地空综合预测方法在该区得到了很好的应用效果。

航空高光谱；地面高光谱；诊断性光谱特征；地空综合找矿模型；红山铜金矿床

前言

遥感地质找矿是遥感信息获取、含矿信息提取以及含矿信息成矿分析与应用的过程，近年来，被广泛应用于地质、矿产资源及相关环境的调查中[1-7]。航空高光谱是当前遥感的前沿技术，通过高光谱成像所获取的地球表面的图像包含光谱维信息融合为一体，即“图谱合一”[8]，具有波段数多，波带窄，对地物能定量分析等优点，是二十世纪末期以来遥感领域最大的技术进展[9]。航空高光谱遥感技术相对于传统的星载高光谱遥感技术，可以获取高空间分辨率的高光谱遥感数据，因而对微小地物具有更强的识别能力[10]。通过高光谱矿物填图，可以大面积、快速提取蚀变矿物[11]。王润生等[12]系统总结了高光谱矿物填图技术流程、工作方法和技术体系。中国国土资源航空物探遥感中心在新疆土屋东-三岔口地区，应用HyMap航空成像光谱数据填绘出了白云母、绿泥石、绿帘石、绿泥石和绿帘石组合、高岭石、蒙脱石、透辉石、透闪石、蛇纹石、褐铁矿、方解石等10多种矿物种类[13]。

本文以新疆卡拉塔格地区获取的HyMap航空高光谱数据、地面准同步定标数据、地面高光谱数据为例，提取并筛选了基于HyMap航空高光谱影像的成矿有利蚀变信息，在此基础上对红山铜金矿床进行了航空-地面高光谱综合剖析，并结合矿区地质背景，建立了红山铜金矿床高光谱遥感地空综合找矿模型，圈定了两处找矿有利地段，找矿效果显著。

1、研究区地质概况

研究区大地构造位于东天山吐哈盆地南缘、大南湖-头苏泉晚古生代岛弧带（图1）。

区内地层整体呈北西向产出，由中奥陶统荒草坡群大柳沟组（O2Hd）、中-上志留统红柳峡组（S2-3h）、下泥盆统大南湖组（D1d）、上石炭统脐山组（C2qs）、中二叠统阿尔巴萨依组（P2a）、上二叠统库莱组（P3k）、下侏罗统三工河组（J1s）、全新统（Qh）组成。

该区经历了多期次的构造变动，区域构造形成类型多样、特征复杂，褶皱、断裂和火山机构发育，区域性的断裂构造控制着该区火山岩和侵入岩的展布，次级断裂控制着矿化蚀变带的分布，更次一级构造裂隙带控制着矿化体的产出。主要发育北西、北北西和北东东向三组断裂构造，同火山断裂发育，含矿火山热液沿断裂活动强烈。

区内火山活动强烈，火山岩分布广泛。岩石类型以中基性火山岩为主，中酸性火山岩次之，从奥陶系至二叠系均有出露，其中奥陶系最为发育，其次为泥盆系、石炭系，二叠系局部出现。侵入岩较为发育，主要为一套古生代侵入体，包括中酸性闪长岩、花岗闪长岩、花岗岩等侵入岩。其中志留纪侵入岩由老到新依次为英云闪长岩、花岗闪长岩、二长花岗岩、花岗斑岩。区内经受了区域埋深变质作用，主要变质矿物为绿泥石、钠长石、葡萄石、石英、阳起石、方解石、绿帘石。

区内已知矿床自北西向南东依次产出有：红山铜金矿床→梅岭铜锌金矿床→红石铜金矿床→红海铜锌金矿床。

2、高光谱遥感地空综合预测方法用于寻找铜金矿床的原理

2.1 高光谱遥感地空综合预测方法的依据

不同矿床类型由于岩石岩性和蚀变矿物不同，它们的光谱特征有明显的差异[14]。铜金矿床属于内生金属矿床，这类矿床与地表蚀变密切相关，尤其地表蚀变矿物直接指示矿床的成因类型和指导找矿。利用航空高光谱遥感影像进行蚀变矿物提取，结合地表岩石光谱测试分析，应用HyMap航空成像光谱仪和美国ASD公司的FieldSpec Pro FR光谱仪对不同类型矿床的光谱特征进行测量，分析其不同类型矿床的光谱特征，进行蚀变矿物遥感信息提取，开展高光谱遥感地空综合预测研究，进而指导找矿。

2.2 蚀变矿物诊断性光谱特征分析

岩矿的光谱包含有一系列特征吸收谱带。每一个特征谱带或谱带组合与岩矿内部微粒的物质属性存在一定的对应关系[15]。在短波红外光谱区域，褐铁矿、黄钾铁矾、白云母、绿泥石等具有可以识别的诊断性光谱特征，特征吸收波段的深度与岩石中这些矿物的含量密切相关[16]。矿物诊断性光谱吸收包括金属阳离子在可见光区域的电子过程以及阴离子基团在近红外区域的振动过程。蚀变矿物的可见光-近红外范围内几种官能团主要包括Fe2+、Fe3+、Al-OH、Mg-OH、OH-、CO32-。

一般阴离子诊断谱带位于2000-2500nm光谱区域。而Fe2+，Fe3+和Mn2+诊断谱带一般位于400-1200nm光谱区域[17]。所以根据400-1100nm光谱区间的Fe吸收谱带和2000-2500nm羟基和碳酸根谱带将矿物分为含铁矿物、含羟基矿物和含碳酸根矿物以及其它矿物[18]。褐铁矿的光谱吸收峰是矿物光谱分析、遥感铁染异常提取常用的波段。其中Fe3+离子的特征吸收峰位主要在600～800nm，中心位置在700nm左右，Fe2+离子的特征吸收峰位主要在900～1500nm，中心位置在1200nm左右；白云母的诊断性光谱特征是Al-OH键在2180～2230nm之间的尖锐而深的吸收特征，以及在2340nm和2440nm 附近的较弱Al-OH特征，Clark等[19]在高光谱分辨率下，检测到白云母矿物在2200nm附近的吸收波段随Al含量的增加向长波方向的移动；绿泥石的光谱具有Fe-OH和Mg-OH的诊断性吸收特征，波长位置分别为2235～2255nm和2320～2360nm，这些吸收特征的波长随着绿泥石中铁离子含量的增加而增大[16]。

3、HyMap航空高光谱数据获取及处理

综合考虑仪器参数、成图精度的要求、研究区自然地理、地形起伏、野外工作条件等因素，进行飞行方案的最优设计（表1）。将HyMap成像光谱仪搭载在Y-12飞机上，利用国外引进的HyMap航空高光谱遥感测量系统、地面ASD光谱测量系统。采集高光谱分辨率的HyMap航空高光谱遥感数据、地面准同步的地物光谱定标数据。具体获取过程分为三个步骤：首先对HyMap成像光谱仪进行实验室定标，保证其性能指标达到项目要求，以获取优质的高光谱影像数据；其次进行HyMap仪器的安装和测试，按照仪器安装规范方法，将HyMap安装在运12飞机上，保证机下的摄影窗口满足数据获取要求，按照航空成像光谱数据获取规范进行地面测试；最后为HyMap航空高光谱数据获取和初步质量检查。

在此基础上，开展HyMap航空高光谱遥感数据辐射定标、大气校正光谱反演、几何校正和地理编码等数据处理工作。通过此项研究，为高光谱遥感地空综合预测铜金矿床的技术方法研究和应用提供关键的数据基础。

4、地面高光谱数据获取及处理

地面光谱测量采用美国A S D公司生产的FieldSpec Pro FR地面光谱测量仪（表2）。利用该仪器开展了地面光谱准同步定标（图2）。在HyMap高光谱数据获取过程中，准同步布设黑白布定标场和明暗地物定标场，利用FieldSpec Pro FR地面光谱测量仪进行光谱定标数据获取（图2-3），并基于经验线性模型法，对高光谱影像进行大气校正和光谱比对分析。

同时开展了典型铜金矿床地表蚀变围岩和其他相关的岩石地物的光谱测量，完成了4条贯穿地层走向的成矿有利地质单元光谱测量和4条典型矿床光谱测量。共完成光谱测量点287个，光谱测量曲线2870条。测量完成后，开展了各岩石类型的岩矿光谱数据处理、光谱库建立、光谱特征分析等研究工作，研究及分析结果为航空高光谱蚀变矿物信息提取、地面光谱验证、铜金矿床围岩蚀变信息筛选等提供重要的地面光谱信息支持。

地面高光谱数据处理采用ViewSpec Pro软件，利用ViewSpec Pro软件挑选符合要求的地面岩矿光谱曲线，并将这些曲线转换为txt格式的文件。由于受到太阳位置、角度条件、大气条件、地形影响及传感器本身性能的影响，传感器所记录的数据与目标的光谱反射率或光谱辐射亮度值并不一致。需先将传感器记录的原始辐射值（DN值）转化为地物反射率[20]。

5、蚀变矿物信息提取处理方法及分布特征

目前基于高光谱遥感影像数据提取矿物的处理方法可以分为两大类：a）通过像元光谱与矿物标准光谱曲线的匹配实现识别和量化目标矿物；b）提取矿物光谱特征吸收波段参数作为识别和量化指标[21]。前者在处理过程中光谱数据信息可能由于重定标而部分丢失。且基于单幅图像统计参数的处理方法难以做到无缝衔接，该方法也涉及大量的运算。本次研究区蚀变矿物提取主要基于第二种方法。该方法利用ENVI软件，基于可编程实现的特征提取处理通道，该通道中脚本设计是基于提取光谱参数，首先对影像光谱数据进行均值标准化、去连续统、卷积平滑等数据预处理，在进行特征深度、特征面积、特征宽度、极小值波长、比值、算术和逻辑运算等特征参数提取。同时利用多个诊断性特征来标记或约束特定的矿物，提取矿物的分布和相对含量，并增强矿物识别的精度。

根据矿物光谱特征吸收波段参数的蚀变矿物提取方法，在研究区提取了褐铁矿-黄钾铁矾、白云母-蒙脱石、绿泥石-绿帘石-碳酸盐、叶腊石等蚀变矿物或组合（图3），其中褐铁矿-黄钾铁矾、白云母-蒙脱石、绿泥石-绿帘石-碳酸盐等蚀变矿物组合在全区范围内分布，且呈团块状、条带状及星散状。叶腊石则仅在研究区北东部呈星点状分布，蚀变不明显，对于指导找矿意义不大。

6、红山铜金矿床高光谱遥感地空综合找矿模型的建立与找矿有利地段优选

6.1 矿区地质特征

红山铜金矿床是高硫型-浅成低温热液型矿床[22]。矿床产于陆相蚀变酸性火山弯窿及其通道中，该次火山机构受近东西与北西向两组断裂控制，且主要控岩控矿构造为北西向断裂。铜矿体受北西向断裂及其次级构造破碎带控制，为浸染状辉铜矿化和黄铜矿化。金矿化主要与酸性次火山岩中的细粒黄铁矿有关[23]（图4）。

矿区地表发育黄钾铁钒化、绢云母化、硅化、褐铁矿化、黄铁矿化、绿泥石化、绿帘石化、青盘岩化。近矿蚀变以火山岩酸淋滤带蚀变体的黄钾铁矾化-明矾石化-高岭土化-石膏和中性斑岩体的绢云母化-硅化为显著特征，远矿蚀变为青盘岩化[24]。矿区中心为泥化、硫酸盐带，向北为硅化-伊利石化-绢云母化-褐铁矿化带。具有找矿标志的矿化蚀变为：硅化-绢云母化、伊利石化、褐铁矿化、硅化、高岭土化-明矾石化等[23]。

红山铜金矿床矿化以铜、金矿为主，赋矿岩性主要为流纹岩、英安岩、安山玢岩、闪长玢岩、火山角砾岩。矿石矿物主要为黄铜矿、蓝铜矿、闪锌矿、磁铁矿、赤铁矿、黄钾铁钒等。脉石矿物为石英、斜长石、绢云母、绿泥石等。矿石呈稠密浸染状、细脉状、角砾状构造[25]。

6.2 红山铜金矿床高光谱遥感地空综合找矿模型的建立

充分研究区域成矿有利的构造、地层、岩浆岩条件、变质条件的基础上，结合区域控矿因素、矿床、矿化蚀变特征，进行航空高光谱蚀变矿物分布特征分析及评价。本文截取了研究区内红山铜金矿床及外围基于HyMap航空高光谱遥感影像的蚀变矿物分布图（图5），图中显示红山铜金矿床中心位置呈黑白相间且颜色突变的色调，影纹粗糙且不均匀，矿床外围呈灰黑色、灰绿色色调，影纹粗糙。与矿床有关的航空高光谱蚀变矿物组合主要为褐铁矿+黄钾铁钒、绢云母+蒙脱石等，呈大面积的团块状分布。在矿床外围分布着呈星散状、星点状的绿泥石+绿帘石蚀变矿物组合，该组蚀变为铜金矿床外围主要发育的蚀变组合，对于寻找该类型矿床有一定的指示作用。

在红山铜金矿区收集了具体的剖线采样和光谱数据，图6为红山铜金矿区的地质光谱综合剖面，出露的地层为奥陶统荒草坡群大柳沟组，岩性主要为英安岩、片理化带、流纹斑岩、花岗岩、花岗闪长岩、构造蚀变岩以及铜金矿化体等。蚀变矿物主要有褐铁矿、黄钾铁钒、绢云母、石膏、高岭石、明矾石及方解石等。通过样品采集和光谱测试结果，研究与该区铜金矿关系密切的几种赋矿岩性的光谱曲线，分别为构造蚀变岩、花岗闪长岩、花岗岩、流纹斑岩、蚀变英安岩、及铜金矿化体。

经研究分析发现：红山铜金矿床及外围不同岩性地面光谱特征各不相同，各类岩性具体曲线特征如下（图7）：

构造蚀变岩（图7-a）地面光谱曲线特征显示：880nm处的Fe3+特征吸收可能为褐铁矿化所致；2210nm处的弱吸收峰为发育弱绢云母化所致，2260nm处的较强吸收峰为岩石发育黄钾铁矾化所致。故构造蚀变岩蚀变矿物组合为：褐铁矿+黄钾铁矾+绢云母。

花岗闪长岩（图7-b）地面光谱曲线特征显示：880nm处较强的对称吸收峰为岩石发育褐铁矿化所致；2210nm处较强的吸收峰，及右肩2250nm处较弱次级吸收峰，是Al-OH的特征吸收，为明矾石化所致；2350nm处强的主吸收峰，左肩2250nm处较强的次级吸收峰，是Mg-OH的特征吸收，为绿泥石化所致。故其蚀变矿物组合为：褐铁矿+明矾石+绿泥石。

花岗岩（图7-c）地面光谱曲线特征显示：910nm处较弱的Fe3+特征吸收峰为褐铁矿化所致；2210nm处的弱Al-OH特征吸收峰为绢云母化引起，其右肩2250nm处的较弱吸收峰可能为明矾石化所致；2265nm处的较弱吸收峰可能为黄钾铁矾化所致；2330nm处可见较强吸收峰，是CO32-的特征吸收峰，为岩石发育碳酸盐化所致，故花岗岩蚀变矿物组合为：褐铁矿+黄钾铁矾+绢云母+明矾石+方解石。

流纹斑岩（图7-d）地面光谱曲线特征显示：880nm处的较强吸收峰为褐铁矿化所致；2210nm处的强Al-OH特征吸收峰为岩石发育绢云母化所致；2265nm处的弱Fe-OH特征吸收峰为岩石发育黄钾铁矾化所致，故流纹斑岩蚀变矿物组合为：褐铁矿+黄钾铁钒+绢云母。

蚀变英安岩（图7-e）地面光谱曲线特征显示：880nm处的弱吸收峰为弱褐铁矿化所致；2210nm处的弱Al-OH特征吸收峰为绢云母化引起；2350nm处的较强主吸收峰及其左肩2250nm处的较强次级吸收峰，是Mg-OH的特征吸收峰，为绿泥石化所致；2330nm处的较强主吸收峰以及1880nm、2060nm、2170nm处较弱次级吸收峰，是CO32-的特征吸收峰，为岩石发育碳酸盐化所致；1836nm、1935nm、2213nm、2270nm处可见多处较强吸收峰，为岩石发育石膏所致。故蚀变英安岩蚀变矿物组合为：褐铁矿+绢云母+绿泥石+方解石+石膏。

铜金矿化体（图7-f）地面光谱曲线特征显示：880nm处的较强Fe3+特征吸收为岩石发育褐铁矿化所致；2210nm处的弱Al-OH特征吸收峰为岩石发育绢云母化；2265nm处的弱Fe-OH特征吸收峰为岩石发育黄钾铁矾化；1540nm、1835nm、1920nm、2000nm、2165nm、2208nm、2250nm、2329nm、2374nm、2385nm处可见多处较强吸收峰，可能为岩石发育黝帘石化所致。故铜金矿化体蚀变矿物组合为：褐铁矿+绢云母+黄钾铁钒+黝帘石。

根据航空高光谱和地面高光谱蚀变矿物或组合的类型和分布特征，结合红山铜金矿床的地质特征，进行典型矿床航空-地面高光谱综合剖析，最终建立了红山铜金矿床高光谱遥感地空综合找矿模型（表3）。

6.3 铜金矿床找矿有利地段优选

在建立红山铜金矿床高光谱遥感地空综合找矿模型的基础上，结合地质、区域化探、野外调查验证等多源地学信息，运用高光谱遥感地空综合预测方法，优选了两处与已知矿床航空-地面高光谱向类似的找矿有利地段，分别为I-1和I-2（图8）。

I-1找矿有利地段分布玉带铜金矿点北西侧一带，呈北北西向展布，面积约5.2km2。区内分布的地层为中-上志留统红柳峡组（S2-3h）火山岩，近东西向断裂构造发育，区内多发育褐铁矿化、硅化、绿帘石化等蚀变。化学分析显示钛含量最高达0.23×10-2，金红石含量最高达0.38×10-2。为铜、金、钛等元素多金属找矿有利地段。

I-2找矿有利地段分布在东二区南东侧，呈近东西向展布，面积约6.8km2。区内分布的地层为中奥陶统荒草坡群大柳沟组（O2Hd）火山岩，北东东向断裂构造发育，花岗闪长岩体呈近东西向展布，区内多发育硅化、褐铁矿化、赤铁矿化、绿帘石化、孔雀石化等蚀变。化学分析显示钛元素含量最高达0.21×10-2，金红石含量最高达0.35×10-2，且区内已有红山铜金矿床分布，为铜、金、钛等元素多金属找矿有利地段。

7、结论

（1）新疆东天山卡拉塔格铜金多金属矿床集中，地表具有指示意义的蚀变矿物分布较为典型，适合利用高光谱遥感地空综合预测方法在该区寻找铜金及多金属矿床，且找矿效果显著，故该方法为寻找铜金矿床的有效技术方法。

（2）在运用HyMap航空成像光谱仪和FieldSpec Pro FR光谱仪进行高光谱数据获取和处理的基础上，结合HyMap航空高光谱蚀变矿物提取效果分析，最终开展典型矿床航空-地面高光谱综合剖析研究，建立了红山铜金矿床高光谱遥感地空综合找矿模型，并综合地质、区域化探及野外调查验证等多源地学信息，圈定了I-1和I-2等两处找矿有利地段。

（3）与常规的地质、区域化探等地质矿产勘查方法相比，高光谱遥感地空综合预测方法具有低成本、快速高效的优势。但由于各类矿床的找矿标志和高光谱异常信息多种多样，从而造成蚀变异常信息的多解性，加上利用高光谱地空综合预测方法开展多金属矿产资源调查研究还比较少，需要今后继续完善成熟。

致谢：核工业航测遥感中心汪冰研究员在成文过程中给予了建设性的指导和帮助，野外工作得到王瑞军、李存金等同事的大力帮助，在此表示衷心的感谢！

（

）

[1]Sabins F F.Remote sensing for mineral exploration[J].OreGeology Reviews，1999，14：157-183.

[2]Everett J R，Jengo C J，Staskowski R J，et al.Remote and GIS enable futureexploration success[J].World Oil，2002，223(11)：59-65.

[3]连长云，章革，元春华，等.短波热红外光谱矿物测量技术在热液蚀变矿物填图中的应用-以土屋斑岩铜矿床为例[J].中国地质，2005，32(3)：483-495.

[4]廖崇高，杨武年，徐凌，等.成矿预测中遥感与地质异常的综合分析-以兰坪盆地为例[J].中国地质，2002，29(4)：416-420.

[5]刘颖璠，陈建平，郝俊峰，等.内蒙古朱拉扎嘎地区成矿遥感信息提取及成矿预测研究[J].中国地质，2012，39(4)：1062-1068.

[6]钱建平，伍贵华，陈宏毅.现代遥感技术在地质找矿中的应用[J].地质找矿论丛，2012，27(3)：355-360.

[7]谭克龙，万余庆，王晓峰，等.基于遥感技术的煤炭勘查方法研究[J].中国地质，2012，39(1)：218-227.

[8]束炯，王强，孙娟.高光谱遥感的应用研究[J].华东师范大学学报，2006(4)：1-10.

[9]Campbell，J.b.，R.H.Wyname.2011.Introduction to remote sensing [M].Fifth Edition.New York：Guilford Press：3-28.

[10]刘德长，邱骏挺，田丰，等.区域控矿断裂带的航空高光谱遥感技术研究-以黑石山-花牛山深大断裂带为例[J].地质与勘探，2015，51(2)：366-375.

[11]Clack，R.N.，T.V.V.King，M.Klejwa，G.A.Swayyze，N.Veego.1990.High spectral resolution reflectance spectroscopy of minerals [J].Journal of Geophysical Research，95：12653-12680.

[12]王润生，甘甫平，闫柏琨，等.高光谱矿物填图技术与应用研究[J].国土资源遥感，2010，21(1)：1-13.

[13]王润生，熊盛青，聂洪峰，等.遥感地质勘查技术与应用研究[J].地质学报，2011，85(11)：1700-1743.

[14]姚佛军.荒漠戈壁区金属矿床波谱特征和蚀变遥感异常信息提取研究及在矿产资源评价中的应用[D].北京：中国地质科学院，2006.

[15]Povarennykh A S.The Use of Infrared Spectra for the Determination of Minerals.Am.Mineral.， 1978，63：956-959.

[16]刘圣伟，甘甫平，闫柏琨，等.成像光谱技术在典型蚀变矿物识别和填图中的应用[J].中国地质，2006(2)：178-186.

[17]甘甫平，王润生，马蔼乃.基于特征谱带的高光谱遥感矿物谱系识别[J].地球科学，2003，10(2)：445-454.

[18]甘甫平，王润生，马蔼乃，等.遥感地质信息提取集成与矿物遥感地质分析模型[J].遥感学报，2003，7(3)：207-213.

[19]Clark R N，King T V V，Klejwa M et al.High spectral resolution reflectance spectroscopy of minerals[J].J. Geophys Res，1990，95:12653-12680.

[20]李志忠，汪大明，刘德长，等.高光谱遥感技术及资源勘查应用进展[J].地球科学，2015，40(8)：1287-1294.

[21]卢燕.机载HyMap高光谱数据获取和处理工作报告[R].2014.12.

[22]董承维.卡拉塔格矿区重磁特征及找矿预测[D].昆明：昆明理工大学，2013.

[23]王瑞军，牛海威，孙永彬，等.新疆重点地区航空高光谱调查与找矿预测技术研究阶段成果报告[R].2014.12.