AVO及叠前反演

第五章油气综合检测技术1第一节叠前流体检测技术近几年，随着地震采集处理技术的进步，尤其叠前偏移技术的发展和推广应用，使得研究人员可以得到来自地下真实反射点的叠前道集（CRP道集），为叠前烃类检测技术的发展奠定了资料基础。目前基于叠前道集的直接烃类检测方法主要有两种：一种是在岩石物理建模的基础上进行叠前道集AVO响应特征分析；一种是利用多个限角叠加数据体进行叠前弹性参数反演，利用纵横波波阻抗、纵横波速度比、泊松比、拉梅系数等敏感属性反映含油气性。一.AVO分析技术1、AVO理论简介AVO（AmplitudeVariationwithOffset），早先也称为AmplitudeVersusOffset，译为振幅随炮检距变化。由此而衍生的有振幅随入射角变化AVO（AmplitudeVariationwithAngle），振幅随方位角变化AVA（AmplitudeVariationwithAzimuth），振幅随炮检距和方位角变化AVOA（AmplitudeVariationwithOffsetandAzimuth）等。AVO作为一种含气砂岩的异常地球物理现象，最早在20世纪80年代初被Ostrander发现。这一现象表现为：当储层砂岩含气后，地震反射振幅随炮检距会发生明显的加大（基于SEG标准极性）。因为AVO现象与含气砂岩的对应关系，从而引起勘探地球物理界广泛的重视。后续的研究表明：这种异常现象并非一种特殊的形式，而是遵循Zoepprittz早先所提出的地震反射波动力学方程式，从而对AVO现象的解释有了完整的理论基础。针对AVO现象继而出现的AVO技术是继亮点之后又一项利用振幅信息研究岩性、检测油气的技术手段。AVO技术具有以下特点：A、直接利用CDP道集资料进行分析，这就充分利用的多次覆盖得到的丰富的原始信息；B、利用振幅随炮检距（入射角）的变化的特点，即利用整条曲线的特点。而亮点技术只是利用了这一特殊情况下曲线的一个数值。所以，AVO技术对岩性的分析比亮点技术更为可靠。C、这几年波动方程对地震剖面的成像有了更大的成果，是对地下构造形态的反演。AVO技术从严格意义上说算不上是利用波动方程进行岩性反演分析的方法，但是其理论和思路是对波动方程得到的结果的比较精确的利用。D、AVO技术是一种研究岩性的比较细致的方法，并且需要有测井资料的配合。2、AVO技术的理论基础振幅随炮检距的变化来自于所谓的“能量分区”。当地震波入射到地层界面时，一部分能量反射，一部分能量透射。如果入射角不等于零度，纵波（P波）能量一部分反射，一部分转化成透射P波和S波。反射和透射波的振幅能量取决于地层边界的物理性质差异。纵波速度Vp、横波速度Vs和密度ρ是非常重要的。同时，需要注意反射振幅也依赖于入射波的入射角（图5-3-1）。因此，当一个平面纵波非垂直入射到两种介质的分界面上，就要产生反射纵、横波和透射纵、横波。在界面上，根据应力连续性和位移连续性，依据边界条件并引入反射系数、透射系数，就可以得出四个相应波的位移振幅应当满足的方程叫做Zoeprritz方程，这个方程是Zoeprritz在1919年解出的。这个方程组比较复杂，不能解出新产生的波的振幅与有关参数明确的函数关系。但是从方程组可以看出，一般反射纵波的反射系数Rpp是入射角界面上部介质的密度ρ1，纵波速度Vp1，横波速度Vs1以及界面以下的介质密度ρ2，纵波速度Vp2，横波速度Vs2等七个参数的函数，可以简单的表示为Rpp（＆，Vp1，Vs1，ρ1，Vp2，Vs2，ρ2），虽然不能直接从方程中解出Rpp与七个参数的具体关系，但是可以假设以物质的六个物性参数为参变量，以＆为变量，仔细分析可以得到，六个参数是以两个参数的比值，例如21ppVV、21等形式出现。这样就可以把21ppVV、21等分别看作一第五章油气综合检测技术2个参数，再加上在同一种介质中，纵波速度Ｖp，横波速度Ｖs，以及泊松比σ之间又有关系，如：21)1(22VsVp，于是有关系式：Rpp=f（＆，21ppVV，21，σ），这样来达到减少参数的目的。从理论上说，在实际地震记录上得到某个界面的反射波的振幅与入射角的变化关系曲线，并且又知道某些参数，就可以利用曲线族作为量板来估算地层参数。图5-3-1：入射到地层边界的地震波示意图通过对地层弹性参数的研究，得到泊松比σ是一个对岩性和含油气情况反应比较敏感，参数，所以就要对方程进行适当的化简，得出以泊松比σ为参数的以σ为变量的简单的近似关系，即Rpp=f（＆，σ），这样通过反算σ，来达到对储层的参数测定和检测。由于Zoeppritz方程过于复杂，因此有许多学者尝试对其进行简化或近似，其中比较著名和实用的主要有：1）Koefoed的试算Koefoed在1955年第一个给出了用Zoeppritz方程，以R为参数计算出的Rpp～＆曲线。他用17组纵横波速度、密度和泊松比参数，较为详细地研究了泊松比对两个各向同性介质之间反射和折射面所产生的反射系数的影响，最大的入射角达到30°。他的研究结果被公认为Koefoed五原则。虽然Koefoed的结论说明了利用Rpp～＆曲线是可以反算出泊松比σ的，但是用未简化的Zoeppritz方程进行计算太复杂，因此反求弹性参数也是很复杂的。2）AkiKI和RichardsPG的纵波近似公式1980年AkiKI和RichardsPG给出了一个纵波反射振幅的一个近似表达式：R（θ）=)sin42secsin41(212222222VsVVVVVVVspsppps（5-3-1）这个公式是有适用条件的，就是界面两边的弹性介质性质的百分比变化小，式中的参数有：ΔVp=Vp2-Vp1、ΔVs=Vs2-Vs1、Vp=（Vp2+Vp1）/2、Vs=（Vs2+Vs1）/2、Δρ=ρ2-ρ1、ρ=（ρ2+ρ1）/2、θ=（θ1+θ2）/2，θ1是入射角，θ2是按斯奈尔定理计算的透射角。这个式子虽然对Zoeppritz方程进行了化简，但是用了Vs做参数而未用σ作参数。3)Bortfeld近似式1961年，Bortfeld提出近似式第五章油气综合检测技术3Rpp=211212122221211211122lnln)ln(2)()sin()coscosln(21spspppssppVVVVVVVVVpVV（5-3-2）4）Hilterman近似式211122211122coscoscoscos=RppppppVVVV21221221211123sinsssssspVVVVVVV（5-3-3）5）Shuey简化公式1985年Shuey利用上面的式子进行了简化，用σ代替Vs，修改得：R（θ）=2222000sintan21sin1VpVpRAR（5-3-4）其中：122122VpVs、VpVpVpVpB、121120BBA、0011RAA、12、212。垂直入射时的反射振幅VpVpR210，σ1、σ2分别为入射介质和透射介质的泊松比。其中，界面两侧泊松比的差Δσ是一个至关重要的因素，这就是振幅与炮检距关系研究的物理基础。Shuey近似式的特点就是其三项都有明确的物理意义：第一，垂直入射时，θ1=θ2=θ=0，R（0）=R0，即R0是垂直入射时的反射振幅。第二，在中等入射的情况(0〈θ〈30°），有近似tanθ≈sinθ，于是有：20200sin11RARR（5-3-5）此时，反射振幅与A有关，前两项起作用。这时的反射系数与介质的泊松比有密切关系，因此，利用此式更能突出油气特征。第三，对于大角度入射情况，反射振幅与速度变化有关中的第三项起主要作用。6）Simth和Gidlow加权叠加方法1987年Smith和Gidlow根据Gardner等人给出的水饱和岩石的密度与速度的4次方根成正比的假设，将Aki和Richards（1980）提出的Zoeppritz方程进行了修改化简：VsVCiVVBRsppipr（5-3-6）第五章油气综合检测技术4式中：iipsitgVVB222221sin2185、isVpVCi22sin24、Vp、Vs分别为纵横波速度，B、C为常数。因此，尽管表示反射系数的简化式有很多，但其最终目的都是为了求得一个最简表达式来表示反射系数随入射角的变化。而上、下介质的泊松比又对反射系数随入射角的变化起重要作用，因此，需要在表达式中包含标志油气特征的参数。Hilterman的简化方程与Shuey的简化方程满足了以上要求，且适用于各种层状模型的AVO模拟，因此得到了广泛的应用。在应用中最常见的一种简化公式为：2sin)(GPR（5-3-7）式中：R为反射系数、θ为入射角；P为近似为零偏移距下的纵波的反射振幅，也称AVO截距，其大小决定于上下层之间的纵波波阻抗差异；G为纵波反射振幅随入射角的变化梯度，也称AVO斜率，其大小取决于泊松比的变化。这个方程不能很好处理入射角较大的情况，但是简单易懂，至今仍然非常实用。对NMO（正常时差校正）道集中每一个时间采样点作sin2θ（θ为入射角）和振幅交汇图（图5-3-2），截距描述了正常入射角时的P波反射率（NIP），同时斜率就是梯度（振幅随入射角怎么变化）。式（5-3-7）表明，在入射角小于中等角度（一般为30°）时，纵波反射系数近似与入射角正弦值的平方成线性关系。基于该方程的叠前反演可以获得AVO属性参数P与G及其各种转换属性参数，对这些属性参数的解释是：a、由AVO截距组成的P剖面是一个真正的法线入射零炮检距剖面。b、由AVO斜率组成的G剖面反映的是岩层弹性参数的综合特征。c、在横波剖面上，当纵、横波速度比近似等于2时，P-G可以反映出横波波阻抗的特征。d、在泊松比剖面上，当纵、横波速度比近似等于2时，P+G反映的是泊松比的特征。e、截距与梯度的乘积（P×G）剖面，也称AVO强度剖面，更有利于识别气层。图5-3-2：sin2θ（θ为入射角）和振幅交汇图3、AVO的地质意义AVO应用的基础是泊松比的变化，而泊松比的变化是不同岩性和不同孔隙流体介质之间存在差异的客观事实。大量的试验研究和实践表明，沉积岩的泊松比值具有如下特点：a、未固结的浅层盐水饱和沉积岩往往具有非常高的泊松比值（0.4以上）；b、泊松比往往随孔隙度的减小及沉积固结程度而减少；c、高孔隙度的盐水饱和砂岩往往具有较高的泊松比值（0.3～0.4）；d、气饱和高孔隙度第五章油气综合检测技术5砂岩往往具有极低的泊松比值（如低到0.1）。一般来讲，不同岩性按泊松比值从高到低排序依次是：石灰岩、白云岩、泥岩和砂岩，在砂岩中由于孔隙流体性质的差异，依次是：水砂岩、油砂岩和气砂岩。从波速到泊松比，人们对地层的研究已逐渐深入到其本质。在研究振幅变化与波速的关系时，发现振幅变化与其他弹性参数的关系更有意义，这是因为波速不是一个独立的弹性参数，而是介质的几个弹性参数组合的结果。在拉梅常数（λ）、切变模量（μ）和泊松比（σ）等弹性参数中，泊松比在Zoeppritz方程中是以一个独立的参数出现的。而且，通过试验研究，发现泊松比是对区分岩性有特殊作用的一个参数。Gregory（1976）分别对各种岩性的泊松比做了大量测量试验，结果表明不同岩性泊松比的差别比速度的差别大。因此，利用泊松比判别岩性更可靠。所以，AVO技术的地质基础在于不同岩石以及含有不同流体的同类岩石之间泊松比存在差别。Domenico（1977）研究了含气、含油、含水砂岩的泊松比随埋藏深度的变化规律，结果发现含不同流体砂岩的泊松比随深度的变化特征是不同的：含气砂岩的泊松比随着深度的增加而增加，