Modélisation cinétique de la réduction de minerais par le monoxyde de carbone utilisés dans le procédé de combustion en boucle chimique

Le procédé de combustion en boucle chimique est un procédé de production d'énergie générant un courant ne contenant que du dioxyde de carbone et de l'eau. Contrairement aux autres procédés de production d'énergie basés sur la combustion, les gaz ne sont pas dilués dans un courant riche en diazote (N2). Pour séparer le CO2 sans la pénalité énergétique de post-séparation, on utilise l'oxygène contenu dans un oxyde métallique (solide), séquentiellement oxydé et réduit dans deux réacteurs interconnectés. La configuration de lit fluidisé s'est graduellement imposée dans la littérature. Dans le réacteur d'oxydation, opéré en colonne montante (riser), le solide est réoxydé en présence d'air, et dans celui de réduction, opéré en régime bullant, on alimente un combustible (solide, liquide ou gazeux) qui servira d'agent réducteur. On récupère de l'énergie essentiellement dans le réacteur d'oxydation (variation d'enthalpie négative), mais certains oxydes sont aussi réduits par le biais d'une réaction exothermique (p. ex., cuivre). La somme des enthalpies d'oxydation et de réduction du solide est égale à la variation d'enthalpie de la combustion du combustible utilisé. Le déploiement du procédé à l'échelle industrielle repose sur l'identification d'un solide (porteur d'oxygène) pouvant transférer une grande quantité d'oxygène sur une base massique (capacité de transfert d'oxygène), ayant une cinétique rapide (minimisation du volume nécessaire pour une conversion donnée), résistant à l'attrition, non toxique, peu cher et chimiquement stable (désactivation lente). Les critères de non-toxicité et de faible coût favorisent l'utilisation de minerais dans le procédé (par opposition aux oxydes de synthèse). Le dimensionnement du procédé repose sur la connaissance des cinétiques d'oxydation et de réduction du solide. Or, la réaction d'oxydation est toujours beaucoup plus rapide (quelques secondes) que sa réduction (quelques minutes). Ainsi, pour le dimensionnement, la détermination de la cinétique de réduction est suffisante. L'ilménite (de formule FeTiO3) s'est graduellement imposée comme porteur d'oxygène d'origine minérale de référence en raison de son abondance et de sa tenue mécanique exemplaire. Malheureusement, celui-ci est sujet à l'agglomération et possède une capacité de transfert d'oxygène limitée. Nous avons préalablement identifié la pyrolusite (de formule MnO2 dans sa forme fraîche) comme porteur d'oxygène potentiel de remplacement pour l'ilménite. La pyrolusite est un minerai amorphe de manganèse, beacoup moins cher que l'ilménite et abondant. Lors de sa calcination, son degré de cristallinité augmente, et on peut identifier la bixbyite (Mn2O3), l'hausmannite (Mn3O4) et la manganosite (MnO) par analyse de diffraction à rayons X. Un des objectifs de cette thèse est d'identifier le meilleur candidat pour le procédé de combustion en boucle chimique en procédant à la comparaison des cinétiques de réduction par le monoxyde de carbone (CO) de l'ilménite et de la pyrolusite. L'ilménite a d'abord été préoxydé à 1000 �C durant 24 h afin de diminuer sa tendance à l'agglomération par frittage. Malheureusement, cette pratique a comme effet de diminuer sa capacité de transfert d'oxygène en raison de l'augmentation de sa cristallinité, et ce, malgré une augmentation de sa capacité de transfert d'oxygène théorique (formation de pseudobrookite ferrique, de formule Fe2TiO5). La capacité de transfert d'oxygène de l'ilménite fraîche est de 1,7 %, en comparaison à 1,2% pour la pseudobrookite ferrique, lorsque réduite durant 2 min avec du monoxyà carbone à 910 �C. Aux mêmes conditions, la pyrolusite transfère 2 %, et cette capacité double si on coalimente de l'hydrogène (sur la base de l'oxygène transféré au CO seulement, c.-à-d. en excluant la quantité d'hydrogène converti en eau). La cinétique de réduction des deux minerais peut être caractérisée de la même façon : soit une réaction de surface entre les molécules de CO et l'oxygène de surface, et une réaction de réoxydation des sites réduits par l'oxygène en provenance du réseau cristallin. Le modèle permet d'expliquer plus de 90 et 95% de la variance dans les données, pour l'ilménite et la pyrolusite, respectivement. Bien entendu, les valeurs numériques des énergies d'activation des deux oxydes diffèrent grandement. Dans le cas de la réaction de surface, elle est de 51 kJ mol−1 pour l'ilménite. Cette valeur est négative dans le cas de la réaction directe de réduction de la pyrolusite, en accord avec la diminution de l'avancement de cette réaction avec la température. La capacité de transfert d'oxygène tend vers une asymptote pour une augmentation de la température. Aussi, la réaction de réduction des sites surfaciques est réversible, c.-à-d. qu'il est possible de réduire le CO2 en CO à une vitesse augmentant avec la température. La caractérisation cinétique de la réduction de l'ilménite par le CO a aussi permis de déterminer que cette réaction est limitée par le transfert de matière externe, soit entre la volume du gaz et la surface de l'oxyde (théorique du film). Dans le cas de la pyrolusite, l'analyse des critères de Mears et de Weisz-Prater a permis de conclure que le régime est contrôlé par la cinétique. La réaction de réoxydation des sites surfaciques réduits par le CO présente une énergie d'activation de 89 kJ mol−1 pour l'ilménite, alors qu'elle est de 55 kJ mol−1 dans le cas de la pyrolusite La réaction moins aisée de réoxydation de la surface dans le cas de l'ilménite peut être expliquée par sa structure cristalline : des couches alternées Ti et Fe, induisant une certaine résistance étant donné que l'oxygène lié au titane n'est pas pas impliqué dans la réaction de réduction. L'effet de la cinétique plus rapide de la réduction de la pyrolusite se reflète sur l'inventaire de solide nécessaire dans le réacteur de réduction : une diminution d'un ordre de grandeur par unité de puissance thermique. En considérant une cinétique moyenne, on trouve des valeurs de 4,3 et 0,41 tMW−1 pour l'ilménite préoxydée et la pyrolusite, respectivement. Il a été question de modéliser la cinétique de ces réactions de réduction à l'aide de l'approche classique mécanistique. Le choix de cette méthode repose sur les nombreuses incongruités retrouvées dans la littérature concernant la cinétique de réduction de l'ilménite par le CO. Des auteurs différents (mais parfois les mêmes auteurs travaillant avec les mêmes données expérimentales) ont caractérisé la cinétique de réduction de l'ilménite par le CO tour à tour avec des modèles morphologiques (ne considérant que l'évolution de la structure du solide) à coeur rétrécissant (shrinking core) en une et trois dimensions (géométrie planaire et sphérique, respectivement), avec des ordres de réaction variant entre 0,67 et 2 pour le gaz, et des énergies d'activation allant du simple au double, révélant du coup la nature de pur ajustement mathématique de cette approche de modélisation. Nous avons montré que presque la moitié des formes canoniques de ces modèles parvenait à caractériser la courbe de conversion du solide, d'où la nécessité d'un modèle plus rigoureux. Malheureusement, la modélisation cinétique selon l'approche mécanistique fait intervenir des paramètres de surface difficilement mesurables. Par exemple, pour exprimer la variation de la concentration de sites oxydés, nous avons besoin de quantifier la densité surfacique de sites actifs (soit des atomes d'oxygène en surface). Nous avons réussi dans un premier temps à justifier la validité de notre approche ayant recours à un paramètre ajusté pour cette quantité. Les valeurs numériques ajustées ont permis de corroborer la diminution dans le temps de la densité surfacique de sites actifs avec le nombre de cycles de réduction. Ce phénomène peut être interprété par la tendance de l'ilménite à l'agglomération. Dans un deuxième temps, nous avons établi une méthode d'estimation de la densité surfacique de sites actifs. On peut interpréter la méthode comme une courbe de titration simulée. En effet, en discriminant entre réaction de surface et de diffusion, il est possible d'estimer la densité surfacique de sites en prenant l'aire sous la courbe de la concentration du pic initial de CO2 formé dans le cas de la réduction par le CO. Nous en avons illustré le principe à l'aide de la cinétique de réduction de l'ilménite : une fois la surface complètement réduite, on observe que la concentration de CO2 diminue jusqu'à zéro, n'étant affecté que par la dispersion axiale dans le réacteur et le transfert de matière. L'hypothèse sous-jacente à la méthode proposée est la supériorité de la vitesse de la réaction surfacique par opposition à celle de diffusion. La méthode respecte également les critères des méthodes de mesures chimiques existantes, soit l'utilisation d'une sonde chimique identique ou similaire à celle du réactif, et des conditions expérimentales de caractérisation de la surface le plus près possibles des conditions où la réaction a lieu. La méthode nécessite par contre un certain comportement chimique afin de fournir une borne supérieure pour la densité surfacique de sites, soit la saturation de la surface par les réactifs. Les profils de gaz obtenus lors de la réduction de la pyrolusite par le CO ont permis de constater une concentration maximale asymptotique de CO2 produite pour des concentrations de CO croissant. Ce faisant, nous avons pu estimer la densité surfacique de sites actifs de la pyrolusite à 40 μmolm−2, une valeur en accord avec les valeurs mesurées publiées d'oxydes de manganèse de synthèse. Nous nous sommes aussi intéressés à des améliorations à apporter au procédé CLC, en particulier à la valorisation des gaz de combustion. Nous avons en effet démontré la possibilité du reformage de gaz de combustion en gaz de synthèse à l'aide de la pyrolusite réduite. En effet, nous avons montré que la pyrolusite soumise à une centaine de cycles d'oxydoréductions sous méthane et oxygène (séquentiellement) à 800 �C devenait apte à réduire le CO2 en CO, de même que la vapeur d'eau en hydrogène. La production spécifique maximale de monoxyde de carbone atteint 5,5 mol kg−1 entre 810 et 960 �C, même pour un temps de contact de seulement 0,01 s. Cette valeur correspond à 100 fois celle des meilleurs catalyseurs de synthèse, les pérovskites à base de Sr, Ce et Fe. De même, la production d'hydrogène a atteint un maximum de 4,8 mol kg−1 à 947 �C, soit le double des valeurs obtenues à l'aide du cycle dit hercynite à l'aide de catalyseurs de CoFe2O4 et Al2O3. Finalement, nous nous sommes intéressés au couplage d'un oxyde de cuivre et d'un oxyde de calcium dans un contexte de production d'hydrogène. Pour pallier les inefficacités du procédé classique d'enlèvement du CO2 des gaz de reformage dans une tour de lavage, nous avons fait appel à un oxyde mixte cuivre-calcium. Cet oxyde mixte, combinant les propriétés du procédé CLC et les cycles à carbonate, permet d'éviter le recours à de l'oxygène pur pour la régénération de l'oxyde de calcium. Dans cette configuration, la chaleur émise lors de la réduction d'un oxyde métallique fournit la chaleur de régénération de l'oxyde de calcium. CanMET nous a fourni un solide de synthèse combinant un oxyde de cuivre et de calcium, additionné d'un agent liant à base de ciment pour une résistance améliorée à la désactivation et à l'attrition. Cet oxyde mixte peut convertir jusqu'à 4,3% en masse de CO2 par masse initiale en carbonate sur la plage de température allant de 425 à 665 �C sous des concentrations de CO2 entre 10 et 20 %. Nous avons également caractérisé sa cinétique de carbonatation à l'aide d'un modèle empirique de premier ordre à l'équilibre avec des énergies d'activation de 220 et 120 kJ mol−1 avec un coefficient de corrélation d'au moins 99 %. De même, à l'aide de mesures répétées de surface spécifique BET avant et après les 15 premiers cycles de carbonatation-régénération, nous avons pu modéliser sa désactivation par frittage à l'aide d'un modèle de puissance d'ordre 0,5 et d'une constante de détérioration de 2,86 min.