Show simple item record

dc.contributor.author
Koene, Erik F. M.
dc.contributor.supervisor
Robertsson, Johan O.A.
dc.contributor.supervisor
Andersson, Fredrik
dc.contributor.supervisor
Arntsen, Børge
dc.contributor.supervisor
Blanch, Joakim O.
dc.date.accessioned
2021-02-24T08:09:01Z
dc.date.available
2021-02-23T10:46:55Z
dc.date.available
2021-02-24T07:49:11Z
dc.date.available
2021-02-24T08:09:01Z
dc.date.issued
2020
dc.identifier.uri
http://hdl.handle.net/20.500.11850/471059
dc.identifier.doi
10.3929/ethz-b-000471059
dc.description.abstract
Seismograms (i.e., recordings of seismic waves that propagate through the earth) can be used to uncover information about the earth's subsurface. Such investigations require accurate numerical wave simulations. One of the most common techniques to carry out these simulations is the finite-difference (FD) method. In the FD method, (1) derivatives are replaced with approximations of limited accuracy, and (2) continuous space and time are discretized into finite steps. The FD method is a fast numerical method, but it also introduces inaccuracies. In this thesis, we propose four procedures that reduce these inaccuracies. The overarching aim is to provide fast FD simulations (using large steps in space and time) while yielding accurate solutions. The first proposed method is the use of a filter pair: the forward and inverse time-dispersion transforms. These transforms must be applied before the simulation (to modify the source wavelet) and after the simulation (to modify the recorded seismic signals). They correct for the inaccuracy induced by the approximation of the temporal derivative in the wave equation. We show that the method applies to acoustic and elastic wave simulations. Furthermore, we show that the method applies to viscoelastic FD simulations if they use standard memory variables. The second proposed method is the use and design of `optimal' FD operators. Such FD operators are highly accurate for a prescribed wavenumber range. We obtain these FD operators using the Remez exchange algorithm, a well-known algorithm in the field of filter design. Our work generalizes the existing literature drastically: (1) we consider arbitrary derivative orders, (2) we consider three cost-functions [the absolute error, the relative error, the group velocity error], (3) we consider arbitrary input locations, (4) we can compute solutions that are optimal in a least-squares or maximum norm sense. Optimal results in FD modeling are obtained with the FD operator designed for the relative error. The third proposed method concerns the implementation of point-sources in FD simulations. These sources are typically modeled by exciting the source on a single FD node. We show that such an implementation leads to wavenumber-varying amplitude errors. In effect, two artifacts are generated: (1) ringing is introduced, and (2) erroneous wave modes may be excited. We show how to correct for this error using a filter in the wavenumber domain. The `FD-consistent' point-source that we propose minimizes these artifacts. The fourth proposed method concerns the use of interface representation schemes in wave simulations. For this, we compare five interface representations from geophysical literature. We find that, in acoustic simulations, optimal results are obtained with anti-aliasing of the fine velocity model. Conversely, in isotropic and anisotropic elastic simulations, optimal results are obtained with the Schoenberg \& Muir (1989) calculus. The proposed methods have two attractive features: (1) they allow FD simulations with large steps in space and time, (2) they must only be applied before and after the simulation to improve the accuracy, and have a negligible computational cost. Hence, they allow for fast FD simulations with a minimal computational cost, while yielding excellent accuracy.
en_US
dc.format
application/pdf
en_US
dc.language.iso
en
en_US
dc.publisher
ETH Zurich
en_US
dc.rights.uri
http://rightsstatements.org/page/InC-NC/1.0/
dc.subject
SEISMIC WAVES/PROPAGATION (GEOPHYSICS)
en_US
dc.subject
GEOPHYSICS
en_US
dc.subject
FINITE DIFFERENCE METHOD (NUMERICAL MATHEMATICS)
en_US
dc.subject
filtering
en_US
dc.subject
anisotropy
en_US
dc.subject
elastic wave propagation
en_US
dc.title
A filtering approach to remove finite-difference errors from wave equation simulations
en_US
dc.type
Doctoral Thesis
dc.rights.license
In Copyright - Non-Commercial Use Permitted
dc.date.published
2021-02-24
ethz.size
267 p.
en_US
ethz.code.ddc
DDC - DDC::5 - Science::550 - Earth sciences
en_US
ethz.notes
This work was supported by SNF grant 2-77220-15.
en_US
ethz.identifier.diss
27218
en_US
ethz.publication.place
Zurich
en_US
ethz.publication.status
published
en_US
ethz.leitzahl
ETH Zürich::00002 - ETH Zürich::00012 - Lehre und Forschung::00007 - Departemente::02330 - Dep. Erdwissenschaften / Dep. of Earth Sciences::02506 - Institut für Geophysik / Institute of Geophysics::03953 - Robertsson, Johan / Robertsson, Johan
en_US
ethz.date.deposited
2021-02-23T10:47:06Z
ethz.source
FORM
ethz.eth
yes
en_US
ethz.availability
Open access
en_US
ethz.rosetta.installDate
2021-02-24T07:49:23Z
ethz.rosetta.lastUpdated
2021-02-24T07:49:23Z
ethz.rosetta.exportRequired
true
ethz.rosetta.versionExported
true
ethz.COinS
ctx_ver=Z39.88-2004&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.atitle=A%20filtering%20approach%20to%20remove%20finite-difference%20errors%20from%20wave%20equation%20simulations&rft.date=2020&rft.au=Koene,%20Erik%20F.%20M.&rft.genre=unknown&rft.btitle=A%20filtering%20approach%20to%20remove%20finite-difference%20errors%20from%20wave%20equation%20simulations
 Search print copy at ETH Library

Files in this item

Thumbnail

Publication type

Show simple item record