基于小波谱能量曲线EWT的面波压制算法

面波提取和压制是地震资料处理的关键问题之一,对地震资料处理的准确性具有重要影响.针对传统地震面波压制算法面波分离精度低、有效波损伤大的问题,提出一种基于小波谱能量曲线(Wavelet Spectrum Energy Curve,WSEC)经验小波变换(Empirical Wavelet Transform,EWT)的面波压制算法,该算法利用连续小波变换(Continuous Wavelet Transform,CWT)精细化的时频分析能力得到地震信号的WSEC,从频率和能量两个维度对信号进行自适应分解,提高面波提取精度.首先对地震记录进行CWT,根据小波谱计算各频率点的能量得到WSEC;基于WSEC极大值点的频率及ε邻域法确定频谱分割边界进行改进EWT,得到依能量和频率分解的各固有模态分量(Intrinsic Mode Function,IMF);根据各子频带的频率和能量确定面波IMF,最后对面波IMF进行带通滤波以保护其中的有效波.仿真实验表明,算法能根据地震信号的频率和能量自适应分解,实现面波的精准分离与提取;实际地震资料处理表明,算法能准确分离面波,能最大限度保护有效波,提高资料信噪比,是一种有效可靠的面波压制算法.     

面波压制的Ridgelet域方法

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

多分量地震数据矢量滤波面波压制方法（英文）

针对陆上三分量地震资料去噪处理中面临的面波噪声能量强、有效信号能量弱、转换波与面波主频叠置、非矢量保真处理等问题,本文基于多分量数据中有效信号与强面波噪声在极化属性、振幅能量、视速度三方面的特征差异,提出了一种时频域压制多分量数据强面波的去噪方法,并设计了相应的矢量滤波器。该方法利用小波变换将地震信号变换到时频域,并基于时频域瞬时极化率公式计算出各点极化率参数,然后根据面波视速度划分出有效信号区和面波噪声区,并据此对有效信号进行能量均值统计,最后联合极化率、能量、频率特征对噪声识别和分频、分区压制。理论模型和实际资料处理验证了本文方法的正确性和有效性,为解决低信噪比陆上多分量资料面波压制提供了又一种有效手段。     

多方向正交多项式变换压制多次波

提出一种基于Radon 变换和正交多项式变换的多方向正交多项式变换压制多次波方法.抛物Radon变换对不同曲率方向的同相轴叠加,根据速度差异区分一次波和多次波,但Radon反变换会损伤振幅特性,不利于AVO分析.多方向正交多项式变换在Radon变换(某一曲率方向的零阶特性)的基础上,利用正交多项式变换进一步分析同相轴的高阶多项式分布特性,用正交多项式谱表征同相轴AVO特性;根据一次波和多次波速度差异和同相轴能量分布特征实现多次波压制.该方法的优点是仅用一个曲率参数就可描述同相轴剩余时差参数,提高了一次波和多次波的剩余时差分辨率.实验结果表明,该方法可以有效压制多次波并保留一次波AVO特性.     

基于第二代Curvelet变换的面波压制(英文)

针对面波能量强、频率低、视速度低等特点,本文提出了基于第二代Curvelet变换的多尺度、多方向分解的面波压制方法,充分利用Curvelet变换的多尺度、多方向功能,分离出含有面波的尺度和方向然后利用面波与有效波在Curvelet域基本不重叠的特性,进行信噪分离处理。这种方法在实际叠前地震数据的处理中,能够有效压制面波,同时较好地保护有效反射波,特别是反射波的低频分量。文中提出的基于第二代Curvelet变换的去面波方法是一种较好的保幅去噪方法。     

基于小波分频与径向道变换的联合压制面波方法（英文）

面波是地震记录上严重的干扰波。基于小波变换的面波衰减方法考虑了面波频率低的特点,但是当有效信号和面波干扰存在频率重叠时,其压制效果不理想。径向道变换考虑了面波与有效信号视速度差异,有利于去噪,但是在去除面波的同时也会损害有效信号的低频分量。本文根据面波和有效信号在视速度以及小波域能量的差别,将小波变换的局部分析能力与径向道变换的去噪优势相结合,提出了基于小波分频与径向道变换的联合面波压制方法。首先应用小波变换,将地震记录分解为不同频段,对出现面波的频段作径向道变换,然后再作低切滤波处理,最后利用径向道反变换后的记录与其他频段的记录进行小波重构得到去除面波的记录。两种方法的联合使用提高了小波分频去除面波的能力,同时也较好的保护了有效信号。通过模型数据实验分析和实际资料数据处理的结果表明,本文提出的方法具有较强的去噪能力和良好的保幅性能。     

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

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

保幅低频面波压制方法与应用

面波是陆地地震勘探中的一种强能量规则干扰,面波的存在对地震数据成像有很大影响,因此需要对面波予以压制.常规面波压制方法只考虑面波的强能量特性,很容易损伤面波频带内的有效信号,难以满足地震资料保幅处理的要求.在传统方法的基础上,本文提出低频弱信号保护的面波压制方法,该方法在有效地消除面波干扰的同时,最大程度的减少了低频有效信息的损失,提高地震资料的信噪比和保幅性,松辽盆地北部西部斜坡地区实际地震资料的处理表明,该方法是一种有效的相对保幅的面波压制方法,可为岩性油气藏勘探提供重要的基础数据,具有很强的实用性.     

基于希尔伯特-黄变换的多分量地震去噪方法研究及应用

针对多分量地震数据具有信噪比低、波场比较复杂等特点,同时为了解决有效波和面波在低频带具有相关性且频带有重叠的部分.本文基于希尔伯特—黄变换(Hilbert—Huang transform,HHT)提出一种多分量地震数据去噪方法.HHT的核心部分是固有模态分解(Empirical mode composition method,EMD),它将地震信号自适应分解为各个尺度的固有模态函数(Intrinsic mode function,IMF);然后根据有效波和噪声分布在不同的分辨尺度范围内进行数据重构;对于一些尺度分布比较宽的,进行更加严格的EMD分解,以及通过瞬时参数阈值大小进行处理以达到更好的去噪效果.本文所做的工作主要针对多波资料中面波的处理,利用面波跟有效波的频率尺度和瞬时属性差异,从仿真数据、模型和实际地震数据进行面波的压制实验,取得了比较满意的压制效果,提高了地震资料的分辨率.也表明了该方法在多分量地震去噪处理中的可行性.     

基于非二次幂Curvelet变换的最小二乘匹配算法及其应用

本文提出了基于非二次幂Curvelet变换的最小二乘匹配算法.首先,根据输入地震信号的频谱和方向等特征进行非二次幂Curvelet变换,根据其特征不同,最大程度地将有效信号和噪声分开;然后,在噪声能量集中的非二次幂Curvelet子记录上对输入数据和预测的噪声模型进行最小二乘匹配滤波处理.本方法提高了常规最小二乘匹配算法在时间空间域内进行信噪分离的稳定性和准确性.对含有面波的实际地震数据进行测试,其结果表明本方法可以有效地压制面波干扰,特别是当面波和有效信号有交叉或重叠等现象出现时,能较好地保护反射同相轴信息.本方法还可用于对含自由表面多次波和层间多次波等地震数据进行自适应信噪分离.     

Using the Wiener–Levinson algorithm to suppress ground-roll

In land seismic surveys, the seismic data are mostly contaminated by ground-roll noise, high amplitude and low frequency. Since the ground-roll is coherent with reflections and depends on the source, the spectral band of seismic signal and ground-roll always overlap, which can be clearly seen in the spectral domain. So, separating them in time or frequency domain commonly causes waveform distortions and information missing due to cut-off effects. Therefore, the combination of these factors leads to search for alternative filtering methods or processes. We applied the conventional Wiener–Levinson algorithm to extract ground-roll from the seismic data. Then, subtracting it from the seismic data arithmetically performs the ground-roll suppression. To set up the algorithm, linear or nonlinear sweep signals are used as reference noise trace. The frequencies needed in creating a reference noise trace using analytical sweep signal can be approximately estimated in spectral domain. The application of the proposed method based on redesigning of Wiener–Levinson algorithm differs from the usual frequency filtering techniques since the ground-roll is suppressed without cutting signal spectrum. The method is firstly tested on synthetics and then is applied to a shot data from the field. The result obtained from both synthetics and field data show that the ground-roll suppression in this way causes no waveform distortion and no reduction of frequency bandwidth of the data.     

基于地震子波支撑和Curvelet域稀疏约束的面波衰减方法

面波噪声衰减是地震数据处理流程中的重要一环,传统的面波衰减方法主要依靠面波与有效信号的几何特征差异,在变换域中将两者进行分离.受复杂近地表因素的影响,面波往往呈现非线性特征,并且在变换域中面波与有效信号存在部分重叠,这都导致面波噪声与有效信号难以彻底分离,消除面波的同时也损伤了有效信号.针对这一问题,本文综合利用Curvelet变换对地震数据的稀疏表征特性以及地震子波支撑来构建方程,通过Curvelet域稀疏约束来恢复压制面波时损失掉的有效信号.文中对该方法进行了模型试算和实际资料处理,处理结果表明:本文方法能够在一定程度上恢复损失的有效信号,提高了面波压制方法的保幅性.     

基于非稳态多项式拟合的地震噪声衰减方法研究(英文)

基于非稳态多项式拟合理论,针对地震数据中同相轴振幅变化这一特征,我们提出了一种地震噪声衰减的新方法。非稳态多项式拟合系数是时变的,通过整形正则化约束多项式拟和系数的光滑性,自适应的估计地震数据的相干分量。基于动校正后的共中心点道集(CMP)中地震信号的相干性,利用非稳态多项式拟合估计有效信号,从而衰减随机噪声。对于线性相干噪声,如地滚波,首先利用径向道变换(RadialTraceTransform,RTT)将地震数据变换到时间一视速度域,在时间—视速度域利用非稳态多项式拟合估计出相干噪声,然后减去相干噪声。该方法可以有效的估计振幅变化的相干分量,不需要相干分量振幅为常量的假设。模拟和实际资料处理结果表明,与传统的稳态多项式拟合和低切滤波相比,该方法可以更为有效的衰减地震噪声,同时保真了地震有效信号。     

基于地震信号波形形态差异的面波噪声稀疏优化分离方法

实际地震信号通常可表示为具有波形特征差异的多种基本波形信号的线性组合,如叠前道集中的工频干扰噪声与有效波信号、面波噪声与体波信号等.选择单一数学变换方法,往往不易实现地震信号的稀疏表示.近年来发展的形态成分分析理论,通过联合多种数学变换,可实现对复杂信号的稀疏表示.本文根据单道地震记录中面波与体波信号波形结构特征的差异性,提出一种基于形态成分分析的面波噪声衰减方法.针对面波的低频、窄带以及频散特性选择一维平稳小波变换作为其稀疏表示字典,而针对体波波形的局部相关特性选择局部离散余弦变换作为其稀疏表示字典,建立基于双波形字典的形态成分分析模型,通过求解该稀疏优化问题获得最终的信噪分离结果.理论模型和实际地震资料处理证实该方法不仅能够衰减单炮地震记录中的强面波干扰噪声,同时能够更好地保护有效信号的波形特征与频谱带宽,为地震资料的后续处理和分析提供良好的数据基础.     

径向道变换压制相干噪声方法研究

径向道变换(Radial Trace Transform)是将地震道集振幅值从偏移距一双程旅行时坐标系变换到视速度-双程旅行时坐标系,通过这种坐标系的变换,使相干噪声与有效信号在视速度和频率方面都有效分离.本文在介绍RT变换基本原理基础上,分析了RT变换中两种常用插值方法及其特点.并利用对模拟地震资料的处理,证明了RT域模拟-相减法较其带通滤波法在相干噪声压制与反射信号保持方面具有明显优势.最后,根据噪声特点,通过选择合理RT滤波参数,对实际地震资料进行处理试验,获得了较好的去噪效果,明显提高了资料信噪比,验证了研究方法的有效性.     

Ground-roll attenuation using a 2D time-derivative filter

We present a new filtering method for the attenuation of ground-roll. The method is based on the application of a bi-dimensional filter for obtaining the time-derivative of the seismograms. Before convolving the filter with the input data matrix, the normal moveout correction is applied to the seismograms with the purpose of flattening the reflections. The method can locally attenuate the amplitude of data of low frequency (in the ground-roll and stretch normal moveout region) and enhance flat events (reflections). The filtered seismograms can reveal horizontal or sub-horizontal reflections while vertical or sub-vertical events, associated with ground-roll, are attenuated. A regular set of samples around each neighbourhood data sample of the seismogram is used to estimate the time-derivative. A numerical approximation of the derivative is computed by taking the difference between the interpolated values calculated in both the positive and the negative neighbourhood of the desired position. The coefficients of the 2D time-derivative filter are obtained by taking the difference between two filters that interpolate at positive and negative times. Numerical results that use real seismic data show that the proposed method is effective and can reveal reflections masked by the ground-roll. Another benefit of the method is that the stretch mute, normally applied after the normal moveout correction, is unnecessary. The new filtering approach provides results of outstanding quality when compared to results obtained from the conventional FK filtering method.     

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

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

Restoration of missing offsets by parabolic Radon transform1

Restoration of missing offsets and trace interpolation is an interesting and important problem in seismic data processing. Based on the parabolic Radon transform, a method is presented for missing offset restoration, resampling and regularization of prestack individual CMP gathers. The method is also valid for resampling spatially aliased seismic data. The method is based on the parabolic assumption of the seismic events which is generally verified after a partial NMO correction in the CMP organization of the data. The essence of the method consists of a band-limited forward parabolic Radon transform of the data containing zero traces at the missing offset locations. The curvature range is chosen to map properly the coherent energy while the zero traces map beyond that range. After inverse transform the originally zero traces are partly filled with information. Several iterations of forward and inverse transform, every time replacing the zero traces in the original gather with the partially reconstructed ones, almost fully restore the zero traces. Efficient and fast algorithms can be built up to process data having a uniform geometry. Examples on synthetic as well as on field data demonstrate clearly the robustness of the method.     

基于混合时频分析技术的地震数据噪声压制(英文)

针对复杂地质结构、陡倾角相干噪声、空间采样不均匀等情况下F-x域反褶积去噪技术的不足,提出首先应用具有时-频聚集性度量准则的广义S变换将时间-空间域的地震数据变换至时间-频率-空间域(t-f-x)的数据,在t-f-x域中对每一个频率切片应用经验模态分解(EMD),移除噪声占主导地位的本征模态函数以压制相干和随机噪声的滤波方法。模型分析表明第一本征模态函数表征的高频信息以噪声为主,移除第一本证模态函数可以达到压制噪声的目的。经广义S变换后形成t-f-x域中EMD滤波方法等效于具有依赖于空间位置、频率、高波数截断特征的自适应f-k滤波。此滤波方法考虑了数据的局部时-频特征,且具有执行简单的特点。与AR预测滤波方法比较,此法滤除的成分包含较少的低波数的信息,滤除的成分非常的局部化,且获得结果没有表现出过度平滑的特征。实际资料的应用表明在经广义S变换后形成t-f-x域中运用EMD滤波方法能够有效地压制随机和陡倾角相干噪声。     

基于反射地震记录变子波模型提高地震记录分辨率

本文给出了地震记录变子波模型的一种近似数学表达式．基于该表达式研究了反射系数序列不满足白噪假设和子波在地下传播时发生变化这两种情况下地震道谱的组成及结构,讨论了谱白化及反褶积方法在这两种情况下效果不佳的原因．然后基于变子波模型,提出了一种新的提高地震记录分辨率的方法：第一步,用自适应于地震记录的Gabor分子窗把地震记录恰当地划分成若干片断,每段内信号近似平稳,然后将地震记录变换到时间-频率域;第二步,在变换域对每个分子窗内信号的振幅谱进行处理以拓宽频带;最后把处理后的时间-频率域函数反变换回时间域得到提高分辨率后的结果．本文提出的方法具有能较好地适用于反射系数不满足白噪假设的情况及提高分辨率后的地震记录能较好地保持原地震记录的相对能量关系等优点,模型和实际资料算例结果均表明,本文方法在拓宽地震资料频带及保持地震记录局部能量相对关系方面均明显优于谱白化方法．