scholarly journals Convolutional sparse coding for noise attenuation in seismic data

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. V23-V30
Author(s):  
Zhaolun Liu ◽  
Kai Lu

We have developed convolutional sparse coding (CSC) to attenuate noise in seismic data. CSC gives a data-driven set of basis functions whose coefficients form a sparse distribution. The noise attenuation method by CSC can be divided into the training and denoising phases. Seismic data with a relatively high signal-to-noise ratio are chosen for training to get the learned basis functions. Then, we use all (or a subset) of the basis functions to attenuate the random or coherent noise in the seismic data. Numerical experiments on synthetic data show that CSC can learn a set of shifted invariant filters, which can reduce the redundancy of learned filters in the traditional sparse-coding denoising method. CSC achieves good denoising performance when training with the noisy data and better performance when training on a similar but noiseless data set. The numerical results from the field data test indicate that CSC can effectively suppress seismic noise in complex field data. By excluding filters with coherent noise features, our method can further attenuate coherent noise and separate ground roll.

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. S65-S70 ◽  
Author(s):  
Lele Zhang ◽  
Evert Slob

Internal multiple reflections have been widely considered as coherent noise in measured seismic data, and many approaches have been developed for their attenuation. The Marchenko multiple elimination (MME) scheme eliminates internal multiple reflections without model information or adaptive subtraction. This scheme was originally derived from coupled Marchenko equations, but it was modified to make it model independent. It filters primary reflections with their two-way traveltimes and physical amplitudes from measured seismic data. The MME scheme is applied to a deepwater field data set from the Norwegian North Sea to evaluate its success in removing internal multiple reflections. The result indicates that most internal multiple reflections are successfully removed and primary reflections masked by overlapping internal multiple reflections are recovered.


Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1408-1416 ◽  
Author(s):  
David C. Henley

Coherent noise is a persistent problem in seismic imaging, and a number of techniques have been developed to attenuate it. The radial trace (RT) transform, a simple seismic data mapping algorithm, can be used as the basis for a particularly flexible and effective method for attenuating coherent noise on both prestack and poststack seismic data. Described here are the principles and some practical application details for attenuating coherent noise in the RT domain. A comparison between frequency–wavenumber (f–k) and RT domain filtering on a synthetic model is presented, and some of the differences and advantages of RT methods are identified. Next, RT coherent noise attenuation is demonstrated using a set of good‐quality field data; it is then applied to a very noisy data set. The results obtained with this last set prove to be as good as, or better than, those produced using f–k filtering.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB153-WB164 ◽  
Author(s):  
William Curry ◽  
Guojian Shan

Reflection seismic data typically are undersampled. Missing near offsets can be interpolated in reflection seismic data with pseudoprimaries, generated by crosscorrelating multiples and primaries in incomplete recorded data. These pseudoprimary data can be generated at the missing near offsets but contain many artifacts, so it is undesirable simply to replace the missing data with the pseudoprimaries. A nonstationary prediction-error filter (PEF) can instead be estimated from the pseudoprimaries and used to interpolate missing data to produce an interpolated output that is superior to direct substitution of the pseudoprimaries into the missing offsets. This approach is applied successfully to 2D synthetic and field data. Limitations in conventional acquisition geometry limit this approach in 3D, which can be illustrated using a synthetic data set.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. U67-U76 ◽  
Author(s):  
Robert J. Ferguson

The possibility of improving regularization/datuming of seismic data is investigated by treating wavefield extrapolation as an inversion problem. Weighted, damped least squares is then used to produce the regularized/datumed wavefield. Regularization/datuming is extremely costly because of computing the Hessian, so an efficient approximation is introduced. Approximation is achieved by computing a limited number of diagonals in the operators involved. Real and synthetic data examples demonstrate the utility of this approach. For synthetic data, regularization/datuming is demonstrated for large extrapolation distances using a highly irregular recording array. Without approximation, regularization/datuming returns a regularized wavefield with reduced operator artifacts when compared to a nonregularizing method such as generalized phase shift plus interpolation (PSPI). Approximate regularization/datuming returns a regularized wavefield for approximately two orders of magnitude less in cost; but it is dip limited, though in a controllable way, compared to the full method. The Foothills structural data set, a freely available data set from the Rocky Mountains of Canada, demonstrates application to real data. The data have highly irregular sampling along the shot coordinate, and they suffer from significant near-surface effects. Approximate regularization/datuming returns common receiver data that are superior in appearance compared to conventional datuming.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. C81-C92 ◽  
Author(s):  
Helene Hafslund Veire ◽  
Hilde Grude Borgos ◽  
Martin Landrø

Effects of pressure and fluid saturation can have the same degree of impact on seismic amplitudes and differential traveltimes in the reservoir interval; thus, they are often inseparable by analysis of a single stacked seismic data set. In such cases, time-lapse AVO analysis offers an opportunity to discriminate between the two effects. We quantify the uncertainty in estimations to utilize information about pressure- and saturation-related changes in reservoir modeling and simulation. One way of analyzing uncertainties is to formulate the problem in a Bayesian framework. Here, the solution of the problem will be represented by a probability density function (PDF), providing estimations of uncertainties as well as direct estimations of the properties. A stochastic model for estimation of pressure and saturation changes from time-lapse seismic AVO data is investigated within a Bayesian framework. Well-known rock physical relationships are used to set up a prior stochastic model. PP reflection coefficient differences are used to establish a likelihood model for linking reservoir variables and time-lapse seismic data. The methodology incorporates correlation between different variables of the model as well as spatial dependencies for each of the variables. In addition, information about possible bottlenecks causing large uncertainties in the estimations can be identified through sensitivity analysis of the system. The method has been tested on 1D synthetic data and on field time-lapse seismic AVO data from the Gullfaks Field in the North Sea.


2010 ◽  
Vol 14 (3) ◽  
pp. 545-556 ◽  
Author(s):  
J. Rings ◽  
J. A. Huisman ◽  
H. Vereecken

Abstract. Coupled hydrogeophysical methods infer hydrological and petrophysical parameters directly from geophysical measurements. Widespread methods do not explicitly recognize uncertainty in parameter estimates. Therefore, we apply a sequential Bayesian framework that provides updates of state, parameters and their uncertainty whenever measurements become available. We have coupled a hydrological and an electrical resistivity tomography (ERT) forward code in a particle filtering framework. First, we analyze a synthetic data set of lysimeter infiltration monitored with ERT. In a second step, we apply the approach to field data measured during an infiltration event on a full-scale dike model. For the synthetic data, the water content distribution and the hydraulic conductivity are accurately estimated after a few time steps. For the field data, hydraulic parameters are successfully estimated from water content measurements made with spatial time domain reflectometry and ERT, and the development of their posterior distributions is shown.


Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. R199-R217 ◽  
Author(s):  
Xintao Chai ◽  
Shangxu Wang ◽  
Genyang Tang

Seismic data are nonstationary due to subsurface anelastic attenuation and dispersion effects. These effects, also referred to as the earth’s [Formula: see text]-filtering effects, can diminish seismic resolution. We previously developed a method of nonstationary sparse reflectivity inversion (NSRI) for resolution enhancement, which avoids the intrinsic instability associated with inverse [Formula: see text] filtering and generates superior [Formula: see text] compensation results. Applying NSRI to data sets that contain multiples (addressing surface-related multiples only) requires a demultiple preprocessing step because NSRI cannot distinguish primaries from multiples and will treat them as interference convolved with incorrect [Formula: see text] values. However, multiples contain information about subsurface properties. To use information carried by multiples, with the feedback model and NSRI theory, we adapt NSRI to the context of nonstationary seismic data with surface-related multiples. Consequently, not only are the benefits of NSRI (e.g., circumventing the intrinsic instability associated with inverse [Formula: see text] filtering) extended, but also multiples are considered. Our method is limited to be a 1D implementation. Theoretical and numerical analyses verify that given a wavelet, the input [Formula: see text] values primarily affect the inverted reflectivities and exert little effect on the estimated multiples; i.e., multiple estimation need not consider [Formula: see text] filtering effects explicitly. However, there are benefits for NSRI considering multiples. The periodicity and amplitude of the multiples imply the position of the reflectivities and amplitude of the wavelet. Multiples assist in overcoming scaling and shifting ambiguities of conventional problems in which multiples are not considered. Experiments using a 1D algorithm on a synthetic data set, the publicly available Pluto 1.5 data set, and a marine data set support the aforementioned findings and reveal the stability, capabilities, and limitations of the proposed method.


Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 191-203 ◽  
Author(s):  
A. Frank Linville ◽  
Robert A. Meek

Primary reflections in seismic records are often obscured by coherent noise making processing and interpretation difficult. Trapped water modes, surface waves, scattered waves, air waves, and tube waves to name a few, must be removed early in the processing sequence to optimize subsequent processing and imaging. We have developed a noise canceling algorithm that effectively removes many of the commonly encountered noise trains in seismic data. All currently available techniques for coherent noise attenuation suffer from limitations that introduce unacceptable signal distortions and artifacts. Also, most of those techniques impose the dual stringent requirements of equal and fine spatial sampling in the field acquisition of seismic data. Our technique takes advantage of characteristics usually found in coherent noise such as being localized in time, highly aliased, nondispersive (or only mildly so), and exhibit a variety of moveout patterns across the seismic records. When coherent noise is localized in time, a window much like a surgical mute is drawn around the noise. The algorithm derives an estimate of the noise in the window, automatically correcting for amplitude and phase differences, and adaptively subtracts this noise from the window of data. This signal estimate is then placed back in the record. In a model and a land data example, the algorithm removes noise more effectively with less signal distortion than does f-k filtering or velocity notch filtering. Downgoing energy in a vertical seismic profile (VSP) with irregular receiver spacing is also removed.


Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. V1-V10
Author(s):  
Julián L. Gómez ◽  
Danilo R. Velis ◽  
Juan I. Sabbione

We have developed an empirical-mode decomposition (EMD) algorithm for effective suppression of random and coherent noise in 2D and 3D seismic amplitude data. Unlike other EMD-based methods for seismic data processing, our approach does not involve the time direction in the computation of the signal envelopes needed for the iterative sifting process. Instead, we apply the sifting algorithm spatially in the inline-crossline plane. At each time slice, we calculate the upper and lower signal envelopes by means of a filter whose length is adapted dynamically at each sifting iteration according to the spatial distribution of the extrema. The denoising of a 3D volume is achieved by removing the most oscillating modes of each time slice from the noisy data. We determine the performance of the algorithm by using three public-domain poststack field data sets: one 2D line of the well-known Alaska 2D data set, available from the US Geological Survey; a subset of the Penobscot 3D volume acquired offshore by the Nova Scotia Department of Energy, Canada; and a subset of the Stratton 3D land data from South Texas, available from the Bureau of Economic Geology at the University of Texas at Austin. The results indicate that random and coherent noise, such as footprint signatures, can be mitigated satisfactorily, enhancing the reflectors with negligible signal leakage in most cases. Our method, called empirical-mode filtering (EMF), yields improved results compared to other 2D and 3D techniques, such as [Formula: see text] EMD filter, [Formula: see text] deconvolution, and [Formula: see text]-[Formula: see text]-[Formula: see text] adaptive prediction filtering. EMF exploits the flexibility of EMD on seismic data and is presented as an efficient and easy-to-apply alternative for denoising seismic data with mild to moderate structural complexity.


Sign in / Sign up

Export Citation Format

Share Document