一种等时地震切片获取方法、储层预测方法和装置

技术领域

本发明涉及石油勘探与开发技术领域，特别涉及一种等时地震切片获取方法、储层预测方法和装置。

背景技术

近年来，我国每年新增探明储量中，岩性地层油气藏储量均超过50％，且呈逐年递增趋势。无论在我国东、西部的大型坳陷盆地(松辽、鄂尔多斯盆地)和断陷型盆地(渤海湾盆地)，还是在西部的前陆盆地(准噶尔盆地)，都发现了类型多样、储量丰富的岩性地层油气藏。然而这些油气藏包含的油气层单层厚度一般都小于地震纵向分辨率，这给地震预测带来了巨大的困难。

针对地震资料分辨率的技术瓶颈，许多学者从提高垂向分辨率角度开展反演算法研究，进而预测薄储层的厚度与分布，并形成了稀疏反射系数反演、谱反演、波形指示反演和相控随机反演等多种反演方法。还有学者以调谐理论为基础，研究了薄层频谱、薄层分辨率、薄层峰值频率与厚度的关系、薄互层调谐与分辨率的关系等，并形成了谱分解、峰值瞬时频率属性、频率与振幅比属性等薄层预测方法。曾洪流、Posamentier等则从利用地震资料的横向分辨率角度分别提出了地震沉积学和地震地貌学的概念，它的核心思想是通过地层切片扫描来检测薄层，所谓的地层切片就是在2个等时标志层间通过等比例内插得到的切片，从其原理可以看出，地层切片主要适合于研究区内横向沉积速率保持稳定的地区。然而，陆相盆地不同区域沉积速率往往变化较大，传统的地层切片容易穿时，从而导致薄砂体预测失败。

发明内容

鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的沥一种等时地震切片获取方法、储层预测方法和装置，能够合理从地震数据体中切取等时地震切片，为实现储层的合理预测提供了数据依据。

第一方面，本发明实施例提供一种等时地震切片获取方法，包括：

获取地震数据体目的层切片中基准点和参考点的位置信息与漂移量，所述基准点的漂移量为0，所述参考点的漂移量不等于0；

针对地震数据体目的层切片中除所述基准点和参考点外的任一点，根据该点的位置信息及所述基准点和参考点的位置信息与漂移量，确定该点的漂移量；

将所述目的层切片中每一点按该点的漂移量在所述地震数据体中沿深度方向进行漂移，得到等时地震切片。

第二方面，本发明实施例提供一种储层预测方法，包括：

从地震数据体中提取设定属性的第二地震属性体，从所述第二地震属性体中提取等时地震切片处的第二地震属性切片，所述等时地震切片是按照上述方法从所述地震数据体中切取的；

根据井目的层储层特征对所述第二地震属性切片中与该井对应的点进行储层特征标记；

根据所述第二地震属性切片和储层特征标记，进行储层特征预测。

第三方面，本发明实施例提供一种等时地震切片获取装置，包括：

参考数据获取模块，用于获取地震数据体目的层切片中基准点和参考点的位置信息与漂移量，所述基准点的漂移量为0，所述参考点的漂移量不等于0；

漂移量确定模块，用于针对地震数据体目的层切片中除所述基准点和参考点外的任一点，根据该点的位置信息及所述基准点和参考点的位置信息与漂移量，确定该点的漂移量；

漂移模块，用于将所述目的层切片中每一点按该点的漂移量在所述地震数据体中沿深度方向进行漂移，得到等时地震切片。

第四方面，本发明实施例提供一种储层预测装置，包括：

提取模块，用于从地震数据体中提取设定属性的第二地震属性体，从所述第二地震属性体中提取等时地震切片处的第二地震属性切片，所述等时地震切片是按照权上述方法从所述地震数据体中切取的；

标记模块，用于根据井目的层储层特征对所述第二地震属性切片中与该井对应的点进行储层特征标记；

预测模块，用于根据所述第二地震属性切片和储层特征标记，进行储层特征预测。

第五方面，本发明实施例提供一种具备储层预测功能的计算机程序产品，包括计算机程序/指令，其中，该计算机程序/指令被处理器执行时实现上述等时地震切片获取方法，或实现上述储层预测方法。

本发明实施例提供的上述技术方案的有益效果至少包括：

本发明实施例提供的等时地震切片获取方法，针对地震数据体目的层切片中除所述基准点和参考点外的任一点，根据该点的位置信息及基准点和参考点的位置信息与漂移量，确定该点的漂移量；将目的层切片中每一点按该点的漂移量在地震数据体中沿深度方向进行漂移，得到等时地震切片。以基准点的时间为基准，同时参考参考点的漂移量来合理确定目的层切片中其他点的漂移量，从而对目的层切片中每个点进行合理漂移，在地震数据体中切取等时地震切片，大大降低了等时地震切片的穿时性，使得通过地震数据进行薄层储层预测成为可能。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：

图1为提取的地震子波波形示意图；

图2为本发明实施例一中等时地震切片获取方法的流程图；

图3为本发明实施例一中基准点确定方法的具体实现流程图；

图4为本发明实施例一中参考点确定方法的具体实现流程图；

图5为本发明实施例中等时地震切片获取方法的具体实现流程图；

图6A为地震剖面示例图；

图6B为常规切片方法得到的目的层地震振幅属性平面图；

图6C为实施例一方法得到的目的层地震振幅属性平面图；

图7为古地形图的示意图；

图8为本发明实施例中等时地震切片获取装置的结构示意图；

图9为本发明实施例中储层预测装置的结构示意图。

具体实施方式

下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。

应理解本发明中所述的术语仅仅是为描述特别的实施方式，并非用于限制本发明。另外，对于本发明中的数值范围，应理解为还具体公开了该范围的上限和下限之间的每个中间值。在任何陈述值或陈述范围内的中间值以及任何其他陈述值或在所述范围内的中间值之间的每个较小的范围也包括在本发明内。这些较小范围的上限和下限可独立地包括或排除在范围内。

除非另有说明，否则本文使用的所有技术和科学术语具有本发明所属领域的常规技术人员通常理解的相同含义。虽然本发明仅描述了优选的方法和材料，但是在本发明的实施或测试中也可以使用与本文所述相似或等同的任何方法和材料。本说明书中提到的所有文献通过引用并入，用以公开和描述与所述文献相关的方法和/或材料。在与任何并入的文献冲突时，以本说明书的内容为准。

在本发明的描述中，需要说明的是，术语“包含”、“包括”、“具有”、“含有”等等，均为开放性的用语，即意指包含但不限于。此外，术语“第一”、“第二”和“第三”等仅用于描述目的，而不能理解为指示或暗示相对重要性。

为了解决现有技术中存在的地震数据体切片容易穿时导致薄砂体预测失败的问题，本发明实施例提供了一种等时地震切片获取方法、储层预测方法和装置，能够合理从地震数据体中切取等时地震切片，为实现储层的合理预测提供了数据依据。

本发明实施例所需的基础数据及其处理包括：

1、原始叠后地震数据体。

2、建立时间和深度对应关系。

对研究区内具有测井曲线的井制作合成记录，将合成记录与井旁叠后地震数据道进行对比，建立地震反射时间与深度对应关系。

3、原始叠后地震数据体的相位转换。

可以是通过希尔伯特变换将原始叠后地震数据体转换成多个相位的叠后地震数据体。通过时间和深度对应关系将测井岩性指示曲线(如自然电位、自然伽马等)和叠后地震数据体相匹配，选定一种匹配吻合度较高的叠后地震数据体。

申请人经过多次实验对比，发现-90度相位叠后地震数据体与测井岩性曲线的吻合度较为，故可以选用-90度相位叠后地震数据体作为本实施例所用的地震数据体。

例如，以某研究区目的层的等时地震切片获取为例，根据测井资料和井旁道提取的子波(图1)可知，原始叠后地震数据体接近零相位，极性为正极性，子波长度约70ms。其与典型井(wellA)的合成地震记录可知，测井曲线显示目标层深度位于3540米左右，对应地震数据时间为3810毫秒。通过希尔伯特变换将原始叠后地震数据体转换成-90度相位叠后地震数据体。

4、上部参考层和下部参考层。

选择等时的标志地震同相轴作为参考层。在叠后地震数据体中根据测井标定结果选择地质时间相同的标志地震同相轴作为参考层，通常情况下，可选择最大洪(湖)泛面或横向稳定分布的煤层等反射连续的层。最终选定目的层的上部参考层和下部参考层。

5、叠后地震数据体的目的层切片。

根据井的合成记录、目的层深度值和井轨迹数据，在地震数据体中确定目的层的标注点，根据多口井得到多个标注点；在上部参考层和下部参考层的约束下，根据多个标注点，通过厚度等比例内插方法在叠后地震数据体中切取目的层切片。

6、设定属性的地震属性体提取。

从叠后地震数据体中提取设定属性的地震属性体。

优选的，上述设定属性可以为地震振幅；也可以是频率或相位等其它属性。可以是一种地震属性的地震属性体，也可以是多种地震属性的多个地震属性体。

7、生成目的层沉积前的古地形图。

确定目的层切片与参考层间的厚度分布数据，参考上部参考层或下部参考层，得到目的层沉积前的古地形图。这里的厚度，可以是时间厚度，也可以是实际厚度。

实施例一

本发明实施例一提供一种等时地震切片获取方法，其流程如图2所示，包括如下步骤：

步骤S21：获取地震数据体目的层切片中基准点和参考点的位置信息与漂移量。

上述基准点的漂移量为0，参考点的漂移量不等于0。

在一些实施例中，参见图3所示，通过下述步骤确定基准点：

步骤S2111：从地震数据体中提取设定属性的第一地震属性体，从第一地震属性体中提取目的层切片处的第一地震属性切片。

设定属性可以是地震振幅、频率和相位中的一种或多种。

步骤S2112：根据井目的层储层特征和对应的标注点，对第一地震属性切片中与该井的标注点对应的点进行储层特征标记。

储层特征可以是岩相特征，也可以是沉积微相特征，也可以是其他特征。

步骤S2113：根据储层特征标记，从第一地震属性切片中确定基准点。

通常选择储层特征明显的点作为基准点。

在一些实施例中，参见图4所示，通过下述步骤确定参考点：

步骤S2121：确定目的层切片与参考层间的厚度分布数据。

参数层为上部参考层或下部参考层。

步骤S2122：根据厚度分布数据从目的层切片中确定参考点和参考点的漂移量。

根据厚度分布数据得到的厚度等值线图可以代表目的层沉积前的古地形图，以基准点处古地形厚度为基准，若参考点处古地形厚度大于基准点处古地形厚度，则参考点处向下漂移，反之，则参考点处向上漂移；根据参考点处古地形厚度与基准点处古地形厚度的差值大小，确定参考点的漂移幅度。根据参考点的漂移方向和漂移幅度确定其漂移量。

在目的层切片基础上，通过古地形指示切片移动方向，通过地震地貌特征的变化确定最佳漂移量，从而得到等时切片并避免穿时性。

步骤S22：针对地震数据体目的层切片中除基准点和参考点外的任一点，根据该点的位置信息及基准点和参考点的位置信息与漂移量，确定该点的漂移量。

具体的，可以根据该点的深度值及基准点和参考点的深度值与漂移量，确定该点的漂移量。

深度值可以为时间深度值，或实际深度值，或相对深度值。

在一些实施例中，可以根据下式确定点的漂移量：

其中，AutoShift为点的漂移量，T

步骤S23：将目的层切片中每一点按该点的漂移量在地震数据体中沿深度方向进行漂移，得到等时地震切片。

以基准点处古地形厚度为基准，其它地区如果古地形厚度大于基准点处古地形厚度，则该点处目的层切片向下漂移，反之，则目的层切片向上漂移；随古地形厚度的不同，其漂移量也不同。最终得到等时的地震切片。

本发明实施例一提供的等时地震切片获取方法，针对地震数据体目的层切片中除所述基准点和参考点外的任一点，根据该点的位置信息及基准点和参考点的位置信息与漂移量，确定该点的漂移量；将目的层切片中每一点按该点的漂移量在地震数据体中沿深度方向进行漂移，得到等时地震切片。以基准点的时间为基准，同时参考参考点的漂移量来合理确定目的层切片中其他点的漂移量，从而对目的层切片中每个点进行合理漂移，在地震数据体中切取等时地震切片，大大降低了等时地震切片的穿时性，使得通过地震数据进行薄层储层预测成为可能。

在一些实施例中，参见图5所示，上述等时地震切片获取方法，可以包括：

1、目的层沉积前古地形的生成；2、等时切片标准点的提取；3、等时切片的生成。其中目的层沉积前古地形的生成1包括：建立时间和深度对应关系101、生成-90度相位地震数据体102、选择等时的标志地震同相轴作为参考层103、生成目的层地层切片104、生成目的层沉积前的古地形图105。等时切片标准点的提取2包括：提取目的层地层切片处的地震属性值201、选取等时切片的基准点201。等时切片的生成3包括：地层切片漂移方向的确定301、确定地层切片漂移量，生成等时切片302。

以某研究区目的层的等时地震切片获取为例，参见图6A～6C所示，图6A为图6B中剖面位置的地震剖面图，图6A中的t1和t2分别是目的层的上部参考层和下部参考层，在目标层的上下选取等时的标志地震同相轴t1、t2作为参考层，并在t1、t2约束下提取目的层切片s1，进一步提取s1对应的地震振幅属性，参见图6B所示。图6B中根据wellA标注的地质特征与地震属性吻合，证明提取的地震振幅属性存在一定的可信度。在图6B中东部可见一条规模较大的呈北东-南西走向的河道，河道的中部和南部显示清晰，但在北部显示模糊，存在河道模糊区。在地震振幅属性上选取河道显示清晰的中部点B作为等时切片标准点。

图7是根据目的层切片s1和上部参考层t1制作的该期沉积时的古地形图，根据古地形图指示切片移动方向，对原始目的层切片作漂移处理，并提取相应的地震振幅值，图6C是漂移处理后，也就是等时切片对应的振幅属性图，由图6C可见北部河道显示清晰。

可见上述方法解决了由于陆相沉积地层不同区域沉积速率差异大，传统的地层切片技术容易穿时，导致薄层预测结果不准的问题，提供了一种新的薄层预测技术，是地震沉积学在薄互层预测方面的重大突破。

基于本发明的发明构思，本发明实施例还提供一种储层预测方法，包括：

从地震数据体中提取设定属性的第二地震属性体，从第二地震属性体中提取等时地震切片处的第二地震属性切片，等时地震切片是按照上述方法从地震数据体中切取的；根据井目的层储层特征对第二地震属性切片中与该井对应的点进行储层特征标记；根据第二地震属性切片和储层特征标记，进行储层特征预测。

上述储层特征可以是岩相特征，也可以是沉积微相特征，通过上述方法例如可以实现薄层砂体的储层预测。

基于本发明的发明构思，本发明实施例还提供一种等时地震切片获取装置，该装置的结构如图8所示，包括：

参考数据获取模块81，用于获取地震数据体目的层切片中基准点和参考点的位置信息与漂移量，所述基准点的漂移量为0，所述参考点的漂移量不等于0；

漂移量确定模块82，用于针对地震数据体目的层切片中除所述基准点和参考点外的任一点，根据该点的位置信息及所述基准点和参考点的位置信息与漂移量，确定该点的漂移量；

漂移模块83，用于将所述目的层切片中每一点按该点的漂移量在所述地震数据体中沿深度方向进行漂移，得到等时地震切片。

在一些实施例中，漂移量确定模块82，所述根据该点的位置信息及所述基准点和参考点的位置信息与漂移量，确定该点的漂移量，具体用于：

根据该点的深度值及所述基准点和参考点的深度值与漂移量，确定该点的漂移量。

在一些实施例中，漂移量确定模块82，所述根据该点的深度值及所述基准点和参考点的深度值与漂移量，确定该点的漂移量，具体用于根据下式确定该点的漂移量：

其中，AutoShift为该点的漂移量，T

在一些实施例中，参考数据获取模块81，还用于通过下述步骤获取所述目的层切片：

根据井的合成记录、目的层深度值和井轨迹数据，在所述地震数据体中确定目的层的标注点，根据多口井得到多个标注点；

在上部参考层和下部参考层的约束下，根据所述多个标注点，通过厚度等比例内插方法在所述地震数据体中切取目的层切片。

在一些实施例中，参考数据获取模块81，具体用于通过下述步骤确定所述基准点：

从所述地震数据体中提取设定属性的第一地震属性体，从所述第一地震属性体中提取所述目的层切片处的第一地震属性切片；根据井目的层储层特征和对应的标注点，对所述第一地震属性切片中与该井的标注点对应的点进行储层特征标记；根据储层特征标记，从所述第一地震属性切片中确定基准点。

在一些实施例中，参考数据获取模块81，具体用于通过下述步骤确定所述参考点：

确定所述目的层切片与参考层间的厚度分布数据，所述参数层为所述上部参考层或下部参考层；根据所述厚度分布数据从所述目的层切片中确定参考点和所述参考点的漂移量。

基于本发明的发明构思，本发明实施例还提供一种储层预测装置，该装置的结构如图9所示，包括：

提取模块91，用于从地震数据体中提取设定属性的第二地震属性体，从所述第二地震属性体中提取等时地震切片处的第二地震属性切片，所述等时地震切片是按照上述方法从所述地震数据体中切取的；

标记模块92，用于根据井目的层储层特征对所述第二地震属性切片中与该井对应的点进行储层特征标记；

预测模块93，用于根据所述第二地震属性切片和储层特征标记，进行储层特征预测。

关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。

基于本发明的发明构思，本发明实施例还提供一种具备储层预测功能的计算机程序产品，包括计算机程序/指令，其中，该计算机程序/指令被处理器执行时实现上述等时地震切片获取方法，或实现上述储层预测方法。

应该明白，公开的过程中的步骤的特定顺序或层次是示例性方法的实例。基于设计偏好，应该理解，过程中的步骤的特定顺序或层次可以在不脱离本公开的保护范围的情况下得到重新安排。所附的方法权利要求以示例性的顺序给出了各种步骤的要素，并且不是要限于所述的特定顺序或层次。

在上述的详细描述中，各种特征一起组合在单个的实施方案中，以简化本公开。不应该将这种公开方法解释为反映了这样的意图，即，所要求保护的主题的实施方案需要清楚地在每个权利要求中所陈述的特征更多的特征。相反，如所附的权利要求书所反映的那样，本发明处于比所公开的单个实施方案的全部特征少的状态。因此，所附的权利要求书特此清楚地被并入详细描述中，其中每项权利要求独自作为本发明单独的优选实施方案。

上文的描述包括一个或多个实施例的举例。当然，为了描述上述实施例而描述部件或方法的所有可能的结合是不可能的，但是本领域普通技术人员应该认识到，各个实施例可以做进一步的组合和排列。因此，本文中描述的实施例旨在涵盖落入所附权利要求书的保护范围内的所有这样的改变、修改和变型。此外，就说明书或权利要求书中使用的术语“包含”，该词的涵盖方式类似于术语“包括”，就如同“包括，”在权利要求中用作衔接词所解释的那样。此外，使用在权利要求书的说明书中的任何一个术语“或者”是要表示“非排它性的或者”。