Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : We designed a thermo-mechanical numerical model for fast-spreading mid-ocean ridge with variable viscosity, hydrothermal cooling, latent heat release, sheeted dyke layer, and variable melt intrusion possibilities. The model allows for modulating several accretion possibilities such as the "gabbro glacier" (G), the "sheeted sills" (S) or the "mixed shallow and MTZ lenses" (M). These three crustal accretion modes have been explored assuming viscosity contrasts of 2 to 3 orders of magnitude between strong and weak phases and various hydrothermal cooling conditions depending on the cracking temperatures value. Mass conservation (stream-function), momentum (vorticity) and temperature equations are solved in 2-D cartesian geometry using 2-D, alternate direction, implicit and semi-implicit finite-difference scheme. In a first step, an Eulerian approach is used solving iteratively the motion and temperature equations until reaching steady states. With this procedure, the temperature patterns and motions that are obtained for the various crustal intrusion modes and hydrothermal cooling hypotheses display significant differences near the mid-ocean ridge axis. In a second step, a Lagrangian approach is used, recording the thermal histories and cooling rates of tracers travelling from the ridge axis to their final emplacements in the crust far from the mid-ocean ridge axis. The results show that the tracer's thermal histories are depending on the temperature patterns and the crustal accretion modes near the mid-ocean ridge axis. The instantaneous cooling rates obtained from these thermal histories betray these discrepancies and might therefore be used to characterize the crustal accretion mode at the ridge axis. These deciphering effects are even more pronounced if we consider the average cooling rates occurring over a prescribed temperature range. Two situations were tested at 1275-1125 °C and 1050-850 °C. The first temperature range covers mainly the crystallization range that is characteristic of the high temperature areas in the model (i.e. the near-mid-oceanic-ridge axis). The second temperature range corresponds to areas in the model where the motion is mainly laminar and the vertical temperature profiles are closer to conductive. Thus, this situation results in less discriminating efficiency among the crustal accretion modes since the thermal and dynamic properties that are described are common to all the crustal accretion modes far from the ridge axis. The results show that numerical modeling of thermo-mechanical properties of the lower crusts may bring useful information to characterize the ridge accretion structure, hydrothermal cooling and thermal state at the fast-spreading ridges and may open discussions with petrological cooling rate results.
ISSN: 1991-959X

hal-00907012
https://hal.archives-ouvertes.fr/hal-00907012
DOI : 10.5194/gmd-6-1659-2013

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
The Alps/Apennines system, as well as many collisional orogens through the world, shows a finite deformation produced during a long geological history which involves numerous superimposed tectonic events. As a result, complex and often contrasted reconstructions for the setting and tectonics of the different stages of the growing and interfering Alps/Apennines system have been proposed. To enlight some of the geometric and kinematic signatures related to past geodynamics in the composite Alps(Corsica)/Apennines system, we analyse the major deformation processes, their superimposition, and resulting finite structural geometries in the currently growing Taiwan orogenic system. The Taiwan island and its offshore northward and southward prolongations represent the emerged and the sub-marine parts of an orogenic wedge related to the interactions between the Eurasian and the Philippine Sea and recording different subduction settings. Using the Taiwan actualistic geodynamic scenarios and the related first-order geometric-kinematic constraints, we will attempt a reconstruction of the tectonic evolution for the Corsica/Apennines orogenic system foregrounding the key role of continental subduction and subduction reversal during such evolution.
ISSN: 1437-3254

hal-00617149
https://hal.archives-ouvertes.fr/hal-00617149
DOI : 10.1007/s00531-010-0598-y

Éditeur(s) : HAL CCSD Institute for Scientific Computing and Information
Résumé : International audience
The production of oil in soft highly compacted reservoirs produces a significant reduction in the pore space which in turn causes the surface to subside. The authors couple two simulators to model the mechanical deformations in the presence of a Darcy flow, both one- and two-phase. The simulators solve for flow in a porous medium using a finite volume method, and solve for the mechanical deformation of the porous media using a finite element method. The flow equations and the deformation equations are coupled using Biot's law. The authors use two different methods for the coupling: the first uses a staggered algorithm; the second is based on a preconditioned conjugate method. The authors use a multiscale method in both approaches. The authors show that the inf-sup condition is not satisfied, because of the mixed discretization which is done by using two separate simulators. The authors give convergence comparisons for the two methods for several examples. Two of the examples are for the linear problem (one-phase) where they know the exact solution and one is for the nonlinear problem (two-phase). In all cases, the preconditioned conjugate algorithm is faster and more robust, in particular when the inf-sup condition is not satisfied. In the nonlinear problem the authors state that the results agree with physical predictions but do not give actual empirical data. Finally, they show the effects of grid orientation on the nonlinear problem.
ISSN: 1705-5105

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

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : International audience
The hydrogen content of nominally anhydrous minerals is of great interest, because it can influence many physical and mechanical properties of mantle rocks. Moreover, the hydrogen diffusion profiles can be used to constrain timescales related to magma eruptions. Here, we report models of ionic diffusion for trace elements in anisotropic crystals and apply them to hydrogen diffusing out of mantle-derived olivine. We first compare and discuss the characteristics of 1D and 3D models and show that only 3D anisotropic diffusion models can lead to diffusion profiles exhibiting non-equilibrium plateau at the center of the solid along the slowest axis, as measured in natural samples. In a second part, we discuss the differences between hydration and dehydration of olivine for diffusion that is linked to two different atomic sites involved in hydrogen mobility. Finally, we apply our 3D anisotropic model to previous results on mantle-derived olivine from Pali-aike to better characterize diffusion coefficients and their anisotropy that could be relevant for dehydration of olivine. Our results show that dehydration has to be strongly anisotropic, with a fast [100] axis and a significantly slower [001] axis.
ISSN: 0342-1791

hal-01115253
https://hal.archives-ouvertes.fr/hal-01115253
DOI : 10.1007/s00269-014-0685-x

Éditeur(s) : HAL CCSD Global Science Press
Résumé : International audience
In this paper, the authors present two different time-discretization schemes combined with a finite-volume space discretization for the incompressible Navier-Stokes equations. They show that the collocated space discretization, usually used with a projection method, can also be extended to other time discretizations over previous projection methods, like the splitting methods. The main advantage of this splitting method is that the pressure approximation is truly of second order. In the first and second parts of this article, two different time-discretization schemes are considered with a collocated space discretization and it is explained how the unknown can be correctly coupled. Numerical simulations are presented in the last part. Tests of spatial and time accuracies are made. The numerical results when the Navier-Stokes problem is solved are shown, in the model problem of the driven cavity.
ISSN: 1815-2406

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

É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 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 Elsevier Masson
Résumé : International audience
The seismological study of recent seismic crises near Oleron Island confirms the coexistence of an extensional deformation and a transtensive regime in the Atlantic margin of France, which is different from the general western European stress field corresponding to a strike-slip regime. We argue that the switch of the principal stress axes σ1/σ2 in a NW-SE vertical plane is linked with the existence of crustal heterogeneities. Events of magnitude larger than 5 sometimes occur along the Atlantic margin of France, such as the 7 September 1972 (ML = 5.2) earthquake near Oleron island and the 30 September 2002 (ML = 5.7) Hennebont event in Brittany. To test the mechanism of local strain localization, we model the deformation of the hypocentral area of the Hennebont earthquake using a 3D thermo-mechanical finite element code. We conclude that the occurrence of moderate earthquakes located in limited parts of the Hercynian shear zones (as the often reactivated swarms near Oleron) could be due to local reactivation of pre-existing faults. These sporadic seismic ruptures are favoured by stress concentration due to rheological heterogeneities.
ISSN: 1631-0713

hal-00936220
https://hal.archives-ouvertes.fr/hal-00936220
DOI : 10.1016/j.crte.2013.07.004

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Many iterative processes can be interpreted as discrete dynamical systems and, in certain cases, they correspond to a time discretisation of differential systems. In this article the authors propose generating numerical methods in numerical linear algebra by modelling the linear system to be solved as a given state of a dynamical system; the solution can be reached asymptotically, as a (asymptotically stable) steady state, but also as a finite time (shooting methods). In that way, any (stable) numerical scheme for the integration of such a problem can be presented as a method for solving linear systems. The authors discuss aspects of this approach, which allows them to recover some known methods but also to introduce new ones. Finally, some convergence results and numerical illustrations are presented.
ISSN: 1017-1398

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

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

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

Éditeur(s) : HAL CCSD
Résumé : Aricle soumis à Solar Energy
In order to qualify the fluctuating nature of solar radiation under tropical climate, we classify daily distributions of the clearness index kt by estimating a finite mixture of Dirichlet distributions without assuming any parametric hypothesis on these daily distributions. The method is applied to solar radiation measurements performed in Guadeloupe (16°2 N, 61 W) where important fluctuations can be observed even within a short period of a few minutes. The results put in evidence four distinct classes of distributions corresponding to different types of days. The sequence of such classes can be of interest for prediction.
https://hal.archives-ouvertes.fr/hal-00289567

hal-00289567
https://hal.archives-ouvertes.fr/hal-00289567
https://hal.archives-ouvertes.fr/hal-00289567/document
https://hal.archives-ouvertes.fr/hal-00289567/file/Article_Classif_soumis.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In the Chalkidiki Peninsula of northern Greece a thrust complex made of a basement (Vertiskos Unit), a cover (Circum–Rhodope belt) and arc/back-arc units (Chortiatis Magmatic Suite and eastern Vardar Ophiolites) is exposed in the Chalkidiki Peninsula of northern Greece. The complex forms the western part of the Rhodope Metamorphic Province and lies on the hanging-wall of the Kerdylion Detachment, the structure responsible for the exhumation of the Southern Rhodope Core Complex and the most prominent and visible ductile structure related to the Tertiary Aegean extension. The Chalkidiki thrust complex arguably preserves a complete record of Cretaceous deformation and related fabrics. In this contribution we describe the geometry of foliation, stretching lineation and shear sense(s) on a regional scale. The regional foliation strikes NW–SE and displays different patterns in the three studied units: (i) dominantly dipping at low angle in the Vertiskos Unit, (ii) affected by upright folding in the Circum–Rhodope belt and (iii) systematically steeply dipping to the NE in the Chortiatis Magmatic Suite. Stretching lineation trend dominantly SW–NE in the three mentioned units. On the basis of new mapping, neglecting local perturbations and deformation related to Tertiary extension, we infer the regional kinematics of Cretaceous syn-metamorphic thrusting and subsequent exhumation of the metamorphic units. Thrusting took place toward the SW (in present-day coordination) and the related fabrics are recorded throughout the metamorphic pile. On the contrary, exhumation-related fabrics are related to shear toward the NE and are preferentially recorded in the uppermost part of the metamorphic pile suggesting that extension was more localised and of less finite intensity compared to thrusting.
ISSN: 0040-1951

insu-01231485
https://hal-insu.archives-ouvertes.fr/insu-01231485
DOI : 10.1016/j.tecto.2015.09.034

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This article is devoted to the presentation of a multilevel method using finite differences that is well adapted for solving Stokes and Navier-Stokes problems in primitive variables. We use Uzawa type algorithms to solve the saddle point problems arising from spatial discretization by staggered grids and a semi-explicit temporal scheme. By means of a new change of basis operator, the two-dimensional velocity and pressure fields of an M.A.C mesh are gathered in a hierarchical order, into several grids preserving the M.A.C property on each of them. These new hierarchical unknowns, called Staggered Incremental Unknowns (SIU), allow us to develop techniques which reduce the cost of the resolution of either Stokes or generalized Stokes problems. An experimental estimation of the condition number of the inner matrix is given, and justifies the preconditioning effect of the staggered incremental unknowns.
ISSN: 0168-9274

hal-00770268
https://hal.univ-antilles.fr/hal-00770268
DOI : 10.1016/S0168-9274(99)00044-6

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In order to qualify the fluctuating nature of solar radiation under tropical climate, we classify daily distributions of the clearness index kt by estimating a finite mixture of Dirichlet distributions without assuming any parametric hypothesis on these daily distributions. The method is applied to solar radiation measurements performed in Guadeloupe (16°2 N, 61 W) where important fluctuations can be observed even within a short period of a few minutes. The results put in evidence four distinct classes of distributions corresponding to different types of days. The sequence of such classes can be of interest for prediction.
ISSN: 0038-092X

hal-01527878
https://hal.archives-ouvertes.fr/hal-01527878
DOI : 10.1016/j.solener.2009.01.010

Éditeur(s) : HAL CCSD
Résumé : International audience
Two opposing visions of the Aegean backarc tectonics implicitly contain contrasting images of the rheological behaviour of the lihosphere, supported by different sets of observations. The propagation of the NAF and extension in the Corinth Rift suggest a strongly localising rheology whereas the formation of MCC in the Cyclades and the Rhodope suggests instead a more viscous behaviour. This paradox was the seminal question addressed by the ANR-EGEO programme (2007-2010) that leads to the following conclusions: (1) The exhumation of the Cycladic MCC's is accommodated by the N-dipping North Cycladic Detachment System (NCDS) that partly reworks the Vardar suture and at the base by a series of thrusts, including the basal contact of the Cycladic Blueschists (CBS) over the Cycladic Basement (CB). The activity of LANF is due to the reactivation of pre-existing discontinuities such as thrusts or earlier detachments and much less to the interaction with granitic plutons. (2) Exhumation in the Cyclades has proceeded in two stages: (a) Eocene syn-orogenic exhumation within the subduction lower/upperplate interface while the post-orogenic Rhodope MCC formed further north, (b) Oligo-Miocene post-orogenic extension in an MCC mode coeval with syn-orogenic exhumation in Crete and the Peloponnese. (3) Both stages were associated with slab retreat, as early as the Eocene with an acceleration at 30-35 Ma. (4) The localisation of the presently active steeply-dipping normal faults on the southern margin of the Corinth Rift may have been preceded by a partial reworking of thrusts and syn-orogenic detachments by shallow-dipping decollements. (5) The localisation of the NAF in the Marmara Sea region and the Northern Aegean Sea is now accurately dated between 5.3 and 5 Ma Ma using the erosion and deposition surfaces that mark the Messinian salinity Crisis and nannofossils ages. (6) The formation of MCCs in the Cyclades resulted in the draining of the low viscosity lower crust from the northern Aegean Sea and led to an easier coupling between upper mantle and upper crust, easing the localisation and propagation of the NAF. (7) A comparison between crustal and mantle finite strain suggests that mantle flow due to slab retreat controls the stretching of the crust. Our model involves progressive slab retreat and mantle flow below the overriding plate inducing the observed succession of tectonic episodes with a progressive localisation of extension due to several slab tearing events. There is no paradox in the tectonic history and rheological behaviour of the Aegean lithosphere but an evolution of the kinematic boundary conditions due to slab deformation at depth.
American Geophysical Union
San Francisco, United States

insu-00860518
https://hal-insu.archives-ouvertes.fr/insu-00860518
BIBCODE : 2010AGUFM.T13G..04J

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

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

Résumé : L’objet de cette thèse est l’étude des comportements mécanique et à la corrosion d’un nouveau système de fondation, dit pieux SPIRMEC, adapté aux structures légères. Le travail s’est focalisé sur deux grands aspects. Le premier aspect concerne l’analyse du comportement sous chargement axiale des pieux SPIRMEC sur la base de l’expérimentation et de la modélisation numérique dans un sol mou non cohésif. Quant au deuxième aspect, il traite du problème de la corrosion de l’acier constituant le pieu dans un environnement tropical humide. Le comportement des pieux SPIRMEC sous une charge axiale de traction, dans un sol fin non cohésif, a montré le développement d’une surface cylindrique de rupture. Celle-ci ne dépend pas de la configuration géométrique des pieux tel que : l’espacement S, le diamètre d’hélice Dh et le diamètre du tube d. En adoptant ce critère de rupture et à l’aide d’une analyse statistique des données, nous avons proposé quatre méthodes analytiques de prédiction de la capacité portante des pieux SPIRMEC à partir des essais in-situ et au laboratoire qui sont: •Méthode CPT basée sur la résistance de pointe pénétromètrique (essai CPT), •Méthode MPT basée sur la pression limite pressiomètrique (essai MPT), •Méthode du couple d’installation, •Méthode analytique basée sur les paramètres de cisaillement de sol notamment l’angle de frottement interne (essais triaxial et boîte de cisaillement). Ces méthodes tiennent compte de la géométrie du pieu, du procédé de mise en oeuvre et de la compacité du sol. L’analyse de performance des méthodes proposées a montré que les méthodes analytique et CPT représentent mieux le comportement des pieux testés. La modélisation numérique du comportement des pieux SPIRMEC sous charge axiale de traction a été également étudiée à l’aide du code de calcul Plaxis 2D. Un modèle axisymétrique a été adopté pour modéliser la géométrie de la structure. Le modèle de comportement hardening soil a été retenu pour modéliser le comportement élastoplastiq des couches de sol. Un modèle élastique linéaire a été retenu pour le pieu. Quant au comportement de l’interface sol-structure, il a été pris en compte par un modèle élastoplastique basé sur les critères de Mohr-Coulomb. L’expansion du sol résultant de l’installation des pieux a été approchée par une cavité cylindrique. Cette dernière a été modélisée par des déformations volumiques anisotropes. Nous avons constaté, qu’appliquer une déformation volumique de l’ordre de 40% permet une bonne concordance entre les résultats numériques et expérimentaux. L’étude paramétrique nous a également permis de constater l’existence d’une variation linéaire entre la capacité portante et l’espacement S. Le comportement de l’acier galvanisé dans le sol a été également étudié au moyen de trois approches: électrochimique, visuelles grâce au microscope électronique et gravimétrique. Les tests électrochimiques de corrosion du zinc ont été réalisés dans la solution de sol artificielle pour deux périodes bien distinctes : saison des pluies et après 8 mois de contact sol-pieu. Une forte dépendance est constatée de la vitesse de corrosion de zinc au pH et à la conductivité de la solution. La vitesse de corrosion est plus importante en saison des pluies qu’après 8 mois de contact de pieu avec le sol. Les fortes concentrations en sulfates semblent augmenter la vitesse de corrosion. Les analyses menées au MEB montrent que la zone de fluctuation de nappe et en surface externe (sol/air) sont les zones les plus corrodées. La comparaison entre les différentes techniques d’évaluation de la vitesse de corrosion a montré que l’eau de nappe peut être utilisée comme un milieu représentatif du profil chimique d’un sol non cohésif et légèrement granulaire.
The purpose of this thesis is to study the both behaviors mechanical and of the corrosion of a new foundation system, called spirec piles, for lightweight structures. The work was based on two main aspects. The first aspect concerns the analysis of behavior of the spirmec piles under axial loading on the basis of experimental and numerical modeling in a non-cohesive soft soil. On the second aspect, it addresses the problem of corrosion of steel forming the pile in a humid tropical environment. The behavior of the tested piles under axial traction load, in a non-cohesive soil, showed the development of a cylindrical surface of failure. This does not depend on the geometrical configuration of the tested piles such that: the spacing s, the helix diameter dh and the diameter of the shaft. By adopting this failure criterion and using statistical data analysis, we proposed four analytical methods for predicting the bearing capacity of the piles from in-situ and laboratory tests, which are: cpt method based on tip resistance penetrometer (cpt test), mpt method based on pressuremeter limit pressure (test mpt), method of installation torque and analytical method based on the shear parameters of soil (triaxial tests and shear box). These methods take into account the geometry of the pile, the method and implementation of soil compaction. Performance analysis of the proposed methods showed that the analytical methods and cpt represent better the behavior of piles tested. An axisymmetric condition was assumed to model this geometry in two-dimensional space. In this study, we proposed to model the pipe screw pile with plates elements available in plaxis 2d. The helixes were modeled as circular disks. Interface elements were incorporated along the pile to simulate the soil-pile interaction. Hardening soil model was chosen to describe the soil with drained behavior was assigned for all layers. Linear elastic for pile and elastic plastic model based on the mohr-coulomb criterion was used to describe the interface behavior. The lateral expansion generated by the pile shaft was modeled by anisotropic volumetric strains (ɛxx=ɛzz, ɛyy=0). The pile behavior under tensile loading at different volumetric strains, it is clear that the best fit was obtained at volumetric strain of 40%. The parametric study also allowed us to establish the existence of a linear variation between the bearing capacity and spacings. The behavior of the galvanized steel in the ground has also been studied using three approaches such as: electrochemical, microscopic imaging and gravimetrically. Electrochemical tests of the corrosion of the zinc have been made in artificial soil solution into two periods, the rainy season and after eight months of contact soil-pile. Is a strong dependence of the observed corrosion rate of the zinc at the ph and conductivity of the solution. The corrosion rate is higher in the rainy season compared to eight months of contact of the pile with soil. High sulfate concentrations increase the rate of corrosion. The meb analysis shows that the fluctuation zone of the water table is the most corroded area. Comparison between the different methods of evaluating the corrosion rate showed that the ground water may be used as a representative profile of the chemical environment of a non-cohesive soil and slightly granular.

http://www.theses.fr/2014AGUY0790/document

É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) : HAL CCSD Springer Verlag
Résumé : International audience
Continuous Global Positioning System (CGPS) position time series from eastern North America constrain the pattern and magnitude of regional crustal deformation. Initial analysis delineates consistent uplift patterns, as expected from glacial isostatic adjustment (GIA) predictions, but the associated horizontal deformation is not definitive, in part due to the short time periods for a significant number of the available stations. We employ an eigenpattern decomposition in order to define a unique, finite set of deformation patterns for this continuous GPS data. Similar in nature to the empirical orthogonal functions historically employed in the analysis of atmospheric and oceanographic phenomena, the method derives the eigenvalues and eigenstates from the diagonalization of the correlation matrix using a Karhunen-Loeve expansion (KLE). The KLE technique is used to identify the important modes in both time and space for the CGPS data, modes that potentially include signals such as horizontal and vertical GIA, tectonic strain, and seasonal effects. Here we filter both the vertical and horizontal velocity patterns on different spatiotemporal scales in order to study the potential geophysical sources, after the removal of correlated and random noise. The method is successful in disaggregating the linear vertical signal from more variable and less spatially correlated signals. The vertical and horizontal results are compared to the predictions of the ICE-3G GIA loading model with a number of plausible mantle viscosity profiles. The horizontal velocity analysis allows for qualitative differentiation between several potential GIA models and suggests that, with longer time series, this technique can be employed to remove correlated noise and improve estimates of crustal strain patterns and their sources.
ISSN: 0033-4553

hal-00795454
https://hal.archives-ouvertes.fr/hal-00795454
DOI : 10.1007/s00024-011-0420-1