确定与感兴趣地下体积有关的地质信息的方法

摘要

一种确定与地下体积有关的地质信息的方法。在一个实施例中，所述方法包括：获得与地下体积有关的地震信息；确定自动地把地震信息映射到展平体积内的层位体积，其中，所述展平体积的一个轴对应于年代地层时间，并且其中，由地震信息表示的层位在层位体积内被自动解释，并且由层位体积移动到展平体积内，从而基本上是平面的并且基本上垂直于展平体积的与年代地层时间相对应的轴；确定层位体积相对于年代地层时间的导数；以及根据层位体积相对于年代地层时间的导数，来确定与地下体积有关的地质信息。

权利要求

1.一种计算机实现的确定与地下体积内的地层形成有关的信息的方法，所述方法包括：获得与地下体积有关的表示地下体积内的地质形成的地震信息，其中，所获得的地震信息的参数包括(i)表面平面内的二维位置，和(ii)地震时间；

分析所获得的地震信息，以自动识别地下体积内存在的由所获得的地震信息表示的一组层位；

从所述地震信息确定层位体积，其中，所述层位体积把所获得的地震信息映射到展平体积内，使得在展平体积内，每个所识别的层位被移动成与由层位体积限定为单个年代地层时间的估计的表面基本上共面，其中，展平体积的坐标是(i)表面平面内的二维位置，和(ii)与年代地层时间有关的度量；

在层位体积内，确定对于给定年代地层时间，地震时间相对于年代地层时间的导数，从而提供地层形成度量，所述地层形成度量与在展平体积内对应于给定年代地层时间的表面存在的地层在给定年代地层时间的形成速率有关；

从确定的导数来确定等时线体积；

根据确定的等时线体积来确定等容线体积；以及

从展平体积和确定的等容线体积，来确定等厚线体积。

2.按照权利要求1所述的方法，其中，在所述展平体积内，所述一组层位之间的间距与所述一组层位的沉积之间的年代地层时间有关。

3.按照权利要求1所述的方法，其中，对于多个年代地层时间，确定用于层位体积的地震时间相对于年代地层时间的导数。

4.按照权利要求1所述的方法，其中，通过把等时线体积内所包括的由层位体积指定的地震时间相对于与年代地层时间有关的度量的导数的值乘以地下体积的局部地震层速度的一半，来确定所述等容线体积。

5.按照权利要求1所述的方法，其中，通过把等容线体积内所包括的等容线值乘以展平体积内的点的局部倾斜的余弦，来确定所述等厚线体积。

6.按照权利要求1所述的方法，其中，给定年代地层时间对应于地下体积内的层位。

7.按照权利要求1所述的方法，进一步包括：把在给定年代地层时间的、用于层位体积的地震时间相对于年代地层时间的导数映射回地下体积的未映射坐标。

8.按照权利要求1所述的方法，进一步包括：实现在给定年代地层时间的、用于层位体积的地震时间相对于年代地层时间的导数，以推断与地下体积有关的地质信息。

9.按照权利要求1所述的方法，进一步包括：沿着展平体积内与给定年代地层时间相对应的表面，确定地震时间相对于年代地层时间的导数的梯度。

10.按照权利要求1所述的方法，进一步包括：生成向观察者图形传递地层形成度量的显示。

11.按照权利要求1所述的方法，其中，对于地下体积内的多个层位中的每个层位，沿该层位确定所述导数。

12.按照权利要求1所述的方法，其中，地震信息是从地震数据确定的，所述地震数据是基于接收从地下体积内形成的层位反射的地震脉冲而生成的。

13.一种计算机实现的确定与地下体积内的地层形成有关的信息的方法，所述方法包括：获得与地下体积有关的表示地下体积内的地质形成的地震信息，其中，所获得的地震信息的坐标是(i)表面平面内的二维位置，和(ii)地震深度；

分析地震信息，以自动识别地下体积内存在的由地震信息表示的一组层位；

从地震信息确定层位体积，其中，所述层位体积把获得的地震信息映射到展平体积内，使得在展平体积内，每个所识别的层位被移动成与由层位体积限定为单个年代地层时间的估计的表面基本上共面，其中，展平体积的坐标是(i)表面平面内的二维位置，和(ii)与年代地层时间有关的度量；

在层位体积内，对于给定年代地层时间，确定地震深度相对于年代地层时间的导数，从而提供地层形成度量，所述地层形成度量与在展平体积内对应于给定年代地层时间的表面存在的地层在给定年代地层时间的形成速率有关；

从确定的导数来确定等容线体积；以及

从展平体积和确定的等容线体积，来确定等厚线体积。

14.按照权利要求13所述的方法，其中，通过把等容线体积内所包括的等容线值乘以展平体积内的点的局部倾斜的余弦，来确定所述等厚线体积。

说明书

确定与感兴趣地下体积有关的地质信息的方法

[0001] 相关申请

[0002] 本申请涉及同此提交的题为“Method for indexing a SubsurfaceVolume for the Purpose of Inferring Geologic Information(用于索引地下体积以推测地质信息的方法)”的美国专利申请(代理人案卷号No.T7206)，其内容全部并入本公开。

技术领域

[0003] 本申请涉及用于确定与感兴趣地下体积有关的地质信息的系统和方法。具体地讲，本发明涉及确定与感兴趣地下体积内存在的地层的形成速度有关的度量，并且在一些情况下，涉及从所得到的度量推导其它地质信息。

背景技术

[0004] 用于确定与感兴趣地下体积内存在的地层形成的沉积速度有关的信息的技术已为人们所知。然而，这些技术通常需要与感兴趣地下体积有关的地震数据的人工分析(例如，地震数据内的层位的人工挑选)和/或不精确的计算，并通常导致被稀疏采样和/或低分辨率的数据。因而，实现所确定信息在其它地质信息的推导中的有用性被减小。

发明内容

[0005] 本发明的一个方面涉及一种计算机实现的确定与地下体积内的地层形成有关的信息的方法。在一个实施例中，所述方法包括：获得与地下体积有关的表示地下体积内的地质形成的地震信息，其中，获得的地震信息的参数包括(i)表面平面内的二维位置，和(ii)地震时间；分析所获得的地震信息，以自动识别地下体积内存在的由所获得的地震信息表示的一组层位；从地震信息确定层位体积，其中，所述层位体积把获得的地震信息映射到展平体积内，使得在展平体积内，每个所识别的层位被移动成与由层位体积限定为单个年代地层时间的估计的表面基本上共面，其中，所述展平体积的坐标是(i)所述表面平面中的二维位置，和(ii)与年代地层时间有关的度量；以及在层位体积内，确定对于给定年代地层时间，地震时间相对于年代地层时间的导数，从而提供地层形成度量，所述地层形成度量与在展平体积内对应于给定年代地层时间的表面存在的地层在给定年代地层时间的形成速率有关。

[0006] 本发明的另一方面涉及一种确定与地下体积有关的地质信息的方法。在一个实施例中，所述方法包括：获得与地下体积有关的地震信息，其中，地震信息表示地下体积内存在的层位；确定自动地把地震信息映射到展平体积内的层位体积，其中，展平体积的一个轴对应于年代地层时间，并且其中，由地震信息表示的层位在层位体积内被自动解释，并被层位体积移动到展平体积内，从而基本上是平面的，并且基本上垂直于展平体积的与年代地层时间相对应的轴；确定层位体积相对于年代地层时间的导数；以及根据层位体积相对于年代地层时间的导数，确定与地下体积有关的地质信息。

[0007] 本发明的另一方面涉及一种确定与地下体积有关的地质信息的方法。在一个实施例中，所述方法包括：获得地震数据，其中，所述地震数据是根据接收从在地下体积内形成的层位反射的地震脉冲而生成的；从地震数据确定地震信息，其中，所述地震信息描述了包含在地下体积内的层位的三维位置；确定把地震信息映射到展平体积内的层位体积，其中，所述展平体积的一个轴对应于年代地层时间，并且其中，层位体积把地震信息映射到展平体积内，使得在地震信息中描述的基本上所有层位都被单独移动成基本上是平面的，并且基本上垂直于展平体积的与年代地层时间相对应的轴；确定层位体积相对于年代地层时间的导数；以及根据层位体积相对于年代地层时间的导数，来确定与地下体积有关的地质信息。

[0008] 本发明的另一方面涉及一种计算机实现的确定与地下体积内的地层形成有关的信息的方法。在一个实施例中，所述方法包括：获得与地下体积有关的、表示地下体积内的地质形成的地震信息，其中，所获得的地震信息的坐标是(i)表面平面内的二维位置，和(ii)地震深度；分析所述地震信息，以自动识别存在于地下体积内的由地震信息表示的一组层位；从所述地震信息确定层位体积，其中，所述层位体积把获得的地震信息映射到展平体积内，使得在展平体积内，每个所识别的层位被移动成与由层位体积限定为单个年代地层时间的估计的表面基本上共面，其中，所述展平体积的坐标是(i)表面平面内的二维位置，和(ii)与年代地层时间有关的度量；以及在层位体积内，对于给定年代地层时间，确定地震深度相对于年代地层时间的导数，从而提供地层形成度量，所述地层形成度量与在展平体积内对应于给定年代地层时间的表面存在的地层在给定年代地层时间的形成速率有关。

[0009] 参考构成本说明书的一部分的附图，考虑下面的说明和附加权利要求，本发明的这些和其它目的、特征及特性，以及结构的相关元件的操作方法和功能，部件的组合和制造的经济性将变得更明显，在附图中，相同的附图标记表示各个图中对应的部件。然而，附图只是用于举例和说明，而不是对本发明的限制。说明书和权利要求中使用的单数形式“一个”包括复数的所指事物，除非上下文明确地另有说明。

附图说明

[0010] 图1示例了按照本发明的一个实施例，与感兴趣地下体积有关的地震信息的图形表示，该图形表示表现了感兴趣地下体积内的地质形成。

[0011] 图2示例了按照本发明的一个实施例，与感兴趣地下体积相对应的展平体积的图形表示的正视图。

[0012] 图3示例了按照本发明的一个实施例，确定感兴趣地下体积的层位体积，并处理层位体积以得到与感兴趣地下体积有关的地质信息的方法。

具体实施方式

[0013] 反射地震学(或者说地震反射)是一种使用地震学的原理，根据被导向感兴趣地下体积并从其反射的地震波，来估计地球的地下的感兴趣地下体积的性质的勘探地球物理学方法。该方法通常需要地震波的源，比如爆炸物(例如，炸药/Tovex)、专门的气枪、震动源和/或其它源。所述源被用于在感兴趣地下体积内引入地震波(通常在地表)，并借助(通常置于地表的)地震计阵列获得数据，所述地震计探测当源生成的地震波到达地下体积的表面时的反射。随后处理(例如，移动等)由地震计获得的数据以形成地震数据立方体，所述地震数据立方体在地下体积内的数据点阵列表示地下体积内存在的地质形成。由于地震数据立方体内的数据点一般是高度空间采样的，因此地震数据立方体实际上是感兴趣地下体积内存在的地层的图像。应意识到，在一些情况下，地震数据能够包括比这里讨论的三维地震数据立方体的维数少的地震数据(例如一维或二维地震数据)，和/或它能够包括时间维度以及空间维度。

[0014] 举例来说，图1示例了与感兴趣地下体积12有关的地震信息的图形表示的正视图，它表示感兴趣地下体积12内的地质形成(例如，从对感兴趣地下体积12进行反射地震学而获得的地震数据立方体)。虽然图1提供的图形表示把地震信息描述成是连续的，但地震信息一般在感兴趣地下体积12内的数据点阵列表示感兴趣地下体积12内的地质形成。然而，由于地震信息的高度空间采样(例如，数据点的接近)，为了便于举例说明，信息可被表示成是连续的。此外，尽管图1只提供了地震信息的图形表示的二维视图，但这只是为了便于举例说明，应意识到，可以提供具有进入或离开图1的平面的另一维的三维视图。例如，在一个实施例中，表现感兴趣地下体积12内的地质形成的地震信息的参数是感兴趣地下体积12的表面平面内的二维位置(例如，x轴14和与图1的平面垂直的y轴16)、以及与感兴趣地下体积12的地震深度有关的度量(例如，t轴18，其中t表示地震时间)。与感兴趣地下体积12的地震深度有关的度量可以是地震时间或地震深度。地震时间与地震波在感兴趣地下体积12内的点和地表之间传播所用的时间有关。

[0015] 从图1可看出，感兴趣地下体积12包括多个层位20。层位是在感兴趣地下体积12的地层内的两层不同组成之间的边界处形成的表面。由于层位20表示感兴趣地下体积
12的地层的组成的边界变化，因此为了分析目的，假定每个层位20表示感兴趣地下体积12内在共同年代地层时间沉积的表面。这是合理的，因为对于沉积在感兴趣地下体积12的一个区域中的地层的组成变化的推动力很可能是沉积在感兴趣地下体积12的另一个区域中的地层的类似组成变化的推动力。特别是在指示两个区域的这种组成变化的边界相连接，和/或在感兴趣地下体积12内具有相似深度的情况下更是如此。

[0016] 在层位(或者感兴趣地下体积12内与单个年代地层时间相对应的某个其它表面)上的给定点处，可用其位置(例如，x，y，t坐标)及其“倾斜”来描述层位(或其它表面)。“倾斜”是层位相对于水平面(或者某个其它平面，通常垂直于“走向”)的方位的测量。例如，在图1中，层位20中的一个给定层位的倾斜可用两个分量表征，即，x-倾斜(或者给定层位20相对于x轴14的角度)和y-倾斜(或者给定层位20相对于y轴16的角度)。

[0017] 如图1中所示，感兴趣地下体积12内的层位20一般具有起伏和/或不连续面。这些起伏和/或不连续面中的一些可能是由与地层沉积相关联的各种现象造成的。例如，这些现象包括：可变沉淀速率、可变侵蚀速率、沉积时的表面起伏、可变岩石组成，和/或其它现象。在一些情况下，所述起伏和/或不连续面可能是影响层位20的形状的各种地下作用力的结果。例如，这些现象可包括：地震活动、火山活动、地下水流动、地下沉积物流动(例如，盐运动)、不均匀负载、可变压实速率、可变成岩固结速率、构造驱动的变形(例如，褶皱和/或断层)，和/或其它现象。

[0018] 通过分析表示感兴趣地下体积12内的地质形成(例如，层位20、层位20之间的多层地层等)的地震信息，能够确定与感兴趣地下体积12有关的地质信息。用于分析这种信息的一种技术通常被称为“展平”。例如，图2示例了与感兴趣地下体积12相对应的展平体积22的图形表示的正视图。在展平体积22内，在感兴趣地下体积12内存在的一些或全部层位20被整形以对应于某个共同形状。所述共同形状与由层位体积限定为沉积的单个年代地层时间的估计的表面共面。在图2中所示的展平体积22内，层位20已被整形成与平坦且水平的表面共面。

[0019] 由于每个层位20被认为表示感兴趣地下体积12内在共同年代地层时间沉积的表面，因此展平体积22的纵轴(图2中所示，并且在下文中被称为年代地层时间轴24)变成与沉积的年代地层时间而不是与地震深度有关的度量。更具体地说，展平体积22的参数是展平体积22的表面平面中的二维位置以及与涉及沉积的年代地层时间的信息有关的度量。

[0020] 尽管在图2中，层位20已被展平为基本上与平坦并且垂直于年代地层时间轴24的表面共面，但这不应被看作是对本发明的限制。在一些情况下，层位20可被展平为基本上与不平坦的表面和/或不垂直于年代地层时间轴24的表面中的一个或两者共面(例如，被展平成与感兴趣地下体积12内的层位20之一的形状相对应的表面)。

[0021] 通过查看获得的展平体积22形式的地震信息，解释者(例如，用户，计算机等)能够像感兴趣地下体积12内的地质特征原先被沉积那样地查看所述地质特征。例如，在展平之后，解释者能够在单个图像中查看具有弯曲河道的整个洪泛平原区的保留特征。

[0022] 例如，在下述文献中说明了用于展平获得的地震信息的一种技术：“Flattening without picking”，Lomask 等 人，StandfordExploration Project，Report 112，2002年11月11日，第141-151 页；“Flattening 3-D data cubes in complex geology”，Lomask，Standford Exploration Project，Report 113，2003年7月8日，第247-261页；“Flattening without picking”，Lomask等人，Geophysics，第71卷，第4期(2006年7月-8月)，第13-20页；以及“Volumetricflattening：an interpretation tool”，Lomask等人，The Leading Edge，2007年7月，第888-897页(在下文中被总称为“Lomask”)。这些出版物全部并入本公开文献中。如在Lomask中所述，当把其中描述的展平技术应用于获得的地震信息时，结果包括“层位体积”。这里使用的层位体积是把获得的地震信息映射到展平体积22内或相反的数据集。从而，在一个实施例中，层位体积可被概念化成一个函数，该函数根据其在展平体积22内的位置(例如，(x，y，τ))，来提供地震信息内与展平体积
22内的给定位置相对应的位置(例如，(x，y，t))。由于表面平面中的二维位置(例如，图
1和2中的x，y坐标)在感兴趣地下体积12和展平体积22之间并不改变，因此，对于展平体积22内的给定位置(x，y，τ)指定t使得来自指定t的信息能够将感兴趣地下体积12内的给定位置的相同二维表面位置(例如，相同的x和y)在展平体积22内的给定位置被映射到展平体积22内。

[0023] 如前所述，地震信息通常将在感兴趣地下体积12内的数据点阵列表示感兴趣地下体积12内的地质形成。因此，概念化层位体积的另一种方式是作为包括在展平体积22内的所获得地震信息中的数据点的索引，其中，所述索引指定与感兴趣地下体积12有关的地震信息中应被映射到展平体积22的给定位置的数据点的位置。

[0024] 图3示例了确定层位体积并处理层位体积以提取与感兴趣地下体积有关的地质信息的方法26。在下面的方法26及其一个或多个操作的说明中，具体参考了在Lomask中说明的展平技术。然而，这不应被视为是对本发明的限制。相反，方法26应被理解为能够与各种不同的展平技术一起使用。此外，图3中示例的和在下文描述的方法26的操作的特定安排并不是限制性的。应理解，在一些实现中，可以按照除所述顺序之外的顺序进行各个操作(或者可伴随其它操作地进行各个操作)，各个操作可以与其它操作相结合和/或完全被省略，和/或可以增加另外的各种操作，而不脱离本发明的范围。

[0025] 在方法26的操作28，获得与感兴趣地下体积有关的地震信息，所述地震信息表示感兴趣地下体积内的地质形成。在一个实施例中，获得的地震信息的参数包括表面平面中的二维位置和与地震深度有关的度量。在一些情况下，在操作28获得的地震信息在地下体积内的数据点阵列表示在感兴趣地下体积内存在的地质形成。在一个实施例中，在操作28获得的地震信息包括表示在感兴趣地下体积内存在的地质形成的三维位置的地震数据立方体。

[0026] 在操作30，分析在操作28获得的地震信息以识别在感兴趣地下体积内存在的层位。在一个实施例中，在操作30识别层位包括：通过将在操作28获得的地震信息中在所识别层位的给定层位上的每个数据点与该给定层位相关联，来解释在获得的地震信息内表示的各个层位。例如，可按照在Lomask中描述的技术自动地解释层位。

[0027] 在一个实施例中，操作30识别存在于感兴趣地下体积内的相对于地震深度/时间(例如，沿图1中所示并且在上面说明的t-轴18)的密度大于预定密度阈值的层位。预定密度阈值可由用户根据与感兴趣地下体积有关的一个或多个参数，和/或根据获得的地震信息的参数来配置。预定密度阈值的例子可以是对于地震时间体积内的每个地震周期至少一个层位。备选层位解释策略可以是把地震道的每个最小值或最大值解释为层位。

[0028] 在操作32，确定层位体积，层位体积把在操作28获得的地震信息映射到展平体积内，使得在展平体积内，在操作30识别的每个层位按照沉积的年代地层时间排序，并被移动成与由层位体积限定为单个年代地层时间的估计的表面基本上共面。因而，展平体积的坐标是表面平面中的二维位置和与沉积的年代地层时间有关的度量。在一个实施例中，通过指定与展平体积内的位置相对应的地震时间/深度值，层位体积把在操作28获得的地震信息映射到展平体积内。

[0029] 在一个实施例中，与沉积的年代地层时间有关的度量是年代地层时间的相对测量。情况通常是这样，因为层位体积按照沉积的年代地层时间对在操作30识别的层位进行排序，但并不指定沉积的年代地层时间的实际值。然而，在一些情况下，用户可以为沉积的年代地层时间指定更绝对的标度，以便应用于最初由层位体积提供的相对度量。

[0030] 在一个实施例中，由层位体积限定为单个年代地层时间的估计的表面具有平面形状，并且与和年代地层时间有关的度量相对应的轴正交。在一个实施例中，由层位体积限定为单个年代地层时间的估计的表面具有与在操作28识别的一个或多个层位的形状相对应的形状。

[0031] 在操作34，按照层位体积，把在操作28获得的地震信息映射到展平空间内。如图2中所示和上面所述，这导致在操作30识别的每个层位被移动成与由层位体积限定为沉积的单个年代地层时间的估计的表面基本上共面。

[0032] 在操作36，取得在操作32确定的层位体积相对于与年代地层时间有关的度量的导数。例如，当用地震深度来参数化在操作28获得的地震信息时，在操作36取得的导数是由层位体积指定的地震深度相对于年代地层时间的导数。作为另一个例子，当用地震时间来参数化在操作28获得的地震信息时，在操作36取得的导数是由层位体积指定的地震时间相对于与年代地层时间有关的度量的导数。在一个实施例中，沿着与沉积的单个年代地层时间相对应的表面，在层位体积内取得在操作36相对于与年代地层时间有关的度量取得的导数。这种表面的一个例子是在操作30识别的层位。然而，在层位体积内也存在与单个年代地层时间相对应的其它表面。

[0033] 在操作38，确定是否已经沿着层位体积内与单个年代地层时间相对应的一组预定表面的每个表面，取得了层位体积相对于与年代地层时间有关的度量的导数。在一个实施例中，所述一组预定表面至少包括在操作30识别的层位。在一个实施例中，所述一组预定表面不包括在操作30识别的一个或多个层位。表面相对于年代地层时间的频率、表面的数量、表面相对于地震深度/时间的频率、和/或所述一组表面的其它参数中的一个或多个可由用户配置，和/或可被自动确定(例如，根据在操作30的层位的识别、根据在操作28获得的地震信息等)。

[0034] 如果在操作38确定还未在层位体积内的一组预定表面中的每个表面确定层位体积相对于与年代地层时间有关的度量的导数，那么方法26返回操作36，并且沿着包括在所述一组预定表面中的另一个表面，取得层位体积相对于与年代地层时间有关的度量的导数。如果在操作38确定已经沿着所述一组预定表面中的每个表面确定了层位体积相对于与年代地层时间有关的度量的导数，那么方法26前进至操作40。

[0035] 在操作36对于层位体积内的给定点取得的导数可被概念化成描述在地下体积内与层位体积内的给定点相对应的点处地下体积变得“更厚”(例如，地下体积沿图1的t轴18增大)的速率。在感兴趣地下体积的形成过程中，当沉积/形成另外的地层时，感兴趣地下体积变得“更厚”。因而，在操作36取得的导数提供了与感兴趣地下体积内存在的地层的形成速率有关的地层形成度量。地层形成度量涉及，例如，沉降速率、同沉积生长速率、一个或多个成岩过程、地震速度、和/或与感兴趣地下体积内的地层的形成有关的其它因素。

[0036] 在操作40，从取得的层位体积相对于与年代地层时间有关的度量的导数，来确定与感兴趣地下体积有关的地质信息。可自动地(例如，按照某种预定算法)、人工地(例如，通过用户分析在操作36确定的一个或多个导数)，和/或借助自动和/或人工分析的某种组合来进行操作40。在一个实施例中，地质信息包括感兴趣地下体积的一个或多个储层特性(例如，孔隙率、净总值比、页岩体积百分数、渗透率、岩性、沉积相等)。从操作42、44、46、48、50、52和/或54的描述应明白，方法26可包括其它操作，以便于从在操作32取得的层位体积的导数来确定地质信息。

[0037] 在一个实施例中，方法26包括操作42，在操作42，确定层位体积相对于与年代地层时间有关的度量的导数的梯度。在一些情况下，在操作42确定的梯度是沿着层位体积内与沉积的各个年代地层时间相对应的表面针对该导数的。在一些情况下，在三维区域(例如，展平体积内的三维区域、感兴趣地下体积内的三维区域等)上，确定在操作42确定的层位体积相对于与年代地层时间有关的度量的导数的梯度。在一些情况下，可垂直于或平行于感兴趣的特定地质形成，来确定在操作42确定的梯度。在一些情况下，在操作42确定的梯度可包括沿着从感兴趣的特定地质形成向外扩展的径向线投射的梯度。感兴趣的地质形成、和/或其边界可由用户人工识别，和/或被自动识别。在操作42的一个或多个梯度的确定可便于确定与古河道梯度、古流向和流速、古地形有关的地质信息，和/或其它地质信息。

[0038] 在一个实施例中，方法26包括操作44，在操作44，取得层位体积相对于与沉积的年代地层时间有关的度量的高阶导数(高于在操作36取得的一阶导数)。在一个实施例中，操作44包括：沿着层位体积内与各个年代地层时间相对应的一个或多个表面取得所述高阶导数。在操作44的高阶导数的确定可便于确定与更高沉降周期的开始和终止的定时有关的地质信息。

[0039] 在一个实施例中，方法26包括操作46，在操作46，从在操作36取得的导数来确定等时线体积。例如，在用地震时间(例如，而不是地震深度)参数化在操作28获得的地震信息的一些情况下，通过将在操作36为其确定层位体积的导数的表面上的各个位置与在操作36确定的层位体积相对于与年代地层时间有关的度量的导数的对应值相关联，来确定等时线体积。因而，在操作46确定的等时线体积把由层位体积指定的地震时间相对于与年代地层时间有关的度量的导数描述成展平体积内的位置的函数。

[0040] 在一个实施例中，方法26包括操作48，在操作48，从在操作36取得的导数来确定等容线体积。例如，在用地震深度(例如，而不是地震时间)参数化在操作28获得的地震信息的一些情况下，通过将在操作36为其确定层位体积的导数的表面上的各个位置与在操作36确定的层位体积相对于与年代地层时间有关的度量的导数的对应值相关联，来确定等容线体积。因而，在操作48确定的等容线体积把由层位体积指定的地震深度相对于与年代地层时间有关的度量的导数描述成展平体积内的位置的函数。作为另一个例子，在用地震时间(例如，而不是地震深度)参数化在操作28获得的地震信息的一些情况下，根据在操作46确定的等时线体积来确定等容线体积。例如，可通过把在操作46确定的等时线体积内所包括的由层位体积指定的地震时间相对于与年代地层时间有关的度量的导数的值乘以地下体积的局部地震层速度的一半(或者，在与年代地层时间有关的度量是相对的非标度的年代地层时间的情况下，乘以局部地震层速度的任意标量倍数)，来确定等容线体积。

[0041] 在一个实施例中，方法26包括操作50，在操作50，从在操作34得到的展平体积和在操作48得到的等容线体积，来确定等厚线体积。为了确定等厚线体积，把在操作48得到的等容线体积内所包括的等容线值乘以局部倾斜的余弦。如上所述，展平体积内的给定层位上的点的局部倾斜(实现以便确定等容线体积)描述了给定层位相对于感兴趣未展平地下体积内的水平面的方位。用于确定地下体积内的局部倾斜的技术已相当成熟。因此，在一个实施例中，从在操作28获得的地震信息，来计算展平体积内的点的局部倾斜，随后将其映射到展平体积内，使得在操作48确定的等容线值能够被乘以局部倾斜的余弦，以确定等厚线体积。

[0042] 在一个实施例中，方法26包括操作52，在操作52，按照层位体积，把在操作46得到的等时线体积、在操作48得到的等容线体积、和/或在操作50得到的等厚线体积内的一个或多个映射回在操作28获得的地震信息的参数(例如，映射回感兴趣未展平地下体积内)。因而，操作52提供分别包括在操作46、48或50得到的等时线、等容线或等厚线值的一个或多个体积，并表示在操作30识别的层位的三维形状和位置。

[0043] 在一个实施例中，方法26包括操作54，在操作54，向用户显示信息(例如，借助电子显示器)。在操作54显示的信息可包括在操作32、34、36、42、44、46、48、50和/或52中的一个或多个操作得到的信息。因而，在操作54显示的信息可以使用户能够执行和/或监视或控制性能操作40。在一些情况下，操作54可包括显示在操作40确定的地质信息(不同于或者除了在操作32、34、36、42、44、46、48、50和/或52中的一个或多个操作确定的地质信息之外)。在这种情况下，操作54实现了操作40的结果(或者部分结果)的显示。

[0044] 尽管基于目前认为最实际并且优选的实施例详细举例说明了本发明，但应理解这样的细节只是用于举例说明，并且本发明并不局限于所公开的实施例，相反，本发明意图覆盖在附加权利要求的精神和范围内的修改和等同安排。例如，应理解，在可能的程度，本发明预期任意实施例的一个或多个特征能够与任何其它实施例的一个或多个实施例相组合。