Multi-trace basis pursuit inversion with spatial regularization

Basis pursuit inversion (BPI) was developed to derive spiky reflection coefficients, which was applied on post-stack seismic data trace-by-trace for subsurface reflectivity. Such a trace-by-trace operation has no control on lateral continuity. To address this, we extend the former BPI to a multi-trace procedure with spatial regularization - a spatial derivative operator - to stabilize the lateral variation. The original objective function of BPI is reformed with a 'Z' shape spatial derivative imposed, which can balance the lateral continuity and vertical resolution. Synthetic tests show improved lateral continuity by spatial regularized BPI in the case of noise contamination. Real data applications demonstrate its effectiveness to enhance the lateral continuity of inverted spiky reflectivity.