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.

Examples of pre-stack reverse-time migration applied to ground penetrate radar synthetic data

Elastic Reverse Time Migration Without Wavefield Decomposition

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

Isotropic Elastic Reverse Time Migration

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

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

Study on Influencing Factors of Pre-Stack Reverse Time Migration Imaging of Ground Penetrating Radar

Revers Time Migration Method of Ground Penetrating Radar Based on Gprmax

Wave Equation Migration for 3-D VSP Data

Implementation strategy of 3D reverse time migration on GPU/CPU clusters

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

An efficient least-squares reverse time migration in image domain

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

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

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

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

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

An efficient 3D wave-equation pre-stack time migration for high-density and wide-azimuth data

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

Acoustic prestack reverse time migration based on the symplectic method

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