理论基础-bgp.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 Reflectivity P波阻抗变化率： S-Wave Impedance Reflectivity S波阻抗变化率： P-Wave Velocity Reflectivity P波速度变化率：ΔVp/Vp S-Wave Velocity Reflectivity S波速度变化率：ΔVs/Vs Pseudo-Poisson Reflectivity 伪泊松比变化率： Fluid Factor 流体因子： Shuey公式 （5）变化率: Mu*Rho Reflectivity 变化率: Elastic Impedance Reflectivity 弹性波阻抗变化率:应用Shuey公式计算, Near Normal Angle Stack，近角叠加；Wide Angle Stack，远角叠加。Medium Angle Stack，中等角度叠加。 远角度、中等角度、近角度的大小由用户自己来定义。 2．参数试验 2.1基于模型的预处理 在应用AVO技术之前，应对地震资料进行振幅保真性预处理。该项工作的目的是保证叠前地震资料的反射振幅与反射系数成正比。这意味着应该尽可能的消除地震波传播效应对反射振幅的影响。这些效应包括： 球面扩散 吸收衰减 排列影响 多次波 激发和接收效应等。 经试验后确定，在本项目AVO处理之前的预处理，主要是对CRP原始道集采用Q吸收补偿（Effective Q）、球面扩散补偿（Geometrical spreading）及带通中值滤波（Band pass median filter）处理； 经计算：本工区地震资料的入射角基本均在60。以下，因此在采用自动切除时定义的最大角度为60。；采用的角度门槛值（Angle Threshold），取临界角的百分比，即90％。 图2是预处理前后CRP道集的对比，预处理后CRP道集的质量明显高于预处理前。 图2 预处理前后CRP道集对比 2.2 泥岩基线 在AVO反演过程中需要用到泥岩基线的参数。泥岩基线用来定义地层水饱和砂岩的纵波速度波和横波速度波的线性关系。它是假定两者具有线性关系。 这就意味着：Vs=A*Vp + B，一般情况下可以从具有横波测井的资料中得到这种关系。 由于本工区没有横波测井资料，结合目的层的沉积岩性，按照经验公式确定纵横波之间