Single frequency 2D acoustic full waveform inversion |

where is the velocity model, is the computed data, and is the observed data. is computed as:

where 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:

We then separate the model into a background and a perturbation:

where is the background component, which is the current model in slowness squared units, and is the perturbation component. After this separation, we can use Taylor expansion on the data around the background component as follows:

By neglecting the higher-order terms in the data series, we can define the linearized modeling operator as:

The first order Born approximation can be used to define the operator:

where the Green's functions now satisfy the acoustic wave equation as follows:

We can now compute the model gradient as follows:

Finally, we can update the model with the gradient:

where is the step size. To estimate the step size, we first evaluate the objective function with the gradient scaled to have a maximum of 2% and 4% of the minimum value of the current model. Using these two points as well as the objective function value at the current model, which is already computed in the gradient calculation, we fit a parabola. If the parabola has positive-side minimum, i.e. both the curvature and the x-axis shift are positive, a new objective function evaluation is performed at the parabola minimum. Then, the two or three evaluations are compared and the scale that resulted in the smallest objective function is used as the step size given that the objective function decreases. Otherwise, the line search is repeated after shrinking the gradient by a factor of 4. The optimization scheme is implemented on the CPU in Fourier domain notations using frequency domain Green's function solutions computed on the GPU.

Single frequency 2D acoustic full waveform inversion |

2012-05-10