Matrix‐free crosshole elliptical‐anisotropy tomography: parametrization analysis and ground‐penetrating radar applications in carbonates. Issue 6 (28th August 2020)

• Record Type: Journal Article
• Authors: Salcedo, Marco; Garambois, Stéphane; Le Bouteiller, Philippe; Li, Yubing; Sénéchal, Guy; Danquigny, Charles; Virieux, Jean
• Is Part Of: Near surface geophysics. Volume 18:Issue 6(2020)
• Page Start: 697
• Page End: 712
• Publication Date: 2020-08-28
• Subjects: Anisotropy -- Multiparameter -- Tomography -- Traveltime -- Velocity; Earth (Planet) -- Surface -- Periodicals; Geophysics -- Technique -- Periodicals; Engineering geology -- Periodicals; Geophysics -- Periodicals; Planets -- Surfaces; Engineering geology; Geophysics -- Technique; Geophysics; Earth (Planet); Periodicals; 550
• DOI: 10.1002/nsg.12116
• Languages: English
• ISSNs: 1569-4445

Abstract

ABSTRACT: A novel traveltime tomographic approach is applied to anisotropic media, limited to 2D geometry at present. A general anisotropic Eikonal solver based on a discontinuous Galerkin method is combined with an efficient adjoint formulation for multiparameter least‐squares inversion. This new approach is tested considering synthetic crosshole ground‐penetrating radar data. The configuration of the ground‐penetrating radar survey is inspired by a real experiment done on layered carbonate media disturbed by the presence of a deep gallery, which induces a localized high‐electromagnetic contrast. This made it possible to define a well‐adapted general workflow in this context. We notably show that, under the elliptical anisotropic assumption, the parametrization based on vertical and horizontal velocities provides less biased results than those obtained by considering the vertical velocity and the relevant Thomsen parameter ε . The initial vertical and horizontal velocity models are identical and built from an isotropic inversion. The presence of the high‐contrast gallery generates a weak diffraction pattern, which is taken into account in our tomography approach. It also creates potential artefacts due to the model discretization, which are mitigated by a model regularization term within the definition of the misfit function. This general workflow is then applied to the real experiment dataset. The vertical and horizontal velocity images provide similar structures as those previously obtained by isotropic full waveform inversion, complemented by an image of a rather weak elliptical anisotropy.