D’importantes découvertes de gaz ont été faites récemment en Méditerranée orientale (www.nobleenergyinc.com), incitant les compagnies pétrolières à s’intéresser de plus près au bassin du Levant, considéré aujourd’hui comme une véritable province pétrolière. Par conséquent, une quantité considérable de données géophysiques a été produite et une série d'études académiques et industrielles ont été réalisées. La compréhension de l’architecture crustale et sédimentaire associée à celle de la thermicité actuelle et passée des marges de ce bassin, notamment la marge continentale du Liban, présente des enjeux industriels et scientifiques majeurs. Cette question a des implications majeures pour l'évolution tectonique, les prévisions des tremblements de terre ainsi que celle des systèmes pétroliers. Malgré les différents travaux géophysiques menés sur la Méditerranée orientale ces dernières années, la configuration crustale profonde du bassin du Levant, connu pour avoir été le siège d’un rifting à la fin du Paléozoïque et au début du Mésozoïque, reste imprécise. La transition d’une croûte continentale épaisse vers une croûte atténuée en mer (peut-être même une croûte océanique) a été invoquée, mais pas encore prouvée. Des approches géophysiques intégrées ainsi qu’un travail de modélisation ont été utilisés dans cette thèse pour étudier la structure profonde de la lithosphère sous la région Est de la Méditerranée.Une modélisation crustale 2D à l’échelle régionale (du delta du Nil au sud à la Turquie au nord, et du bassin Hérodote à l’ouest à la plaque arabe à l’est) a été effectuée dans le but d’étudier l’architecture de la croûte dans cette partie de la méditerranée orientale. L’algorithme utilisé est une méthode d’essai-erreur qui fournit l’épaisseur crustale et la profondeur de la limite lithosphère- asthénosphère (LAB) ainsi que la distribution de la densité crustale par l’intégration du flux de chaleur surfacique, l’anomalie gravimétrique à l’air libre, les données du géoïde et la topographie. La profondeur du Moho et l’épaisseur de la croûte ont été contraintes localement par des données de sismique réfraction là où elles sont disponibles. Les résultats montrent une croûte cristalline progressivement atténuée dans une direction EW. Dans le bassin du Levant, la croûte est interprétée comme continentale et composée de deux croûtes distinctes, une supérieure et une inférieure, contrairement au bassin Hérodote qui repose sur une croûte mince, probablement océanique.Une inversion 3D jointe des données de gravité, du géoïde et de la topographie appliquée sur la même région a confirmé les résultats de la modélisation crustale 2D. A total of 168 simulations ont été réalisées, parmi lesquelles, la simulation avec les erreurs les moins grandes sur les données correspond à l’inversion d’un modèle dans lequel la profondeur du Moho varie entre 23 et 26 km dans le bassin du Levant et devient plus profond dans le bassin Hérodote et aux larges des côtes africaines. La profondeur de la LAB est située entre 100 et 150 km dans le bassin du Levant et atteint plus de 180 km dans le bassin Hérodotes. L’interprétation de cinq lignes de sismique réflexion 2D PSTM à 14 s TWT couvrant la partie nord du bassin du Levant a révélé un total de 10 horizons, dont le plus profond pourrait être une interface croûte-manteau. L’interprétation des paquets sismiques, leurs surfaces de raccord ainsi que l’analyse des facies ont été contraints par les interprétations sismiques 2D publiées de la partie nord de l’offshore Libanais (Hawie et al., 2013b), dans lesquelles les connaissances stratigraphiques et sédimentologiques récentes de la marge libanaise ont été extrapolées jusqu’au bassin. Un total de huit paquets sédimentaires a été identifié dans le bassin aux âges variant du Jurassique Moyen au Quaternaire. / Significant gas discoveries have been made recently in the Eastern Mediterranean (www.nobleenergyinc.com), which turned the attention of oil companies towards the Levant Basin. This region is considered today as a typical hydrocarbon frontier province. Hence, a considerable amount of geophysical data has been produced and a series of academic and industry-based studies have been performed. Understanding the crustal and sedimentary architecture, the actual and past thermicity of this basin, in particular on the Lebanese continental margin, has major academic and economic interests. This has important implications on understanding tectonic evolution and earthquakes generation and on assessing petroleum systems. Despite numerous old and recent geophysical studies in this region, the deep crustal configuration of the Levant Basin, known to be the site of rifting in the Late Paleozoic and Early Mesozoic, remains enigmatic. The transition from a typical thick continental crust to thinner attenuated crust offshore (possibly even oceanic crust) has been invoked, but not yet proven. Integrated geophysical approaches and modeling techniques are used in this thesis to study the deep structure of the lithosphere underlying the easternmost Mediterranean region.A 2D modeling approach was accomplished at a regional scale (1000x1000 km2) extending from the Nile delta in the south, to Turkey in the north, from the Herodotus Basin in the west to the Arabian plate in the east. The algorithm used is a trial and error method that delivers the crustal thickness and the depth of the lithosphere-asthenosphere boundary (LAB) as well as the crustal density distribution by integrating top basement heat flow data, free-air gravity anomaly, Geoid and topography data. Moho depth and crustal thickness were locally constrained by refraction data where available. Three models are presented, two in EW direction (580 and 650 km long) and one in SN direction (570 km long). The models in EW sections show a progressively attenuated crystalline crust from E to W (35 to 8 km). The SN section presents a 12 km thick crust to the south, thinning to 9-7 km towards the Lebanese offshore and reaching 20 km in the north. The crystalline crust is best interpreted as a strongly thinned continental crust under the Levant Basin, represented by two distinct components, an upper and a lower crust. The Herodotus Basin, however, shows a very thin crystalline crust, likely oceanic, with a thickness between 6 and 10 km. The Moho under the Arabian plate is 35-40 km deep and becomes shallower towards the Mediterranean coast. Within the Levant Basin, the Moho appears to be situated between 20 and 23 km, reaching 26 km in the Herodotus Basin. While depth to LAB is around 110 km under the Arabian and the Eurasian plates, it is about 150 km under the Levant Basin and plunges finally to 180 km under the Herodotus Basin.A 3D joint inversion of gravity, geoid and topography data applied on the same region confirmed the results of the 2D modeling. A total of 168 of simulations were run, among which the simulation with the minimal data misfits corresponds to a model where the Moho depth varies between 23 and 26 km in the Levant Basin and becomes deeper in the Herodotus Basin and off the African coast. The LAB is 100 to 150 km deep in the Levant Basin and deepens to more than 180 km in the Herodotus Basin.

Ce mémoire présente le développement d'un dispositif expérimental permettant de mesurer deux composantes du gradient de pesanteur, Γzz et Γzx, ainsi que l'accélération de pesanteur. Ces grandeurs sont déterminées en mesurant l'accélération d'atomes froids de rubidium, en chute libre dans le vide, par interférométrie atomique. Pour la gradiométrie, la mesure différentielle est réalisée entre deux nuages atomiques spatialement distants. Pour la mesure de Γzz, l'utilisation de réseaux optiques mobiles permet d'obtenir deux nuages atomiques à partir d'une unique source atomique. Ce travail présente la mise en place du dispositif complet, avec notamment la réalisation de l'enceinte à vide, et des systèmes laser et micro-onde. Les lasers sont basés sur la technologie télécom doublée, permettant d'obtenir des modules compacts et robustes, afin d'envisager des applications embarquées. L'architecture laser originale permet de réaliser des expériences d'atomes froids combinant interférométrie atomique et réseaux optiques, en réduisant au minimum le nombre de composants. Le bruit du laser a été caractérisé, et il limite la sensibilité gravimétrique à 10-9 g en monocoup, la sensibilité différentielle à 10-10 g en monocoup, et la sensibilité gradiométrique à 38 E, en monocoup. / This thesis presents the development of the experimental setup allowing the measurement of two gravity gradient components, Γzz and Γzx, and the gravity acceleration. These quantities are resulted from the measuring of rubidium cold atoms acceleration, in free fall in vacuum, by atom interferometry. For gradiometry, the differential measurement is realized between two atomic clouds spatially separated. For the measurement of Γzz, the use of mobile optical lattices allows to get two atom clouds from a single atomic source. This work presents the setting up of the complete device, in particular with the built of the vaccum chamber, laser and micro-wave systems. Lasers are based on frequency-doubled telecom technology, which allows to obtain compact and robust systems, dedicated for onboard applications. The innovative laser architecture allows to combine atom interferometry and optical lattices, while minimizing the amount of components. The laser noise has been characterized, and limits the single shot gravimetric sensitivity to 10-9 g, the single shot differential sensitivity to 10-10 g, and the single shot gradiometric sensitivity to 38 E.

Au cours de ce travail, un dispositif de mesure couplé manométrique-calorimétrique a été développé pour déterminer simultanément les isothermes et les enthalpies d'adsorption dans des intervalles élevés de pression et température. Ce dispositif couplé comprend un calorimètre Calvet C80 (de setaram) associé à un dispositif manométrique développé in-situ. Il permet d'étudier les phénomènes d'adsorption de gaz sur des solides poreux pour des températures comprises entre 303.15K à 423.15K et des pressions comprises entre 0 et 2.5MPa. Dans un premier temps, une description complète de l'appareillage et des protocoles expérimentaux sont exposés. Dans un second temps, la fiabilité et la reproductibilité de ce dispositif a été testée en réalisant des mesures sur un échantillon poreux benchmark (Filtrasorb F400) à 318.15K. Les mesures ont été complétées, à différentes températures, avec celles obtenues par méthode gravimétrique afin d'établir des informations à haute température et haute pression. Ces dispositifs ont permis d'étudier l'adsorption de fluide supercritique (azote N2, méthane CH4, dioxyde de carbone CO2) dans des charbons actifs et des silices microporeuse ou mésoporeuse. L'adsorption du méthane dans une roche de type (TGR) a été également étudiée. Ces résultats expérimentaux sont utilisés pour l'étude des interactions fluide/substrat à prendre en compte dans la simulation moléculaire ou la théorie DFT.

Dans la première partie de ce travail de thèse, on a étudié les causes d'anomalie de pesanteur, et plusieurs causes de bruit afin d'en tirer des conclusions sur la faisabilité de certaines applications industrielles qui impliqueraient notamment l'utilisation de gravimètres ou de gradiomètres embarqués. On envisage en particulier la possibilité de constituer un système de prévention des collisions pour la navigation sous-marine, d'utiliser un gravimètre pour détecter des cavités enfouies ou encore d'observer l'anomalie de pesanteur créée par le passage d'une vague de tsunami d'une part, et d'autre part on cherche autant que possible à quantifier en spectre de puissance les bruits classiques rencontrés en gravimétrie embarquée, ainsi que le bruit gravitationnel causé par les vagues. Dans la seconde partie, on décrit la réalisation d'un gravimètre à ondes de matière, qui aura la particularité d'utiliser des atomes piégés au voisinage d'une puce en carbure de silicium. Le développement des gravimètres à ondes de matière est en effet extrêmement prometteur en terme d'exactitude de mesure du champ de pesanteur, mais le principe de réalisation utilisé jusqu'à maintenant implique que la sensibilité limite de l'instrument est proportionnelle à sa taille. D'un autre côté depuis une dizaine d'années des puces constituées de fils conducteurs déposés sur un substrat en silicium ont été développées pour le piégeage et le refroidissement d'atomes. L'utilisation d'une puce à atomes devra permettre de démontrer la possibilité de mesurer le champ de pesanteur avec une sensibilité indépendante de la taille de l'instrument, ce qui mènera à la réalisation d'un gravimètre à atomes froids compact, donc potentiellement utilisable dans un véhicule. Le défi de ce démonstrateur est d'effectuer pour la première fois la séparation spatiale cohérente d'un nuage d'atomes sur une puce atomique, à des fins de métrologie / In the first part of this work causes of gravity anomalies are studied, along with causes of noise. The feasibilities of a few industrial applications involving mobile gravity or gravity gradient meters were hence evaluated. These applications include in particular the realization of a collision avoidance device for underwater navigation, detection of underground void spaces and tsunami wave detection. Classical noises encountered in on-board gravity measurements are studied, as well as the less conventional gravity noise caused by ocean waves. The second part of the work is devoted to the implementation of a matter waves interferometry gravimeter. The particularity of the device is to use atoms trapped in the vicinity of a silicon carbide atom chip. The goal of the project is to reach for the measurement accuracies of current matter waves gravimeters with free falling atoms, with a principle that does not imply a proportionality between the size of the device and its sensitivity limit. We hope to demonstrate a proof of concept that may lead to a new generation of atomic gravimeters that are compact and therefore better suited for mobile uses. Coherent splitting of a non-condensed atom cloud for metrology purposes is probably the main challenge of the project

Nous présentons le développement d'un prototype de gravimètre atomique compact reposant sur l'interférométrie atomique avec des transitions Raman stimulées. Nous démontrons une amélioration importante de la compacité et de la simplicité de chaque élément du dispositif expérimental (tête de senseur, source laser, référence de fréquence micro-onde, système de filtrage des vibrations). Ce travail s'appuie sur l'utilisation d'une pyramide creuse comme miroir de rétro-réflexion, ce qui permet de réaliser toutes les fonctions d'un interféromètre atomique (piégeage et refroidissement des atomes, interféromètre, détection) avec un unique faisceau laser. Nous avons donc développé une tête de senseur très compacte, dont les fonctions clés ont toutes été simplifiées. La source laser met en jeu un unique laser émettant à 1560 nm pour interroger des atomes de Rubidium 87. Elle bénéficie de l'utilisation de composants optiques télécoms fibrés, qui ont déjà démontré leur performance et leur robustesse aux conditions environnementales. Tous les éléments du prototype sont assemblés pour permettre la mise en place de l'interféromètre. Ce type de gravimètre compact est très intéressant pour la gravimétrie de terrain. En parallèle, nous avons développé un système de réjection du bruit de vibration, basé sur l'électronique numérique. La contribution des vibrations sur la phase atomique est pré-compensée avant la fermeture de l'interféromètre, directement sur la phase optique des lasers. Ceci garantit que chaque point de mesure a une sensibilité maximum, malgré un bruit de d'accélération important. Ainsi, pour un gravimètre posé au sol en environnement urbain, nous avons démontré une sensibilité à l’accélération de l’ordre de à 1 seconde, qui atteint après 300 secondes d’intégration. Notre dispositif nous a finalement conduit à l’hybridation complète du gravimètre atomique avec un accéléromètre classique, conduisant à un accéléromètre exact très large bande [DC , 430 Hz]. Ce résultat est très prometteur, notamment pour la navigation inertielle. / We present the development of a compact atomic gravimeter, relying on atom interferometry using stimulated Raman transitions. We demonstrate a significant improvement in terms of compactness and simplicity for each element of the device (sensor head, laser source, micro-wave frequency reference, vibration rejection system). This work relies on the use of a hollow-pyramid in place of the usual retro-reflecting mirror. This component allows realizing each step of measurement (trapping and cooling of the atoms, interferometer, detection) with a single laser beam. We developed a very compact sensor head, for which we have simplified every key element. The laser source features only one single laser diode, emitting at 1560 nm to interrogate Rubidium 87 atoms. It benefits from the use of fibered optical telecom components, which have already demonstrated to be performing and robust to environmental conditions. All the elements of the accelerometer prototype are now gathered to allow for the operation of the interferometer. On-field gravimetry will greatly benefit from such compact absolute gravimeter. In parallel, we have developed a novel vibration rejection method, based on digital electronics. The atomic phase shift induced by vibrations is pre-compensated before the recombination of the wave-packets, directly on the optical phase on the Raman lasers. This ensures that each measurement point stays at maximum sensitivity, even in the presence of great acceleration noise. As a result, for a gravimeter operating directly on the ground in an urban environment, we have reached sensitivity to acceleration at a level of at 1 s, which improves down to after 300 s. Finally, our method lead us to fully hybridize the atom gravimeter with a classical accelerometer, which results in an accurate very large-band accelerometer [DC , 430 Hz]. This demonstration is very promising for applications in inertial navigation.

Le projet balance du watt propose de relier la définition du kilogramme à la constante de Planck h. La pesée de la masse impliquée nécessite une détermination de l'accélération de la pesanteur g avec une exactitude meilleure que 10-8. Cette thèse vise à réaliser cette détermination à l'aide d'un gravimètre atomique et d'un site gravimétrique dédié. Avec un gravimètre relatif caractérisé, une cartographie gravimétrique des deux massifs consacrés à l'expérience du LNE a été réalisée puis un modèle des variations de gravité a été développé. Il permet de déterminer la différence de gravité entre deux points dans un volume de 50 m3 au dessus des massifs avec une incertitude inférieure à 3 µGal (3 * 10-8 m.s-2). Entre les deux points centraux, l'incertitude est inférieure au µGal. La détermination absolue de g est réalisée avec un nouveau gravimètre absolu mobile à ondes de matière de 87Rb. Sa conception repose sur le travail débuté en 2002 avec la réalisation d'un prototype dont les limites ont été identifiées. Les différents éléments de ce nouveau gravimètre sont caractérisés dans cette thèse. Les premiers signaux ont été obtenus en 2009 puis l'instrument a été le premier gravimètre atomique à participer à une comparaison internationale : l'ICAG'09 au BIPM. La caractérisation de l'instrument a été poursuivie sur le site du LNE où la sensibilité atteint un plateau à 0,4 µGal après 100 min de mesure. Le budget d'incertitude obtenu de 5,4 µGal a été éprouvé lors d'une comparaison bilatérale avec un FG5 : l'écart de mesure obtenu est de (4,3 +- 6,2) µGal (k=1).

Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).

La méthode classique de traitement des données aérogravimétriques par filtrage fréquentiel ne prend pas en compte les propriétés statistiques des mesures, telles les corrélations, et ne permet pas une estimation aisée des erreurs sur le champ de gravité. Nous avons donc développé une nouvelle méthode de réduction par inversion régularisée au sens des moindres carrés, basée sur la transformation de l'équation différentielle du gravimètre en équation intégrale, qui permet cette estimation. Nous avons appliqué cette nouvelle méthode avec succès à la modélisation du champ de gravité dans le domaine des moyennes longueurs d'ondes sur la partie Sud-Ouest des Alpes avec une précision interne de 3 mGal. Les données utilisées ont été acquises avec un gravimètre « LaCoste & Romberg modèle Air/Sea » lors d'un levé gravimétrique aéroporté réalisé en 1998 dans le cadre du programme national de recherche « GéoFrance 3D ».

Au cours du dernier siècle, plusieurs séismes majeurs ont affecté l'Himalaya. Cependant, la taille maximale de ces événements et la probabilité d'occurrence de méga-séismes avec des magnitudes proches de 9 sont toujours matière à débat. L'étude de la segmentation de l'arc Himalayen est donc primordiale afin de comprendre les mécanismes qui contrôlent ces séismes ainsi que leur extension spatiale. La compréhension du cycle sismique en Himalaya est aujourd'hui essentiellement basée sur des études menées au Népal central, ce qui limite notre connaissance de son fonctionnement tri-dimensionnel. Ce travail de thèse permet d'apporter de nouvelles contraintes sur les variations latérales de la déformation crustale dans la zone comprise entre l'ouest Népal (~80°E) et l'est Bhoutan (~92°E). La première partie de cette thèse est consacrée à l'étude des variations latérales de la structure lithosphérique. Quatre campagnes gravimétriques ont été réalisées entre 2010 et 2012. Elles ont permis d'étendre le réseau gravimétrique Népalais jusqu'à l'ouest Népal et d'établir le premier réseau gravimétrique Bhoutanais. Ce nouveau jeu de données, combiné avec les données sismologiques disponibles, permet de contraindre des modèles thermo-mécaniques de la flexure de la plaque Indienne sous l'Himalaya. Les résultats suggèrent qu'il n'existe pas de variations latérales majeures du comportement mécanique de la lithosphère Indienne entre le centre et l'est Népal contrairement au Bhoutan où une rigidité flexurale plus faible est nécessaire pour expliquer les données.La deuxième partie de cette thèse est consacrée à l'étude des variations latérales des déformations récentes dans le prisme Himalayen. Les études paléosismologiques menées depuis 15 ans en Himalaya ont permis d'étendre le catalogue des événements majeurs sur le dernier millénaire. En intégrant à la fois les contraintes disponibles sur ces séismes majeurs et le chargement séculaire, nous étudions les variations spatiales et temporelles de l'état de contrainte du chevauchemlent Himalayen principal le long de l'arc. Nos calculs montrent que même dans la cas où un méga-séisme se serait produit en ~1100 en Himalaya, le chargement séculaire a quasiment compensé la chute de contrainte associée. Les modèles montrent aussi que le séisme du Shillong ne joue pas un rôle majeur sur les contraintes accumulées sur le chevauchement Himalayen depuis 1897 au niveau du Bhoutan. Enfin, nous présentons les premières contraintes sur la tectonique active au Bhoutan. L'étude morphotectonique réalisée au centre du Sud-Bhoutan a permis de montrer que la majeure partie de la déformation Holocène est accommodée au niveau du front Himalayen, comme au Népal. Nous montrons aussi des évidences de ruptures de surface avec des décalages verticaux de plusieurs mètres associés à deux séismes majeurs sur le dernier millénaire. Ces évidences remettent en cause l'interprétation d'un méga-séisme en ~1100.

Le développement de l’énergie géothermique a conduit à l’exploitation de ressources établies dans des contextes géologiques et géodynamiques très variés. L’exploration géophysique de ces réservoirs complexes nécessite l’utilisation de plusieurs méthodes d’imagerie complémentaire. Ce travail de thèse porte sur l’exploration d’une ressource géothermique située en contexte de socle fracturé dans le Massif Central français par magnétotellurique, tomographie de bruit ambiant et gravimétrie.La magnétotellurique est une méthode d’imagerie 3D résolvante qui est sensible à la présence d’eau et aux argiles d’altération hydrothermale mais limitée par sa couverture spatiale. La tomographie de bruit sismique présente une bonne résolution verticale mais ne résout pas les variations horizontales de vitesse. Cette méthode est sensible aux variations des propriétés mécaniques des roches et donc aux milieux fracturés. Enfin la gravimétrie apporte une contrainte sur les variations lithologiques et possède une bonne résolution latérale mais une faible résolution verticale.Nous présentons une méthode d’inversion conjointe des données sismiques et gravimétriques sous contrainte d’un modèle de résistivité obtenu par inversion magnétotellurique indépendante. L’inversion conjointe nécessite de définir des couplages entre modèles. Par absence de connaissance a priori de relations pétrophysiques, nous avons couplé les modèles de densité, de résistivité et de vitesse avec une loi qui contraint les paramètres à être corrélés en moyenne. Cette stratégie vise à faire ressortir des relations caractéristiques des objets géologiques de la ressource géothermique.Cette méthodologie d’inversion conjointe a été testée sur des modèles synthétiques. L’application aux données réelles acquises dans le Massif Central a permis de définir une zone en profondeur de forte corrélation interprétée comme la transition ductile fragile. La partie intermédiaire des modèles, plus homogène, permet de distinguer différentes unités géologiques séparées par une zone de faille. Enfin la partie superficielle se distingue par une forte hétérogénéité des paramètres résultants probablement de processus d’altération de surface. / The development of geothermal energy has led to the exploitation of resources established in varied geological and geodynamic contexts. Geophysical exploration of these complex reservoirs requires the use of several complementary imaging methods. This PhD thesis focuses on the exploration of a geothermal resource located within the fractured basement in the French Massif Central using magnetotelluric, ambient noise tomography and gravimetry. Magnetotelluric is a 3D imaging method with a good resolution power that is sensitive to the presence of water and hydrothermal weathering clays but is limited by its spatial coverage. Seismic noise tomography has a good vertical resolution but does not resolve well horizontal velocity variations. This method is sensitive to variations of the mechanical properties of rocks and thus to fractured media. Finally gravimetry brings constraint on the lithological variations and has a good lateral resolution but lacks vertical resolution.We present a method of joint inversion of seismic and gravimetric data under the constraint of a resistivity model obtained by independent magnetotelluric inversion. Joint inversion requires defining model couplings. By lack of prior knowledge of petrophysical relationships, we have coupled the density, resistivity and velocity models with a law that constraints the parameters to be correlated on average.This strategy aims to bring out the characteristic relationships of the geological objects of the geothermal resource. This joint inversion methodology has been tested on synthetic models. The application to the real data acquired in the Massif Central has made it possible to define a deep zone of high correlation interpreted as the fragile ductile transition. The intermediate part of the models, more homogeneous, allows to distinguish different geological units separated by a fault zone. Finally the superficial part is distinguished by strong heterogeneity of the parameters resulting probably from surface alteration process.