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

De grond onder onze voeten is zelden een solide, ononderbroken blok. Op veel plaatsen, van diepe olievoorraden tot de aardkorst zelf, is het gesteente doorkruist door breuken—dunne, vlakke scheuren die kilometers lang kunnen zijn. Wanneer energie door dit gebroken landschap reist, of het nu gaat om een seismische golf van een aardbeving of een gelft die wordt gebruikt om een materiaal te inspecteren, veranderen de scheuren het verhaal. Ze verstrooien de energie, verzwakken het signaal en zorgen er soms zelfs voor dat de rotswanden tegen elkaar aan schuren of uit elkaar glijden. Het begrijpen van de exacte manier waarop deze golven met de scheuren interageren, is essentieel voor geologen die proberen te in kaart te brengen wat er onder de grond ligt, en voor ingenieurs die de veiligheid van constructies waarborgen. De fysica van deze interactie is echter berucht moeilijk vast te leggen. Een eenvoudige scheur is niet slechts een gat; het is een complexe grens waar rotswanden tegen elkaar aan kunnen drukken, met wrijving kunnen glijden of open en dicht kunnen klappen, terwijl de golf er in een fractie van een seconde doorheen trekt.

Decennialang hebben wetenschappers vertrouwd op vereenvoudigde modellen om deze gedragingen te voorspellen, waarbij ze breuken vaak behandelden alsof ze gemaakt waren van onzichtbare veren die uitrekken en samendrukken. Hoewel deze modellen nuttig zijn, schieten ze tekort wanneer de golven sterk genoeg zijn om de rotswanden volledig tegen elkaar te laten drukken, te laten glijden of te laten scheiden. Ze kunnen de 'stick-slip'-beweging niet vatten waarbij wrijving de rots op zijn plaats houdt totdat de kracht te groot wordt, noch kunnen ze de nietlineaire verstijving nauwkeurig beschrijven die optreedt wanneer een scheur dichtgeknepen wordt. Een nieuwe studie door onderzoekers van de Universiteit van Bergen, het Noorse Onderzoekscentrum en de Technische Universiteit van München pakt deze hiaten aan door een uitgebreid computationeel kader te bouwen dat elastische golven die door gebroken gesteente bewegen met ongekende realiteitszin simuleert. In plaats van te vertrouwen op vereenvoudigde aannames, heeft het team een methode ontwikkeld die vier verschillende niveaus van complexiteit kan afhandelen, variërend van basisgedrag als veren tot volledige contactmechanica waarbij wrijving en het fysieke openen en sluiten van scheuren in realtime worden berekend.

De onderzoekers construeerden een verenigde digitale omgeving waarin het gesteente en de breuken worden behandeld als één enkel, onderling verbonden systeem. In hun simulaties wordt het gesteente gerepresenteerd als een driedimensionaal volume, terwijl de breuken worden gemodelleerd als tweedimensionale oppervlakken die erin zijn ingebed. Dit stelt de computer in staat om bij te houden hoe een golf die door het gesteente reist een scheur raakt, hoe de scheur vervormt en hoe die vervorming weer rimpelingen terug de rots in stuurt. Het team testte vier specifieke manieren waarop de scheuren zich konden gedragen. De eerste twee modellen gebruikten veerachtige wetten, waarbij de kracht over de scheur proportioneel is aan de mate waarin het gesteente beweegt, met één versie die een eenvoudige lineaire relatie toestaat en een andere die een complexere, nietlineaire regel gebruikt die rekening houdt met het feit dat scheuren moeilijker samen te drukken zijn naarmate ze sluiten. De andere twee modellen introduceerden de realiteit van wrijving en contact. In deze scenario's mogen de scheurvlakken tegen elkaar aan raken en langs elkaar glijden, gestuurd door wrijvingswetten die vergelijkbaar zijn met die welke voorkomen dat autobanden slippen, terwijl ook een regel wordt afgedwongen dat de rotswanden niet door elkaar heen kunnen gaan.

Om de nauwkeurigheid van hun nieuwe kader te waarborgen, voerde het team een reeks rigoureuze tests uit. Ze simuleerden golven die een enkele scheur raken en vergeleken hun computeresultaten met bekende wiskundige oplossingen voor eenvoudige gevallen. Ze ontdekten dat hun methode het gedrag van de golf met hoge precisie kon voorspellen, waarbij de theoretische verwachtingen voor zowel de snelheid van de golf als de hoeveelheid energie die erdoorheen ging of terugkaatste, werden benaderd. Vervolgens duwden ze het systeem verder door golven te simuleren die sterk genoeg waren om de scheuren te laten openen, plakken en glijden. In deze complexere scenario's, waar de scheurvlakken met een wrijvingscoëfficiënt van 1,0 tegen elkaar aan wreven, lieten de simulaties zien dat de methode stabiel bleef en consistente resultaten produceerde. De onderzoekers verifieerden dat hun benadering de plotselinge veranderingen in de fysica kan afhandelen die optreden wanneer een scheur overgaat van plakken naar glijden, een overgang die eenvoudigere modellen vaak in de war brengt.

Wanneer het team de vier verschillende modellen naast elkaar legde, waren de verschillen in de resulterende golfpatronen opvallend. De simulaties toonden aan dat de keuze van het model er aanzienlijk toe doet. Wanneer een scheur mocht sluiten en verstijven volgens de nietlineaire wet, transmitteerde deze veel meer golfenergie dan een scheur die gemodelleerd werd met een eenvoudige lineaire veer. Dit komt omdat het nietlineaire model rekening houdt met het feit dat naarmate de scheur sluit, de rotswanden harder tegen elkaar aan drukken, waardoor de scheur effectief stijver en transparanter voor de golf wordt. Bovendien veranderde de inclusie van wrijving de vorm van de wavefronten. In de wrijvingsmodellen was de golf die aan de andere kant van de scheur tevoorschijn kwam breder en meer vervormd dan in de wrijvingsloze versies. De onderzoekers observeerden dat de wrijving ervoor zorgde dat de scheurvlakken gedurende langere perioden aan elkaar bleven plakken, wat de vrijgave van energie vertraagde en de manier waarop de golf verstrooid werd, veranderde.

De studie stopte niet bij enkele scheuren. De onderzoekers pasten hun kader toe op een driedimensionaal domein met zes snijdende breuken, waardoor een klein netwerk van scheuren met verschillende eigenschappen ontstond. In deze complexe omgeving vochten ze verschillende maximale sluitingslimieten toe aan verschillende sets breuken, wat betekende dat sommige scheuren strakker konden sluiten dan andere. De simulatie toonde aan dat het golfveld zeer heterogeen werd. Breuken die stijver waren en nauwer konden sluiten, reflecteerden meer energie en transmitteerden minder, wat leidde tot duidelijke schaduwen in het golfpatroon. De onderzoekers waren in staat om exact bij te houden hoeveel elke scheur sloot en hoeveel hij gleed, wat onthulde dat zelfs kleine verschillen in de fysieke eigenschappen van de scheuren leidden tot merkbaar verschillende vervormingsreacties. Bijvoorbeeld, een scheur die de mogelijkheid had om strakker te sluiten, vertoonde meer 'plakgedrag', waarbij de vlakken langer bij elkaar werden gehouden voordat ze gingen glijden, wat op zijn beurt de golf die erdoorheen trok, beïnvloedde.

De implicaties van deze bevindingen reiken verder dan de directe resultaten. Door aan te tonen dat één enkel computationeel kader alles kan afhandelen, van eenvoudig veerachtig gedrag tot complexe wrijvingscontacten, hebben de onderzoekers een instrument geleverd dat gebruikt kan worden om real-world seismische data met grotere getrouwheid te interpreteren. Het vermogen om te simuleren hoe golven verstrooien en dempen in aanwezigheid van wrijvende, openende en sluitende scheuren betekent dat geowetenschappers beter onderscheid kunnen maken tussen verschillende soorten gesteentevormaties en vloeistofgevulde breuken. De studie bevestigt dat het negeren van de nietlineaire en wrijvingsaspecten van breukmechanica kan leiden tot significante fouten in het voorspellen van hoe energie door de ondergrond beweegt. Hoewel het huidige werk zich richt op puur mechanische interacties, is het kader gebouwd op een fundament dat compatibel is met modellen voor doorstroming van vloeistoffen, wat wijst op een pad naar nog realistischere simulaties waarbij water of olie die door de scheuren beweegt, interageert met de passerende golven. Het werk staat als een demonstratie van het feit dat door de volledige complexiteit van hoe gesteentesscheuren zich gedragen te omarmen, we een helderder beeld kunnen vormen van de verborgen wereld onder onze voeten.