CN103217710A

CN103217710A - 获取共偏移距共方位角域道集的方法及系统

Info

Publication number: CN103217710A
Application number: CN2013100815578A
Authority: CN
Inventors: 段文胜; 李飞
Original assignee: China Petroleum and Natural Gas Co Ltd
Current assignee: China Petroleum and Natural Gas Co Ltd
Priority date: 2013-03-14
Filing date: 2013-03-14
Publication date: 2013-07-24
Anticipated expiration: 2033-03-14
Also published as: CN103217710B

Abstract

本发明涉及一种在交叉排列观测系统中获取共偏移距共方位角域道集的方法及系统，包括：将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在交叉排列观测系统中炮线的数目，m为在交叉排列观测系统中检波线的数目；根据炮线距和检波线距等距离，对所述交叉排列地震数据道集进行划分获取共偏移距共方位角单元，并对共偏移距共方位角单元进行编号；将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集。

Description

获取共偏移距共方位角域道集的方法及系统

技术领域

本发明涉及地震数据处理领域，特别涉及一种在交叉排列观测系统中获取共偏移距共方位角域道集的方法及系统。

背景技术

为了提高碳酸盐岩缝洞体成像和裂缝预测精度，目前常用的技术手段是通过采集全方位高密度三维地震资料，虽然通过三维地震资料的处理和裂缝预测可以取得非常好的效果。但是，从所采用的整体技术方案来看仍存在一些不足。其一，整体是采用针对窄方位的处理技术来应对全方位资料，对全方位资料处理应有的技术特点没有充分的研究和挖掘，也未形成全方位资料处理技术流程；其二，整体上看，现行技术方案是一种串联解决思路。首先，用VTI各向异性叠前时间偏移得到最终成果数据体，解决层状各向异性问题；然后，用分方位处理进行裂缝预测来解决方位各向异性问题。这个方案能较好地进行裂缝预测，但有一个关键问题往往被忽略，就是最终成果数据体上只解决层状各向异性问题而没有考虑到方位各向异性效应，不同方位角处理得到的数据体如何重构成最终数据体又是一个技术难点。产生上述问题根本原因在于没有采用宽方位/全方位处理技术。

实际碳酸盐岩地层既有成层性，会产生层状介质具有的VTI各向异性，又发育裂缝，会产生裂缝介质具有的HTI方位各向异性，即总体上表现为正交各向异性特征。传统的共偏移距偏移技术对VTI各向异性问题解决得较好，但因共偏移距偏移不能保存方位角信息，使得后续方位各向异性分析不太可行。目前的做法是通过分方位处理，将数据划分成不同扇区各自重新偏移后再做一些分析拟合的工作。然而这种方法由于扇区内覆盖次数有限往往成像很差，比如扩大扇区则各向异性参数分析又不准，比如扩大面元形成超道集则精度又受影响。而且，不同方位角处理得到的数据体如何重构成最终数据体还没有解决方案。

发明内容

本发明的目的是针对现有常规地震处理技术不适应宽方位地震采集数据的局限性，提供在交叉排列观测系统中获取共偏移距共方位角域道集的方法及系统，基于该技术方案获取的共偏移距共方位角域道集在地震数据处理可以有效保留方位角信息。

为实现上述目的，本发明提供了一种在交叉排列观测系统中获取共偏移距共方位角域道集的方法，包括：

将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在交叉排列观测系统中炮线的数目，m为在交叉排列观测系统中检波线的数目；

根据炮线距和检波线距等距离，对所述交叉排列地震数据道集进行划分获取共偏移距共方位角单元，并对共偏移距共方位角单元进行编号；

将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集。

可选的，在本发明一实施例中，还包括：

对所述共偏移距共方位角域道集内的数据先后进行去噪声、插值以及规则化处理；

获取每个地震数据道集的偏移距和方位角，根据偏移距和方位角来获取处理后的所述共偏移距共方位角域道集的偏移距和方位角；

根据所述共偏移距共方位角域道集的偏移距和方位角，对所述共偏移距共方位角域道集进行偏移处理获取共偏移距共方位角域共成像点道集；

利用剩余方位动校正方法获取消除方位各向异性的共偏移距共方位角共成像点道集。

可选的，在本发明一实施例中，所述共偏移距共方位角单元的个数等于地震数据的覆盖次数。

可选的，在本发明一实施例中，所述对共偏移距共方位角单元进行编号的步骤包括：

按照象限和顺序对共偏移距共方位角单元进行编号。

可选的，在本发明一实施例中，所述偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。

为实现上述目的，本发明还提供了一种在交叉排列观测系统中获取共偏移距共方位角域道集的系统，包括：

地震数据处理单元，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在交叉排列观测系统中炮线的数目，m为在交叉排列观测系统中检波线的数目；

划分编号单元，用于根据炮线距和检波线距等距离，对所述交叉排列地震数据道集进行划分获取共偏移距共方位角单元，并对共偏移距共方位角单元进行编号；

共偏移距共方位角域道集获取单元，用于将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集。

可选的，在本发明一实施例中，还包括：

共偏移距共方位角域数据处理单元，用于对所述共偏移距共方位角域道集内的数据先后进行去噪声、插值以及规则化处理；

偏移距方位角单元，用于获取每个地震数据道集的偏移距和方位角，根据偏移距和方位角来获取处理后的所述共偏移距共方位角域道集的偏移距和方位角；

偏移处理单元，用于根据所述共偏移距共方位角域道集的偏移距和方位角，对所述共偏移距共方位角域道集进行偏移处理获取共偏移距共方位角域共成像点道集；

消除方位各向异性单元，用于利用剩余方位动校正方法获取消除方位各向异性的共偏移距共方位角共成像点道集。

可选的，在本发明一实施例中，所述划分编号单元进一步按照象限和顺序对共偏移距共方位角单元进行编号。

可选的，在本发明一实施例中，所述偏移处理单元进行的偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。

上述技术方案具有如下有益效果：基于本发明提出的技术方案获取的共偏移距共方位角域进行地震数据处理，可以将产生采集脚印的数据空间不连续性稀疏分布的特质，就是将不连续性稀疏地分布在整个工区，幅度小且跳跃性小。因为空间不连续性稀疏分布了，采集脚印也变得稀疏了，也避免像常规处理那样过多使用规则化和内插技术。

另外，基于本发明提出的技术方案获取的共偏移距共方位角域偏移在限定偏移距和方位角的前提下进行偏移，更能提高偏移结果的精度，提高横向分辨率。传统的共偏移距域偏移后，方位角信息丢失。共偏移距共方位角域偏移后的共成像点道集具有方位角信息，有利于开展叠前反演，也使得方位AVO（AZAVO）分析成为可能。共偏移距共方位角偏移后无需再做分方位处理，偏移后的共成像点道集可用于求取方位各向异性和预测裂缝，因此其具有识别油藏的各向异性和裂缝特征的潜力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提出的获取共偏移距共方位角域道集的方法流程图；

图2为基于图1获取的共偏移距共方位角域道集进行的地震处理的方法流程图；

图3为本发明提出的获取共偏移距共方位角域道集的系统结构图；

图4为基于图3获取的共偏移距共方位角域道集进行的地震处理的系统结构图；

图5为实施例中某一十字排列地震数据道集偏移距分布图；

图6为实施例中某一十字排列中共偏移距共方位角数据单元划分以及编号分布图；

图7为实施例中第一象限中编号为1组成的共偏移距共方位角域地震数据道集图；

图8为实施例中第一象限中编号为19组成的共偏移距共方位角域地震数据道集图；

图9为实施例中第一象限中编号为37组成的共偏移距共方位角域地震数据道集图；

图10为实施例中第一象限中编号为55组成的共偏移距共方位角域地震数据道集图；

图11为实施例中某一未去噪的共偏移距共方位角域的地震数据道集图；

图12为实施例中某一去噪后的共偏移距共方位角域的地震数据道集图；

图13为实施例中共偏移距共方位角域地震数据道集叠前时间偏移叠加剖面图；

图14为实施例中剩余方位动校正前的共偏移距共方位角域共成像点道集图；

图15为实施例中剩余方位动校正后的共偏移距共方位角域共成像点道集图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，为本发明提出的获取共偏移距共方位角域道集的方法流程图。包括：

步骤101）：将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在交叉排列观测系统中炮线的数目，m为在交叉排列观测系统中检波线的数目；

步骤102）：根据炮线距和检波线距等距离，对所述交叉排列地震数据道集进行划分获取共偏移距共方位角单元，并对共偏移距共方位角单元进行编号；

步骤103）：将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集。

如图2所示，为基于图1获取的共偏移距共方位角域道集进行的地震处理的方法流程图。包括：

步骤103）：将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集；

步骤104）：对所述共偏移距共方位角域道集内的数据先后进行去噪声、插值以及规则化处理；

步骤105）：获取每个地震数据道集的偏移距和方位角，根据偏移距和方位角来获取处理后的所述共偏移距共方位角域道集的偏移距和方位角；

步骤106）：根据所述共偏移距共方位角域道集的偏移距和方位角，对所述共偏移距共方位角域道集进行偏移处理获取共偏移距共方位角域共成像点道集；

步骤107）：利用剩余方位动校正方法获取消除方位各向异性的共偏移距共方位角共成像点道集。

在图1和图2中，所述共偏移距共方位角单元的个数等于地震数据的覆盖次数。

所述对共偏移距共方位角单元进行编号的步骤包括：按照象限和顺序对共偏移距共方位角单元进行编号。

所述偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。

如图3所示，为本发明提出的获取共偏移距共方位角域道集的系统结构图。包括：

地震数据处理单元301，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在交叉排列观测系统中炮线的数目，m为在交叉排列观测系统中检波线的数目；

划分编号单元302，用于根据炮线距和检波线距等距离，对所述交叉排列地震数据道集进行划分获取共偏移距共方位角单元，并对共偏移距共方位角单元进行编号；

共偏移距共方位角域道集获取单元303，用于将每个交叉排列地震数据道集中同一象限相同编号的共偏移距共方位角单元提取出来组合到一起，即获取共偏移距共方位角域道集。

如图4所示，为基于图3获取的共偏移距共方位角域道集进行的地震处理的系统结构图。包括：

共偏移距共方位角域数据处理单元304，用于对所述共偏移距共方位角域道集内的数据先后进行去噪声、插值以及规则化处理；

偏移距方位角单元305，用于获取每个地震数据道集的偏移距和方位角，根据偏移距和方位角来获取处理后的所述共偏移距共方位角域道集的偏移距和方位角；

偏移处理单元306，用于根据所述共偏移距共方位角域道集的偏移距和方位角，对所述共偏移距共方位角域道集进行偏移处理获取共偏移距共方位角域共成像点道集；

消除方位各向异性单元307，用于利用剩余方位动校正方法获取消除方位各向异性的共偏移距共方位角共成像点道集。

在图3和图4中，所述共偏移距共方位角单元的个数等于地震数据的覆盖次数。

所述划分编号单元302进一步按照象限和顺序对共偏移距共方位角单元进行编号。

所述偏移处理单元306进行的偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。

实施例：

首先进行地震数据准备。包括原始地震数据读入、地震数据编辑、观测系统加载和预处理等。

然后，将三维地震数据分选成十字排列地震数据道集；本实施例中，某实际地震资料观测系统共有炮线46条，检波线53条。所以，十字排列地震数据数目为46×53个，即2438个。如图5所示，为实施例中某一十字排列地震数据道集偏移距分布图。

对分选好的每个十字排列地震数据道集按照炮线距和检波线距等距离划分成若干共偏移距共方位角小矩形数据，并按照象限和顺序对每个矩形进行编号；本实施例是利用某实际地震资料（炮线距与检波线距都是330米）进行，划分的小矩形数据个数等于本三维地震数据的覆盖次数225次。如图6所示，为实施例中某一十字排列中共偏移距共方位角数据单元划分以及编号分布图。图6显示的小矩形个数为16×16个，仔细观察可以看出，边界处的小矩形数据都是半个，所以完整的小矩形数据仍然是15×15个，等于覆盖次数225次。按照象限和一定顺序将小矩形数据进行编号。

将所有十字排列地震数据道集中相同象限相同编号的小矩形数据提取出来组成一个相应的共偏移距共方位角域的地震数据道集；如图7所示，为实施例中第一象限中编号为1组成的共偏移距共方位角域地震数据道集图。如图8所示，为实施例中第一象限中编号为19组成的共偏移距共方位角域地震数据道集图。如图9所示，为实施例中第一象限中编号为37组成的共偏移距共方位角域地震数据道集图。如图10所示，为实施例中第一象限中编号为55组成的共偏移距共方位角域地震数据道集图。本实施例中小矩形数据个数为16×16个，即256个；按照相同象限相同编号抽取组合，我们可以得到共偏移距共方位角域的地震数据道集的数目等于小矩形的编号数目，也是256个；每个相应的共偏移距共方位角域的地震数据道集都代表了一定的偏移距和方位角范围。

在得到的共偏移距共方位角域的地震数据道集上可以进行相应数据域的去噪以及插值、规则化的地震数据处理；如图11所示，为实施例中某一未去噪的共偏移距共方位角域的地震数据道集图。如图12所示，为实施例中某一去噪后的共偏移距共方位角域的地震数据道集图。图11与图12对比可知，共偏移距共方位角域的地震数据道集的去噪处理可以明显改善偏移前道集的信噪比。

计算每个地震数据道集的偏移距和方位角，根据获取的每个地震数据道集的偏移距和方位角来获取共偏移距共方位角道集的偏移距和方位角，之后可以进行共偏移距共方位角的叠前时间偏移或叠前深度偏移；也可以将偏移距相等，方位角方向相反（即不同象限，相同编号）的两个共偏移距共方位角域地震数据作为输入，形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。如图13所示，为实施例中共偏移距共方位角域地震数据道集叠前时间偏移叠加剖面图。共偏移距共方位角域地震数据道集叠前时间偏移叠加剖面显示信噪比高，且碳酸盐岩串珠反射更清晰。

共偏移距共方位角域地震偏移形成共偏移距共方位角域共成像点道集，利用剩余方位动校正方法，形成消除方位各向异性的共偏移距共方位角地震共成像点道集；如图14所示，为实施例中剩余方位动校正前的共偏移距共方位角域共成像点道集图。如图15所示，为实施例中剩余方位动校正后的共偏移距共方位角域共成像点道集图。共偏移距共方位角域共成像点道集剩余方位动校正前后可知，消除了方位各向异性。从方位各向异性提取的这些有用信息可以直接展示为各种反映储层裂缝方向分布和裂缝密度分布的信息。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims

1.在交叉排列观测系统中获取共偏移距共方位角域道集的方法，其特征在于，包括：

2.如权利要求1所述的方法，其特征在于，还包括：

3.如权利要求1或2所述的方法，其特征在于，所述共偏移距共方位角单元的个数等于地震数据的覆盖次数。

4.如权利要求1或2所述的方法，其特征在于，所述对共偏移距共方位角单元进行编号的步骤包括：

按照象限和顺序对共偏移距共方位角单元进行编号。

5.如权利要求2所述的方法，其特征在于，所述偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。

6.在交叉排列观测系统中获取共偏移距共方位角域道集的系统，其特征在于，包括：

7.如权利要求6所述的系统，其特征在于，还包括：

8.如权利要求6或7所述的系统，其特征在于，所述共偏移距共方位角单元的个数等于地震数据的覆盖次数。

9.如权利要求6或7所述的系统，其特征在于，所述划分编号单元进一步按照象限和顺序对共偏移距共方位角单元进行编号。

10.如权利要求7所述的系统，其特征在于，所述偏移处理单元进行的偏移处理包括共偏移距共方位角的叠前时间偏移，共偏移距共方位角的叠前深度偏移，偏移距相等、方位角方向相反的两个共偏移距共方位角域地震数据作为输入来形成互换共偏移距共方位角地震数据进行叠前时间偏移或叠前深度偏移。