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

The method of modeling wavefield propagation with an implicit finite-difference approximation to the two-way acoustic isotropic wave equation, using spectral factorization and helical deconvolution, exhibits instability of the propagating wavefield as the time step is increased. In this study, I test several potential sources of the instability problem: the implicit finite-difference scheme itself, the precision of the floating point representation of the filter coefficients, the number of filter coefficients, and the spectral factorization method. None of these issues is the cause for the apparent instability.

- Introduction
- Review of methodology
- Implementation of methodology with increasing time step size

- Effect of floating point precision
- Effect of number of spectrally factorized coefficients
- Wilson-Burg vs. Kolmogoroff spectral factorization
- conclusion and future work
- Bibliography
- APPENDIX A
- About this document ...

2010-11-26