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

Diep onder het aardoppervlak, in de poreuze gesteenten die olie, gas en geothermische warmte vasthouden, wordt een constante stille strijd gevoerd tussen vast gesteente en de vloeistoffen die gevangen zitten in de kleine poriën ervan. Wanneer ingenieurs deze hulpbronnen proberen te extraheren, pompen ze vaak vloeistoffen in het gesteente om het uit elkaar te drukken, een proces dat bekend staat als hydraulisch breken (hydraulic fracturing). Decennialang hebben wetenschappers dit proces gemodelleerd door het gesteente te behandelen als een eenvoudige, elastische spons die uitrekt en weer terugveert. Echter, echte gesteenten zijn complexer; ze kunnen ook permanent vervormen, verbrijzelen en langs interne scheuren glijden, vooral wanneer ze worden samengeperst door hoge druk. Het begrijpen van exact hoe deze gesteenten breken wanneer zowel het solide skelet als de vloeistof binnenin duwen en trekken, is cruciaal voor het veilig boren van putten, het beheren van afvalopslag en zelfs het voorspellen van natuurlijke geologische gebeurtenissen zoals het barsten van ijsbergen of het omhoog duwen van magma door de korst. De uitdaging ligt in het creëren van een wiskundige beschrijving die de rommelige realiteit van brekend gesteente vastlegt zonder verloren te raken in onmogelijke berekeningen.

Een team onderzoekers heeft de fundamentele vergelijkingen die worden gebruikt om dit breukproces te beschrijven, opnieuw onderzocht en ontdekt dat een veelgebruikte afkorting in eerdere modellen leidt tot een significante fout. In de wereld van de rotsmechanica gebruiken wetenschappers een "phase-field"-benadering, waarbij een scheur niet wordt behandeld als een scherpe, onzichtbare lijn, maar als een wazige zone waar het materiaal geleidelijk verzwakt. Dit stelt computers in staat om te simuleren hoe een scheur groeit en van vorm verandert. Om dit nauwkeurig te doen, moet het model rekening houden met twee verschillende gedragingen: wanneer het gesteente uit elkaar wordt getrokken (trekspanning) en wanneer het wordt samengedrukt (compressie). Bij trekspanning openen minuscule scheurtjes zich; bij compressie sluiten ze en glijden de ruwe oppervlakken binnenin langs elkaar, wat wrijving veroorzaakt. De onderzoekers ontdekten dat veel bestaande modellen, die breed worden gebruikt in de industrie en de academische wereld, de variabelen die de energie van het systeem beschrijven, door elkaar halen. Specifiek ontdekten zij dat het behandelen van de vloeistofdruk als een directe input in de energievergelijking, in plaats van de hoeveelheid vloeistof in de poriën, een wiskundige "sprong" of discontinuïteit creëert. Deze sprong betekent dat het model voorspelt dat de sterkte van het gesteente abrupt en onrealistisch verandert wanneer het overgaat van trek naar compressie, wat leidt tot onnauwkeurige voorspellingen over hoe en waar een breuk zal ontstaan.

Om dit op te lossen, hebben de onderzoekers een nieuwe, consistente formulering ontwikkeld die de relatie tussen de vervorming van het gesteente en de vloeistofinhoud zorgvuldig bijhoudt. Ze bouwden hun model op een microscopisch beeld van het gesteente, waarbij ze het voorstellen als een verzameling van kleine, muntvormige scheurtjes gevuld met vloeistof. Door te analyseren hoe deze micro-scheurtjes openen en sluiten onder verschillende spanningen, hebben ze een reeks regels afgeleid die ervoor zorgen dat het sterkteoppervlak van het gesteente — de wiskundige grens die bepaalt wanneer het zal breken — glad en continu blijft, ongeacht hoe de druk verandert. Ze bewezen dat door een specifiek type stroomregel te gebruiken die de permanente vervorming van het gesteente direct koppelt aan de werkende krachten, en door correct een ontbrekende koppelingsterm in de energievergelijking op te nemen, het model de overgang van trek naar compressie nauwkeurig kan beschrijven. Deze correctie is niet slechts een kleine aanpassing; het verandert fundamenteel hoe het model de drijvende kracht berekent die een scheur voortstuwt.

Het team testte hun nieuwe model tegen twee zeer verschillende scenario's om te zien of het standhield. Eerst simuleerden ze een klassieke benchmark voor hydraulisch breken waarbij vloeistof in een gesteente wordt geïnjecteerd om een scheur te creëren. Ze vergeleken hun resultaten met een bekende analytische oplossing, een precieze wiskundige uitkomst afgeleid voor ideale omstandigheden. Het nieuwe model kwam veel dichter bij deze oplossing dan de oudere, gemengde formuleringen, en voorspelde correct de druk die nodig is om de breuk te starten en de breedte van de scheur terwijl deze groeide. De oudere modellen hadden de neiging om de benodigde vloeistofdruk te overschatten, een gebrek dat zou kunnen leiden tot inefficiënte of onveilige booroperaties. In de tweede test simuleerden ze een gesteenteprofiel dat van twee kanten wordt samengedrukt, een conditie die de diepe ondergrondse omgeving nabootst waar gesteenten onder enorme druk staan. Hierbij brak het gesteente niet simpelweg; het vervormde plastisch, waarbij een afschuivingsband ontstond waarin het materiaal langs zichzelf gleed voordat het brak. Het nieuwe model reproduceerde dit complexe gedrag succesvol, waarbij zowel de breuken door afschuiving (veroorzaakt door mechanische compressie) als de breuken door trek (gedreven door vloeistofinjectie) werden gevangen.

De resultaten van deze simulaties laten zien dat de plastische vervorming van het gesteente fungeert als een significante barrière voor de groei van de breuk. Wanneer het gesteente in staat is om permanent te vervormen, absorbeert het meer energie, wat betekent dat de vloeistofdruk binnen de scheur hoger moet zijn om de breuk open en groeiend te houden. Deze extra weerstand verandert het hele profiel van de breuk, waardoor deze smaller wordt en meer kracht vereist om te voortplanten dan een puur elastisch model zou suggereren. Bovendien toonde het model in de compressietests aan hoe de poriedruk binnen het gesteente verandert terwijl de breuk vormt. Naarmate het gesteente vervormt en ruimte creëert voor de vloeistof om te bewegen, daalt de druk lokaal, waardoor er meer vloeistof naar de beschadigde zone wordt getrokken. Deze interactie creëert een feedbackloop waarbij de vloeistofstroom en het mechanische falen van het gesteente nauw aan elkaar gekoppeld zijn, een dynamiek die het nieuwe model met hoge getrouwheid vastlegt.

Door de wiskundige basis van deze modellen te corrigeren, hebben de onderzoekers een betrouwbaarder instrument ontwikkeld voor ingenieurs en wetenschappers die werken met poreuze media. De nieuwe formulering zorgt ervoor dat de overgang tussen verschillende soorten gesteentefalen vloeiend en fysiek realistisch is, waardoor de kunstmatige sprongen die eerdere pogingen teisterden, worden geëlimineerd. Deze nauwkeurigheid is essentieel voor toepassingen variërend van het optimaliseren van olie- en gasproductie tot het ontwerpen van veilige ondergrondse opslag voor koolstofopslag. De studie toont aan dat zelfs in complexe systemen waar vaste stoffen en vloeistoffen interageren, een zorgvuldige, consistente benadering van de onderliggende fysica de ware aard van hoe materialen breken kan onthullen. Het model staat nu klaar om de bewegingen van de aardkorst met meer vertrouwen te voorspellen, waardoor de kloof tussen theoretische mechanica en de rommelige, plastische realiteit van de grond onder onze voeten wordt overbrugd.