Implicit finite difference in time-space domain with the helix transform |

- forward deconvolution with the spectrally factorized coefficients of the above finite-difference weights.
- adjoint deconvolution.
- adjoint convolution.
- forward convolution.

The sequence and it's results are shown in Figure 5. The final result is not a perfect impulse, however the remaining values in the wavefield are 6 or more orders of magnitude smaller than the central impulse.

imp-test3-in
Input impulse used for helical convolution test
Figure 4. | |
---|---|

imp-test3-panels
Results of convolving and deconvolving with filter coefficients whose convolution yields the finite-difference weights required by the implicit finite-difference approximation in Equation 17
Figure 5. | |
---|---|

exp-vs-imp-1d
1D Explicit Vs. Implicit finite difference with constant velocity =
. Source is a Ricker wavelet with central frequency =
. The explicit time step was
. The implicit time step was
.
.
Figure 6. |
---|

In order to test whether the derived implicit finite-difference coefficients are actually valid for wavefield propagation, the next step was to compare standard explicit finite difference to implicit propagation with these coefficients in one dimension. The implicit solution was done using the SEPlib module `rtris` which recursively solves a tridiagonal equation system.
The results are in Figure 6. Severe dispersion is visible in the implicitly propagated wavefield. So far I have been unable to remove this dispersion without damaging the kinematics of the wavefield.

On the other hand, the wavefield derived by implicit finite difference does not diverge even when increasing the temporal time step beyond the stability limit of the explicit solver.

imp-vs-helimp-snap
2D Implicit Vs. Helical Implicit finite difference with constant velocity =
. Source is a Ricker wavelet with central frequency =
. The time step for the top figures was
, and for the bottom figures
.
.
Figure 7. |
---|

exp-vs-helimp-snap
2D Explicit Vs. Helical Implicit finite difference with constant velocity =
. Source is a Ricker wavelet with central frequency =
. The time step was
- at the stability limit of the explicit scheme for the 2nd order in time and space approximation.
.
Figure 8. |
---|

Implicit finite difference in time-space domain with the helix transform |

2010-05-19