Figure 1 a) The true velocity model, b) the initial velocity model, c) multiscale full wave inversion (FWI) tomogram, and d) multiscale phase inversion (MPI) tomogram.
Objective:
Learn how to do multiscale phase inversion (MPI) on local workstation. Read Lei Fu's MPI paper.
Load PhaInv_SeisF90_Demo.zip, from which extract all the files to your working directory. It includes following directories: src for the source codes (Fortran), bin for executable files, model for the velocity models and acquisition coordinates, results for the observed data and inversion results, working for parameter files for the inversion, visual for visualization matlab scripts.
cd to the src directory, create executable files.
modify the directories (eg. DIR_CORE) in Makefile.config, make sure these directories are consistent with the working directory;
cd to folder core, Make clean, then Make;
cd to folder apps, Make clean, then Make; which will make executable files.
cd to the model directory, create true (synthetic case) and initial velocity models, acquisition coordinates.
cd to folder model → 2D_models → Marmous_51x201x10;
modify parameters in mkcoord.m, and then run it in matlab, the velocity (vel.bin, vsmooth.bin) and coordinate (coord_csg.dat) files are created.
cd to the working → 2D_models directory, obtain observed data, run inversions.
modify the directories (DIR_CURRENT) in Makefile.config,
cd to Marmous_51x201x10, and modify parameters in parfile_a2dmod.sh, parfile_fwi_ms.sh, parfile_phainv.sh,
cd.., make a2dmod to obtain the observed data,
make a2dfwi_ms to run the multiscale fwi,
make a2dPhaInv_Field to run the multiscale phase inversion.
Visualization: cd to the visual directory, then run result_plot.m in matlab, we can see the results in Figure 1.