Transition operator approach to full waveform inversion of anisotropic elastic media
Identifikátory výsledku
Kód výsledku v IS VaVaI
<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F67985530%3A_____%2F20%3A00508959" target="_blank" >RIV/67985530:_____/20:00508959 - isvavai.cz</a>
Výsledek na webu
<a href="http://www.global-sci.com/intro/article_detail/cicp/16838.html" target="_blank" >http://www.global-sci.com/intro/article_detail/cicp/16838.html</a>
DOI - Digital Object Identifier
<a href="http://dx.doi.org/10.4208/cicp.OA-2018-0197" target="_blank" >10.4208/cicp.OA-2018-0197</a>
Alternativní jazyky
Jazyk výsledku
angličtina
Název v původním jazyce
Transition operator approach to full waveform inversion of anisotropic elastic media
Popis výsledku v původním jazyce
We generalize the existing distorted Born iterative T-matrix (DBIT) method to seismic full-waveform inversion (FWI) based on the scalar wave equation, so that it can be used for seismic FWI in arbitrary anisotropic elastic media with variable mass densities and elastic stiffness tensors. The elastodynamic wave equation for an arbitrary anisotropic heterogeneous medium is represented by an integral equation of the Lippmann-Schwinger type, with a 9-dimensional wave state (displacement-strain) vector. We solve this higher-dimensional Lippmann-Schwinger equation using a transition operator formalism used in quantum scattering theory. This allows for domain decomposition and novel variational estimates. The tensorial nonlinear inverse scattering problem is solved iteratively by using an expression for the Frechet derivatives of the scattered wavefield with respect to elastic stiffness tensor fields in terms of modified Green's functions and wave state vectors that are updated after each iteration. Since the generalized DBIT method is consistent with the Gauss-Newton method, it incorporates approximate Hessian information that is essential for the reduction of multi-parameter cross-talk effects. The DBIT method is implemented efficiently using a variant of the Levenberg-Marquard method, with adaptive selection of the regularization parameter after each iteration. In a series of numerical experiments based on synthetic waveform data for transversely isotropic media with vertical symmetry axes, we obtained a very good match between the true and inverted models when using the traditional Voigt parameterization. This suggests that the effects of crosstalk can be sufficiently reduced by the incorporation of Hessian information and the use of suitable regularization methods. Since the generalized DBIT method for FWI in anisotropic elastic media is naturally target-oriented, it may be particularly suitable for applications to seismic reservoir characterization and monitoring.
Název v anglickém jazyce
Transition operator approach to full waveform inversion of anisotropic elastic media
Popis výsledku anglicky
We generalize the existing distorted Born iterative T-matrix (DBIT) method to seismic full-waveform inversion (FWI) based on the scalar wave equation, so that it can be used for seismic FWI in arbitrary anisotropic elastic media with variable mass densities and elastic stiffness tensors. The elastodynamic wave equation for an arbitrary anisotropic heterogeneous medium is represented by an integral equation of the Lippmann-Schwinger type, with a 9-dimensional wave state (displacement-strain) vector. We solve this higher-dimensional Lippmann-Schwinger equation using a transition operator formalism used in quantum scattering theory. This allows for domain decomposition and novel variational estimates. The tensorial nonlinear inverse scattering problem is solved iteratively by using an expression for the Frechet derivatives of the scattered wavefield with respect to elastic stiffness tensor fields in terms of modified Green's functions and wave state vectors that are updated after each iteration. Since the generalized DBIT method is consistent with the Gauss-Newton method, it incorporates approximate Hessian information that is essential for the reduction of multi-parameter cross-talk effects. The DBIT method is implemented efficiently using a variant of the Levenberg-Marquard method, with adaptive selection of the regularization parameter after each iteration. In a series of numerical experiments based on synthetic waveform data for transversely isotropic media with vertical symmetry axes, we obtained a very good match between the true and inverted models when using the traditional Voigt parameterization. This suggests that the effects of crosstalk can be sufficiently reduced by the incorporation of Hessian information and the use of suitable regularization methods. Since the generalized DBIT method for FWI in anisotropic elastic media is naturally target-oriented, it may be particularly suitable for applications to seismic reservoir characterization and monitoring.
Klasifikace
Druh
J<sub>imp</sub> - Článek v periodiku v databázi Web of Science
CEP obor
—
OECD FORD obor
10507 - Volcanology
Návaznosti výsledku
Projekt
—
Návaznosti
I - Institucionalni podpora na dlouhodoby koncepcni rozvoj vyzkumne organizace
Ostatní
Rok uplatnění
2020
Kód důvěrnosti údajů
S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů
Údaje specifické pro druh výsledku
Název periodika
Communications in Computational Physics
ISSN
1815-2406
e-ISSN
—
Svazek periodika
28
Číslo periodika v rámci svazku
1
Stát vydavatele periodika
CN - Čínská lidová republika
Počet stran výsledku
31
Strana od-do
297-327
Kód UT WoS článku
000532318100015
EID výsledku v databázi Scopus
2-s2.0-85080025693