Faculté des sciences

Modeling heat and fluid flow in discrete fracture networks : application to fluid-driven seismicity and engineered geothermal systems

Jansen, Gunnar ; Miller, Stephen A. (Dir.) ; Valley, Benoît (Codir.)

Thèse de doctorat : Université de Neuchâtel, 2018.

Understanding the dynamics of naturally fractured systems and fractured reservoirs in terms of flow, heat transport and fracture stability (e.g. earthquakes) is important for a range of applications including geothermal systems, waste water injection and natural fluid-driven seismicity. Many numerical codes exist that can simulate aspects of these dynamics, but in general require either expensive... More

Add to personal list
    Résumé
    La compréhension de la dynamique des systèmes naturels et des réservoirs fracturés en termes des écoulements, des transferts de la chaleur ou de stabilité de fractures (p. ex. les tremblements de terre) est importante pour diverses applications telles que les systèmes géothermiques et la séismicité naturelle induite par des fluides. Aujourd’hui, beaucoup de codes numériques existent pour simuler certains aspects de ces dynamiques, mais en général, ils sont très couteux ou alors leur utilisation n’est possible que par leurs développeurs eux-mêmes. Dans cette thèse, je présente une solution open-source MATLAB efficace pour des simulations numériques de processus couplés dans les systèmes fracturés. J’utilise le modèle des fractures discrètes intégrées.
    Dans une étude numérique thermoélastique, je cherche à comprendre l’effet de l’injection de fluides froids dans un réservoir plus chaud et plus spécifiquement les changements des contraintes associées aux variations de température. Ces changements induisent des ruptures en cisaillement dans le réservoir. J’ai observé une influence très forte des propriétés hydrauliques du réservoir sur l’évolution des contraintes thermiques. De plus, j’ai trouvé que le changement des contraintes thermiques peut conduire à du cisaillement sur des fractures orientées non optimalement pour la rupture. Ces résultats suggèrent que les changements des contraintes thermiques devraient être pris en compte dans tous les modèles lors d’injections de fluides de longue durée dans des réservoirs fracturés.
    J’utilise le modèle développé pour comprendre une séquence de tremblements de terre qui prit place à l’Ouest de Reno, Mogul, Nevada, aux États-Unis à la fin février 2008. Après 2 mois de séismes précurseurs répétés près d’une structure de faille inconnue préalablement, un tremblement de terre de magnitude 4.9 a eu lieu. Ici, je montre que la séquence de séismes précurseurs a pu être engendrée par l’intrusion de fluides à haute pression dans la structure de faille préexistante et j’ai trouvé une forte corrélation entre le front de pression et l’hypocentre des événements sismiques.
    Finalement, je présente un algorithme simple et efficace pour générer des maillages hexahedraux à partir d’un modèle géologique utilisable par n’importe quels simulateurs numériques. En utilisant la méthode de «binary space partitioning» et d’«octree refinement» sur la géométrie importée, un maillage précis est créé. L’algorithme fournit une nouvelle méthode pour la création de maillages hexahedraux pour toutes configurations géologiques.
    Summary
    Understanding the dynamics of naturally fractured systems and fractured reservoirs in terms of flow, heat transport and fracture stability (e.g. earthquakes) is important for a range of applications including geothermal systems, waste water injection and natural fluid-driven seismicity. Many numerical codes exist that can simulate aspects of these dynamics, but in general require either expensive licences or their utility is limited mostly to the developers. In this thesis I present an open source MATLAB package for efficient numerical simulations of the coupled processes in fractured systems. I take advantage of the embedded discrete fracture model that accounts for discrete fractures that are computationally efficient.
    In a numerical thermo-elasticity study, I investigate the effect of cold fluid injection on the reservoir and the resulting thermal stress change on potential shear failure in the reservoir. I observe a strong influence of the hydraulic reservoir properties on thermal stress propagation. I also find that thermal stress change can lead to induced shear failure on non-optimally oriented fractures. The results suggest that thermal stress changes should be taken into account in all models for long-term fluid injections in fractured reservoirs.
    Furthermore, I use the developed model to investigate an earthquake sequence that occurred in Mogul west of Reno, Nevada, USA in late February 2008. It culminated in a magnitude 4.9 main-shock after a foreshock-rich period of approximately 2 months on a previously unidentified fault structure. Here I show that the foreshock sequence may have been driven by a fluid pressure intrusion into this preexisting structure and find a strong correlation between high fluid pressure fronts and foreshock hypocenters.
    Finally, I present a simple and efficient algorithm to generate hexahedral meshes from a geological model to be used in numerical simulation tools. Using binary space partitioning of the input geometry and octree refinement on the grid, a successive increase in accuracy of the mesh is achieved. The algorithm provides a new method for hexahedral mesh generation in geological settings. It generates high accuracy discretizations with cell counts suitable for state-of-the-art subsurface simulators.