一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法技术

本发明专利技术属于地球科学技术领域，具体涉及一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法。本发明专利技术利用地面成像高光谱数据光谱分辨率的优势和图谱合一的特点，通过采集钻孔岩芯高光谱图像数据，识别并提取出与热液成矿作用关系密切的蚀变矿物，同时定量估算钻孔深度方向上蚀变矿物的空间分带。本发明专利技术既能够准确提取与深部热液成矿关系密切的蚀变矿物信息，直观表现蚀变矿物在岩芯上的发育，又能够实现钻孔尺度上蚀变矿物的定量化表达，对深部找矿勘查具有重要的实用价值。

A method for hyperspectral extraction and quantitative estimation of altered minerals in borehole cores

The invention belongs to the field of Earth Science and technology, in particular to a hyperspectral extraction and quantitative estimation method for drilling core alteration minerals. The invention uses the ground imaging hyperspectral data resolution and map integration features, through the acquisition of core of hyperspectral image data to identify and extract and hydrothermal mineralization is closely related to the altered minerals, and quantitative estimation of drilling depth alteration mineral zoning. The invention can accurately extract and deep hydrothermal mineralization is closely related to the alteration mineral information, visual performance of altered minerals in core development, but also to achieve quantitative drilling scale alteration expression of deep ore prospecting has important practical value.

【技术实现步骤摘要】
一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法
本专利技术属于地球科学
，具体涉及一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法。
技术介绍
深部找矿是当前亟需突破的找矿方向，钻孔岩芯是研究深部成矿环境的重要载体，不同种类的蚀变矿物是热液矿床的重要识别标志和成矿有利因素。与地表或近地表不同，深部岩石中发育的变质类矿物大都是热液成因而非风化成因，因此，如何准确识别和提取钻孔岩芯中发育的与成矿密切相关的各种蚀变矿物及其组合，查明它们在深部空间的分带状况，对深部找矿勘查具有重要的现实意义。当前，对钻孔岩芯蚀变矿物的传统研究方法主要通过地质人员的钻孔岩芯编录实现，编录的同时对岩芯进行取样，然后进行送样化学分析，这种方法费时费力，且花费较大。同时，由于不同地质编录人员的知识储备和经验差异，容易得出不同的结果导致影响下一步的分析和结论。高光谱技术在识别矿物方面具有独特的优势。当前，高光谱遥感技术用于地表矿物填图已较为成熟，已逐渐广泛用于大面积地质调查中，基于不同矿物在不同波长位置的诊断性光谱特征和光谱匹配算法，能够准确、客观地将高光谱图像中的矿物提取出来，目前可识别的矿物种类多达十几种。基于上述技术背景，如何将高光谱技术用于钻孔岩芯蚀变矿物提取以为深部找矿勘查提供新的技术手段和思路，如何将钻孔岩芯蚀变矿物的发育状况和空间分布特征更加直观定量地呈现给地质勘探人员以便更好地研究地下深部热液蚀变分带和流体演化规律，是当今地球科学
亟需解决的技术问题。
技术实现思路
本专利技术需要解决的技术问题为：提出一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，利用地面成像高光谱数据光谱分辨率的优势和图谱合一的特点，通过采集钻孔岩芯高光谱图像数据，识别并提取出与热液成矿作用关系密切的蚀变矿物，同时定量估算钻孔深度方向上蚀变矿物的空间分带，为深部地质环境研究与找矿勘查提供技术支撑。本专利技术的技术方案如下所述：一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：包括以下步骤：步骤1.获取钻孔岩芯高光谱图像数据；步骤2.钻孔岩芯高光谱图像数据预处理；步骤3.建立钻孔岩芯蚀变矿物光谱匹配库；步骤4.提取蚀变矿物信息；步骤5.设定钻孔岩芯蚀变矿物估算尺度；步骤6.定量估算蚀变矿物含量；步骤7.钻孔岩芯蚀变矿物估算结果编录可视化。步骤1中，每隔一段时间，利用地面高光谱传感器设备，同步采集获取钻孔岩芯高光谱图像数据和标准板数据。作为优选方案，步骤1中，数据采集前岩芯箱内各岩芯展布方向平行；采集过程中，保持光照条件稳定，同时保持岩芯箱内岩芯展布方向与地面高光谱传感器设备的图像扫描方向平行。步骤2包括以下步骤：步骤2.1.辐射校正基于步骤1所采用的地面高光谱传感器设备定标参数文件，对采集的钻孔岩芯高光谱图像数据以及同步采集的标准板图像数据进行系统级辐射校正，得到辐射亮度图像数据；步骤2.2.光谱重建获取标准板图像的辐射亮度均值光谱，建立标准板辐射亮度均值光谱与标准板反射率光谱之间的经验线性公式，利用所述经验线性公式对辐射校正后的钻孔岩芯高光谱图像进行反射率计算，得到钻孔岩芯反射率图像数据。步骤3包括以下步骤：采集钻孔岩芯蚀变矿物样本的实测光谱曲线，或从已有标准矿物光谱库中挑选钻孔岩芯中可能存在的蚀变矿物标准光谱曲线，建立钻孔岩芯蚀变矿物光谱匹配库。步骤4包括以下步骤：步骤4.1.波段重采样以步骤3建立的钻孔岩芯蚀变矿物光谱匹配库为标准，对步骤2预处理后的钻孔岩芯高光谱图像数据进行波段重采样，使二者波谱范围一致；步骤4.2.光谱匹配计算利用高光谱匹配算法对波段重采样后的钻孔岩芯高光谱图像数据进行蚀变矿物光谱匹配计算，获得匹配计算后的图像数据；所述高光谱匹配算法可以采用光谱角匹配算法、光谱相关系数匹配算法，以及匹配滤波算法。步骤4.3.设定阈值提取蚀变矿物信息针对每一种蚀变矿物匹配计算后的图像数据，设定各自的统一提取阈值，提取全孔钻孔岩芯图像上发育的各种蚀变矿物信息。步骤5中，以钻孔深度为标准，根据具体应用需求，建立统一的估算尺度。步骤6中，针对每一估算尺度深度段的岩芯，利用数字图像像元统计方法统计计算每种蚀变矿物所占含量。所述数字图像像元统计方法如下所述：对于第m个估算尺度深度段的岩芯，获取第m个估算尺度深度段岩芯内的所有像元数Nm；在第m个估算尺度深度段岩芯内，获取提取的第i种蚀变矿物像元数nmi，则在第m个估算尺度深度段岩芯内第i种蚀变矿物的含量为(nmi/Nm)×D，其中D表示步骤5设定的估算尺度。步骤7中，针对全孔岩芯，根据钻孔深度顺序，对各估算尺度内钻孔岩芯蚀变矿物估算量进行编录，利用统计制图软件对编录结果进行可视化表达。本专利技术的有益效果为：本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，既能够准确提取与深部热液成矿关系密切的蚀变矿物信息，更加直观地表现蚀变矿物在岩芯上的发育，又能够实现钻孔尺度上蚀变矿物的定量化表达，对深部找矿勘查具有重要的实用价值，具体具有以下有益效果：(1)本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，充分利用高光谱图像图谱合一的特点，基于地面成像高光谱数据，通过采集、处理钻孔岩芯的地面高光谱图像数据，利用高光谱匹配填图算法从光谱特征角度识别提取钻孔岩芯中发育的蚀变矿物信息，提取速度快、准确性高；(2)本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，与传统岩芯编录方法相比，通过钻孔岩芯高光谱蚀变信息提取，识别出的矿物信息更加丰富；(3)本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，划定钻孔深度尺度进行蚀变矿物的定量估算，实现了钻孔深度空间上的岩芯蚀变矿物定量化估算，有利于开展进一步的深入分析；(4)本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，提取精度高，估算结果与实际一致性好，对深部找矿勘查具有一定的现实意义。附图说明图1为本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法流程图；图2为某铀矿区钻孔蚀变矿物定量估算编录图。具体实施方式下面结合附图和实施例对本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法进行详细说明。如图1所示，本专利技术的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，包括以下步骤：步骤1.获取钻孔岩芯高光谱图像数据按照统一的数据采集标准，利用Hyspex或其他类似的地面高光谱传感器设备采集获取钻孔岩芯高光谱图像数据。数据采集过程中，每隔一段时间，同步采集标准板数据，用于步骤2进行岩芯高光谱图像数据的反射率光谱计算。作为优选实施方式，数据采集前，整理岩芯箱中的岩芯摆放，使各岩芯展布方向大致平行；采集过程中，保持光照条件稳定，同时尽可能保持岩芯箱中的岩芯展布方向与地面高光谱传感器设备的图像扫描方向平行，为后续数据处理提供便利，提高处理效率；采集的数据文件严格依据钻孔号、钻孔深度信息命名，防止由于数据量较多而产生文件名的混淆，避免给后续数据处理造成不便。步骤2.钻孔岩芯高光谱图像数据预处理步骤2.1.辐射校正基于步骤1所采用的地面高光谱传感器设备定标参数文件，对采集的钻孔岩芯高光谱图像数据以及同步采集的标准板图像数据进行系统级辐射校正，得到辐射亮度图像数据。所述辐射校正为本领域技术人员公知常识，可以通过本文档来自技高网...

【技术保护点】
一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：包括以下步骤：步骤1.获取钻孔岩芯高光谱图像数据；步骤2.钻孔岩芯高光谱图像数据预处理；步骤3.建立钻孔岩芯蚀变矿物光谱匹配库；步骤4.提取蚀变矿物信息；步骤5.设定钻孔岩芯蚀变矿物估算尺度；步骤6.定量估算蚀变矿物含量；步骤7.钻孔岩芯蚀变矿物估算结果编录可视化。

【技术特征摘要】
1.一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：包括以下步骤：步骤1.获取钻孔岩芯高光谱图像数据；步骤2.钻孔岩芯高光谱图像数据预处理；步骤3.建立钻孔岩芯蚀变矿物光谱匹配库；步骤4.提取蚀变矿物信息；步骤5.设定钻孔岩芯蚀变矿物估算尺度；步骤6.定量估算蚀变矿物含量；步骤7.钻孔岩芯蚀变矿物估算结果编录可视化。2.根据权利要求1所述的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：步骤1中，每隔一段时间，利用地面高光谱传感器设备，同步采集获取钻孔岩芯高光谱图像数据和标准板数据。3.根据权利要求2所述的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：步骤1中，数据采集前岩芯箱内各岩芯展布方向平行；采集过程中，保持光照条件稳定，同时保持岩芯箱内岩芯展布方向与地面高光谱传感器设备的图像扫描方向平行。4.根据权利要求2所述的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：步骤2包括以下步骤：步骤2.1.辐射校正基于步骤1所采用的地面高光谱传感器设备定标参数文件，对采集的钻孔岩芯高光谱图像数据以及同步采集的标准板图像数据进行系统级辐射校正，得到辐射亮度图像数据；步骤2.2.光谱重建获取标准板图像的辐射亮度均值光谱，建立标准板辐射亮度均值光谱与标准板反射率光谱之间的经验线性公式，利用所述经验线性公式对辐射校正后的钻孔岩芯高光谱图像进行反射率计算，得到钻孔岩芯反射率图像数据。5.根据权利要求4所述的一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法，其特征在于：步骤3包括以下步骤：采集钻孔岩芯蚀变矿物样本的实测光谱曲线，或从已有标准矿物光谱库中挑选钻孔岩芯中可能存在的蚀变矿物标准光谱曲线，建立钻孔岩芯蚀变矿物光谱匹配库。6.根据权利要求5所述的一种...

【专利技术属性】
技术研发人员：张川，叶发旺，徐清俊，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11