This paper reports the application of a multi-objective evolutionary algorithm to the design of water distribution systems. The objective functions are the minimization of the capital costs and the maximization of the hydraulic reliability of the network so as to identify the tradeoff between cost and reliability. Furthermore, a new deterministic reliability index named fuzzy reliability index is introduced utilizing fuzzy logic and the definition of hydraulic reliability. The multi-objective honey-bee mating optimization (MOHBMO) algorithm is developed to achieve the design of two well-known water distribution networks. Results show that the proposed reliability measure is able to produce appropriate hydraulic reliability in the system and the MOHBMO algorithm calculates desirable Pareto-optimal solutions.
INTRODUCTION
Many recent investigations have been conducted in many aspects of water resources systems such as reservoir operation (Fallah-Mehdipour et al. 2011, 2012a, 2013a; Aboutalebi et al. 2015; Garousi-Nejad et al. 2016a), hydrology (Orouji et al. 2013; Aboutalebi et al. 2016), project management (Bozorg-Haddad et al. 2010a; Fallah-Mehdipour et al. 2012b), cultivation rules (Bozorg-Haddad et al. 2009; Fallah-Mehdipour et al. 2013b), pumping scheduling (Bozorg-Haddad et al. 2011), hydraulic structures (Bozorg-Haddad et al. 2010b), water distribution networks (Solgi et al. 2015, 2016; Bozorg-Haddad et al. 2016), operation of aquifer systems (Bozorg-Haddad & Mariño 2011), site selection of infrastructures (Karimi-Hosseini et al. 2011), and algorithmic developments (Shokri et al. 2013; Garousi-Nejad et al. 2016b). However, only a few of these works dealt with the fuzzy multi-objective design of water distribution systems (WDSs).
The use of optimization methods to design a WDS began about 30 years ago. Before that time, WDS design was usually performed based on trial and error by experienced engineers and, as a result, because of the complexity of these problems and also their large decision space, those designs were far from optimal. Thus, development of design methods with the aid of optimization techniques was inevitable. Alperovits & Shamir (1977) applied linear programming (LP) for the optimal design of WDS. Subsequently, Quindry et al. (1981), Goulter et al. (1986), and Kessler & Shamir (1989) complemented and further developed the method introduced earlier by Alperovits & Shamir (1977). Fujiwara & Kang (1990) suggested a two-phase decomposition method based on a nonlinear programming (NLP) gradient procedure. LP and NLP consider pipe diameters as continuous parameters while they are discrete commercial diameters. Thus, a modification is needed on the final continuous solutions obtained by LP or NLP methods. Furthermore, optimal solutions found by these methods usually contain one or two pipe segments with different diameters between each pair of nodes and this split-pipe design must be altered into only one diameter (Savic & Walters 1997). These modifications may impact the final optimal solution and may introduce constraint violations (Cunha & Sousa 1999). Therefore, it is advantageous to develop methods that can overcome these problems or at least lessen them.
Savic & Walters (1997) suggested the application of evolutionary algorithms to overcome the aforementioned shortcomings in the least-cost design of WDSs. Evolutionary and meta-heuristic methods are efficient tools in solving nonlinear problems and do not involve linearization or calculation of partial derivatives. Moreover, because of the global sampling of these methods, the probability of reporting local optima as the final answer decreases significantly and also the impact of a starting solution on the final solution is negligible. Other investigators suggested the application of simulated annealing (Cunha & Sousa 1999), taboo search (Lippai et al. 1999), shuffled frog leaping algorithm (Eusuff & Lansey 2003; Seifollahi-Aghmiuni et al. 2011), ant colony optimization (Maier et al. 2003), harmony search (Geem 2006), particle swarm optimization (Suribabu & Neelakantan 2006), and honey-bee mating optimization (HBMO) (Soltanjalili et al. 2010; Ghajarnia et al. 2011), in the optimal design, the rehabilitation, and the calibration of WDSs (Bozorg-Haddad et al. 2008; Sabbaghpour et al. 2012). Those algorithms have been able to report more flexible and realistic designs for WDSs in comparison with traditional mathematical methods. Evolutionary algorithms permit the search for discrete commercial diameters easily in the optimization procedure. In addition, because of their global sampling for solutions these methods are more reliable than traditional methods in finding near-global optima of complex problems.
WDS design is beset by the trade-off between cost and reliability. Mays (1996) considered two types of failures in WDSs: mechanical (such as pipe breakage, pump failure, power outages, control valve failure, etc.) and hydraulic (such as changes in demand or in pressure head, aging of pipes, inadequate pipe sizing, insufficient pumping capacity, insufficient storage capability). Several investigators have tried to combine these two types of failures to define a measure of reliability (Todini 2000). However, there are no universally accepted definitions for risk and reliability for WDSs. Moreover, Walski et al. (1987) pointed out that a reliability measure should reflect the way in which water users are affected. There are various viewpoints on how to apply reliability in WDS design (Rowell & Barnes 1982; Kettler & Goulter 1983; Morgan & Goulter 1985; Goulter & Coals 1986; Goulter & Bouchart 1990; Ostfeld & Shamir 1996; Xu & Goulter 1999; Tabesh et al. 2004, 2009; Geem 2015).
The necessity of building cost-effective WDSs while ensuring consumers' satisfaction with supply during any operation period leads us to consider the design of a WDS as a multi-objective model with the following objectives: (1) minimization of the cost and (2) maximization of the efficiency factors. Walski (2001) lamented that in spite of numerous papers on optimal design of WDS in the past 15 years none has managed to take a leading role in the design of water distribution networks. The former author lamented the common emphasis on minimizing the cost of WDS without enough attention paid to the WDSs' efficiency, this being the main reason that these methods have not gained widespread acceptance. Focusing on the minimization of WDS costs neglects the efficiency factors, such as the reliability in delivering demand with a desired pressure. Hence, the optimal design of WDS must aim at cost effectiveness and achieving maximum system reliability. This explains the pre-eminence taken by the multi-objective design of WDS.
Several researchers implemented multi-objective optimization to the design of WDSs (i.e. Farmani et al. 2005; Fu & Kapelan 2011; Wang et al. 2014; Geem 2015). Todini (2000) stated that if the WDS design is performed based on the complete satisfaction of the hydraulic constraints it may not be capable of delivering water with sufficient pressure and discharge during failures. Therefore, the probability of delivering water with a desirable pressure increases if each node is provided with more power (energy per unit time) than required so as to have a sufficient surplus to be dissipated internally in the case of failures (Todini 2000). This surplus was considered as a reliability measure by Todini (2000) and is named resilience (Ir). The latter author considered the two objective functions of cost minimization and maximization of (Ir) for the multi-objective design of two water networks.
According to Prasad & Park (2004), by maximizing Todini's (2000) resilience index some surplus pressure can be provided in the network and this does not have any influence on the redundancy of the network loops. They indicated that by having a high Ir index in a tree-shaped network is not enough to guarantee a high level of reliability since any failure, and, particularly, a pipe failure, in a tree-shaped or cost-optimized network can have severe consequences in terms of reliability. Thus, resiliency alone is not enough for assessing the reliability of a network and a modification is needed to assess the reliability of the network loops. Prasad & Park (2004) modified the resilience index in a way to be capable of ensuring the network with reliable loops and named the improved index as network resilience (In). Neither the Ir index nor the In index consider the maximum allowable pressure in their calculations. In fact, as the Ir or In reliability indices increase, the nodal pressures of the network increase without any restriction, which may lead to uncontrollable leakage of the pipes in failure conditions.
This paper proposes a new deterministic reliability index based on fuzzy logic. In addition, the multi-objective MOHBMO algorithm is presented to perform multi-objective WDS design. The results are assessed based on the final Pareto front obtained after applying the algorithm on two well-known and benchmark networks. The effectiveness of the introduced reliability measure is also examined by assessing the final solutions from the viewpoint of nodal pressures and having reliable loops. Identification of the payoff between cost and reliability is also another aim of this study, which is obtained by applying the MOHBMO algorithm in complex, nonlinear, and discrete networks.
THE FUZZY RELIABILITY INDEX




The importance of high consumption nodes is considered in a reverse order here since the FRI index is defined to be used in an optimization model. In an optimization model intending to maximize an index, solutions with values lower than the index are restrictive. As a result, the model improves solutions with lower values of the desired index. That is, the optimization model focuses on improving the condition of high-consumption nodes over those of other junctions. The sensitivity of high-consumption nodes is increased by decreasing their membership values, which is a desirable feature of WDSs. Thus, by defining the objective function to maximize the FRI index it is possible to improve the value of critical (high-consumption) nodes. This assures that the final optimal answer cannot have low FRI values in high consumption nodes. It is noteworthy that there are other considerations in the design of WDSs. For example, there may be low-demand nodes which are critical and strategically important, e.g. supplying hospitals and other essential facilities. Thus, one must consider site-specific requirements that are implemented in the FRI equations.
THE MULTI-OBJECTIVE OPTIMIZATION MODEL
The aim of this study is to minimize the cost of WDSs and maximize the FRI. Letting Cost and Reliability denote the cost of network pipes and system's reliability, respectively, the following objective function takes into account the minimization of the network capital costs and the maximization of the FRI index.


THE MOHBMO ALGORITHM
MOHBMO is the multi-objective version of HBMO algorithm presented by Bozorg-Haddad et al. (2006). The superiority of the MOHBMO over the multi-objective genetic algorithm and multi-objective particle swarm optimization has been demonstrated in other fields of study (i.e. Niknam et al. 2011), although it has not been applied to design water distribution networks. MOHBMO is inspired by the natural life of honey-bees. Each bee hive can contain one or more queens. The queens exit the hive during the mating season and perform the mating flight by attracting drones. Successful drones which mate with the queen do not necessarily belong to the queen's hive. In other words, the queen might mate with drones belonging to other hives. In addition, two queens might not live in the same hive. Therefore, if a hive contains two queens, one of them leaves the hive and migrates to another bee hive which does not contain a queen. It is possible for the migratory queen to mate also with drones of a new hive, leading to the replacement of honey-bee genes between different bee hives. This mating ritual inspired the development of the multi-objective version of the HBMO algorithm.
In MOHBMO, at first two hives with two (or more) different queens are generated. Each hive performs a single objective optimization with one objective function and attempts to improve its queen after a predefined number of iterations. In a multi-objective problem, the final Pareto solutions found by the algorithm do not extend beyond the area defined by the solutions obtained by solving the single objective problem with each objective function. Thus, to obtain a proper multi-objective Pareto it is necessary for its two extreme points to be as near as possible to the solutions found by the single objective solutions of the problem. To ensure the aforementioned Pareto condition a warm up (WU) period is added to the MOHBMO algorithm (Fallah-Mehdipour et al. 2011). During the WU period, each hive attempts to improve the queen through the global optimum of its own objective function by performing a considerable number of iterations. This produces two different solutions; one having optimal conditions related to the first objective function and the other having optimal conditions related to the second objective function. These solutions are compared with each other. If one of them dominates the other, it is then added to the non-dominated list, otherwise both of them are added to the non-dominated list. The queens are then replaced in the hives. This means that the best solution for first objective function moves to the hive which improves the second objective function, and vice versa. Thereafter, each hive starts evolving once more through its own objective function using a new queen, while at this stage the number of iterations is much less than that in the WU period. After conducting a predefined number of iterations the final solutions are compared and moved to the non-dominated list. Each new member which is added to the non-dominated list is compared with older ones and at the end of the iterations only the non-dominated solutions remain in the list.
Flowchart of the MOHBMO algorithm. NWU1 and NWU2 are the number of iterations in the WU periods for the first and second hives, respectively.
RESULTS
MOHBMO is applied to the design of two well-known benchmark problems for the design of WDSs. The following sections present the results of these applications.
Case 1: The two-loop network


Two extreme points in the Pareto of Figure 4 (points A and B) are obtained by a single objective design of the two-loop network, using each objective function separately (WU period). It is seen in Figure 4 that the MOHBMO algorithm was able to cover the area between these two boundary points.
The main aim of the multi-objective design of a WDS is to find a solution which is optimal with respect to the single objective functions. Hence, it is necessary to select one of the solutions on the Pareto set to report it as the final result of the multi-objective design. In this regard, there are different conflict resolution methods to achieve this selection. This study implements the method of Young (1993), in which a desirability or utility function is distributed for each objective function and one of them is chosen by maximizing a mathematical equation based on the gradient of different points of the Pareto (Young 1993). By using this method, point C (Figure 4) is chosen and differentiated by a triangle in Figure 4. The minimum cost solution, maximum FRI solution, and the selected point by Young's (1993) model have been named A, B, and C, respectively, in Figure 4. Their corresponding pipe diameters, pressure head, and FRI values are reported in Tables 1 and 2. Also, the values of both objective functions are shown in these tables. By comparing these values it is seen that the reliability of solution C is 40 times that of solution A, while its cost is only 2.5 times of the cost associated with solution A. A similar comparison on points C and B indicates that solution C achieves 80% of the reliability of solution B while its cost is only a third (0.33) of that of B.
The pipe diameters for solutions A, C, and B in the two-loop network
Point . | A . | C . | B . |
---|---|---|---|
Pipe no. . | D (cm) . | D (cm) . | D (cm) . |
1 | 45.72 | 45.72 | 45.72 |
2 | 25.40 | 50.80 | 69.96 |
3 | 40.46 | 50.80 | 69.96 |
4 | 10.16 | 35.56 | 69.96 |
5 | 40.64 | 45.72 | 69.96 |
6 | 25.40 | 45.72 | 69.96 |
7 | 25.40 | 50.80 | 69.96 |
8 | 2.54 | 45.72 | 69.96 |
Cost ($) | 419,000 | 1,090,000 | 3,980,000 |
Point . | A . | C . | B . |
---|---|---|---|
Pipe no. . | D (cm) . | D (cm) . | D (cm) . |
1 | 45.72 | 45.72 | 45.72 |
2 | 25.40 | 50.80 | 69.96 |
3 | 40.46 | 50.80 | 69.96 |
4 | 10.16 | 35.56 | 69.96 |
5 | 40.64 | 45.72 | 69.96 |
6 | 25.40 | 45.72 | 69.96 |
7 | 25.40 | 50.80 | 69.96 |
8 | 2.54 | 45.72 | 69.96 |
Cost ($) | 419,000 | 1,090,000 | 3,980,000 |
The nodal pressure heads and FRI indices for solutions A, C, and B in the two-loop network
Point . | A . | C . | B . | |||
---|---|---|---|---|---|---|
Node . | FRI . | Head (m) . | FRI . | Head (m) . | FRI . | Head (m) . |
2 | 0.34 | 53.25 | 0.40 | 53.25 | 0.38 | 53.25 |
3 | 0.04 | 30.46 | 0.76 | 42.41 | 0.79 | 42.93 |
4 | 0.60 | 43.45 | 0.66 | 47.19 | 0.73 | 47.78 |
5 | 0.10 | 33.81 | 0.36 | 51.88 | 0.37 | 52.74 |
6 | 0.02 | 30.44 | 0.31 | 36.55 | 0.36 | 37.64 |
7 | 0.02 | 30.55 | 0.63 | 41.55 | 0.69 | 42.64 |
FRI obj. fun. | 0.0234 | 0.9749 | 1.2055 |
Point . | A . | C . | B . | |||
---|---|---|---|---|---|---|
Node . | FRI . | Head (m) . | FRI . | Head (m) . | FRI . | Head (m) . |
2 | 0.34 | 53.25 | 0.40 | 53.25 | 0.38 | 53.25 |
3 | 0.04 | 30.46 | 0.76 | 42.41 | 0.79 | 42.93 |
4 | 0.60 | 43.45 | 0.66 | 47.19 | 0.73 | 47.78 |
5 | 0.10 | 33.81 | 0.36 | 51.88 | 0.37 | 52.74 |
6 | 0.02 | 30.44 | 0.31 | 36.55 | 0.36 | 37.64 |
7 | 0.02 | 30.55 | 0.63 | 41.55 | 0.69 | 42.64 |
FRI obj. fun. | 0.0234 | 0.9749 | 1.2055 |
Clearly, the best point on a Pareto set is the point whose objective functions' values are close to the optimal values of each single objective function. In other words, the best design for the two-loop network is the one whose cost is equal to $419,000 and its FRI reliability index is 1.2055. Obviously, such a network does not exist. Showing this point with a rectangle in Figure 4 helps evaluating the Pareto set. It can be seen in Figure 4 that the Pareto set obtained by the MOHBMO algorithm has a proper convexity towards the rectangle. This means that some of the non-dominant solutions have close proximity to this hypothetical optimal point, and even the point selected with the method of Young (1993) fall within that area too.


Nodal pressure heads for solutions A, C, and B (meters of water) in the two-loop network.
Pipe diameters of solutions A, C, and B (inches) in the two-loop network.
Case 2: The Hanoi network


Point H was selected in this Pareto set with the method of Young (1993), and is shown by a triangle in Figure 8. Capital cost and FRI values for this solution are $7,244,479 and 0.17, respectively. Solution H is able to supply 80% of the reliability of solution G while its cost is only 0.7 of that. By a similar comparison between points H and F it is concluded that although the cost of solution H is 1.2 times higher than solution F, its reliability is 13 times higher. Furthermore, Table 3 shows the diameters of solution H, while the FRI index data for this solution is reported in Table 4. It is seen in Table 4 that Node 2 appears to have a high head. The model chooses pipes' diameters such that the pressure of the nodes immediately downstream the reservoir is relatively high in order to cope with the head losses throughout the gravity network.
Pipe diameters for solution H in the Hanoi network
Pipe . | D (cm) . | Pipe . | D (cm) . |
---|---|---|---|
1 | 101.60 | 18 | 60.96 |
2 | 101.60 | 19 | 101.60 |
3 | 101.60 | 20 | 50.80 |
4 | 101.60 | 21 | 50.80 |
5 | 101.60 | 22 | 101.60 |
6 | 101.60 | 23 | 101.60 |
7 | 101.60 | 24 | 101.60 |
8 | 101.60 | 25 | 50.80 |
9 | 76.20 | 26 | 50.80 |
10 | 101.60 | 27 | 60.96 |
11 | 101.60 | 28 | 50.80 |
12 | 60.96 | 29 | 50.80 |
13 | 40.64 | 30 | 30.48 |
14 | 30.48 | 31 | 40.64 |
15 | 30.48 | 32 | 40.64 |
16 | 60.96 | 33 | 101.60 |
17 | 60.96 | 34 | 101.60 |
Pipe . | D (cm) . | Pipe . | D (cm) . |
---|---|---|---|
1 | 101.60 | 18 | 60.96 |
2 | 101.60 | 19 | 101.60 |
3 | 101.60 | 20 | 50.80 |
4 | 101.60 | 21 | 50.80 |
5 | 101.60 | 22 | 101.60 |
6 | 101.60 | 23 | 101.60 |
7 | 101.60 | 24 | 101.60 |
8 | 101.60 | 25 | 50.80 |
9 | 76.20 | 26 | 50.80 |
10 | 101.60 | 27 | 60.96 |
11 | 101.60 | 28 | 50.80 |
12 | 60.96 | 29 | 50.80 |
13 | 40.64 | 30 | 30.48 |
14 | 30.48 | 31 | 40.64 |
15 | 30.48 | 32 | 40.64 |
16 | 60.96 | 33 | 101.60 |
17 | 60.96 | 34 | 101.60 |
Nodal pressure heads and FRI indices for solution H in the Hanoi network
Node . | Head (m) . | MemF . | FRI . | Node . | Head (m) . | MemF . | FRI . |
---|---|---|---|---|---|---|---|
2 | 97.21 | 0.0092 | 0.0087 | 18 | 52.23 | 0.5227 | 0.4874 |
3 | 62.55 | 0.0099 | 0.0089 | 19 | 61.01 | 0.0100 | 0.0090 |
4 | 58.56 | 0.1051 | 0.1045 | 20 | 52.19 | 0.5252 | 0.4097 |
5 | 53.63 | 0.4305 | 0.4149 | 21 | 43.06 | 0.8718 | 0.8311 |
6 | 48.51 | 0.7681 | 0.7294 | 22 | 42.64 | 0.8928 | 0.7050 |
7 | 47.35 | 0.8449 | 0.7877 | 23 | 46.62 | 0.8928 | 0.7050 |
8 | 46.04 | 0.9312 | 0.9055 | 24 | 45.47 | 0.9690 | 0.9292 |
9 | 45.04 | 0.9976 | 0.8499 | 25 | 44.77 | 0.9845 | 0.8134 |
10 | 42.16 | 0.8126 | 0.5670 | 26 | 42.16 | 0.8127 | 0.7760 |
11 | 41.78 | 0.7878 | 0.7680 | 27 | 42.15 | 0.8121 | 0.7306 |
12 | 41.51 | 0.7694 | 0.5982 | 28 | 44.18 | 0.9461 | 0.9323 |
13 | 37.39 | 0.4980 | 0.4745 | 29 | 42.99 | 0.8672 | 0.6813 |
14 | 39.50 | 0.6371 | 0.5403 | 30 | 42.80 | 0.8548 | 0.7344 |
15 | 39.57 | 0.6416 | 0.6326 | 31 | 42.97 | 0.8659 | 0.8613 |
16 | 42.26 | 0.8189 | 0.6718 | 32 | 44.61 | 0.9745 | 0.6546 |
17 | 45.41 | 0.9730 | 0.9308 | – | – | – | – |
Node . | Head (m) . | MemF . | FRI . | Node . | Head (m) . | MemF . | FRI . |
---|---|---|---|---|---|---|---|
2 | 97.21 | 0.0092 | 0.0087 | 18 | 52.23 | 0.5227 | 0.4874 |
3 | 62.55 | 0.0099 | 0.0089 | 19 | 61.01 | 0.0100 | 0.0090 |
4 | 58.56 | 0.1051 | 0.1045 | 20 | 52.19 | 0.5252 | 0.4097 |
5 | 53.63 | 0.4305 | 0.4149 | 21 | 43.06 | 0.8718 | 0.8311 |
6 | 48.51 | 0.7681 | 0.7294 | 22 | 42.64 | 0.8928 | 0.7050 |
7 | 47.35 | 0.8449 | 0.7877 | 23 | 46.62 | 0.8928 | 0.7050 |
8 | 46.04 | 0.9312 | 0.9055 | 24 | 45.47 | 0.9690 | 0.9292 |
9 | 45.04 | 0.9976 | 0.8499 | 25 | 44.77 | 0.9845 | 0.8134 |
10 | 42.16 | 0.8126 | 0.5670 | 26 | 42.16 | 0.8127 | 0.7760 |
11 | 41.78 | 0.7878 | 0.7680 | 27 | 42.15 | 0.8121 | 0.7306 |
12 | 41.51 | 0.7694 | 0.5982 | 28 | 44.18 | 0.9461 | 0.9323 |
13 | 37.39 | 0.4980 | 0.4745 | 29 | 42.99 | 0.8672 | 0.6813 |
14 | 39.50 | 0.6371 | 0.5403 | 30 | 42.80 | 0.8548 | 0.7344 |
15 | 39.57 | 0.6416 | 0.6326 | 31 | 42.97 | 0.8659 | 0.8613 |
16 | 42.26 | 0.8189 | 0.6718 | 32 | 44.61 | 0.9745 | 0.6546 |
17 | 45.41 | 0.9730 | 0.9308 | – | – | – | – |
CONCLUSIONS
This study presented a new deterministic fuzzy reliability index named FRI for the design of WDSs. Another goal of this study was the multi-objective design of WDS. The multi-objective version of the HBMO (MOHBMO) algorithm was developed and applied to the bi-objective design of two benchmark problems with one objective function for cost minimization and the other for reliability maximization. This paper's results demonstrated that the FRI index is able to bring the nodal pressures close to their optimal values (those defined as the average of minimum and maximum permissible values in this study), and is capable of producing redundant loops. The results indicated acceptable performance of the FRI index and MOHBMO algorithm that produced a set of optimal non-dominant solutions for decision makers.
ACKNOWLEDGEMENTS
This research was supported by the University of Tehran's Research Institute under Award No. 73148796/1/07. The authors would like to thank the University of Tehran's Research Institute for its support.