Calcul de Vérification Et Optimisation Des.......

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE
MINISTERE DE L'ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITE MOHAMED KHIDER BISKRA
FACULTE DES SCIENCES ET DE LA TECHNOLOGIE

DPARTEMENT DE CHIMIE INDUSTRIELLE
Mmoire de Fin d'tude en Vue de l'Obtention du Diplme de Master

En Gnie des Procds
Option : Gnie chimique
THME
Calcul de vrification et optimisation des

paramtres de fonctionnement du colonne
de Stabilisation PX0-CB-21-01
de l UTBS
Prsent par:
RAHIM Khaled et BEN RAHLA Sohaib
Prsident: Mme ADJEL Fatima

Examinateur: Mme MENASRA Hayet
Encadreur: Mr BEN TEMAM Hachemi
Promotion Juin 2012

Mes remerciements les plus respectueux vont mon promoteur Mr BEN TEMAM.H
qui nous a dirig et qui a t la hauteur de sa noble tche.
Nous tenons aussi remercier profondment tout le

Personnel de Hassi-Messaoud :
Mr. GHORAB , chef de Division Exploitation ;
Messieurs : LASMAR FOUAD, MEDFOUNI SAMIR et BEN AMARA HOSIN et
tous les travailleurs de la Division Exploitation(UTBS) ;
Et avec immense plaisir que nous prsentons nos vives gratitudes tous
ceux qui ont contribu de prs ou de loin au droulement de notre stage pratique.
Nous remercier galement ;
les membres du jury pour lhonneur quil nous fait en acceptant lvaluation
de notre travail.
Finalement nous remercions toutes les enseignants de DPARTEMENT DE
CHIMIE INDUSTRIELLE
nombreuses sont les personnes que nous n`avons pas cit.
RAHIM KHALED
BEN RAHLA SOHAIB.

Ww|vtvx
Je ddie ce modeste travail :
Ma trs chre mre,
Mon trs cher pre,
Mes trs chers frres et soeurs,
Toute ma famille,
A tous les enseignants et toutes les enseignantes de

chimie industrielle,
Et tous mes amis.
RAHIM K
Ww|vtvx
Jai le grand honneur de ddier ce modeste travail :

A ceux qui mont encourag et soutenue moralement et
matriellement pendant les moments les plus difficiles et durant
toute ma vie.
Ma trs chre mre Mon trs cher pre
A mes trs chrs frres et mes surs
Sans oublie rendre un grand hommage
Mes cousins dans nos mmoires et toute la famille

Khaled .R
BEN RAHLA.S
SOMMAIRE
Nomenclature
INTRODUCTION GENERALE..01
CHAPITRE I : Historique de champ de Hassi Messaoud

I.1. Historique...........02
I.2. Situation du champ de Hassi Messaoud.02
I.2.1. Situation gographique.02
I.2.2. Situation gologique.03
I.3. Zones et numrotation des puits.04
I.4. Intrt ptrolier...05
CHAPITRE II : Gnralits sur le ptrole brut

II.1. Gnralits sur le ptrole brut..07
II.2. Dfinition du ptrole brut.07
II.3. Origine du ptrole brut..07
II.4. Classification des ptroles brut.07
II.4.1. Classification industrielle08
II.4.2. Classification chimique...08
II.4.3. Classification conomique (technologique)....08
II.4.3.a. Daprs la teneur en soufre..........08
II.4.3.b. Selon la teneur potentielle en fractions claires PF=350c... .08
II.5. Composition des ptroles brut..09
CHAPITRE III : Description de lunit

III.1. Introduction.10
III.2. Vue gnrale de LUTBS10
III.3. Description de lunit..13
III.3.1. Units de Production..13

III.3.2. Rseau de Collecte Production de Satellite....13
III.3.3. Traitement et Stabilisation dHuile14
3-3.a) Sparation dHuile ..14
3-3.b) Dessalage.....15
3-3.c) Stabilisation......16
III.3.4. Stockage, Comptage et Expdition de lHuile Stabilise..19
III.3.5. Traitement et du stockage de lhuile hors spec..19
III.3.6. Unit de Compression de Gaz....20
III.3.7. Lunit de traitement des eaux huileuses..21
III.3.8. Les systmes utilits...21
CHAPITRE IV : Thorie de rectification

IV.1. Introduction.23
IV.2. Dfinition....23
IV.3. Description dune colonne de rectification.....24
IV.4. Principe de fonctionnement dune colonne de rectification....25
IV.5. Types des colonnes de rectification....26
IV.5.1. Colonnes plateaux..26
IV.5.2. Colonnes garnissage....27
IV.6. Lois fondamentales de la rectification...27
IV.6.1. Lois de DALTON..27
IV.6.2. Lois de RAOULT et dHENRY27
IV.6.3. Bilan matire de la colonne et de ses zones...28
IV.6.4. Utilisation de la volatilit relative..32
IV.6.5. Choix de la pression dans la colonne....33
IV.6.6. Rgime de la temprature dans la colonne34
IV.6.7. Taux de vaporisation et composition des phases liquides et vapeurs de la

charge...34
IV.6.8. Choix du taux de reflux.35
IV.6.8.a. Fonctionnement reflux minimal...36
IV.6.8.b. Taux de reflux optimal36
IV.6.9. Bilan thermique de la colonne....36
IV.6.10. Nombre de plateau thorique reflux opratoire.37
IV.6.11. Dtermination des dimensions de la colonne..39
IV.6.11.a. Hauteur de la colonne...39
IV.6.11.b. Nombre de plateaux.....39
IV.6.11.c. Espacement entre les plateaux...40
IV.6.11.d. Diamtre de la colonne..40
CHAPITRE V : Partie calcul

Prsentation du projet.......41
V.1. Donnes de dpart.....42
V.2. Les dbits de distillat et du rsidu de la colonne PX0-CB-21-01...44
V.3. Pression dans les diffrents zones de la colonne.....48
V.3.1. Pression de la charge froide (reflux)...48
V.3.2. la pression au sommet de la colonne..48
V.3.3. la pression au fond de la colonne...49
V.3.4. la pression de zone dalimentation chaude.49
V.4. Dtermination de la temprature.49
V.4.1. Pour lalimentation chaude et froide......49
V.4.2. Dans le sommet de la colonne50
V.4.3. Dans le fond de la colonne.51
V.5. Volatilits relatives..52
V.6. Taux de vaporisation e0 et composition des phases vapeur et liquide

dalimentation...54
V.7. Taux de rebouillage.....57
V.8. Bilan thermique de la colonne.....58
V.8.1. Dtermination des enthalpies des phases vapeurs et liquides.59
V.8.2. Calcul des charges thermiques...62
V.9. Calcul du diamtre de la colonne62
V.10. Calcul de la hauteur de la colonne..66
CHAPITRE IV : Partie optimisation et simulation

VI.1. Dfinition...68
VI.2. Le choix du modle thermodynamique..68
VI.3. Remarques concernant lutilisation des logiciels de simulation69
VI.4. Le choix des paramtres..71
VI.5. Prsentation des rsultats et interprtation..72
VI.5.1. Influence de taux reflux froid sur le fonctionnement de la colonne de

stabilisation...72
VI.5.2. Prsentation et interprtation des profils de la colonne de stabilisation..........76
VI.8. Tableaux comparatifs.79
CONCLUSION GNRALE....80
RFRENCES BIBLIOGRAPHIQUES
ANNEXES
Nomenclature
Symboles Dsignation Unit
C coefficient qui dpend du type de plateau
Di Dbit molaire du constituant i dans le distillat [Kmol/h]
Dim Dbit massique partiel du constituant i dans le distillat [Kg/h]
D Diamtre de la colonne [m]
d 420 Densit mesure 20C
d 415 Densit mesure 15C
e'0 Taux de vaporisation de la charge

F Dbit molaire de la charge froide [Kmol/h]
HT La hauteur totale [m]
H1 Hauteur de lespace libre au dessus du plateau suprieur [m]
H2 Hauteur de LA zone de rectification [m]
H3 Hauteur de la zone dexpansion [m]
H4 Hauteur de la zone dpuisement [m]
H5 Hauteur de lespace libre au-dessous du plateau infrieur [m]
H6 Lespace entre le plateau 13 et le plateau 14 [m]
Hp Distance entre deux plateaux [m]
hL Enthalpie molaire de la charge chaude [Kcal/kmol]
hF Enthalpie molaire de la charge froide [Kcal/kmol]
hD Enthalpie molaire de distillat [Kcal/kmol]
hR Enthalpie molaire de rsidu [Kcal/kmol]
H0 Enthalpie molaire ltat vapeur [Kcal/kmol]
h0 Enthalpie molaire ltat liquide [Kcal/kmol]
kiL Coefficient dquilibre du constituant i dans la charge chaude
kiF Coefficient dquilibre du constituant i dans la charge froide
kiD Coefficient dquilibre du constituant i dans le distillat
kiR Coefficient dquilibre du constituant i dans le rsidu
Kuop Facteur de caractrisation
L Dbit molaire de la charge chaude [Kmol/h]
Mi Masse molaire du constituant i [Kg/kmol]
Mn Masse molaire moyenne de la charge chaude et froide [Kg/kmol]
MmR Masse molaire moyenne du rsidu [Kg/kmol]
MmD Masse molaire moyenne du distillat [Kg/kmol]
NR Nombre de plateaux dans la zone de rectification
NE Nombre de plateaux de la zone dpuisement
np1 Nombre de plateau entre le sommet et lalimentation chaude
P Taux de rebouillage
Pr La Pression de la charge froide [bars]
PS Pression au sommet de la colonne [bars]
Pf Pression au fond de la colonne [bars]
Pac La pression de zone dalimentation chaude [bars]
PCm Pression critique moyenne [bars]
Pr m Pression rduits moyenne [bars]
QF Capacit thermique de lalimentation froide [Kcal/h]
QL Charge thermique de lalimentation chaude [Kcal/h]
QD Capacit thermique du distillat [Kcal/h]
QR Capacit thermique du rsidu [Kcal/h]
QRb Capacit thermique du rebouilleur [Kcal/h]
R Dbit molaire du constituant i dans le rsidu [Kmol/h]
Rim Dbit massique partiel du constituant i dans le rsidu [Kg/h]
R Dbit molaire de rsidu totale [Kmol/h]
Tac La temprature de lalimentation chaude [C]
Taf La temprature de lalimentation froide [C]
Ts Temprature au sommet de la colonne [C]
Tf Temprature au fond de la colonne [C]
TCm Temprature critique moyenne [C]
Tr m Temprature rduits moyenne [C]
Vi Dbit volumique des vapeurs qui quittent le sommet de la colonne [m3/h]
Vm Dbit molaire de la vapeur qui quitte la colonne [Kmol/h]
W La vitesse des vapeurs [m/s]
X'Fi Fraction molaire du constituant i dans la charge froide [mol]
XiF Fraction massique du constituant i dans la charge froide [mol]
X'Li Fraction molaire du constituant i dans la charge chaude [mol]
XiL Fraction massique du constituant i dans la charge chaude [mol]
XRi Fraction molaire du constituant i dans le rsidu [mol]
X0.i Concentration molaire du constituant (i) dans la phase liquide [mol]
Y0.i Concentration molaire du constituant (i) dans la phase vapeur [mol]
YDi Fraction molaire du constituant i dans le distillat [mol]
Z Facteur de compressibilit
Lettres Grecques
Symboles Dsignation Unit
P1 La perte de charge dues aux rsistances des plateaux [bars]
P2 La perte de charge au niveau de chaque plateau [bars]
Lc,i Volatilit relative du constituant i dans la charge chaude
Lf,i Volatilit relative du constituant i dans la charge froide
Li Volatilit relative du constituant i dans la zone dalimentation
Di Volatilit relative du constituant i dans le distillat
Ri Volatilit relative du constituant i dans le rsidu
iDm Volatilit relative moyenne dans la zone de rectification
iRm Volatilit relative moyenne dans la zone dpuisement
V Le poids spcifique de la vapeur dans les conditions de fonctionnement [Kg/m3]
L Le poids spcifique du liquide dans les conditions de fonctionnement [Kg/m3]
Li La masse volumique du constituant (i) aux conditions opratoires [Kg/m3]
Constante qui varie entre les volatilits relatives des constituants cls
lourds et le plus volatiles.
Abbreviations
API : American Petroleum Institute.
CFPA : The Center for Professional Advancement.
UTBS : unit de traitement de brut sud.
CIS : centre industriel sud.
GPL : Gaz de ptrole liqufi
TVR : tension de vapeur Reid
HEH : Haoud El Hamra

Liste des Fuguer
Fig. I. 1: Situation gographique du champ de Hassi Messaoud...3
Fig. I. 2: Rpartition des zones ptrolifres et gazifres4
Fig. I. 3: Zones de production...5
Fig. III. 1: Plan gnral de lUTBS............12
Fig. III. 2: Rseau de Collecte du Brut..14
Fig. III. 3: Sparation dHuile...........15
Fig. III. 4: Dessalage...16
Fig. III. 5: Stabilisation..18
Fig. III. 6: Systme dHuile Hors-Spec..22
Fig. III. 7: Compression de Gaz.21
Fig.IV. 1: Schma simplifi Rectification..23
Fig.IV. 2: Schma d'une colonne plateaux.25
Fig.IV. 3: Reprsentation dune tour plateaux...........26
Fig.IV. 4: Schma descriptif du fonctionnement de la zone de rectification30
Fig.IV. 5: Schma descriptif du fonctionnement de la zone dpuisement...........31
Fig.IV. 6: Schma descriptif du fonctionnement de la zone dalimentation33
Fig. VI. 1: Schma technologique de la simulation...........70
Fig. VI. 2: Schma synoptique du colonne de stabilisation...........72
Fig. VI. 3: Influence du taux de reflux en fonction la TVR..73
Fig. VI. 4: Influence du taux de reflux en fonction du poids molculaire...........74
Fig. VI. 5: Influence du taux de reflux en fonction de duty four..........74
Fig. VI. 6: Influence du taux de reflux en fonction de la temprature de sommet..75
Fig. VI. 7: Influence du taux de reflux en fonction de la temprature de fond de

colonne..........75
Fig. VI. 8: Influence du taux de reflux en fonction de U.A Prchauffeur du
stabiliseur..............................................................................................................76
Fig. VI. 9: Le profil de temprature au niveau de chaque plateau...77
Fig. VI. 10: Le profil de pression au niveau de chaque plateau77
Fig. VI. 11: Le Profil des dbits molaires du liquide et de la vapeur au niveau de chaque
plateau ..78
Liste des tableaux
Tableau [V-1]: Composition de la charge chaude...42
Tableau [V-2]: Composition de la charge froide.43
Tableau [V-3]: Composition et dbits molaire et massique du distilla.47
Tableau [V-4]: Composition et dbits molaire et massique du rsidu 48
Tableau [V-5]: Coefficients d'quilibre des constituants Dalimentation.50
Tableau [V-6]: Coefficients d'quilibre des constituants du distillat51
Tableau [V-7]: Coefficients d'quilibre des constituants du residu..52
Tableau [V-8]: les volatilits relatives de la charge de la colonne..53
Tableau [V-9]: Volatilits moyennes54
Tableau [V-10]: Fractions des phases liquide et vapeur de l'alimentation
chaude.56
Tableau [V-11]: Verification du taux de vaporisation57
Tableau [V-12]: Calcul Mmoy de la charge Chaud ltat liquide et vapeur59
Tableau [V-13]: Calcul Mmoy du distillat a ltat liquide et vapeur.. 61
Tableau [V-14]: Reprsentant la temprature et la pression critiques moyennes
du Distillat.64
Tableau [V-15]: Calcul de la masse volumique du constituant au sommet de la
Colonne..65
Tableau [V-16]: Tableau rcapitulatif et de comparaison entre les dimensions
relles et calcules.67
Tableau [VI -1]: Influence de taux reflux froid sur le fonctionnement du colonne
de Stabilisation 73
Tableau [VI -2]: Les paramtres de fonctionnement de la colonne
P10-CB-21-01...79
Glossaire
1. TERMES GENERAUX
1.1. Point initial de distillation
Temprature releve au moment o la premire goutte de distillat tombe de

lextrmit du condenseur lors dune distillation effectue selon une mthode
normalise.
1.2. Point final de distillation
Temprature maximale releve au cours de lessai de distillation effectu suivant

une mthode normalise.
1.3. Intervalle de distillation
Intervalle de temprature qui dfinit une coupe par son point initial et son point
final de distillation.
1.4 .Coupe
Distillat obtenu entre deux tempratures dtermines au cours dune opration de

distillation.
1.5. Point de bulle
Temprature laquelle apparat la premire bulle de vapeur.
1.6 .Point de rose
Temprature laquelle apparat la premire goutte de liquide une pression fixe

en refroidissant progressivement.
1.7. Reflux
Quantit dun produit issu de distillation servant contrler la temprature

diffrentes zones de la colonne de distillation.
2. PRODUITS ET MATIERES PREMIERES

2.1. Ptrole brut
Huile minrale existant ltat naturel (avant raffinage et autre traitement),

constitue essentiellement dhydrocarbures.
Le ptrole brut peut tre de nature paraffinique, asphaltique ou mixte selon la

prsence dans le rsidu aprs distillation de paraffines, de bitumes ou dun mlange
des deux.
2.2 .Gaz naturel
Hydrocarbures gazeux (principalement du mthane) provenant de gisements

souterrains et dont la production peut tre associe celle du ptrole brut.
Ce gaz est dsign comme humide ou sec suivant la proportion

dhydrocarbures facilement condensables quil contient. Ce terme sapplique
galement au gaz pur.
2.3. Distillat de tte
Fraction la plus volatile du produit trait passant en tte de distillation.
3. PROPRIETES ET ESSAIS
3.1. Tension de vapeur
Pression exerce par les vapeurs mises par un produit ptrolier dans un appareil
appropri et dans des conditions normalises.
La tension de vapeur, une temprature donne, mesure la tendance des molcules

schapper dune phase liquide pour engendrer une phase vapeur en quilibre
thermodynamique.
Lessai de tension de vapeur est trs important, car il indique dune manire directe, la
teneur en produits trs lgers qui conditionnent la scurit au cours du transport.
4.2. Stabilisation
Procd de sparation des fractions lgres du ptrole brut en vue den abaisser la
tension de vapeur.
4.3. Volatilit
A lquilibre, la concentration de tout composant prsent dans la phase liquide

peut tre relie sa concentration dans la phase vapeur par le cfficient de partage
appel volatilit absolue .
4.4.Engorgement par entranement
Si le dbit du gaz augmente trop, un moment, le mlange de liquide et de gaz

peut former une mulsion. Cette mulsion peut remplir tout lespace entre les
plateaux et atteindre le plateau suprieur. Cela augmente fortement la perte de charge
et diminue lefficacit du plateau suprieur, et donc de la colonne.
4.5. Entranement
On parle de phnomne dentranement lorsque le liquide est transport par le gaz

dun plateau au plateau suprieur. Lentranement des gouttelettes du liquide sur le
plateau suprieur est d une vitesse trop leve du gaz dans la zone de contact des
phases.
Lefficacit du plateau est rduite car il y a une diminution des forces motrices
dabsorption due au versement du liquide du plateau moins volatile sur le plateau de
plus grande volatilit. Dans ce cas, la colonne sapproche du point dengorgement, car
la quantit du liquide sur le plateau suprieur augmente.
Introduction gnrale
INTRODUCTION GENERALE
Lnergie occupe une place prpondrante dans la vie de ltre humain. Les
sciences et les progrs techniques ont permis lhomme de dcouvrir des nouvelles
ressources nergtiques, savoir tous les produits vocation nergtique dont nous
citons :
Le ptrole et le gaz naturel, qui sont les seuls capables de rpondre laccroissement
des besoins en nergie.
LAlgrie est riche en ptrole brut. Dans certain cas, ce ptrole contient des
quantits importantes de compos hydrocarbonique plus lger tel que le mthane. Il
faut alors de sparer ces hydrocarbures lgers et de les employer pour dautres
utilisations.
Lusine de lUTBS (unit de traitement de brut sud) comprend trois units de
traitement dhuile.
Lunit de traitement permet datteindre les spcifications requises pour

lexpdition en termes de TVR (entre 7 et 10 psi selon la temprature de stockage
de lhuile stabilise), de salinit (40 mg/l) et de teneur en eau (0,1% volume).
Le but de ce travail est la vrification et loptimisation des paramtres de

fonctionnement dune colonne de stabilisation du ptrole brut , la rectification permet
datteindre une meilleure rcupration en hydrocarbures intressants (iC5-nC12)
avec une bonne puret et limination des composants lgers haute tension de
vapeur (C1 , C2 , C3 et C4).
Notre travail est compos de 6 chapitres :
Historique de la rgion de Hassi Messaoud.
Gnralits sur le brut.
Description gnrale de LUTBS.
Thorie de la rectification.
Calculs des paramtres de la colonne.
Simulation et optimisation des paramtres de la colonne.
Et finalement en termine par une conclusion gnrale.
Page 1
CHAPITRE I
HISTORIQUE DU CHAMP
DE HASSI MESSAOUD
CHAPITRE I Historique du champ de Hassi Messaoud
I.1.Historique
Aprs la mise en vidence par la sismique rfraction de la structure de Hassi Messaoud
comme tant un vaste dme structural, la socit SN REPAL implant le 16 janvier
1956, le premier puit Md1, pas loin du puits chamelier de Hassi Messaoud. Le 15 juin de la
mme anne, ce forage a recoup des grs cambriens 3338 m comme tant producteurs
dhuile.
En mai 1957, la socit CFPA ralise un puits Om1 environ 7Km au Nord-Nord-
Ouest du puits Md1, ce forage confirmait l'existence d'huile dans les grs du Cambrien.
Par la suite, le champ de Hassi Messaoud fut divis en deux concessions distinctes :
C.F.P.A pour le champ Nord. SN.REPAL pour le champ Sud.
La mise en production avait commenc en 1958 avec 20 puits d'exploitation.
Ds lors, les forages se sont dvelopps et multiplis sans cesse, jusqu' arriver plus de
950 puits. Aprs plusieurs annes de production, la pression du gisement a normment
chute ce qui a incit utiliser les mthodes de rcupration secondaire (injection de gaz,
deau, fracturation, acidification, etc...).
I.2. Situation du champ de Hassi Messaoud
I.2.1. Situation gographique
Le champ de Hassi Messaoud est considr comme l'un des plus grands gisements
dans le monde. Il fait partie d'un ensemble de structures formant la partie Nord de la
province Triasique, et se situe environ 850 km au Sud-Sud-Est d'Alger, 280 km au Sud-
Est du gisement de gaz-condensat de Hassi R'Mel et 350 km l'Ouest de la frontire
tunisienne (Fig. I.1). Il stend sur 2500 km, il a pour coordonnes Lambert :
X = [790.000 - 840.000] Est.
Y = [110.000 - 150.000] Nord.
Et il est encadr par les latitudes 31.30et 32.00 et les longitudes 5.40et 6.20.
Page 2
Fig. I.1 : Situation gographique du champ de Hassi Messaoud.
I.2.2. Situation gologique

Par rapport aux gisements, le champ de Hassi Messaoud est limit (Fig. I.2) :
au Nord-Ouest par les gisements de Ouargla (Guellela, Ben Kahla et Haoud
Berkaoui).
au Sud-Ouest par les gisements d'El Gassi, Zotti et El Agreb.
Page 3
au Sud-Est par les gisements Rhourde El Baguel et Mesdar.
Fig. I.2 : Rpartition des zones ptrolifres et gazifires.
Le mle de Hassi Messaoud est le rsultat dune histoire palotectonique assez

complique, cest le prolongement de la dorsale dAmguid El Biod de plus de 800 km de
long. Sa structure fait partie dun ensemble de structures formant la province triasique
Nord orientale.
Gologiquement, il est limit :

l'Ouest par la dpression d'Oued Mya.
au Sud par le mle d'Amguid El Biod.
au Nord par la structure Djamma-Touggourt.
l'Est par les hauts-fonds de Dahar, Rhourde El Baguel et la dpression de
Ghadames.
I.3. Zones et numrotation des puits
Lvolution des pressions des puits en fonction de la production a permis de subdiviser
le gisement de Hassi Messaoud en 25 zones dites de production, dextension variable. Ces
zones sont relativement indpendantes et correspondent un ensemble de puits
Page 4
communiquant entre eux et non pas avec ceux des zones avoisinantes, et se comportant de
la mme manire du point de vue pression de gisement. Les puits dune mme zone
drainent conjointement une quantit dhuile en place bien tablie. Toutefois, il est
important de souligner que le facteur pression ne peut tre le seul critre de caractrisation
des zones (Fig. I.3).
Fig. I.3 : Zones de production.
I.4. Intrt ptrolier
Le champ de Hassi-Messaoud de par sa superficie et ses rserves est considr

parmi les plus grands gisements du monde avec une pression de gisement variant de 120
400 kgf/cm, une temprature de lordre de 118C 123C et une permabilit trs faible
de 0 1darcy.
Le rservoir est li au grs-quartztiques du Combrien, lhorizon le plus productif

est li aux lithozones Ra et RI dont les qualits ptrophysiques sont assez bonnes. La
couverture est assure par une paisse et tanche assise argilo-salifre du Trias. Le champ
Page 5
de Hassi-Messaoud est considr comme une mosaque de gisement, dlimite par les
barrires de permabilit.
Aprs plus de 40 annes de production et plus dune centaine de puits fors, la logique
gologique de ce champ demeure nigmatique, cest ainsi que certaines zones produisent
en dpltion et dautres sont en rcupration secondaire par injection de gaz et deau.
Page 6
CHAPITRE II
GENERALITES SUR LE
PETROLE BRUT [12]

CHAPITRE II Gnralits sur le ptrole brut
II.1. Gnralits sur le ptrole brut

La dcouverte du ptrole date de la plus haute antiquit; considr sous langle de son
emploi comme produit nergtique, le ptrole est dabord un carburant. Il a permis
lindustrie de se dvelopper et laviation de prendre son essor [1].
II.2. Dfinition du ptrole brut

Le ptrole, vient du mot latin petra pierre et oleum huile (soit huile de
pierre ). Cest un liquide brun plus ou moins visqueux d'origine naturelle, mlange
complexe dhydrocarbures. Son exploitation est lun des piliers de lconomie industrielle
contemporaine, car il fournit la quasi totalit des carburants liquides. Le ptrole est aussi
souvent appel or noir en rfrence sa couleur et son cot lev[2].
II.3. Origine du ptrole brut
Le ptrole, est une roche liquide carbone, ou huile minrale de couleur variant du vert
au brun noir et dgage une odeur caractristique. nergie fossile, son exploitation est lun
des piliers de lconomie industrielle contemporaine, car il fournit le quasi totalit des
carburants liquides.
Il sest form dans des lagunes peu profondes par dcomposition de matires
organiques accumules en milieu anarobique. Celles-ci se sont transformes en gaz et en
liquide plus ou moins visqueux qui est le ptrole ; ces transformations sont trs lentes. En
raison de sa faible densit, il a tendance migrer jusqu tre pig dans une couche
poreuse qui donne alors la roche rservoir, coince entre deux couches impermables.
Cest l o se trouvent les gisements [3].
Il existe une grande diversit de ptroles bruts aussi bien du point de vue physique que
du point de vue de la composition chimique ce qui entrane des conditions spcifiques de
production, de transport, de stockage et de raffinage selon les ptroles et les produits qui en
drivent, do la ncessit dune connaissance, la plus exacte possible, de leur
composition [4].
II.4. Classification des ptroles bruts

Il est possible de distinguer les diffrents types de ptrole selon leur densit, leur
fluidit, leur teneur en soufre et autres impurets (vanadium, mercure et sels) et leur teneur
Page 7
en diffrentes classes dhydrocarbures. Le ptrole est alors paraffinique, naphtnique ou

aromatique[5].
Il est aussi possible de les classifier parfois selon leur provenance (golfe Persique,
mer du Nord, Venezuela, Nigeria), car le ptrole issu de gisements voisins a souvent des
proprits proches.
II.4.1. Classification industrielle
Les ptroles sont rpartis dans la partie industrielle en : lgers, moyens et lourds,
dtermins par simple mesure de densit tels que[5] :
Ptrole lger avec : d415 < 0.828
Ptrole moyen avec : 0.828 < d415 < 0.884
Ptrole lourd avec : d415 > 0.884
II.4.2. Classification chimique
Cette classification divise le ptrole selon la prdominance de lune des familles
dhydrocarbure, ils peuvent tre[5] :
Paraffiniques
Naphtniques
Aromatiques
Mixtes : Paraffino-aromatiques
II.4.3. Classification conomique (technologique)
La mthode est base sur la connaissance de la[5] :
Teneur potentielle en soufre dans le ptrole et les fractions ptrolires
Teneur en fractions claires avec PF=350C.
Teneur potentielle en fraction dhuiles de base et leurs qualits.
Teneur en paraffines.
II.4.3.a. Daprs la teneur en soufre
Le ptrole se divise en trois classes[5] :
A faible teneur, lorsque le ptrole ne contient pas plus de 0.5%
A moyen teneur, lorsque le ptrole contient entre 0.5-2%
Sulfureux lorsque sa teneur >2%
II.4.3.b. Selon la teneur potentielle en fractions claires PF=350c [5]
Teneur leve PF = 350c >45% (ptrole lger)
Teneur moyenne entre 30 et 45% (ptrole moyen)
Teneur faible < 30% (ptrole lourd)
Page 8
II.5. Composition des ptroles bruts

Le ptrole brut est un mlange complexe dhydrocarbures de diffrentes familles
(paraffinique, naphtnique, aromatique) associs des composs organiques (oxygns,
azots, sulfurs) ainsi qu des traces de mtaux particuliers (Vanadium, Molybdne,
Nickel) [3].
Les hydrocarbures constituent les lments essentiels des ptroles ; leurs molcules ne
contiennent que du carbone et de lhydrogne dont les teneurs sont respectivement (83%-
87%) et (11%-14%) ; elles se divisent en plusieurs familles chimiques selon leurs
structures [3].
Toutes ces structures sont bases sur la ttravalence du carbone.
La composition du ptrole dpend du lieu o il est extrait, il y a actuellement une

centaine de bruts diffrents sur le march ptrolier caractriss par leurs densits (lourds,
moyens ou lgers) souvent exprimes en degrs API (American Petroleum Institute). On
distingue trois catgories dhydrocarbures prsents dans les fractions des ptroles bruts :
les alcanes ou paraffines (18 65%), les cycloalcanes ou naphtnes (25 90%) et des
composs aromatiques (jusqu 15%). Il ny a pas dalcnes (olfines) ni dalcynes [3].
Page 9
CHAPITRE III
DESCRIPTION DE
LUNITE [6]
CHAPITRE III Description de lunit
III.1. Introduction
Un gisement produit un mlange dhydrocarbures lourds et lgers sous forme
liquide et gazeuse auxquels peuvent sajouter dautres constituants non hydrocarbons (tels
que de lazote, du gaz carbonique, de lhydrogne sulfur, ), de leau plus ou moins
sale, ventuellement des solides. Cet effluent ne peut tre vendu tel quel.
Les transporteurs ainsi que les clients (raffineurs, gaziers, ) demandent des produits
satisfaisant certaines spcifications, en particulier vis ou nuisibles.
Les installations de production surface ont pour objectif de raliser les diffrentes
oprations de traitement des effluents de puits afin de respecter les spcifications
demandes par le client. Ceci englobe :
- la collecte des effluents depuis les puits jusquau centre de traitement.
- les installations de traitement des effluents des puits proprement dites et leurs liaisons.
- ventuellement le stockage des produits aprs traitement.
- les ouvrages dvacuation des produits finis jusquaux points de livraison.
Pour cela, SONATRACH a dcide de construire le complexe dUTBS [6].
III.2. Vue gnrale de LUTBS

L'UTBS est une nouvelle unit de traitement de brut destine recevoir et traiter
lhuile non stabilise provenant de six champs satellites existants de la rgion dHassi-
Messaoud Sud et dexpdier lhuile stabilise vers le centre de stockage situ Haoud El
Hamra via le pipeline 24 CIS-HEH.
Le traitement consiste dessaler et stabiliser le brut provenant des satellites. Le gaz

associ est utilis comme gaz combustible, lexcdant tant envoy vers lunit de GPL
situe au CIS. Leau ncessaire au lavage du brut provient de deux puits deau dAlbien.
Les eaux huileuses traites sont envoyes vers lunit dinjection deau situe OMN77
ou, en secours, vers le bassin dvaporation [6].
LUTBS comprend les units suivantes :
Un rseau de collecte de brut (6 pipelines) pour acheminer lhuile non stabilise des
satellites existants vers la nouvelle installation.
Page 10
Trois units de traitement de brut afin de dessaler et stabiliser le brut provenant de ces
satellites.
Une unit de compression comportant 4 trains.
Une unit de traitement dhuile hors-spec et un bac de stockage dhuile hors-spec.
Quatre bacs de stockage dhuile stabilise.
Une pomperiez dexpdition de lhuile stabilise.
Un systme de comptage dhuile stabilise.
Une unit de traitement des eaux huileuses.
Des pipelines dexpdition (huile stabilise, gaz associ, eau traite).
Les systmes utilits ncessaires aux units de procds.
Page 11
Fig. III.1 : Plan gnral de lUTBS
Page 12
III.3. Description de lunit

III.3.1. Units de Production
LUTBS comprend trois units de traitement dhuile identiques pouvant produire
chacune 100 000 barils par jour dhuile stabilise partir du brut provenant des satellites
existants. Les units de traitement permettent datteindre les spcifications requises pour
lexpdition en termes de TVR (entre 7 et 10psi selon la temprature de stockage de lhuile
stabilise), de salinit (40mg/l) et de teneur en eau (0.1% volume).
En fonctionnement normal, le brut non stabilis provenant des satellites existants est
rparti vers les units de traitement de brut de lUTBS. Si le dbit dhuile provenant des
satellites est suprieur la capacit des units de traitement en service, lexcdent est
dirig vers lunit de traitement dhuile hors-spec [6].
III.3.2. Rseau de Collecte Production de Satellite

Des puits producteurs alimentent les plateformes satellites existantes o une premire
sparation huile/gaz/eau est ralise. Au niveau de chaque plateforme satellite, le brut non
stabilis est pomp du sparateur triphasique vers la nouvelle unit UTBS via le nouveau
rseau de collecte.
Un ensemble de connexions (tie-ins) est ralis afin de relier les plateformes satellites
existantes au nouveau rseau de collecte [6].
Le nouveau rseau de collecte permet dacheminer le brut non stabilis des plateformes
satellites existantes vers le manifold M01 situ lentre de lUTBS :
le satellite W1C est connect au manifold M01 par un pipeline enterr de 20

le satellite W1A est connect au manifold M01 par un pipeline enterr de 12
le satellite W2A est connect au manifold M01 par un pipeline enterr de 12
le satellite E1C est connect au manifold M01 par un pipeline enterr de 14
le satellite E2A est connect au manifold M01 par un pipeline enterr de 10
le satellite S1A est connect au manifold M01 par un pipeline enterr de 10
Page 13
CHAPITRE III Descripti
tion de lunit
Fig.. III.2
I : Rseau de Collecte du Brut
III.3.3. Traitement et Stabilisat

ation dHuile
3.3.a. Sparation dHuile
La sparation huile/eau/gaz se fait sur deux tages, le premier tage de ssparation est
constitu dun ballon de sparat
ration triphasique et le second tage dun ballon
n biphasique.
b
Le brut non trait arrivant de M01 est envoy dans le premier tage d
de sparation
(sparateur tri-phasique). Les
es gaz sont spars puis envoys vers lee systme de
compression. Leau est envoye
e vers le traitement des eaux huileuses.
Lhuile est ensuite rchauffe ju

jusqu 70C dans un changeur brut non trait/
t/brut stabilis
constitu de deux calandres dis
disposes en srie. La temprature de sortie estt contrle
c par
lintermdiaire dun by-pass cot
co calandre. Lhuile rchauffe alimente le deu
deuxime tage
de sparation (sparateur bi-ph
phasique). Les gaz spars sont galement env
nvoys vers le
systme de compression [6].
Le brut est ensuite renvoy par

ar une pompe (2 x 100%) vers le systme de dessa
ssalage.
Page 14
Vers rseau torche BP

Du prchauffeur de la
colonne de stabilisation
TC
Du manifold de collecte
dhuile brute
Vers compression
PC
FC
PC
PX0-GA20-01B
PC PX0-VA-20-02
LC
PX0-VA-20-01
LC
PX0-GA20-01A
LC
Vers arorfrigrants
dhuile stabilise
FC Vers traitement des

eaux huileuses
Vers Dessaleurs
P10-PA-20-01 A/B
Fig. III.3 : Sparation dHuile
3.3.b. Dessalage
Le systme de dessalage est constitu de deux dessaleurs lectrostatiques en srie avec
injection deau de lavage lentre du deuxime dessaleur. Leau extraite du deuxime
tage sous contrle de niveau dinterface est entirement recycle vers lentre du premier
tage par lintermdiaire dune pompe (2 x100%). Leau extraite du premier tage est
envoye sous contrle de niveau dinterface vers lunit de traitement des eaux huileuses.
Chaque dessaleur lectrostatique est quip dune vanne de mlange son entre.
Leau de lavage provient dun puits deau (albien) for par le Matre de lOuvrage [6].
Page 15
FC
De la distribution deau
de lavage
Du contrle de niveau du
sparateur biphasique
Des pompes %
FY
FC
dalimentation du
dessaleur PDC Vers sparateur hors
spec

FY
%
LC LC FY FC
PDC
PX0-VW-21-01 PX0-VW-21-02 Vers colonne de
stabilisation
FC
Vers prchauffeur de la
colonne de stabilisation
PX0-PA-21-03 A/B
PX0-PA-21-02 A/B Du rseau azote BP
Fig. III.4 : Dessalage
3.3.c. Stabilisation
Lhuile dessale alimente la colonne de stabilisation P10-CB-21-01. 25% du dbit
alimente directement la tte de la colonne (alimentation froide) ; 75% du dbit alimente le
prchauffeur de la colonne de stabilisation P10-GA-21-01 A/B (alimentation chaude de la
colonne). Le prchauffeur de la colonne de stabilisation est constitu de deux calandres en
srie. Chaque calandre peut tre by-passe. La temprature de lhuile provenant du
package de dessalage est de 70C environ. Elle passe ct tubes et elle est chauffe par
lhuile stabilise qui passe ct calandre. Afin doptimiser la rcupration de chaleur, la
temprature de lhuile stabilise en sortie du prchauffeur, ct calandre, est rgule
120C.
La colonne de stabilisation oprant 5 barg permet de retirer les composs les plus lgers
du brut et datteindre en fond de colonne la TVR (tension de vapeur Reid) requise pour le
stockage dans les bacs toit flottant. En dautres termes, la colonne de stabilisation permet
daugmenter la temprature de vaporisation de lhuile stabilise au-dessus de la
temprature de sortie des aro-rfrigrants dhuile. La chaleur en fond de colonne est
fournie par le four de rebouillage P10-FA-21-01.
Page 16
La colonne de stabilisation permet datteindre une TVR de 7 psi dans le cas t (cest--
dire pour obtenir un point de bulle de 61C pression atmosphrique) et de 10 psi dans le
cas hiver (cest--dire pour obtenir un point de bulle de 42C pression atmosphrique).
La colonne de stabilisation comporte 21 plateaux rels. 8 plateaux 2 passes sont installs

dans la partie suprieure de la colonne, 13 plateaux quatre passes dans la partie
infrieure. Lalimentation froide liquide est ralise en tte de colonne, au-dessus du
plateau 21 ; lalimentation chaude mixte est ralise au-dessus du plateau 13.
Afin de limiter la prsence de sel dans la colonne, des soutirages deau sont installs sur les
plateaux 11 et 19. Ils permettent de collecter leau qui scoule par diffrence de densit
dans des pots de rcupration. Leau accumule est draine manuellement vers le package
de traitement des eaux huileuses.
Le liquide descendant du plateau 1 est mlang lhuile stabilise en fond de colonne qui
ne comporte pas de cloison.
Une partie de lhuile en fond de colonne alimente le four de rebouillage P10-FA-21-01 par
les pompes de recirculation du rebouilleur P10-PA-21-01 A/B/C (3 x 50%). Le four
apporte la chaleur ncessaire la stabilisation et permet la vaporisation dune partie de
lhuile stabilise. Le retour vers la colonne du mlange biphasique en sortie du four de
rebouillage est ralis sous le plateau 1. La vapeur alimente le plateau 1 alors que le liquide
est mlang au liquide stabilis en fond de colonne et en vaporise une partie.
Le four de rebouillage est un four tirage naturel quatre passes. Le dbit vers les passes
est quilibr manuellement. Le four dispose de six brleurs et six pilotes. Le gaz provient
du systme de gaz combustible. Un pot de rcupration est install en amont du package de
gaz combustible pour collecter le liquide ventuellement prsent dans le gaz. Chaque four
dispose de son propre ballon dazote haute pression pour permettre ltouffement du four
en cas dincendie.
Lhuile stabilise chaude (entre 135 C dans le cas hiver et 160C dans le cas t) sort de la
colonne et passe ct calandre dans le prchauffeur de la colonne de stabilisation puis ct
calandre dans le rchauffeur dhuile et cde ainsi de sa chaleur lhuile non stabilise. Le
refroidissement final avant stockage est assur par les rfrigrants dhuile stabilise P10-
GC-21 A/B/C/D, constitus de 4 baies en parallle, chaque baie comportant deux
ventilateurs, dont un pales variables [6].
Page 17
TC
Vers compression
Du dessaleur
CHAPITRE III
PX0-CB- 21-01
LC
FC
PX0-GA-21-01 B
PX0-FA-21-01
PX0-GA-21-01 A
Du dessaleur PX0-PA-21-01 A/B/C PC

Du rseau de gaz
combustible
TC
Fig. III.5 : Stabilisation

Vers sparateur
biphasique TC
PX0-GA-20-01 B Vers gaz boot
PX0-GA-20-01 A
TC
Vers stockage huile
PX0-GC-21-01 A/B/C/D 1/2
Du sparateur stabilise
triphasique
Description de lunit
Page 18
III.3.4. Stockage, Comptage et Expdition de lHuile Stabilise

Le brut stabilis est ensuite stock dans 4 bacs toit flottant ayant chacun une capacit
brute de 50 000 m3.
La pomperiez dexpdition est compose de pompes de gavage (4 x 33%) et de pompes

dexpdition (4 x 33%).
Un comptage transactionnel (et non fiscal) est implant avant la ligne dexpdition vers
CIS.
La ligne dexpdition de brut stabilis de 30" partir de lUTBS se raccorde sur la

nouvelle ligne 24 pouces CIS-HEH ralise par le Matre de lOuvrage. Le point de
raccordement est situ lintrieur du CIS, proximit des pompes dexpdition existantes
et proche de la gare de racleurs dpart vers HEH (30 m environ). Le raccordement se fait
sur la partie arienne de la canalisation, avant son enfouissement. A noter que le point de
raccordement sur le pipe dexpdition vers HEH, est quip dune vanne 24" srie 300
existante.
La ligne est quipe de gares de racleur dpart/arrive et pourvu de vannes de

sectionnement [6].
III.3.5. Traitement et du stockage de lhuile hors spec
Le systme dhuile hors-spec est utilis de faon exceptionnelle lors du dmarrage de

linstallation ou du dclenchement dune ou plusieurs units de traitement dhuile, du
dclenchement du systme de stabilisation (colonne), du dclenchement du bac de
stockage en remplissage ou de TVR trop leve. Celui-ci permet de continuer la
production en stockant de manire sre de lhuile mal stabilise. La production dhuile au
niveau des plateformes des satellites peut donc continuer tant que le bac de stockage
dhuile hors spec nest pas compltement rempli.
Le bac de stockage dhuile hors-spec est galement utilis en opration normale. Il reoit
en effet le liquide collect dans le ballon de drains ferms et dans les ballons de torche,
ainsi que les effluents huileux du package de traitement des eaux huileuses et du traitement
API [6].
Page 19
PC
Vers rseau torche

TBP
Du rseau de gaz
combustible
P90-VG-21-01
PC
Du rseau de gaz
combustible
FC
De la sortie des P90-RM-21-01
dessaleurs
PC
Du collecteur dhuile
hors spec basse
pression
P90-VA-21-01 LC
Des arorfrigrants
dhuile stabilise
Vers traitement des

eaux huileuses
P90-PG-21-02
Du collecteur dhuile
hors spec
Vers sparateur
triphasique train 10
FC
P90-PA-21-01 A/B
Vers sparateur
Vers sparateur
Vers pompes
dexpdition de lhuile
stabilise
Fig. III.6 : Systme dHuile Hors-Spec
III.3.6. Unit de Compression de Gaz

Le gaz associ gaz de flash provenant des sparateurs triphasiques, des sparateurs
biphasiques et des colonnes de stabilisation des units de traitement dhuile est envoy
dans le collecteur commun de gaz de flash oprant 4.4 barg. Une partie du gaz de flash
est utilis au sein de lUTBS comme gaz combustible. Lexcs de gaz est comprim et
expdi vers lunit de GPL situe au CIS.
Le systme de compression de gaz de flash consiste en quatre trains de compression

identiques (trois trains en opration, un train en attente). Chaque train comporte deux
compresseurs centrifuges (compresseur de 1er tage et compresseur de 2me tage) installs
sur le mme arbre avec un corps commun et entrains par un moteur lectrique commun.
La pression dans le collecteur commun daspiration est rgule grce aux vannes de
laminage situes lentre de chaque train et au rpartiteur de charge (load sharing) des
trains de compression [6].
Page 20
0G3-VM-27-01 CIS-VM-27-01
Vers Gaz Combustible
de Secours
Vers Systme de rpartition de Charge et de
Contrle de pression dAspiration
Anti pompage et Contrle de Dbit

PC
Injection dinhibiteur
de Corrosion
G51-GC-23-01
G51-VD-23-02
G51-GC-23-02
G51-VD-23-03
Gaz de la Colonne de G51-VD-23-01
Stabilisation
TC LC
Gaz des Sparateurs LC TC
Triphasiques
G51-KA-23-01 M
M
G51-KA-23-02
Gaz des Sparateurs

Biphasiques
Vers Gaz
Combustible
Vers Sparateurs
Vers Drains Ferms Triphasiques
Train B
Train C
Train D
Fig. III.7 : Compression de Gaz
III.3.7. Lunit de traitement des eaux huileuses
LUnit de Traitement des Eaux Huileuses permet de traiter les eaux huileuses de
procd ainsi que les effluents pr-dshuils des drains ouverts avant stockage dans le bac
tampon deau traite puis expdition vers OMN77 pour injection dans un puits [6].
III.3.8. Les systmes utilits [6]
Gaz combustible.
Torches.
Drains ferms.
Eau brute.
Eau potable.
Drains ouverts.
Traitement des eaux huileuses.
Stockage et expdition de leau traite.
Azote.
Traitement de lhuile de lubrification.
Stockage et distribution de diesel.
Injection de produits chimiques.
Page 21
Station dpuration des eaux uses.

Eau incendie.
Emulseur.
Page 22
CHAPITRE IV
THEORIE DE
RECTIFICATION
CHAPITRE IV Thorie
rie de rectification
IV.1. Introduction
Il est bien connu que la rec
ectification consiste simplement dans la format
ation dune phase
vapeur riche en lment le pluss volatile et une phase liquide riche en lmentt le plus lourd que
le mlange initial.
IV.2. Dfinition
La rectification est un pro
rocd physique de sparation dun mlange
ge en constituant
pratiquement pur par des chan
anges successifs de masses et de chaleur entre
re les deux phases
circulant contre courant dans u
un appareil appele : Colonne de Rectification.
La colonne de rectification
n peut tre munie des plateaux successifs ou
o de garnissage
mtallique qui dont le rle cons
nsiste assurer un contact intime entre les pha
hases. Sur chaque
plateau, il y a barbotage de la v
vapeur dans le liquide avec un transfert de mas
asse et de chaleur.
Les plateaux sont garnis dlmeent, appel : Cloches, Clapets, Calottesetc.
Les phases liquides et vapeur
eurs pntrants sur un mme plateau ne sont jam
jamais en quilibre
(tliq<tvap) [7].
A la suite de ce contact intim
me entre les deux phases, le systme volue plus
us ou moins rapide
jusqu atteindre lquilibre. Daans le but de simplifier la conception et le calc
lcul du procd de
rectification, on introduit la notio
tion ou lhypothse du plateau thorique.
Un plateau est dit thorique
ue si les deux phases partantes de ce mmee plateau sont en
quilibre. La rectification seff
effectue toujours dans une colonne munie de plateaux, dun
condenseur partiel et dun reboui
uilleur [7].
F IV.1 : Schma simplifi

Fig.
Rectification
Page 23
CHAPITRE IV Thorie de rectification
IV.3. Description dune colonne de rectification
Chaque colonne de rectification compose toujours de trois zones :
Une zone dalimentation dite (zone de flash) ou dalimentation est introduite ltat
vapeur, liquide ou bi phasique (selon la temprature de la charge) [8].
Une zone situe au dessus de la zone de flash est appele (zone de rectification),
caractrise par un condenseur qui sert condenser partiellement ou totalement la
vapeur sortante de la tte de colonne, une partie de cette vapeur est envoye en tant que
(reflux froid) vers la partie suprieure de la colonne pour garder le gradient de
temprature et corriger la composition des produits de tte. Lautre partie est extraite
comme distillat ltat vapeur [8].
Une zone situe au dessous de la zone de flash est appele (zone dpuisement),
caractrise elle aussi par la prsence dun rebouilleur, qui sert apporter une certaine
quantit de chaleur au liquide soutir du dernier plateau (reflux chaud ou rebouillage)
afin de rcuprer les constituants lgers entrans par les lourds [8].
La quantit de matire change dpend de la concentration initiale et finale, des pressions

partielles, de la surface de contact et de lcoulement (vitesse). Le fonctionnement normal de la
colonne de rectification pour lobtention dun distillat et dun rsidu avec des compositions
voulues dpend principalement de ltat de la charge et par consquent du lieu dalimentation.
Gnralement la charge peut tre envoye dans la colonne sous diffrente forme tats [8] :
Etat liquide (e0=0) :
Lalimentation doit tre introduite dans la partie suprieure de la colonne.
Etat vapeur (e0=1) :
Lalimentation doit tre introduite dans la partie infrieure de la colonne.
Etat bi phasique liquide vapeur (0<e0<1) :

Lalimentation doit tre introduite dans la partie mdiane de la colonne lieu exact de
lalimentation peut tre dtermin en calculant le nombre de plateau de la zone de rectification
et dpuisement.
Page 24
IV.4. Principe de fonctionnement dune colonne de rectification
La colonne de rectification est un appareil qui spare un mlange en deux ou plusieurs

produits qui sont soit des corps purs, soit des mlanges comportant moins des constituants
que le mlange initial [9].
La colonne est une enveloppe cylindrique verticale (la virole) qui renferme des
dispositifs (plateaux ou garnissage) permettant dassurer un contact intime entre les deux
phases (reflux liquide et gazeux)
Afin de permettre un transfert de matire et dnergie entre les deux phases et atteindre
ltat dquilibre. Sur le long de la colonne la phase vapeur senrichit, de bas en haut en
lment le plus volatil et la phase liquide senrichit en lment le moins volatil de haut en
bas. Pour avoir en tte de la colonne un distillat D et en bas un rsidu R.
Fig. IV.2 : Schma d'une colonne plateau
Page 25
IV.5. Types des colonnes de rectification

Les colonnes utilises pour la rectification se divisent en deux types distincts :
les colonnes lments spars (plateaux, grilles, cascades).

Les colonnes garnissage.
IV.5.1. Colonnes plateaux

Le premier moyen de favoriser le contact gaz liquide est de faire buller le gaz dans un
liquide retenue dans une capacit : le plateau. Nous donnons dans la Fig.IV.3 la
reprsentation dune tour plateaux [9].
Fig. IV.3 : Reprsentation dune tour plateaux
Les plateaux doivent tre conus comme des lieux de rencontre optimale du liquide et
de la vapeur. Le liquide introduit est vaporis partiellement dans le bouilleur. La vapeur
remonte jusquau condenseur. Une partie de la vapeur est condense et redescend dans la
colonne contre courant de la vapeur. En distingue plusieurs types plateaux [9] :
Types de plateaux :
a) Plateau calotte (bubble-cap tray)
b) Plateau perfor (sieve trays)
c) Plateau clapets (valve trays)
d) Plateaux sans dversoirs
Page 26
IV.5.2. Colonnes garnissage

Elles sont utilises dans la distillation et l'absorption. Ces colonnes ralisent contre-courant
continu. L'change de matire croit en rgime turbulent, c'est pourquoi elles sont remplies
d'lments solides. On emploi comme garnissage n'importe quels morceaux solides,
gnralement du coke ou de la brique qui ont l'avantage d'tre particulirement bon march [9].
Le garnissage efficace doit tre tel que la surface par unit de volume soit la plus grande
possible; distribution uniforme du liquide et du gaz et faible pertes de charges.
Les garnissages les plus utiliss dans l'industrie sont de (RASHING) : des manchons de
cylindres dont la hauteur est gale au diamtre extrieur .Leur construction est en cramique ou
en porcelaine (bon march et grande rsistance la corrosion). Comme il existe aussi des
anneaux en acier, en carbone et en matire plastique (polythylne) [9].
IV.6. Lois fondamentales de la rectification
IV.6.1. Lois de DALTON
Considrerons une phase vapeur compose de constituants aux quels on peut appliquer la loi
des gaz parfaits, alors la pression partielle dun constituant (pvi) est proportionnelle la pression
du systme (PT) et sa concentration molaire (Yin) [8]:
Pvi = PT Yi (1)
Avec :
pvi : Pression partiel dun constituant i
PT : Pression du systme.
Yi : Concentration molaire du constituant i en phase vapeur.
IV.6.2. Lois de RAOULT et dHENRY
Les lois de Raoult et Henry ne sont applicables que pour les liquides. Considrons une phase
liquide compose de constituants entirement miscibles, On conoit aisment que la tendance
des molcules dun constituant schapper de la phase liquide soit proportionnelle la
concentration molaire de ce constituant cette phase.
Daprs la loi de RAOULT la pression partielle (pLi) dun constituant est gale au produit de
sa tension de vapeur (Pi) ltat pur par sa concentration molaire en phase liquide :
piL = Pi . Xi (2)
Avec :
pLi: Pression partiel dun constituant i
Pi: Tension de vapeur du constituant i
Xi : Concentration molaire du constituant i en phase liquide.
Page 27
Dans le cas des mlanges rels les constituants obissent la loi dHENRY, daprs cette loi,
la pression partielle dun constituant, la temprature tant fixe, est proportionnelle la
concentration molaire du constituant dissout dans la phase liquide :
piL =Ei. Xi (3)
Avec:
Ei : constante dHenry ; dpend de la nature du constituant et de la temprature lquilibre,
les pressions partielles dun mme constituant dans les deux phases liquides et vapeur sont
gaux :
p Li = p Vi (4)
Donc :
Yi = (Pi / PT). X (5)
Cest lquation dquilibre
Avec :
Ki= (Pi / PT) : Coefficient dquilibre du constituant i Pour un mlange de deux
constituants A et B :
yA / yB = (KA / KB). (xA / xB) (6)
i = (KA / KB) (7)
Avec :
i : Coefficient de volatilit relative qui caractrise la capacit de A de passer en phase
vapeur par rapport B [8].
IV.6.3. Bilan matire de la colonne et de ses zones
Le bilan matire de la colonne :
L=D+R (8)
Avec :
L : Dbit massique de la charge.
D : Dbit massique du Distillat.
R : Dbit massique du Rsidu.
Le bilan matire partiel peut scrire comme suit :
L.XL=D.YD+R. XR (9)
La rsolution connue de ces quations, en appliquant la loi de LEVIER donne :
L/(YD XR) = D/( XL XR) = R/( YD XL) (10)
Bilan matire de la zone de rectification :
Page 28
CHAPITRE IV Thorie
Gi= g i-1 + D (11)

Avec :
Gi : Dbit pondral de v
vapeur traversant un sectionnement entre deux
eux plateaux de la
colonne.
gi-1 : Dbit massique du liq
liquide traversant un sectionnement.
D: Dbit massique du dist
istillat.
Donc :
G i.Yi = g i-1. X (i-1) + D. YD (12)
Appliquant la loi de LEVIE
IER, on trouve :
G/ (YD X (i-1)) = g i-1 / (YD Y i) = D/ (Y i X (i-1)) (13)
A partir de cette quation,
n, on trouve :
yi=( r f/( r f+1)) x i-1+ (1
(1/( r f+1)) . yD (14)
On note que : r f =g/D
Avec :
rf : le taux de reflux : Lee reflux froid a pour but la de rguler la temp
prature en tte de
colonne, de corriger composit
sition et de ne permettre quaux lments lgers
l de quitter
lappareillage, le reflux froid prend
pr la chaleur des vapeurs et provoque unee ccondensation des
constituants lourds [8].
Il est gal au rapport de laa quantit
q de reflux renvoy dans la colonne su
sur la quantit du
distillat :
Lquation qui tient compte d

de la constante des dbits massiques G, g d
dans la zone de
rectification de la colonne donne
ne une relation entre la concentration en constit
tituant ( i ) dans le
liquide se trouvant sur le plateau
au (i-1) et la concentration en mme constituant
nt dans la vapeur
provenant du plateau ( i ) et perm
rmet de dterminer le nombre dtages thorique
ues dans la zone de
rectification[8].
Page 29
Q1
D, yD, HD
g0, x0, y0
G1, H1, y1 1 g1, x1, h1

G2, H2, y2 g2, x2, h2
Gi-1, Hi-1, yi-1 2 gi-1, xi-1, hi-1
i-1
Gi, Hi , yi i gi, xi, hi
i+1
Gi+1, Hi+1, yi+1 gi+1, xi+1, hi+1
Fig. IV.4 : Schma descriptif du fonctionnement de la zone de rectification
Bilan matire de la zone dpuisement :

gi+1 = Gi+R (15)
Avec :
gi+1 : Dbit massique du liquide traversant un sectionnement.
Gi : Dbit massique de vapeur traversant un sectionnement entre deux plateaux de la

colonne.
R : Dbit massique du rsidu.
gi+1. x ( i+1) = Gi Y i+ R . XR (16)
En appliquant la loi de Levier, on trouve :
Gi+1/ (Yi - XR) = Gi/ (X ( i+1) - XR) = R/ (Yi - X(i+1)) (17)
A partir de cette quation, on trouve :
Xi+1=( r b/( r b +1))Yi +( 1/ rb+1)XR (18)
On note que : r b =G/R

Avec :
r b : le taux de rebouillage
Page 30
CHAPITRE IV Thorie
Le reflux chaud a pou but dap

apporter la chaleur, de maintenir la tempratur
ture au fond de la
colonne et dassurer la vaporisati
ation des lgers entrans par le produits de fond.
d.
Il est gal au rapport de la quan
antit de reflux chaud provenant du rebouilleur
ur la quantit de
rsidu [8].
Gi+2, Hi+2, yI+2 i+2 gi+2, xi+2, hi+2
Gi+1,Hi+1, yi+1 i+1 gi+1, xi+1, hi+1
Gi, Hi, yi i gi, hi, xi
G2, H2, y2 2 g2, h2, x2
G1, H1, y1 1 g1, h1, x1
R , hR
R, xRi
Q2
Fig. IV.5 : Schma descripti

ptif du fonctionnement de la zone dpuisement
nt
Page 31
Bilan matire de la zone dalimentation :

En se rfrant la figure ci-dessous, on voit quil faut mlanger le liquide gk avec la phase go
pour obtenir le liquide gm alimentant la zone dpuisement, le bilan matire sera effectue juste
au dessus du plateau suprieur de cette zone [8].
L = G0 + g0 (19)
L XL = g0. x0 + G 0 . y0 (20)
Avec :
L, XL: Dbit et concentration molaire de lalimentation.
G0, y0: Dbit et concentration molaire de la vapeur.
g0, X0: Dbit et concentration molaire du liquide.
Zone de rectification
Gm, yD
G0,y0
L, xL, hLi gk,xk
GK,yK g0, x0
gm, xR,
Zone dpuisement
Fig. IV.6 : Schma descriptif du fonctionnement de la zone dalimentation
Le bilan matire de la phase vapeur alimentant la zone de rectification :
G m = G0+GK (21)
G m .y m= G0.y0+GK..yK (22)
Le bilan matire de la phase liquide alimentant la zone dpuisement

g m=g0+g K (23)
g m .Xm = g0 .X0 + g K .XK . (24)
IV.6.4. Utilisation de la volatilit relative

Dans les calculs de la rectification des mlanges complexes, on sera amen effectuer des
sries de calcul dquilibre pour chaque zone.
Pour viter le travail lent et fastidieux des approximations successives, on peut utiliser une
mthode simplifie et rapide donnant une prcision assez suffisante. Cette mthode fait
Page 32
intervenir la notion de la volatilit relative i tant le rapport des coefficients dquilibre du

constituant i celui dun constituant de rfrence r choisit gnralement comme tant
llment le plus lourd (le moins volatile du mlange) [9].
= k i /k r (25)
IV.6.5. Choix de la pression dans la colonne
Quelque soit la complexit du mlange sparer, le choix de la pression est avant tout un
problme conomique. On sait qu faible pression la slectivit de la sparation sera meilleur et
que lon vitera laltration thermique des produits, mais on contrepartie, faible pression il
faut condenser basse temprature les vapeurs de tte pour obtenir le reflux liquide sa
temprature de bulle, donc cest le choix du fluide rfrigrant qui va imposer la pression dans la
colonne. En gnrale, on utilise une source froide, bon march et abondante (eau douce, eau de
mer ou lair) qui permet dobtenir un condenst de 30 450C environ [9].
La pression dans la colonne sera alors au moins gale la pression de bulle ces
tempratures.
Lorsque le distillat est constitu de constituants trs lgers, il est ncessaire de faire appel
une rfrigration artificielle pour travailler basse temprature quelque soit la pression de
service.
La temprature de la condensation totale des vapeurs du distillat doit tre prise moins 15
20oC plus grande que celle du fluide rfrigrant, ainsi la pression dans le ballon de reflux Pb est
dtermine daprs lquation de lisotherme de la phase liquide, par approximation successive.
K i .xi =1 (26)
La pression au sommet, au fond et dans la zone dalimentation de la colonne sera

respectivement :
Ps=Pb+P1 (27)
Avec :
P1 : (0,2 : 0,4) atm, perte de charge due aux rsistances des conduites de condenseur
Pb : cest la pression dans le ballon de reflux.
Pf=Ps+P2 (28)
Avec:
P2 : (0,3 : 0,5) atm, perte de charge due aux rsistances des plateaux
Par consquent, la pression dans la zone dalimentation sera :
PL = (Ps+Pf) / 2 (29)
Page 33
IV.6.6. Rgime de la temprature dans la colonne

A la pression de service Ps, la temprature au sommet de la colonne est dtermine par
approximation successive partir de lquation de lisotherme de la phase vapeur :
(Yi /Ki ) = 1 (30)
La temprature au fond de la colonne (dans le rebouilleur doit tre calcul comme la

temprature de bulle du Rsidu la pression Pf) est dtermin par approximation successive
partir de lquation de lisotherme de la phase liquide.
(Ki . x,i ) = 1 (31)
Pour la temprature de la charge, elle peut tre dtermine en trois cas :
Si elle est ltat liquide, on utilise lquation suivante :
(Ki . xL,i ) = 1 (32)
Si elle est ltat vapeur, on utilise lquation suivante :
( yL i /Ki ) = 1 (33)
Si elle est ltat liquide vapeur, on utilise lquation suivante :
xi = [ x L i / 1 + e(Ki-1) ] = 1 (34)
(yi / xi ) = Ki (35)
Avec:
e : Taux de vaporisation massique.
X0.i : Concentration molaire du constituant (i) dans la phase liquide.
Y0.i : Concentration molaire du constituant (i) dans la phase vapeur.
XL.i : Concentration molaire du constituant (i) dans la charge.
IV.6.7. Taux de vaporisation et composition des phases liquides et vapeurs de la charge

Lors du calcul de la zone dalimentation, il est ncessaire de disposer du dbit et de la
composition des phases liquides et vapeurs de la charge.
Daprs la dfinition de taux de vaporisation et de condensation de la charge,
Page 34
on a :
e = G0 / L (36)
r = g0 / L (37)
Avec :
e: Taux de vaporisation massique de la charge.

r : Taux de condensation massique de la charge.
go: Dbit massique de la phase liquide dans la charge.
Go: Dbit massique de la phase vapeur dans la charge.
A la pression et la temprature donne, le taux de vaporisation ainsi que la composition des
phases de la charge sont dtermines par approximation successive daprs lquation connue
de TREGOUBOV.
X i = [XLi /1+ e.(Ki-1)] = 1 (38)
Ou :
r+e=1 (39)
Cette dernire quation peut tre utilise galement pour dterminer la temprature de la
charge si l a pression et (e) sont connus.
Remarque :
Si e = 0 : la charge se trouve sa temprature de bulle et (go = L), on utilise lquation

disotherme de la phase liquide pour dterminer cette temprature.
Si e = 1 : la charge se trouve sa temprature de ros et (G0 = L), on utilise lquation
disotherme de la phase vapeur pour dterminer cette temprature.
Pour calculer le taux de vaporisation en fonction de la temprature, on peut utiliser la relation
approximative suivante :
e = ln(T / Tb) / ln(Tr / Tb) (40)
O :
T: la temprature considre (k).
Tb : la temprature de bulle (k).
Tr : la temprature de ros de la charge (k).
IV.6.8. Choix du taux de reflux
Le taux de reflux dtermine les dimensions de lquipement et les dpenses nergtiques au
sommet et au fond de la colonne.
Page 35
IV.6.8.a. Fonctionnement reflux minimal
Le fonctionnement reflux minimal correspond un nombre infini de plateaux pour la

sparation spcifie, donc daprs la formule dUNDER WOOD : gnralement, on dtermine
le taux de reflux minimal suivant lquation :
(r f)min = (YD / i,L - ) (41)
(r b)min = ( i,L . XR / i,L-) (42)
e =Li .XLi / (Li . ) (43)
Avec :
i,L : Le coefficient de volatilit du constituant i quelconque par rapport au

constituant cls lourd dtermin la temprature dentre de la charge.
: Paramtre conventionnel dterminer par approximations successives.
La relation entre (r f) min et (r b) min est donne par la formule suivante :
(rb)min = [ (D(rf)min / L) + (1- eo) - (R / L) ] / (R / L) (44)
IV.6.8.b. Taux de reflux optimal
En suite, on dtermine le taux de reflux opratoire ou optimal suivant qui correspond un

nombre de plateaux thoriques (chaque taux de reflux correspond un nombre fini de plateaux
thorique et que ce nombre diminue quand le taux de reflux saccrot).
Pour ces calculs approximatifs, on peut utiliser la formule propose par GILLILAND [9] :
[(r f) opt - (r f) min] / [(r f) opt + 1] = 0,10,33 (44)
(rb)opt = [D (rf)min / R) + (1-e)] - 1 (45)
IV.6.9. Bilan thermique de la colonne
Ltablissement du bilan thermique de la colonne peut tre ncessaire pour dterminer :

La quantit de chaleur apporter au rebouilleur.
La quantit de chaleur retirer au condenseur.
Ainsi que pour vrifier :
Le dbit de fluide chauffant et le dbit de fluide refroidissant.
Le dbit de reflux.
Page 36
La somme des flux enthalpiques entrants doit tre gale la somme des flux enthalpiques
sortants.
Les flux enthalpiques entrant sont gnralement :
La chaleur apporte au rebouilleur.
Le flux enthalpique des produits alimentant la colonne.
Les flux enthalpiques sortant sont gnralement :
La chaleur retire au condenseur.Le flux enthalpique des produits extraits de la
colonne.
Lquation de bilan thermique de la colonne entire scrit :
L . hL + Q r b = D. HD + R .h R + Q Con (46)
L . hL = D(H D+Q con/D)+R(h R Q r b/R) (47)
L. h L = D.H*D +R . h* R (48)
Avec:
Q con /D : la charge thermique de condenseur rapporte lunit de distillat.
Q r b /R : la charge thermique de rebouilleur rapporte lunit de rsidu.
Q r b : Charge thermique du rebouilleur. (Kcal/h).
Q Con : Charge thermique du condenseur. (Kcal/h).
hL : Enthalpie de la charge (Kcal/Kmol).
HD : Enthalpie du distillat (Kcal/Kmol).
hR : Enthalpie du rsidu (Kcal/Kmol.).
IV.6.10. Nombre de plateau thorique reflux opratoire
Mthodes approximatives
GILLILAND (II) a propos une corrlation empirique entre le taux de reflux et le nombre
de plateaux thorique. Si N est le nombre de plateaux thoriques correspond un taux de reflux
opratoire (rf) les deux fonctions sont relies entre elles par une courbe de corrlation[9].
(N N min) / (N+1) = ( rf (rf)min) / ( rf +1) (49)
Page 37
La formule de Fenske- Under Wood , nous donne la relation entre la concentration et

les volatilits relatives de deux constituants dans le distillat et dans le rsidu.
Le nombre de plateaux minimal :
Nmin = log X DV . X RL / log V (50)

X DL X RV L
Avec :
V : Volatil
L : Lourd
Le nombre de plateaux thorique dans la zone de rectification et dpuisement est fixe par
lapplication de lquation de Fenke- Under Wood Aux constituants cls entre
lalimentation et le distillat ou le rsidu [9].
Zone de rectification :
X DV X LL
N Rmin = log . / log V (51)
X DL X LV
L
Zone dpuisement :
N Emin = log X LV . X RL / log V

(52)
X
LL X RV L
Le nombre de plateaux thoriques dans la colonne entire o dans chacune de deux zones,
peut calculer par lapplication de la mthode GILLILAND [9].
NN (rf )op (rf ) (53)

min = f min
N 1 (rf )op + 1

Daprs la valeur de : rf ( rf )min

rf + 1
N N min
On tire du diagramme de GILLIAND la valeur de : qui nous permet de
N +1
calculer N.
Page 38
IV.6.11. Dtermination des dimensions de la colonne
Les conditions de marche dune colonne tour sont : pression, temprature et taux de
reflux en fonction du nombre dtages.
Les dimensions dune colonne sont :
IV.6.11.a. Hauteur de la colonne

Cette hauteur est lie au nombre de plateaux que lon doit installer lespacement que lon
doit prvoir entre chacun d eux [10]. .
La hauteur de la colonne est :
HT = h1 + h2 + h3 + h4 + h5 + h6 (m) (54)
D
h1 = : la distance entre le 1er plateau de la zone de rectification et le sommet de la
2
colonne.
h2 = (NR - 1) h p : la hauteur de la zone de rectification.
h3 = 2.hp : la hauteur de la zone dalimentation.
h4 = (NE 1) h p : la hauteur de la zone dpuisement.
h5 : hauteur libre au fond de la colonne (1-1,5m).
h6 : hauteur de la jupe de la colonne (1,5m).
hp= (0, 3 - 0,8m) : la distance entre deux plateau successives.
IV.6.11.b. Nombre de plateaux

Le nombre dtages est fonction du taux de reflux, ces tages sont galement appels
plateaux thoriques . En effet, un plateau parfois reprsente un quilibre entre les phases
liquides et vapeurs qui peuvent stablir parfaitement [10].
En pratique, il nen est pas ainsi et on doit intervenir la notion defficacit de 100%, le
plateau rel une efficacit infrieur 100%, le rapport entre le nombre de plateaux thoriques
et le nombre de plateaux rels est lefficacit moyenne des plateaux dune colonne [10].
Page 39
N br de plateaux thoriques
E : efficacit moyenne (%) = =
N br de plateaux rel
On peut donc lier le nombre dtages au nombre de plateaux. Il faut un taux de reflux tel que
le nombre de plateaux ne soit pas grand pour que la balance conomique soit bonne. En effet,
les plateaux constituent un investissement et le taux de reflux une dpense dexploitation [10].
IV.6.11.c. Espacement entre les plateaux
Les critres sont simples, il faut que la construction soit possible et que lon puisse pntrer
entre les plateaux pour assurer lentretien, la hauteur de la tour sera gal au produit du nombre
de plateaux par lespacement, plus lespace libre en tte et au fond du tour [10].
IV.6.11.d. Diamtre de la colonne
Si la vitesse des vapeurs dans la tour est trop grande (diamtre trop petit) des gouttelettes de
liquide peuvent tre entranes dun plateaux lautre, en effet il y aura un by-pass partiel des
plateaux, lquilibre nest pas atteint, ce qui amne une baisse fficacit de plateaux.
Le calcul se fait en prenant le dbit de vapeur de tte (distillat + reflux interne) car cest la
zone la plus change [10].
.D 2
Q = V. S = V. (55)
4
4.Q 4.Q
D2 = D= (56)
V . V .
O :
V : La vitesse du fluide (m/s).
S : Section du passage (m2).
Q : Dbit (m3/s).
Page 40
CHAPITRE V
PARTIE CALCUL
Prsentation de ltude
Prsentation de ltude
Notre travail a t ralis dans la perspective doptimisation des paramtres de la colonne de
stabilisation du brut au niveau de lusine du traitement du brut dUTBS, qui se compose de 3 trains
de production, pour rcuprer le maximum de brut stabilis.
Lors de la conception du 1er train, les calculs ont t faits sur la base des concentrations des puits
au moment de lexploitation. Ces puits ont subit un appauvrissement en hydrocarbures liquides.
Avec les dbits dexploitation proposs, la colonne PX0-CB-21-01 sengorgent car par la suite
on a constat que les nouveau puits fors sont riches en hydrocarbures liquides.
-But et intrt du travail
Lors des calculs, on doit tenir compte des limites des changements des valeurs des paramtres
dexploitation. Et ne pas proposer doprer des changements sur les installations dj existantes.
Les objectifs de ce calcul sont :
-Premire Partie
- Dtermination des dbits de la colonne de stabilisation du brut avec les concentrations
relles.
- De vrifier que la colonne de stabilisation du brut dj existante peut traiter (stabilis) la
charge, et davoir des produits finis rpondants aux exigences commerciales.
-Deuxime Partie
- Optimiser les paramtres de fonctionnement (pression, tempratureetc) de la colonne de
stabilisation du brut partir du changement des dbits dalimentation (charge chaude et
froide).
Page 41
CHAPITRE V Partie calcul
V.1. Donnes de dpart

La composition et les caractristiques de la charge de la colonne PX0-CB-21-01:
Tableau [V-1] : Composition de la charge chaude
dbit
Constituants Mi X'Li Mi.X'Li XiL=Mi.X'Li/Mn molaire
(Kmol/h)
H2 O 18 0,0128 0,2304 0,002048036 41,74208
N2 28 0,0001 0,0028 0.000024889 0,32611
CO2 44 0,0011 0,0484 0,00043023 3,58721
CH4 16 0,0031 0,0496 0,000440897 10,10941
C2H6 30 0,0265 0,795 0,007066792 86,41915
C3H8 44 0,0713 3,1372 0,027886718 232,51643
i-C4H10 58 0,0202 1,1716 0,010414407 65,87422
n-C4H10 58 0,0587 3,4046 0,030263649 191,42657
i-C5H12 72 0,0341 2,4552 0,021824388 111,20351
n-C5H12 72 0,064 4,608 0,040960728 208,7104
n-C6H14 86 0,0883 7,5938 0,067501644 287,95513
n-C7H16 100 0,0976 9,76 0,086757098 318,28336
n-C8H18 114 0,0697 7,9458 0,070630589 227,29867
n-C9H20 128 0,0724 9,2672 0,082376576 236,10364
n-C10H22 142 0,0663 9,4146 0,083686821 216,21093
n-C11H24 156 0,0523 8,1588 0,072523956 170,55553
n-C12H26 170 0,1824 31,008 0,275631567 594,82464
n-C12H26 170 0,0791 13,447 0,119531014 257,95301
somme 1506 1,00000 112,498 1 3261,1
La temprature dalimentation est de 114C.

La masse molculaire moyenne de la charge chaude est :
Mn = 112,498Kg/Kmole
Page 42
Tableau [V-2] : Composition de la charge froide

dbit
Constituants Mi X'Fi Mi.X'Fi XiF=Mi.X'Fi/Mn molaire
(Kmol/h)
H2 O 18 0,0128 0,2304 0,002048036 10,43584
N2 28 0,0001 0,0028 0.000024889 0,08153
CO2 44 0,0011 0,0484 0,00043023 0,89683
CH4 16 0,0031 0,0496 0,000440897 2,52743
C2H6 30 0,0265 0,795 0,007066792 21,60545
C3H8 44 0,0713 3,1372 0,027886718 58,13089
i-C4H10 58 0,0202 1,1716 0,010414407 16,46906
n-C4H10 58 0,0587 3,4046 0,030263649 47,85811
i-C5H12 72 0,0341 2,4552 0,021824388 27,80173
n-C5H12 72 0,064 4,608 0,040960728 52,1792
n-C6H14 86 0,0883 7,5938 0,067501644 71,99099
n-C7H16 100 0,0976 9,76 0,086757098 79,57328
n-C8H18 114 0,0697 7,9458 0,070630589 56,82641
n-C9H20 128 0,0724 9,2672 0,082376576 59,02772
n-C10H22 142 0,0663 9,4146 0,083686821 54,05439
n-C11H24 156 0,0523 8,1588 0,072523956 42,64019
n-C12H26 170 0,1824 31,008 0,275631567 148,71072
n-C12H26 170 0,0791 13,447 0,119531014 64,49023
somme 1506 1,00000 112,498 1 815.3
La temprature dalimentation est de 68.2C

La masse molculaire moyenne de la charge froide est :
Mn = 112,498Kg/Kmole
Les caractristiques de la charge de la colonne PX0-CB-21-01
Paramtres la charge froide la charge chaude

Temprature (C) 68,2 114
Pression (bar) 6 6,9
Dbit molaire (kmole/h) 815,3 3261,1
Page 43
V.2. Les dbits de distillat et du rsidu de la colonne PX0-CB-21-01

- Le constituant cl lger est n-C4H10.
- Le constituant cl lourd est i-C5H12.
les concentrations admises de i-C5H12, dans le distillat et de nC4H10 dans le rsidu sont:
-La teneur en cl lger dans le rsidu R est XRnC4=0,0433.
-La teneur en cl lourd dans le distillat D est YDiC5 =0,0322.
On peut estimer que les H2O ,N2,CO2,C2, C3 et i-C4 sont pratiquement absents dans le rsidu
ainsi que les n-C5, n-C6 et n-C7.,n-C12 dans le distillat.
On suppose donc que :
X R H2O = X R N2 = XR CO2 = XR C1=XRC2= XRC 3= XR iC4= 0
YD nC5= YD nC6= YD nC7 = .= YD n-C12 = 0
- Le bilan matire global de la colonne est donn par:
F+L=R+D.(1)
F : Dbit molaire de la charge froide (1ere charge).
L : Dbit molaire de la charge chaude (2eme charge).
R : Dbit molaire du rsidu.
D : Dbit molaire du distillat.
Le bilan de matire partiel (pour llment i) scrit :
F.XFi+L.XLi=R.XRi+D.YDi..(2)
O :
XLi : Concentration molaire du constituant i dans la charge chaude.
XFi : Concentration molaire du constituant i dans la charge froide.
XRi : Concentration molaire du constituant i dans le rsidu.
YDi : Concentration molaire du constituant i dans le distillat.
On applique lquation prcdente pour chaque constituant de la charge :
F. XF H2O +L . XL H2O = R. X R H2O +D. YD H2O...... (3)
F. XF N2 +L . XL N2 = R. X R N2+D. YD N2.......(4)
F. XF CO2 +L . XL CO2 = R. XR CO2+D. YD CO2.... (5)
F. XF C1 +L . XL C1 = R. XR C1+D. YD C1..... (6)
F. XF C2+L . XL C2 = R. XR C2+D. YD C2 ....... (7)
F. X C3 +L . XL C3 = R. XR C3+ D. YD C3.... (8)
F. X i-C4 +L . XL i-C4 = R. X i-C4+ D. YD i-C4....... (9)
F. X n-C4 +L . XL n-C4 = R. X n-C4+ D. YD n-C4... (10)
Page 44
F. X i-C5 +L . XL i-C5 = R. XR i-C5+ D. YD i-C5... (11)

F. X n-C5 +L . XL n-C5 = R. XR n-C5+ D. YD n-C5..... (12)
F. XF n-C6+L . XL n-C6 = R. XR n-C6+ D. YD n-C6.... (13)
F. XF n-C7+L . XL n-C7 = R. XR n-C7+ D. YD n-C7.... (14)
F. XF n-C8 +L . XL n-C8 = R. XR n-C8+ D. YD n-C8.... (15)
F. XF n-C9 +L . XL n-C9= R. XR n-C9+ D. YD n-C9......... (16)
F. XF n-C10+L . XL n-C10= R. XR n-C10+ D. YD n-C10.... (17)
F. XF n-C12+L . XL n-C12 = R. XR n-C12+ D. YD n-C12... (19)
En substituant les valeurs connues dans toutes les quations, obtient :
52.17792 = R. 0 +D. YD H2O..... (1)
0.40764 = R. 0+D. YD N2......(2)
4.48404 = R. 0+D. YD CO2.... (3)
12.63684 = R. 0+D. YD C1.... (4)
108.0246 = R. 0+D. YD C2 ....... (5)
290.64732 = R. 0+ D. YD C3. (6)
82.34328 = R. 0+ D. YD i-C4......... (7)
239.28468 = R. 0.0433+ D. YD n-C4. (8)
139.00524 = R. XR i-C5+ D. 0.0322.. (9)
260.8896 = R. XR n-C5+ D. 0... (10)
359.94612 = R. XR n-C6+ D. 0..... (11)
397.85664 = R. XR n-C7+ D. 0..... (12)
284.12508 = R. XR n-C8+ D. 0..... (13)
295.13136= R. XR n-C9+ D. 0.......... (14)
270.26532= R. XR n-C10+ D. 0.... (15)
213.19572= R. XR n-C11+ D. 0.... (16)
743.53536 = R. XR n-C12+ D. 0... (17)
322.44324= R. XR n-C12+ D. 0.... (18)
En additionnant les quations (1), (2), (3), (4), (5), (6), (7) et (8), on obtient :
790.00632=D. (YD H2O +. +YD nC4) + R . 0,0433.
790.00632=D. (YD H2O +..+YD nC4) + (L+F-D). 0,0433.
790.00632=D. (YD H2O +..+YD nC4-0,0433) + (815,3+3261,1). 0,0433.
Mais, on a :
Page 45
YD H2O + YD N2+ YD CO2+YDiC4+YD nC4 =1- 0.0322 = 0,9678.

790.00632=D.(0,9678-0,0433)+176.50812
D = 663.6 Kmol/h
R = F+L-D = 3261.1+815.3-663.6 = 3412.8 Kmol/h
R = 3412.8 Kmol/h
On remplace les valeurs de D et R dans les quations (1,2,3,4,5,,18) on

obtient :
52.17792 = 663.6 YD H2O..... (1)
0.40764 =663.6. YD N2......(2)
4.48404 = 663.6. YD CO2.. (3)
12.63684 = 663.6. YD C1.. (4)
108.0246 = 663.6. YD C2 ..... (5)
290.64732 = 663.6. YD C3.... (6)
82.34328 = 663.6. YD i-C4......... (7)
239.28468 = 3412.8. 0.0433+ 663.6. YD n-C4... (8)
139.00524 = 3412.8. XR i-C5+ 663.6. 0.0322.... (9)
260.8896 = 3412.8. XR n-C5 ... (10)
359.94612 = 3412.8. XR n-C6 . (11)
397.85664 = 3412.8. XR n-C7 ..... (12)
284.12508 = 3412.8. XR n-C8 ..... (13)
295.13136= 3412.8. XR n-C9 .......... (14)
270.26532= 3412.8. XR n-C10 ..... (15)
213.19572= 3412.8. XR n-C11 ..... (16)
743.53536 = 3412.8. XR n-C12 .... (17)
322.44324= 3412.8. XR n-C12 ..... (18)
A laide des quations prcdentes, on calcule la composition du distillat et du rsidu, les
rsultats figurent sur le tableau qui suit :
Page 46
Tableau [V-3] : Composition et dbits molaire et massique du distillat

dbit dbit
Constituants Mi Yi Yi.Mi Yi=Mi*Yi'/MmD molaire massique
(Kmol/h) (Kg/h)
Di=D.Yi Dim=Di.Mi
H2 O 18 0,07862857 1,41531426 0,03238481 52,1779191 939,202543
N2 28 0,000614 0,01720012 0,00039357 0,40764284 11,4139996
CO2 44 0,00675714 0,29731416 0,00680306 4,4840381 197,297677
CH4 16 0,01904286 0,30468576 0,00697173 12,6368419 202,18947
C2H6 30 0,16278571 4,8835713 0,11174446 108,024597 3240,73791
C3H8 44 0,43798571 19,2713712 0,44096192 290,647317 12788,482
i-C4H10 58 0,12408571 7,19697118 0,16467901 82,3432772 4775,91008
n-C4H10 58 0,1379 7,9982 0,18301249 91,51044 5307,60552
i-C5H12 72 0,0322 2,3184 0,05304896 21,36792 1538,49024
n-C5H12 72 0 0 0 0 0
n-C6H14 86 0 0 0 0 0
n-C7H16 100 0 0 0 0 0
n-C8H18 114 0 0 0 0 0
n-C9H20 128 0 0 0 0 0
n-C10H22 142 0 0 0 0 0
n-C11H24 156 0 0 0 0 0
n-C12H26 170 0 0 0 0 0
n-C12H26 170 0 0 0 0 0
somme 1506 0,99999999 43,703028 1 663,599993 29001,3294
La masse molculaire moyenne du distillat est :

MmD =43,703028Kg/Kmole
Page 47
Tableau [V-4] : Composition et dbits molaire et massique du rsidu

dbit dbit
Constituants Mi Xi Xi.Mi Xi=Mi.Xi/MmR molaire massique
(Kmol/h) (Kg/h)
Ri=R.Xi Rim=Ri.Mi
H2 O 18 0 0 0 0 0
N2 28 0 0 0 0 0
CO2 44 0 0 0 0 0
CH4 16 0 0 0 0 0
C2H6 30 0 0 0 0 0
C3H8 44 0 0 0 0 0
i-C4H10 58 0 0 0 0 0
n-C4H10 58 0,0433 2,5114 0,01995156 147,77424 8570,90592
i-C5H12 72 0,03447 2,48184 0,01971673 117,639216 8470,02355
n-C5H12 72 0,07644444 5,50399968 0,04372597 260,889585 18784,0501
n-C6H14 86 0,10546944 9,07037184 0,07205866 359,946105 30955,365
n-C7H16 100 0,11657778 11,657778 0,09261404 397,856648 39785,6648
n-C8H18 114 0,08325278 9,49081692 0,07539884 284,125088 32390,26
n-C9H20 128 0,08647778 11,0691558 0,08793779 295,131368 37776,8151
n-C10H22 142 0,07919167 11,2452171 0,0893365 270,265331 38377,6771
n-C11H24 156 0,06246944 9,74523264 0,07742002 213,195705 33258,53
n-C12H26 170 0,21786667 37,0373339 0,29423937 743,535371 126401,013
n-C12H26 184 0,09448056 16,0616952 0,12760052 322,443255 54815,3534
somme 1520 1 125,874841 1 3412,80191 429585,658
La masse molculaire moyenne du distillat est :

MmR =125,874841Kg/Kmole
V.3. Pression dans les diffrentes zones de la colonne

V.3.1. Pression de la charge froide (reflux)
la Pression de la charge froide (reflux) est gale :
Pr = 6 bars (donne)
V.3.2. la pression au sommet de la colonne
PS = Pr
Page 48
Ou :
PS : pression au sommet de la colonne
Donc la pression du sommet de la colonne gale :
Ps = 6 bars
V.3.3. la pression au fond de la colonne
La perte de charge entre le fond et le sommet est gale :
Pf = Ps + P1
O :
Pf : pression au fond de la colonne.
P1 : pertes de charges dues aux rsistances des plateaux.
P1 = 0, 2 bar
Pf = 6 + 0.2 =6.2
Pf = 6.2 bars
V.3.4. la pression de zone dalimentation chaude
Pac = Ps + P2.np1
Pac : La pression de zone dalimentation chaude.
P2 : La perte de charge au niveau de chaque plateau.
np1 : Nombre de plateau entre le sommet et lalimentation chaude.
P2 = (Pf Ps)/np1
P2 = (6.2 6)/21
P2 = 0,00952 bar
Pac = Ps + P2.np1
Pac = 6 + 0.00952. 9
Pac = 6.0857 bars
V.4. Dtermination de la temprature

V.4.1. Pour lalimentation chaude et froide
La temprature de lalimentation chaude est gale : Tac =114C
La temprature de lalimentation froide est gale : Taf = 68.2C
Ces tempratures sont dtermines par les approximations successives Annexe 1 :
Page 49
Tableau [V-5] : Coefficients d'quilibre des constituants Dalimentation
Pc =6.9 bar et Tac = 114C Pr =6 bar et Taf = 68.2C

Constituants Alimentation chaude Alimentation froide
KiL X'Li KiF X'Fi
H2 O 19 0,0128 4,2 0,0128
N2 80 0,0001 96 0,0001
CO2 21 0,0011 26 0,0011
CH4 52 0,0031 46 0,0031
C2H6 14 0,0265 10 0,0265
C3H8 5,5 0,0713 3,5 0,0713
i-C4H10 2,8 0,0202 1,7 0,0202
n-C4H10 2,3 0,0587 1,4 0,0587
i-C5H12 1,3 0,0341 0,68 0,0341
n-C5H12 1,3 0,064 0,575 0,064
n-C6H14 0,6 0,0883 0,25 0,0883
n-C7H16 0,29 0,0976 0,1 0,0976
n-C8H18 0,15 0,0697 0,046 0,0697
n-C9H20 0,08 0,0724 0,018 0,0724
n-C10H22 0,043 0,0663 0,007 0,0663
n-C11H24 0,024 0,0523 0 0,0523
n-C12H26 0,015 0,1824 0 0,1824
n-C12H26 0,015 0,0791 0 0,0791
somme 1,00000 1,00000
V.4.2. Dans le sommet de la colonne

A une pression P = 6 bar, pour la vrification de la temprature au sommet de la colonne
on utilise lquation disotherme de la phase vapeur par approximations successives :

1

Si linconnue est la temprature (P connue), on fait le calcul suivant :
- Calcul du point de rose dune vapeur :
Le mlange est entirement vapeur yi = Z
Page 50
On crit que la premire goutte de condensat liquide a une composition molaire Xi, telle
que:
Xi = 1

1

On obtient galement une quation en Ki dont la rsolution s'effectue par approximations

successives sur la pression ou la temprature, comme dans le cas prcdent, et fournit
simultanment la composition de la phase liquide.
- Pour une temprature T = 70C on a Y'Di / Ki = 1,138 1
Pour une temprature T = 76.6C Les rsultats de calcul sont illustrs dans le tableau suivant:
Tableau [V-6] : Coefficients d'quilibre des constituants du distillat
PD = 6 bar et Ts = 76.6C
Constituants Y'Di KiD Y'Di / KiD

H2 O 0,07862857 0,127 0,6191226
N2 0,000614 2210 2,7796*10-7
CO2 0,00675714 924,505 7,3089*10-6
CH4 0,01904286 44 0,00043279
C2H6 0,16278571 9.5 0,01713534
C3H8 0,43798571 3.3 0,13272294
i-C4H10 0,12408571 1.65 0,07520346
n-C4H10 0,1379 1.3 0,10607692
i-C5H12 0,0322 0.65 0,04953846
somme 0,99999999 1,0002401
V.4.3. Dans le fond de la colonne

En utilisant lisotherme de la phase liquide XiR.Ki = 1, pour une pression de fond de 6.2bar.
Linconnue est la temprature (P connue), on fait le calcul suivant :
- Calcul du point de bulle dun liquide :
Le mlange est entirement liquide X i = Z i
On crit que la premire bulle de vapeur mise a une composition molaire yi, telle
que:
Yi = 1
Page 51
Ki . xi = 1
La composition du liquide Xi tant connue, cette quation en K i ne dpend que des

deux paramtres pression et temprature, de sorte que si l'on se fixe l'un, l'autre est dtermin.
Mais, comme les coefficients Ki ne sont pas des fonctions explicites de P et T, la
rsolution de cette quation ncessite un calcul par approximations successives qui fournit la
pression ou la temprature de bulle, ainsi que la composition de la phase vapeur.
- Pour une temprature T = 150 CO on a XiR*Ki = 1,35622
-Pour une temprature T = 160 CO Les rsultats de calcul sont illustrs dans le tableau
suivant:
Tableau [V-7] : Coefficients d'quilibre des constituants du residu
PR = 6.2 bar et Tf = 160C
Constituants XRi KiR X'Ri.KiR
n-C4H10 0,0433 5,81 0,25157467
i-C5H12 0,03447 3,15 0,10884349
n-C5H12 0,07644444 3,03 0,2317279
n-C6H14 0,10546944 1,64 0,17317738
n-C7H16 0,11657778 0,94 0,11043283
n-C8H18 0,08325278 0,54 0,04521561
n-C9H20 0,08647778 0,34 0,029491
n-C10H22 0,07919167 0,20 0,01600371
n-C11H24 0,06246944 0,13 0,00828472
n-C12H26 0,21786667 0,08 0,01761143
n-C12H26 0,09448056 0,08 0,00763741
somme 1,00000056 1,00000016
V.5. Volatilits relatives
i = K i K .. (21)
r
Ki : Coefficient dquilibre de llment (i).

Kr : Coefficient dquilibre de llment le plus lourd.
On admet le constituant cl lourd iC5 H12 comme lment de rfrence.
Page 52
Tableau [V-8] : les volatilits relatives de la charge de la colonne
Alimentation chaude Alimentation froide

Constituants KiL Lc,i= KiL/KiC5 KiF Lf ,i= KiF/KiC5
H2 O 19 14,6153846 4,2 6,17647059
N2 80 61,5384615 96 141,176471
CO2 21 16,1538462 26 38,2352941
CH4 52 40 46 67,6470588
C2H6 14 10,7692308 10 14,7058824
C3H8 5,5 4,23076923 3,5 5,14705882
i-C4H10 2,8 2,15384615 1,7 2,5
n-C4H10 2,3 1,76923077 1,4 2,05882353
i-C5H12 1,3 1 0,68 1
n-C5H12 1,3 1 0,575 0,84558824
n-C6H14 0,6 0,46153846 0,25 0,36764706
n-C7H16 0,29 0,22307692 0,1 0,14705882
n-C8H18 0,15 0,11538462 0,046 0,06764706
n-C9H20 0,08 0,06153846 0,018 0,02647059
n-C10H22 0,043 0,03307692 0,007 0,01029412
n-C11H24 0,024 0,01846154 0 0
n-C12H26 0,015 0,01153846 0 0
n-C12H26 0,015 0,01153846 0 0
-Dtermination de la volatilit moyenne

Les valeurs moyennes des volatilits relatives sont dtermines par les formules suivantes :
Li = 0.5( Lc .i + Lf .i ) (22)
Pour la zone de rectification :
i .m = 0.5( Li + Di ) (23)
Pour la zone dpuisement :
i .m = 0 .5( Li + Ri ) .(24)
Les rsultats sont prsents dans les tableaux suivants :
Page 53
Tableau [V-9] : Volatilits moyennes
Alimentation Rectification Epuisement

Constituants
Li KiD Di iDm KiR Ri iRm
H2 O 10,3959276 0,127 0,19538462 5,29565611 - 0 5,1979638
N2 101,357466 2210 3400 1750,67873 - 0 50,678733
CO2 27,1945701 924,50 1422,31538 724,754977 - 0 13,5972851
CH4 53,8235294 44 67,6923077 60,7579186 - 0 26,9117647
C2H6 12,7375566 9.5 14,6153846 13,6764706 - 0 6,36877828
C3H8 4,68891403 3.3 5,07692308 4,88291855 - 0 2,34445701
i-C4H10 2,32692308 1.65 2,53846154 2,43269231 - 0 1,16346154
n-C4H10 1,91402715 1.3 2 1,95701357 5,81 1,84445666 1,8792419
i-C5H12 1 0.65 1 1 3,15 1,0024221 1,00121105
n-C5H12 0,92279412 - 0 0,46139706 3,03 0,96232521 0,94255967
n-C6H14 0,41459276 - 0 0,20729638 1,64 0,52125949 0,46792612
n-C7H16 0,18506787 - 0 0,09253394 0,94 0,30072663 0,24289725
n-C8H18 0,09151584 - 0 0,04575792 0,54 0,1724166 0,13196622
n-C9H20 0,04400452 - 0 0,02200226 0,34 0,10826159 0,07613306
n-C10H22 0,02168552 - 0 0,01084276 0,20 0,06415501 0,04292027
n-C11H24 0,00923077 - 0 0,00461538 0,13 0,04210173 0,02566625
n-C12H26 0,00576923 - 0 0,00288462 0,08 0,02566217 0,0157157
n-C12H26 0,00576923 - 0 0,00288462 0,08 0,02566217 0,0157157
V.6. Taux de vaporisation e0 et composition des phases vapeur et liquide

dalimentation
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 :
n n
X'
i =1
0 .i = 1 (25) et Y'
i =1
0 .i =1 .. (26)
Ou : X0.i et Y0.i sont calculs de la faon suivante :
X ' Li
X '0.i =
1 + e0 ' ( K i 1) ....(27)
Page 54
Y ' 0 .i
X ' 0 .i = Y '0.i = X ' 0.i .K i (28)
Ki
Le taux de vaporisation de lalimentation et la composition de ses deux phases lquilibre
est dtermin laide de lquation de TREGOUBOV :
n n
X ' L.i

i =1
X ' 0.i =
i =1 1 + e0 ' (Ki 1)
= 1 (29)
Avec : X0.i : Concentration molaire du constituant (i) dans la phase liquide.

Y0.i : Concentration molaire du constituant (i) dans la phase vapeur.
XL.i : Concentration molaire du constituant (i) dans la charge.
Pour la charge chaude.
-Taux de vaporisation e0 et composition des phases vapeur et liquide de la charge
Lapplication de la formule pour lalimentation chaude :
le taux de vaporisation est gale :

e'0 =0,09
Page 55
Tableau [V-10] : Fractions des phases liquide et vapeur de l'alimentation chaude
P = 6.9 bars T= 114C e'0 = 0.09
Constituents XL.i KiL X0.i Y0.i

H2 O 0,0128 19 0,0048855 0,09282443
N2 0,0001 80 1,233*10-5 0,00098644
CO2 0,0011 21 0,00039286 0,00825
CH4 0,0031 52 0,00055456 0,02883721
C2H6 0,0265 14 0,01221198 0,17096774
C3H8 0,0713 5,5 0,05074733 0,27911032
i-C4H10 0,0202 2,8 0,01738382 0,0486747
n-C4H10 0,0587 2,3 0,05255148 0,1208684
i-C5H12 0,0341 1,3 0,03320351 0,04316456
n-C5H12 0,064 1,3 0,06231743 0,08101266
n-C6H14 0,0883 0,6 0,09159751 0,05495851
n-C7H16 0,0976 0,29 0,10426237 0,03023609
n-C8H18 0,0697 0,15 0,07547374 0,01132106
n-C9H20 0,0724 0,08 0,07893589 0,00631487
n-C10H22 0,0663 0,043 0,07254861 0,00311959
n-C11H24 0,0523 0,024 0,05733643 0,00137607
n-C12H26 0,1824 0,015 0,20014265 0,00300214
n-C12H26 0,0791 0,015 0,08679432 0,00130191
somme 1 1,00135231 0,98632669
Pour la charge chaude, le taux de vaporisation est gal :

e'0 =0,09
n
L.i X 'L.i
e'0 = (30)
i=1 L.i
: Constante qui varie entre les volatilits relatives des constituants cls lourds et le plus
volatiles.
1 < < 1,91402715 Le calcul pour = 1.345 nous donne les rsultats dans le tableau
suivant, sachant que :
Page 56
Tableau [V-11] : Verification du taux de vaporisation
= 1.345
Constituants X'Li Li X'Li.Li Li- e'0
H2 O 0,0128 10,3959276 0,13306787 10,4267897 0,01470213
N2 0,0001 101,357466 0,01013575 116,018888 0,00010134
CO2 0,0011 27,1945701 0,02991403 23,3525081 0,00115723
CH4 0,0031 53,8235294 0,16685294 55,7338301 0,00317945
C2H6 0,0265 12,7375566 0,33754525 12,1022207 0,02962858
C3H8 0,0713 4,68891403 0,33431957 3,4062437 0,09997852
i-C4H10 0,0202 2,32692308 0,04700385 1,11170346 0,04786917
n-C4H10 0,0587 1,91402715 0,11235339 0,60308277 0,19744821
i-C5H12 0,0341 1 0,0341 -0,3064 -0,09884058
n-C5H12 0,064 0,92279412 0,05905882 -0,44289425 -0,13988157
n-C6H14 0,0883 0,41459276 0,03660854 -0,92594023 -0,03934679
n-C7H16 0,0976 0,18506787 0,01806262 -1,13326782 -0,01557214
n-C8H18 0,0697 0,09151584 0,00637865 -1,22292299 -0,00508874
n-C9H20 0,0724 0,04400452 0,00318593 -1,26688851 -0,00244884
n-C10H22 0,0663 0,02168552 0,00143775 -1,28819598 -0,00108648
n-C11H24 0,0523 0,00923077 0,00048277 -1,29848333 -0,00036142
n-C12H26 0,1824 0,00576923 0,00105231 -1,3014 -0,00078576
n-C12H26 0,0791 0,00576923 0,00045635 -1,3014 -0,00034075
somme 1 0,09031158
V.7. Taux de rebouillage

P = J / R..........................................................................................................(31)
R = R + J(32)
P : Taux de rebouillage.
R : dbit molaire de rsidu total.
R : dbit molaire de rsidu.
J : Dbit de laspiration de la pompe.
R = 3412.8 + 3066.4 = 6479.2
R = 6479.2 Kmol/h
Page 57
P = 3066.4 / 6479.2 = 0.4732

P = 0.4732
V.8. Bilan thermique de la colonne

Ltablissement du bilan thermique de la colonne est ncessaire pour dterminer les
charges thermiques.
QF + QL + QRb = QD + QR .... .(33)
QL = L.hL ..........(34)
hL = L .H0,L+ (1-L).h0,L.......(35)
QF = L .hF ...... (36)
QR = R.hR......... (37)
QD = D.HD.. (38)
Donc:
L.hL+ F.hF + QRb= R.hR+ D.HD
Avec :
QL : charge thermique de lalimentation chaude Kcal/h
QF : charge thermique de lalimentation froide Kcal/h
QR: charge thermique de rsidu Kcal/h
QD: charge thermique de distillat Kcal/h
QRb charge thermique fourne par le rebouilleur Kcal/h
hL: enthalpie de la charge chaude Kcal/Kmol
H0,L: Enthalpie molaire de la phase vapeur de la charge chaude Kcal/Kmol
h0,L: Enthalpie molaire de la phase liquide de la charge chaude Kcal/Kmol
hF : enthalpie de la charge froide Kcal/Kmol
e0 : taux de vaporisation de la charge
hR : enthalpie de rsidu en tat liquide Kcal/Kmol
hD: enthalpie de distillat Kcal/Kmol
h0D : enthalpie de distillat en tat liquide Kcal/Kmol
H0D : enthalpie de distillat en tat vapeur Kcal/Kmol
L, R, D, F : sont les dbits molaires respectivement de la charge chaude et froide, du
rsidu et de distillat.
Page 58
V.8.1. Dtermination des enthalpies des phases vapeurs et liquides

Les enthalpies sont dtermines par le monogramme de lenthalpie en fonction de la
pression, la temprature et la masse molculaire (voire annexe 2 et 3).
-Enthalpie de la charge Chaude hc
Donnes :
PAc = 6.9 bar; (Pa = 100,076039 psi)
TAc = 114C ; (Ta= 237,20F)
Mn = 112.498 Kg/Kmol
Tableau [V-12] : Calcul Mmoy de la charge Chaude a ltat liquide et vapeur

P = 6.9 bars T = 114C
e'0 = 0,09
Constituants X0.i Mi X0.i . Mi Y0.i Mi Y0.i. Mi
H2 O 0,0048855 18 0,087939 0,09282443 18 1,67083974

N2 1,233*105 28 0,00034524 0,00098644 28 0,02762032
CO2 0,00039286 48 0,01728584 0,00825 48 0,363
CH4 0,00055456 16 0,00887296 0,02883721 16 0,46139536
C2H6 0,01221198 30 0,3663594 0,17096774 30 5,1290322
C3H8 0,05074733 44 2,23288252 0,27911032 44 12,2808541
i-C4H10 0,01738382 58 1,00826156 0,0486747 58 2,8231326
n-C4H10 0,05255148 58 3,04798584 0,1208684 58 7,0103672
i-C5H12 0,03320351 72 2,39065272 0,04316456 72 3,10784832
n-C5H12 0,06231743 72 4,48685496 0,08101266 72 5,83291152
n-C6H14 0,09159751 86 7,87738586 0,05495851 86 4,72643186
n-C7H16 0,10426237 100 10,426237 0,03023609 100 3,023609
n-C8H18 0,07547374 114 8,60400636 0,01132106 114 1,29060084
n-C9H20 0,07893589 128 10,1037939 0,00631487 128 0,80830336
n-C10H22 0,07254861 142 10,3019026 0,00311959 142 0,44298178
n-C11H24 0,05733643 156 8,94448308 0,00137607 156 0,21466692
n-C12H26 0,20014265 170 34,0242505 0,00300214 170 0,5103638
n-C12H26 0,08679432 170 14,7550344 0,00130191 170 0,2213247
somme 1,00135231 1506 118,684534 0,98632669 1506 49,9452836
Page 59
-Enthalpie de la charge Chaude a ltat liquide

Donnes :
PAc = 6,9 bar ; (Pa = 100,076039 psi)
TAc = 114C ; (Ta = 237,20 F)
Mmoy = 118,684534 Kg/Kmol.
Rsultats :
h0L = 4615,33 Kcal/Kmole.
-Enthalpie de la charge Chaude a ltat vapeur
Donnes :
PAc = 6,9 ba r; (Pa = 100,076039 psi)
TAc = 114C ; (Ta = 237,20 F)
Mmoy = 49,9452836 Kg/Kmol.
Rsultats :
H0L = 5868,783 kcal/Kmole.
-Enthalpie de la charge Chaude
hL = eo . HoL + (1 - eo).hoL
hL = 4728,14077 Kcal/Kmole.
-Enthalpie de la charge froide hF
Donnes :
Pa = 6 bars ; (Pa = 87,0226 psi)
Ta = 68,2C ; (Ta=154,76F)
Mmoy = 112,498 Kg/Kmol.
(eo = 0) liquide totale
hF = 1374,975 Kcal/Kmole
Rsultats :
hF = 1374,975 Kcal/Kmol
-Enthalpie du distillat
Page 60
Tableau [V-13] : Calcul Mmoy du distillat a ltat liquide et vapeur
E = 0,001
Constituants X0.i Mi X0.i . Mi Y0.i Mi Y0.i. Mi
H2 O
0,07869727 18 1,41655091 0,61966356 18 11,1539442
N2
0,00019143 28 0,00535996 8,6619*10-8 28 2,425*10-6
CO2
0,00351293 48 0,15456896 3,7998*10-6 48 0,00016719
CH4
0,01825778 16 0,29212441 0,00041495 16 0,00663919
C2H6
0,16141369 30 4,84241081 0,01699092 30 0,50972745
C3H8
0,43698065 44 19,2271488 0,13241838 44 5,82640873
i-C4H10
0,12400511 58 7,19229619 0,07515461 58 4,35896739
n-C4H10
0,13785864 58 7,99580126 0,10604511 58 6,15061635
i-C5H12
0,03221127 72 2,31921172 0,04955581 72 3,56801804
Somme 0,99312878 43,445473 1,00024722 31,5744909
Enthalpie du distillat ltat liquide h0D

Donnes :
Pb = 6 bars ; (Pb= 87,0226 psi)
Tb = 76,6C ; (Tb= 169,88F)
MD = 43,445473 Kg/Kmol
Rsultats :
h0D = 1206,8 Kcal/ Kmol.
-Enthalpie du distillat ltat vapeur H0D
Donnes :
Pb = 6 bars ; (Pb = 87,0226 psi)
Tb = 76,6C ; (Tb = 169,88F)
MD = 31,5744909 Kg/Kmol
Rsultats :
H0D = 3507,79 Kcal/Kmole.
-Enthalpie du distillat
hD = eo . HoD + (1 - eo).hoD
hD = 1101,696 Kcal/Kmole
-Enthalpie de rsidu hR :
Donnes :
Page 61
Pf = 6,2 bars ; (Pf = 89,923 psi)

Tf =160C ; (Tf = 320F)
MR = 125,874841 Kg/Kmol
Rsultats
hR = 8601,459 Kcal/Kmole.
V.8.2. Calcul des charges thermiques
-Charge thermique de lalimentation chaude QL
QL = L.hL
QL = 1,5419 . 107 Kcal/h
-Charge thermique de lalimentation froide QF
QF = F .hF
QF = 1,121017. 106 Kcal/h
-Charge thermique de distillat QD
QD = D.hD
QD =7,31085. 105 Kcal/h
-Charge thermique de rsidu QR
QR = R.hR
QR = 2,9355059.107 Kcal/Kmole
- Charge thermique fourne par le rebouilleur QRb
QRb = QR + QD - (QL+ QF)
QRb = 2,9355059.107 + 7,31085. 105 (1,5419. 107+1,121017. 106)
QRb = 1,35461. 107 Kcal/h
NB:
1 BTU/lb = 0.5564 Kcal/Kg
1(F) = T (C) 9/5 + 32
1 (psi) = 14.5 bars
V.9. Calcul du diamtre de la colonne
On calcule le diamtre de la colonne l'aide de l'quation suivante :
4.Vi
D' = .................................(39)
.W
Vi : dbit des vapeurs qui quittent le sommet de la colonne en m/s.
W : la vitesse des vapeurs en m/s.
-Calcul de la vitesse admissible des vapeurs :
Page 62
Pour la dtermination de W on utilise la formule suivante :
L V
W = 8,47.10 5.C. .......... ........( 40 )
V
C : coefficient qui dpend du type de plateau, de lespacement entre les plateaux et de la
charge sur les plateaux.
L et v : ce sont respectivement les poids spcifique de la vapeur et du liquide dans les
conditions de fonctionnement.
Calcul de v :
M m .Vm
V = .... (41)
Vi
Mm : la masse molculaire moyenne de la vapeur qui quitte la colonne en Kg/Kmol.
Vm : dbit molaire de la vapeur qui quitte la colonne en Kmol/h.
Vi : dbit volumique de la vapeur qui quitte la colonne en m/h.
Vm = D = 663.6 Kmol/h.
Mm = 43.70302 Kg/Kmol.
273 + TS P0
Vi = 22,4.Vm . .Z . (42)
273 PS
Avec :
Ts et Ps : respectivement la temprature et la pression au sommet de la colonne.
Z : facteur de compressibilit.
Vi : dbit volumique des vapeurs qui quittent le sommet de la colonne.
- calcul du facteur de compressibilit (Z) :
Ces facteurs sont dtermins daprs les paramtres P et T critiques
Connaissant ces derniers, on calcul les paramtres rduits daprs les formules :
Tr m = TS / TCm au : TCm = Yi . TCi
Pr m = PS / PCm au : PCm = Yi . PCi
Le calcul dtaill est donn dans le tableau :
Page 63
Tableau [V-14] : Reprsentant la temprature et la pression critiques moyennes du

Distillat
Constituants Y'Di Tci(K) Pci(bar) Tci.Y'Di Pci.Y'Di
H2 O 0,07862857 647,299 221,187 50,8961947 17,3916175
N2 0,000614 126,194 33,9305 0,07751971 0,02084317
CO2 0,00675714 304,1 73,6868 2,05484627 0,49791202
CH4 0,01904286 190,699 46,3936 3,63145436 0,88346683
C2H6 0,16278571 305,428 48,8253 49,7193138 7,94806113
C3H8 0,43798571 369,898 42,5534 162,010038 18,6377811
i-C4H10 0,12408571 408,096 36,463 50,6388819 4,52453724
n-C4H10 0,1379 425,199 37,953 58,6349421 5,2337187
i-C5H12 0,0322 460,398 33,3227 14,8248156 1,07299094
somme 1,00000 392,488007 56,2109287
TCm = 392,488007K
PCm = 56,2109287 bars
Les paramtres rduits sont :
Trm = (76.6 + 273.15) / 392,488007 = 0,99971 0K
Prm = (6 / 56.2109287) = 0.10674bars
On calcul les pressions et les tempratures rduites afin de dterminer le facteur de
compressibilit Z daprs le diagramme Annexe 4
On trouve : Z = 0.94
On remplace Z dans lquation (43)
273.15 + 76.6 1
Vi = 22, 4 663.6 .0,94 = 2981.86 m 3 / h
273.15 6
43.70328 663.6
V = = 9.726.Kg / m 3 .
2981 .86
V = 9.726 kg/m3
Calcul de L
Page 64
M m
L = '
. (43)
M i .X
Li
Di
Li : Masse volumique du constituant (i) aux conditions opratoires (au sommet de la

colonne).
Le calcul
Li se fait partir du diagramme annexe 5 :
Pour les constituants C3H8, i-C4H10 et n-C4H10 on applique la formule suivante :

t
dt = t = 1000.d t .....................( 44)
1000
Pour les constituants CH4 et C2H6 on applique le suivant :
t = 1000 .(d 420 a.(t 20) ).......... .......... ...( 45)

a = 0,001828 0,00132 .d 420 .......... .......... ..( 46)
avec d 420 d 415
Les rsultats sont reprsents dans le tableau suivant :
Tableau [V-15] : Calcul de la masse volumique du constituant au sommet de la
colonne
'
M i . X Di
Constituant M i di Li Y Di' '
X Di '
M i . X Di
Li
H2 O 18 0.998 998 0,07862857 0,6191226 11,1442068 0,01116654
N2 28 0.80637 806,37 0,000614 2,779*10-7 7,782*10-6 9,651*10-9
CO2 46 0.825335 825,335 0,00675714 7,308*10-6 0,00032159 3,896*10-7
CH4 16 0.24954 249,54 0,01904286 0,00043279 0,00692468 2,775*10-5
C2H6 30 0.29847 298,47 0,16278571 0,01713534 0,51406014 0,00172232
C3H6 44 0.39 390 0,43798571 0,13272294 5,83980947 0,01497387
i-C4H10 58 0.48 480 0,12408571 0,07520346 4,36180072 0,00908708
n-C4H10 58 0.502 502 0,1379 0,10607692 6,15246154 0,0122559
i-C5H12 72 0.55 550 0,0322 0,04953846 3,56676923 0,00648503
Somme 1,00000 1,0002401 31,5863619 0,0557189
Page 65
31,58636
L = = 567 , 07 Kg / m 3
0,055718
On prend C = 770 pour Hp = 600(mm) pour les plateaux clapets
5 567 , 07 9 , 726
W = 8 , 47 . 10 770 = 0 , 493 m / s
9 , 726
Donc :
4.Vi 4 2989,86
D' = = = 1,465m
.W 3,14 0,493 3600
V.10. Calcul de la hauteur de la colonne

La hauteur totale de la colonne est calcule par la formule suivante :
HT = H1 + H2 + H3 + H4 + H5+H6. . (47)
H1 : la distance entre le sommet est le premier plateau de la zone de rectification.
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 distance entre le premier plateau de la zone d'puisement et le fond de la colonne.
Puisque la hauteur de la colonne dpend de la distance entre les plateaux mais la
distance entre les plateaux dpend aussi du diamtre de la colonne.
Hauteur de lespace libre au dessus du plateau suprieur :
H1 = D/ 2 = 1,465 / 2 = 0,7325 m.
Hauteur de la zone de rectification :
H 2 = H p (N R 1 ) . (48)
NR : Nombre de plateaux dans la zone de rectification, NR = 8.

Hp : Distance entre deux plateaux, Hp = 0. 61 m.
H 2 = 0 . 61 ( 8 1 ) = 4 . 27 m
Hauteur de la zone dexpansion :
H 3 = 6 H p = 3 . 61 m
Hauteur de la zone dpuisement :

H 4 = H p (N E 1 ) . (49)
NE : Nombre de plateaux de la zone dpuisement, NE = 13.

H 4 = 0 . 61 (13 1 ) = 7 . 32 m
Hauteur de lespace libre au-dessous du plateau infrieur :
Page 66
H 5 = 9 . 335 m
H6 :lespace entre le plateau 13 et le plateau 14 au niveau de rtrcissement ;

H 6 = 1 . 56 m
- La hauteur totale sera :

H T = H 1 + H 2 + H 3 + H 4 + H 5 + H 6 (50)
H T = 26 . 8275 m
La hauteur relle est de :
HT = 24.435 m
Tableau [V-16] : Tableau rcapitulatif et de comparaison entre les dimensions
relles et calcules
calcule relle
Nombre total des
21 21
plateaux
Diamtre (m) 1 ,465 1,5
Hauteur (m) 24,435 26,8275
Conclusion
Daprs les calculs que nous avons fait sur les paramtres de fonctionnement de la colonne
de stabilisation, on a constat que les rsultats des calculs simuls manuellement sont
acceptables avec les valeurs des paramtres actuels.
Page 67
CHAPITRE VI
PARTIE OPTIMISATION ET
SIMULATION
CHAPITR VI Partie optimisation et simulation
VI.1. Dfinition
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 de lavoir travers
lexprience et lintuition seule, en particulier pour les systmes complexes avec plusieurs
variables interactives.
Le logiciel HYSYS est lun des logiciels de simulation des procds de gnie des
procds. Il a t dvelopp par la socit canadienne HYPROTECH
Il a t conu pour permettre le traitement dune vaste gamme de problmes allant des
sparations bi et tri - phasiques simples, de la compression la distillation et la
transformation chimique.
Les diffrentes taches quun simulateur de procds devrait effectuer sont :
Le calcul des flashs
La rsolution des bilans de matire et dnergie.
Loptimisation du procd.
Le dimensionnement des quipements.
VI.2. Le 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 les
proprits de transports.
La russite de la simulation dpend du choix du modle thermodynamique, parce que ce
dernier est tabli pour une classe de fluide et un domaine des conditions P et T
recommandes.
Lquation de Peng-Robinson est recommande pour le systme tudi (traitement de
gaz), car elle a t vrifi pour le cas actuel et a donn des rsultats avec des erreurs
presque ngligeables.
Cette quation est donne par :
RT a
P =
V b V (V + b ) + b (v b ) [1]
Page 68
Z 3
(1 B )Z 2
+ (A 2B 3B 2
)Z (4 B 2 B 3 B )= 0
N
Rt C i
b =
x i bi b i = 0 . 0778
PC
i 1 ; i
N N
a = X I X j (a i a j )0 . 5 (1 K ij )
i 1 j 1
a i = a ci i
a ci = 0 . 457235
(RT ci )2
Pci
i0 .5 = 1 + m i (1 Tri0 .5 )
2
m i = 0 . 37646 + 1 . 54226 i 0 . 2699 i
bP
B =
A,B : constantes lies la temprature et la pression avec et RT
aP
A =
(RT )2
a, b : constantes lies la temprature et la pression critiques.
Z : coefficient de compressibilit.
ni : nombre de moles du constituant i.
R : constante des gaz parfaits.
Tci : temprature critique du constituant i.
Pci : pression critique du constituant i.
Kij : coefficient dinteraction binaire.
xi, xj : Fractions molaires des constituants i et j.
Tri : temprature rduite.
VI.3. Remarques concernant lutilisation des logiciels de simulation
Lemploi des logiciels de simulation offre de nombreuses possibilits pour dfinir les
spcifications dterminant le fonctionnement des colonnes distiller ; il est ncessaire de
fixer autant de valeurs ou spcifications quil existe de produits ; leurs choix doit tre
raisonner de manire faciliter les convergences.
Page 69
Lorsque, pour un nouveau problme, lutilisateur ne possde pas les ordres de

grandeurs ; il lui est conseill de raliser un premier calcul en dfinissant des spcifications
garantissant une solution. Lanalyse de cette premire solution oriente les volutions et les
dcisions prendre.
Fig. VI. 1: SCHEMA TECHNOLOGIQUE DE LA SIMULATION
Page 70
VI.4. Le choix des paramtres

Lors de la simulation des colonnes de stabilisations, le choix des paramtres
opratoires (pression, temprature, taux de reflux) vise :
Diminuer la quantit de chaleur fournir au rebouilleur.
Diminution du taux de reflux, or une augmentation de ce dernier va entraner une
augmentation des quantits de chaleur dans le rebouilleur ainsi que la puissance des
pompes et par consquent la consommation lectrique.
Avoir une colonne avec un faible taux dengorgement et une perte de charge
minimale.
Mais il est important de noter quil faut assurer une bonne sparation avant tout et obtenir
des produits rpondant aux spcifications.
Dans notre cas les paramtres suivants conditionnent le bon fonctionnement du colonne de
stabilisation :
Le taux de reflux.
La temprature de sommet et du fond.
Q de rebouilleur.
Do une tude de linfluence de la variation de ces paramtres sur les performances de la
colonne de stabilisation est mene.
Pendant la simulation de cette colonne deux considrations sont prisent en compte :
La quantit de nC5 dans le produit de tte doit tre minime (poids molculaire).
La rcupration du produit brut stabilis doit tre maximale.
Tout en assurant :
Une quantit de chaleur minime fournir au rebouilleur.
Une colonne avec un faible taux de lengorgement et de perte de charge.
Page 71
Fig. VI. 2 : Schma synoptique du train 10
VI.5. Prsentation des rsultats et interprtation

VI.5.1. Influence de taux reflux froid sur le fonctionnement de la colonne de
stabilisation
Dans cette partie, on tudie l influence du taux de reflux de tte sur les paramtres de
fonctionnement de la colonne (TVR, poids molculaire, duty four, la temprature,
lefficacit de prchauffeur).
Page 72
Tableau [VI -1]: Influence de taux reflux froide sur le fonctionnement du colonne de
stabilisation
Taux de Poid Duty four.107 T0de T0de U.A.105

TVR (psia)
reflux (%) molculaire (Kj/h) sommet(C0) fond(C0) (Kj/C0.h)
10 6.0 44.90 5.12515 83.46 161.9 9.862
15 6.3 44.22 5.26365 79.11 161.5 9.005
18 6.8 43.57 5.35492 77.87 161.2 8.55
20 7.007 42.85 5.41881 77.21 161 8.258
22 7.2 42.25 5.4819 76.62 160.8 7.977
25 7.5 41.90 5.57667 75.91 160.3 7.579
27 7.7 41.50 5.63869 75.59 160 7.326
30 7.9 39.94 5.73373 75.25 159.6 6.963
40 8.0 37.74 6.0441 74.77 157.8 5.818
Fig. VI. 3: Influence du taux de reflux en fonction la TVR
Page 73
Fig. VI. 4 : Influence du taux de reflux en fonction du poids molculaire
Fig. VI. 5 : Influence du taux de reflux en fonction de duty four
Page 74
Fig. VI. 6 : Influence du taux de reflux en fonction de la temprature de sommet
Fig. VI. 7 : Influence du taux de reflux en fonction de la temprature de fond de

colonne
Page 75
Fig. VI. 8 : Influence du taux de reflux en fonction de U.A Prchauffeur du stabiliseur
On note daprs ce tableau et ces profils que :

A Chaque fois on augmente le taux de reflux nous remarquons la diminution de
temprature de fond et temprature de tte.
La temprature en sortie du four de rebouillage est ajuste de telle sorte que la TVR en
sortie de colonne soit de 7 psi.
La quantit de chaleur absorbe dans le rebouilleur qui correspond aux taux de reflux 20
%.
pour un ratio charge froide charge chaude de 20%-80% on obtient un couple Duty four-
UxA Prchauffeur du stabiliseur optimum qui permet davoir une taille du Four et du
Prchauffeur raisonnable.
La masse molculaire de vapeur de tte lintervalle (39- 44).
-On remarque que la bonne sparation sera faite aux taux de reflux de 20 %.
VI.5.2. Prsentation et interprtation des profils de la colonne de stabilisation

Les profils les plus importants prsenter qui sont donns par les figures ci-dessous sont :
Le profil de temprature.
Le profil de pression
Le profil des dbits molaire du liquide et de la vapeur.
Le profil dengorgement.
Page 76
Le profil de la perte de charge.

Profil de temprature
La figure montre le profil de temprature au niveau de chaque plateau le long de la
colonne de stabilisation.
Fig. VI. 9 : Profil de temprature au nivaux de chaque plateaux

La figure ci-dessus montre que 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 13eme), ce pic sexplique par
la diffrence entre la temprature dalimentation et celle du plateau.
Profil de pression
Fig. VI. 10: Profil de pression au nivaux de chaque plateaux
Page 77
La figure ci-dessus montre que le profil de pression actuel prsente une volution
normale le long de la colonne (diminution de la pression du 1er plateau en tte, jusquau
dernier plateau du fond).
Profils des dbits molaires du liquide et de la vapeur
Fig. VI. 11: Profils des dbits molaires du liquide et de la vapeur au nivaux de chaque
plateaux
Les graphes des dbits liquides et vapeur refltent la prsence de deux sections :
Une premire section qui stend du 13me jusquau 21me plateau, prsente une volution
des profils des dbits liquides et vapeurs uniforme et parallle qui est du labsence dune
source de perturbation.
Au niveau du 13me plateau (plateau dalimentation) on constate une augmentation
importante du dbit liquide accompagne dune diminution du dbit vapeur; ceci est d la
condensation du gaz pendant son mouvement ascendant vers le haut de la colonne.
Une deuxime section qui stend du 1ere jusquau 13me plateau, ce stade on remarque
que la stabilit est relativement rtablie.
Page 78
VI.7. Tableaux comparatifs

Aprs loptimisation, on a obtenu les rsultats ci-dessous :
Tableau [VI -2]:les paramtres de fonctionnement de la colonne P10-CB-21-01 :
Les valeurs relles Les rsultats de hysys
Temprature au sommet(C) 76.6 77.21
Temprature au fond(C) 160 161
7
Duty four.10 (Kj/h) 5.662 5.41881
TVR (psia) 7.1 7.007
Taux de reflux (%) 22 20
Poid molculaire 43.25 42.85
Conclusion :
Pour les rsulta smill par HYSYS, la diminution du pourcentage de la charge froide (22%
- 20%) on a trouv que les valeurs des paramtres de fonctionnement de la colonne sont
acceptable condition quil ny pas une perturbation dans le fonctionnement.
On a montrait que l nergie consomme par le four (duty de four) elle est optimise, et
avec une bonne confissions de transfert de chaleur au niveau de lchangeur thermique.
Page 79
CONCLUTION GENERALE
CONCLUSION GENERALE
Au court de notre sjours UTBS, nous avons mis le point sur plusieurs problmes
techniques, parmi lesquels le problme de sparassions des lments [C1,C2,C3,]au niveau de
la section de fractionnement.
Dans notre tude, nous avons essay doptimiser les paramtres de fonctionnement de la
section de stabilisation (taux de reflux), ceci dont le but dobtenir une colonne avec un faible
taux dengorgement et une quantit de chaleur minimale fournir au rebouilleur et avec une
bonne production de brut en qualit et en quantit.
Daprs le calcul et loptimisation on peut dire que : pour un ratio charge froide charge
chaud de 20%-80% on obtient :
Un couple Duty four-UxA Prchauffeur du stabiliseur optimum.
Une colonne avec un faible taux dengorgement et de perte de charge.
Une bonne stabilisation de brut (TVR=7)
La masse molculaire de vapeur de tte lintervalle [39 44].
- on peut conclure que :
la bonne sparation sera faite aux taux de reflux de 20 %.
La nature de la charge actuelle et les spcifications requises des produits finis

dfinissent le choix des paramtres de fonctionnement optimal des quipements.
Loptimisation de reflux froid au sommet de la colonne de stabilisation (diminution du

dbit du reflux froid) permet daugmenter la temprature au fond et la temprature du
sommet, tout en conservant le gradient de temprature et la qualit des produits finis.
Page 82
RFRENCES BIBLIOGRAPHIQUES
[1] KHELLADI, M., Quest-ce que le ptrole . Elaboration et conception

(SONATRACH), (1984).
[2] www.Wikipdia.com , l'encyclopdie libre Documentation Internet
[3] Durousset, M., Le march du ptrole. Paris : Collection HISTEGE, (1999).
[4] WAUQUIER, J, P., Le raffinage du ptrole, produits ptroliers, schmas
de fabrication.
[5] HADJARAB, D., Cours de raffinage 4me anne
[6] UTBS Manuelle opratoire
[7] CILILE (J-C)., Technique de lingnieur Distillation Absorption .1.
Gnralits sur les colonnes de fractionnement [J2 621].Vol. J2 (1994).
[8] Alexandre, R., Bernard, D., Le gaz naturel production, traitement et
transport1994.
[9] WUITHIER, J.P., Raffinage et gnie chimique. Paris : Edition Technip,
publication de linstitut franais, (1972). Tome 1.
[10] WUITHIER, J.P., Raffinage et gnie chimique. Paris : Edition Technip,
publication de linstitut franais, (1972). Tome 2.
[11] HYSYS simulator 6.0 canada 2006
[12] HOCINE, M., SAADI, O., Etude de performances des PDM
utiliss pour le forage des drains horizontaux de Hassi Massaoude.
Mmoire dingnieur, Universit Mhamed Bouguerra, Boumerdes,
(2006).
ANNEXE 1
Diagramme des coefficients dquilibre des hydrocarbures (Scheibel et Jenny) [9]
ANNEXE 2
Diagramme des enthalpies total des paraffines hydrocarbures liquides en fonction de
la temprature et la masse molaire [10]
ANNEXE 3
Diagramme des enthalpies total des paraffines hydrocarbures vapeurs en fonction de
la temprature et la masse molaire [10]
ANNEXE 4
Facteur de compressibilit (Z) on foction TC et PC [8]
ANNEXE 5
Variation de la densit des hydrocarbures ou fraction paraffiniques en ltat liquide en
fonction de la temprature [9]
Rsum
La composition de la charge change avec le temps et devient plus lgre par rapport
au dpart; ce qui perturbe le fonctionnement de lunit de traitement du ptrole brut de
lUTBS rgion de Hassi Messaoud qui pour mission de produire le brut stabilis et gaz
lifts.
Daprs les calculs que nous avons fait sur les paramtres de fonctionnement de la
colonne de stabilisation, on a constat que les rsultats des calculs simuls manuellement
sont acceptables par rapport aux valeurs des paramtres actuels.
Notre tude consiste optimiser de nouveaux paramtres de fonctionnement de la

colonne de stabilisation P10-GA-21-01 afin damliorer son rendement.
Abstract
The composition of the load changes in time and becomes lighter with regard to the
departure; what perturbs the functioning of the unity of treatment of the crude oil of the
UTBS region of Hassi messaoud which has for mission to produce the stabilized crude oil
and the gas lifts.
According to the calculations that we made on the parameters of functioning of the

column of stabilization, we noticed that the results of the calculations feigned manually are
acceptable with regard to the values of the current parameters.
Our study consists in optimizing new parameters of functioning of the column of

stabilization P10-GA-21-01 to improve its efficiency.

" !

# %&* ) (' +,- . / .1
."+86 ,* -2 * - -&& 6- .2 3
4! 6* -
<. =. >* @

?! - -A &8 . +,B- CD3 E =*!- 2/! -3
.+E .=. +,B-* G! =*! - 2/ HD- '+
.K' *+2/ =* ' - -A &8 . +,B-* G! I / E JE !-3

Calcul de Vérification Et Optimisation Des.......

Transféré par

Droits d'auteur :

Formats disponibles

Calcul de Vérification Et Optimisation Des.......

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

Calcul de Vérification Et Optimisation Des.......

Transféré par

Droits d'auteur :

Formats disponibles

REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L'ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

FACULTE DES SCIENCES ET DE LA TECHNOLOGIE

Mmoire de Fin d'tude en Vue de l'Obtention du Diplme de Master

Option : Gnie chimique

Calcul de vrification et optimisation des

RAHIM Khaled et BEN RAHLA Sohaib

Prsident: Mme ADJEL Fatima

Promotion Juin 2012

Nous tenons aussi remercier profondment tout le

BEN RAHLA SOHAIB.

Je ddie ce modeste travail :

Ma trs chre mre,

Mon trs cher pre,

Mes trs chers frres et soeurs,

A tous les enseignants et toutes les enseignantes de

Et tous mes amis.

Jai le grand honneur de ddier ce modeste travail :

Ma trs chre mre Mon trs cher pre

A mes trs chrs frres et mes surs

Sans oublie rendre un grand hommage

Mes cousins dans nos mmoires et toute la famille

CHAPITRE I : Historique de champ de Hassi Messaoud

I.2. Situation du champ de Hassi Messaoud.02

I.2.1. Situation gographique.02

I.2.2. Situation gologique.03

I.3. Zones et numrotation des puits.04

I.4. Intrt ptrolier...05

CHAPITRE II : Gnralits sur le ptrole brut

II.2. Dfinition du ptrole brut.07

II.3. Origine du ptrole brut..07

II.4. Classification des ptroles brut.07

II.4.1. Classification industrielle08

II.4.2. Classification chimique...08

II.4.3. Classification conomique (technologique)....08

II.4.3.a. Daprs la teneur en soufre..........08

II.4.3.b. Selon la teneur potentielle en fractions claires PF=350c... .08

II.5. Composition des ptroles brut..09

CHAPITRE III : Description de lunit

III.2. Vue gnrale de LUTBS10

III.3. Description de lunit..13

III.3.1. Units de Production..13

III.3.3. Traitement et Stabilisation dHuile14

3-3.a) Sparation dHuile ..14

III.3.4. Stockage, Comptage et Expdition de lHuile Stabilise..19

III.3.5. Traitement et du stockage de lhuile hors spec..19

III.3.6. Unit de Compression de Gaz....20

III.3.7. Lunit de traitement des eaux huileuses..21

III.3.8. Les systmes utilits...21

CHAPITRE IV : Thorie de rectification

IV.3. Description dune colonne de rectification.....24

IV.4. Principe de fonctionnement dune colonne de rectification....25

IV.5. Types des colonnes de rectification....26

IV.5.1. Colonnes plateaux..26

IV.5.2. Colonnes garnissage....27

IV.6. Lois fondamentales de la rectification...27

IV.6.1. Lois de DALTON..27

IV.6.2. Lois de RAOULT et dHENRY27

IV.6.3. Bilan matire de la colonne et de ses zones...28