TY - JOUR
T1 - A tool for designing digital filters for the Hankel and Fourier transforms in potential, diffusive, and wavefield modeling
AU - Werthmüller, Dieter
AU - Key, Kerry
AU - Slob, Evert C.
N1 - Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public.
PY - 2019
Y1 - 2019
N2 - The open-source code fdesign makes it possible to design digital linear filters for the Hankel and Fourier transforms used in potential, diffusive, and wavefield modeling. Digital filters can be derived for any electromagnetic (EM) method, such as methods in the diffusive limits (direct current, controlledsource EM [CSEM]) as well as methods using higher frequency content (ground-penetrating radar [GPR], acoustic and elastic wavefields). The direct matrix inversion method is used for the derivation of the filter values, and a brute-force minimization search is carried out over the defined spacing and shifting values of the filter basis. Included or user-provided theoretical transform pairs are used for the inversion. Alternatively, one can provide layered subsurface models that will be computed with a precise quadrature method using the EM modeler empymod to generate numerical transform pairs. The comparison of the presented 201 pt filter with previously presented filters indicates that it performs better for some standard CSEM cases. The derivation of a longer 2001 pt filter for a GPR example with a 250 MHz center frequency proves that the filter method works not only for diffusive EM fields but also for wave phenomena. The presented algorithm provides a tool to create problem specific digital filters. Such purpose-built filters can be made shorter and can speed up consecutive potential, diffusive, and wavefield inversions.
AB - The open-source code fdesign makes it possible to design digital linear filters for the Hankel and Fourier transforms used in potential, diffusive, and wavefield modeling. Digital filters can be derived for any electromagnetic (EM) method, such as methods in the diffusive limits (direct current, controlledsource EM [CSEM]) as well as methods using higher frequency content (ground-penetrating radar [GPR], acoustic and elastic wavefields). The direct matrix inversion method is used for the derivation of the filter values, and a brute-force minimization search is carried out over the defined spacing and shifting values of the filter basis. Included or user-provided theoretical transform pairs are used for the inversion. Alternatively, one can provide layered subsurface models that will be computed with a precise quadrature method using the EM modeler empymod to generate numerical transform pairs. The comparison of the presented 201 pt filter with previously presented filters indicates that it performs better for some standard CSEM cases. The derivation of a longer 2001 pt filter for a GPR example with a 250 MHz center frequency proves that the filter method works not only for diffusive EM fields but also for wave phenomena. The presented algorithm provides a tool to create problem specific digital filters. Such purpose-built filters can be made shorter and can speed up consecutive potential, diffusive, and wavefield inversions.
UR - http://www.scopus.com/inward/record.url?scp=85062694650&partnerID=8YFLogxK
U2 - 10.1190/GEO2018-0069.1
DO - 10.1190/GEO2018-0069.1
M3 - Article
AN - SCOPUS:85062694650
SN - 0016-8033
VL - 84
SP - F47-F56
JO - Geophysics
JF - Geophysics
IS - 2
ER -