É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 Elsevier
Résumé : International audience
Most peridotite massifs and mantle xenoliths show a wide range of isotopic variations, often involving significant decoupling between Hf, Nd and Sr isotopes. These variations are generally ascribed either to mingling of individual components of contrasted isotopic compositions or to time integration of parent-element enrichment by percolating melts/fluids, superimposed onto previous depletion event(s). However, strong isotopic decoupling may also arise during porous flow as a result of daughter-elements fractionation during solid-liquid interaction. Although porous flow is recognized as an important process in mantle rocks, its effects on mantle isotopic variability have been barely investigated so far. The peridotites of the Lherz massif (French Pyrenees) display a frozen melt percolation front separating highly refractory harzburgites from refertilized lherzolites. Isotopic signatures observed at the melt percolation front show a strong decoupling of Hf from Nd and Sr isotopes that cannot be accounted for by simple mixing involving the harzburgite protolith and the percolating melt. Using one dimensional percolation-diffusion and percolation-reaction modeling, we show that these signatures represent transient isotopic compositions generated by porous flow. These signatures are governed by a few critical parameters such as daughter element concentrations in melt and peridotite, element diffusivity, and efficiency of isotopic homogenization rather than by the chromatographic effect of melt transport and the refertilization reaction. Subtle variations in these parameters may generate significant inter-isotopic decoupling and wide isotopic variations in mantle rocks. (C) 2009 Elsevier B.V. All rights reserved.
ISSN: 0012-821X

hal-00413135
https://hal.archives-ouvertes.fr/hal-00413135
DOI : 10.1016/j.epsl.2008.12.033

É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 Oxford University Press (OUP)
Résumé : International audience
We report new structural, microstructural, petrological, and major- and trace-element data on ultramafic rocks from the Kondyor zoned ultramafic complex in Far-East Russia. The ultramafic rocks are subdivided into three subconcentric lithologies, from core to rim: (1) a metasomatic domain where generally phlogopite-rich dykes pervasively intrude dunite; (2) a main dunite core; (3) a pyroxenite rim. The ultramafic rocks have nearly vertical contacts with the surrounding Archaean basement (gneisses, quartzites and marbles) and hornfelsed Riphean sediments. The hornfelsed sediments show a relatively steep ( 60), outward dipping layering, which rapidly flattens to horizontal away from the inner contact. Although the Riphean sediments define a dome-like structure, the inward, shallow dipping foliation of the dunites indicates a synformal structure. Detailed petro-structural investigations indicate that the Kondyor dunites were deformed by solid-state flow under asthenospheric mantle conditions. The outward textural change from coarse- to fine-grained equigranular dunite and the outward-increasing abundance of subgrains and recrystallized olivine grains suggest dynamic recrystallization while fluid circulation was channelized within the core metasomatic zone, with a decreasing melt fraction from core to rim, and also suggest that solid-state deformation induced grain-size reduction towards the cooling border of the Kondyor massif. Based on their geochemistry, the dunites are interpreted as mantle rocks strongly affected by reaction with melts similar to the JurassicCretaceous Aldan Shield lamproites. Rim pyroxenites were formed by a melt-consuming peritectic reaction, implying the existence of at least a small, conductive thermal gradient around the dunite body while the latter was still at near-solidus temperature conditions. This suggests that the zoned structure of Kondyor was initiated at mantle depths, most probably within the subcontinental lithosphere. Upon cooling, the lamproitic melts were progressively focused in the central part of the massif and drained into vein conduits where they reacted with the wall-rock dunite. Two-dimensional numerical modelling based on finite-differences with a marker-in-cell technique incorporates temperature-dependent rheologies for both molten and non-molten host rocks. The modelling consolidates the structural, petrological and geochemical interpretations, which show that the dunites represent the synformal, flat-lying apex of an asthenospheric mantle diapir, triggered by fluid pressure channelized in the core, which nearly reached the Earths surface. We conclude that translithospheric mantle diapirism is an important mode of mass transfer in the Earth.
ISSN: 0022-3530

hal-00424508
https://hal.archives-ouvertes.fr/hal-00424508
DOI : 10.1093/petrology/egn083

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
3D seismic data located in the Gjallar Ridge (Vøring Basin, offshore Norway) reveals a closely-spaced polygonal fault system affecting more than 800 m of homogeneous mud-dominated Quaternary and Tertiary sequences. As some faults reach the modern seafloor, they represent an active polygonal fault system at present day. Even if the processes remain unclear and are still under debate, it is generally agreed that the initiation of polygonal faults is the result of shallow burial dewatering of fine-grained unconsolidated sediments by volumetric compaction. 3D seismic data are commonly interpreted by propagating horizons automatically and by picking faults manually. However, in the case of polygonal fault intervals, this approach is time consuming due to the huge number of faults and because automatic propagation can be misleading. In this study, we applied a new technique of 3D seismic interpretation based on a sequential stratigraphy analysis, using the new PaleoScan© software (Eliis Company). It allowed us to build a 3D geological model computing more than 300 horizons within the faulted intervals. We then used the coherency attribute, depicting anomalies in the shape of seismic waveform like faults, in order to constrain a possible link between fault distribution and stratigraphic levels. Our approach allows fault throws to be calculated in milliseconds on any polygonal fault plane. The result shows that fault segments have been reactivated by dip-linkage. Distribution of faults depends on mechanical units, intervals characterized by different petrophysical properties, which are independent from lithological and diagenetic changes. According to these results, we propose a model showing the evolution of polygonal fault intervals in which faulting stages are separated by a quiescence phase during burial. A first tier of polygonal faults is initiated at a specific depth, according to the Cam-clay model. Then, following a period of quiescence during which mud-rich sediments continued to accumulate, new fault segments are initiated above the first mechanical unit and within this undeformed interval. New nucleated faults then connect downward to pre-existing underlying polygonal fault system, thus progressively increasing the thickness of the faulted interval.
ISSN: 0025-3227

hal-00772470
https://hal.archives-ouvertes.fr/hal-00772470
DOI : 10.1016/j.margeo.2012.07.016

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Electrical resistivity soundings are used by geophysicists to determine the structure and composition of the Earth's crust and mantle and to explore natural resources (ore, oil, gas, water). Their interpretations in terms of composition and in-situ physical conditions depend mainly on laboratory measurements of electrical conductivity of rocks at simulated crustal conditions of temperature, pressure, saturation and pore pressures. These measurements present a numbers of limitations, in particular, in the case where conductive pore fluids are present, as in the case of deep reservoir conditions, where temperature exceeds 250 A degrees C. Here, we present a new cell capable of measuring electrical conductivity of large saturated samples at confining pressure up to 200 MPa, pore pressure up to 50 MPa, and temperature up to 500 A degrees C. The measurement cell has been developed in a commercial, internally heated, gas pressure apparatus (Paterson press). It is based on the concept of "guard ring" electrode, which is adapted to samples that are jacketed by a very conductive, metallic material. Numerical modeling of the current flow in the electrical cell allowed defining the optimal cell geometry. Calibration tests have been performed on Fontainebleau sandstones saturated with electrolytes of different conductivities, up to 350 A degrees C. The resulting electrical formation factor and temperature dependence of electrical conductivity are in very good agreement with previous studies. This new cell will improve the exploration and exploitation of deep fluid reservoirs, as in unconventional, high enthalpy geothermal fields. In particular, the investigations address possible effects of fluid-rock interactions on electrical resistivity of a reservoir host rock.
ISSN: 0169-3913

hal-00681997
https://hal.archives-ouvertes.fr/hal-00681997
DOI : 10.1007/s11242-011-9846-2

É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 Geological Society
Résumé : Polygonal faults, mainly oriented N50, N110 and N170, are abundant in the upper part of the mud-dominated Kai Formation (upper Miocene-lower Pliocene) of the Voting Basin. A second, less-developed tier of polygonal faults, oriented N20, N80 and N140, exists at the base of the overlying Naust Formation (upper Pliocene-Present). The faults abruptly terminate upward below a thick interval of debris flows. We propose a dynamic model in which: (1) the development of polygonal faults discontinues temporarily as a result of a change in regional sedimentation, leading to inactive polygonal faults; (2) rapid emplacement of debris flows in the late Pleistocene creates a new interval of polygonal faults in the lower part of the Naust Formation immediately beneath the debris flow and some faults penetrate into the underlying Kai Formation; (3) some polygonal faults within the Kai Formation are reactivated and propagated upward into the base of the Naust Formation. The high interconnectivity between faulted layers allows the fluids to reach shallower depths, forming well-expressed pipes and pockmarks on the sea floor. The model of cessation/reactivation of polygonal faults constrains the sealing capacity of sedimentary cover over the reservoirs and helps to reconstruct the fluid migration history through the sedimentary column.
ISSN: 0016-7649

hal-00406645
https://hal.archives-ouvertes.fr/hal-00406645
DOI : 10.1144/0016-76492005-178

Éditeur(s) : HAL CCSD
Résumé : International audience
The purpose of our study is to model a geothermal reservoir. When geothermal reservoir are assumed to be composed of pure water, the transfer of mass and energy is classically described by two balance equations: mass balance equation and the energy balance equation. In addition to those equations, fluid velocity ist classically given by the Darcy law while thermodynamic properties, inferred from theoretical or empirical equations of state, are used to close the mathematical system. Once this system is closed, there exist different solutions. The first one is to solve for pressure and temperature with a variable switch to saturation in the two-phase region (e.g. TOUGH2). The second one is to solve for pressure and enthalpy to increase stability of phase transition between single and two-phase states (e.g. Hydrotherm). We adopted the second option and chose te use a splitting method to get rid of the complexity of coupling equations and a finite volume method for the spatial discretization. Selecting object-oriented languages, we developed a multi-language framework, combining Python, Fortran and a C++ implementation of IAPWS (from the freesteam project) including the supercritical equations, in porous media velocity is given by Darcy law and to close the system physical properties are determined by the IAPWS-IF97 thermodynamic formulation. We resolve the equations in pressure and enthalpy instead of pressure and temperature in order to increase stability and to handle easier the passage from a single-phase to a two-phase system. We solve the system by using a splitting method - to get rid of the complexity of coupling equations - and a finite volume method. We offer some freedom to users thanks to the implementation of several methods like explicit or implicit Euler, Runge-Kutta or BDF2 for time solvers or GMRES and BICGSTAB for the linear solver. We can handle several boundary conditions like no-flow - describing a boundary which can not exchange matter with the exterior - or like a mixed-therm condition - a Dirichlet condition to the pressure and a Dirichlet or an outflow condition to the temperature in order to describe a recharge or a discharge zone - ... We're developing a multi-language framework, combining Python, Fortran and the C++ implementation of IAPWS (from the freesteam project). All these languages are object-oriented. We've applied this simulation model to the dogger in Paris, France, to several one-dimensional systems and a two-dimensional one made by Coumou with the CSMP++ platform. The dogger is a reservoir exploited to produce heat by pumping water at 70 and reinjecting it in the reservoir at 40. In the one-dimensional systems we wanted to observe the process of heat transfer from a higher temperature boundary to a smaller one in a high-energy domain. The last simulation shows the natural convection of water in a fault. For every simulation we compared the solutions we found with an other code (TOUGH2 or CSMP++) and they agreed. The next step will be to model the geothermal plant in Guadeloupe, West Indies. It's the only place in France - and in the West Indies so far - producing electricity with the earth power. The temperature can reach up to 1000 and the pressure range is around a few hundreds MPa. In some surface zones we can see two-phase water at atmospheric pressure. In the 1980s Bouillante was a laboratory for France. Since 1995 Bouillante has given 30GWh electricity a year to the Guadeloupeans.
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering
Stanford, United States

Droits : info:eu-repo/semantics/OpenAccess
hal-00944133
https://hal-brgm.archives-ouvertes.fr/hal-00944133
https://hal-brgm.archives-ouvertes.fr/hal-00944133/document
https://hal-brgm.archives-ouvertes.fr/hal-00944133/file/Copol.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The processes active in the deep crust above an oceanic subduction zone during its evolution have been constrained through a detailed geochemical study (major and trace elements and Sr, Nd and Pb isotopes) of representative samples through an 30 km thick exhumed crustal section of the Cretaceous Kohistan oceanic island arc (Northern Pakistan). The use of both trace elements and radiogenic isotopes reveals two distinct geochemical suites (suites A and B) within the JijalPatanKiruKamila (JPKK) complex. Suite A is characterized by a progressive enrichment in Pb-207 and a decrease in Nd-143/Nd-144 with increasing La-N/Sm-N. Suite B has higher Pb-207/Pb-204 and lower Nd-143/Nd-144 ratios with approximately constant La-N/Sm-N. By combining trace elements with different partitioning behaviour it is demonstrated that there is an increasing contribution of the subduction component in the magmas with time. It is also possible to distinguish a slab component imprint carried by aqueous fluids from one corresponding to sediment melts. Intrusive granites are abundant in the upper levels of the JPKK section. All were generated at the arc root level (Jijal crustal section) during dehydration-melting of hornblende-rich plutonic rocks. A three-stage geodynamic model is proposed for the evolution of the arc over a period of 30 Myr. The first stage (117105 Ma) starts with the onset of subduction, which was followed by the building of the volcanic arc. The second stage (10599 Ma to 9691 Ma) corresponds to a major igneous event, which was characterized by abundant magma underplating and granulite-facies metamorphism at the arc base. Recycling of the residualcumulative lower crust into the convective asthenospheric mantle was efficient during this stage, and was probably related to thermo-mechanical erosion of the base of the crust. The last stage (9585 Ma) corresponds to a period of low magmatic activity, which marked the end of the intra-oceanic subduction. This is related to the formation of a cold blanket above the slab surface as a result of thermo-mechanical erosion of the cold walls of the subduction zone (i.e. the upper part of the slab and the base of the overriding plate), and corner flow dragging the cold material into the zone of melt generation. Ultimately, a voluminous magmatic pulse occurred around 85 Ma (forming the Chilas complex), before arccontinent collision.
ISSN: 0022-3530

hal-00420068
https://hal.archives-ouvertes.fr/hal-00420068
DOI : 10.1093/petrology/egp010

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Random walk (RW) methods are recurring Monte Carlo methods used to modelconvective and diffusive transport in complex heterogeneous media. Many applications canbe found, including fluid mechanic, hydrology and chemical reactors modeling. These methodsare easy to implement, very versatile and flexible enough to become appealing for manyapplications because they generally overlook or deeply simplify the building of explicit complexmeshes required by deterministic methods. RWprovides a good physical understandingof the interactions between the space scales of heterogeneities and the transport phenomenaunder consideration. In addition, they can result in efficient upscaling methods, especiallyin the context of flow and transport in fractured media. In the present study, we review theapplications of RW to several situations that cope with diverse spatial scales and differentinsights into upscaling problems. The advantages and downsides of RW are also discussed,thus providing a few avenues for further works and applications.
ISSN: 0169-3913

insu-01308060
https://hal-insu.archives-ouvertes.fr/insu-01308060
DOI : 10.1007/s11242-016-0693-z

Éditeur(s) : HAL CCSD
Résumé : International audience
The Achmmach tin mineralization (NE of theMoroccan Central Massif) is associated with tourmaline-richalteration halos, veins, and faults hosted in sandstones andmetapelites of the Upper Visean-Namurian. These depositsare reported to be late Hercynian in age and related to theemplacement of late-orogenic granite not outcropping in thestudied area. Structural and paragenetic studies of theAchmmach tin deposit were conducted in order to establisha general model of the mineralization. From field constraints,the late Hercynian phase is marked by a transition fromtranspression to extension with deformation conditions evolvingfrom ductile to brittle environments. The transpression(horizontal shortening direction roughly trending E-W) is coevalwith the emplacement of the first tourmaline halos alongseveral conjugated trends (N070, N020, and N120).Thereafter, a tourmaline-rich breccia formed in response tothe fracturing of early tourmaline-altered rocks.Subsequently, during the extensional phase, these structureswere reactivated as normal faults and breccias, allowing theformation of the main tin mineralization (cassiterite) associatedwith a wide variety of sulfides (arsenopyrite, chalcopyrite,sphalerite, galena, pyrrhotite, bismuthinite, pyrite, and stannite). This evolution ends with fluorite and carbonate deposition.The hydrothermal fluid flow velocity, calculated byapplying statistical measures on the tourmaline growth bands,varies with the lithology. Values are lower in metapelites andhigher in breccia. In the general evolution model proposedhere, tourmaline alteration makes the rock more competent,allowing for brittle fracturing and generation of open spacewhere the main Sn mineralization was precipitated.
Mineralium Deposita

insu-01240481
https://hal-insu.archives-ouvertes.fr/insu-01240481
DOI : 10.1007/s00126-015-0613-0