ZO Least Squares Migration Lab
OBJECTIVE: Least squares migrate synthetic data with a MATLAB script.
PROCEDURE:
- Load into your working directory MATLAB scripts for the main
program
(testforw.m),
diffraction stack modeling
(forw.m),
Ricker wavelet generator (ricker.m).
migration program
(mig.m),
and veltime model generator
(TIME.m).
Name each file by their names given above.
- In MATLAB, type "testforw" to generate and migrate the shot gather synthetic data for
a point scatterer model.
- Turn code into steepest descent least squares migration code.
That is, the steepest dscent formula is given by m'=m-alpha L^Tdelta d,
where m is the new migration image, m is the old image, delta d
is the residual data d^predicted - d^observed, L^T is the migration operator and
alpha is the step length. In practice the step length is an empirically
derived fraction of the analytic step length. Here we can try with the analytical optimal step
length (steplength.m) by adding this step length calculation into the iterations.
Also turn the code into conjugate gradient least square migartion code and apply optimal step length in the iterations.
- Same as previous question, except create an LSM conjugate gradient code.
- (Optional) Migrate the ZO radar data in radar lab.