Water Pipe Replacement Scheduling Based on Life Cycle Cost Assessment and Optimization Algorithm

: Water distribution networks (WDNs) comprise a complex network of pipes and are crucial for providing potable water to urban communities. Therefore, WDNs must be carefully managed to avoid problems such as water contamination and service failures; however, this requires a large budget. Because WDN components have different statuses depending on their installation year, location, transmission pressure, and ﬂow rate, it is difﬁcult to plan the rehabilitation schedule within budgetary constraints. This study, therefore, proposes a new pipe replacement scheduling approach to smooth the investment time series based on a life cycle cost (LCC) assessment for a large-scale WDN. The proposed scheduling plan simultaneously considers both the annual budget limitation and the optimum expenditure on the useful life of pipes. A multi-objective optimization problem consisting of three decision-making objectives—minimum imposed LCC on the network, minimum standard deviation of annual investment, and minimum average age of the network— is thus solved using a nondominated sorting genetic algorithm to obtain an optimal plan. Three scenarios with different pipe replacement time spans and different annual budget constraints are considered accordingly. The results indicate that the proposed scheduling framework provides an efﬁcient water pipe replacement scheduling plan with a smooth management budget.


Introduction
A water distribution network (WDN) comprises a complex network of underground pipes and plays a key role in delivering potable water to urban communities [1]. WDN is one of the most expensive components of the overall water infrastructure and directly impacts people's quality of life and economic prosperity [2]. Unfortunately, water infrastructure managers worldwide are now facing the problem of aging pipes in WDNs, and limited historical data have hampered cost-efficient replacements. Furthermore, low budgets and underinvestment in water infrastructure over recent decades have led to the degradation of WDNs. Therefore, replacing, designing, and optimally managing WDNs is expected to be a significant social challenge in the 21st century [3]. Pipelines age because of various internal and external factors that lead to a decrease in their functionality and an increase in their risk of failure; both of these issues can incur large social and economic costs [4]. Thus, pipeline rehabilitation must be scheduled effectively to ensure adequate water quality and structural performance [5]. In practice, suitable pipe replacement/rehabilitation requires a large budget. However, as noted above, budgets for rehabilitating water infrastructure are typically limited, thereby constraining the options available to water infrastructure managers. Thus, although water pipelines are aging, their replacement and rehabilitation rate is lagging. For a considerable proportion of WDN assets, an efficient and practical pipe replacement scheduling plan must satisfy the annual budget and minimize future costs while maintaining WDN functionality. Water infrastructure managers require sustainable, practical, and affordable solutions to address required investment to peak in some years when the number of pipes needing replacement exceeds the budget. Indeed, though water infrastructure governors need to start putting more effort and funding towards fixing the problem of aging pipes in WDNs, over the last decade, most studies have proposed scheduling plans without considering budget restrictions (i.e., expenditure ceilings). In this study, the LCC is therefore used to propose an economical replacement plan for individual pipes in a real large-scale WDN by considering annual budget restrictions to smooth the annual investment time series based on optimization models. This approach considers three objective functions: Imposed LCC minimization, annual investment smoothing, and network age minimization. To the best of our knowledge, no previous study has investigated the use of smoothing to avoid investment peaks by considering all individual assets. In the proposed plan, the replacement time is relaxed from peak years to off-peak years while maintaining WDN reliability. The proposed method can be used to schedule the rehabilitation of a wide range of WDN assets and can determine the optimum number of service years of an individual asset by considering the above-mentioned objective functions. Thus, this method can enable the development of a more realistic budget for WDN rehabilitation.
The remainder of this paper is structured as follows. Section 2 describes the proposed approach and its methodology. Section 3 presents the optimization results obtained under three different scenarios. Finally, Section 4 discusses the application of the proposed method and presents the conclusions of this study.

Optimal Water Pipe Replacement Based on LCC Assessment
LCC is one of the most critical factors in cost-effective decision-making for a system [30]. The LCC evaluation was thus employed in this study to achieve optimal pipe replacement scheduling. The LCC of a pipe is simply defined as the sum of all costs incurred during its entire life span. In general, the LCC for each pipe attribute is therefore calculated as the sum of acquisition, maintenance, replacement, repair, and disposal costs. In this study, the LCC is calculated as LCC(D, t) = CI(D, t) + CR(D, t), where CI is the initial cost ($/km/year), CR is the running cost ($/km/year), D is pipe diameter (mm), and t is pipe replacement interval (year). Here, the CI includes only the pipe replacement cost per year and the CR includes only the pipe repair cost per year. The CI is calculated using: where CP is the pipe replacement cost ($/km). The CR is the product of the repair cost and the average failure rate during the replacement interval, and is therefore expressed as where Cr is the pipe repair cost ($/failure) referred from [31] and expressed by Cr(D) = 1.300 D 304.8 0.62 × 800, (4) and Fr is pipe failure rate (failure/km/year), which is dependent on the pipe diameter and pipe age (A). This study uses the failure rate obtained by [4] as follows: Available records of historical failure data must be used to estimate the WDN failure rate. This failure rate was derived through nonlinear regression analysis using the Statistical Package for the Social Sciences (SPSS, v. 18) and the input model suggested by [16]. This input model explains the relationship between the failure rate and the pipe attributes (i.e., diameter and age). Figure 1 illustrates an example of an LCC curve, indicating that as the replacement age of a pipe increases, the running cost increases, and the initial cost decreases. The optimal time for replacement, or the most economic life span, is that for which the total cost is the minimum (t * ), shown in Figure 1. In this study, the least life cycle cost (LLCC) for individual pipes is defined as the LCC at t * . Available records of historical failure data must be used to estimate the WDN failure rate. This failure rate was derived through nonlinear regression analysis using the Statistical Package for the Social Sciences (SPSS, v. 18) and the input model suggested by [16]. This input model explains the relationship between the failure rate and the pipe attributes (i.e., diameter and age). Figure 1 illustrates an example of an LCC curve, indicating that as the replacement age of a pipe increases, the running cost increases, and the initial cost decreases. The optimal time for replacement, or the most economic life span, is that for which the total cost is the minimum ( * ), shown in Figure 1. In this study, the least life cycle cost (LLCC) for individual pipes is defined as the LCC at * . The summation of costs for all individual pipes at t * forms the LLCC of the network (LLCCN), which represents the cost at which the summation of the initial and running costs is the minimum, and is defined as where n is the number of pipes in the network and Li is the length of pipe i.
The LCC for each pipe attribute (i.e., diameter and material) is calculated to determine the most economic replacement age for each individual pipe. However, to smooth the annual investment time series, the replacement time of each pipe should be relaxed within a time window around the optimal point ( * ) as shown in Figure 2.
Notably, the LCC curve can help engineers and decision-makers choose an appropriate and economical replacement policy [32]. The advantage of using the LCC curve is that it approximates the optimal economic replacement age and shows the change in total cost around this optimal point. If the total cost curve is flat around the optimum, as shown in Figure 2, the engineer need not plan for the replacement to be performed precisely at the optimum age ( * ), thus indicating some flexibility in the replacement scheduling plan. Therefore, as shown in Figure 2, performing the replacement within the feasible replacement time boundary of ± (bounded by the lower boundary (LB) and upper boundary (UB)) has no significant effect on the total cost. By contrast, if the LCC curve is not flat around the optimal point and changes sharply on both sides, the replacement should be performed as close to the optimal age as possible. Note that Figure 2 shows that a change in replacement time will lead to an increase in the LCC whether the replacement time is moved forward or backward; however, moving the replacement time will impact the overall system age (i.e., pipe age). The summation of costs for all individual pipes at t * forms the LLCC of the network (LLCC N ), which represents the cost at which the summation of the initial and running costs is the minimum, and is defined as where n is the number of pipes in the network and L i is the length of pipe i. The LCC for each pipe attribute (i.e., diameter and material) is calculated to determine the most economic replacement age for each individual pipe. However, to smooth the annual investment time series, the replacement time of each pipe should be relaxed within a time window around the optimal point (t * ) as shown in Figure 2.

Smoothing the Investment Time Series
An emerging challenge worldwide is the need to replace the vast numbers and lengths of deteriorating pipes beneath streets, as most pipes are nearing the end of their useful lives. According to the LCC concept, pipe replacements accrue towards the end of the useful life of the pipes. The optimal economic replacement age is when the replacement cost becomes lower than the failure cost (i.e., costs of burst and leaks).
Based on the LCC formulation discussed in Section 2.1, the investment in each year Notably, the LCC curve can help engineers and decision-makers choose an appropriate and economical replacement policy [32]. The advantage of using the LCC curve is that it approximates the optimal economic replacement age and shows the change in total cost around this optimal point. If the total cost curve is flat around the optimum, as shown in Figure 2, the engineer need not plan for the replacement to be performed precisely at the optimum age (t * ), thus indicating some flexibility in the replacement scheduling plan. Therefore, as shown in Figure 2, performing the replacement within the feasible replacement time boundary of ±α (bounded by the lower boundary (LB) and upper boundary (UB)) has no significant effect on the total cost. By contrast, if the LCC curve is not flat around the optimal point and changes sharply on both sides, the replacement should be performed as close to the optimal age as possible. Note that Figure 2 shows that a change in replacement time will lead to an increase in the LCC whether the replacement time is moved forward or backward; however, moving the replacement time will impact the overall system age (i.e., pipe age).

Smoothing the Investment Time Series
An emerging challenge worldwide is the need to replace the vast numbers and lengths of deteriorating pipes beneath streets, as most pipes are nearing the end of their useful lives. According to the LCC concept, pipe replacements accrue towards the end of the useful life of the pipes. The optimal economic replacement age is when the replacement cost becomes lower than the failure cost (i.e., costs of burst and leaks).
Based on the LCC formulation discussed in Section 2.1, the investment in each year (I y ) is the summation of the cost (CI and CR) of each pipe in a WDN; this leads to the formation of an investment time series. This investment time series imposes the minimum LCC on the network as well as the individual pipe if all pipes are replaced at their t * , which is determined based on the attribute (i.e., diameter and material) of each pipe. However, in some years, the required investment may exceed the annual budget, as shown in Figure 3, making the scheduling plan infeasible. To ensure that the scheduling plan remains feasible and close to actual practice, the investment time series needs to be smoothed to avoid investment peaks and meet the available annual budget. In order to do so, the investment time series needs to be changed so that investments in peak years are moved to off-peak years by changing the replacement time points for all individual pipes in the feasible replacement time boundary. The optimization algorithm (NSGA-II) was thus used in this study to obtain a set of replacement intervals for all pipes among those within the feasible replacement time boundary and smooth the investment time series as described in Section 2.3.

Smoothing the Investment Time Series
An emerging challenge worldwide is the need to replace the vast numbers and lengths of deteriorating pipes beneath streets, as most pipes are nearing the end of their useful lives. According to the LCC concept, pipe replacements accrue towards the end of the useful life of the pipes. The optimal economic replacement age is when the replacement cost becomes lower than the failure cost (i.e., costs of burst and leaks).
Based on the LCC formulation discussed in Section 2.1, the investment in each year ( ) is the summation of the cost (CI and CR) of each pipe in a WDN; this leads to the formation of an investment time series. This investment time series imposes the minimum LCC on the network as well as the individual pipe if all pipes are replaced at their * , which is determined based on the attribute (i.e., diameter and material) of each pipe. However, in some years, the required investment may exceed the annual budget, as shown in Figure 3, making the scheduling plan infeasible. To ensure that the scheduling plan remains feasible and close to actual practice, the investment time series needs to be smoothed to avoid investment peaks and meet the available annual budget. In order to do so, the investment time series needs to be changed so that investments in peak years are moved to off-peak years by changing the replacement time points for all individual pipes in the feasible replacement time boundary. The optimization algorithm (NSGA-II) was thus used in this study to obtain a set of replacement intervals for all pipes among those within the feasible replacement time boundary and smooth the investment time series as described in Section 2.3.

Overview of NSGA-II
This study adopted the well-known and widely used multi-objective NSGA-II [33] as the optimization method. Figure 4 shows the optimization flowchart for NSGA-II and the general steps are briefly described below.

1.
NSGA-II randomly generates an initial population consisting of a number of chromosomes. Each chromosome is the value of the design parameter (i.e., the replacement time of an individual pipe).

2.
Each individual in the population generated in the previous step is ranked based on an evaluation of the objective function. Then, the individuals are sorted based on their rank by using the crowding distance.

3.
Parents are randomly selected from a mating pool to create the new generation. The mating pool consists of solutions with a higher crowding distance and rank. 4.
To produce offspring, parents from the previous step undergo crossover and mutation procedures. In the crossover procedure, two parents breed to produce an offspring that inherits its genes from both parents. In the mutation procedure, some values of the genes in each offspring are changed, thus providing the offspring with the opportunity to have at least one different gene value than their parents.
The previous steps are repeated until convergence is reached. In this study, the optimization algorithm was terminated after a defined number of generations. The efficiency of this algorithm can be improved by tuning its control parameters, such as the crossover and mutation rates, number of generations, and population size. To achieve the highest efficiency, various sets of control parameters were examined in this study.
This study adopted the well-known and widely used multi-objective NSGA-II [33] as the optimization method. Figure 4 shows the optimization flowchart for NSGA-II and the general steps are briefly described below. 1. NSGA-II randomly generates an initial population consisting of a number of chromosomes. Each chromosome is the value of the design parameter (i.e., the replacement time of an individual pipe).
2. Each individual in the population generated in the previous step is ranked based on an evaluation of the objective function. Then, the individuals are sorted based on their rank by using the crowding distance.
3. Parents are randomly selected from a mating pool to create the new generation. The mating pool consists of solutions with a higher crowding distance and rank.
4. To produce offspring, parents from the previous step undergo crossover and mutation procedures. In the crossover procedure, two parents breed to produce an offspring that inherits its genes from both parents. In the mutation procedure, some values of the genes in each offspring are changed, thus providing the offspring with the opportunity to have at least one different gene value than their parents.
The previous steps are repeated until convergence is reached. In this study, the optimization algorithm was terminated after a defined number of generations. The efficiency of this algorithm can be improved by tuning its control parameters, such as the crossover and mutation rates, number of generations, and population size. To achieve the highest efficiency, various sets of control parameters were examined in this study.

Construction of Optimization Model
This study focused on developing a smooth, optimized replacement scheduling plan in which the peaks of investment time series are dispersed across a time span to meet the budget limit and use the available budget efficiently, while the LCC of each pipe in a WDN is kept as close as possible to its minimum value. The decision variable of the optimization

Construction of Optimization Model
This study focused on developing a smooth, optimized replacement scheduling plan in which the peaks of investment time series are dispersed across a time span to meet the budget limit and use the available budget efficiently, while the LCC of each pipe in a WDN is kept as close as possible to its minimum value. The decision variable of the optimization model is a vector containing the optimal replacement time for all individual pipes, and the replacement time should be determined within a replacement time boundary (t * ± α). The three objective functions of the problem are given by Min where E[·] represents the expectation operator, I t y is a set of the annual investment, A t y is the annual average age of pipes, and n is the number of pipes in the network. Equation (7) minimizes the imposed LCC of the network, that is, the difference between LCC N (i.e., LCC of the network) and LLCC N . By applying this objective function, the optimized LCC N is kept as close as possible to LLCC N , implying that the imposed cost is minimized. Equation (8) minimizes the standard deviation (SD) of the annual investments to smooth the investment time series. Equation (9) minimizes the average pipe age in the WDN to more reliably maintain the network [34].
The two main constraints considered in this study are the annual budget and the allowable replacement time interval (t * ± α), during which the pipe replacement is performed. Considering the minimization of the three objective functions in Equations (7)-(9), the multi-objective optimization problem can be formulated as: where β is the available annual budget constraint (the total investment in each year must be equal or less than this value) and α is the replacement time boundary.

Case Study and Assumptions
The proposed method was applied to a currently operating WDN in South Korea that consists of 6326 pipes with diameters ranging from 25-500 mm, as shown in Figure 5. Here, the total length of pipes is 812.5 km, and the percentages of pipe lengths in each diameter range are 57.4% (D < 80 mm), 32.6% (80 mm ≤ D < 200 mm), 7.1% (200 mm ≤ D < 300 mm), 2.6% (300 mm ≤ D < 400 mm), and 0.3% (400 mm ≤ D). In this study, only pipes with diameters of 80-500 mm were considered, which constituted 42.6% of the total length of all pipes and counted as 3042 pipes. Note that all pipes are made of the same material (i.e., ductile iron). The pipes can be divided into four groups according to their age in the year 2020, as shown in Figure 5b; 65.4% of the pipes are 10-20 years old.
It should be noted that the following assumptions were made in this case study: 1.
All pipes need to be replaced at least once in the simulation time horizon in consideration of the diameter and the first installation year of the pipes.

2.
Pipes that already passed their replacement time (t * ) will be replaced on priority in the first year of the scheduling plan. It should be noted that the following assumptions were made in this case study: 1. All pipes need to be replaced at least once in the simulation time horizon in consideration of the diameter and the first installation year of the pipes.
2. Pipes that already passed their replacement time ( * ) will be replaced on priority in the first year of the scheduling plan. The initial costs of ductile iron pipe were obtained from the 'Water Facilities Construction Cost Estimation Report', published by K-water [35], and are summarized in Table 1. Table 2 shows the LLCC and corresponding * for different pipe diameters. The LCC curves are shown in Figure 6, in which each color represents a different pipe diameter. Table 1. Initial cost data of ductile iron pipe [35].

Pareto Front Obtained by Multiobjective Optimization
The NSGA-II was used to obtain the optimal replacement scheduling plan. The optimization model was developed in Python using the Pymoo library [36] and compiled with Cython 3 for parallel processing. For better convergence, the population size was set to 2000. The number of surviving offspring was set to 1500 to retain the best solutions in the population. Different generations were used for different scenarios in NSGA-II to obtain an extensive Pareto front and guarantee convergence. A simulated binary crossover function and a polynomial mutation operator with probabilities of 0.9 and 0.1, respectively, were used. The final run of the program with 2000 generations and population size of 2000 required approximately 99 h on an Intel ® Core i9-10920X CPU @ 3.5 GHz with 128 GB memory. Three decision-making scenarios were considered during the optimization procedure as described below.
In the first scenario (Scenario 1), the feasible replacement time boundary ( ) for changing the replacement time was set to 5 years, which means that the new replacement time could be in the interval of * − 5 to * + 5 years. The annual budget constraint was set to 2.5 M$ for the first scenario. In the second scenario (Scenario 2), the boundary was set to * ± 10 years, and the annual budget constraint was 2.2 M$. In the third scenario (Scenario 3), the range was set to * ± 16 years, and the annual budget constraint was 2.0 M$. The replacement time in all three scenarios lies within the LB and UB; this range was chosen based on the obtained LCC curve. As explained in Section 2.1, in the flat area of an LCC curve, the change in annual average cost is negligible, and the cost increases marginally when the replacement time is shifted to a different year to smooth the investment.

Pareto Front Obtained by Multiobjective Optimization
The NSGA-II was used to obtain the optimal replacement scheduling plan. The optimization model was developed in Python using the Pymoo library [36] and compiled with Cython 3 for parallel processing. For better convergence, the population size was set to 2000. The number of surviving offspring was set to 1500 to retain the best solutions in the population. Different generations were used for different scenarios in NSGA-II to obtain an extensive Pareto front and guarantee convergence. A simulated binary crossover function and a polynomial mutation operator with probabilities of 0.9 and 0.1, respectively, were used. The final run of the program with 2000 generations and population size of 2000 required approximately 99 h on an Intel ® Core i9-10920X CPU @ 3.5 GHz with 128 GB memory. Three decision-making scenarios were considered during the optimization procedure as described below.
In the first scenario (Scenario 1), the feasible replacement time boundary (α) for changing the replacement time was set to 5 years, which means that the new replacement time could be in the interval of t * − 5 to t * + 5 years. The annual budget constraint was set to 2.5 M$ for the first scenario. In the second scenario (Scenario 2), the boundary was set to t * ± 10 years, and the annual budget constraint was 2.2 M$. In the third scenario (Scenario 3), the range was set to t * ± 16 years, and the annual budget constraint was 2.0 M$. The replacement time in all three scenarios lies within the LB and UB; this range was chosen based on the obtained LCC curve. As explained in Section 2.1, in the flat area of an LCC curve, the change in annual average cost is negligible, and the cost increases marginally when the replacement time is shifted to a different year to smooth the investment. Moreover, the chosen replacement time boundary depends on the budget constraint. Specifically, a wider time range (i.e., 16 years in the third scenario) leads to a smoother investment time series; by contrast, a narrower range (i.e., 5 years in the first scenario) leads to a small deviation from the LLCC N and, consequently, a less smooth investment time series that is likely to violate the budget constraint. Figure 7 shows the nondominated Pareto solution obtained using NSGA-II for the three evaluated scenarios with different annual budget constraints and replacement time boundaries. The Pareto solutions of the first, second, and third scenarios are indicated in green, blue, and navy, respectively.
To analyze the characteristics of the optimal schedule, four representative solutions, including three corners and one knee-point, were selected from among all solutions for each of the three Pareto front scenarios, and are indicated by the red plus signs in Figure 7. The solutions at the three selected corners provided the minimum imposed LCC, minimum SD of annual investment, and minimum average WDN age among all solutions, and the knee-point provided the balanced solution for all three objectives. The knee-point was extracted after normalizing all three axes, after which it was the closest point to the origin. To increase the benefit of one of the three objectives, the other two objectives must clearly be sacrificed. Figure 7 shows the trade-off relationships among the three evaluated objectives. 7. The solutions at the three selected corners provided the minimum imposed LCC, minimum SD of annual investment, and minimum average WDN age among all solutions, and the knee-point provided the balanced solution for all three objectives. The knee-point was extracted after normalizing all three axes, after which it was the closest point to the origin. To increase the benefit of one of the three objectives, the other two objectives must clearly be sacrificed. Figure 7 shows the trade-off relationships among the three evaluated objectives. To better clarify the relationship among these three objectives, the three-dimensional relationships in Figure 7 are projected in two-dimensional planes in Figure 8. To better clarify the relationship among these three objectives, the three-dimensional relationships in Figure 7 are projected in two-dimensional planes in Figure 8.
Referring to Figures 7 and 8, it can be observed that the Pareto fronts resulting from this multi-objective pipe replacement scheduling problem indicate that a trade-off relationship exists among all three objectives, with the trade-off between the imposed LCC and SD of annual investment being the strongest, followed by that between the average WDN age and the imposed LCC. As the SD of the annual investment increases, the average WDN age increases and imposed LCC decreases, implying that the investment time series becomes less smooth, and the annual investment fluctuation increases. Comparing the Pareto fronts of all three scenarios, with an increasing replacement time boundary, the overall imposed LCC increases, but the average WDN age and variation of the annual investment decrease. Further, it can be observed that the Pareto front widens from Scenarios 1 to 3; as the replacement time boundary is more relaxed, the decision space increases, resulting in a larger number of alternative optimal solutions.
Comparing Scenarios 1 and 3 in Figures 7 and 8, it can be seen that allowing the wider replacement time boundary and reducing the annual budget limit broadens the change in the imposed LCC, causing a wide variety of optimal solutions to appear in the Pareto front. In Scenario 1, which has the narrowest replacement time boundary and the highest annual budget, the replacement time does not change significantly from t * , thus LCC is minimum, but the annual variation and system age are the highest. Whereas in Scenario 3, which has the widest replacement time boundary and the lowest annual budget, the annual investment shows the smallest SD, implying that the investment time series is the smoothest. The results from all three scenarios indicate the presence of a trade-off relationship between the variation of the annual investment and the imposed LCC. A smoother investment time series leads to an increase in the imposed LCC, because the replacement is conducted sooner or later than t * . To better clarify the relationship among these three objectives, the three-dimensional relationships in Figure 7 are projected in two-dimensional planes in Figure 8. Referring to Figures 7 and 8, it can be observed that the Pareto fronts resulting from this multi-objective pipe replacement scheduling problem indicate that a trade-off relationship exists among all three objectives, with the trade-off between the imposed LCC and SD of annual investment being the strongest, followed by that between the average WDN age and the imposed LCC. As the SD of the annual investment increases, the average WDN age increases and imposed LCC decreases, implying that the investment time series becomes less smooth, and the annual investment fluctuation increases. Comparing the Pareto fronts of all three scenarios, with an increasing replacement time boundary, the overall imposed LCC increases, but the average WDN age and variation of the annual investment decrease. Further, it can be observed that the Pareto front widens from Scenarios 1 to 3; as the replacement time boundary is more relaxed, the decision space increases, resulting in a larger number of alternative optimal solutions.   Figure 9 shows the investment time series for the WDN before smoothing. The value of investment for a specific year is calculated as the product of the replacement cost (initial cost) and the length of pipes replaced within that particular year; for the rest of the pipes, the investment value is calculated as the product of the running cost and the pipe length. In Figure 9, the replacement and running costs are indicated in navy and blue, respectively.
In some years in Figure 9, a majority of pipes need to be replaced (e.g., 2053), and the share of replacement cost is high enough to jeopardize the replacement schedule. Moreover, there are some years in Figure 9 (e.g., 2026) in which no replacement is planned. Figure 9 shows the investment time series for the WDN before smoothing. The value of investment for a specific year is calculated as the product of the replacement cost (initial cost) and the length of pipes replaced within that particular year; for the rest of the pipes, the investment value is calculated as the product of the running cost and the pipe length. In Figure 9, the replacement and running costs are indicated in navy and blue, respectively. In some years in Figure 9, a majority of pipes need to be replaced (e.g., 2053), and the share of replacement cost is high enough to jeopardize the replacement schedule. Moreover, there are some years in Figure 9 (e.g., 2026) in which no replacement is planned. In these cases, moving items in the replacement queue from peak years to off-peak years will relieve budgetary pressure and lead to a reasonable and sensible replacement plan. In the next sections (3.3.2-3.3.4), obtaining a smoother investment time series for the replacement schedule is explained in detail for the three scenarios. For a better comparison of the results in each of the three evaluated scenarios, details of the network before smoothing are shown in Table 3. Here, the total annual investment (TAI) is the total cost divided by the planning time horizon (i.e., 108 years). The investment time series obtained by the four representative solutions of the first scenario are shown in Figure 10, and the details are summarized in Table 4. In Figure 10, the red dashed line indicates the annual budget limit. The first scheduling plan (Solution 1, shown in Figure 10a) shows the smoothest time series with the minimum SD of annual investment. The second scheduling plan (Solution 2, shown in Figure 10b) is the closest to LLCCN and therefore has the lowest imposed cost. The third scheduling plan (Solution 3, shown in Figure 10c) yields the youngest network with the lowest pipe average age, which can be considered as a more reliable plan than other solutions. The use of the knee-point (Solution 4) provides a balanced scheduling plan that considers all three objective functions simultaneously, as shown in Figure 10d.

Future Investment Based on LLCCN
In Table 4, MODE represents the pipe replacement year relative to the * with the highest frequency; that is, the year picked most often as the new replacement time for In these cases, moving items in the replacement queue from peak years to off-peak years will relieve budgetary pressure and lead to a reasonable and sensible replacement plan. In the next (Sections 3.3.2-3.3.4), obtaining a smoother investment time series for the replacement schedule is explained in detail for the three scenarios. For a better comparison of the results in each of the three evaluated scenarios, details of the network before smoothing are shown in Table 3. Here, the total annual investment (TAI) is the total cost divided by the planning time horizon (i.e., 108 years). Table 3. Details of the investment before smoothing.  The investment time series obtained by the four representative solutions of the first scenario are shown in Figure 10, and the details are summarized in Table 4. In Figure 10, the red dashed line indicates the annual budget limit. The first scheduling plan (Solution 1, shown in Figure 10a) shows the smoothest time series with the minimum SD of annual investment. The second scheduling plan (Solution 2, shown in Figure 10b) is the closest to LLCC N and therefore has the lowest imposed cost. The third scheduling plan (Solution 3, shown in Figure 10c) yields the youngest network with the lowest pipe average age, which can be considered as a more reliable plan than other solutions. The use of the knee-point (Solution 4) provides a balanced scheduling plan that considers all three objective functions simultaneously, as shown in Figure 10d.
In Table 4, MODE represents the pipe replacement year relative to the t * with the highest frequency; that is, the year picked most often as the new replacement time for pipes in the network. The maximum I t year represents the maximum annual investment value. The objective of this study was to minimize and disperse this maximum annual investment compared to the LLCC N schedule to provide the smallest fluctuation in annual expenditure. The running and initial costs and their summation after changing the replacement time represent the change in expenditure due to the selection of a different budget and replacement boundary. When the replacement time is postponed by 5 years from the expected t * in the first scenario, a 5-year running cost is added to the new plan, therefore 5 years are added to its life cycle. On the other hand, when a pipe is replaced 5 years sooner than the expected t * , 5 years are omitted from its life cycle, which reduces the running cost and pipe age, but increases the replacement cost. The TAI, defined as the total cost divided by the time horizon (i.e., 108 years), provides an estimation of the annual WDN maintenance cost and enables a more realistic budget to be assigned for WDN rehabilitation. Note that the LCC does not depend on the time scale and represents the continuous life cycle of an asset. However, the total cost depends on the time horizon to be investigated; upon changing the time horizon, the total cost changes because the cycle is cut in a particular year. In this case, only the cost up to that specific year is counted. Note that the time horizon of the case study WDN was determined based on the above-mentioned assumption in which all pipes need to be replaced at least once. In Scenario 2, the budget limitation decreased to 2.2 M$, and as a consequence, the replacement boundary was increased to ±10 years. In this scenario, as shown in Table 5, the MODE for Solution 1 was −10 years, and the replacement times were dispersed to the extent possible to smooth the time series. Consequently, Figure 11a shows a wellsmoothed investment time series for Solution 1. The imposed cost for this solution was the highest, while the maximum annual investment was the lowest, leading to a 66.25%  As shown in Table 4, for Scenario 1 (in which the replacement boundary was ±5 years), the replacement times for the majority of pipes at the edge of the boundary were determined to occur 5 years sooner in order to disperse all the peaks and provide a sufficiently smooth time series.
As shown in Figure 10a, Solution 1 provides the smoothest investment time series with the lowest maximum annual investment of 1.99 M$. This solution shows a 59% decrease in the SD of annual investment compared to the before-smoothing-plan. For Solution 2, the majority of pipe replacement times remained unchanged in order to ensure a small deviation from LLCC N , and the MODE is 0. This solution was obtained to prioritize keeping the replacement time the same as t * (MODE = 0) while dispersing the peaks as necessary. As shown in Figure 10b, this time series follows the same pattern as the time series before smoothing ( Figure 9), but in a smoother manner, and imposed the minimum cost to the network compared to the other solutions, with a maximum annual investment of 2.22 M$, and the LCC marginally increased by 0.08% compared to the before-smoothingplan. For Solution 3, shown in Figure 10c, the priority was to keep the network younger; as a consequence, the majority of pipes are replaced 5 years sooner than t * , which results in a more reliable plan than provided by the other three solutions. The average age in this solution decreased by 14.7% compared to the before-smoothing-plan. In Solution 4, which used the knee-point to fairly consider all three objectives, the replacement time can be observed to be dispersed evenly around the provided boundary, as shown in Figure 10d In Scenario 2, the budget limitation decreased to 2.2 M$, and as a consequence, the replacement boundary was increased to ±10 years. In this scenario, as shown in Table 5, the MODE for Solution 1 was −10 years, and the replacement times were dispersed to the extent possible to smooth the time series. Consequently, Figure 11a shows a wellsmoothed investment time series for Solution 1. The imposed cost for this solution was the highest, while the maximum annual investment was the lowest, leading to a 66.25% smaller SD compared to the before-smoothing-plan. Solution 2 was used to obtain the most economical scheduling plan by keeping the imposed LCC near LLCC N (MODE = 0) and consequently a less smooth curve showing the highest maximum annual investment of 2.18 M$ can be observed in Figure 11b. It is observed that the overall investment of Scenario 2 marginally increases compared to Scenario 1, because the replacement time boundary was relaxed by five additional years. The youngest scheduling plan was provided by Solution 3, which showed a 23.53% reduction compared to the before-smoothing-plan and is shown in Figure 11c.  Figure 11b. It is observed that the overall investment of Scenario 2 marginally increases compared to Scenario 1, because the replacement time boundary was relaxed by five additional years. The youngest scheduling plan was provided by Solution 3, which showed a 23.53% reduction compared to the before-smoothing-plan and is shown in Figure 11c. Finally, Figure 11d shows the future investment time series for Solution 4, which used the knee-point showing the balanced values for the applied three objectives as listed in Table 5. Finally, Figure 11d shows the future investment time series for Solution 4, which used the knee-point showing the balanced values for the applied three objectives as listed in Table 5.

Future Investment Based on Scenario 3
Scenario 3 dispersed the replacement time interval the widest (±16 years) with the lowest budget limitation of 2.0 M$. As shown in Figure 12a, the time series for Solution 1 was the smoothest, not only among the four representative solutions in this scenario, but also among all three scenarios, showing a 73.2% reduction in the SD of annual investment compared to the before-smoothing-plan. In Solution 2, the imposed LCC is the minimum (with MODE = 0) while the SD is the maximum as shown in Table 6 and Figure 12b. The most reliable and youngest plan was provided by Solution 3, as shown in Figure 12c. It is interesting to observe that the MODE of Solution 3 is −12 years, which do not reach the LB (i.e., −16 years). This means that the majority of pipes should be replaced later than the LB to refrain from violating the other objectives while minimizing the system age. This indicates that the proposed method acted intelligently in selecting the replacement time and showed reasonable and sensible results. The average age of the network decreased by 26.8% compared to the before-smoothing-plan.
(with MODE = 0) while the SD is the maximum as shown in Table 6 and Figure 12b. The most reliable and youngest plan was provided by Solution 3, as shown in Figure 12c. It is interesting to observe that the MODE of Solution 3 is −12 years, which do not reach the LB (i.e., −16 years). This means that the majority of pipes should be replaced later than the LB to refrain from violating the other objectives while minimizing the system age. This indicates that the proposed method acted intelligently in selecting the replacement time and showed reasonable and sensible results. The average age of the network decreased by 26.8% compared to the before-smoothing-plan.   Table 7 shows a comparison between all three evaluated scenarios in terms of the minimum value of each objective function. First, a notable trade-off can be observed between two constraints, that is, increasing the replacement time boundary could reduce the budget ceiling. Setting the replacement time boundary affects the ability to meet the   Table 7 shows a comparison between all three evaluated scenarios in terms of the minimum value of each objective function. First, a notable trade-off can be observed between two constraints, that is, increasing the replacement time boundary could reduce the budget ceiling. Setting the replacement time boundary affects the ability to meet the budget constraints and vice versa. Next, the reason for evaluating these three scenarios was to assess the effect of different budgets and time boundaries in terms of the three objective functions. Scenario 1, with the highest budget limitation and the smallest replacement boundary, provided the lowest-cost scheduling plan in which the imposed LCC was negligible, while the investment time series was fairly smooth due to the small change in replacement time. Whilst, Scenario 3, with the lowest budget limitation and widest replacement range, provided the smoothest and the youngest scheduling plan while imposing the highest LCC on the WDN. The change in the average age of the network from Scenarios 2 to 3 was negligible, while the change in imposed LCC was sufficiently high to make a competitive difference in selecting between these scenarios. On the other hand, there was no significant change in the imposed LCC from Scenarios 1 to 2, while the average age in Scenario 2 implied that its scheduling plan would maintain a reliable WDN by keeping the assets younger. Therefore, Scenario 2 provided a reasonable case that comprehensively considered the three objectives, balancing the smoothness, imposed LCC, and average age of the WDN. This indicates that the proposed method allows decisionmakers to consider various scenarios with combinations of annual budget and replace time span and select the most appropriate solution as required. The optimal scheduling plans obtained in Scenarios 1, 2, and 3 increased the LCC marginally by 0.08%, 0.27%, and 1.06%, respectively, from the LLCC. In terms of smoothness, a well-prorated investment time series was obtained with small fluctuations in annual investment by dispersing the peak annual investment of 3.36 M$ and smoothing to 1.99, 1.80, and 1.79 M$ for Solutions 1 in Scenarios 1, 2, and 3, respectively. In all three scenarios, moving from Solution 1 to Solution 2 led to an increase in the year-to-year fluctuation in investment but reduced the imposed LCC. Solutions 1 and 3 provided the highest initial cost and, consequently, the lowest running cost.
Comparing Figures 10-12 with Figure 9, the investment time series clearly becomes smoother, and the fluctuation in annual investment decreases due to the optimal scheduling of pipe replacements. In addition, it can be observed that among the four representative solutions evaluated, Solution 1 provided a well-prorated plan that can make the budget allocations easier than in Solutions 2 and 3. Budget allocation is an essential issue in water infrastructure asset management that increases system reliability by considering all the assets in the budget allocation process. Thus, using the proposed method will help water infrastructure managers by improving their ability to assign a realistic budget that better matches the reality of the system and the financial structure by which it is maintained.

Conclusions
The NSGA-II was tuned to develop a multi-objective optimization model used to solve a pipe replacement scheduling problem of WDN. Three scenarios with different replacement time boundaries and budget limitations were defined by targeting the minimum imposed LCC, minimum SD of the annual investment, and minimum average age of the WDN. For individual pipes in a real-life WDN, the results of the simulation proposed the replacement times, and four scheduling plans were evaluated based on four different budgetary viewpoints and decision-maker opinions. We confirmed that the optimization provides a helpful decision-making result to visualize the trade-offs among pipe replacement schedules. The following conclusions were derived based on the proposed scheduling plan.
Using the proposed scheduling framework, a useful number of service years was obtained for all individual pipes, and rehabilitation management was performed more efficiently. For a large WDN with ductile iron pipe diameters of 80-500 mm in an age range of 4-53 years, smoothed investment time series were proposed using three different scenarios. We ranked four different solutions for each scenario based on a real field situation and decision-maker opinions, obtaining twelve alternative scheduling plans that better match reality.
The running cost used in this study plays an essential role in developing the investment time series; therefore, data analytics and techniques for the control and management of a large-scale WDN are needed. By using online monitoring and recording failure data, the pipe failure rate could be continuously updated. Thus, by updating the annual replacement plan over years of monitoring, the accuracy of the failure rate would be improved and the running cost would therefore become more realistic, allowing the scheduling plan to become near-optimal.
The methods employed in this study can be extended in future research to investigate different rehabilitation scenarios (i.e., repair and renovation) and techniques (i.e., sleeves, CIPP) considering more components and other pipe materials.