CN109326002A - 基于钻孔数据的矿体建模方法、装置、***及存储介质 - Google Patents
基于钻孔数据的矿体建模方法、装置、***及存储介质 Download PDFInfo
- Publication number
- CN109326002A CN109326002A CN201811421277.6A CN201811421277A CN109326002A CN 109326002 A CN109326002 A CN 109326002A CN 201811421277 A CN201811421277 A CN 201811421277A CN 109326002 A CN109326002 A CN 109326002A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- ore
- point
- distance field
- constraint
- 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.)
- Granted
Links
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于钻孔数据的矿体建模方法、装置、***及存储介质，该方法包括：对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；根据所述接触点约束和所述离面点约束构建离散点距离场；对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。可以直接基于钻孔数据进行矿体三维建模，省去了人工解译矿体剖面轮廓线导致的交互繁琐的麻烦，有效提高了矿体建模的效率和建模质量。
Description
技术领域
本发明涉及矿体三维建模领域，具体涉及一种基于钻孔数据的矿体建模方法、装置、***及存储介质。
背景技术
现有的矿体三维建模一般采用基于人工解译剖面线框进行轮廓线拼接或三角网重构的矿体显式建模方法。由于可以将矿体限制在一定的形状或岩性中，该矿体显式建模方法在早期一直被认为是“最佳”的行业实践。然而，该矿体显式建模方法具有交互过程繁杂、无法对逐步精细化的模型进行动态更新、生成的模型质量差等先天性缺陷。故亟需设计一种新的矿体建模方法，以有效改善矿体建模的效率和建模质量。
发明内容
有鉴于此，本发明实施例提供了一种基于钻孔数据的矿体建模方法、装置、***及存储介质，旨在有效改善矿体建模的效率和建模质量。
本发明实施例的技术方案是这样实现的：
本发明实施例第一方面，提供一种基于钻孔数据的矿体建模方法，该方法包括：
对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；
根据所述接触点约束和所述离面点约束构建离散点距离场；
对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；
基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
进一步地，所述对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，包括：
按照设定的条件对所述钻孔数据进行预处理得到钻孔数据样品组合；
在所述钻孔数据样品组合的基础上按钻孔轨迹对样段和非样段进行离散化采样得到离散采样点。
进一步地，所述根据所述接触点约束和所述离面点约束构建离散点距离场，包括：
校验并修正所述离面点约束；
根据所述接触点约束和修正后的所述离面点约束构建所述离散点距离场。
进一步地，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，包括：
根据所述离散点距离场和所述离散采样点构建插值约束，将所述插值约束代入插值方程进行插值计算得到与插值计算后的插值方程对应的所述隐式函数距离场。
进一步地，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，还包括：
添加附加插值趋势约束；
根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
进一步地，所述添加插值趋势约束，包括：
添加插值趋势约束点、插值趋势约束线和/或插值趋势面信息；
所述根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场，包括：
根据所述插值趋势约束点、所述插值趋势约束线和/或所述插值趋势面信息，对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
进一步地，所述基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型，包括：
根据所述插值约束提取用于界定模型边界的最小外包；
根据所述最小外包确定模型三维重构的分辨率；
按照所述分辨率对所述最小外包进行空间分割，构造离散化的空间规则数据场；
对所述隐式函数距离场按照所述离散化的空间规则数据场进行采样计算，重构所述三维模型。
本发明实施例的第二方面，还提供一种基于钻孔数据的矿体建模装置，该装置包括：
采样模块，用于对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；
第一距离场构建模块，用于根据所述接触点约束和所述离面点约束构建离散点距离场；
第二距离场构建模块，用于对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；
建模模块，用于基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
本发明实施例的第三方面，还提供一种终端设备，包括存储器、处理器及存储在所述存储器上并能够由所述处理器运行的可执行程序，所述处理器运行所述可执行程序时执行前述任一实施例所述的基于钻孔数据的矿体建模方法。
本发明实施例的第四方面，还提供一种存储介质，其上存储有可执行程序，其特征在于，所述可执行程序被处理器执行时实现前述任一项所述的基于钻孔数据的矿体建模方法。
本发明实施例的技术方案中，通过对地质勘探采样的钻孔数据进行离散化采样，提取样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束，根据所述接触点约束和所述离面点约束构建离散点距离场，对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型，可以直接基于钻孔数据进行矿体三维建模，省去了人工解译矿体剖面轮廓线导致的交互繁琐的麻烦，有效提高了矿体建模的效率和建模质量。
附图说明
图1为本发明实施例基于钻孔数据的矿体建模方法的流程示意图；
图2为本发明实施例中钻孔数据样品组合的结构示意图；
图3为本发明实施例中钻孔数据离散化的结构示意图；
图4为本发明实施例中异常离面点约束的结构示意图；
图5为本发明实施例中修正后的离面点约束的结构示意图；
图6为本发明实施例中1#钻孔数据矿体建模过程示意图；
图7为本发明实施例中2#钻孔数据矿体建模对应的不同边界品位下的品位壳模型示意图；
图8为本发明实施例基于钻孔数据的矿体建模装置的结构示意图；
图9为本发明实施例基于钻孔数据的矿体建模***的硬件结构示意图。
具体实施方式
以下结合说明书附图及具体实施例对本发明技术方案做进一步的详细阐述。应当理解，此处所提供的实施例仅仅用以解释本发明，并不用于限定本发明。另外，以下所提供的实施例是用于实施本发明的部分实施例，而非提供实施本发明的全部实施例，在不冲突的情况下，本发明实施例记载的技术方案可以任意组合的方式实施。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本申请。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
请参阅图1，图1为本发明一实施例中基于钻孔数据的矿体建模方法，该方法包括：
步骤101，对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束。
这里，地质勘探采样的钻孔数据是指地质勘探中，对于地下的矿体通过打钻孔的方式来采样的钻孔数据。在采样过程中，对每个钻孔划分不同的分段，有矿且品位足够高的分段认为是矿体，该分段就是样段。样段以外的分段就是非样段。
在一可选的实施方式中，所述对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，包括：
按照设定的条件对所述钻孔数据进行预处理得到钻孔数据样品组合；
在所述钻孔数据样品组合的基础上按钻孔轨迹对样段和非样段进行离散化采样得到离散采样点。
本实施例中，预处理是指利用采集的钻孔数据创建包括孔口表、测斜表、样品表及岩性表的钻孔数据库。类似于地质解译过程，在钻孔数据库校验的基础上，选定组合样长、边界品位、工业品位和最小可采厚度等参数，应用“穿鞋戴帽”的方法按地质解译要求对钻孔数据进行样品组合(参见图2)，使样长参数达到可采厚度的要求，为离散化采样提供数据基础。
本实施例中，可选地，在所述钻孔数据样品组合的基础上按钻孔轨迹对样段和非样段进行离散化采样得到离散采样点之前，还包括：
确定采样点间距；提取钻孔数据所有样段和非样段信息，以最小样段长度为最大间距，以模型重构精度为最小间距，按每个样段至少采样一次的原则沿着钻孔轨迹方向确定采样间距。
本实施例中，可以采用径向基函数RBF插值方法，通过提取样段的采样点和非样段的采样点将钻孔线框数据转化为接触点约束与离面点约束条件。图3示出了本发明实施例钻孔数据离散化的结构示意图。
其中，提取样段的采样点并转化为接触点约束包括：以所有样段两端边界点为矿体边界点，进而建立接触点约束条件(f(x)＝0)，表示相应约束点的函数值为零。沿着钻孔轨迹方向，对样段按采样点间距从两端向中间进行离散化采样，按距离邻近矿体边界点样段长度赋以函数值，建立模型内部离面点约束条件；为保证模型的可靠性，样段内至少采样一次。
其中，提取非样段的采样点并转化为离面点约束包括：对非样段按采样点间距从两端向中间进行离散化采样，按距离邻近矿体边界点非样段长度赋以函数值，用于建立模型内外离面点约束条件(f(x)＞0或f(x)＜0)，其中，函数值大于零对应模型外离面约束点，函数值小于零对应模型内离面约束点。为保证模型的可靠性，非样段内也至少采样一次。由于缺乏邻接边界信息，对于没有任何样品段的钻孔不再进行离散化采样和约束。
步骤102，根据所述接触点约束和所述离面点约束构建离散点距离场。
本实施例中，所述根据所述接触点约束和所述离面点约束构建离散点距离场，包括：
校验并修正所述离面点约束；
根据所述接触点约束和修正后的所述离面点约束构建所述离散点距离场。
请参阅图4和图5，为了避免重复、异常、二义性的距离函数值，在构建离散点距离场之前，本实施例根据模型内外离面点约束应符合最近距离的特征规则来检查离面点约束的二义性，并对异常的离面约束点按照上述规则进行修正。然后基于钻孔数据构建非欧式距离场，将钻孔轨迹长度作为非欧式距离场的距离，从而可以构造更加符合按钻孔轨迹趋势分布的地质距离的离散点距离场。
步骤103，对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场。
在一可选的实施方式中，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，包括：
根据所述离散点距离场和所述离散采样点构建插值约束，将所述插值约束代入插值方程进行插值计算得到与插值计算后的插值方程对应的所述隐式函数距离场。本实施例中，可以根据样段的采样点和/或非样段的采样点在所述离散点距离场中的距离值构建插值约束，该插值约束为各采样点对应的距离值。
本实施例中，插值方程如下：
其中，x＝(x,y,z)为任意待求插值点；f(x)为径向基函数(简称为RBF)；p(x)为低阶多项式；λi为径向基函数的系数；φ为严格正定或条件正定的核函数；ci为径向基函数的插值中心点。
可选地，在插值计算前，设置一定的距离容差，采用快速空间索引方法来过滤重复约束，以保证插值方程求解的唯一性。对于插值方程的求解，可以采用FMM快速多极化的方法加速求解过程。
在另一可选的实施方式中，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，还包括：
添加附加插值趋势约束；
根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
为了处理复杂条件下矿体各向异性的特性，利用隐式函数便于动态更新的优点，在隐式模型自动构建的基础上可以通过添加插值趋势约束点、插值趋势约束线和/或插值趋势面等信息对隐式模型进行插值约束，动态更新数学模型，建立符合矿化域分布趋势和地质工程师人工解译约束的三维矿体隐式模型。
在一可选的实施方式中，添加附加插值趋势约束点进行动态更新。通过插值趋势约束点可以改变局部稀疏不均匀数据的插值趋势和局部地质界面的连续性，构建符合稀疏钻孔数据地质特征的矿体模型。
在一可选的实施方式中，添加附加插值趋势约束线进行动态更新。采用附加法向约束剖面的约束线，通过限定约束线所在切平面方向构建法向信息。约束线所在切平面限定了约束线的法向所在平面，根据约束线各分段方向，分段确定约束线法向信息；再结合切平面所切割模型的内外关系确定法向的正方向(指向模型外部)。对插值趋势约束线按一定间距进行离散化采样，并由采样点沿采样点法向的正反两个方向偏移一定距离得到离面约束点，将插值趋势约束线转化为插值趋势约束点，以便构造RBF插值约束条件。
在一可选的实施方式中，对于添加的插值趋势约束线，需要采用最近点约束校验方法来处理处理二义性的离面约束点。根据法向偏离附加的离面点约束应符合最近距离特征来检查离面点约束的二义性，以避免造成病态约束矩阵。对于存在二义性的离面点约束，在选定迭代次数内，通过逐步迭代减小法向离面间距的方法来修正离面点约束。对于超过迭代次数的二义性离面点约束，则不再定义法向离面约束点，单独定义接触点约束。
在一可选的实施方式中，为了处理矿体按特定矿化域方向延伸趋势的特征，可以通过构造插值趋势面的方法来处理矿体模型的各向异性插值趋势要求。矿体模型的各向异性通常表现为几何各向异性，其中主轴、次轴、短轴方向趋势面各向异性表现为各向异性椭球体，轴长表现为各向异性的强度大小。利用几何各向异性的线性变换关系，可以通过一定的矩阵变换把几何各向异性转换为各向同性。先按各向异性椭球体建立新的坐标系，再按主轴、次轴、短轴方向各向异性强度E1、E2、E3(E1<E2<E3)构建缩放矩阵S进行缩放，形成各向同性的坐标系空间，来计算原始坐标系下的各向异性距离。因此，各向异性插值问题可以转化为各向同性插值问题。
插值趋势面非常适用于稀疏钻孔数据在较大距离上沿一定平面方向连续的建模情况。为了提高给定方向的插值连续性，可以增大插值权重在该方向上的相关性，而与最强各向异性方向对齐的点的各向异性距离小于其各向同性距离。
步骤104，基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
本实施例中，所述基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型，包括：
根据插值约束提取用于界定模型边界的最小外包；
根据所述最小外包确定模型三维重构的分辨率；
按照所述分辨率对所述最小外包进行空间分割，构造离散化的空间规则数据场；
对所述隐式函数距离场按照所述离散化的空间规则数据场进行采样计算，得到矿体的三维模型。
本实施例中，采用等值面提取(Marching Cubes)算法对隐式函数按移动立方体进行采样计算，重构矿体三维模型。
本实施例矿体建模方法，根据样品组合结果提取钻孔数据，离散化样段与非样段信息，构建插值数据接触点约束和离面点约束，不需要人工解译矿体剖面轮廓线，可以直接基于钻孔数据进行矿体三维建模。省去了人工解译矿体剖面轮廓线导致的交互繁琐的麻烦，有效提高了矿体建模的效率和建模质量。
本实施例矿体建模方法，进一步通过点、线、面约束规则，可以构建具有特定约束规则的几何边界约束条件，可以按照地质规则约束要求对矿体隐式模型进行动态更新，构建更加符合地质趋势的矿体三维模型，改善建模质量。且能实现对逐步精细化建模的复杂矿体隐式模型进行动态更新，提高三维建模的效率和灵活性。
下面结合具体的应用实例对本发明实施例矿体建模方法的应用进行举例说明：
该应用实例根据矿体勘探建立的矿体钻孔数据如表1所示：
表1矿体隐式建模参数
请参阅图6，图6为本发明实施例中1#钻孔数据矿体建模过程示意图。具体地，首先根据边界品位、工业品位和最小可采厚度等信息对钻孔数据库进行样品组合，如图6a所示。选取样段采样间距2、非样段采样间距3，对钻孔组合样段、非样段按采样参数进行离散化采样，如图6b所示。自动建模结果如图6c所示。
通过隐式建模进行自动建模的结果有时并不能反映地质规则的约束，比如矿体模型的外推趋势等。为了使建模结果更加符合矿体的延伸趋势，可以通过添加约束线、趋势面约束对隐式模型进行动态调整，最终的矿体隐式建模结果可以很好地符合地质工程师的解译要求，如图6d所示。不同分辨率尺寸下矿体模型三维重构结果的参数如表1所示。
图7为2#钻孔数据在不同边界品位下所建立的矿体动态品位壳模型示意图，其中图7a为不同边界品位Cu＝2％和Cu＝3％时重构的矿体模型，图7b为不同边界品位Cu＝1％和Cu＝3％时重构的矿体模型，可以很好地反映不同经济指标下矿体模型的空间分布趋势，对隐式模型进行动态更新。
为了实现上述实施例建模方法，本发明实施例还提供一种基于钻孔数据的矿体建模装置，请参阅图8，该装置包括：
采样模块801，用于对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；
第一距离场构建模块802，用于根据所述接触点约束和所述离面点约束构建离散点距离场；
第二距离场构建模块803，用于对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；
建模模块804，用于基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
在一可选的实施方式中，采样模块801具体用于：按照设定的条件对所述钻孔数据进行预处理得到钻孔数据样品组合；在所述钻孔数据样品组合的基础上按钻孔轨迹对样段和非样段进行离散点采样。
在一可选的实施方式中，第一距离场构建模块802具体用于：校验并修正所述离面点约束；根据所述接触点约束和修正后的所述离面点约束构建所述离散点距离场。
在一可选的实施方式中，第二距离场构建模块803具体用于：根据所述离散点距离场和所述离散采样点构建插值约束，将所述插值约束代入插值方程进行插值计算得到与插值计算后的插值方程对应的所述隐式函数距离场。
在一可选的实施方式中，第二距离场构建模块803还用于：
添加附加插值趋势约束；
根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
在一可选的实施方式中，所述添加插值趋势约束，包括：
添加插值趋势约束点、插值趋势约束线和/或插值趋势面信息；
所述根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场，包括：
根据所述插值趋势约束点、所述插值趋势约束线和/或所述插值趋势面信息，对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
在一可选的实施方式中，建模模块804具体用于：根据插值约束提取用于界定模型边界的最小外包；根据所述最小外包确定模型三维重构的分辨率；按照所述分辨率对所述最小外包进行空间分割，构造离散化的空间规则数据场；对所述隐式函数距离场按照所述离散化的空间规则数据场进行采样计算，重构所述三维模型。
需要说明的是：上述实施例提供的矿体建模装置在进行三维建模时，仅以上述各程序模块的划分进行举例说明，实际应用中，可以根据需要而将上述处理分配由不同的程序模块完成，即将矿体建模装置的内部结构划分成不同的程序模块，以完成以上描述的全部或者部分处理。另外，上述实施例提供的矿体建模装置与矿体建模方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
在实际应用中，上述各程序模块均可由服务器上的中央处理器(CPU，CentralProcessing Unit)、微处理器(MPU，Micro Processor Unit)、数字信号处理器(DSP，Digital Signal Processor)、或现场可编程门阵列(FPGA，Field Programmable GateArray)等实现。
为了实现上述矿体建模方法，本发明实施例还提供了一种矿体建模***的硬件结构。现在将参考附图描述实现矿体建模***，该矿体建模***可以以各种类型的处理器，如客户端、服务器等形式来实施。下面对本发明实施例的矿体建模***的硬件结构做进一步说明，可以理解，图9仅仅示出了矿体建模***的示例性结构而非全部结构，根据需要可以实施如图9示出的部分结构或全部结构。
参见图9，图9为本发明实施例提供的一种矿体建模***的硬件结构示意图，实际应用中可以应用于前述运行应用程序的服务器，图9所示的矿体建模***900包括：至少一个处理器901、存储器902、用户接口903和至少一个网络接口904。该矿体建模***900中的各个组件通过总线***905耦合在一起。可以理解，总线***905用于实现这些组件之间的连接通信。总线***905除包括数据总线之外，还可以包括电源总线、控制总线和状态信号总线。但是为了清楚说明起见，在图9中将各种总线都标为总线***905。
其中，用户接口903可以包括显示器、键盘、鼠标、轨迹球、点击轮、按键、按钮、触感板或者触摸屏等。
可以理解，存储器902可以是易失性存储器或非易失性存储器，也可包括易失性非易失性存储器两者。
本发明实施例中的存储器902用于存储各种类型的数据以支持矿体建模***900的操作。这些数据的示例包括：用于在矿体建模***上操作的任何计算机程序，如可执行程序9021和操作***9022，实现本发明实施例的矿体建模方法的程序可以包含在可执行程序9021中。
本发明实施例揭示的矿体建模方法可以应用于处理器901中，或者由处理器901实现。处理器901可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述矿体建模方法的各步骤可以通过处理器901中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器901可以是通用处理器、DSP，或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。处理器901可以实现或者执行本发明实施例中提供的各矿体建模方法、步骤及逻辑框图。通用处理器可以是微处理器或者任何常规的处理器等。本发明实施例所提供的矿体建模方法的步骤，可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于存储介质中，该存储介质位于存储器902，处理器901读取存储器902中的信息，结合其硬件完成本发明实施例提供的矿体建模方法的步骤。
在本实施例中，该矿体建模***900包括存储器902、处理器901及存储在存储器902上并能够由所述处理器901运行的可执行程序9021，所述处理器901运行所述可执行程序9021时实现前述实施例的矿体建模方法。
在示例性实施例中，本发明实施例还提供了一种存储介质，该存储介质可为移动存储设备、只读存储器(Read Only Memory，ROM)、光盘、闪存或磁盘等存储介质，可选为非瞬间存储介质。存储介质上存储有可执行程序9021，该可执行程序9021被处理器运行时，实现前述实施例的矿体建模方法。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或可执行程序产品。因此，本发明可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的可执行程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和可执行程序产品的流程图和/或方框图来描述的。应理解可由可执行程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些可执行程序指令到通用计算机、专用计算机、嵌入式处理机或参考可编程数据处理设备的处理器以产生一个机器，使得通过计算机或参考可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些可执行程序指令也可存储在能引导计算机或参考可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些可执行程序指令也可装载到计算机或参考可编程数据处理设备上，使得在计算机或参考可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或参考可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种基于钻孔数据的矿体建模方法，其特征在于，包括：
对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；
根据所述接触点约束和所述离面点约束构建离散点距离场；
对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；
基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
2.如权利要求1所述的基于钻孔数据的矿体建模方法，其特征在于，所述对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，包括：
按照设定的条件对所述钻孔数据进行预处理得到钻孔数据样品组合；
在所述钻孔数据样品组合的基础上按钻孔轨迹对样段和非样段进行离散化采样得到离散采样点。
3.如权利要求1所述的基于钻孔数据的矿体建模方法，其特征在于，所述根据所述接触点约束和所述离面点约束构建离散点距离场，包括：
校验并修正所述离面点约束；
根据所述接触点约束和修正后的所述离面点约束构建所述离散点距离场。
4.如权利要求1所述的基于钻孔数据的矿体建模方法，其特征在于，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，包括：
根据所述离散点距离场和所述离散采样点构建插值约束，将所述插值约束代入插值方程进行插值计算得到与插值计算后的插值方程对应的所述隐式函数距离场。
5.如权利要求4所述的基于钻孔数据的矿体建模方法，其特征在于，所述对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场，还包括：
添加附加插值趋势约束；
根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
6.如权利要求5所述的基于钻孔数据的矿体建模方法，其特征在于，所述添加插值趋势约束，包括：
添加插值趋势约束点、插值趋势约束线和/或插值趋势面信息；
所述根据所述附加插值趋势约束对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场，包括：
根据所述插值趋势约束点、所述插值趋势约束线和/或所述插值趋势面信息，对所述插值计算后的插值方程进行插值计算，动态更新所述隐式函数距离场。
7.如权利要求4所述的基于钻孔数据的矿体建模方法，其特征在于，所述基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型，包括：
根据所述插值约束提取用于界定模型边界的最小外包；
根据所述最小外包确定模型三维重构的分辨率；
按照所述分辨率对所述最小外包进行空间分割，构造离散化的空间规则数据场；
对所述隐式函数距离场按照所述离散化的空间规则数据场进行采样计算，构建所述三维模型。
8.一种基于钻孔数据的矿体建模装置，其特征在于，包括：
采样模块，用于对地质勘探采样的钻孔数据进行离散化采样得到离散采样点，提取所述离散采样点中样段的采样点并转化为接触点约束，提取非样段的采样点并转化为离面点约束；
第一距离场构建模块，用于根据所述接触点约束和所述离面点约束构建离散点距离场；
第二距离场构建模块，用于对所述离散点距离场进行空间插值，构建用于表征矿体边界距离场分布的隐式函数距离场；
建模模块，用于基于参考等值面提取算法对所述隐式函数距离场进行采样计算，得到矿体的三维模型。
9.一种终端设备，包括存储器、处理器及存储在所述存储器上并能够由所述处理器运行的可执行程序，其特征在于，所述处理器运行所述可执行程序时执行如权利要求1至7任一项所述的基于钻孔数据的矿体建模方法。
10.一种存储介质，其上存储有可执行程序，其特征在于，所述可执行程序被处理器执行时实现如权利要求1至7任一项所述的基于钻孔数据的矿体建模方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421277.6A CN109326002B (zh) | 2018-11-27 | 2018-11-27 | 基于钻孔数据的矿体建模方法、装置、***及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811421277.6A CN109326002B (zh) | 2018-11-27 | 2018-11-27 | 基于钻孔数据的矿体建模方法、装置、***及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109326002A true CN109326002A (zh) | 2019-02-12 |
| CN109326002B CN109326002B (zh) | 2022-11-18 |
Family
ID=65257514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811421277.6A Active CN109326002B (zh) | 2018-11-27 | 2018-11-27 | 基于钻孔数据的矿体建模方法、装置、***及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109326002B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111192359A (zh) * | 2019-12-30 | 2020-05-22 | 核工业北京地质研究院 | 利用地质剖面钻孔空间轨迹形态提取钻孔岩性数据的方法 |
| CN111260783A (zh) * | 2020-01-19 | 2020-06-09 | 中国地质大学（武汉） | 一种基于K近邻和Poisson曲面的矿体三维自动建模方法 |
| CN111459955A (zh) * | 2020-03-13 | 2020-07-28 | 济南轨道交通集团有限公司 | 一种基于gis平台的三维地质结构模型自动更新方法及*** |
| CN111968231A (zh) * | 2020-08-14 | 2020-11-20 | 成都理工大学 | 一种基于地质图切剖面的三维地层建模方法 |
| CN112036030A (zh) * | 2020-08-31 | 2020-12-04 | 中南大学 | 矿体组合约束建模方法、装置、设备及存储介质 |
| CN112115613A (zh) * | 2020-09-21 | 2020-12-22 | 中南大学 | 矿体隐式模型协同构建方法、装置、设备及存储介质 |
| CN112907743A (zh) * | 2021-03-08 | 2021-06-04 | 中南大学 | 一种基于钻孔数据的地质体剖/平面自动成图方法、装置、设备及介质 |
| CN112927350A (zh) * | 2021-02-23 | 2021-06-08 | 中南大学 | 多域地质体模型构建方法、装置、设备及存储介质 |
| CN115163025A (zh) * | 2021-12-14 | 2022-10-11 | 核工业北京化工冶金研究院 | 一种砂岩型铀矿原地浸出精准开采方法 |
| CN118037980A (zh) * | 2024-04-11 | 2024-05-14 | 中国地质大学（武汉） | 基于地质空间插值的地层确定方法及装置、介质、设备 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130262052A1 (en) * | 2012-04-03 | 2013-10-03 | Jean-Laurent Mallet | System and method for generating an implicit model of geological horizons |
-
2018
- 2018-11-27 CN CN201811421277.6A patent/CN109326002B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130262052A1 (en) * | 2012-04-03 | 2013-10-03 | Jean-Laurent Mallet | System and method for generating an implicit model of geological horizons |
Non-Patent Citations (3)
| Title |
|---|
| 刘亚静，米雪玉，陈光等著: "《非层状矿体空间构模与数据存储关键技术研究》", 31 August 2014, 武汉：武汉大学出版社 * |
| 杨云保，唐永虎，徐惠长主编: "《固体矿产勘查技术》", 31 July 2007, 北京：地质出版社 * |
| 郭甲腾等: "基于径向基函数曲面的矿体隐式自动三维建模方法", 《煤炭学报》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111192359B (zh) * | 2019-12-30 | 2023-10-20 | 核工业北京地质研究院 | 利用地质剖面钻孔空间轨迹形态提取钻孔岩性数据的方法 |
| CN111192359A (zh) * | 2019-12-30 | 2020-05-22 | 核工业北京地质研究院 | 利用地质剖面钻孔空间轨迹形态提取钻孔岩性数据的方法 |
| CN111260783B (zh) * | 2020-01-19 | 2022-06-24 | 中国地质大学（武汉） | 一种基于K近邻和Poisson曲面的矿体三维自动建模方法 |
| CN111260783A (zh) * | 2020-01-19 | 2020-06-09 | 中国地质大学（武汉） | 一种基于K近邻和Poisson曲面的矿体三维自动建模方法 |
| CN111459955A (zh) * | 2020-03-13 | 2020-07-28 | 济南轨道交通集团有限公司 | 一种基于gis平台的三维地质结构模型自动更新方法及*** |
| CN111459955B (zh) * | 2020-03-13 | 2023-09-29 | 济南轨道交通集团有限公司 | 一种基于gis平台的三维地质结构模型自动更新方法及*** |
| CN111968231A (zh) * | 2020-08-14 | 2020-11-20 | 成都理工大学 | 一种基于地质图切剖面的三维地层建模方法 |
| CN111968231B (zh) * | 2020-08-14 | 2023-05-30 | 成都理工大学 | 一种基于地质图切剖面的三维地层建模方法 |
| CN112036030A (zh) * | 2020-08-31 | 2020-12-04 | 中南大学 | 矿体组合约束建模方法、装置、设备及存储介质 |
| CN112036030B (zh) * | 2020-08-31 | 2021-06-01 | 中南大学 | 矿体组合约束建模方法、装置、设备及存储介质 |
| CN112115613A (zh) * | 2020-09-21 | 2020-12-22 | 中南大学 | 矿体隐式模型协同构建方法、装置、设备及存储介质 |
| CN112927350A (zh) * | 2021-02-23 | 2021-06-08 | 中南大学 | 多域地质体模型构建方法、装置、设备及存储介质 |
| CN112927350B (zh) * | 2021-02-23 | 2021-11-02 | 中南大学 | 多域地质体模型构建方法、装置、设备及存储介质 |
| CN112907743B (zh) * | 2021-03-08 | 2022-05-06 | 中南大学 | 一种地质体剖/平面自动成图方法、装置、设备及介质 |
| CN112907743A (zh) * | 2021-03-08 | 2021-06-04 | 中南大学 | 一种基于钻孔数据的地质体剖/平面自动成图方法、装置、设备及介质 |
| CN115163025A (zh) * | 2021-12-14 | 2022-10-11 | 核工业北京化工冶金研究院 | 一种砂岩型铀矿原地浸出精准开采方法 |
| CN115163025B (zh) * | 2021-12-14 | 2023-09-08 | 核工业北京化工冶金研究院 | 一种砂岩型铀矿原地浸出精准开采方法 |
| CN118037980A (zh) * | 2024-04-11 | 2024-05-14 | 中国地质大学（武汉） | 基于地质空间插值的地层确定方法及装置、介质、设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109326002B (zh) | 2022-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109326002A (zh) | 基于钻孔数据的矿体建模方法、装置、***及存储介质 | |
| CN109584357B (zh) | 基于多轮廓线的三维建模方法、装置、***及存储介质 | |
| Zhihong et al. | Evaluation of granular particle roundness using digital image processing and computational geometry | |
| CN110033519A (zh) | 基于隐式函数的三维建模方法、装置、***及存储介质 | |
| JP2008527473A (ja) | ３次元モデルの検索方法、検索装置及び検索プログラム | |
| CN105654483A (zh) | 三维点云全自动配准方法 | |
| Liu et al. | Parameterization for molecular Gaussian surface and a comparison study of surface mesh generation | |
| Remelli et al. | Neuralsampler: Euclidean point cloud auto-encoder and sampler | |
| Qin et al. | A comparison of high-level design tools for soc-fpga on disparity map calculation example | |
| US9582626B1 (en) | Using waveform propagation for accurate delay calculation | |
| Xu et al. | Developing an extended IFC data schema and mesh generation framework for finite element modeling | |
| JP2020115338A (ja) | 生メッシュからの構造化３ｄモデルの生成 | |
| CN115984503A (zh) | 一种地质剖面图生成方法、***、电子设备及介质 | |
| Clyne et al. | Physically-based feature tracking for CFD data | |
| JP7056452B2 (ja) | シミュレーション装置、シミュレーション方法およびシミュレーションプログラム | |
| He et al. | Reverse Engineering of Free‐Form Surface Based on the Closed‐Loop Theory | |
| Emiris et al. | Exact Voronoi diagram of smooth convex pseudo-circles: General predicates, and implementation for ellipses | |
| US9745842B2 (en) | Proximity calculation in a geoscience domain | |
| Fang et al. | Volumetric data modeling and analysis based on seven-directional box spline | |
| Krone et al. | Parallel computation and interactive visualization of time-varying solvent excluded surfaces | |
| CN112927350B (zh) | 多域地质体模型构建方法、装置、设备及存储介质 | |
| Wang et al. | Parametric Primitive Analysis of CAD Sketches With Vision Transformer | |
| Zhao et al. | Improved DDA method based on explicitly solving contact constraints | |
| US9754180B2 (en) | Robust automatic computation of ridges and valleys of a height field | |
| Lizier et al. | Comparing techniques for tetrahedral mesh generation |
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 |