Acquisition and processing strategies for 3D georadar surveying a region characterized by rugged topography

Efficiently performing 3D ground-penetrating radar (GPR or georadar) surveys across rugged terrain and then processing the resultant data are challenging tasks. Conventional approaches using unconnected GPR and topographic surveying equipment are excessively time consuming for such environments, and special migration schemes may be required to produce meaningful images. We have collected GPR data across an unstable craggy mountain slope in the Swiss Alps using a novel acquisition system that records GPR and coincident coordinate information simultaneously. Undulating topography (dips of 8 ◦ to 16 ◦ ) and boulders with diameters up to about 2 m complicated the field campaign. After standard processing, the data were found to be plagued by time shifts associated with minor coordinate inaccuracies, uneven antenna-ground coupling, and numerous small gaps in data coverage. These problems were resolved by passing the data sequentially through an adaptive f -xy deconvolution routine and f -kx and f -ky filters. This filtering also reduced incoherent noise. Finally, the data were migrated using a 3D algorithm that accounted for the undulating topography. The nonmigrated and migrated images contained gently and moderately dipping reflections from lithological boundaries and actively opening fracture zones. A suite of prominent diffraction patterns was generated at a steeply dipping fracture zone that projected to the surface. Through this case history we introduce a general strategy for 3D GPR studies of topographically rugged land.

[1]  Rosemary Knight,et al.  Processing Ground Penetrating Radar Data To Improve Resolution Of Near-Surface Targets , 1993 .

[2]  Alan G. Green,et al.  3-D Georadar Mapping: Looking Into the Subsurface , 1996 .

[3]  Samuel D. Stearns,et al.  Signal processing algorithms using Fortran and C , 1992 .

[4]  Alan G. Green,et al.  Mapping the architecture of glaciofluvial sediments with three-dimensional georadar , 1995 .

[5]  S. Spitz.  Seismic trace interpolation in the F-X domain , 1991 .

[6]  H. Willenberg.  Multidisciplinary monitoring of progressive failure processes in brittle rock slopes – Concepts and system design , 2002 .

[7]  B. Widrow,et al.  Adaptive noise cancelling: Principles and applications , 1975 .

[8]  Alan G. Green,et al.  Semiautomated georadar data acquisition in three dimensions , 1999 .

[9]  Michael K. Chase,et al.  Random noise reduction by 3‐D spatial prediction filtering , 1992 .

[10]  Alan G. Green,et al.  Shallow geometry and displacements on the San Andreas Fault near Point Arena based on trenching and 3-D georadar surveying , 2002 .

[11]  Duncan Murray,et al.  Adaptive noise cancelling filter for 60 Hz , 1999 .

[12]  B. Widrow,et al.  Adaptive antenna systems , 1967 .

[13]  L. Canales.  Random Noise Reduction , 1984 .

[14]  Samuel D. Stearns,et al.  Signal processing algorithms , 1988 .

[15]  A. Green,et al.  Location and geometry of the Wellington Fault (New Zealand) defined by detailed three‐dimensional georadar data , 2004 .

[16]  Yanghua Wang.  Seismic trace interpolation in the f‐x‐y domain , 2002 .

[17]  E. Forte,et al.  2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data:a case history from an archaeological site , 1999 .

[18]  Alan G. Green,et al.  Topographic migration of georadar data: Implications for acquisition and processing , 2000 .