叠前深度偏移、速度建模、保幅偏移.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,三 维 叠 前 深 度 偏 移,叠前深度偏移成像因稳健且效率与精度高，已成为人们关注的重点和焦点。包括射线法和波动方程法（,FXFD,SSF,FFD,GS)。,以,粘滞声波、弹性波,波动方程为基础的偏移方法，以及适于,各向异性,介质的偏移方法都得到了一定的应用，并取得了较好的效果。,射线法三维叠前深度偏移,Kirchhoff,积分法的关键是绕射旅行时的计算,Kirchhoff,积分法叠前深度偏移成像公式,（,1,）,绕射旅行时计算方法,可以采用变速射线追踪法、基于费马原理的二维有限差分法、和稳健高效的三维迎风有限差分法,常规方法,存在不足。,三 维 叠 前 深 度 偏 移,波动方程法叠前深度偏移,基于共炮集的波动方程叠前深度偏移概论,偏移思路,引入基于单程波方程的波场传播算子,正向,延拓,反向延拓,三 维 叠 前 深 度偏 移,偏移结果的表示,共炮集数据的叠前深度偏移成像公式：,（,3,）,（,2,）,可见,成像的核心,是波场外推算子，目前基于波动方程的波场延拓算子不外乎：有限差分波场延拓算子和付立叶波场延拓算子。,分析：,三 维 叠 前 深 度 偏 移,频率空间域有限差分,(,FXFD),法叠前深度偏移,基本理论,基于频率波数域波动方程的分解，绕射项解决,绕射波收敛,问题，折射项,校正由于横向变速引起的时差,。因此,该方法不仅能对陡倾角进行成像，而且能适应速度的任意横向变化，故该算法可使复杂地质体精确成像。,（,4,）,三 维 叠 前 深 度 偏 移,（,7,）,（,6,）,（,5,）,相关成像条件,（,8,）,三 维 叠 前 深 度 偏 移,Marmousi,模型及其,FXFD,法偏移剖面,模型试算,下面给出,Marmousi,模型,FXFD,的试算结果,三 维 叠 前 深 度 偏 移,分步傅立叶,(,SSF),法波动方程叠前深度偏移,基本原理,（,10,）,（,9,）,SSF,法叠前深度偏移，把速度场分解为常速背景和变速扰动两部分，对常速背景在频率波数域采用相移处理，对层内的变速扰动，在频率空间域采用时移校正（第二次相移）,（,12,）,（,11,）,（,13,）,（,14,）,试算情况,分别采用雷克子波和提取的子波模拟震源波场时的,Marmousi,模型,SSF,偏移剖面,三 维 叠 前 深 度 偏 移,三 维 叠 前 深 度 偏 移,傅立叶有限差分（,FFD,）,法波动方程叠前深度偏移,（,15,）,（,16,）,三 维 叠 前 深 度 偏 移,外推公式,（,17,）,（,19,）,（,18,）,FFD,法与,SSF,法叠前深度偏移的精度比较,三 维 叠 前 深 度 偏 移,三 维 叠 前 深 度 偏 移,广义屏法波动方程叠前深度偏移,基本理论,广义屏偏移算子是基于波的散射理论，从波动方程,Green,函数解出发，借助,Born,近似、屏近似等一系列数学手段而导出的,。,广义屏偏移算子,（,21,）,（,20,）,三 维 叠 前 深 度 偏 移,（,24,）,（,25,）,（,23,）,（,22,）,x,z,screen,三 维 叠 前 深 度 偏 移,广义屏算子的小角度近似及相屏算子,扩展的局部,Born,近似的广义屏算子,（,4.26,）,（,4.27,）,（,4.30,）,（,4.28,）,（,4.29,）,三 维 叠 前 深 度 偏 移,还有,MRS-ELBF,的广义屏算子、,ELRF,的广义屏算子及,基于,Born/,Rytov,近似的联合广义屏算子。,模型试算,Marmousi,模型相屏偏移剖面,Marmousi,模型,ELBF,法偏移剖面（,dz,=2m,）,a.MRS-ELBF,法，,b.,最优,Born,近似法，,c.Born/,Rytov,联合法,三 维 叠 前 深 度 偏 移,三 维 叠 前 深 度 偏 移,弹性波波动方程叠前深度偏移,方法原理,由全弹性波波动方程出发计算合成数据的垂直分量和水平分量，然后从这些分量中提取出纯的,P,波和,SV,波，然后分别用所提取的纯,P,波和,S,波分别进行波动方程偏移，可以得到两个深度偏移剖面,这两个深度偏移剖面,从理论上讲应该相等,因此可以相互比较,用来检验偏移方法和偏移结果的实用性和可行性,。,模型试算,三 维 叠 前 深 度 偏 移,概述,一般的说来，不管是时间偏移还是深度偏移都存在一个,矛盾,，即偏移的目的是为了求得地下地层的真实结构和形态，但为了进行偏移，我们必须事先给出一个速度模型。即出现了“已知答案，求解问题”的矛盾,。为了更好地进行偏移成像，必须建立正确合理的速度模型。,分类：理论：射线法、波动方程法、层析成像法,道集,：,CMP,CIP,CRP,CFP,CRS,等,偏 移 速 度 建 模,常用：速度谱分析法、相干反演法和层析成像法,射线法偏移速度建模,基本原理,基于,射线理论,和,最大叠加能量,或最大相关系,数准则，,在,CMP,道集通过,层剥离,实现偏移速,度建模,校正偏移速度，这种方法计算,效率高，,但存在速度,误差累计,和由于在,CMP,道集进,行速度分析所带来的,多解性,问题。,偏 移 速 度 建 模,模型试算,偏 移 速 度 建 模,波动方程法偏移速度建模,共成像点道集偏移速度建模,基本原理,基于,Kirchhoff,积分法和波动方程法叠前深,度偏移和,CIP,道集,剩余曲率分析（,RCA,）,偏 移 速 度 建 模,偏 移 速 度 建 模,FFD,法,偏 移 速 度 建 模,SSF,法,速度修正准则和迭代速度反演,模型试算及分析,偏 移 速 度 建 模,1000,0,200,150,100,50,0,2000,(,d),h,/m,TRACE,/(,25,m),(,c),1000,0,0,200,2000,150,100,50,TRACE,/(,25,m),h,/m,0,100,200,1000,300,400,2000,0,(,b),CDP,/(,12.5,m),h,/m,0,1000,2000,100,200,400,300,(,a),CDP,/(,12.5,m),h,/m,共聚焦点,CFP,道集偏移速度建模,基本原理,该方法基于等时原理、,CFP,偏移和射线追踪,/,波,动方程正演模拟,，,借助逆时聚焦算子和,CFP,响,应之间的差异，即通过,DTS,分析来实现偏移速度,建模,偏 移 速 度 建 模,模型试算及分析,偏 移 速 度 建 模,偏 移 速 度 建 模,共,反射,点,CRP,道集偏移速度建模,CRP,道集,偏移速度建模是基于叠前深度偏移、地震层析成像（,ST,）、,和深度聚集分析（,DFA,）。,在,DFA,中，若偏移深度和聚焦深度之间的误差为零，则偏移速度合理；否则，依据误差修正偏移速度,。,偏 移 速 度 建 模,共,反射,点,CRP,道集偏移速度建模,深度聚集分析（,DFA,）,基于实际速度和过大速度的向下延拓示意图,（来自,Jeannot,et al.,1986,）,偏 移 速 度 建 模,共,反射面,CRS,道集偏移速度建模,CRS,道集偏移速度建模是基于射线理论、三参数优化、和CRS,叠加,。,效率较高，速度分析精度取决于优化三参数的精度，适于较复杂地质体的速度建模.,CRS,道集偏移速度建模就是利用优化后的地震三参数来实现的。,地震层析成像(,ST),法偏移速度建模,地震层析成像法是在最小平方意义上，通过迭代反演修正速度。一般的，射线法和波动方程法与速度谱分析方法一样简单实用，而层析成像法主要用来对速度异常区成像，解决较复杂地区的速度建模。,偏 移 速 度 建 模,偏 移 速 度 建 模,小结,除了上述的几种速度建模方法以外，随着转换波受人们重视程度的增加，,转换波,偏移速度分析方法应运而生。由于转换波路径的不对称性，处理起来比,P,波复杂的多，表现出对速度各向异性敏感且衰减快。因此转换波偏移速度建模必须考虑到各向异性及相干噪声的影响。,波 动 方 程 保 幅 偏 移,波动方程法三维叠前深度偏移技术在很大程度上能够实现保幅处理，特别是付立叶有限差分（,FFD）,法，综合了付立叶方法与有限差分法的优点，对陡倾角地层和强横向变速介质有很好的适应性，是目前精度最高的叠前深度偏移方法。下面主要在叠前深度偏移方法的基础上讨论,适用于各向异性介质,的保幅偏移，以及基于,粘滞声波、弹性波和粘弹性,波动方程的保幅偏移。,波 动 方 程 保 幅 偏 移,一、基于散射理论的保幅处理,1,方法原理,基于散射理论（波动理论与,Born,反演）的三维叠前保幅偏移，该方法是将解波动方程的问题转化为解格林函数的问题，其中格林函数的微小扰动可表示为,：,反演公式如下：,波 动 方 程 保 幅 偏 移,，J,为雅可比矩阵,2,实际应用,三维保幅偏移剖面,2.5,维保幅偏移剖面,波 动 方 程 保 幅 偏 移,偏移立体图的各种剖面显示,波 动 方 程 保 幅 偏 移,二、各向异性介质的保幅叠前深度偏移,1,方法原理（,VTI）,相移算子,反相移算子,波 动 方 程 保 幅 偏 移,2,模型试算（脉冲响应）,弱各向异性三维偏移脉冲响应,强各向异性三维偏移脉冲响应,波 动 方 程 保 幅 偏 移,3,实例,实例,1,：,（Oklahoma,地区,）,由各向同性,PreSDM,得到的偏移剖面,最终的各向异性偏移结果，井与目的层有了较好的匹配,波 动 方 程 保 幅 偏 移,实例,2,Husky,数据的各向同性,(,a),及各向异性,(,b),叠前深度偏移剖面对比图,Nordegg,/,Chungo,数据的各向同性,(,a),及各向异性,(,b),叠前深度偏移剖面,波 动 方 程 保 幅 偏 移,实例3,Blackstone,二维测线各种偏移剖面对比图,(,a),叠后时间偏移,(,b),叠前时间偏移,(,c),叠前深度偏移,(,d),各向异性深度偏移剖面,波 动 方 程 保 幅 偏 移,三、基于弹性波波动方程的保幅叠前深度偏移,1,方法原理,该方法从全弹性波波动方程出发计算出合成数据的垂直分量和水平分量，并从这些分量中提取出纯,P,波和纯,SV,波，然后分别用所提取的纯,P,波和纯,SV,波进行波动方程逆时偏移，这样可以得到两个深度偏移剖面。从理论上讲，这两个深度偏移剖面应该相等，因此可以相互比较，并用来检验偏移方法和偏移结果的实用性和可行性,波 动 方 程 保 幅 偏 移,2,模型试算,波 动 方 程 保 幅 偏 移,波 动 方 程 保 幅 偏 移,垂直分量,水平分量,波 动 方 程 保 幅 偏 移,四、基于粘滞声波波动方程的保幅叠前深度偏移,地下岩芯的实验室测试和油田的勘探开发实践，都已经证实地下介质或多或少地表现出非完全弹性，即粘弹性行为,(,衰减和频散,),。,1,基本原理,标准线性体粘滞声波方程表示为：,M,为复模量，为密度,复速度,波 动 方 程 保 幅 偏 移,定义复波数,则上、下行波场的相移算子分别表示为：,采用线性最优与非线性最优相结合的方法设计最佳的褶积延拓算子，该算子是频率、速度及粘滞性品质因子,Q,的函数，算子结合,McClellan,变换（,Hale,1991,）,可以实现三维保幅叠前深度偏移，由于其计算量较大的缘故，可以通过并行算法来实现。,波 动 方 程 保 幅 偏 移,2,模型试算,时间域地下构造模型,时间域地下构造模型的粘滞声波自激自收响应,波 动 方 程 保 幅 偏 移,对粘滞声波自激自收响应做声波偏移后的结果,对粘滞声波自激自收响应做粘滞声波偏移后的结果,波 动 方 程 保 幅 偏 移,中央凹陷地质模型,声波叠前深度偏移的结果,CDP80,处粘滞性声波单炮记录,粘滞声波叠前深度偏移的结果,