Elastic Reverse-Time Migration with Complex Topography

Migration is an important step in seismic data processing for oil and gas exploration. The accuracy of migration directly affects the accuracy of subsequent oil and gas reservoir characterization. Reverse-time migration is one of the most accurate migration methods at present. Multi-wave and multicomponent seismic data contain more P- and S-wave information. Making full use of multi-wave and multicomponent seismic data can offer more information about underground structure and lithology, as well as improve the accuracy of seismic exploration. Elastic reverse-time migration (ERTM) has no dip restriction and can be applied to image multi-wave and multicomponent seismic data in complex structural areas and some special lithology structures. However, the surface topography of complex regions has an influence on wavefield and seriously degrades the quality of ERTM’s migration results. We developed a new ERTM method to migrate multi-wave and multicomponent seismic data in the region with complex surface topography. We first fill the layers between the highest and lowest undulating surface with near-surface elastic parameters in a complex topography model to obtain a new model with a horizontal surface. This allows the finite difference (FD) method based on the regular rectangular grid to be used to numerically solve elastic wave equations in the model with complex topography. The decoupled wave equations are used to generate source P- and S-waves and receiver P- and S-waves to reduce crosstalk artefacts in ERTM. A topography-related filter is further used to remove the influence of surface topography on migration results. The scalar imaging condition is also applied to generate PP and PS migration images. Some numerical examples with different complex topographies demonstrate that our proposed ERTM method can remove the influence of complex topography on ERTM’s images and effectively generate high-quality ERTM images.