Next Article in Journal
Using SWAT Model to Assess the Impacts of Land Use and Climate Changes on Flood in the Upper Weihe River, China
Previous Article in Journal
Hydrological Time Series Clustering: A Case Study of Telemetry Stations in Thailand
Previous Article in Special Issue
Integrating Ecosystem Services into Risk Assessments for Drinking Water Protection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimal Maintenance and Replacement Strategy for Deteriorating Water Mains

Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), School of Civil and Engineering, Sun Yat-sen University, Zhuhai 519082, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(13), 2097; https://doi.org/10.3390/w14132097
Submission received: 29 April 2022 / Revised: 17 June 2022 / Accepted: 27 June 2022 / Published: 30 June 2022
(This article belongs to the Special Issue Risk Assessment and Decision Support in Drinking Water Systems)

Abstract

:
Municipal water mains are built with a target service age of several decades. In such a long life, breaks can occur, even multiple times. Water mains can be maintained before or right at breaks. The former is referred to as the preventive strategy, whereas the latter is the corrective strategy. Depending on the costs of repair, replacement, and failure consequence, different strategies should typically be implemented in order to achieve the optimal watermain management in terms of life cycle costs. This study aims to investigate the optimal scenarios for the two strategies based on a two-time-scale (TTS) point process used to model the deterioration of water mains. The corrective strategy is to determine the optimal number n , where upon the n -th break, implementing a replacement for water main is justified, compared to a minimal repair. The preventive strategy is to determine the optimal replacement time in terms of pipe survival probability P s . Monte Carlo simulations are used to investigate the optimal n and P s considering a number of influential factors, including model parameters of the intensity function and ratios of maintenance, replacement, and consequence costs. Then, the full distributions of the life cycle costs are characterized with the mean of total life cycle costs being the target for optimization. Last, a case study is illustrated to demonstrate the application of both strategies in real water systems. An important finding is that with a typical pipe diameter of 400 mm and length of 200 m, the optimal n is typically less than five, and the optimal P s is below 50%.

1. Introduction

Water mains are key municipal infrastructure assets supporting the daily life of citizens [1]. Though water mains are designed to serve for several decades, they deteriorate with increasing service age, and can break before reaching the target design life. As a matter of fact, unexpected water main breaks are very common in reality. A main reason is that water mains are buried underground, and thus, it is very difficult to inspect and monitor their service conditions nondestructively, non-disruptively, timely, and easily [2,3].
Upon a water main break, municipal water managers have to decide whether to repair or replace the pipe. Even for water mains installed at some critical areas where a failure would cause severe loss, a preventive strategy can be adopted. A corrective strategy means maintaining a water main only after breaks occur, whereas a preventive strategy refers to maintaining a water main before it breaks.
Usually, if the consequence of failure is of little concern, maintenance can be implemented after a failure takes place. On the other end, if the failure would bring about massive loss [4,5], then a preventive strategy is warranted. Other than these two extreme cases, how to optimize the managing strategies depends on a set of influencing factors, for example, deterioration of the mains, repair cost, replacement cost, consequence, readiness of repairing technology, pipe location, ease of access, municipal budget, citizen affordability upon break, flooding risk, etc.
Although advanced managing strategies such as multiple target optimization on water pipe network level have been advocated [6,7], the application seems still limited in local areas due to the complexity. The current practice is still based on a ‘fire-fighting’ strategy coupled with a service age strategy, for example, replacing a pipe when the service age reaches 40 years regardless of its conditions. Moreover, the threshold service age at which a pipe will be replaced is selected somewhat arbitrarily, and with a lack of sound evidence. Obviously, improvements can be made.
This study investigates the optimization of managing strategies for water mains by setting the minimization of life cycle cost as the target. To do this, a proper deterioration model capturing breaks of water mains must be first determined. On the one hand, repairing a broken water main minimally will no doubt resume its functionality; however, the pipe condition remains the same, or at least not significantly improved, and degradation continues as before. A drawback for minimal repair, although it costs less, is that the water main could break again shortly after treatment. If so, another repair is then needed. Continuous repairs add up to the maintenance cost. On the other hand, replacing a broken water main with a new one starts a new life cycle. The new main is in excellent condition, and can be expected to serve for a much longer time before a break takes place. Nevertheless, a replacement typically costs much more than that of a minimal repair. Despite extensive studies on the management and maintenance of infrastructure assets, e.g., van Noortwijk [8], van Noortwijk and Frangopol [9], Jardine and Tsang [10], Bull [11], and Campbell et al. [12], under which scenarios a water main should be replaced rather than repaired remains a question to answer.
As a matter of fact, the estimation of the life cycle cost of water mains depends significantly on the deterioration. There have been extensive deterioration models proposed for water pipes, e.g., both break-history-independent [13,14,15,16,17,18,19,20] and -dependent [21,22,23]. Reviews of the modelling of water pipe deterioration can be seen in [2,24,25,26,27,28]. Obviously, history independence is a strong assumption in the prediction of water main breaks, which will inevitably reduce the prediction accuracy. As such, degradation models considering the breaking histories are superior and advantageous to be adopted in life cycle costing analysis (LCCA). Nevertheless, this is the part that has not yet received due attention in the current literature for the life cycle management of municipal water pipe networks. Existing studies, though they have made significant advances in water pipe management, have focused largely upon network effects and a simplification of methods for LCCA, e.g., Termes-Rifé et al. [29], Roshani and Filion [30], Naoum-Sawaya et al. [31], Lee et al. [32], Ghobadi et al. [33], among other excellent works. A brief review of LCCA for water pipes can be found in Ghobadi et al. [33].
In this paper, we examine two simple maintenance and replacement strategies for water mains. The first is a corrective strategy which is to determine the optimal number of breaks n at which the water main is replaced; before reaching the n -th break, repair is implemented upon a break. The second is a preventive strategy, which is to determine the optimal survival probability at which the water main is replaced. The water main break, along with time, is modelled as a two-time-scale (TTS) point process, proposed by Lin and Yuan [34]. The influences of degradation intensity, ratios of maintenance, replacement, and failure consequence costs are discussed. The Monte Carlo simulation approach is adopted to conduct the analyses. A case study is provided in the end to elaborate the application of the two strategies for the real world. The outcomes of this study are valuable to municipalities for implementing better management strategies for water mains.
The rest of the paper is organized as follows. The next section introduces the TTS degradation model for water mains and the simulation technique for breaks. The definition of life cycle cost is formulated in Section 3 with respect to corrective and preventive strategies. Section 4 shows the results of strategy optimization, followed by a case study in Section 5. The paper is concluded in Section 6, followed by references.

2. Degradation Model

2.1. Modelling Water Main Breaks as a Point Process

Upon a break, a water main is not always replaced by a new one; instead, it is usually repaired, which brings the main back to functioning. Since the working conditions, e.g., surrounding soils, inner pressure, corrosion, external loading, etc., are typically unknown or just partially known, physical modeling of the birth and growth of defects on a pipe turns out to be very challenging, if not impossible [35]. An alternative is to model the occurrence of a water main break phenomenically [36]. If ordering all water main breaks in terms of breaking times, then a series of ‘points’ in the timeline are obtained, as shown in Figure 1. Due to the unknown, the occurrence of these points is uncertain. This can then be simply modelled as a stochastic point process.
For a simple point process, it can be fully characterized by a conditional intensity function describing the potential of an event to take place within the next interval of time given the event history. According to Cox and Isham [37] and Cook and Lawless [38], the conditional intensity λ ( t | H ( t ) ) can be mathematically expressed as:
λ ( t | H ( t ) ) = lim Δ t 0 1 Δ t Pr { N [ t , t + Δ t ) = 1 | H ( t ) }
where H ( t ) = history of the N ( t ) until before time t ; and N [ t , t + Δ t ) is the number of occurrences within the time interval [ t ,   t + Δ t ) . where Δ t > 0 . When Δ t approaches zero, the probability for more than one event to occur within the Δ t period is ignorable. This assumption can be written as Pr { N [ t , t + Δ t ) 2 | H ( t ) } = o ( Δ t ) . Discussion on the general theory of point processes in detail can be found in Cox and Isham [37], Jacobsen [39], Cook and Lawless [38], and Daley and Vere-Jones [40], among others.

2.2. The Two-Time-Scale (TTS) Intensity Function

The break of water mains has been most commonly modelled as a Poisson process. An implicit assumption of this type of modelling is that breaks to occur do not depend on the breaking history. Nevertheless, this is a strong assumption and does not have any solid scientific support. To account for the influence of breaking history upon the future breaks, Lin and Yuan [34] developed a two-time-scale (TTS) intensity function for breaks of water mains. The intensity function includes an overall time scale representing the entire pipe deterioration, as well as a local time scale representing the renewal effect. The TTS function is expressed as:
λ = α t β w γ
where t = current time; w = sojourn time since the last break; α , β , and γ are model parameters. Obviously, the break history in Equation (2) is simplified as inter-occurrence duration w .
Lin and Yuan [34] demonstrated the capability of the TTS model in capturing different deterioration trends, i.e., both monotonic and non-monotonic. In particular, it is able to model a bathtub curve deteriorating trend, which is of great interest to engineers, as shown in Figure 2, where the surface of λ is plotted against model parameter β and the sojourn time w given α = 0.1 , β = 0.6 , γ = 0.3 , and T L = 20 . Note that these values for α , β , and γ are intentionally selected here for λ to be in a bathtub curve shape, which is a highly desirable deterioration trend for engineers to model. With an increase in the sojourn time w from zero, the intensity λ first drops dramatically until reaching a minimal, then increases back as w continues to increase. The parameter α controls the overall scale of λ , whereas β and γ are the shape factors controlling the trend together. With the increasing of w , λ first decreases to a minimum and then increases, showing a bathtub curve. Table 1 summarizes the trends that correspond to the permutations of β and γ . In general, the larger the λ , the higher the breaking potential of a water main. Readers are directed to Lin and Yuan [34] for more explorations on Equation (2).

2.3. Simulating a TTS Point Process

Given λ , the probability of a water main surviving through a time duration [ t 1 ,   t 2 ] can be computed as [38,41]:
P s = Pr { N [ t 1 , t 2 ] = 0 | H ( t 1 ) } = exp [ t 1 t 2 λ ( u | H ( u ) ) d u ]
Equation (3) suggests that once P s , t 1 , and λ are known, t 2 can be determined by solving the equation. Based on the breaking data in the city of Toronto, Canada, Lin and Yuan [34] showed that the TTS is superior to Poisson process, and water mains follow a bathtub deterioration trend. Therefore, the TTS intensity function is adopted here, and the model parameters are taken as α = 0.1 , β = 0.6 , and γ = 0.3 for analyses. Figure 3 shows 100 simulated TTS point processes starting from t 1 = 0 . For the simulations, P s is randomly sampled from a uniform distribution of U ~ [ 0 , 1 ] .
On average, the first two breaks occur within t = 10 , i.e., t 1 = 6.5 , t 2 = 9.8 . A total of eight breaks can be expected within t = 20 . It should be noted that here, t does not have a unit, meaning that it can be a year, month, minute, second, or others, depending on the data used for investigation. Another observation is that the simulated break events tend to occur either one shortly after another, or none within a relatively long period of time. This matches observations from engineering experience.

3. Definition of Life Cycle Cost

Upon a break, a water main pipe can be either maintained or replaced. Maintenance resumes the serviceability of the main while the deterioration continues at the original rate [42,43,44]. On the other hand, a replacement brings in a new main pipe and starts a new life cycle. Suppose the cost of a maintenance is C m and the cost of a replacement is C r ; then, it is reasonable to assume C r > C m , especially when the no-dig technology is employed [45,46,47,48,49]. In addition to C m and C r , failure consequence is another primary source of costs. For example, a water main break could block the traffic for hours, flooding infrastructure assets, and thus, resulting in damage and loss. The cost due to failure consequence is denoted as C c in this study. Depending on λ , C m , C r , and C c , different strategies must be practiced to reach the goal of minimum life cycle cost for water main asset management. Two asset management approaches are considered here: corrective and preventive strategies. The following introduces these two strategies, including the definition of life cycle cost and the steps to compute these costs.

3.1. Corrective Strategy

A corrective strategy is that the maintenance, repair, or replacement is implemented only after water main breaks take place. The cost for each break event is either ( C m + C c ) for maintenance or ( C r + C c ) for replacement. As was pointed out earlier in this paper, a repair or maintenance will only resume the serviceability of the main, but does not improve its condition. The degradation is still accelerating with service age, meaning that the water main will break more frequently as time goes. After some point, minimal repair is no longer cost-effective compared to a renewal.
Consider a simple corrective maintenance strategy: replacing the water main with a new one at the n -th breaks; then, the question becomes: what is the optimal n ? This idea is illustrated in Figure 1. In terms of the concept of net present value, the total cost of the water main due to maintenance and replacement in a sufficient long planning horizon can be written as:
C L | p , n = k = 1 k = p { j = 1 j = n 1 D ( t k , j ) [ C c ( t k , j ) + C m ( t k , j ) ] + D ( t k , n ) [ C c ( t k , n ) + C r ( t k , n ) ] }
where p = number of lifecycles; D = discount factor converting future cost to net present value. Usually, D = ( 1 + d ) t , with d being the discounted rate. The most commonly used discounted rate is 5% annually [50,51]. This study adopts d = 5 % as a reference for analyses.
Assuming C m , C r , and C c are all constant with time, then for convenience, the life cycle cost expressed by Equation (4) can be normalized by C m as:
C L | p , n C m = k = 1 k = p [ j = 1 j = n 1 D ( t k , j ) ( C c C m + 1 ) + D ( t k , n ) ( C c C m + C r C m ) ]  
Equation (5) indicates that the optimal n mainly depends on the cost ratios of C c C m and C r C m . If C m = C r , then obviously the main should be replaced at each break. The breaking event becomes a renewal process. This also applies to cases with severe consequence, i.e., C r C c . If C m C r , then the optimal solution for n will depend on the ratio of C c / C r .
The steps to determine the optimal n value are as follows.
Step 1. Determine the values for the ratios of C c C m and C r C m , discounting factor d , and for model parameters α , β , and γ ;
Step 2. Select an initial n value, e.g., n = 1 ;
Step 3. Generate a realization of the TTS point process given the pre-determined α , β , and γ ;
Step 4. Compute the normalized life cycle cost C L | p , n C m using Equation (5);
Step 5. Repeat Steps 3 to 4 for a sufficient number of times, i.e., 100,000 times in this paper, to obtain the full distribution of life cycle cost C L | p , n C m ;
Step 6. Calculate the statistics of C L | p , n C m , i.e., mean, coefficient of variation;
Step 7. Repeat Steps 2 to 6 to compute the means and COVs of C L | p , n C m for different n values; then, the one that gives minimal mean of C L | p , n C m is determined as the optimal;
Step 8. Change the values of C c C m , C r C m , d , α , β , and γ to investigate their influences on the optimal n .

3.2. Preventive Strategy

As will be demonstrated later in this study, for cases with high failure consequence, a corrective strategy might not be adequate, and the life cycle cost can still be further optimized by adopting a preventive strategy, that is, replacing a water main before it fails.
Since the exact time instance for a break to take place is unknown, this study takes the survival probability P s computed as Equation (3) as a time indicator for taking a replacement action. The idea is illustrated in Figure 4.
Suppose the last break of the water main is at t 1 ; then, the steps to compute the cost of the next life cycle C L are as follows.
Step 1. Select an initial survival probability P s ;
Step 2. Compute the corresponding t P s using Equation (3);
Step 3. Generate a random sample u from a uniform distribution U ~ [ 0 , 1 ] ;
Step 4. Determine t 2 corresponding to the u with Equation (3);
Step 5. Compare t 2 against t P s ; if t 2 > t P s , then the cost is calculated as C L = D ( t P s ) × C r , and t 1 is updated as t 1 = t P s ; otherwise, t 2 t P s , C L = D ( t 2 ) × ( C r + C c ) , and t 1 = t 2 ;
Step 6. Repeat steps 2 to 5 for a large number of times; then sum all the costs up to be the total cost of one simulated path;
Step 7. Repeat Step 6 a large number of times, and take the average as the mean total life cycle cost.
Step 8. Change the initial P s and repeat steps 2 to 7 until the minimum mean total life cycle cost is reached; then, the corresponding P s is determined as the optimal P s .

4. Results of Strategy Optimization

A reference example is firstly set up for analysis purposes: α = 0.1 , β = 0.6 , γ = 0.3 , d = 5 % , C c / C m = 1 , and C r / C m = 10 . This setting assumes that the loss due to failure is the same as the cost for a maintenance, whereas the cost of a replacement is 10 times to that of a maintenance. These ratios will be varied to reflect different scenarios, as well as their influences on establishing optimal strategies.

4.1. Corrective Strategy

Figure 5 shows the distributions and the means of the normalized total life cycle costs C L / C m (in net present value) against the n value for the reference example. In the figure, the mean C L / C m was averaged from 10,000 simulations; each simulation had 1000 cycles. The means, medians, and COVs for the C L / C m in Figure 5 are summarized in Table 2. As shown, all the distributions are more or less symmetric. Applying K-S tests shows that most of them can be considered normal distributions. The mean C L / C m is about 30 at n = 1 , but decreases to about 17.5 as n increases to 7; after that, C L / C m increases to about 22 for n = 50 . These results suggest that the optimal n value is 7, meaning a water main pipe should be replaced by a new one rather than repaired every 7 breaks. The difference of the mean C L / C m reaches 12.5, i.e., n = 7 against n = 1 .
Another observation from Figure 5 is that the coefficient of variation (COV) of C L / C m exhibits a decreasing trend with increasing n value. The C L / C m COV is as large as 32% for the n = 1 case, and continuously drops to about 22% for the n = 50 case. It is worth mentioning here that this study takes the mean C L / C m as the target for life cycle cost optimization. For those taking both mean and COV of C L / C m as the targets, a different optimal n value can be resulted.
It is also worth noting that the mean C L / C m has a very flat bottom for n , varying from 5 to 10. The differences in the mean C L / C m values are practically negligible. Hence, though the optimal value is n = 7 , it also seems acceptable to take any integers between 5 and 10 for n .
Keeping the other two parameters unchanged while varying the third one for the three TTS model parameters, the effects of α , β , and γ on the final outcome of the optimal n value are examined. The results are provided in Figure 6. The increases of α and β lead to a decrease of optimal n value, whereas for γ , the opposite is observed. This can be explained by the fact that both α and β are parameters acting on the overall level of λ . The increases of them will result in a higher deterioration intensity, meaning that the water main tends to break more and more frequently. This causes a quick add-up of the total repair cost, and thus, the need of a replacement is justified so as to yield an economic benefit. After replacement, t is re-zeroed, and the λ dramatically drops. Hence, a lower optimal n matches better the target of minimizing the total life cycle cost. On the other hand, parameter γ is a factor acting on the local level of λ . A smaller γ (negative value) indicates a stronger local renewal effect of the water main degradation, which, in turn, calls for a more frequent replacement of the asset.
The second observation is that the effect of β is much larger than those by α and γ . As β drops from 0.8 to 0.2, the optimal n correspondingly increases from 6 to 63, which is over an order of magnitude of difference. For α and γ , the variation ranges are within 6 to 10 for α decreasing from 0.20 to 0.05, and 4 to 7 for γ increasing from −7 to −0.1. This reveals that β is the most influential parameter in λ .
Despite this, it should be emphasized that, in reality, the sojourn time of a break of water mains could be very long, and the cost, after converted to the net present value, would become very small. An indication of this is that the cost ratios C c / C m and C r / C m could play a more important role than the three model parameters.
Figure 7 shows the influences of ratios of consequence to maintenance cost C c / C m and replacement to maintenance cost C r / C m on the determined optimal n value. For C r / C m = 1 , the optimal n value will always be 1, regardless of C c / C m . This is understandable, as a replacement gives a new water main whose condition is better than a repaired one with a certain service age without extra cost. For C r / C m > 1 , a smaller optimal n is obtained for a larger C c / C m . For example, for C r / C m = 10 , the optimal n is found to be 14 if C c / C m = 0 , but drops to 1 after C c / C m exceeds 9. For C r / C m , its increase will lead to an increase of the optimal n value. Given C c / C m = 10 , the optimal n value is 1 for C r / C m = 10 , and increases to 4 for C r / C m = 20 . Figure 7 demonstrates the importance of C c / C m and C r / C m in the decision-making of maintenance and replacement strategies for water mains. This also implies the need to implement pertinent strategies for water mains located in different areas, e.g., rural, commercial, residential, or industrial areas.
The last observation from Figure 7 is that it is not uncommon for the optimal n value to be 1 under the corrective strategy framework. This suggests implementation of a preventive strategy could be a better solution for those scenarios.

4.2. Preventive Strategy

For preventive strategy, there is no maintenance; only replacement is considered. As a result, C m is removed from the calculation of life cycle cost C L ; only C c and C r are taken into account. Figure 8 shows the normalized life cycle cost C L / C r varying with replacement at different survival probabilities of water mains. The means, medians, and COVs for the C L / C r in Figure 8 are summarized in Table 3. For the scenario of C c / C r = 10 , the optimal survival probability for implementing a replacement is found to be P s = 70%. The corresponding life cycle cost is C L / C r = 28.2 . If P s is increased to 90%, C L / C r would increase to be about 33.0; if P s is decreased to 0%, which is equivalent to the corrective strategy of n = 1 , C L / C r would increase to be about 32.0. This means that turning the strategy from corrective to preventive, in this case, would reduce the life cycle cost by (32 − 28.2)/32 = 12%. For other cases, the improvement could be larger.
Similarly, though the mean C L / C r reaches its minimum at P s = 70 % , the COV does not concurrently reach its minimum. The COV of C L / C r is found to be 0.38 at P s = 70 % , but only 0.29 at P s = 90 % and 0.33 P s = 0 % . Again, if minimizing both means, and COV of C L / C r is the target, then the optimal P s could not be 70%. This observation also calls for the development of a new approach for the optimization of life cycle cost, e.g., multiple targets or information entropy. The means of C L / C r for P s between 60% and 70% are essentially the same; however, the sojourn times corresponding to these two P s values are different. For cases where assets will be decommissioned after a limited planning horizon, the advantage of C L / C r by setting P s = 70 % over P s = 60 % will then be more evident.
Figure 9 shows the influences of C c / C r on the optimal survival probability for replacement P s with consideration of discounted rates at d = 1%, 5%, and 10%. As C c / C r increases, i.e., the consequence due to water main breaks becomes higher, the main should be replaced earlier, i.e., at larger P s . For larger discounting rate with a fixed C c / C r , a higher optimal P s for replacement is reached. This is because a larger d results in less costs than the net present value, which, in turn, gives higher optimal P s . Interestingly, it is found that the effect of d diminishes at a larger C c / C r value, for example, for C c / C r exceeding 20, the optimal P s appears to be independent on d . This can be explained by the fact that for large C c / C r , the optimal P s is also high, i.e., over 80%. The water main will be replaced ‘shortly’ after being put into service. Due to this short time duration, the effect of d is small.
Table 4 summarizes the influences of model parameters, as well as cost ratios on the optimal corrective and preventive strategies. In general, the model parameter β and the C c / C r ratio are the two most influential factors on the decision-making of maintenance and replacement of water mains.

5. Case Study

This section presents the analysis of optimal maintenance and replacement strategies for water mains using model parameters obtained from a real project in a relatively young municipality in Canada. The water mains in the municipality are mainly ductile iron pipes and concrete pressure pipes. This case study focuses on ductile pipes. Most of the pipes were put into service before 1998; the oldest one was installed in 1967. The pipes are typically 400 mm in diameter, and less than 200 m in length. Other information relating to deterioration include pressure zone, surface type, and soil types; however, these pieces of information are not complete. The break records for these pipes are available only after January 1st, 1998. The records are basically the time instances when the breaks took place.
Based on the water main inventory and break records, Lin and Yuan [34] reported that the degradation of ductile iron water mains in the municipality can be modelled by a TTS intensity function including the length and diameter, expressed as:
λ = α ( l l o ) k L ( D D o ) k D t β w γ
where l and D are length and diameter of a water main, respectively; l o = 50 m and D o = 400 mm are constants used to normalize l and D , respectively. The model parameters were reported to be α = 0.0071 , κ L = 0.4344 , κ D = 1.0703 , β = 0.4619 , and γ = 0.5380 .
Based on the lengths and diameters reported for the municipality by Lin and Yuan [34], we consider an example water main with l = 200 m and D = 400 mm. Figure 10 shows the distributions of predictions for the first five breaks for this example water main. The means, medians, and COVs of these predictions are summarized in Table 5. The predictions spread widely, with a COV of 1.055 for the first one, but reducing to 0.546 for the fifth one. The means are 103.0, 148.5, 185.3, 217.0, and 241.4 years. The medians are 64.4, 118.0, 154.9, 194.6, and 218.1 years. Obviously, the distributions of the predictions are not symmetric, since the means and medians are significantly deviated from each other; this deviation reduces for later breaks.
When corrective strategy is implemented for this water main, it is found that for most cases, the optimal n value is under 5. A typical n value of 2 or 3 is obtained. This can be expected, as the breaking times are so large (over 100 years) that the gaps of costs due to maintenance or replacement are small in terms of net present value. Similar to earlier analyses, an increase of C c / C m or a decrease of C r / C m would lead to a smaller optimal n value, as shown in Figure 11.
When a preventive strategy is preferred, the optimal survival probability to carry out a replacement would be 10% for C c / C r = 1 (See Figure 12). If the consequence of failure is much more severe than a replacement, for example, C c / C r > 20 , then the water main is best replaced at a survival probability as early as 40% (or equivalently, the probability of failure reaches 60%). In general, it appears that the best time to implement a replacement of a water main is when its survival probability lies below 50%.
It is worth mentioning that the model parameters and ratios are firmly based on the case study used here. Though the analysis process can be applied to other cases, caution should be practiced to see whether the values of the parameters and ratios are applicable to those cases or not. It is also worth emphasizing here that the hydraulic pressure at which a water main operates is an important factor for pipe ageing, and thus, for life cycle costing analysis. This factor is not explicitly accounted for in this case study due to the lack of data. Nevertheless, it should be formulated into the intensity function whenever possible to enable a more accurate LCC analysis.

6. Conclusions

This paper presented analyses of both corrective and preventive repair and replacement strategies for water mains. The recurrent water main breaks with increasing service age were modelled as a two-time-scale (TTS) stochastic point process. The corrective strategy was simplified as finding the optimal n -th break at which implementing a replacement is more cost-effective than a minimal repair. The preventive strategy is set as determining the optimal survival probability P s for a water main to receive a replacement. Both strategies were investigated from the perspective of life cycle cost (LCC). The Monte Carlo simulation technique was adopted to conduct the analyses. A case study was illustrated to demonstrate the application of the strategies for real water pipes. The main conclusions drawn from this study are the following.
(1)
Both of the optimal corrective and preventive strategies, optimal n and P s , depend highly on the degradation intensity of water mains, as well as the ratios of maintenance, replacement, and failure consequence costs. For corrective strategy, the higher the breaking intensity, the higher the failure consequence loss, and the lower the replacement cost, the lower the optimal n value. For preventive strategy, the higher the ratio of consequence loss to replacement cost, the higher the optimal P s . In addition, the effect of discounting rate diminishes when the ratio of consequence loss to replacement cost becomes large.
(2)
The optimal maintenance and replacement strategies are sensitive to the overall shape parameter of the service age in the TTS intensity function, whereas the effects of the scale parameter and the local shape parameter of the sojourn time since the last break on the optimal n and P s are comparatively secondary.
(3)
A case study using a typical pipe diameter of 400 mm and length of 200 m shows that the optimal n value is typically under 5, and the optimal P s lies below 50%.
Though the approach adopted in this study for carrying out a life cycle cost (LCC) analysis is easy to implement, it has three limitations. First, the TTS deterioration model is used, which is advantageous, as it takes breaking history into account; however, the parameter estimation for the TTS model would be very computationally demanding if the early break history is missing. Moreover, only the mean of LCC is set as the single target for optimization in this study. If the COV is also considered, then the optimization outcomes could be different. Last, the LCC analysis presented in this study is on individual pipe level, rather than pipe network level. Note that the latter two issues exist not only for this study, but also for many other studies in the literature. It is suggested that the extension from pipe to network level LCC analysis based on the TTS model should be carried out considering multiple objective optimizations. Related topics can be referenced to, e.g., [6,52].

Author Contributions

Conceptualization, P.L., S.H. and B.M.; methodology, P.L.; validation, X.C. and S.H.; formal analysis, X.C.; investigation, P.L.; writing—original draft preparation, P.L.; writing—review and editing, S.H.; supervision, B.M.; funding acquisition, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities, Sun Yat-sen University [22hytd06], the National Natural Science Foundation of China [52008408], the Guangdong Basic and Applied Basic Research Foundation [2021A1515012088], Science and Technology Program of Guangzhou, China [202102021017]. And The APC was funded by the Central Universities, Sun Yat-sen University [22hytd06].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sköld, N.-P.; Bergion, V.; Lindhe, A.; Keucken, A.; Rosén, L. Risk-Based Evaluation of Improvements in Drinking Water Treatment Using Cost-Benefit Analysis. Water 2022, 14, 782. [Google Scholar] [CrossRef]
  2. Liu, Z.; Kleiner, Y.; Rajani, B.; Wang, L.; Condit, W. Condition Assessment Technologies for Water Transmission and Distribution Systems; United States Environmental Protection Agency (EPA): Washington, DC, USA, 2012; EPA/600/R-12/017; p. 108. [Google Scholar]
  3. Bergion, V.; Lindhe, A.; Sokolova, E.; Rosén, L. Accounting for Unexpected Risk Events in Drinking Water Systems. Expo. Health 2021, 13, 15–31. [Google Scholar] [CrossRef]
  4. Sjöstrand, K.; Klingberg, J.; Zadeh, N.; Haraldsson, M.; Rosén, L.; Lindhe, A. The Value of Water—Estimating Water-Disruption Impacts on Businesses. Water 2021, 13, 1565. [Google Scholar] [CrossRef]
  5. Sjöstrand, K.; Lindhe, A.; Söderqvist, T.; Rosén, L. Water Supply Delivery Failures—A Scenario-Based Approach to Assess Economic Losses and Risk Reduction Options. Water 2020, 12, 1746. [Google Scholar] [CrossRef]
  6. Scholten, L.; Scheidegger, A.; Reichert, P.; Mauer, M.; Lienert, J. Strategic rehabilitation planning of piped water networks using multi-criteria decision analysis. Water Res. 2014, 49, 124–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Sjöstrand, K.; Lindhe, A.; Söderqvist, T.; Rosén, L. Sustainability assessments of regional water supply interventions—Com-bining cost-benefit and multi-criteria decision analyses. J. Environ. Manag. 2018, 225, 313–324. [Google Scholar] [CrossRef]
  8. van Noortwijk, J.M. Explicit formulas for the variance of discounted life-cycle cost. Reliab. Eng. Syst. Saf. 2003, 80, 185–195. [Google Scholar] [CrossRef]
  9. van Noortwijk, J.M.; Frangopol, D.M. Two probabilistic life-cycle maintenance models for deteriorating civil infrastructures. Probabilistic Eng. Mech. 2003, 19, 345–359. [Google Scholar] [CrossRef]
  10. Jardine, A.K.; Tsang, A.H. Maintenance, Replacement, and Reliability: Theory and Applications; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  11. Bull, J.W. Life Cycle Costing for the Analysis, Management and Maintenance of Civil Engineering Infrastructure; Whittles Publishing: Dunbeath, UK, 2015. [Google Scholar]
  12. Campbell, J.D.; Jardine, A.K.; McGlynn, J. Asset Management Excellence: Optimizing Equipment Life-Cycle Decisions; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  13. Constantine, A.; Darroch, J. Pipeline Reliability. Stochastic Models in Engineering Technology & Management; Osaki, S., Murthy, D.N.P., Eds.; World Scientific Publishing: Singapore, 1993. [Google Scholar]
  14. Constantine, G.; Darroch, J.; Miller, R. Predicting underground pipeline failure. Water-Melb. Artarmon. 1996, 23, 9–10. [Google Scholar]
  15. Economou, T.; Kapelan, Z.; Bailey, T.C. On the prediction of underground water pipe failures: Zero inflation and pipe-specific effects. J. Hydroinform. 2012, 14, 872–883. [Google Scholar] [CrossRef] [Green Version]
  16. Kleiner, Y.; Rajani, B. Using limited data to assess future needs. J. Am. Water Work. Assoc. 1999, 91, 47–61. [Google Scholar] [CrossRef] [Green Version]
  17. Kleiner, Y.; Rajani, B. I-WARP: Individual Water mAin Renewal Planner. Drink. Water Eng. Sci. 2010, 3, 71–77. [Google Scholar] [CrossRef] [Green Version]
  18. Le Gat, Y. Extending the Yule process to model recurrent pipe failures in water supply networks. Urban Water J. 2014, 11, 617–630. [Google Scholar] [CrossRef]
  19. Røstum, J. Statistical Modelling of Pipe Failures in Water Networks. Ph.D. Thesis, Norwegian University of Science and Technology, Trondheim, Norway, February 2000. [Google Scholar]
  20. Watson, T.G.; Christian, C.D.; Mason, A.J.; Smith, M.H.; Meyer, R. Bayesian-based pipe failure model. J. Hydroinform. 2004, 6, 259–264. [Google Scholar] [CrossRef] [Green Version]
  21. Mailhot, A.; Pelletier, G.; Noël, J.-F.; Villeneuve, J.-P. Modeling the evolution of the structural state of water pipe networks with brief recorded pipe break histories: Methodology and application. Water Resour. Res. 2000, 36, 3053–3062. [Google Scholar] [CrossRef]
  22. Pelletier, G.; Mailhot, A.; Villeneuve, J.-P. Modeling water pipe breaks—Three case studies. J. Water Resour. Plan. Manag. 2003, 129, 115–123. [Google Scholar] [CrossRef]
  23. Scheidegger, A.; Scholten, L.; Maurer, M.; Reichert, P. Extension of pipe failure models to consider the absence of data from replaced pipes. Water Res. 2013, 47, 3696–3705. [Google Scholar] [CrossRef] [Green Version]
  24. Clair, A.M.S.; Sinha, S. State-of-the-technology review on water pipe condition, deterioration and failure rate prediction models. Urban Water J. 2012, 9, 85–112. [Google Scholar] [CrossRef]
  25. Kleiner, Y.; Rajani, B. Comprehensive review of structural deterioration of water mains: Statistical models. Urban Water 2001, 3, 131–150. [Google Scholar] [CrossRef] [Green Version]
  26. Nishiyama, M.; Filion, Y. Review of statistical water main break prediction models. Can. J. Civ. Eng. 2013, 40, 972–979. [Google Scholar] [CrossRef]
  27. Rajani, B.; Kleiner, Y. Comprehensive review of structural deterioration of water mains: Physically based models. Urban Water 2001, 3, 151–164. [Google Scholar] [CrossRef] [Green Version]
  28. Scheidegger, A.; Leitão, J.P.; Scholten, L. Statistical failure models for water distribution pipes—A review from a unified perspective. Water Res. 2015, 83, 237–247. [Google Scholar] [CrossRef] [PubMed]
  29. Termes-Rifé, M.; Molinos-Senante, M.; Hernández-Sancho, F.; Sala-Garrido, R. Life Cycle Costing: A tool to manage the urban water cycle. J. Water Supply Res. Technol. 2013, 62, 468–476. [Google Scholar] [CrossRef]
  30. Roshani, E.; Filion, Y.R. Event-Based Approach to Optimize the Timing of Water Main Rehabilitation with Asset Management Strategies. J. Water Resour. Plan. Manag. 2014, 140, 04014004. [Google Scholar] [CrossRef] [Green Version]
  31. Naoum-Sawaya, J.; Ghaddar, B.; Arandia, E.; Eck, B. Simulation-optimization approaches for water pump scheduling and pipe replacement problems. Eur. J. Oper. Res. 2015, 246, 293–306. [Google Scholar] [CrossRef]
  32. Lee, H.; Shin, H.; Rasheed, U.; Kong, M. Establishment of an Inventory for the Life Cycle Cost (LCC) Analysis of a Water Supply System. Water 2017, 9, 592. [Google Scholar] [CrossRef]
  33. Ghobadi, F.; Jeong, G.; Kang, D. Water Pipe Replacement Scheduling Based on Life Cycle Cost Assessment and Optimization Algorithm. Water 2021, 13, 605. [Google Scholar] [CrossRef]
  34. Lin, P.; Yuan, X.-X. A two-time-scale point process model of water main breaks for infrastructure asset management. Water Res. 2018, 150, 296–309. [Google Scholar] [CrossRef]
  35. Lin, P.Y.; Yuan, X.X.; Tovilla, E. Integrative modeling of performance deterioration and maintenance effectiveness for infrastructure assets with missing condition data. Comput. Aided Civ. Infrastruct. Eng. 2019, 34, 677–695. [Google Scholar] [CrossRef]
  36. Yuan, X.X. Principles and guidelines of deterioration modelling for water and waste water assets. Infrastruct. Asset. Manag. 2017, 4, 19–35. [Google Scholar] [CrossRef]
  37. Cox, D.R.; Isham, V. Point Processes; Taylor & Francis: Abingdon, UK, 1980. [Google Scholar]
  38. Cook, R.J.; Lawless, J.F. The Statistical Analysis of Recurrent Events; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar] [CrossRef]
  39. Jacobsen, M. Point Process Theory and Applications: Marked Point and Piecewise Deterministic Processes; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  40. Daley, D.J.; Vere-Jones, D. An Introduction to the Theory of Point Processes: Volume II: General Theory and Structure; Springer: New York, NY, USA, 2007. [Google Scholar]
  41. Lin, P.; Yuan, X.-X. An Advanced Statistical Method for Point Process Modelling with Missing Event Histories. IEEE Access 2020, 8, 210081–210098. [Google Scholar] [CrossRef]
  42. Pietrucha-Urbanik, K. Failure Prediction in Water Supply System—Current Issues; Springer: Berlin/Heidelberg, Germany, 2015; pp. 351–358. [Google Scholar] [CrossRef]
  43. Pietrucha-Urbanik, K. Assessing the Costs of Losses Incurred as a Result of Failure; Springer: Berlin/Heidelberg, Germany, 2016; pp. 355–362. [Google Scholar] [CrossRef]
  44. Pietrucha-Urbanik, K.; Tchórzewska-Cieślak, B. Cost Analysis of Water Pipe Failure; Springer: Berlin/Heidelberg, Germany, 2020; pp. 411–424. [Google Scholar]
  45. Hashemi, B.; Najafi, M.; Mohamed, R. Cost of Underground Infrastructure Renewal: A Comparison of Open-Cut and Trenchless Methods. In Proceedings of the 12th International Pipelines Conference 2008, Atlanta, GA, USA, 22–27 July 2008. [Google Scholar]
  46. Najafi, M.; Gokhale, S.B. Trenchless Technology: Pipeline and Utility Design, Construction, and Renewal; McGraw-Hill Education: New York, NY, USA, 2005. [Google Scholar]
  47. Wu, Y.; Kang, C.; Nojumi, M.M.; Bayat, A.; Bontus, G. Current water main rehabilitation practice using trenchless technology. Water Pract. Technol. 2021, 16, 707–723. [Google Scholar] [CrossRef]
  48. Yang, M.-D.; Su, T.-C.; Wu, T.-Y.; Huang, K.-S. No-Dig Inspection Technologies for Underground Pipelines. J. GeoEng. 2010, 5, 99–104. [Google Scholar]
  49. Zhao, Y.; Ma, B.; Ariaratnam, S.T.; Zeng, C.; Yan, X.; Wang, F.; Wang, T.; Zhu, Z.; He, C.; Shi, G.; et al. Structural performance of damaged rigid pipe rehabilitated by centrifugal spray on mortar liner. Tunn. Undergr. Space Technol. 2021, 116, 104117. [Google Scholar] [CrossRef]
  50. Lemer, A. Life Cycle Costing: For the Analysis, Management and Maintenance of Civil Engineering Infrastructure. Constr. Manag. Econ. 2015, 33, 689–691. [Google Scholar] [CrossRef]
  51. van den Boomen, M.; Schoenmaker, R.; Wolfert, A.R.M. A life cycle costing approach for discounting in age and interval replacement optimisation models for civil infrastructure assets. Struct. Infrastruct. E 2018, 14, 1–13. [Google Scholar] [CrossRef]
  52. Mittal, A.; Scholten, L.; Kapelan, Z. A review of serious games for urban water management decisions: Current gaps and future research directions. Water Res. 2022, 215, 118217. [Google Scholar] [CrossRef]
Figure 1. Concept of modelling failure events as a point process.
Figure 1. Concept of modelling failure events as a point process.
Water 14 02097 g001
Figure 2. Three-dimensional plot of λ against model parameter β and the sojourn time w given α = 0.1 , β = 0.6 , γ = 0.3 , and T L = 20 .
Figure 2. Three-dimensional plot of λ against model parameter β and the sojourn time w given α = 0.1 , β = 0.6 , γ = 0.3 , and T L = 20 .
Water 14 02097 g002
Figure 3. Simulated water main breaking events using the TTS point process model assuming α = 0.1 , β = 0.6 , and γ = 0.3 .
Figure 3. Simulated water main breaking events using the TTS point process model assuming α = 0.1 , β = 0.6 , and γ = 0.3 .
Water 14 02097 g003
Figure 4. Illustration of the preventive strategy in terms of survival probability.
Figure 4. Illustration of the preventive strategy in terms of survival probability.
Water 14 02097 g004
Figure 5. Normalized life cycle cost C L / C m varying with replacement at the n -th failure n ( d = 5 % , C c / C m = 1 , and C r / C m = 10 ).
Figure 5. Normalized life cycle cost C L / C m varying with replacement at the n -th failure n ( d = 5 % , C c / C m = 1 , and C r / C m = 10 ).
Water 14 02097 g005
Figure 6. Effects of the model parameters α , β , and γ in the TTS intensity function on the optimal n value ( d = 5 % , C c / C m = 1 , and C r / C m = 10 ).
Figure 6. Effects of the model parameters α , β , and γ in the TTS intensity function on the optimal n value ( d = 5 % , C c / C m = 1 , and C r / C m = 10 ).
Water 14 02097 g006
Figure 7. Influences of ratios of consequence to maintenance cost C c / C m and replacement to maintenance cost C r / C m on optimal n value.
Figure 7. Influences of ratios of consequence to maintenance cost C c / C m and replacement to maintenance cost C r / C m on optimal n value.
Water 14 02097 g007
Figure 8. Normalized life cycle cost C L / C r varying with replacement at the survival probability of water mains ( d = 5 % , α = 0.1 , β = 0.6 , and γ = 0.3 ).
Figure 8. Normalized life cycle cost C L / C r varying with replacement at the survival probability of water mains ( d = 5 % , α = 0.1 , β = 0.6 , and γ = 0.3 ).
Water 14 02097 g008
Figure 9. Influences of C c / C r on the optimal survival probability for replacement considering different discounted rates.
Figure 9. Influences of C c / C r on the optimal survival probability for replacement considering different discounted rates.
Water 14 02097 g009
Figure 10. Predictions of occurrence times for the first five breaks for the example water main.
Figure 10. Predictions of occurrence times for the first five breaks for the example water main.
Water 14 02097 g010
Figure 11. Influences of C c / C m and C r / C m on optimal n value for the example water main.
Figure 11. Influences of C c / C m and C r / C m on optimal n value for the example water main.
Water 14 02097 g011
Figure 12. Influences of C c / C r on the optimal survival probability for replacement for the example water main.
Figure 12. Influences of C c / C r on the optimal survival probability for replacement for the example water main.
Water 14 02097 g012
Table 1. Summary of deterioration trends captured by the TTS model with respect to permutations of model parameters (after Lin and Yuan [34]).
Table 1. Summary of deterioration trends captured by the TTS model with respect to permutations of model parameters (after Lin and Yuan [34]).
Case #Model ParameterConditional Intensity, λ (Equation (2))
β + γ γ
1 β + γ > 0 γ > 0 Monotonically increasing
2 γ = 0 Monotonically increasing
3 γ < 0 First decreasing, then increasing
4 β + γ = 0 γ > 0 Monotonically increasing
5 γ = 0 Constant
6 γ < 0 Monotonically decreasing
7 β + γ < 0 γ > 0 First increasing, then decreasing
8 γ = 0 Monotonically decreasing
9 γ < 0 Monotonically decreasing
Table 2. Summary of means, medians, and COVs of the C L / C m for the example based on the simulations.
Table 2. Summary of means, medians, and COVs of the C L / C m for the example based on the simulations.
Para.Replacement at the n -th Failure
1234567891020304050
Mean29.923.220.018.817.917.817.517.617.617.819.020.121.021.7
Median29.322.619.518.417.617.617.117.517.217.518.820.120.821.7
COV0.320.300.300.300.280.270.270.260.260.260.230.220.230.23
Table 3. Summary of means, medians, and COVs of the C L / C r for the example based on the simulations.
Table 3. Summary of means, medians, and COVs of the C L / C r for the example based on the simulations.
ParameterReplacement at the Survival Probability, P s (%)
9080706050403020100
Mean32.929.628.228.729.429.830.730.831.831.8
Median31.828.627.228.228.328.229.729.930.731.1
COV0.290.350.380.370.380.380.370.360.350.33
Table 4. Summary of effects of model parameters and cost ratios on determination of the optimal n and P s value.
Table 4. Summary of effects of model parameters and cost ratios on determination of the optimal n and P s value.
ParameterOptimal n Optimal P s Influence
α NegativePositiveMedium
β NegativePositiveHigh
γ PositiveNegativeMedium
C r / C m PositiveN/AHigh
C c / C m NegativeN/AHigh
C c / C r NegativePositiveHigh
d NegativePositiveMedium
Note: N/A = not applicable.
Table 5. Mean, median, and COV of predictions of the first five breaks for the example water main.
Table 5. Mean, median, and COV of predictions of the first five breaks for the example water main.
Parameter1st2nd3rd4th5th
Mean (Year)103.0148.5185.3217.0241.4
Median (Year)64.4118.0154.9194.6218.1
COV1.0550.7960.6780.6050.546
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, P.; Chen, X.; Huang, S.; Ma, B. An Optimal Maintenance and Replacement Strategy for Deteriorating Water Mains. Water 2022, 14, 2097. https://doi.org/10.3390/w14132097

AMA Style

Lin P, Chen X, Huang S, Ma B. An Optimal Maintenance and Replacement Strategy for Deteriorating Water Mains. Water. 2022; 14(13):2097. https://doi.org/10.3390/w14132097

Chicago/Turabian Style

Lin, Peiyuan, Xianying Chen, Sheng Huang, and Baosong Ma. 2022. "An Optimal Maintenance and Replacement Strategy for Deteriorating Water Mains" Water 14, no. 13: 2097. https://doi.org/10.3390/w14132097

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