Multiple attenuation in complex geology with a pattern-based approach

Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. V97-V107 ◽  
Author(s):  
Antoine Guitton
Keyword(s):  
Seismic Data ◽  
Prediction Error ◽  
Input Data ◽  
Multiple Model ◽  
Ideal Case ◽  
Multiple Attenuation ◽  
Time Space ◽  
Multiple Prediction ◽  
Complex Geology ◽  
Space Variant

Primaries (signal) and multiples (noise) often exhibit different kinematics and amplitudes (i.e., patterns) in time and space. Multidimensional prediction-error filters (PEFs) approximate these patterns to separate noise and signal in a least-squares sense. These filters are time-space variant to handle the nonstationarity of multioffset seismic data. PEFs for the primaries and multiples are estimated from pattern models. In an ideal case where accurate pattern models of both noise and signal exist, the pattern-based method recovers the primaries while preserving their amplitudes. In the more general case, the pattern model of the multiples is obtained by using the data as prediction operators. The pattern model of the primaries is obtained by convolving the noise PEFs with the input data. In this situation, 3D PEFs are preferred to separate (in prestack data) the multiples properly and to preserve the primaries. Comparisons of the proposed method with adaptive subtraction with an [Formula: see text] norm demonstrate that for a given multiple model, the pattern-based approach generally attenuates the multiples and recovers the primaries better. In addition, tests on a 2D line from the Gulf of Mexico demonstrate that the proposed technique copes fairly well with modeling inadequacies present in the multiple prediction.

scholarly journals Multiple prediction through inversion: A fully data‐driven concept for surface‐related multiple attenuation

Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 547-553 ◽  
Author(s):  
Yanghua Wang
Keyword(s):  
Seismic Data ◽  
Model Building ◽  
Data Driven ◽  
Multiple Model ◽  
Data Set ◽  
Multiple Attenuation ◽  
Multiple Prediction ◽  
Surface Operator ◽  
Previous Iteration ◽  
Conventional Scheme

This paper introduces a fully data‐driven concept, multiple prediction through inversion (MPI), for surface‐related multiple attenuation (SMA). It builds the multiple model not by spatial convolution, as in a conventional SMA, but by updating the attenuated multiple wavefield in the previous iteration to generate a multiple prediction for the new iteration, as is usually the case in an iterative inverse problem. Because MPI does not use spatial convolution, it is able to minimize the edge effect that appears in conventional SMA multiple prediction and to eliminate the need to synthesize near‐offset traces, required by a conventional scheme, so that it can deal with a seismic data set with missing near‐offset traces. The MPI concept also eliminates the need for an explicit surface operator, which is required by conventional SMA and is comprised of the inverse source signature and other effects. This method accounts implicitly for the spatial variation of the surface operator in multiple‐model building and attempts to predict multiples which are not only accurate kinematically but are also accurate in phase and amplitude.

scholarly journals Efficient 2D multiple attenuation using SRME with curvelet-domain subtraction

2022 ◽  
Vol 43 (1) ◽  
Author(s):  
Szu-Ying Lai ◽  
Yunung Nina Lin ◽  
Ho-Han Hsu
Keyword(s):  
Finite Number ◽  
Least Squares ◽  
Seismic Data ◽  
Parameter Determination ◽  
Prediction Errors ◽  
Geological Condition ◽  
Multiple Model ◽  
Multiple Attenuation ◽  
Time Space ◽  
Multiple Elimination

AbstractSurface Related Multiple Elimination (SRME) usually suffers the issue of either over-attenuation that damages the primaries or under-attenuation that leaves strong residual multiples. This dilemma happens commonly when SRME is combined with least-squares subtraction. Here we introduce a more sophisticated subtraction approach that facilitates better separation of multiples from primaries. Curvelet-domain subtraction transforms both the data and the multiple model into the curvelet domain, where different frequency bands (scales) and event directions (orientations) are represented by a finite number of curvelet coefficients. When combined with adaptive subtraction in the time–space domain, this method can handle model prediction errors to achieve effective subtraction. We demonstrate this method on two 2D surveys from the TAiwan Integrated GEodynamics Research (TAIGER) project. With a careful parameter determination flow, our result shows curvelet-domain subtraction outperforms least-squares subtraction in all geological settings. We also present one failed case where specific geological condition hinders proper multiple subtraction. We further demonstrate that even for data acquired with short cables, curvelet-domain subtraction can still provide better results than least-squares subtraction. We recommend this method as the standard processing flow for multi-channel seismic data.

scholarly journals Improving adaptive subtraction in seismic multiple attenuation

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. V59-V67 ◽  
Author(s):  
Shoudong Huo ◽  
Yanghua Wang
Keyword(s):  
Least Squares ◽  
Seismic Data ◽  
Primary Energy ◽  
Multiple Models ◽  
Multiple Model ◽  
Data Set ◽  
Multiple Attenuation ◽  
Matching Filter ◽  
Theoretical Analyses ◽  
Different Order

In seismic multiple attenuation, once the multiple models have been built, the effectiveness of the processing depends on the subtraction step. Usually the primary energy is partially attenuated during the adaptive subtraction if an [Formula: see text]-norm matching filter is used to solve a least-squares problem. The expanded multichannel matching (EMCM) filter generally is effective, but conservative parameters adopted to preserve the primary could lead to some remaining multiples. We have managed to improve the multiple attenuation result through an iterative application of the EMCM filter to accumulate the effect of subtraction. A Butterworth-type masking filter based on the multiple model can be used to preserve most of the primary energy prior to subtraction, and then subtraction can be performed on the remaining part to better suppress the multiples without affecting the primaries. Meanwhile, subtraction can be performed according to the orders of the multiples, as a single subtraction window usually covers different-order multiples with different amplitudes. Theoretical analyses, and synthetic and real seismic data set demonstrations, proved that a combination of these three strategies is effective in improving the adaptive subtraction during seismic multiple attenuation.

Surface‐related multiple elimination on land seismic data—Strategies via case studies

Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 719-734 ◽  
Author(s):  
Panos G. Kelamis ◽  
D. J. Verschuur
Keyword(s):  
Case Studies ◽  
Seismic Data ◽  
Input Data ◽  
Multiple Attenuation ◽  
Processing Strategies ◽  
Processing Step ◽  
Subsequent Elimination ◽  
Seismic Amplitudes ◽  
Multiple Elimination

Three processing strategies for the estimation and subsequent elimination of surface‐related multiple energy on land seismic data are presented. They can be applied in a prestack mode (to shot and common‐midpoint gathers) or in a poststack mode. The algorithm for the multiple attenuation is based on wave theoretical principles in which the data are used as a prediction operator. The estimated multiples are then adaptively subtracted from the input data to obtain primary‐only data. A processing step prior to applying multiple elimination is an important component of these methodologies, particularly in the prestack analysis. Its aim is to regularize the data, improve the S/N ratio, and balance the seismic amplitudes. This results in smooth prediction operators. The effectiveness of these schemes in suppressing multiples is demonstrated with a number of case studies involving processing land seismic data.

scholarly journals Practical Implementation of Multiple Attenuation Methods on 2D Deepwater Seismic Data : Seram Sea Case Study

2017 ◽  
Vol 32 (1) ◽  
Author(s):  
Tumpal Bernhard Nainggolan ◽  
Deny Setiady
Keyword(s):  
Seismic Data ◽  
Input Data ◽  
Complex Structure ◽  
Least Square ◽  
Inversion Method ◽  
Practical Implementation ◽  
Proper Solution ◽  
Multiple Attenuation ◽  
Multiple Elimination ◽  
Water Bottom

Some deepwater multiple attenuation processing methods have been developed in the past with partial success. The success of surface multiple attenuation relies on good water bottom reflections for most deepwater marine situations. It brings the bigger ability to build an accurate water bottom multiple prediction model. Major challenges on 2D deepwater seismic data processing especially such a geologically complex structure of Seram Sea, West Papua – Indonesia are to attenuate surface related multiple and to preserve the primary data. Many multiple attenuation methods have been developed to remove surface multiple on these seismic data including most common least-squares, prediction-error filtering and more advanced Radon transform.Predictive Deconvolution and Surface Related Multiple Elimination (SRME) method appears to be a proper solution, especially in complex structure where the above methods fail to distinguish interval velocity difference between primaries and multiples. It does not require any subsurface info as long as source signature and surface reflectivity are provided. SRME method consists of 3 major steps: SRME regularization, multiple modeling and least-square adaptive subtraction. Near offset regularization is needed to fill the gaps on near offset due to unrecorded near traces during the acquisition process. Then, isolating primaries from multiples using forward modeling. Inversion method by subtraction of input data with multiple models to a more attenuated multiple seismic section.Results on real 2D deepwater seismic data show that SRME method as the proper solution should be considered as one of the practical implementation steps in geologically complex structure and to give more accurate seismic imaging for the interpretation.Keywords : multiple attenuation, 2D deepwater seismic, Radon transform, Surface Related Multiple Elimination (SRME). Banyak metode atenuasi pengulangan ganda dikembangkan pada pengolahan data seismik dengan tingkat keberhasilan yang rendah pada masa lalu. Keberhasilan dalam atenuasi pengulangan ganda permukaan salah satunya bergantung pada hasil gelombang pantul pada batas dasar laut dan permukaan pada hampir seluruh survei seismik laut. Hal tersebut menentukan keakuratan dalam membuat model prediksi pengulangan ganda dasar laut dan permukaan air. Tantangan utama dalam pemrosesan data seismik 2D laut dalam khususnya struktur geologi kompleks seperti Laut Seram, Papua Barat – Indonesia adalah pada kegiatan menekan pengulangan ganda permukaan sekaligus mempertahankan data primer. Beberapa metode yang dikembangkan untuk menghilangkan pengulangan ganda permukaan pada data seismik seperti least-square, filter prediksi kesalahan dan transformasi Radon.Dekonvolusi Prediktif dan Metode Surface Related Multiple Elimination (SRME) digunakan sebagai solusi yang baik pada struktur kompleks dimana metode-metode lain gagal untuk memisahkan perbedaan kecepatan interval data primer dan pengulangan ganda. Metode tersebut tidak membutuhkan informasi bawah permukaan selain parameter sumber dan reflektivitas permukaan. Metode SRME terdiri dari 3 tahapan utama : regularisasi SRME, pemodelan pengulangan ganda dan pengurangan adaktif least-square. Regularisasi near offset diperlukan untuk mengisi kekosongan pada near offset yang disebabkan oleh adanya sejumlah tras terdekat yang tidak terekam selama akuisisi. Pemodelan maju digunakan untuk memisahkan data primer dan pengulangan ganda kemudian inversi dengan pengurangan input data dengan model multiple.Hasil pada data seismik 2D laut dalam menunjukkan bahwa metode SRME layak diterapkan sebagai salah satu pengembangan metode atenuasi multiple permukaan serta untuk meningkatkan akurasi data seismik terutama untuk struktur geologi kompleks.Kata kunci : peredaman pengulangan ganda (multiple), seismik 2D laut dalam, transformasi Radon, Surface Related Multiple Attenuation (SRME).

Multiple attenuation with 3D high-order high-resolution parabolic Radon transform using lower frequency constraints

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. V317-V328
Author(s):  
Jitao Ma ◽  
Guoyang Xu ◽  
Xiaohong Chen ◽  
Xiaoliu Wang ◽  
Zhenjiang Hao
Keyword(s):  
High Resolution ◽  
Radon Transform ◽  
Seismic Data ◽  
Computational Cost ◽  
Real Data ◽  
Multiple Model ◽  
Multiple Attenuation ◽  
Common Depth Point ◽  
Frequency Constraint ◽  
Parabolic Radon Transform

The parabolic Radon transform is one of the most commonly used multiple attenuation methods in seismic data processing. The 2D Radon transform cannot consider the azimuth effect on seismic data when processing 3D common-depth point gathers; hence, the result of applying this transform is unreliable. Therefore, the 3D Radon transform should be applied. The theory of the 3D Radon transform is first introduced. To address sparse sampling in the crossline direction, a lower frequency constraint is introduced to reduce spatial aliasing and improve the resolution of the Radon transform. An orthogonal polynomial transform, which can fit the amplitude variations in different parabolic directions, is combined with the dealiased 3D high-resolution Radon transform to account for the amplitude variations with offset of seismic data. A multiple model can be estimated with superior accuracy, and improved results can be achieved. Synthetic and real data examples indicate that even though our method comes at a higher computational cost than existing techniques, the developed approach provides better attenuation of multiples for 3D seismic data with amplitude variations.

scholarly journals Multiple subtraction using an expanded multichannel matching filter

Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 346-354 ◽  
Author(s):  
Yanghua Wang
Keyword(s):  
Hilbert Transform ◽  
Spatial Dimension ◽  
Multiple Model ◽  
Space Domain ◽  
First Derivative ◽  
Multiple Attenuation ◽  
Time Space ◽  
Multichannel Filter ◽  
Matching Filter ◽  
The Hilbert Transform

An expanded multichannel matching (EMCM) filter is proposed for the adaptive subtraction in seismic multiple attenuation. For a normal multichannel matching filter where an original seismic trace is matched by a group of multiple‐model traces, the lateral coherency of adjacent traces is likely to be exploited to discriminate the overlapped multiple and primary reflections. In the proposed EMCM filter, a seismic trace is matched by not only a group of the ordinary multiple‐model traces but also their adjoints generated mathematically. The adjoints of a multiple‐model trace include its first derivative, its Hilbert transform, and the derivative of the Hilbert transform. The convolutional coefficients associated with the normal multichannel filter can be represented as a 2D operator in the time‐space domain. This 2D operator is expanded with an additional spatial dimension in the EMCM filter to improve the robustness of the adaptive subtraction. The multiple‐model traces are generated using moveout equations to afford efficiency in the multiple attenuation application.

5D de-aliased seismic data interpolation using non-stationary prediction error filter

Geophysics ◽  
2021 ◽  
pp. 1-92
Author(s):  
Yangkang Chen ◽  
Sergey Fomel ◽  
Hang Wang ◽  
shaohuan zu
Keyword(s):  
Seismic Data ◽  
Prediction Error ◽  
Reduction Method ◽  
Signal To Noise Ratio ◽  
Data Interpolation ◽  
Computationally Efficient ◽  
Rank Reduction ◽  
Frequency Space ◽  
Time Space ◽  
Ill Posed

The prediction error filter (PEF) assumes that the seismic data can be destructed to zero by applying a convolutional operation between the target data and prediction filter in either time-space or frequency-space domain. Here, we extend the commonly known PEF in 2D or 3D problems to its 5D version. To handle the non-stationary property of the seismic data, we formulate the PEF in a non-stationary way, which is called the non-stationary prediction error filter (NPEF). In the NPEF, the coefficients of a fixed-size PEF vary across the whole seismic data. In NPEF, we aim at solving a highly ill-posed inverse problem via the computationally efficient iterative shaping regularization. The NPEF can be used to denoise multi-dimensional seismic data, and more importantly, to restore the highly incomplete aliased 5D seismic data. We compare the proposed NPEF method with the state-of-the-art rank-reduction method for the 5D seismic data interpolation in cases of irregularly and regularly missing traces via several synthetic and real seismic data. Results show that although the proposed NPEF method is less effective than the rank-reduction method in interpolating irregularly missing traces especially in the case of low signal to noise ratio (S/N), it outperforms the rank-reduction method in interpolating aliased 5D dataset with regularly missing traces.

Multiple attenuation in the parabolic τ-p domain using wavefront characteristics of multiple generating primaries

Geophysics ◽  
1999 ◽  
Vol 64 (6) ◽  
pp. 1806-1815 ◽  
Author(s):  
Evgeny Landa ◽  
Igor Belfer ◽  
Shemer Keydar
Keyword(s):  
Radon Transform ◽  
Time Windows ◽  
Real Data ◽  
Multiple Model ◽  
Transform Method ◽  
Multiple Attenuation ◽  
The Common ◽  
P Domain ◽  
Multiple Prediction ◽  
Parabolic Radon Transform

The problem of multiple attenuation has been solved only partially. One of the most common methods of attenuating multiples is an approach based on the Radon transform. It is commonly accepted that the parabolic Radon transform method is only able to attenuate multiples with significant moveouts. We propose a new 2-D method for attenuation of both surface‐related and interbed multiples in the parabolic τ-p domain. The method is based on the prediction of a multiple model from the wavefront characteristics of the primary events. Multiple prediction comprises the following steps: 1) For a given multiple code, the angles of emergence and the radii of wavefront curvatures are estimated for primary reflections for each receiver in the common‐shotpoint gather. 2) The intermediate points which compose a specified multiple event are determined for each shot‐receiver pair. 3) Traveltimes of the multiples are calculated. Wavefields within time windows around the predicted traveltime curves may be considered as multiple model traces which we use for multiple attenuation process. Using the predicted multiple traveltimes, we can define the area in the τ-p domain which contains the main energy of the multiple event. Resolution improvement of the parabolic Radon operator can be achieved through a simple multiplication of each sample in the τ-p space by a nonlinear semblance function. In this work, we follow the idea of defining the multiple reject areas automatically by comparing the energy of the multiple model and the original input data in the τ-p space. We illustrate the usefulness of this algorithm for the attenuation of multiples on both synthetic and real data.

scholarly journals Multiple prediction through inversion: Theoretical advancements and real data application

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. V33-V39 ◽  
Author(s):  
Yanghua Wang
Keyword(s):  
Explicit Knowledge ◽  
Real Data ◽  
Second Phase ◽  
Multiple Model ◽  
Subsurface Structures ◽  
Multiple Attenuation ◽  
Data Application ◽  
Multiple Prediction ◽  
Marine Seismic ◽  
Surface Operator

Wave-equation-based multiple attenuation seismic methods may be divided into the two distinct phases of multiple modeling and multiple subtraction. These two are interrelated and must be optimized in order to produce an optimal final result. The multiple prediction through inversion (MPI) scheme updates the multiple model iteratively, as we usually do in a linearized inverse problem. The scheme models the multiple wavefield without an explicit knowledge of surface and subsurface structures or of the source signature; both are generally unknown in seismic surveys. However, compared to a conventional surface-related multiple attenuation method, the accuracy of the multiple model is improved both kinematically and dynamically. It is because the MPI scheme implicitly takes account of the spatial variation of the surface reflectivity, the source signature, the detector patterns and receiver ghosts, and other effects included in the so-called surface operator. When the MPI scheme is used in the first phase it also significantly reduces the nonlinearity of the problem in the second phase that involves attenuating multiples without removing or altering primaries. The effectiveness of the MPI scheme is demonstrated by examples involving real marine seismic data.

Sign in / Sign up

Export Citation Format

Share Document