Mathematical Modeling of Salt Precipitation and Multi-Phase Flow in High Enthalpy Fractured Geothermal Systems

Cet article présente un nouveau modèle de flux compositionnel en code source ouvert, implémenté dans le cadre PorePy, qui simule l'écoulement multiphasique et non isotherme ainsi que la précipitation de la halite dans les réservoirs géothermiques fracturés à haute enthalpie, en utilisant une formulation robuste des variables primaires et une approche matrice-fracture discrète pour prédire avec précision les dommages de perméabilité et les défis opérationnels.

L'essentiel

Imaginez un réservoir souterrain à haute température comme une gigantesque cocotte-minute naturelle remplie d'eau très chaude et salée. C'est la source de l'énergie géothermique. Pour extraire cette énergie, les ingénieurs forent des puits et pompent de l'eau froide pour repousser l'eau chaude vers le haut. Cependant, ce processus est délicat à cause du « sel » présent dans l'eau.

Considérez le sel (halite) comme le sucre dans une tasse de thé chaud. Si vous refroidissez le thé ou si une partie de l'eau s'évapore, le sucre ne peut plus rester dissous et commence à redevenir des cristaux solides. Dans un puits géothermique, cela se produit lorsque l'eau chaude se refroidit près d'un puits d'injection ou s'évapore près d'un puits de production. Le résultat ? Des cristaux de sel se forment et bouchent les minuscules trous dans les roches et les fissures qui permettent à l'eau de circuler. C'est comme un embouteillage causé par des cristaux de sucre bloquant une autoroute.

Cet article présente un nouvel outil de simulation informatique conçu pour prédire exactement où et comment ces « embouteillages de sucre » se produiront dans des roches souterraines complexes et fissurées.

Voici une décomposition du fonctionnement de l'outil et de ce qu'il a découvert, en utilisant des analogies simples :

1. La Carte : Voir les fissures clairement

Les cartes traditionnelles des roches souterraines lissent souvent tout, traitant la roche comme une éponge solide. Mais en réalité, l'eau circule principalement à travers un réseau de fissures, comme l'eau circulant à travers un trottoir fissuré plutôt qu'à travers le béton lui-même.

• L'Innovation : Ce nouveau modèle utilise une approche de « Matrice à Fractures Discrètes ». Imaginez dessiner les fissures comme des lignes fines et distinctes sur une carte, plutôt que de simplement les fondre dans l'arrière-plan. Cela permet à l'ordinateur de voir exactement comment les fissures se connectent (ou ne se connectent pas) et comment le sel pourrait boucher une fissure spécifique par rapport à la roche environnante.

2. Le Moteur : Une « télécommande universelle » pour la physique

Simuler l'eau bouillante, la vapeur et le sel solide simultanément est incroyablement difficile pour les ordinateurs. Habituellement, l'ordinateur doit constamment changer de « mode » (par exemple, « D'accord, maintenant c'est un liquide ; maintenant c'est un gaz ; maintenant c'est un solide »), ce qui peut provoquer un plantage du calcul ou un blocage.

• L'Innovation : Les auteurs ont créé un système « unifié ». Imaginez une télécommande universelle qui fonctionne pour tous les appareils sans avoir besoin de changer de piles ou de mode. Le modèle utilise trois « cadrans » fixes (Pression, Énergie Thermique et Quantité de Sel) qui restent les mêmes, que l'eau soit liquide, vapeur ou en train de devenir du sel solide. Cela rend la simulation beaucoup plus fluide et stable, lui permettant de gérer les changements d'états chaotiques sans s'interrompre.

3. L'astuce de vitesse : La « fiche de révision »

Calculer la physique exacte de l'eau salée à haute température nécessite généralement que l'ordinateur résolve des puzzles mathématiques complexes encore et encore, ce qui est très lent.

• L'Innovation : L'équipe a créé une « fiche de révision » pré-calculée (une table de correspondance). Avant que la simulation ne commence, ils ont calculé tous les résultats possibles de la façon dont le sel se comporte sous différentes conditions et les ont stockés. Pendant la simulation, au lieu de résoudre les mathématiques difficiles à chaque fois, l'ordinateur consulte simplement la réponse sur la fiche. Cela rend la simulation beaucoup plus rapide tout en restant précise.

4. L'effet d'obstruction : Le « rétrécissement des pores »

Lorsque les cristaux de sel se forment, ils prennent de la place.

• L'Innovation : Le modèle réduit automatiquement la taille des « tuyaux » (porosité et perméabilité) à mesure que le sel s'accumule. Il utilise une règle (Kozeny-Carman) qui stipule : « Si le sel remplit 10 % du trou, le tuyau devient considérablement plus étroit. » Cela permet au modèle de prédire comment le débit ralentira ou s'arrêtera complètement à mesure que l'« embouteillage de sucre » s'aggrave.

Ce que les simulations ont montré

L'équipe a testé cet outil dans deux scénarios principaux :

Scénario A : L'autoroute brisée (Fissures déconnectées)

• Configuration : Imaginez un réservoir où les fissures ne sont pas connectées entre elles ; l'eau doit se frayer un chemin à travers la roche solide entre elles.
• Résultat : Lorsque l'on a pompé de l'eau froide, l'eau chaude près du puits de production a bouilli rapidement. Cela a provoqué la cristallisation du sel et le bouchage de la roche juste autour du puits.
• Le rebondissement : Si l'on pompait l'eau plus rapidement, l'obstruction s'aggravait considérablement et la production d'énergie chutait de manière significative. Le modèle a montré que l'« embouteillage » se produisait principalement dans la roche près du puits, et pas seulement dans les fissures.

Scénario B : L'autoroute connectée (Fissures connectées)

• Configuration : Imaginez un réservoir où les fissures forment une autoroute continue à grande vitesse de l'injection vers le puits de production.
• Résultat : L'eau froide a circulé rapidement à travers les fissures. Comme elle se déplaçait très vite et restait fraîche, elle a en fait dissous le sel près du puits de production au lieu de le boucher !
• Le rebondissement : La précipitation du sel s'est déplacée vers un autre endroit — juste à la limite de la zone d'eau froide — plutôt que de boucher le puits lui-même. Cela suggère qu'un réseau de fissures connecté pourrait en fait protéger le puits contre l'obstruction, même si cela change l'endroit où le sel s'accumule.

L'essentiel

Ce document présente un nouvel outil logiciel en libre accès (open-source) qui aide les ingénieurs à comprendre la danse complexe entre la chaleur, la pression et le sel dans les puits géothermiques. En cartographiant avec précision comment les fissures se connectent et comment le sel les bouche, l'outil peut aider à prédire :

1. Où les puits pourraient être bloqués par le sel.
2. Quelle quantité d'énergie peut être extraite en toute sécurité avant que les « tuyaux » ne soient obstrués.
3. Si la disposition des fissures souterraines aidera ou nuira au processus de production.

Les auteurs ont vérifié leur outil par rapport à un standard établi de l'industrie et ont constaté qu'il correspondait parfaitement, prouant qu'il s'agit d'un moyen fiable de simuler ces environnements de roches fissurées et salines à haute température.

Résumé technique

Résumé Technique : Modèle de Flux Compositionnel Unifié pour la Simulation de Systèmes Géothermiques à Haute Enthalpie et Fracturés

Énoncé du Problème
La simulation des réservoirs géothermiques fracturés à haute enthalpie présente des défis importants en raison du couplage complexe des flux multiphasiques et multicomposants non isothermes, d'une thermodynamique fortement non linéaire et de l'influence dominante des réseaux de fractures. Un problème opérationnel critique dans ces systèmes est la mise en solution minérale, spécifiquement la précipitation de la halite (NaCl). Pendant la production, l'abaissement de la pression induit l'ébullition, concentrant la saumure résiduelle et dépassant les limites de solubilité, tandis que l'injection d'eau froide abaisse la solubilité locale. Ces deux mécanismes conduisent à la précipitation de sel dans les pores et les fractures, dégradant gravement la porosité et la perméabilité, réduisant la productivité des puits et causant potentiellement le colmatage des puits. Les modèles continus existants échouent souvent à capturer la précipitation localisée et l'évolution de la perméabilité à l'échelle des fractures, tandis que les simulateurs compositionnels traditionnels reposant sur des formulations de changement de variables (variable-switching) peinent avec la stabilité numérique près des frontières de phase.

Méthodologie
Les auteurs présentent un nouveau modèle de flux compositionnel unifié implémenté au sein du framework open-source PorePy. La méthodologie intègre trois composantes fondamentales :

1. Formulation à Variables Persistantes : Contrairement aux approches traditionnelles qui changent les variables primaires en fonction de la présence de phases, ce modèle utilise un ensemble fixe de variables primaires — pression ($p$), enthalpie spécifique ($h$) et fraction massique globale de sel ($z$) — à travers toutes les configurations de phase (liquide, vapeur et halite solide). Cette approche gère naturellement les transitions de phase sans changement manuel, améliorant la robustesse et la stabilité numériques.
2. Approche de Fracture Discrète-Matrice (DFM) : Le domaine physique est représenté par une géométrie de dimension mixte où les fractures et leurs intersections sont traitées comme des variétés de dimensions inférieures (1D et 0D) noyées dans la matrice poreuse 2D. Cette représentation explicite préserve la connectivité et l'hétérogénéité des fractures, permettant de résoudre des processus localisés tels que la réduction de la perméabilité induite par la précipitation.
3. Clôture Thermodynamique et Linéarisation : Le modèle emploie les corrélations complètes de Driesner et Christoph (2007) pour la thermodynamique H2O–NaCl, couvrant les conditions magmatiques. Pour réduire le coût computationnel, le modèle utilise une stratégie de Linéarisation Basée sur l'Opérateur (OBL). Les propriétés thermodynamiques et leurs dérivées sont précalculées sur une grille structurée dans l'espace des variables primaires et récupérées via une interpolation multilinéaire pendant la simulation, éliminant ainsi les calculs de séparation de phase (flash) coûteux effectués à la volée.
4. Évolution Dynamique de la Perméabilité et de l'Ouverture : Le modèle incorpore des relations de type Kozeny-Carman pour mettre à jour dynamiquement la porosité et la perméabilité de la matrice en fonction de la saturation en halite. Pour les fractures, l'ouverture évolue selon une relation de loi de puissance dépendant de la saturation en halite, ce qui met à jour la perméabilité de la fracture via la loi cubique.

Les équations directrices (bilans de masse et d'énergie) sont discrétisées à l'aide d'un schéma de volume fini centré sur les cellules avec une approximation de flux multipoints (MPFA) pour les termes diffusifs et un upwinding pour les termes advectifs. Le système non linéaire résultant est résolu à l'aide de la méthode de Newton avec une recherche de ligne Armijo de repli (backtracking) et un pas de temps adaptatif.

Résultats Clés
Le modèle a été vérifié et appliqué à travers plusieurs expériences numériques :

• Vérification (Benchmark 1D) : Le modèle a été vérifié par rapport au simulateur propriétaire CSMP++ en utilisant un benchmark de dissolution de sel en 1D. Les résultats ont montré une forte concordance à travers quatre régions de phases distinctes (vapeur+halite, vapeur+liquide, phase liquide unique, et liquide+halite) ainsi que leurs transitions, validant la formulation à variables persistantes et l'implémentation thermodynamique.
• Réseau de Fractures Déconnecté (2D) : Dans un réservoir avec des fractures isolées, l'écoulement est dominé par la matrice. La simulation a démontré que l'abaissement de la pression au puits de production induit l'ébullition, entraînant une précipitation significative de halite dans la matrice environnante et les segments de fractures adjacents. Cette précipitation a réduit la perméabilité de la matrice jusqu'à 36 % ($K/K_0 \approx 0,64$) et a altéré les lignes de courant. L'augmentation du taux d'injection a intensifié l'ébullition près du puits, conduisant à une saturation en halite plus élevée ($s_{hal} \approx 0,46$) et à une réduction de perméabilité plus sévère ($K/K_0 \approx 0,29$), impactant significativement les taux de production.
• Réseau de Fractures Connecté (2D) : Dans un réservoir présentant une chaîne de fractures continue reliant les puits d'injection et de production, le transport est devenu dominé par les fractures. Le réseau connecté a facilité une communication de pression rapide et le transport de fluide à faible salinité vers le puits de production. Par conséquent, bien que l'ébullition se soit produite près du producteur, le fluide arrivant est resté suffisamment dilué pour empêcher la saturation locale en halite. Au lieu d'une précipitation au niveau du producteur, une dissolution a été observée, et la précipitation a été redistribuée dans la matrice devant le front d'injection et le long de la chaîne de fractures en raison des effets de refroidissement. Cette configuration a permis d'éviter la grave dégradation de la perméabilité observée dans le cas déconnecté.

Signification et Revendications
L'article affirme fournir un outil numérique robuste et open-source pour analyser les transferts complexes de chaleur et de masse avec des effets de mise en solution minérale dans les systèmes géothermiques fracturés à haute enthalpie. Sa principale importance réside dans :

• Stabilité Numérique : La formulation unifiée à variables persistantes évite les difficultés de convergence associées au changement de variables près des frontières de phase.
• Fidélité Géométrique : L'approche DFM résout explicitement la connectivité des fractures, ce que les auteurs démontrent être un facteur critique contrôlant la distribution spatiale de la précipitation et la réduction de transmissivité résultante.
• Aperçu Opérationnel : Les résultats soulignent que la connectivité des fractures dicte les défis opérationnels ; les réseaux déconnectés sont sujets au colmatage près du puits et à la perte de productivité, tandis que les réseaux connectés peuvent atténuer ces problèmes en maintenant la dilution du fluide au puits de production.
• Reproductibilité : Le travail fournit le modèle numérique complet et le flux d'exécution via un conteneur Docker, facilitant la reproductibilité et le développement ultérieur au sein de la communauté géothermique.

Les auteurs concluent que leur modèle prédit efficacement les motifs de précipitation de la halite et leur impact sur la performance du réservoir, offrant un outil précieux pour comprendre l'interaction entre l'écoulement, la thermodynamique et la mise en solution minérale dans les réservoirs géothermiques fracturés à haute enthalpie.