Comparison of scaling methods for waveform inversion

Abstract

Waveform inversion can lead to faint images for later times due to geometrical spreading. The proper scaling of the steepest-descent direction can enhance faint images in waveform inversion results. We compare the effects of different scaling techniques in waveform inversion algorithms using the steepest-descent method. For the scaling method we use the diagonal of the pseudo-Hessian matrix, which can be applied in two different ways. One is to scale the steepest-descent direction at each frequency independently. The other is to scale the steepest-descent direction summed over the entire frequency band. The first method equalizes the steepest-descent directions at different frequencies and minimizes the effects of the band-limited source spectrum in waveform inversion. In the second method, since the steepest-descent direction summed over the entire frequency band is divided by the diagonal of the pseudo-Hessian matrix summed over the entire frequency band, the band-limited property of the source wavelet spectrum still remains in the scaled steepest-descent directions. The two scaling methods were applied to both standard and logarithmic waveform inversion. For standard waveform inversion, the method that scales the steepest-descent direction at every frequency step gives better results than the second method. On the other hand, logarithmic waveform inversion is not sensitive to the scaling method, because taking the logarithm of wavefields automatically means that results for the steepest-descent direction at each frequency are commensurate with each other. If once the steepest-descent directions are equalized by taking the logarithm of wavefields in logarithmic waveform inversion, the additional equalizing effects by the scaling method are not as great as in conventional waveform inversion.