CN102565857A

CN102565857A - 自动剩余动校正方法

Info

Publication number: CN102565857A
Application number: CN2011104241145A
Authority: CN
Inventors: 罗红明; 何光明; 陈爱萍; 曹中林; 张华�; 刘鸿; 金德刚; 陈三平
Original assignee: Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Current assignee: China National Petroleum Corp; BGP Inc
Priority date: 2011-12-16
Filing date: 2011-12-16
Publication date: 2012-07-11
Anticipated expiration: 2031-12-16
Also published as: CN102565857B

Abstract

本发明提供了一种自动剩余动校正方法。所述方法包括如下步骤：从接收的动校正后的多个道集中选择一个道集并对选择的道集所包括的多个地震道进行互相关时差校正；对选择的道集内的时差校正后的地震道进行叠加从而形成一个地震道数据；对叠加得到的地震道数据进行滤波处理从而形成选择的道集的模型道；确定所述模型道的极值点；从所述选择的道集所包括的时差校正前的多个地震道中选择一个地震道，确定所述选择的地震道中与所述模型道的极值点对应的极值点；基于选择的地震道的极值点与模型道的与选择的地震道的极值点对应的极值点的时差，对选择的地震道进行平移；对平移后的地震道进行插值，从而完成对所述选择的地震道的自动剩余动校正。

Description

自动剩余动校正方法

技术领域

本发明涉及地球物理勘探领域，具体地讲，涉及一种基于互相关的自动剩余动校正方法。

背景技术

复杂地区成像方法研究是石油地球物理勘探中长期存在并在目前尚未完全解决的前沿性研究课题，目前研究取得了很多理论和实际应用成果。但目前的研究大都集中于以下两个方面：

(1)研究如何使用各项异性动校正和高阶动校正方法

主要的研究包括两部分内容：一是提高叠加速度的计算精度，主要的方法利用各向异性速度分析得到叠加速度和各项异性参数等，从而实现动校正道集的中、远偏移距正确归位；二是提高动校正计算的精度，如高阶动校正方法等实现动校正道集的中、远偏移距正确归位。由于复杂构造区地震速度的影响因素很多，上述方法虽然在一定程度上提高了速度分析的精度，但并不能保证获得真正精确的叠加速度场。

(2)利用叠前偏移技术完成复杂构造的成像

目前该类方法是主要的研究热点，已取得了许多较好的研究成果，获得了较为广泛的应用。但这方面的研究主要依靠解释人员的经验进行地震速度拾取和调整，人工工作量很大，而且不同的处理解释人员根据他们的经验会得出差异较大的地震速度场，很难做到准确的地下速度场建立与成像。

虽然各向异性动校正和高阶动校正方法具有精度高的优点，但它们不适用于复杂构造成像，叠前偏移成像处理过程复杂，计算周期长，实际操作意义不大。

发明内容

为了解决现有技术中存在的上述问题，本发明提供了一种基于互相关的高效自动剩余动校正方法。

根据本发明的实施例的一种自动剩余动校正方法，包括如下步骤：a)接收动校正后的多个道集，其中，一个道集包括多个地震道；b)从接收的多个道集中选择一个道集并对选择的道集所包括的多个地震道进行互相关时差校正；c)对选择的道集内的时差校正后的地震道进行叠加从而形成一个地震道数据；d)对叠加得到的地震道数据进行滤波处理从而形成选择的道集的模型道；e)确定所述模型道的极值点；f)从所述选择的道集所包括的时差校正前的多个地震道中选择一个地震道，确定所述选择的地震道中与所述模型道的极值点对应的极值点；g)基于选择的地震道的极值点与模型道的与选择的地震道的极值点对应的极值点的时差，对选择的地震道进行平移；h)对平移后的地震道进行插值，从而完成对所述选择的地震道的自动剩余动校正；i)重复步骤f)和h)，直到对选择的道集的所有地震道均完成自动剩余动校正；j)重复步骤b)至i)，直到对所有的道集均完成了自动剩余动校正。

另外，在步骤d)中，所述滤波处理为带通滤波。

另外，在步骤b)中，在选择的道集所包括的多个地震道中选择一个地震道作为基准地震道，计算基准地震道与选择的道集中的另一地震道之间的互相关时差，并基于该互相关时差对该另一地震道进行时差校正。

另外，在步骤e)中，按预定窗口长度将模型道划分为N个子部分并确定每个子部分的极值点，其中，N为正整数。

另外，在步骤f)中，按预定窗口长度将选择的地震道划分为N个子部分并确定每个部分的极值点；基于模型道的第i子部分的极值点与选择的地震道的第i子部分的极值点的时差，平移选择的地震道的第i子部分，其中，i为不大于N的正整数。

另外，在步骤h)中，如果平移后的地震道中存在子部分重叠区域，则对重叠区域进行加权处理；如果平移后的地震道中存在子部分相隔而形成的空白区域，则对空白区域进行插值处理。

根据本发明的基于互相关的高效自动剩余动校正方法，适用于复杂山地的构造处理与成像，且计算步骤简便，计算效率高以及和成像效果好等特点，在目前我国的西部的复杂山地地震资料处理中具有广阔的应用前景

附图说明

图1是示出根据本发明示例性实施例的基于互相关的高效自动剩余动校正方法的流程图。

具体实施方式

现在对本发明实施例进行详细的描述，其示例表示在附图中，其中，相同的标号始终表示相同部件。下面通过参照附图对实施例进行描述以解释本发明。

下面，结合图1来详细描述根据本发明示例性实施例的基于互相关的高效自动剩余动校正方法。

在步骤100，接收动校正后的多个道集。如本领域所公知的，一个道集包括多个地震道。地震道表示的是一系列样点值(地震数据的样点值)，且每个样点值表示的是一个时间系列的某个时间点的振幅(也称为能量)。

在步骤110，进行互相关时差校正。如本领域所公知的，对地震道进行时差校正就是将样点值序列向前或向后移动，并去掉0毫秒前面和最大振幅值记录时间点后面的值，从而得到一个新的样点值序列。

具体地讲，从在步骤100中接收到的动校正后的多个道集中的选择一个道集(以下，为了便于描述，简称为选择的道集)，从选择的道集所包括的多个地震道中选择一个地震道作为基准地震道，计算选择的道集中的另一地震道与该基准地震道之间的互相关时差，并基于计算的互相关时差对所述另一地震道进行时差校正。这里，优选地，可将通过时差校正产生的空白部分插入0值，即，空白部分的振幅均为0。重复上述过程直到对选择的道集所包括的所有地震道均完成时差校正。

这里，计算互相关时差以及进行时差校正是本领域的公知常识，因此不对其进行详细描述。

在步骤120，对选择的道集内的时差校正后的地震道进行叠加从而形成一个地震道数据。

在步骤130，对叠加得到的地震道数据进行滤波处理从而形成选择的道集的模型道。这里，优选以带通滤波方式对叠加得到的地震道数据进行滤波。

在步骤140，获得所述模型道的极值点。具体地讲，按预定窗口长度将模型道划分为N个子部分(N为自然数)。这里，为了便于描述，假设模型道总长度为5秒，且预定窗口长度为500毫秒，则所述N为10个。在每个部分中选择一个极值点，即，选择一个振幅的绝对值最大的点。

在步骤150，从所述选择的道集所包括的时差校正前的地震道中选择一个地震道(以下，为了便于描述，简称为选择的地震道)，确定选择的地震道中与所述模型道的极值点对应的极值点。

对于选择的地震道，也按预定窗口长度对地震道进行划分，并在选择的地震道中的每个划分的部分中确定一个极值点(即，振幅的绝对值最大的点)。

例如，基于上述当假设，将选择的地震道按500毫秒的预定窗口长度划分为10个子部分，并从每个子部分中确定具有最大的振幅绝对值的极值点。

在步骤160，基于选择的地震道的极值点与模型道的与该极值点对应的极值点的时差，对选择的地震道进行平移。

例如，选择的地震道的第一子部分的极值点的时间为360毫秒，而在模型道中第一子部分的极值点的时间为380毫秒，则将选择的地震道的第一部分往后移20毫秒。相反，在模型到道中第一部分的极值点的时间为360毫秒，而选择的地震道的第一部分的极值点的时间为380毫秒，则将选择的地震道的第一部分往前移20毫秒。按这种方式对选择的地震道的所有子部分进行平移。

在步骤170，对平移后的地震道进行插值，从而完成对选择的地震道的自动剩余动校正。具体地讲，对于在平移后的地震道的重叠的区域进行加权处理；对于在平移后的地震道中的空白区域进行插值处理。

这里，为了便于描述，假设选择的地震道的第一子部分向后平移了10毫秒，第二子部分部分向前平移了5毫秒，第三子部分向后平移了5毫秒，则在平移后的地震道中，在0毫秒至10毫秒区域是空白(即，没有振幅)，495毫秒至510毫秒区域有振幅的重叠，在990毫秒至1005毫秒区域是空白。

对于空白区域，即0毫秒至10毫秒区域和990毫秒至1005毫秒区域进行差值；对于重叠区域，即495毫秒至510毫秒区域进行加权处理。

然后，重复步骤150至170直到对选择的道集内所有地震道均完成自动剩余动校正。在对选择的道集的所有地震道完成自动剩余动校正之后，重复步骤110至步骤170，对所有剩余道集完成自动剩余动校正。

根据本发明的基于互相关的高效自动剩余动校正方法，适用于复杂山地的构造处理与成像，且计算步骤简便，计算效率高以及和成像效果好等特点，在目前我国的西部的复杂山地地震资料处理中具有广阔的应用前景。

虽然已表示和描述了本发明的一些实施例，但本领域技术人员应该理解，在不脱离由权利要求及其等同物限定其范围的本发明的原理和精神的情况下，可以对这些实施例进行修改。

Claims

1.一种自动剩余动校正方法，其特征在于，包括如下步骤：

a)接收动校正后的多个道集，其中，一个道集包括多个地震道；

b)从接收的多个道集中选择一个道集并对选择的道集所包括的多个地震道进行互相关时差校正；

c)对选择的道集内的时差校正后的地震道进行叠加从而形成一个地震道数据；

d)对叠加得到的地震道数据进行滤波处理从而形成选择的道集的模型道；

e)确定所述模型道的极值点；

f)从所述选择的道集所包括的时差校正前的多个地震道中选择一个地震道，确定所述选择的地震道中与所述模型道的极值点对应的极值点；

g)基于选择的地震道的极值点与模型道的与选择的地震道的极值点对应的极值点的时差，对选择的地震道进行平移；

h)对平移后的地震道进行插值，从而完成对所述选择的地震道的自动剩余动校正；

i)重复步骤f)和h)，直到对选择的道集的所有地震道均完成自动剩余动校正；

j)重复步骤b)至i)，直到对所有的道集均完成了自动剩余动校正。

2.根据权利要求1所述的基于互相关的自动剩余动校正方法，其特征在于，在步骤d)中，所述滤波处理为带通滤波。

3.根据权利要求1所述的基于互相关的自动剩余动校正方法，其特征在于，在步骤b)中，在选择的道集所包括的多个地震道中选择一个地震道作为基准地震道，

计算基准地震道与选择的道集中的另一地震道之间的互相关时差，并基于该互相关时差对该另一地震道进行时差校正。

4.根据权利要求1所述的基于互相关的自动剩余动校正方法，其特征在于，在步骤e)中，

按预定窗口长度将模型道划分为N个子部分并确定每个子部分的极值点，其中，N为正整数。

5.根据权利要求4所述的基于互相关的自动剩余动校正方法，其特征在于，在步骤f)中，

按预定窗口长度将选择的地震道划分为N个子部分并确定每个部分的极值点；

基于模型道的第i子部分的极值点与选择的地震道的第i子部分的极值点的时差，平移选择的地震道的第i子部分，其中，i为不大于N的正整数。

6.根据权利要求5所述的基于互相关的自动剩余动校正方法，其特征在于，在步骤h)中，

如果平移后的地震道中存在子部分重叠区域，则对重叠区域进行加权处理；

如果平移后的地震道中存在子部分相隔而形成的空白区域，则对空白区域进行插值处理。