Three-parameter Radon transform based on shifted hyperbolas

Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. V39-V48 ◽  
Author(s):  
Ali Gholami ◽  
Toktam Zand

The focusing power of the conventional hyperbolic Radon transform decreases for long-offset seismic data due to the nonhyperbolic behavior of moveout curves at far offsets. Furthermore, conventional Radon transforms are ineffective for processing data sets containing events of different shapes. The shifted hyperbola is a flexible three-parameter (zero-offset traveltime, slowness, and focusing-depth) function, which is capable of generating linear and hyperbolic shapes and improves the accuracy of the seismic traveltime approximation at far offsets. Radon transform based on shifted hyperbolas thus improves the focus of seismic events in the transform domain. We have developed a new method for effective decomposition of seismic data by using such three-parameter Radon transform. A very fast algorithm is constructed for high-resolution calculations of the new Radon transform using the recently proposed generalized Fourier slice theorem (GFST). The GFST establishes an analytic expression between the [Formula: see text] coefficients of the data and the [Formula: see text] coefficients of its Radon transform, with which a very fast switching between the model and data spaces is possible by means of interpolation procedures and fast Fourier transforms. High performance of the new algorithm is demonstrated on synthetic and real data sets for trace interpolation and linear (ground roll) noise attenuation.

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. V197-V206 ◽  
Author(s):  
Ali Gholami ◽  
Milad Farshad

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.


Geophysics ◽  
2014 ◽  
Vol 79 (1) ◽  
pp. V1-V11 ◽  
Author(s):  
Amr Ibrahim ◽  
Mauricio D. Sacchi

We adopted the robust Radon transform to eliminate erratic incoherent noise that arises in common receiver gathers when simultaneous source data are acquired. The proposed robust Radon transform was posed as an inverse problem using an [Formula: see text] misfit that is not sensitive to erratic noise. The latter permitted us to design Radon algorithms that are capable of eliminating incoherent noise in common receiver gathers. We also compared nonrobust and robust Radon transforms that are implemented via a quadratic ([Formula: see text]) or a sparse ([Formula: see text]) penalty term in the cost function. The results demonstrated the importance of incorporating a robust misfit functional in the Radon transform to cope with simultaneous source interferences. Synthetic and real data examples proved that the robust Radon transform produces more accurate data estimates than least-squares and sparse Radon transforms.


Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. V223-V232 ◽  
Author(s):  
Zhicheng Geng ◽  
Xinming Wu ◽  
Sergey Fomel ◽  
Yangkang Chen

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.


Geophysics ◽  
1996 ◽  
Vol 61 (6) ◽  
pp. 1846-1858 ◽  
Author(s):  
Claudio Bagaini ◽  
Umberto Spagnolini

Continuation to zero offset [better known as dip moveout (DMO)] is a standard tool for seismic data processing. In this paper, the concept of DMO is extended by introducing a set of operators: the continuation operators. These operators, which are implemented in integral form with a defined amplitude distribution, perform the mapping between common shot or common offset gathers for a given velocity model. The application of the shot continuation operator for dip‐independent velocity analysis allows a direct implementation in the acquisition domain by exploiting the comparison between real data and data continued in the shot domain. Shot and offset continuation allow the restoration of missing shot or missing offset by using a velocity model provided by common shot velocity analysis or another dip‐independent velocity analysis method.


Geophysics ◽  
2012 ◽  
Vol 77 (1) ◽  
pp. A5-A8 ◽  
Author(s):  
David Bonar ◽  
Mauricio Sacchi

The nonlocal means algorithm is a noise attenuation filter that was originally developed for the purposes of image denoising. This algorithm denoises each sample or pixel within an image by utilizing other similar samples or pixels regardless of their spatial proximity, making the process nonlocal. Such a technique places no assumptions on the data except that structures within the data contain a degree of redundancy. Because this is generally true for reflection seismic data, we propose to adopt the nonlocal means algorithm to attenuate random noise in seismic data. Tests with synthetic and real data sets demonstrate that the nonlocal means algorithm does not smear seismic energy across sharp discontinuities or curved events when compared to seismic denoising methods such as f-x deconvolution.


Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 225-231 ◽  
Author(s):  
Rongfeng Zhang ◽  
Tadeusz J. Ulrych

This paper deals with the design and implementation of a new wavelet frame for noise suppression based on the character of seismic data. In general, wavelet denoising methods widely used in image and acoustic processing use well‐known conventional wavelets which, although versatile, are often not optimal for seismic data. The new approach, physical wavelet frame denoising uses a wavelet frame that takes into account the characteristics of seismic data both in time and space. Synthetic and real data tests show that the approach is effective even for seismic signals contaminated by strong noise which may be random or coherent, such as ground roll or air waves.


Author(s):  
Maxim I. Protasov ◽  
◽  
Vladimir A. Tcheverda ◽  
Valery V. Shilikov ◽  
◽  
...  

The paper deals with a 3D diffraction imaging with the subsequent diffraction attribute calculation. The imaging is based on an asymmetric summation of seismic data and provides three diffraction attributes: structural diffraction attribute, point diffraction attribute, an azimuth of structural diffraction. These attributes provide differentiating fractured and cavernous objects and to determine the fractures orientations. Approbation of the approach was provided on several real data sets.


Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. R59-R67 ◽  
Author(s):  
Igor B. Morozov ◽  
Jinfeng Ma

The seismic-impedance inversion problem is underconstrained inherently and does not allow the use of rigorous joint inversion. In the absence of a true inverse, a reliable solution free from subjective parameters can be obtained by defining a set of physical constraints that should be satisfied by the resulting images. A method for constructing synthetic logs is proposed that explicitly and accurately satisfies (1) the convolutional equation, (2) time-depth constraints of the seismic data, (3) a background low-frequency model from logs or seismic/geologic interpretation, and (4) spectral amplitudes and geostatistical information from spatially interpolated well logs. The resulting synthetic log sections or volumes are interpretable in standard ways. Unlike broadly used joint-inversion algorithms, the method contains no subjectively selected user parameters, utilizes the log data more completely, and assesses intermediate results. The procedure is simple and tolerant to noise, and it leads to higher-resolution images. Separating the seismic and subseismic frequency bands also simplifies data processing for acoustic-impedance (AI) inversion. For example, zero-phase deconvolution and true-amplitude processing of seismic data are not required and are included automatically in this method. The approach is applicable to 2D and 3D data sets and to multiple pre- and poststack seismic attributes. It has been tested on inversions for AI and true-amplitude reflectivity using 2D synthetic and real-data examples.


Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. V69-V78 ◽  
Author(s):  
Jinlin Liu ◽  
Wenkai Lu

Adaptive multiple subtraction is the key step of surface-related multiple elimination methods. The main challenge of this technique resides in removing multiples without distorting primaries. We have developed a new pattern-based method for adaptive multiple subtraction with the consideration that primaries can be better protected if the multiples are differentiated from the primaries. Different from previously proposed methods, our method casts the adaptive multiple subtraction problem as a pattern coding and decoding process. We set out to learn distinguishable patterns from the predicted multiples before estimating the multiples contained in seismic data. Hence, in our method, we first carried out pattern coding of the predicted multiples to learn the special patterns of the multiples within different frequency bands. This coding process aims at exploiting the key patterns contained in the predicted multiples. The learned patterns are then used to decode (extract) the multiples contained in the seismic data, in which process those patterns that are similar to the learned patterns were identified and extracted. Because the learned patterns are obtained from the predicted multiples only, our method avoids the interferences of primaries in nature and shows an impressive capability for removing multiples without distorting the primaries. Our applications on synthetic and real data sets gave some promising results.


Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. V317-V328
Author(s):  
Jitao Ma ◽  
Guoyang Xu ◽  
Xiaohong Chen ◽  
Xiaoliu Wang ◽  
Zhenjiang Hao

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.


Sign in / Sign up

Export Citation Format

Share Document