Thèse--Techniques et économie de l'entreprise minière--Paris--Ecole nationale supérieure des mines de Paris, 1987. / Bibliogr. p. 171-179.

L'effet mécanique des forces de succion, forces exercées par un fluide qui se déplace dans un milieu poreux, a été étudié dans le cadre de la fracturation hydraulique des roches de la croûte terrestre. Cet effet a été étudié par des méthodes analytiques issues de la mécanique des milieux continus, et par des simulations numériques. Ces forces de succion sont des forces de frottement causées par des gradients de pression fluide dans les milieux poreux. Différents modes de fracturation (bandes de cisaillement localisées, fractures en mode I) causés par une augmentation localisée de la pression fluide dans la croûte ont été reproduits dans un milieu poro-élastique grâce à plusieurs codes numériques spécialement développés à cet effet. La valeur de la pression fluide lors de la nucléation de la fracturation est aussi prédite à l'aide d'une nouvelle solution analytique. <br />Dans la solution analytique, une cavité elliptique dans un solide poreux est remplie avec un fluide à une pression non-hydrostatique. On considère aussi que le milieu poreux est soumis à un champ de contrainte externe. Puisque la pression du fluide dans la cavité est différente de la pression de pore dans la roche; le couplage poro-élastique est pris en compte dans le calcul des déformations. A partir de la théorie de Griffith qui donne une condition pour la propagation d'une fracture, et en utilisant la solution analytique obtenue, une équation généralisée a été obtenue pour la contrainte effective dans le milieu. Cette nouvelle loi décrit la fracturation dans un milieu poreux saturé avec un fluide, et dans lequel la distribution de pression fluide n'est pas homogène.

La fracturation des roches concerne plusieurs applications industrielles, telles que les tunnels, l’excavation, l’extraction de blocs de roches et l'industrie minière. Les techniques classiques pour fracturer la roche sont associées à des inconvénients ceux qui restreignent leur utilisation dans certaines circonstances. Par conséquent, la fracturation par spallation thermique a été suggérée comme une solution alternative. Cela produit des contraintes de compression dans la roche grâce à un flux de chaleur intense. Une fois que ces contraintes dépasseraient la résistance à compression de la roche, la rupture aurait lieu en éjectant des petits éclats. Cette mode de fracturation thermique nécessite une source d’énergie puissante ainsi qu’elle influence négativement les conditions du travail sur le site. Profitant d'une faible résistance à la traction du matériau de roche, le présent travail propose une méthode basée sur un chauffage rapide de la matière rocheuse par un rayonnement micro-ondes, suivi par un refroidissement rapide par de l'eau pour fracturer la roche. Au cours de ce processus, une dilatation thermique est produite à l'intérieur de la zone traitée, qui est ensuite suivie d'une contraction thermique. Cette contraction produit des contraintes de traction, ce qui pourrait conduire à la rupture de la matière de roche. Le processus d'échauffement/refroidissement sera analysés grâce à un modèle numérique ce qui décrit les événements physiques et l'influence de différents paramètres. Expérimentalement, le processus sera examiné du point de vue thermique.‎ / The issue of rock fracturing concerns several industrial applications, such as tunneling or excavation ‎in rock material, extraction of rock blocks and the mining industry. The conventional techniques to ‎fracture the rock material are associated with disadvantages restrict their use in certain circumstances. ‎Consequently, the fracture by thermal spallation was suggested as an alternative solution. This ‎solution basically depends on exceeding the compressive strength of the rock material by compressive ‎stresses induced by heating effects. This necessitates high temperature at the rock surface, which ‎requires high power supply and it degrades the work conditions. Profiting from low tensile strength ‎of the rock material, the present work proposes a method based on a rapid heating of the rock ‎material by microwaves radiation followed by a rapid cooling by water to fracture the rock material. ‎During this process, a thermal dilatation is produced within the treated zone, which is then followed ‎by a thermal shrinkage. This shrinkage produces tensile stresses, which could lead to fracturing the ‎rock material. The heating/cooling process will be analyzed through numerical model describes the ‎physical events and the influence of different parameters. Experimentally, the process will be ‎examined from thermal point of view.‎

La déshydratation, l’oxydation et la fracturation sont les trois principaux processus couplés subis par l’argilite callovo-oxfordienne recoupée par les galeries du Laboratoire Souterrain de Meuse/Haute-Marne de l'Andra lors de l'excavation de galeries. Une série de forages à différents temps après l'ouverture de la galerie a permis d'étudier l'évolution des perturbations en fonction du temps, jusqu'à 18 mois. Les profils de teneur en eau et les analyses isotopiques montrent une déshydratation en surface des parois d'environ 50%. Le front de déshydratation est rapidement sub-stationnaire car il résulte d'un équilibre entre la recharge par l'eau de la formation et l'évaporation imposée par la pression partielle de l'eau dans l'air de la galerie. La zone perturbée affecte principalement les premiers 70 cm, et ne se fait plus sentir à partir de 1 m 50. L’oxydation des pyrites provoque la formation de fluides sulfatés acides qui dissolvent les minéraux environnants, libèrent dans la solution des ions Ca2+, K+ et Na+ puis lors de la déshydratation de l’argilite, l'évaporation de l'eau porale induit la précipitation de minéraux sulfatés principalement. La dégradation de la matière organique dans les conditions des galeries montre une dégradation limitée, due soit à une oxydation très discrète, soit à une dégradation lente non détectable aux temps courts. Une série d'expérience de mise en équilibre d'eau avec une roche plus ou moins oxydée a permis d'évaluer l'impact de la présence des phases minérales néoformées sur les caractéristiques chimiques des eaux de réhydratation. Cette approche expérimentale a été confrontée aux résultats de modélisation géochimique. Les résultats obtenus permettent de mieux préciser les caractéristiques géométriques, minéralogiques et physico-chimiques de cette zone perturbée qui est susceptible d'évoluer par la suite en condition de réhydratation et d'élévation de température lors du stockage / Dehydration, oxidation, and fracturing are the three main and coupled processes underwent by Callovian-Oxfordian argillite cut by the galleries of the Laboratoire Souterrain of Meuse / Haute-Marne of Andra during the excavation. Several boreholes drilled at different times after the opening of the gallery were used to study the evolution of disturbances with time up to 18 months. Profiles of water content and isotopic analysis show dehydration at the surface of the walls of about 50%. The area of dehydration is rapidly sub-stationary because it results from a balance between the suction by the pore water and the evaporation imposed by the partial pressure of water in the air of the gallery. The disturbed area mainly affects the first 70 cm, and not exceeds 1 m 50. The oxidation of pyrite causes the formation of acid sulphate fluids which dissolve the surrounding minerals in the solution, release Ca2+, K+ and Na+ in the pore water then during the dehydration of argillite, the evaporation of the pore water leads to the precipitation of sulphate minerals mainly. The degradation of organic matter in the galleries shows a limited deterioration due to a very moderate oxidation or the slow degradation is not detectable at short times. Several experiments of equilibration between more or less oxidized samples and water were used to evaluate the impact of the presence of new-formed mineral phases (sulphates) on the chemical composition of the rehydration water. This experimental approach has been compared to results of geochemical modelling. The results obtained allow us to define geometrical, mineralogical and physico-chemical properties of the disturbed zone which could evolve thereafter during rehydration and temperature rise in storage condition

L'accrétion océanique au niveau des dorsales lentes est alimentée par un apport magmatique et par l'exhumation de péridotites mantelliques. Une partie de la chaleur dégagée lors du refroidissement et de la cristallisation des magmas sert de moteur à la circulation hydrothermale qui peut transporter de l'eau de mer jusqu'aux roches mantelliques. Une réaction d'oxydoréduction impliquant l'eau de mer, l'olivine et le pyroxène des péridotites peut alors avoir lieu pour produire de la serpentine, de la magnétite et de l'hydrogène. Cette réaction, dite de serpentinisation, contrôle les propriétés physico-chimiques de la lithosphère océanique et est, de fait, incorporée dans les modèles géophysiques avec une cinétique encore mal contrainte. Ce travail de thèse est centré sur la cinétique d'hydratation de l'olivine de San Carlos, une réaction de serpentinisation simple. Un dispositif expérimental de suivi de la cinétique d'hydratation a été mis au point, basé sur la mesure du signal magnétique associé à la magnétite formée. Cette méthode, couplée à une modélisation thermodynamique, permet également de quantifier indirectement la production d'hydrogène. L'influence de la taille initiale des grains l'olivine et de la température sur la cinétique d'hydratation a aussi été étudiée. Une augmentation de la vitesse de serpentinisation dans les premiers temps de la réaction est mise en évidence et expliquée par la mise en place de microtextures également observées dans le milieu naturel (puits de corrosion et fractures) contribuant à un accroissement de la surface réactive. De manière générale, la cinétique de serpentinisation obtenue ici est de un a deux ordres de grandeur plus lente que celle utilisée jusqu'alors dans les modèles géophysiques. En plus de cette étude expérimentale, une réaction d'oxydoréduction dans un contexte naturel a été étudiée au niveau d'une zone réactionnelle entre des serpentinites et des marbres en Corse alpine. Sur la base d'une approche pétrographique et thermodynamique, la présence de wollastonite, CaSiO3, est expliquée par l'atteinte de conditions très réductrices associées à la présence de fluides riches en hydrogène en contexte de subduction. Ces fluides sont interprétés comme résultant de la conservation du potentiel réducteur des serpentinites au cours de la subduction. / Slow-spreading ridges are fed by a magmatic input and the exhumation of mantle peridotites. Part of the magmatic heat is evacuated through the hydrothermal circulation which can transport seawater up to the mantle rocks. A RedOx reaction between seawater and the olivine and pyroxene from the peridotite can then take place, producing serpentine, magnetite and hydrogen. This reaction, so-called serpentinization, controls the physical and chemical properties of the oceanic lithosphere and is, thus, considered in geophysical models even if its kinetics are poorly constrained. Therefore, this study focuses on the kinetics of San Carlos olivine hydratation, a simple serpentinization reaction. An experimental method has been developed which allows monitoring magnetically the kinetics of the reaction of interest through the amount of magnetite that it produces. This method provides also an indirect estimate of the hydrogen production when coupled to thermochemical modeling. The influence of initial grain size and temperature on the hydration kinetics has also been investigated. An increase in the serpentinization rate at the beginning of the reaction is explained by the formation of microtextures also observed in natural samples (fractures and etch pits) which contribute to the generation of new reactive surface area. The inferred kinetics of serpentinization are by one to two orders of magnitude lower than the one commonly input in geophysical models. In addition to this experimental study, a natural case of RedOx reaction involving serpentinites has been studied which deals with a decimetre wide reaction zone between serpentinites and marbles from Alpine Corsica. The presence of wollastonite, CaSiO3, is explained on a petrological and thermodynamical basis as related to highly reducing conditions characterized by the presence of H2-rich fluids in this subduction zone environment. These fluids are interpreted as resulting from the conservation of the reducing potential of serpentinites throughout the subduction process.

Une nouvelle approche de modélisation des milieux géologiques fracturés représentés par un modèle conceptuel de fractures discrètes et déterministes est présentée dans cette thèse. L'objectif principal de l'étude est de reproduire l'hétérogénéité et la complexité des milieux poreux fracturés dans un modèle géométrique tridimensionnel afin d'effectuer des simulations numériques dans le but d'améliorer les capacités de modélisation en hydrogéologie. Ceci est réalisé à travers le couplage entre une plateforme de modélisation géologique (GOCAD) et un code numérique (HydroGeoSphere). Les principaux défis à relever sont: la représentation géométrique du réseau de fractures, la sélection d'un maillage approprié pour la discrétisation spatiale du domaine de simulation et l'adaptation du code numérique à ce maillage. La nouvelle approche est basée sur une première phase de modélisation géologique 3D, suivie par la génération d'un maillage tétraédrique 3D et par la simulation numérique de l'écoulement souterrain en conditions saturées et du transport de solutés. En général, le maillage tétraédrique s'avère plus adéquat que les maillages de blocs ou de prismes pour discrétiser les geometries complexes telles que les milieux fracturés. De plus, une définition alternative du maillage "dual", qui est essentiel pour appliquer la méthode numérique élément finis - volume de contrôle utilisée par HydroGeoSphere, est analysée et intégrée dans le code numérique. Le code numérique proposé est d'abord vérifié par l'intermédiaire de simples scénarios de simulation dont les solutions, analytiques et numériques, sont déjà connues. La complexité des simulations est augmentée de façon graduelle. L'approche de modélisation est finalement appliquée au site Olkiluoto (Finlande) où un laboratoire de recherche souterrain est en construction afin d'évaluer la faisabilité du stockage géologique profond de déchets nucléaires à haute activité. Les techniques de modélisation géologique mises au point permettent de modéliser facilement la géométrie des fractures identifiées à travers la caractérisation géologique in situ. De plus, le modèle numérique s'avère adéquat pour la simulation de l'écoulement et du transport de solutés dans ce site complexe. Ce travail de recherche présente une contribution au développement des techniques de modélisation hydrogéologique des milieux fracturés.

Cette étude porte sur les structures et déformations qui les ont engendrées des roches de la trilogie priabonienne du Champsaur. L'étude de terrain est axée sur les affleurements de schiste, roche déformable, particuliérement riche en petites structures engendrées par les différents évènements tectoniques et plus accessoirement sur les affleurements calcaires ou gréseux.

Les affleurements de Bevons, Nyons et Rosans dans le Sud-Est de la France, comme ceux du Numidien (Sicile, Tunisie, Maroc...) ou de la Tourelle (Canada) permettent l'observation de nombreuses injectites sableuses alimentées par des chenaux turbiditiques. Deux types d'injectites sont présents : les sills (horizontaux) et les dykes (verticaux), les dykes étant issus des sills. Leurs processus de mise en place sont, selon les auteurs, per ascensum, post-dépositionnels, et (sans doute le plus souvent) per descensum, contemporains de la mise en place du sable nourricier et objet de la présente étude. L'approche des mécanismes intervenant dans cette fracturation et la modélisation de celle-ci ont notamment pour but de préciser l'évolution de l'étanchéité des matériaux argileux pouvant constituer un site de stockage de déchets. D'un autre point de vue, des injections sableuses étant localisées entre des réservoirs contenant des hydrocarbures, une détermination de leur géométrie est nécessaire pour prédire les circulations potentielles de fluides durant la production. Un modèle géométrique a été établi à partir de ces observations et a permis de faire certaines hypothèses sur les mécanismes à prendre en compte dans les modélisations de cette fracturation et/ou injection. Plusieurs approches sont menées en parallèle : - L'étude de la compaction dans le massif, avant l'injection : pour des concentrations en sable assez élevées (> 35 %), le poids de la colonne de mélange d'injection est supérieure à la contrainte verticale dans la partie superficielle du massif. Cette inversion de densité peut être responsable du fait que l'injection se propage plutôt vers le bas dans le cas des dykes ou le long de discontinuités stratigraphiques (telles que les niveaux de cendres par exemple). Lorsque l'on descend dans les sédiments, la pression dans la colonne de mélange sableux varie linéairement avec la profondeur alors que la contrainte dans les sédiments varie non linéairement (issue d'une courbe de porosité en exponentielle). A partir d'une certaine profondeur, la pression fluide devient inférieure aux contraintes dans la succession sédimentaire : l'injection ne peut pas dépasser cette profondeur. Ainsi, la propagation des dykes est limitée par le poids de son encaissant. - L'utilisation de la fracturation hydraulique comme mécanisme de l'injection : la mécanique des roches et plus précisément les mécanismes d'hydrofracturation ont été testés pour tenter de modéliser ces observations à grande échelle. La force responsable de l'initiation et de la propagation des fractures est la pression due à la mise en place du chenal turbiditique. Le modèle tient compte du fait que l'injection doit soulever les terrains sus-jacents. Les premiers résultats ont permis de montrer que dès qu'une fracture horizontale est initiée, la propagation est rapide. L'ouverture obtenue (épaisseur d'un sill) est décamétrique, soit du même ordre de grandeur que les observations de terrain. - L'utilisation d'un fluide à seuil dans un réseau établi, avec des probabilités cohérentes avec la réalité des phénomènes entrant en jeu, a permis de déterminer les extensions maximales d'injection en accord avec les observations de terrain : - Un sill seul peut se propager sur 2400 m en quelques dizaines d'heures, - Un dyke seul peut descendre jusqu'à 400 m sous le paléofond de mer en quelques secondes, - Un réseau de sills et de dykes connectés (ce réseau étant connecté au chenal par un sill) peut avoir une extension horizontale de 1200 m et verticalement de 400 m, et se mettre en place en quelques heures. La formation des sills est plus longue que celle des dykes, la formation d'un sill seul dure plusieurs heures alors que celle d'un dyke dure quelques secondes. Lorsque les sills et les dykes sont connectés, le temps de propagation du réseau est d'environ 2 heures. Les dykes se forment toujours en quelques secondes et ont un effet sur la propagation des sills qui se propagent deux fois moins loin que s'ils étaient seuls.

Le mouvement des glaciers alpins est fortement conditionné par un grand nombre de processus complexes, alliant à la fois l'hydrologie, la mécanique, la thermodynamique, et dont la taille caractéristique peut varier de l'échelle microscopique à l'échelle macroscopique. Les mécanismes de fracturation peuvent être associés à des causes aussi nombreuses que diverses : frottement basal, ouverture de crevasses, variations de la température (fracturation thermique), variations de hauteur d'eau (fracturation hydraulique). Dans tous les cas, il a été demontré par le passé que la fracturation était génératrice de sismicité, enregistrable à l'aide de sismomètres. Cette sismicité permet par conséquent l'auscultation des conditions mécaniques d'un glacier, et en particulier son champ de contrainte.<br />Nous nous proposons, dans ce travail, d'appliquer des méthodes sismologiques inédites pour l'étude et la caractérisation de la fracturation dans les glaciers alpins, au travers de deux exemples (le glacier d'Argentière, dans le massif du Mont-Blanc, et le glacier du Gorner, dans le massif du Mont Rose, en Suisse). La localisation des sources a été effectuée par le biais de techniques d'antennes. Dans ce cadre, nous avons développé une méthode de recherche sur grille, bien adaptée à ces objets de taille limitée. Cette dernière permet le calcul de la densité de probabilité complète de la position de l'hypocentre. Elle nous a permis de mettre en évidence une activité sismique en essaim sur le glacier d'Argentière, et dont il a été possible de calculer une magnitude locale. Nous avons également montré que les séismes profonds ayant lieu au glacier du Gorner sont localisés à proximité immédiate de l'interface glace - roche. Parallèlement, d'autres techniques (corrélation de bruit, beamforming) nous ont permis de préciser l'information obtenue par la méthode de recherche exhaustive.

Le fonctionnement et la structure hydrogéologiques des aquifères de socle des granites et des gneiss sont aujourd'hui relativement bien compris. En revanche, dans cet ensemble, les aquifères dans les péridotites sont très peu étudiés et mal compris. Dans ce contexte, les massifs obductés de Nouvelle-Calédonie présentent un laboratoire naturel exceptionnel pour améliorer la connaissance de cet hydrosystème original en contexte tropical. Ainsi l'objectif du présent travail de thèse est d'approfondir la connaissance de la structure et du fonctionnement hydrogéologiques de ces massifs. L'analyse porte d’abord sur le manteau d'altération constitué de la cuirasse, des latérites, des saprolites grossières et du saprock. Plus de 60 essais hydrauliques sont menés sur les massifs de péridotites et les résultats sont compilés aux données existantes. La conductivité hydraulique moyenne des latérites est évaluée à 1.10-7 m/s et celle des saprolites grossières et du saprock à 8.10-7 m/s. L'hétérogénéité de cet horizon altéré est marquée par une gamme de variation de la conductivité hydraulique sur six ordres de grandeur et l'analyse piézométrique met en évidence des connexions hydrauliques avec le substratum fracturé profond. Le substratum est ensuite considéré. L'étude de la fracturation est réalisée à partir de mesures structurales sur affleurement et de la description de près de 1000 m de carottes de forages. L'analyse de la fracturation met en évidence l'importance du réseau serpentineux par sa densité d’une part, et par son lien avec l'altération supergène d'autre part. De plus, il est vérifié que la conductivité hydraulique du substratum diminue avec la profondeur. Cette variation est liée à la diminution de la densité de fractures altérées. Ainsi, à l'issue de ces analyses, la structure des massifs de péridotites est définie. Un réseau primaire de fractures d’espacement décimétrique lié au réseau serpentineux préstructure les péridotites. Sur ce réseau se surimpose un réseau de fractures dont l'espacement est décamétrique et caractérisé par une altération supergène. Les fractures altérées présentent localement de fortes conductivités hydrauliques, de l'ordre de 10-5 m/s. En profondeur l'espacement des fractures est hectométrique et les fractures sont majoritairement fermées, scellées par les minéraux néoformés ou par l'effet de la pression lithostatique. Les réseaux de fractures déca et hectométriques, visibles également sur l'effet d’échelle de la conductivité hydraulique, sont majoritairement verticaux, développés par instabilité de dissolution lors des processus d’altération. Cependant, des structures à faible pendage existent également et permettent la percolation du réseau. Enfin, à partir de ces nouveaux résultats et de l’intégration de l'ensemble des données acquises sur les différents massifs, un modèle de structure et de fonctionnement hydrogéologiques est proposé à l'échelle du massif. Ce modèle comprend l'horizon des latérites qui constitue un aquitard homogène sous lequel se développe l'aquifère dont l'épaisseur est de l'ordre de cinquante mètres. Le substratum est discrétisé en trois couches dont la conductivité hydraulique décroît de 2.10-7 à 2.10-8 m/s entre 50 m et 250 m environ sous le mur de l'aquifère. Les modèles numériques construits permettent de valider le modèle conceptuel unitaire et montrent que l'état de saturation des massifs est contraint par leur géomorphologie. Au terme de ce travail, plusieurs aspects doivent encore être approfondis. Le rôle hydrogéologique de la cuirasse doit être précisé et considéré dans le modèle hydrogéologique. Enfin, compte tenu du développement possible de structures très perméables, voire pseudo-karstiques, la connaissance de la distribution des structures drainantes doit être améliorée. Les résultats appliqués de ce travail de recherche sont exposés dans un rapport final et un guide méthodologique livrés dans le cadre du projet CNRT « HYPERK ». / Water resources of hard-rock (gneisses or granites) aquifers have significantly been studied in the past two decades. The hydrogeological behavior and structure of these aquifers are thus relatively well understood. On the other hand, aquifers in mantle-type basements, such as peridotites, are poorly studied and understood, mainly because they are not common and of limited extent. In this context, New Caledonia is a great laboratory offering unique opportunity to improve the knowledge of these original types of hydrosystems in tropical climate. Thus, the objective of this thesis is to improve the knowledge of these aquifer systems within weathered peridotites. Firstly, the study focuses on the characterization of the weathered layers of the peridotites composed of, from top to bottom, iron oxides/ferricrete, laterite, coarse saprolite and saprock (ie. top of the bedrock, with up to 20 % of weathered material). More than 60 hydraulic tests are performed and results were compiled with existing data. Mean hydraulic conductivity (K) of laterites is estimated around 1.10-7 m/s while mean value in coarse saprolites and saprock is around 8.10-7 m/s. Heterogeneity of this altered layer is high; K varies between six orders of magnitude and hydraulic head data analysis reveals a hydraulic connection with the deep fractured bedrock. Secondly, the fresh rock part of peridotites is studied. Fracture network analysis is derived from outcrop structural measurements and from the description of about 1000 m of cumulated borehole cores. This work highlights the importance of serpentine network, because of its high density and its critical impact on weathering. Moreover, the observations reveal that hydraulic conductivity decreases with depth within the substratum, due to the vertical decrease of weathered fractures density. These new results allow defining a structural framework of the massifs. It is characterized by a primary decimetrical fracture network closely related to the serpentine network. This network is overprinted by a secondary weathering network which reveals decametric spacing (ca. 30 m) and in places K values of 10-5 m/s. At depth, spacing is hectometrical and fractures are sealed by lithostatic pressure and/or subsequent mineral precipitations. These deca- and hectometric fracture networks, which are also visible on the scale effect of hydraulic conductivity, are primarily vertical and are the result of dissolution instabilities occurring during weathering processes. However, low- angle fractures do occur and allow the percolation of the network. Finally, on the basis of these new results and the integration of all existing data from different massifs a new hydrogeological conceptual model is proposed at the scale of a massif. The model includes a homogeneous lateritic aquitard and a coarse saprolite and saprock aquifer which is about 50 m thick. The bedrock is subdivided into three layers whose hydraulic conductivity decreases from 2.10-7 m/s to 2.10-8 m/s, respectively 50 and 250 m below the aquifer base. Numerical modelling validates this unitary conceptual model and reveals that the saturation of the massifs depends on their morphology. At last, several aspects require further research. The role of the ferricrete layer must be specified and considered in the hydrological model. Moreover, distribution of the fracture network remains to be fully addressed and should be studied with care given the potential development of highly permeable structures that could conform to pseudo-karstic drains. The applied results of this work are available in a “Technical guide” and a “Technical report” of the “HYPERK” CNRT Project.