Fractal-Constrained Crosshole/Borehole-to-Surface Full-Waveform Inversion for Hydrogeological Applications Using Ground-Penetrating Radar

• Published in: IEEE transactions on geoscience and remote sensing Vol. 60; pp. 1 - 10
• Main Authors: Giannakis, Iraklis, Giannopoulos, Antonios, Warren, Craig, Sofroniou, Anastasia
• Format: Journal Article
• Language: English
• Published: New York IEEE 2022; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; Boreholes; Computational modeling; Computer applications; Computers; Dimensions; Distribution; Finite difference methods; Finite-difference time-domain (FDTD); Fractal geometry; Fractals; full-waveform inversion (FWI); Geology; Ground penetrating radar; ground-penetrating radar (GPR); Hydrogeology; Iterative methods; Optimization; Permittivity; principal component analysis (PCA); Principal components analysis; Radar; Soil; Soil properties; Soils; Time-domain analysis; Transmitters; Uptake; Waveforms
• ISSN: 0196-2892; 1558-0644
• DOI: 10.1109/TGRS.2021.3054173

Full-waveform inversion (FWI) is considered one of the most promising interpretation tools for hydrogeological applications using ground-penetrating radar. However, FWI has had limited practical uptake for several reasons: large computational requirements, an inability to reconstruct loss mechanisms of soil, and the need for a good initial starting model. We aim to address these issues via a novel FWI subject to a fractally correlated distribution of water. Initially, the dispersive properties of the soil are expressed as a function of the water fraction using a semiempirical model. This approach means that the permittivity, conductivity, and relaxation mechanisms are all correlated, and therefore, sensitivity problems between the permittivity and loss mechanisms no longer affect the performance of FWI. Subsequently, the distribution of the water fraction is constrained to follow a fractal geometry. Fractal-correlated noise is then compressed using a principal component analysis (PCA) in order to further reduce the number of the system's unknowns and accelerate FWI. PCA reduces the volume and dimensions of the optimization space, and thus, initialization is no longer necessary. Finally, a novel measurement configuration is suggested that uses superposition with all the individual measurements in order to reduce the number of forward models that need to be executed for every iteration of FWI. These enhancements substantially reduce the computational requirements of FWI and therefore eliminate the need for high-performance computers and time-consuming algorithms. The proposed scheme has been successfully tested with several numerical case-studies, which indicates the potential of this approach to become a commercially appealing interpretation tool for hydrogeology.