Le renforcement des fondations par soil mixing consiste à créer dans le sol des colonnes à module mixte obtenues en mélangeant le sol sur place avec un liant hydraulique. Ces colonnes, qui possèdent des caractéristiques mécaniques entre celles du béton et celles du sol, permettent d’améliorer la portance du sol d’assise et réduire les tassements excessifs des structures surjacentes. Ce travail de recherche consiste en une étude du comportement des fondations superficielles renforcées par des colonnes de soil mixing. Il vise à mettre en exergue les améliorations que cette technique apporte en termes d’augmentation de la capacité portante et de réduction des tassements. Dans une première étape, un matériau de soil mixing est formulé à partir d’un sable d’Hostun et d’un ciment CEM III. Une étude préalable a été réalisée pour déterminer les caractéristiques mécaniques de ce matériau de soil mixing, du sable d’Hostun, ainsi que les propriétés d’interface sable d’Hostun - soil mixing. Dans une seconde étape, une cuve expérimentale a été conçue pour permettre la réalisation de la procédure expérimentale du renforcement du sable par des colonnes de soil mixing. Plusieurs configurations du modèle physique ont été testées: Il s’agit d’essais de chargement sur fondations superficielles reposant sur un massif renforcé ou non. La taille de la semelle, le nombre des colonnes de renforcement, leurs âges, la densité de sol ainsi que le type de renforcement (renforcement type fondation mixte ou inclusion rigide) sont des paramètres dont l’effet a été étudié. Les essais de caractérisation, la conception du modèle expérimentale ainsi que les essais de chargement réalisés ont permis de mieux comprendre le comportement du sol renforcé. La connaissance des paramètres géomecaniques du sol et du matériau de soil mixing a permis de comprendre l’évolution de la capacité portante des semelles en présence des colonnes de renforcement. La présence d’un horizon porteur dans lequel les colonnes sont encastrées fait augmenter significativement la capacité portante de la semelle. Dans le cas du renforcement type inclusions rigides, l’augmentation de l’épaisseur du matelas granulaire réduit la capacité portante de la semelle renforcée. Cependant, elle fait augmenter la contribution des colonnes de renforcement dans la reprise des efforts appliqués sur la semelle. Une confrontation a été réalisée entre les résultats expérimentaux des essais de chargement et les résultats provenant d’un modèle éléments finis mis en œuvre avec le logiciel Abaqus. Les paramètres utilisés proviennent de nos essais de caractérisation réalisés dans le cadre de cette étude. Le modèle a présenté des résultats satisfaisant et a réussi de reproduire le comportement physique des semelles renforcées. Le modèle n’a non seulement réussi à reproduire les essais de chargement à l’échelle d’une colonne, mais aussi à l’échelle des semelles renforcées par une seule colonne et par quatre colonnes. / The reinforcement of shallow foundations by soil mixing consists on mixing the soil beneath with a hydraulic binder in order to obtain columns with mixed modulus. These columns, which have mechanical properties bounded between those of concrete and those of soil, can improve the soil bearing capacity and reduce excessive settlements as well. This research work is a phenomenological and quantitative study of the behavior of shallow foundations reinforced by soil mixing columns. It leads to analyze and highlight improvements that this technology brings in terms of bearing capacity increase and settlement reduction. In a first step, a soil mixing material is formulated from Hostun sand and cement CEM III. A preliminary study is conducted to determine the mechanical characteristics of the soil mixing material, Hostun sand and sand – soil mixing interface. In a second step, an experimental tank is designed to allow the installation of 1 m3 of hostun sand under different density. The experimental procedure allows also the installation of soil mixing columns that reinforce the soil. Finally a device is installed in order to carry on loading tests. Different configurations of the physical model were tested. Loading tests are carried on soil mixing columns and on unreinforced and reinforced shallow foundations. Size of the foundation, number of reinforcing columns, their age, soil density and type of reinforcement (rigid inclusion or mixed foundation) are parameters we change during this study. The characterization tests, the design of the experimental model, and the loading tests helped us to better understand the behavior of the reinforced soil. The knowledge of geomechanical parameters of soil and soil mixing material allowed us to understand the evolution of the bearing capacity of the reinforced foundations. The presence of a carrier horizon, in which the tip columns is embedded, leads to significantly increase the bearing capacity of the foundation. In the rigid inclusions reinforcement case, the increasing of granular mattress thickness reduces the bearing capacity of the reinforced foundation. However, it increases the part of the loading transmitted to soil mixing columns. A comparison was made between the experimental results of loading tests and the results of a finite element model implemented with Abaqus software. Model parameters are taken from our material characterization study. It is proven that the model developed is able to reproduce not only loading tests at column scale, but alsot at reinforced foundations scale.

Dans le cadre de la mise à niveau sismique des bâtiments existants, la technique de renforcement par placage et enrobage de polymères renforcés de fibres (FRP) offre une réponse pertinente. L’objectif de cette thèse est de proposer une modélisation fiable pour la détermination de la réponse quasi-statique et dynamique d’une grande variété d’éléments de structure et d’identifier les apports possibles d’une telle modélisation à la conception. Les stratégies de modélisation s’appuient sur l’utilisation d’éléments finis massifs ou basés sur une cinématique simplifiée (coques multicouches et poutres multifibres), associés à des modèles locaux d’endommagement et de plasticité pour les matériaux en présence (béton, armatures et FRP). L’hypothèse d’adhérence parfaite, d’une part, entre les armatures et le béton, et d’autre part, entre les bandes de FRP et le béton, permet de calculer efficacement les cas des poteaux courants, des poteaux courts, des voiles longs et courts, dans les régimes quasi-statique (pushover) et dynamique. Le renforcement par placage et enrobage FRP pour des structures en béton armé, est reproduit au sein de la modélisation par un apport de matière (éléments de type barre avec les caractéristiques propres des FRP) et des modifications des paramètres de la loi de comportement du béton, justifiées par l’expérimentation et la littérature. La pertinence de l’approche est finalement démontrée en confrontant la réponse temporelle de structures à l’échelle 1, de type ossature ou de type contreventé par des voiles, aux résultats expérimentaux issus de benchmarks internationaux. / In the context of the building’s protection against seismic risk, the strengthening technique by FRP (Fiber Reinforced Polymers) plating and wrapping provides a relevant solution. The objective of this thesis is to propose a reliable modeling for determining the quasi-static and dynamic responses of a wide variety of structural elements and to draw advantage in design. The modeling strategies make use of 3D elements or finite elements based on simplified kinematics (multilayer shells or multifiber beams), associated with local damage and plasticity laws for the constitutive materials (concrete, rebar and FRP). The perfect bond assumption between steel-concrete and FRP-concrete allows efficiently calculating the quasi-static and dynamic behaviors of short and slender columns, as well as short and slender walls. The mechanical contribution of FRP plating and wrapping is reproduced in the modeling by adding material (bar type elements with FRP characteristics) and by changing the parameters of the concrete model on the basis of the experience and analytical formula issued from the literature. The relevance of the approach is finally highlighted by comparing the time-history response of real braced frame or wall structures with the experimental results.

Cette contribution, en s’appuyant sur un travail expérimental et de modélisation numérique à l’aide de la méthode des éléments distincts, vise l’étude de structures en maçonnerie et de structures en pisé. Pour la partie maçonnerie, notre étude traite tout d’abord de modèles réduits testés sur table inclinable, puis dans un second temps nous avons testé des murs en briques à une échelle représentative. Dans la deuxième partie, nous avons mené des essais sur des murs en pisé sous un chargement en cisaillement, avec ou sans précontrainte axiale de confinement qui correspond à une descente de charge. Une modélisation numérique par la méthode des éléments discrets a été réalisée pour ces différents cas d’étude. La comparaison entre les résultats expérimentaux et numériques, nous a permis, d’évaluer les pertinences et limites de la modélisation via la méthode des éléments discrets (DEM). / This contribution, based on experimental work and numerical modeling using the distinct elements method, aims to study masonry structures and rammed earth structures. For the masonry part, our study first deals with reduced models tested on tilting table, then in a second time we tested brick walls on a representative scale. In the second part, we conducted tests on rammed earth walls under a shear loading, with or without axial prestressing of the containment which corresponds to a descent of load. Numerical modeling using the discrete element method has been carried out for these different case studies. The comparison between the experimental and numerical results, allowed us to evaluate the relevance and limits of modeling via the discrete element method (DEM).

L'étude du système aquifère du Chari Baguirmi vise à améliorer sa connaissance en vue d'une gestion durable etprudente de la ressource en eaux souterraines. Dans la zone d'étude qui couvre 70000 Km2, des campagnes demesures et d'échantillonnages des eaux et de sédiments ont été réalisés. Les eaux échantillonnées ont fait l'objetd'analyses chimiques et isotopiques (O18 et 2H) afin de comprendre le comportement hydrogéochimique dusystème aquifère. La synthèse des données géologiques existantes et la mise en évidence de l'hétérogénéité litho-stratigraphiquepar l'étude granulométrique ont confirmé la complexité hydrogéologique de la zone d'étude. La minéralogie desargiles a montré que le dépôt des différents sédiments se fait de la périphérie vers le centre de la dépression. Cemode de sédimentation permet d'avancer l'hypothèse de l'existence d'une dépression vers laquelle les sédimentsseraient entrainés par l'agent de transport. Cette hypothèse de “ dépression structurale ” réfute l'hypothèseadmise jusqu'à présent d'une origine hydrogéologique (forte évaporation) de cette dépression piézométrique.L'étude basée sur la piézométrie et les données chimiques et isotopiques (O18 et 2H) des eaux a permis decomprendre les processus hydrogéochimiques qui gouvernent la minéralisation des eaux. En plus des isotopes(O18 et 2H), l'utilisation des ratios Br/Cl, Sr/Ca et l'indice d'échanges de bases a permis de mettre en évidenced'une part, l'origine de la salinité des eaux de la dépression piézométrique et d'autre part, les phénomènesd'échanges de bases liés au long temps de séjour des eaux au contact des roches encaissantes et les zones derecharge potentielles de la nappe. / The study of the Chari Baguirmi aquifer system aims to improve the knowledge about this system forsustainable groundwater resource management. In the study area, which covers 70,000 km2, measurementcampaigns and sampling of water and sediments were undertaken. The waters have been sampled for chemicaland isotopic analyzes (18O and 2H) to understand the hydrogeochemical behavior of the aquifer system.The synthesis of existing geological data and highlighting the litho-stratigraphic heterogeneity through thegranulometric study confirmed the hydrogeological complexity of the study area. Clay mineralogy showed thatthe deposition of different sediments occurred from the periphery to the center of the depression. This mode ofsedimentation can hypothesize the existence of a depression into which sediments are trained by the transportagent. This assumption of "structural depression" refutes the hypothesis of a hitherto accepted hydrogeologicalorigin (high evaporation) of the piezometric depression. The study based on the piezometry and chemical and isotopic data (18O and 2H) of water enabled to understand the hydrogeochemical processes that govern the mineralization of the water. In addition to isotopes (2H and 18O), the use of ratios Br/Cl, Sr/Ca and base exchanges index has highlighted on the one hand, the origin of the salinity of waters in the piezometric depression and secondly, the phenomena of base exchanges related to longresidence time of water in contact with rocks and areas of potential groundwater recharge. Numerical modeling of the aquifer was performed in steady state and the results confirm the prevalence of the phenomenon of evaporation in the western p

Les submersions marines d’origine météorologique sont des catastrophes naturelles majeures, responsables chaque année de milliers de morts et de milliards d’euros de dégâts. La partie centrale du Golfe de Gascogne est un territoire particulièrement vulnérable à cet aléa, comme nous l’a rappelé la forte submersion engendrée par la tempête Xynthia en 2010. L’objectif de ce travail est d’améliorer la compréhension des surcotes et des submersions marines dans cette zone de France dans une approche pluridisciplinaire mêlant la géomorphologie, l’océanographie et l’analyse d’archives historiques. Afin de juger du caractère exceptionnel de Xynthia, une recherche de l’ensemble des submersions marines qui ont affecté la région d’étude depuis 500 ans a été menée. La modélisation numérique des surcotes des 5 tempêtes engendrant des submersions au 20ème siècle, révélées par ces recherches, montre que des conditions météo-marines variées ont induit des niveaux d’eau et des submersions comparables à ceux provoqués par Xynthia. Ce constat est en désaccord avec les estimations de périodes de retour de niveau marin extrêmes basées sur l’analyse statistique de mesures marégraphiques et met en avant l’apport de l’approche historique dans de telles problématiques. Devant la forte vulnérabilité des Pertuis Charentais aux submersions marines, la modélisation statique de la submersion marine, méthode simple mais néanmoins fréquemment utilisée pour estimer l’extension des zones inondées, a été évaluée. Cette méthode fournit de bonnes estimations de l’extension de l’inondation dans les zones de faibles altitudes caractérisées par une faible distance entre le trait de côte et la limite continentale de la zone inondable, mais mauvaises lorsque cette distance est grande. En effet, lorsque l’inondation se propage loin du trait de côte, la dynamique de l’écoulement ne peut plus être négligée sur ces grandes distances. Afin d’anticiper de futures submersions, deux configurations des digues ont été testées par modélisation numérique au travers de l’exemple de l’estuaire de la Charente. Les hauteurs d’eau et l’inondation de Xynthia sont simulées en augmentant la hauteur des digues de l’estuaire, empêchant toute inondation des zones basses adjacentes ; puis en abaissant les digues bordant l’estuaire au niveau des plus hautes marées astronomiques et en créant une seconde rangée de digues protégeant les zones habitées. Cette seconde configuration permet l’inondation des zones non habitées mais empêche l’inondation des zones à enjeux importants, comme la ville de Rochefort. Il est montré que la rehausse de l’ensemble des digues entraîne des niveaux d’eau supérieurs de 1.2 m à Rochefort par rapport à la simulation sans modification de digues, alors que l’abaissement de celles-ci et la protection des zones à forts enjeux ne modifient pas la hauteur d’eau dans l’estuaire. Ainsi, la rehausse des digues côtières n’est pas une solution systématique car la protection contre l’inondation de toutes les zones côtières peut augmenter la vulnérabilité des zones à forts enjeux. / Storm-induced coastal flooding are major natural disasters, responsible for thousands of deaths and billions of euros of damages each year. The central part of the Bay of Biscay is vulnerable to this hazard, as recently shown by the strong flooding induced by Xynthia in 2010. This study aims to improve the understanding of storm surges and coastal floods in this area of France, using several methods such as numerical or static modeling. To assess the uniqueness of Xynthia, historical researches of coastal floods affecting the study area for 500 years was conducted. Numerical modeling of the storm surges related to 5 storms of the 20th century revealed by these researches shows that various meteo-oceanic settings induced water levels and coastal floods comparable to those caused by Xynthia. This finding challenges return periods estimations of extreme sea levels based on statistical analysis of tide gauges measurements and highlights the contribution of the historical approach to such issues. Given the high vulnerability on coastal floods of Pertuis Charentais, static modeling, a simple but frequently used method to estimate the extension of flooded areas is evaluated. This analysis shows that this method provides good estimations of flood extents in low-lying areas characterized by a small distance between the shoreline and the continental limits of the lower area, but bad estimations when this distance is large. These poor performances when floods spread away from the coastline are explained by the dynamics of the flow, which can no longer be ignored. Two coastal defenses strategies are investigated in the Charente-river Estuary by numerical modeling. Water levels and coastal floods induced by Xynthia are simulated with increased dikes height, preventing flooding of adjacent low-lying areas, and then with dikes lowered to highest high spring tide height and with a second rank of dikes preventing flooding of important issues areas, such as the town of Rochefort. It is shown that raising dikes leads to higher levels of 1.2 m in Rochefort compared to the simulation without changing dikes, while protecting issues do not affect the water level in the estuary. Thus, it is demonstrated that the systematic raising of dikes is not a solution because it can increase the vulnerability of important issues areas.

Cette thèse s'attache à la compréhension des phénomènes de couplages sismo-électromagnétiques survenant dans des milieux poreux contenant des fluides, ainsi qu'au développement de méthodes d'imagerie géophysique basées sur ces conversions. Deux axes de travail sont investigués. Le premier concerne le traitement des données sismo-électromagnétiques et notamment l'extraction des réponses interfaciales, qui peuvent renseigner sur le sous-sol en profondeur. Nous comparons l'impact de plusieurs méthodes de filtrage sur les amplitudes et les formes d'onde de ces réponses, dont une nouvelle technique de filtrage en domaine des curvelets, qui s'appuie sur les similarités entre les données sismiques et cosismiques. D'autre part, nous étendons la théorie de la propagation couplée des ondes sismiques et électromagnétiques à des milieux partiellement saturés pour le cas d'un mélange air/eau. Nous adaptons un programme de simulation numérique des ondes sismo-électromagnétiques à ces conditions. / This thesis aims to improve our understanding of seismo-electromagnetic coupling phenomena occuring in fluid-containing porous media and to develop geophysical imaging techniques based on conversions.Two areas of work are explored. The first encompasses seismoelectromagnetic data processing, especially the extraction of interface responses, which may provide information about the subsurface at depth. We compare the influence of several filtering techniques on the amplitudes and waveforms of the interface responses, among which a new filtering technique in the curvelet domain, taking advantage of the similarity between seismic and coseismic data. Furthermore, we extend the theory for the coupled propagation of seismic and electromagnetic waves here to partially saturated media for the case of a water/air mixture. We adapt a seismoelectromagnetic wave modelling program to these conditions.

Cette thèse porte sur une nouvelle approche pour le contrôle d’écoulement aérodynamique. Cette nouvelleapproche est basée sur l’utilisation d’actionneurs plasma. La modélisation numérique peut être une outilpuissante entre les mains des scientifiques et des ingénieurs pour comprendre, optimiser et ainsi ouvrir lavoie à la commercialisation et l’application de cette technologie. Le couplage entre l’électromagnétisme, leplasma et l’écoulement, nécessite des modèles et des techniques numériques avancées. Le travail présentédans cette thèse, a pour principaux objectifs : le développement et la validation de méthodes numériques poursimuler efficacement le fonctionnement de certains des plus importants types d’actionneurs plasma. Nousnous sommes intéressés à trois types d’actionneurs plasma : les décharges micro-ondes, la décharge à barrièrediélectrique (DBD) et le jet synthétique plasma (JSP).En ce qui concerne les décharges microondes, les objectifs sont plus fondamentaux que pour les autrestypes d’actionneurs. Il s’est agit de mieux comprendre la création du plasma, son évolution et de calculerl’efficacité énergétique de dispositifs microondes par la simulation numérique. Un schéma couplé implicite(ADI) - FDTD avec un modèle de plasma fluide simplifié est présenté. Cette formulation conserve la simplicitéet la robustesse des systèmes de FDTD, tout en dépassant la barrière du critère de stabilité CFL. Elle conduità un temps de calcul réduit et la possibilité de réaliser des simulations tridimensionnelles de la formationdu plasma et de l’évolution d’un plasma dans un champ micro-ondes. Afin d’étudier l’énergie absorbée parle plasma et le transfert vers le gaz sous forme de chaleur ainsi que le changement consécutif de la densitédu gaz, un solveur Euler a été couplé avec le modèle EM-plasma en tenant compte des effets de gaz réel.Diverses validations et applications sont ensuite étudiés. Des simulations tridimensionnelles de formationdu plasma sont réalisée qui montrent la formation de structures dans une décharge micro-ondes librementlocalisée. Les effets de chauffage de gaz sur le développement d’un "streamer" et la durée d’un volumepré-ionisé avec des champs sous-critiques sont également calculés.En ce qui concerne les deux autres groupes d’actionneurs, les objectifs de cette thèse se concentrent sur lamodélisation de leur fonctionnement et sur la production d’écoulement qui en résulte. Le Jet Synthétique Plasma a été numériquement étudié par trois modèles couplés. Les résultats obtenus sont prometteurspour l’optimisation du JSP et une meilleure compréhension des mécanismes qui limitent ses performances.L’actionneur DBD a été modélisée en utilisant deux solveurs différents basés sur des modèles physiquessimilaires - celui développé à l’ONERA et l’autre à LAPLACE. Des études paramétriques ont montré queles modèles donnent une estimation assez précise de la force produite par le DBD par rapport à des mesuresexpérimentales. Des applications aérodynamiques de contrôle d’écoulement ont démontré les effets possiblesde ces actionneurs pour la transition laminaire - turbulente et l’amélioration de la portance. Ces travauxouvrent une perspective nouvelle dans la conception et l’optimisation de ces actionneurs. / As aerodynamic flow control still remains one of the top subjects of research in the aerospace scientific world, new ways to perform such a control are being constantly studied. Plasma actuators based on momentum or energy addition in the flow, have been proven capable of positively modifying the flow aerodynamic features. Nevertheless, the development and optimization of such actuators, require further understanding of the basic multi-scale physics involved. In this thesis, we are interested in the numerical modeling of plasma flow control actuators. Three types of plasma actuators are considered: Microwave Plasma Discharges (MPD), the Dielectric Barrier Discharge (DBD) and the Plasma Synthetic Jet (PSJ). Concerning MPDs, a novel implicit approach has been developed which with have enabled three-dimensional simulations in time domain in reduced CPU time. The microwave breakdown and evolution of the plasma due to the electromagnetic waves has been studied numerically, demonstrating the three-dimensional nature of such discharges. Coupling of the EM-plasma model with an Euler based solver accounting for real gas effects, have revealed the plasma modification due to the intense gas heating. For the PSJ actuator, the numerical solver consists of three coupled numerical models and the obtained results of its operation offer important information of its performance and its limits. The DBD actuator has been numerically studied using 2 different solvers (based on the same physical model). Both solvers were capable to give quite accurate estimations of the induced force due to the plasma and various parametric studies have been conducted. These studies offer new perspectives in the understanding and the optimization of plasma actuators for flow control purposes.

Utilisée depuis des milliers d’années sous ses manifestations naturelles par l’Homme, cette ressource fait l’objet d’une exploitation commerciale depuis seulement le XXe siècle, à destination du chauffage de bâtiments, de certains usages industriels ainsi que de la production d’électricité. La France compte parmi les pionniers concernant l’usage direct de la chaleur alors qu’aucune filière industrielle n’est véritablement effective pour la production d’électricité d’origine géothermique. Le projet sélectionné, intitulé FONGEOSEC, a pour objectif la conception et la réalisation d’un démonstrateur innovant préindustriel d’une centrale géothermique haute enthalpie exploité par cogénération d’électricité et de chaleur. Un travail de recherche et développement, conduit par un consortium composé de partenaires industriels et scientifiques, vise au lancement de la filière industrielle géothermique haute température en France. L’objectif général des travaux de thèse porte sur une meilleure compréhension globale des comportements hydrauliques, massiques et thermiques des formations profondes en réponse à une sollicitation anthropique de longue durée. Il s’agit d’identifier les paramètres clés régissant la réponse du complexe réservoir à la suite d’une exploitation géothermique. Un point particulier sera consacré à caractériser la part de chacun des modes de transport de chaleur en milieu poreux – conduction thermique, convection libre et forcée – dans l’établissement des performances thermiques de l’ouvrage considéré. Plusieurs dispositifs techniques d’exploitation seront proposés afin de réduire les incertitudes associées au système géothermique souterrain et garantir le succès du projet FONGEOSEC. L’impact des mécanismes thermo-convectifs au voisinage des forages d’exploitation géothermique de grande profondeur reste peu documenté, a fortiori dans le cas de dispositifs déviés adoptant une complétion particulière. L’outil retenu pour l’évaluation des performances du dispositif au contact de l’encaissant est la modélisation numérique distribuée. La variabilité des propriétés physiques de l’hydrosystème, de la conception et des modalités d’exploitation du dispositif sur le comportement hydraulique et thermique de l’exploitation est envisagée selon différentes approches développées à partir de modèles numériques 3D. / Used for thousands of years under its natural manifestations, this resource has been commercially exploited since the twentieth century, for the heating of buildings, certain industrial uses and the production of electricity. France is one of the pioneers in the direct use of heat, whereas no industrial cluster is truly effective for the production of geothermal electricity. The selected project, FONGEOSEC, aims to design and produce an innovative pre-industrial demonstrator of a high enthalpy geothermal power plant operated by cogeneration of electricity and heat. A research and development project, led by a consortium of industrial and scientific partners, aims to launch the high-temperature geothermal industrial sector in France. The general objective of this thesis is to improve the understanding of the hydraulic, mass and thermal behavior of deep porous formations in response to long-term anthropogenic stress. The aim is to identify the key parameters governing the response of the reservoir complex related to geothermal operation. A particular point will be devoted to characterize the part of each mode of transport of heat in porous medium – thermal conduction, free and forced convection – in the establishment of the thermal performances of the geothermal power plant. Several technical operating devices will be proposed to reduce the uncertainties associated with the underground geothermal system and guarantee the success of the FONGEOSEC project. The impact of thermo-convective mechanisms in the vicinity of deep geothermal borehole remains poorly documented, especially in the case of deviated wells with a complex inner geometry. The evaluation of the hydraulic and thermal performances of the device, based on 3D numerical modeling, is conducted according to different approaches.

Dans ce travail, nous étudions un schéma élément fini que nous appliquons à un modèle décrivant l'intrusion saline dans les aquifères côtiers confinés et libres. Le modèle est basé sur l'approche hydraulique qui consiste à moyenner verticalement le problème initial 3D, cette approximation repose sur une hypothèse d'écoulement quasi-hydrostatique qui, loin des épontes et des sources, est vérifiée. Pour modéliser les interfaces entre l'eau douce et l'eau salée (respectivement entre la zone saturée et la zone sèche), nous combinons l'approche 'interface nette' à l'approche avec 'interface diffuse' ; cette approche est déduite de la théorie de champ de phase, introduite par Allen-Cahn, pour décrire les phénomènes de transition entre deux zones. Compte tenu de ces approximations, le problème consiste en un système fortement couplé d'edps quasi-linéaires de type parabolique dans le cas des aquifères libres décrivant l'évolution des profondeurs des 2 surfaces libres et de type elliptique-prabolique dans le cas des aquifères confinés, les inconnues étant alors la profondeur de l'interface eau salée/eau douce et la charge hydraulique de l'eau douce. Dans la première partie de la thèse, nous donnons dans le cas d'un aquifère confiné, des résultats d'estimation d'erreur d'un schéma semi-implicite en temps combiné à une discrétisation en espace de type élément fini Pk Lagrange. Ce résultat utilise entre autre un résultat de régularité du gradient de la solution exacte dans l'espace Lr(ΩT), r > 2, ce qui permet de traiter la non-linéarité et d'établir l'estimation d'erreur sous des hypothèses de régularité raisonnables de la solution exacte. Dans la seconde partie de la thèse, nous généralisons l'étude précédente au cas de l'aquifère libre. La difficulté principale est liée à la complexité du système d'edps paraboliques mais à nouveau, grâce au résultat de régularité Lr(ΩT), r > 2 établi pour les gradients des surfaces libres, nous montrons que le schéma est d'ordre 1 en temps et k en espace pour des solutions suffisamment régulières. Nous concluons ce travail par des simulations numériques dans différents contextes (impact de la porosité et de la conductivité hydraulique sur l'évolution de l'interface, pompage et injection d'eau douce, effet des marées) validant ainsi le modèle et le schéma. Puis nous comparons les résultats à ceux obtenus avec un schéma volume fini construit à partir d'un maillage structuré. / In this work, we study a finite element scheme we apply to a model describing saltwater intrusion into coastal aquifers confined and free. The model is based on the hydraulic approach of vertically averaging the 3D original problem, this approximation is based on a quasi-hydrostatic flow hypothesis which, instead of the walls and springs, is checked. To model the interface between freshwater and salt water (respectively between the saturated zone and dry zone), we combine the approach net interface (approach with the diffuse interface) ; This approach is derived from the phase field theory introduced by Allen-Cahn, to describe the phenomena of transition between two zones. Given these approximations, the problem consists of a strongly couple to edps parabolic quasi-linear system in the case of unconfined aquifers describing the evolution of the depths of two free surfaces and elliptical-parabolic type in the case confined aquifer, the unknowns being then the depth of salt water / fresh water and the hydraulic load of fresh water. In the first part of the thesis, we give in the case of a confined aquifer, error estimation results of a semi-implicit scheme in a combined time discretization space finite element type Pk Lagrange. This result among other uses a regularity result of the gradient of the exact solution in the space Lr(ΩT), r > 2, which can handle the non-linearity and to establish the error estimate under assumptions reasonable regularity of the exact solution. In the second part of the thesis, we generalize the previous study to the case of the free aquifer. The main difficulty is related to the complexity of the system of parabolic edps but again, thanks to regularity result Lr(ΩT), r > 2 gradients established for the free surfaces, we show that the scheme is of order 1 time and space k for sufficiently regular solutions. We conclude this work by numerical simulations in different contexts (impact of porosity and hydraulic conductivity of the evolution of the interface, and pumping fresh water injection, tidal effects) thus validating the model and diagram. The we compare the results with those obtained using a finite volume scheme constructed from a structured mesh.

Dans un contexte de changement climatique, il est important d'étudier les processus de distribution et les tendances des gaz à l'état de traces dans l'atmosphère tels que l'ozone et la vapeur d'eau. L'observatoire atmosphérique du Maïdo, situé à 2160 m, héberge un ensemble d'instruments de mesures atmosphériques in-situ et par télédétection, dont le Lidar1200, pour la mesure de profils verticaux de vapeur d'eau. La validation, l'évaluation et une première exploitation des données du Lidar1200 font l'objet de cette thèse. Des comparaisons avec des radiosondages ont permis de valider les profils de vapeur d'eau du lidar depuis la basse troposphère jusqu'à 22 km d'altitude dès lors qu'un nombre suffisant de données est intégré. À courte échelle de temps, le Lidar1200 est capable de suivre temporellement des structures de l'ordre de la centaine de mètres d'épaisseur jusqu'à une dizaine de kilomètres d'altitude. L'analyse des données montre le potentiel du Lidar1200 de surveillance de la vapeur d'eau à long-terme dans la basse stratosphère avec un ou deux points de mesure par mois. La base de données 2013-2015 apporte des informations sur la variabilité saisonnière de la vapeur d'eau au-dessus de La Réunion. Les données de vapeur d'eau permettent de documenter divers processus atmosphériques dont les échanges stratosphère-troposphère. Des signatures d'intrusion stratosphérique ont été observées sur des observations lidar quasi simultanées et co-localisées d'ozone et de vapeur d'eau pendant la campagne de mesure MALICCA (MAïdo LIdar Calibration Campaign, La Réunion), le 4 avril 2013, suggérant un processus d'intrusion stratosphérique avec plusieurs origines. / In a context of climate change, it is important to study distribution processes and the trends of trace atmospheric gases such as ozone and water vapor. The Maïdo atmospheric observatory, located at 2160 m in height, hosts a set of atmospheric in-situ and remote sensing instruments, including the Lidar1200 providing vertical water vapor profiles. This thesis aims to provide a thorough work of validation, evaluation and analysis of the Lidar1200 data. Comparisons with radiosoundings allowed to validate the water vapor profiles of the lidar from the lower troposphere up to an altitude of 22 km if a sufficient number of data is integrated. At a smaller temporal scale, the Lidar1200 is able to follow temporally structures in the range of a hundred of meters thick up to a dozen of kilometers in altitude. The data analysis show the potential of the lidar in monitoring the water vapor in the lower stratosphere on a long term with one to two points of measurements by month. The database 2013-2015 brings information on the seasonal variability of the water vapor over Reunion Island. The water vapor data allow to study various atmospheric processes including stratosphere-troposphere exchanges. Signatures of stratospheric intrusions were observed on quasi-simultaneous and co-localized ozone and water vapor observations realized by two lidars of the Maïdo Observatory during the measurement campaign MALICCA-1 (MAïdo LIdar Calibration CAmpaign, Reunion Island), on April 4, 2013, suggesting stratospheric intrusions with multiple origins.