Kerbouche Ali Kaabouche Boumadiane

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
UNIVERSITE KASDI MERBAH OUARGLA

FACULTE DES SCIENCES APPLIQUEES
DEPARTEMENT DE GENIE DES PROCEDES
Mmoire
Prsent pour lobtention du diplme de
MASTER
Domaine : Science et Technologies
Filire :Gnie des Procds
Spcialit: Ingnierie du Gaz Natural Master Acadmique
Prpar par : KERBOUCHE Ali
KAABOUCHE Boumadiane
Thme
Vrification et optimisation des paramtres de

fonctionnement de la colonne du dbutaniseur T202 a
HASSIRMEL MPP1
Soutenu Publiquement le : 31/05/2016
Devant le jury :
Mr TABCHOUCHE Ahmed Pr Prsident UKM Ouargla
Mr ROUANE Azeddine Mc encadreur UKM Ouargla
Mr KAHOUL Fares Mc examinateur UKM Ouargla
Mr HENNI Abdellah Mc examinateur UKM Ouargla
_
Anne Universitaire : 2015/2016
Avant tout nous remercions "Allah" tout puissant qui nous a donn le courage, la
volont et la force pour accomplir ce modeste travail. Merci de nous avoir clair
le
chemin de la russite.
Nous remercions vivement notre professeur consultant
Mr.RAOUAN Azedin pour son aide prcieuse et ses conseils judicieux
ainsi que tous les enseignants du dpartement de gnie des procdes pour leur
contribution notre formation et leur disponibilit orienter les tudiants.
Nos sincres gratitudes sadressent aussi mes amis et les responsables de
CSTF Hassi Rmel,et particulirement Mrs :
MAHOUR.Dj, ATHAMNIA.F, TITAM .R ,DAIBOUN SAHEL .A ,
BOUNEZRA .R , AIET AHMED .F , EL GOUTNI .F ,KAHOUL . w et
tous les tableautistes et les operateurs sans exception .
Finalement, nous remerciements tous ceux ou celles qui ont contribu de prs
ou de
loin laccomplissement de ce travail.
Ingnierie de Gaze Natural

UNIVERSIT KASDI MERBEH - OUARGLA 2016
Rsum
Suite la chute consquente de la pression du gisement, la composition de la charge tend vers un
changement graduel et significatif avec une diminution de la richesse. L Ce problme a t pris en charge
par la mise en service dune station Boosting, en septembre 2004, pour permettre aux installations de
fonctionner dans de bonnes conditions selon le design.
Nanmoins, ces variations induisent souvent une dstabilisation de la section de fractionnement, ce qui
ncessite une optimisation des paramtres dexploitation et surtout pour le dbutaniseur qui est la partie
la plus sensible du train. La baisse de richesse implique la diminution de la charge dalimentation du
dbutaniseur qui induit un Changement des paramtres de fonctionnement du dbutaniseur.
Mots cls : la pression du gisement, Boosting, le design, une dstabilisation de la section de
fractionnement, le dbutaniseur, train, dalimentation, les paramtres de fonctionnement.
Abstract
Following the substantial fall in reservoir pressure, the composition of the charge tends towards
a gradual and significant changes with a decrease in wealth.
This has been supported by the commissioning of a Boosting station in September 2004, to
allow the equipment to operate in good conditions according to the design.
However, these variations often induce destabilization of the section fractionation, which
requires optimization of operating parameters and especially for debutanizer which is the most
sensitive part of the train. The decline in wealth implies the reduction of the debutanizer
feedstock which induces a changing operating parameters of the debutanizer.
Keywords: reservoir pressure, Boosting, the design, destabilization of the section fractionation,
the debutanizer, the train, feedstock, operating parameters.

4002 .
.

.
( )
.
, , , , , :
,
Nous ddions ce travail :
A nos trs chers parents, sources inpuisables damour et de tendresse, en reconnaissance
des sacrifices consentis avec dvouement pour notre ducation et notre formation .
A nos amis et toute notre famille ainsi qu tous ceux qui nous sont chers veuillez trouver
ici le symbole de lamour et lintime attachement que vous nous portez.
A notre parrain et tous ceux qui taient prs de nous au cours de ce stage, surtout nos
trs chers professeurs qui sont la source de notre inspiration quils trouvent dans ce travail les
sincres tmoignages de notre profonde affectation et de notre haute considration .

Liste de figures et schmas
Liste de figures
Figure I.1 : Production mondiale du GN. 4
Figure I.2 : Consommation mondiale du GN.. 5
Figure II-1: La Carte gomtrique de HASSI RMEL, et les sites des modules. 12
Figure III .1 : Demande nationale en GPL 23
Figure III .2 : Offre prvisionnelle en national de GPL 23
Figure. III .3 : Sphre de stockage 25
Figure V.1 : Simulation du procd de la section. 71
Figure V.2-Profil de temprature.. 72
Figure V.3 : Profil de pression. 72
Figure V.4 : Profil du dbit.. 73
Figure V.5 : Profil de la composition de la charge... 73
Liste de schmas
Schma II .1 : Schma des diffrentes sections du module 1 15
Schma II .2 : schma de Section de stabilisation et fractionnements.. 19
Schma IV.1 : Schma gnral du fonctionnement dune colonne de rectification.. 31
Schma IV.2 : schma de la tte de la colonne. 34
Schma IV.3 : Schma descriptif du fonctionnement de la zone de rectification. 35
Schma IV.4 : schma du fond de la colonne... 35
Schma IV.5 : Schma descriptif du fonctionnement de la zone dpuisement... 36
Liste des tableaux
Liste des tableaux

Tableau I.1 : Composition du gaz naturel de Hassi Rmel 3
Tableau II.1: Capacit de production des diffrents modules... 12
Tableau II.2 : Capacit de stockage des produits finis et rinjection des gaz.. 12
Tableau III.1 : Composition molaire du GPL 21
Tableau V.1 : Composition molaire de la charge dalimentation dbutaniseur 43
Tableau V.2 -Fractions et masse molaire de distillat et de rsidu. 45
Tableau V.3 : Valeurs des Ki du distillat liquide... 46
Tableau V.4 : Valeurs des Ki du distillat vapeur... 47
Tableau V.5 : Valeurs des Ki du rsidu. 47
Tableau V.6 : Valeurs des Ki de la charge 48
Tableau V.7 : Volatilits relatives dans les diffrentes zones de T202. 49
Tableau V.8 : Taux de vaporisation.. 50
Tableau V.9 : Verification du taux de vaporisation.. 51
Tableau V.10 : Calcul de rf. 51
Tableau V.11 : Calcul de rb (taux de rebouillage). 52
Tableau V.12 : Composition du liquide du dernier plateau.. 55
TableauV.13 : Composition massique et molaire du gaz de tte et du rsidu... 55
Tableau V.14 : Calcul de la masse volumique moyenne pour le distillat. 56
Tableau V.15 : Calcul de la masse volumique moyenne pour le rsidu... 56
Tableau V.16 : Calcul de la masse volumique pour le dernier plateau. 57
Tableau V.17 : Calcul de la masse volumique moyenne pour lalimentation... 57
Tableau V.18 : Les masses volumiques 58
Tableau V.19: Les enthalpies 58
Tableau V.20 : Calcul des paramtres critiques 60
Tableau V.21 : Nombre de mole dO2 ainsi que les produits de la combustion(CO2,H2O). 64
Tableau V.22 : Calcul de PCS.. 65
Tableau V.23 : Calcul de la masse volumique de gaz combustible.. 66
Tableau V.24 : Paramtres opratoires du design et calcul de la colonne T202... 69
Tableau V.25 : Composition design de lalimentation du dbutaniseur... 70
Tableau V.26 : la comparaison entre cas design rel et cas design simul... 74
Tableau V .27 : Des diffrentes spcifications des produits. 74
NOMENCLATURE
NOMENCLATURE
GN :gaz naturel
GNL : gaz naturel liqufi
Gpl : gaz de ptrole liqufi
Kwh : kilowatt par hour
Mtep : million de tonnes quivalents ptrole par jour
TVR : Tension de vapeur rad
a : densit de lair Tmoy= 57,61 C
am : viscosit dynamique de lair T= 57 ,61C
a : densit de lair la temprature initial
t : diffrence caractristique de temprature
P1 : Pertes de charges dues aux rsistances des conduites et du condenseur
P2 : Pertes de charges dues aux rsistances des plateaux
Bd : Rsistance thermique de lencrassement externe de la surface des tubes
C : Coefficient qui dpend du type de plateau, de la charge et de la distance entre eux
PC1 ,PC2 : sont les pouvoirs calorifiques moyens de lair dtermins la temprature dentre
et de sortie respectivement, et une pression constante (kcal/ kg).
PCi : Pouvoir calorifique du constituant i
D : Dbit molaire du distillat de la colonne T 202 (dbutaniseur)
Dm : Dbit massique du distillat dfinie
e 0 : Taux de vaporisation molaire de la charge.
e0 : Taux de vaporisation molaire de la charge dalimentation du dbutaniseur.
G : Dbit des vapeurs de GPL entrantes dans larocondenseur E-211.
g 0 : Dbit molaire de la phase liquide de la charge
G 0 : Dbit molaire de la phase vapeur de la charge
Ga : Dbit dair ncessaire pour le refroidissement.
GCOM: Dbit du gaz combustible
GF: Dbit massique du fuel gaz
GR: Dbit du reflux chaud (vapeur provenant du rebouilleur)
gX : Dbit du reflux froid
: Enthalpie de la phase liquide de la charge
: Hauteur du fond suprieur du dbutaniseur
: Enthalpie des vapeurs de GPL 64C.
I
NOMENCLATURE
: Hauteur de la zone de rectification.

: Enthalpie de GPL sortant de E -211 54C.
: Hauteur de la zone dalimentation.
: Hauteur de la zone dpuisement.
: Hauteur du fond infrieur du dbutaniseur [m]
: Hauteur de la jupe.
: Enthalpie du distillat ltat vapeur (vapeur de tte de colonne) [kcal /kg]
: Enthalpie du distillat ltat liquide [kcal /kg]
: Enthalpie du rsidu ltat liquide [kcal /kg]
: Enthalpie de la charge [kcal /kg]
: Enthalpie de la phase vapeur de la charge [kcal /kg]
: Distance entre deux plateaux successives. [m]
: Enthalpie liquide du rsidu. [kcal/kmole]
HR: Enthalpie de leffluent provenant du rebouilleur [kcal /kg]
hR: Enthalpie du rsidu ltat liquide [kcal /kg]
K : Coefficient global de transfert de chaleur [kcal /m2.h.C]
Ki : Constante dquilibre du constituants i une pression et temprature fixe.
N : Nombre de mole de H2O [moles]
P : Pression de tte. [kg/cm2]
PCI : Pouvoir calorifique infrieur [kcal/Nm3]
Pci: Pression critique dun constituant i [kg/cm2]
PCm : Pression critique moyenne [kg/cm2]
PCS : Pouvoir calorifique suprieur [kcal/Nm3]
PCI : Pouvoir calorifique infrieur [kcal/Nm3]
Q : Charge thermique de larocondenseur [kcal / h]
QC: Charge thermique du condenseur [kcal/h]
QR : Charge thermique du rebouilleur H 202 [kcal/h]
QR : Charge thermique du rebouilleur [kcal /h].
R' : Dbit molaire du rsidu de la colonne T 202 (dbutaniseur) [kmole/h]
r 0 : Taux de condensation molaire.
T : Temprature de tte [C]
Tci: Temprature critique dun constituant i [C]
Tcm: Temprature critique moyenne [C]
II
NOMENCLATURE
tmax, tmin : Sont respectivement, la diffrence maximale et minimale des tempratures [C]
V : Reprsente le volume entre le point haut et bas [m3/h]
V1 : Volume correspondant 5 min de soutirage du distillat D. [m3/h]
V2 : Volume correspondant 3min de reflux de tte [m3/h]
4t : Masse volumique de la Charge T= 54C [kg/m3]
Vi : Dbit volumique des vapeurs qui quittent le sommet de la colonne [m3/h]
W : Vitesse admissible des vapeurs [m/s]
Gv : Dbit molaire des vapeurs [kmole/h]
Wa : Vitesse de lair dans la section rtrcie du faisceau tubulaire ailette [m/s]
Wi : Composition molaire massique du fuel gaz.
Xi : Composition molaire des constituants dans le rsidu.
Yi : Composition molaire des constituants dans le distillat.
Z : Facteur de compressibilit
i ,L : Coefficient de volatilit dun constituant i.
L : Masse volumique du GPL liquide T = 54C [kg/m3]
V : Masse volumique des vapeurs [kg/m3]
: Paramtre conventionnel dtermin par approximations successives.
III
Sommaire
Sommaire
Nomenclature... I
Liste des figures et schmas.... IV
Liste des tableaux.... V
Sommaire.VI
INTRODUCTION GENERALE 1
Chapitre I : Gnralit Sur Le Gaz Naturel
I.1-Le gaz naturel dans le monde 2

I.2 - Gaz naturel en Algrie.... 2
I .3 -La production de gaz Naturel.. 3
I.4 -La consommation de gaz naturel...... 4
I.5 - Traitement du gaz naturel 5
I.6 - Caractristiques du gaz naturel 6
I.7 - Types du gaz naturel 8
I.8 - Utilisation du gaz naturel. 8
I.9 - Avantages et inconvnients du gaz naturel 9
Chapitre II : Prsentation de la rgion de Hassi Rmel et MPP1
II.1 -Prsentation de la rgion HASSI RMEL .. 10

II.2 -Situation gographique de HASSI RMEL. 10
II.3 -Les principales tapes du dveloppement de HASSI RMEL. 10
II.4-Organisation du champ de HASSI RMEL.. 11
II.4.1 -Les salles de contrle 13
II.4.2 -Le centre de traitement dhuile (CTH) 13
II.4.3 -Station de rcupration du gaz associ (SRGA).. 13
II.4.4 -CTG-HR- SUD.... 13
II.4.5 -Le stockage et lexpdition des hydrocarbures liquides.. 14
II.4.6 -Centre national de dispatching gaz (CNDG).. 14
II.4.7 -Les Stations Boosting.. 14
II.4.8 -Organisation du laboratoire.... 14
I.5 -Prsentation du module 1 15
II.5.1 -Description du procd de traitement du gaz naturel dans le module1.. 16
II .5.1.1 Zone de rcupration du condenst (Haute pression) 16
II.5. 1.2 Zone de rcupration du GPL et stabilisation du condenst (Basse pression). 17
II.5.2 -Unit de rgnration du glycol 19
VI
Sommaire
II.5.2 .1-Description du Procd 20
Chapitre III : Gnralits sur les produits
III. 1 .GPL... 21
III.1.1 Dfinition du GPL... 21
III 1.2 -Sources du GPL.. 21
III.1 .3 Caractristiques gnrales des GPL.. 22
III .1.4 -GPL en Algrie. 22
III.1.5 La demande nationale de GPL.... 23
III.1.6 -Domaines dutilisation des GPL.... 24
III.1.6.1 -GPL carburant 24
III.1.6.2 -GPL dans la production lectrique. 24
III.1.6.3 -GPL dans les mnages... 24
III.1.6.4 -Le GPL dans la climatisation. 24
III.1.6.5 -GPL dans la ptrochimie 24
III.1.7 -Stockage de GPL 25
III.2 Condenst.... 26
III.2.1 -Caractristiques du condenst 26
III.2.2 -Utilisation du condenst. 27
III.3 -gaz sec (gaz de vente).... 28
III.3.1 -Caractristiques du gaz de vente.. 28
III.3.2 -Utilisation de gaz de vente 28
Chapitre IV : Thorie de la rectification
IV.1-Thorie de la rectification... 30
IV.1.1 Dfinition.. 30
IV.1.2 -Description de la colonne de rectification. 30
IV.2 -Lois fondamentales de distillation..... 32
IV.3 -Calcul dune colonne de rectification 33
IV.3.1 -Bilan matire de la colonne... 33
IV.3.2 -Les paramtres opratoires.... 38
IV.3.3 -Taux de vaporisation et composition des phases liquide et vapeur de la charge... 39
IV.3.4 -Fonctionnement reflux 39
V.3.5 -Taux de reflux optimal... 39
VII
Sommaire
V.3.6 -Bilan thermique de la colonne.... 40

V.3.7 -Calcul du nombre de plateaux par la mthode de GILLILAND.... 41
V.3.8 Calcul des dimensions de la colonne. 41
V.3.9 -Hauteur de la colonne... 42
Chapitre V : Parties de Calculs
V.1 -Bilan de matire de la colonne T202 43

V.2 -Dbit et composition du distillat et du rsidu.. 44
V.3 -Pressions dans les diffrentes zones de la colonne... 46
V.3.1 -Calcul de la pression et les tempratures aux diffrentes zones... 46
V.4 -Volatilits relatives... 46
V.5 -Taux de vaporisation et composition des phases vapeur et liquide dalimentation. 49
V.6 -Calcul du rgime minimal 51
V.6.1 -calcul de taux de reflux et taux de rebouillage. 51
V.6.2 -Calcul de taux de reflux minimal. 51
V.6.3 -Calcul de taux de reflux optimal.. 52
V.6.4 -Calcul de taux de rebouillage minimal. 52
V.6.5 -Calcul de taux de rebouillage optimal. 52
V.6.6 -Calcul du nombre de plateaux thoriques dans la colonne.. 52
V.7 -Bilan thermique de la colonne.. 53
V.7.1 -Calcul des diffrents dbits de reflux de la colonne. 54
V.7.2-Les enthalpies des diffrentes phases liquides. 56
V.7.3 -Calcul de lenthalpie des vapeurs quittant le sommet de la colonne. 58
V.7.4 -Calcul de la charge thermique du condenseur 58
V.7.5 -Calcul de lenthalpie de la charge.. 58
V.8 -Dimensionnement de la colonne 59
V.8.1 -Calcul du diamtre de la colonne. 59
V.8.2 -Dtermination du facteur de compressibilit. 60
V.8.3 -Calcul de la vitesse admissible des vapeurs. 60
V.8.4 -Calcul de la hauteur de la colonne. 62
V.9 -Calcul Du Ballon De Reflux V209 62
V.9.1 -Dbit massique du distillat. 63
V.10 -Calcul du vrification du four 64
V.10.1. -Calcul de la quantit de chaleur apporte par le gaz combustible 64
VIII
Sommaire
V.10.1.1 Dbit du gaz combustible. 64

V.10.2 -Calcul du pouvoir calorifique suprieur PCS. 65
V.11 Partie simulation... 68
V.11. 1 -La simulation dun procd... 68
V.11 .2 Choix du modle thermodynamique... 68
V.11.3. Vrification des paramtres design du dbutaniseur... 68
V.11.3.1. Composition de charge du dbutaniseur. 69
V.11.3.2. Procdure de la simulation.. 69
V.11.4 -Profil de temprature... 72
V.11.5 Profil de pression.. 72
V.11.6-Profil du dbit.. 73
V.11.6.1-Profil de la composition de la charge....... 73
Conclusion Gnrale.................................................................................. 74
Annexe..
Bbliographie
IX
Introduction gnrale
INTRODUCTION GENERALE :
Aujourdhui, lindustrie gazire aborde le nouveau millnaire dans un contexte de

restructuration de la scne nergtique dans son ensemble. Sous limpulsion de nouvelles
opportunits technologiques et commerciales, lindustrie se repositionne et entame une
nouvelle phase de croissance et a fait merger une nouvelle conomie du gaz.
Les stratgies des intervenants, la globalisation, les fusions, les alliances, se sont
combines pour donner au gaz le dynamisme ncessaire la matrialisation des espoirs quil
nourrit pour obtenir. Lnergie du 21eme sicle.
Au cours des dix annes coules, des hausses majeures travers le monde ont
fortement influence la croissance de la demande du gaz naturel.
Cette activit du gaz a enregistre un dveloppement modeste pendant la dcennie 90,
en laugmentant un taux moyen annuel de 1.7% entre 1990 et 2003. Parmi les nombreux
facteurs qui convergent a lexpansion future du gaz ; on cite de nouvelles innovations
technologiques qui permettent rduire les couts et valorises les nouveaux dbouchs pour
le gaz.
La technologie ayant servie convertir le gaz naturel en produits ptroliers nest pas
nouveau. De nos jours, dans un contexte de hausse du prix du ptrole brut et de la baisse des
investissements, cette technologie offre des opportunits prometteuses pour la valorisation, des
rserves de gaz dlaisses.
En fait, la conversion chimique du gaz en produits de qualit suprieure avec son
potentiel de marche des produits ptroliers devient une alternance commerciale que plusieurs
compagnies sont en train dtudier au dbut de sicle.
Les raisons qui nous ont pousses entamer ce travail ont pour objectif de la ralisation
dune unit de stabilisation du condensat aprs sa sortie du procd de sparation.
Cette stabilisation par distillation permet datteindre les spcifications des produits
stabiliss avec une meilleure rcupration en hydrocarbures intressants (C5+) et une
meilleure limination des composants lgers haute tension de vapeur (C2, C3 et C4).
1
Chapitre I Gnralit Sur Le Gaz
Naturel

Chapitre I Gnralit Sur Le Gaz Naturel
I.1-Le gaz naturel dans le monde

Le gaz naturel est une nergie fossile comme la houille, le charbon ou le lignite.
Cest un mlange dont le constituant principal, de 75 % 95 %, est le mthane (CH4). Il
est donc compos majoritairement dhydrogne et de carbone, do son nom
dhydrocarbure.
Pendant longtemps le gaz naturel a t considr comme un parent pauvre
de lexploitation du ptrole, et lon voyait sur les champs ptroliers mondiaux, des
centaines de torchres brlantes lair libre la plus grande partie du gaz naturel
produite en mme temps que le ptrole .
La plus part des pays ont rapidement ragi contre un tel gaspillage
dnergie en interdisant de brler ainsi le gaz lair libre et en imposant soit
une rinjection de ce gaz dans ces formations ptrolires souterraines, soit lutiliser
directement comme source dnergie.
Au cours de ces deux dernires annes la consommation du gaz naturel a
connu un accroissement particulirement rapide ; les besoins mondiales annuels en
nergie devraient passer de 180 millions de tonnes quivalents ptrole par jour
(mtep/j) en 1995 317 mtep/j en 2020 ; ce qui reprsente une augmentation de 75
% avec 3.3 % de croissance annuelle; le gaz naturel aura la plus grande part et
pourrait tre la premire source dnergie en 2050. [1]
I.2 - Gaz naturel en Algrie
Les rserves algriennes prouves en gaz sont de lordre de 4,550

milliards de m 3, lAlgrie est lun des pays pilote exportateur de gaz naturel et de
GNL dans le monde Le domaine minier algrien est particulirement riche en gaz
naturel. Les rserves gazires algriennes reprsentent prs de 3% du total des
rserves mondiales.
Prs de 66% du gaz actuellement commercialis par SONATRACH
proviennent du gisement gazier gant de Hassi Rmel qui possde un gaz
condenst.
Champs exploits par SONATRACH
Les champs gaziers Hassi Rmel, Rhourde Nouss, Gassi Touil, Alghar et
TFT exploits par SONATRACH .
2
Sachant que la composition du gaz naturel ainsi que sa densit diffre dun
gisement un autre, Hassi Rmel les puits exploits par lunit de traitement
module1 ont une densit de 0,6 et une composition illustre dans le tableau ci-
dessous :
Tableau I.1 Composition du gaz naturel de Hassi Rmel. [2]
Constituents Percentage Molaire %
N2 5,56
CO2 0,2
CH4 78,36
C2H6 7,42
C3H8 2,88
i-C4H10 0,62
n-C4 H10 1,1
i-C5H12 0,36
n-C5H12 0,48
C6H14 0,59
C7H16 0,56
C8H18 0,45
C9H20 0,37
C10H22 0,27
C11H24 0,21
C12H26 0,57
Total 100
I .3 -La production de gaz Naturel

La production mondiale totale en 2006 tait de 2865 milliards de mtres
cubes en croissance de 3% par rapport l'anne prcdente. Une croissance
suprieure la moyenne annuelle sur la priode 1990-2006 (2.3%).
La production gazire mondiale enregistre une bonne performance en
2012, qui s'explique par l'Amrique du Nord (33 % de l'accroissement total),
suivie du Moyen- Orient (23 %), l'Asie Ocanie (17 %), et l'Afrique (13 %). Les
volumes extraits en 2012 par les pays producteurs du Moyen-Orient et d'Asie ont
continu d'afficher des taux de progression significatifs, mais ceux-ci sont en
retrait par rapport aux annes prcdentes.
3
Une hausse de la production de gaz naturel dans le monde est attendue en
raison des projets d'exploration et d'expansion planifis en anticipation d'une
demande future haussire.
Figure I.1 :Production mondiale du GN [3]
I.4 -La consommation de gaz naturel

Le gaz naturel reprsente prs du quart de la consommation nergtique
mondiale. Comme le montre clairement le graphique ci-aprs, cette consommation
a considrablement augmente au cours des 40 dernires annes.
Les principaux pays consommateurs de gaz naturel en 2006 taient les
Etats-Unis avec 22% de la consommation totale et la Fdration de Russie, avec
15.1%. L'Amrique du Nord et l'ex Union sovitique ont consomm ensemble
environ 49% de la production de gaz naturel. La part de l'Europe (UE 25) dans
la consommation totale de gaz naturel tait de 16.3%. Ces trois zones ont
compt elles seules pour prs des deux tiers de la consommation globale en
2006.
La consommation gazire aux tats-Unis pourrait grimper de
plus de 4 % en 2012, en raison d'une hausse considrable de l'utilisation
du gaz dans le secteur lectrique, consquence des prix du gaz plus
comptitifs par rapport au charbon. Selon l'agence amricaine EIA, cette
situation pourrait voluer court terme, en raison d'une hausse probable
des prix du gaz.
La Chine devrait voir sa consommation augmenter de plus de 13 % en 2012
et ses importations crotre de prs de 40 % . La Chine est devenue importatrice en
4
2006. Depuis, sa dpendance aux importations a bondi de 2 % presque 30 %
aujourd'hui. En aot 2012, la Chine a commenc importer du gaz en provenance
d'Ouzbkistan.
Les quantits livres cette anne nont pas t spcifies, mais un contrat
sign entre les deux pays en juin de cette anne portait sur la fourniture de 10
Gm3/an partir de 2013.
Faisant exception sur le march mondial, la consommation europenne
devrait, pour la deuxime anne conscutive, enregistrer une baisse, toutefois
moins prononce qu'en 2011. En 2011, la consommation gazire de l'UE-27 avait
chut de 10 % selon CEDIGAZ. En 2012, CEDIGAZ anticipe une baisse de la
consommation gazire europenne (Turquie incluse) denviron 3% par rapport
l'anne prcdente.
Figure I.2 : Consommation mondiale du GN.
I.5 - Traitement du gaz naturel

Il est ncessaire de faire subir au gaz naturel les traitements suivants
afin de lui donner les spcifications requises pour son utilisation :
1. Dtente la pression de traitement et sparation de leau.
2. Dsulfuration pour liminer les gaz acides et les mercaptans.
3. Elimination des fractions lourdes ou dgazolinage.
Ltape de traitement qui suit dpend du mode de transport adopt. Le gaz naturel
ainsi que ses diffrentes fractions peuvent tre transports sous diverses formes :Gaz naturel
comprim (transport par gazoduc), gaz naturel liqufi (GNL), gaz de ptrole liqufi (G.P.L)
et produits chimiques drivs (mthanol, ammoniac).
5
Les procds de traitement de gaz naturel sont caractriss par des moyens de
rfrigration. A Hassi Rmel, deux procds sont utiliss pour le traitement de gaz. Le
premier procd est appel PRITCHARD quant au deuxime, c'est le procd HUDSON. Ce
dernier est le plus utilis dans les modules de la rgion.
Procd PRITCHARD : il consiste refroidir le gaz par change thermique et une

dtente Isenthalpique adiabatique (Polytropique), en utilisant une vanne appel JOULE-
THOMSON pour atteindre des tempratures infrieures ou gale -18C.
Procd HUDSON : il est bas sur le fractionnement par diffrence de temprature
et une dtente en utilisant le TURBO-EXPENDER qui peut atteindre une temprature
infrieure ou gale -40C.
Le procd HUDSON est plus performant que le procd PRITCHARD, il permet
une meilleure rcupration des hydrocarbures liquides (condenst, GPL). Il est utilis dans les
modules 2 ,3et 4 et le procd de PRITCHARD est utilis dans les modules 0 et 1.
Le gaz de Hassi RMel est un gaz condenst, par consquent les produits qui
dcoulent de son traitement sont : le gaz propane liqufi (GPL) et le condenst stabilis .
I.6 - Caractristiques du gaz naturel
Au stade final de son exploitation le gaz naturel rpond aux caractristiques
suivantes :
a) Masse volumique
Cest la masse de lunit de volume du gaz exprim en Kg / m 3, elle est fonction de
la temprature et de la pression, en dehors des conditions spcifies on se rfre des
conditions dites normales ou standards :
Conditions normales : T = 0C, P = 1atm.
Conditions standards : T = 15C, P = 1atm.
1m3(n) = 1.057 m3(st)
b) Volume massique
Reprsente le volume occup par unit de masse de gaz, il est donc linverse de la
masse volumique et sexprime en m3 / Kg.
c) La densit
Elle est dfinie pour un gaz comme tant le rapport de sa masse volumique celle de
lair dans des conditions bien dtermines de temprature et de pression, comme elle peut tre
obtenue partir de sa masse molculaire que lon peut dfinir partir de sa composition
chimique en utilisant la relation :Densit du gaz = masse molculaire / 28.966
6
d) Pouvoir calorifique
Cest la quantit de chaleur dgage par la combustion dune unit de volume du gaz,
mesure dans les conditions de rfrence, le pouvoir calorifique pour le gaz naturel sexprime
en Joules/m3. [4]
On distingue deux pouvoirs calorifiques :
1. Pouvoir calorifique suprieur PCS
Il correspond la chaleur dgag lorsque tous les produits de combustion sont

ramens la temprature ambiante, leau forme tant ltat liquide.
2. Pouvoir calorifique infrieur PCI
Il correspond la combustion dans laquelle leau reste ltat vapeur, le PCI diffre
du PCS dune quantit de chaleur latente de vaporisation de leau.
e) Composition chimique
Indique la nature des composs dhydrocarbures et dautres constituants du gaz, et

leurs importances relatives dans le mlange par lintermdiaire de leurs fractions volumiques
ou molculaires.
La composition chimique dun gaz est utilise pour ltude de vaporisation, et elle
sert aussi calculer certaines de ses proprits en fonction de la pression et de la temprature
(compressibilit, densit) et dfinir les conditions de son traitement lors de son exploitation.
Elle est dtermine par mthode chromatographique en phase gazeuse.
f) Composition du gaz naturel
Les diffrentes compositions du gaz naturel sont dues la diversit de ses origines, qui
peuvent tre :
U n gaz bactrien
Le mode bactrien est d laction de bactries sur les dbris organiques qui
saccumulent dans les sdiments, le gaz ainsi form est appel gaz bactrien ou biochimique.
7
U n gaz thermique
Au cours de lvolution des bassins sdimentaires, les sdiments sont ports des
tempratures et pressions croissantes, subissant ainsi une dgradation thermique qui va donner
ct des hydrocarbures, une large gamme de composs non hydrocarbons.
U n gaz inorganique
Le mode inorganique de formations des gaz dhydrocarbures reste trs secondaire.

Les gaz volcaniques ou les sources hydrothermales contiennent parfois du mthane et les
inclusions fluides des minraux des roches mtamorphiques ou magmatiques renferment
souvent des hydrocarbures lgers, principalement du mthane. [4]
I.7 - Types du gaz naturel
Lapparition dune phase liquide dpend des conditions de temprature et de pression

dans le rservoir et en surface, ce qui conduit distinguer les types suivants
Gaz sec: Ne formant pas de phase liquide dans les conditions de production, le gaz
est concentr en mthane et contient trs peu dhydrocarbures lourds.
Gaz humide: Formant une phase liquide en cours de production dans les
conditions de surface, il est moins concentr en mthane.
Gaz a condenst: Formant une phase liquide dans le rservoir en cours de
production, la phase condense est riche en constituants lourds.
Gaz associ: Coexistant dans le rservoir avec une phase huile (gisement de
ptrole). Le gaz associ comprend le gaz de couverture (phase gazeuse prsente dans le
rservoir) et les gaz dissous.
I.8 - Utilisation du gaz naturel
Le gaz naturel est une source d'nergie polyvalente qui peut tre employe dans des
domaines trs varis, tel que :
Le secteur industriel
A la fois comme combustible et comme matire premire, le gaz naturel est utilis
dans lindustrie chimique, notamment pour la ptrochimie et le raffinage, son utilisation dans
les pays industrialiss est de 25%, titre dexemple : la synthse dammoniac partir du gaz
naturel permet de fabriquer des engrais pour lagriculture. la synthse du mthanol partir du
gaz naturel est utilise en chimie de spcialits et comme base dadditif des essences.
8
La production dlectricit
Depuis une dizaine dannes, le secteur lectrique est devenu le moteur principal de
laugmentation de lutilisation du gaz naturel dans le monde ainsi le gaz naturel offre un bilan
environnemental trs favorable dans la production dlectricit.
Le secteur domestique
Le gaz naturel est traditionnellement utilis comme combustible dans la roduction de
chaleur pour la cuisson ou le chauffage. 30 % au niveau mondial est destine au secteur
rsidentiel, en particulier la production deau chaude et au chauffage des particuliers.
Le secteur dautomobile
Lutilisation du gaz naturel est particulirement intressante car elle ne ncessite pas de
transformation majeure du moteur. Les moteurs au gaz naturel offrent par ailleurs un bon
rendement nergtique combin un potentiel important de rduction dmissions de CO2.
I.9 - Avantages et inconvnients du gaz naturel
Le gaz naturel comme toute source dnergie possde des avantages et des
inconvnients, parmi ces avantages le gaz naturel est une source dnergie relativement
propre et bon march .
Le gaz naturel possde un pouvoir calorifique raisonnable (PCS=10 kwh/Nm3) et
permet dengendrer de nombreux produits de base pour l'industrie chimique et ptrochimique,
comme l'hydrogne, le mthanol et l'ammoniac.
Quand aux inconvnients du gaz naturel hormis le fait que ses rserves sont limites
et son stockage est trs dlicat, le prix du gaz naturel est index sur celui du ptrole.
9
Chapitre II Prsentation de la rgion de
Hassi Rmel et MPP1

Chapitre II Prsentation de la rgion de Hassi Rmel et MPP1
II.1 -Prsentation de la rgion Hassi Rmel

Le champ gazier de Hassi RMel est le ple conomique le plus important
au niveau national et international. Il est lun des plus grands gisements de gaz
lchelle mondiale, on le considr comme le poumon de lAlgrie, il est rpartie
en trois secteurs : Nord, Centre et Sud. [5]
II.2 -Situation gographique de Hassi Rmel
Hassi RMel, port du dsert, se trouve 550 km au sud dAlger, entre les wilayas de
Ghardaa et Laghouat. Dans cette rgion relativement plate du Sahara, laltitude
moyenne est denviron 750m et le climat est caractris par une pluviomtrie
faible (140 mm par ans) et une humidit moyenne de 19% en t et 34% hiver. Les
amplitudes thermiques sont importantes et les tempratures varient entre 5C en
hiver et +45C en t .
II.3 -Les principales tapes du dveloppement de Hassi Rmel
Le dveloppement du gisement de HASSI RMEL a t ralis en plusieurs tapes, rpondant

lvolution conomique du pays et au dveloppement technologique du march du gaz
naturel.
1961 : Ralisation de 2 units de traitement du gaz avec une capacit de 1,3 milliards
m3 / an.
1969 : Ralisation de 4 units supplmentaires pour augmenter la capacit 4
milliards m3/ an.
1972- 1974 : Ralisation de 6 units supplmentaires pour arriver une capacit
de 14 milliards m3/ an.
1975- 1980 : Russir et rendre effectif le plan de dveloppement qui a vis les
objectifs suivants :
Augmenter la capacit de traitement du gaz de 14 94 milliards m3/ an.
Maximiser le GPL et rcupration du condenst.
1985 : Lunit de la rcupration de gaz dvasement et production du GPL
des modules 0 et 1.
1981- 1993 : Ralisation du centre de traitement dhuile.
1987- 2000 : Le dbut en haut des centres de traitement dhuile.
1999 : Unit de rcupration du gaz de lassocies.
10
2001 : Dmarrage du projet Boosting.
2005 : Mise en service du Boosting.
2013 :Dmarrage du projet unit de glycol
2014 : Mise en service dunit de glycol
II.4-Organisation du champ de Hassi Rmel
Le champ de Hassi Rmel comprend trois secteurs :

Secteur Nord :
Le module 3.
L'unit de dpropaniseur.
Une station de compression nord.
CTH (Centre de Traitement dHuile).
Secteur centre :
3 modules 0, 1 et 4.
CSTF (Centre de stockage et transfert facilites).
CNDG (Centre de notionnel dispatching de gaz )
Lunit de la phase B .
Secteur sud
Le module 2.
Une station de compression sud.
Djebel Bissa.
HR sud.
11
Figure II-1: La Carte gomtrique de Hassi Rmel, et les sites des modules. .[6]
Tableau II.1: Capacit de production des diffrents modules.
Produits Modules Modules Modules Module Djebel Bissa

I II, III IV 0
Gaz sec (106 Nm /J) 60 37 60 20 6
Condenst (T/ J) 3200 1500 3800 1200 360
GPL (T / J) 1200 720 1500 500 -
Tableau II.2 : Capacit de stockage des produits finis et rinjection des gaz.
Condenst (m) 224000
GPL (m) 68000

Capacit de rinjection par jour (Nm) 180
12
II.4.1 -Les salles de contrle
Cest partir des salles de contrle que sont donnes toutes les instructions aux oprateurs
sur site en se basant sur les indications des instruments et des diagrammes denregistrement.
Les salles de contrles gres par un systme numrique DCS pour les modules I, II, III et
IV.
En phase de dmarrage, les vannes stratgiques du processe et celles des machines tournantes
sensibles sont actionnes normalement partir des rgulateurs en salle de contrle pour assurer
une volution graduelle des diffrents paramtres (pression, temprature, dbit, vitesse de
rotationetc.).
Cette opration en semi-automatique sera effectue jusqu' la mise en service de tous les
quipements la stabilisation de tous les paramtres du process leur valeur dexploitation
normale. Cest ce moment seulement que tous les rgulateurs de commodes sont disposs
en contrle automatique, en prenant soit dajuster minutieusement les points de consigne, et en
sollicitant le laboratoire pour les diffrentes analyses de contrle de conformit des points avant de
procder leur expdition.
II.4.2 -Le centre de traitement dhuile (CTH)
Cinq units sont implantes HASSI RMEL (CTH1, CTH2, CTH3, CTH4 et CTH SUD elles
sont dotes de manifold dont des sparateurs et de bacs de stockage qui permettent de produire le
ptrole brut et le transfrer au module 0 pour stockage intermdiaire avant quil soit expdi
vers la station de pompage et de transport SPT .
II.4.3 -Station de rcupration du gaz associ (SRGA)
Cette unit a dmarr le 18 avril 1999 avec une capacit de 4000 m par jour, elle comporte
quatre turbocompresseurs.
II.4.4 -CTG-HR- SUD

Ce projet comporte une unit de traitement de gaz naturel destin traiter le gaz brut pour
obtenir du gaz sec et la rcupration des hydrocarbures liquides estims 2830 millions de tonnes
par mois.
II.4.5 -Le stockage et lexpdition des hydrocarbures liquides
Le condenst et le GPL, produits par tous les modules sont achemins vers le centre de
stockage et de transfert CSTF qui se trouve dans la zone centrale de HASSI RMEL, ou ils
13
sont dbarrasss des ventuelles quantits deau rsiduelles et comptabilises avant quils
soient expdis vers SPT , puis vers Arzew.
II.4.6 -Centre national de dispatching gaz (CNDG)
Son rle est la collecte de toutes les quantits du gaz produites au niveau de HASSI RMEL
provenant de lextrme sud via le gazoduc GR1 . Ces quantits sont ensuite distribues
vers les centres de consommation, comme les centres GNL , SONALGAZ
et ltranger ; lItalie via la Tunisie GR1 , GR2 et lEspagne via le Maroc GME .
II.4.7 -Les Stations Boosting
La pression dentre du gaz brut aux modules dcrot avec le temps, ce qui influe sur la
quantit et la qualit des produits de
chaque catgorie, et sur lusine car il est capable de fonctionner a une pression dentre
comprise entre 140 100 bars. Les stations Boosting ont pour but la compression de ces gaz
bruts issus des puits producteurs de 100 bars 120 bars afin davoir une dtente importante, dont
la rcupration optimale des liquides. Les modules 2 et 3 ont leurs propres stations Boosting,
les autres modules 0,1 et 4 ont un Boosting commun.
II.4.8 -Organisation du laboratoire
Les travailleurs de laboratoire sont des techniciens assurant le systme de travail par quipe ; Et
des ingnieurs responsables sur la bonne marche du laboratoire, Ils effectuent des analyses des
produits finis et des autres produits utiliss lors du traitement du gaz naturel.
Ces analyses sont faites pour vrifier la conformit des produits aux spcifications donnes et
dtecter les anomalies de fonctionnement du processe pour optimiser les paramtres.
Pour cela, les laborantins suivent ces tapes :
Echantillonnage : Dans chaque module existent des points d'chantillonnage o les
laborantins ramnent leurs chantillons dans les capacits appropries a chaque produit.
Analyse : Les laborantins suivent les modes opratoires dans les analyses des produits
avec une grande prcaution car ils manipulent des produits inflammables et dangereux capable de
provoquer facilement des incendies.
Enregistrement des rsultats : les rsultats obtenus sont compars aux spcifications et
enregistrs sur un registre et saisies sur un micro-ordinateur en prcisant l'heure d'analyse, les
rsultats sont transmis en salle de contrle.
14
I.5 -Prsentation du module 1
Notre projet de fin dinduction a t accompli au sein du module 1. Ce dernier a t ralis par
la compagnie amricaine Stone and Webster corporation en 1976 et sa capacit de production
totale par jour avoisine 60 millions de standard m 3 de gaz de vente,1700 tonnes de GPL et 4889
tonnes de condenst .[6]
Le module1 est constitue principalement de trois trains identiques : I-II-III chacun deux assure
la production du gaz sec, GPL et condenst.
Le schma ci-dessous reprsente les diffrentes sections du module I
Schma II .1 : Schma des diffrentes sections du module 1.
Boosting : Cest une unit installe rcemment afin daugmenter la pression dentre
du gaz brut allant de 100 105 kg/cm2 pour pouvoir alimenter les modules.
Au niveau du Boosting, le gaz est dbarrass dune grande partie deau libre et le condenst dans
des ballons de sparation, ainsi aprs compression il est envoy vers le diffuseur V 201 (entre du
module 1).
Train: Section du traitement du gaz naturel comporte 3 trains identiques (I.II.III).
Salle de contrle : Elle est munie dun systme numrique DCS (Digital Contrle Systme)
et cest partir de la salle de contrle que sont donnes toutes les instructions aux oprateurs sur
site, sur la base des indications des instruments et des diagrammes denregistrement.
Les utilits : L'exploitation du module '1', comme les autres modules exige la mise en service
des utilits, elle fournit lair instrument, l'air service, l'eau de refroidissement, l'eau potable, l'eau
de service et l'eau incendie ainsi que le gaz inerte .
Bacs de stockage : Trois rservoirs de stockage ( toit fixe) existes (S-204 A/B/C) dans le
module 1 et ils sont destins au stockage intermdiaire du condenst.
Unit de rgnration du Glycol : Le Di-thylne Glycol joue un rle trs important
dans la marche du procd de PRITCHARD, son cot est trs important, donc cest pour cela
15
quil a t mis au point une unit de rgnration du DEG afin den rcuprer le maximum. La
rgnration est base sur la dshydratation par chauffage pression normale.
Phase B : Cette unit a pour but de rcuprer les gaz moyenne pression de 20 et 21kg/cm
du module 1 et de fournir les frigories ncessaires pour le refroidissement du gaz permettant ainsi
une meilleure rcupration du GPL au niveau du module.
Laboratoire danalyse : Le rle dun laboratoire consiste effectuer des analyses des
produits finis (Gaz, GPL et Condenst) et sassurer quils rpondent aux normes et
spcifications demandes.
II.5.1 -Description du procd de traitement du gaz naturel dans le module1
Les installations du MPPI sont composes de deux sections

a) Section de sparation du condenst (Haute pression).
b) Section de rcupration du GPL et stabilisation du condenst (Basse pression), cest la
partie qui nous intresse dans notre projet.
II .5.1.1 Zone de rcupration du condenst (Haute pression)
Le gaz dalimentation du module 1 provenant du Boosting arrive une pression de 102

kg/cm2 et une temprature de 62C sous forme dun mlange monophasique, passe par le
diffuseur darrive V-201 o il est divis dune manire gale et envoy vers les trois trains
parallles. Le gaz pntre dans le sparateur dadmission V-202, sous le contrle de la vanne de
pression PV-7 une pression de 100 bars et une temprature dentre 59C ou il subit une
sparation flash, pour sparer le liquide du mlange dalimentation et limin leau libre qui est
vacue vers le puisard S-210 sous le contrle de la vanne LV-1.
Les hydrocarbures liquides provenant du V-202 ainsi que les liquides envoys de la station
Boosting et qui constituent la plus grande quantit subissent une dtente isotherme jusqu
une pression de 27 kg/cm2 avec une temprature de 54C dans le sparateur de condenst
V-203 et les vapeurs issues de ce dernier sont envoyes vers la colonne T-201 sous le contrle
de la vanne de pression PV-9.
Le liquide provenant du V-203 contrl par la vanne de dbit FV-12, est envoy au 18me plateau
du dthaniseur T-201 aprs avoir t chauffe dans lchangeur E-216 ct calandre
Le gaz brut provenant du sparateur V-202 est dshydrat et refroidi simultanment
une temprature de -1C dans trois batteries dchangeurs (gaz/gaz) E-201, E-202(A/B/C/D) et
(gaz /liquide) E-203(A/B) Le gaz sortant des changeurs passe dans deux Chillers (E-204 et E-205)
avec du propane comme fluide frigorigne pour atteindre une temprature de -9C. La
dshydratation du gaz dalimentation est effectue par linjection directe et de manire uniforme
16
dune solution de di-thylne glycol (DEG) 86% dans la plaque tubulaire de chaque
changeur afin dviter la formation des hydrates.
A la sortie des chillers, le gaz passe par la vanne Joule Thomson PV12 o il subit une dtente
isentropique la pression de 78 bars et une temprature de -18 C, puis entre dans le sparateur
froid V-204 o se produit une sparation entre les deux phases gaz/liquide.
Le gaz rcupr du V-204 qui a une teneur en eau <50ppm, refroidit contre- courant le gaz
brut provenant du V-202 dans les changeurs E-201 et E-202 pour rcuprer les frigories quil
contient, une fois rchauff, il pntre dans le pipeline du gaz de vente avec une pression de 72
kg/cm2 et une temprature de 60C
Le liquide (mlange de condenst et de glycol) issu du V-204 est prchauff dans lchangeur
E-203 (gaz/liquide) ct calandre jusqu une temprature de 30 C et une pression de 25
kg/cm2.
Le liquide est ensuite admis dans le ballon sparateur V-205, ou il va avoir une sparation tri-
phasique (gaz- condenst glycol hydrat). Le condenst et le gaz sont envoys au 18me plateau
du dthaniseur T-201 sous le contrle de la vanne de dbit FV-32 et la vanne de pression PV-
13B respectivement, alors que le glycol hydrat est envoy vers lunit de rgnration du glycol
travers la vanne de niveau LV-8 [6]
II.5. 1.2 Zone de rcupration du GPL et stabilisation du condenst (Basse pression)
Cette zone a pour but de stabiliser le condenst et de produire du GPL partir du liquide
rcupr du gaz brut dans la zone de sparation du condenst, suivant les spcifications et les
normes tablies pour chaque produit.
Le condenst stabilis provenant du dbutaniseur est refroidi dans les Chillers E-208 et E-
209 par du propane et admit ensuite au sommet du dthaniseur T-201 pour une meilleure
extraction des C3+. Les vapeurs de tte du dthaniseur (gaz) sont envoyes vers la PHASE
B pour tre comprims et dshydrat puis rinjects dans le gaz de vente.
Le reflux latral du dthaniseur qui est assur par la pompe P-202A/B, est refroidi dans le
Chiller E-207.
Le rebouillage latral du T-201 est assur par lchangeur E-206 contre courant (ct calandre)
avec les produits de fond du dbutaniseur T-202.
Le rebouillage principal est fourni par le rebouilleur du dthaniseur H-201 laide de la
pompe P-201 A/B/C.
17
produits de fond du dthaniseur T-201 sont prchauffs dans lchangeur E-217 (ct
calandre) puis envoys autant que charge au 27me plateau du dbutaniseur T-202.
Les produits de tte provenant du T-202 sont condenss dans larocondenseur E-211 une
temprature de 51C. Le GPL liquide est ensuite collect dans un ballon de reflux du dbutaniseur
V-209, puis pomp par la pompe de reflux P-205A/B/C.
Une grande partie de ce liquide (GPL) est utilise comme reflux pour le dbutaniseur et le reste
comme production vers CSTF si le produit est On Spec , ou vers les bacs de stockage du
module 4si le produit est Off Spec .
Le rebouillage du T-202 est fourni par le rebouilleur du dbutaniseur H-202 travers la
pompe P-204 A/B.
Les produits de fond condenst stabilis du T-202 sont refroidis partiellement dans le E-
217 pour prchauffer la charge du dbutaniseur(ct calandre) puis dans le E-206 pour le
rebouillage latral du dthaniseur enfin dans larorfrigrant E-210 ou ils sont refroidis
environ 45C.
Le condenst stabilis transite par la V-208 ou une partie est pompe par les pompes de reflux P-
203A/B vers la tte du dthaniseur comme huile dabsorption et le reste est envoy au stockage
intermdiaire de condenst du module 1 vers les bacs de stockage (S-204A/B/C) ou irectement
vers CSTF .
18
Schma II .2 : schma de Section de stabilisation et fractionnements.
II.5.2 -Unit de rgnration du glycol
La formation dhydrates peut tre vite en se plaant en dehors des conditions

thermodynamiques de formation, donc il est ncessaire pour viter la formation des
hydrates soit de rduire la teneur en eau du gaz par une opration de schage, soit dutiliser
des inhibiteurs.
Ces derniers agissent comme des antigels permettant dabaisser la temprature de
formation des hydrates. Les alcools, notamment les glycols sont trs largement utiliss comme
inhibiteurs.
Les glycols prsentent lavantage de pouvoir tre facilement rcuprs en phase liquide et
rgnrs aprs chaque utilisation.
19
II.5.2 .1 : Description du Procd
Le DEG (Dithylene Glycol) hydrat provenant des sparateurs (glycol condenst) des
trois trains avec une concentration de (76%) est dtendu une pression de 1,5 bars dans le
ballon de dtente V-214 sous le contrle de la vanne de niveau LV-8 ,ensuite stock
dans le bac dalimentation de lunit de rgnration S-201.
Le DEG hydrat est refoul par la pompe de charge P206. Le glycol hydrat avant son
entre dans les 2 units de rgnration de glycol travers les Prfiltres de Glycol F-201 A/B
Le DEG hydrat et pr-filtr est prchauff 40C par lchangeur de chaleur
E-301 (cote tube) avec le DEG pauvre chaud est ensuite envoy vers un sparateur trphasique
V-301(P=2.5kg /cm2 T=40C), o les gaz dissouts et
les hydrocarbures liquides entrains sont retirs respectivement par flash et par dcantation.
Le DEG (Dithylene Glycol) hydrat quitte le sparateur et scoule travers un filtre glycol
F-301A/B pour liminer les particules solides qui pourraient causer du moussage et
augmenteraient alors les pertes en DEG.
Avant dentrer dans la colonne de rgnration le DEG hydrat est chauff jusqu
approximativement 92C grce lchangeur de chaleur E-302A/B (cote tube) avec le DEG
pauvre chaud travers un second changeur de chaleur glycol/glycol.
Le DEG riche scoule vers le bas de la colonne T-301 et contrl par la vanne LV102 en
rencontrant le DEG chaud et la vapeur deau qui remontent dans la colonne. Ces vapeurs sont
mises dans le rebouilleur H-301 et sont spares du DEG dans la colonne de rgnration. Les
vapeurs deau sont refroidies et condense dans larorfrigrant EF-301, puis entrent dans le
ballon de reflux V-302. Environ 25% de cette eau est renvoye en tte de la colonne par la pompe
(P-302A/B) et lautre partie de ce distillat envoye vers le puisard des eux glycol 14-S208-B.
La temprature du rebouilleur est maintenue 120C. Cette temprature permet de retirer
leau du glycol riche et dobtenir une puret en glycol de 85% en poids.
Le liquide rgnr scoule ensuite par gravit dans lchangeur de chaleur chaud
glycol/glycol E-302A/B (cot calandre) chaud puis dans lchangeur de chaleur froid
glycol/glycol E-301(cot calandre). Enfin il est pomp par les pompes P-301A/B de glycol vers
larorfrigrant E-215 sous le contrle de la vanne de niveau LV-301 et stock dans le bac de
stockage du glycol rgnr S-202. Linjection du glycol dans les trois trains est assure par la
pompe volumtrique P- 210 3A/1C/2C/3C piston avec une pression de refoulement de 115
bars. [6]
20
Chapitre III Gnralits sur les produits

III. 1 .GPL
III.1 1 Dfinition du GPL

Le GPL gaz du ptrole liqufi est un mlange gazeux constitue essentiellement de
propane, et de butane, et un peu de traces de mthane, de lthane et le pentane.
Le gaz du ptrole liqufi est un produit gazeux la temprature ambiante et pression
atmosphrique, mais il peut demeurer ltat liquide sous des pressions relativement basses (4
-18 bars). Pour faciliter son transport, sans stockage et sa commercialisation, il se gazifie au
moment de son utilisation.
La composition chimique de GPL est variable selon les normes et les utilisations dans
diffrents pays. Il peut contenir le propylne, butne, une faible quantit de mthane, dthylne,
de pentane et exceptionnellement des hydrocarbures tels que le butadine, lactylne et le
mthyle-actylne [7]
III 1.2 -Sources du GPL

Les gaz de ptrole liqufis (GPL) sont produits principalement :
Dans les raffineries de ptrole brut, soit au cours de la distillation du ptrole ,soit
pendant le craquage thermique ou reforming catalytique des produits en vue de la production
des essences.
Dans les units de traitement et de sparation du gaz naturel (GN) qui a pour but de
sparer et de recueillir les condenst (propane, butane, essence lgre etc..).
Par liqufaction du gaz naturel.
Tableau III.1 : Composition molaire du GPL [2]
Les composants du GPL % Molaire

Mthane 0,32
Ethane 1,12
Propane 60,95
Iso butane 15,46
Normal butane 22,14
Iso pentane 0,01
21
III.1 .3 : Caractristiques gnrales des GPL
1 -Odeur : le GPL est inodore ltat naturel, mais pour des raisons de
scurit on doit ajouter des odorants (composs sulfures tel que le
dithylmercaptan ou le dimthylsulfide); afin de dtecter les fuites.
2 -Tension de vapeur : la TVR du GPL est de 8 et 2 bars pour le propane et le
butane respectivement 20C.
3 -Dilatation : ltat liquide, le GPL a un haut coefficient de dilatation dont il faut
tenir compte lors de leur stockage (les sphres ne doivent jamais tre compltement remplies).
4 -Densit : aux conditions normales de temprature et de pression, les GPL sont plus
lourds que lair, elle diminue avec laugmentation de la temprature.
5 -Temprature dbullition : la pression atmosphrique la temprature dbullition de
propane est de ( 42C), celle de butane est de ( 6C).
6 -Pouvoir calorifique : cest la proprit la plus intressante tant donn que le GPL
est traditionnellement utilis pour les besoins domestiques :
Iso-butane: PCI = 29460 (kcal/kg). Normal butane : PCI = 29622 (kcal/kg). Propane :
PCI = 22506 (kcal/ kg).
7 -Impurets : le plus important est le soufre, la teneur en soufre est infrieure ou
gale 0,005 % en masse, aussi leau est lune des impurets importantes.
8- Le GPL est non corrosif lacier Le GPL na aucune proprit de lubrification et cette
caractristique doit tre prise en considration lors de la conception des quipements du GPL
(pompes et compresseurs) Le GPL est incolore, que ce soit phase liquide ou gazeux.
III .1. 4 : GPL en Algrie
Durant la dernire dcennie, lindustrie Algriennes du GPL a connu des changements

profonds, notamment en termes de production, dexportations et le transport maritime.
Le programme de valorisation des ressources gazires, lanc au dbut des annes 90 fait
bnficier aujourdhui SONATRACH des disponibilits importantes de GPL.
Depuis la mise en exploitation des champs HAMRA et OHANET, la production des GPL
a suivi une croissance soutenue, La production est tombe 8,4 millions de m3 en 2006. Avec
la mise en service de la nouvelle usine de sparation ARZEW prvue en 2010, les
disponibilits lexportation devraient atteindre 10 millions de tonnes. [8]
22
III.1.5 : La demande nationale de GPL
Le niveau de la demande nationale de GPL est de lordre 1,4 millions de

tonnes (sources NAFTAL) dont 90% de butane, 5% de propane et 5% de GPL carburant
(GPL/c).
Elle est rpartie comme suit :
336.000 tonnes dans la rgion Ouest (24%)
516.000 tonnes dans la rgion Est (37%)
Et en fin 401.000 tonnes dans la rgion centre (29%)
Figure III .1 Demande nationale en GPL [8]
Figure III .2 Offre prvisionnelle en national de GPL (en MT) [9]
23
III.1.6 -Domaines dutilisation des GPL
III.1.6.1 -GPL carburant
Suite une phase dtude et dexprimentation entame en 1977, la dcision dintroduire

le GPL carburant SIRGHAZ est intervenue en 1983 avec ladoption de la bicarburation et
la mise en place de rglementation lie aux conditions dutilisation et de distribution du
GPL/C.
III.1.6.2 -GPL dans la production lectrique
Il est utilis pour la gnration dlectricit, pourrait tre important si les conditions de
nature conomique venaient tre favorables, les marchs cibls cet effet, et qui sont
attractives pour le GPL.
La production dlectricit en Algrie seffectue uniquement par le gaz naturel. Les GPL
nont pas t introduits jusqu ici dans la gnration dlectricit ou dans la combustion.
III.1.6.3 -GPL dans les mnages

Le butane et le propane qui sont obtenus par la sparation du GPL ont une grande
importance dans les mnages. Le butane est utilis principalement dans le secteur domestique
pour la cuisine et le chauffage, la consommation nationale reprsente lquivalent de 11 millions
de bouteille de 13 Kg en 1996.
Le propane reprsente 4,8% de la consommation nationale du GPL, il est utilis comme
combustible dans les units industrielles et artisanales
en substitution domestique dans la rgion non relies au gaz naturel.
III.1.6.4 -Le GPL dans la climatisation
Il est utilis pour la construction des rfrigrateurs et des climatiseurs (moyenne capacit)
Grce sa dtente dabsorption de chaleur et de crer le froid.
III.1.6.5 -GPL dans la ptrochimie
Les 18% de la production mondiale des GPL sont destins l'industrie ptrochimique
dans le but d'une utilisation dans le vapocraquage pour l'obtention de bases olfiniques et
aromatiques. Il existe cependant d'autres usages ptrochimiques des GPL.
24
III.1.7 -Stockage de GPL
Il est indispensable de penser, dores et dj, la valorisation des excdents de GPL Pour
ne pas bruler dans les torches ; il faut prvoir des moyens efficaces pour le stockage. Les GPL
sont stocks sous pression dans des rservoirs sphriques. Lors de leur stockage, des vapeurs de
GPL (boil off) se dgagent par bullition sous leffet de plusieurs paramtres notamment :
La convection et la radiation de la chaleur atmosphrique.
Lchauffement d la fraction GPL produit dans les canalisations.
La diffrence de temprature entre le liquide stock et celui introduire dans le

rservoir.
Pour viter la perte de ces vapeurs ; les paramtres de stockage doivent tre maintenus
dans une plage bien dtermine, pour cette raison, un systme de rfrigration par
compression, refroidissement et dtente des vapeurs sont utiliss. Ces vapeurs sont
comprimes de 5.3 bars 19 bars et 115C par un compresseur alternatif puis condenses
54C dans un aro-rfrigrant et enfin dtendues 5 bars et 15C. Ceci permet le
refroidissement de la sphre et par consquent la rduction de sa pression.
Il est noter que certaines units utilisent dautres techniques de stockage savoir :
stockage dans des sphres, des pressions variables entre 12.5 et 14.5 bars
stockage comme fuel gaz dans des cigares
Une extension de la capacit de stockage en surface pour recevoir de grandes quantits de
GPL excdentaire est carter cause de ses cots levs et lincertitude des quantits de GPL
stocker. Quant aux techniques de stockage souterraines en usage actuellement dans le Monde,
elles peuvent apparatre, difficiles envisager cause des investissements levs quelles
requirent.
Certains experts de SONATRACH proposent un procd de stockage qui consiste rinjecter le GPL
sous forme liquide au sein du mme gisement du gaz ou il est produit laide dun ou plusieurs puits
injecteurs. Cela devient une possibilit de stocker de grands volumes avec un investissement en plus
la possibilit de dstocker n importe quel moment. [10]
25
Figure. III .3: Sphre de stockage [10]

III.2 : Condenst
Le condenst est constitu dun mlange dhydrocarbures paraffiniques, de liso pentane
jusquau dodcane et plus, gnralement extrait des gisements de gaz condenst, il est sous
forme liquide dans les conditions normales de temprature et pression.
III.2.1 : Caractristiques du condenst
Les principales caractristiques physiques du condenst sont :

Aspect : cest un liquide incolore, avec une odeur dessence.
Poids spcifique : il est compris entre 0.7 0.8.
Point dclair : infrieur 40C.
Limite dinflammabilits : il est trs inflammable, car il a un point dclair infrieur
0 , ses limites dinflammabilit sont approximativement 1.4 7.6%vol (dans lair).
Densit de vapeur : les vapeurs du condenst sont plus lourdes que lair, la densit
de vapeur est de 3 4 (dans lair).
Explosivit et inflammabilit : le condenst est un fluide hautement inflammable et
vaporable temprature et pression normale, du fait quil nest pas lectro-conducteur,
26
prsente un danger dinflammation ou explosion du fait de ltincelage lectrostatique amorc

par coulage, filtration, chute, pulvrisationetc.
On devra faire attention car les vapeurs du condenst constituent un mlange gazeux explosif
se rpondant par terre du fait de sa densit plus leve que celle de lair.
Toxicit physiologique : les vapeurs du condenst sont toxiques, lorsque un homme
sy expose, le premier symptme constat sert lirritation des yeux ce qui sera suivi de
symptmes nvropathies, lorsque la concentration des vapeurs du condenst sont de lordre de
0.025 0.05% vol.
Dans lair, elles ne peuvent pas causer de symptmes srieux, mme aprs plusieurs heures
dinhalation.
Prcautions hyginiques : pour prvenir lintoxication, on doit raliser une
ventilation convenable des locaux de travail et maintenir la concentration des vapeurs du
condenst moins de 300 PPm. [11]
III.2.2 : Utilisation du condenst
La particularit du condenst rside dans son utilisation dans deux secteurs industriels
entirement spars, le raffinage et la ptrochimie.
Cette bivalence confre ce produit une valeur diffrente par secteur dutilisation, qui
reprsente en fait un cot dopportunit par rapport la charge classique, dans une raffinerie
ou un craqueur ptrochimique.
En plus de cette diversification dutilisation finale, le march de destination intervient
galement dans la dfinition de ce cot.
La valorisation de ces condensats, se fait principalement dans le domaine de la
ptrochimie, secteur ou la demande est trs soutenue. Le vapocraquage est ainsi la principale
voie de valorisation pour la production dolfines.
Dans ce secteur, plusieurs lments dominent :
laccentuation du dsquilibre entre les demandes respectives dthylne et de
propylne, avec une croissance moyenne annuelle voisine de 4% pour le premier
et proche de 5% pour le second.
Le fonctionnement (de plus en plus pratiqu) des vapocraquages de charges
liquides la limite de flexibilit du rapport propylne/thylne.
Laccroissement de la contribution du craquage catalytique la fourniture de
propylne.
Lintrt de ce produit dans la ptrochimie est justifi par les considrations
suivantes :
27
Le condenst est riche en hydrocarbures paraffiniques et naphtniques en C5 etC7

ayant un bon potentiel de production dolfines (thylne et propylne).
Les hydrocarbures aromatiques, qui se transforment en polyaromatiques sans
valeur commerciale sont gnralement prsents en faibles quantits.
La teneur en soufre observe en moyenne ne gne pas le vapocraquage, le
condenst est utilis sans prtraitement ce titre.
Et gnralement leur prix de march est favorable en comparaison avec dautres
charges similaires, particulirement dun naphta de distillation atmosphrique.
III.3 -gaz sec (gaz de vente)
Un gaz sec ne forme pas de phase liquide, dans les conditions de production cest dire que
les points reprsentants la condition dans le rservoir et ensurface se trouve tous deux en
dehors du domaine biphasique au cours de la production dun tel gaz il ne se forme jamais de
phase liquide, ceci reste applicable lorsque la pression dcrot dans le rservoir avec le temps.
Pendant lexploitation du gisement ; cette situation implique d un domaine biphasique
relativement troit, le gaz doit tre concentr en mthane et contenir trs peu dhydrocarbures
plus lourd que lthane.
III.3.1 -Caractristiques du gaz de vente
-La pression ; 71 bars 65C.

-Point de rose : -6C 80.5 bars effectifs.
- Teneur en C5 et plus lourd : 0.5%.
-Teneur en humidit 50PPm en volume max.
-Pouvoir calorifique ( 15C et 16 bars) :9350 9450 kcal/sm3.
-Densit : 0.66.
III.3.2 -Utilisation de gaz de vente
1- Conversion chimique :
Il est possible de transformer le mthane en produit liquide aux conditions ambiantes, le but
cest davoir un carburant, essence, krosne ou gas-oil facile transporter et utiliser.
Pour effectuer cette transformation, on dispose 2 mthodes :
- Mthode direct
Elle consiste essayer de transformer le mthane en hydrocarbures liquides sans la
production intermdiaire du gaz de synthse.
28
- Mthode indirect
cest la mthode la plus utilis et qui passe par la production du gaz de synthse, mlange
doxyde de carbone et dhydrogne, aprs lajustement de la composition de ce mlange.
2 -Production du gaz de synthse
CH4 + O2 CO + 2H2 H298 = -35.7 kj/mol.
CH4 + H2O CO + 2H2 H298 = -35.7 kj/mol. [12]
3 -Synthse de mthanol :
CO + 2H2 CH3OH H298 = 90.8 kj/mol.
CO + 2H2 CH3OH H298 = 49.5 kj/mol.
3 -comme gaz combustible
Il est obtenu soit partir du gaz dalimentation, soit partir du gaz de vente, ce dernier est
utilis pour le fonctionnement normal, le gaz dalimentation tant disponible pour le dmarrage,
la distribution du gaz dutilit est assure par un systme de moyenne pression (15 kg/cm2) et
un systme basse pression (5.5 kg/cm2).
- Gaz dutilit moyenne pression alimente le module IV et le CSTF.
- Gaz dutilit basse pression :
Rebouilleur du dethaniseur.
Rebouilleur de dbutaniseur.
Rebouilleur du glycol.
Les torches hautes et basse pression.
29
30
Chapitre IV
Thorie de la rectification

Chapitre IV Thorie de la rectification
IV.1-Thorie de la rectification
IV.1.1 Dfinition
La rectification est un procde physique de sparation dun mlange dhydrocarbures
et cela par des changes de matire et dnergie entre deux phases circulant contre-
courant dans un appareil appel colonne de rectification. [13]
IV.1.2 -Description de la colonne de rectification
La colonne de rectification est un appareil cylindrique vertical, qui peut tre constitu
des plateaux successifs ou de garnissages, dont le rle consiste assurer un contact intime
entre les deux phases.
Dans chaque plateau il ya barbotage de la vapeur dans la phase liquide et ceci
permettra le transfert de matire et dnergie entre les deux phases, pour cela les plateaux
seront munis dlments spciaux comme les calottes. Le liquide se dplace par gravit de
haut en bas, et la vapeur de bas en haut, grce lnergie de pression utilise pour quil
yest un bon barbotage.
Le liquide arrivant au fond de la colonne est partiellement vaporis dans le rebouilleur
qui cre la phase vapeur ncessaire la distillation, le liquide non vapor est extrait du
fond de la colonne et constitue le rsidu.
La vapeur arrivant en tte de la colonne est condense, une partie du produit condense
est envoye en tte de la colonne en qualit de reflux, le reste est soutir en tant que distillat
(produit de tte)
La colonne est devise en trois zones :
Zone dalimentation : cest lentre du mlange fractionner.
Zone de rectification : se situe au-dessus de la zone dalimentation.
Zone dpuisement : se situe au-dessous de la zone dalimentation.
IV.1.3 -Principe de fonctionnement dune colonne de rectification
Le principe de fonctionnement dune colonne distiller est simple, tout en travaillant
une pression constante, il consiste dplacer les quilibres de phases laide dun gradient
de temprature cre par une source froide (condenseur), qui cre un flux de liquide froid
descendant dans la colonne, une source chaude (le rebouilleur), gnrant un flux de vapeur
chaud montant.
30
Schma IV.1 : Schma gnral du fonctionnement dune colonne de

rectification.
Pour assurer les transferts de matire et de chaleur entre ses deux flux, des lments de
contact sont placs lintrieur de la virole.
Ils sont constitus soit par des plateaux, dont le principe de base est de faire barboter la
vapeur dans une rtention liquide provoque par un barrage coupant la phase liquide, soit par
des garnissages constitus le plus souvent par un treillis mtallique dispersant les deux
phases et assurant une bonne surface dchange entre les fluides circulant contre-courant.
Le rebouilleur fonctionne laide dun fluide extrieur (vapeur, huile chaude) le plus
souvent refroidi par lair ou de leau.
Lorsque la colonne fonctionne en continu, la charge est injecte en un point de la virole
dpendant des caractristiques de la charge, celle-ci est alors spare en rsidu au fond de la
colonne et en distillat en tte.
La zone dpuisement a pour objectif dappauvrir le rsidu en produits lgers, la zone de
tte ou zone de rectification, a pour objectif dliminer les produits lourds du distillat, cette
dfinition implique que la zone dpuisement pourra tre moins efficace que la zone de
rectification si les exigences de qualit sont plus faibles pour le rsidu que pour le distillat.
31
IV.2 -Lois fondamentales de distillation

Lois de DALTON et de RAOULT
- La pression dquilibre dun mlange gazeux suppos parfait est la somme des
pressions partielles de ces constituants.
PT=Pi (IV.1).
- En appliquant les lois des gaz parfaits chacun des constituants i dans la phase
gazeuse on peut crire :
yi=Pi/ PT (IV.2).
- Daprs la loi de RAOULT la pression partielle dun constituant i dans une
solution idale est gale au produit de la tension de vapeur de ce constituant par sa
concentration molaire en phase liquide.
Pi=Pi.xi (IV.3).
Avec Pi : La tension de vapeur dun constituant i.
- Dans le cas de mlange rel la relation (IV3) devient comme suit :
Pi=Ei.xi (IV.4).
O : Ei: coefficient dHENRY.
- Equilibre liquide-vapeur:
On dit que les deux phases sont en quilibre lorsque la pression partielle de chacun des
constituants est la mme dans les deux phases (liquide et vapeur), Des relations (VI.2) et
(VI.3) on peut faire apparatre le coefficient dquilibre :
Ki=Pi/ PT =yi/xi (IV.5).
Le coefficient dquilibre ainsi dfini ne peut tre que dun usage restreint tant donc
- La dfinition des pressions partielles dans la phase vapeur implique que chacune de
celle-ci se comporte comme un gaz parfait ce qui ne peut tre vrifi qua basse pression.
- Lorsque la pression augmente, la valeur de Ki calcule comme gale au rapport
(Pi/), ne caractrise pas la distribution relle des constituants entre les deux phases en
quilibre. Dans ce cas, on utilise les valeurs exprimentales de Ki ou celle calcules par la
mthode thermodynamique (diagramme JENNY SHEIBEL).
Dtermination de la temprature de bulle tb et de rose tr :
1) Taux de vaporisation molaire :
Considrons un mlange de n constituants sous une pression P, sachant que :
NL : Nombre de Mole de lalimentation dune colonne.
ND: Nombre de Mole de la phase vapeur (distillat).
NR : Nombre de Mole de la phase liquide (rsidu).
xLi: Concentration molaire du constituant i dans la charge.
32
NLxLi=NDyi+NRxi (IV.6)
Sachant que = ND /NL
Par substitution on aura :
xi =xLi /[1+ (k i- 1)] = 1 (IV.7)
yi =kixLi /[1+ (ki 1)] =1 (IV.8)

Cela en combinaison avec la formule (IV5).
2) Dtermination de la temprature de bulle tb :
A la temprature de bulle, le mlange se trouve compltement ltat liquide (=0),
et substituant la valeur de dans lquation (IV7) nous aurons :
Ki xi= 1 ( IV.9)
La temprature de bulle est dtermine par approximations successives.
3) Dtermination de la temprature de rose tr :
A cette temprature le mlange se trouve comptemment ltat vapeur (=1), nous aurons :
yi/Ki=1 (IV.10)
La temprature tr est dtermine par approximations successives. [14]
IV.3 -Calcul dune colonne de rectification
Pour le calcul dune colonne de rectification on doit passer par les tapes suivantes
IV.3.1 -Bilan matire de la colonne
Le bilan matire pour la colonne entire et pour llment i scrit comme suit :
L=D+R (IV.11).
L.xL,I=D.yD,i+R.xR,I (IV.12).
Lquation obtenue sert calculer le rsidu (R) et distillat (D), si la concentration du
constituant i dans le distillat et dans le rsidu est connue.
Avec :
L : Le dbit molaire de la charge [kmole/h].
R : Le dbit molaire du rsidu [kmole/h].
D : Le dbit molaire du distillat [kmole/h].
xL : concentrations molaires du constituant le plus volatil dans la charge.
yD: concentrations molaires du constituant le plus volatil dans le distillat.
xR: concentrations molaires du constituant le plus volatil dans le rsidu.
33
a) Bilan matire de la zone de rectification

Le taux de reflux rf
Le reflux froid a pour but de rguler la temprature en tte de colonne, de corriger la
composition et de ne permettre quaux lments lgers de quitter lappareillage, le reflux
froid prend la chaleur des vapeurs et provoque une condensation des constituants lourds.
Il est gal au rapport de la quantit de reflux
renvoy dans la colonne sur la quantit du distillat :
rf=g/D (IV.13)
Schma IV.2 : schma du tte de la colonne
On donne les significations suivantes :

G, y : Le dbit et la concentration molaire de la vapeur.
g, x : Le dbit et la concentration molaire du liquide.
D ,y : Le dbit et la concentration molaire du distillat.
Hi,hi : Enthalpie respectives de la vapeur et du distillat (kcal/kg) .
Q1 : Quantit de chaleur soustraite du distillat (kcal/kg).
Le bilan matire pour la zone de rectification pour le constituant le plus volatil i scrit
:
G=g+D (IV.14)
G.yn,i=g.xn-1,i+D.yD (IV.15)
A partir des quations (13), (14) et (15), on trouve lquation de concentration de la
zone de rectification : Yn,i=m.xn-1,i+(1-m).yD,i (IV.16)
Avec : m = g / G = rf/(rf +1)
34
Schma IV.3 : Schma descriptif du fonctionnement de la zone de rectification

.
b) Bilan matire de la zone dpuisement

1 -Le taux de rebouillage rb (reflux chaud)
Le reflux chaud a pour but dapporter la chaleur, de maintenir la temprature au fond
de la colonne et dassurer la vaporisation des lgers entrans par le produits de fond.
Il est gal au rapport de la quantit de reflux chaud provenant du rebouilleur la quantit
de rsidu :
rb=GR/ R (VI.17).
2-Bilan matire de zone dpuisement
Schma IV.4 : schma du fond de la colonne.

Admettant que :
R, xR : dbit et concentration molaire de rsidu.
Q2 : Quantit de chaleur fournie par le rebouillage (Kcal)
35
Le bilan matire pour la zone dpuisement et pour le constituant i le plus volatil scrit :
g=G+R (IV.18)
g.xn+1,i=G.yn,i+R.xR,I (IV.19)
La rsolution de lquation (IV.19) par rapport Xn+1 ,i permet dobtenir lquation de
concentration de la zone dpuisement :
xn+l,i=yn,i/m+((m-1)/m)xR,i (IV.20)
Avec : m=g/G=(rb+1)/rb (IV.21)
Schma IV.5 : Schma descriptif du fonctionnement de la zone dpuisement.
3 -Bilan matire de zone dalimentation

Il y a lieu de procder au bilan de matire de la zone dalimentation. Nous avons donc le
bilan de matire effectu juste au-dessous du plateau suprieur de cette zone.
gm=gk+Gd=GK+R (IV.22)
gmxm,i=GkyK,i+RxR,I (IV23)
36
Schma IV.6 Schma descriptif du fonctionnement de la zone

dalimentation
La vapeur Gm alimentant la zone de rectification est la somme de vapeur GK et de la

phase GD, donc le bilan de matire effectu juste au-dessous du plateau infrieur de la zone
de rectification donne :
Gm=GK+GD=gk+D (IV.24)
Gmym,i=GKyK,i+GDxD,I (IV.25)
Dans le cas de rectification des mlanges complexes, nous serons amen effectuer
des sries de calculs dquilibre sur chaque plateau. Pour viter le travail long et laborieux
des approximations successives, nous pouvons utiliser une mthode simplifie rapide
donnant une prcision souvent suffisante. Cette mthode fait intervenir la notion de volatilit
i tant le rapport des coefficients dquilibre des constituant i et le coefficient
dquilibre du constituant cl r.
i=ki/Kr (IV.26)
En qualit de constituant cl, nous pouvons choisir nimporte quel constituant
(habituellement le constituant le moins volatil du mlange).
En tenant compte de (IV.26), lquation dquilibre peut tre crite sous la forme :
yi=KiX=iKrxI (IV.27)
xi=y/iKr (IV.28)
Dans la zone de rectification : i,m=0.5 (D,i + L,i) (IV.29)
Dans la zone dpuisement : i,m =0.5 (R,i + L,i) (IV.30)
37
IV.3.2 - Les paramtres opratoires

a) Choix de la pression dans la colonne
Le choix de la pression dans la colonne dpend essentiellement de la composition et
de la nature physique et chimique des mlanges sparer quelle que soit leur complexit.
A faible pression, la slectivit de la sparation est meilleure, car on peut viter
laltration thermique des produits, mais en contre partie, faible pression, il faut condenser
basse temprature les vapeurs de tte pour obtenir la concentration voulue, donc cest le
choix du fluide rfrigrant qui va imposer la pression de la colonne.
Gnralement, nous utilisons une source froide bon march et abondante (eau et air)
qui permet dobtenir un condenst de 30 50 C environ. Si la temprature de condensation
des vapeurs de distillat la pression atmosphrique est plus basse que celle de leau ou de
lair, la colonne doit fonctionner sous haute pression pour que la condensation des vapeurs
du distillat puisse tre ralise par rfrigration ordinaire et leur temprature de condensation
doit tre prise au moins de 15 20 C plus grande que celle du fluide rfrigrant quelle que
soit la temprature de service. Ainsi, la pression Pb dans le ballon de reflux est dtermine
par approximations successives daprs lquation disotherme de la phase liquide :
(Ki.XD,i)=1 (IV.31)
Les valeurs de Ki sont dtermines daprs le diagramme de JENNY et SCHEIBEL
(Annexe 1)
Pression au sommet de la colonne
PS=Pb+P1 (IV.32)
Avec : Les pertes de charge au fond de colonne qui sont compris entre
P1 = 0.2 et 1 atm.
Pression au fond de la colonne
Pf=PS+P2 (IV.33)
Avec : Les pertes de charge au fond de colonne qui sont compris entre
P2 = 0.3 et 0.6 atm.
Pression dans la zone dalimentation
PA=0.5 (PS+Pf) (IV.34)
b) Choix de temprature dans la colonne
A la pression de service PS (Pression au sommet de la colonne) et en utilisant
lquation disotherme de la phase vapeur, nous dterminons par approximations successives
la temprature au sommet de la colonne.
(yD,i/Ki)=1 (IV.35)
38
La temprature au fond de la colonne est dtermine par lapplication de lquation de

lisotherme de la phase liquide (par approximations successives).
(Ki.xR,i)=1 (IV.36)
Pour la temprature de la charge, trois cas peuvent se prsenter :
Si la charge est ltat liquide, nous utilisons lquation suivante :
(Ki.xl,i)=1 (IV.37)
Si la charge est ltat vapeur, nous utilisons lquation suivante :
(yL,i/Ki)=1 (IV.38)
Si la charge est ltat liquide vapeur, nous utilisons lquation suivante :
x0,i=(xL,i/(1+e()Ki1))=1 (IV.39)
x0,i=y0,i/Ki (IV.40)
IV.3.3 -Taux de vaporisation et composition des phases liquide et vapeur de la charge
Lors du calcul de la zone dalimentation, il est ncessaire de disposer du dbit et de la
composition des phases liquide et vapeur de la charge.
Daprs la dfinition du taux de vaporisation (ou de condensation) de la charge, nous avons :
e0=G0/L o G0=e'0.L (IV.41)
r0=g0/Ld o g=r0.L (IV.42)
A une pression et temprature donne, le taux de vaporisation ainsi que la composition de
la charge sont dtermins par approximations successives, daprs lquation connue de
TREGOUBOV :
x0,i=(xL,i/(1+e0(Ki1))=1 (IV.43)
O:
e0+r0=1 (IV.44)
IV.3.4 -Fonctionnement reflux
Le fonctionnement reflux minimal correspond un nombre infini de plateaux pour la
sparation spcifie, donc daprs la formule dUNDERWOOD :
(rf)min=(yD,i/(i, L- ) (IV.45)
(rb)min=[(i,LxR,i)/(i,L-)] (IV.46)
V.3.5 : Taux de reflux optimal
A partir du taux de reflux opratoire choisi en tte, ltude du procd de rectification
montre qu chaque valeur du taux de reflux correspond un nombre fini de plateaux
thoriques et que ce nombre diminue quand le taux de reflux saccrot.
Pour ces calculs approximatifs, nous pouvons utiliser la formule propose par
GILLILAND :
[((rf)op(rf)min)/((rf)op+1)]=[0.1-O.33] (IV.47)
39
Pour le rapport moyen, non aurons :

(rf)op=1.3(rf)min+0.30 (IV.48)
V.3.6 - Bilan thermique de la colonne
Ltablissement du bilan thermique de la colonne est ncessaire pour vrifier le dbit du
reflux liquide provenant du condenseur et pour dterminer les charges thermiques du
condenseur et du rebouilleur.
Lquation du bilan thermique de la colonne entire scrit :
L.hL+QR=D.hD+R.hR+QC (IV.49)
hL=e0.H0+(1-e0).h0 (IV.50)
QR=GR.(HRhR)=GF.hF (IV.51)
QC=(D+gX).(HDhD)=Ga.(H2 H) (IV.52)
Avec: hF = CPF. T et CPF= CPi . Wi

Les enthalpies de la phase liquide et la phase vapeur sont dtermins respectivement partir
de lannexe 2 et lannexe 3.
En apposant lquation de la zone dpuisement; nous aboutissons a la composition du
liquide descendant du dernier plateau:
X1,i=(1/m)yR,i+(m-1/m)xR,i (IV.53)
Sachant que :
yR,i=imR.xR,i/(im.XR,i) (IV.54)
Et
m= (rb)opt +(1/ (rb)opt) (IV.55)
Avec :
1515=420moy+5 (IV.57)
=0,0018480.00132.420moy (IV.58)
A partir du bilan thermique de la colonne, nous avons:
gx=(rf)opt.D (IV.59)
GR=(rb)opt.R (IV.60)
Taux de rebouillage optimal dfini par :
(rb)opt=[(D(rf)opt/L)+(1-e0)-R/L]/(R/L) (IV.61)
40
V.3.7 -Calcul du nombre de plateaux par la mthode de GILLILAND

GILLILAND: propose une corrlation empirique unique, entre le taux de reflux et le
nombre dtages thoriques pour les colonnes distillant des mlanges dhydrocarbures
lgers.
Si Nmin est le nombre minimal dtages thoriques reflux total, (rf) min le taux de reflux
minimal le taux de reflux pour un nombre infini dtages.
N : un nombre dtages finis associs un taux de reflux fini, sont runies entre elles par
une corrlation traduite graphiquement . [15]
F(N)= (N Nmin )/(N+1) et F (rf)= (rf- (rf) min) / (rf+1)
Le nombre minimal de plateaux thoriques est dtermin comme suit:
Dans la zone de rectification
N R min =log [(y'D,v x'L,L)/( y'L,L x'R,V)] / log V,L
Dans la zone dpuisement
N E min =log [(x'L,v x'R,L)/( x'L,L x'R,V)] / log V,L
Le nombre thorique de plateaux thoriques est dtermin daprs le formule de
GILLILLAND :
NENE min/ (NE +1) =0,75[1-((r-r min) /(r+1) 0,5688]
V.3.8- Calcul des dimensions de la colonne
a) Calcul de diamtre de la colonne

Le calcul du diamtre de la colonne se fait partir de la formule ci aprs :
b) Calcul de dbit volumique des vapeurs qui quittent la tte de colonne

Le dbit des vapeurs quittant la tte de colonne T202 est valu selon lquation suivante :
Le facteur de compressibilit Z est dtermin partir de lannexe 4, compte aux pressions

et tempratures critique de chaque constituant sont dtermin partir du DATA BOOK .
c) Calcul de la vitesse admissible
La vitesse admissible des vapeurs note W est calcule daprs la formule suivante :
41
V.3.9 -Hauteur de la colonne

La hauteur de la colonne est la somme des distances suivantes :
=1 +2 + 3 +4 + 5 + 6
= ( 1)
3 = 3
4 = ( 1)
Sachant que: = 0,6 m
42
Chapitre v Partie de Calcule
V.1 -Bilan de matire de la colonne T202

V.2 -Dbit et composition du distillat et du rsidu
V.3 -Pressions dans les diffrentes zones de la colonne
V.4 -Volatilits relatives
V.5 -Taux de vaporisation et composition des phases vapeur et liquide
dalimentation
V.6 -Calcul du rgime minimal
V.7 -Bilan thermique de la colonne
V.8 -Dimensionnement de la colonne
V.9 Partie simulation

Chapitre V Parties de Calculs
V.1 -Bilan de matire de la colonne T202

Donnes de dpart
Charge dalimentation =rsidu de dthaniseur L = 974.7 kmole/h.
Masse molculaire M = 77,9736 kg /kmole.
Temprature dentre T=165C.
La composition molaire de la charge est donne dans le tableau .
Tableau V.1 : Composition molaire de la charge dalimentation dbutaniseur.
Constituants xL, i Mi xL, i. Mi

CH4 0 16 0
C2H 6 0,0033 30 0,099
C3H 8 0,2791 44 10,30
iC4H10 0,0673 58 4,4718

nC4H10 0,1228 58 8,8624
IC5H12 0,0773 72 4,4856
n C5H12 0,1051 72 5,7672
C6H14 0,0836 86 6,3296
C7H16 0,0833 100 8 ,33
C8H18 0,0298 114 5,1072
C9H20 0,0324 128 4,1472
C10H22 0,0326 142 4,6292
C11H24 0,0299 156 4,6644
C12H26 0,0534 170 15,878
TOTAL 0,9999 - 77,9736

xL,i = xL,i M= Mi . xLi
Les produits du dbutaniseur (GPL, condenst stabilis) doivent rpondre aux spcifications
suivantes :
- Teneur en C2- 3%.
- Teneur en C5+ 0 ,4%.
- Condenst stabilis : Cest les C5+, dont la TVR10 Psia.
43
V.2 -Dbit et composition du distillat et du rsidu

Pour la sparation du GPL et du condensat, on peut fixer :
- constituant cl volatil nC4H10
- constituant cl lourd iC5H12
Le dbit molaire dalimentation est de 974,7 Kmol/h
La concentration de iC5dans le distillat et la concentration de nC4dans le rsidu sont 0,004
et de 0,02 respectivement.
Equation de bilan de matire : L = D + R
Equation de bilan par rapport un constituant i :
L* xi,L = D * yi,D + R*xi,R
974,7* 0 ,0033= D* y2,D + R * 0.(01)
974,7* 0 ,2791= D* y3,D + R * 0 ....(02)
974,7* 0,0673= D* yi4,D + R * 0 ....(03)
974,7* 0,1228= D* yn4,D + (974,7D) * 0,02 ......(04)
974,7* 0,0773= D* 0 ,004+ (974,7D) * x 5,R ......(05)
974,7* 0,1051= D* 0 + (974,7D) * xn5,R ........(06)
974,7* 0,0836= D* 0 + (974,7D) * x6,R ......(07)
974,7* 0,0833= D* 0 + (974,7D) * x7,R ...(08)
974,7* 0,0298= D* 0 + (974,7D) * x8,R ...(09)
974,7* 0,0324= D* 0 + (974,7D) * x9,R ...(10)
974,7* 0,0326= D* 0 + (974,7D) * x10,R .....(11)
974,7* 0,0299 = D* 0 + (974,7D) * x11,R ....(12)
974,7* 0,0534= D* 0 + (974,7D) * x12,R .....(13)
La somme des quations de la 5me la 13me pour obtenir le distillat puis on injecte le
rsultat dans lquation de bilan de matire et on dduit les fractions, on aboutit aux
rsultats suivants :
>D = 451.8973 kmol/h
> R= 522.8026 kmol/h
- De la relation (1) on a : Y2,D = 0.0071
- De la relation (2) on a : Y3,D = 0.6019
- De la relation (3) on a : Yi4,D= 0.1451
- De la relation (4) on a : Yn4,D=0.2417
44
Des relations (5, 6,7, 8, 9, 10, 11, 12,13) on a :

x5i,R = 0.1406 x9,R = 0.0604
x5n,R= 0.1959 x10,R= 0.0607
x6,R= 0.1558 x11,R = 0.0557
x7,R= 0.1553 x12,R = 0.0995
x8,R=0.0555
La composition molaire du distillat et du rsidu ainsi que leurs masses molaires sont
regroupes dans le tableau suivant ;
Tableau V.2 -Fractions et masse molaire de distillat et de rsidu .
Constituants Distillat Rsidu
YD, i YD, i. Mi XR, i XR, i. Mi
CH4 0 0 0 0
C2 H6 0,0071 0,213 0 0
C3 H8 0,6019 26,4836 0 0
iC4H10 0,1451 8,4158 0 0
nC4H10 0,2417 14,0186 0,02 11,6
IC5H12 0,004 0,2880 0,1406 10,1232
n C5H12 0 0 0,1959 14,1048

C6H14 0 0 0,1558 13,3988
C7H16 0 0 0,1553 15,5300
C8H18 0 0 0,0555 6,327
C9H20 0 0 0,0604 7,7312
C10H22 0 0 0,0607 8,6194
C11H24 0 0 0,0557 8,6892
C12H26 0 0 0,0995 16,915
TOTAL 1.0000 49.419 0,9999 113,0386
45
V.3 -Pressions dans les diffrentes zones de la colonne

V.3.1 -Calcul de la pression et les tempratures aux diffrentes zones
La pression et la temprature au ballon de reflux
Il faut dterminer la pression du ballon de reflux our le GPL aprs sa sortie du
condenseur, en prenant la composition prcdent et out en sachant que toute la vapeur de tte de
colonne passe ltat liquide via le condenseur total.
Pour une temprature critique de la rgion de Hassi Rmel de 54C et via le nomogramme
de Jenny et Scheibel et ( lannexe 1 ) parit rations successives jusqu ce que KiXi soit gale
1,nous avons dtermin les coefficients dquilibre suivants:
Pour les distillats arrivant au ballon de reflux avec P =11,9atm et KiXi
Tableau V.3 :Valeurs des Ki du distillat liquide.

Constituant XD, i K i ( T= 54C) XD, i*K i
C H4 0 0 0
C 2H6 0,0071 3,9 0,0276
C 3 H8 0,6019 1,3 0,7824

C 4H10 0,1451 0,62 0,0986
C 4 H4 0,2417 0,52 0,1256

C 5H12 0,004 0,25 0,0010
Totale 0,9999 1,0352
La pression et la temprature au sommet de la colonne

PS = Pb + P1.
PS : pression au sommet de la colonne.
P1 : perte de charge due aux rsistances des conduites du condenseur.
On prendra P1 = 1,1atm.
PS= 11,9 + 1,1.
Donc : PS =13 atm.
Maintenant, passant la temprature de tte de colonne, les rsultats sont rcapituls
dans le tableau qui suit:
Nous trouvons une temprature de 62 C, et la somme Yi/Ki est gale 1,0084
46
Tableau V.4 :Valeurs des Ki du distillat vapeur .

Constituants: Y K Y i/ Ki
C H4 0 0 -
C2H6 0,0071 3,8 0,0018
C3H8 0,6019 1,52 0,3959
iC4H10 0,1451 0,75 0,1934
nC4H10 0,2417 0,6 0,4028
iC5H12 0,004 0,28 0,0142
Total 0,9999 - 1,0084

La pression et la temprature au fond de la colonne
PF = PS +P2
P2 : pertes de charge dues aux rsistances des plateaux
On prendra P2 = 0,5 atm. PF = 13 + 0,5 >PF = 13,5 atm .
Le fond de colonne est P =13,5atm, T =189 C et o KiXi =1,0004.
Les rsultats sont rsums dans le tableau suivant:
Tableau V.5 : Valeurs des Ki du rsidu .

Constituants: X K X R,i* Ki
nC4H10 0,0200 3,10 0,0600

iC5H12 0,1406 1,90 0,2671
nC5H12 0,1959 1,70 0,3330

C6H14 0,1558 1,15 0,1791
C7H16 0,1553 0,68 0,1056
C8H18 0,0555 0,44 0,0244
C9H20 0,0604 0,30 0,0181
C10H22 0,0607 0,20 0,0121
C11H24 0,0557 0,15 0,0083
C12H26 0,0995 0,08 0,0079
Total 0,9999 - 1,0156
47

La pression et la temprature a la zone dalimentation
PA =1/2 (PS + PF) > PA= 0.5 (13 +13,5) >PA=13,25 atm
Il ne nous reste plus que lalimentation de la colonne qui est une pression de
13,25atm et une temprature de 165 C, et dont les rsultats sont:
Tableau V.6 : Valeurs des Ki de la charge.

Constituants Xi Ki
C2 H6 0,0033 14
C3 H8 0,2791 6
iC4H10 0,0673 3,8
nC4H10 0,1228 3,3
iC5H12 0,0773 2,2
nC5H12 0,1051 2,2
C6H14 0,0836 1,25
C7H16 0,0833 0,8
C8H18 0,0298 0,5
C9H20 0,0324 0,35
C10H22 0,0326 0,25
C11H24 0,0299 0,15
C12H26 0,0534 0,01
V.4 -Volatilits relatives

Calculons maintenant les volatilits relatives des constituants aux diffrentes zones de
la colonne en se basant sur la formule :
i = Ki / Kr
Kr : Coefficient dquilibre de llment le plus lourd.
En qualit de constituant de rfrence, on admet le constituant cl lourd iC5H12. Et le
constituant cl lger nC4H10.
On admet le constituant cl lourd iC5 H12 comme lment de rfrence.
Les valeurs moyennes des volatilits relatives sont dtermines par les formules
suivantes :
Pour la zone de rectification : im =0.5 ( D , i + L, i )
Pour la zone d'puisement : im =0.5 ( R , i + L, i )

Les rsultats sont illustrs dans le tableau (V.7)
48
Tableau V.7 : Volatilits relatives dans les diffrentes zones de T202.

Constituants Alimentation Rectification Epuisement
Ki L ,i Ki D, i i, m Ki R, i i, m
CH4 37 23,12 0 - - 35 18,42 -

C 2H 6 12 6,36 3,8 13,5 10,5 13 6,84 6,6
C 3H 8 5 2,72 1,52 5,42 4,27 6,5 3,42 3,07
iC4H10 3 1,72 0,75 2,67 2,03 3,5 1,84 1,78
nC4H10 2,6 1,5 0,6 2,14 1,88 3,10 1,63 1,63
iC5H12 1,6 1 0,28 1 1 1,90 1 1

n C5H12 1,55 0,97 0 0 0,46 1,70 0,89 0,91
C6H14 0,95 0,56 - - - 1,15 0,60 0,59
C7H16 0,57 0,36 - - - 0,68 0,35 0,35
C8H18 0,3 0,22 - - - 0,44 0,23 0,20
C9H20 0,2 0,15 - - - 0,30 0,15 0,13
C10H22 0,15 0,11 - - - 0,20 0,10 0,09
C11H24 0,9 0,068 - - - 0,15 0,07 0,06
C12H26 0,06 0,004 - - - 0,08 0,04 0,03
V.5 -Taux de vaporisation et composition des phases vapeur et liquide dalimentation

Le taux de vaporisation e0 de l'alimentation et la composition de ses phases l'quilibre
( vapeur- liquide ) l'entre de la colonne sont dtermins laide de l'quation de
TREGUBOV.
X0, i = X'L, i (1+ e0 (Ki 1 ))
Y0,i = X0,i . Ki
Pour le calcul, on procde par ttonnement, cest--dire quon prend des valeurs diffrentes
du taux de vaporisation, et on calcule les concentrations molaires de la phase liquide X0.i et la
phase vapeur Y0.i dans la zone de flash jusqu lobtention de :
X0i =1 et Y0i = 1
Les rsultats des calculs sont donnes dans le tableau.
49
Tableau V.8 : Taux de vaporisation .
Constituants P =13,25atm T = 165C e0= 0,7345

XL, i Ki X0, i Y0, i
CH4 0 37 0 0
C2 H6 0,0033 12 0,0003 0,0043
C3 H 8 0,2791 5 0,0708 0,3543
iC4H10 0,0673 3 0,0272 0,0817

nC4H10 0,1228 2,6 0,0564 0,1467
IC5H12 0,0773 1,6 0,0536 0,0858
n C5H12 0,1051 1,5 0,0768 0,1153
C6H14 0,0836 0,95 0,0867 0,0824
C7H16 0,0833 0,57 0,1217 0,0693
C8H18 0,0298 0,3 0,0615 0,0184
C9H20 0,0324 0,2 0,0785 0,0157
C10H22 0,0326 0,15 0,0869 0,0130
C11H24 0,0299 0,9 0,0903 0,0081
C12H26 0,0534 0,06 0,1724 0,0103
TOTAL 0,9999 - 0,9979 1,0058
Daprs le tableau on trouve que e0= 0.7345 e0

, . ,
e0=

, _
: Constante qui varie entre les volatilits relatives des constituants cls lourds et le
plus
volatiles.
1 < < 1,612 Le calcul pour = 1.314 nous donne les rsultats dans le tableau
suivant, sachant que :
.,. ,
e0 =
.,
50
V.6 -Calcul du rgime minimal

V.6.1 -calcul de taux de reflux et taux de rebouillage
Tableau V.9 : Verification du taux de vaporisation .
.,. ,
Constituants X'L, i L, i L, i . X'L, i L, i- e0 =
.,
C2 H 6 0,0033 6,36 0,020 5,046 0,0040

C3 H 8 0,2791 2,72 0,7591 1,406 0,4816
iC4H10 0,0673 1,72 0,1157 0,406 0,2249
nC4H10 0,1228 1,5 0,1842 0,186 0,6416
iC5H12 0,0773 1 0,0773 -0,3141 -0,2461
nC5H12 0,1051 0 ;97 0,1051 -0,377 -0,2612

C6H14 0,0836 0,56 0,0468 -0,721 -0,0687
C7H16 0,0833 0,36 0,0299 -0,958 -0,0309
C8H18 0,0298 0,22 0,0065 -1,127 -0,0049
C9H20 0,0324 0,15 0,0048 -1,189 -0,0034
C10H22 0,0326 0,11 0,0035 -1,221 -0,0024
C11H24 0,0299 0,068 0,0020 -1,258 -0,0013
C12H26 0,0534 0,004 0,0002 -1,277 -0,0015
Total 0,7313
V.6.2 -Calcul de taux de reflux minimal
Tableau V.10 :Calcul de rf .
Constituants Yi L, i- Yi / ( L, i- )
C2 H6 0,0071 5,046 0,0014

C3 H8 0,6019 1,406 0,4280
Ic4H10 0,1451 0,406 0,3573
nC4H10 0,2417 0,186 1,2994
iC5H12 0,004 -0,3141 -0,0127

Total 2,0734
Donc (rf) min = 2.7245
51
V.6.3 -Calcul de taux de reflux optimal

le taux de reflux optimal est calcul par la formule propose par GILLILAND.
(rf) op = 1,3*(rf) min+0,3.
Donc (rf) op = 1,3 * 2,7245+0,3 > (rf) op = 3,8418
V.6.4 -Calcul de taux de rebouillage minimal

., . ,
(r b) min =
.,
Tableau V.11 :Calcul de rb (taux de rebouillage).
. , . ,
Constituants x'R, i L, i L, i- . ,
nC4H10 0,0200 1,5 0,186 0,3835

iC5H12 0,1406 1 -0,314 -0,5093
nC5H12 0,1959 1 -0,377 -0,5671
C6H14 0,1558 0,56 -0,721 -0,1057
C7H16 0,1553 0,36 -0,958 -0,0515
C8H18 0,0555 0,22 -1,127 -0,0143
C9H20 0,0604 0,15 -1,189 -0,0061
C10H22 0,0607 0,11 -1,221 -0,0038
C11H24 0,0557 0,068 -1,258 -0,0401
C12H26 0,0995 0,004 -1,277 -0,0042
Total 0,9999 -1,0775
Daprs le tableau : (rb) min=1,0775
V.6.5 -Calcul de taux de rebouillage optimal
(rb) opt = [(rf)opt(D/R)+(1-e)L/R] > (rb) opt =1,45

V.6.6 -Calcul du nombre de plateaux thoriques dans la colonne
a) Calcul du nombre de plateaux thoriques dans la zone de rectification
Pour effectuer ce calcul, nous devons calculer dabord le nombre de plateaux minimal
N min R =log[(YD,V XL,L)/( YD,L XL,V)] /logimD,V
52
N min R =log[(0.2417*0.0773 )/( 0.004*0.1228)] /log 1.88

N minR = 5,43
On : ((rf) op- (rf) min)/((rf) op +1) =(3,8418-2,7245)/ (3,8418+1)
> (rf) op- (rf) min/(rf) op +1 =0,23
Et daprs le diagramme de GILLILLAND:
> (NR - NminR) / (NR +1) =0,43
> NR= 10,29 donc NR= 10 plateaux
b) Calcul du nombre de plateaux thoriques dans la zone dpuisement
N minE =log[(X'L,vX'R,L)/( X'L,L X'R,V)] /log Vm,L
N minE =log[(0.1228*0.1406)/(0.0773*0.02)] /log 1.59
Nmin E =5,01
Et daprs le diagramme de GILLILLAND:
NE - NminE / NE +1 =0,43
> NE= 9,49
Le nombre de plateaux thoriques dans la zone dpuisement est de NE = 9 plateaux.
Le nombre total de plateaux thoriques dans la colonne est la somme des plateaux dans
les deux zones : Nth =NR +NE= 19 plateaux
En introduisant lefficacit =0,6 du plateau clapet, le nombre de plateaux rel dans la
colonne entire sera de:
NREEL= Nth / =31,66 NREEL= 32 plateaux
V.7 -Bilan thermique de la colonne
Ltablissement du bilan thermique de la colonne est ncessaire pour vrifier les dbits du
reflux liquide et dterminer les charges thermiques au niveau du condenseur et du rebouilleur.
Lquation de bilan thermique de la colonne entire s'crit :
QR +LhL=D HD +R hR + Qc
Zone dalimentation :
LhL = h0 (1- e0) + e0 H0
Condenseur :
QC = (D + gx) (HD - hD)
La quantit de chaleur Q apporte au fond de la colonne, peut tre dtermine partir du bilan
thermique labor pour la partie situe au dessous du dernier plateau
g1=GR+ R
53
g1 hL + QR = GRhr + R HR.
(GR + R )h1 +QR = GR HR+ R Hr
QR = GR(HR-h1) + R(hR-h1)
QR =D HD+R hR+ QC - L hL
L, D et R : dbits molaires respectivement de la charge, distillat et du rsidu (k moles/h).
QR: charge thermique du four (kcal/h).
Qc: charge thermique du condenseur (kcal/h).
hL: Enthalpie de la charge biphasique (kcal/kmol).
e0: taux de vaporisation de la charge (molaire).
H0: enthalpie de la phase vapeur de la charge.
h0: enthalpie de la phase liquide de la charge.
hD: enthalpie de distillat ltat liquide (kcal/kmol).
HD: enthalpie de distillat ltat vapeur (kcal/kmol).
hR: enthalpie de rsidu (kcal/kmol).
gx: dbit de reflux liquide (kmol/h).
GR: dbit de la vapeur provenant du rebouilleur (kmol/h).
V.7.1 -Calcul des diffrents dbits de reflux de la colonne
gx=(rf )opt. D
gx= 3,8418*451,8973 = 1736,099 >gx = 1736,099 (kmol/h).
Gr=(rb )opt. R
Gr =1,45*522,8026 =758,06377 >Gr=758,06377 (kmol/h).
g1= Gr+gx
g1 = 758,06377 + 1736,099 =2494,1628 > g1= 2494,1628 (kmol/h).
laide de lquation de la zone dpuisement; on calcule la composition du liquide
descendant du dernier plateau :
X'1, i = (1/m) YR, i + (m-1/m) X'R,I YR,i= i m R* X'R,i / ( i m R * X'R,i)
Sachant que : m = (rb) opt + (1/ (rb) opt) =1,45 + (1/ 1,45)
> m=2,14 1/m=0, 46 m-1/m=0,5327
54
Tableau V.12 : Composition du liquide du dernier plateau .
Constituant X'r,i imR X'r,i* i mR YR,i X1,i Mi Mi. X1,i

nC4 H10 0,0200 1,63 0,0326 0,0617 0,03016 58 1,7497
iC5 H12 0,1406 1 0,1406 0,2662 0,1712 72 12,3266

nC5 H12 0,1959 0,91 0,1782 0,3374 0,2303 72 16,5868
C6 H14 0,1558 0,59 0,0919 0,1740 0,1602 86 13,7793
C7 H16 0,1553 0,35 0,0543 0,1028 0,1424 100 14,2495
C8 H18 0,0555 0,20 0,0111 0,0210 0,0470 114 5,3684
C9 H20 0,0604 0,13 0,0078 0,0147 0,0492 128 6,3070
C10 H22 0,0607 0,09 0,0054 0,0102 0,0483 142 6,8718
C11 H24 0,0557 0,06 0,0033 0,0062 0,0436 156 6,8083
C12 H26 0,0995 0,03 0,0029 0,0054 0,0765 170 13,0186
TOTAL 0,9999 0,5281 1,000 0,9994 97,0663
TableauV.13 :Composition massique et molaire du gaz de tte et du rsidu.

Mi Y'i Y'i*Mi Yi X'r,i X'r,i *Mi Xr,i
30 0,0071 0,213 0,0043
44 0,6019 26,4836 0,5358
58 0,1451 8,4158 0,1702
58 0,2417 14,0186 0,2836 0,0200 11,6 0,1026
72 0,004 0,2880 0,0058 0,1406 10,1232 0,0895
72 0 0 0 0,1959 14,1048 0,1247
86 0 0 0 0,1558 13,3988 0,1185
100 0 0 0 0,1553 15,5300 0,1373
114 0 0 0 0,0555 6,327 0,0559
128 0 0 0 0,0604 7,7312 0,0683
142 0 0 0 0,0607 8,6194 0,0762
156 0 0 0 0,0557 8,6892 0,0768
170 0 0 0 0,0995 16,915 0,1496
1,0001 49.419 0,9998 0,9999 113,0386 1
55
V.7.2-Les enthalpies des diffrentes phases liquides

hlt=(1/1515)0,5 (0,403.t+0,000405t2)
Avec : 1515 = 420moy+ 5
Et : = 0.001848 0.00132. 420moy
420moy = i . di420
i= (X i. Mi di415 ) Xi. Mi di415
Tableau V.14 : Calcul de la masse volumique moyenne pour le distillat.
Constituants: Y'i*Mi di415 Y'i*Mi / di415 i.D di420 i.D * di420
C2H6 0,213 0,54 0,3944 0,0046 0,36 0,0016
C3 H 8 26,4836 0,582 45,5044 0,5310 0,5 0,2655

iC4H10 8,4158 0,557 15,1091 0,1763 0,56 0,0987
nC4H10 14,0186 0,579 24,2117 0,2825 0,58 0,1638
iC5H12 0,2880 0,62 0,4645 0,0054 0,619 0,0033
TOTAL 49.419 85 ,6841 1,000 0.5329
Tableau V.15 -Calcul de la masse volumique moyenne pour le rsidu .
15 15 20 20
Constituants: x'i *Mi di4 x'i*Mi / di4 i.R di4 i.R * di4
nC4 H10 11,6 0,579 20,0345 0,1193 0,58 0,0692
iC5 H12 10,1232 0,62 16,3277 0,0972 0,619 0,0602
nC5 H12 14,1048 0,626 22,5316 0,1342 0,62 0,0832
C6 H14 13,3988 0,659 20,3320 0,1211 0,66 0,0799
C7 H16 15,5300 0,684 22,7046 0,1352 0,683 0,0923
C8 H18 6,327 0,703 9 0,0536 0,695 0,0372
C9 H20 7,7312 0,718 10,7676 0,0641 0,718 0,0460
C10 H22 8,6194 0,73 11,8073 0,0703 0,73 0,0513
C11 H24 8,6892 0,74 11,7421 0,0699 0,74 0,0517
C12 H26 16,915 0,748 22,6136 0,1347 0,748 0,1007
TOTAL 113,0386 167,8614 1,000 0,58 0,6721
56
Tableau V.16 :Calcul de la masse volumique pour le dernier plateau.
Constituants
X'li *Mi di415 X'l i*Mi / di415 i.l di420 i.l * di420
nC4 H10 1,7497 0,579 3,0219 0,0210 0,58 0,0122
iC5 H12 12,3266 0,62 19,8817 0,1384 0,619 0,0857
nC5 H12 16,5868 0,626 26,4965 0,1845 0,62 0,1144
C6 H14 13,7793 0,659 20,9094 0,1456 0,66 0,0961
C7 H16 14,2495 0,684 20,8326 0,1450 0,683 0,0990
C8 H18 5,3684 0,703 7,6364 0,0531 0,695 0,0369
C9 H20 6,3070 0,718 8,7842 0,0611 0,718 0,0439
C10 H22 6,8718 0,73 9,4134 0,0655 0,73 0,0478
C11 H24 6,8083 0,74 9,2004 0,0640 0,74 0,0474
C12 H26 13,0186 0,748 17,4045 0,1212 0,748 0,0906
TOTAL 97,0663 143,5815 1,000 0,6743
Tableau V.17 :Calcul de la masse volumique moyenne pour lalimentation.
15 15 20 20
Constituants X'L*Mi di4 X'L*Mi /di4 i.L di4 i.L * di4
C2H6 0,099 0,54 0,0061 0,00005 0,36 0,00001
C3 H8 10,30 0,582 21,1003 0,1758 0,5 0,0879
iC4H10 4,4718 0,557 7,0078 0,0584 0,56 0,0327
nC4H10 8,8624 0,579 12,3012 0,1025 0,58 0,0594

iC5H12 4,4856 0,62 8,9767 0,0748 0,619 0,0463
nC5H12 5,7672 0,626 12,0881 0,1007 0,62 0,0624
C6H14 6,3296 0,659 10,9098 0,0909 0,66 0,0600
C7H16 8 ,33 0,684 12,1783 0,1015 0,683 0,0693
C8H18 5,1072 0,703 4,8324 0,0402 0,695 0,0279
C9H20 4,1472 0,718 5,7760 0,0481 0,718 0,0345
C10H22 4,6292 0,73 6,3414 0,0528 0,73 0,0385
C11H24 4,6644 0,74 6,3032 0,0525 0,74 0,0388
C12H26 15,878 0,748 12,1363 0,1011 0,748 0,0756
TOTAL 77,9736 126,3593 1,000 0,6340
57
Tableau V.18 : Les masses volumiques.
Rsidu Distillat Alimentation Dernier plateau

20
4 moy 0,6721 0.5329 0,6340 0,6743
0,0009 0,0011 0,001 0,0009
15
15 0,6766 0,5384 0,6390 0,6788
Tableau V.19: Les enthalpies.

Enthalpie kcal/kg Enthalpie kcal/kmol
hR 110,86223 12531,71127
hD 31,267798 1545,223309
h0 103,7411 8089,067035
hl 110,0069 10677,97100
V.7.3 -Calcul de lenthalpie des vapeurs quittant le sommet de la colonne
D'aprs le ANNEXE 3 et la temprature T= 143.6 F et P= 191.0473 psia.

HD = 150 BTU / LB = 83,46 kcal / kg= 4124,5097 kcal / kmole
V.7.4 -Calcul de la charge thermique du condenseur
Elle est donne par la formule :

QC = ( D + gx ) (HD - hD )
Qc = (451,8973 +1736,099)*( 4124,5097 -1545,2233)
>Qc = 5643469,0998 kcal/h.
V.7.5 -Calcul de lenthalpie de la charge
Calcul de lenthalpie de la charge

hL = h0 (1- e0)+ e0 H0
D'aprs le ANNEXE 2 la temprature T=347 F et P=200,0528 psia.
H0= 165 BTU/LB = 91,66 kcal/ kg >H0=4659,2469 kcal/kmol
Donc : hL =(0,7345*4659,2469)+((1-0,7345)* 8089,067035)
>hL= 5569,8641 kcal/kmol.
58

Calcul de la charge thermique du four
Elle est calcule selon lquation :
QR =D HD+R hR+ QC - L hL
QR =(451,8973*4124,5097)+(522,8026*12531,71127)+ 5643469,0998 - (974,7*5569,8641)
>QR=8629988,59 kcal/h.
V.8 -Dimensionnement de la colonne
V.8.1 -Calcul du diamtre de la colonne

Le calcul du diamtre se fait par la formule suivante :
4
D=
.
D'o :
Vi: dbit volumique des vapeurs qui quittent le sommet de la colonne en m3 /h.
W : vitesse admissible des vapeurs;

= 8,47 . 105 .

O.
C = 530 : coefficient qui dpend du type de plateau, de la distance entre les plateaux et de
la charge sur les plateaux.
v : Masse volumique des vapeurs.
L: Masse volumique des liquides.
Calcul de dbit volumique des vapeurs qui quittent le sommet
+
Vi = 22.4 . . . .Z

DV : dbit molaire des vapeurs.
T : temprature de tte.
P : pression de tte.
Z : facteur de compressibilit.
DV = g x + D = D + D (rf)op
DV =3,8418*451,8973 +451,8973 > DV = 2187,9963 kmol/h.
59
V.8.2 -Dtermination du facteur de compressibilit
Tr =Tgaz /Tcm Avec Tcm = yi . Tci

Pr =Pgaz/Pcm Avec Pcm = yi . Pci
Tableau V.20 : Calcul des paramtres critiques.
Constituant PC (atm) TC (C) yi,D PC .y 'i,D TC .yi,D

CH4 45,4 -82.4 0 0 0
C 2H 6 48.2 32.4 0,0071 0,3422 0,2300
C 3H 8 41.9 96.4 0,6019 25,2196 58,0231
iC4H10 36 135 0,1451 5,2236 19,5885
nC4H10 37,5 152.2 0,2417 9,0637 36,7867
iC5H12 33,4 187.4 0,004 0,1336 0,7496

nC5H12 33, 3 196,6 0 0 0
TOTAL - - 1 39,9827 115,3779
Do :
62 13
Tr = 115.3779 = 0.53 C Pr = 39.9 = 0.32 atm
D'aprs lannexe 4, nous trouvons : Z = 0,8

Donc le dbit volumique:
+
Vi = 22.4 . . . . Z

Vi = 22.4 . 2187,9963. +

. . 0.8

Vi =1,728 m3/s > Vi = 2759,4577 m3/h
V.8.3 -Calcul de la vitesse admissible des vapeurs
Calcul de v masse volumique des vapeurs
..
v =
..(+)
: pression au sommet
M : Masse volumique des vapeurs
60
....
v =
..(+)
>v =23,3726 kg/m3
Calcul de la masse volumique du liquide L 54C

On calcule d'abord la masse volumique du liquide (reflux):
1
420 = , >420 = 532.9 kg/m3

=1 20
4,
1 = 420 ( 20)
=1.848-(0.00132*532.9) > =1.1445
1=532.9-1.1445*(62-20) > 1=484.831 kg/m3
L= 454= 420-(54-20)
Nous aurons L=501.2kg/m3
484.83123.3726
= 8,47 . 105 . w=8,47 .10-5.530
23.3726
> W=0.1994 m/s
Do la vitesse admissible des vapeurs est de : W = 0,1994 m /s
> D = 3,35 m.
le diamtre est estim : D =3,35 m
61
V.8.4 -Calcul de la hauteur de la colonne

La hauteur de la colonne est la somme des distances suivantes :
hT = h1 +h2 + h3 +h4 + h5 + h6
h1 : la hauteur du fond suprieur
h2 : la hauteur de la zone de rectification.
h3 : la hauteur de la zone d'alimentation.
h4 : la hauteur de la zone d'puisement.
h5 : la hauteur du fond infrieur
h6: la hauteur de la jupe
Donnes de dpart
h1 = 1.695 m
h5 =12.806 m
h6 = 4.62 m
h2 = hp( NR 1 )
hp : distance entre deux plateaux successives = 0,6 m.
Nous avons :
>h2 =0,6*(16-1) = 9,4 m.
>h3 = 2 hp =2*0,61 = 1,2 m.
h4= hp (NE- 1)
NE = Ne / 0,6
>h4 =0,6*(15-1)= 8,4 m.
H = 38,12 m.
V.9 -Calcul Du Ballon De Reflux V209

Pour le cas dune condensation totale du produit de tte de la colonne, nous utilisons la
mthode ci-aprs.
Le volume du distillat (GPL) sur lequel sera bas le calcul du diamtre et la longueur du
ballon
Donnes de dpart
Dbit de la charge D = 451.8973 Kmol/h.
Nature de la charge : liquide.
Temprature dans le ballon : 54 C.
Pression dans le ballon : 11,9 kg/cm2
62
Densit de la charge : d420=0,5367
Pour des valeurs de pression 7 kg/cm2
V.9.1 -Dbit massique du distillat
Dm = D. MD
MD: masse molculaire moyenne du distillat.
Dm : le dbit massique du distillat donn par lquation
On a: M= 49.419 g > Dm = 22332.3126 kg/h.
Dbit volumique du distillat
v= Dm/ 454. Avec : 4t = 420 (t-20)
> 420 = 0.5329
> =0.001848 - 0.00132 420 > 4

54 = 0,493984552 g/cm3
On trouve : V = 45,21626362 m3/h.

La retenue correspondante 5 min de soutirage est de :
V1 = (V.5)/60 >V1= 3,7680 m3.
La retenue correspondante 3min de reflux de tte gD sera gale :
V2 = (3 rf V) /60 >V2 = 8,6855 m3.
Donc la retenue totale est gale : VT = 12,4535 m3.
Dimensions du ballon
D b = 0.75 (V)1/3. >V = VT .2 =24,9070 m3. >Db = 2,20 m.
Donc : >R=1,1 m.
Prenons h = 20 cm pour les hauteurs des segments circulaires suprieurs et infrieurs, le rapport :
>h/R = 0,180.
donne pour la valeur trouve de h/R :
S= 0,1410188 R2. S :Laire de deux segmen >S= 0,17 m2
Laire entre le NB et celui NH est gale
F = R2- S >F= 3,62 m2
On dduit la longueur du ballon L : L= V/F >L= 6,90 m
Pour le cas dune condensation totale du produit de tte de la colonne. Pour une pression du ballon
suprieure 7 bars.
3 < L/D < 5
On vrifie le rapport L/D >L/D = 6,90/2,20 = 3,13 m.
63
V.10 -Calcul du vrification du four
V.10.1. -Calcul de la quantit de chaleur apporte par le gaz combustible
QGC=QR/
= 0.78 Avec (donne pratique du constructeur).
QGC=8629988,59/0,78 >QGC =11064087,9 Kcal/h.
V.10.1.1 Dbit du gaz combustible :
G=QGC / PCI
Avec :
G : dbit de gaz combustible (Nm3/h,Kg/h).
PCI : pouvoir calorifique infrieur (Kcal/Nm3,Kcal/Kg).
Calcul de la combustion
Le calcul est bas sur lquation stoechiomtrique de la combustion des hydrocarbures ayant la
formule chimique gnrale CnHm. CnHm + (n+m/4) O2 n CO2 + m /2 H2O.
On dtermine daprs cette quation la quantit thorique dO2 ncessaire la combustion et la
quantit des fumes obtenue (CO2 ,H2O).
Tableau V.21 :Nombre de mole dO2 ainsi que les produits de la combustion (CO2, H2O).
bre bre bre
Constituants Gaz combustible N de mole N de mole N de mole
(%mole)
De H2O De CO2 De O2
N2 0,0561
CO2 0,0021 0,0034
C1 0,8353 1,6706 0, 8353 1,6706
C2 0,0770 0,2310 0,1540 0,2695
C3 0,0194 0,0776 0,0582 0,0970
iC4 0,0027 0,0108 0,0108 0,0175
nC4 0,0043 0,0215 0,0172 1,8167
iC5 0,0010 0,0050 0,0050 0,0080
nC5 0,0011 0,0066 0,0055 0,0088
C6 0,0010 0,0070 0,0060 0,0950
TOTAL 1 2,0301 0,2601 3,9831
V.10.2 -Calcul du pouvoir calorifique suprieur PCS
PCS = PCSi. YIgc
64
Tableau V.22 : Calcul de PCS.
Constituants Gaz combustible (%mol) PCS(Kcal/Nm) PCSi

N2 0,0561 / /
CO2 0,0021 / /
C1 0,8353 8859 7399,9227
C2 0,0770 15523 1195,2710
C3 0,0194 22101 428,7594
iC4 0,0027 28576 77,1552
nC4 0,0043 28644 123,1692
iC5 0,0010 35123 35,123
nC5 0,0011 35190 38,709
C6 0,0010 41754 41,754
TOTAL 1 9389,8635
3
Donc :>PCS=9389,8635 Kcal/Nm .
Calcul du pouvoir calorifique infrieur PCI

PCI = PCS 10,9N
>PCI = 9389,8635 - (10,9*2,0301). >PCI = 9367,7354 kcal / Nm3.
Calcule de la masse volumique du gaz combustible
com= m / 22,4 m= n.M = M
n= 1 ; parce quon a une mole de gaz combustible.
M : masse molaire du gaz combustible.
65
Tableau V.23 : Calcul de la masse volumique de gaz combustible.
Constituants Gaz combustible Mi Mi*YiGC

(%mol)
N2 0,0561 44 2,4684
CO2 0,0021 28 0,0588
C1 0,8353 16 13,3648
C2 0,0770 30 2,3100
C3 0,0194 44 0,8536
iC4 0,0027 58 0,1566
nC4 0,0043 58 0,2494
iC5 0,0010 72 0,0720
nC5 0,0011 72 0,0792
C6 0,0010 86 0,0860
TOTAL 1 19,6986
com= 19,6986/ 22,4 >com= 0,87941071 Kg/Nm3.

>PCI= 9389,8635 /0,87941071 >PCI=9367,7354 Kcal/kg.
A la fin on aura : G=QGC / PCI =11064087,9 / 9367,7354
>G = 1181,08459 Nm3/ h = 1038,65844 Kg/h.
66
Comparaison entre les rsultats du calcul et les paramtres de design
Tableau V.24 : Paramtres opratoires du design et calcul de la colonne

T202.
Dsigne Calcul
Dbit dalimentation [kmole/h] 2950 974,7
Diamtre de la colonne [m] 3,50 3,35
Hauteur de la colonne [m] 39 38,13
Nombre de plateaux 31 32
Dbit du reflux froid [kmol/h] 2800 1736,099
Dbit du reflux chaud [kmol/h] 3985 758,06377
Temprature du sommet [C] 68 62
Temprature du fond [C] 190 167
Temprature dalimentation [C] 190 165
Temprature du reflux froid [C] 60 54
Temprature du reflux chaud 220 189
[C]
Diamtre du ballon de reflux [m] 2,50 2.20
Longueur du ballon de reflux 7,8 6 ,90
[m]
Interprtation
Daprs les rsultats cits ci-dessus et le calcul de dimensionnement on peut tirer la conclusion
suivante :
On a trouv des valeurs bien inferieures a celle des valeurs dans le cas design donc on propose
daugmenter le dbit de la charge alimentant la colonne pour augmenter la production et les
bnfices.
67
V.11 Partie simulation

V.11. 1 -La simulation dun procd
La simulation peut tre dfinie comme lutilisation dun modle mathmatique pour
gnrer une description du comportement dun systme physique (procd).
Lavantage majeur de la simulation est la fourniture dun bon aperu du comportement du
systme rel. Cet aperu peut tre difficile avoir travers lexprience et lintuition seule, en
particulier pour les systmes complexes avec plusieurs variables interAactives.
Ds que le modle mathmatique rpond aux ajustements de changement des paramtres
comme un procd rel, la simulation peut tre classe comme une mthode commode qui ne
cote pas chre et sans danger pour la comprhension du comportement du procd rel et sans
intervention sur le procd en cours de fonctionnement.
Les simulateurs disponibles pour concevoir de nouvelles units et pour optimiser les
procds industriels sont : ASPEN PLUS, ChemCAD-III, HYSIM, PRO-II et HYSYS.
Les diffrentes taches quun simulateur de procds devrait effectuer sont :
la rsolution des bilans matire et dnergie.
loptimisation du procd.
le dimensionnement des quipements.
lvaluation conomique du procd. [16]
V.11 .2 Choix du modle thermodynamique

Les logiciels de simulation donnent accs plusieurs modles thermodynamiques pour la
prdiction de lquilibre liquide-vapeur, denthalpie et dentropie ainsi que des proprits de
transport.
La russite de la simulation dpend du choix du modle thermodynamique, parce que ce dernier
est tabli pour une classe de fluides et un domaine de conditions de pression et de temprature
recommands en utilisant des hypothses et des suppositions pratiques.
Lquation de Peng - Robinson est recommande pour le systme tudi. [17]
V.11.3. Vrification des paramtres design du dbutaniseur
Dans cette premire tape, notre objectif est de vrifier la validit du modle
thermodynamique choisi (Peng-Robinson), pour son utilisation ultrieure.
68
V.11.3.1. Composition de charge du dbutaniseur

Donnes du design
1- Le dbitdalimentation:2960 Kmol/h
2- La temprature dentre:190 C
3- La pression : 16.2 Kg/ cm2
Tableau V.25 : Composition design de lalimentation du dbutaniseur.
Constituants Xi,
C2 H6 0,0058
C3 H8 0,1190
iC4 H10 0,0450
nC4 H10 0,0941

iC5 H12 0,0535
nC5 H12 0,0770
C6 H14 0,1109
C7 H16 0,1117
C8 H18 0,0916
C9 H20 0,0758
C10 H22 0,0554
C11 H24 0,0431
C12 H26 0,1170
TOTAL 1,0000
V.11.3.2. Procdure de la simulation

La procdure de calcul du dbutaniseur dans lenvironnement HYSYS ncessite les donnes
suivantes :
La composition, dbit, temprature et pression de la charge.
Le design de la colonne tels que : diamtre, type et nombre total de plateaux.

La position du plateau dalimentation
69
Spcification de l efficacit des plateaux

Le simulateur ncessite la spcification de trois dautres paramtres opratoires qui
reprsentent le degr de libert de la colonne (tableau V.5):
Temprature au condenseur (TCondenseur);
Temprature au rebouilleur (TRebouilleur);
Le dbit de rebouilleur (Rb).
Tableau V.26 : la comparaison entre cas design rel et cas design simul.
Paramtre Cas design rel Cas design simul

(HYSIS)
Temprature dalimentation (C) 190.5 185

2
Pression dalimentation (kg /cm g) 16.9 16.9
Temprature du sommet (C) 72 74.1
2
Pression du sommet (kg /cm g) 14.2 14.2
Temprature de condenseur (C) 57.6 58.7
2
Pression de condenseur (kg /cm g) 14.2 13.8
Temprature de fond (C) 205 197.5
2
Pression de fond (kg /cm g) 15.4 15.3
Temprature de rebouilleur (C) 220.9 211
OPTIMISATION
A partir de la vrification de design par simulateur, on dit que ce simulateur est en bon
fonctionnement.
Donnes de dpart:
Le dbit dalimentation : 974,7 Kmol/h
Temprature dalimentation : 165C
Pression dalimentation : 13.25 bar
70
Figure V.1 : Simulation du procd de la section

Le tableau suivant reprsente les dfrentes fractions de la charge et distillat et rsidu
Figure V.2: les fractions de la charge et du produit.
71
V.11.4 -Profil de temprature :
Figure V.2 :Profil de temprature.

Interprtation du graphe
Le profil de temprature dans la colonne de stabilisation prsente une volution normale le long
de la colonne, sauf une prsence dun pic qui est remarqu au niveau du plateau dalimentation
(au niveau du 27eme), ce pic explique par la diffrence entre la temprature dalimentation et
celle du plateau.
V.11.5 Profile de pression
Figure V.3 : Profil de pression.
72
Le profil de pression actuel prsente une volution normale le long de la colonne

(augmentation de pression du 1er plateau en tte, jusquau dernier plateau du fond).
V.11.6-Profil du dbit
Figure V.4 : Profil du dbit .

Daprs la figure on constate que le dbit de vapeur dans la section de rectification est
me
suprieur celui de liquide mais quand on arrive au 27 plateau le dbit de liquide augmente
dune faon considrable par rapport au dbit de vapeur qui est diminu cause de la fraction
liquide vapeur de lalimentation (v/f=0.2208).
V.11.6.1-Profil de la composition de la charge
Figure V.5 : Profil de la composition de la charge.
73
La fraction du GPL le long de la colonne a une augmentation progressive pour atteindre les
o elle se stabilise, par contre la fraction du nC+5 diminue progressivement au top de la colonne.
Les diffrentes spcifications des produits dans le tableau suivant

Tableau V .27 : Des diffrentes spcifications des produits.
Paramtres Produit
Dbit distillat (Kgmol/h) 449,1
Dbit Reflux (Kgmol/h) 1706
Dbit Rsidu (Kgmol/h) 525,6
Dbit alimentation (Kgmol/h) 974,7
Temprature de tte (C) 63,5
Temprature de fond (C) 165
Temprature de reflux (C) 48,4
6
Quantit de chaleur dans le rebouilleur (Kcal/h) 8,478*10
6
Quantit de chaleur dans le condenseur (Kcal/h) 5,222*10
-
Teneur en C2 : 3 % mol (max) 0,002
+
Teneur en C5 : 0.4 % mol (max) 0,00045
Tension de vapeur Reid (TVR) 10 psia 6,05
Taux de reflux 3,8
74
Conclusion Gnrale
Conclusion Gnrale
Daprs la formation effectu au niveau des units de traitement de gaz (module 1) de
Hassi Rmel, nous a permis dacqurir et dapprofondir nos connaissances thoriques
leffet de les mettre en pratique.
Par ailleurs, cette formation nous avons galement permis de visiter les diffrents
quipements industriels tel que : les colonnes, les arocondenseurs, les changeurs, les
pompes.etc.
Nous avions pour objectif de vrifier les paramtres opratoires et les dimensions des
diffrents quipements de la section de stabilisation du condenst et les calculs effectus
ont montr que les paramtres de fonctionnement optimiss analytiquement et
numriquement sont infrieurs aux celles du design.
Aprs avoir effectu une comparaison entre les valeurs calcules, et celles du design des
paramtres opratoires ainsi que les dimensions des diffrents quipements de la section de
stabilisation du condenst. Nous avons montr que le manque de la richesse en GPL et
condensat dans le gaz brut traiter influe beaucoup sur les paramtres de fonctionnement
de dbutaniseur. Ce dernier peut fonctionner avec une plus grande capacit de production.
Donc on a propos daugmenter le dbit de lalimentation et ainsi travailler avec les
paramtres calculs.
Finalement, la meilleur solution que nous pouvons proposer cest de faire une rpartition
des charges liquides dun train, la charge introduite la section fractionnement et
stabilisation vers les deux autres trains.
75
Bibliographie
Bibliographie
[1] Le gaz Natural production traitements transports ; prfrence de Lucien SAJUS directeur
dlgu gnral technique prsident de la communication gaz du comit dtude ptrolires et
Marines (CEPM). Editions technique 27 rue ginoux 75737 Paris cedex 15.
[2] fiches danalyses. Et projet module IV
[3] Rapport de stage du module I, Centre de formation Hassi Rmel.2015.
[4]Procds de ptrochimie: caractristiques techniques et conomiques A Chauvel, G

Lefebvre, L Castex - 1985 - Editions Technip
[5]Rapport de stage du center stockage et transfert facilit, Centre de formation Hassi

Rmel.2010.
[6] Manuel opratoire module 1, Hassi Rmel. Algrie, Stone and Webster corporation
.1976.
[7]site;www.archive.wikiwix.com/cache/?url=https%3A%2F%2Ffanyv88.com%3A443%2Fhttp%2Fwww.horizongpl.com%2F
tout_gpl_cquoi.php
[8]site ;www.memoireonline.com/12/13/8151/m_La-dynamique-des-prix-des-GPL-au-regard-
des-determinants-du-marche-spot-americain12.html
[9]site ;www.memoireonline.com/03/07/392/m_performance-commerciale-reseau-distribution-
sctm-gaz10.html
[10]Stockage des produits ptroliers .destin pour le groupe oil and gas procssing .cours raliss
par ; BENMOUNAH .A , IAP.cu boumerdes .2006
[11]NaturaGasCondensate [http://archive.wikiwix.com/cache/?url=https%3A%2F%2Ffanyv88.com%3A443%2Fhttp%2Fwww.mar
athonpetroleum.com%2Fmsds%2F0245MAR001.pdf ] ; Marathon Oil Company MSDS
[12] Jean-Martin Brault. Conception et synthse de procds. Ecole polytechnique

Montral.2005.
[13] Brochures :R.ABOULLAEV ,V.KOSSIAKOV. Thorie et calcule de la rectification des

mlange complexes
[14]: cours opration unitaires - la rectification FHC (IGCP/08) boumerdes 2011/2012.

Bibliographie
[15] Emilian Koller, Aide-mmoire Gnie chimique, 3e dition .
[16] BELHOCINE. Mokhtaria ; Etude et optimisation des paramtres de fonctionnement de la

section fractionnement. Universit dOran (U.S.T.O). 2010.
[17] Formation professionnalisant PRO/EXP4 groupe 2 .IFP ingnierie de gaz traitement du

gaz naturel ; HMD bas 24 -02-07: Mr ALEXANDER Scondo. mai 2015
ANNEXE 1
Diagramme des coefficients dquilibre des hydrocarbures (Scheibel et Jenny)

ANNEXE 2
Diagramme des enthalpies total des paraffines hydrocarbures liquides en
fonction de la temprature et la masse molaire .
ANNEXE 3
Diagramme des enthalpies total des paraffines hydrocarbures vapeurs en
fonction de la temprature et la masse molaire .
ANNEXE 4
Coefficient de compressibilit en fonction de la pression rduite

Kerbouche Ali Kaabouche Boumadiane

Transféré par

Droits d'auteur :

Formats disponibles

Kerbouche Ali Kaabouche Boumadiane

Transféré par

Informations du document

Description originale:

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Kerbouche Ali Kaabouche Boumadiane

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

UNIVERSITE KASDI MERBAH OUARGLA

Vrification et optimisation des paramtres de

Soutenu Publiquement le : 31/05/2016

Ingnierie de Gaze Natural

A nos trs chers parents, sources inpuisables damour et de tendresse, en reconnaissance

ici le symbole de lamour et lintime attachement que vous nous portez.

sincres tmoignages de notre profonde affectation et de notre haute considration .

Ingnierie de Gaze Natural

Liste des tableaux

: Hauteur de la zone de rectification.

Chapitre I : Gnralit Sur Le Gaz Naturel

I.1-Le gaz naturel dans le monde 2

II.1 -Prsentation de la rgion HASSI RMEL .. 10

II.5.2 .1-Description du Procd 20

Chapitre III : Gnralits sur les produits

Chapitre IV : Thorie de la rectification

V.3.6 -Bilan thermique de la colonne.... 40

Chapitre V : Parties de Calculs

V.1 -Bilan de matire de la colonne T202 43

V.10.1.1 Dbit du gaz combustible. 64

Aujourdhui, lindustrie gazire aborde le nouveau millnaire dans un contexte de

Ingnierie de Gaze Natural

I.1-Le gaz naturel dans le monde

I.2 - Gaz naturel en Algrie

Les rserves algriennes prouves en gaz sont de lordre de 4,550

Champs exploits par SONATRACH

I .3 -La production de gaz Naturel

Figure I.1 :Production mondiale du GN [3]

I.4 -La consommation de gaz naturel

Figure I.2 : Consommation mondiale du GN.

I.5 - Traitement du gaz naturel

Procd PRITCHARD : il consiste refroidir le gaz par change thermique et une

1. Pouvoir calorifique suprieur PCS

Il correspond la chaleur dgag lorsque tous les produits de combustion sont

2. Pouvoir calorifique infrieur PCI

Indique la nature des composs dhydrocarbures et dautres constituants du gaz, et

f) Composition du gaz naturel

Le mode inorganique de formations des gaz dhydrocarbures reste trs secondaire.

I.7 - Types du gaz naturel

Lapparition dune phase liquide dpend des conditions de temprature et de pression

I.9 - Avantages et inconvnients du gaz naturel

Ingnierie de Gaze Natural

II.1 -Prsentation de la rgion Hassi Rmel

II.2 -Situation gographique de Hassi Rmel

II.3 -Les principales tapes du dveloppement de Hassi Rmel

Le dveloppement du gisement de HASSI RMEL a t ralis en plusieurs tapes, rpondant

II.4-Organisation du champ de Hassi Rmel

Le champ de Hassi Rmel comprend trois secteurs :

Tableau II.1: Capacit de production des diffrents modules.

Produits Modules Modules Modules Module Djebel Bissa

Condenst (m) 224000

GPL (m) 68000

II.4.1 -Les salles de contrle

II.4.2 -Le centre de traitement dhuile (CTH)

II.4.3 -Station de rcupration du gaz associ (SRGA)

II.4.4 -CTG-HR- SUD

II.4.6 -Centre national de dispatching gaz (CNDG)

II.4.7 -Les Stations Boosting