Next Article in Journal
Improving the Accuracy of Land Use and Land Cover Classification of Landsat Data in an Agricultural Watershed
Previous Article in Journal
Precision Detection of Dense Litchi Fruit in UAV Images Based on Improved YOLOv5 Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Processing GPR Surveys in Civil Engineering to Locate Buried Structures in Highly Conductive Subsoils

1
Mechanical and Mining Engineering Department and CEACTEMA, Higher Polytechnic School of Linares, Technological Scientific Campus, University of Jaen, 23700 Linares, Spain
2
Andalusian Institute of Geophysics, Campus de Cartuja, University of Granada, 18071 Granada, Spain
3
Geology Department and CEACTEMA, Higher Polytechnic School of Linares, Technological Scientific Campus, University of Jaen, 23700 Linares, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(16), 4019; https://doi.org/10.3390/rs15164019
Submission received: 25 July 2023 / Revised: 11 August 2023 / Accepted: 12 August 2023 / Published: 14 August 2023

Abstract

:
Many studies have illustrated the great benefit of ground-penetrating radar (GPR) in civil engineering. However, in some cases, this geophysical survey method does not produce the desired results due to the electromagnetic characteristics of the subsoil. This study presents the results obtained in two locations near Linares (southern Spain), evaluating the detection of structures buried in conductive host materials (0.02 S/m in site 1 and 0.015 S/m in site 2) characterized by strong signal attenuation. Accounting for the study depth, which was 1.5 m, a 500 MHz shielded GPR antenna was used at both sites. At the first site, a controlled experiment was planned, and it consisted of burying three linear elements. An iron pipe, a PVC pipe, and a series of precast blocks were buried at a depth of 0.5 m in a subsoil composed of highly conductive clayey facies. To eliminate additional multiples caused by other superficial structures and increasing the high-frequency content, the predictive deconvolution flow was applied. In the 3D processing, the cover surfaces technique was used. Once the acquired GPR signals was analyzed and the optimal processing flow established, a second site in which different infrastructures in a conductive host medium formed by marly facies was explored. The 2D flow and 3D processing applied in this work allows to detect and see the continuity of some structures not visible for the default processing.

Graphical Abstract

1. Introduction

1.1. Theoretical Overview

Ground-penetrating radar (GPR) is a geophysical prospecting technique that is based on the study of the propagation of electromagnetic waves (EM) through the subsoil. This wave aspect means that it is a high-resolution method that operates at “high frequencies (>10 MHz) or in a low-loss state when the product between the conductivity (σ) and the dielectric constant ( ε ) is less than 1 (loss factor P << 1, known as loss tangent)” [1]; thus, the velocity of the electromagnetic waves does not depend on the frequency but only on the electromagnetic response of the medium in which they travel [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33]. The GPR signals are affected by a series of characteristics intrinsic to the materials, such as the composition, texture, degree of saturation, and composition of the interstitial fluids [2,3,4,5]. These two aspects (resolution and response to small variations in the environment) make GPR exploration a widely used method in surface geoscience studies, including studies focused on surface geology, hydrology, mining, environment, forensics, and archaeological or heritage studies [5,6,7,8,9,10,11,12,13,14,15].
In civil engineering, this technique has been successfully applied to evaluate the state of pavements, roads, bridges, tunnels, the location of utilities buried in the substrate and hidden structures in buildings as well as in the characterization and diagnosis of construction materials or the detection of humidity [16,17,18,19,20,21,22,23].
An algorithm that automatically detects and segments small cracks in asphalt pavement at the pixel level has been described [24]. The research results represent a further step toward accurately detecting and characterizing internal vertical cracks in asphalt pavement. This application is practical for considering subsequent maintenance work due to the use of the pavement [25]. A novel YOLOv3 model identifying and localizing concealed cracks in asphalt pavements through GPR was proposed [26,27]. An overview from 2015 to the present of the states of research on employing the GPR technique in the civil engineering was provided [28], in which promising directions for future research were identified. More directly related to the line of this research, there are also numerous studies that have applied GPR for the detection of buried services, basically cables and pipes. The technique allows locating the position and depth of these services as well as locating possible leaks or illegal taking of supplies [17,29,30,31].
However, despite the many examples that demonstrate its effectiveness, there are situations, generally unpublished, where this method fails to produce the expected outcomes. Fundamentally, these situations are due to two causes: (i) when there is little contrast between the structure to be detected and the host medium or (ii) when the host material is conductive and/or has high magnetic susceptibility [2,5,22,32].
In an undulatory wave motion and in the first approximation [2,5,33,34,35], the propagation equation of the electric field ( E ) for a plane wave is obtained as follows:
E = f ( β ± v t ) e α β
where f is a multivariable function in which β represents the distance in the direction of propagation, v represents the velocity of propagation, and α represents the attenuation coefficient as the parameter describing the decay of the wave energy with the propagation distance through a host media. For low-loss and diamagnetic materials [33,34,35,36,37,38,39], when P << 1 and considering that the loss factor is equal to the ratio between the conductivity (σ) and the product of the angular frequency (ω) with the permittivity (ε), P = σ/ωε, the above parameters are expressed as follows:
β = r   ·   κ ^
υ = c ε r μ r 1 + P 2 + 1 2 1 2 ,   if   σ < < ω ε   and   μ r 1 ,   then   resulting   in   υ = c ε r
α = ω ε r μ r c 1 + P 2 1 2 1 2 ,   if   σ < < ω ε   and   μ r 1 ,   resulting   in   α = 1 2 σ μ ε
where r is the position vector, k ^ is the direction and magnitude of EM wave propagation, ε r = ε ε 0 is the relative dielectric constant or relative permittivity, c is the speed of light in a vacuum, and σ and μ represent the conductivity and magnetic susceptibility (also known as permeability) of the material, respectively; μr is the relative permeability, and α the attenuation coefficient.
Equation (3) describes the dependence of the EM wave propagation velocity on the relative dielectric constants of the materials, establishing an inverse dependence. This relationship implies variations of up to an order of magnitude, where the extreme values are located between the air (εr = 1) and water (εr = 81). Equation (4) describes the attenuation coefficient as the decay of the EM’s energy as a function of the constitutive parameters. However, in the attenuation of an EM signal, other attenuation phenomena, such as diffraction, scattering, and chromatic dispersion, should be incorporated [5,33,34,35]. In this framework, the penetration of the EM wave is parameterized by the skin depth (δ), which is an important concept in EM theory and is used in GPR to describe the characteristic distance at which the amplitude of an electromagnetic wave significantly decreases as it propagates through a conducting or semiconducting host media [33,34,35], as it is an energy attenuation factor of the EM study due to its inverse relationship with the attenuation coefficient [1,3,33]. Normally, the skin depth refers to the distance at which the amplitudes are reduced by a factor of 1/e, that is, 36.8% of the original signal, and can be obtained as follows:
δ = 1 α resulting   in : δ = 503 1 σ f       f o r   ω ε < < σ     ( f i r s t   c o n d i t i o n ) 0.0053 ε r σ         f o r   σ < < ω ε   ( s e c o n d   c o n d i t i o n )
From the above equation, there are two conditions for the skin depth. The first one is when the product between the irradiated frequency and the dielectric constant is much smaller than the conductivity of the host media. In this case, the skin depth depends on the frequency of the device used, as is the case in studies performed with inductive methods, such as frequency domain EM (FDEM) [1,3]. In contrast, for the second condition, the conductivity is less than the product between the irradiated frequency and the dielectric constant giving rise to an equation that is calculated independently of the frequency of the device used. This is the case for non-inductive methods, such as ground-penetrating radar [1,3]. Therefore, for conductive media the skin depth tends to be low, increasing the attenuation coefficient of the EM signal.
The objective of this study is to evaluate signal detection in host media with a low skin depth, leading to materials such as conductive media near the electromagnetic wave state boundary for GPR surveys.

1.2. Study Areas

Studies were carried out at two sites near Linares (southern Spain, Figure 1A). The first zone (site 1) was selected in the vicinity of the city’s university campus (Figure 1B), which has a very clayey substrate. Three linear elements were buried (a PVC pipe, an iron pipe, and precast blocks), all of them at a depth of 0.5 m. In this case, knowing the buried structure, its position, and its depth is intended to verify the ability of the technique to detect them. The second one (site 2) consists of marly subsoil where urbanization works were carried out years ago (Figure 1C), but they were abandoned before their completion. In this second case, the goal is to determine the technique’s ability to detect the buried utilities associated with this urbanization, as there is no information on the existing structures under the substrate.

2. Materials and Methods

2.1. Experimental Design and Data Acquisition

At the first field experimental site (university campus; Figure 1B), three parallel trenches approximately 0.5 m wide, 4 m long, and 0.7 m deep were constructed (Figure 2A). The top of these three buried elements (Figure 2A–C, respectively) should be detected at approximately 0.5 m depth, and the average thickness of the structures is approximately 0.2 m. A PVC pipe was placed in the first trench to the east, next to one of the pre-existent concrete slabs (Figure 2A). In the second trench, the precast blocks were installed (Figure 2B). Note that these bricks have voids inside them. Finally, in the third trench, the iron pipes were placed (Figure 2C).
In site 1 (university zone), a square of 11 m × 11 m was planned for detecting the three buried structures described above. They also incorporated part of two old, highly degraded concrete slabs. Figure 3A shows the spread GPR chart with the location of buried elements and the cemented portions.
In site 2 (old, abandoned urbanization; Figure 1C and Figure 2D), some utilities were buried, such as the infrastructure for water, gas, sewer, electrical, and telephone facilities, but this information is not available (Figure 2D).
The two surveys were carried out in September 2022 on very dry lithologies. In both, a shielded antenna of 500 MHz (e RAMAC/GPR system, Pro-Ex model, by MALA GEOSCIENCE) was used (Figure 2D). This shielded antenna is a bow-tie design with an approximate maximum penetration depth of 6 meters. It provides a radial resolution of about 0.04 meters in the main radiation lobe for an environment with typical soil conditions where high-conductivity layers are absent. To define the parameters used in this work, several preliminary GPR campaigns were carried out, varying the time window (TW), trace spacing, and profile spacing. This previous work made it possible to define the parameters that offered the best response. Table 1 summarizes the main acquisition parameter in the two sites, whose acquisition geometry is shown in Figure 3.

2.2. Equipment

The experiment was performed with a 500 MHz shielded antenna from MALA Geosciences. This antenna consists of a shielded bow-tie design modified by the manufactured to avoid signal noise from external sources, such as mobile telephones, FM broadcast, television, and radar system. The shielded antenna means that most of the energy is only transmitted in one direction. The shielded antenna comprises both transmitter and receiver antenna elements in one single housing. The antenna is calibrated by the manufacturer to generate emissions below the limits established by Federal Communications Commission (FCC15.509 and FCC15.209) [40,41] and American National Standard Institute (ANSI C63.4) [42]. In this instance, the manufacturer had previously conducted multiple tests under antennas with the same characteristics. Based on these tests, the manufacturer provided information on the radiation patterns [43,44]. This includes a 500 MHz-centered main emission lobe, a pulse emission of 10.40 ns, a 3dB bandwidth of 361 MHz with an approximate maximum penetration depth of 6 m, and a radial resolution of 5 cm in typical soil conditions where high-conductivity layers are absent. Figure 4A,B show an idealized radiation pattern of the 500 MHz antenna. The main lobe is centered at 500 MHz, and the side lobes are around 200 MHz and 700 MHz.
With the above information, the theoretical signal attenuation of the antenna in clay media at 0.5 m could be calculated through the attenuation coefficient (α), considering the 500 MHz central frequency (Equation (6)) [45,46,47], a dielectric constant (ε, permittivity) for clay media of ~7.0 × 10−11 F/m (approximated to the permittivity of free space), and a magnetic susceptibility (μ, permeability) of ~1.0 H/m (approximated to the permeability of free space). The skin depth (δ) for clay is generally very low and can be considered negligible when calculated, such as tan (δ). Finally, the theoretical value of the attenuation coefficient (α) at the average frequency at 0.5 m is 0.026 Np/m, which is approximately equal to 0.226 dB/m.
α = ω μ ε tan δ
Therefore, despite the fact that Equation (6) includes frequency as a parameter for calculating the attenuation coefficient, it is more appropriate to use Equations (4) and (5), which are approximations made within the EM theory for the frequency range from 40 MHz to 5 GHz, which refer to light speeds and low-loss state. As can be inferred from these equations, regardless of the fact that the radiated frequency is a crucial factor in depth penetration and resolution, the signal losses (like the attenuation coefficient) do not depend on the frequency but rather on the properties of the host media [1,2,3,4,5,33].
It is important to notice that attenuation coefficient and losses are related concepts in GPR systems, but they represent different aspects of signal propagation and signal strength reduction in the subsurface medium. On the one hand, the attenuation coefficient (α) in GPR refers to the rate at which the electromagnetic wave’s energy decreases as it propagates through the subsurface medium. It is a measure of how much the signal weakens over a unit distance of propagation. The attenuation coefficient is usually expressed in units of Nepers per meter (Np/m) or decibels per meter (dB/m). The attenuation coefficient accounts for factors such as absorption, scattering, and reflection of the electromagnetic waves as they interact with the subsurface materials. Different materials have varying attenuation coefficients, and they can change with frequency and moisture content. For instance, dry clay soils have a theoretical α-value that spans from 0.06 to 6 dB/m, whereas wet clay produces values between 6 and 50 dB/m. Higher attenuation coefficients indicate stronger signal absorption and scattering, leading to faster signal weakening as the wave travels through the medium [3,33,48,49,50,51].
On the other hand, losses in GPR systems refer to the other mechanisms that cause a reduction of signal strength during its transmission, reception, and processing. These losses are caused by various factors, e.g., geometric expansion of wavefront, transmission and reflection losses, scattering losses, absorption losses, and antenna losses, among others [3,7,33,48,49,50,51].
Therefore, the attenuation coefficient quantifies the rate of signal weakening due to the subsurface medium’s properties, while losses encompass a broader set of factors that result in signal reduction during its propagation and measurement in the GPR system.

2.3. GPR Data Analysis of Site 1 to Stablish the 2D Processing Flow

Raw GPR data were processed using two codes: the Reflexw software, version 7.0 [38] and GPR-SLICE (https://www.gpr-survey.com/) that is accessed on 25 January 2023 [34] (see Table 2 and Figure 5). As a first starting point, in order to establish the processing flow, the effective skin depth (mean) in site 1 was performed utilizing Equation (5) for clay and marl media. These lithologies are characterized by a dielectric constant spanning from 10 to 40 and a host medium velocity ranging between 4.7 and 9.5 cm/ns, reaching a skin depth ranging from 1.0 to 1.67 meters, with attenuation coefficient between 5 and 9 dB/m when employing Equation (4). The values were originally provided in Nepers per meter and subsequently converted into decibels per meter (Figure 6D). In this instance, the velocity of the host media was calculated based on the known fact that there are three buried objects at 50 cm depth, which are detected at around 20 ns in the B-scans. Using the classic equation v = x t , a velocity of 5 cm/ns was calculated, which was then corrected by knowing that the air space between the antenna and soil travels at 4 ns. Hence, the effective time within the substrate is 16 ns, resulting in a medium velocity of 6.3 cm/ns, which corresponds to a dielectric constant of 22 (see Figure 6B). A hyperbolic analysis improved the calculation of the medium velocity and the dielectric constant, resulting in a dielectric constant of 20 and a medium velocity of 6.7 cm/ns. Based on the previously calculated values, a conductivity value of 0.02 S/m was obtained by interpolating between the limits of the dry and wet clay curves depicted in Figure 6C.
This relationship is contingent upon the conductivity of the substratum; in this case, 0.02 S/m gives an attenuation coefficient of 6 dB/m (Figure 6D), where media with higher conductivity (moist conditions) exhibit greater dielectric constants and slower wavefront velocities, whereas drier media display the opposite trend (Figure 6A,B).
Figure 7A shows an example of a raw B-scan where the linear structures are marked with arrows; weak hyperbolas indicate their poor detection. In the same Figure, the red points are the calculated skin depth for each scan; thus, this means an effective signal up to 20 ns on average. Another advisable control test to perform on the raw data is to analyze the band of frequencies where the different GPR signals are located. In this sense, Figure 7B shows the spectral frequency power normalized for same B-scan, where the interesting GPR signals are placed between 200 MHz to 500 MHz. The maximum spectral density in the spectrum is approximately 500 MHz and was calculated for the time window of 50 ns for site 1 and 40 ns for site 2. The maximum values of the spectral density indicate the frequency with the highest concentration. During the processing flow, these spectra are used as reference to configure the bandpass filters. Notice that the amplitudes are normalized between −1 and 1 in all radargrams, and this feature has no specific physical units, such as volts or amperes. Normalization in the context of B-scans is used to adjust the signal amplitudes and emphasize amplitude variations within the image. When amplitudes are normalized between −1 and 1, a relative range of amplitudes is established based on the maximum and minimum amplitudes present in the radargram data. In this case, a value of −1 represents the minimum detected amplitude, while a value of 1 represents the maximum detected amplitude. The normalization of amplitudes in a radargram allows for a more convenient and meaningful representation of amplitude variations in the image display. This facilitates the identification and analysis of subsurface features based on relative amplitude variations.
In conductive media with low detection, signal amplitude restitution must not only be compensated by the geometric expansion of the wavefront: additional gain corrections must also be applied due to the intrinsic absorption of the medium (Equations (4) and (5)). The above analyses indicate that the processing flow should be aimed at increasing the amplitude of the signals of interest and at “enriching” the 200 MHz–400 MHz frequency band, in which the reflections of buried elements are attributed during the experiment, as shown in Figure 7B. Furthermore, both processes must be carried out in such a way that the other reflections contained in the B-scans are affected the least. That is to say, a general increase in the amplitude of the data is not enough where the hyperbolas cannot be distinguished either, but this increase must be “selective”. Another important derived piece of information is that the processing must consider that the average value of the skin depth is located almost below the three structures. Therefore, the algorithms that are applied must be very conservative in the time window of 10 ns to 30 ns in order to avoid artifacts due to the low signal-to-noise ratio (S/N) that prevails in this interval.
Considering the previous aspects, Table 2 shows the processing flow applied to raw data for site 1, while the two images in Figure 8 correspond to the effect of the processing on the previous selected B-scan. The key step of this flow was the predictive deconvolution since with them, most of the multiples caused by other structures were eliminated, such as those located at the ends of the B-scan (path and concrete elements). In other words, the signals of interest were emphasized, minimizing the surrounding signals as much as possible. The data from site 2 (detailed in the following section) were analyzed using the same criteria.

2.4. Analysis of Site 1 Using 3D-GPR to Enhance the Coherence Signal

Figure 9 shows the depth slices for site 1 up to the mean depth of 20 ns, i.e., 0.65 m depth, considering a dielectric constant of 20, obtained from the hyperbola analysis for migration and confirmed by the velocities used to reach the buried features using the plots described in Figure 6. The successive images show the evolution of the buried elements, verifying that the three linear structures still present a low detection. It was also verified that the PVC pipe and the eastern end of the row of bricks were first detected at a depth of 40 cm, possibly caused by irregularities in the base of the trenches (±10 cm). Figure 9A emphasizes the buried elements, and the other reflections observed in the slices are related to concrete slabs with a structural mesh that masks the medium, as they are highly reflective. It can be seen that the buried blocks have very low reflectivity, and the pipes show hardly any amplitude contrast, as the signal was attenuated in the shallower centimeters.
Furthermore, the normalized B-scan amplitudes underwent a multiplication by a factor of 10,000 when generating slices from a three-dimensional volume. The processing software converts decimal numbers to integers in order to facilitate their manipulation. Consequently, the range of amplitudes displayed in the slices, spanning from −10,000 to 10,000, does not possess a corresponding physical unit.
In order to obtain a greater continuity of these three elements, the cover surfaces method [39] was applied considering three slices between 0.4 m and 0.65 m. 3D GPR Cover surfaces is an alternative technique to construct 3D images that improve the visualization of GPR results. Broadly speaking, this technique consists of considering the reflections whose amplitudes exceed a pre-set threshold value, in our case A = 1000, for each slice, which refers to the maximum amplitude to be covered by the algorithm. The threshold value was selected after several test with values between the 50% and 85% of the maximum amplitude. Thus, when the amplitude of a mesh-point in the initial slice is below the selected threshold, it will not be retained in the 3D cover surface. On the other hand, if the amplitude is much higher than the threshold value, it will be retained in the cover surface and may appear as a bright spot or anomaly normalized to the threshold value established.
Thereby, the data mesh (xi, yj) is traversed comparing the amplitudes Aij of each depth slice and considering the deepness of the first slice that meets the above condition. Therefore, an image is obtained that represents the depth at which the roof of the different structures is obtained. Figure 10 shows the result obtained for site 1, where the structure with the lowest detection was the brick row. From our perspective, it can be observed that during the burial process of the blocks, the side holes of the bricks were spontaneously occupied by sediments from the host medium. This resulted in the attainment of nearly identical dielectric constants between the brick–clay and clay materials, primarily due to the presence of low-impedance contrasts. Figure 10 is a result of the cover surface analysis and demonstrates an improvement in the detection of the buried pipes and scattered reflections near the blocks.

3. Results in Site 2 and Discussion

The processing flow and the cover surface algorithms were applied in site 2 (Figure 1C and Figure 3B). In this case, the relative dielectric constant was around 30, and the conductivity of the subsurface was 0.015 S/m (about 60 ohms*m), similar to site 1 (0.02 S/m, equal to 50 ohms*m), where the differences are based on the diverse geological environment (the lutitic subsoil was replaced, in this site, by marly facies), which in turn explains the skin-depth variations from one location to another. However, the GPR signals have similar features (Figure 11).
There is one aspect to comment on, which also occurs at site 1: when the ground surface is highly compacted, strong reflections can be produced in the interface of the subsequent layer. To reduce these multiples as much as possible, an additional algorithm was used, consisting of a horizontal filter that eliminates coherent signals longer than 3 m; the result is that the base of the compact layer is located around 15 ns on average, about 0.25 m deep with a dielectric constant of 15 (Figure 11C). These values give a skin-depth range from 1.5 to 2.1 m and attenuation coefficient from 3 to 7.5 dB/m (lower than previous stage). The conductivity value of 0.015 S/m gives an approach of 4.5 dB/m for the attenuation coefficient.
Two very clear hyperbolas at 3.5 m and 4.5 m from the origin of the GPR profile are shown in Figure 11A (profile position: X = 3 m). The first hyperbola was found at a depth of 0.3 m, and as seen in Figure 11, it can be correlated in all profiles. Its direction, which is parallel to the sidewalk, is detected in Figure 12. Given the strong dielectric contrast that it offers, it could correspond to a metallic pipe, which suggests a drinking water supply. The second of the hyperbolas (at 4.5 m from the origin in Figure 11) is at a similar depth (approximately 0.40 m). This hyperbola is correlated from x = 0 m to x = 6 m, and it could be a branch of the main metallic pipe towards the plots located in that section of sidewalk. It can be deduced from Figure 12. Therefore, it would be associated with the drinking water supply.
After additional processing was performed to remove the multiples associated with the compacted layer, two new hyperbolas were detected in the profile X = 3 (Figure 11C). The first, 10 m from the origin, is related to the trench made to install the central sewer line. According to the direct information taken in the outcrop (from the well that is accessed from one of the old manholes), the pipe should have been found at about 1.5–2 m depth. However, given its depth, it could not be detected in our study (only the excavation trench associated with it could be detected).
Between 14–16 m from the origin, a new hyperbola appeared at a depth of 0.30 m, and it showed good continuity in the rest of the profiles and in the slices made (Figure 12). This hyperbola could also correspond to a pipe, although in this case, given the low dielectric contrast, it may be non-metallic (probably polyethene) and destined for another type of supply (probably gas). On the surface, approximately 18 m from the origin, the profile passes over a small trench in the direction of the “x”-axis, where the remains of corrugated wire loom tubing can be observed. This trench may be the result of the looting, from the nearby old manhole, of the wiring to provide electrical or telephone supply. For this reason, the wiring in the GPR profile was not detected with continuity (Figure 12).
The slices depicted in Figure 12 were utilized as input files to perform the cover surfaces shown in Figure 13A. The algorithm systematically assesses the amplitude of each Aij position within the data mesh grids corresponding to each slice. In this context, Aij indicates the matrix notation representing the data values within each slice, with i and j denoting the rows and columns that define the position of the data amplitude within the mesh. During this evaluation process, the algorithm compares all the Aij positions across the input slices and records the depth at which the highest amplitude is observed. To ensure uniformity in the maximum amplitude values, a threshold of 1000 is established before executing the cover surfaces algorithm.
Consequently, the amplitudes of the four slices for each position are meticulously traversed and compared. As a result, a consolidated representation emerges in the form of a new single slice, possessing the key features exhibited by the preceding four slices. This cover surface slice accentuates the amplitude of signals emanating from the pipes or other detectable elements such as walls while discarding low amplitudes to maintain continuity of the reflected materials. The cover surface representation encapsulates a comprehensive overview of the detected bodies within the substrate, adhering to preselected thresholds. This enables a more precise interpretation of the findings, particularly in regions where signal attenuation is evident, as observed during the investigation process. Figure 13B represents a 3D image of site 2. The X/Y plane is represented at a depth of 0.50 m, showing all the pipes buried in the subsoil. The X/Z plane is at a distance of 18 meters from the origin, and it represents the trenches of each element in depth, as it is able to deduce the lateral limit of each trench.

4. Conclusions

In conductive media, the detection of structures is poor. This factor must be accounted for when defining a GPR survey aimed at detecting buried utilities in civil engineering. In this sense, in clayey substrate, the results may not be as expected.
We carried out a controlled experiment (site 1) to establish an optimal processing flow in order to recover the GPR signal as much as possible without creating artifacts. In this first site, consisting of a subsoil composed of highly conductive clayey facies, different structures were buried at a known position and depth in order to be used for calibration purposes during the experiment.
This optimal processing flow was divided into two stages: 2D processing of the profiles and 3D processing. The key step of the 2D flow was the predictive deconvolution since with them, most of the multiples caused by other structures were eliminated. It was possible to emphasize the signals of interest by reducing the surrounding signals as much as possible. In order to obtain a greater continuity of these three elements (3D processing), the cover surfaces method was applied. This technique involves an analysis of reflections with amplitudes exceeding a predetermined threshold value across a sequence of consecutive depth slices. This approach yields enhanced visual outcomes, thereby facilitating the interpretation process in environments characterized by signal attenuation, such as clay formations in the site 1. Once this processing flow was established, a GPR exploration was carried out in an old, abandoned urbanization, and the same flow was applied. The use of 2D flow defined in site 1 allowed visualizing some structures that were not detected by applying preliminary processing flow. On the other hand, the cover surfaces algorithms were also applied in site 2, where they made it possible to highlight the continuity of some of the structures.
Although the GPR technique has great limitations for studies of conductive subsoils (0.02 S/m in site 1 and 0.015 S/m in site 2), the obtained results show the improvement in detecting and defining the elements buried in the subsoil due to optimized processing determined by knowing the type of structure, its dimensions, and the position and depth at which it is located.

Author Contributions

J.R. conceived the initial ideas of the study; B.M., R.M. and J.R. participated in the field survey; C.A.-P. and J.R. processed and analyzed the data; R.M., C.A.-P., B.M., J.R. and M.C.H. contributed to the data interpretation; R.M., B.M., J.R. and M.C.H. contributed to the preparation of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been financed by the FEDER Andalucía R+D+i Project (reference 1380520), the Ministerio de Ciencia e Innovación Project (reference PID2021-123506OB-I00) and by the University of Jaen’s own funds.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank T. Teixido and J.A. Peña for the help provided during the course of the research and writing. In addition, the GuidelineGeo Company provided the information of the Mala 500 Mhz antenna, which was necessary to carry out the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Telford, W.M.; Geldart, L.P.; Sheriff, R.E. Applied Geophysics; Cambridge University Press: Cambridge, UK, 1990; 770p. [Google Scholar]
  2. Annan, A.P. Ground Penetrating Radar: Principles, Procedures & Applications; Sensors & Software Incorporated: Mississauga, QN, Canada, 2003; 278p. [Google Scholar]
  3. Daniels, D.J. Ground penetrating radar. In The Institution of Electrical Engineers, 2 ed.; Engineers: London, UK, 2004; 726p, Available online: http://refhub.elsevier.com/S0963-8695(17)30017-8/sbref8 (accessed on 15 January 2023).
  4. Gómez-Ortiz, D.; Martín-Crespo, T.; Martín-Velázquez, S.; Martínez-Pagán, P.; Higueras, H.; Manzano, M. Application of ground penetrating radar (GPR) to delineate clay layers in wetlands. A case study in the Soto Grande and Soto Chico watercourses, Doñana (SW Spain). J. Appl. Geophys. 2010, 72, 107–113. [Google Scholar] [CrossRef]
  5. Reynolds, J.M. An Introduction to Applied and Environmental Geophysics; John Wiley and Sons Ltd.: Chichester, UK, 2011; 796p. [Google Scholar]
  6. Neal, A. Ground-penetrating radar and its use in sedimentology: Principles, problems and progress. Earth-Sci. Rev. 2004, 66, 261–330. [Google Scholar] [CrossRef]
  7. Jol, H.M. Ground Penetrating Radar Theory and Applications; Elsevier Science: Amsterdam, The Netherlands, 2009; 544p. [Google Scholar]
  8. Pueyo-Anchuela, O.; Pocoví, A.; Casas-Sainz, A.M. Characterization of karst hazards from the perspective of the doline triangle using GPR—Examples from Central Ebro Basin (Spain). Eng. Geol. 2009, 108, 225–236. [Google Scholar] [CrossRef]
  9. Grangeia, C.; Ávila, P.; Matias, M.; Ferreira da Silva, E. Mine tailings integrated investigations: The case of Rio tailings (Panasqueira Mine, Central Portugal). Eng. Geol. 2011, 123, 359–372. [Google Scholar] [CrossRef]
  10. Martínez-Moreno, F.J.; Galindo-Zaldívar, J.; Pedrera, A.; Teixido, T.; Ruano, P.; Peña, J.A.; González-Castillo, L.; Ruiz-Constán, A.; López-Chicano, M.; Martín-Rosales, W. Integrated geophysical methods for studying the karst system of Gruta de las Maravillas (Aracena, Southwest Spain). J. Appl. Geophys. 2014, 107, 149–162. [Google Scholar] [CrossRef]
  11. Cortada, U.; Martínez, J.; Rey, J.; Hidalgo, M.C.; Sandoval, S. Assessment of tailings pond seals using geophysical and hydrochemical techniques. Eng. Geol. 2017, 223, 59–70. [Google Scholar] [CrossRef]
  12. Galdón, J.M.; Rey, J.; Martínez, J.; Hidalgo, M.C. Application of geophysical prospecting techniques to evaluate geological-mining heritage: The Sinapismo mine (La Carolina, Southern Spain). Eng. Geol. 2017, 218, 152–161. [Google Scholar] [CrossRef]
  13. Duarte, J.; Cunha, M.; Carvalho, J. Modeling and Comparison of Data Obtained by GPR, for Geological/Structural Analysis of a Carbonated Ornamental Rock Quarry—Blocometry Validation—Case Study in Valinho De Fátima, Portugal. In IOP Conference Series: Earth and Environmental Science, Proceedings of the 6th World Multidisciplinary Earth Sciences Symposium, Prague, Czech Republic, 7–11 September 2020; IOP Publishing: Philadelphia, PA, USA, 2020; p. 609. [Google Scholar] [CrossRef]
  14. Rey, J.; Mendoza, M.; Martínez, J.; Hidalgo, M.C.; Florez Rodríguez, C. Combining geophysical methods (DC, IP, TDEM and GPR) to characterize mining waste in the Linares-La Carolina district (southern Spain). J. Environ. Manag. 2022, 322, 104921. [Google Scholar] [CrossRef]
  15. Oliveira, R.J.; Caldeira, B.; Teixidó, T.; Borges, J.F.; Carneiro, A. Increasing the Lateral Resolution of 3D-GPR Datasets through 2D-FFT Interpolation with Application to a Case Study of the Roman Villa of Horta da Torre (Fronteira, Portugal). Remote Sens. 2022, 14, 4069. [Google Scholar] [CrossRef]
  16. Benedetto, A.; Pajewski, L. Civil Engineering Applications of Ground Penetrating Radar; Springer Transactions in Civil and Environmental Engineering; Springer International: New York, NY, USA, 2015. [Google Scholar]
  17. Prego, F.J.; Solla, M.; Puente, I.; Arias, P. Efficient GPR data acquisition to detect underground pipes. NDTE Int. 2017, 91, 22–31. [Google Scholar] [CrossRef]
  18. Wai-Lok Lai, W.; Dérobert, X.; Annan, P.A. Review of Ground Penetrating Radar application in civil engineering: A 30-year journey from Locating and Testing to Imaging and Diagnosis. NDTE Int. 2018, 96, 58–78. [Google Scholar] [CrossRef]
  19. Pajewski, L.; Fontul, S.; Solla, M. Ground-penetrating radar for the evaluation and monitoring of transport infrastructures. In Innovation in Near-Surface Geophysics. Instrumentation, Application, and Data Processing Methods; Elsevier: Amsterdam, The Netherlands, 2019; pp. 341–398. [Google Scholar]
  20. Rathod, H.; Debeck, S.; Gupta, R.; Chow, B. Applicability of GPR and a rebar detector to obtain rebar information of existing concrete structures. Case Stud. Constr. Mater. 2019, 11, e00240. [Google Scholar] [CrossRef]
  21. Solla, M.; Pérez-Gracia, V.; Fontul, S.A. Review of GPR Application on Transport Infrastructures: Troubleshooting and Best Practices. Remote Sens. 2021, 13, 672. [Google Scholar] [CrossRef]
  22. Neto, P.X.; Medeiros, W.E. A practical approach to correct attenuation effects in GPR data. J. Appl. Geophys. 2006, 59, 140–151. [Google Scholar] [CrossRef]
  23. Peng, Z.; Jing, G.; Wang, S.; Li, Y.; Guo, Y. Railway ballast layer inspection with different GPR antenas and frequencies. Transp. Geotech. 2022, 36, 100823. [Google Scholar] [CrossRef]
  24. Liu, Z.; Yeoh, J.K.W.; Gu, X.; Dong, Q.; Chen, Y.; Wu, W.; Wang, L.; Wang, D. Automatic pixel-level detection of vertical cracks in asphalt pavement based on GPR investigation and improved mask R-CNN. Autom. Constr. 2023, 146, 104689. [Google Scholar] [CrossRef]
  25. Liu, Z.; Gu, X. Performance evaluation of full-scale accelerated pavement using NDT and laboratory tests: A case study in Jiangsu, China. Case Stud. Constr. Mater. 2023, 18, e02083. [Google Scholar] [CrossRef]
  26. Liu, Z.; Gu, X.; Yang, H.; Wang, L.; Chen, Y.; Wang, D. Novel YOLOv3 model with structure and hyperparameter optimization for detection of pavement concealed cracks in GPR images. IEEE Trans. Intell. Transp. Syst. 2022, 23, 22258–22268. [Google Scholar] [CrossRef]
  27. Liu, Z.; Gu, X.; Chen, J.; Wang, D.; Chen, Y.; Wang, L. Automatic recognition of pavement cracks from combined GPR B-scan and C-scan images using multiscale feature fusion deep neural networks. Autom. Constr. 2023, 146, 104698. [Google Scholar] [CrossRef]
  28. Hou, F.; Rui, X.; Fan, X.; Zhang, H. Review of GPR Activities in Civil Infrastructures: Data Analysis and Applications. Remote Sens. 2023, 14, 5972. [Google Scholar] [CrossRef]
  29. Borgioli, G.; Capineri, L.; Falorni, P.; Matucci, S.; Windsor, C.G. The Detection of Buried Pipes from Time-of-Flight Radar Data. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2254–2266. [Google Scholar] [CrossRef]
  30. Rashed, M.A.; Al-Garni, M.A. On the application of GPR for locating underground utilities in urban areas. Arab. J. Geosci. 2012, 6, 3505–3511. [Google Scholar] [CrossRef]
  31. De Coster, A.; Pérez Medina, J.L.; Nottebaere, M.; Alkhalifeh, K.; Neyt, X.; Vanderdonckt, J.; Lambot, S. Towards an improvement of GPR-based detection of pipes and leaks in water distribution networks. J. Appl. Geophys. 2019, 162, 138–151. [Google Scholar] [CrossRef]
  32. Takahashi, K.; Preetz, H.; Igel, J. Soil properties and performance of landmine detection by metal detector and ground-penetrating radar—Soil characterisation and its verification by a field test. J. Appl. Geophys. 2011, 73, 368–377. [Google Scholar] [CrossRef]
  33. Cassidy, N.J. Electrical and Magnetic Properties of Rocks, Soils and Fluids. In Ground Penetrating Radar: Theory and Applications; Elsevier: Amsterdam, The Netherlands, 2009; pp. 41–72. [Google Scholar]
  34. Annan, A.P. Transmission dispersion and GPR. J. Environ. Eng. Geophys. 1996, 1, 125–136. [Google Scholar] [CrossRef]
  35. Annan, A.P. Ground-penetrating radar. In Near-Surface Geophysics; Society of Exploration Geophysicists: Tulsa, OK, USA, 2005; pp. 357–438. [Google Scholar]
  36. Ward, S.H.; Hohmann, G.W. Electromagnetic theory for geophysical applications. In Electromagnetic Methods in Applied Geophysics, Vol. 1, Theory; Nabighian, M.N., Ed.; Society of Exploration Geophysicists: Tulsa, OK, USA, 1998; pp. 131–312. [Google Scholar]
  37. Grimm, R.E.; Heggy, E.; Clifford, S.M.; Dinwiddlie, C.; McGinnis, R.; Farrell, D. Absorption and scattering in ground-penetrating radar: Analysis of the Bishop tuff. J. Geophys. Res. 2006, 11, E06S02. [Google Scholar] [CrossRef] [Green Version]
  38. Sandmeier, K.J. Program for the Processing of Seismic, Acoustic or Electromagnetic Reflection, Refraction and Transmission Data; REFLEXW Version 7.0; Software Manual; Sandmeier Geophysical Researcher: Karlsruhre, Germany, 2012; Available online: https://www.sandmeier-geo.de/reflexw.html (accessed on 15 January 2023).
  39. Peña, J.A.; Teixidó, T. Cover Surfaces as a New Technique for 3D GPR Image Enhancement. Archaeological Applications. 2013. Available online: http://hdl.handle.net/10481/22949 (accessed on 25 January 2023).
  40. Federal Communications Commission (FCC). Title 47—Telecommunication: Volume 1, Part 15, Section 15.509. 2021. Available online: https://www.ecfr.gov/current/title-47/chapter-I/subchapter-A/part-15/subpart-F/section-15.509 (accessed on 15 July 2023).
  41. Federal Communications Commission (FCC). Title 47—Telecommunication: Volume 1, Part 15, Section 15.209. 2021. Available online: https://www.ecfr.gov/current/title-47/chapter-I/subchapter-A/part-15/subpart-C/section-15.209 (accessed on 15 July 2023).
  42. American National Standards Institute (ANSI). American National Standard for Methods of Measurement of Radio-Noise Emissions from Low-Voltage Electrical and Electronic Equipment in the Range of 9 kHz to 40 GHz. In ANSI C63.4-2014 (Revision of ANSI C63.4-2009). 2014, pp. 1–170. Available online: https://ieeexplore.ieee.org/document/6840852 (accessed on 25 January 2023).
  43. Aerotech Telub: EMC test, 500MHz Shielded antenna. Ref: E014-TR020223. Mala Geoscience AB. 2002. Available online: https://fcc.report/FCC-ID/QLA500MHZ/290495.pdf (accessed on 15 July 2023).
  44. Pereira, M.; Rial, F.I.; Lorenzo, H.; Arias, P.; Novo, A. Analysis of Characteristics of Ground Penetrating Radar Antennas. In Natural Resources & Environmental Engineering, EUET Forestal; University of Vigo: Vigo, Spain, 2006. [Google Scholar]
  45. Jordan, E.C.; Balmain, K.G. Electromagnetic Waves and Radiating Systems; Prentice-Hall International, INC: Englewood Cliffs, NJ, USA, 1968; 753 p. [Google Scholar]
  46. Balanis, C.A. Antenna Theory: Analysis and Design; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  47. Cheng, D.K. Field and Wave Electromagnetics; Pearson Education India: Noida, India, 1989. [Google Scholar]
  48. Conyers, L.B.; Goodman, D. Ground-Penetrating Radar: An Introduction for Archaeologists; AltaMira Press: Lanham, MD, USA, 1997. [Google Scholar]
  49. Fisher, A.J.; Lowe, D.J. (Eds.) Ground Penetrating Radar in Sediments; Geological Society of London: London, UK, 2004. [Google Scholar]
  50. Matsuoka, T.; Takeda, H. Ground Penetrating Radar (GPR) Surveys for Railway Ballast Evaluation. J. Appl. Geophys. 2012, 78, 48–56. [Google Scholar]
  51. Sato, M.; Tsuchiya, N.; Yonezawa, C.; Sato, S. Attenuation Characteristics of GPR Signals in Pavement and Subsurface Structure. J. Appl. Geophys. 2004, 55, 117–130. [Google Scholar]
Figure 1. (A) Geographic location of the studied areas. (B) Aerial photograph showing the position of site 1, where the red square represents the controlled test, at the university campus. (C) Aerial photograph showing the rectangle placement of site 2 with industrial infrastructures.
Figure 1. (A) Geographic location of the studied areas. (B) Aerial photograph showing the position of site 1, where the red square represents the controlled test, at the university campus. (C) Aerial photograph showing the rectangle placement of site 2 with industrial infrastructures.
Remotesensing 15 04019 g001
Figure 2. (A) Views of the buried elements used in the tests carried out in site 1. PVC pipe was placed in the first trench to the east. (B) Second trench and the precast blocks installed. (C) Third trench, where the iron pipe was installed. (D) General view of the site 2: remains of old manholes and the direction of two perpendicular streets are observed. All data were acquired with a 500 MHz shielded antenna.
Figure 2. (A) Views of the buried elements used in the tests carried out in site 1. PVC pipe was placed in the first trench to the east. (B) Second trench and the precast blocks installed. (C) Third trench, where the iron pipe was installed. (D) General view of the site 2: remains of old manholes and the direction of two perpendicular streets are observed. All data were acquired with a 500 MHz shielded antenna.
Remotesensing 15 04019 g002
Figure 3. (A) Layout of the acquisition geometry in site 1 with the location of the controlled buried elements and the coverage of acquisition GPR profiles. (B) Layout of the acquisition geometry in site 2 with the location of the superficial elements and the GPR profiles coverage.
Figure 3. (A) Layout of the acquisition geometry in site 1 with the location of the controlled buried elements and the coverage of acquisition GPR profiles. (B) Layout of the acquisition geometry in site 2 with the location of the superficial elements and the GPR profiles coverage.
Remotesensing 15 04019 g003
Figure 4. (A,B) Idealized radiation pattern of the 500 MHz antenna constructed from the manufacturer emission test. The main lobe is centered at 500 MHz, and the side lobes are around 200 MHz and 700 MHz.
Figure 4. (A,B) Idealized radiation pattern of the 500 MHz antenna constructed from the manufacturer emission test. The main lobe is centered at 500 MHz, and the side lobes are around 200 MHz and 700 MHz.
Remotesensing 15 04019 g004
Figure 5. General GPR processing work-flow chart used in site 1 and site 2.
Figure 5. General GPR processing work-flow chart used in site 1 and site 2.
Remotesensing 15 04019 g005
Figure 6. (A) Plot of dielectric constant and wave travel time for velocities in wet and dry clays and marls media. (B) Plot of velocity as a function of dielectric constant using Equation (3). (C) Plot of skin depth as a function of dielectric constant for conductivities in clay media using Equation (5). (D) Absorption plot as a function of dielectric constant for conductivities in wet and dry clay media. The calculations for all these graphic parameters consider a centrally irradiated frequency of 500 MHz. The vertical black lines indicate the limits of the dielectric constants for clays and marls in dry and wet environments, whose skin depth and absorption values at the researched sites should fall within the shaded zones.
Figure 6. (A) Plot of dielectric constant and wave travel time for velocities in wet and dry clays and marls media. (B) Plot of velocity as a function of dielectric constant using Equation (3). (C) Plot of skin depth as a function of dielectric constant for conductivities in clay media using Equation (5). (D) Absorption plot as a function of dielectric constant for conductivities in wet and dry clay media. The calculations for all these graphic parameters consider a centrally irradiated frequency of 500 MHz. The vertical black lines indicate the limits of the dielectric constants for clays and marls in dry and wet environments, whose skin depth and absorption values at the researched sites should fall within the shaded zones.
Remotesensing 15 04019 g006
Figure 7. (A) Selected raw B-scan of the site 1 with the three buried structures (arrows). X-axis represents traces, and Y-axis represents depth in meters and nanoseconds. The image’s right side depicts the A-scan (trace) indicated by the light-yellow line on the B-scan (radargram). The corresponding low signals indicate the poor detection. The red line indicates the skin-depth points for any scan; the black rectangle indicates where the structures are mainly located. (B) The spectral power of the same B-scan. X-axis represents traces, and Y-axis represents the frequency in megahertz. On the right site is the color scale of the power-normalized spectrum density: higher-power spectral density values are around 500 MHz (vertical scale on the left side); lower-power spectral density values are around the 300 MHz and 200 MHz bandwidth.
Figure 7. (A) Selected raw B-scan of the site 1 with the three buried structures (arrows). X-axis represents traces, and Y-axis represents depth in meters and nanoseconds. The image’s right side depicts the A-scan (trace) indicated by the light-yellow line on the B-scan (radargram). The corresponding low signals indicate the poor detection. The red line indicates the skin-depth points for any scan; the black rectangle indicates where the structures are mainly located. (B) The spectral power of the same B-scan. X-axis represents traces, and Y-axis represents the frequency in megahertz. On the right site is the color scale of the power-normalized spectrum density: higher-power spectral density values are around 500 MHz (vertical scale on the left side); lower-power spectral density values are around the 300 MHz and 200 MHz bandwidth.
Remotesensing 15 04019 g007
Figure 8. (A) Resulting B-scan after the processing flow of Table 2 was applied. X-axis represents traces, and Y-axis represents depth in meters and nanoseconds. The three buried structures are marked with arrows. The image’s right side depicts the A-scan (trace) indicated by the light-yellow line on the B-scan (radargram). (B) Corresponding spectral normalized power. X-axis represents traces, and Y-axis represents the frequency in megahertz.
Figure 8. (A) Resulting B-scan after the processing flow of Table 2 was applied. X-axis represents traces, and Y-axis represents depth in meters and nanoseconds. The three buried structures are marked with arrows. The image’s right side depicts the A-scan (trace) indicated by the light-yellow line on the B-scan (radargram). (B) Corresponding spectral normalized power. X-axis represents traces, and Y-axis represents the frequency in megahertz.
Remotesensing 15 04019 g008
Figure 9. Slices at different depths to see the evolution of the buried elements in site 1. The top of the PVC pipe and some precast blocks were first detected at a depth of 0.4 m; the grid units are in meters. Rectangles highlight piperlines.
Figure 9. Slices at different depths to see the evolution of the buried elements in site 1. The top of the PVC pipe and some precast blocks were first detected at a depth of 0.4 m; the grid units are in meters. Rectangles highlight piperlines.
Remotesensing 15 04019 g009
Figure 10. Result of applying the cover surface technique between 0.4 and 0.65 m depth at site 1. Using this method, it was possible to increase the continuity of the two pipes. White rectangles highlight piperlines.
Figure 10. Result of applying the cover surface technique between 0.4 and 0.65 m depth at site 1. Using this method, it was possible to increase the continuity of the two pipes. White rectangles highlight piperlines.
Remotesensing 15 04019 g010
Figure 11. Selected raw B-scan of site 2 (X = 3) with two buried structures (arrows). (A) The red rectangle highlights strong reflections produced by compacted path. (B) The spectral power of the same B-scan. (C) Resulting B-scan after the flow processing.
Figure 11. Selected raw B-scan of site 2 (X = 3) with two buried structures (arrows). (A) The red rectangle highlights strong reflections produced by compacted path. (B) The spectral power of the same B-scan. (C) Resulting B-scan after the flow processing.
Remotesensing 15 04019 g011
Figure 12. Selected depth slices at different depths to see the constructive evolution of the buried elements of site 2; the grid units are in meters.
Figure 12. Selected depth slices at different depths to see the constructive evolution of the buried elements of site 2; the grid units are in meters.
Remotesensing 15 04019 g012
Figure 13. (A) Cover surface of site 2 up to the 0.65 m, coinciding with the mean of skin depth. (B) Three-dimensional view with the main pipelines placed at a depth of 0.5 m; the grid units are in meters.
Figure 13. (A) Cover surface of site 2 up to the 0.65 m, coinciding with the mean of skin depth. (B) Three-dimensional view with the main pipelines placed at a depth of 0.5 m; the grid units are in meters.
Remotesensing 15 04019 g013
Table 1. Main GPR acquisition parameters.
Table 1. Main GPR acquisition parameters.
Main GPR Acquisition Parameters
ParametersSite 1Site 2
Profile spacing0.5 m0.5 m
Number of profiles2341
Trace spacing 0.01 m 0.01 m
Time window30 ns40 ns
Samples trace 248312
Table 2. Planned 2D processing flow for site 1—controlled experiment.
Table 2. Planned 2D processing flow for site 1—controlled experiment.
StepsDescription
Subtract DC shiftTo eliminate the residual voltage
Static correctionTo adjust the delay time at time = 0 ns
Remove the background noiseTo eliminate the multiples proceed by coupled wave (air-surface)
Spiking predictive deconvolutionTo eliminate additional multiples caused by other superficial structures and increase the high-frequency content; Operator length of 21 ns, predication lag of 7 ns, pre-whitening 10%
Bandpass filter50 MHz–550 MHz to eliminate the high and low noise frequencies
Geometric amplitude correction To restitute the amplitude content for wave front propagation
Gain control (GC)To increase the amplitude between 10–30 ns where the buried structures are placed; a gain function was created by modifying the AGC
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mendoza, R.; Araque-Perez, C.; Marinho, B.; Rey, J.; Hidalgo, M.C. Processing GPR Surveys in Civil Engineering to Locate Buried Structures in Highly Conductive Subsoils. Remote Sens. 2023, 15, 4019. https://doi.org/10.3390/rs15164019

AMA Style

Mendoza R, Araque-Perez C, Marinho B, Rey J, Hidalgo MC. Processing GPR Surveys in Civil Engineering to Locate Buried Structures in Highly Conductive Subsoils. Remote Sensing. 2023; 15(16):4019. https://doi.org/10.3390/rs15164019

Chicago/Turabian Style

Mendoza, Rosendo, Carlos Araque-Perez, Bruna Marinho, Javier Rey, and Mari Carmen Hidalgo. 2023. "Processing GPR Surveys in Civil Engineering to Locate Buried Structures in Highly Conductive Subsoils" Remote Sensing 15, no. 16: 4019. https://doi.org/10.3390/rs15164019

APA Style

Mendoza, R., Araque-Perez, C., Marinho, B., Rey, J., & Hidalgo, M. C. (2023). Processing GPR Surveys in Civil Engineering to Locate Buried Structures in Highly Conductive Subsoils. Remote Sensing, 15(16), 4019. https://doi.org/10.3390/rs15164019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop