Regularized Levenberg–Marquardt Type Method Applied to the Structural Inverse Gravity Problem in a Multilayer Medium and its Parallel Realization

Aliya F. Skurydina


The structural inverse gravity problem in a multilayer medium is one of the most important geophysics problem. Until recently, the problem was reduced to the separation of gravitational fields and the restoration of unknown layers independently. Now the methods are in demand that allow find unknown layers simultaneously. For solving Urysohn integral equation of the first kind describing the problem regularized algorithms Levenberg–Marquardt type with weight factors are investigated. A new Levenberg–Marquardt type method based on Levenberg–Marquardt scheme is proposed. A regularized Levenberg–Marquardt type method compared with classic Levenberg–Marquardt method. For classic Levenberg–Marquardt method some computational optimizations are offered. The numerical experiments using model gravitational data allow to compare convergence rates, relative errors and program execution times of classic Levenberg–Marquardt algorithm and Levenberg–Marquardt method. The parallel programs implementing the algorithms are developed using CUDA and OpenMP technologies.

Ключевые слова

Tikhonov regularization scheme; integral Urysohn type equation of first kind; regularized Levenberg-Marquardt method; regularized componentwise Levenberg-Marquardt type method; inverse gravimetry multilayer problem

Полный текст:

PDF (English)


Akimova E., Skurydina A.A. Componentwise Newton Type Method for Solving the Structural Inverse Gravity Problem. XIV EAGE International Conference — Geoinformatics: Theoretical and Applied Aspects (Kiev, Ukraine, 11–14 May, 2015). DOI: 10.3997/2214-4609.201412361.

Akimova E., Skurydina A. On Solving the Three-Dimensional Structural Gravity Problem for the Case of a Multilayered Medium by the Componentwise Levenberg–Marquardt Method. XV EAGE International Conference — Geoinformatics: Theoretical and Applied Aspects (Kiev, Ukraine, 10–13 May, 2016). pp. 181–184. DOI: 10.3997/2214-4609.201600505.

Vasin V.V., Perestoronina G.Ya. Levenberg–Marquard Method and Its Modifying Variants for Solutions to Nonlinear Equations with Applications to Inverse Gravimetry Problem. Proceedings of the Institute of Mathematics and Mechanics UrB of RAN. 2011. vol. 11, no. 2. pp. 53–61. DOI: 10.1134/S0081543813020144.

Kaltenbacher B., Neubauer A., Scherzer O. Iterative Regularization Methods for Nonlinear Ill-Posed problems. Berlin, New York. Walter de Gruyter, 2008. 194 p.

Hanke M. A Regularization Levenberg–Marquardt Scheme, with Applications to Inverse Groundwater Filtration Problems. Inverse Problems. 1997. vol. 13, no. 1. pp. 79–95.

Vasin V.V. The Levenberg–Marquardt Method for Approximation of Solutions of Irregular Operator Equations. Automation and Remote Control. 2012. vol. 73, no. 3. pp. 440–449 DOI: 10.1134/S0005117912030034.

Akimova E.N., Martyshko P.S., Misilov V.E. Algorithms for Solving the Structural Gravity Problem in a Multilayer Medium. Doklady Earth Sciences. 2013. vol. 453, no. 2. pp. 1278–1281. DOI: 10.1134/S0081543813020144.

Martyshko P.S., Prutkin I.L. Technology of Separation of Gravitational Field Sources in Depth. Geofizicheskii Zhurnal [Geophysical Journal]. 2003. vol. 25, no. 3. pp. 159–168. (in Russian)

Martyshko P.S., Akimova E.N., Misilov V.E. Solving the Structural Inverse Gravity Problem by the Modified Gradient Methods // Izvestiya, Physics of the Solid Earth. 2016. vol. 52, no. 5. pp. 704–708. DOI: 10.1134/S1069351316050098.

Akimova E.N., Martyshko P.S., Misilov V.E. A Fast Parallel Gradient Algorithm for Solving Structural Inverse Gravity Problem. XIII International Conference on Numerical Analysis and Applied Mathematics: AIP Conference Proceedings (Rodos, Greece, 22–28 September, 2015). 1648, 850063. DOI: 10.1063/1.4913118.

Akimova E.N., Misilov V.E. A Fast Componentwise Gradient Method for Solving Structural Inverse Gravity Problem. International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management (SGEM): Proceedings of 15th Intern. Multidisciplinary Scientific GeoConference SGEM. (Albena, Bulgaria, 18–24 June, 2015). vol. 3. pp. 775–782.

Vasin V.V., Eremin I.I. Operators and Iterative Processes of Fejer Type. Theory and Applications. Walter de Gruyter, Berlin. 2009. 155 p. DOI: 10.1515/9783110218190.

Bakhvalov N.S., Zhidkov N.P., Kobel'kov, G.M. Numerical Methods. Nauka, Moscow. 1987. 600 p.