Laplace-domain wave-equation modeling and full waveform inversion in 3D isotropic elastic media
SCIE
SCOPUS
Cited 8 time in
WEB OF SCIENCE
Cited 9 time in
Scopus
-
Title
- Laplace-domain wave-equation modeling and full waveform inversion in 3D isotropic elastic media
-
Author(s)
- Son, Woohyun; Pyun, Sukjoon; Shin, Changsoo; Kim, Han-Joon
-
Alternative Author(s)
- 손우현; 김한준
-
Publication Year
- 2014-06
-
Abstract
- The 3D elastic problem has not been widely studied because of the computational burden. Over the past few years, 3D elastic full waveform inversion (FWI) techniques in the time and frequency domains have been proposed by some researchers based on developments in computer science. However, these techniques still have the non-uniqueness and high nonlinearity problems. In this paper, we propose a 3D elastic FWI algorithm in the Laplace domain that can mitigate these problems. To efficiently solve the impedance matrix, we adopt a first-order absorbing boundary condition that results in a symmetric system. A conjugate gradient (CG) solver can be used because the Laplace-domain wave equation is naturally positive definite. We apply the Jacobi preconditioner to increase the convergence speed. We identify the permissible range of Laplace damping constants through dispersion analysis and accuracy tests. We perform the Laplace-domain FWI based on a logarithmic objective function, and the inversion examples are designed for a land setting, which means that the source is vertically excited and multi-component data are considered. The inversion results indicate that the inversion that uses only the vertical component performs slightly better than the multi-component inversion. This unexpected result is obtained partly because we use a vertically polarized source. We analyze the residuals and Frechet derivatives for each component to examine the characteristics of the Laplace-domain multi-component FWI. The results indicate that the residuals and Frechet derivatives for the horizontal component have a singularity problem. The numerical examples demonstrate that the singularity problem is related to the directivity of the displacement and to taking the logarithm of Laplace-domain wave fields. To avoid this singularity problem, we use a simple method that excludes the data near the singular region. Although we can use either simultaneous or sequential strategies to invert the Laplace-domain data, we apply a simultaneous inversion strategy in this paper. Nonetheless, the numerical examples demonstrate that our inversion algorithm yields reasonable long-wavelength velocity structures. (C) 2014 Elsevier B.V. All rights reserved.
-
ISSN
- 0926-9851
-
URI
- https://sciwatch.kiost.ac.kr/handle/2020.kiost/2803
-
DOI
- 10.1016/j.jappgeo.2014.03.013
-
Bibliographic Citation
- JOURNAL OF APPLIED GEOPHYSICS, v.105, pp.120 - 132, 2014
-
Publisher
- ELSEVIER
-
Subject
- FINITE-DIFFERENCE; FREQUENCY-DOMAIN; TRAVEL-TIME; ACCURACY; SOLVER
-
Keywords
- 3D; Elastic media; Modeling; Full waveform inversion; Laplace domain
-
Type
- Article
-
Language
- English
-
Document Type
- Article
- Files in This Item:
-
There are no files associated with this item.
Items in ScienceWatch@KIOST are protected by copyright, with all rights reserved, unless otherwise indicated.