Deformation Due to a Pair of Infinite and a Pair of infinite-finite planar faults in Fractional Standard Linear Solid

Technische Samenvatting: Deformatie door een paar oneindige en een paar oneindig-eindige planaire breuken in een Fractionele Standaard Lineaire Vloeistof

Probleemstelling
Seismisch actieve regio's bevatten vaak meerdere interagerende breuken waarbij slip langs één breuk de stressaccumulatie en -verdeling in naburige breuken aanzienlijk beïnvloedt. Hoewel eerdere studies de interactie tussen breuken hebben gemodelleerd met behulp van integer-orde afgeleiden binnen visco-elastische media (bijv. Maxwell- of Standaard Lineaire Vloeistofmodellen), negeren deze traditionele benaderingen vaak de "geheugeneffecten" die inherent zijn aan geologische materialen. Dit artikel behandelt de deformatie van het lithosfeer-asthenosfeer-systeem van de Aarde tijdens de aseismische periode door twee specifieke scenario's te modelleren van interagerende kruipende strike-slip breuken:

1. Geval-I: De interactie tussen twee oneindige begraven breuken.
2. Geval-II: De interactie tussen één oneindige begraven breuk en één eindige oppervlakte-doorbrekende breuk.

De studie stelt dat het medium van de Aarde zich gedraagt als een visco-elastische halve ruimte gekenmerkt door een Fractionele Standaard Lineaire Vloeistof (FSLS) model, wat de tijdafhankelijke rheologische gedrag en geheugeneffecten van crustale gesteenten beter vastlegt dan integer-orde modellen.

Methodologie
De auteurs hebben theoretische wiskundige modellen ontwikkeld om analytische oplossingen af te leiden voor verplaatsing, spanning en rekcomponenten. De methodologie omvat de volgende belangrijke stappen:

• Wiskundige Formulering: Het medium wordt gedefinieerd met een fractionele standaard lineaire vloeistof constitutieve vergelijking die gebruikmaakt van de Caputo fractionele afgeleide operator (${}_{0}D^{\alpha}_{t}$) waarbij de orde $\alpha \in (0, 1]$. De spanning-rek relaties omvatten effectieve stijfheid ($\mu$) en effectieve viscositeit ($\eta$).
• Beherende Vergelijkingen: De analyse gaat uit van quasi-statische deformatie (het verwaarlozen van inertiekrachten). De spanningen van de bewegingsvergelijkingen en randvoorwaarden worden vastgesteld voor zowel oneindige als eindige breukgeometrieën.
• Oplossingstechniek: De auteurs maken gebruik van de Laplace-transformatie om de tijdafhankelijke fractionele differentiaalvergelijkingen te behandelen. De oplossingen worden afgeleid met behulp van:
  • Het Correspondentieprincipe (Rybicki, 1971).
  • Een Gemodificeerde Green-functie techniek (Maruyama, 1966) om de dislocatiecondities over de breukvlakken te verklaren.
  • De Mittag-Leffler functie ($E_{\alpha}$), die van nature voortkomt uit de inverse Laplace-transformatie van de fractionele operatoren.
• Geanalyseerde Scenario's:
  • De oplossingen worden afgeleid voor de staat vóór breukkruip (accumulatie van tektonische spanning) en na breukkruip (dislocatie door slip).
  • Kruipsnelheden ($V_1, V_2$) en diepteafhankelijke kruipfuncties worden toegepast om de beweging van breuken $F_1$ en $F_2$ te simuleren.
  • Numerieke berekeningen worden uitgevoerd met parameters afgeleid van geologische gegevens (bijv. stijfheid $\mu = 3,5 \times 10^{10}$ N/m², viscositeit $\eta = 5 \times 10^{19}$ Pa·s) om grafische resultaten te genereren voor verschillende breukhellingen ($\theta$), fractionele orden ($\alpha$) en kruipsnelheden.

Belangrijkste Bijdragen

1. Fractionele Modellering van Breukinteractie: Dit artikel breidt bestaande breukinteractiemodellen uit door fractionele calculus te integreren, wat een nauwkeurigere representatie mogelijk maakt van het geheugenafhankelijke visco-elastische gedrag van de lithosfeer.
2. Analytische Oplossingen voor Complexe Geometrieën: De studie biedt gesloten vorm analytische expressies voor verplaatsings-, spanning- en rekvelden voor twee verschillende configuraties: een paar oneindige begraven breuken en een gemengd systeem van oneindige en eindige oppervlakte-doorbrekende breuken.
3. Kwantificering van Geheugeneffecten: Het onderzoek toont expliciet aan hoe de fractionele orde $\alpha$ de accumulatie en afgifte van spanning beïnvloedt, waarbij de wiskundige parameter wordt gekoppeld aan het fysieke "geheugen" van het geologische medium.

Resultaten en Numerieke Bevindingen
De numerieke analyse, ondersteund door grafische weergaven, levert verschillende specifieke observaties op:

• Verplaatsingsgrootten: De verandering in oppervlakteverplaatsing nabij de breuken is van de orde van $10^{-3}$ m/jaar. Dit komt goed overeen met waargenomen geofysische gegevens (bijv. USGS-rapporten van 2–6 mm/jaar), terwijl eerdere integer-orde modellen waarden in de orde van $10^{-5}$ voorspelden, wat een discrepantie vertoont met real-world observaties.
• Invloed van Breukgeometrie:
  • In Geval-I (twee oneindige breuken) is de verplaatsingsgrootte nabij de tweede breuk groter dan nabij de eerste.
  • In Geval-II (oneindige en eindige breuken) is de verplaatsingsgrootte nabij de eindige breuk lager dan nabij de oneindige breuk.
  • De verplaatsingspieken bevinden zich nabij de breuklijn en nemen af naarmate de afstand toeneemt, om uiteindelijk naar nul te tenderen.
• Effect van de Fractionele Orde ($\alpha$): De studie vindt dat de snelheid van spanningaccumulatie en -afgifte sterk afhankelijk is van $\alpha$. Lagere waarden van $\alpha$ (wat wijst op sterkere geheugeneffecten) resulteren in hogere magnituden van schuifspanningaccumulatie en een langzamere relaxatie. Naarmate $\alpha$ toeneemt richting 1 (benadering van integer gedrag), neemt de snelheid van spanningsafgifte toe.
• Drijvende Kracht en Spanningsverdeling: De totale drijvende kracht over de breuken vertoont een snelle afname na het begin van de kruip. De magnitude van deze afname varieert met de breukhelling en de kruipsnelheid. Contourkaarten onthullen dat spanningaccumulatie het meest significant is nabij het vrije oppervlak en de breukdislocatiezones, terwijl spanningsreductie prominent aanwezig is op diepten die overeenkomen met de breuklocaties.
• Rek-analyse: Variaties in schuifrek ($e_{12}$) worden gevonden in de orde van $10^{-4}$, wat consistent is met gegevens van real-world gebeurtenissen zoals de 2011 Tohoku-oki aardbeving.

Betekenis en Claims
De auteurs beweren dat hun fractionele model een algemenere en nauwkeurigere representatie biedt van de dynamiek van breukinteractie vergeleken met traditionele integer-orde modellen. Specifiek:

• De fractionele benadering slaagt er succesvol in om de geheugeneffecten en niet-lineariteit van visco-elastische media te vangen, die vaak worden genegeerd in standaardmodellen.
• De voorspellingen van het model voor verplaatsingssnelheden ($10^{-3}$) tonen een sterke overeenstemming met geobserveerde geofysische gegevens, wat suggereert dat fractionele calculus een superieur instrument is voor het voorspellen van seismisch gedrag en substantiële deformatie.
• De bevindingen geven aan dat de interactie tussen breuken de spanningsvelden uitgebreider beïnvloedt dan een enkele breuk, en dat de fractionele orde $\alpha$ een kritieke parameter is bij het bepalen van de snelheid van spanningsafgifte.

Het artikel concludeert dat deze theoretische modellen het begrip van ondergrondse deformatie en de rol daarvan bij het beïnvloeden van breukdynamiek en aardbevinggeneratie verbeteren, en een robuust kader bieden voor het analyseren van aseismische perioden in seismisch actieve zones.