Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • harvard-cite-them-right
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Finite-difference modelling of 2-D wave propagation in the vicinity of dipping interfaces: a comparison of anti-aliasing and equivalent medium approaches
ETH Zürich, Institute of Geophysics, CH-8092 Zürich, Switzerland.
University of Borås, Faculty of Textiles, Engineering and Business. Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund, Sweden.ORCID iD: 0000-0002-0905-6188
ETH Zürich, Institute of Geophysics, CH-8092 Zürich, Switzerland.
2022 (English)In: Geophysical Journal International, ISSN 0956-540X, E-ISSN 1365-246X, Vol. 229, no 1, p. 70-96Article in journal (Refereed) Published
Abstract [en]

Finite-difference (FD) modelling of seismic waves in the vicinity of dipping interfaces gives rise to artefacts. Examples are phase and amplitude errors, as well as staircase diffractions. Such errors can be reduced in two general ways. In the first approach, the interface can be anti-aliased (i.e. with an anti-aliased step-function, or a lowpass filter). Alternatively, the interface may be replaced with an equivalent medium (i.e. using Schoenberg & Muir (SM) calculus or orthorhombic averaging). We test these strategies in acoustic, elastic isotropic, and elastic anisotropic settings. Computed FD solutions are compared to analytical solutions. We find that in acoustic media, anti-aliasing methods lead to the smallest errors. Conversely, in elastic media, the SM calculus provides the best accuracy. The downside of the SM calculus is that it requires an anisotropic FD solver even to model an interface between two isotropic materials. As a result, the computational cost increases compared to when using isotropic FD solvers. However, since coarser grid spacings can be used to represent the dipping interfaces, the two effects (an expensive FD solver on a coarser FD grid) equal out. Hence, the SM calculus can provide an efficient means to reduce errors, also in elastic isotropic media.

Place, publisher, year, edition, pages
2022. Vol. 229, no 1, p. 70-96
Keywords [en]
Numerical modelling, Computational seismology, Wave propagation, NONPERIODIC HOMOGENIZATION, NUMERICAL-SIMULATION, HETEROGENEOUS MEDIA, FORM INVERSION, LEBEDEV SCHEME, DISPERSION, REPRESENTATION
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:hb:diva-27446DOI: 10.1093/gji/ggab444ISI: 000743517400005Scopus ID: 2-s2.0-85130559846OAI: oai:DiVA.org:hb-27446DiVA, id: diva2:1635510
Funder
Swedish Research Council, 2019-04878Swedish Nutrition Foundation (SNF), 2-77220-15Available from: 2022-02-07 Created: 2022-02-07 Last updated: 2023-02-06Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

Wittsten, Jens

Search in DiVA

By author/editor
Wittsten, Jens
By organisation
Faculty of Textiles, Engineering and Business
In the same journal
Geophysical Journal International
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 63 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • harvard-cite-them-right
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf