一种油气检测的频变反演方法

摘要

本发明提供了一种油气检测的频变反演方法，属于油汽地球物理领域。本发明方法包括：(1)输入叠前地震记录；(2)根据步骤(1)输入的叠前地震记录生成叠前角道集；(3)对步骤(2)得到的叠前角道集进行频谱分解获得分频角道集记录；(4)利用步骤(3)的分频角道集记录以及分频角道集记录对应的频率进行频变AVO反演，获得纵波频散梯度；(5)利用步骤(4)得到的纵波频散梯度预测油气储层。本发明的反演结果要比利用已有Wilson‑Wu的频变AVO反演公式进行反演得到的结果更准确；另外，本发明方法中待反演量前的系数是线性无关的，所以利用本发明方法不会出现反演结果不唯一的情况。

主权项

一种油气检测的频变反演方法，其特征在于：所述方法包括以下步骤：(1)输入叠前地震记录；(2)根据步骤(1)输入的叠前地震记录生成叠前角道集；(3)对步骤(2)得到的叠前角道集进行频谱分解获得分频角道集记录；(4)利用步骤(3)的分频角道集记录以及分频角道集记录对应的频率进行频变AVO反演，获得纵波频散梯度；(5)利用步骤(4)得到的纵波频散梯度预测油气储层，所述步骤(3)是这样实现的：利用信号平滑伪Wigner‑Ville分布的频谱分解方法对叠前角道集进行频谱分解获得分频角道集记录，所述步骤(4)中是利用下面的公式进行频变AVO反演：$\begin{array}{c}R(\theta ,f)\approx A_2\left(\theta \right)\frac{{\mathrm{\Delta V}}_p}{V_p}\left(f_0\right)+(f-f_0)A_2\left(\theta \right)I_{a2}\\ +B_2\left(\theta \right)\frac{{V_s}^2}{{V_p}^2}\left(\frac{{\mathrm{\Delta V}}_s}{V_s}+\frac{1}{8}\frac{{\mathrm{\Delta V}}_p}{V_p}\right)+(f-f_0)B_2\left(\theta \right)I_{b2}\end{array}---\left(1\right)$其中I_a2是纵波频散梯度；I_b2是混合剩余频散梯度；R是反射系数，ΔV_p是界面上下层纵波速度差，ΔV_p＝V_p2‑V_p1，V_p是界面上下层的平均纵波速度，f₀是参考处的频率值，f是频率，V_s是界面上下层的平均横波速度，ΔV_s是界面上下层横波速度差，ΔV_s＝V_s2‑V_s1，θ是入射角；下标2表示界面下层对应物理量，下标1表示界面上层对应物理量；系数A₂、B₂的表达形式如下：$A_2\left(\theta \right)=\frac{5}{8}+\frac{1}{2}{\tan }^2\left(\theta \right)---\left(2\right)$ B₂(θ)＝‑4sin²(θ)    (3)。