Inversion of the power spectrum from gravity anomalies of prismatic bodies

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1698-1703 ◽  
Author(s):  
Juan García‐Abdeslem
Keyword(s):  
Inverse Problem ◽  
Power Spectrum ◽  
Ridge Regression ◽  
Field Data ◽  
Gravity Data ◽  
Source Parameters ◽  
Gravity Anomalies ◽  
Independent Functions ◽  
Prismatic Bodies ◽  
And Inversion

I develop methods for modeling and inversion of gravity anomalies produced by prismatic bodies. The forward problem is solved in the wavenumber domain, where the power spectrum of the gravity anomaly is given by the product of independent functions that describe depth, thickness, horizontal dimensions, and density of the source body. The solution of the inverse problem is based on an iterative ridge‐regression algorithm, starting from an initial trial of the geometry and density of the source body. This procedure is assessed in a number of examples with both numeric and field data spectra. The method is first tested on the spectrum obtained from simple gravity anomalies and is found to be capable of recovering the source parameters. Next the method is applied to gravity data from a survey made near Noranda, Quebec, Canada. This interpretation compares favorably with drillhole data and provides an estimation of the mass of the source body which is similar to previous estimates.

Linear inversion of gravity data using the spectral expansion method

Geophysics ◽  
1985 ◽  
Vol 50 (5) ◽  
pp. 820-824 ◽  
Author(s):  
Rene E. Chavez ◽  
George D. Garland
Keyword(s):  
Inverse Problem ◽  
Gravity Data ◽  
Linear Equations ◽  
Thin Sheet ◽  
Gravity Anomalies ◽  
Expansion Method ◽  
Spectral Expansion ◽  
Uniform Density ◽  
Linear Inversion ◽  
Prismatic Structure

Inversion of gravity anomalies in terms of an anomalous mass distribution with irregular outline but uniform density leads to a nonlinear inverse problem. An alternative approach based on the thin‐sheet approximation can, however, be formulated as a linear inverse problem provided the structures are two‐dimensional. The anomalous mass is represented by a thin sheet, which is located at a depth [Formula: see text] and is divided into M strips with N < M data points. Thus, an underdetermined system of linear equations is obtained which is solved by the spectral expansion method for the surface density distribution of each segment. This set of parameters is then transformed into a prismatic structure with variable depth but uniform density. The modeling procedure involves a noniterative method. A gravity problem is investigated, and the solution obtained compares well with previous interpretations.

A method for inversion of 1D vertical soundings of gravity anomalies

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. G15-G23
Author(s):  
Andrea Vitale ◽  
Domenico Di Massa ◽  
Maurizio Fedi ◽  
Giovanni Florio
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Gravity Data ◽  
Finite Size ◽  
Gravity Anomalies ◽  
Real Data ◽  
Gravity Inversion ◽  
Chain Process ◽  
Potential Field Data ◽  
Gamma Density

We have developed a method to interpret potential fields, which obtains 1D models by inverting vertical soundings of potential field data. The vertical soundings are built through upward continuation of potential field data, measured on either a profile or a surface. The method assumes a forward problem consisting of a volume partitioned in layers, each of them homogeneous and horizontally finite, but with the density changing versus depth. The continuation errors, increasing with the altitude, are automatically handled by determining the coefficients of a third-order polynomial function of the altitude. Due to the finite size of the source volume, we need a priori information about the total horizontal extent of the volume, which is estimated by boundary analysis and optimized by a Markov chain process. For each sounding, a 1D inverse problem is independently solved by a nonnegative least-squares algorithm. Merging of the several inverted models finally yields approximate 2D or 3D models that are, however, shown to generate a good fit to the measured data. The method is applied to synthetic models, producing good results for either perfect or continued data. Even for real data, i.e., the gravity data of a sedimentary basin in Nevada, the results are interesting, and they are consistent with previous interpretation, based on 3D gravity inversion constrained by two gamma-gamma density logs.

scholarly journals Interpretação gravimétrica na Borda Norte do Cráton do Congo, Sul de Camarões

2010 ◽  
Vol 33 (1) ◽  
pp. 73-82
Author(s):  
Yves N Shandini ◽  
Jean Marie Tadjou ◽  
Charles T Tabod ◽  
James Derek Fairhead
Keyword(s):  
Gravity Data ◽  
Bouguer Anomaly ◽  
Gravity Anomalies ◽  
Continental Collision ◽  
Tectonic Structures ◽  
Collision Zone ◽  
Granitic Intrusion ◽  
Southern Cameroon ◽  
Northern Edge ◽  
And Inversion

Gravity data in the southern Cameroon are interpreted to better understand the organization of underlying structures throughout the northern edge of the Congo craton. The Bouguer anomaly maps of the region are characterized by an elongated SW-NE trending negative gravity anomaly which correspond to a collapsed structure associated with a granitic intrusion beneath the center of the region and limited by fault systems. We applied 3-D gravity modelling and inversion in order to obtain the 3-D density structure of the area. Our result demonstrated that observed gravity anomalies in the region are associated to tectonic structures in the subsurface. The resulting model agrees with the hypothesis of the existence of a major continental collision zone between the Congo Craton and the Pan-African belt. The presence of deep granulites structures in the northern part of the region expresses a continental collision.

Nonlinear inversion of isostatic residual gravity data from Montage Basin, northern Gulf of California

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. G45-G55 ◽  
Author(s):  
Juan García-Abdeslem
Keyword(s):  
Inverse Problem ◽  
Gravity Anomaly ◽  
Gulf Of California ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Nonlinear Inversion ◽  
Residual Gravity ◽  
Sedimentary Sequences ◽  
Residual Gravity Anomaly ◽  
Northwestern Mexico

The flexural isostatic response to surface loads is used to estimate the crustal thickness in northwestern Mexico and Southwestern USA. This estimate is used to compute an isostatic regional gravity, which, subtracted from Bouguer gravity anomalies, led to the isostatic residual gravity anomaly at Montage Basin. This basin is located between the southern portion of the Mexicali Valley and the northern Gulf of California, it roughly has an extension of [Formula: see text] wide, and it shows a gravity minimum reaching approximately [Formula: see text]. Montage Basin is within the extensional province of the Gulf of California, where rifting is currently an ongoing geologic process, and deep exploratory wells drilled by Petróleos Mexicanos have shown that the basin accommodates thick sedimentary sequences greater than 5 km. The interpretation of the isostatic residual gravity anomaly is considered as a nonlinear inverse problem, constrained using density as a function of depth derived from Gardner’s equation applied to dual time [Formula: see text]-logs, assuming isostatic equilibrium and considering the basin as a subsurface load that is compensated at depth by a mass of unknown shape and density. The outcome of the inverse problem suggests that Montage Basin accommodates as much as 7.5 km thick sedimentary sequences and a compensating mass at a minimum depth of 13 km.

scholarly journals Joint Gravity and Seismic Reflection Methods to Characterize the Deep Aquifers in Arid Ain El Beidha Plain (Central Tunisia, North Africa)

Water ◽  
2021 ◽  
Vol 13 (9) ◽  
pp. 1310
Author(s):  
Hajer Azaiez ◽  
Hakim Gabtni ◽  
Mourad Bédir
Keyword(s):  
Seismic Reflection ◽  
Large Scale ◽  
Seismic Analysis ◽  
Gravity Data ◽  
Geophysical Methods ◽  
Gravity Anomalies ◽  
Integrated Approach ◽  
Central Tunisia ◽  
Deep Aquifers ◽  
Advanced Analysis

Electric resistivity sounding and tomography, as well as electromagnetic sounding, are the classical methods frequently used for hydrogeological studies. In this work, we propose the development and implementation of an original integrated approach using the unconventional hydro–geophysical methods of gravity and seismic reflection for the fast, large–scale characterization of hydrogeological potential using the Ain El Beidha plain (central Tunisia) as an analogue. Extending the values of vintage petroleum seismic reflection profiles and gravity data, in conjunction with available geological and hydrogeological information, we performed an advanced analysis to characterize the geometry of deep tertiary (Oligocene and Eocene) aquifers in this arid area. Residual and tilt angle gravity maps revealed that most gravity anomalies have a short wavelength. The study area was mainly composed of three major areas: the Oued Ben Zitoun and Ain El Beidha basins, which are both related to negative gravity trends corresponding to low–density subsiding depocenters. These basins are separated by an important NE–SW trend called “El Gonna–J. El Mguataa–Kroumet Zemla” gravity high. Evaluation of the superposition of detected lineaments and Euler deconvolution solutions’ maps showed several NE–SW and N–S relay system faults. The 3D density inversion model using a lateral and vertical cutting plane suggested the presence of two different tectonic styles (thin VS thick). Results from the gravity analysis were in concordance with the seismic analysis. The deep Oligocene and Eocene seismic horizons were calibrated to the hydraulic wells and surrounding outcrops. Oligocene and Eocene geological reservoirs appear very fractured and compartmented. The faulting network also plays an important role in enhancing groundwater recharge process of the Oligocene and Eocene aquifers. Finally, generated isochron maps provided an excellent opportunity to develop future comprehensive exploration surveys over smaller and more favorable areas’ sub–basins.

Computer Application in Geosciences: Imaging of the Subsurface by Structural Coupled Inversion of Potential Field Data

2014 ◽  
Vol 644-650 ◽  
pp. 2670-2673
Author(s):  
Jun Wang ◽  
Xiao Hong Meng ◽  
Fang Li ◽  
Jun Jie Zhou
Keyword(s):  
Potential Field ◽  
Field Data ◽  
Gravity Data ◽  
Computer Application ◽  
Magnetic Data ◽  
Imaging Method ◽  
Inversion Algorithm ◽  
Constant Property ◽  
Near Surface ◽  
Potential Field Data

With the continuing growth in influence of near surface geophysics, the research of the subsurface structure is of great significance. Geophysical imaging is one of the efficient computer tools that can be applied. This paper utilize the inversion of potential field data to do the subsurface imaging. Here, gravity data and magnetic data are inverted together with structural coupled inversion algorithm. The subspace (model space) is divided into a set of rectangular cells by an orthogonal 2D mesh and assume a constant property (density and magnetic susceptibility) value within each cell. The inversion matrix equation is solved as an unconstrained optimization problem with conjugate gradient method (CG). This imaging method is applied to synthetic data for typical models of gravity and magnetic anomalies and is tested on field data.

scholarly journals Basin Depth Mapping Beneath Wellington City Based on Residual Gravity Anomalies

2021 ◽  
Author(s):  
◽  
Alistair Stronach
Keyword(s):  
Gravity Anomaly ◽  
Sedimentary Basin ◽  
Gravity Data ◽  
Three Dimensional ◽  
Maximum Amplitude ◽  
Central Business District ◽  
Gravity Anomalies ◽  
Depth Map ◽  
Capital City ◽  
Detailed Knowledge

<p><b>New Zealand’s capital city of Wellington lies in an area of high seismic risk, which is further increased by the sedimentary basin beneath the Central Business District (CBD). Ground motion data and damage patterns from the 2013 Cook Strait and 2016 Kaikōura earthquakes indicate that two- and three-dimensional amplification effects due to the Wellington sedimentary basin may be significant. These effects are not currently accounted for in the New Zealand Building Code. In order for this to be done, three-dimensional simulations of earthquake shaking need to be undertaken, which requires detailed knowledge of basin geometry. This is currently lacking, primarily because of a dearth of deep boreholes in the CBD area, particularly in Thorndon and Pipitea where sediment depths are estimated to be greatest.</b></p> <p>A new basin depth map for the Wellington CBD has been created by conducting a gravity survey using a modern Scintrex CG-6 gravity meter. Across the study area, 519 new high precision gravity measurements were made and a residual anomaly map created, showing a maximum amplitude anomaly of -6.2 mGal with uncertainties better than ±0.1 mGal. Thirteen two-dimensional geological profiles were modelled to fit the anomalies, then combined with existing borehole constraints to construct the basin depth map. </p> <p>Results indicate on average greater depths than in existing models, particularly in Pipitea where depths are interpreted to be as great as 450 m, a difference of 250 m. Within 1 km of shore depths are interpreted to increase further, to 600 m. The recently discovered basin bounding Aotea Fault is resolved in the gravity data, where the basement is offset by up to 13 m, gravity anomaly gradients up to 8 mGal/km are observed, and possible multiple fault strands identified. A secondary strand of the Wellington Fault is also identified in the north of Pipitea, where gravity anomaly gradients up to 18 mGal/km are observed.</p>

scholarly journals Basement Mapping of the Fucino Basin in Central Italy by ITRESC Modeling of Gravity Data

Geosciences ◽  
2021 ◽  
Vol 11 (10) ◽  
pp. 398
Author(s):  
Federico Cella ◽  
Rosa Nappi ◽  
Valeria Paoletti ◽  
Giovanni Florio
Keyword(s):  
Seismic Activity ◽  
Site Response ◽  
Gravity Data ◽  
A Priori ◽  
Central Italy ◽  
Gravity Anomalies ◽  
Gravity Modeling ◽  
Ground Shaking ◽  
Sedimentary Evolution ◽  
Fucino Basin

Sediments infilling in intermontane basins in areas with high seismic activity can strongly affect ground-shaking phenomena at the surface. Estimates of thickness and density distribution within these basin infills are crucial for ground motion amplification analysis, especially where demographic growth in human settlements has implied increasing seismic risk. We employed a 3D gravity modeling technique (ITerative RESCaling—ITRESC) to investigate the Fucino Basin (Apennines, central Italy), a half-graben basin in which intense seismic activity has recently occurred. For the first time in this region, a 3D model of the Meso-Cenozoic carbonate basement morphology was retrieved through the inversion of gravity data. Taking advantage of the ITRESC technique, (1) we were able to (1) perform an integration of geophysical and geological data constraints and (2) determine a density contrast function through a data-driven process. Thus, we avoided assuming a priori information. Finally, we provided a model that honored the gravity anomalies field by integrating many different kinds of depth constraints. Our results confirmed evidence from previous studies concerning the overall shape of the basin; however, we also highlighted several local discrepancies, such as: (a) the position of several fault lines, (b) the position of the main depocenter, and (c) the isopach map. We also pointed out the existence of a new, unknown fault, and of new features concerning known faults. All of these elements provided useful contributions to the study of the tectono-sedimentary evolution of the basin, as well as key information for assessing the local site-response effects, in terms of seismic hazards.

scholarly journals DETERMINATION AND MAPPING OF THE TERRESTRIAL GRAVITY ANOMALIES IN THE MOUNTAINOUS AREAS OF IRAQ

2021 ◽  
Vol 6 (24) ◽  
pp. 213-225
Author(s):  
Shazad Jamal Jalal ◽  
Tajul Ariffin Musa ◽  
Ami Hassan Md Din ◽  
Wan Anom Wan Aris
Keyword(s):  
Gravity Data ◽  
Physical Geodesy ◽  
Gravity Anomalies ◽  
Gravity Survey ◽  
Mountainous Areas ◽  
Entire Study ◽  
Terrestrial Gravity ◽  
The Past ◽  
Regional Gravity ◽  
Entire Study Area

Gravity data and computing gravity anomalies are regarded as vital for both geophysics and physical geodesy fields. The mountainous areas of Iraq are characterized by the lack of regional gravity data because gravity surveys are rarely performed in the past four decades due to the Iraq-Iran war and the internal unstable political situation of this particular region. In addition, the formal map of the available terrestrial gravity which was published by the French Database of Bureau Gravimetrique International (International Gravimetric Bureau-in English) (BGI), introduces Iraq and the study area as a remote area and in white color because of the unavailability of gravity data. However, a dense and local (not regional) gravity data is available which was conducted by geophysics researchers 13 years ago. Therefore, the regional gravity survey of 160 gravity points was performed by the authors at an average 11 km apart, which was covers the whole area of Sulaymaniyah Governorate (part of the mountainous areas of Iraq). In spite of Although the risk of mine fields within the study area, suitable safe routes as well as a helicopter was used for the gravity survey of several points on the top of mountains. The survey was conducted via Lacoste and Romberg geodetic gravimeter and GPS handheld. The objective of the study is to determine and map the gravity anomalies for the entire study area, the data of which would assist different geosciences applications.

scholarly journals Structural model of the Northern Latium volcanic area constrained by MT, gravity and aeromagnetic data

1997 ◽  
Vol 40 (5) ◽  
Author(s):  
P. Capuano ◽  
G. Florio ◽  
P. Gasparini
Keyword(s):  
Structural Model ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Volcanic Area ◽  
Aeromagnetic Data ◽  
Geothermal Exploration ◽  
Magnetotelluric Soundings ◽  
Aeromagnetic Anomalies ◽  
Deep Boreholes

The results of about 120 magnetotelluric soundings carried out in the Vulsini, Vico and Sabatini volcanic areas were modeled along with Bouguer and aeromagnetic anomalies to reconstruct a model of the structure of the shallow (less than 5 km of depth) crust. The interpretations were constrained by the information gathered from the deep boreholes drilled for geothermal exploration. MT and aeromagnetic anomalies allow the depth to the top of the sedimentary basement and the thickness of the volcanic layer to be inferred. Gravity anomalies are strongly affected by the variations of morphology of the top of the sedimentary basement, consisting of a Tertiary flysch, and of the interface with the underlying Mesozoic carbonates. Gravity data have also been used to extrapolate the thickness of the neogenic unit indicated by some boreholes. There is no evidence for other important density and susceptibility heterogeneities and deeper sources of magnetic and/or gravity anomalies in all the surveyed area.

Sign in / Sign up

Export Citation Format

Share Document