Guidelines For Design of Small Hydropower Plants

TECHNICAL CARD
Copyright c 2000.
Published in 2000 through WREAN (Western Regional Energy Agency &
Network) and DED (Department of Economic Development), Belfast, North
Ireland.
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 except under terms
of the Copyright, Designs and Patents Act 1988 without the permission in writing
of the publisher, with the exception of any material supplied specifically for the
purpose of being entered and executed on a computer system for exclusive use by
the purchaser of the publication.
Title:
Guideline for Design of
SMALL HYDROPOWER PLANTS
Helena Ramos
A. Betmio de Almeida
M. Manuela Portela
H. Pires de Almeida
(Chapters: 1, 2, 3, 5, 6, 7 and 9)
(Chapters 1, 2, 3, 5, 7, 9 and reviser)
(Chapters: 4 and 10)
(Chapter 8)
A. Carvalho Quintela
(Foreword and reviser)
Legal deposit: July 1999.

Contents production: Helena Ramos.
Cover production - Designer: Joana Cruz.
This book has the support of CEHIDRO, WREAN and DED.
Typeset in 11/14 pt Times New Roman MS Word 97.
FOREWORD
The book Guidelines for design of SMALL HYDROPOWER PLANTS

comprises the following ten chapters covering the main subjects connected with
the design of such plants:
1 Hydroenergy.
2 General types of small hydropower plants.
3 Hydropower and water uses.
4 Hydrology.
5 Hydraulic design of small power plants.
6 Small hydraulic turbines.
7 Hydraulic transients and dynamic effects.
8 Electrical equipment.
9 Environment.
10 Economic analysis.
The first three chapters contain an introduction to the problems involved in the
subject under consideration.
Chapter 4 deals with hydrologic studies for small hydroelectric schemes with a
run-of-the-river exploitation in order to define the design discharges of water
conveyance systems and of hydraulic turbines, and power productions. The
proposed methodologies may overcome the non-existence of basic flow
information at the river section of the scheme intakes.
Chapter 5 contains relevant information about the hydraulic design of intakes,

bottom outlets, hydraulic conveyance systems (including canals and penstocks)
and powerhouses.
The object of Chapter 6 is to characterise the differences among the types of
small turbines and their behaviour, as well as turbines application ranges and to
estimate main dimensions of the turbines.
Chapter 7 is the longest one (it fulfils with chapters 5 and 6 half of the book) and
presents the main problems concerning hydraulic transients and dynamic effects.
Chapter 8 deals with synchronous and asynchronous generators, electrical
facilities and control systems.
Environment impact assessment is the object of chapter 9.
Chapter 10, after considering capital costs, annual operations costs and benefits,
presents an economic analysis according to two different concepts and defines
economic indexes. It includes an application example of economic analysis.
The chapters authors are as following:
Helena Ramos
- author of chapters 1, 2, 3, 5, 9 (being second author A. Betmio de
Almeida), and 6, and second author of chapter 7 (being first author A.
Betmio de Almeida).
- author of chapter 7 (being second co-author Helena Ramos) and second
author of chapters 1, 2, 3, 5 and 9 (being Helena Ramos the first author).
M. Manuela Portela
- author of chapters 4 and 10.
Pires de Almeida.
- author of chapter 8.
Each author has a good professional experience as consulting engineer in the

field of his participation in this book. The first three mentioned authors also have
a significant teaching and research experience in such field. The very special
background of Prof. A. Betmio de Almeida must be emphasised.
The book identifies the more important subjects of interest for the design of
small hydropower plants, which are treated under a rigorous, synthetic and
understandable way.
The text is completed by appropriated illustrations.
Thus, I think that the book will be a very useful tool for those who have a
participation in one or more of the multi-disciplinary aspects involved in the
design of small power plants and a favourable acceptance will be expected.
Antnio de Carvalho Quintela
ACKNOWLEDGEMENTS
This book is the result of the applied experience in several small hydropower
base and final type projects, since I joined the HIDRO4 (now a part of
SOMAGUE Engineering Group). A research program in this domain was
developed, since the eighties, at Instituto Superior Tcnico (IST) of Technical
University of Lisbon, and was the basis of my PhD. The book also benefited
from the contribution and support from different parts that I would like to
acknowledge hereafter.
Firstly, I would like to give special thanks to my adviser Professor A. Betmio
de Almeida (Full Professor of IST and an expert in hydraulic transients), whose
enthusiasm and scientific revision contributed substantially with their
suggestions. Moreover, the colleague Prof. M. Manuela Portela for chapters 4
and 9 and Eng. H. Pires de Almeida, Manager of ENERSIS a mini-hydro power
company, with his contribution for chapter 7, are greatly appreciated.
I also want to express my gratitude to Professor A. Carvalho Quintela (Full
Professor of IST and a very experienced designer of hydroelectric schemes) for
taking the time out to do relevant remarks to the original and to write the
foreword of this book.
Special thanks go to Dr. Robert Gibson, manager from WREAN (Western

Regional Energy Agency and Network, North Ireland) and Dr. Lawrence Foy
from DED (Department of Economic Development Energy Division) through
their invitation, confidence and interest always proved along the development of
this book.
I am thankful for the co-operation of several companies, such as ALSTHOM,

BOUVIER, CEHIDRO, DED - Energy Division, DUMONT, ENERSIS, EDP,
HIDRO4/SOMAGUE, HIDROCENTRAIS REUNIDAS, LNEC, NEYRPIC,
VOITH and WREAN, which allowed to support the publishing of this document
including some knowledge, photos, project elements and data information.
Helena Ramos
PhD in Hydraulic/Civil Engineering
I S T - Technical University of Lisbon
[email protected]
CONTENT
1- HYDROENERGY
1.1- Renewable energy concept

1.2- Hydroelectricity: energy conversion and hydropower principles
1.3- Environmental considerations
1
3
7
2- GENERAL TYPES OF SMALL

HYDROPOWER SCHEMES
2.1- Classification of hydropower plants

2.2- Hydropower schemes
2.3- Site location
9
10
13
3- HYDROPOWER
AND WATER USES
15
3.1- Power generation and different sector-users

3.2- The benefit of small hydropower plants
15
19
4- HYDROLOGY
21
4.1- Introduction. Scope of the studies

4.2- Basic data required for the hydrologic study
4.3- Mean annual And mean daily flow series
21
22
25
4.3.1- Introduction
4.3.2- Annual precipitation and annual flow. Mean
daily flow series
4.3.3- Mean annual flow duration curve
4.4- Energy Evaluation
25
26
31
34
5- HYDRAULIC DESIGN OF
SMALL POWER PLANTS
39
5.1- Introduction
5.2- Layout schemes
5.2.1- Intakes
5.2.1.2- Minimum submergence
5.2.1.3- Bottom intake
5.2.2- Protection rack
5.2.3- Sedimentation or desilting basin
5.2.4- Weir
5.3- Sluice bottom outlet
5.4- Conveyance system
5.4.1- General layout
5.4.2- Head losses and net head
5.4.3- Canals
5.4.3.1- Uniform and steady state hydraulic regimes
5.4.3.2- Boundary conditions
5.4.3.3- Forebays
5.4.4- General remarks about mixed circuit
39
40
40
40
44
48
51
54
57
58
58
60
66
67
69
69
71
5.4.5- Penstocks
5.5- Powerhouses
5.6- Analysis of hydropower schemes
72
75
78
6- SMALL HYDRAULIC
TURBINES
79
6.1- Types of turbines

6.2- Turbine similarity laws and specific speed (Ns)
6.3- Turbine efficiency
79
84
89
6.4- Dimensions of turbines
90
7- HYDRAULIC TRANSIENTS
AND DYNAMIC EFFECTS
101
7.1- Introduction
7.2- Canal systems
7.3- Pressurised systems
7.3.1- Typical transient regimes
7.3.2- Preliminary analysis
7.3.3- Governing equations
7.4- Overspeed dynamic effects
7.4.1- Overspeed runner control
7.4.2- Overspeed effect on turbine discharge
7.4.3- Turbine overspeed effects on waterhammer
7.5- Special protection devices
7.5.1- Introduction
7.5.2- Surge tanks
7.5.3- Differential surge tank
7.5.4- Air vessel
7.5.5- Synchronised valve or relief valve
7.5.6- Flywheel
7.5.7- Protection devices behaviour
7.5.7.1- Analysis of a surge tank
7.5.7.2- Analysis of an air vessel
7.6- Examples
7.7- Other protection devices and procedures
7.8- Integrated analysis and design
7.9- Case studies
101
104
108
108
110
113
115
115
119
122
123
123
124
125
127
129
130
131
131
136
138
140
142
143
8- ELECTRICAL EQUIPMENT
145
8.1- Generators
8.1.1- Synchronous generators
8.1.2- Asynchronous generators
8.2- Electrical installations
8.2.1- Main transformer
8.2.2- Switchgear
8.2.3- Control equipment
8.2.4- Electrical protection
145
145
148
149
149
150
150
150
8.3- Control system considerations

8.3.1- Introduction
8.3.2- Speed control
8.3.3- Water level regulation
8.3.4- Generator output power control
8.3.5- Voltage and power factor regulation
8.3.6- Flood level control
8.3.7- Blade control of Kaplan turbines
8.3.8- Analytical representation
8.3.9- Switchyard
151
151
152
152
153
153
153
153
154
154
9- ENVIRONMENT
155
9.1- Environmental impact assessment

9.2- Identification of impacts
9.3- Activities associated with hydropower plants
9.4- Reference criteria to support environmental studies
9.5- Safety
155
157
161
163
163
10- ECONOMIC ANALYSIS
167
10.1- Introduction. Costs and benefits

10.2- Economic analysis
10.2.1- Introduction
10.2.2- Constant market prices system concept
10.2.3- Discount rate and present value concept
10.2.4- Economic indexes
10.2.4.1- Basic considerations
10.2.4.2- Net Present Value (NPV)
10.2.4.3- Benefit/cost ratio (B/C)
10.2.4.4- Internal Rate of Return (IRR)
10.2.4.5- Average Price of the kWh (AP)
10.2.4.6- Payback period (T)
10.2.5- Sensitivity analysis
10.2.6- Application example
167
172
172
172
173
181
181
184
184
185
187
188
188
188
BIBLIOGRAPHY
191
SYMBOLOGY
A
a
AF
Af
AGS
AHP
AP
ARAC
As
B/C
Bf
BH
BR
Bs
b
bo
C
flow cross section (m2).

net interval between two bars of a rack of a rack (m).
actualisation factor (-).
horizontal area of the forebay (m2).
area of the watershed of the stream-gauging station (km2).
area of the watershed of the hydropower scheme (km2).
average price of the kWh (cost unit).
cross-section of the air vessel (m2).
superficial area (m2).
benefit/cost ratio (-).
forebay width (m).
hydraulic torque (N.m).
resistant or electrical torque (N.m).
superficial width (m).
length bar of a rack in the flow direction (m).
inlet runner height (m).
discharge coefficient (-); runner speed (m/s); parameter of flow
approximation (-); Hazen-Williams coefficient (m0.37/s); present value
of the capital costs (cost unit).
c
Cd
CE
CEC
Ci
CMP
Cp
Cr
C vd
Cvalv
D
d
thickness or diameter of a rack bar (m); wave celerity in open channel

or in a pressure pipe (m/s); specific runner speed (-).
discharge coefficient of a rack (-).
equivalent tariff of energy (cost unit).
coefficient of energy cost (-).
capital costs in year i (cost unit).
pipe material coefficient (-).
pipe material cost per unit weight (cost unit).
ratio area between free and total or gross rack area (-).
weir discharge coefficient (-).
valve discharge coefficient (-).
pipe diameter (m); average number of days per year during which
mean daily flows are equal or greater than a given flow (days/year).
solid particle density (ratio between specific weight of a particle p and
of the water , d =
D1; D2
D50
DN
Do
E
Eb
E
f
fb
Fr0
Fr1
fs
g
H
h
Ha
Hd
HF
hG
p
(-); opening of the intake (m).
inlet and outlet diameter of a turbine runner (m).

reference dimension of a settling solid particle (m).
nozzle diameter (m).
optimum diameter or economic diameter (m).
Young's modulus of elasticity of the conduit walls (N/m2); Euler
number (-).
modulus of elasticity of rack bars (N/m2).
mean annual energy production (GWh).
Darcy-Weisbach friction factor (-); frequency of the electric grid (Hz).
rack bar oscillation frequency (Hz).
Froude number of the flow (-).
Froude number related to the wave propagation of a bore (-).
swirl frequency (Hz).
gravitational acceleration (m/s2).
annual flow depth (mm); specific energy flow (m); head (m); head over
the rack (m); piezometric head (m).
water depth (canal) (m).
atmospheric piezometric head (m).
downstream water level or outlet level added by kinetic head (m).
final head (energy) (m).
depth of the gravity centre (m).
Hg
HGSi
H g ,c
HHPi
Ho
H RAC
HP
Hr
Hs
I
i
Ii
IRR
Jch
K
k
K ch
K or
K RAC
ku
Ku1, Ku2
K
L
L desc
Lf
available gravity potential head or gross head (m).

annual flow depth over the watershed of the stream-gauging station in
year i (mm).
piezometric head in the junction section of the surge tank (m).
annual flow depth over the watershed of the hydropower scheme in
year i (mm).
neat head (m); head over the spillway (m).
piezometric head in the junction section of the air vessel (m).
piezometric head at any conveyance system section (m).
upstream level variation (m).
freeboard between the Nriver and the nozzle axe for turbine with vertical
shaft, and between water level of the river or rejection canal and the
lowest runner point for turbines with horizontal shaft (m).
rotating mass inertia (kg m2).
time step index (-).
inflow to the scheme water intake during time step i (m3 /s).
internal rate of return (%).
angle between runner blades and rotational speed C direction (degree).
angle between the velocity vector (V) and the rotational speed (C)
(degree).
hydraulic headloss gradient (-).
hydraulic grade line of the interior pipe of the differential surge tank (-).
fluid bulk modulus of elasticity (N/m2); pipe Gauckler-ManningStrickler coefficient (m1/3/s).
absolute pipe roughness (m); trash cleaning coefficient (-); index day
(-); overspeed coefficient (-).
singular head loss coefficient in junction section of the surge tank (-).
singular orifice head loss coefficient (-).
inlet or outlet singular head loss coefficient of the air vessel (-).
peripheral speed coefficient (-).
numeric factors to unit conversion (-).
coefficient that depends upon the angle of the conical needle end (-).
distance between rack bar supports (m); pipe length (m); length of the
hydraulic circuit (m); weir length (m); length of the sedimentation
basin (m); field inductance (-).
perimeter of the internal weir in a differential surge tank (m).
forebay length (m).
N
Nf
nGS
nHS
Nnoz
no
NPV
Nres
Nriver
Ns
O
Oi
Oj
P
p
p atm
PGSi
PGS
PHPi
PHP
Pm
Po
Pt
PV
rotational speed (r.p.m.); period with simultaneous annual records of

precipitation depths and of flows depths over the watershed of the
stream-gauging station (-); time step index (-); polytropic coefficient
for the air behaviour inside a vessel (-); number of the project lifetime
periods (generally measured in years).
number of nozzles of a Pelton turbine(-).
free surface level in a forebay (m).
period with records either of mean daily flows at the stream-gauging
station or of annual precipitation depths in rain gauges located inside
or near the station watershed (years).
period with records of annual precipitation depths in rain gauges
located inside or near the hydropower scheme watershed (years).
nozzle level (m).
nominal rotational speed (r.p.m.).
net present value (cost unit).
reservoir water level (m).
water level at downstream powerhouse (m).
specific speed (r.p.m.).
present value of the annual operation costs (cost unit).
discharge over the spillway during time step i (m3/s).
annual operation cost for year j (cost unit).
annual precipitation depth (mm); installed power (kW); present value
of the reposition costs (cost unit).
number of alternator pairs of polos (-).
atmospheric local head (m).
precipitation depth over the watershed of the stream-gauging station in
year i (mm).
mean annual precipitation depth in the watershed of the
stream-gauging station (mm).
precipitation depth over the watershed of the hydropower scheme in
year i (mm).
mean annual precipitation depth in the watershed of the hydropower
scheme (mm).
reposition cost in year m (cost unit).
hydraulic linear growing power (kW).
installed turbine power (kW).
present value (cost unit).
Q
q
Qd
qd
i
q GS k
turbine discharge (m3/s); mean daily flow (m3/s); spillway discharge

(m3/s).
unit discharge over an inlet bottom rack (m2/s).
design discharge (m3/s).
unit discharge trough the rack (m2/s).
mean daily flow measured in the stream-gauging station during day k
of the year i (m3/s).
q HP k
Qin
Qinlet
Qmax
Qmin
Qmod
QmodHP
Qnoz
Qout
Qoutlet
QP
Q RAC
Qres
Qresi
QRW
QT+V
R
R
r
Re
RH
Rj
S
s
St
estimate, for the water intake of the hydropower scheme, of the mean
daily flow during day k of year i (m3/s).
forebay inflow (m3/s).
inlet discharge (m3/s).
design discharge of the hydropower scheme (m3/s).
minimum discharge compatible with turbine operation (m3/s).
modulus discharge or average of the mean daily flows (m3/s).
modulus discharge at the intake of the hydropower scheme (m3/s).
maximum discharge of each nozzle (m3/s).
forebay outflow, over the weir and for penstock (m3/s).
outlet discharge from the reservoir (m3/s).
discharge at any conveyance system section (m3/s);
3
discharge into or out of the air vessel (m /s).
reserved discharge (m3/s).
reserved discharge during time step i (m3/s).
discharge at runaway conditions (m3/s).
turbine discharge and the valve discharge (m3/s).
present value of the revenues (cost unit).
mean annual benefit (cost unit).
radial distance (m); discount rate (%).
VD
Reynolds number (for sedimentation: R e = s 50 ; for pressurised
UD
flow R e =
) (-).
hydraulic radius (m).

revenues in year j (cost unit).
flow cross section of the nozzle (m2); cross-section of the rack (m2);
submergence (m); elastic wave transmission (%).
bottom slope (-).
Strouhal number (-).
T
t
TC
TE
ti
Tm
To
t0
Ts
Tscour
payback period (years).

time (s); pipe wall thickness (m).
guide vane closing discharge time (s).
elastic time constant (s).
instant i t after instant t0 (s).
machine starting time (s).
average annual time operation (h).
initial instant of numerical simulation (s).
sedimentation time (s).
scour time of a particle (s).
tv
vapour head (m).
TW
U
hydraulic inertia time constant (s).

mean velocity of the flow (m/s).
v
Vo
Vs
V 2 2g
Zd
Zu
z ch
z RAC
unit absolute flow velocity (-).

flow velocity (m/s).
sedimentation velocity (m/s).
kinetic energy head (m).
reservoir volume (m3).
maximum mean annual volume utilised in the energy production
(hm3).
mean annual flow volume at the stream-gauging station (hm3).
mean annual flow volumes at the water intake of the hydropower
scheme (hm3).
air volume inside the air vessel (m3).
relative flow (or meridian) velocity in the turbine runner (m/s).
unit relative velocity (-).
distance along the canal bottom or the pipe axis (m); flow and rack
direction (-).
water level at the tailrace (m).
water level at the intake (m).
water level in the surge tank (m).
water level in the air vessel (m).
angle with horizontal (degree); Coriolis coefficient (-).

coefficient of the rack bar (-).
GS
HP
RAC
W
w
x
absolute flow velocity in the turbine runner (m/s); flow velocity (m/s).
b
p
w
i
Ei
H
h
Hf
hi
HM
HT
p
S
t
i
relative roughness (-).

angular rotating speed (rad/s).
kinematic viscosity of the water (m/s2 ) ( =10-6 m/s2).
specific mass of the fluid (kg/m3).
specific weight of fluid (N/m3) (e.g. for water = 9800 N/m3).
specific weight (N/m3) of the bar material (e.g for steel b=78000
N/m3).
specific weight of the pipe material (N/m3).
specific weight of the water (N/m3).
turbine-generator efficiency (-).
efficiency of the powerhouse during time step i (-).
energy production during time step i (GWh).
head loss (m).
constant net head (m).
pipe or canal friction loss (m).
constant net head during time step i (m).
maximum transient overpressure (m).
total head losses in the hydraulic circuit (m).
pressure difference between two sections (N/m2).
variation of the volume of water stored in the reservoir during a time
step of the simulation program (hm3).
time step of the simulation program (s).
turbined volume during time step i (hm3).
output power angular non-phase (degree).
open degree of the nozzle (degree); angle of the canal with horizontal
(degree).
non-dimensional parameter that depends upon the elastic properties of
the conduit (-).
singular loss coefficient that depends on the geometry of the singularity
(-).
dynamic depression coefficient or Thoma coefficient (-); allowable
stress of pipe wall material (N/cm2).
Subscripts
1 and 2 upstream and downstream of wave front in a canal; inlet and outlet of a
runner.
m; p
model and prototype.
Hydroenergy
HYDROENERGY
Helena Ramos
1.1- Renewable energy concept

The renewable energy concept is basically associated to the following remarks:
Inexhaustible energetic sources, in spite of being limited or

conditioned;
Low polluted energy with small environmental impacts;
Relevant component of a sustainable development.
Nowadays, the policy in most of the countries is devoted to assure additional

generating energy from renewable, in particular with small hydropower
schemes, which can contribute with a cheap source, as well as to encourage
internationally competitive small industries across a wide range of new energy
sources options and technologies.
The hydraulic power is one of the oldest energy sources of the mankind, namely
for irrigation and industry. Nowadays, small hydro is one of the most valuable
answers to the question of how to offer to isolated rural communities the benefits
of electrification and the progress associated with it, as well as to improve the
quality of life. The hydroelectric power plant utilises a natural or artificial fall of
-1-
Guidelines for design of SMALL HYDROPOWER PLANTS
a river. The water flow energy is used to turn the wheel of a turbine and returns
again to the river. This type of electricity production does not consume water,
thus it is usually considered a renewable energy source. The flow will continue
to fall downhill and the water will continue to be available as a resource for men
and environment needs, thanks to the natural hydrologic cycle.
The economic utilisation of renewable energies is now based on new
technologies and on environmental protection techniques. Small hydropower,
with its multiple advantages, as a decentralised, low-cost and reliable form of
energy, is in the forefront of many developing countries to achieve energy selfsufficiency.
Fig. 1.1 - Typical scheme of a renewable energy source based on the waterpower.
For environmental protection it must be considered, in each small hydro project,

the ecological or reserved flow in order to protect downstream the wildlife
habitats and to encourage or maintain the migration through fish-passages.
It will be enhanced the main advantages to develop small hydro comparing with
other electricity sources:
It saves consumption of fossil, fuel, and firewood.

It is self-sufficient without the need of fuel importation.
It does not contribute for environment damages by resettlement, as it
occurs with large dams and reservoirs.
It can be a good private capital investment in developing or developed
countries.
-2-
Hydroenergy
It offers a decentralised electrification at a low running cost and with

long life.
A small-scale project can also induce tourist activities and can benefit both rural
and small urban areas with a friendly water scenario.
1.2- Hydroelectricity:
energy conversion and hydropower principles
The theory of hydropower generation is based on the conversion of the hydraulic
potential energy of a flow into electric energy, which corresponds to a
differential net head. The energy of the flow is associated to the gravity energy
through natural or artificially created topographic water falls in rivers or through
hydraulic conveyance systems, composed by pressurised pipes or penstocks or
by mixed hydraulic conveyance system composed by canal and penstocks.
According to the principle of conservation of energy, the energy balance of a
steady flow from A to B will obey to the following relationship:
ZA +
p A A U 2A
p
U2
+
= Z B + B + B B + H AB
g
2g
g
2g
HA
(1.1)
HB
where ZA and ZB (m) are the elevations above a datum plane, pA and pB (Pa) are
the pressures at the centres of gravity of the flow cross sections at A and B, UA
and UB (m/s) are the average flow velocities respectively at A and B, (kg/m3)
is the water density, g (m/s2) is the gravity acceleration and A and B are
numerical coefficients accounting for the non-uniform velocity distribution.
Equation (1.1) express that the difference between total heads at A (HA) and B
(HB) equals the headloss HAB between the two flow cross sections, where the
head is the total flow energy by the weight of the flowing water.
For free surface flow, equation (1.1) simplifies to the following form:
H AB = N A N B + A U 2A B U 2B / 2g
-3-
(1.2)
where NA and NB (m) are respectively the elevations of the free surface of cross
sections at A and B. Equation (1.2) express that the dissipated head equals the
head difference between A and B. Should the difference A U 2A B U 2B be very
small or equal and the dissipated head equals the difference between elevations
NA and NB.
NA
HAB = Hg
NB
A
B
NA
Headwater Dam
HCB Head loss
Ho
NB
A
C
or
Net head to
be converted
B Tailwater
Powerplant
HC
Ho
Diversion works
Total head
upstream the
powerplant
Fig. 1.2 Different types of hydropower schemes.
The basic hydropower principle is based on the conversion of Ho or net head, the
large part of the naturally dissipated head along a watercourse into mechanical
and electrical energy (Figure 1.2):
H o = H g H CB
(1.3)
-4-
Hydroenergy
implying that the headloss HCB along the hydraulic conveyance circuit of the
scheme will be much less than Hg. The headloss between A and C are artificially
reduced in order to convert Ho into electricity.
The net head of a small hydropower plant built at section B can be created in
quite a number of ways. Two fundamental ways are the following ones:
to built a dam across a stream to increase the water level just above the
powerplant;
to divert part of the stream, with a minimum of headloss, to just above
the powerplant built far away the dam.
The net flow power Po and the corresponding energy Eo over an interval time t
(s) of the hydropower plant will be respectively:
Po = g Q H o
(1.4)
E o = g Q H o t
where Q (m3/s) is the constant discharge diverted to the powerplant.
The final useful head delivered to the electrical network is smaller than the
available gross head:
HF =H g
(1.5)
where is the global efficiency, resulting of the multiplication of partial

efficiencies from the successive transport and conversion phases ( < 1 ).
Along the same watercourse the total gross head can be profited by a multistage
scheme with several powerplants (i.e. cascade system).
The hydroelectricity production is an energy conversion process, in which the
water is the vehicle of transmission and transformation of the gravity potential
energy into mechanical and electric energy by the turbine-generator set installed
in the powerhouse. The water is led through pipes and/or canals to the turbine,
which turn the shaft of the generator to produce electric energy. From the
powerhouse, after a convenient voltage transformation, transmission lines carry
electricity out to communities or to the national grid.
-5-
In a river, the available potential energy or gross head (Hg) will be converted in a
system through the following components:
Reservoir: constitute a storage form of the available potential energy and

creates the conditions for water diversion through the intake.
Conveyance system, including the intake, conveyance canal, penstock,
galleries and tailrace or outlet where part of the available energy is
converted into kinetic energy; another part is transformed into reversible
flow work capacity (pressure head) and another part is dissipated in heat
(by fluid viscosity) resulting in the net or useful head.
Hydraulic turbine: where the net head is converted into rotor speed of the
turbomachine.
Generator rotor: the mechanical energy on the shaft maintain the speed of
the rotor and it is transformed in electric energy according to
electromagnetic laws.
Line to link to the grid: the electric energy is driven and transformed in
order to connect to the grid for transportation to long distances and
distribution.
Powerplant
Hg
Tailwater
Turbine
Penstock
Diversion canal
from the river
Generator
Fig.1.3 Small old hydroelectric scheme

(adapted from MACINTYRE, 1983).
The hydraulic engineering intervention has a major contribution in these types of

projects, namely in the planning, conception, study, design, building and
-6-
Hydroenergy
exploitation phases of a project that involves multidisciplinary teams, with

technicians and experts in several domains: civil, mechanical, electric
engineering and specialists in geologic and environmental sciences, among
others.
Adequate head and flow are necessary requirements for hydro generation. For
the head characterisation is necessary to consult available maps and confirm in
situ trough field visits and surveys about the potential sites. Typically, maps
information need to be complemented by field surveying along the hydraulic
circuit. Such maps are particular important for positioning the structural
components of the system. In what concerns the flow, a hydrologic study must
be carried out. After a first approach to power potential calculation, estimation
of energy output, identification of civil works and other critical issues (e.g.
environmental and social constraints), a technical-economic feasibility study
must be performed.
In small hydroelectricity, the hydraulic structures will be much less complex
than in large hydroelectricity. The hydraulic conveyance circuit can be
integrated in other components for multiple purposes (e.g. irrigation or water
supply schemes). Sometimes, small hydropower plants can be very
unconventional both in design concept and in components (e.g. the turbines can
replace pressure reducing valves or other types of localised dissipation of excess
head). The small hydroelectricity depends mainly upon the local and regional
characteristics. It has low environmental impacts but also a relatively less
guaranteed energy production due to the very small storage volume in the
upstream reservoir.
1.3- Environmental considerations

Hydropower plants produce no carbon dioxide, sulphur oxides or nitrous oxides,
no air emissions and no solid or liquid wastes. Nevertheless there are impacts by
retaining water and inducing sediment to settle down, as well as by obstructing
the fish passage and upsetting the wildlife habitats.
The world concern about global planet warming phenomenon essentially due to
CO2, SO2, NOx emissions in energy generating process with fossil fuels and the
problem for the future of nuclear wastes, more and more will be emphasised the
advantages of energy production trough renewable sources. Hydropower
represents an important environmental benefit to aid sustainable development
-7-
because there is no release of carbon dioxide that contributes to ozone depletion

and global warming. In the future a penalty cost by combustion of fossil sources
should be applied. Being hydroelectric energy 11.7% from total generated
energy in European Union, the small hydropower plants can still occupy a
relevant contribution around the world, being estimated that the energy
production through renewable sources can increase three more times than the
values obtained in ninety decade.
However the abstraction of water from a watercourse must be controlled in order
to avoid serious damages to aquatic biota between the intake and the tailrace. In
fact, a variable residual flow should remain along the year between the diversion
dam and the powerplant. This residual flow, also known by reserved flow,
compensation or ecological flow, must be environmentally acceptable. In
particular, in seasons of low flow the residual flow is very important in order to
keep a steady regime to warrant the aquatic natural development and the water
quality. In this way, a compromise has therefore to be sought to ensure the
maximum energy production with maintenance of the equilibrium of the aquatic
system.
The small hydropower does not require high dams because the majority of them
are run-of-river schemes, meaning simply that the turbine only generates when
there is available water. A minimum daily storage and flow regulation is
typically guaranteed. When the river dries up the generation ceases.
-8-
____________________________________General Types of Small Hydropower Plants
GENERAL TYPES
OF SMALL
HYDROPOWER PLANTS
Helena Ramos
2.1- Classification of hydropower plants

The classification of hydropower plants can be based on different factors:
head: low (less than 50 m); medium (between 50 and 250 m); high
(greater than 250 m);
exploitation and storage: with daily (or seasonal) flow regulation
(reservoir type); without flow regulation (run-of-the-river type);
conveyance system: pressurised (penstock); mixed circuit (canal and
penstock);
powerhouse site: dam or diversion scheme;
energy conversion mode: turbining or reversible pumping-turbining;
type of turbines: impulse, reaction and reversible;
installed power: micro (Pt < 100 kW); mini (100 kW < Pt < 500 kW);
small (500 kW < Pt < 10 MW).
The classification based on the power is very important because is an

institutional and legislate reference (ESHA, 1994).
-9 -
The present concept of small-scale hydropower is neither a miniature of large

hydro nor a simple repetition of old techniques, but it requires advanced and
modern studies and properly adaptations for different solutions.
2.2- Hydropower schemes

In Europe, the development of small hydroelectricity has known a new growing
since the seventies, due to the world energy crisis and concerns of negative
environmental impacts, as well as due to the development of automation and
remote control (i.e. abandoned exploitation) and the standardisation of the
equipment (e.g. for turbines and generators). The small hydropower can fill the
gap of the decentralised production, for instance for private or municipal activity
production for sale to national electric grid or for supply industries, rural or
isolated zones, improving their development.
Small hydropower is typically associated with rivers with catchement areas of
less than 200 km2. The net head and the plant discharge are two important
parameters to be considered in the small hydropower design. The reservoir for
discharge regulation allows, under certain limits, to make independent the
energy production from river flow variations. With an adequate reservoir volume
it is possible to program the power generation in order to be able to satisfy the
demand or to sale energy in rush hours, when its price is valorised. Nevertheless,
a small hydropower plant is rarely compatible with a large reservoir due to
economic constrains. However, it can be shared with other types of water uses
(e.g. irrigation systems, water supply systems, flood protection structures or
discharge regulation structures), and can generate energy whenever excess
discharge exists.
According to the mode of net head characterisation, the following main
hydropower schemes can be identified:
Dam scheme - The dam is used to concentrate the head, which raises the
upstream water level. In this way, the powerhouse can be placed
either at downstream incorporated inside the dam.
- 10 -
Fig. 2.1 Typical scheme of dam used to raise head for intake
diversion and small water storage.
Diversion scheme The utilisation of diversion structures, such as canals,

tunnels or galleries or low-pressure conduits allows the head gain.
A small dam can be used with long hydraulic circuit in order to
obtain the net head at downstream end.
Small dam
Canal
Penstock
Powerhouse
Fig. 2.2 Typical scheme with a mixed hydraulic conveyance

circuit with a canal and a penstock.
Mixed scheme - A dam can partly raise the net head and a long hydraulic
conveyance circuit will raise the other part.
Small dam
Surge tank
Tunnel or
gallery
Penstock
Powerhouse
Fig. 2.3 Typical scheme of a low pressure tunnel

and a penstock as a diversion structure.
-11 -
According to the mode of discharge exploitation, the following hydropower

scheme types can be identified:
Run-of-the-river scheme Power is generated without inflow regulation.
It is a common scheme applied to mini or micro hydropower plants.
Daily regulation scheme Power is generated according to the natural
fluctuation of the daily demand, the water being stored in a
regulating pond or small reservoir at off-peak times and discharging
it at peak hours, resulting a bigger energy output comparatively to
without regulation capability.
Seasonal regulation scheme This scheme is commonly applied in larger
power plants, which needs a reservoir to store water in rainy season
and discharge it in the dry season, enhancing a constant energy all
year (normally it is not a common scheme in small power plants).
Cascade scheme The cascade scheme is a typical exploitation of the
river, in order to make the best use of the river falls.
Fig. 2.4 Example of a cascade scheme of large dams (courtesy of EDP).
In case of existing water intake in differential canal (e.g. irrigation system) it can
be used a siphon, which penstock pass over the dam without affecting it.
Nevertheless, for the turbine start-up it is necessary to use a vacuum pump and
for the stoppage, the system must have installed an air valve. The gross head is
equal to the difference between upstream and downstream water levels
(normally for schemes with heads smaller than 10 m and discharges between 1to
50 m3/s). These solutions can avoid onerous civil works.
- 12 -
2.3- Site location

A hydro generation site is mainly conditioned by head and flow requirements. In
a preliminary stage, the power potential and energy output, the estimation of the
need works, and the economic feasibility, must be studied.
Depending upon the geographical characteristics of the available site, the
hydropower scheme can have high or low head. High head schemes can be less
expensive because, the required hydraulic equipment (e.g. dam, turbines and
valves) will be smaller. Nevertheless, the topographical and geotechnical
characteristics could constrain the type and the alignment for the conveyance
system. The nature of the streambed, bends and the access, stability of the soil,
are important condition factors to take into account in the site choice.
-13 -
- 14 -
Hydropower and Water Uses
HYDROPOWER
AND WATER USES
Helena Ramos
3.1- Power generation and different sector-users

Water resources can be used in different ways to serve the society, taking into
consideration all demands arising from different social and economic sectorusers. Meanwhile, the exploitation and the utilisation of water resources aim to
obtain the maximum benefit that should be controlled in order to reduce natural
hazards and environmental impacts. Small hydropower can be associated with
different water uses:
Power generation and water supply - Water conveyance systems to feed
a water supply of a town through a pressure pipe, from a reservoir
to a treatment plant are, e.g., equipped by a valve system to
dissipate the excess energy. In that case, a turbine can substitute this
dissipation system. However, this solution demands the installation
of a by-pass in order to guarantee the flow continuity in case of
turbine stoppage. The turbo-generator units will work according to
special regulation and control systems under the command of a
water treatment plant. They can be built in canals or on water
supply high pressure pipes (e.g. micro turbines replacing pressurereducing valves to take the advantages of excess head Figure 3.1).
- 15 -
Water supply system

Reservoir
Pressure
reducing valve
Housings
Water distribution
system
Generator
Inlet
control
valve
Daily demand
Turbine
Outlet
Inlet
Fig. 3.1 Example of application a micro-turbine as a pressure reducing

valve in a water supply system to yield power.
Power generation and irrigation It is possible to install a powerplant in

irrigation canal or in a lateral by-pass canal, with a small penstock,
in order to take the advantage of the head created by an upstream
dam. In several irrigation schemes there is the opportunity to insert
a small powerplant and to convert, permanently or seasonally, the
excess head into electric energy:
should an upstream dam imposes an excess head relatively to

the downstream delivery canal, a turbine scheme can be
envisaged in order to replace, totally or partially, the dissipation
structure;
along an irrigation canal, significant topographic uneven can
be used by a diversion scheme out of irrigation period (e.g.
whenever exists uneven between canal branches, with different
water levels).
- 16 -
Fig. 3.2 Diversion canal for a powerplant.
Fig. 3.3 Scheme of a power generation in an irrigation system.

Diversion works for the powerhouse.
- 17 -
Power generation and flood prevention - Dams can be used to prevent

floods through the creation of reservoirs that should be emptied
ahead of a rainy season, although allocating a certain volume in the
reservoir for power generation. An integrated management of
reservoirs operation tends, nowadays, to reduce the energy loss to a
minimum and to benefit water users. Nevertheless in small
hydropower plants this is not enough significant.
Fig. 3.4 Visualisation of the dry zone to absorb floods.
Power generation and fish protection - Damming a river impedes fish

migration and the passage between downstream to upstream sides
and vice-versa, when it is applied, being necessary to built fishpasses and to assure the navigation, in case of it is applied.
However, the reservoir would create a special spot for recreation
and breeding conditions for aqua-culture.
Fig. 3.5 Fish-ladders of pool-type with bottom orifice and top weir.
- 18 -
3.2- The benefit of small hydropower plant

Automatic isolated or autonomous small hydro plants provide an alternative
solution to electric grid extension that serves widely scattered communities, as a
high efficiency power supplement in urban areas, small industries and for
domestic purposes. Namely, there is a great disparity between urban and isolated
rural zones, with a consequent imbalance in the accessibility of energy.
Local conditions limit the adequate alternatives where is feasible to invest. The
global profit of watercourse management based on integrated programs related
to the water uses can indicate the best sites to install these types of plants. It is
convenient to know about the existing or future water uses and the needs
whether for water supply, either irrigation, or energy and for tourist purposes in
order to avoid incompatibilities between different needs.
Fig. 3.6 The effect on the landscape of low dams.
- 19 -
- 20 -
Hydrology
HYDROLOGY
M. Manuela Portela
4.1- Introduction. Scope of the studies

The main objectives of the hydrologic study of a small hydropower scheme are
the characterisation of:
the runoff at the water intake of the scheme in order to allow the
determination of the design discharge, and, thus, the design of the water
intake, of the diversion circuit and of the powerhouse, as well as the
evaluation of the energy production;
the floods or, more precisely, the peak flows, to consider in the design of the
weir, of some of the diversion works and of the powerhouse (for instance, if
the turbines are of the Pelton type they should be located above the water
surface elevation in flood conditions, at the powerhouse outlet.
The flood characterisation involves a domain, which is not specific of small

hydropower schemes. For this reason, this subject will not be treated herein.
However the dam works built along the river need to be prepared for typical
floods according to an accepted risk or design flood criteria based on a return
period.
- 21 -
Most of the time it is not easy to carry out the fulfilment of the previous
objectives, as there is not the required basic hydrologic data in the watershed of
the hydropower scheme.
In this chapter some simple methodologies that can be applied to the hydrologic
studies of small hydropower schemes when the available hydrologic data is
scarce are briefly and systematically presented. The proposed methodologies
partially overcome the non-existence of basic hydrologic information and with
some minor approximations allow the hydrologic and energetic characterisation
of the schemes under consideration.
More detailed analysis in the field either of the hydrologic studies or of theirs
application to small hydropower schemes can be found in books of hydrology
(LINSLEY and al., 1985, CHOW and al., 1988, MAIDMENT, 1992) and in
books specially dealing with the design of small hydropower schemes (ESHA,
1994, and JIANDONG and al., 1997)
4.2- Basic data required for the hydrologic study

The basic hydrologic data required for the evaluation of the energy production in
a small hydropower scheme is the mean daily flow series at the scheme water
intake in a period that has to be long enough in order to represent, in average, the
natural flow regime. By this way, it is reasonable to assume that the errors of the
estimates that result from the variability of the natural flows are minimised.
However, the small hydropower schemes, being frequently located in the upper
zones of the streams, have small drainage areas for which recorded stream flow
series are seldom available. For this reason, the inflows at the hydropower
scheme are usually evaluated by indirect procedures. Some of these procedures
utilise the transposition to the scheme water intake of the flow records in other
watershed, namely in the watershed of a stream-gauging station1.
For this last scenario, the required basic hydrologic data is next systematised. In
the next items, the methodology that allows the evaluation of the energy
production based on the previous data is presented. If the natural stream-flow
A gauging station where records of the discharges of a stream are obtained, LANGBEIN and ISERI,
1960, p. 19.
- 22 -
Hydrology
series at the section of scheme water intake is known, some of the procedures
presented in the next items become less important or, even, not applicable.
It should be pointed out that the decision of installing a stream-gauging station in
an ungauged watershed is mainly a long term decision that will only produce
usable results after several years of measures and, thereby, it is not a suitable
decision for design of any small hydropower scheme. In fact, only after a long
period of measurements it will be possible to achieve a series of flow records
having statistics characteristics that can be considered to represent the variability
of the natural flow regime. Once the small hydropower scheme is built and
started to operate it is advisable to measure the stream-flows in order to confirm
the design assumptions.
As mentioned, when the inflows to the hydropower scheme are not known, the
evaluation of the energy production can be based on the transposition to the
scheme water intake section of a mean daily flows series measured in another
site, namely in a stream-gauging station having a natural flow regime similar to
the one expected in the hydropower scheme.
When selecting the previous station one should guarantee that the corresponding
watershed and the one of the hydropower scheme have similar areas and are
close enough or both located in regions with similar hydrological behaviour, in
terms of similar mean annual values of precipitation and of runoff. The
hydrological comparison between the drainage areas of the stream-gauging
station and of the hydropower scheme can be based on maps of annual average
values of precipitation and runoff, which are published in most countries by the
appropriate water authorities.
The climate and the geological constitution, as well as the vegetal cover and the
human occupation and activity in the watersheds of the hydropower scheme and
of the steam-gauging station should also be similar. The period with stream-flow
records at the station has to be long enough, for saying, at least, twenty years.
In this conditions, it is necessary the collect the following data:
For the stream-gauging station:
series of mean daily flows (in m3/s) in a period that will be denoted by
nGS years;
- 23 -
in each year with known mean daily flows, the annual precipitation in
rain gauges located inside and near the watershed of the station.
For the hydropower scheme:
series of annual precipitation in rain gauges located inside or near the

watershed of the hydropower scheme. These records have to be
collected for a period as long as possible and during which all the rain
gauges have simultaneous data. This period will be denoted by nHS.
The available precipitation data usually allows a detailed description of the

spatial variability of the precipitation as the number of sample points of the
corresponding gauging network is non-comparatively higher than the number of
stream-flow sample points. The annual values of precipitation are generally
related to hydrological years or water years2 and are expressed as the uniform
vertical depths of water that would accumulate on the level surface of the
corresponding watershed if the precipitation remained where it fell (LINSLEY
and al., 1985, p. 55) annual precipitation depths (in mm).
The quality of the collected data should be tested, at least in terms of the general
consistency of the records. The most common test to evaluate the consistency of
the annual precipitation data at each rain gauge is the double-mass analysis
which compares the accumulated annual precipitation in a station with the
concurrent accumulated values of mean precipitation for a group of surroundings
stations (LINSLEY and al., 1985, p. 70). If the records are consistent, the points
thus achieved are displayed along a straight line (Figure 4.1).
Mean of the accumulated annual
precipitations (10 3 mm)
40
30
20
10
0
0
10
20
30
40
50
60
Accumulated annual precipitation (10
70
mm)
Fig. 4.1 Double-mass test.

2
Generally from October 1 to September 30 (LANGBEIN and ISERI, 1960, p. 21 and LINSLEY and
al., 1985, p. 116)
- 24 -
Hydrology
The application of the double-mass test to stream-flow records is, most of the
time, impossible, as there are not stream-gauging stations close enough in order
to allow the comparison of their records. So, the analysis of the consistency of
the annual flow records is usually based on a simple-mass curve analysis: graph
of the cumulative annual flows, generally as ordinates, plotted against time
(namely years), as abscissa. Once more, if the representation thus achieved
follows a straight line, the records are considered to be consistent (Figure 4.2).
6
Accumulated annual flow (10 m )

150
100
50
0
0
10
20
30
40
50
Hydrologic year
Fig. 4.2 Simple-mass test.
There are several other procedures, most of them of statistical nature, to

evaluate, in terms apart from the general consistency, the quality of the recorded
series. These procedures are, however, far beyond the scope of this book.
4.3- Mean annual and mean daily flow series

4.3.1 Introduction
Based on the hydrological basic data systematised in item 4.2, either for the
stream-gauging station or for the hydropower scheme, the mean daily flows
series at the scheme water intake can be performed according to the following
main steps:
Evaluation of the mean annual flow depth in the watershed of the

hydropower scheme, based on the mean annual precipitation depth over the
same watershed and in a correlation between annual flow and annual
- 25 -
precipitation. This correlation has to be previously established for the

watershed of the selected stream-gauging station.
By using the previous estimate of the mean annual flow depth, transposition
of the mean daily flow series measured in the stream-gauging station to the
water intake of the hydropower scheme.
4.3.2 Annual precipitation and annual flow. Mean

daily flow series
Once the annual records of the rain gauges located inside or near the watersheds
of the stream-gauging station and of the small hydropower scheme are collected,
the annual areal average depths of precipitation over the corresponding areas can
be evaluated. For this purpose the Thiessen method can be applied (LINSLEY
and al., 1985, p. 71, SMITH, 1992, p. 3.20).
In this method the stations are plotted on a map, and connecting lines are drawn.
Perpendicular bisectors of these connecting lines form polygons around each
station. The sides of each polygon are the boundaries of the effective area
assumed for the station. Over this area the precipitation is considered to be
constant and equal to the one measured in the station. The area of each polygon
within the basin boundary is determined by planimetry and is expressed as a
percentage of the total area Figure 4.3.
Fig. 4.3 Areal averaging of precipitation by the Thiessen method

(from CHOW et al., 1988, p. 79).
- 26 -
Hydrology
Weighted average rainfall for the total area is computed by multiplying the
precipitation at each station by its assigned percentage of area and totalling.
Besides the Thiessen method, other methods, as the ishoyetal one, can be applied
to evaluated the precipitation over each watershed. The presentation of these
methods and the discussion of their comparative advantages can be found in
SHAW, 1984, 209-216, LINSLEY and al., 1985, p. 71, and CHOW et al., 1988,
p. 78-80.
If the digital elevation models of the watersheds are available, the previous work
can be significantly simplified taking advantages of a GIS capabilities, where
others area averaging methods, besides the Thiessen method, are implemented
and easily manipulated.
Once the series of annual precipitation depth over each watershed is known, the
corresponding mean annual value is obtained: PGS (mm), for the stream-gauging
station, and PHP (mm), for the hydropower scheme
PGS =
PGSi
n GS
and
PHP =
PHP i
n HP
(4.1)
where Pi (mm) denotes the precipitation depth in year i over the watershed
identified by the sub index.
According to the methodology presented herein, the evaluation of the mean
annual flow depth in the watershed of the hydropower scheme requires the
establishment of a relation between annual precipitation and annual flow for the
watershed of the stream-gauging station.
Let PGSi (mm) denote the precipitation depth over the watershed of the
stream-gauging station in year i (i=1, , n) and HGSi (mm) the corresponding
annual flow depth obtained from the records of mean daily flows according to
365
i
H GS =
q GS k
k =1
A GS
86,4
(4.2)
- 27 -
where AGS (km2) is the watershed area and q GS ik (m3/s), the mean daily flow that
was measured during day k of the year i.

If the plotting of the annual precipitation, PGSi, versus annual runoff, HGSi, for a
common period of n years (i = 1, , n, n being as great as possible) displays a
good correlation between those two variables it will be acceptable to establish a
relation expressing, by means of statistics, the observed dependency (Figure
4.4).
Fig. 4.4 Relation between annual precipitation and runoff

(from LINSLEY et al., 1985, p. 255).
If a linear relation is assumed (which provides the simplest relation between

annual precipitation and annual runoff) the following expression can be obtained
by means of a regression analysis:
H=P-
(4.3)
where and are positive parameters that can be estimated from the sample of
annual values of precipitation and flow by the least square method (DRAPPER
- 28 -
Hydrology
and SMITH, 1981). The parameter should be less than one. The sub indexes of
H (annual flow depth) and P (annual precipitation depth) have been omitted as
(4.3) can be applied to a watershed different from the one to which the relation
was established (although both watersheds should have the similarities previous
pointed out).
Other types of relation between P and H can be assumed, for instances a
parabolic relation
H = P2 -
(4.4)
where and should also be positive.

The choice of the most suitable relation between annual values of precipitation
and flow should result from the graphical analysis of the behaviour of those two
variables (Figure 4.4) and from the hydrologic specific knowledge of the region
where the small hydropower scheme is going to be built.
Applying (4.3) or (4.4), or an equivalent expression, to the mean annual
precipitation depth over the watershed of the hydropower scheme, PHP (mm), the
corresponding mean annual flow depth, H HP (mm), is obtained.
The mean annual flow depth over the watershed of the stream-gauging station,
i
HGS (mm), can be evaluated from the series either of the mean daily flows, q GS k
(m3/s), or of the values of HGS i (mm) given by (4.2)
H GS
n GS 365
n GS
i
86,4 q GSik
H GS
i =1 k =1
i =1
=
=
A GS n GS
n GS
(4.5)
Let GS (hm3) and HP (hm3) denote the mean annual flow volumes at the
stream-gauging station and at the water intake of the hydropower scheme,
respectively
GS =
H GS A GS
1000
(4.6)
- 29 -
HP =
H HP A HP
1000
(4.7)
The transposition to the water intake of the hydropower scheme of the mean
daily flow series measured at the stream-gauging station can be finally
accomplished by applying the following relation:
i
q HP k = q GS k
HP
GS
(4.8)
where qHP ik (m3/s) is the estimate, for the water intake of the hydropower
scheme, of the mean daily flow during day k of year i.
Stream-gauging station
60
Mean daily flow (m3 s -1)
Hydropower scheme water intake
40
20
60
0
0
2000
4000
Time (day)
40
20
HP GS
)
0
0
2000
4000
Time (day)
Fig. 4.5 Obtainment of the chronological diagram of the mean daily flows at the
hydropower scheme water intake.
Applying (4.8) to the flow series measured at the stream-gauging station the
mean daily flow series for the water intake of the small hydropower scheme is
obtained. Once this series is known the correspondent chronological diagram can
be established Figure 4.5.
The application of relation (4.8) leads to an approximated idea of the mean daily
flow series at the water intake of the scheme. However, this idea is frequently
the best one, as alternative daily models require much more data and are often
highly complex, as the deterministic models that perform the transformation of
- 30 -
Hydrology
precipitation into runoff (rainfall-runoff process models). It should be noticed

that assuming relation (4.8) is equivalent to consider that the shapes of the
non-dimension chronological diagrams of the mean daily flow series3 at the
stream-gauging station and at the scheme water intake are equal.
4.3.3- Mean annual flow duration curve

Two techniques are available to determine the energy potential of a hydropower
site, namely, the flow duration curve method and the sequential stream-flow
routing.
The evaluation of the mean annual energy production in a small hydropower
scheme with a run-of-river exploitation was traditionally based on a mean annual
daily flow duration curve. This curve gives for each value of the mean daily
flow, Q, the average number of days per year, D, during which occurred mean
daily flows equal or greater than Q (MOSLEY and McKERCHAR, 1992, p.
8.27) Figure 4.6.
25
20
15
10
5
0
0
100
200
300
Time equalled or excedeed (day)
Fig.4.6- Mean annual daily flow duration curve.
Once the mean daily flows at the hydropower scheme are evaluated, the
corresponding mean annual daily flow duration curve can be obtained by
organising those flows by magnitude instead of chronological.
The flow duration curve can also be represented in a dimensionless form in what
concerns the mean daily flow series or both these flows and the time
Figure 4.7.
3
Chronological diagram of the mean daily flows divided by the corresponding average.
- 31 -
In the Figure 4.6, Qmod represents the modulus or the average mean daily flow
derived from the known mean daily flows by summation and averaging. For the
hydropower scheme QmodHP (m3 s-1) can be obtained by one of the following
expressions:
Q mod HP
n 365
i
q HP k
HP
= H HP A HP =
= i =1 k =1
365 n
0.365 24 3.6 0.365 24 3.6
(4.9)
where H HP (mm) is evaluated by applying (4.3) or (4.4) or an equivalent

relation, HP (hm3) is given by (4.7) and the meaning and units of the other
variables have already been presented.
Q/Qmod (-)
15
10
0
0
100
200
300
Time equalled or exceeded (day)
Q/Qmod (-)
15
10
0
0
25
50
75
100
Time equalled or exceeded (%)
Figure 4.7 Dimensionless forms of the mean annual flow duration curve.
Taking into account the procedures that led to the mean daily flow series at the
hydropower scheme water intake, q HP ik (m3/s), it is easy to conclude that the
- 32 -
Hydrology
dimensionless forms of the mean annual daily flow duration curves thus reached
for the scheme and the one relative to the stream-gauging are equal.
Once a dimensionless mean annual daily flow duration curve is put in a
dimension form by multiplying its ordinates by a modulus, Qmod, the curve thus
obtained will become exclusively representative of the watershed to which that
modulus is referred.
The mean annual daily flow duration curve at an ungauged site was traditionally
based on the adoption of a standardised regional non-dimensional duration
curve. The curve selected was next synthesised by multiplying its ordinates by
the hydropower scheme modulus, Qmod, in order to provide its particular flow
duration curve with numerical discharges.
However, the previous traditional methodology only leads to reasonable results
with respect to the energy evaluation when no storage capacity is available in the
hydropower scheme, that is to say, when the scheme has a pure run-of-river
exploitation. At the same time and even for a run-of-river scheme, this
methodology only gives the estimate of the mean annual energy production and
it does not allow any simulation study in order to analyse the variability of the
energy production due to the natural variability of the flows.
Regional standardised flow duration curves are generally available in all
European countries. The criteria that lead to the establishment of those curves
depend, however, on each specific country. When selecting a standardised curve
care should be taken in order to ensure that the watershed of the hydropower
scheme and the one inherent to the regional curve are similar (in terms of
hydrological behaviour, climate, area, occupation and geological constitution).
The procedures proposed in this chapter, namely those that result from relation
(4.8) are also regional transposition methodologies. The main advantages of
these procedures comparatively to the traditional ones based on the adoption of a
regional standardised mean annual flow duration curve can be systematised by:
In what concerns the factors that determine the similarity of the flow
regime at the stream-gauging station and at hydropower scheme, the
methodologies presented herein allow a more rigorous selection of the
daily data to be transposed.
- 33 -
4.4-
They provide not only a mean annual flow duration curve but also a
mean daily flow series allowing the evaluation of the energy
production either in a pure run-of-river scheme or in a scheme having
some storage capacity, in this last case, by means of simulations
algorithms.
In a general sense, whenever mean daily flows exist it is possible to
perform any simulation study. By this way, either the design of the
small hydropower scheme, or the comparison of different solutions for
the same become much easier.
Energy Evaluation
The evaluation of the mean annual production in a small hydropower scheme

with a pure run-of-river exploitation can be accomplished only by using the
mean annual flow duration curve at the scheme water intake. The simple
procedures that can be adopted in this situation are next briefly described.
For this purpose, let Qmax denote the design discharge of the hydropower
scheme. In this situation, the maximum mean annual volume, , that can be
utilised in the energy production is represented by the dashed area of in Figure
4.8 a).
If the efficiency of the powerhouse, , and the net head, h (m), are considered
constants, the mean annual energy production, E (GWh), in the scheme having
the flow duration curve of Figure 4.8 a) is given by
E=
h
3600
g
(4.10)
where g is the gravity acceleration (m s-3) and is expressed in hm3.

Most of the time, the average turbined volume will be less than as no turbine
can operate from zero flow to its rated discharge and, so, a minimum discharge,
Qmin, compatible with the energy production has to be considered.
A reserved discharge (for ecological purpose but also for others purposes as
irrigation or water supply) has also to be taken into account. In this situation the
- 34 -
Hydrology
mean annual turbined volume to be considered in the application of (4.10) will

be the one represented by the dashed area of Figure 4.8 b).
a)
3
-1
Mean daily flow (m s )
4
Qmax
0
0
100
200
300
b)
3
-1
Mean daily flow (m s )
4
Qmax
2
Qmin
Qres
0
0
100
200
300
Fig. 4.8 Maximum turbined volume in a pure run-of-river scheme: a) without

considering and b) considering a minimum turbined flow, Qmin, and
a reserved flow, Qres.
Usually neither nor h will be constant but, instead, will exhibit a

non-negligible variation. In this situation, as well as when there is a storage
capacity, the mean annual flow duration curve does not provide a correct way of
estimate the energy production and other procedures have to be implemented.
- 35 -
These procedures generally simulate the exploitation of the powerhouse by

means of a sequential stream-flow routing methodology. The simulation studies
can be developed according to different points of view, namely:
by adopting the mean daily flow series previously established for the
scheme, that is to say, assuming that the future flows will be equal to the
flows that occurred in the pass;
by adopting a synthetic mean daily flow series obtained from the series
previously establish for the scheme water intake by a model, generally a
stochastic model (for instances, a disagregation model).
As this last option involves models far beyond the scope of this work, only the
first option will be considered.
The simulation studies have to be performed by means of a computational
program. The main data generally required for this program is the following one:
mean daily flow series at the water intake of the scheme;

design discharge, Qmax (m3/s), and minimum discharge, Qmin (m3/s),
compatible with the operation of the powerhouse;
ecological discharge and any other discharge required for the
consumption between the sections of the weir and the outlet of the
powerhouse (reserved flow, Qres, in m3/s).
head losses in the diversion circuit as a function of the diverted flow;
rating curve at outlet section of the powerhouse;
efficiency curve of the equipment of the powerhouse.
If the weir of the scheme creates a lagoon providing a storage capacity and if the
operation of the schemes foresees the use of this capacity in order to improve the
conditions of the energy production, the following additional data is also
required:
the reservoir stage-capacity curve4;

the exploitation rules, that is to say, the rules that translate the way the
stored volume of water is going to be exploited, for instances, in order
A graph showing the relation between the surface elevation of the water in the reservoir, usually
plotted as ordinates, against the volume of water stored below that elevation, plotted as abscissa.
- 36 -
Hydrology
to concentrate the energy production during the periods where the sale
tariffs are higher (or when the energy demands are higher).
The simulation algorithm is almost exclusively based on the application of the
continuity equation.
Let t0 (usually t0 = 0 s) denote the initial instant and ti (s) the instant i t seconds
after t0, that is to say
ti=t0+i t
(4.11)
where t is the computational time step (s).

In these conditions the volume utilised in the energy production during the time
step i, i (hm3), will result from the net inflow also during the time step
inflow to the water intake, Ii (m3/s), minus the discharge over the spillway, Oi
(m3/s), minus the reserved flow, Qresi, (m3/s) and from the variation of the
volume of water stored in the reservoir, Si (hm3)
i = (Ii-Qi-Qresi) t + Si
(4.12)
The variation of volume Si is a consequence of the exploitation rules which,

generally, results into the two following different actions:
priority to the energy production;

priority to the storage.
The former of the previous rules is generally applicable to the peak hour period
and the latter, to the low hour period.
When priority to the energy production prevails, the turbined flow should be as
great as possible. For this purpose, the inflow to the scheme water intake is
increased (if possible, until the limit of the design flow, Qmax), by emptying
totally or partially the reservoir.
If priority to the storage prevails, as much water as possible should be stored. In
this situation the powerhouse will work only if discharges over the spillway are
foreseen.
- 37 -
The computation of the energy produced during the time step i, Ei (GWh), is
accomplished by an expression equivalent to (4.10) but where is replaced by
the volume turbined during the time step, i (hm3) and the variables h and
are replaced by hi and i in order to denote the specific conditions under which
the production occurred
E i =
i h i
3600
g i
(4.13)
The estimate of the energy produced in the hydropower scheme, E (GWh),

during the interval between instants t0 and tn, such as tn = t0 + n t, is obtained by
totalling the production during the successive time steps
n
E = E i
(4.14)
i =1
By omitting the term Si of equation (4.12), the previous methodology is also

applicable to the evaluation of the energy production in a hydropower scheme
with no storage capacity. In this situation, the values of the remaining variables
that appear in (4.12) will only reflect the particular conditions in what
concerns the net inflow, the net head and the efficiency of the equipment under
which the production occurred during each time step.
- 38 -
Hydraulic Design of Small Power Plants
HYDRAULIC DESIGN OF
SMALL POWER PLANTS
Helena Ramos
5.1- Introduction
Civil works are quite significant in Small Hydropower Plants (SHP). The
investment is generally limited and there is not capability to spend so much on
geological or survey exploration, hydraulic and structural design. These
undertakings have the main advantage in the use of local materials and attracting
local people for the construction. They can not be scaled down from large
projects. A SHP design should be the result of the work of a multi-disciplinary
engineering or multi-specialist team including hydrologic, hydraulic, structures,
electric, mechanical, geologic and environmental experts.
Several types of SHP layout schemes are, essentially, characterised by different
intakes and diversion structures depending on the type of the conveyance
system. The powerhouse depends on the type of turbines. The headwork
typically includes a low dam or barrage belonging to low dam risk category (e.g.
depending on the specific legislation of each country).
- 39 -
5.2- Layout schemes

5.2.1- Intakes
There are different types of intakes with diversion flow to the turbine through
the penstock: frontal, lateral, bottom drop and siphon type. For the two first
types it is necessary to design the entrance shape in order to avoid separated
zones of the flow and excessive head loss through wing walls. It is also
necessary for all types to verify the minimum submergence in order to avoid
vortex formation and, consequently, air entrance. The discharge control can be
improved by special structures (e.g. narrow sections, weirs and gates or valves)
depending upon if the flow is pressurised or of free surface type.
The elevation of the water level can be controlled by a weir or by an inflatable or
rubber dam or gate in order to create enough available head (Ho) to divert the
desired discharge (Qo).
5.2.1.2- Minimum submergence

The vortex formation by insufficient intake submergence can induce air
dragging or even solid material to the intake, reducing the turbine efficiency.
The vortex development will depend on the geometry of the intake, the
submergence and the approach flow velocity.
Fig. 5.1 Development vortex phenomenon (ASCE/EPRI, 1989).
The intake design criteria is based on the definition of minimum submergence

in order to enable the guarantee of non-vortex formation, with air dragging into
the hydraulic conveyance system. Several models and field studies were
- 40 -
developed in order to obtain design criteria like those based on Gordon

(ASCE/EPRI, 1989) formula (Figure 5.2).
Fig. 5.2 Intake scheme with submergence. Gordon symbols.
Gordon (in ASCE/EPRI, 1989) developed a criterion for intake design in order
to avoid the vortex formation (Figure 5.3). Two different types of flow
approximation (i.e. symmetric and asymmetric) were considered and the
following dimensionless equation was proposed:
S
V
=C
d
gd
(5.1)
where:
S is the submergence (m); d is the intake opening (m); V is the mean flow
velocity at the inlet (m/s); g is the gravity acceleration (9.8 m/s2) and C=1.7
(symmetric) or C=2.3 (asymmetric).
S/d
6
0,5
0,5
S/d=2,3 V/(gd)
S/d=1+2,3 V/(gd)
Penino and Hecker
0,5
S/d=1,7 V/(gd)
Gordon
3
2
S/d=0,5(V /(gd) -1)
Pressure pipe
1
0
0,5
1,5
2
0,5
V/(gd)
Fig. 5.3 Different criteria for minimum submergence.
- 41 -
Figure 5.3 shows that formula deduced by Pennino and Hecker (in ASCE/EPRI,
1989) is a conservative criterion in what concerns the vortex formation.
Different types of vortexes can occur in the water intakes. The following
classification of situations is proposed (adapted from ASCE/EPRI 1989),
considering main types (Figure 5.4): Type 1 developed vortex with a deep
nucleus and with drape air; Type 2 - superficial depression, without drag bulb air
but with well defined nucleus; Type 3 depression quite negligible with
unstable nucleus; Type 4 rotational movement without free-surface depression
but with superficial circulation.
Type 1
Type 2
Type 3
Type 4
Fig. 5.4 Scheme of different types of vortex.
Some of the procedures to consider, in order to avoid vortex formation beyond

the minimum submergence, consist in creating a good inflow approximation,
with a canal or by a convergent. If there is any singularity that provokes flow
circulation the minimum submergence criteria could be not enough. These types
of vortexes can occur close to partial opening gates, sluice valves or bottom
outlet.
- 42 -
Fig. 5.5 A vortex in a model facility with a visible air dragging.
A pproac h intak e c onditions

1B ad c onditions
Vortex type
2Good c onditions
3V ery good c onditions
4-
0,2
0,4
0,6
E uler num ber
0,8
Fig. 5.6 Relation between Euler number and the type of vortex
(adapted from NEIDERT et al., 1991).
Based on experimental tests and for the intake without type 1 vortex the
following formulation for the minimum submergence should be applied:
S 1 V 2 gd
=
1
2
d 2 E
(5.2)
- 43 -
with E, the Euler Number, that is defined as

E=
p
V 2
(5.3)
where
p - differential pressure between two sections (upstream and downstream of
the vortex);
- specific mass of the water;
V - mean flow velocity at the inlet.
The Euler number (E) is a dimensionless parameter that represents physically the
pressure drop by the velocity arising, which can influence the appearance of
vortex (Figure 5.6). For instance, if an intake has good approach conditions it
can appear type 1 vortex (with air drag) for E > 0.85 and it can be avoided, for
new operational conditions, characterised by E < 0.60. Very good approach
conditions are characterised by the non-existence of separated flow zone or any
type of singularity near the intake. When there are bad conditions (e.g. with
swirl conditions or presence of singularities near the intakes) the conservative
formulation of Penino and Hecker (in ASCE/EPRI, 1989) can be insufficient
requiring lab tests for a complementary analysis.
5.2.1.3- Bottom drop intake

The bottom drop intake is largely used in SHP (Figure 5.7). Due to its simplicity
it can be easily adapted to a low dam. The diverted discharge by this type of
intake structure is a function of the rack characteristics, namely the opening
degree or free area under non-submerged operational conditions. The design
criteria must have into consideration the removing capacity of solid material by
dragging.
The rack is normally located at the top of the low dam, allowing the absorption
of a discharge less or equal to the design one.
Should the turbulence be not significant over the weir where is placed the
bottom intake, the total head of the flow can be considered approximately
constant along the crest.
- 44 -
Ho
Upstream
Metalic rack
Downstream
Spillway
Overflow
spillway
qd
Collect
channel
Collect
channel
Low dam
Fig. 5.7 Bottom drop intake schemes.
The minimum length of the rack for diversion of the turbine design discharge is
determined by the integration of the following equations (NOSEDA, 1956):
dq d
= C d C r 2gh
dx
(5.4)
dh 2C d C r h (H o h )
=
dx
2H o 3h
(5.5)
with
Ho = h +
q2
2g h 2
being
Ho = specific energy over the rack;
h = head flow over the rack;
q = unit discharge over the rack (Q/L with L the total transversal length of the
rack);
qd = unit discharge absorbed by the rack;
Cd = discharge coefficient of the rack;
Cr = ratio area between free and total rack area;
x = flow and rack direction.
- 45 -
The combined integration of equation (5.4) and (5.5) results in the following
solutions:
L=
Ho
[(y 2 ) (y1 )]
Cd Cr
L=
and
Ho
[(u 2 ) (u 1 )]
Cd C r
(5.6)
in which
1
arc cos y(1 y ) and
2
1
1
2
2 cos + 1
(u ) = arc cos
(2 cos + 1)(1 cos )
2
2
3
with
h
q
(for head variation); u =
(for discharge variation).
y=
Ho
q mx
(y ) =
Equations (5.6) have the following positive roots:
y=
2
1
cos arc cos 1 2u 2
3
3
) + 13 - for subcritical flow
2
1
1
cos arc cos 1 2u 2 + 240 o + - for supercritical flow
3
3
3
1
= arc cos 1 u 2 - for subcritical flow
3
1
= arc cos 1 u 2 + 240 o - for supercritical flow
3
y=
These equations, based on NOSEDA (1956) study, are only valid for a free or
not submerged flow through the rack. This method allows the calculation of
head and discharge at upstream and downstream of the rack and the rack
absorbed discharge (Figure 5.7). Based on computational simulations, Figure 5.8
shows the influence of the discharge design on the rack length, for different
values of rack width.
- 46 -
length of the rack (m)
3,0
2,5
2,0
1,5
1,0
w idth = 2 m
w idth = 4 m
0,5
w idth = 6 m
0,0
0,5
1,0
1,5
2,0
design discharge (m3/s)
Fig. 5.8 Estimation of the rack length as a function of the design discharge.
water profile along the rack
8
water level (cm)
rack
lab tests
simulation
2
0
0
10
20
30
40
50
60
-2
-4
rack length (cm)
discharge variation
80
rack discharge
70
discharge (l/s)
60
50
lab tests (l/s)

simulation (l/s)
40
30
20
downstream discharge
10
0
0
10
20
30
rack length (cm )
40
50
60
Fig. 5.9 Comparisons between lab tests and computer simulations

(RAMOS and ALMEIDA, 1991).
- 47 -
In the calculations it was supposed a rack discharge coefficient of 0.63, a net

rack area coefficient of 0.60, a rack slope of 20% inclined to downstream and a
20% of obstruction by solid materials. This assumption is currently applied in
design phase for bottom drop intake type.
Another formulation based on CHARDONNET, MEYNARDY and OTH, 1967,
and posterior developed by a computer code (RAMOS, H. and ALMEIDA, A.
B., 1991) to calculate the discharge and head variation along the rack was
compared with experimental tests (Figure 5.9).
This type of intake (also known by Tyrolian, bottom-grade, streambed and drop
intake) has been applied in Alpine areas, in Europe. Worldwide practice shows it
is applicable to small rivers, in particular, in mountainous or hilly regions with
flows transporting a large amount of debris and big stones, where there are steep
gradients or falls provoking rapids or flash floods to facilitate the automatic rack
cleaning.
5.2.2- Protection rack

The rack is an important protection element against the inflow of solid material
that can provoke damage in turbines. Longitudinal bars lean against crossbars
with or not transversal reinforcement beams characterise each rack element. The
rack is located at inlet entrance and can be of fixed or movable type. The rack
can be defined by its spacebar, a, length in flow direction, b, thickness, c, and the
total cross-section, S.
c
a
b
Visualisation of two bars
Scheme of the rack
Figure 5.10 Characterisation of the rack.
In order to avoid rack obstructions by solid materials, namely when installed in

an approach canal, it will be appropriate to provide it with an auto-trash-rack.
The rack must be specified in order to avoid excessive head loss by grid
obstruction if the spacebars are too small, neither to allow driving solid material
- 48 -
into the turbine flow if the spacebars are very large. According to LENCASTRE,
1991, the spacebar dimension, a, is a function of the turbine type (Table 5.1).
Table 5.1 Spacebar for different turbine types
(LENCASTRE, 1991)
Turbine
type
a (m)
Kaplan
0.10 - 0.15
Fast
Francis
0.08 - 0.10
Slow
Francis
0.06 0.09
Pelton
0.03 0.05
Pumps as
Turbines
0.02
For a submerged rack with an automatic trash rack system, the maximum flow
approach velocity is about 0.80 to 1.00 m/s. This velocity is based on the
kinematic law:
V=
Q
S
(5.7)
where Q is the turbine discharge value (m3/s); S is the cross-section of the rack
(m2).
The flow through the rack can induce severe vibrations due to detached swirls
with a frequency, fs, which should be different of the bars frequency, fb, in order
to avoid resonance phenomena and the collapse of the rack. According to the
stability criteria, fb must obey to the following condition:
fb 1.5 fs
(5.8)
The swirl frequency, fs (Hz) is given by

fs =
St V
c
(5.9)
where
St is the Strouhal number (Table 5.2) that depends upon the cross-section of the
bar, whose value must be majored by the safety factor F, according to Table 5.3.
- 49 -
Table 5.2 Strouhal number for different types of bars (LENCASTRE, 1987).
Type
of bar
b=c
b=2.8 c
Diameter=c
b=2.8 c
b=2.8 c
b=5 c
St
0.130
0.155
0.200
0.255
0.265
0.275
Table 5.3 Safety coefficient for Strouhal number (LENCASTRE, 1987).
(a+c)/c
F
1.50
2.15
2.00
1.70
2.50
1.40
3.00
1.20
4.00
1.05
5.00
1.01
The calculation of the structural frequency of the bar of a rack with fixed
extremities, is based on the following equation:
fb = 3.6
gE b
c
3.46L + a
b
c
(5.10)
being L the distance between bar supports (m); g the gravity acceleration (9.8
m/s2); Eb is the elasticity modulus (N/m2); b the specific weight (N/m3) of the
bar material (for steel bars Eb = 2.1x1011 N/m2 and b = 78000 N/m3); is the
water specific weight (9800 N/m3).
In order to guarantee the rack structural stability, the dimension L must be
reduced until equation (5.8) be verified. Sometimes the solution leads to
consider a transversal bar in order to obtain a lower bar length.
- 50 -
5.2.3- Sedimentation or desilting basin

A sedimentation basin allows settling fine particles in suspension that may wear
parts of the turbine and to decrease its efficiency and lifetime. The settling is due
to gravity action through the reduction of the flow velocity by increasing the
sedimentation basin cross-section. Normally these basins are prismatic tanks,
where the flow has a main horizontal direction. The tank can be considered
divided in four zones:
1234-
inlet zone;
settling zone;
outlet zone;
sludge storage zone.
The inlet zone should assure a uniform flow velocity distribution by imposing
the flow through a perforated baffle. The settling zone allows the sedimentation
of the solid particles. Clean water is driven into the hydraulic conveyance system
through an outflow weir or orifice (Figure 5.11).
For design purposes, the solid particles are assumed as spherical shape, with a
uniform specific weight. The typical dimensions of the sedimentation (or a
desilting) basin are represented in Figure 5.12, for an inflow discharge Q at the
inlet of the basin with a depth D, a width W and a length L.
Outflow weir
Clean water
Target baffle
Q
Inflow with
solids
Settling zone
Inlet zone
Solids
Sludge zone
Outlet
zone
Fig. 5.11 Different zones in a sedimentation basin.
Through each basin cross-section, the free-surface flow is supposed to be

horizontal and with a uniform velocity. Particles tend to settle to the bottom with
a fall velocity Vs. Normally, the median diameter of the particles (D50)
considered admissible for small turbines is about 0.2 mm. This value will depend
- 51 -
on each turbine type (the turbine manufacturer should indicate the maximum
allowable particle dimension).
Q
V=Q/(WD)
Outlet
zone
D
Vs
Settling zone
W
L
Fig. 5.12 Schematic sedimentation basin.
The flow velocity is dependent of the basin main dimensions:

Q
V=
W.D
(5.11)
The sedimentation velocity can be obtained by the following equation:
Vs = 3.61
D 50 (d 1)
CR
(5.12)
where d is the solid particle density (ratio between specific weight of the particle
and the water, d =
p
, normally for silt type d 2.5); Re is the Reynolds
Vs D 50
, being the water viscosity =10-6 m2/s and CR the
shape particle coefficient.

number, R e =
Based on the design discharge (Q) and the sedimentation velocity (Vs) the
superficial area, As, necessary for settling particles of dimension greater than D50
is obtained by
As =
Q
Vs
(5.13)
- 52 -
Fig. 5.13 Spherical shape coefficient, CR, as a function of Reynolds number, Re

(adapted from QUINTELA, 1981).
The sedimentation criterion requires that the scour time Tscour of a particle must
be greater than the sedimentation time Ts, (Tscour Ts) in order to the particle be
settled.
Tscour =
Ts =
L
V
(5.14)
D
Vs
(5.15)
In practice the basin behaviour will be different from the one predicted by the
theoretical calculation. Main factors that influence the basin behaviour are the
flow turbulence and the non-uniform horizontal velocity distribution. One way to
obviate this problem consists in increasing the superficial area, As, of the basin
or the width, W, and removing all angles or unnecessary singularities.
The width and the flow velocity must be controlled in order to avoid the lift up of
settling particles (i.e. verification of the non-drag criterion). The equation (5.16)
gives the relation between the scour velocity and the geometry of the
sedimentation basin (ASCE/EPRI, 1989):
- 53 -
Vscour = 0.5 K R 1H6 C R (d 1) 50

304.8
12
(5.16)
with RH = hydraulic radius of the sedimentation basin, R H =
W.D
, and K
W + 2D
the Gauckler-Manning-Strickler coefficient, for concrete basin K = 75 m1/3/s.

The methodology consists in a set value for W, based on As, and the considered
L. Then D values will be estimated in order to verify that Vscour is greater than the
mean flow velocity, V.
Rack
Forebay
Approach canal
Intake
Fig. 5.14. Scheme of a forebay.
The depth, D, can vary along the length of the basin, as showed in Figure 5.14.
Nevertheless the obtained D will be an average depth.
5.2.4- Weir
The weir crest (in this case it will also be the spillway) is, normally, located in
the middle of the low dam. Sometimes it is convenient to install a gate in order to
increase the storage capacity. The full storage water level (FSWL) coincides with
the spill crest because it is the maximum capacity under normal operating
conditions. The weir will be designed for a discharge flow with a chosen return
period (e.g. 100 years) and should also be verified for greater discharges,
corresponding to the maximum flood water level (MFWL).
The typical spillway profile is defined in order to have atmospheric pressure
along its surface for design conditions. For discharge verification (Qmax>Q)
under-pressures along the spillway profile will be accepted, since the head is
greater than the design head (according to criteria based on Waterways
Experimental Station (WES) Corps of Engineers, US).
The general free-surface weir discharge law is given by
Q = CL 2g H 3 / 2
(5.17)
- 54 -
being:
Q = weir discharge (m3/s);
C = discharge coefficient (-);
L = weir width (m);
Ho = head over the spillway crest (m);
g = gravity acceleration (m/s2).
The discharge coefficient for a WES type weir is around C = 0.50, depending on
the head value (e.g. is superior to 0.50 for heads greater than design one and
conversely for lower heads).
Figure 5.15 shows the variation of discharge coefficient with the head if lateral
entrance wing walls exist in order to avoid flow lateral contraction.
The design head is obtained by equation (5.17) based on the design flood
discharge. The analytic definition of the spillway profile is obtained by the
following equation (WES profile):
y=
x 1.85
2H o0.85
(5.18)
The connection curve between the vertical upstream wall and the spillway profile
is obtained by an ellipse or by several arcs of curves according to Figure 5.16.
For verification discharge (Qmax>Q), the spillway profile will operate with a
pressure lower than atmospheric one and the discharge coefficient will increase.
For a smaller discharge (Qmin<Q) over-pressures (greater than atmospheric
pressure) will occur along the spillway profile with a smaller discharge
coefficient.
In order to avoid flow separation and cavitation risk for higher heads, the
following design criteria should be obeyed (LENCASTRE, 1987):
Hmax/Ho 1.4
According to the value accepted for the maximum or verification head (Hmax) the
wing lateral walls are fixed with a certain freeboard (e.g. 1 m at least) above the
MFWL.
- 55 -
Fig. 5.15 - Discharge variation of a WES spillway as a function of H/Ho

(ABECASIS, 1977).
Fig. 5.16 - Spillway profile for the design head (Ho)

(CORPS OF ENGINEERS, WES referred in QUINTELA, 1981).
In small dams it is very acceptable to choose a stepped spillway where the

downstream side is build with steps, where the flow can dissipate a certain
amount of energy (see Figure 3.6 - at bottom). When the available space for
locating the spillway is limited, labyrinth or semicircular solutions should be
considered, in order to obtain higher discharge.
A dissipation structure should be provided downstream the dam if the
geotechnical conditions are not strong enough.
- 56 -
5.3- Sluice bottom outlet

To empty the reservoir the dam should be provided with a sluice bottom
discharge. To obtain the emptying law of a dam sluice bottom outlet is necessary
to define the hydraulic conveyance system characteristics, namely singular loss
coefficients, the roughness and the cross section. Through the energy equation
between reservoir level, Hr, and the axe of the outlet works, Hd, (flow into the
atmosphere):
Hr H = Hd
(5.19)
where
Hr = is the upstream water level;
H = is the total head losses (as a function of Q - see 5.4 - Conveyance System);
Q2
, where A is the
Href = outlet water level plus kinetic outlet head =
2
2gA
outlet cross-section area.
It is important to know the velocity value at outlet, in order to prevent energy
dissipation structure (e.g. impact blocks), as well as to predict the safe height of
wing walls to direct the flow towards the river.
The bottom discharge will allow the total or partial reservoir emptying. To
estimate the emptying duration, it is common to consider a null inflow discharge
and the full storage water level (FSWL). Based on the discharge law for sluice
bottom outlet and on mass flow balance equation:
Qoutlet = A 2 g ( Hr H ref )
d
= Q inlet Q outlet
dt
(5.20)
where Qinlet is the inlet discharge into the reservoir (= 0 m3/s), Qoutlet is the outlet
discharge from the reservoir (m3/s), is the discharge coefficient and the
volume of the reservoir (m3).
The differential equation can be integrated by a finite difference technique in
order to obtain the total duration of the reservoir emptying.
- 57 -
5.4- Conveyance system

5.4.1- General layout
The conveyance system includes all elements designed to water transport from
the intake to the powerhouse. The conveyance layout can be composed either by
pressure galleries or pipes or by a mixed system composed by free-surface canals
and pressurised pipes. The layout of the conveyance system depends on several
factors according to the following remarks:
Small
dam
Analysis of different sites for low dam and powerhouse installation

and conveyance layout alternatives to seek for the economic net
head, the minimum length that involves smaller volumes of
excavation, and land volume to replace. All alternatives should have
good access and the concern about minimisation of expropriating,
must be performed. The selection of the more adequate final
proposes will be based on technical, economic and social
considerations.
The cross-section of the hydraulic system (e.g. characterised by a
diameter) should be the most economic in order to attend the annual
average energy production and the costs for conveyance system
installation.
Definition of all system accessories and components along the total
conveyance system and the more correct place to install them.
Canal
Small
dam
Forebay
Penstock
Powerhouse
Low pressure
pipe
Surge tank
Penstock
Powerhouse
Fig. 5.17 Schematic layouts for the hydraulic system.
Site selection for the headwork (e.g. a small dam and the intake structure) and
the powerhouse should be carried out on the basis of topographical and
geological conditions and on site visits. The main factors to be considered in the
alignment of canals or low pressure pipes are:
If a reach of a pipe can gain at least 2% of head, H, the experience

shows it is economically feasible.
- 58 -
A long canal could be unfeasible geomorphologically because of the

need to overcome side streams and gullies that imposes the
construction of crossing structures, such as inverted siphon.
Potential slides should be properly considered.
River bends associated to enclose geological conditions are
generally favourable to a straight tunnel or a combined system with
canals and low-pressure pipes.
Plastic pipes could be considered but there are limitations of
available diameters, or sunlight effects, or overpressures, or
temperature variations. These factors impose that the plastic pipes
should be buried or covered.
These constraints should be studied comparatively in order to obtain the best

solution.
For a total pressurised circuit with pipes following the land profile, blow-off
valves should be installed at the lowest points, allowing the pipe emptying, and
air-valves should be installed at higher points of the pipeline to avoid it
collapsing or buckling associated to air pressure problems. The highest point of
the conduit should be positioned below the lowest hydraulic gradient line under
extreme operating conditions (Figure 5.18).
Hmax maximum upsurge
Hydraulic grade lines
Ho steady state regime
Hmin minimum downsurge
Natural land profile
Air valves
Purge or
blow-off
valves
Penstock profile
Powerhouse
Fig. 5.18 Example of the hydraulic grade line and penstock profile.
Due to turbine discharge variations or turbine guide vane manoeuvres the

pressure will vary, along the penstock, with the time. The maximum and
- 59 -
minimum piezometric head envelopes should be indicated and compared with

the pipe resistance pressure classes and with the penstock profile, respectively.
The control of both the maximum and minimum transient pressure variations is
one of design objectives.
5.4.2- Head losses and net head

The net head, Ho, is one of the most important parameters for a feasibility
analysis and design of a hydroelectric system. The material, the cross-section and
the length of the penstock will influence the headloss. Net head calculation will
be based on the knowledge of the gross head (Hg - geometric water level or the
difference between water level at the intake (Zu) and the water level at the
tailrace (Zd)- for reaction turbines or the nozzle axe for action turbines (in
particular for Pelton)) and all head losses (friction and all singularities losses,
Hi, along the hydraulic conveyance system). The headloss will be a function of
the diverted discharge.
H o = Z u Z d H i
(5.21)
The design head is used to define the design power output of a turbine (typically
the maximum power output for the best efficiency head). For a hydropower
station with a high head scheme (e.g. with long penstock), the variation in the
water level can have less effect on power output than, perhaps, the variation of
the head losses with the discharge.
MFWL
FSWL
MEWL
Spillway
Dam
Fig. 5.19 Water levels at the reservoir.
When the reservoir water level changes during the powerplant operation, a
maximum water level may appear when the weir discharges the maximum flood
(MFWL - maximum flood water level). The crest of the weir gives the full
- 60 -
storage water level (FSWL) and the minimum exploitation water level (MEWL)
is determined by requirements of the intake operation (e.g. minimum
submergence criteria).
Friction head loss

The friction loss along a canal or a penstock is expressed by
H f = J L
(5.22)
where
J = hydraulic gradient;
L= length of the canal or the penstock (m).
Reservoir
Hydraulic grade line

Total singular losses
(curves, bends and
valves)
Gross head
Net head
Generator
Conveyance system
(low pressure pipe and penstock)
V2/(2g)
Turbine
Tailrace
Draft tube
Fig. 5.20 Hydraulic grade line of a hydro scheme equipped with a reaction turbine.
By assuming that the flow velocity is negligible in the reservoir, the gross and
net heads are schematically defined in Figure 5.20, for a reaction turbine with
downstream draft tube. In this case the internal draft headlosses are considered in
the turbine efficiency.
- Pressure flow
In closed pressurised pipes the flow is typically turbulent (Re > 2000), thus the
Colebrook-White formula is advisable for friction head loss calculation. This
formula requires an iterative method to solve the friction factor for each
discharge value. The Moody diagram allows a graphical calculation. For this
calculation it is necessary to know the mean velocity of the flow (U), the
- 61 -
diameter of the pipe (D), the absolute roughness (k) and the kinematic viscosity
of the water ( ).

1
2.51
= 2 log
+
f
3.7 R e f
(5.23)
where
k
UD
JD
=
Re =
2
D
U 2g
being f = Darcy-Weisbach factor; = relative roughness; Re = Reynolds number.
f=
In order to obtain the hydraulic gradient (J), it can be possible to establish the
following iterative process based on any initial arbitrary value Jn:
U 2 k
2.51
+
J n +1 =
log
8 g D 3.7 D D 2g D J n
(5.24)
The value of the absolute roughness (k) depends on the material of the conduit.
In Table 5.4 some k values are presented.
Table 5.4 Typical values of absolute roughness, k, for different
type of materials (BUREAU OF RECLAMATION, 1977)
material
cast iron
- new
- strongly rusty
soldier steel
- new
- rusty
concrete
- rough
- smooth
- 62 -
k
(mm)
0.25
1.50
0.10
0.40
0.60
0.18
The friction head loss calculation can also be based on empirical formulae. In
fact there are several empirical equations to estimate the hydraulic gradient (J),
but they are only valid for some flow conditions. Two of the most popular
empirical formulations are the Hazen-Williams (5.25) and Gauckler-ManningStrickler (5.26), which when expressed in SI units, are presented by the
following equations:
Q = 0.849 C S R 0.63 J 0.54
(5.25)
Q = K S R 2 3 J1 2
(5.26)
where K and C depend on the pipe wall material; R is the hydraulic radius
(R = D/4) and S is the pipe cross-section area.
Table 5.5 Typically coefficients of Gauckler-Manning-Strickler
and Hazen-Williams formulas
Materials
K (m1/3/s)
C (m0.37 /s)
steel
cast iron
concrete
PVC
90
80
75
110
130
120
120
140
- Free-surface flow
The Gauckler-Manning-Strickler formula (see equation 5.26) is currently applied
in open channel flows. This formulae due to its simplified application is often
chosen to obtain a quick estimation of the hydraulic gradient or unit headloss (J).
The following relation defines the hydraulic radius
R=
A
P
(5.27)
where A is the cross-section area of the canal and P is the wet perimeter of the
cross-section.
- 63 -
Singular or local losses

Local losses are expressed by a general equation of the following type:
U2
(5.28)
2g
where = coefficient of singular loss that depends mainly on the geometry of
the singularity and Reynolds number.
H L =
Typical values of are next presented (QUINTELA, 1981 and LENCASTRE,

1987):
Contractions
Type of
contraction
sharp edge
rounded edge
conical horn
gradual
<5
20
45
60
75
0.50
0.25
0.10
U2
2g
Piezometric
grade line
0.06
0.20
0.30
0.32
0.34
Hydraulic
grade line
Expansions
A
= 1 1
A2
(5.29)
Piezometric grade line
- 64 -
Bends
Separation
of the flow
30
40
60
80
90
0.20
0.30
0.55
0.99
1.10
Separation of
the flow
Gates and valves
h/D
0.20
0.50
0.70
31.4
3.3
0.8
5
20
40
60
0.24
1.54
10.8
118
gate
h
Racks
The flow through a rack induces a singular head loss that can be calculated as
follows:
4
c 3
= k sin
a
(5.30)
- 65 -
c
Types of
bars
a
2.42
1.79
0.76
0.92
where
= shape coefficient of a rack bar;
c = thickness or diameter of a bar (mm);
a = distance between bars (mm);
= angle of the bar with the horizontal;
k = coefficient rack obstruction.
5.4.3- Canals
The economic feasibility of small powerplants is based, among other factors, on
the available flow energy. The construction of long circuits in order to raise the
total net head is a solution depending on the topographic and geologic
conditions. Typically, a diversion scheme composed by a canal has a reservoir or
a forebay at downstream, where it is located the intake to the penstock. This
forebay will benefit the turbine flow changes.
The length of a canal is an important constraint that depends on the
technical/economic study. The height of the lateral walls of a canal must be
specified in order to avoid the overtopping by the flow including both the
unsteady and hydrostatic conditions.
A lateral weir is a way to control the canal water level (Figure 5.21). The
location of a weir depends on the terrain characteristics and strength against the
erosion provoked by the outflow. The weir will control the hydrostatic water
level and the height of the canal lateral walls. The implementation of automatic
control systems for water level regulation and discharge by gates are also
important solutions for the exploitation conditions in order to avoid an excessive
height of lateral walls.
- 66 -
Top of horizontal walls

Steady-state regime
Intake
A)
Hydrostatic
level
Forebay
Top of walls
Steady-state regime
Intake
B)
Weir crest
Hydrostatic
level
Fig. 5.21 Canal layouts: A) with top of horizontal walls; B) with a lateral weir.
5.4.3.1- Uniform and steady state

hydraulic regimes
The steady state uniform flow regime is a state of equilibrium between the
weight and the resistance force components along a prismatic canal with constant
roughness. Under uniform conditions, the free-surface profile, the hydraulic
grade line and the longitudinal bottom profile are straight and parallel lines.
In uniform regime the discharge, flow velocity, roughness and coefficient are
constants along the canal. For canals with small slope (), the hydraulic gradient,
J, is equal to the bottom slope:
U2/2g
J = sin tan = s
(5.31)
Fig. 5.22 Scheme of uniform hydraulic regime
By assuming that the flow is a rough turbulent flow, the Gauckler-ManningStrickler formula is currently applied due to its simplicity, with enough accuracy.
- 67 -
For non-uniform steady flows the free-surface profile or backwater curve is

obtained in engineering by the following formula:
dH
= sin J
dx
(5.32)
U2/2g
with
H = h cos +
Q2
2gA 2
(5.33)
Fig. 5.23 - Scheme of free surface steady state flow
For small slopes (cos 1), turbulent flow ( 1) and considering no lateral
discharge and a prismatic canal with a slope s = tan sin , the following
equation can be obtained by using a finite difference technique (standard step
method):
J + J n +1
s n
h n +1 h n =
x
2
Q Bsn Bsn +1
1
+
2g A 3n A 3n +1
(5.34)
where
A = cross-section area;
Bs =surface flow width;
g = gravity acceleration;
H = specific energy flow;
h = water depth;
J = hydraulic gradient;
s = bottom slope;
= angle of the canal bottom with the horizontal;
dx = space step discretisation along the canal.
The flow regime in a small powerplant canal is typically controlled by
downstream conditions (the forebay water level), being Q the turbine discharge.
- 68 -
5.4.3.2- Boundary conditions

The analysis of typical boundary conditions of diversion canals in small
hydropower plants is based on the following fundamental basic concepts:
- Upstream
Specific energy law from a reservoir to a subcritical flow canal;
Fixed water level downstream control gate with float;
Gate or valve that imposes a discharge variation law canal inflow
law.
- Downstream
Turbine discharge law a discharge demand law;
Lateral weir and water level regulation control of the forebay
water level, including an automatic control action through the
turbine guide vane or the control by the outflow through the weir.
5.4.3.3- Forebays
Sudden turbine discharge variations will provoke water level oscillations along
the diversion canal. A forebay can be considered as a regulation reservoir
(PINHEIRO, 1989), in order to reduce the water level variations and to improve
the canal response to turbine discharge variations and, can also operate as a
protection against silt or floating particles (see Figure 5.14).
When the plant demands a greater discharge, the water level quickly draws down
while the canal cannot supply enough flow. Otherwise, when the plant shutdown
a hydraulic bore will propagate upstream while the canal is still supplying the
forebay. This last event can induce secondary oscillatory waves and the canal
wall overflow.
A forebay (Figure 5.24) positioned at downstream end of a canal has its size
conditioned by the following factors (PINHEIRO, 1989):
1) To assure conditions to install the penstock intake with its equipment
(e.g. trash-rack, level detectors, sluices, gates, and weirs) always with
the minimum submergence criteria verified.
2) To limit the flow oscillations along the canal by turbine discharge
variations.
3) To assure the regulation function (e.g. to allow the transient turbine
demand satisfaction independent of the flow regime).
- 69 -
In a simple way, the modelling of the flow inside a forebay can be considered as
a reservoir with horizontal free surface in each time step, because the forebay is,
normally, much deeper than the canal.
Intake
Canal
Forebay
Forebay
Canal
Weir
Ice
sluice
Sand
sluice
Intake
Weir Sand sluice
Ice
sluice
Forebay
Intake
Canal
Weir
Sand
sluice
Fig. 5.24 General forebay layouts.
Applying the continuity equation at the forebay, the level oscillations can be
obtained by the following relationship:
dN f (Q in Q out )
=
dt
Af
(5.35)
where
Nf = free surface level in the forebay;
Qin = forebay inflow, from the canal;
Qout = forebay outflow, over the weir and for penstock;
Af = horizontal area of the forebay.
Typically Lf >2.5 Bf (being Lf the forebay length and Bf the forebay width) and
the velocity in the forebay is less than 0.5 m/s, in order to induce settling of the
harmful solid particles.
- 70 -
Lateral
weir
Intake
Maximum level
Nf
Normal level
Minimum level
Trash-rack
Trash-rack
Canal
Penstock
Sand sluice
Lf
Fig. 5.25 A typical profile of a forebay.
5.4.4- General remarks about mixed circuit

Depending on local conditions, a mixed circuit solution can be chosen. The flow
velocity and the shape of its cross section can be defined through practical
considerations. The canal is specified for operation with the maximum turbine
discharge, although it can also operate under much smaller discharges (during
dry seasons), which can induce serious silting problems along the canal or in the
forebay.
Small dam
Canal
Forebay
Small dam
Penstock
Powerhouse
Penstock
Powerhouse
Storage canal
Fig. 5.26 Two different types of mixed circuits.
The main disadvantages of a canal, comparing with a total pressurised circuit, are
the following factors:
The canal length will be greater than a straight pressurised tunnel

(unless a free surface flow tunnel be considered).
Potential construction troubles of the canal along the river
border.
- 71 -
Crossing structures to overcome side streams and gullies.

People and animal passages and the danger to fall into the canal.
Ice troubles in cold regions.
Efficient drainage systems for sediments.
Environmental impacts due to excavations.
Potential slope slides along the canal.
Maintenance costs.
Alternatively, in many cases the tunnelling is preferable. However, tunnels could

present, normally, difficulties during construction, especially caused by geologic
fault or shear zones.
5.4.5- Penstocks
The penstock must be designed to bear the maximum internal pressure due to
water hammer phenomenon during normal and abnormal operational conditions.
It is always laid on a stable site and towards the hill-slope. For small slopes, the
pipe must be buried to avoid temperature effects, or these effects must be
considered in the pipe design.
There are two suitable layout schemes, namely separate penstocks for each unit,
or a common penstock, with branching pipes only close to the powerhouse. The
principle of maximum economy and safe constructions must be taken into
account in the more suitable layout selection.
Economic diameter
The economic diameter is obtained considering the incremental energy benefits
from a lower energy loss associated to a larger diameter, as well as stability of
operational conditions (i.e. water hammer and wear) with the total investment
cost increase.
The selection of the more suitable diameter for each pipe branch can be based on
simplified economic study, with the main objective the minimisation of total
cost (e.g. dependent of unit pipe cost and of the loss energy by friction loss):
H
T
dC
AF
= K u1 p w o C p D o K u 2 C E o2 Q 3 19 / 3 = 0
dD
K
Do
where
AF = actualization factor;
- 72 -
(5.36)
CE = equivalent tariff of energy (cost unit);

Cp = pipe material cost per unit weight (cost unit);
Do = optimum diameter or economic diameter (m);
Ho = net head (m);
K = Gauckler-Manning-Strickler coefficient for the pipe material (m1/3/s);
Ku1, Ku2 = numeric factors to unit conversion;
Q = equivalent turbine discharge (m3/s);
To = average annual time operation (h).
= admissible pressure for pipe material (N/cm2);
p = specific weight of the pipe material (N/m3);
= specific weight of the water (N/m3).

Nevertheless, empirical formulations are recommended to economic diameter
estimation for small hydropower plants by several authors, on the basis of data
from existing penstocks (JIANDONG et al., 1997):
Do = C EC C MP Q 0.43 H o0.24
(5.37)
where
CEC = coefficient of energy cost (zones where the energy cost is low = 1.2,
medium = 1.3 and high or no alternative source = 1.4);
CMP = coefficient for the pipe materials (for steel = 1; or plastic = 0.9);
Ho = net head (m);
Q = design discharge (m3/s).
Another criteria based on the maximum velocity flow in the penstock can give a
first estimation for the diameter: 2-3 m/s low head plants; 3-4 m/s for medium
head and 4-5 m/s for high head plants. These are typical values based on real
cases.
For penstocks different types of materials can be used:
1) steel and cast iron pipe is widely used in hydropower design and is
preferable for very high head plants that allow a long service life. Expansion
joints should be provided at certain intervals, depending on local climatic
conditions;
2) prestressed concrete pipe presents a reduce cost but greater difficulty in
installation;
3) plastic and glassfibre-reinforced plastic pipe induces small friction loss, but
must be buried or wrapped to protect it from sunlight effect;
4) reinforced concrete pipe, should sewage pipe producers exist and it has a
long service life, little maintenance, low cost, but high difficulty in
installation and resistance problems.
- 73 -
Minimum thickness of steel pipe

The critical condition of pipe minimum thickness corresponds to the pipe crush
condition due to sub-atmospheric-pressure by vacuum occurrence because it is
the much warning condition. For this situation, when a long pipe is forced by an
external pressure p2 greater than internal pressure p1, creates a differential
critical pressure pcr that can be obtained by the following equation:
p cr = (p 2 p1 )cr
2E t
=

1 2 D
(5.38)
being E the steel elasticity modulus (E = 2x1011 N/m2), the steel Poisson
coefficient ( = 0.3), t the pipe thickness and D the pipe diameter. By
substitution of the former values on equation (5.38), yields the following result:
t
p cr = 4.53x10
D
11
(5.39)
According to European Convention for Constructional Steelwork (ECCS) in

order to take into account eventual geometric imperfections, it is proposed the
following relation:
p atm =
7.35
p cr
F
(5.40)
where patm is the atmospheric pressure (patm = 1.012x105 N/m2) and F a safety
factor (e.g. F = 2).
Based on equations (5.39) and (5.40) the minimum steel pipe thickness for a
given diameter (D) will be obtained:
t
8.4 x10 3
D
(5.41)
In order to consider corrosion effects it is adopted more 1 mm to the pipe

thickness.
Fig. 5.27 Resultant forces produced by weight, pressure and change in momentum.
- 74 -
Net forces due to pressure and momentum change are computed as follows:
Fpx = p A (cos 1 cos 2 )

Fmx = QU (cos 1 cos 2 )
Fpy = p A (sin 1 sin 2 )

Fmy = QU (sin 1 sin 2 )
(5.42)
where p = pressure in penstock; A = cross-section area of penstock; Q =

discharge; 1 and 2 = angles as shown in Figure 5.27.
At each pipe change direction, the penstock and its supporting structures must be
designed to resist the forces resulting from changes in direction (Figure 5.27).
5.5- Powerhouses
The function of powerhouses consists in housing and protecting turbo-generator
groups and the auxiliary equipment (e.g. safety and protection valves, electric
boards, control equipment, remote controller, switchgear panel and protection
equipment). The powerhouse layouts need to allow an easy installation of the
equipment as well as access for inspection and maintenance of the turbines and
all other equipment. Generally speaking, the dimensions of powerhouses are
mainly determined by the size of the generating unit(s) and equipment
(Figure 5.28).
For small power plants with long hydraulic conveyance circuit (e.g. of diversiontype) the head is generally greater than 40 m, even up to several hundred meters
and the discharge is, normally, smaller than 10 m3/s. As a result, open-flume
Francis, propeller or Kaplan and S-type turbines are not suitable. Experience
shows that horizontal Francis and impulse (e.g. Pelton, Turgo and crossflow
types) type of turbines are the appropriate. For horizontal shafts, the civil work
costs can be reduced 20% due to smaller height, easier inspection and
installation of protection devices, such as flywheels. A hand-operated travelling
bridge crane, with a maximum capacity to move the heaviest part of the
equipment (e.g. 5-15 tons) is used to assemble generators, turbines, protection
valves and other components. The crane type should be considered during the
design through the definition of the layout and dimensions of the powerhouse. In
small hydro an alternative layout should be considered based on the use of
mobile cranes and large opening on the building roof (Figure 5.28). For
economic and environmental reasons the powerhouse should be as compact as
possible, in order to minimize the landscape disturbance.
- 75 -
One of the most important construction constraints is the noise prevention

outside the powerhouse. This is a crucial aspect should the powerhouse be
installed near building or urban areas. To avoid the external noise all the
building openings need to be prepared with special isolation devices and noise
sources should be carefully identified (e.g. generator ventilation or flow through
the turbine wheel).
The generator can be ventilated by natural air. The control panel is housed in an
isolated and easy accessible part of the powerhouse.
1 gate isolation; 2 counting board; 3 6 kV equipment; 4 batery and loader; 5 auxiliary services board; 6
comand and control switchgearr; 7 automation switchgear; 8 protections switchgear; 9 - tension regulation
board; 10 support table; 11 bookcase for technical documents; 12 60 kV equipment; 13 main transformer
6/60 kV 5500 kVA; 14 auxiliary transformer 6/0.4 kV 50 kVA; 15 Pelton turbine with 4 nozzles; 16
synchronous generator of 5500 kVA 6 kV; 17 valve of isolation DN 800; 18 oil hydraulic central; 19 well
of drainaige and pumping; 20 workbench; 21 tools-case; water reservoir for WC
Fig. 5.28 A typical example of a small powerhouse with a Pelton turbine installed.
- 76 -
Different characteristics of the units will impose different powerhouse layouts

and sizes (see Figure 5.29 and 5.30).
The design of a powerhouse needs to take into account the hydraulic constraints
for a good turbine inflow, the weight of the equipment and hydrodynamic forces
that will be transmitted to the structure including the massive concrete (anchor
blocks).
Pelton
turbine
Spiral Francis
turbine
Open-flume Francis turbines
(adapted from MACINTYRE, 1983)
Fig. 5.29 Different types of powerhouse layouts.
Fig. 5.30 Powerhouse floor area required for Pelton and high head Francis turbine
installations (adapted from FRITZ, 1984).
- 77 -
In some powerhouses the control panel and the switchboard need to be properly
arranged at a level above the maximum design flood in the river, in order to be
protected against at least 500-year flood. For impulse turbines it must be
provided that the runner is non-submerged during operation.
5.6- Analysis of hydropower schemes

After the definition of the total characteristic of each component of a
hydropower scheme (see following chapters 6, 7 and 8), it is finally possible to
obtain the final design and the best solution. Figure 5.31 shows the importance
of the identification of different elements to input in a complete computational
model in order to obtain the dynamic response of the system.
Fig. 5.31 Typical inputs for hydropower scheme analysis.
- 78 -
____________________
Small Hydraulic Turbines
SMALL HYDRAULIC
TURBINES
Helena Ramos
6.1- Type of turbines

The choice of standardised turbines for small hydro schemes depends upon the
main system characteristics: net head, unit discharge and unit power. Hydraulic
turbines convert hydropower energy into rotating mechanical energy. The main
different types of turbines depend upon the way the water acts in the runner:
a free jet at atmospheric pressure - impulse turbines;

a pressurised flow - reaction turbines.
Impulse turbines are more efficient for high heads. The Pelton turbine is the
most known model of this type and is composed by a runner and one or more
nozzles. The runner has blades with the shape of a double spoon (Figure 6.1).
The jet coming from the nozzle hits the blades of the runner, transforming the
flow kinetic energy into rotational mechanical energy. Each nozzle has a
movable needle to control the discharge. The maximum number of nozzles is
two, for horizontal shaft, or six for vertical shaft. The nozzle has a deflector,
which is a device to control the flow whenever a load rejection occurs,
provoking a deviation of the jet enabling its slow closing, controlling the
overpressure in the penstock and avoiding the overspeed of the runner.
- 79 -
Fig. 6.1 - The trajectory of a jet from a nozzle in blades

direction of a Pelton (adapted from MACINTYRE, 1983).
The main advantages of these turbines are:
They can be easily adapted to power variation with almost constant

efficiency.
The penstock overpressure and the runner overspeed control are
easier.
The turbine enables an easier maintenance.
Due to the jet, manufacturers of these turbines impose a better solid
particles control inducing, consequently, a lower abrasion effect.
The Pelton turbine wheel diameter is usually 10 20 times the nozzle jet
diameter, depending on the spacing of the buckets. The net head in this type of
impulse turbine is measured between the total head in the penstock just
upstream the nozzle and the axis level of the water jet. Often runners are
installed on both sides of the generator (double- overhang installation). Impulse
turbines are provided with housings to prevent splashing, but the air within this
housing is substantially at atmospheric pressure.
Reaction turbines have, normally, a closed chamber (spiral case), where the
flow takes place in transforming part of pressure energy into rotational
mechanical energy of the runner. A movable guide vane (or wicket gate) guides
the flow around the runner, making, simultaneously, the regulation of the
turbine discharge.
- 80 -
____________________
Spiral case
Runner
Draft tube
Wicket gate
Fig. 6.2 A scheme of a reaction turbine.
A shifting ring to which each gate is attached moves the guide vane. At
downstream of the runner follows the draft tube that due to its shape (section
progressively rising) allows the partial recuperation of the kinetic energy of the
runner. The main advantages of this type of turbine are:
It needs lesser installation space (e.g. the runners are smaller than
Pelton runners).
It provides a greater net head and a better protection against
downstream high flood levels (can run submerged).
It can have greater runner speed.
It can attain higher efficiencies for higher power values.
The arrangement of the turbine shaft as vertical or horizontal and with fixed
(Francis) or adjustable blades (Kaplan) are important factors to take into
account in the classification. Francis turbines can be of single or double effect.
Thus, depending upon the type and the main conditioning factors, a particular
turbine can be classified as a function of:
- installation type;
- number of runners;
- position of the runner shaft;
- net head available.
The turbine types are typically characterised by a well-known parameter: the

specific speed Ns defined in 6.2 section.
In Table 6.1 it is shown a general classification of turbines with typical nominal
head, discharge power and Ns.
- 81 -

Table 6.1 Resume of application range of standard turbines
(based on several manufacturers data)
Hydraulic Turbines
Reaction
Impulse
Bulb
Kaplan and
propeller
axial flow
Francis with
high specific
speed diagonal
flow
Francis with
low specific
speed
radial flow
Pelton
Turgo
Cross-flow
H
(m)
Q
(m3/s)
P
(kW)
Ns
(r.p.m.)
2 - 10
2 - 20
3 - 40
3 - 50
100 - 2500
50 - 5000
(kW, m)
200 - 450
250 - 700
10 - 40
0.7 - 10
100 - 5000
100 - 250
40 - 200
1 - 20
500 - 15000
30 - 100
60 - 1000
30 - 200
2 - 50
0.2 - 5
200 - 15000
100 - 6000
2 - 15
<30
0.01 0.12
Main differences between turbines (see Figure 6.3)
Impulse
turbines
Pelton turbine can have single or multiple jets and is used

for medium or high heads. The flow incises on the blades
and falls into the tailrace canal at atmospheric pressure. It
can not work submerged.
Turgo turbine is also an impulse turbine type but with
different buckets, when compared with Pelton. The water
enters into the runner through one side of the runner disk,
through it and then emerges from the other side.
Cross-flow turbine (also designated by Ossberger, Banki or
Mitchell) is used for a wide range of heads. The water passes
through a guide-vane located at upstream of the runner and
has a double action on the blades of the runner.
- 82 -
____________________
Reaction
turbines
Reverse-pump operates as a reaction turbine with reverse

flow, from the outlet to inlet. Since it has no flow regulation
(guide vane) it can only operate under approximate constant
head and discharge.
Francis turbine is of radial or mixed flow with adjustable
wicket gate. It is used for medium heads. In these turbines
the flow is uniformly distributed through the spiral case, that
is connected to the penstock. The turbine may have vertical
or horizontal shaft and it can be arranged in open flumes too.
Kaplan and propeller turbines are axial flow reaction types,
generally used for low heads. The Kaplan has adjustable
runner blades with adjustable (double regulated) or not guide
vane (or single regulated). Propeller is chosen when both
flow and head remain practically constants. Both types can
be arranged in an open flume or with a spiral case of
concrete or cast iron, similarly to Francis turbines. When the
head is low and the flow is high a Bulb unit and S conduit
can be used.
Fig. 6.3- Typical turbines for small hydropower plants.
In engineering practice, a standard type of turbine can be selected in an early

design stage from charts prepared by the manufacturers like the one presented in
Figure 6.4.
- 83 -
H (m)
Q (m3/s)
Fig. 6.4 Example of hydro turbine range application

(adapted from ALTENER, 1997).
6.2- Turbine similarity laws and specific speed (Ns)

The application of similarity laws are required in order to predict the behaviour
of a full size prototype through the interpretation of model tests. The similitude
theory of turbomachines is used in different applications. This theory requires
the verification of three similarity conditions: geometric, kinematics and
dynamic conditions. For geometric similarity, the turbine dimension and the
flow passage must obey to one geometrical scale. The kinematic similarity
implies equivalent velocity triangles at inlet and outlet of the runner. For the
dynamic similarity there are analogous action forces (e.g. equivalent force
polygon). A complete similarity between runners of different dimensions is
always difficult and scale effects can occur. Reynolds similarity is not valid
- 84 -
____________________
(due to scale effects of internal loss and efficiency), because Reynolds number
has lower value in laboratory model than in prototype (MATAIX, 1975 and
JACOB, 1994). JACOB, 1994 and RAMOS, 1995 verified that the use of
Froude similarity to guarantee the relation between the inertia forces and the
gravity forces is the same for the model and the prototype. It can also guarantee
the pressure gradient similarity for a given average velocity.
Any similarity is related with homologous relationships in model and in
prototype, in particular, to allow the definition of the specific speed of turbines,
as an important parameter of each set of similar turbines that characterises its
dynamic behaviour. Based on similarity laws, a full description of the external
and internal (inertia) forces balance acting on a control volume defined between
inlet and outlet runner sections, through momentum equation, will provide the
discharge variation.
Under similarity operational conditions, the turbine speed, head and power, both
in model and prototype, follow the general equation:
nop
nom
P
= m
Pp
1/ 2
H op
H om
5/ 4
(6.1)
that yields in the following specific speed definition:
Ns = n o
(6.2)
H1o.25
where no (rev/min) is the nominal rotational speed and Ns ((rev/min) in m, kW)

is the speed of a similar turbine with unit head and unit output power during
similar operating conditions. Thus, under the same conditions, the Ns value is
considered as constant for similar turbines.
Impulse turbines have low specific speeds, Francis turbines have medium values
of Ns and propeller or Kaplan turbines have high values.
It is important to analyse the influence of multiple turbines in the specific speed
value. In small hydropower, double Francis wheel, reaction turbines with more
than one-stage and Pelton turbines with several nozzles can be selected. For the
case of wheels installed in parallel the discharge will increase (Q = n . Q and H
- 85 -
= cte, with n the number of wheels or nozzles) and the specific speed
( N s ,n = N s n ) will increase too. For different stages of wheels the head will
increase (Q = cte and H = n H) yielding in a decreasing of Ns ( N s ,n =
Ns
).
n 3/ 4
The relation between real velocity components, absolute (V) and relative (or
meridian) (W) of a water particle and the blade speed (C), at the inlet and outlet
of the runner will depend on Ns value of each turbine:
V=W+C
(6.3)
The ratio between absolute velocities and Torricelli velocity defines the specific
or unit speeds.
v=
V
C
W
; c=
; w=
2gH o
2gH o
2gH o
(6.4)
The increase of the specific speed (Ns) implies variations on unit speeds
(absolute, relative and tangential), on velocity triangles (Figures 6.5 and 6.6)
and on turbine discharge. Two geometrically similar turbines have same unit
speeds and analogous velocity triangles.
The Ns parameter can also be a complement parameter for turbine selection in
an early design stage (Figure 6.7).
- 86 -
____________________
Fig. 6.5 - Velocity diagrams through a Francis runner.
Ns (kW, m)
Pelton
28
80
130
Francis
210
250
300
Kaplan
300 - 700
Fig. 6.6 Variation of inlet triangle velocity with the specific speed (or runner shape)
(adapted from MATAIX, 1975).
- 87 -
A - Courtesy of VOITH HYDRO Power Generation Group, 1997.
Range application
Ns (m, kW)
1000
100
Pelton
Francis
1 nozzle
Kaplan
propeller
bulb
10
8
10
20
40
80
Ho (m )
100
200
500
1000
B - Based on ALTENER, 1997.
Fig. 6.7- Example of turbine selection charts based on Ns parameter.
- 88 -
____________________
6.3- Turbine efficiency

The turbine efficiency is defined as the ratio of power delivered to the shaft to
the power taken from the flow. Hence, in general, the efficiency of turbines is
defined by:
=
BH
P
=
QH o QH o
(6.5)
where BH (N.m) is the torque (and P (W) the power) delivered to the shaft by
the runner, (rad/s) is the angular velocity, Q (m3/s) is the flow rate, Ho is the
net head on the turbine. The efficiency of various types of turbine change with
the discharge, as shown in Figure 6.8. As can be seen the impulse turbine
maintains high efficiency over a wide range of discharges. The efficiency of
propeller turbines is very sensitive these values. The Kaplan turbine (with
movable blades) maintains high efficiency over a wide range of discharge
values.
Turbines can not economically operate from zero flow to rated discharge. The
efficiency decreases rapidly below a certain percentage of the rated discharge.
Many turbines can only operate upward 40% of rated discharge.
The total unit and powerhouse efficiency will be obtained by multiplication of
other efficiencies (e.g. generator efficiency).
1 0 0
efficiency (%)
8 0
6 0
4 0
P e lt o n
C r o s s - f lo w
Tu rg o
F r a n c is
K a p la n
2 0
P r o p e le r
0
0
2 0
4 0
6 0
8 0
1 0 0
Q /Q m a x ( % )
Fig. 6.8- Typical efficiency as a function of a percentage of the rated discharge for
several types of turbines (adapted from ATERNER, 1997).
- 89 -
6.4- Dimensions of turbines

Dimensions of turbines can be obtained by manufacturer information or
estimated through empirical formulas.
Impulse turbines
The main criteria for design of impulse turbines, in particular, for Pelton
turbines, consists in the calculation of the number of nozzles, turbine speed and
the dimensions of the runner and, finally an economic comparison between the
head benefit vs. powerhouse costs.
One of the most important parameter that characterises the turbo-generator
group, allowing the comparison of its behaviour with other machine, is the
specific speed (defined in 6.2). It is important to remember that if a turbine has
more than one nozzle (N) the actual specific speed (Ns,N) relates with the
specific speed of a turbine (Ns) with only one nozzle, through the following
equation:
N s,N = N s N
(6.6)
The generator is an element coupled to the turbine that transforms the

mechanical energy into electric energy. For synchronous generators, its speed
depends upon the frequency of the electric grid (f(Hz) in EU f=50 Hz) and the
pair of poles of the generator (p):
p n = 60 f
(6.7)
The net head is obtained by

H o = N res N noz H T
(6.8)
where
Nres = reservoir (upstream) water level;
Nnoz = nozzle axis level (Nnoz=Nriver+Hs);
Nriver = water level at outlet powerhouse (e.g. for 10 years of return period);
Hs = freeboard between the Nriver and the nozzle axis for turbine with vertical
shaft, and between river and the lowest runner point for turbines with
horizontal shaft (SIERVO and LUGARESI, 1978):
- 90 -
____________________
H s = 1.87 + 2.24
Q
Ns
(e.g. Hs 2m) ;
Q = turbine rated flow;

HT = total head losses in the hydraulic circuit.
and the turbine output capacity is given by
P = Q Ho
(6.9)
In order to avoid turbines with low efficiency or bad designed characteristics it is

convenient to adopt Ns according to manufacturer data information (or by
available literature) based, normally, on the net head and results of turbines
already tested.
Ho
Nnozzle
Fig. 6.9 - Definition of the net head for a Pelton with one nozzle.
The nozzle can have different open/close positions and it can be defined by a
displacement X (of the needle). When the nozzle is completely closed X is null.
X
Needle
2 2
db
Fig. 6.9 Scheme of a Pelton nozzle.
- 91 -
According to experimental tests (VIVIER, 1966 and MATAIX, 1975) the

discharge coefficient can be considered constant (except for very small open
degrees) and the discharge will be proportional to the free section of the nozzle.
Before the simulation of nozzle manoeuvres it is necessary to calculate the initial
conditions (steady state) and the nozzle diameter. Based on the maximum
discharge and the free cross section defined by:
sen 2 2
S = sen D N X
X
2
Q noz = CS 2gH o
(6.10)
DN
1.2
where X is represented in Figure 6.9.
For each nozzle it yields a diameter given by the following equation
X=
DN = K
Q noz
(6.11)
Ho
where
Qnoz = maximum discharge of each nozzle;
= the open degree of the nozzle ( = 1 completely open; = 0 completely
closed);
C = discharge coefficient ( 0.97- VIVIER, 1966 and MATAIX, 1975);
S = flow cross-section area of the nozzle;
K = coefficient that depends upon the angle of the conical needle end (e.g.
= 1 K = 0.5445).
For constructive proposes, according to SIERVO and LUGARESI, 1978, the
external diameter D2 and the diameter related to the centreline of the blades D1
are
D1 =
60 K u 2gH o
(6.12)
n
D 2 = D1 (1.028 + 0.0137 N s ,N )
(6.13)
and the blade dimensions are a function of nozzle diameter (DN)

- 92 -
____________________
DN =
D1 N s , N
(6.14)
250.74 1.796 N s , N
H1 = 3.20D 0N.96
and H 2 = 3.23D1N.02
(6.15)
where
K u = 0.5445 0.0039 N s,N is the periphery speed coefficient;
N s,N =
n
H 5o / 4
P
is the specific speed for each nozzle (N = number of nozzles);
N
H1
H2
D2
D1
Fig. 6.10 A typical Pelton runner

(adapted from SIERVO and LUGARESI, 1978).
Theoretically the choice falls on the greater specific speed turbine solution,
because it corresponds to smaller dimensions. However, on the one hand it must
be based on economic assessment, that includes civil works and equipment costs,
on the other hand the mechanical disposition of the group, maximum peripheral
speed of the rotor, and the minimum dimension for powerhouse, could help to
find the best choice. It is important to be aware that turbines with greater specific
speed lead a smaller runner diameter. However, for Pelton turbines the number
of nozzles can constrain powerhouse dimensions.
The dimensions of a Pelton turbine casing depend on the outlet diameter D2 of
the wheel, where L gives the horizontal size of the casing. For prismatic casings
this value has been assumed equal to the average diameter of the circle inscribed
and circumscribed on the casing.
L = 0.78 + 2.06D 2
(6.16)
- 93 -
The distance G between the wheel centreline and the top of the casing and the
other dimensions for casing and for spiral case are defined by
G = 0.196 + 0.376D 2
F = 1.09 + 0.71L
H = 0.62 + 0.513L
I = 1.28 + 0.37 L
B = 0.595 + 0.694L
C = 0.362 + 0.68L
D = 0.219 + 0.7 L
E = 0.43 + 0.70L
(6.17)
Fig 6.11 Casing dimensions (adapted from

SIERVO and LUGARESI, 1978).
Reaction turbines
Any reaction turbine is composed by the following elements:
Spiral case - with decreasing cross section to downward direction to
transform the pressure energy into kinetic one.
Wicket gate (or guide vane) - guides the inlet of the flow into the runner,
delivers it uniformly and controls the turbine discharge.
Runner - radial or axial with or without movable blades.
Draft tube - pipe with increasing cross section to downward direction.
For reaction turbines, the flow is totally pressurised and the net head is defined
by
H o = N res N river H
(6.18)
- 94 -
____________________
and the rotor speed is similarly calculated as for impulse turbines. The estimation
of specific speed is obtained from manufacturer data as a function of head
(see 6.2)
Defined the turbine speed, the next step consists in calculation the runner
position in order to avoid cavitation. The admissible suction head (hsmax) is
defined by the difference between the characteristic runner section and the
tailrace level. When the suction head is negative, the turbine operates in backpressure.
h s max =
p atm t v
H o
(6.19)
where
p atm
= local barometric head (m);
tv
= vapour pressure (m);
= dynamic depression coefficient or Thoma coefficient;

Ho = net head (m).
Fig. 6.12 - Definition of the net head for a reaction turbine.
The barometric head depends on the local altitude and the vapour pressure from
the local temperature (Table 6.1).
- 95 -

Table 6.1 - Barometric head and vapour pressure
Altitude
(m)
p atm
Temperature
(C)
(m)
tv
(m)
10.35
0.089
500
9.75
10
0.125
1000
9.18
15
0.174
1500
8.64
20
0.239
2000
8.12
25
0.324
(adapted from Almeida, 1995)
Thoma coefficient ( ) gives the susceptibility to the cavitation occurrence of a

reaction turbine. This parameter depends upon Ns. Several curves have been
proposed beeing one of them presented by Bureau of Reclamation, 1976 for
reaction turbines with vertical shaft:
=
N1s.64
38652
(6.20)
with Ns in (m, kW)
Fig. 6.13 Shematic representation of suction turbine head (hs).
Out of normal operating conditions induce the cavitation phenomena

occurrence. As a result of cavitation development, the performance of the
turbine will be affected, since efficiency, appearance of vibrations and noise or
blades erosion:
- 96 -
____________________
whenever exists a discontinuity of the flow;

the surface of a blade with different stages of
vortex development: initially by separating zones,
after by bulb formation and consequent rupture
that provokes vibrations and serious noise;
vortex band cavitation at the inlet of a draft tube
during partial load operating conditions, especially
for Francis, inducing periodical pressure Fig 6.14 Vortex in the draft
fluctuation, noise and hard vibrations in the
tube of an experimental
facility (at LNEC, RAMOS, 1995).
turbine.
In order to predict turbine dimensions and the necessary excavation for the
powerhouse in an early phase of design, SIERVO and LEVA, 1976 present a
practical formula based on several parameters:
D 3 = 84.5K u
Ho
n
K u = 0.31 + 2.5x10 3 N s
94.5
D1 = D 3 0.4 +
N s
D3
D2 =
0.96 + 0.00038N s
H1 = D 3 (0.094 + 0.00025 N s )
Fig. 6.15 Reaction turbine runner

42
H 2 = D 3 0.05 +
N s
(6.21)
In the same way, the steel spiral case and the draft tube dimensions are estimated
by the following relationships (SIERVO and LUGARESI, 1976):
- 97 -
- Francis turbine
- Kaplan turbine
19.56
A = D 3 1.2
N
s
54.8
B = D 3 1.1 +
N s
A = D M 0.40 N s0.2
C = DM
49.25
C = D 3 1.32 +
N s
48.8
D = D 3 1.5 +
N s
)
N )
s
D = D M 1.59 + 5.74 x10 4 N s
63.6
E = D 3 0.98 +
N s
131.4
F = D 3 1 +
N s
E = D M 1.21 + 2.71x10 4 N s
72.17
F = D M 1.45 +
N s
41.63
G = D M 1.29 +
N s
96.5
G = D 3 0.89 +
N s
81.75
H = D 3 0.79 +
N s
(
)
L = D (0.88 + 4.9 x10 N )
M = D (0.60 + 1.5x10 N )
I = D 3 0.1 + 6.5x10 4 N s
4
203.5
N = D 3 1.54 +
N
s
140.7
O = D 3 0.83 +
N s
(
(1.46 + 3.24x10
B = D M 1.26 + 3.79 x10 4 N s
P = D 3 1.37 5.6 x10 4 N s
31.86
H = D M 1.13 +
N s
31.80
I = D M 0.45
N s
L = D M 0.74 + 8.7 x10 4 N s
M = D M / 2.06 1.2 x10 N s
N = D M 2.0 2.14 x10 6 N s

O = D M 1.4 1.67 x10 5 N s
16.35
P = D M 1.26
N s
- 98 -
____________________
22.6
Q = D 3 0.58 +
N s
R = D 3 (1.6 0.0013N s )
18.40
Q = D M 0.66
N s
R = D M 1.25 7.98x10 5 N s
S = N s / ( 9.28 + 0.25 N s )
(
U = D (0.51 7 x10
T = D 3 1.50 + 1.9 x10 4 N s

4
53.7
V = D 3 1.10 +
N s
33.8
Z = D 3 2.63 +
N s
Ns
201.51
S = D M 4.26 +
N s
T = D M (1.20 + 5.12 x10 4 N s )
102.66
Z = D M 2.58 +
N s
(6.56)
D3,M
A)
B)
Fig 6.34 Spiral case dimensions (A) and draft tube dimensions (B)
- 99 -
Through these interpolation functions the dimension of Francis and Kaplan

turbines can be characterised, allowing a better estimate of economic costs for
powerhouses and the definition of the relative position of each component.
- 100 -
Hydraulic Transients and DynamicEffects
de Almeida
HYDRAULIC TRANSIENTS A. Betmio
Helena Ramos
AND DYNAMIC EFFECTS
7.1- Introduction
When the turbined flow changes during the hydropower operation disturbance
will occur along the hydraulic conveyance system. Hydraulic transients are the
regimes caused by these types of disturbances during a change from one steady
state to another. Physically, the hydraulic transients provoke surface gravity
waves along the diversion canals and pressure elastic waves along the
pressurised pipes.
Since the early design phases of the hydropower scheme, the hydraulic transients
should be considered in order to find the technical specifications corresponding
to a more economic and safer layout. The flow changes are inevitable: any turbogenerator unit must, at sometimes, to start-up, to undergo changes of load or to
be switched off. Unpredictable events, like human errors, equipment failures or
environmental hazards, can also cause severe unsteady regimes.
- 101 -
In canal systems a flow stoppage at the powerhouse will induce a transient water
depth wave that propagates upstream. The canal freeboard and all hydraulic
components along the canal need to be prepared for this situation.
In pressurised systems all the components need to support the maximum
transient overpressures and under pressures due to flow changes.
The final optimum solution need to conciliate the hydraulic transients with all
collateral dynamic effects, as well as the interaction with the hydro-electric
equipment:
hydrotransients will be the major factor in definition of operational safety

conditions, in order to avoid pipe rupture, water column separation or air
entrance into the system and the overtopping of canal or forebay walls;
hydrotransients and flow control manoeuvres will influence the overall
system response, in what concerns the turbine runner overspeed, the
automatic control system (e.g. regulators and governors) and the stability
conditions as well as the structural efforts induced by the hydrodynamic
forces.
The type of methods of analysis to be applied will depend on the design phase
and on the characteristics of each hydrosystem.
Accidents due to hydraulic transients can represent a very important risk in what
concerns both economic and life losses and of the powerplant operation,
reliability, as well as the production quality.
Since the sixties, computer methods for hydraulic analysis and simulation were
developed and it is now possible to have a large number of techniques to deal
with hydraulic transients and the global dynamic behaviour of the hydrosystems.
In some cases the diversion system and the hydraulic conveyance circuit is a
mixed one (free-surface diversion canals and pressure circuits or penstocks) but
in some cases it is more economic and environmentally more acceptable to
select a full pressurised hydraulic system.
For small hydro-systems the general methodology of transient analysis can be
the following one:
A- Preliminary and feasibility studies and early design phases.
Preliminary transient (waterhammer) analysis for basic situations and

manoeuvres.
Objective: to guarantee a feasible and economic solution without special

protection devices or to predict operational constraints or the type of protection
to be specified later.
- 102 -
B- Detailed design studies for the all for tenders.
Detailed transient analysis and studies, including the selected protection systems,
in order to obtain the hydraulic response to normal and abnormal turbine
operational conditions and select the main parameters of the equipment.
Objective: to specify the main component characteristics, namely in what

concerns the hydraulic conveyance circuit (canal, forebay and penstocks or
tunnels and conduits) and the flow control equipment (safety and control valves,
time of manoeuvres, unit inertia, among others), based on estimated equipment
characteristics.
C- Final studies for construction and operation.
Detailed transient analysis and computer simulations including the

characteristics of the selected equipment and final specifications of the civil
works.
Objective: to verify the safety level of the hydro-system and to specify operation
rules and to support the software development for special automation systems.
The complete hydraulic analysis, including the transient regimes and the
interactive dynamic effects, is a complex topic that justifies, in large
hydroelectric schemes, research activities related to special phenomena and the
development of advanced computer codes. For small hydropower plants, the
restricted design budget will certainly impose the application of well-known
criteria and already existing operational methods of analysis. However, the
computer codes based on the complete unsteady and transient hydrodynamic
equations are, now, relatively easy to use and the old approximate waterhammer
analysis methods are not justifiable and can be dangerous from both the
economic and the safety point of views.
The hydraulic analysis will strongly depend on the type of electric grid to which
the small hydroplant is connected: 1) a connection with an isolated (or island)
grid will impose more severe constraints, especially in what concerns the
dynamic effects and the stability of regulation of the turbine speed; 2) a
- 103 -
connection to a large national electric grid with a much greater power production
will make easier the powerplant operation in what concerns the dynamic effects.
7.2- Canal systems

In a diversion canal, the water levels will vary during the transient regime
according to the flow changes and the system characteristics.
Normal operations will include the complete powerplant shutdown: for this
condition the water level will increase and after a period of time the wave
fluctuations will dampen and the water level will align itself with the static
headwater level or other situation compatible with the control gates and safety
spillways.
As soon as the flow through the turbines is locked a pressure flow will travel
upstream along the penstock and when it arrives at the pipe intake, the water
level at the canal system begins to rise until the depth, the velocity and, hence,
the discharge reach new values. The sudden increase in depth gives rise to a
gravity wave (or bore) which propagates with a relative velocity Vp:
Vp =
gA 2 (A 2 h G 2 A1h G1 )
A1 (A 2 A1 )
(7.1)
where 1 and 2 correspond to upstream and downstream of the wave front,

respectively; A is the flow cross section area; hG - depth of the gravity centre of
the area A.
In free-surface flow, the complete dynamic model is based on the Saint-Venant
equations, which must be written under the conservative form, in order to be able
to simulate the bore propagation (RAMOS, 1995):
U F(U )
= D (U )
+
x
t
(7.2)
where U, F(U) and D(U) are the following vectors:

Q
; D(U ) = 0
U = ; F(U ) = Q 2
+ gAh
Q
gA (s J )
A
- 104 -
(7.3)
and x = distance along the canal axis (m); t = time (s); A = cross-section flow
area (m2); Q = discharge (m3/s); h = water depth (m); s = canal bottom slope (-);
J = slope of the energy grade line (-); and g = gravity acceleration (m/s2).
A true bore will be formed should the surface elevation of the surge above the
initial level be more than 20% of the average initial depth. When the wave
height is smaller, a train of short waves (or movable undulate jump) will be
created. The computer solution of the system of equation (7.3) based on
adequate numerical techniques (method of characteristics or finite difference
techniques) and on the boundary conditions will give the depth variation along
the canal (see ALMEIDA and KOELLE, 1992).
In a detailed transient analysis and canal design, the secondary oscillations,
known by Favre waves, need to be added to the surge waves, or bores obtained
through the integration of Saint-Venant equations (7.2). PREISSMAN AND
CUNGE, 1967 present a methodology to calculate the amplitude of these waves
(Figure 7.2), based on Froude number for the propagation of a bore:
Fr1 = Fr 0 +
1 h1
1 +
2 h 2
(7.4)
where
Fr 0 =
V1
g h1
c
h1
V1
V2
h2
(A)
h1
Favre waves
V1
V2
h2
(B)
Fig. 7.1 Propagation of bores A) based on Saint-Venant equations;
B) considering Favre waves.
- 105 -
The limits of different types of configurations were obtained: in case of 1 < Fr1 <
1.3 will appear secondary waves overlapping the main wave (Figure 7.1); in case
of Fr1 > 1.7 breaking waves or moving hydraulic jump will appear.
Favre waves constrain the height of the canal walls, especially along the
downstream 2/3 length of the canal. However, it is necessary to consider other
factors that can influence the water level along the canal: sedimentation of solid
particles; wind induced waves; lateral discharge from hillside; variation of canal
wall roughness and a variation on the turbine maximum discharge. Considering
these factors, US Bureau of Reclamation (AISENBREY et al., 1978) and
INVERSIN, 1978 propose a minimum limit of 0.15 m for the wall freeboard.
h*/h'
2.3
2.1
Near the
near
thewalls
walls
1.9
B=h2
B=3h2
1.7
1.5
in the
axe
Along
the canal
canal axis
1.3
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
h'/h1
h
h1
h*
1/1
h2
B
Fig. 7.2 Maximum amplitude of Favre waves in trapezoidal canal with 1/1 lateral wall
slope (adapted from PREISSMAN AND CUNGE, 1967).
For a preliminary analysis and canal freeboard evaluation, the Feifel formula can
be applied in order to calculate the surge depth induced by the complete
downstream flow stoppage:
- 106 -

2
V2
V2
V2 A
h = 1 + 1 + 2 1 1
2g
2g B2
2g
(7.5)
where:
h = surge depth above the initial canal water level;
V1 = initial upstream flow velocity (before the flow stoppage in the
powerhouse);
A1 = flow cross section for the initial discharge and flow velocity V1;
B2 = surface width induced by the bore (this value needs to be evaluated by
iteration).
In powerplants equipped with action or impulse type turbines (e.g. Pelton
turbines), the turbined water will flow downstream turbines by a tailrace canal.
Should this canal be long enough and the transient free surface flow due to
turbine discharge variations should also be considered. For this type of
powerplant unit the turbine wheel should be always placed above the highest
downstream level, including the selected design flood level at the downstream
river, where the turbined flow will be conducted. The canal freeboard or the
ceiling level, in case of a closed canal or tunnel, need to consider these aspects
as well as the need for an air flow circulation during normal turbine operation.
If a forebay is placed between the diversion canal and penstock intake, their
geometric characteristics will modify the canal surge.
The canal-forebay characteristic
parameter PCC is defined by:
PCC =
h (m)
Qo = 5 m3/s
A F H c
Qo
L
Qo
+
gh
H c
gA
Qo = 2.5 m3/s
Qo = 0.5 m3/s
being
A flow cross section area;
AF horizontal forebay area;
h flow depth;
L canal length;
Qo initial discharge;
Hc water level difference
between upstream and
downstream end sections of the
diversion canal.
Fig. 7.3 Example of upsurge (h) wave attenuation due to forebay effect after an
instantaneous flow stoppage at the powerplant (L = 1000 m and s = 0.0005 adapted
from PINHEIRO, 1989)
- 107 -
For PCC=0 (no forebay) the h values will correspond to the no modified surge
or bore height. The increase in the horizontal forebay surface area will attenuate
the bore height h that will propagate upstream the diversion canal.
In order to guarantee the maximum powerplant head and generating power, an
automatic water level is typically placed in the forebay. At upstream side of the
diversion canal the flow changes can be caused either by natural river discharge
variation or by a gate action. At downstream side, the main source of level
variation is the variation of the turbine discharge. In a simplified way, for high
head plants with impulse turbines it can be admitted that penstock discharge is
imposed by the regulator through the turbine gate (nozzle gate) position. In this
case, the water level regulating can be based on a P.I.D. (Proportional, Integral
and Derivative) regulator for a stable and efficient control (ALMEIDA and
KOELLE, 1992 and RAMOS, 1995).
7.3- Pressurised systems

7.3.1- Typical transient regimes
When the flow through a penstock diminishes too rapidly, the pressure along the
pipe upstream the flow control device will rise above the initial level and might
cause the penstock to burst. At downstream side of the control device the
pressure variation will follow the opposite way, when the flow diminishes the
pressure will tend to lower. This phenomenon is typically known as the
waterhammer and is one of the most dramatic aspects of the hydraulic transients
(RAMOS, 1995).
In any pressure transient analysis the following steps need to be considered:
1- The physical origin of the phenomenon and its mathematical
characterisation.
2- The selection of scenarios compatible with the hydropower characteristics
and the evaluation of the transient pressure variations, as accurate as
possible according to the design stage.
3- The selection, analysis and specification of special protection operational
procedures or/and devices or components to control the transient pressure
variations and other harmful dynamic effects on the plant operation.
In what concerns the operational conditions to be considered in small
hydropower schemes the following regimes can be selected:
Normal operating conditions, as expected or specified, that should not

induce any difficulty or problem (maximum safety factors).
- 108 -
Emergency operating conditions, probable but unexpected they may cause

some inconveniences but should not strongly damage the hydrosystem
(average safety factors).
Exceptional operating conditions, very unexpected and highly improbable,
they may cause severe damage to the hydrosystem (minimum safety
factors).
In Table 7.1 some examples of operating conditions are shown.
Table 7.1 - Examples of operational conditions for transient analysis and

simulation (adapted from PEJOVIC et al., 1987).
Normal
Emergency
Exceptional
Steady-state flow regimes for different turbine discharges

and load elements.
Steady-state flow regimes for turbining no-load
conditions.
Transient operation to synchronous regime and loading
conditions.
Load rejection during turbine operation followed by
closure of the flow control equipment.
Penstock filling and emptying situations.
A failure on the closure mechanism of a turbine wicket.
A failure on a protection device against waterhammer.
The closure of the safety guard valve of one unit after a
complete load rejection of the plant.
Flow control under turbine runaway regime.
Rapid start-up of the plant followed by a failure of one
unit.
Unfavourable succession of operations with overlapping
of transient regimes.
Complete failure of the turbine wicket closure
mechanisms with a very fast flow stoppage (in a time less
than 2L/c).
More than one protection or safety device fails.
Complete blockage of all turbine wickets and closure of
all penstock safety valves.
Downstream or draft column separation and closure due
to the reverse flow.
Resonance or oscillation phenomena.
Breakdown of the penstock.
Seismic induced hydrodynamic actions.
The hydraulic transient effects will strongly depend on the overall hydrosystem
characteristics, including the number of units and type of turbines. In all cases
- 109 -
the main objective is to prevent any serious damage to the penstocks or other
pressure conduits, as well as to any other component of the hydrosystem. In
what concerns the structural penstock protection both the maximum and
minimum transient pressures need to be controlled in order to avoid that: 1) the
allowable maximum pressure be exceeded, to prevent a breakdown or pipe burst;
and 2) the sub-atmospheric pressure in order to avoid cavitation phenomena and
water column separation effects, as well as a potential pipe wall buckling event.
Typically, the allowable maximum relative transient head variations H/Ho will
depend on the design head:
Ho (m)
Allowable maximum
H/Ho
500 - 200
200 - 50
<50
0.15 - 0.20
0.25 - 0.35
0.50
In all cases the downsurge or minimum transient pressure should be properly

considered in order the penstocks or conduits be protected against the formation
of under atmospheric pressures and vacuum formation (water column
separation).
7.3.2- Preliminary analysis

Any disturbance induced to the flow will propagate as a wave with a finite
relative velocity known as the elastic wave celerity whose value will depend on
the elastic properties of both the fluid and the pipe or tunnel wall according to
this general formula:
c=
K
[1 + (K E ) ]
(7.6)
where
E - Young's modulus of elasticity of the pipe wall;
K fluid bulk modulus of elasticity (2x103 MPa for water);
- parameter that depends upon the pipe structural constraints;
fluid specific mass (103 kg/m3 for water).
For thin wall pipes =D/t (with D the pipe diameter and t the wall pipe
thickness) and the celerity formula will be for water and S.I. units:
- 110 -
c=
1425
(7.6a)
KD
1+
Et
Should the pipe or tunnel thickness be very large and the elastic wave celerity
will be 1425 m/s. The pipe or tunnel elasticity will make to lower the celerity.
The pressure fluctuations related to these elastic waves will have, in a uniform
pipe of length L, the period 4L/c.
The basic waterhammer theory (ALMEIDA and KOELLE, 1992) teaches us that
the maximum transient pressure variation PJ, due to a flow disturbance, follows
the Joukowsky formula:
p J = c V
(7.7)
where V is the flow velocity variation.

The negative signal is for upstream side of the flow control device and the
positive signal is for downstream side. Thus, for a flow velocity reduction p>0
(the pressure will increase) at upstream side and will be p<0 at downstream
side. For very long pipes with high friction headloss formula (7.7) for upstream
pressure variation needs to be modified:
3
H
1 H o
*
o

p J = p J 1 +

(7.8)
H J 12 H J
where Ho is the initial pipe head loss and HJ the head variation corresponding
to pJ (ALMEIDA, 1985).
The Joukowsky formula will be only valid for instantaneous or fast manoeuvres,
when the duration of flow change TC is less or equal to TE, the time of one round
trip of the first pressure elastic wave or the time for the elastic wave to propagate
from the source of the disturbance to the penstock intake, where it will be
reflected back with the same velocity, and returns to the disturbance source.
This situation can be avoided if TC is carefully chosen for normal operation
conditions in order that TC>>TE (slow manoeuvre). This could be easily obtained
by closing or opening the turbine flow control devices (nozzle or guide vane)
any safety valve or gate very slowly. For slow manoeuvres of complete flow
stoppage or start up), in real and normal design situations the approximate
relative head variation HM/Ho can be given by
- 111 -
H M
T
= KT W
Ho
TC
(7.9)
where KT is a factor that depends on the turbine type and operation (LEIN, 1965)
and TW is the hydraulic inertia time constant defined by
TW =
L Vo
g Ho
(7.10)
where
L = pipe length (m);
Vo = initial or final flow velocity (m/s);
Ho = reference net head (m).
For Pelton turbines KT will vary between 3.7 (closure) to 3.3 (opening) and for
Francis turbines KT will vary for both manoeuvres from 1.2 to 2.0 (being the
majored KT factor between 1.8 and 2.2). For certain theoretically conditions
KT=2.0 (Michaud formula mentioned in RAMOS, 1995).
A similar equation can be applied for the calculation of minimum heads or
pressures near the turbine gate. For pipe design, the transient maximum and
minimum heads need to be obtained along the axis profile (see Figure 5.18). For
a computer this information is easily obtained in each computational section. In
most of the approximate methods, a criterion for the head envelops needs to be
followed: typically, for slow manoeuvres, a straight line is adopted for each
envelope, from the upstream or the downstream side of the turbine to the fixed
reservoir or tailrace levels, respectively.
Based on the two extreme enveloped lines, the maximum and minimum
pressures can be easily found by comparing those lines with the pipe or penstock
axis profile along terrain.
For both normal and abnormal turbine operations it is necessary to be very
careful with the minimum transient pressures downstream the turbine runner (in
the draft tube) in order to avoid the water column separation, specially when
there is a long tailrace pressurised tunnel or pipe. An approximate criterion was
proposed by LEIN, 1965:
p rm
V22
= Hd H m
Zr
2g
(7.11)
where
Hd downstream level (e.g. downstream river level);
- 112 -
Hm head variation giving the minimum transient head downstream the draft
tube;
V2 flow velocity at runner outlet;
Zr turbine runner outlet section level;
prm minimum transient pressure at runner outlet section;
Water column separation will be avoided if prm/ > -6 m w.c.. This means that a
reaction type turbine runner elevation and the powerhouse elevation should be
fixed according, among other factors, to the downstream pressure transient
analysis.
7.3.3- Governing equations

The basic differential equations of the unsteady and transient pressure flows can
be written as follows (RAMOS, 1995):
U F(U )
+
= D (U )
t
x
(7.12)
where U, F(U) and D(U) are the following vectors:

H
U = ;
Q
c2
Q
F(U ) = gA ;
gAH
D(U ) = JgA
Q
Q
Q 2
(7.13)
in which x = distance along the pipe axis (m); t = time (s); A = cross-section flow
area (m2); Q = discharge(m3/s; H = piezometric head (m); J = slope of the energy
grade line; g = gravity acceleration (m/s2); and c = elastic wave celerity (m/s).
Equations (7.12) can be solved by computer codes based on numerical
techniques (e.g. method of characteristics) and on the specified boundary
conditions. Nowadays, the computer modelling is based on the discretisation of
the hydrosystem and the computational model will be composed by several
components (tube and non-tube elements) and nodes (ALMEIDA and KOELLE,
1992).
With these kind of models it is possible to simulate different operation conditions
and study the hydraulic behaviour of complex systems, including pipe networks
and all major hydro-mechanical equipment. Among these ones, the turbines are
the most important source of flow disturbances and their modelling is crucial for
a reliable computational model as a design aid tool.
- 113 -
The turbine modelling, especially in what concerns the reaction turbines (e.g.
Francis and Kaplan turbines), the turbine couples the upstream and downstream
sides of the pressurised hydraulic circuit. In these cases, the turbine modelling
implies the knowledge of the complete characteristic curves based on the
manufactures tests. These curves are very difficult to obtain for small
hydroplants.
In what concerns the action turbines (e.g. Pelton turbines), the pressure transients
will need to be considered along the upstream penstock. The disturbance source
will be the turbine nozzles, a special type valve that can control the flow by
changing the position of an internal part known by the needle. In this case, as in
all pressure transients induced by a flow control device placed downstream a
long pipe, it is very important to know the characteristic equation of such control
device, or the controlled flow for each nozzle (or valve) opening and each head
variation across the upstream and downstream sides or nodes of this specific
component.
Small hydropower schemes when installed in mountainous regions are typically
associated to long penstocks, in order to increase the available head, and to
action or impulse turbines (e.g. Pelton turbines).
With computational models it is possible to include the real nozzle
characteristics and the headloss variation along the penstock during the transient
regimes. The friction headloss can be characterised by the following pipe
parameter (RAMOS, 1995):
Fp =
H
fL
=
2
V 2g D
(7.14)
where H is the flow total head loss; V 2 2g is the kinetic head; f is DarcyWeisbach friction factor; L the pipe or penstock length; and D the pipe or
penstock diameter.
The computer simulations indicate that for high Fp values the flow control can
be very difficult to obtain, because the time duration of the physical manoeuvre
of the nozzle or valve differs very much of the real flow change duration TC. In
these cases, for a nozzle linear time closure, the discharge variation will be only
effective near the full closed position: the duration of the nozzle mechanical
manoeuvre is very different from the effective time of flow discharge variation
and a theoretically slow manoeuvre (TC>TE) can be, in fact, a fast one (TC<TE) as
- 114 -
can be seen in Figure 7.4 (where Qo is the initial turbine discharge for full open
nozzle or turbine maximum discharge):
- the decrease of fL/D allows the discharge variation to be more
favourable;
- for large values of fL/D most of the discharge variation only occurs for
small values of nozzle opening, at the end of the nozzle manoeuvre.
Q/ Qo
H (m)
Pelton with 4 nozzles
1.0
0.8
0.6
0.4
f L/ D=0,1
f L/ D=5
f L/ D=15
0.2
f L/ D=10 0
0.0
0.0
0.2
0.4
0.6
% ope ning
0.8
1.0
Time (s)
Fig. 7.4- Discharge variation for different penstock characteristics and nozzle closures.
Head variation in the powerhouse for a typical manoeuvre (RAMOS, 1995).
The hydraulic conveyance circuit at upstream side of impulse turbines will be

influenced by the nozzle control system of the discharge variation imposed
through the closure law. This law and the hydraulic penstock characteristics will
be very important on the transient pressure analysis.
It is very important to note that it is no sufficient to define the total manoeuvre
time duration without knowing the conveying system characteristics.
7.4- Overspeed dynamic effects

7.4.1- Overspeed runner control
During normal steady-state operating conditions the speed of a turbine is
constant, but it will rise quickly when the load of the turbine is rejected. During
the closure of the nozzle or the turbine guide vane, the speed of the turbine will
rise to a maximum value that can not surpass the allowable one fixed by the
turbine manufacturer.
- 115 -
For impulse turbines the transients overspeed does not influence the pressure
along the hydraulic circuit. With reaction turbines the runner speed will change
the turbine discharge and the overspeed conditions will also influence the
hydraulic circuit (RAMOS, 1995). The transient overspeed can be a
waterhammer additional source, especially in what concerns the slow Francis
turbines (with small Ns values). A complete computer analysis should be
performed based on the turbine characteristic equations, the complete hydraulic
equations (7.12) and the equation of the rotating masses of each unit:
BH BR = I
d
dt
(7.15)
in which BH is the transient hydraulic torque (N m), BR is the resistant or

electrical torque (N m), I is the rotating masses inertia (kg m2) and is the
angular rotating speed (rad/s).
An inertia machine time constant (or start-up time of rotating masses), Tm, can be
defined based on equation (7.16):
Tm =
WD 2 n o2
x10 3
3575 Po
(7.16)
where no is the nominal runner speed (r.p.m.), Po the reference power or full load
turbine power (kW) and WD2 = 4gI (N m2). Tm has the order of magnitude of the
time for the unit to attain the speed no when submitted to a linear increasing
hydraulic power from 0 to Po.
For impulse turbines, and especially in what concerns the Pelton type turbine, the
runner overspeed control can be obtained by a deflector that deviates the flow of
the runner. This means that a slower closure manoeuvre can be specified in order
to better control the transient pressures along the penstock. For reaction turbines
the simultaneous transient pressure and overspeed control is much more difficult
(RAMOS, 1995). In order to evaluate, in a preliminary analysis, what is the
turbine maximum relative overspeed after a full load rejection some approximate
formula based on equation (7.16) can be used (LEIN, 1965 and HADLEY,
1970):
Lein formula
- 116 -
kT
H C
n
1
= 1 + C 1 +
no
Tm
H o
(7.17)
valid for n < 0.5 no. For Pelton turbines with jet deflectors, k=0.9, H is the
upstream waterhammer or transient head variation due to the manoeuvre and TC
is the time required by the deflector to change the flow direction, including the
dead time (TC 1.5 s). For Francis turbines, k=0.8 and TC is the guide vane
closure time.
Hadley formula
kT
H C
n
= 1 + C 1 +
no
Tm
Ho
3/ 2
(7.18)
where k is now related to turbine characteristics as a function of the specific

speed Ns, as indicated in the following table, and TC is the guide vane time of
closure:
Ns
(r.p.m.)
(kW, m)
76
134
286
381
573
672
763
0.98
0.96
0.84
0.77
0.66
0.61
0.57
Typically n/no should not exceed 0.6. For normal runner speed governing,
Hadley estimated that the following relationship between the minimum Tm and
TW (see eq. (7.10)) should be obeyed:
Pelton turbines Tm = 2.5 TW

Francis and Kaplan turbines Tm = 3.0 TW
The overspeed control and speed regulation stability in small hydroplants, during
normal operation with load demand changes, can be improved by increasing the
Tm value through the WD2 value (flywheel effect) or by decreasing TW value
(e.g. by inserting a surge tank or by selecting a larger penstock diameter). The
advanced electric speed-load regulators can now guarantee a better regulation
stability criterion.
- 117 -
Values of WD2 and TC should be guaranteed by the turbine and generator

manufacturers and are agreed upon at design stage. Another effective way to
control the transient overspeed and pressure in hydrosystems equipped with
Francis turbines is to equip the spiral case with special pressure regulators or
water outlets which automatically open and discharge water (as a relief valve)
into the tailwater pool when guide vanes are closed too rapidly (synchronous
relief valves). If none of these procedures is enough to control the transient
overspeed and pressure then another special protection device should be
considered like a surge tank as a last resort.
In small hydropower plants connected to a large national electric grid, the
stability of speed regulation does not present any special problem, because the
grid has the capacity to stabilise any unstable behaviour.
In exceptional conditions the turbine runner can be forced to accelerate under the
hydraulic power until a maximum limit steady state speed or runaway speed nrw
is attained in 3 or 5 seconds. This maximum speed will depend on turbine Ns and
the head values and varies typically between 1.5 no (propeller turbines) and 2.0
no (Pelton turbines) to 3.5 no (Kaplan turbines) for design head Ho (see Table
7.2). For other values of H the runaway speed can be estimated as follows:
*
n rw = n rw
H*
Ho
(7.19)
According to the BUREAU OF RECLAMATION, 1976, the runaway speed can

be estimated by this empirical formula:
n rw = 0.63N1s/ 5 n o
H*
Ho
(7.20)
where H* is the turbine head, Ho is the net head for the best efficiency regime and
no rotational speed.
Approximate formula must be used with care, as well as the estimation of the
WD2 or I values for turbine and generator. For a detailed analysis and final
design stage, computer simulations should be made based on the real selected
machine characteristics in order to be possible to analyse the complete dynamic
behaviour and interaction between the system components.
- 118 -
7.4.2- Overspeed effects on turbine discharge

For impulse turbines (e.g. Pelton turbines) at a constant needle position, the flow
is independent of the runner speed. It is usual to operate the nozzle and the
deflector simultaneously (double regulation turbine) during a closure operation.
The deflector moves the jet away from the runner and the nozzle closure law has
the main influence in the discharge regulation and penstock flow variation. For
Francis turbines with low specific speed, the flow drops with the transient
overspeed a constant gate position. Conversely, for Francis turbines with high
specific speed and for Kaplan turbines the transient discharge tends to increase
(RAMOS, 1995).
On Table 7.2 and Figure 7.5 is shown the turbine discharge reduction for low
specific speed Francis turbines.
Table 7.2 Runaway conditions for different types of turbines

Turbine
Type
(m, kW)
Ns
Pelton
Turgo
Cross-flow
Francis
Kaplan
15-60
20-70
20-80
80-300
200-800
Normal speed
no (r.p.m.)
Runaway
speed
nRW/no
Runaway
discharge
QRW/Qo
500-1500
600-1000
60-1000
500-1500
75-150
1.8-2.0
2.0
1.8-2.0
1.8-2.2
2.0-3.2
0.45-1.10
0.80-1.90
The flow across a runner, as former mentioned, is characterised by three types

of velocities: absolute velocity of the water V, with the direction imposed by the
guide vane blade, relative velocity W and tangential velocity C of the runner.
For uniform velocity distribution assumed at inlet (section 1) and outlet (section
2) of the runner and by application of Euler's theorem it enables obtaining the
equation, which relates the motor binary and the momentum moment between
these two sections:
BH = Q(r1V1 cos 1 r2 V2 cos 2 )
(7.21)
being a and r are the angle and radius (see Figure 6.6).
The output power in the turbine shaft is a result of multiplication the binary by
the angular speed .
- 119 -
P = BH
(7.22)
where
BH is the hydraulic torque and P the output power.
It yields in the following equation, after some transformations (RAMOS, 1995),
Q=A
+Bn
n
(7.23)
with
A=
60gH o
2 r2
1
1
+
2 b o r2 tan o A 2 tan 2
and B =
2 r2
60
1
1
+
2 b o r2 tan o A 2 tan 2
(7.24)
that gives a relation for turbine discharge depending upon the rotational speed
value and the characteristic of the runner.
The subscript o denotes outlet from the wicket gate, 1 and 2 denote,
respectively, inlet to and outlet from the runner, bo is the runner height (or freearea), r is the radial distance, is the angle that the velocity vector (V) makes
with the rotational velocity (C), is the angle that runner blades makes with the
C direction and A2 is the exit flow cross-section area.
Looking to this equation, the discharge regulation can be obtained by variation
of bo, o or 2. The bo variation is not easy to obtain, because this parameter is
a fixed characteristic of a runner (related to the height of the runner). Thus, for a
constant rotational speed, the discharge variation can be obtained by the
following procedures:
1 - variation of o;
2 - simultaneous variation of o and 2;
3 - variation of 2.
- 120 -
Procedure 1 is used for discharge regulation of Francis and propeller turbines

through the guide vane; procedure 2 is used for Kaplan turbines and procedure 3
is used for Kaplan turbines with a fixed wicket gate.
QRW /Qo (%)
120
110
100
90
80
70
60
DUM ONT
50
line a r co rre la tio n(b=22; m =0,245)
B OUVIER
40
70
100
130
160
190
220
N s (m,kW )
250
280
310
A Overspeed effect on discharge variation of reaction turbines (RAMOS, 1995).

Q
Q
% of nozzle opening
Q
% of wicket gate opening

100%
80%
100%
80%
60%
H=ct
Pelton
100%
80%
60%
60%
H=ct
100%
H=ct
n
Francis with low

sepecific speed
80%
60%
H=ct
n
Francis with medium

sepecific speed
Francis with high

specific speed
B Discharge variation with the percentage of guide vane opening.

Pump-Turbine Ns = 230 (m, kW) - Suter parameters
Pump-Turbine Ns = 65 (m, kW) - Suter parameters

-2.0
-2.0
Q/QR
Q/QR
-1.5
-1.5
BR =0
-1.0
-1.0
Turbine zone
Turbine zone
-0.5
0.0
0.5
runaw ay conditions
h=0.7
h=0.8
h=0.9
h=1
h=1.1
h=1.2
h=1.3
h=1.4
runaway
-0.5
runaway
runaw ay conditions
h=0.7
h=0.8
h=0.9
h=1
h=1.1
h=1.2
h=1.3
h=1.4
BR=0
0.0
0.5
1.0
0.0 -0.1 -0.2 -0.3 -0.4 -0.5 -0.6 -0.7 -0.8 -0.9 -1.0 -1.1 -1.2 -1.3 -1.4 -1.5
1.0
0
N/NR
-0.2
-0.4
-0.6
-0.8
-1
-1.2
-1.4
-1.6
-1.8
-2
-2.2
N/NR
C Example of a low and medium specific pump-turbine speeds.

Fig. 7.5 Turbine discharge variation analysis as a function of the runner speed.
- 121 -
-2.4
7.4.3- Turbine overspeed effects on waterhammer

For low specific speed reaction turbines (Francis turbines) or QRW/Qo < 1 (where
QRW is the turbine discharge at runaway conditions and Qo the nominal turbine
discharge), the overspeed effect provoked by runaway conditions will potentially
induce greater overpressures than those induced by just the guide vane closure
effect and it can be obtained through chart of Figure 7.6 based on systematic
computer simulations (RAMOS, 1995 and RAMOS; ALMEIDA, 1996).
In hydropower plants with long penstocks (or large TW / ((Q RW Q o )Tm ) values)
severe waterhammer troubles can occur. As presented on Figure 7.6, the
maximum upstream transient head variation will depend on Ns and two other
parameters: TW/Tm and TC/TE.
Ns
/Ho
____
HH+M0.45
Hb
0.40
T w/ T m
(m, CV)
330
310
290
0.35
Ns
0.30
270
250
0.25
230
T c/ T E=4,5
0.20
T c/ T E=6,7
0.15
210
190
T c/ T E=12,5
0.10
170
0.05
150
0.6
0.7
0.8
0.9
T w/ T m
2,99
1,50
0,99
2,29
1,15
0,76
1,38
0,69
0,46
1
QRw/Qo
variation
with Ns
Fig. 7.6 Maximum overpressure induced by both the overspeed and gate effects
of low specific speed reaction turbines on upstream penstock (RAMOS, 1995).
As indicated in the figure (follow the example through the arrows), the
calculation begins with the Ns turbine value. Moving horizontally the Ns dashed
line is reached and the QRw/Qo can be known. Knowing the TW/Tm value (the
relative water and turbine inertia time constants) and selecting the relative wicket
closure time TC/TE, the relative maximum upsurge variation can be obtained. For
QRw = Qo, the overpressure will only depend on the gate effect. These results
allow an approximate prediction of the maximum overpressure due to a sudden
load rejection in Francis turbines of a small powerplant.
- 122 -
When the powerhouse is equipped with reaction turbines, both the upstream and
downstream water columns must be analysed, in order to avoid excessive
overpressure at penstock and at powerhouse or even water column rupture at
draft tubes. A classic solution to minimise hydraulic transients or to reduce the
total length of the hydraulic conveyance system under waterhammer action
consists to insert a surge tank as near the powerplant as possible. However, this
solution can not be advisable, because it can be a very costly structure and can
cause significant environmental impacts for a small hydropower plant. Thus, a
good estimation of the maximum transient pressure without any special
protection device is very important.
7.5- Special protection devices

7.5.1- Introduction
The protection devices to control transient pressures due to turbine operation
can act by different ways: by supplying or removing water, or by flow energy
accumulation or dissipation. Typical protection devices for small hydropower
plants are among others (RAMOS, 1995):
- surge tanks;
- air vessels;
- synchronised valves, automatic discharge valves or pressure relief
valves;
- flywheels.
The mathematical formulation applied in computational modelling and in the
analysis of protection components of a hydropower scheme are presented. The
selection of the best solution for the protection devices depends upon the
hydraulic characteristics and topography that can influence the conveyance
system profile of the hydraulic circuit. Also important is the simulation of the
integrated system in order to know the best and correct dynamic response and the
influence of different devices always kept up with an economical comparison
study. Suitable protection devices to small hydropower plants are next presented
based on a simplified mathematical formulation and a brief explanation of its
performance.
For abnormal conditions (emergency or exceptional operating conditions)
creating potential excessive overpressures during a full load rejection of very fast
flow stoppage due to a safety valve closure, a pressure fusible section can be
- 123 -
used. This protection device is based on safety discs or membranes that will be
specified to rupture whenever pressure rises above a pre-set value. A battery of
discs can be selected for different pre-set pressures. The disc burst will allow an
outflow as a relief valve.
7.5.2- Surge tanks

Surge tanks allow the attenuation and control of rapid variations of discharge
and pressure, through storing excess energy (water volume) in an open
reservoir. During a transient, a surge tank will act as a large reservoir, where the
elastic pressure waves are supposed to be totally or partially reflected. Thus, the
penstock length submitted to waterhammer can be reduced to the length
between the surge tank and the powerplant. In this way TE and TW values can be
drastically reduced.
The unsteady flow between the upstream intake and the surge tank will have a
mass flow oscillation type, without elastic effects, as a rigid liquid column
oscillation between two open reservoirs. The singular headloss at the surge tank
connection side will dampen the water level oscillations.
Surge tanks can have different shapes and be of different types (Figure 7.7).
Fig. 7.7 Different types of surge tanks.

Ach
Tank
Qg
Tunnel or gallery
Zch
Qch
Qc
Penstock
Fig. 7.8 Typical scheme of a surge tank.
Equations used in a computational modelling of a simple surge tank (Figure 7.8)

are presented as follow:
- 124 -
- Continuity equation at the surge tank connection section or node

Q g + Q ch = Q c
(7.25)
- Kinematic equation (rigid model without considering inertia of the surge tank water)
dZ
(7.26)
A ch ch = Q ch
dt
- Head loss equation in the surge tank connection section
(7.27)
H g ,c = Z ch K ch Q ch Q ch
- Hydrodynamic characteristic equations
C+: H P H A + B(Q P Q A ) + R Q A Q A = 0 (gallery)
C-: H P H B B(Q P Q B ) R Q B Q B = 0 (penstock)
(7.28)
(7.29)
being:
B = parameter of characteristic lines of method of characteristics (MOC)
c
B =
, in which c = wave celerity (m/s) and A = pipe cross-section area
gA
(m2);
H g ,c and HP = piezometric head at surge tank insertion and any conveyance
system section (m);
K ch = singular head loss coefficient at surge tank insertion;
Qc = discharge in the penstock (m3/s);
Qch = discharge in the surge tank (m3/s);
Qg = discharge in the gallery (m3/s);
QP = discharge at any conveyance system section (m3/s);
R = friction coefficient of the gallery and penstock (R = Jx/Q2);
Z ch = water level in the surge tank (m);
A, B, P left, right and device section, respectively.
7.5.3- Differential surge tank

A differential surge tank is similar to a simple one, but has a much better
performance in what concerns the overpressure attenuation capacity (Figure
7.9).
It is provided by an eccentric or symmetric orifice connected to the main pipe
(i.e. gallery or penstock) and an interior pipe or riser that discharges into the
external or outer tank (RAMOS, 1995).
- 125 -
At
Tank
Ach
Zt
Qdes
Qch
Qor
Qc
Tunnel or gallery Qg
Penstock
A
Fig. 7.9 Differential surge tank. A scheme with notation at left

and real example under construction.
The computational modelling of this element is based on the following set of

equations (RAMOS, 1995):
- Continuity equation at the surge tank connection section or node
Q g + Q ch + Q or = Q c
(7.30)
- Continuity equation of the interior pipe or riser
Ach
dZ ch
= Qch + Qdesc
dt
(7.31)
- Discharge law of the weir at the top river section

3
2
Q desc = C vd Ldesc 2g H desc
(7.32)
- Continuity equation of the surge tank

dZ t
At
= (Q or Q desc )
dt
- Head equation of the orifice
H g ,c = Z t K or Q or Q or
- Head equation at the insertion node of the differential surge tank
1 dQ ch
H g ,c = Z ch
+ J ch ( Z ch Z bch ) K ch Q ch Q ch
gA
dt
ch

C-: H H B(Q Q ) R Q Q = 0 (penstock)
P
- 126 -
(7.33)
(7.34)
(7.35)
(7.36)
(7.37)
being:
B = parameter of characteristic lines of the method of characteristics (MOC)
c
B =
which c = wave celerity (m/s) and A = pipe cross-section (m2);
gA
C vd = discharge coefficient of the weir;

g = gravity acceleration (m/s2);
H g ,c = piezometric head at the surge tank connection section (m);
HP = piezometric head at any conveyance system section (m);

Jch = hydraulic grade line of interior pipe of differential surge tank;
K or = singular head loss coefficient in the orifice;
K ch = singular head loss coefficient at the tank connection section;
L desc = perimeter of the riser pipe of the differential surge tank (it works as a weir)
(m);
Qc = discharge in the penstock (m3/s);
Qch = discharge in the interior pipe (m3/s);
Qdesc = discharge over the weir (m3/s);
Qg = gallery discharge (m3/s);
Qor = discharge through the orifice (m3/s);
R = friction coefficient of the gallery and penstock (R = Jx/Q2);
Zch = water level in the riser pipe (m);
Z t = water level in the outer tank (m);
A, B, P left, right and device section.
7.5.4- Air vessel

An air vessel has a similar function to a surge tank, but it is closed and a smaller
reservoir, with entrapped air inside avoiding, thus, too large reservoir
dimensions due to the air compressibility absorption (Figure 7.10).
VRAC
air
ZRAC
QRAC
QA
QB
A
Fig. 7.10 Air vessel scheme.
- 127 -
The volume of air due to its compressibility effect will contribute for the
overpressure attenuation. An energy dissipater can be inserted in the connecting
pipe in order to better control the maximum overpressures and to dampen their
oscillations.
The computational modelling of this element can be based on the following set
of equations:
- Continuity equation at the air vessel connection section or node
Q A = Q RAC + Q B
- Head equation at the air vessel connection section
H RAC = H A + H a z RAC K RAC Q RAC Q RAC
- Polytropic equation for perfect gas behaviour
H RAC nRAC = C
- Continuity equation inside the air vessel
dZ
A RAC rac = Q RAC
dt
C-: H H B(Q Q ) R Q Q = 0 (penstock)
P
(7.38)
(7.39)
(7.40)
(7.41)
(7.42)
(7.43)
being:
2
ARAC = cross-section area of the air vessel RAC (m );

c
B =
which c = wave celerity (m/s) and A = pipe cross-section (m2);
gA
Ha = atmospheric pressure (= 10.33 m);

H RAC = piezometric head at the air vessel insertion (m);
K RAC = inlet and outlet singular head loss coefficient of the air vessel (RAC);
n = polytropic coefficient for the gas inside the air vessel;
3
Q RAC = discharge in the interior of air vessel (m /s);
R = friction coefficient of the penstock (R = Jx/Q2);
Z RAC = water level in the air vessel (m);
- 128 -
RAC
= air volume inside the air vessel RAC (m3);

A, B, P left, right and device section.
7.5.5- Synchronised valve or relief valve

These valves allow the exit of the flow, within a control discharge law and,
consequently, the attenuation of maximum upsurges by increasing the effective
flow stoppage time (Figure 7.11).
Pressure relief valve
Fig. 7.11 Typical installation of a pressure relief valve upstream a turbine.
These valves will be placed near each unit of the powerplant and will
automatically open to the atmosphere or outlet pipe if a fast closure of the
wicket gate occurs. The relief valve will open during the turbine fast gate
operation (synchronous valve RAMOS, 1995). After the end of the wicket
gate closure, the relief valve will slowly close. Other types of relief valves will
open as soon as the penstock pressure exceeds a limit value. A relief valve is an
adequate device for high heads (typically Ho > 50 m).
The computational modelling of this element is based on the following set of
equations:
- Turbine characteristic equations
- Characteristic equation of the turbine and valve (when the valve is coupled to the
turbine)
Q V = ( C Valv ) H o
(7.44)
- Hydrodynamic characteristic equation C+ of the penstock

H P = C1 BQ
and a tailrace at constant water level ZJ:
H P = H o + ZJ
where
- 129 -
(7.45)
(7.46)
c
B =
which c = wave celerity (m/s) and A = pipe cross-section area (m2);
gA
C1 = parameter of characteristic lines (MOC) (C1 = H A + (B R Q A )Q A ) ;
CValv = valve discharge coefficient that will change as a function of the valve closure
or opening laws (-);
Ho = net valve head (m);
QV = valve discharge (m3/s);
R = friction coefficient of the hydraulic circuit (R = Jx/Q2).
7.5.6- Flywheel
The flywheel, whose function is to accumulate energy by increasing the rotating
mass energy of turbo-generators, allows the increase of the stoppage time TC of
the units as well as the time to attain the runaway speed (Figure 7.12).
Flywheel
Generator
Turbine
Fig. 7.12 Typical turbine flywheel installation.
With this device, the turbine runner overspeed and the transient pressure will be
both controlled.
The computational modelling of this element is based on the following equation
that should be solved with the turbine characteristics and the hydraulic
equations:
-Rotating mass equation
- 130 -
BH BR = I
2 dn
60 dt
(7.47)
in which
I = Itur + Ige + Ifw
(7.48)
where
BH = hydraulic torque (N m);
BR = resistent electromagnetic torque (N m);
I = total rotating mass inertia (kg m2);
Itur = rotating mass inertia of the turbine (kg m2);
Ige = rotating mass inertia of the generator (kg m2);
Ifw = rotating mass inertia of the flywheel (kg m2);
n = rotating speed of the unit (r.p.m.).
It can be concluded that the flywheel computational analysis just implies an
increase of the total rotating mass inertia, I, of each unit in the model. With this
increase Tm will also be increased, as well as the starting up time and the
regulation stability conditions will be improved, especially in what concerns the
isolated grid operation. On the other hand the increase in Tm will diminish the
transient maximum runner overspeed.
7.5.7- Protection devices behaviour

7.5.7.1- Analysis of a surge tank
Generally, surge tanks are installed at downstream of a diversion gallery (or
tunnel) or at upstream of a long tailrace tunnel, in order to reduce the total length
of the hydraulic conveyance circuit under waterhammer effects (due to elastic
waves) and to improve the turbine speed regulation stability should the plant be
connected to an isolated electric grid.
In a surge tank analysis (Figure 7.13), it is usual to consider the mass oscillation
(or rigid water column) regime between the upstream reservoir and the surge tank
or between the surge tank and the tailrace when the surge tank is placed
downstream the powerplant (RAMOS, 1995).
In some cases, due to construction and economic constraints, protection devices
are not directly connected to the gallery or pipe and, sometimes, it is necessary to
have a long connection pipe. In practice, at small power plants this pipe does not
have a significant influence, since the pipe length is not too much long. However
- 131 -
a surge tank may be not able to reflect the most severe incident elastic waves
coming from the penstock. In this situation a transmitted wave with reduced
intensity will propagate along the protected tunnel or gallery. To avoid this
wave transmission the surge tank connection should not have a very strong
restricted linking pipe.
Fig. 7.13. Surge tank effect on the reflection and transmission of the pressure elastic
waves induced by turbine discharge variation.
Some sensitivity analysis carried out concerning the gallery and the penstock
cross sections allow concluding that the main influence is obviously imposed
by the gallery cross section area that as smaller as greater is the water level
attained inside the surge tank (RAMOS, 1995).
Fig. 7.14 Two typical surge tank schemes: simple cylindrical tank with a restricted
connection and a differential tank.
Based on the main oscillation theory, analytical solutions can be obtained for
the calculation of the extreme surge tank water levels for simplified conditions
- 132 -
Hydraulic Transients and Dynamic Effects
as, by example, for simple (cylindrical) tanks and the complete and
instantaneous flow stoppage at the powerplant. Neglecting the friction and
singular head losses, the maximum water level oscillation Z* in the surge tank
will be:
Z * = Qo
L
LA
= Vo
gA Ach
g Ach
(7.49)
where
Qo = initial gallery or tunnel discharge (m3/s);
Vo = initial flow velocity corresponding to Qo (m/s);
L = gallery length (m);
A = gallery cross section area (m2);
Ach = surge tank cross section area (m2);
g = gravity acceleration (m/s2).
The water level period of oscillations T*, in these simplified condition, may be
obtained by
T * = 2
LAch
g A
(7.50)
It can be easily concluded that this period is much longer than the pressure
elastic fluctuations (4L/c). The water level oscillations will be damped by the
friction and singular headlosses along the gallery. The initial headlosses will
also modify the maximum water level fluctuation Z*. This influence can be
evaluated by the following approximate formula deduced by Jaeger:
2
2
1 H o
Z1 = Z * H o +
3
9 Z *
(7.51)
where
Z1 - deduced water level oscillation;
Z* - theoretical water level variation without headloss effect;
Ho - gallery total headlosses for the initial discharge Qo.
Another water level oscillation reduction can still be considered as a function
of the time of the turbine discharge change or time of gate closure (e.g. for
TC = 0.4 T* and Z 0.75 Z*).
- 133 -
A longer closure time for the wicket turbine can be feasible for Francis turbines
if a synchronous relief valve is installed. For Pelton turbines the deflector will
deviate the jet and the transient overspeed will be controlled.
In order to reduce the water level variation in the surge tank some special
energy dissipative devices or flow restrictions are sometimes placed in
connecting pipe (e.g. an orifice or diaphragm). However, a strong flow
restriction will increase the transmission of the waterhammer waves for the
gallery.
One of the classic problems to be considered in a surge tank design and
specification is the stability of the water level oscillations, should the
hydroplant be under the action of an automatic regulator with feedback
imposing a turbine discharge variation through the turbine gate position, as a
function of the continuous variation of the load power demand and the net
head. This is a typical situation found in an isolated electric grid.
The surge tank stability analysis depends on the overall hydro and electric
(demand) system behaviour. An approximate analysis can be based on very
simplified conditions and small disturbances and a linear method of analysis of
the set of equations. A minimum cross section area, or Thoma section, is
theoretically obtained for simple cylindrical tanks in order to guarantee the
stability of the water level oscillations:
ATh =
L A V2
H o H o 2 g
(7.52)
where
Ath = Thoma cross section area;
A = gallery cross section area;
L = gallery length;
V = gallery flow velocity;
Ho = net head;
Ho = gallery headloss.
This limit area can impose a very costly surge chamber. Formula (7.52) can be
improved with minimum cross section areas reduced to 50 - 60%, should the
dimensions and cost, as well as the powerplant importance, justify it by taking
into consideration several additional factors (GARDEL, 1956).
- 134 -
Computer simulations with the complete non-linear equations will allow to

consider the influence of the entire main factors, as well the large water level
oscillations, in the stability and operational analysis.
There are a lot of types of special surge tanks. Each one has some advantage in
what concerns the hydraulic response to transient regimes, the local technical
constraints (e.g. geological and topographical), the environmental impacts or
the construction costs. Most of the types of surge tanks are not adequate to
small hydropower systems because they are too complex and costly.
One exception is the differential surge tank (e.g. the Johnson type) with an
internal pipe or riser and an outer tank. When the flow reaches the top of the
internal pipe, after a turbine discharge stoppage, it will spill into the outer tank.
The rapid head increase in the riser will create a faster deceleration head
applied to the gallery flow. In order to prevent excessive surge amplitude, the
height of the riser is limited and the outer tank will store the spilled water.
A load demand increase will induce a rapid drop in the riser with a forcing
acceleration effect upon the gallery flow.
The water required to supply the increased turbine demand is drawn from the
storage outer tank through one or more openings.
The savings in volume of this type of surge tank over that of a single chamber
with a restricted - orifice can be obtained if the critical turbine flow demand is
smaller than the critical rejected flow (e.g. when the powerplant has multiple
units).With a differential surge tank including an eccentric extra orifice (see
Figure 7.14), the elastic wave transmission (S) can considerably be reduced.
Relatively to the Johnson type surge tank (Figure 7.15), especially in what
concerns the wave transmission a reduction of S = 5% can be obtained
depending on the relative dimensions of the orifice diameter (Dorif) and of the
tank diameter (Dtank).
15
14
13
12
S (%) 11
10
9
8
7
Johnson
eccentric orifice
0.1
0.15
0.2
0.25
c=(Dorif/Dtank)2
0.3
0.35
Fig. 7.15 Example of comparison between a Johnson surge tank

and eccentric orifice surge tank (RAMOS, 1995).
- 135 -
Some conclusions can be presented according to this analysis:

1 - the differential surge tank with eccentric orifice seems to be more
efficient in wave transmission;
2 - the internal tank cross-section must be carefully studied because it
can reduce or induce variations in the gallery wave transmission
factor;
3 in general, there is an increasing in the gallery wave transmission
when the outer tank cross-section is reduced;
5 both surge tank types have less environmental impact than a normal
surge tank due to have smaller global volumes.
More information about methods of analysis and the behaviour of different
types of surge tanks can be found in JAEGER, 1977 and MOSONYI, 1991.
7.5.7.2- Analysis of an air vessel

The use of air vessels or air-cushion surge chambers is well known as an
efficient protection device in pumping systems. Air vessel has also been used in
large high head hydroelectric schemes through rock caverns (e.g. in Norway).
In small hydroplants this device is not frequently selected. However, the air
vessel can theoretically replace an open surge tank with much less volume
(Figure 7.16).
Alb u feira
ch a m in
Reservoir
Surge
co n d u ta
tank
fo ra d a
Penstock
tnel
RAC
Tunnel or
gallery
Air
vessel
Restitu i o
Tailrace
Fig. 7.16 - Comparison of surge tank with air vessel layouts.
- 136 -
The range of water level oscillations depends on selection of the air-filled

volume of the vessel or chamber, as well as on the singular headlosses created
at the connecting pipe.
Several simulations were made for a specific system to enable the analysis of
an air vessel behaviour based on the analysis of the singular head loss
coefficient at the inlet (Ke) and outlet (Ks) of the vessel, and its influence in the
surge attenuation (Fig. 7.17). Sensitivity analyses were also developed, in what
concerns the air volume, the coefficients Ks,e and the gallery length variations
(Fig. 7.18). Longer circuits are unfavourable for the overpressures, as well as it
is known, and Ks has more influence on minimum pressure values and Ke on
the maximum overpressure values. It is necessary to have some care with these
values because they can also increase the overpressures when the constraining
is too high.
H/Ho
1.0
0.8
0.6
0.4
0.2
0.0
-0.2
10
12
14
-0.4
being a
characteristic
parameter for
air vessels:
-0.6
-0.8
=
= f (Vo)
oc
nAVoL
Fig. 7.17 Extreme variation of the piezometric head by increasing the air volume
inside the air vessel (RAMOS, 1995)
The air vessel dimensions or volume should also satisfy the stability constraints
similar to open surge tanks. A modified Thoma critical section ATha was
presented by Mosonyi for preliminary studies (Svee formulae):
p
ATha = ATh 1 + n o
ho
(7.53)
where
n - polytropic coefficient for air thermal process;
- 137 -
po - reference air pressure inside the vessel; ho - air height inside the vessel;
- water volumetric weight.
3500
3000
2500
2000
Ks
Ke 1500
Ks
1000
Ke
500
0
0
10
20
Vo (m3)
30
40
50
Fig. 7.18 Optimisation of the head loss coefficient value

as a function of air vessel volume
It is always important to determine the maximum downsurge in consequence of

a sudden shut-down of the electric load demand in a standstill period.
7.6- Examples
In this section a graphically transient comparative analysis will be presented,
for a total pressurised penstock. Different alternative types of protection
devices are considered. The system dynamic response varies with the hydraulic
and equipment characteristics.
Figure 7.19 shows the more or less capability of each protection device to
attenuate the pressure variation along the conveyance system, and the time
duration to dampen the transient regimes and to establish again the operational
capacity system. Only through a technical and economic comparative analysis
will be possible to select the best solution for each case.
A sensitivity analysis was carried out in order to obtain the dynamic response
of a canal connected to a forebay, during the fill up followed by start-up and
stoppage of the turbo-generator units and manoeuvres in valves (Figure 7.20).
The knowledge of the dynamic behaviour of the complete integrated system is
fundamental to establish criteria to specify the walls height, the dimensions of
the canal, forebay and weir, the need of a PID type (or equivalent) regulator,
the characteristic time constants and automation parameters.
- 138 -

I
II
III
IV
V
VI
Fig. 7.19 - Sensitivity analysis regarding different protection devices for overpressure
control due to overspeed effect induced by a Francis turbine: I) without protection device; II)
with flywheel; III) with elastic pipe and expansion chamber; IV) with differential surge tank; V) with relief
valve by pressure control; VI) synchronised relief valve by overspeed effect
- 139 -
(RAMOS, 1995).
I
II
III
IV
Fig. 7.20 - Dynamic behaviour induced by the filling of a canal through the opening of
an upstream gate (with a forebay at downstream): I) without regulation and the turbine discharge
equal to canal discharge; II) without regulation and the turbine discharge greater than canal discharge
inducing the emptiness of the canal; III) with the actuation of a PID regulator to control downstream level at
0.90 m although the turbine discharge is greater than canal discharge; IV) with a PID regulator to control the
water level in 0.90 m although turbine discharge is smaller than canal discharge
(RAMOS, 1995).
7.7- Other protection devices and procedures

In small hydroplants with long hydraulic conveyance systems upstream or
downstream the turbines, or in both sides, the most feasible, reliable and
economic solution should be selected and specified.
- 140 -
For normal operations the maximum and minimum transient pressures along
the tunnels and penstocks should be compatible with specific safety factor in
order to avoid any pipe bursting, buckling effect or cavitation or vacuum
formation phenomena.
The designer should always consider the most adverse conditions. Some of
these situations correspond to the maximum turbine flow, others correspond to
smaller flow velocities as, by example, the full closing of the turbine gate with
speed no load condition: the minimum flow rate for nominal speed can be
nullified by a fast manoeuvre and the overpressure will be the one given by
Joukowsky formula.
To guarantee the global safety condition several alternative solutions are
typically considered during the design phases. The most common solution
involves the specification of the adequate laws of closure and opening of the
turbine flow control devices (e.g nozzles and guide vanes (or wicket gates)) and
of any safety valve (e.g. pressure relief valves).
The filling and emptying of pressurised pipes need also to be considered,
including the position of bottom (or purge) valves and strategically placed air
vents for air inlet and outlet.
In each selected case the extreme transient head envelopes along the hydraulic
circuit should be obtained and compared with the profile of the pipe system.
Abnormal conditions should also be envisaged according to the hydrosystem
characteristics, including the malfunction of any vital device or an operational
failure due to human error.
A large number of these situations involve the operation of safety valves as
the fast closing of the safety valve placed just upstream a turbine should the
normal flow control fail and the turbine runner velocity increase;
the fast closing of a safety valve due to an abnormal flow velocity increase
due to a pipe burst.
In these cases, a convenient closure law needs to be specified in order to

control the transient maximum and minimum pressures. This specification will
mainly depend on the following factors:
- 141 -
the characteristics of the pipe system to be protected, especially the pipe

factor Fp = fL/D as was already explained; in fact, the relative pipe
headloss can adversely modify the system behaviour and the same valve
closure time can induce a slow or a rapid flow change according to Fp
value;
the intrinsic characteristics of the valve; a butterfly valve (for medium
heads) and a spherical valve (for high heads) have different effects on the
flow for the same closure law;
the type of valve actuator.
Most of the fast closing hydroplant safety valves close under an external arm
action with a weight that is liberated by a mechanism as soon as the flow
velocity at a nearby pipe section exceeds a pre-fixed value. Under the weight
action the valve will tend to close very fast. A hydraulic damper will slow the
movement in the closing final phase. A valve closure law with two speeds is
obtained that can be considered as a bi-linear closure law with two fundamental
time parameters: the total time TF and the intermediate time, where the closing
gradient changes, TI.
For these types of valves this bi-linear closure law is also very adequate and
easy to model should the intrinsic valve characteristics or the valve discharge
equation as a function of the valve relative opening be known.
7.8- Integrated analysis and design

The integrated system response with the different components depends on the
disturbance type or excitation induced (Figure 7.21). The interaction of the
different components can induce potential accidents and, at limit, resonance
phenomena. As mentioned before, the type of turbines will strongly influence
the system response.
The hydraulic conveyance system behaviour will depend on several
components and interactions. In each case, especially in what concerns a small
hydroplant the design needs to select the most important disturbance sources to
be considered. A small powerplant connected to the national grid can justify
not to consider most of the stability criteria applied to isolated electric grids. In
all cases the integrated analysis should not forget some special situations that
can be modelled with less components, as by example:
- 142 -
penstock emptying or drainage by bottom valves and selection of air inlet

or outlet valves if it is considered necessary;
stability analysis of pipe anchors based on the extreme hydrodynamic
forces obtained from the transient analysis and the extreme pressure
envelope along the pressurised penstock.
Reservoir
Turbine
Hydraulic Circuit:
conveyance system
and outlet
Protection
Devices
Outlet
Generator
Hydromechanic
Equipment
Electric
Line
National
Grid
Automation
and
Regulation
Fig. 7.21 - Oscillatory characteristics of different components

of an hydropower scheme.
The control of hydrotransients and the dynamic behaviour of hydroelectric

systems are fundamental for the design and the exploitation, in order to
guarantee a safe and reliable solution. The effectiveness of new design criteria
based upon computational techniques should allow the analysis of the global
behaviour in order to identify eventual operational constraints, since the
beginning of the exploitation, as well as particular local and partial situations
that are vital for the system safety.
7.9- Case studies

In the Table 7.3 are presented the main characteristics of eight small
hydropower systems, in the North of Portugal. All schemes have long hydraulic
conveyance systems. One of them has a differential surge tank and two are
equipped with relief valves. The schemes equipped with Pelton have no special
protection devices against transient pressures. Two schemes with Francis
turbines have not also any special protection devices. The S. Pedro do Sul
powerplant is a very special case with TW = 8.5 s and QRw/Qo = 0.9. In the
Torga scheme, field tests (load rejection) were performed. The Francis turbines
- 143 -
have Ns = 130 (m, kW) and and estimated QRw/Qo = 0.6. Comparing with S.
Pedro do Sul this scheme has smaller TW value but greater discharge reduction
by overspeed (runaway condition). In fact, the overpressure at S. Pedro do Sul
is due to the wicket gate (or guide vane) closure and at Torga is due to transient
turbine overspeed.
Table 7.3 Case studies. Comparison between eight small hydropower
plants with different characteristics
Identification
L (m)
Hg
Qo
Po
Ermida
3900
Ovadas
Turbines
Prot.
TW
(m)
(m3/s)
(MW)
395
2.35
7.7
Pelton (1)
None
2.5
3025
334
2.15
5.9
Pelton (1)
None
2.5
Sordo
3600
321
3.6
9.8
Pelton (2)
None
2.6
Torga
1200
60
20
9.1
Francis (2)
None
3.2
Nunes
2500
107
12
9.9
Francis (2)
Relief
3.2
device
valve
V. Viosa
1935
123
3.59
3.6
Francis (2)
Relief
4.3
valve
S. Pedro Sul
4000
74
14
8.2
Francis (2)
None
8.5
Terragido
1575
125
10
10
Francis (3)
Dif. sur.
2.0
tank
100
(%)
80
Overspeed ( n/no)
60
Overpressure ( H/Ho)
40
Two units operating
20
Wicket gate(s) closure (s)
0
2.2
3.4
4.0
4.0
4.5
4.6
5.0
5.0
6.0
7.0
8.0
9.0
9.5
P ot
ncia s (MW )
Power
Fig. 7.22 Full load rejection of Torga units. Field tests of overpressures, overspeeds
and wicket gate(s) time closures (RAMOS, 1995).
- 144 -
__
Electrical Equipment
ELECTRICAL EQUIPMENT
H. Pires de Almeida
8.1- Generators
8.1.1. Synchronous generators
The generator is a rotating machine, with its shaft coupled to the turbine,
providing the conversion from mechanical to electrical power.
There are two main types of generators used for this purpose, the synchronous
and the asynchronous.
The synchronous generators have alternating current in the three-phase armature
windings (normally in the stator) and direct current in the field winding
(normally in the rotor). With the rotor at synchronous speed, its direct current
field flux establishes a rotating field on the stator travelling at equivalent speed,
hence at the nominal system frequency (50 Hz).
Its speed and the number of poles define the generator operating frequency.
n = f/p
(8.1)
with
n - speed in revolutions per second;
f - frequency in hertz;
p - number of pole pairs.
- 145 -
The equivalent circuit shown in Figure 8.1 represents these generators. The
airgapg flux and resulting voltage Vo depend on the field current. If an
automatic voltage regulator (AVR) which establishes the voltage under no-load
conditions and the power factor when the generator is connected to the grid
normally controls this current.
Fig. 8.1 Synchronous generator equivalent circuit.
In approximate terms, the following expression shows the influence of the field
current:
Vo = d / dt = = 2 f L If
(8.2)
with
L - field inductance;
generator angular velocity.
The generator is started-up unexcited and driven by the turbine from standstill to
nominal speed. The field excitation is applied at around 95% nominal speed and
Vo and V1 are then controlled by the AVR, through the If value.
Controlling speed and phase angle via the turbine wicket gate control and with
the AVR providing fine voltage adjustments does the manoeuvre of
synchronisation with the grid. The machine is synchronised by closing of its
circuit breaker when the generator and grid voltages are in phase and with the
same magnitude.
Once the generator is connected to the grid the system voltage imposes a
reference to its behaviour, which can be interpreted from the P/Q diagram shown
in Figure 8.2.
The voltages Vo and V1 are typically separated by an angle , called the load
angle. The angle shows the phase difference between voltage and current. The
cos is normally referred to as power factor. The complex power components
are:
Active power - P=V1 Ia cos
Reactive power - Q= V1 Ia sin
- 146 -
__
The load angle is an indicator of generator stability with respect to the grid.
The active power supplied to the grid is given by P=V1Vo sin / Xg, reaching
the steady state stability limit for = 90 (Figure 8.2). Further increases in load
angle only take place under transient stability conditions. In such cases, a
suitable response of the voltage regulators is necessary to avoid pole slipping.
Fig. 8.2 P/Q Diagram.

d /dt
0.15
0.1
0.05
0
-0.05
10
15
20
25
-0.1
-0.15
(degree)
Fig. 8.3 Typical response of a transitory stability of a generator. Linear variation of the
turbine guide vane opening degree from parallel situation (RAMOS, 1995).
Increases in load angle are due to accelerating power, Pa= Pmec - Pel = M
d/dt, where Pmec is P, turbine power, Pel is generator power, with energy
transfers between the electrical power and the machine rotational kinetic energy
(ER =1/2 I 2), where I is the polar inertia moment.
The AVR fast response provides enhanced transient stability by increasing If and
Vo. This allows for a transient increase in output power, reducing the load angle
and allowing the machine to re-gain stability (see Figure 8.2).
- 147 -
Figure 8.3 shows that transitory generator stability is assured for a final
inferior than 45o. In the same figure is visible the phase plan for an increasing of
opening degree of the guide vane (25% to 30%) in different time increments
(dashed line in 2 s and fill line in 10 s).
Under isolated grid operation Vo and V1 set up the output voltage and a precise
wicket gate control ensures appropriate speed control and approximately
constant frequency. In these operating conditions, the generating units are
provided with flywheels with sufficient polar inertia and kinetic energy to absorb
and provide transient power fluctuations due to load variation, compensating for
the wicket gate response times.
8.1.2 Asynchronous generators

These machines are also called induction generators, being magnetised from the
three-phase winding in the stator and inducing slip currents on the rotor.
The rotor revolves at a speed slightly above synchronism and the currents
induced on it are of a frequency corresponding to the speed difference.
Slip = (synchronous speed - rotating speed) / synchronous speed
Torque results from the magnetic interaction between stator and rotor fluxes and
varies with slip, as indicated on the torque / slip curve of Figure 8.4.
operating point
Fig. 8.4 Torque / Slip curve
These generators can only supply active power to the grid and must import
reactive power for their magnetisation.
The procedure for grid connection requires running up driven by the turbine. The
generator breaker is closed when the actual speed passes over the synchronous
speed, at very reduced acceleration, in order to limit circulating currents when
closing. On the grid connection the generator absorbs a transient magnetising
- 148 -
__
current of short duration. The turbine power is then increased at a suitable rate
by the wicket gate control, while the slip and stator current increase according to
the generator characteristics.
The asynchronous generator control and handling towards the grid is simpler
than the synchronous but its application in small power plants (or mini-hydro)
schemes is normally limited to 2 MW machines due to transient effects of grid
connection and the need to install capacitors for power factor correction.
8.2. Electrical installations

8.2.1. Main transformer
The typical single line diagram of a mini-hydro power station is shown in
Figure 8.5.
The transformer is a static unit with the purpose of stepping-up the generation
voltage to the grid connection voltage level.
The voltage transformation is achieved by having the suitable number of turns in
the primary and secondary windings. The voltage ratio is defined by
N1/N2=V1/V2= I2/I1
(8.3)
with
N1 and N2 - primary and secondary winding number of turns.
V1 and V2 - primary and secondary voltages.
I1 and I2 - primary and secondary currents (inverse ratio).
The transformers are normally of the immersed oil type, but dry types/resin
impregnated are also used for the lower powers.
The high voltage winding of the power stations main transformers are normally
fitted with tap changers to allow for expanded adjustment range to the grid
voltage.
Fig. 8.5 Typical single line diagram.
- 149 -
8.2.2. Switchgear
The mini-hydro power stations normally use the medium voltage switchgear,
switches and circuit breakers, housed in metal cubicles, these being standard
market products satisfying the European electrical regulations (CEI).
The circuit breakers operate in SF6 or vacuum also with standard ranges of
breaking capacity and open / close times.
8.2.3. Control equipment

The typical powerstation is equipped with electrical control and supply boards,
housing the manual and automatic control circuits and auxiliary equipment.
Switches and push buttons support the manual controls, with the operator being
informed about the state of the plant by the indicating instruments.
The automatic controls are based on a Programmable Logical Controller (PLC),
which receives on line information through transducers and digital input signals
and takes the necessary control actions for water utilisation, flow or level
control. The PLC output is processed via suitable relays.
The PLC is also used for data processing and transmission via telephone
modems or radio signalling.
8.2.4. Electrical protection

The main features of electrical protection are the generator, transformer and grid
connection.
Generator protection
The amount of protection relays depends on the machines rating, normally
being installed overcurrent, under-voltage, reverse power and stator earth fault.
The units above 2 MW may have also differential, loss of field and negative
phase sequence relays.
Transformer protection
The overcurrent relays are always installed. The differential relay is used in units
above 2 MW. Transformers are normally supplied with Buchholz, oil and
winding temperature relays and pressure relief valves.
Grid protection
These protections must ensure appropriate separation from the grid in case of a
fault. The minimum protection scheme is normally defined by the utility
standards.
- 150 -
__
Nowadays the tendency is to install the less important protection features, on the
small machines, as part of the PLC software. This option requires however some
careful risk evaluation and circuitry analysis, ensuring minimum reliability
levels.
8.3. Control system considerations

8.3.1. Introduction
The generating unit behaviour shows the effects of the hydraulic circuit, the
turbine and the generator dynamics.
Figure 8.6 shows the interaction of the various blocks.
The hydraulic circuit and turbine features have been dealt in previous chapters.
This chapter describes the way the control actions are carried out for speed,
water level, power output and voltage / power factor control.
Speed, water level regulation and generated power are controlled, as previously
mentioned, by the powerstation PLC using appropriate software. This software
is based on sampled data with a short (20 to 50 ms) scan cycle time. The
controller is continuously acquiring data at a fixed sample rate. This system is
very suitable to mini-hydro stations because the time constants involved in these
processes are very long when compared to the cycle time. Furthermore, the
regulation parameters are very accessible.
Fig. 8.6 Power conversion scheme.
8.3.2. Speed control

Taking a Francis turbine as an example, the wicket gate is operated by a
hydraulic servomotor, controlled by electro-valves activated by the PLC outputs.
On start-up the wicket gate opens to a pre-set value and as the unit gains speed
the wicket gate closes to the value needed to run at no-load around synchronous
speed. Acceleration must be very low as speed reaches the nominal value. A
percentage opening varying between 6% and 12%, depending on the turbine
characteristics ensures the no-load running. The precise speed regulation for
synchronisation takes place via an automatic synchroniser. The total time from
standstill to synchronisation is of the order of 3 to 4 minutes.
- 151 -
The initial stage of speed increase is done with proportional control impulses,
adding derivative and integral control when necessary during the closure
manoeuvre and final regulation to reach synchronous speed.
The automatic synchroniser is an analogue device fitted with proportional and
integral control loops.
8.3.3. Water level regulation

This type of regulation is used mainly on run-of-the-river schemes. For cases of
small forebays, it is necessary to have proportional and derivative control and
sometimes, integral controls if it is relevant to maintain a very precise water
level.
In other cases, with some storage capacity, it is sufficient to have proportional
control and slow actions.
Normally in this type of regulation, maximum and minimum water levels are
established to back up the normal control system and avoid excessive deviations.
8.3.4. Generator output power control

This is normally done when significant storage capacity is available and consists
in operating a powerstation at pre-established fixed power / flow values during
the more relevant hours of the day and making use of the best combination of
turbine and hydraulic efficiencies.
The wicket gate control is done by impulses until the pre-set generator output
power is reached.
As in the previous case, upper and lower limits are established to the water level.
In this type of water utilisation, the volumes of water to be used are calculated
daily to establish the operating power values.
8.3.5. Voltage and power factor regulation

This control is done over the generator excitation system.
Under no-load conditions the AVR controls voltage and assists the automatic
synchroniser, providing the correct voltage for synchronisation. When on-load,
connected to the grid, it changes to power factor regulation and controls the
reactive power exchange with the grid, according to the regulations.
The voltage regulation is still typically done by dedicated AVR`s, fitted with
Proportional, Derivative and Integral control loops (PID).
The parameters on this equipment are available for adjustment in terms of the
gains and time constants of each loop, as well as the set points.
- 152 -
__
The AVR settings are normally pre-set at the factory and adjusted
experimentally on site.
8.3.6. Flood level control

When flood control gates are installed at the dam, these should have a control
system of the water level regulation type. The gates open by the appropriate
amounts, in order to regulate the upstream level. The flood gate manoeuvres are
very slow, avoiding sudden changes of flow downstream.
8.3.7. Blade control of Kaplan turbines

This is also an example of application of Proportional and Integral control since
precision is necessary as it affects the turbine efficiency. For specific cases
derivative control may also be added.
8.3.8. Analytical representation

The PID control may be represented (Figure 8.7) by the following basic
expression:
t
e(t ) = K p Z(t ) + K i Z(t )dt + K d

0
d Z(t )
dt
(8.4)
where
e(t) is the relative variation of any parameter (ex: voltage deviation: Vref - Vg)
with respect to the reference;
Z(t) is the relative variation of other associated variable (ex: exciter field
current Ife) with respect to the regulation variable;
Kp; Ki;Kd are the proportional, integral and differential gains of the regulator.
Fig. 8.7 Typical response of a PID regulator (RAMOS, 1995).
- 153 -
Applying Laplace Transforms this may be represented under the block diagram
transfer function form, as indicated in Figure 8.8.
Fig. 8.8 PID control system. Example of AVR block diagram.
8.3.9- Switchyard
For plants with capacities greater than 600 kW the transformers and switching
gear should be located in a yard outside. This switchyard should be as close to
the powerhouse as possible and should be surrounded by a sturdy, durable fence
for safety (Figure 8.9).
Fig. 8.9 Grid substation to transform the energy produced from 6kV to 60kV and fed it
into the supply network of the national electric grid.
- 154 -
Environment
ENVIRONMENT
Helena Ramos
9.1- Environmental impact assessment

The purpose of an environment impact assessment (EIA) applied to small hydro
power plant consists in the evaluation of the favourable and unfavourable
impacts, in what concerns the two categories of impacts, in natural and social
environmental context. Natural impacts include hydrology and sediment effects,
as well as the water temperature and quality, ecology, engineering construction,
biology, landscape and effects on archaeology and cultural assets, soils and
geology, air, noise and, eventually, climate local change (in case of large
reservoirs). Social impacts involve social, cultural and economic development
inducing local industrialisation and changes in citizens life quality as well as
potential people displacement due to submersion by the reservoir (nevertheless
rarely applied at small hydro plants - Figure 9.1).
For an EIA, firstly, it is necessary to identify potential impacts based on intrinsic
characteristics of the small hydro plant and on the site, taking into consideration
the experience and the knowledge of impacts provoked by similar projects and
contributions from the technical team. According to the magnitude and the
importance of impacts, the prediction specific techniques, based on quantitative
and qualitative analyses, must be carried out, in order to classify them as
- 155 -
positive, significant or not, whenever occurs a violation of standards or rules

legally approved or accepted by the scientific community.
Being the water an essential resource that sustains all life on Earth, it is
important to plan its correct use. Nowadays it is obligatory to provide all
information about watercourses through suitable catchement management plans,
which includes several aspects, such as water abstraction, effluent disposal,
flood defence, amenity and recreation, pollution, land uses, flow variation and
wildlife conservation.
The flow in streams and rivers is, normally, the main outflow component from a
catchement, and the river damming even by a low dam constitutes a major
concern of populations (when not well informed).
Social and
economic impacts
Natural impacts
engineering
construction
employment effects
hydrology
access improvement
agriculture
water quality
expenditure on local
assets
geology and land

uses
development of
monitoring systems
biology
EIA
effects on local
services
sedimentation
human health
air and noise
accommodation
effects during
construction stage
pollution
archaeological assets
sociocultural effects
on life style and
integration into the
community
landscape
Fig. 9.1 - Linkage between natural and socio-economic impacts in an

EIA study for a small hydropower plant (SHP).
- 156 -
Environment
9.2- Identification of impacts

Water quality - Water in natural rivers is never completely pure, varying,
considerably, in the range and in concentrations of dissolved substances present,
such as suspended particles, pH and temperature. Fast-flowing streams and
rivers, normally, have high oxygen levels absorbed from the atmosphere due to
natural turbulence. Still waters, such as reservoirs or forebays have highly
variable levels that may range from a minimum during groups operation to a
maximum during non-operation period. Stratification can occur that in upper
layer is well oxygenated and in the lower layer there is no any contact with the
atmosphere suffering oxygen depletion.
The aquatic life is adapted to temperature regimes presented in freshwater
systems. Changes in temperature can eventually promote oxygen consumption.
Nevertheless, in these type of systems, there is no significant change in water
temperature.
Sedimentation Dams will retain solid material transported by the river. The
reduction of solid flow for the downstream river reach may cause concern due to
the increased potential erosion capacity. Small reservoirs and sediment discharge
outlets provided in the dam will mitigate this impact..
Air - There are no emissions of toxic gaseous, such as oxides of sulphur and
nitrogen that could affect the natural pH and, consequently, the water quality.
Although, during the construction is necessary the installation of a concrete
central and due to excavations and the circulation of heavy vehicles these can
induce particle emissions and gases (e.g. CO, CO2 and SOx).
Pollution During exploitation, in powerhouses the oil of generator lubrication
and transformers must be caught in appropriate site, in order to avoid the water
or land contamination with complex toxicity. Residual materials from the
construction must be removed and located in appropriated sites.
Gas emissions by vehicles and machines in working and siltation by excavations
are a matter of concern.
Soils The implementation of the hydro scheme should take into account the
soil ability, in particular for agricultural purposes. During the construction phase,
soils occupation by civil works and reservoir (in case of existing) can affect the
- 157 -
soils in a non-reversible way. Excavation, land moving and superficial

deforestation are relevant actions in soil impacts.
Noise - The noise produced during hydropower operation caused by the
turbine(s), the generator(s) and the cooling ventilation, is transmitted to the
powerhouse structure and can be harmful for the external or outdoor
environment. The external noise level can be minimised by improving the
acoustic isolation of the powerhouse and turbine, controlling vibrations of the
ventilation system, improving hydrodynamic design of hydraulic structures and
by using non-reflecting and sound absorbing materials such as fibreglass mat,
false ceiling and heavy mass trap doors. The sound produced in a powerhouse
through the generator is about 82 dB(A) at 7 m of distance. A covered wall can
decrease it about 25 50 dB(A) using specific soundproofing materials. The
World Health Organisation (WHO) suggests that daytime outdoor noise levels
should be below 50 dB(A) (i.e. Leq levels must be lesser than 55 dB(A) in
daytime and 45 dB(A) at night, where Leq is the average notorious steady noise
sounds for a given period of time). During the construction, a temporary
increasing of sound levels caused by working machines will be verified. The
noise level study should include also the flowing water, especially at the tailrace
(as it can be higher than 60 70 dB(A)).
Biota Civil works will provoke several disturbances in fish habitats. Any
turbined flow variation due to power generation change can affect these habitats.
It is required that any engineering works provide a minimum flow (or reserved
flow), in order to ensure the life and the free passage of migrating fishes. This
happens more frequently whenever there are diversion works with reduction of
flow between the intake and the powerhouse. The reserved flow depends on the
guarantee of aquatic biota protection and the physical characteristics of the
catchment of the hydropower schemes. Nevertheless, as an approximate value,
this flow can be about 5 to 10% of the long-term average flow. Fish-passages or
fish ladders should be installed to allow migratory fish. A fish-ladder should
allow the movement of the existing species in the ecosystem by providing a well
designed access and flow rates whether by swimming or jumping.
- 158 -
Environment
Fig. 9.2 A fish-ladder of a SHP.
Landscape Civil works and exploration of stone quarries, lend areas, the
transport of materials and sites for waste deposition should not affect areas with
high landscape quality. Deforestation, land removal and the installation of long
conveyance systems (e.g. canal or penstock) may have a harmful visual impact
in higher slope zones or mountainsides. The constructions of the dam and the
powerhouse, as well as the access imply local deforestation, creating a clearing
in the woods with potential significant morphologic alterations.
Dam
Natural falls
Fig. 9.3 A good example of SHP dam in landscape integration.
- 159 -
Fig. 9.4 An example of a SHP powerhouse in landscape integration.
Fig. 9.5 Deforestation for pipe installation.
Social-economics, land arrangement and municipal planning Whenever

plans land arrangement policy and life population standards are affected or
modified special concerns get higher. During the construction will happen an
increase of the local traffic that can disturb the quotidian life of a rural isolated
- 160 -
Environment
region. The construction of a small hydropower plant or mini-hydro can affect

zones that can belong to ecological protected areas (EPA) or agricultural
protected areas (APA). Disturbances related with natural flow river interruption,
dragging solid materials to downstream and excavations are predictable in EPA.
Reservoir filling and vehicles motion can affect APA.
Geomorphology - Excavation and land depositions normally associated to the
dam and powerhouse construction and conveyance system installation can affect
the ground and geologic characteristics. The retention of water and solid
particles by a dam will modify the geomorphologic equilibrium of a river
inducing the erosion phenomenon.
Climate local change During the construction stage a small alteration of the
micro-clime conditions close to the soil can occur, through an increasing of the
temperature and wind, and a decreasing of relative air humidity due to
destruction of natural vegetation in the dam site, powerhouse access and
hydraulic conveyance installation. During exploitation it will occur more
frequently smog at the reservoir zone (if it is large enough).
Water resources The deforestation occurred during the construction induces
the production of organic material that induces the degradation of the water
resources. Land removing also becomes water turbidity.
Archaeological and historic heritage It is advisable to develop surveys in situ
and to have necessary caution for special local heritage materials, namely old
mills and barrages or aqueducts or cultural heritage monuments.
Developers and engineers must be aware of what it is required to be taken into
consideration during the different stages of a project (e.g. during the planning,
construction and exploitation).
9.3- Activities associated with hydropower plants

The use of an area for construction will affect the soil structure, reducing the
porosity and the aeration, diminishing the diversity of the soil biota, and
damaging vegetation through the destruction, reducing plant growth and ability
land to recover. For a low dam and powerhouse these effects are not especially
relevants.
- 161 -
Table 9.1- Activities and impacts due to a small hydropower plant project
Activities
- River engineering
- Pipe installations
- Channelization
- Dam and
Powerhouse
construction
- Reservoir
(even small)
- Access - roads
- Powerhouse
- Excavations
- Deforestation
Potential Impacts
Enhanced erosion and silt production
Destruction of some vegetation
Damage in habitats and biota
Land expropriation and access affectation
Wastes creation and noise production
Different land uses affected
Localised changes in the water level
Loss of terrestrial habitats (at small scale)
Changes in flow regimes, detention time, siltation, depth
and potential stratification
Changes in groundwater recharge areas and in the water
quality
Changes in drainage systems due to landscaping (e.g.
gradient changes and embankments)
Noise production (noise control will be possible)
Increase runoff, affect stream flows and stream sediment
loads and modify soil compaction
Reduced interception with consequent channel erosion,
flood risk and runoff (can be minimised)
A small hydro scheme that comprises different components (i.e. a low dam,
weir, spillway, intake, channel/penstock, powerhouse, tailrace, substation and
transmission electric lines) changes the visual aspect of the valley. Nevertheless,
the most significant impacts occur during the construction phase. Some
architectural creativity for the powerhouse external aspect and integration can
avoid contrasting with the background and landscape colours resulting in an
acceptable visual project (e.g. colours of penstock and powerhouse).
Buried pipes in order to avoid a barrier for wildlife or people passage or
agricultural works, utilising local materials (e.g. stones or rocks to adornment
channel walls and powerhouses or small dams and weirs are appropriated),
reforesting excavation areas, according to local existent vegetation, are examples
to improve the compatibility of artificial works with the environment.
- 162 -
Environment
9.4- Reference criteria to support environmental studies

University of Manchester Centre developed a methodology based on a list of
revision criteria hierarchically organised that includes the following stages: 1)
project description and reference situation; 2) identification and assessment of
main impacts; 3) alternative solutions and mitigation; 4) communication of final
results (see Figure 9.6).
Reference situation
and project
description
Identification and
evaluation of key
impacts
Alternative solutions
Reference situation
Project description and

analyses of solutions
Description of the place to

install the project
Quantification of residues
production
Description of envelop
environment
Alternatives and
mitigation
Definition of impacts:
positive, reversible,
permanent and direct
Identification of environment
factors: water and air quality,
land use, noise, landscape,
public health, employment,
economy, road net,
architectonic heritage
Scoping and
effectiveness of
mitigation measures
Devotion on mitigation
Scoping definition
Foresight the magnitude and

significance
Comunication of results:
format, presentation, emphasis, no
technical resume
Fig. 9.6 - Environmental impact assessment criteria.
Mitigation strategies are based on elimination, reduction or compensation

measurements. The development of an EIA should attend the existing
- 163 -
legislation, the decision of the project proponent or by imposition of responsible

organisms or even for licence purposes.
The content of an EIA should characterise the environmental system (scooping
analysis and definition), the relation between project components and
environment changes (screening), identification of affected environmental
factors, impacts quantification through specific models, environmental effect
assessment and consequences and, finally, a correct and sufficiently clear
presentation of results. This analysis summarises three important components: a)
identification of the type and the cause; b) foresight the quantity resulting from
relationship between cause-effect of each impact; c) consequence assessment.
9.5- Safety
The exploitation of a small hydropower scheme can be dangerous. In fact, the
water volume stored by the dam and the flow towards the intake can be danger
for children and others that fall into the water or want to enjoy the benefits of an
artificial lake. Protection devices should be installed to safe anyone and to
prevent the suction towards the intake. Downstream the powerhouse there is
always the danger associated to a sudden flow variation during a turbine
operation. This hazard should be clearly noticed and sound signals should be
considered, in order to warn fishermen and other people crossing the river.
Low dams do not constitute a great risk for downstream valley, regarding a
failure scenario due to an abnormal hazard. However, the potential damages and
losses along the valley due to a dam-break flood need to be evaluated.
Emergency plans and risk zoning should be considered in certain special cases.
- 164 -
Environment
Fig. 9.7 A life guard grid with float balls.
- 165 -
- 166 -
Economic analysis
ECONOMIC ANALYSIS
M. Manuela Portela
10.1-Introduction. Costs and benefits

The final decision on whether or not a small hydropower scheme should be
constructed, or the selection among alternative design solutions for the same is
generally based on the comparison of the expected costs and benefits for the
useful life of the project, by means of economic analysis criteria. This analysis
should be performed in the first stages of the design (along with the feasibility
study) as nothing ensures that a project suitable from a technical point of view is
also advantageous from an economic point of view.
It should be pointed out that the choice among alternative solutions having
identical benefits is simply accomplished by identifying which one has the lower
cost or the less negative environmental impact. Only the comparison of
alternative projects or design solutions with different costs and benefits requires
the application of economic criteria in order to identify the most desirable
alternative.
The effectiveness of the economic analysis as a decision tool of a small hydro
investor depends on the accuracy of the project cost and benefits estimates.
These estimates are not easy to reach, especially in the early stages of design
where some of the scheme characteristics are only preliminary defined.
- 167 -
The costs of a small hydropower scheme can be grouped in the three following
categories, also systematised in Table 10.1:
The capital cost that may be defined as the sum of all expenditures required
to bring the project to completion. These costs occur during the construction
period (generally, from one to three years).
The annual operation costs resulting from the exploitation and maintenance
of the scheme during its useful live.
The reposition costs concerning to the substitution of the equipment having
a useful live lesser than the one of the scheme.
Table 10.1 Costs of a small hydropower scheme
Capital costs
Annual operation costs
Studies and design

Supervision during the execution
Civil works
Equipment
Land acquisition
Contingencies or unforeseen cost
Exploitation
Maintenance
Spare parts
Grant of permission
Reposition costs
The following items should generally be considered within the capital costs
Table 10.1:
Studies and design.

Supervision during the execution.
Civil works (including access roads and building site facilities).
Equipment (including hydromechanical, electromechanical and electrical
equipment and the interconnection to the electrical grid).
Land acquisition.
It should also be considered an item for contingencies or unforeseen cost. This

item represents expenditures that are possible but not certain or yet foreseen. For
instance, this item intends to overcome the uncertainty that result from the fact
- 168 -
Economic analysis
that the characteristics of the site where the scheme is going to be constructed
are not totally known, especially in the first design phases and, so, they may be
incompletely expressed in the scheme conception. The uncertainty under
consideration is mainly related with the civil works and can result, for instance,
from undervalued excavation or landing works or from bedrock having
characteristics worse than those foreseen in the design.
The studies and design costs and the supervision costs during the scheme
execution result from agreements between the investor and consulting firms. If
the investor has previous experience in the execution of small hydropower
schemes, he should be able to express the costs under consideration as a
percentage of the civil works and equipment costs.
The civil work costs are evaluated from the design, by measuring the work
quantities relative to the different components of the scheme and affecting them
of unit prices, which generally are not difficult to obtain from each country civil
work contractors market. To evaluate the equipment costs, budget prices from
the suppliers should also be obtained.
The cost of the acquisition (or renting) of the land that will be occupied by the
scheme (including accesses and the area that will be submerged by the scheme
reservoir) depends strictly on the land valorisation in each country or region.
If possible, the capital costs estimated for the small hydropower scheme under
study should be compared with the ones of similar schemes, already built or
previously characterised. This procedure is more important in what concerns the
evaluation of the equipment cost as frequently it is not possible to get the budget
prices from the suppliers in a period as short as the design one.
The annual operation costs include the following main components Table 10.1:
Exploitation.
Maintenance.
Spare parts.
Grant of permission or legal permit.
The exploitation costs represent the charges with the staff responsible for the
scheme operation. To reduce these costs the scheme should be made fully
- 169 -
automatic, that is to say, explored in an abandon mode. This kind of exploitation

requires additional equipment (for instance, automation systems, telephone or
satellite lines for telemetry and alarm signals).
The maintenance costs include two parcels, related one with the civil works and
the other, with the equipment. The former is generally the smallest representing
from 0,25% to 0,50% of the capital costs in civil woks, while the latter can reach
about 2,50% of the respective capital costs.
The spare parts costs are the costs with the reposition of the material that is
necessary to keep in stock in order to perform the maintenance of the
hydropower scheme or to execute small repairs in the same. These costs can be
assigned to annual costs or to costs without periodicity, occurring whenever it is
necessary to restore the stock of the material.
The grant of permission costs occurs once the scheme starts to operate and
represent the annual payments due for the scheme license and for the water
utilisation. Although these costs result from the legislation in force in each
country, in terms of economic appraisal they can be treated as fixed percentages
of the energy incomes.
From the investor point of view, the only tangible revenue or benefit in a small
hydropower scheme is the annual income with the energy production sale. This
income depends on the amount of energy produced during the scheme lifetime
and on the specific conditions that rule the hydroelectric sector, namely the
energy sale contract conditions and the tariffs policy, which are specific in each
country.
In the further development of this chapter it will be assumed that the purchase of
the energy produced in a small hydropower scheme is ensured, no matter the
amount and the characteristics of the production. This scenario should be real as
it translates the special rule that the renewable energies are expected to play in
the energetic policy, not only of Europe, but also of the world.
The evaluation of the energy production requires accurate hydrological studies
based on specific methodologies that overcome the non-existence of basic
hydrological data that characterises most of the small hydropower schemes,
generally located in small and ungauged watersheds. The results of the
- 170 -
Economic analysis
hydrological study have to be interpreted taking into account the characteristics

of the watershed that condition the water availability and regime at the water
intake of the scheme: watershed geology, vegetation covering, land use, water
consumption, existing upstream storage reservoirs, ecological requirements.
As the energy production that will result from the future inflows at the scheme
water intake can not be specified, the revenues in the economic analysis are
considered to be fixed during the lifetime of the project and based on a annual
long-term average discharge of the river. By other words, the revenues are
generally defined on the basis of an annual production constant and equal to the
expected mean one, E .
The uncertainty concerning the future inflows and, thus, concerning the future
energy productions, can be regard as a hydrologic risk. Based on the statistic
analysis of the historic flow series, this risk can be expressed, for instances, in
terms of the probability of occurring future periods of dry years or of years with
annual flows bellow a limit. This kind of analysis is most important and should
be performed in countries, as those of south Europe, characterised by an
irregular flow regime where dry periods of several years can occur.
Based on simulation studies, PORTELA and ALMEIDA, 1995, carried out
hydrologic risk analysis for several hypothetical small hydropower schemes
located in the North of Portugal. This analysis clearly shows that the profitability
of those schemes could be drastically reduced if periods of dry years occur,
especially in the first years of exploitation, when the incomes are crucial for the
capital recovery.
Figure 10.1 represents some of the results achieved by those authors. The cost
factor, , was defined as the ratio between the income under real flow conditions
and the income evaluated on the basis of an annual production constant and
equal to the mean annual production. Two exploitation periods where
considered: one with 33 years and the other with only 10 years. Based on real
flow data, simulation studies were performed in order to evaluate the maximum
(best series) and the minimum (worst series) incomes during each of the
previous periods. The results thus achieved were represented as a function of the
ratio between the design flow, Qmax, and the average mean daily flow, Qave.
- 171 -

Cost factor,
1.3
Analysis period of 10 years
1.1
Analysis period of 33 years
0.9
0.7
0.5
0.0
0.5
1.0
1.5
2.0
2.5
Best series
3.0
3.5
4.0
4.5
5.0
Qmax/Qave
Worst series
Figure 10.1 Variation of the income in a hypothetical small hydropower scheme,

located in the North of Portugal, due to the natural variability of the flow regime
(PORTELA and ALMEIDA, 1995).
10.2-Economic analysis
10.2.1-Introduction
In this item some concepts necessary to the economic analysis or appraisal of
small hydropower schemes are presented in a simple and straight way.
Several literature references are available either in general terms of economic
and financial analysis principles and applications (which are far beyond the
scope of this Guide) or in specific terms for small hydropower schemes. Among
these last ones KUIPER, 1981, ESHA, 1994 and 1997, and JIANDONG et al.,
1996, may also have interest in the domain of the present chapter.
10.2.2- Constant market prices system concept

Although the design of a small hydropower scheme generally foresees a
schedule for the investment costs and revenues, the respective estimates are
based on market prices referred to a given year, the one during which the studies
were performed. The evolution of the inflation during the project lifetime and its
effect on each component of the previous cost estimates is practically impossible
to establish.
- 172 -
Economic analysis
So, one could say that one of the main problems of the economic analysis in to
conceive a scenario for the future evolution of the inflation. To overcome this
problem a common and simple economic approach based on a constant market
prices system referred to a given year is generally applied in the comparison of
costs and benefits either of a project or of alternative design solutions for the
same. This approach, which will be adopted in the present chapter, assumes that
it is not necessary to account for the inflation, as it will have the same effect in
any monetary flux. The future costs and benefits are, then, evaluated at present
market prices.
10.2.3- Discount rate and present value concept

A given monetary unit (for instances, a pound or a dollar) is more worthy in the
present than in the future. By other words, the future value of the monetary flux
that has a present value of one unit will be great than one. At the same time, the
present value of a future unitary monetary flux will be lesser than one.
Through the years, this situation generates different appetencies to transfer
money from the present to the future and vice-verse. Theses appetencies can
be expressed in terms of different discount rates, r. The values of these rates
depend, among other factors, on the state of the economy, on the risk that
involves the investment, the capital availability and on the expected future rate
of inflation.
Let n denote a period of n years, from year 1 to year n. According to the discount
rate concept and as represented in Figure 10.2, one monetary unit of today will
be changed in year n by (1+r)n monetary units and one monetary unit of year n
will be changed today by 1/(1+r)n units. If r is defined, not on a yearly basis, but
for a generic period of time , n will denote the total number of periods
between the two monetary fluxes of Figure 10.2.
- 173 -
(1+ r)n
n-2
n-1
Year
(1 + r) n
1
n-2
n-1
Year
Fig. 10.2 Transference of unitary monetary fluxes by means of the

discount rate, r.
The factor
(1 + r )n
that expresses the depreciation suffered by future monetary
fluxes when transferred to the present is call present worth or value factor.
According to this concept, the present value, PV, of a single generic monetary
flux occurring in future year i, Ci, is given by
PV =
(1 + r )i
Ci
(10.1)
The present value for the beginning of year 1 of the continuous sequence of
annual monetary fluxes represented in Figure 10.3 is given by
PV =
1
1
1
C1 +
C2 + +
Ci + +
2
(1 + r )
(1 + r )
(1 + r )i
1
1
++
C n 1 +
Cn
(1 + r )(n 1)
(1 + r )n
n
PV =
(1 + r ) i C i
1
(10.2)
(10.3)
i =1
- 174 -
Economic analysis
PV
C1
C2
C3
C4
Ci
Annual monetary flux
Cn-2
n-2
Cn-1
n-1
Cn
Present value
Fig. 10.3 Continuos sequence of monetary fluxes.
For economical analysis purpose the monetary fluxes are grouped in periods
(generally years) and are considered to occur as represented in Figure 10.3, that
is to say, concentrated in the end of each of those periods. The present value
operations are performed for the beginning of the year adopted as reference
(year one in the previous figure).
If the monetary fluxes of Figure 10.3 are constant and equal to an annuity, C,
uniform series of annual monetary fluxes the following relation is achieved
for the respective present value referred to the beginning of year 1:
n
PV = C
(1 + r ) i
(1 + r )n 1
(1 + r )n r
is the present worth factor for uniform series. The inverse of
i =1
where
(1 + r )n 1
(1 + r )n r
=C
(10.4)
this factor is denoted by capital-recovery factor.

The present worth factor either for a single monetary flux or for a series of
uniform monetary fluxes is generally provided by tables as the ones presented in
next pages, which can be found in most books dealing with economic analysis
concepts.
- 175 -
Year
The factor (1+ r )n represented in the first part of Figure 10.2 is the
capitalisation factor or the future value factor. Table 10.2 also includes its value
either for a single monetary flux or for an uniform series of monetary fluxes. The
formula that provides the future value factor for an uniform series of monetary
(1 + r )n 1
fluxes is given by
. The inverse of the previous ratio is the
r
sinking-fund factor.
Although the computation of the previous factors does not offer special
difficulty it is useful to provide tables as those included herein in order to allow
expedite evaluations or to control results of the economic appraisal.
It should be pointed out that the expressions and concepts presented so far, as
well as their development and application in the rest of this chapter, assume that
the discount rate, r, is constant during the n years period under analysis. If this is
not the situation, different sub periods having different discount rates should be
considered. However, this is not a common procedure in the economic analysis
of small hydropower schemes.
- 176 -
Economic analysis
Table 10.2 (1/4) Values of the present worth factor and of the capitalisation factor for
different time periods, n, and discount rates, r
1
a1) Present worth factor for single monetary fluxes
n
(1 + r )
Time
Period, n
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
35
40
45
50
60
70
80
90
100
Discount rate, r, of:

8%
10%
2%
4%
6%
12%
14%
0,98039
0,96117
0,96154
0,92456
0,94340
0,89000
0,92593
0,85734
0,90909
0,82645
0,89286
0,79719
0,87719
0,76947
0,94232
0,88900
0,83962
0,79383
0,75131
0,71178
0,67497
0,92385
0,85480
0,79209
0,73503
0,68301
0,63552
0,59208
0,90573
0,82193
0,74726
0,68058
0,62092
0,56743
0,51937
0,88797
0,79031
0,70496
0,63017
0,56447
0,50663
0,45559
0,87056
0,75992
0,66506
0,58349
0,51316
0,45235
0,39964
0,85349
0,73069
0,62741
0,54027
0,46651
0,40388
0,35056
0,83676
0,70259
0,59190
0,50025
0,42410
0,36061
0,30751
0,82035
0,67556
0,55839
0,46319
0,38554
0,32197
0,26974
0,80426
0,64958
0,52679
0,42888
0,35049
0,28748
0,23662
0,78849
0,62460
0,49697
0,39711
0,31863
0,25668
0,20756
0,77303
0,60057
0,46884
0,36770
0,28966
0,22917
0,18207
0,75788
0,57748
0,44230
0,34046
0,26333
0,20462
0,15971
0,74301
0,55526
0,41727
0,31524
0,23939
0,18270
0,14010
0,72845
0,53391
0,39365
0,29189
0,21763
0,16312
0,12289
0,71416
0,51337
0,37136
0,27027
0,19784
0,14564
0,10780
0,70016
0,49363
0,35034
0,25025
0,17986
0,13004
0,09456
0,68643
0,47464
0,33051
0,23171
0,16351
0,11611
0,08295
0,67297
0,45639
0,31180
0,21455
0,14864
0,10367
0,07276
0,65978
0,43883
0,29416
0,19866
0,13513
0,09256
0,06383
0,64684
0,42196
0,27751
0,18394
0,12285
0,08264
0,05599
0,63416
0,40573
0,26180
0,17032
0,11168
0,07379
0,04911
0,62172
0,39012
0,24698
0,15770
0,10153
0,06588
0,04308
0,60953
0,37512
0,23300
0,14602
0,09230
0,05882
0,03779
0,59758
0,36069
0,21981
0,13520
0,08391
0,05252
0,03315
0,58586
0,34682
0,20737
0,12519
0,07628
0,04689
0,02908
0,57437
0,33348
0,19563
0,11591
0,06934
0,04187
0,02551
0,56311
0,32065
0,18456
0,10733
0,06304
0,03738
0,02237
0,55207
0,30832
0,17411
0,09938
0,05731
0,03338
0,01963
0,50003
0,25342
0,13011
0,06763
0,03558
0,01894
0,01019
0,45289
0,20829
0,09722
0,04603
0,02209
0,01075
0,00529
0,41020
0,17120
0,07265
0,03133
0,01372
0,00610
0,00275
0,37153
0,14071
0,05429
0,02132
0,00852
0,00346
0,00143
0,30478
0,09506
0,03031
0,00988
0,00328
0,00111
0,00039
0,25003
0,06422
0,01693
0,00457
0,00127
0,00036
0,00010
0,20511
0,04338
0,00945
0,00212
0,00049
0,00012
0,00003
0,16826
0,02931
0,00528
0,00098
0,00019
0,00004
0,00001
0,13803
0,01980
0,00295
0,00045
0,00007
0,00001
0,00000
- 177 -
(1 + r )n 1
a2) Present worth factor for uniform monetary fluxes
n
(1 + r ) r
Time
period, n
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
35
40
45
50
60
70
80
90
100
2%
Uniform series of monetary fluxes and discount rate, r, of:

4%
6%
8%
10%
12%
14%
0,98039
1,94156
0,96154
1,88609
0,94340
1,83339
0,92593
1,78326
0,90909
1,73554
0,89286
1,69005
0,87719
1,64666
2,88388
2,77509
2,67301
2,57710
2,48685
2,40183
2,32163
3,80773
3,62990
3,46511
3,31213
3,16987
3,03735
2,91371
4,71346
4,45182
4,21236
3,99271
3,79079
3,60478
3,43308
5,60143
5,24214
4,91732
4,62288
4,35526
4,11141
3,88867
6,47199
6,00205
5,58238
5,20637
4,86842
4,56376
4,28830
7,32548
6,73274
6,20979
5,74664
5,33493
4,96764
4,63886
8,16224
7,43533
6,80169
6,24689
5,75902
5,32825
4,94637
8,98259
8,11090
7,36009
6,71008
6,14457
5,65022
5,21612
9,78685
8,76048
7,88687
7,13896
6,49506
5,93770
5,45273
10,57534
9,38507
8,38384
7,53608
6,81369
6,19437
5,66029
11,34837
9,98565
8,85268
7,90378
7,10336
6,42355
5,84236
12,10625
10,56312
9,29498
8,24424
7,36669
6,62817
6,00207
12,84926
11,11839
9,71225
8,55948
7,60608
6,81086
6,14217
13,57771
11,65230
10,10590
8,85137
7,82371
6,97399
6,26506
14,29187
12,16567
10,47726
9,12164
8,02155
7,11963
6,37286
14,99203
12,65930
10,82760
9,37189
8,20141
7,24967
6,46742
15,67846
13,13394
11,15812
9,60360
8,36492
7,36578
6,55037
16,35143
13,59033
11,46992
9,81815
8,51356
7,46944
6,62313
17,01121
14,02916
11,76408
10,01680
8,64869
7,56200
6,68696
17,65805
14,45112
12,04158
10,20074
8,77154
7,64465
6,74294
18,29220
14,85684
12,30338
10,37106
8,88322
7,71843
6,79206
18,91393
15,24696
12,55036
10,52876
8,98474
7,78432
6,83514
19,52346
15,62208
12,78336
10,67478
9,07704
7,84314
6,87293
20,12104
15,98277
13,00317
10,80998
9,16095
7,89566
6,90608
20,70690
16,32959
13,21053
10,93516
9,23722
7,94255
6,93515
21,28127
16,66306
13,40616
11,05108
9,30657
7,98442
6,96066
21,84438
16,98371
13,59072
11,15841
9,36961
8,02181
6,98304
22,39646
17,29203
13,76483
11,25778
9,42691
8,05518
7,00266
24,99862
18,66461
14,49825
11,65457
9,64416
8,17550
7,07005
27,35548
19,79277
15,04630
11,92461
9,77905
8,24378
7,10504
29,49016
20,72004
15,45583
12,10840
9,86281
8,28252
7,12322
31,42361
21,48218
15,76186
12,23348
9,91481
8,30450
7,13266
31,72839
21,57725
15,79217
12,24336
9,91810
8,30561
7,13304
31,97842
21,64146
15,80910
12,24793
9,91936
8,30597
7,13315
32,18353
21,68485
15,81855
12,25005
9,91985
8,30609
7,13317
32,35179
21,71416
15,82383
12,25104
9,92004
8,30612
7,13318
32,48982
21,73396
15,82678
12,25149
9,92011
8,30614
7,13318
- 178 -
Economic analysis
b1) Capitalisation factor for single monetary fluxes (1+ r ) n

Time
Period, n
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
35
40
45
50
60
70
80
90
100
Discount rate, r, of:

6%
8%
10%
2%
4%
12%
14%
1,02000
1,04040
1,04000
1,08160
1,06000
1,12360
1,08000
1,16640
1,10000
1,21000
1,12000
1,25440
1,14000
1,29960
1,06121
1,12486
1,19102
1,25971
1,33100
1,40493
1,48154
1,08243
1,16986
1,26248
1,36049
1,46410
1,57352
1,68896
1,10408
1,21665
1,33823
1,46933
1,61051
1,76234
1,92541
1,12616
1,26532
1,41852
1,58687
1,77156
1,97382
2,19497
1,14869
1,31593
1,50363
1,71382
1,94872
2,21068
2,50227
1,17166
1,36857
1,59385
1,85093
2,14359
2,47596
2,85259
1,19509
1,42331
1,68948
1,99900
2,35795
2,77308
3,25195
1,21899
1,48024
1,79085
2,15892
2,59374
3,10585
3,70722
1,24337
1,53945
1,89830
2,33164
2,85312
3,47855
4,22623
1,26824
1,60103
2,01220
2,51817
3,13843
3,89598
4,81790
1,29361
1,66507
2,13293
2,71962
3,45227
4,36349
5,49241
1,31948
1,73168
2,26090
2,93719
3,79750
4,88711
6,26135
1,34587
1,80094
2,39656
3,17217
4,17725
5,47357
7,13794
1,37279
1,87298
2,54035
3,42594
4,59497
6,13039
8,13725
1,40024
1,94790
2,69277
3,70002
5,05447
6,86604
9,27646
1,42825
2,02582
2,85434
3,99602
5,55992
7,68997
10,57517
1,45681
2,10685
3,02560
4,31570
6,11591
8,61276
12,05569
1,48595
2,19112
3,20714
4,66096
6,72750
9,64629
13,74349
1,51567
2,27877
3,39956
5,03383
7,40025
10,80385
15,66758
1,54598
2,36992
3,60354
5,43654
8,14027
12,10031
17,86104
1,57690
2,46472
3,81975
5,87146
8,95430
13,55235
20,36158
1,60844
2,56330
4,04893
6,34118
9,84973
15,17863
23,21221
1,64061
2,66584
4,29187
6,84848
10,83471
17,00006
26,46192
1,67342
2,77247
4,54938
7,39635
11,91818
19,04007
30,16658
1,70689
2,88337
4,82235
7,98806
13,10999
21,32488
34,38991
1,74102
2,99870
5,11169
8,62711
14,42099
23,88387
39,20449
1,77584
3,11865
5,41839
9,31727
15,86309
26,74993
44,69312
1,81136
3,24340
5,74349
10,06266
17,44940
29,95992
50,95016
1,99989
3,94609
7,68609
14,78534
28,10244
52,79962
98,10018
2,20804
4,80102
10,28572
21,72452
45,25926
93,05097
188,88351
2,43785
5,84118
13,76461
31,92045
72,89048
163,98760
363,67907
2,69159
7,10668
18,42015
46,90161
117,39085 289,00219
700,23299
3,28103
10,51963
32,98769
101,25706
304,48164
3,99956
15,57162
59,07593
218,60641
789,74696 2787,79983 9623,64498
4,87544
23,04980
105,79599
471,95483 2048,40021 8658,48310 35676,9818
5,94313
34,11933
189,46451 1018,91509 5313,02261 26891,9342 132262,467
7,24465
50,50495
339,30208 2199,76126 13780,6123 83522,2657 490326,238
- 179 -
897,59693 2595,91866
(1 + r )n 1
b2) Capitalisation factor for uniform monetary fluxes
Time
period, n
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
35
40
45
50
60
70
80
90
100
2%
Uniform series of monetary fluxes and discount rate, r, of:

4%
6%
8%
10%
12%
14%
1,00000
2,02000
1,00000
2,04000
1,00000
2,06000
1,00000
2,08000
1,00000
2,10000
1,00000
2,12000
1,00000
2,14000
3,06040
3,12160
3,18360
3,24640
3,31000
3,37440
3,43960
4,12161
4,24646
4,37462
4,50611
4,64100
4,77933
4,92114
5,20404
5,41632
5,63709
5,86660
6,10510
6,35285
6,61010
6,30812
6,63298
6,97532
7,33593
7,71561
8,11519
8,53552
7,43428
7,89829
8,39384
8,92280
9,48717
10,08901
10,73049
8,58297
9,21423
9,89747
10,63663
11,43589
12,29969
13,23276
9,75463
10,58280
11,49132
12,48756
13,57948
14,77566
16,08535
10,94972
12,00611
13,18079
14,48656
15,93742
17,54874
19,33730
12,16872
13,48635
14,97164
16,64549
18,53117
20,65458
23,04452
13,41209
15,02581
16,86994
18,97713
21,38428
24,13313
27,27075
14,68033
16,62684
18,88214
21,49530
24,52271
28,02911
32,08865
15,97394
18,29191
21,01507
24,21492
27,97498
32,39260
37,58107
17,29342
20,02359
23,27597
27,15211
31,77248
37,27971
43,84241
18,63929
21,82453
25,67253
30,32428
35,94973
42,75328
50,98035
20,01207
23,69751
28,21288
33,75023
40,54470
48,88367
59,11760
21,41231
25,64541
30,90565
37,45024
45,59917
55,74971
68,39407
22,84056
27,67123
33,75999
41,44626
51,15909
63,43968
78,96923
24,29737
29,77808
36,78559
45,76196
57,27500
72,05244
91,02493
25,78332
31,96920
39,99273
50,42292
64,00250
81,69874
104,76842
27,29898
34,24797
43,39229
55,45676
71,40275
92,50258
120,43600
28,84496
36,61789
46,99583
60,89330
79,54302
104,60289
138,29704
30,42186
39,08260
50,81558
66,76476
88,49733
118,15524
158,65862
32,03030
41,64591
54,86451
73,10594
98,34706
133,33387
181,87083
33,67091
44,31174
59,15638
79,95442
109,18177 150,33393
208,33274
35,34432
47,08421
63,70577
87,35077
121,09994 169,37401
238,49933
37,05121
49,96758
68,52811
95,33883
134,20994 190,69889
272,88923
38,79223
52,96629
73,63980
103,96594
148,63093
214,58275
312,09373
40,56808
56,08494
79,05819
113,28321
164,49402
241,33268
356,78685
49,99448
73,65222
111,43478
172,31680
271,02437
431,66350
693,57270
60,40198
95,02552
154,76197
259,05652
442,59256
767,09142 1342,02510
71,89271
121,02939
212,74351
386,50562
718,90484 1358,23003 2590,56480
84,57940
152,66708
290,33590
573,77016 1163,90853 2400,01825 4994,52135
114,05154
237,99069
533,12818 1253,21330 3034,81640 7471,64111 18535,1332
149,97791
364,29046
967,93217 2720,08007 7887,46957 23223,3319 68733,1784
193,77196
551,24498 1746,59989 5886,93543 20474,0021 72145,6925 254828,441
247,15666
827,98333 3141,07519 12723,9386 53120,2261 224091,118 944724,767
312,23231 1237,62370 5638,36806 27484,5157 137796,123 696010,547 3502323,12
- 180 -
Economic analysis
10.2.4- Economic indexes

10.2.4.1-Basic considerations
The evaluation of the profitability of a small hydropower scheme or the
comparison of alternative design solutions for the same can be based on
economic indexes or parameters. Among these parameters the next five ones
will be treated herein: net present value (NPV), benefit/cost ratio (B/C), internal
rate of return (IRR), average price of the kWh (AP) and payback period (T).
The economic parameters do not always identify the same project as being the
most economic amongst a number of alternatives. At the same time, the
preference of one particular parameter may alter with a change, even small, in
the discount rate.
Being the results of the economic analysis so sensitive to the value considered
for the discount rate much attention should be paid in the establishment of this
value. However, in most situations, there is not a well-defined and unique
discount rate but, instead, a range of possible discount rates. Although this range
is often narrow, it is advisable to perform a sensitivity analysis of the project
response to different discount rates. By this way it will be possible to partially
overcome the uncertainty related with this rate and, at the same time, to conclude
if there is any change in the relative preferences indicated by the economic
parameters.
Specially to compare alternative design solutions for a small hydropower it is
usual to evaluate, for each of these solutions, the ratio of the total capital cost to
the installed capacity, that is to say, the cost per installed kW. Although this
index is quite easy to obtain it is poorly informative from an economic point of
view, as it does not take into account neither other costs of the project, nor its
benefits. However it may be useful in the systematisation of the capital costs
required for small hydropower projects or in comparison of alternative projects
having similar benefits.
In the next items the following symbols will be utilised:
n - number of the project lifetime periods (generally measured in years and

equal to the duration of the legal permits);
r - discount rate (constant along the project useful live);
- 181 -
Ci - capital costs in year i;

Oj - annual operation cost for year j;
Rj - revenues in year j;
Pm - reposition cost foreseen for year m ( n 2 < m < n );
E - mean annual production.
The abbreviations C, O, R and P will respectively denote the present values of

the costs, Ci, Oj, Pm and revenues Rj. If
the capital costs occur in the first k years,

the energy production takes place from year (k+1) until year n,
the equipment reposition is foreseen for year m,
the following relations are obtained for the beginning of the n years period:
k
C=
(1 + r )i
Ci
(10.5)
i =1
O=
j
j= k +1 (1 + r )
Oj
R=
P=
(10.6)
(1 + r )k
j
j= k +1 (1 + r )
Rj
(10.7)
(1 + r )k
Pm
(10.8)
(1 + r )m
The numerator of expressions (10.6) and (10.7) provides the present value (for
the beginning of year k+1) of an annuity that will occur during n-k years, from
year k+1 until year n. The denominator performs the transference of the previous
value from year k to the beginning of year 1 Figure 10.4.
- 182 -
Economic analysis
Expression (10.6) assumes that the parcel of the annual operation costs related
with the legal permits will occur only when the scheme starts to operate which is
generally a far enough realistic scenario.
.....
k+1
....
n-2
n-1
n Year
Present value operation applied to:

Present value
an uniform series of monetary fluxes;
Annuity
a single monetary flux.
Fig. 10.4 Compound present value operations.
As the annual benefits in a small hydropower scheme are generally evaluated on

a constant annual production basis (the mean annual production, E ), expression
(10.7) becomes
R
R=
(1 + r )(n k ) 1
(1 + r )(n k ) r
(1 + r )k
(10.9)
where R is the mean annual benefit obtained by multiplying the E by unitary

sale price of the energy in the first year of the period under analysis (year 1).
The period for which the economical analysis is performed should be equal to
the period foreseen in the legal permit of the hydropower exploitation. However,
it is not reasonable to adopt periods greater then twenty five to thirty years, as
the present values of the expected costs and benefits in the last years become so
small that their contribution to the economic indexes is negligible.
- 183 -
If the lifetime of the equipment is greater than the assumed useful life of the
hydropower scheme a residual value of the investment over its lifetime can be
considered and jointed to the benefits. As the estimate of this residual quantity is
quite subjective and fallible and its present value contribution for the scheme
profitability will certainly be very modest, it will not be considered in the further
development of this chapter.
10.2.4.2-Net Present Value (NPV)

The net present value, NPV, represents the cumulative sum of all expected
benefits during the lifetime of the project minus the sum of all its cost during the
same period, both expressed in terms of present values
NPV = R C O P
(10.10)
If NPV is negative, the project should be rejected as it is expected that the

benefits during its lifetime will be insufficient to cover the project costs.
Assuming that there are not any restrictions with respect to the required initial
capital availability, among projects or alternative design solutions with positive
NPV the best ones will be those with greater NPV.
The net present value can also be evaluated by computing the discount
cumulative cash flow. This cash flow gives, through each year of the economic
analysis period, the value of the cumulative sum of the present values of the
costs minus the present values of the benefits. The NPV equals the result
provided by the cash flow for the last year of the period under analysis (year n).
10.2.4.3-Benefit/cost ratio (B/C)

The benefit/cost ratio, B/C, compares the present values of the benefits and costs
of the hydropower scheme on a ratio basis. This index can be defined as the ratio
between present values of the net annual benefits and of the capital and
reposition costs or as the ratio between present values of the benefits and of the
total costs, that is to say
B C=
R O
C+P
(10.11)
- 184 -
Economic analysis
B C=
R
C+O+P
(10.12)
The first definition seems more coherent as it combines together the annual
monetary fluxes that occur during the project lifetime.
The B/C parameter has much popular appeal since it gives an immediate
indication of the degree of desirability of a project. If the benefit/cost is less
than one, the project is evidently undesirable. If it is exactly one the project has a
marginal interest and if it is greater than one, its implementation would seem
justified and as much as B/C is higher.
It should be pointed out that the B/C ratio evaluated by expression (10.12) is
always equal or greater then the B/C ratio that results from expression (10.11).
However, if (10.11) provides a unitary value, the B/C ratio given by (10.12) will
also be one. An unitary B/C ratio implies a NPV equal to zero.
10.2.4.4-Internal Rate of Return (IRR)

The internal rate of return, IRR, is defined as the discount rate that makes the net
present value, NPV, equal to zero.
Let us consider a project having capital costs, Ci, during the first k years and
energy revenues, Rj, and annual operation costs, Oj, from year j=(k+1) until year
n. A reposition cost is foreseen for year m. In this condition, the project NPV for
the beginning of year 1 can be evaluated according to
n
NPV = R C O P =
n
(1 + r ) j Rj
j= k +1
(1 + r )
i
i =1 (1 + r )
Ci
Oj
j
j= k +1 (1 + r )
(1 + r )
Pm
(1 + r )m
The respective IRR will be obtained by
- 185 -
(10.13)
NPV =
(1 + IRR ) j (Rj Oj)

j= k +1
(1 + IRR ) k
k
i =1 (1 + IRR )
Ci
Pm
(1 + IRR )m
=0
(10.14)
A discount rate equal to IRR will imply an unitary B/C ratio and a null NPV.
A process of trial and error can provide the solution of expression (10.14). This
process, which is generally simple to carry out with the software available
nowadays, begins with an arbitrary discount rate. If the NPV thus obtained is
positive, a higher discount rate must be tested; if NPV is negative, the search
must proceed with a lower discount rate Figure 10.5.
NPV
3000
2000
IRR=12,5%
1000
Dis count
0
5
10
15
20
25
r ate (%)
-1000
-2000
Figure 10.5 Internal rate of return, IRR.
The previous figure assumes that the project is well behaved from an economic
point of view. This is a correct assumption for small hydropower projects in
which the capital costs occur during a very short period of time (generally, less
than three years) while the incomes last during the useful life of the scheme (for
saying, more than twenty years).
Among projects or alternative design solutions with different internal rates of
return the best ones will be those with greater IRR. If the rates thus achieved are
- 186 -
Economic analysis
greater than the opportunity costs of the capital, the projects are advantageous
from an economic point of view.
10.2.4.5-Average price of the kWh (AP)

Although this parameter is not generally referred in expert literature, it is an
interesting index as it provides economic information in a way very simple to
understand. It is also useful in the comparison of alternative design solutions.
According to the abbreviations established, the average price of the kWh, AP, is
be defined by
n
(1 + r ) i Ci +
AP = i =1
Oj
j
j= k +1 (1 + r )
(1 + r )
(1 + r )(n k ) 1 E
(1 + r )(n k ) r
=
C+O+P
(1 + r )(n k ) 1 E
(1 + r )( n k ) r
Pm
(1 + r )m
(10.15)
This parameter represents the ratio between the present value of all costs and an
equivalent energy present amount. This last amount is obtained by applying
the present worth factor to the mean annual production. By other words, the AP
index represents the unitary energy sale price that makes the NPV equal to zero
and the B/C ratio equal to one.
If the AP price is lower than the nominal unitary energy-selling price the project
is justified from an economic point of view. The best project will be the one with
lower AP.
10.2.4.6-Payback period (T)

The payback period or recovery period, T, can be defined either on a discounted
or without discount basis and represents the number of years it takes before
- 187 -
cumulative forecasted cash flows equal the initial investment. Its value is
provided by the year when the cumulative cash flow changes from a negative
value to a positive value.
10.2.5-Sensitivity analysis
A small hydropower scheme project is usually characterised by a lack of
certainty about the capital cost estimates, future annual costs and future value of
the energy. So a sensitivity analysis should be performed in order to analyse the
project response capability to different scenarios. These scenarios can be
obtained by increasing the costs and maintaining or decreasing the benefits or, if
the cost estimates are considered to be accurate enough, maintaining the costs
and decreasing the benefits. Variations from 10% to 20% are usually adopted
in the previous analysis.
In what concerns the discount rate, at least two values should be considered
around the most probable discount rate: a pessimist value higher than the
expected one and a optimist value, lower than the expected discount rate.
10.2.6-Application example
Table 10.3 provides an application example for a hypothetical hydropower
scheme having an installed capacity of 2,2 MW and an average energy
production of 8,5 GWh.
The costs were considered to occur during the first two years (execution period)
and the benefits during the next twenty-eight years (total duration of the legal
licence period of 30 years). No residual costs, reposition costs and spare parts
costs were considered. The economic analysis was developed for three discount
rates: 10%, 8% and 6%. The first year is adopted for reference year.
The annual operation costs were estimate on the basis of fixed percentages either
of the capital costs (exploitation and maintenance costs), or of the incomes (legal
permit costs).
As one can conclude from the previous table, if the discount rate is set equal to
the internal rate return then the net present value becomes zero, the benefit/cost
ratio one, the average price of the kWh the nominal unitary sale price of the
energy and the payback period the analysis period.
- 188 -
Economic analysis
Table 10.3 Application example.
- 189 -
- 190 -

Guidelines For Design of Small Hydropower Plants

Uploaded by

Copyright:

Available Formats

Guidelines For Design of Small Hydropower Plants

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Guidelines For Design of Small Hydropower Plants

Uploaded by

Copyright:

Available Formats

What are the main chapters covered in the book?

What are the main chapters covered in the book?

What hydrologic studies are discussed in Chapter 4?

What hydrologic studies are discussed in Chapter 4?

TECHNICAL CARD

All rights reserved. No part of this publication may be reproduced, stored in a

(Foreword and reviser)

Legal deposit: July 1999.

The book Guidelines for design of SMALL HYDROPOWER PLANTS

Chapter 5 contains relevant information about the hydraulic design of intakes,

Each author has a good professional experience as consulting engineer in the

Antnio de Carvalho Quintela

Special thanks go to Dr. Robert Gibson, manager from WREAN (Western

I am thankful for the co-operation of several companies, such as ALSTHOM,

1.1- Renewable energy concept

2- GENERAL TYPES OF SMALL

2.1- Classification of hydropower plants

3.1- Power generation and different sector-users

4.1- Introduction. Scope of the studies

6.1- Types of turbines

6.4- Dimensions of turbines

8.3- Control system considerations

9.1- Environmental impact assessment

10- ECONOMIC ANALYSIS

10.1- Introduction. Costs and benefits

flow cross section (m2).

thickness or diameter of a rack bar (m); wave celerity in open channel

inlet and outlet diameter of a turbine runner (m).

available gravity potential head or gross head (m).

rotational speed (r.p.m.); period with simultaneous annual records of

turbine discharge (m3/s); mean daily flow (m3/s); spillway discharge

hydraulic radius (m).

payback period (years).

vapour head (m).

hydraulic inertia time constant (s).

unit absolute flow velocity (-).

angle with horizontal (degree); Coriolis coefficient (-).

relative roughness (-).

1.1- Renewable energy concept

Inexhaustible energetic sources, in spite of being limited or

Nowadays, the policy in most of the countries is devoted to assure additional

Guidelines for design of SMALL HYDROPOWER PLANTS

For environmental protection it must be considered, in each small hydro project,

It saves consumption of fossil, fuel, and firewood.

It offers a decentralised electrification at a low running cost and with

Guidelines for design of SMALL HYDROPOWER PLANTS

Fig. 1.2 Different types of hydropower schemes.

where is the global efficiency, resulting of the multiplication of partial

Guidelines for design of SMALL HYDROPOWER PLANTS

Reservoir: constitute a storage form of the available potential energy and

Fig.1.3 Small old hydroelectric scheme

The hydraulic engineering intervention has a major contribution in these types of

exploitation phases of a project that involves multidisciplinary teams, with

1.3- Environmental considerations

Guidelines for design of SMALL HYDROPOWER PLANTS

because there is no release of carbon dioxide that contributes to ozone depletion

____________________________________General Types of Small Hydropower Plants

2.1- Classification of hydropower plants

The classification based on the power is very important because is an