Create and apply a source wavelet correction filter for full-waveform
inversion of field data.
Design a filter which minimizes the misfit to the field recordings
when being applied to the raw synthetics. Other criteria are possible.
Field recordings and raw synthetics for one shot.
Optionally additional time series can be passed to be convolved with the
determined source time function.
All time series have identical temporal sampling.
Recorded data and corresponding synthetics must have identical
spatial and temporal sampling. All
recorded traces must have passed the same filters. Raw synthetics have
larger bandwidth than recordings (otherwise just swap synthetics and
Impulse response of the designed filter and synthetics
convolved with filter response.
If additional series were passed, they are convolved with the newly determined
filter and are returned to the user.
The directory src/libs/libstfinv contains the source code for a library which supports the estimation of a source time function in waveform inversion.
The library is designed to provide a common API (application programming
interface) to the library user in C++ as well as in C.
A Fortran API is not yet implemented but could be constructed on top of the C
The library provides several approaches to finding an optimized source wavelet
Each of the approaches is addressed through the same API, such that programs
using the library could immediatly benefit from an new approach after its
implementation without need to modify consumer programs.
The approaches are encapsulated in so-called 'engines'.
coming within the source code for further details.