基于叠前道集小波域Q值求取方法研究.pdf

《基于叠前道集小波域Q值求取方法研究.pdf》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、第50卷第4期石油物探Vo1．50，No．42011年7月GE0PHYSICALPR0SPECTINGFORPETROLEUMJu1．，2011文章编号：1000—1441(2011)04—0355—06基于叠前道集小波域Q值求取方法研究魏文，王小杰，李红梅(1．中国石油化工股份有限公司胜利油田物探研究院，山东东营257022；2．中国石油大学(华东)地球资源与信息学院，山东青岛266555)摘要：地层Q值是用于油气检测的重要参数，针对应用叠后地震资料求取Q值存在多解性的问题，探索了基于叠前道集求取Q值的Qvo(地震品质因子Q随偏移距变

2、化)方法，进行了数值模型试算和实际地震资料的应用，结果证明了QVO方法的正确性和有效性。该方法首先利用小波变换将经过动校正处理的CMP道集进行时频域转换，通过选择不同的时窗和与之相对应的尺度因子，利用小波域谱比法计算非零偏移距频谱比斜率值；然后根据频谱比斜率值随偏移距的变化关系，求得零偏移距处的Q值。模型试算提取的频谱比斜率与偏移距平方近似呈线性关系，与理论吻合，证实了理论的正确性。将Qvo方法应用于实际工区地震资料提取地层Q值，结果表明，含气区与非含气区的Q值差异明显，含气区的Q值较低，表现为强吸收。与叠后地震资料提取的Q值相比，Qv

3、o方法提取的Q值精度高，更能真实地反映地层吸收特征，提高了储层含气性检测的精度。关键词：叠前道集；QVO小波变换；频谱比法；斜率值DOI：10．3969／j．issrL1000—1441．2011．04．007中图分类号：P631．4文献标识码：A在储层描述时，地层Q值对岩性的变化具有的精度。很高的灵敏度，当地层含油气尤其是含气时表现出很大的异常值[1]，因此，如何准确求取地层Q值1基本原理是地球物理学家较重视的问题。现阶段国内外有很多估算Q值的方法[6。。。，广义上说，可以分为两1．1小波变换大类，一类是在时间域计算Q值，如振幅衰减法

4、、小波变换是傅里叶分析的继承和发展，是为了子波模拟法、上升时间法和解析信号法等；另一类克服傅里叶分析不能做局部分析的特点而提出的。是在频率域计算Q值，如中心频率偏移法、频谱模能量有限实信号厂()对于解析小波(￡)的拟法、频谱比法和拟合技术等。其中频谱比法计算连续小波变换为结果最优，上升时间法和解析信号法计算结果要差(，6)一I以l一{仁)(d)一些，说明频率域方法要比时间域方法稳定。但是式中：t为时间；窗函数(￡)称为母小波，其上面的上述求取Q值的方法均基于叠后地震资料，由于横线表示取复数共轭；，b分别为尺度因子和平移叠后地震资料缺乏丰

5、富的振幅和旅行时信息，因而一些细微的地层特征无法在叠后地震资料中得到因子，满足∈R，口≠0，b∈R；(Wcf)(n，6)称为小反映，特别是当油气藏的厚度远远小于地震分辨能波变换系数；ll一专是归一化因子，使不同尺度的力时，用叠后地震资料提取的Q值很难确定含气小波保持相同的能量。储层的准确位置。母小波()满足允许条件目前利用叠前资料求取Q值的方法比较少。Jf。。<+∞(2)1998年，Dasgupta等提出了一种基于CMP道集一。。l(l计算Q值的Qvo方法[1，利用这种方法求取的并且平方可积。实际应用中可根据具体要求选取不同类型的小Q值

6、精度高，油气预测效果显著，改善了基于叠后地震资料求取Q值精度的局限。我们引入QVO方法，利用小波变换进行时频域转换，应用小波域收稿日期：2010一O5—31；改回日期：2011一O3—02。谱比法基于叠前CMP道集估算Q值，并经过线性作者简介：魏文(198O一)，男，工程师，主要从事天然气地震地质综合研究。回归得到零偏移距处地层Q值，一定程度上消除基金项目：国家科技重大专项“渤海湾盆地南部精细勘探关键技术”了炮检距和非固有衰减的影响，提高了计算Q值(2oo8zxo5006—004)资助。356石油物探第5o卷波。地震资料处理时，通常要研

7、究信号的振幅谱和方程，由(6)式得零偏移距处斜率值：相位谱，这就需要选择复数型小波，如Morlet小波，s一T(8)．(￡)一{ee-毒(3)式中：s为零偏移距处的斜率值；T为零偏移距处式中：盯为调制频率，≥5。为了使小波与地震子地震波通过目的层顶、底界面时所用的时间差；Q波的波形特征相接近，可对Morlet小波进行修正，为零偏移距处的地层品质因子。修正后的Morlet小波在不同的偏移距处，斜率值近似为()一丌-{ee。。(4)s≈古+(9)式中：C为调幅因子，它控制了小波函数的长度。Morlet小波和修正的Morlet小波具有一个式中

8、：Q为零偏移距处的地层品质因子。共同的性质，即利用它们对信号厂(￡)进行小波变由(9)式可以看出，当假定传播时间的变化满换得到的小波系数的虚部就是其实部的Hilbert足经典的正常时差方程时，求得的斜率值随