Sparse reflectivity inversion for nonstationary seismic data

Abstract

Conventional reflectivity inversion methods are based on a stationary convolution model and theoretically require stationary seismic traces as input (i.e., those free of attenuation and dispersion effects). Reflectivity inversion for nonstationary data, which is typical for field surveys, requires us to first compensate for the earth’s [Formula: see text]-filtering effects by inverse [Formula: see text] filtering. However, the attenuation compensation for inverse [Formula: see text] filtering is inherently unstable, and offers no perfect solution. Thus, we presented a sparse reflectivity inversion method for nonstationary seismic data. We referred to this method as nonstationary sparse reflectivity inversion (NSRI); it makes the novel contribution of avoiding intrinsic instability associated with inverse [Formula: see text] filtering by integrating the earth’s [Formula: see text]-filtering operator into the stationary convolution model. NSRI also avoids time-variant wavelets that are typically required in time-variant deconvolution. Although NSRI is initially designed for nonstationary signals, it is suitable for stationary signals (i.e., using an infinite [Formula: see text]). The equations for NSRI only use reliable frequencies within the seismic bandwidth, and the basis pursuit optimizes a cost function of mixed [Formula: see text] norms to derive a stable and sparse solution. Synthetic examples show that NSRI can directly retrieve reflectivity from nonstationary data without advance inverse [Formula: see text] filtering. NSRI is satisfactorily stable in the presence of severe noise, and a slight error in the [Formula: see text] value does not greatly disturb the sensitivity of NSRI. A field data example confirmed the effectiveness of NSRI.