2d Frequency-Domain Elastic Full-Waveform Inversion Using Time-Domain Modeling and a Multistep-Length Gradient Approach

DSpace/Manakin Repository

2d Frequency-Domain Elastic Full-Waveform Inversion Using Time-Domain Modeling and a Multistep-Length Gradient Approach

Show simple item record

title 2d Frequency-Domain Elastic Full-Waveform Inversion Using Time-Domain Modeling and a Multistep-Length Gradient Approach
contributor.author Xu, Kun (UT Dallas)
contributor.author McMechan, George A. (UT Dallas)
description.abstract To decouple the parameters in elastic full-waveform inversion (FWI), we evaluated a new multistep-length gradient approach to assign individual weights separately for each parameter gradient and search for an optimal step length along the composite gradient direction. To perform wavefield extrapolations for the inversion, we used parallelized high-precision finite-element (FE) modeling in the time domain. The inversion was implemented in the frequency domain; the data were obtained at every subsurface grid point using the discrete Fourier transform at each time-domain extrapolation step. We also used frequency selection to reduce cycle skipping, time windowing to remove the artifacts associated with different source spatial patterns between the test and predicted data, and source wavelet estimation at the receivers over the full frequency spectrum by using a fast Fourier transform. In the inversion, the velocity and density re-constructions behaved differently; as a low-wavenumber tomography (for velocities) and as a high-wavenumber migration (for density). Because velocities and density were coupled to some extent, variations were usually underestimated (smoothed) for V_P and V_S and correspondingly overestimated (sharpened) for ρ. The impedances I_P and I_S from the products of the velocity and density results compensated for the under-or overestimations of their variations, so the recovered impedances were closer to the correct ones than V_P, V_S, and ρ were separately. Simultaneous reconstruction of V_P, V_S, and ρ was robust on the FE and finite-difference synthetic data (without surface waves) from the elastic Marmousi-2 model; satisfactory results are obtained for V_P, V_S, ρ, and the recovered I_P and I_S from their products. Convergence is fast, needing only a few tens of iterations, rather than a few hundreds of iterations that are typical in most other elastic FWI algorithms.
description.sponsorship "The research in this paper was supported by the sponsors of the UT-Dallas Geophysical Consortium. S. Operto (SEISCOPE) provided a 2D P-SV-wave fourth-order staggered-grid FD modeling algorithm to verify our elastic FWI. K. Xu wishes to express his gratitude to F. Gao (at Total) for constructive discussions on the FWI algorithm. Finally, we would like to thank T. Alkhalifah, M. Sen, and two anonymous reviewers for their fruitful comments on the paper. The computations were performed on Linux clusters in the center for Lithospheric Studies at UT Dallas and at the Texas Advanced Computing Center. This paper is contribution no. 1255 from the Department of Geosciences at the University of Texas at Dallas."
identifier.citation Xu, Kun, and George A. McMechan. 2014. "2D frequency-domain elastic full-waveform inversion using time-domain modeling and a multistep-length gradient approach." Geophysics 79(2), doi:10.1190/GEO2013-0134.1
identifier.issn 0016-8033
identifier.uri http://hdl.handle.net/10735.1/4105
identifier.uri http://dx.doi.org/10.1190/GEO2013-0134.1
identifier.volume 79
identifier.issue 2
identifier.startpage 41
subject Waveform inversion
subject Time-domain analysis
subject Frequency domain
subject Seismic waves--Speed
subject Density
date.issued 2014-02-17
publisher Soc Exploration Geophysicists
rights ©2014 Society of Exploration Geophysicists. All Rights Reserved.
language.iso en
source.journal Geophysics

Files in this item

Files Size Format View
NSM-FR-GAMcMechan-271114.14.pdf 5.435Mb PDF View/Open Article

This item appears in the following Collection(s)


Show simple item record