estructura gravimétrica y magnética de la corteza del suroeste ...

7. MODELOS GRAVIMÉTRICOS Y MAGNÉTICOS
7.1. INTRODUCCIÓN
7. Modelos Gravimétricos y Magnéticos
El siguiente y último paso en el proceso de interpretación gravimétrica y magnética
consiste en el análisis cuantitativo de las anomalías, esto es, la modelización. Para ello, se
parte de un modelo geológico inicial en el que se incorpora toda la información disponible,
geológica, sísmica, profundidad de las superficies que presentan un contraste de densidad y/o
susceptibilidad, así como cualquier otro dato que nos ayude a eliminar indeterminaciones y a
restringir las posibles soluciones de esa distribución de densidades propuesta que, en
realidad, constituye una simplificación del modelo geológico inicial. Se calcula la anomalía
teórica que produce el modelo propuesto y se compara con la anomalía observada. Se
modifica en sucesivas iteraciones el modelo inicial propuesto hasta que la respuesta del
modelo calculado ajuste satisfactoriamente con los datos observados, en un intervalo de error
previamente definido en función de la escala de trabajo.
La modelización ha sido llevada a cabo utilizando el programa GM-SYS v. 2.03P
(GEOSOFT), que opera en modo directo o inverso y permite la interpretación simultánea de
los datos gravimétricos y magnéticos en 2 ½ D.
La modelización se desarrolla mediante la definición de cuerpos de sección
poligonal, con dirección perpendicular al perfil y con extensión infinita en esa dirección,
asignándoles una densidad y susceptibilidad magnética, pues representan estructuras
geológicas. La interpretación es progresivamente mejorada ajustando tanto la forma del
cuerpo como sus propiedades, hasta que se obtiene un buen ajuste con los campos
gravimétrico y magnético observados. El programa tiene también un módulo de inversión
por el que se puede ajustar aún más la respuesta del modelo, fijando bien la densidad, la
susceptibilidad magnética, o bien la posición de los vértices del polígono que representa cada
cuerpo geológico y dejando libres las demás variables para, en sucesivas iteraciones, ajustar
más el resultado de la modelización.
Los métodos utilizados para calcular las respuestas de los modelos gravimétricos y
magnéticos están basados en los métodos de Talwani et al. (1959), Talwani y Heirtzler
(1964), y utiliza los algoritmos descritos en Won y Bevis (1987). Éstos calculan en 2D la
contribución punto a punto, de cada uno de los lados de los polígonos que representan a los
distintos cuerpos geológicos, en cada punto del perfil. La anomalía calculada en cada punto
incluye la contribución de cada lado de los distintos polígonos. Los cálculos en 2 ½ D,
basados en Rasmussen y Pedersen (1979), permiten acotar la extensión lateral de algunos
cuerpos en la dirección perpendicular al perfil en función de los datos geológicos, asignando
una densidad determinada con respecto al “encajante” desde el final del cuerpo hasta +/- ∞.
Esta característica resulta muy útil a la hora de modelizar cuerpos como intrusiones
graníticas o estructuras geológicas de pequeña extensión.
129