Fundamentos de Gravimetría y Magnetometría
Fundamentos de Gravimetría y Magnetometría
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 .
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.
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.
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 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.
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.
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
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).
6
Si la tierra fuera perfectamente homogénea, el geoide coincidiría con el elipsoide.
Elipsoide
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.
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.
9
G teórica: Atracción gravitatoria elipsoide normal.
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.
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.
Las derivas normales pueden ser de 0,1 a 0,5 mGal por día.
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.
Ideal: 1 hora.
Terrestre: 4/8 horas. Depende del área que yo quiera estudiar. Si es muy grande voy a poder
volver menos veces.
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.
12
La gravedad varía con la latitud debido a dos componentes:
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.
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.
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.
16
Corrección topográfica
Es más compleja
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.
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.
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.
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 equilibrio isostático significa que no hay ni exceso ni déficit de masas en una región
determinada.
Es por ello que ciertas regiones que se encuentran en equilibrio dejan de estarlo.
20
Modelo de Airy.
El peso de todas las columnas se iguala por encima de una superficie de compensación.
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.
temperatura).
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.
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?
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
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.
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:
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
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.
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.
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.
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:
26
Densidades
Densidad promedio de la corteza continental 2,7 gr/cm3.
• Rocas metamórficas: tienden a ser más densas que sus protolitos (ej. eclogitas granatíferas).
• 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.
28
Cambios de densidad con el metamorfismo
29
● Señal analítica
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.
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.
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).
También tengo que ver bien el mapa resultante, para ver si tiene o no mucho ruido.
3) Filtrado
El método es distinto.
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).
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.
Con esta ecuación puedo simular a qué distancia quiero medir (el valor de h).
Filtros de continuación:
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.
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.
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.
• Señal analítica
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).
Gradientes: enfatizan variaciones de las anomalías en la dirección en que son calculados (x, y,
N30E, etc.) dA/dx; dA/dy
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.
Señal analítica
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.
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.
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: la isolínea de CERO sigue de forma precisa los límites superiores 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
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í.
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.
Deconvolución de Euler
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
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.
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)
El dipolo genera un campo magnético, y a este campo lo podemos cuantificar con distintas
propiedades.
● 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.
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.
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.
Aparece en el cuerpo una magnetización J, paralela al campo ordenador. Ésta a su vez, genera
un campo magnético.
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.
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).
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π
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.
54
La respuesta del dipolo de abajo está atenuada por la mayor distancia. Genera una respuesta
de menor amplitud y mayor longitud de onda.
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.
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 FIGRF
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:
60
Campo actual:
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).
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
Para el caso del campo magnético. El grado más sencillo es el dipolo. Es el grado 1.
63
El origen puede ser externo o interno:
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.
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 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.
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.
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:
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.
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.
La roca internamente tiene momentos magnéticos. Cada uno tiene que ver con el
desplazamiento de cargas generado por el movimiento electrónico.
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.
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.
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
73
Cualquier mineral que no tenga Ni, Co Fe y Mn, va a ser diamagnético.
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.
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.
76
Los minerales con estructuras más apretadas van a ser más ferromagnéticos, como óxidos y
sulfuros.
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.
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
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.
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.
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.
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.
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.
83
Es más fácil magnetizar en la dirección larga.
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.
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.
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.
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.
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:
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.
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 Jrem
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
De puerta de flujo
Mide la intensidad del campo direccional
92
La frecuencia de esto es proporcional a la intensidad campo terrestre.
Efecto overhauser
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.
96
Relevamiento terrestre
Espaciado entre 10 y 100 metros. Depende de la profundidad de los que queremos buscar.
Procedimiento
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.
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.
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
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.
También hay que eliminar el campo regional. El IGRF es nuestra primera opción como campo
regional a remover.
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.
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.
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.
2:05:00
107