CN112837413A - 面向地质钻孔的虚拟地层推断方法及装置 - Google Patents
面向地质钻孔的虚拟地层推断方法及装置 Download PDFInfo
- Publication number
- CN112837413A CN112837413A CN202110180713.0A CN202110180713A CN112837413A CN 112837413 A CN112837413 A CN 112837413A CN 202110180713 A CN202110180713 A CN 202110180713A CN 112837413 A CN112837413 A CN 112837413A
- Authority
- CN
- China
- Prior art keywords
- stratum
- borehole
- drilling
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005553 drilling Methods 0.000 claims abstract description 97
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 40
- 238000004590 computer program Methods 0.000 claims description 5
- 238000005755 formation reaction Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 5
- 239000007787 solid Substances 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
- 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 (9)
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被遍历完,完成所有钻孔虚拟地层的推断,更新地层信息,保存含有虚拟地层的钻孔数据。
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所述的面向地质钻孔的虚拟地层推断方法,其特征在于:步骤(5)具体包括:
(5-1)从推断钻孔集合ND中获取任一推断钻孔dk,判断推断钻孔dk是否存在虚拟地层;
(5-2)若存在虚拟地层,则根据钻孔集合UD内每个钻孔的地层厚度,采用反距离加权插值法计算虚拟地层的厚度;
(5-3)循环执行步骤(5-1)-(5-2),直至完成推断钻孔集合ND中所有推断钻孔的遍历,得到所有推断钻孔的虚拟地层推断。
8.根据权利要求7所述的面向地质钻孔的虚拟地层推断方法,其特征在于:步骤(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。
9.一种面向地质钻孔的虚拟地层推断装置,包括处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于:所述处理器执行所述程序时实现权利要求1-8中任意一项所述的方法。
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 true CN112837413A (zh) | 2021-05-25 |
CN112837413B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113240771A (zh) * | 2021-05-27 | 2021-08-10 | 广东省地震局 | 一种地层等厚线图的自动生成方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130297274A1 (en) * | 2011-01-27 | 2013-11-07 | Landmark Graphics Corporation | Methods and systems regarding models of underground formations |
CN109003330A (zh) * | 2018-07-02 | 2018-12-14 | 南京师范大学 | 一种基于基岩边界约束的三维地层建模方法 |
CN111429573A (zh) * | 2020-03-06 | 2020-07-17 | 南京师范大学 | 一种三维地质钻孔模型的自动构建方法 |
-
2021
- 2021-02-08 CN CN202110180713.0A patent/CN112837413B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130297274A1 (en) * | 2011-01-27 | 2013-11-07 | Landmark Graphics Corporation | Methods and systems regarding models of underground formations |
CN109003330A (zh) * | 2018-07-02 | 2018-12-14 | 南京师范大学 | 一种基于基岩边界约束的三维地层建模方法 |
CN111429573A (zh) * | 2020-03-06 | 2020-07-17 | 南京师范大学 | 一种三维地质钻孔模型的自动构建方法 |
Non-Patent Citations (1)
Title |
---|
田宜平;吴冲龙;董志;: "引入虚拟钻孔的地层连接推理与地质体建模", 地质科技情报, no. 05, 15 September 2010 (2010-09-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113240771A (zh) * | 2021-05-27 | 2021-08-10 | 广东省地震局 | 一种地层等厚线图的自动生成方法及装置 |
CN113240771B (zh) * | 2021-05-27 | 2021-12-07 | 广东省地震局 | 一种地层等厚线图的自动生成方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN112837413B (zh) | 2024-05-07 |
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 | |
CN111899341B (zh) | 长大带状线性工程三维地质快速建模方法 | |
CN112381937B (zh) | 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法 | |
US10853533B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
CN111079217B (zh) | 一种基于bim的岩土工程综合勘察信息解译方法及*** | |
CN109003330B (zh) | 一种基于基岩边界约束的三维地层建模方法 | |
CN106844871A (zh) | 基于bim的三维岩土工程勘察信息模型构建方法 | |
AU2013402201B2 (en) | In-situ wellbore, core and cuttings information system | |
CN103824329B (zh) | 一种地质勘探三维可视化储量估算方法 | |
US10380793B2 (en) | Geobody surface reconstruction | |
US10387583B2 (en) | Rotations from gradient directions | |
US10113421B2 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
US20180348400A1 (en) | Method and System for Creating and Using a Subsurface Model in Hydrocarbon Operations | |
CN111459955A (zh) | 一种基于gis平台的三维地质结构模型自动更新方法及*** | |
CN112150582B (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
Jia et al. | A triangulated irregular network constrained ordinary Kriging method for three-dimensional modeling of faulted geological surfaces | |
CN116776422A (zh) | 一种盾构隧道三维正向设计方法 | |
CN110197013B (zh) | 基于Morphing的河床基岩面建模方法 | |
CN103871101A (zh) | 一种基于三维布尔运算的基坑与地质体模型无缝集成方法 | |
Mahmoudi et al. | Optimisation of geotechnical surveys using a BIM-based geostatistical analysis | |
CN112837413B (zh) | 面向地质钻孔的虚拟地层推断方法及装置 | |
Hegemann et al. | Hybrid ground data model for interacting simulations in mechanized tunneling | |
US10650107B2 (en) | Three-dimensional subsurface formation evaluation using projection-based area operations | |
CN116721227A (zh) | 一种复杂地质体三维地质模型自动化建模方法 | |
CN112581558B (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 |