Approximate Finite Rate of Innovation Based Seismic Reflectivity Estimation

Reflectivity inversion is an important deconvolution problem in reflection seismology that helps to describe the subsurface structure. Generally, deconvolution techniques iteratively work on the seismic data for estimating reflectivity. Therefore, these techniques are computationally expensive and may be slow to converge. In this paper, a novel method for estimating reflectivity signals in seismic data using an approximate finite rate of innovation (FRI) framework, is proposed. The seismic data is modeled as a convolution between the Ricker wavelet and the FRI signal, a Dirac impulse train. Relaxing the accurate exponential reproduction limitation given by generalised Strang-Fix (GSF) conditions, we develop a suitable sampling kernel utilizing Ricker wavelet which allows us to estimate the reflectivity signal. The experimental results demonstrate that the proposed approximate FRI framework provides a better reflectivity estimation than the deconvolution technique for medium-to-high signal-to-noise ratio (SNR) regimes with nearly 18% of seismic data.