Complex seismic wavefi eld interpolation based on the Bregman iteration method in the sparse transform domain

In seismic prospecting, fi eld conditions and other factors hamper the recording of the complete seismic wavefi eld; thus, data interpolation is critical in seismic data processing. Especially, in complex conditions, prestack missing data affect the subsequent highprecision data processing workfl ow. Compressive sensing is an effective strategy for seismic data interpolation by optimally representing the complex seismic wavefi eld and using fast and accurate iterative algorithms. The seislet transform is a sparse multiscale transform well suited for representing the seismic wavefield, as it can effectively compress seismic events. Furthermore, the Bregman iterative algorithm is an efficient algorithm for sparse representation in compressive sensing. Seismic data interpolation methods can be developed by combining seismic dynamic prediction, image transform, and compressive sensing. In this study, we link seismic data interpolation and constrained optimization. We selected the OC-seislet sparse transform to represent complex wavefields and used the Bregman iteration method to solve the hybrid norm inverse problem under the compressed sensing framework. In addition, we used an H-curve method to choose the threshold parameter in the Bregman iteration method. Thus, we achieved fast and accurate reconstruction of the seismic wavefi eld. Model and fi eld data tests demonstrate that the Bregman iteration method based on the H-curve norm in the sparse transform domain can effectively reconstruct missing complex wavefi eld data.

Study of denoising method for nonhyperbolic prestack seismic reflection data

Removing random noise in seismic data is a key step in seismic data processing. A failed denoising may introduce many artifacts, and lead to the failure of final processing results. Seislet transform is a wavelet-like transform that analyzes seismic data following variable slopes of seismic events. The local slope is the key of seismic data. An earlier work used traditional normal moveout(NMO) equation to construct velocity-dependent(VD) seislet transform, which only adapt to hyperbolic condition. In this work, we use shifted hyperbola NMO equation to obtain more accurate slopes in nonhyperbolic situation. Self-adaptive threshold method was used to remove random noise while preserving useful signal. The synthetic and field data tests demonstrate that this method is more suitable for noise attenuation.