Le vieillissement croissant de la population mondiale augmente les risques d’incidence de maladies dégénératives comme la maladie d’Alzheimer (AD). L’AD représente la forme de démence la plus commune. C’est plus de 40 millions de personnes qui sont affectés à l’échelle mondiale, créant ainsi un lourd fardeau économique annuel de plusieurs centaines de milliards de dollars. L’AD est une maladie dégénérative du cerveau qui se traduit par une perte graduelle de la mémoire et des fonctions cognitives. Cette dernière se caractérise par trois symptômes pathophysiologiques intimement reliés qui sont : (1) un système cholinergique défectueux, (2) une accumulation excessive de plaques séniles toxiques composées de peptides β–amyloïdes et (3) une hyperphosphorylation de la protéine Tau qui résulte en des enchevêtrements neurofibrillaires. Il n’existe à ce jour aucune cure contre l’AD et les traitements actuellement retrouvés sur le marché n’ont qu’un effet transitoire. En revanche, il y a de plus en plus d’indices qui suggèrent que l’utilisation de certains facteurs de croissance comme l’IGF-2, le bFGF et les neurotrophines pourraient agir à titre de traitement. Parmi ces derniers, la protéine morphogénétique osseuse 9 (BMP-9) a montré un fort potentiel thérapeutique. Cette dernière est toutefois coûteuse à produire en plus d’être volumineuse, ce qui limite fortement son acheminement au cerveau en raison de la barrière hématoencéphalique. Dans cette optique, nous avons développé deux peptides de 23 acides aminés dérivés de l’épitope knuckle de la BMP-9, pBMP-9 et SpBMP-9, 300 fois moins chers à produire. Les travaux de ce projet de doctorat portent sur le développement d’un système de libération du SpBMP-9, le plus efficace des deux peptides, comme stratégie de traitement contre l’AD. Tout d’abord, il a été montré que le SpBMP-9 permettait d’induire la différenciation neuronale de cellules humaines SH-SY5Y de façon plus efficace que la protéine native par la présence d’une croissance plus importante des neurites et l’expression de marqueurs de différenciation neuronaux. De plus, le SpBMP-9 pouvait agir sur au moins deux symptômes de l’AD en orientant la différenciation vers le phénotype cholinergique et en inactivant la GSK3β, une Tau kinase. Ensuite, un système de libération sous forme de nanoparticules à base de chitosane et d’alginate a été développé, puis caractérisé. Ce système de libération a permis d’encapsuler efficacement le SpBMP-9. Le suivi et la modélisation des cinétiques de libération à l’aide d’un modèle mathématique prenant en considération la distribution de taille des particules a permis d’identifier que la libération était principalement diffusive, mais qu’il existait de fortes interactions entre le peptide et l’alginate. Les tests de viabilité ont montré que les nanoparticules n’étaient pas toxiques sur des cellules SH-SY5Y. Enfin, le SpBMP-9 libéré à partir des nanoparticules était toujours bioactif et permettait une différentiation neuronale. Finalement, une étude exploratoire a permis de mettre en évidence que le SpBMP-9 pouvait agir en synergie avec le bFGF et le NGF. La combinaison du SpBMP-9 avec le bFGF ou le NGF a montré la présence d’une différenciation neuronale accrue mise en évidence par la présence d’une croissance de neurites supérieure, par l’expression de plusieurs marqueurs de différenciation ainsi que par un niveau de calcium intracellulaire plus important. Pour conclure, ces travaux de recherche suggèrent que le SpBMP-9 possède un fort potentiel thérapeutique dans le contexte de l’AD. Une étude plus approfondie de ses mécanismes d’action est requise ainsi que la validation de son effet in vivo sur modèle animal.

Dans ce travail nous avons développé des modèles de Monte Carlo cinétique d’objets (OKMC) qui ont permis de prédire l'évolution nano-structurelle des amas de lacunes et des auto-interstitiels sous irradiation neutronique, à la température de fonctionnement des réacteurs de génération II dans les alliages Fe-C-MnNi (alliages modèles pour les aciers de cuve) et Fe-Cr-C (matériaux envisagés pour les réacteurs de génération IV). Un véritable acier de cuve venant du programme de surveillance de la centrale nucléaire suédoise de Ringhals a aussi été modélisé. Pour ce faire nous avons développé deux modèles OKMC fondés sur les données les plus actuelles concernant la mobilité et la stabilité des amas de défauts. Les effets des solutés d'intérêt ont été introduits dans nos modèles dans l’hypothèse simplifiée ‘‘d’alliage gris’’, c'est-à-dire que les solutés ne sont pas explicitement introduits dans le modèle, qui ne peut donc pas décrire leur redistribution, mais leur effet est introduit dans les paramètres liés à la mobilité des amas de défauts. A l’aide de cette approche nous avons modélisé diverses conditions de température et de débit de dose ainsi que des études de recuits isochrones d’alliages Fe-C-MnNi. L'origine du durcissement par irradiation neutronique à basse température a également été étudiée et les modèles ont fortement soutenu l'hypothèse selon laquelle les solutés ségrégent sur des boucles interstitielles immobiles, qui agissent donc comme des sites de nucléation hétérogène pour la formation d’amas enrichis en NiSiPCr et MnNi. A chaque fois nos modèles ont été validés par comparaison des résultats obtenus avec les observations expérimentales disponibles dans la littérature. / We developed object kinetic Monte Carlo (OKMC) models that proved able to predict the nanostructure evolution under neutron irradiation in both RPV and F/M steels. These were modelled, respectively, in terms of Fe-C-MnNi and Fe-C-Cr alloys, but the model was also validated against data obtained on a real RPV steel coming from the surveillance programme of the Ringhals Swedish nuclear power plant. The effects of the substitutional solutes of interest were introduced in our OKMC model under the simplifying assumptions of ‘‘grey alloy’’ scheme, i.e. they were not explicitly introduced in the model, which therefore cannot describe their redistribution under irradiation, but their effect was translated into modified parameters for the mobility of defect clusters. The possible origin of low temperature radiation hardening (and subsequent embrittlement) was also investigated and the models strongly supported the hypothesis that solute clusters segregate on immobile interstitial loops, which act therefore as heterogeneous nucleation sites for the formation of the NiSiPCr- and MnNi-enriched cluster populations experimentally, as observed with atom probe tomography in, respectively, F/M and RPV steels. In other words, the so-called matrix damage would be intimately associated with solute atom clusters and precipitates which increase their stability and reduce their mobility: their ultimate effect is reflected in an alteration of the macroscopic mechanical properties of the investigated alloys. Throughout all our work the obtained results have been systematically validated on existing experimental data, in a process of continuous improvement of the physical hypotheses adopted.

Le recours aux modèles et à la modélisation est mentionné dans la documentation scientifique comme un moyen de favoriser la mise en oeuvre de pratiques d'enseignement-apprentissage constructivistes pour pallier les difficultés d'apprentissage en sciences. L'étude préalable du rapport des enseignantes et des enseignants aux modèles et à la modélisation est alors pertinente pour comprendre leurs pratiques d'enseignement et identifier des éléments dont la prise en compte dans les formations initiale et disciplinaire peut contribuer au développement d'un enseignement constructiviste des sciences. Plusieurs recherches ont porté sur ces conceptions sans faire de distinction selon les matières enseignées, telles la physique, la chimie ou la biologie, alors que les modèles ne sont pas forcément utilisés ou compris de la même manière dans ces différentes disciplines. Notre recherche s'est intéressée aux conceptions d'enseignantes et d'enseignants de biologie au secondaire au sujet des modèles scientifiques, de quelques formes de représentations de ces modèles ainsi que de leurs modes d'utilisation en classe. Les résultats, que nous avons obtenus au moyen d'une série d'entrevues semi-dirigées, indiquent que globalement leurs conceptions au sujet des modèles sont compatibles avec celle scientifiquement admise, mais varient quant aux formes de représentations des modèles. L'examen de ces conceptions témoigne d'une connaissance limitée des modèles et variable selon la matière enseignée. Le niveau d'études, la formation préalable, l'expérience en enseignement et un possible cloisonnement des matières pourraient expliquer les différentes conceptions identifiées. En outre, des difficultés temporelles, conceptuelles et techniques peuvent friner leurs tentatives de modélisation avec les élèves. Toutefois, nos résultats accréditent l'hypothèse que les conceptions des enseignantes et des enseignants eux-mêmes au sujet des modèles, de leurs formes de représentation et de leur approche constructiviste en enseignement représentent les plus grands obstacles à la construction des modèles en classe.

La modélisation est essentielle pour comprendre et reproduire les phénomènes physiques dominants ayant lieu dans des écoulements turbulents naturels (atmosphériques, océaniques). En effet, la dynamique des écoulements géophysiques résulte d'interactions complexes à des échelles et intensités variées, et sur des temps différents. La description précise de tels écoulements est pour le moment hors de portée des simulations numériques directes, surtout à cause des limitations en nombre de Reynolds. C'est pourquoi dans cette thèse on s'attaque à la modélisation de la turbulence homogène avec le formalisme spectral de l'approximation EDQNM. Ceci nous permet d'obtenir des résultats rapidement en termes de ressources numériques à très grands nombres de Reynolds, et ainsi d'étudier séparément la plupart des mécanismes en jeu dans les écoulements turbulents naturels, à savoir le cisaillement, le gradient de température, la stratification, l'hélicité, et des combinaisons de ces éléments. On procède en deux étapes: tout d'abord, l'EDQNM permet de fermer les équations des moments d'ordre 2, et ensuite l'anisotropie est modélisée grâce à des tenseurs moyennés sphériquement. Cette méthode est appliquée aux différentes configurations mentionnées ci-dessus, nous permet de proposer de nouveaux résultats et de les valider numériquement à grands nombres de Reynolds. Parmi les points les plus importants, nous nous sommes concentrés sur (i) la prédiction des lois de croissance et décroissance de quantités telles que l'énergie cinétique, la variance scalaire et l'hélicité; (ii) la détermination des comportements spectraux; et (iii) la distribution d'anisotropie échelle par échelle. / Modelling is essential to understand and reproduce the dominant physical mechanisms occurring in natural turbulent flows such as atmospheric and oceanic ones. Indeed, the dynamics of geophysical flows results of multiple complex processes interacting with each others, at various scales, intensities, and on different characteristic times. The fine description of such flows is currently out of reach of direct numerical simulations, notably because of Reynolds numbers limitations. Consequently, we address in this thesis the modelling of homogeneous turbulence, using the spectral formalism of the eddy-damped quasi-normal Markovian (EDQNM) approximation. This first allows us to obtain results rapidly in terms of computational resources at very large Reynolds numbers, and thus to investigate separately some of the fundamental mechanisms at stake in natural turbulent flows, namely shear, mean temperature gradient, stratification, helicity, and combinations of these processes. In this framework, a two-step approach is considered: first, EDQNM is used to close the non-linear terms in the second-order moments equations, and anisotropy is then modelled through spherically-averaged tensors. This methodology is applied to the various configurations mentioned above, permits to propose new theoretical results, and to assess them numerically at large Reynolds numbers. Among the most important findings, we focused on (i) the prediction of the decay and growth laws of crucial one-point statistics such as the kinetic energy, the scalar variance, and helicity; (ii) the determination of spectral scalings; and (iii) the scale by scale distribution of anisotropy.

L’azote (N) est un élément indispensable pour les végétaux mais reste limitant dans la majorité des écosystèmes forestiers. La source principale de N pour les arbres provient de la décomposition des litières. Afin d’étudier le devenir du N des litières dans l’écosystème forestier, une expérience de marquage d’une cohorte de litière enrichie en azote 15 (15N) a été mise en place dans dix hêtraies européennes. Sur le moyen terme (une décennie), les profils 15N dans le sol montrent une faible incorporation du marqueur en profondeur (>10cm). Dans les feuilles, la dynamique d’incorporation du 15N par l’arbre témoigne de la mise à disposition du N des litières puis de sa stabilisation dans les agrégats du sol. L’utilisation du N dérivé des litières par l’arbre dépend des besoins en N des arbres mais aussi du type d’humus. Les arbres sous moder sont capable d’utiliser le N dérivé des litières directement depuis l’horizon organique. Environ 25 % des besoins totaux en N de l’arbre sont couverts par la réallocation du N, et 75 % sont issus du prélèvement du N dans le sol. / Nitrogen ( N) is an essential element for plants but is limited in most forest ecosystems. The main source of N for trees comes from the decomposition of litter. The fate of litter N on the long term (ten years) in beech forest was followed using one cohort of litter enriched in nitrogen-15 (15N) in ten sites in Europe. 15N in the soil profile showed a low incorporation of litter-derived N deeper than 10 cm. Leaf 15N dynamics illustrate the litter N availability and then its stabilization in the soil aggregates.The use of litter derived N in the tree depends on the needs of N by trees but also on the type of humus. At moder humus sites, trees are able to mine the litter-derived N directly from the organic horizon. Approximately 25 % of the total N requirement of the trees are covered by the N reallocation, and 75 % comes from soil N uptake.

La modélisation numérique de la propagation des failles décrochantes est essentielle à notre compréhension du comportement mécanique long terme de la lithosphère. Dans un premier temps, nous avons cherché à reproduire la déformation intersismique autour d'une faille préexistante. Cette modélisation du chargement intersismique a permis de déterminer quelle rhéologie et quel type de chargement, justifiables mécaniquement et géologiquement, étaient nécessaires à l'obtention d'une déformation de surface cohérente avec les mesures géodésiques. Pour un déplacement latéral tel qu'imposé par le mouvement long terme des plaques, il faut imposer une hétérogénéité structurelle à l'aide d'une zone de faible viscosité sous la partie bloquée de la faille. Dans un second temps, nous avons étudié la propagation d'une faille à l'aide de la mécanique de l'endommagement. Cette approche permet de prédire l'évolution de l'endommagement sans prédéterminer le chemin de propagation, et de faire le lien avec la théorie de Griffith. Nous avons amélioré cette approche en introduisant une contrainte tangentielle résiduelle sur la faille qui s'apparente à l'introduction de friction sur les plans de faille, et en permettant l'utilisation, pour la première fois en mécanique de l'endommagement, de méthodes de résolutions mixtes. Ces développements ont permis de décrire dans des cas simples l'évolution d'une faille, ouvrant ainsi la voie à l'utilisation de nouveaux types de modèles de croissance et de propagation des failles. / Numerical modelling of strike-slip faults propagation is essential to understand the long-term mechanical behavior of the lithosphere. Firstly, we aim to reproduce the interseismic deformation of a preexisting fault. The modelling of this interseismic regime enables us to determine what rheology and loading, both mechanically and geologically justifiable, are required to obtain a surface deformation that matches the geodetic measurements. When the lateral displacement is due to the long-term plates motion, a structural heterogeneity has to be imposed by the mean of a weak viscous zone beneath the locked fault. Secondly, we study fault propagation by using damage mechanics. This approach allows the prediction of damage evolution, without previously knowing the propagation path, and makes the link with the Griffith's theory. We improve this approach with the introduction of a residual tangential stress on the fault that is similar to fault planes friction, and we use mixed resolution methods, which represents a novelty in damage mechanics. These developments make possible the description of fault evolution for simple cases, and open the way towards the use of new growth and propagation fault models

Les propriétés statiques, dynamiques, rhéologiques et la cinétique de scissions et recombinaisons de micelles linéaires auto-assemblées sont étudiées à l'équilibre et sous-écoulement par simulations sur ordinateur, en utilisant un modèle mésoscopique nouveau. Nous représentons les micelles comme des séquences linéaires de billes browniennes dont l'évolution spatio-temporelle est gouvernée par la dynamique de Langevin. Un algorithme de Monte-Carlo contrôle l'ouverture des liens ou la fusion de deux chaînes par les bouts. Un paramètre cinétique o, qui modélise l'effet d'une barrière le long d'un chemin de réaction, est introduit dans notre modèle. A l'équilibre, nous nous concentrons sur les mécanismes de scission/recombinaison aux temps long et court. Nos résultats montrent que pour les temps plus grands que le temps de vie d'une chaîne moyenne, la cinétique est en accord avec le modèle champ-moyen de Cates. L'étude de fonctions de relaxation macroscopique confirme que nos constantes cinétiques effectives obtenues aux temps longs sont pertinentes pour ces relaxations. Pour la situation hors équilibre, nous étudions les effets du couplage entre un écoulement de cisaillement et la cinétique de scission et recombinaison sur les propriétés structurales et rhéologiques du système micellaire. Nous nous plaçons dans un régime semi-dilué et dynamiquement 'unentangled'. Le paramètre o est choisi de façon à ce que la durée de vie d'une chaîne moyenne soit plus courte que son temps de relaxation de Rouse le plus long. Nos analyses font apparaître une longueur dynamique A, le fiagrnent de chaîne dont la durée de vie TA est égale à son temps de Rouse. Nous trouvons que les propriétés telles que la rhéo-fluidification, l'orientation des chaînes et l'étirement des liens sont des fonctions du taux de cisaillement réduit PA= YT* , alors que la longueur moyenne des micelles est une fonction décroissante du taux de cisaillement, indépendamment de la barrière du processus scission/recombinaison / Statics, Dynamics, and Rheological properties of Micellar solutions by Computer Simulation Statics, dynamics, rheology and scission-recombination kinetics of self-assembling linear micelles are investigated at equlibrium state and under shear flow by computer simulations using a newly proposed mesoscopic model. We model the micelles as linear sequences of Brownian beads whose space-time evolution is governed by Langevin dynamics. A Monte Carlo algorithm controls the opening of a bond or the chain-end fusion. A kinetic parameter o, modelling the effect of a potential barrier along a kinetic path, is introduced in our model. For equilibrium state we focus on the analysis of short and long time behaviors of the scission and recombination mechanisms. Our results show that at time scales larger than the life time of the average chain length, the kinetics is in agreement with the mean-field kinetics model of Cates. By studying macroscopic relaxation phenomena such as the average micelle length evolution after a T-jump, the monomer diffusion, and the zero shear relaxation function, we confirm that the effective kinetic constants found are indeed the relevant parameters when macroscopic relaxation is coupled to the kinetics of micelles. For the non-equilibrium situation, we study the coupled effects of the shear flow and the scissionrecombination kinetics, on the structural and rheological properties of this micellar system. Our study is performed in semi-dilute and dynamically unentangled regime conditions. The explored parameter o range is chosen in order for the life time of the average size chain to remain shorter than its intrinsic (Rouse) longest relaxation time. Central to our analysis is the concept of dynamical unit of size A, the chain fiagrnent for which the life time TA and the Rouse time are equal. Shear thinning, chain orientation and bond stretching are found to depend upon the reduced shear rate P1\=y~A while the average micelle size is found to decrease with increasing shear rate, independently of the height of the barrier of the scission-recombination process

La réduction catalytique sélective des oxydes d'azote au moyen de l'urée ou de l'ammoniac est une technique établisse, pour une conversion efficace des NOx, qui existent dans le gaz d'échappement des moteurs Diesel mobiles. Parmi les formulations catalytiques proposées, les zéolithes échangées au fer sont considérés comme des catalyseurs SCR prometteuse, grâce à leur performance deNOx élevée, sur une large gamme des températures, ainsi qu'une résistance importante au vieillissement hydrothermique. Dans le cadre de cette thèse, une étude expérimentale détaillée des voies réactionnelles de la SCR sur les zéolithes au Fe a été réalisée de façon à développer un modèle phénoménologique et macro-cinétique, appliqué pour la simulation des systèmes de post-traitement des moteurs Diesel. Des expériences dédiées qui comprenaient des mesures de spectroscopie IR ont été réalisées sur les catalyseurs H- et Fe-BEA, synthétisés en laboratoire pour élucider les aspects mécanistiques liées à la formation et aux interactions de l'ammoniac et d'oxydes d'azote (notamment des espèces nitrates et nitrites) sur des sites catalytiques acides et redox. La formation et la décomposition de NH4NO3 ont été particulièrement préoccupantes. Basé sur des résultats expérimentaux un modèle cinétique multi-site a été développé. Le modèle a été capable de prendre en compte la variation de l'acidité de surface, les interactions des espèces gazeux du NH3 et des NOx avec des sites métalliques, ainsi que des phénomènes de physisorption à basse température. Le modèle a été validé avec succès par la simulation des expériences réalisées sur divers échantillons, CE qui souligne l'intérêt de l'approche suivie. / Selective Catalytic Reduction of nitrogen oxides through urea or ammonia is a well established technique for the efficient abatement of NOx, from mobile diesel engines exhaust stream. Among catalytic formulations proposed, Fe-exchanged zeolites are considered as promising SCR catalysts, due to their enhanced deNOx performance over a broad range of operating temperatures, as well as significant resistance to hydrothermal ageing. In the present work, detailed experimental investigation of the NH3-SCR reactions network over Fe zeolites was performed in order to develop a phenomenological macro-kinetic model oriented to diesel engines aftertreatment system simulation. Dedicated experiments, including IR spectroscopy measurements were performed over lab-synthesized H- and Fe-BEA catalysts in order to unravel mechanistic aspects related to the formation and interactions of ammonia and nitrogen oxides adspecies (namely nitrate and nitrites) on acidic and redox catalytic sites. Formation and decomposition of NH4NO3 was of particular concern. Based on experimental results a multi-site kinetic model was developed, which was able to account for surface acidity variation, interactions of gaseous NH3 and NOx species with iron sites, as well as low temperature physisorption phenomena. The model was successfully validated through simulation of experiments performed over various Fe-zeolite samples, applying a broad range of operating conditions, which showed the interest of the followed approach.

Les polluants atmosphériques, qu’ils soient de sources naturelles ou anthropiques, peuvent être transformés et transportés sur de longues distances. Les retombées atmosphériques constituent des puits définitifs pour ces polluants atmosphériques. Ce travail de recherche s’inscrit dans le cadre du projet SESAME du programme PRIMEQUAL « Pollution longue distance ». Les objectifs sont d’identifier les déterminants de l’évolution des retombées atmosphériques acidifiantes et eutrophisantes et d’évaluer l’efficacité des politiques de réduction en vigueur. Les méthodologies reposent sur des approches statistiques appliquées à une compilation de bases de données issues de trois observatoires MERA, CATAENAT et BAPMoN. L’utilisation de procédures de consolidation et de comparaison a permis d’optimiser la base de données pour envisager des analyses descriptives en toute objectivité. La détermination de la composition chimique des espèces ioniques majeures a permis de distinguer les déterminants des contrastes géographiques et saisonniers de la chimie des retombées atmosphériques. L’étude des changements temporels annuels et saisonniers par la méthode de Sen Mann Kendall a mis en évidence des liens pouvant être linéaires ou non linéaires entre les évolutions de la composition chimique des retombées atmosphériques et leurs déterminants. L’expertise des déterminants a contribué à la création pas à pas d’un modèle statistique en vue de proposer des prévisions des retombées atmosphériques au site de mesure à l’horizon de quelques années (suivant différents scénarios climatiques). Des valeurs ponctuelles de dépassement de charges critiques et leur évolution face aux dépôts mesurés et projetés ont été évaluées. Enfin, la densité de sites a permis d’obtenir une spatialisation des retombées atmosphériques, qui a pu être comparée aux cartes de dépôts modélisées par le modèle EMEP. / Atmospheric pollutants emitted by both natural and anthropogenic sources are precursors of secondary compounds such as acidic gases and aerosols, which can be transported over long distances. Atmospheric deposition is the ultimate sink for these pollutants and can lead to acidification and eutrophication issues. This research was conducted within the SESAME project as part of the PRIMEQUAL long distance pollution program. The objectives seek to identify determinants of atmospheric deposition changes and to evaluate policies of pollutant reductions. Statistical approaches were applied to measurement databases of three French networks: MERA, CATAENAT and BAPMoN. Several procedures were developed to consolidate and to compare the datasets for a descriptive analysis. Geographical and seasonal determinants were identified based on the chemistry of concentrations and deposition fluxes. The trends resulting from Sen Mann Kendall tests highlighted that the changes in chemical compositions of atmospheric deposition can be either linear or nonlinear with the identified determinants. The findings were used to develop a predictive model to forecast the deposition fluxes of sulfur and nitrogen species under different climate scenarios and over several decades. The exceedances of critical loads and their changes were evaluated from current and forecasted deposition fluxes. Finally, the density of measurement sites over the French area was strong enough to create maps of deposition fluxes, which were then compared to the EMEP model outputs.

Comme bien d'autres grandes métropoles construites sur des formations géologiques aquifères, Paris présente une situation hydrogéologique complexe et perturbée.L'imperméabilisation du sol modifie notablement la recharge naturelle des nappes qui se fait alors, principalement, par les fuites des réseaux (voirie, eau non potable, eaux usée …). L'urbanisation intense du sous-sol, qui représente plus de 30 % de la surface parisienne (hors bois de Vincennes et Boulogne), influence les écoulements des aquifères. A cela, s'ajoute le rôle des 20 à 30 millions de mètres cubes d'eau qui son prélevés annuellement dans les nappes phréatiques.Les aquifères parisiens sont surveillés par l'Inspection Générale des Carrières de Paris (Ville de Paris) grâce à un réseau de 305 piézomètres, qui suit en particulier les ouvrages à proximité de la Seine pour anticiper les conséquences d'une crue. La réalisation d'un modèle hydrogéologique du sous-sol parisien permettant d'évaluer les effets de tels évènements et tenant compte des spécificités de l'anthropisation présentait donc un intérêt.Dans cet objectif, un SIG en trois dimensions du sous-sol parisien a, en premier lieu, été construit. Il intègre plusieurs éléments :- La géologie parisienne depuis la craie campanienne jusqu'aux formations tertiaires. 20 000 sondages ont été exploités pour cela.- Les infrastructures souterraines : 20 000 parkings résidentiels et constructions en sous-sol, plus de 140 parcs de stationnement concédés par la ville de Paris, 197 km de galeries du métro souterrain et 46 km du Réseau Express Régional.- Les prélèvements d'eau, en se basant sur les prélèvements annuels déclarés (12 millions de mètres cubes).- Les données piézométriques issues de 669 piézomètres répartis sur Paris et des communes avoisinantes.Avant ce travail, il n'existait aucune représentation synthétique de l'organisation générale du sous-sol de la capitale et de ses aménagements.Un modèle mathématique hydrodynamique s'étendant sur l'emprise de la capitale et de la proche banlieue a ensuite été construit. Ce modèle multicouche comporte 11 couches alternant 6 niveaux aquifères et 5 niveaux semi perméables, discrétisées selon 607 171 mailles carrées de 50 mètres de coté. Au stade actuel, ce modèle a permis la réalisation de simulations exploratoires reproduisant en régime permanent la structure de l'écoulement dans les conditions d'aujourd'hui. / Paris, as many large cities built on aquifers, shows a disturbed hydrogeological environment.Land waterproofing reduces the natural recharge of the water table, then mainly driven by sewer network leakage. Underground intensive urbanisation covering more than 30 % of the Parisian surfaces (without Vincennes and Boulogne woods) influences groundwater flow. In addition, had to be taken in consideration the 20-30 millions cubic meters of water yearly pumped from the aquifer system.Parisian aquifers are surveyed by IGC (Inspection General des Carrières) from the city of Paris, thanks to a 305 piezometers network, in particular close to the Seine river, in order to monitor and anticipate a flood event. In order to appraise such events a hydrogeological model of the Parisian underground had to be setup, including anthropic specification.To achieve this model a 3D GIS of the Parisian underground has been done, integrating several topics:- Parisian geology, including Campanian chalk and tertiary formation. 20 000 geological drillings have been used.- Underground infrastructures: 20 000 residential car parks and cellars more than 140 public car parks conceded by the city of Paris, 197 km of underground subway and 46 km of RER.- Water pumping. Only declared samples (12 millions cubic meters)- Piezometric data from 669 piezometers over Paris and neighbouring townsBefore this thesis, no synthetic representation of Paris geology and underground has been done.In a second time a mathematical model including Paris and suburbs has been built. This multilayer model consists of 11 layers interbedding, 6 aquifers level and 5 semi-permeable discretized through 607 171 square mesh (50 m side). Until now, this mathematical model allowed exploratory simulation in steady state showing flow structure of the present day hydrogeological situation.