High resolution imaging of the unsaturated and saturated zones of a gravel aquifer using full-waveform inversion

Full-waveform inversion can significantly improve imaging results of cross-hole GPR data compared to conventional ray-based inversion schemes. Recently, a 2D full-waveform finite difference time domain (FDTD) approach was used to invert cross-hole GPR data measured in the saturated zone of a gravel aquifer. Due to water table refractions and reflections, the upper part of the aquifer was not reliably imaged. Here, we reconstruct the upper part of the aquifer by imaging both the saturated and unsaturated zones. Estimating one effective wavelet, as was done for the saturated zone inversion alone, is insufficient because the effective wavelet strongly depends on the location of both the transmitter and receiver antennas. Therefore, four different source wavelets were estimated for the different combinations of antennas placement in the two zones, and the full-waveform inversion algorithm adapted accordingly. This resulted in improved images of the aquifer. In general, the observed and the synthetic traces show a good correspondence in both shape and amplitude. For the transmitters in the unsaturated domain, the amplitude fit was not optimum and can probably be improved by adjusting the conductivity starting model.