Shearlet域稀疏约束地震数据重建

在地震数据处理流程中,通常对不规则的、稀疏的或者缺失的地震数据进行插值处理,通过插值方法来避免多次波的预测错误和成像假频等现象,使地震数据处理更加精准。Shearlet变换是一种多尺度变换,具有最佳的稀疏性、方向性以及局部化特性。将Shearlet变换与基于Landweber加速下降迭代方法结合起来对地震数据进行插值,在保证求解精度的同时提高了计算效率。信号和噪声在Shearlet域具有不同的分布特点,通过阈值法压制随机噪声,可提高算法的抗噪性。此外,采用jitter采样的方式,更好地压制了假频信息。理论和实际地震数据验证了该方法的有效性。     

Radon变换多次波压制方法及应用研究

多次波在地震资料中是普遍存在的。Radon变换是现行商业地震资料多次波压制处理软件应用最多的模块之一。这里在详细阐述Radon变换基本原理和离散采样计算的基础上,给出了理论模拟地震数据的多次波压制算法和计算结果。进而将Radon变换多次波压制方法应用于松辽盆地地震资料处理中,在对实际的CRP道集合理设计滤波切除函数后,应用本文的Radon变换算法,得到多次波压制的地震剖面。对比多次波压制前、后的叠前时间偏移地震剖面,压制多次波后地震剖面的信噪比很高,地质构造特征清晰,可应用于后续的地质构造解释,地震属性提取和油气储层预测。该方法的计算效率较高,算法易于实现。理论和实际地震数据的多次波压制结果表明,该算法具有高精度和实用性强的特点。     

海上拖缆随机采样抗混叠机理分析

目前海上拖缆地震采集只有道间距满足采样条件,而炮间距、炮线距和缆间距都不能满足采样条件,造成了地震数据的采样不足,导致频谱产生混叠效应。近年来发展了海上拖缆随机采集技术,为了分析该采集技术抗假频特性,以简谐信号和合成地震记录为研究对象,在空间规则采样的条件下,分别分析了充足采样和不充足采样情况下的频谱特征,分析了空间假频产生的原因。同时,在规则采样分析的基础上,对简谐信号和合成地震记录进行了随机抽样和频谱分析,研究了随机采样对空间假频的压制作用,数据分析结果表明海上拖缆地震随机采样可以降低频谱的混叠效应。     

拉东变换压制多次波在新疆巴楚地区的应用

多次波是地震勘探中一种常见的干扰波,巴楚地区浅层强反射界面导致了多次波异常发育,严重影响了地震数据的成像,甚至产生了构造假象,压制多次波对于新疆巴楚地区地震数据的构造精确成像具有重要的理论意义与实用价值。结合巴楚地区地质结构特征,对地震采集的原始数据进行深入剖析,同时根据多次波自身的特点,在试验了多种技术方法之后,认为拉东变换压制多次波方法能够在该区取得地震数据构造成像满意的效果。给出了在不同数据域准确识别多次波的几种基本方法和判断准则,同时详细阐述了拉东变换压制多次波的基本原理。通过工区实际资料的处理,系统地演示了该技术方法运用的过程、结果及关键点。     

海洋资料多次波组合衰减技术及应用

多次波的压制是海洋地震资料处理中的重要步骤。海洋地震勘探中的多次波一般都异常发育且较难压制,压制多次波的好坏直接决定了处理结果的质量。目前压制多次波的方法有很多种,但不同的方法基于不同的理论,具有不同的针对性。单一的方法压制多次波往往不能针对所有多次波,达不到理想的效果。笔者使用组合衰减多次波技术,综合了SRME、拉东变换、LIFT去噪对海洋实际地震数据的多次波进行压制,取得了较好的效果。     

南海南部深水多次波模拟分析与压制处理

南海南部陆缘构造复杂多变,海底崎岖不平,造成该区多道地震资料多次波十分发育且混杂复杂绕射波。为了识别深层有效反射信息、压制多次波,首先通过正演模拟分析深水条件下多次波的传播规律和特征,然后对复杂理论Pluto模型进行详细处理并试验了几种典型多次波压制方法,寻找到一套有效的多次波压制处理思路:通过表面相关多次波衰减(SRME)法降低多次波能量级别,运用τ-p域预测反褶积压制周期规律的近中道多次波,再在此基础上灵活运用高分辨率抛物线Radon变换对共反射点道集(CRP)进行中远道多次波压制。理论模型和礼乐盆地海区的实际资料处理结果表明,该方法能够较大程度上压制深水资料多次波并增强对深部有效信息的识别。     

基于地震信号相似性的S变换域多次波压制技术

多次波压制是海洋地震数据处理的关键,如何保护一次反射波的同时高效准确地识别并压制多次波是地震资料处理的重点。笔者在共中心点(CMP)叠加压制多次波方法的基础上,基于CMP道集地震信号与叠加道地震信号的相似性特征在S变换域进行多次波压制。首先,利用多次波速度动校正并对CMP道集数据减去叠加道数据完成第一次多次波减除,此时,波场中残余的多次波以随机噪声的形式存在;然后,对减除多次波的地震波场应用一次波速度动校正,并以一次波速度叠加道为初始模型对CMP道集S变换谱进行相似性滤波完成第二次残余多次波去除。针对理论数据和实际数据进行了有效性测试,结果表明本文所提方法能够较好地压制多次波,提高反射波的成像精度。     

面波压制方法与道间距关系的数值分析

野外无检波器组合的采集方式使得面波的压制成为单点高密度地震数据处理的关键。F-K滤波、τ-p变换是两种常用的压制面波的方法,这两种方法压制面波的效果与道间距有关。对面波速度、频率、射线参数与道间距关系的数值分析发现,道间距越小,压制面波的效果越好;大道间距会产生假频而使得面波难以有效地滤除。这从另一侧面也证明这两种方法可以更好地压制高密度地震数据中的面波。      

混合Radon变换地震噪声压制的应用

如何有效地从地震资料中剔除相干噪声一直是地震勘探十分关心的问题,尤其是多次波的压制和消除。笔者采用相似函数压制端点效应和截断效应的混合Radon变换从地震剖面中一次性分离面波等线性噪声,再应用自适应滤波技术在Radon域识别并剔除多次波,并使用双曲Radon反变换计算获得有效反射地震信号。理论模型及实测资料试算结果表明,该方法对多次波压制效果良好,同时,应用混合Radon变换可简化数据中噪声压制处理流程。     

逆数据域表面多次波压制方法

基于波动方程压制表面多次波的反馈迭代方法,推导出压制表面多次波的逆数据域方法。该方法利用了一次波和表面多次波在逆数据域的可分离性质。含表面多次波的地震数据变换到逆数据域后,表面多次波的能量在逆数据域形成能量聚焦区域,可通过简单切除能量聚焦区域的信号实现表面多次波压制,避免了基于波动方程预测减去方法中的自适应相减过程,完全损伤不到一次波的能量。单道地震记录和合成炮数据模型试算表明：逆数据域多次被压制方法简单易行,可有效压制表面多次波,并保持一次波振幅不受损失。与常规的基于波动方程的表面多次被压制方法（SRME方法）对比表明,逆数据域方法能更好地应对有效波和多次波同相轴相交的情况。     

一种加权抛物Radon变换地震波场重建方法

基于带限正反抛物Radon变换的基本原理和计算方法,笔者给出了一种加权抛物Radon变换（PRT）算法,该方法可应用于叠前地震波场恢复和重建,如缺失偏移距地震数据的插值外推和反假频重采样,其权系数完全由数据域来确定,也就是说在迭代计算过程中,该权系数也在变化。对于均匀采样的地震数据,笔者给出了近偏移距外推和波场的反假频重采样的计算结果,试算表明笔者给出的加权PRT方法可有效地应用于地震波场重建中,并具有算法稳健、计算精度高和效率高的特点。     

线性Radon域地震干涉层间多次波预测方法

为避免时空域地震干涉技术预测层间多次波过程中出现假象干扰,影响地震资料处理。笔者结合线性Radon变换与地震干涉技术,在线性Radon域中利用地震干涉技术预测层间多次波。通过数值模拟验证了线性Radon域地震干涉层间多次波预测方法的有效性,并将该方法应用到复杂速度模型中。结果表明,线性Radon域地震干涉技术可以有效减少时空域地震干涉技术预测层间多次波存在的假象,压缩数据体积大小,提高计算效率,并能够适用于复杂地质环境。     

抛物 Ｒadon 域地震干涉技术

当地震干涉技术应用在不完整数据体上时,缺失的地震记录会在合成的虚炮集记录上造成大量的假象。为减少假象的产生,笔者将抛物 Ｒadon 变换与地震干涉技术相结合,在抛物 Ｒadon 域中应用地震干涉技术。本文通过数值模拟验证了抛物 Ｒadon 域地震干涉技术的有效性,并将该方法应用在不完整数据体上。结果表明,抛物 Ｒadon 域地震干涉技术可以有效减少假象的产生,并能够压制缺失的地震数据的影响,提高虚炮集记录的质量。     

干涉波场重建衰减OBC层间多次波

由于OBC(ocean bottom cable)数据的采集方式为接收电缆铺设在海底,炮点在海面,常规的层间多次波压制方法不再适用,这给实际数据的处理带来了巨大的困难。本文利用地震波场干涉重建方法,对采集到的OBC数据进行波场重建,得到炮点和检波点都位于海底的虚拟炮集记录,与常规的在地表激发、地表接收的炮集一样,把OBC数据的层间多次波转化成常规的表面相关多次波和层间多次波,再进行预测。将其结果与共聚焦点(CFP)边界算法去除层间多次波的结果进行了对比,体现干涉重建技术在OBC数据层间多次波预测方面的准确性。     

对比抛物Radon正变换几种矩阵求解方法

抛物Radon变换法(Parabolic Radon Transform)在地震资料处理中有广泛的应用。PRT可对不同频率的地震数据解耦处理,这一特点使得抛物Radon变换的计算效率比双曲Radon变换有数量级上的提高。在频率域求解时,需要对每一个频率成份求解同样大小的线性方程组。求解抛物Radon正变换的计算方法主要有Levinson递推法、共轭梯度法、Cholesky分解法和直接矩阵求逆法。最小平方抛物Radon正变换所形成的矩阵具有Toeplitz结构,可采用Levinson递推法进行计算。高分辨率抛物Radon正变换所形成矩阵的Toeplitz结构被破坏,一般采用共轭梯度法或Cholesky分解法进行求解。这里详细推导了复Toeplitz矩阵的Levinson递推算法,并分别对求解方程的四种方法进行了讨论,最后给出抛物Radon正变换求解的数值算例,并对所给出的四种方程求解方法的计算效率及计算精度进行了对比。     

基于F-K和Radon变换的多次波衰减方法

多次波干扰是地震勘探中最重要的问题之一,其衰减技术也是地震数据处理的难题之一。本文主要介绍了多次波的形成、分类及其特征,并在ProMax系统上对合成地震记录的多次波进行了压制:首先采用合适的速度对时-空域的地震记录进行NMO校正;然后,分别用F-K变换和Radon变换方法将其变换到F-K域和Radon域,对多次波所占据的主要区域进行充零衰减;最后,通过对应的反变换将其变换到时-空域,从而达到多次波衰减的目的。通过在ProMax系统中的应用表明:不同方法对多次波的衰减效果各有优劣,组合滤波方法能达到更好的效果。     

各向异性 Ｒadon 变换在深水地震多次波压制中的应用

在深海多缆、长拖缆地震勘探的各种噪音影响中,多次波对深层的有效反射影响难以消除。笔者介绍了一种高分辨率各向异性 Ｒadon 变换方法,在 Ｒadon 变换积分路径中增加非椭圆率 eta 参数,将传统 Ｒadon 变换积分曲线由慢度、偏移距两参数改为慢度、偏移距及非椭圆率三参数,更好地刻画了长偏移距下反射非双曲时差。通过对某海域实际采集的多缆长偏移距地震数据进行压制多次波处理,与传统方法相比取得了较为理想的效果,为后续深层构造成像提供了基础数据。     

Attenuation of multiple waves using the continuous wavelet transform (CWT)

This paper describes how the continuous wavelet transform is used to filter multiple waveforms in both time and frequency domains. It is well suited to process the stationary signals, and it shows the signal in both time and frequency scales. This new approach was tested first on synthetic data and then on real data. The results obtained on both cases were good. The method consists of identifying the multiples on which we apply a normal move out using the multiple velocity law. The multiples will be aligned and the primary reflections will not be aligned. This operation allows locating the multiples in the time-scale domain. We compute the continuous wavelet transform (CWT for short) in order to focus on the patterns relative to seismic events. To filter the multiples, we define a zone with frequency and time bounds. These bounds are deduced from the projection of the seismic trace. Then an automatic mask is applied to the pattern to be isolated. Filtering in time–frequency domain is done by keeping only the wavelet coefficients that are outside the mask and assigning zero to the coefficients larger than a threshold amplitude inside the defined zone. The mask shape does not matter, which is not the case in classical filtering, where both the window size and shape play a key role. The mask is defined from three parameters: time, frequency, and the wavelet coefficients. To go back to the time domain, one has to compute the wavelet transform inverse of the trace. This procedure is repeated for all traces. To reset the traces to their initial positions, we apply the dynamic correction inverse with the same velocity law as the multiples. It turns out that the attenuation of multiples by the CWT works fine, in particular, the two identified multiples were quasi eliminated (Fig. 10).     

基于三参数小波变换的地震瞬时属性计算方法及应用

为提高地震瞬时属性的计算精度,针对薄层反射地震信号含有快速变化的振幅和频率分量的特点,研究了新的分析小波即三参数小波,基于三参数小波变换,提出了在相空间计算地震瞬时属性的方法。三参数小波有三个可调参数,对信号做小波分析时有很高的自由度,能够很好地匹配地震子波或给定的有效信号,三参数小波与BMSW小波或其他小波相比,具有更好的时域局部化性质。实际地震资料的应用效果表明,三参数小波变换地震瞬时属性比Hilbert变换瞬时属性有更高的信噪比和分辨率,基于三参数小波变换的地震瞬时属性分析方法是识别薄层砂体的有效手段。