地震属性分析技术课件_第1页
地震属性分析技术课件_第2页
地震属性分析技术课件_第3页
地震属性分析技术课件_第4页
地震属性分析技术课件_第5页
已阅读5页,还剩198页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

地震属性分析技术.ppt地震属性分析技术.ppt课程内容第一讲地震属性概述

第二讲地震属性提取与分析第三讲地震属性优化处理第四讲地震属性应用实例第五讲储层预测结果评价与检验课程内容第一讲地震属性概述一、地震属性的发展历程二、地震属性的定义三、地震属性的分类四、地震属性的发展趋势第一讲地震属性概述地震属性分析技术一、地震属性的发展历程第一讲地震属性概述地震属性分析技术“地震属性(Seismicattribute)”一词于20世纪70年代开始引入地球物理界。起初,国内在译名上并不完全统一,类似的译名还有地震特征、地震参数、地震标志等等,直到20世纪90年代初才基本统一称为地震属性。第一讲地震属性概述一、地震属性的发展历程“地震属性(Seismicattribute)”一词于20地震属性的发展大致经历了以下3个阶段。第1阶段:20世纪60年代末到70年代末,为起步阶段。该阶段以“亮点”技术为代表,这时的属性研究既没有考虑地震资料的运动学、动力学特征,也没有特定的地质含义,而仅仅是对地震剖面特征的一种定性描述与分析,如利用亮点标志解释气藏。第一讲地震属性概述地震属性的发展大致经历了以下3个阶段。第1阶段:20世纪第2阶段:20世纪70年代末到80年代末,为迅速发展阶段。该阶段以属性定量提取方法大量出现为主要特征。这是地震属性研究蓬勃发展的阶段,地震属性的应用开始走向各个领域,如储层预测、油气分析、物性研究等。这也是地震属性应用较为混乱的阶段,有多达几十种的地震信息被提取出来,但多数方法仅停留在地震波场的几何学、运动学、动力学等特征的研究上,没有对地震属性所代表的地质意义进行分析与解剖。第一讲地震属性概述第2阶段:20世纪70年代末到80年代末,为迅速发展阶段。该第3阶段:20世纪90年代以后,基本成熟阶段。这一阶段以多维属性的出现为主要标志。90年代初,以相干、倾角、方位角等为代表的一批多维属性开始出现。这类属性的一个显著特点是能够直观地反映地层的结构性信息,如倾角、方位角反映地层的视倾角/倾向,相干异常指示断裂带、裂缝(孔洞)发育带等。因此,它们的开发与应用,很快得到有关各方的广泛认同。与此同时,属性标定与优化的方法也开始大量涌现,地震属性研究开始向规范化、科学化方向健康发展。第一讲地震属性概述第3阶段:20世纪90年代以后,基本成熟阶段。这一阶段以多维第一讲地震属性概述二、地震属性的定义对地震属性的定义主要有以下几种:1、WesternAltlasInternational公司的QuincyChen与SteveSidney所给出的定义:“地震属性是地震资料的几何学、运动学、动力学及统计学特征的一种量度”。这一定义基本属于纯数学定义。第一讲地震属性概述二、地震属性的定义对地震属性的定义主要有2、LandmarkGraphics公司的ArthurE.Barnes对地震属性的定义是:“地震属性是一种描述和量化地震资料的特性,是原始地震资料中所包含全部信息的子集;地震属性的求取是对地震数据进行分解,每一个地震属性都是地震数据的一个子集”。该定义强调了地震属性的提取过程,但未对地震属性在地学中的作用进行阐述。定义中涉及的地震数据分解是一种广义的波场分解,没有统一的计算准则,包罗了所有的地震属性计算方法。第一讲地震属性概述2、LandmarkGraphics公司的ArthurE3、从应用地球物理学的角度看,比较合理的定义是:地震属性是地震数据中反映不同地质特征(信息)的分量或子集,是刻画、描述地层结构、岩性以及物性等地质信息的地震特征量。第一讲地震属性概述在众多的地震属性中,有些对特定的油藏环境比较敏感,有些对不易检测的地下界面异常更有利,还有些直接用于烃类检测或油气预测。

3、从应用地球物理学的角度看,比较合理的定义是:地震属性是地目前,大多数学者按照数学、物理学方法对地震属性进行分类(6家)。1、刘企英将地震属性分为振幅、频率、速度、吸收衰减、波形与时间6大类。2、AlistairR.Brown将地震属性分为时间、振幅、频率及衰减4大类。第一讲地震属性概述三、地震属性的分类目前,大多数学者按照数学、物理学方法对地震属性进行分类(第一讲地震属性概述3、20世纪90年代初Taner等人将地震属性分为几何属性和物理属性。几何属性通常与波形及地震层位的几何形态(如倾角,方位,曲率等)有关;物理属性包括运动学和动力学属性,主要有速度、振幅、频率、衰减等。4、1996年Brown等人将地震属性分为叠后属性与叠前属性。这种分类方法反映出人们对地震信号保真度和地震解释精度的日益重视。第一讲地震属性概述3、20世纪90年代初Taner等第一讲地震属性概述5、到20世纪90年代末,QuincyChen与SteveSidney在上述分类方法的基础上,提出了一种比较完整的分类方法,他们将叠后属性与叠前属性看成属性技术发展的2个阶段,就象叠后偏移成像和叠前偏移成像一样。在这个意义上,将地震属性分为几何学属性、运动学属性、动力学属性和统计学属性4大类。统计学属性大多是采用统计学方法获得的次生属性或导出属性,如相似性、广义主分量(GPC)、边缘平滑(EPS)等。一般没有明确的物理意义,但比其它属性含有更丰富的地质意义。第一讲地震属性概述5、到20世纪90年代末,Quinc6、从地震属性应用的角度出发,这些分类并不能反映地震属性在地球物理学或更广泛的地学中所起的作用。按地震属性在油气勘探开发的应用领域进行分类,较容易为属性的应用者所理解,也更具有地球物理学的特征。根据这一思路,曹辉等将地震属性分为下列2大类:第一讲地震属性概述6、从地震属性应用的角度出发,这些分类并不能反映地震属性第一讲地震属性概述(1)多尺度结构属性:包括相干、倾角、方位角、平行度、连续性、发散度、波度、随机度、曲率和间距等。其中大尺度结构属性用于构造分析;小尺度结构属性用于裂隙(缝)、孔洞分析和地震(沉积)相分析。这类属性的解释大都使用定性分析方法,根据专家知识将属性直接与地层结构特征相联系。第一讲地震属性概述(1)多尺度结构属性:包括相干、倾角、方(2)岩性或物性指示属性:这类属性最多,按动力学、统计学和反演等方法提取的属性以及部分按运动学提取的属性都可归入此类,主要有振幅、波阻抗、AVO、频率等及其它们的衍生参数。这类属性多数使用定量或半定量分析方法,例如通过井孔等资料对属性进行标定,建立定量或半定量模式,再以此模式识别岩性或推算物性参数。第一讲地震属性概述(2)岩性或物性指示属性:这类属性最多,按动力学、统计学综合上述几种分类,我们倾向于如下分类:1、建立在运动学、动力学基础上的地震属性类型,包括振幅、波形、频率、衰减特性、相位、相关分析、能量、比率等。2、以油藏特征为基础的地震属性类型,包括表征亮点、暗点、AVO特性、不整合圈闭或断块隆起异常、含油气异常、薄层油藏、地层间断、构造不连续、岩性尖灭、特殊岩性体等的地震属性。第一讲地震属性概述综合上述几种分类,我们倾向于如下分类:第一讲地震属性概述地震属性分析技术ppt课件地震属性分析技术ppt课件第一讲地震属性概述3、不同数据对象的地震属性类型,包括①以剖面为基础的属性,如传统的瞬时类属性,或经速度、声阻抗等特殊处理后的剖面;②以同相轴为基础的属性,提供了在地质分界面上或分界面之间的地震属性的变化信息,如沿层或层间瞬时属性、单道时窗的沿层或层间属性、多道时窗的沿层或层间属性。③以数据体为基础的属性,由三维地震数据体得到的相关类型的属性体具有很大的研究价值,例如可提供地震信号相似性和连续性方面的最佳信息。第一讲地震属性概述3、不同数据对象的地震属性类型,包括第一讲地震属性概述四、地震属性的发展趋势从20世纪60年代的烃类检测和亮点技术,到70年代基于复数道分析的瞬时属性,80年代的多属性分析(由叠前或叠后数据提取的各类属性),90年代的多维属性(倾角,方位和相干等)分析,直至21世纪的自动地震相分析等,地震属性技术的发展循序渐进,目前已逐渐走向成熟。对地震属性近几年的新进展,分以下8个方面加以简要介绍:第一讲地震属性概述四、地震属性的发展趋势从20世纪601、属性分类呈现几何学、物理学(运动学和动力学)、统计学并存,准属性(metaattribute)倍受关注。

准属性的含义就是对不同类型的地震信息进行目标拟合,生成一些按用户意图定义的属性,如地震相和照明技术等。这类属性虽然不能归入上述几种属性类型,但其地质意义明确,易于直观显示,能获得较好的解释效果。准属性实际上是多属性分析发展的结果,它属于多维属性的范畴。准属性分析通过智能化手段实现多属性的目标拟合,不仅能避免多属性解释的混乱,而且降低了对计算机性能的依赖性,使解释更为准确、便捷。第一讲地震属性概述1、属性分类呈现几何学、物理学(运动学和动力学)、统计学第一讲地震属性概述

地震相通常是由多种地震属性加波形分析确定的,它有助于精细地划分沉积单元,并揭示储层内部的细节。

照明技术则是通过对地震数据体中所有数据点进行同步照明(垂直于照明方向的特征被加强,平行于照明方向的特征被抑制)而形成的属性,包含反射系数、局部倾角和方位等信息,能很好地揭示各种构造细节,并弥补其他属性的不足。第一讲地震属性概述地震相通常是由多种地震属性加波形分析2、属性提取从线性(根据时窗数据提取的各种属性,如振幅、频率、相关函数等)向非线性(如小波变换的多尺度属性、应用照明技术生成的属性等)方向发展。3、属性分析智能化程度提高,多属性联合分析流行。4、属性反演(如波阻抗、弹性阻抗、AVO等)对岩石物理测试和正演模拟的依赖性增加,且以多维属性为主。第一讲地震属性概述2、属性提取从线性(根据时窗数据提取的各种属性,如振幅、第一讲地震属性概述5、属性研究从叠后向叠前推进,从纵波向弹性波扩展。6、属性应用的重点进一步转向开发,尤其是向精细储层描述、动态监测等方面转移。7、属性成果的可视化程度更高,地质意义更明了,使专家获取更丰富的地下信息。8、属性技术的定量化进程加快,解释的置信度提高。第一讲地震属性概述5、属性研究从叠后向叠前推进,从纵波第一讲地震属性概述总之,地震属性分析技术的发展已进入三维体属性分析时代,迫切需要立体显示、动态显示和仿真显示等现代化显示手段。恰好,三维可视化技术和虚拟现实技术的迅速发展基本满足了这一需求。目前,许多地震资料解释系统大都配上了现代化的可视化软件,虚拟现实技术也在一些大型石油公司得到应用。三维可视化和虚拟现实技术使地震属性分析结果更为直观、清晰,更易于反映油藏内部的地质意义,有利于启发解释人员的想象力,有利于从不同角度开展不同形式的分析与解释。第一讲地震属性概述总之,地震属性分析技术的发展已进入三以三维体可视化方式检测强振幅区域的工作步骤:①扫描并识别地质目标;②应用遮光度分析;③分离要确定的目标体;④振幅外延与构造凸现;⑤内部振幅变化;⑥沉积标志。第一讲地震属性概述以三维体可视化方式检测强振幅区域的工作步骤:①扫描并识别地第一讲地震属性概述在地震属性的显示方面,国外许多公司研发了新颖的显示技术,如Landmark公司的研究人员利用光照原理,形象地展现地震属性分析结果。休斯敦MagicEarth公司的研究人员提出了地质体可视化显示技术,该技术源于20世纪80年代的层位自动追踪,90年代的象素追踪,而属性可视化技术涉及到用三维象素追踪实现多属性和多样本的有机结合。地质体可视化不仅能使地学家检验各种数据集的多种解释结果,而且能提高解释速度和精度。第一讲地震属性概述在地震属性的显示方面,国外许多公司研第一讲地震属性概述地震属性分析结果的另一个发展趋势是地质意义更为明了。一方面,地震地层学和层序地层学的广泛应用丰富了地震属性技术的地质内涵,它不仅给出了地震属性的总体科学构思,而且激励人们利用地震属性建立地震相,再将地震相转化为地质相(层序体系、沉积相等);另一方面,多维属性分析使人们淡化了属性的具体类别,而更加注重围绕地质解释目标的有机结合,使得形成的属性体本身就具有直观的地质意义,而3D可视化和虚拟现实技术则使这种地质意义更加凸现。因此,未来地震属性总体将更加地质化而不是地球物理化。第一讲地震属性概述地震属性分析结果的另一个发展趋势是地第一讲地震属性概述目前,地震属性分析技术在正演模拟、多分量研究、多属性联合分析与反演、强化地质含义和改善显示手段等方面取得了明显的进展,但在叠前深度域、弹性波分析、属性解释的定量化及风险评估方面的进步不大。未来高难度油气勘探开发对地球物理技术的需求中除精确成像外,还有待在下列7方面进行不懈的研发:(1)用定量属性分析、多分量反演和非常重要的岩石物理及井资料标定,使地下储层中物性的空间展布和流体特性随时间的变化得到可靠的估计。第一讲地震属性概述目前,地震属性分析技术在正演模拟、多第一讲地震属性概述(2)不确定性量化和风险预测将与地球物理成果输出联成一体,被用于整体风险评估。(3)必须注重大数据量和多种类型数据体的实际交互,促进各种类型专业软件的无缝集成。(4)开展动态储集性能(如相对孔渗饱参数)的地震预测研究。第一讲地震属性概述(2)不确定性量化和风险预测将与地球物理第一讲地震属性概述(5)加强单点万道、弹性(矢量)地震处理方法以及解释分析方面的技术开发。(6)加速解决与相干噪声及多次波、各向异性、衰减、空间或时间分辨率以及保真度等有关的地震波传播、记录和处理问题。(7)加强时移地震属性分析和高分辨率拟地质统计反演结果与井孔资料综合,建立油藏精细模型方面的研究。第一讲地震属性概述(5)加强单点万道、弹性(矢量)地震处理课程内容第一讲地震属性概述

第二讲地震属性提取与分析第三讲地震属性优化处理第四讲地震属性应用实例第五讲储层预测结果评价与检验课程内容第一讲地震属性概述第二讲地震属性提取与分析一、提取地震属性的分析方法二、地震属性的描述与适用性三、相干属性体及其应用四、地震属性的物理意义五、地震属性间的相关性地震属性分析技术第二讲地震属性提取与分析一、提取地震属性的分析方法地震属性一、提取地震属性的分析方法1、自相关分析2、傅立叶谱、功率谱分析3、振幅特征分析4、复地震道分析5、地震记录的信噪比分析6、地震记录的相对分辨率分析7、地震记录的自回归分析8、地层平均吸收衰减特性分析9、相干体分析10、AVO分析11、各种变换和反演方法等第二讲地震属性提取与分析一、提取地震属性的分析方法第二讲地震属性提取与分析1、自相关分析设f(n)为离散地震信号,N1和N2为分析时窗的起始、终止时间所对应的离散时间序号,时窗内离散时间样点数为N=N2-N1+1,f(n)自相关函数ACF(τ)表示为:自相关函数如右图所示。第二讲地震属性提取与分析1、自相关分析第二讲地震属性提取与分析自相关函数是地震记录特征的反映,是地震记录中重复性的标志。在一般情况下,自相关函数主极值幅度代表着记录段的能量;主极值宽度与记录的视周期有关,频率低的信号,主极值宽度大;频率高的信号,主极值宽度窄;旁极值的幅值和面积表示地震记录的重复性及延续时间的长短。当反射层具有薄互层结构,反射记录出现干涉现象时,自相关函数幅值和面积增大。第二讲地震属性提取与分析自相关函数是地震记录特征的反映,是地震记录中重复性的标志。在2、傅立叶谱、功率谱分析谱分析是在频率域内描述地震记录特征的重要方法,它有两种形式,一是傅立叶谱分析,用于确定性函数;二是功率谱分析,用于随机过程。当用于分析的地震数据是一个均值为零的随机过程,功率谱为它的一个统计特性,可以较好表示反射波的特征;当用于分析的地震数据是一个确定的时间函数,或记录的信噪比较高、分析时窗中有稳定的反射波出现,使用傅立叶谱分析描述反射波特征较为适宜。第二讲地震属性提取与分析2、傅立叶谱、功率谱分析第二讲地震属性提取与分析地震波的频谱与其波形的关系是互为正、反傅立叶变换的关系,其数学表达式为:第二讲地震属性提取与分析地震波的频谱与其波形的关系是互为正、反傅立叶变换的关系,其数获取功率谱P(f)的方法:1、利用振幅谱与功率谱的关系,即P(f)

C2A2(f)2、利用自相关函数的FFT得到P(f)3、伯格谱方法-利用自回归分析的相关系数定义伯格谱。第二讲地震属性提取与分析获取功率谱P(f)的方法:3、伯格谱方法-利用自回归分析的相3、振幅特征分析反射波振幅或能量特征是地震岩性解释和储层预测常用的动力学属性,利用振幅特征分析可提取下列一些属性:记录能量、最大振幅值、均方根振幅等,这类属性反映了目标层内波阻抗、地层厚度、岩石成分、孔隙度及含流体成分的变化,可用来识别振幅异常或用于层序特征分析;也可用来追踪地层学特征,如三角洲、河道、各种扇或特殊岩性体;还可用于识别岩性变化、不整合、气体以及流体的聚集等。第二讲地震属性提取与分析3、振幅特征分析第二讲地震属性提取与分析4、复地震道分析Taner等(1979)解释的复数道由实部(常规地震道f(t))和虚部(正交道h(t))组成,即式中虚数i是-1的平方根。正交道h(t)可以通过Hilbert转换的常规地震道独立地确定,即

式中*表示褶积。实际上,h(t)是f(t)的900相移形式,所有其他的复数道都可以从正交道和实数道中求得。第二讲地震属性提取与分析4、复地震道分析第二讲地震属性提取与分析瞬时频率:三瞬剖面属性定义为:瞬时振幅强度:瞬时相位:第二讲地震属性提取与分析瞬时频率:三瞬剖面属性定义为:瞬时振幅强度:瞬时相位:第二5、地震记录的信噪比分析在使用地震资料进行储层预测前,对资料的信噪比作定量评价是必要的。在一定的地震地质条件下,记录的干扰背景很强;在另外的地质条件下,记录的干扰背景可能较平静,如油藏部分经常有地震记录规律性变差、信噪比低这一特点。因此,测定地震记录的信噪比,也具有重要的地质解释意义。从技术上看,信噪比经常用来评价野外工作的施工质量,考核资料处理效果。建立信噪分离方法,定量评价地震记录的信噪比,具有很好的现实意义。第二讲地震属性提取与分析5、地震记录的信噪比分析第二讲地震属性提取与分析设地震记录f(t)由信号q(t)和噪音n(t)迭加而成,且噪音n(t)为随机的,与信号q(t)不相关,是一个满足正态概率分布的稳定的随机过程:f(t)=q(t)+n(t)在这种假设条件下,用多道记录可以分别获得信号和干扰的自相关函数的可靠估计。第二讲地震属性提取与分析设地震记录f(t)由信号q(t)和噪音n(t)迭加而成,且噪考虑到信号与干扰不相关,其间互相关函数为零,从而有:式中,k=1,2,…,K为多道处理时使用的各地震道道序号,K为参与处理的总道数,T为各道分析时窗长度。第二讲地震属性提取与分析考虑到信号与干扰不相关,其间互相关函数为零,从而有:6、地震记录的相对分辨率分析所谓分辨率,指的是分辨相距很近的分界面而来的反射波信号的能力。主要决定于地震记录上反射脉冲的延续时间、地层厚度以及反射波的频率成分等。由于不可能在地震记录上直接测定反射脉冲的延续时间,因此分辨率的量度按地震记录自相关参数来定义:ACF(

)

0S1S4S3S2

第二讲地震属性提取与分析6、地震记录的相对分辨率分析ACF()0S1S4S3S2其中θ为自相关函数主极值半周期宽度,S1为主极值半周期内自相关曲线所包含的面积,S2、S3、S4,简记S2-4为自相关曲线依次三个旁极值面积之和。θ与反射波视周期有关,视周期越长,分辨率越低,所以参数Rλ与θ成反比;自相关函数的面积表示地震波能量的分布情况。相位个数越少、脉冲越窄,则地震波分辨率越高,对这类脉冲主极值半周期面积S1越大;相反,脉冲延续时间长、相位个数多,则地震波分辨率低,而自相关函数旁极值面积S2-4越大;所以参数Rλ与S1成正比,与S2-4成比。第二讲地震属性提取与分析其中θ为自相关函数主极值半周期宽度,S1为主极值半周期内7、地震记录的自回归分析区别于上述地震记录反射特征分析,自回归分析是一个数学分析方法,使用自回归模型,或AR模型,将地震记录表示为本身样点的线性组合,用多项式的系数来表达地震记录的反射特征。自回归系数作为反射特征参数,它们不具有象其它参数那样明显的物理意义,而是一个抽象的数学系数,描述着地震记录的面貌和特征。第二讲地震属性提取与分析7、地震记录的自回归分析第二讲地震属性提取与分析作为随机过程xt的一段记录,可以用M阶自回归模型表示,即以xt过去时样点线性组合来表示该过程现在时的估计值:其中k=M+1,M+2,…,M0;也就是说,如果地震记录的前M个样点是已知的,则M+1以后的样点值可以用前面的样点值来预测。这是一个线性滤波问题,使用最小平方法解解线性方程组:第二讲地震属性提取与分析作为随机过程xt的一段记录,可以用M阶自回归模型表示,即可以求得滤波因子{aj,j=1,2,…,M}自回归分析的应用包括两方面内容,一方面是计算伯格功率谱,另一方面,自回归分析用于油藏描述,可以预测油气藏的存在、类型及油藏边界。第二讲地震属性提取与分析可以求得滤波因子{aj,j=1,2,…,M}第二讲地震属8、地层平均吸收衰减特性分析地震波传播的实际介质并非是理想的完全弹性介质,在许多情况下具有非弹性性质。研究弹性波吸收与岩相的关系,结果表明这种关系在不同频率上是不同的。对油气藏模型研究结果表明,其物理性质分布是复杂的、不均匀的。对弹性波吸收性质主要决定于岩石骨架的弹性性质,影响较大的因素还有孔隙度、孔隙可压缩性及孔隙饱和液成分,特别是孔隙中含有天然气成分对吸收性质的影响更为明显。第二讲地震属性提取与分析8、地层平均吸收衰减特性分析第二讲地震属性提取与分析对地震波传播的实际介质的吸收性质作一个归纳,以下6点事实是可以肯定的。(1)岩石吸收性质从地表到地球深处变化范围很大,且对纵波和横波各不相同。横波的吸收衰减高于纵波的吸收衰减。(2)各类波的吸收系数与频率有关,随频率增加而增大,接近于线性关系。(3)吸收系数与地震速度之间存在明显的对比关系,高速的岩石,吸收系数低;而低速的岩石,吸收系数高。吸收系数如同速度一样,频散异常现象较弱。第二讲地震属性提取与分析对地震波传播的实际介质的吸收性质作一个归纳,以下6点事实是可(4)地层岩石吸收性质首先决定于岩石保存状态和内部结构。矿物颗粒和粒度对吸收性质影响不大,地层静压力随深度而加大,使岩石压紧,结构致密,引起吸收系数降低。受到破坏的岩石结构,将使它的吸收系数提高。(5)由固、液、气构成的多相介质中,对吸收性质影响最显著的是气态物质。在岩石孔隙饱和液中渗入少量气态物质,可以明显提高对纵波能量的吸收。第二讲地震属性提取与分析(4)地层岩石吸收性质首先决定于岩石保存状态和内部结构(6)对大多数地区,泥岩中平均吸收系数比砂岩高。岩石的吸收系数与砂岩系数有关,随砂岩系数增加,吸收系数降低。岩石的吸收系数还与其埋藏深度有关,随埋藏深度增加而减少。地震勘探中用来表示地层吸收性质的参数有吸收系数α、衰减因子h、对数衰减率δ、品质因子Q等。第二讲地震属性提取与分析(6)对大多数地区,泥岩中平均吸收系数比砂岩高。岩石

吸收系数

用来表示地震波振幅A沿传播距离的衰减,即其中x为波的传播距离,A0为起始振幅,t为波的传播时间。衰减因子h:对数衰减率

表示地震波振幅在一个波长λ距离上的衰减量或与之相应的在一个周期T内的衰减量,即:第二讲地震属性提取与分析吸收系数用来表示地震波振幅A沿传播距离的衰减,即其中从δ=αλ=hT关系看,对数衰减率是一个无量纲参数。根据大量实验数据证明,在地震频带内,吸收系数α与频率f成正比。正是表达了这种函数关系,比例系数应为一常数。忽略地震速度的频散现象,它是一个与频率无关的常量,则对数衰减率δ亦将是与频率无关的常量,单一地反映着传播介质的吸收性质。第二讲地震属性提取与分析从δ=αλ=hT关系看,对数衰减率是一个无量纲参数。根据大量

品质因子Q是用来比较地震波能量在一个波长范围内的相对变化。若地震波起始振幅为A0,在一个波长λ范围内衰减到A,则能量相对变化为:定义品质因子:在以上四个吸收参数中,实际工作常用的是对数衰减率δ和品质因子Q。第二讲地震属性提取与分析品质因子Q是用来比较地震波能量在一个波长范围内的相对变化求取地层平均吸收衰减的具体方法有:(1)傅立叶谱算法——这是一种直接根据地震反射波记录的振幅谱求对数衰减率的方法,本方法适用于分析时窗中反射波单一、稳定,记录的信噪比较高的情况。(2)功率谱算法基本公式:t为两个时间段记录的中心时间之差;S0(f)为吸收层顶面反射波的振幅谱,Si(f)为该层底面反射波的振幅谱。第二讲地震属性提取与分析求取地层平均吸收衰减的具体方法有:t为两个时间段记录的(3)复赛谱算法——复赛谱分析是一种非线性滤波,用以分解两个序列相褶积的过程。地震反射波可视为地震子波与脉冲地震记录的褶积。脉冲地震记录具有白噪性质,变化快;而地震子波一般较为稳定,变化慢。当介质无吸收作用时,子波无变化。这些性质符合两个过程可分辨的条件。第二讲地震属性提取与分析(3)复赛谱算法——复赛谱分析是一种非线性滤波,用以分解设地震记录f(t)为子波w(t)与脉冲地震记录k(t)的褶积:f(t)=w(t)*k(t),对此做傅立叶变换,得S(ω)=W(ω)K(ω);取对数得:lnS(ω)=lnW(ω)+lnK(ω),经分解后再取傅立叶逆变换,获得时间域复谱,称为复赛谱。第二讲地震属性提取与分析设地震记录f(t)为子波w(t)与脉冲地震记录k(t)的褶积(4)特征值法——这是利用地震记录相邻时窗的自相关函数特征值、振幅特征值的变化来判断目的层吸收性质的相对变化。①极大振幅比和均方根振幅比——取两个时窗中的极大振幅或均方振幅,求它们的比值。比值大小与界面性质、地层的吸收性质甚至含流体情况有关。第二讲地震属性提取与分析(4)特征值法——这是利用地震记录相邻时窗的自相关函数特②测定记录段自相关函数的主极值幅度变化——由于自相关函数主极值ACF(0)相应于振幅的平方和,则由主极值幅度变化可以确定介质对地震信号能量的吸收参数。③自相关函数的主瓣面积比或旁瓣面积比④低频能量之比⑤波峰数之比⑥主频之比第二讲地震属性提取与分析②测定记录段自相关函数的主极值幅度变化——由于自相关9、相干体分析相干体技术由相干技术公司(CTC)和Amoco公司发明,1997年获美国专利。该技术被称为是近几十年来三维地震解释方面最重要的突破。与原来揭示地下异常体的方法相比,相干体技术更能清楚地识别断层和地层特征。相干体技术的特有算法是通过三维数据体来比较局部地震波形的相似性。相干值较低的点与地质不连续性如断层和地层、特殊岩性体边界密切相关。对相干数据体作水平切片图,可揭示断层、岩性体边缘、不整合等地质现象,为油藏描述提供了识别油藏特征的有利证据。第二讲地震属性提取与分析9、相干体分析第二讲地震属性提取与分析10、AVO分析

AVO技术是利用CMP道集资料,分析反射波振幅随偏移距(也即入射角)的变化规律,估算界面两侧的弹性参数泊松比,进一步推断地层的岩性和含油气性。AVO技术的理论基础是Zoeppritz方程及其各种形式的简化关系。第二讲地震属性提取与分析10、AVO分析第二讲地震属性提取与分析Zoeppritz方程的简化公式很多,目前使用较为广泛的是Shuey简化式:

式中包含三部分:(1)当垂直入射时的反射振幅:

(2)当广角即

>临界角入射时;R(

)主要与速度变化有关,即:

第二讲地震属性提取与分析Zoeppritz方程的简化公式很多,目前使用较为广泛的是S(3)当中等入射角时,即0<

<300,近似有sin=tan,此时上式可写为:

或式中P为截距,反映垂直入射时的反射振幅;G为梯度,反映振幅随偏移距的变化率。利用P和G的线性组合表示AVO的属性W,即:

第二讲地震属性提取与分析(3)当中等入射角时,即0<<300,近似有sin=tan

a,b取不同的值,赋予W不同的物理意义(见下表)。比较AVO属性P、G和岩石背景趋势的差异有助于提高烃类识别能力。

abW意义10PRp0,垂直入射的P波反射系数-截距01G反射振幅随偏移距的变化率-梯度1/21/2(P+G)/2RP0-RS0,纵横波差异1/2-1/2(P-G)/2RS0,横波垂直入射的反射系数G0指示振幅随偏移距增加而增加AVO属性及物理意义

第二讲地震属性提取与分析a,b取不同的值,赋予W不同的物理意义(见下表)。比较上述分析与理论试算表明:在

<300情况下,入射角对反射系数的影响不是很大,这正是CMP叠加的理论基础,也是AVO的应用区间。下面通过应用实例来说明AVO属性在烃类检测或油气预测中的能力。

第二讲地震属性提取与分析上述分析与理论试算表明:在<300情况下,入射角对利用AVO资料进行烃类检测的实例

在进行AVO分析之前,必定要对AVO资料进行相应的处理。简要说明AVO资料处理中两个阶段应该注意的问题。一是AVO资料的预处理阶段。预处理阶段的主要任务是得到高质量的CMP道集记录,这是AVO处理成功的关键。结合野外数据及地下地质情况,AVO处理时应做好以下工作:①精细的波前扩散补偿;②震源组合与检波器组合效应的校正;③反Q滤波;④地表一致性处理;⑤叠前去噪处理;⑥叠前剩余振幅补偿;⑦精细的初至切除、速度分析和高精度的动校正;⑧叠前时间偏移。第二讲地震属性提取与分析利用AVO资料进行烃类检测的实例第二讲地震属性提二是AVO属性处理阶段。AVO属性处理是以CMP道集为输入,经角度道转换后,不同角度道组成一个角度道道集,从角道集上便可观测和分析地震反射振幅随入射角的变化。以角道集为输入,根据相应的公式即可拟合出垂直入射时的纵波反射、梯度及各种AVO属性剖面。AVO分析的基本方法包括:AVO的正演方法(如反射系数曲线分析法、非零偏移距地震模型法、流体替代模型法等)和AVO的反演方法。第二讲地震属性提取与分析二是AVO属性处理阶段。AVO属性处理是以CMP道集为输入,根据反演目标的不同可分为AVO属性参数反演、储层参数反演和叠前多波地震反演。例如,广义振幅属性反演可用于AVO交汇图分析;岩性参数变化量反演主要用于岩性和油气的定性预测;弹性参数反演主要用于定量储层预测;波型反演(如利用纵波AVO提取横波剖面或转换波剖面)常用于改善烃类检测能力。根据反演算法的不同,AVO反演又可分为迭代正演模拟法、线性反演和非线性反演。第二讲地震属性提取与分析根据反演目标的不同可分为AVO属性参数反演、储层参数图1声阻抗AI与弹性阻抗EI(300)反演曲线比较图1声阻抗AI与弹性阻抗EI(300)反演曲线比较图2(a)Sonic、AI反演图2(b)Sonic、EI曲线和入射角等于100时的反演曲线和入射角等于角道集叠加剖面300时的角道集叠加剖面图2(a)Sonic、AI反演图3墨西哥湾MC619-1井的典型岩石物性显示图3墨西哥湾MC619-1井的典型岩石物性显示AI与含油饱和度EI(300)与含油饱和度图4含油饱和度与AI、EI(300)间的关系第二讲地震属性提取与分析AI与含油饱和度EI(300)关于AVO烃类指示因子

常规的AVO分析方法(P,G两参数)是从地震资料中提取振幅,并通过截距P和斜率G两种属性将振幅随偏移距的变化与岩石物性联系起来。由于截距是纵波速度α和密度ρ的函数,斜率是横波速度β和纵波速度α的函数,所以三个主要弹性参数(α、β和ρ)彼此依赖,互有联系,不能分离。这在岩性变化剧烈、产层含气饱和度低的情况下就可能导致两参数AVO反演的失效。为此,不少研究者提出了一些改进方法,如Kabir建议利用密度差异作为含气饱和度的指示因子,Skidmore等采用三参数()的AVO岩性属性反演提高烃类检测能力。第二讲地震属性提取与分析关于AVO烃类指示因子第二讲地震属性提取与分析图5(a)产油区的振幅层位切片在振幅异常中,除了A5井外所有井都含气。A5井为含水井,B1井位于邻近含水区的薄气层,在振幅异常中的灰白区域则没有发现气藏图5(a)产油区的振幅层位切片在振幅异常中,除了A5图5(b)产油区的AVO层位切片在AVO异常中,除了A5井外所有井都含气。A5井为含水井,B1井位于邻近含水区的薄气层,在AVO异常中的灰白区域则没有发现气藏图5(b)产油区的AVO层位切片在AVO异常中,除了图5(c)产油区的密度异常层位切片在密度异常平面图上,A5井处在密度异常之外,B1井处在密度异常边缘上。图5(c)产油区的密度异常层位切片在密度异常平面图上图6(a)墨西哥湾某开发区的振幅异常分布图右图为墨西哥湾一个开发区的振幅异常分布图。研究区内4个AVO异常区块中的3个已经钻探并出气,第四区块还没有测试。图6(a)墨西哥湾某开发区的振幅异常分布图右图为墨西图6(b)墨西哥湾某开发区的AVO异常分布图在第一、第二断块中的3口井附近,AVO异常明显。在第三断块中的AVO异常已被钻井证实仅发现少量的气,且在开采几个月后被关闭。图6(b)墨西哥湾某开发区的AVO异常分布图在第一、第图6(c)墨西哥湾某开发区的密度异常分布图在第一至第三断块中的密度变化比较明显。第四振幅异常区还未钻探,它在密度异常图上具有与振幅和AVO相同的形状,表明具有良好的开采前景。图6(c)墨西哥湾某开发区的密度异常分布图在第一至第Rutherford和Williams的四类气层分类法:Ⅰ类为高阻抗含气砂岩(C1),其AVO特征为:零偏移距振幅强,且为正极性,AVO曲线呈减小趋势,当入射角足够大时可看到极性反转。Ⅱ类为近零阻抗差的含气砂岩(C2),其AVO特征为:零偏移距振幅很小,趋于零。故在零偏移距附近不易检测,随着偏移距的增大,其AVO特征变化较大。此类气层又可分为C2P与C2N。前者AVO曲线开始大于零,随着偏移距的增大振幅会减小,并出现极性反转;后者AVO曲线都小于零,随着偏移距的增大振幅绝对值在增加。Ⅲ类和Ⅳ类为低阻抗含气砂岩(C3和C4)。它比上覆介质的阻抗低,其AVO特征为:零偏移距振幅很强,呈负极性,AVO(指振幅绝对值)呈增加趋势的为第Ⅲ类;AVO呈减小趋势的为第Ⅳ类。Rutherford和Williams的四类气层分类法:Ⅰ类图7气层AVO分类模型图7气层AVO分类模型图8AVO的截距A和梯度B交汇图图8AVO的截距A和梯度B交汇图ⅠⅡⅡⅢⅣⅠⅡⅣⅢⅠⅡⅡⅢⅣⅠⅡⅣⅢ11、各种变换和反演方法小波变换是90年代比较活跃的地震属性提取方法,它不仅具有提高地震属性分辨率的潜力,而且能优化属性提取的时窗长度。与小波变换有关的地震属性包括:利用连续子波变换提取地震时频属性,时频分解能以低频成分提供较高的频率分辨率,以高频成分提供较高的时间分辨率;低频成分可用来识别储层中的油气,高频成分则有利于分辨薄层。利用离散小波变换提取振幅、瞬时频率和相位属性。小波变换为属性计算中采用多尺度分析和去噪技术提供了潜力。第二讲地震属性提取与分析11、各种变换和反演方法第二讲地震属性提取与分析二、地震属性的描述与适用性通过以上讨论我们知道:提取地震属性常用的分析方法包括复地震道分析、相关分析、富立叶谱分析、功率谱分析、自回归分析、数理统计分析等。经过相应的分析计算后,可得到一系列地震属性参数,包括:

1.瞬时属性(InstantaneousAttributes)-这是根据一定位置上地震数据的复数道分析技术导出的一系列地震属性参数,主要包括:第二讲地震属性提取与分析二、地震属性的描述与适用性第二讲地震属性提取与分析(1).瞬时真振幅f(t)

(2).瞬时积分振幅(正交道)h(t)

(3).瞬时相位

(t)=tg-1[h(t)/f(t)]

(4).瞬时相位的余弦cos(

(t))

(5).瞬时真振幅乘瞬时相位的余弦,f(t)·cos(

(t))

(6).瞬时频率ω(t)=d

(t)/dt

(7).振幅加权瞬时频率fa=(8).能量加权的瞬时频率fe=(9).瞬时频率的斜率,d[d

(t)/dt]/dt

(10).反射强度A(t)=

(11).分贝表示的反射强度,20·lgA(t)

第二讲地震属性提取与分析(1).瞬时真振幅f(t)第二讲地震属性提取与分析(12).反射强度的中值滤波能量(13).反射强度的斜率dA(t)/dt

(14).视极性-定义为实际地震道在反射强度波峰处的极性,其数值等于反射强度×sign;sign为反射系数的正负符号。常与反射强度一起用来检查沿层极性的横向变化。(15).响应相位-是从围绕反射强度叶瓣的瞬时相位中导出的;这是追踪地震子波的相位特性在时间、空间上改变的一种替代方法。可用于识别地震地层层序,检测岩相或流体含量的变化,当振幅特征相似时,可以作为区域性特征的识别标准。第二讲地震属性提取与分析(12).反射强度的中值滤波能量第二讲地震属性提取与分析(16).响应频率-是从围绕反射强度叶瓣的瞬时频率中导出的;这是追踪地震子波的优势频率在时间、空间上改变的一种替代方法。可用来识别与气聚集而引起的频率异常;与瞬时频率所起作用类似。

2.单道时窗属性(SingleTraceWindowedAttributes)-这是对单道、一定时窗的地震数据经振幅特征分析、相关分析、富立叶谱分析、功率谱分析、自回归分析、数理统计分析和计算后得到的一系列地震属性参数,包括:

第二讲地震属性提取与分析(16).响应频率-是从围绕反射强度叶瓣的瞬时频率中导出的;(1).平均振动能量(2).平均振动路径长度(3).最大波峰、波谷振幅(4).累积绝对振幅(5).复合绝对振幅(6).均方根振幅(7).复合包络差(8).平均过零点(9).带宽比(10).优势频率比(11).中心频率比(12).形心频率(功率谱面积的中心所对应的频率)比第二讲地震属性提取与分析(1).平均振动能量第二讲地震属性提取与分析(13).第一、第二、第三峰值谱频率(14).优势功率谱(15).优势功率谱的集中程度(16).有限的或特定的带宽能量(17).特定的与有限的能量比(18).衰减敏感带宽(有限带宽能量值除以优势频率值)(19).功率谱对称性(频谱相对中心频率的形状对称性)(20).功率谱斜率(21).相邻峰值比(22).自相关函数的峰值比(23).地震记录分辨率的评价参数

第二讲地震属性提取与分析(13).第一、第二、第三峰值谱频率第二讲地震属性提取与(24).样点的数学均值(25).振幅斜率(26).振幅峭度(27).正负振动比(28).样点值小于、大于门槛值的百分比(29).半能量时间(30).正负振幅个数比(31).复赛谱分析算法求取对数衰减率δc

(32).付立叶谱分析算法求取对数衰减率δf

(33).功率谱分析算法求取对数衰减率δp

(34).低频能量比LER=LE2/LE1

第二讲地震属性提取与分析(24).样点的数学均值第二讲地震属性提取与分析3.多道时窗属性(Multi-TraceWindowedAttributes)-这是对多道(地震道的空间组合模式如下图所示)、一定时窗的地震数据经相关分析、主元素分析或K-L变换等处理后得到的一系列地震属性参数,包括:(1).相关KLPC1,KLPC表示主元素分析法(2).相关KLPC2

(3).相关KLPC3(4).相关KLPC比,主分量间差值之比,即(PC1-PC3)/(PC1-PC2)(5).相关长度-定义为所有道互相关值降到0.5时用道数表示的平均距离。第二讲地震属性提取与分析3.多道时窗属性(Multi-TraceWindowed(6).平均相关(多道时窗内互相关的平均值)(7).加权相关(做互相关时乘以不同权系数)(8).最小相关值(多道时窗内互相关的最小值)(9).最大相关值(多道时窗内互相关的最大值)(10).相关峭度(多道时窗内平均相关的4次幂)

(11).相似系数第二讲地震属性提取与分析(6).平均相关(多道时窗内互相关的平均值)第二讲地震属上述属性中,最小相关强调局部地震道间的不相似性;最大相关用于检测大的地层边界;均值、中值、加权相关给出道间相似性的趋势,对噪音有压制作用,但有平滑效应;最大相关与最小相关之差、均值与中值之差、方差给出局部相似性之间的差异。(12).最大相关-最小相关(13).中值(14).均值减中值(15).方差第二讲地震属性提取与分析上述属性中,最小相关强调局部地震道间的不相似性;最大相

4.沿层构造属性(EventObjectStructureAttributes)-包括倾角分析,方位分析,边缘检测,差异检测等沿层构造属性。这些沿层构造属性可以揭示细小断层、地层学特征等,另外,通过这些属性的分布情况,可以较为直观的展示目的层的构造分布。(1).倾角分析-倾角定义为时间梯度的大小:(2).方位分析-定义为y方向时间梯度与x方向时间梯度之比的反正切:第二讲地震属性提取与分析4.沿层构造属性(EventObjectStruc

倾角/方位分析示意图

135-224315-4546-134225-314红黄绿蓝第二讲地震属性提取与分析

夏一段倾角方位分布图N夏一段倾角方位分布图N(3).倾角/方位分析(4).边缘检测从上面的讨论可知,从地震资料尤其是三维地震数据体中确实可以提取很多很多的地震属性参数,这些属性参数都是地下地层、岩性、物性特征的具体反映,例如有反映储层含油气特征的地震属性参数;有反映局部高振幅带的地震属性参数;有反映油藏厚度或断层特征变化的地震属性参数;也有反映储层频率吸收衰减的地震属性参数;还有研究储层裂缝及其发育带的地震属性参数,如反映储层吸收衰减特性的各种属性参数、沿层倾角方位检测的有关属性参数,如倾角分析,方位分析,边缘检测等沿层构造属性参数。第二讲地震属性提取与分析(3).倾角/方位分析第二讲地震属性提取与分析

1、相干体的概念计算地震相干数据体的目的主要是对地震数据进行求同存异,以突出那些不相干的数据。通过计算纵向和横向上局部的波形相似性,可以得到三维地震相关性的估计值。在出现断层、地层岩性突变、特殊地质体的小范围内,地震道之间的波形特征发生变化,进而导致局部的道与道之间相关性的突变。通过三维相关属性体的提取,就可以把三维反射振幅数据体转换成三维相似系数或相关值的数据体。三、相干属性体及其应用第二讲地震属性提取与分析1、相干体的概念三、相干属性体及其应用第二讲地震属性互相关分析算法:1995年由Bahorich,M.S.和Farmer,S.L.提出的,它是基于互相关的相干体算法;相似性分析算法:1998年由Marfurt,K.J.,KirlinR.L.和Farmer,S.L.提出的,主要利用地震道的相似性来得到相干属性;特征值分析算法:1999年由Gersztenkorn,A.和Marfurt,K.J.提出的,即通过计算协方差矩阵的特征值来得到相干属性。下表给出了三种相干体方法的基本特征比较。

2、相干体的主要算法互相关分析算法:1995年由Bahorich,M.S.和F算法名称算法描述贡献者计算公式互相关分析算法将纵、横测线的各道数据与中心道进行简单的互相关,得到相关系数。1995年,由Bahorich和Farmer引入到地震属性提取中,早前Taner和Koehler在速度分析中已有应用。相似性分析算法通过构建协方差矩阵来计算多道数据的相似系数,计算过程采用求和运算。1998年,由Marfurt、Kirlin等提出,其等价的计算公式可追踪到1984年Sheriff的求和公式。特征值分析算法计算协方差矩阵的特征值,将其中最大值与矩阵迹的比值作为相干值。1999年,由Gersztenkorn和Marfurt给出了详细的讨论,此前已有多人对这种基于特征值的方法进行了研究。相干体算法的基本特征比较算法名称算法描述贡献者计算公式互相关分析算法将纵、横测线的(1)InitialAlgorithmC1相干体分析技术的发展(2)SemblanceC2(3)EigenStructureC3C1是互相关分析法。下图展示了盐丘周围的放射状断层以及峡谷,用相干体的等时切片展示其形态十分有利。(1)InitialAlgorithmC1相干体分析C3标准的(左)与高分辨率本征结构(右)的相干体切片比较。C2大大改善了信噪比和稳定性;C3对断层和地层特征的成像可用相同参数序列且对其形态无任何损害。C3标准的(左)与高分辨率本征结构(右)的相干体切片比较。(a)为原始数据体切片(b)是互相关分析算法(c)是相似性分析算法(d)是特征值分析算法相干体切片(a)为原始数据体切片(b)是互相关分析算法(e)是局部结构熵(LSE)算法(f)是超级地震道-高阶统计学相干(ST-HOSC)算法多种相干算法的相干体切片(t=2000ms)比较相干体的算法简介(e)是局部结构熵(LSE)算法(f)是超级地震道-高阶普通三维切片相干体切片3、相干体的应用概述(1)展示断层发育细节普通三维切片相干体切片3、相干体的应用概述(1)展示断层发育目标层Ng相干体沿层属性平面图

北东向断层十分发育,并收敛于西南部目标层Ng相干体沿层属性平面图北东向断层十分发育,并收胜利油田潜山油藏顶面沿层相干切片+信噪比迭合图胜利油田潜山油藏顶面沿层相干切片+信噪比迭合图(a)常规地震数据体的相干切片;(b)波阻抗体的时间切片;(c)由波阻抗数据体计算的相干切片展示断层的常规数据与波阻抗数据的对比(a)常规地震数据体的相干切片;(b)波阻抗体的时间切片;(常规数据体与相干数据体t=2.4s的水平切片比较(据GeoEnergy公司资料)常规数据体与相干数据体t=2.4s的水平切片比较(据GeoE(2)研究礁体结构(2)研究礁体结构(3)检测裂缝发育带(3)检测裂缝发育带胜利油田罗家地区沿层相干切片

胜利油田罗家地区沿层相干切片地面地震时间-振幅切片相干体切片地面地震+相干切片展示断层和裂缝发育带的相干体切片2、相干体的应用概述地面地震时间-振幅切片相干体切片全方位和不同方位角的相干体切片

全方位和不同方位角的相干体切片(4)控制三维资料处理质量(4)控制三维资料处理质量(a)粗叠加处理(b)剩余静校正(c)偏移处理经过不同处理过程的常规数据体时间-振幅切片与相干体时间切片的比较

2、相干体的应用概述(a)粗叠加处理(b)剩余静校正(5)估计偏移速度场(5)估计偏移速度场(6)展示可能含气砂岩位置(6)展示可能含气砂岩位置(7)地质灾害预测2、相干体的应用概述指明浅水流动沙和直线状河道特征的相干切片(408ms)北西-南东向直线状河道,物源来自西北角的岩丘(7)地质灾害预测2、相干体的应用概述指明浅水流动沙和直线表明坍塌特征的相干切片(725ms)图中箭头所指处为具有冲刷痕迹的坍塌特征

表明坍塌特征的相干切片(725ms)首先在相干数据体上进行浏览,作小断层以及特殊岩性体的调查,了解其空间分布,这项工作不需要进行地震反射层位的解释就可实现。一般地说,高连续性数据对应连续的地层;中等宽连续性数据对应层序特征,如海侵/海退序列;窄条带低连续性对应断层、岩性的变化或特殊岩性体的边界;宽条带低连续性对应数据质量不好或无反射层位。然后对相干数据体切片进行解释,这种解释与常规解释思路不同,不需要先观察垂直剖面,只需在相干数据体切片上对不相干数据带进行解释。最后,要进行地质分析,搞清地层关系,分析工区内影响地震反射波连续性的因素,并结合地震纵测线、地质、测井资料对相干体数据进行综合解释。4、相干数据体解释的三大步骤首先在相干数据体上进行浏览,作小断层以及特殊岩性体的很显然,断面上一般没有反射同相轴,它是以不连续或无资料区为特征的,而这些不连续性或无资料的区域追踪是比较困难的,但也正是数据体因为断层的存在而产生的不连续性形成了追踪不连续性的相干体技术。5、断面的自动拾取第二讲地震属性提取与分析很显然,断面上一般没有反射同相轴,它是以不连续或无资料三维相干数据体的静态显示

三维相干数据体的静态显示相干体技术在追踪不连续性的断面时,一般有两种方法:其一是利用初始解释结果作为被追踪断面的种子拾取点的方法;其二是利用预先处理来创建突出不连续性的数据体,如差异法,微商法和多属性分析法等,以实现对不连续性进行追踪形成断面。后者中的差异法主要是利用振幅、三瞬、波阻抗等数据体来进行连续性的比较,具有很高的分辨能力和边界检测功能。微商法实际上是对差异法的补充,如对瞬时相位作空间微商计算,即可寻找相位横向上的突变点,而这些突变点常常与断层、地层边界和同相轴的转折有关联。第二讲地震属性提取与分析相干体技术在追踪不连续性的断面时,一般有两种方法:其一是断块静态立体显示第二讲地震属性提取与分析断块静态立体显示第二讲地震属性提取与分析

1、油气藏的地震异常特征(1)振幅异常—亮点,暗点,平点及AVO特征等;(2)能量异常—含油气储层与围岩相比具有较高的吸收系数,导致能量的变化;(3)频率异常—油气藏存在使地层频率响应发生变化,表现为低频趋势;(4)速度异常—储层岩石弹性性质的改变引起油藏部位上出现速度异常和纵横波速度比的改变;四、地震属性的物理意义第二讲地震属性提取与分析1、油气藏的地震异常特征四、地震属性的物理意义第二讲地震(5)时间厚度异常—储层部位速度的降低导致油藏底面反射波旅行时的延长;(6)地震记录特征异常—与油藏物性复杂分布相联系,在油藏部位可能出现记录面貌和反射特征的异常。据此,可采用一系列分析方法获取众多的地震属性参数。第二讲地震属性提取与分析(5)时间厚度异常—储层部位速度的降低导致油藏底面反射波沿层属性是沿目的层段开一时窗,对窗内的记录作自相关、功率谱、傅立叶谱、自回归及其它统计特征分析,来提取出相关的地震属性,下面从以下7类阐述这些属性(32种)的地质含义。(1)地震波频率属性;(2)地震波振幅或能量属性;(3)地震波波形属性;(4)地震记录自相关函数的属性;(5)分辨率属性;(6)信噪比属性;(7)其它属性。2、单时窗沿层地震属性的地质含义第二讲地震属性提取与分析沿层属性是沿目的层段开一时窗,对窗内的记录作自相关序号地震属性序号地震属性1功率谱上25%对应的频率23整波形面积2功率谱上50%对应的频率24整波形能量3功率谱上75%对应的频率25波形正半周面积4平均中心频率26波形正半周能量5振幅谱主频27波形正半周平均振幅6频带宽度28有效段均方根频率7记录能量29全谱平均频率8自相关函数第一零值点位置30全谱均方根频率9自相关函数第二零值点位置31有效段带宽频率10自相关函数极小值与主极值之比32全谱带宽频率11自相关函数旁极值面积与总面积之比33最大振幅之比12频谱二阶距34均方根振幅之比13时窗内波峰数35能量之比14相对分辨率36主频之比15信噪比37波峰数之比16L2模38自相关函数主瓣面积之比17峭度39自相关函数旁瓣面积之比18自回归系数比a2/a1;40复赛谱算法对数衰减率19自回归系数比a3/a1;41傅立叶谱算法对数衰减20自回归系数比a4/a1;42功率谱算法对数衰减率21自回归系数比a5/a1;43低频能量之比22最大振幅值

序号地震属性序号地震属性1功率谱上25%对应的频率23整波形(1)地震波频率属性

功率谱上25%对应的频率f25;功率谱上50%对应的频率f50;

功率谱上75%对应的频率f75;

平均中心频率fAV;

振幅谱主频f0;

有效段均方根频率fR;

全谱平均频率fCQ;

全谱均方根频率fRQ。上面8个属性主要反映了地层厚度、岩性及含流体成分的变化,常用来检测由于上覆地层异常如气饱和或裂缝存在所致的选频吸收,也可识别由于地层学特征、岩相等改变而引起的细小的频率变化。(1)地震波频率属性(2)地震波振幅或能量属性

记录能量ACF(0);

最大振幅值Amax;

整波形能量ENZ;

波形正半周能量EN;

波形正半周平均振幅AVA。上面5个属性反映了目标层内波阻抗、地层厚度、岩石成分、地层压力、孔隙度及含流体成分的变化。可用来识别振幅异常或用于层序特征分析;也可用来追踪地层学特征,如三角洲、河道、各种扇或特殊岩性体;还可用于识别岩性变化、不整合、气体以及流体的聚集等。第二讲地震属性提取与分析(2)地震波振幅或能量属性第二讲地震属性提取与分析(3)地震波波形属性

时窗内波峰数;

整波形面积SGT,Z;

波形正半周面积SGT;上面3个属性反映了目标层内波阻抗的变化规律、沉积层序、地层层理特征、古代剥蚀面、古构造特征、沉积过程及其连续性、沉积盆地的大小等。第二讲地震属性提取与分析(3)地震波波形属性第二讲地震属性提取与分析

(4)地震记录自相关函数的属性在一般情况下,自相关函数主极值幅度代表着记录段的能量;主极值宽度与记录的视周期有关,频率低的信号,主极值宽度大;频率高的信号,主极值宽度窄;旁极值的幅值和面积表示地震记录的重复性及延续时间的长短。当反射层具有薄互层结构,反射记录出现干涉现象时,自相关函数幅值和面积增大。这类属性有:

自相关函数第一零值点位置τ1;

自相关函数第二零值点位置τ2;

自相关函数极小值振幅与主极值振幅之比ACF(min)/ACF(0);

自相关函数旁极值面积与总面积之比S234/St。上面4个属性反映了沉积条件的稳定性、地层分界面的光滑度。(4)地震记录自相关函数的属性(5)分辨率属性

频带宽度fband;

相对分辨率Rλ;

有效段带宽频率fb;

全谱带宽频率fbQ。上面4个属性反映了地震记录分辨地层厚度的能力,通常在大套均匀岩层中夹有反射性能强的夹层,则在地震记录上可以见到在较平静的记录背景上,有少量清晰反射波出现,此时记录分辨率参数较大;当反射层为薄互层结构,相邻反射层反射波相互干涉,形成复杂波组时,地震记录分辨率参数较小,频带宽度较窄。第二讲地震属性提取与分析(5)分辨率属性第二讲地震属性提取与分析(6)信噪比属性信噪比也经常取作有效波平均振幅Aq与干扰的均方振幅之比:信噪比属性反映了地震记录上的干扰背景强弱,也可反映地质条件的变化,如油藏部分经常有地震记录规律性变差、信噪比降低这一特点。第二讲地震属性提取与分析(6)信噪比属性第二讲地震属性提取与分析(7)其它属性这是用数学方法定义的属性参数,无明显的物理意义,这类属性有:

频谱二阶距FBD2;

L2模;

峭度表征地层岩性纯度,是一种用数学方法定义的属性参数;

自回归分析使用AR模型,将地震记录表示为本身样点的线性组合,用多项式的系数来表达地震记录的反射特征。第二讲地震属性提取与分析(7)其它属性第二讲地震属性提取与分析3、地层平均吸收属性分析的方法有傅立叶谱、功率谱和复赛谱等。除了谱比法外,还可用特征值法,利用地震记录相邻时窗的自相关函数特征值、振幅特征值的变化来判断目的层吸收性质的相对变化。利用这些方法,最终提取了以下11个属性:(1)最大振幅Amax比;(2)均方根振幅Arms比;(3)ACF(0)比;(4)主频f0之比;(5)波峰数之比;(6)自相关函数主瓣面积S1之比;3、地层平均吸收属性分析的方法有傅立叶谱、功率谱和复赛谱等。(7)自相关函数旁瓣面积S234之比;(8)傅立叶谱算法对数衰减率δf;(9)复赛谱算法对数衰减率δK;(10)功率谱算法对数衰减率δP;(11)低频能量比。根据地层吸收性质与岩相、孔隙度、含油气成分等的密切关系,可以用它来预测岩性、预测砂泥岩分布,在有利的条件下可以用来预测石油和天然气的存在。第二讲地震属性提取与分析(7)自相关函数旁瓣面积S234之比;第二讲地震属性提取与4、地震属性体沿层属性的地质含义我们从众多的地震属性中选取以下9种属性体:(1)瞬时真振幅;(2)瞬时相位;(3)瞬时频率;(4)道积分;(5)道微分;(6)能量半时;(7)线积分;(8)噪信比;(9)均方根振幅比。第二讲地震属性提取与分析4、地震属性体沿层属性的地质含义第二讲地震属性提取与分析(1)瞬时真振幅:这是所选样点上各道时间域振动幅值,即为地震道数据的隐含表示。广泛用于地震资料的构造和地层岩性解释,常与其它振幅属性一起用于分离高幅区或低幅区,如亮点和暗点技术,也可用于地震相划分。(2)瞬时频率:定义为瞬时相位对时间的导数,用度/ms或弧度/ms表示。经常用来估计地震波在地下介质中的衰减变化;往往油气的存在引起高频成分的衰减,故可用这一属性检测油气。第二讲地震属性提取与分析(1)瞬时真振幅:这是所选样点上各道时间域振动幅值,即为地震(3)瞬时相位:表示在所选样点上各道的相位值,以度或弧度表示。主要用于地震资料地质解释中的同相位追踪对比,同时也增强了油藏内弱同相轴,不利之处是对噪音也有放大作用;最终成图的彩色色标应考虑到结果的周期性,即由于油气的存在经常引起相位的局部变化,所以这一属性常和其它属性一起用作油气检测的指标之一;也可用于测定薄层的相位特征,其横向变化与流体含量变化及薄层组合有关。第二讲地震属性提取与分析(3)瞬时相位:表示在所选样点上各道的相位值,以度或弧度表示(4)道积分:方法的计算简单快捷,计算的结果反映了波阻抗相对变化,而相对波阻抗数值的横向变化反映了储层特性、岩石物性、流体性质等的空间变化,所以该参数在储层横向预测和油藏描述中具有较好的应用价值。(5)能量半时:它定义为在给定的分析时窗内,计算能量达到1/2时的相对时间位置。能量半时可用来测定分析时窗内能量累积的速度,可以指示岩相和岩性的相对变化。第二讲地震属性提取与分析(4)道积分:方法的计算简单快捷,计算的结果反映了波阻抗相对(6)道微分:其实是地震道的另一种描述方式,它记录了相邻采样点之间的幅值之差。因而,可以用它作为提取其他一些属性的原始数据。道微分剖面则不存在受浅层干扰的问题。道微分处理后剖面的高频成分增多,可以更加细致地反映地层内部的细节,有利于地震相和沉积相研究。第二讲地震属性提取与分析(6)道微分:其实是地震道的另一种描述方式,它记录了相邻采样(7)线积分:这是将时窗内的地震波包络线的长度作为一种属性。这种属性的计算公式为:,式中,N是采样点数,T是采样间隔(ms),a是离散振幅值。S的大小体现了地震剖面上振幅和频率的变化,它可以反映岩相、岩性和地震波衰减的特征。通常S值的增大对应着高水位期的细岩性,而较小的S值对应着低水位期的粗岩性。剥蚀面或断层面上通常对应着高S值。第二讲地震属性提取与分析(7)线积分:这是将时窗内的地震波包络线的长度作为一种属性。(8)噪信比:计算三维地震数据体的噪信比三维属性,再按照目标层提取沿层噪信比属性,用于储层研究。噪信比主要反映分析数据的资料品质,沿层噪信比反映了分析时窗内目标层反射信息的稳定性,同相轴的连续性以及噪音所占的比值。(9)均方根振幅比:其沿层切片反映了特定时窗内地震波振幅的平均变化水平,其数

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论