1393-Article Text PDF-5095-1-10-20130303
1393-Article Text PDF-5095-1-10-20130303
1393-Article Text PDF-5095-1-10-20130303
75–90 (2002)
Abstract
The paper approaches the optimization of looped networks supplied by direct pumping from one or
more node sources, according to demand variation. Traditionally, in pipe optimization, the objective
function is always focused on the cost criteria of network components. In this study an improved
nonlinear model is developed, which has the advantage of using not only cost criteria, but also energy
consumption, consumption of scarce resources, operating expenses etc. Discharge continuity at
nodes, energy conservation in loops and energy conservation along some paths between the pump
stations and the adequate ‘critical nodes’ are considered as constraints. This problem of nonlinear
programming with equality constraints finally turns into a system of nonlinear equations to be solved
by the ‘gradient method’. The nonlinear optimization model considers head losses or discharges
through pipes as variables to be optimized in order to establish the optimal diameters of pipes and
is coupled with a hydraulic analysis. The paper compares a nonlinear optimization model to some
others, such as the classic model of average economical velocities and Moshnin optimization model.
This shows the good performance of the new model. For different analyzed networks, the saving
of electrical energy, due to diminishing pressure losses and operation costs when applying the new
model, represents about 10 . . . 30%.
Keywords: water supply, distribution, looped networks, optimal design, nonlinear model.
1. Introduction
Distribution networks are an essential part of all water supply systems. The relia-
bility of supply is much greater in the case of looped networks. Distribution system
costs within any water supply scheme may be equal to or greater than 60% of the
entire cost of the project. Also, the energy consumed in a distribution network sup-
plied by pumping may exceed 60% of the total energy consumption of the system
[16].
Attempts should be made to reduce the cost and energy consumption of the
distribution system through optimization in analysis and design.
A water distribution network that includes pumps mounted in the pipes, pres-
sure reducing valves, and check-valves can be analyzed by several common methods
such as Hardy–Cross, linear theory, and Newton–Raphson.
Traditionally, pipe diameters are chosen according to the average economical
velocities (Hardy–Cross method) [6]. This procedure is cumbersome, uneconomi-
cal, and requires trials, seldom leading to an economical and technical optimum.
76 I. SÂRBU and F. KALMÁR
Linear programming is one of the common methods that has been used to
design water distribution systems, specially, in branched pipes systems. S ÂRBU
and B ORZA [15] approach distribution networks with concentrated outflows or
uniform outflow along the length of each pipe. In this method, pipe lengths are
considered variables to be optimized.
For optimizing the design of pipe network with closed loops that is a nonlinear
problem, the bulk transport function can be used as an objective function. Strictly,
this will not be the optimum for nonlinear flow rate – cost relationships, since
economy of scale is not introduced [17].
D IXIT and R AO [7] have used a method in which only the cost of pipes is
minimized. There are other analytical and numerical models which make use of
optimization of cost criteria [1], [3], [5], [12]. Some of these methods either require
more feasible variants, or do not include the case of looped networks supplied by
more sources and having pumps mounted in the pipes. On the other hand, all of
these optimization models consider quadratic turbulence regime of water flow.
In this study it is thought that water pumping is to be performed directly in
the network according to demand variation, by means of a complex automation
of pump stations, and that the distribution network does not have any impounding
reservoir (Fig. 1). This water pumping system is used especially for large dis-
tribution networks, situated on flat ground. The solution, even if it creates some
difficulty in operation, is flexible; if, after a while, the distribution system needs to
be further developed (new higher buildings are built), the pressure or water flow in
the network can be increased by changing the pumps and, subsequently, modifying
the network.
The present paper develops a nonlinear model for optimal design of looped
networks supplied by direct pumping from one or more node sources, which has the
advantage of using not only cost criteria, but also energy consumption, consump-
tion of scarce resources, operating expenses etc. Also, this new model considers
the transitory or quadratic turbulence regime of water flow. The discharge conti-
nuity at nodes, energy conservation in loops, and energy conservation along some
paths between the pump stations and the adequate ‘critical nodes’ are considered as
constraints. The nonlinear optimization model considers head losses or discharges
through pipes as variables to be optimized in order to establish the optimal dia-
LOOPED WATER SUPPLY NETWORKS 77
meters of pipes and is coupled with a hydraulic analysis. This model can serve as
guidelines to supplement existing procedures of network design.
9.81 12
C e = We e = 730eτ k Q p h i j + Ho , (3)
η 1
(1 + βo )t − 1
Cac = Cc + C p + Cex (6)
βo (1 + βo )t
and is considered during the whole operation period (t = Tr ).
Network included energy Wc is defined by the binomial objective function of
the form (1), where a, b, α parameters have statistically corresponding determined
values [16].
Energetic consumption Wt represents the energy included in the pipes of the
network and the energy consumed in network exploitation during one year and is
expressed by:
Wt = (βo + p1 )Wc + We , (7)
where We is the pumping energy, having the expression determined from (3).
Taking into account relations (1). . .(7) and denoting:
(1 + βo )t − 1
ra = , (8)
βo (1 + βo )t
t t
ξ1 = ra p1 + ; ξ2 = ra p2 + , (9)
Tr Tr
9.81
12
ψ= f σ ξ2 + 730ra eτ k , (10)
η 1
T
NP
Fc = ξ1 (α + bDiαj )L i j +ψ Q p, j h i j + Ho , (11)
j
i j =1 j =1
where: t is the period for which the optimization criterion expressed by the objective
function is applied, having the value 1 or Tr ; N P – number of pump stations.
The general function (11) enables us to obtain a particular objective function
by particularization of the time parameter t and of the other economic and energetic
parameters, characteristic of the distribution system. For example, from t = 1,
ra = 1, e = 1, f = 0 the minimum energetic consumption criterion is obtained.
For networks supplied by pumping, the literature [1], [5], [12], [17] suggests
the use of minimum annual total expenses criterion (CAN), but choosing the optimal
diameters obtained in this way, the networks become uneconomical at some time
after construction, due to inflation.
Therefore, it is recommended the fore-mentioned criterion to be subject to
dynamization by using the criterion of total updated minimum expenses (CTA), the
former being in fact a specific case of the latter when the investment is realised
LOOPED WATER SUPPLY NETWORKS 79
within a year; the operating expenses are the same from one year to another and
the expected life-time of the distribution system is high. In particular, the use of
energetical criteria different from cost criteria is recommendable.
Thus, another way to approach the problem which has a better validity in
time and the homogenization of the objective function is network dimensioning
according to minimum energetic consumption (WT).
8 Lij
hij = λi j r Q 2i j , (12)
π g Di j
2
in which:
1
= −2 lg + 1.138, (14)
λp D
C = Re λ p + 8 λ p − 4, (15)
D
where: Re is Reynolds number; D – pipe diameter; – absolute roughness of the
pipe wall; λ p – friction factor for quadratic turbulence regime of water flow.
Eq. (12) is difficult to use in the case of pipe networks and therefore it is
convenient to write it similar to the Chézy–Manning formula:
β
h i j = Ri j Q i j , (16)
The nonlinear optimization model (MON) allows the optimal designing of looped
networks by using one of the CAN, CTA or WT optimization criteria, expressed by
the objective function (11).
If the diameter Di j is expressed in relation (17) through the discharge and
head losses: β
1 −1 1
Di j = K r Q irj h i j r L irj , (18)
and in the objective function (11) the resulting expression is replaced, we have:
T
βα
NP
α − αr αr
Fc = ξ1 a + bK r Q i j h i j L i j L i j + ψ
r
Q p, j h i j + Ho (19)
j
i j =1 j =1
• energy conservation along some paths between the pump stations and ade-
quate ‘critical nodes’:
NTj
Z S P, j − εi j (h i j − H p,i j ) − Z o, j = 0 ( j = 1, . . . , N P) (22)
i j =1
LOOPED WATER SUPPLY NETWORKS 81
in which: Q i j is the discharge through pipe i j , with the sign (+) when entering
node j and (−) when leaving it; qj – concentrated discharge at node j with the sign
(+) for node inflow and (−) for node outflow; hi j – head loss of the pipe i j ; εi j –
orientation of flow through the pipe, having the values (+1) or (−1) as the water
flow sense is the same as or opposite to the path sense of the loop m and (0) value
if i j ∈
/ m; f m – pressure head introduced by the potential elements of the loop m
[15]; M – number of independent loops (closed-loops and pseudo-loops); ZS P, j –
available piezometric head at the pump station S Pj ; H p,i j – pumping head of the
pump mounted in the pipe i j , for the discharge Qi j , approximated by parabolic
interpolation on the pump curve given by points [16]; Z o, j – piezometric head at
the critical node O j ; N T j – pipe number of a path S Pj − O j .
The optimization model (19). . .(22) represents a nonlinear programming prob-
lem, which results in a system of non-linear equations by applying the Lagrange
non-determined coefficients method. The Lagrange function can be written as:
N−N P
N
M
T
= Fc + n Qi j + q j + m εi j h i j − f m
n=1 i= j m=1 i j ∈m
i=1 i j =1
NP
NTj
+ j Z S P, j − εi j (h i j − H p,i j ) − Z o, j (23)
j =1 i j =1
T
Hi∗j = 0 (m = 1, . . . , M) (27)
i j ∈m
i j =1
in which: βα−r
−α α+r
Hi∗j = h i j r L i jr Q i j r . (28)
Eqs. (24) can be expressed similarly with the discharge continuity equations, by
giving Q ∗i j the same sign as Qi j . Eqs. (27) are similar to the energy conservation
equations in the loops, by giving Hi∗j the same sign as hi j .
Generally, the (20), (21), (22), (24) and (27) equation system allows the
determination of variables Qi j and h i j , but we must also consider the existence of
the objective function’s extreme.
Second order derivatives of function Fc by h i j and Q i j are:
∂ 2 Fc α+r βα
− α+2r α+r
= AQ r
h r
L j ,
r
(29)
∂h i j
2 r i j i j i
∂ 2 Fc βα − r βα−2r − αr α+r
= β A Q r
h L ij .
r
(30)
∂ Q 2i j r ij ij
As Q i j ≥ 0 and h i j ≥ 0, and considering that for usual values of α [16], (α +r)/r >
0, it results that ∂ 2 Fc /∂h 2i j > 0. For practical values of α and β, (βα − r)/r < 0,
so it results that ∂2 Fc /∂ Q 2i j < 0.
Consequently, in all cases the objective function Fc has a convex-concave
form for its definition range, and, therefore, has no extreme. In order to establish
an extreme, we should specify a set of variables (Qi j or h i j ). Thus, if the flow
discharges in pipes are known, the values hi j are to be determined by minimizing
the objective function Fc . If only the head losses are the given values, the variables
Q i j are to be determined by maximizing the objective function Fc .
Considering variables hi j to be unknown, pipes discharges could be calcu-
lated in a variety of ways for Eqs. (20) to be satisfied; this, however, affects the
reliability and technical and economic-energetical conditions of the system. That
is why optimization of the flow discharges in pipes must be performed according
to the minimum bulk transport criterion [13], which takes into account the network
reliability.
In this case, computation of the optimal design of looped networks must be
performed in the following stages:
• Establishment of optimal distribution for discharges through pipes, Qi j [15].
• Determination of head losses through pipes (hi j ) and piezometric heads at the
supply nodes (Z S P, j ), by solving the nonlinear equation system (21), (22),
(24) using the gradients method [11].
LOOPED WATER SUPPLY NETWORKS 83
4. Numerical Application
The looped distribution network with the topology from Fig.2 is considered. It is
made of cast iron and is supplied with a discharge of 1.22 m3 /s, provided from two
pump stations (Q p,1 = 0.806 m3 /s, Q p,2 = 0.404 m3 /s). The following data are
known: pipe length Li j , in m, elevation head Z Tj , in m, and necessary pressure
H N j = 24 m H2 O.
84
Table 1. Hydraulic characteristics of the pipes
Pipe Li j Classic model (MVE) Moshnin model (MOM) Nonlinear model (MON)
i−j [m] Qi j Di j hi j Vi j Qi j Di j hi j Vi j Qi j Di j hi j Vi j
[m3 /s] [mm] [m] [m/s] [m3 /s] [mm] [m] [m/s] [m3 /s] [mm] [m] [m/s]
0 1 2 3 4 5 6 7 8 9 10 11 12 13
2–1 480 0.01372 150 2.370 0.78 0.01488 300 0.097 0.21 0.01715 200 0.838 0.55
3–2 560 0.04476 250 1.992 0.91 0.04467 400 0.221 0.36 0.04557 250 2.062 0.93
4–3 450 0.07719 300 1.812 1.09 0.07195 400 0.460 0.57 0.07005 350 0.682 0.73
5–4 410 0.10431 350 1.345 1.08 0.08931 350 1.316 0.93 0.08692 400 0.477 0.69
6–5 470 0.12616 350 2.235 1.31 0.07299 350 1.008 0.76 0.08514 350 1.039 0.89
12–6 380 0.15178 400 1.306 1.21 0.09861 350 1.487 1.03 0.11076 400 0.707 0.88
8–7 470 0.02729 200 2.010 0.87 0.01984 300 0.169 0.28 0.01934 200 1.034 0.62
9–8 540 0.06848 300 1.722 0.97 0.04405 350 0.422 0.46 0.03863 250 1.444 0.79
85
86 I. SÂRBU and F. KALMÁR
• all the pipes of the network are functioning in a transitory turbulence regime
of water flow;
• there is a general increase of pipe diameters obtained by optimization models
(MOM, MON) with respect to MVE, because the classical model does not
take into account the minimum consumption of energy and the diversity of
economical parameters;
• in comparison with the results obtained by MVE, those obtained by opti-
mization models are more economical, a substantial reduction of specific
energy consumption for water distribution is achieved (MOM - 29.9%, MON
88 I. SÂRBU and F. KALMÁR
5. Conclusions
Notations
References
[1] A BRAMOV, N. N., Teoria i metodica rasceta sistem podaci i raspredelenia vodî, Stroizdat,
Moskva, 1972.
[2] A RSENIE , D., O formulǎ pentru calculul coeficientului de rezistenţǎ Darcy-Weisbach de utilizat
la proiectarea conductelor sub presiune, Hidrotehnica, nr. 12, 1983.
[3] A STRAHAN ŢEV, D. – A STRAHAN ŢEV, V., Optimizarea multicriteriala a reţelelor de distribuţie
a apei din centre populate, Simp. Aplicaţii ale informaticii în proiectarea şi cercetarea de
construcţii, Sibiu, 1986.
[4] B OULOS , P. F. – W OOD , D. J., An Explicit Algorithm for Calculating Operating Parameters
for Water Networks, Journal of Civil Engineering Systems, 8, 1991.
90 I. SÂRBU and F. KALMÁR
[5] C ENEDES , A. – M ELE , P., Optimal Design of Water Distribution Networks, Journal of the
Hydraulics Division, ASCE, no. HY2, 1978.
[6] C ROSS , H., Analysis of Flow in Network of Conduits or Conductors, Bulletin no. 286, Univ.
of Illinois Engrg. Experiment Station, III, 1936.
[7] D IXIT, M. – R AO , B. V., A Simple Method in Design of Water Distribution Networks, Afro-
Asian Conference on Integrated Water Management in Urban Areas, Bombay, 1987.
[8] G OFMAN , E. – RODEH , M., Loop Equation with Unknown Pipe Characteristics, Journal of
Hydraulics Division, ASCE, no. HY9, 1981.
[9] J EPPSON , R., Equivalent Hydraulic Pipe for Parallel Pipes, Journal of Hydraulics Division,
ASCE, no. HY1, 1982.
[10] O RMSBEE , L. E. – W OOD , D. J., Hydraulic Design Algorithms for Pipe Networks, Journal
of Hydraulic Engineering, ASCE, no. HY2, 1986.
[11] P CHÉNITCHNY, B. – DANILINE , Y., Méthodes numériques dans les problèmes d’extremum,
Edition Mir, Moscou, 1977.
[12] P ÎSL ǍRA ŞU , I. – ROTARU , N. – T EODORESCU , M., Alimentǎri cu apǎ, Editura Tehnicǎ,
Bucureşti, 1981.
[13] L’optimisation de la répartition des débits dans les réseaux maillés de distribution d’eau, Bulet-
inul Ştiinţific al U.P. Timişoara, Tom 36, 1991.
[14] S ÂRBU , I. – B ORZA , I., Optimal Design of Water Distribution Networks, Journal of Hydraulic
Research, no. 1, 1997.
[15] S ÂRBU , I., Energetical Optimization of Water Distribution Systems, Editura Academiei,
Bucureşti, 1997.
[16] S TEPHENSON , D., Pipe Flow Analysis, Elsevier Science Publishers B.V., 1984.
[17] T HAWAT, W., Least-Cost Design of Water Distribution Systems, Journal of the Hydraulics
Division, ASCE, no. HY9, 1973.