Least squares migration can eliminate the artifacts introduced by the direct imaging of irregular seismic data but is computationally costly and of slow convergence. In order to suppress the migration noise, we propos...Least squares migration can eliminate the artifacts introduced by the direct imaging of irregular seismic data but is computationally costly and of slow convergence. In order to suppress the migration noise, we propose the preconditioned prestack plane-wave least squares reverse time migration （PLSRTM） method with singular spectrum constraint. Singular spectrum analysis （SSA） is used in the preconditioning of the take-off angle-domain common-image gathers （TADCIGs）. In addition, we adopt randomized singular value decomposition （RSVD） to calculate the singular values. RSVD reduces the computational cost of SSA by replacing the singular value decomposition （SVD） of one large matrix with the SVD of two small matrices. We incorporate a regularization term into the preconditioned PLSRTM method that penalizes misfits between the migration images from the plane waves with adjacent angles to reduce the migration noise because the stacking of the migration results cannot effectively suppress the migration noise when the migration velocity contains errors. The regularization imposes smoothness constraints on the TADCIGs that favor differential semblance optimization constraints. Numerical analysis of synthetic data using the Marmousi model suggests that the proposed method can efficiently suppress the artifacts introduced by plane-wave gathers or irregular seismic data and improve the imaging quality of PLSRTM. Furthermore, it produces better images with less noise and more continuous structures even for inaccurate migration velocities.展开更多
基金supported by the National Science and Technology Major Project(No.2016ZX05014-001-008)the National Key Basic Research Program of China(No.2014CB239006)+2 种基金the National Natural Science Foundation of China(Nos.41104069 and 41274124)the Open foundation of SINOPEC Key Laboratory of Geophysics(No.33550006-15-FW2099-0033)the Fundamental Research Funds for Central Universities(No.16CX06046A)
文摘Least squares migration can eliminate the artifacts introduced by the direct imaging of irregular seismic data but is computationally costly and of slow convergence. In order to suppress the migration noise, we propose the preconditioned prestack plane-wave least squares reverse time migration （PLSRTM） method with singular spectrum constraint. Singular spectrum analysis （SSA） is used in the preconditioning of the take-off angle-domain common-image gathers （TADCIGs）. In addition, we adopt randomized singular value decomposition （RSVD） to calculate the singular values. RSVD reduces the computational cost of SSA by replacing the singular value decomposition （SVD） of one large matrix with the SVD of two small matrices. We incorporate a regularization term into the preconditioned PLSRTM method that penalizes misfits between the migration images from the plane waves with adjacent angles to reduce the migration noise because the stacking of the migration results cannot effectively suppress the migration noise when the migration velocity contains errors. The regularization imposes smoothness constraints on the TADCIGs that favor differential semblance optimization constraints. Numerical analysis of synthetic data using the Marmousi model suggests that the proposed method can efficiently suppress the artifacts introduced by plane-wave gathers or irregular seismic data and improve the imaging quality of PLSRTM. Furthermore, it produces better images with less noise and more continuous structures even for inaccurate migration velocities.
