基于中高分辨率数据的内蒙古...遥感找矿信息提取及找矿预测_连琛芹.pdf

第 37 卷第 1 期2023 年 2 月现代地质GEOSCIENCEVol.37No.1Feb.，2023DOI:10.19657/j.geoscience.1000 8527.2023.01.24基于中高分辨率数据的内蒙古红石山地区遥感找矿信息提取及找矿预测连琛芹1，2，封志明1，2，刘永新3，雷聪聪3(1.中国科学院 地理科学与资源研究所，北京100101;2.中国科学院大学，北京100049;3.中国地质调查局 呼和浩特自然资源综合调查中心，内蒙古 呼和浩特010010)收稿日期:2022-01-15;改回日期:2022-11-22;责任编辑:楼亚儿。基金项目:科技部第二次青藏高原综合科学考察研究项目(2019QZKK1006)。作者简介:连琛芹，男，博士研究生，1992 年出生，自然资源学专业，主要从事资源遥感研究工作。Email:lianchenqin1911 。通信作者:封志明，男，研究员，1963 年出生，自然地理学专业，主要从事资源开发与区域可持续发展的研究工作。Email:fengzm 。摘要:内蒙古红石山地区位于东天山兴安成矿带内，具有良好的成矿条件，但其位于荒漠戈壁，自然条件恶劣，基础地质研究程度低，因此利用遥感手段进行多金属找矿勘查具有重要意义。采用 OLI 数据，选择掩膜+无损线性拉伸+特征向量主成分变换+异常切割+异常滤波的方法，提取铁染、羟基等矿化异常信息;选择假彩色合成及线性拉伸，提取断裂构造信息。综合多元信息进行找矿预测，共圈定2 处找矿远景区。本次研究对其中1 处区域进行进一步重点勘查。重点勘查区域采用高分辨率 PLEIADES 数据，选择图像融合、彩色合成、直方图拉伸等方法，进一步圈定巨型石英脉体。在重点勘查区域进行野外查证，发现了多处矿化点。此次研究表明综合使用中分辨率 OLI 数据和高分辨率 PLEI-ADES 数据，可以大大提升地质找矿效率，对未来该地区进一步找矿勘查具有重要的指导意义。关键词:遥感找矿;红石山;OLI;PLEIADES;多金属矿中图分类号:P627;TP79文献标志码:A文章编号:1000 8527(2023)01 0227 06emote Sensing Prospecting Information Extraction and Mineral PredictionBased on Medium and High esolution Data in Hongshishan，Inner MongoliaLIAN Chenqin1，2，FENG Zhiming1，2，LIU Yongxin3，LEI Congcong3(1.Institute of Geographic Sciences and Natural esources esearch，Chinese Academy of Sciences，Beijing100101，China;2.University of Chinese Academy of Sciences，Beijing100049，China;3.Hohhot Natural esources Comprehensive Survey Center，China Geological Survey，Hohhot，Inner Mongolia010010，China)Abstract:The Hongshishan area of Inner Mongolia is located in the East Tianshan-Xingan metallogenic beltand has good metallogenic conditions However，poor natural conditions and low level of basic geological surveyrestrict mineral exploration in this area，so it is of great significance to use remote sensing to conduct polymetal-lic prospecting OLI data were used in this study，and the mineralization anomalies such as iron stain andhydroxyl group were extracted by mask，lossless linear stretching，eigenvector principal component transforma-tion，density segmentation，anomaly filtering，while the fracture information was extracted by pseudo-color syn-thesis and linear stretching Based on the multivariate information，two target areas were delineated and one ofthem was further explored High-resolution PLEIADES data，image fusion，color synthesis，histogram stretchingand other methods were used to delineate the giant quartz veins in key exploration areas On the basis of remotesensing interpretation，the field investigation proves that there are several mineralization points in the key explo-ration area This study shows that the combination of medium-resolution OLI data and high-resolution PLEIA-DES data can greatly improve the efficiency of geological prospecting，which has important guiding significancefor further prospecting in this areaKey words:remote sensing prospecting;Hongshishan;OLI;PLEIADES;polymetallic deposit0引言遥感技术已经广泛应用于地质找矿工作中1 5。国内外学者在利用遥感数据进行找矿勘查方面进行了大量而深入的研究，并取得了许多重要的成果6。不同的遥感数据都有其自身的优势和特点。Landsat 和 ASTE 卫星数据是地质找矿中应用最广泛的中分辨率数据，适合区域性遥感勘查。QuickBird、WorldView 及 Spot 等高分辨率卫星数据，适合大比例尺地质对象的解译和信息提取。综合使用中高分辨率遥感数据，可以最大限度地提取找矿信息，为矿产勘查提供有力支撑。内蒙古红石山地区自然条件恶劣，位于大漠戈壁，研究程度低，遥感找矿研究处于空白状态。因此利用遥感手段进行多金属找矿勘查，对本地区具有重要的意义。本文采用 OLI 数据，选择掩膜+无损线性拉伸+特征向量主成分变换+异常切割+异常滤波的方法，提取铁染、羟基等矿化异常信息;选择假彩色合成及线性拉伸，提取断裂构造信息。综合多元信息进行找矿预测，共圈定了 2 处找矿远景区，并对其中 1 处区域进行进一步重点勘查。重点勘查区采用高分辨率 PLEIA-DES 数据，选择图像融合、彩色合成、直方图拉伸等方法，进一步圈定巨型石英脉体。在重点勘查区域进行野外查证，发现了多处矿化点。此次研究获得了良好的找矿效果，表明综合使用中高分辨率遥感数据，可以大大提升本地区地质找矿效率，对未来该地区进一步找矿勘查具有重要的指导意义。1研究区概况研究区地理位置位于内蒙古自治区境内，行政区划隶属内蒙古阿拉善盟额济纳旗苏泊淖尔苏木管辖，地处中蒙边界，研究区完全位于大漠戈壁，遥感影像如图 1 所示。该区大地构造位于阴山天山东西向复杂构造带中段，处在构造带北侧近边缘部位。研究区内地层出露不全，主要以奥陶系咸水湖组为主，其次为蓟县系青白口系圆藻山群以及白垩系巴音戈壁组。除白垩系下统不整合覆于下伏地层之上外，其余各岩石地层单位之间皆为构造接触关系。区内地质构造的总体特点是，北西近东西向断裂构造为本区的主要构造，并伴随有北东向的褶皱构造。区内侵入岩广泛出露且岩石类型较多，从超基性岩到酸性岩均有出露。其中酸性岩分布广泛，出露面积大，岩石类型以二长花岗岩、花岗闪长岩为主，多以岩基、岩株产出。其次为中性岩，主要岩石类型为变质闪长岩、变质石英闪长岩、变质石英二长闪长岩，以岩株和小岩体产出。基性岩出露面积比较小，主要岩石类型有变质橄榄辉长岩、变质辉长岩，以岩株或岩瘤产出。区内的脉岩也较发育。图 1研究区遥感影像图Fig.1emote sensing image of the study area2遥感数据基本特征2.1OLI 数据Landsat 8 卫星 OLI 数据与 Landsat 7 卫星数据相比，Landsat 8 卫星数据很好地延续了 Landsat 系列卫星的一贯特点，同时 Landsat 8 的波段更多，波段划分更加精细7。Landsat 8 卫星装载了陆地成像仪(OLI)和热红外传感器(TIS)。本次使用的是 OLI 数据，该数据与 ETM 数据相比，增加了光谱范围和信噪比，具有 12 位的量化级数，因此OLI 传感器展现出了更加良好的性能8。2.2PLEIADES 数据Pleiades 1 和 Pleiades 2 相同的两颗卫星组成Pleiades 高分辨率卫星星座。两颗卫星分别于 2011822现代地质2023 年年和 2012 年先后发射，并开始获取遥感数据。Pleiades 卫星属于军民两用光学成像卫星，其最高空间分辨率已达 0.5 m。双星配合每日可实现全球的重访，可以满足不同客户的要求。该数据包含 4个多光谱波段和 1 个全色波段，空间分辨率分别为 2 m 和 0.5 m，幅宽为 20 km。PLEIADES 遥感数据凭借其高分辨率的优势，在地球科学领域的应用潜力十分巨大，但 PLEIADES 数据在地质找矿领域的应用鲜有报道。本文充分发挥 PLEIADES数据的独特优势，服务于地质找矿工作。3矿化蚀变信息提取3.1波谱特征分析通常情况下，在遥感获得的地表信息中，只要有一定强度的蚀变岩石出露，就可以在遥感影像上有所反映6。根据研究区所提供的地质矿产资料，与成矿有关的围岩蚀变主要为硅化、绿泥石化、绿帘石化、绢云母化，次为褐铁矿化、黄钾铁矾化等。这些蚀变矿物多为含有羟基和铁离子的矿物，在 Landsat 8 OLI 的波段上具有诊断性波谱特征，为遥感蚀变信息的提取提供了理论支撑。特定的波段形成了识别特定矿物的特征波谱，蚀变矿物的特征波谱是遥感蚀变信息提取的波谱前提9。在本次蚀变矿物信息提取研究中，将USGS 波谱库中的蚀变矿物波谱以 OLI 遥感数据为标准进行了波谱重采样(图 2)。根据含铁离子的蚀变矿物的波谱曲线在 OLI 影像数据的各通道对应关系可知，在 OLI 数据的 2 波段反射率很低，具有强吸收特征，在 4 波段的反射率相对 2 波段高，在第 5 波段具有吸收谷，第 6 波段具有强反射峰的诊断性波谱特征。根据含羟基的蚀变矿物的波谱曲线在 OLI 影像数据的各通道对应关系可知，在 OLI 数据的第 5 波段形成吸收谷，在第 6图 2蚀变矿物波谱曲线Fig.2Spectral curves of the altered minerals波段具有反射峰，第 7 波段具有吸收谷的诊断性波谱特征。3.2矿化蚀变信息提取一般常见的遥感蚀变信息提取的技术方法有比值变换、主成分变换、光谱角、混合像元分解和 MPH 法等，其中主成分变换技术是应用最为广泛，具有对影像大气校正质量要求不高、实现简单、提取效果好、效果稳健等优点，在地质找矿实践过程中，广受地质工作者的青睐10 11。主成分变换又称 K L 变换，此方法变换前后遵循能量守恒定律，即运算前和运算后的信息总量保持不变，其对影像进行集中和压缩，将各通道的高度相关的信息集中在少数几个波段中，使变换后的各主分量相互独立，互不相关。OLI 遥感影像在进行蚀变信息提取前的预处理主要包括研究区裁剪和干扰地物剔除等。植被、云雪、水体、阴影等干扰，对异常的提取有重要影响，会引起大量的假异常，因此对干扰地物的去除十分重要。在本次研究中，我们选择了干扰地物很少的遥感数据，只有少量的阴影，采取掩膜的方法将其剔除。根据与研究区多金属成矿有关的蚀变矿物的光谱诊断特征，选择 Landsat 8 OLI 数据的 2、4、5、6 波段进行主成分分析，提取研究区与铁染有关的遥感蚀变信息。根据含铁矿物的光谱特征，可知构成铁染主分量的波段 2、波段 4 的系数符号通常是相反的，且波段 2 的值为负值。表 1 列出了主成分分析之后的铁染各主分量的物理参数。根据表中铁染的 4 个波段的主成分分析后的各自的特征向量矩阵可知，遥感蚀变信息位于 PC4 中。表 1OLI 数据主成分分析特征向量矩阵表Table 1The eigenvectors of covariance of PCA of OLI主成分B2B4B5B6PC10.143370.412590.412590.70186铁染PC20.374730.447160.43340.68686PC30.873720.00970.44890.18711PC40.275020.793550.542260.02441主成分B2B5B6B7PC10.125410.502380.651710.55422羟基PC20.332840.731830.135840.57896PC30.436870.132350.677030.57729PC40.826220.441050.313790.15614选择 Landsat 8 OLI 数据的 2、5、6、7 波段进行主成分分析，提取研究区与羟基有关的遥感蚀922第 1 期连琛芹等:基于中高分辨率数据的内蒙古红石山地区遥感找矿信息提取及找矿预测图 3蚀变异常信息图Fig.3Extraction results of the alteration anomalies information变信息。根据表 1 中含羟基矿物的各主分量的判断标准和提取蚀变信息的主分量的物理参数可知，PC4 为羟基信息的主分量。根据含羟基矿物的光谱特征，构成羟基主分量的波段 2、波段 6 的系数符号是相同的，且均为正值。因此需将 PC4 乘以(1)处理，使得波段 6 的值为正值。3.3蚀变信息成图研究区的遥感蚀变信息需进行异常切割，通过对“均值+n 标准差”中 n 值的调整，进行蚀变信息的量化分级。在本次遥感铁染蚀变信息级别的划分过程中，n 值取 1.3、1.5、1.7，分别对应一、二、三级蚀变异常。在本次遥感羟基蚀变信息级别的划分过程中，n 值取 1.4、1.7、2.0，分别对应一、二、三级蚀变异常。对异常切割后的图像进行滤波以消除因噪声引起的孤立的“伪”异常。最后将铁染和羟基信息的不同等级赋予不同的颜色，得到研究区遥感矿化蚀变信息图(图 3)。从遥感矿化蚀变信息分布图发现，提取的遥感蚀变信息中，羟基异常信息占主导地位，铁染异常信息较少，与已知地质矿产资料吻合，同时，其主要呈北西向和近东西向展布，且部分呈大面积块状或条带状分布，可能与构造活动有关，这些地方可能具有一定的找矿意义。4遥感构造信息提取地质构造的解译是根据各类地物波谱特征，以及所表现的颜色和纹理特征，通过增强图像中的地物信息的差异，突出感兴趣区的地物信息，进行地质构造解译12 13。本次研究通过 ENVI 软件对 OLI 数据进行假彩色合成及图像线性拉伸等图像增强处理，增大不同地质体之间的差异，从而更加容易进行地质构造的解译工作(图 4)。构造活动有利于成矿元素的活化和成矿流体的运移聚集，对矿床形成具有重要作用。构造密集和交汇部位是成矿作用的有利部位，可为找矿勘查提供构造基础。从解译出的构造可以发现，研究区以北西向构造为主，近东西向、北东向构造为辅，且在其中部发现大型北西向构造与遥感蚀变带吻合好，说明该区域的构造和蚀变信息具有高度相关性，进一步印证蚀变信息和构造信息提取的有效性，为后续找矿远景区圈定奠定基础。5多元信息综合分析及找矿远景区优选矿产勘查就是找矿信息逐渐获取和找矿风险逐渐减小的过程，在地质找矿过程中，各类信息都有其独特优势，应充分发挥各类信息针对性的作用14。本次研究将多元信息叠合于 GIS 平台，进行分析验证。在 GIS 平台上，综合遥感、地质和地球化学等多元信息，进行找矿远景区优选，共圈定 2 处找矿远景区(图 5)。其中，1 号找矿远032现代地质2023 年图 4构造解译图Fig.4Structural interpretation map图 5多元信息叠加图Fig.5Multi-source information overlay map景区位于工作区的西北部，北西向构造发育，化探数据显示该区位于综合化探异常区，遥感蚀变信息以羟基信息为主，呈北西向条带状沿构造展布，各类找矿信息吻合度高。2 号找矿远景区位于工作区西南部，北东向及近东西向构造发育，化探数据显示其位于综合化探异常区附近，遥感蚀变信息为羟基信息，呈带状和团块状展布于构造附近，遥感蚀变信息受构造控制明显。两处找矿远景区均位于遥感蚀变信息和构造信息的周围，同时与综合化探异常吻合较好，具有一定的找矿意义。6重点区域找矿勘查及野外查证本次研究对其中的 2 号远景区进行进一步勘查，使用了该区域的高分辨率 PLEIADES 数据。本次研究为了充分利用高分辨率数据的优势，将全色和多光谱波段进行融合处理，得到可见光和近红外波段的 0.5 m 分辨率的遥感数据。采用彩色合成和直方图拉伸等方法进行图像增强，以解译石英脉体。本次研究在重点勘查区进行野外查证，发现了多处矿化点，它们具有不同强度的硅化、绿泥石化、绿帘石化、绢云母化、褐铁矿化等特征。图 6a 是重点勘查区的 1 处巨型石英脉体露头，其位于蚀变信息图中北西向和近东西向的蚀变带中。在整个重点勘查区中，展布了许多大小不一的巨型石英脉体。图 6b 是蚀变样品，可以看见褐铁矿化和绿泥石化等现象。野外验证工作进一步表明本次研究是有效的，对未来该区及其附近区域的找矿勘查具有重要的指导意义。132第 1 期连琛芹等:基于中高分辨率数据的内蒙古红石山地区遥感找矿信息提取及找矿预测图 6巨型石英脉体与蚀变样品Fig.6Giant quartz vein body and altered sample7结语本文在内蒙古红石山地区，利用中分辨率Landsat 8 OLI 数据和高分辨率 PLEIADES 数据、各类影像处理方法、地学数据综合分析，进行遥感综合找矿工作，得到以下结论:(1)选择 Landsat 8 OLI 数据，采用掩膜+无损线性拉伸+特征向量主成分变换+异常切割+异常滤波的方法，成功提取了内蒙古红石山地区的铁染和羟基遥感蚀变信息。提取的遥感矿化蚀变信息以羟基蚀变信息为主导，铁染为辅。(2)选择 Landsat 8 OLI 数据，进行假彩色合成及线性拉伸，有效地解译了区内北西向大型主断裂和北西向及北东向次级断裂，确定了区内构造格架，对构造与矿化关系研究具有重要意义。(3)综合分析遥感找矿信息，发现呈北西向和近东西向的羟基信息展布于大型断裂构造周围，说明该异常信息的展布明显受控于区内断裂构造，同时与综合化探异常吻合较好。各类信息相互吻合度高，更进一步说明本次提取的各类信息的有效性。(4)综合多元信息进行找矿远景区优选，共圈定2 处找矿远景区。利用高分辨率 PLEIADES 数据对其中 1 处作为重点区进行进一步研究，并进行了野外查证，发现多处矿化点，表明本次研究对该地区未来找矿具有重要的指导意义，为后续持续找矿勘查提供了方向。参考文献:1LIU L，ZHOU J，JIANG D，et al Mineral resources prospectingby synthetic application of TM/ETM+，quickbird and hyperiondata in the Hatu area，West Junggar，Xinjiang，China JScientific eports，2016，6(1):21851 2尹芳，刘磊，张继荣，等 新疆谢米斯台地区小岩体型矿化遥感探测 J 地球学报，2014，35(5):561 566 3赵玉梅，刘磊，周军，等 新疆乌什塔拉地区遥感、化探、地质综合找矿研究 J 地质找矿论丛，2016，31(1):135 141 4王军年，刘磊，周军，等 基于遥感的新疆西准哈图地区宝贝断裂的发现及其找矿意义 J 地质找矿论丛，2018，33(1):100 107 5刘磊，周军，冯敏，等 甘肃北山辉铜山地区镁铁岩体遥感识别方 法 研 究 J遥 感 技 术 与 应 用，2013，28(3):520 525 6连琛芹，姚佛军，杨建民，等 半裸露区遥感蚀变信息提取研究 以甘肃玛曲地区为例 J 现代地质，2019，33(5):1079 1085 7初庆伟，张洪群，吴业炜，等 Landsat 8 卫星数据应用探讨 J 遥感信息，2013，28(4):110 114 8唐超，周可法，张楠楠，等 基于 Landsat 8 OLI 和 ASTE数据集成和融合的矿化蚀变信息提取:以包古图斑岩型铜矿为例 J 地质科技情报，2018，37(6):217 223 9刘知，周萍 基于 ASTE 数据的甘肃花牛山地区蚀变信息提取 J 中北大学学报(自然科学版)，2017，38(2):124 132 10 吴志春，叶发旺，郭福生，等 主成分分析技术在遥感蚀变信息提取中的应用研究综述 J 地球信息科学，2018，20(11):1644 1656 11 张远飞，袁继明，杨自安，等 基于物理意义的二维散点图类型划分与遥感蚀变信息提取 J 国土资源遥感，2013，25(2):57 62 12 赵小星，钱建平，张渊，等 河南祁雨沟金矿及外围高植被区遥感找矿预测 J 大地构造与成矿学，2012，36(4):642 651 13 耿新霞，杨建民，姚佛军，等 新疆阿巴宫铅锌矿遥感找矿信息提取及地化多元信息综合分析 J 地质与勘探，2010，46(5):942 952 14 耿新霞，杨建民，姚佛军，等 新疆阿勒泰阿巴宫铁矿遥感找矿综合信息研究 J 地质论评，2010，56(3):365 373232现代地质2023 年