一种方位各向异性速度反演的方法及装置的制造方法
【技术领域】
[0001] 本发明涉及地球物理宽方位地震数据处理技术领域，尤其涉及一种方位各向异性 速度反演的方法及装置。
【背景技术】
[0002] 根据油气开发生产实践的经验总结，油气储层往往伴随着裂缝发育，所述裂缝可 以对油气的运移产生重要的影响，因此裂缝检测是油气开发中的重要课题。宽方位地震勘 探可以为裂缝检测提供有效的手段，裂缝构造在宽方位地震数据上可以表现为方位各向异 性，即地震波特征随着观测方位的变化而变化。方位各向异性速度作为方位各向异性中的 参数，不仅可以影响宽方位地震数据处理的效果，还可以为裂缝发育的解释提供有用的定 量信息，因此，求取准确的方位各向异性速度是宽方位地震数据处理的关键环节。
[0003] 现有技术中，方位各向异性速度可以通过分扇区的方式进行获取，上述方法中可 以将地震数据按方位角分布情况划分成多个方位扇区，再对各个扇区进行速度分析，最后 将所有扇区的速度拟合成方位各向异性速度。
[0004] 但是，现有技术中对扇区中的数据进行人工速度分析，人工速度分析不仅消耗大 量人力，还不可避免地带来拾取误差，影响最终的拟合结果。一般情况下，地震数据的方 位分布是不规则的，分扇区后，某些扇区的覆盖次数较低，导致这些扇区的速度分析精度不 够，影响最终的拟合结果，同时，分扇区一般将扇区内的方位信息进行平均，一个扇区对应 着一个方位，因此，现有技术方法不能根据方位分布均衡覆盖次数，还降低方位精度，从而 降低拟合后的方位各向异性速度精度。
【发明内容】
[0005] 本发明的目的在于提供一种方位各向异性速度反演的方法及装置，可以充分利用 宽方位的成像点道集数据，提高方位各向异性速度的反演精度。
[0006] 本发明提供的一种方位各向异性速度反演的方法及装置是这样实现的：
[0007] -种方位各向异性速度反演的方法，所述方法包括：
[0008] 获取炮检距矢量片叠前时间偏移后的成像点道集数据以及炮检距矢量片叠前时 间偏移后的叠加数据；
[0009] 在预设的处理时窗内对所述成像点道集数据和所述叠加数据作相关计算，生成成 像点道集的剩余时差数据和相关系数数据；
[0010] 根据所述剩余时差数据和相关系数数据计算所述成像点道集中零炮检距数据样 点的方位各向异性速度。
[0011] 可选的，在本发明的一个实施例中，所述在预设的处理时窗内对所述成像点道集 数据和所述叠加数据作相关计算，包括：
[0012] 获取当前成像点道集的线号和点号，从所述叠加数据中提取与所述线号和点号相 对应的叠加数据道，将所述叠加数据道作为模型数据道；
[0013] 根据预设的处理时窗计算得到互相关时窗；
[0014] 在所述互相关时窗内，从所述模型数据道中截取与所述互相关时窗相对应的模型 数据道序列，并从所述成像点道集数据的第一数据道中提取与所述互相关时窗相对应的第 一数据道序列；
[0015] 在所述互相关时窗内对所述模型数据道序列和所述第一数据道序列作相关计算， 生成所述互相关时窗的相关系数序列。
[0016] 可选的，在本发明的一个实施例中，所述生成成像点道集的剩余时差数据和相关 系数数据，包括：
[0017] 选取所述相关系数序列中的最大相关系数，计算与所述最大相关系数对应的序 号；
[0018] 根据所述序号、所述成像点道集数据的采样间隔计算所述互相关时窗中点的剩余 时差值以及相关系数；
[0019] 对所述互相关时窗中点的剩余时差值以及相关系数进行线性插值，得到所述第一 数据道中每个采样点的剩余时差值以及相关系数。
[0020] 可选的，在本发明的一个实施例中，所述计算与所述最大的相关系数对应的序号， 包括：
[0021] 根据所述最大相关系数以及所述最大相关系数的两个相邻相关系数，采用三点逆 抛物线内插法计算得到所述最大相关系数对应的序号。
[0022] 可选的，在本发明的一个实施例中，所述根据所述剩余时差数据和相关系数数据 计算所述成像点道集中零炮检距数据样点的方位各向异性速度，包括：
[0023] 根据所述成像点道集的剩余时差数据、相关系数数据以及叠前时间偏移数据，采 用加权最小平方算法计算成像点道集在零炮检距数据样点的方位各向异性速度。
[0024] -种方位各向异性速度反演的装置，所述装置包括：
[0025] 数据获取单元，用于获取炮检距矢量片叠前时间偏移后的成像点道集数据以及炮 检距矢量片叠前时间偏移后的叠加数据；
[0026] 剩余时差计算单元，用于在预设的处理时窗内对所述成像点道集数据和所述叠加 数据作相关计算，生成成像点道集的剩余时差数据和相关系数数据；
[0027] 各向异性速度计算单元，用于根据所述剩余时差数据和相关系数数据计算所述成 像点道集中零炮检距数据样点的方位各向异性速度。
[0028] 可选的，在本发明的一个实施例中，所述剩余时差计算单元还包括：
[0029] 模型数据道确定单元，用于获取当前成像点道集的线号和点号，从所述叠加数据 中提取与所述线号和点号相对应的叠加数据道，将所述叠加数据道作为模型数据道；
[0030] 互相关时窗计算单元，用于根据预设的处理时窗计算得到互相关时窗；
[0031] 数据道序列提取单元，用于在所述互相关时窗内，从所述模型数据道中截取与所 述互相关时窗相对应的模型数据道序列，并从所述成像点道集数据的第一数据道中提取与 所述互相关时窗相对应的第一数据道序列；
[0032] 相关计算单元，用于在所述互相关时窗内对所述模型数据道序列和所述第一数据 道序列作相关计算，生成所述互相关时窗的相关系数序列。
[0033] 可选的，在本发明的一个实施例中，所述剩余时差计算单元还包括：
[0034] 最大相关系数选取单元，用于选取所述相关系数序列中的最大相关系数，计算与 所述最大相关系数对应的序号；
[0035] 中点剩余时差计算单元，用于根据所述序号、所述成像点道集数据的采样间隔计 算所述互相关时窗中点的剩余时差值以及相关系数；
[0036] 插值单元，用于对所述互相关时窗中点的剩余时差值以及相关系数进行线性插 值，得到所述第一数据道中每个采样点的剩余时差值以及相关系数。
[0037] 可选的，在本发明的一个实施例中，所述最大相关系数选取单元还包括：
[0038] 三点逆抛物线内插计算单元，用于根据所述最大相关系数以及所述最大相关系数 的两个相邻相关系数，采用三点逆抛物线内插法计算得到所述最大相关系数对应的序号。
[0039] 可选的，在本发明的一个实施例中，所述各向异性速度计算单元还包括：
[0040] 加权最小平方法计算单元，用于根据所述成像点道集的剩余时差数据、相关系数 数据以及叠前时间偏移数据，采用加权最小平方算法计算成像点道集在零炮检距数据样点 的方位各向异性速度。
[0041] 由此可见，本发明一种方位各向异性速度反演方法及装置的技术方案可以通过充 分利用宽方位成像点道集数据，在进行剩余时差校正后计算得到方位各向异性速度，不仅 提高了方位各向异性速度的反演精度，使用校正后的道集还可以提高同相轴振幅的拾取精 度，从而进一步得到较高的成像精度。
【附图说明】
[0042] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可以根据这些附图获得其他的附图。
[0043] 图1是本发明提供的方位各向异性速度反演方法的一种实施例的方法流程图；
[0044] 图2是本发明提供的方位各向异性速度反演装置的一种实施例的模块结构示意 图；
[0045] 图3是本发明提供的剩余时差计算单元的一种实施例的模块结构示意图；
[0046] 图4是本发明提供的剩余时差计算单元的另一种实施例的模块结构示意图。
【具体实施方式】
[0047] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实 施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术