SlideShare a Scribd company logo
3
Most read
4
Most read
GRAVITY DATA INTERPRETATION
Group – I
Shimul Al-razi
Muhid Hossain Peal
Sadia Islam
Lovely Khatun
Khalid Hasan
Mahmood Mosaddeque
Sajib Mahmud
INDEX
 The inverse problem
 Regional fields and residual anomalies
 Direct interpretation
 Indirect interpretation
 Summary
THE INVERSE PROBLEM
The interpretation of potential field anomalies (gravity, magnetic and electrical) is
inherently ambiguous. The ambiguity arises because any given anomaly could be
caused by an infinite number of possible sources. For example, concentric spheres of
constant mass but differing density and radius would all produce the same anomaly,
since their mass acts as though located at the center of the sphere. This ambiguity
represents the inverse problem of potential field interpretation, which states that,
although the anomaly of a given body may be calculated uniquely, there are an
infinite number of bodies that could give rise to any specified anomaly. An important
task in interpretation is to decrease this ambiguity by using all available external
constraints on the nature and form of the anomalous body. Such constraints include
geological information derived from surface outcrops, boreholes and mines, and
from other, complementary, geophysical techniques.
An Introduction to Geophysical Exploration. page(139)
Bouguer anomaly fields are often
characterized by a broad, gently varying,
regional anomaly on which may be
superimposed shorter wavelength local
anomalies. Usually in gravity surveying it is
the local anomalies that are of prime interest
and the first step in interpretation is the
removal of the regional field to isolate the
residual anomalies.
REGIONAL FIELDS AND RESIDUAL ANOMALIES
An Introduction to Geophysical Exploration. page(139, 140)
The separation of regional and residual gravity
anomalies from the observed Bouguer anomaly.
DIRECT INTERPRETATION
Direct interpretation provides, directly from the gravity
anomalies, information on the anomalous body which is
largely independent of the true shape of the body.
 Limiting depth
a. Half-width method
b. Gradient–amplitude ratio method
c. Second derivative methods
 Excess mass
a. Inflection point
b. Approximate thickness
An Introduction to Geophysical Exploration. page(140, 141, 142)
a. Half-width method
The depth to a line mass or to the center of a horizontal
cylinder with the same mass distribution is given by
𝑧 =
𝑥1/2
3
4 − 1
For any three-dimensional body, the limiting depth is then
given by
𝑧 <
𝑥1/2
3
4 − 1
.The depth to a line mass or to the center of a horizontal
cylinder with the same mass distribution is given by
𝑧 = 𝑥1/2
For any two-dimensional body, the limiting depth is then given
by
𝑧 < 𝑥1/2
LIMITING DEPTH
Limiting depth calculations using (a) the half-width method and (b) the
gradient–amplitude ratio
An Introduction to Geophysical Exploration. page(140, 141)
(b) Gradient–amplitude ratio method
By differentiation of the relevant
formulae, for any three-dimensional
body
𝑧 < 0.86
𝐴𝑚𝑎𝑥
𝐴𝑚𝑎𝑥
′
and for any two-dimensional body
𝑧 < 0.65
𝐴𝑚𝑎𝑥
𝐴𝑚𝑎𝑥
′
(c) Second derivative methods.
LIMITING DEPTH
Limiting depth calculations using (a) the half-width method and (b) the
gradient–amplitude ratio
An Introduction to Geophysical Exploration. page(140, 141)
EXCESS MASS
The survey area is divided into 𝑛 grid
squares of area ∆𝑎 and the mean residual anomaly ∆𝑔 found for each
square. The excess mass 𝑀𝑒 is then given by
𝑀𝑒 =
1
2𝜋𝐺
𝑖=1
𝑛
∆𝑔𝑖∆𝑎𝑖
To compute the actual mass 𝑀 of the body, the densities of both
anomalous body 𝜌1 and country rock 𝜌2 must be known:
𝑀 =
𝜌1𝑀𝑒
𝜌1 − 𝜌2
An Introduction to Geophysical Exploration. page(141, 142)
(a) Inflection point
The locations of inflection points on gravity profiles,
i.e. positions where the horizontal gravity gradient
changes most rapidly, can provide useful
information on the nature of the edge of an
anomalous body.
(b) Approximate thickness
If the density contrast ∆𝜌 of an anomalous body is
known, its thickness 𝑡 may be crudely estimated
from its maximum gravity anomaly ∆𝑔 by making
use of the Bouguer slab formula
𝑡 ≈
∆𝑔
2𝜋𝐺∆𝜌
EXCESS MASS
Bouguer anomaly profiles across (a) a
granite body, and (b) a sedimentary
basin.
An Introduction to Geophysical Exploration. page(141, 142)
INDIRECT INTERPRETATION
(a) The circular gravity anomaly at Darnley Bay, NWT, Canada. Contour interval 250 gu. (b) Two possible
interpretations of the anomaly in terms of a model constructed from a suite of coaxial vertical cylinders.
(After Stacey 1971.)
An Introduction to Geophysical Exploration. page(142, 143)
The gravity anomaly of this slab ∆𝑔 is given by
∆𝑔 = 2𝐺∆𝜌[− 𝑥1𝑠𝑖𝑛𝜃 + 𝑧1𝑐𝑜𝑠𝜃 × 𝑠𝑖𝑛𝜃𝑙𝑜𝑔𝑒
𝑟2
𝑟1
+ 𝑐𝑜𝑠𝜃 𝜑2 − 𝜑2 + 𝑧2𝜑2 − 𝑧1𝜑1]
INDIRECT INTERPRETATION
Parameters used in defining the gravity anomaly of a
semi-infinite slab with a sloping edge.
The computation of gravity anomalies of two
dimensional bodies of irregular cross-section.
An Introduction to Geophysical Exploration. page(143)
However a model calculation is performed, indirect
interpretation involves four steps:
1. Construction of a reasonable model.
2. Computation of its gravity anomaly.
3. Comparison of computed with observed anomaly.
4. Alteration of model to improve correspondence of
observed and calculated anomalies and return to step 2.
The most flexible of such methods is non-linear optimization
(Al Chalabi 1972). All variables (body points, density contrasts,
regional field) may be allowed to vary within defined limits.
The method then attempts to minimize some function F
which defines the goodness of fit, for example
𝐹 =
𝑖=1
𝑛
∆𝑔𝑜𝑏𝑠𝑖
− ∆𝑔𝑐𝑎𝑙𝑐𝑖
2
INDIRECT INTERPRETATION
A two-dimensional interpretation of the gravity
anomaly of the Bodmin Moor granite, southwest
England. (After Bott & Scott 1964.)
An Introduction to Geophysical Exploration. page(143, 144)
SUMMARY
Gravity data interpretation is useful in exploring regions that have
different geological structures, which contain minerals, ores and oil
deposits. There are different numerical methods for the model
parameters (depth (z), origin location (xo), shape parameter (q) and
amplitude coefficient (A)) evaluation of a covered structure such as
gradient method, particle swarm optimization technique and Werner
deconvolution method. In this study, application of these methods is
utilized to appraise the model parametric quantity of the covered
structures.
Geophysics-Gravity Method

More Related Content

PDF
Magnetic Survey
PPT
Gravity, Expl.ravity
PPTX
Geophysical survey
PPTX
Exploration Geology- Radioactive method of exploration
PPTX
Geophysics. GRAVIMETERY METHOD
PPTX
Kutch
PPTX
Applications Gravity survey Magnetic survey Electrical resistivity survey Sei...
Magnetic Survey
Gravity, Expl.ravity
Geophysical survey
Exploration Geology- Radioactive method of exploration
Geophysics. GRAVIMETERY METHOD
Kutch
Applications Gravity survey Magnetic survey Electrical resistivity survey Sei...

What's hot (20)

PPTX
Gravity Method
PPT
Seismic survey
PPTX
Geophysical methods brief summary
PPTX
Geophysical studies
PPT
Bundelkhand craton
PDF
Seismic Refraction Surveying
PDF
2 d and 3d land seismic data acquisition and seismic data processing
PPTX
Gravity method report
PDF
Introduction to seismic interpretation
PPT
Lecture 13-electrical method -field procedure
PPT
Magnetic prospecting
PPTX
Gravity field separation
PPTX
Shear zones.
PPT
Resistivity method
PPTX
Fmi log presentation
PPTX
Seismic interpretation - Fluvial Deltaic System
PPTX
Geological mapping in Exploration Geology( surface and subsurface)
PPT
Structur Alanalysis
PPT
Anderson s-theory-of-faulting (1)
PPTX
Density log and its uses
Gravity Method
Seismic survey
Geophysical methods brief summary
Geophysical studies
Bundelkhand craton
Seismic Refraction Surveying
2 d and 3d land seismic data acquisition and seismic data processing
Gravity method report
Introduction to seismic interpretation
Lecture 13-electrical method -field procedure
Magnetic prospecting
Gravity field separation
Shear zones.
Resistivity method
Fmi log presentation
Seismic interpretation - Fluvial Deltaic System
Geological mapping in Exploration Geology( surface and subsurface)
Structur Alanalysis
Anderson s-theory-of-faulting (1)
Density log and its uses
Ad

Similar to Geophysics-Gravity Method (20)

PPT
Geotermal 1
PPT
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
PPTX
our igu poster
PDF
EGU2016-2988
PPTX
dsadfsafasdsadasdasdasdsadsadasdasdsadsad
PPT
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
PPTX
minor project ppt (2)
PDF
Gravity Method | Geophysics | Geology
PPT
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
PPTX
GEL 513 Geophysics PowerPoint material for 500 leve
PDF
A Rapid Location Independent Full Tensor Gravity Algorithm
PDF
Crustal Structure from Gravity and Magnetic Anomalies in the Southern Part of...
PDF
Crustal Structure from Gravity and Magnetic Anomalies in the Southern Part of...
PPTX
Grav-nghjghgh jhj jhjjj jhjh jhjhjMag.pptx
PPTX
Grav-Mag-2019-Rajjj kjkjk kjkjjib Sinharay.pptx
PDF
eguPrutkin
PDF
20320130406015 2-3-4
PDF
ambiguity in geophysical interpretation_dileep p allavarapu
PPTX
Gravity anomaly across reagional structures
PDF
Radial Derivative and Radial Inversion for Interpreting 4D Gravity Anomaly Du...
Geotermal 1
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
our igu poster
EGU2016-2988
dsadfsafasdsadasdasdasdsadsadasdasdsadsad
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
minor project ppt (2)
Gravity Method | Geophysics | Geology
ÖNCEL AKADEMİ: INTRODUCTION TO GEOPHYSICS
GEL 513 Geophysics PowerPoint material for 500 leve
A Rapid Location Independent Full Tensor Gravity Algorithm
Crustal Structure from Gravity and Magnetic Anomalies in the Southern Part of...
Crustal Structure from Gravity and Magnetic Anomalies in the Southern Part of...
Grav-nghjghgh jhj jhjjj jhjh jhjhjMag.pptx
Grav-Mag-2019-Rajjj kjkjk kjkjjib Sinharay.pptx
eguPrutkin
20320130406015 2-3-4
ambiguity in geophysical interpretation_dileep p allavarapu
Gravity anomaly across reagional structures
Radial Derivative and Radial Inversion for Interpreting 4D Gravity Anomaly Du...
Ad

Recently uploaded (20)

PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PPTX
OOP with Java - Java Introduction (Basics)
PDF
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
PDF
Automation-in-Manufacturing-Chapter-Introduction.pdf
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PDF
Operating System & Kernel Study Guide-1 - converted.pdf
PDF
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
PDF
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
PDF
Model Code of Practice - Construction Work - 21102022 .pdf
PDF
Embodied AI: Ushering in the Next Era of Intelligent Systems
PPTX
Geodesy 1.pptx...............................................
PPTX
Construction Project Organization Group 2.pptx
PPTX
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PPTX
Sustainable Sites - Green Building Construction
PPTX
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PPTX
Lecture Notes Electrical Wiring System Components
PDF
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
OOP with Java - Java Introduction (Basics)
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
Automation-in-Manufacturing-Chapter-Introduction.pdf
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
Operating System & Kernel Study Guide-1 - converted.pdf
keyrequirementskkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
Model Code of Practice - Construction Work - 21102022 .pdf
Embodied AI: Ushering in the Next Era of Intelligent Systems
Geodesy 1.pptx...............................................
Construction Project Organization Group 2.pptx
Engineering Ethics, Safety and Environment [Autosaved] (1).pptx
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
CYBER-CRIMES AND SECURITY A guide to understanding
Sustainable Sites - Green Building Construction
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
Lecture Notes Electrical Wiring System Components
Mitigating Risks through Effective Management for Enhancing Organizational Pe...

Geophysics-Gravity Method

  • 1. GRAVITY DATA INTERPRETATION Group – I Shimul Al-razi Muhid Hossain Peal Sadia Islam Lovely Khatun Khalid Hasan Mahmood Mosaddeque Sajib Mahmud
  • 2. INDEX  The inverse problem  Regional fields and residual anomalies  Direct interpretation  Indirect interpretation  Summary
  • 3. THE INVERSE PROBLEM The interpretation of potential field anomalies (gravity, magnetic and electrical) is inherently ambiguous. The ambiguity arises because any given anomaly could be caused by an infinite number of possible sources. For example, concentric spheres of constant mass but differing density and radius would all produce the same anomaly, since their mass acts as though located at the center of the sphere. This ambiguity represents the inverse problem of potential field interpretation, which states that, although the anomaly of a given body may be calculated uniquely, there are an infinite number of bodies that could give rise to any specified anomaly. An important task in interpretation is to decrease this ambiguity by using all available external constraints on the nature and form of the anomalous body. Such constraints include geological information derived from surface outcrops, boreholes and mines, and from other, complementary, geophysical techniques. An Introduction to Geophysical Exploration. page(139)
  • 4. Bouguer anomaly fields are often characterized by a broad, gently varying, regional anomaly on which may be superimposed shorter wavelength local anomalies. Usually in gravity surveying it is the local anomalies that are of prime interest and the first step in interpretation is the removal of the regional field to isolate the residual anomalies. REGIONAL FIELDS AND RESIDUAL ANOMALIES An Introduction to Geophysical Exploration. page(139, 140) The separation of regional and residual gravity anomalies from the observed Bouguer anomaly.
  • 5. DIRECT INTERPRETATION Direct interpretation provides, directly from the gravity anomalies, information on the anomalous body which is largely independent of the true shape of the body.  Limiting depth a. Half-width method b. Gradient–amplitude ratio method c. Second derivative methods  Excess mass a. Inflection point b. Approximate thickness An Introduction to Geophysical Exploration. page(140, 141, 142)
  • 6. a. Half-width method The depth to a line mass or to the center of a horizontal cylinder with the same mass distribution is given by 𝑧 = 𝑥1/2 3 4 − 1 For any three-dimensional body, the limiting depth is then given by 𝑧 < 𝑥1/2 3 4 − 1 .The depth to a line mass or to the center of a horizontal cylinder with the same mass distribution is given by 𝑧 = 𝑥1/2 For any two-dimensional body, the limiting depth is then given by 𝑧 < 𝑥1/2 LIMITING DEPTH Limiting depth calculations using (a) the half-width method and (b) the gradient–amplitude ratio An Introduction to Geophysical Exploration. page(140, 141)
  • 7. (b) Gradient–amplitude ratio method By differentiation of the relevant formulae, for any three-dimensional body 𝑧 < 0.86 𝐴𝑚𝑎𝑥 𝐴𝑚𝑎𝑥 ′ and for any two-dimensional body 𝑧 < 0.65 𝐴𝑚𝑎𝑥 𝐴𝑚𝑎𝑥 ′ (c) Second derivative methods. LIMITING DEPTH Limiting depth calculations using (a) the half-width method and (b) the gradient–amplitude ratio An Introduction to Geophysical Exploration. page(140, 141)
  • 8. EXCESS MASS The survey area is divided into 𝑛 grid squares of area ∆𝑎 and the mean residual anomaly ∆𝑔 found for each square. The excess mass 𝑀𝑒 is then given by 𝑀𝑒 = 1 2𝜋𝐺 𝑖=1 𝑛 ∆𝑔𝑖∆𝑎𝑖 To compute the actual mass 𝑀 of the body, the densities of both anomalous body 𝜌1 and country rock 𝜌2 must be known: 𝑀 = 𝜌1𝑀𝑒 𝜌1 − 𝜌2 An Introduction to Geophysical Exploration. page(141, 142)
  • 9. (a) Inflection point The locations of inflection points on gravity profiles, i.e. positions where the horizontal gravity gradient changes most rapidly, can provide useful information on the nature of the edge of an anomalous body. (b) Approximate thickness If the density contrast ∆𝜌 of an anomalous body is known, its thickness 𝑡 may be crudely estimated from its maximum gravity anomaly ∆𝑔 by making use of the Bouguer slab formula 𝑡 ≈ ∆𝑔 2𝜋𝐺∆𝜌 EXCESS MASS Bouguer anomaly profiles across (a) a granite body, and (b) a sedimentary basin. An Introduction to Geophysical Exploration. page(141, 142)
  • 10. INDIRECT INTERPRETATION (a) The circular gravity anomaly at Darnley Bay, NWT, Canada. Contour interval 250 gu. (b) Two possible interpretations of the anomaly in terms of a model constructed from a suite of coaxial vertical cylinders. (After Stacey 1971.) An Introduction to Geophysical Exploration. page(142, 143)
  • 11. The gravity anomaly of this slab ∆𝑔 is given by ∆𝑔 = 2𝐺∆𝜌[− 𝑥1𝑠𝑖𝑛𝜃 + 𝑧1𝑐𝑜𝑠𝜃 × 𝑠𝑖𝑛𝜃𝑙𝑜𝑔𝑒 𝑟2 𝑟1 + 𝑐𝑜𝑠𝜃 𝜑2 − 𝜑2 + 𝑧2𝜑2 − 𝑧1𝜑1] INDIRECT INTERPRETATION Parameters used in defining the gravity anomaly of a semi-infinite slab with a sloping edge. The computation of gravity anomalies of two dimensional bodies of irregular cross-section. An Introduction to Geophysical Exploration. page(143)
  • 12. However a model calculation is performed, indirect interpretation involves four steps: 1. Construction of a reasonable model. 2. Computation of its gravity anomaly. 3. Comparison of computed with observed anomaly. 4. Alteration of model to improve correspondence of observed and calculated anomalies and return to step 2. The most flexible of such methods is non-linear optimization (Al Chalabi 1972). All variables (body points, density contrasts, regional field) may be allowed to vary within defined limits. The method then attempts to minimize some function F which defines the goodness of fit, for example 𝐹 = 𝑖=1 𝑛 ∆𝑔𝑜𝑏𝑠𝑖 − ∆𝑔𝑐𝑎𝑙𝑐𝑖 2 INDIRECT INTERPRETATION A two-dimensional interpretation of the gravity anomaly of the Bodmin Moor granite, southwest England. (After Bott & Scott 1964.) An Introduction to Geophysical Exploration. page(143, 144)
  • 13. SUMMARY Gravity data interpretation is useful in exploring regions that have different geological structures, which contain minerals, ores and oil deposits. There are different numerical methods for the model parameters (depth (z), origin location (xo), shape parameter (q) and amplitude coefficient (A)) evaluation of a covered structure such as gradient method, particle swarm optimization technique and Werner deconvolution method. In this study, application of these methods is utilized to appraise the model parametric quantity of the covered structures.