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

Questo articolo presenta un nuovo modello di flusso composizionale open-source implementato nel framework PorePy che simula il flusso multifase non isotermico e la precipitazione di alite in serbatoi geotermici fratturati ad alta entalpia, utilizzando una formulazione robusta delle variabili primarie e un approccio matrice-frattura discreta per prevedere accuratamente il danno di permeabilità e le sfide operative.

L'essenziale

Immaginate un serbatoio sotterraneo ad alta temperatura come una gigantesca pentola a pressione naturale piena di acqua salata caldissima. Questa è la fonte dell'energia geotermica. Per estrarre quell'energia, gli ingegneri scavano pozzi e pompano acqua fredda all'interno per spingere l'acqua calda verso l'alto. Tuttavia, questo processo è complicato a causa del "sale" presente nell'acqua.

Pensate al sale (alite) come allo zucchero in una tazza di tè caldo. Se raffreddate il tè o lasciate che parte dell'acqua evapori, lo zucchero non può più rimanere sciolto e inizia a trasformarsi nuovamente in cristalli solidi. In un pozzo geotermico, questo accade quando l'acqua calda si raffredda vicino a un pozzo di iniezione o evapora vicino a un pozzo di produzione. Il risultato? Cristalli di sale solidi si formano e intasano i piccoli fori nelle rocce e le crepe (fratture) che permettono il passaggio dell'acqua. È come un ingorgo causato da cristalli di zucchero che bloccano un'autostrada.

Questo articolo presenta un nuovo strumento di simulazione al computer progettato per prevedere esattamente dove e come questi "ingorghi di zucchero" avverranno in complessi scogli sotterranei fratturati.

Ecco una ripartizione di come funziona lo strumento e di ciò che ha scoperto, utilizzando analogie semplici:

1. La Mappa: Vedere le Crepe chiaramente

Le mappe tradizionali delle rocce sotterranee spesso levigano tutto, trattando la roccia come una spugna solida. Ma nella realtà, l'acqua scorre principalmente attraverso una rete di crepe, come l'acqua che scorre attraverso un marciapiede crepato piuttosto che attraverso il cemento stesso.

• L'Innovazione: Questo nuovo modello utilizza un approccio "Discrete Fracture-Matrix" (Matrice a Fratture Discrete). Immaginate di disegnare le crepe come linee sottili e distinte su una mappa, invece di limitarsi a sfumarle nello sfondo. Ciò consente al computer di vedere esattamente come le crepe si collegano (o non si collegano) e come il sale potrebbe intasare una specifica crepa rispetto alla roccia circostante.

2. Il Motore: Un "Telecomando Universale" per la Fisica

Simulare acqua che bolle, vapore e sale solido tutto insieme è incredibilmente difficile per i computer. Di solito, il computer deve cambiare costantemente la sua "modalità" (ad esempio, "Ok, ora è liquido; ora è gas; ora è solido"), il che può causare il crash o il blocco dei calcoli.

• L'Innovazione: Gli autori hanno creato un sistema "unificato". Pensate a un telecomando universale che funziona per ogni dispositivo senza dover cambiare batterie o modalità. Il modello utilizza tre "manopole" fisse (Pressione, Energia Termica e Quantità di Sale) che rimangono invariate indipendentemente dal fatto che l'acqua sia liquida, vapore o si stia trasformando in sale solido. Questo rende la simulazione molto più fluida e stabile, permettendole di gestire il caotico passaggio tra gli stati senza interrompersi.

3. Il Trucco della Velocità: Il "Foglietto di Ripasso"

Calcolare l'esatta fisica dell'acqua salata ad alte temperature richiede solitamente che il computer risolva complessi enigmi matematici ripetutamente, il che è molto lento.

• L'Innovazione: Il team ha creato un "foglietto di ripasso" pre-calcolato (una tabella di consultazione). Prima dell'inizio della simulazione, hanno calcolato tutti i possibili risultati su come il sale si comporta in diverse condizioni e li hanno memorizzati. Durante la simulazione, invece di risolvere il difficile calcolo matematico ogni volta, il computer consulta semplicemente le risposte sul foglietto. Questo rende la simulazione molto più veloce pur mantenendo l'accuratezza.

4. L'Effetto di Intasamento: Il "Rimpicciolimento dei Pori"

Quando si formano i cristalli di sale, essi occupano spazio.

• L'Innovazione: Il modello restringe automaticamente i "tubi" (porosità e permeabilità) man mano che il sale si accumula. Utilizza una regola (Kozeny-Carman) che dice: "Se il sale riempie il 10% del buco, il tubo diventa significativamente più stretto". Ciò consente al modello di prevedere come il flusso rallenterà o si fermerà completamente man mano che l'ingorgo di zucchero peggiora.

Cosa hanno mostrato le Simulazioni

Il team ha testato questo strumento in due scenari principali:

Scenario A: L'Autostrada Interrotta (Crepe Disconnesse)

• Configurazione: Immaginate un serbatoio in cui le crepe non sono connesse tra loro; l'acqua deve incunearsi attraverso la roccia solida tra di esse.
• Risultato: Quando hanno pompato acqua fredda, l'acqua calda vicino al pozzo di produzione è evaporata rapidamente. Ciò ha causato la cristallizzazione del sale che ha intasato la roccia proprio intorno al pozzo.
• Il Colpo di Scena: Se hanno pompato l'acqua più velocemente, l'intasamento è peggiorato notevolmente e la produzione di energia è calata drasticamente. Il modello ha mostato che l'ingorgo è avvenuto principalmente nella roccia vicino al pozzo, non solo nelle crepe.

Scenario B: L'Autostrada Connessa (Crepe Connesse)

• Configazione: Immaginate un serbatoio in cui le crepe formano un'autostrada continua ad alta velocità dal pozzo di iniezione al pozzo di produzione.
• Risultato: L'acqua fredda è corsa attraverso le crepe rapidamente. Poiché si muoveva velocemente e rimaneva fresca, ha effettivamente sciolto il sale vicino al pozzo di produzione invece di intasarlo!
• Il Colpo di Scena: La precipitazione del sale si è spostata in un punto diverso, proprio al limite della zona d'acqua fredda, invece di intasare il pozzo stesso. Ciò suggerisce che avere una rete di crepe connessa potrebbe effettivamente proteggere il pozzo dall'intasamento, anche se cambia il luogo in cui il sale si accumula.

In Sintesi

Questo articolo presenta un nuovo strumento software open-source che aiuta gli ingegneri a comprendere la complessa danza tra calore, pressione e sale nei pozzi geotermici. Mappando accuratamente come le crepe si collegano e come il sale le intasa, lo strumento può aiutare a prevedere:

1. Dove i pozzi potrebbero essere bloccati dal sale.
2. Quanta energia può essere estratta in sicurezza prima che i "tubi" si intasino.
3. Se la disposizione delle crepe sotterranee aiuterà o ostacolerà il processo di produzione.

Gli autori hanno verificato il loro strumento rispetto a uno standard consolidato del settore e hanno scoperto che corrispondeva perfettamente, dimostrando che si tratta di un modo affidabile per simulare questi ambienti di roccia fratturata ad alta temperatura e salinità.

Sommario tecnico

Sintesi Tecnica: Modello di Flusso Composizionale Unificato per la Simulazione di Sistemi Multifase ad Alta Entalpia con Fratture

Problematica
La simulazione di serbatoi geotermici fratturati ad alta entalpia presenta sfide significative a causa del complesso accoppiamento tra il flusso multifase, multicomponente e non isotermico, la termodinamica fortemente non lineare e l'influenza dominante delle reti di fratture. Un problema operativo critico in questi sistemi è la formazione di incrostazioni minerali, specificamente la precipitazione di alite (NaCl). Durante la produzione, il calo di pressione induce l'ebollizione, concentrando la salamoia residua e superando i limiti di solubilità, mentre l'iniezione di acqua fredda abbassa la solubilità locale. Entrambi i meccanismi portano alla precipitazione di sale all'interno dei pori e delle fratture, compromettendo gravemente porosità e permeabilità, riducendo la produttività dei pozzi e potenzialmente causando il blocco del pozzo stesso. Gli esistenti modelli a continuo spesso non riescono a catturare la precipitazione localizzata e l'evoluzione della permeabilità su scala di frattura, mentre i tradizionali simulatori composizionali che si affidano a formulazioni di commutazione delle variabili (variable-switching) incontrano difficoltà di stabilità numerica in prossimità dei confini di fase.

Metodologia
Gli autori presentano un nuovo modello di flusso composizionale unificato implementato all'interno del framework open-source PorePy. La metodologia integra tre componenti fondamentali:

1. Formulazione a Variabili Persistenti: A differenza degli approcci tradizionali che scambiano le variabili primarie in base alla presenza di fase, questo modello utilizza un set fisso di variabili primarie — pressione ($p$), entalpia specifica ($h$) e frazione di massa salina complessiva ($z$) — attraverso tutte le configurazioni di fase (liquido, vapore e alite solida). Questo approccio gestisce naturalmente le transizioni di fase senza necessità di commutazione manuale, migliorando la robustezza e la stabilità numerica.
2. Approccio Discrete Fracture-Matrix (DFM): Il dominio fisico è rappresentato mediante una geometria a dimensioni miste in cui le fratture e le loro intersezioni sono trattate come varietà a dimensione inferiore (1D e 0D) immerse nella matrice porosa 2D. Questa rappresentazione esplicita preserva la connettività e l'eterogeneità delle fratture, permettendo di risolvere processi localizzati come la riduzione della permeabilità indotta dalla precipitazione.
3. Chiusura Termodinamica e Linearizzazione: Il modello impiega le correlazioni complete di Driesner e Christoph (2007) per la termodinamica H2O–NaCl, coprendo condizioni magmatiche. Per ridurre il costo computazionale, il modello utilizza una strategia di Linearizzazione Basata su Operatore (OBL). Le proprietà termodinamiche e le loro derivate vengono precomputate su una griglia strutturata nello spazio delle variabili primarie e recuperate tramite interpolazione multilineare durante la simulazione, eliminando i costosi calcoli di separazione di fase (flash) effettuati on-the-fly.
4. Evoluzione Dinamica di Permeabilità e Apertura: Il modello incorpora relazioni di tipo Kozeny-Carman per aggiornare dinamicamente la porosità e la permeabilità della matrice in base alla saturazione di alite. Per le fratture, l'apertura evolve secondo una relazione di legge di potenza dipendente dalla saturazione di alite, che a sua volta aggiorna la permeabilità della frattura tramite la legge cubica.

Le equazioni governanti (bilancio di massa ed energia) sono discretizzate utilizzando uno schema di volume finito a celle centrate con l'approssimazione di flusso multi-punto (MPFA) per i termini diffusivi e l'upwinding per i termini avvettivi. Il sistema non lineare risultante viene risolto utilizzando il metodo di Newton con una ricerca di linea Armijo con backtracking e passi temporali adattivi.

Risultati Chiave
Il modello è stato verificato e applicato attraverso diversi esperimenti numerici:

• Verifica (Benchmark 1D): Il modello è stato verificato rispetto al simulatore a codice chiuso CSMP++ utilizzando un benchmark di dissoluzione del sale 1D. I risultati hanno mostrato un forte accordo attraverso quattro distinte regioni di fase (vapore+alite, vapore+liquido, fase liquida singola e liquido+alite) e le loro transizioni, validando la formulazione a variabili persistenti e l'implementazione termodinamica.
• Rete di Fratture Disconnesse (2D): In un serbatoio con fratture isolate, il flusso è dominato dalla matrice. La simulazione ha dimostrato che il calo di pressione al pozzo di produzione induce l'ebollizione, portando a una significativa precipitazione di alite nella matrice circostante e nei segmenti di frattura adiacenti. Questa precipitazione ha ridotto la permeabilità della matrice fino al 36% ($K/K_0 \approx 0,64$) e ha alterato le linee di flusso. L'aumento della velocità di iniezione ha intensificato l'ebollizione vicino al pozzo, portando a una maggiore saturazione di alite ($s_{hal} \approx 0,46$) e a una riduzione più severa della permeabilità ($K/K_0 \approx 0,29$), impattando significativamente i tassi di produzione.
• Rete di Fratture Connesse (2D): In un serbatoio con una catena di fratture continua che collega i pozzi di iniezione e produzione, il trasporto è diventato dominato dalle fratture. La rete connessa ha facilitato la rapida comunicazione della pressione e il trasporto di fluido a bassa salinità verso il pozzo di produzione. Di conseguenza, mentre l'ebollizione si è verificata vicino al produttore, il fluido arrivato è rimasto sufficientemente diluito da prevenire la saturazione locale di alite. Inveve della precipitazione al produttore, è stata osservata la dissoluzione, e la precipitazione è stata ridistribuita nella matrice davanti al fronte di iniezione e lungo la catena di fratture a causa degli effetti di raffreddamento. Questa configurazione ha evitato la severa compromissione della permeabilità osservata nel caso disconnesso.

Significatività e Rivendicazioni
L'articolo sostiene di fornire uno strumento numerico robusto e open-source per l'analisi del complesso trasporto di calore e massa con effetti di incrostazione minerale in sistemi geotermici fratturati ad alta entalpia. La sua principale significatività risiede in:

• Stabilità Numerica: La formulazione unificata a variabili persistenti evita le difficoltà di convergenza associate alla commutazione delle variabili in prossimità dei confini di fase.
• Fedeltà Geometrica: L'approccio DFM risolve esplicitamente la connettività delle fratture, il che gli autori dimostrano essere un fattore critico che controlla la distribuzione spaziale della precipitazione e la conseguente riduzione della trasmissibilità.
• Intuizione Operativa: I risultati evidenziano come la connettività delle fratture determini le sfide operative; le reti disconnesse sono soggette a intasamenti vicino al pozzo e perdita di produttività, mentre le reti connesse possono mitigare tali problemi mantenendo la diluizione del fluido al pozzo di produzione.
• Riproducibilità: Il lavoro fornisce il modello numerico completo e il workflow di esecuzione tramite un contenitore Docker, facilitando la riproducibilità e l'ulteriore sviluppo all'interno della comunità geotermica.

Gli autori concludono che il loro modello predice efficacemente i pattern di precipitazione dell'alite e il loro impatto sulle prestazioni del serbatoio, offrendo uno strumento prezioso per comprendere l'interazione tra flusso, termodinamica e incrostazione minerale nei serbatoi geotermici fratturati.