Prestack imaging of seismic data using Lp iterative reweighted least-squares wavefield extrapolation filters in the frequency-space domain

A stable explicit depth wavefield extrapolation is obtained using [Formula: see text] iterative reweighted least-squares (IRLS) frequency-space ([Formula: see text]-[Formula: see text]) finite-impulse response digital filters. The problem of designing such filters to obtain stable images of challenging seismic data is formulated as an [Formula: see text] IRLS minimization. Prestack depth imaging of the challenging Marmousi model data set was then performed using the explicit depth wavefield extrapolation with the proposed [Formula: see text] IRLS-based algorithm. Considering the extrapolation filter design accuracy, the [Formula: see text] IRLS minimization method resulted in an image with higher quality when compared with the weighted least-squares method. The method can, therefore, be used to design high-accuracy extrapolation filters.

Sparse signal recovery with prior information by iterative reweighted least squares algorithm

In this paper, the iterative reweighted least squares (IRLS) algorithm for sparse signal recovery with partially known support is studied. We establish a theoretical analysis of the IRLS algorithm by incorporating some known part of support information as a prior, and obtain the error estimate and convergence result of this algorithm. Our results show that the error bound depends on the best {(s+k)}-term approximation and the regularization parameter λ, and convergence result depends only on the regularization parameter λ. Finally, a series of numerical experiments are carried out to demonstrate the effectiveness of the algorithm for sparse signal recovery with partially known support, which shows that an appropriate q ({0<q<1}) can lead to a better recovery performance than that of the case {q=1}.