Les systèmes turbiditiques chenalisés permettent le transfert de sédiments depuis le plateau continental vers le domaine marin profond de la même façon que les cours d’eau sur les continents. Les dépôts sédimentaires turbiditiques forment parmi les plus grands réservoirs d’hydrocarbures actuellement exploités. MINES ParisTech développe depuis une vingtaine d’années Flumy, un modèle numérique simulant l’évolution d’un système fluviatile méandriforme pour la simulation de réservoirs. L’approche proposée dans ce travail est de transposer Flumy aux systèmes sous-marins profitant de l’analogie entre ces environnements. La géomorphologie des systèmes chenalisés méandriformes fluviatiles et sous-marins est d'abord comparé à titre de données naturelles, et des expériences en laboratoire mettent en évidence les différences dans les écoulements aériens et sous-aquatiques. L'intégration dans Flumy de ces résultats permet de reproduire l'architecture stratigraphique des réservoirs turbiditiques. / Channelized turbidite systems act as sediment transfer routes through the submarine realm similarly to rivers on the continents. The turbidite sedimentary deposits create among the biggest hydrocarbon reservoirs currently produced. MINES ParisTech has developed in the last twenty years Flumy, a numerical model that simulates the evolution of a meandering fluvial system to simulate reservoirs. In this work, we propose to transpose Flumy to submarine systems taking advantage of the analogy between these environments. Using natural data, the geomorphology of channelized meandering fluvial and submarine systems is compared, and laboratory experiments highlight the differences in terms of aerial and subaqueous flow behavior. The integration of the results of these studies allows Flumy to accurately simulate the stratigraphic architecture of channelized turbidite reservoirs.

Alors que l’on utilise couramment de nombreux modèles numériques différents pour la description statique des réservoirs souterrains et des incertitudes associées, les incertitudes sur les écoulements des fluides à l’intérieur de ces réservoirs ne peuvent, pour des raisons de performance, que s’appuyer que sur quelques simulations d’écoulements. Les travaux de cette thèse ont donc pour objectif d’améliorer la transmission de l’impact des incertitudes statiques sur les réponses du simulateur dynamique d’écoulements sans augmenter le temps de calcul, grâce à des modèles approchés (proxy). Pour cela deux axes de recherche ont été menés : - L’implémentation de nouveaux proxys basés sur le Fast Marching, afin de modéliser la propagation d’un fluide dans un réservoir avec seulement quelques paramètres. Cela permet d’obtenir des courbes de réponse similaires à celles fournit par le simulateur d’écoulement pour un temps de calcul très court ; - La mise en place d’une procédure de minimisation mathématique afin de prédire les courbes de réponses du simulateur d’écoulement à partir d’un modèle analytique et des distances entre les modèles calculées avec les réponses des proxys. Les méthodes développées ont été appliquées sur deux cas d’études réels afin de les valider face aux données disponibles dans l’industrie. Les résultats ont montrés que les proxys que nous avons implémentés apportent de meilleures informations que les proxys disponibles bien que les nôtres soient toujours perfectibles. Nous avons aussi mis en évidence l’intérêt de notre procédure de minimisation pour mieux évaluer les incertitudes dynamiques à partir du moment où le proxy utilisé est suffisamment fiable / Although it is common to use many different numerical models for the static description of underground reservoirs and their associated uncertainties, for fluid flow uncertainties through these reservoirs only few dynamic simulations can be used due to performance reasons. The objective of this thesis’ work is to better transmit the impact of static uncertainties on flow simulator’s responses without increasing computation time, using approximated models (proxies). Research has been undertaken in 2 directions: - Implementation of new proxies based on Fast Marching in order to better approach fluid propagation behavior in a reservoir using only a few parameters. This allows to obtain response curves close to those provided by the flow simulator in a very short period of time - Set up a mathematical minimization’s procedure in order to predict flow simulator’s response curves using an analytical model and distances between proxy responses computed on every model. The methods developed during this PhD have been applied on two different real cases in order to validate them with industry data. Results have shown that our new proxy improve the quality of the information about fluid behavior compared to the available proxy even though ours can still be improved. We also highlight that our minimization procedure better assesses dynamic uncertainties if the proxy used is reliable enough

Ce travail propose des solutions innovantes en vue de l'exploitation des modèles d'ingénierie hétérogènes. Il prend pour exemple le domaine de la prospection pétrolière. Les stratégies de prospection sont élaborées à partir de représentations tridimensionnelles du sous-sol appelées modèles géologiques. Ceux-ci reposent sur un grand nombre de données hétérogènes générées au fur et à mesure de la conduite de l'exploration par des activités telles que la prospection séismique, les forages, l'interprétation des logs de puits. A fin d'optimisation, les utilisateurs finaux souhaitent, pouvoir retrouver et réutiliser à tout moment les données et les interprétations attachés aux différents modèles successivement générés. Les approches d' intégration des connaissances susceptibles d'être mises en œuvre pour résoudre ce défi, doivent être dissociées aussi bien des sources et des formats de données que des outils logiciels en constante évolution. Pour cela, nous proposons d'utiliser l'annotation sémantique, technique courante du Web sémantique permettant d'associer la connaissance à des ressources au moyen d' "étiquettes sémantiques". La sémantique ainsi explicitée est définie par un certain nombre d' ontologies de domaine, qui, selon la définition classique, correspondent à autant "de spécifications formelles de la conceptualisation" des domaines considérés. En vue d'intégrer les modèles d'ingénierie considérés, nous proposons une architecture, qui permet de relier des concepts appartenant respectivement à des ontologies locales et à une ontologie globale. Les utilisateurs peuvent ainsi avoir une vision globale, intégrée et partagée de chacun des domaines impliqués dans chaîne de modélisation géologique. Un prototype a été développé qui concerne la première étape de la chaîne de modélisation (interprétation séismique). Les expérimentations réalisées prouvent que, grâce à l'approche proposée, les experts peuvent, en utilisant le vocabulaire de leur domaine d'expertise, formuler des questions et obtenir des réponses appropriées.

La plupart des sites industriels abritent des équipements et des réservoirs sous pression. Pour des raisons diverses (suppression, impact mécanique, surchauffe ou autre), ils peuvent être endommagés et même éclater. Cette explosion peut engendrer de nombreux projectiles. Au cours de leur vol, ces derniers peuvent impacter d'autres équipements, tels que des réservoirs sous pression ou d'autres installations sensibles (poste de commande, par exemple). Si une des cibles impactées subit une ruine mécanique, elle peut exploser et générer une nouvelle série de projectiles. Ces projectiles menacent, à leur tour, d'autres installations et ainsi de suite. Ce type d'enchaînement accidentel catastrophique est connu sous le nom d'effet domino ou de suraccident. Dans ce document, l'effet domino pouvant se produire sur des sites industriels est analysé au travers des projections produites par l'accident initial. Une approche probabiliste globale est ainsi développée dans laquelle le calcul de la probabilité d'occurrence du phénomène requiert le passage par quatre étapes : - Analyse des termes sources : les projectiles générés par l'explosion d'un réservoir ont différentes caractéristiques, à savoir le nombre de projectiles, la forme, la masse, la vitesse de départ et les angles de départ. Toutes ces grandeurs sont modélisées par des variables aléatoires. A l'aide du principe du maximum d'entropie et des données existantes, des distributions probabilistes sont développées pour toutes ces variables. On se limite, cependant, au cas de l'explosion d'un réservoir cylindrique ou sphérique. - Analyse de l'impact ou analyse du mouvement : la trajectoire d'un projectile (ou fragment de structure), en fonction de ses caractéristiques de départ, est décrite par une combinaison des effets d'inertie, de gravitation et d'aérodynamique. Une approche simplifiée faisant l'hypothèse de constance des coefficients aérodynamiques permet d'identifier analytiquement la trajectoire du projectile tandis qu'une solution numérique est obtenue par une approche complète où toutes les valeurs de ces coefficients sont prises en compte.Les mouvements de translation et de rotation sont également étudiés. A l'aide de l'analyse complète et des conditions d'impact, la probabilité d'impact est déterminée. L'étude est restreinte à des projectiles en forme de fond de réservoir, fond oblong de réservoir et plaque. Les formes des cibles de l'étude sont restreintes au cas ellipsoïdal, cylindrique et cubique. - Analyse de l'état des cibles impactées : dans un premier temps, des modèles simplifiés d'impact sont utilisés afin d'étudier l'interaction mécanique entre les projectiles et les réservoirs impactés. Un modèle mécanique complet comprenant une loi de comportement élasto-plastique et un modèle de rupture est également proposé. Ce modèle est ensuite implémenté dans un code de calcul sans maillage de type SPH, i.e. Smoothed Particle Hydrodynamics. Afin d'estimer la probabilité de rupture des réservoirs impactés, les modèles simplifiés sont mis en œuvre, ce qui permet de réduire le coût de calcul. - Occurrence du sur-accident : selon l'état mécanique résiduel de la cible et son état physique (conditions thermodynamiques, niveau de remplissage, etc.), l'impact de projectiles peut conduire à la poursuite du pnénomène. Ce dernier point n'est pas traité dans le document présenté

Les propriétés pétrophysiques des réservoirs silicoclastiques sont influencées par de nombreux facteurs sédimentaires et diagénétiques. Les principaux phénomènes diagénétiques affectant les réservoirs sont généralement la cimentation de quartz et les compactions mécanique et chimique. Des réservoirs gréseux ayant des contextes géologiques différents ont été considérés dans cette étude ; les réservoirs carbonifères d'origine deltaïque-marine du bassin des Palmyrides-Sud en Syrie et les réservoirs cambro-ordoviciens du bassin de Sbaa en Algérie. Cette thèse consiste à établir l'histoire diagénétique, déterminer les contrôles sédimentaires et structuraux influençant l'évolution des phénomènes diagénétiques, caractériser l'habitus des cristaux authigènes de quartz formés autour les grains détritiques et aussi mettre en relation les différentes contrôles sédimentaires et diagénétiques sur les caractéristiques des pores et ainsi que sur la variation de la perméabilité. L'histoire diagénétique entre le réservoir du bassin de Sbaa se caractérise par une forte cimentation de quartz composée de trois phases Q1, Q2 /et Q3, par tapissage illitique et ainsi par une importante compaction chimique liée à certains faciès glaciaires et également une cimentation d'argiles principalement en illite mais surtout dans les champs d'Oued Zine et de Bou Hadid. A l'exception du champ de Hassi Ilatou, où une faible cimentation de quartz composée de Q1 a eu lieu. Alors que la diagenèse des réservoirs gréseux du bassin des Palmyrides-Sud est représentée par une faible cimentation de quartz composée d'une seule phase Q1, une absence de compaction chimique, ainsi qu'une cimentation d'argiles dominée par la chlorite et les kaolins. Les analyses microthermométriques des inclusions fluides dans les surcroissances de quartz mettent en évidence une silicification se déroulant principalement entre 100 et 160°C dans les deux bassins. D'après la reconstitution de l'histoire thermique de bassin, cet intervalle de température a été atteint entre le Viséen et la fin du Namurien dans le bassin de Sbaa et au Crétacé supérieur-Paléocène dans le bassin des Palmyrides-Sud. Les analyses isotopiques indiquent une eau originelle météorique et marine, progressivement réchauffée lors de l'enfouissement, et s'enrichissant au fur et à mesure en ¹⁸O dans les pores intergranulaires et des fluides évolués et chauds à l'origine des filonnets. L'habitus des cristaux authigènes de quartz et la forme de croissance montrent une relation avec les phases de ciment de quartz, son taux et la présence/absence de gaz. En effet, des cristaux à prisme court, tronqués par des facettes additionnelles, et des cristaux trapus caractérisent les grès cimentés par une seule phase de quartz authigène, et une fréquence importante des cristaux de quartz à multiples nucléas est constatée dans ces grès. Des cristaux à prisme développé et rarement des cristaux à prisme court caractérisent les grès contenant deux phases du ciment de quartz. Des cristaux peu développés et limités à quelques faces sont présents dans les grès cimentés par trois phases du ciment de quartz dans la paléozone à eau du réservoir dans le champ ODZ. Une forme de croissance en escalier est présente uniquement dans ces derniers grès. La présence des inclusions à hydrocarbures dans les surcroissances de quartz dans la partie supérieure du réservoir ordovicien du champ de Oued Zine indique que la mise en place des hydrocarbures dans le réservoir a été contemporaine à la cimentation de quartz à des températures 100-140°C en raison de la paléostructure anticlinale dans ce champ. Un deuxième épisode a eu lieu suite à la fracturation hercynienne à des températures comprises entre 117-185°C qui augmente vers le nord-ouest du bassin. La composition du gaz dans les inclusions monophasées (92 ± 5 mole %) est comparable à la composition actuelle du gaz dans le réservoir.

Les systèmes méandriformes construisent, à long terme, des architectures sédimentaires composées de corps poreux disséminés dans un volume de sédiments peu perméables. Ces corps sont susceptibles de contenir des ressources naturelles. Afin d'optimiser leur exploitation, il importe d'estimer la répartition et la connectivité de ces corps. À cet effet, Mines-Paristech développe Flumy, un logiciel de modélisation par processus de ces architectures . Cette thèse vise à améliorer la simulation de la migration dans Flumy en introduisant une influence.Pour cela, trois différents modèles ont été reliés entre eux et comparés. Le premier (modèle à pente constante) est celui actuellement en place dans Flumy, héritier des travaux de [Ikeda 1981] et dans lequel la pente de la surface du cours d'eau est constante. Le deuxième (modèle à pente variable) fut proposé par [Lopez 2003]. Il attribue à la surface libre la pente locale de la topographie. Enfin, un dernier modèle (modèle Saint Venant) a été élaboré en modifiant le précédent. L'écoulement moyen, initialement calculé sous une surface libre connue, est désormais obtenu par résolution des équations de Saint Venant sur un fond connu. Ces trois modèles ont été appliqués à la migration à long terme de méandres libres. Les modèles à pente constante et à pente variable ont également été utilisés pour reproduire les méandres confinés de deux rivières canadiennes, la Beaver et la Red Deer.Les résultats montrent un développement plus réaliste des méandres simulés avec le modèle à pente variable que de ceux simulés avec le modèle à pente constante. Cette amélioration est observable dans les méandres isolés (ralentissement de la vitesse d'extension) et dans les trains de méandres (confinement en ceintures de méandres). Le modèle à pente variable permet également une meilleure reproduction de la morphologie des méandres confinés. L'utilisation du modèle Saint Venant apporte les mêmes améliorations dans la simulation des méandres libres, quoique moins prononcées. Il permet surtout de construire une surface libre physique sur une grande variété de fonds, résolvant ainsi une limite identifiée dans le modèle à pente variable.

Les interfaces cerveau-machine visent à établir un lien de communication entre le cerveau et un système externe à ce dernier. Les électroencéphalogrammes (EEG), dans ce contexte, ont l’avantage d’être non invasifs. Par contre, l’information sensorielle qui se retrouve dans un signal EEG est beaucoup moins ciblée que dans un signal neuronal acquis par une méthode invasive. De plus, étant donné que le cortex auditif est situé dans des repliements du tissu cortical, les neurones qui déchargent, suite à un stimulus auditif, sont parallèles à la surface corticale sur laquelle les EEG sont enregistrés. Par conséquent, l’information auditive qui se retrouve dans le canal EEG situé vis-à-vis du cortex auditif est faible. L’objectif principal de ce projet de recherche consiste donc à étudier la répartition de l’information auditive dans l’ensemble des canaux EEG. Pour ce faire, nous utilisons deux approches. Dans la première, nous tenterons d’estimer l’activité corticale sous-jacente à partir des signaux EEG en utilisant un modèle de couplage bande fréquence. En effet, certaines bandes de fréquences sont des bons prédicteurs des décharges neuronales. Cependant, cette approche n’a pas été validée pour le système auditif, nous confronterons donc l’estimation obtenue à une autre estimation en ayant recours à un modèle spécialisé pour l’encodage du signal de parole faisant appel aux processus ponctuels. Ce modèle prend en compte les dynamiques intrasèques des neurones et également des propriétés spectrotemporelles du stimulus d’entrée. Dans la seconde approche, nous étudierons la possibilité de classifier 3 voyelles (a, i et u) en fonction du nombre de canaux EEG utilisés ainsi que leur répartition sur le cuir chevelu. Nous aurons recours, pour cela, à un réservoir de neurone à décharge récurrent activé en entrée par les données EEG. Les résultats démontrent que l’information auditive se retrouve en fait dans l’ensemble des canaux EEG et qu’elle n’est pas confinée à un nombre restreint d’électrodes. Il est également montré que lorsque l’on utilise les 64 électrodes que comporte l’EEG pour classifier les 3 voyelles, on obtient une classification de l’ordre de 80%, mais aussi qu’un nombre limité de 10 électrodes suffit pour obtenir une classification satisfaisante et, qu’en plus, la position de ces électrodes sur le cuir chevelu est peu importante.

Les matériaux composites constitués de fibres de carbone et de résine époxy présentent des propriétés spécifiques remarquables qui les destinent à une utilisation à grande échelle dans de nombreux domaines où des gains de masse sont recherchés, comme le transport. Les réservoirs d’hydrogène de type IV présents dans les véhicules automobiles en sont un exemple. Dans le cadre de la sécurité des biens et des personnes, le risque incendie doit alors être considéré : lorsqu’un réservoir d’hydrogène subit une agression thermique telle qu’un feu, son enveloppe en matériau composite est sujette à décomposition thermique qui, couplée aux transferts de chaleur et à l’endommagement dû au chargement mécanique, peut conduire à l’explosion de la structure. Pour prédire numériquement le comportement des réservoirs dans des situations de charges couplées, thermiques et mécaniques, un modèle construit sur la prise en compte des phénomènes ayant une conséquence majeure sur le comportement du matériau a été développé. Il intègre, dans un cadre thermodynamique, la microfissuration matricielle, la rupture des fibres et des interfaces fibres/matrice ainsi que le délaminage et représente les effets de la température sur les propriétés mécaniques. A cet endommagement mécanique est adjoint la décomposition thermique due aux températures élevées (>350°C). Elle induit des changements de structure du matériau due à la gazéification de la résine époxy, une modification des paramètres thermiques (ce qui a une influence sur les transferts de chaleur) et une perte des propriétés mécaniques.A l’échelle de l’éprouvette, des calculs sont menés afin de déterminer les paramètres des différents sous-modèles. Cela comporte les critères d’initiation et les lois d’évolution des endommagements, les paramètres réactionnels pour la décomposition thermique et les paramètres thermiques (masse volumique, capacité thermique et conductivité) pour chaque état de décomposition. Une méthode est proposée pour fixer les paramètres clés du couplage thermomécanique, à savoir l’influence de la décomposition thermique sur le comportement mécanique. Des calculs entièrement couplés sont également réalisés afin de déterminer le poids de chaque phénomène (température, décomposition thermique et endommagement mécanique) sur la rupture finale du matériau en condition d’exposition au feu.A l’échelle du réservoir hyperbare, des prédictions de pressions d’éclatement, en situation d’incendie et à température ambiante, sont menées. A température ambiante, le rôle de chaque processus d’endommagement dans la pression ultime est évalué afin de déterminer leur nécessité dans la modélisation. Dans des conditions d’incendie, le temps à éclatement est évalué lorsque le réservoir est soumis à différentes pressions internes. Le modèle est capable de prédire correctement la transition entre un mode d’éclatement à haute pression et un mode de fuite par le corps à plus faible pression, dû à la fusion du liner avant que le niveau de contraintes soit critique dans la coque composite et conduise à l’éclatement. Cette approche, mise en place à l’échelle du réservoir, permet donc d’établir la succession des événements menant à son éclatement (ordre et plis d’apparition des endommagements, champ de température dans l’épaisseur l’enveloppe composite, temps à explosion). / Composite materials made of carbon fibres and epoxy resin have remarkable specific properties that make them suitable for large-scale use in many areas where mass savings are required, such as transport. An example is the type IV hydrogen tank in motor vehicles. In the context of the safety of persons, the fire risk must then be considered: when a hydrogen tank undergoes thermal aggression such as a fire, its composite shell is subject to thermal decomposition which, coupled with heat transfers and damage due to mechanical loading, can lead to burst of the structure. To numerically predict the behaviour of tanks subjected to coupled load (thermal and mechanical), a model based on the phenomena having a major impact on the behaviour of the material has been developed. It involves, in a thermodynamic framework, matrix microcracking, fibre failure and fibre/matrix interfaces decohesion as well as delamination and represents the effects of temperature on mechanical properties. To this mechanical damage is added the thermal decomposition due to high temperatures (>350°C). It induces structural changes in the material due to the gasification of the epoxy resin, a change in thermal parameters (which has an influence on heat transfers) and a loss of mechanical properties.At the specimen scale, calculations are carried out to determine the parameters of the different sub-models. This includes initiation criteria and damage evolution laws, reaction parameters for thermal decomposition and thermal parameters (density, thermal capacity and conductivity) for each decomposition state. A method is proposed to determine the key parameters of thermomechanical coupling, namely the influence of thermal decomposition on mechanical behaviour. Fully coupled calculations are also performed to determine the weight of each phenomenon (temperature, thermal decomposition and mechanical damage) on the final failure of the material under fire exposure conditions.At the hyperbaric tank scale, burst pressure predictions, in a fire situation and at room temperature, are carried out. At room temperature, the role of each damage process in the ultimate pressure is evaluated to determine their relative weight in the modelling. Under fire conditions, the time to burst is evaluated when the tank is subjected to different internal pressures. The model is able to correctly predict the transition from a burst mode at high pressure to leak mode at lower pressure due to melting of the liner before the stress level is critical in the composite shell and leads to burst. This approach, implemented at the scale of the tank, therefore makes it possible to establish the sequence of events leading to burst (plies in which damage occurs, temperature field in the composite wall, burst time).

Trois magmas (haplotonalite synthétique, basalte et dacite de Santorin) ont été juxtaposés et déformés à l’aide d’une presse de type Paterson pour contraindre les conditions de mélange de magmas à haute pression (300 MPa). Ces trois magmas ont été utilisés en conditions sèches ou hydratées et dans une gamme de température comprise entre 600 et 1200°C, permettant d’obtenir plusieurs fractions cristallines et contrastes de viscosité. Les textures de mélange produites lors des expériences dépendent de la fraction cristalline et sont similaires aux textures rencontrées dans la nature. Les textures de mélange mécanique (mingling) produites sont la ségrégation de cristaux depuis leur magma source, la formation d’enclaves par détachement et de filaments par étirement de parcelles de magma. Le mélange chimique (mixing) est illustré par des zones d’interactions comportant une large variété de liquides intermédiaires et la cristallisation de nouvelles phases. Le mélange est produit en l’absence d’une charpente cristalline dans l’un ou l’autre des magmas, si le contraste de viscosité est faible (< 0,3 unité log) et en dessous d’un seuil de viscosité absolue, déterminé entre 107 et 108 Pa.s. Par ailleurs, l’eau joue un rôle important sur les propriétés rhéologiques des magmas, et donc sur leurs capacités de mélanges : elle abaisse le seuil de mélange de près de 200°C entre des magmas saturés en eau, et son exsolution (présence de bulles) entraîne une réduction significative des viscosités promouvant le mélange à des fractions cristallines plus faibles. Les conditions rhéologiques favorables aux mélanges entre le magma d’un réservoir et un magma plus chaud qui le recharge ont été déterminées selon la fraction de magma injecté. Cette fraction est de 0,5 minimum pour la majorité des réservoirs typiques du contexte d’arc. La comparaison avec des systèmes plutoniques et volcaniques exposant des figures de mélange montre que les conditions favorables aux mélanges sont atteintes avec des fractions inférieures. Ceci suggère que le fonctionnement d’un réservoir magmatique est influencé par ses styles et taux de recharge, ainsi que la quantité de volatils. / Three magmas (a synthetic haplotonalite, a basalt and a dacite from Santorini Volcano) have been deformed together at P=300 MPa in a Paterson type apparatus to constrain the mixing capacities between chemically contrasted magmas. Both dry and hydrous conditions were investigated at temperatures varying between 600 and 1200°C to vary magma crystallinity and visc osity contrasts. Mingling and mixing features were produced, which closely resemble those encountered in nature. Mingling features are xenocrysts, enclaves, and melt filaments. Diverse intermediate liquids and newly grown crystals characterize the mixing textures. Both mixing and mingling features were produced after crystal network breakage, at low viscosity contrasts (less than 0.3 log unit) and bulk viscosities lower than 107 – 108 Pa.s. The presence of water considerably affects the rheological properties of magmas and their mixing capacities: water-saturated magmas are able to mix at temperatures 200°C colder than their dry equivalent, and bubbles may influence the viscosity contrast leading to mixing at lower crystal fractions. Rheological conditions favourable to magma mixing between the resident magma in the reservoir and a hot replenishing one were determined according to the injected mafic magma fraction. The minimum fraction of mafic magma for mixing is 0.5 for typical arc magma chambers. The comparison with plutonic and volcanic rocks displaying mixing features shows that mixing may, however, occur at lower fractions. This suggests that reservoir dynamism may depend on the style and the frequency of the injection rate, and the amount of volatiles.

Dans le cadre de monitoring des puits et des réservoirs pétroliers, les études de thermométries se sont développées. Suite à l'émergence de nouvelles technologies d'acquisition avec l'apparition des fibres optiques depuis une vingtaine d'années, il est maintenant possible d'obtenir un profil de température continu à la fois sur la hauteur du puits et dans le temps. Dans ce travail, on étudie d'un point de vue thermomécanique deux modèles d'écoulements afin d'interpréter ces mesures de température. <br /><br />On couple d'abord un modèle réservoir axisymétrique avec un modèle puits pseudo 1D pour l'écoulement d'un fluide monophasique. La modélisation du réservoir intègre une équation d'énergie non-standard qui prend en compte les phénomènes de décompression du fluide (Joule-Thomson) ainsi que ceux liés à la friction dans la formation (dissipation visqueuse). Le puits est décrit par un modèle couplant classiquement les équations de Navier-Stokes compressibles avec un bilan d'énergie. Des conditions de transmission adéquates sont imposées à l'interface entre les deux domaines et ensuite dualisées par des multiplicateurs de Lagrange. On obtient ainsi une formulation variationnelle mixte non standard pour le problème couplé. Les flux sont discrétisés par des éléments finis de Raviart-Thomas. Enfin, on montre que les problèmes couplés continu et discret sont bien posés.<br /><br />On développe ensuite un modèle réservoir multi-phasique (huile, gaz et eau) et multi-composant. Compte tenu de la complexité du problème, on a choisi d'étendre un simulateur non-thermique existant GPRS (General Purpose Reservoir Simulator) en lui ajoutant une équation d'énergie adéquate et la thermodynamique correspondante. L'écoulement est régi par les équations de conservation de la masse pour chaque composant et une équation de conservation d'énergie exhaustive, couplées avec la loi de Darcy géenéralisée appliquée à chacune des phases. L'ensemble des propriétés thermodynamiques des fluides et les équilibres des phases sont calculés par l'équation d'état de Peng-Robinson. La discétisation en espace est faite à l'aide des volumes finis et le système non-linéaire obtenu est résolu par la méthode de Newton-Raphson.<br /> <br />Des tests numériques avec des données réelles validant les codes développés sont présentés.