1887
Volume 23, Issue 3
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

The ray tomography problem can be reduced to a minimum norm, least-squares constrained problem. As is well known, the generalized inverse matrix methods, such as singular value decomposition (SVD) and damped least squares residual inversion (DLSQR, WLSQR), which work well for dense matrix problems, become expensive for a large, sparse problem in ray tomography. Some iterative algorithms (ART, SIRT) are directly applied to solve the system without considering a information and the non-uniqueness of the problem. The conjugate gradient (CG) method and the Paige-Saunders' algorithm are popular solutions only for the (undamped) LSQR problem. In this paper, the CG iterative algorithm is developed to solve the damped minimum norm, least squares constrained problem and the dual inversion problem in non-linear traveltime tomography. The suggested algorithms exhibit the following desirable features: (1) consideration is given to errors in both the primary estimation of the model and the traveltime data; (2) solutions are less sensitive to data error; (3) there is a reduction in the non-uniqueness of the solution using a information ('hard' and 'soft' bounds); and (4) the algorithms take advantage of the sparsitv of the matrix and consume less computer time and memory space for a large problem. The quantitative comparisons between most of the popular ray tomography methods and our algorithms were based on three measures: the misfit with the true model; solution stability; and image quality. Studies were made under different levels of noise. The results indicate that our algorithm is an accurate, flexible and efficient method for non-linear traveltime tomography and has greater tolerance to noise than any standard tomography algorithm.

© 1992 ASEG
Loading

Article metrics loading...

/content/journals/10.1071/EG992497
1992-06-01
2026-01-24

  • From This Site

    /content/journals/10.1071/EG992497
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Bregman, N., Bailey, R. and Chapman, C. (1989). ‘Crosshole seismic tomography’. Geophysics 54, 200-215.
  2. Carrion, P. M. (1989). ‘Generalized non-linear elastic inversion with constraints in model and data space’. Geophys. J. 96, 151-162.
  3. Carrion, P. M. (1991). ‘Dual tomography for imaging complex structure’. Geophysics 56, 1395-1404.
  4. Carrion, P. M. and Carneiro, D. J. (1989). ‘Generalized approach to linearized inversion with constraints’. Geophys. Res. Lett. 16, 1039-1042.
  5. Humphreys, E. and Clayton, R. W. (1988). ‘Adaptation of the tomographic method of inversion to seismic traveltime problem’. J. Geophys. Res. 93, 1073-1086.
  6. Lo, T W. (1991). ‘Fuzzy logic approach for solving a traveltime tomography problem’. SEG 1991 Annual Meeting, Houston Expanded Abstract, 925-927.
  7. McMechan, G. A., Harris, J. M. and Anderson, L. M. (1987). ‘Crosshole tomography for strongly variable media with applications to scale model data’. Bull. Seis. Soc. Am. 77, 1945-1960.
  8. Menke, W. (1984). ‘Geophysical data analysis: Discrete inverse theory’. Academic Press Inc., San Diego, California., 40-126.
  9. Meyerholtz, K. A. and Szpakowski, S. A. (1989). ‘Convolutional quelling in seismic tomography’.Geophysics 54, 570-580
  10. Moser, T. J. (1991). ‘Shortest path calculation of seismic rays’. Geophysics 56, 59-67.
  11. Nolet, G. (1987). ‘Seismic tomography’. D. Reidel Publishing Company, Dordrecht, Holland, 49-83.
  12. 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.
  13. Peterson, J. E., Paulsson, B. N. P. and McEvilly, T. V. (1985). ‘Application of algebraic reconstruction techniques to crosshole seismic data’. Geophysics 50, 1566-1580.
  14. Phillips, W. S. and Fehler, M. C. (1991). ‘Traveltime tomography: A comparison of popular methods’. Geophysics 56, 1639-1649.
  15. Scales, J. (1987). ‘Tomographic inversion via the conjugate gradient method’. Geophysics 52, 179-185.
  16. Shaw, P. R. and Orcutt, J. A. (1985). ‘Waveform inversion of seismic refraction data and applications to young Pacific crust’. Geophys. J. R. Astr. Soc. 82, 375-414.
  17. Singh, R. P. and Singh, Y. P. (1991). ‘RAYPT-a new inversion technique for geotomographic data’. Geophysics 56, 1215-1227.
  18. Tarantola, A. and Valette, B. (1982). ‘Generalized non-linear inverse problems solved using the least-squares criterion’. Rev. Geophys. Space. Phys. 20, 219-232.
  19. Zhou, B., Greenhaigh, S. A. and Sinadinovski, C. (1992). ‘Non-linear traveltime tomography: imaging high contrast inhomogeneities’.
/content/journals/10.1071/EG992497
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