Une brèche est un ensemble de blocs anguleux noyés dans un ciment de nature variable. Les brèches hydrothermales sont générées par un processus de fracturation, de dissolution des fragments, ainsi que des changements de composition causés par des eaux souterraines sous pression à haute température. La nature de la majorité des processus impliqués dans la formation des brèches hydrothermales est bien comprise d'un point de vue géochimique et plusieurs modèles basés sur cette perspective existent. Par contre, il n'existe pas de modèles approchant ces processus d'un point de vue géométrique. Dans ce mémoire, nous proposons un modèle basé sur les automates cellulaires, capable de simuler les principaux processus qui interviennent dans la formation des brèches. Un automate cellulaire est un modèle discret qui consiste en une grille de cellules pouvant chacune prendre à un instant donné un nombre fini d'états. Le temps est également discret et l'état d'une cellule au temps t est fonction de l'état au temps t -1 d'un nombre fini de cellules appelé son voisinage. À chaque nouvelle unité de temps, les mêmes règles sont appliquées pour toutes les cellules de la grille, produisant une nouvelle génération de cellules dépendant entièrement de la génération précédente. Cette approche est compatible avec l'aspect discret de la dissolution des minéraux et permet l'étude de l'évolution géométrique de fragments de roche virtuelle. Plus spécifiquement, on veut mesurer la complexité morphologique des fragments par leur dimension fractale de bordure, une méthode de mesure utilisée sur des échantillons réels et permettant de valider notre modèle avec des données analogiques. Un simulateur a été conçu pour mettre en oeuvre un tel modèle. Celui-ci est codé en Java et l'interface graphique est en HTML. Des expériences sur le simulateur ont mis en évidence deux régimes de dissolution: l'un limité par la diffusion (Diffusion Limited Regime -DLR), l'autre cinétique. Le premier régime dépend de la surface exposée et on y observe l'arrondissement et le lissage progressif des fragments. Le second régime est indépendant de la surface et on observe la formation de cavités dendritiques et une augmentation progressive de la complexité morphologique. D'un point de vue géochimique, le régime DLR est dit «contrôlé par la surface» alors que le régime cinétique est dit «contrôlé par le transport». Les extensions possibles au modèle sont variées et nombreuses. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Brèche hydrothermale, Automate cellulaire, Modélisation, Dissolution, Dimension fractale.

La démarche d’intégration des hétérogénéités dans les modèles réservoirs en est à ses prémices dans le domaine du stockage géologique de CO2. C’est dans ce contexte que s’inscrivent ces travaux de thèse. Un protocole d’analyse depuis l’étude de terrain jusqu’aux simulations réservoirs a été établi. La caractérisation du Minjur Sandstone (formation Triasique d’Arabie Centrale) met en avant le caractère crucial de la connectivité des corps dans l’architecture du réservoir, notamment en liant génétiquement leur nature, leur connectivité et leur position dans la séquence de dépôt. S’appuyant sur la connaissance de cette formation, un modèle conceptuel est construit, puis reproduit stochastiquement par un algorithme permettant l’élaboration de modèles conditionnés par une histoire sédimentaire. Le protocole prévoit la création de 50 scénarios illustrant divers degrés de connectivité ; chaque scénario étant composé de deux modèles de même architecture mais à remplissage sédimentaire différent. Cette approche permet d’appréhender (a) l’impact de la connectivité et (b) des hétérogénéités sédimentaires sur les performances réservoirs. L’estimation de capacité par l’approche statique des volumes disponibles estime une capacité moyenne d’environ 13Mt (aquifère semi-infini de 25 km par 25 km et 60m d’épaisseur à 1000 m de profondeur). Les hétérogénéités internes (sédiments argileux appelés oxbow lakes) entraînent une différence de capacité de 30%. Les simulations dynamiques confirment ces résultats et révèle une variabilité de capacité de 23% liée la connectivité des corps. De plus les hétérogénéités réduisent la migration verticale du gaz ce qui peut augmenter l’intégrité du stockage. / In the CO2 storage context, heterogeneity has only been rarely considered in reservoir models to date. To address this key issue, the project aims at developing a workflow that manages the heterogeneity from the field observations up to the reservoir simulation. The characterisation of the Minjur Sandstone (a Triassic formation from Central Saudi Arabia) shows the crucial role of connectivity in the reservoir architecture, and the genetic link between the nature, location and connectivity of the sedimentary bodies in the sequence. Stemming from this study, a conceptual model was established and stochastically reproduced through an algorithm simulating models conditioned to a sedimentary history. Fifty scenarios were simulated, representing various connectivity degrees. Each of these scenarios is composed of two models, identical by their architecture but different in their internal sedimentary fill. This approach allows the study of the impact of the (a) reservoir bodies’ connectivity and (b) their internal sedimentary heterogeneity on the reservoir’s performances. The capacity estimates using a static calculation based on the available pore volumes reveals a mean capacity of 13 Mt (for a 25 x 25 km x 60 m semi finite aquifer at 1000m deep). The sedimentary heterogeneity (shaly deposits called oxbow lakes) is responsible for a 30% difference of capacity. The flow simulations confirm these results and show that the connectivity of the reservoir bodies creates a 23% capacity variation. Moreover, the heterogeneities tend to reduce the amount of CO2 able to reach the uppermost reservoir which may enhance the storage integrity.

Résumé indisponible

Cette thèse s'articule autour de la reconstruction géométrique de formes. Plusieurs méthodes, basées sur les diagrammes de Voronoï, sont proposées pour la reconstruction automatique d'objets naturels. L'application principale est la modélisation et l'imagerie géologique. Une première méthode permet la reconstruction de volumes et surfaces géologiqu- es à partir de données incomplètes et hétérogènes : données ponctuelles sur des affleurements, portions de contours cartographiques, sondages, coupes incomplètes ou interprétées, modèles numériques de terrains... L'idée majeure de la méthode consiste à assembler les objets différents selon leurs proximités, en utilisant le diagramme de Voronoï de ces objets. Les diagrammes de Voronoï sont des structures géométriques permettant de partitionner l'espace en régions d'influence. En pratique toutes les données sont discrétisées en un ensemble de points colorés, les couleurs représentant ici les caractéristiques géologiques ou géophysiques des données, que nous souhaitons imager. La partition "colorée" de ces points nous donne une première solution topologique au problème de reconstruction. Elle nous fournit en outre, une représentation du bord de l'objet géologique et de son intérieur. L'utilisation de courbes et de surfaces déformables sous contraintes (tension, courbure et respect de la topologie initiale) permet ensuite d'obtenir des interfaces plus lisses et plus conformes. Une étape particulière permet de prendre en compte des surfaces de discontinui- té comme les failles. Afin de représenter un objet S, non plus par des éléments discrets (polyèdres de Voronoi), mais par les valeurs positives d'une fonction continue, nous avons introduit une nouvelle méthode. L'objectif de la méthode est de définir une fonction interpolante s telle que l'ensemble des zéros de s passe exactement par les données de départ et soit une approximation cohérente et lisse de S par ailleurs. Dans un premier temps nous définissons, une fonction caractéristique locale en chaque donnée (point, contour...) et l'objet volumique final résulte alors d'une interpolation de ces fonctions.

Cette thèse doctorale aborde divers aspects méthodologiques de l'analyse de levés électromagnétiques aéroportés en domaine temporel (TDEM) pour une interprétation détaillée à finalités géologique et hydrogéologique. Ce travail s'est appuyé sur un levé réalisé dans la région de Courtenay (Nord-Est de la région Centre) caractérisée par un plateau de craie karstifié (karst des Trois Fontaines) recouvert par des argiles d'altération et des alluvions. Tout d'abord, une méthode de filtrage des données TDEM utilisant la Décomposition en Valeurs Singulières (SVD) a été développée. L'adaptation rigoureuse de cette technique aux mesures TDEM a permis de séparer avec succès les bruits, qui ont pu être cartographiés, et le " signal géologique ", diminuant grandement le temps nécessaire à leur traitement. De plus, la méthode s'est avérée efficace pour obtenir, rapidement, des informations géologiques préliminaires sur la zone. Ensuite, une analyse croisée entre le modèle de résistivité obtenu en inversant les données filtrées et les forages disponibles a été effectuée. Celle-ci a mené à une amélioration de la connaissance géologique et hydrogéologique de la zone. Une figure d'ondulation, séparant deux dépôts de craie, et le réseau de failles en subsurface ont pu être imagés, apportant un cadre géologique au karst des Trois Fontaines. Enfin, une nouvelle méthode combinant l'information aux forages et les pentes issues du modèle de résistivité EM a permis d'obtenir un modèle d‟une précision inégalée du toit de la craie. L'ensemble de ces travaux fournit un cadre solide pour de futures études géo-environnementales utilisant des données TDEM aéroportées, et ce, même en zone anthropisée.

Cette étude s'inscrit dans un contexte de demande croissante en ressources en eau, conduisant à mieux contraindre et estimer la distribution des aquifères karstifiés qui représentent une ressource importante, notamment sur le pourtour méditerranéen. Le but de ce travail est d'élaborer une méthode prédictive de la distribution de la karstification et de l'épaisseur de massifs carbonatés karstifiés, notamment pour les karsts profonds. La méthodologie est développée à partir de l'étude des karsts de la plateforme carbonatée téthysienne du Languedoc montpelliérain; région qui a subi une histoire tectonique complexe et polyphasée, entrainant des phases de karstification successives. Elle est basée sur la compréhension des relations entre eustatisme, géodynamique et niveau de base, enregistrées par les paléokarsts de la région. Le remplissage sédimentaire des paléokarsts permet d'identifier les épisodes de karstification et de les dater ; c'est un marqueur des variations -que l'on cherche à quantifier- du niveau de base local responsable de la formation et du comblement des karsts. On analyse 3 paléokarsts d'extension régionale dont on peut étudier : le remplissage karstique, la distribution le long de profils amont-aval et l'extension verticale indicatrice de l'amplitude de l'abaissement du niveau de base dont il est issu. On en déduit les causes géodynamiques à l'origine de leur formation: - Le " paléokarst des bauxites " à remplissage caractéristique, résulte de l'évolution en 2 phases de " l'Isthme Durancien " : altération bauxitique de la couverture néocomienne sur le sud du Massif Central, liée à une inversion au Barrémien, suivi d'une surrection ≥ 1600m à l'Aptien-Albien, causant l'érosion des altérites en amont et le dépôt des bauxites détritiques en aval. - Le " paléokarst à remplissage marin paléocène inférieur " indique des variations du niveau de base local ≥350m dans la zone amont sud-cévenole. On propose que ces variations de niveau de base d'amplitude très supérieures aux variations eustatiques et dans une zone faiblement déformée, résultent d'assèchements et d'ennoiements d'un bassin endoréique à seuil. - Le " paléokarst à remplissage détritique d'origine cévenole ", enregistre la chute du niveau de base au Serravalien-Tortonien, qui se traduit par l'incision de canyons dans l'arrière-pays de la marge du Golfe du Lion. On détermine une surrection de plus de 400 mètres à l'amont, près de 250m dans les blocs faillés intermédiaires et quelques dizaines de mètres dans la zone avale. En parallèle, on établit la géométrie et la distribution actuelles des massifs carbonatés dans la zone d'étude (y compris sous couverture), en construisant un modèle géologique 3D, des Cévennes à la Méditerranée. On intègre les données de sismique réflexion, de forages et les données gravimétriques. La restauration tectonique séquentielle d'une coupe Nord-Sud extraite du modèle géologique permet d'appréhender la position des paléo-niveaux de base par rapport aux massifs carbonatés, pour chaque époque de karstification. On déduit l'organisation 3D des paléo-écoulements à travers la géométrie complexe des massifs carbonatés languedociens, notamment pendant l'évènement messinien. En terme d'évaluation de la ressource en eau, on déduit que: - Les surrections de l'arrière-pays (Crétacé " moyen ", Miocène) produisent des karsts de grande amplitude, mais ne constituent pas de réservoirs significatifs. - La karstification à l'origine de réservoirs exploitables importants et profonds est créée par l'assèchement de bassins endoréiques, qui ont affecté la région au Paléocène inférieur et au Messinien . - Ces deux évènements sont à l'origine de deux principaux réservoirs karstiques profonds qui ont une signature gravimétrique : (1) En amont du chevauchement de Montpellier, entre le Pic Saint Loup, la faille des Matelles et le bassin de l'Hérault, le karst profond se développe sur une épaisseur atteignant 900m, et est attribuée à l'événement du Paléocène. (2) Au Sud du chevauchement de Montpellier les karsts côtiers s'étendent jusqu'à plus de 10 km au large, scellés par une épaisse couverture Miocène et Plio-Quaternaire, et résultent de l'événement Messinien.

La vallée de la Komadougou Yobé constitue un site privilégié de recharge de l'aquifère quaternaire du bassin du lac Tchad. Depuis les années 1980, le développement des cultures de rente (p. ex. poivron) a conduit à un doublement des surfaces irriguées (1995-2005) dans sa partie aval. Afin de quantifier la recharge supplémentaire induite, une modélisation des flux hydriques verticaux dans la zone non-saturée (0-8 m) a été engagée. La première étape, objet de ce travail de thèse, consiste à décrire l'organisation spatiale et les propriétés hydrodynamiques des corps sédimentaires constitutifs de l'alluvium.Des forages à travers les formations superficielles mettent en évidence des alternances sable-sable argileux dans la totalité de la zone non-saturée, conséquence des migrations du lit mineur de la Komadougou Yobé. Le calcul de la corrélation entre les forages indique une forte analogie avec les hétérogénéités identifiées en surface à partir de données de télédétection (Landsat 7 binarisée). Ces données sont donc utilisées pour l'apprentissage de statistiques multipoints représentatives des hétérogénéités et un modèle géologique 3D de l'alluvium est généré via l'algorithme snesim.Des suivis neutroniques de l'humidité le long de plusieurs profils caractéristiques des différentes unités sédimentaires associées à ce modèle sont réalisés pour des conditions de flux contrôlées en surface. A partir de simulations numériques 1D, des jeux de paramètres hydrodynamiques permettant de reproduire les humidités mesurées sont déterminés par une approche de type Monte-Carlo. Des densités de probabilité intégrant l'incertitude sur les mesures sont obtenues pour les paramètres de Mualem - van Genuchten décrivant les courbes de rétention et de conductivité hydraulique des sédiments.Une procédure 1D-distribuée est utilisée pour simuler les écoulements non-saturés verticaux au sein de plusieurs réalisations du modèle géologique et pour différents jeux de paramètres hydrodynamiques probables. La recharge diffuse calculée se montre particulièrement sensible au paramètre de pression d'entrée d'air attribué aux dépôts superficiels, siège des principales interactions sol-plante-atmosphère, ainsi qu'aux contrastes verticaux de conductivité hydraulique. / The downstream part of the Komadugu Yobe River is an important recharge area for the Lake Chad Quaternary aquifer. Since the 1980s, cash crop development (e.g. sweet pepper) has led to the doubling (1995-2005) of irrigated surfaces in the vicinity of the river. A modeling approach of vertical water fluxes through the vadose zone (0-8 m) was designed to quantify the related increase in groundwater recharge. The first step, which is the main topic of this PhD thesis, consisted in describing both spatial arrangement and hydrodynamic properties of the sedimentary bodies that make up the alluvium.Boreholes in surficial deposits highlighted sandy to clayey alternations within the whole unsaturated zone; this was interpreted as the result of frequent migrations of the River channel. Spatial correlation between bore logs showed strong similarities with heterogeneities depicted on ground by means of remote sensing data (binarized Landsat 7 image). This data were therefore used to train multiple-point statistics representative of heterogeneities, and a 3D geological model was generated through the snesim algorithm.For each representative sedimentary unit, soil moisture under controlled hydraulic surface conditions was monitored by vertical neutron probe soundings. Using 1D numerical simulations, different data sets of hydrodynamic properties that reproduced moisture measurements were determined by a Monte-Carlo approach. Probability density functions including measurement uncertainties were deduced for the Mualem - van Genuchten parameters which describe both retention and hydraulic conductivity curves.A 1D-distributed procedure was applied for modeling vertical flows in the unsaturated zone within several geological model realizations with different probable sets of hydrodynamics parameters. The simulated diffuse recharge was shown to be particularly sensitive to two main parameters: air-entry pressure linked to superficial deposits, where soil-plant-atmosphere interactions do occur, and vertical hydraulic conductivity contrasts within the alluvium.

Cette thèse a pour but de contribuer à l'étude des réservoirs carbonatés karstiques par une approche intégrant modélisation géologique numérique 3D et hydrogéologie. Une méthodologie originale a été ainsi développée pour être appliquée à la caractérisation des ressources en eau souterraine de l'Unité du Beausset (Bouches-du-Rhône et Var, SE France). Cette méthodologie lie l'architecture géologique 3D du réservoir à ses propriétés hydrodynamiques extraites de son étude hydrogéologique. Cette zone est un synclinal de roches sédimentaires mésozoïques d'une épaisseur pouvant atteindre 3000m. Ces roches affleurent sur environ 650 km², et ont subi de nombreuses phases de fracturation et karstification. A partir de la synthèse des données existantes géologiques, hydrogéologiques et karstologiques, sept intervalles lithostratigraphiques ont été définis. Ces intervalles ont également été caractérisés hydrodynamiquement par l'analyse des séries temporelles de hauteur, température et conductivité électrique enregistrées en continu en sept sources représentatives. Cette caractérisation a été complétée par des campagnes exhaustives d'hydrochimie et de mesures de débit. L'ensemble des données et des résultats obtenus ont permis de définir où sont les unités de stockages et d'écoulement de l'eau souterraine. Ces dernières ont été intégrées de manière qualitative en 3D à l'aide d'un modèle surfacique de l'architecture géologique de l'Unité du Beausset. Cette intégration a permis de définir les meilleures zones de captage potentiel pour de nouvelles ressources en eau souterraine. / The aim of this PhD thesis is to provide an original approach for the characterization of karstic carbonate reservoirs that integrates 3D geological modeling and hydrodynamics. The 3D geological model stands for the structural compartmentalization of the reservoir and the hydrodynamics provide dynamic properties for the understanding of the fluid flow. This methodology is applied to the characterization of the groundwater resources of the Beausset Unit (SE France). The study area is a syncline composed of nearly 3000 m of Mesozoic sedimentary rocks, which age from late Triassic to Upper Cretaceous. These rocks, which are mainly carbonates, crop out over an area of 650 km² and have known several phases of fracturation and karstification. This study is divided in four parts: (1) synthesis of existing data, (2) 3D geological modeling, (3) study of hydrodynamics, and (4) integration of static and dynamic data in order to discuss the characterization of the groundwater resources in the Beausset Unit. The existing data on stratigraphy, lithology, petrophysical properties, karstology and hydrogeology allow the definition of seven lithostratigraphic intervals that have been chosen for modeling. These intervals represent the vertical compartmentalization of the reservoir: each interval has a unique combination of karst and rocks properties. The integration of all this data end results leads to the recognition of the flow and reservoir units in each lithostratigraphic interval. A qualitative link has been made between the 3D architecture of the basin and the reservoir properties; it allows the characterization of the best potential targets for new groundwater tapping zones.

Cette thèse doctorale aborde divers aspects méthodologiques de l’analyse de levés électromagnétiques aéroportés en domaine temporel (TDEM) pour une interprétation détaillée à finalités géologique et hydrogéologique. Ce travail s’est appuyé sur un levé réalisé dans la région de Courtenay (Nord-Est de la région Centre) caractérisée par un plateau de craie karstifié (karst des Trois Fontaines) recouvert par des argiles d’altération et des alluvions. Tout d’abord, une méthode de filtrage des données TDEM utilisant la Décomposition en Valeurs Singulières (SVD) a été développée. L’adaptation rigoureuse de cette technique aux mesures TDEM a permis de séparer avec succès les bruits, qui ont pu être cartographiés, et le « signal géologique », diminuant grandement le temps nécessaire à leur traitement. De plus, la méthode s’est avérée efficace pour obtenir, rapidement, des informations géologiques préliminaires sur la zone. Ensuite, une analyse croisée entre le modèle de résistivité obtenu en inversant les données filtrées et les forages disponibles a été effectuée. Celle-ci a mené à une amélioration de la connaissance géologique et hydrogéologique de la zone. Une figure d’ondulation, séparant deux dépôts de craie, et le réseau de failles en subsurface ont pu être imagés, apportant un cadre géologique au karst des Trois Fontaines. Enfin, une nouvelle méthode combinant l’information aux forages et les pentes issues du modèle de résistivité EM a permis d’obtenir un modèle d‟une précision inégalée du toit de la craie. L’ensemble de ces travaux fournit un cadre solide pour de futures études géo-environnementales utilisant des données TDEM aéroportées, et ce, même en zone anthropisée. / This PhD thesis addresses various methodological aspects of the analysis of airborne Time Domain ElectroMagnetic (TDEM) surveys for a detailed interpretation in geological and hydrogeological purposes. This work was based on a survey conducted in the region of Courtenay (north-east of the Région Centre, France) characterized by a plateau of karstified chalk (karst des Trois Fontaines) covered by weathering clays and alluvium. First, a TDEM data filtering method using the Singular Value Decomposition (SVD) was developed. The rigorous adaptation of this technique to TDEM data has successfully separated the noise, which was mapped, and the “geological signal”, greatly reducing the time required for processing. Furthermore, the method has proved to be effective to obtain quick preliminary geological information on the area. Then, a cross analysis between the resistivity model obtained by inverting the filtered data and the available boreholes was conducted. This has led to the improvement of the geological and hydrogeological knowledge of the area. An undulating feature, separating two chalk deposits, and a fault network were imaged in subsurface, providing a geological framework for the Trois Fontaines karst. Finally, a new 3D modelling method combining the information at boreholes and the slopes derived from the EM resistivity model yielded an accurate model of the top of the chalk. All of this work provides a solid framework for future geo-environmental studies using airborne TDEM data, even in anthropized area.

Cette étude s'inscrit dans un contexte de demande croissante en ressources en eau, conduisant à mieux contraindre et estimer la distribution des aquifères karstifiés qui représentent une ressource importante, notamment sur le pourtour Méditerranéen. Le but de ce travail est d'élaborer une méthode prédictive de la distribution de la karstification et de l'épaisseur de massifs carbonatés karstifiés, notamment pour les karsts profonds. La méthodologie est développée à partir de l'étude des karsts de la plateforme carbonatée téthysienne du Languedoc montpelliérain; région qui a subi une histoire tectonique complexe et polyphasée, entrainant des phases de karstification successives. Elle est basée sur la compréhension des relations entre eustatisme, géodynamique et niveau de base, enregistrées par les paléokarsts de la région. Le remplissage sédimentaire des paléokarsts permet d'identifier les épisodes de karstification et de les dater ; c'est un marqueur des variations -que l'on cherche à quantifier- du niveau de base local responsable de la formation et du comblement des karsts. On analyse 3 paléokarsts d'extension régionale dont on peut étudier : le remplissage karstique, la distribution le long de profils amont-aval et l'extension verticale indicatrice de l'amplitude de l'abaissement du niveau de base dont il est issu. On en déduit les causes géodynamiques à l'origine de leur formation: - Le « Paléokarst des bauxites » à remplissage caractéristique, résulte de l'évolution en 2 phases de « l'isthme Durancien » : altération bauxitique de la couverture Néocomienne sur le sud du Massif Central, liée à une inversion au Barrémien, suivi d'une surrection ≥ 1600m à l'Aptien-Albien, causant l'érosion des altérites en amont et le dépôt des bauxites détritiques en aval..- Le « Paléokarst à remplissage marin Paléocène inférieur » indiquent des variations du niveau de base local ≥350m dans la zone amont sud-cévenole. On propose que ces variations de niveau de base d'amplitude très supérieures aux variations eustatiques et dans une zone faiblement déformée, résultent d'assèchements et d'ennoiements d'un bassin endoréique à seuil.- Le « paléokarst à remplissage détritique d'origine cévenole », enregistre la chute du niveau de base au Serravalien-Tortonien, qui se traduit par l'incision de canyons dans l'arrière-pays de la marge du Golfe du Lion. On détermine une surrection de plus de 400 mètres à l'amont, près de 250m dans les blocs faillés intermédiaires et quelques dizaines de mètres dans la zone aval. En parallèle, on établit la géométrie et la distribution actuelles des massifs carbonatés dans la zone d'étude (y compris sous couverture), en construisant un modèle géologique 3D, des Cévennes à la Méditerranée. On intègre les données de sismique réflexion, de forages et les données gravimétriques. La restauration tectonique séquentielle d'une coupe Nord-Sud extraite du modèle géologique permet d'appréhender la position des paléo-niveaux de base par rapport aux massifs carbonatés, pour chaque époque de karstification. On déduit l'organisation 3D des paléo-écoulements à travers la géométrie complexe des massifs carbonatés languedociens, notamment pendant l'évènement Messinien. En terme d'évaluation de la ressource en eau, on déduit que: - Les surrections de l'arrière-pays (Crétacé « moyen », Miocène) produisent des karsts de grande amplitude, mais ne constituent pas de réservoirs significatifs.- La karstification à l'origine de réservoirs exploitables importants et profonds est créée par l'assèchement de bassins endoréiques, qui ont affecté la région au Paléocène inférieur et au Messinien. / The present study addresses the question of increasing need in water supply. One line of research consists in investigating karst aquifers that represent significant water supplies for the areas around the Mediterranean. This work aims at developing a methodology of predictive evaluation of karst distribution and karst thickness within the carbonate massifs. It is based on the study of the Tethysian carbonate platform from the area of Montpellier, South of France. This area underwent a polyphased, complex, tectonic evolution, involving successive karst forming events. The approach is based on the relationships between eustacy and geodynamics, on one hand, and base-level variation recorded within palaeokarsts, on the other hand. Sedimentary fill of palaeokarsts allows to identify distinct karstification phases and sometimes to date them. It is a marker of base-level variations, which we attempt to quantify, that is responsible for karst formation and karst filling and sealing. For three palaeokarsts of regional extent, we analyze : karst sedimentary-fill, the distribution along upstream-downstream profile, the vertical extension of palaeokarsts, indicative of the base-level lowering.- The « Bauxite palaeokarst » results from the two-stage evolution of the « Isthme Durancien » : i) Bauxitic weathering of the Neocomian cover over the southern Massif Central , related to an inversion during Barremian ; ii) ≥1600 m uplift and exhumation during Aptian-Albian, inducing erosion of alterites, upstream, and deposition of detrital bauxites, downstream.- The « Early Paleocene palaeokarsts with marine sedimentary-fill » indicate base level variations of more than 350m in the upstream area. It is suggested here that such base-level variations, exceeding eustacy and recorded in tabular undeformed series, result from desiccation-flooding events in an endoreic silled basin.- The « Miocene palaeokarst with Cevennes-sourced sediment-fill » records a Serravalian-Tortonian base-level drop, depicted by canyon incision in the hinterland of the Gulf of Lion Margin. Uplift values of ≥ 400m upstream, about 250m in an intermediate area and 10's of m downstream, are measured. The present 3D geometry and distribution of carbonate massifs (exposed and covered) are modeled in the study area, from the Cevennes to the Mediterranean. The 3D geological model integrates seismic reflection, borehole and gravimetric data. Tectonic restoration of a N-S section, extracted from the 3D model, allows to reconstruct paleo-base-levels with respect with the carbonate massifs, at different stages of deformation and successive phases of karstification. The 3D paleo-flows through carbonates across the study area, are thus constrained, especially during the Messinian desiccation stage. Application of the results to hydrogeology leads to the following :- Uplift of upstream hinterland (Cretaceous, Miocene) produces large amplitude karsts, but they do not make significant reservoirs.- Karstification responsible for the development of important and deep karst reservoirs is related to desiccation of endoreic basins (early Paleocene and Messinian).- Theses two events have produced two main deep karst reservoirs, with a gravimetric signature : (1) Upstream of Montpellier Thrust, between Pic St Loup, Les Matelles Fault and Hérault Basin, a deep karsts developed over 900m thickness in the Malm limestone, during Paleocene. (2) Downstream (south) of Montpellier Thrust, the coastal karsts extend offshore over more than 10km, they are sealed by a thick Miocene to Quaternary cover, and they result from the Messinian event.