Professional Documents
Culture Documents
Final Book
Final Book
Faculty of Engineering
Mechanical Department
Following is a list of the team members with a brief description about their work
and contacts:
Ahmed Bahaa El-Deen Negmy
+201120840171
ahmedbahaa3519@gmail.com
Team Leader
Ahmed Tarek Mohamed Ahmed
+201200707769
ad_tk_2008@yahoo.com
Solar systems.
Desalination systems.
Electrolysis.
Hydrogen blending.
Techno-economic
analysis.
Hydraulic modeling
using Synergi Gas.
Regulator station.
Unsteady state
modeling.
Solar systems.
Desalination systems.
Electrolysis.
Hydrogen blending.
Techno-economic
analysis.
+201111967235
shaza_hady1993@yahoo.com
Hydraulic modeling
using Synergi Gas.
Compressor station
modeling using
Synergi Gas
+201096722029
mohamedmostafa9310@gmail.com
Team (HR)
Acknowledgment
First, Thank to Allah, Most Merciful and Most Compassionate. We take this
opportunity to express our deepest and sincere gratitude to our supervisors
Prof. Dr. Hassan Warda, Eng. Sherif Hadara and Prof.Dr. Essam Wahba for
their insightful advice, motivating suggestions, invaluable guidance, help and
support in successful completion of this project and also for their constant
encouragement and advice throughout this work. We feel greatly honored to
have been working under their supervision.
Our deepest thanks are further extended to Eng/ ahmed abd elmoety from
natgas, Eng/ hesham Khalil from Gasco , Gasser Elfaham from Gasco, Dr.Ahmed
Hegazy, Dr.Osama El-Masry, Eng/ Ayman El-Morshedy from Sidi Krir and abd ElMoneam.
We would like to convey our thanks to the teaching staff of the Mechanical
Engineering Department, for their invaluable help and support throughout the
period of Bachelors Degree.
Finally, yet more importantly, we would like to express our deep appreciation to
our families for their prompt help and continuous care which has a great
influence throughout our lives.
Preface
Natural gas is one of the most widely used forms of energy today, a cleanerburning fuel than either gasoline or diesel, used by industrial, residential,
electricity generation, commercial, and transportation sectors in varying
amounts.
Egypts proved reserves of natural gas are the third largest in Africa, after Nigeria
and Algeria, despite that Egypt faces a shortage in Natural gas.
in recent years. For power generation, old power stations which used to run on
oil crude has been converted to natural gas increasingly, this produces less
carbon dioxide than crude oil. Its also more economic than crude oil.
Due to increasing demand of natural gas markets it became a must to optimize
natural gas network operation and to find new energy sources to balance the
network.
Our project target is to optimize the Natural Gas National Grid (NGNG)
Operation. And Also to find new power generation application power to gas.
In this book all the algorithms procedures and considerations will be discussed,
also the final results and recommendations will be stated.
The screenshots from a licensed software package which was used to perform
the Modelling and results will be presented.
Table of contents
Abstract ........................................................................................ 18
Natural Gas Network Components .................................................... 20
1.1 Introduction ............................................................................. 21
1.2Components of network .............................................................. 24
1.2.1 Transmission pipelines ......................................................... 24
1.2.2 Compressor Stations: .......................................................... 50
1.2.3 Regulation and Gas Metering stations..................................... 50
1.2.4 City Gate Stations ............................................................... 51
1.2.5 Valves ................................................................................ 52
1.2.6 Flow Measurements ............................................................. 56
1.2.7 Remote Control System & Telecommunications ....................... 59
Power to Gas ................................................................................. 61
2.1 Introduction: ............................................................................ 62
2.1.1 What is the problem of energy in Egypt? ................................ 62
2.1.2 energy sources available in Egypt:......................................... 62
2.1.3 Power to Gas: ..................................................................... 62
2.2 What is PtG? ............................................................................ 63
2.2.1 Definition: .......................................................................... 63
2.2.2 Advantages of P2G: ............................................................. 63
2.2.3 Components of P2G plant ..................................................... 64
2.3 Solar systems ........................................................................... 64
2.3.1 Solar system in Egypt .......................................................... 64
2.3.2 History ............................................................................... 65
2.3.3 Places of solar plants in Egypt ............................................... 65
2.3.4 Types of solar system .......................................................... 67
2.4 Desalination process .................................................................. 80
2.4.1 Water Classification: ............................................................ 81
2.4.2 Desalination techniques: ...................................................... 81
2.5 Electrolysis ............................................................................... 94
2.5.1 Theory of water electrolysis operation: ................................... 94
2.5.2 Factors affecting on electrical efficiency: ................................ 95
2.5.3 Alkaline electrolysis ............................................................. 96
10
11
12
13
14
15
16
17
Abstract
Egypt has large potential of natural gas. The
consumption has grown unsteadily in recent
years until it has reached the first largest in
Africa which led to the increase in the gap
between production and consumption. One
of the main demands is power station which
is the major consumer of natural gas in Egypt
as they consume about 57 % of the natural
gas produced in Egypt, and 85 % of the
energy consumption in power stations is
from natural gas. Natural gas shortage and
low network pressure are the major
problems which need to be solved and this
can be achieved by studying the gas network
components and power stations and their
effect on the network and searching for new
solutions to operate the network at high
performance.
18
19
Chapter 1
Natural Gas Network
Components
20
1.1 Introduction
Natural gas, like oil and coal, is a fossil fuel
energy source, situated deep under the
earth. The main component of natural gas is
methane (CH4), a colorless, odorless,
although flammable gas. When natural gas is
distributed commercially, it is supplemented
with an artificial odor so as to prevent leaks.
Natural gas is different than the well-known
cooking gas (known as LPG Liquefied
Petroleum Gas), whereas LPG is a product of
petroleum distillation, natural gas is pumped
from the underground, undergoes cleaning
and treatment and then is ready for use.
Similarly, to other fuels and energy sources,
like oil, diesel, fuel oil, coal and others,
natural gas also constitutes a primary energy
source, from which electricity may be
manufactured and mechanical energy and
heat produced. Natural gas is considered as
an environmentally friendly clean fuel,
offering important environmental benefits
when compared to other fossil fuels. The
superior environmental qualities over coal or
crude oil are that emissions of sulfur dioxide
are negligible or that the levels of nitrous
oxide and carbon dioxide emissions are
lower. This helps to reduce problems of acid
rain, ozone layer, or greenhouse gases.
Natural gas is also a very safe source of
energy when transported, stored, and used.
21
Table 1-1
22
23
1.2Components of network
The natural gas pipeline network is a system
that moves gas from the areas of production
to the areas of consumption. Natural gas
pipelines provide three transportation
functions,
gathering,
transmission
(interstate and intrastate), and distribution.
As a system, the pipeline network comprises
a number of components that help move the
gas from production to consumption.
24
25
hourly fluctuations
efficiently.
in
demand
very
26
27
28
29
30
(1.1)
(1.2)
Where:
Where
=density of gas
g= density of gas
m=mass of gas
V=volume of gas
29
(1. 1)
31
(.91)+(.082)+(.023)
29
(1.5)
(1.4)
= kinematic viscosity, St
Where:
=dynamic viscosity, P
= density, kg/m3
Kinematic viscosity is expressed in ft2/s in
USCS units and Stokes (St) in SI units. Other
units of viscosity in SI units include
centipoise (cp) for dynamic viscosity and
centistokes (cst) for kinematic viscosity. The
viscosity of a gas depends on its temperature
and pressure. Unlike liquids, the viscosity of
a gas increases with increase in temperature.
Since viscosity represents resistance to flow,
as the gas temperature increases, the
quantity of gas flow through a pipeline will
decrease; hence, more throughput is
possible in a gas pipeline at lower
temperatures.
1.2.1.8.5 Viscosity
The viscosity of a fluid represents its
resistance to flow. The higher the viscosity,
the more difficult it is to flow. Lowerviscosity fluids flow easily in pipes and cause
less pressure drop. Liquids have much larger
values of viscosity compared to gases. For
example, water has a viscosity of 1.0
centipoise (cp), whereas viscosity of natural
gas is approximately 0.0008 cp. Even though
the gas viscosity is a small number, it has an
important function in determining the type
of flow in pipelines.
32
Where
P = absolute pressure of gas, psia
(1.13)
Where
P = absolute pressure of gas, psia
T = absolute temperature of gas, R
Tr = reduced temperature, dimensionless
(1.12)
(1.2)
33
Is defined as follows:
(1.16)
34
2
2
2
2
+ = +
(1.18)
Where
Hp is the equivalent head added to the fluid
by a compressor at A and
35
= 1.1494 103 ( )
12 22
.5
2.5
Where
Q=gas flow rate, measured at standard
conditions, m3/day
f=friction factor, dimensionless
2. Colebrook-White equation
4. AGA equation
5. Weymouth equation
6. Panhandle A equation
7. Panhandle B equation
8. IGT equation
9. Spitzglass equation
36
= 5.747 104 ( )
12 22
.5
) 2.5 (1.20)
(1.21)
= 5.747 104 ( )
12 22
.5
) 2.5 (1.22)
Where
( 1)
(1.23)
37
= .0684 (
2 1
) (1.24)
Where
s= elevation
dimensionless
adjustment
parameter,
(1.25)
= 1 1+ 2 2 +
(1.26)
38
(1.27)
= Z1RT1 (1.32)
1 = ( ) (1) 1
(1.28)
(1.33)
= 14.7349 ( 2 ) ( ) ( ) (1.34)
Where
39
P = pressure, kPa
= 10029
(1.36)
Where
100
(1.35)
Where
max = maximum or erosional velocity, ft/s
= gas density at flowing temperature,
lb/ft3
Since the gas density r may be expressed in
terms of pressure and temperature,
usingthe gas law Equation, the maximum
(1.37)
Where
Re = Reynolds number, dimensionless
40
= .0004778 ( ) ( ) (1.39)
Where
Pb = base pressure, psia
= .5134 ( ) ( ) (1.40)
Where
Pb = base pressure, kPa
(1.41)
Where
=Fanning friction factor
41
64
(1.42)
(1.43)
Where
e = absolute or internal roughness of pipe,
in.
D = pipe inside diameter, in.
1.2.1.8.12.1 Colebrook-White Equation
The Colebrook-White equation, sometimes
referred to simply as the Colebrook
equation, is a relationship between the
friction factor and the Reynolds number,
pipe roughness, and inside diameter of pipe.
The following form of the Colebrook
equation is used to calculate the friction
factor in gas pipelines in turbulent flow.
(1.44)
Where
f = friction factor, dimensionless
42
(1.45)
For turbulent flow in smooth pipes similarly,
for turbulent flow in fully rough pipes, with
Re being a large number, f depends mostly
on the roughness e and, therefore, the
friction factor equation reduces to
(1.46)
43
F = transmission factor
(1.49)
Since Re, f, and F are dimensionless, as long
as consistent units are used for both e and D,
the modified Colebrook equation is the
same, regardless of the units employed.
(1.47)
Therefore,
(1.48)
Where
f = friction factor
44
length of the
consideration.
pipe
section
under
(1.53)
Note: The drag factors above are based on
40-ft joints of pipelines and mainline valves
at 10-mile spacing.
(1.50)
Equation 1.2.50 is also known as the Von
Karman rough pipe flow equation.
For the partially turbulent zone, F is
calculated from the following equations
using the Reynolds number, a parameter Df
known as the pipe drag factor, and the Von
Karman smooth pipe transmission factor Ft:
(1.51)
And,
Where
(1.52)
= 433.5 ( )
45
12 .22
0.961 .
) . 2.667 (1.54)
Where
1.0788 2
(1 22 )
( )
= 435.87 0.8539
. .
2.2.618
(1.56)
Where
Q = volume flow rate, standard ft3/day
(SCFD)
( ) (12 22 )
= 3.7435 103
. .
2.667
(1.55)
Where
Q = gas flow rate, standard m3/day
46
= 4.5969 103
1.0788 2
(1 22 )
( )
0.8539 . .
2.2.618
(1.57)
Where
1.02
= 1.002 102 ( )
0.51
12 . 22
( 0.961
)
. 2.53
. .
(1.59)
Where
Q = gas flow rate, standard m3/day
E = pipeline efficiency, a decimal value less
than 1.0
Tb = base temperature, K (273 + C)
1.02
= 737 ( )
12 . 22
( 0.961
)
. .
. 2.53
(1.58)
Where
47
(IGT) Equation
Tb
= 1.28822 103 E ( )
Pb
0.555
2
s 2
P1 e . P2
( 0.961
)
. D2.667
G
. Tf . Le .2
(1.61)
Where
Q = gas flow rate, standard m3/day
E = pipeline efficiency, a decimal value less
than 1.0
Tb = base temperature, K (273 + C)
Pb = base pressure, kPa
Tf = average gas flow temperature, K (273 +
C)
P1 = upstream pressure, kPa (absolute)
P2 = downstream pressure, kPa (absolute)
Le = equivalent length of pipe segment, km
= gas viscosity, Poise
Tb
= 136.9 E ( )
Pb
0.555
P12 es . P22
( 0.961
)
G
. Tf . LeZ
. D2.667
(1.60)
Where
Q = volume flow rate, standard ft3/day
(SCFD)
E = pipeline efficiency, a decimal value less
than 1.0
48
49
50
Gas Filtration
Flow metering
Heating of the fuel
Reduction the fuel pressure
51
1.2.5 Valves
Valves and flow measurements used on gas
pipelines. The design and construction codes
for valves, materials of construction, and
application of the different types of valves
and their performance characteristics will be
explained. The importance of flow
measurement in a gas pipeline and the
accuracy of available instruments, codes,
and standards used will be discussed.
Various American National Standards
Institute (ANSI), American Petroleum
Institute (API), and American Gas Association
(AGA) formulas used in connection with
orifice meters will be reviewed. Since a small
error in mea surement in gas flow in a
pipeline can translate to several thousand
dollars of loss of revenue, it is important that
industry strives to improve upon
measurement methods. Accordingly, gas
transportation companies and related
industries have been researching better
ways to improve flow measurement
accuracy.
52
Gate valve
Ball valve
Plug valve
Butterfly valve
Globe valve
Check valve
Control valve
Relief valve
Pressure regulating valve
53
54
55
56
57
58
59
60
Chapter 2
Power to Gas
61
2.1 Introduction:
2.1.1 What is the problem of energy
in Egypt?
In Egypt, as everywhere, energy plays a
substantial role in countrys economic
development.
On the other hand, Egypt as a developing
fast growing country suffers from rapid
annual population growth currently at a
rate of 1.68%. Egypts population is
expected to reach 110 million by 2030 and
128 million by 2050.
Such a fast population growth along with
other environmental challenges is overstraining the limited energy resources of
the country.
62
63
Photovoltaic plants.
Concentrating solar power plants (CSP)
(Parabolic trough).
64
2.3.2 History
The first solar power plant in Egypt was
built
in
MAADI
in
1913!
it is not a hoax because if you check the
history of the solar energy you will find that
solar energy is much older than we think, it
is not a 21th century idea after. It is 20th
century idea.
65
2) NORTH COAST
3) TIBA-QENA
66
2.3.4.1.1 Components of
photovoltaic plant
1) Solar cell: is the basic component of
any photovoltaic, or PV, system. A
device like a calculator requires only a
single cell, but most devices or
commercial
applications
needs
multiple cells.
A solar cell is a solid-state electrical
device (p-n junction) that converts the
energy of light directly into electricity
(DC) using the photovoltaic effect.
67
(p-n) junctions:
Depletion region:
A depletion region forms instantaneously
across a p-n junction. Electrons and holes
diffuse into regions with lower
concentrations of electrons and holes,
much as ink diffuses into water until it is
uniformly distributed.
By definition, N-type semi-conductor has
an excess of free electrons compared to
the P-type region, and P-type has an excess
of holes compared to the N-type region.
When N-doped and P-doped pieces of
semiconductor are placed together to
form a junction, electrons migrate into the
P-side and holes migrate into the N-side.
Departure of an electron from the N-side
to the P-side leaves a positive donor ion
behind on the N-side, and likewise the hole
leaves a negative acceptor ion on the Pside.
68
Photovoltaic effect:
When a photon of light is absorbed by
a valence electron of an atom, the
energy of the electron is increased by
the amount of energy of the photon. If
69
4) Tracking system
To increase solar yield and electricity
production, solar trackers and
concentrator systems are sometimes
used. Trackers direct solar panels or
modules toward the sun.
These devices change their orientation
throughout the day to follow the suns
path to maximize energy capture.
In photovoltaic systems, trackers help
to minimize the angle of incidence (the
angle that a ray of light makes with a
line perpendicular to the surface)
between the incoming light and the
panel, which increases the amount of
energy the installation produces.
70
2.3.4.1.3 PV technologies
There are a wide range of PV cell
technologies on the market today, using
different types of materials, and an even
larger number will be available in the
future.
PV cell technologies are usually classified
into three generations, depending on the
basic material used and the level of
commercial maturity:
71
2) Cadmium-Telluride (CdTe).
Figure
(CIGS)
2.15 Copper-Indium
Gallium-Diselenide
72
73
2.3.4.1.5 Advantages of PV
system
1. PV panels provide clean, green energy.
During electricity generation with PV
panels there is no harmful greenhouse
gas emissions thus solar PV is
environmentally friendly.
2. Solar energy is energy supplied by
nature, so its free and abundant. Solar
energy can be made anywhere there is
sunlight.
3. Solar Panels cost is currently on a fast
reducing track and is expected to
continue reducing for the next years.
Solar PV panels has a highly promising
future both for economic viability and
environmental sustainability.
4. Compared to other renewable green
energy technologies, PV systems
require minimum maintenance; only
minor checking of cable connections
and a basic-regular cleaning of the
panel surfaces is adequate to keep
them operational for several years.
5. PV panels have no mechanically
moving parts, except in cases of
sun-tracking
mechanical
bases;
consequently they have far less
breakages or require less maintenance
than other renewable energy systems
(e.g. wind turbines).
6. PV panels are totally silent, producing
no noise at all; consequently, they are
- Methods of cleaning:
- cleaning warning:
74
4.
5.
6.
7.
2.3.4.1.7 Photovoltaics
calculations:
Power calculations
The global formula to estimate the
electricity generated in output of a
photovoltaic system is:
E = A. r. H. PR
Where,
75
( )
= (
PV array sizing:
To calculate the PV area:
PV Area =
Where,
EL: Average
(KWh/day)
daily
load
energy
76
77
2. Mirrors.
The parabolic mirrors that reflect solar
energy to the HCE use in this study will
be low-iron laminated or monolithic
glass mirrors. These mirrors are known
to be reliable components that show
no long-term degradation in reflective
quality. Reflectivity of the mirrors is
closely monitored to determine the
frequency and method of mirror
washing. Regular mirror washing has
been proven to maintain mirror
reflectivity at installed conditions.
78
2.3.4.2.2 Configurations of
parabolic trough systems
Using only parabolic trough panels to
get heat energy, some of that heat
79
2.3.4.2.4 Disadvantages of
parabolic trough systems
- High Costs, Solar Thermal Energy costs at
least 3.5 Euro /watt. However, these costs
are too high as Solar PV already costs 2.5
Euro /watt and it will have its costs
reduced by 5% in the next 10 years making
it attain half the cost of Solar Thermal
Technology
by
2020.
- Water Issue, Solar Thermal Plants use lots
of Water which is Major Problem in Desert
Areas which raises the cost of CSP projects
(concentrating solar power) too much.
While using Sea Water has been proposed
it. Remains to be seen if it possible to
implement this solution as this would
imply building Plants very near the Coastline.
- Limited Locations and Size Limitations,
Solar thermal energy can only be built in
places which have the high amount of solar
radiation. They can be built in deserts
mostly and require a large land area. This
means it is not possible to build them in
80
81
82
One method used recently for the pretreatment of seawater before directed to
RO modules is nano-filtration (NF). NF is
primarily developed as a membrane
softening process which offers an
alternative to chemical softening.
The main objectives of NF pre-treatment
are:
1. Minimize particulate fouling of the RO
membranes by removal of turbidity and
bacteria.
2. Prevent scaling by removal of the
hardness
ions.
3. Lower the operating pressure of the RO
process by reducing the feed-water total
dissolved solids (TDS) concentration.
83
Figure. 2.21. Membrane module types, a) plate and frame, b) tubular, c) spiral wound and d) hollow fiber modules
2.4.2.1.4 RO calculations
- To determine the rate of water transport
across the membrane:
NAw = L ( )
Where,
= CRT
where,
p:
the
difference
transmembrane
pressure
- L is given by:
L=
where,
D: the water diffusivity in the membrane
S: the water solubility in the membrane
V: the molar volume of water
R: the ideal gas constant
T: the ambient temperature
84
2.4.2.2.1.2 Disadvantages of
MSF:
A disadvantage of MSF is that precise
pressure levels are required in the
different stages and therefore some
transient time is required to establish the
normal running operation of the plant.
2.4.2.2.1.3 Configurations of
MSF:
The MSF desalination plants are divided
into two models, Once-through MSF (MSFOT) and Brine recirculation MSF designs
(MSF-BC).
85
formation.
- The feed seawater (Mf) enters the
flashing stages. The flashing process takes
place due to a decrease in the stage
saturation temperature that causes the
reduction in the stage pressure.
- In each stage, the flashed off vapor
condenses on the outside surface of the
condenser tubes, where the feed seawater
(Mf) flows inside the tubes from the cold to
the
hot
side
of
the
plant.
- The condensed fresh water vapor outside
the condenser tubes accumulates across
the stages and forms the distillate product
stream (Md). This stream cascades in the
same direction of the flashing brine from
stage to stage and is withdrawn from the
last
stage.
- The vapor flashed off in the flash chamber
is separated from the condenser by a
demister which intercepts brine droplets
entrained in the flashing vapor resulting in
86
87
2.4.2.2.1.3.3 Comparison
between MSF processes:
The most common multistage flash
desalination process is the brine recycle
system
(MSF-BC).
As for the once through system (MSF-OT)
its use is found on a very limited scale.
88
89
Stage Model:
- The mass conservation equation for
the brine in the flashing pool is given
by:
()
() ()
()
= () () ()
+ () ()
() () ()
= (1) () ()
, () = () () ()
- The balance of the dissolved solids
reads:
()
()
= (1) (1)
() ()
- The energy balance for brine pool is
given by:
()
() ()
()
(() () ) (() (+1) )
=
(() () )
ln[
]
(() (+1) )
- The mass balance for the condenser tubes
can be written as:
() () () () ()
() = (+1)
()
= (1) + () ()
()
() ()
()
()
()
() ()
- The mass balance for the vapor space
produces:
+
() ()
()
= () + () ()
()
- The energy balance for the vapor space is
as follows:
90
() () ()
() ()
Brine Heater Model:
- The mass and energy balance
equations are given by:
91
92
2.4.2.2.2.3 Disadvantages of
MED:
- They are expensive to build and operate
as energy consumption is particularly high.
- The product water is at an elevated
temperature and can require cooling
before using as a potable water.
- The recovery ratio is low, although not as
low as for MSF.
2.4.2.2.3 Comparison between
MSF and MED:
MSF and MED are both thermal
desalination
processes.
The
way
evaporation and heat transfer take place in
these processes are completely different.
In MED process, vapor is originated from
the falling lm of seawater which brought
in contact with the heat transfer surface. In
MSF process, seawater is rstly heated
within the tubes by convection then vapor
is initiated by ashing as it streams in each
stage of reduced pressure. MSF system
normally consists of three sections: heat
recovery, heat rejection and brine heater.
93
2.5 Electrolysis
Water electrolysis is currently one of the
most promising methods for hydrogen
production. Hydrogen is produced via
electrolysis by passing direct current (DC)
through two electrodes in water. The
water molecule is split, producing oxygen
at the anode (positive electrode) and
hydrogen at the cathode (negative
electrode). The technique can be used to
make hydrogen fuel and breathable
oxygen; though currently most industrial
methods make hydrogen fuel from natural
gas instead.
Electrolysis reactions:
Cathode Electrode reaction:
2 2 + 2 2 + 2
Anode Electrode reaction:
2 2 + 2 +
Hence,
the
2 2 2 +
total
1
2
1
2 2
reaction
is:
2 + 2
94
95
96
97
2 0.5 2 + 2 + 2
Hydrogen evolution reaction (HER):
2 2 + 2
2 + 2 ()
98
99
Alkaline
electrolysis
is
the
cheapest
and
most
reliable
technology.
Alkaline
electrolyser
can be operated between 20 and
100% of the design capacity, and
100
101
Cathode reaction:
High pressure
12000 psi).
operation
(300-
102
Cathode reaction:
2 +2 = 2 + 2-
2- = 2 + 2
103
other types:
In the best low-temperature electrolytic
cells operating under normal pressure, the
average cell voltage is 1.7 V. It
corresponds to a power input of 4.07 KWh
/ Nm3H2. However, this number does not
reflect all the power inputs for hydrogen
production
in
a
low-temperature
electrolytic cell. Electrolyzers based on
low temperature electrolytic cells require
very pure water (distilled or even bidistilled), so for a proper analysis of energy
consumption of such devices it is necessary
to take into account the power inputs for
preparation of pure water. To prepare a
sufficient amount of water for production
of 1 m3 of hydrogen, it is necessary to
spend no less than 0.73 kWh in the case of
one step distillation and about 1.5 kWh in
the case of double distillation. Electrical
energy is generally used for water
distillation, so, the real power input in a
low temperature electrolytic cell is 4.8 to
6 KWh / Nm3H2.
104
1 H2O + Electricity = 1 H2 + 2 O2
This reaction is endothermic and the
required energy by the process is:
= +
where, (J/mol) is the change of
enthalpy, (J/mol) is the Gibbs free
energy and represents the minimum
96485.3365
ln
2
Where,
0 : standard equilibrium potential
105
0
2
Where,
0 is increment of Gibbs free energy
F is Faraday constant (96519.4
coulomb/mol)
EO is 1.226 V under 298 K and 1 atm. This
is slightly higher than the decomposition
voltage of pure water under the same
condition.
Since is theoretical equilibrium
potential, actual cell voltage needed to
continue the reaction is higher with the
addition of ohmic loss of electrolyte and
diaphragm and overvoltage by electrode
reaction.
_ =
Where,
141.86 (
= 1.4819 V
E is electrolysis voltage
is reversible voltage
is overvoltage
_ =
is ohmic loss.
119.93 (
1.253
of hydrogen.
- The amount of hydrogen produced
can be calculated by:
= ( )
2
where p (MWh) is the available energy
and is the electrolyser efficiency
2 =
106
hydrogen
2 =
is
calculated
by:
W=
2 ()
nRTln( 2 )
) is the
.
2
(MW)
107
108
Table 2.8 Summary of the key operational parameters of Alkaline, PEM, and Solid Oxide Electrolysis.
109
110
111
2.7 Implementation of
Power to Gas (P2G) around
the world:
2.7.1 Germany experience:
Hydrogen production (70 percent energy
conversion efficiency) is just the first step
in a process which will result in the
production of synthetic gas (i.e. methane 55 to 65 percent energy conversion
efficiency) which can replace natural gas.
Existing energy storage technologies like
compressed-air energy Storage (CAES) and
pumped storage hydroelectricity certainly
have a role to play in the future, but alone
they are not able to meet the expected
surpluses of solar and wind energy or
expected energy shortfalls in periods of
low sunshine and wind (future renewable
energy-based domestic power supply
requires a storage capacity between 20 to
40 terawatt hours, current energy storage
capacity is approximately 0.04 terawatt
hours). At present, a compensating
capacity level of 17 GW is already required
to maintain domestic supply and demand
balance. This figure is set to rise to 28GW
in 2025. To date, balancing power has
largely been provided by hydro pumped
storage
and
gas power
plants.
112
113
114
2.8 Techno-Economic
Analysis
2.8.1 Technical analysis:
Based on the Egyptian national natural gas
network data, the capacity of NG is
1083000 cubic meter per day in WDGC,
and a 3% Hydrogen blending is determined
as it is a safe injection blending into the gas
network. Hence, the rate of hydrogen
production needed from the P2G is 0.03
multiplied by 1083000 which equals to
32490 m3/day.
115
Mass Balance:
2 =
1
2
2 .
: pump efficiency 80 %
Where,
2 : The molar mass of hydrogen (2.015
).
).
(E8) =
Energy Balance:
1 2
2 .
where:
and
panels
Where:
E=
Energy
(kWh
/year)
A = Total solar panel Area (m)
H = Annual average solar radiation on
tilted panels (3200 KWh /m2. year)
: Photovoltaics panel efficiency
: assumed be 0.765
system:
2: Number of electrons
power (E6) =
For
Photovoltaic
output from PV = E6 + E7+ E8
where:
: The osmotic pressure of solution (Pa)
C: The concentration of solution
(mole/3 )
116
117
Mass Balance:
For
PEM
electrolyser:
1
2 =
2
2 .
. 78 %
For
MSF
process:
5 = 2 + 3 + 4
- For CSP system:
= +
+ 5 5 5 = 4 4 4
+ 3 3 3 + 2 2 2
= + +
9 =
Energy Balance:
For
PEM
For
, assume =0.23
7 =
( )
electrolyser:
2
2
2 .
Assume = 2 volt
-
=
MSF
process:
- Input Electricity = 6 + 8
118
Mass Balance:
- For PEM electrolyser:
0
2 =
2
2 .
Output power = 1 + 8
8 =
5 = 2 + 3 + 4
2 + 3 = 0.7 5
Energy Balance:
- For PEM electrolyser:
=
& =
Input power =
0 2
2 .
Input
Assume = 2 volt
power
Assume = 0.23
Output power ( 6 ) =
(
)
Assume = 300
= 240
= 5 5 5 +
( )
119
Mass Balance:
- For SOEC electrolyser:
0
1 =
2
2 .
.
For
MSF
4 = 1 + 2 + 3
1 + 2 = 0.7 4
Energy Balance:
=
.
-
54
%
process:
Output power = 7 + 8
7 =
Input power =
0 2
& =
2 .
10
=
()
10
panels
= hS/H hsat.
Input power =
Output power ( 5 ) =
(
)
99 %
P.P =
All other things being equal, the better
investment is the one with the shorter
payback period. For example, if a project
costs $100,000 and is expected to return
$20,000 annually, the payback period will
be $100,000/$20,000, or five years.
1. It
121
and,
therefore,
does
not
measure profitability.
2. It ignores the time value of money.
NPV = =1 (1+)
Where,
Ct = net cash inflow during the period t
Co = total initial investment costs
r
=
discount
rate
t = number of time periods
A positive net present value indicates that
the projected earnings generated by a
project or investment (in present dollars)
exceeds the anticipated costs (also in
present dollars).
122
123
2.8.3.6 Inflation:
A sustained, rapid increase in prices, as
measured by some broad index (such
as Consumer Price Index) over months or
years, and mirrored in the correspondingly
decreasing purchasing
power of
the currency. When a dollar buys less than
you would expect it to, we call that
inflation. Inflation is caused by a variety
of factors, but most of them are related to
interest and debt. When the Federal
Reserve bank raises interest rates, it
causes the dollar to inflate. There is
more money in the system, so every dollar
is worth just a little bit less.
124
125
finance
government
activities.
In economics, taxes fall on whomever pays
the burden of the tax, whether this is the
entity being taxed, like a business, or the
end consumers of the business's goods.
2.8.3.13 Depreciation:
Depreciation is the process by which a
company allocates an asset's cost over the
duration of its useful life. Each time a
company prepares its financial statements,
it records a depreciation expense to
allocate a portion of the cost of the
buildings, machines or equipment it has
purchased to the current fiscal year. The
purpose of recording depreciation as an
expense is to spread the initial price of the
asset over its useful life.
2.8.3.14 Loan or Debt:
Any investment project takes a loan from a
bank which afford the total cost of the
project with the investor. Generally, that
loan is a percentage of the total cost and
that percentage is used to be 70 %.
2.8.3.15 Taxes:
126
127
128
Sensitivity
of
CSP
Panel
price:
When the price of PEM electrolyser is
1,000,000,000 $
129
130
Sensitivity of PV price:
131
132
Sensitivity of PV price:
CSP Panels price dont have a significant
effect, hence assume the worst case of its
price which is 3.5 $ / watt.
133
2.8.5 Conclusion:
The most effective parameters of the
Internal Rate of Return (IRR) are the
electrolyser price and the hydrogen price
per Kg.
134
Hydrogen price
110
$20
$80
PV Price
-20%
800
$3
-10%
0%
$0.5
10%
20%
135
30%
40%
50%
0.4
Alkaline
Hydrogen price
PV cost
0.3
0.2
0.1
-0.1
-0.2
40%
60%
80%
100%
120%
140%
136
160%
180%
Hydrogen price
1000
$50
$120
CSP
-20%
300
$3.5
-10%
0%
10%
$1.0
20%
137
30%
40%
50%
0.4
PEM
Hydrogen price
0.3
CSP
0.2
0.1
-0.1
-0.2
40%
60%
80%
100%
120%
140%
160%
180%
Hydrogen price
1000
$44
$95
PV
-20%
$3
-10%
0%
$0.5
10%
20%
138
30%
40%
50%
0.4
PEM
Hydrogen price
0.3
PV
0.2
0.1
-0.1
-0.2
40%
60%
80%
100%
120%
% Change in Input Value
140%
139
160%
180%
800
Hydrogen price
$50
280
$95
PV
$3
-20%
-10%
0%
$0.5
10%
20%
30%
40%
0.4
PEM
Hydrogen price
PV
0.3
0.2
0.1
-0.1
-0.2
40%
60%
80%
100%
120%
% Change in Input Value
140%
140
160%
180%
50%
141
Chapter 3
Power Stations
142
3.1 Abstract
Egypt has been an industrialized nation.
But, energy shortages have been a
challenge to Egyptian industry.
143
144
1. Cycle.
2. Application.
3. Shaft Configuration.
145
146
147
1. The
system is dependent on
external means as considerable
quantity of cooling water is
required in the pre-cooler.
2. Higher internal pressures involve
complicated
design
of
all
components and high quality
material is required which
increases the cost of the plant.
3. The response to the load variations
is poor compared to the open-cycle
plant.
4. It requires very big heat exchangers
as the heating of workings fluid is
done indirectly.
is
3.4.2.5 Microturbines
turbines
The frame units are the large power
generation units ranging from 3 to 480
MW in a simple-cycle configuration,
with efficiencies ranging from 30% to
48%.
148
149
02
01
3.7 Thermodynamic
Analysis of Gas Turbines
=(
02 1
01
Considering
efficiency:
02 = 01 +
= ( )
compressor
isentropic
02s 01
= 12 (02 -01 )
03 =
. + . 23 . 02
. 34
34 03 23 02
.
= 34 03 23 02
(3-4) Isentropic expansion process
= ( )
= ( )
150
= g 34 (03-04 )
=
(4-1) Constant-pressure heat rejection
= 41 (04 01 )
Fig. 3-8 Brayton cycle with Regeneration
Cycle Efficiency:
=
3.8.1 Regeneration
The most common cycle improvement is
regeneration where exhaust heat is used
to increase the temperature of
compressed air before combustion. This is
accomplished by installing a heat
exchanger in the exhaust to preheat the air
between the compressor and the
combustors.
3.8.2 Intercooling
In some gas turbine arrangements, the
compression of the inlet air is done in two
stages with the air being cooled between
the stages in a heat exchanger
(intercooler).
Since
isothermal
compression (compression without an
increase in temperature) takes less work
than adiabatic compression (compression
where no heat is removed so that the air
temperature increases), more of the
151
3.8.3 Reheat
In addition to compressing the air in two
stages and intercooling between these
stages, the gas turbine plant may also be
arranged to expand the hot gases in two
stages, with the gases being reheated
between the stages. The gases are
expanded first in a high-pressure turbine
and then reheated before entering a lowpressure turbine.
152
However,
increased
the
ambient
temperature an increased the air to fuel
ratio caused increased of specific fuel
consumption and heat rate as shown in
Figures 3-15 and 3-16.
153
154
3.9.4 Conclusion
The simulation result from the modeling of
the influence of parameter showed that
compression ratio, ambient temperature,
air to fuel ratio and turbine inlet
temperature effect on performance of gas
turbine power plant. The results were
summarized as follows:
The
compression
ratios,
ambient
temperature, air to fuel ratio as well as the
isentropic efficiencies are strongly
influence on the thermal efficiency of the
gas turbine power plant
The variation of thermal efficiency at
higher compression ratio, turbine inlet
temperature and ambient temperature are
very important.
The thermal efficiency and power output
decreases linearly with increase of
ambient temperature as well as the air to
fuel ratio.
The specific fuel consumption and heat
rate increase linearly with increase of both
ambient temperature and air to fuel ratio.
155
2. Programmable
(PLC)
3. Distribution control system (DCS)
Governor control
Lubrication system
Fuel system
Starting system
Cooling system
Instrumentation and control
1.
Sequencing
2. Control
3. Protection
4. Display
156
2. The
3. The
4. Fire
occurred
enclosure
in
the
turbine
157
using accelerometers
mounted on the engine case if antifriction bearings are used (applies
to most aero-derivatives)
3.10.1.5 Instrumentation
The following instrumentation is normally
provided for control and monitoring of a
gas turbine:
1.
7. Vibration
using
eddy-current
displacement probes if journal or
tilt-pad bearings are used (applies
to most heavy-duty gas turbines)
2. Air
inlet temperature in C,
preferably after the intake filter but
often ambient air temperature.
8. Bearing temperatures in C if
journal
or tilt-pad bearings are used
(applies to most heavy-duty gas
turbines)
3. Differential
9. Fuel
gas flow,
temperature
4. Compressor
discharge pressure
measured at the exit of the
compressor
and
before
combustion, in kPa
5. Exhaust
gas temperatures in C,
usually measured at multiple
circumferential points and as an
average after the first or second
158
pressure
pressures
and
10.
Oil
system
temperatures
and
11.
1.
Sequencing
The startup and shutdown of a gas turbine
may occur automatically if predetermined
conditions occur. For example, a backup
power generation unit may start if there is
an increase in demand and a compressor
may start if the pressure drops in a process.
Often, operators monitoring the overall
process will initiate a manual start. Once
the start or shutdown is initiated, the
sequencing is almost always automatic.
3.10.2.2 Lubricants
Synthetic oils with low viscosity and high
viscosity index are recommended for gas
3.10.2.3 Functions lube oil systems
1.
2.
159
2. Radial
160
161
162
Hydraulic System
7. Taking
8. Calibrating
instrumentation
testing protective devices
2. Monitoring
oil pressures
temperatures (usually daily)
3. Checking
and
alarm
4. Topping
5. Changing
and
the
163
Alcohol
Kerosene
Gasoline
Stove Oil
Diesel
Heavy residual oil
164
165
Compressor data
Type: centrifugal compressor
Actual inlet flow: 364 to 4338 m3/h
Inlet temperature: 300C
Discharge: 37.5 bar
Main motor speed: 2980 RPM
Maximum speed of pinion shaft: 27000 RPM
Lube system: 2.2 bar & 540C
Oil filling: 1250 liter
Cooling water supply pressure: 5 bar
166
2. Monitoring
pressures
and
temperatures (usually daily)
3. Changing fuel filters at the required
differential pressure
4. Calibrating
instrumentation
and
testing protective devices
167
5.
6.
7.
8.
9.
168
3.11 Protection
2. Exhaust
3. Excessive vibration
4. Excessive
oil
temperatures
and
bearing
6. Excessive
169
170
6. Bore
On-Line Waterwash
1.
4. Perform
7. Check
TROUBLE SHOOTING
5. Check
1.
171
Check
the
pressure
instrumentation (gauge or
transmitter)
o The
pressure regulator
need adjusting
o The lube oil filter need
replacement due to hi
deferential pressure
o There is blockage in main
lube oil header
o There is major leak in the
lube oil line
o The lube oil pump is not
pumping
due
to
(cavitations- failure- pump
gear box failure)
2. High oil temp. alarm
o Check
the
temp.
Instrumentation
(gaugeRTD)
o Oil cooler is not operating
o Oil cooler is not clean
o The temp. Control valve is
internal leaking (by-passing
the cooler)
o Problem in one or more
lube oil bearing
o The lube oil needs to be
replaced due to aging
3. Ignition failure (fail to ignite - flame
failure)
o The fuel pressure is not
enough
o The igniter is bad or not
cleaned
o The flame detector is not
working
o Bad fuel metering valve
4. The engine fails to reach the self
sustain speed
o The fuel pressure is not
enough
o The fuel metering valve is
not working good
172
1.
2.
Start initiation
3.
4.
Warm-up
5.
Loading
173
174
2.
3.
175
3.14.5 Loading
3.14.4 Warm-up
176
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
177
Hours
Load (MW)
44133
149
27
43479
150
27
43966
150
28
43868
149
28
43786
150
28
43431
149
27
43255
150
27
44219
149
27
44127
150
27
10
44380
149
26
11
44105
149
26
12
43629
149
26
13
43960
149
27
14
44283
148
27
15
44101
149
27
16
43951
149
27
17
44463
149
27
18
44103
149
27
19
44255
149
26
20
59435
219
25
21
59690
218
25
22
59417
219
25
23
59129
218
26
24
59248
219
26
178
30
Gas Pressure
25
20
15
10
5
0
0
10
15
20
25
Hour
Fig. 3-41 gas pressure (bar) with time (hours)
10
15
179
20
25
200
150
100
50
0
0
10
15
180
20
25
Chapter 4 | Compressors
Chapter 4
Compressor Stations
181
Chapter 4 | Compressors
4.1 Introduction
Compressor station is one of the main
components of the natural gas national
grid, as metering stations, pipes, valves,
and reduction station. Compressor
stations have two main purposes, they can
act as a lifting station or booster station.
The lifting station is located near the fields
e.g. Rashpetco compressor station, or it is
located at the processing plants e.g. WDGC
(western desert gas complex), the lifting
stations act as a source of gas to the
national grid.
The boosting compressor station is located
through the network e.g. Dahshour, or it is
located at the consumers e.g. Al-Akaba
compressor station to rise the gas pressure
in order to provide the gas to the
consumers with the required pressure.
4.2 Compressor
4.2.1 Introduction
A compressor is a device used to increase
the pressure of a compressible fluid. The
inlet pressure level can be any value from
a deep vacuum to a high positive pressure.
The discharge pressure can range from sub
atmospheric levels to high values in the
tens of thousands of pounds per square
inch. The inlet and outlet pressures are
related, corresponding with the type of
compressor and its configuration. The fluid
can be any compressible fluid, either gas or
vapor, and can have a wide molecular
weight range. Recorded molecular weights
of compressed gases range from 2 for
hydrogen to 352 for uranium hexafluoride.
Applications of compressed gas vary from
consumer products, such as home
refrigerators,
to
large
complex
petrochemical plant installations.
4.1.1 Definition
A compressor station is a facility which
provides the motive power for the
transportation of natural gas from one
location to another. Natural gas, while
being transported through a gas pipeline,
needs to be constantly pressurized at
certain distance intervals (usually
compressor stations are located 40 to 100
miles apart along the pipeline, but are
sometimes even farther apart). The gas
enters the compressor station with low
pressure and velocity, looking for an
energy source to move it further down the
pipeline.
4.1.2 Components
1) Compressor: Re-pressurize the volume
of gas flowing through the pipeline.
2) Lubrication system: Lubricate, cool and
protect the moving parts.
3) Cooling system: When natural gas is
compressed, its pressure and temperature
increase. The gas may be cooled before its
return to the pipeline to protect the
182
Chapter 4 | Compressors
Compression Methods
the
family
of
183
Chapter 4 | Compressors
184
Chapter 4 | Compressors
185
Chapter 4 | Compressors
186
Chapter 4 | Compressors
187
Chapter 4
compressors
number of stages.
188
Chapter 4
compressors
189
Chapter 4
compressors
4.2.3.5 SR Compressors
These compressors are suitable for
relatively low pressure services. They have
the feature of having several shafts with
overhung impellers. The impellers are
normally open type, i.e., shroudless, to
achieve high tip speeds with low stress
levels and high pressure ratios per stage.
Each impeller inlet is coaxial whereas the
outlet is tangential. These compressors are
generally employed for air or steam
compression, geothermal applications etc.
(see Fig. 4.20).
190
Chapter 4
compressors
Ep,1-2=h2-h1
(4.1)
Ek,1-2 =
(4.2)
Turbocompressors
(more
briefly,
compressors) constitute a special category
of process machines. They operate with
compressible fluids and are characterized
by an appreciable increase in the density of
the fluid between the first and the last
compression stages. The compression
process is frequently distributed among
several stages, a term used to indicate an
elementary system composed of mobile
blading, in which the fluid acquires energy,
h0,4 = h4 +
(4.3)
191
Chapter 4
compressors
HP =
(4.7)
p = Cos tan te
Hp=
Zo R Too [( )
With
(4.10)
p = Cos tan te
(4.11)
(4.4)
H s=
Zo R Too [ ( )
We also have
Hs =
]
(4.9)
HR =
(4.8)
(4.5)
(4.12)
HR = ( h04 h01)
(4.6)
192
Chapter 4
compressors
p =
s =
[ ( )
p =
() ( )
( )
(4.15)
= ( ) (4.18)
The work transferred through the blading
per mass unit of fluid processed is thus
given by:
Wx = /m = ( )
s =
(4.19)
The first principle of thermodynamics
establishes that the work per mass unit is
equal to, for an adiabatic flow, the
variation in total enthalpy. We thus obtain:
, = h2 h1
(4.17)
(4.14)
()
( )
(4.13)
( )
p =
( )(/ )
[ ( )
= ( )
=
( )
(4.20)
(4.16)
193
Chapter 4
compressors
(4.21)
flow
(4.24)
, = Cp ( T04 T02 )
= , = Cp ( T02 T01 )
(4.23)
(4.22)
Mu =
(4.25)
194
Chapter 4
particular
question.
compressors
fluid-dynamic
problem
in
(4.30)
The two quantities defined above can be
interpreted as dimensionless work per
Re =
(4.26)
Re =
p =
(4.27)
( )
(4.29)
[ ( )
(4.31)
(4.28)
195
Chapter 4
compressors
(4.32)
= [
( + ) ]
= [ ( + ) ]
(4.33)
196
Chapter 4
machines
deviation
too
the
compressors
phenomenon
of
(4.35)
(4.36)
(4.37)
4.2.5.3 Conventional
Representation of Pressure Drop in
Compressors
(4.34)
197
Chapter 4
compressors
WA= KA
WA = KA Q002
(4.39)
(4.38)
198
Chapter 4
compressors
WT = WI+WA
=KAQ002 + KU(Q00-Q00*)2 + Ko
(4.42)
or in dimensionless form:
WT = KA + KU (1-1*)2 + Ko
(4.43)
FIGURE 4.28 Losses breakdown as a function of the flow
coefficient.
Fig. 4.28.
199
Chapter 4
compressors
200
Chapter 4
compressors
201
Chapter 4
compressors
202
Chapter 4
compressors
Q1= CQ 2
CQ =
Hp = KN2
(4.45)
Q1 =
(4.49)
Hp = KQ2
(4.48)
(4.44)
()
Q1 = K N
(4.47)
(Kg/m3)
(4.50)
(4.46)
Q1 =
(4.51)
203
Chapter 4
compressors
II. Position.
1) Pre- cooler.
2) Inter cooler.
3) After cooler.
(/ )=
Cooler function
1- Protect the equipments from
overheating due to material temperature
limitation.
Texit = Tinlet
(4.53)
4.3.2Cooler classification
Coolers classified according to:
I. The refrigerant Fluid.
1) Water-cooled Heat Exchangers.
2) Air-cooled Heat Exchangers.
204
Chapter 4
compressors
% Uncooled Hp =
(4.56)
Therefore, the minimum %Uncooled HP
could be satisfied is:
% Uncooled Hpmin =
(4.54)
(4.57)
WI.T = m Rgas T ln ( )
PV = constant
(4.55)
205
Chapter 4
compressors
(4.58)
()
(4.59)
(4.60)
206
Chapter 4
compressors
1) Electrical motors.
207
Chapter 4
compressors
2) Steam turbines.
Mechanical drive steam turbines are prime
movers for compressor, blower, and pump
applications. Steam turbines are available
for a wide range of steam conditions,
horsepower, and speeds. Typical ranges
for each design parameter are:
Inlet Pressure, psig 30 2000
Inlet Temperature, F saturated 1000
Exhaust Pressure, psig saturated 700
Horsepower 5 100,000
Speed, rpm 1800 14,000
4) Gas turbines.
The gas turbine was first widely used as an
aircraft power plant. However, as they
became more efficient and durable, they
were adapted to the industrial
marketplace. Over the years the gas
turbine has evolved into two basic types
for high power stationary applications: the
industrial or heavy-duty design and the
aircraft derivative design.
208
Chapter 4
compressors
209
Chapter 4
compressors
210
Chapter 4
compressors
compression
(in
(4.61)
211
(4.62)
Chapter 4
compressors
= 1
= 1
(4.63)
( )
( )
= 1
( )
( )
(4.64)
=( )
= ( )
th,regen = 1 ( ) ( )
(4.68)
(4.65)
th,brayton = 1
Where
(4.66)
Figure 4.44 A gas-turbine engine with regenerator.
(4.67)
CYCLE
WITH
Figure 4.45 T-S diagram of a Brayton cycle with
regeneration.
212
Chapter 4
compressors
THE
BRAYTON
CYCLE
INTERCOOLING,
REHEATING,
REGENERATION
WITH
AND
Figure 4.48 T-S diagram of an ideal gas-turbine cycle
with intercooling, reheating, and regeneration.
213
Chapter 4
compressors
2) Compressor Section.
In the compressor, air is confined to the
space between the rotor and stator blading
where it is compressed in stages by a series
of alternate rotating (rotor) and stationary
(stator) airfoil shaped blades. The rotor
blades supply the force needed to
compress the air in each stage and the
stator blades guide the air so that it enters
in the following rotor stage at the proper
angle. The compressed air exits through
214
Chapter 4
compressors
3) Diffuser.
The diffuser is located between the
compressor section and the main
components of the combustion section.
The purpose of the diffuser is to prepare t
the air for entry into the combustion
section. The front end of the diffuser is
215
Chapter 4
compressors
6) Exhaust Section.
5) Turbine Section.
2- Humidity.
3- Inlet and Exhaust Losses.
4- Fuels.
5-Diluent Injection.
6- Air Extraction.
7- Inlet Cooling.
4.5 Scrubber
4.5.1 Definition
A scrubber is a two-phase separator that is
designed to recover liquids carried over
from the gas outlets of production
separators or to catch liquids condensed
due to cooling or pressure drops. Liquid
loading is much lower than that in a
separator.
Figure 4.56 turbine section details.
216
Chapter 4
compressors
4.5.2 Applications
1) large liquid
accommodated.
slugs
have
to
be
217
Chapter 4
compressors
I. Feed inlet.
6) filter candles.
Options include:
1) wiremesh
3)
Schoepentoeter-mistmat-swirldeck
separator
4) Schoepentoeter-vane pack-swirldeck
separator
3) (multi)cyclones.
218
Chapter 4
compressors
5)
Schoepentoeter-mistmat-swirldeckmistmat separator.
6) Cyclone with straight inlet.
7) Vertical separator with reversed-flow
multicyclone
bundle
(conventional
multicyclone).
Vertical wiremesh demister is the most
commonly used in compressor station
then vane-type. Others types stated above
not used in this case because they have
high pressure drop.
Figure 4.59 Horizontal gas scrubber
219
Chapter 4
compressors
To
safeguard
the
environment,
atmospheric discharges shall be limited to
nonflammable and non-toxic gases. Any
requirement to discharge hydrocarbons to
atmosphere shall be highlighted at the
earliest opportunity .The discharge point
shall be at a safe location and the flow
directed away from areas where personnel
may be present.
220
Chapter 4
compressors
221
Chapter 4
compressors
Design Feature
3- Liquid seal
Process vent streams are usually passed
through a liquid seal before going to the
flare stack. The liquid seal can be
downstream of the knockout drum or
incorporated into the same vessel. This
prevents possible flame flashbacks, caused
when air is inadvertently introduced into
222
Chapter 4
compressors
4- Flare stack
6- Burner tip
7- Pilot burners
223
Chapter 4
compressors
hydrocarbon
molecule,
thereby
minimizing polymerization, and forms
oxygen compounds that burn at a reduced
rate and temperature not conducive to
cracking and polymerization. Another
theory claims that water vapor reacts with
the carbon particles to form CO, CO2, and
H2, thereby removing the carbon before it
cools and forms smoke. An additional
effect of the steam is to reduce the
temperature in the core of the flame and
suppress thermal cracking.
The physical limitation on the quantity of
steam that can be delivered and injected
into the flare flame determines the
smokeless capacity of the flare. Smokeless
capacity refers to the volume of gas that
can be combusted in a flare without smoke
generation. The smokeless capacity is
usually less than the stable flame capacity
of the burner tip. Significant disadvantages
of steam usage are the increased noise and
cost. Steam aggravates the flare noise
problem by producing high-frequency jet
noise. The jet noise can be reduced by the
use of small multiple steam jets and, if
necessary, by acoustical shrouding. Steam
injection is usually controlled manually
with the operator observing the flare and
adding steam as required to maintain
smokeless operation. To optimize steam
usage infrared sensors are available that
sense flare flame characteristics and adjust
the steam flow rate automatically to
maintain smokeless operation. Automatic
control, based on flare gas flow and flame
radiation, gives a faster response to the
need for steam and a better adjustment of
the quantity required. If a manual system
is used, steam metering should be installed
to
significantly
increase
operator
awareness
and
reduce
steam
consumption.
8- Steam jets
A diffusion flame receives its combustion
oxygen by diffusion of air into the flame
from the surrounding atmosphere. The
high volume of fuel flow in a flare may
require more combustion air at a faster
rate than simple gas diffusion can supply.
High velocity steam injection nozzles,
positioned around the outer perimeter of
the flare tip, increase gas turbulence in the
flame boundary zones, drawing in more
combustion air and improving combustion
efficiency. For the larger flares, steam can
also be injected concentrically into the
flare tip.
The injection of steam into a flare flame
can produce other results in addition to air
entrainment and turbulence. Briefly, one
theory suggests that steam separates the
224
Chapter 4
compressors
9- Controls
4.6.3 Operation
1) Video camera
2) Pilot gas flow and temperature Relief
gas flow, pressure, temperature and
composition
225
Chapter 4
compressors
! WARNING
226
Chapter 4
compressors
227
Chapter 4
compressors
228
Chapter 4
compressors
Figure 4.68 first operating scenario (two booster compressor and two main compressor)
229
Chapter 4
compressors
Second
scenario:
COMPRESSOR OFF)
(BOOSTER
Figure 4.69 second operating scenario (booster compressor off and three main compressors on)
Fig. (4.70) third operating scenario (booster compressor on and main compressor off)
230
Chapter 4
compressors
Fig. (4.71) fourth operating scenario (booster compressor off and three main compressors on and two main compressor
stand by)
231
Chapter 4
compressors
C2/C3.
Commercial propane(exporting)
Condensate (C5+).
232
Chapter 4
compressors
233
Chapter 4
compressors
234
Chapter 4
compressors
235
Chapter 4
compressors
236
Chapter 4
compressors
237
Chapter 4
compressors
238
Chapter 4
compressors
4.7.3.4Compressor suction
description (k.o drum)
2. Compressor 200-K-11x.
239
Chapter 4
compressors
1. Compressor 200-K-11x.
2. Scrubber 200-V-81.
3. Filters 200-FL-81 A/B.
240
Chapter 4
compressors
241
Chapter 4
compressors
1. inlet SDV-x915.
242
Chapter 5
Modeling Using ArcGIS
243
5.1 introduction:
We have two reasons for drawing the Gas
network model on ArcGIS, Due to the
abilities ArcGIS provides to draw accurate
map of Egyptian national grid of natural
gas we use this map as a model on
SynerGEE package to perform hydraulic
analysis and optimize its operation.
After
performing
analysis
and
optimization, generated data from
SynerGEE can be added to the attribute
table of each layer this could make us
supported with accurate map with
accurate coordinates of each feature and
hydraulic data of each point on the map to
be displayed.
244
5.2.2.4 People:
GIS technology is of limited value without
the people who manage the system and
develop plans for applying it to real world
problems. GIS users range from technical
specialists who design and maintain the
system to those who use it to help them
perform their everyday work.
5.2.2.5 Procedures or methods:
5.2.2.1 Hardware:
Hardware is the computer system on
which a GIS operates. Today, GIS software
runs on a wide range of hardware types,
245
For example, you can use GIS to map pipecorrosion patterns, allowing you to predict
risks and meet DIMP requirements. The
technology lets you view patterns and
trends to predict the future and, more
importantly, prepare for it.
246
essential:
One step beyond data capture is data
editing and maintenance. Editing involves
functions such as adding, duplicating and
deleting spatial and attribute data as well
as creation and management of metadata.
Using tools to detect errors and to perform
transformations is another aspect of data
management.
247
248
249
Measurements,
computations
and
computer applications are more difficult to
manage using angular units, thus, latitude
and longitude values are commonly
converted to a rectangular coordinate
system of northings (Y-axis) and
eastings (X-axis) that can be expressed
and easily understood in linear units such
as meters or feet.
250
Figure5.8
Projection.
251
252
253
5.5.2 ArcGlobe:
ArcGlobe is used for seamless 3D
visualization of geographic data using a
continuous global view.
5.5 Mapping:
5.5.3 ArcScene:
ArcScene is used for 3D visualization of
focused scenes or areas.
5.5.1 ArcMap:
254
ArcMap
represents
geographic
information as a collection of layers and
other elements in a map.
Common map elements include the data
frame containing map layers for a given
extent plus a scale bar, north arrow, title,
descriptive text, a symbol legend, and so
on.
Figure5.12 ArcScene visualization.
5.5.4 ArcCatalog:
255
Data view
Layout view
Each view lets you look at and interact with
the map in a specific way.
256
257
258
259
260
261
262
263
264
Chapter 6
Chapter 6
Modeling Using Synergi
265
Chapter 6
6.1 introduction.
6.1.1 about synergi
Synergi Gas is the leading natural gas
network modeling software and it is the
commercial package that is used by GASCO
& NATA for modeling the national grid of
natural gas.
1. Solved
266
Chapter 6
2. Unsolved
3. Infeasible
If an unknown variable has reached a
limiting
value
and
constraint
interchange does not apply, an
infeasible solution is reached. An
infeasible solution is not a broken
simulation or failure; often, an
infeasible solution contains results that
are meaningful. Possible infeasible
solutions are pressures below absolute
zero, or regulators closing due to sonic
conditions.
Node flow.
Node pressure.
Node pressure of an adjacent
connected node or a facility
parameter of a facility connected to
this node as a control node
267
Chapter 6
6.1.6.1.1 Introduction
Optimizer simulates elements of a gas
pipeline using Synergis modeling
techniques.
268
Chapter 6
Figure 6- 2
Figure 6- 1
269
Chapter 6
Figure 6- 3
minimization
During
a
steady-state
analysis,
assumptions must be made about the
systems operating parameters. Usually,
either the suction or discharge pressure of
a compressor is set. This pressure is
maintained for as long as the compressor
engine has the horsepower to do so When
you run an optimization analysis, the
following takes place:
solution
technique
for
fuel
minimization
Performing a fuel minimization analysis
270
Chapter 6
Figure 6- 4
create an equation
edit an equation
rename an equation
delete an equation
Figure 6- 5
Figure 6- 6
271
Chapter 6
- Total Deliveries
272
Chapter 6
273
Chapter 6
274
Chapter 6
increased
system
reliability
and
significantly reduced operating cost.
Transient analysis helps you to minimize
compressor fuel usage, determine proper
emergency operating procedures, and
assess system survival time under certain
conditions.
Gas blow down or lost gas simulations. If a
pipeline must be blown down, the
engineer can determine the safest and
fastest way to perform this operation. USM
is an excellent tool for designing pipeline
venting systems and determining proper
evacuation times. In addition, USM is often
used to perform lost gas calculations in
case of a pipeline rupture. How unsteadystate modeling complements steady-state
modeling
275
Chapter 6
=
Transient analyses, in contrast, are not
constrained by the Q relationship. This
category uses completely different
numeric algorithms and techniques to
solve differential equations that relate
mass flow rate, momentum, velocity,
acceleration, and friction, etc., with
respect to time and horizontal pipeline
distance.
276
Chapter 6
6.1.7 script
Synergies scripts are a very powerful tool
that allows you to automate tasks that you
want to repeatedly perform in your model.
Synergi provides for several different types
of scripts, each with a distinct purpose:
277
Chapter 6
Figure 6- 8
Figure 6- 9
278
Chapter 6
Figure 6- 12
Figure 6- 10
Figure 6- 13
Figure 6- 11
279
Chapter 6
Figure 6- 17
Figure 6- 15
Figure 6- 16
Figure 6- 19
280
Chapter 6
geodatabases,
SDE
(multi-user)
geodatabases, and query layers. It is not
possible to import both CSV files and nonCSV (spatial DSD) files at the same time.
For CSV import, the entire set of CSV files
in a specified directory will be converted.
The types of CSV files in the specified
directory control the set of imported
facilities. Each set of CSV files must have at
least a node file and a pipe file (which is the
same requirement in MiddleLink 3.x).
For spatial DSD import, you can select DSDs
that represent combinations of conversion
options.
You can:
Create Synergi pipe facilities, non-pipe
facilities, and isolation valves.
Create only Synergi pipe facilities.
Create only Synergi non-pipe facilities.
You can also define a DSD to create Synergi
nodes. Such a DSD only has meaning in the
context of also importing one of the above
combinations of pipe and non-pipe
facilities. Any nodes remaining at the end
of the conversion process that do not have
at least one attached facility are
automatically deleted from the model. An
alternative form of a node DSD is a node
insert DSD. When imported as a node
insert DSD, Model Builder can attempt to
locate a pipe at the node points XY
coordinate, split the pipe, and insert the
node, and also to update existing nodes
with data from incoming nodes.
281
Chapter 6
6.2.3.3 Attributes
Source file attributes can be mapped to
Synergi facility attributes, including
dynamic attributes (such as multiple flow
categories) and custom attributes. These
attribute field mappings, along with other
user selections, can be saved in a
configuration file so they can be reused in
future
Model Builder conversions.
NOTE: Custom attributes must exist in the
object model before they can be mapped
properly. If you have an existing model
with the necessary custom attributes, you
may find it easier to import them using
Model Builders supplemental data import
feature. These custom attributes must
have been assigned to model features.
Figure 6- 20
282
Chapter 6
283
Chapter 6
2. Click
.
3. In the DSD list, type the new name for
the DSD and then press Enter.
To delete a DSD
1. From the Model Builder Configuration
dialog box, under the list of DSDs (located
on the right side of the dialog box), select
the DSD that you want to delete.
2. Click
3. Click
DSD.
284
Chapter 6
285
Chapter 6
286
Chapter 6
Figure 6- 21
Figure 6- 22
287
Chapter 6
NAME
Flow rate
(m3/h)
WDGC
934 500
Abu Qir
268 500
Badre El-Din
630 000
Rashpetco
1 085 000
UDGC
1 657 000
Abo Mady
630 000
El.Mansoura
126 000
Dahshour
123 000
Ras bakr
25000
Ras Shouqir
25 000
North Sinai
90 000
288
Chapter 6
Delta
Alexandria
REGION
NAME
FLOW (m3/hr)
TYPE
Abo Qir
87 000
Steam
El Mahmoudia
68 000
Combined
Damanhour
46 000
Combined
El Seoof
50 000
Gas
Kafr El Dawar
74 000
Steam
42 000
steam
90 000
Combined
Intergen (3&4)
138 000
Steam
Atf new
138 000
Combined
51 000
Steam
Damietta
208 000
Combined
Talkha gas
64 000
Gas
Talkha steam
51 000
Steam
Nubaria
380 000
Combined
Talkha Comb.
125 000
Combined
New damietta
138 000
Gas
West damietta
140 000
Gas
Banha
95 000
Gas
289
Cairo
Chapter 6
Shobra El Khema
33 000
Steam
South Cairo
56 000
Combined
West Cairo
63 000
Steam
North Cairo
168 000
Combined
Wady hoof
11 000
Gas
Teppin
14 000
Steam
Korymat
127 000
Steam
Korymat (comb.)
123 000
Combined
Abou Rawash
150 000
Gas
Ataka
46 000
Steam
Abou Sultan
28 000
Steam
Port Said
14 000
Gas
El Shabab
30 000
Gas
Mousa
156 000
Steam
West Suez
153 000
Steam
180 000
steam
El Arish
19 000
Steam
The data shown in the table are collected from confirmed data sheets.
290
Chapter 6
Alexandria
REGION
INDUSTRY
Petroleum ,
Derivatives
And
petrochemicals
NAME
FLOW RATE(m3/h)
Alex petroleum
19 700
Ameriya petroleum
26 600
petrochemical
27 000
64 500
Amoic
10 500
Sidpic
7 000
Poly probelin
21 800
Misr Chemical
6 100
7 500
46 000
Misr petroleum
160
ANRPC
6 770
5000
Star Nex
2000
129 800
Alexandria Fertilizer
40 000
Dekheila Steel
80 000
Cement
MAX Cement
17 550
Fertilizers
291
Chapter 6
Others
Petroleum ,
Derivatives
And
petrochemicals
Fertilizers
2 720
1 870
Alex IA
18 420
Alex sugar
31200
Nobaria Sugar
18030
Methanics
23 000
Fenosa Facilities
4 000
Uoria
72 400
Natarat
Delta
Ghazl el mahla
Others
Fertilizers
Cairo
America Cement
Cement
Others
200
19 500
Quesna(Egypt gas)
2700
Repco GAS
2800
Gastec Talkha
800
Dakhlia Sugar
26 000
4 400
Cebaght Elmahala
2800
Helwan Fertilizer
50 000
Helwan Cement
20 000
National Cement
37 900
Ktamya cement
7 400
Kama en El Tob
75 000
6 000
35000
10th of Ram
95 000
292
Chapter 6
Petroleum ,
Derivatives
And
petrochemicals
Fertilizers
Cement
Hawmdya sugar
13 800
Fayoum Gas
22300
14000
EL Nasr Petroleum
7 500
Suez Petroleum
11 500
Gas Derivatives
40 000
Morgan
31 000
Gapco
7 000
42 000
Suez Fertilizer
15 000
Suez Cement
1 000
Arabiya Cement
1 000
Egypt Cement
5 000
Sweedy Cement
1 000
Gastec Suez
Others
400
Sinai IA
15000
Epec
45000
5000
Upper Egypt IA
5000
Abo Qurqas
6000
El Minia IA
1000
Qift
17700
5000
El Minia
10000
Assuit
10000
293
Chapter 6
Upper
Egypt
Industries
Qena
22800
Naga hamadi
11500
Kima
185 800
Aswan
30000
Bani swaif
5000
10000
The data shown in the table are collected from confirmed data sheet
6.2.4.3.2.3 Residential sector
The residential category is the less Priority as it represents 3-5 % from the total gas
consumption.
table 6- 4 residential flowrates
sector
Residential
Sector
NAME
14400
9920
34620
5600
800
9000
4420
2410
2810
1970
15000
Trans Gas-Delta
13000
Egypt Gas-Delta
40000
15400
95000
294
Chapter 6
40000
15000
72000
35000
5000
16000
City Gas
124000
The data shown in the table are collected from confirmed data sheets.
295
Chapter 6
Figure 6- 23
296
Chapter 6
Synergi facilities:
Pipes
Storage facilities
Processing plants
Loss elements
Forced flow elements
Valves
Regulators
Regulator stations
Compressor stations
297
Chapter 6
6.3.3 Pipes
Pipes are the most basic component of any
system, connecting all of the facilities in
the system into a functional network. In
Synergi, you can select from different pipe
equations to represent the different
characteristics of each pipe in the system.
Different pipe types can be mixed together
in a single modeled system.
298
Chapter 6
Pipe diameter
Pipe length
Pipe roughness
pipe equations
6.3.4 valves
A valve is a facility that has a fixed valve
stem position with varying inlet and outlet
pressures. A valve is not a pressure
regulating or sustaining device.
299
General valve
Check valve
Fisher valve.
Mokveld valve.
Reliance valve.
Valtek valve
Chapter 6
6.3.5
Regulators
regulating stations
and
300
Chapter 6
Figure 6- 28
6.3.5.3 REGULATORS
Adding a regulator: To add a regulator to
your model, select a regulator in your
equipment warehouse, drag it onto the
map, and click the points where you want
the from-node and the to-node to be
located. For detailed information on
adding facilities
Editing a regulator: (Service state Reference flow-Minimum pressure DropInlet temperature-Travel time-Regulator
equation-Pressure
Profile-Maximum
flowrate)
regulator stations
Regulator stations provide a number of
advantages over modeling individual
regulators directly in the Synergi model.
301
Chapter 6
302
Chapter 6
From
Warehouse>Regulators>Sample
(fisher, General, Reliance) regulator
To
add
regulator
station
Warehouse>Regulator
station>Sample
regulator station
Setting pressure is the same as regulators
but to edit the entire station click on
Figure 6- 32
303
Chapter 6
304
Chapter 6
gas (
(
)
$
gas (
) = (
Figure 6- 35
) = (
( )
electric (
(
) = 24 (
) ()
In
ambiguous
situations,
Synergi
determines which path or series of paths
should be operating. For example, when
305
Chapter 6
6.3.6.6 Drivers
Drivers are one of two required
components for a compressor station unit.
(The other required component is the
compressor itself).
Editing drivers: TO EDIT A DRIVER
Theoretical
power/flow
(TH)
compressor
Actual power/flow theoretical form
(HP) compressor
Actual ratio/flow polynomial form (C2)
compressor
N-stage
reciprocating
(RN)
compressors
3.0303
[( )
306
1]
Chapter 6
= [1 ( )3 2 ]( )
ec = compressor
(dimensionless)
efficiency
Figure 6- 36
[( )
= 3.0303
1
1]
Ec : Compression efficiency
Em : Mechanical efficiency
= (0 + 1 + 2 2 + 3 3 + 4 4 )
( ) 1
307
Chapter 6
6.3.6.8 Coolers
Coolers are optional and are not required
to perform an analysis. Synergi offers two
types of coolers for compressor station
modeling: generic coolers (the traditional
cooling that has been available in Synergi)
and detailed aerial coolers. Detailed aerial
coolers are only available in Enhanced
compressor station mode when you have
selected either Aerial Cooler Only or Mixed
Cooler as the Cooler Preference in the
compression controls.
1. Cooling
method
(Delta
temperature & Target temperature
& Delta heat removed)
2. Heat removal rate or Degrees of
cooling
3. 3. Use Suction Temperature : This
option is available if you selected
Target Temperature as the Cooling
method
4. 4.
Maximum
Temperature:
Specifies the limiting output
temperature. This option is
available if you selected Target
Temperature as the Cooling
method. The temperature of the
gas discharged from an individual
compressor is given by:
= +
1
( ) 1
= [
]
Generic coolers
Ps = Pin S drop
-Flow dependent pressure drop
308
Chapter 6
=
= 2 +
2
(. 000024 )2
= ( , 11.2 )
7. Wall thickness
8. Outside diameter
Use the Design Conditions tab to specify
rated conditions that help determine how
the aerial cooler will operate under
specified conditions.
= .000024 2 2
2. Pressure drop
3. Enable constraints
4. Minimum loss
5. Maximum loss
1.
2.
3.
4.
5.
6.
7.
6. Diameter
2. Thermal conductivity
3. Inside Heat Transfer coefficient
4 .Outside Heat Transfer coefficient
Click OK
Aerial coolers
1. Method
309
Chapter 6
2. Valve constant
Figure 6- 37
Recycle Lines
Figure 6- 38
Figure 6- 39
310
Chapter 6
Figure 6- 40
Figure 6- 41
6.3.6.16 Digitizer
How to define Compressors performance
curve?
The default compressor station is using a
theoretical compressor unit
The required compressor unit is centrifugal
not theoretical
311
Chapter 6
Figure 6- 42
6.3.6.16.1 Terminologies
Figure 6- 43
312
Chapter 6
Figure 6- 45
Figure 6- 47
Figure 6- 46
Figure 6- 48
313
Chapter 6
6.3.6.16.3 Data
To get the points of variables
Figure 6- 51
6.3.6.16.5 Importing
Figure 6- 49
6.3.6.16.4 Exporting
Export CSV file for the performance curve
and the data obtained
Figure 6- 52
6.4.1 Reports
Reporting provides detailed information
about your model in a format that is easy
to read and use. Reports are generated in
either an HTML or a table-based format
and can be saved to different formats,
based on the type of report. In addition to
the different report types that can be
generated in the Report window, Synergi
can also create certain analysis based
reports, which are saved as text files in
your personal Synergi directory.
Figure 6- 50
314
Chapter 6
Compressor report
Demand Contribution report
Dew Point report
Exception report
Full report
Hydrate report
Length and Diameter report
Model report
Regulator report
Storage Field report
315
Chapter 6
Nodes
Pipes
Looped pipes
Valves
Regulators
Compressors
Storage fields
Forced flow elements
Loss elements
Processing plants
Demand
Supply
Fuel and Power
Status Summary
316
Chapter 6
6.4.2 Charts
Chart Viewer allows you to view different
charts relating to a Synergi model. Charts
can be created and viewed using the Chart
tab within Synergi and also using the
standalone Chart Viewer application.
Distance charts
Time charts
Centrifugal compressor maps
Distance charts
317
Chapter 6
318
Chapter 6
319
Chapter 6
320
Chapter 6
321
Chapter 6
322
Chapter 6
Name
Result pressure
Ras Shouqir
37.3
North Sinai
41.1
WDGC
46.9
Dahshour
27.5
Ras bakr
37.5
LPG(badreldin)
45.1
Aboqir
35.1
Rashpetco
39.5
Abo madi
35.4
UGDC(port fouad)
44.6
Elmansoura
25.5
Ras shouqir
37.3
North sinai
41.1
323
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
ps1
talkha steam
-51000.0
9.00
ps2
shobra el khema
steam
krir 1,2 steam
-33000.0
9.00
-42000.0
9.00
-138000.0
9.00
ps5
intergen 3,4
steam
abo qir steam
-87000.0
9.00
ps6
-51000.0
9.00
ps7
-74000.0
9.00
ps8
kafr el dawar
steam
west cario steam
-63000.0
9.00
ps9
teppin steam
-14000.0
9.00
ps10
korymat steam
-127000.0
9.00
ps11
-180000.0
9.00
ps12
-156000.0
9.00
ps13
ataka steam
-46000.0
9.00
ps14
-153000.0
9.00
ps15
al arish steam
-19000.0
9.00
ps16
-28000.0
9.00
psg1
el attf new
combined
damietta
combined
new damietta gas
-138000.0
27.00
-208000.0
18.49
-138000.0
18.62
ps3
ps4
psg2
psg3
324
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
psg4
-140000.0
27.00
psg5
talkha combined
-125000.0
23.32
psg6
talkha gas
-64000.0
23.70
psg7
nubaria combined
-380000.0
17.14
psg8
banha gas
-95000.0
23.31
psg9
-90000.0
27.00
psg10
noth cairo
combined
-168000.0
22.39
psg11
6 october gas0
-150000.0
21.79
psg12
south cairo
combined
-56000.0
23.31
psg13
-11000.0
23.44
psg14
korymat comb
-123000.0
27.00
psg15
elshabab gas
-30000.0
27.00
psg16
-14000.0
27.00
psg17
elsyof gas
-50000.0
27.00
-46000.0
27.00
-68000.0
27.00
psg18
psg19
damanhor
combined
el mahmoudia
combined
325
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
in1
natgas borg el
arab
-34620.0
7.00
in2
natgas agamy
-9920.0
7.00
in3
-14400.0
7.00
in4
-15000.0
35.26
in5
natgas kafr
eldawar local
-5600.0
4.00
in6
-16000.0
42.76
in7
transgas delta
-13000.0
25.23
in8
-40000.0
22.54
in9
quesna(egyptgas)
-2700.0
23.60
in10
egyptgas delta(2)
-40000.0
21.46
in11
egyptgas abo
zabal
-15400.0
23.62
in12
egyptgas cairo 1
-45000.0
23.11
in13
natgas cario
-15000.0
23.93
in14
towngas cairo
-72000.0
21.39
in15
natgas 6 october
-35000.0
25.99
in16
egyptgas cairo 2
-45000.0
32.25
in17
towngas el
ismailia
-5000.0
38.78
in18
city gas
-124000.0
28.03
in19
midor and
medalic
-64500.0
38.69
in20
star nex
-2000.0
7.00
in21
misr chemical
-6100.0
7.00
in22
alex petrolum
-19700.0
7.00
in23
amoic
-10500.0
42.92
in24
anrpc
-6770.0
43.12
in25
polyprobelin
-21800.0
43.37
326
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
in26
elmogammaa for
petrol
-46000.0
40.49
in27
amryia petroleum
-26600.0
43.17
in28
misr petroluem
-160.0
44.76
in29
sidpec
-7000.0
46.70
in30
petrochemical
-27000.0
45.46
in31
carbon black
factory
-7500.0
46.45
in32
fenosa facilities
-4000.0
40.11
in33
methanics
-23000.0
39.90
in34
morgan
-31000.0
36.54
in35
gapco
-7000.0
36.85
in36
suez petroluem
-11500.0
28.17
in37
elnasr petroluem
-7500.0
26.86
in38
gas derivatives
-40000.0
26.55
in39
-129800.0
33.61
in40
alexandria fert
-40000.0
31.26
in41
-72600.0
24.28
in42
helwan fert
-50000.0
24.53
in43
-42000.0
33.37
in44
-15000.0
33.28
in45
amrya cement
-2720.0
42.14
in46
dkhela steel
-80000.0
7.00
327
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
in47
max cement
-17550.0
7.00
in48
ganob elwadi
cement
-6000.0
25.10
in49
katmya cement
-7400.0
23.21
in50
helwan cement
-20000.0
23.47
in51
national cement
-37900.0
23.42
in52
kamaen el top
-75000.0
22.31
in53
-35000.0
24.32
in54
swedy cement
-1000.0
34.12
in55
suez cement
-1000.0
34.12
in56
arabya cement
-1000.0
33.70
in57
egypt cement
-5000.0
33.62
in58
-1870.0
42.21
in59
alex sugar
-31200.0
7.00
in60
nubaria sugar
-18030.0
7.00
in61
alex IA
-18420.0
32.20
in62
dakhlia sugar
-26000.0
39.14
in63
repco gas
-2800.0
24.79
in64
damietta pump
station
-4400.0
40.24
in65
gastec talkha
-800.0
25.43
in66
gazl el mahala
-19500.0
25.00
in67
cebaghet el
mahla
-2800.0
25.24
328
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
in68
10 of ram
-95000.0
24.64
in69
hawamdya sugar
-13800.0
22.99
in71
fayoum gas
oftake
-22300.0
27.76
in72
-14000.0
41.93
in74
epec
-45000.0
33.23
in75
gastec suez
-400.0
27.35
-15000.0
41.08
-5000.0
41.57
in76
in78
in79
bani swaif
-5000.0
41.12
in80
upper egypt IA
-5000.0
41.10
in81
-5000.0
41.07
in82
elmania
-10000.0
37.22
in83
elmania AI
-1000.0
36.94
in84
abo qurqas
-6000.0
36.62
in85
asyout
-10000.0
33.08
in86
naga hamadi
-11500.0
27.37
in87
qena
-22800.0
24.44
in88
qift
-17700.0
21.13
in89
aswan
-30000.0
14.24
in90
egyptgas aswan
-10000.0
12.80
in91
kema fert
-185000.0
10.00
329
Chapter 6
No
Name
Flow
SM^3/h
Result pressure
in101
-800.0
4.00
in102
local damanhur
-9000.0
4.00
in103
local
elmahmoudia
-4420.0
6.00
in104
-2410.0
7.00
in105
local rashid
-2810.0
7.00
-1970.0
7.00
-5000.0
7.00
local wadi
elnatron
sodum
carbonat(solvay)
in106
in107
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi1
Pi2
Pi3
Pi5
Pi4
Pi7
Pi6
Pi8
Pi9
Pi10
Pi12
Pi13
Pi15
Pi14
Pi16
Pi17
Pi18
Pi19
Pi20
24
24
24
24
24
24
16
8
24
8
24
24
24
8
24
10
10
24
6
14
18
13
1
1
1
1
2
2
1
3
1
9
1
11
1
1
3
1
-2720
37340
270000
138000
90000
228000
42000
64500
101840
1870
113768
373710
103848
-9920
-269862
-2000
-80000
195768
-27000
in45
No3
No5
No7
No7
No9
No9
No3
No12
No12
No14
No15
No17
No19
No15
No29
No23
No22
in30
42.14
42.18
42.22
40.85
40.85
40.93
40.93
42.18
42.21
42.21
42.6
42.41
42.54
42.5
42.41
42.65
41.79
42.8
45.46
No2
No2
No9
No329
No348
No7
No328
in19
No3
in58
No17
No5
No15
No17
No18
No20
No14
No20
No24
42.14
42.14
40.93
40.83
40.84
40.85
40.91
38.69
42.18
42.21
42.54
42.22
42.41
42.54
43.47
42.65
42.6
42.65
46.66
330
Result
Downstre
am
Velocity
(m/s)
0.1
0.8
5.6
2.9
1.9
4.8
2
12.9
2.1
0.3
2.3
7.6
2.1
1.8
5.4
0.2
9.4
3.9
8.1
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi21
Pi22
Pi23
Pi24
Pi25
Pi26
Pi27
Pi29
Pi28
Pi30
Pi31
Pi32
Pi33
Pi34
Pi35
Pi37
Pi36
Pi39
Pi38
Pi40
Pi41
Pi42
Pi43
Pi45
Pi44
Pi47
Pi46
Pi49
Pi50
Pi51
Pi52
Pi54
Pi53
Pi55
Pi56
Pi58
Pi57
Pi59
6
18
32
8
6
24
24
32
6
6
6
6
24
24
24
18
6
18
6
18
24
6
24
18
6
32
10
18
16
8
18
18
8
16
8
8
4
28
1
1
1
1
1
1
2
17
1
1
1
1
3
1
1
12
1
2
1
10
2
1
2
1
1
39
1
10
5
9
3
4
2
4
4
9
1
8
-7000
-108304
690227
-6100
-17550
224418
206868
648727
-7500
-19700
-10500
-6770
261388
254618
244118
-181064
-160
-180904
-26600
248463
630000
21800
-531250
-154304
-46000
617527
-31200
234063
50000
-15000
184063
145043
-18420
-129800
-40000
-24020
-5600
268500
in29
d1
s14
No25
No30
No33
No32
No35
in31
No27
in23
in24
No18
No40
No38
No42
in28
No44
in27
d1
s2
d1
No18
No49
in26
No51
No28
No53
No45
in4
No45
No59
in61
in39
in40
No64
No13
s3
46.7
44.19
46.91
42.79
42.42
43.04
42.97
46.54
46.45
42.34
42.92
43.12
43.47
43.2
43.12
44.76
44.76
44.4
43.17
44.19
45.19
44.19
43.47
44.26
40.49
44.55
44.44
40.54
36.97
35.26
36.97
35.24
32.2
33.61
31.26
32.57
31.94
35.02
No26
No49
No26
No22
No32
No32
No41
No51
No35
No33
No38
No40
No40
No38
No33
s14
No42
No42
No44
No53
d1
in25
d1
No44
No49
No75
No51
No45
No316
No56
No56
No57
No64
No61
No57
No59
No64
d2
46.78
44.26
46.78
42.8
42.97
42.97
42.86
44.55
46.54
43.04
43.12
43.2
43.2
43.12
43.04
46.91
44.76
44.76
44.4
40.54
44.19
43.37
44.19
44.4
44.26
40.08
44.55
36.97
36.8
36.26
36.26
34.62
32.57
34.54
34.62
35.24
32.57
34.59
331
Result
Downstre
am
Velocity
(m/s)
2
3.7
7
1.1
5.7
4.5
4.1
7
2.2
6.4
3.3
2.1
5.2
5
4.8
6.1
0
6.2
8.4
9.3
12.1
6.9
10.4
5.3
15.6
7.4
3.4
9.7
2.6
3.3
7.8
6.4
4.5
7.5
10
5.8
5.5
4.9
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi60
Pi61
Pi63
Pi62
Pi64
Pi67
Pi66
Pi68
Pi70
Pi69
Pi72
Pi71
Pi73
Pi74
Pi75
Pi76
Pi77
Pi79
Pi78
Pi81
Pi80
Pi82
Pi83
Pi84
Pi85
Pi86
Pi87
Pi89
Pi88
Pi90
Pi92
Pi91
Pi94
Pi93
Pi96
Pi95
Pi97
Pi99
18
18
18
16
20
16
16
28
32
10
28
16
28
30
30
42
32
28
24
32
10
24
24
28
32
32
32
32
10
32
32
10
32
24
32
6
36
24
1
3
2
1
18
1
1
17
25
5.5
14
2
6.8
26
30
165
9.4
5.4
3
59
9
2
18
28.2
40
9
110
10
7
79
23
3
6
2
38
60
12
9
-105043
-24757
-48478
-51000
114279
-138000
-87000
220022
599497
-18030
146022
-74000
-516542
677971
362572
614223
-69823
-571542
-46000
552527
-45000
-68000
210420
306122
-380000
294030
304030
293030
-1000
287030
287030
-6000
333854
14210
35492
0
298361
-196374
No61
No61
No68
No331
s4
No72
No330
d2
No75
No31
No77
No332
d6
s5
s16
d7
s5
No84
No317
No86
in16
No318
d6
d6
No346
in82
c2
No97
in83
No94
No100
in84
in71
d8
No103
No103
No103
No106
34.54
34.54
34.56
34.25
35.55
34.29
34.18
34.59
40.08
39.83
33.96
33.81
29.6
39.57
35.79
33.73
39.57
32.1
32.08
35.2
32.25
27.87
29.6
29.6
19.97
37.22
40.62
36.94
36.94
35.96
36.64
36.62
27.76
25.03
27.44
27.44
27.44
27.02
No57
No68
d2
No72
No70
No70
No72
No77
No383
No75
d7
No77
No372
d7
d7
c1
No366
No368
No84
in71
No86
No88
No374
No376
No92
No97
in82
No100
No97
No98
No94
No100
No103
No378
d9
in70
c1
in71
34.62
34.56
34.59
34.29
34.55
34.55
34.29
33.96
37.09
40.08
33.73
33.96
31.14
33.73
33.73
27.16
39.59
33.48
32.1
27.76
35.2
27.89
27.99
27.16
23.48
36.94
37.22
36.64
36.94
33.49
35.96
36.64
27.44
25.03
27.42
27.44
27.16
27.76
332
Result
Downstre
am
Velocity
(m/s)
4.7
1.1
2.2
2.9
4.1
7.8
5
4.1
7.8
2.2
2.8
4.3
11.2
11.1
6
6.5
0.9
11.4
1.2
9.8
7
2.1
6.6
7.2
9.4
3.9
4
3.9
0.1
4.2
3.9
0.8
6
0.5
0.6
0
4.3
6.4
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi138
Pi140
Pi139
Pi142
Pi141
Pi143
Pi144
Pi145
Pi146
Pi149
Pi148
Pi151
Pi150
Pi153
Pi152
Pi155
Pi154
Pi157
Pi156
Pi159
Pi158
Pi161
Pi160
Pi163
Pi162
Pi165
Pi164
Pi167
Pi166
Pi169
Pi168
Pi171
Pi170
Pi173
Pi172
Pi175
Pi174
Pi177
18
30
24
28
10
32
16
16
24
28
24
24
10
24
24
24
12
24
16
24
24
24
10
24
10
24
10
24
24
24
10
24
12
24
12
16
8
28
159
3
1
28
6
6
14
10
10
6
3
2
10
4
2
3
3
6
2
8
4
2
1
1
2
4
4
5
3
3
2
11
4
2
4
3
8
18
68360
338030
-123000
-195488
40000
811608
4579
33000
238158
-155488
-95000
-68042
-13800
25158
168000
193158
45000
-366942
-14000
-135042
67000
-352942
-50000
-302942
-35000
-155042
-20000
56000
11000
-192942
-37900
-46842
-72000
-267942
-75000
-10821
15400
93800
No105
No146
No340
No147
No147
d11
d8
d8
d8
No154
No347
No156
in69
No157
No157
No159
No159
No162
No334
No163
No163
No166
in42
No168
in53
No170
in50
No171
No171
No174
in51
No176
in14
No178
in52
No179
No179
No182
42.2
42.11
42.09
23.8
23.8
41.24
25.03
25.03
25.03
23.64
23.57
23.46
22.99
23.39
23.39
23.66
23.66
25.66
25.65
23.49
23.49
25.09
24.53
24.87
24.32
23.85
23.47
23.44
23.44
24.09
23.42
23.37
21.39
24.36
22.31
25.02
25.02
25.43
No272
No143
No146
No377
in10
No227
No179
No327
No159
No147
No154
No163
No156
No371
No336
No157
in12
d9
No162
No170
No171
No162
No166
No166
No168
No174
No170
No338
No339
No178
No174
No183
No176
d10
No178
No367
in11
No136
37.34
41.96
42.11
24.97
21.46
40.01
25.02
24.81
23.66
23.8
23.64
23.49
23.46
23.38
23.24
23.39
23.11
27.42
25.66
23.85
23.44
25.66
25.09
25.09
24.87
24.09
23.85
23.4
23.44
24.36
24.09
23.43
23.37
24.7
24.36
25.03
23.62
25.26
333
Result
Downstre
am
Velocity
(m/s)
2.8
4.4
2.5
5.3
9.4
9.8
0.4
2.6
8.8
4.2
3.5
2.5
3
0.9
6.3
7.2
6.8
12.5
1.1
5
2.5
12.3
10.3
10.7
7.3
5.7
4.3
2.1
0.4
7
8.2
1.8
11.8
9.6
11.8
0.9
5.1
2.4
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
Pi176
Pi179
Pi178
Pi181
Pi180
Pi182
Pi184
Pi183
Pi185
Pi187
Pi186
Pi189
Pi188
Pi191
Pi190
Pi192
Pi193
Pi194
Pi195
Pi197
Pi196
Pi199
Pi201
Pi200
Pi203
Pi202
Pi204
Pi206
Pi205
Pi208
Pi207
Pi210
Pi209
Pi212
Pi211
Pi213
Pi214
Pi215
8
24
10
28
28
24
28
10
22
22
22
22
22
32
10
22
12
12
24
16
8
24
32
16
24
10
24
32
10
32
10
16
12
32
10
24
16
8
1
4
16
1
4
7
2
1
7
1
1
1
1
22
7
40
28
25
18
18
7
5
14
1
7
2
1
3
6
3
4
19
1
3
4
4
9
25
-800
-54242
7400
94600
126000
630000
-31400
-72600
240000
125000
64000
189000
51000
-98633
-26000
199065
41067
41067
348800
-25821
15000
348800
-238633
-140000
346000
-2800
208000
-242633
-4000
-265633
-23000
-120821
-95000
-270033
-4400
-270033
-30000
6000
in65
No183
No183
No186
s9
s6
No190
in41
d4
No192
No192
No194
No194
No197
in62
d3
d3
d3
d3
No199
No199
No202
No204
No343
No206
in63
psg3
No209
in32
No211
in33
No213
in68
No215
in64
c4
No314
No219
334
FromNode
Pressure
To-Node
Name
To-Node
Pressure
25.43
23.43
23.43
25.46
25.52
35.43
25.45
24.28
25.46
23.84
23.84
23.97
23.97
39.79
39.14
30.63
30.63
30.63
30.63
25.08
25.08
26.18
39.87
39.64
24.79
24.79
18.62
40.13
40.11
40.19
39.9
25.4
24.64
40.26
40.24
40.33
30.02
34.19
No182
No156
in49
No369
No186
d3
No186
No190
No194
No344
No345
No192
No319
No204
No197
d4
No361
No357
No202
No365
in13
No206
No209
No204
No342
No206
psg2
No211
No209
No215
No211
No218
No213
c4
No215
No360
No218
No224
25.43
23.46
23.21
25.45
25.46
30.63
25.46
25.45
23.97
23.77
23.82
23.84
23.96
39.87
39.79
25.46
27.12
27.52
26.18
25.32
23.93
24.79
40.13
39.87
22.75
24.79
18.49
40.19
40.13
40.26
40.19
30.16
25.4
40.33
40.26
40.74
30.16
33.7
Result
Downstre
am
Velocity
(m/s)
0.2
2
1.6
2.4
3.2
17.9
0.8
15.1
10.4
5.5
2.8
8.3
2.2
1.2
3.3
8.1
5.3
5.2
11.7
2
4.9
12.3
2.9
6.8
13.3
0.6
9.9
2.9
0.5
3.2
2.8
9.4
13.5
3.2
0.5
5.7
2
1.4
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi217
Pi219
Pi218
Pi221
Pi220
Pi222
Pi223
Pi224
Pi226
Pi225
Pi227
Pi228
Pi229
Pi230
Pi231
Pi233
Pi232
Pi235
Pi234
Pi236
Pi238
Pi237
Pi240
Pi239
Pi242
Pi241
Pi244
Pi243
Pi245
Pi246
Pi248
Pi247
Pi250
Pi249
Pi251
Pi252
Pi253
Pi254
24
8
8
24
24
42
16
24
32
32
8
8
32
30
30
16
10
32
10
16
16
16
8
16
16
12
16
10
16
24
24
16
24
16
24
36
36
16
9
6
4
6
8
40
16
165
82
55
6
12
144
14
16
69
2
18
2
4
4
2
1
1
2
2
77
1
76
3
29
2
1
2
124
31
22
77
-150248
5000
-1000
119785
-119785
-1091608
96920
198439
481106
330503
1000
2000
151681
144788
135212
61920
-5000
265530
-11500
102000
82920
14000
15000
-42000
57000
-45000
66920
16000
-20685
301594
-190806
153000
-37806
0
156000
23090
100910
-20550
No222
No224
in56
No225
No225
d11
s8
s8
No227
No227
No229
No231
c5
d11
d11
No236
in17
No238
in86
No239
No241
No241
No244
in43
No246
in74
No247
No247
d13
d13
No250
No250
No252
No252
d12
d12
d12
d13
44.43
33.7
33.7
44.55
44.55
41.24
44.7
44.7
40.01
40.01
34.13
34.16
37.96
41.24
41.24
38.79
38.78
27.42
27.37
34.08
43.1
43.1
33.4
33.37
33.49
33.23
42.8
42.8
35.37
35.37
34.08
34.08
34.08
34.08
41.11
41.11
41.11
35.37
s8
in57
No224
No222
s8
s8
No241
d13
No363
No379
in54
No229
No233
d12
d12
d13
No236
No129
No238
No246
No247
No315
in44
No244
No244
No246
No236
in6
No249
No257
No249
No324
No250
ps35
No322
c3
No265
No249
44.7
33.62
33.7
44.43
44.7
44.7
43.1
35.37
33.32
38.01
34.12
34.13
36.78
41.11
41.11
35.37
38.79
26.8
27.42
33.49
42.8
43.09
33.28
33.4
33.4
33.49
38.79
42.76
35.84
34.92
35.84
33.43
34.08
34.08
36.61
41.11
41.07
35.84
335
Result
Downstre
am
Velocity
(m/s)
2.9
1.2
0.2
2.3
2.3
7.4
4.3
4.8
7
4.2
0.2
0.5
2
1.9
1.8
3.4
0.6
4.9
2.1
5.9
3.7
0.6
3.5
2.5
3.3
4.7
3.3
1.8
1.1
7.5
4.8
8.9
1
0
3.7
0.2
0.9
1.1
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi256
Pi255
Pi258
Pi257
Pi260
Pi259
Pi261
Pi263
Pi262
Pi265
Pi264
Pi267
Pi266
Pi269
Pi268
Pi271
Pi270
Pi272
Pi273
Pi274
Pi275
Pi276
Pi277
Pi279
Pi278
Pi281
Pi280
Pi283
Pi282
Pi284
Pi286
Pi285
Pi288
Pi287
Pi290
Pi289
Pi293
Pi292
24
10
24
12
24
16
32
36
16
10
10
10
10
10
10
18
16
30
30
10
30
16
16
36
36
30
8
36
36
30
16
16
24
18
16
16
16
10
5
11
5
6
3
7
18
26
7
2
3
2
1
1
2
3
25
28
204
15
12
75
98
19
81
3
19
79
29
6
7
54
2
2
2
2
32
5
242194
59400
72194
124000
196194
46000
265530
-79090
-180000
47900
-11500
7500
40000
47500
-400
37676
41549
225030
225030
-22800
242730
-10865
-1135
34000
23090
225030
-17700
10910
90000
225030
-25000
14135
30000
195030
185030
10000
-1135
0
No257
No257
No259
No259
No261
No261
No129
No265
No321
No267
in36
No268
No268
No271
in75
No272
No272
No273
No274
in87
No263
No272
No144
No279
c3
No281
in88
No283
s12
No275
No285
No285
No287
No287
No289
No289
No293
No292
34.92
34.92
34.24
34.24
34.44
34.44
26.8
41.07
38.42
28.25
28.17
26.88
26.88
27.35
27.35
37.34
37.34
25.55
24.49
24.44
26.16
37.34
37.31
41.1
41.11
25.66
21.13
41.1
41.14
14.64
37.47
37.47
14.25
14.25
12.81
12.81
37.31
37.31
No261
No267
No219
in18
No259
No323
No263
No283
No265
No271
No267
in37
in38
No268
No271
No144
No233
No274
No275
No263
No281
No285
No293
No307
No279
No273
No281
No279
No283
No287
s18
No286
in89
No289
in91
in90
No278
No293
34.44
28.25
34.19
28.03
34.24
34.23
26.16
41.1
41.07
27.35
28.25
26.86
26.55
26.88
27.35
37.31
36.78
24.49
14.64
26.16
25.66
37.47
37.31
41.09
41.1
25.55
25.66
41.1
41.1
14.25
37.53
37.32
14.24
12.81
10
12.8
37.32
37.31
336
Result
Downstre
am
Velocity
(m/s)
6.1
10.6
1.8
15.4
5
2.6
5
0.7
9.1
8.8
2.1
1.4
7.6
8.9
0.1
1.5
2.2
5.2
8.6
4.7
5.3
0.6
0.1
0.3
0.2
4.9
6.6
0.1
0.8
8.8
1.3
0.7
1.8
23.5
35.6
1.5
0.1
0
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi294
Pi295
Pi297
Pi296
Pi299
Pi298
Pi301
Pi300
Pi302
Pi304
Pi303
Pi306
Pi305
Pi307
Pi310
Pi309
Pi312
Pi311
Pi313
Pi314
Pi315
Pi316
Pi317
Pi318
Pi319
Pi320
Pi321
Pi322
Pi323
Pi324
Pi325
Pi326
Pi327
Pi328
Pi329
Pi330
Pi406
Pi407
16
24
16
12
16
12
16
12
16
16
12
36
16
24
36
36
36
20
24
24
32
32
18
24
28
28
16
24
8
24
24
24
24
12
17.25
36
22
24
5
4
13
4
21
3
18
4
23
31
6
8
5
45
32
220
44
208
3.3
1
1
1
5
2
31
5
25
5
4
10
2
9
115
15
1
1
87
1
0
-13000
0
0
0
0
0
0
38000
7000
31000
15000
-19000
-15000
0
0
0
0
34620
193768
656227
683227
169063
138000
-57788
322212
-150821
66194
1000
64194
-37806
-14135
0
-63209
127000
912584
-256194
64194
s13
s10
No298
in94
No300
in93
No302
ps36
s10
No304
No304
No307
No325
in76
No310
No310
No313
ps39
No2
No20
No24
No26
No56
No88
No92
No111
No218
No219
No229
No231
No239
No278
No286
d1
No105
c1
d9
No358
37.53
37.31
37.53
37.53
37.53
37.53
37.53
37.53
37.31
36.88
36.88
41.09
41.07
41.08
41.09
41.09
41.09
41.09
42.14
42.65
46.66
46.78
36.26
27.89
23.48
24.06
30.16
34.19
34.13
34.16
34.08
37.32
37.32
44.19
42.2
27.16
27.42
34.08
No300
No278
in95
No298
No302
No300
No298
No302
No304
in35
in34
No309
No307
No309
in77
No313
s11
No313
No21
No14
No35
No24
No59
No335
No140
No92
No379
No231
in55
No358
No252
No286
in73
s14
No320
No381
No105
No239
37.53
37.32
37.53
37.53
37.53
37.53
37.53
37.53
36.88
36.85
36.54
41.09
41.09
41.09
41.09
41.09
41.09
41.09
42.14
42.6
46.54
46.66
35.24
27.81
23.61
23.48
38.01
34.16
34.12
34.08
34.08
37.32
37.32
46.91
42.07
26.94
42.2
34.08
337
Result
Downstre
am
Velocity
(m/s)
0
0.3
0
0
0
0
0
0
2
0.4
2.9
0.1
0.9
0.3
0
0
0
0
0.7
3.9
6.7
7
7.4
4.3
1.6
8.8
9.8
1.7
0.2
1.6
1
0.3
0
4.9
5
13.2
9.7
1.6
Chapter 6
Name
Diameter
(in)
Length
(km)
Flow
(sm^3/hr)
FromNode
Name
FromNode
Pressure
To-Node
Name
To-Node
Pressure
Pi408
Pi412
Pi417
Pi418
Pi419
Pi420
Pi421
Pi422
Pi423
Pi424
Pi425
Pi426
Pi410
Pi411
Pi415
Pi416
Pi351
Pi352
Pi353
Pi354
Pi355
Pi356
Pi357
Pi358
Pi360
Pi361
Pi345
Pi346
Pi347
Pi348
Pi349
Pi350
Pi359
Pi362
Pi363
Pi364
Pi365
Pi366
Pi367
32
24
32
6
16
6
16
6.065
28
6.065
24
10
12
23.25
12
23.25
16
24
28
24
23.25
16
32
36
17.25
17.25
32
6.065
28
6.065
28
6.065
24
6.065
28
6.065
32
7.981
24
6
38
80
6
6
3
20
6
2
2
6
3
15
7
12
8
27
11
3
18
4
10
5
63
2
2
112.6
1
1
1
4.2
2
1
1
16.8
23
14.9
9.8
1
0
-270033
481106
0
-25821
0
-10821
0
94600
0
25158
0
41067
0
41067
0
-38811
-232041
-128159
-53119
-67329
28000
151681
912584
-114279
-476852
-72633
-2810
-572342
800
-525542
9000
206000
4420
303712
2410
597527
1970
201868
s1
No360
No363
No363
No365
No365
No367
No367
No369
No369
No371
No371
No357
No357
No361
No361
No233
No249
No377
No378
No377
No379
No379
No382
s4
s16
No366
No10
No368
No368
No372
No372
No374
No374
No376
No376
No383
No383
No41
34.08
40.74
33.32
33.32
25.32
25.32
25.03
25.03
25.45
25.45
23.38
23.38
27.52
27.52
27.12
27.12
36.78
35.84
24.97
25.03
24.97
38.01
38.01
50
35.55
35.79
39.59
39.57
33.48
33.48
31.14
31.14
27.99
27.99
27.16
27.16
37.09
37.09
42.86
No358
No222
d8
in96
No213
in97
No199
in98
No182
in99
No176
in100
d4
s19
d4
s20
No144
No233
d8
No123
No378
No326
c5
No105
s16
s5
No197
No366
d7
No1
No84
No4
No88
No6
No114
No8
No86
No11
No22
34.08
44.43
25.03
33.32
25.4
25.32
25.08
25.03
25.43
25.45
23.37
23.38
25.46
27.52
25.46
27.12
37.31
36.78
25.03
25.15
25.03
37.9
37.96
42.2
35.79
39.57
39.79
39.59
33.73
33.47
32.1
30.75
27.89
27.93
25.61
26.75
35.2
37.07
42.8
338
Result
Downstre
am
Velocity
(m/s)
0
5.7
9.5
0
2
0
0.9
0
2.4
0
0.9
0
5.6
0
5.6
0
2
5.6
3.3
1.9
2.5
1.4
1.9
8.2
5.4
22.3
0.9
1
10.9
0.3
10.8
4
6.4
2.2
7.6
1.2
8.3
0.4
4
Chapter 6
339
Chapter 6
340
Chapter 6
Suction
Pressure
(bar)
Discharge
Differential
Ratio
Power (KW)
Flow
(sm^3/h)
minimum
20
Current
26.94
Current
Maximum
50
65
Current
23.06
Maximum
45
Minimum
1.2
Current
1.82
Maximum
Current
20000
Installed
40000
Max. Available
40000
Fractional Capacity
.5931
Minimum
400000
Current
912584
Maximum
2000000
341
Chapter 6
Actual in
Count
Total km
4.000
1
1.00
6.000
14
80.04
6.065
7
36.00
7.981
2
10.77
8.000
21
134.28
10.000
45
205.56
12.000
15
132.00
16.000
57
997.97
17.250
3
5.05
18.000
15
218.00
20.000
3
256.00
22.000
7
138.00
23.250
3
19.00
24.000
86
909.44
28.000
25
240.37
30.000
16
374.00
32.000
35
1235.89
36.000
15
677.96
42.000
2
205.00
6.4.3.3.7 Pipe Length Summary
Model Reliability
Pipes
Greater
Than Or
Equal To
0.0
20.0
40.0
60.0
80.0
100.0
120.0
140.0
160.0
200.0
220.0
Pipes
Less
Than
Count
Total
km
20.0
40.0
60.0
80.0
100.0
120.0
140.0
160.0
180.0
220.0
240.0
306
32
8
9
5
3
1
2
3
2
1
1624.60
890.18
377.00
654.98
428.00
337.57
124.00
303.00
505.00
412.00
220.00
342
Chapter 6
Alpha:
343
Chapter 6
344
Chapter 6
6.4.4.4 PROFILES
A profile is a curve that defines a functional
relationship between two variables. These
variables are typically boundary condition
variables with respect to time. Each profile
is defined by sets of X and Y coordinates.
The variable represented on the X axis is
the independent variable, while the
variable represented on the Y axis is the
dependent variable.
345
Chapter 6
b) Pressure
g) Utilization factor
c) Gas temperature
2) Pipes:
a) Pipeline ambient temperature
b) Elevation
3) Regulators:
a) Setpoint
b) Cgmax
c) Maximum flow rate
4) Valves:
a) Dimensionless valve position
5) Compressor Stations:
a) Ambient temperature
b) Gas cost
c) Electric cost
d) Other cost
e) Set pressure
f) Maximum flow rate
346
Chapter 6
HOURS
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
FLOW
EVERY
HOUR
43.69
43.13
44.05
42.98
43.99
43.65
34.4
35.67
44.17
38.92
38.49
38.5
39.4
38.85
39.44
40.04
39.52
39.4
42.94
42.7
44.84
45.95
44.96
45.16
FACTOR
0.996
0.983
1.004
0.980
1.003
0.995
0.784
0.813
1.007
0.887
0.877
0.878
0.898
0.886
0.899
0.913
0.901
0.898
0.979
0.973
1.022
1.047
1.025
1.029
347
Chapter 6
Figure 6- 65
348
Chapter 6
349
Chapter 6
System linepack
350
Chapter 6
This model is ready to be used in optimization algorithms to find the best operation
scenario for the network.
351
Chapter 7 | Optimization
Chapter 7
Optimization
352
Chapter 7 | Optimization
7.1 Introduction
Optimization is the act of obtaining the
best result under given circumstances. In
design, construction, and maintenance of
any engineering system, engineers have to
take many technological and managerial
decisions at several stages. The ultimate
goal of all such decisions is either to
minimize the effort required or to
maximize the desired benefit.
353
Chapter 7 | Optimization
354
Chapter 7 | Optimization
355
Chapter 7 | Optimization
planning,
356
Chapter 7 | Optimization
7.2.1 Objectives
Maximize profit;
Minimize cost;
be
an
Maximize utility;
Maximize turnover;
Maximize return on investment;
Maximize net present value;
Maximize number of employees;
Minimize number of employees;
Minimize redundancy;
357
Chapter 7 | Optimization
358
Chapter 7 | Optimization
359
Chapter 7 | Optimization
Subjected
to
constraints.
conventional
linear
360
Chapter 7 | Optimization
by a variable t .
(ii) Represent the products xj t by variables
wj The objective now becomes
(iii)
Introduce
constraint
361
Chapter 7 | Optimization
7.2.2 Constraints
7.2.2.1 Productive capacity constraints
If resources to be used in a productive
operation are in limited supply, then the
relationship between this supply and the
possible uses made of it by the different
activities give rise to such a constraint.
362
Chapter 7 | Optimization
limitations
If there is a limitation on the amount of a
product that can be sold, which could well
result in less of the product being
manufactured than would be allowed by
the other constraints, this should be
modelled. Such constraints may be of the
form
xM
constraints
363
Chapter 7 | Optimization
364
Chapter 7 | Optimization
z ui 0, for all i,
z vi 0, for all i.
bounds
365
Chapter 7 | Optimization
written as
Optimizaiton
Problem types
Linear
Quadratic
Integer
Stochastic
Dynamic
366
Non-linear
Chapter 7 | Optimization
367
Chapter 7 | Optimization
X2
to:
X2 0
368
Chapter 7 | Optimization
369
Chapter 7 | Optimization
Convex Sets
We now define the important notion of a
convex set .
370
Chapter 7 | Optimization
Figure 7.5:
Maximize 5 X1 + 3 X2
Subject
to:
2 X1 + X2 40 labor constraint
X1 + 2 X2 50 material constraint
and both X1, X2 are non-negative.
371
Chapter 7 | Optimization
372
Chapter 7 | Optimization
373
Chapter 7 | Optimization
constraints:
{3X1 + X2 8, 4X1+ 3X2 19, X1+ 3X2 7,
X1 0, X2 0};
Next we plot the solution set of the
inequalities to produce a feasible region of
possibilities.
c) The 2:1 cost ratio of steak to potatoes
dictates that the solution must be here
since, as a whole, we can see that one unit
of steak is slightly less nutritious than one
unit of potatoes. Plus, in the one category
where steak beats potatoes in healthiness
(proteins), only 7 total units are necessary.
Thus it is easier to fulfill these units without
buying a significant amout of steak. Since
steak is more expensive, buying more
potatoes to fulfill these nutritional
requirements is more logical.
Standard Form
An integer linear program in canonical
form
is
expressed
as:[
374
Chapter 7 | Optimization
Example
The graph on the right shows the following
problem.
7.3.2.2 Applications
There are two main reasons for using
integer variables when modeling problems
as a linear program:
The integer variables represent quantities
that can only be integer. For example, it is
not possible to build 3.7 cars.
The integer variables represent decisions
and so should only take on the value 0 or 1
.
points
and
which both have an
objective value of 2. The unique optimum
of the relaxation is
with objective
375
Chapter 7 | Optimization
7.3.2.3 Algorithms
Exact
there are a variety of algorithms that can
be used to solve integer linear programs
exactly. One class of algorithms are cutting
plane methods which work by solving the
LP relaxation and then adding linear
constraints that drive the solution towards
being integer without excluding any
integer feasible points.
Scheduling
These problems involve service and vehicle
scheduling in transportation networks. For
example, a problem may involve assigning
buses or subways to individual routes so
that a timetable can be met, and also to
equip them with drivers. Here binary
decision variables indicate whether a bus
or subway is assigned to a route and
whether a driver is assigned to a particular
train or subway.
Telecommunications networks
The goal of these problems is to design a
network of lines to install so that a
predefined set of communication
requirements are met and the total cost of
the network is minimal. This requires
optimizing both the topology of the
network along with the setting the
capacities of the various lines. In many
cases, the capacities are constrained to be
integer quantities. Usually there are,
depending on the technology used,
additional restrictions that can be modeled
as a linear inequalities with integer or
binary variables.
Heuristic
many problem instances are intractable
and so heuristic methods must be used
instead. For example, tabu search can be
used to search for solutions to ILPs.To use
tabu search to solve ILPs, moves can be
defined as incrementing or decrementing
an integer constrained variable of a
feasible solution, while keeping all other
integer-constrained variables constant.
The unrestricted variables are then solved
for. Short term memory can consist of
Cellular networks
The task of frequency planning in GSM
mobile networks involves distributing
available frequencies across the antennas
so that users can be served and
interference is minimized between the
antennas. This problem can be formulated
as an integer linear program in which
binary variables indicate whether a
frequency is assigned to an antenna.
376
Chapter 7 | Optimization
minimizes
subject to
where xT denotes
the
vector transpose of X.
The
notation Ax b means that every entry of
the vector Ax is less than or equal to the
corresponding entry of the vector b.
neural networks
interior point,
active set
gradient projection,
Example
Solve a simple quadratic programming
problem: find values of x that minimize
f(x)=(1/2)x21+x22x1x22x16x2
subject to:
x1 + x2 2
x1 + 2x2 2
2x1 + x2 3
0 x1 , 0 x2.
Given:
a
real-valued, n-dimensional
vector c,
an n n-dimensional
real Symmetric matrix Q,
an m n-dimensional
matrix A, and
real
Q = [1 -1; -1 2]
f = [-2; -6]
A = [1 1; -1 2; 2 1]
b = [2; 2; 3]
377
Chapter 7 | Optimization
lb = zeros(2,1)
And
378
Chapter 7 | Optimization
Example
Portfolio Selection
An investor has $5000 and two potential
investments. Let xj for j = 1 and j = 2
denote his allocation to investment j in
thousands of dollars.
From historical data, investments 1 and 2
have an expected annual return of 20 and
16 percent, respectively. Also, the total risk
involved with investments 1 and 2, as
measured by the variance of total return, is
given by 2x1 2+x 2 2 +(x1+x2) 2 ,
Maximize f(x),
subject to:
gi(x) bi (i = 1, 2, . . . , m).
As in linear programming, we are not
restricted to this formulation. To minimize
f(x), we can of course maximize f(x).
Equality constraints h(x) = b can be written
as two inequality constraints h(x) b and
h(x) b. In addition, if we introduce a
slack variable, each inequality constraint is
transformed to an equality constraint.
Thus sometimes we will consider an
alternative equality form:
Maximize f(x),
subject to:
hi(x) = bi (i = 1, 2, . . . , m) xj 0 (j = 1, 2, .
. . , n).
subject to:
g1(x) = x1 + x2 5, x1 0, x2 0,
379
Chapter 7 | Optimization
are nondeterministic
380
Chapter 7 | Optimization
381
Chapter 7 | Optimization
382
Chapter 7 | Optimization
383
Chapter 7 | Optimization
Si=ti(Si+1 , Xi)
Ri=ri(si+1 , Xi)
T = t(X, S)
Since the input state of the system
inuences the decisions we make, the
return function can be represented as
R = r(X, S)
A serial multistage decision process can be
represented schematically as shown in Fig.
384
Chapter 7 | Optimization
Fig 7.13
si , i = 1, 2, . . . , n + 1.
385
Chapter 7 | Optimization
386
Chapter 7 | Optimization
7.4.2.1 Introduction
Evolutionary Algorithms
Simulated
Anneling
Ant Colony
Neural
Network
Fuzzy
Genetic
Algorithm
387
Chapter 7 | Optimization
(1)
7.4.2.2 Idea of SA
In (Dowsland, 1995) table 6.1 is presented
which shows how physical annealing can
be mapped to simulated annealing. It is
repeated here
System States
Feasible Solutions
Energy
Cost
Change of State
Neighbouring
Solutions
Where
Temperature
Control Parameter
Frozen State
Heuristic Solution
Thermodynamic Combinatorial
Simulation
Optimisation
(2)
Table 6. 1
388
Chapter 7 | Optimization
Initialize the
solution(T=Th)
For N (number of
iteration)
Calculate
(its the difference of
energy between
current and new state)
YES
If
<0
NO
Accept new state
with probability
Pr(accepted)=exp(E/KB.T
Decrease the
temperature
End
389
Chapter 7 | Optimization
Starting temperature
390
Chapter 7 | Optimization
Temperature Decrement
Once we have our starting and stopping
temperature we need to get from one to
the other. That is, we need to decrement
our temperature so that we eventually
arrive at the stopping criterion.
Final temperature
t = t
Where < 1
Experience has shown that should be
between 0.8 and 0.99, with better results
being found in the higher end of the range.
Of course, the higher the value of , the
longer it will take to decrement the
temperature to the stopping criterion.
391
Chapter 7 | Optimization
t = t / (1 + t)
Where is a suitably small value. I have
entered this formula on a spreadsheet
(available from the web site) so that you
can play around with the parameters, if
you are interested.
Hard Constraints
:
these
constraints cannot be violated in a feasible
solution. For example, in designing the
Jubilee Campus the space allocation cost
function could define a hard constraint as
not allowing a professor to occupy an
office that is smaller than a certain size.
Cost Function
Presented with a solution to a problem,
there must be some way of measuring the
quality of the solution.
Soft Constraints
:
these
constraints should, ideally, not be violated
but, if they are, the solution is still feasible.
Delta Evaluation
:
the
difference between the current solution
and the neighbourhood solution is
evaluated.
Partial Evaluation
:
a simplified
evaluation function is used that does not
give an exact result but gives a good
indication as to the quality of the solution.
Cooling
If you plot a typical cooling schedule you
are likely to find that at high temperatures
392
Chapter 7 | Optimization
393
Chapter 7 | Optimization
real ants
394
Chapter 7 | Optimization
395
Chapter 7 | Optimization
396
Chapter 7 | Optimization
397
Chapter 7 | Optimization
swimming.
Figure 7.15. Two membership functions relating to
swimming water temperature. Set A is the set defining
the fraction of all individuals who think the water
temperature is too cold, and Set B defines the fraction of
all individuals who think the water temperature is too
hot.
398
Chapter 7 | Optimization
7.4.5.1 Introduction
In machine learning, artificial neural
networks (ANNs) are a family of models
inspired by biological neural networks and
are used to approximate functions that can
depend on a large number of inputs and
are generally unknown. ANNs are
presented as systems of interconnected
"neurons" which exchange messages
between each other. The connections have
weights that can be tuned, making neural
nets adaptive to inputs and capable of
learning. The ability of the ANNs to
accurately
approximate
unknown
functions , underlies this growing
popularity. Neural networks have been
used to solve a wide variety of tasks that
are hard to solve using ordinary rule-based
programming, including computer vision
and speech recognition, Character
Recognition, Image Compression, Stock
Market Prediction, Traveling Saleman's
Problem, and Medicine, Electronic Nose,
Security, and Loan Applications . The word
network in the term 'artificial neural
network' refers to the interconnections
between the neurons in the different
layers of each system. An example system
has three layers (see figure 7.17). The first
layer has input neurons which send data
via synapses to the second layer of neurons
(hidden layer), and then via more synapses
to the third layer of output neurons. More
complex systems will have more layers of
399
Chapter 7 | Optimization
7.4.5.2.1 Back-propagation
Back-propagation, an abbreviation for
"backward propagation of errors", is a
common method of training artificial
neural networks used in conjunction with
an optimization method such as gradient
descent. The method calculates the
gradient of a loss function with respect to
all the weights in the network. The
gradient is fed to the optimization method
which in turn uses it to update the weights,
in an attempt to minimize the loss
function. Back-propagation requires a
known, desired output for each input value
in order to calculate the loss function
gradient.
400
Chapter 7 | Optimization
401
Chapter 7 | Optimization
In security
There are lots of applications for neural
networks in security too. Suppose you're
running a bank with many thousands of
credit-card transactions passing through
your computer system every single
minute. You need a quick automated way
of identifying any transactions that might
be fraudulentand that's something for
which a neural network is perfectly suited.
Your inputs would be things like
1. Is the cardholder actually present?
2. Has a valid PIN number been used?
3. Have five or more transactions
been presented with this card in
the last 10 minutes?
In airplanes
7.4.6.1 Introduction
Genetic Algorithms (GA) is a search
heuristic that
mimics
the
process
of natural selection. This heuristic is
In industry
402
Chapter 7 | Optimization
Optimization problems
In a genetic algorithm, a population of
candidate solutions (called individuals,
creatures, or phenotypes) to an
optimization problem is evolved toward
better solutions. Each candidate solution
has a set of properties (its chromosomes or
genotype) which can be mutated and
altered; traditionally, solutions are
represented in binary as strings of 0s and
1s, but other encodings are also possible.
Initialization
The population size depends on the nature
of the problem, but typically contains
several hundreds or thousands of possible
solutions. Often, the initial population is
generated randomly, allowing the entire
403
Chapter 7 | Optimization
Selection
Genetic operators
Main
articles:
algorithm) and
algorithm)
Crossover
Mutation
(genetic
(genetic
404
Chapter 7 | Optimization
Population:
The number of chromosomes available to
test.
Adaptation:
process of adjustment to given (outer)
conditions.
Fitness:
evaluation of an individual with respect to
its reproduction capability. Selection in EA
is based on the fitness. Generally, it is
determined on the basis of the objective
value(s) of the individual in comparison
with all other individuals in the selection
pool. The fitness function may additionally
depend
on
different
side
conditions/constraints and stochastic
influences. The term ``fitness function'' is
often used as a synonym for objective
function.
Representation:
Also coding or problem representation,
the mapping of the optimization problem
to the genome level.
Convergence:
Approach to a stationary system state,
where the variation rate of the considered
state variables goes to zero. It is necessary
to
distinguish gene
convergence (or
premature convergence) from the
convergence to optimal states. The latter is
divided into local and global convergence
depending on whether a local or a global
optimum is approached.
Fitness-proportionate selection:
Stochastic selection method where the
selection probability of an individual is
proportional to its fitness. Therefore,
fitness must be non-negative (scaling
function). Only fitness maximization is a
possible evolution goal.
Gene pool:
Total amount of all genes in a population.
Crossover rate:
Probability of a CO execution.
Generation:
(Natural) time unit of the EA, an iteration
step of the EA, a complete cycle comprising
the procreation and evaluation (fitness
determination and selection) of one or
more individuals.
405
Chapter 7 | Optimization
Mutation/ mutation operator:
Variation operator, which changes the
information contained in the genome of a
parent according to a given probability
distribution. In the case of bit strings, this
is realized by random negation (bit
mutation, mutation rate) of single bits. For
real-valued object parameters a random
number from a suitable distribution is
added (mutation strength).
Individual:
Carrier of the genetic information
(chromosome). It is characterized by its
state
in
the
search
space,
its fitness (objective function), and, where
applicable, a set of strategy parameters.
The individual is the unit of selection.
Mutation rate:
Mutation probability of a single
gene/object parameter of an individual.
With respect to binary representation, the
mutation rate is the probability of flipping
a single bit position. The mutation rate
does not directly relate to the mutation
strength/step size.
Intermediate recombination:
Production of one recombinant in realvalued search spaces by calculating the
centroid (center of mass) of parents, a
special case of (weighted) arithmetic
recombination.
Inversion:
Mutation operator, which extracts a
randomly
chosen
fraction
of
a chromosome and reinserts the genes in
reverse order. In ordering problems.
Multi-criteria optimization:
Also multi-objective optimization, vector
optimization, or multi-criteria decision
making. Optimization with regard to
multiple objective functions aiming at a
simultaneous improvement of the
objectives.
Multi-recombination:
Variation operator, which recombines
more than two parents in order to produce
one or more offspring.
406
Chapter 7 | Optimization
Object parameter:
Also search space parameter, genotype or
phenotype variables which directly
influence the fitness. The objective
function directly depends on the object
parameters. In general, EA aims at
identifying those object parameters which
maximize (or minimize) the objective
function.
Optimization:
Main application domain of EAs.
Optimization is the process of determining
or approximating object parameter sets
that maximize resp. minimize a
given objective function(s). A distinction is
made between global and local
optimization (local search). Optimal states
are always defined with respect to a
given neighborhood in the search space.
For global optima the neighborhood
includes the whole search space.
Depending on the structure of the search
space parameter, combinatorial, or
structure optimization are distinguished.
Phenotype:
Expression of the properties coded by the
individual's genotype.
The
expression/development of the phenotype
can additionally be influenced by
(stochastic) constraints. The precise
definition is mostly problem-dependent.
For parameter optimization the phenotype
is usually identical with the object
parameters, whereas for structure
optimization (e.g. of neural networks) the
phenotype represents a specific structure.
Population:
Pool of individuals exhibiting equal or
similar genome structures, which allows
the application of genetic operators.
Parent:
Individual which has been selected for
reproduction due to its fitness.
Performance:
General term concerning the evaluation of
the optimization behavior/quality of EAs. It
is often used synonymously for monitoring
the
evolution
of fitness or objective
function values over the generations
(fitness dynamics).
Population size:
Number of individuals in a population.
Problem size:
Number
of object
parameters,
components, or degrees of freedom of an
optimization problem.
Performance measures:
Local and global performance measures
are differentiated.
Progress rate:
Theoretical
performance
measure,
expected value of the distance change to a
predefined goal (usually, the optimal
solution) in one generation
407
Chapter 7 | Optimization
Rank-based selection:
Also ranking selection. The individuals'
selection probabilities are assigned
according to the individuals' rank which is
based on the objective function values.
Real-coded GA:
GA tailored for the optimization in realvalued search spaces. In contrast to the
canonical GA, the genome consists of (realvalued) object parameters, i.e., evolution
operates on the ``natural'' representation.
The real-coded GA employs special
recombination operators, which are hybrid
constructs
of intermediate
recombination and mutation. Real-coded
GAs can exhibit self-adaptive behavior.
Selection pool:
The part of the population containing all
individuals which constitute mutual
reproduction partners (parents).
Selection pressure:
Ratio of the best individual's selection
probability to the average selection
probability
of
all
individuals
in
the selection pool.
Recombination:
Variation operator, which merges the
genetic information of two or more parent
individuals producing one or more
descendants, the so-called recombinants.
In
canonical GAs,
normally
two
recombinants
are
produced
by crossover using two parents.
Roulette selection:
Randomized
variant
proportionate selection.
Stopping condition:
Also termination criterion or stop rule,
conditions
which
determine
the
termination of the evolutionary process
(examples: number of objective function
evaluations, maximum run time)
of fitness-
408
Chapter 7 | Optimization
409
Chapter 7 | Optimization
7.4.6.3 GA Operators
7.4.6.3.1 Selection
In selection the offspring producing
individuals are chosen. The first step is
fitness assignment. Each individual in the
selection pool receives a reproduction
probability depending on the own
objective value and the objective value of
all other individuals in the selection pool.
This fitness is used for the actual selection
step afterwards.
7.4.6.3.1.1 Selection operator types:
= ()
following equation:
2
_( +1)
The
Steps:
1. =
.
2. For
the
1
()
410
first
parent:
1 =
Chapter 7 | Optimization
3. For
=
the
remaining
parents:
()
fitness:
2. Generate
[0,1]
random
number:
()
(1)
reproduced
Steps:
1. Order the N individuals of P (t)
according to their fitness.
With,
411
Chapter 7 | Optimization
Individual 1
Individual 2
Crossover point
Child 1
Child 2
Double-point crossover:
In double-point crossover two crossover
positions are selected uniformly at random
and the variables exchanged between the
individuals between these points. Then
two new offspring are produced.
Single-point crossover:
Individual 1
Individual 2
Crossover points
Child 1
412
Chapter 7 | Optimization
Multi-point crossover:
For multi-point crossover, m crossover
positions ki[1,2,...,Nvar-1], i=1:m, Nvar:
number of variables of an individual, are
chosen at random with no duplicates and
sorted into ascending order. Then, the
variables between successive crossover
points are exchanged between the two
parents to produce two new offspring.
Individual 1
Individual 2
Sample/Mask
After crossover the new individuals are
created:
Individual 1
Child 1
Child 2
7.4.6.3.3 Mutation
By mutation individuals are randomly
altered. These variations (mutation steps)
are mostly small. They will be applied to
the variables of the individuals with a low
probability (mutation probability or
mutation rate). Normally, offspring are
mutated after being created by recombination.
Child 1
Child 2
413
Chapter 7 | Optimization
Ex:
Proto-child
operators:
Single flip mutation:
Child
Steps:
1. Select one gene randomly.
New =
2.
Ex:
New =
3.
Proto-child
Example:
=70%
Child
Min value = 70 %
value = 83 %
100%
Old
Max value =
New =
+
10083
83+
= 91.5%.
New =
Ex:
= 83 -
8370
2
= 76.5%
Proto-child
Child
414
Chapter 7 | Optimization
Thrors mutation:
Fitness-based reinsertion:
Produce more offspring than needed for
reinsertion and reinsert only the best
offspring.
Elitist combined
reinsertion:
with
fitness-based
7.4.6.3.4 Reinsertion
Once the offspring have been produced by
selection, recombination and mutation of
individuals from the old population, the
fitness of the offspring may be determined.
If less offspring are produced than the size
of the original population then to maintain
the size of the original population, the
offspring have to be reinserted into the old
population. Similarly, if not all offspring are
to be used at each generation or if more
offspring are generated than the size of the
old population then a reinsertion scheme
must be used to determine which
individuals are to exist in the new
population.
Reinsertion methods
Pure reinsertion:
Produce as many offspring as parents and
replace all parents by the offspring.
Uniform reinsertion:
Produce less offspring than parents and
replace parents uniformly at random.
Elitist reinsertion:
415
Chapter 7 | Optimization
7.5.1 GA toolbox:
416
Fitness Functions.
Variables.
Populations.
Crossover.
Mutation.
Plot.
Main GA code.
Chapter 7 | Optimization
7.5.5 Crossover
Crossover types:
Crossover heuristic.
Crossover scattered.
Crossover intermediate.
Crossover single point.
Crossover two point.
Crossover arithmetic.
You cant use the all types for your case
because every type has its method for
crossover.
7.5.4 Population
Structure:
Function pop=population(Nvar,FitnessFcn,options)
totalPopulationSize = sum(options.PopulationSize);
pop = cell(totalPopulationSize,1);
for i=1:totalPopulationSize
% main code
end
figure 7.20
417
Chapter 7 | Optimization
Figure 7.21
figure 7.22
418
Chapter 7 | Optimization
Structure
function xoverKids =
crossover_pop (parents,options,NVARS, ...
FitnessFcn,thisScore,thisPopulation)
function mutationChildren =
mutate_pop(parents ,options,NVARS, ...
FitnessFcn, state, thisScore,thisPopulation)
nKids =length(parents)/2;
xoverKids = cell(nKids,1);
index=1;
for I = 1 : nKids
r1 = thisPopulation { parents (index) };
r2 = thisPopulation { parents (index+1) };
index=index+2;
t=randi (NVARS,1,2);
t=sort(t);
child =
[r1(1:t(1)), r2(t(1)+1:t(2)), r1(t(2)+1:end)];
xoverKids { i } = child;
end
mutationRate=0.1;
mutationChildren = cell (length(parents),1);
for i=1:length (parents)
child = thisPopulation {parents (i) , :} ;
mutationPoints =
find(rand(1,length(child)) < mutationRate);
for j = 1:length(mutationPoints)
b= mutationPoints(j);
if child(b)==0
child(b)=1;
else
child(b)=0;
end
end
mutationChildren{i} = child;
end
Figure 7.23
Figure 7.24
7.5.6 Mutation
Mutation types:
Mutation for binary system.
Mutation for real number.
Every type has its own way to mutate
419
Chapter 7 | Optimization
Structure
function mutationChildren =
mutate_pop(parents ,options,NVARS, ...
FitnessFcn, state, thisScore,thisPopulation)
mutationRate=0.1;
mutationChildren = cell(length(parents),1);
for I = 1 : length (parents)
child = thisPopulation { parents(i), : };
mutationPoints =
find ( rand (1 , length(child)) < mutationRate) ;
for m=1:length(mutationPoints)
p=mutationPoints(m);
if p<=5 %% limit number for first variable
mut1=1:5;
child(p) = mut1 (randi(length(mut1),1,1));
else
mut2= zeros(1,5);
child(p) = mut2(randi(length(mut2),1,1));
end
end
mutationChildren{i} = child;
end
Figure 7.26
Figure 7.25
420
Chapter 7 | Optimization
'MutationFcn',@mutate_pop,'PlotFcn',@gaplotbestf ,
'Generations',30,'PopulationSize',20,
'SelectionFcn',@selectionroulette, ...
'StallGenLimit',6,'Vectorized','on');
numberOfVariables = 9;
[x,fval,reason,output] = ...
ga(FitnessFcn,numberOfVariables,options)
scores = zeros(size(x,1),1);
for i=1:size(x,1)
p=x{i}
% main code
scores(i)=V;
end
Figure 7.27
The function
( x , inputs)
contain
Figure 7.28
two
variables
Options of GA:
Population type is set custom because
you create the population
In the option you define every
function you have made:
Population
Crossover
Mutation
Plot: its set as default
421
Chapter 7 | Optimization
PopulationType: 'doubleVector'
PopInitRange: [2x1 double]
PopulationSize: 20
EliteCount: 2
CrossoverFraction: 0.8000
MigrationDirection: 'forward'
MigrationInterval: 20
MigrationFraction: 0.2000
Generations: 100
TimeLimit: Inf
FitnessLimit: -Inf
StallLimitG: 50
StallLimitS: 20
InitialPopulation: []
InitialScores: []
PlotInterval: 1
CreationFcn: @gacreationuniform
FitnessScalingFcn: @fitscalingrank
SelectionFcn: @selectionstochunif
CrossoverFcn: @crossoverscattered
MutationFcn: @mutationgaussian
HybridFcn: []
PlotFcns: []
OutputFcns: []
Vectorized: 'off'
Figure 7.29
422
Chapter 7 | Optimization
423
Chapter 7 | Optimization
Figure 7.30
Figure 7.33
Figure 7.31
424
Chapter 7 | Optimization
figure 7.34
figure 7.35
425
Chapter 7 | Optimization
figure 7.36
figure 7.37
figure 7.38
426
Chapter 7 | Optimization
Note:
figure 7.39
427
Chapter 7 | Optimization
Number of variables
The options fields
Diversity
Population Options
Fitness Scaling Options
Selection Options
Reproduction Options
Mutation Options
The Crossover Fraction
Setting the Maximum Number of
Generations
7.6.9.1 Diversity
Population
options
control
the
characteristics of the individuals in the
population. This section describes the
following options:
428
Chapter 7 | Optimization
429
Chapter 7 | Optimization
Figure 7.40
figure 7.42
Figure 7.41
430
Chapter 7 | Optimization
431
Chapter 7 | Optimization
Figure 7.44
Figure 7.45
432
Chapter 7 | Optimization
433
Chapter 7 | Optimization
7.8 Optimization
Implementation
Medium constraints:
constraints that can be accepted if you
have no choice and the hard
constraints arent broken.
Soft constraints:
constraints that if the solution broke it,
the solution will be accepted but it will
not be the best.
7.8.1.2 Constraints
Constraints differ from case to another
depending on what you want.
You can divide it to many section just like:
Hard constraints:
constraint that cant be broken and if
the solution breaks any one of them,
the solution must be ignored.
434
Chapter 7 | Optimization
Initial Population
Simulator
synergi script
Fitness Function
matlab sign
Scaling
simulation by synergi
Selection
Crossover
Mutation
synergi sign
Simulator
Fitness Function
Stopping Criteria
Yes
Optimum
Results
435
End
No
Chapter 7 | Optimization
436
Chapter 7 | Optimization
7.9.2 Example
Consider the following set of cities: [A, B, C,
D, E and F], The problem lies in finding a
minimal path passing from all vertices
once.
437
Chapter 7 | Optimization
0 0 0 0 0 0 0 | 1 1 1 1 1 1 1
Initial Population
Fitness Function
Selection
Crossover
Mutation
1 1 1 1 1 1 1 0 0 0 0 0 0 0
1 0 1 1 1 1 1 0 0 1 0 0 0 0
Stopping Criteria
End
7.9.3.4.1 initial population
each chromosome will have a different
sequence to pass through all cities.
Figure 7.46
1 1 1 1 1 1 1 | 0 0 0 0 0 0 0
438
Chapter 7 | Optimization
439
Chapter 7 | Optimization
7.9.4.2 Crossover
In this case we use single point crossover
due to the chromosome length.
flow direction
7.9.4.1 Population
Chromosome length will be three, contains
three diameters from the range that you
define:
440
Chapter 7 | Optimization
7.9.4.3 Mutation
h =
0.81 1 2
0.82 2 2
2
0.83 3 2
3 5
= + * h
= h
1
S.P = 1000
cost = S.P *
1
m
(hr/year)
Figure 7.51 mutation code
Initial Cost
V= 4 *( d 2 di 2 )* L
Total mass
m=*V
Cost =
m
1000
kg
* 590 $
7.9.4.5 Results
Pipes diameters
d1 = 24 in
d2 = 21 in
d3 = 18 in
441
Chapter 7 | Optimization
7.9.5.3 solutions
7.9.5.1 Problems
442
Chapter 7 | Optimization
443
Chapter 7 | Optimization
Dahshour (Co1)
Shabshir (Co4)
Fayoum(Co5)
UGDC(Co6)
444
Chapter 7 | Optimization
7.9.5.3.1.2 Unsteady state analysis
Pressures before optimization
The pressures of power stations are out of allowable range, so we have to choose which
compressors will make the pressures stable at the allowable range.
After we propose the eight locations, we use the GA to find the optimum solution for the
operation for 24 hr.
after optimization
After optimization the code choose three compressors to operate in the grid with minimum
required power for compressors operation.
By using these compressors, the pressures at all power stations are within the range.
But what about the power of compressor stations for 24hr, it will be illustrated in the figure,
every compressor has its own power according to the location of the compressor in the grid.
445
Chapter 7 | Optimization
7.9.5.3.1.3 Result pressures from optimization
pressures of compressor stations are stable for 24hr which increase the efficiency of the
operation of the national grid.
The natural gas reaches gas power stations at pressure 27 bar and steam power stations at 9
bar which Achieve our target from optimization.
446
Chapter 7 | Optimization
The power differs from compressor to other according to the pressure of the grid.
The power of compressor 4 is constant for 24hr but compressor 2 is working for few hours.
Linepack
Linepack in the national grid for 24hr is varying between 83.5 and 97 mmsm 3, this range is
suitable for our model.
This linpack results from the variation in power stations load and load of industries.
447
Chapter 7 | Optimization
7.9.5.3.2 Flow Reduction Strategy
7.9.5.3.2.1 Understanding the
notion
Here in this case study we aim to solve a
problem that already exist in the national
grid which is shortage of the amount of gas
in the grid which will deteriorate the
pressures at some power station
decreasing its efficiency and threatening
the national security .
Total Supply
(sm^3/hr)
5559499.961
5559499.961
448
Description
Pressure
krir
5
,
6
combined
elsyof gas
damanhor
combined
el
mahmoudia
combined
nubaria
combined
6 october gas0
korymat comb
banha gas
noth
cairo
combined
27
Flow
(sm^3/hr)
-90000
27
27
-50000
-46000
27
-68000
27
-380000
27
27
27
27
-150000
-123000
-95000
-168000
Chapter 7 | Optimization
south
cairo
combined
wadi hoof gas
talkha combined
talkha gas
west
damietta
gas
new damietta gas
damietta
combined
elshabab gas
port said gas
el
attf
new
combined
27
-56000
27
27
27
27
-11000
-125000
-64000
-140000
27
26.96
-138000
-208000
27
27
27
-30000
-14000
-138000
Pressure Flow
(sm^3/hr)
intergen 3
9
-138000
krir 1
9
-42000
abo qir new
9
-51000
abo qir steam
9
-87000
kafr el dawar 9
-74000
steam
west
cario 9
-63000
steam
shobra el khema 9
-33000
steam
teppin steam
9
-14000
talkha steam
9
-51000
west suez steam 9
-153000
mousa steam
9
156000.016
ataka steam
9
-46000.004
east port said 9
-180000
steam
al arish steam
9
-19000
abo
sultan 9
-28000.002
steam
korymat steam
9
126999.992
449
Description
krir 5,6 combined
elsyof gas
damanhor combined
el mahmoudia combined
nubaria combined
Pressure
27
27
27
27
18.44
6 october gas0
25.15
korymat comb
27
banha gas
noth cairo combined
26.22
25.43
26.31
26.43
26.53
talkha gas
west damietta gas
26.85
27
22.96
damietta combined
22.91
elshabab gas
port said gas
el attf new combined
27
27
27
Flow
-90000
-50000
-46000
-68000
380000
150000
123000
-95000
168000
-56000
-11000
125000
-64000
140000
138000
208000
-30000
-14000
138000
Chapter 7 | Optimization
7.9.5.3.2.3 Discussion of the code
Now we will run both codes
simultaneously ( matlab code and synergi)
450
1911232
2016100
Chapter 7 | Optimization
Before Optimization
No
1
2
3
NodeDescripti
on
natgas borg el
arab
natgas agamy
After Optimization
NodeFlow
(sm^3/hr)
-34620
-3790
natgas
west -7300
alex
town gas alex
-15000
-7300
-15000
natgas
kafr -8540
eldawar
town
gas -16000
portsaid
transgas delta -13000
5
6
natgas
kafr -8540
eldawar
town gas portsaid -16000
transgas delta
-40000
-2700
quesna(egyptgas) -2700
-40000
10
egyptgas delta(2)
-15400
11
-45000
12
egyptgas
abo -15400
zabal
egyptgas cairo 1
-45000
13
egypt
gas
delta(1)
quesna(egyptg
as)
egyptgas
delta(2)
egyptgas abo
zabal
egyptgas cairo
1
natgas cario
-15000
13
natgas cario
-15000
14
towngas cairo
-72000
14
towngas cairo
-72000
15
natgas
6 -35000
october
egyptgas cairo -45000
2
towngas
el -5000
ismailia
city gas
-124000
15
natgas 6 october
-35000
16
egyptgas cairo 2
-45000
17
towngas
ismailia
city gas
-124000
midor
medalic
star nex
19
4
5
6
7
8
9
10
11
12
16
17
18
19
20
18
and -64500
-1830
20
451
-13000
-40000
el -5000
midor
and -62249
medalic
star nex
-1825
Chapter 7 | Optimization
NO
Description
Flow
NO
Description
Flow
21
misr chemical
-6100
21
misr chemical
-5542
22
alex petrolum
-19700
22
alex petrolum
-18450
23
amoic
-10500
23
amoic
-10438
24
anrpc
-6770
24
anrpc
-6534
25
polyprobelin
-21800
25
polyprobelin
-20669
26
26
27
elmogammaa
for -46000
petrol
amryia petroleum
-26600
27
elmogammaa
for -42528
petrol
amryia petroleum
-25458
28
misr petroluem
-160
28
misr petroluem
-148
29
sidpec
-7000
29
sidpec
-6986
30
petrochemical
-27000
30
petrochemical
-25516
31
-7500
31
-6906
32
fenosa facilities
-4000
32
fenosa facilities
-3742
33
methanics
-23000
33
methanics
-22978
34
morgan
-31000
34
morgan
-29667
35
gapco
-7000
35
gapco
-6600
36
suez petroluem
-11500
36
suez petroluem
-11380
37
elnasr petroluem
-7500
37
elnasr petroluem
-7467
38
gas derivatives
-40000
38
gas derivatives
-37901
39
-129800
39
-124803
40
alexandria fert
-40000
40
alexandria fert
-37301
41
-72600
41
-68085
42
helwan fert
-50000
42
helwan fert
-49450
43
-42000
43
-36919
44
-15000
44
-13403
452
Chapter 7 | Optimization
NO
Description
Flow
No
Descripiton
Flow
45
amrya cement
-2720
45
amrya cement
-2420
46
dkhela steel
-80000
46
dkhela steel
-69161
47
max cement
-17550
47
max cement
-17432
48
48
49
ganob
elwadi -6000
cement
katmya cement
-7400
49
ganob
elwadi -5983
cement
katmya cement
-6427
50
helwan cement
-20000
50
helwan cement
-19401
51
national cement
-37900
51
national cement
-35627
52
kamaen el top
-75000
52
kamaen el top
-65232
53
53
54
54
55
suez cement
-1000
55
suez cement
-955
56
arabya cement
-1000
56
arabya cement
-893
57
egypt cement
-5000
57
egypt cement
-4578
58
-1870
58
-1787
59
alex sugar
-31200
59
alex sugar
-29469
60
nobaria sugar
-18030
60
nobaria sugar
-17895
61
alex IA
-18420
61
alex IA
-17702
62
dakhlia sugar
-26000
62
dakhlia sugar
-25948
63
repco gas
-2800
63
repco gas
-2706
64
64
65
65
damietta
pump -3964
station
gastec talkha
-645
66
gazl el mahala
66
gazl el mahala
-18551
67
67
cebaghet el mahla
-2467
68
10 of ram
-95000
68
10 of ram
-74481
69
hawamdya sugar
-13800
69
hawamdya sugar
-13062
70
egyptgas aswan
-10000
70
egyptgas aswan
-8500
-19500
453
Chapter 7 | Optimization
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
-22300
-14000
-17700
-45000
-400
-15000
71
72
73
74
75
76
-30000
-5000
-5000
-5000
-5000
-10000
-1000
-6000
-10000
-11500
-22800
77
78
79
80
81
82
83
84
85
86
87
454
-21330
-12811
-17134
-44978
-394
-11955
-27825
-4613
-4756
-4938
-4593
-8151
-952
-5316
-8296
-10966
-21786
Conclusion
Egyptian natural gas grid has problem in the power stations pressures, due to
low pressure the power stations works at low efficiency and consume more flow
than the required in normal operation.
So to solve this problem we have two solutions:
1. Converting power to gas
This solution will increase the calorific value of natural gas and increase the
amount of gas in the national grid by blending hydrogen in the grid at WDGC
with percentage equal to 3% from the total output flow from WDGC.
Increasing the calorific value of natural gas will decrease the total
consumption of power stations and that will increase the linepack and
pressure in the national grid.
2. Optimization
We have two scenarios first one is demand reduction
In this scenario we reduce the amount of the required flow for industrial
demands according to their priorities, it differs from industrial demand to
other.
The reduction of flow will increase the amount of gas flowing in the
national grid which increase the pressure in the grid and the linepack.
Our target from this scenario is to increase the pressure at all power
stations gas or steam to the required level and minimize the amount of
flow reduction.
Second scenario is to install new compressor stations
Install new compressors stations in the grid will increase the pressure in
the grid without the need to reduce the flow at any demand.
But increasing the number of new compressor stations in the grid above
the optimum limit will decrease the pressure in the grid instead of
increasing it, so to reach the optimum solution GA will search for the
optimum solution to increase the pressure at power stations to the
required level.
We installed eight compressor stations beside Dahshour compressor and
the locations are:
WDGC, Abohoumos, Shabshir, Cairo, Fayoum, UDGC, East El-Tina and
Suez.
455
Our target from this scenario is to increase the pressure at all power stations gas
or steam to the required level and minimize the total power of compressor
stations.
So after running the GA, the optimum solution was including three compressor
stations where are in Abohoumos, Shabshir and East El-Tina.
456
References
Chapter 1
Chapter 2
Links:
http://www.pv-magazine.com/news/details/beitrag/power-to-gas-platform-established-ineurope_100011076/#axzz4BgHhXjsY
http://www.powermag.com/progress-for-germanys-power-to-gas-drive/
http://blogs.worldwatch.org/from-wind-and-sun-to-gas-fraunhofer%E2%80%99s%E2%80%9Crenewable-methane%E2%80%9D-energy-storage-technology/
http://www.solar-energy-at-home.com/photovoltaic-solar-power.html
http://conceptselectronics.com/diodes/structure-pn-junction-diode/
http://www.odec.ca/projects/2008/full8e2/background.html
http://www.ontariosolarfarms.com/solar-panels.html
http://www.renewablegreenenergypower.com/solar-energy-pros-and-cons-pvsystems/#.UFtfyVEfo_d
https://www.environment.gov.au/system/files/resources/ef2c1cc7-07d8-4ed8-8f79816d36fb959e/files/desalination-summary.pdf
http://www.solidcell.com/techElectrolyzer.htm
http://www.nrel.gov/docs/fy13osti/51995.pdf
http://www.gtai.de/GTAI/Navigation/EN/Invest/Industries/Smarter-business/Smartenergy/power-to-gas.html
http://www.gtai.de/GTAI/Navigation/EN/Invest/Industries/Smarter-business/Smartenergy/power-to-gas,t=model-commercialization-projects-in-germany,did=547394.html
http://www.businessdictionary.com/definition/inflation.html
http://www.investopedia.com/terms/n/npv.asp
457
Chapter 3
Chapter 4
Chapter 5
Help of ArcGIS
User Guide of ArcGIS DESK 1, 2, 3
GIS and science.com
Bgis introduction to GIS .com
Chapter 6
Chapter 7
[Randy_L._Haupt,_Sue_Ellen_Haupt]_Pratical_Genetic
[Singiresu_S._Rao]_Engineering_Optimization_Theory
[H._Paul_Williams]_Model_Building_in_Mathematical
Genetic Algorithms and Engineering Optimization [Mitsuo Gen and Runwei Cheng]
web.mit.edu/15.053/www/AMP-Chapter-13.pdf
www.tik.ee.ethz.ch/sop/publicationListFiles/zitz1999a.pdf
www.utwente.nl/ctw/wem/education/afstuderen/Loucks_VanBeek/05_chapter05.pdf
courses.mai.liu.se/FU/MAI0083/Report_Mina_Nikanfs.pdf
Bioinformatics - Genetic Algorithm And Direct Search Toolbox - Matlab
458
Appendix
Compressors code steady state
Population
function pop=population(Nvar,FitnessFcn,options)
totalPopulationSize = sum(options.PopulationSize);
x=[0 1]; % dim. in "m"
n=Nvar;
pop = cell(totalPopulationSize,1);
for i=1:totalPopulationSize
pop{i}= randi(x,1,n);
end
Crossover
function xoverKids = crossover_pop(parents,options,NVARS, ...
FitnessFcn,thisScore,thisPopulation)
nKids =length(parents)/2;
xoverKids = cell(nKids,1);
index=1;
for i=1:nKids
r1 = parents(index);
index = index + 1;
r2 = parents(index);
index = index + 1;
v=zeros(1,NVARS);
for j = 1:NVARS
if(rand > 0.5)
b= thisPopulation{r1};
v(j)= b(j);
else
b= thisPopulation{r2};
v(j)= b(j);
end
end
xoverKids{i} = v;
end
mutation
function mutationChildren = mutate_pop(parents ,options,NVARS, ...
FitnessFcn, state, thisScore,thisPopulation)
mutationRate=0.1;
mutationChildren = cell(length(parents),1);
for i=1:length(parents)
child = thisPopulation{parents(i),:};
mutationPoints = find(rand(1,length(child)) < mutationRate);
mut=[0 1];
child(mutationPoints) = mut(randi([1 length(mut)],
length(mutationPoints),1));
mutationChildren{i} = child;
end
459
fitness function
function scores=cost_fitness(x,inputs)
%%%%%%%%%%%%%%%%%%%%%% no.of components of nodes and valves
ID3=fopen('C:\Users\Ahmed Bahaa\Desktop\comp file\text
file\components.txt','wt');
fprintf(ID3,'%f\n',inputs);
fclose(ID3);
iter=0;
Syn_sign=0.5;
ddddd=55;
flow=0;
power=10000000;
prs=0;
it=1;
for i=1:size(x,1)
p=x{i};
ID1=fopen('C:\Users\Ahmed Bahaa\Desktop\comp file\text
file\syn_input.txt','wt');
fprintf(ID1,'%f\n',p);
fclose(ID1);
%%%%%%%%%%%%%% matlab sign to make synergee start simulation
ID2=fopen('C:\Users\Ahmed Bahaa\Desktop\comp file\text
file\mat_firstsign.txt','wt');
fprintf(ID2,'%f\n',it);
fclose(ID2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%% waiting for simulation end in synergee
for wait=1:100000;
pause(0.1);
try
Syn_sign = textscanu('C:\Users\Ahmed Bahaa\Desktop\comp
file\text file\mat_sign.txt');
ddddd = cellfun(@str2double,Syn_sign );
end
if ddddd ==it+1;
break
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Reading pressure values of gas power stations
try
Syn_sign = textscanu('C:\Users\Ahmed Bahaa\Desktop\comp
file\text file\syn_output.txt');
pr = cellfun(@str2double,Syn_sign );
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%pressure of steam power stations %%%%%%%%
Try
Syn_sign = textscanu('C:\Users\Ahmed Bahaa\Desktop\comp file\text
file\steam.txt');
prs = cellfun(@str2double,Syn_sign );
end
460
461
GA main code
clc,clear
winopen('C:\Users\Ahmed Bahaa\Desktop\comp file\network model.lnk');
inputs=[9 19 16]; %%[comp psg pss]
FitnessFcn = @(x)cost_fitness(x,inputs);
options = gaoptimset('PopulationType', 'custom');
options = gaoptimset(options,'CreationFcn',@population, ...
'CrossoverFcn',@crossover_pop, ...
'MutationFcn',@mutate_pop,'PlotFcn',@gaplotbestf , ...
'Generations',20,'PopulationSize',20,'SelectionFcn',@selectionroulet
te, ...
'StallGenLimit',6,'Vectorized','on');
numberOfVariables = 9;
[x,fval,reason,output] = ...
ga(FitnessFcn,numberOfVariables,options)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% optimum simulation data %%%%%%%%%%%%%%%%%%%%
p=x{1};
ID1=fopen('C:\Users\Ahmed Bahaa\Desktop\comp file\text
file\syn_input.txt','wt');
fprintf(ID1,'%f\n',p);
fclose(ID1);
winopen('C:\Users\Ahmed Bahaa\Desktop\comp file\optimum.lnk');
462
fitness function
function scores=cost_fitness(x,inputs)
%%%%%%%%%%%%%%%%%%%%%% no.of components of nodes and valves
ID3=fopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\text
file\components.txt','wt');
fprintf(ID3,'%f\n',inputs);
fclose(ID3);
iter=0;
ddddd=55;
it=1;
for i=1:size(x,1)
p=x{i};
ID1=fopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\text
file\syn_input.txt','wt');
fprintf(ID1,'%f\n',p);
fclose(ID1);
%%%%%%%%%%%%%% matlab sign to make synergee start simulation
ID2=fopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\text
file\mat_firstsign.txt','wt');
fprintf(ID2,'%f\n',it);
fclose(ID2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%% waiting for simulation end in synergee
for wait=1:100000;
pause(0.1);
try
Syn_sign = textscanu('C:\Users\Ahmed Bahaa\Desktop\comp
unsteady\text file\mat_sign.txt');
ddddd = cellfun(@str2double,Syn_sign );
end
if ddddd ==it+1;
break
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Reading pressure values of gas power stations
T1 = readtable('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\press
excel\pressuresg.csv');
B1= table2array(T1);
A1=B1(:,2:20);
u1=size(A1,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%pressure of steam power stations %%%%%%%%%
T2 = readtable('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\press
excel\pressuress.csv');
B2= table2array(T2);
A2=B2(:,2:17);
u2=size(A2,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% power of compressors %%%%%%%%%%%%%%%%%%%%%
463
464
GA main code
clc,clear
winopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\usmodel.lnk');
inputs=[9]; %%[compressors]
FitnessFcn = @(x)cost_fitness(x,inputs);
options = gaoptimset('PopulationType', 'custom');
options = gaoptimset(options,'CreationFcn',@population, ...
'CrossoverFcn',@crossover_pop, ...
'MutationFcn',@mutate_pop,'PlotFcn',@gaplotbestf , ...
'Generations',30,'PopulationSize',20,'SelectionFcn',@selectionroulet
te, ...
'StallGenLimit',6,'Vectorized','on');
numberOfVariables = 9;
[x,fval,reason,output] = ...
ga(FitnessFcn,numberOfVariables,options)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% optimum simulation data %%%%%%%%%%%%%%%%%%%
p=x{1};
ID1=fopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\text
file\syn_input.txt','wt');
fprintf(ID1,'%f\n',p);
fclose(ID1);
winopen('C:\Users\Ahmed Bahaa\Desktop\comp unsteady\optimum.lnk');
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479