Open AccessNon-orthogonal beam coordinate system wave propagation and reverse time migration
Author(s) -
Zhaolin Zhu,
Danping Cao,
Bangyu Wu,
Xingyao Yin,
Ying Wang
Publication year - 2019
Publication title -
journal of geophysics and engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.623
H-Index - 38
eISSN - 1742-2140
pISSN - 1742-2132
DOI - 10.1093/jge/gxz067
Subject(s) - coordinate system , seismic migration , paraxial approximation , coordinate space , grid , computation , wave equation , gaussian beam , ray tracing (physics) , finite difference , representation (politics) , beam (structure) , computer science , coordinate descent , algorithm , mathematical analysis , geometry , mathematics , optics , physics , politics , law , political science , geophysics
Grid size has a significant influence on the computation efficiency and accuracy of finite-difference seismic modeling and can change the workload of reverse time migration (RTM) remarkably. This paper proposes a non-orthogonal analytical coordinate system, beam coordinate system (BCS), for the solution of seismic wave propagation and RTM. Starting with an optical Gaussian beam width equation, we expand the representation on vertically variable velocity media, which is the most common scenario in seismic exploration. The BCS based on this representation can be used to implement an irregular-grid increment finite-difference that improves the efficiency of RTM. Based on the Laplacian expression in Riemannian space, we derive the wave equation in the BCS. The new coordinate system can generate an irregular grid with increment increasing vertically along depth. Through paraxial ray tracing, it can be extended to non-analytical beam coordinate system (NBCS). Experiments for the RTM on the Marmousi model with the BCS demonstrate that the proposed method improves the efficiency about 52% while maintaining good image quality.