A Model Order Reduction Method for Seismic Applications Using the Laplace Transform

Titel: Een Model Order Reduction-methode voor Seismische Toepassingen met behulp van de Laplace-transformatie

Auteurs: Fernando Henríquez en Matthias Schlottbom
Datum: Februari 2026 (voorgestelde versie)

1. Probleemstelling en Motivatie

De numerieke modellering van seismische verschijnselen, zoals seismische tomografie en full waveform inversion (FWI), vereist het oplossen van golfvergelijkingen (akoestisch of elastisch) voor duizenden of miljoenen queries. Deze herhaalde berekeningen zijn computationally zeer intensief en maken real-time simulaties vaak onmogelijk.

De bronterm in deze modellen wordt vaak gemodelleerd met een Ricker-golfje (een tijdsafhankelijke functie die de spectrale eigenschappen van seismische bronnen nabootst). De uitdaging ligt in het versnellen van de oplossing van deze evolutieproblemen zonder de nauwkeurigheid te verliezen. Traditionele Model Order Reduction (MOR) methoden, zoals Reduced Basis (RB) methoden, stuiten bij golfvergelijkingen vaak op het "Kolmogorov-barrière"-probleem: de oplossingsmanifold heeft een trage afname van de N-width, wat betekent dat zeer grote reductieruimtes nodig zijn voor hoge nauwkeurigheid.

2. Methodologie: LT-MOR (Laplace Transform Model Order Reduction)

De auteurs stellen een nieuwe MOR-strategie voor die werkt in het Laplace-domein in plaats van direct in de tijdsdomein. De kern van de methode is als volgt:

• Transformatie: De tijdsafhankelijke golfvergelijking (met homogene beginvoorwaarden en een bronterm bestaande uit een Ricker-golfje vermenigvuldigd met een ruimtelijk profiel) wordt getransformeerd naar het Laplace-domein. Dit resulteert in een familie van tijds-onafhankelijke, stationaire problemen (vergelijkbaar met een gedempte Helmholtz-vergelijking) die afhankelijk zijn van de complexe Laplace-variabele $s$.
• Parametrisatie: Hoewel het oorspronkelijke probleem niet-parametrisch is, fungeert de Laplace-variabele $s$ als een parameter in het getransformeerde domein. Dit maakt het mogelijk om gevestigde RB-technieken voor parametrische PDE's toe te passen.
• Offline-fase:
  1. Er wordt een set steunpunten (samples) geselecteerd in het Laplace-domein, specifiek langs een lijn $s = \mu + i\tau$.
  2. Voor elk steunpunt wordt de bijbehorende stationaire probleem opgelost met een hoge precisie (High-Fidelity, HF) discretisatie (bijv. eindige elementen).
  3. De oplossingen (snapshots) worden verzameld en gecomprimeerd via Proper Orthogonal Decomposition (POD) om een laag-dimensionale reductieruimte $V_R^{(rb)}$ te construeren.
  4. De auteurs gebruiken sinc-kwadratuur voor de selectie van de steunpunten, wat zorgt voor een exponentiële convergentie van de discretisatiefout in het Laplace-domein.
• Online-fase:
  1. De oorspronkelijke tijdsafhankelijke probleem wordt geprojecteerd op de gereduceerde ruimte.
  2. Dit resulteert in een klein stelsel gewone differentiaalvergelijkingen (ODE's) dat efficiënt in de tijd kan worden geïntegreerd (bijv. met een Newmark-beta schema).
  3. Cruciaal: De tijdsoplossing wordt niet verkregen door een numerieke inversie van de Laplace-transformatie (wat vaak instabiel is), maar door direct te integreren in de gereduceerde ruimte.

3. Belangrijkste Bijdragen

De paper levert zowel theoretische als praktische inzichten:

1. Exponentiële Convergentie: De auteurs bewijzen dat de oplossing van de abstracte golfvergelijking met exponentiële nauwkeurigheid kan worden benaderd door eindig-dimensionale deelruimtes. De convergentiesnelheid is onafhankelijk van de parameters van het Ricker-golfje (piektijd $t_0$ en breedte $\alpha$), zolang de fout groter is dan een intrinsieke nauwkeurigheidsdrempel.
2. Robuustheid: De afgeleide foutgrenzen zijn expliciet en robuust ten opzichte van de golfje-parameters. Dit is een significant voordeel ten opzichte van eerdere methoden die gevoelig kunnen zijn voor deze parameters.
3. Intrinsieke Nauwkeurigheidslimiet: De analyse identificeert een fundamentele limiet in de nauwkeurigheid die wordt bepaald door de waarde van het Ricker-golfje en zijn afgeleide op het tijdstip $t=0$. Deze waarden vervallen exponentieel ($e^{-(\alpha t_0)^2}$), wat betekent dat de fout niet verder kan dalen dan deze drempel, ongeacht de grootte van de reductieruimte.
4. Onafhankelijkheid van Discretisatie: De theoretische resultaten zijn onafhankelijk van de onderliggende Galerkin-discretisatie (ruimtelijke mesh), wat de methode geldig maakt in een oneindig-dimensionale setting.
5. Praktische Implementatie: Er wordt een volledig werkend algoritme gepresenteerd dat gebruikmaakt van sinc-kwadratuur voor steunpuntselectie, wat zorgt voor exponentiële afname van de steekproeffout.

4. Numerieke Resultaten

De methode is getest op de scalare golfvergelijking in een 2D domein:

• Nauwkeurigheid: De relatieve fout tussen de hoge precisie-oplossing en de gereduceerde oplossing neemt exponentieel af naarmate de dimensie $R$ van de reductieruimte toeneemt. De resultaten bevestigen de theoretische voorspellingen.
• Invloed van Samples: De nauwkeurigheid hangt ook exponentieel af van het aantal Laplace-domein samples ($M$). Zodra $M$ groot genoeg is, wordt de fout bepaald door de reductie van de ruimte ($R$).
• Snelheidswinst (Speed-up): De methode toont aanzienlijke snelheidswinsten. Voor de geteste configuraties varieerde de snelheidswinst tussen 6x en 14x ten opzichte van de volledige tijdsstap-methode. De offline-fase (het oplossen van de Laplace-problemen) is paralleliseerbaar, terwijl de online-fase zeer efficiënt is omdat deze slechts werkt met een klein aantal vrijheidsgraden ($R$).
• Visualisatie: De gereduceerde basisfuncties en de oplossing op verschillende tijdstippen tonen een hoge overeenkomst met de referentieoplossing, zelfs met een relatief kleine $R$ (bijv. $R=50$).

5. Significatie en Toekomstperspectief

Deze studie is significant omdat het een brug slaat tussen de theorie van Laplace-transformaties en praktische seismische modellering. Het biedt een stabiel en nauwkeurig alternatief voor traditionele tijdsstap-methoden en parametrische MOR-methoden die worstelen met de Kolmogorov-barrière bij golfproblemen.

Toekomstige uitbreidingen die door de auteurs worden genoemd, omvatten:

• Het incorporeren van extra fysieke of geometrische parameters (bijv. materiaaleigenschappen of bronlocaties) voor een echt parametrisch probleem.
• Het toepassen van de methode op complexere golfsystemen, zoals volledig elastische golfvergelijkingen.
• Het uitbreiden naar andere seismische bronmodellen (bijv. Kupper-golfjes) en het gebruik van adaptieve steekproefstrategieën in het Laplace-domein.

Kortom, de paper presenteert een wiskundig onderbouwde en computationally efficiënte methode om seismische simulaties aanzienlijk te versnellen, met gegarandeerde foutgrenzen die robuust zijn voor de kenmerken van de bron.