L'objectif de cette thèse était de modéliser à partir des expériences de percolation-réactive, les processus hydrodynamique et réactionnel qui gouvernent l'altération des ciments de puits. Différentes expériences ont été réalisées dans des conditions représentatives de celles du stockage du CO2. Des échantillons fracturés ont été utilisés pour injecter une saumure enrichie en CO2 à 60°C et 10MPa et à différentes pressions partielles de CO2. Le débit d'injection variait en fonction des propriétés hydrauliques de l'échantillon exposé. L'injection d'une saumure enrichie en CO2 à débit constant à travers une fracture supposée plane a permis d'étudier les modifications des propriétés hydrodynamiques et ces conséquences sur la géochimie et la microstructure du ciment altérée. L'impact dynamique de l'évolution microstructurale a été mis en évidence. Les expériences conduites sur des échantillons présentant de larges ouvertures, réalisées sur une durée de 5 h, ont montré que la perméabilité était maintenue constante le long de l'expérience. Trois couches d'altération se sont développées consécutivement à la dissolution de la portlandite et la décalcification des CSH. L'altération a entrainé la précipitation des carbonates et de la silice amorphe à proximité de la fracture. Dans le cas d'une expérience longue durée appliquant les mêmes conditions que précédemment on a observé que la croissance de la silice amorphe a entrainé la diminution de la perméabilité. Par ailleurs les expériences effectuées sur des échantillons présentant de faibles ouvertures, ont indiqué que la conversion de la portlandite en calcite conduit au colmatage de la fracture. L'évolution des assemblages de phases conduisant à la formation des carbonates et de la silice amorphe ont été modélisés à partir du code géochimique GEMS. Les mécanismes de diffusion et les processus de mise en place des couches d'altération ont été étudiés à partir d'un modèle analytique et d'un modèle de transport réactif à partir du code géochimique PHREEQC. Mots-clés Percolation-réactive, processus hydrodynamique et thermochimique, altération, ciments de puits, stockage du CO2, modélisation. / Title: “Experimental characterization and modelling of the alteration of fractured cement under CO2 storage conditions.”The main purpose of this thesis was to characterize and to model the hydrodynamic and thermochemical processes leading to the alteration of the wellbore cement materials under borehole conditions. Percolation experiments were performed on fractured cement samples under CO2 storage conditions (60°C and 10MPa). Injection flow rate was dictated by the fracture aperture of each sample. CO2 enriched brine was flowed along the fracture aperture, and permeability changes as well as chemical evolution of major cations were continuously acquired during the experiment time. Reaction paths developed by the alteration of the cement were characterized using microtomography and ESEM images. The experiments conducted using samples presenting large fracture apertures during 5h showed that permeability was maintained constant during the experiment time. Three reacted layers were displaying by the alteration of portlandite and CSH. Long term experiment conducted with large initial fracture aperture showed a decrease of the permeability after 15hours of CO2 exposure. Otherwise, experiments performed on samples presenting narrow apertures indicated the conversion of portlandite and CSH to calcite leading to the permeability reduction and the fracture clogging. Assemblages of phases and chemical changes were modelled using GEMS-PSI speciation code. We studied also using a coupled transport-reactive model the conditions leading to the cement alteration and the formation of associated layers.Key words: Hydrodynamic and thermochemical processes, alteration, wellbore cement, CO2 storage, percolation experiments, numerical modeling.

Cette thèse concerne l'étude des risques associés aux projets de CTSC (Captage, Transport et Stockage de CO2) dont le développement est prévu à l'échelle industrielle. Les projets de CTSC sont des systèmes sociotechniques complexes pour lesquels une approche systémique de management des risques est nécessaire. L'approche doit couvrir les différents aspects du risque pour analyser l'influence de la dynamique des risques sur la dynamique des projets. Une méthodologie systémique de management des risques est proposée. Cette méthodologie est fondée sur les concepts de la pensée systémique, de l'approche STAMP (Systems-Theoretic Accident Model and Processes) développée au sein du Massachusetts Institute of Technology, et de la dynamique des systèmes. L'objectif est de modéliser et d'analyser la structure de contrôle de sécurité impliquée dans un projet de CTSC. La structure de contrôle de sécurité est la structure organisationnelle des parties prenantes (contrôleurs) qui sont responsables de maintenir les contraintes de la sécurité. L'objectif de cette structure de contrôle de sécurité dans cette thèse est d'éviter le retard ou l'échec des projets de CTSC. Cet objectif a été reformulé comme étant la définition et le management des risques majeurs qui pourraient empêcher ou limiter le maintien des contraintes de sécurité. Les risques ont été d'abord identifiés et classés selon huit catégories : Technique, SSE (Santé, Sécurité et Environnement), Politique/Stratégie, Réglementation, Organisationnel/Humain, Financier/Economique, Social et Projet. Les risques majeurs liés aux phases amont ont été extraits et modélisés en utilisant la méthodologie proposée. Les rétroactions affectant la propagation et l'amplification de chaque risque ont été étudiées. Les structures de contrôle de sécurité, le contexte et les risques associés des projets de Barendrecht, de Lacq et de Weyburn ont été analysés. L'application de la méthodologie sur ces trois retours d'expériences permet de proposer un modèle générique de contrôle de sécurité pour les projets de CTSC. L'accent est mis sur le rôle majeur des facteurs endogènes conduisant à l'échec des projets de CTSC. Ce modèle met en évidence les flux d'information et de communication entre les parties prenantes qui conduisent à améliorer leurs modèles mentaux et leurs décisions.

Le stockage géologique du CO2 et/ou la récupération assistée de pétrole par injection de CO2 dans des réservoirs pétroliers, pourraient permettre de limiter le CO2 atmosphérique. Cependant, le CO2 peut être associé à de l’oxygène. Prédire l’évolution des hydrocarbures dans ces conditions, implique d’étudier les mécanismes de l’oxydation. Des expériences d’oxydation et des modélisations cinétiques détaillées ont été réalisées avec des composés modèles purs ou en mélange. La comparaison des résultats expérimentaux et de modélisation a permis la construction d’un mécanisme d’oxydation d’hydrocarbures, et a souligné les paramètres influençant l’apparition d’une auto-inflammation. La bonne cohérence des expérimentations et des modélisations, est prometteuse pour le développement d’un outil de prédiction afin de déterminer la limite d’auto-inflammation ainsi que l’évolution de la composition des hydrocarbures, pour estimer la stabilité d’un système pétrolier en contexte d’injection de CO2 / The geological storage of CO2 (CO2 Capture-Storage – CCS) and the Enhanced Oil Recovery (EOR) by CO2 injection into petroleum reservoirs could limit CO2 atmospheric accumulation. However, CO2 can be associated with oxygen. To predict the hydrocarbon evolution under these conditions involves the study of oxidation mechanisms. Oxidation experiment and kinetic detailed modeling were carried out with pure compounds. The comparison between experimental and modeling results led to the construction of a hydrocarbon oxidation kinetic model and emphasized the parameters leading to auto ignition. The good agreement between our experiments and modeling are promising for the development of a tool predicting the critical temperature leading to auto-ignition and the evolution of hydrocarbon composition, to estimate the stability of a petroleum system in CO2 injection context

Un modèle de transport réactif est proposé pour simuler la réactivité des matériaux à base de ciment en contact avec une saumure saturée en CO2 et/ou le CO2 supercritique (CO2sc) dans les conditions de stockage géologique du CO2. Un code a été développé pour résoudre simultanément le transport et la chimie par une approche globale couplée, compte tenu de l'effet de la température et de la pression. La variabilité des propriétés du CO2sc avec la pression et la température, telles que la solubilité dans l'eau, la densité et la viscosité sont pris en compte. On suppose que tous les processus chimiques sont en équilibre thermodynamique. Les réactions de dissolution et de précipitation de la portlandite (CH) et de calcite (CC) sont décrites par des lois d'action de masse et des seuils de produit d'activité ioniques. Une cinétique de dissolution de CH est introduite pour faciliter la convergence numérique. La définition d'une variable principale permet de capturer la précipitation et la dissolution des phases solides à base de calcium. Une généralisation de la loi d'action de masse est développée et appliquée aux silicates de calcium hydratés (CSH) pour tenir compte de la variation continue (diminution) du rapport Ca/Si au cours de la dissolution des CSH. Les variations de porosité et de la microstructure induites par les réactions de précipitation et de dissolution sont également prises en compte. Le couplage entre le transport et la chimie est modélisé par cinq équations de bilan de masse écrites pour chaque atome (Ca, Si, C, K, Cl), ainsi que par une équation de conservation de la masse totale et celle de la charge électrique. Les lois de Darcy et de Nernst-Planck sont utilisées pour décrire le transport de masse et d'ions. Les propriétés de transport dépendent du degré de saturation et de la porosité. Le modèle est implémenté dans le code de volumes finis, Bil. Les principes de cette méthode et l'approche de modélisation sont discutés et illustrés sur un exemple simple. Ce modèle est en mesure de simuler les processus de carbonatation des matériaux à base de ciment, dans des conditions à la fois saturés et insaturés, dans une large plage de concentration de CO2, de température et de pression. Plusieurs expériences, rapportées dans la littérature, sont simulées en utilisant divers types de conditions aux limites: (i) solutions saturées ou non en CO2 et carbonate de calcium, (ii) gas supercritique de CO2. Les prédictions sont comparées avec les observations expérimentales. Certains phénomènes observés expérimentalement peuvent être également expliqués par le modèle / A reactive transport model is proposed to simulate the reactivity of cement based material in contact with CO2-saturated brine and supercritical CO2 (scCO2) under CO2 geological storage conditions. This code is developed to solve simultaneously transport and chemistry by a global coupled approach, considering the effect of temperature and pressure. The variability of scCO2 properties with pressure and temperature, such as solubility in water, density and viscosity are taken into account. It is assumed that all chemical processes are in thermodynamical equilibrium. Dissolution and precipitation reactions for portlandite (CH) and calcite (CC) are described by mass action laws and threshold of ion activity products in order to account for complete dissolved minerals. A chemical kinetics for the dissolution and precipitation of CH and CC is introduced to facilitate numerical convergence. One properly chosen variable is able to capture the precipitation and dissolution of the relevant phase. A generalization of the mass action law is developed and applied to calcium silicate hydrates (C-S-H) to take into account the continuous variation (decrease) of the Ca/Si ratio during the dissolution reaction of C-S-H. The changes in porosity and microstructure induced by the precipitation and dissolution reactions are also taken into account. Couplings between transport equations and chemical reactions are treated thanks to five mass balance equations written for each atom (Ca, Si, C, K, Cl) as well as one equation for charge balance and one for the total mass. Ion transport is described by using the Nernst-Plank equation as well as advection, while gas and liquid mass flows are governed by advection. Effect of the microstructure and saturation change during carbonation to transport properties is also considered. The model is implemented within a finite-volume code, Bil. Principles of this method and modeling approach are discussed and illustrated with the help of a simple example. This model, with all the efforts above, is able to simulate the carbonation processes for cement based materials, at both saturated and unsaturated conditions, in a wide CO2 concentration, temperature and pressure range. Several sets of experiments, including sandstone-like conditions, limestone-like conditions, supercritical CO2 boundary and unsaturated conditions reported in the literature are simulated. Good predictions are provided by the code when compared with experimental observations. Some experimental observed phenomena are also explained by the model in terms of calcite precipitation front, CH dissolution front, porosity profile, etc

Le stockage géologique du carbone dans les veines de charbon est une solution transitoire pour lutter contre le réchauffement climatique. La faisabilité de ce stockage à un coût abordable reste incertaine, en particulier parce que l'injection de dioxide de carbone dans les veines de charbon est lente. Les projets pilotes ont montré que la perméabilité du réservoir diminue lors de l'injection, suite au gonflement du charbon induit par l'adsorption préférentielle du dioxyde de carbone par rapport au méthane présent naturellement. Ce mémoire de thèse est consacré à l'étude de ce gonflement. Un premier travail théorique a consisté à étendre les équations constitutives de poromécanique classique dans les cas où l'adsorption sur des surfaces ou dans des micropores devient significative. Nous avons montré que le comportement poromécanique du solide ne peut être compris que si la dépendance de l'adsorption en fonction de la déformation du milieu poreux est connue. Le couplage entre adsorption et déformation est peu étudié dans la littérature et difficile à mesurer expérimentalement. Dans ce travail, nous avons utilisé la simulation moléculaire qui permet facilement de contrôler indépendamment la pression du fluide adsorbé et la déformation du milieu poreux. A l'aide de simulations moléculaires d'adsorption dans des systèmes modèles unidimensionnels, nous avons validé les nouvelles équations consitutives. Nous avons montré également que l'adsorption peut dépendre de la déformation de façon complexe et qu'elle est très sensible à la structure des micropores. Les résultats de simulations moléculaires d'adsorption dans un modèle moléculaire réaliste de la matrice organique du charbon nous ont permis de montrer que le gonflement du charbon en présence de fluide peut être expliqué par l'adsorption dans les micropores, mais pas dans les mésopores. Nous avons étudié numériquement le couplage entre adsorption et déformation dans le charbon. Le gonflement estimé en associant les simulations moléculaires d'adsorption aux nouvelles équations constitutives de poromécaniques est en bon accord avec les mesures expérimentales. De même, nous avons simulé l'adsorption de mélanges de dioxide de carbone et de méthane dans le charbon à des températures et pressions représentatives des conditions souterraines. Le résultat de ces simulations a permis d'estimer le gonflement différentiel durant l'injection de carbone pour des veines de charbon à différentes profondeurs. / Pas de résumé en anglais

Les émissions de CO2 dans l'atmosphère sont reconnues comme ayant un effet significatif sur le réchauffement climatique. Le stockage géologique de CO2 est largement considéré comme une approche essentielle pour réduire l'impact de telles émissions sur l'environnement. De plus, injecter du dioxyde de carbone dans les veines de charbon remplies de méthane présent naturellement facilite la récupération de ce méthane, un processus connu sous le nom de récupération assistée du méthane des veines de charbon (ECBM en anglais). Mais le gonflement de la matrice de charbon induite par l'adsorption préférentielle de dioxyde de carbone par rapport au méthane conduit à la fermeture du système de cleats (un ensemble de petites fractures naturelles) du réservoir et donc à une perte d'injectivité. Cette thèse de doctorat est consacrée à l'étude de comment cet injectivité évolue en présence de fluides. Nous dérivons deux modèles poromécaniques à double porosité pour une veine de charbon saturée par un liquide pur. Les équations constitutives obtenues permettent de mieux comprendre et modéliser le lien entre injectivité de la veine de charbon et gonflement du charbon induit par l'adsorption. Pour les deux modèles, on considère l'espace poreux du réservoir comme divisé en les cleats macroporeuses et les pores de la matrice de charbon. Les deux modèles diffèrent dans la manière dont l'adsorption de fluide est prise en compte : le premier modèle est limité à une adsorption surfacique, tandis que le deuxième modèle peut être appliqué à l'adsorption dans un milieu possédant un réseau poreux générique, et donc dans un milieu microporeux comme le charbon, pour lequel l'adsorption se déroule principalement par remplissage de micropores. Le second modèle est calibré sur deux charbons avec des propriétés de sorption et de gonflement différentes. Nous effectuons alors des simulations à différentes échelles (du Volume Élémentaire Représentatif, de l'échantillon de charbon, la veine de charbon). En particulier, nous validons notre modèle sur des données expérimentales de variations de perméabilité de charbon induites par l'adsorption. Nous effectuons aussi des simulations de veines dont le méthane serait produit (un processus connu sous le nom de CBM en anglais) ou de veines sans méthane dans lesquelles du CO2 serait injecté. Nous étudions l'effet de différents paramètres tels que les conditions aux limites, la compressibilité de la matrice de charbon, ou la cinétique de transfert de liquide entre les cleats et la matrice de charbon. Dans une dernière partie, le modèle dérivé est étendu aux cas pour lesquels le charbon est en présence de mélanges fluides binaires tels que les mélanges de méthane et le dioxyde de carbone. Nous calibrons entièrement calibré ce modèle étendu sur des données disponibles obtenues expérimentalement et par simulations moléculaires. Des calculs sont alors effectués à l'échelle d'un Volume Élémentaire Représentatif pour prévoir comment sa porosité et sa perméabilité varient en présence de mélanges fluides de méthane et de dioxyde de carbone

Ce mémoire fait un point sur seize années de recherche dans le domaine du transport réactif, depuis mes travaux de thèse en 1997. Ces recherches s'articulent autour de deux axes équilibrés : d'une part le développement du programme de transport réactif Hytec, d'autre part l'application de ce modèle dans différents domaines. Les deux premières parties du rapport détaillent certains points marquants, pour la plupart publiés ou en cours de publication, de ces travaux de recherche en suivant la double logique développement et applications. La dernière partie propose une réflexion sur des pistes de prolongement de ces recherches. Développement d'outil de transport réactif La première partie, très numérique, détaille le fonctionnement du code Hytec auquel j'ai participé depuis mon arrivée dans le laboratoire. Les équations sous-jacentes du modèle sont décrites. Les méthodes de résolution reposent sur une discrétisation volumes finis sur des éléments de Voronoï pour toute la partie hydrodynamique (écoulement, transport et thermique). Le couplage entre la chimie et le transport est résolu de manière séquentielle itérative. Des développements spécifiques sont ensuite présentés. La prise en compte de la rétroaction de la chimie sur le transport exige des méthodes de couplage spécifiques pour assurer la convergence vers la solution : son effet doit être intégré au coeur du couplage entre chimie et transport. La simulation de milieux à double porosité peut être simulée de manière très élégante en dédoublant les mailles chimiques. L'intégration des gaz a des implications sur la résolution de l'écoulement (et la nécessité de résoudre simultanément les équations en pression et en teneur en gaz), le transport (de la phase gaz indépendamment de la phase aqueuse) et le couplage avec la chimie et l'équilibre eau-gaz. Applications Le code Hytec est utilisé dans de nombreux domaines, dans l'équipe et à l'extérieur no- tamment par les membres du consortium de développement Pôle Géochimie Transport. Deux familles d'applications sur lesquelles j'ai particulièrement travaillé ces dernières an- nées sont détaillées. Le stockage géologique de CO2 est une technologie intéressante pour limiter les émissions de gaz à effet de serre. Avant que la technologie soit déployée, il est indispensable de démontrer son efficacité. Les méthodes détaillées s'attachent à vérifier la durabilité des barrières de stockage : puits et couverture. L'exploitation d'uranium par ISR est une technique utilisée pour les gisements de fort tonnage mais faible teneur. Elle consiste à faire circuler des fluides dans le gisement pour dissoudre et extraire l'élément visé. La simulation par transport réactif est particulière- ment bien adaptée pour décrire les phénomènes. Les travaux détaillés visent à simuler la formation des gisements et leur exploitation.

Le stockage géologique du carbone dans les veines de charbon est une solution transitoire pour lutter contre le réchauffement climatique. La faisabilité de ce stockage à un coût abordable reste incertaine, en particulier parce que l'injection de dioxide de carbone dans les veines de charbon est lente. Les projets pilotes ont montré que la perméabilité du réservoir diminue lors de l'injection, suite au gonflement du charbon induit par l'adsorption préférentielle du dioxyde de carbone par rapport au méthane présent naturellement. Ce mémoire de thèse est consacré à l'étude de ce gonflement. Un premier travail théorique a consisté à étendre les équations constitutives de poromécanique classique dans les cas où l'adsorption sur des surfaces ou dans des micropores devient significative. Nous avons montré que le comportement poromécanique du solide ne peut être compris que si la dépendance de l'adsorption en fonction de la déformation du milieu poreux est connue. Le couplage entre adsorption et déformation est peu étudié dans la littérature et difficile à mesurer expérimentalement. Dans ce travail, nous avons utilisé la simulation moléculaire qui permet facilement de contrôler indépendamment la pression du fluide adsorbé et la déformation du milieu poreux. A l'aide de simulations moléculaires d'adsorption dans des systèmes modèles unidimensionnels, nous avons validé les nouvelles équations consitutives. Nous avons montré également que l'adsorption peut dépendre de la déformation de façon complexe et qu'elle est très sensible à la structure des micropores. Les résultats de simulations moléculaires d'adsorption dans un modèle moléculaire réaliste de la matrice organique du charbon nous ont permis de montrer que le gonflement du charbon en présence de fluide peut être expliqué par l'adsorption dans les micropores, mais pas dans les mésopores. Nous avons étudié numériquement le couplage entre adsorption et déformation dans le charbon. Le gonflement estimé en associant les simulations moléculaires d'adsorption aux nouvelles équations constitutives de poromécaniques est en bon accord avec les mesures expérimentales. De même, nous avons simulé l'adsorption de mélanges de dioxide de carbone et de méthane dans le charbon à des températures et pressions représentatives des conditions souterraines. Le résultat de ces simulations a permis d'estimer le gonflement différentiel durant l'injection de carbone pour des veines de charbon à différentes profondeurs.

La détermination des paramètres d'un milieu poreux, notamment de la porosité, de la perméabilité et des propriétés du fluide saturant est un enjeu important pour des problèmes hydrologiques, pétroliers ou de risques naturels. L'objectif principal de ce travail est d'estimer ces propriétés à partir des ondes sismiques réfléchies.<br />Les théories poro-élastiques (Biot,1956) nécessitent de nombreux paramètres pour caractériser les milieux poreux et conduisent à des propriétés particulières des ondes sismiques (deux ondes de compression, atténuation intrinsèquement définies,...). Ces équations sont résolues pour un milieu poreux stratifié plan saturé par un fluide homogène par une méthode de réflectivité associée à une intégration en nombre d'ondes discrets. Ce programme de simulation est tout d'abord utilisé pour estimer la sensibilité des ondes réfléchies à la localisation et à la concentration du dioxyde de carbone dans le cas d'un stockage dans un aquifère marin profond. La sensibilité de la réponse sismique aux différents paramètres du milieu poreux est ensuite établie de manière plus systématique par le calcul analytique des dérivées de Fréchet des sismogrammes et leur mise en oeuvre numérique. Les applications numériques réalisées indiquent que les paramètres primordiaux à déterminer sont la porosité et la consolidation.<br />Ces opérateurs de sensibilité ont ensuite été intégrés dans un code d'inversion de formes d'ondes complètes (algorithme de Quasi-Newton). Les calculs d'inversion réalisés à partir de données synthétiques indiquent que les distributions de porosité et les paramètres caractérisant le solide et le fluide (densité et modules mécaniques) peuvent être correctement reconstruits lorsque les autres paramètres sont bien déterminés.<br />Cependant, l'inversion de plusieurs paramètres reste un problème difficile du fait des couplages sismiques existant entre eux. Il est cependant possible de résoudre des problèmes complexes en ne considérant qu'un seul paramètre pour le fluide (saturation) et un pour les minéraux (lithologie), ou en effectuant des inversions différentielles pour suivre des variations du fluide.<br />La méthode d'inversion est finalement appliquée à un jeu de données réelles acquis sur le site côtier de Maguelonne dans l'Hérault. Les variations du milieu peuvent être reconstruites en utilisant de l'information a priori venant de forages.

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.