一种小波域叠前地震道集吸收衰减参数估计算法制造技术

本发明专利技术利用地层对不同频率地震波衰减程度的差异，提出了一种时频域地震道集吸收衰减参数估计算法。具体算法是：对每个共中心点道集，在小波域提取某一时刻的高低频衰减曲线，以这两个单频衰减曲线围成的面积来表示该CMP点在这一时刻地震波通过地层的吸收程度。通过对整个时间窗逐时间进行计算，得到整个分析区域完整的吸收特性参数数据体。在吸收特性剖面上，以高吸收特性来预测烃类异常。实际资料算例结果验证了该方法的有效性。与叠后吸收衰减方法相比，该方法不受叠加处理环节对信号的影响，提高了地震资料估计地层衰减参数的精度。

【技术实现步骤摘要】
【专利摘要】本专利技术利用地层对不同频率地震波衰减程度的差异，提出了一种时频域地震道集吸收衰减参数估计算法。具体算法是：对每个共中心点道集，在小波域提取某一时刻的高低频衰减曲线，以这两个单频衰减曲线围成的面积来表示该CMP点在这一时刻地震波通过地层的吸收程度。通过对整个时间窗逐时间进行计算，得到整个分析区域完整的吸收特性参数数据体。在吸收特性剖面上，以高吸收特性来预测烃类异常。实际资料算例结果验证了该方法的有效性。与叠后吸收衰减方法相比，该方法不受叠加处理环节对信号的影响，提高了地震资料估计地层衰减参数的精度。【专利说明】
本专利技术涉及石油地球物理勘探
，属于地震资料解释范畴，是一种时频域叠前地震道集吸收衰减参数估计算法具体地说是。
技术介绍
地震资料解释是地震数据处理的重要环节，它是将地震资料转化成我们对勘探区地下地质情况的认识，并从中找到与油气有关的信息。在进行储层含油气性预测过程中，地层的吸收特性是一个非常重要的参数，它对岩性变化有着很高的灵敏度，特别当地层中含有油气时，它会表现出很大的异常突变特征。地震波在地层中传播时，由于散射、地层吸收等原因，能量会发生衰减，这种介质本身所固有的衰减特性通常用品质因子Q来描述。它是地震波传播一个波长距离后储能与耗散能的比率，与地层内部结构以及岩性和流体特征密切相关，一般来讲，Q值越大表明吸收越小，Q值越小吸收越大。由波传播理论可知，地震波的这种衰减特性还与频率密切相关，频率越高，弹性波衰减越快，导致地震波主频向低频方向移动。当储层岩石中含有油气时，储层具有低Q的特征，地震波在聚集了石油、天然气的储层中传播时，会发生非弹性衰减，对高频成分的吸收衰减更强，低频能量相对增强，因此利用高、低频率段的地震波吸收衰减特征可以间接预测油气存在以及分布范围。当前对于地震波衰减的研究主要集中在Q估算上:频谱比率法，升时法，时间域振幅衰减法，质心频率偏移法等。以上每一种估算方法都是建立在一系列的假设之上，而且都是基于叠后地震资料的。叠加模糊了振幅和旅行时信息，掩盖了道集中某些岩性和流体特征，一些细微的地层特征在叠后资料上是检测不到的，特别是当油藏本身的厚度远远小于地震的分辨能力时，由叠后地震资料提取的数据很难对产层做准确定位。相较于叠后资料，叠前地震道集含有更丰富的地层信息，对叠前资料的处理和解释，可以提高烃类预测精度。小波分析是上世纪80年代,Grossman和Morlet首次提出的一种信号时间_尺度或时间-频率分析方法，具有多分辨率的特点，能根据高低频信号特点自适应调整时-频窗宽度，在时域和频域都具有良好的局部化性质，完全解决了傅里叶变换局部特征差，短时傅里叶变换对高低频信号处理缺乏选择性的缺陷，自问世以来，在地震数据处理中得到了广泛应用。
技术实现思路
本专利技术的目的是提供一种利用小波变换对叠前地震道集吸收衰减参数进行估计的算法，并将该算法应用于叠前实际地震资料的数据处理中，根据吸收特性剖面的高吸收特性来预测烃类异常。本专利技术解决其技术问题所采用的技术方案是:基于地层对不同频率地震波衰减程度的差异，将叠前地震数据通过Morlet小波变换转换到时间频率域，即将时间域的全频地震信号分解为各个不同频率段的时间域地震信号，提取每个共中心点道集在某一时刻的高频和低频衰减曲线，以这两个单频衰减曲线围城的面积来表示该CMP点在该时刻地震波通过地层时的吸收程度。通过对整个时间窗逐时间进行计算，得到整个分析区域完整的吸收特性参数数据体，在吸收特性剖面上，以高吸收特性来预测烃类异常。本专利技术主要利用小波域分频计算瞬时振幅，然后利用某一时刻某个CMP道集高低频曲线围成的面积作为该CMP在该时刻的吸收衰减参数，利用了小波分析局部化优势，具有一定的抗噪性能，可准确刻画目的层高频衰减特征。与叠后吸收衰减方法相比，该方法不受叠加处理环节对信号的影响，提高了地震资料估计地层衰减参数的精度。该技术的实现过程为: I)采用人工激发并由检波器记录地震波，取得原始叠前地震数据，去噪保幅 处理，抽取叠前地震道集并对其做同相轴拉平处理； 将某个CMP道集s(.% t)经小波变换，得到其时频域数据体:TiCti /; t),其中， ?为接收器空间坐标(即偏移距)，/为尺度或频率。3)利用以下公式求得小波域CMP道集地层吸收衰减参数 固定时刻t = T0,从数据体『rjx, /;r)中抽取两个单频曲线p(x;/J:,)和,其中和= 分别表示I5时刻频率为厶和的高、低频曲线，则小波域CMP道集地层吸收衰减参数由下式得到:【权利要求】1.，其特征在于基于地层对不同频率地震波衰减程度的差异，将叠前地震数据通过Morlet小波变换转换到时间频率域，即将时间域的全频地震信号分解为各个不同频率段的时间域地震信号，提取每个共中心点道集在某一时刻的高频和低频衰减曲线，以这两个单频衰减曲线围城的面积来表示该CMP点在该时刻地震波通过地层时的吸收程度；通过对整个时间窗逐时间进行计算，得到整个分析区域完整的吸收特性参数数据体，在吸收特性剖面上，以高吸收特性来预测烃类异常。2.根据权利要求1所述的方法，其特征在于主要利用小波域分频计算瞬时振幅，然后利用某一时刻某个CMP道集高低频曲线围成的面积作为该CMP在该时刻的吸收衰减参数，利用了小波分析局部化优势，具有一定的抗噪性能，可准确刻画目的层高频衰减特征。3.根据权利要求1所述的方法，其特征在于该技术的实现过程为: 1)采用人工激发并由检波器记录地震波，取得原始叠前地震数据，去噪保幅处理，抽取叠前地震道集并对其做同相轴拉平处理； 2)将某个CMP道集经小波变换，得到其时频域数据体pa/d,其中，X为接收器空间坐标(即偏移距)，/为尺度或频率； 3)利用以下公式求得小波域CMP道集地层吸收衰减参数 固定时刻t = ,从数据体中抽取两个单频曲线和為)，其中和分别表示&时刻频率为厶和Z的高、低频曲线，则小波域CMP道集地层吸收衰减参数由下式得到:ρ: (r0) = V ((L ΑΛ)|-f:, h )|) 其中，i力CMP道集数据随偏移距变化的序号，M为该CMP道集中的地震道个数； 4)逐CMP道集进行上式的运算，会形成一条在t= L时刻的吸收衰减特 性曲线，由地层的吸收特性可知，当储层中含有油气时，该曲线会有明显的异常凸起，由此可检测烃类异常，确定油气藏范围； 5)最后，通过对整个时间窗逐时间进行计算，得到对整个分析区域完 整的吸收特性参数数据体。【文档编号】G01V1/28GK103984013SQ201410167091【公开日】2014年8月13日 申请日期:2014年4月24日 优先权日:2014年4月24日 【专利技术者】高亚力 申请人:浪潮电子信息产业股份有限公司本文档来自技高网...

【技术保护点】
一种小波域叠前地震道集吸收衰减参数估计算法，其特征在于基于地层对不同频率地震波衰减程度的差异，将叠前地震数据通过Morlet小波变换转换到时间频率域，即将时间域的全频地震信号分解为各个不同频率段的时间域地震信号，提取每个共中心点道集在某一时刻的高频和低频衰减曲线，以这两个单频衰减曲线围城的面积来表示该CMP点在该时刻地震波通过地层时的吸收程度；通过对整个时间窗逐时间进行计算，得到整个分析区域完整的吸收特性参数数据体，在吸收特性剖面上，以高吸收特性来预测烃类异常。

【技术特征摘要】

【专利技术属性】
技术研发人员：高亚力，
申请(专利权)人：浪潮电子信息产业股份有限公司，
类型：发明
国别省市：山东;37