CN112837413B - 面向地质钻孔的虚拟地层推断方法及装置 - Google Patents

面向地质钻孔的虚拟地层推断方法及装置 Download PDF

Info

Publication number
CN112837413B
CN112837413B CN202110180713.0A CN202110180713A CN112837413B CN 112837413 B CN112837413 B CN 112837413B CN 202110180713 A CN202110180713 A CN 202110180713A CN 112837413 B CN112837413 B CN 112837413B
Authority
CN
China
Prior art keywords
drilling
stratum
borehole
virtual
inferred
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
CN202110180713.0A
Other languages
English (en)
Other versions
CN112837413A (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.)
Nanjing Normal University
Original Assignee
Nanjing Normal 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 Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN202110180713.0A priority Critical patent/CN112837413B/zh
Publication of CN112837413A publication Critical patent/CN112837413A/zh
Application granted granted Critical
Publication of CN112837413B publication Critical patent/CN112837413B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00—Computing arrangements using knowledge-based models
    • G06N5/04—Inference or reasoning models
    • G06N5/041—Abduction
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Artificial Intelligence (AREA)
  • Remote Sensing (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种面向地质钻孔的虚拟地层推断方法及装置，方法具体包括步骤：(1)加载钻孔基本信息及钻孔地层信息，计算已知地层厚度；(2)基于约束边界数据，划分钻孔类型，并确定需推断地层编号区间C；(3)读取地层编号c，构建推断钻孔集合ND与非推断钻孔集合ID；(4)基于钻孔集合ID，构建尖灭钻孔集合AD，二者合并创建三角网M；(5)基于三角网M，推断钻孔集合ND中每钻孔缺失的虚拟地层，反距离加权插值计算地层厚度。本发明实现了一种快速推断钻孔地层缺失信息方法。本发明与现有技术相比，既提高了地质钻孔地层完整程度，又具有较高的自动化程度。

Description

面向地质钻孔的虚拟地层推断方法及装置
技术领域
本发明涉及地理信息和地质学领域，具体涉及一种面向地质钻孔的虚拟地层推断方法及装置。
背景技术
地质钻孔是获取地质信息的主要勘察手段，也是进行三维建模的重要数据来源。根据钻孔数据构建的三维地质模型，能够准确地表达复杂地质现象的边界条件、地质体内发育的各种地质构造，逼真地呈现三维动态显示效果。
然而，由于钻探成本昂贵、钻探目标不同等原因，同一研究区域的不同钻孔，钻探深度和钻遇地层有较大不同。部分较浅钻孔的存在，不可避免地会影响基于相关钻孔所构建三维模型的质量。为了弱化这一问题，针对较浅钻孔，基于地质规律和地质知识进行相关虚拟地层的自动推断，具有重要的研究意义和应用价值。
发明内容
发明目的：本发明针对现有技术存在的问题，提供一种面向地质钻孔的虚拟地层推断方法及装置。
技术方案：本发明所述的面向地质钻孔的虚拟地层推断方法包括：
(1)读取钻孔数据、钻孔地层信息和约束边界矢量数据，生成钻孔集合D、地层集合S、约束边界集合Boundary和地层厚度集合H；
(2)基于约束边界集合Boundary，将钻孔集合D划分为边界钻孔集合BD和非边界钻孔集合FBD，并根据地层集合S确定需推断地层的地层编号区间C；
(3)从地层编号区间C中获取任一地层编号c，构建推断钻孔集合ND与非推断钻孔集合ID；
(4)基于钻孔集合ID，构建尖灭钻孔集合AD，合并集合ID和集合AD为钻孔集合UD，并根据集合UD创建三角网M；
(5)基于三角网M、边界钻孔集合BD和非边界钻孔集合FBD，获取推断钻孔集合ND中所有钻孔缺失的虚拟底层，并采用反距离加权插值计算虚拟地层的厚度；
(6)循环执行步骤(3)-(5)，直至地层编号区间C被遍历完，完成所有钻孔虚拟地层的推断，更新地层信息，保存含有虚拟地层的钻孔数据。
进一步的，步骤(1)具体包括：
(1-1)读取钻孔数据到钻孔集合D＝{di|i＝1,2,…,DN}；其中，i表示钻孔序号，DN表示钻孔数量；
(1-2)读取所有钻孔的地层信息，存入地层集合S＝{sij|i＝1,2,…,DN,j＝1,2,…,SNi}，其中，sij表示钻孔di的第j个地层，SNi表示钻孔di的地层数量；
(1-3)读取约束边界矢量数据到约束边界集合Boundary；
(1-4)计算每一地层的厚度，存入地层厚度集合H＝{sij.h|i＝1,2,…，DN，j＝1，2,…,SNi},其中，sij.h＝sij.et-sij.eb表示sij的地层厚度，sij.et为sij的层顶深度，sij.eb为sij的层底深度。
进一步的，步骤(2)具体包括：
(2-1)获取钻孔集合D内每一钻孔的坐标，判断钻孔是否位于约束边界集合Boundary中的约束边界上，若是则标记钻孔为边界钻孔，并存入边界钻孔集合BD中，否则标记为非边界钻孔，并存入非边界钻孔集合FBD中；
(2-2)从地层集合S中获取钻孔集合D内所有钻孔的末尾地层的地层编号，得到地层编号的最小值α与最大值β，构建需推断地层的地层编号区间C＝[α+1，β],区间内取正整数。
进一步的，步骤(3)具体包括：
(3-1)从地层编号区间C中获取任一地层编号c；
(3-2)从地层集合S中获取钻孔集合D内每一钻孔的末尾地层的地层编号，若该地层编号小于c，则代表该地层编号对应的钻孔为需推断钻孔，将其保存到需推断钻孔集合ND，否则代表该地层编号对应的钻孔为非推断钻孔，将其保存到非推断钻孔集合ID中；
(3-3)若非推断钻孔集合ID中的非推断钻孔数量LN≥3，则返回执行步骤(3-1)，否则执行步骤(4)。
进一步的，步骤(4)具体包括：
(4-1)根据非推断钻孔集合ID，采用1/2尖灭规则构建尖灭钻孔集合AD；
(4-2)合并非推断钻孔集合ID与尖灭钻孔集合AD为钻孔集合UD；
(4-3)根据钻孔集合UD内的钻孔坐标信息，创建Delaunary三角网M＝{mt|t＝1,2,…，TM}，其中，mt表示三角网内第t个三角形，TM表示三角网内三角形数量。
进一步的，步骤(4-1)具体包括：
(4-1-1)根据非推断钻孔集合ID中每一钻孔坐标信息，构建Delaunary三角网T0，并获取T0内的所有三角边，存入三角边集合E；
(4-1-2)从三角边集合E中读取任一三角边e，获取位于三角边e两端的钻孔dl、dl+1，在两钻孔之间的1/2处创建一尖灭钻孔do，若钻孔dl、dl+1的地层厚度slc.h、s(l+1)c.h均不为0，则将尖灭钻孔do的地层厚度soc.h赋值为(slc.h+s(l+1)c.h)/2，否则地层soc.h＝0，
(4-1-3)重复执行步骤(4-1-2)，直至实现集合E中所有三角边的遍历，得到所有尖灭钻孔，并将所有尖灭钻孔存入尖灭钻孔集合AD中。
进一步的，步骤(5)具体包括：
(5-1)从推断钻孔集合ND中获取任一推断钻孔dk，判断推断钻孔dk是否存在虚拟地层；
(5-2)若存在虚拟地层，则根据钻孔集合UD内每个钻孔的地层厚度，采用反距离加权插值法计算虚拟地层的厚度；
(5-3)循环执行步骤(5-1)-(5-2)，直至完成推断钻孔集合ND中所有推断钻孔的遍历，得到所有推断钻孔的虚拟地层推断。
进一步的，步骤(5-1)具体包括：
(5-1-1)从推断钻孔集合ND中获取任一推断钻孔dk，根据钻孔dk与三角网M内三角形的位置推断是否存在虚拟地层，具体方法如下：(a)读取三角网M内任一三角形mt，根据三角形mt三个顶点与钻孔dk点坐标采用同向法判断钻孔dk是否位于三角形mt内；(b)若钻孔dk位于三角形mt内，则获取mt三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示推断钻孔dk存在虚拟地层，执行步骤(5-2)，反之表示推断钻孔dk不存在虚拟地层，则设置推断钻孔dk对应的虚拟地层厚度为0；若钻孔dk不位于三角形mt内，则返回执行步骤(a)；(c)若三角网M被遍历完，钻孔dk不位于所有三角形mt内，则执行步骤(5-1-2)；
(5-1-2)根据推断钻孔dk类型推断虚拟地层，具体方法如下：(a)若dk位于非边界钻孔集合FBD内，则表示dk不存在虚拟地层，则设置dk对应的虚拟地层厚度为0；(b)若dk位于边界钻孔集合BD内，则计算三角网M内每一三角形的中心与dk点的距离，根据距离大小获取离dk最近的三角形mp；(c)以钻孔dk点坐标为中心，dk与三角形mp中心的距离为半径构建几何圆；若几何圆内存在非边界钻孔，则表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0；若几何圆内不存在非边界钻孔，则获取mp三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示dk存在虚拟地层，执行步骤(5-2)，反之表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0。
本发明所述的面向地质钻孔的虚拟地层推断装置，包括处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述方法。
有益效果：本发明与现有技术相比，其显著优点是：本发明通过数据加载、虚拟地层推断区间确定，尖灭钻孔集合构建、虚拟地层推断与厚度计算等环节，提供了一种面向地质钻孔虚拟地层的推断方法，其既提高了地质钻孔地层完整程度，又具有较高的自动化程度，采用本方法可以使得钻孔地层信息得到完善，以实现高效精确的地质体三维建模。
附图说明
图1是本实施例中采用的钻孔基本信息(部分)；
图2是本实施例中采用的地层基本信息(部分)；
图3是本发明提供的面向地质钻孔的虚拟地层推断方法的流程图；
图4是钻孔分布与约束边界示意图；
图5是钻孔类型分布图；
图6是读取地层编号13划分的推断钻孔与非推断钻孔分布图；
图7是Delaunary三角网T0示意图；
图8是尖灭钻孔分布图；
图9是Delaunary三角网M示意图；
图10是一工程地质钻孔虚拟地层推断前后对比表(a为推断前，b为推断后)；
图11是一工程地质钻孔虚拟地层推断前后三维实体模型示意图(a为推断前全局图，b为推断前局部图，c为推断后全局图，d为推断后局部图)。
具体实施方式
下面对本发明技术方案作进一步详细的说明，本实施例的实验数据采用的是某研究区域工程地质钻孔数据(图1)与地层数据(图2)，该实验数据采用的投影坐标系为WGS_1984_UTM_Zone_50N。下面结合附图，并通过描述一个具体的实施例，来进一步说明。
如图3所示，本发明提供的面向地质钻孔的虚拟地层推断方法包括：
(1)读取钻孔数据、钻孔地层信息和约束边界矢量数据，生成钻孔集合D、地层集合S、约束边界集合Boundary和地层厚度集合H。在本实施例中，钻孔分布与约束边界示意图如图4所示。
该步骤具体包括：
(1-1)读取钻孔数据到钻孔集合D＝{di|i＝1,2,…,DN}；其中，i表示钻孔序号，DN表示钻孔数量；在本实施例中，DN＝69；
(1-2)读取所有钻孔的地层信息，存入地层集合S＝{sij|i＝1,2,…,DN,j＝1,2,…,SNi}，其中，sij表示钻孔di的第j个地层，SNi表示钻孔di的地层数量；
(1-3)读取约束边界矢量数据到约束边界集合Boundary；
(1-4)计算每一地层的厚度，存入地层厚度集合H＝{sij.h|i＝1,2,…,DN,j＝1，2，…,SNi},其中，sij.h＝sij.et-sij.eb表示sij的地层厚度，sij.et为sij的层顶深度，sij.eb为sij的层底深度。本实施例中，取i＝1，j＝1时，s11.et＝0，sij.eb＝16.74，sij.h＝16.74。
(2)基于约束边界集合Boundary，将钻孔集合D划分为边界钻孔集合BD和非边界钻孔集合FBD，并根据地层集合S确定需推断地层的地层编号区间C。
该步骤具体包括：
(2-1)获取钻孔集合D内每一钻孔的坐标，判断钻孔是否位于约束边界集合Boundary中的约束边界上，若是则标记钻孔为边界钻孔，并存入边界钻孔集合BD中，否则标记为非边界钻孔，并存入非边界钻孔集合FBD中；在本实施例中，钻孔类型分布图如图5所示；
(2-2)从地层集合S中获取钻孔集合D内所有钻孔的末尾地层的地层编号，得到地层编号的最小值α与最大值β，构建需推断地层的地层编号区间C＝[α+1，β],区间内取正整数。在本实施例中，α＝12，β＝35。
(3)从地层编号区间C中获取任一地层编号c，构建推断钻孔集合ND与非推断钻孔集合ID。
该步骤具体包括：
(3-1)从地层编号区间C中获取任一地层编号c；
(3-2)从地层集合S中获取钻孔集合D内每一钻孔的末尾地层的地层编号，若该地层编号小于c，则代表该地层编号对应的钻孔为需推断钻孔，将其保存到需推断钻孔集合ND，否则代表该地层编号对应的钻孔为非推断钻孔，将其保存到非推断钻孔集合ID中；
(3-3)若非推断钻孔集合ID中的非推断钻孔数量LN≥3，则返回执行步骤(3-1)，否则执行步骤(4)。在本实施例中，取c＝13时，LN＝68需推断钻孔与非推断钻孔分布图如图6所示。
(4)基于钻孔集合ID，构建尖灭钻孔集合AD，合并集合ID和集合AD为钻孔集合UD，并根据集合UD创建三角网M。
该步骤具体包括：
(4-1)根据非推断钻孔集合ID，采用1/2尖灭规则构建尖灭钻孔集合AD；构建方法包括：
(4-1-1)根据非推断钻孔集合ID中每一钻孔坐标信息，构建Delaunary三角网T0，并获取T0内的所有三角边，存入三角边集合E；在本实施例中，三角网T0如图7所示；
(4-1-2)从三角边集合E中读取任一三角边e，获取位于三角边e两端的钻孔dl、dl+1，在两钻孔之间的1/2处创建一尖灭钻孔do，若钻孔dl、dl+1的地层厚度slc.h、s(l+1)c.h均不为0，则将尖灭钻孔do的地层厚度soc.h赋值为(slc.h+s(l+1)c.h)/2，否则地层soc.h＝0，
(4-1-3)重复执行步骤(4-1-2)，直至实现集合E中所有三角边的遍历，得到所有尖灭钻孔，并将所有尖灭钻孔存入尖灭钻孔集合AD中；尖灭钻孔分布图如图8所示；
(4-2)合并非推断钻孔集合ID与尖灭钻孔集合AD为钻孔集合UD；
(4-3)根据钻孔集合UD内的钻孔坐标信息，创建Delaunary三角网M＝{mt|t＝1,2,…,TM}，其中，mt表示三角网内第t个三角形，TM表示三角网内三角形数量。三角网M如图9所示。
(5)基于三角网M、边界钻孔集合BD和非边界钻孔集合FBD，获取推断钻孔集合ND中所有钻孔缺失的虚拟底层，并采用反距离加权插值计算虚拟地层的厚度。
该步骤具体包括：
(5-1)从推断钻孔集合ND中获取任一推断钻孔dk，判断推断钻孔dk是否存在虚拟地层；具体判断方法为：
(5-1-1)根据钻孔dk与三角网M内三角形的位置推断是否存在虚拟地层，具体方法如下：(a)读取三角网M内任一三角形mt，根据三角形mt三个顶点与钻孔dk点坐标采用同向法判断钻孔dk是否位于三角形mt内；(b)若钻孔dk位于三角形mt内，则获取mt三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示推断钻孔dk存在虚拟地层，执行步骤(5-2)，反之表示推断钻孔dk不存在虚拟地层，则设置推断钻孔dk对应的虚拟地层厚度为0；若钻孔dk不位于三角形mt内，则返回执行步骤(a)；(c)若三角网M被遍历完，钻孔dk不位于所有三角形mt内，则执行步骤(5-1-2)；(所述同向法，参考在下列文件:Philip J.S.,David H.E.,Geometric Tools for Computer Graphics,Elsevier Science,2003,p64-65)；
(5-1-2)根据推断钻孔dk类型推断虚拟地层，具体方法如下：(a)若dk位于非边界钻孔集合FBD内，则表示dk不存在虚拟地层，则设置dk对应的虚拟地层厚度为0；(b)若dk位于边界钻孔集合BD内，则计算三角网M内每一三角形的中心与dk点的距离，根据距离大小获取离dk最近的三角形mp；(c)以钻孔dk点坐标为中心，dk与三角形mp中心的距离为半径构建几何圆；若几何圆内存在非边界钻孔，则表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0；若几何圆内不存在非边界钻孔，则获取mp三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示dk存在虚拟地层，执行步骤(5-2)，反之表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0；
(5-2)若存在虚拟地层，则根据钻孔集合UD内每个钻孔的地层厚度，采用反距离加权插值法计算虚拟地层的厚度；
所述反距离加权插值法，参考于下列文件：[1]李安波,周良辰,闾国年.2013.地质信息系统.北京:科学出版社.P63.
(5-3)循环执行步骤(5-1)-(5-2)，直至完成推断钻孔集合ND中所有推断钻孔的遍历，得到所有推断钻孔的虚拟地层推断。
(6)循环执行步骤(3)-(5)，直至地层编号区间C被遍历完，完成所有钻孔虚拟地层的推断，更新地层信息，保存含有虚拟地层的钻孔数据。
本发明提供的面向地质钻孔的虚拟地层推断装置包括处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述方法。
在本实施例中，研究区域内一地质工程钻孔的虚拟地层推断前后对比表如图10所示。基于推断前后钻孔构建该研究区域三维实体模型，模型侧视图如图11所示，相比于推断前三维模型，推断后的实体模型地层连续，符合地层发育的特点，建模效果更加合理。本发明实施例中基于GDAL/OGR提供部分GIS操作，相关步骤也可以使用SuperMap、ArcgisObject等软件的API进行相应GIS操作。
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (7)

1.一种面向地质钻孔的虚拟地层推断方法，其特征在于该方法包括：
(1)读取钻孔数据、钻孔地层信息和约束边界矢量数据，生成钻孔集合D、地层集合S、约束边界集合Boundary和地层厚度集合H；
(2)基于约束边界集合Boundary，将钻孔集合D划分为边界钻孔集合BD和非边界钻孔集合FBD，并根据地层集合S确定需推断地层的地层编号区间C；
(3)从地层编号区间C中获取任一地层编号c，构建推断钻孔集合ND与非推断钻孔集合ID；
(4)基于钻孔集合ID，构建尖灭钻孔集合AD，合并集合ID和集合AD为钻孔集合UD，并根据集合UD创建三角网M；
(5)基于三角网M、边界钻孔集合BD和非边界钻孔集合FBD，获取推断钻孔集合ND中所有钻孔缺失的虚拟底层，并采用反距离加权插值计算虚拟地层的厚度；
(6)循环执行步骤(3)-(5)，直至地层编号区间C被遍历完，完成所有钻孔虚拟地层的推断，更新地层信息，保存含有虚拟地层的钻孔数据；
其中，步骤(5)具体包括：
(5-1)从推断钻孔集合ND中获取任一推断钻孔dk，判断推断钻孔dk是否存在虚拟地层；
(5-2)若存在虚拟地层，则根据钻孔集合UD内每个钻孔的地层厚度，采用反距离加权插值法计算虚拟地层的厚度；
(5-3)循环执行步骤(5-1)-(5-2)，直至完成推断钻孔集合ND中所有推断钻孔的遍历，得到所有推断钻孔的虚拟地层推断；
步骤(5-1)具体包括：
(5-1-1)从推断钻孔集合ND中获取任一推断钻孔dk，根据钻孔dk与三角网M内三角形的位置推断是否存在虚拟地层，具体方法如下：(a)读取三角网M内任一三角形mt，根据三角形mt三个顶点与钻孔dk点坐标采用同向法判断钻孔dk是否位于三角形mt内；(b)若钻孔dk位于三角形mt内，则获取mt三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示推断钻孔dk存在虚拟地层，执行步骤(5-2)，反之表示推断钻孔dk不存在虚拟地层，则设置推断钻孔dk对应的虚拟地层厚度为0；若钻孔dk不位于三角形mt内，则返回执行步骤(a)；(c)若三角网M被遍历完，钻孔dk不位于所有三角形mt内，则执行步骤(5-1-2)；
(5-1-2)根据推断钻孔dk类型推断虚拟地层，具体方法如下：(a)若dk位于非边界钻孔集合FBD内，则表示dk不存在虚拟地层，则设置dk对应的虚拟地层厚度为0；(b)若dk位于边界钻孔集合BD内，则计算三角网M内每一三角形的中心与dk点的距离，根据距离大小获取离dk最近的三角形mp；(c)以钻孔dk点坐标为中心，dk与三角形mp中心的距离为半径构建几何圆；若几何圆内存在非边界钻孔，则表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0；若几何圆内不存在非边界钻孔，则获取mp三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示dk存在虚拟地层，执行步骤(5-2)，反之表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0。
2.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(1)具体包括：
(1-1)读取钻孔数据到钻孔集合D＝{di|i＝1,2,…,DN}；其中，i表示钻孔序号，DN表示钻孔数量；
(1-2)读取所有钻孔的地层信息，存入地层集合S＝{sij|i＝1,2,…,DN,j＝1,2,…,SNi}，其中，sij表示钻孔di的第j个地层，SNi表示钻孔di的地层数量；
(1-3)读取约束边界矢量数据到约束边界集合Boundary；
(1-4)计算每一地层的厚度，存入地层厚度集合H＝{sij.h|i＝1,2,…,DN,j＝1,2,…,SNi},其中，sij.h＝sij.et-sij.eb表示sij的地层厚度，sij.et为sij的层顶深度，sij.eb为sij的层底深度。
3.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(2)具体包括：
(2-1)获取钻孔集合D内每一钻孔的坐标，判断钻孔是否位于约束边界集合Boundary中的约束边界上，若是则标记钻孔为边界钻孔，并存入边界钻孔集合BD中，否则标记为非边界钻孔，并存入非边界钻孔集合FBD中；
(2-2)从地层集合S中获取钻孔集合D内所有钻孔的末尾地层的地层编号，得到地层编号的最小值α与最大值β，构建需推断地层的地层编号区间C＝[α+1，β],区间内取正整数。
4.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(3)具体包括：
(3-1)从地层编号区间C中获取任一地层编号c；
(3-2)从地层集合S中获取钻孔集合D内每一钻孔的末尾地层的地层编号，若该地层编号小于c，则代表该地层编号对应的钻孔为需推断钻孔，将其保存到需推断钻孔集合ND，否则代表该地层编号对应的钻孔为非推断钻孔，将其保存到非推断钻孔集合ID中；
(3-3)若非推断钻孔集合ID中的非推断钻孔数量LN≥3，则返回执行步骤(3-1)，否则执行步骤(4)。
5.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(4)具体包括：
(4-1)根据非推断钻孔集合ID，采用1/2尖灭规则构建尖灭钻孔集合AD；
(4-2)合并非推断钻孔集合ID与尖灭钻孔集合AD为钻孔集合UD；
(4-3)根据钻孔集合UD内的钻孔坐标信息，创建Delaunary三角网M＝{mt|t＝1,2,…,TM}，其中，mt表示三角网内第t个三角形，TM表示三角网内三角形数量。
6.根据权利要求5所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(4-1)具体包括：
(4-1-1)根据非推断钻孔集合ID中每一钻孔坐标信息，构建Delaunary三角网T0，并获取T0内的所有三角边，存入三角边集合E；
(4-1-2)从三角边集合E中读取任一三角边e，获取位于三角边e两端的钻孔dl、dl+1，在两钻孔之间的1/2处创建一尖灭钻孔do，若钻孔dl、dl+1的地层厚度slc.h、s(l+1)c.h均不为0，则将尖灭钻孔do的地层厚度soc.h赋值为(slc.h+s(l+1)c.h)/2，否则soc.h＝0，
(4-1-3)重复执行步骤(4-1-2)，直至实现集合E中所有三角边的遍历，得到所有尖灭钻孔，并将所有尖灭钻孔存入尖灭钻孔集合AD中。
7.一种面向地质钻孔的虚拟地层推断装置，包括处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现权利要求1-6中任意一项所述的方法。
CN202110180713.0A 2021-02-08 2021-02-08 面向地质钻孔的虚拟地层推断方法及装置 Active CN112837413B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110180713.0A CN112837413B (zh) 2021-02-08 2021-02-08 面向地质钻孔的虚拟地层推断方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110180713.0A CN112837413B (zh) 2021-02-08 2021-02-08 面向地质钻孔的虚拟地层推断方法及装置

Publications (2)

Publication Number Publication Date
CN112837413A CN112837413A (zh) 2021-05-25
CN112837413B true CN112837413B (zh) 2024-05-07

Family

ID=75933274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110180713.0A Active CN112837413B (zh) 2021-02-08 2021-02-08 面向地质钻孔的虚拟地层推断方法及装置

Country Status (1)

Country Link
CN (1) CN112837413B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113240771B (zh) * 2021-05-27 2021-12-07 广东省地震局 一种地层等厚线图的自动生成方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003330A (zh) * 2018-07-02 2018-12-14 南京师范大学 一种基于基岩边界约束的三维地层建模方法
CN111429573A (zh) * 2020-03-06 2020-07-17 南京师范大学 一种三维地质钻孔模型的自动构建方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2823710C (en) * 2011-01-27 2019-08-13 Landmark Graphics Corporation Methods and systems regarding models of underground formations

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109003330A (zh) * 2018-07-02 2018-12-14 南京师范大学 一种基于基岩边界约束的三维地层建模方法
CN111429573A (zh) * 2020-03-06 2020-07-17 南京师范大学 一种三维地质钻孔模型的自动构建方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
引入虚拟钻孔的地层连接推理与地质体建模;田宜平;吴冲龙;董志;;地质科技情报;20100915(05);全文 *
田宜平 ; 吴冲龙 ; 董志 ; .引入虚拟钻孔的地层连接推理与地质体建模.地质科技情报.2010,(05),全文. *

Also Published As

Publication number Publication date
CN112837413A (zh) 2021-05-25

Similar Documents

Publication Publication Date Title
Guo et al. Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm
CN112381937B (zh) 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法
Wu et al. An approach to computer modeling and visualization of geological faults in 3D
CN106844871A (zh) 基于bim的三维岩土工程勘察信息模型构建方法
CN106981093B (zh) 一种分区约束耦合的三维地层并行建模方法
AU2013402201B2 (en) In-situ wellbore, core and cuttings information system
US10853533B2 (en) Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof
CN109003330B (zh) 一种基于基岩边界约束的三维地层建模方法
CA3065179A1 (en) Constructing structural models of the subsurface
CN112150582B (zh) 一种面向多模态数据的地质剖面图近似表达方法
CN111459955A (zh) 一种基于gis平台的三维地质结构模型自动更新方法及系统
US10113421B2 (en) Three-dimensional fracture abundance evaluation of subsurface formations
CN103069460B (zh) 用于呈现钻井记录值的方法和系统
Jia et al. A triangulated irregular network constrained ordinary Kriging method for three-dimensional modeling of faulted geological surfaces
CN112837413B (zh) 面向地质钻孔的虚拟地层推断方法及装置
CN117197377A (zh) 矿山开采区三维动态可视化建模方法
CN110349262B (zh) 一种基于bim的隧道超前地质预报三维建模方法及系统
CN110197013B (zh) 基于Morphing的河床基岩面建模方法
CN105719342B (zh) 一种地裂缝地质体的三维建模可视化方法及装置
Hegemann et al. Hybrid ground data model for interacting simulations in mechanized tunneling
CN111862326A (zh) 地质模型的存储方法和计算方法、电子设备及存储介质
Carrell Tools and techniques for 3D geologic mapping in Arc Scene: Boreholes, cross sections, and block diagrams
CN113160399B (zh) 构建三维地质模型的方法、系统及可读存储介质
Gao et al. Mathematical Interpolation and Correction of Three-Dimensional Modelling of High-Speed Railway.
CN112581558A (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