CN114663627B - 一种基于三维点云数据库的矿山数字模型建立方法 - Google Patents

一种基于三维点云数据库的矿山数字模型建立方法 Download PDF

Info

Publication number
CN114663627B
CN114663627B CN202210312210.9A CN202210312210A CN114663627B CN 114663627 B CN114663627 B CN 114663627B CN 202210312210 A CN202210312210 A CN 202210312210A CN 114663627 B CN114663627 B CN 114663627B
Authority
CN
China
Prior art keywords
mine
model
point cloud
rock stratum
coordinate set
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
CN202210312210.9A
Other languages
English (en)
Other versions
CN114663627A (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202210312210.9A priority Critical patent/CN114663627B/zh
Publication of CN114663627A publication Critical patent/CN114663627A/zh
Priority to PCT/CN2023/084075 priority patent/WO2023185734A1/zh
Priority to GB2400439.2A priority patent/GB2623664A/en
Priority to US18/580,031 priority patent/US20240346758A1/en
Application granted granted Critical
Publication of CN114663627B publication Critical patent/CN114663627B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21C—MINING OR QUARRYING
    • E21C41/00—Methods of underground or surface mining; Layouts therefor
    • E21C41/16—Methods of underground mining; Layouts therefor
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22—Indexing; Data structures therefor; Storage structures
    • G06F16/2291—User-Defined Types; Storage management thereof
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/28—Databases characterised by their database models, e.g. relational or object models
    • G06F16/283—Multi-dimensional databases or data warehouses, e.g. MOLAP or ROLAP
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G06F16/903—Querying
    • G06F16/90335—Query processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V99/00—Subject matter not provided for in other groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10028—Range image; Depth image; 3D point clouds
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/56—Particle system, point based geometry or rendering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Computer Graphics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Mining & Mineral Resources (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于三维点云数据库的矿山数字模型建立方法，属于矿山三维数字建模领域方法。收集矿山数据资料，依据点云数据提取标准，提取有效坐标点位置信息，建立矿山三维点云数据库；构建岩层分界面、采区及巷道模型，构建地表三维地质模型，进而构建矿山三维地质模型；划分网格，得到矿山采掘反演数字模型。依据该方法可获得矿山三维点云数据库，并能够还原各标识岩层地貌、采区和巷道位置、断层形态及地表形态，根据仿真模拟需求，利用本方法能够建立不同精度的矿山采掘反演数字模型，为地应力场反演奠定基础，为矿井灾害防治提供可靠的依据。

Description

一种基于三维点云数据库的矿山数字模型建立方法
技术领域
本发明涉及一种基于三维点云数据库的矿山数字模型建立方法，属于矿山三维数字建模技术领域。
背景技术
随着地下矿产资源的不断开采，开采深度和开采强度逐渐加大，矿山动力灾害一旦发生将对人员安全和矿井生产造成严重影响。矿压观测是目前解决采动问题的主要决策依据，影响巷道矿压控制、开采程序设计、巷道位置合理选择与维护、冲击地压和突出的预测和防治、承压水上煤体安全开采等重大技术问题的决策部署。矿压观测受限于测试位置和成本高等问题，很难在矿区进行大规模的监测，仿真模拟技术的发展解决了实测受限的问题，可以为现场实践提供全面科学的指导。
仿真模拟的关键之一是构建模型，模型的准确性是模拟结果可靠性的重要依据，其中对于大型复杂模型的建立一直都是学者们关注的重中之重。前人所建模型有三点不足：1)模型较小，多为局部采区模型；2)模型简单，忽略了实际煤层褶曲、断层等地质构造；3)前人建立的矿山模型很少采用专业的建模软件，而多数模拟软件内置的建模模块较为简单，无法满足复杂地质条件下的建模需要。本发明选用专业的3D建模软件Rhino，可以较高的还原煤层实际样貌，提高与现场实际的贴合性，满足复杂大型矿山建模的需要。
点云模型往往由测量直接得到，每个点对应一个测量点，未经过其他处理手段，故包含了最大的信息量。然而，这些信息隐藏在点云中需要以其他提取手段将其萃取出来，提取点云中信息的过程则为三维图像处理。
三维点云建模思想已经在多个领域应用，多采用非接触的方式采集目标对象的表面属性点信息，获得三维点云数据，而对目标对象的内部结构不做细致建模。目前，在采矿领域，三维激光扫描技术应用于露天矿山、边坡等地上模型构建，对地下复杂模型构建未有详细研究。为此本发明详细介绍了一种复杂矿山采掘反演数字模型的建立方法，其中重点介绍三维点云建模思想在地下模型构建的应用。
发明内容
根据现有技术的不足之处，提供一种基于三维点云数据库的矿山数字模型建立方法，通过建立矿山三维点云数据库，进而利用三维建模软件，建立能够反映岩层起伏和延伸情况、采区和巷道位置、地质构造的矿山采掘反演数字模型。
为实现本发明的目的，本发明基于三维点云数据库的矿山数字模型建立方法，其特征在于：首先收集矿山数据资料，依据点云数据提取标准，提取有效岩层界面的坐标点位置信息，建立矿山三维点云数据库；然后利用矿山的三维点云数据构建岩层分界面、采区及巷道模型，构建地表的三维地质模型，进而构建矿山三维地质模型；在矿山三维地质模型上划分网格，得到还原了各标识岩层地貌、采区和巷道位置、断层形态及地表形态的矿山采掘反演数字模型，提高了模型的贴合性和准确性。
具体步骤如下：
S1、通过收集目标矿区地质资料，确定模型尺寸，并分析获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、地表等高线空间坐标集、断层位置、断层断距及断层角度；
S2、处理各岩层上下分界面空间坐标集，提取有效坐标点位置信息，建立矿山三维点云数据库：
S2.1、根据模拟需求确定模型精度ε(0＜ε＜1)；
S2.2、确定矿山三维点云数据提取标准和岩层分界面空间坐标集被提取的最小重要系数；
设共有m组岩层分界面空间坐标集和n条矿山三维点云数据提取标准，利用下式计算第k组岩层分界面空间坐标集的重要系数Pk
Figure BDA0003567516010000021
Figure BDA0003567516010000022
其中，Xi为n条矿山三维点云数据提取标准中标准i的重要系数，wi为标准i的权重，l为修正系数；
根据模拟需求确定岩层分界面空间坐标集被提取的最小重要系数Pmin且满足
Figure BDA0003567516010000023
Pmax为最大重要系数；
S2.3、将满足Pk≥Pmin的各岩层分界面空间坐标集直接提取到矿山三维点云数据库中并分别命名；
S2.4、将剩余岩层分界面空间坐标集分为两类，即满足
Figure BDA0003567516010000024
或
Figure BDA0003567516010000025
的坐标集，其中α为可调系数且α∈(0,1)，然后将符合
Figure BDA0003567516010000026
的岩层删去，将其与附近满足
Figure BDA0003567516010000027
的岩层合并形成新的岩层分界面坐标集，并将其提取至矿山三维点云数据库中并分别命名；
S3、构建还原断层形态、采区及巷道的三维地质模型的地下部分：
S3.1、将包含各标识岩层上下分界面空间坐标集导入Rhino软件，利用Rhino软件由下到上的分层拟合岩层分界面；由模型尺寸建立长方体模型，并由各拟合得到的岩层分界面进行切割，对各岩层进行命名，此时初步得到长方体结构的地下三维地质模型；
S3.2、根据断层位置，建立断层分界线并根据断层角度，切割煤层；然后根据断层断距，调整被分割的煤层还原断层形态；
S3.3、将包含回采区段拐点空间坐标集、关键巷道拐点空间坐标集输入Rhino软件，分别拟合出工作面各回采区段和关键巷道并命名；
S3.4、对岩层空间模型、回采区段及关键巷道进行布尔运算，去除岩层空间模型中的多余部分，形成全矿区三维地质模型的地下部分；
S4、根据地表等高线空间坐标集，构建地表三维地质模型并命名，将两部分模型拼接在一起，由此全矿区三维地质模型建立完成；
S5、利用Rhino软件自带的网格划分工具对全矿区三维地质模型进行初步网格划分，网格尺寸根据实际需要设置，然后利用Griddle插件对模型再次进行网格划分，通过限制网格形状和最大、最小边缘尺寸来修正网格，使网格变化均匀；最后将模型导出为FLAC3D软件可识别的网格文件，在该模拟软件中导入网格文件，就得到全矿井采掘反演数字模型。
在步骤S1中收集的地质资料包括目标矿区采掘工程平面图、工程地质剖面图、地勘钻孔柱状图、煤层底板等高线图。通过分析目标矿区的地质资料，获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、地表等高线空间坐标集、模型尺寸、断层位置、断层角度、断层断距。关键巷道主要指区段回风平巷、区段运输平巷。
在步骤S1中所述的各岩层上下分界面空间坐标集是由地勘钻孔柱状图分析获得。每个钻孔的纬距和经距一定，可以直接获得该组钻孔坐标集的横坐标X和纵坐标Y，根据岩层数量p和埋藏深度hj可以获得该组钻孔坐标集，由上到下分别为(X,Y,h1)、(X,Y,h2)、……(X,Y,hj)、……(X,Y,hp+1)，其中hj≤hp+1。将获得的所有坐标先按照地质单位“段”重新进行分类，进一步的在每“段”中，再按照不同岩层进行分类，这样便得到各岩层上下分界面空间坐标集。
在步骤S2中所述的模型精度ε是根据模拟需求确定的仿真模型还原精细程度，既可以是局部某些连续岩层还原精度，也可以是整体地层还原精度。地质资料越丰富，地勘钻孔打钻越密集，模型最大还原精度越高，由于无法做到100％还原矿山地下样貌，所以ε＜1。
在步骤S2中所述的矿山三维点云数据提取标准为根据仿真模拟需求、矿山地下岩层形态、地勘钻孔柱状图和岩层分界面空间坐标集确定的详细要求，该要求的精细程度直接影响矿山三维点云数据库的建立。
在步骤S3中所述的标识岩层指矿山三维点云数据库中确定的岩层，模型精度越高且矿山三维点云数据提取标准越详细，标识岩层数量越多。在漫长的地质时期，由老到新逐渐形成的各标识岩层，可以在地质考查、钻孔勘探和矿井采掘过程中用于判别层位信息和地质年代。
在步骤S3中，由于回采区段和关键巷道分别与岩层空间模型存在空间重复，故利用软件布尔运算功能中的差集运算，对岩层空间模型、回采区段及关键巷道进行差集运算，目的是在岩层空间模型中切除回采区段及关键巷道，避免空间位置重复，这样岩层空间模型和回采区段及关键巷道就构成了全矿区三维地质模型的地下部分。
有益效果：1)本方法将三维点云建模思想应用于地层，依据该方法可获得矿山三维点云数据库；2)本方法既可以最高限度的还原煤层实际地质构造、岩层起伏和延伸情况及地表形态，又可以根据仿真模拟需求建立不同精度的矿山采掘反演数字模型。采用本方法建模并进行地应力场反演，可以为矿井开拓、巷道布置和支护等现场实践提供更加全面科学的指导。
附图说明
图1为本发明基于三维点云数据库的矿山数字模型建立方法的流程图。
图2为本发明实施例通过点云数据库创建流程图。
图3为本发明实施例等比例三维地质模型示意图。
具体实施方式
下面结合附图对本发明的实施例做进一步说明；
实施例一、
如图1所示，一种基于三维点云数据库的矿山数字模型建立方法，包括以下步骤：
S1收集数据资料
收集目标矿区采掘工程平面图、工程地质剖面图、地勘钻孔柱状图、煤层底板等高线图，确定模型尺寸，并分析获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、断层位置、断层断距及断层角度。
其中，各岩层上下分界面空间坐标集是由地勘钻孔柱状图分析获得。每个钻孔的纬距和经距一定，可以直接获得该组钻孔坐标集的横坐标X和纵坐标Y，根据岩层数量p和埋藏深度hj可以获得该组钻孔坐标集，由上到下分别为(X,Y,h1)、(X,Y,h2)、…(X,Y,hj)、…(X,Y,hp+1)，其中hj≤hp+1。将获得的所有坐标先按照地质单位“段”重新进行分类，进一步的在每“段”中，再按照不同岩层进行分类，这样便得到各岩层上下分界面空间坐标集。
S2建立矿山三维点云数据库；
1)根据模拟需求确定模型精度ε(ε＜1)；
2)确定矿山三维点云数据提取标准和岩层分界面空间坐标集被提取的最小重要系数；
设共有m组岩层分界面空间坐标集和n条矿山三维点云数据提取标准，Pk为第k组岩层分界面空间坐标集的重要系数，Xi为n条矿山三维点云数据提取标准中标准i的重要系数，wi为标准i的权重，l为修正系数，则有：
Figure BDA0003567516010000051
其中，
Figure BDA0003567516010000052
根据模拟需求确定岩层分界面空间坐标集被提取的最小重要系数Pmin，且满足
Figure BDA0003567516010000053
3)将满足Pk≥Pmin的各岩层分界面空间坐标集直接提取到矿山三维点云数据库中并分别命名；
4)将剩余岩层分界面空间坐标集分为两类，即满足
Figure BDA0003567516010000054
或
Figure BDA0003567516010000055
的坐标集，其中α为可调系数且α∈(0,1)，然后将符合
Figure BDA0003567516010000056
的岩层删去，将其与附近满足
Figure BDA0003567516010000057
的岩层合并，重新获取岩层分界面坐标集，并将其提取至矿山三维点云数据库中并分别命名。
S3构建还原断层形态、采区及巷道的三维地质模型的地下部分；
1)将包含各标识岩层上下分界面空间坐标集导入Rhino软件，利用Rhino软件由下到上的分层拟合岩层分界面；由模型尺寸建立长方体模型，并由各拟合得到的岩层分界面进行切割，对各岩层进行命名，此时初步得到地下三维地质模型；
2)根据断层位置，建立断层分界线并根据断层角度，切割煤层；然后根据断层断距，调整被分割的煤层，还原断层形态；
3)将包含回采区段拐点空间坐标集、关键巷道拐点空间坐标集输入Rhino软件，分别拟合出工作面各回采区段和关键巷道并命名；
4)对岩层空间模型、回采区段及关键巷道进行布尔运算，去除岩层空间模型中的多余部分，形成全矿区三维地质模型的地下部分；
S4构建矿山三维地质模型
根据地表等高线空间坐标集，构建地表三维地质模型并命名，将两部分模型拼接在一起，由此全矿区三维地质模型建立完成；
S5划分网格，得到矿山采掘反演数字模型
利用Rhino软件自带的网格划分工具对全矿区三维地质模型进行初步网格划分，网格尺寸根据实际需要设置，然后利用Griddle插件对模型再次进行网格划分，通过限制网格形状和最大、最小边缘尺寸来修正网格，使网格变化均匀；最后将模型导出为某模拟软件可识别的网格文件，在该模拟软件中导入网格文件，就得到全矿井采掘反演数字模型。
本发明提供了一种基于三维点云数据库的矿山数字模型建立方法，下面将以某矿山为实施例，详细说明本发明提供的方法。
S1收集数据资料
收集目标矿区采掘工程平面图、工程地质剖面图、地勘钻孔柱状图、煤层底板等高线图，确定模型尺寸，并分析获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、断层位置、断层断距及断层角度。
其中，各岩层上下分界面空间坐标集是由地勘钻孔柱状图分析获得。每个钻孔的纬距和经距一定，可以直接获得该组钻孔坐标集的横坐标X和纵坐标Y，根据岩层数量p和埋藏深度hj可以获得该组钻孔坐标集，由上到下分别为(X,Y,h1)、(X,Y,h2)、……(X,Y,hj)、……(X,Y,hp+1)，其中hj≤hp+1。将获得的所有坐标先按照地质单位“段”重新进行分类，进一步的在每“段”中，再按照不同岩层进行分类，最终获得m组(m≥p+1)岩层分界面空间坐标集。
S2图2为本实施例三维点云数据库创建流程图，根据图中流程建立该矿山三维点云数据库；
1)根据模拟需求确定模型精度为0.4，可调系数且α为0.8；
2)确定矿山三维点云数据提取标准和岩层分界面空间坐标集被提取的最小重要系数；
设共有4条矿山三维点云数据提取标准，Pk为第k组岩层分界面空间坐标集的重要系数，Xi为4条矿山三维点云数据提取标准中标准i的重要系数，标准1为岩层厚度情况，w1＝0.4；标准2为层位稳定情况，w2＝0.3；标准3为分布范围，w3＝0.2；标准4为钻孔密集程度，w4＝0.1，具体X1、X2、X3、X4的取值如表1、表2、表3和表4所示，修正系数l为0.5，则有：
Figure BDA0003567516010000071
根据模拟需求确定岩层分界面空间坐标集被提取的最小重要系数Pmin为7，其中，Pmax＝10.5，满足
Figure BDA0003567516010000072
Figure BDA0003567516010000073
即Pk≥3.36；
Figure BDA0003567516010000074
即Pk＜3.36。
表1岩层厚度情况和X1取值
岩层厚度 0.5m及以下 0.5m～5m 5m～10m 10m～20m 20m以上
<![CDATA[X<sub>1</sub>]]> 1 2 5 7 10
表2层位稳定情况和X2取值
层位稳定情况 非常不稳定 较不稳定 一般 较稳定 非常稳定
<![CDATA[X<sub>2</sub>]]> 1 2 5 7 10
表3分布范围和X3取值
分布范围 30％以下 30％～50％ 50％～70％ 70％～90％ 90％以上
<![CDATA[X<sub>3</sub>]]> 1 3 5 7 10
表4钻孔密度和X4取值
钻孔密集程度 很分散 较分散 一般 较密集 很密集
<![CDATA[X<sub>4</sub>]]> 1 3 5 7 10
3)将满足Pk≥7的各岩层分界面空间坐标集直接提取到矿山三维点云数据库中并分别命名；
4)将剩余岩层分界面空间坐标集分为两类，即满足Pk≥3.36或Pk＜3.36的坐标集，然后将符合Pk＜3.36的岩层删去，将其与附近满足Pk≥3.36的岩层合并，重新获取岩层分界面坐标集，并将其提取至矿山三维点云数据库中并分别命名。
S3构建还原断层形态、采区及巷道的三维地质模型的地下部分；
1)将包含各标识岩层上下分界面空间坐标集导入Rhino软件，利用Rhino软件由下到上的分层拟合岩层分界面；由模型尺寸建立长方体模型，并由各拟合得到的岩层分界面进行切割，对各岩层进行命名，此时初步得到地下三维地质模型；
2)根据断层位置，建立断层分界线并根据断层角度，切割煤层；然后根据断层断距，调整被分割的煤层，还原断层形态；
3)将包含回采区段拐点空间坐标集、关键巷道拐点空间坐标集输入Rhino软件，分别拟合出工作面各回采区段和关键巷道并命名；
4)对岩层空间模型、回采区段及关键巷道进行布尔运算，去除岩层空间模型中的多余部分，形成全矿区三维地质模型的地下部分；
S4构建矿山三维地质模型
根据地表等高线空间坐标集，构建地表三维地质模型并命名，将两部分模型拼接在一起，由此全矿区三维地质模型建立完成；
S5划分网格，得到矿山采掘反演数字模型
利用Rhino软件自带的网格划分工具对全矿区三维地质模型进行初步网格划分，网格尺寸根据实际需要设置，然后利用Griddle插件对模型再次进行网格划分，通过限制网格形状和最大、最小边缘尺寸来修正网格，使网格变化均匀；最后将模型导出为某模拟软件可识别的网格文件，在该模拟软件中导入网格文件，就得到如图3所示的全矿井采掘反演数字模型。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (7)

1.一种基于三维点云数据库的矿山数字模型建立方法，其特征在于：首先收集矿山数据资料，依据点云数据提取标准，提取有效岩层界面的坐标点位置信息，建立矿山三维点云数据库；然后利用矿山的三维点云数据构建岩层分界面、采区及巷道模型，构建地表的三维地质模型，进而构建矿山三维地质模型；在矿山三维地质模型上划分网格，得到还原了各标识岩层地貌、采区和巷道位置、断层形态及地表形态的矿山采掘反演数字模型，提高了模型的贴合性和准确性；
具体步骤如下：
S1、通过收集目标矿区地质资料，确定模型尺寸，并分析获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、地表等高线空间坐标集、断层位置、断层断距及断层角度；
S2、处理各岩层上下分界面空间坐标集，提取有效坐标点位置信息，建立矿山三维点云数据库：
S2.1、根据模拟需求确定模型精度ε，0＜ε＜1；
S2.2、确定矿山三维点云数据提取标准和岩层分界面空间坐标集被提取的最小重要系数；
设共有m组岩层分界面空间坐标集和n条矿山三维点云数据提取标准，利用下式计算第k组岩层分界面空间坐标集的重要系数Pk
Figure FDA0004155246200000011
Figure FDA0004155246200000012
其中，Xi为n条矿山三维点云数据提取标准中标准i的重要系数，wi为标准i的权重，l为修正系数；
根据模拟需求确定岩层分界面空间坐标集被提取的最小重要系数Pmin且满足
Figure FDA0004155246200000013
Pmax为最大重要系数；
S2.3、将满足Pk≥Pmin的各岩层分界面空间坐标集直接提取到矿山三维点云数据库中并分别命名；
S2.4、将剩余岩层分界面空间坐标集分为两类，即满足
Figure FDA0004155246200000014
或
Figure FDA0004155246200000015
的坐标集，其中α为可调系数且α∈(0,1)，然后将符合
Figure FDA0004155246200000016
的岩层删去，将其与附近满足
Figure FDA0004155246200000021
的岩层合并形成新的岩层分界面坐标集，并将其提取至矿山三维点云数据库中并分别命名；
S3、构建还原断层形态、采区及巷道的三维地质模型的地下部分：
S3.1、将包含各标识岩层上下分界面空间坐标集导入Rhino软件，利用Rhino软件由下到上的分层拟合岩层分界面；由模型尺寸建立长方体模型，并由各拟合得到的岩层分界面进行切割，对各岩层进行命名，此时初步得到长方体结构的地下三维地质模型；
S3.2、根据断层位置，建立断层分界线并根据断层角度，切割煤层；然后根据断层断距，调整被分割的煤层还原断层形态；
S3.3、将包含回采区段拐点空间坐标集、关键巷道拐点空间坐标集输入Rhino软件，分别拟合出工作面各回采区段和关键巷道并命名；
S3.4、对岩层空间模型、回采区段及关键巷道进行布尔运算，去除岩层空间模型中的多余部分，形成全矿区三维地质模型的地下部分；
S4、根据地表等高线空间坐标集，构建地表三维地质模型并命名，将两部分模型拼接在一起，由此全矿区三维地质模型建立完成；
S5、利用Rhino软件自带的网格划分工具对全矿区三维地质模型进行初步网格划分，网格尺寸根据实际需要设置，然后利用Griddle插件对模型再次进行网格划分，通过限制网格形状和最大、最小边缘尺寸来修正网格，使网格变化均匀；最后将模型导出为FLAC 3D软件可识别的网格文件，在该模拟软件中导入网格文件，就得到全矿井采掘反演数字模型。
2.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S1中收集的地质资料包括目标矿区采掘工程平面图、工程地质剖面图、地勘钻孔柱状图、煤层底板等高线图；通过分析目标矿区的地质资料，获取各岩层上下分界面空间坐标集、回采区段拐点空间坐标集、关键巷道拐点空间坐标集、地表等高线空间坐标集、模型尺寸、断层位置、断层角度、断层断距；关键巷道指区段回风平巷、区段运输平巷。
3.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S1中所述的各岩层上下分界面空间坐标集是由地勘钻孔柱状图分析获得；每个钻孔的纬距和经距一定，可以直接获得该组钻孔坐标集的横坐标X和纵坐标Y，根据岩层数量p和埋藏深度hj可以获得该组钻孔坐标集，由上到下分别为(X,Y,h1)、(X,Y,h2)、……(X,Y,hj)、……(X,Y,hp+1)，其中hj≤hp+1；将获得的所有坐标先按照地质单位“段”重新进行分类，进一步的在每“段”中，再按照不同岩层进行分类，这样便得到各岩层上下分界面空间坐标集。
4.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S2中所述的模型精度ε是根据模拟需求确定的仿真模型还原精细程度，既是局部某些连续岩层还原精度，也是整体地层还原精度；地质资料越丰富，地勘钻孔打钻越密集，模型最大还原精度越高，由于无法做到100％还原矿山地下样貌，所以ε＜1。
5.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S2中所述的矿山三维点云数据提取标准为根据仿真模拟需求、矿山地下岩层形态、地勘钻孔柱状图和岩层分界面空间坐标集确定的详细要求，该要求的精细程度直接影响矿山三维点云数据库的建立。
6.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S3中所述的标识岩层指矿山三维点云数据库中确定的岩层，模型精度越高且矿山三维点云数据提取标准越详细，标识岩层数量越多；在漫长的地质时期，由老到新逐渐形成的各标识岩层，可以在地质考查、钻孔勘探和矿井采掘过程中用于判别层位信息和地质年代。
7.根据权利要求1所述的基于三维点云数据库的矿山数字模型建立方法，其特征在于：在步骤S3中，由于回采区段和关键巷道分别与岩层空间模型存在空间重复，故利用软件布尔运算功能中的差集运算，对岩层空间模型、回采区段及关键巷道进行差集运算，目的是在岩层空间模型中切除回采区段及关键巷道，避免空间位置重复，这样岩层空间模型和回采区段及关键巷道就构成了全矿区三维地质模型的地下部分。
CN202210312210.9A 2022-03-28 2022-03-28 一种基于三维点云数据库的矿山数字模型建立方法 Active CN114663627B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202210312210.9A CN114663627B (zh) 2022-03-28 2022-03-28 一种基于三维点云数据库的矿山数字模型建立方法
PCT/CN2023/084075 WO2023185734A1 (zh) 2022-03-28 2023-03-27 一种基于三维点云数据库的矿山数字模型建立方法
GB2400439.2A GB2623664A (en) 2022-03-28 2023-03-27 Mine digital model construction method based on three-dimensional point cloud database
US18/580,031 US20240346758A1 (en) 2022-03-28 2023-03-27 Mine digital model construction method based on three-dimensional point cloud database

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210312210.9A CN114663627B (zh) 2022-03-28 2022-03-28 一种基于三维点云数据库的矿山数字模型建立方法

Publications (2)

Publication Number Publication Date
CN114663627A CN114663627A (zh) 2022-06-24
CN114663627B true CN114663627B (zh) 2023-05-16

Family

ID=82034098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210312210.9A Active CN114663627B (zh) 2022-03-28 2022-03-28 一种基于三维点云数据库的矿山数字模型建立方法

Country Status (4)

Country Link
US (1) US20240346758A1 (zh)
CN (1) CN114663627B (zh)
GB (1) GB2623664A (zh)
WO (1) WO2023185734A1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114663627B (zh) * 2022-03-28 2023-05-16 中国矿业大学 一种基于三维点云数据库的矿山数字模型建立方法
CN117197398B (zh) * 2023-11-07 2024-02-09 中国电建集团昆明勘测设计研究院有限公司 一种矿山资源模型的统计方法、装置、设备及存储介质
CN117371068A (zh) * 2023-11-09 2024-01-09 山西工程技术学院 一种对bif铁矿顺层边坡变形破坏的数值模拟分析方法
CN117456122B (zh) * 2023-12-22 2024-03-08 江西省地质博物馆 一种三维地质模型生成方法及系统
CN117788740B (zh) * 2023-12-29 2024-08-30 武汉永业赛博能规划勘测有限公司 基于三维模型的计算方法、系统和存储介质
CN117574518B (zh) * 2024-01-15 2024-04-09 国网湖北省电力有限公司经济技术研究院 一种抽水蓄能电站地下厂房三维地质模型建模方法及系统
CN117574687B (zh) * 2024-01-15 2024-04-02 垒知科技集团有限公司 一种煤气柜施工建设中活塞运行仿真模拟方法
CN117576335B (zh) * 2024-01-19 2024-04-26 长安大学 面向矿产区域勘查的三维空间模型数据处理方法及系统
CN117593477B (zh) * 2024-01-19 2024-04-19 中国煤炭地质总局勘查研究总院 融合无人机三维模型和地质填图的三维地质图填绘方法
CN118036125A (zh) * 2024-01-19 2024-05-14 北京科技大学 一种地下矿山透明采场的构建方法和三维可视化模型
CN117671160B (zh) * 2024-01-30 2024-04-23 北京星天地信息科技有限公司 多源数据协同的煤层建模方法、装置、设备及存储介质
CN118036287B (zh) * 2024-02-05 2024-08-23 赣州有色冶金研究所有限公司 一种矿山开采区域模拟方法及系统
CN117830586B (zh) * 2024-03-05 2024-06-11 煤炭科学研究总院有限公司 一种井巷三维建模增量更新方法及系统
CN118096703B (zh) * 2024-03-08 2024-10-18 中国地质科学院 一种基于重磁电三维模型的找矿靶区优选方法及相关装置
CN118128527B (zh) * 2024-05-06 2024-09-17 中国科学技术大学 深部岩体随钻激光测量地应力的测量系统及测量方法
CN118505921B (zh) * 2024-07-17 2024-10-11 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) 一种地质档案的模型三维可视化方法及其系统
CN118657902B (zh) * 2024-08-20 2024-10-18 临沂大学 新近系地层结构与薄基岩开采工程地质模型构建方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622526A (zh) * 2012-04-06 2012-08-01 中国矿业大学(北京) 数字矿山中的巷道开掘搜索预测方法
CN109712241A (zh) * 2018-12-17 2019-05-03 长沙矿山研究院有限责任公司 一种包含采空区的三维矿山实体模型的建立方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6608913B1 (en) * 2000-07-17 2003-08-19 Inco Limited Self-contained mapping and positioning system utilizing point cloud data
US8224631B2 (en) * 2008-08-18 2012-07-17 Fci Holdings Delaware, Inc. Stress, geologic, and support analysis methodology for underground openings
CN110118550B (zh) * 2019-05-14 2021-02-05 中国矿业大学 一种综采工作面高精度三维导航地图的生成系统及方法
CN110163966B (zh) * 2019-05-16 2022-12-09 中煤科工集团北京华宇工程有限公司 利用点云数据自动建立煤炭矿井三维地质信息模型的方法
CN113379909B (zh) * 2021-05-08 2023-10-27 华北科技学院(中国煤矿安全技术培训中心) 一种透明工作面智能开采大数据分析决策方法和系统
CN114005051A (zh) * 2021-09-27 2022-02-01 辽宁工程技术大学 一种基于图像识别的井巷通风阻力系数计算方法
CN113962008A (zh) * 2021-11-09 2022-01-21 华北科技学院(中国煤矿安全技术培训中心) 一种矿山三维井巷模型生成及透明工作面构建方法
CN114663627B (zh) * 2022-03-28 2023-05-16 中国矿业大学 一种基于三维点云数据库的矿山数字模型建立方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622526A (zh) * 2012-04-06 2012-08-01 中国矿业大学(北京) 数字矿山中的巷道开掘搜索预测方法
CN109712241A (zh) * 2018-12-17 2019-05-03 长沙矿山研究院有限责任公司 一种包含采空区的三维矿山实体模型的建立方法

Also Published As

Publication number Publication date
WO2023185734A1 (zh) 2023-10-05
US20240346758A1 (en) 2024-10-17
GB2623664A (en) 2024-04-24
GB202400439D0 (en) 2024-02-28
CN114663627A (zh) 2022-06-24

Similar Documents

Publication Publication Date Title
CN114663627B (zh) 一种基于三维点云数据库的矿山数字模型建立方法
Guo et al. Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm
Kaufmann et al. 3D geological modelling from boreholes, cross-sections and geological maps, application over former natural gas storages in coal mines
CN113821953B (zh) 基于动态地应力和剩余油特征的暂堵重复压裂设计方法
CN103700141B (zh) 煤矿底板突水通道的多级尺度地质建模分析方法
CN107515423B (zh) 一种高陡构造的三维地质模型构建方法
CN113313825B (zh) 一种油藏储层三维地质建模方法和装置
CN112446956B (zh) 一种基于钻孔数据的地质三维模型构建方法
Zhang et al. A case study on integrated modeling of spatial information of a complex geological body
CN113689559B (zh) 提高地质三维模型内容表达精度的网格状图切地质剖面布线方法
CN105321204A (zh) 一种三维地质栅格体模型构建方法
CN109712239A (zh) 一种矿床精细三维地质建模方法
CN113742827A (zh) 一种基于有限差分分析的公路边坡监测网络体系构建方法
CN113536414B (zh) 基于三维建模的岩质边坡稳定性分析方法、系统及介质
CN110197013B (zh) 基于Morphing的河床基岩面建模方法
CN110705168A (zh) 构造应力场的模拟方法
CN108256259B (zh) 基于模糊层次分析确定煤层底板注浆靶区的方法及装置
Kaufmann et al. Reprint of “3D geological modelling from boreholes, cross-sections and geological maps, application over former natural gas storages in coal mines”[Comput. Geosci. 34 (2008) 278–290]
CN111859687A (zh) 一种刻画含铀砂层地质结构的混合地质建模方法及系统
CN114820959B (zh) 一种冲击地压矿井多地质要素数字化建模方法
CN115713599A (zh) 一种地下采煤矿山精确监测防治方法
CN115035258A (zh) 一种高效的基于cad钻孔柱状图的城市三维地质建模方法
CN111612900B (zh) 一种回采工作面煤层透明化三维地质模型构建方法及装置
CN116129063A (zh) 一种煤矿地质三维模型的构建方法及系统
CN112949106A (zh) 一种岩土工程地质地表移动变形状态的检测方法

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