Implementing implicit finite-difference in the time-space domain using spectral factorization and helical deconvolution |

An example of the lack of the effect of number of coefficients on the propagation is shown in Figure 9. This 1D example shows how propagation using 2 spectrally factorized filter coefficients is basically identical to propagation when using 50 filter coefficients. Another indication comes from observing the filter coefficients themselves. This example was produced by a 2nd order scheme, which means that there are only 2 finite-difference weights. When factorizing using only 2 filter coefficients, the Wilson-Burg algorithm (for the propagation parameters used in Figure 9) yielded:

The coefficients are displayed in order of lags, so the
is the zero-lag filter coefficient. Correlating these coefficients, we get:

at lag 0 and lag 1, which are equal to the floating point representations of the finite-difference weights for Figure 9.

Factorizing using 50 filter coefficients produced (only the first four coefficients are shown, the rest are in A-5 and A-6):

After lag 20, the coefficients are all zeros. Note that the first two coefficients are identical to the ones produced by the factorizer when requesting only two coefficients.The correlation of this filter is:

This correlation again shows the accurate representation of the finite-difference coefficients at lag 0 and lag 1. In addition, the correlation produces a set of other values at later lags, which are much smaller than the weights themselves.

The fact that two filter coefficients were sufficient to produce the finite-difference weights by correlation was interesting, but what is more important is to test what effect the change in the number of coefficients might have on the deconvolution process. Correlating the coefficients is like convolving them over a spike, once in the forward direction and once in reverse order. In order to test the exact effect that a change in the number of coefficients had on the deconvolution, I tested the result of deconvolving the coefficients over a spike. Here as well, the result was identical. I shall spare displaying the numbers themselves for this case.

In summary, I could not a find a combination of parameters (of propagation or of factorization) for which wavefield propagation was more stable if more filter coefficients were used.

helimp-ncoeffs
1D helical implicit finite-difference with 2 spectrally factorized filter coefficients (top), and 50 coefficients (bottom). Velocity =
,
,
,
.[Figure 9. ER]
| |
---|---|

Implementing implicit finite-difference in the time-space domain using spectral factorization and helical deconvolution |

2010-11-26