[发明专利]一种横波偏移速度建模方法

说明书

技术领域

本发明涉及石油地球物理勘探技术领域，尤其涉及一种应用于转换波叠前时间偏移处理的横波偏移速度建模方法。

背景技术

在现有的转换波偏移速度分析和建模技术中，最具有代表性的方法是首先基于精确纵波偏移速度模型和初始横波偏移速度模型对整条测线逐个道集进行转换波叠前时间偏移处理，在偏移处理的成像道集上逐点进行反射波组同相轴偏移效果作横波速度解释拾取，然后根据解释拾取的新的横波速度模型再次进行转换波叠前时间偏移处理。这其中需要再次检查偏移道集中反射波组同相轴是否过偏(速度小)或欠偏(速度大)，从而判断是否需要再次修改速度。由于对速度的评价取决于人眼视觉对成像道集校平的判断以及相关经验，在人机交互解释过程中人为因素较多，因此工作量繁琐而巨大。图1显示了目前大多数偏移速度分析所采用的一种交互界面。其中，图1(a)显示了偏移处理的成像道集(CIP)，图1(b)显示了人机交互解释拾取的横波速度谱。在具体实施时，通常是综合速度谱上的相似性能量团大小和成像道集的偏移效果进行人机交互解释，只有当完成整条测线上所有成像控制点的分析，才能最终得到横波偏移速模型。

另外一种经典的方法是由英国地质调查局各向异性研究项目(EAP)所建立的等效转换波速度建模的方法。该方法不需要纵波偏移速度模型，而是将下行纵波和上行横波的转换波等效为一种C波，应用多参数各向异性速度分析技术，获得转换波等效波速度、有效速度比和各向异性等三个参数。同前述第一种方法一样，这种方法也存在人机交互解释过程，需要根据人眼视觉对成像道集校平的判断和相关经验评价速度和其它参数，人为因素较多。图2显示了C波多参数各向异性偏移速度分析的一种交互界面。同样地，在具体实施时只有完成整条测线上所有成像控制点的分析，才能最终得到C波偏移速模型。

目前所有的转换波偏移速度分析和建模技术，都必须在初始纵、横波或者C波各向异性速度模型的基础上，进行分析点上的叠前时间偏移处理，输出成像道集(CIP)和计算相应的相似性能量的速度谱。由于工作量巨大，因此往往必须借助于计算机的GPU或CPU并行加速技术，才能将偏移速度分析实用化。此外，这其中还需要检查偏移修正的成像道集中反射波组同相轴是否过偏(速度小)或欠偏(速度大)，人机交互解释工作量巨大，是资料处理中最繁重的工作环节。有鉴于此，本发明的技术人员一直致力于寻找一种快速有效的方法建立横波偏移速度模型，以提高转换波地震资料叠前时间偏移处理的效率。

发明内容

针对上述问题，本发明的目的是提供一种快速的横波偏移速度建模方法，为时间一致性转换波地震资料叠前时间偏移处理提供有效的、合理的横波偏移速度模型。

一种横波偏移速度建模方法，包括以下步骤：

S100，分别对纵波地震资料和转换波地震资料进行叠加处理，获取纵波和转换波的等效速度终叠剖面，并建立纵波偏移速度模型以及获取纵波叠前时间偏移剖面；

S200，对纵波终叠剖面和转换波终叠剖面进行层位匹配处理，求取两个终叠剖面之间的时移序列，根据两个终叠剖面之间的时移序列反演纵横波速度比；

S300，根据纵波偏移速度模型和纵横波速度比，建立横波偏移速度模型；

S400，根据纵波偏移速度模型和横波偏移速度模型，对转换波地震资料进行叠前时间偏移处理，获取其时间与纵波旅行时一致的转换波叠前时间偏移剖面；

S500，对纵波叠前时间偏移剖面和转换波叠前时间偏移剖面进行层位匹配处理，求取两个偏移剖面之间的时移序列；

S600，判断两个偏移剖面之间的时移序列是否大于给定的精度阈值：

如果是，进入步骤S700；

如果否，进入步骤S900；

S700，根据两个偏移剖面之间的时移序列反演纵横波速度比微量变化，根据纵横波速度比微量变化修正纵横波速度比；

S800，根据纵波偏移速度模型和修正后的纵横波速度比，更新横波偏移速度模型，返回步骤S400；

S900，输出横波偏移速度模型。

根据本发明的实施例，所述步骤S200中，优选采用动态图像变形最优化技术对两个终叠剖面进行层位匹配处理，求取两个终叠剖面之间的时移序列。

根据本发明的实施例，所述步骤S500中，优选采用动态图像变形最优化技术对两个偏移剖面进行层位匹配处理，求取两个偏移剖面之间的时移序列。

根据本发明的实施例，所述步骤S200中，在反演纵横波速度比之后，还以纵波速度场的时间-速度对为控制点，对纵横波速度比做层内的趋势平滑处理。