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.

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 ATENUASI WATER-BOTTOM MULTIPLE DENGAN METODE TRANSFORMASI PARABOLIC RADON

2016 ◽  
Vol 12 (3) ◽  
pp. 145
Author(s):  
Subarsyah Subarsyah ◽  
Tumpal Benhard Nainggolan
Keyword(s):  
Radon Transform ◽  
Seismic Data ◽  
Wide Distribution ◽  
Seismic Section ◽  
Optimal Separation ◽  
Multiple Attenuation ◽  
Parabolic Radon Transform ◽  
Water Bottom

Interferensi water-bottom multipel terhadap reflektor primer menimbulkan efek bersifat destruktif yang menyebabkan penampang seismik menjadi tidak tepat akibat kehadiran reflektor semu. Teknik demultiple perlu diaplikasikan untuk mengatenuasi multipel. Transformasi parabolic radon merupakan teknik atenuasi multipel dengan metode pemisahan dalam domain radon. Multipel sering teridentifikasi pada penampang seismik. Untuk memperbaiki penampang seismik akan dilakukan dengan metode transformasi parabolic radon. Penerapan metode ini mengakibatkan reflektor multipel melemah dan tereduksi setelah dilakukan muting dalam domain radon terhadap zona multipel. Beberapa reflektor primer juga ikut melemah akibat pemisahan dalam domain radon yang kurang optimal, pemisahan akan optimal membutuhkan distribusi offset yang lebar. Kata kunci: Parabolic radon, multipel, atenuasi Water-bottom mutiple interference often destructively interfere with primary reflection that led to incorrect seismic section due to presence apparent reflector. Demultiple techniques need to be applied to attenuate the multiple. Parabolic Radon transform is demultiple attenuation technique that separate multiple and primary in radon domain. Water-bottom mutiple ussualy appear and easly identified on seismic data, parabolic radon transform applied to improve the seismic section. Application of this method to data showing multiple reflectors weakened and reduced after muting multiple zones in the radon domain. Some of the primary reflector also weakened due to bad separation in radon domain, optimal separation will require a wide distribution of offsets. Keywords: Parabolic radon, multiple, attenuation

Relative time seislet transform

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. V223-V232 ◽  
Author(s):  
Zhicheng Geng ◽  
Xinming Wu ◽  
Sergey Fomel ◽  
Yangkang Chen
Keyword(s):  
Seismic Data ◽  
Computational Cost ◽  
Lifting Scheme ◽  
Real Data ◽  
Data Sets ◽  
Relative Time ◽  
New Approach ◽  
New Formulation ◽  
Wavelet Lifting ◽  
Seismic Images

The seislet transform uses the wavelet-lifting scheme and local slopes to analyze the seismic data. In its definition, the designing of prediction operators specifically for seismic images and data is an important issue. We have developed a new formulation of the seislet transform based on the relative time (RT) attribute. This method uses the RT volume to construct multiscale prediction operators. With the new prediction operators, the seislet transform gets accelerated because distant traces get predicted directly. We apply our method to synthetic and real data to demonstrate that the new approach reduces computational cost and obtains excellent sparse representation on test data sets.

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.

Multiple attenuation using λ-f domain high-order and high-resolution Radon transform based on SL0 norm

2019 ◽  
Vol 16 (4) ◽  
pp. 473-482
Author(s):  
Wen-Zhi Sun ◽  
Zhen-Chun Li ◽  
Ying-Ming Qu ◽  
Zhi-Na Li
Keyword(s):  
High Resolution ◽  
Radon Transform ◽  
High Order ◽  
Multiple Attenuation

Revisiting levees in southern Texas using Love-wave multichannel analysis of surface waves with the high-resolution linear Radon transform

2017 ◽  
Vol 5 (3) ◽  
pp. T287-T298 ◽  
Author(s):  
Julian Ivanov ◽  
Richard D. Miller ◽  
Daniel Feigenbaum ◽  
Sarah L. C. Morton ◽  
Shelby L. Peterie ◽  
...  
Keyword(s):  
High Resolution ◽  
Rayleigh Wave ◽  
Surface Waves ◽  
Radon Transform ◽  
Seismic Data ◽  
Love Wave ◽  
S Wave ◽  
Multichannel Analysis ◽  
Masw Method ◽  
Wave Methods

Shear-wave velocities were estimated at a levee site by inverting Love waves using the multichannel analysis of surface waves (MASW) method augmented with the high-resolution linear Radon transform (HRLRT). The selected site was one of five levee sites in southern Texas chosen for the evaluation of several seismic data-analysis techniques readily available in 2004. The methods included P- and S-wave refraction tomography, Rayleigh- and Love-wave surface-wave analysis using MASW, and P- and S-wave cross-levee tomography. The results from the 2004 analysis revealed that although the P-wave methods provided reasonable and stable results, the S-wave methods produced surprisingly inconsistent shear-wave velocity [Formula: see text] estimates and trends compared with previous studies and borehole investigations. In addition, the Rayleigh-wave MASW method was nearly useless within the levee due to the sparsity of high frequencies in fundamental-mode surface waves and complexities associated with inverting higher modes. This prevented any reliable [Formula: see text] estimates for the levee core. Recent advances in methodology, such as the HRLRT for obtaining higher resolution dispersion-curve images with the MASW method and the use of Love-wave inversion routines specific to Love waves as part of the MASW method, provided the motivation to extend the 2004 original study by using horizontal-component seismic data for characterizing the geologic properties of levees. Contributions from the above-mentioned techniques were instrumental in obtaining [Formula: see text] estimates from within these levees that were very comparable with the measured borehole samples. A Love-wave approach can be a viable alternative to Rayleigh-wave MASW surveys at sites where complications associated with material or levee geometries inhibit reliable [Formula: see text] results from Rayleigh waves.

The Shuey-Radon transform

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. V197-V206 ◽  
Author(s):  
Ali Gholami ◽  
Milad Farshad
Keyword(s):  
Radon Transform ◽  
Seismic Data ◽  
Matching Pursuit ◽  
Real Data ◽  
Transform Domain ◽  
Amplitude Variation ◽  
Amplitude Variation With Offset ◽  
Linear System Of Equations ◽  
Matching Pursuit Algorithm ◽  
Zero Offset

The traditional hyperbolic Radon transform (RT) decomposes seismic data into a sum of constant amplitude basis functions. This limits the performance of the transform when dealing with real data in which the reflection amplitudes include the amplitude variation with offset (AVO) variations. We adopted the Shuey-Radon transform as a combination of the RT and Shuey’s approximation of reflectivity to accurately model reflections including AVO effects. The new transform splits the seismic gather into three Radon panels: The first models the reflections at zero offset, and the other two panels add capability to model the AVO gradient and curvature. There are two main advantages of the Shuey-Radon transform over similar algorithms, which are based on a polynomial expansion of the AVO response. (1) It is able to model reflections more accurately. This leads to more focused coefficients in the transform domain and hence provides more accurate processing results. (2) Unlike polynomial-based approaches, the coefficients of the Shuey-Radon transform are directly connected to the classic AVO parameters (intercept, gradient, and curvature). Therefore, the resulting coefficients can further be used for interpretation purposes. The solution of the new transform is defined via an underdetermined linear system of equations. It is formulated as a sparsity-promoting optimization, and it is solved efficiently using an orthogonal matching pursuit algorithm. Applications to different numerical experiments indicate that the Shuey-Radon transform outperforms the polynomial and conventional RTs.

Direct estimation of shear‐wave interval velocities from seismic data

Geophysics ◽  
1985 ◽  
Vol 50 (4) ◽  
pp. 530-538 ◽  
Author(s):  
P. M. Carrion ◽  
S. Hassanzadeh
Keyword(s):  
Shear Wave ◽  
Small Angle ◽  
Seismic Data ◽  
Mode Conversion ◽  
Real Data ◽  
Compressional Wave ◽  
Compressional Waves ◽  
Common Depth Point ◽  
Interval Velocities ◽  
Wave Decomposition

Conventional velocity analysis of seismic data is based on normal moveout of common‐depth‐point (CDP) traveltime curves. Analysis is done in a hyperbolic framework and, therefore, is limited to using the small‐angle reflections only (muted data). Hence, it can estimate the interval velocities of compressional waves only, since mode conversion is negligible when small‐angle arrivals are concerned. We propose a new method which can estimate the interval velocities of compressional and mode‐converted waves separately. The method is based on slant stacking or plane‐wave decomposition (PWD) of the observed data (seismogram), which transforms the data from the conventional T-X domain into the intercept time‐ray parameter domain. Since PWD places most of the compressional energy into the precritical region of the slant‐stacked seismogram, the compressional‐wave interval velocities can be estimated using the “best ellipse” approximation on the assumption that the elliptic array velocity (stacking velocity) is approximately equal to the root‐mean‐square (rms) velocity. Similarly, shear‐wave interval velocities can be estimated by inverting the traveltime curves in the region of the PWD seismogram, where compressional waves decay exponentially (postcritical region). The method is illustrated by examples using synthetic and real data.

Sign in / Sign up

Export Citation Format

Share Document