scholarly journals Imitation modelling technology for gravity inversion cases

2019 ◽  
Keyword(s):  
Gravity Field ◽  
Gravity Data ◽  
A Priori ◽  
Mineral Resources ◽  
Structural Type ◽  
Practical Significance ◽  
Gravity Inversion ◽  
Density Structure ◽  
Research Areas ◽  
Computer Technologies

Formulation of the problem. A gravity method is aimed at prospecting and exploration of mineral resources which are based on the study of the geological section structure. The task of quantitative interpretation of the gravimetric materials, which uses methods for solving direct and inverse gravity problems, is the modelling of a gravity field (direct problem) and geological media’s density structure (inverse problem). The important features of methods for density structure modelling of complex geological media are geological content, consistency with a priori data and its subordination to geological hypotheses. It is proposed to analyze these properties by a imitation technique. The purpose of the article is to describe the imitation gravimetric modelling method, based on the construction of an informal sequence of equivalent solutions. The purpose of imitation modelling is to study the properties of gravity inversion in general formulation as well as to assess the degree of detail and reliability of the methodology and technologies of gravity modelling, which is claimed to be an effective solution to geological problems. Methods. Imitation modelling technology and methods of solving gravity direct and inverse problems for geodensity model of complex geological environment. Results. The examples of density and structural simulation testing of the informal sequence of equivalent solutions and its computer technologies show that complex interpretation of wells, seismic and gravity data enables to create detailed density models of geological medium. Studies have also been conducted of ways to increase the reliability of gravitational modelling. Scientific novelty and practical significance. It is revealed that the best approximation of the regional background is an inclined plane, which approximates the observed gravity field along characteristic pickets over the research areas that are better studied. Also, an increase in the reliability of modelling can be achieved by rebuilding near side zones in structural type models in an interactive process of solving structural inverse gravity problems. Substantive modelling depends primarily on the experience of the interpreter, since computer technologies for solving direct and inverse gravity problems are only an interpretation tool.

scholarly journals SIMULATION MODELLING IN THE STRUCTURAL GRAVITY PROSPECTING

2019 ◽  
pp. 38-48
Author(s):  
S. H. Anikeyev ◽  
S. M. Bahriy ◽  
B. B. Hablovskiy
Keyword(s):  
Cross Sections ◽  
Gravity Data ◽  
A Priori ◽  
Structural Type ◽  
Simulation Modelling ◽  
Gravity Modeling ◽  
Gravity Inversion ◽  
Density Structure ◽  
Data Modelling ◽  
Computer Technologies

In accordance with the purpose of geophysical exploration, the gravity data interpretation is aimed at prospecting mineral resources which is based on the study of the geological cross-section structure. The task of quantitative interpretation, which uses methods of gravity modeling and gravity inversion, is the modelling of a gravity field (gravity modeling) and of a density structure of geological environments (gravity inversion). The article presents the definition and steps of the gravity data modelling technique. This technique is based on the construction of an informal sequence of equivalent solutions. The technological and geological features of methods for modelling the density structure of complex geological environments are given; among them geological content, consistency with a priori data and the subordination of modelling to geological hypotheses are important. The topicality and methods of simulation modelling are outlined. The purpose of simulation modelling is to study the properties of gravity inversion in the general formulation, as well as to evaluate the degree of detail and reliability of the methods and technologies of gravity modelling, which claim to be an effective solution to geological problems. The example of structural simulation testing of the methods of informal sequence of equivalent solutions and its computer technologies shows that a complex interpretation of seismic and gravity measurements data enables the creation of detailed density models of structural cross-sections. The ways of increasing the veracity of gravity data modelling of structural cross-sections have been studied. It is revealed that the best approximation of the regional background is an inclined plane which approximates the observed field of gravity according to characteristic pickets over the research areas that are better studied. The increase in the veracity of modeling can also be achieved by rebuilding the near side zones in the structural type models in an interactive process of solving structural gravity inversion problems. Substantive modeling depends primarily on the experience of the interpreter since computer technologies for gravity modeling and gravity inversion are merely an interpretation tool.

scholarly journals P-wave velocity and density structure beneath Mt. Vesuvius: a magma body in the upper edifice?

2013 ◽  
Vol 56 (4) ◽  
Author(s):  
Paolo Capuano ◽  
Guido Russo ◽  
Roberto Scarpa
Keyword(s):  
High Resolution ◽  
Wave Velocity ◽  
Gravity Data ◽  
High Energy ◽  
P Wave ◽  
Gravity Inversion ◽  
Density Structure ◽  
Resolution Image ◽  
High Resolution Image ◽  
P Wave Velocity

<p>A high-resolution image of the compressional wave velocity and density structure in the shallow edifice of Mount Vesuvius has been derived from simultaneous inversion of travel times and hypocentral parameters of local earthquakes and from gravity inversion. The robustness of the tomography solution has been improved by adding to the earthquake data a set of land based shots, used for constraining the travel time residuals. The results give a high resolution image of the P-wave velocity structure with details down to 300-500 m. The relocated local seismicity appears to extend down to 5 km depth below the central crater, distributed into two clusters, and separated by an anomalously high Vp region positioned at around 1 km depth. A zone with high Vp/Vs ratio in the upper layers is interpreted as produced by the presence of intense fluid circulation alternatively to the interpretation in terms of a small magma chamber inferred by petrologic studies. In this shallower zone the seismicity has the minimum energy, whilst most of the high-energy quakes (up to Magnitude 3.6) occur in the cluster located at greater depth. The seismicity appears to be located along almost vertical cracks, delimited by a high velocity body located along past intrusive body, corresponding to remnants of Mt. Somma. In this framework a gravity data inversion has been performed to study the shallower part of the volcano. Gravity data have been inverted using a method suitable for the application to scattered data in presence of relevant topography based on a discretization of the investigated medium performed by establishing an approximation of the topography by a triangular mesh. The tomography results, the retrieved density distribution, and the pattern of relocated seismicity exclude the presence of significant shallow magma reservoirs close to the central conduit. These should be located at depth higher than that of the base of the hypocenter volume, as evidenced by previous studies.</p>

Three‐dimensional gravity inversion using simulated annealing: Constraints on the diapiric roots of allochthonous salt structures

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1438-1449 ◽  
Author(s):  
Seiichi Nagihara ◽  
Stuart A. Hall
Keyword(s):  
Simulated Annealing ◽  
Gulf Of Mexico ◽  
Gravity Data ◽  
A Priori ◽  
Gravity Inversion ◽  
Inversion Method ◽  
A Priori Information ◽  
Final Density ◽  
Smoothing Effects ◽  
Priori Information

In the northern continental slope of the Gulf of Mexico, large oil and gas reservoirs are often found beneath sheetlike, allochthonous salt structures that are laterally extensive. Some of these salt structures retain their diapiric feeders or roots beneath them. These hidden roots are difficult to image seismically. In this study, we develop a method to locate and constrain the geometry of such roots through 3‐D inverse modeling of the gravity anomalies observed over the salt structures. This inversion method utilizes a priori information such as the upper surface topography of the salt, which can be delineated by a limited coverage of 2‐D seismic data; the sediment compaction curve in the region; and the continuity of the salt body. The inversion computation is based on the simulated annealing (SA) global optimization algorithm. The SA‐based gravity inversion has some advantages over the approach based on damped least‐squares inversion. It is computationally efficient, can solve underdetermined inverse problems, can more easily implement complex a priori information, and does not introduce smoothing effects in the final density structure model. We test this inversion method using synthetic gravity data for a type of salt geometry that is common among the allochthonous salt structures in the Gulf of Mexico and show that it is highly effective in constraining the diapiric root. We also show that carrying out multiple inversion runs helps reduce the uncertainty in the final density model.

scholarly journals Seismic Constrained Gravity Inversion: A Reliable Tool to Improve Geophysical Models Away from Seismic Information

Geosciences ◽  
2021 ◽  
Vol 11 (11) ◽  
pp. 467
Author(s):  
Daniele Sampietro ◽  
Martina Capponi
Keyword(s):  
High Resolution ◽  
Gravity Field ◽  
A Priori ◽  
Geophysical Methods ◽  
Real Life ◽  
Gravity Inversion ◽  
Seismic Surveys ◽  
Real Model ◽  
Gravity Fields ◽  
Priori Information

The exploitation of gravity fields in order to retrieve information about subsurface geological structures is sometimes considered a second rank method, in favour of other geophysical methods, such as seismic, able to provide a high resolution detailed picture of the main geological horizons. Within the current work we prove, through a realistic synthetic case study, that the gravity field, thanks to the availability of freely of charge high resolution global models and to the improvements in the gravity inversion methods, can represent a valid and cheap tool to complete and enhance geophysical modelling of the Earth’s crust. Three tests were carried out: In the first one a simple two-layer problem was considered, while in tests two and three we considered two more realistic scenarios in which the availability on the study area of constraints derived from 3D or 2D seismic surveys were simulated. In all the considered test cases, in which we try to simulate real-life scenarios, the gravity field, inverted by means of an advanced Bayesian technique, was able to obtain a final solution closer to the (simulated) real model than the assumed a priori information, typically halving the uncertainties in the geometries of the main geological horizons with respect to the initial model.

Isostatic constraint for 2D nonlinear gravity inversion on rifted margins

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. G17-G34
Author(s):  
B. Marcela S. Bastos ◽  
Vanderlei C. Oliveira Jr.
Keyword(s):  
Gravity Data ◽  
A Priori ◽  
Petroleum Industry ◽  
Synthetic Data ◽  
Gravity Inversion ◽  
Rifted Margins ◽  
Isostatic Equilibrium ◽  
Priori Information ◽  
Pelotas Basin ◽  
Nonlinear Gravity

We have developed a nonlinear gravity inversion for simultaneously estimating the basement and Moho geometries, as well as the depth of the reference Moho along a profile crossing a passive rifted margin. To obtain stable solutions, we impose smoothness on basement and Moho, force them to be close to previously estimated depths along the profile and also impose local isostatic equilibrium. Different from previous methods, we evaluate the information of local isostatic equilibrium by imposing smoothness on the lithostatic stress exerted at depth. Our method delimits regions that deviate and those that can be considered in local isostatic equilibrium by varying the weight of the isostatic constraint along the profile. It also allows controlling the degree of equilibrium along the profile, so that the interpreter can obtain a set of candidate models that fit the observed data and exhibit different degrees of isostatic equilibrium. Our method also differs from earlier studies because it attempts to use isostasy for exploring (but not necessarily reducing) the inherent ambiguity of gravity methods. Tests with synthetic data illustrate the effect of our isostatic constraint on the estimated basement and Moho reliefs, especially at regions with pronounced crustal thinning, which are typical of passive volcanic margins. Results obtained by inverting satellite data over the Pelotas Basin, a passive volcanic margin in southern Brazil, agree with previous interpretations obtained independently by combining gravity, magnetic, and seismic data available to the petroleum industry. These results indicate that combined with a priori information, simple isostatic assumptions can be very useful for interpreting gravity data on passive rifted margins.

scholarly journals NEW ALGORITHM FOR INVERSION OF GRAVITATIONAL DATA

2018 ◽  
pp. 51-56
Author(s):  
P. Gryshchuk
Keyword(s):  
Gravity Data ◽  
Gravity Anomalies ◽  
Practical Significance ◽  
Gravity Inversion ◽  
Section Data ◽  
Number Of Layers ◽  
Coal Deposits ◽  
Mass Displacement ◽  
Definition Of ◽  
Gravity Acceleration

The article is devoted to the description of the new algorithm for inversion of gravitational anomalies.The principle of the geological objects mass displacement is used for definition of its geometry. In the grid model, the movement of blocks is performed in four directions for a two-dimensional model. The motion and changing of density is determined in a random style. This algorithm defined the geometry of three rectangular bodies of model. The conditions for determining the exact form of gravity sources by the proposed approach were clarified. An analysis of the inversion of gravity acceleration anomalies from three bodies for sections, which consisted of three, four and five layers, was performed. The bodies of identical density were determined exactly for three and four layers, and with errors for the five layers. The model with two densities had some errors in determining the geometry of the bodies. Scientific novelty is related to the development of new algorithm, which is based on the principle of movingfor density between blocks. The new approach fulfils gravity inversion much faster than the genetic algorithm. The practical significance is determination of the number of layers needed for a proper inversion of gravity anomalies. The main factors influencing the accuracy of geometry are the data of body density and the number of layers. The method was used for the interpretation of gravity data over the coal deposits of the Donbas. The model consisted of five layers and a range of minimum to their maximum densities for the geological section. Data of the inversion was determined by the depth to the surface of the limestone, which is confirmed by the data of the drilling. The developed approach determines the geometry of section from density, which is important for geological, search and environmental tasks.

Geophysical subsurface modelling based on the updated, enhanced regional gravity field solution in Antarctica

2021 ◽  
Author(s):  
Theresa Schaller ◽  
Mirko Scheinert ◽  
Philipp Zingerle ◽  
Roland Pail ◽  
Martin Willberg
Keyword(s):  
Gravity Field ◽  
Gravity Data ◽  
Average Density ◽  
Weddell Sea ◽  
Moho Depth ◽  
Gravity Inversion ◽  
Satellite Gravity ◽  
German Research Foundation ◽  
Field Solution ◽  
Gravity Disturbances

&lt;p&gt;The gravity field reflects mass inhomogeneities (mostly) inside the Earth. Therefore, gravity inversion and geophysical gravity field modelling are important tools to study the Earth's inner structure and tectonic evolution. In Antarctica, it is extremely challenging to carry out geoscientific studies due to its harsh environment and difficult logistics. Additionally, the continent is covered by an up to 5 km thick ice sheet. However, in the framework of IAG Subcommission 2.4f &amp;#8220;Gravity and Geoid in Antarctica&amp;#8221; (AntGG) a large database of airborne, shipborne and ground based gravity data has been compiled. Especially airborne data have been acquired during recent years, among others in the area of the polar gap of satellite gravity data. Now, in a joint project funded by the German Research Foundation (DFG) all existing and new gravity data were processed to infer an enhanced gravity field solution for Antarctica (see contribution by Scheinert et al., session G1.5). Processed data e.g. gravity disturbances on the ground or a constant height and other functionals will be provided on a regular grid with 5 km grid spacing. Subsequently, the new Antarctic gravity field solution can now be used for further geophysical and tectonic studies. We use the newly calculated gravity disturbances to study subglacial topography, sediment thickness and Moho depth and to improve respective existing models in Antarctica. For this, we apply 2D Parker-Oldenburg inversion in combination with results from other gravity based studies and further constraining data (e.g. seismic data and ice penetrating radar). We investigate how the higher resolution (5 km) of the new Antarctic gravity field solution facilitates the study of smaller regions in more detail, specifically parts of Wilkes Land, Dronning Maud Land and the Weddell Sea. Additionally, we will infer accuracy estimates for the resulting boundaries in terms of the used inversion parameters (density contrast, average density and filter wavelengths) and their respective gravity signal. Thus, the challenges of gravity field inversion in Antarctica will be discussed in detail and first results of the subsurface modelling will be presented.&lt;/p&gt;

An integrated geophysical approach for imaging subbasalt sedimentary basins: Case study of Jam River Basin, India

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. B141-B147 ◽  
Author(s):  
V. Chakravarthi ◽  
G. B. K. Shankar ◽  
D. Muralidharan ◽  
T. Harinarayana ◽  
N. Sundararajan
Keyword(s):  
Gravity Field ◽  
River Basin ◽  
Gravity Data ◽  
Sedimentary Basins ◽  
Gravity Inversion ◽  
Data Points ◽  
3D Gravity ◽  
Regional Gravity ◽  
Godavari Basin

An integrated geophysical strategy comprising deep electrical resistivity and gravity data was devised to image subbasalt sedimentary basins. A 3D gravity inversion was used to determine the basement structure of the Permian sediments underlying the Cretaceous formation of the Jam River Basin in India. The thickness of the Cretaceous formation above the Permian sediments estimated from modeling 60 deep-electric-sounding data points agrees well with drilling information. The gravity effect of mass deficit between the Cretaceous and Permian formations was found using 3D forward modeling and subsequently removed from the Bouguer gravity anomaly along with the regional gravity field. The modified residual gravity field was then subjected to3D inversion to map the variations in depth of the basement beneath the Permian sediments. Inversion of gravity data resulted in two basement ridges, running almost east to west, dividing the basin into three independent depressions. It was found that the Katol and Kondhali faults were active even during post-Cretaceous time and were responsible for the development of the subsurface basement ridges in the basin. The inferred 3D basement configuration of the basin clearly brought out the listric nature of these two faults. Further, the extension of the Godavari Basin into the Deccan syneclise and the fact that the source-rock studies show the presence of hydrocarbons in the Sironcha block in the northern part of the Godavari Basin also shed some light on the hydrocarbon potential of the Jam River Basin.

Three-dimensional gravity inversion using graph theory to delineate the skeleton of homogeneous sources

Geophysics ◽  
2015 ◽  
Vol 80 (2) ◽  
pp. G53-G66 ◽  
Author(s):  
Rodrigo Bijani ◽  
Cosme F. Ponte-Neto ◽  
Dionisio U. Carlos ◽  
Fernando J. S. Silva Dias
Keyword(s):  
Genetic Algorithm ◽  
Graph Theory ◽  
Minimum Spanning Tree ◽  
Gravity Data ◽  
A Priori ◽  
Three Dimensional ◽  
Real Data ◽  
Gravity Inversion ◽  
Inversion Method ◽  
Dimensional Gravity

We developed a new strategy, based on graph theory concepts, to invert gravity data using an ensemble of simple point masses. Our method consisted of a genetic algorithm with elitism to generate a set of possible solutions. Each estimate was associated to a graph to solve the minimum spanning tree (MST) problem. To produce unique and stable estimates, we restricted the position of the point masses by minimizing the statistical variance of the distances of an MST jointly with the data-misfit function during the iterations of the genetic algorithm. Hence, the 3D spatial distribution of the point masses identified the skeleton of homogeneous gravity sources. In addition, our method also gave an estimation of the anomalous mass of the source. So, together with the anomalous mass, the skeleton could aid other 3D methods with promising geometric a priori parameters. Several tests with different values of regularizing parameter were made to bespeak this new regularizing strategy. The inversion results applied to noise-corrupted synthetic gravity data revealed that, regardless of promising starting models, the estimated distribution of point masses and the anomalous mass offered valuable information about the homogeneous sources in the subsurface. Tests on real data from a portion of Quadrilátero Ferrífero, Minas Gerais state, Brazil, were performed for complementary analysis of the proposed inversion method.

Antarctica crustal model by means of the Bayesian gravity inversion

2020 ◽  
Author(s):  
Martina Capponi ◽  
Daniele Sampietro
Keyword(s):  
Gravity Field ◽  
Gravity Data ◽  
Mass Density ◽  
Geophysical Methods ◽  
Gravity Inversion ◽  
Inversion Algorithm ◽  
Global Gravity ◽  
Uniqueness Of The Solution ◽  
Definition Of ◽  
Regional Gravity

&lt;p&gt;The Antarctica crustal structure is still not completely unveiled due to the presence of thick ice shields all over the continent which prevent direct in situ measurements. In the last decades, various geophysical methods have been used to retrieve information of the upper crust and sediments distribution however there are still regions, especially in central Antarctica, where our knowledge is limited. For these kind of situations, in which the indirect investigation of the subsurface is the only valuable solution, the gravity data are an important source of information. After the recent dedicated satellite missions, like GRACE and GOCE, it has been possible to obtain global gravity field data with spatial resolution and accuracy almost comparable to those of local/regional gravity acquisitions, paving the way to new geophysical applications. The new challenge today is the capability to invert such gravity data on large areas with the aim to obtain a 3D density model of the Earth crust. This is in fact a problem characterized by intrinsic instability and non-uniqueness of the solution that to be solved requires the definition of strong constrains and numerical regularization.&lt;/p&gt;&lt;p&gt;In this work the authors propose the application of a Bayesian inversion algorithm to the Antarctica continent to infer a model of mass density distribution. The first operation is the definition of an initial geological model in terms of geological horizons and density. These two variables are considered as random variables and, within the iterative procedure based on Markov Chain Monte Carlo methods, they are adjusted in such a way to fit the gravity field on the surface. The test performed show that the method is capable of retrieving an estimated model consistent with the prior information and fitting the gravity observations according to their accuracy.&lt;/p&gt;

Sign in / Sign up

Export Citation Format

Share Document