Projected Barzilai-Borwein method for the acceleration of gravity field data inversion

2020 ◽  
pp. 1-12
Author(s):  
Zhaohai Meng ◽  
Wei Li ◽  
Fengting Li ◽  
Hongyu Li
2021 ◽  
Author(s):  
Bart Root ◽  
Javier Fullea ◽  
Jörg Ebbing ◽  
Zdenek Martinec

<p>Global gravity field data obtained by dedicated satellite missions is used to study the density distribution of the lithosphere. Different multi-data joint inversions are using this dataset together with other geophysical data to determine the physical characteristics of the lithosphere. The gravitational signal from the deep Earth is usually removed by high-pass filtering of the model and data, or by appropriately selecting insensitive gravity components in the inversion. However, this will remove any long-wavelength signal inherent to lithosphere. A clear choice on the best-suited approach to remove the sub-lithospheric gravity signal is missing. </p><p>Another alternative is to forward model the gravitational signal of these deep situated mass anomalies and subtract it from the observed data, before the inversion. Global tomography provides shear-wave velocity distribution of the mantle, which can be transformed into density anomalies. There are difficulties in constructing a density model from this data. Tomography relies on regularisation which smoothens the image of the mantle anomalies. Also, the shear-wave anomalies need to be converted to density anomalies, with uncertain conversion factors related to temperature and composition. Understanding the sensitivity of these effects could help determining the interaction of the deep Earth and the lithosphere.</p><p>In our study the density anomalies of the mantle, as well as the effect of CMB undulations, are forward modelled into their gravitational potential field, such that they can be subtracted from gravity observations. The reduction in magnitude of the density anomalies due to the regularisation of the global tomography models is taken into account. The long-wavelength region of the density estimates is less affected by the regularisation and can be used to fix the mean conversion factor to transform shear wave velocity to density. We present different modelling approaches to add the remaining dynamic topography effect in lithosphere models. This results in new solutions of the density structure of the lithosphere that both explain seismic observations and gravimetric measurements. The introduction of these dynamic forces is a step forward in understanding how to properly use global gravity field data in joint inversions of lithosphere models.</p>


Geophysics ◽  
1971 ◽  
Vol 36 (3) ◽  
pp. 605-608 ◽  
Author(s):  
Edwin S. Robinson

Investigation of geological structure by gravimetric and magnetic field surveys requires consideration of relationships between gravity anomaly and magnetic anomaly generating sources. The possibility of using Poisson’s Relation to examine magnetic and gravity fields related to a common source is intriguing. This relation is expressed as follows: [Formula: see text] (1) where A (x, y, z) is the magnetic field potential and U (x, y, z) is the gravity field potential at a point in space due to a source of uniform density ρ and uniform magnetization I in the direction α. This expression has been used to derive magnetic anomalies over idealized forms (Nettleton, 1940) and, by Baranov (1957), to extract pseudogravity fields from magnetic field data. The purpose of this paper is to develop an expression for extracting a pseudomagnetic field from gravity field data and to examine the practical applications of this expression.


2020 ◽  
Vol 222 (2) ◽  
pp. 1074-1089 ◽  
Author(s):  
Yajun Liu ◽  
Pritam Yogeshwar ◽  
Xiangyun Hu ◽  
Ronghua Peng ◽  
Bülent Tezkan ◽  
...  

SUMMARY Electrical anisotropy of formations has been long recognized by field and laboratory evidence. However, most interpretations of long-offset transient electromagnetic (LOTEM) data are based on the assumption of an electrical isotropic earth. Neglecting electrical anisotropy of formations may cause severe misleading interpretations in regions with strong electrical anisotropy. During a large scale LOTEM survey in a former mining area in Eastern Germany, data was acquired over black shale formations. These black shales are expected to produce a pronounced bulk anisotropy. Here, we investigate the effects of electrical anisotropy on LOTEM responses through numerical simulation using a finite-volume time-domain (FVTD) algorithm. On the basis of isotropic models obtained from LOTEM field data, various anisotropic models are developed and analysed. Numerical results demonstrate that the presence of electrical anisotropy has a significant influence on LOTEM responses. Based on the numerical modelling results, an isolated deep conductive anomaly presented in the 2-D isotropic LOTEM electric field data inversion result is identified as a possible artifact introduced by using an isotropic inversion scheme. Trial-and-error forward modelling of the LOTEM electric field data using an anisotropic conductivity model can explain the data and results in a reasonable quantitative data fit. The derived anisotropic 2-D model is consistent with the prior geological information.


Geophysics ◽  
1953 ◽  
Vol 18 (4) ◽  
pp. 907-909 ◽  
Author(s):  
L. J. Peters ◽  
T. A. Elkins

We would like to call attention to a point of considerable practical importance which is neglected in this interesting and ingenious paper on the computation of the second derivative. This is the fact that gravity field data inevitably contain errors so that the second derivative values computed by coefficients from this gravity data also will contain errors, which may be of such magnitude as to mask the real effects caused by geologic structure, the finding of which was the purpose of the gravity survey.


1987 ◽  
Author(s):  
Roman Alvarez ◽  
Alfredo Cortés
Keyword(s):  

Geophysics ◽  
2012 ◽  
Vol 77 (4) ◽  
pp. WB219-WB231 ◽  
Author(s):  
P. Kaikkonen ◽  
S. P. Sharma ◽  
S. Mittal

Three-dimensional linearized nonlinear electromagnetic inversion is developed for revealing the subsurface conductivity structure using isolated very low frequency (VLF) and VLF-resistivity anomalies due to conductors that may be arbitrarily directed towards the measuring profiles and the VLF transmitter. We described the 3D model using a set of variables in terms of geometric and physical parameters. These model parameters were then optimized (parametric inversion) to obtain their best estimates to fit the observations. Two VLF transmitters, i.e., the [Formula: see text], [Formula: see text] (“E”) and the [Formula: see text], [Formula: see text] (“H”) polarizations, respectively, can be considered jointly in inversion. After inverting several noise-free and noisy synthetic data, the results revealed that the estimated model parameters and the functionality of the approach were very good and reliable. The inversion procedure also worked well for the field data. The reliability and validity of the results after the field data inversion have been checked using data from a shear zone associated with uranium mineralization.


2018 ◽  
Vol 2018 (1) ◽  
pp. 1-4
Author(s):  
Zhaohai Meng ◽  
Fengting Li ◽  
Hao Yu ◽  
Lin Ma ◽  
Zhongli Li

Geosciences ◽  
2018 ◽  
Vol 8 (10) ◽  
pp. 373 ◽  
Author(s):  
Petr Martyshko ◽  
Igor Ladovskii ◽  
Denis Byzov ◽  
Alexander Tsidaev

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.


Sign in / Sign up

Export Citation Format

Share Document