Higher accuracy linearizations

As noted earlier, the approximation  

 150#150 (60)

corresponds to an explicit numerical solution to the differential equation (). However, this is neither the only possible solution, nor the most accurate, and furthermore it is only conditionally stable.

We can, however, solve Equation () using other numerical schemes. Two possibilities are implicit numerical solutions, where we approximate  

 151#151 (61)

or bilinear numerical solutions, where we approximate  

 152#152 (62)

Equations () and () are first order, but Equation () is second order accurate as a function of the phase 5#5.Numerical schemes based on Equation () are conditionally stable, but numerical schemes based on Equations () and () are unconditionally stable.

In the context of partial differential equations, the bilinear approximation () is known under the name of Crank-Nicolson and has been extensively used in migration by downward-continuation using the paraxial wave equation (). Figures  and compare the approximations in Equations (), () and () as a function of phase. Both the explicit and implicit solutions lead to errors in amplitude and phase, while the bilinear solution leads just to errors in phase (Figure ).

 

unit Figure 1 Explicit, bilinear and implicit approximations plotted on the unit circle. The solid line corresponds to the exact exponential solution.

 

exap Figure 2 Amplitude and phase errors for the explicit, bilinear and implicit approximations.

If, for notation simplicity, we define  

 153#153 (63)

the WEMVA equation () can be written as  

 154#154 (64)

and so the linearizations corresponding to the explicit, bilinear and implicit solutions respectively become Aparently, just the first equation in () provides a linear relationship between 146#146 and 141#141. However, a simple re-arrangement of terms leads to

For MVA, both the background (157#157) and perturbation wavefields (146#146) are known, so it is not a problem to incorporate them in the linear operator. In any of the cases described in Equation (), the approximations can be symbolically written using the fitting goal

158#158 (67)

where the data 9#9 is the wavefield perturbation, and the model 10#10 is the slowness perturbation. The same operator 11#11 is used for inversion in all situations, the only change being in the wavefield that is fed into the linear operator. Therefore, the new operators are not more expensive than the Born operator.

All linear relationships in Equation () belong to a family of approximations of the general form

   159#159 (68)

The various approximations can be obtained using appropriate values for the parameter 160#160.All forms of Equation (), however, are approximations to the exact non-linear relation (), therefore they are all likely to break for large values of the phase, or equivalently large values of the slowness perturbation or frequency. Nevertheless, these approximations enable us to achieve higher accuray in slowness estimation as compared to the simple Born approximation.

An interesting comparison can be made between the extreme members of the sequence given by Equation (): for 161#161 we use the background wavefield 157#157, and for 162#162 we use the full wavefield 163#163. The physics of scattering would recommend that we use the later form, since the scattered wavefield (146#146) is generated by the total wavefield (164#164), and not by an approximation of it (157#157), thus naturally accounting for multiple scattering effects. The later situation also corresponds to what is known in the scattering literature as wavefield renormalization (). The details of these ideas and their implications remain open for future research.

Finally, we note that Equation () cannot be used for forward modeling of the wavefield perturbations 146#146, except for the particular case 161#161, since the output quantity is contained in the operator itself. However, we can use this equation for inversion for any choice of the parameter 165#165.