Performance des Tours de Refroidissement
Performance des Tours de Refroidissement
Deux méthodes de calcul ont été développées pour le calcul de la performance de la tour
de refroidissement. La première méthode consiste à faire des bilans thermique et massique afin
de déterminer le rendement de la tour de refroidissement 33E101, pour la deuxième méthode
nous avons utilisé l’équation du Merkel pour la détermination du NUT (nombre d’unité de
transfert) ce qui nous a permis par la suite de calculer l’efficacité de la tour de refroidissement.
Les résultats obtenus ont montré que la deuxième méthode est plus représentative grâce
à sa précision. Néanmoins, les résultats des deux autres méthodes restent assez proches.
Une évaluation économique a été effectuée dans le but de déterminer la quantité d’eau
qui peut être économisée afin de réduire les pertes dues à l’évaporation de l’eau durant le
processus de refroidissement.
Les résultats ont montrés qu’il été avantageux de travailler avec un écart de température
de 10°C au lieu d’un écart de température plus grand, car cela nous permet de réaliser un gain
considérable.
Mots clés : température sèche, température humide, humidité relative, humidité absolue, NUT,
efficacité d’une tour de refroidissement.
I
DEDICACE
Je dédie ce projet :
A ma chère mère,
Qui n’ont jamais cessé, de formuler des prières, à mon égard, de me soutenir et de m’épauler
pour que je puisse atteindre mes objectifs.
II
REMERCIMENTS
III
TABLE DES MATIERES
RESUME ..................................................................................................................................... I
DEDICACE ............................................................................................................................... II
REMERCIMENTS................................................................................................................... III
TABLE DES MATIERES ....................................................................................................... IV
LISTE DES TABLEAUX .......................................................................................................VII
LISTE DES FIGURES .......................................................................................................... VIII
LISTE D’ABREVIATIONS .................................................................................................... IX
INTRODUCTION GENERALE ................................................................................................ 1
I. PRESENTATION DE LA RAFFINERIE D’ARZEW ...................................................... 2
I.1 INTRODUCTION ....................................................................................................... 3
I.2 CAPACITE DE TRAITEMENT ................................................................................. 3
I.3 PRESENTATION DES INSTALLATIONS DE PRODUCTION ............................. 4
I.3.1 DEPARTEMENT DE PRODUCTION P1 .......................................................... 4
I.3.2 DEPARTEMENT DE PRODUCTION P2 .......................................................... 4
I.3.3 DEPARTEMENT DE PRODUCTION P3 .......................................................... 4
I.3.4 DEPARTEMENT DE PRODUCTION P4 .......................................................... 4
I.3.5 LABORATOIRE DE CONTROLE ..................................................................... 5
I.4 PRESENTATION DE LA ZONE 03 (UTILITES) ..................................................... 5
I.4.1 UNITE 31 (PRODUCTION DE VAPEUR ET D’ELECTRICITE) .................... 5
I.4.2 UNITE 32 (PRODUCTION DE L’EAU DISTILLEE) ....................................... 6
I.4.3 UNITE 33 (TOUR DE REFROIDISSEMENT) .................................................. 6
I.4.4 UNITE 35 (RECEPTION ET DISTRIBUTION DU FUEL GAZ) ..................... 7
I.4.5 UNITE 36 (PRODUCTION D’AIR SERVICE ET D’AIR INSTRUMENT) ..... 8
I.4.6 UNITE 37 (PRODUCTION DU GAZ INERTE) ................................................ 8
II. THEORIES SUR TOURS DE REFROIDISSEMENT .................................................... 10
II.1 INTRODUCTION ..................................................................................................... 11
II.2 PRINCIPE DE TOUR DE REFROIDISSEMENT ................................................... 11
II.3 CLASSIFICATION DES TOURS DE REFROIDISSEMENT ................................ 12
II.3.1 TOURS DE REFROIDISSEMENT A TIRAGE NATUREL............................ 12
II.3.2 TOURS DE REFROIDISSEMENT A TIRAGE MECANIQUE ...................... 13
II.4 DIFFERENT TYPE DES CIRCUITS DE REFROIDISSMENT ............................. 15
II.4.1 CIRCUIT OUVERT ........................................................................................... 15
II.4.2 CIRCUIT SEMI-OUVERT ................................................................................ 16
IV
II.4.3 CIRCUIT FERME.............................................................................................. 16
II.5 APERÇUS SUR QUELQUE ELEMENTS DE CALCULS ..................................... 17
II.5.1 PUISSANCE D’UNE TOUR DE REFROIDISSEMENT ................................. 17
II.5.2 DEBIT DE RECIRCULATION ......................................................................... 17
II.5.3 ECART DE TEMPERATURE .......................................................................... 18
II.5.4 TAUX DE CONCENTRATION ....................................................................... 18
II.5.5 DEBIT DE DECONCENTRATION D .............................................................. 18
II.5.6 DEBIT D’EVAPORATION 𝑬 ........................................................................... 18
II.5.7 ENTRAINEMENT VESICULAIRE 𝑬V ............................................................ 19
II.5.8 FUITES SUR L’INSTALLATION (F) .............................................................. 19
II.5.9 PERTES TOTALE ............................................................................................. 19
II.5.10 APPOINT A ................................................................................................... 19
II.5.11 SURFACE D’ECHANGE 𝒂 .......................................................................... 19
II.6 PARAMETRES DE PERFORMANCE DE LA TOUR DE REFROIDISSEMENT 20
II.7 CONDITIONS POUR EVALUATION DES PERFORMANCES ........................... 21
II.7.1 CONDITIONS MECANIQUES ........................................................................ 21
II.7.2 CONDITIONS ATMOSPHERIQUES .............................................................. 21
II.8 THEORIE DE LA TOUR DE REFROIDISSEMENT .............................................. 22
II.8.1 INTEGRATION A CONTRE-COURANT ....................................................... 25
III. CALCUL DE PERFORMANCE DE LA TOUR DE REFROIDISSEMENT (RA1Z) 27
III.1 INTRODUCTION ................................................................................................. 28
III.2 LA MESURE DE LA TEMPERATURE ET DE L’HUMIDITE RELATIVE ..... 28
III.3 PARAMETRES INFLUENÇANT LA PERFORMANCE DE LA TOUR DE
REFROIDISSEMENT ......................................................................................................... 30
III.4 CONDITIONS POUR L’EVALUATION DES PERFORMANCES ................... 31
III.5 CALCULS DE PERFORMANCE DE LA TOUR DE REFROIDISSEMENT .... 31
III.5.1 PREMIERE METHODE DE CALCUL ............................................................ 31
III.5.2 ETUDE DE L’EFFET DES VARIABLES TEMPERATURES ET HUMIDITE
SUR LES PARAMETRES DE PERFORMANCE .......................................................... 38
III.5.3 DEUXIEME METHODE DE CALCUL ........................................................... 41
III.5.4 CONSERVATION DE L’EAU DANS LA TOUR DE REFROIDISSEMENT 51
III.6 CONCLUSION ...................................................................................................... 56
CONCLUSION GENERAL .................................................................................................... 58
RECOMMANDATIONS ......................................................................................................... 60
BIBLIOGRAPHIE ................................................................................................................... 62
ANNEXE 01 ........................................................................................................................ 64
V
ANNEXE 02 ........................................................................................................................ 65
ANNEXE 03 ........................................................................................................................ 66
ANNEXE 04 ........................................................................................................................ 66
ANNEXE 05 ........................................................................................................................ 67
VI
LISTE DES TABLEAUX
Tableau III.1 : L’humidité relative et les températures de l’eau, de l’air à l’entré et la sortie de
la TR.
Tableau III.4 : Effet de la température d’entrée d’air sec sur le ratio (Qc/G).
Tableau III.10 : Résultat de calcul (effet de l’écart de température sur la quantité d’eau
évaporée.
Tableau III.12 : Résultat de calcul (effet de l’humidité relative sur la quantité d’eau évaporée).
VII
LISTE DES FIGURES
Figure I.1 : Tour de refroidissement 33E101.
Figure III.2 : Effet de la température d’entrée d’air sec sur le ratio (Qc/G).
Figure III.11 : Effet de la température de l’air sec sur la quantité d’eau évaporée.
VIII
LISTE D’ABREVIATIONS
𝜂 Efficacité d’une tour de refroidissement.
Ω Puissance d’une tour de refroidissement.
𝑎 Surface d’échange.
𝐴 Appoint d’eau.
BRI Brut réduit importé.
𝐶 Taux de concentration.
C0, C1, C2 Coefficients d’Antoine.
𝐶𝑠 Chaleur humide.
𝐷 Débit de déconcentration.
𝐸 Débit d’évaporation.
𝐸𝑣 Entrainement vésiculaire.
𝐹 Fuites sur installation.
𝐺 Débit d’air sec.
Ha Enthalpie spécifique de l’air sec.
Hah Enthalpie spécifique de l’air humide.
Hah e Enthalpie spécifique de l’air humide à l’entrée de la tour.
Hah s Enthalpie spécifique de l’air humide à la sortie de la tour.
Hv Enthalpie spécifique de la vapeur d’eau.
𝐻𝑅 Humidité relative de l’air.
HRe Humidité relative de l’air à l’entrée de la tour.
HRs Humidité relative de l’air à la sortie de la tour.
𝑁𝑈𝑇 Nombre d’unités de transfert.
𝑃 Pertes totales.
p Pression atmosphérique.
pv Pression de vapeur.
𝑝sat Pression de vapeur saturante.
𝑇𝑎𝑒 Température de l’air à l’entrée de la tour.
𝑇𝑎𝑠 Température de l’air à la sortie de la tour.
𝑡𝑒 Température de l’eau chaude.
𝑡𝑠 Température de l’eau froide.
𝑡ℎ Température du bulbe humide.
∆𝑇 Ecart de température.
TR Tour de refroidissement.
𝑄𝑐 Débit de circulation d’eau.
we Humidité absolue de l’air à l’entrée de la tour.
ws Humidité absolue de l’air à la sortie de la tour.
yv Fraction molaire de la vapeur.
yvéq Fraction molaire de la vapeur dans l’air saturé.
Introduction générale
Introduction générale
INTRODUCTION GENERALE
La plupart des processus dans l’industrie chimique implique un transfert d’énergie
thermique d’une source à travers le processus et finalement au radiateur. Qu’il s’agisse d’un
procédé de séparation tel que la distillation ou la cristallisation, ou d’un procédé de conversion
de la chaleur en énergie mécanique ou électrique, ou d’une réaction chimique exothermique, il
doit toujours exister un bon moyen d’évacuer la chaleur perdue. L’eau froide est généralement
préférable.
Etant donné que très peu d’usines chimiques ont accès à une quantité d’eau froide
abondante et illimitée, le recyclage de l’eau est une nécessité économique. Il doit donc exister
un moyen de dissiper la chaleur acquise par l’eau ayant traversé le processus et de la renvoyer
à une température où elle est à nouveau utile pour le refroidissement.
Dans le premier chapitre, nous décrirons la raffinerie d’Arzew ainsi que la zone d’utilités 03.
En plus, dans le troisième chapitre nous discuterons les résultats obtenus à base des valeurs
collectés durant notre stage pratique.
1
Chapitre I – présentation de la
raffinerie d’Arzew
2
Chapitre I – présentation de la raffinerie d’Arzew
I.1 INTRODUCTION
La raffinerie d’Arzew est considérée comme un complexe industriel important, réalisé
dans le cadre du premier plan dans les années 70. Elle occupe une superficie de 150 jusqu’à
170 hectares dans la zone industrielle d’Arzew, sur le plateau d’El Mohgoun, à environ 40
kilomètres d’Oran et 5 km de la mer méditerranée. La raffinerie a été construite dans le cadre
du premier plan quinquennal 1970-1973 par un groupe japonais (Société Japonaise JAPAN
GASOLINE CORPORATION). En 1984, d’autres unités ont été rajoutées [1].
Zone 03 (utilités).
Zone 04 (carburants).
Zone 07 (production des huiles de base).
Zone 10 (production des bitumes).
Zone 19 (utilités).
Zone 05 (production des huiles de base).
Le laboratoire contrôle aussi la qualité des eaux, les rejets et procède aux analyses des
différents échantillons spéciaux.
L’eau utilisée au niveau de la zone d’utilités est une eau brute qui provienne de
la station de dessalement KAHRAMA où cette eau sera stockée dans le bac T602 d’une
capacité de 12000 m3, en cas de coupure de cette dernière une réserve est assurée par le bac
T603. Elle comporte :
L’unité 31 comporte 3 chaudières H1, H2, H3 dont deux sont identiques H1 et H2, elles
produisent chacune 62 t/h max de vapeur à une pression de 50 Kg/cm2 à une température
de 430 °C, la troisième chaudière H3 produit une capacité maximale de 124 t/h, elle est prévue
pour remplacer les deux anciennes chaudières. Les chaudières utilisent le fioul gaz comme
combustible [3, 4].
Le débit de circulation est assuré par deux pompes. Cette eau rencontre en tombant le
garnissage, ce qui augmente la surface de contact entre l’eau et l’air aspiré cela permet un
meilleur échange thermique , ensuite l’eau refroidit est refoulée à l’aide de la pompe 33G102
vers un ensemble de trois filtres (33k102 A/B/C) à sable ou seulement deux filtres sont en
fonctionnement , ce filtrage concerne 10% du volume de bassin et il permet d’éliminer la boue
au fond [3], l’eau qui retourne au bassin subit une injection des produits suivants :
Le gaz naturel arrive à la raffinerie à une pression d’environ 28 bars après avoir subi
une détente de 40 bars à 30 bars. Avant son introduction dans le ballon 35D1 le gaz est porté à
une température 40°C dans un réchauffeur à vapeur le 35E1 afin d’éviter son givrage lors de sa
détente jusqu’à une pression de 4 bars. Le gaz en provenance de la zone 4 riche en hydrogène
est introduit directement dans le ballon. La pression dans le ballon 35D1 est maintenue à 4 bars
par les 35PIC 1-2, l’excédent de gaz est envoyer vers torche.
Un autre ballon 35D2 a été installé afin d’alimenter les besoins de la 31GTG01 en gaz
naturel après avoir subi une détente de 30 bars à 25 bars [3].
10
Chapitre II – Théories sur les tours de refroidissement
II.1 INTRODUCTION
L’amélioration du processus de dissipation de l’énergie thermique indésirable est l’une
des tâches les plus importantes des centrales. Dans ce but, les tours de refroidissement sont
utilisées presque partout, elles sont utilisées pour le refroidissement de grandes quantités d’eau
dans les industries de traitement chimique, les centrales thermiques, les centrales nucléaires et
les industries pétrolières. Ces dispositifs de transfert de chaleur et de masse sont basés sur le
refroidissement par évaporation de l’eau en contact avec l’air ambiant. Le volume de travail de
la tour est rempli de matériaux d’emballage pour augmenter le contact entre les deux phases [6,
7].
La théorie des tours de refroidissement a été étudiée en profondeur depuis les premiers
travaux de Merkel en 1925. C’est une description mathématique raisonnablement précise et
relativement simple de la chaleur et des phénomènes de transfert de masse dans une tour à
contre-courant. Jaber et Webb ont présenté un nombre d'efficacité d'unités de transfert méthode
d’analyse particulièrement utile pour les échanges flux de refroidissement des tours. Les deux
méthodes font des hypothèses simplificatrices afin de faciliter le calcul de la hauteur de
remplissage [8, 9].
Les tours de refroidissement à tirage naturel, également connues sous le nom de tours
de refroidissement hyperboliques, sont habituellement privilégiées lors du traitement de grands
volumes d'eau, que l’on retrouve principalement dans les centrales thermiques et les centrales
nucléaires et plus rarement dans l’industrie [11, 12].
Le principe de fonctionnement est le suivant : l’air chaud se déplace vers le haut tandis
que l’air froid est entrainé à l’intérieur de la tour via une entrée d’air située dans le fond.
Les principaux avantages des tours de refroidissement à tirage naturel résident dans le
fait qu'aucun ventilateur n'est utilisé, permettant ainsi d'économiser de l'énergie et de maintenir
un faible niveau sonore. L'absence de pièces en mouvement réduit les opérations de
maintenance et augmente le cycle de vie opérationnel [14].
La circulation d'air dans l'une ou l'autre classe peut être à contre-courant ou à courant
croisé en ce qui concerne l'eau de chute. Le croisement de flux indique que le flux d'air est
horizontal dans la partie remplie de la tour tandis que le contre-courant signifie que le flux d'air
est dans la direction opposée de l'eau de chute.
Toutes les tours mécaniques doivent être localisées de sorte que l'air de décharge se
répande librement sans passer par la tour de recirculation, et ainsi les admissions d'air ne sont
pas restreintes. Les tours de refroidissements devraient être localisées aussi près que possible
des systèmes de réfrigération qu'elles servent, mais ne devraient jamais être localisées au-
dessous d'eux afin de permettre à l'eau du condensateur de sortir hors du système par la tour de
bassin quand le système est arrêté [13].
C’est un système très répondu lorsque l’on dispose d’eau en assez grande quantité. Ce
type de système est de moins en moins utilisé en raison des prélèvements d’eau importants qu’il
nécessite et de l’augmentation de chaleur du milieu récepteur [17].
1000
𝜔 = 3600 𝑄𝐶 × 𝐶𝑝𝑒𝑎𝑢 × ∆𝑡 (II.1)
Avec :
𝑠𝑎𝑙𝑖𝑛𝑖𝑡é 𝑑𝑢 𝑐𝑖𝑟𝑐𝑢𝑖𝑡
𝐶 = 𝑠𝑎𝑙𝑖𝑛𝑖𝑡é 𝑑′𝑎𝑝𝑝𝑜𝑖𝑛𝑡 (II.6)
Habituellement, ce taux de concentration est déterminé par la mesure des chlorures, facilement
dosable qui du fait de leur solubilité sont les plus stables dans le circuit.
𝐸
𝐷 = 𝐶−1 (II.5)
Avec :
𝑄𝐶⁄
𝐸𝑉 = 1000 (II.3)
𝑃 = 𝐸 + 𝐷 + 𝐸𝑉 + 𝐹 (II.4)
II.5.10 APPOINT A
C’est le débit d’eau nécessaire pour compenser celle qui est perdue par évaporation, plus
celle qui est perdue par la purge de déconcentration, par entraînement vésiculaire sur la tour et
par d’autres pertes. Il se mesure généralement par un débitmètre.
𝐴 = 𝐷 + 𝐸𝑉 + 𝐸 + 𝐹 (II.7)
𝑎 = 𝑆⁄𝑉 (II.8)
Ou bien :
𝑣𝑎𝑟𝑖𝑎𝑡𝑖𝑜𝑛
𝜂 = 𝑣𝑎𝑟𝑖𝑎𝑡𝑖𝑜𝑛+𝑎𝑝𝑝𝑟𝑜𝑐ℎ𝑒 (II.10)
La vitesse maximale du vent ne dépassant pas 5 m/s, elle est mesurée à une hauteur de
1,5 à 2 m du sol.
La température du bulbe humide de l’air à l’entrée ne dépassant pas ± 5 °C par rapport
au design.
La quantité de l’eau circulante ne doit pas contenir des matières dissoutes plus
que 5000 ppm, et pas plus de 10 ppm de traces d’hydrocarbures.
Le niveau d’eau refroidie dans le bassin ne doit dépasser pas la hauteur normale.
Le débit d’appoint doit être suffisant pour récupérer toutes les pertes qui sont
rencontrées durant le fonctionnement.
La quantité d’eau circulante ne doit pas être supérieure ou inférieure à 10% par rapport
au design.
La goutte d’eau à la température T est entourée par l’air à la température du bulbe sec
Ta, ayant une enthalpie Ha et une humidité absolue Wa. L’interface est assumée d’être un
film d’air saturé avec une température intermédiaire T", une enthalpie h" et une humidité
absolue W [20].
En prenant une valeur constante de 4,18 kЈ/ (kg.K°) pour la chaleur spécifique de l'eau
Cp, le transfert d'énergie total de l'eau vers l'interface est ;
Ou :
Ou :
Le processus atteint l’équilibre lorsque ta=t, et l'air devient saturé d'humidité à cette
température dans des conditions adiabatiques, l'équilibre est atteint à la température de
saturation adiabatique ou à la température du bulbe humide de l'air.
C'est la température la plus basse possible dans une tour de refroidissement. L'eau en
circulation approche rapidement cette température lors une tour fonctionne sans charge
thermique. Le processus est le même lorsqu'une charge de chaleur est appliquée, mais
l'enthalpie de l'air augmente à mesure qu'elle se déplace à travers la tour afin que la température
d'équilibre augmente progressivement.
Merkel (1925) a supposé que la relation de Lewis était égale à un en combinant le transfert
de masse et de chaleur sensible en un coefficient global basé sur la différence d'enthalpie,
𝐾𝐺
=1 (II.16)
𝐾′ 𝐶𝑝𝑚
L'équation considère le transfert de l'interface vers le flux d'air, mais les conditions inter-
faciales sont indéterminées. Si la résistance du film est négligée et le coefficient K’ ce base sur
la force motrice de l’enthalpie h’ à la température de la gouttelette d’eau t, l’équation (II.16)
devient
𝐾′ 𝑎 𝑉 𝑡 𝐶𝑝
= ∫𝑡 2 ℎ′ −ℎ 𝑑𝑡 (II.19)
𝐿 1 𝑎
𝐾′ 𝑎 𝑉 ℎ 𝑑ℎ
= ∫ℎ 2 ℎ′ −ℎ (II.20)
𝐺 1 𝑎
𝐶𝑝 ∆𝑡
=1 (II.21)
ℎ′ −ℎ𝑎
Les équations ne sont pas autosuffisantes et ne sont pas soumises à une solution
mathématique directe. Elles reflètent l'équilibre de masse et d'énergie à n'importe quel
point dans une tour de refroidissement et elles sont indépendantes du mouvement relatif
des deux courants de fluide. L'intégration de l'équation (III.19) donne la NUT pour un
ensemble donné de conditions [11].
27
Chapitre III – calcul de performance de la tour de refroidissment (RA1Z)
III.1 INTRODUCTION
Les bilans massique et thermique donnent les quantités d’air à mettre en jeu et les surfaces
de transfert de chaleur nécessaire au refroidissement de l’eau. Cela permet un dimensionnement
optimal de l’agent de refroidissement atmosphérique.
Tableau III.1 : L’humidité relative et les températures de l’eau, de l’air à l’entré et la sortie
de la TR.
Paramètres
l'entré
de l'air à
température
sortie
de l'air à la
température
l'entré
de l'eau à
température
sortie
de l'eau à la
température
l'air à l'entré
relative de
humidité
Symboles T ae T as te ts HR
Unités °C °C °C °C %
01-mars 18 22,1 34,8 25,1 54
02-mars 20,4 23,1 37 24,1 51
03-mars 18,7 21,3 36,4 26,2 63
04-mars 20 23,8 36,7 26,4 54
05-mars 22,4 25,6 37 27 54
06-mars 25,4 28,6 37 27 29
07-mars 17 20,2 34,3 24,6 57
08-mars 18,6 22,5 35 25 52
09-mars 16,4 19,3 35 25,4 76
10-mars 17,1 20 36 26,4 76
11-mars 19,5 22,5 37,2 26,8 60
12-mars 18,5 21,4 38 28 66
13-mars 22,4 26 36,4 25,5 47
14-mars 17,1 19,8 35,9 25,6 65
15-mars 19,3 23,2 36,4 25,5 57
16-mars 19,2 23 36,7 26,2 60
17-mars 24,2 28 36,7 26,2 23
18-mars 19,6 21,7 36,5 27 69
19-mars 15,8 18 35,3 25,5 68
20-mars 14,7 16,9 34,9 25 70
21-mars 15,4 17,7 35 25.2 72
22-mars 16,1 18 36.2 26 74
23-mars 16,3 18,1 36,8 26,6 79
24-mars 18,6 20,8 35,7 26,4 72
25-mars 18,4 20,5 34,8 26 78
26-mars 16,9 18,6 34 25,6 86
Les relevés ci-dessus ont été effectués au cours de la période de suivi (mois de mars 2019).
L’objectif de ces mesures consiste à déterminer l’influence des conditions météorologiques sur
le fonctionnement de la tour de refroidissement.
Ces performances sont reliées aux conditions de l’air à l’entrée, en particulier le taux
d’humidité qu’il contient. On peut classer les paramètres influents sur ces performances comme
suit :
La vitesse du vent.
La température du bulbe humide de l’air.
La température de l’air sec.
La température d’eau chaude.
La température d’eau refroidie.
La température du débit d’appoint.
L’écart de température de l’eau entre l’entré et la sortie
La quantité de la chaleur évacuée.
La quantité eau évaporée.
Rapport liquide gaz L/G.
Tae=18.41 °C
HRe=63.16%
Tas=21.36 °C
HRs= 100% car supposant l’air à la sortie de la tour est saturé.
La détermination expérimentale de cette valeur est souillée par l’erreur : par mesure de
sécurité, cette température ne peut pas être mesurée directement à la sortie immédiate de la tour
de refroidissement. Elle est mesurée à une hauteur de 2 mètres de ce point. Il s’agit d’un
abaissement de cette température.
La détermination de cette valeur par calcul passe à travers les étapes suivantes :
𝐶1
ln 𝑝 𝑠𝑎𝑡 = 𝐶0 − (III.1)
𝑇−𝐶2
C0 C1 C2
Valeurs standards 23.1963 3816.44 46.13
du P en Pa
A.N :
𝑝 𝑠𝑎𝑡 = 2093.24 𝑃𝑎
é𝑞
Ou 𝑦𝑣 représente la fraction molaire de la vapeur dans l’air et 𝑦𝑣 représente la fraction molaire
de la vapeur dans l’air saturé
𝑝𝑣 𝐻𝑅 . 𝑝𝑣 𝑠𝑎𝑡
𝑊𝑒 = 0.622 𝑝−𝑝 = 0.622 (III.3)
𝑣 𝑝 – 𝐻𝑅 . 𝑝𝑣 𝑠𝑎𝑡
L’équation (III.3) permet de calculer l’humidité absolue qui est en fonction de l’humidité
relative, la pression atmosphérique 𝑝 et la pression de vapeur saturante𝑝𝑣 .
A.N :
𝐻𝑎 (𝑇) = 𝐶𝑎 𝑇 (III.4)
𝐶𝑎 : Capacité calorifiquede l’air pour des températures entre -20°C et 50°C et égale
1.0006kj/kg.K°.
D’autre part l’enthalpie spécifique de la vapeur d’eau est exprimée par la relation suivante
Avec :
A.N:
∆𝐻 = 22.63 𝑘𝐽/𝑘𝑔
Le débit massique d’air sec est constant, tandis que le débit massique d’air humide augmente
par suite de l’évaporation d’une partie de l’eau.
𝐸 = 𝐺(𝑤𝑠 − 𝑤𝑒 ) (III.9)
𝑄 (ℎ −ℎ𝑒 )
𝐺 = (𝐻 −𝐻 𝑐)−ℎ𝑠 (III.11)
𝑠 𝑒 𝑠 (𝑤𝑠 −𝑤𝑒 )
A.N :
𝐺 = 7.21 . 106 𝑚3 /ℎ
Donc la quantité d’eau évaporée durant le refroidissement d’après l’équation (III.9) est :
𝐸 = 68.02 𝑚3 /ℎ
𝜔 = 𝑄𝑐 𝐶𝑝𝐻2 𝑂 . ∆𝑇 (III.12)
A.N:
𝑐 𝑄
𝐸𝑣 = 1000 (III.13)
A.N :
𝐸𝑣 = 4.5 𝑚3 ⁄ℎ
𝐸
𝐷 = 𝐶−1 (III.14)
Dans notre cas au niveau de la zone 03 le taux de concentration est égal à 2,078.
D’où :
𝐷 = 63.10 𝑚3 ⁄ℎ
𝐴 = 𝐷 + 𝐸𝑣 + 𝐸 (III.15)
A.N:
𝐴 = 135.62 𝑚3 ⁄ℎ
𝑣𝑎𝑟𝑖𝑎𝑡𝑖𝑜𝑛 𝑡𝑒 − 𝑡𝑠
𝜂= =
𝑣𝑎𝑟𝑖𝑎𝑡𝑖𝑜𝑛 + 𝐴𝑝𝑝𝑟𝑜𝑐ℎ𝑒 𝑡𝑒 − 𝑇ℎ
Tels que :
𝑡ℎ est la température humide, donnée par le diagramme d’air humide connaissant la température
d’air sec et son humidité (relative ou absolue).
Dans notre cas pour une température d’air sec est de l’ordre de 18.41 °C et une humidité relative
de 63.16% on a :
𝑡ℎ = 14.31°𝐶
A.N :
𝜂 = 46.09 %
Le tableau suivant (III.3) présente toutes les valeurs trouvées dans la 1re méthode de calcul.
Tableau III.4 : Effet de la température d’entrée d’air sec sur le ratio (Qc/G).
Tae HR we ws H ah e H ah s Rapport
(°C) (%) (kJ/kg) (kJ/kg) Qc/G
8 63.16 0.00416 0.01606 18.50 62.28 1.06
10 63.16 0.00477 0.01606 22.08 62.28 0.97
12 63.16 0.00547 0.01606 25.86 62.28 0.88
14 63.16 0.00625 0.01606 29.87 62.28 0.78
16 63.16 0.00713 0.01606 34.13 62.28 0.68
18 63.16 0.00812 0.01606 38.67 62.28 0.57
20 63.16 0.00922 0.01606 43.52 62.28 0.45
La figure (III.2) montre l’effet de la température d’entrée d’air sec sur le ratio (Qc/G).
Figure III.2 : Effet de la température d’entrée d’air sec sur le ratio (Qc/G).
D’après ce graphe nous n’entrevoyons que l’augmentation de la température de l’air sec
à une influence directe sur la performance de la tour de refroidissement. L’enthalpie de l’air
humide augmente d’une façon considérable, cela diminue la capacité de l’air a absorbé énergie
qui se trouve emmagasiner dans l’eau chaude afin de la refroidir. Par résultante, le débit d’air
augmente d’où la diminution du rapport (Qc/G) qui passe de 1,06 pour une température de 8 °C
à 0,45 pour une température de 20 °C.
Tae HR we ws H ah e H ah s Rapport
(°C) (%) (kJ/kg) (kJ/kg) Qc/G
18.41 20 0.00261 0.01606 25.14 62.28 0.90
18.41 30 0.00393 0.01606 28.48 62.28 0.82
18.41 40 0.00525 0.01606 31.83 62.28 0.74
18.41 50 0.00658 0.01606 35.19 62.28 0.66
18.41 60 0.00791 0.01606 38.57 62.28 0.57
18.41 70 0.00925 0.01606 41.86 62.28 0.49
18.41 80 0.01059 0.01606 45.36 62.28 0.41
18.41 90 0.01194 0.01606 48.78 62.28 0.33
La figure (III.3) illustre l’effet de l’humidité relative de l’air sur le ratio (Qc/G) :
Selon ce graphe, la variation de l’humidité relative a aussi une grande influence sur la
performance de fonctionnement de la tour de refroidissement. Quand l’humidité s’accroit,
l’enthalpie de l’air humide accroisse aussi, ce qui diminue la capacité de l’air a absorbé de
l’énergie à partir de l’eau chaude. Par conséquent, le débit d’air augmente, d’où la diminution
du rapport (Qc/G) qui passe de 0.90 pour une humidité relative de 20% à 0,33 pour une humidité
de 90% à la même température 18.41 °C.
Les itérations commencent du bas de la tour (entrée d'air et également sortie de l'eau
fraîche) et continuent jusqu'au sommet de la tour (entrée eau chaude et sortie de l’air
chaud).Connaissant la température d’air sec et son humidité relative, on calcul l’enthalpie
de l’air humide à l’entrée et sortie de la tour de refroidissement. On calcule ensuite la
différence entre ces deux enthalpies, puis l'inverse de la différence.
𝑄𝑐 𝐶𝑝 ∆𝑡 = 𝐺 ∆𝐻
𝑄𝑐 𝐶𝑝 ∆𝑡
𝐻𝑠 = 𝐻𝑒 +
𝐺
Données :
T ae = 18.41°C
HR % = 63.16
Rapport Qc/G =0.55
Après les calculs on trouve que le NUT = 1.842cela nous a permis de déterminer
l’efficacité de la tour de refroidissement. Pour cette méthode de calcul, nous utilisons la formule
suivante :
𝑄𝑐⁄
𝜂%= 𝐺 . 100
𝑁𝑈𝑇
𝜂 % = 29.85 %
Les figures (III.4) et (III.5) montrent respectivement, l’effet de la température sèche sur le
nombre d’unité de transfert (NUT) et l’efficacité de la TR.
Ces graphes montrent nettement que lorsque la température de l’air sec augmente,
l’efficacité de la tour de refroidissement diminue, contrairement au NUT, l’augmentation de la
Cette évolution affirme qu’une température élevée de l’air sec provoque une augmentation
considérable de l’enthalpie de l’air humide ce qui diminuera par la suite la capacité de l’air à
absorber l’énergie de l’eau chaude.
Les figures (III.6) et (III.7) montrent respectivement, l’effet de l’humidité relative sur le
nombre d’unités de transfert (NUT) et l’efficacité de la TR.
On voie bien que d’après ces graphes, l’augmentation de l’humidité relative provoque
l’augmentation du NUT et cela engendre une diminution de l’efficacité de la TR, le contraire
est correct.
Cela peut-être expliqué par l’élévation d’enthalpie humide à l’entrée de la TR qui a une
influence directe sur la capacité de l’air à refroidir l’eau chaude.
Les figures (III.8) et (III.9) illustrent respectivement, l’effet du rapport (Qc/G) sur le nombre
d’unité de transfert (NUT) et l’efficacité de la TR.
D’après la figure (III.9), on voit bien que la variation de l’efficacité en fonction du ratio
Qc/G admet un extremum à Qc/G = 1,5. Ce rapport est en accord avec les études faites par
certains fabricants et qui recommande un rapport Qc / G compris entre [0,75 et 1,5]. [22]
Nous allons s’intéresser seulement aux pertes dues à l’évaporation, car on ce qui
concerne l’entrainement vésiculaire c’est des pertes involontaire et pour le cas de la purge de
déconcentration, un appoint des eaux saumâtres est assuré pour compenser ces pertes vues que
ces eaux présentent une qualité dans les normes des eaux de refroidissement, mais elles ont
comme inconvénient leur température élevée.
Tableau III.10 : Résultat de calcul (effet de l’écart de température sur la quantité d’eau
évaporée).
∆T (°C) E (m3/h)
7 48.20
8 55.08
9 61.97
10 68.85
11 75.74
12 82.62
13 89.51
14 96.39
15 103.28
La figure (III.10) représente l’effet de l’écart de température sur la quantité d’eau évaporée.
Selon ce graphe on voit bien que lorsque l’écart de température augmente, la quantité
d’eau évaporée augmente elle aussi, captivant l’exemple ou la différence de température est de
7° C, la quantité d’eau évaporée est égale à 48.20 m3/h. Pour un écart de 15 °C, la quantité d’eau
évaporée est de 103.28 m3/h.
Cette évolution peut être expliquée par le fait que la quantité d’eau évaporée est en
fonction de l’état de l’air c’est-à-dire la température d’entrée et le taux d’humidité ainsi que le
débit d’air requit pour le refroidissement de l’eau. Comme dans notre cas nous fixons les
caractéristiques de l’air à l’entrée de la tour, le refroidissement de l’eau sera en fonction du
débit d’air cela veut dire que l’écart de température augmente, le débit d’air nécessaire au
refroidissement de l’eau augmente aussi ce qui engendre une augmentation de la quantité d’eau
évaporée car le refroidissement se fais par évaporation et par convection.
Tae HR we ws H ah e H ah s E
(°C) (%) (kJ/kg) (kJ/kg) (m3/h)
10 63.16 0.00477 0.01606 22.08 62.28 53.91
12 63.16 0.00547 0.01606 25.86 62.28 55.91
14 63.16 0.00625 0.01606 29.87 62.28 58.26
16 63.16 0.00713 0.01606 34.13 62.28 61.16
18 63.16 0.00812 0.01606 38.67 62.28 65
20 63.16 0.00922 0.01606 43.52 62.28 70.63
22 63.16 0.01046 0.01606 48.71 62.28 80
La figure (III.11) présente l’effet de la température de l’air sec sur la quantité d’eau évaporée.
Figure III.11 : Effet de la température de l’air sec sur la quantité d’eau évaporée.
Ce graphe montre nettement que lorsque la température sèche augmente, la quantité
d’eau évaporée augmente, prenant l’exemple du jour où la température de l’air sec est de 10°C,
la quantité d’eau évaporée est égale à 53.91 m3/h. pour une température de l’ordre de 22°C, la
quantité évaporée est égale à 80 m3/h.
Par résultante l’action de la température sèche sur la quantité d’eau évaporée est
observée. Une élévation de cette température provoque une augmentation d’enthalpie de l’air
humide à l’entrée de la tour de refroidissement ce qui provoque plus l’évaporation d’eau par
conséquent une augmentation du débit eau d’appoint.
Tableau III.12 : Résultat de calcul (effet de l’humidité relative sur la quantité d’eau
évaporée).
Tae HR we ws H ah e H ah s E
(°C) (%) (kJ/kg) (kJ/kg) (m3/h)
18.41 20 0.00261 0.01606 25.15 62.28 70.14
18.41 30 0.00393 0.01606 28.48 62.28 69.50
18.41 40 0.00525 0.01606 31.83 62.28 68.71
18.41 50 0.00658 0.01606 35.19 62.28 67.72
18.41 60 0.00791 0.01606 38.57 62.28 66.45
18.41 70 00.00925 0.01606 41.96 62.28 64.75
18.41 80 0.01059 0.01606 45.36 62.28 62.36
18.41 90 0.01194 0.01606 48.78 62.28 58.76
La figure (III.12) montre l’influence de l’humidité relative sur la quantité d’eau évaporée
L’humidité relative procrée une enthalpie de l’air humide élevée à l’entrée de la tour de
refroidissement, par conséquent cela diminue l’échange massique et thermique entre l’air et
l’eau ce qui affaiblie la quantité d’eau évaporer car l’air entre à la tour pré-saturer.
La raffinerie utilise des quantités très importantes d’eau soit 12000 m³/j pour une
capacité de traitement 3 750 000 t/an. Les circuits de refroidissement, gros consommateurs
d’eau (45 à 65%de la consommation totale d’eau), sont de plus en plus souvent alimentés en
eau d’appoint.
Dans cette partie, nous allons étudier la possibilité d’économiser une quantité d’eau
évaporée ce qui permet de réaliser un gain concernant les quantités d’eau utilisée comme
appoint de la tour de refroidissement.
D’après les calculs effectués dans ce chapitre, pour un écart de température de 10°C, la
quantité d’eau évaporée est de l’ordre de 68.85m3/h et pour un écart de 15°C la quantité d’eau
perdue est égale à 103.28 m3/h.
Si on opte pour un écart de température de 10°C, la quantité d’eau économisée sera égale à :
Soit
1 m3 d’eau brute (KAHRAMA) coûte les environs de 86 DA, donc le gain réalisé sera égale à :
Donc une tour de refroidissement qui fonctionne avec une ∆t = 10°C réalise un gain de
25 955 950.68 𝐷𝐴.
III.6 CONCLUSION
Les calculs effectués dans ce chapitre nous ont permet d’avoir une estimation de
la quantité d’eau économisée ainsi que le gain réalisé.
Une tour de refroidissement opérant avec un écart de température de 10°C sera mieux
convenable qu’une tour avec un écart plus grand, car cela nous a permet de réalisé un gain de
25 955 950.68 𝐷𝐴, on utilisant seulement trois aéroréfrigérants au lieu de cinq donc on a
économisé aussi de l’énergie électrique nécessaire pour faire fonctionner les deux autres
aéroréférgérants.
Il ne faut pas oublier que les résultats trouvés dans ce chapitre ont été obtenus sur la base
des conditions météorologiques du mois de mars 2019, c’est-à-dire la quantité d’eau évaporée
va changer selon ces conditions ainsi que la quantité d’eau économisée
58
Conclusion générale
CONCLUSION GENERAL
Dans le présent travail, nous avons étudié la performance de la tour de refroidissement
en faisant appel à des concepts approfondis de psychrométrie. Afin d’octroyer notre travail,
nous avons tenté d’étudier, d’évaluer le circonvient des différents paramètres en distinctive
la température sèche, la température humide, l’humidité relative et absolue ainsi que
l’enthalpie de l’air humide sur le bon fonctionnement d’une tour de refroidissement, cela nous
a permis d’effectuer des bilans énergétiques et massiques.
La première méthode est une application directe des bilans massiques et thermiques. La
température de l'air à la sortie de la tour de refroidissement doit être connue pour pouvoir
déterminer l’enthalpie de l’air humide afin d’établir un bilan énergétique.
La deuxième méthode utilisée est celle développée par Merkel. Il a étudié le transfert
massique et thermique à l'interface du film formé entre l'eau et l'air. Les équations obtenues
sont ensuite intégrées numériquement de la température de l'eau froide jusqu'à la valeur de la
température de l'eau chaude pour déterminer le nombre d’unités de transfert (NUT).
Une évaluation économique a été effectuée dans le but de déterminer la quantité d’eau
qui peut être économisée afin de réduire les pertes dues à l’évaporation de l’eau durant le
processus de refroidissement. Les résultats ont montré qu’il été avantageux de travailler avec
un écart de température de 10°C au lieu d’un écart de température plus grand, car cela nous
permet de réaliser un gain considérable.
RECOMMANDATIONS
La réduction des prélèvements d’eau et de la consommation d’eau est la principale
exigence du système de refroidissement de la centrale électrique. Dans les systèmes de
recirculation par voie humide, l’eau de refroidissement chaude est généralement pompée du
condenseur vers une tour de refroidissement ou la chaleur est dissipée directement dans l’air
ambiant par évaporation de l’eau et chauffage de l’air.
L’échangeur de chaleur air-air intégré à une tour de refroidissement crée une résistance
accrue, due à l’échange de la chaleur, nécessitera une puissance supplémentaire du ventilateur
afin de maintenir le même débit dans la section de refroidissement.
La troisième méthode est plus favorable car les résultats ont montré que dans un modèle de
laboratoire, on pouvait s’attendre à une économie d’eau d’environ 35% de l’évaporation, selon
les conditions de l’air ambiant [23].
Entretenir régulièrement le garnissage ainsi que les diffuseurs d’eau qui ont une
influence primordiale sur le rendement de la tour en augmentant le temps de séjour et le
temps de contact air-eau.
Prétraiter l’eau des circuits de refroidissement afin de prolonger la durée de vie des
installations et des échangeurs thermiques.
Suivre régulièrement la qualité d’eau de refroidissement pour éviter les pertes d’eau et
produits chimiques dues à la purge de déconcentration.
Régulier le débit d’air en employant un moteur de ventilateurs de transmission variable
de vitesse et ajuster également l’inclinaison des pales.
Réutiliser les eaux de rejets parvenant de la zone 27 comme appoint au niveau de la tour
de refroidissement afin de minimiser l’utilisation d’eau d’appoint provenant de
ressources externes à la raffinerie.
Réutiliser les eaux saumâtres comme une eau d’appoint des deux évaporateurs 32B101
et 32B102 vue la bonne qualité qu’elles présentent qui répond aux spécifications de la
norme d’eau bruts dans le but de réduire les dépenses en eau d’appoint.
[2] OTMANE Fatima. SI-TAYEB Atika, valorisation des effluents traités à l’unité 27 comme
eau d’appoint de la tour de refroidissement de la raffinerie d’Arzew, magistère raffinage et
pétrochimie, institut algérien du pétrole Boumerdes, 2014.
[6] J.R Picardo. J.E Variyar, the merkel equation revisited: A novel method to compute the
packed height of cooling tower, article, School of Mechanical and Building Sciences, VIT
University, Vellore 632 014, Tamil Nadu, India.
[7] MEHDI Rahmatia. SEYED Rashid Alavi. Mohammad Reza Tavakoli, Experimental
investigation on performance enhancement of forced draft wet cooling towers with special
emphasis on the role of stage numbers, article, Department of Mechanical Engineering, Isfahan
University of Technology, Isfahan 84156-83111, Iran.
[9] Jaber H, Webb RL. Design of cooling towers by the effectiveness-NTU method. ASME J
Heat Transfer 1989 ; 111:837–43.
62
[12] Iffis Tayeb. Amimeur Silia, étude numérique d’une tour de refroidissement humide à
contre-courant, mémoire master, université A. MIRA-BEJAIA, 2016/2017.
[13] D. Hantz, Les différents procèdes de refroidissement d’eau dans les installations
industrielles et tertiaires, Guide technique, CETEAT, France, 2015.
[17] Rouane Islam, Les eaux de refroidissement industriel cas E.N.I.C.A.B (Wilaya de Biskra)
et impact sur l’environnement, mémoire de master, Université Mohamed Khider, Biskra, juin
2016.
[18] Duvivier. L, (2006). Traitement des eaux de refroidissement, Ed TEC & DOC, Paris.
[21] [Link], Raffinage et Génie chimique Tome 2, PARIS : Edition Technip, 1972.
[23] M. Deziania, Kh. Rahmani, S.J. Mirrezaei Roudaki, [Link], Feasibility study for
reduce water evaporative loss in a power plant cooling tower by using air to Air heat exchanger
with auxiliary Fan, Energy and Mechanical Engineering Department, Abbaspour College of
Engineering, Shahid Beheshti University, P.O. Box.16765-1719, Tehran, Iran.
63
ANNEXES
64
ANNEXE 01
Tableau 01 : Caractéristiques physico-chimiques des eaux d’appoint. [5]
CIRCUIT 1 (zone03) CIRCUIT 2 (zone19)
EAU DE SAUMURE EAU SAUMURE ISSUE
ISSUE D’EAU D’EAU
BARRAGE DESSALEE EAU DE EAU EAU EAU DE
2mois/12 10mois/12 BARRAGE DESSALE DESSALEE BARRAGE
E 2mois/12 10mois/12
pH à 25 6 - 7.2 7 – 7.5 6.8 – 7.8 8.2 – 8.6 4–8 3–8
°C
Conducti 2500 – 3300 250 - 500 1400 - 1800 100 - 200 150 – 550 2500 - 4500
vité
𝝁S/cm
TH °F 65 - 80 10 - 15 45 - 50 6 - 10 4 – 20 20 - 130
TAC °F 2–4 3-5 10 - 15 6 - 10 2max 2-6
Calcium 40 - 50 10 – 15 28 - 35 5 - 10 7 – 20 40 – 70
°F
Chlorure 40 - 55 5 – 10 25 - 32 0-5 1–3 20 - 80
°F
ANNEXE 02
Tableau 02 : Données calculé par SONATRACH et NALCO Champion. [5]
Données C/C SONATRACH Données calcul NALCO
Rapport Volume Débit Rapport Volume Débit
concentration (m3) appoint concentration (m3) appoint
(m3/h) (m3/h)
CIRCUIT 1.5 à 1.8 3000 150 à 200 1.8 3000 150
1
CIRCUIT 2.8 à 3.2 1500 30 à 60 3.0 1500 35
2
65
ANNEXE 03
Tableau 03 : Synthèse des programmes de traitement des circuits 1 et 2. [5]
Eau de saumure issue Eau de saumure issue
d’eau d’eau
dessalée ou eau dessalée de barrage ou eau de
(10 mois/12) barrage (2 mois/12)
Rapport de concentration <8 < 2,5 (CT 2) - < 2,0 (CT 1)
Consigne de régulation du 8,0 – 8,3 7,9 – 8,2
pH
Nalco 3DT 227C 50 ppm en continu (pompe doseuse)
Eau de Javel 12% (non 2x/semaine en injection manuelle choc
fournie)
Acide sulfurique En continu pour régulation En continu pour régulation
(non fourni par NALCO) pH (pompe doseuse) pH (pompe doseuse)
Biocide 7330 1 X/2 semaine en injection manuelle choc
Bio détergent 73550 1 X/semaine en injection manuelle choc
Dispersant 7313 En injection manuelle choc en cas de fuite d’hydrocarbure
(10 chocs/an inclus)
ANNEXE 04
Figure 01 : application numérique diagramme psychrométrie
66
ANNEXE 05
67