[发明专利]一种转换波角道集抽取方法

说明书

技术领域

本发明涉及石油地球物理勘探领域，具体是一种用于转换波AVA分析、转换波振幅分析和地震属性反演中的转换波角道集抽取方法。

背景技术

在多波多分量地震勘探中，有效利用转换波信息能够获得更加可靠的岩石弹性参数。从转换波AVA分析中能够提取有用的岩性信息，用来预测地层岩性和裂缝等介质信息。在转换波AVA方程中，自变量是转换波入射角。转换波AVA分析中的一个基本步骤就是进行将炮检距域道集转换为角度域道集，即进行转换波入射角的计算。目前，有两种方法可以计算转换波入射角：一是直射线法；二是射线追踪法。

前者在已知转换波的转换点位置和反射层深度的条件下，用直射线法求取转换波入射角；后者在已知层速度和层厚度的情况下，用射线追踪法求取转换波入射角。

直射线法的计算速度快，但是其计算的角度总是偏小，误差较大；射线追踪法的计算结果准确，但是其运算量大，计算速度慢。

发明内容

本发明目的在于提供一种提高计算转换波入射角的速度和精度的转换波角道集抽取方法。

本发明采用以下步骤实现：

1)在野外利用纵波震源激发地震波并利用三分量检波器记录地震波，按照多波多分量地震资料处理流程对采集的纵波入射和横波反射（PS）转换波数据进行相对振幅保持的高保真处理，形成振幅随偏移距变化（AVO）的正常时差校正（NMO）后的道集；

2)在横波反射（PS）转换波数据处理过程中通过速度分析得到横波反射（PS）转换波时差校正（NMO）速度V_C2、垂直速度比γ₀、有效速度比γ_eff和PS转换波各向异性参数χ_eff；

3)由横波反射（PS）转换波时差校正（NMO）速度V_C2、垂直速度比γ₀、有效速度比γ_eff和PS转换波各向异性参数χ_eff，计算纵波速度V_P2和纵波各向异性参数η；

4)由横波反射（PS）转换波时间t_C、横波反射（PS）转换波炮检距x、横波反射（PS）转换波时差校正（NMO）速度V_C2、垂直速度比γ₀和有效速度比γ_eff，计算与横波反射（PS）转换波炮检距x对应的纵波炮检距x'；

所述的横波反射（PS）转换波炮检距x与纵波炮检距x'具有相同的震源O和相同的转换点水平位置x_C；

5)由横波反射（PS）转换波时间t_C、横波反射（PS）转换波炮检距x、横波反射（PS）转换波时差校正（NMO）速度V_C2、垂直速度比γ₀、有效速度比γ_eff和PS转换波各向异性参数χ_eff，计算零炮检距的转换波旅行时t_C0；

6)由零炮检距的转换波旅行时t_C0和垂直速度比γ₀计算零炮检距的纵波旅行时t_P0；

7)由零炮检距的纵波旅行时t_P0、纵波速度V_P2、纵波各向异性参数η和纵波炮检距x'得到纵波时差方程t_P；

8)由零炮检距的纵波旅行时t_P0和纵波速度V_P2得到层间纵波速度V_int；

9)由纵波时差方程t_P和层间纵波速度V_int，用射线参数法计算PS转换波入射角θ；

10)计算横波反射（PS）转换波数据中每一样点的入射角，将此样点放置在入射角的角度道上，对所有样点重复此过程，得到横波反射（PS）转换波角道集。

本发明容易实现，实验结果证明提高了抽取转换波角道集的速度和精度。

附图说明

图1是转换波入射角计算示意图，图中纵波反射：P波入射，P波反射；PS转换波反射：P波入射，S波反射；

图2是转换波炮检距道集，图中纵轴为时间，单位秒，横轴为道号，道间距为100米，偏移距为0米；

图3是不同方法计算的入射角结果对比，根据表1中的转换波速度数据，分别用本发明、射线追踪方法和直射线法计算的转换波炮检距为1000米的不同深度界面的入射角的结果；