Flared gas analysis
The analysis of the composition of the flared gas is carried out by using gas chromatography (G.C.). In order to obtain the best result that can give us the composition of the flared gas, the selection of the samples was based on :
Following the analysis of the different gas samples, we obtained the results reported in Table 1.
The operating conditions of the two flares during normal train operation are summarized in Table 2.
Table 1
Composition of flared gas.
Components
|
Hot torch (% molar)
|
Cold torch (% molar)
|
CO2
|
0.102
|
0.980
|
Nitrogen
|
1.221
|
4.000
|
Methane
|
77.080
|
60.230
|
Ethane
|
9.621
|
16.140
|
Propane
|
5.450
|
8.810
|
n-Butane
|
2.930
|
4.430
|
i-Butane
|
2.335
|
4.110
|
n-Pentane
|
0.195
|
0.255
|
i-Pentane
|
0.960
|
0.329
|
n-Hexane
|
0.106
|
0.725
|
Table 2
Operating conditions of the two flares.
|
Hot torch
|
Cold torch
|
Flow rate Kg/h
|
15954.15
|
23798.4
|
Pressure Bar
|
1
|
1
|
Temperature oC
|
25
|
-70
|
Simulation of flare gas recovery system
The simulation of flare gas recovery methods is carried out using ASPEN-HYSYS V12 software. This software is extremely powerful for steady state simulation. In the present research, the Peng-Robinson thermodynamic model is used for the simulations, which is the most widely used model in the field of oil and gas (hydrocarbon). In this study we have simulated three methods of flared gas recovery at the GL1Z complex, namely electricity, LPG and GTL production (Fig. 5).
Production of electricity
Flared gases are used to generate electricity. This production is used for distribution or for internal consumption, generally using a turbine. The Brayton cycle is one of the most efficient cycles for converting gaseous fuels into mechanical energy and electricity using a turbine.
The role of the combustion chamber is to combust a mixture of air and flared gas coming from the compressor, where the pressure is increased from atmospheric pressure to 23 bar. The gases produced by combustion, where the potential chemical energy contained in the fuel is transformed into heat energy, are directed towards the turbine. The turbines of turbo engines are the site of an adiabatic detent which transforms the energy available in the active fluid into mechanical energy, which is converted to electricity using the generator. The exhausted gases then leave the gas turbine. Fig. 6 shows the Brayton cycle simulated in this study.
Figure 6. Production of electricity (Rahimpour et al. 2012).
Using the ASPEN HYSYS v12.0 simulator, a gas turbine cannot be simulated directly, but a system representing a gas turbine can be simulated. Figure 7 shows the simulation of the transformation of flared gases to produce electricity.
The combustion reactions that take place in the combustion chamber are grouped together in the following Table 3:
Table 3
Reactions
|
ΔHo298 (*10+ 6 KJ/Kmol)
|
CH4 + 2O2 → CO2 + 2H2O
|
-0.8
|
C2H6 + 3.5O2→ 2CO2 + 3H2O
|
-1.4
|
C3H8 + 5O2→ 3CO2 + 4H2O
|
-2.0
|
C4H10 + 6.5O2→ 4CO2 + 5H2O
|
-2.6 for iC4
-2.7 for nC4
|
C5H12 + 8O2 → 5CO2 + 6H2O
|
-3.3
|
C6H14 + 9.5O2→ 6CO2 + 7H2O
|
-3.9
|
To ensure the complete combustion of the entire gas charge, a sufficient quantity of oxygen must be given for the complete combustion of each component of the charge. In the ASPEN-HYSYS simulator we have the "Set" function to establish a relationship between the gas molar flow rate QGas and the oxygen molar flow rate Q O2 :
𝑄O2 = 𝛼 ∙ 𝑄𝐺𝑎s (1)
To calculate the α coefficient, we need the stoichiometric coefficients Ki of each reaction, then multiplying by the molar fraction Xi of the components between the fuels, eliminating the non-combustible components (CO2 and N2). To do this we need to recalculate the X'i of New.
Xi: molar fraction of each compound in the gas charge.
Ki: stoichiometric coefficient of oxygen in the reaction with component i.
X’i: molar fraction of compound in combustible compounds.
$${X}_{i}^{{\prime }}=\frac{{X}_{i}}{{\sum }_{C1}^{C6}Xi}$$ 2
αi : partial coefficient of flow regulation.
𝛼𝑖 = 𝐾𝑖 ∙ 𝑋′𝑖 (3)
α : flow regulation coefficient.
$$\alpha =\sum _{1}^{6}{\alpha }_{i}$$ 4
⟹𝛼 ≥ 2.6
The relation (1) becomes :
𝑄𝑂2 = 2.6 𝑄g𝑎s (5)
The molar composition of air is 78.05% nitrogen N2, 20.95% oxygen O2 and 1% other gases.
The relation (5) becomes
Q air = 12.41 Qgas (6)
- For the compressor « gas comp » :
Consumed power E1 = 4297 kWH.
- For the compressor « Air comp » :
Consumed power E2 = 9.527 104 kWH.
- For the turbine « Elecrical Power » :
Power generated E3 = 2.590 105 kWH
Production of electricity E = E3-E2-E1 = 1.13 105 kWH
Production per year = 1.13 105*24*365 = 9.9108 kWH.
Simulation of LPG unit
One of the methods used to recover flared gas using a treatment plant is to produce:
Recovery of the LPG contained in the flared gas is a method of increasing the quantity of LPG produced (Fig. 8).
The system for transforming flared gases into LPG, condensate and Fuel gas comprises four steps:
- Demethanisation: recovery of methane gas.
- Deethanisation: recovery of ethane.
- Depropanisation: recovery of propane.
- Debutanisation: recovery of butane.
The mixture of gases at the outlet of the demethaniser and dethaniser is the Fuel gas. The LPG is the product of the mixture of gases at the top of the depropaniser and debutaniser. The condensate is recovered at the bottom of the debutaniser column. Figure 9 shows the simulation of the flared gas recovery unit in the form of LPG, condensate and Fuel gas.
Simulation of a process to recover 39752.55 Kg/h of flared gas results in :
27663.79 Kg/h Fuel gas.
10900.01 Kg/h LPG.
1188.75 Kg/h condensate.
The composition of the Fuel gas, LPG and condensate is given in Table 4.
Table 4
Compositions of finished products.
Components | Fuel gas (%molar) | LPG (%molar) | Condensate (%molar) |
---|
CO2 | 3.211 | 0.000 | 0.000 |
Nitrogen | 0.686 | 0.000 | 0.000 |
Methane | 78.323 | 0.000 | 0.000 |
Ethane | 15.289 | 0.431 | 0.000 |
Propane | 2.152 | 53.116 | 0.000 |
n-Butane | 0.133 | 28.729 | 0.999 |
i-Butane | 0.201 | 14.858 | 0.071 |
n-Pentane | 0.002 | 0.555 | 17.055 |
i-Pentane | 0.003 | 2.294 | 34.091 |
n-Hexane | 0.001 | 0.015 | 47.784 |
The molar composition of the LPG produced by the simulation is compared with the composition of the LPG produced at the GP1Z1 plant. This comparison is reported in Table 5.
Table 5
Comparison of LPG composition.
Components | LPG (GP1Z) % molar | LPG (simulation) % molar |
---|
Methane | 0.23 | 0.0 |
Ethane | 1.59 | 0.431 |
Propane | 59.30 | 59.118 |
i-Butane | 14.10 | 14.858 |
n-Butane | 24.54 | 24.729 |
i-Pentane | 0.20 | 0.294 |
n-Pentane | 0.03 | 0.555 |
n-Hexane | / | 0.015 |
H2O | 0.002 | / |
Simulation of GTL unit
GTL (Gas to Liquids), a technology that enables natural gas to be converted directly into various liquid synthetic petroleum products, has been the subject of growing interest in recent years from oil companies, governments, industry, universities and research and technology development centres. The aim is to use natural gas, which is made up of small molecules, to produce products with much larger molecules. The main component of natural gas is methane, whose molecule comprises 4 hydrogen atoms solidly arranged around a carbon atom, and is characterised by particular stability.
Its transformation into a liquid first requires the chemical bonds between atoms to be destroyed by a major external input of energy, in the form of heat and at high pressures. Carefully selected catalysts are involved in the chemical reaction, without actually being altered or consumed. The conventional "indirect" approach used to break these bonds is based on pure power; the chemical bonds of the methane molecule are ruptured using vapour, heat and a catalyst, resulting in the formation of a mixture of CO and hydrogen, generally known as "SynGas" or synthetic gas.
The next step is to convert the syngas into various liquid fuels using the Fisher Tropsch process, invented in 1923. Catalysts are essential to this process. The liquid fuels formed are refined to obtain a variety of products using conventional methods (Fig. 10).
The objective is to transform flared gases containing methane and ethane of a majority composition into components with longer molecular chains. The GTL process is very complex compared with other methods.
The GTL process is composed of different units, some of which are principal, such as the synthesis gas production unit, the FT (Fisher Tropsch) synthesis unit and the product valorization unit.
✓ Syngas production unit: the reaction section comprises two reactors. The Pre- reformer is feed by a charge with a flow ratio of H2O/CH4 = 2.3.
All hydrocarbons heavier than methane react completely in the first reformer. As a result, the feed stream to the secondary steam reformer contains no hydrocarbons heavier than methane.
The next Table 6 shows the various reactions:
Table 6
| Reactions | ΔH298 0 KJ/mol |
---|
Reactions of conversion | C2H6 + 2H2O→ 2CO + 5H2 | + 350 |
C3H8 + 3H2O → 3CO + 7H2 | + 500 |
iC4H10 + 4H2O→ 4CO + 9H2 | + 660 |
nC4H10 + 4H2O→ 4CO + 9H2 | + 650 |
iC5H12 + 5H2O→ 5CO + 11H2 | + 810 |
nC5H12 + 5H2O→ 5CO + 11H2 | + 800 |
nC6H14 + 6H2O→ 6CO + 13H2 | + 950 |
Equilibrium reactions | CO + 3H2 → CH4 + H2O | -210 |
CO + H2O → CO2 + H2 | -41 |
A third Auto Thermal Reformer Reactor (ATR) fed by a charge with an O2/CH4 = 0.5 ratio given by Otaraku and Vincent (2015). Table 7 shows the equations reactions and their corresponding reaction enthalpies.
Table 7
Auto Thermal Reformer Reactor (ATR) reactions.
Reactions | ΔH298 0 KJ/mol |
---|
CH4 + 1.5O2 → CO + 2H2O | −519 |
CH4 + H2O → CO + 3H2 | 210 |
CO + H2O → CO2 + H2 | -41 |
✓ The Fisher -Tropsch (FT) synthesis unit: the reaction section is charged with H2/CO = 2.1. The FT reactions are given in Table 8.
The stoichiometric coefficients for the FT reactions are based on the ASF distribution determined by (Kerron et al. 2014) and Otaraku & Vincent (2015) :
Xn = (1-ά)2άn-1 ά = 0.9 (7)
ά : chain growth probability, which is a direct measure of the probability of an FT catalyst to catalyse a chain.
Xn : mole fraction of each carbon number (n)
Table 8
Fisher -Tropsch reactions Otaraku and Vincent (2015).
Reactions | ΔH298 0 KJ/mol |
---|
\(CO+2.1{H_2} \to ~\sum\limits_{{n=1}}^{{20}} {{{(ASFcoefficient)}_n}{C_n}{H_{2n+2}}+{{(ASFcoefficient)}_{30}}{C_{30}}{H_{62}}+{H_2}O}\)
|
-160
|
CO + 3H2→CH4 + H2O
|
-210
|
The catalyst used is Cobalt, and the following reaction given by Otaraku and Vincent (2015) gives a better representation for the H2/CO = 2.1 ratio.
CO + 2.1H2 → 0.01CH4 + 0.009C2H6 + 0.008C3H8 + 0.007C4H10 + 0.007C5H12 + 0.006C6H14 + 0.005C7H16 + 0.005C8H18 + 0.004C9H20 + 0.004C10H22 + 0.003C11H24 + 0.003C12H26 + 0.003C13H28 + 0.003C14H30 + 0.002C15H32 + 0.002C16H34 + 0.002C17H36 + 0.002C18H38 + 0.002C19H40 + 0.001C20H42 + 0.012C30H62 + H2O. (8)
Fig .11 below shows the simulation of flared gas recovery using the GTL process
The simulation of a flared gas recovery process enables :