Finite Difference Redatuming Lab




Figure 1. Comparison between analytic Green's function and numerical one calculated by the reciprocity theorem of the correlation type.

Objective: Execute program fd1.m to reproduce above figure. This result is produced by forward modeling the acoustic wave equation to solve for the pressure field along the enclosing boundary of the 1-scatterer model, which is a single scatterer embedded in a homogeneous medium. The reciprocity equation is implemented by backprojecting the time-reversed recorded fields at the boundary, and the result is that the Green's function G(A|B) is recovered for A and B somewhere inside the model volume.

Procedure:

  1. Download the following programs: backprop.m, ricker.m, indata2.m, indata1.m, chap3b.m, mkdamp.m, plot1.m, plot2.m, plot4.m, plot5.m, plot6.m, plot7.m, plot9.m, plot33.m, and store1.m.
  2. Type chap3b.m to generate the above picture for a single scatterer model.