Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The mantle above a subducting slab is the site of complex interactions between deformation, partial melting, fluid migration and magma transport. To constrain these interactions and their effects on olivine deformation, we analyze microstructures, crystal preferred orientations, and water contents of peridotite xenoliths entrained by andesites of the Avacha volcano, southern Kamchatka arc. These xenoliths are refractory spinel harzburgites that have coarse-grained microstructures with widely spaced subgrain boundaries and sinuous grain boundaries in olivine, consistent with deformation by dislocation creep under low deviatoric stress (< 13 MPa) and with a significant contribution from diffusional processes. Analysis of crystal preferred orientations (CPO) indicates dominant activation of high-temperature, low-stress {0kl}[100] slip systems in olivine and of(100)[001] in orthopyroxene. In most samples, coarse opx crystals, elongated parallel to the lineation, enclose small olivine grains in crystallographic continuity with neighbouring crystals, indicating crystallization of orthopyroxene at the expense of olivine as a result of reactive percolation of Si-rich fluids coeval with the high-temperature deformation. Secondary crystallization of interstitial orthopyroxene led locally to development of opx-rich lenses parallel to the foliation, characterized by a decrease in olivine grain size and dispersion of the olivine CPO without changing the dominant slip systems. Half of the samples also show acicular orthopyroxene aggregates with which is associated a fine-grained matrix composed of rounded strain-free olivine, orthopyroxene, spinel, and rare amphibole crystals. This matrix occurs pervasively along grain boundaries or forms millimeter-scale irregular lenses and anastomosing veinlets that crosscut the coarse crystals and their ductile deformation structures. Both acicular orthopyroxene and the fine-grained matrix are interpreted as resulting from reactive transport of H2O-rich fluids under static conditions, probably in the lithospheric mantle. Infrared analyses show that olivine contains 1-8 center dot 6 ppm by weight of water. These low water contents are similar to those observed in spinel peridotites from other subduction zones and probably record both the low solubility of water in olivine at low pressure and dehydration during exhumation of the xenoliths. Water contents in orthopyroxene are highly variable (25-506 ppm H2O), probably recording spatially heterogeneous interaction with fluids or melts and compositional disequilibrium in the studied samples. Change in the dominant percolation mechanism from porous flow to fracturing suggests cooling, consistent with the low temperatures estimated from pyroxene thermometry (< 800-900 degrees C). The Avacha xenoliths therefore record pervasive deformation of a region of the mantle under asthenospheric conditions, followed by its accretion to the base of the lithosphere, probably as a result of cooling of the mantle wedge. Percolation of Si-rich fluids or hydrous melts is recorded at all stages; this probably enhanced diffusion and lowered deviatoric stresses during ductile deformation, but did not change the dominant slip direction in olivine from [100] to [001].
ISSN: 0022-3530

hal-00456053
https://hal.archives-ouvertes.fr/hal-00456053
DOI : 10.1093/petrology/egp085

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Production rates in horizontal shale gas wells display declines that are controlled by the low permeability and the intrinsic heterogeneity of the shale matrix. We present an original multi-continuum approach that yields a physical model able to reproduce the complexity of the decreasing gas rates. The model describes the dynamics of gas rate as function of the physical reservoir parameters and geometry, while the shale matrix heterogeneity is accounted for by a stochastic description of transmissivity field. From the 3D (Dimensional) problem setting, including the heterogeneous shale matrix, the fractures generated by the hydrofracking operations, as well as the production well characteristics, we establish an effective upscaled 1D model for the gas pressures in fracture and matrix as well as the volumetric flux. We analyse the decline curves behaviour, and we identify the time scales that characterize the dynamics of the gas rate decline using explicit analytical Laplace space solutions of the upscaled process model. Asymptotically, the flux curves decrease exponentially, while in an intermediate regime we find a power-law behaviour, in which the flux scales with a power law in time as (t^{-beta }), where (beta ) reflects the medium heterogeneity. We use this solution to fit a set of real data displaying distinctly different decline trends and study the sensitivity of the model to the reservoir parameters in order to identify their respective controls at the different stages of the decline curve dynamics. Results indicate that the initial value of the gas rate is determined by the transmissivity of the fractures and the initial pressure of the gas in the shale matrix. The latter causes mainly a shift in the entire decline curve. The early time of decline curve shape is primarily controlled by the fracture properties (compressibility and transmissivity). During the main part of the economically valuable production times, i.e. before the production rate drops exponentially, the decline curve is strongly controlled by the properties of the shale rocks including their heterogeneity, which is modelled by two parameters describing the non-Fickian pressure diffusion effects in a stochastic framework.
ISSN: 0169-3913

hal-01217247
https://hal.archives-ouvertes.fr/hal-01217247
DOI : 10.1007/s11242-015-0504-y

É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é : International audience
This study addresses the mechanical and chemical behavior of allanite during shear zone formation under high-pressure metamorphism. Understanding physico-chemical processes related to the retention or resetting of Pb isotopes in allanite during geological processes is essential for robust petrochronology. Dating of allanite in meta-granodiorite showing variable amounts of strain (from an undeformed protolith to mylonite) at Monte Mucrone (Sesia Zone, NW Italy) gave surprising results. Based on structural and petrographic observations the shear zones at Mt Mucrone are Alpine, yet allanite located within an eclogite facies mylonite yielded Permian ages (208Pb/232Th average age: 287 ± 7 Ma). These mm-sized allanite grains are rimmed by an aggregate of coarse-grained garnet + phengite, thought to derive from former epidote. These aggregates were immersed in a weak matrix that experienced granular flow, and they were thus chemically and mechanically shielded during Alpine mylonitization. In undeformed samples (8a and 8b), two populations of epidote group minerals were found. Allanite forms either coronas around Permian monazite or individual grains with patchy zoning. Both types yield Permian ages (208Pb/232Th age: 291 ± 5 Ma). On the other hand, grains of REE-rich clinozoisite of Cretaceous age are found in undeformed rocks. These grains appear as small fragments with embayed surface outlines and minute satellites or rims around Permian allanite. These (re)crystallized grains are Sr-rich and show mosaic zoning. These results indicate that allanite crystals retained their chemical and isotopic characteristics, and thus their Permian age, as a result of strong strain partitioning between the epidote group porphyroclasts and the eclogite facies matrix in HP-mylonites. The observed partial mobilization of Pb isotopes, which lead to the Cretaceous-aged rims or grains in undeformed samples was facilitated by (re)crystallization of allanite and not by mere Pb diffusion alone under the HP conditions.
ISSN: 0024-4937

insu-00679730
https://hal-insu.archives-ouvertes.fr/insu-00679730
DOI : 10.1016/j.lithos.2011.01.005

Éditeur(s) : HAL CCSD Royal Society of Chemistry
Résumé : International audience
Abstract View references (60)Here we show the long term performance at high temperatures of a multi-tube module containing 8 membranes in 4 parallel lines with a total of 545 cm 2 area. The membranes were prepared via thin film dip coating of cobalt oxide silica (CoO xSi) sol-gel on tubular alumina supports. A preliminary study found that the sol-gel containing 20 mol% cobalt oxide formed the best microporous structure with the highest surface area and pore volume. All resulting membranes delivered permeances of ∼1 × 10 -7 mol m -2 s -1 Pa -1 at 500 °C, indicating a high repeatability for the membrane fabrication process. The permselectivities of helium (He) and hydrogen (H 2) over carbon dioxide (CO 2) and nitrogen (N 2) increased from 10-20 at 100 °C to values close to 1000 at 500 °C. Additionally, the apparent energies of activation (E act) for the smaller kinetic diameter gases He and H 2 at 12.2 and 19.5 kJ mol -1 were high and contrary to the negative values for larger gases N 2 and CO 2 at -1.8 and -7.4 kJ mol -1. These remarkable results were attributed to the molecular sieving mechanism of the microporous silica which was enhanced by the embedding of cobalt oxide into the matrix, delivering structural control with an average pore size of 3 Å. The E act for H 2 permeance was higher than that of He, indicating that the cobalt oxide played an important role in H 2 transport. Two membrane lines performed exceptionally well for binary gas mixture processing with H 2 purity reaching values close to 100% in the permeate stream for argon (Ar) concentrations of up to 80% in the retentate stream. A major finding here is that the binary gas selectivity was independent of temperature, contrary to the permselectivity observed for single gas permeance. Further, the H 2 flow rate was greatly affected by the concentration of Ar in the mixture, while the temperature dependency played only a marginal role. In particular, competitive adsorption in the percolative pathways containing pore constrictions or bottlenecks of the anisotropic CoO xSi matrix allowed Ar to impede H 2 diffusion. Finally, the CoO xSi membranes proved thermally stable and robust for 2000 h of testing for various thermal cycles up to 500 °C. © The Royal Society of Chemistry 2012.
Energy & Environmental Science

hal-01683199
https://hal.univ-antilles.fr/hal-01683199
DOI : 10.1039/c2ee03247c

Éditeur(s) : HAL CCSD
Résumé : Differentiation of the lithospheric mantle occurred principally through partial melting and extraction of melts. Harzburgites are generally considered as melting residues whereas lherzolites are regarded as pristine mantle weakly affected by melting. However, some orogenic peridotites show evidence of igneous refertilization. In this context, this work re-investigates the nature of the Lherz lherzolites (Pyrenees), type-locality of lherzolites, described as a piece of preserved fertile mantle. Structural and geochemical data show that these lherzolites are not pristine but formed through a refertilization reaction between MORB-like melts and refractory lithosphere. Moreover, the Lherz peridotites were partly used to infer the composition of the primitive upper mantle and these results may have important implications for the nature of the late veneer. Additionally, crystal-preferred orientations of minerals (CPO) highlight a strong feedback between melt percolation and finite strain in the percolated rocks. CPO variations are ruled by a subtle balance between instantaneous melt fraction and local strain rate. This work also investigated the effect of melt percolation on Hf, Nd and Sr isotopes. Isotope systematics in Lherz shows that strong isotopic decoupling may arise in a percolation front. The modelling suggests that decoupled isotopic signatures are generated during porous flow and governed by the melt/matrix elements concentrations, chemical diffusivities or efficiency of isotopic homogenization. Melt-rock interactions can generate “intraplate-like” isotopic signatures. This suggests that a part of isotopic signatures of mantle-related rocks could be generated by diffusional processes associated with melt transport.
La différenciation du manteau lithosphérique s'est principalement effectuée par fusion partielle et extraction de magmas. Les harzburgites sont généralement considérées comme des résidus de fusion tandis que les lherzolites représentent un manteau originel peu affecté par la fusion. Cependant, certaines péridotites orogéniques montrent des évidences de refertilisation magmatique. Dans ce contexte, ce travail remet en cause la nature primitive des lherzolites de Lherz (Pyrénées), lithotype de la lherzolite. Les données structurales et géochimiques montrent que ces lherzolites ne représentent pas du manteau préservé mais qu'elles se sont formées par une réaction de refertilisation entre un liquide de type MORB et une lithosphère réfractaire. De plus, les orientations préférentielles de réseau (OPR) des minéraux montrent une forte interdépendance entre la percolation magmatique et la déformation finie observée dans les roches percolées. Les variations d'OPR semblent contrôlées par un équilibre subtil entre fraction de liquide instantanée et taux de déformation local. Ce travail a également abordé les effets de la percolation magmatique sur la variabilité isotopique des roches mantelliques. Les résultats montrent qu'un fort découplage isotopique peut s'observer dans un front de percolation. Les modèles suggèrent que le découplage des signatures isotopiques se produit pendant la percolation, gouverné par le rapport des concentrations en éléments (magma/matrice), la diffusivité et le processus d'homogénéisation isotopique. Les réactions magma/roche peuvent générer des signatures isotopiques de type intraplaque. Ceci suggère qu'une part des signatures isotopiques des roches ou laves dérivées du manteau pourrait être produite par des processus diffusionnels associés à du transport de magma.
https://tel.archives-ouvertes.fr/tel-00431325

tel-00431325
https://tel.archives-ouvertes.fr/tel-00431325
https://tel.archives-ouvertes.fr/tel-00431325/document
https://tel.archives-ouvertes.fr/tel-00431325/file/theseleroux08.pdf

É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