Attenuation analysis of real GPR wavelets: The equivalent amplitude spectrum (EAS)

Abstract

Absorption of a Ground Penetrating Radar (GPR) pulse is a frequency dependent attenuation mechanism which causes a spectral shift on the dominant frequency of GPR data. Both energy variation of GPR amplitude spectrum and spectral shift were used for the estimation of Quality Factor (Q*) and subsequently the characterization of the subsurface material properties. The variation of the amplitude spectrum energy has been studied by Spectral Ratio (SR) method and the frequency shift by the estimation of the Frequency Centroid Shift (FCS) or the Frequency Peak Shift (FPS) methods. The FPS method is more automatic, less robust. This work aims to increase the robustness of the FPS method by fitting a part of the amplitude spectrum of GPR data with Ricker, Gaussian, Sigmoid–Gaussian or Ricker–Gaussian functions. These functions fit different parts of the spectrum of a GPR reference wavelet and the Equivalent Amplitude Spectrum (EAS) is selected, reproducing Q* values used in forward Q* modeling analysis. Then, only the peak frequencies and the time differences between the reference wavelet and the subsequent reflected wavelets are used to estimate Q*. As long as the EAS is estimated, it is used for Q* evaluation in all the GPR section, under the assumption that the selected reference wavelet is representative. De-phasing and constant phase shift, for obtaining symmetrical wavelets, proved useful in the sufficiency of the horizons picking. Synthetic, experimental and real GPR data were examined in order to demonstrate the effectiveness of the proposed methodology.