Seismic imaging of simultaneous-source data using constrained least-squares reverse time migration

Simultaneous-source acquisition can help obtain better-sampled seismic data with a tremendously faster efficiency. However, the interference caused between simultaneous sources will cause strong artifacts in the migrated image. We propose to directly migrate the blended data without deblending. We treat the imaging for blended seismic data as a least-squares inversion problem, and use preconditioned conjugate gradient algorithm to iteratively solve it. The constraining operator is chosen as the structural smoothing operator, which helps to smooth the seismic image along the local structure. Compared with other available approaches, the proposed approach can obtain a true-amplitude and much cleaner image.