Premium
Matrix‐free crosshole elliptical‐anisotropy tomography: parametrization analysis and ground‐penetrating radar applications in carbonates
Author(s) -
Salcedo Marco,
Garambois Stéphane,
Le Bouteiller Philippe,
Li Yubing,
Sénéchal Guy,
Danquigny Charles,
Virieux Jean
Publication year - 2020
Publication title -
near surface geophysics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.639
H-Index - 39
eISSN - 1873-0604
pISSN - 1569-4445
DOI - 10.1002/nsg.12116
Subject(s) - geology , anisotropy , ground penetrating radar , geometry , isotropy , solver , discretization , tomography , radar , mathematical analysis , physics , optics , mathematics , computer science , mathematical optimization , telecommunications
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.