An implicit stabilization strategy for Q-compensated reverse time migration

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. S169-S183
Author(s):  
Hanming Chen ◽  
Hui Zhou ◽  
Ying Rao
Keyword(s):  
Robust Stabilization ◽  
Synthetic Data ◽  
Frequency Noise ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Time Migration ◽  
Wavefield Extrapolation ◽  
Amplitude Compensation ◽  
Seismic Images

Reverse time migration with [Formula: see text] compensation ([Formula: see text]-RTM) is an effective approach to enhance the resolution of seismic images because it retrieves the amplitude loss and phase distortion induced by the viscosity of media. According to the crosscorrelation imaging condition, [Formula: see text]-RTM requires compensation for the amplitude loss in the propagation paths of source and receiver wavefields, which can be realized by solving an amplitude-boosted wave equation. However, the amplitude-boosted simulations suffer from numerical instability due to the amplification of high-frequency noise. We have developed a robust stabilization strategy for [Formula: see text]-RTM by incorporating a time-variant filter into the amplitude-boosted wavefield extrapolation step. We modify the Fourier spectrum of the operator that controls the amplitude compensation to be time variant, and we add to the spectrum a stabilization factor. Doing so, we integrate the time-variant filter into the viscoacoustic wave propagator implicitly, and we avoid any explicit filtering operation in [Formula: see text]-RTM. We verify the robustness of this stabilized [Formula: see text]-RTM with two synthetic data examples. We also apply this technique to a field data set to demonstrate the imaging improvements compared to an acoustic RTM and a more traditional [Formula: see text]-RTM method.

Reverse time migration of multiples: Reducing migration artifacts using the wavefield decomposition imaging condition

Geophysics ◽  
2017 ◽  
Vol 82 (4) ◽  
pp. S307-S314 ◽  
Author(s):  
Yibo Wang ◽  
Yikang Zheng ◽  
Qingfeng Xue ◽  
Xu Chang ◽  
Tong W. Fei ◽  
...  
Keyword(s):  
Complex Function ◽  
Synthetic Data ◽  
Velocity Model ◽  
High Amplitude ◽  
Frequency Noise ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Time Migration ◽  
Wavefield Decomposition

In the implementation of migration of multiples, reverse time migration (RTM) is superior to other migration algorithms because it can handle steeply dipping structures and offer high-resolution images of the complex subsurface. However, the RTM results using two-way wave equation contain high-amplitude, low-frequency noise and false images generated by improper wave paths in migration velocity model with sharp velocity interfaces or strong velocity gradients. To improve the imaging quality in RTM of multiples, we separate the upgoing and downgoing waves in the propagation of source and receiver wavefields. A complex function involved with the Hilbert transform is used in wavefield decomposition. Our approach is cost effective and avoids the large storage of wavefield snapshots required by the conventional wavefield separation technique. We applied migration of multiples with wavefield decomposition on a simple two-layer model and the Sigsbee 2B synthetic data set. Our results demonstrate that the proposed approach can improve the image generated by migration of multiples significantly.

Acoustic and elastic numerical wave simulations by recursive spatial derivative operators

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. T167-T174 ◽  
Author(s):  
Dan Kosloff ◽  
Reynam C. Pestana ◽  
Hillel Tal-Ezer
Keyword(s):  
Synthetic Data ◽  
Analytic Solutions ◽  
Numerical Dispersion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Lamb’S Problem ◽  
Dynamic Elasticity ◽  
Time Migration ◽  
Recursive Operators

A new scheme for the calculation of spatial derivatives has been developed. The technique is based on recursive derivative operators that are generated by an [Formula: see text] fit in the spectral domain. The use of recursive operators enables us to extend acoustic and elastic wave simulations to shorter wavelengths. The method is applied to the numerical solution of the 2D acoustic wave equation and to the solution of the equations of 2D dynamic elasticity in an isotropic medium. An example of reverse-time migration of a synthetic data set shows that the numerical dispersion can be significantly reduced with respect to schemes that are based on finite differences. The method is tested for the solutions of the equations of dynamic elasticity by comparing numerical and analytic solutions to Lamb’s problem.

scholarly journals LAGUERRE-GAUSS FILTERS IN REVERSE TIME MIGRATION IMAGE RECONSTRUCTION

2018 ◽  
Vol 35 (2) ◽  
Author(s):  
Juan Guillermo Paniagua Castrillón ◽  
Olga Lucia Quintero Montoya ◽  
Daniel Sierra-Sosa
Keyword(s):  
Cross Correlation ◽  
Synthetic Data ◽  
Low Frequency ◽  
Frequency Noise ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Migration Imaging ◽  
Gauss Transform

ABSTRACT. Reverse time migration (RTM) solves the acoustic or elastic wave equation by means of the extrapolation from source and receiver wavefield in time. A migrated image is obtained by applying a criteria known as imaging condition. The cross-correlation between source and receiver wavefields is the commonly used imaging condition. However, this imaging condition produces spatial low-frequency noise, called artifacts, due to the unwanted correlation of the diving, head and backscattered waves. Several techniques have been proposed to reduce the artifacts occurrence. Derivative operators as Laplacian are the most frequently used. In this work, we propose a technique based on a spiral phase filter ranging from 0 to 2π, and a toroidal amplitude bandpass filter, known as Laguerre-Gauss transform. Through numerical experiments we present the application of this particular filter on three synthetic data sets. In addition, we present a comparative spectral study of images obtained by the zero-lag cross-correlation imaging condition, the Laplacian filtering and the Laguerre-Gauss filtering, showing their frequency features. We also present evidences not only with simulated noisy velocity fields but also by comparison with the model velocity field gradients that this method improves the RTM images by reducing the artifacts and notably enhance the reflective events. Keywords: Laguerre-Gauss transform, zero-lag cross-correlation, seismic migration, imaging condition. RESUMO. A migração reversa no tempo (RTM) resolve a equação de onda acústica ou elástica por meio da extrapolação a partir do campo de onda da fonte e do receptor no tempo. Uma imagem migrada é obtida aplicando um critério conhecido como condição de imagem. A correlação cruzada entre campos de onda de fonte e receptor é a condição de imagem comumente usada. No entanto, esta condição de imagem produz ruído espacial de baixa frequência, chamados artefatos, devido à correlação indesejada das ondas de mergulho, cabeça e retrodifusão. Várias técnicas têm sido propostas para reduzir a ocorrência de artefatos. Operadores derivados como Laplaciano são os mais utilizados. Neste trabalho, propomos uma técnica baseada em um filtro de fase espiral que varia de 0 a 2π, e um filtro passabanda de amplitude toroidal, conhecido como transformada de Laguerre-Gauss. Através de experimentos numéricos, apresentamos a aplicação deste filtro particular em três conjuntos de dados sintéticos. Além disso, apresentamos um estudo comparativo espectral de imagens obtidas pela condição de imagem de correlação cruzada atraso zero, a filtragem de Laplaciano e a filtragem Laguerre-Gauss, mostrando suas características de frequência. Apresentamos evidências não somente com campos simulados de velocidade ruidosa, mas também por comparação com os gradientes de campo de velocidade do modelo que este método melhora as imagens RTM, reduzindo os artefatos e aumentando notavelmente os eventos reflexivos. Palavras-chave: Transformação de Laguerre-Gauss, correlação cruzada atraso zero, migração sísmica, condição de imagem.

Viscoacoustic reverse-time migration with a robust space-wavenumber domain attenuation compensation operator

Geophysics ◽  
2021 ◽  
pp. 1-95
Author(s):  
Jidong Yang ◽  
Jianping Huang ◽  
Hejun Zhu ◽  
Zhenchun Li ◽  
Nanxun Dai
Keyword(s):  
Time Window ◽  
Back Propagation ◽  
Homogeneous Medium ◽  
Frequency Noise ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Frequency ◽  
Attenuation Compensation ◽  
Time Migration ◽  
Wavefield Extrapolation

Intrinsic attenuation gives rise to phase dispersion and amplitude loss during seismic wave propagation. Not correcting these effects in seismic imaging can result in inaccurate reflector locations, dimmed amplitudes and degraded spatial resolution. In reverse-time migration (RTM), attenuation compensation can be implemented by reversing the sign of the dissipation term and keeping the dispersion term unchanged for backward wavefield extrapolation. Although this Q-compensated RTM scheme can effectively correct attenuation effects, amplitude amplification during back-propagation might lead to numerical instabilities, especially for field data with strong high-frequency noise. To mitigate this problem, we develop a robust space-wavenumber compensation operator, and apply it to viscoacoustic RTM. By analyzing the dispersion-only and viscoacoustic Green’s functions, we obtain an analytical solution for the attenuation compensation operator in a homogeneous medium. Because it is a time-frequency operator, to apply it directly in viscoacoustic RTM requires access to the extrapolated wavefields within a certain time window. To avoid storing the wavefields and improve computational efficiency, we use an approximated dispersion relation and convert the time-frequency operator to an equivalent space-wavenumber operator, which allows us to implement attenuation compensation on the fly during wavefield extrapolation. The hybrid-domain property of the operator enables us to account for the wavenumber-dependent compensation. A similar strategy can also be applied to the migrated images as a poststack processing approach, which is more efficient than the prestack compensation. Two synthetic and one land field dataset examples demonstrate the feasibility and adaptability of the proposed method.

Fast 3D target-oriented reverse-time datuming

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA141-WCA151 ◽  
Author(s):  
Shuqian Dong ◽  
Yi Luo ◽  
Xiang Xiao ◽  
Sergio Chávez-Pérez ◽  
Gerard T. Schuster
Keyword(s):  
Synthetic Data ◽  
Data Sets ◽  
Target Area ◽  
Complex Structures ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Salt Domes ◽  
Time Migration ◽  
Subsalt Imaging

Imaging of subsalt sediments is a challenge for traditional migration methods such as Kirchhoff and one-way wave-equation migration. Consequently, the more accurate two-way method of reverse-time migration (RTM) is preferred for subsalt imaging, but its use can be limited by high computation cost. To overcome this problem, a 3D target-oriented reverse-time datuming (RTD) method is presented, which can generate redatumed data economically in target areas beneath complex structures such as salt domes. The redatumed data in the target area then can be migrated inexpensively using a traditional migration method. If the target area is much smaller than the acquisition area, computation costs are reduced significantly by the use of a novel bottom-up strategy to calculate the extrapolated Green’s functions. Target-oriented RTD is tested on 2D and 3D SEG/EAGE synthetic data sets and a 3D field data set from the Gulf of Mexico. Results show that target-oriented RTD combined with standard migration can image sediments beneath complex structures accurately with much less calculation effort than full volume RTM. The requirement is that the area over the target zone is smaller than that of the acquisition survey.

Reverse time migration of multiples: Eliminating migration artifacts in angle domain common image gathers

Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. S263-S270 ◽  
Author(s):  
Yibo Wang ◽  
Yikang Zheng ◽  
Lele Zhang ◽  
Xu Chang ◽  
Zhenxing Yao
Keyword(s):  
Free Surface ◽  
Radon Transform ◽  
Synthetic Data ◽  
Velocity Model ◽  
Migration Velocity ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Time Migration ◽  
Nonzero Curvature

Free-surface-related multiples are usually regarded as noise in conventional seismic processing. However, they can provide extra illumination of the subsurface and thus have been used in migration procedures, e.g., in one- and two-way wave-equation migrations. The disadvantage of the migration of multiples is the migration artifacts generated by the crosscorrelation of different seismic events, e.g., primaries and second-order free-surface-related multiples, so the effective elimination of migration artifacts is crucial for migration of multiples. The angle domain common image gather (ADCIG) is a suitable domain for testing the correctness of a migration velocity model. When the migration velocity model is correct, all the events in ADCIGs should be flat, and this provides a criterion for removing the migration artifacts. Our approach first obtains ADCIGs during reverse time migration and then applies a high-resolution parabolic Radon transform to all ADCIGs. By doing so, most migration artifacts will reside in the nonzero curvature regions in the Radon domain, and then a muting procedure can be implemented to remove the data components outside the vicinity of zero curvature. After the application of an adjoint Radon transform, the filtered ADCIGs are obtained and the final denoised migration result is generated by stacking all filtered ADCIGs. A three-flat-layer velocity model and the Marmousi synthetic data set are used for numerical experiments. The numerical results revealed that the proposed approach can eliminate most artifacts generated by migration of multiples when the migration velocity model is correct.

Visibility analysis for optimal imaging of target areas and its application to a Gulf of Mexico deep-water data set

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB119-WB126 ◽  
Author(s):  
Elive Menyoli ◽  
Shengwen Jin ◽  
Shiyong Xu ◽  
Stuart Graber
Keyword(s):  
Gulf Of Mexico ◽  
Synthetic Data ◽  
Target Area ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Visibility Analysis ◽  
Velocity Models ◽  
Time Migration ◽  
Subsalt Imaging

Marine wide-azimuth data in the Gulf of Mexico, reverse time migration (RTM) and anisotropic velocity models have led to significant improvement in subsalt imaging. However, imaging of some steeply dipping subsalt targets such as three-way closures against salt is still difficult. This can be attributed to poor illumination and noise contaminations from various shot records. We apply the visibility analysis method that quantitatively determines which shot records contribute most energy on a specific subsalt prospect area. As a result we selectively migrate only those shot records thereby reducing noise contamination from low energy contributing shot records, improving signal continuity and better trap definition in the target area. Like conventional illumination analysis, the computation takes into account the overburden velocity distribution, acquisition geometry, target reflectivity and dip angle. We used 2D and 3D synthetic data examples to test the concepts and applicability of the method. A Gulf of Mexico case study example using wide-azimuth data demonstrated its use in an industry scale project. It is shown that for the particular 60°–65° subsalt target of interest only 30% of the wide-azimuth shot records are sufficient for the imaging. By reducing noise, the image results show significant improvement in the subsalt area compared to the full shot record RTM volume.

A stable and efficient approach of Q reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S557-S567 ◽  
Author(s):  
Yan Zhao ◽  
Ningbo Mao ◽  
Zhiming Ren
Keyword(s):  
Wave Equation ◽  
Computational Efficiency ◽  
High Frequency ◽  
Frequency Noise ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
High Frequency Noise ◽  
Time Migration ◽  
Low Pass ◽  
Wavefield Extrapolation

Amplitude energy attenuation and phase distortion of seismic waves caused by formation viscoelasticity reduce the resolution of reverse time migration (RTM) images. Q-RTM often is used to compensate the attenuation effects and improve the resolution of seismic imaging. However, serious high-frequency noise and tremendous amplitude will be produced during the wavefield extrapolation of Q-RTM, resulting in its inability to be imaged. Many Q-RTM algorithms solve the problem of instability through low-pass filtering in the wavenumber domain, but the method is less efficient in computation and has a truncation effect in the wavefield. We have developed a stable and efficient Q-RTM method, in which a regularization term was introduced into the viscoacoustic wave equation to suppress the high-frequency noise, and the finite-difference method was used to solve the viscoacoustic wave equation with a regularization term. We used the model example to visually demonstrate the instability of wavefield extrapolation in Q-RTM and compared the effect and computational efficiency of the two stabilization processing methods, low-pass filtering and regularization. Meanwhile, our method is not involved in solving the fractional derivatives by using the pseudo-spectral method, the computational efficiency also can be improved. We tested the Q-RTM approach on a simple layered model, Marmousi model, and real seismic data. The results of numerical examples demonstrated that the Q-RTM method can solve the problem of instability effectively and obtain a higher resolution image with lower computational cost.

Implementation of prestack reverse time migration using frequency-domain extrapolation

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. S61-S72 ◽  
Author(s):  
Kun Xu ◽  
Bing Zhou ◽  
George A. McMechan
Keyword(s):  
Frequency Domain ◽  
Synthetic Data ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Disk Storage ◽  
Time Migration ◽  
Seismic Experiment ◽  
Excitation Time ◽  
Parallel Iterative ◽  
Wavefield Extrapolation

We implement prestack reverse time migration (RTM) using frequency-domain extrapolation. We assume the observed data are the scattered field from the seismic experiment; then we time reverse them and derive the corresponding multipoint virtual sources to place along the receiver line rather than treat them as boundary conditions, which are common in time-domain reverse time migration (TRTM). Because the number of virtual sources is equal to the number of independent true sources in prestack frequency-domain reverse-time migration (FRTM), wavefield extrapolation is efficient in the frequency domain using a direct lower/upper (LU) triangular solver for 2D solutions and is still feasible using a parallel iterative or hybrid direct/iterative solver for 3D solutions. Each frequency implicitly contains information for all times, so we can implement the excitation-time and crosscorrelation imaging conditions straightforwardly without disk storage or input/output (I/O). Comparing costs, prestack shot-record FRTM is competitive with prestack shot-record TRTM. In our tests, prestack acoustic-wave FRTM accurately recovers the correct images of subsurface reflectors using 2D synthetic data.

Dirty salt velocity inversion: The road to a clearer subsalt image

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB169-WB174 ◽  
Author(s):  
Shuo Ji ◽  
Tony Huang ◽  
Kang Fu ◽  
Zhengxue Li
Keyword(s):  
Model Building ◽  
Synthetic Data ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Velocity Inversion ◽  
Time Migration ◽  
And Migration ◽  
Subsalt Imaging

For deep-water Gulf of Mexico, accurate salt geometry is critical to subsalt imaging. This requires the definition of both external and internal salt geometries. In recent years, external salt geometry (i.e., boundaries between allochthonous salt and background sediment) has improved a great deal due to advances in acquisition, velocity model building, and migration algorithms. But when it comes to defining internal salt geometry (i.e., intrasalt inclusions or dirty salt), no efficient method has yet been developed. In common industry practices, intrasalt inclusions (and thus their velocity anomalies) are generally ignored during the model building stages. However, as external salt geometries reach higher levels of accuracy, it becomes more important to consider the once-ignored effects of dirty salt. We have developed a reflectivity-based approach for dirty salt velocity inversion. This method takes true-amplitude reverse time migration stack volumes as input, then estimates the dirty salt velocity based on reflectivity under a 1D assumption. Results from a 2D synthetic data set and a real 3D Wide Azimuth data set demonstrated that the reflectivity inversion scheme significantly improves the subsalt image for certain areas. In general, we believe that this method produces a better salt model than the traditional clean salt velocity approach.

Sign in / Sign up

Export Citation Format

Share Document