République Algérienne Démocratique et Populaire
Université M’Hamed Bougara Boumerdes
Faculté des Hydrocarbures et de la Chimie
Département Gisement Minier et pétrolier
Module : Thermodynamique des
fluides Pétroliers
TD Exercice 03 & 04
Application des Equation d’état en utilisant les méthodes de calcule
Encadreur Réalisé Par
Madame MEFTAH Ahmed Merzoug
Abderraouf Chemmakh
Année Universitaire 2018-2019
MAPH15
TD Thermodynamique Exercice 3 & 4 Ahmed Merzoug
Chemmakh Abderraouf
Exercice 03
- Calcule de Pression
Pour que le propane existe à l’état Vapeur humide P=Pv
Pv= 185 D’après la projection dans le diagramme de Cox
D’après la loi de Van Der Waals
Tel que
P : pression du system.
R : Constante des gaz parfaits 10.73 psi.ft3/lb-mole, OR.
T : Température du System.
V : volume, ft3/mole.
a : Paramètre d’attraction
b : Paramètre de répulsion.
On a
Z3 − (1+ B) Z2 + AZ − AB = 0 ________________________________________(1)
- Calcule des paramétré a et b :
𝑅2 ∗ 𝑇𝑐2
𝑎 = 0.421875 ∗ = 34951,17216
𝑃𝑐
𝑅 ∗ 𝑇𝑐
𝑏 = 0.125 ∗ = 1,4491522
𝑃𝑐
Tc=666 K
Pc= 616 psi
- Calcule de A et B
𝑎𝑃
𝐴= = 0,179122
𝑅2 𝑇 2
𝑏𝑃
𝐵= = 0,044625
𝑅𝑇
MAPH15 Année Universitaire 2018-2019
TD Thermodynamique Exercice 3 & 4 Ahmed Merzoug
Chemmakh Abderraouf
Tel que P=185 psi et T= 100 F Alors T= 560 R
On substitue A et B dans l’équation 1
On aura :
Z3 − 1.044625Z2 + 0.179122Z − 0.007993= 0
On utilise la calculatrice Kajib fx991ES Plus suivant les étapes suivante :
1- Clique sur Mode.
2- Clique sur 5 (Mode EQN).
3- On choisit 4 Mode d’équation d’un variable d’ordre 3.
4- On remplit les facteurs a, b, c, d selon la puissance d’ordre décroissant
(1, 1+B, A, -AB) et on clique sur égale.
5- On clique sur égale 3 fois pour avoir les résultats à chaque fois qu’on
clique on trouve une nouvelle solution
Z1=0.075 Z2=0.84 Z3=0.125
On choisit les résultats plus grand et plus petit d’où
Zv=0.84 Zl= 0.075
Calcule de la masse Volumique
𝑃𝑀
ρ= ; ρg=1.61 lb/ft3 ; ρl=17.98 lb/ft3
𝑍𝑅𝑇
Tel que : P= 185 psi M=44g/mole R=10.73 psi.ft3/lb-mole, OR. T=560 R
MAPH15 Année Universitaire 2018-2019
TD Thermodynamique Exercice 3 & 4 Ahmed Merzoug
Chemmakh Abderraouf
Exercice 04
Du tableau on peut avoir
Loi des états correspondants
Mi Ppc M
g/mole Yi Tci R Pci kPa Tpc R kPa g/mole
16 0.9443 343.08 4600 323.970444 4343.78 15.1088
30.1 0.0158 549.72 4880 8.685576 77.104 0.47558
44.1 0.0132 665.64 4250 8.786448 56.1 0.58212
58.1 0.0051 734.58 3650 3.746358 18.615 0.29631
58.1 0.0053 765.36 3790 4.056408 20.087 0.30793
72.2 0.0023 830.52 3380 1.910196 7.774 0.16606
72.2 0.014 845.46 3370 11.83644 47.18 1.0108
362.99187 4570.64 17.9476
Tpc=∑ 𝑦𝑖𝑇𝑐𝑖=362.9918 R
Ppc=∑ 𝑦𝑖𝑃𝑐𝑖= 4570,64 Psi
Tpr=T/Tpc= 1,875165964
Ppr=P/Ppc= 4,92272417
Tel que P=22.5 Mpa T=680 R
On projetant sur le graph de Standing on trouve Z=0.93
Calcule de la masse volumique :
𝑃𝑀
ρ= ρ=8.63 lb/ft3
𝑍𝑅𝑇
MAPH15 Année Universitaire 2018-2019
TD Thermodynamique Exercice 3 & 4 Ahmed Merzoug
Chemmakh Abderraouf
En utilisant EOS
Redlich Knwong Equation of State
P : pression du system.
R : Constante des gaz parfaits 10.73 psi.ft3/lb-mole, R.
T : Température du System.
V : volume, ft3/mole.
a : Paramètre d’attraction
b : Paramètre de répulsion.
Z3 – Z2 + (A – B – B2)Z – AB = 0__________________________________(2)
Mi Yi Tci (R) Pci Kpa Tpc (R) Ppc kPa M lb/mole ai bi
16 0.9443 343.08 4600 323.970444 4343.78 15.1088 745.4504974 0.07747
30.1 0.0158 549.72 4880 8.685576 77.104 0.47558 1804.052166 0.11701
44.1 0.0132 665.64 4250 8.786448 56.1 0.58212 3037.216367 0.16269
58.1 0.0051 734.58 3650 3.746358 18.615 0.29631 4306.963319 0.20905
58.1 0.0053 765.36 3790 4.056408 20.087 0.30793 4502.753264 0.20977
72.2 0.0023 830.52 3380 1.910196 7.774 0.16606 5945.239137 0.25524
72.2 0.014 845.46 3370 11.83644 47.18 1.0108 6179.339656 0.26060
362.99187 4570.64 17.9476
Pc =662.82 psi
On calcule les valeurs am bm pour les mélanges :
am = 262,240783 bm = 0,046426
MAPH15 Année Universitaire 2018-2019
TD Thermodynamique Exercice 3 & 4 Ahmed Merzoug
Chemmakh Abderraouf
- Calcule de A et de B
P= 3263.34 Psi T=680 R R=10.73 psi.ft3/lb-mole, R.
A=0,00058022 B=0,020264244
En remplaçant A et B dans 2 on trouve :
Z3 – Z2 + (-0,020094664)Z – 1,17577E-05 = 0
On utilise la calculatrice Kajib fx991ES Plus suivant les étapes suivante :
1. Clique sur Mode.
2. Clique sur 5 (Mode EQN).
3. On choisit 4 Mode d’équation d’un variable d’ordre 3.
4. On remplit les facteurs a, b, c, d selon la puissance d’ordre
décroissant (1, -1, A – B – B2, -AB) et on clique sur égale.
5. On clique sur égale 3 fois pour avoir les résultats à chaque fois qu’on
clique on trouve une nouvelle solution
Z1=-0.019 Z2=1.019 Z3=-6.032*10^-4
On choisit
Zv=1.019
Calcule de la masse Volumique
ρ = PM/ZRT ; ρg=7.87 lb/ft3
MAPH15 Année Universitaire 2018-2019
