KAUST DepartmentEarth Science and Engineering Program
Physical Science and Engineering (PSE) Division
Seismic Wave Analysis Group
Online Publication Date2018-08-27
Print Publication Date2018-08-27
Permanent link to this recordhttp://hdl.handle.net/10754/631148
MetadataShow full item record
AbstractNumerical solutions of the acoustic wave equation, especially in anisotropic media, is crucial to seismic modeling, imaging and inversion as it provides efficient, practical, and stable approximate representation of the medium. However, a clean implementation (free of shear wave artifacts and dispersion) of wave propagation, especially in anisotropic media, requires an integral operator, the direct evaluation of which is extremely expensive. Recently, the low-rank method was proposed to provide a good approximation to the integral operator utilizing Fourier transforms. Thus, we propose to split the integral operator into two terms. The first term provides a differential operator that approximates that can be approximated with a standard finite-difference method. We, then, apply the lowrank approximation on the residual term of the finite-difference operator. We implement the two terms in two complementing steps, in which the spectral step corrects for any errors admitted by the finite difference step. Even though we utilize finite-difference approximations, the resulting algorithm admits spectral accuracy. Also, through the finite difference step, the method can deal approximately with the free surface and absorbing boundary conditions in a straight forward manner. Numerical examples show that the method is of high accuracy and efficiency.
CitationWu Z, Alkhalifah T, Zhang Z (2018) A partial-low-rank method for solving acoustic wave equation. SEG Technical Program Expanded Abstracts 2018. Available: http://dx.doi.org/10.1190/segam2018-2991046.1.
SponsorsWe thank KAUST for its support and the SWAG group for the collaborative environment. We also thank BP for providing the benchmark dataset. The numerical examples are reproducible using the open source package Madagascar(http://www. ahay.org/wiki/Main_Page). The research reported in this publication is supported by funding from King Abdullah University of Science and Technology (KAUST). We thank the computing support from the Supercomputing Laboratory at King Abdullah University of Science and Technology (KAUST) in Thuwal, Saudi Arabia.
PublisherSociety of Exploration Geophysicists
Conference/Event name88th Society of Exploration Geophysicists International Exposition and Annual Meeting, SEG 2018