Multi-trace semi-blind nonstationary deconvolution

We proposed a multitrace semi-blind nonstationary deconvolution method. The proposed method estimates reflectivity and source wavelet simultaneously for pursuing high resolution seismic processing. The mathematical framework is derived based on convolution exchange law and Fourier transform property. In this framework, seismic records are treated as the convolution of a time-varying wavelet and non-attenuated reflectivity or the convolution of a constant wavelet and attenuated reflectivity. Using these two equivalence relations, we devise an objective function containing two variables, the reflectivity and wavelet. In addition, we add the 2D total variation constraint to the cost function, which preserves lateral and vertical continuity of the estimated reflectivity. The cost function is solved by alternating iteration and proximal splitting methods, under the assumptions of a known attenuation model and sparse reflectivity. Besides, the mathematical framework is extended to implement semi-blind deconvolution in an approximate layered earth model. To demonstrate the effectiveness of the proposed method, we apply the proposed method to synthetic data and field data, and confirm that the proposed method can achieve better reflectivity and source wavelet.