Abstract
The pre-stack depth migration of reflection seismic data can be expressed, in the framework of waveform inversion, as a linear least squares problem. Together with the precise definition of this operator, we detail additional main characteristics of the forward model, like its huge size, its sparsity and the composition with convolution. It ends up with a so-called discrete ill-posed problem, whose acceptable solutions have to undergo a regularization procedure. Both direct and iterative methods have been implemented with specific attention to the convolution, and then applied to a given data set: a synthetic 2-dimensional profile of revealing size with some added noise. The efficiency with regard to computational effort and storage requirements is evaluated. The needed regularization of the solution is thoroughly studied in both cases. From the point of the global inverse problem, the extra feature of providing a solution that can be differentiated with respect to a parameter such as background velocity is also discussed.
Similar content being viewed by others
References
P.R. Amestoy, I.S. Duff and C. Puglisi, Multifrontal QR factorization in a multiprocessor environment, Numer. Linear Algebra Appl. 3(4) (1996) 275–300.
A. Björck, Numerical Methods for Least Squares Problems (SIAM, Philadelphia, PA, 1996).
N. Bleistein, Mathematical Methods for Wave Phenomena (Academic Press, New York, 1984).
G. Böhm and A.L. Vesnaver, In quest of the grid, Geophysics 64(4) (1999) 1116–1125.
D. Calvetti, L. Reichel and Q. Zhang, Conjugate gradient algorithms for symmetric inconsistent linear systems, in: Cornelius Lanczos Internat. Centenary Conference, eds. J.D. Brown, M.T. Chu, D.C. El-lison and R.J. Plemons (SIAM, Philadelphia, PA, 1994) pp. 267–272.
T.F. Chan, Rank revealing QR factorizations, Linear Algebra Appl. 88/89 (1987) 67–82.
G. Chavent, Duality method for waveform inversion, Technical Report 2975, INRIA, Le Chesnay, France (1996).
G. Chavent, F. Clément and S. Gómez, Waveform inversion by migration based travel time formula-tion, in: 3rd Internat. Conf. on Mathematical and Numerical Aspects of Wave Propagation (SIAM, Philadelphia, PA, 1995) pp. 1179–1182.
J.F. Claerbout, Towards a unified theory of reflector mapping, Geophysics 36 (1971) 467–481.
Y.-H. De Roeck, Sparse linear algebra and geophysical migration, Technical Report 3876, INRIA, Le Chesnay, France (2000).
I.S. Duff, A.M. Erisman and J.K. Reid, Direct Methods for Sparse Matrices, Monographs on Numer-ical Analysis (Clarendon Press, Oxford, 1986).
R.D. Fierro and P.C. Hansen, Accuracy of TSVD solutions computed from rank-revealing decompositions, Numer. Math. 70 (1995) 453–471.
G. Golub and C. Van Loan, Matrix Computations(North Oxford Academic, 1983).
M. Hanke, Conjugate Gradient Type Methods for Ill-Posed Problems, Pitman Research Notes in Mathematics (Longman Scientific and Technical, Harlow, 1995).
M. Hanke and C.R. Vogel, Two-level preconditionioners for regularized inverse problems I: Theory, Numer. Math. 83 (1999) 385–402.
P.C. Hansen, Regularization, GSVD and truncated GSVD, BIT 29 (1989) 491–504.
P.C. Hansen, Truncated singular value decomposition solutions to discrete ill-posed problems with ill-determined numerical rank, SIAM J. Sci. Statist. Comput. 11(3) (1990) 503–518.
P.C. Hansen, Analysis of discrete ill-posed problems by means of the l-curve, SIAM Rev. 34(4) (1992) 561–580.
P.C. Hansen, Rank-Deficient and Discrete Ill-Posed Problems, Monographs on Mathematical Modeling and Computation (SIAM, Philadelphia, PA, 1998).
M.T. Heath, Some extensions of an algorithm for sparse linear least squares problems, SIAM J. Sci. Statist. Comput. 3(2) (1982) 223–237.
C.C. Paige and M.A. Saunders, LSQR: An algorithm for sparse linear equations and sparse least squares, ACM Trans. Math. Software 8(1) (1982) 43–71.
R.-E. Plessix, G. Chavent and Y.-H. De Roeck, A quantitative Kirchhoff migration to esimate the 2D velocity distribution, in: 3rd Internat. Conf. on Mathematical and Numerical Aspects of Wave Propagation (SIAM, Philadelphia, PA, 1995) pp. 704–712.
R.-E. Plessix, Y.-H. De Roeck and G. Chavent, Waveform inversion of reflection seismic data for kinematic parameters by local optimization, SIAM J. Sci. Comput. 20(3) (1999) 1033–1052.
G. Quintana-Ortí, X. Sun and C.H. Bischof, A BLAS-3 version of the QR factorisation with column pivoting, SIAM J. Sci. Comput. 19(5) (1998) 1486–1494.
G.W. Stewart, Four algorithms for the efficient computation of truncated pivoted QR approximations to a sparse matrix, Numer. Math. (1999).
W.W. Symes, The reflection inversion problem for acoustic waves, in: Mathematical and Numerical Aspect of the Wave Propagation Phenomena (SIAM, Philadelphia, PA, 1991) pp. 423–433.
W.W. Symes and J.J. Carazzone, Velocity inversion by differential semblance optimization, Geophysics 56(5) (1991) 654–663.
A. Tarantola, Linearized inversion of seismic reflection data, Geophysical Prospecting 32 (1984) 998–1015.
G. Wahba, Pratical approximate solutions to linear operator equations when the data are noisy, J. Nu-mer.Anal. 14 (1977) 651–667.
Ö. Yilmaz, Seismic Data Processing, Investigations in Geophysics, No. 2 (Society of ExplorationGeophysicists, Tulsa, 1987).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
De Roeck, YH. Sparse Linear Algebra and Geophysical Migration: A Review of Direct and Iterative Methods. Numerical Algorithms 29, 283–322 (2002). https://doi.org/10.1023/A:1015210828153
Issue Date:
DOI: https://doi.org/10.1023/A:1015210828153