GPR clutter amplitude processing to detect shallow geological targets

The analysis of clutter in A-scans produced by energy randomly scattered in some specific geological structures, provides information about changes in the shallow sedimentary geology. The A-scans are composed by the coherent energy received from reflections on electromagnetic discontinuities and the incoherent waves from the scattering in small heterogeneities. The reflected waves are attenuated as consequence of absorption, geometrical spreading and losses due to reflections and scattering. Therefore, the amplitude of those waves diminishes and at certain two-way travel times becomes on the same magnitude as the background noise in the radargram, mainly produced by the scattering. The amplitude of the mean background noise is higher when the dispersion of the energy increases. Then, the mean amplitude measured in a properly selected time window is a measurement of the amount of the scattered energy and, therefore, a measurement of the increase of scatterers in the ground. This paper presents a simple processing that allows determining the Mean Amplitude of Incoherent Energy (MAEI) for each A-scan, which is represented in front of the position of the trace. This procedure is tested in a field study, in a city built on a sedimentary basin. The basin is crossed by a large number of hidden subterranean streams and paleochannels. The sedimentary structures due to alluvial deposits produce an amount of the random backscattering of the energy that is measured in a time window. The results are compared along the entire radar line, allowing the location of streams and paleochannels. Numerical models were also used in order to compare the synthetic traces with the field radargrams and to test the proposed processing methodology. The results underscore the amount of the MAEI over the streams and also the existence of a surrounding zone where the amplitude is increasing from the average value to the maximum obtained over the structure. Simulations show that this zone does not correspond to any particular geological change but is consequence of the path of the antenna that receives the scattered energy before arriving to the alluvial depositsPeer ReviewedPostprint (published version

Near-field imaging of locally perturbed periodic surfaces

This paper concerns the inverse scattering problem to reconstruct a locally perturbed periodic surface. Different from scattering problems with quasi-periodic incident fields and periodic surfaces, the scattered fields are no longer quasi-periodic. Thus the classical method for quasi-periodic scattering problems no longer works. In this paper, we apply a Floquet-Bloch transform based numerical method to reconstruct both the unknown periodic part and the unknown local perturbation from the near-field data. By transforming the original scattering problem into one defined in an infinite rectangle, the information of the surface is included in the coefficients. The numerical scheme contains two steps. The first step is to obtain an initial guess, i.e., the locations of both the periodic surfaces and the local perturbations, from a sampling method. The second step is to reconstruct the surface. As is proved in this paper, for some incident fields, the corresponding scattered fields carry little information of the perturbation. In this case, we use this scattered field to reconstruct the periodic surface. Then we could apply the data that carries more information of the perturbation to reconstruct the local perturbation. The Newton-CG method is applied to solve the associated optimization problems. Numerical examples are given at the end of this paper to show the efficiency of the numerical method