基于Gabor变换反褶积技术在渤海某工区的应用研究

使用常规的Wiener反褶积必须假设震源子波在地层旅行过程中是平稳的即一成不变的,这个前提条件与实际野外地震资料采集差别较大,而基于Gabor变换反褶积技术考虑到地震能量的衰减、子波的形变等非平稳性特征.地震道在Gabor域可因式分解成三项即震源子波、衰减函数和反射系数,该技术设计POU窗函数,并利用此函数在Gabor域对地震信号进行局部时频分解.Gabor域反褶积算法在Gabor域通过除以衰减函数和震源子波的乘积来估算地层反射系数,然后再做Gabor反变换可求得时间域的地层反射系数.理论模型的测试和实际地震资料的应用均表明,与Wiener反褶积相比较,基于Gabor变换反褶积可补偿中深层的能量衰减并因此拓宽有效频带和提高时间分辨率.     

几种时频分析方法的比较和实际应用

地震信号往往是非线性非平稳的,时频分析技术能同时展示信号在时间域和频率域的局部化特征。本文研究4种时频分析方法:短时傅里叶变换(STFT)、小波变换、广义S变换和Wigner-Ville分布。在理论上运用雷克子波模拟地震记录进行时频分析对比,发现广义S变换具有相对较好的时频聚焦性以及较好的交叉项抑制性。为了验证这一理论,我们分别在两个工区做了实际的地震实验,分析并阐述4种时频分析方法的优缺点,同时也验证广义S变换的优越性。      

基于物理小波的频谱分解方法及应用研究

在地震资料频谱分解中,采用匹配地震子波的物理小波,依据地震信号的特征,用振幅、能量衰减率、能量延迟时间及地震子波的中心频率等四类参数构造基本小波,把地震信号分解在小波域,高频分量能够得到精细的刻画.本文以物理小波变换为工具, 给出了该变换中的核函数的选择方法,进而提出了基于物理小波变换的频谱成像方法.我们将此方法用于海上某油田河流相储层的描述,并与常规软件中的小波变换频谱成像结果进行了对比, 结果表明,本文提出的方法更能精细地刻画地质事件.     

基于零相位子波地层Q值估计

地震子波一般并非脉冲,而地层品质因子的提取公式大都是假设地震子波为脉冲时推导得到的,这在实际应用中误差较大,并且在进行频谱比时不容易确定所选择的频率区间的范围,从而影响方法的适用性.由此,本文讨论了地震子波为一般零相位子波时,基于地震波的传播理论和S变换理论,推导了S域中地层品质因子的估算方法,并通过Marmousi2...     

深度域叠前弹性波阻抗反演（英文）

叠前深度偏移的广泛应用产生了大量的深度域地震数据,提供了深度域地震处理的数据基础,并为探索直接基于深度域地震数据的储层弹性参数反演提供了可能性。由于深度域地震子波存在与速度密切相关的特征,直接提取深度域地震子波存在困难,本文将深度域速度变换与脊回归方法结合,实现了深度域地震子波的确定性提取,将提取的深度域常速地震子波构建成依赖速度的空变地震子波,基于加权叠加思想,提出了深度域地震记录合成方法。再结合非线性共轭梯度算法,建立了基于深度域地震数据和模型驱动的深度域叠前波阻抗反演流程,不再进行额外的时间域和深度域之间的数据转换。深度域合成地震数据与实际地震资料的叠前波阻抗反演取得了较好的结果,表明该方法是可行的。该方法可作为深度域正演分析和叠前弹性参数反演的基础方法,为深度域反演方法的稳定性、精度和效率的进一步研究创造了条件。     

基于S域谱模拟技术的时变子波提取方法研究

反褶积是叠前地震数据处理中的重要环节,反褶积效果的好坏很大程度上依赖于地震子波的准确性.早期的反褶积处理大多数都是基于Robinson提出的平稳褶积模型,即地震子波是时不变的,但实际上由于地下介质的吸收衰减作用,地震子波是随时间不断变化的,这说明要进一步改善反褶积,使用时变的地震子波是必要的,因此本文提出了一种从地震资料中直接提取时变子波的方法.具体地,首先对单道地震数据做S变换求出其时频谱,进而得到其时变振幅谱,然后利用谱模拟技术从求得的地震记录振幅谱中拟合出每一时刻的子波振幅谱,在子波是零相位假设的前提下,拟合出的时变子波振幅谱即是所求频率域的时变子波,本文最后利用正演的单道地震记录和实际资料分别验证了所求频率域时变子波的准确性.     

基于小波变换的薄层地震信号奇点的检测

通过小波变换的局部极大模可检测出奇点的位置,本文对反射地震信号小波变换进行了数值计算和聚焦处理,发现小波变换后奇点的位置仅受奇点实际位置及地震子波长度两因素制约,与子波形状无关．子波长度可以通过信号的小波变换本身求取,并进一步加以消除,从而实现了薄层地震信息的检测,得出薄层的地震分辨率可达1/32波长．实际资料处理表明,有效信号通过奇点分析可得到突出,剖面分辨率可得到提高．      

基于广义S变换的地震资料振幅谱补偿和相位谱校正方法研究（英文）

由于地下介质对地震波振幅的影响和地震波频散因素,地震波振幅和相位随时间、空间及频率的变化而发生改变,本文提出一种基于广义S变换的振幅谱补偿和相位谱校正新方法。该方法在S域中分为振幅谱补偿和相位谱校正两个步骤进行处理:振幅谱补偿是在地震记录可靠频带范围内恢复反射系数的振幅谱,其具体实现是在S域中利用谱模拟技术来拟合时变子波振幅谱,从而补偿由地层吸收所引起的振幅衰减;相位谱校正是消除子波剩余相位的影响,其具体实现是在S域中利用相位扫描来拾取随时间、空间和频率而变化的相位校正量,并由Parsimony准则来进行最佳相位判别。本文方法不需要直接求取Q值,能够适用于变Q值情况。理论模型和实际资料处理表明,该方法不仅能恢复地层反射系数的振幅谱,还可以有效消除子波剩余相位的影响,使子波接近或达到零相位,从而提高地震资料分辨率。     

时频域时变反褶积方法

常规反褶积方法具有很多局限性,往往要求地震子波是最小相位并且是平稳的,反射系数序列为白噪等等.本文将常规褶积模型扩展成时变褶积模型,通过S变换对地震记录进行谱分解,估计出震源子波和衰减因子,求出时变子波,在S域进行时变反褶积,再反S变换到时间域输出结果.此方法完全打破了常规反褶积方法的局限性,数值试验和实际资料处理均证明了此方法的有效性.     

地震资料处理中小波函数的选取研究

本文给出了常见地震子波的一个模拟公式，可以很好地模拟零相位及混合相位子波，在一定意义上也可以近似模拟最大相位及最小相位子波．模拟出的子波加上适当的修正项后满足允许条件，可用作小波函数．与Ｍｏｒｌｅｔ小波类似，在实际应用中这些修正项在一定条件下可以略去，文中对Ｍｏｒｌｅｔ小波作了改造，使其能更好地适应于地震资料处理．研究了反射波能量及噪声等干扰波在时间－尺度域的分布特征与所选基本小波的关系．提出用地震子波（或与地震子波相近的函数）作为基本小波，对地震资料进行去噪及分频解释的方法．最后用实例证明方法的有效性．     

面波压制的Ridgelet域方法

面波压制是地震数据处理中的一个重要问题. 常规的处理方法虽然能在一定程度上压制面波,但是在处理过程中只是单一的利用面波的一种特性,例如频率域滤波中利用面波与有效信号频率之间的差别,因此难以有效地压制面波. 利用Ridgelet变换可将原始地震记录拓展到（a,τ,p）三维空间,从而可以同时利用地震记录的视速度、时间和尺度域特性差别,实现有效信号与面波的分离. 文中通过理论合成记录及实际地震记录的算例,证实了基于Ridgelet变换的面波压制方法是有效且可行的.     

S变换谱分解技术在深反射地震弱信号提取中的应用

在深反射地震资料处理中,当来自深部的有效弱信号和噪声干扰频带差异较小且难以区分时,传统滤波方法的应用会受到限制.谱分解方法是一种使用离散傅里叶变换,基于信号的频率-振幅谱等信息生成高分辨率地震图像的方法,通常用来识别介质物性横向分布特征,处理复杂介质内频谱变化和局部相位的不稳定性等问题,包括定位复杂断层和小尺度断裂等.S变换作为一种新的时频分析方法,具有自动调节分辨率的能力,近些年来被广泛应用到勘探地震、大地电磁等数据处理中,逐渐成为地球物理方法中噪声压制的有效方法之一.与常规石油反射地震资料相比,深反射主动源地震为了探测深部结构信息,常采用大药量激发方式、长排列观测系统等,导致深部有效信号基本湮灭在噪声干扰之中.针对深反射数据特点,本文结合谱分解和S变换技术,首先设计了简单的脉冲函数实验数据,证实S变换方法的有效性,同时说明谱分解方法的效果受所用时频分析方法影响较大,而其中决定分辨能力的变换窗函数的选取尤为重要.在此基础上,分别应用到深反射地震资料的单道和叠加剖面实际数据上,对比分析了传统变换谱分解和S变换谱分解的应用效果,单道资料对比结果表明:相比传统谱分解,S变换谱分解方法具有自动调节分辨率的能力,能够精确的标定深反射地震资料中弱信号不同时刻的频率分量;叠加剖面资料应用结果表明:由S变换谱分解得到的剖面结果与其他谱分解方法结果整体上具有较高的一致性,同时清晰地刻画出原叠加剖面上被噪声湮灭的低频细节特征,提高了剖面的分辨率及同相轴连续性;对比结果明显看出,Gabor变换谱分解方法得到的结果同相轴较为破碎,分析原因认为这是由Gabor变换的时频分解方法的定长窗函数所致,窗口大小不会随着信号频率的变化来调节长度,只能在处理的过程中根据一定的记录长度范围选取窗函数参数,而S变换谱分解方法在窗函数的选取时,通过时变信号的局部频率特征自动调节窗口长度,能够更好的刻画各个频段的细节特征,在深反射剖面成像应用中效果尤为明显.本文结果表明S变换谱分解技术在深地震叠加剖面上的应用有效地提高了来自深部弱反射信号的信噪比和分辨率,并刻画出了叠加剖面上所不具有的低频细节特征,在实际深反射地震资料处理中能有效保护低频弱信号获得更好的成像效果.本文为深地震反射资料中弱信号的保护处理找到一种有效的方法.     

基于Dreamlet变换的地震数据压缩理论与方法

为达到更加有效地表示地震数据的目的,仅仅将地震数据当作普通的图像数据处理是远远不够的,地震数据中蕴含的地震波的运动学特性也应作为重要因素而被考虑到.本文讨论了利用Dreamlet变换方法实现地震数据压缩的方法,并针对地震数据本身所蕴含的频散关系特性进一步提出了多尺度Dreamlet变换压缩方法.Dreamlet变换由2个一维局部谐波变换的张量积构成,它在提供地震波场时间-空间局部化性质的同时可以保留波场的运动学特性.通过对二维SEG/EAGE叠前、叠后数据的算例说明了Dreamlet变换用于地震数据压缩的有效性.利用压缩后的数据进行成像的结果更表明,与Curvelet变换方法相比,Dreamlet与多尺度Dreamlet方法可以提供更高的压缩比;在相同压缩比的条件下,使用Dreamlet与多尺度Dreamlet方法压缩重建后的数据进行成像能更好地保留成像结果中的重要结构.     

多尺度奇异性检测及其在地震相分析中的应用

本文从地震道的奇异属性出发,利用连续小波变换求取地震道的小波变换模极大值曲线,并沿此曲线提取地震道的小波变换系数,称之为小波变换模极大值连线振幅.此属性不仅可以代表信号本身,而且可以最大程度的区别于相邻道,并且具有地震道多尺度的特征,即兼具时频域的特征.由此结合自组织神经网络,我们提出了一种新的地震相分析方法.通过模型合成地震记录实验分析,证明此方法是可行的,且对地质层位的解释误差具有一定的容许度.最后,将此方法应用于了实际资料,取得了良好的效果.     

井间地震资料精细解释方法研究与应用

井间地震资料分辨率很高,但是缺乏成熟的解释方法,本文针对井间地震资料的特点,借鉴地面地震解释技术对井间地震资料开展了精细解释方法研究.首先在深度域资料特点分析的基础上,利用区域速度和层析成像速度将深度域井间地震资料转换到时间域进行研究,并通过粗化降频建立起井间地震资料与地面地震资料之间的相似性联系,得到对井间地震低频反射特征的认识,结合井旁道合成地震记录标定建立起井间地震资料的反射特征,准确的解释了井间地震的反射层位.采用地震瞬时属性分析和时频属性分析得到了对沉积环境的认识,进一步通过波阻抗反演更好的认识井间储层特征.在实际应用中综合多种方法进行解释,并结合测井资料对比落实了井间砂体展布特征,得到对井间储层的精细描述结果,解决了两口井之间储层分布范围、横向连通性及厚度变化等问题,取得了较好的应用效果.     

基于保幅拉东变换的多次波衰减

为在去除多次波时有效保护地震一次反射波数据的AVO现象,给后续反演、解释提供准确的地震数据,本文提出了一种基于保幅拉东变换的多次波衰减方法,该方法是对常规抛物拉东变换的修改,把常规的稀疏拉东变换在拉东域分成两部分：一部分用于模拟零偏移距处的反射波能量,增加的另一部分用于模拟反射波振幅的AVO特性.该方法不仅考虑了反射波同相轴的形状,还考虑了反射波同相轴振幅幅度的变化,从而可把反射波信息进行有效转换,进而有利于多次波的消除,更好地恢复有效波的能量.在把地震数据由时间域转换到拉东域时,本文采用了IRLS算法实现保幅拉东算子的反演.模型数据和实际地震道集的试算分析表明,与常规拉东变换相比,保幅拉东变换在去除多次波的同时可有效保护一次反射波的AVO现象.     

基于S变换的软阈值滤波在深地震反射数据处理中的应用

深地震反射原始单炮数据是非平稳的弱能量反射信号,信噪比较低.如何提高信噪比一直是深地震反射数据前处理中的一大难题.S变换是一种适用于分析非平稳信号的时频变换方法.同其他分析时变信号的方法相比,S变换的基本小波不必满足小波在时间域均值为零的容许性条件,它的时频分辨率与分析信号的频率有关,且其在时间域的积分可以得到傅里叶频谱,其反变换也简单.因此,S变换容易表示深地震反射信号复杂的时频特性.本文在S变换的基础上,利用软阈值滤波方法对深地震反射数据进行处理,实验结果表明,该方法有效地提高了信噪比,压制了有效频带范围内的混频干扰,突出了弱反射信号,使得波组信息更加丰富,有利于连续追踪有效反射波组和识别薄地层,特别是提高了深部Moho界面反射层位的分辨率,为深地震反射剖面后续处理和准确解释奠定了基础.     

Fast hyperbolic deconvolutive Radon transform using generalized Fourier slice theorem

The hyperbolic Radon transform has a long history of applications in seismic data processing because of its ability to focus/sparsify the data in the transform domain. Recently, deconvolutive Radon transform has also been proposed with an improved time resolution which provides improved processing results. The basis functions of the (deconvolutive) Radon transform, however, are time-variant, making the classical Fourier based algorithms ineffective to carry out the required computations. A direct implementation of the associated summations in the time–space domain is also computationally expensive, thus limiting the application of the transform on large data sets. In this paper, we present a new method for fast computation of the hyperbolic (deconvolutive) Radon transform. The method is based on the recently proposed generalized Fourier slice theorem which establishes an analytic expression between the Fourier transforms associated with the data and Radon plane. This allows very fast computations of the forward and inverse transforms simply using fast Fourier transform and interpolation procedures. These canonical transforms are used within an efficient iterative method for sparse solution of (deconvolutive) Radon transform. Numerical examples from synthetic and field seismic data confirm high performance of the proposed fast algorithm for filling in the large gaps in seismic data, separating primaries from multiple reflections, and performing high-quality stretch-free stacking.     

基于Curvelet变换的地震资料信噪分离技术

在地震资料中,噪声干扰严重影响了有效信号的提取,为此必须进行信噪分离处理.本文提出一种基于Curvelet变换和KL变换相结合的软硬阈值折衷处理方法.首先对地震数据进行Curvelet变换,然后对各尺度系数选取适当阈值压制噪声干扰,再利用KL变换提取数据中的相干有效信号,最后重构得到去噪后的记录.经合成记录和实际地震资料处理实验证明,该方法与小波变换法相比较,更能有效进行信噪分离,提高地震剖面信噪比和分辨率.