Solution of the gravimetric inverse problem using multidimensional grids

Мартышко, П. С.; Ladovskii, I. V.; Бызов, Д. Д.

doi:10.1134/s1028334x13060172

Cited by 14 publications

(4 citation statements)

References 0 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…This 1D density distribution is continued outside the study volume and the excess density is calculated relative to it. As we seek 3D density distribution as a product of this function, ρ 0 (z) and some 2-dimensional corrective addition Φ(x, y), the inverse problem for lateral density in a flat layer will be sufficiently stable [23,26]. The described interpretation algorithm of the gravity data was applied to the difference of the observed field and the initial model field (Figure 5).…”

Section: Experimental Results Three-dimensional Density Model Of Midd...mentioning

confidence: 99%

Parallel Algorithms for Solving Inverse Gravimetry Problems: Application for Earth’s Crust Density Models Creation

2021

View full text Add to dashboard Cite

The paper describes a method of gravity data inversion, which is based on parallel algorithms. The choice of the density model of the initial approximation and the set on which the solution is sought guarantees the stability of the algorithms. We offer a new upward and downward continuation algorithm for separating the effects of shallow and deep sources. Using separated field of layers, the density distribution is restored in a form of 3D grid. We use the iterative parallel algorithms for the downward continuation and restoration of the density values (by solving the inverse linear gravity problem). The algorithms are based on the ideas of local minimization; they do not require a nonlinear minimization; they are easier to implement and have better stability. We also suggest an optimization of the gravity field calculation, which speeds up the inversion. A practical example of interpretation is presented for the gravity data of the Urals region, Russia.

show abstract

Section: Experimental Results Three-dimensional Density Model Of Midd...mentioning

confidence: 99%

Parallel Algorithms for Solving Inverse Gravimetry Problems: Application for Earth’s Crust Density Models Creation

2021

View full text Add to dashboard Cite

show abstract

“…Input data for initial 3D density model contains 3 parts: a) profile hodographs, time fields and appropriate 2d velocity models, b) empiric correlation between density and velocity and c) gravity field anomalies digital maps in Bougier reduction. The main stages of proposed methods are [Martyshko et al, 2010[Martyshko et al, , 2013[Martyshko et al, , 2016 construction of velocity sections of the Earth's crust, refinement of the coefficients of the regression "velocity-density" dependence from the results of 2D gravitational modeling for the given region, construction of 3D model of the initial approximation, calculation of difference between observed and modelling field, extracting field from layers, density values calculating by method of local corrections with adaptive regularization in each layer. It is known that correlation between density and longitudinal wave velocity for Earth crust rocks can be represented as piece-wise linear regression dependence.…”

Section: Methodsmentioning

confidence: 99%

“…Layered velocity distribution in a form of grid arrays fits initial model ideally. This provides stability of inverse gravity problem in a class of weak-unique solutions for models with inhomogeneous layers [Martyshko et al, 2013]. The stable algorithm of layer-by-layer inversion use 2D corrective additions with zero average value [Martyshko et al, 2016].…”

Section: Methodsmentioning

confidence: 99%

See 1 more Smart Citation

Structural Tectonic Scheme Creation Based on Seismic-Gravity Models and Isostasy Usage: Ural Case

Martyshko,

Ladovskii,

Byzov

et al. 2024

Russian Journal of Earth Sciences

View full text Add to dashboard Cite

Process of Earth’s density models creation leads to the solution of direct and inverse gravimetry problems. The inverse problem of gravimetry is a classic example of an ill-posed problem: in the common statement, its solution is not unique and unstably depends on input data. Therefore, it is necessary to determine solutions belonging substantial sets of correctness, choosing reasonable models of an initial approximation. In this paper the application of complex interpretation methods of seismic and gravitational data for the creation of three-dimensional models of crust and the upper mantle are presented. Original algorithms and programs were developed for implementation of these methods. They contain solution of non-linear (structural) inverse problem and the solution of the linear three-dimensional inverse problem taking note of the side sources. Coefficients of the “density-velocity” correlation formulas for a number of geo-traverses were defined. Also, we suggest a technique of tectonic maps construction, which is based on the lithostatic pressure calculation. Its idea can be applied to both two- and three-dimensional cases. In the 2D case we show the way to split the mantle to blocks with vertical boundaries. If lithostatic compensation hypothesis is adopted, the method also allows one to calculate density value for each block. Such separation of the mantle helps to diminish discrepancy between model and observed fields. In 3D case we suggest a method, which can be used to construct tectonic structure maps with information about approximate depth and height of each tectonic block.

show abstract