Laboratory Catalysis

Spectrochemistry
Laboratoire
Catalyse&&
Spectrochimie

Introduction la Distillation
2CAD2: CM = 10

Jean-Pierre GILSON
Tl.: [Link].15 email: [Link]@[Link]

ENSICAEN, Universit de Caen, CNRS

6 bd du Marchal Juin, 14050 Caen, France
http:/[Link]

Introduction
Procds de Sparation
But du Cours

E. Sorel, La Rectification de lAlcool, 1893

Industrie
(Batch + Continu)
Taille: : ~15 m x H: 50 m
Dbit: ~160 000 b/j, i.e. 25 400 m3/j

Laboratoire (Batch)

JPG, 2CAD2, 2012-13

Plan du Cours
Procds de Sparation
I.
II.

Introduction: Ncessit, types de sparations, Labo vs. Industrie,

Aspects de la distillation: Bilan matire, nergie, phnomnes de transport

quilibres Liquide - Vapeur

I.
II.

Thermodynamique
Simulations des quilibres

Distillation Flash & Continue

I.
II.
III.
IV.

Distillation Flash
Distillation Continue
Phnomnes de Transport
quipement

Conception & Utilisation de Colonnes

I.
II.
III.
IV.
V.

Bilans Matires
Temprature de Tte
Temprature de Fond
Temprature dAlimentation
Taux de Reflux & Nombre de Plateaux (MacCabe-Thiele)

Applications en Raffinage Ptrolier

I.
II.
3

Distillation atmosphrique
Distillation sous vide
JPG, 2CAD2, 2012-13

Bibliographie
1.

E. Sorel, La Rectification de lAlcool, (1894), Gauthier-Villars et fils, Paris

2.

J-P. Wauquier, Le Raffinage du Ptrole 2. Procds de Sparation, (1998), ditions Technip, Paris

3.

J. Vidal, Thermodynamique: Application au Gnie Chimique & lIndustrie Ptrolire, (1997), ditions
Technip, Paris

4.

D. Morvan, Les Oprations Unitaires, (2009), Ellipses Marketing, Paris

5.

J. Gmehling, B. Kolbe, M. Kleiber & J. Rarey, Chemical Thermodynamics for Process Simulation, (2012),
Wiley-VCH, Weinheim

6.

H. Kister, Distillation Design, (1992), McGraw-Hill, New York

7.

P.C. Wankat, Separation Process Engineering, 3rd ed., (2011), Prentice-Hall, Upper Saddle River

8.

J. D. Seader & E. J. Henley, Separation Process Principles, (1998), J. Wiley & Sons, Chichester

9.

M. Chabanel & B. Illien, Thermodynamique Chimique, (2011), Ellipses, Paris

10. J.M. Prausnitz, R.N. Lichtenthaler & E. Gomez de Azevedo, Molecular Thermodynamics of Fluid-Phase
Equilibria, 3rd ed., (1999), Prentice-Hall, Upper Saddle River
11. W. McCabe, J. Smith & [Link], Unit Operations of Chemical Engineering, 7th ed., (2004), McGraw-Hill,
New York
12. R.H. Perry & D.W. Green, Perry's Chemical Engineers' Handbook, 8th ed., (2008), McGraw-Hill, New York
4

JPG, 2CAD2, 2012-13

Logiciel Base de Donnes

1. ChemSep ([Link]): Gratuit en version lite. Effectue toutes les
oprations dcrites dans ce cours

2. Dortmund Data Bank ([Link]):Base de donnes thermodynamiques et

logiciel accessible gratuitement (versions explorer et acadmique lite)

JPG, 2CAD2, 2012-13

Procds de Sparation
I. Introduction
Sparations: Pourquoi ?
a. Fractionner: nenvoyer vers des units industrielles que les coupes ncessaires
b. Recycler: les procds de conversion (catalytique) nont pas toujours une conversion
maximum et la slectivit dsire; on recycle ce qui nest pas aux spcifications
c. Purifier: les oprations en aval ont des exigences (contaminants, polluants) sur
leurs charges

b.

a. & c.
6

JPG, 2CAD2, 2012-13

Procds de Sparation
I. Introduction
Sparations: Comment ?

L = Liquide, V = Vapeur, G = Gaz, S = Solide

7

JPG, 2CAD2, 2012-13

Procds de Sparation
I. Introduction
Sparations: Utilisation & Maturit

JPG, 2CAD2, 2012-13

Procds de Sparation
I. Introduction
Ptrole Brut (Fractionationner)

Procds (Recycler)

JPG, 2CAD2, 2012-13

Produits (Purifier)

Procds de Sparation
I. Introduction
Capacits de Distillation de Ptrole Brut
Raffinerie

10

Location

Barrels per Day*

Jamnagar (Reliance Industries

Limited)

Jamnagar, India

Paraguana (PDVSA)

Paraguan,
Falcon, Venezuela

940 000

Ulsan (SK Energy)

Ulsan, South Korea

850 000

GS Caltex Yeosu Refinery (GS

Caltex)

Yeosu, South Korea

730 000

ExxonMobil

Singapore

605 000

1 240 000

Port Arthur (Motiva Enterprises, i.e.

Port Arthur, Texas, USA
RD Shell & Saudi Aramco)

600 000

Baytown (Exxon Mobil)

Baytown, TX, USA

572 500

Ras Tanura (Saudi Aramco)

Saudi Arabia

550 000

S-Oil Ulsan (S-Oil)

Ulsan, South Korea

503 000

Marathon Oil (Marathon

Petroleum)

Garyville, LA

490 000

Pernis (RD Shell)

Pernis, Pays-Bas

416 000

Antwerp (Total)

Anvers, Belgique

360 000

Normandie (Total)

Gonfreville, France

205 000

JPG, 2CAD2, 2012-13

*: 1 barrel = 159 l; par 6,25 pour convertir en m3

Procds de Sparation
II. Aspects de la Distillations
Distillation: Paramtres ?
Sparation par EQUILIBRE: Thermodynamique & Transfert de Matire
Sparation par CHAUFFAGE: Thermodynamique & Transfert de Chaleur
Fonctionnement en CONTINU: Bilan & Transfert de Matire
Dimension INDUSTRIELLE: Intgration sur site, Contrle, conomie
Complexit !

11

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
1. Rappel de Formules:
Corps Pur: Soit la transformation T & P constantes dans un systme ferm

+dni

Evolution ( ) si dGTot < 0, soit i > i

Equilibre entre phases si dGTot = 0, soit i = i, i.e. galit des i de
chaque composant i dans toutes les phases.

-dni

dP H
=
dT T V0
P H 1 1
ln 0 =

P
R T0 T

Relation de Clapeyron:
Clausius-Clapeyron:

Mlanges:
Degrs de Libert

Variance : v = (
Nombre de variables - contraintes

1
c
Constituants

0
r

0
k

Equations chimiques indpendantes

Autres conditions imposes

2
=1

2 pour T & P

Phases

PTot = Pi et Pi = yi .PTot

Loi de Dalton (Gaz Parfaits):

12

2
n

) +

Loi de Raoult (Solution idale):

PA = [Link]*

Loi de Henry (Solution idale dilue):

PA = xA .K A

JPG, 1MA1C1, 2012-13

quilibres L - V
I. Thermodynamique
2.

Corps Pur (ex. CO2)

Volume Molaire (intensif)

13

JPG, 2CAD2, 2012-13

quilibres L - V
2.

14

I. Thermodynamique

Corps Pur (ex. H2O)

JPG, 2CAD2, 2012-13

quilibres L - V
I. Thermodynamique
2.

Corps Purs: H2O

2
3

Courbe de tension de vapeur

PointTriple PointCritique

Chaleur Latente
Chaleur Sensible

Dpenses & Pertes nergtiques

Distiller est Energivore !
15

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
3.

Mlanges Binaires:
xi yi: sauf dans le cas dun Azotrope (+ ou -) avec:
xi fraction molaire de i en phase liquide
P = PB0 + ( PA0 + PB0 ) .xA
PA0 .PB0
yi fraction molaire de i en phase vapeur
P= 0
0
0
P
+
P

P
(
). y A
A
B
A
zi fraction molaire de i (phase liquide + vapeur)

Bulle
Ligne Conodale

Rose

16

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
3.

Mlanges Binaires:
Diagramme T-H: les lignes isothermes courbes des T-H des constituants
Evolution de T en fonction de la Chaleur absorbe

Corps Pur : n-Butane

OA (A, A): le liquide chauffe (Chaleur Sensible)

AB (AB, AB): vaporisation (Chaleur Latente)
17

Mlange Binaire

Remplacement des paliers par des courbes

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Principe de la Distillation (Batch)

T1: T1 < TbA&B, x1 = x2

T2:
bullition
commence,
composition de la vapeur
(distillat), enrichie en B est y2

18

T3: T3 > TbB & T3 < TbA, x1 = x2

La vapeur (distillat) est enrichie en
B et les compositions x3 & y3 se
dplacent vers la droite p/ T2

T4: T4 = TbA, x1 = x2
La dernire goutte de liquide se
vaporise et sa composition, y4, est A
pur

A & B purs: une dtapes (Plateaux) est ncessaire. Plus facile de mlanger que de sparer !
La composition de la vapeur (distillat) change goutte goutte, de y2 y4
A la fin de la distillation, la composition du distillat est identique la solution de dpart
JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Rgle du Levier

Systme binaire de composition 1

A la Tconodale, la composition est:

liquide: 2
vapeur: 3

- Rgle du levier:

Liquide = a 0,25 ici

Vapeur b

19

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Principe de la Distillation Fractionne
Dpart:
xA = 0,8 & xB = 0,2

R = Rsidu
D = Distillat
20

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Principe de la Distillation Fractionne - Rsum: E = A+B+C+D

Distillation en 4 Plateaux
B concentr en D4: xB = 0,2 0,8
A concentr en R1: xa = 0,8 0,95
R2 & R3 intermdiaire, peu utiles
Sparation totale: Plateaux ! (Plus
facile de mlanger [S ] que de
sparer [S ])
Rendements + Succession dtapes
de Chauffe/Refroidissement !

21

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
En pratique au Laboratoire: Colonne de Fractionnement (Vigreux)

22

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
En pratique: Fractionnement (Rectification) Discontinu
DFINITIONS
Distillat
Vapeur condense enrichie en llment le plus
volatil
Rsidu (Fond)
Reste du bouilleur, appauvri en compos le plus
volatil
BILANS MATIRE
Global
F = W +D1 +D2 +D3 + Pertes
Composant Volatil
Fxf = Wxw + D1xd1 + D2 xd2 + Dxd3 +Pertes
VARIABLES
Taux de Reflux

23

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Cas des Mlanges Azotropiques:
On butte sur lAzotrope !
Formellement: de 2 Mlanges: A-Az + B-Az

24

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
AZEOTROPES: Composition Liquide (x) vs. Vapeur (y)
Diagramme de McCabe-Thiele

25

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Rsum: Idal Dviation - Azotropie

26

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Rsum: Azotropie

27

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
4.

Mlanges Ternaires:
Donnes dquilibre reportes sur un Triangle quilatral
Chaque point dans le triangle (R, S) y/c artes (L, P, Q) &
sommets (A, B, C) reprsente une composition du mlange
La des longueurs des abaisses dun point intrieur (eg. R)
sur les 3 cts est gale la hauteur du triangle = Cste
Pour chaque point: xA + xB + xC = Cste
LMNO-OPQ : Courbe (binodale) de Saturation ou Envelope de
Solubilit de ABC, i.e. dlimite la zone dimmiscibilit
R: lintrieur de la zone dimmiscibilit; se scinde en 2
phases de composition M & P (ligne dquilibre ou
conodale)
O: Point Critique ou Point de Plissement (Plait Point) de la
binodale

Calculs de VARIANCE !

28

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
5.

Mlanges Complexes: les charges & produits ptroliers contiennent une grande

varit dhydrocarbures (P, N, A, O) de poids molculaire trs variable (C1C100+).

Lapparition de charges renouvelables (biomasse) va encore complexifier la situation
(oxygns). Il est inconcevable de connaitre la nature chimique de chaque constituant.
On recueille des coupes entre Tx-Ty*

On peut recueillir a & b ~purs

PF = TRose

Tx
Ty

PI = TBulle

*: on ne peut plus recueillir a, b, c, d, purs

29

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Mthodes dvaluation des Mlanges complexes: Mthodes standardises (ASTM,
ISO, API) qui sont soit rapides et peu prcises ou lentes et (plus) prcises
Mthode Distillation Flash (FC*): sparation la plus simple, un plateau
(routine) - *:Flash Curve
Discontinu (Batch) ou Continu (F, z)
La charge est chauffe sous PZ, une TZ1
Dtente (isenthalpique)dans le Ballon de Flash (Flash Drum)
Vapeur enrichie en volatils & Liquide enrichi en lourds
On rpte une nouvelle TZ2 et sous PZ
On peut recommencer une P diffrente

30

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Mthode Distillation Flash (FC): On peut mettre en CASCADE (Srie) des Ballons
de Flash P (ou T) constante pour amliorer les sparations Consommation dEnergie !

P = Constante
31

Diagramme T, x dj rencontr
Meilleure sparation car plus de plateaux
Cot important: Chauffer & Refroidir
Retenons le principe !
JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Mthode True Boiling Point (TBP): Sparation plus longue, plus prcise

32

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Comparaison True Boiling Point (TBP): mlange 4 Constituants vs. Ptrole Brut

33

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Slectivit (S) de Sparation: Overlap & Gap des Mlanges complexes
a: courbe TBP dun mlange binaire

b: courbe ASTM dun mlange binaire

S=1

S=

34

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Slectivit (S) de Sparation: Overlap & Gap des Mlanges complexes
c: Courbe ASTM dun Brut Ptrolier
S = Overlap

35

d: Courbe ASTM dun Brut Ptrolier

S = Gap

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Comparaisons FC, ASTM & TBP

TBP: Sparation plus prononce (TTBP > TFC ) et donc de meilleure qualit

36

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

I. Thermodynamique
Distillation Simule: Analyse Chromatographique en Phase Gaz

talonnage par n-paraffines de Tb connues

La colonne GC spare par Tb
Mesurer laire entre 2 paraffines
Reporter sur un diagramme T = f(% Distill)
ASTM D 2887: fractions dont Tb < 538C
ASTM D 3710: fractions dont Tb < 250C

Remarque:
- Il nest pas encore possible danalyser compltement
(molculairement) les fractions ptrolires
- Mme si cela tait possible, la quantit disomres
obtenus ncessiterait un regroupement par familles
(lumping)

37

JPG, 2CAD2, 2012-13

quilibres Liquide - Vapeur

II. Simulation

38

JPG, 2CAD2, 2012-13

Chemistry Data Series, Dechema

Distillation Continue
I. Distillation Flash
1. Rsultat Dsir: en continu !
FD (Distillate)
x (Composition)
PDi (Pressure)
TDi (Temperature)
F (Feed)
z (Composition)
Pi (Pressure)
Ti (Temperature)

39

JPG, 2CAD2, 2012-13

Distillation Continue
I. Distillation Flash
2.

Distillations Flash en cascade rarranges:

40

JPG, 2CAD2, 2012-13

Distillation Continue
I. Distillation Flash

+ Volatil

3. Cascades de Flash P Constante

Botes Bleues Ballons de Flash

Liquide : doit tre rchauff
Vapeur : doit tre refroidie

41

L & V changent Chaleur

changeur de Chaleur

JPG, 2CAD2, 2012-13

- Volatil
L & V changent Matire
Exit changeur de Chaleur
Systme en mouvement:
Condenser une partie de V1 & Recycler
Vaporiser une partie de L5 & Recycler

Distillation Continue
I. Distillation Flash
4. Distillations Batch: Rectification & Stripping

Rectification
Seuls les composs lgers sont rcuprs
eg.: Production dalcools (Calvados)

42

Stripping
Seuls les composs lourds sont rcuprs
eg.: Enlever de (petites quantit de) solvant
Peu utilis industriellement (casser des
azotropes par ajout dun 3ime composant
eg.: MeOH/Tol + Et3N)
JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
1. Description: Nous venons de dcrire le Principe de la Distillation Continue !

TColonne

[Link]
43

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
1. Description: Instrumentation & Rgulation de la colonne

44

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
2. Bilans Matire & Thermique: Tour de Stabilisation d essence
Charge: Essence (Gasoline) ~ C5-10 - Essence Brute: contient encore du C3-4
Stabilisation: Distiller Essence Brute (14,34 t/h)
Distillat: C3-4 , Essence Sauvage (4,32 t/h)
Rsidu: C5-10 , Essence Stabilise (10,02 t/h) Base pour Carburant (Spcification RVP)

Bilan Thermique
-

Chaleur reue
Enthalpie dalimentation
Chaleur fournie par le rebouilleur
Enthalpie du reflux froid

Chaleur extraite
Enthalpie des vapeurs de tte (D + Reflux)
Enthalpie du Rsidu

45

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
Schma:

d & e: Corps cls de la Sparation

C
o
r
p
s
C
l

46

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
3. Rglage dune colonne: Tour de Stabilisation d essence
a)

But de la Stabilisation:
Ajuster la tension de vapeur de lessence stabilise la norme en vigueur (Reid
Vapor Pressure)
Veiller ce que lessence sauvage ne contienne pas de C5

b) Moyens disponibles:
Tour de Stabilisation permettant de ttaiter un certain Dbit horaire de charge
Pression peut tre ajuste
Temprature peut tre ajuste
Taux de Reflux peut tre ajust
c) Influence de P
PMIN : Le Distillat doit pouvoir tre condens totalement la T de leau de
refroidissement disponible dans le condenseur (Tension de vapeur de l essence
sauvage ~30C : 10 kg/cm2, i.e. 9,8 bar
Si la pression dopration est > PMIN, les TBULLE augmentent et donc la TTOUR
donc, peu de libert au niveau de P
47

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
3. Rglage dune colonne: Tour de Stabilisation d essence
d) Influence de T: - TTETE, T en tte de tour rgle la qualit du produit de tte, le Distillat (essence
sauvage). Ces vapeurs sont en quilibre avec le liquide du premier plateau, donc la TROSEE. En
agissant sur cette T, on agit sur la composition du Distillat

T : le Distillat salourdit (senrichit en produits moins volatils); son point sec

T
: le Distillat sallge (sappauvrit en produits plus volatils); son point sec
Rglage de la TTETE: en agissant sur le dbit de reflux froid qui, extrayant plus ou
moins de calories au systme, fait varier TTETE Qualit de lessence sauvage
- TFOND, T en fond de tour, rgle la qualit du produit de fond , le Rsidu (essence
stabilise). Ce liquide est en quilibre avec la vapeur du dernier plateau, donc la TBULLE. En
agissant sur cette T, on agit sur la composition du Rsidu

T : le Rsidu salourdit, sa teneur en butane diminue; sa tension de vapeur

T
: le Rsidu sallge, sa teneur en butane augmente; sa tension de vapeur
Rglage de la TFOND: en agissant sur le dbit de produit chauffant (eau en
gnral) dans le rebouilleur

48

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
3. Rglage dune colonne: Tour de Stabilisation d essence
e) Influence du Reflux: La chaleur fournie au systme provient du rebouilleur. Elle sert chauffer
(chaleur sensible), vaporiser (chaleur latente) le produit distiller et vaporiser le reflux froid inject en
tte de tour.

Reflux: extrait du systme une quantit de chaleur excdentaire p/ quantit

ncessaire au chauffage & vaporisation du distillat. Comme TTETE (= TROSEE) &
TFOND (= TBULLE) sont constantes, le dbit du reflux est variable variation du flux
de calories fourni par le rebouilleur
Reflux externe & interne: permet le transfert de matire (lavage des vapeurs) et
de chaleur. Cest la cl du Fractionnement.
Taux de Reflux: L/D permet de contrler la Slectivit et le nombre dtages de
la Distillation

Normand, Fig. 26 & 27, p. 149

49

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
3. Rglage dune colonne: Tour de Stabilisation d essence
e) Influence du Reflux:

L0: Reflux externe, dbit molaire

L: Reflux interne (Rectification)
L: Reflux interne (Strippage: L + LA)

+ riche en V

+ riche en L

50

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
Exemple industriel: Dbutaniseur

51

JPG, 2CAD2, 2012-13

Distillation Continue
II. Distillation Continue
Exemple industriel: Dbutaniseur

52

JPG, 2CAD2, 2012-13

Distillation Continue
III. Phnomnes de Transport
1.

Types de Phnomnes - Gnralisation

Matire - Loi de Fick (2ime) :
Chaleur - Loi de Fourier:
Impulsion - Loi de Newton:

C
= D 2C
t
T
= 2T
t

u
= 2u
t

C: Concentration
D: Diffusivit
T: Temprature
: Conductivit thermique
u: Vitesse linaire
: Viscosit
t: Temps

2
2
2
+
+
=

x 2 y 2 z 2 Laplacien

En commun: Transport dune quantit physique (Matire, Chaleur, Impulsion) travers un

gradient (Concentration, Temprature, Vitesse linaire) afin dliminer ce dernier.

2. Transport de Matire
2 mcanismes sont responsables du transport de matire:
Diffusion Molculaire: mouvement spontan scopique dans une phase donne (g, l, s); il
est assez lent. Voir ci-dessus
Diffusion Turbulente: mouvement macroscopique alatoire dun fluide
La Diffusion Molculaire entre 2 phases (g, l) joue un rle important en distillation. Elle se passe
linterface des 2 phases et conditionne les choix dquipements.
Les coefficients de Diffusion Molculaire [cm2/sec] sont de lordre de:
Gaz:
Di = 10-1
Liquides: Di = 10-5
Solides: Di = 10-10 - 10-13
53

JPG, 2CAD2, 2012-13

Distillation Continue
III. Phnomnes de Transport
3. Combinaison Transport de Matire & Chaleur
Phnomnes complexes (Turbulences): Fick & Newton ( linaires ) se compliquent
Introduction de Flux (J), coefficients de transfert (k, Matire & Chaleur) et corrlations
(nombres sans dimensions, e.g. Re), plus utiles Gnie Chimique
J

M ,q
Matire, Chaleur

k M ,q
Coefficient de Transfert

C ,T

.(Interface).(

Gradient de Concentration, Temprature

Qualitativement: en distillation, on va favoriser les changes, contacts (Matire Ji &

Chaleur JqL, JqV) dans les films (L & V) lInterface L/V

TCol

54

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
1. Rle
Effectuer le plus efficacement possible le transfert de matire et de chaleur entre les
phases liquide & vapeur; ces transferts se passent aux interfaces et les quipement
devront favoriser ltendue de celles-ci. En particulier, les quipements favoriseront:
temps de contact entre les phases
Qualitatif ! Mtier de lIngnieur:
mlange appropri
rendre ces termes QUANTITATIFS

2. Types de Colonnes
Colonnes Plateau (TRAY):
L & G scoulent contre-courant ou courant crois
Colonne peut travailler sous Pression
changes de Matire (& Chaleur) seffectuent sur des plateaux (Trays)
Contre-Courant

55

Courant Crois

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
Colonnes TRAY: Types de plateaux
Valve Plate

Sieve Plate

Ouverture variable
Fonction de la vitesse ascendante du Gaz
Se ferme si le flux gazeux est trop faible
Empche le liquide de descendre trop rapidement
vers ltage infrieur

Bubble (Bell) Cap Plate

Ouverture fixe
Turbulence accrue

56

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements

57

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
Colonnes TRAY: Comparaisons

58

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
Colonnes Remplies (Packed):

L & G scoulent contre-courant

Colonne peut travailler sous Pression
changes de Matire (& Chaleur) seffectuent
sur les structures de remplissage (Packings)
59

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
PACKED Columns: Comparaisons

60

JPG, 2CAD2, 2012-13

Distillation Continue
IV. Equipements
TRAYS vs. PACKED Columns: Comparaisons

TRAYS

PACKED
Moins cher pour de petites colonnes ( < 0,6 m)
Packings en cramique rsistent mieux aux
milieux corrosifs
Moussage plus controllable
Fentre dopration plus troite
Nettoyage difficile et coteux si contamination
par des particules solides
lments assez fragiles lors de stress en T ou P

Consulter, par exemple:

- SULZER: [Link]
- KOCH-GLITSCH: [Link]

61

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

Cahier des Charges
Pour:
Concevoir (Design) de nouvelles colonnes
Conduire (Operate) des units existantes
Rparer/Corriger (Troubleshooting) des dysfonctionnements
Nous avons besoin de:
Bilans Matires
Temprature de Tte de Colonne
Temprature de Fond de Colonne
Temprature dAlimentation
Taux de Reflux & Nombre de Plateaux Mthode de McCabe-Thiele
62

Pour:

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

I. Bilans Matires
Colonne: Rectification dEssence Sauvage (cf. Distillation Continue)
Charge: Essence Sauvage , i.e., Propane & Butane sous P = 20 bar

Propane: 35 vol %

Butane: 65 vol %
Spcifications de la Sparation:
Propane: < 1 vol % Butane
Butane: < 2 vol % Propane
Soient:
F: nombre de molcules dAlimentation (posons 100 molcules)
D: nombre de molcules de Distillat (Propane < 1 vol % Butane)
W: nombre de molcules de Rsidu (Butane < 2 vol % Propane)
zF: fraction molaire dun composant dan l Alimentation
xD: fraction molaire dun composant dan le Distillat
xW: fraction molaire dun composant dan le Rsidu
Bilan Global:
D+W=F
Bilan par compos:
[Link] + [Link] = [Link] = [Link]
o Propane:
D.0,99 + W.0,02 = 100 x 0,35
o Butane:
D.0,01 + W.0,99 = 100 x 0,65
ce qui donne: D = 34,02 & W = 65,98
63

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

I. Bilans Matires
Colonne: Rectification dEssence Sauvage
On dresse ainsi le tableau suivant:

Alimentation
zF

[Link]

[Link]

Rsidu
xD

[Link]

xW

Propane

35

0,35

33,68

0,99

1,32

0,02

Butane

65

0,65

0,34

0,01

64,66

0,98

100

1,00

34,02

1,00

65,98

1,00

Bilan

64

Distillat

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

II. Temprature de Tte
Colonne: Rectification dEssence Sauvage: TTte = TRose du Distillat sous
liquide du 1ier Plateau; ce Liquide= D aliment

P = 20 bar ( cet endroit, vapeurs de D

en reflux aprs condensation)
P
= PA
x
Raoult
Pression du mlange

TRose : v

+ PB (1 x )

Fraction Molaire du compos le + lger, A

l, donc: PA . x = P. y
liquide, A

ou x =

vapeur, A

P. y
P
et PB (1 x) = P(1 y ) ou (1 x ) = (1 y )
PA
PB
I

II

P. y P
y
et: I + II = 1 =
+
(1 y ) ou encore: i = 1 au point de Rose, TRose
PA
PB
i Ki
1

KA

KB

que lon vrifie par approximations successives (rappel: P = 20 bar)

TRose = 60C
y

TRose = 55C
y/K

y/K

Propane

0,99

1,1

0,9

0,99

1,01

0,98

Butane

0,01

0,54

0,02

0,01

0,51

0,02

0,92
65

JPG, 2CAD2, 2012-13

1,00

Conception & Opration dune Colonne

III. Temprature de Fond
Colonne: Rectification dEssence Sauvage: TFond = TBulle du Rsidu sous
P = 20 bar
TBULLE: l

v; la composition de cette vapeur: y =

liquide, A

et: I + II = 1 ou

PA
P
. x et 1 y = B .(1 x )
P vapeur, A
P
I

II

PA
P
.x + B (1 x ) = 1 et en posant K: K A .x + K B .(1 x ) = 1 ou
P
P

K .x
i

=1

On procde de mme que prcdemment:

TBulle = 100C
x

TBulle = 98C
x

K.x

K.x

Propane

0,02

1,75

0,035

0,02

1,70

0,034

Butane

0,98

1,00

0,98

0,98

0,99

0,970

1,015

66

JPG, 2CAD2, 2012-13

1,004

Conception & Opration dune Colonne

IV. Temprature dAlimentation
Colonne: Rectification dEssence Sauvage: linstallation est telle que la
TAlimentation = 75C. Cette T est-elle > ou > TBulle de lAlimentation ?
On essaie TBulle = 85C qui est la TTte

K.x

Propane

0,35

1,45

0,505

Butane

0,65

0,80

0,520
1,025

K.x = 1,025 est proche de 1, donc TBulle de lAlimentation ~85C

A 75C et sous P = 20 bar, la Charge est Liquide

67

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
1. Introduction: la Slectivit de la Sparation est lie au Taux de Reflux & au Nombre de
Plateaux. De nombreuses mthodes existent. Pour les Mlanges binaires, celle de McCabeThiele est la plus utilise
Warren Lee McCabe (1899 - 1982) & Ernest W. Thiele (1895 1993): 2 tudiants
en Master & Doctorat au MIT (Chemical Engineering) en 1925
Autres Mthodes:
Ponchon-Savarit (plus rigoureuse mais plus difficile)
Mthodes numriques lquilibre (MESH): rsolution simultane des quations de:
M: Bilans Matires
E: Relations dquilibre entre lments finis en quilibre ( Stages )
S: Sommes sur les fractions molaires, i.e. forcer: xi = 1
i
H: Bilans nergtiques & Enthalpiques (H)
Mthodes hors-quilibre (MERSHQ):
R: Transferts de Matire & Energie
Q: quations dquilibres

Stages

68

hors-quilibre

JPG, 2CAD2, 2012-13

quilibre

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
1. Introduction:
Caractristiques & Limitations de la Mthode de McCabe-Thiele:
Calcul du nombre de plateaux thoriques pour la sparation de mlanges binaires
Solution graphique base sur les diagrammes dquilibre (x-y), trs visuelle
Acceptable pour de nombreuses applications, populaire par sa simplicit
(~couteau Suisse)
Valeur pdagogique, en brainstorming (remue-mninge) & on vite la Black-Box !
Manque de riguer d, inter alia, la non-utilisation des bilans enthalpies
P est suppose constante le long de toute la colonne
Base sur les hypothses de Lewis:

69

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
2. Hypothses de Lewis:

- suppose Dbits Molaires des Flux Saturs Constants

- implique Taux de Reflux & Rebouillage Constants
- valides si volatilits proches & mme famille chimique
- permet dliminer les quations de bilan enthalpiques
Rectification

70

Epuisement

Vo

JPG, 2CAD2, 2012-13

L L0
L L + A
Taux de Reflux = L0/D = Cste

V V0
V V +VA
Taux de Rebouillage = V0/R = Cste

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
a) Donnes ncessaires:
Caractristiques de la Charge A, zA et son tat de vaporisation: VA & LA
Spcifications de la Sparation: xD & xR
quilibres l,v du Mlange trait sous la Pression opratoire (suppose constante)
par le diagramme y x
Le taux de Reflux rf = L0/D
Ensuite on procde par tapes:
Zone de Rectification
Zone dpuisement (Stripping)
Combinaison des zones de Rectification & puisement
Zone dAlimentation

71

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
b) Bilan en Zone de Rectification: en Tte de Colonne: y1 = xD si condenseur total

les quations de bilans liant les flux liquide et

vapeur entre 2 plateaux p & p+1 scrivent:
V = L + D (avec L = L0)
Vyp+1= Lxp+ DxD
soit: y p +1 =
et

L
Dx
L
L
DxD
x p + D avec rf = , y p +1 =
xp +
,
V
V
D
L+D
L+D
r
x
y p +1 = f x p + D
rf + 1
rf + 1

Rectification Operating Line - Droite Bilan de Rectification, DBR

si rf est fix, les dbits liquide et vapeur sont

dtermins

72

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
b) Bilan en Zone de Rectification: Rectification Operating Line Droite Bilan de
Rectification

rf = 0

La DBR:

aucun Reflux

rf

Reflux ~ total
(Distillation des livres de CP)

73

JPG, 2CAD2, 2012-13

Passe par x = y = xD

Pente:

Ordonne lorigine:

rf
rf +1

ou

L
V

xD
rf +1

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
c) Bilan en Zone dpuisement (Stripping):

les quations de bilans liant les flux liquide et

vapeur entre 2 plateaux q-1 & q et p+1 scrivent:
V = L + R
Lxq = Vyq-1+ RxR
soit: yq 1 =
et,

L'
Rx
V'
V '+ R
Rx
xq R ,avec rb = , yq 1 =
xq R
V'
V'
R
V'
V'
r +1
x
yq 1 = b
xq R
rb
rb

Stripping Operating Line - Droite Bilan d'Epuisement, DBE

si rb est fix, les dbits liquide et vapeur sont

dtermins

74

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
c) Bilan en Zone dpuisement (Stripping): Stripping Operating Line Droite
Bilan d puisement

La DBE:

75

JPG, 2CAD2, 2012-13

Passe par x = y = xR

Pente:

Ordonne lorigine:

rb + 1
rb

xR
rb

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
d)Raccordement des 2 Droites de Bilan: lintersection des droites bilans de
Rectification (DBR) & dpuisement (DBE) possde des proprits remarquables
DBR: Vy = Lx + DxD
DBE: V ' y = L ' x RxR

pour yp+1, yq1 et xp, xq on a:

(1)
(2)
PENTE

et: (1) - (2) (V V ') y = ( L L ') x + DxD + RxR , soit: VA y = LA x + Az A ou encore: y =

DBR & DBE ne dpend que de lAlimentation !
passe par x = zA & y = zA
sa pente, -LA/VA , dpend de l tat de vaporisation de
la charge
si la charge est
o TBulle: VA 0 & LA A, droite VERTICALE
(pente )
o TRose : LA 0 & VA A, droite HORIZONTALE
(pente 0)
76

JPG, 2CAD2, 2012-13

LA
Az
x + A (q-line)
VA
VA

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
4.

Construction de McCabe-Thiele en pratique:

a) Donnes du Problme:
Charge: A, zA, LA & VA (q = LA/LA +VA = f(T), fixe la Droite dAlimentation ou q-line)
Spcifications des Produits: xD & xR qui donnent D & R
Taux de Reflux (rf = L/D): fixer les trafics l (L, L) & v (V, V) dans la colonne + DB
b) Mise en Place des Droites de Bilan (DBR & DBE)
Zone de Rectification

Zone dpuisement

(xm+1, ym): Spcification

calculable, vide infra

Spcification

77

xn > xA > xn+1

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
c) Remarques:
Dernier Plateau de la Rectification

G
A
P

Position Correcte dAlimentation

xn > xA > xn+1

xn+1: 1 Plateau aprs n, xn+1 < LA, Alimentation

TBulle (Charge): vapeur pas modifie par la
charge (V = V & yn+1 = ym) GAP 0
TRose (Charge): liquide pas modifi par la
charge (L = L) & (xn, ym) sur la DBE

Premier Plateau de lpuisement

(xn, yn+1) DBR & (xm+1, ym) DBE

Calcul de ym

on a en faisant le Bilan Matire des Flux LA & L:

L ' = LA + L
L x + Lxn
soit : xm+1 = A A
et le point ( xm+1, ym ) DBE, qui dtermine ym

L
'
x
=
L
x
+
Lx
L
L
+
A A
n
m+1
A

78

JPG, 2CAD2, 2012-13

Effet de la Composition (l, v) de la Charge

a) q > 1 (subcooled l, <TB)

d) q = 0 (saturated v, TR)

q = fraction molaire de l dans la Charge

b) q = 1 (saturated l, TB)

e) q < 0 (superheated v, >TR)

c) 0 < q < 1 (mix of l & v)

79

JPG, 2CAD2, 2012-13

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
- En rsum

!: IN =

OUT =

Slope = R/(R+1)
Slope = q/(1-q)
Slope = L/ V
-

80

JPG, 2CAD2, 2012-13

Toutes ces oprations sont effectues

par les logiciels courants (ChemSep)
On entre les donnes du pb et les
approximations
!: effet black-box

Conception & Opration dune Colonne

V. Taux de Reflux & Nombre de Plateaux McCabe-Thiele
Rvision

Voir:
Tutorial Distillation Basics dans: [Link]
Raffinage Ptrolier (2 vidos):

&
Distillation [Link]

81

Distillation [Link]

JPG, 2CAD2, 2012-13

Applications en Raffinage Ptolier

I. Distillation Atmosphrique
Caractristiques du Raffinage
Grande complexit des Charges
Gamme de Tebullition
Nature des charges (P, O, N, A)
Ractivit des Charges vs. T
Quantits traiter
Nombreuses coupes Soutirages multiples

82

JPG, 2CAD2, 2012-13

Applications en Raffinage Ptolier

I. Distillation Atmosphrique
Schma Gnral dune raffinerie de Ptrole

83

JPG, 2CAD2, 2012-13

Applications en Raffinage Ptrolier

I. Distillation Atmosphrique
Vue Gnrale

84

JPG, 2CAD2, 2012-13

Applications en Raffinage Ptrolier

I. Distillation Atmosphrique

85

JPG, 2CAD2, 2012-13

Applications en Raffinage
II. Distillation sous Vide
Ncessit de Distillation sous vide:
Ractivit des HC haute Tbullition
Si PHC diminue, TEbullition diminue
Augmentation des carts de volatilit
Diminution de la consommation nergtique

86

JPG, 2CAD2, 2012-13

Symboles Graphiques

87

JPG, 2CAD2, 2012-13

Symboles Graphiques

88

JPG, 2CAD2, 2012-13

Equipements & Instruments

89

JPG, 2CAD2, 2012-13