常用地震属性的意义.docx

1、常用地震属性的意义常用地震属性的意义地震反射波来自地下地层，地下地层特征的横向变化，将导致地震反射波特征 的横向变化，进而影响地震属性的变化，因此，地震属性中携带有地下地层信息， 这是利用地震属性预测油气储层参数的物理基础。随着地震属性处理及提取技术的 大量涌现，属性种类多达 几百种，实际应用人员应用起来遇到了很大困难，迫切需要 按实用的角度，总结各地震属性参数与储层特征参数间的内在联系，为进一步研 究建立地震信息与储层参数之间的关系提供可靠的前提条件， 做到信息提取有方向、有目标。为了达到这一目的，首先按类别较全面总结了目前常用地震属性，从算法开 始，分析了各属性所表达的在地震波波形上的意义

2、，从正向上分析地震属性变化 与油气储层特征变化的关系，进而探讨总结了它的潜在地质应用。1、属性体、属性剖面这类属性是按剖面（或体）处理的，是一个体文件（或剖面文件），属性值对应 空间位置，即（x、y、to.属性值），可以用于常规地震剖面的方式显示与使用，常 用的属性有：相干体（方差体、相似体等）、波阻抗、道积分数据体，经希尔伯特变 换得到的瞬时属性体、倾角、倾向数据体等，这些属性体可以直接应用于解释，也可以 用解释层位提取出来转变为属性层，下表为常用属性体属性意义及潜在地质应用一览表。属性名称（BCM同义词）反射强度定又在解释屮旳应用属性特征XX kz 11V/ V L1V/ 11Streng

3、th ) 振幅包络( Amplitude Enve lope ) 瞬时振幅(Instaneous Amplitude ) REFLSTAN (缩 写)瞬时相彳立At() f2(t) h2(t)1用于振幅异常的品质分 析；用于检测断层、河 道、地下矿床、薄层调 谐效应；从复合波屮分 辨出厚层反射进行地震地层层序和提供声阻抗差的信息。横向 变化常与岩性及油气聚集有关。 值总是正的描述了复相位图中实部和(InstaneousPhase)INS TP HAS (缩写)、 1 h( t)(t) tanf(t)特征的识别；加强同 相轴的连续性，因此 使得断层、尖灭、河 道更易被发现。可对 相位反转成图，有

4、可虚部之间的角度。它的值总 在180。之间。瞬时相位 是不连续的，从+ 180到 -180的反转可引起锯齿状 波形能指示含气否瞬时频率(InstaneousFrequency )INSTFREQ (缩 写)d (t)W(t) dt用于气体聚集带和低 频带的识别；确定沉 积厚度；显不尖灭、 怪水界面边界等突变 现象嶙时相位对时间的变化率。值域为(一 f W, + f w)。然而， 大多数瞬时相位都为正。 可提供同相轴的有效频率吸 收效应及裂缝影响和储层 厚度的信息正 交 道(QuadratureTrace ) 希尔伯特变换( HilbertTransform) QUADRATR (缩 写)h(t

5、)是f(t)的希 尔伯特变换，也是 f(t)的90相移 h(t) 1 * f (t )t用于复数道分析的品质控制当实地震道代表地震响应 屮质点位移的动能时， 正交道相当于质点位移的势能视极性( ApparentPolarity) APPAPOLA (缩 写)在振幅包络峰值处实地震道的极性用于振幅异常的品质 分析为实地震道的符号位，假设零相位子波、视极性与反射 系数的极性相同响应相位( ResponsePhase )RESPPHAS (缩写)在振幅包络峰值处的瞬时相位值地震地层层序的识 另I、检测。由于流体 含量或岩性引起的横 向变化，在具有相似 的振幅响应时，用来 区分有利和不利带强调反射界面

6、的主相位特 征。与瞬时相位的应用相同响应频率( ResponseFrequency )RESPFREQ(缩写)在振幅包络峰值 处的瞬时频率值识别与气藏聚集有关的可能区带相应频率在区域上更具可 解释性。与瞬时频率的应用 相同反射强度交流 分量(Perigram )PERIGRAM (缩 写)消除了反射强度中 的均值(直流分量) 部分后的偏差用于振幅异常的品质 分析。与反射强度的 应用相同，但更适合 于分析和处理，因为 它有正负这种显示使能量最大值的 定位比在地震剖面上更明 显、更清晰相位余弦(cosine ofPhase)瞬时相位 余弓玄 (cosine of Instaneous Phase)

7、道归一 化(Normalized Trace )f(t)C0S )A(t)用于地震地层层序和 特征的识别，与瞬时 相位的应用相同，但 克服了相位反转的跳 断，可用于数据加强 处理它在正值和负值之间平滑 地振荡。它可能影响地震显 示屮同相轴的外观，更便于 用传统的彩色图进行分析反射强度交流分量相位余弦( Perigra mcosine 八酋Phase)刀里 GRPXPER1 (缩当 Perigram0时，反射强度交流亍相位余弦乘积；否则为0写)强振幅、连续相位成 图，用于振幅异常分 析，与反射强度应用 相同将实际资料分离成振幅 (Perigram)和相位(cosine of phase)两部分,

8、消除小 于振幅能量一半的数据相干体计算相邻地震道的互相关系数识别断层、裂缝带、 河道和砂体边界等时窗长度可以选择， 还可以 选连续度处理和非连续度 处理。另外还有相干系数的 平均、均方、屮值等选项。相似体计算相邻地震道 的相似系数它将地震资料、测同上不但可以对三维体数据作 不连续分析，还可以对基于 层位的二维数据作相似性 预测，以及倾角、方位角 ，边界检测和图象增强。 还可 以沿层解释的层位作相似 性分析波阻抗井数据、地质解释 相结合，利用测井 资料具有较高的 垂向分辨率和地 震剖面有较好的横向连续性的特 点，将地震剖面“转换成”波阻抗 剖面用于储集层的研究，识别砂体的分布特征和范围将地震资料

9、与测井资料连 接对比，能有效地对地层物 性参数的变化进行研究，对 储层特征进行描述道积分对地震道进行积 分识别砂体、岩性尖灭点等相对对数波阻抗倾角倾向数据体计算同相轴的倾角识别尖灭点、不整合、了解地层产状2、沿层地震属性这种属性是用解释层位在地震数据体(剖面)中提取出来的属性，它的数值对应一个层位或一套地层，每个属性值对应一个x、y坐标。提取方式有两类：沿一个解释 层开一个常数时窗，在此时窗内提取地震属性，提取方式有4种(图2-la ) o用两个解释层提取某一段地层对应的地震属性，提取方式也有 4种(图2lb ) o常用地震属性的计算方法总结如下：(1)、均方根振幅 (RMS Amplitud

10、e )均方根振幅是将振幅平方的平均值开平方。由于振幅值在平均前平方了，因此， 它对特别大的振幅非常敏感。(2)、平均绝对值振幅KMS =丄(32 +奸+ H72 + 462) 16=烏945)二 J5246.56=72.43(Average Absolute Amplitude )Average Absolute = 口131 日bsoluiE rrp讦udc?Amplitude number of samples二 1045/16=65.31平均绝对值振幅没有均方根振幅那样，对特别大的振幅敏感。(3)、最大波峰振幅(Maximum Peak Amplitude )最大波峰振幅的求取方法是，对

11、于每一道，PAL在分析时窗里做一抛物线，恰好通过最大正的振幅值和它两边的两个采样点，沿着这曲线内插叮得到最大波峰值 振幅值。PAL画一个使这三个采样点适合曲线并且 07 沿这一曲线确定出最大值。imum Peak Amplitude = 125(4)、平均波峰振幅(Average Peak Amplitude)平均峰值振幅是对每一道在分析时窗里的所有正振幅值相加，得到总数除以时 窗里的正振幅值采样数得到的。(5 )、最大波谷振幅(Maximum Trough Amplitude)最大波谷振幅的求取方法是，对于每一道，PAL在分析时窗里做一抛物线，恰好 通过最大负的振幅值和它两边的两个采样点，沿

12、着这曲线内插可得到最大波谷振幅值。1276Q7imum Trough Amplitude = I-90I = 90117波谷振幅(Average Trough Amplitude)平均波谷振幅是对每一道在分析时窗里的所有负振幅值相加，得到总数除以时 窗里的负振幅值采样数得到的。(Maximum Absolute Amplitude)计算每道的最大绝对值振幅的求取方法是，首先在分析时窗内计算岀波峰和波谷 的值，得出最大的波峰或波谷值，然后，PAL画一抛物线，恰好通过最大波峰或波谷 振幅值和它两边的两个采样点，沿着这曲线内插可得到最大绝对值振幅值。画一个适合这三个采样点的曲线 并且沿着这一曲线确定

13、出最大值。Maximum Absolute Amplitude = 123.6(8 )、总绝对值振幅(Total Absolute Amplitude)总绝对值振幅是计算确定时窗内的所有道的绝对值振幅值。Total Absolute Amplitude = sum of absolute values of(9)、总振幅(Total Amplitude)每一道的总振幅是，在层内对采样点求取总的振幅值。Total Amplitude = sum of amplitudes=559(10)、平均能量(Average Energy)对于每一道的平均能量的求取方法是，对分析时窗内的振幅值平方相加，对总

14、数除以时窗内的采样数求得。*9482sum of squared amplitudesAverage Energy = number of samples、7687= 83,945/1638-86-17(11)、总能量(Total Energy)对于每一道总能量的求取方法是，对分析时窗内的振幅值平方相加求和得到的。Total Energy = sum of squared amplitudes,945(12)、平均振幅(Mean Amplitude)对于每一道的平均振幅的求取方法是，对分析时窗内的振幅值相加，总数除以 非零采样点数得到的。、122 A . . sum of amplitudes

15、Mean Ainphiude = ! number of non-zero values= 559/16七6-94 1225広12 7687J-1176S =召工(“-对=范(叫 _ 34.94 )5=亠】,74& 392.31 16=109, 274.52(15 )、振幅的峰态(Kurtosis in Amplitude)对于每一道的振幅的峰态的求取方法是，对分析时窗内的所有采样点求取平均 值，然后减去每道的平均值，计算差值的四次方，求出这些值的总和，除以采样点 数就可得到。(16 )、有效带宽(Effective Bandwidth )数据体时窗的有数带宽是曲数据体的零延时的自相关除以采样

16、/所有自相关的总和之积而求得的。r(t) = the two-sided autocorrelation of the data in the windowT = sample periodER =广(0)+ t 工 r(fOWirfiiBwAengti = M+l有效带宽被看作是定量化的相似数据体。狭窄的带宽就是比较相似的数据体； 而较宽的带宽是不太相似的数据体。因此，宽的带宽表示不均质的反射特征，被认 为是复杂的地层；窄的带宽表示的是较简单的或平滑的反射特征，认为是均质的地 层模式。带宽能帮助我们在数据体中识别噪声区，有噪声的数据体比没有噪声的纯 数据体有很明显宽的带宽。应用地震地层学的方

17、法，可以从与其它属性相配合的有效带宽中推断出一系列 地震反射所代表的沉积环境。如一个狭窄的带宽区域，低振幅，高频，连续的平行 反射代表了低能量沉积环境，认为是深海页岩。(17)、弧长(Arc Length )弧长是作为地震道的波形长度来定义的，它是在时窗内对所有地震道的变化范围的比例测量。假想，用道的波形样式绘制地震道曲线，然后想象一根绳子放在地 震道上跟着每个波形波动。地震道的弧长就是当绳子伸展开的总长度。N s =馬工皿+1)-曲+ T2 r = 1N = number of samples in the window弧长是用于高振幅高频率和高振幅低频率之间与低振幅高频率和低振幅低频率 之

18、间的区别。如因为页岩和砂岩的界面，一般有一些突变和高阻抗的反差，弧长就 用于页岩层序和含砂量较高的层序之间的识别，带宽越小，弧长就越接近总绝对值 振幅。这一属性相似于反射的非均质性。(18)、过零值平均频数(Average Zero Crossings Frequency )过零值平均频率的计算方法是通过数据体时窗中的过零点的个数(Nzc),和求出 第一个通过零值的反射时间和最后一个通过零值的反射时间，根据下式计算出过零点平 均次数(fzc ) O一】 t 1 = time of first zero crossingf .= * J 2(&=-t4ifie of last zero cros

19、sing对于过零值平均频数的用途相似于瞬时频率，由于它不会有尖脉冲，并且它的 值不会为负值或无穷大，因此它是一个比较稳定的量。当时窗比较小时，过零值平 均频数对波形中较小的变化比平均瞬时频率敏感。(19)、Dominant Frequency Series Fl、F2、F3 (主频系列 Fl、F2 F3)对于所确定时窗的每一个输入道的估算值是由能量谱中的三个最主要频率分量组成，如下图中的Fl、F2、F3o其中F1是低频段中的峰值，F2是中间频段中的峰 值，F3是高频段中的峰值。运行这些属性，PAL就会用最大嫡方法，对每道进行谱 分析，六次多项式是用于能量谱模式和识别它的三个峰值。它应用的优点是

20、能够输入有限的数据得到可靠的估算值。对于一定的输出格式必须由40ms的数据，当分析时 窗在40ms以下时，PAL将会输出无效值。上图所绘的能量谱图是通过对所有道进行快速傅立叶变换得到的，主频估算值是在50ms的分析时窗中得到的。最大嫡方法是在有限的时间时窗内得到可靠的估算 值，但这些是对三个主频的数学方法估算值，并且这些估算值可能不总是于与你在实 际能谱上看到的峰值一样。这三组属性帮助你在数据时窗内来确定主频特征，在任意或所有主频系列属性 里的侧向变化可能有由油气饱和度或断裂导致的频率吸收效应的特征。例如，油气 饱和的砂体削弱了较高的频率，这样你就会看到较低的一个或所有的主频。虽然同样的是计算

21、峰值谱频率，因为它可以显示在振幅谱中的最重要的三个点, 所以主频系列有更多的信息。通过更多的振幅谱特征，主频系列可以揭示与地层或岩 性有关的频率趋势。(20)、峰值谱频率(Peak Spectral Frequency )对于所确定时窗内的每一输入道，峰值谱频率的估算值是由能量谱中单一的最 主要的频率组分组成。峰值谱频率相似于主频系列，主频系列估算值是由能量谱中 的三个最主要的频率段组成。大体上，峰值谱频率将描述的是主频系列(Fl、F2、F3) 中所给任意道的最主要的谱组分。运行这些属性，PAL就会用最大炳方法，对每道进行谱分析，多系数多项式是 用于能量谱模式和识别它的最重要的峰值谱频率。它的

22、应用的优点是能够输入有限 的数据得到可靠的估算值。对于一定的输出格式必须由40ms的数据，当分析时窗在 40ms以下时，PAL将会输出无效值。上图所绘的能量谱图是通过对所有道进行快速傅立叶变换得到的，主频估算值 是在50ms的分析时窗中得到的。最大嫡方法是在有限的时间时窗内得到可靠的估算 值，但这些是峰值谱频率数学方法估算值，并且这些估算值可能不总是于与你在实际 能谱上看到的峰值一样。峰值谱频率提供了一种追踪主频特征的方法，主频特征可能由油气饱和度、断 裂、岩性、地层的变化有关现象导致的频率吸收效应所带来的的特征。例如，油气 饱和砂体吸收了较高的地震频率，这样你可能看到较低的峰值谱频率值。任何

23、大于门槛值的频率都将从峰值谱频率分析中被排除的。在数据体中设定门 槛值为最大有效频率，一般来说，这个值是信噪比为 1的频率值。在这个频率值以(21 )、从谱的峰值到最高频率的斜率 (SpectmlSlopefiomPeakto MaximumFrequency)这个属性表明了在分析时窗内高频成分被吸收的特点。你确定了一个感兴趣的 最大值，PAL就计算出在谱中的峰值频率到你设定的门槛值衰减比率。如果斜率 是一个高值，高频成分很快被吸收；如果斜率是一个低值，就没有信息被吸收。对每一输入道,PAL会用最大爛方法计算峰值谱频率，多系数多项式是用于能 量谱模式和识别它的最主要的峰值。这个过程是用最小二次

24、方回归法确定一个线性关系，适合于在峰值频率和对于 谱估算的最大频率之间的所有能量谱模式的线性关系。斜率用 db/HZ表示。下图实例中，所绘的能量谱图是通过对所有道进行快速傅立叶计算得到的，波 峰谱频率和从波峰到70HZ最大频率的斜率是用PAL在50ms的分析时窗中得到的。这个属性想通过对能量谱的衰减的估算，用频率在典型的能量谱内定量表示频 率的吸收效应。谱斜率的侧向变化可能由于油气饱和度或断裂或与岩性或地层的变 化有关现象导致的频率吸收效应所带来的特征。例如，油气饱和砂体衰减了较高的 地震频率，谱斜率就会比较陡峭。谱估算的最大频率值用于规定了峰值谱频率的上限，也就是，当峰值谱频率计算 出来时，

25、任何高于这个门槛值的频率都会被排除。在数据体中设定门槛值为最大有效 频率，一般来说，这个值是信噪比为1的频率值。在这个频率值以上，存在的噪声比 信息多。因此这个数据不会对整个时间道有建设性的作用。(22)、大于门槛值百分比(Percent Greater than Threshold )对于每一道来说，在分析时窗中，大于设定的门槛值的采样个数除以总采样个 数，乘以100o当在时窗内振幅属性采样率为平均采样率时，为求得的振幅大值或小值所占比 例，大于门槛值的百分比就决定了大于设定的振幅门槛值的采样数的多少。在某种 意义上，你所计算的主要是在时窗中的相对高振幅部分。这一方法的优点是，它是 对某一层

26、统计计算，并且对数据体特征中的侧向变化是非常敏感的。小于门槛值的百分比主要用于分析地层的延伸，海进和海退垂直序列层序会在 高振幅砂岩面和低振幅页岩面之间产生。通过计算大于门槛值的振幅百分比，你可 以确定这些垂直变化和绘出横向变化的范围图。同样的，这一属性可以帮助你区分出整合基底（高振幅）、丘状起伏基底（较低 振幅）和杂乱反射基底（低振幅）之间的不同。另一个应用是在层序或沿确定的反射层内可以画出异常振幅图，例如由油气或 流体的聚集，不整合和调谐效应所导致的异常。（23 ）、小于门槛值百分比（Percent Less than Threshold ）对于每一道来说，在分析时窗中，小于设定的门槛值的

27、采样个数除以总采样 个数，结果乘以100o当在时窗内振幅属性采样率为平均采样率时，为求得的振幅的大值或小值所占 比例，小于门槛值的百分比就决定了，小于设定的振幅门槛值的采样数的多少。在 某种意义上，你所计算的是在时窗中的相对低振幅部分。这一方法的优点是，它是 对某一层的统计计算，并且对数据体特征中的侧向变化是非常敏感的。小于门槛值的百分比主要是用于地层走向方面的。在特定的第三纪盆地内，三 角洲层序是从富含砂，高均方根振幅，到富含页岩前三角洲或深海平原里面的低振 幅来划分的。这些油页岩比率的变化通过看图中的小于门槛值的百分比就可以很容 易确定。振幅）和杂乱反射基底（低振幅）之间的不同。另一个应用

28、是在层序或沿确定的反射层内可以画出异常振幅图，例如由油气或 流体的聚集，不整合和调谐效应所导致的异常。（24）、能量半衰时（Energy Half-Time ）在研究的时窗内，从上到下根据样点数求能量累加之和。当能量之和达到计算时 窗内总能量的一半时，到这点的样点个数除以总的样点个数为这点的能量半衰时。能量半衰时是在一个周期内时间达总时间的一半测量时间所需要的能量。它用这个周期的时间域的百分数来表示。如果在分析时窗内振幅是相对一致的，那么总能量的一半就会在时窗中心附近（能量半衰时二40%-60%）；如果在时窗中较浅的部分是强振幅，那么它就会用较少的 时间到达总能量的一半（能量半衰时=10%-4

29、0%）；相反的，如果在时窗中较深的部分 是强振幅，那么能量半衰时就会较长（能量半衰时 二60%90%）。下面是一道的例子:8394583945/241972.5 ( between 8th & 9th samples表示的是地层的变化或由流体含量、不整合或岩性有关所造成的振幅异常。例如，海进和海退层序常常具有高能的砂岩的反射和低能的页岩的分布变化特征。当从页岩向下到砂岩层序分层时，能量半衰时将大于50%o当从砂岩向下到页岩 层序分层时，能量半衰时将小于50%o能量半衰时中的横向变化图可以帮助整个地层 解释。能量半衰时也能对振幅异常描述由帮助。例如，亮点和暗点与油气含量有关， 当这些异常在分析时窗内改变了能量的分布时，你可以看到能量半衰时中的变化。 对于可以被检测出来的在中心的能量分布的时移，时窗必须包括最前或尾部的数据 体作为异常振幅的参考。（25 ）、能量半衰时斜率（Slope at Energy Half-Time ）能量半衰时斜率所计算的是当所累计的能量是总能量的一半时所需时间的能量曲线的斜率。Slope at Energy Half-Time= E(nhalf) E(nhalf-l)E (energy) = amplitude squared of the