Computational analysis of extended full waveform inversion |

where is the modeled data, is the velocity model, is the source function, is frequency, and are the source and receiver coordinates, and is the model coordinate. In the acoustic, constant-density case, the Green's function satisfies:

The propagation can be done in the time domain by convolving each model point with a finite-difference stencil. However, the time marching requires the time axis sampling to satisfy dispersion and stability conditions (Marfurt, 1984), generally much finer than the data sampling. Moreover, each time step requires multiplying the time slice by the velocity squared. Therefore, the cost of forward modeling can be written as:

where , and , are the number of points along the three spatial axes, is the number of sources, is the cost of convolving one model location by the time-domain finite-difference stencil and is the number of time samples for propagation. By linearizing equation 1 over the squared slowness, we can compute the adjoint as:

where is the perturbation in squared slowness and denotes the complex conjugate. For the adjoint, the imaging time sampling can be much larger than that of propagation since it does not need to satisfy the dispersion and stability conditions. Hence, the cost of computing the adjoint of FWI can be written as:

where is the number of time samples for imaging. The total cost of one iteration of FWI becomes

Computational analysis of extended full waveform inversion |

2012-10-29