Advanced Process Analysis For Source Red
Advanced Process Analysis For Source Red
' ( ' ) * #
* ) ' ( # ) '+ (
, #
. / #
.
.
0 ) #
1 +# 2 0 3 4
0 5
6 78 # 9 (+
0 * 5
': : 4
0 ) #
.
) 2 8 * - #
5 * )
* - $ ;%%$
) 2 8 * - #
0 5
6 , < ===$%
'6*,#' ,
0 *, > ,'60 *
0 *, > 5# *
2'(, # ; 0 , #',5# # ? 9 $;
' ( ' * $;
' )> 0 > 8 $
/ )> 0 > 8 $
$&
( ( $"
* ) ' '+ ( ;;
'+ ( ;
* ) ' '+ ( ;
, )'+ # ;"
# . ) '+ ) -
4 ( %
# . ) '+ ) -
4 6 ( %
) )( ? -
+ "
( "
=
*
# # * $$
'(( *
' *, .', > > > ,#', > * ?'# >5*
*( * ,2 ','03*, (2'* '$
Advanced Process Analysis System was successfully applied to the 15,000 BPD
alkylation plant at the Motiva Enterprises Refinery in Convent, Louisiana. Using the
of the Advanced Process Analysis System an average increase in the profit of 127% can
be achieved. Energy savings attained through reduced steam usage in the distillation
obtainable through pinch analysis. Also one of the main sources of the pollution from
Besides economic savings and waste reductions for the alkylation process,
which is one of the most important refinery processes for producing conventional
optimization of the process has been demonstrated. Application to other processes could
The alkylation plant model had 1,579 equality constraints for material and
energy balances, rate equations and equilibrium relations. There were 50 inequality
constraints mainly accounting for the thermodynamic feasibility of the system. The
equality and inequality constraints described 76 process units and 110 streams in the
plant. Verification that the simulation described the performance of the plant was shown
using data validation with 125 plant measurements from the distributed control system.
v
List of Tables
Table 5.4 Alkylation Plant Raw Material/Utility Costs and Product Prices 102
vi
List of Figures
Figure 1.4 Composite Curves for Hot Streams and Cold Streams 8
vii
Figure 5.2 Network Grid Diagram for Alkylation Process 107
Figure 5.3 Integrating columns (5C-601 and 5C-603) with the process:
Pressure shift for column 5C-601 only (left),
for both columns (right). 108
viii
! "
# $
$ $ "
& "
' ( )
' ( $ * + # !
! , " !
! , )
& " - . #
/ # #
$ " #
" $
$ 0 "
&
# $ # # 1 " "
# )
" #
!
$%
& !'
) "
!" #
! ("
! ! )
* + %
#
#
!"!
* + %
#
, % $ $
' ' % &
- ' $
' % & %
'
' $
1 " 31 #
2
$ " 5
$ "
"
$ #
6 7 &"
# 48 '
) "
# $ $
"4 " #
# # $ & #
$ $ & $
& $ 2:8"
" ? #
@ <<="
*) !"
# 1 #
# $ "
# # #
$ # " #
4
" # & 3
0 0 $ $
" #
) 1 #
$ $ "
# # "1 &
& $
0 #
$ & "
1 #
# "
BC $ ? 7" 1 &
$ $ "
7 & D" )
+ $ & !'
5 # %
$ &
$ $ "
1 "2" $
" $
A
# $ "
$ $ "
# $
"
1 "23 5 5
. /
! ) 0
/ 0
0 /
:
, !
) ) 0 $
# 1 "4" $ $
&
1 ) "
$ $
) 0 #
# ) 0
0 ) 0
" 0
" " # )# 0
# " F "
* *
C !0 C )! 0
C !0 C )! 0
. $$
1 .
.
1& . .
1 . 1& . $$ .
. $$
4) 1 .
1 "43 ? 5
E
-
# <>;F
& # $ &
& $
$ & # " 6
G "
6 # 1 "A
1 ": "
$ $
$ " 1
&
$ # $ #
0 &
# "
G #
"
# #
$ " /
# $ "
) $ $
$ $
" $
>
# $ "
* H#
"
461 461
451 451
* :=-< 31 * :=-< 31 I 2
* I*2
21 21
1 1
1 411 511 711 211 811 1 411 511 711 211 811
9 :;< 9 :;<
* ,
*2 * *
* 2
* * & !
1 ":3 C ?
$ $
"
=
. " !
B 6 & 7" 6 $ /
/ $ $ #
"
J 6 I
C ) 5 ' (
6 0 F 6 &F #
" & 3
& $ "
# " /
$ # # #
& "
/ " 6 $
" $
"
0 $ ".
$ $ 0 "
<
/ 0
$ " 6
$ # 4) :
) ' <<4("
$ $
$ 0 $ "
'*2 5A(
'*1(
$ 0 " /
# & " /
# " & $
'*2 5A(
'*1( " $ 2; #
*2 5A *1"
*1 $
*2 5A *1 $ $
" *1
;
*1 0
1 "!!
# #
# $ $ "
# "
6 & 0
" $
"
"
"
In this chapter the current status of the methodology and literature is reviewed
for the methods used in the Advanced Process Analysis System. Also the current
1.1, includes chemical reactor analysis, process flowsheeting, pinch analysis and on!line
optimization. All of these programs use the same information of chemical processes
(material and energy balances, rate equations and equilibrium relations). Consequently,
The need of an integrated approach to process analysis has been given by Van
Reeuwijk et al. (1993) who proposed having a team of computer aided process
described by Ballinger et al. (1994) called ‘ whose goal is have to a user interface
to create and manipulate objects such as processes, streams and components with
environment. The Clean Process Advisory System (CPAS) has been described by Baker
university team. This includes technology selection and sizing, potential and designs,
article by Shaney (1995) describes the various modeling software and databases
12
available for process analysis and design. A review of computer aided process
engineering by Winter (1992) predicts linking various applications will result in better
quality of process design, better plant operations and increased productivity. It also
multiple sources into a common multi!user framework for analysis, process definition
and process engineering documentation rather than the original concept of a common
user interface and datastore linked with a range of applications computing tools.
!! "# $ #! %
the new computing power for process optimization and control that leads to higher
product qualities and better processes, which are cleaner, safer, more efficient, and less
costly.
reported in the literature from 1983 to 1991 and cited nine applications – five ethylene
plants, a refinery, a gas plant, a crude unit and a power station. The results showed a
sulfuric acid plant of IMC Agrico at Convent, Louisiana. Economic optimization can
achieve a 17% increase in plant profit and 25% reduction in sulfur dioxide emission.
The plant was studied by Chen (1998), in developing the optimal way to conduct on!
13
line optimization. Also, Chen (1998) reported a number of other successful applications
& "# $ #! %
based on plant’s current operating and economic conditions. As shown in Figure 1.2,
! Data Reconciliation
! Parameter Estimation
! Optimization Algorithm
The relationship between these key elements is outlined in Figure 2.1. Both
plant model and optimization algorithms are required in the three steps of on!line
which are solved by the optimization algorithm. For data validation, errors in plant
and a test statistic is used to detect gross errors in the measurements. For parameter
14
constraints in the plant model. For economic optimization, the plant model is used with
economic model to maximize plant profit and provide optimal setpoints for the
' (
Heat Exchanger Network Synthesis (HENS) for maximum heat recovery is the
key to energy conservation in a chemical plant. The problem of design and optimization
of heat exchanger networks has received considerable attention over the last two
decades (Ahmad et al., 1990, Duran et al., 1986, Linhoff et al. 1978, 1979, 1982).
network to exchange heat among a set of process streams where any heating and
cooling that is not satisfied by exchange among these streams must be provided by
external utilities (Shenoy, 1995). Attempts at solving this problem have been based on
Rudd (1969). At that time, the design methods were generally based on heuristic
15
approaches. One of the commonly used rules was to match the hottest stream with the
coldest stream. Several other methods were based on tree search techniques (Lee et al.,
Hindmarsh (Linnhoff et al., 1979) first introduced Pinch Analysis; a method based on
formed the basis for further research. These concepts were reviewed by Gundersen et al.
Paoulias and Grossman (1983) formulated Maximum Energy Recovery (MER) problem
as a linear programming (LP) model based on the transshipment model, which is widely
used in operations research. This model was expanded to make restricted hot and cold
16
( ) * Many different approaches (Telang, 1998) have
been suggested to deal with this problem. One of these is the Waste Reduction
Algorithm (WAR) (Hilaly, 1994). The WAR algorithm is based on the generic pollution
algorithm can be used to minimize waste in the design of new processes as well as
potential environmental impacts, and it can be used in the design and modification of
caused by the energy and material that the process takes from and emits to the
environment. The potential environmental impact is a conceptual quantity that can not
The generic pollution balance equation of the WAR algorithm is now applied to
the conservation of the Potential Environmental Impact in a process. The flow of impact
, in and out of the process is related to mass and energy flows but is not equivalent to
17
= − + (2.3)
where is the potential environmental impact content inside the process, is the
input rate of impact, is the output rate of impact and is the rate of impact
generation inside the process by chemical reactions or other means. At steady state,
0= − + (2.4)
relates the conceptual impact quantity to measurable quantities. The input rate of
= ∑ = ∑ ∑ Ψ (2.5)
where the subscript ‘in’ stands for input streams. The sum over j is taken over all the
input streams. For each input stream j, a sum is taken over all the chemical species
present in that stream. Mj is the mass flow rate of the stream j and the xkj is the mass
chemical k.
The output streams are further divided into two different types: Product and
potential impact and all product streams are considered to have zero potential impact.
18
The output rate of impact can be written as:
= ∑ = ∑ ∑ Ψ
(2.6)
where the subscript ‘out’ stands for non!product streams. The sum over j is taken over
all the non!product streams. For each stream j, a sum is taken over all the chemical
species.
Knowing the input and output rate of impact from the equations 2.5 and 2.6, the
generation rate can be calculated using equation 2.4. Equations 2.5 and 2.6 need values
Ψ = ∑α Ψ (2.7)
,
where, the sum is taken over the categories of environmental impact. αl is the relative
impact of chemical k for impact of type l. Values of Ψsk,l for a number of chemical
species can be obtained from the report on environmental life cycle assessment of
products (Heijungs, 1992). Some non!zero values of Ψsk,l for the components used in
There are nine different categories of impact. These can be subdivided into four
photochemical oxidant formation), three human toxicity effects (air, water and soil) and
two ecotoxicity effects (aquatic and terrestrial). The relative weighting factor αl allows
the above expression for the impact to be customized to specific or local conditions.
The suggested procedure is to initially set values of all relative weighting factors αl to
19
one, and then allow the user to vary them according to local needs. More information on
impact types and choice of weighting factors can be obtained from the report on
equation is used to define two categories of Impact Indexes. The first category is based
on generation of potential impact within the process. These are useful in addressing the
questions related to the internal environmental efficiency of the process plant, i.e., the
environmental impact. The second category measures the emission of potential impact
by the process. This is a measure of the external environmental efficiency of the process
i.e. the ability to produce the desired products while inflicting on the environment a
minimum of impact.
Within each of these categories, three types of indexes are defined which can be
used for comparison of different processes. In the first category (generation), the three
20
1) This measures the the total rate at which the process generates potential
subtracting the input rate of impact ( ) from the output rate of impact
( ), i.e. = ! .
manufacturing a unit mass of all the products. This can be obtained from
= −
.
∑
3) This is a measure of the mass efficiency of the process, i.e., the ratio of
− ( ) − ( )
∑ ∑ −∑ ∑
= −
.
∑
4) This measures the the total rate at which the process outputs potential
− ( )
2.6, i.e. = ∑ ∑ Ψ .
21
5) This measures the potential impact emitted in manufacturing a unit mass of
all the products. This is obtained from dividing by the rate at which the
− ( )
∑ ∑
= −
.
∑
basis whereas indices 2, 3, 5 and 6 can be used to compare them independent of the
plant size. Higher values of indices mean higher pollution impact and suggest that the
plant design is inefficient from environmental safety point of view. Negative values
mean that the input streams are actually more harmful to the environment than the non!
" 1
average quality of all components available to the gasoline pool, increased amounts of
gasoline per volume of crude oil and high heats of combustion. Alkylates permit use of
internal combustion engines with higher compression ratios and hence the potential for
22
increased miles per gallon. Alkylates burn freely, promote long engine life, and have
low levels of undesired emissions (Albright 1990a, Corma and Martinez 1993).
containing tertiary hydrogen to an olefin. The process is used by the petroleum industry
to prepare highly branched paraffins mainly in the C7 to C9 range for use as high!quality
fuels for spark ignition engines. The overall process is a composite of complex
1 )
distill in the gasoline range (up to ca. 200 oC). Commercial refinery plants operate with
the C3 and C4 hydrocarbon streams; alkylation involving high molecular weight olefin
or isoparaffins (over C5) are not attractive, partly because of numerous side reactions
temperatures are reduced to the range of 0!10 oC. Cooling requirements are obtained by
flashing of unreacted isobutane. Some form of heat removal is essential because the
heat of reaction is approximately 14 x 105 J/kg (600 Btu/lb.) for butenes. In order to
zone. Isobutane!to!olefin molar ratios of 6:1 to 14:1 are common. More effective
23
The alkylation reaction system is a two!phase system with a low solubility of
isobutane in the catalyst phase. In order to ensure intimate contact of the reactant and
the catalyst, efficient mixing with fine subdivision must be provided. Presence of
unsaturated organic diluent in the acid catalyst favors the alkylation reaction. The
organic diluent has been considered to be a source of carbonium ions that promote the
" 1
More than 60% of the worldwide production of alkylate using a sulfuric acid
catalyst is obtained from effluent refrigeration process of Stratco Inc. A typical process
24
The Stratco reactor or contactor, shown in Figure 2.3, is a horizontal pressure
vessel containing a mixing impeller, an inner circulation tube and a tube bundle to
remove the heat generated by the alkylation reaction. The hydrocarbon and acid feeds
are injected into the suction side of the impeller inside the circulation tube. The impeller
rapidly disperses the hydrocarbon feed with the acid catalyst to form an emulsion. The
increase the overall heat transfer coefficients; improved pump/agitator system and
injection devices for introducing the hydrocarbon feed and the acid into the contactor,
which in turn improves alkylate quality and lowers refrigeration costs. The
pump/agitator has been positioned below the centerline in order to minimize partial
A part of the emulsion is continuously removed and sent to the acid settler or
decanter, where the acid and hydrocarbon phases separate. Then in a flash drum the
hydrocarbon phase is flashed to separate C4 and lighter components from the heavier
25
hydrocarbons. The lighter components are sent to the contactor as refrigerants and are
then recycled. The heavier components are sent to the deisobutanizer column, where
alkylate is separated from unreacted butane and isobutane. The C3 content in the system
is decreased by sending the lighter components through the depropanizer column. The
' ) " 1
tertiary carbon atom are found to undergo catalytic alkylation with olefins. Reactions
The principal reactions that occur in alkylation are the combinations of olefins
CH3 CH3
| |
CH2 = CH ! CH3 + CH3 ! CH ! CH3 → CH3 ! CH ! CH2 ! CH2 ! CH3
|
CH3
26
Steps in the alkylation reaction mechanism involving carbonium ions are shown
below, with typical examples to illustrate each reaction step (Cupit 1961). (X is OSO3H
or F):
1. The first step is the addition of proton to olefin molecule to form a tertiary butyl
cation:
C X
| | +
C ! C = C + HX C! C ! C C ! C ! C + X! (2.8)
| |
C C
2. Then, the tertiary butyl cation is added to the olefin :
C
| +
C=C!C → C!C!C!C!C (2.9)
| | |
C C C
C C
| | +
C – C+ + C!C=C!C → C!C!C!C!C (2.10)
| | |
C C C
C
| +
C=C!C!C → C!C!C!C!C!C (2.11)
|
C
27
3. Carbonium ions may isomerize via hydride and methyl shifts to form more stable
carbonium ions:
C C C
| |
~CH3! |
→ C!C!C!C!C (2.12)
+
C C C C C C
| | + ~H!
| | ~CH3!| |
C!C!C!C!C → C!C!C!C!C → C!C! C!C!C (2.13)
| | + +
|
C C C
C C C C
|
~CH3! | ~H!
| |
→ C!C!C!C!C → C!C!C!C!C (2.14)
| + | +
C C
4. These carbonium ions suffer rapid hydride transfer from isobutane, leading to the
C C C C C C C C
| | | | | | | |
C!C!C!C!C + C!C!C → C!C!C!C!C + C!C!C (2.15)
+ +
Unfortunately, these are not the only reactions occurring during alkylation. There
are a number of secondary reactions that in general reduce the quality of the
5. Polymerization results from the addition of a second olefin to the carbonium ion
formed in the primary reaction, as was seen in step 2 of the above mechanism:
28
The iC12H25+ can continue to react with an olefin to form a larger isoalkyl cation:
lower and a higher molecular weight isoparaffin than the initial one:
7. Larger isoalkyl cations can crack, leading to smaller isoalkyl cations and olefins:
iC12H25+
alkylated with olefins other than butanes. At the same time, saturated paraffin of the
same carbon number as the olefin is obtained. The reaction scheme for pentene, for
example is :
+ +
C!C!C!C +C!C!C C!C!C!C + C!C!C (2.22)
| | | |
C C C C
+
C!C!C C ! C = C + H+ (2.23)
| |
C C
29
C
+ | +
C!C!C+C!C=C → C!C!C!C!C (2.24)
| | | |
C C C C
The above reactions are believed to be fundamental to the alkylation process and are
used to explain the formation of both primary and secondary products. Although
isobutane and butenes were used as examples, these reactions also applied to other
- , ) " 1 " 2 3 ) !
Langley and Pike (1972) studied the sulfuric acid alkylation of isobutane with
Schmering carbonium ion mechanism with modification introduced to account for iC9
continuous flow stirred tank reactor at the temperature range 18!57 oC. The model was
found to be valid in the range of 27!57 oC using 95% sulfuric acid catalyst. Lower
/ , ) " 1 " 2 3 ) 4
Vichailak (1995) extended the mechanism given by Langley and Pike (1972) to
a nineteen!reaction mechanism model for isobutane with butylene alkylation and twenty
one!reaction mechanism model for isobutane with pentylene alkylation. The reaction
mechanism model for alkylation of isobutane with butylene is shown in Table 2.3. The
reaction rate constants (frequency factor and activation energy) for each reaction of
30
these mechanisms that were taken to be identical to the reactions in the propylene
mechanism. A list of these rate constants, k values, are given in Appendix C.8. The
results from this reaction model showed a fair agreement (90%) with the published data
(Albright, 1992, Srichanachaikul, 1996). The reaction mechanism and rate constants
shown in Table 2.3 to describe the alkylation reaction were used for the STRATCO
31
Table 2.2. Reaction Mechanism and Material Balances for Sulfuric Acid Alkylation of Isobutane with
Propylene (Langley, 1969)
+ − + − + −
8
+ 4
→
6
8
+ 4 6
= 4
[ 6
][ 4
]
+ −
7
= 5
[ 7
][ 4
]
+ −
Destructive alkylation reactions 8
= 6
[ 8
][ 4
]
+ − = + −
7
→
12
7
+ 9
= 7
[ 9
][ 4
]
= + − = + − + −
7
+ 4
→
13
5
+ 6 10
= 8
[ 10
][ 4
]
= + −
5
+ →
14
5
(c) Olefinic intermediate rate equations.
− + − + − = + −
5
+ + 4
→
3
5
+ 4 4=
=0= 9
[ 4
]− 10
[ 4
][ 4
]
+ − + − = + − + −
6
+ 4
→
4
6
+ 4 5=
=0= 13
[ 7
][ 4
]+ 17
[ 10
]−
+ − = + − = = + −
7
+ 3
→
15
10 14
[ 5
][ ]− 16
[ 5
][ 4
]
+ − + − + − = + −
10
+ 4
→
8
10
+ 4 7=
=0= 12
[ 7
]− 13
[ 7
][ 4
]
= + − + −
5
+ 4
→
16
9
(d) Carbonium ion rate equations.
+ − + − = + −
9
+ 4
→
7
9
+ 4 3+ −
= 0 = 1[ 3
][ ]− 2
[ 3
][ 4
]
+ − = + − + − = + −
10
→
17
5
+ 5 4+ −
=0= 4
− 9
[ 4
]− 10
[ 4
][ 4
]−
+ − = = + −
11
[ 4
][ 3
]− 13
[ 7
][ 4
]−
= + −
16
[ 5
][ 4
]
= + −
5+ −
=0= 14
[ 5
][ ]+ 17
[ 10
]−
+ −
3
[ 5
][ 4
]
= + − + −
6+ −
=0= 13
[ 7
][ ]− 4
[ 6
][ 4
]
+ − = + −
7+ −
=0= 11
[ 4
][ 3
]− 5
[ 7
][ 4
]−
+ − = + −
15
[ 7
][ 3
]− 12
[ 7
]
= + − + −
8+ −
=0= 10
[ 4
][ 4
]− 6
[ 8
][ 4
]
= + − + −
9+ −
=0= 16
[ 5
][ 4
]− 7
[ 9 [ 4
]
+ − = + −
10 + −
=0= 15
[ 7
][ 3
]− 17
[ 10
]−
+ −
32 8
[ 10
][ 4
]
Table 2.3. Reaction Mechanism and Material Balances for Sulfuric Acid Alkylation of Isobutane
with Butylene (Vichailak, 1995)
+ − + − + − =
6+ −
=0= 17
[ 11
]− 4
[ 6
][ 4
]− 19
[ 6
][ 4
]
+ − + − + − = + − + − =
6
+ 4
→
4
6
+ 4 7+ −
=0= 13
[ 4
][ 8
]− 5
[ 7
][ 4
]− 15
[ 7
][ 4
]
+ − = + − = + − + − =
6
+ 4
→
19
10 8+ −
=0= 11
[ 4
][ 4
]+ 10
[ 4
][ 4
]−
+ − + − + − + −
10
+ 4
→
8
10
+ 4 6
[ 8
][ 4
]− 12
[ 8
]
= + − + −
9+ −
=0= 16
[ 5
][ 4
]− 7
[ 9 [ 4
]
+ − = + −
10 + −
=0= 19
[ 6
][ 4
] − 8[ 10
][ 4
]
+ − = + − + −
11+ −
= 033
= 15
[ 7
][ 4
]− 18
[ 11
][ 4
]− 17
[ 11
]
0 " " 2
The most important process variables are reaction temperature, acid strength,
isobutane concentration, and olefin space velocity. Changes in these variables affect
Reaction Temperature
The reaction temperature in sulfuric acid alkylation is usually in the range of 32!
o
50 F. At higher temperatures, oxidation reactions become important and acid
polymerization of the olefins becomes significant and yields are decreased. If the
operation takes place at lower temperatures, the effectiveness decreases due to the
increase in the acid viscosity and the decreased solubility of hydrocarbons in the acid
Acid Strength
effectiveness of the reactor mixing and the water content of the acid. In sulfuric acid
alkylation, the best quality and highest yields are obtained with acid strengths of 93!
95% by weight of acid, 1!2% of water and the remainder hydrocarbon diluents. The
water content in the acid lowers its catalytic activity by about 3!5 times as much as
hydrocarbon diluents, thus, an 88% acid containing 5% water is much less effective
catalyst than the same strength acid containing 2% water. At concentrations higher than
99%, isobutane reacts with SO3, and below 85!88% concentration, the catalyst becomes
higher acid strength necessary to keep acid dilution down. Increasing the acid strength
34
from 89% to 93% increases the alkylate quality by 1!2 octane numbers (Cupit et.al.,
1962).
Isobutane Concentration
of isobutane/olefin ratio. This is one of the most important process variables that
controls acid consumption, yield and quality of the alkylate. When the isobutane/olefin
ratio in the feed is high (15:1), the olefin is more likely to react with an isobutane
Thus undesired reactions are minimized. If this ratio is kept low (<5:1), production of
heavy polymers increases due to polymerization reactions, and the acid consumption
increases (Corma and Martinez, 1993). An increase in this ratio produces a decrease in
acid consumption while increasing both the yield and the quality of the alkylate. The
external isobutane to olefin ratio in sulfuric acid plants are usually in the range of 5:1 to
Olefin space velocity is defined as the volume of olefin charged per hour
! " $#
." * !
Lowering the olefin space velocity reduces the amount of high boiling hydrocarbons
produced, increase the product octane and lowers acid consumption (Gary and
Handwerk, 1984). Olefins space velocity is one way of expressing space!time, another
35
56 1
Olefins and isobutane are used as alkylation feedstocks. The chief sources of
olefins are catalytic cracking and coking operations. Butenes and propenes are the most
common olefins used but ethylene and pentenes are included in some cases. Olefins can
In addition to the alkylate stream, the products leaving the alkylation unit
include the propanes and normal butane that enter with the saturated and unsaturated
3. C5+ Alkylate
Only about 0.1 % by volume of olefin feed is converted into tar. This is not truly a tar
but a thick dark brown oil containing complex mixtures of conjugated cyclopentadienes
36
8
Concentrated sulfuric and hydroflouric acids are the only catalysts used
commercially today for the production of high octane alkylate gasoline but other
catalysts are used to produce ethylbenzene, cumene and long chain (C12 to C16)
The desirable reactions are the formation of C8 carbonium ions and the
olefins. Only strong acids can catalyze the alkylation reaction but weaker acids can
cause polymerization to take place. Therefore the acid strengths must be kept above 88
containing free SO3 also causes undesired side reactions and concentrations greater than
Isobutane is soluble in the acid phase only to the extent of about 0.1 % by
weight in sulfuric acid and about 3 % in hydroflouric acid. Olefins are more soluble in
the acid phase and a slight amount of polymerization of the olefins is desirable as the
polymerization products dissolve in the acid and increase the solubility of isobutane in
If the concentration of the acid becomes less than 88 %, some of the acid must
be removed and replaced with stronger acid. In hydroflouric acid units, the acid
removed is redistilled and the polymerization products removed as thick dark oil. The
concentrated HF is recycled in the unit and the net consumption is about 0.3 lb per
37
The sulfuric acid removed must be regenerated in a sulfuric acid plant which is
generally not part of the alkylation unit, and the acid consumption ranges from 18 to 30
lb per barrel of alkylate produced. Makeup acid is usually 99.3 % by weight H2SO4.
'
This chapter reviewed the literature for chemical process analysis and for the
Motiva’s Alkylation process and the results of applying the system to this process.
38
! " # $
" " %
&
& .
, /% -
' () *+ ) $ &
,
$ " !
! %
-
% " $
% $ $
"
3 "
! " 3 3 !
! ! " # $ !
3 4 ! 5"
3 !
6 " ' 3
3 /3 !
% 3 "
" ! '
! 6 ! !
! 4+& 5
12
! ! ! "
! ! 3
! '
&
" % #
( *
! % 7
" ! %
" 18
! !
% 9 3 "
% "
% " !
1
3 " "
! %
7 /% .
(: ! + * % "
! ! !
17
% ! !
% " % !
" # $
% ; * " #! #
( ) *
) ! " %
# $ ) # $
< !
% $ $ !
! 8" '
! "
$ !
"
% * + ! & 4+& 5
$ $
1
!
/ ! # $ ! # $
'
3 4 5
$ 3 !
! " ! "
" 3 ! "
" ! "
" !
! ! ! "
3 4 . " =>5
!" #
" ! ! !
! " ! ! ! ! !
! " ! ! ( ! !
! !< "
11
: ! !
& !
! #
$ ' ! "
< #? * : ! " !
! ! & " =
"
#3
! !
! ?
! ? & * "
3 !
"
18
" " ) * ' !< $
% $ "
" 3 !
" ) $ !
! % $ 4 !< 5 !<
!< !
$ ! # " $ @ !
!< $
! 4+& 5
) # $
! 9 "! !
! 9 !<
$ ! 8 : 6 4: !
6 5 3 . ! !
%;
1>
" " ,
, $
" %
" %
% 0 %
3 "
" !
"& . !
. " 3 "
3 " # 3 " !
# !
" % !
! "
# $ "
1A
.
" "
# ! ! !
* %; #! #
"
% "
" ! % !
- ." / 0
$ ! %
$ 6
3 !
3 !
! ! %
% "
, /% - $
% !
$ % !
! $
1=
! ! , /% -
4 !! ,/-5 !
,/-
"
% "
% $
+ %
$ " "
/ !
" !
% ! ,/-
$ !
"
% $
0 "
1
% ! 3
, ! " %
$ ,/-
% :
% " "
- : + %
1 + ! !
* %;
#! #
! ,
/% - !
"
$ !
,/- ! %
% " !
! ,/-
82
1
( % * ! B .
/ ! & 7
3 % !
/ !
% ) "
" % % "
# !
% %
! 3
" "
# #
# $
8
!
B . "
! ! "
% $
% " %
, ! 3 " %
6 !
" 3 /
% %
! % , "
! !
/ # $ ! !
% %
87
" % ! #
2+ ,
$ "
) #C ) $ % $
$ "
# & .
$ , /% - $
$ % "
! 3 ! $ 6
8
CHAPTER 4
MOTIVA ALKYLATION PROCESS
Advanced Process Analysis System was applied to a commercial alkylation
process at Motiva Refinery Complex in Convent, Louisiana. This chapter describes the
Motiva alkylation process and the simulation of the process developed with APAS.
Also, the validation of the process simulation is given, demonstrating that the
Alkylation Plant. The heart of the process is the STRATCO reactor or contactor, which
contacts the reactants in a high velocity propeller stream and removes heat from the
exothermic reaction.
(propylene, butylenes) are reacted with isobutane in the presence of sulfuric acid
catalyst to form hydrocarbon, mainly in the iC7 to iC8 range, called alkylate. The
blended with the refinery gasoline pool to increase the gasoline octane.
Figures 4.1, 4.2 and 4.3. These were prepared from P&ID’s of the plant. The reaction
section of the process is shown in Figure 4.1. The refrigeration, depropanizer and
deisobutanizer sections are shown in Figure 4.2. The saturate deisobutanizer section of
The process has four reactor pairs and four acid settlers as shown in Figure 4.1.
In the reaction section there are three feed streams, the olefin feed (HC01), the
54
isobutane feed (HC03) and the recycled olefin/isobutane mixture (HC32). The olefin
feed contains the light olefins that are reacted with isobutane in the alkylation unit's
STRATCO stirred reactors. The isobutane stream is in excess to fully react with all of
the olefins being charged to the unit. The process units are described in this section.
The two feed streams (HC01 and HC03) are cooled by exchanging heat with the
net Contactor effluent stream (5E-628, 5E-628, 5E-630) from the suction trap/flash
drum (5C-614) as shown in Figure 4.1. The two streams are then combined and fed to
pressure vessel containing a mixing impeller, an inner circulation tube and a tube
bundle to remove the heat generated by the alkylation reaction as shown in Figure 4.4.
The feed is injected into the suction side of the impeller inside the circulation tube. The
impeller rapidly disperses the hydrocarbon feed with acid catalyst to form an emulsion.
The emulsion is circulated by the impeller at high rates within the Contactor. The
sulfuric acid present within the reaction zone serves as a catalyst to the alkylation
reaction. The acid absorbs a small amount of hydrocarbon from side reactions and feed
55
1
5 M-3
Olefins Feed HC28 M-4
M-2 HC24
HC31 HC26
C301
HC27 HC25
C401
5C-614 Acid Settler Acid Settler
HC01 5C-631 5C-632
5E-628
STFD
HC30 AC15
AC05
C402 AC23
5E-629, HC32 AC12
630 HC02 Fresh Acid
M-24 M-1 M-7 M-11
AC02
3 HC03 HC04 5E-633 S-5 S-7
HC08 HC14
C403 HC06
AC07 AC18
S-2 HC07 AC09 HC11 AC20
3’ 5C-623 5C-625
2 HC29
Isobutane R1
R29
STRATCO
4 Reactor
R20
R3 R2 R7
R6
S-19 HC34 HC38
HC33 M-15
S-23
HC23
R10 HC22
56
Alkaline
Compressor 5E-634-639 ,
5 Water
C301 5E-655, 656 5E-641-644 Fuel Gas
C306
1 C313A
5K-601 C307 C310 C312
C312A
C303
C313 C315
Compressor KO Refrigerant
Drum C310 Accumulator
5C-615 5C-618 5C-658
5C-622
Economizer Depropanizer
Depropanizer
C302 Drum Charge
C308 Charge Drum
4 Caustic Settler
5C-616 C309
C312B C314
C314A
5E-640
C311 Discharge
Water
5
C315
5E-610-653 5E-613 To FCCU,
Fuel Gas
C325 6
C316 C325
C318
5C-604 7
5E-611 C329
Depropanizer
S-23
C329A C327 5E-621-624 C417
5C-603 IC4 from
C319 storage
5E-657, 658 C415
i-butane to
Storage C332
C317
Propane C414
S-34 5C-607
C320 C328 LPG
Alkylate
C321
Deisobutanizer
C331 5E-612
3
C322 M-33 5E-627, 647
Fresh Caustic C406 C412 C415
S-41
C404 n-butane
2 5C-606 C410
C403 C407 C411 C404 C418
5E-616 M-34
C411
S-42
5C-660 5C-608
C405A To IC4 Flush Pump C427 C406 5E-696
C405 C409
C409
C410 C419 5E-626
C405 5E-695
Discharge Mixer
C408
5E-617-620
C425
C407
C408 3’
Wash Water Wash Water
Alkylate
Figure 4.2. Refrigeration, Depropanizer and Deisobutanizer Sections of the Alkylation process
57
SC406
5E-605 - 608
5E-609A
SC404
SC413
5E-602
7
5E-603
n-butane to storage
SC403
SC405
58
In order to maintain the desired spent acid strength (89.0 wt% H2SO4) in the last two of
the eight Contactors, fresh acid (98.5 wt% H2SO4) is continuously charged to the first
two Contactors, and an equivalent amount of 89.0 wt% spent acid is withdrawn from
The hydrocarbon and acid are separated in four settlers 5C-631 to 5C-634 and
the acid-free hydrocarbon phase (HC22, HC23, HC25, HC27) flows from the top of the
acid settler through a backpressure control valve. Part of it goes into the tube side of the
contactor tube bundle (R1). The backpressure control valve is set at a pressure (60.0
psig) to maintain the contents of the settler in the liquid phase. As the hydrocarbon
stream passes through the control valve, its pressure is reduced to about 5.0 psig,
flashing a portion of the stream's lighter components, thereby, cooling the stream to a
temperature of about 30 °F. When the two-phase stream passes through the tube bundle,
reaction.
As shown in Figure 4.1, the second portion (HC29) of the hydrocarbon effluent
from the settler goes to the isobutane chiller (5E-633) and is used to reduce the
temperature of the isobutane feed stream. It then combines with the Contactor effluent
59
stream leaving the tube bundle and flows into the suction trap side of the suction
trap/flash drum (5C-614), where the liquid and vapor portions of this stream are
vapor space. The net Contactor effluent is accumulated on one side of a separation
baffle and is pumped to the effluent treating section. The cold refrigerant condensate is
accumulated on the other side of the baffle. This effluent recycle stream consists mostly
of isobutane and is returned to the Contactor by the effluent recycle pump. The vapor
portions of both of these streams combine and flow to the suction of the refrigeration
refrigerant vapor to allow condensing by cooling water. The vapor is condensed in two
from the propane-rich vapor in the refrigerant accumulator (5C-615), further cooled by
the refrigerant cooler, and then flashed across a control valve. Vapor and liquid are
stage of the compressor (5K-601) to be recompressed, while the liquid is returned to the
The remaining vapor portion of the compressor discharge stream (C310) flows
from the separator to the "total condenser" (5E-641-644). Liquid from the total
fed to the Depropanizer (5C-603) for stripping of propane before being sent to the
60
4.1.3 Depropanizer (5C-603)
The depropanizer feed stream is pumped from the charge drum to the
The overhead propane vapor from the depropanizer (C325) is totally condensed
by cooling water and collected in the depropanizer reflux accumulator drum (5C-604).
A part of the liquid is refluxed to the column's top tray (C329), and the rest is further
cooled and sent to the LPG section (C328). The depropanizer bottoms product (C317)
enters a water-cooled exchanger (5E-611) before it is split into three portions, one each
The net effluent from the suction trap/flash drum (5C-614) is heated in the
feed/effluent exchangers (5E-628; 5E-629, 5E-630) as shown in Figure 4.1 and then
treated with caustic to remove any traces of acid (5C-660 and 5C-608 in Figure 4.2).
The treated net effluent is fed to the alkylate deisobutanizer (DIB) tower (5C-606) to
strip isobutane from alkylate as shown in Figure 4.2. Two steam-fired thermosiphon
reboilers (5E-695, 5E-696) supply the majority of the heat to the lower portion of the
column by vaporizing liquid side draws (C411) from the column as well as the DIB
bottoms liquid (C408). A portion of the depropanizer bottoms is fed to the column as
and cooled by cooling water (5E-627, 5E-647). The normal butane product (C415) can
61
be used for vapor pressure blending or sold. The rate of the normal butane product will
vary depending on the desired deisobutanizer alkylate product Reid Vapor Pressure
(RVP).
in the isobutane accumulator pot (5C-607). Make-up isobutane (C417) from storage is
added to the pot and the combined stream (C418) is sent to the Contactors to react with
the olefins ( 3’ ).
deisobutanizer feed/effluent exchanger (5E-616) and also a heat exchanger with cooling
water (5E-617-5E-620). The cooled alkylate (C407) is removed as a product and sent to
storage.
Another source for isobutane is the saturate deisobutanizer (SatDIB) column (Figure
4.3), which strips isobutane from saturate feed coming from reformer unit of the
refinery (SC401). The saturate feed, containing mainly C4’s, flows to the SatDIB
feed/bottoms exchanger (5E-601) before entering the tower, as shown in Figure 4.3.
The overhead isobutane vapor is condensed by cooling water and a fraction of it (stream
SC413) is accumulated in the isobutane accumulator pot (5C-607, Figure 4.2). The
the SatDIB feed/bottoms exchanger (5E-601) and also by a heat exchanger with cooling
This section describes the models of the various process equipments in the
alkylation process. The process flow diagram, as developed with Flowsim, the
62
flowsheet simulation tool of Advanced Process Analysis System, is shown in Figure
4.5.
The equipment in the plant can be categorized according to their functions as follows:
Reaction Zone:
• STRATCO Contactor
• Acid settler
Separation Zone:
• Depropanizer
• Alkylate Deisobutanizer
• Saturate Deisobutanizer
• Economizer
• Compressor
• Heat Exchangers
• Condensers
• Reboilers
Miscellaneous:
• Mixers
• Splitters
• Reflux Accumulators
The models are developed for each of these units and are given in this chapter.
63
The alkylation process involves mixtures of a variety of components. The
components and compositions of streams vary widely from one section to the other. For
record keeping, all the hydrocarbon streams in the plant are composed of the following
components:
2. C4= Butenes
3. iC4 iso-butane
4. nC4 normal-butane
5. iC5 iso-pentane
6. nC5 normal-pentane
7. iC6 iso-hexane
8. iC7 iso-heptane
9. iC8 iso-octane
All the sulfuric acid streams in the plant are composed of:
The streams carrying the reaction products from the contactors to the acid settlers
64
Figure 4.5. Process flow diagram, as developed with the Flowsheet Simulation tool of Advanced Process Analysis System
65
4.2.1 STRATCO Contactor (5C-623)
The sulfuric acid alkylation reaction mechanism for pure propylene and pure
butylene feeds were shown in Tables 2.2 and 2.3. These are based on the Schmering
carbonium ion mechanism with modification introduced to account for iC9 and iC10
formation. The material balance on reactants and the product formation equations are
shown in the same tables. The parameters that affect the performance of the contactor
are: AC07
1. Operating Temperature
R3 R2
4. Olefin Space Velocity
HC34
with heat exchange. The intense mixing provided by the propellers ensures that there is
The composition and temperature of the exit stream were the same as those inside the
reactor.
Referring to Figures 4.1 and 4.6, the material and energy balance equations for
Description
Material Balance
Overall:
FR 2 − FR 3 = 0 (4.3)
Component:
i
for the hydrocarbon components. Here, r is the reaction rate of component i, in terms
of moles of i produced per unit time, per unit volume of the acid catalyst (metric ton
moles/(m3⋅min)), and MW
i
is the molecular weight of component i. The quantity
V5aC 623 , shown in the material balance equations, is the volume of acid in the contactor,
Considering the degradation of the catalyst, the material balance equation for the acid
07 − FAC 09 − R5 C 623 = 0
i i a
FAC (4.6)
a
where, R5C 623 is the rate of degradation of the acid (metric ton/min). The acid
degradation rate can be related to the volumetric flowrate of butylenes into the contactor
67
R5aC 623 = 0.121QHC
C 4=
07 (4.7)
The total mass flow rates are calculated as the sum of the component flow rates.
∑ FXi − FX = 0 (4.8)
The subscripts to the process variables like mass flow rates, specific enthalpies and
stream temperatures indicate the stream names, and their superscripts indicate the
component. The subscripts to the process parameters indicate the process unit they
represent.
hydrocarbon phase and an acid phase. Since the reaction model assumes that all the
reactions occur in the acid phase, the concentrations shown in the reaction model
The olefin feed to the process is composed of about 45% lv butylenes (1-BUT +
C2B + T2B) and about 3.5% lv propylene among other paraffins. Because of the low
the reaction mechanism of alkylation of isobutane with pure butylene (Table 2.3) was
adopted. The reaction rate terms in Equation 4.5 is evaluated using corresponding
equations from Table 2.3. Appendix A shows the procedure followed for calculating the
obtain the desired yields. Higher operating temperatures result in oxidation and
68
effectiveness due to the increase in acid viscosity and decrease in solubility of
maintained by recycling the cold refrigerant condensate from the Suction Trap Flash
Drum, which consists mostly of isobutane, and by circulating a part of the flashed, acid-
free effluent from the settlers, through the contactor tube bundle. When the two-phase
The tube bundle in the contactor, is effectively, a partial boiler, which vaporizes
a portion of the coolant. This can be modeled by assuming that the two-phases in the
coolant mixture are in equilibrium with each other throughout the tube bundle. The
component mass balance for the two-phase coolant stream can be written as
Ki = ∑ xi = 1 , ∑ y i = 1
yi
where, (4.10)
xi
In the above equation, y i is the mole fraction of component i in the vapor phase, and
The model can be simplified by assuming that all the vapor generation in the
tube bundle, due to transfer of heat from the contactor, occurs at a point outside the tube
bundle; i.e. after the coolant has passed through it. This is the same as considering the
FRi 3v FRi 3l
yi = xi =
∑ ∑
i
MW and MW i (4.11)
FRi 3v FRi 3l
MW i MW i
69
Energy Balance
Overall:
where, h X is the enthalpy of stream X in (MJ/min), and H 5C 623 is the heat generated in
The stream enthalpies are calculated from the component specific enthalpies, as
i
where, h X is the specific enthalpy of the component i in stream X in (MJ/metric ton).
Heat Transfer:
where, U 5C 623 is the overall heat transfer coefficient (MJ/m2⋅°C⋅min), A5C 623 is the
area of heat transfer (m2), and ΔTlm is the log mean temperature difference in the
(T5C 623 − TR 3 )
ln 5C 623
70
The mass and energy balance equations for contactor 5C-623 are summarized in Table
4.1.
Material Balances
FR 2 − FR 3 = 0
07 − FAC 09 − R5 C 623 = 0
i' i' a
FAC
Energy Balances
The equations given in Table 4.1 describe each of the four pairs of contactors.
However, the kinetic model is not a function of the catalyst concentration, and plant
data are not available for sulfuric acid concentration in the contactors. Consequently,
71
4.2.2 Acid Settler (5C-631)
The acid settler separates the emulsion from the contactor into the acid and
hydrocarbon phases. One acid settler is provided to each pair of contactors. The inputs
to the settler are the effluents from the contactor pair and the outputs are the separated
The material and energy balance equations for the acid settler (5C-631) are as
shown in Table 4.2. The mass flow rate of each component coming in is equal to the
mass flow rate of each component leaving of the settler. Since there is no heat
transferred from or into the settler, the temperatures of the various input and output
streams can be assumed to be equal. We also assume that the separation of the two
phases in the settler is complete, so that there is no acid component in the output
Material Balances
09 − FHC 27 = 0
i i
FAC i=1,2,3,4,5,7
09 + F AC 09 − F AC 05 − F AC 12 = 0
11 12
FAC
05 − s C 631 F AC 09 = 0
11 11
FAC
05 − s C 631 F AC 09 = 0
12 12
FAC
AC 05 − x AC 12 = 0
x 11 11
Energy Balances
T AC 09 = THC 27 ; T AC 09 = T AC 12 ; T AC 09 = T AC 05
72
4.2.3 Depropanizer (5C-603)
1963; Perry et al., 1984). This method can be applied to absorption, extraction
processes and distillation. The equation, which applies to the distillation process, is for
(1 + S ) + R(1 − S )
each component i,
(1 + S ) + R(1 − S ) + hS (1 − S )
N −M
f = N −M
n,i n
N −M M +1
(4.19)
n,i n n m
leaving in the bottoms product. The quantity S is the stripping factor (for each
component) and is defined for each group of stages in the column by Sn,i=KiV/L and
Sm,i=K'iV'/L'. K is the equilibrium ratio y/x for each component and V, L are the vapor
and liquid molar flow rates in the column. The Ki, V and L values are the effective
values for the top column section; K'i, V' and L' are the effective values for the section
below the feed stage. The quantity h depends on whether K or K' is used for the feed
stage. If the feed is mostly liquid, the feed stage is grouped with the lower stages and
K 'i L ⎛ 1 − S n ⎞
hi = ⎜ ⎟
K i L' ⎜⎝ 1 − S m ⎟⎠i
(4.20)
L ⎛ 1 − Sn ⎞
hi = ⎜ ⎟
L' ⎜⎝ 1 − S m ⎟⎠i
(4.21)
3. The reflux rate or the maximum vapor rate at some point in the column
73
4. The distillate rate D
Specification of D and L (or V), along with specifications of the feed stream
variables fixes the assumed phase rates in both sections of the column. Determination of
ideal solutions are assumed, K values are functions of only temperature and column
specified pressure. Estimation of K values and, in turn, Sn and Sm values for each
group of stages.
Tn =
(TN + TM +1 )
(4.22)
2
(
Tm = TM +1 + T1
) (4.23)
2
The feed enters the column at the (M+1)th stage. The material and energy balance
equations for the depropanizer are shown in Table 4.3. The Saturate Deisobutanizer
74
Table 4.3. Summary of the depropanizer model
Material Balances
where i=1,3,4,5
K i = f ( PC 603 , Tn )
K i' = f ( PC 603 , Tm )
S n ,i =
K i FC 325
FC 329
S m ,i =
K i'V '
(1 + S ) + R(1 − S )
L'
(1 + S ) + R(1 − S ) + h S (1 − S )
N −M
fi =
n ,i n ,i
N −M N −M M +1
n ,i n ,i i n ,i m ,i
fi = i
FCi 317
FC 316
K 'i L ⎛ 1 − S n ⎞
hi = ⎜ ⎟
K i L' ⎜⎝ 1 − S m ⎟⎠i
where i=1,3,4,5,7
TC 317 − TC 323 = 0
75
4.2.4 Alkylate Deisobutanizer (5C-606)
The alkylate deisobutanizer column has two feed streams, two product streams,
two side streams and two reboilers, one of which is used as a side reboiler. One of the
feed streams is fed to the column, at the top tray, in place of the reflux. There was no
readily available short cut method to model such a distillation column as it is, therefore,
the column was divided into three reasonable sections (C-606A, C-606C and C-606D)
which can be modeled using Smith Brinkley Method and by employing equilibrium
column with recycle and reboiler streams, however the model of deisobutanizer requires
the derivation of the method for the rectifying and stripping sections separately.
The equations for the alkylation deisobutanizer are given in Table 4.4.
76
Table 4.4. Summary of the deisobutanizer model
Material Balances
FC 427 − FC 431 = 0
FC 425 − FC 430 = 0
Species FCi 404 + FCi 432 + FCi 322 − FCi 414 − FCi 430 = 0
where i=1,3,4,5
S m ,i ,C 606 A =
K i',C 606 AVC' 606 A
L'C 606 A
S n ,i ,C 606 A =
K i ,C 606 A FC 414
FC 322
S m ,i ,C 606 D =
K i ,C 428 FC 428
FC 426
77
Table 4.4. Summary of the deisobutanizer model (Cont’d)
(1 − S N −M
)+ q ( S nN,i−,CM606 A − S n ,i ,C 606 A ) +
( )
S
f i ,C 606 A =
(1 − S N −M
n ,i ,C 606 A )+ h i ,C 606 A S N −M
n ,i ,C 606 A (1 − S M +1
m ,i ,C 606 A )
f i ,C 606 A =
FCi 430
( FCi 404 + FCi 322 + FCi 432 )
⎛ 1 − S n ,C 606 A ⎞
hi ,C 606 A = ⎜ ⎟
⎜1− S ⎟
FC 322
L'C 606 A ⎝ ⎠i
( ) ( )
m , C 606 A
where i=1,3,4,5,7
TC 431 − TC 425 = 0
78
4.2.5 Suction Trap/Flash Drum (5C-614)
The suction trap/flash drum (STFD) splits the hydrocarbon effluent from the
The net contactor effluent is accumulated on one side of the baffle and the cold
refrigerant condensate on the other. The vapors from the two compartments combine,
C401 HC32
combining in a mixer, as shown in Figure
Suction Trap Flash Drum
4.7. The suction trap acts as the separating
Figure 4.7: STFD 5C-614
vessel for the two-phase contactor from the
contactor tube bundles and from the isobutane chiller (5E-633). The liquid and vapor
streams leave the vessel as two separate streams. In the Flash drum, the cold condensate
from the economizer is flashed and separated. This unit is modeled as an equilibrium
stage with adiabatic flash. The material and energy balance equations for the two units
are as follows:
Suction Trap:
Overall material balance: (FHC 31l + FHC 31v ) − (FC 401 + FC 301 ) = 0 (4.24)
31l − FC 401 = 0
i i
FHC (4.27)
79
Energy balance: THC 31 = TC 301 = TC 401 (4.28)
Flash Drum
Ki = ∑ xi = 1 , ∑ y i = 1
yi
, where, (4.31)
xi
The flash operation has been described, in detail, in the contactor modeling section. The
Suction Trap Flash Drum model described above is summarized in Table 4.5.
80
Table 4.5. Summary of the suction trap/flash drum model
Material Balances
Species (Fi
HC 31l + FHC
i
31v ) − (FC 401 + FC 301 ) = 0
i i
31l − FC 401 = 0
i i
FHC
Ki = ∑ xi = 1 , ∑ y i = 1
yi
, where,
xi
Energy Balances
81
4.2.6 Economizer (5C-616)
flashed. This unit is modeled the same way as the flash C311
drum of STFD, i.e. assuming an adiabatic flash Figure 4.8. Economizer 5C-616
operation, where the two phases generated are in equilibrium with each other. The mass
Ki = ∑ xi = 1 , ∑ y i = 1
yi
, where, (4.35)
xi
In the above equation, K i is the distribution coefficient or the K-factor, y i is the mole
fraction of component i in the vapor phase, and xi is that in the liquid phase.
Enthalpy balance:
The combined vapor stream from Suction Trap and the Flash Drum is rich in
propane, and is sent to the depropanizer section. The vapor is compressed in the
82
Refrigeration compressor (5K-601), to allow condensing by cooling water. The
compressor was assumed to operate under adiabatic conditions, and hence the
⎡ ( k −1)
⎤
⎢ ⎛ ⎞
= RT1 ⎜⎜ ⎟⎟ − 1⎥
⎢⎝ p1 ⎠ ⎥
k p k
k −1
2
H ad (4.37)
⎣⎢ ⎥⎦
The work expended on the gas during compression is equal to the product of the
adiabatic head and the weight of the gas handled. Therefore, the adiabatic power is as
follows:
⎡ ( k −1)
⎤
k WRT1 ⎢⎛ p2 ⎞ k
hpad = = ⎜⎜ ⎟⎟ − 1⎥
⎢ ⎥
WH ad
k − 1 550 ⎝ p1 ⎠
(4.38)
⎢⎣ ⎥⎦
550
⎡ ( k −1)
⎤
⎢ ⎛ ⎞
hpad = 4.36 × 10 −3
Q1 p1 ⎜⎜ ⎟⎟ − 1⎥
⎢⎝ p1 ⎠ ⎥
k p k
k −1
2
or (4.39)
⎣⎢ ⎦⎥
where hp ad = power, hp; W = mass flow, lb/s; and Q1 = volume rate of gas flow,
ft3/min.
83
( k −1)
⎛p ⎞
T2 = T1 ⎜⎜ 2 ⎟⎟
k
⎝ p1 ⎠
(4.40)
C303 C304
C310
The constraint equations, for the compressor, are developed using the above formulas.
Enthalpy balance:
( k −1)
⎛p ⎞
T2 = T1 ⎜⎜ 2 ⎟⎟
k
⎝ p1 ⎠
(4.44)
where hp ad is the compressor as given in Equations (4.38, 4.39). The constraints of the
84
Material Balances
Energy Balances
( k −1)
⎛p ⎞
T2 = T1 ⎜⎜ 2 ⎟⎟
k
⎝ p1 ⎠
⎡ ( k −1)
⎤
⎛ ⎞
⎢⎜ 2 ⎟ k − 1⎥
=
⎢⎜⎝ p1 ⎟⎠ ⎥
k 8.314 FC 306TC 303 p
k −1
hp ad
⎢⎣ ⎥⎦
55
The olefin feed effluent exchanger cools the olefin feed by exchanging heat with
the alkylate effluent stream from the suction trap/flash drum. The material and energy
balance equations for this exchanger are shown in Table 4.7. Since no mass transfer is
involved in this unit, the masses of the two streams are conserved. The heat transferred
between the streams can be accounted for as shown in the energy balance equations.
85
Material Balances
FC 401 + FC 402 = 0
01 + x HC 02 = 0
i i
Species x HC
where i= 1,2,3,4,5
Energy Balances
86
4.3 Model Validation
The accuracy of the alkylation process model to predict the performance of the
plant has to be established. This model validation was performed by solving the data
validation problem using a set of plant data (steady state operation point #1) and then
comparing the results of the model with the plant data. For the 125 measured plant
variables, 88 were within the accuracy of the measurements. The remaining 37 variables
are shown in Table 4.8 for the plant data from the DCS, their reconciled value and the
standard measurement error ∈i (∈i = |yi -xi |/σi, where is σi the standard deviation of the
process variable). In consulting with the alkylation plant process engineers, it was
concluded that these 37 variables were within the range of possible process values. This
includes FC316, FC322, x5C417, xx5C412, xx7C414. The discrepancy in the variable
xx1C322 is a result of the model assumption that the propane purge from the drum 5C-
618 occurs infrequently. In summary, the model of the process accurately predicts its
(∈i)
(xi) Error
87
x1C417 0.0372 0.0255 3.1309
x2SC402 0.0136 0.0084 3.7929
x2SC408 0.0221 0.0002 9.9048
x3C325 0.0017 0.0000 10.0000
x3SC403 0.0103 0.0212 10.5665
x4C316 0.0580 0.0796 3.7155
x4SC408 0.0331 0.0088 7.3475
x5C316 0.0020 0.0060 19.8000
x5C417 0.0009 0.0295 286.2300
x5HC32 0.0096 0.0306 22.0134
x6SC402 0.0167 0.0666 29.8204
x6SC403 0.0250 0.0950 27.9946
x7HC32 0.0197 0.0497 15.2312
x7SC402 0.0022 0.0032 4.3956
x7SC408 0.0022 0.0000 10.0000
xx1C322 0.0027 0.1167 428.5338
xx1C414 0.0330 0.0800 14.2498
xx2HC01 0.4525 0.1291 7.1481
xx3C407 0.0003 0.0000 7.4194
xx3HC01 0.3558 0.0125 9.6498
xx4C407 0.1124 0.0853 2.4068
xx5C407 0.0803 0.1506 8.7555
xx5C412 0.0022 0.0581 255.6751
xx5C414 0.0021 0.0011 4.8325
xx7C414 0.0015 0.0080 44.4218
4.4 Summary
In this chapter, Motiva alkylation process and the system model for the process
was described. Also the validity of the model was demonstrated by comparing plant
data with model predictions. The complete model equations can be found in Appendix
B. The next chapter describes the results from the Advanced Process Analysis System.
88
! "
# $
% &' % &
( $
) # $ $
$ *
+ , $
,-
/ 0$
$ $ , $
$ $
$ $ 1 $
$ +
$ *
$ $ # $
$ # "2 * $
$ # "3 "4 # *
0 ,- * *
& 1 6 2 7 48
! 4 $ 9 * *
+* $ 6 8 *
1 *
$ +*
$ ! " # , "3 $ *
2!&:74 !
2!&:74
* *
( 0 ,
; < 3 2= * 25 *
6 , .8 =: 335 $
, # $
5
"# $% &
! 4 $
0 #
$ * 1
1 $
+, $
! 4 % &' %
! )0
,1 , 6:8 $
,1 $
# $ 6 , ! 48 372
$ 375 334
# *
* ? 6?
% $
3
# $
,!4
' $% &
# $ $ $
* %
3 25
1 # $
,!"
! 4 * $
$ $ *
* # $
$ 0
' )
, , $ 62!&:53 2!&:548
7
# $
,!2
! * $ # ,
, $ 1 # $
/ 0
! # , ,
/ A 0 $
/B +, 0 0
,!
) *+ &
# $
% &' %
9
7 4 " 2
+ C 5 D 3E D 7E D 4E D "E D 2E
$ 5 2
0 +
* # $
!: !=
4
, +
# $
23 22
23
# F
@ =:
335
+* ! 32=
1 * ! 25
G 372
@ G 325
:"
- . -* /
! 7 4 % &' %
# $
$ $ # 37
0 $
! &
! 7 4 % &' %
& $
"
? 0
# ( # $ H
!%I% 7 $ ?
= 33 7 J3 ,
27 J3 # 2K
4 24 6∈ C L &, L>σ 8 4 24
@ 43
$ 27 ,)3 )7
, B3
? ) G
27 G J3
L &, L>σ
! ! ! " #
$ ! ! !
# # ! ! !
#" #%$ ! ! !
" " ! ! ! #%"!
# % $$ # $& # ! ! # # "
# ! ! % & % % % &#
# "!% " $ "$ "$&
# "" & !%" !"" ! $ %
# % $" " "& #
# $ && " &!# &! &"$
2
" # # %&!! #!&& # # %#" " $&
" & & $#"% $ $ "
" # " " " " %
" & &$$ "#$ & " &"%
' %$% % !$ % #
' # % $! $ !$ %" %!$
" # " %& "$ ! "$"#
" ##"" # "% #" #"!#
" & %& # !& #"$$ # "#!
" # "" "#$ "!" "%
! " %%% %%$ %%$ # &! &
# #! %$$! !#%
# !" ##$ ! !" #!! !" & %
! ! %" ! ! !
! # !# #& !$ #&# !$" # % !$"
(' & #"# $"!$ &#$ %%"
(' %"$! &#$ &&! &#$
(' " ! &#$ $ "$ &#$
(' ! #"" &#$ &#$ &#$
(' #" $ "!" $# " $
(' #% & # &! &$
(' " %#"" %!$$ %"! %"%$
(' " %#$ % % % %&
) $ %" % % % !!!$ % & #
) %" &%# % % !!!$ % & #
) # %" &" ! % % % "&$#
) # % $! % % % #" % " !&
) #" %" $ % % #" % " !&
) " %& &## % %# #&%# %" $"
) " %" $#& % %# #&%# %" $"
) # # !" !"" % !%! % !$%& % "$$!
) # ! #" $ &% #"$ " & #"$ &%#$ #"$ %%$&
) # & # % ! # % " ## # $ #! % # %$ $
) # % # %"! # % " ## # $ #! % # %$ $
) # # % #&% # & #! # % #&% # %! !
) # ! #"" & #" #" " #" %
) # & #!# " ! # $ !$ # $ # $
) # # " # # $"% # $#
) # " #! # $ !$ # $ # $
) # # $ !! & # $$& # !%! # !% $
) " " # % !& # "$ # ! %% # & $"%
) " " ! " " ## "
) " & # $ # # & # # $ # $
) " % " ! " " "
) " # ! " #!# # ! #!# #$ & #!# #!!
) " " # "" ! ##& % $ ##! &% " ##! &&&"
:
) " % # $ % # & "## # & "%!& # % "#"
) " $ # !#!& # " " # % # $& &
)' # !# "# ! $ && ! $&% ! ""
) " # " !$%$ # %%%! # "! ! # " ! !&
) " # ##& $ % ##! ! "! ## &" " ## &"$
) " # ! # # ! # !
) " % # $ $ # % &"" # & $ &! # &$ %
) " # $% !$ & # # #
* $!& $!%& $& % $!$
* # $#% $" " $" $"#
* #" $ $ ! $ ! $ !$
* " $ ! %$ %$ %$ !
* # ! " % !! %" &%
* # $&
* " & #& #
* ' # % % # ##
* " "# !# " "
* " #
* " % "& "&& "&"
* " #! %" %" # &$ $
* " # $ $ $
* " % $ $ "%
*# # ! &&$ &% &%% &%%!
*# # &
*# " & %!$$ $&#% % #& % "
*#' # % %& &! " && % && %
*# " $ # ! $ $ $&#
*# " # # #" !!
*# " % %&! $&#$ $"# $"#!
*" # ! % %$" &$ &$! #&
*" " & ! $ ! " $
*"' # " " &
*" " $#" % ! %&
*" " # % " &!!! &$" &$#!
*" " % ## ! %! %% & #"&
* # ! ! ! ! $%
* " & $ $ $ %! #
* ' # $! %% # ! # ! #"
* " # !%! ! !
* " # & % &#" &#!
* " %
*! " !& !%! !!! !!! $% "
*! " # $"$ $ & $$"!
*! " %
*&' # $& "&$ "$! "$& #
*& " !$ # # " #$ !
=
*& " # "! "! "!
*& " %
** # & ## !& !& " % ##%
** " " ## &! &$% % " "$%
** ' %% $&! $$
** ' " $! $ $ & "%
**# # & $"# &$ &$# &$#
**# # $"# &$ &$# &$#
**# " & # & " $"
**# " & & &
**# " " %&& % % % $! % $!
**#' # % $ !"$%
**" # & !& %$$ % %
**" # !& %$$ % %
**" " & " % % % #
**" " $ " %&& %!!% % "$
**" " " %! !& $ ! $ "
**"' % % !" !&
** " & % # " ! %&
** " !% % % !&
** " " " %#
**& " " % % % "" " %
334 3" 5 J3 ,
$ 27
24 4: :" $
H $
$ ,B7
? +
24 J3
+ , ,
!#" & &
!"
!
! #
! $
! $#"
! %
! !
! &
! & &#$&
! -
! !
! & !
! &- &
! %
! $
!#"
!" & !$
!"
. ! " "
. ! " "
. !$ $ $
. !$! " "
! ! $ $
! ! %$ %$
! ! $ $
!## " "
/ !
/ ! # $$ "
/ ! !
! $ " $
! # " "
0 !# $&& $&&"
0 !# $% $%
0 !## $$ $$ #
0 !#" $$ $%$#
) 1 $ $
2 ! % %
2 ! ! !
2 ! #
2 ! " "$%
2 ! $ " "
2 ! %# %$%
2 ! $$
2 ! # #
2 ! # !
2 ! !
2 ! & !
2 ! " "
2 ! - &% & $
2 ! !
2 ! &
2 ! &-
2 ! % ""
2 ! $
2 !## ! &&
2 !#" !
2 !"
2 !" %" %&"
2 !$ ## # $
2 !$ - #$ #%%
2 !$!
2 !$!-
2 !33 # %!
" -* /
, $ H
$ $
2" H $ $ $
$ $9
C &! &@
! CΣ1 E!
@ CΣ1 E@ !
355
$ 1 ! 6 8 !" 0 6
) 8 1 @
$ $ J3 ,1
:4 $ M7 33>
$ 3""K 6J38
2K 7 7K
2 2K 355K
6 $ 8 = 7K
) 7 7K *
M2:55 >
, $ 22$ $ 72K
73:K $ ,
$ , $
455K "555K $
)3 )7
,B4 ?
! $ &
353
2" ( $ >@ !
# ! 6M> 8
I N
$
B!53 33 = 35 ""
. B!53 3 55 3: 2:
1 & !"3" 3: 3= 4
$
I& !"52 !"34 34 7 37 =3
!4 !"5= 7" " 77 45
!" !"5= 7: 47 7" 5:
!" (
!"5= 7: 4" 7" 3
! @ !
357
22 ! ) G 6) G 8 + 6 +8 +
% 6+ % 8 ,) % 6 8
% % - 12
J3 33 37 3 7 3 3""
J4 73 " 77 3 7: 7:
J" =5 =5 77 3 73:
J2 35 3 74 4 7: 4 3:5
J: 77 5 74 : 7= : 72
K 37=
" 3 4 5 -* /
B +, I $ % 6 B+I8 $
, 1
$ &
* $ & ,
, *
354
7 , , $
$ $
! 4 B+I
,+
$ $ 0
, - +
0 # !: !=
* *
& * $9
7 4 " 2
BC 5 D 3 D 7 D 4 D " D 2 -> 62 38
CΣ 1 ,1 62 78
$ 1 ,1 1
! 4
$ 0 #
$ 34 O 6 0 8
35"
$ $
#%( ( I B+I
% ? ! ! 6?!!8 $ # 23
?!! $ $ * $
$ H
, *
* 3="7 ->
, "5"4 -> , #
* 3K 3 -> $
= =K * = "K
"4:= ->
3 :K * "455 ->
# 23? ! !
352
, + ( 6 +(8
$ $ ,
?!! # 27 $
3: , " 32 $ : , "
37 *
B !"5: 6 )
6) 8 !"37 6 ) 8
, , 6 +&:3:8 B $
* & 6) 3 8 2!&
:53 6 = O8 *
35:
# , $
, $
1 ,
# 27 I $ ? )
35=
" "
# #
!
!
# #
% %
" % ! " % !
# 24 1 62!&:53 2!&:548 $ 9
2!&:53 6 8 6 8
35
'
$ , 1
$ 0
,-
35 2:
# $ & #
% &' % 0
2: 1 %
) 1 ,
!57 # # 1 5 5
B!53 % # 1 3 :77
!"3" & 1 1 3 :33
!"53 &) # 1 3=
!"2 I & 3 54"
!475 ' ? 5
!47 # ? 5
!"5= 5
!"34 I& 5
!"52 I& 5
35
! 4
6 +1 8 +1
,? ( $ $
@ +1 $
&
$ 2= 6.+%8 6 +%8
&
2= G 6.+%8 6 +%8
1 , G
6.+%8 6 +%8
&" 375 &" = :: >
&4 7 :5 &4 "2 " >
&5 === &5 ="7 >
3 5472 3 544= >
5: = 5 ="24 >
5 35: 5 332" >
0 $
0 6 7 7K8
335
+
$ # $ % &' % $
333
!
"
# "
$% &
"
! ,
' +. ' +. , /
" 0 ( 1 ++ %.
23 4 " 5 '.
. ' '.
&&. " $ .
6 " ! .
) 7
( +, &( 8 #1 & %.
0'%&& 1 "
! .
, ) (&$ 9:1
+*&& 9:1 ! )
) " "
*
REFERENCES
Albright, L. F., 1990a, Modern Alkylate - 1:Alkylation will be key process in
reformulated gasoline era, Oil and Gas Journal, 88, no. 46, 79-92
Brook, A., D. Kendrick, A. Meeraus, 1998, GAMS- A User’s Guide, Release 2.25, The
Scientific Press
Cabezas, H., J. C. Bare and S. K. Mallick, 1997, Pollution Prevention with Chemical
Process Simulators: The Generalized Waste Reduction Algorithm, Computers
Chem. Engng., Vol. 21, Supp. p. S305-S310
Chen, X., 1998, The Optimal Implementation of On-Line Optimization for Chemical
and Refinery Processes, Doctoral Dissertation, Louisiana State University,
Baton Rouge, LA
Corma, A. and A. Martinez, 1993, Chemistry, Catalysts and Processes for Isoparaffin -
Olefin Alkyalation: Actual situation and Future trends, Catalysis Reviews:
Science and Engineering, 35, no. 4, 485-523
Gary, H. James and H. E. Glenn, 1984, Petroleum Refining, Second Edition, Marcel
Dekker, Inc.
Graves, D. C., 1999, Acid consumption and acid runaway on a sulfuric acid alkylation
unit, Alkylation Seminar in Scottsdale, Arizona.
114
Gruse, W. A. and D. R. Stevens, 1960, Chemical Technology of Petroleum, Third
Edition, McGraw-Hill, New York
Gundersen, T. and Naess, L., 1988, The Synthesis Of Cost Optimal Heat Exchanger
Networks, Computers Chem. Engng., Vol. 12, No. 6, p. 503-530
Hegde, S. C., 1999, Simulations for Pollution Prevention, Master’s Thesis, Lamar
University, Beaumont, TX
Heijungs, R. (Final Ed.); Guinee, J. B., Huppes, G., Lankreijer, R. M., Udo de Haes, H.
A. and Sleeswijk, A. W., 1992, Environmental Life Cycle Assessment of
Products Guide - October 1992, Centre of Environmental Science, Leiden
Hilaly, A. K. and Sikdar, S. K., 1994, Pollution Balance: A New Methodology for
Minimizing Waste Production in Manufacturing Processes, J. Air and Waste
Manage. Assoc., Vol. 44, p. 1303-1308
Hohmann, E. C., 1971, Optimum Networks for Heat Exchange, Ph.D. Dissertation,
University of Southern California.
Langley, J. R., 1969, The Kinetics of Alkylation of Isobutane with Propylene Using
Sulfuric Acid Catalyst, Doctoral Dissertation, Louisiana State University, Baton
Rouge, LA
Lee, K.F., A. H. Masso and D. F. Rudd, 1970, Branch and Bound Synthesis of
Integrated Process Design, Ind. Engng. Chem. Fundam., Vol. 9, p. 48-58
Masso, A. H. and D. F. Rudd, 1969, The Synthesis of System Designs. II. Heuristic
Structuring, AIChE Journal, Vol. 15, p. 10-17
Perry, R. H., Green, D., Malony, J. O. (editors), 1984, Perry’s Chemical Engineers’
Handbook, Sixth Edition, McGraw-Hill International Publishers
115
Rosenwald, R. H., 1978, Alkylation, Kirk-Othmer Encyclopedia of Chemical
Technology, Third Edition, John Wiley & Sons, New York
Shaney, Agnes, 1995, Cleaner by Design, Chemical Engineering, Vol. 102, No. 7, p.
32.
Telang, K., 1998, Advanced Process Analysis System, Master’s Thesis, Louisiana State
University, Baton Rouge, LA
Templeton, P. C. and B. H. King, 1956, paper presented before the Western Petroleum
Refiners Association, June 21-22
Thomas, C. L., 1970, Catalytic Processes and Proven Catalysts (McGraw Hill Book
Company, New York), pp. 87-96
Van Reeuwijk, P. A., R. Zwart and M. E. J. Dreissen, 1994, CAPE in the Development
of an Energy Efficient Chemical Process, Computers Chem. Engg., Vol. 18,
Suppl. P. S161-S164
Warren L. McCabe, Julian C. Smith, Peter Harriott, 1993, Unit Operations of Chemical
Engineering, Fifth Edition, McGraw-Hill, Inc.
Winter, Peter, 1992, Computer Aided Process Engineering: The Evolution Continues,
Chem. Engng. Prog., Vol. 88, No. 2, p. 76-83
Yongkeat, S, 1996, The sulfuric Alkylation of isobutane with light olefins, Master’s
Dissertation, Lamar University, Beaumont, TX
116
117