Éditeur(s) : Yale University
Résumé : Sea-surface temperature images of the coastal upwelling regions off Northwest Africa show that the core of upwelling is sometimes located far from the coast. This has been documented in three regions that share a common feature, namely a wide and shallow continental shelf. This upwelling feature plays a key role in the ecology of the Canary Current System. It creates an innerfront which provides retention for biological material, e.g. fish eggs and larvae, in the highly productive nearshore environment. An analytical model has been developed based on a two dimensional extension of Ekman's solution. The linear and steady response of a homogeneous ocean forced by an upwelling-favorable wind provides a mechanism for the upwelling separation from the coast. The merging of the surface and bottom Ekman layers induces a very weak cross-shore circulation and a "kinematic barrier" for the Ekman transport divergence. In the case of an alongshore wind, the barrier is located near the isobath h ≈ 0.4D, where D is the thickness of Ekman layers. This yields an upwelling cell which is essentially concentrated in the region 0.5D < h < 1.25D, with upwelling occurring preferentially near the isobath h ≈ 0.6D. It turns out that the cross-shore width of upwelling scales with D/S, the ratio of Ekman depth to bottom topographic slope. The application of this solution to real bathymetric profiles rationalizes, not only the offshore upwelling observations in Northwest Africa, but also the influence of topography on the cross-shelf structure of a wind-driven coastal upwelling. The model also quantifies the effect of the cross-shore wind component showing how it drives the nearshore pressure gradient adjustment and how it affects the upwelling. A linear numerical experiment reproduces the theoretical steady solution, thereby allowing investigation of the transient regime. Relaxation of the hypothesis in the numerical model validates the linear assumption of the theory and then allows investigation of the sensitivity to friction parameterizations and the influence of stratification. The latter leads to an "oscillation" of the upwelling cell with seaward migration driven by outcropping and homogeneization of the water column, and, coastal incursion driven by a "boundary layers splitting" process caused by shoreward advection of the isopycnal dome and stratification of the inner shelf.
Journal of Marine Research (Yale University), 2008-09 , Vol. 66 , N. 5 , P. 589-616

Droits : 2008 Yale University
http://archimer.ifremer.fr/doc/2008/publication-6334.pdf
DOI:10.1357/002224008787536790
http://archimer.ifremer.fr/doc/00000/6334/

Éditeur(s) : HAL CCSD American Association of Petroleum Geologists
Résumé : International audience
We describe the structure, microstructure, and petrophysical properties of fault rocks from two normal fault zones formed in low-porosity turbiditic arkosic sandstones, in deep diagenesis conditions similar to those of deeply buried reservoirs. These fault rocks are characterized by a foliated fabric and quartz-calcite sealed veins, which formation resulted from the combination of the (1) pressure solution of quartz, (2) intense fracturing sealed by quartz and calcite cements, and (3) neoformation of synkinematic white micas derived from the alteration of feldspars and chlorite. Fluid inclusion microthermometry in quartz and calcite cements demonstrates fault activity at temperatures of 195degC to 268degC. Permeability measurements on plugs oriented parallel with the principal axes of the finite strain ellipsoid show that the Y axis (parallel with the foliation and veins) is the direction of highest permeability in the foliated sandstone (10-2 md for Y against 10-3 md for X, Z, and the protolith, measured at a confining pressure of 20 bars). Microstructural observations document the localization of the preferential fluid path between the phyllosilicate particles forming the foliation. Hence, the direction of highest permeability in these fault rocks would be parallel with the fault and subhorizontal, that is, perpendicular to the slickenlines representing the local slip direction on the fault surface. We suggest that a similar relationship between kinematic markers and fault rock permeability anisotropy may be found in other fault zone types (reverse or strike-slip) affecting feldspar-rich lithologies in deep diagenesis conditions.
ISSN: 0149-1423

hal-00903428
https://hal.archives-ouvertes.fr/hal-00903428
DOI : 10.1306/03071312127

Résumé : Cette thèse à pour objet l'étude théorique et numérique de solutions dans les équations aux dérivées partielles non linéaires de la physique, en particulier en dynamique des fluides. La présence de discontinuités dans les solutions de ces équations complique la compréhension mathématique des phénomènes mis enjeu et leur traitement numérique, notamment en vue de simulations informatiques . Nous étudions ces équations par une méthode de régularisation dans un espace fonctionnel approprié. Lorsque des schémas numériques construits par des méthodes différentes conduisent à des résultats identiques, ceci jusque dans leurs moindres détails, il semble alors naturel de s'interroger dans quelle mesure ces suites de solutions numériques constituent une approximation d'une solution des équations étudiées. Nous construisons des suites de solutions approchées à partir d'un schéma numérique original,stable et suffisamment simple pour démontrer que ses suites constituent une méthode asymptotique de Maslov au sens des distributions en dimension trois d'espèce. La technique de régularisation employée consiste à étendre les variables réelles du problème ne des variables complexes, ce qui nous permet de construire des familles de solutions particulières que l'on ramène au cas réel en faisant tendre un petit paramètre vers O. Les solutions physiques recherchées apparaissent alors comme valeurs au bord de fonction holomorphes. Nous illustrons les résultats obtenus par des applications en cosmologie dans les cadres Newtoniens et relativistes pour des systèmes sans pression, puis avec pression et auto-gravitation, ainsi que pour le système des gaz parfaits.
This thesis is devoted to the theoretical and numerical study of singular solutions appearing in nonlinear partial differential complicates the mathematical understanding of the phenomena under concem as well as their numerical treatment, in particular in view of computation. These equations are studied by a regularization method in an appropriate functional space. When completely different numerical methods give the same results up to the smallest details one can reasonably expect that these numerical results suggest the existence of a mathematical solution of theses equations. We construct sequences of approximate solutions from an original numerical scheme, which is stable and simple enough to prove that these sequences constitute a Maslov asymptotic method in three space dimension. The regularization technique in use consits in extending the real variables of the problem into complex ones, which perrnits to construct families of particular equations that we bring back to the real case by letting a small paramater tend to zero. The expected physical solutions appear as boundary values of holomorphie functions . Illustrations are given by applications to cosmology in the Newtorian and re1ativistic settings for pressure1ess fluid dynamics, then in presence of self-gravitation and pressure as weil as for the systemof ideal gases

http://www.theses.fr/2011AGUY0478/document

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

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

Éditeur(s) : HAL CCSD 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 Elsevier
Résumé : International audience
The kinetics of hydrogenation of dry iron-bearing olivine single crystals was determined by performing hydration experiments under hydrothermal conditions at high pressure. The experiments were performed in a multi-anvil press at 3 GPa, for a temperature range between 900 and 1200 °C and for various durations. The oxygen fugacity was buffered along Ni-NiO joint. Polarized Fourier transform infrared spectroscopy and recent empirical calibration were used to quantify the hydroxyl distributions in the samples along crystallographic axes after the experiments. The chemical diffusion coefficients are similar (barely slower) than in olivine hydrated at lower pressure (0.2 and 0.3 GPa) for the same diffusion mechanism. Under the given experimental conditions, the anisotropy of diffusion is the same as for proton-vacancy mechanism, with diffusion along the [0 0 1] axis faster than along the [1 0 0]. However, the anisotropy at 3 GPa is weaker compared to measurements at lower pressures and the analysis of concentration profiles using 3D models shows that an isotropic solution could also be relevant. Fits of the diffusion data to an Arrhenius law yield activation energies for the slightly faster [0 0 1] axis of the crystallographic axes around 198 ± 5 kJ mol−1, a value only slightly lower than the results from previous experimental studies for natural iron-bearing olivine hydrogenated at lower confining pressure. At 3 GPa, hydrogenation can be well approximated by a single mechanism controlled by coupled diffusion of protons and octahedral vacancies (di- and tri-valent ions). The diffusion rates are fast enough to alter hydrogen concentration within olivine in xenoliths ascending from the mantle or experiencing hydrogen-rich metasomatism events, but too slow to permit complete homogenization of hydrogen in olivine-rich rocks at kilometer scale in less than one My.
ISSN: 0031-9201

hal-01412663
https://hal.archives-ouvertes.fr/hal-01412663
DOI : 10.1016/j.pepi.2016.08.005

Éditeur(s) : HAL CCSD
Résumé : Thrusting fault zone in foreland basins are characterized by highly foliated zones generally enriched in phyllosilicates which can play a major role on the mechanical behaviour of the fault. In this context, investigations of synkinematic clay minerals permit to determine the origin of the fluid from which they precipitated as well as the mechanisms of deformation. Our study is focused on clay mineral assemblages (illite and chlorite) in a major thrust fault located in the Monte Perdido massif (southern Pyrenees), a shallow thrust that affects upper cretaceous-paleocene platform carbonates and lower Eocene marls and turbidites. It implied 3 km of displacement of the Monte Perdido thrust unit with respect to the underlying Gavarnie unit. In this area the cleavage development by pressure-solution is linked to the Monte Perdido and Gavarnie thrust activity. The core zone of the fault, about 6 m thick, consists of an interval of intensely deformed clay-bearing rocks bounded by major shear surfaces. The deformed sediment is markedly darker than the protolith. Calcite-quartz shear veins along the shear planes are abundant. Detailed SEM and TEM observations of highly deformed fault zone samples indicate that clay mineral enrichment in the core zone of the fault is not only related to passive increase by pressure-solution mechanism but that dissolution-recrystallization of phyllosilicates occurs during deformation. A mineral segregation is observed in the highly deformed zone. Newly formed 2M 1 muscovite is present along the cleavage whereas IIb chlorite crystals fill SV2 shear veins suggesting syntectonic growth of phyllosilicates in the presence of fluids in low-grade metamorphic conditions. These mineralogical reactions act as weakening processes and would favour Monte Perdido fault creeping.
Swiss Journal of Geosciences

hal-00795553
https://hal.archives-ouvertes.fr/hal-00795553
DOI : 10.1007/s00015-012-0098-0

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We describe an experimental system including monitoring of temperature, pressure, pH, oxidation reduction potential and optical density at 600 nm, designed for studying the role of microorganisms on the geological sequestration Of CO2 and its transformation into solid carbonate phases. Measurements were performed in an artificial ground water (AGW) supplemented with urea (2 g.l(-1)) and equilibrated at controlled temperatures with a gaseous phase before bacterial inoculation. We used the ureolytic strain Bacillus pasteurii as a model carbonate precipitating bacteria and showed that it can successfully promote strong pH increases by ureolysis in the AGW equilibrated with CO2 pressures of up to I bar. Increasing salinities (5.8,13.5 and 35.0 g.l(-1)) have a positive effect on the rate of pH increase, whereas the effect of increasing temperatures (30,35 and 38 degrees C) is less important. Calcium is also shown to have a specific positive influence on the rate of ureolysis. The number of viable cells present in solution decreases greatly during the carbonate precipitation event but the population partially recovers once precipitation is over.
ISSN: 0009-2541

hal-00445282
https://hal.archives-ouvertes.fr/hal-00445282
DOI : 10.1016/j.chemgeo.2008.12.012

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
People are exposed to multiple pollutants, especially indoors, in particular through ingestion of indoor settled dust. In the perspective of a cumulative risk assessment, a multi-residue analytical method based on pressurized liquid extraction (PLE) and gas chromatography/tandem mass spectrometry (GC/MS/MS) was developed for the simultaneous analysis in indoor dust of several classes of semi-volatile organic compounds (SVOCs) of health concern, from trace to highly concentrated compounds, including musk fragrances, organochlorines (OCs), organophosphates (OPs), polycyclic aromatic hydrocarbons (PAHs), polybromodiphenylethers (PBDEs), polychlorobiphenyls (PCBs), phthalates and pyrethroids. The method was validated in terms of limits of quantification (LOQ), and accuracy and precision via spiking experiments on an inert material (Celite(®) 545) and replicate analysis of the standard reference material SRM 2585 supplied by the National Institute of Standards and Technology (NIST). Method LOQs for 200 mg samples of sieved dust were 26 ng g(-1) for PCBs and some OCs, 65 ng g(-1) for musks, OPs, PAHs, PBDEs, pyrethroids, other OCs and some phthalates, 132 ng g(-1) for butylbenzylphthalate (BBP), 197 ng g(-1) for tributylphosphate and 1579 ng g(-1) for other phthalates. Quadratic calibration curves were established for each compound by analyzing at least five calibration solutions and exhibited determination coefficients higher than 0.999. The method was successfully applied to the SRM 2585 and seven real indoor dust samples. The results obtained on SRM 2585 demonstrate both excellent reproducibility and agreement with the indicative, reference or certified values and provide, for the first time, indicative concentrations for chlorpyrifos, diazinon, diisononylphthalate (DiNP) and tetramethrin. The results obtained on real dust samples illustrate the ability of the proposed method to quantify a wide range of SVOCs in a single analysis, making it appropriate for environmental monitoring programs or large-scale studies with a large number of samples.
ISSN: 0021-9673

hal-00994493
https://hal-univ-rennes1.archives-ouvertes.fr/hal-00994493
DOI : 10.1016/j.chroma.2014.02.004
PUBMED : 24598454

É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 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 Wiley
Résumé : International audience
In this paper, the generalized Stokes problem αu−νΔu+∇p=f, ∇u=0, is considered, where ν denotes the kinematic viscosity, α is a positive parameter proportional to the inverse of the time-step, and the source term f gathers the external forcing, the nonlinear term, and the explicit part of the linear term. The system is derived from a time discretization for the Navier-Stokes equation governing the motion of a two-dimensional incompressible fluid. A domain decomposition numerical method is developed for solving the system. The original domain is divided into several nonoverlapping subdomains, and in each subdomain the concerned equations are discretized by certain finite element methods. A new feature of the method is that the pressure p is required to be piecewise linear rather than piecewise constant, so that the nonphysical oscillations in the solution using piecewise linear pressure are eliminated. However, a difficulty in using piecewise linear pressure is that the Stokes problem in each subdomain involves Dirichlet boundary conditions on pressure at the interface, which is not a real Stokes problem. To avoid the difficulty, a modification is proposed which removes the continuity conditions on the pressure across the interface. In this way, the Stokes problem in each subdomain has only Dirichlet boundary conditions on the velocity. Formulations and discussions for both the methods are presented, which shows that the modified method is simpler than the original one. Numerical experiments are displayed to show the efficiency of the methods, especially that the modified method has the same accuracy as the original one.
ISSN: 0749-159X

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

Éditeur(s) : Université des Antilles Études caribéennes
Résumé : The tropical coast of Queensland, Australia is experiencing rapid population and tourism growth. Heavy development pressure is being placed on the biophysical environment, which includes the largest sand island in the world (World Heritage listed), Fraser Island. Despite the existence of State planning legislation oriented towards environmental sustainability, it is apparent that little, if any, regard is being paid to the long-term effects of resort and canal-estate development on the underlying biophysical environment of the coast and islands. While the remarkable surface features of the tropical coastal environment such as the surf beaches, the unique dune lakes, and the coastal vegetation, as well as the prolific wildlife, are used as selling points for new settlement and above all tourism, little is known or apparently of major concern with respect to the long term environmental impact of coastal development. While the crowding effect of large numbers of visitors and their impact on the natural environment - through significant degradation, site hardening or lack of proper on-site management – has been identified as a particular problem in certain areas, there is little knowledge of impacts on the underlying coastal geomorphology. In particular, developmental choices are often made in ignorance of, or disregard of, the widespread existence of tropical acid sulfate soils in that coastal environment. Reactive soils such as these can, and do, greatly affect the biophysical outcomes of a development if they are exposed. This paper documents such impacts in the form of a case study of the acid sulfate soil problem in relation to the development of tourist resorts, and suggests that part of the solution is for the State Government to insist on enforceable standards relating to their development for tourism in the coastal zone. These would reduce the use by many resort owners of fancy golf course developments as an environmental cover up to make it look as if they put a lot of effort into protecting the environment….
La côte tropicale du Queensland est confrontée à un accroissement rapide de la population. Les fortes pressions anthropiques pèsent sur le milieu biophysique, qui comprend notamment a plus grande île de sable du monde, un site classé Patrimoine de l’Humanité, Port Hinchinbrook et Fraser Island. Malgré l’existence d’une législation nationale orientée vers la préservation de l’environnement, il est manifeste que peu d’attention est accordée aux effets à long terme des sites touristiques et du développement de l’urbanisation sur ces franges côtières sensibles. L’aspect massif des flux de visiteurs et leurs impacts sur le milieu naturel a été identifié comme un sujet d’inquiétude particulier dans certaines régions. La connaissance des impacts sur la géomorphologie côtière souterraine est très réduite et les choix de développement sont souvent faits dans l’ignorance ou au mépris de la présence répandue de terrains de sulfate d’acide dans ce milieu côtier. Des sols réactifs comme ceux-là peuvent affecter, et affectent, grandement les résultats géophysiques d’un développement lorsqu’ils sont connus. Des études de cas sur le problème des sols de sulfate d’acide suggèrent qu’une partie de la solution pour le gouvernement de l’Etat du Queensland est d’imposer des normes obligatoires pour le développement du tourisme dans la zone côtière.
Australie
Queensland

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

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In orogenic systems, thrust faults play a major role in stacking different tectonic units and may act as conduits for the expulsion of large amounts of fluid of different origins (metamorphic, diagenetic, meteoric). This study focuses on the Monte Perdido thrust unit emplaced in the Paleogene Jaca thrust-sheet-top basin, in the SW-central Pyrenees. We aim to decipher the mechanisms and P-T conditions of deformation in fault zones and characterize the related fluid involvement, through combined microstructural, geochemical and microthermometry analyses. Two thrust faults cutting platform limestones, marls and siliciclastic turbidites of the lower part of the basin-fill (Paleocene-lower Eocene) have been studied. The fault zones are characterized by metre-thick shear zones with highly deformed, foliated clay-rich sediments. Foliation is underlined by preferentially oriented phyllosilicates. Several generations of shear and extension calcite, quartz and chlorite-bearing veins attest to fluid-rock interactions during a multi-stage deformation. Microstructural observations and stable isotope analyses on calcite from veins and host sediments suggest that deformation was aseismic and dominated by diffusive mass transfer from pressure solution sites along cleavage and stylolites to the precipitation sites in veins, with mineralizing fluids in equilibrium with the host sediments. Our results suggest an essentially closed hydrologic system, and imply the absence of significant fluid flow along the studied fault zones. Microthermometric study on fluid inclusions present in calcite and quartz veins, and calcite-quartz oxygen isotopic fractionation determined for the first generation shear veins, allow a geothermal gradient of 34 degrees C/km to be estimated. Analytical results demonstrate an evolution of the fault zones in three stages. The first stage was related to the emplacement of the Monte Perdido thrust unit during the middle Eocene at a temperature of similar to 208 degrees C and a burial depth of similar to 5.7 km. The second stage corresponds to a fault reactivation at a temperature of similar to 240 degrees C and a burial depth of similar to 6.5 km. The latter deformation may have been related to folding of the Monte Perdido thrust unit during the emplacement of the underlying Gavarnie thrust unit during the late Eocene-early Oligocene, with deeper burial resulting from aggradation of the thrust-sheet-top basin-fill. The last event corresponds to the formation of a dilatant vein system likely related to the exhumation of the massif.
ISSN: 0191-8141

hal-00644828
https://hal.archives-ouvertes.fr/hal-00644828
DOI : 10.1016/j.jsg.2011.06.006

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Despite many studies, our understanding of silica precipitation from natural waters remains limited, in particular for geothermal waters. Here we present a detailed study on the kinetics of silica polymerization as a function of fluid temperature and pH using high-temperature (250-260 °C) seawater-derived geothermal fluids as those discharged from the Bouillante geothermal site.We monitored the on-site decrease in monomeric silica concentration(initial SiO2 concentration of about 600 mg/l)with time using the molybdenum blue spectrophotometric method on samples of separated water collected from the high-pressure separator at 167 °C and cooled to 25, 50, 75, and 90 °C and for pH values ranging from 5 to 8. During all these experiments, only silica polymerization was observed, with the formation of colloidal particles in suspension in the solutions. No scaling of amorphous silica was formed. The collected data were after modeled in order to determine the useful kinetic parameters for predicting and preventing amorphous silica precipitation in the production wells during fluid exploitation in the specific context of Bouillante. Results under the investigated experimental conditions show that the kinetics of silica precipitation is affected more strongly by pH than by temperature change. The reaction in the acidic Bouillante solution begins with a transition period that significantly decreases the kinetics of silica polymerization. Modeling the experimental data indicates that the silica polymerization up to the state of equilibrium is characterized by a 2nd-order kinetic law relative to dissolved silica; this could indicate that the polymerization is controlled mainly by the formation of dimers. However, the first hour of the experiments is better characterized by a 4th-order kinetic law, suggesting more complex polymerization reactions in the initial stage with the formation of nanocolloidal particles containing 3 to 4 monomers. In both cases, the corresponding kinetics rate constant k is linearly dependent on pH for pH between 5 and 7. The activation energy Ea for the overall reaction, calculated using the Arrhenius equation under the considered pH and temperature conditions, ranges between 41±9.8 and 54±9.6 kJ/mol. To complete this work, the colloid particles formed at the end of the kinetics experiments were extracted and analyzed by SEM and TEM microscopy, X-ray Diffraction and the BET method. Results show that the size of the colloids increase and their specific surface decrease with increasing pH and are thus dependent on pH. Globally, our work provides a reliable database for understanding silica polymerization kinetics in natural geothermal brines or geologic waters characterized by a near neutral pH and moderate dissolved silica concentrations.
ISSN: 0009-2541

hal-01403092
https://hal-brgm.archives-ouvertes.fr/hal-01403092
DOI : 10.1016/j.chemgeo.2016.08.031

Éditeur(s) : HAL CCSD Wiley-Blackwell
Résumé : This study describes normal fault zones formed in foreland arkosic turbidites (the Grès d'Annot Formation, SW French Alps) under deep diagenesis conditions (~200 °C) and highlights the occurrence of two markedly different fault-rock types: (1) the foliated fault rocks of the Moutière-Restefond area; and (2) the dilatant fault rocks of the Estrop area. The deformation of (1) is dominated by intra- and transgranular fracturing, pressure solution of quartz and feldspar grains and syn-kinematic phyllosilicate precipitation resulting from feldspar alteration. The combination of these mechanisms results in a strongly anisotropic strain with intense shortening normal to the foliation (pressure solution) and extension parallel to the foliation (quartz- and calcite-sealed extension veins). This deformation implies local mass transfer that may be achieved without (or with limited) volume change. The deformation of (2) is expressed as dilatant quartz-sealed veins and breccia textures in which the main mechanisms are transgranular fracturing and quartz precipitation. Type (2) implies fault volume increase, isotropy of deformation and mass transfer at distances larger than in type (1). This study discusses the origins of (1) and (2) and shows that the permeability of (1) is anisotropic, with higher values than the host rocks parallel to the Y main deformation axis (i.e. perpendicular to the slip vector), whereas the permeability of (2) is isotropic and equivalent to that of the host rocks.
ISSN: 0954-4879

hal-01115177
https://hal.archives-ouvertes.fr/hal-01115177
DOI : 10.1111/ter.12100

Éditeur(s) : HAL CCSD Elsevier
Résumé : Le titre a été changé Ancien titre: Roles for local crowding and pressure in counteracting density fluctuations at equilibrium
International audience
What is the probability that all the gas in a box accumulates in the same half of this box? Though amusing, this question underlies the fundamental problem of density fluctuations at equilibrium, which has profound implementations in many physical fields. The currently accepted solutions are derived from the studies of Brownian motion by Smoluchowski, but they are not appropriate for the directly colliding particles of gases. Two alternative theories are proposed here using self-regulatory Bernoulli distributions. A discretization of space is first introduced to develop a mechanism of matter congestion holding for high densities. In a second mechanism valid in ordinary conditions, the influence of local pressure on the location of every particle is examined using classical laws of ideal gases. This approach reveals that a negative feedback results from the reciprocal influences between individual particles and the population of particles, which strongly reduces the probability of atypical microstates. Finally, a thermodynamic quantum of time is defined to compare the recurrence times of improbable macrostates predicted through these different approaches.
Physica A: Statistical Mechanics and its Applications

hal-00920985
https://hal-univ-rennes1.archives-ouvertes.fr/hal-00920985
https://hal-univ-rennes1.archives-ouvertes.fr/hal-00920985v2/document
ARXIV : 1404.6223
DOI : 10.1016/j.physa.2014.03.067

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In this paper we solve the Navier Stokes’ equation using finite difference method, on a simulated porous rock structure in 2-D, 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. Depending on the rock composition, both dissolution and precipitation reactions may occur. However precipitation occurs only through the cations that are released in the solution due to dissolution. These combined dissolution–precipitation reactions change the porosity, permeability and pore geometry of the sedimentary rock. We study the temporal changes of these properties as functions of Peclet number, concentration of the reactive solute and ratio of Damkholer numbers of dissolution to precipitation. The final flow property is decided by a combination of these parameters.
ISSN: 0022-1694

hal-01115291
https://hal.archives-ouvertes.fr/hal-01115291
DOI : 10.1016/j.jhydrol.2014.10.019

Éditeur(s) : American Geophysical Union
Résumé : Two time series of deep ocean bottom pressure records (BPRs) in between the Crozet Islands and Kerguelen are compared with GRACE (Gravity Recovery And Climate Experiment) equivalent water heights. An analysis of the correlation is performed for four time series: 1) monthly averages of the equivalent water height at the Crozet Islands, 2) the same near the Kerguelen Islands, 3) the mean of the two preceding series and 4) the difference between the two locations expressed in terms of geostrophic transport. We find that smoothed GRACE solutions are strongly correlated with the BPR data with correlation coefficients in the order of 0.7-0.8. Consequently GRACE measures real oceanic mass variations in this region.
Geophysical Research Letters ( GRL ) (0094-8276) (American Geophysical Union), 2006-11 , Vol. 33 , N. 21 , P. NIL_10-NIL_13

Droits : 2006 AGU
http://archimer.ifremer.fr/doc/2006/publication-2088.pdf
DOI:10.1029/2006GL027452
http://archimer.ifremer.fr/doc/00000/2088/

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
An analytical solution in the Laplace domain is derived for modeling anomalous pressure diffusion during pumping tests in aquifer displaying hierarchical fractal fracture networks. The proposed solution generalizes all of the analytical models for fractal flow published previously by combining multifractal diffusion and nested multiporosity with transient exchanges, interface skin effects, and well storage effect. Solutions are derived for fracture-delimited blocks with planar, cylindrical, and spherical shapes, as well as with any fractional dimensional shapes. Any combinations of these shapes can be defined in order to model a large range of situations. Within each permeability level, the fractal properties of the fracture network can be specified and the fractal dimension can be distinct from the shape dimension of the block.
ISSN: 0043-1397

hal-00411549
https://hal.archives-ouvertes.fr/hal-00411549
DOI : 10.1029/2008WR007125