Gravity Data Inversion with Method of Local Corrections for Finite Elements Models
We present a new method for gravity data inversion for the linear problem (reconstruction of density distribution by given gravity field). This is an iteration algorithm based on the ideas of local minimization (also known as local corrections method). Unlike the gradient methods, it does not requir...
Published in: | Geosciences |
---|---|
Main Authors: | , , , |
Format: | Text |
Language: | English |
Published: |
Multidisciplinary Digital Publishing Institute
2018
|
Subjects: | |
Online Access: | https://doi.org/10.3390/geosciences8100373 |
id |
ftmdpi:oai:mdpi.com:/2076-3263/8/10/373/ |
---|---|
record_format |
openpolar |
spelling |
ftmdpi:oai:mdpi.com:/2076-3263/8/10/373/ 2023-08-20T04:09:07+02:00 Gravity Data Inversion with Method of Local Corrections for Finite Elements Models Petr S. Martyshko Igor V. Ladovskii Denis D. Byzov Alexander G. Tsidaev agris 2018-10-10 application/pdf https://doi.org/10.3390/geosciences8100373 EN eng Multidisciplinary Digital Publishing Institute Geophysics https://dx.doi.org/10.3390/geosciences8100373 https://creativecommons.org/licenses/by/4.0/ Geosciences; Volume 8; Issue 10; Pages: 373 inverse problem 3D density model joint interpretation Text 2018 ftmdpi https://doi.org/10.3390/geosciences8100373 2023-07-31T21:46:17Z We present a new method for gravity data inversion for the linear problem (reconstruction of density distribution by given gravity field). This is an iteration algorithm based on the ideas of local minimization (also known as local corrections method). Unlike the gradient methods, it does not require a nonlinear minimization, is easier to implement and has better stability. The algorithm is based on the finite element method. The finite element approach in our study means that the medium (part of a lithosphere) is represented as a set of equal rectangular prisms, each with constant density. We also suggest a time-efficient optimization, which speeds up the inversion process. This optimization is applied on the gravity field calculation stage, which is a part of every inversion iteration. Its idea is to replace multiple calculations of the gravity field for all finite elements in all observation points with a pre-calculated set of uniform fields for all distances between finite element and observation point, which is possible for the current data set. Method is demonstrated on synthetic data and real-world cases. The case study area is located on the Timan-Pechora plate. This region is one of the promising oil- and gas-producing areas in Russia. Note that in this case we create a 3D density model using joint interpretation of seismic and gravity data. Text Pechora MDPI Open Access Publishing Geosciences 8 10 373 |
institution |
Open Polar |
collection |
MDPI Open Access Publishing |
op_collection_id |
ftmdpi |
language |
English |
topic |
inverse problem 3D density model joint interpretation |
spellingShingle |
inverse problem 3D density model joint interpretation Petr S. Martyshko Igor V. Ladovskii Denis D. Byzov Alexander G. Tsidaev Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
topic_facet |
inverse problem 3D density model joint interpretation |
description |
We present a new method for gravity data inversion for the linear problem (reconstruction of density distribution by given gravity field). This is an iteration algorithm based on the ideas of local minimization (also known as local corrections method). Unlike the gradient methods, it does not require a nonlinear minimization, is easier to implement and has better stability. The algorithm is based on the finite element method. The finite element approach in our study means that the medium (part of a lithosphere) is represented as a set of equal rectangular prisms, each with constant density. We also suggest a time-efficient optimization, which speeds up the inversion process. This optimization is applied on the gravity field calculation stage, which is a part of every inversion iteration. Its idea is to replace multiple calculations of the gravity field for all finite elements in all observation points with a pre-calculated set of uniform fields for all distances between finite element and observation point, which is possible for the current data set. Method is demonstrated on synthetic data and real-world cases. The case study area is located on the Timan-Pechora plate. This region is one of the promising oil- and gas-producing areas in Russia. Note that in this case we create a 3D density model using joint interpretation of seismic and gravity data. |
format |
Text |
author |
Petr S. Martyshko Igor V. Ladovskii Denis D. Byzov Alexander G. Tsidaev |
author_facet |
Petr S. Martyshko Igor V. Ladovskii Denis D. Byzov Alexander G. Tsidaev |
author_sort |
Petr S. Martyshko |
title |
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
title_short |
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
title_full |
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
title_fullStr |
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
title_full_unstemmed |
Gravity Data Inversion with Method of Local Corrections for Finite Elements Models |
title_sort |
gravity data inversion with method of local corrections for finite elements models |
publisher |
Multidisciplinary Digital Publishing Institute |
publishDate |
2018 |
url |
https://doi.org/10.3390/geosciences8100373 |
op_coverage |
agris |
genre |
Pechora |
genre_facet |
Pechora |
op_source |
Geosciences; Volume 8; Issue 10; Pages: 373 |
op_relation |
Geophysics https://dx.doi.org/10.3390/geosciences8100373 |
op_rights |
https://creativecommons.org/licenses/by/4.0/ |
op_doi |
https://doi.org/10.3390/geosciences8100373 |
container_title |
Geosciences |
container_volume |
8 |
container_issue |
10 |
container_start_page |
373 |
_version_ |
1774721850941636608 |