[发明专利]一种基于薄层反射理论的AVA反演方法

摘要

本发明公开了一种基于薄层反射理论的AVA反演方法，首先，实际地震AVA数据为薄层实地震道，利用Hilbert变换，获得薄层虚地震道，由此构建薄层复地震道。利用FFT变换，将时间域的薄层复地震道及地震子波变换到频率域，两者做除法运算由此获得频率域中实际的薄层PP波反射系数。给定初始薄层模型参数，并利用薄层PP波反射系数近似式计算该薄层的PP波反射系数。模型更新利用高斯‑牛顿算法获得，当实际PP波反射系数与模型PP波反射系数之间的差异达到精度要求，输出当前模型参数即为薄层AVA反演结果。相比常规AVA反演方法，本发明给出的基于薄层反射理论的AVA反演方法可以更好地指示薄层的存在，同时获得更高精度的薄层弹性参数及厚度的反演结果。

主权项

1.一种基于薄层反射理论的AVA反演方法，其特征在于，包括如下步骤：S1、根据实际地震记录，获得时间域的薄层实AVA道集及地震子波；S2、利用Hilbert变换，由薄层实AVA道集求得薄层虚AVA道集，并结合薄层实AVA道集和薄层虚AVA道构造时间域的薄层复AVA道集；S3、针对时间域的薄层复AVA道集及地震子波，利用FFT变换，分别获得频率域的薄层复AVA道集及地震子波，两者通过除法运算，即可获得频率域中实际的薄层PP波复反射系数矩阵R^obs；S4、建立初始的薄层模型m；所述薄层模型m包含薄层厚度以及薄层三层介质的纵波速度、横波速度、密度；S5、利用薄层PP波反射系数近似式获得对应薄层模型m的PP波反射系数矩阵R；S6、求解目标函数Q(m)＝‖R^obs‑R‖²,并判断Q(m)是否满足精度要求；S7、当Q(m)小于阈值ε时，Q(m)满足精度要求，输出反演结果m；当Q(m)不小于ε时，Q(m)不满足精度要求，利用高斯‑牛顿算法，计算更新量Δm，令m＝m+Δm作为新的薄层模型返回步骤S5。