Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

Il terreno sotto i nostri piedi è raramente un blocco solido e ininterrotto. In molti luoghi, dai profondi giacimenti petroliferi alla crosta stessa della Terra, la roccia è attraversata da fratture: crepe sottili e planari che possono estendersi per chilometri. Quando l'energia viaggia attraverso questo paesaggio fratturato, che si tratti di un'onda sismica causata da un terremoto o di un impulso sonoro usato per ispezionare un materiale, le crepe cambiano la storia. Esse disperdono l'energia, indeboliscono il segnale e talvolta causano anche lo sfregamento o lo scivolamento delle facce della roccia l'una contro l'altra. Comprendere esattamente come queste onde interagiscano con le crepe è vitale per i geologi che cercano di mappare ciò che si trova nel sottosuolo e per gli ingegneri che garantiscono la sicurezza delle strutture. Tuttavia, la fisica di questa interazione è notoriamente difficile da definire con precisione. Una semplice crepa non è solo un vuoto; è un confine complesso dove le facce della roccia possono premere l'una contro l'altra, scorrere con attrito o aprirsi e chiudersi bruscamente, il tutto mentre l'onda passa in una frazione di secondo.

Per decenni, gli scienziati si sono affidati a modelli semplificati per prevedere questi comportamenti, trattando spesso le crepe come se fossero fatte di molle invisibili che si tendono e si comprimono. Sebbene utili, questi modelli falliscono quando le onde sono abbastanza forti da costringere le facce della roccia a toccarsi, scivolare o separarsi completamente. Non riescono a catturare il movimento di "stick-slip" (attacca-scivola) in cui l'attrito tiene ferma la roccia finché la forza non diventa troppo grande, né possono descrivere accuratamente l'irrigidimento non lineare che si verifica quando una crepa viene schiacciata. Un nuovo studio condotto da ricercatori dell'Università di Bergen, del Centro di Ricerca Norvegese e dell'Università Tecnica di Monaco affronta queste lacune costruendo un quadro computazionale completo che simula onde elastiche che si muovono attraverso rocce fratturate con un realismo senza precedenti. Invece di affidarsi a ipotesi semplificate, il team ha sviluppato un metodo in grado di gestire quattro distinti livelli di complessità, che vanno dal comportamento basico simile a una molla fino alla piena meccanica del contatto, dove l'attrito e l'apertura e chiusura fisica delle crepe vengono calcolati in tempo reale.

I ricercatori hanno costruito un ambiente digitale unificato in cui la roccia e le crepe sono trattate come un sistema unico e interconnesso. Nelle loro simulazioni, la roccia è rappresentata come un volume tridimensionale, mentre le fratture sono modellate come superfici bidimensionali incorporate in esso. Ciò consente al computer di tracciare come un'onda che viaggia attraverso la roccia colpisca una crepa, come la crepa si deformi e come tale deformazione invii increspature nuovamente nella roccia. Il team ha testato quattro modi specifici in cui le crepe potrebbero comportarsi. I primi due modelli utilizzavano leggi basate su molle, dove la forza attraverso la crepa è proporzionale a quanto la roccia si muove, con una versione che permetteva una semplice relazione lineare e un'altra che utilizzava una regola non lineare più complessa che tiene conto del fatto che le crepe diventano più difficili da schiacciare man mano che si chiudono. Gli altri due modelli introducevano la realtà dell'attrito e del contatto. In questi scenari, le facce della crepa possono toccarsi e scivolare l'una contro l'altra, governate da leggi di attrito simili a quelle che impediscono agli pneumatici di un'auto di slittare, applicando anche una regola per cui le facce della roccia non possono attraversarsi l'un l'altra.

Per garantire l'accuratezza del loro nuovo framework, il team ha eseguito una serie di test rigorosi. Hanno simulato onde che colpivano una singola crepa e hanno confrontato i risultati del computer con soluzioni matematiche note per casi semplici. Hanno scoperto che il loro metodo poteva prevedere il comportamento dell'onda con alta precisione, corrispondendo alle aspettative teoriche sia per la velocità dell'onda che per la quantità di energia che passava attraverso o rimbalzava indietro. Hanno poi spinto il sistema oltre, simulando onde abbastanza forti da far aprire, aderire e scivolare le crepe. In questi scenari più complessi, dove le facce della crepa strinano tra loro con un coefficiente d'attrito di 1,0, le simulazioni hanno mostrato che il metodo rimaneva stabile e produceva risultati coerenti. I ricercatori hanno verificato che il loro approccio poteva gestire i cambiamenti improvvisi della fisica che si verificano quando una crepa passa dallo stato di adesione a quello di scivolamento, una transizione che spesso mette in crisi i modelli più semplici.

Quando il team ha confrontato i quattro diversi modelli fianco a fianco, le differenze nei modelli d'onda risultanti sono state sorprendenti. Le simulazioni hanno rivelato che la scelta del modello conta significamente. Quando una crepa veniva lasciata chiudere e irrigidire secondo la legge non lineare, trasmetteva molta più energia dell'onda rispetto a una crepa modellata con una semplice molla lineare. Questo perché il modello non lineare tiene conto del fatto che, man mano che la crepa si chiude, le facce della roccia premono più strettamente, rendendo la crepa effettivamente più rigida e trasparente all'onda. Inoltre, l'inclusione dell'attrito ha cambiato la forma dei fronti d'onda. Nei modelli con attrito, l'onda che emergeva dall'altro lato della crepa era più ampia e distorta rispetto alle versioni senza attrito. I ricercatori hanno osservato che l'attrito causava l'adesione delle facce della crepa per periodi più lunghi, ritardando il rilascio di energia e alterando il modo in cui l'onda si disperdeva.

Lo studio non si è fermato alle singole crepe. I ricercatori hanno applicato il loro framework a un dominio tridimensionale contenente sei fratture intersecanti, creando una piccola rete di crepe con proprietà differenti. In questo contesto complesso, hanno assegnato diversi limiti di chiusura massima a diversi set di fratture, il che significa che alcune crepe potevano chiudersi più strettamente di altre. La simulazione ha mostrato che il campo d'onda è diventato altamente eterogeneo. Le fratture che erano più rigide e potevano chiudersi più strettamente riflettevano più energia e trasmettevano meno, creando ombre distinte nel pattern dell'onda. I ricercatori sono stati in grado di tracciare esattamente quanto ogni crepa si fosse chiusa e quanto avesse scivolato, rivelando che anche piccole differenze nelle proprietà fisiche delle crepe portavano a risposte di deformazione sensibilmente diverse. Ad esempio, una crepa che poteva chiudersi più strettamente mostrava un comportamento di "adesione" maggiore, mantenendo le proprie facce unite più a lungo prima di scivolare, il che a sua volta influenzava l'onda che la attraversava.

Le implicazioni di queste scoperte si estendono oltre i risultati immediati. Dimostrando che un singolo framework computazionale può gestire tutto, dal semplice comportamento a molla alla complessa meccanica del contatto con attrito, i ricercatori hanno fornito uno strumento che può essere utilizzato per interpretare i dati sismici reali con maggiore fedeltà. La capacità di simulare come le onde si disperdono e si attenuano in presenza di crepe soggette a attrito, apertura e chiusura significa che i geoscientifici possono distinguere meglio tra diversi tipi di formazioni rocciose e fratture riempite di fluidi. Lo studio conferma che ignorare gli aspetti non lineari e di attrito della meccanica delle fratture può portare a errori significativi nella previsione di come l'energia si muove nel sottosuolo. Sebbene l'attuale lavoro si concentri sulle interazioni puramente meccaniche, il framework è costruito su una base compatibile con i modelli di flusso dei fluidi, suggerendo una strada verso simulazioni ancora più realistiche in cui l'acqua o il petrolio che si muove attraverso le crepe interagisce con l'onda passante. Il lavoro è una dimostrazione del fatto che, abbracciando la piena complessità del comportamento delle crepe nella roccia, possiamo costruire un'immagine più chiara del mondo nascosto sotto i nostri piedi.