En géosciences, la majeure partie du sous-sol est inaccessible à toute observation directe. Seules des informations parcellaires ou imprécises sont donc disponibles lors de la construction ou de la mise à jour de modèles géologiques ; de ce fait, les incertitudes jouent un rôle fondamental en géomodélisation. La théorie des problèmes inverses et les méthodes de simulations stochastiques fournissent un cadre théorique permettant de générer un ensemble de représentations plausibles du sous-sol, également nommées réalisations. En pratique, la forte cardinalité de l'ensemble des réalisations limite significativement tout traitement ou interprétation sur le modèle géologique.L'objectif de cette thèse est de fournir au géologue des algorithmes de visualisation permettant d'explorer, d'analyser et de communiquer les incertitudes spatiales associées à de larges ensembles de réalisations. Nos contributions sont les suivantes : (1) Nous proposons un ensemble de techniques dédiées à la visualisation des incertitudes pétrophysiques. Ces techniques reposent sur une programmation sur carte graphique (GPU) et utilisent une architecture garantissant leur interopérabilité ; (2) Nous proposons deux techniques dédiées à la visualisation des incertitudes structurales, traitant aussi bien les incertitudes géométriques que les incertitudes topologiques (existence de la surface ou interactions avec d'autres surfaces) ; (3) Nous évaluons la qualité des algorithmes de visualisation des incertitudes par le biais de deux études sur utilisateurs, portant respectivement sur la perception des méthodes statiques et par animation. Ces études apportent un éclairage nouveau sur la manière selon laquelle l'incertitude doit être représentée / Most of the subsurface is inaccessible to direct observation in geosciences. Consequently, only local or imprecise data are available when building or updating a geological model; uncertainties are therefore central to geomodeling. The inverse problem theory and the stochastic simulation methods provide a framework for the generation of large sets of likely representations of the subsurface, also termed realizations. In practice, however, the size of the set of realizations severely impacts further interpretation or processing of the geological model.This thesis aims at providing visualization algorithms to expert geologists that allow them to explore, analyze and communicate on spatial uncertainties associated to large sets of realizations. Our contributions are: (1) We propose a set of techniques dedicated to petrophysical uncertainty visualization, based on a GPU programming approach that maintains their interoperability; (2) We propose two techniques dedicated to structural uncertainty visualization that can handle both geometrical and topological uncertainties (e.g., the existence of the surface or its relationships with other surfaces); (3) We assess the quality of our uncertainty visualization algorithms through two user studies, which respectively focus on the perception of static and animated methods. These studies bring new elements on how uncertainty should be represented

La construction du réseau de failles d'un domaine géologique à l'aide d'un logiciel de géomodélisation peut être longue et fastidieuse pour des géométries et des contacts de failles complexes et peu marqués sur les données de subsurface. En utilisant une représentation implicite de ces surfaces de faille par des champs de potentiel calculés sur un maillage tétraédrique, le processus de création d'un modèle structural a été totalement automatisé. Les contacts de failles sont détectés et des surfaces avec une géométrie cohérente au niveau des contacts sont construites. Les champs de potentiel permettent la modification du réseau de failles tout en préservant la cohérence géologique du modèle. Une technique développée ici permet également de corriger la géométrie d'un maillage tétraédrique à proximité des failles. Ainsi, de nouvelles informations ou interprétations ou des perturbations représentant les incertitudes sur la position des failles sont intégrées au modèle géologique / One of the first steps when a geological study area is modelled in three-dimensional geomodelling software is to build the fault network. This can be tedious and time-consuming when fault geometry and branching are complex and difficult to locate from sub-surface data. The process of creating a three-dimensional structural model from various data types was entirely automated by using an implicit representation of fault surfaces by level sets computed on a tetrahedral mesh. Fault branching is detected automatically and surfaces with coherent contact geometry are built. Using level sets for fault surfaces also enables easy editing of the fault network while maintaining the geological consistency of the model. A further tetrahedral mesh editing technique was developed during this PhD in order to modify mesh geometry close to faults. New information or interpretations or perturbations which represent geometrical uncertainty on faults can thus be integrated to the geological model

Les réseaux de fractures et les karsts constituent des discontinuités au sein de la roche qui affectent considérablement les écoulements de fluides, ce qui engendre des problèmes spécifiques dans divers domaines des géosciences. La problématique générale consiste à déterminer les caractéristiques géométriques et hydrauliques des réseaux de fractures ou de karsts. La caractérisation et la modélisation de ces structures se heurtent cependant à leur complexité géométrique et à leur distribution spatiale hétérogène. De plus, les observations et données directes concernant aussi bien les fractures et karsts que leur encaissant rocheux restent largement insuffisantes pour décrire avec certitudes leurs caractéristiques. Pour ces raisons, la modélisation de réseaux de fractures ou de karsts est le plus souvent réalisée dans un cadre probabiliste. Des simulations stochastiques de type objet ou pixel sont généralement mise en œuvre pour générer des modèles 3D de fractures ou karsts. Cependant les mécanismes sur lesquels repose ce type d'approche ne permet pas de reproduire toutes la complexité de ces objets naturels et fournit des modèles manquant de réalisme géologique.Dans ces travaux de thèse, nous proposons d'aborder la problématique de la modélisation des fractures et des karsts suivant une approche pseudo-génétique. Il s'agit de contraindre le processus de simulation stochastique de fractures et karsts par des règles géométriques et heuristiques qui imitent les processus physiques gouvernant leur formation. Deux méthodes poursuivant cet objectif ont été développées, l'une adressant la simulation des fractures et la seconde celle des karsts. Les modèles ainsi générés exposent des caractéristiques similaires à celles des réseaux de fractures et karsts naturels. / Fractures and karstic networks are known to significantly affect flow paths and therefore raise specific issues in a wide variety of geoscience fields. The common question beyond these problems isto determine whether there is a network of fractures and/or karstic conduits and if yes what are itsgeometrical and hydraulic characteristics. Characterization and modeling of these features is a challenge for it usually displays complex geometries and heterogeneous spatial distribution. Moreover, in most cases, neither fracture and karst nor their host environment can be observed or described with certainty at all scales and location of relevance. For these reasons, fractures and karstic networks are usually integrated into 3D geological model through a probabilistic framework. Stochastic object- or pixel-based simulations are commonly performed to generate 3D models of fractures and karst but failed to reproduce the whole complexity of these natural objects and 3D models often lack geological realism.To address the issues related to fracture and karst modeling, we present two genetic-like approaches. The motivation of this work is to constrain the stochastic simulation of fractures and karsts by geometrical and heuristic rules which mimic the physical processes governing their formation. The resulting fracture and karst models display similar characteristics as those of natural pattern.

La modélisation du sous-sol est un outil indispensable pour décrire, comprendre et quantifier les processus géologiques. L'accès au sous-sol et son observation étant limités aux moyens d'acquisition, la construction de modèles tridimensionnels du sous-sol repose sur l'interprétation de données éparses à résolution limitée. Dans ce contexte, de nombreuses incertitudes affectent la construction de tels modèles, dues aux possibles biais humains cognitifs lors de l'interprétation, à la variabilité naturelle des objets géologiques et aux incertitudes intrinsèques des données utilisées. Ces incertitudes altèrent la prédictibilité des modèles et leur évaluation est donc nécessaire afin de réduire les risques économiques et humains liés à l'utilisation des modèles. Le travail de thèse s'est déroulé dans le cadre plus spécifique des incertitudes sur les structures géologiques. Les réponses apportées sont multiples : (1) une méthode stochastique de génération de modèles structuraux à géométrie et topologie changeantes, combinant une connaissance a priori des structures géologiques aux données interprétées, a été développée ; (2) le réalisme géologique des structures modélisées est garanti grâce à la modélisation implicite, représentant une surface par une équipotentielle d'un champ scalaire volumique ; (3) la description des failles en un nombre restreint de paramètres incertains a permis d'aborder la modélisation inverse, ce qui ouvre la voie vers l'assimilation de données géophysiques ou d'écoulement fluides grâce à des méthodes bayesiennes / Subsurface modeling is a key tool to describe, understand and quantify geological processes. As the subsurface is inaccessible and its observation is limited by acquisition methods, 3D models of the subsurface are usually built from the interpretation of sparse data with limited resolution. Therefore, uncertainties occur during the model building process, due to possible cognitive human bias, natural variability of geological objects and intrinsic uncertainties of data. In such context, the predictability of models is limited by uncertainties, which must be assessed in order to reduce economical and human risks linked to the use of models. This thesis focuses more specifically on uncertainties about geological structures. Our contributions are : (1) a stochastic method for generating structural models with various fault and horizon geometries as well as fault connections, combining prior information and interpreted data, has been developped ; (2) realistic geological objects are obtained using implicit modeling that represents a surface by an equipotential of a volumetric scalar field ; (3) faults have been described by a reduced set of uncertain parameters, which opens the way to the inversion of structural objects using geophysical or fluid flow data by baysian methods

La caractérisation sismique des réservoirs pétroliers nécessite l'intégration de plusieurs techniques telles que la lithosismique, la géomodélisation, la géostatistique, l'utilisation des algorithmes évolutionnaires et la pétrophysique. L'information sismique est d'abord utilisée pour la description de l'architecture externe des réservoirs car son utilisation pour la description des faciès ne se fait pas sans difficultés. L'objectif de cette thèse est d'apporter des outils nouveaux pour aider à l'utilisation de l'information sismique pour caractériser les réservoirs.<br />Un premier travail a consisté à évaluer l'impact des incertitudes structurales sur les inversions pétroélastiques et les conséquences en terme de classification de faciès. Ensuite, nous considérons la modélisation sismique comme aide à l'évaluation du modèle réservoir. Cette modélisation permettra de faire le lien entre les simulateurs réservoir ou les géomodeleurs et la réponse sismique du réservoir. <br />Nous développons ensuite deux approches alternatives aux méthodes traditionnelles en inversion pétroélastique et pétrophysique. La première utilise la méthode géostatistique des déformations graduelles pour créer des réalisations de propriétés réservoirs. Elle permet de créer des propriétés à l'échelle réservoir, conditionnées aux puits, tout en respectant une fonction coût basée sur la comparaison des données sismiques réelles et issues de ces réalisations. <br />La seconde méthode repose sur le principe de la classification supervisée et utilise des réseaux de neurones pour analyser la forme des traces sismiques. Une première étape consiste à générer un volume d'apprentissage contenant tous les modèles pétrophysiques envisageables pour un champ donné. Ces modèles sont analysés par les réseaux de neurones. Les neurones ainsi identifiés sont appliqués aux données réelles, pour identifier des relations pétrophysique/sismique identiques aux données d'apprentissage.<br />Toutes les méthodologies sont validées sur plusieurs réservoirs choisis pour leurs particularités géologiques (complexité structurale, lithologie du réservoir).

La corrélation entre les différents puits disponibles des unités stratigraphiques identifiées correspond à l'une des premières étapes de la construction d'un modèle de sous-sol. Les horizons stratigraphiques construits définissent, de part leur géométrie et leur topologie, l'architecture stratigraphique du sous-sol guidant simulations géostatistiques des propriétés réservoirs et des faciès. Cependant, la faible qualité et/ou quantité des données disponibles ainsi que la complexité de la mise en oeuvre des concepts sédimentologiques utilisés pour construire un modèle de corrélation rendent cette étape incertaine. Afin de gérer et d'échantillonner ces incertitudes, une méthode stochastique de corrélation stratigraphique est proposée. Cette méthode se base sur le développement de règles de corrélation, assises sur des concepts sédimentologiques ainsi que sur les connaissances a priori de la zone étudiée. Cette méthodologie est appliquée à la corrélation de séquences stratigraphiques de la de marge carbonatée, d'âge Crétacé Supérieur, du bassin Sud-Provençale. Dans cette optique, deux règles de corrélations, basées sur la cohérence des paléoangles et sur une représentation de la paléotopographie construite a priori sont définies. L'étude de la répartition des faciès dans les différents modèles de corrélations générés indique que les incertitudes sur les corrélations stratigraphiques conduisent à envisager plusieurs scenarios de compartimentation de réservoir. Une étude liant corrélations stratigraphiques, simulations d'écoulement et géophysique est menée sur le champ de gaz de Malampaya. Dans ce cas, des unités diagénétiques sont corrélées afin de construire différents modèles statiques et dynamiques de ce réservoir. Cette étude permet de montrer un fort contrôle stratigraphique sur la répartition de l'écoulement alors que l'imagerie séismique semble moins affectée. L'étude de dépôts fluviaux deltaïques au niveau du bassin de la mer du Nord sert de base au développement de règles de corrélations intégrant des informations issues de l'imagerie séismique. Ces nouvelles règles de corrélations permettent de mieux contraindre le problème de corrélations stratigraphique lorsque les sédiments enregistrés au niveau des puits montrent une forte cyclicité. La magnétostratigraphie est une autre source d'information permettant d'étudier l'histoire du remplissage et de la déformation de bassins sédimentaires. Le développement de règles dédiées à cette discipline, et l'utilisation de celles-ci sur des séries sédimentaires himalayennes, permet d'appréhender les incertitudes sur l'âge des sédiments ainsi que sur les paléo-taux d'accumulation de sédiments / Stratigraphic correlation consists in linking boundaries of correlative units between wells or outcrops over a given study area, and is therefore one of the first steps of the characterization of the subsurface geometry, supporting geostatistical modeling of static reservoir properties. However, this early step is subject to uncertainties, since stratigraphic well correlation is constrained only by sparse observations (wells and outcrops), low resolution information coming from geophysics, regional knowledge and geological concepts. The Dynamic Time Warping (DTW) algorithm serves as a base for the development of a generic method stochastically performing stratigraphic correlation between units identified on available wells. The proposed method relies on the definition of correlation rules that are applied using the available data and some regional knowledge, according to the way stratigraphic units are defined. An application to the Cretaceous southern Provence carbonate basin has been performed using correlation rules based on paleo-angles and the theoretical architecture of the depositional environment. The computation of vertical proportions of facies on numerous models generated from the stochastic correlations of sequence stratigraphic units indicates that the uncertainties on the stratigraphic correlation impact the compartmentalization of the modeled reservoirs. The impact of stratigraphic correlation uncertainties on fluid flow behavior is assessed through the example of the Malampaya diagenetic carbonate reservoir. Diagenetic units are correlated on the basis of their wireline log signature and diagenetic types. Different models are generated from the stochastic well correlations, and the corresponding water saturation profiles are computed. They show different displacement patterns, indicating a stratigraphic control of the dynamic property, which contrasts with the synthetic seismic model constructed from the corresponding geomodel. Magnetostratigraphic correlation is another way to study sedimentary basin deposition and deformation history. Adapting the DTW algorithm to magnetostratigraphic data, we generate dating models of Himalayan deposits, for which conflicting interpretations are proposed in the literature. This allows managing the associated accumulation rates uncertainties

Les orogènes à double vergence sont classiquement définis comme deux prismes critiques opposés (pro et retro) qui évoluent ensemble. Les études récentes montrent que les rétro-prismes et leurs bassins d’avant-pays associés se comportent différemment des pro-prismes. Cependant, ni les facteurs qui mènent un orogène à devenir doublement vergent, ni la relation entre le pro- et rétro-prisme ne sont bien compris. Le but de cette étude est d'améliorer notre connaissance 1) de la relation entre le pro- et le rétro-prisme pendant l'orogénèse, 2) des facteurs contrôlant l'évolution d'un orogène à double vergence, et 3) d’un lien dynamique possible entre le pro- et le rétro-prisme. Répondre à ces questions nécessite une connaissance améliorée de l'évolution d'un orogène à double vergence. Nous nous sommes concentrés sur les Pyrénées Orientales, en raison de la grande quantité de données disponibles. Nous avons effectué une étude de terrain tectono-stratigraphique détaillée à l’est du Massif de Saint Barthelemy et dans l’avant-pays autour de Lavelanet (plaque Européenne). Notre interprétation d’une coupe restaurée intègre une configuration crustale pré-orogenique en tant qu'une marge hyper-amincie. Nous relions l'évolution détaillée du rétro-prisme à celle du pro-prisme (plaque Ibérique), afin de mieux contraindre la dynamique à l'échelle crustale. Nous subdivisons l'évolution des Pyrénées Orientales en quatre phases. La première phase (Crétacé Supérieur) est caractérisée par la fermeture d'un domaine de manteau exhumé entre les plaques et l'inversion synchrone d'un système de rift riche en sel et thermiquement déséquilibré. Le raccourcissement était distribué de façon égale entre les deux marges pendant cette première phase d’inversion. Une phase de quiescence (Paléocène), limitée au rétro-prisme, enregistre la transition entre l'inversion et la phase de collision. La phase de collision principale (Éocène) enregistre le taux de raccourcissement le plus élevé, et était principalement accommodé dans le pro-prisme. Pendant la phase finale (Oligocène) le rétro-prisme était largement inactif et le raccourcissement du pro-prisme a ralenti. Cela démontre que la relation entre le pro- et rétro-prisme change avec le temps. Nous avons utilisé des modèles numériques 2D thermomécaniques à l'échelle lithosphérique pour simuler l'évolution d'un orogène à double vergence s'initie après avec un rift. Nos résultats montrent un modèle évolutif similaire à celui observé dans les Pyrénées Orientales avec une phase d'inversion du rift approximativement symétrique suivie d'une phase de collision asymétrique. L'héritage du rift est essentiel pour permettre le développement d’un orogène à double vergence. Des autres facteurs, comme les processus de surface et la déformation de la couverture, ont un effet significatif sur la structure crustale et la répartition du raccourcissement entre les deux prismes. Un niveau de décollement (sel) à la base de la couverture favorise la formation d'un empilement antiformal d’écailles crustales, similaire à la géométrie observée dans la Zone Axiale des Pyrénées, en formant un prisme à faible pente qui force la déformation crustale à se concentrer dans l'arrière-pays. Enfin, nous montrons que l'évolution des pro- et rétro-prismes est inextricablement liée : des événements ou des conditions d'un côté de l'orogène ont un effet direct sur l'autre côté de l'orogène. Ceci est clairement démontré dans nos modèles par des variations constantes des taux de raccourcissement du pro- et rétro-prisme en réponse à l'accrétion dans le pro-prisme. Le Haut Atlas (Maroc) et Pyrénées peuvent être respectivement considérés comme des exemples d'inversion de rift symétrique et de phases de collision asymétrique ultérieures / The doubly vergent nature of some natural orogens is classically understood as two opposing thrust wedges (pro and retro) that comply with critical taper theory. The evidence that retro-wedges and their associated basins behave differently from their pro-wedge counterparts has been steadily increasing over the past few decades. However, what causes an orogen to become doubly vergent is currently not well understood. Nor is the relationship between the pro- and retro-wedge during the evolution of a doubly vergent orogen. It is the aim of this work to improve our understanding of: 1) how the pro- and retro-wedges relate to each other during the orogenic process, 2) what factors control the evolution of a doubly vergent orogen and 3) a possible link between the pro- and retro-wedge. Answering these questions requires an improved knowledge of the evolution of a doubly vergent orogen. We focussed on the Eastern Pyrenees as a type example of a doubly vergent orogen, due to the large amount of available data. We performed a detailed tectonostratigraphic study of the retro-foreland of the Eastern Pyrenees (European plate), updating the interpretation based on recent insights into its hyperextended rift origins. We link the evolution of the retro-foreland to that of the pro-foreland (Iberian plate) in order to derive insight into the crustal scale dynamics. Based on cross section restoration, reconstructed shortening rates and subsidence analysis, we subdivide the East Pyrenean evolution into four phases. The first (Late Cretaceous) phase is characterised by closure of an exhumed mantle domain between the European and Iberian rifted margins, and simultaneous inversion of a salt-rich, thermally unequilibrated rift system. Shortening was distributed roughly equally between both margins during this early inversion phase. Following inversion, a quiescent phase (Paleocene) was apparently restricted to the retro-foreland. This phase may record the period of transition between inversion and full collision in the Eastern Pyrenees. The main collision phase (Eocene) records the highest shortening rates, which was predominantly accommodated in the pro-wedge. Retro-wedge shortening rates were lower than during the rift inversion phase. During the final phase (Oligocene) the retro-wedge was apparently inactive and shortening of the pro-wedge slowed. This demonstrates that the relationship between the pro- and retro-wedges changes through time. We used lithosphere-scale thermo-mechanical numerical models to simulate the evolution of a doubly vergent orogen. Our results show a similar evolutionary pattern as observed in the Pyrenees: A roughly symmetrical rift inversion phase is followed by an asymmetric collision phase. Rift inheritance was found to be essential for enabling double vergence. Other factors, such as surface processes and thin-skinned deformation, were found to have a significant effect on the crustal structure and strain partitioning between both wedges. A salt décollement layer in the sedimentary cover promotes the formation of a crustal antiformal stack such as observed in the Pyrenees and Alps by forming a wide and low-taper thin-skinned fold-and-thrust belt that forces crustal deformation to focus in the hinterland. Finally, we show that the evolution of the pro- and retro-wedges is inextricably linked: events or conditions on one side of the doubly vergent orogen have an immediate effect on the other side of the orogen. This is clearly demonstrated in our models by constant variations in shortening rates of the pro- and retro-wedge in response to accretion of new pro-wedge thrust sheets. The High Atlas (Morocco) and Pyrenees can be seen as examples of symmetric rift inversion and later asymmetric collision phases, respectively

La Zone Non Saturée (ZNS) peut jouer un rôle majeur pour la circulation et le stockage de l'eau dans les hydrosystèmes karstiques mais ces mécanismes restent peu connus. L'hydrogéophysique peut-elle permettre d'évaluer les variations spatiales et temporelles du stock de l'eau de la ZNS pour améliorer la compréhension de son fonctionnement ? En 2011 et 2012 différents outils ont été utilisés sur un site expérimental à l'aplomb du Laboratoire Souterrain à Bas Bruit de Rustrel (LSBB) dans le bassin d'alimentation de la Fontaine de Vaucluse afin d'en étudier la structure et le fonctionnement hydrogéologique. Ces investigations se sont déroulées selon une méthodologie comprenant quatre phases : (1) phase préliminaire d'étude géologique et hydrogéologique (2) caractérisation géophysique de la stucture géologique (3) caractérisation hydrogéophysique du fonctionnement hydrodynamique (4) synthèse par géomodélisation 3 D. Les résultats obtenus sur le site expérimental soulignent le potentiel de l'hydrogéophysique en milieu karstique et valident la méthodologie proposée. Les investigations ont permis d'identifier rapidement des structures géologiques pouvant jouer un rôle majeur dans la régulation des flux d'eau. Des variations saisonnières du stock d'eau ont été mesurées jusqu'à plusieurs dizaines de mètres dans la ZNS. Le suivi géophysique évènementiel a permis d'imager l'hétérogénéité spatiale et temporelle de la recharge. Enfin, la géomodélisation s'avère essentielle dans la méthologie proposée pour synthétiser les connaissances du site, finaliser le schéma de fonctionnement et poser les bases de futurs travaux de modélisation hydrodynamique dans la ZNS du karst. / Karst Unstatured Zone can play a major hydrodynamic role but water flow and storage mechanisms remains poorly known. Can hydrogeophysics help assessing the spatial and temporal waterstorage variation in UZ and improve understanding karst UZ functioning ? In 2011 and 2012, various tools have been employed on an experimental site above the low noise underground research laboratory of Rustrel (LSBB) within the Fontaine de Vaucluse catchment, in order to study the hydrogeological structure and functioning. These investigations were conducted using a methodology including four phases : (1) Preliminary study of geological and hydrogeological contexte, (2) Geophysical characterization of the geological features, (3) Hydrogeophysical characterization of the hydrodynamic functioning, (4) 3D geomodeling of the karst UZ. The acquired results highligth the potential of hydrogeophysical investigation in karst media and validated the proposed methology. These investigations allowed a fast identification of geological structures which may play a major role in water flow regulation. Seasonal variations in water storage were measured until several tens of meters deep within the UZ. Hydrogeophysical monitoring allowed imaging spatial and temporal heterogeneityof recharge; Finally, the geomodeling is essential within the proposed methology to synthesize site knowledge, finalize the conceptual functioning schema and lay the foundation for future work of hydrodynamicmodeling in karst UZ / La Zona No Saturada (ZNS) puede jugar un papel importante en la circulación y el almacenamiento del agua en los hidrosistemas kársticos, sin embargo sus mecanismos son poco conocidos. ¿La hidrogeofísica podría evaluar las variaciones espaciales y temporales del agua subterránea de la ZNS para mejorar la comprensión de su comportamiento?En el año 2011 y 2012, diferentes herramientas fueron utilizadas en un sitio experimental del Laboratorio Subterráneo de Bajo nivel de Ruido (LSBB, siglas en francés) ubicada en la cuenca de alimentación de la Fontaine de Vaucluse (Francia) con el propósito de estudiar su estructura y comportamiento hidrogeológico. Estas investigaciones fueron desarrolladas según una metodología compuesta de cuatro fases: (1) una fase preliminar de estudio geológico e hidrogeológico, (2) una caracterización geofísica de la estructura geológica, (3), una caracterización hidrogeofísica del comportamiento hidrodinámico, y (4) una síntesis por geomodelización 3D.Los resultados obtenidos en el campo validan la metodología utilizada y muestran la importancia del uso de la hidrogeofísica en el estudio del medio kárstico. Las investigaciones permitieron identificar rápidamente las estructuras geológicas que juegan un rol importante en la regulación de los flujos de agua. Las variaciones temporales del agua subterránea fueron medidas hasta varias decenas de metros de profundidad en la ZNS. El monitoreo geofísico de eventos puntuales permitió realizar una representación grafica de la heterogeneidad espacial y temporal de la recarga. Finalmente, la geomodelización es considerada como parte esencial en la metodología propuesta para sintetizar los conocimientos del área de estudio, finalizar el esquema del funcionamiento y proponer las bases de los futuros trabajos en modelización hidrodinámica de la ZNS del karst. / يمكن للطبقة غير المشبعة٬ في نظم المياه الجوفية ذات الطبيعة الجيرية (الكارست) أن تلعب دورا مهما في حركة و تخزين المياه الجوفية٬ لكن آليات عمل هذه الطبقة تبقى إلى حد ما مبهمة. إن السؤال الذي يبقى مطروحا : هل بإمكان علم الهيدروجيوفيزياء أن يمكننا من تقييم التغيرات في البعدين المكاني و الزمني لمخزون المياه في الطبقة غير المشبعة بغية تطوير فهمنا لآلية عملها؟خلال الفترة الممتدة بين سنتي 2011 – 2012 م٬ تم إستخدام العديد من الأدوات في حقل تجريبي على مستوى المختبر التحت أرضي عديم التشويش بروستال بفرنسا (LSBB de Rustel)٬ والذي يتواجد في حوض التغذية لمنبع الفوكلوز(Fontaine de Vaucluse) بغية دراسة بنية الطبقة غير المشبعة للكارست و كيفية أدائها فيما يتعلق بالمياه الجوفية. إن هذه الأبحاث نفذت بمنهجية مكونة من أربعة مراحل٬ و هن كالآتي: (1) المرحلة الابتدائية و تخص الدراسة الجيولوجية و الهيدروجيولوجية٬ (2) التوصيف الجيوفيزيائي للبنية الجيولوجية٬ (3) التوصيف الهيدروجيوفيزيائي لحركة الطبقة المائية٬ (4) الخلاصة من خلال إنشاء نموذج جيولوجي ثلاثي الأبعاد.إن النتائج المتحصل عليها في الموقع التجريبي أظهرت أهمية الهيدروجيوفيزياء في الوسط الكارستي و أكدت نجاعة المنهجية المقترحة. إن الأبحاث المنجزة٬ مكنت في وقت وجيز٬ من تحديد البنيات الجيولوجية التي يمكن أن تلعب دورا مهما في معادلة التدفقات المائية. من منظور آخر٬ مكنتنا الطريقة المقترحة من قياس التغيرات الموسمية للمخزون المائي في الطبقة غير المشبعة (ZNS) و ذلك على عمق عشرات الأمتار. إن الرصد الجيوفيزيائي الحدثي مكن من تصوير التغيرات المكانية و الزمنية لإعادة التغذية. أخيرا يظهر جليا أن النمذجة الجيولوجية ثلاثية الأبعاد ذات الأسس الجيولوجية٬ الهيدروجيولوجية و الجيوفيزيائية٬ أمر أساسي بالنسبة للمنهجية المقترحة من أجل حوصلة معلوماتنا بخصوص الموقع المدروس٬ تكملة مخطط التشغيل و وضع أسس الأعمال المستقبلية لنماذج حركة المياه خلال الطبقة غير المشبعة لأحواض المياه الجوفية ذات الطبيعة الجيرية (الكارست Karst).المفردات الرئيسية:الهيدروجيوفيزياء (علم جيوفيزياء المياه الجوفية)، الكارست (حوض مائي جوفي ذو طبيعة جيرية)، إعادة تغذية الطّبقة المائيّة، الطّبقة غير المشبّعة (ZNS)، النّمذجة الجيولوجيةGéomodélisation) )، منبع الفوكلوز (ontaine deF Vaucluse)، المختبر التحت أرضي عديم التشويش بروستال بفرنسا (LSBB de Rustel). / 非饱和带对喀斯特地区水分的运动和存储具有重要作用，但是这一机制仍然不太清楚。水文地球物理学方法是否有助于评估非饱和带水分存储的时空变异，并深化我们对喀斯特非饱和带功能的认识？为了研究沃克吕兹泉流域的水文地质结构和功能，分别于2011和2012年，在吕斯特勒地区地下低噪声试验室(lsbb)上方的试验点进行了各种试验。研究内容包括四个部分：(1)：初步的地质和水文地质调查；(2)：地质结构的地球物理方法勘探；(3)：水动力过程的水文地球物理监测；(4) 喀斯特非饱和带结构的三维模拟。研究结果表明，水文地球物理学方法在喀斯特地区的研究中具有很大应用潜力，并检验了我们提出的方法的可靠性。利用这种方法，对对水流调节起主要作用的地质结构进行了快速识别；对深达几十米的非饱和带储水量的季节性变化进行了监测；对地下补给过程的时空异质性进行了成像。最后，综合研究区的各种信息，建立地学模型，确定了水文地质结构和功能的概念模式，为后续喀斯特非饱和带水动力学过程模拟研究奠定了基础。

La géomodélisation du sous-sol est considérée comme stratégique par les compagnies pétrolières. Chacune a en effet pour objectif de construire les meilleurs modèles géologiques sur lesquels seront évalués les réserves de gaz et d'huile des territoires étudiés, ainsi que le placement optimal des puits de forage qui exploiteront ces réserves. Ces modèles sont bâtis sur l'interprétation de la sismique par les géologues et les géophysiciens. Le processus de modélisation consiste à structurer le sous-sol en un ensemble d'horizons (strates géologiques), de failles (fractures de la roche) et de corps (amas de roches de même nature). Ces objets géologiques sont au sein de deux problématiques complémentaires abordées dans cette thèse chez Total : la construction du modèle structural à proprement parler constitué d'un réseau de failles et d'horizons mis en cohérence et la constructions des corps selon des attributs sismiques élaborés. La construction du modèle structural part des données issues d'interprétations manuelles et/ou semi-automatiques ne décrivant que ponctuellement les objets géologiques : un ensemble de lignes brisées pour les failles et un nuage de points pour les horizons. Une première partie de nos travaux porte sur l'élaboration d'une méthodologie complète de construction de surfaces modélisant les failles et les horizons adaptée à la nature et la densité particulières des données initiales. Dans un premier temps les failles sont reconstruites indépendamment les unes des autres à l'aide d'une triangulation de Delaunay contrainte des données projetées dans leur plan de régression, puis raffinées par subdivision. Ensuite un réseau de failles global est créé et mis en cohérence par des opérations d'extrapolation et de clipping. Enfin, les trous dans les horizons sont bouchés par des méthodes d'interpolation spatiale étendus par nos soins à l'aide des "voisinages partiels" contraints par les failles, pour lesquelles nous montrons l'influence des stratégies de remplissage et leur robustesse au bruit. Les corps sont construits par extraction d'isosurfaces sur des volumes d'attributs sismiques pouvant atteindre des tailles de plusieurs dizaines de giga-octets. Ces tailles de données nous ont imposées de développer un algorithme parallèle d'extraction d'isosurfaces simplifiées. Le volume est d'abord récursivement découpé pour former une structure hiérarchique de blocs voisins. Des isosurfaces partielles sont extraites de ces blocs par des processeurs indépendants exécutant une version étendue de l'algorithme tandem. Ainsi, les noeuds de calcul extraient, simplifient, s'échangent et assemblent les morceaux de surfaces ainsi obtenues jusqu'à l'obtention d'une isosurface simplifiée globale. Cette dernière est organisée et stockée par construction en un flux de composantes connexes qui permettra une interaction ultérieure de haut niveau.

Les calcaires hémipélagiques du Sud-Est de la France sont caractérisés par de faibles valeurs de porosité et de perméabilité. Ces caractéristiques permettent de considérer ces roches comme de potentiels analogues des réservoirs pétroliers non-conventionnels ou des aquifères fissurés. Ces réservoirs sont d’une importance capitale dans le contexte économique actuel. Néanmoins, l’exploitation de ce type de réservoir représente des risques industriels et des challenges techniques qui doivent être gérés pour réduire leur empreinte sur l’environnement.Ce travail apporte une quantité significative de nouvelles données et des méthodes pour caractériser les réservoirs carbonatés non-conventionnels. Il démontre l’importance des figures de compaction (stylolites, microstylolites et faisceaux de pression-dissolution) sur la capacité de stockage des fluides dans ce type de réservoir. Il montre aussi que l’évolution diagénétique des sédiments carbonatés hémipélagiques dépend de processus précoces (cimentations précoces, compaction différentielle) et éventuellement de réactivations plus tardives (inversion de contraintes tectoniques). Le couplage des modèles 3D avec les géostatistiques apparaît comme un protocole valide et rapide pour identifier les éléments clef (faciès, fractures et stylolites) qui jouent un rôle dans la distribution des hétérogénéités pétrophysiques des réservoirs non-conventionnels. Cette méthode permet aussi de générer des modèles équiprobable de corps géologiques fracturés et perméables permettant de possible investigations sur la distribution des zones d’écoulement préférentiels. / Hemipelagic limestones from the South-East of France are characterized by low values of porosity and permeability that allow considering these rocks as potential analogs to unconventional oil and gas reservoirs or fissured aquifers. These reservoirs are of primary importance in today’s economical context where the energy and the water’s demands grow with the increase in the world population. Nevertheless, the industrial exploitation of such reservoirs faces series of environmental risks and technical challenges that require to be managed in order to reduce their footprint on the environment.This work brought a significant amount of new data on unconventional carbonate reservoirs. It highlights how important compaction features (stylolites, microstylolites, solution seams) may be to fluids storage in such reservoirs. It shows that the diagenetic evolution of hemipelagic carbonate sediments depends on both early processes (sorting, early cementation, differential compaction) and eventual late reactivations (under tectonic stress inversion for example). Coupling three-dimensional modeling with geostatistics appear as an accurate and fast protocol to identify key elements (facies, fractures and stylolite) that play a role in the petrophysical heterogeneity of such unconventional reservoirs. This method also allows identifying equiprobable models of fractures geobodies and investigating the distribution of flow preferential pathways.