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.

Les liaisons bimétalliques entre acier faiblement allié et acier inoxydable sont nombreuses au sein des réacteurs nucléaires français où elles assurent la connexion entre les principaux composants et les tuyauteries du circuit primaire. Le revêtement interne de cuve réalisé par soudage feuillard-flux est un autre cas de soudure dissimilaire dont le rôle est d'assurer une bonne protection contre la corrosion. Ce travail de thèse a notamment pour but de comprendre la genèse des microstructures complexes se formant à l'interface entre les deux aciers pendant le soudage. Il étudie d'autre part leur évolution durant le traitement thermique post-soudage à 610°C ainsi que les conséquences de ces transformations sur le comportement mécanique du joint soudé. Partant du métal de base, on rencontre successivement une fine bande de martensite, une zone purement austénitique puis la microstructure biphasée δ/γ de l'acier inoxydable. Des techniques de microscopie (MEB, EDS, EBSD) combinées à des calculs thermo-cinétiques (modèle de Scheil-Gulliver, surfusion en pointe de dendrite) ont permis d'expliquer le gradient de microstructure et les raisons des transitions de phases observées. Au cours du traitement thermique de détensionnement à 610°C, le gradient de potentiel chimique du carbone à travers l'interface de fusion cause la dissolution de la cémentite et la croissance des grains du côté faiblement allié. On observe également la diffusion du carbone à travers l'interface et la précipitation de carbures riches en Cr dans le liseré martensitique et la zone austénitique. Une caractérisation détaillée des profils de composition et de la précipitation a été réalisée à différentes échelles (depuis le millimètre jusqu'au niveau atomique). Un modèle mésoscopique, s'appuyant sur des bases de données thermodynamiques et cinétiques, a été développé pour coupler la diffusion à longue distance dans un milieu multi-constitué à la germination et croissance des précipités (approche de type Kampmann-Wagner). Il a permis de prévoir la teneur en carbone ainsi que la fraction volumique et la distribution de taille des précipités en fonction de la distance à l'interface. Les conséquences de la forte variabilité de microstructure sur le comportement mécanique local ont été analysées dans la dernière partie de ce travail, en particulier les aspects de déformation localisée et de rupture ductile. Des lois de comportement élasto-plastique ont été déterminées pour chacune des régions de l'assemblage à l'état détensionné. L'étude des mécanismes de rupture ductile dans les zones les plus faibles de la soudure, c'est-à-dire le métal de base décarburé et les couches de revêtement austénitique a donné lieu à des observations in-situ et une modélisation de l'endommagement. / Dissimilar welds between low-alloy steel and stainless steel are numerous within the French nuclear power plants where they enable connecting the main components to the primary circuit pipes. The internal cladding (in stainless steel) of the pressure vessel (in bainitic steel) made by submerged arc welding is another case of dissimilar weld whose goal is the protection against corrosion. This PhD work aims at understanding the complex microstructures which form at the interface between both steels during welding together with their evolution during the post-weld heat-treatment at 610°C and their consequences on the mechanical behavior of the welded assembly. Starting from the base metal, one meets successively a thin layer of martensite, a fully austenitic zone and the two-phase δ/γ microstructure of the stainless steel. Microscopy techniques (SEM, EDS, EBSD) combined with thermo-kinetics calculations (Scheil-Gulliver model, dendrite tip undercooling) have allowed explaining the graded microstructure and the reasons for the observed phase transitions. During the post-weld heat-treatment, the large gradient of carbon chemical potential across the fusion line leads to cementite dissolution and grain growth on the low-alloyed side. Carbon diffusion through the interface and Cr-rich carbides precipitation in both the martensitic layer and the austenitic weld have also been observed. An in-depth characterization has been performed at different scales (from the millimeter to the atomic level) to quantify the extent of carbon diffusion and carbides precipitation during the phase transformations. A mesoscopic thermodynamic and kinetic model based on Calphad databases has been developed to fully couple long-range diffusion in a multi-component system with precipitates nucleation and growth (Numerical Kampmann-Wagner approach). It allowed a prediction of the carbon content, volume fraction and size distribution of the precipitates at any distance from the fusion line. The consequences of the high variability of microstructures on the local mechanical behavior have been examined in the last part of this work, in particular the localization of deformation and the ductile failure. Elasto-plastic constitutive laws were determined for each region of the dissimilar weld in the heat-treated state. Ductile failure mechanisms in the weak zones of the weld, namely the decarburized base metal and the stainless steel cladding layers, were investigated through in-situ observations and damage modeling.

Ces travaux ont permis d'approfondir les mécanismes de (dé)lithiation et de vieillissement dans des électrodes à base de silicium et de LiFePO4 pour accumulateurs Li-ion à partir d'observations multi-échelles. Des cartographies de phases, autant à l'échelle de la particule qu'à l'échelle de l'électrode, ont été menées par microscopie électronique mettant en évidence de fortes hétérogénéités. Pour le silicium, la mise en place de cartographie unique par STEM/EELS, s'appuyant sur une base de données des pertes faibles d'alliages sensibles à l'air et au faisceau d'électrons, a permis de comprendre les mécanismes de lithiation à l'échelle du nanomètre. L'étude de la première lithiation a montré des différences de mécanismes de réaction avec le lithium suivant deux facteurs : la taille des particules et les défauts au sein de celles-ci. Il a été observé une composition d'alliage LixSi plus faible pour les nanoparticules que pour les microparticules. Les défauts dus notamment au broyage constituent des sites préférentiels de lithiation. En vieillissement, les nanoparticules subissent de profonds changements structuraux et morphologiques, passant d'un état sphérique cristallin (50 nm) à un réseau de fils amorphe (5-10 nm d'épaisseur) contenu dans une matrice de SEI. Pour le LiFePO4, il a été clairement montré, par la combinaison de plusieurs techniques de microscopies électroniques (diffraction des électrons en précession, EFSD : Electron Forward Scattering Diffraction, EFTEM), que les particules de taille nanométrique (100-200 nm) étaient soit entièrement lithiées soit entièrement délithiées à l'équilibre thermodynamique. De fortes hétérogénéités ont été observées dans les électrodes fines comme dans les électrodes épaisses. A l'échelle des particules, l'analyse statistique de plus de 64000 particules a montré que les plus petites particules se délithient en premier. A l'échelle de l'agglomérat, les cartographies de phases ont révélé un mécanisme « cœur-coquille » : la réaction débute de la surface vers le centre des agglomérats. A l'échelle de l'électrode, le front de propagation de phase se déplace suivant des chemins préférentiels de plus grandes porosités de la surface de l'électrode vers le collecteur de courant. La conductivité ionique au sein de nos électrodes est le facteur limitant. / This work aimed at better understanding the (de)lithiation and aging mechanisms in LiFePO4 and silicon-based electrodes for Li-ion batteries from multiscale investigations. Phase mapping was performed by electron microscopy at the particle scale and at the electrode scale. This highlights some strong heterogeneities. The silicon study has shown some different lithium reaction mechanisms following two effects: particle size and crystalline defects. A smaller lithium amount in LixSi alloy was highlighted for the nanoparticles rather than for the microparticles. The defects mainly due to milling are preferential sites for the lithiation. In aging, the nanoparticles have undergone structural and morphological changes. The pristine crystalline spherical shape (50 nm) was transformed into an amorphous wire network (5-10 nm of thickness) contained in a SEI matrix. Thanks to a combination of electron microscopy techniques (precession electron diffraction, Electron Forward Scattering Diffraction, EFTEM), it was clearly shown that the LiFePO4 particles (100-200 nm) are either fully lithiated or fully delithiated at the thermodynamic equilibrium. Strong heterogeneities were observed in the thin and thick electrodes. At the nanoscale, the statistical analysis of 64000 particles unambiguously shows that the small particles delithiate in first. At the mesoscale, the phase maps reveal a core-shell mechanism at the scale of the agglomerates, from the surface to the center of these agglomerates. At the electrode scale, the phase front would move following preferential paths into the higher porosity from the surface in contact with electrolyte toward the current collector. The electrode ionic conductivity is the limiting parameter.

L'étude a pour but de caractériser les facteurs de contrôles de la distribution des propriétés pétrophysiques et mécaniques dans les réservoirs carbonatés. Pour y répondre, il a été entrepris une démarche intégrant des analyses sédimento-diagénétiques (sédimentologiques, pétrographiques, géochimiques…) et de la fracturation (stratigraphie mécanique et stratigraphie de fracture, hiérarchisation, chronologie…). L'acquisition des données a été réalisée sur un analogue de terrain d'un réservoir carbonaté fracturé : la Formation Madison, d'âge Carbonifère inférieur (357-340 Ma), affleurant dans le Bassin de Bighorn (Wyoming, USA) et qui est aussi un réservoir exploité en subsurface. Cette acquisition a été réalisée dans cinq sites : Wind River Canyon, Shell Canyon, Sheep Mountain, Shoshone Canyon, et Clark's Fork Canyon (selon une polarité paléogéographique proximale - distale).Dans le Wyoming (USA), la Formation Madison est une série carbonatée, atteignant une épaisseur de 340 m, formée de sept séquences de dépôt basse fréquence (SBF). Les deux premières (SBF1 et 2) se sont formées sur une rampe qui, en conséquence de progradations extensives, s'est aplanie progressivement pour former une plate-forme sur laquelle se sont déposées les cinq autres séquences (SBF3 à 7). SBF1 à 3, d'âge Kinderhookien à Osagéen inférieur, se sont formées sous climat aride à la faveur d'une subsidence tectonique générale compensant les chutes eustatiques épisodiques. Les cortèges de rétrogradation se caractérisent par le développement étendu des milieux inter- à supratidaux et des processus de précipitation d'évaporites et dolomitisation. / The purpose of this study is to characterize factors controlling the distribution of the petrophysical and mechanical properties in carbonate reservoirs. To do so, a pluridisciplinary approach integrating sedimento-diagenetic and fracturing analyses on a carbonate reservoir analogue: he Madison Formation, (Lower Carboniferous, 357-340 My), outcropping in the Bighorn Basin (Wyoming, USA) and that is also a subsurface reservoir.The Madison Formation is a 340 m thick carbonate series composed of seven low frequency depositional sequences (LFS). The first two (LFS 1 and 2) formed on a prograding ramp passing upward into a vast platform on which the other five LFS deposited (LFS3 à 7). LFS 1 to 3 (Kinderhookian to lower Osagean) deposited under arid conditions during general subsidence that balanced eustatic falls. Retrograding system tracts are characterized by the development of supratidal to intertidal environment dominated by evaporite precipitations and carbonate dolomitization. Prograding system tracts are mainly mainly by early-lithified grainstones. LFS 4 to 7 (Osagean) deposited under humid conditions (glacial conditions in high-latitudes) that contributed to a decrease in evaporite precipitations and carbonate dolomitization but that involved karstifications on tops of LFS 4 to 7 during uplift episodes and eustatic falls.

Les injections de CO2 à l'échelle industrielle, dans les aquifères salins profonds, vont affecter le système naturel en induisant des perturbations des gradients de pression à court et à moyen terme.Des études prévisionnelles sont nécessaires pour évaluer les risques de contamination (fuites de CO2, déplacement des fluides natifs) et les risques d'interférences entre projets d'injection ou entre utilisations du sous-sol.Les aquifères salins étant généralement peu caractérisés, les incertitudes géologiques sont à considérer lors de l'étude de faisabilité du stockage et des risques associés puisque les paramètres géologiques influencent la réponse du système à l'injection.Nous nous sommes intéressés aux incertitudes résultantes en termes de prévisions de perturbations de pression et de migration de CO2 et à leurs conséquences sur la faisabilité des projets de stockage.Dans un premier temps, les incertitudes de modélisation (changement d'échelle, résolution de la variabilité spatiale des propriétés pétrophysiques) et géologiques (propriétés pétrophysiques de la formation d'injection et de la couverture) ont été étudiées sur des modèles 2D conceptuels. L'objectif étant de balayer les champ d'incertitudes pour des modèles peu coûteux en temps de calcul, pour ensuite, réduire les évaluations à effectuer dans le cadre de la modélisation 3D régionale du système souterrain et y appliquer des méthodes simplifiées, validées en deux dimensions.Des centaines de réalisations stochastiques sont utilisées pour évaluer l'influence de la variabilité spatiale de la perméabilité. Pour limiter le nombre de simulations d'écoulement à effectuer, des méthodes de sélection de réalisations, à partir de "proxy-response" (i.e. approximation de la réponse par une méthode de calcul simplifiée) ont été testées et validées.Ensuite, les modèles 3D sont construits à partir des données d'un modèle hydrogéologique du bassin parisien. Différents scénarios d'injection sont envisagés. La sensibilité de la réponse est étudiée principalement par rapport à la variabilité spatiale de la perméabilité et à la compressibilité des pores. Cette dernière étape permet de mieux appréhender les risques d'interférences en fonction des incertitudes majeures, d'une part des paramètres géologiques, et d'autre part des paramètres physiques liés à l'injection. / This thesis studies the regional-scale response of an aquifer system to a massive CO$_2$ injection. Industrial-scale CO$_2$ injections into deep saline aquifers affect natural groundwater systems by generating short-term to medium-term pressure gradient perturbations. To evaluate contamination risks and interference risks between injection projects or other uses of underground space, modelling studies become necessary. The geological parameters of underground formations are also to be taken into consideration as they certainly influence the injection reponse. But, saline aquifers are generally poorly-characterized which adds uncertainties to an already complex system. This thesis aims to explore uncertainties in pressure perturbations and CO$_2$ migration predictions, and their consequences in terms of CO$_2$ storage feasibility studies. Firstly, modelling and geological uncertainties have been tested on 2D conceptual models. This step, based on simpler models than 3D ones, allows a fast uncertainties discrimination and save computational time. Hundreds of stochastic realizations are generated to define the influence of permeability spatial variability. To limit the number of flow simulations, selection procedures of realizations are applied and tested. Selections are derived from fast-calculations methods called "`proxy-response"'. Secondly, once these methods have been 2D tested and validated, and once a number of uncertainties have been eliminated, these methods and related ones are applied to the underground system 3D modelling. The 3D models have been built based on available data from an existing Paris Basin hydrogeological model. Several injection scenarios have been considered and tested. Permeability spatial variability and pore compressibility are the two main parameters chosen to evaluate the injection response. This last step allows a better definition of interference risks between the major uncertainties from geological parameters and injection-related physical parameters.

Les polymères amorphes présentent des hétérogénéités dynamiques à l'échelle de 3 à 5nm proche de la température de transition vitreuse. Leur contribution aux propriétés mécaniques est encore mal connue. Pour ce faire, nous considérons un modèle simple dans lequel chaque hétérogénéité dynamique possède un unique temps de relaxation distribué aléatoirement selon une loi log-normale. Un modèle coarse-grained à l’échelle d’une hétérogénéité dynamique est alors développé dans le cadre de l'approximation des milieux continus. Ainsi, des simulations par éléments finis nous permettent de calculer les propriétés mécaniques macroscopiques des polymères en tenant compte des effets de couplage mécaniques en hétérogénéités. Nous montrons que la transition vitreuse est pilotée par un réseau de percolation des domaines lents. Nous mettons également en évidence les couplages mécaniques intervenant entre hétérogénéités dynamiques. Dans le cas des films minces de polymères confinés à l’échelle de quelques dizaines de nanomètres, nous montrons que les hétérogénéités dynamiques permettent de rendre compte d’un ralentissement de la dynamique macroscopique du film. En considérant que pour de fortes contraintes, la dynamique locale des polymères massifs est accélérée, nous démontrons que les hétérogénéités dynamiques permettent de rendre compte d’un phénomène de contrainte interne à un système, celle-ci étant la cause d’une partie non récupérable de l’énergie élastique stockée un système lors d’une déformation plastique. / Amorphous polymers present dynamical heterogeneities at the scale of 3 to 5nm near Tg. Their contributions to mechanical properties are still not well known. We thus consider a simple model where each dynamical heterogeneities has its own relaxation time randomly drawn is a log-normal time distribution. A coarse –grained model at the dynamical heterogeneity is then developed in the continuous medium approximation. Finite element simulations are performed to calculate macroscopic mechanical properties of amorphous polymers taking account of mechanical couplings between heterogeneities. We show that the glass transition is controlled by a percolation of slow domains. Mechanical couplings result in viscoelastic spectrum highly narrowed as compared to the microscopic one. In thin films of polymers, we evidence a slowdown of the dynamics of relaxation as compared to the bulk one. Considering that a high applied stress increase the dynamics of polymers, we show that dynamical heterogeneities result in an internal stress network. The latter is a consequence a plastic deformation and result in an internal energy which is not released after an unloading.

L'objet de cette thèse est la construction, l'étude mathématique et numérique de modèles déterministes pour des systèmes Proie-Prédateur en milieu insulaire hétérogène. Il s'agit d'évaluer les effets de l'introduction d'espèces invasives, prédateurs et compétiteurs, sur une population de proies natives. La première partie présente l'étude de modèles faiblement structurés, bas´es sur des systèmes d'E.D.O. singuliers, le dénominateur d'un des termes de réaction pouvant s'annuler. L'analyse mathématique permet d'isoler des conditions d'extinction en temps fini ou de persistance. Dans ce second cas, le comportement en temps long dépend d'hypothèses supplémentaires. Une étude similaire est menée dans le cadre d'une population de proies natives structurée en 2 classes d'âge : juvéniles et adultes. Dans la seconde partie, on étend les modèles précédents au cadre avec structuration en espace, pour prendre en compte les hétérogénéités spatiales du milieu. On obtient des systèmes d'E.D.P. du type Réaction-Diffusion singuliers. Une analyse approfondie donne des critères d'existence globale en temps et d'existence sur un intervalle de temps fini des solutions. Parallèlement, nous mettons en place une méthode numérique du type splitting d'opérateurs dans un but double : valider les modèles spatiaux et étudier des processus d'invasion. Les simulations numériques permettent d'établir le rôle fondamental des proies introduites dans le succès de l'invasion par les prédateurs de colonies isolées de proies natives. Enfin, la structuration discrète en âge pour les proies natives permet d'exhiber des dynamiques oscillatoires.

L’objet de cette thèse est de développer des méthodologies permettant d’identifier la distribution spatiale des valeurs de perméabilité dans des échantillons de roches. Nous avons tout d’abord développé en laboratoire des expériences d’injection de fluide miscible très visqueux dans des échantillons initialement saturés par une saumure peu visqueuse. Pendant l’injection, l’évolution au cours du temps de la pression différentielle entre les deux faces de l’échantillon a été enregistrée par des capteurs de pression. En outre, des mesures scanner ont fourni une carte 3D de la porosité ainsi que des cartes 3D décrivant la distribution spatiale des concentrations dans l’échantillon à différents temps. Nous avons mis en place une méthode d’interprétation donnant directement le profil 1D de la perméabilité le long de la direction d’écoulement à partir de la pression différentielle mesurée au cours du temps. Cette méthode a été validée numériquement et expérimentalement. Puis, afin d’affiner la description de l’agencement des valeurs de perméabilité dans l’échantillon, c’est à dire d’obtenir un modèle 3D de perméabilité représentatif de l’échantillon, nous avons développé une méthodologie itérative de calage des pressions et des concentrations. Cette méthode passe par deux étapes : une optimisation simple pour capturer l’hétérogénéité dans la direction de l’écoulement et une optimisation complexe pour capturer l’hétérogénéité transverse. Cette méthode a été validée à partir de tests numériques. La méthode a été appliquée à deux des expériences d’injection de fluide visqueux. Nous avons pu alors déterminer des modèles de perméabilité capables de reproduire assez bien les données de pression et de concentration acquises pendant l’injection / The objective of this study is to develop new methodologies to identify the spatial distribution of permeability values inside the heterogeneous core samples. We developed laboratory viscous miscible displacements by injecting high viscosity glycerin into the core samples initially saturated by low viscosity brine. The pressure drop across the samples was measured as a function of time until breakthrough. Meanwhile, CT scan measurements provided a 3D porosity map plus several 3D maps of concentration distribution inside the core samples at different times. A simple permeability mapping technique was developed deducing a one-dimensional permeability profile along the flow direction from the measured pressure drop data. The method was validated with both numerical and laboratory experiments. To go beyond one-dimensional characterization of permeability into cores, we developed an iterative process for matching pressure and concentration data. This method consisted of two steps: a simple optimization for capturing the permeability heterogeneity along the flow direction axis and a complex optimization for capturing transversal permeability heterogeneities. The methodology was validated by numerical data. It was also applied to the data collected from two laboratory viscous miscible displacements. We showed that the final 3D permeability models reproduce well the measured pressure drop and concentration data

La présence d'H2S dans un gisement est un important facteur de dévalorisation économique, c'est pourquoi connaitre sa concentration et sa distribution est capital pour optimiser toutes les étapes de l'exploration à la production d'un champ. Dans les conditions de pression et température des réservoirs, l'H2S est beaucoup plus soluble que les hydrocarbures et autres gaz. Le lessivage préférentiel de l'H2S par un aquifère actif pourrait ainsi entrainer la création d'hétérogénéités compositionnelles au cours du temps. L'objectif de la thèse est d'illustrer et quantifier ce processus à l'aide de simulations numériques réalisées avec Hytec, logiciel couplé géochimie-transport diphasique développé par MINES ParisTech. Le lessivage préférentiel de l'H2S est contrôlé par : la solubilité différentielle des gaz qui modifie les quantités relatives de chacun des gaz à proximité de l'aquifère ; le transport aqueux qui exporte les gaz dissous et améliore ainsi la dissolution des gaz sur le long terme ; le transport gazeux qui renouvelle les gaz dissous à l'interface et étend le lessivage du gaz à toute la structure ; les caractéristiques d'un réservoir (type d'aquifère, hétérogénéité des perméabilités). / The H2S occurrence in gas reservoirs is an important factor of economic depreciation. Thus, the knowledge of its content and distribution is a critical parameter when planning field development. Under typical reservoirs conditions of pressure and temperature, H2S is far more soluble than hydrocarbons and other gases. The preferential leaching of H2S by an active aquifer over time could explain the creation of compositional heterogeneities. The thesis aims at illustrating and quantifying this process based on numerical simulations performed with the two-phase transport and geochemical software Hytec, developed by Mines ParisTech (France). This mechanism may be controlled by: Differential solubility of gases, which changes the relative amounts of each gas near the gas-watercontact; Contact with an active aquifer, which can export the dissolved gases thus enhancing dissolution on the long-term; Diffusional transport in the gas phase, which transfers the compositional anomalies farther from the gas-water contact; Geological parameters (type of aquifer, permeability heterogeneities) which can modify the transport scenario.

L’objectif de ce travail est double : d’une part, il s’agit de développer une méthodologie intégrée pour la construction d’un modèle géomécanique ainsi que la représentation des incertitudes associées aux propriétés poro-élastiques des roches constitutives, en exploitant l’ensemble des données disponibles et en s’appuyant de façon cohérente sur les modèles de réservoir statique et dynamique classiquement utilisés par les géologues et les ingénieurs réservoir ; d’autre part, il s’agit d’analyser quel est l’impact des hétérogénéités géologiques, souvent négligées, dans la réponse mécanique du réservoir sollicité par son exploitation, et d’aboutir à des incertitudes sur les champs de contraintes et de déplacements, issues des incertitudes sur ces hétérogénéités et leurs paramètres mécaniques. Pour ce faire, une méthodologie intégrée s’appuyant sur des simulations géostatistiques a été développée. Après la construction du cadre géométrique 3D, le remplissage des propriétés au sein du réservoir suit une démarche de simulations géostatistiques 3D emboîtées, dans laquelle la représentation des hétérogénéités lithologiques conditionne la génération des propriétés poro-élastiques. La démarche consiste ensuite à représenter les incertitudes sur le modèle géomécanique par des ensembles de réalisations géostatistiques dont la réponse mécanique est alors calculée avec un simulateur mécanique aux éléments finis. Les incertitudes sur les champs de contraintes et de déformations sont déduites ensuite des différentes réponses mécaniques obtenues. La démarche a été mise en œuvre sur un réservoir réel, dans un environnement fluvio-deltaïque, produisant en Mer du Nord. Dans ce cadre, il a été démontré que les hétérogénéités du réservoir et leurs incertitudes influencent significativement les calculs des champs de contraintes et de déformations, ainsi que les risques mécaniques de rupture. Des incertitudes sur les quantités mécaniques analysées (premier invariant du tenseur des contraintes et subsidence) ont été aussi estimées / This work has two main objectives. The first one is to develop an integrated methodology allowing to build a 3D geomechanical model and also to image the uncertainties attached to the poro-mechanical properties of the constitutive rocks. This geomechanical model should be based on all related available data and should be consistent with the static and dynamic models, currently built by reservoir geologists and engineers. The second objective is to analyse the impact of geological heterogeneities, which are often neglected, in the mechanical response of the reservoir induced by its exploitation, and furthermore to derive uncertainties on the stress and deformation fields related to the uncertainties on the input properties of the geomechanical model. An integrated methodology based on geostatistical simulations is developed. First, the geometric frame is built; then an approach of embedded stochastic simulations is carried out to infill the different reservoir properties, the lithological description constraining the petrophysical and poro-elastic descriptions. The next step is to generate the mechanical responses of the stochastic realisations, using a finite-element mechanical simulator. The uncertainties on the resulting stress and displacement fields are then deduced from the multiple mechanical responses which are computed. This approach is demonstrated on a real field case, a fluvio-deltaic reservoir in North Sea. It is shown on this example that the reservoir heterogeneities and their uncertainties significantly influence the calculations of stress and strain fields, and also the risks of mechanical failure. Uncertainties on the mechanical quantities under analysis (first invariant of the stress tensor and subsidence) are also derived