Abstract:Ground-penetrating radar (GPR) profiling is the primary tool to provide detail information of the internal structure and characteristics inside a sand dune in the desert area. However, with the severe elevation change in a short horizontal distance and on the rugged and complex surface topography of the sand dunes, getting clear imaging of the internal structure of sand dunes from the GPR profiling is still a big challenge. The classic imaging technique such as the Kirchhoff migration has been applied to process GPR data on sand dunes, with limited success. The reverse-time migration (RTM) technique is the most advanced imaging technique that can handles GPR data acquired on rugged surface with severe topographic relief to generate subsurface structural images with high fidelity and has been tried to process GPR data on sand dunes. The results are encouraging, and the imaging quality is significantly improved. The finite-difference time-domain (FDTD) method is the major numerical tool for forward and backward continuations of the wave field during the RTM process. There are two aspects for using FDTD in RTM of GPR data on sand dunes still need improvement: the numerical scattering caused by the staircase approximation of the ground surface by using gridding in the Cartesian coordinate, and the negligence of the possible anisotropy of the electromagnetic material properties due to the calcareous cementation bedding inside a sand dune. In this paper, we develop the RTM algorithm based on the staggered grid FDTD that handles the rugged topographic surface by using the curvilinear coordinate, and the possible anisotropic radar wave velocity of the sand dune media. We first demonstrate the equivalency of the nonuniform, isotropic medium and the uniform, anisotropic medium for justifying using the uniform, anisotropic velocity in the RTM by synthetic modeling. Next, we validate our approach of using the synthetic data with the comparison of using the Cartesian coordinate and the curvilinear coordinate in an isotropic medium. The results indicate that the RTM algorithm using the curvilinear coordinate can efficiently eliminate the adverse effect of the staircase approximated boundary of the topography surface. Finally, we processed the real GPR data collected on a sand dune in the Badain Jaran desert by using the curvilinear coordinate and the uniform, anisotropic velocity in FDTD forward and backward wave field continuation. Comparison of the results indicates that the RTM imaging using the boundary-conforming curvilinear coordinate and anisotropic velocity gains more coherent and higher resolution images for the calcareous cementation layers and the water table.

Reverse Time Migration Method of Linear Frequency Modulation Ground Penetrating Radar Based on gprMax

Revers Time Migration Method of Ground Penetrating Radar Based on Gprmax

Three-Dimensional Reverse Time Migration of Ground-Penetrating Radar Signals

Improving reverse time migration of ground-penetrating radar under zero-time imaging conditions using cross-correlation superposition projection

GPR Least-Squares Reverse Time Migration Based on the Improved Cross Correlation Window

Reverse-time migration of GPR data for imaging cavities behind a reinforced shield tunnel

Frequency Compensation Technology in Ground Penetrating Radar for Improving Resolution

3D prestack reverse time migration of ground penetrating radar data based on the normalized correlation imaging condition

Attenuation compensated reverse time migration method of ground penetrating radar signals

Prestack Least-Squares Reverse Time Migration With an Exact Adjoint Operator for Ground-Penetrating Radar

Application of Reverse Time Migration on GPR Data for Detecting Internal Structures in a Sand Dune

Reverse Time Migration for Extended Obstacles: Electromagnetic Waves

Three-dimensional reverse-time migration applied to a MIMO GPR system for subsurface imaging

GPR detection localization of underground structures based on deep learning and reverse time migration

Imaging Shallow Fault Structures by Three-Dimensional Reverse Time Migration of Ground Penetration Radar Data

Reverse time migration for imaging periodic obstacles with electromagnetic plane wave

Reverse time migration of ground‐penetrating radar with full wavefield decomposition based on the Hilbert transform

3D Migration Depth Focus Velocity Analysis of Hand-Held Ground Penetrating Radar

Anisotropic Reverse-Time Migration of Ground-Penetrating Radar Data Collected on the Sand Dunes in the Badain Jaran Desert

An interpolating scaling functions method with low-storage five-stage fourth-order explicit Runge-Kutta schemes for 3D ground penetrating radar simulation

Elastic reverse time migration based on nested triangular mesh in 2D isotropic media