[发明专利]一种方位各向异性速度反演的方法及装置

说明书

技术领域

本发明涉及地球物理宽方位地震数据处理技术领域，尤其涉及一种方位各向异性速度反演的方法及装置。

背景技术

根据油气开发生产实践的经验总结，油气储层往往伴随着裂缝发育，所述裂缝可以对油气的运移产生重要的影响，因此裂缝检测是油气开发中的重要课题。宽方位地震勘探可以为裂缝检测提供有效的手段，裂缝构造在宽方位地震数据上可以表现为方位各向异性，即地震波特征随着观测方位的变化而变化。方位各向异性速度作为方位各向异性中的参数，不仅可以影响宽方位地震数据处理的效果，还可以为裂缝发育的解释提供有用的定量信息，因此，求取准确的方位各向异性速度是宽方位地震数据处理的关键环节。

现有技术中，方位各向异性速度可以通过分扇区的方式进行获取，上述方法中可以将地震数据按方位角分布情况划分成多个方位扇区，再对各个扇区进行速度分析，最后将所有扇区的速度拟合成方位各向异性速度。

但是，现有技术中对扇区中的数据进行人工速度分析，人工速度分析不仅消耗大量人力，还不可避免地带来拾取误差，影响最终的拟合结果。一般情况下，地震数据的方位分布是不规则的，分扇区后，某些扇区的覆盖次数较低，导致这些扇区的速度分析精度不够，影响最终的拟合结果，同时，分扇区一般将扇区内的方位信息进行平均，一个扇区对应着一个方位，因此，现有技术方法不能根据方位分布均衡覆盖次数，还降低方位精度，从而降低拟合后的方位各向异性速度精度。

发明内容

本发明的目的在于提供一种方位各向异性速度反演的方法及装置，可以充分利用宽方位的成像点道集数据，提高方位各向异性速度的反演精度。

本发明提供的一种方位各向异性速度反演的方法及装置是这样实现的：

一种方位各向异性速度反演的方法，所述方法包括：

获取炮检距矢量片叠前时间偏移后的成像点道集数据以及炮检距矢量片叠前时间偏移后的叠加数据；

在预设的处理时窗内对所述成像点道集数据和所述叠加数据作相关计算，生成成像点道集的剩余时差数据和相关系数数据；

根据所述剩余时差数据和相关系数数据计算所述成像点道集中零炮检距数据样点的方位各向异性速度。

可选的，在本发明的一个实施例中，所述在预设的处理时窗内对所述成像点道集数据和所述叠加数据作相关计算，包括：

获取当前成像点道集的线号和点号，从所述叠加数据中提取与所述线号和点号相对应的叠加数据道，将所述叠加数据道作为模型数据道；

根据预设的处理时窗计算得到互相关时窗；

在所述互相关时窗内，从所述模型数据道中截取与所述互相关时窗相对应的模型数据道序列，并从所述成像点道集数据的第一数据道中提取与所述互相关时窗相对应的第一数据道序列；

在所述互相关时窗内对所述模型数据道序列和所述第一数据道序列作相关计算，生成所述互相关时窗的相关系数序列。

可选的，在本发明的一个实施例中，所述生成成像点道集的剩余时差数据和相关系数数据，包括：

选取所述相关系数序列中的最大相关系数，计算与所述最大相关系数对应的序号；

根据所述序号、所述成像点道集数据的采样间隔计算所述互相关时窗中点的剩余时差值以及相关系数；

对所述互相关时窗中点的剩余时差值以及相关系数进行线性插值，得到所述第一数据道中每个采样点的剩余时差值以及相关系数。

可选的，在本发明的一个实施例中，所述计算与所述最大的相关系数对应的序号，包括：

根据所述最大相关系数以及所述最大相关系数的两个相邻相关系数，采用三点逆抛物线内插法计算得到所述最大相关系数对应的序号。

可选的，在本发明的一个实施例中，所述根据所述剩余时差数据和相关系数数据计算所述成像点道集中零炮检距数据样点的方位各向异性速度，包括：

根据所述成像点道集的剩余时差数据、相关系数数据以及叠前时间偏移数据，采用加权最小平方算法计算成像点道集在零炮检距数据样点的方位各向异性速度。

一种方位各向异性速度反演的装置，所述装置包括：

数据获取单元，用于获取炮检距矢量片叠前时间偏移后的成像点道集数据以及炮检距矢量片叠前时间偏移后的叠加数据；

剩余时差计算单元，用于在预设的处理时窗内对所述成像点道集数据和所述叠加数据作相关计算，生成成像点道集的剩余时差数据和相关系数数据；

各向异性速度计算单元，用于根据所述剩余时差数据和相关系数数据计算所述成像点道集中零炮检距数据样点的方位各向异性速度。