sfzomig3 performs 3-D zero-offset modeling or migration using one-way wave extrapolation and the exploding reflector concept.

The following example from rsf/su/rsflab10 shows the result of migrating the benchmark Viking Graben dataset.

The algorithm used by sfzomig3 is known as extended split-step or phase-shift plus interpolation (PSPI). It works in the frequency domain and mixed space-wavenumber domain.

  • Gazdag, J., and Sguazzero, P., 1984, Migration of seismic data by phase shift plus interpolation: Geophysics, 49, 124-131.
  • Kessinger, W., 1992, Extended split-step Fourier migration: 62nd Ann.
    Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 917-920.

Depending on inv= flag, sfzomig3 performs modeling or migration (if mode=m) or datuming (upward or downward) wavefield continuation (if mode=d). The default is migration (mode=m inv=n).

The algorithm efficiency depends on the number of reference velocities, which corresponds to the number of spatial Fourier transforms. By default, the algorithm is trying to estimate this number automatically at each depth step. The maximum number can be set by nrmax=.

The pmx=, pmy=, and tmx=, tmy= parameters control padding and tapering in space, needed to avoid boundary reflections and Fourier wrap-around artifacts.

10 previous programs of the month: