Next Article in Journal
Wideband and Channel Switchable Mode Division Multiplexing (MDM) Optical Power Divider Supporting 7.682 Tbit/s for On-Chip Optical Interconnects
Next Article in Special Issue
Broadband Measurements of Soil Complex Permittivity
Previous Article in Journal
A Low-Profile Antenna for On-Body and Off-Body Applications in the Lower and Upper ISM and WLAN Bands
Previous Article in Special Issue
A Method for Extracting Debye Parameters as a Tool for Monitoring Watered and Contaminated Soils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accurate Detection and Localization of Water Pipe Leaks through Model-Based TDR Inversion

1
Department of Electrical and Information Engineering, Polytechnic University of Bari, Via E. Orabona 4, 70125 Bari, Italy
2
Department of Engineering for Innovation, Complesso Ecotekne-Corpo O, University of Salento, 73100 Lecce, Italy
3
Department of Information Engineering, Electronics and Telecommunications (DIET), Sapienza University of Rome, 00184 Rome, Italy
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(2), 710; https://doi.org/10.3390/s23020710
Submission received: 16 November 2022 / Revised: 15 December 2022 / Accepted: 4 January 2023 / Published: 8 January 2023
(This article belongs to the Collection Dielectric Sensing-Based Systems and Applications)

Abstract

:
The problem of water scarcity affects many areas of the world due to water mismanagement and overconsumption and, more recently, to climate change. Monitoring the integrity of distribution systems is, therefore, increasingly important to avoid the waste of clean water. This paper presents a new signal processing technique for enhancing the performance of the methodology of leak detection in water distribution pipes based on time domain reflectometry (TDR). The new technique is based on a particular kind of TDR inversion (spatial TDR) based on a “gray-box” lumped parameter model of the system. The model does not include, e.g., radiative phenomena, non-TEM (transverse electromagnetic) modes etc. but is capable of reproducing accurately the complicated reflectograms obtained by a TDR leak detection system assuming a proper profile of capacitance per unit length along the sensing element. Even more importantly, the model is identified using only the reflectograms taken by the system with very little prior information about the system components. The developed technique is able to estimate with good accuracy, from reflectograms with unclear or ambiguous interpretation, the position and the extension of a region where water is located. The measurement is obtained without prior electromagnetic characterization of the TDR system components and without the need of modeling or quantifying a number of electromagnetic effects typical of on-site measurements.

1. Introduction

Pipeline leakage is a significant problem in the public water supply industry. According to EurEau (European Federation of National Associations of Water Services), the mean leakage rate of drinking water distribution networks throughout Europe is about 25% of all water provided [1]. This means that approximately 2700 m3 of water are wasted yearly per km of pipe. Financial losses for suppliers are the most direct consequences of water waste due to leaks, but this problem has environmental and social consequences too. Leaks can indeed exacerbate the water scarcity issue, which, at certain times of year, is affecting a growing number of people and bigger areas of Europe as a result of climate change [2]. For these reasons, the problem of water leakage has been addressed by a recent European Union (EU) Directive [3], which states that water leakage levels should be monitored and reduced. Reliable and accurate methods for locating leaks in pipes are, therefore, of great interest nowadays.
Many different approaches have been proposed and developed over the years for leak detection in water pipes [4]. The most direct method is based on visual inspection, which relies on the usage of endoscopes or robots equipped with cameras and guided inside pipes. This is a very simple method conceptually, but it also requires a big specialized human intervention for interpreting images and data. Electromagnetic (EM) inspection methods such as magnetic flux leakage, remote field eddy current and ground penetrating radar are also used for status assessment of metallic pipes. These techniques evaluate the metal conditions based on the interaction between transmitted EM signals and the metallic pipes. Other methods are, instead, based on the use of acoustic and ultrasonic waves [5]. The effectiveness of all these methods depends on the kind of material and on the geometry of the pipes, on the working conditions of the water distribution network (e.g., water pressure) and on the environmental conditions. Even more importantly, expert personnel are required for measurement operations and for result interpretation, which prevents them from being considered for continuous monitoring of pipes.
An alternative approach proposed in 2012 [6,7,8] is based on the use of time domain reflectometry (TDR). The method exploits EM waves propagating in a distributed sensing element (SE) placed near the pipe (typically less than 1 cm from the pipe) to detect and localize water leaks, as shown in Figure 1. The SE is a transmission line (TL), e.g., a bi-wire cable, whose primary parameters are affected by changes in the dielectric constant of the surrounding medium. The presence of water leaking from the pipe increases locally the dielectric constant of the medium and can, therefore, be detected by a reflectogram, as a reflection caused by the locally increased capacitance per unit length. The reflectogram is obtained by transmitting a stimulus signal from one end of the SE and measuring the reflected signal at the same point. In this way, long segments of the pipe can be analyzed with just one measurement. As a matter of fact, SEs up to 100 m long have been successfully used for water leak detection [6]. As of the current date, systems of this kind have been installed in more than 50 km of water pipes by AQP (public water operator of the Apulia region, Italy, and the largest water facility in Europe), proving the effectiveness of the methodology.
With respect to other available techniques for water leak detection in pipes, the TDR-based method has the advantage of not imposing requirements on the construction material of the pipes (because water in the soil around the pipe is detected) or on their operating conditions, leading to greater ease of use. TDR-based techniques also drastically reduce the amount of time needed for inspection compared to conventional leak detection techniques; in fact, hundreds-of-meters-long pipes can be inspected with a single measurement operation. Besides, TDR is a well-established technique already used for many applications such as detection and localization of cable defects [9,10,11,12], monitoring of landslide [13,14], characterization of antennas [15] and measurement of dielectric properties of materials [16,17], of liquid level in tanks [18], of soil moisture [19], of water diffusion in irrigation [20] among many others, even in the biomedical field [21]. A possible disadvantage of TDR-based water leak detection is that the SE must be installed near the water pipe, and, hence, it is difficult to use the technique in existing pipelines. This is not very limiting, however, since the SE can be installed during usual maintenance operations, which provide for the periodic replacement of sections of pipes.
A limitation of the TDR method is that the local increase in the dielectric constant caused by the leak must be clearly visible in the reflectogram in terms of a clear difference [7] between the “reference reflectogram” (taken without leak), and the “measurement reflectogram” (taken in presence of leak). As better shown in the next sections, this is not always the case: quite often, an incipient water leak causes a comparatively small variation of dielectric constant, which cannot be accurately detected visually, since the corresponding variation in the reflectogram is small or masked by other disturbance factors and unwanted effects which typically influence on-field measurement. It is worth noting that, in this case, evaluating a simple “difference reflectogram” between the measurement and the reference reflectogram is not a solution: a local change of dielectric constant causes changes in the whole reflectogram due to reflections, changes in the propagation velocity, etc. This is clearly observable in the examples reported in Section 3.
As a matter of fact, when the measurement reflectogram is not clearly interpretable, a TDR inversion–or spatial TDR–is a key issue for localizing, more accurately than an intuitive visual interpretation, the change in the spatial profile of dielectric constant. This operation is also useful for an accurate localization in addition to a reliable detection. Indeed, it must be noted that having a reasonably accurate measurement of the leak position is very important for the repair team, who is required to carry out excavations.
TDR inversion algorithms have already been presented and used successfully in the literature but in experimental situations quite different from the one here considered. Usually, the inversion problem is reduced to the minimization of an objective function whose gradient has a known analytical expression. In [22] and [23], the TDR inversion works well but for simulated ideal models. In [24], the TDR inversion is performed for an experimental setup, which uses a soil moisture TDR probe with accurately known electrical parameters (a classic flat-ribbon cable with three conductors) and in conditions allowing the use of a very simple model of its primary parameters. The same kind of TDR probe is used in [25], and, in this case too, the electrical parameters of the probe are accurately determined and validated through laboratory experiments and full-wave simulations. A satisfactory TDR inversion is also obtained in [26]: in this case, the probe is carefully characterized with frequency domain measurements, performed using a properly calibrated VNA (Vector Network Analyzer). In general, in all TDR inversions known to the authors, determining the parameters of the employed probe and of the other components of the measurement system is a preliminary step, necessary to performing the actual inversion.
In the set-up depicted in Figure 1, an accurate characterization of the measurement system is unavailable and cannot be obtained a priori with, e.g., VNA measurements. First of all, the probe, i.e., the SE, is a buried bi-wire cable. Therefore, it is tens of meters long, inevitably twisted and stressed, and it is surrounded by a medium with an uneven and unknown profile of dielectric constant. Characterizing a short section of the SE in air is not simple (it is a balanced line, not matched with the ports of a VNA), but, more importantly, the characterization would be of no use to model the cable when it is buried near the pipe.
In the experimental situation of Figure 1, other parameters are also not known with sufficient accuracy: namely, the electromagnetic parameters of the coaxial cable and of its connection with the TDR unit are not necessarily the same that could be determined in laboratory making use of dedicated instrumentation. These parameters should be determined for the actual setup on the field. The same applies, of course, to the SE termination.
Consequently, classic TDR inversion methods are substantially inapplicable to the problem at hand: a method that does not require prior characterization of the measurement system components is needed. It must also be underlined that classic methods often use reflectograms taken from both ends of the probe (in early works like [22,23], and in recent works like [27,28]) since TDR inversion with a single reflectogram, under the hypotheses considered in these papers, is mathematically an ill-posed problem. Instead, a TDR inversion with a single reflectogram, like in [29], is highly desirable for the considered problem of water-leak detection, but it is necessary not to limit too severely the class of identifiable capacitance profiles.
In this paper, a TDR inversion method suitable for water-leak localization in underground pipes, using a setup similar to that described in Figure 1, is proposed. The basic idea is to identify a circuital model of the on-field measurement system, able to reproduce very accurately both the “reference reflectogram” and the “measurement reflectogram” even if:
(1)
the reflectograms are very irregular, a condition inherent to the situation of a long SE buried along a pipe;
(2)
the electromagnetic parameters of the SE, and of other components of the measurement system, are not known in advance, and must be determined directly on the field.
This model is not designed to represent accurately the physical reality of the system (which includes non-TEM (transverse electromagnetic) modes, radiative phenomena, interfering fields, non-ideality of connections, etc.) but simply to reproduce the observed reflectograms including the multiple reflections (voltage recorded after the first reflected pulse).
The circuital model, and the identification procedure, are then used for leak detection and localization, based on measuring the capacitance profile needed to reproduce the reflectogram as it derives from the presence of a leak. As shown in the next sections, the proposed technique has clear advantages with respect to traditional techniques when analyzing reflectograms measured in the presence of water leaks which produce small variations of dielectric constant. In these cases, indeed, results obtained with traditional techniques based on the observation of the reflectogram are very ambiguous while the proposed algorithm is able to effectively detect and localize the water leak.
It is worth noting that the proposed method can be easily extended to a number of other TDR-based applications, such as those for detecting soil moisture, cable faults, defects in transmission lines, etc. The advantage of the method is that almost no prior knowledge of the electrical parameters of the probes and of the measurement system components is required, which makes it more flexible and general, and much easier to apply. Investigating such different application of the method, of course, is the subject for separate works.

2. Materials and Methods

2.1. Reflectograms in TDR Leaks Detection in Water Pipes

The on-field measurement setup for TDR-based leak localization in underground water pipes (Figure 1) has been already illustrated and discussed with some detail in the Introduction. We add here that the TDR unit can be a dedicated commercial device or can be realized with general purpose instruments (e.g., waveform generator + oscilloscope).
A typical signal obtained using the described measurement setup, with the SE terminated in open circuit, is shown in Figure 2. An impulsive voltage waveform is transmitted to the SE and reflected at different points. The first significant reflection occurs at the interface between the coaxial cable and the SE and is due to the quite different characteristic impedance of the connected TLs. The reflection at the SE end is also very strong, due to the open termination. The reflection at the water leak is also visible in the example, but, in this case, its amplitude is comparable to other small reflections due to non-idealities of the SE (e.g., twists and stress on the bi-wire) and to the inhomogeneity of the medium around it. In this TDR signal, the water leak is hardly detectable and localizable. Detection and localization in this situation is exactly the goal of the method illustrated in this work.

2.2. Laboratory Experimental Setup

For the purpose of the present investigation, an “ad-hoc” laboratory set-up has been arranged, described in Figure 3. The TDR unit is implemented using an Agilent 33250A arbitrary waveform generator (AWG) for the stimulus signal and a LeCroy Waverunner-2 LT262 oscilloscope for signal acquisition. The AWG and the SE are connected to the oscilloscope through a T-junction and a short RG58-CU coaxial cable. The SE is a bi-wire cable having the cross-section shown in Figure 3. The two conductors of the SE are composed of 24 copper-clad aluminum wires coated with polyvinyl chloride (PVC).
The whole setup is designed to perform several laboratory experiments with many different accurately known distributions of water near the SE. Nonetheless, it is not so easy to obtain an accurate electromagnetic model of the system, and this situation mimics an on-field measurement setup.
First of all, the SE is 15 m long and lies freely on the floor: it is, therefore, unavoidably twisted at random points. The parameters of the TL constituted by the SE are influenced by the floor and are not the same as would be observed in air. The SE is not shielded from electromagnetic fields in the laboratory. The connection between the (unbalanced) coaxial cable and the (balanced) SE is realized with simple crocodile clips, and, therefore, its effect on the reflectogram, although certainly not very large, is not known.
It is useful to note that, in practical on-field leak measurements, the SE has more distance between the wires and, therefore, higher sensitivity. In laboratory experiments it is preferable to use an SE with low sensitivity, since it makes it simpler to obtain realistic reflections associated with leaks using only water, instead of a mix of soil and water.
More specifically, in our experiments water leaks are simulated by (method 1) submerging the SE in water or (method 2) by placing a water container on the SE. Method 1 produces a visible reflection, like those observed in the field with a more sensitive SE, when a big leak soaks the soil surrounding it. Method 2 produces a much weaker reflection, like those observed in the field when the leak is smaller and the surrounding soil is moist but not soaked. It is worth noting that the placement of the water container above the SE is a practical device to obtain a weaker reflection in laboratory tests and does not replicate an actual on-field condition. In installed systems, indeed, all the soil around the SE is moist due to water leaking from the pipe.
Typical measured signals that are analyzed in the paper are depicted in Figure 4. Signal (a) is acquired in normal (i.e., dry) conditions. It can be a “reference reflectogram”. Signal (b) is acquired with method 1, with a water container length of 33.5 cm (water height is, for this and the other experiments, about 2 cm). The pulse reflected by water, highlighted by a red box, is clearly observed. Signal (c) is acquired with method 2, with the same 33.5 cm long water container. Due to the distance between the SE and the water (separated by the bottom of the 2 mm thick container), the parameters of the TL are less influenced, and, hence, a smaller reflection is produced. The last signal, in Figure 4d, is also acquired with method 2 but, in this case, the water container is 2 m long. Therefore, this case is similar to (c), but the water is distributed in a wider region.

2.3. Measurement System Modeling

The first step for locating the water leaks via TDR inversion is modeling the measurement system. The adopted model is represented in Figure 5.
For the coaxial cable and the SE, a classic lumped-element model with pi cells is used, terminated on the load Z L . The two TLs are composed of N C and N S E cells, respectively, of length d z . The total lengths of the cables are therefore, respectively,
l c = N c   d z
l S E = N S E   d z
The following frequency-dependent models are used for the primary parameters R , L , G , C , for both the coaxial cable and the SE:
R ( ω ) = R 0 + R 1 ω
L ( ω ) = L 0 + L 1 / ω
G ( ω ) = G 0   ω
C = C 0
where ω is the angular frequency and R 0 ,   R 1 ,   L 0 ,     L 1 ,   G 0 ,     C 0 are six adjustable parameters. The employed functions of the frequency are justified by the dielectric dispersion and the skin effect phenomena. It can be easily verified that, neglecting these effects, the model would be able to reproduce only very coarsely the actual observed reflectograms.
The capacitance of the SE, which is the parameter of interest in the analysis (influenced by the permittivity of the medium around the SE), is independent of the frequency, but dependent on the position z = n d z . Thus, the model of the capacitance per unit length of the n -th cell is
C 0 ( z ) = C 0 , n = C S E p n
where p n is a multiplicative parameter applied the n -th cell. The circuital model is completed by that of the waveform generator and of the oscilloscope, also depicted in Figure 5, and by the impedance Z L of the SE termination.
The described circuital model is clearly only an approximation of the physical reality and is not intended to provide a faithful representation of the electromagnetic behavior of the measurement system, as a full-wave simulation can do. Non-TEM modes and radiative phenomena, for example, are not modeled. However, the model can produce reflectograms that fit on-field measurements very well, as shown in Section 3. Moreover, the model is much simpler and computationally cheaper than full-wave models and can, therefore, be effectively integrated in the simulation-based optimization routine needed for its identification.

2.4. Measurement System Simulation

The model is used to actually simulate the measurement system and generate the reflectograms corresponding to given values of all the parameters. To this purpose, it is implemented in the simulator described in [30], which works in the frequency domain and is able to compute efficiently the output, handling frequency-dependent primary parameters. In order to obtain good simulation results, the cell length must be chosen conveniently small with respect to the minimum wavelength λ m i n propagating in the cables. The cell length is given by:
d z = λ m i n / α λ λ m i n = min { min f [ 0 , f max ]   1 f L C ( f ) C C , min f [ 0 , f max ] , z 1 f L S E ( f ) C S E ( z ) }
where f m a x is the maximum frequency considered in the simulation, and α λ is a tunable parameter. The frequency f m a x is determined by the impulsive signal of the generator; a rule of thumb for the choice of α λ is α λ 10 [30]. The value of this parameter determines the density of the TL model and, hence, the accuracy in simulating a continuous system. In this work, the value α λ = 20 has been used; greater values lead to negligible improvements and higher computation time. Since during the optimization the values of the parameters and, therefore, of λ m i n , change at each iteration, the number of cells is recomputed each time, according to Equations (8), (1) and (2).

2.5. Measurement System Identification

The system model is identified by minimizing the root mean square error (RMSE) between the measured reflectogram and the simulated one. Minimizing this objective function is conceptually very simple but requires some attention for its actual implementation.
It is important to underline that, in principle, all the parameters of the model described in Section 2.3 can be involved in the optimization process, including the impedance of the generator and of the oscilloscope. We kept these impedances fixed to their nominal values reported in Figure 5, since optimizing them proved, in our experiments, not useful to improve the fitting. We also kept fixed the length l C of the coaxial cable used for the connection and the length l S E of the SE, which are easy to be known in advance with sufficient accuracy.
It is convenient for the optimization process to introduce a functional representation with a fixed number of parameters of the position-dependent capacitance profile p n , since the number of cells depends on λ m i n and is not fixed (besides being quite high, N S E 1000 in our tests). A possibility is the multi-sine signal, which represents well a general band-limited function. In our case, however, we opted for representing p n as a sum of N g equispaced Gaussian pulses, with standard deviation equal to the distance between the peaks of two consecutive pulses, and different amplitudes, which are, therefore, the free parameters. This model is more effective for the optimization process, because each parameter influences a limited region of the profile, and not the whole profile as in the case of the multi-sine signal. It was determined that N g = 100 parameters are able to describe the capacitance profile with sufficient generality for typical measured reflectograms. This parameter can be increased, however, when working with more irregular reflectograms.
As a final note, the termination impedance Z L is modeled as a capacitance to consider the fringing at the end of the SE. This model proved to lead to better simulation results with respect to an ideal infinite impedance, which causes the presence of a peak in the final part of p n (an effect that clearly compensates the unmodeled fringing capacitance).
The complete list of the free parameters of the model, subject to the optimization process, is listed here.
  • Primary parameters of the coaxial cable: R 0 , c ,   R 1 , c ,   L 0 , c ,   L 1 , c ,   G 0 , c ,   C 0 , c .
  • Primary parameters of the SE: R 0 , S E ,   R 1 , S E ,   L 0 , S E ,   L 1 , S E ,   G 0 , S E ,   C 0 , S E .
  • Vector of 100 amplitudes of the Gaussian pulses composing the capacitance profile p n .
  • Termination capacitance: C L .
Determining and handling an analytical expression of the gradient of the objective function is an unrealistic task due to the complexity of the frequency-dependent model. Given also the large dimensionality of the parameters vector, a random search method with adaptive step size is effectively used for the optimization [31]. Discussing in detail the numerical optimization algorithm is beyond the scope of this work. The results, instead, deserve a careful analysis, presented in the next section.

2.6. Water-Leak Detection and Localization

There are two ways to detect and localize water leaks from the estimated capacitance profile p n . The first does not use any reference reflectogram but is more prone to errors and ambiguities in interpreting the profile. The second must use the reference reflectogram but gives automatic and univocal results.
When reflections due to the presence of water are strong, like in Figure 4b, water leaks can be directly localized as prominent peaks in the estimated capacitance profile p n . We refer conventionally to this procedure as “direct water-leak detection”.
When reflections due to the presence of water are weak, like in Figure 4c,d, the estimated profile does not provide a solid identification of leaks. In these cases, the reference signal measured with the same setup but without water (Figure 4a) is exploited in the following way. We refer conventionally to this procedure as “water-leak detection after calibration”.
  • A first model identification is performed on the reference reflectogram. All the model parameters are determined. This is the calibration phase, which gives the profile p n = p c a l ( z ) (with z = n d z ).
  • A second model identification is performed on measurement reflectogram. In this case, the formerly determined parameters are fixed. The fitting is obtained by acting only on the parameters of an additional capacitance profile p l e a k ( z ) , i.e., assuming the capacitance profile p ( z ) = p c a l ( z ) + p l e a k ( z ) . This is the measurement phase.
The additional capacitance profile associated with a water leak is modeled as a function of the form:
p l e a k ( z ) = k β ( z z 0 w , α , α )
where β ( x , a , b ) is the beta probability density function computed in x , with form parameters a and b (see Figure 6). In this way, p l e a k ( z ) is non-zero in the interval [ z 0 , z 0 + w ] (it begins at z 0 and has a width w ); k is a magnitude factor; α is a form factor. Depending on the value of α , model (9) assumes different shapes, from a rectangular pulse to a bell-shaped pulse, as depicted in Figure 6.

3. Results

The water-leaks estimation technique has been applied to signals measured with the experimental setup described in Section 2.2.

3.1. Direct Water-Leak Detection

When the signal reflected at the water leak is strong (case (b) of Figure 4), the region of the SE submerged in water is accurately identified, as shown in Figure 7.
When the signal reflected at the water leak is weak (case (c) of Figure 4), results are of the kind illustrated in Figure 8. The region where water is present is still recognizable from the fitted capacitance profile, since it corresponds to its highest peak. However, the magnitude of the peak is not so high with respect to the fluctuations in the profile due to the non-idealities of the SE; the detection and the localization are clearly prone to subjective interpretations. In this case, the procedure “water-leak detection after calibration” is justified.

3.2. Water Leak Detection after Calibration

3.2.1. Calibration Phase

The result of the calibration phase of the procedure is shown in Figure 9. The figure shows that the obtained fitting follows the uneven oscillations of the measured reflectogram very well, even after the principal reflected pulse (i.e., even multiple reflections are reproduced quite well). It is important, however, to point out that the estimated profile must not be interpreted as the actual physical capacitance profile along the SE. The model greatly simplifies the electromagnetic behavior of the system, and, even more important, the transmitted pulse is quite large: to accurately estimate the capacitance profile, a greater spatial resolution is needed, and, hence, a shorter transmitted pulse must be used (this kind of measurement needs costly instrumentation not suitable for on-field use; it has been obtained in some preliminary tests, not of interest here).
The obtained profile, therefore, must be interpreted simply as one of the profiles that reproduces well the input/output behavior of the system. This is sufficient, however, to detect and locate quite accurately the water leak.

3.2.2. Measurement Phase—Case (b)

Using the calibration described above, the measurement phase of the procedure “water-leak detection after calibration”, as described in 2.6, is applied to the reflectogram of Figure 7, obtaining the results reported in Figure 10. The estimated center and length of the water region are 4.50 m and 51 cm, respectively, while the actual ones are 4.52 m and 46 cm, respectively.

3.2.3. Measurement Phase—Case (c)

Using the same calibration, the measurement phase is applied to the reflectogram of Figure 8. The results are depicted in Figure 11. Despite the very small reflection due to the water, the water region is correctly and accurately identified.
The same test is repeated by placing the water container in other positions, and results are reported in Table 1. The estimated water region always intersects the actual water region, which can be classified as a successful localization. The standard deviation of the error in the leak-center estimation and that in the leak-length estimation are both approximatively equal to 13 cm. In Table 1 is also reported the total estimated capacitance of the leak, i.e., the area below the estimated additional capacitance profile. This figure can be useful as an index of the size of the leak, in terms of quantity of water sensed by the SE.
The signal in Figure 11 has also been analyzed through traditional water leak detection techniques, for comparison. Traditional techniques are mainly based on direct observation of the measured reflectogram which leads to the identification of reflected signals. In this case, in which the reflected pulse is very small, the “difference reflectogram” between the measurement reflectogram and the reference reflectogram (depicted in Figure 9) is more informative. The difference reflectogram has been, therefore, computed and is depicted in Figure 12. The Figure also contains an indicative distance scale that correlates time instants to distances from the beginning of the SE. The origin of the scale corresponds to the peak of the pulse reflected at the coaxial/SE interface, while the length of the SE (15 m) corresponds to the peak of the pulse reflected at SE termination. In the difference reflectogram, the reflection caused by the presence of water is clearly visible. Oscillations in the signal are, however, present in the whole region from ~ 1 m to ~ 5 m, making it difficult to accurately localize the leak and reducing, therefore, the practical applicability of the method.

3.2.4. Measurement Phase—Case (d)

Water leak detection is then carried out for 2 m long simulated leaks, i.e., case (d) described in Section 2.2 and in Figure 4; the calibration is that of Section 3.2.1. The results are reported in Table 2. In this case, water leaks are correctly localized in each measurement (estimated and actual water regions always intersect). The standard deviations of the error in the leak center estimation and in the leak length estimation are approximately equal to 33 cm and 43 cm, respectively.

4. Conclusions

In this paper, a TDR inversion method for enhanced TDR-based water-leak detection and localization is proposed. The detection and localization of water leaks is performed by processing reflectograms taken on a bi-wire sensing element placed alongside the buried pipes. TDR inversion uses a model designed to reproduce accurately the uneven and oscillating reflectograms obtained in a leak detection campaign and to be identified even if very little prior information on the model parameters is available, as is unavoidable in this particular kind of measurements. The TDR inversion can localize water leaks working on a single signal, when the reflection caused by the water is strong enough, and exploits an initial calibration carried out in the absence of water, in case of weaker reflections. Experimental results carried out in different conditions show that the method is able to produce accurate estimates of the position and the extension of water leaks. The proposed method has, therefore, the potential to be used successfully also for continuous monitoring of water pipes and has the potential of being employed for other TDR-based applications, since it basically is able to perform a TDR inversion with no prior characterization of the TDR probes and of other components of the measurement system.

Author Contributions

Conceptualization, A.C., N.G., E.P., M.S. (Marco Scarpetta) and M.S. (Maurizio Spadavecchia); methodology, N.G., E.P. and M.S. (Marco Scarpetta); software, M.S. (Marco Scarpetta); validation, A.C., A.M., E.P., M.S. (Marco Scarpetta) and M.S. (Maurizio Spadavecchia); formal analysis, N.G., E.P. and M.S. (Marco Scarpetta); investigation, M.S. (Marco Scarpetta); resources, N.G. and M.S. (Maurizio Spadavecchia); data curation, A.C., N.G., A.M., E.P., M.S. (Marco Scarpetta) and M.S. (Maurizio Spadavecchia); writing—original draft preparation, M.S. (Marco Scarpetta); writing—review and editing, A.C., N.G., A.M., E.P. and M.S. (Maurizio Spadavecchia); visualization, A.C., A.M., E.P., M.S. (Marco Scarpetta) and M.S. (Maurizio Spadavecchia); supervision, A.C., N.G. and M.S. (Maurizio Spadavecchia); project administration, N.G.; funding acquisition, N.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Academic Research Fund (ARF 2021)—Polytechnic University of Bari: “Processing of reflectometric signals for the characterization of cables, materials and sound absorbing fabrics”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. EurEau. Europe’s Water in Figures—An Overview of the European Drinking Water and Waste Water Sectors; EurEau: Brussels, Belgium, 2021. [Google Scholar]
  2. European Environment Agency. Water Resources across Europe—Confronting Water Stress: An Updated Assessment; EEA: Brussels, Belgium, 2021. [Google Scholar]
  3. EU Directive (EU) 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the Quality of Water Intended for Human Consumption 2020. Available online: https://eur-lex.europa.eu/eli/dir/2020/2184 (accessed on 3 January 2023).
  4. Liu, Z.; Kleiner, Y. State of the Art Review of Inspection Technologies for Condition Assessment of Water Pipes. Measurement 2013, 46, 1–15. [Google Scholar] [CrossRef] [Green Version]
  5. Yu, Y.; Safari, A.; Niu, X.; Drinkwater, B.; Horoshenkov, K.V. Acoustic and Ultrasonic Techniques for Defect Detection and Condition Monitoring in Water and Sewerage Pipes: A Review. Appl. Acoust. 2021, 183, 108282. [Google Scholar] [CrossRef]
  6. Cataldo, A.; De Benedetto, E.; Cannazza, G.; Masciullo, A.; Giaquinto, N.; D’Aucelli, G.M.; Costantino, N.; De Leo, A.; Miraglia, M. Recent Advances in the TDR-Based Leak Detection System for Pipeline Inspection. Measurement 2017, 98, 347–354. [Google Scholar] [CrossRef]
  7. Cataldo, A.; Cannazza, G.; Benedetto, E.D.; Giaquinto, N. A TDR-Based System for the Localization of Leaks in Newly Installed, Underground Pipes Made of Any Material. Meas. Sci. Technol. 2012, 23, 105010. [Google Scholar] [CrossRef]
  8. Cataldo, A.; Cannazza, G.; De Benedetto, E.; Giaquinto, N. Experimental Validation of a TDR-Based System for Measuring Leak Distances in Buried Metal Pipes. Prog. Electromagn. Res. 2012, 132, 71–90. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, J.; Zhang, Y.; Guan, Y. Analysis of Time-Domain Reflectometry Combined with Wavelet Transform for Fault Detection in Aircraft Shielded Cables. IEEE Sens. J. 2016, 16, 4579–4586. [Google Scholar] [CrossRef]
  10. Giaquinto, N.; Scarpetta, M.; Spadavecchia, M. Algorithms for Locating and Characterizing Cable Faults via Stepped-Frequency Waveform Reflectometry. IEEE Trans. Instrum. Meas. 2020, 69, 7271–7280. [Google Scholar] [CrossRef] [Green Version]
  11. Scarpetta, M.; Spadavecchia, M.; Adamo, F.; Ragolia, M.A.; Giaquinto, N. Detection and Characterization of Multiple Discontinuities in Cables with Time-Domain Reflectometry and Convolutional Neural Networks. Sensors 2021, 21, 8032. [Google Scholar] [CrossRef]
  12. Scarpetta, M.; Spadavecchia, M.; Andria, G.; Ragolia, M.A.; Giaquinto, N. Analysis of TDR Signals with Convolutional Neural Networks. In Proceedings of the 2021 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Glasgow, UK, 17–20 May 2021; pp. 1–6. [Google Scholar]
  13. Ho, S.-C.; Chen, I.-H.; Lin, Y.-S.; Chen, J.-Y.; Su, M.-B. Slope Deformation Monitoring in the Jiufenershan Landslide Using Time Domain Reflectometry Technology. Landslides 2019, 16, 1141–1151. [Google Scholar] [CrossRef]
  14. Lin, Y.-S.; Chen, I.-H.; Ho, S.-C.; Chen, J.-Y.; Su, M.-B. Applying Time Domain Reflectometry to Quantification of Slope Deformation by Shear Failure in a Landslide. Environ. Earth. Sci. 2019, 78, 123. [Google Scholar] [CrossRef]
  15. Cataldo, A.; Monti, G.; De Benedetto, E.; Cannazza, G.; Tarricone, L.; Catarinucci, L. Assessment of a TD-Based Method for Characterization of Antennas. IEEE Trans. Instrum. Meas. 2009, 58, 1412–1419. [Google Scholar] [CrossRef]
  16. Wilczek, A.; Szypłowska, A.; Kafarski, M.; Skierucha, W. A Time-Domain Reflectometry Method with Variable Needle Pulse Width for Measuring the Dielectric Properties of Materials. Sensors 2016, 16, 191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Ingole, S.A.; Deshmukh, A.R.; Shinde, R.V.; Kumbharkhane, A.C. Dielectric Relaxation Study of Aqueous Tetraethylene Glycol Using Time Domain Reflectometry Technique in the Frequency Range 10 MHz to 50 GHz. J. Mol. Liq. 2018, 272, 450–455. [Google Scholar] [CrossRef]
  18. Giaquinto, N.; D’Aucelli, G.M.; De Benedetto, E.; Cannazza, G.; Cataldo, A.; Piuzzi, E.; Masciullo, A. Criteria for Automated Estimation of Time of Flight in TDR Analysis. IEEE Trans. Instrum. Meas. 2016, 65, 1215–1224. [Google Scholar] [CrossRef]
  19. Piuzzi, E.; Cataldo, A.; Cannazza, G.; De Benedetto, E. An Improved Reflectometric Method for Soil Moisture Measurement Exploiting an Innovative Triple-Short Calibration. IEEE Trans. Instrum. Meas. 2010, 59, 2747–2754. [Google Scholar] [CrossRef]
  20. Cataldo, A.; De Benedetto, E.; Cannazza, G.; Piuzzi, E.; Giaquinto, N. Embedded TDR Wire-like Sensing Elements for Monitoring Applications. Measurement 2015, 68, 236–245. [Google Scholar] [CrossRef]
  21. Schiavoni, R.; Monti, G.; Piuzzi, E.; Tarricone, L.; Tedesco, A.; De Benedetto, E.; Cataldo, A. Feasibility of a Wearable Reflectometric System for Sensing Skin Hydration. Sensors 2020, 20, 2833. [Google Scholar] [CrossRef]
  22. Lundstedt, J.; He, S. A Time-Domain Optimization Technique for the Simultaneous Reconstruction of the Characteristic Impedance, Resistance and Conductance of a Transmission Line. J. Electromagn. Waves Appl. 1996, 10, 581–601. [Google Scholar] [CrossRef]
  23. Norgren, M.; He, S. An Optimization Approach to the Frequency-Domain Inverse Problem for a Nonuniform LCRG Transmission Line. IEEE Trans. Microw. Theory Tech. 1996, 44, 1503–1507. [Google Scholar] [CrossRef]
  24. Schlaeger, S. A Fast TDR-Inversion Technique for the Reconstruction of Spatial Soil Moisture Content. Hydrol. Earth Syst. Sci. 2005, 9, 481–492. [Google Scholar] [CrossRef]
  25. Scheuermann, A.; Huebner, C.; Schlaeger, S.; Wagner, N.; Becker, R.; Bieberstein, A. Spatial Time Domain Reflectometry and Its Application for the Measurement of Water Content Distributions along Flat Ribbon Cables in a Full-Scale Levee Model: Spatial Time Domain Reflectometry. Water Resour. Res. 2009, 45, W00D24. [Google Scholar] [CrossRef] [Green Version]
  26. Minet, J.; Lambot, S.; Delaide, G.; Huisman, J.A.; Vereecken, H.; Vanclooster, M. A Generalized Frequency Domain Reflectometry Modeling Technique for Soil Electrical Properties Determination. Vadose Zone J. 2010, 9, 1063–1072. [Google Scholar] [CrossRef]
  27. Scheuermann, A.; Yermán, L.; Bittner, T.; Serna, M.L.; Zárate, S.; Bajodek, M.; Bore, T. Design of Two Large Coaxial Cells for Studying Physical Process Using Spatial TDR. In Proceedings of the 2018 12th International Conference on Electromagnetic Wave Interaction with Water and Moist Substances (ISEMA), Lublin, Poland, 4–7 June 2018; pp. 1–3. [Google Scholar]
  28. Yan, G.; Bore, T.; Li, Z.; Schlaeger, S.; Scheuermann, A.; Li, L. Application of Spatial Time Domain Reflectometry for Investigating Moisture Content Dynamics in Unsaturated Loamy Sand for Gravitational Drainage. Appl. Sci. 2021, 11, 2994. [Google Scholar] [CrossRef]
  29. Greco, R. Soil Water Content Inverse Profiling from Single TDR Waveforms. J. Hydrol. 2006, 317, 325–339. [Google Scholar] [CrossRef]
  30. D’Aucelli, G.M.; Giaquinto, N.; Andria, G. LineLab-A Transmission Line Simulator for Distributed Sensing Systems: Open-Source MATLAB Code for Simulating Real-World Transmission Lines. IEEE Antennas Propag. Mag. 2018, 60, 22–30. [Google Scholar] [CrossRef]
  31. Schumer, M.; Steiglitz, K. Adaptive Step Size Random Search. IEEE Trans. Autom. Control. 1968, 13, 270–276. [Google Scholar] [CrossRef]
Figure 1. Representation of the on-field measurement setup for the leak localization in underground pipes based on time domain reflectometry (TDR). The sensing element (SE) is placed in the ground near the water pipe and is connected to a TDR unit through a short coaxial cable.
Figure 1. Representation of the on-field measurement setup for the leak localization in underground pipes based on time domain reflectometry (TDR). The sensing element (SE) is placed in the ground near the water pipe and is connected to a TDR unit through a short coaxial cable.
Sensors 23 00710 g001
Figure 2. Typical reflectogram associated to the on-field measurement setup schematized in Figure 1, when an impulsive transmitted signal is used by the TDR unit.
Figure 2. Typical reflectogram associated to the on-field measurement setup schematized in Figure 1, when an impulsive transmitted signal is used by the TDR unit.
Sensors 23 00710 g002
Figure 3. Representation of the laboratory measurement setup. The SE is placed on the floor and is connected to the TDR unit through a short coaxial cable.
Figure 3. Representation of the laboratory measurement setup. The SE is placed on the floor and is connected to the TDR unit through a short coaxial cable.
Sensors 23 00710 g003
Figure 4. Examples of measured signals, SE length: 15 m. The depicted signals have been obtained in the following conditions: (a) SE in absence of water; (b) method 1, 33.5 cm long water container; (c) method 2, 33.5 cm long water container; (d) method 2, 2 m long water container. In each case, the position of the water container is 4 m from the first end of the SE, and the water height is 2 cm in all cases. Reflected pulses corresponding to water-leak conditions are marked with red boxes.
Figure 4. Examples of measured signals, SE length: 15 m. The depicted signals have been obtained in the following conditions: (a) SE in absence of water; (b) method 1, 33.5 cm long water container; (c) method 2, 33.5 cm long water container; (d) method 2, 2 m long water container. In each case, the position of the water container is 4 m from the first end of the SE, and the water height is 2 cm in all cases. Reflected pulses corresponding to water-leak conditions are marked with red boxes.
Sensors 23 00710 g004
Figure 5. Model of the measurement system.
Figure 5. Model of the measurement system.
Sensors 23 00710 g005
Figure 6. Examples of profiles p l e a k ( z ) obtained using model (9).
Figure 6. Examples of profiles p l e a k ( z ) obtained using model (9).
Sensors 23 00710 g006
Figure 7. Results of the “direct water-leak estimation” in the case of a strong reflected pulse: (a) measured and fitting signal; (b) estimated capacitance profile.
Figure 7. Results of the “direct water-leak estimation” in the case of a strong reflected pulse: (a) measured and fitting signal; (b) estimated capacitance profile.
Sensors 23 00710 g007
Figure 8. Results of the “direct water leak estimation” in the case of a weak reflected pulse: (a) measured and fitting signal; (b) estimated capacitance profile.
Figure 8. Results of the “direct water leak estimation” in the case of a weak reflected pulse: (a) measured and fitting signal; (b) estimated capacitance profile.
Sensors 23 00710 g008
Figure 9. Fitting of the measured signal in the absence of water leaks: (a) measured and fitting signal; (b) estimated capacitance profile.
Figure 9. Fitting of the measured signal in the absence of water leaks: (a) measured and fitting signal; (b) estimated capacitance profile.
Sensors 23 00710 g009
Figure 10. Results of the measurement phase for the reflectogram of Figure 7: (a) measured and fitting signal; (b) estimated additional capacitance profile.
Figure 10. Results of the measurement phase for the reflectogram of Figure 7: (a) measured and fitting signal; (b) estimated additional capacitance profile.
Sensors 23 00710 g010
Figure 11. Results of the measurement phase for the reflectogram of Figure 8: (a) measured and fitting signal; (b) estimated additional capacitance profile.
Figure 11. Results of the measurement phase for the reflectogram of Figure 8: (a) measured and fitting signal; (b) estimated additional capacitance profile.
Sensors 23 00710 g011
Figure 12. Difference between the reflectogram measured in the presence of water (Figure 11) and the reference reflectogram (Figure 9). The indicative distance scale depicted in the Figure starts from the position of the pulse reflected at the coaxial/SE interface and ends at the position of the pulse reflected at the SE termination.
Figure 12. Difference between the reflectogram measured in the presence of water (Figure 11) and the reference reflectogram (Figure 9). The indicative distance scale depicted in the Figure starts from the position of the pulse reflected at the coaxial/SE interface and ends at the position of the pulse reflected at the SE termination.
Sensors 23 00710 g012
Table 1. Estimation results for the water leak position and entity estimation in the case of leak length 33.5 cm.
Table 1. Estimation results for the water leak position and entity estimation in the case of leak length 33.5 cm.
Nominal Center
of Water Leak (m)
Measured Center
of Water Leak (m)
Measured Length
of Water Leak (cm)
Measured Total Capacity (% of First Measure)
21.9733.2-
44.1121.664.40
66.0230.4109.63
88.1419.877.54
1010.4162.2136.93
1212.2522.0109.95
1414.1437.5117.07
Table 2. Estimation results for the water leak position and entity estimation in the case of leak length 2 m.
Table 2. Estimation results for the water leak position and entity estimation in the case of leak length 2 m.
Nominal Position
of Water Leak (m)
Measured Position
of Water Leak (m)
Measured Length
of Water Leak (cm)
Measured Total Capacity (% of First Measure)
21.64132.20-
43.59114.5686.85
66.22183.06211.06
87.90241.84187.61
1010.60118.45121.66
1212.22124.75182.85
1414.18147.00115.07
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

Scarpetta, M.; Cataldo, A.; Spadavecchia, M.; Piuzzi, E.; Masciullo, A.; Giaquinto, N. Accurate Detection and Localization of Water Pipe Leaks through Model-Based TDR Inversion. Sensors 2023, 23, 710. https://doi.org/10.3390/s23020710

AMA Style

Scarpetta M, Cataldo A, Spadavecchia M, Piuzzi E, Masciullo A, Giaquinto N. Accurate Detection and Localization of Water Pipe Leaks through Model-Based TDR Inversion. Sensors. 2023; 23(2):710. https://doi.org/10.3390/s23020710

Chicago/Turabian Style

Scarpetta, Marco, Andrea Cataldo, Maurizio Spadavecchia, Emanuele Piuzzi, Antonio Masciullo, and Nicola Giaquinto. 2023. "Accurate Detection and Localization of Water Pipe Leaks through Model-Based TDR Inversion" Sensors 23, no. 2: 710. https://doi.org/10.3390/s23020710

APA Style

Scarpetta, M., Cataldo, A., Spadavecchia, M., Piuzzi, E., Masciullo, A., & Giaquinto, N. (2023). Accurate Detection and Localization of Water Pipe Leaks through Model-Based TDR Inversion. Sensors, 23(2), 710. https://doi.org/10.3390/s23020710

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