On the Generalization of Seismic Tomography Algorithms


The seismic tomography problem often leads to underdetermined and inconsistent system of equations. Solving these problems, care must be taken to keep the propagation of data errors under control. Especially, the non-Gaussian nature of the noise distribution (for example outliers in the data sets) can cause appreciable distortions in the tomographic imaging. In order to reduce the sensitivity to outlier, some generalized tomography algorithms are proposed in the paper. The weighted version of the Conjugate Gradient method is combined with the Iteratively Reweighted Least Squares (IRLS) procedure leading to a robust tomography method (W-CGRAD). The generalized version of the SIRT method is introduced in which the (Cauchy-Steiner) weighted average of the ART corrections is used. The proposed tomography algorithms are tested for a small sized tomography example by using synthetic traveltime data. It is proved that—compared to their traditional versions—the outlier sensitivities of the generalized tomography methods are sufficiently reduced.

Share and Cite:

M. Dobróka and H. Szegedi, "On the Generalization of Seismic Tomography Algorithms," American Journal of Computational Mathematics, Vol. 4 No. 1, 2014, pp. 37-46. doi: 10.4236/ajcm.2014.41004.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] J. A. Scales, “Tomographic Inversion via the Conjugate Gradient Method,” Geophysics, Vol. 52, No. 2, 1987, pp. 179-185.
[2] M. Dobroka, L. Dresen, C. Gelbke and H. Rüter, “Tomographic Inversion of Normalized Data: Double-Trace Tomography Algorithms,” Geophysical Prospecting, Vol. 40, No. 1, 1992, pp. 1-14. http://dx.doi.org/10.1111/j.1365-2478.1992.tb00361.x
[3] L. Amundsen, “Comparison of the Least-Squares Criterion and the Cauchy Criterion in Frequency-Wave Number Inversion,” Geophysics, Vol. 56, No. 12, 1991, pp. 2027-2038. http://dx.doi.org/10.1190/1.1443015
[4] F. Steiner, “Most Frequent Value Procedures,” Geophysical Transactions, Vol. 34, No. 2-3, 1988, pp. 139-260.
[5] A. Gersztenkorn and J. A. Scales, “Smoothing Seismic Tomograms with Alpha-Trimmed Means,” Geophysical Journal, Vol. 92, No. 1, 1988, pp. 67-72. http://dx.doi.org/10.1111/j.1365-246X.1988.tb01121.x
[6] J. A. Scales, A. Gersztenkorn and S. Treitel, “Fast Lp Solution of Large, Sparse, Linear Systems: Application to Seismic Travel Time Tomography,” Journal of Computational Physics, Vol. 75, No. 2, 1988, pp. 314-333.
[7] M. Dobroka, A. Gyulai, T. Ormos, J. Csókás and L. Dresen, “Joint Inversion of Seismic and Geoelectric Data Recorded in an Under-Ground Coal Mine,” Geophysical Prospecting, Vol. 39, No. 5, 1991, pp. 643-665.
[8] H. Szegedi and M. Dobroka, “On the Use of Steiner’s Weights in Inversion-Based Fourier Transformation,” Acta Geodaetica et Geophysica, 2014, pp. 1-10. http://dx.doi.org/10.1007/s40328-014-0041-0

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.