Wavefield compression for adjoint methods in full-waveform inversion

Document typeArticle
Defense date2016-09-13
PublisherSociety of Exploration Geophysicists
Rights accessOpen Access
All rights reserved. This work is protected by the corresponding intellectual and industrial
property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public
communication or transformation of this work are prohibited without permission of the copyright holder
ProjectHPC4E - HPC for Energy (EC-H2020-689772)
GEAGAM - Geophysical Exploration using Advanced GAlerkin Methods (EC-H2020-644202)
GEAGAM - Geophysical Exploration using Advanced GAlerkin Methods (EC-H2020-644202)
Abstract
Adjoint methods are a key ingredient of gradient-based full-waveform inversion schemes. While being conceptually elegant, they face the challenge of massive memory requirements caused by the opposite time directions of forward and adjoint simulations and the necessity to access both wavefields simultaneously for the computation of the sensitivity kernel. To overcome this bottleneck, we have developed lossy compression techniques that significantly reduce the memory requirements with only a small computational overhead. Our approach is tailored to adjoint methods and uses the fact that the computation of a sufficiently accurate sensitivity kernel does not require the fully resolved forward wavefield. The collection of methods comprises reinterpolation with a coarse temporal grid as well as adaptively chosen polynomial degree and floating-point precision to represent spatial snapshots of the forward wavefield on hierarchical grids. Furthermore, the first arrivals of adjoint waves are used to identify “shadow zones” that do not contribute to the sensitivity kernel. Numerical experiments show the high potential of this approach achieving an effective compression factor of three orders of magnitude with only a minor reduction in the rate of convergence. Moreover, it is computationally cheap and straightforward to integrate in finite-element wave propagation codes with possible extensions to finite-difference methods.
CitationBoehm, C. [et al.]. Wavefield compression for adjoint methods in full-waveform inversion. "Geophysics", 13 Setembre 2016, vol. 81, núm. 6, p. R385-R397.
ISSN0016-8033
Publisher versionhttp://library.seg.org/doi/abs/10.1190/geo2015-0653.1
Collections
Files | Description | Size | Format | View |
---|---|---|---|---|
Wavefield Compression for Adjoint Methods in.pdf | 6,807Mb | View/Open |