Wavefield compression for adjoint methods in full-waveform inversion
Visualitza/Obre
Cita com:
hdl:2117/103434
Tipus de documentArticle
Data publicació2016-09-13
EditorSociety of Exploration Geophysicists
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
ProjecteHPC4E - 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.
CitacióBoehm, 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
Versió de l'editorhttp://library.seg.org/doi/abs/10.1190/geo2015-0653.1
Col·leccions
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
Wavefield Compression for Adjoint Methods in.pdf | 6,807Mb | Visualitza/Obre |