Block row recursive least squares migration
This work addresses the computational bottleneck of least squares migration in geophysics, offering a recursive approach for real-time or large-scale seismic imaging.
The authors propose a block row recursive least squares method for wave equation least squares migration, enabling efficient processing with sliding windows. Numerical experiments demonstrate the technique's effectiveness for geophysical imaging.
Recursive estimates of large systems of equations in the context of least squares fitting is a common practice in different fields of study. For example, recursive adaptive filtering is extensively used in signal processing and control applications. The necessity of solving least squares problem via recursive algorithms comes from the need of fast real-time signal processing strategies. Computational cost of using least squares algorithm could also limits the applicability of this technique in geophysical problems. In this paper, we consider recursive least squares solution for wave equation least squares migration with sliding windows involving several rank K downdating and updating computations. This technique can be applied for dynamic and stationary processes. One can show that in the case of stationary processes, the spectrum of the preconditioned system is clustered around one and the method will converge superlinearly with probability one, if we use enough data in each windowed setup. Numerical experiments are reported in order to illustrate the effectiveness of the technique for least squares migration.