Full-Waveform Inversion of Time-Lapse Crosshole GPR Data Using Markov Chain Monte Carlo Method

Abstract

Crosshole ground-penetrating radar (GPR) is an important tool for a wide range of geoscientific and engineering investigations, and the Markov chain Monte Carlo (MCMC) method is a heuristic global optimization method that can be used to solve the inversion problem. In this paper, we use time-lapse GPR full-waveform data to invert the dielectric permittivity. An inversion based on the MCMC method does not rely on an accurate initial model and can introduce any complex prior information. Time-lapse ground-penetrating radar has great potential to monitor the properties of a subsurface. For the time-lapse inversion, we used the double difference method to invert the time-lapse target area accurately and full-waveform data. We propose a local sampling strategy taking advantage of the a priori information in the Monte Carlo method, which can sample only the target area with a sequential Gibbs sampler. This method reduces the calculation and improves the inversion accuracy of the target area. We have provided inversion results of the synthetic time-lapse waveform data that show that the proposed method significantly improves accuracy in the target area.