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.