An efficient Helmholtz solver for acoustic transversely isotropic media

Handle URI:
http://hdl.handle.net/10754/626223
Title:
An efficient Helmholtz solver for acoustic transversely isotropic media
Authors:
Wu, Zedong; Alkhalifah, Tariq Ali ( 0000-0002-9363-9799 )
Abstract:
The acoustic approximation, even for anisotropic media, is widely used in current industry imaging and inversion algorithms mainly because P-waves constitute the majority of the energy recorded in seismic exploration. The resulting acoustic formulas tend to be simpler, resulting in more efficient implementations, and depend on less medium parameters. However, conventional solutions of the acoustic wave equation with higher-order derivatives suffer from S-wave artifacts. Thus, we propose to separate the quasi-P wave propagation in anisotropic media into the elliptic anisotropic operator (free of the artifacts) and the non-elliptic-anisotropic components, which form a pseudo-differential operator. We, then, develop a separable approximation of the dispersion relation of non-elliptic-anisotropic components, specifically for transversely isotropic (TI) media. Finally, we iteratively solve the simpler lower-order elliptical wave equation for a modified source function that includes the non-elliptical terms represented in the Fourier domain. A frequency domain Helmholtz formulation of the approach renders the iterative implementation efficient as the cost is dominated by the Lower-Upper (LU) decomposition of the impedance matrix for the simpler elliptical anisotropic model. Also, the resulting wavefield is free of S-wave artifacts and has balanced amplitude. Numerical examples show that the method is reasonably accurate and efficient.
KAUST Department:
Physical Sciences and Engineering (PSE) Division
Citation:
Wu Z, Alkhalifah T (2017) An efficient Helmholtz solver for acoustic transversely isotropic media. GEOPHYSICS: 1–20. Available: http://dx.doi.org/10.1190/geo2017-0618.1.
Publisher:
Society of Exploration Geophysicists
Journal:
GEOPHYSICS
Issue Date:
11-Nov-2017
DOI:
10.1190/geo2017-0618.1
Type:
Article
ISSN:
0016-8033; 1942-2156
Sponsors:
We thank KAUST for its support and the SWAG group for the collaborative environment. Especially, we thank Zhendong Zhang for useful discussions. We also thank Hemang Shah, Faqi Liu, Scott Morton, Hess Corporation and BP Exploration Operation for providing the benchmark model. The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST). For computer time, this research used the resources of the Supercomputing Laboratory at King Abdullah University of Science and Technology (KAUST) in Thuwal, Saudi Arabia. We also thank the associate editor Dimitri Komatitsch, assistant editor Arthur Cheng, Jiubing Cheng and another anonymous reviewer for their fruitful suggestions and comments.
Additional Links:
https://library.seg.org/doi/10.1190/geo2017-0618.1
Appears in Collections:
Articles; Physical Sciences and Engineering (PSE) Division

Full metadata record

DC FieldValue Language
dc.contributor.authorWu, Zedongen
dc.contributor.authorAlkhalifah, Tariq Alien
dc.date.accessioned2017-11-29T11:13:54Z-
dc.date.available2017-11-29T11:13:54Z-
dc.date.issued2017-11-11en
dc.identifier.citationWu Z, Alkhalifah T (2017) An efficient Helmholtz solver for acoustic transversely isotropic media. GEOPHYSICS: 1–20. Available: http://dx.doi.org/10.1190/geo2017-0618.1.en
dc.identifier.issn0016-8033en
dc.identifier.issn1942-2156en
dc.identifier.doi10.1190/geo2017-0618.1en
dc.identifier.urihttp://hdl.handle.net/10754/626223-
dc.description.abstractThe acoustic approximation, even for anisotropic media, is widely used in current industry imaging and inversion algorithms mainly because P-waves constitute the majority of the energy recorded in seismic exploration. The resulting acoustic formulas tend to be simpler, resulting in more efficient implementations, and depend on less medium parameters. However, conventional solutions of the acoustic wave equation with higher-order derivatives suffer from S-wave artifacts. Thus, we propose to separate the quasi-P wave propagation in anisotropic media into the elliptic anisotropic operator (free of the artifacts) and the non-elliptic-anisotropic components, which form a pseudo-differential operator. We, then, develop a separable approximation of the dispersion relation of non-elliptic-anisotropic components, specifically for transversely isotropic (TI) media. Finally, we iteratively solve the simpler lower-order elliptical wave equation for a modified source function that includes the non-elliptical terms represented in the Fourier domain. A frequency domain Helmholtz formulation of the approach renders the iterative implementation efficient as the cost is dominated by the Lower-Upper (LU) decomposition of the impedance matrix for the simpler elliptical anisotropic model. Also, the resulting wavefield is free of S-wave artifacts and has balanced amplitude. Numerical examples show that the method is reasonably accurate and efficient.en
dc.description.sponsorshipWe thank KAUST for its support and the SWAG group for the collaborative environment. Especially, we thank Zhendong Zhang for useful discussions. We also thank Hemang Shah, Faqi Liu, Scott Morton, Hess Corporation and BP Exploration Operation for providing the benchmark model. The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST). For computer time, this research used the resources of the Supercomputing Laboratory at King Abdullah University of Science and Technology (KAUST) in Thuwal, Saudi Arabia. We also thank the associate editor Dimitri Komatitsch, assistant editor Arthur Cheng, Jiubing Cheng and another anonymous reviewer for their fruitful suggestions and comments.en
dc.publisherSociety of Exploration Geophysicistsen
dc.relation.urlhttps://library.seg.org/doi/10.1190/geo2017-0618.1en
dc.rightsArchived with thanks to GEOPHYSICSen
dc.titleAn efficient Helmholtz solver for acoustic transversely isotropic mediaen
dc.typeArticleen
dc.contributor.departmentPhysical Sciences and Engineering (PSE) Divisionen
dc.identifier.journalGEOPHYSICSen
dc.eprint.versionPost-printen
kaust.authorWu, Zedongen
kaust.authorAlkhalifah, Tariq Alien
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.