1887
Volume 28, Issue 4
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

Tomographic estimation of near-surface velocity anomalies has direct application in prestack depth migration. Seismic velocity information extracted from traveltime measurements often suggests unrealistic geological settings. The linear system constructed by ray tracing is typically very large, sparse and poorly conditioned. In most cases it is too large to invert using a generalised inversion algorithm. Instead, a least squares method to solve the system is employed.

Limited angular ray coverage traveltime measures using three synthetic reflection data sets have been shown to allow reliable velocity image reconstruction of anomalous zones that vary up to 50 per cent from the background.

Tomographic imaging using four summation expansion methods, namely, the algebraic reconstruction techniques (ART), and the simultaneous iterative reconstruction technique (SIRT), a conjugate gradient (CG) method and a singular value decomposition (SVD) method, were used to image the synthetic data sets, and the result of each of these tomographic algorithms was evaluated based on their rms misfit and residual, solution stability and sensitively to data noise and algorithmic operation time.

Modifications to the algorithms that stabilise the inversion are shown to improve the velocity images obtained and enhance solution convergence rate. These modifications include the application of damping parameters, weighting schemes, smoothing filters including average and median, implemented both pre- and post-solution as well as convolutional quelling. The effect of noisy traveltime data on tomographic images, and useful methods for dealing with this, are also studied.

Loading

Article metrics loading...

/content/journals/10.1071/EG997369
1997-09-01
2026-01-20
Loading full text...

Full text loading...

References

  1. Campbell, S., and Meyer, C, 1979, Generalized inverses of linear transformations: Pitman Publ. Ltd.
  2. Censor, Y., 1981, Row-action methods for huge and sparse systems and their applications: SIAM Rev. 23, 444-466.
  3. Dines, K.A., and Lytle, R.J., 1979, Computerized geophysical tomography: Proc. IEEE 67, 1065-1073.
  4. Gilbert, P.F.C., 1972, Iterative methods for three dimensional reconstruction from projections: J. Theor. Biol. 36, 105-117.
  5. Golub, G. and Reinsch, C, 1970. Singular value decomposition and least squares solution, Num. Math.14, 403-420.
  6. Hageman, L., and Young, D., 1981, Applied iterative methods: Academic Press Inc.
  7. Hestenes, M., and Stiefel, E., 1952, Methods of conjugate gradients for solving linear systems: Nat. Bur. Standards J. Res. 49, 409-436.
  8. Ivansson, S., 1983, Remark on an earlier proposed iterative tomographic algorithm: Geophys. J. Roy. Asfr. Soc. 75, 885-860.
  9. Kaczmarz, S., 1937, Angenahert auflosung von systemen linearen gleichungen: Bull. Acad. Polon. Sci. Lett. A, 355.
  10. McMechan, G.A., 1983, Seismic tomography in boreholes: Geophys. J. Roy. Astr. Soc. 74, 601-612.
  11. Menke, W., 1989, Geophysical data analysis: Discrete inverse theory: Academic Press Inc.
  12. Meyerholtz, K.A., Pavlis, G.L., and Szpakowski, S.A., 1989, Convolutional quelling in seismic tomography: Geophysics 54, 570-580.
  13. Paige, C.C., and Saunders, M.A., 1982, LSQR: An algorithm for sparse linear equations and sparse least squares: ACM Trans. Math. Softw. 8, 43-71.
  14. Pissanetsky, S., 1984, Sparse matrix technology: Academic Press Inc.
  15. Tarantola, A., 1987, Inverse Problem Theory: Methods for Data Fitting and Model Parameter Estimation, Elsevier, Amsterdam.
  16. Trampert, J., and Leveque, J.J.. 1990, Simultaneous iterative reconstruction technique: physical interpretation based on the generalized least squares solution: J. Geophys. Rev. 95, 12553-12559.
/content/journals/10.1071/EG997369
Loading
  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error