Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Direct numerical simulations of flow through two millimeter-scale rock samples of limestone and sandstone are performed using three diverse fluid dynamic simulators. The resulting steady-state velocity fields are compared in terms of the associated empirical probability density functions (PDFs) and key statistics of the velocity fields. The pore space geometry of each sample is imaged at 5.06−μm voxel size resolution using X-ray microtomography. The samples offer contrasting characteristics in terms of total connected porosity (about 0.31 for the limestone and 0.07 for the sandstone) and are typical of several applications in hydrogeology and petroleum engineering. The three-dimensional fluid velocity fields within the explicit pore spaces are simulated using ANSYS® FLUENT® ANSYS Inc. (2009), EULAG Prusa et al. (Comput. Fluids 37, 1193–1207 2008), and SSTOKES Sarkar et al. (2002). These computational approaches are highly disperse in terms of algorithmic complexity, differ in terms of their governing equations, the adopted numerical methodologies, the enforcement of internal no-slip boundary conditions at the fluid-solid interface, and the computational mesh structure. As metrics of comparison to probe in a statistical sense the internal similarities/differences across sample populations of velocities obtained through the computational systems, we consider (i) integral quantities, such as the Darcy flux and (ii) main statistical moments of local velocity distributions including local correlations between velocity fields. Comparison of simulation results indicates that mutually consistent estimates of the state of flow are obtained in the analyzed samples of natural pore spaces despite the considerable differences associated with the three computational approaches. We note that in the higher porosity limestone sample, the structures of the velocity fields obtained using ANSYS FLUENT and EULAG are more alike than either compared against the results obtained using SSTOKES. In the low-porosity sample, the structures of the velocity fields obtained by EULAG and SSTOKES are more similar than either is to the fields obtained using ANSYS FLUENT. With respect to macroscopic quantities, ANSYS FLUENT and SSTOKES provide similar results in terms of the average vertical velocity for both of the complex microscale geometries considered, while EULAG tends to render the largest velocity values. The influence of the pore space structure on fluid velocity field characteristics is also discussed.
ISSN: 1420-0597

hal-01174152
https://hal.archives-ouvertes.fr/hal-01174152
DOI : 10.1007/s10596-015-9486-7

Éditeur(s) : HAL CCSD Oxford University Press (OUP): Policy B - Oxford Open Option A Oxford University Press (OUP)
Résumé : International audience
Peridotite xenoliths exhumed by Quaternary alkaline magmatism in the Tahalgha district, southern Hoggar, represent fragments of the subcontinental lithospheric mantle beneath the boundary between the two major structural domains of the Tuareg Shield: the ‘Polycyclic Central Hoggar’ to the east and the ‘Western Hoggar’, or ‘Pharusian Belt’, to the west. Samples were collected from volcanic centres located on both sides of a major lithospheric shear zone at 4°35' separating these two domains. Although showing substantial variations in their deformation microstructures, equilibrium temperatures and modal and chemical compositions, the studied samples do not display any systematic changes of these features across the 4°35' fault. The observed variations rather record small-scale heterogeneities distributed throughout the study area and reflecting the widespread occurrence of vein conduits and metasomatized wall-rocks related to trans-lithospheric melt circulation during the Cenozoic. These features include partial annealing of pre-existing deformation microstructures, post-deformation metasomatic reactions, and trace-element enrichment, coupled with heating from 750–900°C (low-temperature lherzolites) to 900–1150°C (intermediate-T lherzolites and high-T harzburgites and wehrlites). Trace-element modelling confirms that the range of rare earth element (REE) variations observed in the Tahalgha clinopyroxenes may be accounted for by reactive porous flow involving a single stage of basaltic melt infiltration into a light REE (LREE)-depleted protolith. Whole-rock compositions record the final entrapment of disequilibrium metasomatic melts upon thermal relaxation of the veins–wall-rock system. The striking correlations between equilibrium temperatures and trace-element enrichment favor a scenario in which the high-temperature peridotites record advective heat transport along melt conduits, whereas the intermediate- and low-temperature lherzolites reflect conductive heating of the host Mechanical Boundary Layer. This indicates that the lithosphere did not reach thermal equilibrium, suggesting that the inferred heating event was transient and was rapidly erased by thermal relaxation down to the relatively low-temperature present-day geotherm. The low-T (<900°C) deformed lherzolites (porphyroclastic to equigranular) are characterized by only incipient annealing and LREE-depleted clinopyroxene compositions. They were only weakly affected by the Cenozoic events and could represent relatively well-preserved samples from rejuvenated Pan-African lithosphere. Extensive lithospheric rejuvenation occurred either regionally during the Pan-African orogeny, as a result of lithospheric delamination or thermomechanical erosion after thickening, or more locally along the meridional shear zones. The low-T Tahalgha lherzolites are comparable with lherzolites from Etang de Lherz, southern France, interpreted as lithospheric mantle rejuvenated by melt-induced refertilization during a late stage of the Variscan orogeny.
ISSN: 0022-3530

hal-01107410
https://hal.archives-ouvertes.fr/hal-01107410
DOI : 10.1093/petrology/egu056

Éditeur(s) : HAL CCSD American Association of Petroleum Geologists
Résumé : In reservoir engineering, hydrodynamic properties can be estimated from downhole electrical data using heuristic models (e.g., Archie and Kozeny-Carman's equations) relating electrical conductivity to porosity and permeability. Although proven to be predictive for many sandstone reservoirs, the models mostly fail when applied to carbonate reservoirs that generally display extremely complex pore network structures. In this article, we investigate the control of the three-dimensional (3-D) geometry and morphology of the pore network on the electrical and flow properties, comparing core-scale laboratory measurements and 3-D x-ray microtomography image analysis of samples from a Miocene reefal carbonate platform located in Mallorca (Spain). The results show that micrometer- to centimeter-scale heterogeneities strongly influence the measured macroscopic physical parameters that are then used to evaluate the hydrodynamic properties of the rock, and therefore, existing models might not provide accurate descriptions because these heterogeneities occur at scales smaller than those of the integration volume of the borehole geophysical methods. However, associated with specific data processing, 3-D imagery techniques are a useful and probably unique mean to characterize the rock heterogeneity and, thus, the properties variability.
ISSN: 0149-1423

hal-00950066
https://hal.archives-ouvertes.fr/hal-00950066
DOI : 10.1306/05221312134

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
A reactive percolation experiment was conducted by injecting seawater into a permeable olivine aggregate at 190 °C and 19 MPa to explore the relationships between hydration reactions and hydrodynamic properties during the onset of serpentinization in the ultramafic lithosphere exposed at mid-ocean ridges. The experiment was stopped after 23 days when the sample became impermeable. The initial flow rate was 0.2 mL/h and then was decreased down to 0.06 mL/h after 8 days. Permeability decreased continuously throughout the experiment. The analyses of fluid chemistry time series and of the mineralogy and structure of the reacted sample showed olivine dissolution and precipitation of proto-serpentine, brucite and Fe-oxides. These reactions are controlled by coupled hydrodynamic and chemical processes interacting at different time and spatial scales. Their first characteristic is the production of silica rich outlet fluids in disequilibrium with the observed reacted mineral assemblages. These compositions are interpreted as resulting from a suite of coupled dissolution-precipitation reactions, controlled at the pore scale by surface kinetic processes, that rapidly reaches steady-state (constant fluid composition independent of fluid flux). Hence, the effective rate of serpentinization was controlled, to the first order, by the transport of reactants during the experiment. Mass balance calculations show that the rate of olivine conversion was fast (0.2-0.5 wt%/day), yet only ∼8 wt% of the olivine sample reacted, because the permeability drop limited fluid circulation. Porosity varied little during the experiment compared to permeability changes: the decrease of permeability was explained by the structure of the newly formed serpentine favouring the clogging of fluid paths. The rate at which permeability decreased was the fastest at low flux conditions. This suggests that permeability changes were not controlled simply by the kinetics of the serpentinization reactions, but mainly reflected the development of low Peclet pore scale zones forming micro-environments that control the serpentinization reaction-paths. These results outline the need to take into account the coupling of hydrodynamic and chemical processes when modelling the effects of reactive transport within the hydrating oceanic lithosphere.
ISSN: 0012-821X

hal-00855043
https://hal.archives-ouvertes.fr/hal-00855043
DOI : 10.1016/j.epsl.2013.03.052

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Large sand intrusions often exhibit conical morphologies analogous to magmatic intrusions such as saucer-shaped or cup-shaped sills. Whereas some physical processes may be similar, we show with scaled experiments that the formation of conical sand intrusions may be favoured by the pore-pressure gradients prevailing in the host rock before sand injection. Our experiments involve injecting air into a permeable and cohesive analogue material to produce hydraulic fractures while controlling the pore pressure field. We control the state of overpressure in the overburden by applying homogeneous basal pore pressure, and then adding a second local pore pressure field by injecting air via a central injector to initiate hydraulic fractures near the injection point. In experiments involving small vertical effective stresses (small overburden, or high pore fluid overpressure), the fracturing pressure (lambda(fract)) is supralithostatic and two dipping fractures are initiated at the injection point forming a conical structure. From theoretical considerations, we predict that high values of lambda(fract), are due to strong cohesion or high pore fluid overpressure distributed in the overburden. Such conditions are favoured by the pore pressure/stress coupling induced by both pore pressure fields. The dips of cones can be accounted for elastic-stress rotation occurring around the source. Contrary to magmatic chamber models, the aqueous fluid overpressure developed in a parent sandbody (and prevailing before the formation of injectites) may diffuse into the surrounding overburden, thus favouring stress rotation and the formation of inclined sheets far from the parent source. For experiments involving higher vertical effective stresses (thick overburden or low pore fluid overpressure), the fracturing pressure is lower than the lithostatic stress, and a single fracture is opened in mode I which then grows vertically. At a critical depth, the fracture separates into two dilatant branches forming a flat cone. We make use of a P.I.V. (Particle Imaging Velocimetiy) technique to analyse plastic deformation, showing that these inclined fractures are opened in mixed modes. Close to the surface, they change into steep shear bands where fluids can infiltrate. The final morphology of the fracture network is very similar to the common tripartite architecture of various injection complexes, indicating that different mechanisms may be involved in the formation of dykes. Feeder dykes under the sill zones may open as tensile fractures, while overlying dykes may be guided by the deformation induced by the growth of sills. These deformation conditions may also favour the formation of fluid escape structures and pockmarks. (C) 2011 Elsevier B.V. All rights reserved.
ISSN: 0012-821X

hal-00745640
https://hal.archives-ouvertes.fr/hal-00745640
DOI : 10.1016/j.epsl.2011.10.029

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Mid-ocean ridge basalt (MORB) is the most abundant magma on Earth, and provides a geochemical window into the mantle. Deriving mantle composition and melting processes from the erupted lavas requires correction to be made for their evolution as they pass through and generate the oceanic crust. This is typically done by assuming that modification of melts in crustal magma chambers occurs exclusively by fractional crystallisation. However, extensive mineral major- and trace element data from a full section of fast-spread lower crustal rocks exposed in Hess Deep (equatorial Pacific Ocean) demonstrate that their evolution is instead controlled by reactive porous flow. These reactions lead to a strong enrichment in, and fractionation of, incompatible trace elements in the melt (as recorded by clinopyroxene compositions), leading to melt compositions far outside of the compositional realm of MORB both in terms of trace element abundances and ratios. The reactive signature increases in strength up section, peaking in varitextured gabbros interpreted to represent the fossilised axial melt lens, indicating that reactive porous flow occurred on the scale of the entire lower crust. The enrichment of the melt is coupled with a strong trace element depletion of plagioclase, olivine, and, to a lesser extent, clinopyroxene cores, suggesting that these phases represent the residues of the reactions from which trace elements have been removed. The dominant role of reactive porous flow, and the resulting deviations from fractional crystallisation predictions, suggest that the lower oceanic crust plays a much more complex and significant role in modifying the compositions of MORB than previously expected, with consequent implications for models of mantle processes.
ISSN: 0012-821X

hal-00814075
https://hal.archives-ouvertes.fr/hal-00814075
DOI : 10.1016/j.epsl.2012.11.012

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We present a set of single-well injection withdrawal tracer tests in a paleoreef porous reservoir displaying important small-scale heterogeneity. An improved dual-packer probe was designed to perform dirac-like tracer injection and accurate downhole automatic measurements of the tracer concentration during the recovery phase. By flushing the tracer, at constant flow rate, for increasing time duration, we can probe distinctly different reservoir volumes and test the multiscale predictability of the (non-Fickian) dispersion models. First we describe the characteristics, from microscale to meter scale, of the reservoir rock. Second, the specificity of the tracer test setup and the results obtained using two different tracers and measurement methods (salinity-conductivity and fluorescent dye­optical measurement, respectively) are presented. All the tracer tests display strongly tailed breakthrough curves (BTC) consistent with diffusion in immobile regions. Conductivity results, measured over 3 orders of magnitude only, could have been easily interpreted by the conventional mobile-immobile (MIM) diffusive mass transfer model of asymptotic log-log slope of 2. However, the fluorescent dye sensor, which allows exploring much lower concentration values, shows that a change in the log-log slope occurs at larger time with an asymptotic value of 1.5, corresponding to the double-porosity model. These results suggest that the conventional, one-slope MIM transfer rate model is too simplistic to account for the real multiscale heterogeneity of the diffusion-dominant fraction of the reservoir.
ISSN: 0043-1397

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

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Energy dissipation is observed on seismic data when a wave propagates through a porous medium, involving different frequency regimes depending on the nature of rock and fluid types. We focus here on the role of partial fluid saturation in unconsolidated porous media, looking in particular at P-wave phase velocity and attenuation. The study consists in running an experiment in a sand-filled tank partially saturated with water. Seismic propagation in the tank is generated in the kHz range by hitting a steel ball on a granite plate. Seismic data are recorded by buried accelerometers and injecting or extracting water controls the partial saturation. Several imbibition/drainage cycles were performed between the water and gas residual saturations. A Continuous Wavelet Transform applied on seismic records allowed us to extract the direct P-wave at each receiver.We observe an hysteresis in phase velocities and inverse quality factors between imbibition and drainage. Phase velocities and inverse quality factors are then jointly inverted to get a final poro-visco-elastic model of the partially saturated sand that satisfactorily reproduces the data. The model formulation consists in generalizing the Biot theory to effective properties of the fluid and medium (permeability and bulk modulus) in order to properly explain the phase velocity variation as a function of the saturation. The strong level of attenuation measured experimentally is further explained by an anelastic effect due to grain to grain sliding, adding to Biot's losses. This study shows that fluid distribution at microscopic scale has strong influence on the attenuation of direct P-waves at macroscopic scale and confirms that seismic prospection may be a powerful tool for the characterization of transport phenomena in porous media.
ISSN: 0956-540X

hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920/document
https://hal.archives-ouvertes.fr/hal-00779920/file/Barriere_et_al_HAL.pdf
DOI : 10.1111/j.1365-246X.2012.05691.x

Éditeur(s) : Pergamon-elsevier Science Ltd
Résumé : Society's needs for a network of in situ ocean observing systems cross many areas of earth and marine science. Here we review the science themes that benefit from data supplied from ocean observatories. Understanding from existing studies is fragmented to the extent that it lacks the coherent long-term monitoring needed to address questions at the scales essential to understand climate change and improve geo-hazard early warning. Data sets from the deep sea are particularly rare with long-term data available from only a few locations worldwide. These science areas have impacts on societal health and well-being and our awareness of ocean function in a shifting climate. Substantial efforts are underway to realise a network of open-ocean observatories around European Seas that will operate over multiple decades. Some systems are already collecting high-resolution data from surface, water column, seafloor, and sub-seafloor sensors linked to shore by satellite or cable connection in real or near-real time, along with samples and other data collected in a delayed mode. We expect that such observatories will contribute to answering major ocean science questions including: How can monitoring of factors such as seismic activity, pore fluid chemistry and pressure, and gas hydrate stability improve seismic, slope failure, and tsunami warning? What aspects of physical oceanography, biogeochemical cycling, and ecosystems will be most sensitive to climatic and anthropogenic change? What are natural versus anthropogenic changes? Most fundamentally, how are marine processes that occur at differing scales related? The development of ocean observatories provides a substantial opportunity for ocean science to evolve in Europe. Here we also describe some basic attributes of network design. Observatory networks provide the means to coordinate and integrate the collection of standardised data capable of bridging measurement scales across a dispersed area in European Seas adding needed certainty to estimates of future oceanic conditions. Observatory data can be analysed along with other data such as those from satellites, drifting floats, autonomous underwater vehicles, model analysis, and the known distribution and abundances of marine fauna in order to address some of the questions posed above. Standardised methods for information management are also becoming established to ensure better accessibility and traceability of these data sets and ultimately to increase their use for societal benefit. The connection of ocean observatory effort into larger frameworks including the Global Earth Observation System of Systems (GEOSS) and the Global Monitoring of Environment and Security (GMES) is integral to its success. It is in a greater integrated framework that the full potential of the component systems will be realised. Crown Copyright (C) 2011 Published by Elsevier Ltd. All rights reserved.
Progress In Oceanography (0079-6611) (Pergamon-elsevier Science Ltd), 2011-10 , Vol. 91 , N. 1 , P. 1-33

Droits : 2011 Published by Elsevier Ltd. All rights reserved.
http://archimer.ifremer.fr/doc/00044/15557/15923.pdf
DOI:10.1016/j.pocean.2011.05.001
http://archimer.ifremer.fr/doc/00044/15557/

Éditeur(s) : HAL CCSD
Résumé : This study is designed to improve the understanding of structural, mechanical and petrophysical characteristics of deformation bands in porous sandstone such as their potential influence on fluid storage or migration in reservoir setting. Structural and petrophysical analysis of band networks from Provence show that permeability of band is mainly controlled by the degree of cataclasis and diagenetic processes at depth. Map-scale analysis of band networks underlines the major influence of tectonic regime on band distribution, organization and permeability. Band associated to normal-fault regime (cataclasis and low-permeability) are located around map-sclae-fault whereas band associated to thurst-fault regime (moderate cataclasis and permeability) are pervasively distributed in all the area. Reverse-sense band are cataclastic and low-permeabiilty around the Roquemaure thrust, which show the potential role of large-scale fault on structural and petrophysical characteristics of bands. This analysis also shows the influence of fault propagation on formation of low permeability catacastic bands. This analysis suggests also granulometry as factor controlling organization of band networks.A mechanical analysis of deformation bands is realized based on results from Provence. Yield envelopes of these materials are calculated from analytical solution. Triaxial tests done on sandstone (L'Etang quarry, Orange) confirm the shape of yield envelope theoretically calculated. Triaxial tests done on the poorly-lithified sands (Boncavaï quarry, Uchaux) show grain size, sorting and packing influencing yield envelope and cataclastic band localization at shallow burial depth. Stress paths are calculated for burial and tectonic event (extension and contraction) using the relationship between main stress and secondary stresses (K0) calibrated from well data. These different results are integrated into a model which allows to estimate type and structural characteristics of bands susceptible to form for a known geological setting.A structural, mechanical and petrophysical model is established form these results. This model is calibrated by permeability data from deformation bands of Provence and results from the mechanical model of strain localization. This model is also confronted to a synthesis of permeability data from literature, which confirms the influence of the different controlling parameters integrated to the model. This model is applied for two sites of uranium exploration/exploitation in porous sandstone reservoirs and validated by field observations.
Cette étude est destinée à améliorer la compréhension des caractéristiques structurales, pétrophysiques et mécaniques des bandes de déformation dans les grès poreux ainsi que de leur rôle potentiel sur la migration ou le piégeage des fluides en milieu réservoir. L'analyse structurale et pétrophysique des réseaux de bandes de Provence montre que leur perméabilité découle principalement de l'intensité de la cataclase dans leur microstructure et des processus diagénétiques qui peuvent s'y localiser avec la profondeur. L'analyse régionale de ces réseaux de bandes souligne l'influence majeure du régime tectonique sur leur distribution, leur organisation et leur perméabilité. Les bandes associées au régime normal (cataclactiques et peu perméables) sont localisées autour des failles cartographiques, tandis que les bandes associées au régime inverse (modérement cataclastiques et perméables) sont distribuées dans toute la région. Les bandes inverses sont fortement cataclastiques et imperméables uniquement autour du chevauchement de Roquemaure ce qui montre le rôle potentiel de la présence de grande faille sur les caractéristiques structurales et pétrophysiques de ces bandes. D'autre part, cette analyse suggère que la granulométrie de l'encaissant peut influencer l'initiation et l'organisation des réseaux. L'analyse mécanique des bandes de déformation a ensuite été réalisée à partir des résultats obtenus en Provence. L'enveloppe de plasticité de ces matériaux est calculée à l'aide d'une solution analytique. Les essais triaxiaux menés sur le grès de la carrière de l'Etang (orange) confirment la forme de ces enveloppes estimées théoriquement. Les essais réalisés sur les sables de Boncavaï (Uchaux) montrent que la taille des grains, le tri et de la compaction dans ces matériaux peu lithifiés influent sur leurs enveloppes de plasticité et la localisation des bandes cataclastiques à faibles profondeurs. Les trajets de contraintes sont estimés pour les phases d'enfouissement et de chargement tectonique (extension et compression) à partir de la relation entre contrainte principale et contraintes secondaires (K0), calibrée à partir de données de forage. Ces résultats sont intégrés à un modèle permettant d'estimer la nature et les caractéristiques structurales des bandes susceptibles de se former dans un contexte géologique donné. Un modèle structural, mécanique et pétrophysique a été établir à partir de ces différents résultats. Ce modèle est calibré à partir des données de perméabilité mesurées sur les bandes de Provence et des résultats du modèle mécanique de localisation cité précédemment. Ce modèle a ensuite été confronté à une synthèse des données de perméabilité disponibles dans la littérature et confirme l'influence des différents paramètres de contrôle pris en compte dans le modèle. Ce modèle a enfin pu être appliqué sur deux sites d'explorations/exploitations d'uranium en réservoir gréseux poreux et validé par les résultats obtenus, en concordance avec les observations de terrain.
https://tel.archives-ouvertes.fr/tel-01001389

NNT : 2013MON20067
tel-01001389
https://tel.archives-ouvertes.fr/tel-01001389
https://tel.archives-ouvertes.fr/tel-01001389/document
https://tel.archives-ouvertes.fr/tel-01001389/file/36072_BALLAS_2013_archivage-1_cor.pdf

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2-D studies. We analyze seismicity data from a 13 month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3-D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87±0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early-stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies.
ISSN: 1525-2027

hal-01667198
https://hal.archives-ouvertes.fr/hal-01667198
DOI : 10.1002/2017GC007027

É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

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

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

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We aim at testing and validating a mobile-immobile mass transfer model from a set of single-well injection withdrawal tracer tests in a heterogeneous porous aquifer. By varying the duration of the injection phase, different volumes of aquifer are investigated by the tracer. Hence, we focus the transport model validation not only on reproducing a single breakthrough curve (BTC) but also on the model's capacity to predict the amount of mixing as a function of the volume visited by the tracer. All the BTCs are strongly asymmetric, as expected when dispersion is controlled by diffusive mass transfers between the mobile water and the immobile water part of the porosity. However, the BTC cannot be modeled by a conventional mobile-immobile mass transfer model with a simple power law memory function. To account for that, we implement a continuous time random walk model in which the transition time distribution y (t), which is related to the excursion time probability of the tracer in the immobile domain, is a dual-slope power law distribution. The model best fits the BTC data set with a transitional regime controlled by y(t) t2 and an asymptotic regime characteristic of the conventional double-porosity model with y(t) t1.5 . This work emphasizes that high-resolution concentration measurement and multiple-scale tracer tests are required for assessing solute dispersion models in heterogeneous reservoirs and for subsequently obtaining reliable predictions.
ISSN: 0043-1397

insu-00373698
https://hal-insu.archives-ouvertes.fr/insu-00373698
https://hal-insu.archives-ouvertes.fr/insu-00373698/document
https://hal-insu.archives-ouvertes.fr/insu-00373698/file/Borgne_et_al-2008-Water_Resources_Research.pdf
DOI : 10.1029/2007WR006279

Éditeur(s) : HAL CCSD Institute for Scientific Computing and Information
Résumé : International audience
In this paper, we solve a problem describing the mechanical deformations of a porous medium in the presence of a monophasic linear flow or a two phase nonlinear flow with the purpose of modelizing subsidence of hydrocarbon reservoirs. An essential characteristics of this problem is that the mechanical deformation and the flow have different time scales. In petroleum industry, one uses different very efficient simulators for the flow problem and the mechanical deformations, which enables to handle complex models. Therefore it is necessary to be able to combine as efficiently as possible the exploitation of these simulators. We propose two alternative splitting approaches. The first one is the staggered algorithm used by engineers, which amounts to a Gauss-Seidel method in the one phase linear case. The second approach is based upon the preconditioned conjugate gradient method. We use a numerical multi-scale method in both of these algorithms. We compare these two approaches and we show that the preconditioned conjugate gradient algorithm is faster and more robust than the staggered algorithm. Applying the preconditioned conjugate gradient algorithm therefore seems to compensate for the fact that the inf-sup condition for the mixed discretization method is not satisfied when combining the simulators for the mechanical deformations and for the flow computations.
ISSN: 1705-5105

Droits : http://creativecommons.org/licenses/by/
hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039/document
https://hal.archives-ouvertes.fr/hal-01449039/file/DHLE.pdf

É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
Résumé : International audience
The carbon capture and storage is considered as one of the solutions of reducing the effects of global warming. It consists in injecting large quantities of CO 2 in supercritical form directly into deeply located geological formations. One of the methods involves sequestration in the carbonate sedimentary rocks (sandstones and limestones) which are considered as a promising solution. The morphology of the pore network and solid skeleton defines important macroscopic prop-erties of the rock (permeability, stiffness). Modelling at micro scale of random heterogeneous materials (porous media) should be based on morphological and statistical characterisation 1,2 of the investigated microstructure. The work is focused on a reconstruction method by grains deposit 3 , for various range of pa-rameters, to obtain equivalent model of natural rock samples (oolitic limestone and fontainebleau sandstone). Generated samples satisfy morphological and statistical informations 4 which were obtained by 3D image analysis of X-ray tomography of the natural rock samples. For validation of the modelling numerical computations of mechanical and transfer properties are also provided.
14th European Mechanics of Materials Conference - EMMC14
Gothembourg, Sweden

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

Éditeur(s) : Université des Antilles Études caribéennes
Résumé : En développant, à partir de 1984, un Centre touristique Intégralement Planifié dans une localité de la côte Pacifique de l’état d’Oaxaca, l’État mexicain entendait impulser un développement économique dans une région périphérique et pauvre du pays. Basé sur l’expropriation d’une portion du territoire communal, le processus engendre un certain nombre de résistances au sein des communautés locales. Le centre touristique est aujourd’hui un espace fragmenté, reflet et conséquence des inégalités sociales et des conflits sociaux. Depuis une perspective de géographie sociale, cet article cherche à mettre en lumière les conflits et compromis sociaux qui traversent la société locale. L’appropriation de l’espace, les luttes de pouvoir et les stratégies développées par les différents acteurs seront les entrées privilégiées pour interroger un modèle de planification, « par le haut », ignorant trop souvent les réalités sociales à l’échelle locale.
Since 1984, the Mexican government has tried to stimulate the economy in a poor and peripheral area of the country. They aimed to do so by developing an Integrally Planed Resort (CIP, in Spanish) in a community, located in the Oaxaca state on the Pacific coast. Being based on the expropriation of some parts of the communal territory, this process created some remittances among the local communities. Today, this resort is still a fragmented space and reflects social inequalities and social conflicts. From a social geographical perspective, this article aims to show the social conflicts and social compromises among the local society. Space appropriation, power battles and strategies developed by the different actors, will be the main tools to investigate this planning model. In this model, the government too often has ignored social realities on the local scale.
Al crear, a partir de 1984, un Centro turístico Integralmente Planeado en un poblado de la costa del estado de Oaxaca, el gobierno mexicano pretendió impulsar el desarrollo económico de una región periférica y marginada del país. Basado en la expropiación de una parte del territorio municipal, dicho proceso genera resistencias entre las comunidades locales. Hoy en día, el centro turístico es un espacio fragmentado, reflejo y consecuencia de las desigualdades sociales y conflictos. Desde una perspectiva de geografía social, el presente artículo buscará destacar los conflictos y compromisos sociales que atraviesan la sociedad local. La apropiación del espacio, las luchas de poder y las estrategias desarrolladas por los diferentes actores serán los enfoques privilegiados para cuestionar un modelo de planeación “vertical”, que suele ignorar las realidades sociales a nivel local.
Mexique

Droits : info:eu-repo/semantics/openAccess
urn:doi:10.4000/etudescaribeennes.4069
http://journals.openedition.org/etudescaribeennes/4069

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Crystal growth driven by a flowing solution is modeled for a flow with low Reynolds number using a computational dynamic software. Considering equivalent crystallographic faces, the chemical flux is calculated along upstream and downstream faces. Upstream flux is higher compared to downstream and leads to a symmetry breakdown of the crystal shape and develops mirror symmetry parallel to the flow velocity. Moreover the ratio of these two fluxes (upstream / downstream) gives a quantitative relation between the relative crystal growth rate and the flow velocity. Thus, using an inverse method, the flow direction and velocity can be deduced by the study of the variation of the growth band thicknesses of equivalent crystallographic faces. This new method was applied to the formation of metasomatic tourmalinite associated with a leucogranite sill. The approach is complemented by a study of the chemistry of the tourmaline. In the studied case, the application of the new method gives the high fluid velocities in pores (10-3 - 10-4 m/s) during metasomatism. Equivalent Darcy velocities are estimated and discussed accounting for the major role played by the regional deformation. Finally, a two-stage tectono-hydrodynamic model is proposed for the metsomatism. The first stage is genetically linked to the sill injection, and the second is characterized by a wider event with hydrothermal flow passing along the leucogranite sills.
ISSN: 0012-821X

Droits : info:eu-repo/semantics/OpenAccess
insu-00352307
https://hal-insu.archives-ouvertes.fr/insu-00352307
https://hal-insu.archives-ouvertes.fr/insu-00352307/document
https://hal-insu.archives-ouvertes.fr/insu-00352307/file/Sizaret_et_al_final.pdf
DOI : 10.1016/j.epsl.2009.01.013

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Two in situ concrete/claystone interfaces were sampled at the laboratory level in the Andra Meuse/Haute Marne (France) Underground Research Laboratory (URL) in order to study five years of interactions between Callovian-Oxfordian (COx) claystone and two cementitious materials (concrete bottom slab and shotcrete on the walls of the main gallery), with a specific focus on sulfur. Combined mineralogical, chemical and sulfur isotopic investigations were carried out to define the degree of the perturbation of the sulfur system in the claystone and in both the cementitious materials. At both interfaces, results show that the main perturbation on the claystone side is the formation of scarce μm-sized gypsum, the sulfur content of which is essentially derived from pyrite oxidation. The distribution of gypsum is highly correlated with the fissure network of the damaged zone due to excavation of the gallery. Its presence is also often associated with a loss of cohesion of the concrete/claystone interface. Due to the small amounts of gypsum and its μm-size, measurements were performed by ion microprobe. Adaptations were needed on account of the reactivity of gypsum and sulfates in general under the beam. The use of ion microprobe analysis provided evidence of high local isotopic heterogeneity that could be attributed to kinetic fractionation effects. Some analyses suggest a minor contribution of dissolved sulfates in pore water of claystone and possibly of concrete. The perturbation on the concrete side is marked by a significant increase in the bulk sulfur content within three millimeters of the interface with the claystone, showing a sulfur gradient from claystone to concrete. The main objective of this work was to define the extent of the chemical and mineralogical perturbations, taking into account in situ URL conditions, i.e. hydrodynamic conditions (shotcrete sprayed on the gallery walls and subjected to ventilation of the galleries), damaged zone of claystone induced by excavation and multi-scale heterogeneity. This identifies the relevant parameters to better constrain predictive reactive transport models of concrete/claystone interactions at a realistic scale.
ISSN: 0079-1946

hal-01055790
https://hal-brgm.archives-ouvertes.fr/hal-01055790
DOI : 10.1016/j.pce.2014.01.003