2016/04/23 by V. S. Sizikov, Sizikov, Valery, Vadim Evseev +1
Earth and Planetary Sciences · #41A29 #65C20 #65F22 #86A22 #FOS: Mathematics #FOS: Physical sciences #G.1 #Geophysical and Geoelectrical Methods #Geophysics (physics.geo-ph) #Geophysics and Gravity Measurements #I.6 #J.2 #Numerical Analysis (math.NA) #Seismic Imaging and Inversion Techniques
paper · pdf · doi:10.48550/arxiv.1604.06927
openalex publication_date 2016/04/23 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/01
The direct gravimetry problem is solved by dividing each deposit body into a set of vertical adjoining bars, whereas in the inverse problem, each deposit body is modelled by a homogeneous ellipsoid of revolution (spheroid). Well-known formulae for the z-component of gravitational intensity for a spheroid are transformed to a convenient form. Parameters of a spheroid are determined by minimizing the Tikhonov smoothing functional with constraints on the parameters, which makes the ill-posed inverse problem by unique and stable. The Bulakh algorithm for initial estimating the depth and mass of a deposit is modified. The proposed technique is illustrated by numerical model examples of deposits in the form of two and five bodies. The inverse gravimetry problem is interpreted as a gravitational tomography problem or, in other words, as "introscopy" of Earth's crust and mantle.