Book of Design Water System
Book of Design Water System
Book of Design Water System
ELSEVl ER
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system
or transmitted in any form or by any means, electronic, mechanical, photocopying, recording
or otherwise, without the prior written permission of the publisher, Elsevier Science Publish-
ers B.V./Physical Sciences & Engineering Division, P.O. Box 330, 1000 AH Amsterdam, The
Netherlands.
Special regulations for readers in the USA - This publication has been registered with the
Copyright Clearance Center Inc. (CCC), Salem, Massachusetts. Information can be obtained
from the CCC about conditions under which photocopies of parts of this publication may be
made in the USA. All other copyright questions, including photocopying outside of the USA,
should be referred to the publisher.
No responsibility is assumed by the publisher for any injury and/or damage to persons or
property as a matter of products liability, negligence or otherwise, or from any use or opera-
tion of any methods, products, instructions or ideas contained in the material herein.
PREFACE
CONTENTS
Geochemical source o f p o l l u t i o n 1
Effect o f e v a p o r a t i o n on c o n c e n t r a t i o n s 2
Effects o f poor q u a l i t y w a t e r 2
Scaling 3
P r e d i c t i o n of s c a l i n g a n d corrosion 3
P r e v e n t i o n of s c a l i n g 3
Calcium c a r b o n a t e s c a l i n g 5
Sulphate s c a l i n g 6
A d d i t i v e s f o r t h e p r e v e n t i o n of s c a l i n g 6
Fouling 8
Control o f f o u l i n g 9
O i l emulsion b r e a k d o w n 10
Corrosion 10
Types of corrosion 13
Corrosion p r e v e n t ion 14
P o t a b l e water s t a n d a r d s 15
Agriculture and irrigation 17
Mass Balances 20
M i x e d a n d p l u g f l o w systems 21
Systems a n a l y s i s 24
T e r m i n a l concentration in a w a t e r c i r c u i t 24
A p p l i c a t i o n to a m i n e w a t e r c i r c u i t 26
Computer s i m u l a t i o n model 31
Mathematical b a s i s of model 31
Introduction 35
B a s i c mass b a l a n c e e q u a t i o n 35
Oxygen b a l a n c e in r i v e r s 37
Coupled equations f o r DO a n d BOD 37
Analytical solution 39
C a l i b r a t i o n of a m o v i n g BOD model 40
Oxygen b a l a n c e 40
Fie1d measurements 45
S i m u l a t i o n o f H y d r a u l i c Systems 51
Two-step method 52
Demonstration o f n u m e r i c a l i n a c c u r a c y 52
I m p l i c i t f i n i t e d i f f e r e n c e schemes 55
Comments on f i n i t e d i f f e r e n c e methods 56
Numerical methods f o r t h e s o l u t i o n o f s i n g l e
differential equations 57
The E u l e r method 57
The m o d i f i e d E u l e r method 59
Runge-Kutta methods 60
M u l t i s t e p methods 61
F i n i t e elements 62
Boundaries f o r n u m e r i c a l methods 62
vi i
Introduction 64
Catchmen t d e s c r i p t i o n 64
Q u a l i t y Observations 66
Fa1lout measurement 66
Relationship between t o t a l p o l l u t a n t load a n d
r u n o f f volume 67
Chemical constituents 67
Mass b a l a n c e f o r event of 18 January 1985
on H i l l b r o w catchment 72
Mass b a l a n c e f o r event of 7 M a r c h 1983
on Montgomery P a r k catchment 73
Conclusions 77
I n t roduc t ion 99
A l t e r n a t i v e s f o r optimal reuse of waste water 99
Selection of optimum d e s a l i n a t i o n methods 101
Relevant d e s l i n a t i o n methods 103
I n d u s t r i a l wastewater treatment 1 04
Reverse osmosis 104
Membrane d e s c r i p t i o n 105
EIect r o d i a I y s i s 105
Ion exchange 105
Cost a n a l y s i s 107
C a p i t a l costs 107
I n d i r e c t c a p i t a l costs 108
Running costs 108
L a b o u r costs 108
Membrane replacement 108
Conc Ius ions 111
A p p e n d i x 8.2 136
MINOP p r o g r a m f o r o p t i m i z i n g d i s t r i b u t i o n 136
MINOP l i s t o f symbols 136
Introduction 141
Cost a n a l y s i s 146
Mathematical formu tat ion 149
Results 151
Summary a n d conclusions 153
Introduction 155
T h e m a t h e m a t i c a l model 158
Optimization method 162
Int r o d u ct i o n 166
Hydraulic analysis 167
F low measurements 167
H i g h e r income r e s i d e n t i a l 169
L o w income r e s i d e n t i a I 170
Apartment b u i l d i n g s 171
Commercial a r e a s 171
Industrial 172
Conclusions 172
Appendix 174
P r o g r a m SEWSIM 174
Effect of local p e a k s 1 74
Routing effect 175
Non-Circul a r Conduits 175
I nf low components 1 76
Data 177
Program output 177
Sample d a t a f i l e 186
AUTHOR INDEX 21 7
SUBJECT INDEX 21 9
This Page Intentionally Left Blank
1
CHAPTER 1
SCAL I NG
P r e d i c t i o n of S c a l i n g and Corrosion
Prevention of S c a l i n g
3000
2800
2600
2400
2200
zoo0
1800
1600
-E
-I
1400
1200
1000
800
600
400
200
100
50
0
0 20 40 60 80 100 120 140 160
TEMPERATURE O C
Calcium carbonate s c a l i n g
Procedure :
glven temp. OC
TDS mgll
Ca mgll
Alkallnlty
proceed
1-2-3-4-5
PHS
Ryzner Stablllty Index RSI=lpH,-pH
Calclum Carbonate S a l l n g Ilkely If LSI>O
and R S l c 6
Colrorlon Ilkely If RSI>O
Sulphate s c a l i n g
FOUL I NG
Control o f f o u l i n g
CORROS ION
Metal Electrode r e a c t i o n s E q u i I ib r i u r n p o t e n t i a l
(volts)
Cathodic
area area
a-
Iron
L
\
\
\
1 - \
\
%
\
\
\
\
Potential \
\
o f Metal - Oxidation
%
%
Eh - \
Corrosion
relative 0 - %
t o hydrogen
- Corrosion
-
-1 - Immunity due t o
- low i r o n p o t e n t i a l
Types of Corrosion
There are many ways i n which corrosion can occur in the presence of
water. Corrosion i s commonly an electro-chemical phenomenon which occurs
a t an anode when electrons flow from an anode to a cathode, leaving a
positively charged anode to react with oxygen. The cathode does not
corrode. Ways i n which the electrons m i g r a t e for corrosion to occur, are
described below ( U h l i g , 1963).
G a l v a n i c Corrosion:
When e l e c t r i c a l l y dissimilar metals are in contact in o r through an
electrolyte, a p o t e n t i a l difference i s established. The more a c t i v e ( a n o d i c )
metal corrodes, as i t i s least noble.
Pitting:
A shell of permeable magnetite o r ferric h y d r o x i d e may form over an
iron surface. Salts may concentrate under the shell and the resulting
env ironmen t becomes i n c r e a s i n g I y corrosive.
Stress Corrosion :
Metals in stress may exhibit abnormal corrosive properties in a
corrosive environment. Once a c r a c k i s formed i t r a p i d l y deteriorates due
to s a l t b u i l d - u p s i m i l a r to p i t t i n g . Chlorides a n d amonia appear to be the
chief aggressors i n t h i s type of corrosion. Welding may also induce l i n e s
of corrosion unless stress r e l i e v e d .
A c i d Corrosion :
Acids, o r even carbon d i o x i d e i n solution, can increase the hydrogen
ion concentration. T h i s r e s u l t s i n general loss of metal b y corrosion. Some
chelants, e.g. NTA (nitrilotriacetic acid) may also become c o r r o s i v e as
they concentrate.
14
B a c t e r i a l Corrosion :
B a c t e r i a can cause biochemical a c t i o n which r e s u l t s i n corrosion. This
t y p e of corrosion i s often encountered i n s u l p h u r i c c o n d i t i o n s .
E l e c t r i c a l Corrosion
E l e c t r i c c u r r e n t s , d.c. i n p a r t i c u l a r , may cause severe corrosion. I f an
anode i s formed where the c u r r e n t leaves the conductor, corrosion occurs
there.
Reagent Corrosion :
Certain scale preventing agents such as acids and chelants and
complexing agents can promote corrosion
Corrosion p r e v e n t i o n
Substance Concentration,
mg/t
Lead ( P b ) 0.05
Selenium (Se) 0.01
Silver (Ag) 0.05
Lower l i m i t Upper l i m i t
\ Lucerne\ I
I I \ \ I
C o n d u c t i v i t y of g r o u n d w a t e r (mS/m)
REFERENCES
MASS BALANCES
Examples
Flow Balance a, + a2 = a3
Mass Balance ale, + a2c2 = a3
.*.C = alcl+a2c2
3
Q1 +Q2
22
Volume S
Conc. C
QiCi =
d (SC)
Q C + -
e dt
= Q C + SdC
x f o r constant S
.*. dt =
SdC
Qi C i-QeC
e.g. a t t = 0, C = C s , and a t t = m , o r Be = - o r S = 0,
c = (ai/ae)ci
Observe that i f Q . does not equal Qe, there must be i n t e r n a l gains or
losses, e.g. due to evaporation.
The previous example could be studied numerically. Although this
r e q u i r e s specific numbers, it is often the o n l y practical way of solving
more complex problems.
Assume S = 1000 m 3 , Q . = lm’/s = Q C s = 0, Ci = 500 mg/P.
e’
Choose ~t = 100 5. The choice of ~t can affect the speed of solution,
the accuracy of r e s u l t s a n d the numerical s t a b i l i t y of the computations. It
must be determined by trial, from experience or from theoretical
considerations.
c 2 -c 1
NOW Q.C. - Q C = 5- (2.10)
I I At
.’. C2 = C
1
t
5
8 . ( C . -C
1 1 1
) = C1 + 0.1(500-C1) (2.11)
t 500-C1 xo.l
c1 c2
0 0 500 50 50
100 50 450 45 95
200 95 405 40 135
300 135 365 37 172
0
SYSTEMS ANALYSIS
F I :~ ai ~ = aD + ae (2.12)
25
c I c,
a. P l u g flow
c. D i f f u s e systea
(2.14)
#
Pure w a t e r
Remaining
concentrated
S a l i n e water
,Qi
Average
salinity I
Evaporation
Groundwater
100 -
00-
E
s,E
-
-
.-f>
.-
"
a
60
c
P
40-
d0VS
dissolved solids of the o r d e r of 100 to 300 mg/e per cycle as the water
r a n from the workings back to the s h a f t .
I t was therefore not possible to i n s e r t the complete chemical process in
equation form into the computer model of the system. Empirical
r e l a t i o n s h i p s were therefore used and these will have to be verified for
each mine and each o r e mined
(2.16)
and -
dM
dt
= Qi.Cl - Q2.C2 (2.17)
12000
10000
8000
5000
4000
2000
REFERENCES
CHAPTER 3
INTRODUCTION
- _ - KC
at
ISdtAdx
direct i o n
Decay Kc Adxdt
F i g . 3.1 Mass b a l a n c e
36
... ac
-a +
t k C + v s- a2c
E,,I-~=O (3.3)
1 ) i s r a t e of increase in concentration of p o l l u t a n t
2) i s decay r a t e
3) i s advection
4) is diffusion
5) i s source
where U=
, v'(T/P) = shear velocity
a n d k i s the von Karman constant, 0.4.
But i t i s not t h a t simple i n channels as not o n l y molecular diffusion
b u t macro turbulence, tracking, dead water, s t r a t i f i c a t i o n etc. complicate
the action, therefore one needs to c a l i b r a t e models.
E l d e r ( D e i n i g e r , 1973) suggests E = A h J(ghS) (3.7)
where h = depth a n d A = coefficient ( a v e r a g i n g 0.07).
Normally d i f f u s i o n i s n e g l i g i b l e in r i v e r s , except estuaries.
Thus one gets the Streeter-Phelps equation
-a_t - - v ac-_
ax
KC
( o m i t t i n g sources)
(3.8)
(3.9)
(3.10)
or
OXYGEN BALANCE I N R I V E R S
Dissolved Oxygen
Dernond Sag Curve
Dissolved
Oxygen
P-.L._-I n-:-A
Lriilcui ruin1
R e oxygenot ion C u r v e
Deoxygenalion Curve
Distance Downstream
E f f luen t or Time
Outtall
J,
Olcygen
0 Carbonaceous plus
aL a2L aL (3.12)
and - = E 7- v- - K,L, t S
a t ax ax L
charac t e r i s t i C
At
n - l
F i g . 3.5 Solution g r i d
dC
I f - = - K L + K (C - C ) (3.15)
dt 1 2 s
a n d oxygen d e f i c i t D = C - C (3.16)
2:
_ - - KIL - K2D (3.17)
Integrating gives
KILO -K,t -K,t -K, t
D = - (e -e 1 + Doe
Kz-Ki (3.18)
One can also e v a l u a t e K, and K2 a t t C ( D e i n i n g e r 1972 p 126).
40
OXYGEN BALANCE
(3.20)
(3.21)
18h00
. . a . A
24h00
l . . . . . I I1 .....
06h00 06h00 12h00
I . . , . . . ..... .....
18h00
I
24h00
I ,
06h00
1-1 I 1+1 x
F i g . 3.7 x - t grid
- (L 1
2At i,n + L i + l , n - Li,n-l + Li+l,n-l
-- -UAt
-
2 Ax
(Li+l,n + Li+l,n-l - L .i,n - L i , n - l 1
+ At Si - At T i (3.24)
T h i s can be r e w r i t t e n as
- - *
K1 iAt
(bi+l,n + bi+l,n-l + bi,n + bi,n-l 1
1 UAt
- ci+l,n (- + -126x
2
+ A t Pi - At R i = o (3.27)
C. + M. - N. = d. (3.28)
1,n 1,n 1,n 1,n
45
Minim i se
{; (Ui,n; + V.
1,n
+ Mi,n + N.
1,n
subject to the c o n s t r a i n t s g i v e n b y equations 3.27 a n d 3.28.
(3.29)
F I ELD MEASUREMENTS
The length of stream modelled was 6 km. It was divided into four
reaches and two sets of samples were taken a s representative, one set in
mid w i n t e r and one i n mid summer. Samples were taken e v e r y hour f o r 24
hours of each section, which was p r o b a b l y a b i t sparse. DO was measured
w i t h a p o r t a b l e meter. The samples were tested f o r 5-day and 20-day BOD,
COD a n d pH, conductivity, ammonia, nitrate, nitrite, chloride, alkalinity
and suspended sol i d s were determined. Photosynthetic oxygen release was
estimated from light and dark bottle tests, and time of passage and
dispersion were determined w i t h f l u o r e s c i n dye,
Various methods were employed to calibrate the simulation model :
linear programming was used to minimize the absolute value of the
differences between observed and s i m u l a t i o n concentrations of BOD a n d DO.
The method is described elsewhere. In order to render the equations
I inear, the theoretical concentrations were approximated by observed
values whenever p r o d u c t s of two unknowns appeared in the equations. This
may have been the result of often apparently high decay rates and
unaccounted f o r sources a l o n g some of the reaches. The methods a r e b e i n g
extended to non-l i n e a r equations (McPherson and Sharland, 1979) with
encouraging results.
REFERENCES
COD C a l i b r a t i o n BOD C a l i b r a t i o n
Value Value
Parameter jyrnbo Units ?each 1 3each 2 leach 3 teach 1 Reach Z Reach : Method of determination
Reaches 1 a n d 2 -
Tracer studies
Dispersion coeff i c i e n E 10.4 0.4 10.0 10.0 0.4 10.0 Reach 3 - C a l i b r a t i o n
Value
-- - Method of
Parameter ymbol each each each Det erm i n a t ion
1 2 3
--
D i spers ion E 10.4 0.4 10.0 Assumed same a s
coefficient for March survey
BOD source/
sink ( 1 ) S -32 175 50 Model f i t t i n g
DO source/sinb
(2) -1 3 13 -1 5 Model f i t t i n g
R1
Photosynthesis
(2)
DO 3 5 5 Bottle tests
p1
~ n e r tsource/
(1)
sink -48 20 50 Model f i t t i n g
- - -
Notes: (1) A n e g a t i v e v a l u e i n d i c a t e s a source ( p o s i t i v e b e i n g a s i n k )
(2) A p o s i t i v e v a l u e i n d i c a t e s a source ( n e g a t i v e b e i n g a s i n k )
ENn 1.7s bn
-.
I
Vn = 23.40 kruleay
DFFN l0.00kruld.y
aa*
. 0.050 1/*y
.-.
KFm 2.00 I/&Y
SINt -190.0 n p / l / d a Y
91RE 8.0 n p / l / d a Y
Po0
SINt I . 3.0 np/l/dmy
-150.0 np/l/dmy
5.30 bn
20.03 krulday
10.00 krulday
3.m I/&y
1.00 I/&Y
100.0 -/1/dmy
0.0 n p / l / d a Y
0.0 np/l/daY
m.0 np/i/ea~
--
O a X I l l * 0.203 bn
OELXIII 0.241 bn
OaXl3l 0,174 km
NO of pace 1n1ew.I~ NX
No of Ilm Inl.walsNT
- 27
123
( c l T y p i c a l Time V a r i a t i o n Output
KLlPRlVER SIMULATION - RUN
.
TIME SIMULATED OBSERVED
N0.5 m a , Kn 00 5 10 15 20 25 30 3 5 44 4 2 50 55 60 6: ?$ :7 80
.
0
0.0 57.56 3.38 55.0 3.52 ox 0.
1.04 59.69 3.45 59.0 3.51 on
2.09 57.55 3.52 56.00 4.00 ON 0.
2.49 53.24 3.65 54.00 4.02 on 0 .
.. ..
3.06 52.68 3.91 59.00 4.26 0% 0 .
5.00 52.24 4.21 59.00 4.50 a .
6.04 52.55 4.47 60.97 4.50 XO 0 .
7.07 M.55 4.56 60.02 4.57 SO
7.92 49.75 4.71 55.00 4.62 NO
8.% 49.32 4.76 61.23 4.62 NO 0 .
... . ....
10.00 47.41 4.73 57.62 2.37 LO 0 .
11.04 51.56 4.65 58.80 4.25 SO 0 .
12.05 55.93 4.48 60.95 4.13 *O 0.
12.92 61.39 4.19 60.17 3.66 SO
13.96 66.30 4.04 60.24 3.81 X
..
15.00 69.45 3.76 56.00 3.78 no
16.04 70.40 3.65 56.90 3.56 NO
17.09 69.30 3.53 56.36 3.34 no . I I
0
.. .
17.92 53.32 1.54 60.46 3.33
18.96 53.95 3.42 66.55 3.52 0 0 .
20.00 60.98 3.43 57.07 3.60 0 0 .
.
21.36 49.00 3.48 39.95 3.60 0
22.09 49.12 3.50 47.13 3.62 ON 0.
22.92 51.33 3.46 53.00 3.60 0 6
23.94 54.72 3.47 37.50 3.35 0 0
25.0 54.91 3.46 27.70 3.39 0
I
5 10 15 20
51
CHAPTER 4
F i n i t e elements
Characteristics
F i n i t e difference - Implicit - Four p o i n t 0-0
Leap f r o g
Diffusive
Backward centred
L a x - Wendroff = d i f f u s i v e / l e a p frog
Ax2/At 2 2E (4.1
or At 5 A x2/2 E (4.2)
Two-step method
a‘c
_ _ EaxZ ac kC
- v - (4.3)
at ax
c. - Cidl
use aC
Ax
=
Ax
(4.4)
then C .
i,n+l
= c.
i,n
+ EAt Ci-l,n+Ci+l,n-2Ci,n - vAt ‘i+l,n-‘i,n - kC.
1,n
Ax2 A X
(4.6)
The f i r s t and last two terms on the right hand side in the above
equation for advection and decay can be used to get the first
a p p r o x i m a t i o n to C . a n d then the d i f f u s i o n term.
I, n+l
Demonstration of n u m e r i c a l i n a c c u r a c y
i I l u s t r a t e problems a n d i n a c c u r a c i e s d u e to a n i n c o r r e c t n u m e r i c a l scheme.
Neglecting the d i f f u s i o n a n d decay term, we h a v e
‘i ,n+l = C.
1,n
- v ~ t‘i+l,n - ‘i,n
A X (4.7)
53
C i ,n
I f Ax = vAt
then u s i n g a f o r w a r d difference e x p l i c i t method
= c. - c.i , n )
‘i ,n + l 1,n - (‘i+l,n
= 1 - (0-1)
0.5
F i g . 4.3 O s c i l l a t i n g scheme
54
1 -
F i g . 4.4 Numerical d i f f u s i o n
C o n t i n u i n g so, a n o s c i l l a t i n g c u r v e occurs:
I \
F i g . 4.5 Instability
55
I m p l i c i t f i n i t e d i f f e r e n c e schemes
i-1 i i +1 X
F i g . 4.6 I m p l i c i t scheme
(4.10)
E x p l i c i t method:
F i g . 4 . 7 V a r y i n g g r i d s p a c i n g (zooming)
57
NUMER I CAL METHODS FOR THE SOLUT ION OF SINGLE D IFFERENT IAL EQUATIONS
I
Y ( t ) = y (tO)+Y ( t o ) ( t - t O ) + y
II
-
( t o ) (t-t0)2+y
I1
-
( t o )(t-t0)3+ ... (4.13)
2! 3!
L e t t i n g n represent the p r e v i o u s step at time to a n d n + l represent the
next step a t t +h, the series can be w r i t t e n as:
0
-yn I I + c y n
Yn+l=~n+hyn l+h2
Ill+
*.. (4.14)
2 6
(4.15)
w i t h i n i t i a l conditions
Y(0) = 1 (4.16)
This is a linear time variant 1st order differential equation. The
a n a l y t i c a l solution to the problem,
y = 2e-t-1 w i l l be used to compare the
t
numerical r e s u l t s o f some of the methods a n d t o i l l u s t r a t e the e r r o r a t a n y
step.
The E u l e r Method
Anal y t i c a1
solution
. c
t
The numerical solution a f t e r 5 steps i s y(0.10)=1 .lo81 whereas y=2e -t-1
g i v e s the exact a n a l y t i c a l solution as y(0.10)=1.1103. Hence the a b s o l u t e
g l o b a l e r r o r i s 0.0022, i.e. two-decimal-place accuracy. Since the global
59
The M o d i f i e d E u l e r Method
1 1
yn+l = Yn + h'n +',+I (4.24)
2
I
The E u l e r a l g o r i t h m must f i r s t be used to predict yn+l so that y
n+l
can be estimated. Applying the same example (4.15) as before a n d
s u b s t i t u t i n g y1 = x+t i n t o (4.24) gives
(4.27)
= 1.0204 (4.29)
60
y 2=
1.0204 + 0.02(1 .0204+0.02)+(1.0204+0.02(1.0204+0.02)+0.04)
2 (4.30)
= 1.0416 (4.31)
Runge-Kutta Methods
= y + I ( kl +2k2+2k3+k4) (4.32)
Yn+l n6
(4.33)
= h f ( t n + i h , yn+gkl ) (4.34)
k2
k3 = hf(t,+ih,yn+$k2) (4.35)
k4 = hf(tn+l,yn+k3) (4.36)
61
kl =h(tn+yn)
=o. 1 (0+1 = 0.10000 (4.37)
M u l t i s t e p Methods
The simple Euler, Modified Euler and Runge-Kutta methods are called
s i n g l e step methods because they use o n l y the i n f o r m a t i o n from the last
step computed. I n t h i s they have the a b i l i t y to perform the next step w i t h
a d i f f e r e n t step size a n d a r e i d e a l f o r b e g i n n i n g the s o l u t i o n where o n l y
the initial conditions are available. The principle behind a multistep
method is to utilize the past values of y and/or yl to construct a
polynomial that approximates the derivative function and to extrapolate
this into the next time interval. Most multistep methods have the
d i s a d v a n t a g e that they use a constant step size h to make the c o n s t r u c t i o n
of the polynomial easier. Another d i s a d v a n t a g e of m u l t i s t e p methods i s t h a t
62
F I N I TE ELEMENTS
F i g . 4.9 F i n i t e elements
One sets up equations giving balance for each element and solve
simu I taneousl y .
Boundaries for numerical methods
REFERENCES
CHAPTER 5
I NTRODUCT ION
F i g . 5.2 H i l l b r o w Catchment
66
Q U A L I T Y OBSERVATIONS
Fa1lout measurement
- ~~
I l l OM conduct1V l t Y TDI
taman
ma/. .PI 1
--
2Ohl 1 6.20 14.31 138 1010 0.2 12 5.1 36
20hl4 6.30 13.12 112 242 0.3 10 4.1 31
20118 6.05 11.81 134 160 0.8 10 3.0 36
2ohz3 6.00 9.69 102 770 4.1 10 3.0 24
.?Oh26 5.55 9.91 100 512 8.6 11 5.1 10
2013 1 5.85 10.88 126 232 5.3 13 5.1 24
20150 5.45 13.37 120 110 15.0 16 1.1 10
2lho1 5.90 15.43 170 102 12.9 21 8.2 27
I
*/A 5.55 6.60 18 2.7 4 3.8 6
- 63
TABLE 5.4 R e s u l t s o f c h e m i c a l a n a l y s e s on d r y w e a t h e r f l o w s a m p l e s
from Montgomery P a r k
TOI 8usp.na.a
soilas
-11 -1 I
mar
PH
(mpnl
18 10
18 160- 14 6.5
14 140- 12 0.0
12 120- 10 5.5
10 100- 8 5.0
21
Fig. 5.3 Plot of p o l l u t a n t c o n c e n t r a t i o n VS. time f o r r a i n f a l l - r u n o f f
event on H i l l b r o w on 03/01/85
45 - \
-
----Suop.
HW+I+WI
-.-
Nllrole
Sulohale
TO9
Solids
-
--.-.
-
Conducllvlly
pn
cc CMorldoo
Flowale -.18
30- 5.e
25 - 5.t
20 -
15-
10-
5-
14h25 30 35 40 45 50 55 15hoo 5 10 15 T h
,/-•
r--- L*'
--.----.
-
14hOO 16h00
TIME
lEhoO
I
1,20
- /--'
E 0,90
"
-
E
W
I-
U
d 0,60
s -
s
LL
0,30 -
-
wation lainfall U l p h t of U i p h t 01 Ratio O f ?Oi I U t l O l
nlllb-
flow for the Hillbrow catchment and about 5.3 times that due to dry
weather flow f o r the Montgomery P a r k catchment.
The pollutant loading rates derived from the different sources are
summarized i n Table 5.7.
CONCLUSIONS
REFERENCES
CHAPTER 6
OPT I MUM ALLOCAT ION OF WATER RESOURCES SUBJECT TO QUAL I TY CONSTRA I NTS
I NTROOUCT ION
THE SYSTEM
Demand
8
aAW + aBW + aCW=
Q~~ + aBX + aCX = 12
aAy + aBy + aCY = 16
The q u a l i t y c o n s t r a i n t s may be w r i t t e n :
6QAW + l l Q B W + 8QCW 1 1 0 x 8
6QAx + llQBx + 8QCX 511 x 12,
6QAy + llQBy + 8QBx 5 8 x 16
SOLUTION METHOD
Y
83
TABLE 6. 3 T r a n s p o r t a t i o n m a t r i x optimized
Dunand :
D iscussionl
All water users do not require the same high quality water. Where
poorer quality is tolerable, allocation of alternative sources may be
considered. O v e r a l l economy of p u r i f i c a t i o n a n d d i s t r i b u t i o n r e s u l t .
Qua1i ty Quantity
Sources :-
Water Board 500 mg/P 100.0 Me/d
Groundwater 600 mg/e 11.5 Me/d
Wastewater 1750 mg/l Variable
Desalinated wastewater 175 mg/P Variable
Demands:-
Transfer 1750 mg/e 7.0 Me/d
System 1 700 mg/e 9.5 Me/d
System 2 700 mg/e 0.7 Me/d
System 3 700 mg/e 0.5 Me/d
Waste ( S l a c k ) 1750 mg/P unused
U + 1.45 D = 10 (6.10)
where U = Used MSW i n Me/d
D = Desalinated Wastewater in MP/d
87
Two new variables (U and D) and a new constraint Eq. 6.10 are
i n t r o d u c e d to c a t e r f o r t h e d e s a l i n a t i o n of a v a r i a b l e p e r c e n t a g e o f waste.
U was assumed t o b e 7.5 Mk‘/d, h e n c e f r o m Eq. 6.10, D w a s 1.73 Me/d. The
m a g n i t u d e o f t h e s l a c k a l l o c a t i o n t o w a s t e (W) w i l l c o n s e q u e n t l y v a r y . The
v a r i a t i o n i s a c c o r d i n g to E q . 6.11 where the sources a r e b a l a n c e d a g a i n s t
the demands:
W + 7.0 + 9.5 + 0.7 + 0.5 = 100 + 11.5 + U + D
but D = (10 - U ) 0.69 f r o m Eq. 6.10 (6.11)
hence W = 93.8 + U + ( 1 0 - U ) 0.69
W = 100.7 + 0.31 U (6.12)
From a T r a n s p o r t a t i o n E x t e n d e d a n a l y s i s U w a s 7.5 Me/d, a n d t h e 100
Mt/d was o n l y 1 MP/d, hence W was 4.03 MP/d. However, U varies now
with a maximum value of 10 Me/d, (that is without desalination).
T h e r e f o r e from Eq. 6.12
W L 103.8 Me/d (6.13)
The a c c e p t a b l e q u a l i t y assumed f o r t h e System 1 (S), System 2 ( V ) a n d
System 3 ( M ) i s s t i l l 700 m g / Q . T h e n e x t s e c t i o n r e v i e w s a n a n a l y s i s of a
r a n g e of acceptable q u a l i t i e s .
The only unrealistic assumption necessary in this analysis, using
linear programming, i s that the total costs are linearly r e l a t e d to feed
flow. The cost coefficients are summarized below in Table 6.4.
Ab b r e via tio n s f o r the sources a n d demands a r e a l s o i n d i c a te d .
DEMANDS
WB
1
R
w l O / S
0.0 0.0 24.0
I
32.5
v l M 39.5
-_
GROUND
cn WATER F 5.0 2.5 3.0 11.5 14.0
W
u WASTE
K
3 WATER u 7.5 6.0 0.0 1 .o 0.0
SI
DESAL.
WASTE D 44.0 43.0 38.0 46.0 51 .O
WATER
88
n
z = c c. x (6.14)
J J
Hence
Source C o n s t r a i n t s :-
RW + RB + RS t R V t RM 5 100.0 (6.16)
FW + FB + FS + FV + FM = 11.5 (6.17)
UW + UB + US + UV + UM = U (6.18)
DW + DB + DS + DV + DM = D (6.19)
Demand C o n s t r a i n t s : -
R W + FW + UW + DW 103.8 ( f r o m 6.14) (6.20)
RB + FB + UB + DB = 7.0 (6.21)
RS + FS + US + DS = 9.5 (6.22)
R V + FV + UV + DV = 0.7 (6.23)
RM + FM + UM + DM = 0.5 (6.24)
Quality Constraints:-
500 R W + 600 + FW + 1750 UW + 175 DW I 103.8 (1750) (6.25)
500 RB + 600 + FB + 1750 UB + 175 DB 5 7.0 (1750) (6.26)
500 RS + 600 + FS + 1750 US + 175 DS 5 9.5 ( 700) (6.27)
500 RV + 600 + FV + 1750 UV + 175 DV 5 0.7 ( 700) (6.28)
500 RM + 600 + FM + 1750 UM + 175 DM 5 0.5 ( 700) (6.29)
89
B I end i n g Constraints:-
U + 1.45 D = 10 ( from 6 . 1 0 ) ( 6 .3 0 )
Non-nega t i v i t y constra i n ts:-
A l l unknowns 0 (6.31 )
It is necessary to convert Eq. 6.18 and 6.19 to a form with no
unknowns on the r i g h t h a n d side. Hence:
UW + UB + US + UV + UM - U = 0 (6.32)
DW + DB + DS + DV + DM - D = 0 (6.33)
The set of constraints (Eq. 6.16, 6.17 and 6.20 to 6.33 can be
optimized, subject to the objective f u n c t i o n of Eq. 6.15, u s i n g the manual
Simplex Technique. Since the r e s u l t a n t m a t r i x i s r a t h e r l a r g e , a computer
programme to solve I i n e a r programming problems b y the Simplex Technique
i s preferred. The f o l l o w i n g i s noted :
a ) D i a g r a m m a t i c a l l y the solutions i s as i n d i c a t e d i n F i g . 6 . 3
b ) The water obtained from the water board is not required, and is
consequently a l l o c a t e d to waste as slack. This is interpreted to mean
that i s i s not necessary to o b t a i n water from the water b o a r d .
c ) Most of the ground water i s allocated to system 1 and most of the
wastewater i s a l l o c a t e d to the t r a n s f e r .
I WASTE J
1750mg/L
11 .5MP/d /
/f
Ground-
water System 1 3.0 3.0 8.7 1.6 47%
Waste-
water Waste 2.0 7.5 3.8 5.2 31 %
Waste-
water Transfer 2.0 6.0 5.3 4.5 25%
Ground-
water Transfer 5.0 2.5 1.7 3.1 -24%
(6.35)
subject to source c o n s t r a i n t s 6.16, 6.17, 6.32 a n d 6.33,
demand c o n s t r a i n t s 6.20 to 6.24,
qua1 i ty c o n s t r a i n t s 6.25 to 6.29,
b I end in g c o n s t r a i n t s 6.30,
U + 1.45 D = 10 (6.30)
Non-nega t i v i t y c o n s t r a i n t s ,
adjacent c o n s t r a i n t s f o r separable v a r i a b l e s : -
93
RW
RB 0.651 4.339 3.409 7.000
RS
RV
RM
FW 3.318 4.522 3.800
FB 2.059 0.800 1.730 2.661 3;591 7,000
FS 7.265 8.382 9.500 8.674 7.848 7.022 6.196 0.700
FV 0.535 0.618 0.700 0.639 0.578 0.517 0.467
FM 0.382 0.441 0.500 0.457 0.413 0.370 0.326
uw 3.326 3.800 3.800 8.139 7.209 6.278
UB 6.349 4.941 6.200 5.270
us 0.826 1.651 2.478 3.304 8.800
uv 0.061 0.122 0.183 0.243 0.700
UM 0.043 0.087 0.130 0.174 0.500
DW
DB
DS 2.235 1.118
DV 0.165 0.082
DM 0.118 0.059
U 6.349 8.175 10.000 10.000 10.000 10 .ooo 10.000
D 2.518 1.259
OBJ.
FN .1573 1213 698 691 633 61 2 64 9 359
(R/d)
97
I1
--1
LLL
1 6
u
7 800 900 1000 1100 1200
TOTAL D I S S O L V E D S O L I D S (mg/e)
The a b r u p t change in slope a t 600 mg/e i s a r e s u l t of the necessity f o r
p r o c u r i n g Desalinated wastewater for lower acceptable TDS values. As
water of improved q u a l i t y i s required, the amount of water requiring
d e s a l i n a t i o n incrases, a n d consequently the a l l o c a t i o n s to System 1, 2
and 3 increase.
The p r e v i o u s argument v a l i d a t e s the selection of a quality within the
reach, slightly to the r i g h t of 600 m g / t . The exact quality would be
determined a f t e r a g r a p h was d r a w n .
Each successive increase i n acceptable TDS should cause a decrease in
the t o t a l least-cost.
Not a l l the successive TDS increases manifest a decreased cost f o r the
p r e v i o u s TDS. This occurs as a result of the non-convex separable
o b j e c t i v e function. However, on inspection of Fig. 6.5, it is evident
t h a t 700 mg/4, 1000 mg/e a n d 1100 mg/e a r e in f a c t local optima, since
they do not follow the expected trend. Consequently the allocations
produced i n Table 6.7 f o r these TDS v a l u e s a r e also not true optima.
These problems may be overcome be p e r f o r m i n g a sensitivity analysis,
o r p o s s i b l y b y r e v i s i n g the g r i d a n d f u n c t i o n a l equations.
The absolute lowest cost occurs when all quality aspects a r e ignored
completely. The resulting cost is $359/d, a decrease of $34/d in
comparison w i t h the s o l u t i o n u s i n g t r a n s p o r t a t ion programming a lone.
Trends of increases, decreases a n d changes i n a l l o c a t i o n s as the TDS
v a r i e s a r e e v i d e n t from Table 6.6. Two t y p i c a l forms are:
1) An increase in allocation from 500 mg/Q to 600 mg/e and a
decrease thereafter - Groundwater to System 1 .
2) No allocation until after 600 mg/e and a steady increase
thereafter - Wastewater to System 2 .
Most of the groundwater is allocated to System 1, and most of the
wastewater t r a n s f e r r e d a t low TDS values d i s c h a r a g e d to waste.
REFERENCES
CHAPTER 7
I NTRODUCT ION
Clarl-
flcatloi
Grit
500
Sand
50
10
1
0.1
Silt
0,Ol
SAND FILTERS
Clay
0.00
EVAPORATION
Collolc
R E V E R S E OSMOSIS
,001 -
0001 I I
1 10 100 1000 10000 too
RELEVANT D E S A L I N A T I O N METHODS
I n d u s t r i a l Wastewater treatment
The increased energy costs during the last decade have directed
research and development work for al I desalination methods towards
r e d u c i n g energy consumption. Different methods of energy recovery have
been investigated and their applicability depends on the costs and the
size of the p l a n t (Binnies, 1981; Larson, 1979).
Some examples a r e :
1) RO p l a n t energy can be recovered b y installing a turbine on line in
the ( h i g h p r e s s u r e ) b r i n e stream.
2) Underground installation of RO plant can be justified based on
u t i l i s a t i o n of the s t a t i c pressure instead of h i g h p r e s s u r e pumps.
T h i s can be a p p l i c a b l e to the m i n i n g industry f o r f r e s h water production
underground. The energy consumption costs w i l I n o r m a l l y be h i g h .
In the USA it has been suggested that advanced treatment methods
(demineral i z a t i o n ) for domestic and municipal wastewater is the best
a l t e r n a t i v e f o r s o l v i n g the problems of water supply. This is applied iri
Denver (USA) where the RO for demineralization is included in a single
plant, a n d i s now b e i n g contemplated elsewhere.
Reverse Osmosis
Membrane Description
E Iectrod ia I y s i s
I o n Exchange
The ion exchange process has been used f o r many years for softening
106
iter
ter
COST ANALYSIS
C a p i t a l Costs
I n d i r e c t C a p i t a l Costs
R u n n i n g Costs
Running costs a r e g e n e r a l l y d i r e c t l y p r o p o r t i o n a l to p r o d u c t t h r o u g h p u t
and include energy costs, chemical costs, labour for operation and
ma i n tenance, membrane r e p lacemen t , o p e r a t i n g a n d maintenance costs.
For membrane plants it is reasonable to assume electricity cost with
100% load factor of 2c/kWh. Chemical cost and treatment costs vary with
feed water characteristics, the process used and the plant's recovery
ratio.
L a b o u r Costs
Membrane Replacement
D E S R L I N R T I O N PROCESSES
CFTXTRL L RUNNING co- ro11 nLcmtornmysxs L IINLRSL osnos~s
c
8
Fig
n
\
I
Fig. 7.4 Running cost versus P l a n t size (Feed) f o r ED & RO
ED process ( -.-.d
RO process ( )
I10
D E S A L I N A T I O N PROCESSES
ENERGY R f O U I R M W S (hl4hk) FOR CLCCTRODIRLYSIS h REVERSE OSMOSIS
ED process (-.-.- )
RO. p r o c c s d )
111
CONCLUSIONS
REFERENCES
CHAPTER 8
INTRODUCTION
No.3 Shaft
Flsrure water
Mlm u r v l c a water
-
-L
... . .....
C b r water return Ilne
-
-. -
Carcmd. rvstem overflow
F a d rate;
Sludge llmr
----- MIM wwklngs u r d water
6 nrsure water
SYSTEMS ANALYSIS
IDENTIFY SYSTEM
READ CONDUIT DATA
IDENTlfY NODE CONS
SEEK ML LOOPS
HYDRAULIC NETWORK
. I
I I
ANALYSIS FOR CONDUITS
I
I I I
I TRYMTERNATNEDESM
PUNT LOCATIONS AND-SCALE I
I
I
SOLVE SIMULTANEOUS
MASS W C E TDS EOS
AT ALL NODES XPT 0
ADJUST FLOW IN LOOPS
IN BEST WAY UNTIL
AU TDS's WITHIN LIMITS
RETC COST
General o p t i m i z a t i o n problem
0 ,DAY
1 2 3 4 5 6 /
RESULT OF ANALYSIS
Table 8.1 was obtained from the simulation program, where water
quality was arbitrarily assumed to start at 500 mg/P before increasing
over l i t t l e more than a week to e q u i l i b r i u m TDS's of over 4000 mg/l at
p I aces.
~~ ~
Pipe Node 1 Node 2 Water Flow Increase TDS Cost Equilibrium TDS
number (W (will (CPU (WlU
1 0 1 23 800 80 i91
2 1 3 23 0 0 193
3 3 6 8 0 0 183
4 6 I 8 0 0 3296
6 0 11 -8 0 0 4133
6 0 I -8 0 0 3296
I 2 4 36 0 0 2419
8 4 6 61 0 0 2298
9 0 6 18 1800 0 2298
10 5 8 69 400 0 8180
11 8 7 62 0 6 8295
12 8 9 67 0 0 3174
13 9 2 67 0 20 8188
14 I 10 36 0 0 8138
16 10 12 40 200 0 8826
16 I 11 26 0 0 4183
17 11 12 18 200 0 3826
18 0 10 4 1800 0 3138
19 12 8 60 0 10 3180
20 12 0 8 0 60 0
21 2 0 21 0 16 0
22 3 4 16 0 0 2419
~~
Pipe Node 1 Node 2 Water Flow Increase TDS Cost Equilibrium TDS
number (116) ( W l ) (c/W (mg/l)
1 0 1 176 800 30 800
2 1 3 176 0 0 799
3 3 6 118 0 0 798
4 6 7 118 0 0 1239
6 0 11 -8 0 0 1303
6 0 7 -8 0 0 1233
7 2 4 36 0 0 1106
8 4 6 94 0 0 1216
9 0 6 18 1800 0 1216
10 6 8 112 400 0 1604
11 8 7 113 0 6 1233
12 8 9 212 0 0 1609
13 9 2 212 0 20 1603
14 7 10 78 0 0 1489
16 10 12 194 200 0 1697
16 I 11 146 0 0 1303
17 11 12 138 200 0 1697
18 0 10 116 1800 0 1489
19 12 8 213 0 10 1604
20 12 0 118 0 60 0
21 2 0 176 0 16 0
22 3 4 68 0 0 1106
127
REFERENCES
APPENDIX 8.1
The programme MlNSlM can be copied onto new tapes. Data files
129
MlNSlM L i s t of Symbols
A1 a n g l e of v i e w i n g p l a n e from X a x i s , degrees
A2 a n g l e of v i e w i n g p l a n e from 2 a x i s , degrees
B1-9 dummy i n p u t f o r a l t e r a t i o n s
C p r i c e c/kP
c2 t o t a l cost/rands p e r annum
E( i n i t i a l volume, m3. Use n e g a t i v e o r zero v a l u e to s i g n i f y
constant outflow over 24h. Must then b a l a n c e flows over 24h
not a t peaks.
F( ) new mg/P a t node
G( p o l l u t a n t concentration mg/P a t node
G1 G counter f o r p l o t s
GO used in c a l c u l a t i o n of TDS a t nodes
H( head, m not used
H1 size of device
I( top node of p i p e
J( bottom of p i p e
K ( I ,MI number of p i p e s connecting i n t o node ( u p to 5 p e r m i t t e d )
L( 1 length, m
L1 distance to device
MO p i p e counter
M Pipe
M1 number of nodes
M2 connecting p i p e counter
M3 counter f o r i n i t i a l flow c a l c u l a t i o n s
M5 p i p e no. of device
N node counter
NO device p e r p i p e
N1 node a t w h i c h TDS i s to b e p l o t t e d
N2 1 = TDS concentration p l o t r e q u i r e d
2 = volume a t node, m3
N3 0 = o l d data, 1 = new, 2 = r e v i s e d
N4 1 = g r a p h i c d i s p l a y , 0 = none, 2 = r e c o r d d a t a a n d stop
P( 1 i n p u t p o l l u t a n t concentration a t node 2, mg/e i n t o p i p e
PO in it i a I concentration
Q( flow P / s
91 ZQ
R( design flow i n p i p e P / s
S( 1 volume m3
s1 S coun'ter f o r p l o t s
T1 d u r a t i o n of simuln. d a y s
T2 s i m u l a t i o n i n t e r v a l , hours
T3 d r a w o f f hours/day (1st hours of d a y )
T4 time i n t e r v a l s f o r d r a w o f f p e r i o d
1 30
Data Input
Ll System name
L2 Simulation p e r i o d , d a y s
L3 Time increment, d a y s
L4 Period in h o u r s p e r d a y d u r i n g which d r a w o f f occurs. The
b a l a n c e of the time water may flow to r e f i l l r e s e r v o i r s .
L5 I n i t i a l TDS of the e n t i r e systems, mg/P
L6 Node no. a t which a p l o t o f TDS versus time i s r e q u i r e d
L7 Type 0, 1 o r 2 depending on whether the o l d d a t a f i l e , a
new one o r a r e v i s i o n of the o l d one i s r e q u i r e d
L B et seq Type the f o l l o w i n g separated b y commas, w i t h one l i n e p e r
pipe o r conduit;
Program listing
1 I RE-STORE"MINS1M'
10 I MINE WATER RETIC SIMULATION I"M1NSIM"
1 1 GRAPHICS OFF
12 DUMP DEVICE I S 707,EXPANDED
1 4 PRINTER I S 707
2 0 ASSIGN O P a t h l TO "DATMIN" I CREATE "FILNAM",100,88 ag"DATM1N"
2 1 D I S P "SYSTEM NAME"i
22 INPUT NS
2 4 PRINT NE
2 5 INTEGER N ,NO ,N1 ,N2 ,N3 ,N4 ,M ,M0 ,M1 ,M2 ,M5 ,S I ,I( 99 ) , J ( 99 ) ,K( 99.5 )
2 7 D I M X ( 9 9 ) , Y ( 9 9 ) .Z( 9 9 ) ,U( 9 9 ) ,W(99) ,P( 9 9 ) , H ( 9 9 ) ,6( 9 9 ) , S ( 9 9 )
30 D I M E( 9 9 ) ,R( 9 9 ) ,P( 9 9 ) ,L( 9 9 ) .F( 9 9 ) ,C( 99 )
31 D I S P "DURN OF SIMULATION,DAYS"I
3 2 INPUT 1 1
33 D I S P "TIME 1NCREHENT.HOURS"t
3 4 INPUT 1 2
3 5 D I S P "FLO OVER HOURS/DAY"i
36 INPUT T3
37 D I S P " I N I T I A L TDS,mQ/l"(
3 8 INPUT P0
39 D I S P "PLOT TDS AT N0DE"i
4 0 INPUT N1
4 2 DEG
4 3 D I S P "OLD OR NEW OR REV D A T A ( 0 / 1 / Z ) " i
45 INPUT N3
61 NZ-1
62 T4-T3/T2 I I T S / d
66 1 7 - 2 4 / 1 2
70 G(O )=0
80 X ( 0 )-0
90 Y(0j-0
1 0 0 Z(0)-0
110 E ( 0 ) - 0
111 FOR N-1 TO 9 9
112 S ( N ) = 0
113 E ( N ) - 0
114 F ( N ) = 0
115 X ( N ) = 0
116 Y(N)=0
117 Z(N)=0
118 C(N)-0
1 1 9 NEXT N
120 c2-0
122 G( 1 )=P0
1 2 5 M1=0
130 FOR M=1 TO 9 9
1 4 0 I F N3C>1 THEN 190
1 4 5 I NEW PIPE DATA
150 DISP "N1 ,N2,X2,YZ,ZZ,U2,l/s,tmg/l,c/"i
1 6 0 INPUT I ( M ) , J ( M ) , X ( J ( M ) ) ,Y( J ( M ) ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) ,P(M) ,C( M )
170 OUTPUT @ P a t h 1.Mi I ( M ) , J ( M ) , X ( J ( M ) ) .Y( J ( M ) ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) , P ( M ) ,C(M)
1 8 0 GOT0 2 10
1 8 5 I OLD PIPE DATA
1 9 0 ENTER @ P a t h 1, M i I ( M ) , J ( M ) , X ( J ( M ) ) , Y ( J ( M ) ) , Z ( J ( M ) ) ,E( J ( M ) ) ,R(M),P(M),C(M)
2 1 0 I F I ( M ) + J ( M ) = 0 THEN 228
212 S(J(M))=E(J(M))
218 G(J(M))-P0
2 2 0 C2=CZtC( M )rR( M )*315
225 M l = M l + l
132
226 NEXT M
228 I F N3<2 THEN 256
230 FOR M0=l TO 99 I REV P I P E DATA
231 DISP "PIPE N o . " i
232 INPUT M
233 C2=CZ-C( M ) * R ( M ) * 3 1 5
234 DISP "N1 , N 2 , X 2 , Y 2 , ~ 2 , U 2 , l / s , t m g / l , c / " i
235 INPUT I ( M ) , J ( M ) , X ( J ( M ) ) , Y ( J ( M ) ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) , P ( M ) , C ( M )
245 OUTPUT B P a t h l , M i I ( M ) , J ( M ) , X I J ( M ) 1 , Y ( J ( M ) ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) , P ( M ) ,C(M)
246 I F I ( M ) t J ( M ) = 0 THEN 256
?48 I F MCMl THEN 251
249 Ml=MI+l
251 S ( J ( M ) )=E( J ( M ) )
253 G( J ( M ) )-P0
254 CZ=C2+C(M)*R(M)*315
255 NEXT M 0
256 FOR M - l TO M 1
257 L ( M ) = S Q R ( ( X ( J(M))-X(I(M)) ) ^ 2 + ( Y ( J ( M ) ) - Y ( I ( M ) ) ) * 2 + ( 2 ( J(M))-Z( I ( M ) ) ) ^ 2 1
258 NEXT M
262 PRINT "N1 N2 X2 Y 2 22 U2 Q t m g c / "
263 FOR M=1 TO M 1
264 PRINT USING 2651 I ( M ) , J ( M ) , X ( J ( M ) ) , V ( J ( M 1 ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) . P ( M ) , C ( M )
265 IMAGE 20,2D ,5D, 4D,5D ,4D ,3D, 4D,3D
267 NEXT M
268 D I S P "LAYOUT GRAPHICS(0=NO,l=YES,2=RECORD DATA & STOP ) " I
269 INPUT N4
270 I F N4<2 THEN 2 8 0
271 ASSIGN O P a t h l TO "DATMIN"
272 FOR M=1 TO 99
274 OUTPUT @ P a t h 1, M i I ( M ) , J ( M ) ,XC J ( M ) ) , Y ( J ( M ) ) ,Z( J ( M ) ) ,E( J ( M ) ) , R ( M ) ,P(M) ,C(M 1
276 I F I ( M ) t J ( M ) = 0 THEN 2910
278 NEXT M
286 I F N4-1 THEN 1570
706 ALPHA OFF
705 GINIT
707 GRAPHICS ON
708 I F N2>1 THEN 8 5 0
710 WINDOW - . 5 , T I . - 2 0 0 , 2 * G ( I )t3000
71 1 C L I P 0 ,T1 ,0 ,2*Gc 1 )+3006
720 AXES 1 ,100
730 C L I P OFF
740 MOUE T1-1,10
750 LABEL " O A Y "
760 FOR T0=1 TO T I
770 MOVE TO-.5,-200
780 LC\BEL UALO(T0)
790 NEXT T0
800 FOR G I - 0 TO 2 * G ( l ) + 2 5 0 0 STEP 500
810 MOUE -.5,G1
820 LABEL UALS(G1 )
830 NEXT 61
835 MOUE .5,2*G( 1 )+600
840 LABEL "TDSmg/l NODE"LUALB(N1 )
845 60TO 890
850 WINDOW -.5,Tl ,-10,E(NI ) t l U
851 C L I P 0 , T l .-10,E(N1 ) + I 0
854 AXES 1 , l
855 C L I P OFF
133
857 MOUE T l - 1 ,0
868 LABEL "DAY"
860 FOR T0=l TO T1
862 MOUE T0-.5,-10
864 LABEL UALS(T0 )
866 NEXT 1 0
8 6 8 FOR S l - 0 TO E ( N 1 ) t l O STEP 10
870 MOUE -.5,S1
872 LABEL UCILL(SI )
874 NEXT S1
876 MOUE B,E(Nl ) + I 0
878 LABEL "UOLn3,NODE"bUALO(Nl)
890 I END OF LABELING
900 FOR N=1 TO M l t l I NODES CONS
905 M2=0
906 FOR M0=l TO 5
907 K ( N ,M0 )-0
908 I ( K ( N , M 0 ) ) = 0
909 J(K(N,M0))=0
910 NEXT MO
915 FOR M0=1 TO M 1 1 PIPES
920 I F I ( M 0 ) < > N THEN 940
925 MZ-M2tl
927 K ( N ,M2 )-M0
930 GOTO 960
940 I F J(M0)<>N THEN 968
945 M2=M2+1
950 K(N,MZ)=MB
960 I F M2-5 THEN 970
968 NEXT M0
970 NEXT N
975 I F NZ>l THEN 990
980 MOUE 0,G(N1 )
985 GOTO 1090
990 MOUE 0,S(Nl )
1090 FOR T=1 TO T l ! DAYS
1100 FOR T9=1 TO T7
1105 I F T9>T4 THEN 1160
1110 FOR M=1 TO M1
1145 Q ( M ) = R ( M )
1155 NEXT M
1157 GOTO 1200
1160 I F E ( I ( M ) ) < = 0 THEN 1200
1195 Q(M)-0
1200 TDS a STORAGE TINS
1210 FOR N-1 TO M 1 + 1
1220 I F E ( N ) > 0 THEN 1340
1225 Gl=Q
1228 Q1=.001
1230 FOR M0-1 TO 5
I 2 3 5 I F J ( K ( N , M B ) ) O N THEN 1280
1236 I F G( I ( K ( N , M 0 ) ) )+P(K(N ,M0) )<=0 THEN 1255
1240 G1=( G( I ( K ( N ,M0) ) ) t P ( K ( N , M I ) ) )*Q(K(N .M0) )+G1
1255 01 =Ol + Q ( K ( N ,M0 ) )
1280 NEXT MO
1282 G l = G l / Q l
1285 GOTO 1440
1340 G l m 0 1 S>0
1345 G0=G( N )*S( N )
134
2262 H l = E ( J ( N ) ) / 2 5
2264 C0=0
2266 GOT0 2320
2268 ALPHA ON
2 2 7 0 DISP "PIPEn ,X ,TYPE ,SIZE .COST/' 1
2280 INPUT M5,Ll ,T5,H1 ,C0
2 3 2 0 I F M5-0 THEN 2420
2340 X5=X( I ( M 5 ) ) + L l / L ( M 5 ) * ( X ( J(M5 1 )-X( I ( % ) ) )
2350 Y5=Y( I ( M 5 ) )+Ll /L(M5 ) * ( Y ( J( M5 ) )-Y( I(M5 ) ) )
2360 Z5=Z( I ( M 5 ) )+L1 /L(M5 )+(Z( J(M5) )-Z( I ( M 5 ) ) )
2 3 7 0 US=XS*COS(AI ) t Y 5 * S I N ( A 1 )
2 3 8 0 W5=Z5*COS(A2 )t(YS*CDS(Al )-XS*SSN(Al ) ) * S I N ( h Z )
2390 ON T5 60TO 2460,2490,2540,2590,2850
2400 I l=FLANGE,Z=UALUE,3=TANK,4=~RROW,S=S~UARE
2 4 1 0 NEXT N
2415 NEXT N0
2420 MOUE U0,WB
2430 C2=1NT(C2 )
2440 LABEL " R/s="LUALO(C2)
2445 DUMP GRAPHICS
2 4 5 0 60T0 7 0 0
2460 MOUE US ,W5tHI / 2
2 4 7 0 DRAW US ,W5-H1/2
2480 GOT0 2410
2 4 9 0 HOVE U5-H1/2,WStH1/2
2500 DRAW U 5 t H I 12 ,WS-HI / 2
2510 MOUE U5+Hl/2,WS+H1/2
2520 ORAW U5-H1/2,U5-H1/2
2530 GOTO 2410
2540 MOUE U5-HI / 2 ,WStWI
2550 DRAW U 5 - H l / 2 ,W5
2560 DRAW U S H 1 / 2 ,WE
2 5 7 0 DRAW U5+Hl/Z,W5+Hl
2580 60TO 2410
2590 I F U ( J ( M S ) ) < > U ( I ( M S ) ) THEN 2601
2591 IF W ( J ( f l S ) ) > W ( I ( M 5 ) ) THEN 2594
2592 U8=270
2 5 9 3 GOTO 2608
2594 U8=90
2595 GOT0 2608
2601 UE=ATN( (U( J( M 5 ) )-W( I( M5 ) ) ) / (U( J ( M 5 ) )-U( I( M5 ) ) ) )
2602 I F U8>=0 THEN 2606
2603 I F W ( J(M5) ) < W ( I(M5) 1 THEN 2608
2604 GOTO 2607
2606 I F U ( J ( M S ) ) > W ( I ( M S ) ) THEN 2608
2607 U8=U8t180
2608 UG=U5-Hl*COS(U8-45)
2610 W6=W5-Ht*SIN(U8-45)
2 6 2 0 U7-U5-Hl*COS(U8t45)
2630 W7-W5-Hl*SIN(UEt45)
2 8 1 0 MOUE U6,W6
2820 DRAW U5,W5
2830 DRAW U7,W7
2 8 4 0 6QTO 2410
2850 MOUE US-H1/2,WStHl
2860 DRAW U5-H1/2 ,W5
2a70 DRAW u 5 + ~ 12 1 ,w5
2880 DRAW U5tH1/2,W5tHI
2890 DRAW U5-H1/2,W5tHl
2900 60TO 2410
2920 END
136
p r i c e c/ke
cost
t o t a l cost
dummy
TDS
max. TDS desired, o r G of node w i t h m a x . T D S
increment i n TDS
max. increment i n TDS
t o t a l TDS - mg/s i n t o node
t o t a l flow i n t o node
TDS
G(I) - H(I)
max. TDS
top node
bottom node
p i p e no. connecting to node ( u p to 5 p e r m i t t e d )
number of p i p e connecting
number loops
best loop no.
branches i n loop
p o s i t i v e loop
p i p e s out node
loop counter
number loops a n d b e g i n number loop
number connecting p i p e s to node
p i p e number
p i p e counter
number of nodes
number of p i p e s in loop
number of connecting p i p e s from node
dummy
p i p e s out node
number of p i p e s in loop
r e d u c t i o n i n no. p i p e s in loop, o r , p i p e to node w i t h max.
TDS
p i p e number i n loop
b e g i n p i p e f o r loops
name
node counter
no. nodes
p i p e no.
i n p u t TDS, mg/e
flow e/s
dQ/dC
dQ
co-ord. not used i n MINOP
0,
I1
137
Notes on program
The program i s in BASIC for an HP 9816 series 200 micro computer. The
data f i l e i s obtained from the M l N S l M program in appendix 8.1.
138
450 NEXT M0
4 6 0 NEXT N
4 7 0 G( 0 1-0
4 8 0 L1=01LOOFS
4 9 0 FOR.M9=1 TO M1 IBEGINPIPE FOR LOOPS
500 L 0 - L l t l ITRY LOOP
51 0 L8=0
520 M G ( L 0 ) = l I N O . P I P E S I N LOOP
530 M2(L0,1 )=M91PIFES I N LOOPI
540 LE=LBIPOS LOOP
550 FOR L 3 = 1 TO MlIBRANCH ROUTINE
560 L8=0
570 L4=M3( J ( M 2 ( L0 ,M6( L0 ) 1 ) )
580 FOR M5=1 TO L 4 I P I P E S OUT NODE
139
1190 62=62tQ(K(N,O))
1200 NEXT D
1210 63=6(N)
1220 6(N )=GI 162
1230 64=AES( 6( N )-63)/G(N )
1240 I F 64<65 THEN 1260
1250 6544
1260 NEXT N
1270 I F 65<.00I THEN 12901 MAX FC
1280 NEXT L 4
1290 RETURN
1300 FOR M=1 TO N1 IITNS
1310 RI=BIDQ/DC
1320 R3=0IDQ
1330 FOR L 2 = l TO LllEEST LOOP
1340 C1=0
1350 H1=0
1360 M0=6
1370 FOR M8=1 TO MSfL2 )
1380 C1 =C1tC(MZ(L2 ,M8 ) )
1390 I F Q ( M 2 ( L Z 3 M 8 ) ) < 0THEN 1540 INEXT LOOP
1400 M0=M8
14 10 Q( M2( L2 ,M8 ) )=Q( M2( L2 ,M8 ) ) t1
1420 I F J(M2(L2,M8))=0 THEN 1470 !NEXT PIPE
1430 I F G( J ( M 2 ( L2 ,M8 ) ) )-H( J( M 2 ( L Z ,M8 ) ) )<=HI THEN 1470
1440 HI =G( J( M 2 ( L2 ,M8 ) 1 )-H( J ( M 2 ( L2 ,M8 ) ) )
1450 M7=MZ(L2,MB)IPIPE TO NODEWITH MAX TDS
1460 60=G(J(MZ(LZ,M8)))
1470 NEXT ME
1480 I F H l i = l THEN 1540
1490 GOSUE 1120
1500 I F ( G Q - G ( J ( M 7 ) ) ) / C I ~ ' = R 1THEN 1540
1510 Rl=(G0-G( J(M7) ) ) / C l
1520 R3=(G0-H( J(M7)) ) / t G Q - G ( J ( M 7 ) ) )
1536 L3=L2
1540 FOR M8=1 TO MQ
1550 Q( MZ(L2 ,ME) )=Q(MZ(L? , M a ) ) - 1
1560 NEXT M8
1570 NEXT L t
158@ FOR M7=1 TO M6(L3i
1590 Q ( M ~ ~ L ~ , M ~ ) ) I Q ( M ~ ( L )~t ,RM3 ~ )
1600 NEXT M7
1610 NEXT M
1620 C 2 = 0
1630 PRINT " P n N l N2 1 / s tTDSng1 c / L 1 T D S 2 '
1640 FOR M-1 TO M1
1650 PRINT USING 1660;M,I(M) , J ( M ) , Q ( M ) ,P(M) , C ( M ) ,G( J(M))
1660 IMAGE 2D ,4D ,4D ,4D ,6D ,50 ,SO
1670 C2=C2tC(fl)*Q(M)*315
1680 NEXT M
1690 C2=INT(C2 )
1706 PRINT "COST ,R/a="i C Z
1710 ASSIGN 0Pathl TO
1720 EN0
141
CHAPTER 9
INTRODUCTION
The international growth i n water demand over the last few decades
has been p e r s i s t e n t l y h i g h . T h i s r a t e of growth i s l i k e l y to c o n t i n u e a s a
l a r g e p r o p o r t i o n of the p o p u l a t i o n is increasing r a p i d l y in standard of
living. The a v a i l a b i l i t y of new sources of water h a s tended to l a g b e h i n d
demand. Even i f new sources were a v a i l a b l e the cost of p r o v i d i n g to meet
any possible drought extreme in developing areas could be high, on
account of the u n r e l i a b i l i t y of r i v e r flow. Surface water can o n l y be used
t o i t s f u l l e s t extent i f a l t e r n a t i v e sources a r e a v a i l a b l e to meet essential
demands i n times of d r o u g h t . For t h i s reason the w o r l d i s now l o o k i n g to
groundwater a n d wastewater to meet s h o r t f a l IS.
A scheme i s i n v e s t i g a t e d here to s u p p l y from groundwater o n l y a t the
r a t e a t which i t can be n a t u r a l l y replenished. Separate studies a r e b e i n g
conducted on t e r t i a r y treatment of wastewater but the idea of using the
wastewater to a r t i f i c i a l l y recharge groundwater with wastewater is only
now r e c e i v i n g consideration. Such research is a long term project and
cannot be expected to re1 ieve c u r r e n t droughts, w h i c h however p r e c i p i t a t e d
research i n t o a l t e r n a t i v e sources of water.
It i s proposed to use groundwater i n conjunction w i t h surface water
resources in such a way that deficiencies in surface water can be
supplemented b y groundwater, resulting in a higher overall availability.
Surface water resources can then be u t i l i z e d to a g r e a t e r degree since
groundwater reserves can be d r a w n on in times of shortfalls in surface
r i v e r s (Paling, 1984). The r a t e of recharge w i l l also be r e l a t i v e l y slow
owing to limited suitable wastewater being available, the possibility of
natural purification and the limited permeability of the soil. The
h y d r a u l i c s of the recharge process should be i n v e s t i g a t e d w i t h s i t e tests.
The case study a n a l y z e d i s the Witwatersrand area, a h i g h growth r a t e
conurbation based o r i g i n a l l y on m i n i n g . Groundwater constitutes a t present
o n l y one percent of the average d a i l y supply to t h e Rand Water Board of
2400 M l / d to the Witwatersrand area. Privately owned boreholes for
farming and gardening purposes are however in common use as a
I
COST ANALYSIS
20 207.5 * Q + 4625
40 217.5 * Q + 4725
60 230.0 * a + 5800
The costs for the pump motors were based on the shaft power,
increased b y a 20% safety m a r g i n (see Table 9.3).
0 - 250 400 60
250 - 1900 1000 - 3300 70
1900 - 6570 6600 - 11000 80
-
a p p l i e d to solve t h i s problem i n the f o l l o w i n g way.
Rynfield 0---------
--------ODaveyton
,'; \, -0' i
MC Comb
Maple1:on
RESULTS
+' nodes 1 , 2
REFERENCES
Bouwer, H., Rice, R.C., Lance., J.C., and Gilbert, R.G., 1980.
R a p i d - i n f i l t r a t i o n Research a t F l u s h i n g Meadows Project, Arizona. J.
Water Polut. Control Fed., Vol. 52, No. 10, p 2457.
I d e l o v i t c h , E. and M i c h a i l , M., 1984. S o i l - a q u i f e r Treatment - A New
Approach to an O l d Method of Wastewater Reuse. J. Water P o l l u t .
Control Fed., V o l . 56, No. 8, p 936.
Leighton, J.P. and Shoemaker, C.A., 1984. An I n t e g e r Programming
Analysis of the Regionalization of L a r g e Wastewater Treatment a n d
Collection Systems. Water Resources Research, Vol. 20, No. 6, p 671.
154
CHAPTER 10
I NTRODUCT ION
SUE-SKTEM I
'. \
\
Fig. 10.2 Water c i r c u l a t i o n d i a g r a m
‘1 ‘1-2
‘2’2-3
+ ’ ‘3‘3-6 + ‘4‘4-3 (10.1)
Note that Q6-7 is fixed so its cost is not variable i.e. need not be
considered. The cost of conveyance i n n a t u r a l channels i s zero.
P u r i f i c a t i o n costs comprise a component p r o p o r t i o n a l to flow Q . . and a
I-J
component proportional to pollutant load Q. .P.
I-] I-J
. or proportional to
pollutant load removed, Q . . ( P . .-P.
I-J I-J
. )
l-J/2
where subscript 2 refers to
conditions after treatment. The l a t t e r component, i.e. cost p r o p o r t i o n a l to
load removed, i s d i f f i c u l t to e s t a b l i s h a n d i n fact i t i s often assumed that
works a r e designed to produce a n effluent of reasonable standard with
treatment costs p r o p o r t i o n a l to flow r a t e . Vie w i l l consider the general case
a n d designate the coefficients of Q
3-6 P3-4 a n d Q2-6 P2-3 as C5 and C6
respectively. Since Q3-6 + Q3-4 is a constant, the flow-proportional cost
can b e omitted.
The equations o r c o n s t r a i n t s d e s c r i b i n g the system a r e formulated next.
Hall (1977) formulated the system w i t h similar constraints, but at that
time was unaware of a simple method of solution.
For flow b a l a n c e a t the v a r i o u s nodes:
At some source nodes such as (1 1, the yield may be limited but in our
case we consider u n l i m i t e d augmentation possible, a t a cost. At consumer
nodes the s u p p l y must be s u f f i c i e n t :
(10.16)
(10.17)
Q3-4P3-4+Q3-6P3-4-Q2-3P2-3-Q4-3p5-3 = al
2 (10.18)
Note P2-6 e q u a l s P2-3, P4-5 equals P4-2 a n d P3-9 e q u a l s P3-4 so these
s u b s t i t u t i o n s a r e made f o r s i m p l i f i c a t i o n .
It is implicit in the optimization program that all variables are
non-negative a n d r e a l so these c o n s t r a i n t s a r e not s t a t e d e x p l i c i t l y .
The general model then is to minimize (10.1) subject to constraints
(10.2)-(10.17). A method of solution i s o u t l i n e d in the n e x t section.
162
OPTIMIZATION METHOD
QP = M‘ - NZ (10.20)
and
Q = M + N , P = M - N (10.21)
M and N a r e u n r e s t r i c t e d i n s i g n b u t t h i s i s p e r m i t t e d in the d e l t a method
of separable programming (IBM, 1976).
The separable programming algorithm is based on the fact that a
separable f u n c t i o n can be approximated b y piecewise I i n e a r functions. The
p o l y g o n a l a p p r o x i m a t i o n i s represented b y a set of special variables, so
a n y v a l u e of M can be represented a s follows:
M = Mo + GIDl + G2D2 + ...Gk D k +...GRDR (10.22)
where the Ds represent i n t e r v a l s of M a n d the special v a r i a b l e s GI, ...,GR
a r e defined as follows:
f o r M in i n t e r v a l k ,
G , = G 2 = G k-1 = 1 (10.23)
O(Gk( 1 ( 10.24)
Gk+l = Gk+2,..GR = 0 (10.25)
i.e. M comprises a set o f i n t e g r a l i n t e r v a l s D u p to k - 1 p l u s a f r a c t i o n
of i n t e r v a l k .
Note that
M’ = Ma
0
G E+
1 1
+...
G E +... GRER
k k (10.26)
where each i n t e r v a l Ek corresponds to an interval Dk. Thus for the
approximations i n Fig. 10.3, Mo = 0, k = 4 a n d Gk = 0.3. The v a l u e of M
can be confined to a known range.
Although there i s a p o s s i b i l i t y of a t t a i n i n g a local optimum t h i s chance
i s reduced i f the problem i s solved w i t h the special v a r i a b l e s set initially
at their upper bound, a n d then w i t h them set at their lower bounds, to
verify the r e s u l t s . It will be observed that 2R variables are introduced
i n t o the model f o r each p r o d u c t in the o r i g i n a l c o n s t r a i n t s . It i s therefore
d e s i r a b l e to s i m p l i f y the o r i g i n a l system a s much a s p o s s i b l e to minimize
163
subject to:
0.28)
0.29)
0.30)
0.31)
0.32)
(10.33)
(10.34)
( 10.35)
( 10.36
(10.37)
(10.38)
Put
Q3-4p3-4
= M Z - Nt (10.42)
Then (10.37) a n d (10.38) c a n b e r e p l a c e d b y e q u a t i o n s (10.43)-(10.52):
164
(10.43
(10.44
(10.45
(10.46
(10.47)
(10.48)
(10.49)
(10.50)
(10.51)
(10.52)
t /
REFERENCES
CHAPTER 1 1
I NTRODUCT ION
HYDRAULIC ANALYSIS
FLOW MEASUREMENTS
H i g h e r income r e s i d e n t i a l
I I
0 3 6 9 12 15 18 21 24
Time of day h
Detailed land use d a t a were not available for the study of the low
income residential area (Fig. 11.2). Details of houses and flats were
a b s t r a c t e d from c o n s t r u c t i o n d r a w i n g s a n d checked on s i t e before the sewer
d a t a were used. The a r e a chosen embraced most of the newer sections of
Lenasia - a town 25 km south o f Johannesburg - a n d the sewage flow was
monitored with a flume. A time l a g of one hour was allowed for when
comparing the h y d r o g r a p h measured with the hydrograph at the p o i n t of
origin.
15
-
-- - Calculaled l b w
Monday 30 January 1984
... Tuesday 31 January 1984
........... Wednesday 1 February 1984
5
a
E
l 1
0 3 6 9 12
Time of day h
15 18
. 21 24
Apartment buiIdi n g s
10 - - Calculated (low
---
-.-
Monday 21 November 1983
Tuesday 22 N o v e m k 1983
Time 01 day h
Commercial a r e a s
I I
0 3 6 9 12 15 18 21 24
Time 01 day h
F i g . 11.4 Comparison of c a l c u l a t e d a n d metered flows in a commercial area
Industrial
CONCLUSIONS
300 -
--
-..-.-.-
Calculaled llow
Monday 5 July 1982
Tuesday 6 J U I ~1982
Wednesday 7 July 1982
....... Thursaay 8 July 1982
,7 E FE R EN C ES
APPENDIX 11.1
PROGRAM SEWS IM
Ey
0 15*x
~0.013a100 - = 15m
0.003'7
.'. x = t2m
Non-Circular Conduits
A
Q = A v =
- De2 Qe
4
I n f l o w Components
DATA
Program Output
16 DlJMP O E V I C E I S Nd
18 PRINTER I S Nd
20 D I M M ( 4 0 0 ) . M d ( 4 0 0 ) . D ( 4 0 0 ) .S(408). X ( 4 0 0 ) .Hq(400) . Y ( 4 0 0 ) .H(40GI). I t (400) . T x (400)
. T I (400) . Q ( 4 0 0 ) . Q c ( 4 0 0 ) . Q w ( 4 0 0 ) , G ! s ( 4 0 0 ) ,Qi (400) .Ql (400) . Q m ( 4 0 0 ) .Qv(4P)0) . J d ( 4 0 B )
3v1 DJM I21 (99) . Q 2 ( 9 9 ) , Q 3 ( 9 9 ) . T 1 8 ( 9 9 ) . T 1 1 (99) . T 2 0 ( 9 9 ) . T 2 1 (99) . T 3 0 ( 9 9 ) , T 3 1 (99) . M h ( 9
9 ) , Jh (99) ,at (29.99) ,He(99) .G1 (400)
31 COM NSCZ03
40 I N P U T "NAME '?".NS
60 READ N s . T s . T i . N h , O l l !NO.SECNS. S I M L N h . T I N C h .Nhydqphs.GLmBOTMH
70 FOR Jn=l TO Nh ! HYDROGHAPHS
80 HEAD M h ( J n ) !PIPE NO OF HYGPHS
90 NEXT Jn
100 J2=B
110 FDR K = l TO Ns ' SECT DATA
I20 READ Mb,Itk,Np,Am,Fw,Fs,Fi,Fl ! HOTM MH, ZONE TYPE,Npipes,MANNINljn.I-/MIN/HEQ
,mm/h STORM,INFIl/MIN/M/M.LK/MIN/HEQ
130 J2=J2+Np
140 M d ( J 2 ) = M b
150 FOR J=J2-Np+l TO J2 ! P I P E DATA
160 READ M ( J ) , D ( J ) . X ( J ) , S ( J ) , H q ( J ) ,Y(J) , H ( J ) ,G1 (;I) !NO. ,DIClmm,I..m,Sm/m.H~USE EQU
VS(l00m2),DETTHm,DROP BOTMENDmm,GLmMH
180 I F J < = J Z - N p + l THEN 200
190 M d ( J - l ) = M ( J )
200 I t . ( J ) = I t k
205 Q v ( J ) = K
260 T x ( J ) =X ( J) +Am*4". ( 2 / 3 /) ( D ( J )/ 1000) ." ( 2 / 3 )/S( J)". 5 ! L..AG, 5
*
270 Qc ( J) =. 7 8 5 / A m / 4 . " ( 2 1 3 ) (D( J) / 1000) ( 8 / 3 +S
,*% ) ( J) 5+ 1000 !CAPAC I TY ,L./s
280 Qw (J)=Fw+Hq ( J ) /h0/ 1000 !PEAK M3/S INFLOW
290 Q s ( J ) = F s * X ( J ) * l 0 0 / 1 0 0 0 / ~ 6 0 0 ! DO. STORM
300 Qi ( J ) = F i + X (,I) +D ( J )/ 1 0 0 0 / 6 0 / 1000 ! DO. I N F I L T N
310 Q1 ( J ) = F l + H q ( J ) /60/1000 ! DO. L E A K S
320 NEXT J
330 0 (K)=J2-Np+l ! TEMP. TOPMH FOR PLAN P1.OT
340 9rn (F:) =Np
350 NEXT K
400 T40=12 !STORM START h
410 T41-17 !PEAK h
420 91(1)=.7 ! RESID4UPCLAS. 1S T P E A K = l
430 92(1)=.5 ! 2 N D PEAK
440 93(1)=.6 ! 3 R D PEW
450 T10(1)=5 !START h 1ST PEW MUST BE 1st HG TO START
460 Tll(i)=9 ' P E A K h 1ST PEAK
470 T20(1)=7 ! S T A R T h ZND PEAK
480 T21(1)=13 !PEAK h 2ND PEAK
490 T30(1)=14 !START h 3RD PEAK
500 T 3 1 ( 1 ) =20 !PEAK h 3RD PEAK
510 91(2)=.8
520 92(2)=.7 !RESIDLPOOR
530 c-13(2)=.6
540 T10(2)=6
550 T11 (2)=8
560 T20(2)=7
570 T Z 1 ( 2 )=13
580 T30(2)=16
590 T 3 1 ( 2 ) = 2 0
4500 01(3)=.4 ! INDUST
6llb 02(3)=.S
620 93 ( 3 ) =. 45
630 T10(3)=6
640 T 1 I ( 3 )= I 0
650 T 2 0 ( 3 ) = 4
660 T 2 1 (3)=13
670 T 3 0 ( 3 ) = 1 3
675 T 3 1 ( 3 ) = 1 6
680 01 (4)=.3 !COMMERCIAL
690 02(4)=.7
700 03(4)=.45
710 T 1 @ ( 4 ) = 6
--.
720 1 1 1 ( 4 ) = 1 0
/d T 2 0 ( 4 ) = 6
740 T 2 1 ( 4 ) = 1 3 d
W
w
t);
i
B 61
t.4 4
m IXI
2 z
W iii
I I
!- L
w
a
+ r.l
a -
1030 NEXT E
1032 GINIT
1033 GRAPHICS ON ! DRAW LAYOUT
1034 GCLEAR
1035 WINDOW - 1 0 0 0 . T l ( N l ) + 1 0 0 0 , 0 , N s + l
1039 FOR E = l TO Ns
1040 FOR J = B ( E ) TO O m ( K ) + O ( # ) - l
1041 MOVE T 1 ( N l ) - T l ( J ) . K
1042 DRAW T 1 ( N l ) - T l ( J d ( J ) ) - 2 0 . K
1043 L l = T l ( N l ) - T l (J)-20
1044 C S I Z E 3..4
1045 MOVE L1 .K
1046 ! D I R l+PT/b
1047 LABEL M ( J )
1048 LDIR 0
104.9 NEXT J
1050 Jmm=Om ( K ) +O ( K ) - 1
1051 I F J d ( J m m ) = 0 THEN 1055
1052 MOVE T 1 ( N l ) - T l ( J d ( J m m ) ) - 1 0 . K
1053 DRAW 1 1 (N1 ) - T 1 ( J d ( J m m ) ) , G ! v ( J d ( J m m ) )
1055 NEXT K
1060 PAUSE
1070 FOR J=1 TO J2
1071 (;7m(J)=0 ! MAX. FLOW
1072 Qv(J)=0 !SPILL VOL.
1073 NEXT J
1074 Ti=0 !COUNTER FOR T H
1075 FOR T h = T i TO T s STEP T i ! T I M E h AT BOTTM
1076 FOR J=1 TO 52 I N I T I C I L I Z E FLOWS
1077 R(J)=0
1078 NEXT J
1080 FOR J=1 TO 52
J 090 T S T h - T I ( J ) /3h00 !PIPE T I M E FOR REACHNG E X I T A T T h
1092 I F T : : . T 1 0 ( I t ( J ) ) THEN 1100
1093 T=T+24
1100 IF T > = T l B ( I t ( J ) ) THEN 1110 !HYOROGRAPH ORDINATE PER PIPE
1290 NEXT J 5
1300 NEXT J
1 3 1 0 FOR J = 1 TO J Z
1320 I F R(J)<=Qm(J) THEN 1340
1 3 3 0 Qm(J)=Q(J)
1 3 4 0 I F Q ( J ) < P c ( J ) / l 0 0 0 THEN 1360
1350 R V ( J )=Qv (J)+ ( Q ( J )-Qc ( J )/ 1000)+3600+Ti
1360 NEXT J
1370 T i = T J + l
1380 FOR Jn=l TO Nh !HYDRAPH POTNTS
1 3 Y 0 R t (Jn,T j )=a ( J h ( J n )) + l B 0 0
1400 NEXT J n
1.410 NEXT T h
1420 PRINT "SEWER NETWORK ANALYSIS , "
1430 PRINT N8
1440 PRINT 'I P I P E D I A SLOPE HUNITS LAGh M A X L s QCAP XFIJL.1. OFL.Om3"
1450 IMAGE DDDDDD.DDDDD.D. DDDD.DDDDD.DDD. DD.DDDI)D.T)I)DDD.onnI>. D.DDDD. D
146R FnR J=l TO 52
1470 T 1 ( J ) = T l (J)/360cI
1400 Qm (J ) =Qm (J)+1800
1490 P f = a m (J) /Qc (J ) + 1 0 0
1500 PRINT USIND 1450:M(J) ,D(J) . S ( J ) .Hq(J) .T1 ( . 7 ) .Qm(,l) . R c ( , l ) . P f . R v ( J )
1 5 1 0 NEXT 3
1520 IMAGE DDDDDDDD.#
1525 IMAGE DDDDD.DD.#
1526 IMAGE DDDDD.DD
1 5 2 7 IMAGE DDDDDDDD
1530 PRINT
1540 PRINT "SELECTED HYDROORAPHS AT":Ti ;"h 1NTVL.S STCIHTTNG A T - T l FKIR PTPF
1 5 4 5 PRINT USING 1520;Ti
1550 FOR Jn=l. 7'0 Nh-1
1551 PRINT USING 1520;M(Jh(Jn))
1552 NEXT J n
1 5 5 3 PRINT USING 1 5 2 7 : M ( J h ( N h ) )
1 5 6 0 FOR T k = l TO T i
1565 PRINT USING 15'20:Tk
1570 FOR Jn=l TO Nh-1
1571 PRINT I.JSING 1525:Qt (-7n.Tk)
1 5 7 2 NEXT J n
1573 PRINT IJSING 1526:Qt ( N h . T k )
1575 NEXT TIC
1380 GCLECIR
1600 FOR Jn=l TO Nh
1602 QCLEAR ! P L O T HYDROORCIPHS ON SCRN- DUMP GRAPHICS & / O R CONT
1603 CSIZE 4
1605 Omi=Qm(Jh(Jn))
1610 WINDOW -3,24.5,-0.Qmj*l. 1
1620 CIXES 1,1,0,0,12,10
1630 MOVE l , Q m ( J h ( J n ) )
1700 *
Qm (Jh (Jn ) ) = I N T (Qm ( Jh (Jn ) ) 1000)/ 1000
1701 L A B E L Qm (Jh (Jn) ) ;" L / s PIPE" :Mh (Jn)
1710 MOVE 6,0
1711 LCIBEL. N S
1720 MOVE 22.0
1721 L A B E L "h 24"
1724 MOVE -2.10
1725 LhBEL 10
1730 FOR T=2 TO T i
1740 Th=(T)*Ti-Tl (Jh(Jn))
1750 MOVE T h - T i .Qt ( J n , T - l )
1755 DRAW Th,Qt ( J n . T )
1760 NEXT T
1765 MOVE -T1 ( J h ( J n ) ). Q t ( J n , T - 1 )
1766 DRAW T i - T l ( J h ( J n ) ) , Q t ! J n , l )
1770 PAUSE ! T Y P E CONT (El-) TO DO NEXT t i Y B 13R DUMPFH TO DRAW EX
1780 NEXT Jn
1790 FOR 1 = 1 TO 6
1800 He(I)=0
1810 NEXT I
1820 FOR J = 1 TO 52
1830 H e ( 1 t ( J ) ) = H e ( I t (J)) + H q ( J )
1840 XleX1+X(J)
1850 NEXT J
1855 PRINT
SEWSIM DATA FORM
--
NAME
2 W ! DATA ..........
NO. OF SECTIONS, SIMULATION DURATION h, T I M E INCREMENT.h, NO. HYDROGRAPHS REQUIRED, G.L.(M) BOTM M i .
BOTTOM MH NO., ZONE TYPE (1-41, NO. PIPES I N SECN, MANNINGN, SW/min/HE, STORM Chnm/h, I N F I L T N L/mm/m/m, EAKG L/mi /HE.
2030 DATA ...........................................................................................................................
P I P E DATA: ( I L PER P I P E )
P I P E NO.(=TOP MH NO.), m
-, LENGTH m , SLOPE m/m, HE=HDWE E a U I V S ( l W m ’ ) , (DEPTH TOP MHm, DROP BOTM mm, GLm
SAMPLE DATA F I L E
7000 DATA 4 , 2 4 , 1 . 4 , 5 0
2001 DATA lil,1'21,211,.311., 1 1 . 3 , 1 , 2 ,. 0 1 J , 1 1 , 0 ,. , 0
ZQCI? WTA 111,100,100 ,.0100,100,1.30,49
2061.3 DATA 1 1 2 , 3 0 0 , 3 0 0 , . 6 1 1 0 0 , l ~ 0 , 1 , 0 , 4 8
2804 DATA 112,2,9, . 0 2 , I , 0 , . 1 , . 0
288'3 DATA 1 2 1 , 1 5 0 , 1 B 0 , . 0 1 5 , 1 0 0 , 1 , 0 , 4 7
2806 DATA 1 2 2 , 1 5 0 , 3 0 0 ,.004,75,1,C!Il46
2807 D(?ITA 123,200,i00,.802,95,1.5,30,45
2 0 0 8 DATA 122,3,1,.02,1,0,a..1
2089 DATA 211,150,100,.81,100,1,0,48
2 0 1 0 DATA 123,4,1,.02,1,0,0,0
2 0 1 1 DATA 311,100,iQ0,.002,100,1,0,47
100pI0 END
:YF'B2'mX. 78B
VOLUME LABEL: B982.5
F I L E NAME PRO TYPE FiEC/FILE HYTE/HEC ADDPESS
Note _-__
SAME "SEWDAT", don't STORE ,I...
188
- h 24
190
CHAPTER 12
A w o r k i n g memory w i t h a l l d a t a
A r u l e base
A n i n t e r p r e t e r to choose and a p p l y p r o d u c t i o n s
Improved control i s achieved b y a l t e r i n g the r u l e base o r a d d i n g new
ones. For each u n s a t i s f a c t o r y production a l i s t i s created.
A meta production systems was f u r t h e r added. Meta p r o d u c t i o n s do not
a f f e c t the w o r k i n g memory b u t can change the content of the r u l e base.
The meta system i s e v a l u a t e d b y the control interpreter. A simple example
demonstates the technique:
Stormflow could be stored in a detention basin freely, while street
f l o o d i n g would be an u n s a t i s f a c t o r y state. T h i s r u l e could be described b y
a meta f u n c t i o n as flows:
( W E > 1.0) + (pump too L O W ) (Value = - 1 )
Whenever the water level in the sewers i s h i g h e r t h a n manhole level which
may cause street f l o o d i n g t h i s r u l e i s a p p l i e d . At a n y selected time i f the
meta production rule is applicable the decision PUMP = OFF is counter
r u l e d i.e. the corresponding productions a r e decreased b y 1 .
The facts i n the w o r k i n g memory a t the time may h a v e been
W E = 0.4 where W E = water e l e v a t i o n
R I = 10 where R I = r a i n f a l l i n t e n s i t y .
Another s i t u a t i o n may also have been stored in the experience memory,
e.g.
WE = low
R I = 10
I t i s possible b e t t e r productions c o u l d h a v e been a p p l i e d .
The total l i s t i n memory may now be
Va Iue WE
- -
RI
-3 LOW 9
-1 LOW 10
-6 LOW a
-10 LOW 11
The same problem was simplified into a linear system for direct
optimization at discrete times. Sewers were lumped into three
subcatchments.
Wser I1
Hbsssliisc
I2
-1
t PM PRl Po1 93
_--- ,---------------- -----
----
C 0.0 1.0
.---------------- 10
1280
0.26
-----
1730
1 1630 0 0
2 14268 1277 3800 1277 0 1280 1730
3 'I4268 10000 3600 9414 691 1260 1730
4 14268 10000 3600 Be68 9000 0 1730
6 ill88 10000 3600 6346 6346 1260 1730
6 14268 10000 3600 2946 2946 1260 1730
7 14268 10000 3600 1373 1373 1260 1730
8 14268 10000 3800 603 603 1260 1730
9 14264 10000 $600 a 0 1260 1730
10 13814 10000 3600 0 0 1260 1730
11 13049 10000 3600 0 0 1260 1262
12 13088 10000 3600 0 0 1260 193
13 11037 10000 3600 0 0 1260 0
14 10108 10000 3800 0 0 1260 0
16 9125 10000 3600 0 0 0
16 a i r 4 10000 3800 0 0 0 0 0 1260 0
17 7163 10000 3600 0 0 0 0 0 1260 0
18 6182 10000 3800 0 0 0 0 0 1260 0
19 5201 10000 3600 0 0 0 0 0 1260 0
20 4220 10000 3800 0 0 0 0 0 1280 0
21 3231 10000 3600 0 0 0 0 0 1260 0
22 1612
-~~~ 9549 3600 0 0 0 0 451 171 0
23 531 6560 3600 0 0 0 0 989 171 0
24 531 6490' 3600 0 0 0 0 2070 171 0
25 531 4420 3800 0 0 0 0 2070 171 0
28 531 2350 3600 0 0 0 0 2070 0
27 531 260 3600 0 0 0 0 2070 171 0
28 531 0 1610 0 0 0 0 280 171 0
29 531 0 1530 0 0 0 0 0 171 0
t R3 PR3 F3
.---__----- ------------
c 0.3 1.0 1001
.---------- --__----_
1 4944 4944 0
2 9600 6588 5710 1932 12283 13558 3483
3 9600 6588 270 0 6588 9972 6118 1762
4 9600 6032 0 0 6032 12715 6032 752
._
5 9600 3656 0 0 3656 8410 4916 275
6 9600 2047 0 0 2047 5080 3307 205
7 9600 2011 0 0 2011 3519 3271 194
8 9600 1577 0 0 1577 2752 2837 191
9 9600 997 0 0 997 2147 2257 189
10 9600 396 0 0 396 1701 1656 189
11 9600 0 0 0 0 1386 812 189
12 9600 0 0 0 0 1170 171 i 89
13 6704 0 0 896 0 1170 171 189
14 7615 0 0 1069 0 1170 171 189
15 6626 0 0 1089 0 1170 171 ias
16 5437 0 0 1089 0 1170 171 189
17 4848 0 0 1089 0 1170 171 189
~~~
subject to the c a p a c i t y c o n s t r a i n t s
V3 5 1730 m3
R3 5 9600 rn’
V12 5 14268 m’
R1 5 10000 m3
P2 5 0.70 m’/s
PR3 5 3.66 m’/s
PKA 5 2.00 m’/s
PRI 5 8.20 m’/s
PO1 5 5.00 m’/s
The flow time from the inflow site 12 to the downstream pump was
taken as one time step (i.e. 30 m i n ) and the flow time between the two
pumping stations as two time steps.
The problem was solved w i t h s t a n d a r d software. A typical result is
presented i n Table 12.1 f o r a 210 m i n storm and inflow forecasts of one
time step o n l y (0 to 30 m i n from a c t u a l time). The table includes u n i t
costs c of the o b j e c t i v e function. Sensitivity a n a l y s e s showed that these
c o u l d be specified q u i t e a r b i t r a r i l y , provided that u n i t costs o f different
o r d e r s of magnitude a r e a l l o c a t e d to objectives of d i f f e r e n t p r i o r i t y .
I S E U E R PlCIINlEWNCE R M R D S
1 S E U L R OLOCKAGC R C C O R D S
R L P U E S I E D S T A C T DATE:- 87ClUS
*L:UC:TED CND D A T i :- 1175107
EC’PLAIhTS
--
Z O n P L A I h T : R E F i R R E O TO b A T E R S R A N C M
R L F E R R L O 1 0 R O A D S AND YORKS
0
D
I S E Y i l l BLOCKAGE R C C O R D S
R L C U E S l r O S l 8 R T DATE:- 117OlOf
D L i U C S T L D E k J C A T i :- 6 7 0 1 0 7
l E O U E S T E D S T l R l DATE:- 070105
SEOUESTED END DATE :- 870107
€AM: 1 SIZE: 4
CHAPTER 13
1. Sample Collection
2. Laboratory Analysis
3. Data H a n d l i n g
4. Data A n a l y s i s
5. I nformation U t i I i z a t i o n
3. Sampling Frequency
(J 2u
- Y
LY - 2oy (13.1)
SAMPLING FREQUENCY
n = [ aizS] (13.3)
0.9
0.8
0.4
0.3
0.2
I 1 I I I
0. I
10 20 30 40 50
Number of Somples per Yeor
D ISCUSS ION
The above section has pointed out many problems due to not d e s i g n i n g
a m o n i t o r i n g system in a systems context. Perhaps the major concern is
that all aspects of a monitoring program should match in terms of
accuracy. For example, i t would not be wise to use time series analysis
on nonrepresentative, grab sample d a t a . The system would be providing
excessive accuracy i n one segment compared to the accuracy in another
segment .
I n a s i m i l a r manner, i t may be u n r e a l i s t i c to encourage use of more
sophisticated sample collection and laboratory a n a l y s i s techniques if the
d a t a i s not to receive a thorough s t a t i s t i c a l a n a l y s i s .
It i s difficult to test hypotheses, make decisions and i n i t i a t e action
u s i n g water q u a l i t y d a t a which a r e collected o n l y in the daytime a n d not
flow weighted, several times a year, from locations which are not
completely mixed a n d u s i n g l a b analyses procedures which may h a v e more
variation in their results when analyzing the same sample than the
ambiant v a r i a t i o n of the water q u a l i t y v a r i a b l e in the r i v e r .
Perhaps an even l a r g e r concern to those in m o n i t o r i n g network design
i s the use of water quality standards that generally ignore statistics.
T h i s lowers the v a l u e of a n y information from a compliance v i e w p o i n t , to
t h a t of spot checks. I n c o r p o r a t i n g water q u a l i t y means a n d v a r i a t i o n into
standards would greatly facilitate incorporating more statistics into
216
REFERENCES
AUTHOR INDEX
Kemp, P.H. 64
Kim, J.I. 70
Abulnour, A.M. 116
Kleinecke, D. 41
Adarns, B.J. 190
Adarnson, P.T. 76
Lance, J.C. 143
A d r i a n , D.D. 209,212
Larnbert, J.L. 66
A g a r d y , F.J. 66
Larnbourne, J.J. 66
American Water Works Association 37
L a n g e l i e r , W.F. 3, 5, 6
A r n o l d , R.W. 36
L a n y o n , R. 75
L a r s o n , T.J. 104
Baker-Duly, H.L.G. 123
L a u r i a , D.T. 165
B a l l , J.M. 70, 77
L e i g h t o n , J.P. 146, 149
B a r e n b r u g , A.W.T. 2
Lettenmaier, D.P. 212
Bauer, C.S. 143, 146, 149
Lewis, R.W. 119, 149
Beck, M.B. 202
L l o y d , P.J. 1
Bedient, P.B. 66
L o f t i s , J.C. 209, 212
Betz, 3
Loucks, D.P. 116
Bishop, A.B. 165
L u d w i g , L. 9
Boyd, G.B. 66
.
B r a d f o r d , W. J 64
L y n n , W.R. 116
B r e b b i a , C.A. 62
M a d i s h a , J.L. 75
Brownlow, A.H. 1
Mathew, K. 143
Bungaard-Nielsen, M. 210
McDonell, D.M. 56
Chan, W.Y.W. 167
McPherson, D.R. 41, 45
Chiang, C.H. 165
M i c h a i l , M. 143
CIRIA. 164
M i k a l s e n , K.T. 75
C o l w i l l , D.M. 66
Mrost, M. 1
Connor, J.J. 62
MOller, D. 190
Corbetis, S. 116
Cordery, I. 70
Neilsen, K.5. 210
Crabtree, P.R. 167
Neurnann, A. 190
Newrnan, P.W.G. 143
D e i n i n g e r , R.A. 36, 39, 51
D a n t z i g , G.B. 82, 163
O'Conner, B.A. 56
Orlob, G.T. 212
F r i e d , J.J. 55
Fuchs, L. 190 P a l i n g , W.A.J. 141, 143, 145
P e l l e t i e r , R.A. 1
G i l b e r t , R.G. 143 P e r r y , R. 66
Goodier, J.M. 63
Peters, C.J. 66
Green , I. R.A. 64
Petersen, 5.0. 190
G r i z z a r d , T.J. 70
P o l l s , I. 75
Grosman, D.D. 86
Pomeroy, R.D. 212
Porges, J. 2
H a d l e y , G. 162
P r a t i s h t h a n a n d a , S. 165
H a l l , G.C. 160
Helsel, D.R. 70
Quimpo, R.G. 211
Henderson-Sel l e r s , B. 24
H i l t o n , E. 27, 119
R a n d a l l , C.W. 70
Hine, A.E. 197
Rand Water B o a r d 155
Hinton, E. 149
Revelle, C.S. 116
Ho, G.E. 143
Rice, R.C. 143
Hoehn, R.C. 70
Rinaldi, S. 116
Holton, M.C. 75
Ryzner, J.W. 36
Hunter, J.V.I. 66
S h a w , V.A. 167
Shoemaker, C.A. 146, 149
Simpson, D.E. 64
Smeers, Y. 116
Smith, A.A. 119, 149
Soncini-Sessa, R. 127
South A f r i c a n Bureau of S t a n d a r d s 72
S p r i n g e r , N.K. 66
Steel, T.D. 209, 212
Stehfest, H. 127
Stephenson, D. 27, 66, 80, 81, 82,
115, 116, 117, 163, 175,
197, 200
T e r s t r i e p , 66
Thomann, R.V. 39
Timoshenko, 5 . 55
Tyteca, D. 116
U h l i g , H.H. 13
Van Staden, C.M.V.H. 2
Velz, C.J. 41
Zukovs, G. 190
21 9
SUBJECT INDEX
C o r r e l a t i o n 66
Acid 1 Corrosion 3, 13
Additives 6 Cost 79, 107, 146
Advection 21, 52 C r i t e r i a 211
Aerobic 9 Crop 17
A g r i c u l t u r e 17 Crump w e i r 65
Antecedent moisture 66 Crystal 7
Air 1 C y a n i d e 16
Alkalinity 3 Cycle 214
A l l o c a t i o n 79
A l l o y 13 D a t a 177, 204
Ammonia 9 Dead water 24
Anaerobic 9 Decomposition p r i n c i p l e 163
Analyses 195 D e s a l i n a t i o n 99, 115
A n a l y t i c a l 39 D e t e r i o r a t i o n 116
Apartments 167 D i f f u s i o n 36
A q u i f e r 141 Disc 128
Arsenic 17 Dispersants 7
A r t i f i c i a l r e c h a r g e 141 Dispersion 21, 166
D i s t i l l a t i o n 101
Backwater 193 DO ( d i s s o l v e d o x y g e n ) 37, 206
B a c t e r i a 9, 206, 207 Dissolved s o l i d s 206
B a r i u m 16 Downstream 193
B a s i n 209 D r y d a y s 66
Benefits 126 D r y weather 77
Bicarbonate 67
Biocide 9 Economics 99
Blend 89 E l e c t r i c a l corrosion 14
Blowdown 2 E Iect r o d ia I v s i s 105
BOD (biochemical oxygen demand) 37 Emulsion 10
Booster 152 Env ironmenta I 193
Bottleneck 173
Boundaries 62 Equipment 107
Bremen 193 E r r o r 91
B r i n e 104, 122 E s t u a r i e s 37
E u l e r 57, 59
Calcium carbonate 4 Evaporation 2
C a l i b r a t i o n 40 E x p l i c i t 39, 51
C a p i t a l 107, 157
Carbonaceous 38 F a l l o u t 66
Cathode 10 F a r a d a y s l a w 14
Catchment 64 Feedback 205
Cellulose acetate 104 F i e l d 45
C h a r a c t e r i s t i c 39 F i n i t e d i f f e r e n c e 55
Chelant 7 F i n i t e elements 62
Chemical 67 F i r s t f l u s h 70
C h l o r i d e 2, 67 F l o o d i n g 193
Chlorine 9 Flow 166
C i v i l e n g i n e e r i n g 107 Foam 8
Commerci a I 1 70 F o r m u l a t i o n 88
Cleaning 115 Fouling 9
Computer 20, 115, 128 Four p o i n t 51
Concentration 2, 71, 159, 212 F o u r i e r series 54, 169
C o n d u c t i v i t y 18 Freezing 103
Conduit 175 Frequency 21 1
Confidence 214
C o n s t r a i n t s 41, 86 Gain 23
Conveyance 141 G a l v a n i c corrosion 13
Cooling 20 Geochemical 1
Geohydrology 41
220
Tape 128
Taste 16
T a y l o r series 53
TDS ( t o t a l d i s s o l v e d s o l i d s ) 2, 95
Temperature 3, 206, 107
Terminal concentration 24
Time l a g 166
Topography 1%
Toxic 16
T r a f f i c 69
T r a n s p o r t a t i o n p r o g r a m m i n g 80
Treatment 141, 155, 157
T u r b i d i t y 206
Turbulence 8
Two step 39, 52
U n p r e d i c t a b l e 64
Upstream 193
V a a l r i v e r 155
Vapour compression 102
Vegetables 18
Ventilation 2
This Page Intentionally Left Blank