# Application of the Cylindrical Wave Approach to  the Simulation of Buried Utilities

Authors: Fabrizio Frezza; Lara Pajewski; Cristina Ponti; Giuseppe Schettini
Journal: International Journal of Geophysics
 (2011)

---

## Abstract

The cylindrical wave approach is resumed. It is a full-wave technique for the solution of the two-dimensional plane-wave scattering problem by a set of perfectly conducting and dielectric cylinders buried in a dielectric half-space, or in a finite-thickness slab. This method can be applied for simulating scenarios with cylindrical inhomogeneities buried in the earth. New numerical results are presented, concerning the simulation of buried utilities. In particular, metallic pipes and air cavities are simulated.


---

## Body

## 1. Introduction

The Ground Penetrating Radar (GPR) [1, 2] is an established and routine method for the inspection of civil engineering structures. It can provide high-resolution images of the subsurface, typically from 0 to 10 m depth, through wide-band nonsinusoidal electromagnetic waves. The employed frequencies range from 10 MHz to 4 GHz. This technique is effective, rapid, nondestructive and non-invasive. It is employed for surveying of roads and highway pavements [3], bridges [4], tunnels [5], and for detecting subsurface cavities and voids [6]. It can also be used for utility sensing [7], for example, to map all the buried structures in a region, enabling rapid installation of a new plant with the minimum disruption and damage to the existing one. Gas, water, sewage, electricity, telephone and cable utilities can be localized. Moreover, the GPR can be used to perform detailed quality controls of reinforced concrete [8, 9]. An analysis of geological structures can be performed, for the mapping of soil, rock or fill layers in geotechnical investigations and for foundation design [10].An important factor, among those limiting the GPR surveying of buried pipes and cables, is the density of plant in urban areas: if a lot of underground infrastructures are present, it is difficult to interpret the measured data and to clearly image the scenario, but this is just the situation where clarity is more needed.The majority of buried plant is within 0.5 m to 2 m of the ground surface and it may have a wide variation in its size, may be metallic or nonmetallic, in close proximity to other plant. It may be buried in a wide range of soil types, involving large differences in electromagnetic propagation velocity and absorption. Therefore, obtaining adequate penetration of the emitted radiation together with good resolution is not straightforward, and some compromise has to be accepted.The block diagram of a radar system is shown in Figure1. A few nanosecond short impulse of electromagnetic energy is launched by a transmitting antenna. The antenna is mounted on a mobile trolley that is moved forward over the soil, at a very close distance from the ground surface. The energy scattered by the target is gathered by the receiving antenna, which is usually identical to the transmitting antenna and then processed by the receiver, to display the signal in a suitable form for the operator.Figure 1
Block diagram of a GPR system.Radars solve inverse problems, to estimate the electromagnetic properties of a target, or of a complex scenario, from field measurements. At present, different algorithms are employed in the post-processing of collected GPR data; most of them need a fast and accurate forward-scattering solver, to perform repeated evaluations of the scattered electromagnetic field due to known targets and to be used in combination with some optimization techniques.Almost all the objects sought in utility detection are long and thin, so in these applications two-dimensional scattering techniques are employed in the post-processing, being more effective and faster than three-dimensional ones.The direct two-dimensional electromagnetic-scattering problem of buried objects is,  therefore, a theme of great interest in GPR context, and it has been faced by several authors [1], both from a theoretical and a numerical point of view. Different methods have been developed for its solution, both in the frequency and time domains [11–13] (and references therein). Due to the complexity of the problem, many methods present limitations; for example, in [14] the obstacle size and the distance between buried scatterer and air-ground interface have to be much larger than the wavelength, the method proposed in [15] is valid only for small objects, no more than one cylinder is considered and only far-field results are calculated in [16]; some techniques can be applied only when the dielectric contrast between the obstacle and the hosting medium is low [17, 18]; some approaches suppose the ground to be lossless [11]. The finite-element method (FEM) [19], the finite-difference time-domain technique (FDTD) [20], and the method of moments (MoM) [21] can treat more general configurations and are often used.The cylindrical-wave approach (CWA) [11–13, 22, 23] is an efficient spectral-domain technique, developed for the rigorous solution of the two-dimensional electromagnetic forward-scattering problem by a finite set of perfectly conducting or dielectric targets, buried in a dielectric half-space or in a finite-thickness slab. In this method, the field scattered by underground objects is represented in terms of a superposition of cylindrical waves. Use is made of the plane-wave spectrum [24] to take into account the interaction of such waves with the planar interfaces between air and soil and between different layers in the ground. Suitable reflected [25] and transmitted [11] cylindrical functions are defined. Adaptive integration procedures of Gaussian type, together with acceleration algorithms, are employed for the numerical solution of the relevant spectral integrals [12, 26, 27]. All the multiple-reflection and -diffraction phenomena are taken into account.The CWA may deal with both TM and TE polarization fields. It can be applied for arbitrary values of permittivity, size, and position of the targets. Obstacles of general shape can be simulated, by means of a suitable set of small circular-section cylinders. Since the CWA is implemented in the frequency domain, dispersive soils can be modelled. The technique has been extended to study the scattering of an incident pulsed wave, with a rather general time-domain shape, a sampling of the incident-field spectrum can be performed, and the scattering problem can be solved in the frequency domain by using the CWA [13, 28]. The method is accurate and fast; therefore, it may be exploited in iterative algorithms for the solution of inverse problems.In Section2, we briefly resume the theoretical basis of the CWA. In Section 3, new results are presented, showing the effectiveness of the method for the sensing of cylindrical inhomogeneities buried in the earth. In particular, an electromagnetic simulation is performed of suitable scenarios in the context of civil engineering applications.


## 2. Theoretical Basis of the CWA

### 2.1. Perfectly Conducting Cylinders in a Dielectric Half-Space

The application of the CWA to the solution of the monochromatic plane wave-scattering problem by a finite set of perfectly conducting cylinders buried in a dielectric half-space is described in detail in [11]. In this subsection, the method is briefly resumed.Let us considerN perfectly conducting circular cylinders with possibly different radii, buried in a linear, isotropic, homogeneous, dielectric, lossless half-space, as schematized in Figure 2. Each cylinder is parallel to the y axis; the structure is assumed to be infinite along y direction. A monochromatic plane wave, with wavevector ki lying in the xz plane, impinges at an angle φi from medium 0 (a vacuum) on the planar interface with medium 1. We introduce a normalized reference frame (O,ξ,ζ), with coordinates ξ=k0x and ζ=k0z,k0=2π/λ0 being the vacuum wavenumber and λ0 the vacuum wavelength. In such reference frame, the qth cylinder has axis located in (χq,ηq), and its dimensionless radius is called αq=k0aq. The time dependence of the field is assumed to be e-iωt, where ω is the angular frequency. The solution to the scattering problem is carried out in terms of V(ξ,ζ), representing the y-component of the electric/magnetic field: V=Ey(ξ,ζ) for TM(y) polarization, and V=Hy(ξ,ζ) for TE(y) polarization.Figure 2
Geometry of the scattering problem, with cylinders buried in a dielectric half-space.In order to obtain a rigorous solution forV(ξ,ζ), the total field is expressed as the superposition of a set of terms, produced by the interaction between the incident field, the interface and the cylinders: the incident plane wave, the reflected field (due to the reflection in medium 0 of the incident plane wave by the interface), and the transmitted field (due to the transmission in medium 1 of the incident plane wave by the interface); the field scattered by the cylinders in medium 1, the scattered-reflected field (due to the reflection in medium 1 of the scattered field by the interface), the scattered-transmitted field (due to the transmission in medium 0 of the scattered field by the interface).The scattered field is expressed as the sum of the fields scattered by each buried obstacle. The field scattered by a cylinder is represented in terms of a superposition of cylindrical waves, by means of unknown coefficients, and use is made of the plane-wave spectrum of such waves [24] to treat their reflection and transmission by the interface. The presence of the planar interface, in fact, is taken into account by means of its complex plane-wave reflection and transmission coefficients. The final expression for the scattered field turns out to be, in a reference system centered on the qth cylinder axis (q=1,…,N),(1)Vs=V0∑m=-∞+∞Jm(n1ρq)eimθq∑h=1N∑l=-∞+∞ileilφtChl×[CWl-m(n1ξhq,n1ζhq)(1-δhq)+Hm(1)(n1ρq)Jm(n1ρq)δhqδml],
where V0 is the incident plane-wave amplitude, Jm is the first-kind Bessel function of integer order m,n1 is the refractive index of the medium hosting the scatterers, φt is the angle of the transmitted plane wave (according to Snell's law), chℓ are unknown coefficients, CWℓ(ξ,ζ)=Hℓ(1)(ρ)eiℓθ is the cylindrical wave, and Hℓ(1) is the first-kind Hankel function of integer order ℓ;  moreover, ξhq and ζhq are the distances, along ξ and ζ directions, of the axes of the hth and qth cylinders; finally, δ is the Kronecker symbol.To express the scattered-reflected and scattered-transmitted fields, the presence of the planar surface is taken into account by considering the reflection and transmission of each elementary plane wave constituting the Fourier spectrum of a cylindrical function. The scattered-reflected field is given by the sum of the fields scattered by each buried cylinder and reflected by the interface; by introducing the reflected cylindrical functionsRWℓ [11], the final expression is(2)Vsr⁡=V0∑m=-∞+∞Jm(n1ρq)eimθq∑h=1N∑l=-∞+∞Chlile-ilφt×RWm+l[-n1(χh+χq),n1(ηq-ηh)].
The scattered-transmitted field is given by the sum of the fields scattered by each buried cylinder and transmitted to medium 0. Introducing the transmitted cylindrical functions TWℓ [11], it can be found that(3)Vst=V0∑h=1Nchlile-ilφtTWl(ξh,ζh,χh),
Given the expressions of all the field terms, the boundary conditions on the cylinder surfaces are imposed; due to the hypothesis of perfectly conducting cylinders, the tangential component of the total electric field must vanish. A linear system for the unknown coefficients of the cylindrical wave expansions is obtained. Once such system is solved, the V(ξ,ζ) field is determined in any point of the space and for both polarizations. It is possible to derive all the other components of the electromagnetic field by using Maxwell's equations.


### 2.2. Dielectric Cylinders in a Dielectric Half-Space

When the buried scatterers are dielectric cylinders, the field is present also inside them. This problem is treated in detail in [12]. The expression of the field internal to the qth cylinder is given by an expansion in terms of first-kind Bessel functions with unknown coefficients dqm, as follows:(4)Vcq=V0∑m=-∞+∞dqmime-imφtJm(ncqρq)eimθq,
where ncq is the refractive index of the cylinder.On the cylinder surfaces, the boundary conditions have to be imposed. A linear system is thus obtained, for the unknown coefficients of the cylindrical wave expansions and for thedqm coefficients. Once the system has been solved, the total electromagnetic field is completely determined in any point of the space.


### 2.3. Cylinders in a Finite-Thickness Slab

When the scatterers are inside a finite-thickness slab, among two different half-spaces, additional reflected-transmitted, multiple-reflected and multiple-reflected-transmitted cylindrical functions are introduced, to consider the complicated interaction between the cylinders and the planar interfaces delimiting the slab. The case of perfectly conducting cylinders is solved in [22]; the solution to the scattering problem by dielectric cylinders in a slab is described in [23].It is important to observe that this extension of the method allows to simulate several interesting scenarios, as sketched in Figure3: objects in a soil layer, above a ground of different permittivity (a), objects inside a slab among two identical half-spaces, for example, cylinders in a wall (b), or else scatterers in a slab terminated on a perfectly conducting surface (c), in a layer of a stratified medium (d), or, finally, cylinders above (e) or in (f) a dielectric half-space.Cylinders in a finite-thickness slab, various scenarios.
(a)(b)(c)(d)(e)(f)


### 2.4. Pulsed-Wave Analysis

By using the CWA technique, it is possible to solve the scattering problem by buried perfectly conducting or dielectric cylinders of an incident pulsed plane wave, with a rather general time-domain shape. We have to perform a sampling of the incident-field spectrum and of the spectra of the various field terms involved (reflected, transmitted, scattered, scattered-reflected, scattered-transmitted, etc.). The scattering problem can be solved in the frequency domain for any sample through the CWA. Finally, by means of an inverse transform, the solution in the time domain can be calculated. This procedure is described in detail in [13, 28].


## 2.1. Perfectly Conducting Cylinders in a Dielectric Half-Space

The application of the CWA to the solution of the monochromatic plane wave-scattering problem by a finite set of perfectly conducting cylinders buried in a dielectric half-space is described in detail in [11]. In this subsection, the method is briefly resumed.Let us considerN perfectly conducting circular cylinders with possibly different radii, buried in a linear, isotropic, homogeneous, dielectric, lossless half-space, as schematized in Figure 2. Each cylinder is parallel to the y axis; the structure is assumed to be infinite along y direction. A monochromatic plane wave, with wavevector ki lying in the xz plane, impinges at an angle φi from medium 0 (a vacuum) on the planar interface with medium 1. We introduce a normalized reference frame (O,ξ,ζ), with coordinates ξ=k0x and ζ=k0z,k0=2π/λ0 being the vacuum wavenumber and λ0 the vacuum wavelength. In such reference frame, the qth cylinder has axis located in (χq,ηq), and its dimensionless radius is called αq=k0aq. The time dependence of the field is assumed to be e-iωt, where ω is the angular frequency. The solution to the scattering problem is carried out in terms of V(ξ,ζ), representing the y-component of the electric/magnetic field: V=Ey(ξ,ζ) for TM(y) polarization, and V=Hy(ξ,ζ) for TE(y) polarization.Figure 2
Geometry of the scattering problem, with cylinders buried in a dielectric half-space.In order to obtain a rigorous solution forV(ξ,ζ), the total field is expressed as the superposition of a set of terms, produced by the interaction between the incident field, the interface and the cylinders: the incident plane wave, the reflected field (due to the reflection in medium 0 of the incident plane wave by the interface), and the transmitted field (due to the transmission in medium 1 of the incident plane wave by the interface); the field scattered by the cylinders in medium 1, the scattered-reflected field (due to the reflection in medium 1 of the scattered field by the interface), the scattered-transmitted field (due to the transmission in medium 0 of the scattered field by the interface).The scattered field is expressed as the sum of the fields scattered by each buried obstacle. The field scattered by a cylinder is represented in terms of a superposition of cylindrical waves, by means of unknown coefficients, and use is made of the plane-wave spectrum of such waves [24] to treat their reflection and transmission by the interface. The presence of the planar interface, in fact, is taken into account by means of its complex plane-wave reflection and transmission coefficients. The final expression for the scattered field turns out to be, in a reference system centered on the qth cylinder axis (q=1,…,N),(1)Vs=V0∑m=-∞+∞Jm(n1ρq)eimθq∑h=1N∑l=-∞+∞ileilφtChl×[CWl-m(n1ξhq,n1ζhq)(1-δhq)+Hm(1)(n1ρq)Jm(n1ρq)δhqδml],
where V0 is the incident plane-wave amplitude, Jm is the first-kind Bessel function of integer order m,n1 is the refractive index of the medium hosting the scatterers, φt is the angle of the transmitted plane wave (according to Snell's law), chℓ are unknown coefficients, CWℓ(ξ,ζ)=Hℓ(1)(ρ)eiℓθ is the cylindrical wave, and Hℓ(1) is the first-kind Hankel function of integer order ℓ;  moreover, ξhq and ζhq are the distances, along ξ and ζ directions, of the axes of the hth and qth cylinders; finally, δ is the Kronecker symbol.To express the scattered-reflected and scattered-transmitted fields, the presence of the planar surface is taken into account by considering the reflection and transmission of each elementary plane wave constituting the Fourier spectrum of a cylindrical function. The scattered-reflected field is given by the sum of the fields scattered by each buried cylinder and reflected by the interface; by introducing the reflected cylindrical functionsRWℓ [11], the final expression is(2)Vsr⁡=V0∑m=-∞+∞Jm(n1ρq)eimθq∑h=1N∑l=-∞+∞Chlile-ilφt×RWm+l[-n1(χh+χq),n1(ηq-ηh)].
The scattered-transmitted field is given by the sum of the fields scattered by each buried cylinder and transmitted to medium 0. Introducing the transmitted cylindrical functions TWℓ [11], it can be found that(3)Vst=V0∑h=1Nchlile-ilφtTWl(ξh,ζh,χh),
Given the expressions of all the field terms, the boundary conditions on the cylinder surfaces are imposed; due to the hypothesis of perfectly conducting cylinders, the tangential component of the total electric field must vanish. A linear system for the unknown coefficients of the cylindrical wave expansions is obtained. Once such system is solved, the V(ξ,ζ) field is determined in any point of the space and for both polarizations. It is possible to derive all the other components of the electromagnetic field by using Maxwell's equations.


## 2.2. Dielectric Cylinders in a Dielectric Half-Space

When the buried scatterers are dielectric cylinders, the field is present also inside them. This problem is treated in detail in [12]. The expression of the field internal to the qth cylinder is given by an expansion in terms of first-kind Bessel functions with unknown coefficients dqm, as follows:(4)Vcq=V0∑m=-∞+∞dqmime-imφtJm(ncqρq)eimθq,
where ncq is the refractive index of the cylinder.On the cylinder surfaces, the boundary conditions have to be imposed. A linear system is thus obtained, for the unknown coefficients of the cylindrical wave expansions and for thedqm coefficients. Once the system has been solved, the total electromagnetic field is completely determined in any point of the space.


## 2.3. Cylinders in a Finite-Thickness Slab

When the scatterers are inside a finite-thickness slab, among two different half-spaces, additional reflected-transmitted, multiple-reflected and multiple-reflected-transmitted cylindrical functions are introduced, to consider the complicated interaction between the cylinders and the planar interfaces delimiting the slab. The case of perfectly conducting cylinders is solved in [22]; the solution to the scattering problem by dielectric cylinders in a slab is described in [23].It is important to observe that this extension of the method allows to simulate several interesting scenarios, as sketched in Figure3: objects in a soil layer, above a ground of different permittivity (a), objects inside a slab among two identical half-spaces, for example, cylinders in a wall (b), or else scatterers in a slab terminated on a perfectly conducting surface (c), in a layer of a stratified medium (d), or, finally, cylinders above (e) or in (f) a dielectric half-space.Cylinders in a finite-thickness slab, various scenarios.
(a)(b)(c)(d)(e)(f)


## 2.4. Pulsed-Wave Analysis

By using the CWA technique, it is possible to solve the scattering problem by buried perfectly conducting or dielectric cylinders of an incident pulsed plane wave, with a rather general time-domain shape. We have to perform a sampling of the incident-field spectrum and of the spectra of the various field terms involved (reflected, transmitted, scattered, scattered-reflected, scattered-transmitted, etc.). The scattering problem can be solved in the frequency domain for any sample through the CWA. Finally, by means of an inverse transform, the solution in the time domain can be calculated. This procedure is described in detail in [13, 28].


## 3. Numerical Results

Let us consider an underground metallic utility at a standard burial depth of 50 cm: this suggests to operate in a frequency range between 1 and 1.5 GHz. A circular section with a 6 cm diameter is assumed for a target, embedded in soil layer of relative permittivityɛ1=4. In the CWA analysis with normalized quantities, these geometrical values correspond to a depth χ=5π and a radius α=0.6π, at a frequency of 1.5 GHz. The layout is sketched in Figure 4. The total scattered field is evaluated in the upper medium, at the near-field distance of 5 cm, for a normally incident plane wave in TM polarization (electric field parallel to the cylinder axis).Figure 4
Geometrical layout of a buried utility.In Figure5, the scattered field is plotted for a layer of thickness L=120 cm and for different permittivity values of the lower half-space. The relative permittivity of dry soil, such as dry sand, clay, and rock, is assumed as being comprised between 4 and 7. The hypothesis of dry materials is made to better meet the approximation of lossless materials of our analysis;   moreover, in practical surveys the application of GPR to utility detection is limited by attenuation due to wet soil. When ɛ1=ɛ2=4, the scatterer is buried in a semi-infinite medium. The scattered-field shows a behaviour, along a line parallel to the interface, that strongly depends on the presence of the second medium: when the object is in a finite-thickness slab, three peaks can be appreciated in the main lobe of the scattered field, while there are not oscillations in the main lobe when the cylinder is buried in a semi-infinite medium.Figure 5
Total scattered field in air, along a line parallel to the interface, for the layout of Figure4, with L = 120 cm.In Figures6 and 7, the metallic utility is buried in a dielectric half-space with ɛ1=4. In particular, in Figure 6 different values of the cylinder radius are considered, and of course for a bigger utility, a stronger scattered field is observed. In Figure 7, an electric utility with radius a=6 cm is buried at different depths; it can be observed that when it is deeper, the scattered field is less directive; since the soil is considered lossless, we do not obtain a significant reduction of the scattered field amplitude.Total scattered field in air, along a line parallel to the interface, for different values of the radius of the metallic utility.
(a)(b)Total scattered field in air, along a line parallel to the interface, for different values of the burial depth.
(a)(b)Let us now consider an underground air cavity in a dielectric half-space of permittivityɛ1=4. As in previous simulations, the frequency is 1.5 GHz, and the total scattered field is evaluated in the air, at the near-field distance of 5 cm, for a normally incident plane wave in TM polarization. In Figure 8, different values of the cavity radius are considered, the distance between air-soil interface and the target is fixed and it is h-a=10 cm; the electromagnetic effect of a bigger cavity is stronger and the main beam of the scattered field has a more jagged behaviour. In Figure 9 the radius is fixed, a=10 cm, and different values of the burial depth are simulated; considerations similar to those of Figure 7 apply.Total scattered field in air, along a line parallel to the interface, for air cavities of different size.
(a)(b)Total scattered field in air, along a line parallel to the interface, for air cavities located at various depths.
(a)(b)We finally present, in Figure10, time-domain results for three buried metallic utilities. The geometry is sketched in the figure. The cylinder radii are a1=a3=2 cm, a2=4 cm, the burial depths are h1=30 cm, h2=40 cm, h3=35 cm, and the distances along z axis are z1=0,z2=21 cm and z3=42 cm. The polarization of the incident field is TM, the frequency is 1.433 GHz, and the time-domain shape is the one proposed in [13]. The observation point is 15 cm above the ground, aligned with the bigger cylinder axis. We plotted the early-time response of the structure, that is, the total scattered field as a function of time, over a time interval that is of the order of the incident pulse duration and its propagation from the observation point to the targets. Several wavefronts can be identified due to combined effects of the diffraction from the cylinders, the reverberations between the cylinders and the air-ground interface, the creeping-wave circumnavigation of the cylinders.Time-domain results for three buried metallic utilities.
(a)(b)


## 4. Conclusions

In this paper, our work on the development of the cylindrical-wave approach for the solution of the two-dimensional plane wave scattering problem by a set of perfectly conducting and dielectric cylinders buried in a dielectric half-space or in a finite-thickness slab, is shortly reviewed. The proposed full-wave technique deals with both TM and TE polarization cases and yields results in both the near- and the far-field zones. Moreover, it may be applied for any value of the scatterers size and of the distance between the obstacles and the interface between air and soil.The method can be used for the characterization of suitable scenarios in the context of ground-penetrating radar (GPR) applications, which usually employ purely-numerical finite-difference (FD) techniques. The CWA could be exploited in iterative algorithms for the solution of inverse problems, where fast, efficient, and accurate forward solvers are needed to perform repeated evaluations. The approach can be employed to study the scattering of an incident pulsed plane wave, with a rather general time-domain shape.We have presented new numerical results concerning the simulation of buried utilities. In particular, metallic pipes and air cavities have been considered, and the effectiveness of the method for the sensing of cylindrical inhomogeneities buried in the earth has been shown.



---

*Source: 974518-2011-07-25.xml*