CN117828397B - Tunnel rock mineral identification method and system with cooperative multi-element spectrum - Google Patents

Tunnel rock mineral identification method and system with cooperative multi-element spectrum Download PDF

Info

Publication number
CN117828397B
CN117828397B CN202311816245.7A CN202311816245A CN117828397B CN 117828397 B CN117828397 B CN 117828397B CN 202311816245 A CN202311816245 A CN 202311816245A CN 117828397 B CN117828397 B CN 117828397B
Authority
CN
China
Prior art keywords
rock
tunnel
mineral
identification
spectrum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202311816245.7A
Other languages
Chinese (zh)
Other versions
CN117828397A (en
Inventor
许振浩
韩涛
周振发
林鹏
余腾飞
刘福民
邵瑞琦
李珊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN202311816245.7A priority Critical patent/CN117828397B/en
Publication of CN117828397A publication Critical patent/CN117828397A/en
Application granted granted Critical
Publication of CN117828397B publication Critical patent/CN117828397B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3563Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提出协同多元光谱的隧道岩石矿物识别方法及系统,涉及岩石矿物定性定量识别技术领域。包括采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库;获取隧道岩石光谱数据;利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配;利用岩石光谱的特征峰谱带位置和特征峰特征,基于端元矿物光谱库对矿物种类进行精细判识;建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;选取强相关性波段,并不断进行解混模型的迭代实现矿物的准确定量反演,完成定量识别。本发明实现了多元光谱协同的矿物定性定量识别。

The present invention proposes a tunnel rock and mineral identification method and system based on coordinated multivariate spectroscopy, which relates to the technical field of rock and mineral qualitative and quantitative identification. It includes collecting rock samples in the tunnel site area, obtaining accurate mineral identification results, establishing a rock spectrum library and an end-member mineral spectrum library in the tunnel site area; obtaining tunnel rock spectrum data; using the overall waveform of the fused spectrum data to perform a preliminary match with the data in the rock spectrum library in the tunnel site area; using the characteristic peak band position and characteristic peak characteristics of the rock spectrum to perform a fine identification of the mineral type based on the end-member mineral spectrum library; establishing an association between spectral variation characteristics and mineral types, and performing auxiliary identification of the tunnel rock mineral types that cannot be determined by fine identification, and completing the hierarchical qualitative identification of the tunnel rock minerals; selecting strongly correlated bands, and continuously iterating the unmixing model to achieve accurate quantitative inversion of minerals, and completing quantitative identification. The present invention realizes the qualitative and quantitative identification of minerals based on coordinated multivariate spectroscopy.

Description

协同多元光谱的隧道岩石矿物识别方法及系统Tunnel rock and mineral identification method and system based on collaborative multivariate spectroscopy

技术领域Technical Field

本发明属于隧道内岩石矿物定性定量识别技术领域,尤其涉及协同多元光谱的隧道岩石矿物识别方法及系统。The present invention belongs to the technical field of qualitative and quantitative identification of rocks and minerals in tunnels, and in particular to a method and system for identifying rocks and minerals in tunnels using a coordinated multivariate spectrum.

背景技术Background technique

本部分的陈述仅仅是提供了与本发明相关的背景技术信息,不必然构成在先技术。The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.

近些年,我国的隧道建设逐步向地质情况复杂的地区挺进。发明人发现,复杂的地质情况常导致断层、岩溶、蚀变带等不良地质的发育,若不加以及时判识,将极有可能导致地质灾害的产生,如突水突泥、塌方、围岩大变形等,严重影响现场施工安全和进度。In recent years, my country's tunnel construction has gradually moved to areas with complex geological conditions. The inventors have found that complex geological conditions often lead to the development of unfavorable geological conditions such as faults, karst, and alteration zones. If they are not identified in time, they are likely to cause geological disasters, such as sudden water and mud, landslides, and large deformation of surrounding rocks, which seriously affect the safety and progress of on-site construction.

隧道围岩的矿物异常现象常被用作断层等不良地质识别的标志,并可以识别不良地质的性质,已在相关领域中得到广泛应用。然而现有的矿物测试方法主要有镜下薄片法、X射线衍射、红外光谱等方法,镜下薄片和X射线衍射方法需要复杂的样品制备(磨粉或切片),且镜下薄片测试对环境要求高,隧道内的灰尘和黑暗极大影响判识结果准确性。因此,上述方法无法适应现场快速施工要求。The mineral anomalies of tunnel surrounding rocks are often used as signs of faults and other unfavorable geological conditions, and can identify the nature of unfavorable geology, and have been widely used in related fields. However, existing mineral testing methods mainly include microscopic thin section method, X-ray diffraction, infrared spectroscopy and other methods. Microscopic thin section and X-ray diffraction methods require complex sample preparation (grinding or slicing), and microscopic thin section testing has high environmental requirements. Dust and darkness in the tunnel greatly affect the accuracy of the identification results. Therefore, the above methods cannot meet the requirements of rapid on-site construction.

红外光谱法测试速度快、环境适应性强,可在隧道恶劣环境中获取围岩的光谱信息,但现有光谱的解译方法精度低、效率差,无法通过围岩的光谱特征准确识别矿物种类和反演矿物含量,难以对实际工程施工提供技术指导。Infrared spectroscopy has fast testing speed and strong environmental adaptability, and can obtain spectral information of surrounding rocks in the harsh environment of tunnels. However, the existing spectral interpretation methods have low accuracy and poor efficiency. They cannot accurately identify mineral types and invert mineral content through the spectral characteristics of surrounding rocks, making it difficult to provide technical guidance for actual engineering construction.

发明内容Summary of the invention

为克服上述现有技术的不足,本发明提供了协同多元光谱的隧道岩石矿物识别方法及系统,实现了多元光谱协同的矿物定性定量识别,对于提高隧道围岩矿物识别的效率和精度以及保障隧道快速安全施工具有重要意义。In order to overcome the shortcomings of the above-mentioned prior art, the present invention provides a tunnel rock mineral identification method and system based on collaborative multivariate spectroscopy, which realizes the qualitative and quantitative identification of minerals based on collaborative multivariate spectroscopy, and is of great significance for improving the efficiency and accuracy of tunnel surrounding rock mineral identification and ensuring rapid and safe tunnel construction.

为实现上述目的,本发明的一个或多个实施例提供了如下技术方案:To achieve the above objectives, one or more embodiments of the present invention provide the following technical solutions:

本发明第一方面提供了协同多元光谱的隧道岩石矿物识别方法。The first aspect of the present invention provides a tunnel rock mineral identification method based on collaborative multivariate spectroscopy.

协同多元光谱的隧道岩石矿物识别方法,包括以下步骤:The tunnel rock mineral identification method based on collaborative multivariate spectroscopy includes the following steps:

采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库,并获取岩石中不同矿物混合后的光谱变异特征;Collect rock samples in the tunnel site area, obtain accurate mineral identification results, establish a rock spectral library and an end-member mineral spectral library in the tunnel site area, and obtain the spectral variation characteristics of different minerals mixed in the rock;

沿隧道掘进方向获取隧道岩石光谱数据,将两个波段的光谱数据进行融合,利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配,确定隧道岩石的大致矿物组合情况;Acquire tunnel rock spectral data along the tunnel excavation direction, fuse the spectral data of the two bands, and use the overall waveform of the fused spectral data to make a preliminary match with the data in the rock spectral library of the tunnel site to determine the approximate mineral composition of the tunnel rock;

利用岩石光谱的特征峰谱带位置和特征峰精细特征与组合特征,基于端元矿物光谱库对矿物种类进行精细判识;Using the characteristic peak band position of rock spectra and the fine and combined characteristics of characteristic peaks, the mineral types are precisely identified based on the end-member mineral spectral library.

建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;Establish the association between spectral variation characteristics and mineral types, assist in identifying the types of tunnel rock minerals that cannot be determined through fine identification, and complete the hierarchical qualitative identification of tunnel rock minerals;

对隧道岩石光谱解混,利用端元矿物的纯净光谱逆向反推,选取强相关性波段,利用强相关性波段重新进行解混迭代,直至强相关性波段数量不变,完成对矿物的定量识别。The tunnel rock spectrum is unmixed, and the pure spectrum of the end-member mineral is used for reverse inference. The strongly correlated bands are selected and the unmixing is iterated again using the strongly correlated bands until the number of strongly correlated bands remains unchanged, thus completing the quantitative identification of the minerals.

本发明第二方面提供了协同多元光谱的隧道岩石矿物识别系统。The second aspect of the present invention provides a tunnel rock mineral identification system based on collaborative multivariate spectroscopy.

协同多元光谱的隧道岩石矿物识别系统,包括:The tunnel rock and mineral identification system based on collaborative multi-spectral analysis includes:

光谱库建立模块,被配置为:采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库,并获取岩石中不同矿物混合后的光谱变异特征;The spectrum library establishment module is configured to: collect rock samples in the tunnel site area, obtain accurate mineral identification results, establish a rock spectrum library and an end-member mineral spectrum library in the tunnel site area, and obtain the spectrum variation characteristics of different minerals mixed in the rock;

定性识别模块,被配置为:沿隧道掘进方向获取隧道岩石光谱数据,将两个波段的光谱数据进行融合,利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配,确定隧道岩石的大致矿物组合情况;The qualitative identification module is configured to: obtain tunnel rock spectral data along the tunnel excavation direction, fuse the spectral data of the two bands, use the overall waveform of the fused spectral data to perform preliminary matching with the data in the rock spectral library of the tunnel site area, and determine the approximate mineral composition of the tunnel rock;

利用岩石光谱的特征峰谱带位置和特征峰精细特征与组合特征,基于端元矿物光谱库对矿物种类进行精细判识;Using the characteristic peak band position of rock spectra and the fine and combined characteristics of characteristic peaks, the mineral types are precisely identified based on the end-member mineral spectral library.

建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;Establish the association between spectral variation characteristics and mineral types, assist in identifying the types of tunnel rock minerals that cannot be determined through fine identification, and complete the hierarchical qualitative identification of tunnel rock minerals;

定量识别模块,被配置为:对隧道岩石光谱解混,利用端元矿物的纯净光谱逆向反推,选取强相关性波段,利用强相关性波段重新进行解混迭代,直至强相关性波段数量不变,完成对矿物的定量识别。The quantitative identification module is configured to: unmix the tunnel rock spectrum, reversely infer using the pure spectrum of the end-member minerals, select the strongly correlated bands, and re-iterate the unmixing using the strongly correlated bands until the number of strongly correlated bands remains unchanged, thus completing the quantitative identification of the minerals.

本发明第三方面提供了计算机可读存储介质,其上存储有程序,该程序被处理器执行时实现如本发明第一方面所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。The third aspect of the present invention provides a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps of the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in the first aspect of the present invention.

本发明第四方面提供了电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如本发明第一方面所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。The fourth aspect of the present invention provides an electronic device, including a memory, a processor, and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps in the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in the first aspect of the present invention are implemented.

以上一个或多个技术方案存在以下有益效果:One or more of the above technical solutions have the following beneficial effects:

(1)本发明提供了一种协同多元光谱的隧道岩石矿物识别方法及系统,通过联合可见光-近红外波段与热红外波段实现了矿物识别种类的全覆盖,实现了多元光谱协同的矿物定性定量识别,为不良地质识别提供更丰富的矿物异常标志,提高了不良地质识别的准确性,提高了隧道内矿物识别的效率,并通过定量化的数据消除了主观因素对不良地质识别结果的影响。(1) The present invention provides a tunnel rock mineral identification method and system based on collaborative multi-spectral analysis. By combining the visible light-near infrared band and the thermal infrared band, full coverage of mineral identification types is achieved, and multi-spectral collaborative qualitative and quantitative identification of minerals is realized. It provides more abundant mineral anomaly signs for poor geological identification, improves the accuracy of poor geological identification, improves the efficiency of mineral identification in tunnels, and eliminates the influence of subjective factors on the poor geological identification results through quantitative data.

(2)本发明通过充分挖掘围岩光谱特征,将多特征进行深层次联合进行矿物识别,矿物定性的光谱判据层层递进,各判据互相验证,将矿物光谱知识充分融入到识别过程中,并考虑了混合矿物光谱特征的变异性,有效避免光谱的多解性,提高了岩石光谱数据的可解性和准确性。(2) The present invention fully exploits the spectral characteristics of the surrounding rock and combines multiple characteristics in depth for mineral identification. The spectral criteria for mineral qualitative identification are progressive and mutually verified. The mineral spectral knowledge is fully integrated into the identification process. The variability of the spectral characteristics of mixed minerals is taken into account, which effectively avoids the multi-solution of the spectrum and improves the solvability and accuracy of rock spectral data.

(3)本发明提出的定量方法基于逆向思维,推理过程具有扎实的理论基础,通过定量化的误差降低波段选取的主观性,保证模型的客观准确性,并通过解混的强相关波段的不断迭代达到最优的解混效果,实现了隧道内矿物的快速定量化识别以及不良地质识别由定性到定量的跨越。(3) The quantitative method proposed in the present invention is based on reverse thinking, and the reasoning process has a solid theoretical foundation. It reduces the subjectivity of band selection through quantitative errors, ensures the objective accuracy of the model, and achieves the optimal unmixing effect through continuous iteration of unmixed strongly correlated bands, thus realizing the rapid quantitative identification of minerals in the tunnel and the transition from qualitative to quantitative identification of adverse geological conditions.

本发明附加方面的优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Advantages of additional aspects of the present invention will be given in part in the following description, and in part will become obvious from the following description, or will be learned through practice of the present invention.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

图1为第一个实施例的方法流程图。FIG. 1 is a flow chart of a method according to a first embodiment.

图2为钠长石与高岭石混合后热红外波段的变异特征图。Figure 2 is a characteristic diagram of the variation of the thermal infrared band after mixing albite and kaolinite.

图3为地物波谱仪与傅里叶变换红外光谱仪掌子面测点图。Figure 3 is a map of measuring points on the tunnel face using the ground feature spectrometer and the Fourier transform infrared spectrometer.

图4为隧道围岩的光谱实测数据图。Figure 4 is the spectrum measured data of the tunnel surrounding rock.

图5为实测岩石光谱与计算岩石光谱的波段误差图。Figure 5 is a band error diagram of the measured rock spectrum and the calculated rock spectrum.

具体实施方式Detailed ways

应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。It should be noted that the following detailed descriptions are exemplary and are intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。It should be noted that the terms used herein are for describing specific embodiments only and are not intended to be limiting of exemplary embodiments according to the present invention.

在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。In the absence of conflict, the embodiments of the present invention and the features of the embodiments may be combined with each other.

本发明提出的总体思路:The overall idea proposed by the present invention is:

本发明提供了一种协同多元光谱的隧道岩石矿物识别方法及系统,通过前期的地质勘察为围岩光谱的定性解译提供地质约束,充分利用光谱的特征信息,选取“光谱整体波形-特征峰吸收谱带-特征峰精细特征-特征峰组合特征-混合光谱变异特征”的指标,形成层级式矿物定性方法实现矿物种类识别,进而通过选取强相关性波段,并不断进行解混模型的迭代实现矿物的准确定量反演,为不良地质识别提供矿物异常标志。The present invention provides a tunnel rock mineral identification method and system based on collaborative multivariate spectroscopy. The method provides geological constraints for the qualitative interpretation of surrounding rock spectra through preliminary geological surveys, makes full use of the characteristic information of the spectrum, selects indicators such as "overall spectrum waveform-characteristic peak absorption band-characteristic peak fine characteristics-characteristic peak combination characteristics-mixed spectrum variation characteristics", forms a hierarchical mineral qualitative method to realize mineral type identification, and then realizes accurate quantitative inversion of minerals by selecting strongly correlated bands and continuously iterating the unmixing model, thereby providing mineral anomaly marks for poor geological identification.

实施例一Embodiment 1

本实施例公开了协同多元光谱的隧道岩石矿物识别方法。This embodiment discloses a tunnel rock and mineral identification method using collaborative multivariate spectroscopy.

如图1所示,协同多元光谱的隧道岩石矿物识别方法,包括:As shown in FIG1 , the tunnel rock mineral identification method based on collaborative multivariate spectroscopy includes:

(1)进行前期地质勘察,充分调查隧址区的地质特征,为矿物定性识别提供地质信息约束;(1) Conduct preliminary geological surveys to fully investigate the geological characteristics of the tunnel site and provide geological information constraints for qualitative identification of minerals;

(2)采集隧址区岩石样品,采用实验室内精确矿物识别方法,获取准确的矿物组成信息,并通过地物波谱仪和傅里叶变化红外光谱仪进行光谱测试,建立隧址区岩石光谱库,并根据岩石矿物识别结果建立端元矿物光谱库;(2) Collect rock samples from the tunnel site, use precise mineral identification methods in the laboratory to obtain accurate mineral composition information, and perform spectral testing using a ground feature spectrometer and a Fourier transform infrared spectrometer to establish a rock spectral library for the tunnel site. Based on the rock mineral identification results, an end-member mineral spectral library is established;

(3)通过步骤(2)获取的岩石矿物信息、光谱数据以及光谱库,探明岩石中不同矿物混合后的光谱变异特征,辅助判识矿物种类;(3) Using the rock mineral information, spectral data and spectral library obtained in step (2), the spectral variation characteristics of different minerals mixed in the rock are explored to assist in identifying the mineral types;

(4)采用地物波谱仪和傅里叶变换红外光谱仪沿隧道掘进方向进行大量测试,获取隧道岩石的光谱数据;(4) Using a ground object spectrometer and a Fourier transform infrared spectrometer, a large number of tests were carried out along the tunnel excavation direction to obtain spectral data of the tunnel rock;

(5)对获取的岩石光谱数据进行预处理,提升数据的质量,突出光谱特征,提高光谱分析精度;(5) Preprocess the acquired rock spectral data to improve the data quality, highlight the spectral features, and improve the accuracy of spectral analysis;

(6)对获取的数据进行矿物定性识别:a)首先对两个波段的光谱数据进行数据集融合,利用光谱的整体波形,根据前期建立已知岩石矿物成分的岩石光谱库进行匹配,进行矿物共生组合的初步判识;b)再分别根据岩石光谱的其他特征,如“特征峰谱带位置(峰值处的波长)→特征峰的精细特征(特征峰形态、特征峰对称性、特征峰是否为双峰或多峰等特征)→特征峰的组合特征(多个特征峰的组合形式)”进行矿物种类的精细判识;c)对于难以通过上述特征进行判识的矿物,通过步骤(3)中探明的矿物混合变异规律知识(某种矿物发生混合后特征峰谱带发生偏移或特征峰形态发生变化等)进行辅助判识,使变异后非典型的光谱特征(发生偏移的特征峰位置或特征峰的形态等)与矿物种类建立关系,进行辅助判识;(6) Performing qualitative identification of minerals on the acquired data: a) First, the spectral data of the two bands are fused, and the overall waveform of the spectrum is used to match the rock spectrum library with known rock mineral components established in advance to perform preliminary identification of the mineral symbiosis; b) Then, based on other characteristics of the rock spectrum, such as "characteristic peak band position (wavelength at the peak) → fine characteristics of the characteristic peak (characteristic peak shape, characteristic peak symmetry, whether the characteristic peak is double-peaked or multi-peaked, etc.) → combination characteristics of the characteristic peak (combination form of multiple characteristic peaks)", the mineral type is finely identified; c) For minerals that are difficult to identify based on the above characteristics, auxiliary identification is performed based on the knowledge of mineral mixing variation rules discovered in step (3) (the characteristic peak band shifts or the characteristic peak shape changes after a certain mineral is mixed), so that the atypical spectral characteristics after variation (the shifted characteristic peak position or the characteristic peak shape, etc.) are related to the mineral type for auxiliary identification;

(7)通过步骤(6)获得隧道围岩的矿物组成信息后,进行矿物的定量识别:a)首先进行解混以初步获取混合光谱中每种矿物光谱的贡献比例,即矿物的含量;b)利用端元矿物的纯净光谱,按照步骤a)中解混后的含量按照混合光谱解混的规则进行逆向反推,计算获取岩石光谱,称为计算岩石光谱;c)求解实测岩石光谱和计算岩石光谱每个波段的误差值;d)根据误差大小情况,选取满足误差范围的波段作为强相关性波段;e)再利用强相关波段进行解混;f)以上方法不断进行迭代,直至强相关性波段的数量不变结束。(7) After obtaining the mineral composition information of the tunnel surrounding rock through step (6), quantitative identification of minerals is performed: a) first unmixing is performed to preliminarily obtain the contribution ratio of each mineral spectrum in the mixed spectrum, that is, the mineral content; b) using the pure spectrum of the end-member mineral, according to the content after unmixing in step a), reverse inference is performed according to the mixed spectrum unmixing rules to calculate and obtain the rock spectrum, which is called the calculated rock spectrum; c) solving the error value of each band of the measured rock spectrum and the calculated rock spectrum; d) according to the error size, the band that meets the error range is selected as the strong correlation band; e) the strong correlation band is used again for unmixing; f) the above method is iterated continuously until the number of strong correlation bands remains unchanged.

更为具体的,包括:More specifically, they include:

(1)对某隧道沿线进行前期地质勘察,调查隧道地表以及洞内的地质信息,为矿物定性识别提供地质信息约束;(1) Conduct preliminary geological survey along a tunnel to investigate the geological information on the tunnel surface and inside the tunnel to provide geological information constraints for qualitative identification of minerals;

所述地质信息主要包括隧道沿线地表穿越地层的岩性、历史发生的地质作用类型、不良地质的位置等信息;The geological information mainly includes the lithology of the ground strata along the tunnel, the types of geological actions that occurred in the past, the location of unfavorable geology, etc.

所述进行前期地质勘察的方法包括纵向钻进、地表采样、航空电磁物探等方法;The methods for conducting preliminary geological surveys include longitudinal drilling, surface sampling, airborne electromagnetic geophysical exploration, etc.;

其中,对隧道沿线进行勘察,探明该隧道沿线地层为花岗岩地层,历史未发生重大的地质构造运动,但在隧道洞口前方200m左右位置有一规模较小断层,有粘土矿物产生。Among them, an investigation was conducted along the tunnel, and it was found that the strata along the tunnel were granite strata, and no major geological tectonic movements had occurred in history. However, there was a small-scale fault about 200m in front of the tunnel entrance, where clay minerals were produced.

(2)采集隧址区岩石样品,采用实验室内精确矿物识别方法,获取准确的矿物组成信息,并通过地物波谱仪和傅里叶变化红外光谱仪进行光谱测试,建立隧址区岩石光谱库,并根据岩石矿物识别结果建立端元矿物光谱库;(2) Collect rock samples from the tunnel site, use precise mineral identification methods in the laboratory to obtain accurate mineral composition information, and perform spectral testing using a ground feature spectrometer and a Fourier transform infrared spectrometer to establish a rock spectral library for the tunnel site. Based on the rock mineral identification results, an end-member mineral spectral library is established;

采集隧址区演示样品,即将步骤(1)中地质勘查所采集的岩芯、岩块等具有证明性的样品收纳编号;Collect demonstration samples from the tunnel site, i.e., collect and number the cores, rock blocks and other demonstrative samples collected in the geological survey in step (1);

采用实验室内精确矿物识别方法,包括但不限于镜下薄片法、X射线衍射法等其他方法;Use accurate laboratory mineral identification methods, including but not limited to thin section microscopy, X-ray diffraction and other methods;

通过地物波谱仪及傅里叶变换红外光谱仪测试,测试时需选择样品平整面,防止漏光等一系列降低数据质量的可能;When testing with a ground object spectrometer and a Fourier transform infrared spectrometer, the sample surface must be flat to prevent light leakage and other possible factors that may reduce data quality.

建立端元光谱库,即将获取精确矿物识别结果后,将这些矿物的可见光-近红外及热红外波段的光谱收集,建立隧道现场的端元矿物光谱库;Establish an end-member spectral library, that is, after obtaining accurate mineral identification results, collect the spectra of these minerals in the visible light-near infrared and thermal infrared bands to establish an end-member mineral spectral library at the tunnel site;

建立端元矿物光谱库的方法包括从现有的USGS波谱库、JHU波谱库、ASU波谱库等光谱库中搜集,或通过购置纯净矿物,然后用光谱仪采集,建立端元矿物光谱库。Methods for establishing an end-member mineral spectral library include collecting from existing spectral libraries such as the USGS spectral library, the JHU spectral library, and the ASU spectral library, or purchasing pure minerals and then collecting them with a spectrometer to establish an end-member mineral spectral library.

本实施例中,将地质勘查所采集的样品按里程、坐标等信息进行编号,采用X射线衍射的方法对矿物进行识别,结果表明大部分岩石内部含有钾长石、钠长石、白云母、高岭石、伊利石等矿物,并通过两种光谱技术对上述纯净矿物进行光谱测试,建立端元矿物光谱库。In this embodiment, the samples collected by geological survey are numbered according to mileage, coordinates and other information, and the minerals are identified by X-ray diffraction. The results show that most rocks contain minerals such as potassium feldspar, albite, muscovite, kaolinite, illite, etc., and the above pure minerals are spectrally tested using two spectral techniques to establish an end-member mineral spectral library.

(3)通过步骤(2)获取的岩石矿物信息、光谱数据以及光谱库,探明岩石中不同矿物混合后的光谱变异特征,辅助判识矿物种类;(3) Using the rock mineral information, spectral data and spectral library obtained in step (2), the spectral variation characteristics of different minerals mixed in the rock are explored to assist in identifying the mineral types;

矿物混合变异特征,指多种矿物混合后,它们的混合光谱特征会出现偏移、湮没等一系列与纯净矿物光谱特征不符的现象;Mineral mixture variation characteristics refer to the phenomenon that after multiple minerals are mixed, their mixed spectral characteristics will be offset, annihilated, and a series of other phenomena that are inconsistent with the spectral characteristics of pure minerals;

探明不同矿物混合的光谱变异特征对于将矿物混合后的非典型特征与矿物种类建立联系具有积极作用。Exploring the spectral variation characteristics of different mineral mixtures plays a positive role in establishing a connection between the atypical characteristics of mineral mixtures and mineral types.

本实施例中,以钠长石与高岭石混合为例,在近红外波段无明显变异现象,在热红外波段,钠长石的特征峰原位于9616nm,混合后向短波方向发生偏移,其余部分特征被湮没,高岭石可从9899nm以及10950nm/10972nm处的特征证实,如图2所示。In this embodiment, taking the mixture of albite and kaolinite as an example, there is no obvious variation phenomenon in the near-infrared band. In the thermal infrared band, the characteristic peak of albite is originally located at 9616nm, and it shifts to the short-wave direction after mixing, and the remaining characteristics are obliterated. Kaolinite can be confirmed from the characteristics at 9899nm and 10950nm/10972nm, as shown in Figure 2.

(4)采用地物波谱仪和傅里叶变换红外光谱仪沿隧道掘进方向进行大量测试,获取隧道岩石的光谱数据;(4) Using a ground object spectrometer and a Fourier transform infrared spectrometer, a large number of tests were carried out along the tunnel excavation direction to obtain spectral data of the tunnel rock;

采用地物波谱仪和傅里叶变换红外光谱仪沿隧道掘进方向测试,测试时仪器探头需与岩石或围岩的平整面完全接触,并压紧,防止漏光等影响数据精度的现象发生。A ground feature spectrometer and a Fourier transform infrared spectrometer are used to test along the tunnel excavation direction. During the test, the instrument probe needs to be in full contact with the flat surface of the rock or surrounding rock and pressed tightly to prevent light leakage and other phenomena that affect data accuracy.

本实施例中,选择人工或将仪器搭载到机器人上对隧道沿线洞壁及掌子面进行全面测试,测试选点应选择具有代表性点位,测试范围覆盖掌子面,如图3所示。In this embodiment, a comprehensive test is performed on the tunnel wall and the tunnel face along the tunnel by manual labor or by placing the instrument on a robot. The test points should be representative points, and the test range covers the tunnel face, as shown in FIG3 .

(5)对获取的岩石光谱数据进行预处理,提升数据的质量,突出光谱特征,提高光谱分析精度;(5) Preprocess the acquired rock spectral data to improve the data quality, highlight the spectral features, and improve the accuracy of spectral analysis;

对获取的岩石光谱数据进行预处理,其中预处理方法包括但不仅限于多元散射校正、微分、平滑、包络线去除等方法;Preprocessing the acquired rock spectrum data, wherein the preprocessing methods include but are not limited to multivariate scattering correction, differentiation, smoothing, envelope removal and other methods;

本实施例中,步骤(4)获取围岩的全波段光谱数据后,首先进行SG平滑,去除噪声,然后通过多元散校正排除体积散射的影响,同时对近红外数据进行包络线去除,对热红外数据做基线校正,突出其吸收反射特征。In this embodiment, after step (4) obtains the full-band spectral data of the surrounding rock, SG smoothing is first performed to remove noise, and then multivariate dispersion correction is used to eliminate the influence of volume scattering. At the same time, the envelope of the near-infrared data is removed, and the baseline correction is performed on the thermal infrared data to highlight its absorption and reflection characteristics.

(6)对获取的数据进行矿物定性识别:(6) Qualitative identification of minerals on the acquired data:

a)首先对两个波段的光谱数据进行数据集融合,利用光谱的整体波形,根据前期建立已知岩石矿物成分的岩石光谱库进行匹配,进行矿物共生组合的初步判识;a) First, the spectral data of the two bands are fused, and the overall waveform of the spectrum is used to match the rock spectrum library with known rock mineral components established in the early stage to make a preliminary identification of the mineral paragenesis combination;

b)再分别根据岩石光谱的其他特征,如“特征峰谱带位置(峰值处的波长)→特征峰的精细特征(特征峰形态、特征峰对称性、特征峰是否为双峰或多峰等特征)→特征峰的组合特征(多个特征峰的组合形式)”进行矿物种类的精细判识;b) Then, according to other characteristics of the rock spectrum, such as "characteristic peak band position (wavelength at the peak) → fine characteristics of characteristic peaks (characteristic peak shape, characteristic peak symmetry, whether the characteristic peak is double-peaked or multi-peaked, etc.) → combination characteristics of characteristic peaks (combination form of multiple characteristic peaks)", the mineral type is finely identified;

c)对于难以通过上述特征进行判识的矿物,通过步骤(3)中探明的矿物混合变异规律知识(某种矿物发生混合后特征峰谱带发生偏移或特征峰形态发生变化等)进行辅助判识,使变异后非典型的光谱特征(发生偏移的特征峰位置或特征峰的形态等)与矿物种类建立关系,进行辅助判识。c) For minerals that are difficult to identify using the above characteristics, auxiliary identification is performed using the knowledge of the mineral mixing variation rules explored in step (3) (such as the shift of the characteristic peak spectral band or the change of the characteristic peak morphology after a certain mineral is mixed), so that a relationship is established between the atypical spectral characteristics after variation (such as the shifted characteristic peak position or the characteristic peak morphology) and the mineral type for auxiliary identification.

步骤a)中,对两个波段的光谱数据进行融合,融合的方法包括但不仅限于累加融合(CF)、等权融合(ERF)、外积融合(OPF)等。In step a), the spectral data of the two bands are fused, and the fusion methods include but are not limited to cumulative fusion (CF), equal weight fusion (ERF), outer product fusion (OPF) and the like.

本实施例中,我们:In this example, we:

a)采用等权融合的方法将预处理后的两个波段的岩石光谱数据进行数据集融合,相对于其他方法,在减少数据维度的同时不遗漏吸收反射特征,利用融合后的数据与建立的隧址区岩石光谱库进行初步匹配,确定大致矿物组合情况,初步判断为钾长石、高岭石与蒙脱石混合,光谱数据如图4;a) The rock spectral data of the two bands after preprocessing were fused by equal weight fusion method. Compared with other methods, the absorption and reflection features were not missed while reducing the data dimension. The fused data were preliminarily matched with the established rock spectral library of the tunnel site to determine the approximate mineral combination. It was preliminarily judged to be a mixture of potassium feldspar, kaolinite and montmorillonite. The spectral data is shown in Figure 4;

b)首先通过特征峰吸收位置与端元光谱库进行综合判识,通过9480nm反射峰可以判断含有钾长石,10972nm反射峰无法识别高岭石与蒙脱石,三种矿物光谱在热红外波段反射特征位置如表1;b) First, comprehensive identification is performed through the characteristic peak absorption position and the end member spectrum library. The 9480nm reflection peak can be used to determine the presence of potassium feldspar, while the 10972nm reflection peak cannot identify kaolinite and montmorillonite. The reflection characteristic positions of the three mineral spectra in the thermal infrared band are shown in Table 1;

表1:三种矿物光谱在热红外波段反射特征位置Table 1: Reflection characteristic positions of the three mineral spectra in the thermal infrared band

通过近红外光谱中2200nm附近的非对称双峰,可以判断含有高岭石,并且在1400与1900nm处没有呈现斜坡状吸收特征以及没有2210nm处的吸收特征,结合热红外波段没有蒙脱石的特征反射峰位置,判断岩石中不存在蒙脱石;The asymmetric double peaks near 2200nm in the near-infrared spectrum indicate the presence of kaolinite. The absence of slope-shaped absorption features at 1400 and 1900nm and absorption features at 2210nm, combined with the absence of the characteristic reflection peak position of montmorillonite in the thermal infrared band, indicates that montmorillonite does not exist in the rock.

c)通过变异规律,9881nm处的反射峰确认为高岭石,验证步骤a)、b)中的判识结果。c) Through the variation law, the reflection peak at 9881nm is confirmed to be kaolinite, verifying the identification results in steps a) and b).

(7)通过步骤(6)获得隧道围岩的矿物组成信息后,进行矿物的定量识别:(7) After obtaining the mineral composition information of the tunnel surrounding rock through step (6), quantitative identification of minerals is performed:

a)首先进行解混以初步获取混合光谱中每种矿物光谱的贡献比例,即矿物的含量;a) firstly perform unmixing to preliminarily obtain the contribution ratio of each mineral spectrum in the mixed spectrum, that is, the content of the mineral;

b)利用端元矿物的纯净光谱,按照步骤a)中解混后的含量按照混合光谱解混的规则进行逆向反推,计算获取岩石光谱,称为计算岩石光谱;b) using the pure spectrum of the end-member minerals, according to the content after unmixing in step a), reverse inference is performed according to the rules of mixed spectrum unmixing, and the rock spectrum is calculated, which is called calculated rock spectrum;

c)求解实测岩石光谱和计算岩石光谱每个波段的误差值;c) Solve the measured rock spectrum and calculate the error value of each band of the rock spectrum;

d)根据误差大小情况,选取满足误差范围的波段作为强相关性波段;d) According to the error size, select the band that meets the error range as the strong correlation band;

e)再利用强相关波段进行解混;e) Reusing strongly correlated bands for unmixing;

f)以上方法不断进行迭代,直至强相关性波段的数量不变结束。f) The above method is iterated continuously until the number of strongly correlated bands remains unchanged.

所述步骤a)及e)的解混方法,包含线性解混与非线性解混方法。The unmixing methods of steps a) and e) include linear unmixing and nonlinear unmixing methods.

本实施例中:In this embodiment:

我们通过(6)确定岩石中含有钾长石和高岭石两种矿物;We determined through (6) that the rock contains two minerals: potassium feldspar and kaolinite;

遵循(7)中a)步骤从端元矿物光谱库中选取钾长石和高岭石两种矿物的光谱,对岩石光谱进行线性解混,获取两种矿物的相对含量;Following step a) in (7), spectra of two minerals, potassium feldspar and kaolinite, are selected from the end-member mineral spectrum library, and the rock spectra are linearly unmixed to obtain the relative contents of the two minerals;

b)将两种端元矿物的光谱与其相对含量按照下列公式进行线性叠加,获取岩石的计算光谱其中W代表岩石光谱,an代表第n种端元矿物的含量,fn代表第n种端元矿物的纯净光谱。b) The spectra of the two end-member minerals and their relative contents are linearly superimposed according to the following formula to obtain the calculated spectrum of the rock: where W represents the rock spectrum, a n represents the content of the nth end-member mineral, and f n represents the pure spectrum of the nth end-member mineral.

W=a1f1+a2f2+......+anfn W= a1f1 + a2f2 + ......+a n f n

(c)获取岩石的计算光谱后,计算其与岩石实测光谱每个波段的误差值(图5);(c) After obtaining the calculated spectrum of the rock, the error value of each band between it and the measured spectrum of the rock is calculated (Figure 5);

(d)综合考虑每个波段的误差值,我们发现在8um~12um之间误差较高,因此不适宜作为强相关性波段,初步选取5um~8um作为强相关性波段;(d) Considering the error value of each band, we found that the error between 8um and 12um is relatively high, so it is not suitable as a strong correlation band. We initially selected 5um to 8um as a strong correlation band;

e)将5um~8um波段的光谱数据进行线性解混,再次获取每种端元矿物含量;e) Linearly unmixing the spectral data in the 5um to 8um band to obtain the content of each end-member mineral again;

(f)以上流程不断迭代,直至选取的相关性波段数量不变为止。(f) The above process is iterated continuously until the number of selected correlation bands remains unchanged.

实施例二Embodiment 2

本实施例公开了协同多元光谱的隧道岩石矿物识别系统。This embodiment discloses a tunnel rock and mineral identification system based on collaborative multivariate spectroscopy.

协同多元光谱的隧道岩石矿物识别系统,包括:The tunnel rock and mineral identification system based on collaborative multi-spectral analysis includes:

光谱库建立模块,被配置为:采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库,并获取岩石中不同矿物混合后的光谱变异特征;The spectrum library establishment module is configured to: collect rock samples in the tunnel site area, obtain accurate mineral identification results, establish a rock spectrum library and an end-member mineral spectrum library in the tunnel site area, and obtain the spectrum variation characteristics of different minerals mixed in the rock;

定性识别模块,被配置为:沿隧道掘进方向获取隧道岩石光谱数据,将两个波段的光谱数据进行融合,利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配,确定隧道岩石的大致矿物组合情况;The qualitative identification module is configured to: obtain tunnel rock spectral data along the tunnel excavation direction, fuse the spectral data of the two bands, use the overall waveform of the fused spectral data to perform preliminary matching with the data in the rock spectral library of the tunnel site area, and determine the approximate mineral composition of the tunnel rock;

利用岩石光谱的特征峰谱带位置和特征峰精细特征与组合特征,基于端元矿物光谱库对矿物种类进行精细判识;Using the characteristic peak band position of rock spectra and the fine and combined characteristics of characteristic peaks, the mineral types are precisely identified based on the end-member mineral spectral library.

建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;Establish the association between spectral variation characteristics and mineral types, assist in identifying the types of tunnel rock minerals that cannot be determined through fine identification, and complete the hierarchical qualitative identification of tunnel rock minerals;

定量识别模块,被配置为:对隧道岩石光谱解混,利用端元矿物的纯净光谱逆向反推,选取强相关性波段,利用强相关性波段重新进行解混迭代,直至强相关性波段数量不变,完成对矿物的定量识别。The quantitative identification module is configured to: unmix the tunnel rock spectrum, reversely infer using the pure spectrum of the end-member minerals, select the strongly correlated bands, and re-iterate the unmixing using the strongly correlated bands until the number of strongly correlated bands remains unchanged, thus completing the quantitative identification of the minerals.

实施例三Embodiment 3

本实施例的目的是提供计算机可读存储介质。The purpose of this embodiment is to provide a computer-readable storage medium.

计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本公开实施例1所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。A computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in Example 1 of the present disclosure.

实施例四Embodiment 4

本实施例的目的是提供电子设备。The purpose of this embodiment is to provide an electronic device.

电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如本公开实施例1所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。An electronic device includes a memory, a processor, and a program stored in the memory and executable on the processor. When the processor executes the program, the steps in the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in Example 1 of the present disclosure are implemented.

以上实施例二、三和四的装置中涉及的各步骤与方法实施例一相对应,具体实施方式可参见实施例一的相关说明部分。术语“计算机可读存储介质”应该理解为包括一个或多个指令集的单个介质或多个介质;还应当被理解为包括任何介质,所述任何介质能够存储、编码或承载用于由处理器执行的指令集并使处理器执行本发明中的任一方法。The steps involved in the apparatuses of the above embodiments 2, 3 and 4 correspond to the method embodiment 1, and the specific implementation methods can refer to the relevant description part of embodiment 1. The term "computer-readable storage medium" should be understood as a single medium or multiple media including one or more instruction sets; it should also be understood to include any medium that can store, encode or carry an instruction set for execution by a processor and enable the processor to execute any method of the present invention.

本领域技术人员应该明白,上述本发明的各模块或各步骤可以用通用的计算机装置来实现,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。本发明不限制于任何特定的硬件和软件的结合。Those skilled in the art should understand that the modules or steps of the present invention described above can be implemented by a general-purpose computer device, or alternatively, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. The present invention is not limited to any specific combination of hardware and software.

上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。Although the above describes the specific implementation mode of the present invention in conjunction with the accompanying drawings, it is not intended to limit the scope of protection of the present invention. Technical personnel in the relevant field should understand that various modifications or variations that can be made by technical personnel in the field without creative work on the basis of the technical solution of the present invention are still within the scope of protection of the present invention.

Claims (10)

1.协同多元光谱的隧道岩石矿物识别方法,其特征在于,包括以下步骤:1. A tunnel rock mineral identification method based on collaborative multivariate spectroscopy, characterized in that it comprises the following steps: 采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库,并获取岩石中不同矿物混合后的光谱变异特征;Collect rock samples in the tunnel site area, obtain accurate mineral identification results, establish a rock spectral library and an end-member mineral spectral library in the tunnel site area, and obtain the spectral variation characteristics of different minerals mixed in the rock; 沿隧道掘进方向获取隧道岩石光谱数据,将两个波段的光谱数据进行融合,利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配,确定隧道岩石的大致矿物组合情况;Acquire tunnel rock spectral data along the tunnel excavation direction, fuse the spectral data of the two bands, and use the overall waveform of the fused spectral data to make a preliminary match with the data in the rock spectral library of the tunnel site to determine the approximate mineral composition of the tunnel rock; 所述将两个波段的光谱数据进行融合是通过联合可见光-近红外波段与热红外波段实现对矿物识别种类的全覆盖;The fusion of the spectral data of the two bands is to achieve full coverage of mineral identification types by combining the visible light-near infrared band and the thermal infrared band; 利用岩石光谱的特征峰谱带位置和特征峰精细特征与组合特征,基于端元矿物光谱库对矿物种类进行精细判识;Using the characteristic peak band position of rock spectra and the fine and combined characteristics of characteristic peaks, the mineral types are precisely identified based on the end-member mineral spectral library. 建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;Establish the association between spectral variation characteristics and mineral types, assist in identifying the types of tunnel rock minerals that cannot be determined through fine identification, and complete the hierarchical qualitative identification of tunnel rock minerals; 对隧道岩石光谱解混,利用端元矿物的纯净光谱逆向反推,选取强相关性波段,利用强相关性波段重新进行解混迭代,直至强相关性波段数量不变,完成对矿物的定量识别。The tunnel rock spectrum is unmixed, and the pure spectrum of the end-member mineral is used for reverse inference. The strongly correlated bands are selected and the unmixing is iterated again using the strongly correlated bands until the number of strongly correlated bands remains unchanged, thus completing the quantitative identification of the minerals. 2.如权利要求1所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于:2. The tunnel rock and mineral identification method based on collaborative multivariate spectroscopy as claimed in claim 1, characterized in that: 对获取的隧址区岩石样品采用实验室内精确矿物识别方法进行识别,获取精确矿物识别结果;The rock samples obtained from the tunnel site are identified using precise mineral identification methods in the laboratory to obtain precise mineral identification results; 利用地物波谱仪和傅里叶变化红外光谱仪对隧址区岩石样品进行光谱测试,建立隧址区岩石光谱库;Use ground object spectrometer and Fourier transform infrared spectrometer to carry out spectrum test on rock samples in tunnel site, and establish rock spectrum library in tunnel site; 从现有的光谱库中搜集精确矿物识别结果中纯净矿物的光谱数据,或通过购置纯净矿物、用光谱仪采集光谱数据,建立端元矿物光谱库。Collect the spectral data of pure minerals from the accurate mineral identification results from the existing spectral library, or establish an end-member mineral spectral library by purchasing pure minerals and collecting spectral data with a spectrometer. 3.如权利要求1所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于,沿隧道掘进方向获取隧道岩石光谱数据,具体为:3. The tunnel rock mineral identification method based on collaborative multivariate spectroscopy as claimed in claim 1 is characterized in that the tunnel rock spectral data is obtained along the tunnel excavation direction, specifically: 对隧道沿线洞壁及掌子面进行全面测试,选择具有代表性点位作为测试选点,测试范围覆盖掌子面,获取隧道岩石光谱数据。Comprehensive tests are conducted on the tunnel walls and faces along the tunnel, representative points are selected as test points, the test range covers the face, and tunnel rock spectral data is obtained. 4.如权利要求1所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于,所述特征峰精细特征包括特征峰形态、特征峰对称性、特征峰是否为双峰或多峰,所述特征峰组合特征为多个特征峰的组合形式。4. The collaborative multivariate spectroscopy tunnel rock mineral identification method as described in claim 1 is characterized in that the characteristic peak fine features include characteristic peak morphology, characteristic peak symmetry, whether the characteristic peak is bimodal or multimodal, and the characteristic peak combination features are the combination form of multiple characteristic peaks. 5.如权利要求1所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于,对矿物进行定量识别,具体为:5. The tunnel rock mineral identification method based on collaborative multivariate spectroscopy as claimed in claim 1 is characterized in that the quantitative identification of minerals is specifically: a)对岩石光谱进行解混,初步获取矿物含量;a) Unmix the rock spectra and obtain the preliminary mineral content; b)利用端元矿物光谱库获取端元矿物的纯净光谱,基于获取的矿物含量,按照混合光谱解混的规则进行逆向反推,获取计算岩石光谱;b) Use the end-member mineral spectral library to obtain the pure spectrum of the end-member minerals, and based on the obtained mineral content, perform reverse inference according to the rules of mixed spectrum unmixing to obtain the calculated rock spectrum; c)求解实测岩石光谱和计算岩石光谱每个波段的误差值;c) Solve the measured rock spectrum and calculate the error value of each band of the rock spectrum; d)根据误差大小情况,选取满足误差范围的波段作为强相关性波段;d) According to the error size, select the band that meets the error range as the strong correlation band; e)利用强相关波段进行解混,获取新的矿物含量;e) Use strongly correlated bands for unmixing to obtain new mineral content; f)迭代上述步骤b)- e),直至强相关性波段的数量不变结束。f) Iterate the above steps b) - e) until the number of strongly correlated bands remains unchanged. 6.如权利要求5所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于,所述解混方法包括线性解混方法与非线性解混方法。6. The tunnel rock mineral identification method based on collaborative multivariate spectroscopy as described in claim 5 is characterized in that the unmixing method includes a linear unmixing method and a nonlinear unmixing method. 7.如权利要求1所述的协同多元光谱的隧道岩石矿物识别方法,其特征在于,在采集隧址区岩石样品之前,还包括:7. The tunnel rock and mineral identification method based on collaborative multivariate spectroscopy as claimed in claim 1, characterized in that before collecting rock samples in the tunnel site area, it also includes: 对隧道沿线进行前期地质勘察,调查隧道地表以及洞内的地质信息,为矿物定性识别提供地质信息约束。Conduct preliminary geological surveys along the tunnel to investigate the geological information on the tunnel surface and inside the tunnel, and provide geological information constraints for qualitative identification of minerals. 8.协同多元光谱的隧道岩石矿物识别系统,其特征在于:执行时实现如权利要求1-7任一项所述的协同多元光谱的隧道岩石矿物识别方法,包括:8. A collaborative multi-spectral tunnel rock and mineral identification system, characterized in that: when executed, the collaborative multi-spectral tunnel rock and mineral identification method according to any one of claims 1 to 7 is implemented, comprising: 光谱库建立模块,被配置为:采集隧址区岩石样品,获取精确矿物识别结果,建立隧址区岩石光谱库和端元矿物光谱库,并获取岩石中不同矿物混合后的光谱变异特征;The spectrum library establishment module is configured to: collect rock samples in the tunnel site area, obtain accurate mineral identification results, establish a rock spectrum library and an end-member mineral spectrum library in the tunnel site area, and obtain the spectrum variation characteristics of different minerals mixed in the rock; 定性识别模块,被配置为:沿隧道掘进方向获取隧道岩石光谱数据,将两个波段的光谱数据进行融合,利用融合后光谱数据的整体波形与隧址区岩石光谱库中的数据进行初步匹配,确定隧道岩石的大致矿物组合情况;The qualitative identification module is configured to: obtain tunnel rock spectral data along the tunnel excavation direction, fuse the spectral data of the two bands, use the overall waveform of the fused spectral data to perform a preliminary match with the data in the rock spectral library of the tunnel site area, and determine the approximate mineral composition of the tunnel rock; 利用岩石光谱的特征峰谱带位置和特征峰精细特征与组合特征,基于端元矿物光谱库对矿物种类进行精细判识;Using the characteristic peak band positions, fine features and combination features of the rock spectra, the mineral types are precisely identified based on the end-member mineral spectral library. 建立光谱变异特征与矿物种类的关联,对无法通过精细判识确定的隧道岩石矿物种类进行辅助判识,完成隧道岩石矿物的层级式定性;Establish the association between spectral variation characteristics and mineral types, assist in identifying the types of tunnel rock minerals that cannot be determined through fine identification, and complete the hierarchical qualitative identification of tunnel rock minerals; 定量识别模块,被配置为:对隧道岩石光谱解混,利用端元矿物的纯净光谱逆向反推,选取强相关性波段,利用强相关性波段重新进行解混迭代,直至强相关性波段数量不变,完成对矿物的定量识别。The quantitative identification module is configured to: unmix the tunnel rock spectrum, reversely infer using the pure spectrum of the end-member minerals, select the strongly correlated bands, and re-iterate the unmixing using the strongly correlated bands until the number of strongly correlated bands remains unchanged, thus completing the quantitative identification of the minerals. 9.计算机可读存储介质,其上存储有程序,其特征在于,该程序被处理器执行时实现如权利要求1-7任一项所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。9. A computer-readable storage medium having a program stored thereon, wherein when the program is executed by a processor, the program implements the steps of the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in any one of claims 1 to 7. 10.电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-7任一项所述的协同多元光谱的隧道岩石矿物识别方法中的步骤。10. An electronic device comprising a memory, a processor, and a program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the collaborative multivariate spectroscopy tunnel rock mineral identification method as described in any one of claims 1 to 7 are implemented.
CN202311816245.7A 2023-12-26 2023-12-26 Tunnel rock mineral identification method and system with cooperative multi-element spectrum Active CN117828397B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311816245.7A CN117828397B (en) 2023-12-26 2023-12-26 Tunnel rock mineral identification method and system with cooperative multi-element spectrum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311816245.7A CN117828397B (en) 2023-12-26 2023-12-26 Tunnel rock mineral identification method and system with cooperative multi-element spectrum

Publications (2)

Publication Number Publication Date
CN117828397A CN117828397A (en) 2024-04-05
CN117828397B true CN117828397B (en) 2024-07-12

Family

ID=90520498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311816245.7A Active CN117828397B (en) 2023-12-26 2023-12-26 Tunnel rock mineral identification method and system with cooperative multi-element spectrum

Country Status (1)

Country Link
CN (1) CN117828397B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118730927B (en) * 2024-09-02 2024-12-27 洛阳市灿业建设工程有限公司 Highway tunnel detecting system for highway engineering
CN118940844B (en) * 2024-10-12 2025-05-09 山东大学 Tunnel surrounding rock spectrum interpretation method and system based on large language model
CN119164901B (en) * 2024-11-21 2025-11-07 山东大学 Element and spectrum characteristic-based coal seam exploitation in-situ identification method and system
CN119804377B (en) * 2025-03-11 2025-05-16 永春县农业科学研究所(永春县农业检测中心、永春县作物良种场) Agricultural product quality detection method and system
CN121068531B (en) * 2025-11-05 2026-02-10 中国矿业大学 Geological sample component identification and anomaly detection system based on machine learning

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551471A (en) * 2009-05-19 2009-10-07 中国国土资源航空物探遥感中心 High-spectrum remote-sensing mineral content quantitative inversion method
CN101887012A (en) * 2010-06-28 2010-11-17 中国国土资源航空物探遥感中心 Spectral reflectance peak decomposition based quantitative inversion method of hyperspectral remote sensing mineral content

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009017487B4 (en) * 2009-04-15 2013-04-18 Bundesrepublik Deutschland, vertreten durch den Präsidenten der Bundesanstalt für Geowissenschaften und Rohstoffe Method and device for determining a mineral content of a clay mineral, in particular swellable clay mineral in a rock
US11821857B2 (en) * 2020-10-06 2023-11-21 Baker Hughes Oilfield Operations Llc Data-driven solutions for inverse elemental modeling
CN115753632B (en) * 2022-10-19 2024-05-31 山东大学 Real-time identification method and system of unfavorable geological bodies in tunnels based on image spectrum
CN115656053B (en) * 2022-10-19 2024-05-31 山东大学 Rock mineral content testing method and system
CN115641453A (en) * 2022-10-19 2023-01-24 山东大学 A method and system for identifying weathering degree of tunnel surrounding rock based on hyperspectral technology
CN115761038B (en) * 2022-10-19 2023-06-30 山东大学 A Method and System for Geological Sketching of Tunnel Face Based on Image Spectrum Technology

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551471A (en) * 2009-05-19 2009-10-07 中国国土资源航空物探遥感中心 High-spectrum remote-sensing mineral content quantitative inversion method
CN101887012A (en) * 2010-06-28 2010-11-17 中国国土资源航空物探遥感中心 Spectral reflectance peak decomposition based quantitative inversion method of hyperspectral remote sensing mineral content

Also Published As

Publication number Publication date
CN117828397A (en) 2024-04-05

Similar Documents

Publication Publication Date Title
CN117828397B (en) Tunnel rock mineral identification method and system with cooperative multi-element spectrum
CN110967311B (en) Porphyry deposit alteration zone identification method based on infrared spectrum and magnetic susceptibility measurement
CN110333200B (en) Method for delineating mineralization center based on short-wave infrared spectrum
CN105257286B (en) A method and device for obtaining the content of formation rock components
CN103513270B (en) A kind of gas-bearing formation identification and evaluation method and device based on rock acoustic property
CN101975767B (en) Method for measuring bonding potential between soluble organisms and metal ions in water environment
EP3809133B1 (en) A method for characterizing underground metallic mineral deposits based on rock coatings and fracture fills
CN101368481A (en) Quantitative Analysis Method of NMR Rock Oil Abundance in Oil Drilling
CN107589094A (en) Anshan type iron mine stone type based on spectral signature determines method
Hassan et al. The effect of mineralogical composition on strength and drillability of granitic rocks in Hulu Langat, Selangor Malaysia
CN115046941A (en) Hydrothermal uranium deposit alteration zonal identification method based on multi-scale spectral measurement
CN109738955A (en) A kind of metamorphic rock lithology comprehensive distinguishing method based under ingredient-textural classification
Liu et al. A New Quantitative Approach for Element‐Mineral Determination Based on “EDS (Energy Dispersive Spectroscopy) Method”
CN119104580A (en) A high-precision recovery method and system for judging the sedimentary environment of marine mixed shale
CN109001422B (en) A kind of geochemical discrimination method of granite choritization alteration extent
CN117348071A (en) A lithology identification method for diamictite formations based on fusion of multiple well logging information
CN102330553B (en) A Dynamic Spectral Fluid Identification Method for MDT Testing
Huang et al. Conversion models between elements and mineral compositions and their applications for rapid determination of the brittleness of shale
CN106525818A (en) LIBS Mars substance analysis method based on multi-database support and multi-link correction
CN114441581B (en) Discrimination analysis method for multi-stage cause dense dolomite
CN110135298A (en) A Method for Cataloging Drilling Core Alteration Information
Zhenhao et al. Four typical variation patterns of mid-infrared spectra of the mafic mineral anomalies for fault zone identification
CN118817637B (en) Method for discriminating temperature and pressure conditions of pegmatite based on muscovite spectral characteristics
CN108240999A (en) An identification method of clastic rock lithology
CN111487393A (en) Method for measuring paleo-salinity value, single-well and multi-well paleo-salinity plan and method for drawing sedimentary facies development plan

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant