Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Upon subduction, the oceanic crust transforms to blueschists and eclogites, with seismic properties that gradually become similar to those of the surrounding mantle. In order to evaluate the anisotropy of blueschists and glaucophane-bearing eclogites, the elastic constants of glaucophane single-crystal plates from the Sesia-Lanzo Zone (Aosta Valley, Western Alps) were measured using Brillouin spectroscopy at ambient conditions. The mean P- and S-wave velocities are 7.8 and 4.6 km s−1 respectively, and the anisotropy is high (38.1% (AVP) and 27.3% (AVS)). Glaucophane develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the {110} poles concentrated sub-perpendicular to the foliation in both blueschist and eclogite rocks. The measured LPO is in good agreement with viscoplastic selfconsistent numerical models. Seismic properties of glaucophane-bearing blueschist and eclogite are calculated by combining measured LPO and the single-crystal elastic moduli of glaucophane with the other main mineral constituents of the rock: mostly epidote for blueschist, and omphacite and garnet for eclogite. Blueschists present stronger anisotropies (AVP=16.1% and AVS=10.3%) than eclogites (AVP=2.9% and AVS=1.7%). The shear-wave splitting and resulting delay times for a 7-km thick layer of eclogite or blueschist are low for the eclogite (b0.03 s), but signi!cant for blueschist (0.16 s). Application to the subducted oceanic crust yields a decrease of velocity contrast with the surrounding mantle and of anisotropy at depth depending on the temperature of the slab. The details of the velocity variations in subducted oceanic crust are dif!cult to reconcile with the blueschist-eclogite transition as probed by exhumed metamorphic rocks and may require additional phases such as lawsonite or chemical variations such as oxidation state.
ISSN: 0040-1951

hal-00564518
https://hal.archives-ouvertes.fr/hal-00564518
DOI : 10.1016/j.tecto.2010.09.011

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Peridotite xenoliths brought up to the surface by the volcanism of the Kerguelen Islands represent a mantle that has been affected by a high degree of partial melting followed by intense melt percolation above the Kerguelen plume. These xenoliths are therefore particularly suitable to investigate effects of melt-rock interaction on crystallographic fabrics (lattice-preferred orientation (LPO)) of peridotite minerals and on the LPO-induced seismic properties of peridotites above a mantle plume. We have studied a suite of 16 ultramafic samples representative of different degrees of partial melting and magma-rock interaction among which the protogranular harzburgites are the least metasomatised xenoliths and dunites are the ultimate stage of metasomatism. Olivine LPO is characterized by high concentration of [010] axes perpendicular to the foliation and [100] axes close to the lineation or distributed in the foliation plane in harzburgites, whereas the high concentration of [100] axes is parallel to the lineation and [010] axes is perpendicular to the assumed foliation in dunites. Olivine LPO in harzburgites is interpreted as being due to a deformation regime in axial compression or transpression. The fabric strength of olivine decreases progressively from protogranular to poikilitic harzburgites and finally to dunites, for which it remains nevertheless significant (J index 3.8). Seismic properties calculated from LPO of minerals indicate that metasomatism at higher melt/rock ratio lowers the P wave velocities. The most significant difference between harzburgites and dunites corresponds to the distribution of S wave anisotropy. Harzburgites display the maximum of anisotropy within the foliation plane and the minimum of anisotropy perpendicular to the foliation plane, whereas the lowest anisotropy is parallel to the lineation for dunites. These modifications of seismic properties as a result of metasomatic processes may induce seismic heterogeneities in the mantle above the Kerguelen plume. In addition, assuming a lithospheric mantle primarily harzburgitic and structured with a horizontal foliation, the seismic properties calculated for the Kerguelen xenoliths reconcile the rather high anisotropy evidenced by the horizontally propagating surface waves with the apparent isotropy revealed by the absence of splitting of vertically propagating teleseismic SKS waves recorded by the GEOSCOPE Kerguelen station.
ISSN: 1525-2027

hal-00366222
https://hal.archives-ouvertes.fr/hal-00366222
https://hal.archives-ouvertes.fr/hal-00366222/document
https://hal.archives-ouvertes.fr/hal-00366222/file/bascou_Kerguelen_2007GC001879_hal.pdf
DOI : 10.1029/2007GC001879

Éditeur(s) : HAL CCSD
Résumé : International audience
In orogens, shortening is mainly accommodated by thrusts, which constitute preferential zones for fluid–rock interactions. Fluid flow, mass transfer, and mineralogical reactions taking place along thrusts have been intensely investigated, especially in sedimentary basins for petroleum and uranium research. This study combines petrological investigations, mineralogical quantifications, and geochemical characterizations with a wide range of analytical tools with the aim of defining the fluid properties (nature, origin, temperature, and redox) and fluid–host rock interactions (mass transfers, recrystallization mechanisms, and newly formed synkinematic mineralization) in the Pic-de-Port-Vieux thrust fault zone (Pyrenees, Spain). We demonstrate that two geochemically contrasted rocks have been transformed by fluid flow under low-grade metamorphism conditions during thrusting. The hanging-wall Triassic red pelite was locally bleached, while the footwall Cretaceous dolomitic limestone was mylonitized. The results suggest that thrusting was accompanied by a dynamic calcite recrystallization in the dolomitic limestone as well as by leaching of iron via destabilization of iron oxides and phyllosilicate crystallization in the pelite. Geochemical and physical changes highlighted in this study have strong implications on the understanding of the thrust behavior (tectonic and hydraulic), and improve our knowledge of fluid–rock interactions in open fluid systems in the crust.
Contributions to Mineralogy and Petrology

hal-01575548
https://hal-univ-fcomte.archives-ouvertes.fr/hal-01575548
DOI : 10.1007/s00410-017-1394-5

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The hydraulic and electrical properties of pyroclastic rocks have been investigated in laboratory on a representative sampling of Montagne Pelée (Martinique, France) deposits with renewed interest in geophysical applications. This sampling covers all the lithologic units of this volcano: lava dome and lava flows, pumices from ash-and-pumice fall and flow deposits, lava blocks from block-and-ash flow and Peléean “nuées ardentes” deposits, scoriae from scoria flow deposits. The connected porosity varies over a wide range from 3 to 62%. The unconnected porosity is important only on pumices where it can reach 15%. The permeability covers more than 5 orders of magnitude, ranging from 10−16 to 35 × 10−12 m2. The higher values are obtained on lava blocks and the scoriae, even if these rocks are less porous than the pumices. The formation factor ranges from 7 to 1139. The transport properties of these rocks are slightly correlated with porosity. This indicates that these properties are not only controlled by the connected porosity. To connect the transport properties to the textural characteristics of the pore network of pyroclastic rocks, different models, based on geometrical considerations or percolation theory, were tested. The pore access radius distribution and the tortuosity control the transport properties of pyroclastic rocks. Consequently, the models (electric and hydraulic) based on the concept of percolation (e.g., the models of Katz and Thompson), apply better than the equivalent channel model of Kozeny-Carman. In addition, the difference in transport properties observed on lava blocks and pumices confirms that the mechanisms of degassing and vesiculation are different for these two types of rock.
ISSN: 2169-9313

hal-00310237
https://hal.archives-ouvertes.fr/hal-00310237
DOI : 10.1029/2006JB004385, 2007

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Thermal diffusivity of synthetic quartz single-crystals and natural polycrystalline aggregates were measured using a modified Angstrom method at ambient pressure and up to 800 degrees C in order to investigate the effect of crystal preferred orientation on thermal diffusivity of quartz rich rocks. The thermal diffusivity tensor of quartz was determined by measuring the thermal diffusivity of an oriented synthetic quartz single-crystal along the [100] and [001] directions. A naturally deformed, fine-grained quartzite mylonite, was used to compare the single crystal properties with aggregate properties. Using the crystal preferred orientations (CPO) measured by the Electron Back Scattered Diffraction (EBSD) method and the measured thermal diffusivity tensor of quartz, the three-dimensional aggregate thermal diffusivity was modelled using volume averaging models. In order to test the validity of these models that are only based on the properties of the constituent minerals, the thermal diffusivity was then measured on oriented cores as a function of temperature, along three orthogonal directions, parallel to the (X-lineation, Y, Z-pole to foliation) structural frame. Both modelled and direct measurement of thermal diffusivity show an anisotropy of thermal diffusivity, heat transport being more effective parallel to Y-direction and less effective parallel to the lineation X-direction. Thermal diffusivities measured in a given direction on the natural rock are lower by less than 10% than the modelled thermal diffusivity in the same direction, showing that rock imperfections have a weak effect on heat transport. In addition, the modelled thermal diffusivity anisotropy based on the characterization of crystal orientations and thermal diffusivity tensor shows a good agreement with measurements in both alpha and beta stability field. This validation of the petrophysical simulation allows for evaluating the 3D thermal diffusivity of quartz rich rocks in different crustal conditions. The 3D thermal diffusivity of five naturally deformed quartzites displaying quartz fabrics expected in different crustal pressure and temperature ranges were simulated. The resulting thermal diffusivity anisotropies of quartzites range between 40% and 10%, depending on temperature, crystallographic fabric patterns and fabric strength. Under crustal conditions, quartzites deformed in the granulite facies display a high thermal diffusivity in the (XZ) foliation plane, whereas quartzites deformed in greenschist facies have highest thermal diffusivity in the (Z) direction normal to the foliation plane.
ISSN: 0040-1951

hal-00413568
https://hal.archives-ouvertes.fr/hal-00413568
DOI : 10.1016/j.tecto.2008.11.006

É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 Elsevier
Résumé : International audience
The present study aims to apply the AMS method (Anisotropy of Magnetic Susceptibility) at a regional scale to track the fluid circulation direction that has produced an iron metasomatism within pre-existing dolomite host rock. The Urgonian formations hosting the Zn-Pb mineralizations in La Florida (Cantabria, northern Spain) have been taken as target for this purpose. Sampling was carried out, in addition to ferroan dolomite host rock enclosing the Zn-Pb mineralizations, in dolomite host rock and limestone to make the comparison possible between magnetic signals from mineralized rocks, where fluid circulation occurred, and their surrounding formations. AMS study was coupled with petrofabric analysis carried out by texture goniometry, Scanning Electron Microscopy (SEM) observations and also Shape Preferred Orientation (SPO) statistics. SEM observations of ferroan dolomite host rock illustrate both bright and dark grey ribbons corresponding respectively to Fe enriched and pure dolomites. SPO statistics applied on four images from ferroan dolomite host rock give a well defined orientation of ribbons related to the intermediate axis of magnetic susceptibility K2. For AMS data, two magnetic fabrics are observed. The first one is observed in ferroan dolomite host rock and characterised by a prolate ellipsoid of magnetic susceptibility with a vertical magnetic lineation. The magnetic susceptibility carrier is Fe-rich dolomite. These features are probably acquired during metasomatic fluid circulations. In Fe-rich dolomite host rock, ‹c› axes are vertical. As a rule, (0001) planes (i.e. planes perpendicular to ‹c› axes) are isotropic with respect to crystallographic properties. So, the magnetic anisotropy measured in this plane should reflect crystallographic modification due to fluid circulation. This is confirmed by the texture observed using the SEM. Consequently, AMS results show a dominant NE-SW elongation interpreted as the global circulation direction and a NW-SE secondary elongation that we have considered as sinuosities of the fluid trajectory. The second type of magnetic fabric is essentially observed in the limestone and characterised by an oblate form of the ellipsoid of magnetic susceptibility, a horizontal magnetic foliation and mixed magnetic susceptibility carriers. It is interpreted as a sedimentary fabric.
ISSN: 0012-821X

insu-00342499
https://hal-insu.archives-ouvertes.fr/insu-00342499
https://hal-insu.archives-ouvertes.fr/insu-00342499/document
https://hal-insu.archives-ouvertes.fr/insu-00342499/file/Essalhi-EPSL-2008.pdf
DOI : 10.1016/j.epsl.2008.11.011

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Predicting the impact of underground engineering on the environment requires the knowledge of natural media at different scales. In particular, understanding basic phenomena controlling the properties of rocks in the presence of complex fluids necessitates a detailed atomic description of the solid/fluid/solid contacts, the subject of Part I of the present study. First, building the solid interspace between two different crystals in a non-periodic situation is achieved using the ab initio and molecular mechanics code GenMol (TM). A description of the fluid confined within the interspace is then derived from the original genetic iterative multi-species (GIMS) algorithm implemented in the same code. This approach consists of equilibrating chemical potentials, cycle after cycle and species after species, between the confined fluid and the free natural fluid. An elementary iteration for a species k consists of different steps incrementing the number N-k of particles k, with other numbers N-k remaining constant. At each step, an optimum fluid composition is obtained by a genetic process distributing the fluid particles on a grid by stochastic shots, followed in fine by a refining process. The effectiveness of the GIMS approach is demonstrated in the case study of a fluid confined between two (0 0 1) kaolinite faces, with apertures h varying between 4 and 10 A. connected to a 9-species external solution [H2O, Cl-, Na+ CO2(aq), NaCl(aq), Ca2+, Mg2+, HCO3-, H3O+] where concentrations are ranging from 55 to 10(-4) mol/L. The results show a drastic variation in the solute/solvent and cations/ions ratios in the confined fluid when aperture h is lowered to less than 1 nm. These results obtained with a very rapid convergence of the iterative algorithm combined with a very competitive genetic optimizer are validated with high precision on a free solution. This description of contacts between crystals is original and unattainable by standard crystal interface approaches. It opens the way to understanding and/or predicting phenomena between crystals in a complex natural environment, such as adhesion or repulsion investigated in the Part II of this study.
ISSN: 0022-0248

hal-00544569
https://hal.archives-ouvertes.fr/hal-00544569
DOI : 10.1016/j.jcrysgro.2010.08.009

Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
Fluid-rock interactions can control earthquake nucleation and the evolution of earthquake sequences. Experimental studies of fault frictional properties in the presence of fluid can provide unique insights into these interactions. We report the first results from experiments performed on cohesive silicate-bearing rocks (microgabbro) in the presence of pressurized pore fluids (H2O, drained conditions) at realistic seismic deformation conditions. The experimental data are compared with those recently obtained from carbonate-bearing rocks (Carrara marble). Contrary to theoretical arguments, and consistent with the interpretation of some field observations, we show that frictional melting of a microgabbro develops in the presence of water. In microgabbro, the initial weakening mechanism (flash melting of the asperities) is delayed in the presence of water; conversely, in calcite marble the weakening mechanism (brittle failure of the asperities) is favored. This opposite behavior highlights the importance of host-rock composition in controlling dynamic (frictional) weakening in the presence of water: cohesive carbonate-bearing rocks are more prone to slip in the presence of water, whereas the presence of water might delay or inhibit the rupture nucleation and propagation in cohesive silicate-bearing rocks.
ISSN: 0091-7613

hal-01054327
https://hal.archives-ouvertes.fr/hal-01054327
DOI : 10.1130/G34916.1

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Seismic tomography and high surface heat flow models indicate an anomalously thin and hot lithospheric mantle beneath the Canadian Cordillera. To understand the thermal and deformational history of the lithospheric mantle under these conditions, we performed a study of the microstructures, crystal preferred orientations (CPO), and seismic anisotropy of mantle xenoliths collected in the Summit Lake volcano in the southern Canadian Cordillera. The xenoliths vary in composition between spinel lherzolites, cpx-rich lherzolites, and wehrlites. All samples are coarse-grained and show clear evidence for annealing, but deformation microstructures are well-preserved in the spinel lherzolites. Olivine CPO is characterized by a strong point-maximum of [100] and by a girdle distribution of [010[ and [001] axes normal to it; CPO strength is variable, spinel-lherzolites have the strongest CPO. Some samples show a weak orthopyroxene CPO consistent with the olivine one. However most pyroxenes have random CPO, interstitial shapes and undeformed appearance, which suggest a secondary origin by reactive melt percolation under static conditions, leading to crystallization of cpx + sp +/- opx at the expenses of olivine, followed by crystallization of cpx + sp +/- ol at the expenses of opx to produce the wehrlites. The studied xenoliths have typical upper mantle seismic properties patterns with fastest propagation of P-waves and polarization of S-waves parallel to the olivine [100] maxima. Maximum P- and S-waves anisotropy for lherzolites are 8 and 5% respectively. Dispersion of the olivine CPO due to neocrystallization of olivine results in weaker anisotropies (6.5 and 3.5%) in wehrlites. Seismic velocities calculated using a geotherm consistent with the equilibrium temperatures of the Summit Lake xenoliths are 3-4% lower than those expected in a normal, 100 km-thick lithosphere, consistent with seismic tomography models for the region. This high geothermal gradient may be partially explained by considering the anisotropy of thermal conductivity of olivine and subhorizontal foliations in the lithospheric mantle, which result in a decrease by up to 15% of the thermal conductivity in the vertical direction.
ISSN: 0040-1951

hal-00618332
https://hal.archives-ouvertes.fr/hal-00618332
DOI : 10.1016/j.tecto.2011.04.014

É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 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
Résumé : Fault zones are major discontinuities in sedimentary basins. Understanding their role on fluid migrations is an essential issue to (i) characterize the mechanisms and kinematics of deformation and (ii) to determine the parameters which control the distribution of energetic or mineral resources. This work applies to faulting under a temperature range of 200-250°C, corresponding to that of deeply buried reservoirs as well as potentially seismogenic fault zones. The studied faults are normal faults affecting the “Grès d’Annot” Formation, a Priabonian-Rupelian siliclastic turbidite succession of the Alpine foreland basin. They are located in two distinct areas: the Moutière-Restefond area in the eastern part of the basin and the Estrop area in the western part. Vitrinite reflectance outside fault zones indicates maximal temperatures of 240-260°C at Moutière-Restefond and 170-200°C at Estrop, i.e. burial depths around 8 km and 6 km, respectively, assuming a mean geothermal gradient of 30°C/km. Burial was due to underthrusting below the Embrunais-Ubaye Alpine nappes, the difference in burial between the two areas traducing the westward taper of the nappe front.The studied faults affect alternating arkosic sandstone beds and pelite layers with offsets from centimeters to decameters. Two types of core zones are recognized: (1) the foliated core zones of the Moutière-Restefond area, where deformation dominated by pressure solution and synkinematic phyllosilicate neoformation resulted in a foliated fabric, and (2) the core zones of the Estrop area, characterized by dilatant quartz-mineralized veins and breccia. Fault and host rock petrophysical properties were measured on drill plugs, using the water porosity technique and nitrogen permeability technique under 2 Mpa of confining pressure. In the Moutière-Restefond fault zones, plug axes were oriented in three orthogonal directions corresponding to the main deformation axes (X, Y, Z) of the foliated arkose to characterize the effect of structural anisotropy. In the Estrop fault zones, where deformation does not show equivalent anisotropy, the plug axes were oriented according to the main fault plane orientation (i.e., parallel to azimuth, parallel to dip and normal). In each fault zone, the porosity of fault rocks is equivalent to that of host rocks. In the Moutière-Restefond fault zones, the plugs parallel to the X and Z axes of deformation have a permeability similar to that of the host rocks, whereas permeability parallel to the Y axis (i.e. parallel to the foliation and the main fault plane) is one order of magnitude higher. By contrast, the permeability values in the Estrop fault zones are isotropic and equivalent to the host rock values. This study highlights the importance of various types of deformation mechanisms and related fluid-rock interactions in determining different fault rock types within a similar protolith, ultimately controlling distinct petrophysical properties of fault zones.
Les zones de failles constituent des discontinuités majeures dans la croûte supérieure. Ce mémoire s’intéresse aux failles actives à des profondeurs de 5-8 km, i.e. dans des conditions pouvant jouer un rôle important (i) sur le comportement mécanique de de la croûte, en particulier sur son potentiel sismogénique et (ii) sur la compartimentation des réservoirs géologiques profondément enfouis.Les failles étudiées sont des failles normales (rejet jusqu’à 50 mètres) affectant des alternances de bancs de grès arkosiques et de pélites dans les Grès d’Annot, une formation turbiditique d’âge priabonien-rupélien du bassin d’avant-chaine alpin. Les Grès d’Annot ont été enfouis sous les nappes de l’Embrunais-Ubaye peu de temps après leur dépôt et exhumés au Miocène moyen-supérieur. L’étude a été réalisée dans deux secteurs des Alpes de Haute-Provence, (1) le secteur de l’Estrop (chaine des Trois Evéchés, dans la partie occidentale du bassin) et (2) le secteur de Point Vert-Restefond (alentours du col de la Bonette, dans la partie orientale du bassin). La réflectance de la vitrinite dans les Grès d’Annot indique des températures maximales de 160-190°C dans le secteur (1) et de 240-260 °C dans le secteur (2) , soit des profondeurs d’enfouissement de 5-6 km et 8-10 km, respectivement, dans l’hypothèse d’un gradient géothermique moyen de 25-30°C/km. A partir des observations microstructurales, de données d’inclusions de fluides, de modélisations thermodynamiques sur chlorites syncinématiques et des indices de cristallinité de l’illite, on a établi que les failles du secteur (1) ont fonctionné à une température d’environ 200°C pour une pression de 0.6Kbar et celles du secteur (2) dans une gamme de température de 200-220°C et de pression comprise entre 0.6 et 1.1 Kba. L’activité de la faille de Point-Vert a été datée par la méthode 40Ar/39Ar au Miocène inférieur (21 Ma). Ces failles sont interprétées comme associées aux premiers stades du soulèvement résultant de l’activité des chevauchements dans le socle sous-jacent.Deux types principaux de roches de failles caractérisent les failles étudiées. Les brèches dilatantes sont majoritaires dans les zones de failles du secteur (1) alors que les arkoses foliées sont majoritaires dans les zones de failles du secteur (2). Les brèches dilatantes témoignent de mécanismes de déformation dominés par la fracturation transgranulaire et la pression-solution des silicates. Ce dernier mécanisme se manifeste principalement par d’abondantes précipitations cycliques de quartz, les stylolites restant peu développés. Les géométries structurales sont associées à l’activation de plans T de Riedel déterminant l’ouverture de relais extensifs et un régime localement dilatant. Nous associons ce type de roche de faille à une faible pression effective au cours de la déformation.La déformation cisaillante exprimée dans les arkoses foliées témoigne de mécanismes de déformation dominés par (i) la fracturation intragranulaire et transgranulaire, (ii) la pressionsolution/ précipitation des silicates, principalement le quartz, (iii) la néoformation de micas blancs liée à l’altération des feldspaths et (iv) la précipitation de chlorite. Ce travail montre également que les transformations minéralogiques liées aux interactions fluide-roche dans la zone de faille (forte augmentation de la proportion des phyllosilicates au dépend des feldspaths) changent radicalementla composition minéralogique de la roche. Nous associons ce type de roche de faille à une forte pression effective au cours de la déformation.La pétrophysique montre que les brèches de failles (secteur 1) ont une perméabilité isotrope et équivalente à celle du protolithe (10-3mD sous 2 MPa de confinement) tandis que les roches foliées (secteur 2) montrent une augmentation de perméabilité (10-2mD) selon l’axe Y de la déformation, d’un ordre de grandeur supérieure aux perméabilités suivant les axes X et Z et dans la roche hôte (10-3mD). Ce drain préférentiel potentiel est associé à l’orientation des phyllosilicates dans la foliation.
https://hal.archives-ouvertes.fr/tel-01310382

tel-01310382
https://hal.archives-ouvertes.fr/tel-01310382
https://hal.archives-ouvertes.fr/tel-01310382/document
https://hal.archives-ouvertes.fr/tel-01310382/file/Th%C3%A8seCAVAILHES2012.pdf

É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
Homogenization models are widely used to predict the evolution of texture (crystal preferred orientations) and resulting anisotropy of physical properties in metals, rocks, and ice. They fail, however, in predicting two main features of texture evolution in simple shear (the dominant deformation regime on Earth) for highly anisotropic crystals, like olivine: (1) the fast rotation of the CPO towards a stable position characterized by parallelism of the dominant slip system and the macroscopic shear and (2) the asymptotical evolution towards a constant intensity. To better predict CPO-induced anisotropy in the mantle, but limiting computational costs and use of poorly-constrained physical parameters, we modified a viscoplastic self-consistent code to simulate the effects of subgrain rotation recrystallization. To each crystal is associated a finite number of fragments (possible subgrains). Formation of a subgrain corresponds to introduction of a disorientation (relative to the parent) and resetting of the fragment strain and internal energy. The probability of formation of a subgrain is controlled by comparison between the local internal energy and the average value in the polycrystal. A two-level mechanical interaction scheme is applied for simulating the intracrystalline strain heterogeneity allowed by the formation of low-angle grain boundaries. Within a crystal, interactions between subgrains follow a constant stress scheme. The interactions between grains are simulated by a tangent viscoplastic self-consistent approach. This two-level approach better reproduces the evolution of olivine CPO in simple shear in experiments and nature. It also predicts a marked weakening at low shear strains, consistently with experimental data.
ISSN: 0012-821X

hal-01243503
https://hal.archives-ouvertes.fr/hal-01243503
DOI : 10.1016/j.epsl.2015.08.018

Éditeur(s) : HAL CCSD
American Geophysical Union, Fall Meeting 2009
San Francisco, United States

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

Éditeur(s) : HAL CCSD Copernicus Publications
Résumé : International audience
During Ocean Drilling Program Leg 134, a mixed sequence of volcaniclastic and carbonate sediments was cored together with basaltic breccias and a single basalt sill in Hole 833B. In this study, data from Schlumberger's Geochemical Logging Tool, together with physical property logs, characterise part of the section and establish a complete log-derived lithostratigraphy. The major transition between the sedimentary rocks and a basaltic sill at 830 mbsf (metres below sea floor) is clearly seen as well as the distinction between the silicate-rich sediments and the carbonate-rich part of the section. The transition between the volcanic sandstone interbedded with basaltic and sed-lithic breccia of Unit III, and the volcanic sandstone of Unit IV is clearly marked by the physical property logs. The accuracy of the logging data was evaluated and discussed with an improved log versus core correlation. The cluster analysis of the complete data set shows a correlation with the data provided by the mechanical caliper device, which in turn indicates the quality of the borehole.
ISSN: 1816-8957

hal-00642165
http://hal.upmc.fr/hal-00642165
http://hal.upmc.fr/hal-00642165/document
http://hal.upmc.fr/hal-00642165/file/Rabaute_et_al-SciDrilling-1997.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Both oedometric and plane-strain tests were performed with parallelepipedic samples made of synthetic granular, cohesive, frictional and dilatant rock analogue material GRAM2. For the first time parallel sets of fractures that have all the characteristics of natural joints were reproduced in the laboratory. The fractures are regularly spaced, normal to σ3, and have plumose morphology very similar to that of natural joints. These fractures can form at tensile stress σ3 much smaller in magnitude than the tensile strength of material and even at slightly compressive σ3. When mean stress σ exceeds a certain value, the fractures become oblique to σ1 (the obliquity increases with σ), forming networks of conjugate shear bands/fractures. These results of plane-strain experiments are in good agreement with those of better controlled conventional axisymmetric tests on a similar material in Chemenda et al. (2011b) and are closer to real geological situations. Both types of experiments are complementary. Their results lead to the conclusion that at least certain categories of natural fractures (including joints, and conjugate shear fractures/bands) were initiated as deformation localization bands. The band orientation is defined by the constitutive properties/parameters (notably the dilatancy factor) that are sensitive to σ.
ISSN: 0040-1951

hal-00795555
https://hal.archives-ouvertes.fr/hal-00795555
DOI : 10.1016/j.tecto.2011.11.021

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Rheological properties of mantle minerals and rocks at temperatures (T) appropriate to much of Earth's lithosphere have remained poorly constrained, even though past experimental studies on olivine single crystals and polycrystalline aggregates have quantified the high-temperature creep mechanisms (T > 1200 degrees C). Consequently, we have performed deformation experiments on crystals of San Carlos olivine at lower temperatures, from 900 degrees to 1200 degrees C, in triaxial compression along the [101](c) direction. The experiments were carried out at a confining pressure of 300 MPa in a high-resolution gas-medium mechanical testing apparatus at differential stresses of 100 to 500 MPa. Several samples were deformed at constant displacement rate and others at constant load, in order to provide insight into possible effects of work-hardening. Under the deformation conditions investigated, little evidence of work-hardening was observed. The data follow a power-law dependence on stress, as in previous high-temperature deformation studies. The samples were, however, considerably weaker than predicted by the experimentally determined high-temperature constitutive equation for olivine crystals of this orientation from the study of Bai et al. (1991). The mechanical behavior correlates instead with the weaker of the two mechanisms (flow laws) that contribute to the high-temperature constitutive equation. Thus, our experiments demonstrate that published high-temperature constitutive equations overestimate the strength of lithospheric mantle and that the transition to low-temperature creep occurs at lower temperatures than previously inferred. Citation: Demouchy, S., S. E. Schneider, S. J. Mackwell, M. E. Zimmerman, and D. L. Kohlstedt (2009), Experimental deformation of olivine single crystals at lithospheric temperatures, Geophys. Res. Lett., 36, L04304, doi: 10.1029/2008GL036611.
ISSN: 0094-8276

hal-00420070
https://hal.archives-ouvertes.fr/hal-00420070
DOI : 10.1029/2008GL036611

Éditeur(s) : HAL CCSD American Chemical Society
Résumé : International audience
Carbonation of ultramafic rocks in geological reservoirs is, in theory, the most efficient way to trap CO2 irreversibly; however, possible feedback effects between carbonation reactions and changes in the reservoir permeability must be considered to realistically assess the efficiency and sustainability of this process. We investigated changes in the hydrodynamic properties of sintered dunite samples by means of percolation experiments, under conditions analogous to that of in situ carbonation. Our results show that carbonation efficiency is controlled by the local renewal of the reactants and the heterogeneity of the pore structure. Preferential flow zones are characterized by the formation of magnetite and of a silica-rich layer at the olivine surfaces, which eventually inhibits olivine dissolution. Conversely, sustainable olivine dissolution together with coprecipitation of magnesite, siderite, and minor Mg-TOT-phyllosilicates, occur in reduced-flow zones. Thus carbonate precipitation only decreases porosity in zones where diffusion-controlled transport is dominant. Consequently, while high flow rates will decrease the carbonation efficiency of the reservoir and low flow rates may reduce the permeability irreversibly close to the injection point, moderate injection rates will ensure a partial carbonation of the rock and maintain the reservoir permeability.
ISSN: 0013-936X

hal-00424540
https://hal.archives-ouvertes.fr/hal-00424540
DOI : 10.1021/es8018429