Reverse-time prestack depth migration of GPR data from topography for amplitude reconstruction in complex environments

Abstract

With increased computational power, reverse-time prestack depth migration (RT-PSDM) has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar (GPR) applications. Complex topography alters the wavefield kinematics making for a challenging imaging problem. Model simulations show that topographic variation can substantially distort reflection amplitudes due to irregular wavefield spreading, attenuation anomalies due to irregular path lengths, and focusing and defocusing effects at the surface. The effects are magnified when the topographic variations are on the same order as the depth of investigation—a situation that is often encountered in GPR investigations. Here, I use a full wave-equation RT-PSDM algorithm to image GPR data in the presence of large topographic variability relative to the depth of investigation. The source and receiver wavefields are propagated directly from the topographic surface and this approach inherently corrects for irregular kinematics, spreading and attenuation. The results show that when GPR data are acquired in areas of extreme topography, RT-PSDM can accurately reconstruct reflector geometry as well as reflection amplitude.

Keywords

reverse-time prestack depth migration / ground-penetrating radar / topography / wavefield / reflector geometry / reflection amplitude

References

[1] Botelho M. A., Mufti I. R., Neto V. P. Multishot Prestack Depth Migration: An Application on Wide- Angle Reflection and Refraction GPR Data. SEG Technical Program Expanded Abstracts 1998, 1998, 1393-1396.

[2] Bradford J. H. Applying Reflection Tomography in the Postmigration Domain to Multifold Ground-Penetrating Radar Data. Geophysics, 2006, 71(1): K1-K8.

[3] Bradford J. H. Measuring Water Content Heterogeneity Using Multifold GPR with Reflection Tomography. Vadose Zone Journal, 2008, 7(1): 184-193.

[4] Bradford J. H. GPR Prestack Amplitude Recovery for Radiation Patterns Using a Full Wave-Equation, Reverse-Time Migration Algorithm. SEG Technical Program Expanded Abstracts 2012, 2012, 1-5.

[5] Chattopadhyay S., McMechan G. A. Imaging Conditions for Prestack Reverse-Time Migration. Geophysics, 2008, 73(3): S81-S89.

[6] Deng F., McMechan G. A. True-Amplitude Prestack Depth Migration. Geophysics, 2007, 72(3): S155-S166.

[7] Engheta N., Papas C. H., Elachi C. Radiation Patterns of Interfacial Dipole Antennas. Radio Science, 1982, 17(6): 1557-1566.

[8] Fisher E., McMechan G. A., Annan A. P., . Examples of Reverse-Time Migration of Single-Channel, Ground-Penetrating Radar Profiles. Geophysics, 1992, 57(4): 577-586.

[9] Lehmann F., Green A. G. Topographic Migration of Georadar Data: Implications for Acquisition and Processing. Geophysics, 2000, 65(3): 836-848.

[10] Leuschen C. J., Plumb R. G. A Matched-Filter-Based Reverse-Time Migration Algorithm for Ground-Penetrating Radar Data. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(5): 929-936.

[11] Sanada Y., Ashida Y. An Imaging Algorithm for GPR Data. Symposium on the Application of Geophysics to Engineering and Environmental Problems 1999, 1999, 565-573.

[12] Shragge J., Irving J., Artman B. Shot-Profile Migration of GPR Data. Proceedings of the 10th International Conference on Ground Penetrating Radar. Delft, 2004, 337-340.

[13] Topp G. C., Davis J. L., Annan A. P. Electromagnetic Determination of Soil Water Content: Measurements in Coaxial Transmission Lines. Water Resour. Res., 1980, 16: 574-582.

[14] Zhou D., Huang W. P., Xu C. L., . The Perfectly Matched Layer Boundary Condition for Scalar Finite-Difference Time-Domain Method. IEEE Photonics Technology Letters, 2001, 13(5): 454-456.

[15] Zhou H., Sato M., Liu H. J. Migration Velocity Analysis and Prestack Migration of Common-Transmitter GPR Data. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(1): 86-91.