On the Generalization of Seismic Tomography Algorithms

Mihály Dobróka, Hajnalka Szegedi

Department of Geophysics, University of Miskolc, Miskolc, Hungary.
MTA-ME Applied Geosciences Research Group, Miskolc, Hungary.

DOI: 10.4236/ajcm.2014.41004                

Abstract

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.

Keywords

Robust Tomography; Conjugated Gradients; SIRT; Cauchy-Steiner Weights

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	J. A. Scales, “Tomographic Inversion via the Conjugate Gradient Method,” Geophysics, Vol. 52, No. 2, 1987, pp. 179-185. http://dx.doi.org/10.1190/1.1442293
[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. http://dx.doi.org/10.1016/0021-9991(88)90115-5
[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. http://dx.doi.org/10.1111/j.1365-2478.1991.tb00334.x
[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