An efficient and accurate finite-difference operator using adaptively discretized grids and its application for 3D least-squares reverse-time migration

Wang, Ziying; Huang, Jianping; Li, Zhenchun; Liu, Dingjin; Yong, Peng

doi:10.1007/s12517-020-05417-4

Cited by 7 publications

(5 citation statements)

References 28 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…There are two schemes can be considered: First, we can use finite‐difference method to replace the hybrid finite‐difference/pseudo‐spectral method for forward modelling. Second, we can adopt dual‐variable grid forward modelling method to reduce computation cost, which can be realized by matching the size variability of grid spacing using grid spacing of different sizes in different regions and locally variable time steps (Huang et al., 2015; Wang et al., 2020).…”

Section: Discussionmentioning

confidence: 99%

Least‐squares reverse time migration based on an efficient pure qP‐wave equation

Huang

Mao

et al. 2023

Geophysical Prospecting

Self Cite

View full text Add to dashboard Cite

Seismic waves propagating in anisotropic earth media suffer from waveform distortion, which leads to degraded resolution in seismic data imaging if this adverse effect is not corrected. Additionally, the wavefields simulated by the traditional coupled pseudo‐acoustic wave equations contain undesired shear wave artefacts and are unstable when Thomsen's anisotropy parameter ε is less than δ. To address these issues, many pure qP‐wave equations are proposed to describe wave propagation in anisotropic media. However, these equations are either low precision or high accuracy but computationally expensive. Considering the trade‐off between the computation efficiency and simulation accuracy, we derive a pure qP‐wave equation from the exact dispersion formula in tilted transverse isotropic media. The newly derived equation has only two higher order partial derivatives, which can better balance the accuracy and efficiency than the previous pure qP‐wave equations. The least‐squares reverse time migration method can balance amplitudes, suppress migration artefacts and improve imaging resolution. Then, based on the newly derived pure qP‐wave equation, we derive the Born modelling operator and adjoint migration operator and develop an anisotropic least‐squares reverse time migration method. To achieve the accelerated convergence, we introduce a generalized minimum residual algorithm to implement the anisotropic least‐squares reverse time migration method. Numerical simulation results show that the proposed pure qP‐wave equation outperforms the previous equations in balancing the trade‐off between accuracy and efficiency. In addition, synthetic examples demonstrate the effectiveness and robustness of our proposed anisotropic least‐squares reverse time migration method in correcting the deviations due to anisotropic effects.

show abstract

Section: Discussionmentioning

confidence: 99%

Least‐squares reverse time migration based on an efficient pure qP‐wave equation

Huang

Mao

et al. 2023

Geophysical Prospecting

Self Cite

View full text Add to dashboard Cite

show abstract

“…where v min (z) denotes the minimum velocity along the depth axis z, f m is the dominant frequency of the source wavelet, and dz(z) denotes the optimal vertical grid spacing. We use a rectangular sampling method (Wang et al, 2020) to resample the initial migration model. Firstly, we set a small trial step from z 0 and increase it continuously to get the first grid point η 1 , where η 1 dz(η 1 ).…”

Section: Adaptive Sampling Strategymentioning

confidence: 99%

“…Li et al (2014) developed a dual-variable grid algorithm and applied it to RTM. Li et al (2017) introduced the idea of pseudo-time domain (Alkhalifah, 2003;Ma and Alkhalifah, 2013) into LSRTM, which improves imaging efficiency by reducing vertical grid points (Wang et al, 2020). Proposed an adaptive grid discretization strategy and applied it to 3D LSRTM.…”

Section: Introductionmentioning

confidence: 99%

Diffraction Extraction and Least-Squares Reverse Time Migration Imaging for the Fault-Karst Structure With Adaptive Sampling Strategy

et al. 2022

Self Cite

View full text Add to dashboard Cite

The ultra-deep fault-karst structure discovered in the Tarim Basin in Western China is a fractured-vuggy carbonate reservoir with great potential for development. The diffraction generated by fractures, small-scale caves and vugs is often used for reservoir identification and seismic interpretation. Since the diffraction is much weaker than the reflection, it is difficult to separate the diffraction from the full wavefield. We use a plane-wave destruction (PWD) filter to extract the diffraction from the full data. In order to obtain a high-accuracy and amplitude-preserving imaging profile, we use the least-squares reverse time migration (LSRTM) method to image the separated diffraction. The large amount of calculation is the most challenging problem of the LSRTM algorithm. In view of this, we develop an adaptive sampling strategy to improve computing efficiency and reduce memory requirement. We use a fault model, a vugs-fractures model, and a fault-karst model to demonstrate the effectiveness and practicability of the proposed method. The numerical examples show that the proposed method can enhance the imaging resolution of the fault-karst structure and save computing cost without losing accuracy. In addition, a test on field data processing demonstrates the advantages of our algorithm.

show abstract

“…Gu et al (2015) developed a modified imaging condition of excitation amplitude prestack RTM, which saves memory consumption on the premise of ensuring accuracy. Wang et al (2019); Wang et al (2020) combined the variable grid method with the GPU parallel strategy and introduced it into FWI and LSRTM, which effectively improves the efficiency of the computation.…”

Section: Introductionmentioning

confidence: 99%

Pseudo-depth domain reverse time migration in VTI medium based on GPU parallel strategy

Xinwen

Huang²,

Jiang³

et al. 2023

Front. Earth Sci.

View full text Add to dashboard Cite

As one of the main imaging methods, a lot of attention has been directed to reverse time migration anisotropic media. However, reverse time migration still faces cost and memory limitations, especially in imaging high frequencies. A coordinate transformation strategy in transversely isotropic media with a vertical axis of symmetry (VTI) was utilized to obtain high-resolution images useful for stratigraphic interpretation. The coordinate transformation converted traditional depth domain imaging into pseudo-depth domain imaging. In addition, this paper used the Graphics Processing Unit (GPU) parallel strategy and compared the imaging approach with the conventional acoustic isotropy method. By testing several models and field data, it concludes that Graphics Processing Unit parallel strategy can significantly improve the computational efficiency of imaging. The pseudo-depth domain method also has obvious advantages in computational time and memory consumption. Moreover, it is more reasonable when considering anisotropy for complex laminated layers.

show abstract

An efficient and accurate finite-difference operator using adaptively discretized grids and its application for 3D least-squares reverse-time migration

Cited by 7 publications

References 28 publications

Least‐squares reverse time migration based on an efficient pure qP‐wave equation

Least‐squares reverse time migration based on an efficient pure qP‐wave equation

Diffraction Extraction and Least-Squares Reverse Time Migration Imaging for the Fault-Karst Structure With Adaptive Sampling Strategy

Pseudo-depth domain reverse time migration in VTI medium based on GPU parallel strategy

Contact Info

Product

Resources

About