Abstract：
Tikhonov regularization (TR) method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method (EXTR) is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.
WANG Zhu-Wen,XU Shi,LIU Yin-Ping et al. Extrapolated Tikhonov method and inversion of 3D density images of gravity data[J]. APPLIED GEOPHYSICS, 2014, 11(2): 139-148.
[1]
Barbosa, V. C. F., and Silva, J. B. C., 1994, Generalized compact gravity inversion: Geophysics, 59(1), 57-68.
[2]
Blakely, R. J., 1995, Potential theory in gravity and magnetic applications: Cambridge University Press.
[3]
Chai, Y., and Hinze, W. J., 1988, Gravity inversion of an interface above which the density contrast varies exponentially with depth: Geophysics, 53(6), 837-845.
[4]
Chenot, D., and Debeglia, N., 1990, Three-dimensional gravity or magnetic constrained depth inversion with lateral and vertical variation of contrast: Geophysics, 55(3), 327-335.
[5]
Commer, M., 2011, Three-dimensional gravity modeling and focusing inversion using rectangular meshes: Geophysical Prospecting, 59, 966-979.
[6]
Cordell, L., and Henderson, R. G., 1968, Iterative three-dimensional solution of gravity anomaly data using a digital computer: Geophysics, 33(4), 596-601.
[7]
Guillen, A,, and Menichetti, V., 1984, Gravity and magnetic inversion with minimization of a specific functional: Geophysics, 49(8), 1354-1360.
[8]
Guo, L. H., Meng, X. H., and Shi, L., 2009, 3D correlation imaging for gravity and gravity gradiometry data: Chinese J. Geophys. (in Chinese), 52(4), 1092-1106.
[9]
Hamarik, U., Palm, R., and Raus, T., 2008, Extrapolation of Tikhonov and Lavrentiev regularization methods: Journal of Physics: Conference Series, 135, 1-8.
[10]
Hamarik, U., Palm, R., and Raus, T., 2007, Use of extrapolation in regularization methods: Journal of Inverse and Ill-Posed Problems, 15, 277-294.
[11]
Joerg L., and Heinrich V., 2013, Large-scale Tikhonov regularization of total least squares: Journal of Computational and Applied Mathematics, 238, 95-108
[12]
Last, B. J., and Kubik, K.,1983, Compact gravity inversion: Geophysics, 48(6), 713-721.
[13]
Lawson, C. J., and Hanson, R. J., 1974, Solving least squares problems: Prentice-Hall Inc.
[14]
Li, Y. G., and Oldenburg, D. W., 2003, Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method: Geophysical Journal International, 152, 251-265.
[15]
Li, Y. G., and Oldenburg, D. W., 1996, 3-D inversion of magnetic data: Geophysics, 61(2), 394-408.
[16]
Li, Y. G., and Oldenburg, D. W., 1998, 3-D inversion of gravity data: Geophysics, 63(1), 109-119.
[17]
Liu, Y. P., Wang, Z. W., Du, X. J., et al., 2013, 3D constrained inversion of gravity data based on Extrapolation Tikhonov regularization algorithm: Chinese J. Geophys. (in Chinese), 56(5), 1650-1659.
[18]
Medeiros, W. O., and Silva, J. B. C., 1996, Geophysical inversion using approximate equality constraints: Geophysics, 61(6), 1678-1688.
[19]
Meng, X. H., Liu, G. F., and Chen, Z. X., et al., 2012, 3-D gravity and magnetic inversion for physical properties based on residual anomaly correlation: Chinese J. Geophys. (in Chinses), 55(1), 304-309.
[20]
Pedersen, L. B., 1977, Interpretation of potential field data, a generalized inverse approach: Geophys. Pros ., 25, 199-230.
[21]
Silva, J. B. C., and Barbosa, V., 2006, Interactive gravity inversion: Geophysics, 71, J1-J9.
[22]
Silvia, G., and Paolo, N., 2014, Automatic parameter setting for Arnoldi-Tikhonov methods: Journal of Computational and Applied Mathematics, 256, 180-195
[23]
Tikhonov, A. N., 1963, Regularization of ill-posed problems: Doklady Akad. Nauk, SSSR, 153, 49-52.
[24]
Tikhonov, A. N., and Arsenin, V. Y., 1977, Solution of ill-posed problems: V. H. Winston &Sons.
[25]
Yao, C. L., Zheng, Y. M., and Zhang, Y. W., 2007, 3D gravity and magnetic inversion for physical properties using stochastic subspaces: Chinese J. Geophys. (in Chinese), 50(5), 1576-1583.
[26]
Yao, C. L., Hao, T. Y., and Guan, Z. N., et al, 2003, High-speed computation and effective storage in 3-D gravity and magnetic inversion based on genetic algorithms: Chinese J. Geophys. (in Chinses), 46(2), 252-258.
[27]
Zhdanov, M., 1988, Integral Transforms in Geophysics: Springer Verlag.
[28]
Zheng H. L., 2005, Gravity field and gravity exploration: Geological Publishing House, BeiJing.
2014, Vol. 11 
