Éditeur(s) : HAL CCSD Elsevier Masson
Résumé : Contaminant migration in aquifers is one of the most debated issues in hydrogeology, as most experimental results display large deviations from the standard (asymptotic) Fickian dispersion theories. Multi-scale investigation and high-resolution sensors are required to determine the origin of non-asymptotic dispersion and validate models. For this, a set of multi-scale Single-Well Injection-Withdrawal (SWIW) tracer tests using a new dual-packer probe CoFIS, including a high resolution optical sensor TELog, are presented. When compared to standard techniques such as salinity measurements, it is shown that high-resolution optical measurements allow an improved characterization of the long-lasting non-asymptotic dispersion mechanisms.
ISSN: 1631-0713

hal-00445250
https://hal.archives-ouvertes.fr/hal-00445250
DOI : 10.1016/j.crte.2009.07.012

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Diffusion parameters for HTO, 36Cl−, and 125I− were determined on Upper Toarcian argillite samples from the Tournemire Underground Research Laboratory (Aveyron, France) using the through diffusion technique. The direction of diffusion was parallel to the bedding plane. The purpose of the present study was 3-fold; it was intended (i) to confirm the I− interaction with Upper Toarcian argillite and to verify the effects of initial I− concentration on this affinity, as previously observed by means of radial diffusion experiments, (ii) to highlight any discrepancy between Cl− and I− diffusivity, and (iii) to investigate the effect of an increase of the ionic strength of the solution on the anionic tracers' diffusive behaviour. The results show that the effective diffusion coefficient (De) and diffusion accessible porosity (εa) values obtained with an ionic strength (I.S.) synthetic pore water of 0.01 eq L−1 are: De = 2.35-2.50 × 10−11 m2 s−1 and εa = 12.0-15.0% for HTO, and De = 14.5-15.5 × 10−13 m2 s−1 and εa = 2.5-2.9% for 36Cl−. Because of anionic exclusion effects, anions diffuse slower and exhibit smaller diffusion accessible porosities than HTO, taken as a water tracer. The associated effective diffusion coefficient (De) and rock capacity factor (α) obtained for 125I are: De = 7.00-8.60 × 10−13 m2 s−1 and α = 4.3-7.2%. Such values make it possible to calculate low 125I distribution ratios (0.0057 < RD < 0.0192 mL g−1) which confirm the trend indicating that the 125I rock capacity factor increases with the decrease of the initial I− concentration. Additional through-diffusion experiments were carried out with a higher ionic strength synthetic pore water (I.S. = 0.11 eq L−1). No evolution of HTO diffusion parameters was observed. The anionic tracers' effective diffusion coefficient increased by a factor of two but no clear evolution of their accessible porosity was observed. Such a paradox could be related to the particularly small mean pore size of the Upper Toarcian argillite of Tournemire. The most significant finding of this study is the large discrepancy (factor of two) between the values of the effective diffusion coefficient for 125I and 36Cl. Whatever the ionic strength of the synthetic solution used, 125I exhibited De values two times lower than those of 36Cl. A detailed explanation for this difference cannot be given at present even if a hypothesis based on ion-pairing or on steric-exclusion cannot be excluded. This makes questionable the assumption usually made for quantifying 125I sorption and postulating that 36Cl and 125I would diffuse in the same porosity. In other terms, at Tournemire, 125I sorption could be more pronounced than previously indicated.
ISSN: 0883-2927

hal-00757749
https://hal.archives-ouvertes.fr/hal-00757749
DOI : 10.1016/j.apgeochem.2011.12.017

Éditeur(s) : HAL CCSD
Résumé : Geological storage of CO2 is one of diverse technologies being explored to reduce atmospheric carbon fromindustrial processes (i.e. fossil fuel combustion). One of the specific features of CO2 injection is the possibility ofgeochemical reactions (dissolution – precipitation) between mobile reactive brine (e.g. formation water enrichedin CO2) and the host rock during the spatial and temporal evolution of CO2. That leads to modifications in thepore structure, which in turn change the flow dynamics of the reservoir (e.g. the permeability k). Then, thesesstructural modifications can largely control the injectivity, so that the pressure field in the reservoir and also theCO2 propagation. Accordingly, it is crucial to explore the changes in the reservoir properties (e.g. structural andhydrodynamic) induced during a CO2 injection and specially the relationships between them (e.g. k or reactivesurface-Sr versus porosity-φ , k versus rock heterogeneity), for developing predictive modelling tools of thetransport and reaction processes occurring during a CO2 injection and reliable risk assessment. In the case ofcarbonate rocks, the application of the predictive models of transport and reaction is still challenging, because oftheir high heterogeneity so that the incertitude in the reaction kinetics of carbonate minerals.From this perspective, we realized brine-enriched in CO2 percolation experiments through carbonate rocksamples in thermodynamic conditions expected during CO2 injection in deep reservoirs (T = 100°C et P =12MPa). The permeability changes k(t) is monitored during the experiments and the porosity variationφ (t) iscalculated from chemical analyses of the sampled outlet fluids, using ICP-EAS. The pore structure modificationsare investigated from high resolution X ray micro tomography images acquired from the synchrotron ofGrenoble (ESRF). Depending to the dissolution regime, controlled by the reservoir rock fabric and the chemicalcomposition of the brine (e.g. PCO2), we observed that a modification of pore structure can either improve(atypical result in dissolution context) or impair the value of the permeability k.
Le stockage géologique du CO2 est l’une des diverses technologies étant explorées afin deréduire les émissions de carbone atmosphérique des processus industriels (i.e. combustion del’énergie fossile). L’une des spécifiques caractéristiques de l’injection du CO2 en profondeurreste la possibilité de réactions géochimiques (dissolution-précipitation) entre la saumureréactive mobile (e.g. eau de formation enrichie en CO2) et la roche encaissante durantl’évolution spatiale et temporelle du CO2, conduisant à des modifications dans la structure despores et par conséquent dans les propriétés d’écoulement du réservoir (e.g. la perméabilité k).Donc, ces changements structuraux peuvent largement contrôler l’injectivité, ainsi que lechamp de pression dans le réservoir et aussi la propagation du CO2. Il demeure ainsi cruciald’explorer les changement dans les propriétés de réservoirs (e.g. structurales ethydrodynamiques) induits durant une injection de CO2 et explicitement les relations existantesentre eux (e.g. k ou surface réactive-Sr versus porosité-φ , k versus hétérogénéité de la roche),afin de développer des outils de modélisation prédictive des processus de transport etréactionnels se produisant durant une injection de CO2 et d’évaluer de façon fiable les risques.Dans le cas des réservoirs carbonatés, l’application des modèles prédictifs de transport réactifdemeure toujours un enjeu, car contrainte par la forte hétérogénéité en leur sein ainsi que parl’incertitude dans la cinétique de réactions des minéraux carbonatés dans ce contexte. Danscette optique, nous avons réalisé des expériences de percolation à travers des échantillons deroches carbonatées dans les conditions thermodynamiques de stockage en profondeur (T =100°C et P =12 MPa). L’évolution de la perméabilité k(t) est suivie au cours desexpériences ; et la variation de la porosité φ (t) est calculée à partir des résultats d’analyseschimiques au ICP-AES des fluides de sortie échantillonnés. L’investigation des modificationsapportées à la structure des pores est réalisée par le biais de la Micro-Tomographie hauterésolution à rayon X, acquise au synchrotron de Grenoble (e.g. ESRF).Dépendant du régime de dissolution, contrôlé par la fabrique de la roche réservoir et lacomposition chimique de la saumuré chargée en CO2 (e.g. PCO2 engagée), on a observéqu’une modification de la structure de la roche peut soit améliorer soit détériorer (résultatatypique en contexte de dissolution) la valeur de la perméabilité k.
https://hal.archives-ouvertes.fr/tel-01142136

tel-01142136
https://hal.archives-ouvertes.fr/tel-01142136
https://hal.archives-ouvertes.fr/tel-01142136/document
https://hal.archives-ouvertes.fr/tel-01142136/file/thesemangane2013.pdf

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : International audience
The dissolution of carbonate rocks usually leads to both porosity (ϕ) and permeability (k) increase. We present experimental evidences and physical-based models of positive and anti-correlated dynamics of k and ϕ observed during dissolution experiments of carbonate rocks. We study the way the rate of change of ϕ and k is controlled by the degree of undersaturation of the percolating solution for two different types of carbonate rocks. We document the occurrence of an anti-correlated k−ϕ trend when the flowing fluid (deionized water) has a weak capacity of dissolution. A positive correlation is found when CO2 is added to the deionized water to increase the potential dissolution rate. Detailed analyses of the microstructures of the rock performed by X-ray microtomography reveal that low dissolution rate favors detachment of solid particles and their subsequent accumulation at the pore-throat inlet. Particles are detached from the rock matrix due to the differential dissolution rate of the indurated grains and the microporous cement. We then propose a simple phenomenological model to interpret the effect of the pore-throat clogging by the accumulation of partially dissolved carbonate particles. We conjecture that permeability is controlled by the decrease of the effective hydraulic radius and the increase of the tortuosity due to partial and localized obstruction of the pore network. Conversely, increasing the level of undersaturation of the flowing solution leads to an augmented potential of dissolving most of the transported particles before they reach the throats. In this case, both k and ϕ increase and display power-law correlations.
ISSN: 0169-3913

hal-01174145
https://hal.archives-ouvertes.fr/hal-01174145
DOI : 10.1007/s11242-015-0456-2

Éditeur(s) : HAL CCSD American Chemical Society
Résumé : The distribution of water and ions in nanoporous membranes, in particular close to their surfaces with external reservoirs, is investigated by means of molecular dynamics (MD) simulations using classical polarizable force fields. A sodium iodide (NaI) aqueous solution is considered in uncharged weakly hydrophilic (WH) and highly hydrophilic (HH) nanoporous materials modeled by a single truncated nanopore in contact with two reservoirs. In the inner part of the nanopore, the spatial distribution of water and ions is obtained by coupling semi-infinitesimal longitudinal and semi-infinitesimal radial descriptions. In the HH material, in contrast with the WH material, species are no longer uniformly distributed in the inner nanopore and water molecules are distributed in a nearly frozen axially periodic corona-like structure which leads to the formation of ionic tunnels. In the reservoirs, water molecules and ions attracted by the membrane surface create a necking of the confined solution which acts as a barrier at the entrance of the nanopore. In the WH nanopore subsurface, close to the surface of the membrane, the species distribution is analogous to the interface distribution between an aqueous electrolyte solution and air. In the HH nanopore subsurface, this distribution is modulated by the water corona-like structure. These results provide useful information for predicting properties of nanoporous membranes, in particular the drastic reduction of diffusion coefficients in HH materials, and give a guide to designing synthetic membranes for applications in nanofiltration, etc.
ISSN: 1932-7447

hal-00682155
https://hal.archives-ouvertes.fr/hal-00682155
DOI : 10.1021/jp2078924

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Mobile-immobile mass transfer is widely used to model non-Fickian dispersion in porous media. Nevertheless, the memory function, implemented in the sink/source term of the transport equation to characterize diffusion in the matrix (i.e., the immobile domain), is rarely measured directly. Therefore, the question can be posed as to whether the memory function is just a practical way of increasing the degrees of freedom for fitting tracer test breakthrough curves or whether it actually models the physics of tracer transport. In this paper we first present a technique to measure the memory function of aquifer samples and then compare the results with the memory function fitted from a set of field-scale tracer tests performed in the same aquifer. The memory function is computed by solving the matrix diffusion equation using a random walk approach. The properties that control diffusion (i.e., mobile-immobile interface and immobile domain cluster shapes, porosity, and tortuosity) are investigated by X-ray microtomography. Once the geometry of the matrix clusters is measured, the shape of the memory function is controlled by the value of the porosity at the percolation threshold and of the tortuosity of the diffusion path. These parameters can be evaluated from microtomographic images. The computed memory function compares well with the memory function deduced from the field-scale tracer tests. We conclude that for the reservoir rock studied here, the atypical non-Fickian dispersion measured from the tracer test is well explained by microscale diffusion processes in the immobile domain. A diffusion-controlled mobileimmobile mass transfer model therefore appears to be valid for this specific case.
ISSN: 0043-1397

insu-00373705
https://hal-insu.archives-ouvertes.fr/insu-00373705
https://hal-insu.archives-ouvertes.fr/insu-00373705/document
https://hal-insu.archives-ouvertes.fr/insu-00373705/file/Gouze_et_al-2008-Water_Resources_Research.pdf
DOI : 10.1029/2007WR006690

Éditeur(s) : HAL CCSD Elsevier
Résumé : Among the radionuclides that can migrate from radioactive wastes buried within consolidated argillaceous rocks, I-129 poses a particular problem due to its high mobility in geomedia and its very long half-life of 1.6 10(7) years. Previous studies carried out on such rocks showed the existence of a low iodide retardation regarding the other anion species. The purpose of the present study was two-fold. It was hoped (1) to confirm this retardation by means of radial diffusion experiments and (2) to estimate the influence of the initial iodide concentrations towards the iodide uptake. The results obtained from in-diffusion stage show that (1) no iodide uptake occurs at high initial concentrations and (2) an uptake at low initial concentrations cannot be clearly evidenced. Finally. the most singular behaviour of iodide was highlighted for all the initial iodide concentrations, from the out-diffusion stage. during which some amount of iodide, available for out-diffusion, were not involved in the diffusion process. Although the origin of such phenomenon remains unclear, a kinetic effect for the iodide uptake was proposed, assuming that the period of 55days (corresponding to the in-diffusion stage duration) would be necessary for inducing the iodide uptake. This kinetic effect, already shown elsewhere by (Descostes. M., Blin, V., Bazer-Bachi, F., Meier, P., Grenut, B., Radwan. J., Schlegel, M.L., Buschaert, S.. Coelho, D, Tevissen, E. Diffusion of anionic species in Callovo-Oxfordian argillites and Oxfordian limestones. In press for Appl. Geochem.) on similar materials, needs to be deepened.
ISSN: 0079-1946

hal-00412393
https://hal.archives-ouvertes.fr/hal-00412393
DOI : 10.1016/j.pce.2008.05.020

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This study explores the dynamics of porosity and reactive surface area changes during porous limestone dissolution by CO2-rich water. The Sr and Ca concentrations in both the rock and the outlet solution are used to evaluate the reactive surface area changes of the two rock-forming calcites, i.e. micrite grains and sparite crystals, which have different trace element signatures. The geometric surface area measured with X-ray microtomography decreases slightly whereas the reactive surface area increases continuously with increasing porosity from 20.3 to 30.2%. Surprisingly, changes in reactive surface areas are very different between the two calcites. The reactive surface area changes in the micrite are parabolic while the reactive surface area of sparite increases greatly. The numerical model HYTEC is used to model the change in reactive surface areas during the experiment. Different geometrical models are tested. The model based on spherical-grain dissolution and spherical-pore growth fails to reproduce the experimental results, while the sugar-lump model provides reasonable agreement with the experiment.
ISSN: 0009-2541

hal-00445274
https://hal.archives-ouvertes.fr/hal-00445274
DOI : 10.1016/j.chemgeo.2009.01.032

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
One of the main difficulties encountered when characterizing the hydrodynamic properties of a fractured aquifer is to identify the preferential flow paths within it. Different methods may be applied to determine the variability of the permeability at the borehole scale and to image the structure of the main flow zones between boreholes. In this paper, we compare the information obtained from different measurement techniques performed in a set of three 100 m depth wells (well-to-well spacing: 5­10 m) in a fractured crystalline rock setting. Geophysical logging and borehole-wall imaging are used to identify open and closed fractures intersecting the boreholes and their orientation. The comparison with flowmeter and single packer tests shows that few of the fractures interpreted as open from geophysical logs are significantly transmissive. Cross-borehole connectivity is first investigated from single packer tests with pressure monitoring in adjacent boreholes. To determine fracture zone connectivity, we propose a methodology simply based on the variation with packer depth of the ratio of the drawdown in the observation well and the drawdown in the pumping well. The results are compared to the analysis of cross-borehole flowmeter tests. We show that both methods provide consistent results with a similar level of information on connectivity.
ISSN: 0022-1694

insu-00221380
https://hal-insu.archives-ouvertes.fr/insu-00221380
DOI : 10.1016/j.jhydrol.2007.07.007

Éditeur(s) : HAL CCSD Elsevier
Résumé : Deficient sealing in the well cement plugs and annulus due to mechanical fracturing is an important risk of CO2 leakage from the reservoir to others permeable layers and the surface. Such situation was reproduced at laboratory scale in order to determine the hydro-chemical control on the fracture permeability. Specifically, we investigated the effect of CO2 rich-brine flowing through fractured Portland cements at T = 60 °C and P = 10 MPa and variable flow rates. We showed that carbonation process is dominant and induces permeability decrease and leakage mitigation for extremely low flow rate whereas for high flow rate injection the permeability remains constant due to the precipitation of a low density secondary Si-rich phase which maintains the initial fracture aperture. At intermediate flow rate the hydraulic aperture can increase due to the densification of the material triggered by the net precipitation of low porosity calcite. These results emphasize that more complex behaviors than those considered from batch experiments may take place in the vicinity of flowing fractures. Specifically, the redissolution of the neoformed calcite as well as the development of amorphous phases, both controlled by the CO2-rich brine renewing rate in the fracture may prevent the healing fracture.
ISSN: 1750-5836

hal-00857555
https://hal.archives-ouvertes.fr/hal-00857555
DOI : 10.1016/j.ijggc.2013.04.002

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Random walk (RW) or Continuous Time Random Walk (CTRW) are recurring Monte Carlo methods used to model convective and diffusive transport in complex heterogeneous media. Many applications can be found, including fluid mechanic, hydrology, and chemical reactors modeling. These methods are easy to implement, very versatile and flexible enough to become appealing for many applications because they generally overlook of deeply simplify the building of explicit complex meshes required by deterministic methods. RW and CTRW provide a good physical understanding of the interactions between the space scales of het-erogeneities and the transport phenomena under consideration. In addition, they can result in efficient up-scaling methods, especially in context of flow and transport in fractured media. In the present study, we review the applications of RW or CTRW for several situations coping with various spatial scales, and different insights into up-scaling applications. RW and CTRW advantages and downsides are also discussed, thus providing a few avenues for further works and applications.
ISSN: 0169-3913

hal-01449131
https://hal.archives-ouvertes.fr/hal-01449131
https://hal.archives-ouvertes.fr/hal-01449131/document
https://hal.archives-ouvertes.fr/hal-01449131/file/0025972-02.pdf
DOI : 10.1007/s11242-016-0693-z

Éditeur(s) : HAL CCSD
Résumé : International audience
The porosity (from values of <1% up to more than 50%, with a mean value of 36.5%) and permeability (from very low values to 20.4 D, with a mean value around 800 mD) of 38 plugs, with a diameter 40 mm, have been measured. These plugs were taken approximately every 1.5 m along a 100 m-long core, sampling a series of reefal carbonate platform lithofacies. A non-unique relationship between porosity and permeability is obtained. In addition, the porosity was measured from processing X-ray micro-tomography images with resolutions ranging from 0.42 to 190.0 µm and sampling volumes ranging from less than 1 µm3 to few cm3. Depending on the structure of the pore network, the computed porosity is, in some cases, controlled by the X-ray image resolution and the sampled volume. While high resolution is required to image micro-porous material, large samples with lower resolution images are necessary for identifying distributed vugs and millimetre-scale structures. This study shows that multi-resolution X-ray micro-tomography is an effective tool for characterizing the multi-scale pore structure of carbonate rocks and understanding how it may control key petrophysical parameters such as porosity and permeability.
Geological Society Special Publication

hal-01174157
https://hal.archives-ouvertes.fr/hal-01174157
DOI : 10.1144/SP406.12

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Both flow field heterogeneity and mass transfer between mobile and immobile domains have been studied separately for explaining observed anomalous transport. Here we investigate non-Fickian transport using high-resolution 3-D X-ray microtomographic images of Berea sandstone containing microporous cement with pore size below the setup resolution. Transport is computed for a set of representative elementary volumes and results from advection and diffusion in the resolved macroporosity (mobile domain) and diffusion in the microporous phase (immobile domain) where the effective diffusion coefficient is calculated from the measured local porosity using a phenomenological model that includes a porosity threshold ( inline image) below which diffusion is null and the exponent n that characterizes tortuosity-porosity power-law relationship. We show that both flow field heterogeneity and microporosity trigger anomalous transport. Breakthrough curve (BTC) tailing is positively correlated to microporosity volume and mobile-immobile interface area. The sensitivity analysis showed that the BTC tailing increases with the value of inline image, due to the increase of the diffusion path tortuosity until the volume of the microporosity becomes negligible. Furthermore, increasing the value of n leads to an increase in the standard deviation of the distribution of effective diffusion coefficients, which in turn results in an increase of the BTC tailing. Finally, we propose a continuous time random walk upscaled model where the transition time is the sum of independently distributed random variables characterized by specific distributions. It allows modeling a 1-D equivalent macroscopic transport honoring both the control of the flow field heterogeneity and the multirate mass transfer between mobile and immobile domains.
ISSN: 0043-1397

hal-01278285
https://hal.archives-ouvertes.fr/hal-01278285
DOI : 10.1002/2015WR017645

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : Limestone dissolution by CO2 -rich brine induces critical changes of the pore network geometrical parameters such as the pore size distribution, the connectivity, and the tortuosity which govern the macroscopic transport properties (permeability and dispersivity) that are required to parameterize the models, simulating the injection and the fate of CO2 . A set of four reactive core-flood experiments reproducing underground conditions ( T=100∘C and P=12 MPa) has been conducted for different CO2 partial pressures (0.034
ISSN: 0169-3913

hal-00963597
https://hal.archives-ouvertes.fr/hal-00963597
DOI : 10.1007/s11242-013-0257-4

Éditeur(s) : HAL CCSD
Résumé : The geological storage of CO 2 is considered as one of the solutions to the problem of global warming. The most possible scenarios of geological behavior of soil in the long term conditions has to be precisely defined to avoid catastrophic damages in the future. Chemical properties of saline water in the aquifer after the CO 2 injection becomes more acidic with time. Therefore, chemical degradation of rock occurs, which leads to modifications of its mechanical properties (chemo-mechanical coupling) [1]. In this paper an incremental model of chemo-mechanical coupling using homogenization approach (asymptotic expansion method [2,3]) is proposed. The paper is focused on the numerical computations for the X-ray tomography image of carbonate rock microstructure [4]. The chemical degradation of the microstructure is obtained as a homogeneous dissolution in each direction (Fig.1. – Fig.2.) For every porosity state numerical computations to obtain full, orthotropic elasticity tensor are performed. Therefore, the degradations functions for all components of the rigidity tensor were defined (Fig.3.). Moreover, results of numerical calculations of diffusion, tortuosity and permeability tensors are presented [5,6]. The link between geometrical properties and mechanical and transport parameters are discussed. The numerical analysis of the carbonate rock sample in long term conditions shows the increase of porosity and associated mechanical weakening of the rock. Moreover, the increase of diffusion, tortuosity and permeability tensors elements with time were observed. Finally, the solution of the one-dimensional macroscopic problem of soil settlement due to chemical dissolution of rock at constant load is presented. The method is universal, in the sense that the non-dimensional time scale is used.
ALERT Workshop 2013
Aussois, France

hal-01103175
https://hal.archives-ouvertes.fr/hal-01103175

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Claystone caprocks are often the ultimate seal for CO2 underground storage when residual CO2 gas reaches the reservoir top due to buoyancy. Permeability changes of a fractured claystone due to seepage of CO2-enriched brine and water vapor-saturated CO2 gas are investigated. Results show that brine flow induces a large porosity increase (up to 50%) in the vicinity of the fracture due to dissolution of calcite and quartz, while permeability remains unchanged. Conversely, cyclic flows of CO2-brine and CO2-gas increase the fracture aperture abruptly after each gas flow period, producing a progressive decrease of the caprock seal capacity. Aperture increase is controlled by decohesion of the clay framework within a micrometer-scale-thick layer induced by CO2-gas acidification. Results show that hydraulic aperture increases linearly with duration of the preceding CO2-brine flow period, emphasizing the kinetic control of the quartz grains dissolution during the brine flow periods.
ISSN: 0094-8276

hal-00412332
https://hal.archives-ouvertes.fr/hal-00412332
DOI : 10.1029/2008GL034467

Éditeur(s) : HAL CCSD CRC Press Taylor & Francis Group
Résumé : aucun
Water-Rock Interaction XIII

hal-00857370
https://hal.archives-ouvertes.fr/hal-00857370

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Numerical programs for simulating flow and reactive transport in porous media is essential for predicting reservoir properties related to CO2 sequestration performance, subsurface storage and risk assessment. In this paper we solve the Navier Stokes' equation using finite difference method, on a simulated porous rock structure, to study the velocity distribution of fluid flowing through it under a constant pressure gradient. A reactive solute carried through the fluid is allowed to interact with the minerals in the rock. This chemical reaction dissolves the mineral which changes the rock structure thus affecting its flow properties. These changes of flow properties are studied with variation in reactive solute concentration and pressure field. The different mechanisms of dissolution responsible for the variation of flow properties for the different parameters is predicted. Before the onset of homogeneous dissolution, variation in porosity follows a power-law behaviour with change in permeability when the latter is scaled by the concentration of the reactive species. The simulation results are compared with available experimental data and found to give a reasonable match.
ISSN: 0022-1694

hal-00757760
https://hal.archives-ouvertes.fr/hal-00757760
DOI : 10.1016/j.jhydrol.2012.05.024

Éditeur(s) : HAL CCSD CRC Press Taylor & Francis Group
Résumé : In the chapter : Role of water-rock-gas interaction in sequestration of CO2
Water-Rock Interaction XIII

ISBN : 978-0415604260
hal-00855782
https://hal.archives-ouvertes.fr/hal-00855782

Éditeur(s) : HAL CCSD Elsevier
Résumé : The correct characterization of the effective reactive transport dynamics is an important issue for modeling reactive transport on the Darcy scale, specifically in situations in which reactions are localized, that is when different reactions occur in different portions of the porous medium. Under such conditions the conventional approach of homogenizing only the porous medium chemistry is not appropriate. We consider here reactive transport in a porous medium that is characterized by mass transfer between a mobile and a distribution of immobile regions. Chemical and physical heterogeneities are reflected by distributions of kinetic reaction rate constants and residence times in the immobile zones. We derive an effective reactive transport equation for the mobile solute that is characterized by non-local physical mass transfer and reaction terms. Specifically, chemical heterogeneity is upscaled in terms of a reactive memory function that integrates both chemical and physical heterogeneity. Mass transfer limitations due to physical heterogeneity yield effective kinetic rate coefficients that can be much smaller than the volumetric average of the local scale coefficients. These results help to explain and quantify the often reported discrepancy between observed field reaction rate constants and the ones obtained under well mixed laboratory conditions. Furthermore, these results indicate that transport under physical and chemical heterogeneity cannot be upscaled separately.
ISSN: 0169-7722

hal-00617674
https://hal.archives-ouvertes.fr/hal-00617674
DOI : 10.1016/j.jconhyd.2010.06.002