OJOp  Vol.7 No.2 , June 2018
Optimal Water Pipe Replacement Policy
ABSTRACT
Water scarcity is the major problem confronting both urban and rural dwellers in Enugu State. This scarcity emanated from indiscriminate pipe failure, lack of adequate maintenance, uncertainty on the time of repair or replacement of pipes etc. There is no systematic approach to determining replacement or repair time of the pipes. Hence, the rule of thumb is used in making such a vital decision. The population is increasing, houses are built but the network is not expanded and the existing ones that were installed for no less than two to three decades ago are not maintained. These compounded the problem of scarcity of water in the state. Replacement or repair of water pipes when they are seen spilling water cannot solve this lingering problem. The solution can be achieved by developing an adequate predictive model for water pipe replacement. Hence, this research is aimed at providing a solution to this problem of water scarcity by suggesting a policy that will be used for better planning. The interests in this paper were to obtain a water pipe failure model, the intensity function λ(t) [failure rate], the reliability R(t) and the optimal time of replacement and they were achieved. It was observed that the failure rate of the pipes increases with time while their reliability deteriorates with time. Hence, the Optimal replacement policy is that each pipe should be replaced after 4th break when the reliability = 0.0011.

1. Introduction

Water distribution network is a network consisting of pipes, valves and pumps that distribute water to the end users [1] [2] . In water distribution some components are made redundant which makes water supply failure not to be noticed immediately when pipe failure occurs. Pipes are of different types and sizes. In this paper, our interest is on three types of pipes namely:

1) Asbestos-concrete (AC) pipes,

2) Unplastisized polyvinyl chloride (UPVC) or plastic pipes and

3) Galvanized Iron pipes (GVI).

These pipes are further divided into mains, branches and terminals. Mains are of different lengths and diameters ranging from 13 ft to 20 ft and 200 mm to 600 mm. They carry water directly from the pumping station and the distribution storage. Branches have a diameter ranging from 75 mm to 150 mm, but of the same length with the mains. They carry water from the mains to various terminals. Terminals also have the same length as others but have a diameter ranging from 25 mm to 50 mm. They are called service lines in the network. The terminals convey water to the end-users and they are the smallest in size but the largest in number. Table 1 shows the type of water distribution pipes.

Each of the types of pipes in Table 1 has its advantages and shortcomings as presented in Table 2.

The major problems of the GVI pipes are bulkiness, rusting and expensiveness compared to the other two types. Majority of their failure result from rusting. UPVC is fast replacing GVI pipes because of their portability and durability.

Many factors are responsible for water pipe breaks. For the purpose of this paper, these factors are divided into internal and external factors, which are collectively called covariates. Internal factors: Internal factors include the pressure exerted by the flowing water in the pipe; the action of chemicals used in the water treatment; the diameter of the pipes and the intrinsic material of the pipes etc. Pressure exerted by water over a long period of time on pipes leads to their break and this pressure is inversely related to the diameter of pipes. The pressure

Table 1. Types of water distribution pipes.

Table 2. Water distribution pipes, advantages and disadvantages.

is also intensified, as the pipes are located nearer to the source and reduces as the length of the pipes increases. External factors: External factors include the soil type, topography, lay-dept, shock and rust. These factors are called external because they occur outside the pipes. The nature of soil in which the pipe is laid determines the rate of pipe’s break etc. The rates at which these factors occur influence the maintenance decisions.

Enugu State Water Corporation Enugu

Water generation and distribution are the sole concern of the water Generation unit, Enugu State Water Corporation Enugu. The replacement and repair functions are the sole responsibilities of the engineering and plumbing units. The pressure under which the pumps distribute water ranges from 6 bars to 16 bars. Frequency of pipe breaks depend on the following factors namely; lay depth, topography, intrinsic material of the pipes, pressure, diameter of the pipes, lay location etc. Period of replacement: pipes are replaced as they break and the rule of thumb is extensively used. Maintenance actions include periodic greasing of nuts and bolts, replacement or repair of failed pipes. In most cases, the maintenances are done on demand, as there is no regularity on the time limit for the maintenance services. But since pipes constitute the bulk of the money cost of the water network distribution system, there is a need to develop a failure model to detect the probability of failure at any given time, the failure rate, the reliability of the pipes and the optimal replacement time of the pipes to avoid waste of water resources.

[1] [2] used non-homogeneous Poisson process to model the pipe failure. They were of opinion that power law process model (λ(t) = abtb−1) can be used to describe the intensity of a NHPP. They define intensity as the mean number of failure, N(t), per unit time and that once the intensity has been estimated the number of failures in a given time period is Poisson.

[3] [4] observed that if the parameter b is between 0 and 1, then the network is improving, but if b > 1, it shows a deteriorating network. The model parameters can be determined using either graphical or analytical analysis. We can linearized the mean value function M(t) by taking the logarithm of both sides to obtained ln M ( t ) = ln ( a ) + b ln ( t ) , and a plot lnM(t) versus ln(t) should be a straight line with intercept ln(a) and slop b.

Water utilities are concerned with large number of main breaks and the resulting direct and indirect costs due to deterioration. These contribute to frequent breaks and resulting repair/replacement and rehabilitation costs [5] . They observed that replacement represents large capital investment. They suggested that if an estimator that can track the replacement time will be developed; such an estimator will be of great use to practicing engineers. [6] did a work on optimal replacement policy for stochastically failing equipment inaccessible to inspection. The policy was characterized by a single parameter, N. If equipment age is less than N, the appropriate action is to do nothing; if equal to N, the appropriate action is to replace the equipment. [7] obtained the optimal replacement policy by minimizing the average annual cost of equipment whose maintenance cost is a function increasing with time and whose scrap value is known, and their calculations were based on average annual cost incurred on the item. They assert that the item should be replaced when the average annual cost is minimal.

Water distribution pipes are repairable system because some maintenance action of repair can be done on them to restore the operation of the pipes after failures. [8] [9] [10] observed that a repairable system is a system which, after failing to perform its functions satisfactorily, can be restored to fully satisfactory performance by any method other than replacement of the entire system. Repairable systems shared the property that they could be repaired by replacing some failed components and returned to regular operation. Sometimes, minimal repairs bring the system reliability back to the same level it had before the failure, while ‘perfect’ repairs bring the reliability back to the state of the system at the start of the operation [11] [12] . The intensity and mean value functions of a power law process PLP are λ(t; M, β) = MβtB−1 and V(t; M, β) = Mtβ; M, β > 0. β determines how reliability decays or grows. If β < 1, the intensity function decreases over time and the reliability improves. If β = 1, the reliability remains the same over time, but if β > 1, the reliability declines and intensity increases. Again, if β > 1, we have NHPP and for β = 1, we have HPP. λ(t) is an increasing concave (straight, convex) curve if 1 < β < 2 [13] [14] [15] .

2. Overview of Poisson Process/Reliability

The Poisson Process provides a realistic model for many random phenomena. Since values of a Poisson random variable are the non-negative integers, any random phenomenon for which count of some sort is of interest is a candidate for modelling by assuming a Poisson distribution [16] , thus, the first event counted by Ns(.) takes place after an exponential amount of time with parameter λ; the rest of the inter event times are iid exponential with parameter λ [17] . A Poisson process is a renewal process for which the underlying distribution is exponential. A renewal process is a sequence of independent identically distributed non-negative random variables, X1, X2, …, which, with probability 1, are not all zeros [18] . They assert that the Poisson process arises quite naturally in reliability and life testing situations when the underlying life distribution is exponential. [19] observed that the original use of the term “reliability” was purely qualitative and that was why aerospace engineers recognized the need to have more than one engine on an airplane and drivers keep spare tyres in their vehicles without any precise measurement of the failure rate. But [20] were of opinion that reliability is the probability that a component or a system will perform its intended function under specified environmental conditions with time. They assert that the problem of assuring and maintaining reliability has many facets, which include: original equipment design; control of quality of products during production, life testing and design modifications.

Data Presentation

In this section, we present the data that is used for this work in Table 3.

3. Data Analysis

t = 2005 - 2015 = 10 years (1)

λ = x i N = 18 20 = 0.9 ; i = 1 , 2 , , 18 ; N = 1 , 2 , , 20 (2)

λ c = 0.9 10 = 0.09 ; λ c = failureperunittime (3)

We present how the parameters of the model were determined in Table 4.

The graphic method of determining the parameters of the water pipe failure model is presented in Figure 1.

ln a = 2.4 a = exp ( 2.4 ) = 0.09 (4)

b = 1.60 1.575 = 1.02 (5)

Table 3. Type of pipe, Installation period, Replacement and number of Breaks.

Table 4. Determination of the Parameters of the Model.

Figure 1. Graph of lnL(t) versus ln(t).

Hence

p n ( t ) = exp ( 0.09 t 1.02 ) ( 0.09 t 1.02 n ! ) n ; ( n , t ) = 0 , 1 , 2 , , 10 (6)

p 0 ( 1 ) = exp ( 0.09 ) = 0.9139 (7)

p 1 ( 1 ) = exp ( 0.09 ) ( 0.09 ) = 0.0823 (8)

p 2 ( 2 ) = exp ( 0.1825 ) ( 0.1825 ) 2 2 ! = 0.0139 (9)

p 10 ( 10 ) = exp ( 0.9424 ) ( 0.9424 ) 10 10 ! = 0.000 (10)

3.1. Water Pipe Failure Model

Since b > 1, then, our model assumes NHPP (λ(t)).

Equation (6) is the water pipe failure model, where p = probability; n = the number of failures; t = the time in years; a = 0.09 and b = 1.02 are the parameters of the model and a, b, t > 0.

3.2. The Reliability Function R(t)

R ( t ) = p n ( t ) λ ( t ) (11)

3.3. Determination of Optimal Time of Replacement

A reliability approach is used here to determine the optimal time of replacement of pipes. Reliability approach has never been used to determine the optimal water pipe replacement policy. The pipe should be replaced a little time before the reliability becomes zero. Continuous used of the pipe when the reliability become zero add more to the cost because the failure rate λ(t) increases rapidly as the reliability R(t) becomes zero. Hence, t (time in years) is optimal at the least reliability R(t).

4. Interpretation of Results

Form the data collected from Enugu State Water Corporation Enugu we obtain the parameters of our model “a” and “b” and use them to generate a table for the probability of number of failure at any given time t. Also failure rate and reliability of the pipes were determined and presented in Tables 4-7 respectively.

Result A:

p n ( t ) = exp ( 0.09 t 1.02 ) ( 0.09 t 1.02 n ! ) n ; ( t , n ) = 1 , , 10 (12)

The probability of the number of failure of pipes at a given time is presented in Table 5 below.

Result B:

λ ( t ) = ( 0.0918 ) t 0.02 (13)

The computation of the intensity function λ(t) is presented in Table 6.

Table 6 shows that the intensity function, λ(t), increases with time.

Result C:

R ( t ) = p n ( t ) λ ( t ) (14)

In Table 7, we present the computation of reliability, R(t), of the water pipes.

Table 5. Probability of number of failure at Time t.

Table 6. Computation of the Intensity Function λ(t).

Table 7. Computation of the Reliability function R(t).

Table 7 above shows that the reliability R(t) of the pipes deteriorates with time while the failure rate of the pipes increases with time. The pipe should be replaced at the 4th break when the reliability of the pipes is 0.0011. Hence, it is most reasonable to replace the pipe when the reliability is very close to zero.

5. Conclusion

In this paper, we developed water pipe failure model and the optimal water pipe replacement time (policy) of the pipes. The water pipe failure model was a Non-homogenous Poisson process and reliability has never been used to determine optimal replacement policy on a repairable system. We also obtained the intensity function λ(t) [failure rate] of the pipes. From this work, it was observed that the failure rate of the pipes increases with time while their reliability deteriorates with time. Hence, the Optimal replacement policy is that the pipe should be replaced after the 4th break when the reliability = 0.0011.

Cite this paper
Amuji, H. , Ogbonna, C. , Ugwuanyim, G. , Iwu, H. and Nwanyibuife, O. (2018) Optimal Water Pipe Replacement Policy. Open Journal of Optimization, 7, 41-49. doi: 10.4236/ojop.2018.72002.
References
[1]   Watson, T., Colin, C., Mason, A. and Smith, M. (2001) Maintenance of Water Distribution System. The 36th Annual Conference of American Water Works Association.

[2]   Ozger, S. and Mays, L.W. (2004) Optimal Location of Isolation Valves in Water Distribution Systems: A Reliability/Optimization Approach. Urban Water Supply Management Tools. Mcgraw-Hill, New York.

[3]   Tobias, P.A. and Trindade, D.C. (1995) Applied Reliability. 2nd Edition, Chapman and Hall/CRC, New York.

[4]   Gregor, V.L., Kumar, M., Ranji, R., Downey, B., Jim, V., Ken, H., Sarat, S. and Emily, Z. (2006) An Adaptive Cyber Infrastructure for Threat Management in Urban Water Distribution Systems. Secure Water Resource Publications.

[5]   Park, S. and Loganathan, G.V. (2000) An Optimal Replacement Scheduling for Water Distribution Systems.

[6]   Jorgenson, D.W., McCall, J.J. and Rudner, R. (1967) Optimal Replacement Policy. Rand McNally and Co., Chicago.

[7]   Prem, K.G and Hira, D.S. (2008) Operation Research. S. Chand & Co. Ltd, Ram Nagar, New Delhi.

[8]   Lindquist, B.H. (2006). Maintenance of Repairable Systems. Journal of Statistical Planning and Inference, 136, 1701-1717.

[9]   Wu, S. (2017) The Model for Failure Process of a Repairable System. The 10th international Conference on Mathematical Methods in Reliability, Goenoble, 3-6 July 2017.

[10]   Giogio, M., Guida, M. and Pulcini, G. (2014) Repairable System Analysis in Presence of Covariates and Random Effects. Elsevier, 131, 271-281.

[11]   Masini, L., Antonio, P., Ruggeri, F. and Emmanuela, S. (2006) On Bayesian Models Incorporating Covariates in Reliability Analysis of Repairable Systems. Applied Statistics Models in Business and Industry, 20, 253-264.

[12]   Karol, A. (2015) Stochastic Modeling of the Repairable System. Journal of KONBiN, 3, 5-14.

[13]   Ruggeri, F. and Siva, S. (2005) On Modelling Change Points in Non-Homogeneous Poisson Process. Statistical Inference for Stochastic Processes, 8, 311-329.

[14]   Argiento, R., Pievatolo, A., Ruggeri, F., Cagno, E. and Macini, M. (2003) Seasonal Pattern and Double Measurement Scale in Modelling Failures in Underground Trains. 3rd International Conference on Mathematical Methods in Reliability.

[15]   Guo, R., Ascher, H. and Love, E. (2001) Towards Practical and Synthetical Modelling of Repairable Systems. Economic Quality Control, 16, 147-182.
https://doi.org/10.1515/EQC.2001.147

[16]   Graybill, F.A., Boes, D.C. and Mood, A.M. (1974) Introduction to the Theory of Statistics. McGraw-Hill Book Company, New York.

[17]   Chat Field, C. and Zidek, J.V. (1995) Modelling and Analysis of Stochastic Systems. Chapman and Hall, London.

[18]   Barlow, R.E. and Proschan, F. (1975) Statistical Theory of Reliability and Life Testing Probability Models. Holt, Rinehart and Winston, Inc., Boston.

[19]   Smith, R.L. (1991) Statistical Analysis of Reliability Data. Chapman and Hall, London.

[20]   Irwin, M. and John, E.F. (1987) Probability and Statistics for Engineers. 3rd Edition, Prentice-Hall of India Private Ltd., Delhi.

 
 
Top