WEMVA is closer to conventional MVA than other wave-equation tomography methods Bunks et al. (1995); Forgues et al. (1998); Noble et al. (1991) because it tries to maximize the quality of the migrated image instead of trying to match the recorded data. In this respect, our method is related to Differential Semblance Optimization Symes and Carazzone (1991) and Multiple Migration Fitting Chavent and Jacewitz (1995). However, with respect to these two methods, our method has the advantage of exploiting the power of residual prestack migration to speed up the convergence, and it also gives us the ability to guide the inversion by geologic interpretation.
WEMVA benefits from the same advantages over traveltime estimation methods as wave-equation migration benefits over Kirchhoff migration. The most important among them are the accurate handling of complex wavefields which are characterized by multipathing, and the band-limited nature of the imaging process, which can handle sharp velocity variations much better than traveltime-based methods Woodward (1992). Complex geology, therefore, is where WEMVA is expected to provide the largest benefits.
WEMVA is based on a linearization of the downward-continuation operator using the Born approximation. This approximation leads to severe limitations on the magnitude and size of the anomalies that can be handled. It, therefore, cannot operate successfully exactly in the regions of highest complexity. Other methods of linearization are possible Sava and Fomel (2002), but neither one allows for arbitrarily large anomalies.
In our early tests Biondi and Sava (1999), we construct the image perturbation using Prestack Stolt Residual Migration (PSRM) Sava (2000); Stolt (1996). In summary, this residual migration method provides updated images for new velocity maps that correspond to a fixed ratio () of the new velocity with respect to the original velocity map. Residual migration is run for various ratio parameters, and finally we pick the best image by selecting the flattest gathers at every point.
The main disadvantage of building the image perturbation using PSRM is that, for large velocity ratio parameters (), the background and improved images can get more than out of phase. Therefore, the image perturbation computed by the forward operator and the one computed by residual migration are fundamentally different, and can have contradictory behaviors when using the Born-linearized WEMVA operator for inversion.
Alternative methods can be used to create image perturbations for WEMVA, in compliance with the Born approximation and computed directly from the background image Sava and Biondi (2001). We use an analytic differential procedure starting from the background image and leading to image perturbations similar to the ones created using the forward WEMVA operator.