[发明专利]一种三维叠前地震数据层间多次波预测和压制方法

权利要求书

1.一种三维叠前地震数据层间多次波预测和压制方法，通过三维叠前地震数据层间多次波预测和三维叠前地震数据层间多次波匹配压制实现三维叠前地震数据层间多次波衰减；包括：首先通过对含有层间多次波的叠前数据进行预处理，获得规则化后的只含有一次反射波和层间多次波的叠前数据；然后在时间域根据产生层间多次波层位界面确定一次反射波同相轴，对同相轴的时间线进行记录；利用虚同相轴方法在频率-空间域构造虚一次反射波信号，再将虚一次反射波信号与原始的地震数据褶积构造层间多次波；进一步地，通过有效匹配滤波方法对预测的层间多次波进行匹配处理，达到压制三维叠前地震数据层间多次波，突出有效波的目的；其中：

A.进行三维叠前地震数据层间多次波预测包括如下步骤：

A1)规则化三维叠前地震数据；包括：

定义在表面接收的二维矩形观测系统，确定炮点和检波点的位置关系，根据野外采集的道头信息，提出每一道地震数据，将野外采集的地震数据进行规则化处理；

A2)准备层间多次波压制数据；包括：

将地震数据进行静校正；将规则化的地震数据校正到基准面上；同时拾取直达波信号进行切除，直达波不参与叠前地震数据层间多次波的预测；其中，所使用的地震数据为已经进行包括压制表面多次波的预处理的含有一次反射波和层间多次波的三维叠前地震数据；

A3)读取经过预处理的炮集数据，准备三维叠前地震数据层间多次波预测所需数据；

确定产生层间多次波层位的一次反射波信号，用时间框将有效信号包围，并记录包围有效信号同相轴的时间框信息；具体包括：

A31)在一个炮集数据中从上到下确定产生层间多次波层位的一次反射波，采用一个变化的时间框包围一次反射波同相轴，记录时间框数据信息，同时遍历该反射层对应一次反射波在所有炮集上的信息，并记录；将炮集上产生层间多次波的所有层位进行记录，作为层间多次波预测的输入数据；

A32)再规则化地震数据；

将已处理的三维叠前地震数据和一次反射波时间线数据按照传播-反射-传播WRW模型进行组装，得到规则化的WRW模型数据，便于层间多次波预测计算；所述WRW模型为二维数据矩阵，其中一个维度代表共炮点地震数据，另一个维度代表共检波点地震数据；

A4)预测三维叠前地震数据层间多次波；

根据时间线数据将炮集数据分成上和下两个部分的数据体，直接对三维叠前地震数据层间多次波进行准确预测；预测的具体过程包括：

A41)根据一次反射波的时间线数据将规则化的WRW模型数据分成上部分和下部分两个部分；

A42)进行计算获得虚同相轴数据；

A43)将虚同相轴数据与下部分地震数据执行褶积运算，预测出该层位对应的所有阶次的层间多次波；其中，对产生层间多次波地层的顶底界面倾角没有要求；

B.进行三维叠前地震数据层间多次波匹配压制，将三维叠前地震数据中的层间多次波进行多次迭代预测和匹配压制；包括：

B1)利用一范数和二范数匹配方法设计匹配滤波器；

匹配滤波器包括：

单道滤波器模块，只对地震数据中的一道进行范数匹配，获得单道滤波因子；

多道滤波器模块，以地震数据中某一道为中心道，其In-line方向的相邻地震道作为辅助道，获得多道匹配滤波因子；

方型滤波器模块，以地震数据中某一道为中心道，四周的地震道作为辅助数据，获得滤波因子；其中四周的地震道包括In-line和Cross-line方向的地震道；

B2)利用匹配滤波器对预测数据与实际数据进行匹配，生成匹配后叠前层间多次波数据；

进行匹配具体是：根据获得的滤波因子，与预测的叠前层间多次波相乘计算出匹配的层间多次波，

B3)最后获得层间多次波衰减后的叠前地震数据，得到压制层间多次波后的地震数据；并对层间多次波的预测和匹配压制过程进行多次迭代循环，最后获得层间多次波全部压制后的有效波数据；

所述迭代循环包括外循环和内循环；其中，移动一次反射波产生的迭代过程为外循环，用于压制不同层位产生的层间多次波；在外循环计算过程中的多次预测和匹配压制迭代过程为内循环，用于压制某一特定层位产生的所有阶次层间多次波；

通过上述步骤，实现维叠前地震数据层间多次波衰减。