首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 375 毫秒
1.
地震信号随机噪声压制的双树复小波域双变量方法   总被引:2,自引:2,他引:0       下载免费PDF全文
有效地压制地震信号中的噪声是地震信号解释和后续处理的重要环节之一.本文建立两种双树复小波域双变量模型对地震信号中的随机噪声进行压制.地震信号经双树复小波变换后,同一方向实部与虚部系数、实部(或虚部)系数与对应的模之间存在较强的相关性.鉴于此,对同一方向实部与虚部小波系数建立双变量模型,从含噪地震信号小波系数中估计原始信号的小波系数,再基于双树复小波逆变换重构得到降噪后的地震信号.进一步对同一方向实部(或虚部)系数与对应的模建立双变量模型,得到地震信号随机噪声压制的第二种双树复小波域双变量方法.最后对合成地震记录和实际地震资料中的随机噪声进行压制的实验结果证实本文两种方法都能够有效地压制地震信号中的随机噪声.  相似文献   

2.
估计地震数据的信噪比对于地震数据的处理和解释具有重要作用.以往估计地震数据信噪比的方法都需要分离数据中的有效信号和噪声,然后再估计相应的信噪比.这些估计方法的精度严重依赖信号估计方法或噪声压制方法的有效性,往往存在偏差.本文提出一种估计地震数据局部信噪比的深度卷积神经网络模型,通过迭代训练优化参数,构建从含噪地震数据到其信噪比的特征映射.然后使用该神经网络完成信噪比的推理预测,不需要分离地震数据中的有效信号和噪声.模拟数据和实际资料的处理结果都表明,本文的方法可以准确而高效地估计局部地震数据的信噪比,为地震数据质量的定量评价提供依据.  相似文献   

3.
地震信号中的随机噪声是一种干扰波,严重降低了地震信号的信噪比,并影响着资料的后续处理和分析.本文根据地震信号中有效信号和随机噪声的差异,结合分数阶B样条小波变换与高斯尺度混合模型提出了一种地震信号随机噪声压制方法.首先利用分数阶B样条小波变换将含噪地震信号映射到最优分数阶小波时频域内,然后对各小波子带系数分别建立高斯尺度混合模型,由贝叶斯方法估计出源地震信号小波系数,最后使用分数阶B样条小波逆变换重构得到降噪后的地震信号.利用本文方法对合成地震记录和实际地震信号进行降噪处理,实验结果表明本文方法能够有效地压制地震信号中的随机噪声,并且较好地保留了有效信号.  相似文献   

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

5.
基于奇异谱分析的联合去噪及规则化方法   总被引:1,自引:0,他引:1  
在地震数据处理中,噪音压制和数据规则化研究直接影响后期的地震处理及解释效果.本文提出一种基于奇异谱分析的联合去噪及规则化方法,在迭代时自动在地震道缺失的位置进行插值,在含有地震数据的位置压制噪声,并通过分步插值改善数据相干性.通过合成地震记录和实际资料的联合去噪及规则化处理结果表明:联合方法能够在补全地震道的同时有效地压制噪声.  相似文献   

6.
GNMF小波谱分离在地震勘探噪声压制中的应用   总被引:2,自引:1,他引:1       下载免费PDF全文
田雅男  李月  林红波  吴宁 《地球物理学报》2015,58(12):4568-4575
地震勘探资料噪声压制及信噪比提高是整个地震勘探信号处理过程中的重要任务,随着地震勘探深度的增加及其复杂性,人们对地震数据质量的要求越来越高.勘探环境的复杂化使得采集到的地震资料中有效信号被大量噪声淹没,无法清晰辨识,严重影响后续的数据处理与解释.小波去噪是地震勘探中常用且发展较成熟的一种方法,但是其涉及到的阈值函数选取问题一直令人困扰,虽然已有多种阈值函数被提出,但仍存在各自的缺陷.本文利用小波分解在时域及频域良好的信号细节体现特性,引入模式识别中的非负矩阵分解(NMF)谱分离思想,针对小波系数阈值优化问题,提出了一种小波域图非负矩阵分解(GNMF)消噪算法.该方法首先在小波分解基础上,利用GNMF算法实现小波分解系数谱中信号分量与噪声分量的谱分离,然后通过反变换重构各分离子谱对应的子信号,最后利用K均值聚类算法将得到的多个子信号划分为信号类及噪声类,最终得到重构信号及分离噪声.合成记录和实际地震资料的消噪结果验证了新方法在提高信号与噪声分离准确性和精度方面的有效性,同时新方法避免了阈值选取造成的噪声压制不理想或有效成分损失问题.与小波消噪结果的对比及数值分析也说明了新方法在噪声压制及有效成分保持方面的优势.  相似文献   

7.
时频峰值滤波去噪技术及其应用   总被引:3,自引:0,他引:3       下载免费PDF全文
本文将时频峰值滤波(TFPF)去噪技术应用于共炮点地震资料的随机噪声压制.时频峰值滤波技术是通过频率调制将信号调制成解析信号的瞬时频率,利用解析信号的Wigner-Ville分布的峰值进行瞬时频率估计,恢复有效信号,与其它去噪方法相比,TFPF具有在较少的约束条件下压制强随机噪声的优点.本文针对实际地震资料的非线性特性,利用加窗的Wigner-Ville分布实现TFPF,使得地震信号在一个窗长内近似满足线性瞬时频率条件,减小由地震信号非线性引起的偏差.本文对共炮点地震记录做时频峰值滤波处理,滤波结果表明在地震勘探资料中存在强随机噪声的情况下,利用局部线性化处理的时频峰值滤波技术可以有效地压制地震资料中的随机噪声,恢复出湮没在随机噪声中的地震反射信号.信噪比提高3~6 dB.  相似文献   

8.
海上地震勘探中为提高采集效率会采取连续采集方案,导致地震记录尾端与下一炮起始记录重合,这种现象对浅层地震处理影响不大,但是会完全掩盖深层地震信号.由于混叠噪声与真实信号属于同方向及同震源采集,其频谱完全重叠,常规的噪声压制方法无法完全压制混叠噪声.本文提出了一种基于同相轴追踪的混叠记录分离方法,利用时域上混叠噪声与真实地震信号时距曲线的曲率差异得到混叠噪声的模型,再利用最小平方约束反演方法进行自适应相减,使有效信号误差最小,最终完成混叠炮集分离.通过理论分析与正演测试,混叠衰减的关键在于求取可靠的混叠模型,而自适应衰减则控制细微振幅与相位差异.通过对实际地震资料进行处理,并与F-K滤波方法进行了对比,证实了该方法的有效性.  相似文献   

9.
常规的时间-空间域和频率-空间域预测滤波方法假设地震记录由地震信号和随机噪声两部分构成,即所谓的加噪声模型,但是,在对随机噪声进行估算时,又假设随机噪声可以通过预测误差滤波器由地震记录中进行预测,即所谓的源噪声模型。这种前后不一致的噪声模型降低了该类方法的去噪能力和保幅性能。为此,本文提出了一种基于反演的时空域随机噪声衰减方法。它首先从地震数据中估算预测滤波算子,该算子表征了地震信号的可预测性,自适应地描述了地震信号的空间结构。在得到预测误差算子之后,将该算子作为正则化约束引入到地震信号反演系统,由含有随机噪声的地震数据直接反演地震信号。不同于常规随机噪声衰减方法,该方法将随机噪声衰减问题归结为正则化约束下的地震信号反演问题,克服了常规方法噪声模型的不一致性问题。我们采用模型数据和实际数据进行了实验分析,并与常规方法进行了效果对比。实验结果表明:与常规方法相比,本文方法在噪声压制的同时,没有对有效信号产生明显伤害,具有更好的振幅保持能力。  相似文献   

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

11.
This article utilizes Savitzky–Golay (SG) filter to eliminate seismic random noise. This is a novel method for seismic random noise reduction in which SG filter adopts piecewise weighted polynomial via leastsquares estimation. Therefore, effective smoothing is achieved in extracting the original signal from noise environment while retaining the shape of the signal as close as possible to the original one. Although there are lots of classical methods such as Wiener filtering and wavelet denoising applied to eliminate seismic random noise, the SG filter outperforms them in approximating the true signal. SG filter will obtain a good tradeoff in waveform smoothing and valid signal preservation under suitable conditions. These are the appropriate window size and the polynomial degree. Through examples from synthetic seismic signals and field seismic data, we demonstrate the good performance of SG filter by comparing it with the Wiener filtering and wavelet denoising methods.  相似文献   

12.
Conventional time-space domain and frequency-space domain prediction filtering methods assume that seismic data consists of two parts, signal and random noise. That is, the so-called additive noise model. However, when estimating random noise, it is assumed that random noise can be predicted from the seismic data by convolving with a prediction error filter. That is, the source-noise model. Model inconsistencies, before and after denoising, compromise the noise attenuation and signal-preservation performances of prediction filtering methods. Therefore, this study presents an inversion-based time-space domain random noise attenuation method to overcome the model inconsistencies. In this method, a prediction error filter (PEF), is first estimated from seismic data; the filter characterizes the predictability of the seismic data and adaptively describes the seismic data’s space structure. After calculating PEF, it can be applied as a regularized constraint in the inversion process for seismic signal from noisy data. Unlike conventional random noise attenuation methods, the proposed method solves a seismic data inversion problem using regularization constraint; this overcomes the model inconsistency of the prediction filtering method. The proposed method was tested on both synthetic and real seismic data, and results from the prediction filtering method and the proposed method are compared. The testing demonstrated that the proposed method suppresses noise effectively and provides better signal-preservation performance.  相似文献   

13.
随着数字化地震计的应用和发展,数字观测资料中不仅记录了宽频带地震信息,也包含了大量的干扰信号。在实际工作中,对地震信息的分析处理影响很大,有时还会造成地震波形记录不完整或地震信号淹没在干扰信号中,很难再进行深入的分析和应用。因此,利用MATLAB软件设计IIR数字滤波器对不同记录中的干扰波进行滤波器设计,达到消除干扰波的目的。实例表明,此方法在台站应用效果良好。  相似文献   

14.
陈天  易远元 《地震学报》2021,43(4):474-482
本文以提高地震数据的成像质量为目标,提出一种智能的卷积神经网络降噪框架,从带有噪声的地震数据中自适应地学习地震信号。为了加速网络训练和避免训练时出现梯度消失现象,我们在网络中加入残差学习和批标准化的方法,并采用了ReLU激活函数和Adam优化算法优化网络。此外,Marmousi和F3数据集被用来对网络进行训练和测试,经过充分训练的网络不仅能在学习中保留地震数据特征,而且能去除随机噪声。首先充分地训练网络,从中提取出随机噪声,并保留学习到的地震数据特征,之后通过重建地震数据估算测试集中的波形特征。合成记录和实际数据的处理结果显示了深度卷积神经网络在随机噪声压制任务中的潜力,并通过实验验证表明了深度卷积神经网络框架有很好的去噪效果。   相似文献   

15.
空间光滑地震活动性模型中光滑函数的比较研究   总被引:2,自引:1,他引:1       下载免费PDF全文
徐伟进  高孟潭 《地震学报》2012,34(2):244-256
使用Frankel提出的基于空间光滑地震活动性模型的地震危险性分析方法,选择华南、华北、川滇3个地区的地震记录,比较分析了高斯、幂律和地震分形分布光滑函数3种光滑函数在不同地区的适用性.结果表明,使用交叉验证法可以为高斯光滑函数选取合适的相关距离c值,光滑得到的地震活动性模型能够真实反映研究区域的地震活动特征,根据活动性模型计算得出的峰值加速度(PGA)分布也符合人们对研究区域地震危险性的认识.幂律光滑函数适用于地震活动性较强的地区,且具有容易求取光滑参数的优点.光滑程度较低的幂律光滑函数不适用于地震活动性弱的地区,在该类地区应选择光滑程度较高的高斯光滑函数.地震分形分布光滑函数不适用于地震活动较强且地震活动强度差异较大的地区,其容易过分高估高震级地震对地震危险性的影响,而忽略了低震级地震对地震危险性的贡献.但对于地震活动较弱且地震活动强度差异较小的地区,可使用地震分形分布光滑函数,且同样具有容易求取光滑参数的优点.   相似文献   

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

17.
An accurate estimate of the seismic wavelet on a seismic section is extremely important for interpretation of fine details on the section and for estimation of acoustic impedance. In the absence of well-control, the recognized best approach to wavelet estimation is to use the technique of multiple coherence analysis to estimate the coherent signal and its amplitude spectrum, and thence construct the seismic wavelet under the minimum-phase assumption. The construction of the minimum-phase wavelet is critically dependent on the decay of the spectrum at the low-frequency end. Traditional methods of cross-spectral estimation, such as frequency smoothing using a Papoulis window, suffer from substantial side-lobe leakage in the areas of the spectrum where there is a large change of power over a relatively small frequency range. The low-frequency end of the seismic spectrum (less than 4 Hz) decays rapidly to zero. Side-lobe leakage causes poor estimates of the low-frequency decay, resulting in degraded wavelet estimates. Thomson's multitaper method of cross-spectral estimation which suffers little from side-lobe leakage is applied here, and compared with the result of using frequency smoothing with the Papoulis window. The multitaper method seems much less prone to estimating spuriously high coherences at very low frequencies. The wavelet estimated by the multitaper approach from the data used here is equivalent to imposing a low-frequency roll-off of some 48 dB/oct (below 3.91 Hz) on the amplitude spectrum. Using Papoulis smoothing the equivalent roll-off is only about 36 dB/oct. Thus the multitaper method gives a low-frequency decay rate of the amplitude spectrum which is some 4 times greater than for Papoulis smoothing. It also gives more consistent results across the section. Furthermore, the wavelet obtained using the multi-taper method and seismic data only (with no reference to well data) has more attractive physical characteristics when compared with a wavelet extracted using well data, than does an estimate using traditional smoothing.  相似文献   

18.
A new seismic interpolation and denoising method with a curvelet transform matching filter, employing the fast iterative shrinkage thresholding algorithm (FISTA), is proposed. The approach treats the matching filter, seismic interpolation, and denoising all as the same inverse problem using an inversion iteration algorithm. The curvelet transform has a high sparseness and is useful for separating signal from noise, meaning that it can accurately solve the matching problem using FISTA. When applying the new method to a synthetic noisy data sets and a data sets with missing traces, the optimum matching result is obtained, noise is greatly suppressed, missing seismic data are filled by interpolation, and the waveform is highly consistent. We then verified the method by applying it to real data, yielding satisfactory results. The results show that the method can reconstruct missing traces in the case of low SNR (signal-to-noise ratio). The above three problems can be simultaneously solved via FISTA algorithm, and it will not only increase the processing efficiency but also improve SNR of the seismic data.  相似文献   

19.
Optimum pilot sweep   总被引:3,自引:0,他引:3  
The successful application of high-resolution seismic methods requires evaluating each element in the seismic system and ensuring that each part of the system contributes optimally to the success of the method. Unfortunately, unlike data processing, seismic signal generation is not carefully optimized. The purpose of our study was to optimize the source signal in order to better coordinate field operations with subsequent data processing to achieve their common objective. We developed an iterative method for a rational frequency distribution of the energy of a seismic source. The method allows the optimum amplitude spectrum of a source signal to be calculated, thus providing the best data quality at the end of the processing. We assume that the source signal is affected by a total transfer function, by the reflectivity function of a target interval, and by ambient noise, whose characteristics, if not known, can be estimated or measured in practice. The transfer function includes data processing other than the correlation stage and the final trace-optimizing filter. The variance of a reflectivity estimate is considered to be a measure of the data quality and improvement of the characteristic corresponds to a decrease in the variance. For this reason, a constrained Wiener deconvolution filter is used as the final trace-optimizing filter. It not only minimizes the variance of a reflectivity estimate but also ensures a specific signal-to-noise ratio. The method is made feasible by following the Vibroseis technique, primarily because of the versatility of the technique in controlling the signal spectrum. With the optimum amplitude spectrum obtained, the corresponding optimum pilot sweep can be readily calculated. Examples using synthetic data are presented to illustrate the method.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号