A Frequency-Domain Quasi-Newton-Based Biparameter Synchronous Imaging Scheme for Ground-Penetrating Radar With Applications in Full Waveform Inversion

Abstract

Full waveform inversion (FWI) of ground-penetrating radar (GPR) data is becoming a promising technique to facilitate the interpretation of surface-GPR data and the mapping of the subsurface. However, more general FWIs still require a sufficient amount of RAM memory, and it is difficult to produce an accurate and representative reconstruction result due to a large amount of the Hessian matrix calculations and singular value decomposition (SVD). In this article, we developed a novel full-waveform approach of the limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) algorithm for surface-GPR data that is based on a quasi-Newton framework and the total variation (TV) regularization. The proposed approach uses an L-BFGS algorithm and combines a scale-transformation regularization technique to mitigate the ill-posed problem of inversion, which can impose a biparameter preinformation constraint to ensure the stability of inversion, and adaptive regularization weights are applied to improve the convergence efficiency of inversion. To demonstrate the novelty and effectiveness of the proposed scheme, we tested our FWI algorithm using synthetic data and in-site field data. In the testing, we focus on analyzing the influence of different aspects of the FWI results, including different scale factors, regularization weights, inversion strategies, acquisition configurations, initial models, and the noisy data set. In particular, the FWI experiment is performed to demonstrate the applicability of the proposed algorithm. The results show that the proposed algorithm can effectively reconstruct the biparameter near the subsurface with high accuracy, which makes our approach very attractive for attribute analysis applications and makes the surface-GPR FWI commercially viable.