AVO及叠前反演讲解
- 格式:docx
- 大小:171.52 KB
- 文档页数:8
叠前反演技术介绍宋长青 帕拉代姆北京办事处内容提要一、 反演的技术思路 二、地震资料识别岩性技术的理论基础 三、AVO技术介绍2反演的技术思路3反演工作中面临的问题基础数据的质量分析 反演算法的选择 反演结果的质量监控 反演属性定量分析 当某一项目涉及多项技术的综合,如岩石物理技术、 AVO技术、振幅反演技术、地震模拟技术和其它辅助技 术,如地震相分类,往往使研究者望而生畏。
4帕拉代姆解决方案使油藏描述的数据转换、数据分析和数据解释等多项 工作在交互的环境下同步实现。
Interpretation岩石到流体解释技术Rock & Fluid Interpretationd ui Fl & g ck gin Ro ImaRock & Fluid CharacterizationPetrophysicsGeophysicsConnecting Disciplines5反演及分析解释技术• • • • • • 基于井资料的岩石物性分析技术 – 地层评价 属性建模技术 – 背景模型 AVO 反演, 正演和解释技术 – AVO 属性 AI/EI 反演, 模型建立与解释技术 – AI/EI, 波松比, 密度, 拉枚常数, 剪切摸量 基于神经网络油藏物性反演技术 – 空隙度, 泥质含量, 饱和度 反演属性的标定及解释 – 油藏属性空间解释6结论• 关键技术Petrophysics GeologVolum interpretation 3D CanvasAI/EI VanguardFacies Classification StratimagicVisualization& Interpretation VoxelGeoAVO Probe7地震资料识别岩性技术的理论基础k + 3P-wave Velocity:Vp =Vs =ρμ ρμ 4ρ: Density, Mass per unit Volume.k: Bulk Modulus, IncompressibilityS-wave Velocity:μ: Shear Modulus, Rigidity⎛Vp ⎜ ⎜V ⎝ s⎞ K 4 ⎟ = + ⎟ μ 3 ⎠2岩层中地震波的速度决定于弹性模量和密度,岩石的弹性模量 又首先决定于岩石的矿物成分,其次是孔隙度、孔隙流体性质以及 压力、温度等环境因素,而孔隙度、孔隙流体及环境因素是通过影 响岩石的弹性模量和密度而影响速度的,所以决定岩石速度的最重 要因素是岩石成分,因此我们自然想到用速度来判别岩性。
叠前反演技术,与叠前弹性反演技术、叠前地震反演技术和定量AVO都是指同一概念。
该技术是利用叠前CRP道集数据(或部分叠加数据)、速度数据(一般为偏移速度)和井数据(横波速度、纵波速度、密度及其他弹性参数资料),通过使用不同的近似式反演求解得到与岩性、含油气性相关的多种弹性参数并进一步,用来预测储层岩性、储层物性及含油气性。
为什么要进行地震资料的叠前反演呢?首先,由于地震资料野外采集是多炮多道的观测系统,每个CDP点或CMP点记录的不同道具有不同的炮检距,每道上的反射振幅随炮检距不同而变化。
叠后反演基于常规处理的水平叠加数据,以自激自收为假设条件,即每个CDP或CMP道集经动校正后,把不同炮检距的记录道动校正为零炮检距位置,之后进行水平叠加。
这样,叠加剖面无法反应野外采集所记录的振幅随炮检距变化的特性,并损失了与炮检距关系密切的大量横波信息。
其次是叠后波阻抗反演是不随入射角发生变化,仅与纵波速度、密度有关,而叠前反演的弹性阻抗与入射角密切相关并与纵波、横波速度、密度4项参数有关。
由于同时利用了纵横波速度,其计算产生的弹性参数远较叠后反演丰富,可区别岩性与含油气性,为钻探提供更丰富、更准确的依据。
技术人员在研究中发现:进行叠前反演时,要注意资料条件及处理解释的结合。
一是地震资料的采集必须针对目的层深度,有足够大的炮检距来记录大量信息,并在处理中,对振幅进行补偿,严格保持相对振幅关系,避免虚假振幅信息的产生。
二是在地震资料道集进行部分叠加时,炮检距或角度范围的选择要针对目的层深度,使不同炮检距范围能明显反应用振幅的变化。
三是至少需要3个以上的子波,子波振幅谱对应于不同炮检距部分叠加数据。
四是在纵横波资料分析中,当岩石中含有油气时,纵波速度会降低,有时会出现含油气砂岩的速度接近泥岩速度,在声阻抗上无法区分岩性,但横波阻抗受油气影响很少,因此,两者的交汇图分析对划分岩性及含油气意义深远。
五是弹性参数综合分析,其物理意义不同,有的反应弹性模量,有的反应剪断模量,必须综合分析,才能做出合理解释。
一引言岩层中地震波的速度决定于弹性模量和密度,岩石的弹性模量又首先决定于岩石的矿物成分,其次是孔隙度、孔隙流体性质以及压力、温度等环境因素,而孔隙度、孔隙流体及环境因素是通过影响岩石的弹性模量和密度而影响速度的,所以决定岩石速度的最重要因素是岩石成分,因此我们自然想到用速度来判别岩性。
然而,各种岩石的速度范围太宽,互相重叠,我们很难仅仅根据速度来作岩性判别。
因此利用AVO和叠前弹性反演进行岩性预测越来越引起人们的重视。
二叠前反演方法原理(一)AVO的方法原理AVO分析技术是利用地震反射振幅与炮检距变化的关系 (Amplitude—Versus Offset,简称AVO),即:通过分析CDP道集中不同炮检距的地震反射,来识别岩性及检测含气性的一种地震技术。
其物理意义是:在两种不同岩层之间的界面上,当一种岩层的纵、横波速度之比Vp/Vs与另一种岩层的速度之比明显不同时,其反射系数随入射角(炮检距)而变化。
AVO反演主要利用不同岩性泊松比差异所形成的AVO特征响应,来区分波阻抗相近的储层与非储层。
当地震纵波P1非垂直入射到两种介质分界面上时,会产生反射波和透射波,其中反射波包括反射纵波P2和反射横波S1,透射波包括透射纵波P3和透射横波S2(图1-1图1-1纵波倾斜入射到界面产生的反射波和透射波示意图纵波非垂直入射,反射系数和透射系数满足Zeoppritz方程:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∙⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡------11112112222122111122111222221112221112221122112cos 2sin cos sin 2sin 2cos 2sin 2cos 2cos 2sin 2cos 2sin sin cos sin cos cos sin cos sin φθθθφρρφρρφφφρρθρρφθφθφθφθφθPS PP PS PP P S P P P S S P S P S P S S P T T R R V V V V V V V V V V V V V V V (1-1)其中介质1表示入射波和反射波所在的介质,介质2表示透射波所在的介质,VP1、VS1分别表示介质1中纵、横波的速度;VP2和VS2分别表示介质2中纵、横波的速度;ρ1和ρ2分别表示介质1和介质2的密度;θ1表示纵波的入射角和反射角;φ1表示横波的反射角,θ2和φ2分别表示透射纵波、横波的透射角;RPP 、RPS 、TPP 、TPS 分别表示P 波反射系数、SV 波反射系数、P 波透射系数和SV 波透射系数。
叠前AVO地震反演方法概述(刘文劼 095211068)AVO是一项利用振幅随偏移距变化特征分析和识别岩性及油气藏的地震勘探技术。
理论分析表明:振幅系数随入射角变化与分界面两侧岩石的弹性参数有关,它是通过非常复杂的非线性关系与介质的密度p1和p2、纵波速度a1和a2、横波速度b1和b2及入射角联系起来。
振幅系数随入射角变化本身隐含了岩性参数的信息,利用AVO关系可以直接反演岩石的密度p、纵波速度a和横波速度b,定量进行地震油藏描述。
波阻抗反演是零炮检距数据模型反演,那么AVO分析就是非零炮检距数据反演。
由于通过叠加得到零炮检距剖面,一方面丢失与炮检距有关的信息,另一方面叠加道又不是真正的零炮检距道,致使反演结果的稳定性以及它的应用均受到了一定程度的制约,AVO分析是在具有不同炮检距道集上进行分析,充分应用了叠前各种信息。
因此,它有相对好的应用前景。
目前叠前AVO反演方法主要有以下几种:(1)基于Powell算法的AVO非线性反演采用朱向阳和熊有伦提出的改进的Powell共扼方向算法,以Aki-Richard近似式为基础, 充分利用叠前地震数据丰富的振幅和旅行时信息,模拟平面波在层状弹性半空间传播时形成的地面反射记录,并使其与实际数据间的差异最小,从而获得地层的密度、纵波速度和横波速度分布。
(2)基于贝叶斯理论的AVO非线性反演基于贝叶斯理论, 结合似然函数与先验地质信息反演纵横波阻抗及密度。
先验模型参数的分布采用的是Huber分布。
Huber分布对于小的模型参数值进行一致性加权, 对于大的模型参数值采用拉普拉斯分布产生的权函数进行加权, 使之更能准确地反映模型参数的分布规律。
(3)点约束稀疏脉冲叠前反演基于贝叶斯参数估计的理论,假设似然函数服从高斯分布,待反演的参数服从改进的Cauchy分布,从而得到稀疏的反射稀疏序列,然后用已知点的纵波阻抗、横波阻抗和密度对反演结果进行点约束,从而使反演的结果更加准确可靠。
第6章 AVO技术详解AVO技术是利用反射系数随入射角变化的原理,在叠前道集上分析振幅随偏移距变化的规律,估求岩石的弹性参数、研究岩性、检测油气的重要技术。
AVO是振幅随偏移距变化(Amplitude Variation with Offset)的英文缩写或振幅与随偏移距关系(Amplitude Versus Offset) 的英文缩写,AVA 是振幅随入射角变化(Amplitude Variation with Incident Angle)的英文缩写。
在地震勘探中,共中心点道集记录的偏移距可以等价地用入射角表示,故AVO与AVA等价。
该技术自20世纪80年代提出以来,在油气勘探中不断发展,并得到迅速推广与广泛应用。
尤其是在天然气勘探中指导寻找天然气藏发挥了重要作用,对提高天然气勘探成功率受到了很好的效果。
从近几年的技术发展情况看,P波方位AVO已作为一种预测油气藏各向异性的有效方法而受到青睐。
6.1 AVO技术的理论基础根据地震波动力学中反射与透射的相关理论,反射系数(或振幅)随入射角的变化与分界面两侧介质的地质参数有关。
这一事实包含两层意思:一是不同的岩性参数组合,反射系数(或振幅)随入射角变化的特性不同,称为AVO正演方法;二是反射系数(或振幅)随入射角变化本身隐含了岩性参数的信息,利用AVO 关系可以反演岩石的密度、纵波速度与横波速度,称为AVO 反演方法。
6.1.1 Zoeppritz 方程AVO 技术的理论基础就是Zoeppritz 方程及其简化的思路。
设有两层水平各向同性介质,当地震纵波非垂直入射(即非零偏移距)时,在弹性分界面上会产生反射纵波、反射横波、透射纵波与透射横波,见图6—1。
各种波型之间的运动学关系服从斯奈尔定理22221111sin sin sin sin S P S P V V V V ϕθϕθ=== (6-1)图6—1 入射波、反射波与透射波的关系式中 1θ、1ϕ——纵波、横波的反射角;2θ、2ϕ——纵波、横波的透射角;1P V 、2P V ——反射界面上下介质的纵波速度;1S V 、2S V ——反射界面上下介质的横波速度。
个参数是以两个参数的比值,例如V p1V p211等形式出现。
这样就可以把V p1 £1等分别看作一第一节叠前流体检测技术近几年,随着地震采集处理技术的进步,尤其叠前偏移技术的发展和推广应用,使得研究人员可以得到来自地下真实反射点的叠前道集(CRP道集),为叠前烃类检测技术的发展奠定了资料基础。
目前基于叠前道集的直接烃类检测方法主要有两种:一种是在岩石物理建模的基础上进行叠前道集 AVO响应特征分析;一种是利用多个限角叠加数据体进行叠前弹性参数反演,利用纵横波波阻抗、纵横波速度比、泊松比、拉梅系数等敏感属性反映含油气性。
一.AVO分析技术1、 AVO理论简介AVO ( Amplitude Variation with Offset ),早先也称为 Amplitude Versus Offset,译为振幅随炮检距变化。
由此而衍生的有振幅随入射角变化 AVO( Amplitude Variation with Angle ),振幅随方位角变化 AVA( Amplitude Variation with Azimuth ),振幅随炮检距和方位角变化 AVOA( Amplitude Variation with Offset and Azimuth )等。
AVO作为一种含气砂岩的异常地球物理现象,最早在 20世纪80年代初被Ostrander发现。
这一现象表现为:当储层砂岩含气后,地震反射振幅随炮检距会发生明显的加大(基于SEG标准极性)。
因为AVO现象与含气砂岩的对应关系,从而引起勘探地球物理界广泛的重视。
后续的研究表明:这种异常现象并非一种特殊的形式,而是遵循Zoepprittz早先所提出的地震反射波动力学方程式,从而对AVO现象的解释有了完整的理论基础。
针对AVO现象继而出现的 AVO技术是继亮点之后又一项利用振幅信息研究岩性、检测油气的技术手段。
AVO技术具有以下特点:A、直接利用CDP道集资料进行分析,这就充分利用的多次覆盖得到的丰富的原始信息;B、禾U用振幅随炮检距(入射角)的变化的特点,即利用整条曲线的特点。
而亮点技术只是利用了这一特殊情况下曲线的一个数值。
所以,AVO技术对岩性的分析比亮点技术更为可靠。
C、这几年波动方程对地震剖面的成像有了更大的成果,是对地下构造形态的反演。
AVO技术从严格意义上说算不上是利用波动方程进行岩性反演分析的方法,但是其理论和思路是对波动方程得到的结果的比较精确的利用。
D、 AVO技术是一种研究岩性的比较细致的方法,并且需要有测井资料的配合。
2、 AVO技术的理论基础振幅随炮检距的变化来自于所谓的“能量分区”。
当地震波入射到地层界面时,一部分能量反射,一部分能量透射。
如果入射角不等于零度,纵波( P波)能量一部分反射,一部分转化成透射P波和S波。
反射和透射波的振幅能量取决于地层边界的物理性质差异。
纵波速度Vp、横波速度Vs和密度p是非常重要的。
同时,需要注意反射振幅也依赖于入射波的入射角(图5-3-1 )。
因此,当一个平面纵波非垂直入射到两种介质的分界面上,就要产生反射纵、横波和透射纵、横波。
在界面上,根据应力连续性和位移连续性,依据边界条件并引入反射系数、透射系数,就可以得出四个相应波的位移振幅应当满足的方程叫做Zoeprritz方程,这个方程是 Zoeprritz在1919年解出的。
这个方程组比较复杂,不能解出新产生的波的振幅与有关参数明确的函数关系。
但是从方程组可以看出,一般反射纵波的反射系数Rpp是入射角界面上部介质的密度p,纵波速度Vp1,横波速度VS1以及界面以下的介质密度 p 2,纵波速度VP2,横波速度VS2等七个参数的函数,可以简单的表示为Rpp ( &, VP1, VS1, p 1, Vp2, VS2, p 2),虽然不能直接从方程中解出Rpp与七个参数的具体关系,但是可以假设以物质的六个物性参数为参变量,以&为变量,仔细分析可以得到,六图5-3-1 :入射到地层边界的地震波示意图通过对地层弹性参数的研究,得到泊松比 6是一个对岩性和含油气情况反应比较敏感,参数,所以就要对方程进行适当的化简, 得出以泊松比6为参数的以6为变量的简单的近似关系,即Rpp=f(&,6,这样通过反算 6来达到对储层的参数测定和检测。
由于Zoeppritz 方程过于复杂,因此有许多学者尝试对其进行简化或近似,其中比较著名和实用 的主要有: 1) Koefoed 的试算Koefoed 在1955年第一个给出了用 Zoeppritz 方程,以R 为参数计算出的 Rpp 〜&曲线。
他用 17组纵横波速度、密度和泊松比参数,较为详细地研究了泊松比对两个各向同性介质之间反射和折 射面所产生的反射系数的影响,最大的入射角达到30°。
他的研究结果被公认为Koefoed 五原则。
虽然Koefoed 的结论说明了利用Rpp 〜&曲线是可以反算出泊松比6的,但是用未简化的 Zoeppritz方程进行计算太复杂,因此反求弹性参数也是很复杂的。
2) AkiKI 和RichardsPG 的纵波近似公式1980年AkiKI 和RichardsPG 给出了一个纵波反射振幅的一个近似表达式:2 2 ■ 21 “ Ns • 2T拧 sec 二 Vp ,V s • 2TV s 、=(1-42 sin 4 2 sin )2 V p2;2 V p V p 2 Vs这个公式是有适用条件的,就是界面两边的弹性介质性质的百分比变化小,式中的参数有: A/p=V p2-V p1、△Vs=V s2-V s1、Vp=(V p2+V p1)/2、Vs=(V s2+V s1)/2、Ap = p2- p 1、p = (p 2+p 1) /2、0= ( 0+ 02) /2, 0是入射角,0是按斯奈尔定理计算的透射角。
这个式子虽然对 Zoeppritz 方程进行了化简,但是用了 Vs 做参数而未用 6作参数。
3) Bortfeld 近似式 1961年,Bortfeld 提出近似式个参数,再加上在同一种介质中,纵波速度 v p ,横波速度V S,以及泊松比 6之间又有关系,如:于是有关系式:Rpp=f (&,,卫,这样来达到减少参数的目的。
从V p 2P 2理论上说,在实际地震记录上得到某个界面的反射波的振幅与入射角的变化关系曲线,并且又知道 某些参数,就可以利用曲线族作为量板来估算地层参数。
R ( 0) (5-3-1)4) Hilterman 近似式V p2T 2COS 円一V p1C cos * V p2 ?2 COSK V p1 q cos —'si n d) 一 Vs2R — V s2P 21+ k Ws1 -V s2 3(V s1—V s2)+2 s2 1 s2 2(5-3-3)I Vp1 丿-M + 12」5) Shuey 简化公式1985年Shuey 利用上面的式子进行了简化,用 代替Vs,修改得:R (e)=R o + |A o R o+—”2"sin 2。
+l^yP(tan 2 日— sin 2日) (5-3-4)(1 —<! ) 一 2 Vp.c - c - c 1、厂-「2.「 2 。
1 ^A Vp 也 P "垂直入射时的反射振幅 R 。
= ----------- +—— ,d 、(2分别为入射介质和透射介质的泊松比。
2JVp P 丿其中,界面两侧泊松比的差△(是一个至关重要的因素,这就是振幅与炮检距关系研究的物理基础。
Shuey 近似式的特点就是其三项都有明确的物理意义: 第一,垂直入射时,01= 62=0=0, R ( 0) =R,即R 0是垂直入射时的反射振幅。
第二,在中等入射的情况 (0〈 0〈 30°),有近似tan 銓sin 0,于是有:r ( “ 、2 1只(日)=巳 1+ A 0+ ---------- 3— bin 2日」(5-3-5)L v (1-丿此时,反射振幅与A 有关,前两项起作用。
这时的反射系数与介质的泊松比有密切关系, 因此,利用此式更能突出油气特征。
第三,对于大角度入射情况,反射振幅与速度变化有关中的第三项起主要作用。
6) Simth 和Gidlow 加权叠加方法1987年Smith 和Gidlow 根据Gardner 等人给出的水饱和岩石的密度与速度的 4次方根成正比的假设,将 Aki 和Richards (1980)提出的Zoeppritz 方程进行了修改化简:△V pzR pr …B i 匕 Ci s(5-3-6 )V pVsRpp=1ln(V p 2认却2 V pi A cos 日2Vp i--V s22) 2 +Rpp =其中:Vs 2 =Vp 2— —2(1 )△Vp Vp• : Vp Vp - - A 二 B 一2 1 B 上2r1 _aP(5-3-2)纵横波速度,B 、C 为常数。
因此,尽管表示反射系数的简化式有很多,但其最终目的都是为了求得一个最简表达式来表示 反射系数随入射角的变化。
而上、下介质的泊松比又对反射系数随入射角的变化起重要作用,因此, 需要在表达式中包含标志油气特征的参数。
Hilterman 的简化方程与Shuey 的简化方程满足了以上要求,且适用于各种层状模型的 AVO 模拟,因此得到了广泛的应用。
在应用中最常见的一种简化公式 为:2RU ) =P Gsin 2 二(5-3-7 )式中:R 为反射系数、B 为入射角;P 为近似为零偏移距下的纵波的反射振幅,也称 AVO 截距, 其大小决定于上下层之间的纵波波阻抗差异;G 为纵波反射振幅随入射角的变化梯度,也称AVO 斜率,其大小取决于泊松比的变化。
这个方程不能很好处理入射角较大的情况,但是简单易懂,至今 仍然非常实用。
对 NMO (正常时差校正)道集中每一个时间采样点作 sin20 ( B 为入射角)和振幅交汇图(图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 反映的是泊松比的特征。