Acta Geophysica | 2021
Variable density acoustic RTM of VSP data based on the time–space domain LS-based SFD method
Abstract
Vertical seismic profiling (VSP) can provide more abundant seismic wavefield information and better seismic data with high resolution and high quality for the complex underground geological structures compared with surface seismic data. Reverse time migration (RTM) method possesses significant advantages for the accurate identification of complex geological structures, and it’s considered to be the most accurate imaging method at present. Therefore, we develop a variable density acoustic RTM method which is applicable for VSP data to enhance the recognition capability of complex geological structures, and we also discuss different aspects of this proposed imaging method. Firstly, to effectively improve the modeling precision of seismic wavefields, the wavefield extrapolation of our VSP RTM method is realized by using an optimal staggered-grid finite difference (SFD) method to solve the variable density acoustic wave equation, because this optimal SFD method uses the least square (LS) method to optimize the objective function established by the time–space domain dispersion relation to estimate its difference coefficients. In other words, the time–space domain LS-based SFD method has higher numerical simulation accuracy for seismic modeling. Secondly, to effectively reduce the boundary reflections and storage requirements of our VSP RTM method, we adopt the PML absorbing boundary and the effective boundary storage strategy in the process of wavefield extrapolation. Finally, to strengthen the quality and precision of VSP RTM results, the depth imaging profile of a shot is calculated by the normalized cross-correlation imaging condition of sources which can effectively eliminate the source effects on RTM results, and Laplace filtering is applied to eliminate the imaging noises in final RTM results effectively. The imaging results of different models show the effectiveness of our RTM method for VSP data, and it can more accurately identify the complex underground geological structures compared with the RTM method for conventional surface seismic data.