Waveform Reconstruction Inversion
This software release includes a demonstration of the parallel Wavefield reconstruction Imaging algorithm in Matlab that has been developed at SLIM.
Author: Bas Peters Date : September, 2014
Contents
Downloading & Dependencies
The code can be found in the SLIM sofware release under /applications/WaveformInversion/2DWRI/.
The code has been tested with Matlab R2013a and requires the Parallel Computing Toolbox.
This code uses the following packages, also found in the tools part of the SLIM software release.
- utilities/SPOT - object oriented framework for matrix-free linear algebra.
- utilities/pSPOT - parallel extension of SPOT.
- algorithms/2DFreqModeling - 2D constant-density acoustic modeling
- algorithms/WRI - Functions for the Wavefield Reconstruction Inversion algorithm
- operators/misc - Misc. operators such as interpolation, smoothing and splines
- functions/misc - Misc. functions.
- tools/solvers/QuasiNewton/minConf_mod/ - nonlinear optimization
Running & Parallelism
All the examples and results are produced by the scripts found in this software release under /applications/WaveformInversion/2DWRI/. Start matlab from that directory or run startup in that directory to add the appropriate paths.
To run the scripts follow the instrictions in the README file enclosed with the code.
Read more about SLIM's approach to parallel computing in Matlab.
Examples and results
References
[1] Tristan van Leeuwen, Felix J. Herrmann, Geophysical Journal International,2013. Mitigating local minima in full-waveform inversion by expanding the search space.
[2] Tristan van Leeuwen, Felix J. Herrmann. 2013. A penalty method for PDE-constrained optimization.
[3] Bas Peters, Felix J. Herrmann, Tristan van Leeuwen. EAGE, 2014. Wave-equation based inversion with the penalty method: adjoint-state versus wavefield-reconstruction inversion.