CN102903149A - 地质模型的成块成体方法以及装置 - Google Patents

地质模型的成块成体方法以及装置 Download PDF

Info

Publication number
CN102903149A
CN102903149A CN2012104051234A CN201210405123A CN102903149A CN 102903149 A CN102903149 A CN 102903149A CN 2012104051234 A CN2012104051234 A CN 2012104051234A CN 201210405123 A CN201210405123 A CN 201210405123A CN 102903149 A CN102903149 A CN 102903149A
Authority
CN
China
Prior art keywords
block
sub
piece
tomography
intersection
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
Application number
CN2012104051234A
Other languages
English (en)
Other versions
CN102903149B (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201210405123.4A priority Critical patent/CN102903149B/zh
Publication of CN102903149A publication Critical patent/CN102903149A/zh
Application granted granted Critical
Publication of CN102903149B publication Critical patent/CN102903149B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Tires In General (AREA)
  • Component Parts Of Construction Machinery (AREA)
  • Processing Or Creating Images (AREA)

Abstract

提供一种地质模型的成块成体方法和装置。对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面;根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层面,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线;基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体;根据相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线等信息,自动搜索父块、子块以及子体。

Description

地质模型的成块成体方法以及装置
技术领域
本申请涉及一种地质模型的成块成体方法以及装置,尤其涉及一种用于石油地震勘探中在执行二三维地质模型速度建模时自动成块成体的方法。
背景技术
在石油地震勘探领域中,叠前深度偏移是一种实现地质构造空间归位的有效技术之一,用于解决速度存在剧烈的横向变化、速度分界面不是水平层状时共反射点的叠加和绕射点的归位问题,能够使复杂构造或速度横向变化较大的地震资料正确成像。在叠前深度偏移过程中,速度建模的正确与精度的高低直接影响了偏移成像的效果。二三维的地质模型的成块成体又是叠前深度偏移速度模型建立的重要步骤。
传统的二三维地质模型成块成体对层状模型的描述存在如下特点:
1、描述复杂地质体,需要添加大量虚界面;
2、层位需进行复杂编号;
3、断层也采用层位线表示,不合一般习惯;
4、难以描述复杂地质结构。
近年来提出的块状模型采用封闭多边形描述地质体,具有逻辑清晰、操作简单的特定,可方便的构建任意复杂的地质模型。在块状模型构建中,多边形裁剪算法是其核心,它适用于凸的、凹的和带孔的多边形的裁剪。
然而,现有的模型描述方法存在以下的不足:现有技术处理的结果均为一个多边形,而在地质模型的处理中,多边形相互裁减的结果是多个多边形,即多个块。已有技术中的多边形裁剪算法不能很好地满足地质模型处理中的多边形相互裁减的需要。
现有的二三维地质模型成块成体在地质模型描述与模型生成中的流程如下:
1、对在地质剖面上拾取的层位/断层进行插值,然后对插值点进行三角网格化,形成空间的层位/断层面;
2、手工编辑或人为干预勾画出层位与断层,层位与层位,断层与断层的相交线;
3、通过相交线人工指定或采用多边形裁剪方法进行成块处理;
4、对块人工指定成体或采用复杂的算法对块封闭成体;
5、创建一个地质模型的实体,生成正确的地质模型。在模型生成方面,现有的模型生成方法多为人工干预的模型生成,有时需要人为的勾画出断层与层位的相交线,或者是勾画出断层的上下盘,缺乏自动成块成体的生成方法。
由以上的流程介绍可以看出,在现有的地质模型成块成体以及模型生成的各个步骤中,需要大量的手工编辑、勾画处理以及人为指定,工作量巨大、耗时长,并且从准确性上依赖于操作人员的经验和技能。
发明内容
本发明的目的在于提供一种地质模型的成块成体方法以及装置,在执行二三维地质模型速度建模时,自动地执行成块成体,从而实现地质模型建模的自动化或半自动化。
根据本发明的一方面,提供一种地质模型的成块成体方法,所述成块成体方法包括:对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面;根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层面,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线;基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体;将全部地质面上的边线和交线片段都投影至基面上并对边界求交,将每个外推并生成完整的地质面在基面上的投影识别为父块,以构建包含识别的各父块的父块集;为父块集中的每个父块构建子块集,其中,对于每个父块,自父块边界上的每个交点,根据以下算法搜索子块:自以其作为一个端点的交线片段开始遍历父块,从以所述交线片段的另一端点作为其他交线片段的端点的交线片段当中选取顺时针角度最小的交线片段,直到所述遍历回到开始的交点为止,将选取的交线片段形成的曲面确定为子块;基于构建的全部子块集搜索闭合子体,以构建子体集,其中,对每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块添加到子体候选子块集,再对子体候选子块集当中的每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块,继续添加到子体候选子块集,直到再无相互连接的子块为止,然后从子体候选子块集中选取闭合的体作为子体。
根据本发明的可选实施例,在所述外推操作中,根据用户手工编辑指定的封闭层位面/断层面以及边界面,勾画出相交层位与断层、层位与层位以及断层与断层之间的交线片段。
根据本发明的另一方面,提供一种地质模型的成块成体装置,包括:第一模块,用于对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面;第二模块,用于根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层面,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线;第三模块,用于基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体;第四模块,用于将全部地质面上的边线和交线片段都投影至基面上并对边界求交,将每个外推并生成完整的地质面在基面上的投影识别为父块,以构建包含识别的各父块的父块集;第五模块,用于为父块集中的每个父块构建子块集,其中,对于每个父块,自父块边界上的每个交点,根据以下算法搜索子块:自以其作为一个端点的交线片段开始遍历父块,从以所述交线片段的另一端点作为其他交线片段的端点的交线片段当中选取顺时针角度最小的交线片段,直到所述遍历回到开始的交点为止,将选取的交线片段形成的曲面确定为子块;第六模块,用于基于构建的全部子块集搜索闭合子体,以构建子体集,其中,对每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块添加到子体候选子块集,再对子体候选子块集当中的每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块,继续添加到子体候选子块集,直到再无相互连接的子块为止,然后从子体候选子块集中选取闭合的体作为子体。
根据本发明的可选实施例,第二模块在执行外推时,根据用户手工编辑指定的封闭层位面/断层面以及边界面,勾画出相交层位与断层、层位与层位以及断层与断层之间的交线片段。
附图说明
通过下面结合附图进行的描述,本发明的上述和其他目的和特点将会变得更加清楚,其中:
图1是示出根据本发明的示例性实施例的地质模型的成块成体方法的流程图;
图2是示出根据本发明的成块成体方法中执行父块搜索的示意图;
图3是示出根据本发明的成块成体方法中执行子块搜索的示意图;
图4是示出根据本发明的成块成体方法中执行子体搜索的示意图。
具体实施方式
以下,将参照附图来详细说明本发明的实施例。
可通过计算机软件模块实现本发明的成块成体方法。在本申请中描述的各个步骤不限制为上述步骤,其中的一些步骤可被进一步拆分为更多的步骤,并且一些步骤可合并为较少的步骤。
本发明还提供一种实现地质模型的成块成体的装置。
图1是示出根据本发明的示例性实施例的地质模型的成块成体方法的流程图。在所述实施例中,假设在通用计算机系统上通过实现各步骤的功能的软件模块执行所述成块成体方法。
参照图1,在步骤S110,计算机系统对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面。
在步骤S120,计算机系统根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线。
根据本发明的可选实施例,可半手工地执行步骤S120,即在所述外推操作中,手工编辑指定封闭层位面/断层面以及边界面,并勾画出相交层位与断层、层位与层位以及断层与断层之间的交线片段。
此后,在步骤S130,计算机系统基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体。也就是说,当完成步骤S130时,已形成具有多个块和体、面和线的三维立体。
在步骤S140,计算机系统执行父块搜索。具体为,将全部地质面上的边线和交线片段都投影至基面上并对边界求交,将每个外推并生成完整的地质面在基面上的投影识别为父块,以构建包含识别的各父块的父块集。图2示出通过执行步骤S140的操作搜索到的一个父块(由4个圆圈标注的点连接的交线片段形成的曲面)。
在步骤S150,计算机系统执行子块搜索。具体为,为父块集中的每个父块构建子块集,其中,对于每个父块,自父块边界上的每个交点,根据以下算法搜索子块:自以其作为一个端点的交线片段开始遍历父块,从以所述交线片段的另一端点作为其他交线片段的端点的交线片段当中选取顺时针角度最小的交线片段,直到所述遍历回到开始的交点为止,将选取的交线片段形成的曲面确定为子块。图3示出根据步骤S150遍历父块,搜索顺时针交角最小的交线片段形成的曲面作为子块的操作。
在步骤S160,计算机系统执行子体搜索。具体为,基于构建的全部子块集搜索闭合子体,以构建子体集;其中,对每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块添加到子体候选子块集,再对子体候选子块集当中的每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块,继续添加到子体候选子块集,直到再无相互连接的子块为止,然后从子体候选子块集中选取闭合的体作为子体。图4虽然没有提供三维的展示,但是示意性地示出根据步骤S150遍历全部子块集,搜索顺时针夹角最小的子块形成的曲面体作为子体的操作。
在通过步骤S140~S160执行的父块搜索、子块搜索以及子体搜索操作以后,确定了三维地质模型中各基面上的相互关系,取得关于各种交线/切割线、相交层位、父块、子块以及子体等各种信息,作为构建三维地质模型的重要信息。
从以上参照附图对本发明的示例性实施例的详细描述可以看出,本发明的地质模型的成块成体方法能够在对完成外推的立体空间自动地执行父块搜索、子块搜索以及子体搜索,从而使提高了实现地质模型的建模过程的自动化处理。
尽管已参照优选实施例表示和描述了本发明,但本领域技术人员应该理解,在不脱离由权利要求限定的本发明的精神和范围的情况下,可以对这些实施例进行各种修改和变换。

Claims (4)

1.一种地质模型的成块成体方法,所述成块成体方法包括:
对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面;
根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层面,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线;
基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体;
将全部地质面上的边线和交线片段都投影至基面上并对边界求交,将每个外推并生成完整的地质面在基面上的投影识别为父块,以构建包含识别的各父块的父块集;
为父块集中的每个父块构建子块集,其中,对于每个父块,自父块边界上的每个交点,根据以下算法搜索子块:自以其作为一个端点的交线片段开始遍历父块,从以所述交线片段的另一端点作为其他交线片段的端点的交线片段当中选取顺时针角度最小的交线片段,直到所述遍历回到开始的交点为止,将选取的交线片段形成的曲面确定为子块;
基于构建的全部子块集搜索闭合子体,以构建子体集,其中,对每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块添加到子体候选子块集,再对子体候选子块集当中的每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块,继续添加到子体候选子块集,直到再无相互连接的子块为止,然后从子体候选子块集中选取闭合的体作为子体。
2.如权利要求1所述的成块成体方法,其中,在所述外推操作中,根据用户手工编辑指定的封闭层位面/断层面以及边界面,勾画出相交层位与断层、层位与层位以及断层与断层之间的交线片段。
3.一种地质模型的成块成体装置,包括:
第一模块,用于对在地质剖面上拾取的层位/断层进行插值,并且对插值的点进行三角网格化,以形成立体空间的层位/断层面;
第二模块,用于根据预定的外推算法对形成的层位/断层逐个进行外推,指定其外推的封闭层位面/断层面,外推到的边界面作为封闭边界面,并且通过外推层位/断层面与其外推的封闭层位/断层面求交并截断形成交线片段,以生成完整的地质面并形成相交层位与断层、层位与层位以及断层与断层之间的交线片段或切割线;
第三模块,用于基于形成的切割线以及外推得到的地质面对立体空间进行切割,以形成多个块和体;
第四模块,用于将全部地质面上的边线和交线片段都投影至基面上并对边界求交,将每个外推并生成完整的地质面在基面上的投影识别为父块,以构建包含识别的各父块的父块集;
第五模块,用于为父块集中的每个父块构建子块集,其中,对于每个父块,自父块边界上的每个交点,根据以下算法搜索子块:自以其作为一个端点的交线片段开始遍历父块,从以所述交线片段的另一端点作为其他交线片段的端点的交线片段当中选取顺时针角度最小的交线片段,直到所述遍历回到开始的交点为止,将选取的交线片段形成的曲面确定为子块;
第六模块,用于基于构建的全部子块集搜索闭合子体,以构建子体集,其中,对每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块添加到子体候选子块集,再对子体候选子块集当中的每个子块,获取相互连接的子块,从所述相互连接的子块中选取与当前子块夹角最小的子块,继续添加到子体候选子块集,直到再无相互连接的子块为止,然后从子体候选子块集中选取闭合的体作为子体。
4.如权利要求3所述的成块成体装置,其中,第二模块在执行外推时,根据用户手工编辑指定的封闭层位面/断层面以及边界面,勾画出相交层位与断层、层位与层位以及断层与断层之间的交线片段。
CN201210405123.4A 2012-10-22 2012-10-22 地质模型的成块成体方法以及装置 Active CN102903149B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210405123.4A CN102903149B (zh) 2012-10-22 2012-10-22 地质模型的成块成体方法以及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210405123.4A CN102903149B (zh) 2012-10-22 2012-10-22 地质模型的成块成体方法以及装置

Publications (2)

Publication Number Publication Date
CN102903149A true CN102903149A (zh) 2013-01-30
CN102903149B CN102903149B (zh) 2015-05-27

Family

ID=47575361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210405123.4A Active CN102903149B (zh) 2012-10-22 2012-10-22 地质模型的成块成体方法以及装置

Country Status (1)

Country Link
CN (1) CN102903149B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977628A (zh) * 2014-04-11 2015-10-14 中国石油化工股份有限公司 地质层位片段自动组合方法
CN105612561A (zh) * 2013-08-16 2016-05-25 界标制图有限公司 识别并提取表示地质构造的一个或多个地质体中的断块
CN106461810A (zh) * 2014-06-05 2017-02-22 雪佛龙美国公司 构建复杂地球模型的系统和方法
CN106548513A (zh) * 2015-09-22 2017-03-29 中国石油化工股份有限公司 网格模型数据的生成方法
CN111966732A (zh) * 2020-08-12 2020-11-20 河北省地矿局第二地质大队 全智能化地质填图方法、装置、设备及可读存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279980A (zh) * 2010-06-13 2011-12-14 中国地质科学院矿产资源研究所 地质勘探矿体三维建模方法及其装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279980A (zh) * 2010-06-13 2011-12-14 中国地质科学院矿产资源研究所 地质勘探矿体三维建模方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孟宪海 等: "基于层面结构的三维闭合地质区块构造算法", 《北京航空航天大学学报》, vol. 31, no. 02, 28 February 2005 (2005-02-28), pages 182 - 185 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105612561A (zh) * 2013-08-16 2016-05-25 界标制图有限公司 识别并提取表示地质构造的一个或多个地质体中的断块
CN105684047A (zh) * 2013-08-16 2016-06-15 界标制图有限公司 动态更新表示一个或多个地质构造的区划
CN104977628A (zh) * 2014-04-11 2015-10-14 中国石油化工股份有限公司 地质层位片段自动组合方法
CN104977628B (zh) * 2014-04-11 2017-08-18 中国石油化工股份有限公司 地质层位片段自动组合方法
CN106461810A (zh) * 2014-06-05 2017-02-22 雪佛龙美国公司 构建复杂地球模型的系统和方法
CN106548513A (zh) * 2015-09-22 2017-03-29 中国石油化工股份有限公司 网格模型数据的生成方法
CN106548513B (zh) * 2015-09-22 2019-06-18 中国石油化工股份有限公司 网格模型数据的生成方法
CN111966732A (zh) * 2020-08-12 2020-11-20 河北省地矿局第二地质大队 全智能化地质填图方法、装置、设备及可读存储介质
CN111966732B (zh) * 2020-08-12 2022-08-05 河北省地矿局第二地质大队 全智能化地质填图方法、装置、设备及可读存储介质

Also Published As

Publication number Publication date
CN102903149B (zh) 2015-05-27

Similar Documents

Publication Publication Date Title
CN112381937B (zh) 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法
CN104574511B (zh) 一种快速递进的三维地质建模方法
CN101582173B (zh) 复杂地质构造块状模型构建方法
CN103514631B (zh) 基于钻孔数据的海量三维地质模型网格式并行构建方法
US10795053B2 (en) Systems and methods of multi-scale meshing for geologic time modeling
CN103646423B (zh) 一种三维地质建模的方法及装置
CN102867330B (zh) 基于区域划分的空间复杂层位重构方法
CN103824329B (zh) 一种地质勘探三维可视化储量估算方法
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
US20140136171A1 (en) Unstructured Grids For Modeling Reservoirs
CN102903149B (zh) 地质模型的成块成体方法以及装置
CN102385067B (zh) 含逆断层等值线绘图方法
CN104635262A (zh) 一种基于增强型矩形网格的正逆断层等值线自动生成方法
CN110189409B (zh) 一种基于plaxis的快速真三维地质建模方法及系统
CN113593030A (zh) 基于三维模型的地质剖面图生成方法、系统、终端及介质
CN105719342B (zh) 一种地裂缝地质体的三维建模可视化方法及装置
US20240185522A1 (en) Method for generating a hexahedral mesh
CN111951394A (zh) 基于地质图的断层构造单元三维模型构建方法及装置
CN103632396B (zh) 一种对复杂地质构造层序建模的方法
EP3371627B1 (en) Three-dimensional, stratigraphically-consistent seismic attributes
Jun-Qi et al. Study on a computing technique suitable for true 3D modeling of complex geologic bodies
CN110473270B (zh) 一种智能重建2d地质层序模型的方法
Guo et al. An automatic method for generating curvilinear geological section with stratal pinch-out
Yao et al. A method of geological surface reconstruction with reverse fault
Che et al. Study on 3D modeling method of faults based on GTP volume

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20180208

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200918

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.