L'injection de tensioactifs est une méthode très appliquée dans le domaine de la récupération améliorée des hydrocarbures. Cependant, son efficacité repose sur la capacité de ces agents chimiques à mobiliser l'huile résiduelle en diminuant la tension interfaciale entre l'huile et l'eau. Des modèles à l'échelle du réservoir calculent l'efficacité de la récupération d'huile résiduelle par injection de solutions contenant des tensioactifs. Les mécanismes physiques pris en compte dans les modélisations font intervenir la physico-chimie du système roche-fluide et une courbe globale donnant la saturation résiduelle en huile en fonction du nombre capillaire (courbe de désaturation capillaire). Cette donnée est majeure dans le calcul de l'efficacité de récupération d'huile par injection de solutions de tensioactifs. En effet la mobilisation de l'huile résiduelle laissée en place après injection d'eau n'est possible qu'en augmentant considérablement le nombre capillaire. La prédiction de l'efficacité d'un procédé chimique de récupération passe par la compréhension, à l'échelle du pore, du processus de mobilisation des ganglions d'huile suivant la structure poreuse et le nombre capillaire. L'objet de cette thèse est de caractériser la récupération d'huile tertiaire en fonction du nombre capillaire dans diverses roches mouillables à l'eau. Ces courbes permettront de quantifier l'effet de la microstructure, les hétérogénéités du milieu poreux et diverses propriétés pétrophysiques sur la récupération d'huile. Cette thèse permettra aussi de caractériser les différents mécanismes d'action de tensioactifs sur la mobilisation d'huile résiduelle dans le milieu poreux. L'expérimentation par tomographie RX est utilisée. La tomographie RX permettra de caractériser les courbes de désaturation capillaire à l'échelle de Darcy et visualiser localement le déplacement d'huile résiduelle à travers les milieux poreux. Des essais d'écoulement diphasique sous micro-CT permettront d'observer in-situ et d'étudier les interfaces eau/huile et leurs évolutions en 3D au sein du milieu poreux en fonction du nombre capillaire. / Oil recovery by surfactant injection is related to oil-water interfacial tension and rock properties through the capillary number. In the modeling of oil recovery by surfactant injection, fluid flow physical mechanisms are represented through the capillary desaturation curve (CDC). This curve is central in the evaluation of oil recovery efficiency. In order to mobilize residual oil trapped after waterflooding by capillary forces, chemical EOR rely on increasing capillary number to extremely high values. The mechanisms governing oil release can be described at the pore scale where the balance of capillary and viscous forces is achieved. This description will help to predict the efficiency of surfactant based EOR processes by taking into account the porous geometry and topology, the physico-chemical properties of the fluids and the different phase interaction. The objective of this work is to characterize capillary desaturation curves for various strongly water-wet sandstones. These curves will be used to study the relationship between tertiary oil recovery and the pore structure, porous media heterogeneity and petrophysicals properties. The other aim of this work is to map the different mechanisms of oil recovery by surfactant injection. Experiments under X-Ray tomography are proposed. X-Ray tomography will be applied to characterize capillary desaturation curve at Darcy scale and to visualise the two phase flow saturation after injection. Pore scale experiments based on X-Ray micro-tomography imaging are performed to describe the different mechanisms of oil mobilization.

Ce travail a porté sur l'application de la microtomographie à rayons X à la caractérisation de la géométrie de fractures et milieux poreux soumis à des réactions de dissolution, afin de mieux comprendre les relations existant entre les réactions chimiques, les modifications de géométrie, et l'évolution des écoulements et du transport réactif. Les observations ont permis de mettre en évidence la nature instable des phénomènes de dissolution, qui peuvent conduire à une augmentation rapide de la perméabilité. L'implication de la minéralogie et de la texture des roches naturelles paraît également déterminante dans l'évolution de la rugosité des surfaces, du transport des éléments à l'interface eau-minéral, et des motifs de dissolution. Le modèle de transport réactif HYTEC a ensuite été utilisé afin de mettre en évidence le rôle de la géométrie, de la relation porosité - diffusivité, et de la valeur des nombres de Damköhler et de Péclet dans l'évolution du transport des espèces en solution.

Afin d'étudier le rôle de l'humidité dans la déformation du sel gemme, des essais utilisant la microcellule développée au Centre de Géologie de l'Ingénieur ont été conduits sur des échantillons de sel laiteux et de sel phénoblastique du Bassin Bressan ainsi que sur des échantillons de sel de dôme provenant de la mine d'Asse (Allemagne). Le principe de la méthode utilisée est de suivre visuellement, lors des essais effectués, l'évolution texturale des échantillons sous des conditions de contrainte, température et humidité imposées. Une augmentation même faible du degré d'humidité relative se traduit toujours par une phase de fluage transitoire rapide suivie d'une phase de fluage stationnaire à vitesse de déformation constante. La vitesse de déformation lors de cette phase de fluage stationnaire est généralement supérieure à celle de la phase de fluage stationnaire qui précède l'élévation du degré d'humidité. La texture et la composition de l'échantillon ainsi que sa réaction vis-à-vis du couple température-humidité imposé déterminent la prédominance d'un des micromécanismes de déformation possibles. D'un point de vue microstructural, le contenu des cristaux de halite en inclusions fluides et/ou solides exerce une grande influence sur leur ductilité. On montre que le sel laiteux riche en inclusions fluides est plus ductile que le sel de dôme dont les cristaux contiennent des inclusions solides. Au niveau des joints de grains, l'humidité s'adsorbe ou se condense permettant le déclenchement de certains mécanismes de déformation spécifiques. Selon la morphologie et l'état de surface de ces joints on assiste à des mécanismes de lubrification et de glissement (cas du sel de dôme) ou des mécanismes de cisaillement et de microstylolitisation (cas du sel laiteux). Par ailleurs, les films de saumures qui se reconstituent au niveau des joints favorisent la thermomigration des inclusions fluides d'un cristal de halite à un autre à travers ces joints saturés. Dans le sel phénoblastique 1'anhydrite contenue dans la matrice carbonato-sulfatée microporeuse s'hydrate en gypse sous l'action du couple température-humidité. Cette transformation s'accompagne d'un gonflement important qui se traduit par la fissuration prononcée des plages non halitiques ce qui peut entraîner une dégradation notable du matériau. En prenant en compte les trois paramètres essentiels qui conditionnent le fluage du sel gemme: la texture, la température (avec et sans gradient) et l'humidité, l'étude expérimentale tend à reconstituer de manière réaliste l'état des conditions susceptibles de se réunir à proximité d'une galerie de stockage de déchets radioactifs. Par conséquent elle permet d'avancer des arguments nécessaires pour aider à la prévision du comportement d'un massif salifère à long terme.

L'influence des effets de site sur la réponse sismique des barrages-voûtes a été étudiée récemment par différents auteurs. Dans les travaux effectués jusqu'à présent l'excitation sismique est utilisée de différentes façons, mais toujours en se limitant à quelques angles d'incidence d'ondes sismiques, définis à l'avance. Les réponses des barrages-voûtes dans ces conditions ne représentent que des cas particuliers sur lesquels sont basées les conclusions. La méthode des surfaces d'énergie présentée dans ce travail, représente une méthode générale (et en même temps économique) qui se sert du balayage tridimensionnel fin sur tous les angles d'incidence dans l'espace, pour différents types d'ondes sismiques et leurs différentes célérités de propagation. L'excitation sismique est représentée par les ondes P, SV, SH et ondes de Rayleigh. Leurs expressions ont été développées en séries limitées et remplacées par différentes combinaisons de fonctions de forme dont les coefficients sont les fonctions d'angles d'incidence et de direction. Cette méthode nous permet de définir les intervalles d'angles d'incidence les plus dangereux et ceux les moins dangereux pour un barrage ainsi que l'influence et l'importance des effets de site sur la réponse sismique des barrages-voûtes. Dans les conditions générales d'angles d'incidence, plusieurs paramètres ont été variés et analysés à l'aide de 12 différents modèles qui représentent des systèmes complexes du barrage-voûte avec son réservoir et la masse rocheuse qui les entoure. Ces paramètres sont : les angles d'incidence , la position du noeud dans le barrage, les célérités de propagation d'ondes du champ proche dans la roche, les célérités de propagation d'ondes du champ lointain, les effets de site, les types d'ondes, la compressibilité de l'eau, la présence de la retenue pleine, le changement de la taille et de la forme du réservoir, la direction de la vallée, des petites modifications de la forme de la frontière du modèle numérique, l'utilisation des éléments fluides pour modéliser le réservoir. La méthode d'énergie maximale traite le problème extrémal de l'énergie dans le barrage sur l'espace des fonctions de forme. L'excitation sismique est représentée par un signal plus complexe que celui obtenu par des ondes P, S V, SH et ondes de Rayleigh. Dans les analyses sismiques sur les effets maxima on détermine également des mouvements "éventuels" de la roche qui donnent l'énergie maximale dans le barrage. Les conclusions donnent un apport théorique basé sur les analyses paramétriques effectuées dans ces conditions générales d'angles d'incidence, ainsi qu'un apport pratique qui montre le moyen de détermination de la position du barrage la plus favorable par rapport aux intervalles d'angles prédominants de séismes (en remplissant en même temps les autres conditions nécessaires: statiques, géologiques etc.). Dans le domaine théorique cette thèse élargit et classifie nos connaissances sur les phénomènes des effets de site, en premier lieu l'influence sur le barrage (structure étendue) du massif rocheux, dont les propriétés mécaniques sont exprimées par la célérité de propagation d'ondes sismiques et de l'interaction barrage-fluide du lac de retenue. La méthode présentée dans ce travail ainsi que le logiciel entièrement créé pour le même travail peuvent être applicables, avec peu de modifications, sur d'autres structures étendues dont le comportement sismique est décrit dans le domaine linéaire de façon suffisamment correcte.

Les déformations élastoplastiques (instantanées) et les mécanismes d'endommagement du sel gemme sont abordés par une étude expérimentale et théorique. Une étude bibliographique présente les mécanismes de déformation du sel gemme, et permet de dégager les différents facteurs physiques régissant le caractère ductile ou fragile du comportement du sel gemme. On situe ainsi le sel gemme par rapport aux autres géomatériaux du point de vue des propriétés mécaniques. On présente également différentes méthodes d'identification de l'endommagement des géomatériaux. Des essais de compression et d'extension axisymétriques sont présentés. La déformation volumique globale est mesurée au cours du chargement déviatorique. Cette mesure est un indicateur d'endommagement et permet en particulier d'identifier les conditions d'initiation de l'endommagement. Le chargement hydrostatique et les phases déviatoriques renseignent sur la rhéologie du matériau et sur la détérioration du sel gemme consécutive au chargement déviatorique. Les variations des paramètres élastiques sont décrites mais, étant de faible amplitude, elles ne sont pas prises en compte dans l'écriture de la loi de comportement. Une analyse des contraintes internes par le truchement d'un modèle micro-macro autocohérent permet d'éclairer les causes de la précocité de l'endommagement en extension. L'usage d'un tel modèle permet aussi de discuter la nature de l'écrouissage. A l'échelle macroscopique, un modèle de comportement élastoplastique est proposé. Il suit un écrouissage isotrope, l'écoulement est non-associé, et sa détermination s'appuie sur l'analyse de l'évolution de l'angle de dilatance en fonction de la contrainte moyenne et du paramètre d'écrouissage (la distorsion plastique). Une généralisation à l'espace tridimensionnel des contraintes (en faisant intervenir trois invariants du tenseur des contraintes) est développé et permet de présenter un résultat de calcul de creusement de galerie axisymétrique.

Cette thèse examine la réponse de différents types de lois de comportement en chargement et en déchargement et modélise numériquement les problèmes du creusement dans un massif du sol. Deux méthodes de modélisation d'excavation par éléments finis, l'une proposée par Ghaboussi et Pecknold (1984) et Brown et Booker (1985) et l'autre par Borja et al. (1989) ont été programmées dans le code CESAR-LCPC. Ces deux méthodes ont été ensuite comparées pour la simulation par étapes de l'excavation de massifs de sol élastique et élastoplastique. Nous avons également développé la notion de chargement-déchargement-rechargement multiaxial basée sur le signe du travail incrémental. L'application de ce principe (qui englobe le cas unidirectionnel) à une excavation pour différents types de lois de comportement, nous permet de connaître l'étendue des zones en chargement et/ou en déchargement autour de l'ouvrage. Un modèle élastique non-linéaire (Wong et Duncan), un modèle élastoplastique parfait à critère de Mohr-Coulomb et un modèle élastoplastique avec écrouissage (Nova) ont été étudiés plus en détail dans cette étude. Ces modèles ont été d'abord appliqués à la simulation d'essais triaxiaux, puis à un problème d'excavation. Ce travail s'achève par la modélisation d'un ouvrage réel : une excavation réalisée devant un rideau de palplanche butonné dans un milieu sableux (Université de Karlsruhe). Dans cette partie, nous avons comparé avec les observations sur site, les résultats des modélisations en éléments finis avec le modèle élastoplastique à critère de Mohr-Coulomb et avec le modèle de Nova.

La caractérisation géochimique et isotopique des eaux souterraines de la " zone minéralisée de l'Entre-Deux-Mers " indique une origine commune de la minéralisation, directement liée à la minéralogie des formations captées par les forages.La géochimie montre que les interactions eau-roche sont majoritairement influencées par la présence d'évaporites, mais que d'autres interactions mettant en jeu des carbonates, des silicates et des argiles existent. Un modèle géochimique d'acquisition de la minéralisation reconstitue parfaitement la chimie des eaux souterraines à l'échelle de la zone d'étude. Ce modèle, construit en se basant sur la géochimie des eaux et sur la minéralogie des formations tertiaires du nord du Bassin aquitain, met à l'équilibre des eaux avec des formations carbonatées et évaporitiques. Afin de mieux comprendre la distribution latérale et verticale des formations tertiaires et leur minéralogie, une approche paléogéographique et sédimentologique a permis de localiser les différents horizons riches en sulfates et/ou en fluorures, mais aussi de comprendre leur origine de dépôt. En se basant sur l'hydrogéologie, la paléogéographie, la minéralogie et la géochimie, des hypothèses de répartition de la minéralisation à l'échelle du forage ont pu être testées. Les résultats de la modélisation couplée hydrodynamique-transport reconstituent la chimie des eaux prélevées par les forages de la " zone minéralisée de l'Entre-Deux-Mers ". Au vu de ces résultats, un modèle avec obturation des horizons riches en sulfates et en fluorures a été testé et les résultats obtenus ouvrent des perspectives pour des futures recherches. Ce travail a donc permis de comprendre l'origine de la minéralisation des eaux de " la zone minéralisée de l'Entre-Deux-Mers ", mais aussi de proposer des améliorations et des perspectives pour une meilleure gestion d'une des principales ressources en eau potable de la Gironde.

L'énergie solaire est connue pour sa nature intermittente par rapport aux ressources d’énergie fossile. Cette observation souligne la nécessité d'utilisation d’un système de stockage d'énergie thermique. Le système de stockage thermocline est considéré comme un système de stockage rentable. La présente thèse vise à étudier le potentiel des roches basaltiques et siliceuses comme des candidates matériaux de stockage pour les centrales solaires concentrées. Les études expérimentales des propriétés thermo-physiques et thermomécaniques de ces roches à des températures allant jusqu'à 1000°C montrent que ces roches offrent de bonnes propriétés thermiques par rapport aux matériaux classiques de stockage. L'analyse du système de stockage thermocline sur un lit de roches à air direct est réalisée par une approche numérique. En outre, cette recherche vise également à évaluer l’impact environnementale de ce type de système de stockage en effectuant une analyse comparative de son cycle de vie. Enfin, une étude complémentaire réalisée dans le but de produire une carte d'indice de pertinence a permis d’identifier les zones les plus appropriées pour la construction des centrales solaires en Egypte. L'originalité de cette approche alternative pour le stockage d'énergie thermique est qu’elle combine la performance et la disponibilité des matériaux de stockage tout en réduisant leurs impacts environnementaux et financiers. / Compare to fossil fuel energy resources, solar energy is known for its intermittent nature. This observation highlights the need for the use of a thermal energy storage system. The thermocline storage system is considered as a cost-effective storage system. This thesis aims to study the potential of basalt and silex rocks as candidate storage materials for concentrated solar power plants. Experimental studies of the thermo-physical and thermo-mechanical properties of these rocks at temperatures up to 1000°C show that these rocks offer good thermal properties compared with conventional storage materials. The analysis of the thermocline storage system of air rock-packed bed is carried out using a numerical approach. This research also aims to assess the environmental impact of this type of storage system by conducting a comparative analysis of its life cycle. Finally, a complementary study carried out with the aim of producing a relevance index map made it possible to identify the most suitable areas for the construction of solar power plants in Egypt. The originality of this alternative approach for thermal energy storage is that it combines the performance and availability of storage materials while reducing their environmental and financial impacts.

Les structures en enrochements sont parmi les ouvrages les plus usuels de génie civil (barrages, murs de soutènement,. . . ). Des tassements importants peuvent apparaître tout au long de leur durée de vie et sont principalement dus à la rupture des blocs rocheux. Cette thèse propose un modèle numérique permettant de simuler le comportement de matériaux granulaires présentant des ruptures de grains. Afin de prendre en compte la nature discontinue de ces milieux, la méthode des éléments discrets est utilisée. La modélisation adoptée est de type "Non-Smooth Contact Dynamics", où les grains et particules sont supposés rigides. Afin de générer des blocs ayant des formes complexes, un modèle de grain 3D est proposé. Ce modèle de grain est ensuite discrétisé en sous-éléments de forme tétraédrique liés par des liaisons cohésives afin de pouvoir représenter la rupture. Un critère de rupture de Mohr-Coulomb est utilisé. Le modèle est implémenté sur la plateforme logicielle LMGC90. Le modèle est d’abord éprouvé lors de simulations d’écrasement de blocs cassables entre deux plaques. Plusieurs paramètres contrôlant la résistance du grain sont étudiés : cohésion intergranulaire, taille, discrétisation, forme et orientation du grain. L’effet d’échelle observé sur ce type de matériau est vérifié. Le modèle est ensuite testé lors de simulations numériques de compression œdométrique d’enrochements. L’effet des paramètres du modèle et de l’assemblage du milieu granulaire est également étudié. Les simulations œdométriques sont confrontées à des résultats expérimentaux et présentent une bonne concordance. Enfin, des expérimentations numériques sont menées afin d’étudier les énergies mises en jeu dans ces essais. L’énergie de création de surface est estimée pour ce type de matériau. Les résultats sont proches des données de la littérature. / Rockfill structures are very popular among civil engineering structures (dams, retaining walls, . . . ). Important settlements can take place during the lifetime of these structures, settlements mainly caused by the breakage of rockfill grains. This thesis proposes a numerical model that allows the simulation of the behavior of granular materials exhibiting grain breakage. To take into account the discrete nature of these media, the discrete element method is chosen. The adopted strategy is the Non-Smooth Contact Dynamics method, where grains are considered to be rigid. To generate blocks having complex shapes, a 3D grain model is suggested. This grain model is then discretized into tetrahedral subgrains, joined together using cohesive bonds so that breakage can be simulated. A Mohr-Coulomb failure criterion is used for the cohesive bonds. The model is implemented into the LMGC90 software platform. At first, the model is tested in single grain crushing simulations between two plates. Multiple parameters controling the strength of the grain are studied : the intra-granular cohesion, the size, the discretization and the orientation of the grain. The scale effect that characterizes this type of material is verified. Then the model is tested in numerical simulations of œdometric compression of rockfill. The influence of the parameters of the model and of those of the granular medium are studied. The results of œdometric simulations are compared to experimental results, and present a good agreement. Lastly, numerical experimentations are conducted in order to study the energies that are brought into play in the simulations. The surface creation energy is estimated for this type of material. Results are close to the data provided in the literature.

Dans cette étude, nous cherchons à dresser le bilan de l’altération à l’échelle de l’île de Mayotte (en surface et en souterrain), à contraindre l’impact des facteurs de contrôle dominant l’altération, puis à caractériser les processus impliqués dans les fractionnements des isotopes du silicium à l’échelle du profil d’altération.Deux campagnes d’échantillonnage incluant des eaux de rivières et des eaux souterraines ont été réalisées en saison humide et sèche, complétées par le prélèvement mensuel de 5 rivières. La prise en compte des crues (3 % de l’année) dans le calcul des taux d’altération moyens annuels en rivière a entraîné une augmentation de ≈32 % du bilan annuel. Le taux d’altération chimique global de l’île de Mayotte s’élève à 94 t/km²/an (81 t/km²/an en surface & 131 t/km²/an en souterrain). Nos résultats mettent en avant le rôle prépondérant des écoulements souterrains dans le transport de matériel dissous directement à l’océan. La contribution du domaine souterrain au bilan de l’altération diminue avec l’âge des formations, soulignant l’implication de l’âge des roches du bassin versant comme paramètre clé dans le contrôle des taux d’altération.Les isotopes du Si ont été analysés sur des roches totales le long de profils d’altération associés à différentes conditions d'altération du régolite : météorique (basse température) ou hydrothermale (haute température). À l’échelle du profil, les deux types de régimes ont montré un appauvrissement en 30Si en fonction du degré d’altération. À l’échelle du minéral, le fractionnement des isotopes du Si s’est révélé plus négatif pendant la précipitation de kaolinite secondaire à haute température qu’à basse température. / In this study, we seek to establish the weathering budget at the scale of the island of Mayotte (rivers and groundwater), to constrain the impact of dominant control factors to the weathering both locally and at a global scale, and then to characterize the processes involved in the fractionation of silicon isotopes across the weathering profile.Two field campaigns, in order to sample river water and groundwater, were carried out during wet and dry season, supplemented by the monthly monitoring of 5 rivers. By taking account of river floods (3% of the year) in the calculation of average annual weathering rates in the river has increased the annual weathering budget by ≈32%. The overall rate of Mayotte Island chemical weathering is 94 t/km²/yr (81 t/km²/yr from surface & 131 t/km²/yr from underground). Our results highlight the important role of groundwater flow to the dissolved material export directly to the ocean. The contribution of groundwater to the weathering budget decreases with the age of the geological formations, highlighting the involvement of the age of the rocks of the watershed as a key parameter in the weathering rates control.Si isotopes were analyzed for whole rock along two weathering profiles associated with different alteration conditions of the regolith: meteoric (low temperature) or hydrothermal (high temperature). At the weathering profile scale, the two types of alteration regimes showed 30Si depletion as a function of the degree of weathering. At the mineral scale, Si isotope fractionation was more negative during the secondary kaolinite precipitation at high temperatures than at low temperatures.