A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

Profondamente nel sottosuolo, nella complessa rete di rocce che contengono olio, gas e acqua, i fluidi non si muovono sempre in linee semplici e prevedibili. La roccia stessa è raramente uniforme; è un mosaico di materiali con capacità di far passare i fluidi molto diverse tra loro. Alcuni strati sono come autostrade aperte, mentre altri sono come muri densi. Quando gli ingegneri cercano di modellare il modo in cui questi fluidi scorrono, devono tenere conto di questi cambiamenti improvvisi nelle proprietà della roccia. Inoltre, i fluidi stessi possono essere complicati. Mentre l'acqua scorre facilmente, molti fluidi industriali — come le sospensioni utilizzate nella produzione cartaria o i fanghi di perforazione usati nei pozzi petroliferi — si comportano diversamente. Sono "non newtoniani", il che significa che la loro viscosità e la loro resistenza al flusso cambiano a seconda di quanto vengono spinti. Quando questi fluidi complessi incontrano una roccia che cambia carattere bruscamente, la matematica che ne descrive il movimento diventa incredibilmente difficile, spesso interrompendosi proprio nei punti in cui le proprietà della roccia cambiano.

Per decenni, gli scienziati si sono affidati a un modo standard di calcolare questi flussi, un metodo che essenzialmente media le proprietà della roccia tra due punti. Questo approccio funziona abbastanza bene quando la roccia è relativamente costante, ma fallisce clamorosamente quando la roccia è altamente variabile. In quei casi, il metodo standard sottostima costantemente quanto fluido possa effettivamente muoversi, portando a errori che possono essere enormi. Per risolvere questo problema, i ricercatori della University of Science and Technology of China hanno sviluppato una nuova tecnica computazionale. Invece di affidarsi a semplici medie, questo nuovo metodo utilizza una descrizione matematica precisa di come il fluido si comporta proprio al bordo dove incontrano due tipi diversi di roccia. Catturando il vero, netto comportamento del fluido in questi confini critici, il nuovo metodo fornisce un quadro molto più chiaro e accurato del flusso sotterraneo, anche quando la roccia è estremamente irregolare e il fluido è complesso.

I ricercatori si sono concentrati su un tipo specifico di fluido noto come fluido a legge di potenza (power-law fluid). Questi sono fluidi in cui la relazione tra la forza che li spinge e la velocità con cui si muovono segue una regola specifica e non lineare. Nel mondo reale, questo copre una vasta gamma di sostanze, dalle soluzioni polimeriche dense a certi tipi di petrolio greggio. La sfida sorge quando questi fluidi scorrono attraverso mezzi porosi — rocce piene di minuscoli fori — dove la permeabilità, ovvero la facilità di flusso, cambia improvvisamente da un punto all'altro. In tali ambienti, la pressione del fluido non cambia in modo fluido; invece, crea un picco acuto o una "singolarità" proprio al confine tra la roccia a flusso facile e quella a flusso difficile. I modelli computerizzati tradizionali, che dividono la roccia in una griglia di piccoli triangoli, faticano con questi picchi. Essi assumono che la pressione cambi in linea retta tra i punti, il che è fisicamente errato in questi confini. Di conseguenza, i modelli producono portate troppo basse e richiedono un numero enorme di minuscoli punti di griglia per ottenere anche solo un'approssimazione approssimativa della verità.

Per risolvere il problema, il team ha sviluppato un Metodo Analitico Finito. L'idea centrale è smettere di indovinare come si comporta il fluido al confine e utilizzare invece una descrizione matematica esatta di tale comportamento. I ricercatori si sono resi conto che, vicino al punto in cui le proprietà della roccia cambiano, la pressione del fluido segue un particolare schema curvo, salendo o scendendo in un modo che dipende dalla natura a legge di potenza del fluido. Hanno calcolato la forma esatta di questa curva, nota come autofunzione, per la parte più critica del flusso. Questo calcolo ha comportato la risoluzione di un insieme complesso di equazioni che descrivono come il fluido si muove in un'area a forma di cuneo dove le proprietà della roccia differiscono. Trovando la soluzione più piccola e significativa di queste equazioni, hanno catturato il modo dominante in cui il fluido si comporta proprio nel punto critico.

Una volta compreso questo comportamento locale, i ricercatori hanno costruito un nuovo schema numerico attorno ad esso. Invece di utilizzare il vecchio e difettoso metodo di media per collegare i punti della griglia, hanno utilizzato la loro nuova soluzione analitica per calcolare il flusso attraverso i confini di ogni cella triangolare della loro griglia. Ciò ha permesso al computer di "sapere" esattamente come il fluido dovrebbe comportarsi mentre si infiltra attraverso il collo stretto tra zone ad alta e bassa permeabilità. Il risultato è una simulazione che rispetta la fisica del problema nei punti più difficili, invece di attenuarli con un'approssimazione.

Il team ha testato questo nuovo metodo contro diversi scenari impegnativi per vedere come si fosse comportato. Hanno iniziato con un classico schema a "scacchiera", dove rocce ad alta permeabilità e a bassa permeabilità erano disposte in quadrati alternati. In questa configurazione, le proprietà della roccia cambiavano di un fattore di cento. Quando hanno confrontato il loro nuovo metodo con quello tradizionale, la differenza è stata netta. Il metodo tradizionale sottostimava il flusso per margini enormi, mancando talvolta il valore reale di quasi il novanta per cento su una griglia grossolana. Al contrario, il nuovo metodo rimaneva sorprendentemente vicino alla risposta corretta, anche quando la griglia non era molto fine. Man mano che rendevano la griglia più fine, il nuovo metodo convergevano verso la risposta corretta molto rapidamente, mentre il metodo tradizionale faticava a stare al passo, anche con griglie estremamente fini.

Hanno anche testato il metodo su un dominio esagonale regolare e su un dominio con una distribuzione log-normale casuale della permeabilità, che imita la natura caotica dei veri serbatoi sotterranei. In ogni caso, il nuovo metodo si è dimostrato superiore. Ha gestito con facilità le variazioni estreme delle proprietà della roccia e ha mantenuto l'accuratezza indipendentmente dal tipo specifico di fluido a legge di potenza simulato. Uno dei risultati più sorprendenti è stato che l'accuratezza del nuovo metodo non dipendeva pesantemente dalle proprietà specifiche del fluido o della roccia. Che il fluido fosse molto denso o molto fluido, o che il contrasto della roccia fosse lieve o estremo, il metodo forniva risultati coerenti. I ricercatori hanno scoperto che un livello moderato di raffinamento della griglia era sufficiente per ottenere un'elevata accuratezza, con errori che scendevano a meno dell'otto per cento anche negli scenari peggiori.

Questo lavoro rappresenta un passo avanti significativo nel modo in cui simuliamo il flusso di fluidi in ambienti sotterranei complessi. Sostituendo una tecnica di media difettosa con una soluzione che rispetta la vera fisica del confine, i ricercatori hanno creato uno strumento che è sia più accurato che più efficiente. Per le industrie che dipendono dalla comprensione di come i fluidi si muovono attraverso la terra, dall'estrazione energetica alla gestione ambientale, ciò significa che le simulazioni possono essere affidabili e più veloci. Il metodo non richiede la massiccia potenza di calcolo che era precedentemente necessaria per ottenere una risposta decente, rendendo la modellazione ad alta fedeltà di sistemi complessi ed eterogenei molto più accessibile. I risultati suggeriscono che, guardando da vicino i punti singolari dove il mondo cambia, possiamo costruire una migliore comprensione dell'insieme.