Optimal Pipe Size Design for Looped Irrigation Water Supply System Using Harmony Search: Saemangeum Project Area

Water supply systems are mainly classified into branched and looped network systems. The main difference between these two systems is that, in a branched network system, the flow within each pipe is a known value, whereas in a looped network system, the flow in each pipe is considered an unknown value. Therefore, an analysis of a looped network system is a more complex task. This study aims to develop a technique for estimating the optimal pipe diameter for a looped agricultural irrigation water supply system using a harmony search algorithm, which is an optimization technique. This study mainly serves two purposes. The first is to develop an algorithm and a program for estimating a cost-effective pipe diameter for agricultural irrigation water supply systems using optimization techniques. The second is to validate the developed program by applying the proposed optimized cost-effective pipe diameter to an actual study region (Saemangeum project area, zone 6). The results suggest that the optimal design program, which applies an optimization theory and enhances user convenience, can be effectively applied for the real systems of a looped agricultural irrigation water supply.


Introduction
Water supply systems are mainly classified into branched and looped network systems. The main difference between the two is that, in a branched network system, the flow within each pipe is a known value, whereas in a looped network system, the flow within each pipe is considered an unknown value. Therefore, an analysis of a looped network system can be a more complex endeavor.
Water supply systems form part of a larger social infrastructure of an industrial society; their objective is the effective supply of water from a water source to an area in demand. The analysis of a water supply system can be one of the more complex mathematical problems. A significant fraction of the entire set of equations consists of nonlinear equations, and a large number of these equations must be solved simultaneously.
This process requires sufficient consideration of the law of conservation of energy and a continuity equation of mass. In this regard, over the past few decades, many methods have been developed to analyze water supply systems and perform hydraulic simulations of their steady state conditions. Commercial hydraulic analysis programs such as EPANET [1] and WaterGEMS [2] have been developed to analyze the hydraulic simulations of large water supply systems, an achievement that could not even be dreamed of in past years.
The development of such models has played an important role in the design and operation of water supply systems. However, problems related to the selection of the pipe diameter for configuring low-cost water supply systems have emerged as important issues that need to be resolved.
In recent years, many optimization methods have been used for the design of low-cost water supply systems. The process for obtaining an optimal water supply system and pipe diameter is considered important because it helps in determining the final operational costs. However, because the aforementioned problems are extremely complex, they are constrained by the types of methods selected for defining the problems, as well as by analysis methods; thus far, only a minimization of the construction costs has been experimentally applied [3][4][5][6]. 2 The Scientific World Journal If the structure of a water supply system and its constraints (water pressure and velocity) are known, the optimal design of the water supply system can be expressed in terms of the selection of a pipe diameter that minimizes the total cost. The mathematical optimization methods described earlier can easily be used to find the optimal solutions in small systems within an ideal environment.
However, if existing mathematical optimization methods are applied to actual civil engineering problems, the limitations of these methods are revealed. For example, in linear programming, because all functions applied to such problems are linear, simplified assumptions lower the accuracy of the final solution. On the other hand, in dynamic programming, too many combinations have to be considered to obtain the optimal solutions, thereby requiring a considerable amount of computational effort and storage space.
In nonlinear programming, if the initial solution is not located at a good position within the solution zone, the global optimum is unobtainable, and the initial solution cannot escape from the local optima. To overcome these disadvantages, during the last 20 years, researchers have attempted to apply new approaches to optimization technologies that do not use an existing mathematical methodology.
Rather than relying completely on conventional differential derivatives, the technologies currently being developed apply a natural evolution phenomenon, that is, the principle of the survival of the fittest, and artificial imitations of this phenomena in the optimal system design; these have yielded better results than those obtained using an existing optimal mathematical design.
However, most of these studies have disadvantages in that they were applied to small benchmark problems and were not reflected in the actual plans [17]. The present study aims to develop an optimal pipe diameter estimation technique of an actual agricultural looped irrigation water supply system using an HS algorithm. This study has two main purposes. The first is to develop an economic pipe diameter estimation algorithm and program using the optimization techniques for agricultural irrigation water supply systems. The second is to validate the developed program by applying the proposed optimized economic pipe diameter to an actual target region (Saemangeum business area, zone 6).

Saemangeum Project
The Saemangeum Project is a reclamation project intended to create land from a mud flat and sea waters along the western coast of South Korea by constructing a 33.9 km long seawall. Under the Saemangeum Project, which was started on November 16, 1991, the construction of a cofferdam was completed on April 21, 2006, and the reinforcement and embankment projects were completed on April 27, 2010. The Saemangeum seawall is listed in the Guinness Book of World Records as the longest seawall on record and is 1.4 km longer than the Zuider seawall (32.5 km) of the Netherlands, which was earlier regarded to be the longest. The construction of the seawall has resulted in the reclamation of a new region with an area of 401 km 2 , of which land and a fresh water lake account for 283 km 2 and 118 km 2 , respectively.
Upon completion of the seawall construction, the South Korean Government newly formulated its Comprehensive Saemangeum Development Plan. According to this plan, the reclaimed land is to be developed as the central agricultural and economic sector of Northeast Asia, with the reclaimed land mainly divided into nine areas, as shown in Figure 1 and Table 1.
Agricultural lands account for the largest share (30.3%) among the divided areas. These agricultural lands are used for ensuring national competitiveness, producing high-valueadded agricultural products, and developing food-industry facilities through mixed environment-friendly agriculture and ecological crop cultivation. Traditionally, agricultural irrigation water supply systems have been designed as branched water supply networks, which incur less initial costs. However, such networks are disadvantageous in that they do not ensure the reliability of the water supply.
In recent years, new cultivation methods such as greenhouse crop cultivation, high-value-added crop cultivation, and perennial cultivation have been adopted. Accordingly, agricultural irrigation water supply systems also need to provide a stable and reliable water supply, as achieved by urban water supply systems.

Harmony Search
Algorithm. The HS algorithm proposed by Geem et al. [11] is an optimization technique used in pipe design. HS is a solution-finding technique that considers an optimal solution in engineering to correspond to an optimal sound in music. Generally, heuristic search methods involve the observation of natural phenomena, but the HS method is an algorithm based on the artificial phenomenon of harmony.
When sounds are produced by various sources, they together create a single harmony. Some of these created harmonies sound pleasant, whereas others sound dissonant. Eventually, the discordant harmonies disappear through practice, and among the more appropriate harmonies (local optimum), those that are aesthetically the most beautiful (global optimum) are achieved.
In other words, the HS algorithm considers an optimal solution to be an optimal harmony found through practice. The principle of the HS algorithm can be explained in detail by first comparing how music improvisation and optimization calculations correspond to each other.
Improvisation is the spontaneous creation of notes by performers without relying on sheet music (score). The ability of the performers improves the more they perform together, and ultimately a top-level harmony is created. In such an improvisation, each performer (e.g., a saxophonist, guitarist, and double bass player, as shown in Figure 2) can be referred to as a decision variable or design variable ( 1 , 2 , and 3 in Figure 2).
The musical range of each instrument (in the case of the saxophonist, e.g., one of the notes among Do, Re, and Mi, can be created) made by the corresponding performer can be referred to as the range of each variable (in the case of 1 in Figure 2, its pipe diameter may be 100, 200, or 300 mm). Moreover, when each performer plays a different note, the harmony they create (e.g., the harmony in the figure) (i.e., saxophone, Do; double bass, Mi; and guitar, Sol) corresponds to the overall solution vector obtained (the solution vector for Figure 2 is 1 = 100 mm, 2 = 300 mm, and 3 = 500 mm) by substituting the value of each variable.
Whether the harmony played at a point in time is of high quality is judged aesthetically by the performers or audience through auditory stimuli. If the harmony is very pleasant for the performers or audience, it will often be replayed in the memories. Likewise, during an optimization, whether a solution vector is good or bad can be determined by substituting the vector in an objective function; if this yields a better functional value than the existing one, the solution vector will be preserved.
Moreover, in an improvisation, as the performance is repeated, better harmonies are created, and ultimately, a high level of ability is reached; likewise, in an optimization operation, as additional iterations are carried out, better functional values are increasingly developed, and ultimately, the optimum value is obtained.
The harmony memory (HM), harmony memory considering rate (HMCR), and pitch adjustment rate (PAR) are important factors in the HS method for finding an optimal solution. First, each musical performer should have a memory space to preserve a good harmony; before starting the important process of the HS algorithm, a harmony This is called the HM, and the maximum number of harmonies that can be stored in this storage space is called the harmony memory size (HMS). Next, to produce better solutions from the harmony storage space, which is initially filled by as many random vectors as the HMS, the HS algorithm employs three types of operators.

Random Selection.
In the random selection technique, the value of a variable is randomly selected from all values of the playable note range. If is the total number of all possible variable values, one of them is randomly selected, and the probability of this technique being adopted is 1-HMCR.

Memory Consideration.
The memory consideration technique picks the value of a variable from the existing highquality notes. In other words, a single value is picked from all values possessed by a variable within the storage space.
Its probability is HMCR, and although it can have a value between 0 and 1, a value between 0.7 and 0.95 is usually used; nevertheless, the value is changeable.

Pitch Adjustment.
For a pitch adjustment, a note obtained through a memory recall technique is considered a basic note, and its pitch is trimmed by adjusting the note based on the surrounding upper and lower notes. In an actual calculation, when a single value is obtained using a memory recall technique, it is adjusted by a one-step higher or lower value. The PAR is the probability of this technique actually being applied, and it can attain a value between 0 and 1. Generally, the PAR has a value of around 0.01 to 0.3, but this can vary. Pseudocode 1 shows the pseudocode of the HS algorithm.

Objective
Function. An objective function minimizes the design cost of an irrigation system. The algorithm developed in the present study was applied to an optimization; The Scientific World Journal 5 Figure 3: A looped water supply network applied to the target zone of Saemangeum. the construction costs, pipe material costs, and maintenance costs are considered as the design costs according to the pipe diameter. Therefore, the equation for the objective function is as follows: where ( ) is cost function (construction cost) per unit length (m) for each pipe diameter, ( ) is cost function (maintenance cost) per unit length (m) for each pipe diameter, ( ) is cost function (pipe material cost) per unit length (m) for each pipe diameter, is length of the pipe (m), is pipe diameter (mm), and is total number of pipes. Hydraulic constraint equations are considered in optimization problems. Therefore, a penalty function method is introduced to convert the optimization problem subject to constraint conditions into an optimization that is free from the constraint conditions. The final objective function, which is applied using a penalty function, can be defined as follows: where ℎ is pressure head of each node (m), ℎ min is minimum pressure head (m), ℎ max is maximum pressure head (m), V is velocity of each pipe (m/s), V min is minimum pipe velocity (m/s), V max is maximum pipe velocity (m/s), , are penalty functions with regard to the pressure and pipe velocity, and is total number of nodes. The above penalty function is applied only when the pressure of each node and the velocity of the pipe exceed either the minimum or maximum value; the equation below represents the penalty function equation applied to the present model. In the target water supply system, the minimum and maximum nodal pressures were set to 10 and 35 m, respectively, and the minimum and maximum pipe velocities were set to 0.01 and 2.5 m/s, respectively: where , are penalty constants.
When running an optimization model, if the pressure head of each node and the velocity of the pipe do not satisfy the minimum and maximum values, which are the design conditions, the penalty cost is increased by assigning a significantly greater value to so that the solution will not be selected. To prepare for a case in which the pressure head and pipe velocity fall short of the design conditions by a small margin, a model that largely satisfies all of the design conditions was implemented by assigning a large value to . A trial-and-error analysis was conducted using the and values for Saemangeum, which is the target area of the present project. The results indicate that an effective optimal design is possible when and are assigned values of 10,000,000 and 100,000,000, respectively. But, detailed studies about constraint handling techniques and determination of their parameters should be tackled to improve model efficiency and reliability in future.

Target Water Supply Network.
In the present study, proposal data on the loop-type design of the six zones of Saemangeum were obtained and applied to one of the zones. A diagram of the corresponding water supply network is shown in Figure 3. The target water supply network comprises 356 pipelines, and as mentioned earlier, some of the network consists of a circuit-type water supply. The data on the cost incurred per unit of pipe length for the different diameter pipes used in this study are listed in Table 2. For optimization, 18 types of commercial pipes with 6 The Scientific World Journal  different diameters were considered. Data on the construction and pipe material costs corresponding to the different pipe diameters were obtained from the "Water Facilities Construction Cost Estimation Report" from K-water [18], which provides estimated data on the construction costs for different steel pipe diameters. The task of optimization was carried out on Intel(R) Core(TM) i5-3570 CPU at 3.4 GHz with 4 GB RAM. EPANET [1] was used as a hydraulic analysis program.

Parameter
Settings. The number of decision variables, which should be determined through optimization, is 356 because there are 356 pipelines in the target water supply network. As indicated in Table 2, 18 pipe diameters were considered for the target water supply network. Hence, the number of possible solutions that can be considered during the design period is infinite, as mentioned in Table 3.
The parameters applied in the present program for the Saemangeum target water supply network are listed in Table 4. The size of the harmony memory (HMS), the value of the HMCR parameter, and the value of PAR were set to 30, 0.97, and 0.01, respectively.
These values, which correspond to the optimum results, are adjusted; therefore, the convergence time and efficiency of the optimal solution vary. However, when there are many   the factors in evaluating the mathematical stability of an irrigation system. The minimum nodal pressure head is mostly stable at a value greater than 10 m. In the present initial design, a looped network irrigation system is implemented by installing an additional pipeline to a branched network system. In this case, the supply path up to the demand node is determined to be a branched network, that is, only a single type.
However, in the initial design, because various supply paths are possible, the head loss is slight, and a water supply is possible through the hydraulically satisfied supply paths, a system that is more hydraulically stable than a branched network system that can be implemented. Thus, because various supply paths are possible in a looped irrigation water supply system, a looped system provides a better water supply than a branched system during abnormal operating conditions such as during an irrigation path failure or closure.

Optimal Pipe Diameter Design
Results. The pipe diameter was optimized by considering the pressure and pipe velocity constraint conditions and the HS parameters, which were explained earlier in this study. The optimization results from a cost-effective pipe diameter are shown in Figure 4.
The statistical values of the nodal pressure head and pipeline velocities, which are the results of a hydraulic analysis based on cost-effective pipe diameter and the optimal cost results, are shown in Tables 5 and 6. Overall, the pressure head and pipe velocities were confirmed to be stable, and a comparison based on the hydraulic stability and economic feasibility of the initial design was conducted.
The application results indicate that the cost reduction rate of the optimal design was considerably greater (9.08%) than that of the initial design. These results were further analyzed from the viewpoint of current practices that do not employ optimization techniques; this analysis indicates that even without using any optimization technique, branched network systems that do not significantly differ from the optimal designs can be created using the current techniques.
However, in the case of a looped network system, such as the water supply network applied in this study, the differences in the results were significant; therefore, it is necessary to determine an cost-effective pipe diameter for the optimization technique based on the results obtained when employing current practices. The hydraulic analysis results indicate that the minimum pressure head (more than 10 m) was mostly satisfied, as observed in the initial design. Furthermore, the statistical values of the nodal pressure head and pipe velocity indicate that the minimum pressure head, allowable pipe velocity, and average pipe velocity all satisfy the economical pipe velocity requirements.

Differences from Other Existing Plans
In the present study, optimal design reviews of two other design plans in addition to the proposed looped network design plan were conducted. These two design plans are of a branch type and a pump type, as shown in Figures 5 and 6, respectively.
The branch-type water supply network comprises 335 pipelines, with a total length of 37.88 km. The pump-type water supply network comprises 345 pipelines; for the water supplied by the pumping of this irrigation network, the entire area encompassing the six zones was reclassified into four new areas. The total length of the pipelines is approximately 41.39 km.
To compare and evaluate the estimation results for the optimal pipe diameter of the three water supply network systems, that is, the loop type (plan 1), branch type (plan 2), 8 The Scientific World Journal  and pump type (plan 3), the cost results according to the final optimum pipe diameter and the pipe diameters of the initial plan of each of the three networks are listed in Table 7.
The results indicate that the cost of applying the optimal design was at a minimum for plan 2 and at a maximum for plan 3. This is similar to the trends found in the initial plan. However, an examination of the varying cost rate shows that the cost reduction of the optimal design for plan 2 was 4.19% less than that of the initial plan. On the other hand, the cost increased by 0.72% for plan 3, whereas in the case of plan 1, the cost reduction rate was very high (9.08%). The results for plan 1 show that the reduction rate between the optimal cost and the total length of the pipes is inversely proportional when the pressure head and velocity conditions remain constant. Moreover, a looped irrigation system has many nodes and pipes, which vary hydraulically because pipes of different diameters are used in a pipe system; this proves that it is difficult to design a looped irrigation system economically without using an optimization technique.
These results are attributed to the fact that the self-nodal pressure head of the initial version of plan 1 is relatively greater than that of the initial version of plan 2. However, from the viewpoint of current practices, which do not employ optimization techniques, branch-type systems such as plans 2 and 3, which do not differ greatly from optimal systems, can be designed by applying current techniques. In the case of a looped network system such as plan 1, the differences between the results corresponding to the initial and optimal designs were considerable. Therefore, based on the results from current practices, it is necessary to determine a costeffective pipe diameter using an optimization technique.
The results of a hydraulic analysis in which the optimal pipe diameters for plans 1, 2, and 3 were considered are shown in Table 8. The statistical values of the nodal pressure head and pipe velocity indicate that the minimum pressure head, allowable pipe velocity, and average pipe velocity for all three plans satisfy the economical pipe velocity requirements. An examination of the nodal pressure head confirms that the minimum pressure head (10 m) is mostly stable in plans 1 and 2, as is the case of the initial plan. In the case of plan 3, the minimum pressure for the initial plan was very low (0.5 m); however, the cost increases if the minimum pressure of the initial plan (0.5 m) exceeds the minimum pressure standards (10 m) during the optimal design process.
A comparison of the three optimal design types shows that plan 2 (branch type) is the most economic optimal design based only on the criterion of minimum costs. However, because plan 2 does not differ greatly from plan 1 in terms of costs, it is necessary to derive the final design results by considering the hydraulic and maintenance aspects. Plan 1 is a case in which a looped network irrigation system is implemented by installing additional pipelines to plan 2, which is a branched system.
If the pipelines supplied up to the demand node correspond to plan 2 (branch network type), the supply path is determined to be of only one type. However, in the case of plan 1, many supply paths are present; the water supply is made possible through the supply paths, which are  hydraulically satisfactory. Therefore, plan 1, which is a more hydraulically stable system than plan 2, can be implemented.
Thus, the supply of a looped irrigation water supply system during abnormal situations such as an irrigation path failure or closure is better than that of a branched irrigation water supply system because the former has various supply paths. Unlike plans 1 and 2, plan 3 was designed by reclassifying the target pipeline system into four hydraulically independent sections, and water was supplied to each section through pumping heads. By dividing the target pipeline system into four hydraulically independent sections, the fluctuations in the water quantities by each area can be more effectively and reasonably handled, and plan 3 can respond to future pipeline maintenance and expansion plans. However, the increased use of pumps can cause maintenance difficulties and an increase in maintenance costs.

Conclusions
In the present study, the HS algorithm, which is one of the latest optimization techniques, was introduced in the design of an agricultural irrigation system, and a corresponding program was developed. The developed program was applied to the actual target area (Saemangeum business area, zone 6), and the results were presented in this paper. Currently used methods have disadvantages in that the pipe diameter has to be adjusted through a hydraulic calculation of the given water supply network, and this process has to be repeated until satisfactory results are obtained. Unlike this calculation method, the model presented herein yields results that automatically meet the hydraulic conditions through the combined use of the HS algorithm and a hydraulic analysis. Hence, a comparative analysis is simple and effective. The results obtained by applying this method to an actual largescale water supply network are better than those obtained using existing mathematical algorithms even after considering the nonlinearity, which is inevitable during the analysis. The calculation results of the optimal construction costs and the pipe diameter when applying the proposed model to the actual target region (Saemangeum business area zone 6) indicate that the optimal design results obtained using HS yield much better results (9%) in terms of cost than those of the presently utilized economic pipe diameter calculation techniques. In particular, the optimization technique was found to be more necessary in the optimal design of a looped network irrigation system than for a branched network irrigation system. Furthermore, an examination of the hydrological factors of a pipeline system in which cost-effective pipe diameters were applied showed that based on the statistical values of the head and pipe velocity, the minimum pressure head, the allowable pipe velocity, and the average pipe velocity all satisfy the requirements of an economical pipe velocity. Therefore, if the benefits of the proposed model are proven through application in future systems, it will show the model to be a useful decision-making tool for designing looped network water supply systems.