Contents
Constrained FWI
Illustrates the use of constrained FWI in time-domain and frequency-domain. Example of constraints for a simple 2D model in a well-to-well setting. Receivers in one well, sources in the other.
The examples can be run in serial mode, however you have to make sure matlab won't automatically open a parallel pool. Use 2 workers with parpool(2) if you can.
% Author: Bas Peters and Mathias Louboutin % % April 2016 % Seismic Laboratory for Imaging and Modeling % Department of Earch & Ocean Sciences % The University of British Columbia
Dependencies
The modeling code uses the following packages, found in the tools part of the software release.
- SPOT - object oriented framework for matrix-free linear algebra.
- pSPOT - parallel extension of SPOT.
Running & Parallelism
All the examples can be reproduced by running the scripts found in the software release under applications/WaveformInversion/ConstrainedFWI/FreqDomain and applications/WaveformInversion/ConstrainedFWI/TimeDomain. Start matlab from that directory or run startup in that directory to add the appropriate paths.
The scripts can be run in serial mode but parallel mode is advised for the modeling and imaging examples. Use parpool to start parallel pool with the appropriate configuration (2 or 4 workers should be enough for both examples).
Functions
Depending on what modelling kernel you want to use (Frequency domain or Time domain) you will need different packages set up in the startup.m file in each example folder.
Examples
A few examples are included here
- An examples in time domain is shown in Example_Time_Domain.m.
- An examples in Frequency domain is shown in exp_2D_freqFWI_simple_constrained.m.
References
https://www.slim.eos.ubc.ca/Publications/Public/TechReport/2015/peters2015EAGErwi/peters2015EAGErwi.html.[1] Bas Peters*, Brendan Smithyman, Felix J. Herrmann, Regularizing waveform inversion by projection onto intersections of convex sets Technical report