APLICACIÓN DE UNA METODOLOGÍA DE PROCESAMIENTO E INTERPRETACIÓN DE INFORMACIÓN MAGNETOTELÚRICA SOBRE VARIOS PERFILES LOCALIZADOS EN LA ZONA OCCIDENTAL DEL VOLCÁN NEVADO DEL RUÍZ - COLOMBIA

Resumen

Este trabajo está enmarcado dentro de un proyecto de investigación geotérmica cofinanciado por ISAGEN y COLCIENCIAS denominado “Modelación de la Estructura Resistiva del Subsuelo a Partir de Sondeos Magnetotelúricos en Áreas Geotérmicas De Colombia, Caso Nevado Del Ruíz”, en el cual a través de un acuerdo de cooperación entre la Universidad Autónoma de México (UNAM), el Servicio Geológico Colombiano (SGC) y el Centro Internacional de Física (CIF), se ha adquirido información magnetotelúrica compuesta de 105 sondeos localizados al occidente del Nevado del Ruiz. En este trabajo se ilustra la aplicación de una metodología de procesamiento y modelamiento con base en 58 sondeos asociados a cinco (5) perfiles.

La selección de los perfiles se realizó siguiendo criterios como alta densidad de estaciones, proximidad a un pozo perforado en la zona y la posibilidad de atravesar algunas de las estructuras geológicas de superficie más relevantes. En el procesamiento de la información de campo se han utilizado los software SSMT2000 y MTEDITOR de la casa Phoenix Geophysics, en los cuales están programados los algoritmos básicos de procesamiento tales como la transformación al dominio de las frecuencias, el cálculo de los espectros cruzados y auto espectros y la función de transferencia magnetotelúrica; de esta manera es posible obtener como un primer resultado los sondeos magnetotelúricos (MT) expresados en valores de resistividades aparentes y fases para cada frecuencia.

Para la etapa de modelamiento se utilizó el software WinGLink de la casa Geosystem. Como primera aproximación se realizaron inversiones 1D. Se utilizó el invariante del tensor de impedancias con el algoritmo de Bostick (transformada de Bostick), con el algoritmo de Occam y adicionalmente con el modelo discreto de capas de resistividades. Tras un análisis de la dimensionalidad se pasó a la inversión 2D. El procedimiento consistió en rotar el tensor de impedancias a los ejes principales (como primera aproximación) para los perfiles orientados NW (MT 1 y MT 2) y SW (MT 3) y rotar el tensor en dirección E-W para los perfiles orientados N-S (MT 4 y MT 5), selección hecha a partir de la direccionalidad del strike regional encontrado. Para los modelos 2D se usaron los dos modos electromagnéticos (TE y TM) simultáneamente en el proceso de inversión. Finalmente se obtuvieron modelos de resistividad eléctrica a lo largo de los cinco perfiles magnetotelúricos en donde se puede observar que son más precisos los modelos 2D debido a la alta complejidad de la zona (puesto que se aproximan mejor a la información presente en los sondeos).