Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 107

Gravimetría

Para determinar cómo es la configuración del subsuelo. La información que se obtiene no es


directa.
La gravimetría es la medición y análisis del campo gravitatorio terrestre y de sus variaciones en
el espacio y tiempo.

Tiene muchas aplicaciones​:


● Prospección de hidrocarburos, minerales de mena, agua.
● Estructura interna de la Tierra (desde la superficie hasta el núcleo)
● Isostasia y propiedades mecánicas de la litósfera
● Mareas terrestres
● Transferencia de fluidos entre reservorios: agua, magma, hielo (variaciones temporales
de la gravedad).
Nos interesan las variaciones en los valores del campo gravitatorio. Las variaciones indican
contrastes de densidad​. Puede haber dos zonas en el subsuelo con densidades diferentes.
Puede ser por variaciones de litología o por estructuras artificiales.

En estudios geológicos, medioambientales y en ingeniería civil presentan un alto interés las


variaciones de gravedad, las cuales reflejan contrastes en densidad generados por:
● la existencia de estructuras artificiales enterradas a poca profundidad.
● la existencia de distintas formaciones geológicas

La gravimetría se basa en la ​Ley de la gravitación de Newton​. Dos partículas materiales se


atraen con una fuerza de dirección coincidente con la de la recta que las une, de intensidad
directamente proporcional al producto de sus masas e inversamente proporcional al cuadrado
de la distancia que las separa. Un cuerpo esférico atrae a otro en la misma forma que si la
masa estuviera concentrada en su centro.

1
El signo negativo de la ecuación expresa el hecho de que la interacción es atractiva, y la fuerza
aplicada en la masa m1 está por tanto dirigida hacia la masa m2 .

G es la constante de gravitación universal 6,673.10​-11​ N.m​2​ /kg​2


La constante de gravitación G describe la fuerza expresada en N (Newton) ejercida entre dos
cuerpos de masas de 1kg, cuyos centros distan 1m entre sí y cuyas masas están concentradas
en sus centros. Se la mide en laboratorios. En el año 1797, Cavendish realizó por primera vez la
medición de G resultando un valor de G = 6,754 .10​-11​ N.m​2​ /kg​2​ .

Campo Gravitatorio
Puede considerarse que la interacción es fruto de un “cambio" o condición que la masa m1
genera en el espacio que la rodea. La manifestación de este cambio en la región que rodea la
masa recibe el nombre de campo gravitatorio.
Cuando en un punto de este espacio colocamos otra masa de prueba m2, es el campo, y no
directamente la masa m1 el que causa la fuerza de atracción.

Para cuantificar el campo gravitatorio se usa el vector intensidad de campo gravitatorio, que se
define como la fuerza sobre la unidad de masa, de forma que una masa puntual m1=M crea un
campo gravitatorio cuya intensidad es:

La dimensión de esta magnitud es fuerza por unidad de masa, es decir la misma que la
aceleración. Cualquier cuerpo de cualquier masa, experimenta exactamente la misma
aceleración en un determinado campo gravitatorio. Puede pensarse entonces que el campo
gravitatorio no es un campo de fuerzas (por unidad de masa), sino un campo de aceleraciones.
En la definición de aceleración intervienen tanto el espacio como el tiempo, de forma que el
campo gravitatorio es capaz de distorsionar estas dos magnitudes. Esta es la concepción de
partida de la Relatividad General, teoría de la gravitación que corrige algunos aspectos de la
teoría de la gravitación debida a Newton.
Para el caso específico de una masa sobre la Tierra, podemos reescribir la ecuación previa
usando ME como la masa de la tierra ms, la masa del objeto y Re la masa entre los centros de
los objetos (o la distancia al centro de la Tierra).

La constante de proporcionalidad entre la fuerza que ejerce la Tierra sobre una masa y esa
masa es la aceleración de la gravedad.
Para el caso de una esfera uniforme podemos considerar que toda la masa terrestre se ubica
en un punto.

En una tierra simétricamente esférica el valor de g es constante a través de su superficie.

2
La distribución de masas en su interior podría ser variable (una tierra hueca, una tierra con
núcleo hueco, una tierra de densidad homogénea, heterogénea pero concéntrica.

Sin embargo la tierra es asimétricamente heterogénea.

Un relevamiento gravimétrico detecta diferencias de densidad en rocas del subsuelo, por ej


intrusiones ígneas y fallas.

3
Un gabro tiene contraste positivo con rocas sedimentarias (es más denso), pero un domo
salino tiene un contraste negativo.
Puede haber litologías distintas pero que tengan la misma densidad, entonces no hay
anomalía.
Anomalía gravimétrica: resulta de la diferencia en la densidad entre un cuerpo y las rocas del
medio que lo rodean (contraste de densidad).

Las fórmulas previas solo valen apara una tierra irrotacioanl. Al rotar se genera una aceleración
centrífuga que genera achatamiento polar (forma elipsoidal de la tierra)

La gravedad varía con la latitud debido a dos componentes:

La tierra no es esférica debido a que las fuerzas centrífugas han producido abultamiento
ecuatorial. Un punto en el ecuador está relativamente más lejos del centro de la tierra por lo
que la gravedad allí es menor. Sin embargo la amplitud de este efecto es menor a lo esperable
ya que sobre el ecuador la columna de masa terrestre infrayacente es mayor que en los polos
(menos atracción gravitatoria).

Para calcular la aceleración centrífuga tengo que tener en cuenta el ángulo (la latitud).

4
Esta interacción hace que la tierra en rotación se aparte de la forma esférica, tomando la
forma de un esferoide. La tierra, en primera aproximación, es bien definida por un elipsoide
con abultamiento ecuatorial y aplastamiento polar.

Sumamos otro término, tiene en cuenta que la forma de la tierra no es una esfera.

La velocidad angular en un punto sobre la superficie de la tierra decrece desde un máximo en


el ecuador a cero en los polos.

La gravedad de los polos excede a la del ecuador en -5,186 mGal.


La gravedad de la tierra es la suma de la atracción gravitatoria + aceleración centrífuga
´corrección por achatamiento. De esta manera la g solo depende de la latitud.

La gravedad varía con la latitud debido a dos componentes:

1. La tierra no es esférica debido a que las fuerzas centrífugas han producido un


abultamiento ecuatorial. Un punto en el ecuador está relativamente más lejos del
centro de la tierra por lo que la gravedad allí es menor. Sin embargo la amplitud de
este efecto es menor a los esperable ya que sobre el ecuador la columna de masa
terrestre infrayeciente es mayor que en los polos (menos atracción gravitatoria).
2. La velocidad angular en un punto sobre la superficie de la tierra decrece desde un
máximo en el ecuador a cero en los polos.

Fórmula de newton original, corregida por el achatamiento y por la aceleración centrífuga:

5
Es la gravedad teórica, no tiene en cuenta el contraste de densidades.

Elipsoide
Es la superficie matemática que más se aproxima a la forma de la tierra. Calculado para una
tierra con densidad homogénea.

Primera evidencia de que la tierra es redonda: Eratóstenes (275-195 B.C.)

Primera hipótesis de que la tierra está achatada en los polos: Newton

Primeras mediciones del achatamiento terrestre: Clairaut (1736) y Bouguer (1743) • La forma
de la Tierra puede representarse matemáticamente mediante un elipsoide: – Radio Ecuatorial
= a – Radio Polar = c – Achatamiento f = (a-c)/a

Ecuación del achatamiento f = (a-c)/a

Geoide
El ​elipsoide ​es una superficie matemática y el ​geoide ​es una superficie física. El elipsoide
representa la forma de la tierra, pero el geoide una superficie equipotencial, donde el campo
de gravedad tiene potencia constante, se ven todos las inhomogeneidades que tiene la tierra.
El elipsoide ​no tiene en cuenta las variaciones internas de la tierra. El geoide es medido, a
diferencia del elipsoide.
El ​geoide ​es una superficie que representa al nivel medio del mar, la cual se prolonga por
debajo de los continentes y cubre a la tierra en su totalidad. No representa la topografía.
Puede ser imaginada como la superficie del mar en condiciones ideales de quietud y es un todo
punto perpendicular a la línea de plomada o dirección de gravedad. El geoide es un modelos
físico que corresponde a una superficie equipotencial del campo de gravedad (potencial
constante).
Sobre los ​océanos ​el geoide coincide con la superficie del
mar (asumiendo la inexistencia de corrientes oceánicas,
olas, etc).

Sobre los ​continentes​: el geoide no es la superficie


topográfica. Se prolonga por debajo de los continentes.
Ondulaciones del geoide​: son el resultado de la inhomogeneidad de la tierra (distribución de
masas).

Ondulaciones del geoide = ​diferencias, en metros, entre el geoide y el elipsoide de referencia​:


– Exceso de masa ⇒ elevación del geoide (geoide = sup equipotencial ⇒ debe
permanecer perpendicular a la dirección del campo gravitatorio)

– Déficit de masa⇒ depresión del geoide

6
Si la tierra fuera perfectamente homogénea, el geoide coincidiría con el elipsoide.

Elipsoide

La superficie de la tierra es matemáticamente


representada como un elipsoide.
Eje principal coincide con un radio ecuatorial (a).

Eje menor coincide con radio polar (c )

Achatamiento polar es definido como f = (a-c) / a

Se usa el elipsoide que mejor se aproxime


aunque sea localmente al GEOIDE
La utilidad principal del geoide es establecer la
superficie de referencia de la altura ortométrica,
conocida también como la altura sobre el nivel
medio del mar y se aplica en trabajos de ingeniería topográfica, cartográfica, apoyo terrestre
para fotografía aérea, y como insumo para la generación de modelos digitales de elevación.
Combinando información de un modelo de alturas geoidales con alturas geodésicas obtenidas
mediante técnicas de posicionamiento satelital, es posible obtener alturas ortométricas de
cualquier punto sobre el terreno.
La manera de transformar el valor de la altura ​geodésica ​(h) que proporciona un receptor GPS
en un valor de altura ​ortométrica ​(H), es mediante la resta del valor de la altura geoidal ​(N)
dada por un modelo digital de elevación geoidal. Puedo transformar una altura en otra si
conozco la ondulación del geoide. Las diferencias de altura generan variaciones muy grandes
(más grandes que las variaciones de latitud).

H=h-N

7
Un datum es un elipsoide de referencia que es definido y orientado de tal manera de lograr la
mejor aproximación al geoide tanto para su totalidad como para un sector determinado.

La Tierra es redonda. Su achatamiento, 21 km de diferencia entre el radio polar y el ecuatorial,


es imperceptible para un observador situado en el espacio.
El geoide es la representación de una superficie equipotencial de su campo gravitatorio,
próxima al nivel medio del mar, tiene irregularidades que no superan una diferencia máxima
de 100m por encima o por debajo del nivel medio. Si estas diferencias pudieran ser hechas
visibles serían por lo menos 100 veces más difíciles de ver que el achatamiento polar.
El error de considerar a la tierra “con forma de papa” surge de no observar detenidamente la
magnitud de las exageraciones necesarias para hacer “visibles” los valores de superficies
inmateriales representadas en las ilustraciones.

8
Correcciones gravimétricas
Solo me interesa lo que pasa en el subsuelo. A los datos que medimos hay que hacerle
correcciones. Le tengo que restar distintas variaciones (topografía, latitud etc). Quiero
relacionar la gravedad con las estructuras geológicas en el subsuelo.

Le tengo que restar efectos para quedarnos sólo con lo que me interesa.

Usamos como referencia al elipsoide o datum. Tengo que tener bien medida la altura, porque
es lo que más hay que corregir. El elipsoide es el nivel de referencia para corregir la altura.

La ​gravedad observada es la que se mide en el estación, con un gravímetro. Incluyen el efecto


de todas las masas internas en el subsuelo (contraste de densidad), de la rotación de la tierra
(latitud), variaciones de altura. Pero solo me interesa lo del subsuelo (las anomalías).
Entonces la g​obs se mide con un gravímetro. Incluye efectos que responden a variaciones de
altura, rotación de la tierra, variaciones de la latitud. Pero solo nos interesan las ​variaciones de
densidad​.

Anomalía: g​obs​ - g​teórica


La g​teorica es la que calculamos con la fórmula sobre un elipsoide, que asume que la superficie es
homogénea. Esta gravedad incluye efectos de rotación y toma en cuenta una tierra uniforme.
Tiene en cuenta los efectos de rotación. Incluye los efectos de las masas uniformes dentro de
la tierra y de su rotación, porque la forma tiene en cuenta el achatamiento y la velocidad
angular.

Al medir la gravedad y restarle la g teórica se obtiene el valor de la anomalía.

La anomalía se aparta del valor de la gravedad teórica para esa zona.

Anomalías < g​absoluta​ (del orden de los 980000mGal).


Las ​anomalías de la gravedad marcan la componente ​vertical del efecto
gravitatorio de las masas en la tierra (subsuelo) no incluidas en el
modelo teórico. (La gravedad teórica no tiene en cuenta las variaciones
de densidad). Se mide la componente vertical de la aceleración de la
gravedad porque es la que puede medir el resorte.

9
G teórica: Atracción gravitatoria elipsoide normal.

G observada: Masas anómalas, posición, altura.

Hay una masa puntual, un cuerpo más denso que lo del alrededor. Hay un vector de la
aceleración de la gravedad que apunta hacia el cuerpo, pero el resorte en ese punto en la
superficie solo puede medir la componente vertical de la gravedad.
Mientras que la gravedad teórica depende únicamente de la elipsoide, la observada es función
de más cosas (masas anómalas, posición y altura). Si quiero que me quede solo sobre las masas
anómalas, la tengo que corregir.

Entonces a la gravedad observada la corrijo y después le resto la teórica, entonces me dan las
anomalías del subsuelo.

Hay distintas correcciones


● Temporales​ (no me muevo de un punto, pero el valor de g puede ir cambiando)
o Deriva instrumental
o Efecto luni-solar (mareas)
● Estaciales
o Latitud
o Altura
o Masa interpuesta
o Topografía cercana
o Isostasia

Corrección por mareas


La gravedad medida en un punto determinado variará en el tiempo debido a los efectos
gravimétricos del sol y la luna. La tierra se deforma por las mareas.

Estos efectos gravimétricos además cambian la geometría de la tierra alterando el esferoide.

Las mareas en tierra firme son de menor amplitud que las oceánicas.
Las mareas continentales producen variaciones de elevación del orden de los pocos cm. (De
esta manera varía la distancia al centro de la tierra).

La máxima amplitud que este fenómeno produce es de 0,3 mGal. Es muy poquito, pero lo
tengo que corregir si necesito tener mucha precisión.
Estas variaciones debidas a efectos de marea son fácilmente predecibles, calculables y
descontadas. Si se sabe la ubicación, fecha y hora ya se puede saber el efecto de las mareas.

10
Depende de la atracción de la Luna y el sol. Si la luna y el sol están juntos, su efecto se suma,
aunque el efecto del sol es mucho menor que el de la Luna porque está muy lejos.

Registros para un mismo punto durante un día. Las variaciones con luna llena son más grandes:

Correcciones de deriva
Durante el día, el resorte va sufriendo una fatiga, se va a ir estirando por más que tenga el
mismo peso. El valor que voy midiendo entonces varía.
La deriva es el cambio gradual e involuntario en el valor de referencia respecto al cual se
realizan las mediciones.
Los gravímetros son aparatos termostatizados, pero aún así sufren la deformación por fatiga.
Cuanto más viejo el resorte mejor, porque la deriva cada vez es menor, llega hasta un punto
donde ya no se estira más.

Resortes poco sensibles a los cambios de temperatura.

Las derivas normales pueden ser de 0,1 a 0,5 mGal por día.

Mediciones en el mismo punto, en tiempos distintos dan resultados diferentes.

11
Hay dos maneras

1. Estación base, monitorear la variación. Requiere dos instrumentos uno fijo (base) y
otro móvil. Siempre hacia adelante, rápido. Tiene inconvenientes
a. Mayor costo
b. por su construcción, no hay dos gravímetros que tengan la misma deriva
c. elimina los efectos de marea pero no los de las derivas de cada aparato
d. calcularla para cada aparato
2. Recuperar una estación (volver periódicamente). Un instrumento.. Volver hacia atrás.
Más lento.

Entonces:

Para mareas​: se usan dos o se vuelve hacia atrás cada tanto. (ambas)
Para derivar​: solo volviendo hacia atrás. Usar dos gravímetros no me sirve porque cada uno
tiene una deriva diferente.

Tiempo entre lecturas:

Ideal: 1 hora.
Terrestre: 4/8 horas. Depende del área que yo quiera estudiar. Si es muy grande voy a poder
volver menos veces.

Aéreo: al inicio y final de cada vuelo.

Marino: cada vez que el barco vuelve al mismo puerto. Considero que la deriva fue constante.

Está superpuesto el efecto de las mareas (​rojo​) más el cambio de la deriva del resorte (​verde​).
La deriva es constante.

Consideramos que entre 2 puntos usamos una estación base. Primero medimos el valor en esa
base y en varias estaciones. Lo ideal sería volver cada hora y volver a medir en esa base.
Porque en el mismo punto, en tiempos distintos dan resultados diferentes. También se podría
usar dos gravímetros, dejo uno funcionando en la primera y voy usando el otro para las demás
estaciones. Sin embargo, no me serviría para corregir la deriva porque cada resorte funciona
diferente.

Corrección de latitud
Es espacial.

La gravedad varía con la latitud por el achatamiento y por la velocidad angular.

12
La gravedad varía con la latitud debido a dos componentes​:

1. La tierra no es esférica debido a que las fuerzas centrífugas han producido un


abultamiento ecuatorial. Un punto en el ecuador está relativamente más lejos del
centro de la tierra por lo que la gravedad allí es menor. Sin embargo la amplitud de
este efecto es menor a los esperable ya que sobre el ecuador la columna de masa
terrestre infrayeciente es mayor que en los polos (menos atracción gravitatoria).
2. La velocidad angular en un punto sobre la superficie de la tierra decrece desde un
máximo en el ecuador a cero en los polos.

La diferencia entre los polos excede a la del Ecuador en 5,186 mGal.

Para tener una precisión de 0.01 mgal, tengo que estar en una posición de 12 metros.
Si estoy calculando anomalías, y a la gravedad observada le resto la teórica, ya tengo la
corrección de latitud hecha (está incluida en la fórmula de la gravedad teórica).

Hay una fórmula que incluye términos en varias potencias pares, hasta el grado 8.

Todos los coeficientes de las fórmulas se obtienen ajustando valores medidos en tierra y por
satélites.
Para calcular la gravedad teórica, en la fórmula ya está incluido el efecto de la latitud. Entonces
ya no tengo que hacer la corrección cuando le resto la teórica a la observada.
La gravedad teórica va cambiando la fórmula para calcularla, porque cada vez se conoce mejor
el achatamiento.

En promedio, la gravedad varía 0,81 mgal/km en la dirección norte-sur. Si quiero que las
mediciones tengan una precisión con la centésima de mgal, tengo que conocer la posición en
latitud de la estación con un error menor a 12 metros.

13
Si la región en estudio no está demasiado extendida en sentido meridiano, podemos calcular el
monto de variación de la gravedad en función de la latitud media de la región.

En efecto, la variación de gravedad que buscamos es el gradiente de dirección norte.


Si trabajo en una zona no muy extensa, puedo calcular el monto de variación de la gravedad en
función de la latitud media de la región. Admitiendo que la tierra es esférica, tendremos:

14
En los polos es donde menos varía.

Corrección de altitud
La lectura que se hace más arriba es menor que la de abajo. (Inversamente proporcional al
cuadrado de la diferencia de altura).

Se hace la corrección de aire libre. La gravedad varía 0,3086 mgal/m. Con pocos metros de
altura ya hay muchas variaciones. Para seguir obteniéndose una precisión de 0.01 mgal, en
altura la precisión tiene que ser de 3 cm nada más.

Hay 3 correcciones asociadas a la elevación topográfica.

Corrección de aire libre


Es la primera de una serie de 3 correcciones asociadas a la elevación topográfica. La corrección
de aire libre (FAC) corrige la disminución de gravedad por altura, debida a una mayor distancia
al centro de la tierra.
Lleva la medición al nivel del elipsoide de referencia (uno que pase localmente por el nivel del
mar por ejemplo).

Para reducir una medición tomada a una altura h (sobre el nivel del mar) a un elipsoide de
referencia. Lleva a todas los puntos tomados a una sola altura.
Si mido algo arriba del punto P, le tengo que sumar la corrección, porque arriba del elipsoide
es menor el valor. Lo que mido es menor a lo que realmente es, entonces le tengo que sumar
la corrección.
La corrección de aire libre contempla la altura de la estación, sin embargo la atracción
gravitatoria de la masa de rocas interpuesta entre el punto de medición y el elipsoide es
ignorada.

15
Corrección de Bouguer
Esta tiene en cuenta la masa, con densidad definida.

Si bien estoy arriba y hay menos atracción, hay una masa de roca que es más atrayente. Esta
corrección considera que por debajo de la estación hay una losa de roca de altura h, y que se
extiende hacia el infinito en todas direcciones.
La corrección de Bouguer remueve aquel efecto aproximando este volumen de la roca ubicado
por debajo del punto de medición, a través de una losa horizontal infinita con un ancho
equivalente a la distancia entre la estación y el elipsoide.

Sobre la tierra la corrección de Bouguer debe ser sustraída, ya que la atracción extra de la losa
infinita necesita ser removida. Lo que calculo es mayor que lo real, entonces le tengo que
restar la corrección.

La corrección de Bouguer considera que la topografía en torno al punto de muestreo es plana,


pero raramente es este el caso.

La gravedad en A es diferente que en B, por diferencia de altura y topografía.

Va con menos porque se resta a la de aire libre.

16
Corrección topográfica
Es más compleja

Esta corrección es ​siempre positiva​.

Calcula el efecto de las masas por encima y por debajo de la superficie de medición.

Las regiones indicadas con A fueron contempladas en la corrección de Bouguer, sin embargo
las mismas no existen. Por lo tanto la corrección de Bouguer ha “restado” de más, entonces se
debe sumar para corregir ello.
La región indicada con la letra B consiste en material que no ha
sido contemplado en la corrección de Bouguer, pero que genera
un extra de atracción gravitatoria hacia arriba disminuyendo la
gravedad medida, y por lo tanto debe sumarse al valor.

La corrección topográfica se ha realizado clásicamente usando la cartilla de Hammer.


● Esta cartilla está dividida en líneas radiales y concéntricas que definen 130
compartimientos.
● La zona más externa se extiende a 22 km, distancia a partir de la cual el efecto
topográfico se considera despreciable.
● La carilla se desplaza sobre el mapa con el punto de medición en su centro. Se calcula
el promedio de elevación de cada compartimiento a los cuales se resta la elevación del
punto de medición.

17
● Se determina la atracción gravitatoria de cada compartimiento usando tablas de
referencia.
● La corrección topográfica es la suma de cada componente de cada compartimiento.

También se puede trabajar con anomalías. A la gravedad observada le resto la gravedad


teórica. ​(sin restarle la teórica es gravedad corregida, cuando se la restamos es anomalía).
Para las correcciones es mejor usar la altura geodésica (referida al elipsoide), no la altura
ortométrica (referida al geoide). Porque en la geodésica se basa la teórica.

El cálculo de la anomalía de bouguer constituye las bases de la interpretación de datos de


gravedad continentales.

La anomalía de aire libre es típicamente usada en relevamientos marinos.

altura ortométrica: referida al geoide (Estación Total)

altura geodésica: referida al elipsoide (GPS)

Para las correcciones del altura se debe usar la altura geodésica.

Método de Nettleton.
Cuando no puedo saber la densidad, la estimo haciendo mediciones de gravedad.

18
Uso un rasgo topográfico positivo, voy probando varias densidades. Calculo las anomalías para
los distintos valores.
Cada curva del gráfico es calculada con valores distintos. La mejor curva es la más recta (la del
medio), la que refleje menos la topografía. Con eso estimo la densidad.

La densidad más verdadera en este caso es la de 2,3 que


es la que menos reproduce el rasgo topográfico.

Tiene una limitación: no tiene en cuenta las raíces de las


montañas. Por eso se usa para elevaciones pequeñas, no
para cordilleras.

Isostasia
Es una corrección que no siempre se hace. No se usa tanto en prospección sino es más bien
académico
En las expediciones a los Andes, tenían que horizontalizar sus instrumentos (se usaba una
plomada para medir la vertical). Esperaban que se desviaran hacia los andes, pero la
desviación observada era menor a la esperada.

19
Entonces pensaron que debajo de las montañas había rocas menos densas (por eso había
menos atracción).
Entonces hay cadenas montañosas que tienen menos densidad que el medio que los rodea.
Toda la raíz de abajo es menos densa que lo de alrededor.
Las cordilleras están soportadas por raíces de materiales de baja densidad proporcionales a su
volumen.

El término isostasia describe el estado de la tierra en la ausencia de fuerzas dinámicas que


alteran ese estado.

Los materiales menos densos flotan sobre materiales más densos.


El equilibrio isostático significa que no hay ni exceso ni déficit de masas en una región
determinada.
Vemos que la distribución de masas es variable en el tiempo en la tierra. Es por ello que
ciertas regiones que se encuentran en equilibrio dejan de estarlo. La isostasia permite
cuantificar dichas desviaciones al estado de equilibrio.
Entonces: hay una montaña, que tiene su raíz que es menos densa (es como un iceberg, que la
parte de abajo es mucho más grande y hace que flote)

Cuanto más alta la montaña, menos densa.

Si es más densa, se hunde más en el manto.


La existencia de masa por debajo de la
superficie de menor densidad contrarrestaría
el efecto de la mayor masa aflorante en las
regiones montañosas ya sea por:
● Rocas de diferente densidad
● Rocas de igual densidad pero amyor
extensión en profundidad
● Ambos procesos
Las cordilleras están soportadas por raíces de materiales de baja densidad proporcionales a su
volumen.
La isostasia describe el estado de la tierra en la ausencia de fuerzas dinámicas que alteren ese
estado.

Los materiales menos densos flotan sobre materiales más densos.

El equilibrio isostático significa que no hay ni exceso ni déficit de masas en una región
determinada.

Distribución de masa es variable en el tiempo en la tierra

Es por ello que ciertas regiones que se encuentran en equilibrio dejan de estarlo.

La isostasia permite cuantificar dichas desviaciones al estado de equilibrio.

20
Modelo de Airy.

Los bloques corticales tienen igual densidad.

Cuanto más alta una montaña, existe una


raíz más profunda.
En este modelo hay igual densidad en todos
los bloques.
supongo que esta todo en equilibrio, no hay
ninguna fuerza externa.

Para cierta elevación, puedo calcular el


espesor de la raíz.

El peso de todas las columnas se iguala por encima de una superficie de compensación.

Las montañas tienen raíces corticales que soportan el relieve.


● Isostasia y flexión litosférica: La carga topográfica es compensada por fuerzas de
flotabilidad, lo cual resulta en variaciones laterales de ​espesor ​cortical.

Modelo de Pratt.
Las raíces de todos los bloques penetran a la misma profundidad. Los bloques más altos
poseen más baja densidad. Igual
profundidad pero distinta densidad.
O simplemente, el peso de todas las
columnas sobre el nivel de compensación
es el mismo.

La densidad varía lateralmente (debido a


variaciones en composición o

temperatura).

● Isostasia y flexión litosférica: carga topográfica es


compensada por fuerzas de flotabilidad que son
producidas por variaciones de ​densidad ​laterales en la
corteza.
Sin embargo, una carga topográfica también es resistida por la
rigidez de la litosfera.

Si la carga es lo suficientemente grande en comparación con la


rigidez de la litosfera, esta se inflexiona hacia abajo.
Ambos modelos consideran que las montañas “flotan”
en un medio más denso.

21
Condición isostática: el peso de las columnas de roca, a cierta profundidad denominada
profundidad de compensación, es igual.
Tanto Pratt como Airy tratan a la litósfera como si se comportara como una serie de columnas
conectadas. Al margen del modelo, la carga se acomoda bajo la carga localizada en la litosfera,
lo que sólo sería factible en caso de que la misma tuviera rigidez nula. En general, la litosfera
posee algún grado de rigidez lo cual produce una distribución de la carga a nivel más regional.

La anomalía isostática es la anomalía de Bouguer menos la respuesta gravimétrica calculada


para el modelo isostático elegido.

Hay cosas que pueden generar perturbaciones o apartamientos del estado de equilibrio:
● Si hay una erosión muy rápida, el orógeno ya no crece y se puede desequilibrar. A
medida que se va perdiendo altura debido a la erosión, la raíz tendría que ser menor,
pero le lleva más tiempo.
● Rebote glaciario. Cuando se funden las grandes calotas de hielo, todo ese peso que
estaba, desaparece muy rápido entonces la corteza tiene que rebotar para
acomodarse, pero le lleva mucho más tiempo, entonces tardan en recuperar el
equilibrio.

Mediciones de gravedad
¿Cómo medimos la gravedad?

Hay dos formas:


Mediciones absolutas​: Se mide directamente la aceleración de la g en un punto cualquiera.
Antiguamente se utilizaban péndulos. Ahora se trabaja con el método de caída libre. Se deja
caer un objeto, conozco la distancia y mido el tiempo. Período de un péndulo.

22
Mediciones relativas​: se usa generalmente en relevamientos. Se determina entre dos puntos
distintos cuánto cambia la aceleración. Se necesitan dos sitios de referencia donde la gravedad
absoluta es conocida. Se miden con gravímetros que trabajan con un resorte.
● Extensión de un resorte
● Levitación de una masa metálica en un campo electromagnético

Mediciones absolutas de la gravedad


Antes se usaban péndulos, deteminándose el período del péndulo.

Ahora se utiliza el método de caída libre. Son muy precisas, pero para eso se debe medir muy
bien el tiempo y las distancias. Precisión 1 microGal.
Las notables precisiones en las mediciones de tiempos (10​-12 seg) y distancias (10​-10 m) que se
logran en laboratorios de física han hecho posible determinar rápidamente g con precisiones
de +- 0.01 mGal y aún mejores.

Las distancias se miden con interferómetros de luz láser. Los tiempos se miden con relojes de
rubidio.
Algunos microsismos y vibraciones generan aceleraciones adicionales (errores). Se toman
varias lecturas y se promedian para disminuir el error.

Una masa en caída libre, en una cámara de vacío, recorre una trayectoria s en un tiempo t.

Los gravímetros absolutos son muy lentos, muy caros y muy grandes.

Mediciones relativas de gravedad


Se usan gravímetros con resortes.
No se puede saber el valor absoluto de la gravedad, solo la diferencia entre dos puntos en dos
estaciones.
La mayoría de este tipo de gravímetro se basa en una masa de prueba que pende de un
resorte o muelle. Son más rápidos que los absolutos.

La masa m está suspendida del muelle.

Al liberarse de su posición de equilibrio, la masa produce una


elongación del muelle deltaX
Esta elongación es producto de la acción resultante de la
aceleración de la gravedad (mg).
La fuerza ejercida por la gravedad se ve compensada por una
fuerza de reacción kx. (k es la constante del resorte).

Se tiene la siguiente relación:

23
Donde K es la constante del muelle. Cuando el gravímetro es trasladado de un lugar a otro,
este movimiento implica un cambio en la gravedad (delta g) lo que rendida en un cambio en el
desplazamiento del muelle (delta x), por lo tanto, podemos escribir la siguiente relación:

Los primeros gravímetros (​gravímetros estables​) medían con un


dispositivo como el que se muestra, en el cual el resorte, del cual
pende una masa, se estira diferencialmente en respuesta a un
peso que varía según cómo varía g. Sin embargo se trata de un
dispositivo primitivo, ya que el resorte es el instrumento de
medición y al mismo tiempo soporte de la masa, lo cual le quita
precisión a la medición.
El resorte sufría más deriva.

Fueron evolucionando a los


gravímetros inestables​. En estos el resorte ya no soporta el peso
de la masa. Una vara fija sostiene en un extremo la masa y la vara
se sostiene con un resorte. Ambas se adosan a una misma
superficie vertical.
La magnitud del momento que ejerce sobre la vara es proporcional
a la extensión que sufre el resorte y el seno del ángulo tita.
Tienen dos sistemas, hay otro resorte que lo que hace es
compensar la atracción gravitatoria.

El resorte tiene que tener un módulo de elasticidad elevado, coeficiente de dilatación térmica
pequeño, histéresis elástica pequeña, influencia magnética y eléctrica pequeña.
Se aíslan de variaciones de presión (vacío), temperatura (termostatización,
termocompensación), golpes, vibraciones, campos magnéticos y eléctricos.
El gravímetro debe estar alineado con la vertical por lo cual incluye un nivel. Antes de medir
hay que horizontalizarlo para medir la componente vertical correctamente

Miden con una precisión de 0,01 mGal en un rango de 5000 mGal.

Sistemas electrónicos “feedback”

Cambio en “g” produce cambio en longitud de resorte. Circuito feedback aplica CC a placas de
capacitor, genera Fz electrostática que lleva masa a posición nula. Voltaje “feedback”: medida
del valor relativo de “g”, señal digital, procesada.

24
Sistemas electrónicos​: ​hay dos placas capacitoras, se genera una fuerza electrostática que
lleva a la masa a una horizontalización perfecta.

Relevamientos terrestres
Para calcular la g absoluta necesito tener una estación
con valor de gravedad establecida.
El espacio entre estaciones debe ser coherente con la
escala de la anomalía​.

Todos los puntos del relevamiento deben tener la


posición planialtimétrica determinada (X, Y, Z). Las
cotas (Z) deberán asegurarse a los 3 metros si se
trabaja al mGal, a los 30 cm para 0,1 mGal, 3 cm para
0,01 mGal, y con error menor a 12 m para el
centésimo de mGal.

Repetir medición en base cada 2 horas ideal.

25
Gravimetría marina
Inevitablemente miden con menor precisión (alrededor de 1 mGal) debido a severas
aceleraciones verticales y horizontales.
Aceleraciones ​horizontales​, debidas al oleaje, movimiento del barco en el eje vertical, cambios
de velocidad, son en gran medida eliminados montando al gravímetro en un estabilizador
(giroscopio).

Las variaciones de gravedad por movimientos ​verticales debidos al oleaje pueden ser
canceladas promediando lecturas en un período mayor al de movimiento de las olas.

Puede cubrir áreas muy grandes de forma rápida.

Se usa el giroscopio para compensar los movimientos que sufre el barco y reducir los errores.

Gravimetría aérea
Se puede hacer con aviones o helicópteros. Pero también tiene errores mayores, si bien hay
plataformas estabilizadoras, no siempre logran compensar la aceleración.
En aerogravimetría las precisiones terrestres no se alcanzan debido a que no se puede
distinguir entre aceleración inercial (aeronave) de los cambios de g, por debajo de 1 a 10 mGal.

Estas limitaciones se deben a la resolución del método GPS, a deficiencias en las plataformas
estabilizadoras y a la falta de resolución de los gravímetros aerotransportados.

Los métodos magnetométricos aéreos tienen precisiones semejantes a los terrestres.


Hay una corrección extra que se tiene que hacer para la gravimetría marina y
aerotransportada: el gravímetro se está moviendo mientras se mueve.
El movimiento del gravímetro hace que la aceleración centrífuga y por lo tanto que el valor de
la gravedad varíe. Esta variación está relacionada linealmente con la componente E-O de la
velocidad del gravímetro. La corrección de este efecto debe su nombre al geofísico húngaro R.
Eötvös y es positiva cuando el gravímetro se mueve hacia el O y negativa cuando se mueve
hacia el E.

En helicóptero, volando a muy baja velocidad (56 km/h) y con espaciamiento de líneas de
vuelo de 50 m, se obtuvieron precisiones de 0.4 mGal y resoluciones de 300 m. No es posible
con aviones, donde las velocidades típicas son 180-260 km/h.

Relevamientos satelitales
Variaciones en el campo gravitatorio, en función de variaciones en la órbita de satélites:

La órbita de los satélites artificiales es perturbada por variaciones en el campo gravitatorio.


Por lo cual su medición precisa permite detectar variaciones temporales en el campo
gravitatorio. Se pueden obtener las anomalías gravimétricas de toda la tierra.

Determinación del campo gravitatorio desde el espacio:

● Poca precisión. Cosas muy regionales.


● Relevamiento continuo. Se pueden repetir las mediciones y ver cómo varía en el
tiempo.
● Áreas de difícil a nulo acceso como hielo polares.

26
Densidades
Densidad promedio de la corteza continental 2,7 gr/cm​3​.

La densidad no es diagnóstica del tipo de roca.


• Rocas ígneas frescas: de menor a mayor densidad: félsicas –máficas – ultramáficas, debido a
minerales más o menos densos (feldespatos y cuarzo vs. anfíboles, piroxenos y olivina).
• Rocas ígneas máficas y ultramáficas mayor rango de variación que félsicas por mayor
susceptibilidad a alteración (serpentinización y meteorización) y por vesículas y amígdalas.

• Rocas metamórficas: tienden a ser más densas que sus protolitos (ej. eclogitas granatíferas).

• Rocas sedimentarias carbonáticas más densas que siliciclásticas.


• Evaporitas muy baja densidad promedio, pero amplio rango de variación debido a variación
en su mineralogía.

• BIFs (hierro bandeado), muy densas por óxidos de hierro.

• Carbón: una de las densidades más bajas. La porosidad influye un montón en la densidad de
las rocas. Según qué material ocupe los poros, la densidad varia más (agua, aire,
hidrocarburos).

27
La densidad está dada por la sumatoria de la relación entre el volumen de los componentes
sobre el volumen total multiplicado por la porosidad.
Densidad depende de composición mineralógica, de la porosidad y del material ocupando los
poros

Cuanto más alta es la porosidad, más baja es la densidad. Sobre todo cerca de la superficie.
En roca porosa el fluido poral puede controlar la densidad. Fluidos porales tienen densidad
mucho menor que minerales.
El cambio de porosidad afecta mucho más la densidad que el cambio mineralógico de matriz.
Sobre todo en la superficie, porque hay compactación.

Cambios en densidad por compactación:


• Rocas Sedimentarias, con profundidad aumenta compactación, disminuye porosidad,
aumenta densidad

• Efecto disminuye con la presión por colapso de porosidad

• Para modelar anomalías gravimétricas: necesario tener en cuenta variaciones de densidad


con profundidad
La densidad aumenta rápido con la profundidad al principio. Al principio hay mucha
compactación, pero después la densidad queda estable.

28
Cambios de densidad con el metamorfismo

Se produce un aumento de la densidad, porque se producen minerales más densos, como el


granate. Además disminuye la porosidad.

Cuanto más grande es el grado metamórfico, mayor es la densidad.


Pero si hay serpentinización porque hay muchos fluidos disponibles, la densidad puede
disminuir.
Rocas metamórficas: en general mayor densidad que protolito por minerales de alta densidad
(ej. granate) y disminución de porosidad en protolitos sedimentarios (mármol, cuarcitas, etc.)
Serpentinización (olivina y ortopiroxeno a serpentina, brucita y magnetita): disminución de
densidad de 0.9 gr/cm3 en rocas ultramáficas

Anomalías regionales y residuales


La anomalía gravimétrica de interés puede estar superpuesta a otra de extensión regional. La
anomalía residual es de longitud más corta y corresponde a cosas más someras.

Hay una tendencia regional, se puede deber a un basamento cristalino.


Con la anomalía residual puedo ver mejor el efecto del cuerpo local. Las anomalías más
profundas tienen longitud de onda más largas, y disminuye la amplitud con la profundidad.

Promedio datos, suavizo la curva y me quedo con el regional.

Se basan en los gradientes


● Primera derivada vertical. Indica cuantos miligales varía la anomalía verticalmente.
● Derivada horizontal total. Cuantos metros varía la gravedad horizontalmente. Puede
ser a lo largo de un perfil
● Segunda derivada vertical
● Tilt
● Gradiente de Tilt

29
● Señal analítica

Anomalías regionales y residuales


Anomalía regional. Tiene una mayor
longitud de onda. (puede ser un basamento
cristalino en profundidad)
Anomalía residual (Local). Long de onda
corta. (puede ser un domo de sale en una
cuenca sedimentaria)
La anomalía gravimétrica de interés puede
estar superpuesta a otra de extensión
regional.

Para separarlas hay varios métodos:


● Estimaciones directas
● Estaminaciones gráficas. Se basa en
mapas de anomalías que cumplan
ciertos requisitos (son difíciles de encontrar)
● Estimaciones matemáticas:
○ promedios móviles
○ ajuste por funciones
○ filtrado y continuación analitica
○ modelos isostáticos y de densidades
○ gravedad de satélites.
Cuerpos profundos: mayor longitud de onda. (poca amplitud) Cuerpos superficiales: long corta.
(más amplitud).

Estimaciones matemáticas
1) Promedio móviles
Primero tengo todas las estaciones que medí. Después tomo ventanas de datos y voy
promediado todo los puntos que entran en esa ventana. (defino una ventana de 15 datos,
hago un promedio de cada una). Me queda una curva más suave, sin tantos picos.
Cuando más grande sea la ventana, más suave es la curva. Elimino las longitudes de onda
corta. Me queda una señal suavizada con longitudes de onda largas. (Me representa el campo
regional).

Entonces si al original le resto lo que calcule, me quedan las anomalías residuales. Me quedan
anomalías de longitud de onda cortas.

Este método es subjetivo. La cantidad de datos en la venta lo decide la persona.

30
2) Ajuste por funciones / aproximación por polinomios / superficies de tendencia.

Tengo que ajustar con algún ajuste a mis datos. Uso polinomios de distintos grados.

A mi gravedad medida le resto algún polinomio, según cuál vaya mejor.

Si el polinomio ajusta demasiado bien, corro el riesgo de que me quede ruido al restarselo,
porque también le restaría efectos locales (ajusta bien al regional y al local). (Me quedan solo
los errores de medición).

¿Cómo defino el grado polinómico?


Para poder determinar bien qué superficie de tendencia representa mejor, ​tengo que conocer
la geología de la zona.

También tengo que ver bien el mapa resultante, para ver si tiene o no mucho ruido.

3) Filtrado

Acá no hay que restar ninguna superficie.

El método es distinto.

Aplico filtros basados en la longitud de onda.

Aplico un filtro y me quiero quedar con ciertas longitudes de onda.

31
Los objetos poco profundos y/o de menor tamaño presentarán anomalías de longitudes de
onda menores que los objetos o estructuras regionales localizadas a grandes profundidades
y/o de gran tamaño.
Los filtros de longitudes de onda pueden ser utilizados para separar las anomalías locales a
expensas de las regionales o viceversa.

Los errores en los datos pueden provocar un aumento del ruido en las longitudes de onda
corta. Se puede eliminar mediante un filtrado.
La selección del filtro se realizará en función de la longitud de onda (o el número de onda) que
se quiere conservar. P.e.: si las longitudes de onda ,λx , menores que un valor específico, λs ,
se consideran irrelevantes, se pueden filtrar: las longitudes de onda λx (o los números de
onda mayores a 2π/ λs ) mayores que λs son pasadas sin modificación.

tengo longitudes de ondas menores y mayores. Elimino unas y me quedo con otras.

Paso bajo​: elimina las longitudes de onda cortas (altas frecuencias). Resalta las long de onda
largas (regional).
Paso alto​: elimina las longitudes de onda largas (frecuencias bajas). Resalta las long de onda
cortas (local).

Si quiero un local, le hago un filtro de paso alto.

frecuencia alta: longitud de onda corta.


También tiene una parte subjetiva, porque hay que elegir la longitud de onda a partir de la cual
hacer el filtro.

32
5) Continuación analítica o prolongaciones

Simulamos que las mediciones fueron hechas a cierta altura de la tierra. Cuando medimos
mayores distancias, nos van a dar mayores long de onda. (A mayor profundidad, mayor
longitud de onda).
Funciona como un low pass filter. Se simula un mapa como si el relevamiento de datos se
hubiera hecho a más altura por encima del terreno. De esta manera se minimizan los efectos
de las fuentes superficiales. Al alejarme tanto, pierdo de vista los detalles.
En el caso de una fuente que provoca una anomalía de forma sinusoidal, esta anomalía puede
ser extrapolada a cualquier nivel mediante la ecuación.

Siendo: A​0 la amplitud de la anomalía medida en la superficie (h=0) y k=(2π/λ) el número de


onda de la anomalía con una longitud de onda λ, siendo x la distancia en la dirección del
perfil.

Con esta ecuación puedo simular a qué distancia quiero medir (el valor de h).

Filtros de continuación:

Proyectando las anomalías observadas a un nivel inferior o superior logramos obtener el


mismo efecto que un filtro.
La ​continuación ascendente minimiza los efectos de las longitudes de onda más cortas (los
números de onda más altos), que son aquellas asociadas a fuentes locales poco profundas. La
continuación ascendente facilita la obtención de las anomalías provocadas por fuentes
profundas.

Suaviza la cosa. La longitud de onda


medida desde arriba es mayor.
Continuación descendente. ​En lugar de simular que me alejo, simulo que me acerco. Es
resolver o discriminar el número de fuentes que provocan la anomalía, desdoblar
solapamientos. Si el nivel de proyección es mayor al de profundidad del objeto de más
entidad: se obtendrá un campo con fluctuaciones extremas. La continuación descendente hay
que llevarla a cabo con extrema precaución, ya que las frecuencias altas (ruido y errores) son
amplificadas de forma exponencial. Me sirve para discriminar fuentes que son difíciles de
discriminar (si están muy juntas).

33
La prolongación ascendente separa grandes longitudes de onda. No todas ellas son debidas a
efectos profundos. Puede haberlas debidas a amplios cuerpos someros.

¿Por qué para obtener la anomalía local, restamos a la AB (anomalía de bouguer) la


prolongación ascendente en lugar de utilizar directamente la prolongación descendente?
Porque puede amplificar ruidos y además porque necesito saber antes la profundidad de esos
cuerpos que producen las anomalías.

6) Modelos isostáticos y de densidades


Ahora no se basan en longitudes de onda, sino de modelos. Acá también hay subjetividad,
porque uso modelos que asumo correctos para la zona que estoy estudiando.

Residual = anomalía de bouguer - efecto gravimétrico de la corteza inferior del modelo 3D.
Cómo obtener un modelo isostático​: tengo la topografía, y usando las fórmulas de Airy para
calcular la raíz, calculo la profundidad del moho. Tengo que asumir que está todo compensado,
tengo que asumir una densidad y también tengo que asumir que en esta zona se cumple el
modelo de Airy. Después calculo la anomalía gravimétrica que me genera ese moho. Esta
anomalía representa el regional (porque es profundo), se lo resto a mis datos y me queda una
anomalía isostática de bouguer.
Hay otra forma de hacerlo, con modelos directos. Se construye un modelo en 3 dimensiones
con distintos cuerpos. Después se toma como es la geometría del moho, y se calcula la
respuesta gravimétrica para después restársela a la gravedad observada.

34
El efecto gravimétrico de la corteza interior (moho) teórica corresponde a un modelos de
densidad 3D.

Residual = anomalía de bouguer - efecto gravimétrico de la corteza inferior del modelo 3D.

7) Gravedad de satélites.
Tienen poca resolución, porque están lejos. No detectan cosas de poco tamaño. Registran
longitudes de onda más largas.
¿En qué casos me interesaría en un relevamiento de prospección la anomalía regional?, ¿Y la
local?

Con las anomalías regionales puedo buscar hidrocarburos. Puedo detectar cuencas
sedimentarias.

Con anomalías locales puedo buscar intrusivos.


La elección de una tendencia regional implica criterio geológico, es una decisión interpretativa,
no es inocua no es única. Pueden perderse datos útiles.
La longitud de onda que para un estudio es “regional”, para otro puede ser la “residual”
buscada, dependiendo de los alcances y objetivos de relevamiento.

Resalto de bordes-estimación de profundidades


Hay distintas metodologías.
Se quiere mejorar la resolución de rasgos geológicos particulares, para resaltar los bordes y
detectarlos con mayor facilidad. Hay técnicas que ayudan a detectar estos bordes.

Resalto de anomalías
Hay varios procesamientos. Todos se basan en los gradientes (derivadas). Los gradientes me
indican cuántas unidades (mgal) varía la gravedad con cada metro de distancia. La derivada
vertical me dice cuánto varía la gravedad con cada metro en que yo me desplazo
verticalmente.

• Primera derivada vertical

• Derivada horizontal total

• Segunda derivada vertical

• Tilt (ángulo que surge de los gradientes)

• Gradiente del Tilt

• Señal analítica

Primera derivada vertical


● Enfatiza anomalías sobre los cuerpos.
● Permite una localización más precisa de las estructuras que las causan.
● Permite identificar contactos y bordes de cuerpos.
● Puede ser ruidosa, amplifica el ruido de corta longitud de onda. Así como enfatiza los
bordes, puede enfatizar el ruido. Conviene primero aplicar un filtro y eliminar las
frecuencias muy altas.

35
● Se puede calcular partir de prolongación ascendente de la anomalía a una altura z y a
una altura z+Dz. Hago la diferencia de estas dos anomalías calculadas a diferentes
alturas dividido la diferencia de altura.
● T: amplitud de la anomalía a las alturas z y z+Dz
● Primera derivada vertical (gradiente vertical):


● maximos sobre los cuerpos y minimos hacia los bordes

La azul es la que mido en el campo (puede ser cualquiera, bouguer completa, aire libre,
isostática residual).

Derivada horizontal total

La pendiente de la curva de la anomalía obtenida en el campo es mayor en donde está más


vertical el cuerpo.

La derivada horizontal la puedo calcular en la dirección x e y. La trato de calcular en la


dirección perpendicular al rumbo de las estructuras. La derivada horizontal total, calcular la
derivada a lo largo de x y de y. Se eleva al cuadrado, se suman y se hace la raíz. Saco la
hipotenusa del triángulo y esa es la total.

Gradientes: enfatizan variaciones de las anomalías en la dirección en que son calculados (x, y,
N30E, etc.) dA/dx; dA/dy

Estructuras geológicas perpendiculares a la dirección de cálculo: gradientes más altos​.

Podemos combinar gradientes horizontales usando Pitágoras: Gradiente Horizontal Total.


La Derivada Horizontal Total es perpendicular a las isolíneas (máximo gradiente alrededor del
punto de observación).

Permite determinar ubicación de fallas y contactos geológicos.


Máximos indican bordes de las fuentes. Para cada contacto hay un máximo. Da cero si el
contacto es horizontal.

Complementario a primera derivada vertical, pero permite una ubicación más exacta de fallas.

36
Se pueden calcular derivadas horizontales direccionales (no totales) para resaltar rasgos de
rumbo conocido.

Segunda derivada vertical

Enfatiza características locales, remueve el efecto de anomalías regionales.


Principal ventaja: el valor CERO sigue de manera muy precisa bordes subverticales de bloques
o fallas.
Al igual que con otras derivadas, puede ser usada para resaltar la presencia de ruido o de
errores de nivelación.
Para obtener la segunda derivada vertical, se pueden utilizar valores de la primera derivada
vertical y continuarlos hacia arriba a dos elevaciones diferentes:

Se pueden calcular derivadas verticales de cualquier orden. Ventaja: respuestas más


localizadas (más angostas). Desventaja: mayor sensibilidad a ruido de datos.

Señal analítica

También se basa en las derivadas (gradientes). Para detectar bordes.

Es la combinación de la derivada según x, según y y según z. Es siempre positivo.

Máximos sobre cuerpos discretos y sus bordes. (Porque combina la derivada horizontal y
vertical).
Ancho del máximo: puede indicar profundidad al contacto, si la señal de un único contacto
puede ser aislada y resuelta.

Forma de cresta sobre contactos verticales.

Efectiva para delinear contactos geológicos y resolver cuerpos muy próximos.


La señal analítica genera un positivo tanto para un depocentro sedimentario (cuerpo
anómalamente poco denso) como para un cuerpo anómalamente denso (valor absoluto).

37
Las cartas de señal analítica son matemáticamente análogas a las cartas de gradiente
horizontal, pero menos intuitivas, debido a que incorporan la derivada vertical. La señal
analítica la calculamos mediante la siguiente expresión. No puedo saber si tengo más o menos
densidad.

Método del tilt.

El resultado es un ángulo.

La línea roja es la señal analítica. En ángulo tilt es como si fuese la “inclinación” de la señal
analítica.

38
La aplicación de estas técnicas de resaltos de anomalías, aplicado a la carta de anomalías de
bouguer, permite mapear estructuras del basamento con ventajas distintivas sobre las
derivadas convencionales.
Resalta las medianas y cortas longitudes de ondas, producidas por inhomogeneidades en los
primeros kilómetros de la corteza ocultando las heterogeneidades del basamento, tales como
por ejemplo fallas, resaltos o discontinuidades.

Puedo ver variaciones más pequeñas en los gradientes que antes no veía.

Tilt: positivo sobre la fuente, negativo fuera de la fuente.

Tilt: la isolínea de CERO sigue de forma precisa los límites superiores de la fuente.

Tilt: copia la forma en planta de la fuente

39
40
Determinación de las profundidades de las masas que producen la anomalía
Hay 3 métodos.

•Espectro de Potencia

•Deconvolución de Euler
•Señal Analítica

Hay una ambigüedad: diferentes posibilidades pueden producir la misma anomalía


gravimétrica. Una anomalía observada puede ser explicada por una variedad de distribuciones
de masa a diferentes profundidades (aun con igual densidad).

Espectro de potencia
Se descompone en una serie de funciones senos y cosenos. Tienen distintas longitudes de
onda y distintas amplitudes. La sumatoria de n funciones de fase, longitud de onda y amplitud
distintas entre sí.

Me permite graficar el espectro de potencia.


Espectro de potencia del efecto gravimétrico de una masa puntual (m) ubicada en la posición
(x0 , z0 ): cuadrado del módulo de la transformada de Fourier de la anomalía.
Análisis de Fourier: representación de una función de una variable real mediante una serie de
senos y cosenos.

Espectro: variación de la energía o amplitud en función de la frecuencia.

Le aplico el logaritmo y obtengo una recta.

41
Puedo estimar la profundidad a partir de la pendiente de la recta que mejor ajusta.

Puedo tener el cuerpo a distintas profundidades, entonces ajusto cada tramo de la recta con
distinto ajuste.

Dos masas a distintas profundidades.

• Pendientes de las rectas de mejor ajuste a distintas porciones de la curva: profundidades

• Recta de mejor ajuste próxima a la horizontal: ruido

Deconvolución de Euler

También se basa en la longitud de onda. Tengo que usar el índice estructural.

Este método está basado en la ecuación de homogeneidad de Euler y adiciona un “índice


estructural” para producir las estimaciones de profundidad. Con esta metodología se pueden
identificar y hacer estimaciones de profundidad para una variedad de estructuras geológicas
tales como fallas, contactos, diques intrusivos, etc.

Las derivadas o gradientes en la ecuación de Euler pueden medirse, o más comúnmente,


calcularse a partir de los datos. Las únicas cantidades que no se conocen son x0, z0 y N.

Las coordenadas (x0 y z0) representan la profundidad y ubicación a lo largo del perfil de la
fuente puntual y el N representa el tipo de fuente que mejor representa la anomalía.

Es fácil verificar que modelos simples han prescrito valores de N. Hay tablas para eso.

42
Los resultados de la deconvolución de Euler consisten de varios cientos de miles de soluciones
con un amplio rango de profundidades.
El promedio de profundidades para un agrupamiento (cluster) de soluciones provee una
estimación mucho más confiable de la profundidad de la fuente.

Señal analítica

La amplitud de la señal analítica sobre un contacto 2D, ubicado en x a una profundidad h es:

Modelos 2D y 3D

Modelado en 2 dimensiones

Modelado directo interactivo


Software calcula respuesta gravimétrica de modelo geológico definido por el usuario.
Diferencias entre respuesta del modelo y gravedad observada es eliminada variando la
estructura o la densidad de los diferentes cuerpos que forman el modelo.

Magnetometría
El método magnético mide las variaciones en el campo magnético de la tierra. Se interpreta
según la magnetización de las rocas (los contrastes de magnetización).
El campo de potencial natural responsable de las anomalías magnéticas es el campo magnético
terrestre (CMT), que tiene características más complejas que el campo gravitatorio (ej, Polos
de signo opuesto, líneas de fuerza distribuidas de un modo no homogéneo en la superficie
terrestre).

Las líneas de fuerza están más concentradas en los polos. Hay menor densidad en el ecuador.

Hay dos polos, las líneas de fuerza no apuntan todas a un punto como en el caso gravitatorio.

43
El campo potencial es bipolar. Hay un máximo y un min de potencial. Van del máximo al
mínimo. Las superficies equipotenciales ahora son esferas distorsionadas.
Una carga magnética de prueba de signo positivo experimenta una fuerza de atracción o
repulsión de acuerdo con la ley de Coulomb:

Pensamos la flecha roja como un imán, formado por dos polos iguales pero opuestos. Quiero
ver cual es la atracción que el dipolo ejerce en cualquier punto del espacio. Ponemos una carga
aislada de signo positivo para per si se atrae o si se repele. Va a depender de la ecuación de
arriba.
Pongo una carga arriba, se atrae porque es positiva y se atrae hacia la negativa. Pero también
hay efecto de la carga positiva que queda abajo (hay mehyor distancia, entonces la repulsión

44
es menor que la atracción. La resultante es atracción. Abajo de todo pasa lo contrario, la que
está más cerca es la carga positiva, entonces la resultante va a ser repulsión.

A los costados, la resultante es tangencial.


Los vectores magnéticos son tangentes a las líneas de fuerza del campo magnético (figura de la
derecha).
Ahora, si en lugar de poner una carga aislada, en presencia de un campo magnético ponemos
un cuerpo, al tener masa va a sufrir magnetización inducida.

Magnetización inducida paralela a la dirección local del campo magnético inductor. En un


punto del espacio, el cuerpo magnético se magnetiza.
Al cuerpo le aparece una mag
inducida que es paralela al campo
magnético presente en ese punto.
El fenómeno se llama inducción.
Un dipolo es la unidad magnética
más pequeña. Tiene dos cargas o
masas magnéticas con signo
puesto, que generan el campo
magnético.
La expresión de Coulomb tiene
signo positivo para polos de igual
signo (repulsión), y signo negativo
para polos de distinto signo (atracción),
El cuerpo tiene unos “pequeños imanes microscópicos” situados dentro del material. En
presencia de un campo magnético, las líneas de fuerza que atraviesan el material alinean a los
vectorcitos. Este fenómeno se llama ​inducción​. La magnetización que le aparece al cuerpo es
la ​magnetización inducida​. Es paralela al campo magnético inductor.

Teoría clásica
La ​teoría clásica o “modelo Gilbert” considera que la fuente de un campo magnético es un
dipolo con dos “cargas” o “masas” magnéticas de igual intensidad y signo opuesto. Un dipolo
es la unidad más pequeña.

45
La expresión adopta signo positivo para polos de igual signo (repulsión) Y signo negativo para
polos de distinto signo (atracción)

Convención​: polo positivo es el que es atraído hacia el polo Norte terrestre.

Campo magnetizante, o intensidad de campo magnético:

El dipolo genera un campo magnético, y a este campo lo podemos cuantificar con distintas
propiedades.

Propiedades del campo magnético:


● H. ​Intensidad del campo magnético o campo magnetizante​. Es la fuerza magnética de
coulomb pero normalizada a un polo.

● M. Momento magnético del dipolo​. Apunta desde su polo sur (-) hacia el polo norte
(+). Es la flechita roja, une a las cargas negativa y positiva. Es proporcional a ala carga
de los polos y a la distancia que los separa. Cuanto mayor es la distancia entre los
polos, mayor es el momento.

Magnetización​. Momento por unidad de volumen. Normalizamos el momento porque el


momento magnético es una propiedad extensiva.

46
La subdivisión de un dipolo da origen a dipolos más pequeños, no es posible aislar los polos de
distinto signo. Esto sugiere que el campo magnético tiene un origen a nivel atómico, que no
involucra “polos” independientes.
Experiencia de Oersted. Cuando se gira un cable o en una esperia se genera un campo
magnetico.Es lo que ocurre a novel atomico, hay cargas en movimiento que giran alrededor de
sí mismos y en órbitas que genera un campo magnético, hay un momento magnético.

Campo magnético generado por un anillo. Los puntos son la entrada y salida del anillo. Se
genera un campo magnético axial al centro del anillo

Teoría electromagnética
Ley de Biot-Savart. Toda corriente eléctrica genera un campo magnético. Es un conductor
rectilíneo, las líneas de campo son concéntricas al conductor.

La atracción o repulsión magnética experimentada va a tener relación con la corriente que


circula por la espira.
las unidades de H en la teoría electromagnética son unidades de corriente sobre unidades de
distancia (A/m).

Para una espira circular de radio r, por la que circula una corriente i, el campo magnético h en
el centro de la espira es:

47
M y J también cambian en la teoría electromagnética

Dos concepciones del fenómenos magnético. Para la teoría clásica es un dipolo y para la teoría
electromagnética es un campo solenoidal.
Si tapamos la parte central son bastante parecidos, entonces para describir el fenómeno
magnético en el espacio las líneas de fuerza y equipotenciales son parecidas en ambos
modelos. Solo se diferencian en el centro.
Tengo un sólido sin campo H. Hay partículas magnéticas con sus momentos dispuestos
aleatoriamente. La suma de todos los momentos es cero. El material no tiene una
magnetización neta.

48
Si aparece un campo H (líneas rojas al costado). los momentos se alinean paralelos, ahora la
suma de todos sus momento da una magnetización (vector resultante J). (es un momento
magnético por unidad de volumen, una magnetización). Es una magnetización inducida,
paralela al campo inductor H.

En presencia de un campo ordenador, los momentos tienden a orientarse. Se produce


inducción magnética en el interior del material.

Aparece en el cuerpo una magnetización J, paralela al campo ordenador. Ésta a su vez, genera
un campo magnético.

Dentro del material entonces está H y J.


La suma de las dos cosas constituye el flujo magnético o inducción magnética. Es todo lo que
pasa dentro del material que incluye al campo magnetizante (H) y la magnetización (J)

49
Dentro del cuerpo, las líneas de fuerza van para arriba, al igual que el campo inductor (lo
refuerza). Por eso hay más flujo dentro del cuerpo. El campo magnético tienden a fluir por ese
cuerpo (como si el cuerpo fuera más permeable magnéticamente)
Si una sustancia se somete a la acción de un campo magnético externo H, se originará en ella
un magnetismo inducido (J), paralelo (o antiparalelo) al campo inductor. En el caso del que el
cuerpo se magnetiza anti paralelamente, dentro del cuerpo el flujo va a ser menor porque se
magnetiza en dirección contraria al campo magnetizante.
Se establece en la sustancia un flujo magnético, que puede ser de mayor, menor o igual
intensidad que el flujo exterior, dependiendo de la magnitud y dirección de J.

La inducción magnética B se cuantifica como la densidad de flujo magnético.

50
La inducción magnética B en el interior de un material estará compuesta por el campo externo
aplicado H y por la magnetización J. Se multiplica por la constante de proporcionalidad
(permeabilidad magnética).

La intensidad de flujo es mayo cuanto más permeable sea el material.

En el sistema clásico, H y B son intercambiables en el vacío, pero no son lo mismo.

51
En los dos casos de abajo, el flujo magnético dentro del material es mayor que el de afuera.
Para que B sea menor que H la permeabilidad magnética tiene que ser menor que 1.
Susceptibilidad magnética: facilidad con la que se magnetiza un material en presencia de un
campo magnético H. (La permeabilidad es cuanto facilita el flujo de un cierto campo H,
mientras que la susceptibilidad es con qué facilidad se magnetiza).

La susceptibilidad volumétrica es adimensional tanto en sistema cgs como SI, pero sus
magnitudes difieren en 4π

Hay una relación entre susceptibilidad y permeabilidad


La inducción magnética B en el interior de un material estará compuesta por el campo externo
aplicado H, y por la magnetización J.

52
En el primer caso el flujo es menor adentro que afuera porque la permeabilidad es menor que
uno. La susceptibilidad tiene que ser también negativa. Si k es negativa, la dirección de J con
respecto a H es antiparalela. (La magnetización es opuesta al campo aplicado).

El campo magnético puede ser descripto mediante el modelo de polos (gilbert), o el de bucles
de corriente eléctrica (ampere). El segundo es más realista, pero el primero tiene la virtud de
ser más sencillo matemáticamente.

Las unidades de los sistemas cgs y SI se basan en los modelos de GIlbert y Ampere,
respectivamente. Las unidades cgs son más intuitivas y evitan el conflicto en el tratamiento
B/H, por lo que todavía no fueron abandonadas en el ámbito académico dedicado al
magnetismo.

53
El modelo de polos permite una aproximación intuitiva sencilla al problema de determinar
cómo será la respuesta de un cuerpo (geológico) en presencia de un campo magnetizante.

Modelo de polos para calcular la respuesta de cuerpos sencillos:

54
La respuesta del dipolo de abajo está atenuada por la mayor distancia. Genera una respuesta
de menor amplitud y mayor longitud de onda.

- La respuesta observada en superficie, se debe a la contribución sumada de dos polos, y no de


uno. Se atenúa con la distancia como 1/r3. Recordar q gg ue la aceleración gravitatoria g de un
cuerpo de masa concentrada en su centro variaba como 1/r2 - Por esta razón, suponiendo una

55
fuente de anomalia gravitatoria que también genere anomalía magnética, la anomalía
magnética tendrá una forma más apretada que la gravimétrica, más abierta.

Se anulan las cargas y quedan solo las de los bordes.


Los momentos magnéticos dipolares del material tienden a alinearse con el campo inductor H.
La integración de esos momentos atómicos da lugar a la magnetización neta del material J
magnetización neta del material, J.
Utilizando el modelo de polos (izquierda), podemos simplificar el cálculo de J, ya que los polos
N-S de dipolos sucesivos se anulan entre sí, y J es debida solamente al efecto de los polos no
apareados, en las caras libres.

Los vectores son proporcionales a la distancia que los separa del cuerpo y tangentes a las
líneas de fuerza. Estos vectores son la respuesta del cuerpo. En cada punto medimos el vector
rojo más el violeta.

56
ΔF es el vector violeta y F es el rojo. Mido la suma de los dos F​IGRF

Los magnetómetros de uso corriente en prospección


están diseñados para medir el módulo del campo
resultante.
Por lo tanto, se determina la componente del campo
anómalo en la dirección del campo dirección del campo
ambiente (siempre y cuando F​IGRF​ >> ΔF).
Es decir, se determina en cada punto si F tiene un
módulo mayor o menor que el esperado de no existir
anomalías.

El magnetómetro mide si es mayor o menor que el de la tierra. En el primer caso, como el


vector violeta ayuda un poquito al rojo, se grafica el punto a penas arriba del cero. Más hacia la
derecha, cada vez el vector violeta es más importante y “ayuda” más al rojo, entonces hay
maysos intensidad y el punto rojo se grafica más arriba. Cuando los vectores se contrarrestan,
el valor de la intensidad se vuelve negativo.

57
58
Si el cuerpo no es puntual, su respuesta se calcula integrando la de infinitos dipolos. como los
polos de dipolos sucesivos se cancelan mutuamente, sólo tendrán efecto los polos libres en los
límites del cuerpo.

59
Campo magnético local del cuerpo depende de:

● Forma y tamaño del cuerpo, profundidad


● Orientación, posición con respecto al CMT
● Susceptibilidad magnética
La forma de la anomalía depende de cómo interactúa el CMT principal con el campo local del
cuerpo.
Valen generalizaciones sobre relación profundidad/longitud de onda, contraste/amplitud, pero
la forma de la anomalía es compleja debido a los múltiples factores que la controlan.
Ambigüedad!

Para describir anomalías:


● simetría
● ubicación del pico de anomalía con respecto al cuero
● relación de las amplitudes entre máximos y mínimos
● variaciones de amplitud y/o longitud de onda al cambiar parámetros

Campo magnético de la tierra


El campo magnético de la tierra no es como un dipolo centrado en el centro de la tierra.
No es homogéneo en la superficie de la tierra ni constante en el tiempo. Por eso lo
necesitamos conocer con detalle.

En presencia del CMT, los cuerpos magnetizables experimentarán inducción. Es decir,


adquirirán una magnetización paralela a la del CMT.
En prospección geofísica, la señal es la magnetización que puedan adquirir cuerpos geológicos,
en presencial del CMT.

60
Campo actual:

Se asimila a un dipolo centrado en el centro de la tierra. Es casi paralelo al eje de rotación de la


tierra. En la superficie, el campo magnético tiene intensidades variables, depende de la
intensidad del flujo. En los polos el flujo magnético es mayor.

En lugares donde hay más líneas de fuerza, hay mayor densidad de flujo magnético (polos).
Tanto B, como H y J son vectores. Se caracterizan por su dirección y sentido, además de su
módulo. Podemos usar x y z o el ángulo de la declinación (D) y la inclinación (I)
Vector campo magnético. Los podemos caracterizar angularmente con el ángulo de la
declinación y con el ángulo de la inclinación.

61
En el hemisferio norte el campo
apunta hacia abajo, la inclinación es positiva.
Ángulo de declinación. se mida en la horizontal. Es el ángulo que forma la proyección
horizontal con el norte.
Ángulo de inclinación. se mide en el plano vertical que contiene el vector de magnetización. Si
el vector va hacia arriba es negativo (lo que pasa en el hemisferio sur).

Fuentes de información sobre el campo magnético terrestre

Variaciones en espacio y tiempo registradas por:

62
● Datos en observatorios geomagnéticos
● Datos en relevamientos terrestres, con magnetómetro móvil.
● Datos de satélites
● Análisis matemático de las observaciones
● Magnetismo remanente de las rocas

Campo geomagnético de referencia internacional (IGRF)


Es una descripción matemática del campo magnético terrestre a nivel global, basada en los
datos medidos a través de las diversas fuentes.

El modelo se va revisando cada 5 años.


Gauss. Desarrolló el análisis matemático de los
armónicos esféricos.

Toda señal o función se puede representar como la


superposición de ondas básicas (armónicos).
La curva de mi respuesta la descompongo y la
represento como la suma de sinusoides. A medida que
voy para arriba, hay menor longitud de onda. La suma
de todas, me da la respuesta original. Cada una tiene
una amplitud, fase y longitud diferente.
Los armónicos esféricos es eso pero aplicado a una
esfera.

Existe una combinación apropiada de amplitud y fase


para cada sinusoide, cuya duma reproduce la función objeto de estudio.

Para el caso del campo magnético. El grado más sencillo es el dipolo. Es el grado 1.

Hay variaciones más pequeñas también: cuadripolo, octupolo y así.


Cada grado puede tener varios órdenes, puede haber distintas orientaciones. A medida que
aumenta el grado, aumenta el orden. Cada vez hay “arruguitas” más pequeñas.

A mayor grado del armónico, más compleja es la función.

63
El origen puede ser externo o interno:

r es la distancia al centro de la tierra. Entonces para la parte de la ecuación que es de origen


interno, cuando más lejos estoy del centro de la tierra (r está como divisor), menor sería el
potencial magnético. Mientras que para el campo externo, a mayor r, mayor es el campo
magnético.

El análisis de los valores obtenidos empíricamente para los coeficientes de gauss muestra que
en los armónicos de menor grado, la contribución externa despreciable. Para los armónicos de
menor grado, el campo magnético tiene origen dentro de la tierra. Es decir, se demuestra
matemáticamente que el origen del CMT principal está dentro de la tierra y no es en espacio
exterior.

64
El dipolo me indica fuentes más profundas que los grados más altos. (A mayor profundidad,
mayor longitud de onda).

Los grados 2 a 13 vienen de la parte más superficial del núcleo y del 13 para arriba son fuentes
corticales.

La forma del campo magnético se asemeja a la del dipolo, porque domina.

Tenemos el campo dipolar y el campo no dipolar (es el total menos el dipolo).

65
El CMT se describe matemáticamente a través de análisis de armónicos esféricos, cuyos
coeficientes se establecen a partir de datos empíricos (de observaciones, satélites, etc.) El
resultado es el IGRF, calculado con los primeros 13 grados de armónicos.
El cálculo de los coeficientes de los distintos armónicos permitió determinar que el CMT de
origen externo es casi despreciable, especialmente para los armónicos de menor grado.

El armónico de grado 1 (dipolo) es dominante.

El polo geomagnético es solo del dipolo (es el punto de intersección del dipolo con la superficie
terrestre), el magnético incluye todo los demás órdenes.

El polo magnético no coincide con el magnético, por efecto del campo no dipolar.

66
Variaciones temporales del CMT.
El examen del modo en que varían los coeficientes del IGRF a lo largo de lo saños permite
separar dos rangos de variación:

1. De origen interno: en general lentes y continuas a lo largo de los años (en períodos de
polaridad estable) variaciones seculares.
2. De origen externo: de plazo más corto (usualmente menor a 5 años). Variaciones muy
rápidas. Ocurren durante el día.

Variaciones seculares:
La componente no dipolar del CMT no está fija, tiene una variación secular.

Se puede ver como fue migrando la línea agónica (línea de desviación cero) a lo largo de los
años.

La intensidad del dipolo también ha ido decayendo a lo largo de los años, y sigue decayendo.
Deriva del polo geomagnético (bamboleo del dipolo) y del polo magnético. No varió lo mismo
en el norte que en el sur.

Entonces varían tanto el campo dipolar como el campo no dipolar.

67
En relevamientos magnetométricos, la variación secular (lenta) no será observable.

Por otro lado, el relevamiento magnetométrico sí será condicionado por la variación regional
del CMT.
La configuración regional del CMT no es estática, sino que varía a lo largo del tiempo (variación
secular).
Por lo tanto, la variación secular nos obliga a verificar la configuración regional del CMT vigente
para el momento en que se realiza el relevamiento.

Variaciones de origen externo


Sí influyen en el relevamiento que hago en el campo, la tengo que conocer y corregir.

Está relacionada con la actividad solar, cambia a lo largo del día.

La emisión solar de radiación electromagnética es acompañada por gas ionizado (plasma). La


magnetopausa es el paraboloide en el que el viento solar se equilibra con la presión magnética
de la magnetosfera terrestre. Dentro de la magnetosfera se generan fajas de radiación de
partículas.
En el límite inferior de la magnetosfera se sitúa la ionosfera, entre 50 y 1500 km de la
superficie terrestre. Capa ionizada por la radiación solar, de la que dependen las
telecomunicaciones.
El viento solar, plasma, llega a la magnetosfera y alcanza la ionosfera. Dentro de la
magnetosfera hay partículas cargadas.

68
La ionosfera es donde los iones consiguen concentrarse, interactuar entre sí y formar nubes de
electrones que generan corriente horizontal en anillos. Las moléculas ionizadas liberan nubes
de electrones que forman corrientes horizontales en anillo, más intensas en el lado diurno,
fuente de las variaciones de origen externo del CMT
Entonces la ionosfera está cargada por la radiación solar. La máxima carga ocurre cuando está
solo directo. En la parte nocturna, las corrientes cesan.

Si ponemos un magnetómetro fijo a lo largo del día, vamos a ver como va cambiando la
intensidad del CMT, por la migración de corrientes horizontales en la ionosfera.
Se genera una variación diurna. Genera variaciones de intensidad de una decena de nT. En el
Ecuador pasa lo opuesto a lo que pasa en las latitudes medias.
Tormentas magnéticas​. Ocurren esporádicamente. Son muy rápidas y muy cambiantes. No se
puede registrar, por ende no se puede corregir. No es conveniente hacer relevamientos
cuando hay tormentas magnéticas.

69
El modelo que explique el origen del CMT, debe tener en cuenta los siguientes datos:

● Su configuración espacial en la superficie terrestre, coincide aproximadamente con la


de un dipolo (predominante), situado en el interior de la Tierra y alineado
aproximadamente con su eje de rotación.
● Tiene una variación temporal.
● Experimenta cambios de polaridad sin periodicidad fija.
● Parece haber existido desde el Arqueano (magnetismo fósil en rocas de todas las
edades)
El núcleo externo es fundido, es el material conductor; está en movimiento, hay conducción .
Genera la realimentación constente del campo.

En el núcleo externo líquido hay convección vigorosa, con velocidad del orden de 10 km/año
(medible a escala de siglo) (Comparar con la convección en el manto, donde los ciclos de
convección duran 50-200 m.a -según el tamaño de la celda-)

Modelo de geodínamo
La rotación mecánica se convierte en energía, por el movimiento de un cuerpo conductor en
un campo magnético estático, al modo de una dínamo.

Núcleo externo fluido compuesto de hierro (con Ni y otros elementos livianos), posiblemente a
más de 3000°C. A esas temperaturas aún es un buen conductor de electricidad.

Convección térmica y composicional (los min más pesados van bajando) en el núcleo externo.

70
Además la tierra gira, entonces las celdas
convectivas se organizan por efecto coriolis.
Las líneas de fuerza del campo magnético
creado por el conductor en movimiento, se
mueven junto con el conductor.

El efecto coriolis organiza el flujo helicoidal en


columnas tangenciales.
La reconexión de las líneas de fuerza
magnéticas retorcidas genera un campo
alineado con el eje de rotación terrestre. Son
arrastradas por el material conductor (lineas
azules).

En la interfase manto-núcleo se producen turbulencias, que serían responsables de la


componente no dipolar y sus variaciones temporales.

Las reversiones de polaridad podrían estar relacionadas con períodos de gran inestabilidad en
la estructura del CMT (debilitamiento de la componente dipolar, predominio de la
componente no dipolar).

Nos interesan las variaciones en el CMT que se puedan registrar en el relevamiento geofísico:
● Variación regional debida a la inhomogeneidad del campo principal (verificar con IGRF
de la época correspondiente)
● Variación temporal rápida del campo magnético terrestre (ocurrida durante el
relevamiento)
● Variación local debida a la presencia de cuerpos magnetizados en el subsuelo (para
aislar esta, se deben calcular y remover las otra fuentes de variación).

De las variaciones temporales, sólo influirá la variación diurna, que se debe corregir, para lo
cual es necesario registrarla (con un segundo magnetómetro fijo, volviendo repetidamente a
estación base, o recurriendo a datos de observatorio situado a menos de 150 km).

71
Las variaciones regionales del CMT que no sean debidas a fuentes corticales, están
contempladas en el IGRF. El cálculo del IGRF para cada punto de la superficie terrestre es
posible en numerosos sitios de acceso público.

Materiales magnéticos naturales


La inducción magnética B en el interior de un material estará compuesta por el campo externo
aplicado H, y por la magnetización J.

La roca internamente tiene momentos magnéticos. Cada uno tiene que ver con el
desplazamiento de cargas generado por el movimiento electrónico.

¿Qué hace que un material tenga más o menos susceptibilidad magnética?


La causa última del magnetismo de los materiales es que en ellos existen cargas en
movimiento: los electrones.

72
Los distintos arreglos electrónicos dentro de la materia condicionan sus propiedades
magnéticas.
Para cada electrón existe un momento magnético. Los átomos tienen momentos magnéticos
atómicos debido al movimiento orbital y al movimiento de spin de sus electrones.

Los momentos están cuantizados. La unidad menor de momento magnético es el magnetón de


Bohr.

Diamagnetismo
Si tengo un orbital completo, el momento magnético neto es cero. Un orbital completo no
tiene magnetización (en ausencia de campo externo). Ahora si aplico un campo externo, en la
misma dirección de uno de los momentos, se genera una fuerza de lorentz. La dirección de la
fuerza de lorentz se ve con la regla de la mano derecha.
Entonces cuento yo aplico un campo externo, se genera una fuerza
de lorentz que frena un poco el electrón con el momento en esa
dirección, y acelera al electrón con el momento desfavorable. El
momento magnético en contra se hace más grande. Ahora aparece
un momento magnético neto, ​opuesto al campo magnético
aplicado. Esto se llama ​diamagnetismo​. Esto pasa en materiales
que tienen los orbitales completos.
Todos los materiales experimentan diamagnetismo en presencia de
un campo magnético, todos tienen algún orbital completa.

En estas sustancias, la susceptibilidad magnética es negativa.

Podemos calcular cual es la susceptibilidad estándar con la formulita:

N: número de átomos por unidad de volumen; Z: número de electrones del átomo; r radio
medio de la órbita electrónica; e y m carga y masa del electrón.

Cualquier material que esté compuesto exclusivamente por átomos diamagnéticos, tendrá una
susceptibilidad aproximadamente constante (por unidad de volumen) de:

-14 x 10 -6​​ SI

El porcentaje en volumen de minerales diamagnéticos en la roca permitirá aproximar la


contribución diamagnética en la κ de una roca. Sin embargo esta contribución es fácilmente
enmascarada en presencia de otros efectos más fuertes.

73
Cualquier mineral que no tenga Ni, Co Fe y Mn, va a ser diamagnético.

No varía con la temperatura.

Paramagnetismo
Es propio de materiales que contienen átomos con un momento magnético neto, por la
presencia de electrones desapareados. Hay orbitales incompletos. Si estos momentos
magnéticos tienen escasa interacción mutua, ocurre paramagnetismo.

En este caso la magnetización se genera en el mismo sentido que el campo aplicado.

Los metales de transición tienen electrones desapareados. Fe, Mn, Co


y Ni tienen un momento magnético neto. Algunos elementos tienen
electrones desapareados pero se van en las uniones (estado de
valencia).

Energía térmica​. Es un elemento aleatorizador. Provoca agitación


térmica en la estructura cristalina.
Agitación térmica: movimiento aleatorio de átomos y moléculas en
cualquier objeto que está a temperaturas por encima de cero
absoluto.

Debido a la energía térmica, las direcciones de los momentos magnéticos oscilan rápida y
aleatoriamente.

74
Entonces en los materiales paramagnéticos hay momentos netos en las estructuras de los
minerales, si aplico un campo externo se alinean pero por otro lado está la energía térmica que
lo intenta hacer aleatorio.
Debido a la competencia con la omnipresente energía térmica, la capacidad de alineamiento
de los momento magnéticos es muy baja.

Todo esto se aplica si los momentos magnéticos dentro del átomo están relativamente
separados.

75
Ferromagnetismo
Si dentro de la estructura cristalina del Fe situó a los átomos lo suficientemente cerca como
para que interactúen entre sí cambia la cosa.

Si los átomos están muy juntos, en algún punto van a compartir la obrita. Se genera una
magnetización espontanea, por la interacción de permuta.

No es una propiedad del átomo sino de cómo es la estructura cristalina.


El ferromagnetismo es causado por los mismos iones paramagnéticos, organizados en la
estructura cristalina de modo tal, que experimentan permuta

76
Los minerales con estructuras más apretadas van a ser más ferromagnéticos, como óxidos y
sulfuros.

Interacción de permuta. Puede aparecer entre dos átomos:

Superintercambio:

77
Los materiales ferromagnéticos tienen una magnetización espontánea, aún en ausencia del
campo magnético. Pero esta magnetización espontánea desaparece a cierta temperatura. Esto
es porque si aumenta la temperatura, aumenta la agitación térmica. A cierta temperatura deja
de haber interacción de permuta y el material pasa a comportarse como paramagnético.

Cada mineral tiene su temperatura crítica.


El aumento de la temperatura disminuye la posibilidad de interacción entre átomos vecinos y
por lo tanto disminuye la magnetización espontánea. El punto crítico en que la sustancia deja
de ser ferromagnética se denomina ​temperatura de Curie.
La magnetización es muy intensa, cuanto más intensa, más cargas libres. Se genera un gran
campo magnético. Es mucha energía magnetostática. Para minimizar la energía, el material
ferromagnético se divide en dominios. Dentro de cada dominio hay interacción de permuta y
magnetización espontánea.

La magnetización total tiende a cero, porque se anulan simultáneamente todos los dominios.

78
Cuando aplico un campo externo, los dominio con momentos magnéticos orientados
favorablemente al campo magnético se agrandan (las paredes se mueven) y los que están
orientados desfavorables se achican. Entonces la magnetización adquirida es paralela al campo
aplicado.
Saturación del material ferromagnético. Cuando el campo ya es muy intenso llega a un máximo
de magnetización.

En presencia de campo externo, las paredes de dominio se desplazan de modo que los
dominios orientados favorablemente aumentan de tamaño. Por lo tanto, el material adquiere
una magnetización neta paralela al campo aplicado. Al remover el campo, las paredes de
dominio vuelven a sus posiciones iniciales y la mayor parte de la magnetización neta J
desaparece.

79
Minerales ferromagnéticos más abundantes: magnetita (sensu lato), hematita, pirrotina

● Hematita presenta usualmente magnetismo débil


● Pirrotina restringida a ciertos ambientes geológicos
● Magnetita es considerada responsable de la mayoría de las anomalías magnéticas en la
corteza terrestre

La susceptibilidad magnética de la roca estará determinada por la de sus minerales


constituyentes, y sus propiedades relativas:

En ausencia de minerales paramagnéticos y ferromagnéticos, k será la de los minerales


diamagnéticos (aprox -14x10​-6​)
La presencia 5-10% en peso de Fe (en silicatos y otros minerales no ferromagnéticos) aporta
una susceptibilidad que es más de un orden mayor que la de los minerales diamagnéticos.
A su vez la presencia de 0.1% en volumen de magnetita, aporta una susceptibilidad que es más
que un orden mayor que la de los paramagnéticos.

80
Siempre y cuando hablemos de rocas frescas, el gabro será más denso que el granito debido a
su mayor contenido de minerales ferromagnesianos.
Como el primer control de la susceptibilidad magnética es la abundancia de magnetita (mineral
accesorio), la pregunta no se puede responder sin saber cuales son los minerales accesorios
presentes en cada una de las rocas.

La susceptibilidad magnética de las rocas tiene un amplio rango de variación (a diferencia de la


densidad) y es relativamente difícil de predecir sólo a partir del tipo litológico.

Remanencia magnética
Materiales ferromagnéticos:
Fuertemente magnéticos. Capaces de comportarse como imanes en determinadas
circunstancias. Ej magnetita. Se caracterizan por tener un ordenamiento magnético
(magnetización espontánea) que persiste aún en ausencia de un campo externo. Se ve en el
ciclo de histéresis.

Pierde la magnetización espontánea si la T es suficientemente alta. Temperatura de curie

Los minerales opacos, con luz reflejada se ven brillantes, porque reflejan la luz en lugar de
dejarla pasar como los minerales transparentes.
La magnetita que vemos en el microscopio está en un estado multidominio. Internamente está
dividida en dominios. Esta división ocurre microscópicamente, no es visible. Ópticamente es
homogéneo.

La magnetita que se identifica con el microscopio óptico tiene tamaños de varias decenas de
micrones, y por lo tanto se encuentra en un estado multidominio.

La magnetización espontánea está presente dentro de cada dominio, a lo largo de direcciones


energéticamente favorables. Pero las magnetizaciones de los diversos dominios se anulan
mutuamente, resultando en una magnetización neta prácticamente nula.

81
En presencia de un campo externo (como el terrestre), los dominios cuya dirección de
magnetización es más parecida a la de ese campo, aumentan de tamaño en desmedro de los
otros, generando una magnetización neta.
Esa magnetización prácticamente desaparece al remover el campo; es una magnetización
inducida.

La división ocurre para que el cristal disminuya su magnetización porque es muy grande. Pero
si el cristal es muy chiquito, la energía magnetostática es muy baja y no necesita subdividirse
en dominios. Este grano es de dominio sobre. Está uniformemente magnetizada de manera
espontánea, en una dirección que no tiene que necesariamente coincidir con el campo
aplicado, no es modificada por este.

El tamaño crítico por debajo del cual el grano no se subdivide, depende del mineral. En el caso
de la magnetita, es de alrededor de 1 micrón. Se habla de partículas de “dominio simple”.
A diferencia de lo que ocurre en los granos MD, cada grano SD tiene una magnetización
espontánea neta.
En ciertas condiciones, esta magnetización puede ser muy estable y perdurar por tiempos
geológicos, la llamamos ​magnetización remanente​. Es una dirección fósil del CMT.

La magnetización espontánea ocurrirá en


direcciones favorables energéticamente.
Cambiar la dirección de magnetización
requeriría superar barreras energéticas.
Sólo un aporte energético extra permitiría
ese cambio; mientras no ocurra ese
aporte, la magnetización permanece en la
misma dirección.
La barrera energética es consecuencia de
algún tipo de anisotropía, que hace que la
magnetización en una dirección sea más
difícil que en otra.

Hay dos causas principales de la anisotropía magnética:


● Anisotropía magnetocristalina
● Anisotropía de forma

82
Anisotroía magnetocristalina
Tiene que ver con la estructura cristalina de los minerales. Hay direcciones que son más fáciles
para magnetizarse.

Anisotropía de forma
Para la magnetita es la más importante.

Es un efecto de la auto-desmagnetización.
El grano se magnetiza hacia abajo, pero entonces adentro se genera otra magnetización. No la
anula, pero al contrarresta parcialmente.A mayor campo desmagnetizante, menor es la
magnetización.

¿Cuál es el efecto de la forma en la autodesmagnetización?

83
Es más fácil magnetizar en la dirección larga.

El factor desmagnetizante N en cierta dirección, es proporcional al porcentaje de superficie


cubierto por cargas magnéticas cuando se magnetiza en esa dirección.
Esta diferencia en energía magnetoestática es una barrera a la rotación de J. El grano se va a
magnetizar siempre en su eje mayor. En la magnetita domina la anisotropía de forma.

Cada grano tiene más de una dirección fácil:


cada partícula se magnetiza en alguna de esas direcciones fáciles. Elige la más cercana al
campo aplicado.
La suma vectorial de todas las magnetizaciones será un vector paralelo a la dirección de ese
campo (magnetización remanente natural MRN, en la dirección del paleocampo).

En el momento en que el material se formó (o se volvió ferromagnético) las condiciones


ambientales eran otras, tales que, en aquella situación, el CMT sí era capaz de influir en el
balance energético.
Gracias a eso, la dirección de magnetización espontánea de los materiales ferromagnéticos
tiene un significado geológico (es la dirección fósil del CMT).

Para un conjunto de partículas, el estado magnetizado no es un estado de equilibrio. Su estado


“natural” sería el de partículas magnetizadas en direcciones aleatorias, con magnetización neta
resultante cero.
El sistema llegará al estado de equilibrio con el tiempo, por la rotación estocástica de la
magnetización de las partículas, debida a la agitación térmica.

Consideremos un conjunto de partículas uniaxiales alineadas, magnetizadas en la dirección del


eje fácil.
La fluctuación térmica puede hacer que la dirección de cada partícula experimente una especie
de rotación browniana en contra de la barrera energética anisotrópica.
Para todo material magnetizado existe un tiempo en el cual la existe un tiempo en el cual la
magnetización remanente se habrá reducido exponencialmente

84
Relajación térmica: el tiempo de relajación depende de la altura de la barrera energética, la
cantidad de energía térmica y del mecanismo por el cual la energía térmica se transforma en
energía rotacional.
Néel (1949) demostró que el mecanismo es la deformación elástica de la partícula debido a la
agitación térmica.

Cuando más anisótropo el mineral mayor es la altura de la barrera y el tiempo de relajación es


más largo.

Superparamagnetismo
El tiempo de relajación de la remanencia depende del volumen. Si el volumen es demasiado
pequeño, la barrera energética se vuelve demasiado débil y el grano cambia espontáneamente
su dirección de magnetización (tiempos de relajación del orden de los tiempos de laboratorio).
Se alinea con el campo apenas se le aplica uno y lo pierde al toque. Es una magnetización sin
memoria.

85
Sólo las partículas cuyo tamaño esté en el rango de “dominio simple estable” son capaces de
portar una magnetización espontánea cuya dirección perdure en el tiempo (remanente).
La dirección fácil de magnetización será determinada por la anisotropía magnetocristalina o
por la anisotropía de forma.
Cada grano/cristal de material ferromagnético DOMINIO SIMPLE se magnetizará en su propia
dirección fácil. No se produce inducción paralela al campo externo.

En el caso de la magnetita en dominio simple, la dirección fácil de magnetización es controlada


por la anisotropía de forma. Es decir, la dirección de magnetización fácil suele ser la dirección
larga del cristal.
De todas sus orientaciones fáciles posibles, el grano se magnetizará en aquella que sea más
cercana al campo presente en el momento en que se forma. La magnetización neta resultante
por lo tanto será paralela a ese campo (en el caso de las rocas, paralela al CMT).
Cambiar la dirección/sentido de magnetización requiere aplicar campos magnéticos intensos,
mucho más intensos que el terrestre.
Se empieza enfriando una roca desde 600 °C. Pasa primero por la primera curva y se vuelve
ferromagnético, pero al principio acompaña a la orientación del campo (superparamagnético)
y después ya queda fija esa dirección.

86
El tiempo de relajación
también depende del
volumen de las partículas.
Cuando las partículas son
muy pequeñas, el tiempo
de relajación es muy
corto. Por ejemplo si está
cristalizando magnetita,
crecen los cristales desde
el volumen cero.

Si los granos de la roca ígnea


son meteorizados y puesto a
circular en el sistema exógeno,
al depositarse funcionan como
la aguja de una brújula.

Magnetización remanente detrítica (post - Depositacional)

Partículas magnéticas pequeñas tienen a orientarse por acción de


B con suficiente espacio poral.

87
Algunos de los mecanismos más comunes de adquisición de magnetización remanente estable
por tiempos geológicos (es decir, con tiempo de relajación tan largo, que podemos
considerarla “permanente”):
TRM ​(magnetización termorremanente): el material se enfría desde temperaturas superiores a
su temperatura de Curie, en presencia de campo magnético.

CRM (magnetización remanente química): los granos magnéticos crecen a temperatura


constante en presencia de un campo magnético.
DRM ​(magnetización remanente detrítica): los granos magnéticos se depositan en un medio
acuoso en presencia de campo magnético, y se reacomodan mientras tienen libertad de
movimiento.

En el momento en que se forman, los granos ferromagnéticos se magnetizan


espontáneamente a lo largo de su dirección fácil; cada partícula eligirá su dirección fácil más
cercana a la dirección del CMT presente en ese momento. La resultante es el m​agnetismo
remanente natural​.

Ese MRN fue adquirido por los granos ferromagnéticos, en circunstancias especiales en las que
su tiempo de relajación era corto y por lo tanto podían adaptar su dirección de magnetización
a la del campo ambiente.
El MRN adquirido así por la roca es “permanente”, la roca conserva el registro direccional del
CMT antiguo.

En una roca probablemente habrá minerales de tamaños diversos, algunos con magnetización
remanente significativa. La magnetización total de la roca estará compuesta por:

Los granos de tamaño más grueso no tienen remanente, solo inducida.

Paleomagnetismo

La posición original de las rocas fue modificada por traslación/rotación del continente/placa , y
por deformación tectónica. El estudio de la orientación de la remanencia magnética de las
rocas permite reconstruir esos movimientos: base del Paleomagnetismo.

88
Desde el punto de vista de la prospección magnética, la existencia de magnetización
remanente en las rocas es un problema, porque su dirección es difícilmente predecible. Sin
embargo, sólo las partículas ferromagnéticas de tamaño dominio simple son capaces de portar
remanencia significativa. Si la roca tiene muchas partículas de dominio simples es un
problema, si no tiene, toda la magnetización es inducida.

La magnetización total va a ser la suma vectorial

La magnetización neta del


cuerpo es la suma vectorial de
la magnetización inducida
(paralela al campo inductor) y
de la magnetización
remanente. De existir una
magnetización remanente
significativa y si ésta tiene una
dirección diferente, entonces
la magnetización neta no será
paralela al campo inductor.

Coeficiente de Königsberger
Un coeficiente Q mayor que 1 indica que la magnetización
remanente puede controlar la dirección de magnetización neta.
En ese caso, será importante conocer cuál es la dirección de J​rem

La magnetización inducida depende de CUÁNTA magnetita hay en la roca

La intensidad de la magnetización remanente depende de CÓMO es la magnetita (ej. tamaño


de grano); su dirección de CUÁNDO se formó el mineral dirección, de CUÁNDO se formó el
mineral.

89
En los basaltos domina la remanente frente a los granitos porque hay más granos de dominio
simple.

Técnicas de adquisición
Medimos la intensidad del campo, que es la suma del campo principal de la tierra más la
contribución local de la tierra (que este va a tener una magnetización paralela al campo más
una magnetización remanente).

Espaciado de muestreo
Hay variaciones locales del campo terrestre (como varía el campo). Cuando la frecuencia de
muestreo es lo suficientemente grande, alcanza para cubrir todo.

90
Tengo que elegir un especiado de muestreo que sea acorde con la variación. Si la intensidad
del muestreo es muy ruidosa, puedo usar un espaciado de muestreo muy chico y la capto toda,
pero en realidad debería quedarme con la forma general (cuando es muy ruidosa). Tendría
entonces que aplicar un filtrado de paso bajo.

Si hay señal magnética con longitud de onda más corta que el intervalo de muestreo, ésta va a
ser sub-muestreada, y se comportará como “ruido magnético”.

La longitud de onda está relacionada con la profundidad de la fuente. Por lo tanto, las
longitudes de onda más cortas corresponden a las fuentes más someras (ej. Objetos
metálicos/magnéticos en la superficie).

91
Estas fuentes someras no suelen ser las de interés.

Teniendo en cuenta que la señal magnética disminuye rápidamente, como r3, un modo de
eliminar el ruido superficial es alejar el sensor de la superficie una cierta distancia (2 metros).
Esto reduce el ruido.
Si hay una relación entre profundidad de la fuente y longitud de onda, entonces la elección del
intervalo de muestreo condiciona la mínima profundidad desde la que llega señal: todo lo que
llegue desde profundidades menores al intervalo de muestreo estará submuestrado (ruido).

Magnetómetros

Miden el módulo del vector, del campo magnético total.


● De precesión protónica (miden la intensidad de campo total)
● De bombeo óptico (vapor de Rb, Cs, K) -(miden la intensidad de campototal)
● De puerta de flujo (fluxgate) -(miden intensidad de una componente del campo)
● Gradiómetros

De puerta de flujo
Mide la intensidad del campo direccional

Se usan típicamente en laboratorios o en perforaciones

El sensor está compuesto por 2 bobinas que están


montadas e n sentidos opuestos. Cuando circula la
corriente, una se magnetiza en una dirección y otra
en otra. Entonces la magnetización total es cero
porque se compensa. Mide la componente del
campo en la dirección del eje de las bobinas.
Permite medir en forma continua.

Cuando se aplica un campo hay un desfase, porque


una se magnetiza más rápido, entonces hay una
resultante.
En presencia de un campo externo, la saturación se
alcanza más rápido en el sensor alineado con el campo que en el
opuesto. El desfasaje entre ambos genera un voltaje neto
resultante, que es proporcional a la intensidad del campo
paralelo al eje del instrumento.

Magnetómetro de precesión prótica


Miden la intensidad del campo total. Se usan en el campo.
Es un recipiente con líquido rico en protones (aguarrás, alcohol),
rodeado por un solenoide. Los protones tienen un momento
magnético. Se rodea de una bobina y se aplica un pulso
eléctrico. Se genera un pulso axial en el eje de la bobina, todos
los momentos de los protones se alinean.

Cuando el pulso desaparece, los momentos se vuelven a alinear


al campo terrestre. Pero en el acto de volver ocurre precesión.

92
La frecuencia de esto es proporcional a la intensidad campo terrestre.

Mide la frecuencia de precesión de los protones al retirar el campo polarizante. No permite la


medición continua. El sensor se orienta aproximadamente transversal al campo magnético.
Precesión: ocurre cuando los momentos magnéticos unitarios son forzados a desalinearse del
campo magnético ambiente. Cada partícula precede (rota) alrededor de la dirección del campo
en un plano de precesión perpendicular a éste.

La frecuencia de precesión es la frecuencia de Larmor, y es proporcional a la intensidad del


campo ambiente.

Efecto overhauser

Para mejorar el método. Es una modificación impuesta al magnetómetro protónico tradicional,


para conseguir mayor señal, en menos tiempo, y con menos consumo de energía.
Consiste en agregar un fluido con alto porcentaje de electrones libres. En lugar de un campo
polarizante fuerte, se aplica un campo alterno débil de radiofrecuencia. Este polariza a los
electrones, que a su vez transfieren la excitación (“polarizan”) a los protones. La polarización
es unas 5000 veces mayor, con menos consumo.

El fenómeno de precesión electrónica es demasiado débil para ser medido. Se lo puede


registrar mediante la manipulación óptica o “bombeo” de electrones para llevarlos a un estado
excitado de manera selectiva.

Magnetómetro de bombeo óptico (o de vapor de Cs, Rb o K)


En este caso se excitan electrones, no protones.
Hace uso del comportamiento del único electrón desapareado en la capa más externa de cada
átomo de los metales alcalinos en estado vapor.

Se energiza con luz del mismo elemento, circularmente polarizada (“bombeo óptico )” y
aplicando una radiofrecuencia de despolarización. Se mide la frecuencia de oscilación del
sistema, como variación de intensidad de la luz transmitida.

93
Medición prácticamente continua.

Alta sensibilidad (10 a 100 veces más sensible que un magnetómetro protónico). El sensor se
orienta oblicuamente a las líneas de campo magnético para que éstas no coincidan con las
“zonas muertas” del instrumento (no habrá bombeo óptico si el instrumento está paralelo o
perpendicular al campo).

94
Gradiómetro
Consisten en dos sensores (de cualquiera de los magnetómetros) separados por una distancia
fija. Se mide el gradiente, es decir, la variación en la intensidad de campo en función de la
distancia, en nT/m. Puede medir gradiente horizontal o vertical.
Enfatiza la respuesta de fuentes cercanas a la superficie. Si la fuente es somera la respuesta es
de gran longitud de onda y menor amplitud. Si la fuente tuviera 100 metros de profundidad,
medir a 3 metros o a 1,5 metros de distancia no va a ver diferencia.

95
Plataformas de relevamiento
Cambian los resultados según la altura en que sea relevado.

El ruido de alta frecuencia es el primero que se atenúa al subir. Además disminuye la amplitud.

En el relevamiento en helicóptero, la amplitud es menor que en el terreste.

96
Relevamiento terrestre

Espaciado entre 10 y 100 metros. Depende de la profundidad de los que queremos buscar.

Procedimiento

● El observador no debe portar objetos metálicos de ningún tipo.


● Los sitios de observación deben elegirse lo más lejos posible de vehículos, líneas
eléctricas, alambrados, puentes, etc.
● Posicionamiento de los puntos. En cada punto de medición registrar: hora, valor del
campo, condiciones meteorológicas especiales, referencias (excavaciones mineras, por
ejemplo), situación del punto con respecto a las referencias, relieve.
● Complementar con muestreo de las litologías aflorantes y mediciones de
susceptibilidad magnética.

Relevamiento aéreo

97
Hay puntos de intersección, que se relevan más de una vez. Como el campo va variando con el
tiempo, la intensidad medida no van a ser iguales. Por eso se hace la corrección diurna.
Las líneas de control tienen un espaciamiento mayor. Están separadas 10 veces las líneas de
vuelo. La altura de vuelo está entre 60 a 150 metros.

Si las líneas están muy separadas con respecto a la altura, al hacer un grillado, se va a generar
un estiramiento de los datos.

Conviene usar el de bombeo óptico, porque es continuo.

98
Hay que ir tomando la distancia entre la aeronave y la superficie, y con eso se sabe a
topografía.
Si dos rasgos geológicos están separados menor a la altura de vuelo, se verá como un rasgo
solo.

Planificación
● Hay que saber cual es el problema geológico concreto: por ejemplo, búsqueda de
depósitos minerales determinados, especificando tamaños, profundidades y
propiedades magnéticas esperadas.
● Evaluar si el objeto es detectable; contraste de susceptibilidades que exista con la roca
encajante.

99
● Estimar hasta qué profundidad/tamaño podrá detectarse el blanco.
● Circunscribir áreas de interés, recopilar información geológica regional.
● Establecer densidad de puntos, resolución necesaria, métodos de enlace de itinerarios.

En el relevamiento terrestre, lo ideal es la grilla cuadrada, pero a veces está limitado, depende
el terreno.

A la hora de hacer el mapa, se debe grillar en una grilla regular.

Métodos de relevamiento
Relevamientos terrestres​: normalmente desarrollados como perfiles o serie de perfiles. Se
puede ajustar la ubicación con alta precisión
● Útiles en áreas pequeñas con alta resolución (Minería, arqueología, ambiental)
● En áreas más grandes, como apoyo a otros estudios (ej. gravedad)
Relevamientos aéreos: utilizados mayormente en relevamientos regionales (avión) o locales
(helicóptero, drones). Serie paralela de líneas de vuelo con enlaces. Posición del avión es crítica

● Útiles en minería, relevamientos sistemáticos oficiales, exploración de HC


Relevamientos acuáticos​: mayormente como complemento de otros estudios (ej. sísmica) en
relevamientos regionales o específicos en cuerpos de agua más pequeños (variante de los
terrestres)

● Útiles en estudios académicos, ambientales, arqueológicos, exploración de HC

100
Procesamiento de los datos
Lo primero que hay que hacer es la
corrección diurna. El campo magnético va
variando a lo largo del día.
Se hace con un magnetómetro fijo en la
base.

Si no se tiene un segundo magnetómetro, se


puede contar con los datos de los
observatorios cercanos.
La máxima variación del campo magnético
ocurre en las horas del mediodía.
Un día magnéticamente activo es cuando
hay mucha variación (casi tormenta
magnética). Los datos que se toman esos
días no son tan confiables.
La corrección me sirve más cuando la
variación del campo regional es baja (si
tiengo anomalías más chiquitas)

También hay que eliminar el campo regional. El IGRF es nuestra primera opción como campo
regional a remover.

Los datos se pueden presentar en perfiles o en mapas.

101
Para hacer el mapa, tengo que tener los puntos igualmente espaciados. Hay que hacer una
interpolación y generar una grilla regular. El espaciado de la grilla tiene que tener un espaciado
relacionado al espacio real de los datos tomados.
La mayoría de los métodos de procesamiento y
presentación de datos requieren puntos
regularmente espaciados. Esto exige algún tipo de
interpolación. La interpolación se basa en el análisis
de una ventana de datos alrededor del nodo de la
grilla. Normalmente se asume que la propiedad que
se mide varía espacialmente de modo continuo.
El espaciado más adecuado para la grilla está entre
1/3 y 1/5 del espaciado de las líneas del
relevamiento.

Síntesis de los procedimientos de adquisición y reducción de datos

● Realizar las mediciones a lo largo de perfiles o grillas en el área seleccionada, en lo


posible con puntos de control.
● Editar los datos para localizarlos correctamente y para remover picos espurios.
● Corregir las mediciones por el efecto de la variación diurna (y por variación secular si la
adquisición de datos se hizo en varias etapas separadas en el tiempo).
● Sustraer el valor debido al campo principal, nivelar usando líneas de control, grillar.
● Presentar en perfiles o mapas.

Realce y presentación de datos

Filtros​: operaciones matemáticas que cambian los datos de un modo conocido. Se usan para
enfatizar características particulares, o para suprimir rasgos indeseables. Implican
transformación o convolución.

102
Objetivo: transformar respuestas a formas más simples que se relacionan más directamente
con sus fuentes y/o mejorar la resolución de rasgos particulares.
La mayor parte de ellos son similares a los aplicados en gravimetría, ya que se basan en la física
de los campos de potencial:
● Gradientes: derivadas direccionales, primera y segunda derivada, señal analítica
● Continuación ascendente y descendente
● Reducción al polo (o al ecuador)
● Pseudogravedad
● Susceptibilidad aparente
● Filtros de frecuencias, paso alto, paso bajo, paso de banda
● Suavizado: ej. ventanas móviles
● Escalamiento de amplitud: control de ganancia automática (AGC)
● Filtros direccionales

Pseudogravedad

Tanto el potencial gravitatorio como el potencial magnético dependen de 1/r. Esto permite
calcular uno a partir del otro con la ecuación de Poisson, asumiendo que la relación de
intensidad de magnetización con densidad permanece constante.
Se calcula la pseudogravedad, con el objetivo de simplificar la forma de la anomalía y centrarla
sobre la fuente.

103
A partir de la anomalía magnética
podemos calcular cuál hubiera sido la
respuesta gravimétrica. Entonces
transforma una anomalía magnética
(que puede ser compleja) en una
anomalía más sencilla y centrada sobre
el cuerpo, es decir la anomalía
gravitatoria.
Transforma la anomalía magnética en la
respuesta gravitatoria que se habría
obtenido si la fuente magnética hubiera
tenido una distribución similar de
masas. Es decir, asumiendo que la densidad está relacionada con la magnetización.

Se asume que la magnetización es inducida

La respuesta se suaviza, aumenta la longitud de onda, la respuesta se sitúa sobre la fuente.

Es útil comparar pseudogravedad con gravedad.

No todas las masas anómalas son magnéticas

Reducción al polo
EL objetivo es parecido. Es transformar una
anomalía compleja en una sencilla centrada en el
cuerpo.
Transforma a las anomalías producidas por
magnetización inclinada, en aquellas que se
hubieran producido si la magnetización hubiera
sido vertical (como en los polos).
Simplifica las anomalías y las sitúa sobre la
fuente.

Susceptibilidad aparente
Surge de asumir que la respuesta magnética se origina en un conjunto de prismas que se
extienden infinitamente en la vertical. La magnetización en cada prisma es inducida. La
profundidad al techo de los prismas es constante.
Resultados más realistas cuando puede asumirse que las fuentes magnéticas se encuentran a
cierta profundidad constante, ej. Roca cristalina cubierta por sedimentos no magnéticos, hielo,
etc.

104
Interpretación geológica
Es más difícil que en gravimetría, entre otras razones porque:

● Las anomalías magnéticas tienen formas complejas por el carácter dipolar del CMT.
● La densidad es un escalar y la magnetización es un vector
● La magnetización depende de la abundancia, tamaño y forma de los minerales
ferromagnéticos, que son accesorios en la roca.
● Cualitativa: identificación de zonas tranquilas vs. ruidosas. En mapas: identificación de
texturas, longitudes de onda, rasgos lineales.
● Cuantitativa: determinación de profundidad, forma y tamaño de fuentes.
○ Modelos 2D, 2,5D y 3D
○ Análisis espectral Análisis espectral
○ Deconvolución de Euler
La magnetización (controlada por la susceptibilidad magnética) depende de la abundancia de
un mineral accesorio (magnetita), y pequeñas variaciones en esta abundancia provocan
cambios de órdenes de magnitud en la susceptibilidad.

105
106
La “señal magnética” de los distintos tipos de roca se refleja en las texturas del mapa
magnetométrico.

Una textura ruidosa indica rocas magnéticas en superficie (izquierda).


Si las rocas magnéticas están a profundidad, su señal tendrá una longitud de onda acorde con
la profundidad a la que se encuentran. Si las rocas magnéticas están a gran profundidad, sólo
se verá señal de las rocas no magnéticas (cobertura sedimentaria) en superficie (derecha),
como una textura suave.

En la forma y distribución de las anomalías magnéticas influye:

● La orientación de los cuerpos con respecto al campo inductor


● La abundancia de magnetita en las rocas, y su variabilidad espacial
● La presencia de remanencia (muy común en volcanitas jóvenes)
● La interferencia con otras anomalías cercanas

2:05:00

107

También podría gustarte