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
No Thumbnail Available
Date
Authors
Advisors
Evaluators
Journal Title
Journal ISSN
Volume Title
Publisher
Universidad Industrial de Santander
Abstract
Description
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).
This work is framed within a geothermal research project cofinanced by ISAGEN and COLCIENCIAS called “Modeling Subsurface Structure Resistive From Surveys Magnetotelluric In Geothermal Areas Of Colombia, Case Nevado Del Ruiz” In which through a cooperative agreement between the Universidad Autónoma de México (UNAM) the Servicio Geológico Colombiano (SGC) and the Centro Internacional de Física (CIF) has been acquired magnetotelluric information composed of 105 soundings located to west of the Nevado del Ruiz. This paper illustrates the application of a methodology for processing and modeling based on 58 soundings associated with five (5) profiles. The selection of the profiles was performed using criteria such as high density of stations, proximity to a well drilled in the area and the possibility of crossing some of the surface geological structures most relevant. In the data field processing, the SSMT2000 and MTEDITOR software were used, both from the Phoenix Geophysics House, which are programmed in the basic processing algorithms such as transformation to the frequency domain programmed, the cross spectra and auto spectra, and magnetotelluric transfer function; thus is possible to obtain as a first result the MT (magnetotelluric) soundings values expressed in apparent resistivities and phases for each frequency. For the modeling state we used the winGlink software from the Geosystem house. As a first approximation 1D nversion were made. The invariant curve of the impedance tensor with Bostick algorithm (Bostick Transform), with Occam algorithm and additionally with the model resistivities discrete layers were used. After an analysis of the dimensionality is passed to the 2D inversion. The procedure consisted of rotating the impedances tensor o the main axes (as a first approximation) for profiles oriented NW (MT 1 and MT 2) and SW (MT 3) and rotate the tensor in EW direction for profiles oriented NS (MT 4 and MT 5), selection made from regional strike directionality found. For the 2 models, two electromagnetic modes (TE and TM) simultaneously in the inversion process were used. Finally, electrical resistivity models were obtained along five magnetotelluric profiles where we can observe that 2D models are more accurated due to the high complexity of the area (as best approach to information present n soundings).
This work is framed within a geothermal research project cofinanced by ISAGEN and COLCIENCIAS called “Modeling Subsurface Structure Resistive From Surveys Magnetotelluric In Geothermal Areas Of Colombia, Case Nevado Del Ruiz” In which through a cooperative agreement between the Universidad Autónoma de México (UNAM) the Servicio Geológico Colombiano (SGC) and the Centro Internacional de Física (CIF) has been acquired magnetotelluric information composed of 105 soundings located to west of the Nevado del Ruiz. This paper illustrates the application of a methodology for processing and modeling based on 58 soundings associated with five (5) profiles. The selection of the profiles was performed using criteria such as high density of stations, proximity to a well drilled in the area and the possibility of crossing some of the surface geological structures most relevant. In the data field processing, the SSMT2000 and MTEDITOR software were used, both from the Phoenix Geophysics House, which are programmed in the basic processing algorithms such as transformation to the frequency domain programmed, the cross spectra and auto spectra, and magnetotelluric transfer function; thus is possible to obtain as a first result the MT (magnetotelluric) soundings values expressed in apparent resistivities and phases for each frequency. For the modeling state we used the winGlink software from the Geosystem house. As a first approximation 1D nversion were made. The invariant curve of the impedance tensor with Bostick algorithm (Bostick Transform), with Occam algorithm and additionally with the model resistivities discrete layers were used. After an analysis of the dimensionality is passed to the 2D inversion. The procedure consisted of rotating the impedances tensor o the main axes (as a first approximation) for profiles oriented NW (MT 1 and MT 2) and SW (MT 3) and rotate the tensor in EW direction for profiles oriented NS (MT 4 and MT 5), selection made from regional strike directionality found. For the 2 models, two electromagnetic modes (TE and TM) simultaneously in the inversion process were used. Finally, electrical resistivity models were obtained along five magnetotelluric profiles where we can observe that 2D models are more accurated due to the high complexity of the area (as best approach to information present n soundings).
Keywords
Structure resistive, Processing, Modeling, Profiles MT, Rotation MT tensor, Estructura Resistiva, Procesamiento, Modelamiento, Perfiles MT, Rotación del tensor MT