Abstract:Reverse time migration (RTM) is used for subsurface imaging to handle complex velocity models including steeply dipping interfaces and dramatic lateral variations, and promises better imaging results compared to traditional migration methods such as Kirchhoff migration algorithm. RTM has been increasingly applied to seismic surveys for hydrocarbon resource explorations. Based on the similarity of kinematics and dynamics between electromagnetic waves and elastic waves, we develop a pre-stack RTM method and apply it to processing ground penetrating radar (GPR) data. The finite-difference time domain (FDTD) numerical method is used to simulate the electromagnetic wave propagation including forward and backward extrapolation. The cross-correlation imaging condition is used to obtain the final image. In order to provide a velocity model with relatively higher accuracy as the initial velocity model for RTM, we apply a full waveform inversion (FWI) in the time domain to estimate the subsurface velocity structure based on reflection radar data. For testing the effectiveness of the algorithm, we have constructed a complex geological model; and synthesized common-offset radar data and common-shot profile (CSP) radar reflection data. All data are migrated with the traditional Kirchhoff migration method and pre-stack RTM method separately. The migration results from pre-stack RTM show better coincidence with the true model. Furthermore, we have performed a physical experiment in a sandbox where a polyvinyl chloride (PVC) box is buried in the sand. The obtained common-offset radar data and common-shot radar data are migrated by using the Kirchhoff migration method and per-stack RTM algorithm separately. The per-stack RTM result shows that the RTM algorithm can obtain better imaging results. In the numerical experiments, velocity variation exists in the geological model due to the irregular interface between soil and sand. The imaging result from pre-stack RTM is better than that of the traditional Kirchhoff migration algorithm because pre-stack RTM can handle the model with horizontal velocity variation. In the RTM imaging, the shape and position of the interface and their anomalies match well with the true model, while the imaging of Kirchhoff migration cannot do so. As for the physical experiment, the Kirchhoff migration algorithm and pre-stack RTM have been applied to process the measured radar data. The imaging result from RTM show its advantage compared with the Kirchhoff migration method.

Pre-Stack Reverse-Time Migration Applied To Ground Penetrate Radar Synthetic Data

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

Application of Pre-Stack Reverse Time Migration Based on FWI Velocity Estimation to Ground Penetrating Radar Data

Reverse time migration applied to GPR data based on full wave inversion

Reverse Time Migration of Crosswell GPR Data Based on Wavefield Decomposition

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

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

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

Isotropic Elastic Reverse Time Migration

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

Elastic Reverse Time Migration Without Wavefield Decomposition

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

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

Frequency-Domain Reverse-Time Migration of Ground Penetrating Radar Based on Layered Medium Green's Functions

Two-Dimensional Reverse-Time Migration Applied to GPR with a 3-D-to-2-d Data Conversion

The Application Of Pre-Stack Migration To Sar-Gpr System For Imaging Of Obliquely Buried Landmine

Reverse time migration of the pseudo-space-domain elastic wave equation based on staggered-grid finite-difference scheme

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

Mesh-free Radial-Basis-function-generated Finite Differences and Their Application in Reverse Time Migration

Reverse Time Migration Based on the Pseudo-Space-Domain First-Order Velocity-Stress Acoustic Wave Equation