CN1481511A - 用于分析和成像增强三维体数据集合的系统和方法 - Google Patents

用于分析和成像增强三维体数据集合的系统和方法 Download PDF

Info

Publication number
CN1481511A
CN1481511A CNA018207677A CN01820767A CN1481511A CN 1481511 A CN1481511 A CN 1481511A CN A018207677 A CNA018207677 A CN A018207677A CN 01820767 A CN01820767 A CN 01820767A CN 1481511 A CN1481511 A CN 1481511A
Authority
CN
China
Prior art keywords
data value
data sets
volume data
volume
data
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
CNA018207677A
Other languages
English (en)
Other versions
CN1214254C (zh
Inventor
杰克・李
杰克·李
・M・谢菲尔德
塔特姆·M·谢菲尔德
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.)
Landmark Graphics Corp
Original Assignee
Magic Earth Inc
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 Magic Earth Inc filed Critical Magic Earth Inc
Publication of CN1481511A publication Critical patent/CN1481511A/zh
Application granted granted Critical
Publication of CN1214254C publication Critical patent/CN1214254C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32—Transforming one recording into another or one representation into another
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • 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/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/64—Geostructures, e.g. in 3D data cubes
    • G01V2210/641—Continuity of geobodies

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Image Generation (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Processing Or Creating Images (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

公开了一种方法(10)和系统，用于通过将一个或多个属性体(12)结合成单一体来创建结合属性体或组合体。例如，可以使用地震数据体用于由多个地震属性体创建结合地震属性体。通过将地震数据体中标准8位数据值(12)的某几个以表示其他空间位置一致的地震属性数据的某几个其他数值的标记数值(18)进行替换来实现这一点。然后可以显示(22)组合体，以及在如地质体等感兴趣的事件上定位种子拾取。然后自动检出器(24)函数或程序将找到快速进一步确定事件的全部连接点。然后可以显示并译码事件。

Description

用于分析和成像增强三维体数据集合的系统和方法
技术领域
本发明一般涉及三维(“3D”)体数据集合(3D volume data sets)的成像。更具体地，本发明涉及代表3D体数据集合中的事件的物理属性的改进的成像和分析。
背景技术
很多研究领域需要3D体数据集合的分析和成像。例如，在医学领域，使用计算机轴向断层扫描(“CAT”)扫描仪或核磁共振成像(“MRI”)设备来产生患者身体一些部位的图像或诊断图像。扫描仪或MRI设备产生要被成像或显示的3D体数据集合，从而医学人员可以分析图像并做出诊断。
在涉及地球科学的很多研究领域中也应用三维体数据集合。地震探测是一种用于探索地球地下地质的方法。地下爆破或地震激发类似于低频声波、在地球表面以下传播并由地震仪检测的地震波。地震仪记录地震波的到达时间，包括直射波和反射波。已知爆破或地震的时间和地点，就可以计算通过地球内部的波的传播时间，并用以测量在地球内部的波的速度。类似的技术可以被用于近海石油和天然气的探索。在近海探索中，海船拖着声源和水下水听器。例如，由像气球爆裂一样工作的充气设备产生低频(如，50Hz)声波。声波由海底下面的岩层反弹，并由水听器拾音。以这种方式，反射波“映射”了诸如断层、地层褶曲和穹隆等储存石油的地下沉积结构。处理数据来产生包括在地理空间中特定(x、y、z)位置处的反射或地震振幅数据值的3D体数据集合。
3D体数据集合由具有x、y、z坐标的“体素”或体元素构成。在给定位置，每个体素具有用于一些测量或计算的物理特性的数据值。例如，数据值可以是可以包括256个可能值的八位数据字。地质数据值的示例包括振幅、相位、频率和相似度。在不同的3D体数据集合中存储不同的数据值，其中每个3D体数据代表不同的数据值。为了分析被称为“事件”的某些地质结构，必须译码来自不同3D体数据集合的信息，然后用以分析不同的事件。
显示多3D体数据集合的一种传统方法要求将体素重定比例，以便来自每个3D体数据集合的数据值适合用于彩色显示的256个数据值范围，这引起每个3D体数据集合准确度的下降。另一传统方法显示每个3D体数据集合，但是，通过调整透明度来控制每个3D体数据集合的光学显示。
在杰克·李所著、1999年3月在《前沿》上发表的题为《来自通过体素跟踪的种子拾取的结构缺陷》文章中，在单一显示中结合了每个只利用256个数据值范围中的128个数据值的两个3D体数据集合。显著地降低了显示分辨率，因此限制了准确译码某些事件的能力。
因此，在本领域中需要以改进的准确性和分辨率可视化一个或多个3D体数据集合的系统和方法。因此，本领域中的技术人员进行了长期地探索，而且将非常清楚致力于解决这些和其他问题的本发明。为了描述本发明起见，可能交替使用术语“成像”和“可视化”。
发明内容
因此，本发明的一个目的是提供一种用于可视化并译码单一结合3D体数据集合中的多个3D体数据集合的改进系统和方法。
本发明的另一目的是提供一种用于可视化并译码单一增强3D体数据集合中的单一3D体数据集合的改进系统和方法。
本发明的另一目的是提供一种利用在256个点上成比例的基本3D体数据集合用于可视化并译码一个或多个3D体数据集合的改进系统和方法，其中从可以被插入到基本3D体数据集合中而不改变基本3D体数据集合的定标的一个或多个3D体数据集合中选择数据值。
本发明的一个优势是选中的事件的改进分辨率。
本发明的另一优势是能够准确并有效地显示来自多于两个3D体数据集合的与事件相关的选中数据值。
本发明的另一优势是能够同时显示来自多个3D体数据集合的数据值。
本发明的另一优势是具有比透明度显示更好的准确性。
本发明的另一优势是能够在较低质量的数据值范围内集中关键事件。
本发明的另一优势是缩短译码周期时间。
通过附图、这里所给出地描述以及所附的权利要求，本发明的这些和其他目的、特征和优势将变得清楚。
因此，为了更准确而有效地分析和译码不同的选中事件，本发明提供了一种用于成像一个或多个3D体数据集合的系统和方法。每个3D体数据集合包括多个体素，其中每个体素包括定位在各个3D体数据集合中的3D位置的数据值。一个优选实施例包括一种结合多个3D体数据集合的方法，所述方法通过选择代表第一属性的第一3D体数据集合，选择代表第二属性的第二3D体数据集合，以及通过将第一3D体数据集合和第二3D体数据集合二者至少之一中的每个数据值与预选数据值范围或判别式进行比较而返回输出3D体数据集合来结合多个3D体数据集合。针对满足判别式的每个数据值，所述方法还包括在输出3D体数据集合中与相应数据值对应的位置插入第一选中数据值。针对不满足判别式的每个数据值，所述方法还包括在输出3D体数据集合中与相应数据值对应的位置插入第二选中数据值。所述方法还可以包括显示输出3D体数据集合的至少一个截面，以及通过在显示中插入种子拾取来选择数据值，从而用于可视化和译码事件。
第一选中数据值可以与第一属性相关，以及第二选中数据值可以与第二属性相关。利用输出3D体数据集合的显示，将种子拾取视觉定位在选中数据值。最好使用计算机和软件程序，用于确定或“自动拾取”与种子拾取相连、具有与相应的种子拾取相同或相似数据值的全部数据值。因此，本发明可以包括机器可读的程序存储设备，包含机器可执行的指令程序以最终成像输出3D体数据集合。
在优选实施例中，第一3D体数据集合和第二3D体数据集合每个都包括地震数据。所述方法也允许结合额外的3D体数据集合，以及从而，可以包括产生代表第三属性的第三3D体数据集合，并将其中的每个数据值与第二预选数据值范围进行比较。
在本发明的另一实施例中，可以使用与多个属性之一相关的增强3D体数据集合以可视化并译码不同的选中事件。在这个实施例中，所述方法包括确定来自代表具体属性的3D体数据集合的每个数据值。然后，通过选择数据值范围或判别式并将每个数据值与判别式相比较，从而创建增强3D体数据集合。如果满足判别式，那么所述方法还包括在增强3D体数据集合中与相应数据值对应的位置插入第一选中数据值。如果不满足判别式，那么所述方法包括在增强3D体数据集合中保持数据值不变。额外的步骤可以包括至少显示增强3D体数据集合的截面，通过在显示中插入种子拾取来选择数据值，以及自动拾取与种子拾取相连、具有与种子拾取的数据值相等的数据值的多个数据值。
在本发明的另一实施例中，提供了一种方法用于创建从多个3D体数据集合得到的结合3D体数据集合。所述方法包括选择基本3D体数据集合，其中基本3D体数据集合可以包括具有3D坐标和基本数据字的数据值。基本数据字最好可以与第一属性相关。此外，所述方法包括选择第二三维体数据集合，其中第二3D体数据集合可以包括数据值，所述数据值具有与基本3D体数据集合空间位置一致的坐标和与第二属性相关的第二数据字的。所述方法还包括通过在第二3D体数据集合中根据预选数据值范围或判别式选择数据值来返回结合3D体数据集合。如果满足判别式，那么所述方法还包括用与第二属性相关的数据值代替在基本3D体数据集合中相应坐标的基本数据字。如果不满足判别式，那么所述方法包括保持在基本3D体数据集合中的相应坐标与第一属性相关的基本数据字不变。在创建结合3D体数据集合之后，所述方法还可以包括至少显示结合3D体数据集合的一部分，以及利用显示在事件上定位种子拾取。在本发明的一个实施例中，所述事件是地质结构。
继续以这种方式，额外的方法步骤可以包括选择第三3D体数据集合，其中第三3D体数据集合可以包括具有与基本3D体数据集合空间位置一致的坐标和与第三属性相关的第三数据字的数据值。然后，所述方法包括通过在第三3D体数据集合中根据第二预选数据值范围或判别式选择数据值来返回修订的组合3D体数据集合。如果满足第二判别式，那么所述方法还包括用与第三属性相关的第二选中数据值代替在基本3D体数据集合中相应坐标的基本数据字。如果不满足第二判别式，那么所述方法还包括保持在基本3D体数据集合中相应位置与第一属性相关的基本数据字不变。在优选实施例中，所述第一属性、第二属性和第三属性每个都与地震数据相关。
附图说明
图1是描述了用于实现本发明的一个实施例的流程图；
图2是描述了典型地震轨迹与数据值或体素之间的关系的示意图；
图3是以相关的数据值直方图描述了给定范围-128到127之间的地震振幅数据值(八位数据值)的示例的示意图；
图4是描述了峰值(正相位)事件与相应数据值之间的关系的示意图；
图5是描述了用于自动拾取确定的数据值范围内的全部连接点的种子拾取的示意图；
图6是描述通过自动拾取处理画轮廓得到的地球体(geobody)的示意图；以及
图7示出了实现本发明的软件程序或系统结构的一个实施例。
尽管将连同目前的优选实施例描述本发明，可以理解的是，并不想要把本发明限制于这些实施例。相反，想要含盖包括在本发明的精神之内的全部选择对象、修改和等价物。
具体实施方式
通过将一个或多个属性结合到单一结合属性体或组合体中来创建结合体(“组合体”)。在与3D体数据集合内的地震属性相关的一个示例中，通过将某些数据值(在地震属性3D体数据集合中通常具有256个可能的数值)以代表其他空间位置一致的地震属性数据值的某些数值的另一数据值(如8位地震标记数据值)进行替换来实现这一点。组合体对于增强基于体素的自动跟踪装置的性能特别有用。针对在译码地震数据中的应用突出组合体的效用的示例包括用于自动跟踪低振幅不连续事件的地震振幅/瞬时相位组合体。事件意味着如地震数据描述的地质结构等地球体(geobody)。在另一示例中，地震振幅/相似度组合体可以用于在如断层或其他地质边界等地质间断处停止自动跟踪装置。在另一示例中，地震振幅/瞬时频率组合体可以用于突出地震描述的具体地质特征，如在峰值事件上的上超或在波谷事件上的上超。用于石油和天然气探索和开发的商业用途可以包括事件测绘、模型建立、多属性显示以及自动拾取增强。
现在参照图1，描述了与本发明一致的方法10用于确定与用于创建输出体数据集合的体素相关的数据值。正如本领域的技术人员所知的那样，体素包括3D坐标位置和数据值，如256个数值的数据字，即，8位字。在步骤12中，确定用于在相同坐标或点P的每个空间位置一致体的数据值VS、V1、…、VN，其中VS可以是点P处原始地震体的数据值，V1可以是点P处属性体1(VOL 1)的数据值，以及依此类推，VN是点P处属性体N(VOLN)的数据值。然后大体上如步骤14、16和18所述确定点P处组合体(VC)的数据值。
在步骤14中，确定VS、V1、…、VN是否处于用户指定的判别式中；例如，这些数据值是否在可以为VOL S、VOL 1、…、VOL N中的每一个选择的用户指定的数据值范围内。如步骤16所示，如果VS、V1、…、VN不处于判别式中，则设置VC的数据值；例如，设置为与点P处VS相同的数据值。可以使用其他用户指定的数据值，或者可以使用从相应点P处的任意3D体数据集合取得的数据值。如步骤18所示，如果VS、V1、…、VN处于判别式中，则将VC的数据值设置为与属性体VS、V1、…、VN之一相关的用户指定的标记数据值。如步骤20所示，然后将判别式应用于每个点P，即，用于要产生的组合体的每个体素。一旦这种反复处理完成，如步骤22所示，可以显示组合体的截面、切片或其他视图。然后，可以从视觉可辨别的事件中选择种子拾取，以及如步骤24中“自动拾取”处理所示将自动确认与种子拾取相连具有与种子拾取相同数据值的体素。这种方法快速而准确地确定如地质结构等事件的外延。然后可以显示事件的外延用于分析和解释。
用于增强自动跟踪或自动拾取子表面盐边界的组合体的一个实施例利用这些3D体数据集合，每个3D体数据集合代表如瞬时频率、地震振幅和瞬时振幅等不同的属性。可以利用传统方法来确定地震振幅数据(SA)的瞬时振幅属性(IA)和瞬时振幅属性的瞬时频率(IFIA)，配置这个用于检测和显示盐边界的组合体的实施例。因此，可以如下配置新的盐检测(SD)组合体：
如果f1≤IFIA(x，y，z)≤f2而且a1≤SA(x，y，z)≤a2，
则SD(x，y，z)＝v1，
否则SD(x，y，z)＝IA(x，y，z)。
f1、f2、a1、a2和v1的数值是用户选定的。
可以构建由地震振幅数据和瞬时相位数据组成的组合体的另一实施例，以增强自动跟踪另一地震事件的能力。利用传统方法从地震振幅数据中得到瞬时相位数据。结果得到具有与每个地震振幅数据值或体素相对应的瞬时相位数据值的3D体数据集合。在每个体素，瞬时相位数据值与用户指定的数据值或判别式进行比较。如果瞬时相位数据值满足判别式，则在新的组合体中以用户指定的标记数据值替代这个数据值。如果瞬时相位数据值不满足判别式，则以来自代表地震振幅的3D体数据集合的相应样本数据替代新的组合体中的这个数据值。
可以使用多于两个的3D体数据集合，以及为每一个3D体数据集合设置约束条件，考虑空间位置一致的数据体A(x，y，z)、B(x，y，z)和C(x，y，z)，其中如下配置数据体D(x，y，z)：
如果
a1≤A(x，y，z)≤a2而且/或b1≤B(x，y，z)≤b2而且/或c1≤C(x，y，z)≤c2，
则D(x，y，z)＝某些特定的用户选择数据值，
否则D(x，y，z)＝另一数值。
现在参照图2，示出了典型地震轨迹26与多个体素28之间的关系。沿地震轨迹26以预定间隔测量样本数据值比率30。体素28是围绕地震轨迹26的样本25的3D表示或盒子。为了地震数据的目的，体素通常可以具有可以从0到255或从-128到127标注的256个可能的数据值，或者由所用的数据类型确定的数据值范围。在体素28之间将用于任意3D体数据集合的测量数据值分开。
图3以相关数据直方图描述了给定范围-128到127的地震振幅数据值的示例。
图4描述了峰值(正相位)事件27与相关体素29之间的关系。
本发明的另一实施例包含增强3D体数据集合的创建。增强3D体数据集合用以增强自动跟踪装置创建表面、事件和或地球体(geobody)的能力。尽管本发明可以应用于任意类型的3D体数据集合，根据如下的描述使用地震相位数据值。增强3D体数据集合将8位地震数据体中256个数据值中的某些和与这些特定数据值相关的标记数值相结合。
为了创建代表地震相位的增强3D体数据集合，在最大峰值相位数据值附近选择数据值的范围。例如，当使用标度从0到255的峰值数据值时，其中峰值将是127的数据值，可能的范围可以从大约115到139。将测试3D体数据集合中的全部体素，而且可以用用户选择的数据值如127等重新定义任何具有在这个范围内的数据值的体素。如图5所示，然后可以在选中的事件中视觉定位种子拾取32。然后，计算机产生程序可以确定或自动拾取，并显示任何在上边界34和下边界36定义的范围内与种子拾取32相连的点38，可以包括或者排除如图5具体所示的数据值，画出作为应用于增强地震相位3D体数据集合的种子拾取和自动拾取处理的结果形成的地球体(geobody)的轮廓。
在另一实施例中，从多于一个的3D体数据集合中产生组合体。重要的是应当注意利用3D体数据集合的任意结合都可以应用本发明的方法，但是为了解释的目的，利用地震属性、地震相位和地震振幅数据体的结合提供地震示例。选择基本3D体数据集合。基本3D体数据集合将保持其直方图分布跨越振幅范围(-128到127)如图3中示例所示。在0到255的标度上，0代表-128波谷，而255代表127峰值。选择第二3D体数据集合，从中选择关键体素与基本3D体数据集合相结合。
在这个实施例中，参考如下描述使用地震相位数据。在最大峰值地震相位数据值附近选择数据值范围。当使用标度从0到255的峰值数据值时，其中峰值是数据值127，可能的范围可以从大约115到139。将在0到255的标度上，以数据值127(最大峰值)或255重新定义在这个范围内的全部体素。显示得到的组合体，以及然后在关键事件上定位种子拾取。如图5所示，然后自动检出器处理将找出全部连接点。
可以利用硬件、软件或者它们的组合实现本发明，以及可以在计算机系统和其他处理系统中实现本发明。在图7中示出了用于实现本发明的软件或程序结构100的一个实施例。在程序结构100底部的是操作系统102。合适的操作系统102包括，例如，UNIX操作系统、或者微软公司的Windows NT、或者作为相关领域的技术人员所清楚的其他操作系统。
菜单和视窗软件104位于操作系统102之上。菜单和视窗软件104用以提供多种菜单和窗口以便利与用户的交互作用，并获得用户的输入和指令。菜单和视窗软件104可以包括，例如，Microsoft WindowsTM、X WindowSystemTM(麻省理工学院注册的商标)、以及MOTIFTM(Open SoftwareFoundation公司的注册商标)。正如相关领域的技术人员所清楚的那样，也可以使用其他菜单和视窗软件。
基本图形库106位于菜单和视窗软件104之上。基本图形库106是用于3D计算机绘图的应用程序接口(API)。由基本图形库106执行的功能包括，例如，几何与屏面(raster)图元、RGBA或彩色索引模式、显示列表或直接模式、观察与建模变换、照明和阴影、隐藏表面移除、α合成(半透明)、防混淆、材质贴图、模糊效果(雾化、羽化、模糊)、反馈与选择、蜡版、以及累积缓冲。
特别优选的基本图形库106为加利福尼亚州，芒廷维尤市，SiliconGraphics公司(“SGI”)提供的OpenGL。OpenGLAPI是不依赖于硬件、窗口和操作系统的多平台工业标准。OpenGL被设计为可以从C、C++、FORTRAN、Ada和Java程序设计语言调入。OpenGL执行上述列出的用于基本图形库106的每项功能。OpenGL中的一些命令指定要绘出的几何对象，而其他命令控制如何处理对象。OpenGL状态的全部元素，甚至材质存储器和帧缓存器的内容，都可以通过使用OpenGL的客户应用程序获得。OpenGL和客户应用程序可以在相同或不同的机器上操作，因为OpenGL是网络透明的。在这里作为参考的《OpenGL程序设计指南》(ISBN：0-201-63274-8)和《OpenGL参考手册》(ISBN：0-201-63276-4)中更详细地描述了OpenGL。
视觉模拟图形库108位于基本图形库106之上。视觉模拟图形库108是用于创建实时、多重处理的3D视觉模拟图形应用的API。视觉模拟图形库108提供将如照明、材料、材质和透明度等图形库状态控制函数捆绑在一起的功能。这些函数跟踪状态以及随后可以进行绘制的显示列表的创建。
特别优选的视觉模拟图形库108为加利福尼亚州，芒廷维尤市，SGI提供的IRIS Performer。IRIS Performer支持上述的OpenGL图形库。IRISPerformer包括两个主库，libpf和libpr，和四个附属库，libpfdu、libpfdb、libpfui和libpfutil。
IRIS Performer的基础是性能绘制库libpr，基于GeoSet提供高速绘制功能和利用GeoState提供图形状态控制的低级库。GeoSet是将相同类型的图形图元(如，三角形或四边形)分组在一个数据对象中的可绘制几何形状的集合。GeoSet本身不包含几何图形，只包含指向数据阵列和索引阵列的指针。因为GeoSet中的全部图元类型相同，并且具有相同的属性，所以以最大硬件速度执行大多数数据库的绘制。GeoState提供用于GeoSet的图形状态定义(如，材质或材料)。
位于libpr之上的是libpf，提供优化利用多重处理硬件的高性能多重处理数据库绘制系统的实时视觉模拟环境。数据库实用程序库libpfdu提供用于定义3D对象的几何和外观属性的函数，共享状态和材料，以及从独立多边形输入产生三角带。数据库程序库libpfdb利用libpfdu、libpf和libpr的程序以输入大量工业标准数据库格式的数据库文件。libpfui是用户接口库，为用户接口(C和C++程序设计语言)提供用于写入处理组件的构件块。最后，libpfutil是实用程序库，提供用于实现如多通道选项支持等任务和图形用户接口(GUI)工具的例行程序。
使用IRIS Performer和OpenGLAPI的应用程序通常在为实时3D视觉模拟进行的准备中执行如下步骤：
1.初始化IRIS Performer；
2.指定图形管道的数目，选择多重处理配置，以及按照需要指定
  硬件模式；
3.初始化所选多重处理模式；
4.初始化帧速率，以及设置帧扩展策略；
5.按照要求创建、配置并打开窗口；以及
6.按照要求创建并配置显示通道。
一旦应用程序通过执行上面的步骤1到6已经创建了图形绘制环境，则应用程序通常每帧一次地重复主模拟循环。
7.计算动态性能，更新模型矩阵，等；
8.延迟到下一帧时间；
9.执行等待临界视点更新；
10.绘出帧。
本发明的组合/增强体程序110位于视觉模拟图形库108之上。程序110以相关领域的技术人员所知的方式与由视觉模拟和图形库108、基本图形库106、菜单和视窗软件104以及操作系统102中的每一个执行的函数相互作用，并利用这些函数。
更好地以面向对象的程序设计语言编写本发明的程序110以允许创建并使用对象和对象的功能性。特别优选的面向对象的程序设计语言是C++。
在一个实施例中，程序110以相关领域的技术人员众所周知的方式存储3D体数据集合。例如，用于数据体的格式可以由两部分组成，体头(volume header)，其后跟随与数据集合的大小一样长的数据体。体头通常包含在指定序列中的信息，如数据集合的文件路径(位置)、大小、x、y和z方向的尺寸、用于x、y和z轴的注释、用于数据值的注释等。数据体是字节的二进制序列，可以按照如下方式排序每数据值的一个或多个字节。第一字节是在体位置(x，y，z)＝(0，0，0)处的数据值。第二字节是在体位置(1，0，0)处的数据值，第三字节是在(2，0，0)处的数据值，等等。当取尽x维度时，然后增加y维度，以及最后增加z维度。本发明并不以任何方式限制具体数据格式。
程序110利用来自用户的输入来标识一个和多个3D体数据集合以用于成像和分析。当使用多个数据体时，用于每个数据体的数据值代表用于相同地理空间的不同物理参数和属性。作为示例，多个数据体可以包括地质体、温度体和水饱和度体。可以以(x，y，z，地震振幅)的形式表示地质体中的体素。可以以(x，y，z，℃)的形式表示温度体中的体素。可以以(x，y，z，％饱和度)的形式表示水饱和体中的体素。由这些体中的每个中的体素确定的物理或地理空间是相同的。但是，对于任意特定空间位置(x0，y0，z0)，将在地质体中包含地震振幅，在温度体中包含温度，以及在水饱和度体中包含水饱和度。
上述本发明的公开和描述是例证性和说明性的，而且本领域的技术人员应该清楚的是，在不偏离本发明的精神的前提下可以在大小、形状和材料、机制等价物的应用上，以及在所描述的结构或多种元件的特征的组合的细节上进行多种修改。

Claims (19)

1、一种用于在单一输出3D体数据集合中结合多个3D体数据集合的方法，每个3D体数据集合包括多个体素，每个体素包括数据值而且被定位在相应3D体数据集合的3D位置，所述方法包括：
选择代表第一属性的第一3D体数据集合；
选择代表第二属性的第二3D体数据集合；
将所述第一3D体数据集合和所述第二3D体数据集合二者至少之一中的每个所述数据值与预选数据值范围进行比较；
针对在所述数据值范围内的每个数据值，在所述输出3D体数据集合中与所述各个数据值对应的位置插入第一选中数据值；
针对在所述数据值范围之外的每个数据值，在所述输出3D体数据集合中与所述各个数据值对应的位置插入第二选中数据值；以及
显示所述输出3D体数据集合的至少一个截面。
2、按照权利要求1所述的方法，其特征在于所述第一选中数据值与所述第一属性相关。
3、按照权利要求1所述的方法，其特征在于所述第二选中数据值与所述第二属性相关。
4、按照权利要求1所述的方法，其特征在于还包括：
在所述显示中插入种子拾取，用于确定与物理现象相关的事件，利用所述显示在相应数据值位置定位所述种子拾取。
5、按照权利要求4所述的方法，其特征在于还包括：
自动拾取与所述种子拾取相连、在定位所述种子拾取处具有与所述相应数据值相同的数据值的全部数据值。
6、按照权利要求1所述的方法，其特征在于所述第一3D体数据集合和所述第二3D体数据集合包括地震数据，以及所述事件是地质结构。
7、按照权利要求1所述的方法，其特征在于还包括：
选择代表第三属性的第三3D体数据集合；
将所述第一3D体数据集合、所述第二3D体数据集合和所述第三3D体数据集合三者至少之一中的每个所述数据值与所述预选数据值范围进行比较；
针对在所述数据值范围内的每个数据值，在所述输出3D体数据集合中与所述各个数据值对应的位置插入所述第一选中数据值；
针对在所述数据值范围之外的每个数据值，在所述输出3D体数据集合中与所述各个数据值对应的位置插入所述第二选中数据值；以及
显示所述输出3D体数据集合的至少一个截面。
8、一种机器可读程序存储设备，它包含机器可执行指令程序以执行用于成像输出3D体数据集合的方法步骤，所述3D体数据集合包括多个体素，每个体素由3D坐标和数据值确定，所述方法包括：
选择多个空间位置一致的3D体数据集合，其中每个所述3D体数据集合至少与相应多个属性之一相关；
选择一个或多个可选择的判别式，用于与所述多个空间位置一致的数据集合中的一个或多个中的数据值相比较；
当所述数据值满足所述一个或多个可选择的判别式时，在所述输出3D体数据集合中与所述相应数据值相对应的坐标处插入第一预选数据值；以及
显示所述输出3D体数据集合的至少一个截面。
9、按照权利要求8所述的程序存储设备，其特征在于所述第一预选数据值至少与所述多个属性之一相关。
10、按照权利要求9所述的程序存储设备，其特征在于还包括：
当所述数据值满足所述一个或多个可选择的判别式时，在所述输出3D体数据集合中与所述相应数据值相对应的坐标处插入第二预选数据值；
11、按照权利要求10所述的程序存储设备，其特征在于所述第二预选数据值至少与所述多个属性中的另一个相关。
12、按照权利要求8所述的程序存储设备，其特征在于所述多个属性每个都与地震数据相关。
13、一种用于显示增强3D体数据集合的方法，所述增强3D体数据集合与利用包含多个体素的3D体数据集合的多个属性之一相关，每个体素包括数据值而且被定位在所述3D体数据集合中的3D位置，所述方法包括：
为每个数据值指定数据值范围，所述数据值范围与所述多个属性之一相关；
选择判别式，以及将每个数据值与所述判别式相比较；
当所述数据值满足所述判别式时，在所述增强3D体数据集合中与所述数据值相对应的位置插入预选数据值；
当所述数据值不满足所述判别式时，在所述增强3D体数据集合中与所述相应数据值相对应的位置插入所述数据值；以及
显示所述增强3D体数据集合的至少一个截面。
14、按照权利要求13所述的方法，其特征在于还包括：
利用所述显示用于在所述显示中示出的选中事件处插入种子拾取；以及
自动拾取与所述种子拾取相连、在定位所述种子拾取处具有与所述种子拾取相等的数据值的多个数据值，以用于自动确定所述事件的外延。
15、一种用于返回从多个3D体数据集合中得出的组合体的方法，包括：
选择基本3D体数据集合，所述基本3D体数据集合包括具有3D坐标和基本数据字的体素，所述基本数据字与第一属性相关；
选择第二3D体数据集合，所述第二3D体数据集合包括具有与基本3D体数据集合空间位置一致的坐标和与第二属性相关的第二数据字的体素；
根据第一预选数据值范围在所述第二3D体数据集合中选择体素；
当所述第二3D体数据集合中的所述相应体素在所述第一预选数据值范围内时，通过用与所述第二属性相关的第一预选数据值代替在所述基本3D体数据集合中的所述基本数据字来返回所述结合体；以及
至少显示部分所述结合体。
16、按照权利要求15所述的方法，其特征在于还包括：
利用所述显示在事件上定位种子拾取，以及
确定与所述种子拾取相连、具有与所述种子拾取相同的数据字的点。
17、按照权利要求16所述的方法，其特征在于所述事件是地质结构。
18、按照权利要求15所述的方法，其特征在于还包括：
选择第三3D体数据集合，所述第三3D体数据集合包括具有与基本3D体数据集合空间位置一致的坐标和与第三属性相关的第三数据字的体素；
根据第二预选数据值范围在所述第三3D体数据集合中选择体素；
当所述第三3D体数据集合中的所述相应体素在所述第二预选数据值范围内时，通过用与所述第三属性相关的第二预选数据值代替在所述基本3D体数据集合中的所述基本数据字来返回修订的结合体；以及
至少显示部分所述修订的结合体。
19、按照权利要求18所述的方法，其特征在于所述第一属性、所述第二属性和所述第三属性每个都与地震数据相关。
CNB018207677A 2001-01-31 2001-01-31 用于分析和成像增强三维体数据集合的系统和方法 Expired - Fee Related CN1214254C (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2001/003227 WO2002061463A1 (en) 2001-01-31 2001-01-31 System and method for analyzing and imaging an enhanced three-dimensional volume data set using one or more attributes

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100794416A Division CN100523869C (zh) 2001-01-31 2001-01-31 用于分析和成像增强三维体数据集合的系统和方法

Publications (2)

Publication Number Publication Date
CN1481511A true CN1481511A (zh) 2004-03-10
CN1214254C CN1214254C (zh) 2005-08-10

Family

ID=21742297

Family Applications (2)

Application Number Title Priority Date Filing Date
CNB2005100794416A Expired - Fee Related CN100523869C (zh) 2001-01-31 2001-01-31 用于分析和成像增强三维体数据集合的系统和方法
CNB018207677A Expired - Fee Related CN1214254C (zh) 2001-01-31 2001-01-31 用于分析和成像增强三维体数据集合的系统和方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CNB2005100794416A Expired - Fee Related CN100523869C (zh) 2001-01-31 2001-01-31 用于分析和成像增强三维体数据集合的系统和方法

Country Status (9)

Country Link
EP (1) EP1358510B1 (zh)
CN (2) CN100523869C (zh)
AU (1) AU2001234706B2 (zh)
BR (1) BR0116765C1 (zh)
CA (1) CA2428743C (zh)
DE (1) DE60121819T2 (zh)
MX (1) MXPA03004692A (zh)
NO (1) NO337451B1 (zh)
WO (1) WO2002061463A1 (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100349003C (zh) * 2004-12-29 2007-11-14 中国石油天然气集团公司 一种消除起伏地表成像误差的方法
CN101266299B (zh) * 2008-04-14 2011-03-30 林昌荣 利用地震数据体结构特征预测油气的方法
CN102160087A (zh) * 2007-01-05 2011-08-17 兰德马克绘图国际公司，哈里伯顿公司 用于实时可视化多卷数据集的方法及系统
CN103185896A (zh) * 2006-09-01 2013-07-03 哈利伯顿兰德马克绘图公司 用于对波形体进行成像的系统和方法
CN103221844A (zh) * 2010-09-01 2013-07-24 界标制图有限公司 绘制指示地下地层的特性的值的方法和系统
CN106056669A (zh) * 2016-06-07 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种海量地震数据的地质曲面重建方法
CN114333244A (zh) * 2021-12-31 2022-04-12 成都市美幻科技有限公司 多地震预警方法、装置及设备

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1550001B1 (en) * 2002-09-26 2013-03-27 ExxonMobil Upstream Research Company Method for performing stratigraphically-based seed detection in a 3-d seismic data volume
US6901332B2 (en) * 2002-11-22 2005-05-31 Western Geco, L.L.C. Technique for velocity analysis
CA2742850C (en) * 2008-11-07 2018-01-16 Landmark Graphics Corporation Systems and methods for computing and validating a variogram model
GB2503507B (en) 2012-06-29 2015-04-15 Foster Findlay Ass Ltd Adaptive fault tracking
US10663609B2 (en) * 2013-09-30 2020-05-26 Saudi Arabian Oil Company Combining multiple geophysical attributes using extended quantization
CN105676290B (zh) * 2016-04-03 2017-10-13 北京工业大学 基于曲面细分的地震数据三维显示方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4467461A (en) * 1981-01-05 1984-08-21 Conoco Inc. Interactive color analysis of geophysical data
JP3059495B2 (ja) * 1991-01-17 2000-07-04 株式会社リコー 三次元立体形状処理方式
US5332968A (en) * 1992-04-21 1994-07-26 University Of South Florida Magnetic resonance imaging color composites
US5782762A (en) * 1994-10-27 1998-07-21 Wake Forest University Method and system for producing interactive, three-dimensional renderings of selected body organs having hollow lumens to enable simulated movement through the lumen
US5563949A (en) * 1994-12-12 1996-10-08 Amoco Corporation Method of seismic signal processing and exploration
US5930730A (en) * 1994-12-12 1999-07-27 Amoco Corporation Method and apparatus for seismic signal processing and exploration
US5586082A (en) * 1995-03-02 1996-12-17 The Trustees Of Columbia University In The City Of New York Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging
US5798982A (en) * 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
US6144383A (en) * 1997-05-30 2000-11-07 Hewlett-Packard Company Volumetric data organization method that allows for cache efficient rendering speedups and efficient graphics hardware design
US5966672A (en) * 1997-07-28 1999-10-12 Knupp; Daniel F. Visualization technology method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100349003C (zh) * 2004-12-29 2007-11-14 中国石油天然气集团公司 一种消除起伏地表成像误差的方法
CN103185896A (zh) * 2006-09-01 2013-07-03 哈利伯顿兰德马克绘图公司 用于对波形体进行成像的系统和方法
CN103185896B (zh) * 2006-09-01 2016-08-10 哈利伯顿兰德马克绘图公司 用于对波形体进行成像的系统和方法
CN102160087A (zh) * 2007-01-05 2011-08-17 兰德马克绘图国际公司，哈里伯顿公司 用于实时可视化多卷数据集的方法及系统
CN101266299B (zh) * 2008-04-14 2011-03-30 林昌荣 利用地震数据体结构特征预测油气的方法
CN103221844A (zh) * 2010-09-01 2013-07-24 界标制图有限公司 绘制指示地下地层的特性的值的方法和系统
CN103221844B (zh) * 2010-09-01 2016-09-07 界标制图有限公司 绘制指示地下地层的特性的值的方法和系统
CN106056669A (zh) * 2016-06-07 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种海量地震数据的地质曲面重建方法
CN114333244A (zh) * 2021-12-31 2022-04-12 成都市美幻科技有限公司 多地震预警方法、装置及设备

Also Published As

Publication number Publication date
NO20032925L (no) 2003-09-15
NO20032925D0 (no) 2003-06-25
DE60121819D1 (de) 2006-09-07
DE60121819T2 (de) 2007-08-23
EP1358510A4 (en) 2005-07-27
CA2428743A1 (en) 2002-08-02
CA2428743C (en) 2009-05-12
CN100523869C (zh) 2009-08-05
BR0116765C1 (pt) 2004-07-06
WO2002061463A1 (en) 2002-08-08
EP1358510B1 (en) 2006-07-26
MXPA03004692A (es) 2005-07-01
BR0116765A (pt) 2004-02-25
CN1877364A (zh) 2006-12-13
AU2001234706B2 (en) 2006-11-09
NO337451B1 (no) 2016-04-18
EP1358510A1 (en) 2003-11-05
CN1214254C (zh) 2005-08-10

Similar Documents

Publication Publication Date Title
US6690820B2 (en) System and method for analyzing and imaging and enhanced three-dimensional volume data set using one or more attributes
US9171391B2 (en) Systems and methods for imaging a volume-of-interest
AU2008205064B2 (en) Systems and methods for selectively imaging objects in a display of multiple three-dimensional data-objects
US7702463B2 (en) Systems and methods for enhancing a seismic data image
US8022947B2 (en) Systems and methods for imaging waveform volumes
EP2102823B1 (en) Systems and methods for visualizing multiple volumetric data sets in real time
CN1214254C (zh) 用于分析和成像增强三维体数据集合的系统和方法
CN1860522A (zh) 实时共同显现多个属性的系统和方法
AU2001234706A1 (en) System and method for analyzing and imaging an enchanced three-dimensional volume data set using one or more attributes
CA2546458C (en) System and method for analyzing a region of interest relative to a predetermined event
EP1696247B1 (en) System and method for analysing and imaging an enhanced three-dimensional volume data set using one or more attributes

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
ASS Succession or assignment of patent right

Owner name: LANDLER MARK DRAWING CO.,LTD.

Free format text: FORMER OWNER: MAGIC EARTH INC.

Effective date: 20060908

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20060908

Address after: American Texas

Patentee after: Landmark Graphics Corp.

Address before: American Texas

Patentee before: Magic Earth Inc.

C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20050810

Termination date: 20100131