技术领域
[0001] 本
发明属于地球物理勘探领域，涉及一种叠前
地震记录AVA分析方法，特别涉及一种基于弹性波逆时偏移ADCIGs的叠前地震记录AVA分析方法。
背景技术
[0002] 地震资料振幅解释与反演经过多年发展，已经从精确的构造解释到“亮点”识别，再到利用AVO直接进行
碳氢检测。所谓“亮点”，是指地震剖面上由于地下油气藏的存在所引起的地震反射波振幅相对增强的“点”。地下
地层界面的反射系数对反射波振幅有直接的影响，其大小取决于界面两侧的波阻抗差，即取决于界面两侧的速度和
密度，而速度和密度与
岩石的孔隙度及孔隙中
流体性质有着密切的关系。
[0003] 亮点资料处理是地震资料处理中重要的分支之一，其主要目的是使得反射波的振幅尽可能接近于反射界面的反射系数之比，或至少使反射波振幅确定性地反映反射界面的反射系数，即得到相对保幅的剖面。此外，反射波振幅异常是指示油气藏存在的重要标志，但并非唯一标志。综合考虑各种异常，才能较为可靠地确定油气藏的空间
位置与分布，提高解释的成功率。
[0004] AVO/AVA技术是继亮点技术之后又一项利用振幅信息来研究岩性、检测油气的
地震勘探技术。经过多年发展和实践，该技术在寻找非背斜油气藏方面取得了较好的效果。在地震勘探中，
共中心点道集(CMP)记录内的偏移距经过必要的转换可以等价地用入射
角表示，因此AVO和AVA是等价的概念。所谓AVO技术，是利用
共中心点道集或共反射点(CRP)资料分析反射波振幅随偏移距(或入射角)的变化规律，估算界面两侧的弹性参数，进一步推断储层的岩性以及含油气性的一种资料处理和分析技术。该技术直接利用偏移距道集或角道集进行处理和分析，充分利用了多次
覆盖得到的原始地震信息，对岩性和含油气性的解释比亮点技术更为可靠，是解释
烃类异常的重要工具，已广泛应用于声学、弹性介质以及粘弹性介质中烃类饱和储层预测。
[0005] 然而，AVO/AVA分析需要高
质量、高保真度和大角度的角道集数据。常规的AVO分析大多是基于Kirchhoff偏移或单程波偏移生成偏移距域共成像点道集(ODCIGs)，然后利用射线追踪或射线参数逼近法将偏移距道集转换为角道集。当地下介质结构复杂时，基于Kirchhoff或单程波偏移生成的ODCIGs易受到波传播多路径的影响，偏移假象严重，特别是对大偏移距道集。因此，基于该类方法生成的角道集数据角度
分辨率低、AVA相对保幅性差，特别是大角度道集，为后续AVA分析与反演带来不确定性。
发明内容
[0006] 本发明的目的在于克服上述
现有技术的缺点，提供一种基于弹性波逆时偏移ADCIGs的叠前地震记录AVA分析方法。首先，利用弹性波方程对波场进行延拓，充分考虑波场的矢量特征以及丰富的转换波信息；然后，从逆时偏移中基于平面波分解法提取AVA幅度保持的角道集；然后，基于提取的保幅角道集进行AVA分析，实现振幅解释和烃类
异常检测。与传统的基于基于Kirchhoff或单程波偏移方法首先得到的ODCIGs进而转成角道集方法相比，本发明方法提取角道集没有任何运动学和动
力学假象，具有较高的
信噪比，减少了由
地震波在复杂介质多路径传播带来的偏移假象干扰和AVA不保幅问题，从而大大提高了AVA分析
精度与可靠性。
[0007] 本发明的目的是通过以下技术方案来解决的：
[0008] 该种基于弹性波逆时偏移ADCIGs的叠前地震记录AVA分析方法，包括以下步骤：
[0009] 1)采集原始共炮点地震记录，然后对采集到的地震资料做预处理，包括
切除直达波，去除折射波、鬼波或地滚波，称处理后得到共炮地震记录为偏移炮集记录，记为Qx(xR,xS；t)和Qz(xR,xS；t)，其中xR表示
检波器坐标，xS表示
震源坐标，t表示时间，Qx为位移
水平分量，Qz为位移垂直分量；
[0010] 2)基于速度分析构建纵波偏移速度Vp(x)和横波偏移速度Vs(x)，其中x＝(x,z)为空间坐标，不考虑密度，设密度ρ≡1；
[0011] 3)对偏移炮集记录Qx(xR,xS；t)和Qz(xR,xS；t)进行
频谱分析，构建宽频带的偏移子波函数w(t)；
[0012] 4)以偏移子波w(t)作为弹性波方程的震源项进行正向波场延拓，并对t时刻的位移矢量求取散度，得到正传P波波场pF(x；t)数据体；
[0013] 5)以偏移炮集记录Qx(xR,xS；t)和Qz(xR,xS；t)作为弹性波方程的边界条件进行逆时波场延拓，并对t时刻的位移矢量求取散度，得到反传P波波场pB(x；t)数据体；
[0014] 6)对震源波场pF(x；t)和接收波场pB(x；t)分别做三维傅里叶变换，得到
频率-
波数域震源波场 和接收波场
[0015] 7)计算震源波场 和接收波场 之间的夹角θ，然后在频率-波数域应用互相关成像条件，得到波数域角度域共成像点道集R(k；θ)；
[0016] 8)对波数域角道集R(k；θ)进行二维逆傅里叶变换，得到空间域角道集R(x；θ)；
[0017] 9)利用纵波偏移速度Vp(x)，将深度域角道集R(x；θ)转换到时间域R(x,t；θ)；
[0018] 10)利用Aki-Richards三项公式对R(x,t；θ)进行AVA分析，分别提取截距、梯度和
曲率属性，然后生成截距-梯度乘积属性、比例化泊松比变化属性、比例化密度变化属性；
[0019] 11)利用截距与梯度属性数据体进行AVA交汇分析，得到解释后地震剖面。
[0020] 进一步，上述步骤4)中，以偏移子波w(t)作为震源项进行正向波场延拓所用的弹性波方程如下所示
[0021]
[0022] 式中 ，λ和μ分别为
拉梅常数，与速度、密度之间关系为由于假设密度恒定为1，故 ux和uz分别为位移x
分量和z分量；σxx、σxz和σzz分别为
应力分量，对位移矢量u(x；t)＝(ux(x；t),uz(x；t))T求取散度可得纯P波波场
[0023]
[0024] 上述步骤5)中，以偏移炮集记录Qx(xR,xS；t)和Qz(xR,xS；t)作为边界条件进行逆时波场延拓所用的弹性波方程如下
[0025]
[0026] 对位移矢量u(x；t)＝(ux(x；t),uz(x；t))T求取散度获得纯P波波场，即[0027]
[0028] 上述步骤7)中，两个波场 和 之间的夹角θ可由下式计算得出
[0029]
[0030] 式中，kF和kB分别为 和 的波数矢量。然后在频率-波数域应用如下所示互相关成像条件
[0031]
[0032] 式中，δ为狄拉克函数，ω为频率，k为波数矢量，R(k；θ)为波数域角道集。
[0033] 上述步骤8)中，波数域角道集R(k；θ)的二维逆傅里叶变换公式为
[0034]
[0035] 式中，R(x；θ)为空间域角道集。
[0036] 上述步骤10)中，Aki-Richards三项公式表达式如下所示
[0037] R(θ)＝A+B sin2θ+C tan2θsin2θ
[0038] 式中，θ为反射角，A、B和C分别为截距、梯度和曲率。
[0039] 本发明具有以下有益效果：
[0040] 本发明提出一种基于弹性波逆时偏移角道集的叠前地震记录AVA分析方法。首先，利用弹性波方程对波场进行延拓，充分考虑波场的矢量特征以及丰富的转换波信息；然后，从逆时偏移中基于平面波分解法提取角道集；然后，基于提出的角道集进行AVA分析，实现振幅解释和烃类异常检测。与传统的基于基于Kirchhoff或单程波偏移方法得到的ODCIGs相比，基于弹性波逆时偏移生成的角道集没有任何运动学和动力学假象，受多路径干扰小，具有较高的信噪比，可真正反映地下介质的AVA属性，显著优于传统的ODCIGs。因此，ADCIGs可为AVO/AVA分析以及叠前弹性参数反演提供必要的高质量、高保真度叠前角度域数据体，有助于提高含油气性检测、储层刻画和流体识别的准确度和精度。
附图说明
[0041] 图1是本发明流程示意图；
[0042] 图2是弹性介质中波传播示意图；
[0043] 图3是波数-频率域计算波传播方向示意图。
[0044] 图4是Marmousi模型的P
波速度、S波速度和密度。
[0045] (a)纵波速度图；(b)横波速度图；(c)密度图；
[0046] 图5是偏移P波速度和偏移S波速度模型。
[0047] (a)纵波偏移速度图；(b)横波偏移速度图；
[0048] 图6是炮点位于3.0km处接收的位移x分量和z分量共炮地震记录。
[0049] (a)位移x分量图；(b)位移z分量图；
[0050] 图7是CDP 600处提取的角道集与AVA分析结果图。
[0051] 图8是基于Aki-Richards三项公式提取的截距与梯度乘积属性、比例化泊松比变化属性和比例化密度变化属性剖面。
[0052] (a)乘积属性图；(b)比例化泊松变化属性图；(c)比例化密度变化属性图；
[0053] 图9是截距与梯度属性体交汇分析解释后的地震剖面。
具体实施方式
[0054] 下面结合附图对本发明做进一步详细描述：
[0055] AVA分析的目的是利用反射地震振幅信息来研究地层岩性、检测并解释烃类异常的地震勘探技术。本发明采用一种新的叠前地震记录AVA分析流程与方法。首先，基于弹性波逆时偏移，利用波场分解
算法提取相对保持AVA振幅属性的深度域ADCIGs，然后利用偏移速度将其转换到时间域；然后，利用Aki-Richards三项公式提取截距、梯度和曲率属性，同时生成衍生属性，包括截距-梯度乘积属性、比例化泊松比变化属性、比例化密度变化属性；最后，对截距与梯度属性数据体进行交汇分析，解释烃类异常。
[0056] 本发明的物质
基础是通过野外高分辨率地震采集设备采集到的大量共炮地震数据。本发明的基于弹性波逆时偏移ADCIGs的叠前地震记录AVA分析
框架如图1所示，具体步骤分别为：
[0057] 1)采集原始共炮点地震记录，然后对采集到的地震资料做预处理，包括切除直达波，去除折射波、鬼波或地滚波，称处理后得到地震记录为偏移炮集记录，记为Qx(xR,xS；t)和Qz(xR,xS；t)，其中xR表示检波器坐标，xS表示震源坐标，t表示时间，Qx为位移水平分量，Qz为位移垂直分量，设最大记录时间为Tmax。
[0058] 2)基于速度分析构建纵波偏移速度Vp(x)和横波偏移速度Vs(x)，其中x＝(x,z)为空间坐标，设定密度恒为常数。
[0059] 3)对偏移炮集记录Qx(xR,xS；t)和Qz(xR,xS；t)进行频谱分析，构建宽频带的偏移子波函数w(t)。
[0060] 对共炮地震记录逐道沿时间方向做一维傅里叶变换，计算多道平均振幅谱，确定地震记录的有效频带范围为[ω1,ω2]，然后设计
通带为[ω1,ω2]的
滤波器，并通过逆傅里叶变换到时间域，得到的时间序列作为偏移子波w(t)。
[0061] 4)以偏移子波w(t)作为弹性波方程的震源项进行正向波场延拓，并对t时刻的位移矢量求取散度，得到正传P波波场pF(x；t)数据体。
[0062] 震源波场从0时刻到最大时刻Tmax进行正向延拓所用的弹性波方程如下所示[0063]
[0064] 式中，λ和μ分别为拉梅常数，假定密度恒为1，则拉梅常数与纵横波速度之间的关系为 ux和uz分别为位移x分量和z分量；σxx、σxz和σzz分别为应力分量。对位移矢量u(x；t)＝(ux(x；t)，uz(x；t))T求取散度可得纯P波波场
[0065]
[0066] 5)以偏移炮集记录Qx(xR,xS；t)和Qz(xR,xS；t)作为弹性波方程的边界条件进行逆时波场延拓，并对t时刻的位移矢量求取散度，得到反传P波波场pB(x；t)数据体。
[0067] 接收炮集Tmax时刻到0时刻进行逆时延拓所用的弹性波方程如下所示
[0068]
[0069] 对位移矢量u(x；t)＝(ux(x；t),uz(x；t))T求取散度可得纯P波波场，即[0070]
[0071] 6)对震源波场数据体pF(x；t)和接收波场数据体pB(x；t)分别做三维傅里叶变换，得到频率-波数域震源波场 和接收波场 首先对pF(x；t)和pB(x；t)沿时间方向逐道做一维傅里叶变换到空间-频率域
[0072]
[0073] 式中，ω为频率， 和 分别为频域震源波场和接收波场，然后对 和 分别对空间做二维傅里叶变换
[0074]
[0075] 式中，kF和kB分别为正传波场和反传波场的波数向量， 和 分别为频率-波数域震源波场和接收波场。
[0076] 7)计算震源波场 和接收波场 之间的夹角θ，然后在频率-波数域应用互相关成像条件，得到波数域角度域共成像点道集R(k；θ)，如图3所示。
[0077] 两个频率-波数域波场 和 之间的夹角θ可由波数矢量计算得出
[0078]
[0079] 然后在频率-波数域应用互相关成像条件
[0080]
[0081] 式中，δ为狄拉克函数，ω为频率，k为波数矢量，R(k；θ)为波数域角道集。
[0082] 8)对波数域角道集R(k；θ)进行二维逆傅里叶变换，得到空间域角道集，即[0083]
[0084] 式中，R(x；θ)为空间域角道集数据体。
[0085] 9)利用纵波偏移速度Vp(x)，
选定最大时间和时间
采样间隔，将深度域角道集R(x；θ)转换到时间域R(x,t；θ)。
[0086] 10)利用Aki-Richards三项公式对R(x,t；θ)进行AVA分析，分别提取截距和梯度乘积属性、比例化泊松比变化属性、比例化密度变化属性。
[0087] 考虑单反射界面弹性介质，如图2所示，上下层介质的纵波速度分别为VP1和VP2，横波速度分别为VS1和VS2，密度分别为ρ1和ρ2。对入射角为θ1的入射P波，遇到界面将会产生反射P波、反射SV波、透射P波、透射SV波，与法线的夹角分别为θ1、θ2、φ1和φ2。Wiggins推导出一种更直观的反射系数与上下界面弹性参数之间的关系式，即Aki-Richards三项公式[0088] R(θ)＝A+B sin2θ+C tan2θsin2θ
[0089] 式中，θ为反射角，A、B和C的具体表达式为
[0090]
[0091] A、B和C常被称之为截距属性、梯度属性和曲率属性。在AVA属性分析中，由上述三种属性数据体衍生的AVA属性被经常使用，包括乘积属性(A*B)，比例化的泊松比变化属性(A+B)，比例化横波反射系数(A-B)等。
[0092] 11)利用截距与梯度属性数据体进行AVA交汇分析，得到解释后地震剖面。截距与梯度属性交汇图基于Rutherford-Williams分类模型和泥岩线，在交汇图上对含油气
砂岩与含水砂岩进行分离，然后投影到叠后地震剖面，实现对含油气异常区的解释。
[0093] 数值仿真结果
[0094] 本节利用弹性波RTM中提取的PP波ADCIGs数据体，进行AVA分析，检测含油气异常区。图4给出了纵波速度、横波速度和密度模型。模型的横向宽度和深度分别为7.5km和2.5km，网格间距均为5.0m，网格大小为1501*501。震源从模型左侧50.0m处开始，炮间距
75.0m，共100炮覆盖。采用双边接收观测系统，横向每个网格均布置检波器，共1501个。震源和检波器的深度均为10.0m。震源选用主频为45Hz的Ricker子波。ADCIGs的角度范围为0～
60°，角度采样间隔2°。在P波速度模型中，存在两个较为显著的低速异常体；在密度模型中，存在多个低密度区域，如图中箭头所指示区域。图4给出了纵波偏移速度和横波偏移速度模型。图5给出了炮点位于3.0km处接收到的位移x分量和z分量共炮地震记录。
[0095] 首先利用平面波分解法从弹性波逆时偏移中提取深度域PP波ADCIGs，再利用偏移纵波速度将深度域ADCIGs数据体转换到时间域。图6为CDP 600处提取的角道集与AVA分析结果图，可见AVA显著，与
测井资料匹配程度高。然后利用Aki-Richards三项式分别提取梯度属性数据体A，截距属性数据体B和曲率属性数据体C。最后，利用上述属性数据体生成AVA衍生属性，包括截距和梯度乘积属性(A*B)、比例化泊松比变化属性(A+B)、和比例化密度变化属性(A-C)，如图8所示。在A*B属性中，可清楚看到异常体顶部和底部出现强反射；在A+B属性中，异常体顶部对应负变化，底部对应正变化；在A-C属性中，异常体顶部对应负变化，底部对应正变化。这些特征均符合第三类AVA异常的特点，并与实际模型参数一致。
[0096] 进一步，利用截距与梯度属性数据体进行AVA交汇分析，解释结果如图9所示。白色区域截距和梯度均为负，且振幅较大，解释为气砂顶部；黑色区域截距和梯度均为正，且振幅较大，解释为气砂底部。含油气砂顶部和底部分别被清晰地指示出来。
