Search

Aiming at solving the clutter flooding problem in the traditional cross correlation reverse time migration (RTM) of ground penetrating radar (GPR), we proposed an improved RTM method based on multitraces cross correlation window (MCW) imaging condition. The main difference between the proposed method and the traditional direct stacking is that it can effectively enhance the effective signal while weakening the clutter by performing MC calculation on the single trace imaging results, avoiding the enhancement of both clutter interference and effective signal concurrently by direct stacking. Second, the window threshold is set according to the GPR observation accuracy, and the effective signal in the MC result is retained as the abnormal region window, while the imaging results in the nonabnormal region are discarded, so as to suppress the clutter and retain the abnormal region information. Numerical experiments show that compared with the traditional RTM and total variation (TV) denoising method with cross correlation imaging conditions, the MCW imaging condition can accurately locate abnormal region, suppress clutter interference, and have the advantage of no loss of effective information, which greatly improves the imaging quality. Finally, the proposed method is applied to the measured data to verify the practicability and effectiveness in practical engineering applications. [ABSTRACT FROM AUTHOR]

The full-waveform inversion (FWI) of ground-penetrating radar (GPR) data yields promise for quantitatively characterizing the parameters of the Earth’s shallow subsurface. However, conventional FWI is highly nonlinear and suffers from cycle skipping once the low-frequency data are missed or the initial model is poor. Furthermore, having limited prior knowledge of the subsurface in GPR measurements increases the ill-posedness of the inverse problem. Wavefield reconstruction inversion (WRI), which mitigates cycle skipping, extends the FWI search space by relaxing the wave equation constraint, reduces the nonlinearity, and is less sensitive to the initial model. In this article, we extend WRI to the 2-D frequency-domain imaging of on-ground GPR data. To improve the inversion stability and mitigate the ill-posedness, we utilize modified total variation (MTV) regularization to constrain the inverted models. With a simple numerical example, we first investigate the effects of the penalty parameter, initial models, and MTV regularization on WRI and further discuss the differences between WRI and FWI. Then, we analyze the sensitivity of the proposed approach to the frequency component and noise in a multiple-targets model. Furthermore, we assess our method with a complex synthetic example containing noise-contaminated data, showing that our proposed approach works efficiently even given noisy GPR data. Therefore, the reasonable combination of WRI and MTV regularization can improve the accuracy and efficiency of imaging for on-ground GPR data. This joint approach for the multiparameter quantitative reconstruction of GPR data ultimately exhibits good applicability and strong robustness and is worthy of promotion. [ABSTRACT FROM AUTHOR]

Ground-penetrating radar (GPR) full-waveform inversion (FWI) is an emerging near-surface high-resolution imaging technology that has been widely used in GPR 2-D imaging to better quantitatively evaluate the permittivity and conductivity of the subsurface. However, the 2-D imaging method limits the applicability and accuracy of FWI for the 3-D object cases. In this letter, we develop a 3-D FWI algorithm involving permittivity and conductivity for GPR in the frequency domain. The method employs the edge-based finite element method to solve the forward problem and the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method to solve the inverse problem, which avoids the calculation of the Hessian matrix. In particular, by using the parameter scaling and scaling factor, it is convenient to carry out the simultaneous inversion of permittivity and conductivity with high efficiency. In addition, a modified total variation regularization scheme is utilized for ensuring the stability of the inversion, as well as identifying the abnormal body boundary more effectively. As a demonstration, a synthetic 3-D model experiment is presented to test the proposed algorithm. The results show that the proposed inversion strategy can accurately reconstruct the permittivity structure together with the conductivity distribution from the 3-D on-ground GPR data. [ABSTRACT FROM AUTHOR]

To address problems, such as the computationally intensive inversion requirements, low inversion efficiency, and inadequate inversion accuracy caused by multiparameter crosstalk in a synchronous inversion, a new dynamic stochastic source encoding strategy combined with minmax-concave total variation (MCTV) regularization model constraints was proposed. This strategy avoids crosstalk noise between shots caused by the algorithm and greatly improves the inversion efficiency without affecting the inversion accuracy. By comparing a “cross”-shaped model with the multiparameter inversion results, we found that the MCTV regularization strategy boasts the best inversion effect. We further showed that dynamic stochastic source encoding can increase the inversion efficiency threefold by applying the 1994 British Petroleum (BP) migration international standard topography model and establishing a function to evaluate the most efficient inversion strategy from among seven options. Compared with the traditional stochastic source encoding strategies, dynamic stochastic source encoding was shown to better suppress crosstalk noise. The proposed strategy also presented a higher acceleration ratio; additionally, combining this strategy with MCTV regularization model constraints provided the clearest reconstructed image with the highest inversion precision and obtained the best evaluation score among the considered inversion strategies, albeit with a slight reduction in the total elapsed time-acceleration ratio. [ABSTRACT FROM AUTHOR]

Full-waveform inversion (FWI) in the time domain of ground-penetrating radar (GPR) data involves a vast number of calculations; thus, it requires a large amount of memory and is difficult to calculate on a personal computer (PC). In this paper, GPR data are analyzed with multiscale FWI using two parameters (permittivity and conductivity) based on total variation (TV) regularization, which is implemented on a PC using a graphics processing unit (GPU) parallel acceleration strategy. The inverse problem is considered to be a nonlinear optimization problem and is solved with limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) process, which is a quasi-Newton method. The gradient of the objective function is calculated using the adjoint-state method, and the finite-difference method is required to solve the forward problem many times. A multiscale serial inversion strategy is applied to optimize the inversion algorithm and to decompose the inversion problem into 2–3 frequency bands to search in the direction of the global minimum point instead of local minimums. Taking the complex model as an example, experiments are carried out to assess the parameter adjustment factor and regularization parameter. The appropriate parameter adjustment factor and regularization parameter can effectively guarantee the convergence speed and stability of dual-parameter inversion method and improve the accuracy of GPR data inversion. Finally, FWI of the noise-free and 25-dB signal-to-noise ratio (SNR) noise data of the overthrust model is performed. The results show that the multiscale and dual-parameter inversion method proposed in this paper can provide reliable constraints, has better adaptability to noisy data, and can reliably and accurately reconstruct the dielectric properties distribution of the subsurface. [ABSTRACT FROM AUTHOR]