Next Article in Journal
Modelling the Hydraulic Behaviour of Growing Media with the Explicit Finite Volume Solution
Next Article in Special Issue
Quantifying the Fecal Coliform Loads in Urban Watersheds by Hydrologic/Hydraulic Modeling: Case Study of the Beauport River Watershed in Quebec
Previous Article in Journal
Application of the Multimodel Ensemble Kalman Filter Method in Groundwater System
Previous Article in Special Issue
Evaluating the Infiltration Performance of Eight Dutch Permeable Pavements Using a New Full-Scale Infiltration Testing Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Innovative Approach for Drainage Network Sizing

1
Department of Engineering, University of Naples Parthenope, Centro Direzionale di Napoli, Isola C4, Napoli 80143, Italy
2
Department of Civil, Architectural and Environmental Engineering, University of Naples Federico II, Via Claudio n.21, Napoli 80125, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2015, 7(2), 546-567; https://doi.org/10.3390/w7020546
Submission received: 30 July 2014 / Revised: 9 December 2014 / Accepted: 27 January 2015 / Published: 6 February 2015
(This article belongs to the Special Issue Sustainable Drainage Systems)

Abstract

:
In this paper, a procedure for the optimal design of rural drainage networks is presented and demonstrated. The suggested approach, exploring the potentialities offered by heuristic methods for the solution of complex optimization problems, is based on the use of a Genetic Algorithm (GA), coupled with a steady and uniform flow hydraulic module. In particular, this work has focused: on one hand, on the problems of a technical nature posed by the correct sizing of a drainage network; on the other hand, on the possibility to use a simple but nevertheless efficient GA to reach the minimal cost solution very quickly. The suitability of the approach is tested with reference to small and large scale drainage networks, already considered in the literature.

1. Introduction

The problem of the optimal design of rural drainage channels can be approached from two distinct points of view, namely the optimal design of a single channel and the optimal design of an entire channel network. Historically, due to the lack of computers and adequate numerical techniques, the optimization of the single channel’s shape and design has been considered first, and useful analytic solutions can be found in classic hydraulic engineering texts [1]. Despite the precocious availability of these results, researchers have also considered this theme recently. Guo and Hughes [2] presented an analytical procedure for the determination of the best configuration for a trapezoidal cross section of a single channel, able to minimize both frictional resistance and construction cost, taking into account the freeboard and bank slope. Mironenko et al. [3] studied the design of channels with parabolic cross-section. Loganathan [4] presented optimal conditions for a parabolic channel cross section accounting for freeboard and limitations on the velocity and channel sizes. Froehlich [5] used the Langrange’s multiplier method to determine optimal channel cross sections, incorporating in his formulation of the optimization problem, as additional constraints, both limited flow top width and depth. Monadjemi [6] used Langrange’s multipliers method to find the best hydraulic cross section area for different channel shapes. In particular, he solved the problem of optimizing the lining costs, and found that the minimization of the wetted perimeter and the minimization of the cross section area are mathematically equivalent. Swamee et al. [7,8] proposed an approach for optimal open channel design where seepage losses were also considered. Das [9] proposed an optimization model for the design of trapezoidal channels, which considers the flooding probability; the same author [10] proposed an optimization strategy to design open channels with composite lining along the perimeter. Jain et al. [11] considered spatial variations of the velocity across a proposed composite channel cross section, and approximated the solution to this problem using a Genetic Algorithm (GA). Chahar [12] faced the design of parabolic cross section channels using a nonlinear unconstrained optimization method. More recently, Reddy and Adarsh [13] used a Genetic Algorithm (GA) as Particle Swarm Optimization (PSO) to optimally design a composite trapezoidal irrigation channel.
Of course, in practical applications it is important to consider the optimal design of an entire drainage network consisting of multiple channels. With reference to this topic, few studies about the optimization of free surface rural drainage networks are available, while interest of researchers has been focused mainly on the optimal design of drainage networks. Despite their specific characteristics, there is an obvious conceptual link between these two problems. For this reason, and due to the scarcity of contributions on the topic of rural drainage network optimization, the literature available in the field of urban drainage networks will be also considered here. While numerous works focus on the optimal layout of urban drainage networks [14,15,16,17,18,19,20,21,22,23], the majority of research results concerns the optimal channel sizing of a network whose layout is already known. In other cases, the optimization procedures were oriented to solve more general problems. For instance, Lee et al. [24] proposed a methodology for efficient rehabilitation of sewer systems; Chill and Mays [25] and Zhang et al. [26] proposed different procedures to determine the optimal locations to place various types of developments in a watershed to reduce the negative impacts of urbanization on watershed stormwater systems, and then changes in flow rates and volume from natural to developed conditions; Oxley and Mays [27] proposed an optimization model, based upon the simulated annealing method, to optimize the size and location of detention basin systems including the outlet structures subject to design constraints. An interesting review of the optimal design procedures available for sewer networks has been made by Guo et al. [28].
Generally speaking, the techniques proposed for the optimal sizing of drainage networks differ by:
-
the choice of the decision variables (longitudinal slopes, ground elevations, crown elevations, etc.);
-
the constraints used during the optimization procedure;
-
one or more Objective Functions (OF) considered within the optimization procedure;
-
the optimization algorithm used;
-
the hydraulic model used to evaluate the performances of the drainage network;
-
the model used to evaluate the discharges through the network.
Classical nonlinear optimization methods, based on gradient techniques, are not satisfactory when applied to the optimal drainage network design problem, because they have a tendency to get stuck in local optima while searching for global solutions in a non-convex discrete search space. As a result of developments in Artificial Intelligence and Operation Research, different alternative optimization techniques, such as the Evolutionary Computation approaches, have emerged during the last 30 years. With reference to the ability to achieve fast results, Wang et al. [29] made a comparison between GA [30,31,32,33,34], Particle Swarm Optimization [35] and Ant Colony Algorithm [36,37,38], showing that the Ant Colony methods require minor computational burden. Afshar et al. [39] used Cellular Automata approaches, obtaining results comparable to other methods but with higher computational efficiency. Conversely, GA allows obtaining the most accurate solution [32]: this class of algorithms is very robust in handling complex problems that display large variability and intermittency in input parameters and a large degree of nonlinearity in functional relationships [40,41].
In this paper, we propose a GA procedure aiming at the optimal design of rural drainage networks, which enables the network channels to convey the required discharges with minimum construction and maintenance costs, achieving the best compromise between the numerous technical conflicting requirements. In order to develop the main structure of the optimization procedure, the network hydraulic performance is evaluated by means of a very simple hydraulic model, based on a uniform and steady state stage discharge formula, and the a priori knowledge of discharges flowing through each link of the network. However, these assumptions can be easily relaxed, considering realistic hydraulic simulators, coupled with hydrological models able to evaluate the surface runoff to the channel network [42,43].
Besides the main objective of providing a general methodology for the optimal sizing of rural drainage network channels, additional objectives are considered in this paper, namely:
-
exploring the influence, on the optimal design of the network, of the value assigned to the invert elevation of the network ending node;
-
the analysis of the influence of the technical constraint which imposes, at each junction node of the network, that the size of the channel downstream is not smaller than that of the channels upstream;
-
exploring the influence of the mutation probability, which is a GA parameter to be tuned in order to achieve good solutions [44,45,46].
In the following sections, the problem of the optimal rural drainage network design is formulated, the assumptions made are described, and the optimization model is briefly recalled. Then, two case studies are presented and analyzed. Finally, a discussion of the results obtained is carried out, and general conclusions are drawn.

2. Methods

2.1. Problem Formulation

In practical cases, the problem of the rural drainage network design can have many competing solutions, and a criterion should be defined in order to choose a solution that is optimal. In the present case, we define the optimal network that minimizes the construction cost, and the OF is defined accordingly. The optimization process needs much input data, such as the layout of the system, the ground elevation at the network nodes, the location of the network outlet, the unit costs for construction, the shape of the cross sections, the range of variability of the decision variables, and the flow discharges through the network channels. Feasible solutions should satisfy a set of constraints, in order to take into account physical limitations, technical standards and good engineering practices.
With reference to Figure 1, the constraints that can be considered are summarized as follows:
c1:
if h is the water depth corresponding to the design discharge Q, the design filling degree is defined as δ = h/(Hexcc), where Hexc is the excavation depth and c is the ground subsidence. Overflow of the channels should be avoided: this constraint is represented by the condition δ ≤ δmax, where δ = 1 − fb/(Hexcc), and fb is a convenient freeboard. The design discharge is defined as Q = Q T 2 , where T2 = 10 ÷ 20 years is the design return period.
c2:
a maximum excavation depth Hexc,max has to be considered in order to limit the excavation costs and to avoid excessive drainage of sub-surface flow, with subsequent need for irrigation.
c3:
in order to reduce the construction costs, the erosion of non-lined channels bottom and banks should be controlled, taking into account the effects of moderate return period flows Qf. A criterion based on the definition of a threshold velocity Ver can be used to evaluate the start of erosion: if Vf is the velocity corresponding to the frequent flow discharge Qf, the constraint is expressed as VfVer. For the evaluation of Ver, the approach proposed by USDA [47] can be used, while Q f = Q T 1 is the flow corresponding to a moderate return period T1T2.
c4:
sediment deposition should be avoided during flow conditions that have a frequency higher than 3 ÷ 6 times per year. If Vvf is the velocity corresponding to the very frequent flow discharge Qvf, the constraint is expressed as VvfVdep. The limit velocity Vdep is a function of the diameter of the particles carried by flow, while Q v f ( 1 15 ÷ 1 10 ) Q T 2 .
c5:
a sufficient freeboard fcr, equal to the thickness of the crop-roots layer, has to be considered in order to protect crop even during flow conditions that have a frequency higher than 3 ÷ 6 times per year. If hvf is the water depth corresponding to Qvf, and δvf = hvf/(Hexcc) is the filling degree corresponding to Qvf, this constraint is expressed as δvf ≤ δcr, where δcr = 1 − fcr/(Hexcc).
c6:
at each node of the network, the dimensions of the channel downstream should not be smaller than those of the channels upstream [48,49]).
With reference to a network made up of Nr reaches and Nn nodes, let Ωr be the set of the Nr reaches, Ωn the set of the Nn nodes, and Ωup(j) the set of the reaches whose downstream end coincides with the upstream end of the generic reach j ϵ Ωr. For first order channels, the set Ωup(j) is empty. The problem of the optimal rural network design is formulated as the minimization of the following OF:
O F = j Ω r C j ( C S ( j ) , Z e x c ( j , u p ( j ) ) , Z e x c ( j , u p ( j ) ) )
where Cj is the construction cost of the channel j, C S ( j ) = [ C 1 ( j ) C 2 ( j ) ... C S N s p ( j ) ] is the vector of the channel’s geometric characteristics, up(j) and dw(j) are the upstream and downstream end nodes of the channel j, Zexc(j,n) is the bottom elevation at the end n of the channel j. In particular, the cost OF of the network is the sum of Cexc and Clin, where Cexc refers to the cost of excavation, waste transport and landfill, while Clin refers to the lining cost. In order to evaluate Cexc, the scheme of the trench considered in the calculations is shown in Figure 1.
Figure 1. Rural drainage networks: definition sketch of the symbols used.
Figure 1. Rural drainage networks: definition sketch of the symbols used.
Water 07 00546 g001
The OF is subject to the following constraints:
c1:
δ   ( j , n ) δ max   ( j , n )     j Ω r , n = u p ( j ) , d w ( j )
c2:
H e x c ( j , n ) H e x c , max ( j , n )     j Ω r , n = u p ( j ) , d w ( j )
c3:
V f   ( j , n ) V e r   ( j , n )     j Ω r , n = u p ( j ) , d w ( j )
c4:
V v f   ( j , n ) V d e p   ( j , n )     j Ω r , n = u p ( j ) , d w ( j )
c5:
δ v f   ( j , n ) δ c r   ( j , n )     j Ω r , n = u p ( j ) , d w ( j )
c6:
C S i   ( j ) C S i   ( k )     i = 1 , 2 , ... , N s p , j Ω r , k Ω u p ( j )
Though more general approaches and numerical models may be applied [50,51,52,53,54,55,56,57,58,59], in this work, for the sake of simplicity, in order to show the potential of the approach proposed for the optimal sizing of the drainage network, the actual hydraulic behavior of the whole network is neglected, and the performance of each channel is evaluated only by means of an appropriate state stage-discharge formula corresponding to uniform and steady state conditions. In particular, the Manning’s equation V = n M 1 R 2 / 3 i 1 / 2 is adopted, where nM is the Manning coefficient, R is the hydraulic radius, i = sin [tan−1 (s)], and s is the channel’s longitudinal slope.

2.2. The Genetic Algorithm

The Genetic Algorithm implemented by the authors has been described in Palumbo et al. [60]. For this reason, it will be only briefly depicted in this section. GAs are a class of heuristic techniques, inspired by the biological concepts of natural evolution and selection of individuals, which are used to sample the search space, in order to approximate the optimal solution. The candidate solutions of the optimization problem, called individuals, differ by their appearance (phenotype), i.e., by the value of the decision variables. The phenotype is coded as a genotype string, which is in turn formed by sub-strings, each representing the binary Gray coding of the decision variables. The individual characteristics determine the individual’s Fitness Function (FF) value, which depends both on the OF value related to the phenotype and on the degree of satisfaction of constraints.
At the beginning, an initial population of N individuals is randomly generated. The individuals are ranked in increasing order, according to their fitness, and a selection probability, which decreases with the ranking order, is assigned to each individual. Finally, the individuals are picked, according to their selection probability, and accumulated in a “mating pool”, in order to form couples of parents of the subsequent generation individuals. In this work, “exponential ranking” is used to select the individuals to be inserted in the mating pool for the subsequent steps of the GA processes. After “selection”, other operators can be introduced, namely “crossover”, “mutation”, and “elitism”. When the decision variables satisfy the problem constraints, the FF value coincides with the OF. Conversely, the FF value is calculated by adding penalization terms to the OF value when one or more constraints are not satisfied. This mechanism biases the selection in favor of those individuals that satisfy the constraints.
In this work, trapezoidal cross sections with fixed bank slope are adopted, and then the vector CS(j) degenerates to the bottom width B(j) of the channel j. The trench bottom elevation continuity is considered at the nodes of the rural drainage network:
H e x c ( j , u p ( j ) ) = H e x c ( k , d w ( k ) )     j Ω r , k Ω u p ( j )
Under these hypotheses, the phenotype of a candidate network is completely characterized by a vector containing the height of the trench H e x c n e n at the downstream end of the network and, for each reach, the slope s of the channel together with the bottom width B.
The actual form of the FF adopted is the following:
F F = O F + p f b j Ω r max { 0 , max n = u p ( j ) , d w ( j ) { δ ( j , n ) δ max ( j , n ) } } + p e x c j Ω r max { 0 , max n = u p ( j ) , d w ( j ) { H e x c ( j , n ) H e x c , max ( j , n ) } } + p e r j Ω r max { 0 , max n = u p ( j ) , d w ( j ) { V f ( j , n ) V e r ( j , n ) } } + p d e p j Ω r max { 0 , max n = u p ( j ) , d w ( j ) { V d e p ( j , n ) V v f ( j , n ) } } + p c r j Ω r max { 0 , max n = u p ( j ) , d w ( j ) { δ v f ( j , n ) δ c r ( j , n ) } } + p s z j Ω r max { 0 , max k Ω d w ( u p ( j ) ) { B ( l ) B ( k ) } }
In Equation (9), the symbols pfb, pexc, per, pdep, pcr and psz represent the unit penalties corresponding to the constraints of Equations (2)–(7), respectively.
The following GA parameters have been used during the numerical experiments: N = 300 individuals during each generation; I = 5000 generations; crossover probability cp = 1; Ne = 5 individuals preserved by the elitism operator; the values of the unit penalization coefficients pfb, pexc, pser, pdep, pcr, psz may vary from an application to another, and usually fall in the interval (106, 1015) when the relevant constraint is activated, while the value zero is used if the constraint is discarded. The mutation probability mp is variable in the range (0.01 ÷ 6.0) %.

3. Results and Discussion

The optimization procedure discussed in this work is applied to two case studies. The first application is taken from the existing literature about the drainage networks’ optimal design, and is used to test the GA adopted for the optimization. The second application is used to demonstrate the feasibility of the approach for real world applications.

3.1. Genetic Algorithm Verification: World Bank Network (1991)

In the literature, there is a general lack of case studies referring to the optimization of rural drainage networks, while many case studies are available for urban drainage networks. For this reason, the model implemented is easily adapted to solve the problem of the optimal urban drainage network, and then is applied to an urban drainage network with circular pipes taken from the literature [33,61,62].
The network layout is shown in Figure 2. The characteristics of this test case (network geometry, pipe diameters allowed, pipes costs, excavation costs) are summarized by Afshar and Zamani [62], and they are not repeated here. The following constraints are assumed: the maximum filling degree of the pipes is δmax = 0.82; the maximum excavation depth considered is Hexc,max = 4.5 m; the maximum allowed flow velocity is Vmax = 2.5 m/s; the minimum allowed flow velocity is Vmin = 0.5 m/s; the minimum soil cover depth is Hcov,min = 1.5 m. A set of 29 = 512 longitudinal slopes is considered in the range (0.01 ÷ 0.08) m/m, with a step equal to 1.36986 × 10−4 m/m. Finally, the diameters considered in the calculations are 24 = 16.
Figure 2. World Bank (1991) [61] case study. Layout.
Figure 2. World Bank (1991) [61] case study. Layout.
Water 07 00546 g002
The mutation probability mp must be intended here as the number Nbm of bits involved in the mutation process, divided by the total number Nbt of bits which constitute the genotype of the generic individual. Different analyses are performed in order to evaluate how the optimization process is influenced by the values assigned to the network ending node excavation H e x c n e n and to the mutation probability mp. Aiming at this, two sets of runs are considered:
-
Case WB-1: H e x c n e n is not a decision variable, and its value is taken equal to 2.00 m;
-
Case WB-2: H e x c n e n is left free to vary in the range (0.45 ÷ 2.00) m with step 0.05 m.
For each set of runs, the algorithm is restarted using different initial populations, in order to assess the robustness of the optimization model outcome, and considering variable values of the mutation probability mp.
The results obtained for the case WB-1 are summarized in Table 1.
Table 1. World Bank (1991) [61] case study. Optimal results for the case WB-1.
Table 1. World Bank (1991) [61] case study. Optimal results for the case WB-1.
NbmPop 1Pop 2Pop 3Pop 4Pop 5MinMaxMaxRMS
1199,381.54208,480.70221,530.28199,337.83199,288.43199,288.43221,530.28205,603.764866.32
2199,088.63199,108.37199,125.11199,097.66199,097.79199,088.63199,125.11199,103.518.68
3199,095.89199,108.37199,166.08199,118.22199,105.76199,095.89199,166.08199,118.8717.45
4199,109.00199,105.76199,108.50199,097.52199,111.85199,097.52199,111.85199,106.538.30
5199,098.83199,124.74199,128.56199,245.57199,169.35199,098.83199,245.57199,153.4136.96
6199,158.12199,213.63199,235.26199,242.11199,154.04199,154.04199,242.11199,200.6352.83
7199,324.87199,383.58199,599.30199,287.05199,247.62199,247.62199,599.30199,368.48136.85
In particular, the information reported in the generic row are as follows: the number Nbm of bits involved in the mutation process, the optimal cost obtained for different initial populations (Pop1, Pop2, …) with fixed Nbm, the minimum cost obtained (Min), the maximum cost (Max), the average cost (Ave), and the Root Mean Square error (RMS) of the costs. Note that the solutions are not penalized: the constraints are satisfied, and OF coincides with FF. The best solution is OF = 199,088.63, and it is obtained for Nbm = 2, corresponding to mp = 0.017. It is interesting to observe that the average optimal cost Ave attains its minimum value for Nbm = 2 as well, while the maximum cost Max and the root mean square error RMS of the costs are close to their minimum for Nbm = 2. This ensures that, for the present application, the most important numerical parameter is mp: a good choice of mp leads to reliable solutions.
The results obtained for the case WB-2 are summarized in Table 2. Again, no optimal solution is penalized: the best value for the objective function is OF = 199,088.63 and it is found for Nbm ranging between 2 and 4, corresponding to mp ϵ (0.013 ÷ 0.027). The functions Ave, Max and RMS attain their minimum values in the same range.
Table 2. World Bank (1991) [61] case study. Optimal results for the case WB-2.
Table 2. World Bank (1991) [61] case study. Optimal results for the case WB-2.
NbmPop 1Pop 2Pop 3Pop 4Pop 5MinMaxAveRMS
1202,802.76199,320.22199,289.02199,312.14199,299.74199,299.74202,802.76200,004.78747.88
2199,088.63199,098.47199,183.11199,088.63199,098.47199,088.63199,183.11199,111.4619.10
3199,105.76199,088.63199,088.63199,128.11199,135.48199,088.63199,128.11199,109.3212.72
4199,095.89199,088.63199,129.97199,111.93199,118.18199,088.63199,129.97199,108.9211.27
5199,136.47199,240.47199,139.26199,202.66199,089.27199,089.27199,240.47199,161.6240.45
6199,199.19199,206.59199,133.30199,123.84199,220.07199,199.19199,220.07199,176.6043.20
7199,180.98199,145.15199,170.09199,114.00199,227.68199,145.15199,227.68199,167.5839.16
8199,198.38199,201.65199,264.49199,203.22199,822.13199,198.38199,822.13199,337.97155.81
9199,260.81199,304.59199,396.04199,297.51199,267.17199,260.81199,396.04199,305.2399.26
10199,258.31199,326.11199,963.77199,972.76199,318.42199,258.31199,963.77199,567.88259.66
In Table 3, the results obtained for this set of runs are compared with those obtained by other authors.
Table 3. World Bank (1991) case study. Optimal results obtained by various researchers.
Table 3. World Bank (1991) case study. Optimal results obtained by various researchers.
ModelCost ($)
SEWER (World Bank 1991) [62]199,480
Afshar and Zamani (2002) [63]199,320
Afshar et al. (GA-TRANS2, 2006) [36]199,244
Proposed Model199,088.63
By inspection of the results listed in the Table 1, Table 2 and Table 3, it is possible to state that:
-
the best result obtained for this test case is better than those found by previous authors (Table 3);
-
for this test case, there is no difference between the best results obtained considering H e x c n e n fixed and equal to 2.00 m, or left free to vary in the range (0.45 ÷ 2.00) m;
-
the best solutions for OF are obtained for Nbm ranging in the interval (2 ÷ 4), which corresponds to mp ranging approximately in the interval (0.013 ÷ 0.027). This result is in agreement with the values of mp often suggested in the GA literature, with reference to hydraulic engineering applications [28,63];
-
the functions Ave, Max and RMS attain their minimum values in the same range of mp where OF is minimized. This fact ensures the reliability of the optimal solution found.
The characteristics of the optimal network obtained with the proposed approach are listed in Table 4. It is interesting to observe that, in the case under examination, the constraint c6 (no decreasing size of the channel in the downstream direction) is automatically satisfied and then superfluous.
Table 4. World Bank (1991) [61] case study. Optimal decision variables and hydraulic characteristics.
Table 4. World Bank (1991) [61] case study. Optimal decision variables and hydraulic characteristics.
BranchCrown Elevation (m)Diameter (mm)Slope (m/m)Velocity (m/s)Filling Degree (m/m)
UpstreamDownstream
1–31394.59631387.08841500.0722.0630.456
2–31393.89381387.08842500.0282.0570.624
3–51385.48551380.27673000.0272.3070.684
4–51376.60601374.46581500.0762.4990.739
5–301387.08841380.27673000.0302.4530.674
30–311380.27671378.31784500.0182.4960.711
31–251378.31781377.49864500.0182.4960.711
24–251377.49861374.46584500.0172.4370.727
25–261374.46581371.00005000.0162.4940.681

3.2. Case Study: Biggiero and Pianese Network (1996)

The model is applied to a case study available in the literature [64,65], which is used to demonstrate the feasibility of the approach for real world applications. The test considered is a rural drainage network consisting of 37 reaches, whose total length is 8310 m, and 38 nodes (Figure 3). The characteristics of the network are reported in Table 5. For the sake of simplicity, though without loss of generality, the value of the frequent discharge Qf has been taken equal to the value of the very frequent discharge Qvf.
Figure 3. Biggiero and Pianese (1996) [64] case study. Layout.
Figure 3. Biggiero and Pianese (1996) [64] case study. Layout.
Water 07 00546 g003
Table 5. Biggiero and Pianese (1996) [64] case study. Geometric and hydraulic characteristics of the problem.
Table 5. Biggiero and Pianese (1996) [64] case study. Geometric and hydraulic characteristics of the problem.
BranchGround Elevation (m)Horizontal Length (m)Q (m3/s)QfQvf (m3/s)
UpstreamDownstream
1–213.60413.2042000.103730.010373
2–1113.20412.2044000.199770.019977
10–1112.65412.2042500.143100.014310
11–1212.20411.6943000.445350.044535
3–1212.45411.6944000.157540.015754
4–612.81912.5341500.0956070.0095607
5–613.12912.5343500.153820.015382
6–812.53412.1602200.309890.030989
7–812.32012.1601000.0514180.0051418
8–1512.16011.8402000.410000.041000
18–1712.28512.173700.0498720.0049872
9–1712.51512.1731900.0968210.0096821
17–1612.17312.0081100.169840.016984
24–2312.40812.1381800.0799930.0079993
23–1612.13812.0082600.122760.012276
16–1512.00811.8401200.327310.032731
15–1411.84011.6451500.767480.076748
19–1411.70511.6451500.0598840.0059884
14–1311.64511.4052000.853560.085356
12–1311.69411.4051700.641890.064189
13–2211.40510.9253001.54060.15406
21–2211.86010.9255500.238690.023869
22–2510.92510.6452001.82850.18285
20–2611.44111.0412500.0952210.0095221
27–2611.52111.0413200.146600.014660
26–2511.04110.6453300.321100.032110
25–3310.64510.3702502.17740.21774
31–3211.24510.8202500.121710.012171
28–3211.06710.8201300.0932660.0093266
32–3310.82010.3703000.327670.032767
37–3611.01110.5953200.148740.014874
30–3610.79110.5951400.0625990.0062599
36–3510.59510.3911700.278800.027880
29–3510.54710.3911200.0819490.0081949
35–3410.39110.2701100.374670.037467
33–3410.37010.2701002.46750.24675
34–3810.27010.0003002.82550.28255
The cross section shape is assumed trapezoidal, with bottom width B, while the angle between the banks and the horizontal plane is α = 45°. The values allowed for B range from 0.30 to 4.00 m, and are reported in Table 6.
Table 6. Biggiero and Pianese (1996) [64] case study. Bottom width B and network ending node excavation. H e x c n e n : the values.
Table 6. Biggiero and Pianese (1996) [64] case study. Bottom width B and network ending node excavation. H e x c n e n : the values.
ID12345678910111213141516
B (m)0.300.500.801.001.502.002.503.003.504.00------
H e x c n e n (m)0.400.450.500.550.600.650.700.800.901.001.101.201.301.401.451.50
In order to evaluate the network construction cost, the waste transport and landfill are neglected, while only excavation costs are considered. In particular, the unit excavation costs are equal to 9.97 €/m3 for Hexc ≤ 2.00 m, and are equal to 10.29 €/m3 for Hexc > 2.00 m.
The parameters used for the evaluation of Equations (2), (4) and (5), corresponding to constraints c1, c3 and c5, are chosen as follows: fb = 0 m (and then δmax = 1), c = 0 m, fcr = 0.30 m. Without loss of generality, the constraints c2 and c4 about the maximum excavation and the deposition velocity, respectively, have been discarded. The limit velocity Ver is evaluated considering silt gravels, characterized by Plastic Index value PI = 16 and porosity p = 0.35, while the sediment concentration in the water flowing through the channels is assumed to be equal to 0.7%. Under these assumptions, the approach proposed in USDA [47] allows evaluation of the erosion velocity Ver as a function of the water depth hvf corresponding to the very frequent discharge Qvf, using the formula V e r = 2.44 h v f 0.19 .
Four different series of tests are performed:
-
Case BP-1A: H e x c n e n is not a decision variable, and its value is taken equal to 1.50 m, while the constraint c6 is effective;
-
Case BP-1B: H e x c n e n is not a decision variable, and its value is taken equal to 1.50 m, while the constraint c6 is discarded;
-
Case BP-2A: H e x c n e n is considered as a decision variable, and it is left free to vary in the range (0.40 ÷ 1.50), while the constraint c6 is effective;
-
Case BP-2B: H e x c n e n is considered as a decision variable, and it is left free to vary in the range (0.40 ÷ 1.50), while the constraint c6 is discarded.
In each reach, a set of 29 = 512 longitudinal slopes is considered, variable in the range (0.0001 ÷ 0.0064) m/m with step equal to 0.00001233 m/m, while the 24 values allowed for the decision variable H e x c n e n are reported in Table 6. In order to evaluate the FF in Equation (9), the unit penalization coefficients are chosen as follows: pfb = per = pcr = 109, and pexc = pdep = 0. The value used for the unit penalty coefficient psz is 109 for the cases BP-1A and BP-2A, while it is zero for the cases BP-1B and BP-2B. For each case, the algorithm is restarted from different initial populations (Pop1, Pop2, …), and considering variable mutation probability values mp.
The results obtained for the cases BP-1A and BP-1B are reported in Table 7. With reference to the case BP-1A, the best solution is OF = 98,972.09€, and it is obtained for Nbm = 5, corresponding to mp = 0.0075. For the same case, the average optimal cost Ave attains its minimum value for Nbm = 9, corresponding to mp = 0.0150, together with the maximum cost Max and the root mean square RMS of the costs. With reference to the case BP-1B, the best solution is OF = 85,539.03€, and it is obtained for Nbm = 5, corresponding to mp = 0.0075: due to the absence of the constraint about the channel width, a degree of freedom is added, and the best result obtained for the case BP-1B is not greater than the best result for BP-1A. The optimal values for Ave, Max and RMS are obtained for mp ranging in the interval (0.0075 ÷ 0.0225).
The results for the cases BP-2A and BP-2B are reported in Table 8.
Table 7. Biggiero and Pianese (1996) [64] case study. Optimal results for the cases BP-1A and BP-1B.
Table 7. Biggiero and Pianese (1996) [64] case study. Optimal results for the cases BP-1A and BP-1B.
CasempNbmPop1Pop2Pop3Pop4Pop5MinMaxAveRMS
BP-1A0.0011150,114.97173,070141,479.11145,121.82135,643.23135,643.23173,070149,085.8228,989.15
0.00755122,515.75100,911.57109,885.81108,321.0498,972.0998,972.09122,515.75108,121.2510,754.46
0.0159121,415.61108,932.41101,410.799,967.7102,357.299,967.7121,415.61106,816.7210,145.37
0.022514131,916.7147,795.94191,307.27141,376.08143,637.42131,916.7191,307.27151,206.6830,785.86
0.0319153,819.02145,978.72170,365.79209,552.47134,984.9134,984.9209,552.47162,940.1836,507.51
0.037524209,401.27152,818.34221,438.39216,116.22157,821.18152,818.34221,438.39191,519.0849,230.73
BP-1B0.0011104,690.9996,109.18132,357.17110,152.17122,785.9796,109.18132,357.17113,219.112,986.66
0.00755135,483.19100,798.791,897.1785,539.0393,541.785,539.03135,483.19101,451.9610,159.75
0.015987,205.04106,330.03101,287.47103,712.54135,992.8587,205.04135,992.85106,905.5811,343.42
0.022514106,298.8490,893.8109,810.0999,335.5196,656.0790,893.8109,810.09100,598.866710.11
0.0319109,417.35100,446.6108,456.51103,556.1899,369.4999,369.49109,417.35104,249.237837.56
0.037524126,757.81104,155.58101,745.82102,967.66104,154.92101,745.82126,757.81107,956.3610,195.95
Table 8. Biggiero and Pianese (1996) [64] case study. Optimal results for the cases BP-2A and BP-2B.
Table 8. Biggiero and Pianese (1996) [64] case study. Optimal results for the cases BP-2A and BP-2B.
CasempNbmPop1Pop2Pop3Pop4Pop5MinMaxAveRMS
BP-2A0.0011114,994.54129,701.43123,222.31121,269.94110,152.71110,152.71129,701.43119,868.1911,874.28
0.0075599,147.22111,655.6998,806.23104,279.0894,343.2294,343.22111,655.69101,646.294255.66
0.0159104,551.4299,804.07108,322.84106,931.04104,584.2699,804.07108,322.84104,838.734935.65
0.022514131,596.06142,020.01126,446.31137,005.42148,466.35126,446.31148,466.35137,106.8319,500.35
0.0319149,737.04132,467.33197,257.66157,626.5198,950.1132,467.33198,950.1167,207.7234,740.54
0.037524205,444.33204,105.65176,471.38181,740.54264,992.84176,471.38264,992.84206,550.9552,179.77
BP-2B0.0011105,278.1888,998.65100,389.65109,368.9116,342.2788,998.65116,342.27104,075.5314,338.91
0.0075584,640.1577,488.2179,381.0977,382.4590,499.75477,382.4590,499.7581,878.334431.97
0.015982,917.7373,353.3274,360.8692,763.8287,478.17173,353.3292,763.8282,174.785172.12
0.02251488,965.4486,126.2282,616.31101,479.54125,238.1682,616.31125,238.1696,885.1412,610.8
0.031985,571.2996,805.8491,319.6183,952.8297,789.0983,952.8297,789.0991,087.738322.36
0.03752499,792.81112,730.9294,426.2999,286.95111,468.2294,426.29112,730.92103,541.0413,883.82
With reference to the case BP-2A, the best solution is OF = 94,343.22€, and it is obtained for Nbm = 5, corresponding to mp = 0.0075: due to the absence of the constraint about the excavation at the network ending node of the network, a degree of freedom is added, and the optimal solution is not greater than that obtained for the case BP-1A. For the same case, Ave and RMS attain their minimum values for Nbm = 5, corresponding to mp = 0.075, while Max is minimized using mp = 0.015. With reference to the case BP-2B, the best solution is OF = 73,353.32€, and it is obtained for Nbm = 9, corresponding to mp = 0.015: as expected, the best result obtained for the case BP-2B is not greater than the best results for BP-1B and BP-2A. The optimal values for Ave, Max and RMS are obtained for mp = 0.0075.
The optimal network characteristics are reported in Table 9 for all the cases examined. From the inspection of this Table, it is clear that the optimal decision variables are strongly sensitive to the constraints applied. For instance, with reference to the network ending reach 34–38, its bottom width B lies in the range (1.00 ÷ 1.50) m, depending on the case examined. The same is true for the first order channels. For example, the bottom width B of reach 1–2 lies in the range (0.30 ÷ 0.50) m, while the slope lies in the range (0.00145 ÷ 0.00247) m/m.
By exploring the results listed in the Tables above, it is possible to draw the following observations:
-
the optimal results depend strongly on the constraints that are applied. In particular, the optimal result of the most constrained case (BP-1A) is 35% greater than that of less constrained case (BP-2B);
-
when the constraint c6 is not explicitly enforced (cases BP-1B and BP-2B), it may happen (Table 9) that the channel bottom width decreases downstream, despite the increase of the design discharge Q. This is true when the decrease of the channel width is sufficient to compensate, from an economical point of view, the increase of the channel longitudinal slope;
-
differently from the World Bank case study, there is a significant difference between the cases of H e x c n e n fixed or variable in a range. As expected, the optimal results for the cases BP-2A and BP-2B are not greater than those related to the cases BP-1A and BP-1B;
-
the best solutions for OF, Ave, Max and RMS are obtained for mp ranging in the interval (0.0075 ÷ 0.0225), and again this result is in agreement with the values of mp often suggested in the GA literature.
Comparing the best solution cost obtained, in this work, for the case BP-2A, in which the technical constraint c6 is effective, with the cost of the network considered in [64], obtained using the same unit costs and value of H e x c n e n ( H e x c n e n = 1.4 m) (see the following Table 10 and Figure 4, in which the geometric characteristics reported in [64] and the geometric characteristics obtained for the case BP-2A have been reported), it is possible to observe that the minimum cost network obtained by the proposed optimization procedure is € 94,343.22/€ 275,339.25 = 34.3% of the cost of original network, designed just to be effective from a technical point of view, but without considering the need to reduce the intervention costs. In order to show the convergence properties of the presented approach, the behavior of the fitness function for the case BP-2A has been reported in Figure 5.
Table 9. Biggiero and Pianese (1996) [64] case study. Optimal decision variables.
Table 9. Biggiero and Pianese (1996) [64] case study. Optimal decision variables.
ReachCase BP-1ACase BP-1BCase BP-2ACase BP-2B
BsBsBsBs
(m)(m/m)(m)(m/m)(m)(m/m)(m)(m/m)
1–20.50.001950.30.002470.30.001790.30.00145
2–110.80.003080.30.002370.50.002480.50.00267
10–110.30.002540.30.00180.80.001880.50.00227
11–120.80.003110.80.0031510.002620.70.00177
3–120.50.003540.30.003030.30.002570.30.00194
4–60.30.001950.30.003820.80.003491.30.00334
5–60.30.001720.30.002470.30.002150.30.00207
6–80.50.002740.30.002790.80.001160.50.00154
7–80.30.004690.30.006290.80.002760.40.00111
8–150.50.002620.30.00130.80.005910.50.00246
18–170.80.002810.30.003280.30.003650.50.00023
9–170.80.002320.30.002150.30.003790.40.00131
17–160.80.003430.30.004550.30.002570.50.00277
24–230.80.001810.80.001570.80.001210.40.0018
23–160.80.001540.80.002230.80.002490.30.00111
16–150.80.002910.80.001740.80.005150.30.00188
15–141.50.000470.30.0030310.000190.60.00228
19–140.80.002970.30.00550.30.005760.30.00278
14–131.50.002370.30.001210.001230.50.00149
12–130.80.001320.30.0031910.004560.70.0039
13–221.50.001470.80.0013910.001580.80.00203
21–220.30.002530.80.002920.30.003010.30.00274
22–251.50.003060.80.0021110.001121.50.00091
20–260.30.001580.30.00170.30.00250.30.00145
27–260.50.001510.30.001490.30.002110.40.00127
26–250.50.003760.80.003540.80.002680.50.00264
25–331.50.001550.80.0015910.001671.10.00217
31–320.30.001650.80.002260.50.001740.40.00196
28–320.80.002070.30.002760.30.001920.50.0027
32–330.80.004710.30.004210.80.004320.30.00388
37–360.50.001370.30.001410.30.002220.40.00122
30–360.30.00180.30.002230.30.003870.30.00223
36–350.50.005010.30.004720.30.001640.30.00443
29–350.80.006290.30.006230.80.004820.30.00399
35–340.80.004830.80.004620.80.006390.30.00281
33–341.50.002160.30.0025210.0022310.00137
34–381.50.000940.80.00110.001011.10.00111
H e x c n e n (m)1.51.51.41.3
Table 10. Geometric characteristics reported in Biggiero and Pianese (1996) [64] vs. geometric characteristics obtained for the case BP-2A.
Table 10. Geometric characteristics reported in Biggiero and Pianese (1996) [64] vs. geometric characteristics obtained for the case BP-2A.
ReachBiggiero&Pianese (1996)Case BP-2A
BsBs
(m)(m/m)(m)(m/m)
1–20.50.002000.30.00179
2–110.50.002500.50.00248
10–110.50.001800.80.00188
11–120.80.001701.00.00262
3–120.80.001900.30.00257
4–60.50.001900.80.00349
5–60.80.001700.30.00215
6–80.80.001700.80.00116
7–80.50.001600.80.00276
8–151.00.001600.80.00591
18–170.50.001600.30.00365
9–170.50.001800.30.00379
17–160.50.001500.30.00257
24–230.50.001500.80.00121
23–160.80.000500.80.00249
16–150.80.001400.80.00515
15–141.50.001301.00.00019
19–140.50.000400.30.00576
14–131.50.001201.00.00123
12–131.50.001701.00.00456
13–222.00.001601.00.00158
21–221.00.001700.30.00301
22–252.00.001401.00.00112
20–260.50.001600.30.0025
27–260.80.001500.30.00211
26–251.00.001200.80.00268
25–332.00.001101.00.00167
31–320.50.001700.50.00174
28–320.50.001900.30.00192
32–331.00.001500.80.00432
37–360.80.001300.30.00222
30–360.50.001400.30.00387
36–350.80.001200.30.00164
29–350.50.001300.80.00482
35–341.00.001100.80.00639
33–342.50.001001.00.00223
34–382.50.000901.00.00101
H e x c n e n (m)1.41.4
Figure 4. Biggiero and Pianese (1996) [64] vs. case BP-2A. Layout.
Figure 4. Biggiero and Pianese (1996) [64] vs. case BP-2A. Layout.
Water 07 00546 g004
Figure 5. The behavior of fitness function for the case BP-2A.
Figure 5. The behavior of fitness function for the case BP-2A.
Water 07 00546 g005

4. Conclusions

In this work, an automated tool for the optimal design of rural drainage networks is proposed and its application and effectiveness are demonstrated. The optimization procedure makes use of a GA for the choice of the channels’ geometric characteristics that minimize the construction cost, while a uniform flow stage–discharge formula is used to evaluate the hydraulic performance of the channels and the degree of satisfaction of constraints.
Two case studies are considered. The first application, taken from the literature about the optimal design of urban drainage networks, is used to demonstrate the ability of the GA to approximate the optimal solution of the drainage network problem. The second application refers to a realistic large rural drainage network. The results of this application show that:
-
the cost of the optimal rural drainage network can be very sensitive to the choice of the value to assign to the ending node excavation depth. In particular, the optimal solution obtained fixing the ending node elevation can be much more expansive than the optimal solution obtained with the ending node excavation left free to vary in a given interval. For this reason, fixing a priori the network outlet elevation should be avoided, when possible, technically valid solutions could be obtained by exploiting the possibility that the network outlet channel leaps into the receiving water body;
-
in many cases, the optimization procedure tries to find the optimal solution by increasing the channels slope and reducing the channel width; consequently, the channels’ width may decrease in the downstream direction, despite the fact that the design discharges increase downstream. Of course, the solutions with decreasing channels’ cross section in the downstream direction are not desirable, because they are inefficient when backwater effects are present during on-stationary conditions. For this reason, the constraint c6 should be always enforced in practical cases;
-
the optimal values of the mutation probability mp fall in the range (0.0075, 0.0225) for the cases examined. This result is in good agreement with the values of mp often suggested in the GA literature, with reference to hydraulic engineering applications.
The approach proposed in this work is based on the preventive knowledge of the discharges flowing through each channel of the drainage network, and on the hypotheses of steady and uniform flow conditions. These limitations, though unable to help in establishing very different minimum cost solutions (Cimorelli et al. [43]), can be removed considering a hydrologic model for the evaluation of the discharges, and using a hydraulic model (De Saint Venant Equations or their parabolic approximation) in order to evaluate the hydraulic performance of the channels.

Acknowledgments

The present work was developed with financial contributions from the Campania Region, L.R. n.5/2002-year 2008—within the Project “Methods for the evaluation of security of pressurized water supply and distribution systems towards water contamination, also intentional, to guarantee to the users, and the optimization design of water systems”, Prot. 2014.0293987, dated 29.04.2014, CUP: E66D08000060002.

Author Contributions

All authors contributed equally to this work. In particular:
Luca Cozzolino contributed to the article writing.
Luigi Cimorelli wrote the computer programs used for the computations.
Carmela Mucherino performed the computations and the analysis of results.
Carmine Covelli contributed to the article organization and with the elaboration of figures and tables.
Domenico Pianese had the basic idea of the present work and coordinated the research group.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chow, V.T. Open Channel Hydraulics; McGraw-Hill: New York, NY, USA, 1959. [Google Scholar]
  2. Guo, C.; Hughes, W.C. Optimal channel cross section with freeboard. J. Irrig. Drain. Eng. 1984, 110, 304–314. [Google Scholar] [CrossRef]
  3. Mironenko, P.A.; Willardson, L.S.; Jenab, S.A. Parabolic canal design and analysis. J. Irrig. Drain. Eng. 1984, 110, 241–246. [Google Scholar] [CrossRef]
  4. Loganathan, G.V. Optimal design of parabolic canals. J. Irrig. Drain. Eng. 1991, 117, 716–735. [Google Scholar] [CrossRef]
  5. Froehlich, D.C. Width and depth constrained best trapezoidal section. J. Irrig. Drain. Eng. 1994, 120, 828–835. [Google Scholar] [CrossRef]
  6. Monadjemi, P. General formation of best hydraulic channel section. J. Irrig. Drain. Eng. 1994, 120, 27–35. [Google Scholar] [CrossRef]
  7. Swamee, P.K.; Mishra, G.C.; Chahar, B.R. Design of minimum seepage loss canal sections. J. Irrig. Drain. Eng. 2000, 126, 28–32. [Google Scholar] [CrossRef]
  8. Swamee, P.K.; Mishra, G.C.; Chahar, B.R. Design of minimum seepage loss canal sections with drainage layer at shallow depth. J. Irrig. Drain. Eng. 2001, 127, 287–294. [Google Scholar] [CrossRef]
  9. Das, A. Optimal channel cross section with composite roughness. J. Irrig. Drain. Eng. 2000, 126, 68–72. [Google Scholar] [CrossRef]
  10. Das, A. Optimal design of channel having horizontal bottom and parabolic sides. J. Irrig. Drain. Eng. 2007, 133, 192–197. [Google Scholar] [CrossRef]
  11. Jain, A.; Bhattacharjya, R.K.; Srinivasulu, S. Optimal design of composite channels using genetic algorithm. J. Irrig. Drain. Eng. 2004, 130, 286–295. [Google Scholar] [CrossRef]
  12. Chahar, B.R. Optimal design of parabolic canal section. J. Irrig. Drain. Eng. 2005, 131, 546–554. [Google Scholar] [CrossRef]
  13. Reddy, M.J.; Adarsh, S. Chance constrained optimal design of composite channels using meta-heuristic techniques. Water Resour. Manag. 2010, 24, 2221–2235. [Google Scholar] [CrossRef]
  14. Liebman, J.C. A heuristic aid for design of sewer networks. J. Sanit. Eng. Div. 1967, 93, 81–90. [Google Scholar]
  15. Walters, G.A.; Smith, D.K. Evolutionary design algorithm for optimal layout of tree networks. Eng. Optim. 1995, 24, 261–281. [Google Scholar] [CrossRef]
  16. Mays, L.W.; Wenzel, H.G. Optimal design of multi-level branching sewer systems. Water Resour. Res. 1976, 12, 913–917. [Google Scholar] [CrossRef]
  17. Diogo, A.F.; Walters, G.A.; de Sousa, E.R.; Graveto, V.M. Three-dimensional optimization of urban drainage systems. In Computer-Aided Civil and Infrastructure Engineering; Blackwell Publishers: Malden, MA, USA; Oxford, UK, 2000. [Google Scholar]
  18. Burch, N.; Holte, R.; Müller, M.; O’Connell, D.; Schaeffer, V. Automating Layouts of Sewers in Subdivisions; IOS Press: Amsterdam, The Netherlands, 2010; pp. 655–660. [Google Scholar]
  19. Haghighi, A.; Bakhshipour, A.E. Optimization of sewer networks using an adaptive genetic algorithm. Water Resour. Manag. 2012, 26, 3441–3456. [Google Scholar] [CrossRef]
  20. Muschalla, D.; Schrφter, K.; Schütze, M. Multiobjective evolutionary algorithms in the field of urban drainage. In Proceedings of the Seventh International Conference on Hydroinformatics, Nice, France, 4–8 September 2006.
  21. Fu, G.; Khu, S.T.; Butler, D. Optimal distribution of storage tank volume to mitigate the impact of new developments on receiving water quality. ASCE J. Environ. Eng. 2010, 136, 335–342. [Google Scholar] [CrossRef]
  22. Fu, G.; Butler, D.; Khu, S.T. Multiple objective optimal control of integrated urban wastewater systems. Environ. Model. Softw. 2008, 23, 225–234. [Google Scholar] [CrossRef]
  23. Zhu, X.; Zhang, C.; Yin, J.; Zhou, H.; Jiang, Y. Optimization of water diversion based on reservoir operating rules: Analysis of the Biliu River Reservoir, China. J. Hydrol. Eng. ASCE 2014, 19, 411–421. [Google Scholar] [CrossRef]
  24. Lee, J.H.; Baek, C.W.; Kim, J.H.; Jun, H.D.; Jo, D.J. Development of a decision making support system for efficient rehabilitation of sewer systems. Water Resour. Manag. 2009, 23, 1725–1742. [Google Scholar] [CrossRef]
  25. Chill, J.; Mays, L.W. Determination of the optimal location for developments to minimize detention requirements. Water Resour. Manag. 2013, 27, 5089–5100. [Google Scholar]
  26. Zhang, G.; Guhathakurta, S.; Lee, S.; Moore, A.; Yan, L. Grid-based land-use composition and configuration optimization for watershed stormwater management. Water Resour. Manag. 2014, 28, 2867–2883. [Google Scholar] [CrossRef]
  27. Oxley, R.L.; Mays, L.W. Optimization—Simulation model for detention basin system design. Water Resour. Manag. 2014, 28, 1157–1171. [Google Scholar] [CrossRef]
  28. Guo, Y.; Walters, G.; Savic, D. Optimal design of storm sewer networks: Past, present and future. In Proceedings of the 11th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 31 August–5 September 2008.
  29. Wang, L.; Zhou, Y.; Zhao, W. Comparative study on bionic optimization algorithms for sewer optimal design. In Natural Computation, Proceedings of the Fifth International Conference on Natural Computation, Tianjin, China, 14–16 August 2009; IEEE: Piscataway, NJ, USA; pp. 24–29.
  30. Liang, L.Y.; Thompson, R.G.; Young, D.M. Designing wastewater collection systems using genetic algorithms. In PRICAI 2000 Topics in Artificial Intelligence, Proceedings of the 6th Pacific Rim International Conference on Artificial Intelligence, Melbourne, Australia, 28 August–1 September 2000; Mizoguchi, R., Slaney, J., Eds.; Springer: Berlin/Heidelberg, Germany; pp. 415–423.
  31. Weng, H.T.; Liaw, S.L. Establishing an optimization model for sewer system layout with applied genetic algorithm. J. Environ. Inform. 2005, 5, 26–35. [Google Scholar] [CrossRef]
  32. Afshar, M.H. Application of a genetic algorithm to storm sewer network optimization. Sci. Iran. 2006, 13, 234–244. [Google Scholar]
  33. Afshar, M.H.; Afshar, A.; Marino, M.A.; Darbandi, A.A.S. Hydrograph-based storm sewer design optimization by genetic algorithm. Can. J. Civil. Eng. 2006, 33, 319–325. [Google Scholar] [CrossRef]
  34. Pan, T.C.; Kao, J.J. GA-QP model to optimize sewer system design. J. Environ. Eng. 2009, 135, 17–24. [Google Scholar] [CrossRef]
  35. Izquierdo, J.; Montalvo, I.; Pérez, R.; Fuertes, V.S. Design optimization of wastewater collection networks by PS. Comput. Math. Appl. 2008, 56, 777–784. [Google Scholar] [CrossRef]
  36. Afshar, M.H. Improving the efficiency of ant algorithms using adaptive refinement: Application to storm water network design. Adv. Water Resour. 2006, 29, 1371–1382. [Google Scholar] [CrossRef]
  37. Afshar, M.H. Partially constrained ant colony optimization algorithm for the solution of constrained optimization problems: Application to storm water network design. Adv. Water Resour. 2007, 30, 954–965. [Google Scholar] [CrossRef]
  38. Afshar, M.H. A parameter free continuous ant colony optimization algorithm for the optimal design of storm sewer networks: Constrained and unconstrained approach. Adv. Eng. Softw. 2010, 41, 188–195. [Google Scholar] [CrossRef]
  39. Afshar, M.H.; Shahidi, M.; Rohani, M.; Sargolzaei, M. Application of cellular automata to sewer network optimization problems. Sci. Iran. 2011, 18, 304–312. [Google Scholar] [CrossRef]
  40. Mitchell, T.M. Machine Learning; McGraw-Hill Inc.: New York, NY, USA, 1997; pp. 249–270. [Google Scholar]
  41. Nicklow, J.; Reed, P.; Savic, D.; Dessalegne, T.; Harrel, L.; Chan-Hilton, A.; Karamouz, M.; Minsker, B.; Ostfels, A.; Singh, A.; et al. State of the art for genetic algorithms and beyound in water resources planning management. J. Water Res. Plan. Manag. 2010, 136, 412–432. [Google Scholar] [CrossRef]
  42. Cimorelli, L.; Cozzolino, L.; Covelli, C.; Mucherino, C.; Palumbo, A.; Pianese, D. Optimal design of rural drainage networks. J. Irrig. Drain. Eng. 2013, 139, 137–144. [Google Scholar]
  43. Cimorelli, L.; Cozzolino, L.; Covelli, C.; Della Morte, R.; Pianese, D. Enhancing the efficiency of the automatic design of rural drainage networks. J. Irrig. Drain. Eng. 2014, 140, 137–144. [Google Scholar] [CrossRef]
  44. Hesser, J.; Männer, R. Towards an optimal mutation probability for genetic algorithms. Lect. Notes Comput. Sci. 1991, 496, 23–32. [Google Scholar]
  45. Srinivas, M.; Patnaik, L.M. Adaptive probabilities of crossover and mutation in genetic algorithms. IEEE Trans. Syst. Man Cybern. 1994, 24, 656–667. [Google Scholar] [CrossRef]
  46. Korejo, I.; Yang, S.; Li, C. A comparative study of adaptive mutation operators for genetic algorithms design of composite channels using genetic algorithm. In Lecture Notes in Computer Science, Proceedings of the MIC 2009 VIII Metaheuristic International Conference, Hamburg, Germany, 13–16 July 2009; Voss, S., Caserta, M., Eds.; Springer: Verlag, Berlin, Germany, 2010; pp. 1–10. [Google Scholar]
  47. Design of Open Channels; Technical Release No. 25; United States Department of Agriculture Soil Conservation Service (USDA SCS): Washington, DC, USA, 1977.
  48. Walters, G.A.; Templeman, A.B. Non-optimal dynamic programming algorithms in the design of minimum cost drainage systems. Eng. Optimiz. 1979, 4, 139–148. [Google Scholar] [CrossRef]
  49. Peng, W.; Jia, R. Improved genetic algorithms for optimal design of drainage systems. In Control, Automation, Robotics and Vision Conference, Proceedings of the 8th International Conference on Control, Automation, Robotics and Vision, Kunming, China, 6–9 December 2004; IEEE: Piscataway, NJ, USA; Volume 1, pp. 227–231.
  50. Cozzolino, L.; Della Morte, R.; Covelli, C.; del Giudice, G.; Pianese, D. Numerical solution of the discontinuous bottom shallow-water equations with hydrostatic pressure distribution at the step. Adv. Water Resour. 2011, 34, 1413–1426. [Google Scholar] [CrossRef]
  51. Cozzolino, L.; Della Morte, R.; Del Giudice, G.; Palumbo, A.; Pianese, D. A well-balanced spectral volume scheme with the wetting-drying property for the shallow-water equations. J. Hydroinform. 2012, 14, 745–760. [Google Scholar] [CrossRef]
  52. Cozzolino, L.; Della Morte, R.; Cimorelli, L.; Covelli, C.; Pianese, D. A broad-crested weir boundary condition in finite volume shallow-water numerical models. Procedia Eng. 2014, 70, 353–362. [Google Scholar] [CrossRef]
  53. Cozzolino, L.; Cimorelli, L.; Covelli, C.; Della Morte, R.; Pianese, D. Boundary conditions in finite volume schemes for the solution of shallow-water equations: The non-submerged broad-crested weir. J. Hydroinform. 2014, 16, 1235–1249. [Google Scholar] [CrossRef]
  54. Cozzolino, L.; Cimorelli, L.; Covelli, C.; Della Morte, R.; Pianese, D. Novel numerical approach for 1D variable density shallow flows over uneven rigid and erodible beds. J. Hydraul. Eng. 2013, 140, 254–268. [Google Scholar] [CrossRef]
  55. Cozzolino, L.; Della Morte, R.; Covelli, C.; Pianese, D. A well-balanced spectral volume model for constituents trasport in one-dimensional-flows. In Proceedings of the 18th International Conference on Computational Methods in Water Resources (CMWR 2010), Barcelona, Spain, 21–24 June 2010; pp. 88–95.
  56. Cimorelli, L.; Cozzolino, L.; Della Morte, R.; Pianese, D. An improved numerical scheme for the approximate solution of the Parabolic Wave model. J. Hydroinform. 2013, 15, 913–925. [Google Scholar] [CrossRef]
  57. Cimorelli, L.; Cozzolino, L.; Della Morte, R.; Pianese, D. Analytical solutions of the linearized parabolic wave accounting for downstream boundary condition and uniform lateral inflows. Adv. Water Resour. 2014, 63, 57–76. [Google Scholar] [CrossRef]
  58. D’Aniello, A.; Cozzolino, L.; Cimorelli, L.; Covelli, C.; Della Morte, R.; Pianese, D. One-dimensional simulation of debris-flow inception and propagation. In Procedia Earth and Planetary Science, Proceedings of the 3rd Italian Workshop on Landslides (IWL): Hydrological Response of Slopes through Physical Experiments, Field Monitoring and Mathematical Modeling, Naples, Italy, 23–24 October 2013; Volume 9, pp. 112–121.
  59. D’Aniello, A.; Cozzolino, L.; Cimorelli, L.; Della Morte, R.; Pianese, D. A numerical model for the simulation of debris flow triggering, propagation and arrest. Nat. Hazards 2015, 75, 1403–1433. [Google Scholar] [CrossRef]
  60. Palumbo, A.; Cimorelli, L.; Covelli, C.; Cozzolino, L.; Mucherino, C.; Pianese, D. Optimal design of urban drainage networks. Civ. Eng. Environ. Syst. 2014, 31, 79–96. [Google Scholar] [CrossRef]
  61. World Bank SEWER Design Package; World Bank: New York, NY, USA, 1991.
  62. Afshar, M.H.; Zamani, H. An improved storm water network design model in spreadsheet template. Int. J. Eng. Sci. 2002, 13, 135–148. [Google Scholar]
  63. Savic, D.A.; Walters, G.A. Genetic algorithms for least-cost design of water distribution networks. J. Water Res. Plan. Manag. 1997, 123, 67–77. [Google Scholar] [CrossRef]
  64. Biggiero, V.; Pianese, D. Sul proporzionamento delle reti di bonifica. In Scritti in onore di Mario Ippolito, Proceedings of the Scritti in onore di Mario Ippolito, Napoli, Italy, 16–17 May 1996; CUEN srl: Napoli, Italy, 1996; pp. 643–656. (In Italian)[Google Scholar]
  65. Cimorelli, L.; Covelli, C.; Mucherino, C.; Palumbo, A.; Pianese, D. Optimal design of wastewater and rainstorm drainage networks. In Proceedings of the Indo-Italian Workshop on “Advances in Fluvial Hydraulics and Water Resources Development and Management”, Pune, India, 15–16 September 2011; Chandra, P., Ranganath, L.R., Eds.; Central Water & Power Research Station: Khadakwasla, Pune, India; pp. 164–204.

Share and Cite

MDPI and ACS Style

Cozzolino, L.; Cimorelli, L.; Covelli, C.; Mucherino, C.; Pianese, D. An Innovative Approach for Drainage Network Sizing. Water 2015, 7, 546-567. https://doi.org/10.3390/w7020546

AMA Style

Cozzolino L, Cimorelli L, Covelli C, Mucherino C, Pianese D. An Innovative Approach for Drainage Network Sizing. Water. 2015; 7(2):546-567. https://doi.org/10.3390/w7020546

Chicago/Turabian Style

Cozzolino, Luca, Luigi Cimorelli, Carmine Covelli, Carmela Mucherino, and Domenico Pianese. 2015. "An Innovative Approach for Drainage Network Sizing" Water 7, no. 2: 546-567. https://doi.org/10.3390/w7020546

Article Metrics

Back to TopTop