CN110673227B - 地层不整合交切的处理方法及处理装置 - Google Patents
地层不整合交切的处理方法及处理装置 Download PDFInfo
- Publication number
- CN110673227B CN110673227B CN201911050735.4A CN201911050735A CN110673227B CN 110673227 B CN110673227 B CN 110673227B CN 201911050735 A CN201911050735 A CN 201911050735A CN 110673227 B CN110673227 B CN 110673227B
- Authority
- CN
- China
- Prior art keywords
- top layer
- determining
- geological
- area
- positive
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本公开公开了一种地层不整合交切的处理方法及处理装置,属于地球物理勘探数据处理技术领域。所述处理方法包括:建立地质年代沉积顺序表;确定所有地层的顶层面的地质界面类型;计算得到每个顶层面的距离场,并根据每个顶层面的距离场确定该顶层面的正区域和负区域;获取所有顶层面之间的交线段;在交线段的约束下将各顶层面限定剖分为多个三角网格单元;确定顶层面的正区域和负区域;裁剪所有无效区域中的三角网格单元。本公开通过确定所有顶层面的正区域和负区域,从而可以根据地质年代沉积顺序表高效裁剪各顶层面的无效区域的三角网格单元,大大提高了复杂交切处理的工作效率。
Description
技术领域
本公开属于地球物理勘探数据处理技术领域,特别涉及一种地层不整合交切的处理方法及处理装置。
背景技术
在地质研究工作中,常需要建立有效的三维构造框架模型,以表征地层中的岩石物理属性分布。
相关技术中,三维构造框架模型的各地层之间往往会出现不整合交切的情况。在这种情况下,通常是通过空间三角网格曲面之间的两两求交经典算法获得交线,然后基于交线将地层的层面分成独立的几片,并重新剖分插值,最后按照交切在一起的两个地层的主辅关系交互剔除无效部分。
然而,上述方法,对于复杂的交切情况无法自动判断舍弃无效面,需要人工交互指定,从而导致工作效率低。
发明内容
本公开实施例提供了一种地层不整合交切的处理方法及处理装置,通过层面正负区判别模型及地质年代顺序表,可以高效裁剪掉无效层面,提高了不整合交切处理的工作效率。所述技术方案如下:
第一方面,本公开实施例提供了一种地层不整合交切的处理方法,所述处理方法包括:
根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表;
确定所有所述地层的顶层面的地质界面类型,所述地质界面类型包括剥蚀界面类型和自然沉积界面类型;
计算得到每个所述顶层面的距离场,所述顶层面的距离场指的是所有所述地层中所有点到所述顶层面的最短距离的合集;
根据每个所述顶层面的距离场确定每个所述顶层面的正区域和负区域;
获取所有所述顶层面之间的交线段,所述交线段为各顶层面之间相交所形成的线段;
在所述交线段的约束下将各所述顶层面限定剖分为多个三角网格单元;
当所述顶层面为剥蚀界面类型时,所述顶层面的正区域为无效区域,当所述顶层面为自然沉积界面类型时,所述顶层面的负区域为无效区域;
根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元。
可选地,所述计算得到每个所述顶层面的距离场,包括:
根据各所述顶层面的尺寸,确定规则笛卡尔网格单元的尺寸;
建立所有所述地层的笛卡尔网格系统,所述笛卡尔网格系统由多个所述规则笛卡尔网格单元组成;
根据所述地质年代沉积顺序表,将所有所述顶层面嵌入所述笛卡尔网格系统中;
计算被所述顶层面穿过的所有所述规则笛卡尔网格单元的每个顶点,到对应的所述顶层面的最近距离;
根据所述最近距离,计算得到每个所述顶层面的距离场。
可选地,所述根据每个所述顶层面的距离场确定每个所述顶层面的正区域和负区域,包括:
在所述笛卡尔网格系统中选取任一点作为识别点;
基于三线性插值,根据每个所述顶层面的距离场,分别计算得到所述识别点和每个所述顶层面之间的距离值;
确定所述识别点处于对应的所述顶层面的正区域或负区域;
根据每个识别点所处于各所述顶层面的位置,确定每个所述顶层面的正区域和负区域。
可选地,所述确定所述识别点处于对应的所述顶层面的正区域或负区域,包括:
如果所述识别点的距离值为正值,则所述识别点处于对应的所述顶层面的正区域,如果所述识别点的距离值为负值,则所述识别点处于对应的所述顶层面的负区域。
可选地,所述根据每个识别点所处于各所述顶层面的位置,确定每个所述顶层面的正区域和负区域,包括:
分别收集各个所述顶层面的距离值为正值的所有所述识别点,将所有距离值为正值的所述识别点所在的区域确定为所述顶层面的正区域;
分别收集各个所述顶层面的距离值为负值的所有所述识别点,将所有距离值为负值的所述识别点所在的区域确定为所述顶层面的负区域。
可选地,所述获取所有所述顶层面之间的交线段,包括:
确定每个所述顶层面的相交顶层面;
根据每个所述顶层面的相交顶层面,确定每个所述顶层面的交线段;
如果存在相交的所述交线段,则将相交的所述交线段的交点打断,并形成新的所述交线段。
可选地,所述在所述交线段的约束下将各所述顶层面限定剖分为多个三角网格单元,包括:
将每个所述顶层面划分为多个初级三角网格单元;
将具有所述交线段的所述初级三角网格单元标记为待裁剪网格单元;
根据所述待裁剪网格单元的顶点和交线段的端点,将所述待裁剪网格单元限定剖分为多个次级三角网格单元。
可选地,所述根据所述待裁剪网格单元的顶点和交线段的端点,将所述待裁剪网格单元限定剖分为多个次级三角网格单元,包括:
在所述待裁剪网格单元上建立局部平面坐标系;
根据所述笛卡尔网格系统的三维坐标系,将所述待裁剪网格单元的顶点和交线段的端点的三维坐标变换到所述局部平面坐标系中,以得到所述待裁剪网格单元的顶点和交线段的端点的平面坐标;
根据所述待裁剪网格单元的顶点和交线段的端点的平面坐标,在所述局部平面坐标系中对所述待裁剪网格单元进行限定剖分,得到所有所述次级三角网格单元的顶点的平面坐标;
将所有所述次级三角网格单元的顶点的平面坐标变换到所述笛卡尔网格系统的三维坐标系中,以得到所有所述次级三角网格单元的顶点的三维坐标。
可选地,所述根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元,包括:
根据所述地质年代沉积顺序表,确定各所述地层的年代排序;
根据各所述地层的年代排序,由古代到近代的顺序依次裁剪所有所述无效区域中的所述三角网格单元。
第二方面,本公开实施例提供了一种地层不整合交切的处理装置,所述处理装置包括:
地层排序模块,用于根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表;
地质界面类型确定模块,用于确定所有所述地层的顶层面的地质界面类型,所述地质界面类型包括剥蚀界面类型和自然沉积界面类型;
距离场计算模块,用于计算得到每个所述顶层面的距离场;
正负区确定模块,用于根据每个所述顶层面的距离场确定每个所述顶层面的正区域和负区域;
交线段获取模块,获取所有所述顶层面之间的交线段;
限定剖分模块,用于在所述交线段的约束下将各所述顶层面限定剖分为多个三角网格单元;
区域确定模块,用于当所述顶层面为剥蚀界面类型时,将所述顶层面的正区域确定为无效区域,当所述顶层面为自然沉积界面类型时,将所述顶层面的负区域确定为无效区域;
裁剪模块,用于根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元。
本公开实施例提供的技术方案带来的有益效果是:
通过本公开实施例提供的地层不整合交切的处理方法在处理层面不整合交切时,首先,根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表,地质年代沉积顺序表符合地层的沉积发展规律,通过地质年代沉积顺序表便于后续按照地层的沉积发展规律对各顶层面进行裁剪。然后,确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型,从而利用剥蚀界面类型和自然沉积界面类型可以模拟所有种类的地层顶层面不整合交切情况。再然后,计算得到每个顶层面的距离场,并根据每个顶层面的距离场确定每个顶层面的正区域和负区域,从而通过距离场确定不同地质界面类型的顶层面的正区域或负区域,便于后续对无效区域进行裁剪。
接着,获取所有顶层面之间的交线段,从而确定各顶层面之间交线段的位置,进而便于后续对交线段进行限定剖分。再接着,在交线段的约束下将各顶层面限定剖分为多个三角网格单元,从而实现对各顶面的限定剖分,进而使得交线段限定剖分至三角网格单元中,也就实现了对交线段的限定剖分。然后,确定无效区域,当顶层面为剥蚀界面类型时,顶层面的正区域为无效区域,当顶层面为自然沉积界面类型时,顶层面的负区域为无效区域,根据顶层面的不同地质界面类型,分别确定出所有的无效面,以供后续步骤裁剪。最后,根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元,从而得到裁剪后的适于地层的沉积发展规律的顶层面交切模型。
也就是说,通过确定所有顶层面的正区域和负区域,从而可以根据地质年代沉积顺序表高效裁剪各顶层面的无效区域的三角网格单元,大大提高了复杂交切处理的工作效率。
附图说明
为了更清楚地说明本公开实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本公开实施例提供的一种地层不整合交切的处理方法的方法流程图;
图2是本公开实施例提供的另一种地层不整合交切的处理方法的方法流程图;
图3是本公开实施例提供的地质年代沉积顺序表;
图4是本公开实施例提供的一种地层的顶层面的地质界面类型的示意图;
图5是本公开实施例提供的交线段的结构示意图;
图6是本公开实施例提供的在局部坐标系对交线段的三角单元进行重新限定剖分的示意图;
图7是本公开实施例提供的地层的顶层面交切处理后的结构示意图;
图8是本公开实施例提供的地层的顶层面为自然沉积界面类型时的裁剪示意图;
图9是本公开实施例提供的地层的顶层面为剥蚀界面类型时的裁剪示意图;
图10是本发明实施例提供的正逆断层及多值地质体的不整合交切示意图;
图11是本公开实施例提供的一种地层不整合交切的处理装置的结构示意图。
图中各符号表示含义如下:
10、初级三角网格单元;20、次级三角网格单元;100、待裁减顶层面;a、交线段;b、交点;1、地层排序模块;2、地质界面类型确定模块;3、距离场计算模块;4、正负区确定模块;5、交线段获取模块;6、限定剖分模块;7、区域确定模块;8、裁剪模块。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚,下面将结合附图对本公开实施方式作进一步地详细描述。
本公开实施例提供了一种地层不整合交切的处理方法,如图1所示,该处理方法包括:
S101、根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表。
S102、确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型。
S103、计算得到每个顶层面的距离场。
S104、根据每个顶层面的距离场确定每个顶层面的正区域和负区域。
S105、获取所有顶层面之间的交线段a。
S106、在交线段a的约束下将各顶层面限定剖分为多个三角网格单元。
S107、根据每个顶层面的地质界面类型,确定每个顶层面的无效区域。
步骤S107可以通过以下方式实现:当顶层面为剥蚀界面类型时,顶层面的正区域为无效区域,当顶层面为自然沉积界面类型时,顶层面的负区域为无效区域。
S108、根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元。
通过本公开实施例提供的地层不整合交切的处理方法在处理层面不整合交切时,首先,根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表,地质年代沉积顺序表符合地层的沉积发展规律,通过地质年代沉积顺序表便于后续按照地层的沉积发展规律对各顶层面进行裁剪。然后,确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型,从而利用剥蚀界面类型和自然沉积界面类型可以模拟所有种类的地层顶层面不整合交切情况。再然后,计算得到每个顶层面的距离场,并根据每个顶层面的距离场确定每个顶层面的正区域和负区域,从而通过距离场确定不同地质界面类型的顶层面的正区域或负区域,便于后续确定需要裁剪的无效区域。
接着,获取所有顶层面之间的交线段,从而确定各顶层面之间交线段的位置,进而便于后续对交线段进行限定剖分。再接着,在交线段的约束下将各顶层面限定剖分为多个三角网格单元,从而便于后续对无效区域的准确裁剪。然后,确定无效区域,当顶层面为剥蚀界面类型时,顶层面的正区域为无效区域,当顶层面为自然沉积界面类型时,顶层面的负区域为无效区域,根据顶层面的不同地质界面类型,分别确定出所有的无效面,以供后续步骤裁剪。最后,根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元,从而得到裁剪后的适于地层的沉积发展规律的顶层面交切模型。
也就是说,通过确定所有顶层面的正区域和负区域,从而可以根据地质年代沉积顺序表高效裁剪各顶层面的无效区域的三角网格单元,大大提高了复杂交切处理的工作效率。
图2是本公开实施例提供的另一种地层不整合交切的处理方法的方法流程图,如图2所示,该处理方法包括:
S201、根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表。
步骤S201包括:
a、确定每个地层发育的地质年代时间。
b、根据每个地层发育的地质年代时间,对各地层进行排序。
c、根据各地层的排序结果,建立地质年代沉积顺序表。
在上述实施方式中,建立地质年代沉积顺序表,符合地层的沉积发展规律,从而便于后续对无效区域三角单元进行裁剪。
图3是本公开实施例提供的地质年代沉积顺序表,如图3所示,地层Layer1对应的顶层面为H1,地层Layer2对应的顶层面为H2,以此类推,地层Layer4的顶层面为H4,另外,地层Layer4的底层面为H5。由古代到近代的顺序从下至上排列各地层。
S202、确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型。
在上述实施方式中,利用剥蚀界面类型和自然沉积界面类型可以模拟所有种类的地层的顶层面不整合情况。
图4是本公开实施例提供的一种地层的顶层面的地质界面类型的示意图,如图4所示,顶层面H2是一个剥蚀界面类型(EROSION),该顶层面之下的地层的顶层面H3、顶层面H4、顶层面H5超出顶层面H2,将这样的顶层面视为剥蚀界面类型(EROSION)。顶层面H1、顶层面H3、顶层面H4、顶层面H5、顶层面H6、顶层面H7、顶层面H8及顶层面H9这样自然沉积的顶层面视为自然沉积界面类型(ONLAP)。
S203、根据各顶层面的尺寸,确定规则笛卡尔网格单元的尺寸。
在上述实施方式中,通过笛卡尔网格单元便于模拟各地层。
需要说明的是,规则笛卡尔网格单元的尺寸为各顶层面的三角网格单元的平均尺寸。
S204、建立所有地层的笛卡尔网格系统,笛卡尔网格系统由多个规则笛卡尔网格单元组成。
在上述实施方式中,建立所有地层的笛卡尔网格系统,从而便于确定各顶层面的位置。
S205、根据地质年代沉积顺序表,将所有顶层面嵌入笛卡尔网格系统中。
S206、计算被顶层面穿过的所有规则笛卡尔网格单元的每个顶点,到对应的顶层面的最近距离。
在上述实施方式中,通过穿过顶层面的笛卡尔网格单元顶点到顶层面的最近距离,便于后续计算出不穿过地层的顶层面的笛卡尔网格单元到顶层面的最近距离,从而确定每个顶层面的距离场。
S207、根据最近距离,计算得到每个顶层面的距离场。
可选地,采用欧式距离变换方法计算得到每个顶层面的距离场。
在上述实施方式中,采用欧式距离变换方法计算得到每个顶层面的距离场,便于快速求得所有顶层面的距离场。
S208、在笛卡尔网格系统中选取任一点作为识别点。
在上述实施方式中,通过计算一点的位置,从而区分各顶层面的正区域和负区域。
S209、基于三线性插值,根据每个顶层面的距离场,分别计算得到识别点和每个顶层面之间的距离值。
S210、确定识别点处于对应的顶层面的正区域或负区域。
步骤S210可以通过以下方式实现:如果识别点的距离值为正值,则识别点处于对应的顶层面的正区域,如果识别点的距离值为负值,则识别点处于对应的顶层面的负区域。
在上述实施方式中,通过三线性插值可以确定笛卡尔网格系统中任意一点相对于各顶层面所在的区域(正区域或负区域)。
示例性地,首先,判断识别点落在哪个笛卡尔网格中,再利用笛卡尔网格单元8个顶点的最近距离,基于三线性插值计算任意一点距顶层面的距离值,距离值大于零则视为该点在顶层面的正区域,距离值小于零则视为该点在顶层面的正区域,等于零视为该点在顶层面上。
S211、分别收集各个顶层面的距离值为正值的所有识别点,将所有距离值为正值的识别点所在的区域确定为顶层面的正区域,分别收集各个顶层面的距离值为负值的所有识别点,将所有距离值为负值的识别点所在的区域确定为顶层面的负区域。
在上述实施方式中,通过将各层面的视为该层面正区域的识别点集合在一起,从而构成该顶层面的正区域;通过将各层面的视为该层面负区域的识别点集合在一起,从而构成该顶层面的负区域。
S212、确定每个顶层面的相交顶层面。
在上述实施方式中,某个顶层面的相交顶层面,即为与该顶层面存在相交关系的其他顶层面。
S213、根据每个顶层面的相交顶层面,确定每个顶层面的交线段a(参加图5)。
在上述实施方式中,顶层面与相交顶层面相交的部分,即形成交线段a。而存在交线段a的顶层面,则为和其他顶层面存在相交关系的顶层面,这种顶层面均存在层面不整合交切的情况,需要待裁剪。
S214、如果存在相交的交线段a,则将相交的交线段a的交点b打断,并形成新的交线段a(参见图6)。
在上述实施方式中,通过打断相交的交线段a,可以使得交线段a更为贴近实际的顶层面的不整合交切情况,从而便于后续将相交的交线段a重构。
S215、在交线段a的约束下将各顶层面限定剖分为多个三角网格单元(参见图6)。
步骤S215包括:
a、将每个顶层面划分为多个初级三角网格单元10。
b、将具有交线段a的初级三角网格单元10标记为待裁剪网格单元。
c、根据待裁剪网格单元的顶点和交线段a的端点,将待裁剪网格单元限定剖分为多个次级三角网格单元20。
在上述实施方式中,将每个顶层面划分为多个初级三角网格单元10,可以更好的模拟实际的顶层面,以便于后续步骤中的限定剖分。
根据待裁剪网格单元的顶点和交线段a的端点,将待裁剪网格单元限定剖分为多个次级三角网格单元20,可以理解为将待裁剪网格单元的顶点和交线段a的端点进行两两连接,以构成新的线段,通过这些线段,将待裁剪网格单元剖分为多个次级三角网格单元20。
可以理解的是,在对待裁剪网格单元进行限定剖分后,顶层面将由初级三角网格单元10和次级三角网格单元20组成,而初级三角网格单元10和次级三角网格单元20共同组成三角网格单元(参见图7)。
对待裁剪网格单元进行限定剖分的具体步骤包括:
(1)、在待裁剪网格单元上建立局部平面坐标系。
(2)、根据笛卡尔网格系统的三维坐标系,将待裁剪网格单元的顶点和交线段a的端点的三维坐标变换到局部平面坐标系中,以得到待裁剪网格单元的顶点和交线段a的端点的平面坐标。
(3)、根据待裁剪网格单元的顶点和交线段a的端点的平面坐标,在局部平面坐标系中对待裁剪网格单元进行限定剖分,得到所有次级三角网格单元20的顶点的平面坐标。
(4)、将所有次级三角网格单元20的顶点的平面坐标变换到笛卡尔网格系统的三维坐标系中,以得到所有次级三角网格单元20的顶点的三维坐标。
在上述实施方式中,由于待裁剪网格单元位于笛卡尔网格系统中,所有顶层面的坐标均为三维坐标。而三维坐标不利于对对待裁剪网格单元平面进行限定剖分,因此需要将待裁剪网格单元的三维坐标,变换为平面坐标。
由于每一个待裁剪单元可能不在同一平面上,所以需要针对每一个待裁剪网格单元,分别建立局部平面坐标系。具体的建立方式为,选取待裁剪网格单元的任一顶点作为坐标原点,在该待裁剪网格单元上建立局部平面坐标系。在建立了局部平面坐标系后,即可得到该待裁剪网格单元在局部平面坐标系中的平面坐标,从而便于对该待裁剪网格单元进行限定剖分。在进行限定剖分后,得到所有次级三角网格单元20的顶点的平面坐标,再将这些平面坐标转变为笛卡尔网格系统中的三维坐标,从而将次级三角网格单元20体现到笛卡尔网格系统中。
S216、根据每个顶层面的地质界面类型,确定每个顶层面的无效区域。
步骤S216可以通过以下方式实现:当顶层面为剥蚀界面类型时,顶层面的正区域为无效区域,当顶层面为自然沉积界面类型时,顶层面的负区域为无效区域。
图8是本公开实施例提供的地层的顶层面为自然沉积界面类型时的裁剪示意图,如图8所示,顶层面H2为自然沉积界面类型ONLAP,顶层面H2处于顶层面H3的负区域的三角网格单元设置为无效区域,即为图8中的待裁减顶层面100;顶层面H1为剥蚀界面类型EROSION,则基于顶层面H1的正区域和负区域,将顶层面H2处于顶层面H1的正区域的三角网格单元设置为无效区域,即为图8中的待裁减顶层面100。
图9是本公开实施例提供的地层的顶层面为剥蚀界面类型时的裁剪示意图,如图9所示,顶层面H1为剥蚀界面类型,则基于顶层面H1的正区域和负区域,将顶层面H2处于顶层面H1的正区域的三角网格单元设置为无效区域,即为图9中的待裁减顶层面100。紧接着再向上选择一个剥蚀界面类型的新地层的顶层面H0,重复将层面顶层面H2处于顶层面H0的正区域的三角网格单元设置为无效区域,即为图9中的待裁减顶层面100。
S217、根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元。
步骤S217包括:
a、根据地质年代沉积顺序表,确定各地层的年代排序。
b、根据各地层的年代排序,由古代到近代的顺序依次裁剪所有无效区域中的三角网格单元。
在上述实施方式中,由古代到近代的顺序裁剪每一个顶层面的无效区域中的三角网格单元,即可以将地质年代沉积顺序表与裁剪相结合,适应地层的发展沉积规律。
并且,由于在步骤S215中,已经在交线段a的约束下,将待裁剪网格单元限定剖分为多个次级三角网格单元20。所以只需裁剪掉无效区域中的所有次级三角网格单元20,即可完成无效区域的完整裁剪。
需要说明的是,本公开提供的层面不整合处理方法也可以处理正逆断层及多值地质体的不整合交切(见图10)。
图11是本公开实施例提供的一种地层不整合交切的处理装置的结构示意图,如图11所示,该处理装置包括:
地层排序模块1,用于根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表。
地质界面类型确定模块2,用于确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型。
距离场计算模块3,用于计算得到每个顶层面的距离场。
正负区确定模块4,用于根据每个顶层面的距离场确定每个顶层面的正区域和负区域。
交线段获取模块5,获取所有顶层面之间的交线段a。
限定剖分模块6,用于在交线段a的约束下将各顶层面限定剖分为多个三角网格单元。
区域确定模块7,用于当顶层面为剥蚀界面类型时,将顶层面的正区域确定为无效区域,当顶层面为自然沉积界面类型时,将顶层面的负区域确定为无效区域。
裁剪模块8,用于根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元。
通过本公开实施例提供的地层不整合交切的处理方法在处理层面不整合交切时,首先,根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表,地质年代沉积顺序表符合地层的沉积发展规律,通过地质年代沉积顺序表便于后续按照地层的沉积发展规律对各顶层面进行裁剪。然后,确定所有地层的顶层面的地质界面类型,地质界面类型包括剥蚀界面类型和自然沉积界面类型,从而利用剥蚀界面类型和自然沉积界面类型可以模拟所有种类的地层顶层面不整合交切情况。再然后,计算得到每个顶层面的距离场,并根据每个顶层面的距离场确定每个顶层面的正区域和负区域,从而通过距离场确定不同地质界面类型的顶层面的正区域或负区域,便于后续对无效区域进行裁剪。
接着,获取所有顶层面之间的交线段,从而确定各顶层面之间交线段的位置,进而便于后续对交线段进行限定剖分。再接着,在交线段的约束下将各顶层面限定剖分为多个三角网格单元,从而实现对各顶面的限定剖分,进而使得交线段限定剖分至三角网格单元中,也就实现了对交线段的限定剖分。然后,确定无效区域,当顶层面为剥蚀界面类型时,顶层面的正区域为无效区域,当顶层面为自然沉积界面类型时,顶层面的负区域为无效区域,根据顶层面的不同地质界面类型,分别确定出所有的无效面,以供后续步骤裁剪。最后,根据地质年代沉积顺序表,裁剪所有无效区域中的三角网格单元,从而得到裁剪后的适于地层的沉积发展规律的顶层面交切模型。
也就是说,通过确定所有顶层面的正区域和负区域,从而可以根据地质年代沉积顺序表高效裁剪各顶层面的无效区域的三角网格单元,大大提高了复杂交切处理的工作效率。
以上所述仅为本公开的较佳实施例,并不用以限制本公开,凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
Claims (8)
1.一种地层不整合交切的处理方法,其特征在于,所述处理方法包括:
根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表;
确定所有所述地层的顶层面的地质界面类型,所述地质界面类型包括剥蚀界面类型和自然沉积界面类型;
根据各所述顶层面的尺寸,确定规则笛卡尔网格单元的尺寸;
建立所有所述地层的笛卡尔网格系统,所述笛卡尔网格系统由多个所述规则笛卡尔网格单元组成;
根据所述地质年代沉积顺序表,将所有所述顶层面嵌入所述笛卡尔网格系统中;
计算被所述顶层面穿过的所有所述规则笛卡尔网格单元的每个顶点,到对应的所述顶层面的最近距离;
根据所述最近距离,计算得到每个所述顶层面的距离场,所述顶层面的距离场指的是所有所述地层中所有点到所述顶层面的最短距离的合集;
在所述笛卡尔网格系统中选取任一点作为识别点;
基于三线性插值,根据每个所述顶层面的距离场,分别计算得到所述识别点和每个所述顶层面之间的距离值;
确定所述识别点处于对应的所述顶层面的正区域或负区域;
根据每个识别点所处于各所述顶层面的位置,确定每个所述顶层面的正区域和负区域;
获取所有所述顶层面之间的交线段(a),所述交线段(a)为各顶层面之间相交所形成的线段;
在所述交线段(a)的约束下将各所述顶层面限定剖分为多个三角网格单元;
当所述顶层面为剥蚀界面类型时,所述顶层面的正区域为无效区域,当所述顶层面为自然沉积界面类型时,所述顶层面的负区域为无效区域;
根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元。
2.根据权利要求1所述的处理方法,其特征在于,所述确定所述识别点处于对应的所述顶层面的正区域或负区域,包括:
如果所述识别点的距离值为正值,则所述识别点处于对应的所述顶层面的正区域,如果所述识别点的距离值为负值,则所述识别点处于对应的所述顶层面的负区域。
3.根据权利要求2所述的处理方法,其特征在于,所述根据每个识别点所处于各所述顶层面的位置,确定每个所述顶层面的正区域和负区域,包括:
分别收集各个所述顶层面的距离值为正值的所有所述识别点,将所有距离值为正值的所述识别点所在的区域确定为所述顶层面的正区域;
分别收集各个所述顶层面的距离值为负值的所有所述识别点,将所有距离值为负值的所述识别点所在的区域确定为所述顶层面的负区域。
4.根据权利要求1所述的处理方法,其特征在于,所述获取所有所述顶层面之间的交线段(a),包括:
确定每个所述顶层面的相交顶层面;
根据每个所述顶层面的相交顶层面,确定每个所述顶层面的交线段(a);
如果存在相交的所述交线段(a),则将相交的所述交线段(a)的交点(b)打断,并形成新的所述交线段(a)。
5.根据权利要求1所述的处理方法,其特征在于,所述在所述交线段(a)的约束下将各所述顶层面限定剖分为多个三角网格单元,包括:
将每个所述顶层面划分为多个初级三角网格单元(10);
将具有所述交线段(a)的所述初级三角网格单元(10)标记为待裁剪网格单元;
根据所述待裁剪网格单元的顶点和交线段(a)的端点,将所述待裁剪网格单元限定剖分为多个次级三角网格单元(20)。
6.根据权利要求5所述的处理方法,其特征在于,所述根据所述待裁剪网格单元的顶点和交线段(a)的端点,将所述待裁剪网格单元限定剖分为多个次级三角网格单元(20),包括:
在所述待裁剪网格单元上建立局部平面坐标系;
根据所述笛卡尔网格系统的三维坐标系,将所述待裁剪网格单元的顶点和交线段(a)的端点的三维坐标变换到所述局部平面坐标系中,以得到所述待裁剪网格单元的顶点和交线段(a)的端点的平面坐标;
根据所述待裁剪网格单元的顶点和交线段(a)的端点的平面坐标,在所述局部平面坐标系中对所述待裁剪网格单元进行限定剖分,得到所有所述次级三角网格单元(20)的顶点的平面坐标;
将所有所述次级三角网格单元(20)的顶点的平面坐标变换到所述笛卡尔网格系统的三维坐标系中,以得到所有所述次级三角网格单元(20)的顶点的三维坐标。
7.根据权利要求1所述的处理方法,其特征在于,所述根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元,包括:
根据所述地质年代沉积顺序表,确定各所述地层的年代排序;
根据各所述地层的年代排序,由古代到近代的顺序依次裁剪所有所述无效区域中的所述三角网格单元。
8.一种地层不整合交切的处理装置,其特征在于,所述处理装置包括:
地层排序模块(1),用于根据各地层发育的地质年代时间顺序对所有地层进行排序,建立地质年代沉积顺序表;
地质界面类型确定模块(2),用于确定所有所述地层的顶层面的地质界面类型,所述地质界面类型包括剥蚀界面类型和自然沉积界面类型;
距离场计算模块(3),用于根据各所述顶层面的尺寸,确定规则笛卡尔网格单元的尺寸;
建立所有所述地层的笛卡尔网格系统,所述笛卡尔网格系统由多个所述规则笛卡尔网格单元组成;
根据所述地质年代沉积顺序表,将所有所述顶层面嵌入所述笛卡尔网格系统中;
计算被所述顶层面穿过的所有所述规则笛卡尔网格单元的每个顶点,到对应的所述顶层面的最近距离;
根据所述最近距离,计算得到每个所述顶层面的距离场;
正负区确定模块(4),用于在所述笛卡尔网格系统中选取任一点作为识别点;
基于三线性插值,根据每个所述顶层面的距离场,分别计算得到所述识别点和每个所述顶层面之间的距离值;
确定所述识别点处于对应的所述顶层面的正区域或负区域;
根据每个识别点所处于各所述顶层面的位置,确定每个所述顶层面的正区域和负区域;
交线段获取模块(5),获取所有所述顶层面之间的交线段(a);
限定剖分模块(6),用于在所述交线段(a)的约束下将各所述顶层面限定剖分为多个三角网格单元;
区域确定模块(7),用于当所述顶层面为剥蚀界面类型时,将所述顶层面的正区域确定为无效区域,当所述顶层面为自然沉积界面类型时,将所述顶层面的负区域确定为无效区域;
裁剪模块(8),用于根据所述地质年代沉积顺序表,裁剪所有所述无效区域中的所述三角网格单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911050735.4A CN110673227B (zh) | 2019-10-31 | 2019-10-31 | 地层不整合交切的处理方法及处理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911050735.4A CN110673227B (zh) | 2019-10-31 | 2019-10-31 | 地层不整合交切的处理方法及处理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110673227A CN110673227A (zh) | 2020-01-10 |
CN110673227B true CN110673227B (zh) | 2021-07-09 |
Family
ID=69085210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911050735.4A Active CN110673227B (zh) | 2019-10-31 | 2019-10-31 | 地层不整合交切的处理方法及处理装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110673227B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113625344A (zh) * | 2020-05-08 | 2021-11-09 | 中国石油天然气集团有限公司 | 地震数据可视化方法及装置 |
CN112464479A (zh) * | 2020-11-30 | 2021-03-09 | 核工业北京地质研究院 | 一种基于gocad的三维成矿预测方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010142659A2 (en) * | 2009-06-09 | 2010-12-16 | Shell Internationale Research Maatschappij B.V. | Method for stratigraphic analysis of seismic data |
CN102147934A (zh) * | 2011-04-19 | 2011-08-10 | 中南大学 | 基于tin的地质界面三维形态分析方法 |
CN102819865A (zh) * | 2012-08-09 | 2012-12-12 | 成都理工大学 | 一种大地电磁三维地质结构模型的建模方法 |
CN103632397A (zh) * | 2013-12-12 | 2014-03-12 | 薛林福 | 分块分单元三维地质调查方法 |
CN107886575A (zh) * | 2017-11-14 | 2018-04-06 | 辽宁工程技术大学 | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 |
CN108984919A (zh) * | 2018-07-20 | 2018-12-11 | 武汉中地数码科技有限公司 | 一种基于断层约束的三维地层面构建方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850845B2 (en) * | 2001-07-20 | 2005-02-01 | Tracy Joseph Stark | System for multi-dimensional data analysis |
US8861309B2 (en) * | 2011-01-31 | 2014-10-14 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
US11255994B2 (en) * | 2016-10-27 | 2022-02-22 | Schlumberger Technology Corporation | Earth model generation via measurements |
-
2019
- 2019-10-31 CN CN201911050735.4A patent/CN110673227B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010142659A2 (en) * | 2009-06-09 | 2010-12-16 | Shell Internationale Research Maatschappij B.V. | Method for stratigraphic analysis of seismic data |
CN102147934A (zh) * | 2011-04-19 | 2011-08-10 | 中南大学 | 基于tin的地质界面三维形态分析方法 |
CN102819865A (zh) * | 2012-08-09 | 2012-12-12 | 成都理工大学 | 一种大地电磁三维地质结构模型的建模方法 |
CN103632397A (zh) * | 2013-12-12 | 2014-03-12 | 薛林福 | 分块分单元三维地质调查方法 |
CN107886575A (zh) * | 2017-11-14 | 2018-04-06 | 辽宁工程技术大学 | 一种露天矿采场三角形网格裁剪煤层四边形网格的方法 |
CN108984919A (zh) * | 2018-07-20 | 2018-12-11 | 武汉中地数码科技有限公司 | 一种基于断层约束的三维地层面构建方法及系统 |
Non-Patent Citations (1)
Title |
---|
"地震DNA地层超剥点线识别技术及应用";罗红梅 等;《石油物探》;20160531;第55卷(第3期);第414-424页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110673227A (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104635262B (zh) | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 | |
US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
EP1397703B1 (en) | Automated system for modeling faulted multi-valued horizons | |
RU2600944C1 (ru) | Формирование моделей распознаваемых геологических структур на основании набора узловых точек | |
CN101158985B (zh) | 一种超维度河流动力学自适应并行监测的方法 | |
CN110673227B (zh) | 地层不整合交切的处理方法及处理装置 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN105652323B (zh) | 一种储层预测方法 | |
WO2014078355A1 (en) | Unstructured grids for modeling reservoirs | |
CN104317886A (zh) | 断层约束下网格节点插值时近邻条件数据点的搜索选取方法 | |
CN110263791B (zh) | 一种识别功能区的方法和装置 | |
CN114254960A (zh) | 一种低勘探程度区油气资源的评价方法 | |
CN103995944A (zh) | 一种空间异常极值无偏等值图绘制方法 | |
CN110910500B (zh) | 地层的提取方法及提取装置 | |
CN114036608A (zh) | 一种将地质剖面转化为极限平衡计算模型的方法 | |
CN109584369B (zh) | 实际地层全六面体网格生成方法及装置 | |
CN112035430B (zh) | 一种边界约束的自适应海底地形图水深点检测与筛选方法 | |
Dekkers et al. | Resource assessment: estimating the potential of an african rift geothermal reservoir | |
CN115469361B (zh) | 一种碎屑岩地层三维地质建模方法 | |
Liu et al. | Study on a computing technique suitable for true 3D modeling of complex geologic bodies | |
CN107730586A (zh) | 一种地层建模的方法和系统 | |
CN115903033A (zh) | 一种断点数据智能组合及曲面重建的方法 | |
CN110851930A (zh) | 基于WebGL的地表布线设计方法及系统 | |
CN105894580A (zh) | 三维地质表面模型中的曲面延伸数据处理方法 | |
CN118037990A (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 |