地震迭前反演技术印兴耀课件

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,二级,三级,四级,五级,最新课件,*,叠前地震反演,印兴耀,中国石油大学,山东东营,叠前地震反演印兴耀中国石油大学,反演,叠前反演,叠后反演,EI,反演,AVO,反演,利用叠后资料在井的约束下得到,AI,利用叠前资料获取,AVO,属性,利用叠前资料,及弹性波阻抗获取岩石弹性参数,反演分类,反演叠前反演叠后反演EI反演AVO反演利用叠后资料在井的约,水平叠加,角道集叠加,(8,度,),角道集叠加,(16,度,),角道集叠加,(24,度,),水平叠加与角度部分叠加,水平叠加角道集叠加(8度)角道集叠加(16度)角道集叠加(2,页岩、含油、气、水砂岩的速度随深度变化关系,页岩、含油、气、水砂岩的速度随深度变化关系,只利用纵波速度，难以区分砂岩和泥岩。在泊松比和纵波速度交汇图中，砂泥岩、和油气可以区分开。,泊松比与岩性及流体成分的关系,Vp,Vp,只利用纵波速度，难以区分砂岩和泥岩。在泊松比和纵波速,1,、弹性模量,岩石弹性参数及其相互关系,1、弹性模量 岩石弹性参数及其相互关系,2,、纵波速度和横波速度,岩石弹性参数及其相互关系,2、纵波速度和横波速度岩石弹性参数及其相互关系,3,、泊松比,岩石弹性参数及其相互关系,泊松比的物理意义,3、泊松比岩石弹性参数及其相互关系 泊松比的物理意义,是流体的体积模量，所以可以叫做流体的不可压缩性度量（,fluid incompressibility,）。,是抵抗剪切应变的度量，它常被称为刚性模量,（,modulus of rigidity,），不可压缩性的度量,（,imcompressibility,）或剪切模量（,shear modulus,）。,(,泊松比,),：,是反映岩性与含油气性的重要参数，它是用岩石纵向拉伸和横向压缩的比值来表示,常用弹性参数的意义,是流体的体积模量，所以可以叫做流体的不可压缩性度量（,含油气地层一般都具有较低的,和,/,与流体无关，仅反映岩石骨架的性质,和,能够灵敏地反映储层属性，并且物理意义明确，在岩性和流体预测方面具有重大的意义,和,含油气地层一般都具有较低的和/和,描述岩石弹性特征的主要地震参数有岩石的弹性模量,、,密度 、纵波速度 和横波速度 等。弹性模量反映岩石应力应变关系的特征，密度反映岩石的比重，速度则反映地震波在岩石中的传播的特征，它是弹性模量的函数,。,拉梅常数,剪切模量,体积模量，,描述岩石弹性特征的主要地震参数有岩石的弹性模量、密度 、,地震迭前反演技术印兴耀课件,Vp,Vs,LamdaRho,MiuRho,井的砂泥岩纵横波速度交汇图,LamdaRho,和,MiuRho,交汇图,砂岩,泥岩,砂岩,泥岩,0.4,0.8,VpVsLamdaRhoMiuRho井的砂泥岩纵横波速度交汇,Vp,Vs,交会图,交会图,Vp Vs交会图 交会图,/,交汇图,在,/,交会图中,可以看出,，油层（红色）基本都分布在低泊松比、低,/,的区域。说明该区油层的特点之一是比水层具有更低的泊松比。当泊松比大于,0.32,时，砂岩为水层。,/交汇图在 /交会图中可以看出，油层（红,水层,油水同层,油层,水层,砂岩,泥岩,砂岩 泥岩,Nml,砂泥岩交会图,砂岩,泥岩,Nml 砂泥岩交会图砂岩 泥岩,Nml,砂泥岩交会图,砂岩,泥岩,Nml 砂泥岩交会图砂岩 泥岩,Nml,砂泥岩交会图,砂岩,泥岩,Nml 砂泥岩交会图砂岩 泥岩,亮点技术的出现，使地震烃类检测能力有了很大提高，但实践表明，该技术也有其局限性，并不是所有的亮点都与气层有关，某些特殊的岩性体也可以形成地震剖面上的“亮点”,八十年代初，,Ostrander,将反射系数随入射角变化应用于亮点型含气砂岩的识别，发现含气砂岩反射振幅随偏移距增加而增大，含水砂岩反射振幅随偏移距增加而减小，这一发现极大地改善了烃类检测的能力,AVO,分析是一项利用振幅随偏移距变化特征来分析岩性和油气藏的地震勘探技术，其基本思想是振幅随入射角变化(,AVA),的规律与岩性参数有关。,AVO,也有其陷阱，并不是所有气层的,AVO,特征都呈增加现象，特殊岩性体也会引起,AVO,增加现象,亮点技术的出现，使地震烃类检测能力有了很大提高，但实践表明，,充分利用地震反射振幅随偏移距不同或入射角不同而变化的特征，通过对这些数据的纵横波反演技术研究，能得到,纵横波波阻抗，纵横波速度，密度，泊松比,等多种参数；提供了研究岩性和流体变化规律的更多、更有效的反演数据体成果,叠前地震反演较叠后反演推进一步,叠前反演,充分利用地震反射振幅随偏移距不同或入射角不同而变化的特征，通,从含有丰富地质信息的叠前资料入手，深入研究叠前记录的正演模拟、叠前波形反演技术，采用弹性波反演方法由叠前道集中获取地层纵横波速度、密度等参数；然后利用不同角度的角道集数据得到弹性波阻抗数据体。这些技术相互补充，使叠前地震技术的研究进一步得到深入,叠前弹性波反演和,EI,有机结合，可以建立一套较为完善的联合反演技术，从中得到的地层速度、密度、泊松比、弹性模量等岩性参数数据体，将为研究储层的变化规律提供新的手段,从含有丰富地质信息的叠前资料入手，深入研究叠前记录的正演模,P,垂直入射反射,Zoeppritz,方程及其近似,P垂直入射反射Zoeppritz方程及其近似,Zoeppritz,方程,Zoeppritz方程,当地震波垂直入射到界面上时有,按斯奈尔定理有：,由,Zoeppritz,方程解得,当地震波垂直入射到界面上时有按斯奈尔定理有：由Zoeppri,简单模型,简单模型,合成,CMP,道集，共60道，最大炮检距1500米,合成CMP道集，共60道，最大炮检距1500米,近偏移距、中偏移距和远偏移距各20道的叠加结果,反射2缺失,(第二和第三层阻抗相等),反射3微弱,近偏移距、中偏移距和远偏移距各20道的叠加结果反射2缺失反射,某区典型含油砂岩的地震反射道集,某区典型含水砂岩的地震反射道集,泥岩,含油砂岩,泥岩,含水砂岩,含油性研究,Vp,、,Vs,、,某区典型含油砂岩的地震反射道集某区典型含水砂岩的地震反射道集,叠前弹性波,波形反演,弹性参数,叠前地震数据,角道集数据,弹性波阻抗,约 束,纵横波速度、密度、泊松比、,弹性模量等地层弹性参数,叠前弹性波,波阻抗反演,叠前弹性波弹性参数叠前地震数据角道集数据弹性波阻抗约 束纵横,叠前波动方程反演,算法基本流程,随机产生一族弹性地层模型,计算每一地层模型的合成地震道集,合成道集与实际观测数据的比较,遗传,算法,修改,地层模型,收敛？,否,是,计算,适应度值,保存结果,结束,叠前波动方程反演算法基本流程随机产生一族弹性地层模型计算每一,叠前弹性波反演,叠前反演首先需要利用波场模拟合成叠前记录，选择合适的正演模拟方法是叠前反演中极为重要的一步,射线追踪法运算速度较快，但地层数目很多时，对其进行射线追踪极其困难,有限差分或有限元方法虽然精确，但运算量极大,反射率法（谱方法或慢度法）用精确的波动方程方法进行波场模拟，能够计算包括,转换波,和,多次反射,在内的地震响应，运算速度介于有限差分和射线追踪之间,叠前弹性波反演叠前反演首先需要利用波场模拟合成叠前记录，选择,叠前弹性波反演,叠前地震反演能够从地震数据中得到,P,波和,S,波信息。由于,P,波对孔隙流体的变化较敏感，而,S,波主要与岩石骨架有关，流体的变化对其影响较小。为了检测流体类型，需要同时提供岩石的,P,波和,S,波信息,在叠前波形反演的文献中，绝大多数是平面波反演，而实际地震记录是球面波记录,基于声波理论的反演，得不到横波速度,局部线性化反演方法需要地层参数的导数构成雅可比矩阵，而对地下进行精确模拟需要几百层地层，求导非常困难,Why?,叠前弹性波反演 叠前地震反演能够从地震数据中得到P波和S波信,叠前弹性波反演,从左向右分别为初始群体最佳个体的,P,波速度、,S,波速度、密度及其叠前合成记录。图中,红色曲线为真实井曲线,，蓝色曲线为初始估计模型，可见初始模型参数与井数据差别较大。图中中间的黑色细线是地层参数的低频趋势，两侧的黑线是参数的约束范围。,叠前弹性波反演从左向右分别为初始群体最佳个体的P波速度、S波,叠前弹性波反演,初始群体最佳个体的叠前合成记录(,b),与井曲线叠前道集(,a),的对比，可见二者存在很大差别。,(,a),(,b),叠前弹性波反演初始群体最佳个体的叠前合成记录(b)与井曲线叠,从左向右分别为反演的,P,波速度、,S,波速度、密度及其合成记录。图中红色曲线为真实井数据，蓝色曲线为反演得到的地层参数，可见反演结果与真实模型有很好的吻合。,叠前弹性波反演,从左向右分别为反演的P波速度、S波速度、密度及其合成记录。图,从残差剖面上看出，两个地震道集有很好的相似性,叠前弹性波反演,从残差剖面上看出，两个地震道集有很好的相似性叠前弹性波反,叠前弹性波反演,时间(,s),蓝色曲线是实际测井曲线,P,波速度,S,波速度密度,叠前弹性波反演时间(s)蓝色曲线是实际测井曲线 P波速度,Vp Vs,密度,（,蓝色曲线为实际井数据，红色曲线为叠前波形反演结果，纵轴表示深度，单位：米,）,Vp Vs,从实际井曲线与反演结果的对比来看，,P-,波速度和,S-,波速度匹配的很好，密度曲线匹配的较好,密度曲线能较好地反映地层的变化，含气砂岩的密度较低,与泥岩相比，砂岩表现高,P-,波速度和高,S-,波速度。含气后，,P-,波速度降低，而,S-,波速度基本不变。充分表明,P-,波速度的降低是含油气所致而不是由岩性变化引起,密度曲线(深度域),从实际井曲线与反演结果的对比来看，P-波速度和S-波速度匹配,合成道集,(,左,),与实际地震记录,(,右,),对,LS36-1-1,井附近的,CDP,道集进行了遗传算法反演，下图分别是合成道集与实际地震记录的比较，箭头所示的同相轴同样显示了气砂岩反射引起的振幅随偏移距增大现象,合成道集(左)与实际地震记录(右) 对LS36-1-1井附近,反演结果与井曲线的对比，左图是纵波速度曲线，右图是密度曲线；其中红色是反演曲线，兰色是实际井曲线,经过对比认为，反演结果是可靠的，从而为弹性波阻抗反演提供了先验信息。,反演结果与井曲线的对比，左图是纵波速度曲线，右图是密度曲线；,地震迭前反演技术印兴耀课件,地震迭前反演技术印兴耀课件,合成的叠前记录,LS36-1-2,井附近的实际,CMP,道集,合成的叠前记录LS36-1-2井附近的实际CMP道集,根据,AVO,理论，零炮检距（或小炮检距）剖面可以近似为声阻抗,AI,（,Acoustic Impedance,）的函数，它与岩石的密度和纵波速度有关。为了充分利用大炮检距地震振幅信息，需要利用一种,与入射角有关的弹性阻抗,EI,（,Elastic Impedance,），,这种弹性阻抗形式上,类似于声阻抗,AI,，,不但与岩石的密度和纵波速度，而且与,横波速度、入射角,有关,。,弹性波阻抗,(EI),根据AVO理论，零炮检距（或小炮检距）剖面可,EI,最早出现在,Connolly,于,1999,年发表的论文中，随后在,2000,年的,SEG,年会上同时出现了,4,篇论文对,EI,进行了研究，,ARCO,公司介织了他们申请了专利的弹性波阻抗反演方法，认为在求取的反射系数的稳定性方面要好于,Conno11y,方法，而且计算的,EI,和,AI,数值在同一个尺度下：,CGG,公司发表论文认为,EI,反演效果要好于,AV0,反演，因为,EI,在抗噪能力方面比,AV0,方法具有优势：同时,BP Amoco,公司在会上又提出了扩充弹性波阻抗方法，可以用于流体和岩性预测,这些进展说明：弹性波阻抗已经成为波阻抗反演进一步发展的方向之一，地震反演的发展正在走向,AI,和,EI,相结合、,AI,和,AVO,相结合的道路。,EI最早出现在Connolly于1999年发表的论文中，随后,Zoeppritz,方程及其近似,P-P、,P-S,反射和透射,P,P,SV,SV,P,Zoeppritz方程及其近似P-P、 P-S 反射和透射P,Zoeppritz,方程及其近似,Aki, Richards,和,Frasier,近似,第一项包含,P-,波速度,第二项包含密度,第三项包含,S,波速度,Zoeppritz方程及其近似Aki, Richards和F,Shuey,近似,Zoeppritz,方程及其近似,Hilterman,通过设下面的假设进一步简化了,Shuey,方程,(1) 仅用前两项,(2) 设,Shuey近似Zoeppritz方程及其近似Hilterma,Zoeppritz,方程及其近似,Smith,和,Gidlow,近似,用这一近似在地震道集上做加权叠加可得到储层岩石性质的信息,Zoeppritz方程及其近似Smith和Gidlow近似用,其中,P, 垂直入射时的纵波反射系数,G, 反射系数函数的梯度,Zoeppritz,方程的近似均可统一为,标准的斜率截距式,AVO,Zoeppritz,方程及其近似,从中可以得到许多属性剖面,其中P 垂直入射时的纵波反射系数Zoeppritz,依据的,Zeoppritz,近似方程,Aki_Richards(1980),近似,Fatti(1994),近似,Gray(1999),近似,依据的Zeoppritz近似方程Aki_Richards(1,弹性波阻抗反演,弹性波阻抗,Elastic Impedance,地震波垂直入射时的褶积模型，对于与角度相关的数据，有：,S(),R()*W(),S(),为角度地震道，,R(),是与角度有关的反射系数，它是由,P,波速度、,S,波速度和密度，根据,Zoeppritz,方程计算得到的；,W(),是角度子波，由角度地震数据得到。就像对反射系数进行积分而得到声阻抗,(,Acoustic Impedance),一样，角度反射系数用来计算弹性波阻抗。,EI,并不是一个可以进行物理测量的属性，它是由,Zoeppritz,方程或其近似式得出的用来解释地震数据的属性，弹性波阻抗反演是目前为解决岩性解释而提出的比较适用的方法 。,弹性波阻抗反演 弹性波阻抗地震波垂直入射时的褶积模型，对于与,弹性波阻抗反演,寻求一个函数，使其与声阻抗有相似的性质，以便于以任意角入射时反射系数均可以用下面的公式表示：,函数 称作,EI,（弹性波阻抗），当地层弹性参数变化不大时，用,EI,的对数值来表示反射系数是准确的,弹性波阻抗反演 寻求一个函数，使其与声阻抗有相似的性质，以便,弹性波阻抗反演,去掉两边的微分和对数项，得到下式：,该式即为,P,波的弹性波阻抗,弹性波阻抗反演 去掉两边的微分和对数项，得到下式： 该式即为,弹性波阻抗反演,弹性波阻抗反演流程,用测井曲线合成横波速度，进而计算井中的,EI,曲线，对于同一个,CMP,道集来说，,EI,是入射角的函数,将偏移距数据转化为角道集数据，,估算每个角度数据所对应的子波,利用与迭后反演相类似的反演算法计算出各个角度的弹性波阻抗剖面,弹性波阻抗反演 弹性波阻抗反演流程,波阻抗反演,弹性波阻抗反演,波阻抗反演弹性波阻抗反演,子波提取,子波的振幅和频率特性随角度的变化而变化，这保证了反演处理的稳定性。各角度的子波在形状上很相似，振幅和相位也都比较一致，这有利于各角度资料反演效果的对比。,子波提取 子波的振幅和频率特性随角度的变化而变化，这,在不同角度的道集上进行合成记录标定,5-15,度的,y551,井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线,5-15,度,/y551,井,5-15,度合成记录标定,在不同角度的道集上进行合成记录标定5-15度的y551井的合,15-25,度的,y551,井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线,15-25,度,/y551,井,15-25,度合成记录标定,15-25度的y551井的合成记录层位标定图，从左至右分别为,25-35,度的,y551,井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线,25-35,度,/y551,井,25-35,度合成记录标定,25-35度的y551井的合成记录层位标定图，从左至右分别为,三、弹性阻抗反演,角道集资料,5-15,度角道集剖面,Ang5-15,度,三、弹性阻抗反演角道集资料5-15度角道集剖面 Ang5-1,15-25,度角道集剖面,Ang15-25,度,15-25度角道集剖面 Ang15-25度,25-35,度角道集剖面,Ang25-35,度,25-35度角道集剖面 Ang25-35度,四、弹性阻抗反演,Ang5-15,度,Complete_impedance_5-15,度,四、弹性阻抗反演Ang5-15度Complete_imped,“,AVO,反演”，是从地震数据中抽取,AVO,属性，然后在合适的测井曲线的约束下将其转化。 “弹性波阻抗”反演，首先是在模拟的测井曲线约束下转化一定角度范围的角道集，再抽取所需的弹性参数。这两种方法只在子波不随偏移距变化时是相同的，当子波由于噪音等因素发生变化时，会在不同的,AVO,属性之间产生“泄漏”，在通常使用的入射角范围内，这种“泄漏”压制了一些抗噪属性。这使得,AV0,反演不能提供满意的解决办法，然而，弹性波阻抗能够避免这个问题，因此弹性波阻抗具有优越性。,AVO,与,EI,反演,“AVO反演”，是从地震数据中抽取AVO属性，然后在合适,I,P,R,P,(,P-P,),声阻抗,I,P,R,P,q,(,P-P,),弹性波阻抗,I,P,R,S,q,(,P-S,),弹性波阻抗,I,S,R,S,(,S-S,),剪切阻抗,IPRP(P-P)声阻抗IPRPq(P-P)弹性波阻抗IPR,角道集叠加过程可以简单描述为：把固定偏移距的记录转换成固定入射角,(,或一定角度范围内的叠加,),记录，将振幅随偏移距变化的关系转换为振幅随入射角变化的关系,计算时间,T()j,对应的炮检距,X()j,在点,(T()j-1,，,X()j-1),( T()j,，,X()j),之间建立线性关系,大角度对应时间,T(min),小角度对应时间,T(max),提取各道的,T(min),到,T(max),之间的振幅值做角度部分叠加,对应于角度,(i),的角度道,A(),处理后的,CMP,道集资料,确定,CMP,道集的中心角,序列,(1),、,(2),、,(i),确定中心角,(i),的,范围,(i)min,到,(i)max,计算层速度,V,int,均方根速度,V,rms,角道集部分叠加,角道集叠加过程可以简单描述为：把固定偏移距的记录转换成固,角度道数据沿图中曲线拾取。可以看出，一般情况下，同一时间（深度）处远偏移距对应大角度，小偏移距对应小角度；同一偏移距角度随深度减小。,零,偏,移,距,双,程,旅,行,时,(ms),偏移距,(,m,),偏移距和角度的关系,角度道数据沿图中曲线拾取。可以看出，一般情况下，同一时间（深,角叠加剖面可以为弹性阻抗反演提供资料,角度的估算是一个复杂的过程，其准确性受地震资料前期处理的影响，为了或得更准确的角道集叠加资料，应加强地震资料的前期处理，同时在估算角度方面应采用更复杂的数学模型,角部分叠加资料比叠后资料包含了更丰富的信息；,AVO,和弹性波阻抗反演要求地震资料具有,AVA,特征，因此将,CMP,道集资料转化为角道集部分叠加资料是非常必要的,地震迭前反演技术印兴耀课件,弹性波阻抗的近似算法有多种，在此我们用,BP,近似来计算弹性波阻抗，根据,Connolly,在1999年提出的弹性波阻抗近似公式可以得到弹性波阻抗的如下表达式：,进而由上式可以得到如下表达式：,对于同一角度来说,A，B，C,均为常数：,从弹性波阻抗体中提取岩性参数体,弹性波阻抗的近似算法有多种，在此我们用BP近似来计算弹性波阻,对于多个不同的角度而言有：,这是一个典型的,型方程组。写成矩阵形式：,上面表达式中：左边第一项是一个,n,3的常数矩阵；第二项是所要求的纵、横波速度及密度的对数形式；右边是一个,n,1的常数矩阵，它是由反演得到的各个角度的弹性波阻抗体分别取对数得到的。,从弹性波阻抗体中提取岩性参数体,对于多个不同的角度而言有：这是一个典型的型方程组。写成矩阵形,通过最小平方法得到各样点的,可以进一步得到：,Ip、Is,、泊松比、纵波模量、横波模量及拉梅参数等。,对其分别取指数就得到各采样点的岩性参数体：,从弹性波阻抗体中提取岩性参数体,通过最小平方法得到各样点的可以进一步得到：Ip、Is、泊松比,参数,Lamda*Rho,和,Mu*Rho,的获取,由于,Lamda*Rho,和,Mu*Rho,能够灵敏地反映储层属性，并且物理意义明确，在岩性和流体预测方面具有重大的意义，因此我们有必要从反演中拾取这两个属性。有两种方法可以获取,Lambda*Rho,和,Mu*Rho,，在此称其为直接法和间接法,间接法,在,Connolly,（,1999,）方程的基础上，采用回归算法求取,Alpha,、,Beta,和,Rho,，根据各岩性参数间的关系进一步求得,Ip,、,Is,、,Lambda*Rho,和,Mu*Rho,直接法,用,Ip,和,Is,及,Lambda*Rho,和,Mu*Rho,的显式近似方程,Fatti,近似和,Gray,近似来求取它们，这种方法是直接得到,Ip,和,Is,及,Lambda*Rho,和,Mu*Rho,参数Lamda*Rho和Mu*Rho的获取由于Lamda*R,感谢亲观看此幻灯片，此课件部分内容来源于网络，,如有侵权请及时联系我们删除，谢谢配合！,77,感谢亲观看此幻灯片，此课件部分内容来源于网络，77,