High resolution Kirchhoff seismic migration via 1-norm regularized least-squares


1 Ph.D. Student, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran

2 Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran

3 Associate Professor, Department of Earth Physics, Institute of Geophysics, University of Tehran, Iran


For decades Kirchhoff migration has been one of the simplest migration algorithms and also the most frequently used method of migration in industry. This is due to its relatively low computational cost and its flexibility in handling acquisition and topography irregularities. The standard seismic migration operator can be regarded as the adjoint of a seismic forward modeling operator, which acts on a set of subsurface parameters to generate the observed data. Such adjoint operators are able to provide an approximate inverse of the forward modeling operator and only recover the time of the events (Claerbout, 1992). They cannot retrieve the amplitude of reflections, thus leading to a decrease in the resolution of the final migrated image. The standard seismic migration (adjoint) operators can be modified to better approximate the inverse operators. Least-squares migration (LSM) techniques have been developed to fully inverse the forward modeling procedures by minimizing the difference between observed and modeled data in a least-squares sense. An LSM is able to reduce the (Kirchhoff) migration artifacts, enhance the resolution and retrieve seismic amplitudes. Although implementing LSM instead of conventional migration, leads to resolution enhancement. It also brings some new numerical and computational challenges which need to be addressed properly. Due to the ill-conditioned nature of the inverse operator and also incompleteness of the data, the method generates unavoidable artifacts which severely degrade the resolution of the migrated image obtained by the non-regularized LSM method. The instability of LSM methods suggests developing a regularized algorithm capable of including reasonable physical constraints. Including the seismic wavelet into the migration operator, migration will generate the earth reflectivity image which can be considered as a sparse image, so applying the sparseness constraint, e.g., via the minimization of the 1-norm of reflectivity model, can help to regularize the model and prevent it from getting noisy artifacts (Gholami and Sacchi, 2013).
In this article, based on the Bregmanized operator splitting (BOS), we propose a high resolution migration algorithm by applying sparseness constraints to the solution of least-squares Kirchhoff migration (LSKM). The Bregmanized operator splitting is employed as a solver of the generated sparsity-promoting LSKM for its simplicity, efficiency, stability and fast convergence. Independence of matrix inversion and fast convergence rate are two main properties of the proposed algorithm. Numerical results from field and synthetic seismic data show that migrated sections generated by this 1-norm regularized Kirchhoff migration method are more focused than those generated by the conventional Kirchhoff/LS migration.
Regular spatial sampling of the data at Nyquist rate is another major challenge which may not be achieved in practice due to the coarse source-receiver distributions and presence of possible gaps in the recording lines. The proposed model-based migration algorithm is able to handle the incompleteness issues and is stable in the presence of noise in the data. In this article, we tested the performance of our proposed method on synthetic data in the presence of coarse sampling and also acquisition gaps. The results confirmed that the proposed sparsity-promoting migration is able to generate accurate migrated images from incomplete and inaccurate data.


Main Subjects

Beck‎, ‎A. ‎and ‎Teboulle‎, M‎., 2009, A fast iterative sherinkage-thresholding algorithm for linear inverse problems, SIAM Journal on Imaging Sciences,  2(1)‎, ‎183-202‎.
Claerbout, J. F., 1992, Earth soundings analysis: processing versus inversion, Blackwell scientific publications.‎
Combettes, ‎‎‎P. L. and Wajs, V. R., 2005, Signal recovery by proximal forward-backward splitting, Multiscale Modeling and Simulation, 4(4), 1168-1200.
Dai, W., 2012, Multisource least-squares reverse time migration: Ph.D. thesis, King Abdullah University of Science and Technology.
Daubechies‎, ‎I., Defriese‎, M. and ‎De Mol‎, C‎., 2004‎, An iterative thresholding algorithm for linear inverse problems with a sparsity constraint, ‎Communications on Pure and Applied Mathematics‎‎, 57(11), ‎1413–1457‎.
‎Daubechies‎, ‎I.‎, DeVore‎, R‎., Fornasier‎, ‎M‎. and‎ ‎Gunturk‎, S‎., ‎2010‎, Iteratively re-weighted least squares minimization for sparse recovery, ‎Communications on Pure and Applied Mathematics‎, 63(1), ‎1–38‎.
Figueiredo‎, ‎M‎. ‎A‎. ‎T.‎, Nowak‎, R‎. ‎D‎.‎ ‎and Wright‎, S‎. ‎J‎., 2007, Gradient projection for sparse reconstruction‎: ‎Application to compressed sensing and other inverse problems, IEEE Journal of selected topics In Signal Processing, 1(4), 586-597.
Gholami, A. and Sacchi, M. D., 2013, Fast 3D blind seismic deconvolution via constrained total variation and GCV, SIAM Journal on Imaging Sciences, 6(4), 2350-2369.
Gholami A., 2017, Deconvolutive Radon transform, Geophysics, 82(2), V117-V125.
‎Goldstein‎, ‎T.‎, ‎and Osher‎, S‎., 2009, The split bregman method for l1 regularized problems‎, ‎SIAM Journal on ‎Imaging Science‎, 2(2)‎, ‎323-343‎.
Gray, S. H., Etgen, J., Dellinger, J. and Whitmore, D., 2001, Seismic migration problems and solutions, Geophysics, 66(5), 1622-1640.
Ji, J. 1997, Least squares imaging, datuming, and interpolation using the wave equation, Stanford Exploration Project, 75, 121-135.
Kühl, H. and Sacchi, M. D., 2003, Least-squares wave equation migration for AVP/AVA inversion, Geophysics, 68, 262–273.
Keys, R. and Foster, D., 1998, Comparison of seismic inversion methods on a single real data set: Society of Exploration Geophysics.
Lari, H. H. and Gholami, A., 2014, Curvelet-TV regularized Bregman iteration for seismic random noise attenuation, Journal of Applied Geophysics, 109(1), 233-241.
Loewenthal, D., Lu, L., Roberson, R. and Sherwood, J.W.C., 1976, The wave equation applied to migration, Geophysics Prospecting, 24, 380-399.
Luo, S. and Hale, D., 2014, Least-squares migration in the presence of velocity errors, Geophysics, 79(4), 153-161.
Miller, D., Oristaglio, M. and Beylkin, G., 1987, A new slant on seismic imaging: migration and integral geometry, Geophysics, 52, 943-964.
Nemeth, T., Wu, C. and Schuster, G. T., 1999, Least-squares migration of incomplete reflection data, Geophysics, 64, 208-221.
Osher, S., Burger, M., Goldfarb, D., Xu, J. and Yin, W., 2005, An iterative regularization method for total variation-based image restoration, Multiscale Modeling and Simulation, 4, 460-489.‎‏‎
Østmo, S., Plessix, R.E., 2002, Finite-difference iterative migration by linearized waveform inversion in the frequency domain, 72nd Annual International Meeting, Society of Exploration Geophysics, Expanded Abstracts, 1384–1387.
Plessix, R. E. and Mulder, W. A., 2004, Frequency-domain finite-difference amplitude-preserving migration, Geophysical Journal International, 157, 975–987.
Schuster, G. T., 2017, Seismic Inversion, SEG.
Tarantola, A., 1984, Inversion of seismic reflection data in the acoustic approximation, Geophysics, 49, 1259-1266.
Yilmaz, O., 2001, Seismic data analysis: processing, inversion, and interpretation of seismic data, Society of Exploration Geophysics.
Yousefzadeh, A., 2012, ‎High resolution seismic imaging using least squares migration, Ph.D. thesis, Calgary, Alberta.‎
Zand, T. and Gholami, A., 2018, Total-variation based velocity inversion with Bregmanized operator splitting algorithm, Journal of Applied Geophysics, 151, 1-10.
Zhang, X.‎, ‎Burger, ‎M., Bresson, X. and Osher, S., 2010, Bregmanized nonlocal regularization for deconvolution and sparse reconstruction, SIAM Journal on Imaging Sciences, 3(3), 253-276.