0% found this document useful (0 votes)
153 views

L03 GravCorrAnalysis PDF

The document discusses various corrections that must be applied to gravity measurements, including corrections for latitude, elevation, topography, and density, in order to isolate the gravitational effect of subsurface structures. It also covers analyzing gravity anomalies by calculating theoretical anomalies of simple shapes and distributions and comparing to observed anomalies to interpret subsurface structures.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
153 views

L03 GravCorrAnalysis PDF

The document discusses various corrections that must be applied to gravity measurements, including corrections for latitude, elevation, topography, and density, in order to isolate the gravitational effect of subsurface structures. It also covers analyzing gravity anomalies by calculating theoretical anomalies of simple shapes and distributions and comparing to observed anomalies to interpret subsurface structures.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 11

Gravity:

Corrections and analysis

Reading:
Today: p22-39
Next Lecture: p39-64

Applied Geophysics – Corrections and analysis

Gravity corrections

Measure gravity variations at


stations around loop
Correct for drift

Observations still subject to extraneous effects unrelated to


subsurface geology
Î Must make corrections…
1. Latitude correction
2. Free-air correction
3. Bouguer correction
4. Terrain correction

Applied Geophysics – Corrections and analysis

1
Latitude correction
…correct for the spheroid

Geodetic Reference System (GRS-1967) formula


gφ = 9.780318(1 + 0.0053024 sin 2 φ − 0.0000059 sin 2 2φ ) m/s2

For smaller scale studies we can simplify:

Cφ = 8.12 sin 2φ g.u./km-N-S

Correction is greatest at mid-latitudes.


For 0.1 g.u. accuracy need relative N-S
distance to 12 m
Add or subtract?

Applied Geophysics – Corrections and analysis

Free-air correction
Accounts for the 1/r2 decrease in gravity with GM E
distance from the center of the Earth, recall: g= 2
RE
We calculate:
C F = 3.086h
where h is the elevation in meters

For 0.1 g.u. accuracy stations elevation is


needed to 4 cm.

Add or subtract?

Applied Geophysics – Corrections and analysis

2
Bouguer correction
Accounts for rock thickness between current and base station elevation

∆h

Treat the rock as an infinite horizontal slab:

C B = 0.000419∆hρ where ∆h is in m and ρ is in km/m3

For 0.1 g.u. accuracy station elevation is needed to 9 cm


Add or subtract?

Applied Geophysics – Corrections and analysis

Terrain correction
When Bouguer correction is inadequate, also use terrain correction

Add or subtract?
hill

valley

Approaches to correction
Rectangular grid
Hammer segments … both use elevation differences
between station and surround

Applied Geophysics – Corrections and analysis

3
Terrain correction
Hammer correction

Applied Geophysics – Corrections and analysis

Free-air anomaly
A gravity “anomaly” suggests the difference between a
theoretical and observed value.
Î Tie survey to base station where absolute gravity is known

Then the free-air anomaly is:


∆g F = g obs + C F − gφ

Have not corrected for topography


• Onshore: anomaly map similar to topography
• Free-air anomaly mainly used for offshore

Applied Geophysics – Corrections and analysis

4
Bouguer anomaly

Apply all the corrections:

∆g B = ∆g obs + Cφ + C F − C B + CT
…watch the signs!

Smaller scale engineering/environmental surveys:


• Not tied to absolute gravity
• Use corrections with accuracy necessary
…determined by the size of the target signal

Applied Geophysics – Corrections and analysis

Field determination of density


Why do we need rock density?

Nettleton method:
density profile

Other methods:
• The same problem can
be formulated
mathematically and
solved by least squares
• Borehole logging

Applied Geophysics – Corrections and analysis

5
Analysis and interpretation
Once we have made our gravity observations,
corrected for surface effects,
we attempt to deduce sub-surface structure

Considerations:
• Anomaly profile (2D structure) or map (3D structure)?
If anomaly length > twice the width a 2D interpretation is OK
• Ambiguity
There are an infinite number of structures that could generate the
observations
• Forward calculation
“Guess” at structure, calculate the anomaly and compare. Simple
formula for depth and size of geometric shapes.
• Inverse modeling
Directly invert for structure, choose constraints on the geometry
…don’t forget that ambiguity.
Applied Geophysics – Corrections and analysis

Buried sphere
Analytic expressions for simple
geometric shapes
e.g. a buried sphere

4πGR 3 ∆ρ 1
∆g z =
3z 2 [ (
1+ x z2
2
)]
32

Depth rule
z = 1.302 x1 2

Note: it is only the density


contrast that is important

Applied Geophysics – Corrections and analysis

6
Gravity anomaly map

Applied Geophysics – Corrections and analysis

Simple shape anomalies

Applied Geophysics – Corrections and analysis

7
2D vertical column
Gravity anomaly due to a
two dimensional vertical
column:

∆g z = 2G∆ρb ln r2 
 r1 

Applied Geophysics – Corrections and analysis

2D vertical columns
Gravity anomaly due to a
two dimensional vertical
column:

∆g z = 2G∆ρb ln r2 
 r1 
Gravity anomaly due to a
series of vertical columns:

r2 
∆g z = 2G ∑ ∆ρ i bi ln i 1 
i  ri 

Applied Geophysics – Corrections and analysis

8
Gravity Anomalies: 2D forward calculation
for rectangular parallelepipeds
Spreadsheet: with greater vertical extent than horizontal
Grav2Dcolumn
see Dobrin and Savit eq 12-34

Define density structure Calculated gravity anom aly


4.50
Adjust bold numbers…
4.00
coulum density
center contrast top bottom error 3.50
(km) (g/cm3) (km) (km) check 3.00

dgz (mGal)
2.50
0 0.5 0 0 OK
2.00
1 0.5 0 0 OK
2 0.5 0 0 OK 1.50
3 0.5 0 0 OK 1.00
4 0.5 2 8 OK 0.50
5 0.5 0 0 OK
0.00
6 0.5 0 0 OK
0 2 4 6 8 10 12 14 16 18 20
7 0.5 0 0 OK
distance (km )
8 0.5 0 0 OK
9 0.5 0 0 OK
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
10 0.5 0 0 OK
0
11 0.5 0 0 OK
12 0.5 0 0 OK 1
13 0.5 0 0 OK 2
14 0.5 0 0 OK 3
depth (km)

15 0.5 0 0 OK
4
16 0.5 0 0 OK
17 1 0 0 OK 5
18 0.5 0 0 OK 6
19 0.5 0 0 OK
7
20 0.5 0 0 OK
8
9

Applied Geophysics – Corrections and analysis

Gravity Anomalies: 2D forward calculation


for rectangular parallelepipeds
Spreadsheet: with greater vertical extent than horizontal
Grav2Dcolumn
see Dobrin and Savit eq 12-34

Define density structure Calculated gravity anom aly


10.00
Adjust bold numbers… 9.00
coulum density
8.00
center contrast top bottom error
7.00
(km) (g/cm3) (km) (km) check
dgz (mGal)

6.00
0 0.5 0 0 OK 5.00
1 0.5 0 0 OK 4.00
2 0.5 7 9 OK 3.00
3 0.5 6 10 OK 2.00
4 0.5 5.5 9.5 OK
1.00
5 0.5 5 9 OK
0.00
6 0.5 4.7 8 OK
0 2 4 6 8 10 12 14 16 18 20
7 0.5 4.5 7 OK
distance (km )
8 0.5 4.4 6 OK
9 0.5 4.3 5.5 OK
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
10 0.5 0 0 OK
0
11 0.5 0 0 OK
12 0.5 0 0 OK 2
13 0.5 0 0 OK
14 0.5 0 0 OK 4
depth (km)

15 0.5 0 0 OK
16 0.5 0 0 OK 6
17 1 1 2 OK
18 0.5 0 0 OK 8
19 0.5 0 0 OK
20 0.5 0 0 OK 10

12

Applied Geophysics – Corrections and analysis

9
Ambiguity - I

Applied Geophysics – Corrections and analysis

Ambiguity - II

4πGR 3 ∆ρ 1
∆g z =
3z 2
[ (
1+ x z2
2
)]
32

Applied Geophysics – Corrections and analysis

10
Gravity Anomalies: 2D forward calculation
for rectangular parallelepipeds with greater vertical extent than horizontal
Gravity anom aly
Spreadsheet:
see Dobrin and Grav2Dcolumn
Savit eq 12-34
2.50
Define density structure Profile 1
2.00 Profile 2

Profile 1 Profile2

dgz (mGal)
1.50
Adjust bold numbers… Adjust bold numbers…
coulum density density
1.00
center contrast top bottom error contrast bottom error
(km) (g/cm3) (km) (km) check (g/cm3) top (km) (km) check
0.50

0 0.3 0 0 OK 0.9 0 0 OK
0.00
1 0.3 4 4.4 OK 0.9 0 0 OK
0 2 4 6 8 10 12 14 16 18 20
2 0.3 4 4.4 OK 0.9 0 0 OK distance (km)
3 0.3 4 4.3 OK 0.9 0 0 OK
4 0.3 4 4.3 OK 0.9 0 0 OK 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
5 0.3 4 4.4 OK 0.9 0 0 OK 0
6 0.3 4 4.4 OK 0.9 0 0 OK
2
7 0.3 4 4.5 OK 0.9 0 0 OK
8 0.3 4 4.6 OK 0.9 0 0 OK 4

depth (km)
9 0.3 4 4.7 OK 0.9 0 0 OK 6
10 0.3 4 4.8 OK 0.9 8 12 OK
8
11 0.3 4 4.7 OK 0.9 0 0 OK
12 0.3 4 4.6 OK 0.9 0 0 OK 10
13 0.3 4 4.5 OK 0.9 0 0 OK 12 Profile 1
14 0.3 4 4.4 OK 0.9 0 0 OK
14
15 0.3 4 4.4 OK 0.9 0 0 OK
16 0.3 4 4.3 OK 0.9 0 0 OK 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
0
17 0.3 4 4.3 OK 0.9 0 0 OK
18 0.3 4 4.4 OK 0.9 0 0 OK 2

19 0.3 4 4.4 OK 0.9 0 0 OK 4


20 0.3 0 0 OK 0.9 0 0 OK
6

10

12 Profile 2
14

Applied Geophysics – Corrections and analysis

11

You might also like