Influence of the spatial variability of the thermodynamic and chemical kinetics parameters on the precipitation of carbonate minerals in porous media (CO2 mineral storage)

Ce travail entre dans le cadre de l’étude des interactions eau-roche dans le cas du stockage du CO2 en milieu géologique. Un intérêt particulier est accordé aux hétérogénéités des paramètres associés aux phénomènes géochimiques. Ces hétérogénéités peuvent s’observer à différentes échelles: celle des grains (les minéraux présentent des défauts de cristallinité et des impuretés), et l’échelle centimétrique/pluri-décamétrique. En particulier, les paramètres thermodynamiques (logK) et de cinétique chimique (dans ce travail nous avons considéré le produit de la constante cinétique k par la surface spécifique S soit kS comme "paramètre de cinétique chimique") sont connus à partir des expériences de laboratoire pour des échantillons de quelques centimètres de dimension, alors que l’on s’intéresse aux réactions minéralogiques à l’échelle des réservoirs.Nous avons évalué les caractéristiques géostatistiques de la variabilité spatiale après réaction à travers des simulations de transport réactif dans lesquelles différents paramètres (logK et kS) sont perturbés avec une première variabilité imposée. Une combinaison de deux approches est ainsi abordée : déterministe et géostatistique. Le code du transport-réactif COORES (IFP-EN et Ecole nationale supérieure des mines de Saint-Etienne) a été utilisé pour les simulations déterministes et le système géochimique étudié concerne la dissolution du diopside avec précipitation de minéraux secondaires comme la calcite et la magnésite.Après analyse par la méthode des plans d’expériences, les résultats montrent qu’une corrélation spatiale élevée combinée avec une grande variance de dispersion des minéraux favorise une réactivité importante des minéraux lorsqu’on perturbe le paramètre de cinétique chimique kS. Par ailleurs une vitesse d’injection élevée accélère le processus de dissolution du minéral étudié. La variabilité spatiale du paramètre thermodynamique n’a cependant pas d’effet significatif sur les résultats, le système se comporte comme dans le cas homogène. Du point de vue de l’homogénéisation du paramètre kS, on retrouve l’influence de l’historique de dissolution.The present work is based on the study of water-rock interactions in the case of CO2 storage in geological media. Particular attention is devoted to heterogeneities at different observation scales geochemical phenomena. These heterogeneities can be observed at different scales: the grain (mineral crystallinity present defects and impurities), and the centimeter scale / multi- decametric (rocks are heterogeneous at different scales). In particular, the thermodynamic parameters logK and chemical kinetics kS (in this work we considered the product of the rate constant k by the specific surface area S is kS as "chemical kinetics parameter") are known from laboratory experiments to a few centimeters in size, while we are interested in mineralogical reactions across tanks.We propose to evaluate the geostatistical characteristics of the local variability after reaction through simulations of reactive transport on a small scale in which various parameters (logK and kS) are perturbed with a first spatial variability imposed. A combination of both approaches is discussed: deterministic and geostatistical for the study of geochemical problems at different scales. The reactive transport code - COORES (IFP - EN and Ecole nationale supérieure des mines de Saint -Etienne) was used for deterministic simulations and the geochemical system studied concerns the dissolution of diopside with precipitation of secondary minerals such as calcite and magnesite.After analysis by the method of design of experiments, the results show that high spatial correlation variance combined with high dispersion of minerals promotes a high reactivity when minerals chemically disturbing is the kinetic parameter kS. In addition, a high velocity injection accelerates the dissolution of the mineral studied. However, the effect of spatial variability of the thermodynamic parameter, did not significantly affect the results, the system behaves as in the homogeneous case. From the standpoint of homogenizing the parameter kS, include the influence of the history of dissolution