![]() |
![]() |
![]() |
![]() | A preconditioning scheme for full waveform inversion | ![]() |
![]() |
To make a compelling case, we kept only four shots, 2.5 Km. apart. First, we show in Figure 3 a comparison between the gradient without preconditioning
and the gradient with preconditioning back in the velocity space
.
Because only 4 shots are present, the unpreconditioned gradient looks noisy and resemble geology in very few locations only. Some authors suggest attenuating the high wavenumbers in the gradient by smoothing it after each iteration (Ravaut et al., 2004), where the size of the smoothing operator in the horizontal and vertical directions is a function of an average wavelength at a given frequency. This bears a resemblance with our proposed scheme but
doesn't allow for directional smoothing.
On the contrary, thanks to the preconditioning with directional Laplacian filters, the reprojected gradient in Figure 3b shows the geology captured in the dip field of Figure 1b very well.
![]() |
---|
MarmGRADSPARSESEG
Figure 3. (a) Gradient ![]() ![]() |
![]() ![]() |
![]() |
---|
MarmSPARSESEG
Figure 4. (a) Inversion result for the unpreconditioned inversion. (b) Inversion result for the preconditioned inversion with directional Laplacian filters. The geology at the reservoir level is recovered very well in (b). |
![]() ![]() |
Now, we show in Figure 4 the inversion results for the unpreconditioned (Figure 4a) and preconditioned inversion (Figure 4b). Although quite noisy, the unpreconditioned result shows the geology very well: the presence of low frequencies in the data, along with the multi-scale approach, act as regularization operators. This effect will be less pronounced with real data where low frequencies are often missing. The preconditioned inversion result in Figure 4b is much cleaner: the geology at the reservoir level is recovered very well.
![]() |
![]() |
![]() |
![]() | A preconditioning scheme for full waveform inversion | ![]() |
![]() |