Efficient elastic reverse-time migration for the decomposed P-wavefield using stress tensor in the time domain

2015 ◽  
Vol 116 ◽  
pp. 121-134 ◽  
Author(s):  
Jiho Ha ◽  
Sungryul Shin ◽  
Changsoo Shin ◽  
Wookeen Chung
Keyword(s):  
Stress Tensor ◽  
Time Domain ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
The Time Domain

Q-compensated reverse time migration in tilted transversely isotropic media

Geophysics ◽  
2020 ◽  
pp. 1-79
Author(s):  
Ali Fathalian ◽  
Daniel O. Trad ◽  
Kristopher A. Innanen
Keyword(s):  
Time Domain ◽  
Transversely Isotropic ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Modeling And Analysis ◽  
Seismic Amplitude ◽  
Time Migration ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
The Time Domain

Anisotropy and absorption are critical to the modeling and analysis of seismic amplitude,phase, and traveltime data. To neglect any of these phenomena, which are often bothoperating simultaneously, degrades the resolution and interpretability of migrated images.However, a full accounting of anisotropy and anelasticity is computationally complex andexpensive. One strategy for accommodating these aspects of wave propagation, while keepingcost and complexity under control, is to do so within an acoustic approximation. Weset up a procedure for solving the time-domain viscoacoustic wave equation for tilted transverselyisotropic (TTI) media, based on a standard linear solid model and, from this, developa viscoacoustic reverse time migration (Q-RTM) algorithm. In this approach, amplitudecompensation occurs within the migration process through a manipulation of attenuationand phase dispersion terms in the time domain differential equations. Specifically, theback-propagation operator is constructed by reversing the sign only of the amplitude lossoperators, but not the dispersion-related operators, a step made possible by reformulatingthe absorptive TTI equations such that the loss and dispersion operators appear separately.The scheme is tested on synthetic examples to examine the capacity of viscoacoustic RTM to correct for attenuation, and the overall stability of the procedure.

Angle gathers from reverse time migration using analytic wavefield propagation and decomposition in the time domain

Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. S1-S9 ◽  
Author(s):  
Jiangtao Hu ◽  
Huazhong Wang ◽  
Xiongwen Wang
Keyword(s):  
Plane Wave ◽  
Time Domain ◽  
Computational Cost ◽  
Imaging Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Imaging Condition ◽  
Time Migration ◽  
Local Plane ◽  
The Time Domain

Angle-domain common imaging gathers (ADCIGs) are important input data for migration velocity analysis and amplitude variation with angle analysis. Compared with Kirchhoff migration and one-way wave equation migration, reverse time migration (RTM) is the most accurate imaging method in complex areas, such as the subsalt area. We have developed a method to generate ADCIGs from RTM using analytic wavefield propagation and decomposition. To estimate the wave-propagation direction and angle by spatial Fourier transform during the time domain wave extrapolation, we have developed an analytic wavefield extrapolation method. Then, we decomposed the extrapolated source and receiver wavefields into their local angle components (i.e., local plane-wave components) and applied the angle-domain imaging condition to form ADCIGs. Because the angle-domain imaging condition is a convolution imaging condition about the source and receiver propagation angles, it is costly. To increase the efficiency of the angle-domain imaging condition, we have developed a local plane-wave decomposition method using matching pursuit. Numerical examples of synthetic and real data found that this method could generate high-quality ADCIGs. And these examples also found that the computational cost of this approach was related to the complexity of the source and receiver wavefields.

Fast plane-wave reverse time migration

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S549-S556 ◽  
Author(s):  
Xiongwen Wang ◽  
Xu Ji ◽  
Hongwei Liu ◽  
Yi Luo
Keyword(s):  
Time Delay ◽  
Plane Wave ◽  
Green’S Function ◽  
Green's Function ◽  
Time Domain ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Wave Source ◽  
Computation Cost ◽  
Time Migration

Plane-wave reverse time migration (RTM) could potentially provide quick subsurface images by migrating fewer plane-wave gathers than shot gathers. However, the time delay between the first and the last excitation sources in the plane-wave source largely increases the computation cost and decreases the practical value of this method. Although the time delay problem is easily overcome by periodical phase shifting in the frequency domain for one-way wave-equation migration, it remains a challenge for time-domain RTM. We have developed a novel method, referred as to fast plane-wave RTM (FP-RTM), to eliminate unnecessary computation burden and significantly reduce the computational cost. In the proposed FP-RTM, we assume that the Green’s function has finite-length support; thus, the plane-wave source function and its responding data can be wrapped periodically in the time domain. The wrapping length is the assumed total duration length of Green’s function. We also determine that only two period plane-wave source and data after the wrapping process are required for generating the outcome with adequate accuracy. Although the computation time for one plane-wave gather is twice as long as a normal shot gather migration, a large amount of computation cost is saved because the total number of plane-wave gathers to be migrated is usually much less than the total number of shot gathers. Our FP-RTM can be used to rapidly generate RTM images and plane-wave domain common-image gathers for velocity model building. The synthetic and field data examples are evaluated to validate the efficiency and accuracy of our method.

Reverse Time Migration of Elastic Waves Using the Pseudospectral Time-Domain Method

2018 ◽  
Vol 26 (01) ◽  
pp. 1750033 ◽  
Author(s):  
Jiangang Xie ◽  
Mingwei Zhuang ◽  
Zichao Guo ◽  
Hai Liu ◽  
Qing Huo Liu
Keyword(s):  
Elastic Wave ◽  
Time Domain ◽  
Elastic Waves ◽  
Sampling Rate ◽  
Fdtd Method ◽  
Time Domain Method ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Domain Method

Reverse time migration (RTM), especially that for elastic waves, consumes massive computation resources which limit its wide applications in industry. We suggest to use the pseudospectral time-domain (PSTD) method in elastic wave RTM. RTM using PSTD can significantly reduce the computational requirements compared with RTM using the traditional finite difference time domain method (FDTD). In addition to the advantage of low sampling rate with high accuracy, the PSTD method also eliminates the periodicity (or wraparound) limitation caused by fast Fourier transform in the conventional pseudospectral method. To achieve accurate results, the PSTD method needs only about half the spatial sampling rate of the twelfth-order FDTD method. Thus, the PSTD method can save up to 87.5% storage memory and 90% computation time over the twelfth-order FDTD method. We implement RTM using PSTD for elastic wave equations and accelerate it by Open Multi-Processing technology. To keep the computational load balance in parallel computation, we design a new PML layout which merges the PML in both ends of an axis together. The efficiency and imaging quality of the proposed RTM is verified by imaging on 2D and 3D models.

Reverse Time Migration Using the Pseudospectral Time-Domain Algorithm

2016 ◽  
Vol 24 (02) ◽  
pp. 1650005 ◽  
Author(s):  
Jiangang Xie ◽  
Zichao Guo ◽  
Hai Liu ◽  
Qing Huo Liu
Keyword(s):  
Time Domain ◽  
Acoustic Waves ◽  
Fdtd Method ◽  
Pseudospectral Method ◽  
Imaging Method ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Eighth Order ◽  
Time Migration

We propose a pre-stack reverse time migration (RTM) seismic imaging method using the pseudospectral time-domain (PSTD) algorithm. Traditional pseudospectral method uses the fast Fourier transform (FFT) algorithm to calculate the spatial derivatives, but is limited by the wraparound effect due to the periodicity assumed in the FFT. The PSTD algorithm combines the pseudospectral method with a perfectly matched layer (PML) for acoustic waves. PML is a highly effective absorbing boundary condition that can eliminate the wraparound effect. It enables a wide application of the pseudospectral method to complex models. RTM based on the PSTD algorithm has advantages in the computational efficiency compared to traditional methods such as the second-order and high order finite difference time-domain (FDTD) methods. In this work, we implement the PSTD algorithm for acoustic wave equation based RTM. By applying the PSTD-RTM method to various seismic models and comparing it with RTM based on the eighth-order FDTD method, we find that PSTD-RTM method has better performance and saves more than 50% memory. The method is suitable for parallel computation, and has been accelerated by general purpose graphics processing unit.

Sign in / Sign up

Export Citation Format

Share Document