资源描述
Slide Title,储层地球物理实验室,Reservoir Geophysics Laboratory,2011,年,3,月,Slide Title,2011,年,3,月,地震相干体技术概述,4,主要内容,1,相干体的概念,2,相干体技术的历史,3,相干体算法介绍,相干体技术的应用注意事项,5,相干体技术在油气勘探中的应用,相干体的概念,相干体技术,是,利用地震信号相干值的变化来描述地层、岩性等的横向非均匀性,进而研究断层、微断裂的空间分布,地质构造异常及岩性的整体空间展布特征。,相干体,是指由,三维,地震数据体经过相干处理而得到的一个新的数据体,其基本原理是在三维数据体中,求每一道每一样点处小时窗内分析点所在道与相邻道波形的相似性,形成一个表征相干性的三维数据体,即计算时窗内的数据相干性,把这一结果赋予时窗中心样点。,(,Amoco,公司,),相干体,地震数据,参数测试,:,相干体算法,时间孔径,倾角扫描间隔,空间孔径,用相干体算法求每一道,每一样点处的相干值,相干体的概念,由三维地震数据得到相干切片的方法,三维地震数据,相干体,沿解释层位提取的地震数据子体,层拉平,相干体计算,相干水平切片,相干时间切片,提取给定时间的相干值,沿解释层位提取相干值,相干体计算,4,主要内容,1,相干体的概念,2,相干体技术的历史,3,相干体算法介绍,相干体技术的应用注意事项,5,相干体技术在油气勘探中的应用,相干体技术的历史,相干体技术原来是,Amoco,公司的专利,1995年由,Amoco,公司的,Mike Bahorich,引进石油勘探业。,相干体分析技术是当时地球物理界最具突破性的奇思妙想,,,Amoco,公司把这些算法简称为“,C1,”,,,“,C,2,”和“,C,3,”。,该技术揭示了波场的空间变化情况,直接从3,D,地震数据体中定量地得到断层和地层特征,不受任何解释误差的影响,极大地提高了解释精度,并能得到很多通常被忽略的重要信息,因而很快得到了广泛认可。,1996,,相干技术公司(,CTC),成功地将相干体技术商业化并拥有该技术唯一的许可证;,1999,年,,CoreLab,公司收购,CTC,公司;,2000 Core Lab,从,BP Amoco,公司购得相干体技术全套专利。,此后,相干体技术在地震油气勘探领域得到了广泛 的应用和发展,近年来又发展了,基于,几何结构张量的相干体技术等,,在与可视化结合,、与小波变换结合、,相干体的自动化解释研究等方面,都有较大发展,并向叠前地震数据应用的方向延伸。,4,主要内容,1,相干体的概念,2,相干体技术的历史,3,相干体算法介绍,相干体技术的应用注意事项,5,相干体技术在油气勘探中的应用,相干体算法介绍,第一代相干体技术:基于互相关的相干体技术,(Correlation),第二代相干体技术:基于相似的相干体技术,(Semblance),(,1,)基于多道相似的相干体技术,(,2,)基于曼哈顿距离的波形相似性算法,(,3,)基于复地震道的相干体技术,第三代相干体技术:基于本征结构的相干体技术,(Eigenstructure),新一代相干体技术,(,1,)基于几何结构张量的相干体技术,(,2,)基于高阶统计量的相干体技术,(,3,)基于小波变换的多尺度相干体分析技术,(,4,)基于曲波变换的相干体分析技术,第一代相干体技术:,基于互相关的相干体技术,相关算法是根据随机过程的互相关分析,计算相邻地震道的互相关函数来反映同相轴的不连续性。这种算法只能有三道参与计算。,在纵测线方向,t,时刻计算归一化互相关:,在横测线方向,t,时刻计算归一化互相关:,地震道的空间组合模式,第一代相干体技术:,基于互相关的相干体技术,相关运算法计算速度相对较快,但由于参与计算的地震道数少,对于有相干噪声的资料,仅用两道数据确定视倾角会有很大误差。再者,三点互相关算法假设地震道是零平均信号,当相关时窗长度超过地震子波长度时,这种假设才基本成立,即要求窗口大于地震反射的最长周期,显然,这样降低了计算得到的相干体数据的垂向分辨率。,振幅数据切片,C1,算法相干体切片,式中,,C,ii,(,i,=1,,,2,)为第,i,道的自相关量;,C,ij,(,i,=1,,,2,)为第,i,道和第,j,道的互相关量。视倾角(,p,q,)中,p,和,q,分别为,x,方向和,y,方向上的地震道之间的时移量。,沿视倾角(,p,q,)的相干值,C1,为:,定义一个以分析点为中心的包含,J,道的椭圆或矩形分析时窗,,,则,相干值,C2,为分析时窗内平均道的能量与所有道的能量,比,即,:,式中:下标,j,表示落在分析时窗内的第,j,道;,x,j,和,y,j,表示第,j,道与分析时窗内中心点,t,在,x,和,y,方向的距离,;,p,和,q,分别表示分析时窗内中心点,,,所在局部反射界面,x,和,y,方向的视倾角,。,以分析点为中心的椭圆分析窗口,以分析点为中心的矩形分析窗口,第二代相干体技术,(,1,)基于多道相似的相干体技术,振幅数据切片,C2,算法相干体切片,第二代相干体技术,(,1,)基于多道相似的相干体技术,应用基于多道相似的算法可以对任意多道地震数据计算相干性,对地震资料的质量限制不是很严,能够较精确地计算有噪声数据的相干性、倾角和方位角。用一个适当大小的分析窗口,能够较好地解决提高分辨率和提高信噪比之间的矛盾,而用相关算法则不能解决这个问题。因此,该算法,具有较好的适用性和分辨率,而且具有相当快的计算速度,。,对于地震数据体中的相干计算道,,设,G,为中心道,,,H,为相邻道,,定义归一化的,Manhattan 距离计算公式为:,式中,:,M,d,为归一化的Manhattan距离,,n,为相关时窗长度(样点数),,d,为地层倾斜时的时间延迟值,,N,为中心道,G,相关时窗中点的时间样点值。在式中当中心道,G,与相邻道,H,相同时,,M,d,等于零;当中心道,G,与相邻道,H,完全不相似即极性相反时,,M,d,等于1;在其它情况下,,M,d,的值介于0与1之间。,基于归一化的,Manhattan,距离的相干算法只涉及两个随机信号的加减运算,不涉及普通相干计算中的两个随机信号的乘积运算,因此,具有较高的运算效率,。,北偏东,45,度,3,点直线型道组合相干体切片,北偏西,45,度,3,点直线型道组合相干体切片,拐角型,3,道组合相干体切片,第二代相干体技术,(,2,)基于曼哈顿距离的波形相似性算法,如果在小的垂向时窗内应用,基于多道相似的相干体算法,,在计算零值附近同相轴的相干体时会造成一些假象,,,一旦地震信号的振幅低于背景噪音,,,相干算法将识别为地震数据不相干,,,将产生低相干的假象。应用复地震道技术就可以解决这个问题,,,对应的相干体计算公式为,:,式中:,上标,H,表示地震数据的希尔伯特变换或正交地震道。,当地震道的数据为零时,,,其对应的正交地震道数据的量值最大,;,当正交地震道的数据为零时,,,其实际地震数据的量值最大,,,这样相干体的算法就比较稳定,,,保真性好。,振幅数据切片,相干体切片,第二代相干体技术,(,3,)基于复地震道的相干体技术,该矩阵对应的协方差矩阵为:,第三代相干体是通过计算地震数据体的本征值获得的。在算法分析中,首先从给定的分析时窗内提取多道地震数据生成样点矢量,,,由这些样点矢量构成矩阵:,该协方差矩阵是一个对称的、半正定矩阵,其所有的本征值大于或等于,0,。计算协方差矩阵的本征值和本征向量,那么基于本征结构相干性估计可定义为:,第三代相干体技术:,基于本征结构的相干体技术,基于本征结构分析的相干体,算法,应用了主元素分析的思想,在主元素分析中,第一主元素总是通过协方差矩阵本征向量构成的立体角与发散椭球面长轴吻合。椭球面的主轴是由矩阵的本征向量确定的,本征值等于椭球面半轴长度。最大本征值反映了原始观测点信息的公共部分,即相干性。,相同的两道及通过协方差矩阵本征向量构成的椭圆,trace1 trace2,trace1,trace2,minor axis,major axis,“1”,“0”,不同的两道及通过协方差矩阵本征向量构成的椭圆,trace1 trace2,trace1,trace2,minor axis,major axis,“1”,“0”,第三代相干体技术:,基于本征结构的相干体技术,(,a,)振幅数据切片;(,b,),C1,算法切片;(,c,),C2,算法切片;(,c,),C3,算法切片;,几种相干算法比较,第三代相干体技术:,基于本征结构的相干体技术,本征算法是通过多道本征分解处理来计算波形相似性的一种方法,,虽然该,算法计算速度较低,但它,具有比相似系数算法更高的分辨率,。,接着,,使用方向导数构建梯度结构张量,:,Randen等提出用几何结构张量方法(Geometric Structural Tensor)进行地震三维结构属性的研究,这种,几何结构张量包含了反射界面的倾角和方位角信息,可以稳健地估算时窗内分析点的反射界面的倾角和方位角,。,首先,计算三维数据体每一点梯度矢量:,第三步,平滑张量矩阵,T,的元素:,第四步,计算平滑后的张量矩阵的特征值及相干值,:,新一代相干体技术,(,1,)基于几何结构张量的相干体技术,断层的自动识别(,a,)原始地震剖面及(,b,),GST,处理后的相干剖面,GST,相干切片解释(,a,)振幅数据切片及(,b,),GST,相干切片,新一代相干体技术,(,1,)基于几何结构张量的相干体技术,高阶统计量,具有许多优点,,从,提高相干体算法抑制噪声的能力和减小计算量,的角度出发,,陆文凯等人,将高阶统计量方法与相干体技术相结合,,,以 C1 算法为基础,,,提出了一种新的基于高阶统计量的相干体(CHOS)算法,,该算法仅需三道地震记录同时参与运算:,式中:是地震道 与 之间的时间延迟,是地震道,与 之间的时间延迟。最后得到的某一样点处的相干值为:,新一代相干体技术,(,2,)基于高阶统计量的相干体技术,振幅数据切片,CHOS,算法相干体切片,C2,算法相干体切片,C3,算法相干体切片,新一代相干体技术,(,2,)基于高阶统计量的相干体技术,多分辨率分析又叫多尺度分析,其中小波变换是其主要的分析手段。同傅氏变换相比较,小波变换具有更好的局部化特性,可以任意调节空间的分辨率,即有“变焦”功能,被誉为数字上的显微镜。,利用,小波变换可以很方便地进行多分辨率分析,从而进行精细的油藏描述和解释。将小波变换多尺度思想应用到相干体中,使小波多分辨率相干体在断层和裂缝的分频解释中见到更加明显的效果,。,新一代相干体技术,(,3,)基于小波变换的多尺度相干体分析技术,(a)(b)(c)(d)(e)(f),(a)原始地震剖面段,(b)-(f),不同频,段的小波分解,王西文等将小波分析引人到相干计算中,首先利用小波域分频方法计算地震数据各个频带内的瞬时特征参数,然后用互相关算法计算各个频带内的地震相干数据体,最后通过重构系数,对一定频带内的相干体放大或缩小主要突出特定频段的相干体,分频重构的相干体易于突出被忽略的小断层信息。,新一代相干体技术,(,3,)基于小波变换的多尺度相干体分析技术,(,a,)低频相干切片 (,b,)中频相干切片 (,c,)高频相干切片,将小波变换多尺度思想应用到相干体中,使小波多分辨率相干体在断层和裂缝的分频解释中见到更加明显的效果,。但小波变换方向性提取不足,只能在水平、垂直、对角线等几个有限的方向进行表示和检测,不能有效检测和表示信号中的线状变化特征。,Curvelet,变换结合了,Ridgelet,变换的各向异性,(anisotropy),特点和小波变换的多尺度,(multiscale),特点,具有很强的方向性,能为信号处理提供更多的信息。将曲波变换与相干算法相结合,可得到裂缝发育带的强度及走向分布图。,裂缝描述,岩石尺度,微米至厘米级,地质尺度,米级至千米级,多尺度性,矿物尺度,纳米至微米级,地层尺度 厘米至,米,级,裂缝张
展开阅读全文