UNIVERSIDAD DE EL
SALVADOR
FACULTAD MULTIDISCIPLINARIA DE
OCCIDENTE
“Aplicación de métodos geofísicos (Gravimetría, Magnetometría,
MT y Tomografía Eléctrica) en el campo geotérmico de Berlín,
Usulután”
DOCENTE: LIC. CRISTIAN MIRÓN
PRESENTA
Tania Elizabeth Lemus Mendoza LM19039
Madeline Beatriz Morán Gil MG19050
Marcos Javier Mendoza Lorenzana ML20028
DEPARTAMENTO DE FÍSICA,
LICENCIATURA EN GEOFÍSICA
1
INDICE
INTRODUCCIÓN ....................................................................................................................... 2
GRAVIMETRÍA .......................................................................................................................... 2
Caso de estudio:Caracterización de estructuras geológicas de un área geotérmica
mediante la aplicación integrada de las técnicas de mapeo térmico, magnetometría y
gravimetría ..........................................................................................................................2
Metodología .................................................................................................................. 2
Resultados .....................................................................................................................2
Caso de estudio: Caracterización del acuífero intermedio de un sistema geotérmico de El
Salvador mediante el uso de datos geofísicos .......................................................................... 3
Metodología .................................................................................................................. 3
Resultados .....................................................................................................................4
MAGNETOMETRÍA ................................................................................................................... 5
Caso de estudio:Caracterización de estructuras geológicas de un área geotérmica mediante la
aplicación integrada de las técnicas de mapeo térmico, magnetometría y gravimetría .................5
Metodología .................................................................................................................. 5
Para poder realizar el método se utilizaron los siguientes materiales: ................................. 5
Resultados .....................................................................................................................6
MAGNETOTELURICO (MT) ...................................................................................................... 6
Caso de estudio: Modelo conceptual de un sistema geotérmico de alta entalpía; integración de
metodologías geológicas y geofísicas ..................................................................................... 6
Metodología .................................................................................................................. 6
Resultados .....................................................................................................................7
Caso de estudio: Modelo conceptual geofísico de la zona geotérmica de Chinameca, San Miguel,
El Salvador ........................................................................................................................... 8
Metodología .................................................................................................................. 8
Resultados .....................................................................................................................8
TOMOGRAFÍA ELÉCTRICA .................................................................................................... 10
Caso de estudio: Evaluación de la amenaza geológica por deslizamientos de ladera en la
plataforma del pozo TR-7 del campo geotérmico de Berlín, Usulután ..................................... 10
Metodología ................................................................................................................ 10
Resultados ...................................................................................................................10
CONCLUSIÓN ..........................................................................................................................11
BIBLIOGRAFÍA ........................................................................................................................12
2
INTRODUCCIÓN
El campo geotérmico de Berlín, en el departamento de Usulután, ha sido la base de diversos
estudios geológicos, geoquímicos y geofísicos a lo largo de los años los cuales han plasmado
la evolución del sistema geotérmico y para fines exploratorios o de explotación para el
aprovechamiento del recurso de alta entalpía en la generación de la energía eléctrica. El
presente trabajo, resume un conjunto de casos de estudio realizados en el área antes
mencionada, haciendo énfasis en la aplicación de métodos geofísicos, la metodología
implementada y los resultados obtenidos.
GRAVIMETRÍA
Caso de estudio:Caracterización de estructuras geológicas de un área geotérmica
mediante la aplicación integrada de las técnicas de mapeo térmico, magnetometría y
gravimetría
Metodología
Se utilizó el gravímetro digital Autograv Modelo CG-5 de la empresa SCINTREX.
El procedimiento de medición es el siguiente:
● Limpiar el terreno para colocar el trípode de nivelación, teniendo el cuidado de que
este apunte al norte magnético del lado de los tornillos XY, el trípode tiene un
indicador de nivelación por burbuja en una de sus patas.
● A la hora de ubicar el gravímetro debe de estar siempre vertical y con mucho cuidado
por su gran sensibilidad.
● Encender el gravímetro, en la pantalla principal SET UP se establecieron los
siguientes parámetros: nombre del levantamiento, nombre del responsable que realiza
el levantamiento, latitud geográfica del punto gravimétrico, longitud geográfica del
punto gravimétrico, acimut y elevación.
● Nivelar el equipo observando la pantalla del menú de nivelación moviendo los
tornillos de nivelación, que están en el trípode, hecho esto se da «READ GRAV» para
tomar mediciones.
● Este instrumento es programado para realizar dos ciclos de medida, una vez
terminado mostrará un resumen de las mediciones, si los datos son útiles se guardan.
A la hora de tomar las mediciones es de tener en cuenta la posición de la personas para que
no infiere en la calibración del equipo, además, que el instrumento tiene tiempo de retardo
para que la persona pueda alejarse y no afectar la medida tomada.
Resultados
Se observa en la Fig 1a la disminución de la gravedad en dirección hacia el SurOeste,esto
por el aumento de elevación por el volcán Chinchontepec.
3
La Figura 1b presenta el residual de la anomalía de Bouguer completa, donde las líneas
azules son la representación de la Quebrada Infiernillo Ciego, la línea de color verde del
perfil uno, dos y tres es la curva modelada y la línea de color rojo representa el perfil real. De
la Fig 1b los perfiles dos y tres se observa al lado izquierdo del lineamiento en superficie un
cuerpo de menor densidad que puede corresponder con la estructura geológica o una
alteración hidrotermal. En el perfil uno se observa que al lado derecho ocurre un caso
contrario que parece que el cuerpo se encuentra al lado Este de la falla inferida.
Figura 1. a) Perfiles de gravimetría. b) Perfiles 2D de anomalía de bouguer completa.
a b
Caso de estudio: Caracterización del acuífero intermedio de un sistema geotérmico de
El Salvador mediante el uso de datos geofísicos
Metodología
En este estudio se ha utilizado la gravimetría para poder determinar estructuras acuíferas
donde puede haber un flujo ascendente que sea considerable , El set de datos utilizados para
el análisis gravimétrico corresponde a la grilla de exploración del sistema geotérmico de
interés, con alrededor de 687 puntos de gravedad obtenidos en diferentes campañas.
A los datos se les realizó un filtrado por primera derivada horizontal (Figura que aparece a
continuación (derecha)), con el fin de determinar lineamientos o estructuras que estuvieran
influenciando el sistema geotérmico de interés. A partir de esto, se identificaron ciertas
4
anomalías que pueden asociarse a lineamientos tales como fallas o bordes caldericos, los
cuales fueron posteriormente cotejados con los mapas estructurales generados durante
exploraciones geológicas en la zona.
Figura 2. Perfiles de gravimetría, anomalía de bouguer, y primera derivada horizontal.
Resultados
Las anomalías gravimétricas en la región presentadas por el mapa de anomalía de Bouguer
del área de estudio aporta muy poca información para la determinación de estructuras es por
ello que la elaboración del mapa de primera derivada horizontal aporta de gran manera una
idea sobre la distribución de lineamientos o estructuras en la zona.
La interpretación del mapa de primera derivada permite apreciar al menos tres segmentos de
lineamientos con orientación NW – SE en la parte central del área del sistema geotérmico de
interés; además, en el segmento SW se puede interpretar en conjunto con información de
elevación del área una estructura similar a una caldera, la cual realiza un cierre cercano a la
región donde se tiene el primer lineamiento de la zona central del sistema geotérmico. En la
zona Norte del área de estudio se aprecia un lineamiento con orientación casi E – W que
supondría el borde de la región del graben central de El Salvador
5
Figura 3. Perfil de la primera derivada horizontal (a) y (b)perfil del la primera derivada
horizontal con la topografía del terreno.
MAGNETOMETRÍA
Caso de estudio:Caracterización de estructuras geológicas de un área geotérmica
mediante la aplicación integrada de las técnicas de mapeo térmico, magnetometría y
gravimetría
Metodología
Para poder realizar el método se utilizaron los siguientes materiales:
● Materiales y equipos
● Magnetómetro protónico.
● Brújula.
● Antena GPS.
● Bastón metálico.
Teniendo el equipo y los materiales se tienen que tener en consideración lo siguiente.
Se tiene que tener un reconocimiento total de la zona donde se realizará el estudio para poder
ver la factibilidad de acceso al lugar.
Además, poder realizar una visualización detallada de la zona para poder determinar algunos
factores que pueden dañar las mediciones como son metales, si en el sitio de mediciones, si se
encuentran metales grandes entonces se realiza una limpieza del terreno donde se realicen las
mediciones para que estas puedan ser aún mejores y no exista ruido en las mediciones.
Se debe establecer un punto base de medición para registrar las variaciones temporales del
campo magnético durante el tiempo en que se realice la medición.
Luego empezar con las medidas, tomando los puntos de GPS por punto.
6
Resultados
En la siguiente imagen se pueden observar dos áreas especialmente resaltadas por estos
valores, al sur de la zona de estudio se obtuvo una anomalía negativa, la cual responde con la
zona de mayor alteración hidrotermal observada en campo. En la zona central del área, se
muestra el valor más alto de la anomalía magnética.
Figura 4. Perfil de magnetometría, con sus variaciones en el campo magnético.
Al final se concluye que el método de magnetometría cumple con lo que se necesita el cual
era identificar la zona hidrotermal del área de estudio, y se logró los datos fueron bastante
satisfactorios, además respaldo lo visto en campo.
MAGNETOTELURICO (MT)
Caso de estudio: Modelo conceptual de un sistema geotérmico de alta entalpía;
integración de metodologías geológicas y geofísicas
Metodología
Para la toma de los datos de resistividad se realizó un perfil 1D en dirección oeste-este,
ubicado perpendicularmente a las estructuras geológicas (fallas identificadas). La distribución
de las estaciones MT fue homogénea a través del área de estudio haciendo un total de 70
sondeos MT (figura 5 a).
La metodología para el procesamiento se llevó a cabo de la siguiente manera:
7
1) procesamiento de datos magnetotelúricos, donde se obtuvo una buena calidad de los datos
con presencia de poco ruido. Para el análisis e inversión de las curvas de resistividad aparente
y fases, se implementó el programa WinGLink.
2) Procesamiento de modelamiento 1D de MT mediante WinGlink, creando una base de
datos con las coordenadas y las curvas de resistividad aparente y fase. Debido a que las
mediciones se tomaron a distintas horas del día, fue necesario realizar un filtrado de los datos
para eliminar valores obtenidos en horarios con presencia de ruido.
3) Inversión Unidimensional 1D. WinGLink utiliza un algoritmo que obtiene un modelo
suavizado y realiza un promedio de las curvas de resistividad y fase (figura 5 b), interpolada a
partir de las mediciones de las estaciones y el ajuste a un modelo discreto que permite estimar
los espesores de las capas según los valores de resistividad.
Figura 5. a) (Izquierda) Distribución de los sondeos MT y perfil. b) (derecha) Modelo 1D de
sondeos MT.
Resultados
De los resultados obtenidos en el área se lograron distinguir dos anomalías importantes
(figura 6) : Una capa conductiva superficial (de baja resistividad) con valores inferiores a 10
Ω*m que se extiende hasta los 800 mbnm aproximadamente asociada a la capa sello del
sistema y por debajo de la capa conductiva se observó la presencia de una capa resistiva con
valores entre entre 10 Ω*m y 47 Ω*m ubicado entre los 800 y 2500 mbnm y un espesor
aproximado de 2000 m, interpretado como el posible reservorio. Por último, se obtuvo una
anomalía conductiva muy profunda, a partir de los 900 mbnm, que podría asociarse a la
fuente de calor del sistema geotérmico constituido por una cámara magmática a través de la
cual fluye calor y fluidos hidrotermales.
8
Figura 6. Sección de resistividad obtenido.
Caso de estudio: Modelo conceptual geofísico de la zona geotérmica de Chinameca, San
Miguel, El Salvador
En el periodo de agosto de 2005 a marzo de 2006, La Geo realizó una campaña geofísica en
un área geotérmica de Chinameca con una extensión aproximada de 100Km2 donde se
aplicaron métodos de gravimetría, magnetometría, magnetotelúrico (MT) en combinación
con el método electromagnético en el dominio del tiempo (TDEM), con el fin de generar un
modelo conceptual del yacimiento, logrando una correcta correlación de los datos obtenidos.
A continuación, se resume la metodología empleada y resultados obtenidos en el método MT
de dicho caso de estudio:
Metodología
La metodología empleada se resume en tres etapas principales: La primera que consiste en la
determinación del área de interés, impulsada por el deseo de realizar estudios
complementarios en la zona en años anteriores (años 60’ y 90’). Una vez definida el área, se
determinó el número de estaciones MT haciendo un total de 65 y una estación remota de
referencia ubicada a 20km. La segunda etapa que contempla obtención de los datos, en los
meses de febrero y marzo del 2006, los cuales presentaron una buena calidad, salvo en las
cercanías a zonas pobladas donde se planteó que la distribución de la red eléctrica pudo
influir en las mediciones; por último, se realizó la etapa del procesamiento de los datos
considerando una rejilla de puntos de 250 x 250 m, donde con un total de 50 iteraciones se
obtuvo un RMS final de 1.58. En base al procedimiento anterior el archivo de salida se
exportó al software WingLink para generar los mapas y secciones para una mejor
visualización de los resultados e interpretaciones.
Resultados
La representación visual del modelo obtenido con las mediciones de MT se muestra en la
figura x1. De este modelo, se observaron con claridad la disminución de la resistividad
9
partiendo de la dirección Oeste-Este y capas conductivas en los primeros 500 m
aproximadamente, dicho resultado describe el comportamiento típico de un sistema de alta
entalpía. Por otro lado, la discontinuidad en los valores de resistividad en la parte central del
perfil se asocia a la presencia del sistema de fallas NS en el cerro El Limbo, dicha afirmación
se correlaciona con la información obtenida de estudios magnéticos del área de estudio.
Figura 7. a) Distribución de los sitios MT en el área de estudio (izquierda) y b) modelo de
resistividad (derecha).
La distribución de los valores de resistividad muestra el comportamiento característico de un
volcán geotérmico andesítico, figura 8 (izquierda), donde aprecia la similitud entre el área
geotérmica de Chinameca (centro) y Berlín (derecha).
Figura 8. Comparación de un modelo teórico de resistividades de MT en un sistema
volcánico geotérmico andesítico (izquierda) y la secuencia resistiva en el área de Chinameca
(centro) y el ejemplo de Berlín (derecha).
Las estructuras geológicas más representativas del área geotérmica de Chinameca son: A) una
capa resistiva en los primeros metros de profundidad asociado a la presencia de rocas con
alteración hidrotermal debido a la presencia de fluidos geotérmicos distribuidas hacia el
Norte del perfil. B) Una capa conductora con valores entre 2 y 12 ohm-m, situada por debajo
de los 500 m de profundidad asociada a la presencia de arcillas alteradas correspondiente a la
capa sello del sistema, con un espesor estimado entre 400 y 700 m y por último C) una zona
10
de muy alta resistividad en forma de cúpula, correspondiente a la zona del yacimiento de alta
temperatura con presencia de minerales alterados.
TOMOGRAFÍA ELÉCTRICA
Caso de estudio: Evaluación de la amenaza geológica por deslizamientos de ladera en la
plataforma del pozo TR-7 del campo geotérmico de Berlín, Usulután
Metodología
Para la adquisición de los datos se utilizó el arreglo “Wenner” para distancia de 3 y 5m,
donde se registran todos los valores de resistividad aparente. Dichas mediciones se realizaron
empleando el equipo “Terrameter SAS 100” el cual obtiene perfiles 2D y modelos de
resistividad eléctrica (2D), utilizando algoritmos matemáticos para la inversión de datos, Este
equipo contiene: 4 canales de entrada, incluyendo uno para batería, cable de entrada para
corriente directa (CD), cable RS 232 y electrodos de corriente y potencial.
Para la interpretación de las resistividades aparentes se utilizó RES2DINV y posteriormente
sometidos en el programa MATLAB donde se realizaron dos códigos, uno con el objetivo de
ajustar los datos de elevación a coordenadas no reales que mejor se ajustaran a la topografía
del modelo de resistividad eléctrica sin alterar la topografía y el segundo para leer los los
nuevos datos generados de topografía y valores de resistividad eléctrica producto de la
inversión realizada en el programa RES2DINV.
Resultados
En el perfil SO se observa una anomalía de baja resistividad Fig #a, donde se asocia a la falla
san jos|e, donde existe ascenso de fluidos y la depositación de material con alto grado de
alteración hidrotermal, además, de identificar dos capas: una profunda de 50 ohm.m y una
superficial con valores de 3 om.m, en el perfil B1 se identifica las mismas capas que el
modelo A2, donde este alcanzó valores de 300 ohm.m y un valor mínimo de 0.3 ohm.m.
En la Figura 9, hacia NO se muestra una anomalía triangular, disminuyendo a medida que se
acerca a la superficie confirmando que existe una zona de fracturamiento o debilidad que
facilita el ascenso de fluidos a la superficie. Las capas predominantes son arcilla ( 0-2 ohm.m)
y rocas con alto grado de alteración debido a las manifestaciones hidrotermales. El perfil C1
muestra en NO tipos de rocas son lavas con menor grado de alteración hidrotermal, es el
único perfil que presenta una primera capa de mayor resistividad por lo que se cree que
contiene roca dura fragmentada y hacia SE cambia los valores considerable lo que significa
que hay mayor alteración hidrotermal. El perfil B2 presenta valores de 0.2 ohm.m en la parte
central (posiblemente una canalización).
11
Figura 9. Modelos de resistividad eléctrica: a) Perfiles SO-NE (B1 y A2), b) Perfiles NO-SE
(A1 y C1) y c) Perfil B2.
CONCLUSIÓN
Los diferentes métodos Geofísicos son esenciales para la exploración de la Geotermia, donde
dependerá la elección del método a partir de los objetivos de la investigación, los alcances y
limitaciones, cabe destacar que la integración de dos o más métodos para un mismo estudio
es importante para obtener una mayor precisión de lo que yace en el subsuelo y asociarlo a la
información previa, con ésto se debe tener en cuenta que los resultados obtenidos deben ser
geológicamente y geofisicamente admisibles, para la posterior generación de modelos
conceptuales claves en la toma de decisiones referentes a la exploración y explotación del
recurso.
12
BIBLIOGRAFÍA
José, A. G. R., Daniel, E. G. J., Beatriz, M. S. I., Mario, R. M., Arturo, T. L. J., Santos, P., &
Rivas, M. J. Caracterización de estructuras geológicas de un área geotérmica mediante
la aplicación integrada de las técnicas de mapeo térmico, magnetometría y gravimetría.
Mayra, R. H. C., & Alvaro, J. C, R., (2013), Evaluación de la amenaza geológica por
deslizamientos de ladera en la plataforma del pozo TR-7 del campo geotérmico de
Berlín, Usulután
Mayen, L. A. (2019). Caracterización del acuífero intermedio de un sistema geotérmico de El
Salvador mediante el uso de datos geofísicos. Diplomado en geotermia para América
Latina.
Santos, P., & Rivas, J. (2009). Modelo conceptual geofísico de la zona geotérmica de
Chinameca, San Miguel, El Salvador . UNU-GTP y LaGeo.
Vargas, N. B., & Miraval, C. V. (20019). Modelo conceptual de un sistema geotérmico de
alta entalpía; integración de metodologías geológicas y geofísicas. Diplomado en
geotermia para América Latina.
