CN105654551A - 一种构建地质局部断层三维模型的方法及装置 - Google Patents
一种构建地质局部断层三维模型的方法及装置 Download PDFInfo
- Publication number
- CN105654551A CN105654551A CN201511009770.3A CN201511009770A CN105654551A CN 105654551 A CN105654551 A CN 105654551A CN 201511009770 A CN201511009770 A CN 201511009770A CN 105654551 A CN105654551 A CN 105654551A
- Authority
- CN
- China
- Prior art keywords
- triangle
- subset
- local
- topology
- true edge
- 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
- G06T17/05—Geographic models
-
- 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/30—Polynomial surface description
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例提供的一种构建地质局部断层三维模型的方法及装置。所述方法包括:根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;将所述局部拓扑子集合并至所述全局网格拓扑中。利用本发明实施例技术方案,对目标地层中每个离散点进行局部拓扑重建,将所有离散点重建的局部拓扑合并至全局网格拓扑中,生成目标地层的局部断层的三维模型。所述三维模型的生成避免人工干预,生成速度较快,精度较高,局部断层的接缝处结构平滑,更接近于实际地层结构。
Description
技术领域
本发明涉及地球物理勘探数据处理技术领域,尤其涉及一种构建地质局部断层三维模型的方法及装置。
背景技术
构建地质局部断层三维模型的是地球物理勘探中的重要技术领域,地质解释人员可以根据构建的地质层面信息,对目标层位有更加准确的认识。构建地质局部断层三维模型通常是根据地质层位解释数据建立三维地质层位曲面图。一般地,构建地质局部断层三维模型的难点在于构建方案的确定,如果没有局部逆断层的存在,每个层位在一个地震道上一般只有一个投影值,构建处理相对简单。但是由于局部逆断层的存在,导致地质解释数据的处理比较复杂。
现有技术中,对于包含局部逆断层的地质层面构建通常采用分块构建方法。所述分块构建方法通常根据给定的离散数据,人为地根据地质构造将所述离散数据分为多个数据块,数据块与数据块之间的边界点相同,所述边界点可以由人工确定。地质层面构建过程中,可以针对每个数据块进行三角网构建,最终将每个数据块的三角网拼接成一个整体,生成目标地层的地质层位曲面图。
现有技术方法涉及到人工干预,生成速度较慢,在数据变化较快的地方存在一定的误差。具体表现在生成的层位曲面图中,所述逆断层的接缝处不够平滑,与实际的逆断层结构有误差。因此,现有技术中亟需一种高速、高精度的地质层面构建方法。
发明内容
本申请的目的在于提供一种构建地质局部断层三维模型的方法和装置,可以提高局部断层三维模型的生成速度和成像精度。
为了实现上述目的,本发明提供了一种构建地质局部断层三维模型的方法及装置,所述方法及装置具体是这样实现的:
一种构建地质局部断层三维模型的方法,所述方法包括:
根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;
根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;
基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;
将所述局部拓扑子集合并至所述全局网格拓扑中。
可选的,在本发明的一个实施例中,所述基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集,包括:
从所述地质解释层位数据中确定单个离散点,从所述四面体网中选取包含所述离散点的三角形,生成三角形子集;
从所述三角形子集中获取与所述离散点距离最短的三角形顶点,将所述离散点和所述三角形顶点相连的边作为第一基准边;
从所述三角形子集中获取以所述第一基准边为共边的第一三角形队列;
从所述第一三角形队列中获取所述第一基准边的对角最大的第一三角形,并将所述第一三角形加入局部拓扑子集中;
从所述三角形子集中选取与所述对角最大的三角形共边所述第一基准边且夹角大于预设值的第二三角形队列;
从所述第二三角形队列中选取所述第一基准边的对角最大的第二三角形,并将所述第二三角形加入所述局部拓扑子集中。
可选的,在本发明的一个实施例中,所述方法还包括:
将所述第二三角形中包含所述离散点的另一条边作为第二基准边;
从所述三角形子集中选取与所述对角最大的三角形共边所述第二基准边且夹角大于预设值的第三三角形队列;
从所述第三三角形队列中选取所述第二基准边的对角最大的第三三角形,并将所述第三三角形加入所述局部拓扑子集中。
可选的,在本发明的一个实施例中,所述预设值的取值范围包括140度到160度。
可选的,在本发明的一个实施例中,所述预设值的取值为150度。
可选的,在本发明的一个实施例中,所述将所述局部拓扑子集合并至所述全局网格拓扑中包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。
可选的,在本发明的一个实施例中,所述方法还包括:
检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。
一种构建地质局部断层三维模型的装置,所述装置包括:
包围盒获取单元,用于根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;
四面体网构建单元,用于根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;
局部拓扑子集生成单元,用于基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;
合并单元,用于将所述局部拓扑子集合并至所述全局网格拓扑中。
可选的,在本发明的一个实施例中,所述局部拓扑子集生成单元包括:
三角形子集生成单元,用于从所述地质解释层位数据中确定单个离散点,从所述四面体网中选取包含所述离散点的三角形,生成三角形子集;
第一基准边确定单元,用于从所述三角形子集中获取与所述离散点距离最短的三角形顶点,将所述离散点和所述三角形顶点相连的边作为第一基准边;
第一三角形队列获取单元,用于从所述三角形子集中获取以所述第一基准边为共边的第一三角形队列;
第一三角形确定单元,用于从所述第一三角形队列中获取所述第一基准边的对角最大的第一三角形,并将所述第一三角形加入局部拓扑子集中;
第二三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第一基准边且夹角大于预设值的第二三角形队列;
第二三角形确定单元,用于从所述第二三角形队列中选取所述第一基准边的对角最大的第二三角形,并将所述第二三角形加入所述局部拓扑子集中。
可选的,在本发明的一个实施例中,所述局部拓扑子集生成单元还包括:
第二基准边确定单元,用于将所述第二三角形中包含所述离散点的另一条边作为第二基准边;
第三三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第二基准边且夹角大于预设值的第三三角形队列;
第三三角形确定单元,用于从所述第三三角形队列中选取所述第二基准边的对角最大的第三三角形,并将所述第三三角形加入所述局部拓扑子集中。
可选的,在本发明的一个实施例中,所述预设值的取值范围包括140度到160度。
可选的,在本发明的一个实施例中,所述预设值的取值为150度。
可选的,在本发明的一个实施例中,所述合并单元还包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。
可选的,在本发明的一个实施例中,所述装置还包括:
不相容性检查单元,用于检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。
由此可见,本发明一种构建地质局部断层三维模型的方法及装置的技术方案通过对目标地层中每个离散点进行局部拓扑重建,将所有离散点重建的局部拓扑合并至全局网格拓扑中,生成目标地层的局部断层的三维模型。所述三维模型的生成避免人工干预,生成速度较快,成像精度较高,局部断层的接缝处结构平滑,更接近于实际地层结构。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的构建地质局部断层三维模型方法的一种实施例的方法流程图;
图2是本实施例某个目标工区包含局部逆断层的地质解释层位数据;
图3是本发明实施例对图2所述目标工区的地质解释层位数据四面体话的拓扑网络;
图4是对图3生成的四面体网格中的单个离散点的局部拓扑重建示意图;
图5是对图2生成的四面体网格中的离散点进行合并后的全局网格拓扑示意图;
图6是对图5中的全局网格拓扑进行不相容检查后形成的拓扑结构示意图;
图7是本发明提供的构建地质局部断层三维模型装置的一种实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本发明提供的构建地质局部断层三维模型方法的一种实施例的方法流程图,如图1所示,所述方法包括:
S1:根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒。
三维地质模型的构建与所采用的三维地质解释层位数据密切相关,所述三维地质解释层位数据可以由一系列的离散点组成,所述三维地质解释层位数据可以包括所述离散点的三维坐标。图2是本实施例某个目标工区包含局部逆断层的地质解释层位数据。
地质解释层位数据的数据体庞大,结构复杂,本发明实施例中,可以通过对所述三维地质解释层位数据的分析处理,获取所述目标地层的三维数据包围盒以及全局网格拓扑。所述三维数据包围盒可以用于表示复杂的三维地质解释层位数据,可以实现简化运算的目的。在本发明实施例中,所述三维数据包围盒可以包括矩形包围盒,所述矩形包围盒可以表示所述三维地质解释层位数据的有效范围,如最大、最小三维坐标值。
S2:根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑。
根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑。在本发明实施例中,可以采用三维约束Delaunay四面体法建立构建所述目标地层的四面体网以及全局网格拓扑。在三维空间中,Delaunay网格所采用的单元为四面体,一般地,在三维空间中,Delaunay四面体是由一系列离散点组成,所述四面体的外接球均不包含除此四面体的顶点以外的网格点。Delaunay四面体法具有较大的灵活性,能够更好地逼近三维数据的边界。所述Delaunay四面体法具有较严格的数学理论证明作为基础,可以从初始的散点生成性态优化的四面体网格。
本发明实施例中,所述四面体网格可以通过计算机语言实现,具体的算法流程包括:
(1)遍历所有离散数据点,构建初始Delaunay四面体,所述初始Delaunay四面体可以包含离散点。构建空腔的方法可以包括:搜索所述三维数据包围盒中所有四面体,找出满足预设条件的四面体。
(2)将所述空腔内的四面体删除,将所述空腔的外表面顶点与离散点相连形成新的四面体;
(3)依次插入每个离散点,确定所述离散点的影响范围,重复步骤(2),生成Delaunay四面体。
图3是本发明实施例对图2所述目标工区的地质解释层位数据四面体话的拓扑网络。
S3:基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集。
本发明中的所述目标地层包含局部断层,相对于整个层面,所述局部断层的结构相对比较细微,尤其是断层处的接缝处。在本实施例中,对每个离散点进行局部拓扑重建,再将每个离散点的局部拓扑加入到全局网格拓扑中。具体生成局部拓扑子集的方法包括S21-S26:
S21:从所述地质解释层位数据中确定单个离散点di(x,y,z),从所述四面体网中选取包含所述离散点di(x,y,z)的三角形,生成三角形子集Ti;
S22:从所述三角形子集Ti中获取与所述离散点di(x,y,z)距离最短的三角形顶点dj(x,y,z),将所述离散点di(x,y,z)和所述三角形顶点dj(x,y,z)相连的边ix作为第一基准边;
S23:从所述三角形子集Ti中获取以所述第一基准边ix为共边的第一三角形队列Tx1;
S24:从所述第一三角形队列Tx1中获取所述第一基准边ix的对角最大的第一三角形ixy,并将所述第一三角形ixy加入局部拓扑子集Di(P)中;
S25:从所述三角形子集Ti中选取与所述对角最大的三角形共边所述第一基准边ix且夹角大于预设值的第二三角形队列Tx2;
S26:从所述第二三角形队列Tx2中选取所述第一基准边ix的对角最大的第二三角形ixy2,并将所述第二三角形ixy2加入所述局部拓扑子集Di(P)中。
在此之后,还包括S31-S33:
S31:将所述第二三角形ixy2中包含所述离散点di(x,y,z)的另一条边iz作为第二基准边;
S32:从所述三角形子集Ti中选取与所述对角最大的三角形共边所述第二基准边iz且夹角大于预设值的第三三角形队列Tx3;
S33:从所述第三三角形队列Tx3中选取所述第二基准边iz的对角最大的第三三角形ixy3,并将所述第三三角形ixy3加入所述局部拓扑子集Di(P)中。
本实施例中的预设值的取值范围包括140度到160度,优选的,所述预设值的取值可以为150度,此时,获得的全局网格拓扑的一致性较高,不相容的三角形最少。
图4是对图3生成的四面体网格中的单个离散点的局部拓扑重建示意图。本实施例步骤利用地质曲面特征对局部断层曲面的生成进行约束,约束条件可以包括选取基准边对角最大的三角形等,从而提高了局部断层的生成速度和成像精度。
S4:将所述局部拓扑子集合并至所述全局网格拓扑中。
所述将所述局部拓扑子集合并至所述全局网格拓扑中包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。图5是对图2生成的四面体网格中的离散点进行合并后的全局网格拓扑示意图。如图5所示,生成的全局网格拓扑的一致性较高,不相容的三角形较少。
在本发明的一个实施例中,所述方法后续还包:检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。图6是对图5中的全局网格拓扑进行不相容检查后形成的拓扑结构示意图。如图6所示,利用本发明实施例方法生成的全局网格拓扑中局部断层的接缝处结构平滑,精度较高。
本发明另一方面还提供一种构建地质局部断层三维模型的装置,图7是本发明提供的构建地质局部断层三维模型装置的一种实施例的模块结构示意图,结合附图7,该装置70包括:
包围盒获取单元71,用于根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;
四面体网构建单元72,用于根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;
局部拓扑子集生成单元73,用于基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;
合并单元74,用于将所述局部拓扑子集合并至所述全局网格拓扑中。
在本发明的一个实施例中,所述所述局部拓扑子集生成单元73,可以包括:
三角形子集生成单元,用于从所述地质解释层位数据中确定单个离散点,从所述四面体网中选取包含所述离散点的三角形,生成三角形子集;
第一基准边确定单元,用于从所述三角形子集中获取与所述离散点距离最短的三角形顶点,将所述离散点和所述三角形顶点相连的边作为第一基准边;
第一三角形队列获取单元,用于从所述三角形子集中获取以所述第一基准边为共边的第一三角形队列;
第一三角形确定单元,用于从所述第一三角形队列中获取所述第一基准边的对角最大的第一三角形,并将所述第一三角形加入局部拓扑子集中;
第二三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第一基准边且夹角大于预设值的第二三角形队列;
第二三角形确定单元,用于从所述第二三角形队列中选取所述第一基准边的对角最大的第二三角形,并将所述第二三角形加入所述局部拓扑子集中。
在本发明的另一个实施例中,所述所述局部拓扑子集生成单元73,还可以包括:
第二基准边确定单元,用于将所述第二三角形中包含所述离散点的另一条边作为第二基准边;
第三三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第二基准边且夹角大于预设值的第三三角形队列;
第三三角形确定单元,用于从所述第三三角形队列中选取所述第二基准边的对角最大的第三三角形,并将所述第三三角形加入所述局部拓扑子集中。
在本发明的另一个实施例中在本发明的另一个实施例中,所述合并单元74,还可以包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。
所述装置70,还可以包括:
不相容性检查单75,用于检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。
由此可见,本发明一种构建地质局部断层三维模型的方法及装置的技术方案通过对目标地层中每个离散点进行局部拓扑重建,将所有离散点重建的局部拓扑合并至全局网格拓扑中,生成目标地层的局部断层的三维模型。所述三维模型的生成避免了人工干预,生成速度较快,成像精度较高,局部断层的接缝处结构平滑,更接近于实际地层结构。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
本申请中各个实施例所涉及的上述描述仅是本申请中的一些实施例中的应用,在某些标准、模型、方法的基础上略加修改后的实施方式也可以实行上述本申请各实施例的方案。当然,在符合本申请上述各实施例的中所述的处理方法步骤的其他无创造性的变形,仍然可以实现相同的申请,在此不再赘述。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的装置或模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然,也可以将实现某功能的模块由多个子模块或子单元组合实现。
本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(ApplicationSpecificIntegratedCircuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC625D、AtmelAT91SAM、MicrochipPIC18F26K20以及SiliconeLabsC8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的硬件的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
Claims (14)
1.一种构建地质局部断层三维模型的方法,其特征在于,所述方法包括:
根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;
根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;
基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;
将所述局部拓扑子集合并至所述全局网格拓扑中。
2.根据权利要求1所述的一种构建地质局部断层三维模型的方法,其特征在于,所述基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集,包括:
从所述地质解释层位数据中确定单个离散点,从所述四面体网中选取包含所述离散点的三角形,生成三角形子集;
从所述三角形子集中获取与所述离散点距离最短的三角形顶点,将所述离散点和所述三角形顶点相连的边作为第一基准边;
从所述三角形子集中获取以所述第一基准边为共边的第一三角形队列;
从所述第一三角形队列中获取所述第一基准边的对角最大的第一三角形,并将所述第一三角形加入局部拓扑子集中;
从所述三角形子集中选取与所述对角最大的三角形共边所述第一基准边且夹角大于预设值的第二三角形队列;
从所述第二三角形队列中选取所述第一基准边的对角最大的第二三角形,并将所述第二三角形加入所述局部拓扑子集中。
3.根据权利要求2所述的一种构建地质局部断层三维模型的方法,其特征在于,还包括:
将所述第二三角形中包含所述离散点的另一条边作为第二基准边;
从所述三角形子集中选取与所述对角最大的三角形共边所述第二基准边且夹角大于预设值的第三三角形队列;
从所述第三三角形队列中选取所述第二基准边的对角最大的第三三角形,并将所述第三三角形加入所述局部拓扑子集中。
4.根据权利要求2或3所述的一种构建地质局部断层三维模型的方法,其特征在于,所述预设值的取值范围包括140度到160度。
5.根据权利要求4所述的一种构建地质局部断层三维模型的方法,其特征在于,所述预设值的取值为150度。
6.根据权利要求1所述的一种构建地质局部断层三维模型的方法,其特征在于,所述将所述局部拓扑子集合并至所述全局网格拓扑中包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。
7.根据权利要求1所述的一种构建地质局部断层三维模型的方法,其特征在于,还包括:
检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。
8.一种构建地质局部断层三维模型的装置,其特征在于,所述装置包括:
包围盒获取单元,用于根据对目标地层的地质解释层位数据的分析,获取所述目标地层的三维数据包围盒;
四面体网构建单元,用于根据所述三维数据包围盒构建所述目标地层的四面体网以及全局网格拓扑;
局部拓扑子集生成单元,用于基于所述四面体网,对所述地质解释层位数据中的离散点进行局部拓扑重建,生成局部拓扑子集;
合并单元,用于将所述局部拓扑子集合并至所述全局网格拓扑中。
9.根据权利要求8所述的一种构建地质局部断层三维模型的装置,其特征在于,所述局部拓扑子集生成单元包括:
三角形子集生成单元,用于从所述地质解释层位数据中确定单个离散点,从所述四面体网中选取包含所述离散点的三角形,生成三角形子集;
第一基准边确定单元,用于从所述三角形子集中获取与所述离散点距离最短的三角形顶点,将所述离散点和所述三角形顶点相连的边作为第一基准边;
第一三角形队列获取单元,用于从所述三角形子集中获取以所述第一基准边为共边的第一三角形队列;
第一三角形确定单元,用于从所述第一三角形队列中获取所述第一基准边的对角最大的第一三角形,并将所述第一三角形加入局部拓扑子集中;
第二三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第一基准边且夹角大于预设值的第二三角形队列;
第二三角形确定单元,用于从所述第二三角形队列中选取所述第一基准边的对角最大的第二三角形,并将所述第二三角形加入所述局部拓扑子集中。
10.根据权利要求8所述的一种构建地质局部断层三维模型的装置,其特征在于,所述局部拓扑子集生成单元还包括:
第二基准边确定单元,用于将所述第二三角形中包含所述离散点的另一条边作为第二基准边;
第三三角形队列获取单元,用于从所述三角形子集中选取与所述对角最大的三角形共边所述第二基准边且夹角大于预设值的第三三角形队列;
第三三角形确定单元,用于从所述第三三角形队列中选取所述第二基准边的对角最大的第三三角形,并将所述第三三角形加入所述局部拓扑子集中。
11.根据权利要求9或10所述的一种构建地质局部断层三维模型的装置,其特征在于,所述预设值的取值范围包括140度到160度。
12.根据权利要求11所述的一种构建地质局部断层三维模型的装置,其特征在于,所述预设值的取值为150度。
13.根据权利要求8所述的一种构建地质局部断层三维模型的装置,其特征在于,所述合并单元还包括:
对比所述局部拓扑子集中的三角形和所述全局网格拓扑中的三角形,将所述局部拓扑子集中不与所述全局网格拓扑中的三角形重复的三角形添加至所述全局网格拓扑中。
14.根据权利要求8所述的一种构建地质局部断层三维模型的装置,其特征在于,所述装置还包括:
不相容性检查单元,用于检查所述全局网格拓扑是否符合设置的相容性要求,删除所述全局网格拓扑中不符合相容性要求的三角形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009770.3A CN105654551B (zh) | 2015-12-29 | 2015-12-29 | 一种构建地质局部断层三维模型的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009770.3A CN105654551B (zh) | 2015-12-29 | 2015-12-29 | 一种构建地质局部断层三维模型的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105654551A true CN105654551A (zh) | 2016-06-08 |
CN105654551B CN105654551B (zh) | 2019-07-09 |
Family
ID=56477351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511009770.3A Active CN105654551B (zh) | 2015-12-29 | 2015-12-29 | 一种构建地质局部断层三维模型的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105654551B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108279445A (zh) * | 2017-12-18 | 2018-07-13 | 中国石油天然气集团公司 | 一种确定地质模型的方法和装置 |
CN109658505A (zh) * | 2018-11-05 | 2019-04-19 | 武汉中地数码科技有限公司 | 一种带断层面的三维空间构体方法及系统 |
CN111125942A (zh) * | 2018-10-31 | 2020-05-08 | 香港科技大学 | 用于三维单元结构建模和拓扑优化的b样条高清晰度单元水平集方法和计算机存储介质 |
CN112989698A (zh) * | 2021-03-11 | 2021-06-18 | 上海城勘信息科技有限公司 | 一种地层分类方法、装置、电子设备及存储介质 |
CN113938667A (zh) * | 2021-10-25 | 2022-01-14 | 深圳普罗米修斯视觉技术有限公司 | 基于视频流数据的视频数据传输方法、装置及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040246249A1 (en) * | 2001-12-10 | 2004-12-09 | Laurent Deny | Method, device and programme for three-dimensional modelling of a geological volume by 3d parametering of the geological domain |
US20110205844A1 (en) * | 2010-02-22 | 2011-08-25 | Landmark Graphics Corporation, A Haliburton Company | Systems and Methods for Modeling 3D Geological Structures |
CN103514630A (zh) * | 2013-10-16 | 2014-01-15 | 北京石油化工学院 | 断层构造三维建模方法 |
CN103646423A (zh) * | 2013-12-24 | 2014-03-19 | 中国科学院地质与地球物理研究所 | 一种三维地质建模的方法及装置 |
CN104299269A (zh) * | 2014-09-30 | 2015-01-21 | 长江岩土工程总公司(武汉) | 一种三维地质模型的地质剖面输出方法 |
CN104331931A (zh) * | 2014-10-10 | 2015-02-04 | 北京航空航天大学 | 一种面向油藏数值模拟的地质断裂层面pebi网格处理方法 |
-
2015
- 2015-12-29 CN CN201511009770.3A patent/CN105654551B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040246249A1 (en) * | 2001-12-10 | 2004-12-09 | Laurent Deny | Method, device and programme for three-dimensional modelling of a geological volume by 3d parametering of the geological domain |
US20110205844A1 (en) * | 2010-02-22 | 2011-08-25 | Landmark Graphics Corporation, A Haliburton Company | Systems and Methods for Modeling 3D Geological Structures |
CN103514630A (zh) * | 2013-10-16 | 2014-01-15 | 北京石油化工学院 | 断层构造三维建模方法 |
CN103646423A (zh) * | 2013-12-24 | 2014-03-19 | 中国科学院地质与地球物理研究所 | 一种三维地质建模的方法及装置 |
CN104299269A (zh) * | 2014-09-30 | 2015-01-21 | 长江岩土工程总公司(武汉) | 一种三维地质模型的地质剖面输出方法 |
CN104331931A (zh) * | 2014-10-10 | 2015-02-04 | 北京航空航天大学 | 一种面向油藏数值模拟的地质断裂层面pebi网格处理方法 |
Non-Patent Citations (2)
Title |
---|
刘刚: "基于断层轮廓数据的三维形体重构技术的研究", 《中国优秀硕士学位论文数据库》 * |
朱合华: "基于钻孔数据重构地层周围表面模型算法", 《计算机工程与应用》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108279445A (zh) * | 2017-12-18 | 2018-07-13 | 中国石油天然气集团公司 | 一种确定地质模型的方法和装置 |
CN108279445B (zh) * | 2017-12-18 | 2020-08-07 | 中国石油天然气集团公司 | 一种确定地质模型的方法和装置 |
CN111125942A (zh) * | 2018-10-31 | 2020-05-08 | 香港科技大学 | 用于三维单元结构建模和拓扑优化的b样条高清晰度单元水平集方法和计算机存储介质 |
CN111125942B (zh) * | 2018-10-31 | 2023-07-28 | 香港科技大学 | 用于三维单元结构建模和拓扑优化的b样条高清晰度单元水平集方法和计算机存储介质 |
CN109658505A (zh) * | 2018-11-05 | 2019-04-19 | 武汉中地数码科技有限公司 | 一种带断层面的三维空间构体方法及系统 |
CN112989698A (zh) * | 2021-03-11 | 2021-06-18 | 上海城勘信息科技有限公司 | 一种地层分类方法、装置、电子设备及存储介质 |
CN113938667A (zh) * | 2021-10-25 | 2022-01-14 | 深圳普罗米修斯视觉技术有限公司 | 基于视频流数据的视频数据传输方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105654551B (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105654551A (zh) | 一种构建地质局部断层三维模型的方法及装置 | |
US10254440B2 (en) | Process for constructing a volume mesh for modeling geological structures | |
EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
US10901118B2 (en) | Method and system for enhancing meshes for a subsurface model | |
US9182511B2 (en) | System and method for reproducibly extracting consistent horizons from seismic images | |
CN112381937B (zh) | 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法 | |
US20110310101A1 (en) | Pillar grid conversion | |
US11042676B2 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
US20030132934A1 (en) | Three dimensional geological model construction | |
US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
US20180031721A1 (en) | Method and System for Generating a Subsurface Model | |
CA3023655A1 (en) | Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof | |
US11536862B2 (en) | Automatic feature extraction from seismic cubes | |
KR20180056086A (ko) | 터널 정보 삼차원 가시화 모델 구성방법 및 시스템 | |
US20180031720A1 (en) | Method and System for Generating a Subsurface Model | |
MX2010013221A (es) | Distribucion de propiedades en un modelo volumetrico tridimensional utilizando un campo de continuidad maxima. | |
CN112395745A (zh) | 地下暗河储集体地质模型建立方法、处理设备 | |
CA3023651A1 (en) | Three-dimensional fracture abundance evaluation of subsurface formations | |
US20190025460A1 (en) | Reservoir Simulations with Fracture Networks | |
CN104200528A (zh) | 基于矢量闭合的三维建模方法 | |
US20160025877A1 (en) | Methods and systems for determining well drilling paths in a hydrocarbon field | |
CN103698810A (zh) | 混合网最小走时射线追踪层析成像方法 | |
US10422925B2 (en) | 2.5D stadia meshing | |
CN114384584A (zh) | 裂缝建模方法及装置 | |
CN110968930B (zh) | 一种地质体变属性插值方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |