A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

Nel mondo nascosto sotto i nostri piedi, l'acqua e le sostanze chimiche disciolte si muovono attraverso un complesso labirinto di roccia. Questa roccia è raramente uniforme; è spesso un blocco solido, chiamato matrice, attraversato da sottili crepe note come fratture. Queste fratture sono le autostrade del sottosuolo. Possono agire come canali superveloci che fanno scorrere rapidamente i fluidi attraverso la roccia, o come barriere che bloccano completamente il movimento, a seconda della loro natura. Poiché queste crepe sono spesso microscopiche in larghezza ma possono estendersi per chilometri in lunghezza, modellare il modo in cui le sostanze viaggiano attraverso di esse è una sfida enorme per gli scienziati. Se provano a costruire un modello informatico che catturi ogni minimo dettaglio dell'apertura della crepa, il calcolo diventa così pesante che nemmeno i supercomputer più potenti riescono a risolverlo. Per aggirare questo problema, i ricercatori utilizzano una scorciatoia ingegnosa chiamata modello ibrido-dimensionale. Invece di trattare la crepa come una sottile fetta tridimensionale, la appiattiscono in una linea bidimensionale o in una linea monodimensionale che si trova all'interno della roccia. Ciò consente loro di simulare il flusso senza dover risolvere l'infinitesima larghezza della crepa stessa.

La difficoltà sorge quando gli scienziati cercano di collegare la matematica della roccia solida con la matematica di queste fratture appiattite. Nel mondo reale, la roccia e la crepa sono in costante conversazione, scambiando fluidi e sostanze chimiche attraverso il loro confine condiviso. Per simulare questo, i computer di solito dividono la roccia e la crepa in piccoli pezzi, come un mosaico, per calcolare il flusso. Per anni, l'approccio standard ha richiesto che i pezzi della roccia e i pezzi della crepa si allineassero perfettamente lungo il loro bordo comune, come due pezzi di un puzzle che si incastrano. Se la roccia era divisa in grandi quadrati e la crepa in segmenti minuscoli, il computer doveva costringerli a coincidere, spesso creando uno strato artificiale e complesso di dati per colmare il divario. Questo requisito rendeva i modelli rigidi e difficili da usare quando diverse parti di un sistema dovevano essere raffinate a scale differenti.

Un nuovo studio di Alvaro H. Salas e Carlos D. Acosta dell'Universidad Nacional de Colombia affronta questa rigidità ponendo una domanda semplice ma profonda: e se la roccia e la crepa non dovessero allinearsi affatto? Hanno sviluppato un metodo in cui la roccia e la frattura possono essere divise in pezzi di dimensioni e forme completamente diverse, generati indipendentemente l'uno dall'altro, eppure scambiare informazioni perfettamente. Immaginate di cercare di versare l'acqua da un secchio grande e grossolano in un setaccio fine e dettagliato. Di solito, avreste bisogno di un imbuto speciale per far sì che l'acqua si trasferisca senza versarsi. I ricercatori hanno trovato un modo per versare l'acqua direttamente, calcolando l'esatta quantità che passa attraverso ogni punto in cui il bordo del secchio tocca il setaccio, indipendentemente da quanto siano disallineati i fori.

Il team ha creato un quadro matematico che tratta lo scambio di sostanze chimiche tra la roccia e la frattura come un'interazione fisica diretta. Invece di costringere le due griglie a coincidere o di usare un intermediario per tradurre i dati, hanno costruito un sistema che calcola il trasferimento sulle minuscole sezioni sovrapposte dove il bordo della roccia e la superficie della frattura si incontrano effettivamente. Questo approccio è chiamato metodo a griglia non corrispondente (nonmatching-grid). Nelle loro simulazioni, hanno testato questo sistema con reazioni chimiche complesse, in cui le sostanze cambiano mentre si muovono, uno scenario notoriamente difficile da modellare. Hanno dimostrato matematicamente che questa connessione diretta è stabile e accurata, anche quando la roccia è divisa in blocchi molto grandi e la frattura in segmenti molto fini, o viceversa. La stabilità del loro metodo non dipende dal rapporto tra le dimensioni di questi blocchi, il che significa che il computer non va in crash o produce errori selvaggi solo perché una griglia è molto più fine dell'altra.

Uno dei risultati più significativi del loro lavoro è che questo metodo diretto è superiore ai vecchi approcci indiretti che tentano di proiettare i dati da una griglia all'altra. Quando lo scambio di sostanze chimiche è forte, i vecchi metodi di proiezione perdono informazioni. Essenzialmente attenuano i dettagli della frattura, scartando le variazioni su piccola scala che sono cruciali per una previsione accurata. I ricercatori hanno dimostrato che il loro metodo diretto conserva tutte le informazioni necessarie. In test in cui il trasferimento di sostanze chimiche era intenso, il vecchio metodo ha prodotto errori del ventuno per cento superiori rispetto al loro nuovo metodo. Inoltre, quando hanno raffinato la griglia della frattura per renderla ancora più dettagliata mantenendo fissa la griglia della roccia, il vecchio metodo ha smesso di migliorare, restando in uno stato di stagnazione. Il loro nuovo metodo, invece, continuava a diventare più accurato, dimostrando di catturare la vera fisica dell'interazione senza limitazioni artificiali.

I ricercatori hanno anche dimostrato che questo metodo funziona efficientemente sui computer. Hanno scritto un programma software per testare la loro teoria, eseguendo simulazioni su un normale computer desktop. Hanno confrontato il tempo necessario per impostare i calcoli con griglie corrispondenti rispetto a griglie non corrispondenti. I risultati hanno mostato che il lavoro extra richiesto per gestire i pezzi disallineati era trascurabile. Il tempo necessario per assemblare le equazioni era quasi identico, il che significa che gli scienziati possono ora utilizzare griglie molto più flessibili e realistiche senza pagare un prezzo in termini di tempo di calcolo. Questa flessibilità è vitale per applicazioni del mondo reale, dove diverse parti di una formazione geologica potrebbero essere modellate da team diversi utilizzando software differenti, o dove la rete di fratture è così complessa che una singola griglia perfettamente allineata è impossibile da creare.

Lo studio ha anche esplorato come il metodo si comporti in condizioni estreme. Hanno testato scenari in cui le reazioni chimiche erano molto forti e il trasferimento di fluidi era rapido. In questi casi difficili, il metodo è rimasto robusto, convergendo verso una soluzione senza fallire. Hanno persino testato la geometria delle fratture, simulando crepe che si piegano o si intersecano tra loro. Sebbene i test principali abbiano utilizzato una crepa dritta, la logica sottostante del loro metodo ha retto quando applicata a queste forme più complesse, suggerendo che l'approccio potrebbe essere esteso alle reti aggrovigliate presenti in natura. I ricercatori sono stati attenti a sottolineare che il loro lavoro attuale si concentra su tipi specifici di rocce e geometrie di frattura, e che estendere il metodo a reti tridimensionali, curve e intersecanti è un compito per il futuro.

Ciò che rende questo lavoro particolarmente importante è che rimuove una barriera di lunga data nella simulazione di mezzi fratturati. Per decenni, la necessità di griglie corrispondenti ha limitato il modo in cui gli scienziati potevano modellare questi sistemi, spesso costringendoli a scendere a compromessi sulla precisione della frattura o della roccia per far funzionare la matematica. Dimostrando che i due possono essere accoppiati direttamente e in modo conservativo senza una corrispondenza perfetta, i ricercatori hanno aperto la porta a simulazioni più realistiche e flessibili. Ciò significa che in futuro, i modelli di contaminazione delle acque sotterranee, di stoccaggio del carbonio o di estrazione di energia geotermica possono essere costruiti con maggiore libertà, permettendo alla roccia e alle crepe di essere rappresentate esattamente come sono, piuttosto che come devono essere per adattarsi a un rigido vincolo matematico. Il metodo fornisce una nuova, stabile base per comprendere come i fluidi e le sostanze chimiche si muovono attraverso la terra fratturata, assicurando che la connessione tra la roccia solida e le crepe sia tanto fluida nel computer quanto lo è nel terreno.