Problématique :

Le stockage de GPL à Guellala se fait dans des sphères sous une pression de (10.4 -
10.7) bar, et une température de (23 - 25°C), mais en été la température ambiante augmente,
elle peut aller parfois jusqu'à (55 - 60°C), ce qui provoque l’augmentation de la pression dans
la sphère, et par conséquent l’évaporation du GPL suivit de l’ouverture de la vanne installer
au sommet des sphères, on aura donc le torchage de GPL.
Le but de ce mémoire est d’estimation de la quantité du GPL qui se perd, et la proposition
d’une moyen de récupération (calcul de conception) ; nous allons répartir notre travail comme
suit :
-Estimation de la quantité du Boil-off.
-Elaboration d’un schéma d’une section de récupération de Boil-off.
-Calculons des principaux équipements du circuit.
Stockage et expédition :
Les trois principaux paramètres pour la bonne exploitation d’un réservoir de stockage
GPL sont les suivants :
-La température du stockage.
-La pression de la phase gaz du réservoir.
-Le niveau de produit contenu.
Gaz torchés :

Ce sont les gaz brûlés au niveau des torches avant d’être renvoyé vers l’atmosphère,
ces gaz proviennent de l’ensemble des unîtes de liquéfaction de gaz naturel du complexe
GLA.

Estimation du Boil-off :

Etant donné que la quantité du Boil-off est instable, ce qui explique la variation de la
quantité des gaz torchés, cette quantité qui varié d’un mois à un autre et d’une année à une
autre, ceci dépend de :
-Nombre d’unités mises en production ou service.
-Nombre de déclenchements mensuels et annuels.
-La quantité de GPL produite et de la température de coulée..
-Nombre de chargement méthaniers.
-Condition météorologiques…etc.
Notre travail tient en considération l’estimation du Boil-off de l’année 2004 (du 04 /08/2004).
L’estimation du Boil-off de l’année 2004 est basée sur :
-L’analyse de l’autoconsommation an 2004.
-La détermination de la quantité du gaz consommé dans les sphères année 2004.
-La détermination de la quantité des gaz torchés année 2004.

1
Analyse de l’autoconsommation : (consommation de l’usine plus pertes) :

Définition :

-On appelle autoconsommation de l’usine, la quantité de l’énergie PCS consommées

ou perdue durant le fonctionnement et rapportée à la quantité d’énergie PCS entrant de l’usine
de liquéfaction.
-PCS (Pouvoir Calorifique Supérieur), c’est la quantité de chaleur dégagée lors de la
combustion complète et sous pression constate d’une unité de volume de gaz. Les éléments de
la combustion étant pris initialement à une température de référence et les produits de
combustion ramenés à cette même température, l’eau produite par la combustion est
totalement condensée. La température de référence sera fixée d’un commun accord par le
client et le vendeur. Les valeurs courantes sont 0°C, 15°C, 15.6°C et 25°C. Les mesures
calorimétriques de bases sont généralement faites à 25°C. Les valeurs correspondant sont
obtenues par calculs.
-PCI (Pouvoir Calorifique Inférieur) ou LHV (Low Heating Value), sa définition est la
même que celle de PCS. Toutes fois, L’eau reste totalement à l’état de vapeur a la température
de référence.
-Normal mètre cube (Nm3), C’est la quantité de gaz contenue dans un mètre cube dans
les conditions normales, soit 0°C et 1.013 millibars absolue (1.013 millibars = 1 atm = 760
mmHg).
-Standard mètre cube (Cm3), Il est appelé CONTRAT METRE CUBE. Quantité de gaz
contenue dans un mètre cube dans les conditions standard 15°C et 1bar absolu.
1Cm3 = 0.93544 Nm3
-Tonne métrique (TM), Vaut 984.206 Kg par rapport à la tonne ordinaire qui vaut
1000Kg.
-Calorie (Cal), C’est la quantité de chaleur nécessaire pour élever la température d’un
gramme de matière de 1°C sous pression atmosphérique normale.

1Kilcalorie = 1000 calories

1 thermie (th) = 1000Kilcalories = 106 calories
1Kilcalorie = 3.968 BTU (British Thermal Unit)

Bilan année 2004

Quantification du Boil-off:

Le Boil-off est calculé comme étant 1.4% vaporisé de la quantité de GPL produite, au
niveau des sphères des stockages et pendant la coulée. Etant donnée que le phénomène du
Boil-off est inévitable par suite de l’échange thermique avec le milieu ambiante par les trois
mécanismes de transfert de chaleur (rayonnement, convection, conduction) ce coefficient de
1.4% représente le pourcentage évaporé de GPL, et il est donné par le constructeur lors de la
conception des réservoirs.
Comme la quantité de GPL produite est donnée en (m 3), on doit avoir son équivalent
en Cm3 en GN. Pour cela on doit faire appel à un rapport de masses volumiques.

2
L’équation de Boil-off

Q .24.100
T
% 
Boil  off jour  .V .H
G sph

Q   : Quantité de flux thermique totale Kcal/h.

T
G   : Masse volumique totale du mélange g/Cm3, Kg/m3.
.V
sph  : Chaleur de vaporisation de latente Kcal/Kg.
H : Volume de sphère m3.
Le tableau suivant représente la composition moyenne en pourcentage molaire du Boil-off
(04-08-2004) de l’année 2004.

Constituant X’Li (%)

C2 0.75
C3 35.02
iC4 12.61
nC4 51.11
iC5 0.02
nC5 0.49
Total 100

La masse moléculaire moyenne du Boil-off :

M  M i X i' .
D’où : X’ = le pourcentage molaire de chaque constituant.
Mi = Poids molaire de chaque constituant.

Le tableau suivant présente le calcul de la masse moléculaire du Boil-off

Constituant X’Li (%) Mi (g/mol) M i X i' .

3
C2 0.75 30 0.225
C3 35.02 44 15.4088
iC4 12.61 58 7.3138
nC4 51.11 58 29.6438
iC5 0.02 72 0.0144
nC5 0.49 72 0.3528
M  M i X i' . 53.0895

La masse volumique du Boil-off

On pourra toujours considérer qu’un gaz ou des vapeurs d’hydrocarbures sous une
pression inférieure à 03 atm, ont un comportement d’un gaz parfait.

Le débit massique du Boil-off :

Le débit massique est égal au produit du débit volumique et la masse volumique. Le

volume du Boil-off est obtenu en équivalant gaz naturel aux conditions standards (15°C et 1
bar) ; donc il est nécessaire de calculer le débit volumique aux conditions de pression et
température donné.

PARTIE CALCUL

I)-1- Détermination de la composition du GPL :

4
Le tableau suivant permet de présenter de différent composition du GPL, cas hiver et cas
d’été.

Densité à Masse Cas

Constituant X’Li (%) TVR
15°C moléculaire Été
C2 0.44
C3 34.40
0.5570 8.732 53.29 Juillet 2004
iC4 11.90
nC4 52.48
iC5 0.04
nC5 0.74
Constituant X’Li (%) Densité à TVR Masse Cas
15°C moléculaire Été
C2 0.48
C3 34.37
12.28 0.5569 8.738 53.26 Août 2004
iC4
nC4 52.18
iC5 0.04
nC5 0.66
Constituant X’Li (%) Densité à TVR Masse Cas
15°C moléculaire Eté
C2 0.75
C3 35.02
12.61 0.5559 8.735 53.09 04 Août 2004
iC4
nC4 51.11
iC5 0.02
nC5 0.49
Constituant X’Li (%) Densité à TVR Masse Cas
15°C moléculaire Hiver
C2 0.22
C3 33.00
13.27 0.5579 8.630 53.91 Décembre
iC4 2004
nC4 51.58
iC5 1.46
nC5 0.46
nC5 0.56

5
Constituant X’Li (%) Densité a TVR Masse Cas
15°C moléculaire Hiver
C2 0.39
C3 36.61
12.58 0.5557 8.661 52.99 Janvier 2004
iC4
nC4 49.87
iC5 0.02
nC5 0.56

I-2-a Détermination des conditions optimales de stockage du GPL :

Le GPL, à la sortie de l’réréorefregirant dans le cas d’été à T = 45°C, est stocké dans la
sphère à une température T = 23-25°C avec une pression P = 10.4-10.7 bar.
D’après la construction du sphère de stockage de GPL (spécification de la sphère) on peut
conclue que :
 Pression prévue 12.3 bar (au sommet)
 Température prévue 70°C
 Pression de tarage de la soupape de sécurité : PTS = P prévue + 1.15 = 13.45 bar.

I-2-b Pratiquement, la température ambiante atteint 45°C, on aura que la pression dans la
sphère augmente, le tableautiste qui donne le consigne de vanne de torche ouvrir a P =11.4
bar.
A P =11.4 bar  la température dans la sphère augmente jusqu'à 27°C.
I-2-C Vérification que a P =11.4 bar et T =27°C sont la température d’équilibre des phases
(phase liquide et phases vapeur).
La relation I au dessous permet de vérifier qu’ à T =27°C la températures d’équilibre des
phases à P = 11.4 bar.

n x'

Li  1 ………………. (1)
'
i  1 1  e ( K i  1)

X’Li = La composition du mélange a l’état liquide.

e’ = Le taux d’évaporation.
Ki = Coefficient d’équilibre des phases.

Le Ki sont déterminées d’après le diagramme de Jenny Scheibel = f (p, t).

N.B On dit que T est la température d’équilibre des phases  e’ <= 0.

6
 n x'
Li
Constituant X’Li Ki e’ 
'
i  1 1  e ( K i  1)
C2 0.00750 2.680 0.00750
C3 0.35020 0.875 0.35021
iC4 0.12610 0.385 0.12613
0.00010
nC4 0.51110 0.285 0.51116
iC5 0.00020 0.120 0.00010
nC5 0.00490 0.850 0.00490
Total 1.00000 1.00000

A P =11.4 bar et T = 27°C sont des pressions et températures d’équilibres des phases.

 On fixe une température et une pression maximum.

 P = 11.7 bar et T = 28°C, tel que T+ ΔT (ΔT = 0.5°C) on aura une évaporation du
GPL.
 Pratiquement à P = 11.7 bar la température qui existe dans la sphère est T = 28°C.
 On fixe la température ambiante Tamb = 55°C.

Te = température d’équilibre des phases = 28°C.

Pe = pression d’équilibre des phases = 11.4 bar.

II)- Détermination de la quantité de GPL évaporée (torchée)

7
 Tamb = 55°C

Vanne automatique

Aéroréfrégirant

Te = 28°C
T1 = 40-45°C
T2 = ?
GPL

*- La relation de % Boil-off permet de déterminé la quantité évaporée par jour en fonction de

la quantité de la chaleur totale et le volume qui existe dans la sphère, c'est-à-dire une sphère
pleine qui attend l’expédition.

II-1 -a Calculons la température T2 cotée paroi (métal de la sphère).

Pour déterminer la température T2 coté métal, on appliqué les lois de transfert de
chaleur par convection et par conduction.
** Par conduction : Se fait entre le milieu ambiant et la paroi

dT
q   A …………………………….(2)
CN dx

q
CN = flux de chaleur fait par conduction Kcal / h.

 = conductivité thermique de métal Kcal / h.m.C°.

A = la surface d’échange m2.
dT = gradient de température C°.

A = 4πR2 et dx = dR

dT dR
q  .4R 2  q  4dT
CN CN
dR R2

8
 R T
ex dR 2
q   4  dT
CN
R R2 T
in amb

R
 1  ex  T2
q    4 
T 
CN 
 R R  
 Tamb
in

4  T  T 
q   amb 2
CN R  Rex … (Kcal/h) …………… (3)
in
R .Rex
in

Rin.Rex  : sont les rayons de la sphère intérieur et extérieur respectivement.

** Par conduction : Se fait entre le paroi (le métal) et le fluide.

q  h. A.T …………………………(4)
CV

q
CV = le flux thermique par conduction Kcal/h.
h = coefficient de transfert de chaleur par convection Kcal/h.m2.C°.
ΔT = gradient de température C°.
A = la surface d’échange m2.

A = 4 π R2 et R = D/2

D = diamètre de la sphère m.

2
A = 4πR = 4π (D/2)2  A = πD2………… (5)

1 1 1
  f  e,   
h h 
métal … (Kcal/h.m .°C)…………………. (6)
2
air
D

q
CV 2 e
 
 h. .D 2 T  T ….. (Kcal/h)……………………………….. (7)

On a l’épaisseur de la sphère qui est sont très petit, et la sphère n’est pas équipée d’une
isolation ou à un calorifugeage, c’est à dire la sphère utilise pour stockés le GPL, la sphère un
réservoir place dans un milieu ambiant, alors la température T2 coté paroi ≈ la température
ambiante.

9
Rin = Rex + e  Rin ≈ Rex = D/2 tel que Rex - Rin = e donc

4  T  T  2
 amb 2 D 
3 q   4   T
  T  
CN e  2   amb 2
R2

q  D 2  T  T .......... .......... .......... .(8)

CN  amb 2

* Calculons h :

1 1 1
  f  e,   
h h 
air mét ……. (6)
D

e e e e
moy
f (e,  )  1  2  3 
   

1 1 1
 
Donc (6) sera écrit comme suit h h 
air mét ……………………..(9)
D

hair = coefficient de transfert de chaleur par convection de l’air.

hair = 3BTU/h.ft2.°F = 14.68 Kcal/h.m2.°C.

mét = 260BTU/h.sq.ft.°F/in = 32.22 Kcal/m.h.°C.

1 1 1
   h  1.9 Kcal / h.m2.C
h 14.68 32.22
14.68

T .  T .h
(7) = (8)  T  amb e ……………………. (10)
2 h

55.32.22  28.1.9
T   T  53.5C
2 32.22  1.9 2

II-1-b

10
 Q
 T .24.100
% Boil  off jour  .V
G sph
.H ..................................... (11)

Q = quantité de flux thermique total Kcal/h.

T
 = masse volumique totale du mélange g/Cm3, Kg/m3.
G
∆H = chaleur de vaporisation latente Kcal/Kg.
Vsph = volume de sphère m3.
Pour le volume de la sphère on utilisé à 80% de la volume totale.

1697.4 m3 100%
X 80%  X = Vsph = 1357.92 m3

*Calculons la masse volumique total de mélange a T = 28°C.

15
On a T = 15°C 15 = 0.5559 g/Cm3.
D’après la formule de spécifie gravity on peut calculé 420 . Tel que

15   20  5 et   f   20  ……………………. (12)

15 4 4 

α = 0.0011 donc

 420 = 0.5559 -5.0.0011   420 = 0.5504 g/Cm3.

4t  420    t  20 428  0.5504  0.0011(28  20)

  28  0.5416 g 541.6 Kg 3
4 cm 3 m
*Calculons une chaleur de vaporisation latente ΔH
   53.5   28 …………………………..(13)
2 1

 (T  53.9C  128.28F )  142 BTU LB  78.81 Kcal Kg

2

 (T  28C  82.4F )  129 BTU LB  71.60 Kcal Kg

1

  53.5   28  78.81  71.60  7.215 Kcal Kg

2 1

*Calculons la quantité de chaleur totale

Q  h. A.(T  Te )  h. .D 2 (53.5  28)

T 2
Q  33846 Kcal h
T

11
 33846.24.100
(11)  % Boil  off jour  541.6.1357.92.7.215  15.3 % par jour

80% 1357.92 m3
15.3% X  X = Vi = 259.7022 m3/jour

3 3 Kg
Vi  259.7022 m jour  10.82 m h  5860.62 h

II-1-c Calculons le taux d’évaporation « e’ »

N V (le volume évaporé par jour ) 259.7022

e'  D  i   e'  0.19125
N G (80%)volumede sphére 1357.92
L

II-1-d Calcule de composition de GPL torché :

La formule suivante permet de calculer la concentration de GPL évaporé :

X'
X i'  Li ……………………………..(14)
'
1  e ( K  1)
i
A T =28°C et P = 11.7 bar on a les Ki suivantes :

Ki e’ X'
Li
Composition ’
X Li X i' 
1  e ' ( K  1)
i
C2 0.00750 2.7 0.00566
C3 0.35020 0.9 0.35703
iC4 0.12610 0.39 0.14275
nC4 0.51110 0.29 0.19125 0.59141
iC5 0.00020 0.12 0.00024
nC5 0.00490 0.865 0.00503
Totale 1.00000 1.00005

A P = 11.7 bar et T = 53.5°C

La formule suivante permet de calculer la composition de GPL évaporé

y '  Ki . X i' ………………………………………….(15)

Constituant X’Li Ki e’ Xi’ y '  Ki . X i'

12
 C2 0.00750 4.1 0.00475 0.01948
C3 0.35020 1.52 0.32870 050000
iC4 0.12610 0.74 0.14280 0.10567
0.19125
nC4 0.51110 0.58 0.56568 0.32809
iC5 0.00020 0.285 0.00025 0.00007
nC5 0.00490 0.245 0.00575 0.00141
Totale 1.00000 0.95472

III- Dimensionnement d’un condenseur-refrigerant à « eau »

III- Détermination de la quantité d’eau nécessaire au refroidissement :

G (qV  qL )  G
eau (q40  q25 ) ……………………. (16)
GPL 53.5 28

G
GPL  :- quantité de GPL évaporé…………….Kg/h.

Geau  :- quantité d’eau de refroidissement……….Kg/h.

qV q L  :-L’enthalpie du GPL vapeur et liquide respectivement a l’entrer et la sortie

53.5 et 28
de la condenseur …………Kcal/Kg.

q et q  :- L’enthalpie de l’eau a 40°C et 25°C…….Kcal/Kg.

40 25

L’enthalpie de l’eau et calculé par la méthode suivante :-

q L40 (eau)  C p .t 1.40  q L  40 Kcal Kg

40

13
q L (eau)  C p .t  1.25  q L 25 Kcal Kg
25 25

qV  78.81 Kcal Kg et q L  71.60 Kcal Kg

53.5 28

qL qL
 G  G . 53.5 28
D’après l’équation (16) eau Gpl eau eau
q q
40 25

78.81  71.60
 Geau  5860.62.
40  25

 Geau  2817 Kg h

III-1-b-Calculons de la quantité de chaleur cédée par le GPL :

Qcédé G (qV  q L )  5860.62(78.81  71.60) 

GPL GPL 53.5 28

Qcédé  42255.1 Kcal Kg

GPL

III-1-c Détermination de la quantité d’eau consommée par le condenseur :

Q cédé 42255.1
Geau  GPL   Geau 2817 Kg h
q eau  q eau 40  25
40 25

III-1-d Détermination de la surface d’échange du condenseur :

Q  F .K .TLM …………………………………………………..(17)

Q  :- quantité de chaleur que donne le condenseur chaud………………….Kcal/Kg.

F :- surface d’échange m2.
K :- coefficient global de la transfert de chaleur Kcal/h.m2.°C.
K = 120 Kcal/h.m2.°C.

*-Calculs de la TLM

14
 53.5

28

40

25

Tmax  53.5  40  13.5C

T  28  25  03C
min
T T  T
max  13.5  4.51.8  TLM  max min  13.5  03 
T 03 T 13.5
min 2.3Log max 2.3Log
T 03
min
TLM  07C

Q 42255.1
F  
K .TLM 120.7
D’après l’équation (17)
F 50 m2

15
VI- Dimensionnement d’un ballon de reflux :

G= 5860.62 Kg/h

T = 53.5°C et P = 11.7 bar.

*- La pression au sommet de la colonne et calculer par la méthode suivante

Ps = Pb + ∆P1…………………………….(18)
Ps : la pression de sortie de l’évaporation du GPL.
Pb : la pression de ballon.
∆P1 : la correction de pression.

Pb = Ps - ∆P1 = 11.7 - 0.35

Pb = 11.35 bar

A cette pression on déterminé la température de la colonne, on suppose que une condensation

est totale.

Constituant X’Li Ki Yi' yi'   K Yi ' T

i
C2 0.00750 3.40 0.01948 0.06823
C3 0.35020 1.28 0.50000 0.68400
iC4 0.12610 0.66 0.10567 0.07000
nC4 0.51110 0.46 0.32809 0.16010
iC5 0.00020 0.185 0.00007 0.00001 38°C
nC5 0.00490 0.160 0.00141 0.00025
Totale 1.00000 0.95472 0.98260

T = 38°C et P = 11.35 bar sont la température et la pression de la colonne (ballon de reflux).

- pour le cas d’une condensation totale de tête de colonne :

16
P > 7 bar  3 < L/D < 5

Le volume de distillat (GPL évaporé) sur lequel sera basé le calcul de D et L.

V = 2( V1 + V2 ) tel que
L = la hauteur de la colonne m.
D = le diamètre de la colonne m.
V1 = volume retenu correspondant à cinq (05) minutes de soutirage du distillat.
V2 = Volume retenu correspondant à trois (03) minutes de soutirage de reflux.

 Le volume retenue est utilisé comme reflux  le taux de reflux = 1, c'est-à-dire V1 =

0, donc V = 2V2.
* Calculons la masse volumique de distillat à la température de bulle (Tb =38°C), en utilisant
l’équation suivante :

4t  420    t  20  438  0.5504  0.0011(38  20)

438  0.5306 g 3  530.6 Kg 3
cm m

G
V  GPL  5860.62  3  V  0.56 m3  V  1.12 m3
2 2
 38 530.6 60
4

*L’hypothèse : L =3D.
Donne une première approximation du diamètre :

1 1
D  0.75V 3  0.75(1.12) 3  D  0.78 m

h = 0.2m : hauteur des segments circulaires supérieur et inférieur :

- Calculons le rapport « h/R »

R = D/2 = 0.78/2  R = 0.39 m

h/R = 0.2/0.39  h/R = 0.52.
Le tableau d’évaluation donne pour cette valeur h/R l’aire de segment :

h/R = 0.52  S/R2 = 0.60  S = 0.60 R2 = 0.60 (0.39)2

S = 0.09126

L’aire des deux segments sera

S’ = 2.S  S’ = 2*0.09126 = 0.018252.

F = π.R2 - S’  F = 3.14 (0.39)2 - 0.018252  F = 0.335074.

17
L = V/F = 1.12/.335074  L =3.34 m.

On prend comme longueur de la ballon L = 3m, et de diamètre D = 0.78 m.

Ou L/D = 3/0.78 = 3.

Vers torche

Vanne automatique

h
NH
Te = 28°C
NB
h
GPL
L

18