[发明专利]地震数据处理方法和装置

说明书

技术领域

本发明涉及地质勘测领域，特别涉及一种地震数据处理方法和装置。

背景技术

地震勘探技术是利用人工激发的地震波在地层内传播来勘探地下的地质情况。在地面某处激发的地震波向地下传播时，遇到不同弹性的地层分界面就会产生反射波或折射波返回地面，用专门的仪器可记录这些波，分析所记录的波的特点(如波的传播时间、振动形状等)，能较准确地测定这些界面的深度和形态，判断地层的岩性。

由于地震记录是有噪音的，因此通常会通过不同位置进行激发、接收地下同一点的地震反射波，将这些来自同一点的地震反射波的集合称为地震道集。通过多次地震信号叠加，可以达到消除噪音，提高信噪比的目的。

由于野外采集的地震资料统一反射点的地震反射时间是不一样的，其反射时间与偏移距、地震传播速度有关，因此，需要将不同偏移距地震记录时间校正到自激自收地震反射时间位置，这个过程叫动校正，其目的是为了便于对地震信号叠加。

但是，因为地震波的传播速度是未知的，因此动校正的实现过程通常都设计一系列地震传播速度对地震道集进行动校正处理。

由于野外地震采集数据量非常庞大，人们不可能对每一个地震反射点都进行速度分析。往往都是间隔几十道才对一个地震反射点进行速度分析，而中间未进行速度分析的道集采用已作速度分析道集的速度进行插值而得，由于地震速度的非均质性，这些插值的速度或多或少都会存在一定的误差，动校正结果往往很难达到理想效果，很难将所有地震道集同相轴都校正到水平状态。

发明内容

为了解决现有技术的问题，本发明实施例提供了一种地震数据处理方法和装置。所述技术方案如下：

一方面，本发明实施例提供了一种地震数据处理方法，所述方法包括：

对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，S为整数；

计算各个地震道集与所述模型道的相关系数，并选出相关系数最大的第K地震道集；

在预设时间范围内滑动时窗，计算第K-1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K-1地震道集的剩余动校正量，K＞1且K为整数；

在所述预设时间范围内滑动时窗，计算第M地震道集与第M+1地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第M地震道集的剩余动校正量，K-2≥M≥1且M为整数；

在所述预设时间范围内滑动时窗，计算第K+1地震道集与所述第K地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第K+1地震道集的剩余动校正量；

在所述预设时间范围内滑动时窗，计算第N+1地震道集与第N地震道集的相关系数，并选取相关系数最大时的时窗作为最优点，将最优点对应的时移量作为第N+1地震道集的剩余动校正量，S-1≥N≥K+1且N为整数，所述最优点对应的时移量为时窗内的一个地震道集的地震反射时间与偏移距为0的地震道集的地震反射时间的差值；

采用计算出的所述S个地震道集的剩余动校正量对所述S个地震道集进行剩余动校正。

在本发明实施例的一种实现方式中，所述对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道，包括：

采用如下公式对多个地震道集进行叠加得到模型道：

其中，m1、m2为地震道集序号，a_ij为第i道的第j采样点的振幅值，m1、m2、i和j均为正整数，1≤m1＜m2≤S，i∈[m1，m2]。

在本发明实施例的另一种实现方式中，所述计算各个地震道集与所述模型道的相关系数，包括：

采用以下公式计算模型道与地震道集的相关系数：

其中，为第x道与模型道的相关系数，a_j为模型道第j个采样点的振幅值，Aa为模型道振幅平均值，b_j为地震道集第j个采样点的振幅值，Ba为地震道集振幅平均值，n为地震道集中采样点的总数，n、x和j均为正整数，且1≤x≤S。

在本发明实施例的另一种实现方式中，所述方法还包括：

在对动校正处理后的S个地震道集中预设范围内的地震道集进行叠加得到模型道之前，对所述各个地震道集进行振幅能量均衡化处理。

在本发明实施例的另一种实现方式中，所述对所述各个地震道集进行振幅能量均衡化处理，包括：

在设定时窗范围内统计振幅能量：

Q为统计道数，P为采样点数，a_j为振幅值；