[发明专利]一种基于弹性波逆时偏移ADCIGs的AVA分析方法

权利要求书

1.基于弹性波逆时偏移ADCIGs的AVA分析方法，其特征在于，包括以下步骤：

1)采集原始共炮点地震记录，然后对采集到的地震资料做预处理，包括切除直达波，去除折射波、鬼波或地滚波，称处理后得到共炮地震记录为偏移炮集记录，记为Q_x(x_R,x_S；t)和Q_z(x_R,x_S；t)，其中x_R表示检波器坐标，x_S表示震源坐标，t表示时间，Q_x为位移水平分量，Q_z为位移垂直分量；

2)基于速度分析构建纵波偏移速度V_p(x)和横波偏移速度V_s(x)，其中x＝(x,z)为空间坐标，不考虑密度，设密度ρ≡1；

3)对偏移炮集记录Q_x(x_R,x_S；t)和Q_z(x_R,x_S；t)进行频谱分析，构建宽频带的偏移子波函数w(t)；

4)以偏移子波w(t)作为弹性波方程的震源项进行正向波场延拓，并对t时刻的位移矢量求取散度，得到震源波场p_F(x；t)数据体；

5)以偏移炮集记录Q_x(x_R,x_S；t)和Q_z(x_R,x_S；t)作为弹性波方程的边界条件进行逆时波场延拓，并对t时刻的位移矢量求取散度，得到接收波场p_B(x；t)数据体；

6)对震源波场p_F(x；t)和接收波场p_B(x；t)分别做三维傅里叶变换，得到频率-波数域震源波场和接收波场

7)计算震源波场和接收波场之间的夹角θ，然后在频率-波数域应用互相关成像条件，得到波数域角度域共成像点道集R(k；θ)；

8)对波数域角道集R(k；θ)进行二维逆傅里叶变换，得到深度域角道集R(x；θ)；

9)利用纵波偏移速度V_p(x)，将深度域角道集R(x；θ)转换到时间域R(x,t；θ)；

10)利用Aki-Richards三项公式对R(x,t；θ)进行AVA分析，分别提取截距、梯度和曲率属性，然后生成截距-梯度乘积属性、比例化泊松比变化属性、比例化密度变化属性；

11)利用截距与梯度属性数据体进行AVA分析，得到解释后地震剖面。

2.根据权利要求1所述的基于弹性波逆时偏移ADCIGs的AVA分析方法，其特征在于，上述步骤4)中，以偏移子波w(t)作为震源项进行正向波场延拓所用的弹性波方程如下所示：

式中，λ和μ分别为拉梅常数，与速度、密度之间关系为由于假设密度恒定为1，故u_x和u_z分别为位移x分量和z分量；σ_xx、σ_xz和σ_zz分别为应力分量，对位移矢量u(x；t)＝(u_x(x；t),u_z(x；t))^T求取散度可得纯P波波场

3.根据权利要求1所述的基于弹性波逆时偏移ADCIGs的AVA分析方法，其特征在于，上述步骤5)中，以偏移炮集记录Q_x(x_R,x_S；t)和Q_z(x_R,x_S；t)作为边界条件进行逆时波场延拓所用的弹性波方程如下

式中，λ和μ分别为拉梅常数，与速度、密度之间关系为由于假设密度恒定为1，故u_x和u_z分别为位移x分量和z分量；σ_xx、σ_xz和σ_zz分别为应力分量，对位移矢量u(x；t)＝(u_x(x；t),u_z(x；t))^T求取散度获得纯P波波场，即