# Introduction a colonne d'absorption est une unité de séparation physico-chimique largement utilisée en chimie et en biologie. Il s'agit en général d'un tube dans lequel on passe un ou plusieurs mélanges gazeux et qui permet de séparer un ou plusieurs composés du mélange principal. Elle est largement utilisée pour la séparation des gaz acides (CO 2 , H 2 S) du gaz naturel. Le modèle de la colonne d'absorption industrielle présenté dans ce papier est un modèle dynamique qui consiste en un ensemble d'équations non linéaires aux dérivées partielles obtenu à partir de considérations sur les bilans de matière du CO 2 et du MDEA dans les phases gazeuse et liquide [1,2], il prend en compte le gradient thermique le long de la colonne d'absorption. Le bilan énergétique de la colonne industrielle est également posé. Quelques études ont été publiées sur la modélisation et la commande de la colonne d'absorption. Crosby and Durbin [3] étudient les performances d'un régulateur d'état. Roffel [4] développe un régulateur sous optimal avec contraintes d'état. Darwish and Fantin [5] utilisent la commande décentralisée avec placement de pôles. Petrovsky [6] développe un régulateur PI multivariable. Najim [7,8] développe un régulateur auto-ajustable dans le cas de l'absorption du CO 2 par une solution de diéthanolamine. Il reprendra le problème par la suite avec une commande prédictive [9]. Peu d'études ont par contre été publiées concernant la modélisation et la simulation de l'absorption du CO 2 par une solution aqueuses de MDEA ou de MEA en milieu industriel [10][11][12]. Pour le modèle développé dans notre étude, il nous a semblé intéressant de considérer les réseaux de neurones pour la commande de la colonne car elle permet d'obtenir de bonnes performances pour des processus à dynamique complexe. Nous développerons en première étape une régulation PID pour comparer les performances des techniques de régulation classiques avec celles de la commande par réseaux de neurones. # II. # Modelisation et Simulation en Boucle Ouverte de la Colonne Industrielle La colonne d'absorption présentée ici se situe à Khrechba et fait partie du projet In Salah Gaz, elle permet de prélever du CO 2 du gaz naturel en utilisant une solution aqueuse de methydiéthanolamine (MDEA) comme liquide de lavage. Elle est du type colonne à garnissage mesurant 8 mètre de hauteur et 4 m de diamètre. Le garnissage disposé en vrac est du type anneaux de Pall et est destiné à améliorer la surface de contact entre phases. Pour une meilleure élimination du CO 2 du mélange gazeux, le liquide de lavage (eau+ MDEA) circule à contre courant du flux de gaz. La pression et la température de travail sont respectivement de 71.5 bar et 55°C [2]. Lors du contact entre phase liquide et gazeuse sur la surface des anneaux de Pall, le CO 2 passe de la phase gazeuse vers la phase liquide. Cette diffusion est accélérée par réaction chimique du CO 2 avec la MDEA dans la phase liquide. Le débit de MDEA et la concentration du CO 2 en sortie dans le gaz naturel sont respectivement sélectionnés comme variables de commande et de sortie. # Motsclés: colonned'absorption; modélisation; methydiéthanolamine (MDEA); régulation pid; réseaux de neurones, commande non linéaire prédictive neuronale (nnpc). a) Equations Du Modèle Dans le but de simplifier le modèle, nous avons adopté les hypothèses suivantes [2,13]: ? Il n'y a pas de résistance en phase gazeuse Quantité de soluté à l'entrée z = quantité de soluté à la sortie (z+dz) + quantité de soluté transféré de la phase liquide à la phase gazeuse + accumulation. Ce qui donne : (1) où G (m 3 /s) est le débit vol umique du gaz, le flux de CO 2 transféré de la phase gazeuse vers la phase liquide, S la section de la colonne et C Ag (mol/m 3 ) la concentration du CO 2 dans le gaz. Soit UG=G/S (m/s) la vitesse moyenne d'écoulement du gaz, on obtient alors : (2) La réaction chimique entre le CO 2 et le MDEA est [10][11][12]: (3) Où k est la constante de la vitesse de réaction [14,15]: (5) C AL est la concentration du CO 2 dans la phase liquide et C BL la concentration du MEA dans la phase liquide. Le bilan de matière pour le CO 2 dans la phase liquide donne finalement: (6) Ce qui signifie que la totalité du CO 2 transféré dans la phase liquide réagit avec le MDEA. Le bilan de matière pour le MDEA dans la phase liquide donne : La procédure pour calculer le flux ? est donnée en [2] d'après [14][15][16]. En dernier lieu, nous considérons les conditions aux limites qui sont pour la phase gazeuse la concentration du CO 2 en bas de la colonne C Age et pour la phase liquide la concentration du MDEA en haut de la colonne C BLe . (10) La réaction au sein de la colonne industrielle induit un fort dégagement de chaleur et l'apparition d'un La vitesse de réaction r A a la forme suivante [14,15]: dz dt dC S dz S C G C G Ag dz z Ag z Ag + + = + ? ) ( ) ( ? dt dC dz dC U Ag Ag g ? = + ? ? + + ? + + 3 4 2 1 2 3 2 1 2 HCO NCH R R O H NCH R R CO ?? = 2,9610 5 exp( ? 5332 .8 ?? ) ] [ BL AL C kC = ? dz dt dC S dz S C C k C L C L Bl Bl Al dz z Bl z Bl ? ? = + ] [ ) ( ) ( dt dC dz dC U Bl Bl l = ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? = + ? ? t C z C U t C z C U BL BL L Ag Ag g ? ? ? ? ? ? ? = = = = BLe h z BL Age z Ag C C C C 0 © 2015 Global Journals Inc. (US) Global Journal of Researches in Engineering ( ) Volum(4) BL AL A C C k r = gradient de température tout au long de la colonne ; la variation de température est d'environ 5°C entre l'entrée et la sortie de la colonne, ce qui nous amène à établir un bilan énergétique permettant de décrire l'évolution de la température et son effet sur les différentes concentrations le long de la colonne [16] : C i g : la concentration en phase gazeuse à l'interface (mol/m 3 ) C i l : la concentration en phase liquide à l'interface (mol/m 3 ) cp i g : Chaleur spécifique dans la phase gazeuse à l'interface (J/mol.K) h g/l : coefficient de transfert de chaleur (convection)(J/m 2 .K.s) T l : température du liquide (K) T g : température du gaz (K) Î?"H r : l'enthalpie de la réaction (J/mol) cp i l : la chaleur spécifique dans la phase liquide à l'interface (J/mol.K) On prend finalement en compte les conditions aux limites pour la température qui sont les températures respectives du gaz et du liquide en entrée de la colonne. (12) # b) Validation du modèle Un test a pu être effectué sur notre colonne d'absorption industrielle pour comparer la concentration du CO 2 en sortie donnée par le modèle avec celle du système et ce pour une variation de l'entrée du type échelon de 10 t/h. Les données ont pu être recueillies sur un horizon de temps de 6800 secondes. Les résultats sont groupés dans la figure 2 où on représente respectivement, les débits de MDEA, de gaz et les concentrations du CO 2 en sortie de la colonne et celles données par le modèle [2]. On note que la dynamique des concentrations de CO 2 sortie de la colonne obtenue théoriquement agrée avec les résultats expérimentaux. # c) Simulation en boucle ouverte de la colonne industrielle En considérant les équations ( 2) et ( 8), le modèle dynamique de la colonne d'absorption est un système aux paramètres répartis non linéaires. Les résultants de la simulation en boucle ouverte sont présentés en figures 3 et 4. ( Les résultats obtenus en simulation sont satisfaisants, le régulateur PID annule l'erreur en régime permanent et assure une réponse rapide grâce à l'action dérivée. Le système présente en régulation un certain comportement asymétrique dû a la forte non linéarité de la sortie par rapport à la commande choisie qui est le débit du liquide de lavage. ( ) [ ] ? ? ? ? ? ? ? ? ? ? ? = ? ? + ? ? ? ? ? = ? ? ? ? ? ? ? + ? ? ? ? t T T T a.h r Î?"H C cp z T U t T C cp T T a.h z T U l g l l g A r i i l i l l L g i i g i g g l l g g g 1 ? ? ? ? ? ? ? = ? ? = = ? ? = = = = = 0 0 0 0 h z g le h z l z l ge z g z T , T T z T ,) IV. # Commande par Reseaux de Neurones de la Colonne D'absorption Industrielle Les réseaux de neurones artificiels (RNA) sont des approximateurs universels de fonctions et permettent ainsi d'approcher n'importe quelle fonction non linéaire. Cette propriété motive leur utilisation pour la réalisation de commandes non linéaires par identification d'un modèle donné. . Nous envisageons dans notre étude un modèle général pour la synthèse de la commande non linéaire prédictive neuronale (NNPC) [18]- [19]. # Commande non linéaire prédictive neuronale (NNPC) Nous prenons pour la commande NNPC une structure d'identification simple série-parallèle. Le réseau de neurones est choisi en essayant de minimiser l'erreur de rétropropagation du gradient. Après plusieurs essais, notre choix s'est arrêté sur un réseau à quatre entrées u(k),u(k-1),y(k), y(k-1) et une seule couche cachée à 15 neurones. L'apprentissage de notre réseau critère de performance. Le critère J est minimisé en utilisant un algorithme d'optimisation itératif en général l'algorithme BFGS qui est une méthode quasi Newton [18,22]. Le schéma de la figure 9 A cause des erreurs de modélisation on a normalement une erreur en régime permanent. Pour remédier à ce problème, on a introduit un régulateur PI sur lequel on commute dès que la sortie approche le voisinage de la consigne (10%) (figure 13). Toutes fois la réponse du système présente un comportement asymétrique pour les perturbations sur le débit gaz. V. ? ? = = ? + ? ? ? + ? + + ? + = 2 1 1 2 2 ) 2 ( ) 1 ( ) ( ) ( N N j N j d u j t u j t u j t y j t y J ? [ ] T u N t u t u t U ) 1 ( .... .......... ) ( ) ( ? + ? ? = ? # Conclusion Les simulations montrent que nous obtenons des performances satisfaisantes avec les deux techniques de réglage que ce soit pour les échelons de perturbation sur la concentration du CO 2 ou sur le débit de gaz. Il faut noter que la régulation PID donne de meilleurs temps de réponse que la commande neuronale du fait de sa simplicité. 78![où L est le débit volumique du liquide. En tenant compte de (5) et en notant par U L = L/S (m/s) la vitesse moyenne d'écoulement du flux liquide on obtient : (Notre colonne d'absorption est finalement décrite par le système d'équations aux dérivées partielles suivantes: (9)](image-2.png "( 7 ) 8 )") 1![Figure 1 : La colonne d'absorption de In Salah Gas (ISG)](image-3.png "Figure 1 :") 234![Figure 2 : Concentration du CO 2 en sortie du système et du modèle](image-4.png "Figure 2 :Figure 3 :Figure 4 :") 6![Figure 6 : Concentration du CO2 en sortie et débit du liquide de lavage pour un échelon de perturbation de ±10 % sur le débit du gaz](image-5.png "Figure 6 :") ![Concentration du CO2 en sortie (mole %) © 2015 Global Journals Inc. (US) Year 2015 38 de neurones et le test de validation du modèle neuronal sont présentés dans les figures 7 et 8. Simulation Et Régulation Par Réseaux De Neurones D'une Colonne d'Absorption Industrielle](image-6.png "") 58![Figure 5 : Concentration du CO 2 en sortie et débit du liquide de lavage pour un échelon de perturbation de ± 5% sur la concentration de gaz](image-7.png "Figure 5 :Figure 8 :") 910![Figure 9 : Commande NNPC de la colonne d'absorption Industrielle a) Résultats de simulation de la commande NNPC](image-8.png "Figure 9 :Figure 10 :") ![Débit de la MDEA (t/h) Simulation Et Régulation Par Réseaux De Neurones D'une Colonne d'Absorption Industrielle la commande prédictive neuronale. Le contrôleur est composé d'un réseau de neurones pour estimer ? et d'un bloc d'optimisation pour déterminer U'(t).](image-9.png "") 12![Figure 12 : Concentration du CO 2 en sortie et débit du liquide de lavage pour un échelon de perturbation de + /-5% sur la concentration de CO 2 et + /-10% sur le débit de gaz](image-10.png "Figure 12 :") 711![Figure 7 : Identification neuronale de la colonne industrielle](image-11.png "Figure 7 :Figure 11 :") 0.22sortie systeme bruitée[mol%]0.2sortie modeleCgs0.180.1601000200030004000500060007000erreur de prediction0.01erreur [mol%]-0.005 0 0.005-0.0101000200030004000500060007000temps [sec]Year 2015370.2520.25TTdu CO2 en sortie (mole %)0.242 0.244 0.246 0.248Concentration0.238 0.240.2360.2340100200300400500600700800900temps (secondes)0.270.268Concentration du CO2 en sortie (mole %)0.254 0.256 0.258 0.26 0.262 0.264 0.2660.2520.250100200300400500600700800900temps (secondes)debit MDEA [ton/h]1000 1020 104001000200030004000500060007000temps[sec]Gaz [ton/h]415 420debit41001000200030004000500060007000temps[sec][mol%]0.2Cgs0.1501000200030004000500060007000temps[sec]© 20 15 Global Journals Inc. (US) 0.2650.35Concentration du CO2 en (%)0.245 0.25 0.255 0.26Réf Rég Nnpc+PIConcentration du CO2 en (%)0.2 0.25 0.3Réf Rég Nnpc+PI0.24050100150200250300350400050100150200250300350400temps (secondes)temps (secondes)10201050Débit de la MDEA (t/h)1014 1016 1018950 10001012050100150200250300350400900050100150200250300350400temps (secondes)temps (secondes) © 20 15 Global Journals Inc. (US) © 2015 Global Journals Inc. (US) Year 2015 * Design and application of state variable feedback controller for a packed trickle bed column with carbon dioxide absorption JECrosby LDDurbin Joint Automatic Control Conference West Lafayette, IN, USA 1976 * Control of processes with a state variable inequality constraint B Automatica 17 1976 * Stabilization and control of absorber tower chemical process MDarwish JFantin Third IFAC/IFIP/IFORS Conference on System Approach for Development Rabat 1980 * A design scheme for multivariable output constrained PI regulators with applications to absorption column DPetrovski Automatic Control Theory Applications 1980 8 * Modelling and self-adjusting control of an absorption column KNajim International Newspaper of Adaptive Control and Processing Signal 1991 5 * Multilevel learning control of an absorption column KNajim MChtourou Optimal control applications & Methods 12 1991 * Long range predictive control of an absorption packed column KNajim VRuiz Applied References Références Referencias * Simulation of CO 2 absorption in aqueous solution of monoethanolamine (MEA) on a pilot absorption column ASelatnia SSeddari Distillation and absorption conference London, England IChemE 2006 * Identification multimodèle et commande par logique floue d'une colonne d'absorption industrielle. Projet de fin d'études IZidane HZebiri ENP 2010 Alger * Mathematical Modelling Journal 19 1995 * Pilot Plant Studies of the CO2 Capture Performances of Aqueous MEA and Mixed MEA/MDEA Solvents at the University of Regina CO 2 Capture Technology Development Plant and the Boundary Dam CO 2 Capture Demonstration Plant RIdem MWilson PTontiwachwuthikul AChakma Ind. Eng. Chem. Res 45 2006 * CO 2 Absorption by NaOH, MEA and AMP solutions to packed columns PTontiwatchwuthikul AMeisen CJLim Chemical Engineering Science 47 2 1992 * Simulation of pilot plant and industrial CO 2 -MEA absorbers TPintola PTontiwatchwuthikul AMeisen Gas separation and purification 1993 7 * Simulation d'une colonne à garnissage d'absorption du CO 2 par une solution aqueuse de monoéthanolamine à 25°C et 40 °C SSeddari ENP 2004 Alger Thèse de Magister * Rigorous simulation and design of columns for gas absorption and chemical reactions packed columns LDeleye GFFroment Comp. Chem. Eng 10 1986 * Model based control of absorption tower for CO 2 capturing ABedelbayev TGreer BLie 2008 11 Porsgrunn, Norvège Telemark University College * Réseaux de neurones et traitement du signal JHerault CJutien 1994 Paris, Edition Hermès * « Identification et commande par réseaux de neurones d'une colonne d'absorption MSelatnia ENP 2010 Alger Thèse de Magistère * Identification and Control of Dynamical Systems using Neural Networks KSNarendra KParthasarathy IEEE Trans. Neural Networks 1 1 1990 * «Le perceptron Multicouche et Son Algorithme de Rétropropagation des Erreurs PMarc 2004 Université de Laval * Data Sheet ISG. CO2 absorber.C02 control . Ain Salah Gaz 22. Matlab 6 2008 7