Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Profondamente sotto la superficie terrestre, nelle rocce porose che contengono petrolio, gas e calore geotermico, una battaglia silenziosa viene combattuta costantemente tra la pietra solida e i fluidi intrappolati nei suoi minuscoli pori. Quando gli ingegneri cercano di estrarre queste risorse, spesso pompano fluidi nella roccia per separarla, un processo noto come fratturazione idraulica. Per decenni, gli scienziati hanno modellato questo processo trattando la roccia come una semplice spugna elastica che si allunga e torna in forma. Tuttavia, le rocce reali sono più complesse; possono anche piegarsi permanentemente, sbriciolarsi e scivolare lungo crepe interne, specialmente quando vengono schiacciate da alte pressioni. Comprendere esattamente come queste rocce si rompono quando sia lo scheletro solido che il fluido all'interno spingono e tirano è fondamentale per trivellare pozzi in sicurezza, gestire lo stoccaggio dei rifiuti e persino prevedere eventi geologici naturali come la rottura degli iceberg o il passaggio del magma attraverso la crosta. La sfida consiste nel creare una descrizione matematica che catturi la realtà disordinata della rottura delle rocce senza perdersi in calcoli impossibili.

Un team di ricercatori ha rivisitato le equazioni fondamentali utilizzate per descrivere questo processo di rottura, scoprendo che una scorciatoia comune utilizzata nei modelli precedenti porta a un errore significativo. Nel mondo della meccanica delle rocce, gli scienziati utilizzano un approccio a "campo di fase" (phase-field), che tratta una crepa non come una linea netta e invisibile, ma come una zona sfocata in cui il materiale si sta gradualmente indebolendo. Ciò consente ai computer di simulare come una crepa cresca e cambi forma. Per farlo accuratamente, il modello deve tenere conto di due comportamenti distinti: quando la roccia viene tirata (tensione) e quando viene schiacciata (compressione). In tensione, le piccole crepe si aprono; in compressione, si chiudono e le superfici ruvide all'interno di esse scivolano l'una contro l'altra, creando attrito. I ricercatori hanno scoperto che molti modelli esistenti, ampiamente utilizzati nell'industria e nell'accademia, confondono le variabili utilizzate per descrivere l'energia del sistema. Nello specifico, hanno scoperto che trattare la pressione del fluido come un input diretto nell'equazione dell'energia, piuttosto che la quantità di fluido nei pori, crea un "salto" o una discontinuità matematica. Questo salto significa che il modello prevede che la resistenza della roccia cambi bruscamente e in modo irrealistico quando passa dal subire una trazione al subire una compressione, portando a previsioni imprecise di come e dove si formerà una frattura.

Per risolvere questo problema, i ricercatori hanno sviluppato una nuova formulazione coerente che traccia attentamente la relazione tra la deformazione della roccia e il contenuto di fluido. Hanno costruito il loro modello su una visione microscopica della roccia, immaginandola come una collezione di minuscole crepe a forma di moneta riempite di fluido. Analizzando come queste micro-crepe si aprono e si chiudono sotto diverse sollecitazioni, hanno derivato un insieme di regole che garantiscono che la superficie di resistenza della roccia — il confine matematico che definisce quando si romperà — rimanga fluida e continua, indipendentemente dalle variazioni di pressione. Hanno dimostrato che utilizzando un tipo specifico di regola di flusso che collega la deformazione permanente della roccia direttamente alle forze che agiscono su di essa, e includendo correttamente un termine di accoppiamento mancante nell'equazione dell'energia, il modello può descrivere accuratamente la transizione dalla tensione alla compressione. Questa correzione non è un semplice ritocco minore; cambia fondamentalmente il modo in cui il modello calcola la forza motrice che spinge una crepa in avanti.

Il team ha testato il nuovo modello contro due scenari molto diversi per vedere se reggeva. Per prima cosa, hanno simulato un classico benchmark di fratturazione idraulica in cui un fluido viene iniettato in una roccia per creare una crepa. Hanno confrontato i loro risultati con una soluzione analitica nota, che è una risposta matematica precisa derivata per condizioni ideali. Il nuovo modello corrispondeva a questa soluzione molto più da vicino rispetto alle vecchie formulazioni miste, prevedendo correttamente la pressione necessaria per iniziare la frattura e la larghezza della crepa durante la sua crescita. I modelli precedenti tendevano a sovrastimare la pressione del fluido necessaria, un difetto che potrebbe portare a operazioni di trivellazione inefficienti o insicure. Nel secondo test, hanno simulato un campione di roccia schiacciato da due lati, una condizione che imita l'ambiente sotterraneo profondo dove le rocce sono sottoposte a una pressione immensa. In questo caso, la roccia non si è limitata a spezzarsi; si è deformata plasticamente, formando una banda di taglio in cui il materiale è scivolato su se stesso prima di rompersi. Il nuovo modello ha riprodotto con successo questo comportamento complesso, catturando sia le fratture da taglio causate dallo schiacciamento meccanico che le fratture da tensione guidate dall'iniezione di fluido.

I risultati di queste simulazioni rivelano che la deformazione plastica della roccia agisce come una barriera significativa alla crescita della frattura. Quando alla roccia è permesso deformarsi permanentemente, essa assorbe più energia, il che significa che la pressione del fluido all'interno della crepa deve essere più alta per mantenerla aperta e in crescita. Questa resistenza extra cambia l'intero profilo della frattura, rendendola più stretta e richiedendo più forza per propagarsi rispetto a un modello puramente elastico. Inoltre, nei test di compressione, il modello ha mostrato come la pressione dei pori all'interno della roccia cambi durante la formazione della frattura. Mentre la roccia si deforma e crea spazio per il movimento del fluido, la pressione scende localmente, attirando più fluido nella zona danneggiata. Questa interazione crea un ciclo di feedback in cui il flusso del fluido e il cedimento meccanico della roccia sono strettamente accoppiati, una dinamica che il nuovo modello cattura con alta fedeltà.

Correggendo la base matematica di questi modelli, i ricercatori hanno fornito uno strumento più affidabile per gli ingegneri e gli scienziati che lavorano con i mezzi porosi. La nuova formulazione assicura che la transizione tra i diversi tipi di cedimento della roccia sia fluida e fisicamente realistica, eliminando i salti artificiali che affliggevano i tentativi precedenti. Questa precisione è vitale per applicazioni che vanno dall'ottimizzazione della produzione di petrolio e gas alla progettazione di uno stoccaggio sotterraneo sicuro per la cattura del carbonio. Lo studio dimostra che anche in sistemi complessi dove i solidi e i fluidi interagiscono, un approccio attento e coerente alla fisica sottostante può rivelare la vera natura di come i materiali si rompono. Il modello è ora pronto ad aiutare a prevedere il comportamento della crosta terrestre con maggiore fiducia, colmando il divario tra la meccanica teorica e la realtà plastica e disordinata del terreno sotto i nostri piedi.