Práctica de Gravimetrı́a

Félix Ayllón Muñoz

Fı́sica de la Tierra

17 de marzo de 2022

Resumen
El objetivo principal de esta práctica es el cálculo de las anomalı́as de aire libre y Bouger
mediante los datos aportados de distancia, altura y gravedad. Se realizarán representaciones gráficas
de las distintas anomalı́as en función de la distancia, por último supondremos que la anomalı́a
esta producidad por una esfera de diferente densidad situada a una determinada profundidad que
calcularemos de dos maneras distintas ası́ como su radio.

1. Introducción teórica
La superficie equipotencial de la Tierra sufre desviaciones debido a las fluctuaciones de densidad
del terreno, los valores conocidos teóricamente de la gravedad son estimados en la superficie como un
elipsoide, debemos corregir las fluctuaciones en el terreno para adecuarlo a la forma del elipsoide.
El valor que difiere entre la gravedad que medimos corregida con estas fluctuaciones y el valor teórico
del elipsoide es lo que llamamos anomalı́a gravitatoria:

∆g = gcorregido − γelipsoide

Recordamos que el valor de la gravedad normal del elipsoide viene dado por la siguiente expresión:

γ = γe (1 + β1 sin(ϕ)2 + β2 sin(2ϕ)2 )

Donde ϕ es la latitud del lugar, γe es la gravedad en el ecuador y β1 y β2 son constantes.

β1 = 5, 30244 · 10− 3
β2 = −5, 8 · 10− 6
En esta práctica nos centraremos en dos tipos de anomalı́as: anomalı́a de aire libre y anomalı́a de la
placa Bouger.

1.1. Anomalı́a de aire libre

Si el punto donde queremos medir la gravedad se encuentra a una diferente altura del elipsoide, ya
sea por encima o por debajo, necesitaremos una corrección que nos corriga este efecto, la corrección de
aire libre corrige la diferencia de alutra entre un punto P y el geoide. Si P se encuentra por encima del
elipsoide, la correción será positiva, ya que el valor de g será menor al estar situado en un radio mayor.
Si por el contrario P se encuentra por debajo del elipsoide, la correción será negativa, ya que el valor de
g será mayor al estar en un radio menor.
La correción de aire libre (F) viene dada por la siguiente expesión:

F = 0, 3086h(mgal)

1
 Siendo h la altura en metros.
Para calcular la anomalı́a de aire libre (∆gF ) bastará con usar la siguiente relacción:

∆gF = gobs + F − γelip

1.2. Anomalı́a de la placa Bouger

Una vez hemos corregido el efecto de aire libre, es decir la diferencia de altura entre los puntos,
debemos eliminar la influencia gravitatoria de las capas de tierra que se encuentran entre el punto P
y la superficie del elipsioide. Dependiendo de donde estuviera situado el punto P respecto al elipsoide
tendremos un deficit de masa o exceso de masa que marcará el signo de la correción de la placa Bouger.
La correción de placa Bouger (A) viene dada por la siguiente expresión:

AB = 2πGρh = 1, 119h(mgal)

Donde h es la altura en metros.

Podemos calcular la anomalı́a de la placa Bouger de la siguiente manera:

∆gB = gobs + F − γelip − AB = ∆gF − AB

1.3. Anomalı́a producida por una esfera

Podemos suponer que la anomalı́a está producida por una esfera enterrada a una cierta profundidad z y
de un cierto radio R, esta esfera representaria un exceso de densidad ∆ρ respecto al material circundante.
Como hemos visto en clase, si calculamos la gravedad creada por dicha esfera obtenemos la siguiente
expesión:

4 z
∆gz = πG∆ρR3 3
3 (x + z 2 ) 2
2

Donde x es la distancia, podemos observar que el máximo se produce en x=0,

4 ∆ρR3
(∆gz )max = πG 2
3 z

2
 Figura 1: Anomalı́a producida por una esfera.

2. Datos y resultados
2.1. Anomalı́as
La mayorı́a de los calculos se han realizado con el programa Matlab, tanto la tabla con los resultados
obtenidos como el código aparecerá en el anexo.
Se nos ha proporcionado datos de gravedad a cada distancia, lo primero que hemos hecho es
transformar gravedad de gal a mgal que será la unidad en la que trabajaremos, despues calcularemos
la gravedad teórica del elipsoide en el punto que nos dan por latitud y calcularemos las diferentes
correcciones de aire libre y placa Bouger como hemos mencionado en la introducción teórica, para más
adelante obtener las anomalı́as, todo esto realizado para la contribución a cada altura.
Para una mejor visualización hemos representado algunos de los datos más representativos.

Figura 2: Altura vs distancia.

3
 Figura 3: Anomalı́a de aire libre vs distancia.

Figura 4: Anomalı́a de placa Bouger vs distancia.

Figura 5: Comparación entre las dos anomalias.

4
 Como podemos comprobar los valores de la anomalı́a de aire libre son mayores que los de anomalı́a
Bouger, lo cual tiene sentido ya que para aire libre suponemos que no existe masa entre los dos puntos.
Por otra parte observamos como el valor de la anomalı́a Bouger es positiva, lo cual nos indica que es
producida por un exceso de masa en la profundidad.

2.2. Anomalı́a producida por una esfera

En este apartado se nos pide suponer que la anomalı́a es producida por una esfera de radio R y
g
contraste de densidad ∆ρ = 1 cm 3 que está enterrada una profundidad z, debemos calcular los parámetros

R y z de dos maneras diferentes: mediante el dato dado de la relacción entre semianchura y profundidad
y mediante un ajuste de mı́nimos cuadrados.

2.2.1. Relacción de la semianchura

En este caso se nos dice que la semianchura w guarda la siguiente relacción con la profundidad:

z = 0, 652w

Hemos representado la anomalı́a de Bouger frente a la distancia para después realizar un ajuste
Gaussiano que nos permita obtener la altura a media anchura y de este modo poder obtener la
profundidad z. (Este ajuste ha sido realizado mediante el programa Origin por facilidad de manejo)

Figura 6: Ajuste gaussiano de la anomalı́a Bouger.

Obtenemos un valore para w=2566m, por lo tanto hallamos una profundidad

z=0,652*2566=1673m.
Sabiendo que el valore de la anomalı́a máxima es de (∆g)max = 6, 9mgal podemos despejar el radio
R de la siguiente expresión:

4 ∆ρR3
(∆gz )max = πG 2
3 z
Obteniendo un valor R=884,25m

2.2.2. Mı́nimos cuadrados

Este método consistirá en realizar un ajuste de mı́nimos cuadrados con los parámetros libres z y R,
de esta manera con ayuda del comando lsqcurvefit de Matlab nos realizará dicho ajuste de la ecuación:

5
 4 z
∆gz = πG∆ρR3 3
3 (x + z 2 ) 2
2

Debemos introducir unos valores iniciales de R y z, por lo cual usaremos los obtenidos en el apartado
anterior.
De esta manera hemos obtenido unos valores que minimizan la función de:
R=894,1m
z=1671,3m

Figura 7: Ajuste mı́nimos cuadrados de la anomalı́a Bouger.

6
3. Anexo
Distancia (m) Altura (m) Gravedad (gal)
0 740 980.0625
200 770 980.0570
400 800 980.0515
600 820 980.0481
800 870 980.0388
1000 920 980.0296
1200 930 980.0282
1400 980 980.0189
1600 990 980.0174
1800 1000 980.0157
2000 990 980.0178
2200 950 980.0256
2400 900 980.0351
2600 850 980.0445
2800 850 980.0439
3000 810 980.0518
3200 680 980.0762
3400 620 980.0875
3600 560 980.0988
3800 540 980.1022
4000 550 980.0999

Tabla 1: Datos proporcionados de gravedad a distintas alturas

Distancia (m) Altura (m) Anomalı́a aire libre (mgal) Anomalı́a placa Bouger (mgal)
0 740 85,241 2,435
200 770 88,999 2,836
400 800 92,757 3,237
600 820 95,529 3,771
800 870 101,659 4,306
1000 920 107,889 4,941
1200 930 109,575 5,508
1400 980 115,705 6,043
1600 990 117,291 6,510
1800 1000 118,677 6,777
2000 990 117,691 6,910
2200 950 113,147 6,842
2400 900 107,217 6,507
2600 850 101,187 6,072
2800 850 100,587 5,472
3000 810 96,143 5,504
3200 680 80,425 4,333
3400 620 73,209 3,831
3600 560 65,993 3,329
3800 540 63,221 2,795
4000 550 64,007 2,462

Tabla 2: Valores de las diferentes anomalı́as a las diferentes alturas.

7
%PRÁCTICA GRAVIMETRÍA--Félix Ayllón Muñoz
%(Física de la Tierra)

load datos.txt
z=datos(:,3);
altura=datos(:,2);
distancia=datos(:,1);

gravedad_obs=z*1000 %pasar a mgal

latitud=40.4083; %grados

beta1=5.30244*10^-3;

beta2=-5.8e-6;

gamma_e=9.78032; %gravedad elipsoide en el ecuador

gamma=gamma_e*(1+beta1.*(sind(latitud)).^2+beta2.*(sind(2*latitud)).^2);%gravedad en ms-2

gammaa=gamma*1e5 %gravedad en mgal

%Correccion aire libre

F=0.3086.*altura %mgal

%Anomalia aire libre

deltaf=gravedad_obs+F-gammaa %mgal

%Correcion Bouger

A=0.1119*altura;

%Anomalia Bouger

deltab=gravedad_obs+F-gammaa-A

%Graficas

hold on

figure(1)
plot(distancia,altura,'r.');
xlabel('distancia (m)')
ylabel('altura (m)')
title('altura vs distancia')

hold off

hold on

figure(2)

1
plot(distancia,deltaf,'g.');
xlabel('distancia (m)')
ylabel('anomalía de aire libre (mgal)')
title('anomalía de aire libre vs distancia')

hold off

hold on

figure(3)
plot(distancia,deltab,'b.');
xlabel('distancia (m)')
ylabel('anomalía de Bouger (mgal)')
title('anomalía de Bouger vs distancia')

hold off

figure(4)
plot(distancia,deltaf,'g.');

hold on

plot(distancia,deltab,'b.');
xlabel('distancia (m)')
ylabel('anomalías (mgal)')
title('Comparativa contribución anomalías')
legend('anomalía de aire libre','anomalía de Bouger')

hold off

%Mínimos cuadrados

G=6.67e-11;
rho=1000;

R0=884.25;%Valores iniciales
Z0=1673;
u0=[R0,Z0];

fun=@(u,distancia)(4/3*pi*G*rho.*(u(1).^3).*(u(2)./(distancia.^2+u(2).^2).^(3/2)));

u=lsqcurvefit(fun,u0,distancia,deltab);

Radio=u(1)

profundidad=u(2)

ff=[Radio,profundidad]

dcentro=distancia-2000; %calculamos la distancia centrada en x=0

figure(5)

2
plot(dcentro,deltab,'r.')
hold on
plot(dcentro,fun(ff,dcentro)*10^5,'b') %hemos pasado a mgal multiplicando por 10^5
xlabel('distancia (m)')
ylabel('anomalía Bouger (mgal)')
title('Comparativa anomalía Bouger y ajuste mínimos cuadrados')
legend('anomalía de Bouger','ajuste mínimos cuadrados')
hold off