Services on Demand
Journal
Article
Indicators
-
Cited by SciELO
-
Access statistics
Related links
-
Similars in SciELO
Share
Geofísica internacional
On-line version ISSN 2954-436XPrint version ISSN 0016-7169
Geofís. Intl vol.51 n.2 Ciudad de México Apr./Jun. 2012
Original paper
Application of the WignerVille distribution to interpret ground penetrating radar anomalies
Martha Angélica Elizondo1*, René E. Chávez2, María Encarnación Cámara3 and Andrés Tejero4
1 Posgrado en Ciencias de la Tierra Universidad Nacional Autónoma de México Ciudad Universitaria Delegación Coyoacán, 04510 México D.F., México *Corresponding author: maes112@yahoo.com.mx.
2 Instituto de Geofísica, Universidad Nacional Autónoma de México Ciudad Universitaria Delegación Coyoacán, 04510 México D.F., México.
3 Escuela de Ingenieros Industriales Universidad Politécnica de Madrid Madrid, España.
4 Facultad de Ingeniería, Universidad Nacional Autónoma de México, Ciudad Universitaria, Delegación Coyoacán, 04510, México D.F., México.
Received: November 18, 2010.
Accepted: October 30, 2011.
Published on line: March 30, 2012.
Resumen
Con base en la Distribución de WignerVille (WVD) se realizó un análisis en tiempo y frecuencia de datos obtenidos con el Radar de Penetración Terrestre (GPR), basado en el estudio de la descomposición de la señal espectral. Se calcula una correlación entre la señal original y las componentes de tiempofrecuencia para obtener anomalías estructurales de la información contenida en el radargrama relacionándola con la geología disponible. En primer lugar se describe la aplicación de un ejemplo teórico constituido por lo que representaría un túnel (tubería). Se obtuvieron las firmas correspondientes en el dominio del tiempo y en el dominio de la frecuencia. Finalmente se analiza esta metodología en un sitio de prueba en la detección de un tambo enterrado donde son conocidas la geometría y su profundidad. Este especial sitio fue facilitado por la Universidad Nacional Autónoma de México, en los terrenos del Observatorio Magnético de Teoloyucan, Estado de México. Los resultados obtenidos son bastante alentadores, ya que la WVD es capaz de definir los rasgos morfológicos relacionados con el tambo y abre la posibilidad de localizar este tipo de estructuras.
Palabras clave: Radar de Penetración Terrestre (GPR), análisis tiempofrecuencia, distribución de WignerVille, proceso de señales.
Abstract
Time frequency analysis, based on signal spectral analysis decomposition, is performed by the WignerVille distribution (WVD) for Ground Penetrating Radar (GPR) data. The cross correlation between the original signal and the timefrequency components is obtained to get the structural anomalies available in the information provided by the GPR, related to the geology. We describe the application to a theoretical example representing a tunnel (pipe). Corresponding signatures are obtained in the timefrequency domain. An actual application at a test site is presented, where a drum of known geometry has been buried. This is a special test site at the National Autonomous University of México, Magnetic Observatory at Teoloyucan, State of Mexico. Depth and dimensions are well controlled, and the results are encouraging, since the WVD is capable to define morphological features related to the drum and opens the possibility to locate these structures.
Key words: Ground Penetration Radar (GPR); timefrequency analysis; WignerVille distribution; signal processing.
Introduction
GPR (Ground Penetrating Radar) constitutes an important methodology to reveal subsurface features. Shallow buried features such as fractures, mines, caves and anthropogenic structures (pipes, connecting lines, etc.) can be successfully detected. Traditionally the processing of GPR data may be done in the Time Domain (TD) or the Frequency Domain (FD) by employing the Fourier transform. However, it is not always possible to obtain a reasonable resolution from GPR data, due to variations in the relative dielectric permittivity (RDP). Therefore, GPR interpretation requires the application of TD or FD filters to separate important information from noise. Unfortunately, such filtering may be contaminated with spurious information. Also, the estimated events or anomalies may be welldefined in FD, but not in TD. A timefrequency analysis (TF) can reveal the time variation of the frequencies contained in a 1D GPR signal.
In the TF linear distributions as well as in the Wavelet transform, the signal is decomposed in the time frequency domain based on amplitude, over a family of functions derived from a function base to determine the scale on which depends the resolution (Vera, 2003; Auger et al., 1995).
In the case of discrete Wavelet transform, it corresponds to a dyadic sampling of the timefrequency plane in which the function family constitutes an orthonormal base with a wavelet well localized in time and frequency (Auger et al., 1995; Rosado, 2000; RiveraRecillas et al., 2005).
The WignerVille distribution is a timefrequency bilinear transform, which computes the distribution of signal energy in two variables: time and frequency, as energy density. Its spectrum depends on time, preserving mathematical properties such as the marginal and instantaneous frequency. It employs the analytic signal from a real signal and a parametric function of constant value (kernel = 1).Thus, the interference terms appearing in the Wigner distribution due to negative frequency components of the real signal can be suppressed, (Ramos, 1997).
The idea is to extract the time series profile that make up the radargram and define a function in TF domain, which shows features related to anomalous bodies or subsurface structures. In this case, the WignerVille Distribution (WVD) is applied to define these structural objectives (Lopera et al., 2008). It is proposed to obtain the frequency information contained in time, for a 1D GPR signal processed by the WVD, in order to identify an event such as a stratigraphic change, a cavity or fracture. First, a 1D theoretical model is presented to show the ability of the WVD, and then a real example is shown from a test site where the characteristics are known. The test site is located in the facilities of the National Autonomous University of México, Magnetic Observatory at Teoloyucan town, in the State of Mexico, approximately 60 km NE of Mexico City.
Methodology
The WVD of a discrete signal s(t) is given by (Cohen, 1989):
If a signal is sampled at an interval T, the maximum cutoff frequency where the signal is folded in the Fourier spectrum is fc = 1/2T. However, the WVD timefrequency spectrum is folded at fc/2. It means that the original signal must be sampled at an interval T*=T/2 to avoid aleasing effects. Therefore, the new cutoff frequency can be expressed as fc*=1/2T*=1/T=2fc.
For the theoretical analysis a scan synthetic radargram was obtained for a 1D stratigraphy land, where the electric field at the Earth's surface can be obtained. Following the method of Weng (1971), based on the fact that a plane wave enters the Earth's surface and the total field obtained on the surface is the incident field plus the reflected, considering a dissipative dielectric Nlayered earth model and a generalized reflexion coefficient given by Stokes (Bellman and Wing, 1975) and solved in a recursive form as:
Later in the frequency domain, the Ricker's wavelet, given by (Ricker, 1977), is used:
Where f = frequency
fc = central frequency of oscillation
Which is characterized by a slow construction and a decline of energy at the extremes of the function, in which the standard of the wavelet is smooth and continuous in its 1st and 2nd derivatives and can be optimized to meet the center frequency of a GPR antenna, which is designed to match the input impedance at the central feed point and the origin of the wavelet is used to force the electric field vector in the central feed point (Daniels, 2004; Cassidy and Murdie, 2000) Figure 1.
This is done through the Convolution Theorem radagram obtaining the synthetic trace in the frequency domain (LázaroMancilla and GómezTreviño, 1996).
Where
E (f) = the electric field in the frequency domain
S (f) = the Ricker pulse in the frequency domain
Finally, the inverse Fourier transform is computed to go into the time domain and produce a physical meaning to the estimated radagram (Annan and Chua, 1992; Díaz, 2003).
It is expected that reflections produced by each layer represent the double travel time done by the GPR signal, given by:
The pulse travel time along the layered medium are calculated independently, the speed of propagation of electromagnetic waves is obtained through the formula.
Where c = speed of light
εr = relative permitivity
h = layer thickness
The simple travel time of each layer is
And the depth of the objects known can be estimated following, Peniche (2008).
Finally, we calculated the WVD to obtain the timefrequency plane, which is represented in graphical form by his absolute value.
Analysis and results of the GPR traces
Table (1) presents the model used to generate the synthetic trace (Figure 2). The synthetic trace model (Figure 2a) is represented in a sampling window of 60 ns, that would be obtained by an antenna of 270 MHz with 512 samples per record.
First we calculated the electric field, and the synthetic trace. Then, the timefrequency representation for this trace is obtained. However, small features such as a layer change can be diminished by the direct wave and its reflection surface, so the direct wave and its reflection can be removed by applying a filter (a simple average moving autoregressive to eliminate the unwanted echoes). The final result displays only important reflection components (Lopera, 2007).
The position of high amplitudes indicate events as at the times t1 = 18 and t2 = 39 ns. It is interesting to note that the amplitudes are compacted into the centre of the plane in a frequency range from 100 to 300 MHz. These represent the energy spectrum of the signal, which has a maximum value of 270 MHz (Figure 2).
Finally, a practical example is presented. A GPR profile at the test site, obtainid using a monostatic antenna, whose offset is zero, in which no lateral effects have a significant influence in giving a 1D representation; in the experimental site a cylindrical drum is buried at depth of 1.0 m, with the following dimensions: diameter of 0.60 m. An area of 1.20 × 8.00 m was scanned. GPR data were analyzed following the methodology already described, employing the WVD. Each trace presents a similar behavior (Figure 3), providing the possibility to locate with better resolution the position at depth of the drum.
We have taken a trace of the GPR profile. Then, the timefrequency representation for this trace is obtained; the position of high amplitudes can indicate the location where the cylindrical drum was buried; this event is indicated at the time t1 = 18 ns. These plots show us different events and are represented in a timefrequency plot as high energy values. It is interesting to note that the amplitudes are compacted into the centre of the plane in a frequency range from 100 to 300 MHz. These represent the energy spectrum of the signal, which has a maximum value of 270 MHz (Figure 4).
Conclusions
WVD timefrequency representation may provide another view of the GPR signal that complements what other methods may achieve. We provide a more accurate representation of the evolution of a nonstationary signal, which enables a more complete analysis of the signal providing classification, modeling and pattern identification parameters.
Nonstationary properties of the GPR signal allowed the application of WVD based analysis. A timefrequency representation was obtained, which allowed visualization of the frequency range contributing to the signal in the time scale studied. It allows a better visualization of the position of high amplitude that we can indicate the location of events of interest.
The spectral components represented in the timefrequency plane and depth, can be correlated with the characteristic events of the temporal signals or indirectly allow the characterization of the information contained in the GPR signal in terms of depths.
The study objectives from wich the were to allow characterization of the information contained in the GPR signal, applying the WignerVille transform to nonstationary signals, to visualize the evolution over time of the signal in a frequency range, to help in the classification, modelling and identification of patterns with their local geological correlation.
One of the goals of developed GPR methodologies was to recognize with adequate resolution structures lying at depths of 5.00 m at the most. New numerical approximations should be applied in the future to deal with noise and to enhance useful GPR signals to detect anomalies in the subsoil.
Bibliography
Annan A.P., Chua L.T., 1992, Ground Penetrating Radar performance prediction. Ground Penetrating Radar. Geological Survey of Canada. Canada 513. [ Links ]
Auger F., Flandrin P., GonÇalves P., Lemoine O., 19951996, TimeFrequency Toolbox. CNRS (France). Rice University (USA). [ Links ]
Bellman R., Wing M.G., 1975, An introduction to invariant imbedding. Wiley. New York. pp. 250. [ Links ]
CassandeMottin E., Pai A., 2005, Discrete Time and Frequency WignerVille Distribution: Moyal's Formula and Alliasing. IEEE. Signal Processing Letters. Vol. 12, No. 7. [ Links ]
Cassidy N.J., Murdie R.E., 2000, The application of mathematical modeling in the interpretation of nearsurface groundpenetrating radar sections. Proc. 8th Int. Conf. on GroundPenetrating radar (GPR2000). Australia, 2000. pp. 842847. [ Links ]
Classen T., Mecklenbräuker W., 1980, The Wigner Distribution: A Tool for TimeFrequency Signal Analysis. Philps J. Res. Vol 35, No.6. [ Links ]
Cohen L., 1989, Timefrequency distributions. A review. Proc. IEEE, 77(7), 941980. [ Links ]
Daniels D.J., 2004, Ground Penetrating Radar. 2nd Edition. Institution of Electrical Engineers, London United Kingdom. pp 723. [ Links ]
Díaz M.O., 2003, 1D modeling data forward and reverse ground penetrating radar. MA Thesis. Universidad Nacional Autónoma de México. UNAM. [ Links ]
Flandrin P., 1984, Some Features of timefrequency representations of multicomponent signal. IEEE Int. Conf. Acoust., Speech, Signal Processing. [ Links ]
Kumar P., FoufolaGeorgiou E., 1994, Wavelet analysis in geophysics: an introduction. Geophysics Wavelet, 4, 143. [ Links ]
LázaroMancilla O., Gómez Treviño, E. 1996. Synthetic radagrams from electrical conductivivity and magnetic permeability variations. Journal of applied Geophysics, 34. p. 283290. [ Links ]
Lopera O, Milisavljevic N., Daniels D., Macq B., 2007, Timefrequency domain signature analysis of GPR data for landmine identification. IEEE Transactions on Geosciences and Remote Sensing, pp. 159162. [ Links ]
Quian S., Chen D., 1996, Joint TimeFrequency Analysis. Prentice Hall. [ Links ]
Peniche J.F., 2008, Shaped of GPR's response of a stratified land using the method WKBJ. Thesis. Universidad Nacional Autónoma de México. Facultad de Ingeniería. UNAM. Thesis. [ Links ]
Ramos J., 1997, Dtección of earphones micropotentials of high frequency. Doctoral thesis. Catalunya Technical University, Department of Electronic Engineering, to Barcelona, Spain. [ Links ]
Ricker N.H., 1977, Transient waves in viscoelastic media. Elsevier Scientific Pub. Co. Amsterdam and New York. pp. 275. [ Links ]
RiveraRecillas D.E., LozadaZumaeta M.M., RonquilloJarillo G., CamposEnríquez J.O., 2005, Multiresolution analysis applied to interpretation of seismic reflection data. Geofísica Internacional 44, 355368. [ Links ]
Roth P., van Genderen F., Verhaegen M., 2001, Analysis of the influence of mine and soil properties on features extracted from GPR Data. Detections and Remediation Technologies for Mines IV, 4,394, 428439. [ Links ]
Vera De Payer E., 2003, Joint analysis Time frequency: Quadratic Representations: In processing of Signs and Images. Theory and Applications. Ed. Buenos Aires. Center of processing and signs. UTN BsAs.,110161. [ Links ]
Weng C.C. 1971, Waves and fields in inhomogeneous half. University of Illinois UrbanaChampaign. IEEE Press. The Institute of Electrical and Electronics Engineers, Inc. New York. pp. 608. [ Links ]