一种地震层位的可视化显示方法、装置、电子设备及介质制造方法及图纸

本发明专利技术实施例提供一种地震层位的可视化显示方法、装置、电子设备及介质，用以提升石油勘探中储层预测时的3D可视化显示效率。地震层位的可视化显示方法，包括：根据地震层位平面点集建立地震层位的初始三角剖分面；在x维度、y维度上分别建立初始的索引链表；逐点插入各地震层位数据点，针对每一个新插入的地震层位数据点，基于被分割的三角剖分面及其纵向连接关系在x维度、y维度上分别横向更新索引链表；在x维度、y维度上分别纵向更新索引链表，将同一维度内各链表元素中相同名称的三角剖分面之间建立纵向连接关系；直至所有地震层位数据点均插入到索引链表中；基于最终的双向索引链表，生成地震层位的三角网格并进行可视化显示。示。示。

【技术实现步骤摘要】
一种地震层位的可视化显示方法、装置、电子设备及介质


[0001]本专利技术涉及石油勘探
，尤指一种地震层位的可视化显示方法、装置、电子设备及介质。

技术介绍

[0002]三角剖分算法在计算机可视化研究中一直处于非常重要的地位，它广泛地应用在地球物理勘探、地理系统、地质建模、图形图像分析等领域，能够提高可视化显示的效率，对石油勘探中开展储层预测具有十分重要的意义。
[0003]如何把一个散点集合剖分成不均匀的三角网格，这就是散点集的三角剖分问题，散点集的三角剖分，对数值分析以及图形学来说，都是极为重要的一项预处理技术。现有三角剖分算法的工作流程如图1a、1b所示。请参见图1a，假设
△
ABC为初始三角剖分面T0，初始的索引链表包含T0，在插入第一个数据点P1后，P1将初始三角剖分面T0分割成三个新的三角剖分面T1、T2、T3，插入P1后索引链表更新为T1
→
T2
→
T3；请参见图1b，在插入第二个数据点P2后，P2将三角剖分面T3分割成三个新的三角剖分面T4、T5、T6，插入P2后索引链表更新为T1
→
T2
→
T4
→
T5
→
T6；以此类推，逐点插入各数据点，从而将所有数据点都插入索引链表中。现有三角剖分算法，在大规模数据点的应用场景下，三角剖分的效率较低。
[0004]现阶段研究人员的关注点主要集中在现有三角剖分算法基础上进行的优化和改进。本专利技术人在专利技术过程中发现，现有技术中对于在数据预处理阶段通过建立索引网格对平面点集的数据范围进行划分并没有有效的办法，并不能索引出在网格单元内计算点、线、面之间的位置关系，导致3D可视化的显示速度慢。

技术实现思路

[0005]本专利技术实施例提供一种地震层位的可视化显示方法、装置、电子设备及介质，用以提升石油勘探中储层预测时的3D可视化显示效率。
[0006]为了解决上述技术问题，本说明书是这样实现的：
[0007]第一方面，提供了一种地震层位的可视化显示方法，包括：
[0008]根据地震层位平面点集建立地震层位的初始三角剖分面，所述地震层位平面点集中包括大量待三角剖分的地震层位数据点，所述初始三角剖分面是覆盖所有地震层位数据点的最大外围三角形；
[0009]在x维度、y维度上分别建立初始的索引链表，每个维度上初始的索引链表中包括一个初始的链表元素，所述初始的链表元素对应的数据范围为所述初始三角剖分面在该维度上的数据范围，所述初始的链表元素包含所述初始三角剖分面；
[0010]逐点插入各地震层位数据点，针对每一个新插入的地震层位数据点，确定该地震层位数据点所在的三角剖分面，使用该地震层位数据点将其所在的三角剖分面分割为三个新的三角剖分面，基于该地震层位数据点的坐标在x维度、y维度上分别增加一个链表元素，并基于被分割的三角剖分面及其纵向连接关系在x维度、y维度上分别横向更新索引链表；
其中，每个维度上的链表元素对应该维度上特定的数据范围、并在该链表元素包含的各三角剖分面之间建立横向连接关系，同一维度内各链表元素对应的数据范围之间无交集、且所有链表元素对应的数据范围的合集为该维度上初始的链表元素对应的数据范围，有至少一个点位于某个链表元素对应的数据范围内的三角剖分面属于该链表元素；以及
[0011]在x维度、y维度上分别纵向更新索引链表，将同一维度内各链表元素中相同名称的三角剖分面之间建立纵向连接关系；
[0012]直至地震层位平面点集中所有地震层位数据点均插入到索引链表中得到最终的双向索引链表；
[0013]基于最终的双向索引链表，生成地震层位的三角网格并进行可视化显示。
[0014]可选的，所述根据地震层位平面点集建立地震层位的初始三角剖分面，具体包括：
[0015]确定地震层位平面点集中所有地震层位数据点的数据范围；
[0016]基于所有地震层位数据点的数据范围建立覆盖所有地震层位数据点的最大外围三角形，作为所述初始三角剖分面；所述初始三角剖分面的三个顶点坐标分别表示为各顶点坐标具体通过如下公式确定：
[0017][0018][0019][0020][0021][0022][0023]其中，x
min
、x
max
分别表示x维度上所有地震层位数据点的数据范围的最小值、最大值，y
min
、y
max
分别表示y维度上所有地震层位数据点的数据范围的最小值、最大值，max(
·
)为取最大值操作，abs(
·
)为取绝对值操作。
[0024]可选的，所述基于该地震层位数据点的坐标在x维度、y维度上分别增加一个链表元素，具体包括：
[0025]在x维度上，基于各链表元素对应的数据范围，确定该地震层位数据点的横坐标所在的目标链表元素，使用所述横坐标将横坐标所在的目标链表元素分割为两个新的链表元素，并确定两个新的链表元素各自对应的数据范围；
[0026]在y维度上，基于各链表元素对应的数据范围，确定该地震层位数据点的纵坐标所在的目标链表元素，使用所述纵坐标将纵坐标所在的目标链表元素分割为两个新的链表元素，并确定两个新的链表元素各自对应的数据范围。
[0027]可选的，基于被分割的三角剖分面及其纵向连接关系在x维度、y维度上分别横向更新索引链表，具体包括：
[0028]在x维度上，将所述横坐标所在的目标链表元素包含的各三角剖分面作为两个新的链表元素初始包含的各三角剖分面，基于被分割的三角剖分面、分割后得到的三个新的
三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新；以及
[0029]基于被分割的三角剖分面的纵向连接关系确定x维度上包含所述被分割的三角剖分面的其他链表元素，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、其他链表元素对应的数据范围，对其他链表元素当前包含的各三角剖分面进行更新；
[0030]在y维度上，将所述纵坐标所在的目标链表元素包含的各三角剖分面作为两个新的链表元素初始包含的各三角剖分面，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新；以及
[0031]基于被分割的三角剖分面的纵向连接关系确定y维度上包含所述被分割的三角剖分面的其他链表元素，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、其他链表元素对应的数据范围，对其他链表元素当前包含的各三角剖分面进行更新。
[0032]可选的，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新，包括：
[0033]针对每个新的链表元素，将被分割的三角剖分面从其初始包含的各本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种地震层位的可视化显示方法，其特征在于，包括：根据地震层位平面点集建立地震层位的初始三角剖分面，所述地震层位平面点集中包括大量待三角剖分的地震层位数据点，所述初始三角剖分面是覆盖所有地震层位数据点的最大外围三角形；在x维度、y维度上分别建立初始的索引链表，每个维度上初始的索引链表中包括一个初始的链表元素，所述初始的链表元素对应的数据范围为所述初始三角剖分面在该维度上的数据范围，所述初始的链表元素包含所述初始三角剖分面；逐点插入各地震层位数据点，针对每一个新插入的地震层位数据点，确定该地震层位数据点所在的三角剖分面，使用该地震层位数据点将其所在的三角剖分面分割为三个新的三角剖分面，基于该地震层位数据点的坐标在x维度、y维度上分别增加一个链表元素，并基于被分割的三角剖分面及其纵向连接关系在x维度、y维度上分别横向更新索引链表；其中，每个维度上的链表元素对应该维度上特定的数据范围、并在该链表元素包含的各三角剖分面之间建立横向连接关系，同一维度内各链表元素对应的数据范围之间无交集、且所有链表元素对应的数据范围的合集为该维度上初始的链表元素对应的数据范围，有至少一个点位于某个链表元素对应的数据范围内的三角剖分面属于该链表元素；以及在x维度、y维度上分别纵向更新索引链表，将同一维度内各链表元素中相同名称的三角剖分面之间建立纵向连接关系；直至地震层位平面点集中所有地震层位数据点均插入到索引链表中得到最终的双向索引链表；基于最终的双向索引链表，生成地震层位的三角网格并进行可视化显示。2.如权利要求1所述的方法，其特征在于，所述根据地震层位平面点集建立地震层位的初始三角剖分面，具体包括：确定地震层位平面点集中所有地震层位数据点的数据范围；基于所有地震层位数据点的数据范围建立覆盖所有地震层位数据点的最大外围三角形，作为所述初始三角剖分面；所述初始三角剖分面的三个顶点坐标分别表示为各顶点坐标具体通过如下公式确定：各顶点坐标具体通过如下公式确定：各顶点坐标具体通过如下公式确定：各顶点坐标具体通过如下公式确定：各顶点坐标具体通过如下公式确定：各顶点坐标具体通过如下公式确定：其中，x
min
、x
max
分别表示x维度上所有地震层位数据点的数据范围的最小值、最大值，y
min
、y
max
分别表示y维度上所有地震层位数据点的数据范围的最小值、最大值，max(
·
)为取最大值操作，abs(
·
)为取绝对值操作。
3.如权利要求1所述的方法，其特征在于，所述基于该地震层位数据点的坐标在x维度、y维度上分别增加一个链表元素，具体包括：在x维度上，基于各链表元素对应的数据范围，确定该地震层位数据点的横坐标所在的目标链表元素，使用所述横坐标将横坐标所在的目标链表元素分割为两个新的链表元素，并确定两个新的链表元素各自对应的数据范围；在y维度上，基于各链表元素对应的数据范围，确定该地震层位数据点的纵坐标所在的目标链表元素，使用所述纵坐标将纵坐标所在的目标链表元素分割为两个新的链表元素，并确定两个新的链表元素各自对应的数据范围。4.如权利要求3所述的方法，其特征在于，基于被分割的三角剖分面及其纵向连接关系在x维度、y维度上分别横向更新索引链表，具体包括：在x维度上，将所述横坐标所在的目标链表元素包含的各三角剖分面作为两个新的链表元素初始包含的各三角剖分面，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新；以及基于被分割的三角剖分面的纵向连接关系确定x维度上包含所述被分割的三角剖分面的其他链表元素，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、其他链表元素对应的数据范围，对其他链表元素当前包含的各三角剖分面进行更新；在y维度上，将所述纵坐标所在的目标链表元素包含的各三角剖分面作为两个新的链表元素初始包含的各三角剖分面，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新；以及基于被分割的三角剖分面的纵向连接关系确定y维度上包含所述被分割的三角剖分面的其他链表元素，基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、其他链表元素对应的数据范围，对其他链表元素当前包含的各三角剖分面进行更新。5.如权利要求4所述的方法，其特征在于，所述基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、两个新的链表元素各自对应的数据范围，对两个新的链表元素初始包含的各三角剖分面分别进行更新，包括：针对每个新的链表元素，将被分割的三角剖分面从其初始包含的各三角剖分面中删除，根据该新的链表元素对应的数据范围判断分割后得到的三个新的三角剖分面是否属于该新的链表元素，并将属于该新的链表元素的新的三角剖分面加入被分割的三角剖分面的相应位置处，建立该新的链表元素最新包含的各三角剖分面之间的横向连接关系；所述基于被分割的三角剖分面、分割后得到的三个新的三角剖分面、其他链表元素对应的数据范围，对其他链表元素当前包含的各三角剖分面进行更新，包括：针对其他链表元素，将被分割的三角剖分面从其当前包含的各三角剖分面中删除，根据该其他链表元素对应的数据范围判断分割后得到的三个新的三角剖分面是否属于该其他链表元素，并将属于该其他链表元素的新的三角剖分面加入被分割的三角剖分面的相应位置处，建立该其他链表元素最新包含的各三角剖分面之间的横向连接关系。6.如权利要求1所述的方法，其特征在于，所述方法还包括：对于新插入的地震层位数据点，如果该地震层位数据点位于任一三角剖分面的任一条
边上，则不执行任何处理，继续插入下一个地震层位数据点。7.如权利要求1所述的方法，其特征在于，所述方法还包括：对地震资料解释的地震层位数据进行处理，忽略地震层位数据中的时间维度，得到所述地震层位平面点集。8.一种地震层位的可视化显示装置，其特征在于，包括：剖分面初始化单元，用于根据地震层位平面点集建立地震层位的初始三角剖分面，所述地震层位平面...

【专利技术属性】
技术研发人员：黄鑫鹏，宋强功，李磊，李炎，袁世洪，鞠鹏，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：