3D elastic full‐waveform inversion of small‐scale heterogeneities in transmission geometry

Abstract

Three‐dimensional elastic full‐waveform inversion aims to reconstruct elastic material properties of 3D structures in the subsurface with high resolution. Here we present an implementation of 3D elastic full‐waveform inversion based on the adjoint‐state method. The code is optimized regarding runtime and storage costs by using a time‐frequency approach. The gradient is computed from monochromatic frequency‐domain particle‐velocity wavefields calculated with a time‐domain velocity‐stress finite‐difference scheme. The 3D full‐waveform inversion was applied to data of a complex random medium model, which resembles a realistic crystalline rock environment. We show synthetic inversion results of P‐wave and S‐wave velocities for two transmission geometries: (1) a 3D acquisition geometry with planes of sources and receivers and (2) a 2D geometry with two lines of sources and receivers, resembling a realistic two‐borehole geometry. The 3D inversion of data acquired with 3D source‐receiver geometry is capable to reconstruct differently sized 3D structures of shear and compressional velocities with resolution of about a wavelength. The 3D random medium data recorded with 2D acquisition geometry were inverted using 3D inversion and 2D full‐waveform inversion for comparison. The 2D inversion suffers from strong artefacts that are caused by 3D scattering. The multiparameter 3D inversion, by contrast, is capable to invert the 3D scattered waves and to reconstruct 3D structures up to about 1–2 wavelengths adjacent to the plane between sources and receivers. The resolution is lower compared to the 3D acquisition geometry result. Still, a 3D inversion of cross‐hole data can be beneficial compared to a 2D inversion in the presence of complex 3D small‐scale heterogeneities, as it is capable to resolve 3D structures next to the source‐receiver plane.

The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves

“…In most applications to field data, scattering of surface waves due to near-surface 3D structure might be significant. Unfortunately, 3D elastic FWI is still computationally too expensive for frequent repeat runs (Butzer et al, 2013) in a development environment. In this respect, we start with a 2D FWI that allows us to develop appropriate algorithms in an efficient way until sufficient computational power for routine full 3D FWI is available to all of us.…”

Acceleration strategies for elastic full waveform inversion workflows in 2D and 3D

“…Previous studies have already presented results of 3D elastic full waveform inversion: [6,10,17,18,34,42,43,46]. In order to validate our approach in 3D, the SEG/EAGE Overhrust model is used.…”

Pseudo-spectral full-waveform inversion

“…Метод сопряженных гра-диентов [Hestenes, Stiefel, 1952] существенно со-кращает число необходимых направленных спусков. В приложениях чаще всего применяется либо нели-нейная версия метода сопряженных градиентов [Polak, Ribiere, 1969;Butzer et al, 2013], либо алгоритм Брой-дена-Флетчера-Гольдфарба-Шанно с ограничен-ной памятью (в англоязычной литературе -метод L-BFGS) [Nocedal, 1980;Brossier et al, 2009].…”