## Abstract

This paper presents a new burst detection and location technique for pressurized pipelines based on an extension of the Differential Evolution (DE) algorithm. The proposed approach addresses the burst location problem as an optimization task, by considering the dynamic model that describes the behavior of a fluid through a pipeline and the presence of fluid losses produced by a burst. The optimization problem relies on finding suitable estimations related to the burst parameters, i.e. magnitude, pressure and position of a burst, while a defined cost function is minimized. In order to deal with this problem, three strategies are proposed to extend and adapt the DE algorithm: (i) an informed definition of the physical restrictions of the problem according to the pipeline characteristics; (ii) a training stage of the algorithm that allows to find the appropriate synthesis parameters; (iii) a multi-start structure, in order to track dynamical variations of the problem. Experiments on a pipeline prototype illustrate the results obtained by the proposed algorithm on the estimation of the burst parameters, comparing its performance with an algorithm based on the Extended Kalman Filter, which is widely used in the literature.

## INTRODUCTION

Water loss management is one of the most important issues facing water suppliers around the world. In many water pipelines, a significant percentage of the total water volume supplied is lost due to bursts and leaks, resulting in discontinuous, expensive and inefficient services. Furthermore, the presence of a burst can lead to the rapid deterioration of the pipeline, which might allow contaminant ingress and affect the urban infrastructure (Xie *et al.* 2017). Therefore, in the current literature, different works are focused on developing leak and burst diagnosis algorithms. These algorithms have a common objective to detect and locate leaks or bursts automatically in transmission pipelines, with good precision and minimum invasion of the pipeline (e.g. Verde 2001; Wang *et al.* 2002; Kowalczuk & Gunawickrama 2004; Torres *et al.* 2014; Navarro *et al.* 2017; Brunone *et al.* 2018; Lizarraga-Raygoza *et al.* 2018). In general, there are two widely used approaches to address the burst diagnosis problem, the Transient Test Based (TTB) and the Fault Model Approach (FMA). The idea of the TTB is to exploit the information of the transient produced by burst (Wang *et al.* 2002; Brunone *et al.* 2015, 2018), where the governing equations are expressed in terms of a Fourier series, and it is distinguished the deviation of the damping characteristics to the burst occurrence from the damping characteristics of the normal operation conditions. On the other hand, the FMA is based on the use of an analytical representation of the dynamic behavior of the flow rate and the pressure head in a pipeline, where the effect of a possible burst that affects the pipeline dynamics is considered within the model. This makes possible a direct estimation of the burst parameters (Santos-Ruiz *et al.* 2018). The extended Differential Evolution (DE) algorithm proposed to address the Burst Diagnosis and Location (BDL) problem in the present paper is focused in the FMA.

In general, burst detection and localization algorithms based on the FMA are developed under different design considerations and mathematical assumptions. However, in real applications, these considerations are not always fulfilled, causing a reduction in the algorithms reliability and precision. Then, it is convenient to use algorithms developed under different estimation techniques that can work in parallel (analytical redundancy) in order to offer a more informed and reliable diagnosis (Carvajal-Rubio & Begovich 2016; Santos-Ruiz *et al.* 2018). This fact motivates the development of new techniques to address the burst detection and location problem in transmission pipelines. In view of the above, the main contribution of this paper is *the development of a new BDL algorithm based on an extended DE algorithm and under FMA, which presents characteristics such as reliability, easy implementation and simple configuration*. As far as we know, the use of heuristic techniques to address the BDL problem has not been explored before.

The DE algorithm has been widely applied in the water management and optimal design of water distribution systems (Suribabu 2010; Zheng *et al.* 2011). In this paper, the DE algorithm is proposed to address the burst detection and localization problem as an optimization task that relies on finding suitable estimations related to the magnitude, pressure and position of a burst, while a proposed cost function is minimized. This algorithm addresses the optimization problem through a sequence of evolutive operations (mutation, recombination and selection) in order to find the candidate solution with the best cost function evaluation.

In order to achieve feasible solutions with a low number of iterations, the present paper also proposes a three-step procedure to improve the adaptation of the classical DE algorithm in the BDL problem. These steps are: first, the boundaries of the search space are carefully proposed as a set of restrictions of the optimization problem according to the physical characteristics and mechanical properties of the pipeline. The boundaries of the search space play a fundamental role in the performance of the DE algorithm. Thus, by providing informed bounds, as a priori knowledge, the DE algorithm can reach a good feasible solution in a short time. Second, to select the appropriate synthesis parameters for the DE algorithm, a training step is performed using the method presented by Pérez-González *et al.* (2019). This method consists in evaluating the performance of the algorithm in a particular case, under an iterative execution that uses different synthesis parameters. The information acquired from the iterative execution allows to find the values of the parameter settings where few iterations are needed to approximate a solution to the problem. On the other hand, the BDL problem as an optimization task requires a special consideration due to its dynamic characteristics. Therefore, the proposed algorithm includes, as its third strategy, a multi-start structure that allows to track the pressure dynamical variations in the burst point. The proposed burst diagnosis method based on the extended DE algorithm is implemented and evaluated in a pipeline prototype installed at Cinvestav, Guadalajara, México. As it will be shown, the obtained results illustrate the effectiveness of the proposed algorithm compared with the Extended Kalman Filter (EKF), which is widely used for the burst diagnosis problem.

The paper is organized as follows: Preliminaries section presents the pipeline dynamic model, introduces the BDL problem, and describes the classical DE algorithm. In the Burst parameter identification strategy section, the improved DE algorithm to tackle the BDL problem solution is presented. The Experiments and results section presents and discusses the results obtained. Finally, the respective conclusions are stated.

## PRELIMINARIES

This section introduces the burst detection and localization problem, starting with a brief description of the pipeline fluid dynamic equations, and afterwards the model considered in this work to address the BDL problem is introduced. Also, a brief description of the classical DE algorithm is presented.

### Pipeline model

*et al.*1998). In general, the transient flow is considered unidirectional since the terms of the convective acceleration are smaller than the axial terms. Therefore, by neglecting these terms and the spatial density variation, besides considering a fluid slightly compressible, pipes slightly deformable, and a leveled pipeline, the WHE can be expressed as (Chaudhry 2014):where

*a*is the pressure wave speed in the fluid [],

*A*is the cross-sectional area of the pipeline [],

*Q*is the flow rate [],

*z*is the spatial coordinate along the pipeline [],

*g*is the gravitational acceleration [],

*H*represents the pressure head [],

*t*is the temporal coordinate [], represents the friction factor head losses [], and

*D*is the diameter of the pipe [].

The pressure wave speed can be approximated by different formulations; however, the precision of the estimations highly depends on the mechanical properties of the material, which change with the pipeline deterioration. Furthermore, in plastic materials, the pressure wave speed estimations present variations in comparison with the values suggested in classical textbooks (Evangelista *et al.* 2015). Therefore, different authors as Meniconi *et al.* (2015) and Soares *et al.* (2008) propose experimental procedures in order to obtain more suitable pressure wave speed estimations.

*et al.*1995; Duan

*et al.*2017) have been proposed to compute the unsteady part . In this way, Dulhoste

*et al.*(2011, 2017) present a study that evaluates the influence of the friction unstable term in the implementation of leak location algorithms based on observers and developed under the FMA. Based on their results, the authors conclude that the unsteady formulation is complex and it only presents a small improvement in the burst diagnosis. In view of the above, the present paper neglects the unsteady friction term.

### Burst localization model

In order to develop a model based on FMA, it is necessary to consider the burst effect into a model. When a pipe burst occurs, a sudden pressure wave is generated and propagated in both directions from the burst point. This behavior induces spatial and temporal changes in the velocity flow rate and the pressure, i.e. transients, until a new stable state, is reached. This new stable state depends on the burst size and magnitude. A widely used model in practical and real scenarios (e.g. Torres *et al.* 2014; Delgado-Aguiñaga & Begovich 2017; Verde & Torres 2017; Santos-Ruiz *et al.* 2018), that takes into account the burst effects, is deduced from the spatial discretization of the WHE through the explicit Finite Differences method and then incorporating the effect of the burst.

*et al.*2014):where is a factor associated with the burst magnitude, is the orifice area [], is the pressure head at the burst point [] and is the discharge coefficient [].

*et al.*(2015), as follows:whereand

*k*is the index of the discrete time and is the time step. Model (8) can be written in compact form as follows:where , and is a nonlinear function.

In a practical scenario, the typical instrumentation of a transmission pipeline is composed of pressure and flow rate sensors located at the ends of the pipeline. The main objective of this sensor configuration is to evaluate the supplied volumes to the distribution zones and to monitor the efficiency of the system. This pipeline instrumentation scheme is suitable for the use of model (9) in the implementation of a BDL algorithm. Also, it is necessary that the sampling rate satisfies the Courant–Friedrichs–Lewy (CFL) criterion in order to guarantee the numerical stability of model (9). The CFL condition can be expressed as: , where Δ*t* is the time step.

### Longitudinal adjustment

*local loss*, which is added to the friction produced in the straight pipe sections. The WHE are developed considering a straight pipeline, without accessories. Thus, in order to obtain a dynamical model that includes fittings, in terms of the WHE, it is necessary to express the length of the pipeline in an

*equivalent straight length*(). Hence, each accessory in the studied pipeline will be seen as a virtual straight pipe section that generates the same local loss associated with its corresponding accessory (Mataix 1982). Following this, is computed by using the Darcy–Weisbach equation:where is the differential head pressure at the ends of the pipeline in a burst-free condition []. Finally, by letting in Equation (8), it is possible to introduce the effects of the fittings into the pipeline model.

Using the longitudinal compensation, the burst is localized in a virtual equivalent straight length, however, in a practical situation, this does not represent a complete solution. In addition to the burst localization, it is necessary to establish a conversion from the obtained virtual position to a real-length value. In this work, we use the procedure described by Badillo-Olvera *et al.* (2017) in order to perform this conversion.

### Differential Evolution algorithm

The DE is a powerful optimization algorithm, which has proved to be very robust in a wide variety of problems (Price *et al.* 2006; Coello *et al.* 2007); furthermore, it only needs two control synthesis parameters, which makes it easy to configure. In general, the DE algorithm is used to address problems in which the decision variables are real numbers and this feature is suitable to address the BDL problem. In the DE algorithm, the candidate solutions play the role of individuals, being classified according to an objective function . This classification ensures the survival of the elements with the best performance, driving the evolution in a set of greedy directions.

The DE algorithm maintains the general scheme of an evolutionary method considering a population, of the form , where is an integer index that distinguishes each individual from the other members of the -size population, *g* is the generation index and *d* is the dimension of each individual, given by the number of parameters to be optimized, that is, each element of each individual corresponds to an unknown parameter of the optimization problem. The initial population is generated by a uniform random process, settled in the interval {} for . After the initialization, the DE process becomes an iterative sequence of evolutionary operations based on mutation, recombination, and selection (Zhang & Sanderson 2009). This process will be described as follows.

*Mutation:*In this step, mutant vectors are calculated at each*g*generation, based on the current population . In this work, the mutation strategy denominated by Price*et al.*(2006) as*DE current-to-best/1*is used. This technique considers the candidate solution with the best performance found at the moment, , as part of the mutation operation, in addition to a randomized contribution. The equation that describes this process is as follows:where and are positive integers from a uniform random distribution in the interval [],*F*is the mutation factor, which usually ranges in the interval [0,2]. In general, the ability to perform local searches (exploitation) is achieved through values of*F*near to 0 and, on the contrary, the global search (exploration) is given by fixing*F*near to 2.*Recombination:*In order to improve the diversity of the population, a recombination process is defined after the mutation step, performing a crossover operation in each individual , by using as its partner the mutant vector . This operation generates a trial vector by means of:where is the*crossover rate*, which is a constant fixed in the interval [], and it controls the admittance of the entries from the mutant vector in the formation of the trial vector.

The processes of mutation, recombination and selection are iteratively executed until the optimum parameters are obtained or a pre-specified stop criterion is satisfied.

## BURST PARAMETER IDENTIFICATION STRATEGY

The general formulation of identification techniques is based on the comparison between the response of a real system *y* (sensor measurements), against the outputs given by a related mathematical model , while a performance criterion that quantifies the value of the approximation error is used. Therefore, the optimization task relies on finding a set of values of the burst parameters minimizing a criterion function of the estimation error.

*y*, for each candidate solution. So, model (9) is completed by using the values obtained by the proposed algorithm and then, the performance of each is evaluated by the cost function (14) in order to select the best candidate solution:where and , in order to select the best candidate solution.

The following subsections present three strategies that allow to extend the classical DE algorithm in order to tackle the BDL problem efficiently.

### Initial bounds of the search space

Before the population of the DE algorithm can be initialized, both upper and lower bounds for all the candidate solutions must be specified. In this case, for each characteristic of the burst, there is a certain range within which the value of the bounds should be restricted, due to its physical correspondence. A good definition of the feasible search space and their boundaries provides extra knowledge about the optimization problem, making the deployment of the algorithm faster and efficient (avoiding the consideration of unfeasible candidate solutions). The BDL problem is defined by the unknown values of the burst position, pressure and magnitude of model (9), whose physical boundaries are included in Table 1.

Parameter . | Symbol . | Upper bound . | Lower bound . | Units . |
---|---|---|---|---|

Burst position | [] | |||

Burst pressure | [] | |||

Magnitude factor | [] |

Parameter . | Symbol . | Upper bound . | Lower bound . | Units . |
---|---|---|---|---|

Burst position | [] | |||

Burst pressure | [] | |||

Magnitude factor | [] |

*Leak-Before-Break*(LBB) criterion. A basic analysis, based on Linear Elastic Fracture Mechanics (LEFM), to establish the LBB criterion is presented by Anderson (2017). Considering this, the crack size will remain stable if:where is the fracture toughness of the pipe material [], and is the operating hoop stress [], calculated as , where

*P*is the operating pressure [].

### Calibration zone

One of the main difficulties that a user faces when a heuristic algorithm is applied is to select an appropriate set of synthesis parameters (Lobo *et al.* 2007). In the case of the DE algorithm, the main synthesis parameters that are required are the *Mutation Factor**F* and the , whose selection is crucial for a correct algorithm performance.

*et al.*2016)

*Comfort Zone*(CZ). This methodology consists in evaluating the performance of the DE algorithm, in a particular case, under an iterative execution that uses different parameter settings , selected by means of the expression:where the lower bounds in both cases are 0.0001 in order to guarantee numerical stability and the upper bounds of

*F*and are 2 and 1, respectively, in agreement with Price

*et al.*(2006). The information obtained from the CZ allows a relationship to be established between the calibration and performance of an algorithm, providing a statistical measurement of the probability to obtain a high-performance execution by using a calibration setting that lies in a specific region. Of course, the intention is to use the calibration settings that present the higher probabilities of providing good performance for the algorithm. In this case, the number of iterations needed to reach a specific error bound is defined as a performance measure of the algorithm.

### Multi-start scheme

In Equation (9), the pressure at the burst point corresponds to a *state* of the model, which turns the optimization task into a dynamic optimization problem. That is, the proposed algorithm must track the dynamic behavior of at each instant of the sampling time. In order to do that, the use of a multi-start scheme is proposed (Branke 2012), that allows consideration of the evaluation of model (9), in each sampling time *k*, as an independent optimization problem. So, for the -th sample, there exists a set of values that minimizes (14), therefore, the DE algorithm estimates these values. Before the dataset that corresponds to the sampling time arrives, the DE algorithm is initialized and iterated in order to solve the optimization problem represented by the dataset sampled in the -th interval of time. With this solution, the arrival of the sample set is allowed and the DE algorithm starts again over it, taking this new configuration as a new optimization problem. Of course, tackling a single optimization problem in each sampling time is computationally expensive, nevertheless, the set CZ offers a statistical guarantee of reaching a good solution in a low number of iterations, which makes it easier to execute the optimization scheme as a set of independent optimization tasks.

Finally, the proposed DE algorithm to address the BDL problem can be illustrated by the pseudocode presented in the Appendix section (available with the online version of this paper).

## EXPERIMENTS AND RESULTS

### Prototype description

The BDL methodology based on the extended DE algorithm is evaluated using datasets acquired from the prototype built at the Center for Research and Advanced Studies (Cinvestav), located in Guadalajara, México. Figure 2 presents the prototype diagram and Table 2 contains the physical parameters of the prototype. This pipeline consists of a Polypropylene Random Copolymer circuit of pipes, supplied by a water tank (T1) through a centrifugal pump (FB1), which is controlled by a variable frequency driver of . The flow rate and pressure head at the ends of the pipeline are measured by means of two ultrasonic flow sensors *Promag Propiline 10P* (FT1 and FT2) and two pressure sensors *PMP 41* (PT1 and PT2), respectively. The temperature is measured by means of a resistive sensor *PT100* (TT1). The prototype also has one electrovalve and a manual ball-type valve, whose diameter is 0.0127 m, which allow to emulate a burst in their locations: (valve 1) and (valve 2). Each valve also has a pressure sensor (PT3 and PT4) that allows validation of the pressure estimations in the points of each emulated burst. Pictures of some components of the described prototype can be found in https://www.gdl.cinvestav.mx/ofelia/index.php?page=fdg (accessed March 17, 2019).

Parameter . | Symbol . | Value . | Units . |
---|---|---|---|

Length between sensors | 68.147 | [] | |

Roughness | 7.0 × 10^{−6} | [] | |

Internal diameter | 6.271 × 10^{−2} | [] | |

Wall thickness | 13.095 × 10^{−3} | [] | |

Young modulus | 7.008 × 10^{8} | [] | |

Slope | 0 | [] |

Parameter . | Symbol . | Value . | Units . |
---|---|---|---|

Length between sensors | 68.147 | [] | |

Roughness | 7.0 × 10^{−6} | [] | |

Internal diameter | 6.271 × 10^{−2} | [] | |

Wall thickness | 13.095 × 10^{−3} | [] | |

Young modulus | 7.008 × 10^{8} | [] | |

Slope | 0 | [] |

In this research, the data acquisition is performed by a *NI-myRIO*-embedded device. This system is equipped with a *Xilinx Zynq-7010 FPGA* and a dual-core ARM CortexTM processor. The computational processing of the acquired data is carried out by means of a desktop computer equipped with an IntelCore^{TM} i7-2600 K, processor, of , and Windows^{TM} 10 PRO edition OS, including LabView^{TM} 2016 myRIO-oriented software.

In the diagram, the sensors FT1 and FT2 acquire the flow measurements and , respectively, and the sensors PT1 and PT2 acquire the pressure measurements and , respectively.

### Setup and settings

*et al.*(2015), Delgado-Aguiñaga & Begovich (2017), and Lizarraga-Raygoza

*et al.*(2018), in this work attempting to find improved results, the pressure wave speed is experimentally determined as in Pérez-González

*et al.*(2019). In the experiment, a burst is produced by a valve with known position and then the difference between instants of time where the pressure waves arrive at the pressure sensors located at the ends of the pipeline are registered. Once the time instants are obtained, the wave speed is estimated as:where is the distance between pressure sensors, is the known valve position,

*V*is the fluid velocity and is the difference of times and in which the pressure drop, generated by the burst, arrives at the upstream and downstream sensors, respectively. In order to estimate the arrival time of the pressure wave to the sensors, a Discrete Wavelet Transform (DWT) method is used as in Romero-Delgado & Begovich (2017). With the only purpose to calibrate the pressure wave speed, a burst is produced by means of the valve 1 (located at 17.045 m) under the experimental conditions of the Table 3. As result of this experiment, a pressure wave speed value is obtained.

Steady flow [] . | Temperature [] . | . | . | . | . | . |
---|---|---|---|---|---|---|

8.892 | 24.31 | 68.147 | 2.8906 | 17.045 | 30.540 | 30.631 |

Steady flow [] . | Temperature [] . | . | . | . | . | . |
---|---|---|---|---|---|---|

8.892 | 24.31 | 68.147 | 2.8906 | 17.045 | 30.540 | 30.631 |

*F*and of the algorithm to address the BDL problem, the CZ procedure is performed. In this way, the algorithm is executed in an iterative way on the acquired data from the pipeline prototype by producing a burst by the valve 1. To define the CZ, a cost function based on the total mean square error is considered for each execution:where

*K*is the total number of samples in a database. It is expected that the proposed algorithm reaches a good estimation of the burst parameters at the end of each execution, allowing a maximum number of 500 iterations in each execution. In addition to this, an execution is considered

*successful*when the algorithm reaches an error bound defined as in less than 20 iterations, as

*poor execution*when the algorithm reaches the error bound in more than 20 iterations, but less than 500 and as

*failed execution*when the algorithm does not reach the error bound, even with 500 iterations. The results obtained after 1000 executions are shown in Figure 3. The circles correspond to the successful executions and the squares the failed executions. The successful executions represent 18% of the total of the executions. As Figure 3 shows, when the parameters

*F*and are specified arbitrarily, there is a great probability to obtain a bad performance. On the other hand, it can also be seen that there is an interval of

*F*and where the ED algorithm presents a good performance.

The successful executions are distributed between the regions: and . The minimum value of iterations needed to reach the error bound is 14, obtained by using and .

Once the CZ is determined, the parameters setting for the implementation of the DE algorithm is taken as , since for many applications, the population size is taken as , as it is reported in Price (1997). The control parameters , obtained from the CZ are used, and a total of 20 iterations per sample are considered. Finally, in order to detect the presence of a burst, and in consequence, start the algorithm, the burst detection alarm is established when the expression | becomes true, in order to avoid the uncertainty in the flow measurement introduced by the noise level of the sensors.

It is important to point out that the model-based burst diagnosis methods present high dependence to some pipeline parameters, such as the friction and diameter, then it is necessary to periodically calibrate the pipeline parameter.

## RESULTS AND DISCUSSION

In order to show the performance of the proposed algorithm in the BDL problem, this subsection presents two different experiments using real data from the prototype described above. In addition to this, the proposed algorithm is compared against the EKF, as one of the most widely used algorithms that presents good performance in the BDL problem (Delgado-Aguiñaga *et al.* 2016).

In *Experiment 1*, a single burst is emulated under constant flow conditions, by means of the valve 2, located at at time . The mean temperature during the experiment is 24.31°C and the mean viscosity is . The pressure head and flow rate at the ends of the pipeline are shown in Figure 4 and the experimental conditions are presented in Table 4.

Experiment parameter/signal . | Symbol . | Value before the burst occurrence . | Value after the burst occurrence . | Units . |
---|---|---|---|---|

Input pressure head | [] | |||

Output pressure head | [] | |||

Input flow rate | [] | |||

Output flow rate | [] | |||

Reynolds number | See text below | [] | ||

Friction factor | See text below | [] |

Experiment parameter/signal . | Symbol . | Value before the burst occurrence . | Value after the burst occurrence . | Units . |
---|---|---|---|---|

Input pressure head | [] | |||

Output pressure head | [] | |||

Input flow rate | [] | |||

Output flow rate | [] | |||

Reynolds number | See text below | [] | ||

Friction factor | See text below | [] |

After the burst occurrences, two friction factors and Reynolds numbers are estimated with the input and output flows. In this experiment, and . On the other hand, we have and .

The results obtained from Experiment 1 with the extended DE algorithm and with the EKF are shown in Figure 5: (a) presents the real burst position and its estimation with the proposed DE algorithm and the EKF; the real pressure head and its estimation are shown in (b), and finally, the magnitude factor and its estimation are presented in (c). As it is shown in Figure 5, the burst value estimations obtained from the extended DE algorithm and the EKF are close to the real values.

In order to validate the obtained estimations, Figure 6 shows the comparison between the measured flows and their estimations given by model (9), obtained by the incorporation of the best parameters found by the extended DE algorithm at each sample.

Table 5 presents the average, root mean square error (RMSE) and standard deviation (STD) of each estimated value obtained with the extended DE and EKF algorithms, considering a time interval between 130 s and 300 s.

Method . | Parameter . | Real value . | Average estimation . | Max. STD . | RMSE . |
---|---|---|---|---|---|

Differential Evolution | 49.895 | 47.856 | 1.94 | 2.474 | |

11.147 | 12.107 | 0.162 | 0.975 | ||

1.65 × 10^{−4} | 1.49 × 10^{−4} | 1.17 × 10^{−6} | 1.63 × 10^{−5} | ||

Extended Kalman Filter | 49.895 | 47.306 | 0.813 | 2.277 | |

11.147 | 12.222 | 0.0932 | 1.060 | ||

1.65 × 10^{−4} | 1.56 × 10^{−4} | 7.85 × 10^{−6} | 7.37 × 10^{−6} |

Method . | Parameter . | Real value . | Average estimation . | Max. STD . | RMSE . |
---|---|---|---|---|---|

Differential Evolution | 49.895 | 47.856 | 1.94 | 2.474 | |

11.147 | 12.107 | 0.162 | 0.975 | ||

1.65 × 10^{−4} | 1.49 × 10^{−4} | 1.17 × 10^{−6} | 1.63 × 10^{−5} | ||

Extended Kalman Filter | 49.895 | 47.306 | 0.813 | 2.277 | |

11.147 | 12.222 | 0.0932 | 1.060 | ||

1.65 × 10^{−4} | 1.56 × 10^{−4} | 7.85 × 10^{−6} | 7.37 × 10^{−6} |

A *second experiment* is performed, in order to highlight the effectiveness of the proposed method. A new burst is emulated using again the valve 2 (49.895 m) at time , but also, producing variations in the operation point of the pump, by means of the frequency driver. Figure 7 shows the dataset used during this new experiment and Table 6 presents the experimental scenario.

Experiment parameter/signal . | Symbol . | Value before the burst occurrence . | Value after the burst occurrence (variation range) . | Units . |
---|---|---|---|---|

Input pressure head | 14.116–19.9137 | [] | ||

Output pressure head | 7.0705–10.3411 | [] | ||

Input flow rate | 0.0077–0.0093 | [] | ||

Output flow rate | 0.0073–0.0088 | [] | ||

Reynolds number | See text below | [] | ||

Friction factor | See text below | [] |

Experiment parameter/signal . | Symbol . | Value before the burst occurrence . | Value after the burst occurrence (variation range) . | Units . |
---|---|---|---|---|

Input pressure head | 14.116–19.9137 | [] | ||

Output pressure head | 7.0705–10.3411 | [] | ||

Input flow rate | 0.0077–0.0093 | [] | ||

Output flow rate | 0.0073–0.0088 | [] | ||

Reynolds number | See text below | [] | ||

Friction factor | See text below | [] |

In this experiment, the variation range of the friction factors and are and , respectively. On the other hand, the variation ranges of the Reynolds numbers and are [] and [2.3701 × 10^{5} − 1.9610 × 10^{5}], respectively.

Figure 8 shows the estimation results during the second experiment. The plot (a) presents the estimation of , where the estimated position of the burst point does not change in spite of the variations in the operating point of the pump. In this way, the estimation of maintains its behavior in the neighborhood of the real value with a displacement of 1.45%. The same behavior is presented by the magnitude factor, presenting a displacement of 8.23% with respect to the real value of as Figure 8(c) shows. The case of the pressure head at the leak point is different from the previous two, since changes with respect to the variations of the operation point, despite this, the proposed extended DE algorithm can follow the dynamic variations.

In order to validate the obtained estimations, Figure 9 shows the comparison between the measured flows and their estimations given by Equation (9) considering the best parameters found by the extended DE algorithm at each sample.

Table 7 summarizes the obtained results and their corresponding statistical values, considering a time interval between 130 s and 300 s.

Method . | Parameter . | Real value . | Average estimation . | max STD . | RMSE . |
---|---|---|---|---|---|

Differential Evolution | 49.895 | 49.386 | 0.610 | 2.25 | |

See text below | See text below | 0.430 | 0.722 | ||

1.32 × 10^{−4} | 1.431 × 10^{−4} | 1.19 × 10^{−4} | 1.255 × 10^{−5} | ||

Extended Kalman Filter | 49.895 | 51.863 | 1.016 | 3.224 | |

See text below | See text below | 0.387 | 0.523 | ||

1.327 × 10^{−4} | 1.316 × 10^{−4} | 5.793 × 10^{−6} | 3.201 × 10^{−6} |

Method . | Parameter . | Real value . | Average estimation . | max STD . | RMSE . |
---|---|---|---|---|---|

Differential Evolution | 49.895 | 49.386 | 0.610 | 2.25 | |

See text below | See text below | 0.430 | 0.722 | ||

1.32 × 10^{−4} | 1.431 × 10^{−4} | 1.19 × 10^{−4} | 1.255 × 10^{−5} | ||

Extended Kalman Filter | 49.895 | 51.863 | 1.016 | 3.224 | |

See text below | See text below | 0.387 | 0.523 | ||

1.327 × 10^{−4} | 1.316 × 10^{−4} | 5.793 × 10^{−6} | 3.201 × 10^{−6} |

Remember that in Experiment 2, the head pressure after the burst occurrence is not constant and considering this condition the average value is not expressed. However, in Figure 8(c), it is possible to see that the pressure variations are well estimated.

The results obtained from the performed experiments show that the proposed extended DE algorithm offers similar results than the EKF. Nevertheless, in the first experiment, the burst position presents a lower STD with the EKF, in contrast, in the second experiment, the STD is less with the DE algorithm. This fact shows that the data fusion of multiple algorithms could reduce the uncertainty degree provided by an individual estimation.

The extended DE algorithm presents additional advantages since it requires only a few samples to approximate a solution and it does not need a linearization of the mathematical model, as the EKF does, simplifying in this way the implementation. However, like most evolutive algorithms, its real-time implementation is limited by the computational cost of each iteration in relation to the available time between samples, in this case, the mean computational time per iteration of the ED algorithm is 0.00366 s, while the EKF only takes . Furthermore, it is necessary to emphasize that the correct definition of the search space and the adequate selection of the synthesis parameters of the DE algorithm allow to reach a good estimation of the burst values and also the multi-start scheme provides a powerful tool in order to deal with dynamic optimization.

## CONCLUSIONS

In this paper, an extended version of the DE algorithm is proposed to address the burst detection and localization problem, where the problem relies on estimating the parameters related to a burst occurrence. The experimental results obtained from the used pipeline prototype show that the proposed DE algorithm is a suitable method to estimate the burst parameters under different practical scenarios. In addition to this, it is important to point out that the correct definition of the search space, as well as the synthesis parameters of the DE algorithm, play a fundamental role to achieve an adequate estimation of the burst parameters. Moreover, the estimation of the burst parameters provided by the DE algorithm is similar to the ones provided by the EKF. This fact allows to consider the DE algorithm as a useful element for the conformation of a bank of burst diagnosers, as a strategy that provides more informed and reliable burst detection and localization. As future work, the heuristic techniques can be included in more complex BDL problems, such as multiple burst diagnosis and branched pipes.

## REFERENCES

**119**, 243–252.

**15**(8), 760–868.

*13th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE)*

*23rd Mediterranean Conference on Control and Automation (MED)*

*50th IEEE Conference on Decision and Control and European Control Conference*

*Paper presented at the Proceedings of the Latin American Conference on Automatic Control, XVII CLCA*

**141**(5), in press.

*IEEE International Conference on Evolutionary Computation*

*14th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE)*

**20**(6), 1286–1295.

783(1), 012012.