Fondations du transport réactif dans les systèmes géologiques

Le transport réactif dans les environnements souterrains implique le déplacement simultané des fluides et les transformations chimiques qui se produisent lorsqu'ils interagissent avec les minéraux solides, la matière organique et d'autres fluides.Cette interaction régit une large gamme de processus naturels et techniques, y compris la migration des contaminants des eaux souterraines, la formation de dépôts minéraux et la réponse des réservoirs à l'injection de fluides.

La modélisation thermodynamique fournit la base essentielle pour comprendre quelles réactions sont possibles, jusqu'où elles vont se dérouler et quelles phases minérales apparaîtront ou disparaîtront au fil du temps. Sans cette base, les simulations de transport réactives reposeraient sur des hypothèses arbitraires ou des ensembles de données incomplets, conduisant à des prédictions peu fiables.

Principaux moteurs du transport sous-marin réactif

Mécanismes de transport en flux fluidique et en masse

Le mouvement des solutions aqueuses à travers des milieux poreux et fracturés est régi par l'advection, la dispersion et la diffusion moléculaire. L'avis porte les espèces dissoutes à la vitesse du fluide en vrac, tandis que la dispersion les propage en raison de l'hétérogénéité des voies d'écoulement. La diffusion devient significative dans les zones à faible perméabilité ou les régions stagnantes.

Dans la pratique, les modèles de transport réactifs résolvent les équations de conservation pour la masse, l'élan et l'énergie, souvent en utilisant la loi Darcy’s pour le débit. Le couplage entre le débit et la chimie se produit parce que les réactions peuvent modifier la porosité et la perméabilité par dissolution minérale ou précipitation, se réalimentant ainsi dans le champ de débit.

Réactions chimiques : Équilibre versus cinétique

La modélisation thermodynamique suppose généralement que certaines réactions, comme la complexation aqueuse, les équilibres acides-bases et l'échange d'ions, et #8211; atteignent un équilibre instantané par rapport aux échelles de temps de transport. Cette hypothèse simplifie le système en réduisant le nombre d'équations différentielles nécessaires. Cependant, la dissolution minérale et les précipitations sont souvent limitées cinétiquement, surtout à basse température ou lorsque des espèces inhibitrices sont présentes.

Pour gérer les deux régimes, les modèles modernes combinent des calculs d'équilibre thermodynamique pour des réactions rapides avec des lois de vitesse cinétiques pour les plus lentes. La loi de vitesse elle-même dépend de quantités thermodynamiques telles que l'indice de saturation, qui mesure la distance entre la solution et l'équilibre par rapport à un minéral donné.

Thermodynamique : la base énergétique

Constantes d'équilibre et loi de l'action de masse

Chaque réaction chimique est caractérisée par une constante d'équilibre, K, qui est dérivée du changement d'énergie libre standard de Gibbs, ΔG° = -RT ln K[. Pour les réactions dans des milieux poreux, la constante d'équilibre relie les activités des produits et des réactifs lorsque le système est au repos.

Dans une simulation de transport réactif, l'hypothèse d'équilibre local invoque ces constantes à chaque cellule ou nœud, en résolvant pour la spéciation qui minimise l'énergie libre globale Gibbs du système. Cette approche est efficace par calcul mais suppose que la phase du fluide est bien mélangée et que les vitesses de réaction sont assez rapides pour maintenir l'équilibre localement.

Modèles d'activité et comportement non idéal

Les eaux naturelles sont rarement diluées; elles contiennent des concentrations élevées d'ions dissous qui interagissent électrostatiquement. Le coefficient d'activité γ corrige pour ces interactions non-idéales, convertissant les concentrations en activités thermodynamiques significatives. Plusieurs modèles d'activité sont utilisés, chacun avec sa propre gamme d'applicabilité:

  • Théomatique Debye-Hückel – valable pour une faible résistance ionique (typiquement < 0,1 M).
  • Extended Debye-Hückel (p. ex., équation de Davies) – adapté jusqu'à ~0,5 M.
  • Modèle d'interaction spécifique aux ions du pitzer – précis des salinités diluées aux salinités saumâtres.
  • Modèle de la HKF – utilisé pour les espèces aqueuses à haute température et pression.

Le choix du modèle d'activité correct est essentiel car les erreurs de coefficients d'activité se propagent directement dans les indices de saturation et les prévisions d'équilibre. Une saumure dont les concentrations en Na+ et Cl- dépassent 4 M, par exemple, nécessite une approche Pitzer pour éviter de grandes inexactitudes.

Solubilité minérale et diagrammes de phase

La solubilité d'un minéral est définie par son produit de solubilité Ksp. Lorsque le produit d'activité ionique (PAI) dépasse Ksp, la solution est sursaturée et la précipitation est thermodynamiquement favorisée; lorsque IAP < K]sp], la dissolution se produit. L'indice de saturation et#937; = IAP / ]K]sp [ou log Q/K est un diagnostic central dans tous les codes de transport

Les diagrammes de stabilité, comme les diagrammes de Eh-pH (Pourbaix) pour les éléments sensibles au redox ou les diagrammes d'activité pour les minéraux argileux, aident à visualiser les conditions dans lesquelles les différentes phases sont stables. Ces diagrammes sont dérivés de données thermodynamiques et permettent aux modélistes d'identifier rapidement les assemblages minéraux dominants pour une chimie de l'eau donnée.

Approches numériques pour la modélisation couplée

Plateformes logicielles de pointe

Différents codes ont été élaborés pour résoudre le système étroitement couplé d'équations de débit, de transport et de réaction, notamment :

  • PHREEQC – un code géochimique polyvalent de l'USGS qui peut être utilisé comme moteur de spéciation autonome ou relié à des simulateurs de transport.
  • Geochemist’s Workbench (GWB) – une suite commerciale offrant un transport réactif en 1D et 2D, avec des bases de données thermodynamiques robustes et une modélisation réaction-chemin.
  • TOUGHREACT – un code de transport réactif multicomposants non isotherme élaboré au Laboratoire national Lawrence Berkeley, particulièrement adapté aux problèmes de géothermie et de séquestration du CO2.
  • OpenFOAM avec des résolveurs personnalisés – une plate-forme de dynamique de fluide informatique open source qui peut être étendue avec des modules de réaction chimique.

Chaque code a des forces : PHREEQC excelle dans la géochimie par lots et le couplage avec le simple transport ; TOUGHREACT gère le flux multiphasé à grande échelle ; GWB fournit un pré-traitement et un post-traitement convivial. Le choix dépend de l'échelle de problème, dimensionnalité et complexité de la thermodynamique requise.

Bases de données thermodynamiques et leur qualité

La précision de tout modèle thermodynamique dépend de la base de données sous-jacente. Les principales bases de données de référence comprennent thermo.dat (utilisé par PHREEQC), llnl.dat (Lawrence Livermore National Laboratory), wateq4f, et MINTEQ[. Ces bases de données contiennent des valeurs log K, des données enthalpie et des paramètres de modèles d'activité pour les minéraux, les gaz et les espèces aqueuses.

Les efforts récents ont porté sur la cohérence entre les bases de données et l'élargissement de la couverture aux températures et pressions élevées applicables aux dépôts géologiques profonds et aux réservoirs géothermiques. La base de données ThermoChimie, par exemple, a été conçue spécifiquement pour les applications d'élimination des déchets nucléaires et comprend de nombreuses données sur la solubilité des radionucléides.

Les utilisateurs doivent être conscients des incertitudes dans les données thermodynamiques, en particulier pour les oligo-éléments ou les minéraux avec un comportement complexe de résolution solide.

Stratégies de solution séquentielle et entièrement couplée

Les codes de transport réactifs résolvent le couplage des processus physiques et chimiques en utilisant soit une approche de partage de l'opérateur (séquentiel), soit une méthode entièrement couplée (globale implicite).

  • Approche séquentiel non-itérative (SNIA) ou itérative (SIA) : Le débit et le transport sont d'abord résolus, puis des calculs de réaction sont effectués pour chaque cellule en utilisant les concentrations mises à jour. Ceci est moins cher sur le plan calcul, mais peut introduire des erreurs de bilan massique si le temps est trop long ou si les réactions sont rapides.
  • Méthode implicite globale: Toutes les équations qui régissent (droits de débit, de transport, d'action de masse et de vitesse cinétique) sont assemblées en un seul système d'équations non linéaires et résolues simultanément.

Pour de nombreuses applications pratiques, l'approche itérative séquentielle permet de trouver un équilibre entre la précision et le coût, surtout lorsque les délais de réaction sont comparables aux délais de transport.

Applications critiques dans les domaines des sciences de la terre et de l'ingénierie

Remise en état des eaux souterraines

Les aquifères contaminés exigent souvent des stratégies d'assainissement in situ qui reposent sur la manipulation des conditions chimiques pour immobiliser les polluants. Par exemple, des barrières réactives perméables contenant du fer zéro-valent peuvent réduire les solvants chlorés. Les modèles thermodynamiques prédisent la séquence des produits de corrosion du fer (ferrihydrite, magnétite, rouille verte) qui se forment, leur réactivité avec les contaminants et l'évolution à long terme de la conductivité hydraulique.

Récupération améliorée du pétrole (ROE)

Dans les réservoirs de pétrole et de gaz, l'inondation de l'eau et les inondations chimiques modifient la composition ionique des brines de formation, ce qui modifie la capacité d'humidité minérale et peut précipiter les minéraux formant des échelles comme la calcite ou la barite. Les modèles de transport réactif thermodynamique aident les ingénieurs à concevoir des fluides d'injection qui réduisent les dommages à la formation et optimisent le déplacement du pétrole.

Stockage géologique du carbone

L'injection de CO2 dans des aquifères salins profonds ou dans des champs d'huile appauvris déclenche une cascade de réactions acides-basiques et minérales. Le CO2 dissous forme de l'acide carbonique, abaissant le pH et dissolvant les minéraux carbonés. Cela augmente d'abord la porosité, mais à plus long terme, le rejet de cations peut entraîner des précipitations de carbonates secondaires tels que la dawsonite ou la sidérite, qui piègent en permanence le CO2. Des simulations de transport réactifs sont utilisées pour évaluer la capacité de stockage, l'intégrité du confinement et le potentiel de dissolution ou d'autoscellement du caprock.

Élimination des déchets nucléaires

Les modèles thermodynamiques traitent la dissolution des déchets vitrifiés, la formation de phases secondaires (p. ex., minéraux argileux, oxydes de fer, silicates d'uranyl) et la spéciation des radionucléides dans les eaux souterraines. Comme de nombreux radionucléides (p. ex., Tc-99, Np-237) ont une très faible solubilité dans des conditions de réduction, la limite de solubilité thermodynamique constitue une barrière essentielle au rejet. Les modèles doivent également tenir compte de l'influence de la température de la décomposition radioactive et de l'évolution des gradients chimiques après la fermeture du dépôt.

Défis actuels et frontières

Incertitude et lacunes dans les données

Les constantes thermodynamiques de nombreux minéraux, en particulier les phases métastables et les solutions solides, sont mal limitées. Par exemple, les argiles et les zéolites présentent des compositions variables, ce qui rend difficile la définition d'un produit de solubilité unique. De plus, l'extrapolation des données thermodynamiques à des pressions élevées (au-delà de 1000 bar) et des températures (au-dessus de 300 °C) nécessite une validation attentive par rapport aux mesures expérimentales.

Scalabilité informatique

Les simulations de transport réactif à l'échelle du réservoir avec des domaines tridimensionnels et une discrétisation spatiale fine peuvent impliquer des millions de cellules, chacune avec des centaines d'espèces chimiques et de réactions.Les méthodes implicites entièrement couplées deviennent prohibitivement coûteuses.Les chercheurs développent des maillages adaptatifs, une réduction du modèle (p. ex., décomposition orthogonale adéquate) et des émulateurs d'apprentissage de machine qui approximent le sous-système chimique tout en conservant la cohérence thermodynamique essentielle.

Intégration de la biogéochimie

L'activité microbienne peut accélérer ou inhiber les réactions géochimiques par des transformations de redox, la formation de biofilms et la production de ligands organiques. L'incorporation de cinétique microbienne et d'énergie thermodynamique (p. ex., Gibbs, rendements énergétiques libres pour les voies métaboliques) est une frontière émergente. Par exemple, la réduction du sulfate en sulfure par les bactéries réducteurs de sulfates peut déclencher la précipitation de sulfures métalliques, en éliminant les métaux de la solution.

Perspectives d'avenir : le rôle de l'assimilation des données et de l'apprentissage automatique

À mesure que les modèles thermodynamiques s'intègrent davantage aux données de surveillance sur le terrain, les techniques d'assimilation des données telles que le filtrage de l'ensemble Kalman permettent de mettre à jour les paramètres du modèle (y compris les constantes thermodynamiques) en temps réel au fur et à mesure que de nouvelles observations arrivent, ce qui réduit l'incertitude prédictive et améliore la fiabilité des simulations utilisées pour les décisions à haute prise.

Les réseaux neuraux peuvent être formés à de grands ensembles de calculs thermodynamiques pour servir de substituts rapides au moteur géochimique dans un code de transport réactif. Bien que ces modèles reposent toujours sur des données thermodynamiques de haute qualité pour l'entraînement, ils accélèrent considérablement les simulations d'ensembles et les analyses de sensibilité. Le défi reste de veiller à ce que ces substituts fondés sur les données respectent les lois thermodynamiques et ne produisent pas d'extrapolations non physiques.

La modélisation thermodynamique du transport réactif dans les formations géologiques est une discipline mature mais en évolution rapide. En fournissant les contraintes énergétiques et d'équilibre qui régissent les interactions minéralo-fluidiques, la thermodynamique reste le cadre essentiel pour prédire le comportement subsurface à travers les échelles de temps des heures à des millénaires.