CN114266179B - 基于有限单元的矿床钻探信息处理及分析方法和装置 - Google Patents

基于有限单元的矿床钻探信息处理及分析方法和装置 Download PDF

Info

Publication number
CN114266179B
CN114266179B CN202111586366.8A CN202111586366A CN114266179B CN 114266179 B CN114266179 B CN 114266179B CN 202111586366 A CN202111586366 A CN 202111586366A CN 114266179 B CN114266179 B CN 114266179B
Authority
CN
China
Prior art keywords
finite element
information
drilling
section
analysis
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
CN202111586366.8A
Other languages
English (en)
Other versions
CN114266179A (zh
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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202111586366.8A priority Critical patent/CN114266179B/zh
Publication of CN114266179A publication Critical patent/CN114266179A/zh
Priority to US18/066,605 priority patent/US20230204810A1/en
Application granted granted Critical
Publication of CN114266179B publication Critical patent/CN114266179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G01V1/50—Analysing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/46—Data acquisition
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了基于有限单元的矿床钻探信息处理及分析方法和装置，其中，所述信息处理方法包括：基于勘查钻孔的相关信息进行有限单元划分，并构建其对应的分段有限单元平面图；进一步的，将分段有限单元平面图进行三维空间同位堆叠，并通过空间分析，构建扫描解析有限单元剖面图。本发明可实现对地质钻孔信息分段分层地可视化解析呈现，对精准圈定找矿靶区、精准采矿选冶、构造成矿学研究等方面具有很好的实用价值和指示意义。

Description

基于有限单元的矿床钻探信息处理及分析方法和装置
技术领域
本发明涉及矿床钻探信息分析方法的技术领域。
背景技术
如何解译和呈现矿床勘查信息，并藉此指导矿产资源勘探开发与综合利用，是矿床学研究的重要研究内容。相较于在地表和浅近地下开展的基础地质调查、地球化学勘查、地球物理勘查、遥感、探矿工程等，勘探钻孔信息能够直接、客观地反映钻孔范围内地下岩性、产状、结构、构造、矿化、蚀变、隐爆破碎等地质信息，为地质勘查工作提供直接、真实、可靠的数据支撑，且近年来，随着分析技术的发展，钻孔岩心元素分析数据也得到了丰富地积累。
现有技术中主要通过勘探线地质剖面与地球化学插值图件对钻孔勘探资料进行处理呈现，其通常遵循单工程钻孔柱状图→勘探线剖面图→分段平面图的处理顺序。该方法存在以下明显缺陷：
(1)在构建勘探线地质/地球化学剖面时，由于人工勾连的主观性和插值方法的畸变或者过拟合，往往偏离地质实际；(2)不同类型的地质和地化图件边界差别较大，往往难以进行空间综合匹配分析，只能进行粗略的叠加分析；(3)由于剖面在两处钻孔之间缺少竖直方向上的约束，其常用的处理顺序存在先天性的误差，而剖面绘制阶段产生的误差会在绘制平面图时进一步累积，导致失真。
另一方面，近年来，研究者逐渐认识到，地质系统，尤其是岩浆、构造、流体复合作用的体系，是一个具有自组织临界性的开放、复杂、模糊的非线性系统。此外，内生矿床，尤其是流体相关矿床，往往形成自不同地质历史时期多个系统的复杂叠加。试图将其作为整体进行解译在现有技术条件下难以达成。不论是根据经验进行推断，还是根据现有数据进行插值，均有可能产生误导。
发明内容
本发明的目的在于提出一种改进的矿床钻探信息处理方法和装置，其可利用计算机技术，将分段分层的地质钻孔信息等进行可视化地解析呈现，该处理方法可解决现有方法中图件成果难以进行多元地质信息与地球化学综合对比分析的问题，并可进一步解决解决传统的通过勘探线剖面图获得分段平面图方法中，由于钻孔之间缺少竖直方向上的约束而存在系统性误差的问题及勘探剖面地质体界线圈定人为因素影响较大，地化图件绘制插值畸变明显等问题。
本发明的技术方案如下：
基于有限单元的矿床钻探信息处理方法，其包括：
获得待分析矿床的钻孔信息、其采样样品信息、及其钻孔对应的地质信息，其中，所述采样样品信息包括通过钻孔进行采样的样品的位置信息及其化学分析信息；
基于矿产地质勘查规范行业标准所规定的伴生品位、边界品位、工业品位三个标准值，根据所述采样样品信息对不同位置样品对应的矿床进行品级划分；
根据所述钻孔信息、其采样样品信息及钻孔对应的地质信息对待分析矿床进行垂直分段解析，所述解析包括：确定其垂直分段的间距和数量，计算各分段的矿化、隐爆、蚀变参数，对计算得到的各参数进行分级；
获得不同垂直段内勘查钻孔控制点的平面投影图，并根据各勘查点控制范围进行图件区划，形成分段有限单元区划图；
根据所述垂直分段解析结果，为不同垂直段的所述分段有限单元区划图内的各有限单元进行所述不同参数的赋值，形成各参数的分段有限单元平面图；
将所述各参数的分段有限平面图按照其空间位置进行三维叠加，叠加中，同位置的参数进行组合，得到三维复合数据图；
对所述三维复合数据图进行竖直剖面上的扫描解析，获得其对应的含有复合数据的有限单元剖面图。
根据本发明的一些优选实施方式，所述钻孔信息包括钻孔的位置坐标，钻孔的高程，钻孔的方位角、倾角和孔深。
根据本发明的一些优选实施方式，所述采样样品信息包括采样的起始位置坐标，采样的终止位置坐标，样品的化学分析项目和其化学分析结果数据。
根据本发明的一些优选实施方式，所述地质信息包括钻孔对应的岩矿石名称、地质描述、蚀变类型和隐爆破碎程度。
根据本发明的一些优选实施方式，所述品级划分包括：A级：≥伴生品位～＜边界品位；B级：≥边界品位～＜工业品位；C级：≥工业品位。
根据本发明的一些优选实施方式，所述垂直分段解析具体包括：
设置分段解析的垂直间距
根据所述钻孔信息，获得钻孔的最大开孔的高程H(max)、最小终孔高程H(min)，并结合所述垂直间距确定分段数量[N]+1，其中，[]表示取整，N计算如下：
根据钻孔对应的地质信息，计算各垂直段内钻孔的不同所述品级的矿化厚度百分比参数、隐爆段厚百分比参数、蚀变段厚百分比参数，并将各百分比参数从高到低划分为若干等级；
其中，各百分比参数计算如下：
矿化厚度百分比＝(该垂直分段内矿化厚度/样品段厚度)×100％；
隐爆段厚百分比＝(该垂直分段内隐爆段长度/样品段厚度)×100％；
蚀变段厚百分比＝(该垂直分段内蚀变段长度/样品段厚度)×100％。
根据本发明的一些优选实施方式，所述垂直间距为20～25m。
根据本发明的一些优选实施方式，所述等级包括如下的5级：Ⅰ级为100～80％，Ⅱ级为80～50％，Ⅲ级为50～30％，Ⅳ级为30～10％，Ⅴ级为10～0％。
根据本发明的一些优选实施方式，所述分段有限单元区划图的获得具体包括：
在任一垂直分段内，将其所述勘查钻孔控制点投影到平面图上；
将最外围各控制点进行有限外推，并连接外推圆的包络线，确定有限单元区划图边界；
依据邻点共边原则，沿相邻工程之间的中点划分有限单元的界线，据此划分各有限单元，形成分段有限单元区划图；
根据区域地质构造、工程偏斜程度与不连续面分布情况，结合地质专家知识库，调整有限单元边界划分情况；
将所述垂直分段解析的结果逐类赋值于各有限单元，形成各参数的分段有限单元平面图。
根据本发明的一些优选实施方式，所述外推的间距为矿区工程平均间距的1/2。
根据本发明的一些优选实施方式，所述分段有限单元平面图中，不同点的数据组成如下：
对任一点Dm，
当其位于有限单元内时，有：
Dm＝{Xm，Ym，P，αP，βP，γP…}；
当其位于有限单元边界上时，有：
Dm＝{Xm，Ym，Hn，B}；
其中，Xm与Ym表示该点平面坐标，P代表该点在平面上所属的有限单元，αP，βP，γP…代表该有限单元被赋值的待进行分析的各类属性，B代表该点的边界点性质，Hn表示该点所在分段的高程。
根据本发明的一些优选实施方式，所述含有复合数据的有限单元剖面图的获得具体包括：
在所述三维复合数据图的各水平面上，布设水平间距为的平行竖直剖面，与所述分段有限单元区划图进行区-线空间判别分析；
当所述竖直剖面与所述分段有限单元区划图的分段平面相交时，根据所述同位置的参数的组合对其交点进行赋值；
基于所述分段有限单元平面图中赋值的边界点将所述竖直剖面与所述分段平面的交线进行分割，并将该交线上各边界点沿竖直剖面向上或向下延拓后，水平封口，划分剖面上的有限单元；
基于交线上的赋值为(P，αP，βP，γP)的属性点，对剖面有限单元赋值，构建解析的有限单元剖面图和其对应的数据库。
根据上述信息处理方法，可进一步获得一种分析方法，其包括：对所述各参数的分段有限单元平面图和/或解析剖面图进行叠合分析。
根据本发明的一些优选实施方式，所述分析方法还包括：根据叠合分析结果，判定多个对象在空间上的相互关系，结合矿化等级、隐爆等级和蚀变等级的发展趋势与地质规律，圈定潜在成矿的远景区带，对矿区周边及深部矿体资源潜力进行预测评价。
根据上述信息处理方法和/或分析方法，可进一步获得一种进行矿床钻探信息分析的装置，其包括存储介质和运算器，其中，所述存储介质存储有上述任一信息处理方法和/或分析方法的程序和/或结构数据。
本发明具备以下有益效果：
通过本发明的处理或分析方法可获得成套的分段面状投影解析平面图和分级面状扫描解析剖面图，并实现：
(1)规范地处理钻孔勘探地质和地球化学信息，减少人为解译的主观性和插值分析的畸变，更加客观、真实地反映钻孔范围内地下的地层岩性、构造、矿化、蚀变、隐爆破碎等信息，所呈现的图件更为贴合地质实际；
(2)过构建仅与钻孔控制的空间分布相关的有限单元，并将多元信息赋值于外形规范的有限单元，便于针对多期叠加的复杂系统进行地质、地球化学多维分析；
(3)通过不同深度分段有限单元平面图→空间三维同位叠合→扫描解析有限单元剖面图，利用平面上钻孔控制点分布较剖面图上更为均匀的特点，减少传统方法中单工程钻孔柱状图→勘探线剖面图→分段平面图解译方法中，由于钻孔之间缺少竖直方向的约束而带来的误差。
本发明系统性解决了当前传统地质钻孔信息表达方法内容少、用途单一、成矿预测效果差的问题，对精准圈定找矿靶区、精准采矿选冶、构造成矿学研究等方面具有优异的实用价值和指示意义。
附图说明
图1为本发明的一种具体实施方式的流程示意图。
图2为本发明具体实施方式中扫描解析有限单元剖面图的流程示意图。
图3为本发明具体实施方式步骤5中采用全角半距法对分段平面图上的斜孔控制点坐标进行校正的示意图。
图4为本发明实施例1中所得某矿区铜、钼矿化特征分段有限单元平面图。
图5为本发明实施例1中所得某矿区铜、钼矿化特征扫描解析有限单元南北向剖面图。
图6为本发明实施例1中所得某矿区铜、钼矿化特征扫描解析有限单元东西向剖面图。
具体实施方式
以下结合实施例和附图对本发明进行详细描述，但需要理解的是，所述实施例和附图仅用于对本发明进行示例性的描述，而并不能对本发明的保护范围构成任何限制。所有包含在本发明的发明宗旨范围内的合理的变换和组合均落入本发明的保护范围。
以基于钻孔矿化、隐爆(在地下岩层中发生的，与高压流体相关的爆破现象，通常与流体活动和矿化关联)和蚀变信息圈定有利成矿区为例，对本发明的一些具体实施方式进行说明。
参照说明书附图1～2，根据本发明的技术方案，一种具体的实施过程可包括如：
步骤1以一个地质工作项目或矿床为单位，全面收集与地质钻孔资料相关的图件资料、文字资料和表格资料，对收集的地质钻孔资料进行系统的检查和整理。
其中，收集的各资料可主要包括：地质工作报告、工程布置图、钻孔取样品位(基本)分析数据表、钻孔工程参数(如位置参数及其测斜参数)、钻孔柱状图和勘探线剖面图。
检查内容可包括：
(1)完整性检查：检查地质钻孔坐标是否齐全，钻孔数据量是否足够，单工程各分段是否连续，钻孔取样分析数据是否完整。
(2)有效性检查：检查钻孔基本化学分析数据是否有未检出和错误数值数据，分析数据项之间的逻辑关系是否正确。
(3)规范性检查：检查岩矿石名称、地质描述、蚀变类型、破碎情况数据是否规范和科学。
步骤2完成步骤1后，对钻孔的基本信息、化学分析信息、编录信息进行整理。
其中，钻孔的基本信息可包括钻孔施工参数，如开孔的位置坐标X、Y，开孔的高程H；及测斜参数，如钻孔的方位角、倾角和孔深。
钻孔的基本化学分析信息可包括钻孔中的采样位置，如采样起始位置和采样终止位置；分析项名称和分析结果。
编录信息可包括岩矿石名称、地质描述、蚀变类型、隐爆破碎程度。
以上各基本信息可进一步通过如下的钻孔基本信息表(表1)，钻孔基本化学分析信息表(表2)，和钻孔编录信息表(表3)进行填写整理：
表1钻孔基本信息表
表2钻孔基本化学分析信息表
表3钻孔编录信息表
序号 钻孔编号 分层号 起孔深 止孔深 岩矿石名称 地质描述 蚀变类型 隐爆破碎程度
1
2
3
步骤3完成步骤2后，依据矿产地质勘查规范行业标准所规定的伴生品位、边界品位、工业品位三个标准值，将钻孔基本化学分析数据从高到低划分三个品级。
更具体的，所述三个品级可包括：A级：≥伴生品位～＜边界品位；B级：≥边界品位～＜工业品位；C级：≥工业品位。
步骤4确定分段解析的垂直间距和分段数量，计算各分段的矿化、隐爆、蚀变分级参数。
更具体的，其可包括：
(1)设置分段解析的垂直间距一般可优选为25m。
(2)根据钻孔最大开孔高程H(max)、最小终孔高程H(min)和分段解析垂直间距确定分段数量为[N]+1，其中，[]表示取整，N计算如下：
(4)根据钻孔编录信息和步骤3所划分的三个品级，计算各分段钻孔分级(A级或B级或C级)矿化厚度百分比、隐爆段厚百分比、蚀变段厚百分比，并将百分比数据从高到低划分Ⅰ级～Ⅴ级若干个等级，按下表进行填写：
表4分段矿化、隐爆、蚀变分级统计表
其中，具体参数的计算如下：
矿化厚度百分比＝(段内矿化厚度/样品段厚)×100％；
隐爆段厚百分比＝(段内隐爆段长度/样品段厚)×100％；
蚀变段厚百分比＝(段内蚀变段长度/样品段厚)×100％。
分级可设置如：
Ⅰ级为100～80％，Ⅱ级为80～60％，Ⅲ级为60～40％，Ⅳ级为40～10％，Ⅴ级为10～0％。
步骤5在投有控制点位置的勘查控制平面图上，依相邻点控制区共边原则，确定反映各勘查点控制范围区划的图件，获得分段有限单元区划图，并根据步骤4分段量化统计得到的各项成果，为区划图内各有限单元赋值，形成分段有限单元平面图(以下简称分段平面图)。
其可更具体的，包括：
(1)将勘查钻孔控制点投影到各分段平面图上。
向下钻进的勘查钻孔可分为斜孔和垂直孔，其中，垂直钻孔在分段平面图上的控制点坐标Xi、Yi即为步骤2所述的开孔位置坐标X、Y；斜孔在分段平面图上的控制点坐标，可依据步骤2中的钻孔施工参数、测斜参数和各分段平面图的标高不同，以开孔位置坐标X、Y为已知起始点，采用全角半距法对分段平面图上的控制点坐标Xi、Yi进行校正和分深度段赋值，如附图3所示。
全角半距法可更具体的，为：
式中，Xi、Yi为斜孔在各分段平面图上的控制点坐标(i＝0，1，2，…，N)；Li为各分段平面图上的控制点到开孔位置X、Y的长度；Zi为各分段平面图上控制点处的倾角；Mi为各分段平面图上控制点处的方位角；dM为磁偏角。Li、Zi、Mi可在步骤2中测斜参数中查取；dM可在钻孔所在区域1/5万地形图内查取。
(2)将最外围各控制点进行有限外推，并连接外推圆的包络线，确定有限单元区划图边界，其中外推间距可选如矿区工程平均间距的1/2。
(3)依据邻点共边原则，在沿相邻工程之间的中点划分有限单元的界线，据此划分各有限单元，形成分段有限单元区划图；分段有限单元区划图最小比例尺可一般采用如1/25000。
(4)根据区域地质构造、工程偏斜程度与不连续面分布情况，结合地质专家知识库，适度调整有限单元边界划分情况。
(5)将分段统计表中的结果逐类赋值于各有限单元，形成分段有限单元平面图与相应数据库。
进一步的，由于分段平面图有限单元划分只与钻孔控制点位置相关，而不受其他变量影响，故此一处有限单元可以容纳多组信息，以分段平面图上某一点点Dm为例，说明该数据库的结构，如下：
对有限单元内部任一点Dm，有：
Dm＝{Xm，Ym，p，αP，βP，γP…}
其中，Xm与Ym为该点平面坐标；P代表该点在平面上所属的有限单元，在分段有限单元划定之后，可以视为固定值，αP，βP，γP…等代表该有限单元被赋值的欲进行分析的各类属性。
对有限单元边界任一点，则有：
Dm＝{Xm，Ym，Hn，B}
其中B代表该点的边界点性质，Hn为该点所在分段的高程。
根据所得分段有限单元平面图及数据库，即可对多类地质信息进行单独或综合分析呈现，例如：分段分级矿化厚度有限单元平面图、分段隐爆段厚-蚀变矿化叠加有限单元平面图、分段矿化厚度-隐爆段厚-蚀变段厚叠加有限单元平面图等。
步骤6将各类数据的各分段平面图按照其空间位置输入三维空间数据库后，构成同位叠合的复合数据库。
更具体的，下面以三维空间数据库中自上而下第n分段某一点Dm为例，说明该数据库的结构,若该点位于某一有限单元内部，则：
Dm＝{Xm，Ym，Hn，P，αP，βP，γP…}
其中，Hn为该点所在分段的高程(以该分段的垂直间距中点为准)其满足：
H(max)≥Hn≥H(min)；
若该点位于有限单元的边界，则：
Dm＝{Xm，Ym，Hn，B}。
步骤7使用Mapgis等GIS软件，构建相对应的扫描解析有限单元剖面图(以下简称解析剖面图)。
更具体的，其可包括：
(1)布设水平间距为(如一般可选择为25m或50m)的平行竖直剖面，与步骤6中同位叠合的分段有限单元区划图进行区-线空间判别分析。
(2)当竖直剖面与分段平面相交时，则基于步骤6中建立的同位叠合数据库对其进行赋值，以某一竖直剖面Ll(Xl，Yl，Hl)为例，当其与步骤6中提及的第n分段Dm点重合时，即当(Xl＝Xm，Yl＝Ym，Hl＝Hn)时，则为剖面上重合点赋值为(P，αP，βP，γP)或(B)。
(3)基于上一步骤赋值的边界点(即所有赋值为(B)的点)将竖直剖面与分段平面的交线分割。此后，将该交线上各边界点沿竖直剖面向上或向下延拓后，水平封口，划分剖面上的有限单元。
(4)基于交线上的赋值为(P，αP，βP，γP)的属性点，对剖面有限单元赋值，构建解析剖面数据库。
步骤8使用Mapgis等GIS软件多元面状空间分析与叠合功能，对不同类型数据的分段平面图或解析剖面图进行叠合分析，判定多个对象在空间上的相互关系，结合矿化等级、隐爆等级和蚀变等级的发展趋势与一般地质规律，圈定潜在成矿的远景区带，对矿区周边及深部矿体资源潜力进行预测评价。
上述过程中，步骤5-7扫描解析有限单元剖面图的过程如附图2所示。
实施例1
根据以上具体实施方式，以某铜钼矿床为例，全面收集与地质钻孔相关的资料，对收集的地质钻孔资料进行系统的检查和整理，矿床共施工钻孔78个，钻孔总进尺数米数为25427.72m，共采集铅基本化学分析样2796件，其中采样段进尺共13074.93m，占钻孔总进尺米数的51.42％。
为精细刻划Cu、Mo的单元素矿化分级，依据矿产地质勘查规范行业标准所规定的伴生品位、边界品位、工业品位三个标准值，将钻孔Cu、Mo两种矿化元素的基本化学分析数据从高到低划分三个品级(表5)。
表5某铜钼矿床矿化元素分级简表
通过对钻孔工程参数和基本分析数据的处理，确定分段解析的垂直间距为25m，将78个勘查钻孔控制点投影到25m间距的各分段平面图上，并计算不同分段的矿化分级参数，并将矿化厚度百分比从高到低划分为Ⅰ级～Ⅴ级(表6)。
表6某铜钼矿床分段矿化分级统计表
在投有控制点位置的勘查控制平面图上，依据相邻点控制区共边原则，确定反映各勘查点控制范围区划的图件，获得分段有限单元区划图，并根据分段统计得到的各项成果，为区划图内各有限单元赋值，得到如附图4所示的Cu、Mo元素A级矿化分段平面图。
将各类数据的各分段平面图按照其空间位置输入三维空间数据库后，构成同位叠合的复合数据库，使用Mapgis等GIS软件，得到如附图5、6所示的Cu、Mo元素A级矿化解析剖面图。
使用Mapgis等GIS软件多元面状空间分析与叠合功能，对分段平面图或解析剖面图进行叠合分析。
其中，附图4分段平面图显示，矿体中段平面形态多呈星点状、网状、条带状、团块状，Cu、Mo元素矿化厚度等级呈现随机性，中部矿化向网状、团块状聚拢，南部边界出现逐渐离散的现象，局部包裹无矿单元。总体来看，Cu、Mo元素分级矿化平面形态多样，反映成矿作用具有多期次、多类型热液的参与。
附图5、6解析剖面图显示，在南部向剖面上，由西向东的矿体形态呈叉状或纺锤状，垂向上矿化由条带状-团块状-网格状变化；在东西向剖面上，矿床铜钼矿化在南边呈明显的梳状、网状形态，由南至北，矿化体的规模逐渐变小，南部、西部矿化比北部、东部矿化更好，形态也逐渐变成规模不大的梯状，反映了横向上的热液叠加或浸漫作用。由此可见，在垂向上，矿化在深部更为集中，范围更广，强度更大。
以上实施例仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例。凡属于本发明思路下的技术方案均属于本发明的保护范围。应该指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下的改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (11)

1.基于有限单元的矿床钻探信息处理方法，其特征在于，其包括：
获得待分析矿床的钻孔信息、其采样样品信息、及其钻孔对应的地质信息，其中，所述采样样品信息包括通过钻孔进行采样的样品的位置信息及其化学分析信息；
基于矿产地质勘查规范行业标准所规定的伴生品位、边界品位、工业品位三个标准值，根据所述采样样品信息对不同位置样品对应的矿床进行品级划分；
根据所述钻孔信息、其采样样品信息及钻孔对应的地质信息对待分析矿床进行垂直分段解析，所述解析包括：确定其垂直分段的间距和数量，计算各分段的矿化、隐爆、蚀变参数，对计算得到的各参数进行分级；
获得不同垂直分段内勘查钻孔控制点的平面投影图，并根据各勘查点控制范围进行图件区划，形成分段有限单元区划图；
根据所述垂直分段解析结果，为不同垂直分段的所述分段有限单元区划图内的各有限单元进行不同参数的赋值，形成各参数的分段有限单元平面图；
将所述各参数的分段有限平面图按照其空间位置进行三维叠加，叠加中，同位置的参数进行组合，得到三维复合数据图；
对所述三维复合数据图进行竖直剖面上的扫描解析，获得其对应的含有复合数据的有限单元剖面图；
其中，所述垂直分段解析具体包括：
设置分段解析的垂直间距
根据所述钻孔信息，获得钻孔的最大开孔的高程H(max)、最小终孔高程H(min)，并结合所述垂直间距确定分段数量[N]+1，其中，[]表示取整，N计算如下：
根据钻孔对应的地质信息，计算各垂直分段内钻孔的不同所述品级的矿化厚度百分比参数、隐爆段厚百分比参数、蚀变段厚百分比参数，并将各百分比参数从高到低划分为若干等级；
其中，各百分比参数计算如下：
矿化厚度百分比＝(该垂直分段内矿化厚度/样品段厚度)×100％；
隐爆段厚百分比＝(该垂直分段内隐爆段长度/样品段厚度)×100％；
蚀变段厚百分比＝(该垂直分段内蚀变段长度/样品段厚度)×100％。
2.根据权利要求1所述的信息处理方法，其特征在于，其中，所述钻孔信息包括钻孔的位置坐标，钻孔的高程，钻孔的方位角、倾角和孔深；所述采样样品信息包括采样的起始位置坐标，采样的终止位置坐标，样品的化学分析项目和其化学分析结果数据；所述地质信息包括钻孔对应的岩矿石名称、地质描述、蚀变类型和隐爆破碎程度。
3.根据权利要求1所述的信息处理方法，其特征在于，所述品级划分包括：A级：≥伴生品位～＜边界品位；B级：≥边界品位～＜工业品位；C级：≥工业品位。
4.根据权利要求1所述的信息处理方法，其特征在于，其中，所述垂直间距为20～25m；和/或，所述等级包括如下的5级：Ⅰ级为100～80％，Ⅱ级为80～60％，Ⅲ级为60～40％，Ⅳ级为40～10％，Ⅴ级为10～0％。
5.根据权利要求1所述的信息处理方法，其特征在于，所述分段有限单元区划图的获得具体包括：
在任一垂直分段内，将其所述勘查钻孔控制点投影到平面图上；
将最外围各控制点进行有限的外推，并连接外推圆的包络线，确定有限单元区划图边界；
依据邻点共边原则，沿相邻工程之间的中点划分有限单元的界线，据此划分各有限单元，形成分段有限单元区划图；
根据区域地质构造、工程偏斜程度与不连续面分布情况，结合地质专家知识库，调整有限单元边界划分情况；
将所述垂直分段解析的结果逐类赋值于各有限单元，形成各参数的分段有限单元平面图。
6.根据权利要求5所述的信息处理方法，其特征在于，所述外推的间距为矿区工程平均间距的1/2。
7.根据权利要求5所述的信息处理方法，其特征在于，所述分段有限单元平面图中，不同点的数据组成如下：
对任一点Dm，
当其位于有限单元内时，有：
Dm＝{xm，Ym，P，αP，βP，γP…}；
当其位于有限单元边界上时，有：
Dm＝{Xm，Ym，Hn，B}；
其中，Xm与Ym表示该点平面坐标，P代表该点在平面上所属的有限单元，αP，βP，γP…代表该有限单元被赋值的待进行分析的各类属性，B代表该点的边界点性质，Hn表示该点所在分段的高程。
8.根据权利要求7所述的信息处理方法，其特征在于，所述含有复合数据的有限单元剖面图的获得具体包括：
在所述三维复合数据图的各水平面上，布设水平间距为的平行竖直剖面，与所述分段有限单元区划图进行区-线空间判别分析；
当所述竖直剖面与所述分段有限单元区划图的分段平面相交时，根据所述同位置的参数的组合对其交点进行赋值；
基于所述分段有限单元平面图中赋值的边界点将所述竖直剖面与所述分段平面的交线进行分割，并将该交线上各边界点沿竖直剖面向上或向下延拓后，水平封口，划分剖面上的有限单元；
基于交线上的赋值为(P，αP，βP，γP)的属性点，对剖面有限单元赋值，构建解析的有限单元剖面图和其对应的数据库。
9.应用权利要求1-8中任一项所述的信息处理方法进行矿床钻探分析的分析方法，其包括：对所述各参数的分段有限单元平面图和/或解析剖面图进行叠合分析。
10.根据权利要求9所述的分析方法，其特征在于，其还包括：根据叠合分析结果，判定多个对象在空间上的相互关系，结合矿化等级、隐爆等级和蚀变等级的发展趋势与地质规律，圈定潜在成矿的远景区带，对矿区周边及深部矿体资源潜力进行预测评价。
11.一种进行矿床钻探信息分析的装置，其包括存储介质和运算器，所述存储介质存储有实现权利要求1-8中任一项所述的信息处理方法或权利要求9或10所述的分析方法的程序和/或结构数据。
CN202111586366.8A 2021-12-23 2021-12-23 基于有限单元的矿床钻探信息处理及分析方法和装置 Active CN114266179B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111586366.8A CN114266179B (zh) 2021-12-23 2021-12-23 基于有限单元的矿床钻探信息处理及分析方法和装置
US18/066,605 US20230204810A1 (en) 2021-12-23 2022-12-15 Finite element-based ore deposit drilling information processing and analysis method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111586366.8A CN114266179B (zh) 2021-12-23 2021-12-23 基于有限单元的矿床钻探信息处理及分析方法和装置

Publications (2)

Publication Number Publication Date
CN114266179A CN114266179A (zh) 2022-04-01
CN114266179B true CN114266179B (zh) 2023-08-04

Family

ID=80828944

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111586366.8A Active CN114266179B (zh) 2021-12-23 2021-12-23 基于有限单元的矿床钻探信息处理及分析方法和装置

Country Status (2)

Country Link
US (1) US20230204810A1 (zh)
CN (1) CN114266179B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117057633B (zh) * 2023-10-13 2024-01-02 北京科技大学 一种确定低品位多金属矿品位指标的方法和系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6480790B1 (en) * 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
CN103700141A (zh) * 2013-11-19 2014-04-02 中国矿业大学(北京) 煤矿底板突水通道的多级尺度地质建模分析方法
CN104102972A (zh) * 2014-06-30 2014-10-15 南京梅山冶金发展有限公司 一种基于爆破单元体的矿山资源储量核销方法
CN104536054A (zh) * 2015-01-08 2015-04-22 中国科学院地质与地球物理研究所 一种条带状铁矿评估方法和装置
CN107346038A (zh) * 2017-06-08 2017-11-14 昆明理工大学 “四步式”大比例尺定位探测深部热液矿床或矿体的方法
CN110685737A (zh) * 2019-09-02 2020-01-14 河北省矾山磷矿有限公司 平底结构分段凿岩阶段出矿嗣后充填采矿法
CN113338935A (zh) * 2021-04-21 2021-09-03 铜陵有色金属集团股份有限公司 一种圈定矿体的方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102713142B (zh) * 2009-08-14 2015-12-16 Bp北美公司 储层构型和连通性分析
US8453764B2 (en) * 2010-02-01 2013-06-04 Aps Technology, Inc. System and method for monitoring and controlling underground drilling
DE102018130569B4 (de) * 2018-11-30 2021-04-29 Airbus Operations Gmbh System und Verfahren zur Darstellung eines Kabinenlayouts im Originalmaßstab

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6480790B1 (en) * 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
CN103700141A (zh) * 2013-11-19 2014-04-02 中国矿业大学(北京) 煤矿底板突水通道的多级尺度地质建模分析方法
CN104102972A (zh) * 2014-06-30 2014-10-15 南京梅山冶金发展有限公司 一种基于爆破单元体的矿山资源储量核销方法
CN104536054A (zh) * 2015-01-08 2015-04-22 中国科学院地质与地球物理研究所 一种条带状铁矿评估方法和装置
CN107346038A (zh) * 2017-06-08 2017-11-14 昆明理工大学 “四步式”大比例尺定位探测深部热液矿床或矿体的方法
CN110685737A (zh) * 2019-09-02 2020-01-14 河北省矾山磷矿有限公司 平底结构分段凿岩阶段出矿嗣后充填采矿法
CN113338935A (zh) * 2021-04-21 2021-09-03 铜陵有色金属集团股份有限公司 一种圈定矿体的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
采掘工作面底板突水和防治原则的基本理论研究;曹玉清;华北地质矿产杂志(第03期);全文 *

Also Published As

Publication number Publication date
CN114266179A (zh) 2022-04-01
US20230204810A1 (en) 2023-06-29

Similar Documents

Publication Publication Date Title
US7853045B2 (en) Geostatistical analysis and classification of core data
EP2795527B1 (en) Systems and methods for selecting facies model realizations
CN101738639A (zh) 提高岩石裂缝参数计算精度的方法
CN111062544A (zh) 一种铀成矿远景区的预测方法
Bertotti et al. Toward a quantitative definition of mechanical units: New techniques and results from an outcropping deep-water turbidite succession (Tanqua-Karoo Basin, South Africa)
CN114266179B (zh) 基于有限单元的矿床钻探信息处理及分析方法和装置
CN115983505A (zh) 一种固体矿产三维立体成矿预测方法及装置
Shamseddin Meigooni et al. Application of multivariate geostatistical simulation and fractal analysis for detection of rare-earth element geochemical anomalies in the Esfordi phosphate mine, Central Iran
Abzalov et al. Exploratory data analysis
Leung et al. Sample truncation strategies for outlier removal in geochemical data: the MCD robust distance approach versus t-SNE ensemble clustering
Ovaskainen et al. A new subsampling methodology to optimize the characterization of two-dimensional bedrock fracture networks
Guo et al. Identification of homogeneous region boundaries of fractured rock masses in candidate sites for Chinese HLW repository
Sadeghi SPS Model: a significant algorithm to reduce the time and computer memory required in geostatistical simulations
Poniewierski Block model knowledge for mining engineers–an introduction
Carvalho Probabilistic resource modeling of vein deposits
CN110702632B (zh) 一种深部岩矿高光谱信息三维建模方法
Hu A drawing system for pole diagram and rose diagram of rock mass structural surface
La Pointe et al. Geological discrete fracture network model for the Laxemar site. Site Descriptive Modelling. SDM-Site Laxemar
CN107085650B (zh) 三维地震项目评价方法及装置
Yunsel Risk quantification in grade variability of gold deposits using sequential Gaussian simulation
CN117471546B (zh) 一种黑色岩系型金矿找矿方法
Tokoglu Comparative Analysis of 3D Domain Modelling Alternatives: Implications for Mineral Resource Estimates
AU2019373096B2 (en) A method and system for sample classification
Amwaama Comparison of kriging neighbourhood analysis (KNA) results from two different software packages: A case study
Du Toit A Refined Resource Model for Tshepong Mine

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