Laboratory Catalysis

Laboratoire Catalyse&&
Spectrochemistry
Spectrochimie

Introduction à la Distillation
2CAD2: CM = 10

Jean-Pierre GILSON
Tél.: [Link].15 – email: [Link]@[Link]

ENSICAEN, Université de Caen, CNRS

6 bd du Maréchal Juin, 14050 Caen, France
http:/[Link]
 Introduction
Procédés de Séparation
But du Cours

E. Sorel, La Rectification de l’Alcool, 1893

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

Laboratoire (Batch)

2 JPG, 2CAD2, 2012-13

 Plan du Cours
Procédés de Séparation
I. Introduction: Nécessité, types de séparations, Labo vs. Industrie,…
II. Aspects de la distillation: Bilan matière, énergie, phénomènes de transport…

Équilibres Liquide - Vapeur

I. Thermodynamique
II. Simulations des équilibres

Distillation Flash & Continue

I. Distillation Flash
II. Distillation Continue
III. Phénomènes de Transport
IV. Équipement

Conception & Utilisation de Colonnes

I. Bilans Matières
II. Température de Tête
III. Température de Fond
IV. Température d’Alimentation
V. Taux de Reflux & Nombre de Plateaux (MacCabe-Thiele)

Applications en Raffinage Pétrolier

I. Distillation atmosphérique
II. Distillation sous vide

3 JPG, 2CAD2, 2012-13

 Bibliographie
1. E. Sorel, La Rectification de l‘Alcool, (1894), Gauthier-Villars et fils, Paris

2. J-P. Wauquier, Le Raffinage du Pétrole – 2. Procédés de Séparation, (1998), Éditions Technip, Paris

3. J. Vidal, Thermodynamique: Application au Génie Chimique & à l’Industrie Pétrolière, (1997), Éditions
Technip, Paris

4. D. Morvan, Les Opérations 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 Données

1. ChemSep ([Link]): Gratuit en version lite. Effectue toutes les

opérations décrites dans ce cours

2. Dortmund Data Bank ([Link]):Base de données thermodynamiques et

logiciel accessible gratuitement (versions explorer et académique lite)

5 JPG, 2CAD2, 2012-13

 Procédés de Séparation
I. Introduction
Séparations: Pourquoi ?
a. Fractionner: n’envoyer vers des unités industrielles que les « coupes » nécessaires
b. Recycler: les procédés de conversion (catalytique) n’ont pas toujours une conversion
maximum et la sélectivité désirée; on recycle ce qui n’est pas aux spécifications
c. Purifier: les opérations en aval ont des exigences (contaminants, polluants…) sur
leurs « charges »

a. & c. b.
6 JPG, 2CAD2, 2012-13
 Procédés de Séparation
I. Introduction
Séparations: Comment ?

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

7 JPG, 2CAD2, 2012-13
 Procédés de Séparation
I. Introduction
Séparations: Utilisation & Maturité

8 JPG, 2CAD2, 2012-13

 Procédés de Séparation
I. Introduction
Pétrole Brut (Fractionationner) Procédés (Recycler) Produits (Purifier)

9 JPG, 2CAD2, 2012-13

 Procédés de Séparation
I. Introduction
Capacités de Distillation de Pétrole Brut
Raffinerie Location Barrels per Day*
Jamnagar (Reliance Industries
Jamnagar, India 1 240 000
Limited)
Paraguaná,
Paraguana (PDVSA) 940 000
Falcon, Venezuela
Ulsan (SK Energy) Ulsan, South Korea 850 000
GS Caltex Yeosu Refinery (GS
Yeosu, South Korea 730 000
Caltex)
ExxonMobil Singapore 605 000
Port Arthur (Motiva Enterprises, i.e.
Port Arthur, Texas, USA 600 000
RD Shell & Saudi Aramco)
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
Garyville, LA 490 000
Petroleum)
Pernis (RD Shell) Pernis, Pays-Bas 416 000
Antwerp (Total) Anvers, Belgique 360 000
Normandie (Total) Gonfreville, France 205 000
10 JPG, 2CAD2, 2012-13
*: 1 barrel = 159 l; ÷ par 6,25 pour convertir en m3
 Procédés de Séparation
II. Aspects de la Distillations
Distillation: Paramètres ?
Séparation par EQUILIBRE: Thermodynamique & Transfert de Matière
Séparation par CHAUFFAGE: Thermodynamique & Transfert de Chaleur
Fonctionnement en CONTINU: Bilan & Transfert de Matière
Dimension INDUSTRIELLE: Intégration sur site, Contrôle, É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 système fermé

+dni Evolution (α → β) si dGTot < 0, soit µiα > µiβ

β µiβ
Equilibre entre phases si dGTot = 0, soit µiα = µiβ, i.e. égalité des µi de
chaque composant i dans toutes les phases.
α µiα dP ∆H
-dni Relation de Clapeyron: =
dT T ∆V0
P ∆H  1 1 
Clausius-Clapeyron: ln 0 =  − 
P R  T0 T 
Mélanges:
Degrés de Liberté
1 0 0 2 2
Variance : v = ( c − r − k ) + n − φ =1
Nombre de variables - contraintes Constituants Equations chimiques indépendantes Autres conditions imposées 2 pour T & P Phases

Loi de Dalton (Gaz Parfaits): PTot = ∑ Pi et Pi = yi .PTot

i

Loi de Raoult (Solution idéale): PA = [Link]*

Loi de Henry (Solution idéale diluée): PA = xA .K A

12 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
I. Thermodynamique
2. Corps Pur (ex. H2O)

14 JPG, 2CAD2, 2012-13

 Équilibres L - V
I. Thermodynamique
2. Corps Purs: H2O

Courbe de tension de vapeur

1 2 PointTriple – PointCritique
3 4

Chaleur Latente

Chaleur Sensible

Dépenses & Pertes énergétiques

Distiller est Energivore !
15 JPG, 2CAD2, 2012-13
 Équilibres Liquide - Vapeur
I. Thermodynamique
3. Mélanges Binaires:
xi ≠ yi: sauf dans le cas d’un Azéotrope (+ ou -) avec:
xi ≡ fraction molaire de i en phase liquide P = PB0 + ( PA0 + PB0 ) .xA
yi ≡ fraction molaire de i en phase vapeur PA0 .PB0
P= 0
zi ≡ fraction molaire de i (phase liquide + vapeur) P A + ( PB
0
− P 0
A ). y A

Bulle

Ligne Conodale

Rosée

16 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
3. Mélanges Binaires:
Diagramme T-H: les lignes isothermes ⇒ courbes des T-H des constituants
Evolution de T en fonction de la Chaleur absorbée

Corps Pur : n-Butane Mélange Binaire

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

Remplacement des paliers par des courbes
AB (A’B’, A’’B’’…): vaporisation (Chaleur Latente)

17 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Principe de la Distillation (Batch)

T1: T1 < TébA&B, x1 = x2 T3: T3 > TébB & T3 < TébA, x1 = x2 T4: T4 = TébA, x1 = x2
T2: Ébullition commence, La vapeur (distillat) est enrichie en La dernière goutte de liquide se
composition de la vapeur B et les compositions x3 & y3 se vaporise et sa composition, y4, est A
(distillat), enrichie en B est y2 déplacent vers la droite p/γ T2 pur

A & B purs: une ∞ d’étapes (Plateaux) est nécessaire. Plus facile de mélanger que de séparer !
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 départ
18 JPG, 2CAD2, 2012-13
 Équilibres Liquide - Vapeur
I. Thermodynamique
Règle du Levier

- Système binaire de composition 1

- A la Tconodale, la composition est:

liquide: 2
vapeur: 3

- Règle du levier:

Liquide = a ∼ 0,25 ici

Vapeur b

19 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Principe de la Distillation Fractionnée

Départ:
xA = 0,8 & xB = 0,2

R = Résidu
D = Distillat
20 JPG, 2CAD2, 2012-13
 Équilibres Liquide - Vapeur
I. Thermodynamique
Principe de la Distillation Fractionnée - Résumé: 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 intermédiaire, peu utiles…

Séparation totale: ∞ Plateaux ! (Plus

facile de mélanger [S ↑] que de
séparer [S ↓])

Rendements + Succession d’étapes

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

DÉFINITIONS
Distillat
Vapeur condensée enrichie en l’élément le plus
volatil

Résidu (Fond)
Reste du bouilleur, appauvri en composé le plus
volatil

BILANS MATIÈRE
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 Mélanges Azéotropiques:
On butte sur l’Azéotrope !
Formellement: Σ de 2 Mélanges: 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
Résumé: Idéal – Déviation - Azéotropie

26 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Résumé: Azéotropie

27 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
4. Mélanges Ternaires:

Données d’Équilibre reportées sur un Triangle Équilatéral

Chaque point dans le triangle (R, S) y/c arêtes (L, P, Q) &

sommets (A, B, C) représente une composition du mélange
La Σ des longueurs des⊥ abaissées d’un point intérieur (eg. R)
sur les 3 côtés 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. délimite la zone d’immiscibilité
R: à l’intérieur de la zone d’immiscibilité; se scinde en 2
phases de composition M & P (ligne d’équilibre 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. Mélanges Complexes: les charges & produits pétroliers contiennent une grande
variété d’hydrocarbures (P, N, A, O) de poids moléculaire très variable (C1–C100+).
L’apparition de charges renouvelables (biomasse) va encore complexifier la situation
(oxygénés). 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 = TRosée

Tx
Ty

PI = TBulle

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

29 JPG, 2CAD2, 2012-13
 Équilibres Liquide - Vapeur
I. Thermodynamique
Méthodes d’évaluation des Mélanges complexes: Méthodes standardisées (ASTM,
ISO, API…) qui sont soit rapides et peu précises ou lentes et (plus) précises
Méthode Distillation Flash (FC*): séparation la plus simple, à un plateau
(routine) - *:Flash Curve

Discontinu (Batch) ou Continu (F, z)

La charge est chauffée sous PZ, à une TZ1
Détente (isenthalpique)dans le Ballon de Flash (Flash Drum)
Vapeur enrichie en volatils & Liquide enrichi en lourds
On répète à une nouvelle TZ2 et sous PZ
On peut recommencer à une P différente

30 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Méthode Distillation Flash (FC): On peut mettre en CASCADE (Série) des Ballons
de Flash à P (ou T) constante pour améliorer les séparations – Consommation d’Energie !

Diagramme T, x déjà rencontré

Meilleure séparation car plus de plateaux
Coût important: Chauffer & Refroidir
Retenons le principe !
P = Constante
31 JPG, 2CAD2, 2012-13
 Équilibres Liquide - Vapeur
I. Thermodynamique
Méthode True Boiling Point (TBP): Séparation plus longue, plus précise

32 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Comparaison True Boiling Point (TBP): mélange 4 Constituants vs. Pétrole Brut

33 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Sélectivité (S) de Séparation: Overlap & Gap des Mélanges complexes

a: courbe TBP d’un mélange binaire b: courbe ASTM d’un mélange binaire
S=1 S=α

34 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Sélectivité (S) de Séparation: Overlap & Gap des Mélanges complexes
c: Courbe ASTM d’un Brut Pétrolier d: Courbe ASTM d’un Brut Pétrolier
S = Overlap S = Gap

35 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Comparaisons FC, ASTM & TBP

TBP: Séparation plus prononcée (∆TTBP > ∆TFC ) et donc de meilleure qualité

36 JPG, 2CAD2, 2012-13

 Équilibres Liquide - Vapeur
I. Thermodynamique
Distillation Simulée: Analyse Chromatographique en Phase Gaz

Étalonnage par n-paraffines de Téb connues

La colonne GC sépare par Téb
Mesurer l’aire entre 2 paraffines
Reporter sur un diagramme T = f(% Distillé)
ASTM D 2887: fractions dont Téb < 538°C
ASTM D 3710: fractions dont Téb < 250°C

Remarque:
- Il n’est pas encore possible d’analyser complètement
(moléculairement) les fractions pétrolières
- Même si cela était possible, la quantité d’isomères
obtenus nécessiterait 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. Résultat Désiré: 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 réarrangées:

40 JPG, 2CAD2, 2012-13

 Distillation Continue
I. Distillation Flash + Volatil

3. Cascades de Flash à P Constante

- Volatil
Boîtes Bleues ≡ Ballons de Flash L & V échangent Chaleur L & V échangent Matière
Liquide ↓: doit être réchauffé ⇒ Échangeur de Chaleur ⇒ Exit Échangeur de Chaleur
Vapeur ↑: doit être refroidie Système en mouvement:
⇒ Condenser une partie de V1 & Recycler
⇒ Vaporiser une partie de L5 & Recycler
41 JPG, 2CAD2, 2012-13
 Distillation Continue
I. Distillation Flash
4. Distillations Batch: Rectification & Stripping

Rectification Stripping
Seuls les composés légers sont récupérés Seuls les composés lourds sont récupérés
eg.: Production d’alcools (Calvados…) eg.: Enlever de (petites quantité de) solvant
Peu utilisé industriellement (casser des
azéotropes par ajout d’un 3ième composant –
eg.: MeOH/Tol + Et3N)
42 JPG, 2CAD2, 2012-13
 Distillation Continue
II. Distillation Continue
1. Description: Nous venons de décrire le Principe de la Distillation Continue !

TColonne

[Link]
43 JPG, 2CAD2, 2012-13
 Distillation Continue
II. Distillation Continue
1. Description: Instrumentation & Régulation de la colonne

44 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
2. Bilans Matière & 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)
Résidu: C5-10 , Essence “Stabilisée” (10,02 t/h) ⇒ Base pour Carburant (Spécification RVP)

Bilan Thermique
- Chaleur reçue
Enthalpie d’alimentation
Chaleur fournie par le rebouilleur
Enthalpie du reflux froid
- Chaleur extraite
Enthalpie des vapeurs de tête (D + Reflux)
Enthalpie du Résidu

45 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
Schéma:
d & e: Corps clés de la Séparation

C
o
r
p
s

C
l
é
s

46 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
3. Réglage d’une colonne: Tour de Stabilisation d’ essence
a) But de la Stabilisation:
Ajuster la tension de vapeur de l’essence stabilisée à la norme en vigueur (Reid
Vapor Pressure)
Veiller à ce que l’essence sauvage ne contienne pas de C5

b) Moyens disponibles:
Tour de Stabilisation permettant de ttaiter un certain Débit horaire de charge
Pression peut être ajustée
Température peut être ajustée
Taux de Reflux peut être ajusté

c) Influence de P
PMIN : Le Distillat doit pouvoir être condensé totalement à la T de l’eau de
refroidissement disponible dans le condenseur (Tension de vapeur de l’ essence
sauvage ~30°C : 10 kg/cm2, i.e. 9,8 bar
Si la pression d’opération 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. Règlage d’une colonne: Tour de Stabilisation d’ essence
d) Influence de T: - TTETE, T en tête de tour règle la qualité du produit de tête, 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 s’alourdit (s’enrichit en produits moins volatils); son point sec ↑

↓: le Distillat s’allège (s’appauvrit en produits plus volatils); son point sec ↓
T↓
Réglage de la TTETE: en agissant sur le débit de reflux froid qui, extrayant plus ou
moins de calories au système, fait varier TTETE ⇒ Qualité de l’essence sauvage

- TFOND, T en fond de tour, règle la qualité du produit de fond , le Résidu (essence

stabilisée). 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 Résidu

T ↑: le Résidu s’alourdit, sa teneur en butane diminue; sa tension de vapeur ↓

↓: le Résidu s’allège, sa teneur en butane augmente; sa tension de vapeur ↑
T↓
Réglage de la TFOND: en agissant sur le débit de produit chauffant (eau en
général) dans le rebouilleur

48 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
3. Réglage d’une colonne: Tour de Stabilisation d’ essence
e) Influence du Reflux: La chaleur fournie au système provient du rebouilleur. Elle sert à chauffer
(chaleur sensible), vaporiser (chaleur latente) le produit à distiller et vaporiser le reflux froid injecté en
tête de tour.
Reflux: extrait du système une quantité de chaleur excédentaire p/γ quantité
nécessaire au chauffage & vaporisation du distillat. Comme TTETE (= TROSEE) &
TFOND (= TBULLE) sont constantes, le débit du reflux est variable ⇒ variation du flux
de calories fourni par le rebouilleur
Reflux externe & interne: permet le transfert de matière (lavage des vapeurs) et
de chaleur. C’est la clé du Fractionnement.
Taux de Reflux: L/D permet de contrôler la Sélectivité et le nombre d’Étages de
la Distillation

Normand, Fig. 26 & 27, p. 149

49 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
3. Règlage d’une colonne: Tour de Stabilisation d’ essence
e) Influence du Reflux: L0: Reflux externe, débit 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: Débutaniseur

51 JPG, 2CAD2, 2012-13

 Distillation Continue
II. Distillation Continue
Exemple industriel: Débutaniseur

52 JPG, 2CAD2, 2012-13

 Distillation Continue
III. Phénomènes de Transport
1. Types de Phénomènes - Généralisation
∂C C: Concentration
Matière - Loi de Fick (2ième) : = D∇ 2C D: Diffusivité
∂t
∂T T: Température
∂2 ∂2 ∂2
Chaleur - Loi de Fourier: = κ∇ 2T κ: Conductivité thermique ∇
2
= + + = ∆
∂t u: Vitesse linéaire
∂x 2 ∂y 2 ∂z 2 Laplacien

∂u μ: Viscosité
Impulsion - Loi de Newton: = µ∇ 2u
∂t t: Temps
⇒ En commun: Transport d’une quantité physique (Matière, Chaleur, Impulsion) à travers un
gradient (Concentration, Température, Vitesse linéaire) afin d’éliminer ce dernier.
2. Transport de Matière
2 mécanismes sont responsables du transport de matière:
Diffusion Moléculaire: mouvement spontané µscopique dans une phase donnée (g, l, s); il
est assez lent. Voir ci-dessus
Diffusion Turbulente: mouvement macroscopique aléatoire d’un fluide
La Diffusion Moléculaire entre 2 phases (g, l) joue un rôle important en distillation. Elle se passe à
l’interface des 2 phases et conditionne les choix d’équipements.

Les coefficients de Diffusion Moléculaire [cm2/sec] sont de l’ordre de:

Gaz: Di = 10-1
Liquides: Di = 10-5
Solides: Di = 10-10 - 10-13

53 JPG, 2CAD2, 2012-13

 Distillation Continue
III. Phénomènes de Transport
3. Combinaison Transport de Matière & Chaleur
Phénomènes complexes (Turbulences): Fick & Newton (« linéaires ») se compliquent

Introduction de Flux (J), coefficients de transfert (k, Matière & Chaleur) et corrélations
(nombres sans dimensions, e.g. Re), plus utiles ⇒ Génie Chimique
J M ,q = k M ,q .(Interface).( ∆C ,T )
Matière, Chaleur Gradient de Concentration, Température
Coefficient de Transfert

Qualitativement: en distillation, on va favoriser les échanges, contacts (Matière Ji &

Chaleur JqL, JqV) dans les films (L & V) à l’Interface L/V

TCol

54 JPG, 2CAD2, 2012-13

 Distillation Continue
IV. Equipements
1. Rôle
Effectuer le plus efficacement possible le transfert de matière et de chaleur entre les
phases liquide & vapeur; ces transferts se passent aux interfaces et les équipement
devront favoriser l’étendue de celles-ci. En particulier, les équipements favoriseront:
temps de contact entre les phases
Qualitatif ! Métier de l’Ingénieur:
mélange approprié rendre ces termes QUANTITATIFS

2. Types de Colonnes
Colonnes à Plateau (TRAY):
L & G s’écoulent à contre-courant ou courant croisé
Colonne peut travailler sous Pression
Échanges de Matière (& Chaleur) s’effectuent sur des plateaux (Trays)

Contre-Courant Courant Croisé

55 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
Empêche le liquide de descendre trop rapidement
vers l’étage inférieur

√ 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 s’écoulent à contre-courant

Colonne peut travailler sous Pression
Échanges de Matière (& Chaleur) s’effectuent
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 céramique résistent mieux aux
milieux corrosifs
Moussage plus controllable
Fenêtre d’opération plus étroite
Nettoyage difficile et coûteux si contamination
par des particules solides
Éléments assez fragiles lors de stress en T ou P

Consulter, par exemple:

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

61 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
Cahier des Charges
Pour:
Concevoir (Design) de nouvelles colonnes

Conduire (Operate) des unités existantes

Réparer/Corriger (Troubleshooting) des dysfonctionnements

Nous avons besoin de:

Bilans Matières

Température de Tête de Colonne

Température de Fond de Colonne

Température d’Alimentation

Taux de Reflux & Nombre de Plateaux – Méthode de McCabe-Thiele

62 JPG, 2CAD2, 2012-13

Pour:
 Conception & Opération d’une Colonne
I. Bilans Matières
Colonne: Rectification d’Essence Sauvage (cf. Distillation Continue)
Charge:« Essence Sauvage », i.e., Propane & Butane sous P = 20 bar
√ Propane: 35 vol %
√ Butane: 65 vol %
Spécifications de la Séparation:
Propane: < 1 vol % Butane
Butane: < 2 vol % Propane
Soient:
F: nombre de molécules d’Alimentation (posons 100 molécules)
D: nombre de molécules de Distillat (Propane à < 1 vol % Butane)
W: nombre de molécules de Résidu (Butane à < 2 vol % Propane)
zF: fraction molaire d’un composant dan l’ Alimentation
xD: fraction molaire d’un composant dan le Distillat
xW: fraction molaire d’un composant dan le Résidu
⇒ 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 & Opération d’une Colonne
I. Bilans Matières
Colonne: Rectification d’Essence Sauvage
On dresse ainsi le tableau suivant:

Alimentation Distillat Résidu

[Link] zF [Link] 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

Bilan 100 1,00 34,02 1,00 65,98 1,00

64 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
II. Température de Tête
Colonne: Rectification d’Essence Sauvage: TTête = TRosée du Distillat sous
P = 20 bar (à cet endroit, vapeurs de D liquide du 1ier Plateau; ce Liquide= D alimenté
en reflux après condensation)
Raoult ⇒ P = PA x + PB (1 − x )
Pression du mélange Fraction Molaire du composé le + léger, A

P. y P
à TRosée : v l, donc: PA . x = P. y ou x = et PB (1− x) = P(1 − y ) ou (1− x ) = (1− y )
liquide, A vapeur, A
PA PB
I II
P. y P y
et: I + II = 1 = + (1− y ) ou encore: ∑ i = 1 au point de Rosée, TRosée
PA PB i Ki
1 1
KA KB
que l’on vérifie par approximations successives (rappel: P = 20 bar)

TRosée = 60°C TRosée = 55°C

y K y/K 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 1,00
65 JPG, 2CAD2, 2012-13
 Conception & Opération d’une Colonne
III. Température de Fond
Colonne: Rectification d’Essence Sauvage: TFond = TBulle du Résidu sous
P = 20 bar
PA P
à TBULLE: l v; la composition de cette vapeur: y = . x et 1− y = B .(1− x )
liquide, A
P vapeur, A P
II
I
PA P
et: I + II = 1 ou .x + B (1− x ) = 1 et en posant K: K A .x + K B .(1− x ) = 1 ou ∑ K .x i i =1
P P i

On procède de même que précédemment:

TBulle = 100°C TBulle = 98°C

x K K.x x K 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 1,004

66 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
IV. Température d’Alimentation
Colonne: Rectification d’Essence Sauvage: l’installation est telle que la
TAlimentation = 75°C. Cette T est-elle > ou > TBulle de l’Alimentation ?

On essaie TBulle = 85°C qui est la TTête…

x K 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 l’Alimentation ~85°C

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

67 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
1. Introduction: la Sélectivité de la Séparation est liée au Taux de Reflux & au Nombre de
Plateaux. De nombreuses méthodes existent. Pour les Mélanges binaires, celle de McCabe-
Thiele est la plus utilisée
Warren Lee McCabe (1899 - 1982) & Ernest W. Thiele (1895 – 1993): 2 étudiants
en Master & Doctorat au MIT (Chemical Engineering) en 1925
Autres Méthodes:
Ponchon-Savarit (plus rigoureuse mais plus difficile)
Méthodes numériques à l’Équilibre (MESH): résolution simultanée des équations de:
M: Bilans Matières
E: Relations d’Équilibre entre éléments finis en équilibre (« Stages »)
S: Sommes sur les fractions molaires, i.e. forcer: ∑ xi = 1
i
H: Bilans Énergétiques & Enthalpiques (H)
Méthodes hors-Équilibre (MERSHQ):
R: Transferts de Matière & Energie
Q: Équations d’Équilibres

« Stages »

68 hors-Équilibre JPG, 2CAD2, 2012-13 Équilibre

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
1. Introduction:
Caractéristiques & Limitations de la Méthode de McCabe-Thiele:
Calcul du nombre de plateaux théoriques pour la séparation de mélanges binaires

Solution graphique basée sur les diagrammes d’équilibre (x-y), très visuelle

Acceptable pour de nombreuses applications, populaire par sa “simplicité”

(~couteau Suisse)

Valeur pédagogique, en brainstorming (remue-méninge) & on évite la Black-Box !

Manque de riguer dû, inter alia, à la non-utilisation des bilans enthalpies

P est supposée constante le long de toute la colonne

Basée sur les hypothèses de Lewis:

69 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
2. Hypothèses de Lewis: - suppose Débits Molaires des Flux Saturés Constants
- implique Taux de Reflux & Rebouillage Constants
- valides si volatilités proches & même famille chimique
- permet d’éliminer les équations de bilan enthalpiques

Rectification
L ≡ L0
L’ ≡ L + A
⇒ Taux de Reflux = L0/D = Cste

Epuisement
V’ ≡ V’0
V ≡ V’ +VA
⇒ Taux de Rebouillage = V’0/R = Cste
V’o

70 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
a) Données nécessaires:
Caractéristiques de la Charge A, zA et son État de vaporisation: VA & LA

Spécifications de la Séparation: xD & xR

Équilibres l,v du Mélange traité sous la Pression opératoire (supposée constante)

par le diagramme y – x

Le taux de Reflux rf = L0/D

Ensuite on procède par étapes:

Zone de Rectification
Zone d’Épuisement (Stripping)
Combinaison des zones de Rectification & Épuisement
Zone d’Alimentation

71 JPG, 2CAD2, 2012-13

 Conception & Opération d’une 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 Tête de Colonne: y1 = xD si condenseur total

les équations de bilans liant les flux liquide et

vapeur entre 2 plateaux p & p+1 s’écrivent:

V = L + D (avec L = L0)
Vyp+1= Lxp+ DxD
L Dx L L DxD
soit: y p +1 = x p + D avec rf = , y p +1 = xp + ,
V V D L+D L+D
r x
et y p +1 = f x p + D
rf + 1 rf + 1
Rectification Operating Line - Droite Bilan de Rectification, DBR

⇒ si rf est fixé, les débits liquide et vapeur sont

déterminés

72 JPG, 2CAD2, 2012-13

 Conception & Opération d’une 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
aucun Reflux
La DBR:

- Passe par x = y = xD
rf L
- Pente: ou
rf +1 V
xD
rf → ∞ - Ordonnée à l’origine:
rf +1
Reflux ~ total
(Distillation des livres de CP)

73 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
c) Bilan en Zone d’Épuisement (Stripping):

les équations de bilans liant les flux liquide et

vapeur entre 2 plateaux q-1 & q et p+1 s’écrivent:

V’ = L’ + R
L’xq = V’yq-1+ RxR
L' Rx V' V '+ R Rx
soit: yq −1 = xq − R ,avec rb = , yq −1 = xq − R
V' V' R V' V'
r +1 x
et, yq −1 = b xq − R
rb rb
Stripping Operating Line - Droite Bilan d'Epuisement, DBE

⇒ si rb est fixé, les débits liquide et vapeur sont

déterminés

74 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
3. Principe de la Construction Graphique de McCabe-Thiele:
c) Bilan en Zone d’Épuisement (Stripping): Stripping Operating Line – Droite
Bilan d’ Épuisement

La DBE:

- Passe par x = y = xR
rb + 1
- Pente:
rb
− xR
- Ordonnée à l’origine:
rb

75 JPG, 2CAD2, 2012-13

 Conception & Opération d’une 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: l’intersection des droites bilans de
Rectification (DBR) & d’Épuisement (DBE) possède des propriétés remarquables

DBR: Vy = Lx + DxD (1)

pour yp+1, yq–1 et xp, xq on a: 
DBE: V ' y = L ' x − RxR (2)

PENTE

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

⇒ ∩ DBR & DBE ne dépend que de l’Alimentation !

passe par x = zA & y = zA

sa pente, -LA/VA , dépend de l’ État de vaporisation de

la charge

si la charge est à
o TBulle: VA ≅ 0 & LA ≅ A, droite VERTICALE
(pente ∞)
o TRosée : LA ≅ 0 & VA ≅ A, droite HORIZONTALE
(pente 0)

76 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
4. Construction de McCabe-Thiele – en pratique:
a) Données du Problème:
Charge: A, zA, LA & VA (q = LA/LA +VA = f(T), fixe la Droite d’Alimentation ou q-line)
Spécifications 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 d’Épuisement

(xm+1, ym): Spécification

calculable, vide infra

Spécification

xn > xA > xn+1

77 JPG, 2CAD2, 2012-13
 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
c) Remarques:

Dernier Plateau de la Rectification

Position Correcte d’Alimentation

xn > xA > xn+1
G
xn+1: 1 Plateau après n, xn+1 < LA, Alimentation
A TBulle (Charge): vapeur pas modifiée par la
P charge (V = V’ & yn+1 = ym) ⇒ GAP ≡ 0
TRosée (Charge): liquide pas modifié par la
charge (L = L’) & (xn, ym) sur la DBE

Premier Plateau de l’Épuisement

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

on a en faisant le Bilan Matière des Flux LA & L:

 L ' = LA + L L x + Lxn
 soit : xm+1 = A A et le point ( xm+1, ym ) ∈ DBE, qui détermine ym
L '
 m+1 x = L x
A A + Lxn L A + L

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 & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
- En résumé
!: IN = ⇒ OUT =

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

- Toutes ces opérations sont effectuées

par les logiciels courants (ChemSep…)
- On entre les données du pb et les
approximations…
- !: effet black-box

80 JPG, 2CAD2, 2012-13

 Conception & Opération d’une Colonne
V. Taux de Reflux & Nombre de Plateaux – McCabe-Thiele
Révision

Voir:
Tutorial « Distillation Basics » dans: [Link]

Raffinage Pétrolier (2 vidéos): &

Distillation [Link] Distillation [Link]

81 JPG, 2CAD2, 2012-13

 Applications en Raffinage Pétolier
I. Distillation Atmosphérique
Caractéristiques du Raffinage
Grande complexité des Charges
√ Gamme de Tebullition
√ Nature des charges (P, O, N, A)
√ Réactivité des Charges vs. T
√ Quantités à traiter

Nombreuses « coupes » ⇒ Soutirages multiples

82 JPG, 2CAD2, 2012-13

 Applications en Raffinage Pétolier
I. Distillation Atmosphérique
Schéma Général d’une raffinerie de Pétrole

83 JPG, 2CAD2, 2012-13

 Applications en Raffinage Pétrolier
I. Distillation Atmosphérique
Vue Générale

84 JPG, 2CAD2, 2012-13

 Applications en Raffinage Pétrolier
I. Distillation Atmosphérique

85 JPG, 2CAD2, 2012-13

 Applications en Raffinage
II. Distillation sous Vide
Nécessité de Distillation sous vide:
Réactivité des HC à haute TÉbullition
Si PHC diminue, TEbullition diminue
Augmentation des écarts de volatilité
Diminution de la consommation énergétique

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