AVO处理基本原理.doc

1．AVO处理的基本原理 AVO 即振幅与偏移距的关系，它是用来分析叠前地震资料的振幅特征以获得地下岩石和流体信息的一项重要技术。 AVO技术的理论基础是：地震波反射振幅与入射角有关。当地震波遇到一个界面时，部分能量返回地表，另一部分透射到界面以下。反射波的振幅与界面的反射系数成正比。当地震波非垂直入射时，部分能量转化为横波。转化的能量与反射的能量由下列参数确定： ΔVp— 纵波穿过地层界面时的速度变化 ΔVs— 横波穿过地层界面时的速度变化 Δρ— 地层界面两侧的密度差 θ — 入射角。入射角越大，转化为横波速度的能量就会越多。当垂直入射（θ=0）时，反射系数为R(θ)=f(ΔVp,ΔVs, Δρ) 此时不发生能量转化。 反射界面 转换纵波 反射纵波 入射波 图1 基本模型 Zoeppritz方程描述了地震反射率与入射角（炮检距）的关系，是一个很复杂的公式，一般不易求解。在该项目的AVO处理中使用了Zoeppritz方程的两个近似式，一个是Aki Richards公式，另一个是Shuey公式。 Aki Richards公式 （1） 根据泥岩基线公式（Mudrock line relation）， Vs=α1Vp+α0 （2） 或Greenberg Castagna 公式 Vs=α2 Vp2+α1Vp+α0 （3） 可以得到Vs/Vp。这样Aki Richards反演问题就变成了求解Δρ/ρ、ΔVp/Vp、ΔVs/Vs三个弹性变量。 原则上，直接求解Δρ/ρ、ΔVp/Vp、ΔVs/Vs是可能的，但考虑到解方程时矩阵的稳定性、有限的孔径、噪声和不完全的振幅恢复等原因，Δρ/ρ是不稳定的。因此，在实际的求解中将（1）式进一步简化。 用γ*ΔVp/Vp代替Δρ/ρ，（1）式可简化为， （4） 这样反演问题就变成了求解两个变量ΔVp/Vp和ΔVs/Vs。本项目反演的软件就是使用公式（4）求解ΔVp/Vp和ΔVs/Vs。 基于Aki Richards近似的AVO属性反射率意义 P-Wave Impedance ReflectivityP波阻抗变化率： S-Wave Impedance ReflectivityS波阻抗变化率：P-Wave Velocity ReflectivityP波速度变化率：ΔVp/VpS-Wave Velocity ReflectivityS波速度变化率：ΔVs/VsPseudo-Poisson Reflectivity伪泊松比变化率：Fluid Factor流体因子：Shuey公式 （5） 对于一个小的入射角，G表示相对于炮检距的反射率。尽管这个关系是对Zoeppritz方程很不准确的近似，但在描述AVO异常方面常常是有效的。 对于常规的叠加剖面，每一道的振幅值是所有炮检距道集振幅值的平均，而这种平均损失了包含在炮检距中的相对振幅信息。对于R0道，AVO数据是计算垂直入射的振幅，因此与叠加剖面相对，R0是比较好的零炮检距剖面。 在Shuey公式中，系数G是直接与弹性参数相关的参数。在Vs/Vp变化剧烈的地方都显示出绝对值高的G，或者说G绝对值高的点都是在Vs/Vp变化剧烈的地方。因为岩石孔隙中的气体会强烈影响Vs/Vp的比值（Ostrander，1984）。梯度剖面常常对这类含气储层有明显的指示。 基于Shuey近似的AVO属性反射率意义 Normal Incidence ReflectivityNI或R0，入射角为零时的反射系数，所以又称零炮检距剖面Gradient （G）梯度剖面Sign(NI)*GRNI与梯度（GR）的乘积剖面Poisson Reflectivity泊松比变化率：R0+G（Verm Hiltman, 1995）Lambda*Rho Reflectivity变化率: Mu*Rho Reflectivity变化率:Elastic Impedance Reflectivity弹性波阻抗变化率:应用Shuey公式计算,Near Normal Angle Stack，近角叠加；Wide Angle Stack，远角叠加。Medium Angle Stack，中等角度叠加。 远角度、中等角度、近角度的大小由用户自己来定义。 2．参数试验 2.1基于模型的预处理 在应用AVO技术之