6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.

htm

6. Gravimetría
1. Introducción

El método gravimétrico hace uso de campos de potencial natural igual al método magnético y a algunos
métodos eléctricos. El campo de potencial natural observado se compone de los contribuyentes de las
formaciones geológicas, que construyen la corteza terrestre hasta cierta profundidad determinada por el
alcance del método gravimétrico (o magnético respectivamente). Generalmente no se puede distinguir las
contribuciones a este campo proveniente de una formación o una estructura geológica de aquellas de las
otras formaciones o estructuras geológicas por el método gravimétrico, solo en casos especiales se puede
lograr una separación de los efectos causados por una formación o estructura geológica individual. Se
realiza mediciones relativas o es decir se mide las variaciones laterales de la atracción gravitatoria de un
lugar al otro puesto que en estas mediciones se pueden lograr una precisión satisfactoria más fácilmente
en comparación con las mediciones del campo gravitatorio absoluto. Los datos reducidos apropiadamente
entregan las variaciones en la gravedad, que solo dependen de variaciones laterales en la densidad del
material ubicado en la vecindad de la estación de observación.

2. Historia

El método gravimétrico fue aplicado inicialmente en la prospección petrolífera en los Estados Unidos y en
el golfo de México con el objetivo de localizar domos de sales, que potencialmente albergan petroleo.
Luego se buscaron estructuras anticlinales con este método. El fin del siglo 19 el húngaro Roland von
EÖTVÖS desarrolló la balanza de torsión llamada según él, que mide las distorsiones del campo
gravitatorio causadas de cuerpos de densidades anómalas enterrados en el subsuelo como de domos de
sal o cuerpos de cromita por ejemplo. En 1915 y 1916 se emplearon la balanza de torsión de EÖTVÖS en
el levantamiento de la estructura de un campo petrolífero ubicado en Egbell en la Checoslovaquia antigua.
En 1917 SCHWEIDAR levantó un domo de sal ya conocido ubicado cerca de Hanigsen en Alemania por
medio de una balanza de torsión y la estructura deducida y predicha a partir de esos estudios fue
confirmada luego por sondeos.

3. Principio

3.1 Ley de gravitación de NEWTON

Si cualquier cuerpo inicialmente estando en reposo cae sin ser estorbado después un segundo tendrá una
velocidad de 9,80m/s en la dirección vertical. Después de un segundo más su velocidad será: 9,80m/s +
9,80m/s = 19,60m/s. El aumento de la velocidad vertical de 9,80m/s de un cuerpo cayendo sin ser
estorbado durante cada segundo se denomina aceleración de gravedad o sólo gravedad y se la expresa
como 9,80m por segundo por segundo o es decir 9,80m/s 2. El primero término por segundo indica la
velocidad medida como distancia pasada durante un segundo, el otro por segundo indica la variación de la
velocidad de 9,80m/s, que corresponde a un intervalo de 1s. La aceleración de la gravedad g se debe a la
aceleración gravitatoria, que la tierra ejerce en cada cuerpo, menos la fuerza centrífuga causada por la
rotación de la tierra y dirigida en dirección perpendicular al eje de rotación de la tierra y hacia fuera. La
fuerza total, que actúa en el cuerpo, es igual al producto de su masa m y de la aceleración de gravedad g.
Por consiguiente la atracción gravitatoria en cualquier lugar de la superficie terrestre tiene numéricamente
el mismo valor como la fuerza gravitatoria ejercida a una masa unitaria en el mismo lugar.

Según la ley de gravitación de NEWTON los cuerpos de las masas m 1 y m2 separados por una distancia r
se influyen mutuamente por la fuerza F:

F = f ×((m1 × m2)/r2),

donde m1, m2 = masa del cuerpo 1 o 2 respectivamente,

1 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

r = distancia entre los centros de los cuerpos de masa m 1 y m2.

f = constante de gravitación = 6,67 × 10-8cm3g-1s-2 = 6,67 × 10-11Nm2/kg2 (N = kgm/s2). La constante

de gravitación f describe la fuerza expresada en N (Newton) ejercida entre dos cuerpos de masas 1kg,
cuyos centros distan 1m entre sí y cuyas masas están concentradas en sus centros. Se la mide en el
laboratorio. En el año 1797 la primera vez CAVENDISH realizó una medición de f resultando en un valor de
f = 6,754 × 10-11Nm2/kg2.

F = m1 × a,

donde m1 = masa del cuerpo 1 en consideración

a = aceleración producida por la masa m1 en su vecindad.

La aceleración debida a un cuerpo de masa m1 en un punto de masa m2 en distancia r con respecto al

centre del cuerpo de masa m1 se obtiene por división de la ecuación 'F = m1 × a = f × (m1 × m2)/r2' con
m2. Por consiguiente: a = f × (m1/r2).

La unidad de la aceleración a es 1cm/s2 = 1 Gal (según Galilei) y 0,001cm/s2 = 1mgal = 10gu (unidades
de gravedad).

La unidad de la variación de la aceleración o es decir del gradiente de la aceleración es 1s -2, 10-8s-2 =

1mgal/km y 10-9s-2 = 1E (Eötvös).

3.2 El potencial y el campo gravitatorio de la Tierra

El potencial en un punto de un campo dado se define como el trabajo rendido por la fuerza al mover una
masa unitaria desde un punto arbitrario de referencia - usualmente ubicándose en una distancia infinita -
hacia el punto en cuestión.

El potencial correspondiente al cuerpo de la masa m1 se calcula: P = -f × m1/r.

La diferencia en los potenciales P2 - P1 describe el trabajo rendido en contra de la masa m1 al mover una
masa unitaria desde el centro del cuerpo m1 al centro del cuerpo m2.

Las superficies equipotenciales (superficies, que unen todos los puntos del mismo valor potencial) referidas
a este cuerpo de masa m1 son superficies esféricas. El potencial correspondiente al espacio exterior de una
esfera de estructura de estratos es igual al potencial correspondiente al punto material central, en que
está concentrado la masa total de esta esfera. Este hecho se aplica para describir y cuantificar el campo
potencial gravitatorio de la Tierra.

Dos fuerzas distintas contribuyen al campo gravitatorio de la Tierra. En un lugar de la superficie terrestre
la fuerza gravitatoria neta GN ejercida se constituye de la fuerza gravitatoria dirigida hacia el centro de la
Tierra GT y la fuerza centrífuga GC dirigida perpendicularmente al eje rotativo y afuera referente a la
Tierra. Por consiguiente GN = GT + GC. La fuerza centrífuga se calcula de la manera siguiente:

GC = mT × aC = mT × V2 × rT × senb,

donde b = 90º-ß , ß = latitud geográfica,

V = velocidad angular de la rotación de la Tierra = 7,29 × 10 5s-1, rT = radio de la Tierra,

mT = masa de la Tierra.

2 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

Salvo a los polos, donde aC = 0 debido a b = 0º, la fuerza centrífuga actúa en todos los demás lugares de
la superficie terrestre e es apreciadamente menor en comparación a G T. Por esto se abrevia la fuerza
gravitatoria neta solo con 'g'. En la medición de la fuerza gravitatoria neta no se puede distinguir entre G T
y G C.

La aceleración gravitatoria presente en una dirección definida se obtiene por diferenciación del potencial
con respecto a la distancia en esta dirección. La superficie caracterizada por valores del potencial
constantes se denomina superficie equipotencial. A lo largo de una superficie equipotencial se puede
mover un cuerpo de un lugar al otro sin esforzarse en o en dirección opuesta a la gravedad. Una superficie
equipotencial es la superficie del mar, aun la fuerza gravitatoria varía a lo largo de esta superficie mas que
0,5% entre el ecuador y los polos.

3.3 La forma teórica y la forma geométrica de la Tierra

La forma teórica de la Tierra se describe por medio de la superficie equipotencial normal de la Tierra
coincidente con la superficie del mar y denominada geoide. El geoide involucra las variaciones del
potencial, que originan entre otro en la distribución irregular de las masas en y encima de la corteza
terrestre. El geoide se puede describir solo aproximadamente. La aproximación más sencilla es el esferoide
definido por la función esférica, que se interrumpe usualmente después los términos cuadrados, puesto
que los resultados ya se vuelven satisfactorios para su aplicación en la gravimetría.

La figura geométrica de la Tierra se aproxima gruesamente por una esfera y con suficientemente exactitud
por un elipsoide de rotación.

Las variaciones entre el geoide (forma teórica) y el elipsoide de rotación se llama las ondulaciones del
geoide y son una medida para la distribución irregular de las masas con respecto al elipsoide de rotación.
Una ondulación de geoide positivo indica un exceso de masa, una ondulación de geoide negativo implica
un déficit de masa.

3.4 Gravedad normal g0

La gravedad normal g0 o es decir el campo gravitacional normal de la Tierra se refiere al elipsoide de

rotación, se la calcula con la formula siguiente:

g0 = 978, 049 (1 + 0,0052884sen2b - 0,0000059sen22b), donde b = latitud geográfica.

Esta formula, llamada formula internacional de gravedad se basa en un valor absoluto de g =

981,274cm/s2 (Gal) medido por KÜHNEN y FURTWÄNGLER en Potsdam en 1906. La formula fue adoptada
por la Unión Internacional de Geodesía y Geofísica en 1930.

En la tabla siguiente se presenta algunos valores de la gravedad normal g 0 y de la variación de la

aceleración de la gravedad correspondientes a distintas latitudes (b).

Aceleración de la
Latitud geográfica ß Gravedad normal g0
gravedad en
en º en mgal
mgal/km
0 978049,0 0
15 978394,0 0,406
30 979337,8 0,704
45 980629,4 0,812
60 981923,9 0,704
75 982873,4 0,406
90 983221,3 0

3 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

Valores según GASSMAN & WEBER (1960).

Los valores máximos de la gravedad normal se observa en los polos debido a la ausencia de la fuerza
centrifuga en estos puntos y debido al aplanamiento de la Tierra. Un cuerpo cayendo sin ser estorbado
encima de uno de los polos aumenta su velocidad en la dirección vertical más rápidamente que el mismo
cuerpo cayendo encima del ecuador hacia el suelo. Expresado en variaciones de masa un cuerpo de 1g de
masa pesa casi 5mg más en los polos que en el ecuador.

4. Reducciones

Con el fin de comparar los valores de gravedad medidos a la superficie terrestre con la gravedad normal se
corrige los valores de gravedad observados. En lo siguiente se introduce las reducciones comunmente
aplicadas a los datos gravimétricos tomados en terreno. Un valor reducido es igual al valor observado de la
gravedad menos el valor previsto de la gravedad basándose en el modelo terrestre elegido. En
consecuencia una anomalía es la diferencia entre lo observado y lo previsto de acuerdo con el modelo
terrestre aplicado.

Para cada reducción se entrega una explicación general y unos cálculos ejercidos con valores procedentes
del siguiente caso de modelo.

Se ha realizado mediciones gravimétricas a lo largo de un perfil orientado Norte-Sur. La estación de base

se ubica en una latitud geográfica ß = 27º30' y en el extremo Norte del perfil. Todas las demás estaciones
están más al Sur. La variación de la gravedad con la distancia Norte-Sur es 0,653mgal/km. El factor de
conversión para la lectura del gravímetro es 1,0255mgal/subdivisión de escala. La densidad media del área
es 2,70g/cm3. Se llevara a cabo la conversión a mgal, la reducción para la deriva del instrumento y para
las mareas, las reducciones para la latitud y para la altura y la reducción con la losa de Bouguer. Se finaliza
los cálculos con un perfil de la anomalía de Bouguer en función de la distancia.

Distancia con
No.Estación respecto a EB Cota sndm en m
1 (EB) 0 250
2 1 240
3 2 220
4 3 200
5 4 195
6 5 190
7 6 205
8 7 245
9 8 275
10 9 280

4 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

4.1 Calibración

Las lecturas de los gravímetros se presentan en subdivisiones arbitrarias de escala, lo que requiere una
calibración para expresarlas en mgal. La calibración se lleva a cabo realizando mediciones en dos lugares
de gravedad absoluta o relativa precisamente conocida por ejemplo obtenida por observaciones de la
gravedad a partir de un péndulo. En el caso de dos lugares de gravedad precisamente conocida se supone
una repuesta linear y se calibra la escala entera en base de los dos valores conocidos. Teniendo en cuenta
más lugares con valores de gravedad precisamente conocida se aumenta la precisión de las mediciones
gravimétricas.

Las reducciones de ejemplo se basan en un factor de conversión de 1,0255mgal/subdivisión de escala

(sdde).

Tiempo en g1convertida =
Lectura en
No. Estación horas y g1 = gEO-gEB g1×1,0255 en
sdde
minutos mgal
1 (EB) 8:30 3575 0 0
2 8:45 3576,57953 1,5795344 1,61981253
3 9:15 3580,024 5,024 5,152112
4 9:45 3581,5211 6,5211 6,68738805
5 10:15 3579,2018 4,2018 4,3089459
6 10:45 3576,9135 1,9135 1,96229425
7 11:15 3570,8194 -4,1806 -4,2872053
8 11:45 3562,8804 -12,1196 -12,4286498
9 12:15 3557,3069 -17,6931 -18,1442741
10 12:45 3556,8697 -18,1303 -18,5926226

5 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

4.2 Reducción para la deriva del gravímetro

Cuando se repite una medición por un gravímetro en el mismo lugar y bajo las mismas circunstancias
varaias veces en el día se obtiene distintas lecturas variando unos 0,001mgal. Estas variaciones se deben
a los resortes y fibras de torsión no perfectamente elásticos de los gravímetros, al efecto de la
temperatura y a la influencia de las mareas. La llamativa deriva del instrumento se corrige observando los
valores de gravedad correspondientes a la estación de base en intervalos discretos de tiempo, por ejemplo
cada media hora y durante el intervalo de tiempo, en que se realiza todas las mediciones. La primera y la
ultima medición se debe realizar en esta estación.

Se delinea las variaciones observadas en la estación de base construyendo una curva de las variaciones de
gravedad en la estación de base en función del tiempo. Se elige un valor como valor de referencia, por
ejemplo el primero valor observado en la estación de base. Los valores para corregir los valores
observados son las variaciones de la gravedad con respecto al valor de referencia en los instantes de
medición en las estaciones de observación. Se los encuentra por interpolación. Se toma la variación 'valor
interpolado - valor de referencia ' correspondiente a la estación de base en el instante de medición en la
estación de observación.

• En el caso, que el valor interpolado es mayor al valor de referencia se resta la variación 'valor
interpolado - valor de referencia' del valor observado en la respectiva estación de observación.

• En el caso, que el valor interpolado es inferior al valor de referencia se suma la variación encontrada
'valor interpolado - valor de referencia' al valor observado en la respectiva estación de observación.

En el ejemplo graficado el valor de referencia es 3757 subdivisiones de escala (sdde). Las variaciones de
gravedad con respecto al valor de referencia en los instantes de medición en las estaciones de observación
están presentadas por barras. Una barra dirigida hacia arriba es un valor de gravedad mayor al valor de
referencia, que se midiera en la estación de base en el instante de medición en la estación de observación.
A las 8:45 horas en la estación de base se midiera un valor de 3575,015sdde (véase la tabla de valores
puesta abajo). La diferencia 'valor de referencia - valor interpolado' = 3575sdde - 3575,015sdde =
-0,015sdde es negativa. En consecuencia se la resta del valor observado en la estación de observación en
cuestión. Una barra dirigida hacia abajo indica un valor menor al valor de referencia, que se midería en la
estación de base en el instante de medición en una otra estación de observación en cuestión. A las 10:15
horas en la estación de base se tomaría la lectura 3574,965sdde. La diferencia 'valor de referencia - valor
interpolado' = 3575sdde - 3574,965sdde = 0,035sdde se suma al valor observado en la estación de

6 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

observación en cuestión.

Variación
c1convertido g2 =
Tiempo de c1 =
= Valor g1convertida
(MEZ) gravedad Valor de valor
No. de para
en en la referencia para g1convertida +
estación corregir
horas y estación en sdde corregir c1convertido
convertido
minutos de base en en mgal
en mgal
en sdde sdde
1(EB) 8:30 3575 3575 0 0 0 0
2 8:45 3575,015 3575 -0,015 0,0153825 1,61981253 1,60443003
3 9:15 3574,992 3575 0,008 0,008204 5,152112 5,160316
4 9:45 3574,972 3575 0,028 0,028714 6,68738805 6,71610205
5 10:15 3574,965 3575 0,035 0,0358925 4,3089459 4,3448384
6 10:45 3574,989 3575 0,011 0,0112805 1,96229425 1,97357475
7 11:15 3575,013 3575 -0,013 -0,0133315 -4,2872053 -4,3005368
8 11:45 3575,024 3575 -0,024 -0,024612 -12,4286498 -12,4532618
9 12:15 3575,01 3575 -0,01 -0,010255 -18,1442741 -18,1545291
10 12:45 3574,985 3575 0,015 0,0153825 -18,5926226 -18,5772401

Para obtener mediciones gravimétricas de alta precisión en cada estación de observación se realiza varias
mediciones repetitivas.

4.3 Reducción para las mareas

Debido a su alta sensibilidad los instrumentos de medición de la gravedad aptos para la prospección
responden a la atracción gravitatoria del sol y de la luna y registran las variaciones periódicas en la
atracción causadas por los movimientos de la Tierra con respecto al sol y a la luna. Estas fuerzas
regularmente alzan y bajan las aguas terrestres no rígidas en ciclos mareales previsibles. Las mismas
fuerzas actúan en la superficie terrestre sólida deformándola en la misma manera como el agua, aún en
dimensiones más pequeñas. El movimiento actual mareal de un punto ubicado en la superficie terrestre
firme alcanza solo unos centímetros. Este desplazamiento mismo causa pequeñas variaciones en la
gravedad debido a la variación de la distancia con respecto al centro de la Tierra y debido a la distribución
nueva de las masas en el interior de la Tierra. Estos cambios están superpuestos a las variaciones
causadas por las fuerzas atractivas de los cuerpos del sistema Tierra - Sol - Luna.

La magnitud de estos cambios varía con la latitud, con el día y el mes. El ciclo entero de las mareas está
caracterizado por variaciones iguales y menores a 0,2 a 0,3mgal. Su variación media es +/-0,1mgal.

El efecto gravitatorio de las mareas calculado para una Tierra rígida y amplificado con un factor de 1,16
aproximará valores actuales de la variación gravitatoria causada por las mareas para la mayoría de los
lugares del mundo en un rango de unos porcientos.

La corrección para la influencia de las mareas se puede realizar recurriéndose a tablas publicadas
regularmente y con anticipación para cada año por ejemplo en el servicio francés llamado 'Service
Hydrographique de la Marine y Compagnie General de Géophysique: Tidal Gravity corrections.

En el ejemplo la corrección para la influencia de las mareas está incorporada en la corrección para la
deriva del instrumento.

4.4 Corrección para la latitud

7 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

La formula internacional de gravedad toma en cuenta la variación de la gravedad normal con la latitud
(véase capítulo 3.4). Otra manera de corregir la variación de la gravedad con la latitud consiste en referir
los valores medidos en las estaciones de observación a la estación de base de coordenadas geográficas
conocidas. Se aplica la formula siguiente:

dgß = (0,0081 × sen2ß)mgal/10m de distancia N-S,

donde ß = latitud geográfica de la estación de base.

En el caso de una estación de base más cercana al polo en comparación con la estación de base se resta el
valor de corrección del valor observado: Dg = gobserv - dgß.

En el caso de una estación de observación ubicada más cerca al ecuador en comparación con la estación
de base se suma el valor de corrección al valor medido: Dg observ + dgß.

El perfil se ubica en el hemisferio Sur y todas las estaciones de observación están al Sur de la estación de
base. Por consiguiente para cada estación de observación el valor de corrección dg ß o c2 se resta del valor
corregido para la deriva y convertida en mgal y referido al valor medido en EB a las 8:30horas (g 2). La
variación de la gravedad con la distancia Norte-Sur es 0,653mgal/km en el área, en que se han trazado el
perfil. Para la estación de observación No.5 por ejemplo se calcula c 2 = (0,653 × 4)mgal = 2,612mgal. La
gravedad corregida para la deriva, las mareas y para la latitud g 3 correspondiente a la estación de
observación No.5 es: g3 = g2 - c2 = 4,3448384mgal - 2,612mgal = 1,7328384mgal. El valor de corrección
c2 se resta de la gravedad corregida para la deriva y de las mareas g 2, puesto que la estación de
observación está más cerca al polo Sur en comparación con la estación de base (hemisferio Sur).

c2 valor de g2 =
Distancia
corrección g1convertida + g3 = g2 - c2
No. de N-S con
para la
estación respecto a c1convertido en en mgal
latitud en
EB mgal
mgal
1(EB) 0 0 0 0
2 1 0,653 1,60443003 0,95143003
3 2 1,306 5,160316 3,854316
4 3 1,959 6,71610205 4,75710205
5 4 2,612 4,3448384 1,7328384
6 5 3,262 1,97357475 -1,29142525
7 6 3,918 -4,3005368 -8,2185368
8 7 4,571 -12,4532618 -17,0242618
9 8 5,224 -18,1545291 -23,3785291
10 9 5,877 -18,5772401 -24,4542401

4.4 Corrección para la altura

La intensidad de la gravedad varía en relación inversa al cuadrado de la distancia 'centro de la Tierra -

estación de observación'. Refiriendo la variación de la gravedad al modelo esférico de la Tierra en reposo el
nivel de referencia dista 6367,5km igual (al radio de la Tierra) con respecto al centro de la Tierra.

Dgalt = 2 × f ×Mtierra/(rtierra)3 × a = a × 0,3083mgal/m, donde

• f = constante de gravitación = 6,67 × 10-8cm3g-1s-2,

8 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

• Mtierra = 5,977 × 1027g

• rtierra = 6367,5km

• a = altura sndm
Dgalt recibe un signo positivo para estaciones de observación situadas encima del nivel de referencia. Dg alt
lleva un signo negativo para estaciones de observación ubicadas debajo del nivel de referencia.

Aplicando el modelo del elipsoide terrestre resultaría un valor de corrección, que coincide con el valor
introducido las tres primeras decimales. Para lograr una precisión de 0,01mgal para las observaciones
gravimétricas las diferencias de altura con respecto al nivel de mar debería ser conocidas con un error
menor a 4cm.

Al efecto de la altura se denomina efecto de aire libre, puesto que las masas de rocas presentes o no
presentes entre el nivel de referencia y el nivel de la estación de observación no lo influyen.

En el ejemplo se refiere la variación de la gravedad con la altura al nivel de mar. Todas las estaciones de
observación se ubican arriba del nivel de mar. Por consiguiente se suma dgalt al valor de gravedad
corregida para la deriva, las mareas y la latitud g3 y se calcula: g4 = g3 + Dgalt. Para la estación de
observación No.6 por ejemplo la variación por la altura con respecto al nivel de mar es Dg alt = 190m ×
0,3083mgal/m = 58,577mgal. La anomalía de aire libre correspondiente a la estación de observación No.6
es g4 = -1,29142525mgal + 58,577mgal = 57,2855748mgal.

Cota g4 = g 3 +
No. de g3 = g2 - c2
sndm en Dgalt en mgal Dgalt en
estación en mgal
m mgal
1(EB) 0 250 77,075 77,075
2 0,95143003 240 73,992 74,94343
3 3,854316 220 67,826 71,680316
4 4,75710205 200 61,66 66,417102
5 1,7328384 195 60,1185 61,8513384
6 -1,29142525 190 58,577 57,2855748
7 -8,2185368 205 63,2015 54,9829632
8 -17,0242618 245 75,5335 58,5092382
9 -23,3785291 275 84,7825 61,403971
10 -24,4542401 280 86,324 61,8697599

4.5 Corrección topográfica

En preparación

4.6 Corrección con la losa de Bouguer

La corrección con la losa de Bouguer elimina el efecto de las masas de rocas ubicadas entre el nivel de
referencia y la estación de observación. Las masas de rocas ejercen una atracción gravitatoria extra a una
estación de observación situada en una altura mayor a aquella de la estación de base. El incremento
esperado de la atracción gravitatoria debido a las rocas ubicándose entre el nivel de referencia y el nivel de
la estación de observación a menudo se modela utilizando el hipótesis que se puede aproximar la roca con
una losa horizontal de dimensiones infinitas y de densidad uniforme, cuyo piso coincide con el nivel de
referencia y cuyo techo está en el nivel de la estación de observación. Se emplea la formula siguiente para

9 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

cuantificar el efecto de Bouguer:

DgBouguer = 2 × (3,14159) × f × d × a = 0,04191 ×d × a [en mgal], donde

• f = constante de gravitación = 6,67 × 10-8cm3g-1s-2,

• d = densidad de las masas de rocas ubicándose entre el nivel de referencia y el nivel de la estación
de observación en g/cm3,

• a = diferencia de altura entre el nivel de referencia y el nivel de la estación de observación en m.

La corrección con la losa de Bouguer DgBouguer se resta del valor observado en una estación de
observación en el caso que la estación de observación está encima de la estación de base. Se la suma al
valor observado en el caso que la estación de observación se sitúa debajo del nivel de referencia.
Referente a los signos, la corrección para la altura o es decir de aire libre siempre lleva el signo opuesto al
signo aplicado para la corrección con la losa de Bouguer.

La corrección con la losa de Bouguer se basa en una densidad uniforme, que se supone para las masas
rocosas ubicadas entre el nivel de referencia y el nivel de la estación de observación. La mayoría de los
casos reales no cumple esta condición. Por ejemplo en un área sedimentaria construida por una secuencia
de estratos horizontales de distintas densidades se puede modificar la formula para la corrección de
Bouguer de la manera siguiente:

DgBouguer = 0,04191 (d1 × a1 + d2 × a2 + d3 × a3 +... +... + di × ai) mgal, donde

• d1 a di = densidades de los estratos 1 a i,

• a1 a ai = anchos (potencias) de los estratos 1 a i.

La densidad media del área sea 2,7g/cm3. En el primer paso se realiza la corrección con la losa de
Bouguer con respecto al nivel de mar. Por consiguiente la potencia de la losa de Bouguer coincide con las
cotas correspondientes a las estaciones de observación. Como cada una de las estaciones se ubica sobre el
nivel de mar se debe restar la corrección con la losa de Bouguer del valor de gravedad corregida para la
deriva, las mareas, la latitud y la altura g4 y la anomalía de Bouguer con respecto al nivel de mar es g 5 =
g4 - DgBouguer. Para la estación de observación 7 se obtiene por ejemplo:

DgBouguer = 0,04191 × 2,7g/cm3 × 205m = 23,197185mgal.

g5 = 54,9829632mgal - 23,197185mgal = 31,7857782mgal

La variación de la anomalía de Bouguer con respecto a la estación de base g 6 se calcula como sigue: g6 =
g5EO - g5EB, donde EO = estación de observación y EB = estación de base.

g6 =
DgBouguer
Cota g4 = g 3 + g5 = g4 - g5EO -
No. de en mgal
sndm en Dgalt en DgBouguer g5EB
estación (respecto
m mgal en mgal en
al ndm)
mgal
1(EB) 250 77,075 28,28925 48,78575 0,00
2 240 74,94343 27,15768 47,78575 -1,00
3 220 71,680316 24,89454 46,785776 -2,00
4 200 66,417102 22,6314 43,785702 -5,00
5 195 61,8513384 22,065615 39,7857234 -9,00

10 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

6 190 57,2855748 21,49983 35,7857448 -13,00

7 205 54,9829632 23,197185 31,7857782 -17,00
8 245 58,5092382 27,723465 30,7857732 -18,00
9 275 61,403971 31,118175 30,285796 -18,50
10 280 61,8697599 31,68396 30,1857999 -18,60

El grafico ilustra la anomalía de Bouguer en función de la distancia con respecto al nivel de mar.
Comparándolo con el perfil morfológico correspondiente no se nota ninguna correlación. Para una curva de
tal forma existen varias interpretaciones (en preparación). Cuál de estas interpretaciones semeja más a la
realidad se deberá comprobar por medio de otros estudios geofísicos y/o geológicos.

4.7 Anomalías de gravedad

Una anomalía de gravedad se define como la variación de los valores medidos de la gravedad con respecto
a la gravedad normal después de haber aplicado las correcciones necesarias.

La anomalía de aire libre resulta de las correcciones de la influencia de las mareas, de la derive del
instrumento de medición, de la latitud y de la altura.

La anomalía de Bouguer se obtiene aplicando todas las correcciones mencionadas.

5. Densidad

5.1 Determinación de la densidad del área

La densidad media del área en consideración entra en las formulas, que corrigen el efecto topográfico y el
efecto de las masas ubicadas entre el nivel de referencia y el nivel de observación (corrección con la losa
de Bouguer). Por consiguiente el conocimiento de la densidad media del área en consideración contribuye
a la reducción o eliminación de dichos efectos, además el conocimiento de la distribución de la densidad en

11 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

el área de interés es uno de los fundamentos de la interpretación de los perfiles o mapas gravimétricos
resultantes de las mediciones. Los resultados de las mediciones gravimétricas pueden ser ambiguos como
muestra el ejemplo siguiente. El perfil gravimétrico de la figura (en preparación) característico para una
flexura o falla puede ser causado por diferentes distribuciones de la densidad en la profundidad.

La determinación directa de la densidad de muestras representativas provenientes de afloramientos, minas

o sondeos se realiza en el laboratorio por medio de un picnómetro o una balanza de SCHWARZ o JOLLY.
Precisamente se mide la muestra de roca en aire y en agua y se calcula su densidad 'd' de la manera
siguiente:

dmuestra = peso de la muestra en aire/(peso de muestra en aire - peso de la muestra en agua).

De tal manera se puede determinar la densidad de muestras de rocas compactas, no porosas, como de
rocas plutónicas y metamórficas.

En el caso de las rocas sedimentarias su densidad depende del grado, en que sus poros están llenados con
agua. Una muestra porosa se procesa de la manera siguiente: se la satura con agua y se la pesa en aire y
sumergida en agua. Luego se la seca totalmente en un armario de secado y se la pesa de nuevo en aire y
sumergida en agua. Los distintos pesos encontrados se insertan en la formula ya mencionada. La densidad
verdadera de la muestra se ubica en el intervalo limitado por los dos valores extremos calculados
correspondientes a la muestra saturada con agua y a la muestra totalmente secada.

Estas determinaciones de densidad carecen de que las muestras de algunos afloramientos puntuales no
necesariamente son representativos para toda el área. Además las muestras superficiales pueden variar
apreciadamente en su humedad y en su grado de meteorización en comparación a las muestras ubicadas
en una profundidad más alta, en el caso de rocas sueltas como arcillas, margas, depósitos de morrenas las
rocas superficiales pueden ser menos compactadas en comparación a aquellas ubicadas en una
profundidad más alta.

NETTLETON propuso el siguiente método indirecto. Se considera un perfil gravimétrico trazado sobre un
accidente morfológico pronunciado en el área de interés como una colina o un valle pequeño. Para cada
estación de observación a lo largo del perfil se calcula la gravedad corregida insertando distintos valores de
densidad en las formulas, que corrigen el efecto topográfico y el efecto de las masas ubicadas entre el
nivel de referencia y el nivel de observación. El valor de densidad, que genera el perfil gravimétrico de
menor correlación con el perfil morfológico, es el valor más apropiado y él, que se acerca lo más posible al
valor real.

JUNG ha transferido el método gráfico de NETTLETON al lenguaje matemático. La densidad, que genera un
perfil gravimétrico de menor correlación con el perfil morfológico también se puede hallar suponiendo, que
no existe ninguna correlación entre la morfología y los valores de gravedad. Insertando 0 para el cociente
de correlación resulta la formula siguiente:

d = d0 + [S(DgB0i -DgB)× (ai - a)] /[2 × 3,1415927 × f × (S(ai - a)2)] ,donde

• d = densidad buscada.

• d0 = densidad estimada para el área en consideración.

• f = constante de gravitación = 6,67×10-8 cm3g-1s-2

• S = suma de i=1 hasta n.

• DgB0i = anomalía de Bouguer correspondiente a la estación de observación i.

• ai = altura de la estación de observación i

• D gB = promedio aritmético de las anomalías de Bouguer de todas las estaciones de observación del
perfil

12 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

• a = promedio aritmético de las alturas correspondientes a todas las estaciones de observación del
perfil.

Del método de NETTLETON y de su modificación por JUNG resulta un promedio del efecto de la variación
de densidad más preciso en comparación a la determinación de densidad de muestras superficiales en el
laboratorio. Sin embargo el método de NETTLETON está limitado a profundidades relativamente someras y
a litologías homogéneas.

Hoy día en pozos de sondeos se aplican instrumentos de medición denominados 'density loggers' o sondas
de rayos gamma, que entregan una diagrafía de densidad de las distintas formaciones geológicas. La
sonda se constituye de una fuente radiante (rayos gamma), normalmente de cobalto 60, ubicada en el
extremo inferior de la sonda y un detector, normalmente un contador de Geiger, instalado en el extremo
superior en una distancia de aproximadamente 45cm con respecto a la fuente radiante. La sonda está
envuelta por una capa de plomo con dos orificios posicionados en los niveles de la fuente y del detector de
tal modo, que la única radiación, que puede llegar al detector es aquella reflejada de la formación
geológica por la dispersión del tipo Compton. La amplitud de la radiación dispersada depende de la
concentración de electrones de la formación geológica, la cual es aproximadamente proporcional a la
densidad de la formación geológica. La máxima penetración de los rayos gamma tiene un alcance de 15
cm hacia las rocas adyacentes y el volumen efectivo, que capta la sonda mediante un intervalo de
medición es 0,03m3. Para mantener un contacto estrecho entre un lado de la sonda y uno de las paredes
rocosas del pozo el otro lado de la sonda lleva un resorte. La comparación de los resultados de una sonda
de rayos gamma ('density logger') con las determinaciones de densidad de los testigos correspondientes a
los mismos niveles muestra una coincidencia hasta unas centésimas partes de un gramo por centímetro
cúbico para todas las formaciones geológicas excepto las rocas arcillosas, ricas en minerales arcillosas y
otras rocas muy blandas. Estas rocas tienden a ser socavado por el agua o el lodo de perforación lo que
influye desfavorablemente las mediciones. Las mediciones en pozos con la sonda de rayos gammas son
limitadas a volúmenes rocosos pequeños, solo representativos en el caso de formaciones litológicas
homogéneas y deberían ser comprobadas mediante determinaciones de densidad en el laboratorio.

Un gravímetro especialmente apropiado para pozos permite realizar mediciones de densidad versus la
profundidad para un volumen rocoso mayor en comparación al volumen rocoso cubierto por la sonda de
rayos gamma. La densidad se obtiene a través de la diferencia en gravedad medida en dos niveles del
pozo. En general el espaciamiento (distancia entre los dos niveles, donde se toma la lectura) es alrededor
de 3m. Las densidades obtenidas con este método son representativas para un volumen rocoso mayor en
comparación con aquel captado por la sonda de rayos gamma y pueden ser incorporados en mediciones
gravimétricas realizadas en la superficie.

5.2 Densidades de rocas y minerales

Tipo de roca Densidad media en g/cm3 Rango de densidad en

g/cm3

Densidades medias de rocas plutónicas1

Granito 2,667 2,516 - 2,809

Granodiorita 2,716 2,668 - 2,785

Sienita 2,757 2,630 - 2,899

Diorita cuarcífera 2,806 2,680 - 2,960

Diorita 2,839 2,721 - 2,960

Norita 2,984 2,720 - 3,020

Gabro, con gabro de 2,976 2,850 - 3,120

olivino

13 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

Diabasa 2,965 2,804 - 3,110

Peridotita no alterada 3,234 3,152 - 3,276

Dunita2 3,277 3,204 - 3,314

Piroxenita 3,231 3,100 - 3,318

Anortosita 2,734 2,640 - 2,920

Densidades medias de rocas sedimentarias3

Arenisca 2,32 1,61 - 2,76

Pizarra arcillosa (shale) 2,42 1,77 - 2,45

Caliza 2,54 1,93 - 2,90

Dolomia 2,70 2,36 - 2,90

Densidades medias de rocas metamórficas4

Gneis, Chester, Vermont 2,69 2,66 - 2,73

Gneis granítico 2,61 2,59 - 2,63

Gneis, Grenville 2,84 2,70 - 3,06

Gneis con oligoclasa 2,67

Esquisto de cuarzo y 2,82 2,70 - 2,96

mica

Esquisto de muscovita 2,76

y biotita

Esquistos de estaurolita 2,76

y granate y de biotita y
muscovita

Esquistos de clorita y 2,82 2,73 - 3,03

sericita

Slate 2,81 2,72 - 2,84

Anfibolita 2,99 2,79- 3,14

Granulita con 2,93 2,67 - 3,10

hiperstena

Granulita sin hiperstena 2,73 2,63 - 2,85

Eclogita 3,392 3,338 - 3,452

1: Valores de DALY, R.A. (1966): Handbook of Physical Constants. - Ed.: Geol. Soc. Am. Mem., 97. en
DOBRIN (1988): p.610.

2: Valores de BIRCH (1960), J. Geophys. Res., col.65, p.1083 en DOBRIN (1988): p.610.

3: Valores de CLARK, S.P. (1966): Handbook of Physical Constants. - Ed.: Geol. Soc. Am. Mem., 97. en
DOBRIN (1988): p.610.

4: Mobil Oil Co. en DOBRIN (1988): p.610.

5.3 Isostasia

14 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

Si las montañas estuviesen colocadas sobre una capa rígida de la Tierra, las montañas más altas se
derrumbarían debido a su alto peso generando en altas profundidades variaciones de presión mayores a la
resistencia a la presión del material rocoso. Para explicar, por qué las variaciones topográficas de la
superficie terrestre persisten un intervalo de tiempo largo se asume, que la porción rígida de la Tierra,
llamada litosfera es una capa exterior delgada en comparación con el radio de la Tierra y que la litosfera
flota sobre un interior altamente viscoso. Este concepto requiere, que una masa sobresaliente del nivel de
mar tiene que ser compensada por un déficit de masa debajo del nivel de mar y que las cuencas oceánicas
anormalmente livianas deben ser acompañadas por masas de alta densidad ubicadas en la profundidad. Si
no existiese esta compensación las montañas se estarían hundiendo y las cuencas marinas alzarían. El
peso total por cada unidad de área debajo de la litosfera debe ser uniforme, si el equilibrio isostático está
realizado.

Las montañas grandes como por ejemplo los Andes o la Himalaya no muestran una anomalía de gravedad
tan grande como se podría esperar debido a sus dimensiones. Considerando la estructura estratificada de
la Tierra (corteza terrestre, manto, núcleo) se nota, que a las montañas de altura grande corresponden
raíces, que se hunden profundamente en el manto superior. Las montañas altas casi flotan en el material
más denso del manto superior y sus raíces menos densas hundidas en el manto superior más denso
compensan con sus fuerzas ascendentes sus partes superiores. La porción rígida de la Tierra (litosfera) es
una capa muy delgada en comparación con el radio de la Tierra y flota en un interior efectivamente liquido
y altamente viscoso. Este concepto requiere que el exceso de masas debido a las masas ubicadas entre la
superficie terrestre y el nivel de mar y el déficit de masas presente entre el nivel de mar y el fondo de mar
están compensados por masas de signo reverso ubicadas en la profundidad de la Tierra. En equilibrio en
cualquier lugar de cualquier profundidad debajo de la litosfera el peso total por área unitaria debería ser
igual. A este estado de equilibrio se llama isostasía (DUTTON).

Dos modelos distintos explican la isostasía y la forma de la compensación isostática: El modelo de PRATT
(1809 - 1871), el modelo de AIRY, G.B. (1801 - 1892).

Según PRATT el material menos denso de las raíces de las montañas tiene una base uniforme y las
variaciones en la altura de las partes montañosas se basan en variaciones de densidad del material que la
constituye.

AIRY supone una densidad uniforme para el material, que constituye las montañas y los diferentes niveles
sobresalientes de las montañas se compensan con partes distintamente profundas hacia abajo.

15 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

6. Métodos e instrumentos de medición de la gravedad

1. Péndulo

La medición de gravedad por medio de un péndulo (péndulo de reversión) es un método absoluto. Para el
péndulo físico vale:

T = 2 ´ 3,1415927 [Qc/ (m × g × h)]1/2] , donde

• Qc = momento de inercia del eje de rotación c

• m = masa total del péndulo

• h = distancia desde el centro de gravedad al centro de rotación.

• g = aceleración de gravedad.
Con todos los demás parámetros conocidos de puede deducir la aceleración de gravedad.

2. Experimento de caída

Por el experimento de caída se recibe valores absolutos de la aceleración de gravedad. Para un intervalo de
tiempo T y el intervalo de altura, por lo cual pasa un cuerpo durante el intervalo de tiempo T: x = x 0 + u ×
T + 1/2 × g× T2, donde

• u = velocidad inicial en el nivel x0

• g = aceleración de gravedad
Con configuraciones adecuadas de este experimento de caída se logra determinaciones de los parámetros
x, x0, u y T suficientemente exactas para el calculo de la aceleración de gravedad.

En la prospección gravimétrica se han utilizado tres tipos de instrumentos:

1. Péndulo

2. Balanza de torsión

3. Gravímetro
1. Péndulo

2. La balanza de torsión

La balanza de torsión mide gradientes y curvaturas en lugar de aceleraciones gravitatorias.

Balanzas de distintas configuraciones son:

1. La balanza de CAVENDISH

2. La balanza de torsión estándar de EÖTVÖS

3. La variante de la balanza de torsión estándar

4. La balanza de barra inclinada

La balanza de EÖTVÖS

La balanza de EÖTVÖS está equipada con dos pesos iguales situados a distintas alturas y unidos

16 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

solidariamente. Este conjunto está suspendido de un hilo de torsión de tal manera que la construcción
puede girar libremente en torno del hilo en el plano horizontal. En su disposición más común el soporte es
una barra ligera. Una de las masas reposa en uno de los extremos de la barra, la otra masa suspende del
otro extremo de la barra. La barra gira solamente cuando actúa una fuerza diferencial horizontal en ella o
es decir cuando el campo gravitatorio terrestre de las proximidades del instrumento está distorsionado de
tal manera que las componentes horizontales en los extremos de la barra difieren.

En un campo gravitatorio que pudiera representarse por superficies equipotenciales planas y paralelas no
habría ninguna rotación de la balanza puesto que las componentes horizontales actuando en los dos
extremos de la balanza serían iguales. Una deformación de las superficies equipotenciales como puede
originarse por la atracción de una masa enterrada haría girar la balanza con una magnitud de rotación,
que depende de la magnitud de la fuerza horizontal no compensada y de la rigidez del hilo de torsión. El
par de fuerzas opuestas ejercido por el hilo de torsión es proporcional al ángulo de rotación medido con
respecto a su posición no torsionada. Las masas de la balanza de torsión se desplazan paralelamente a la
superficie terrestre y mediante un movimiento giratorio desde una zona de potencial gravitatorio alto hacia
una zona de potencial gravitatorio más bajo.

La variante de la balanza de torsión estándar de EÖTVÖS empleada en trabajos de campo se constituye de

dos barras paralelas de unos 40cm de largo con dos pesos sujetados en los dos extremos de cada una de
las barras. Los dos pesos superiores están apoyados en una de las barras, los dos pesos inferiores
suspenden a unos 60 a 70cm por debajo de la otra barra. Cada peso tiene una masa de 25g. Las
rotaciones pequeñas causadas por el efecto de una fuerza diferencial horizontal se amplifican por medio de
un sistema óptico.

En una estación de observación se debe tomar por lo menos seis lecturas: se orienta el par de barras en
tres direcciones separadas entre sí 120º y se realiza por lo menos una lectura para cada barra en cada una
de las tres distintas direcciones. Estas lecturas posibilitan la determinación del gradiente de gravedad y la
curvatura.

En áreas favorables la balanza de torsión puede alcanzar una precisión semejante a aquella de gravímetros
modernos. No obstante hoy día la balanza de torsión no está más en uso por la cantidad de lecturas
necesarias para una estación de observación y por el tiempo gastado en estas lecturas.

Gravímetros estables

Gravímetro de HARTLEY

El gravímetro de HARTLEY de constituye de un peso suspendido de un resorte. Por variaciones en la

aceleración gravitatoria de un lugar al otro el resorte principal se mueve y puede ser vuelto a su posición
de referencia por medio de un movimiento compensatorio de un resorte auxiliar o de regulación manejable
por un tornillo micrométrico. El giro del tornillo micrométrico se lee en un dial, que da una medida de la
desviación del valor de la gravedad con respecto a su valor de referencia. Por la posición del espejo en el
extremo de la barra, su desplazamiento es mayor que el desplazamiento del resorte principal y como el
recorrido del haz luminoso es grande, se puede realizar medidas de precisión cercanas al miligal.

Se sigue preparando este capitulo.

7. Interpretación

8. Aplicación y ejemplos

Por medio del método gravimétrico se puede detectar contrastes de densidad existentes en la corteza
terrestre debido a los distintos tipos de rocas, que constituyen la corteza terrestre. La mayoría de las rocas
sedimentarias por ejemplo es menos densa en comparación a las rocas, que forman el basamento. En
consecuencia con el método gravimétrico se puede delinear la interfase o el limite entre las rocas
sedimentarias y las rocas del basamento subyacentes o las dimensiones de cuencas sedimentarias
formando lechos o otras depresiones en las rocas del basamento. El método gravimétrico es muy útil en la
exploración inicial de áreas cubiertas por una capa uniforme, que esconde los afloramientos y la estructura
del subsuelo. La cubierta puede componerse de vegetación densa, de agua somera o de aluviones por
ejemplo. Los fines de la década sesenta un gravímetro portado por barco fue desarrollado, lo que empujó

17 de 18 3/9/2021 09:39
6 https://www.medellin.unal.edu.co/~rrodriguez/geologia/gravimetria.htm

la exploración gravimétrica de los márgenes continentales de cubierta somera de agua. Combinando los
datos gravimétricos con los resultados de la exploración sísmica el geofísico puede identificar más
claramente estructuras y formaciones geológicas como domos de sal o de rocas ígneas por ejemplo en
comparación con la aplicación de solo uno de estos métodos geofísicos de exploración.

En la exploración minera se aplica el método gravimétrico en la búsqueda de minerales pesados como la

cromita por ejemplo. Debido al contraste alto de densidad entre los minerales pesados y las rocas
adyacentes más livianas se puede delinear la distribución y dimensión de las rocas de diferentes
densidades por medio del método gravimétrico. Los canales antiguos son prometedores para
acumulaciones de menas de oro y de uranio. Frecuentemente ellos están hundidos y escondidos debajo de
una cubierta de otras rocas. Debido al contraste de densidad entre el relleno menos denso de estos
canales, que hacen incisiones en rocas de mayor densidad, el método gravimétrico está capaz de delinear
la forma de estos canales.

Los estudios de reconocimiento regional por medio del método gravimétrico pueden resultar en el
levantamiento de estructuras geológicas de importancia regional tales como fallas o lineamientos, que son
prometedores para acumulaciones de minerales y mineralizaciones.

Además se emplea el método gravimétrico para distinguir anomalías electromagnéticas causadas por
sulfuros macizos de aquellas causadas por grafitos de densidad relativamente pequeña, por ejemplo en el
escudo canadiense se realizaron tales estudios.

18 de 18 3/9/2021 09:39