GEO6069v1

GOBIERNO DE CHILE
MINISTERIO DE OBRAS PÚBLICAS

DIRECCIÓN GENERAL DE AGUAS
REGIÓN DE ÑUBLE
ANÁLISIS Y CARACTERIZACIÓN DEL

ACUÍFERO DE LA CUENCA DEL RÍO
ITATA, ESTUDIO GEOFÍSICO VALLE
CENTRAL PROFUNDA
INFORME FINAL
REALIZADO POR:
GEO RESOURCE LTDA.
S.I.T. Nº 545
Chillán, marzo de 2024

MINISTERIO DE OBRAS PÚBLICAS
Director General de Aguas
Sr. Rodrigo Sanhueza Bravo
Director Regional de Aguas de Ñuble

Sr. Marcelo Godoy Toro
Inspector Fiscal:
Dr. Héctor Flores Moraga
Geo Resource LTDA.

Jefe de Proyecto
Geofísico Felipe Daza Hernández
Profesionales Especialistas
Geólogo Pablo Alarcón Pinto
Geofísico Felipe Reyes Córdova (MSc)
Geofísico Andrés Stuardo Espinoza
Geofísica Daniela Montecinos Cuadros (MSc)
Geofísico Daniel González Parada
Geólogo Ian Espinoza Gutiérrez
Lic. en Geología Javiera Venegas González
Lic. en Geología Cristian Roncagliolo Belmar
Estudio financiado por el Gobierno Regional de Ñuble

2022-2023
Índice
1. INTRODUCCIÓN 2
1.1. UBICACIÓN 4
1.2. OBJETIVOS 6
1.3. ANTECEDENTES 6
1.4. CONCLUSIONES DEL ANÁLISIS DE ANTECEDENTES 18

2. CONTEXTO GEOLÓGICO E HIDROGEOLÓGICO 20
2.1. UNIDADES LITOLÓGICAS 20
2.2. UNIDADES HIDROGEOLÓGICAS 22
2.3. CAPTACIONES PROFUNDAS EN LA CUENCA DEL ITATA 24

3. METODOLOGÍA 34
3.1. GRAVIMETRÍA 34
3.1.1. Terreno 35
3.1.2. Procesamiento 37
3.1.3. Modelamiento 46
3.2. MAGNETOTELÚRICA 52
3.2.1. Terreno 53
3.3. TRANSIENTE ELECTROMAGNÉTICO 62

3.3.1. Terreno 64
4. RESULTADOS E INTERPRETACIÓN 73
4.1. MAGNETOTELÚRICA 73
4.1.1. Funciones de transferencia 73
4.1.2. Análisis de dimensionalidad 75
4.1.3. Modelos de resistividad eléctrica 77
4.2. TRANSIENTE ELECTROMAGNÉTICO 131

4.2.1. Modelos de resistividad 1D 131
4.3. GRAVIMETRÍA 137
4.3.1. Perfiles E-W (L1-L7) 140
4.3.2. Perfiles auxiliares (D1-D14) 150
4.3.3. Mapas de espesor sedimentario y profundidad de basamento 157
4.3.4. Comparación con estudios anteriores 160
5. CONCLUSIONES Y RECOMENDACIONES 170
6. REFERENCIAS 173
ÍNDICE DE FIGURAS
Figura 1.1: Ubicación de la cuenca del Itata en la Región del Ñuble. Fuente: Elaboración
propia. .................................................................................................................. 5
Figura 1.2: Localización de estaciones geofísicas en estudios previos en las regiones del
Ñuble y Biobío, relevantes para este estudio. Fuente: Elaboración propia. ................... 17
Figura 1.3: Estructura de las fichas bibliográficas presentadas en el Anexo H ............... 18
Figura 2.1: Mapa geológico de la Cuenca del Río Itata. Modificado de DGA, 2011. ........ 21
Figura 2.2: Mapa hidrogeológico de la Cuenca del Río Itata. Modificado de DGA, 2015. . 23
Figura 2.3: Mapa de captaciones subterráneas registradas en la Dirección General de
Aguas Subterráneas en la Cuenca del Río Itata. Elaboración Propia ............................. 25
Figura 2.4: Columnas estratigráficas de las captaciones que detectan roca (Perfiles
gravimétricos L1 y L2) ubicadas en la cuenca del Río Itata. ........................................ 27
Figura 2.5: Columnas estratigráficas de las captaciones que detectan roca (Perfil
gravimétrico L3) ubicadas en la cuenca del Río Itata. ................................................ 28
gravimétricos L6 y L7) ubicadas en la cuenca del Río Itata. ........................................ 31
gravimétricos D7A, D10 y D9B) ubicadas en la cuenca del Río Itata. ........................... 32
Figura 2.10: Plano de ubicación de los perfiles elaborados a partir de las captaciones
seleccionadas. ...................................................................................................... 33
Figura 3.1: Flujo de adquisición, procesamiento y modelamiento de datos gravimétricos.
.......................................................................................................................... 35
Figura 3.2: Mapa de ubicación de estaciones gravimétricas medidas y estación de
gravedad absoluta utilizada para procesamiento. ...................................................... 36
Figura 3.3: Gravedad Absoluta calculada para la zona de estudio. .............................. 38
Figura 3.4: Figura esquemática de la corrección por altura o aire libre. ....................... 39
Figura 3.5: Figura esquemática de la corrección de Bouguer o de placa. ...................... 40
Figura 3.6: Figura esquemática de la corrección topográfica....................................... 40
Figura 3.7: Modelo de elevación topográfica (DEM) utilizado para corrección topográfica.
.......................................................................................................................... 42
Figura 3.8: Anomalía de Bouguer Completa (ABC) calculada para la zona de estudio. ... 43
Figura 3.9: Campo regional calculado para la zona de estudio. En diamantes negros los
puntos de medición de gravedad sobre roca. ............................................................ 44
Figura 3.10: Anomalía de Bouguer Completa reducida calculada para la zona de estudio.
.......................................................................................................................... 45
Figura 3.11: Ejemplo de cuerpo en 2D (Shuey y Pasquale 1973). ............................... 46
Figura 3.12: Mapa de ubicación de perfiles gravimétricos modelados. ......................... 47
Figura 3.13: Ejemplo de dato outlayer, perfil L7, punto de medición 304. .................... 48
Figura 3.14: Ejemplo de modelamiento e integración de información de otras
metodologías para el perfil L1. Elaboración propia. ................................................... 50
Figura 3.15: Modelación del Perfil L1 con distintas densidades de sedimentos (1,9; 2,0 y
2,1 kg/m3). Elaboración propia. .............................................................................. 51
Figura 3.16: Instalación de una estación de magnetotelúrica. .................................... 52
Figura 3.17: Flujo de adquisición, procesamiento y modelamiento de datos MT. ........... 53
Figura 3.18: Mapa de distribución de estaciones MT. ................................................. 54
Figura 3.19: Series de tiempo de una estación de magnetotelúrica (TRAP5). ................ 56
Figura 3.20: Curvas de resistividad aparente y fase de la estación TRAP5. Eje x de la
figura corresponde al periodo en unidades de segundos. ........................................... 58
Figura 3.21: Visualización de las curvas de resistividad aparente (arriba) y de la fase
(abajo). Eje x de las figuras corresponde al periodo en unidades de segundos Derecha:
Visualización modelo 2D en Winglink. Paleta de colores a la derecha indica valores de
resistividad. ......................................................................................................... 60
Figura 3.22: Figura esquemática de cómo se realiza la adquisición de datos de una
estación TEM y sus fundamentos físicos de interacción con el subsuelo. Modificado de
Dentith & Mudge, 2014. ........................................................................................ 63
Figura 3.23: Flujo de adquisición, procesamiento y modelamiento de datos TEM. ......... 63
Figura 3.24: Mapa de distribución de estaciones TEM. ............................................... 65
Figura 3.25: Esquema de medición. (a) Pulso transmisor. (b) Fuerza electromagnética
inducida. (c) Campo magnético secundario medido en distintas gates o ventanas de
tiempo. ............................................................................................................... 66
Figura 3.26: Ejemplo de stacking de estación ITA44. En color rojo, azul y gris los datos
correspondientes a LM, HM y ruido respectivamente. En líneas con trasparencia se
muestran todas las mediciones realizadas, y en línea sólida los valores apilados (stacking)
de cada momento. ................................................................................................ 67
Figura 3.27: Ejemplo de cálculo de resistividad aparente para estación ITA44. ............. 69
Figura 3.28: Ejemplo de curva final de decaimiento de voltaje de la estación ITA44,
obtenida de las curvas de LM y HM traslapadas. ....................................................... 70
Figura 3.29: Ejemplo de limpieza de datos de curva de decaimiento de voltaje de la
estación ITA44. En gris los puntos que se han eliminado para la inversión. .................. 71
Figura 3.30: Modelo de capas y ajuste respecto a la curva de decaimiento de voltaje de la
estación de TEM ITA44. ......................................................................................... 72
Figura 4.1: Curvas de resistividad y fase de las estaciones TRAP1 y YUN4. .................. 74
Figura 4.2: Sección de tensor de fase y parámetro skew de uno de los perfiles del sector
de Maravilla (MAR2 y MAR3). ................................................................................. 75
Figura 4.3: Cálculo de strike para el perfil MAR2-3 ubicado en el sector de Maravilla. .... 76
Figura 4.4: Ubicación de las estaciones de magnetotelúrica en el sector de Trapiche. .... 78
Figura 4.5:Modelo de resistividad eléctrica y sección interpretada de los puntos MT TRAP1
a TRAP6. ............................................................................................................. 80
Figura 4.6: Ubicación de las estaciones de magnetotelúrica en el sector de Maravilla. ... 81
Figura 4.7: Modelo de resistividad eléctrica y sección interpretada de los puntos MT MAR2
– MAR3. .............................................................................................................. 83
Figura 4.8: Ubicación de las estaciones de magnetotelúrica en el sector de Cato. ......... 84
Figura 4.9: Modelo de resistividad eléctrica y sección interpretada de los puntos MT
CATO5 – CATO6. .................................................................................................. 86
Figura 4.10: Ubicación de las estaciones de magnetotelúrica en el sector de Cachapoal. 87
CACHA1 – CACHA2. .............................................................................................. 89
Figura 4.12: Ubicación de las estaciones de magnetotelúrica en el sector Camino
Cachapoal............................................................................................................ 91
CCACHA12 – CCACHA11. ....................................................................................... 92
Figura 4.14: Ubicación de las estaciones de magnetotelúrica en el sector de Ninquihue. 93
Figura 4.15:Modelo de resistividad eléctrica y sección interpretada de los puntos MT QUI3
– QUI4. ............................................................................................................... 95
Figura 4.16: Ubicación de las estaciones de magnetotelúrica en el sector El Ala. .......... 96
Figura 4.17: Modelo de resistividad eléctrica y sección interpretada del punto MT ALA3. 97
Figura 4.18: Ubicación de las estaciones de magnetotelúrica en el sector de Coihueco. . 98
Figura 4.19: Modelo de resistividad eléctrica y sección interpretada de los puntos MT COI4
– COI3. ............................................................................................................. 100
Figura 4.20: Ubicación de las estaciones de magnetotelúrica en el sector de Viña Ribera.
........................................................................................................................ 101
VRC1 – VRC2. .................................................................................................... 103
Figura 4.22: Ubicación de las mediciones de magnetotelúrica en el sector de Bellavista.
........................................................................................................................ 104
Figura 4.23: Modelo de resistividad eléctrica y sección interpretada de los puntos MT RC1
– RC2 – RC3. ..................................................................................................... 106
Figura 4.24: Ubicación de las estaciones de magnetotelúrica en el sector de Liucura. .. 107
LIU11 – LIU8 – LIU6. .......................................................................................... 109
Figura 4.26: Ubicación de las estaciones de magnetotelúrica en el sector de Pinto. ..... 110
Figura 4.27: Modelo de resistividad eléctrica y sección interpretada de los puntos MT PIN2
– PIN1. ............................................................................................................. 112
Figura 4.28: Ubicación de la medición de magnetotelúrica en el sector de Cartago. ..... 113
Figura 4.29: Modelo de resistividad eléctrica y sección interpretada de la estación de MT
CAR1. ............................................................................................................... 115
Figura 4.30: Ubicación de las estaciones de magnetotelúrica en el sector de Atacalco. 116
ATA8 – ATA7. ..................................................................................................... 118
Figura 4.32: Ubicación de las estaciones de magnetotelúrica en el sector de San Antonio.
........................................................................................................................ 119
SAN7 – SAN6. .................................................................................................... 121
Figura 4.34: Ubicación de las estaciones de magnetotelúrica en el sector de Trehualemu.
........................................................................................................................ 122
Figura 4.35:Modelo de resistividad eléctrica y sección interpretada de los puntos MT TRE3
– TRE4. ............................................................................................................. 124
Figura 4.36: Ubicación de las estaciones de magnetotelúrica en el sector de Yungay. .. 125
YUN1 – YUN2. .................................................................................................... 127
Figura 4.38: Ubicación de las estaciones de magnetotelúrica en el sector de Huepil. ... 128
Figura 4.39: Modelo de resistividad y secciones interpretadas de los puntos MT HUE1 –
HUE2. ............................................................................................................... 130
Figura 4.40: Modelos de inversión 1D de las mediciones que se encuentran en el extremo
este cercanos a la Cordillera Principal .................................................................... 133
Figura 4.41: Izquierda: Modelos de inversión 1D de las mediciones que se encuentran Sur
de Maravilla y en Trapiche. Derecha: Modelo de inversión 1D al Oeste de Trapiche. .... 134
Figura 4.42: Modelos de inversión 1D de las mediciones que se encuentran en El Ala y al
Sur de El Ala. ..................................................................................................... 135
Figura 4.43: Modelos de inversión 1D de las mediciones que se encuentran en la franja
entre Bulnes y Cabrero. ....................................................................................... 136
Figura 4.44: Localización de puntos de control utilizados para modelamiento gravimétrico
2D, obtenidos de este estudio, capaciones subterráneas y estudios previos. Fuente:
Elaboración propia. ............................................................................................. 139
Figura 4.45: Modelamiento e integración de información de otras metodologías del perfil
L1..................................................................................................................... 141
Figura 4.46: Modelamiento e integración de información de otras metodologías de perfil
L2..................................................................................................................... 142
L3..................................................................................................................... 143
L4..................................................................................................................... 145
L5..................................................................................................................... 146
L6..................................................................................................................... 147
L7..................................................................................................................... 149
D1. ................................................................................................................... 151
D2. ................................................................................................................... 152
D3. ................................................................................................................... 153
D4. ................................................................................................................... 154
D5. ................................................................................................................... 156
Figura 4.57: Visión general en planta del espesor sedimentario. ............................... 158
Figura 4.58: Visión general en planta de la profundidad del basamento. .................... 159
Figura 4.59: Visión en planta de espesor sedimentario estimado por DGA (2011). Imagen
de elaboración propia, a partir de datos de DGA (2011). ......................................... 161
Figura 4.60: Visión en planta de espesor sedimentario estimado por DGA (2012). Imagen
modificada de DGA (2012). .................................................................................. 162
Figura 4.61: Modelo de inversión de relleno sedimentario obtenido por DGA (2015).
Imagen obtenida directamente de DGA (2015). ...................................................... 164
Figura 4.62: Anomalía de Bouguer completa residual obtenida por DGA (2015). Obtenida
directamente de DGA (2015). .............................................................................. 165
Figura 4.63: Modelo de profundidad del basamento estimado por Aguirre et al., 2022.
Imagen obtenida directamente de Aguirre et al., 2022. ........................................... 167
Figura 4.64: Anomalía de Bouguer Completa residual de la base de datos gravimétrica
conjunta del presente trabajo, DGA (2011) y DGA (2012). Círculos negros: presente
estudio; Triángulos blancos: DGA (2011); Triángulos azules: DGA (2012). ................ 169
ÍNDICE DE TABLAS
Tabla 2.1: Tabla de información de los pozos profundos recopilados que detectan roca en
la cuenca del Rio Itata. Base cartográfica en WGS84 H18S. ....................................... 26
Tabla 3.1: Esquema de adquisición de las estaciones de magnetotelúrica .................... 55
Tabla 4.1: Ubicación de estaciones geofísicas MT realizadas en el sector de Trapiche
(WGS84 - 18S). ................................................................................................... 77
Tabla 4.2: Ubicación de estaciones geofísicas MT realizadas en el sector de Maravilla
(WGS84 – 18S). ................................................................................................... 81
Tabla 4.3: Ubicación de estaciones geofísicas MT realizadas en el sector de Cato (WGS84
– 18S). ............................................................................................................... 84
Tabla 4.4: Ubicación de estaciones geofísicas MT realizadas en el sector de Cachapoal
(WGS84- 18S). .................................................................................................... 87
Tabla 4.5 Ubicación de estaciones geofísicas MT realizadas en el sector de Camino
Cachapoal (WGS84 – 18S)..................................................................................... 90
Tabla 4.6: Ubicación de estaciones geofísicas MT realizadas en el sector de Ninquihue
(WGS84 – 18S). ................................................................................................... 93
Tabla 4.7: Ubicación de estaciones geofísicas MT realizadas en el sector de El Ala (WGS84
– 18S). ............................................................................................................... 96
Tabla 4.8: Ubicación de estaciones geofísicas MT realizadas en el sector de Coihueco
(WGS84 – 18S). ................................................................................................... 98
Tabla 4.9: Ubicación de estaciones geofísicas MT realizadas en el sector de Viña Ribera
(WGS84 – 18S). ................................................................................................. 101
Tabla 4.10: Ubicación de estaciones geofísicas MT realizadas en el sector de Bellavista
(WGS84 – 18S). ................................................................................................. 104
Tabla 4.11: Ubicación de estaciones geofísicas MT realizadas en el sector de Liucura
(WGS84 – 18S). ................................................................................................. 107
Tabla 4.12: Ubicación de estaciones geofísicas MT realizadas en el sector de Pinto
(WGS84 – 18S). ................................................................................................. 110
Tabla 4.13: Ubicación de estaciones geofísicas MT realizadas en el sector de Cartago
(WGS84 – 18S). ................................................................................................. 113
Tabla 4.14: Ubicación de estaciones geofísicas MT realizadas en el sector de Atacalco
(WGS84 – 18S). ................................................................................................. 116
Tabla 4.15: Ubicación de estaciones geofísicas MT realizadas en el sector de San Antonio
(WGS84 – 18S). ................................................................................................. 119
Tabla 4.16: Ubicación de estaciones geofísicas MT realizadas en el sector de Trehualemu
(WGS84 – 18S). ................................................................................................. 122
Tabla 4.17: Ubicación de estaciones geofísicas MT realizadas en el sector de Yungay
(WGS84 – 18S). ................................................................................................. 125
Tabla 4.18: Ubicación de estaciones geofísicas MT realizadas en el sector de Huepil
(WGS84 – 18S). ................................................................................................. 128
Tabla 4.19: Ubicación de estaciones geofísicas TEM de mejor calidad para la
interpretación geológica y geofísica en la cuenca Itata (WGS84 – 18S). .................... 131
RESUMEN
El presente informe hidrogeológico desarrollado por GEO RESOURCE LTDA. a solicitud de
la Dirección General de Aguas (DGA), presenta en detalle la búsqueda de investigaciones
previas, la recopilación de información bibliográfica, contexto geológico e hidrogeológico de
la cuenca y el resultado de las diferentes metodologías geofísicas realizadas para
determinar la profundidad del basamento rocoso presente en la zona acuífera de la Cuenca
del Río Itata en la región del Ñuble.
Para esto, se realizaron un total de 535 mediciones de gravimetría, 100 mediciones de

magnetotelúrica (MT) y 100 mediciones de transiente electromagnético (TEM), con las
cuales se pudo caracterizar la profundidad del basamento y la disposición de las unidades
que rellenan el Valle central.
La información obtenida con TEM en la Cordillera de la Costa, permitió reconocer la

profundidad del basamento rocoso, así como también conocer los valores de resistividad
en superficie.
Los modelos de MT, ubicados principalmente en la Depresión Intermedia muestran en

general la existencia de una zona conductiva en profundidad. Esta zona podría interpretarse
como el miembro inferior de Fm. Mininco, así como también podría corresponder a alguna
unidad geológica no descrita, o parte del basamento muy alterado y/o fracturado. Por otro
lado, en la vertiente occidental de la Cordillera Principal, fue interpretada la ocurrencia de
unidades como Fm. La Montaña (Malleco?), Fm. Mininco y Fm. Cura-Mallín.
Finalmente, el modelo de espesor sedimentario obtenido a partir de la gravimetría sugiere

la existencia de un depocentro ubicado entre las ciudades de San Carlos y Coihueco que
alcanzaría un espesor máximo de 2.080 m al oeste de San Carlos. La cuenca hacia el sur
se vuelve menos profunda, alcanzando un espesor de 800 m aproximadamente, en las
proximidades de Pemuco.
Estos resultados son consistentes con algunos modelos de espesor sedimentario

presentados previamente en la zona de interés, pero que tenían una menor extensión y
resolución espacial que la desarrollada en este estudio.
I
1. INTRODUCCIÓN
La Dirección General de Aguas (DGA) es el organismo estatal encargado de gestionar,

verificar y difundir información sobre la cantidad y calidad del agua en el país. Entre sus
funciones destacan el estudio constante de las fuentes de recursos hídricos, con el objetivo
principal de asegurar su aprovechamiento sustentable tanto en el presente como en el
futuro.
El sostenido crecimiento económico y desarrollo social de las últimas décadas ha generado

demandas cada vez mayores sobre los recursos hídricos por parte de los diferentes tipos
de usuarios. Este crecimiento sostenido, se produce en un contexto en el cual los recursos
hídricos, en una buena medida, ya están comprometidos en el abastecimiento de los usos
actuales. La relación extracción/disponibilidad de agua manifiesta una separación
importante debido a la fuerte presión sobre los recursos hídricos.
La escasez hídrica se ha agudizado en los últimos años. En efecto, se observa una

significativa disminución en las precipitaciones de la región de Ñuble, aumento de solicitud
de aguas subterráneas y una fuerte presión de los diferentes actores sobre los recursos
hídricos existente. Por lo cual, se requiere información fidedigna de la capacidad del acuífero
y sistemas naturales para evitar impactos antrópicos en el acuífero por sobre explotación,
ante eventos extremos de déficit hídrico, aumento de la superficie de riego, crecimiento de
la población, entre otros.
La Dirección General de Aguas actualmente posee una Red de Medición de Aguas

Subterráneas para la región de Ñuble, y considerando que la necesidades de información
son cada vez más exigentes respecto a tiempo y oportunidad de su acceso, esta debe ser
fortalecida respecto a la densidad de las redes de monitoreo, para ser capaz de medir las
variaciones espaciales de los niveles piezométricos y calidad, y de esta forma, generar
información fundamental para la creación y retroalimentar a los modelos numéricos, que
aporte a la toma de decisión.
El desarrollo de modelos de simulación que permitan planificar y gestionar los recursos

hídricos, el que, mediante intervención programada e introducción directa, indirecta o
inducida de agua en un acuífero, incremente el grado de ganancia y disponibilidad de los
recursos hídricos, permitirá desarrollar protocolos, reglas, planes, proyectos, entre otros.
2
Esto hace necesario conocer la línea base de los reservorios subterráneos, sus interacciones
con el medio y sus limitaciones.
En la actualidad, para iniciar el proceso de búsqueda y generación de respuestas y

herramientas que permitan mejorar la administración del recurso hídrico en la región de
Ñuble, se encuentra en ejecución el estudio “Análisis y caracterización del acuífero de la
cuenca del río Itata” ID: 1019-125-LR21, que busca caracterizar el comportamiento y la
disponibilidad de los recursos subterráneos de la región de Ñuble, información necesaria
para realizar una gestión adecuada de los recursos hídricos, tanto superficiales como
subterráneos y sus interacciones. El actual escenario de crisis climática ha afectado de
manera importante la disponibilidad de agua, teniendo como efectos: el aumento de
conflictos sociales, impactos en la economía agraria, entre otras. Lo anterior hace urgente
contar con herramientas de gestión basadas en información a una escala adecuada que
permita planificar, gestionar y evaluar el comportamiento de los sistemas hídricos y así
evitar una sobreexplotación de los recursos y asegurar la disponibilidad para la vida y salud
de las personas.
Sin embargo, para el correcto desarrollo e implementación de dicha herramienta y/o

modelo es necesario contar con una geometría detallada en profundidad que permita
optimizar los análisis y disminuir la incertidumbre al momento de desarrollar y ajustar los
modelos tanto conceptuales como numéricos de la cuenca del Río Itata. Dado lo anterior
es que surge la necesidad de llevar a cabo el presente estudio que permitirá generar
información geofísica clave y de esta forma, obtener un modelo en profundidad del relleno
sedimentario y el basamento rocoso de la cuenca.
El presente informe hidrogeológico desarrollado por GEO RESOURCE LTDA. a solicitud de

la Dirección General de Aguas (DGA), presenta en detalle la búsqueda de investigaciones
previas, la recopilación de información bibliográfica, contexto geológico e hidrogeológico de
la cuenca y el desarrollo de las diferentes metodologías geofísicas realizadas.
Para el presente estudio, en la cuenca del Itata Valle Central profundo, se realizaron
campañas geofísicas mediante los métodos: gravimétrico, magnetotelúrico (MT) y
transiente electromagnético (TEM). Los modelos generados a partir de este levantamiento
geofísico presentan la variación de los valores de resistividad eléctrica y de densidades a
través de la cuenca, lo que permite la interpretación tanto de las características del relleno
sedimentario como de la profundidad del basamento rocoso.
3
Para una mejor compresión y representación de los datos, se presentan tanto los resultados
obtenidos utilizando cada método geofísico como la interpretación conjunta entre las
distintas metodologías utilizadas. Los trabajos de terreno se ejecutaron durante el periodo
comprendido entre octubre del año 2023 y febrero de 2024.
Los Anexos comprenden una variedad de documentos, que incluyen planos y mapas
elaborados, captaciones revisadas, modelos de inversión 1D para las estaciones de TEM,
2D para las estaciones de magnetotelúrica y 2D para los perfiles gravimétricos. Además,
se incluye el registro fotográfico desarrollado durante las campañas geofísicas que
documenta el progreso de las labores realizados en terreno.
1.1. Ubicación
La Región de Ñuble comprende una superficie de 13.104,41 km 2, de los cuales 11.326,37

km2 (86%) están representados por la cuenca del río Itata, la cual cruza la región desde el
límite internacional por el este hasta el Océano Pacífico hacia el oeste.
La cuenca del río Itata presenta una intrincada red hidrográfica superficial que incluye 3
cursos de agua principales, correspondientes a los ríos Itata, Ñuble y Diguillín. La cabecera
de la cuenca se ubica en la zona de la Cordillera de los Andes, conocida a estas latitudes
como Cordillera Principal, la cual representa las mayores alturas del sector, como es el caso
los volcanes Chillán y Nevados de Chillán, con aproximadamente 3.200 m s.n.m. Esta
cordillera se compone principalmente de rocas volcánicas y volcanosedimentarias, cuyo
origen se asocia a procesos eruptivos pasados, y rocas de origen intrusivo. A medida que
disminuye la altura geográfica hacia el oeste inicia la transición hacia la precordillera, con
alturas geográficas que fluctúan entre los 300 a 800 m s.n.m. Esta zona se caracteriza por
la presencia de una diversidad de litologías, entre ellas rocas plutónicas, volcánicas y
volcanosedimentarias.
Al continuar hacia el oeste, inicia la una amplia explanada que representa el Valle Central,
el cual conforma una zona deprimida entre la Cordillera de la Costa y la Cordillera Principal.
En este sector se han acumulado grandes potencias de relleno sedimentario y
volcanosedimentario, cuyo origen se asocia a procesos erosivos y de transporte productor
de la densa red de drenaje que surca el valle, y también por la acción de procesos eruptivos
pasados.
4
La vertiente oriental de la Cordillera de la Costa marca el aumento de las elevaciones que
llegan hasta los 900 m s.n.m. aproximadamente. Esta unidad se compone de un antiguo
basamento rocoso altamente erosionado, representado por rocas intrusivas y en menor
medida metamórficas. Finalmente, hacia la desembocadura del río Itata disminuyen
abruptamente las elevaciones con el desarrollo de una restringida planicie costera.
La Figura 1.1 representa el área de estudio, definida por el mandante a partir de un archivo
.kmz en las bases de la licitación (Itata_profunda.kmz). En esta área del Valle Central
profundo, se desarrollaron las campañas geofísicas utilizando las metodologías
magnetotelúricas (MT), transiente electromagnético (TEM) y gravimetría.
Figura 1.1: Ubicación de la cuenca del Itata en la Región del Ñuble. Fuente: Elaboración propia.
5
1.2. Objetivos
Objetivo General
El objetivo general del presente estudio es generar nueva información geofísica que permita
caracterizar en profundidad el relleno sedimentario y el basamento rocoso presente en la
zona acuífera de la cuenca del río Itata.
Objetivos Específicos
• Recopilar y analizar antecedentes geológicos y geofísicos de la zona de estudio.

• Compilar los resultados de los estudios geofísicos ejecutados con anterioridad.
• Adquirir nuevos datos de gravimetría, transiente electromagnético y
magnetotelúrica.
• Identificar las características de los estratos sedimentarios que componen el relleno
de la cuenca del río Itata.
• Definir la geometría del basamento rocoso de la cuenca del río Itata.
• Integrar los resultados obtenidos con los datos de prospecciones geofísicas
anteriores.
1.3. Antecedentes
En este subcapítulo se exponen los trabajos e informes recopilados que aportan a los
objetivos de este estudio con respecto a la Geología, Hidrogeología y Geofísica. Se resume
la información que se consideró relevante de cada estudio para la realización de este
informe.
• Mapas geológicos preliminares de Chile, Hoja Los Ángeles – Angol. Escala 1:250.000
(Instituto de Investigaciones Geológicas, 1981)
Carta de las unidades geológicas a una escala de 1:250.000 en el área comprendida entre
Los Ángeles y Angol. Se identifican los afloramientos geológicos más relevantes y
características morfoestructurales que se presentan en el extremo sur de la cuenca de
Itata.
6
Se presentan formaciones geológicas presentes en el área de estudio, tales como Fm.
Mininco y Fm. La Montaña, junto a una descripción de su origen, composición, correlaciones
estratigráficas y quiebres estructurales.
La Fm. Mininco es definida en este documento como sedimentitas lagunares y fluviales

periglaciales, constituidas mayoritariamente por areniscas medias finas y limolitas de
colores claros, las cuales sobreyacen en discordancia a las Lavas de Hueleheico.
Por otro lado, la Fm. La Montaña, constituyen los lomajes antepuestos a la Cordillera
Principal. Está constituida por una serie de depósitos morrénicos y aluviales, de composición
relativamente uniforme y con bancos lenticulares de arcillas y arenas finas.
• Mapas geológicos preliminares de Chile, Hoja Concepción – Chillán. Escala

1:250.000 (Instituto de Investigaciones Geológicas, 1981)
El área estudiada se ubica entre los 36° y 37° de latitud sur y al oeste de los 71°31’. Esta
carta presenta las unidades geológicas a una escala de 1:250.000 en el área comprendida
entre Concepción y Chillán. Se identifican los afloramientos geológicos más relevantes y
características morfoestructurales que se presentan en el extremo oeste de la cuenca de
Itata.
Estratigráficamente, las rocas más antiguas están representadas por rocas metamórficas
paleozoicas, que se distribuyen principalmente en el sector occidental de la cordillera de la
costa.
La unidades geológicas y sistemas estructurales presentes en este estudio permiten

correlacionar de mejor manera los espesores de la porción oeste de la cuenca en el presente
proyecto.
• Mapa Hidrogeológico Nacional (DGA, 1989)
Mapa hidrogeológico a nivel nacional, donde se identifica de forma general la importancia

acuífera de las unidades hidrogeológicas, y se relaciona estas a las distintas formaciones
geológicas presentes en Chile. Se señala rangos de productividad de pozos.
Para el área de estudio se identificaron principalmente unidades de alto potencial

hidrogeológico asociados a material que rellena la Depresión Central. Estas unidades se
componen mayoritariamente por sedimentos fluviales, aluviales, glaciales, lacustres,
7
aluvionales y eólicos, que conformarían acuíferos libres a semiconfinados de una extensión
variable.
Hacia los extremos este y oeste, la presencia de la Cordillera Principal y Cordillera de la

Costa repercute en un menor potencial hidrogeológico asociado a rocas volcánicas,
plutónicas, hipabísales y metamórficas.
• Mapa Geológico de Chile: Versión Digital (SERNAGEOMIN, 2003)
Mapa geológico de escala 1:1.000.000 que proporciona una representación conceptual de

las distintas unidades geológicas presentes en Chile. Se focaliza en agrupar las unidades
de acuerdo con los tiempos y procesos de mayor relevancia a una escala regional.
Adicionalmente cuenta con información de los principales dominios morfoestructurales a

nivel nacional.
Para la zona de estudio se destaca la presencia de secuencias sedimentarias y volcánicas

en la Depresión Central, y una predominancia de rocas intrusivas y metamórficas en los
extremos occidentales y orientales donde se encontrarían la Cordillera de la Costa y la
Cordillera Principal respectivamente.
• Estudio hidrogeológico cuencas del Biobío e Itata (DGA, 2011)
El objetivo general del estudio consistió en el levantamiento de información hidrogeológica

con el fin de desarrollar un modelo conceptual del acuífero de la cuenca del río Itata.
En el estudio se presenta una caracterización y un modelo conceptual a partir del análisis

de antecedentes y de derechos de agua, considerando también la medición de niveles de
209 pozos para la cuenca del Biobío, generando de información hidrogeológica preliminar.
Se realizó una estimación de la recarga para el acuífero de la cuenca del río Itata, a partir
de un modelo de simulación hidrológica simplificado, el cual se realizó considerando la
interacción existente entre los sectores de riego, los cauces superficiales, los embalses y
los pozos de bombeo.
Dentro de los resultados de interés, se destaca la campaña geofísica realizada donde se

ejecutaron un total de 353 estaciones gravimétricas (Figura 1.2), obteniendo 9 perfiles
8
longitudinales que muestran la profundidad de la roca. Se concluye que existe un solo gran
acuífero subdividido en tres zonas diferenciadas por el tipo de relleno permeable.
• Estudio hidrogeológico cuenca del Biobío (DGA, 2012)
El objetivo general del estudio fue el levantamiento de la información hidrogeológica y, a

partir de ello, el desarrollo de un modelo conceptual del acuífero de la cuenca del río Bío-
Bío.
Se recopiló la información de la geología y geomorfología del área de interés, estructurando

un mapa geológico a escala 1:250.000 que permitiera el trazado detallado del contacto
roca-relleno.
Los resultados de mayor interés incluyen los resultados del levantamiento gravimétrico
realizado en la cuenca del río Biobío, abarcando las regiones del Biobío y La Araucanía.
Mediante la ubicación de estaciones gravimétricas, se pudo proporcionan información
general sobre la morfología del basamento subterráneo del Valle Central en latitudes
aproximadas de 37º S a 38º S. Se llevaron a cabo un total de 15 perfiles gravimétricos,
con estaciones distribuidas aproximadamente cada 2 kilómetros.
Adicionalmente el estudio describe las unidades hidrogeológicas presentes,

categorizándolas según su importancia que va desde Alta a Nula, basándose en su potencial
hidrogeológico. Los resultados presentan una Unidad de Alta Importancia Hidrogeológica
comprende los depósitos no consolidados del cuaternario, incluyendo depósitos fluviales y
fluvio-aluviales, tanto actuales como antiguos, así como depósitos coluviales y sedimentos
de terrazas fluviales y terrazas marinas. La Unidad de Media a Baja Importancia
Hidrogeológica engloba sedimentos con un origen glacio-lacustre, caracterizados por tener
porcentajes de finos superiores al 30%. Estos sedimentos están vinculados genéticamente
con procesos de glaciación y deglaciación de masas de hielo situadas en elevaciones
mayores.
A partir del relevamiento de información hidrogeológica, DGA elabora un modelo conceptual

del acuífero presente en la cuenca del río Biobío. El proceso involucró un análisis exhaustivo
del catastro de demanda, la definición de unidades geológicas y acuíferas, la estimación del
flujo de aguas subterráneas, la identificación de zonas de recarga y descarga, la obtención
de parámetros hidrogeológicos del acuífero y la definición geométrica mediante geofísica,
con 15 perfiles gravimétricos.
9
• Modelación Hidrogeológica cuenca Itata - Bajo, Región del Biobío (DGA, 2013)
El principal objetivo del estudio corresponde a desarrollar una modelación conceptual y

numérica del acuífero de la cuenca del río Itata bajo, utilizando datos antiguos y propios
del levantamiento geofísico que realizaron.
La recopilación de antecedentes y levantamiento de información realizada en el estudio

tuvo como objetivo la incorporación de cada uno de los elementos necesarios para la
generación de un modelo conceptual y posteriormente, un modelo numérico.
En la caracterización hidrológica de detalle de la cuenca, se incluye la descripción de la red

de drenaje, los afluentes del río Itata, la definición de sus cuencas y subcuencas, y los
análisis de los datos fluviométricos y de precipitaciones brutas.
La campaña geofísica realizada utilizo el método de sondeos eléctricos verticales (SEV). La

interpretación de la campaña geofísica permitió asignar el comportamiento geoeléctrico a
distintos materiales presentes en la zona:
− Depósitos detríticos no consolidados cuaternarios: la resistividad eléctrica decrece

con la disminución del tamaño de grano medio de la litología, pudiéndose dar
valores, en ausencia de agua, desde pocas decenas de Ωm (arcillas y limos) hasta
varios cientos (arenas gruesas y gravas) o incluso algunos miles (bloque). La
presencia de agua y su salinidad condiciona, principalmente en el caso de agua
salada, de manera importante los valores resistividad eléctrica que presentan los
materiales detríticos. De esta forma y para una misma litología, ésta será menos
resistiva si tiene presencia de agua dulce y significativamente menos resistiva si
alberga agua salada.
− Depósitos plio-cuaternarios consolidados: en general presentan resistividades más
altas que los depósitos no consolidados.
− Basamento rocoso: Los valores de resistividad eléctrica esperados para el
basamento intrusivo (granodioritas y tonalitas) serán mucho más altos que el de los
materiales suprayacentes (detríticos no consolidados) mientras que el del
metamórfico (esquistos, filitas y metareniscas) será sensiblemente inferior.
• Levantamiento hidrogeológico en cuencas pluviales costeras en la Región del

Libertador Bernardo O'Higgins, del Maule y del Biobío, Etapa 2 (DGA, 2014)
10
El objetivo general de este estudio consistió en el desarrollo de modelos conceptuales de
las distintas cuencas y acuíferos estudiados, permitiendo generar una herramienta de
gestión basada en función de la hidrología e hidráulica, así como en la disponibilidad de
recursos hídricos.
Este estudio consistió en la modelación conceptual de 13 acuíferos, correspondientes a las

tres cuencas costeras de la región de O’Higgins (Río Rapel, Costeras entre estero Nilahue
y límite región y Costeras entre estero Topocalma y estero Nilahue), de las 6 cuencas
costeras de la región del Biobío (Costeras entre río Itata y río Pingueral, Costeras entre río
Pangue y río Paicaví, Costeras entre rio Biobío y río Manco, Río Biobío bajo, Río Paicaví y
Río Andalién) y de las tres comunas costeras de la Región del Maule (Constitución,
Empedrado y San Javier).
Se realizó una caracterización hidrológica e hidrogeológica, donde en cada una de las

unidades se efectuó un balance hídrico considerando las entradas y salidas a cada uno de
los sistemas.
Con el objetivo de servir de apoyo para la caracterización geológica y geométrica, que

permita una mejor elaboración de los modelos conceptuales de los acuíferos de las cuencas
de estudio, se llevó a cabo un estudio de prospección geofísica con los métodos transiente
electromagnético (TEM) y gravimétrico. Los resultados de las 45 estaciones TEM y 180
estaciones gravimétricas permitieron interpretar las unidades hidrogeológicas presentes
junto a la geometría de los potenciales acuíferos.
• Estimación Preliminar de las Recargas de Aguas Subterránea y Determinación de los

Sectores Hidrogeológicos de Aprovechamiento Común en la Cuenca de las Regiones
del Maule, Biobío, La Araucanía, Los Ríos y Los Lagos (DGA, 2014)
El objetivo principal de este informe busca delimitar y definir los distintos sectores
hidrogeológicos de aprovechamiento común de las cuencas de las regiones del Maule,
Biobío, La Araucanía, Los Ríos y Los Lagos.
En el presente estudio se realizó una sectorización de acuíferos y delimitación de 82 SHAC

de las cuencas de las regiones del Maule, Biobío, La Araucanía, Los Ríos y Los Lagos. Para
esto, se delimitaron las cuencas hidrográficas estudiadas, por medio de la definición de
unidades geomorfológicas como la Cordillera de los Andes, Depresión Intermedia y
Cordillera de la Costa para posteriormente, continuar con la subdivisión en unidades de
11
cuencas menores. Una vez realizada la delimitación de cuencas, se estimó la precipitación
media de éstas utilizando diferentes herramientas de cálculo, basándose en datos obtenidos
de las isoyetas del Balance Hídrico de Chile (DGA, 1987). Posterior al cálculo se
precipitaciones, se aplicó un coeficiente de infiltración a los valores obtenidos tomando
como base estudios previos de la DGA. De esta manera, se obtuvo una representación
preliminar de recarga media anual que se produciría en cada sector hidrogeológico
delimitado.
En la región del Biobío, se delimitan 18 sectores hidrogeológicos para ser estudiados, de

los cuales el Río Andalién corresponde al sector de interés. El estudio estima los factores
de infiltración, oferta de los recursos hídricos y demanda comprometida permitiendo
establecer la recarga sustentable de cada sector.
Se concluye que el sector hidrogeológico Río Andalién presenta una condición en la que la
explotación actual y el otorgamiento de nuevos derechos de aprovechamiento no generan
riesgo en el acuífero.
• Análisis y aplicación metodología para la delimitación y sectorización de acuíferos

en la provincia de Ñuble, VIII región (DGA, 2015)
El objetivo general de este estudio consistió en el análisis y determinación de la aplicabilidad

de la metodología de delimitación y sectorización de acuíferos en el acuífero de Ñuble,
cercano al estero Changaral, a partir de la elaboración de un análisis crítico y
complementación de la metodología de delimitación y sectorización de acuíferos a nivel
nacional desarrollada en el año 2014 por Guillermo Cabrera para la DGA.
En este estudio se establece una propuesta de sectorización para el acuífero asociado al río
Ñuble con 14 sectores. El trabajo presenta una gran cantidad de datos con respecto a la
caracterización hidrológica e hidrogeológica en el sector estudiado de la cuenca del río Itata,
generando un mapa con las unidades hidrogeológicas presentes en el área de estudio.
La campaña geofísica realizada incluye un total de total de 81 estaciones de TEM y 141

estaciones gravimétricas, de las cuales 81 estaciones gravimétricas fueron localizando en
una ubicación similar a la de las estaciones TEM. Las restantes estaciones gravimétricas se
ubicaron en un perfil longitudinal atravesando la Depresión Central en sentido O-E. Las
estaciones de TEM, en general, se ubicaron en perfiles de cuatro a seis y hasta ocho
12
estaciones TEM contiguas, abarcando longitudes de aproximadamente 1 km con estaciones
separadas cada ~200 m de cada punto central y loops individuales de 200x200 m.
Los resultados del estudio geofísico indican que el basamento alcanza predominantemente
profundidades menores a 300 m, con extremos de hasta 500 m en el área de estudio.
• Prospección geofísica sector rural déficit hídrico Valle del Itata (DGA, 2018)
El objetivo del estudio consistió en la definición de una serie de puntos de prospección para
la construcción de captaciones de agua, enmarcado dentro de un proyecto para paliar la
escasez hídrica que afecta la zona del valle del Itata.
Estudio geofísico en el cual se realizó una campaña de prospección geofísica con el método
del transiente electromagnético (TEM), con un total de 259 prospecciones (Figura 1.2) en
46 sectores distribuidas en 9 comunas del valle del Itata, con diferentes loop de prospección
(10x10 m, 25x25 m, 50x50 m y 100x100 m).
Se definieron siete unidades hidrogeológicas:
- Unidad I – Unidad saturada integrada por una secuencia de sedimentos finos con
presencia de arcillas.
- Unidad II – Unidad saturada constituida por arenas medias, arenas gruesas y

gravas.
- Unidad III – Unidad parcialmente saturada integrada por arenas medias, arenas
gruesas y gravas.
- Unidad IV – Unidad parcialmente saturada constituida por arenas medias y arenas

finas.
- Unidad V– Unidad con saturación baja a seca integrada por sedimentos finos.
- Unidad VI – Unidad constituida por sedimentos o roca fracturada.
- Unidad VII – Unidad constituida por sedimentos compactos o roca fracturada.
13
Además, en el documento se recopilaron los parámetros hidráulicos de los diferentes
sectores donde se realizaron las prospecciones (transmisividad y coeficiente de
almacenamiento), a partir de la consulta de los expedientes de la DGA.
• Sectorización Hidrogeológica del acuífero del Río Laja, Región del Biobío (DGA,
2018)
El objetivo general del estudio fue realizar una propuesta de sectorización hidrogeológica y
una estimación de la recarga (oferta) para el acuífero del valle del río Laja
En el estudio se determina la disponibilidad de aguas subterráneas susceptibles de

explotación a nivel del Sistema Hidrogeológico de Áreas de Captación (SHAC) en el valle
del río Laja.
Se llevó a cabo una sectorización hidrogeológica y una estimación de la recarga del acuífero,
con el objetivo de identificar los flujos de agua subterránea, así como las zonas de recarga
y descarga. Este proceso implico una recopilación de antecedentes, que abarca información
hidrológica, geológica e hidrogeológica, así como datos de captaciones registradas por la
Dirección General de Aguas.
Además, llevan a cabo análisis hidrogeoquímicos y una exploración geofísica, que incluye
la instalación de 55 estaciones de transiente electromagnético (TEM) y 50 estaciones
gravimétricas (Figura 1.2).
Los resultados obtenidos en base a los ensayos geofísicos ejecutados, define que el acuífero
del río Laja se emplaza en una depresión sedimentaria que mantiene una marcada
tendencia de profundización hacia el este. Se identificaron rellenos saturados de hasta 400
m de espesor, donde el basamento se ubica de manera más somera al oeste y flancos del
área de estudio. Asimismo, se distingue gran variedad granulométrica en la composición
de las unidades que conforman el subsuelo, desde bolones y ripios húmedos que presentan
medianas a altas resistividades, hasta materiales más arcillosos que arrojan menores
resistividades y representan las zonas menos permeables dentro del sistema acuífero. Se
confirma de esta manera la existencia de una formación acuífera de gran continuidad tanto
vertical como horizontal, de importante espesor, donde el basamento se hallaría por lo
menos a 1 km de profundidad desde la cota de terreno.
14
Finalmente, formulan un modelo conceptual del acuífero, proporcionando una base sólida
para la gestión sostenible de este recurso hídrico subterráneo.
• El Rol de la exploración geofísica en acuíferos profundos en ambientes semiurbanos

y rurales en cuencas andinas de ante arco: caso de estudio en acuífero del río Ñuble,
Valle Central de Chile (Aguirre et al., 2022).
El presente estudio propone una metodología para la detección y caracterización de

acuíferos profundos, utilizando una combinación de técnicas geofísicas. Proponen aplicar
métodos eléctricos, electromagnéticos y gravimétricos (Figura 1.2), complementados con
magnetometría, y considerando información hidrogeológica y propiedades físicas de la
cuenca.
Con el objetivo de determinar la profundidad hasta el basamento y el espesor del relleno

de la cuenca se utilizaron metodologías de resistividad eléctrica y gravimetría, las cuales,
además, permitieron establecer las variaciones de resistividad del subsuelo, y definir la
geometría de los acuíferos y las particularidades hidrogeológicas del área estudiada.
Para la interpretación geofísica se utilizaron 467 estaciones de gravedad, de las cuales 343
estaciones corresponden a mediciones realizadas por DGA (2011) y 124 estaciones de DGA
(2015), además de la información de estaciones TEM y MT que corresponden a datos
adquiridos por DGA (2015). El estudio realiza 81 estaciones TEM de loop 200 m (TEM-
200m), agrupadas en 16 perfiles de 5 estaciones cada uno, y 14 estaciones TEM de loop
50 m, ubicadas en torno a los sitios de los perfiles TEM-200m. Las 5 estaciones MT medidas
en ese mismo estudio se ubican en la parte rural-semiurbana de la cuenca del Ñuble (parte
central y oriental), dada la imposibilidad de medir en la parte urbana por dificultades de
ruido cultural. La metodología propuesta tiene el potencial de identificar un recurso vital
para la sustentabilidad en las próximas décadas.
La estimación regional de la profundidad de la cuenca se realizó en conjunto con un estudio

integrado TEM y MT en puntos de control. Estos han sido identificados para aquellos puntos
de medición TEM o MT en los cuales se reconoce un incremento significativo de la
resistividad. De esta forma fue posible definir la profundidad hasta el basamento en estos
puntos de control, lo que permitió calibrar el campo regional del estudio gravimétrico.
Los resultados permitieron caracterizar el acuífero entre los 50 y 500 m de profundidad,

como parte de una cuenca sedimentaria de más de 1.000 m de espesor.
15
Adicionalmente se generó un perfil hidrogeológico con orientación SO-NE en la que se
sintetiza la geometría del acuífero y las características hidrogeológicas del área estudiada.
• Análisis y caracterización del Acuífero de la cuenca del rio Itata, Estudio geofísico
Itata Bajo (DGA, 2023)
Estudio con el objetivo de generar nueva información geofísica en la zona acuífera de la

cuenca del río Itata bajo que permita caracterizar en profundidad el relleno sedimentario y
el basamento rocoso.
Se recopilo información bibliográfica referente a la geomorfología, geología e hidrogeología

del área de estudio, con esta información, se realizaron mediciones geofísicas utilizando el
método de TEM y ERT.
Con los resultados de los modelos de resistividad eléctrica obtenidos se interpretaron 8

unidades hidrogeológicas para el área de estudio:
- Unidad A: Capa somera, resistividad variable debido a los sedimentos y humedad.
- Unidad B1: Resistividades entre 10 a 300 Ωm, se asocia a depósitos de origen

sedimentario fluvioaluvial generalmente saturados.
- Unidad B2: Resistividades, en general, menores a 10 Ωm, unidades saturadas con

mayor proporción de material fino, hacia la costa se asocia a mayor salinidad.
- Unidad AB: Resistividades entre 20 a 80 Ωm. Asociables a la unidad B1 donde no se

reconoce la unidad A.
- Unidad BD: Unidad basal, resistividades variables en torno a 100 Ωm y sobre 2000
Ωm en algunos sectores. Se asocia a rocas metamórficas alteradas con un mayor
grado de saturación, salinidad y/o alteración en niveles con menores valores de
resistividad.
- Unidad BC: Resistividad entre 40 a 350 Ωm. Se identificaron en el sector de

Confluencia y se asociaron a secuencias sedimentarias y volcanosedimentarias.
- Unidad D: Unidad asociada al basamento rocoso con resistividades mayores a 350

Ωm.
16
Figura 1.2: Localización de estaciones geofísicas en estudios previos en las regiones del Ñuble y Biobío, relevantes
para este estudio. Fuente: Elaboración propia.
17
1.4. Conclusiones del análisis de antecedentes
El análisis de los antecedentes recopilados en el subcapítulo anterior permitió identificar las

unidades geológicas presentes en el área, junto a determinar la estratigrafía de pozos
profundos, la cual es de suma importancia para realizar una correcta interpretación de las
unidades geológicas e hidrogeológicas presentes.
Por otra parte, la búsqueda de estudios DGA que incluyen labores geofísicas previas en la
cuenca, permitió planificar de mejor forma las localidades donde se realizaron las
mediciones geofísicas del presente informe, además de darle una mayor robustez a los
modelos generados.
Para cada estudio revisado en el subcapítulo anterior se generaron fichas resúmenes, las
cuales se presentan en el Anexo H, cuya estructura base se presenta a continuación:
Figura 1.3: Estructura de las fichas bibliográficas presentadas en el Anexo H
18
Los estudios recopilados de mayor relevancia corresponden a:
• Estudio hidrogeológico cuencas Biobío e Itata (DGA, 2011).
• Estudio hidrogeológico cuenca Biobío (DGA, 2012).
• Prospección geofísica sector rural déficit hídrico valle del Itata (DGA, 2018).
• El Rol de la exploración geofísica en acuíferos profundos en ambientes semiurbanos

y rurales en cuencas andinas de ante arco: caso de estudio en acuífero del río Ñuble,
Valle central de Chile (Aguirre et al., 2022).
Los antecedentes del estudio DGA (2011) presentan una descripción geológica del área de
interés donde se indican las relaciones estratigráficas de las unidades junto a sus
composiciones y espesores aproximados. En términos geofísicos presentan modelos
gravimétricos de profundidad de basamento de algunos perfiles a lo largo de la cuenca, que
permitieron comprender la geometría general de la cuenca y con esto definir zonas de
interés para el presente estudio.
De manera similar, el estudio DGA (2012) presenta un plano de la geología presente al sur
de la zona de estudio y en términos geofísicos, presentan un modelo de profundidad de la
cuenca sedimentaria desde Bulnes hacia el sur.
En el estudio DGA (2018) se realizó una campaña de prospección geofísica con el método
del transiente electromagnético (TEM), con un total de 259 estaciones en 46 sectores
distribuidos en 9 comunas del valle del Itata. De estas mediciones, las de mayor
importancia corresponden a las emplazadas al este de su sector de estudio, debido a su
cercanía a la Cordillera de la Costa. La interpretación de los valores de resistividad y
unidades identificadas en este sector permitieron constreñir de mejor manera los modelos
geofísicos generados en este trabajo.
El estudio de Aguirre et al. (2022) complementa las mediciones geofísicas de DGA (2011)
y DGA (2015) con nuevas estaciones de gravimetría, TEM y MT, ubicadas en sectores
específicos para el desarrollo de un modelo mejorado de la profundidad de basamento y
espesor sedimentario de la cuenca. Este trabajo, aportó con puntos de control importantes
de profundidad de basamento en algunos sectores (principalmente de estaciones MT), para
el modelamiento de los perfiles gravimétricos del presente trabajo.
19
2. CONTEXTO GEOLÓGICO E
HIDROGEOLÓGICO
La identificación y caracterización de las rocas presentes en el área de estudio es

fundamental para conocer y entender la interacción del agua con el subsuelo, ya que la
formación de acuíferos está estrechamente relacionada con las características químicas,
físicas e hidráulicas de las rocas y los sedimentos. En este capítulo se indica y detalla la
información geológica e hidrogeológica recopilada, con el objetivo de entender y
caracterizar las unidades hidrogeológicas presentes.
2.1. Unidades litológicas
La geología del sector (Figura 2.1) indica un predominio por material sedimentario que se
encuentra rellenando valles y de manera subordinada depósitos volcánicos y
volcanosedimentarios hacia el extremo oriental y occidental, donde geomorfológicamente
se emplazarían la Cordillera Principal y la Cordillera de la Costa respectivamente.
La unidad más antigua corresponde a una franja intrusiva de edad Carbonífero – Pérmico
(CPg en Figura 2.1) con orientación NNE-SSO, la cual aflora hacia el oeste del área de
estudio. Estas rocas conforman gran parte de la Cordillera de la Costa y su composición
petrográfica varía entre tonalitas y granodioritas.
Los depósitos de mayor extensión en la cuenca corresponden a los de Formación Mininco

(PPlm en Figura 2.1) de edad Plioceno – Cuaternario, y comprenden sedimentitas lagunares
y fluviales en forma de areniscas tobáceas, limonitas y conglomerados con intercalaciones
de arcillolita y tobas.
20
Figura 2.1: Mapa geológico de la Cuenca del Río Itata. Modificado de DGA, 2011.
21
La Formación La Montaña (PlHlm en Figura 2.1), de edad Pleistoceno – Holoceno, aflora
principalmente al extremo noreste del área de estudio y está conformada por una serie de
depósitos morrénicos y aluviales compuestos principalmente por arcillas y arenas finas con
poca compactación y sin cementación. Estudios más recientes, al sur del área de estudio,
determinan que la parte superior del relleno sedimentario del Valle Central correspondería
principalmente a la Formación Malleco (antes definida como Fm. La Montaña), ya que no
se reconoce diferencias en las facies que permitan diferenciar estas dos unidades
(Anavalón, 2018). Fm. Malleco estaría constituida por una sucesión de conglomerados,
areniscas y en menor medida, limolitas y tobas, las cuales estarían cubriendo a Fm. Mininco.
Los Depósitos piroclásticos (Pl3t en Figura 2.1) están en el sureste del área y corresponden
a depósitos no consolidados de edad cuaternaria, asociados a flujos piroclásticos de
composición riolítica por el colapso de calderas volcánicas.
Los Depósitos de avalancha volcánica (Q3av en Figura 2.1) se encuentran al suroeste del
área y de manera similar a los depósitos descritos anteriormente, se relacionan a colapsos
parciales de edificios volcánicos. Presenta clastos de composición riolítica a andesítica con
tamaños que varían entre grava gruesa y sedimentos finos.
Los Depósitos fluvio–aluviales (Q1fa en Figura 2.1) corresponden a materiales no

consolidados, compuestos por gravas, ripios y arenas que se encuentran rellenando cauces
de ríos y esteros. Estos depósitos se asocian a paleocauces presentes en el área, a
diferencia de los Depósitos actuales (Q1) que abarcarían el material de los cauces actuales.
2.2. Unidades hidrogeológicas
Según DGA (2015), hidrogeológicamente la cuenca de estudio presenta 6 unidades

hidrogeológicas (Figura 2.2), con un mayor predominio de la Unidad B al centro del valle.
La Unidad A se extiende a lo largo de cauces de ríos y esteros, y las unidades C1 y D hacia
la porción este del área donde se emplaza la Cordillera Principal.
La Unidad A se asocia a acuíferos de alta a mediana importancia hidrogeológica en depósitos

no consolidados. Acuíferos intergranulares, productivos, amplia a medianamente
extendidos en forma longitudinal en los depósitos fluviales y fluvio-aluviales en las unidades
geológicas Q1 y Q1fa.
22
Figura 2.2: Mapa hidrogeológico de la Cuenca del Río Itata. Modificado de DGA, 2015.
23
La Unidad B corresponde a acuíferos de mediana importancia hidrogeológica en rocas o
depósitos no consolidados que forman acuíferos de amplia distribución con intercalaciones
de litologías de menor permeabilidad, asociado principalmente a las formaciones geológicas
PPlm y PlHlm.
La Unidad C1 conforma acuíferos de baja importancia hidrogeológica en rocas fisuradas. A

pesar de que no se conocen muchos antecedentes, se deduce de las características
geológicas de esta unidad podría conformar acuíferos locales y con productividad irregular.
Esta unidad se correlaciona con las secuencias volcano-sedimentarias OM2c.
La Unidad D presentaría condiciones de acuicludos y/o acuífugos de nula a muy baja

importancia hidrogeológica con rocas esencialmente sin recursos de agua subterránea. Se
asocia a unidades intrusivas y volcánicas como Kiag, Mg, PPl3.
Las unidades hidrogeológicas que predominan al interior del área a sectorizar corresponden
a las Unidades A, B, C1 y D. Respecto de estas dos últimas, no se conocen antecedentes
hidrogeológicos que den cuenta de la presencia de acuíferos relevantes. En esta área, las
Unidades C1 y D se localizan en zonas elevadas o de mayor cota respecto de las Unidades
A y B, por lo que, se asume que corresponden a áreas de baja relevancia para la
sectorización acuífera.
El basamento de la cuenca está conformado por las Unidades C, C1 y D, los cuales con los
antecedentes de DGA (2015) se consideraron impermeables.
2.3. Captaciones profundas en la cuenca del Itata
Para una mejor caracterización de las unidades hidrogeológicas presentes en la cuenca del
río Itata se recurrió a la búsqueda en la base de datos de la Dirección General de Aguas de
captaciones profundas cercanas. Debido a la densidad de las captaciones presentes en el
área (Figura 2.3), se revisó y analizó las captaciones en distintas áreas de interés, donde
se planifica, realizar las campañas geofísicas con el fin de identificar pozos profundos que
alcanzan a detectar roca y así poder ajustar los modelos geofísicos.
24
Figura 2.3: Mapa de captaciones subterráneas registradas en la Dirección General de Aguas Subterráneas en la
Cuenca del Río Itata. Elaboración Propia
25
Se recopilaron y analizaron un total de 856 captaciones (ANEXO A), siendo la mayoría de
estas norias, drenajes y pozos zanjas de poca profundidad. Sin embargo, se han
identificado 24 pozos profundos que resultan fundamentales para el ajuste y la corrección
de los modelos geofísicos (Tabla 2.1). De estos pozos, se han realizado columnas
estratigráficas (Figura 2.4 a Figura 2.9) que muestran la presencia de roca del basamento
en profundidad. Además, en la Figura 2.10 se observa la relación de estos datos con los
perfiles gravimétricos realizados.
Tabla 2.1: Tabla de información de los pozos profundos recopilados que detectan roca en la cuenca
del Rio Itata. Base cartográfica en WGS84 H18S.
ID Código de Profundidad Profundidad

Tipo UTM E (m) UTM N (m)
Expediente (m) Roca (m b.n.t.)
1 ND-1601-16 Pozo 730.605 5.966.483 60 22
2 ND-0801-1102 Pozo 732.350 5.911.700 58 31
3 ND-0801-1255 Pozo 738.632 5.921.319 76 70
4 ND-0801-1917 Pozo 739.434 5.912.904 40 24
5 ND-0801-1928 Pozo 730.524 5.930.098 40 25
6 ND-0801-2278 Pozo 784.548 5.913.436 104 82
7 ND-0801-2312 Pozo 739.061 5.960.133 87 39
8 ND-0801-5665 Pozo 734.720 5.920.941 100 28
9 ND-0801-6272 Pozo 734.320 5.920.641 40 36
10 ND-0801-6680 Pozo 728.568 5.941.406 42 31
11 ND-0801-7125 Pozo 742.348 5.937.766 84 84
12 ND-0801-8231 Pozo 724.955 5.945.197 62 53
13 ND-0802-3678 Pozo 777.485 5.876.327 63,5 30
14 ND-0802-6154 Pozo 718.439 5.901.456 60 30
15 ND-0802-6759 Pozo 794.248 5.860.328 50 40
16 ND-1601-43 Pozo 725.488 5.944.081 44 43
17 ND-1602-125 Pozo 790.840 5.921.179 104 98
18 ND-1602-299 Pozo 741.455 5.941.160 38 36
19 ND-1602-327 Pozo 724.143 5.933.282 24 23
20 ND-1602-440 Pozo 746.175 5.941.642 53 24
21 ND-1602-521 Pozo 743.567 5.937.821 74 72
22 ND-1603-319 Pozo 750.580 5.954.964 151 150
23 ND-1603-811 Pozo 809.517 5.950.104 100 48
24 ND-1602-354 Pozo 724.496 5.932.842 42 17
26
Figura 2.4: Columnas estratigráficas de las captaciones que detectan roca (Perfiles gravimétricos L1
y L2) ubicadas en la cuenca del Río Itata.
27
Figura 2.5: Columnas estratigráficas de las captaciones que detectan roca (Perfil gravimétrico L3) ubicadas en la cuenca del Río Itata.
28
29
30
Figura 2.8: Columnas estratigráficas de las captaciones que detectan roca (Perfiles gravimétricos L6 y L7) ubicadas en la cuenca del Río
Itata.
31
Figura 2.9: Columnas estratigráficas de las captaciones que detectan roca (Perfiles gravimétricos D7A, D10 y D9B) ubicadas en la cuenca del
Río Itata.
32
Figura 2.10: Plano de ubicación de los perfiles elaborados a partir de las captaciones seleccionadas.
33
3. METODOLOGÍA
3.1. Gravimetría
El método gravimétrico es una técnica de exploración geofísica que utiliza mediciones de la

componente vertical de la gravedad en la superficie de la Tierra para determinar la densidad
de las rocas presentes en el subsuelo.
Esta metodología es una técnica pasiva (remota) que consiste en determinar las diferencias
de gravedad entre distintos puntos distribuidos en la zona de estudio. Las variaciones del
campo gravitatorio son generadas por la presencia de cuerpos de distinta densidad en
profundidad. Dado esto, debido a las diferencias de densidad entre los rellenos
sedimentarios de una cuenca y su roca basal, la interpretación de las anomalías
gravimétricas son una gran herramienta para el estudio y modelamiento de cuencas
sedimentarias.
La adquisición de datos consiste en la medición simultánea de la componente vertical del

campo gravitatorio, utilizando un instrumento llamado gravímetro, y la posición exacta del
punto de medición mediante un GPS diferencial, que permite obtener una resolución
centimétrica en la componente vertical. En general, para este tipo de estudios se utilizan
gravímetros relativos, por lo que durante la campaña de terreno es necesario medir en
algún punto de referencia donde exista una medición absoluta de gravedad. Además,
también es necesaria la medición diaria en una estación base, al comenzar y finalizar el
día, para poder realizar lo que se conoce como corrección por deriva instrumental, que se
detalla en la sección de procesamiento.
El flujo de adquisición, procesamiento y modelamiento de datos gravimétricos consiste en

distintas etapas que son detalladas a continuación (Figura 3.1).
34
Figura 3.1: Flujo de adquisición, procesamiento y modelamiento de datos gravimétricos.
3.1.1. Terreno
En base a los requerimientos y objetivos del estudio, se realizó la exploración geofísica

mediante el método gravimétrico, midiendo un total de 535 puntos de gravedad relativa
en la zona de estudio (Figura 3.2), de los cuales se utilizaron 519 para la creación y
modelamiento de 24 perfiles: 7 perfiles principales en dirección E-O, que cruzan el valle
central desde la Cordillera de la Costa a la Cordillera Principal, y 17 perfiles auxiliares, en
distintas direcciones. En total, se logró obtener alrededor de 1.194 km de perfiles lineales.
Las mediciones se realizaron con un gravímetro relativo G-370 marca Lacoste & Romberg.
Los datos fueron adquiridos bajo un protocolo de loop diario mediante la medición de una
estación base en distintos momentos del día, para poder hacer la corrección por deriva
instrumental. Y, para la correcta conversión de gravedad relativa a gravedad absoluta, se
utilizó como estación base de gravedad absoluta la estación CHILLAN, medida por el
Instituto Geográfico Militar. En cuanto a los datos topográficos se realizaron con un GPS
diferencial marca Sanding, modelo GNSS RTK Sanding T7, que fueron procesados en el
software SurvX 4.0.
35
Figura 3.2: Mapa de ubicación de estaciones gravimétricas medidas y estación de gravedad absoluta utilizada para procesamiento.
36
3.1.2. Procesamiento
El procesamiento de datos gravimétricos se separa en 2 grandes etapas: las correcciones

dinámicas y el cálculo de la anomalía de Bouguer completa. Todo este procesamiento se
lleva a cabo mediante el software Oasis Montaj de Geosoft.
3.1.2.1. Corrección de variaciones dinámicas
• Corrección por marea: Corresponde a la corrección de las variaciones de gravedad

en un punto debido a la posición relativa del Sol y la Luna respecto a la Tierra. Para
su cálculo, se utilizó el algoritmo matemático de Longman et al. (1959), utilizando
la fecha, ubicación y hora exacta de cada medición.
• Corrección por deriva instrumental: La deriva instrumental corresponde a las

variaciones en las mediciones de gravedad debido a errores instrumentales,
principalmente por el desgaste y/o uso de los equipos de medición. Esta deriva se
suele corregir diariamente y se aproxima a una relación lineal en el tiempo. Se
corrige realizando más de una medición de gravedad en un mismo punto durante el
día de medición, con lo que se estima la relación de variación lineal diaria que
permite calcular la deriva para cada momento del día.
• Cálculo de gravedad absoluta: Dado que durante las campañas gravimétricas se

utilizan gravímetros relativos, es necesario hacer una medición de un punto de
referencia en que la gravedad absoluta sea conocida. De esta manera, conociendo
la diferencia relativa entre un punto de la campaña y el punto de gravedad absoluta,
es posible calcular la gravedad absoluta para todos los puntos de medición (Figura
3.3). Para este estudio se utilizó como referencia la estación CHILLÁN, ubicada en
la Plaza de Armas de la ciudad de Chillan, exactamente en las coordenadas 36º 36’
24” S, 72º 06’ 12” O, que posee un valor de gravedad conocido de 979.815 mGal
(ANEXO D).
37
Figura 3.3: Gravedad Absoluta calculada para la zona de estudio.
38
3.1.2.2. Cálculo de anomalía de Bouguer Completa
Una vez obtenida la gravedad absoluta para cada punto de medición, se realiza el cálculo
de la anomalía de Bouguer Completa, según el proceso de corrección estándar de datos
gravimétricos (Blakely, 1995; Lowrie, 2007) que se presenta a continuación.
❖ Gravedad normal (𝑔0 ): Corresponde a la gravedad teórica sobre el elipsoide

equipotencial de referencia, que depende de la latitud geográfica. Está definida
mediante la siguiente relación:
1 + 𝑘 𝑠𝑖𝑛2 𝜑
𝑔0 = 𝑔𝑒 .
√1 + ⅇ 2 𝑠𝑖𝑛2 𝜑
Donde 𝜑 corresponde a la latitud, 𝑔𝑒 a la gravedad en el ecuador y k, e son valores

adimensionales que dependen de 𝑎, 𝑏, 𝑔𝑒 y 𝑔𝑝 (gravedad en el polo). Para este estudio se
utilizan las constantes definidas según el sistema de referencia WGS84 (World Geodetic
System 1984), detalladas en Blakely (1995).
❖ Aire Libre (ⅆg A ): Corresponde a la variación de gravedad dada por la diferencia de

altura entre el punto P de medición en superficie y el datum P’ (Figura 3.4).
Figura 3.4: Figura esquemática de la corrección por altura o aire libre.
Se estima de la siguiente manera (Lowrie, 2007):
ⅆ𝑔𝐴 [𝑚𝐺𝑎𝑙] ≈ 0.3086 ∙ ℎ [𝑚]
Donde h es la elevación del punto de medición, medida mediante un GPS diferencial.
39
❖ Bouguer Simple (ⅆ𝑔𝐵 ): Corresponde a la variación de gravedad generada por la
masa que se encuentra entre el datum y el punto de medición (Figura 3.5).
Figura 3.5: Figura esquemática de la corrección de Bouguer o de placa.
Se estima de la siguiente manera:
𝑔
ⅆ𝑔𝐵 [𝑚𝐺𝑎𝑙] ≈ 0.0419 ∙ 𝜌[ ] ∙ ℎ [𝑚]
𝑐𝑐
Donde h es la elevación del punto de medición y en este estudio se considera una densidad
de reducción clásica de 2,67 kg/m3.
❖ Corrección topográfica (ⅆ𝑔𝑇 ): Corresponde a la variación de gravedad por el efecto

gravitatorio que genera la topografía de la zona de estudio (Figura 3.6).
Figura 3.6: Figura esquemática de la corrección topográfica.
Esta corrección fue calculada siguiendo la combinación de algoritmos propuestos por Kane
(1962) y Nagy (1966) y una grilla topográfica ASTER GDEM, de resolución de 30mx30m
(Figura 3.7).
40
❖ Cálculo de Anomalía de Bouguer Completa (ABC): Una vez calculadas todas las
correcciones para cada punto de medición, se calcula la anomalía de Bouguer
completa (Figura 3.8) de la siguiente manera:
𝐴𝐵𝐶 = 𝑔𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑎 − 𝑔𝑜 + ⅆ𝑔𝐴 − ⅆ𝑔𝐵 + ⅆ𝑔𝑇
• Cálculo de anomalía de Bouguer Completa residual (extracción de regional): Las

anomalías gravitatorias se ven afectadas por las distintas dimensiones,
profundidades y contrastes de densidad de los cuerpos anómalos presenten bajo
superficie. En general un mapa de anomalías de Bouguer contiene anomalías
superpuestas de distintas fuentes, donde las anomalías de larga longitud de onda,
asociadas a cuerpos grandes y profundos, son denominadas anomalías regionales,
y las de longitud de onda corta, asociada a cuerpos más someros y pequeños que
suelen ser los de interés, se denominan anomalías residuales. Por esta razón, la
separación de las anomalías regionales de las anomalías de interés es fundamental
para el análisis gravimétrico.
En Chile, la tendencia regional de gran escala se debe principalmente a la diferencia
de densidad entre las raíces corticales de la Cordillera Principal y la roca basal. Para
su análisis y determinación existen diversas técnicas de descomposición, dentro de
ellas el análisis visual, representación polinomial, filtros y el uso de puntos
referenciales. En este estudio, a través de puntos de referencia de gravedad
medidos sobre afloramientos de roca a ambos extremos de la zona de estudio, se
determina un plano de tendencia regional representativo de la zona (Figura 3.9).
Luego, para separar esta tendencia de gran escala de los datos observados, basta
con hacer una resta directa entre la ABC y la tendencia regional calculada, con lo
que se obtiene la anomalía de Bouguer completa residual (Figura 3.10), mapa a
partir del cual se puede entender en términos generales la estructura de la cuenca
sedimentaria.
41
Figura 3.7: Modelo de elevación topográfica (DEM) utilizado para corrección topográfica.
42
Figura 3.8: Anomalía de Bouguer Completa (ABC) calculada para la zona de estudio.
43
Figura 3.9: Campo regional calculado para la zona de estudio. En diamantes negros los puntos de
medición de gravedad sobre roca.
44
Figura 3.10: Anomalía de Bouguer Completa reducida calculada para la zona de estudio.
45
3.1.3. Modelamiento
El modelamiento de los datos gravimétricos se realiza mediante modelamiento directo

(forward) 2D utilizando Model Vision de Encom basado en el algoritmo de Won y Bevis
(1987), el cual calcula la señal gravimétrica de un polígono de n lados con una densidad.
Los cuerpos construidos son de extensión limitada y simétricos en la tercera coordenada
(Figura 3.11).
Figura 3.11: Ejemplo de cuerpo en 2D (Shuey y Pasquale 1973).
Considerando que la finalidad de este estudio es conocer la geometría de la cuenca

sedimentaria de la zona de estudio, los perfiles se modelaron considerando un modelo
simple de relleno sedimentario de baja densidad, sobre una unidad de roca basal con una
densidad mayor. Para esto, todos los datos gravimétricos han sido medidos y dispuestos
en perfiles como se muestra en la Figura 3.12. De esta manera es posible modelar perfil a
perfil de forma bidimensional y determinar la profundidad de basamento bajo cada estación
gravimétrica.
46
Figura 3.12: Mapa de ubicación de perfiles gravimétricos modelados.
47
❖ Limpieza de datos
El procesamiento de los datos de cada perfil comienza con un control de calidad de la

anomalía de Bouguer Completa, eliminando todos aquellos puntos que se consideren
outlayer o fuera de tendencia respecto al resto del perfil. En la Figura 3.13 se muestra un
ejemplo de un punto gravimétrico (G304) que se muestra fuera de tendencia dentro del
perfil L7, que tuvo que ser eliminando para un correcto modelamiento.
Figura 3.13: Ejemplo de dato outlayer, perfil L7, punto de medición 304.
❖ Cálculo de tendencia regional individual para cada perfil e integración de información
Para hacer el modelamiento de cada perfil se utilizó la anomalía de Bouguer completa y la

tendencia regional se calculó de manera individual y específica para cada uno de los perfiles
dentro del software de modelamiento. Para esto último, se integró toda la información de
profundidad de basamento en puntos cercanos al perfil mediante gravimetría y otras
metodologías (MT, TEM, estratigrafía de pozos, estudios anteriores, etc.), permitiendo
obtener distintos puntos de control a lo largo de la línea. Esta información se incluye en el
modelo, constriñendo específicamente la profundidad de basamento entorno a estos
puntos, y por lo tanto permite obtener la ubicación de la tendencia regional en ese sector.
El modelamiento se realizó primero para los perfiles principales (L1-L7), y la profundidad
de basamento de estos fue utilizada también para constreñir los modelos de los perfiles
auxiliares (D1-D14). En la Figura 3.14, se muestra el ejemplo del modelamiento del perfil
L1, donde se pueden observar todos los puntos de profundidad de basamento utilizados
48
(líneas verticales moradas) y los puntos de control de la tendencia regional (puntos
rosados).
❖ Modelamiento 2D
Como se mencionó anteriormente, en este estudio se realizaron modelamientos forward 2D

de relleno sedimentario simple. Para esto se asume una roca basal (background) de 2,67
kg/m3, densidad promedio del basamento rocoso continental (Blakely, 1996) y se modela
el espesor de la cuenca sedimentaria asumiendo una densidad única de relleno
sedimentario, que en general puede variar entre 1,8 y 2,2 kg/m3. La profundidad de relleno
sedimentario del modelo que se obtiene depende directamente del contraste de densidad
entre basamento y sedimento, lo cual es muy importante de considerar en la interpretación
de los datos. El relleno sedimentario de una cuenca de escala regional no suele poseer
densidades laterales homogéneas ni tampoco en profundidad, dado que las unidades
sedimentarias varían a lo largo de la cuenca y, además, la compactación con la profundidad
juega un rol importante en el aumento de la densidad.
Para este estudio se realizaron pruebas de modelamiento forward del perfil L1 con distintas
densidades típicas de relleno sedimentario (1,9, 2,0 y 2,1 kg/m3) que corresponden a
contrastes de 0,77, 0,67 y 0,57 kg/m3 respectivamente, con el objetivo de poder escoger
el valor más adecuado para modelar la cuenca regional y poder analizar la variabilidad de
los espesores sedimentarios en función de este parámetro. Como se puede observar en la
Figura 3.15, los modelos realizados con las distintas densidades ajustan muy bien los datos
observados, pero existe una variación importante en el espesor de la cuenca sedimentaria
obtenida. No obstante, la forma de la cuenca se conserva y es muy similar en todos los
casos. En base a este análisis, se respalda la importancia de considerar las posibles
variaciones en los contrastes de densidad durante la interpretación, a pesar de que los
modelos restantes son modelados con una densidad única de 2,0 kg/m3 (contraste de 0,67
kg/m3), escogida como un valor promedio coherente, que se ajusta también a información
de otros estudios respecto a la profundidad de la cuenca sedimentaria en la región.
49
Figura 3.14: Ejemplo de modelamiento e integración de información de otras metodologías para el perfil L1. Elaboración propia.
50
Figura 3.15: Modelación del Perfil L1 con distintas densidades de sedimentos (1,9; 2,0 y 2,1 kg/m3). Elaboración propia.
51
3.2. Magnetotelúrica
El método magnetotelúrico (MT) es una técnica de exploración geofísica perteneciente a los

métodos electromagnéticos pasivos que utiliza un amplio espectro de variaciones
geomagnéticas naturales como fuente de energía, tales como tormentas eléctricas,
corrientes ionosféricas y actividad solar. Esto permite la obtención de información sobre la
distribución de la resistividad eléctrica en el subsuelo.
La adquisición de datos consiste en la medición de forma simultánea de las variaciones

temporales de los campos eléctricos y magnéticos en direcciones ortogonales en la
superficie de la Tierra (Figura 3.16). El campo eléctrico es medido a través de dos dipolos
eléctricos, formados por electrodos impolarizables, distribuidos en las direcciones norte-sur
(Ex) y este-oeste (Ey), respectivamente. El campo magnético es medido a través de tres
bobinas de inducción, orientadas en las direcciones norte (Hx), este (Hy) y vertical (Hz).
Figura 3.16: Instalación de una estación de magnetotelúrica. Obtenida de Metronix (2017).
El flujo de adquisición, procesamiento y modelamiento de datos MT consiste en distintas

etapas que son detalladas a continuación (Figura 3.17).
52
Figura 3.17: Flujo de adquisición, procesamiento y modelamiento de datos MT.
3.2.1. Terreno
❖ Adquisición
En base a los requerimientos y objetivos del estudio, se realizó la exploración geofísica de

Magnetotelúrica (MT), realizando en esta campaña un total de 100 estaciones de medición
(Figura 3.18). Las estaciones se distribuyeron con énfasis en sectores cercanos a la
Cordillera de la Costa y a la Cordillera Principal, donde la posibilidad de observar presencia
de basamento mediante esta metodología es mayor, y considerando también las zonas con
altos gravimétricos donde es importante tener más información para el modelamiento de
la cuenca sedimentaria.
53
Figura 3.18: Mapa de distribución de estaciones MT.
54
Para llevar a cabo la adquisición de datos, fue utilizado el equipo de medición ADU08e, de
la compañía Metronix (para más información referirse a Metronix (2017)). Este dispositivo
está equipado con tres bobinas de inducción, modelo MF07e, que se instalan en terreno
con orientación norte-sur, oeste-este y vertical, respectivamente, con el propósito de
registrar las variaciones naturales en el tiempo del campo magnético. Además, cuentan con
cuatro electrodos impolarizables EFP06, que se instalan distribuidos en dos dipolos de
longitud entre 10 m y 80 m, con orientación norte-sur y este-oeste, y se utilizan para
capturar información sobre las fluctuaciones temporales del campo eléctrico terrestre.
Todos estos sensores se conectan a la unidad de adquisición (ADU08e), que almacenan
simultáneamente las series temporales de los campos eléctricos y magnéticos medidos. El
esquema de medición de las estaciones de AMT para este estudio se resume en la Tabla
3.1.
Tabla 3.1: Esquema de adquisición de las estaciones de magnetotelúrica
Frecuencia de muestreo (Hz) Tiempo de medición (s)
65.536 90
16.384 240
4.096 360
1.024 480
512* 2.700 – 18.000
* No todas las estaciones de AMT fueron adquiridas con la banda de 512 Hz.
❖ Control de calidad en terreno
Durante la adquisición, es posible realizar un control de calidad de los datos analizando

tanto las series de tiempo, como las funciones de transferencia, las cuales son
determinadas preliminarmente en terreno mediante el software ProcMT, proporcionado por
Metronix. Esto posibilita el análisis la calidad de datos obtenidos para las diferentes
estaciones in situ, permitiendo repetir algunas mediciones en caso de considerarse
necesario.
55
❖ Control de calidad de series de tiempo
El procesamiento de datos comienza con la inspección visual de las series de tiempo

registradas por los cinco sensores desplegados en terreno (Ex, Ey, Hx, Hy y Hz) (Figura
3.19). Para esto se utiliza el módulo Tsplotter del software ProcMT. En esta etapa es posible
identificar el ruido antrópico existente en las series de tiempo u otros problemas de
medición, pudiéndose excluir del procesamiento los tramos con una menor razón señal a
ruido de la serie de tiempo.
Figura 3.19: Series de tiempo de una estación de magnetotelúrica (TRAP5).
❖ Estimación de funciones de transferencia
Las series temporales son transformadas desde el dominio temporal al dominio de la

frecuencia usando la Transformada de Fourier. Esto se realiza dividiendo la serie en
segmentos de igual tamaño (denominado Ancho de Ventana de Fourier), para luego
multiplicar cada segmento por un tipo de función ventana definido (e.g.: Hamming,
Blackman Rectangular, Gaussian, etc.), y posteriormente aplicar la Transformada de
Fourier para cada segmento. Estando en el dominio de la frecuencia, se busca la solución
para el sistema de ecuaciones definidos a partir de la definición del Tensor de impedancias
56
(Z) (Ecuación 1), el cual relaciona las componentes horizontales del campo eléctrico (𝐸𝑥 , 𝐸𝑦 )
y magnético (𝐻𝑥 , 𝐻𝑦 ).
𝑍𝑥𝑦 (𝜔)
(𝐸𝐸𝑥 (𝜔)
(𝜔)
) = ( 𝑍𝑥𝑥 (𝜔)
𝑍
) 𝐻𝑥 (𝜔)
(𝐻 (𝜔)
) Ecuación 1
𝑦 𝑦𝑥 (𝜔) 𝑍𝑦𝑦 (𝜔) 𝑦
Por otra parte, se realiza el cálculo de los vectores de inducción (𝑇𝑥 , 𝑇𝑦 ) (Ecuación 2), que
corresponde a la función de transferencia que relaciona las componentes horizontales del
campo magnético con su componente vertical (𝐻𝑧 ):
𝐻𝑥 (𝜔)
𝐻𝑧 (𝜔) = (𝑇𝑥 (𝜔) 𝑇𝑦 (𝜔)) ( ) Ecuación 2
𝐻𝑦 (𝜔)
Para más información sobre la formulación teórica de método referirse a Chave y Jones
(2012).
La estimación del tensor de impedancias y los vectores de inducción se lleva a cabo a partir
de dos algoritmos independientes, la primera mediante el software ProcMT (Metronix,
2002) que utiliza una regresión por mínimos cuadrado, y la segunda con la librería SigMT
(Ajithabh y Patro, 2023) que implementa una estimación robusta para estimar las funciones
de transferencia.
❖ Obtención de curvas de resistividad aparente y fase
A partir del tensor de impedancias es posible calcular las curvas de resistividad aparente y
fase en función de la frecuencia. En la Figura 3.20 se muestran a modo de ejemplo las
curvas de resistividad aparente y fase de la estación TRAP5.
3.2.3. Modelamiento
Para la inversión de los datos se utiliza el software Winglink, el cual permite realizar el
control de calidad de los datos y posteriormente efectuar el modelamiento a través de
inversiones bidimensionales. Para esto, al software son importados tanto el tensor de
impedancia, como los vectores de inducción, pudiéndose visualizar las curvas de
resistividad aparente y fase.
57
Winglink posibilita la realización de un proceso de control de calidad de las funciones de

transferencia previo a la ejecución del modelamiento. Para esto, es posible eliminar aquellos
puntos que sean definidos como fuera de tendencia. Además, al disponer de información
de transiente electromagnético, es posible identificar y remediar la presencia de distorsión
galvánica, como el static shift.
Figura 3.20: Curvas de resistividad aparente y fase de la estación TRAP5. Eje x de la figura
corresponde al periodo en unidades de segundos.
58
❖ Inversión 2D
El algoritmo utilizado por Winglink para la inversión bidimensional es el de gradientes

conjugados no lineales (Rodi & Mackie, 2001). Este algoritmo calcula soluciones
regularizadas para el problema inverso de datos magnetotelúricos 2D (Figura 3.21). El
proceso de inversión se lleva a cabo mediante procesos iterativos a partir de una grilla
compuesta por celdas con su resistividad eléctrica determinada, a modo de obtener el
menor desajuste entre los datos y la respuesta sintética del modelo (RMS).
Para llevar a cabo las inversiones se consideraron los datos originales de cada estación,
utilizando ambas componentes de resistividad aparente y fase obtenidas (XY e YX). Se
definió un piso de error de 10% para la resistividad aparente y de 5% para la fase. La grilla
se creó en función del rango de frecuencias obtenidas para cada sitio.
❖ Análisis de dimensionalidad
A partir de las funciones de transferencia (impedancia y vectores de inducción) es posible

calcular distintos parámetros que permiten evaluar la dimensionalidad de la estructura de
resistividad eléctrica, que es relevante para valorar la validez de un modelo bidimensional.
En este estudio se utilizan las siguientes metodologías:
o Elipses del tensor de fase
El tensor de fase, definido por Caldwell et al. (2004), es un tensor de segundo orden que
no se encuentra afectado por la distorsión galvánica del campo eléctrico, y que se
representa gráficamente mediante elipses, denominadas como elipses de fase. La forma y
orientación de éstas entregan información sobre la dimensionalidad del medio. En los casos
unidimensionales, donde la resistividad varía sólo en profundidad, las elipses tienden a ser
circulares. Por otro lado, en situaciones 2D donde la resistividad varía en profundidad y en
una dirección preferencial, las elipses se deforman y se alinean en la dirección del strike
geoeléctrico. En situaciones tridimensionales suelen tener un comportamiento más errático.
59
Figura 3.21: Visualización de las curvas de resistividad aparente (arriba) y de la fase (abajo). Eje x de las figuras corresponde al periodo en
unidades de segundos Derecha: Visualización modelo 2D en Winglink. Paleta de colores a la derecha indica valores de resistividad.
60
Un parámetro que se desprende del cálculo del tensor de fase es el ángulo de distorsión β,
el que es utilizado como un índice de la dimensionalidad, y que representa el grado de
simetría del tensor de fase. Este invariante rotacional tiende a cero en situaciones 2D y
aumenta su valor en situaciones 3D. En general un valor mayor a 3° es indicativo de una
situación tridimensional.
o Strike geo-eléctrico
El strike geo-eléctrico es un parámetro dimensional que permite determinar la existencia

de un rumbo eléctrico predominante que controla la distribución de resistividad eléctrica
del medio, y por lo tanto la circulación de corrientes en subsuperficie. En caso de existir
una estructura de resistividad bidimensional es posible rotar los ejes de coordenadas del
tensor de impedancia para la correcta determinación de los modos TE y TM.
El cálculo del strike se realiza mediante la librería implementada en Python MTpy (Krieger
y Peacock, 2014), que permite la estimación de un rumbo eléctrico a través de tres
metodologías diferentes: mediante el algoritmo de Zhang et al. (1987), a partir de la
definición del tensor de fase y, finalmente, a partir de la orientación de los vectores de
inducción.
61
3.3. Transiente electromagnético
El método transiente electromagnético (TEM) o método electromagnético en el dominio del

tiempo (TDEM), es uno de los métodos más utilizados en Geofísica Aplicada. Corresponde
a una técnica de exploración geofísica perteneciente al grupo de los métodos
electromagnéticos (EM), que son utilizados para mapear las estructuras de resistividad
subsuperficial basándose en el fenómeno de inducción, que permite relacionar la variación
de campos electromagnéticos en el tiempo con las propiedades eléctricas de la Tierra. Los
métodos EM incluyen una gran variedad de técnicas, aplicaciones, alcances e
interpretaciones, sin embargo, los principios son los mismos, y cada técnica involucra la
medición de una o más componentes del campo eléctrico o magnético, ya sea de una fuente
natural o artificial. En particular, el método TEM corresponde a un método activo, ya que
utiliza como campo primario inductor una fuente artificial, a diferencia de otros, como el
método magnetotelúrico (MT) que utiliza fuentes naturales.
Existen diversas configuraciones de medición en el método transiente electromagnético,

pero la más utilizada, y la considerada en este estudio, consiste en la utilización de un loop
cuadrado horizontal (loop transmisor o bobina transmisora) sobre el suelo, por el que se
hace circular una corriente directa (turn-on) que induce un campo magnético primario. En
un tiempo t=0 esta corriente es interrumpida de manera abrupta (turn-off), lo que produce
una variación del campo magnético primario, que es capaz de inducir corrientes (corrientes
de Eddie) en los cuerpos conductores existentes en el subsuelo. Luego, las corrientes que
han sido inducidas en los conductores generan un campo magnético secundario que
decaerá en el tiempo. Este campo magnético secundario es medido como un decaimiento
de voltaje por un segundo loop (loop receptor o bobina receptora), en ausencia del campo
magnético primario (Figura 3.22). Este decaimiento es el que contiene información sobre
la resistividad del subsuelo, ya que la magnitud y distribución de las corrientes inducidas
dependen de las propiedades conductoras del medio y determinará cómo es el proceso de
difusión de los campos electromagnéticos.
62
Figura 3.22: Figura esquemática de cómo se realiza la adquisición de datos de una estación TEM y
sus fundamentos físicos de interacción con el subsuelo. Obtenida de Dentith & Mudge, 2014.
El flujo de adquisición, procesamiento y modelamiento de datos TEM consiste en distintas

etapas que son detalladas a continuación (Figura 3.23).
Figura 3.23: Flujo de adquisición, procesamiento y modelamiento de datos TEM.
63
3.3.1. Terreno
❖ Adquisición
En base a los requerimientos y objetivos del estudio, se realizó la exploración geofísica

mediante el método transiente electromagnético realizando un total de 100 estaciones de
medición (Figura 3.24), distribuidas en distintos sectores de la zona de estudio, con énfasis
en sectores cercanos a la Cordillera de la Costa y a la Cordillera Principal, donde la
posibilidad observar presencia de basamento mediante esta metodología es mayor, y
considerando también las zonas con altos gravimétricos donde es importante tener mayor
información para el modelamiento de cuenca sedimentaria.
El equipo utilizado para la adquisición de datos de este estudio corresponde al WalkTEM2

RX Estándar, de la empresa ABEM (para más información referirse a ABEM (2021)). La
configuración de medición de cada una de las estaciones fue tipo INLOOP, utilizando el
transmisor Tx8 de 40 x 40 m y la antena receptora RC-5 de 0,5 x 0,5 m. El script de
medición escogido, que determina la configuración de la adquisición de datos fue
10ms_40gates, que mide 2 momentos (o inyecciones de corriente): Low Moment (LM) que
inyecta una corriente de 1 A y High Moment (HM) que inyecta una corriente de 7 A. Por
cada ciclo de medición realiza un total de 5.340 mediciones para LM, 2.500 mediciones
para HM, además de 500 mediciones de ruido ambiental. Este script mide durante una
ventana de tiempo total de 10 ms posterior al apagado de la inyección de corriente (turn-
off) y las mediciones del campo magnético secundario se realizan en un total de 40
ventanas de tiempo o gates (Figura 3.25). Para más información sobre la formulación
teórica del método referirse a Ward y Hohmann (1987).
64
Figura 3.24: Mapa de distribución de estaciones TEM.
65
Figura 3.25: Esquema de medición. (a) Pulso transmisor. (b) Fuerza electromagnética inducida. (c)
Campo magnético secundario medido en distintas gates o ventanas de tiempo. Obtenida de ABEM
(2021).
❖ Control de calidad en terreno
Durante el terreno es posible hacer un control de calidad de datos visualizando la tendencia

y forma de las curvas de caída de voltajes para LM, HM y ruido, durante y después de la
medición. Esto brinda un primer acercamiento a la calidad de los datos según la suavidad,
dispersión de las curvas, y polaridad, y permite repetir las mediciones en caso de ser
necesario.
Para realizar un correcto procesamiento de los datos es necesario hacer una inspección
inicial de ellos, para de esta manera evaluar su calidad, lo que permite posteriormente
hacer una limpieza acertada de estos, que luego es considerada al momento de efectuar la
interpretación de los modelos. Este análisis considera la dispersión de los datos (en función
del stacking realizado), el análisis del nivel de ruido, análisis de relación de High/Low
Moment, obtención de una curva única de decaimiento de voltaje, y cálculo de la curva de
resistividad aparente.
66
❖ Stacking de mediciones
El stacking de las mediciones consiste en el apilamiento o unión de todas las mediciones

realizadas en una sola. Esto se realiza mediante el cálculo de un promedio de las más de
2.000 repeticiones de las mediciones para las distintas inyecciones de corriente (LM y HM),
y con esto se obtiene una curva única de decaimiento de voltaje para cada una de ellas. El
análisis del stacking permite tener una idea de cómo es la dispersión y calidad de los datos
registrados para posteriormente hacer una limpieza de datos e interpretación coherente. A
continuación, se presenta a modo de ejemplo, el stacking de la estación ITA44 (Figura
3.26).
Figura 3.26: Ejemplo de stacking de estación ITA44. En color rojo, azul y gris los datos
correspondientes a LM, HM y ruido respectivamente. En líneas con trasparencia se muestran todas
las mediciones realizadas, y en línea sólida los valores apilados (stacking) de cada momento.
67
❖ Análisis de nivel de ruido
El análisis del nivel de ruido consiste en comparar las curvas de decaimiento de voltaje con
las curvas de ruido electromagnético ambiental de la estación. Esto nos permite
comprender y evaluar hasta qué tiempos la curva de decaimiento de voltaje se puede
considerar válida, entendiendo que en los tiempos en que ambas curvas se traslapan
posiblemente el voltaje observado estará siendo afectado por el ruido natural del lugar. En
la Figura 3.26 se puede observar la comparación de la curva de decaimiento (rojo y azul)
de voltaje versus la curva de ruido ambiental (gris) de la estación ITA44, de la cual se
puede deducir que los datos de mejor calidad para utilizar en el modelamiento en este caso,
por encontrarse sobre el nivel de ruido ambiental, son hasta 100 ms aproximadamente.
❖ Cálculo de resistividad aparente
A partir de las curvas de decaimiento de voltaje, que es el dato directo de medición del
equipo, es posible calcular lo que se conoce como resistividad aparente para ambas
inyecciones de corriente (Figura 3.27). Si bien, también es posible interpretar las curvas
de decaimiento, este parámetro nos entrega de forma más directa información respecto a
cómo es el cambio o contraste de resistividad eléctrica en profundidad.
❖ Relación High Moment/Low Moment
En general, las mediciones de Low Moment deberían ser más confiables en los tiempos
tempranos, mientras que las de High Moment lo serían para tiempos tardíos. En condiciones
ideales, las curvas de ambas mediciones se deberían traslapar de buena manera para
tiempos intermedios. El análisis de cómo se traslapan ambas curvas nos permite escoger
qué datos utilizar para cada tiempo de medición y obtener una única curva final de
decaimiento de voltaje y resistividad aparente para utilizar en la inversión (Figura 3.28).
68
Figura 3.27: Ejemplo de cálculo de resistividad aparente para estación ITA44.
69
Figura 3.28: Ejemplo de la curva final de decaimiento de voltaje de la estación ITA44, obtenida a partir del traslape de las curvas de decaimiento de LM y HM.
70
3.3.3. Modelamiento
En base a todos los análisis anteriores, antes de realizar las inversiones TEM 1D, se realiza
una limpieza de las curvas de decaimiento de voltaje y resistividad aparente, eliminando
todos aquellos puntos que se consideran anómalos, fuera de rango o sobre el nivel de ruido
ambiental. En la Figura 3.29, se muestra un ejemplo de la limpieza de la estación ITA44,
donde se eliminan todos los puntos que se encuentran en periodos de tiempo donde la
curva de decaimiento de voltaje se traslapa con el nivel de ruido ambiental y los valores
tienen un comportamiento dispersivo (Figura 3.28).
Figura 3.29: Ejemplo de limpieza de datos de curva de decaimiento de voltaje de la estación ITA44.
En gris los puntos que se han eliminado para la inversión.
71
❖ Inversión 1D
El software utilizado para hacer inversiones 1D de TEM corresponde a Winglink, que permite
realizar 2 tipos de modelos, un modelo de capas simple (layered model) y un modelo
suavizado de varias capas (smooth model). En este estudio se decidió trabajar con modelos
de capas simples, que consiste en un modelo capaz de interpretar los datos con la menor
cantidad de capas posibles (Figura 3.30). De esta manera se evita una sobre interpretación
de los datos.
Figura 3.30: Modelo de capas y ajuste respecto a la curva de decaimiento de voltaje de la estación de
TEM ITA44.
72
4. RESULTADOS E INTERPRETACIÓN
4.1. Magnetotelúrica
4.1.1. Funciones de transferencia
Al evaluar los resultados obtenidos para el tensor de impedancia y vectores de inducción,

se puede observar una calidad muy variable en los datos recopilados a lo largo de la zona
de estudio, producido principalmente por la distribución de fuentes de ruido
electromagnético (Mapa de ruido, Anexo B), siendo las más relevantes las líneas de alta
tensión, antenas y rutas con un gran tráfico automovilístico, como lo es la Ruta 5,
llegándose a obtener tanto estaciones con una alta calidad como las ubicadas en el sector
de Trapiche – TRAP1 y otras como YUN 4 (Yungay), en la que se observa una gran
dispersión en la información obtenida. Lo anterior es fácilmente apreciable al observar la
Figura 4.1, en donde existe una clara diferencia en la calidad de la información recopilada
por estas estaciones. Por otro lado, cabe destacar que de todas las componentes medidas,
la que fue principalmente afectada en los sectores con un mayor nivel de ruido, fue la
componente Hz dificultando la interpretación de los vectores de inducción.
Por otro lado, los datos recopilados en la zona de estudio abarcaron un rango promedio de
frecuencias entre 10.000 Hz hasta 1 Hz, en el cual es posible apreciar que las estaciones
desplegadas hacia los sectores cordilleranos, presentan en frecuencias altas (> 10.000 Hz)
valores de resistividad aparente que oscilan entre los 30 Ωm – 300 Ωm, llegando hasta los
3.000 Ωm en frecuencias más bajas (10 Hz – 1 Hz) sugiriendo la presencia de estructuras
resistivas, mientras que las estaciones ubicadas en el valle registran una tendencia a
valores mayores de resistividad aparente 30 Ωm – 200 Ωm en frecuencias altas,
disminuyendo hacia frecuencias bajas, lo que sugiere la presencia de zonas conductivas en
profundidad.
73
Figura 4.1: Curvas de resistividad y fase de las estaciones TRAP1 (izquierda) y YUN4 (derecha).
74
4.1.2. Análisis de dimensionalidad
❖ Secciones de tensor de fase y skew:
El análisis del tensor de fase de las estaciones desplegadas en la zona de estudio permite
evaluar la dimensionalidad de las estructuras presentes en el subsuelo, mostrando
diferencias claras entre los sectores estudiados, teniendo que en algunos predomina la
unidimensionalidad y/o bidimensionalidad del medio, sin embargo, en otros se ve una clara
incidencia de estructuras tridimensionales.
Figura 4.2: Sección de tensor de fase y parámetro skew de uno de los perfiles del sector de Maravilla
(MAR2 y MAR3).
75
En el caso de la Figura 4.2 se logran apreciar elipses bastante circulares hasta los 0.01 s
(100 Hz) en el caso de la estación MAR3 y hasta los 0,001 s (1.000 Hz) en el caso de la
estación MAR2, con valores de skew (β) menores a 3°, lo que podría indicar la presencia
de una estructura de resistividad bidimensional al menos hasta los 50 Hz. Por otro lado, las
frecuencias más bajas muestran elipses mucho más achatadas y con skew altos (β > 3°)
que podrían ser indicadores de la presencia de una estructura de resistividad tridimensional
en profundidad.
Para mayor información consultar los anexos de la sección MT, en donde se presenta el
resto de los análisis de tensor de fase y skew (Anexo B).
❖ Strike geo-eléctrico
Al igual que en el caso anterior, el análisis del strike geo-eléctrico entrega diferentes
configuraciones a lo largo de toda la zona de estudio, teniendo casos en los que el strike
geo-eléctrico es prácticamente nulo a otros en que el ángulo es de más de 45°. Es
importante mencionar que la orientación del rumbo geo-eléctrico posee una ambigüedad
intrínseca de 90°, por ende, durante el análisis de dimensionalidad se debe discernir cuál
de los ángulos determinados mediante diferentes algoritmos se debe seleccionar para
posteriormente ser utilizado en la rotación de los modelos bidimensionales, de acuerdo a
la coherencia que este ángulo posea con un determinado perfil de estudio.
En el siguiente caso se presenta el strike geo-eléctrico determinado para el perfil MAR2-3

ubicado en el sector de Maravilla (Figura 4.3), el cual posee una orientación N55°E.
Figura 4.3: Cálculo de strike para el perfil MAR2-3 ubicado en el sector de Maravilla.
Para más información consultar los anexos de la sección MT, donde se presentan los
análisis de strike geo-eléctrico restantes (Anexo B).
76
4.1.3. Modelos de resistividad eléctrica
Los modelos de resistividad eléctrica de las estaciones MT realizados se han separado en 3

grandes sectores para su descripción: zona norte, centro y sur, tal y como se muestra en
la (Figura 3.18).
Para una mejor visualización y comprensión de los modelos obtenidos, se presentan los
modelos de resistividad v/s profundidad con la mejor calidad de datos obtenidos en cada
zona junto a su interpretación hidrogeológica respectiva. La totalidad de los perfiles
generados se presenta en el apartado de ANEXO B.
❖ Sector Norte
o Trapiche (TRAP 1-7)
En el sector de Trapiche se realizaron un total de 7 mediciones MT cuya información

geográfica se presenta en la Tabla 4.1 y en la Figura 4.4.
Tabla 4.1: Ubicación de estaciones geofísicas MT realizadas en el sector de Trapiche (WGS84 - 18S).
Localidad Punto MT Este UTM (m) Norte UTM (m)

Trapiche TRAP1 761.547 5.977.160
Trapiche TRAP2 761.400 5.976.895
Trapiche TRAP3 761.277 5.976.620
Trapiche TRAP4 761.157 5.976.350
Trapiche TRAP5 761.659 5.977.443
Trapiche TRAP6 760.936 5.976.249
Trapiche TRAP7 761.668 5.975.950
77
Figura 4.4: Ubicación de las estaciones de magnetotelúrica en el sector de Trapiche.
El perfil representado en la Figura 4.5, con orientación NE-SO, indica valores de resistividad
menores a 50 Ωm (tonalidades lilas) los que podrían asociarse a la Fm. Mininco que se
correlacionan de buena forma con los materiales finos como sedimentitas lagunares y
fluviales, y limolitas que componen a esta formación. Se identifica que esta unidad presenta
un mayor espesor hacia la porción norte del perfil donde alcanza profundidades cercanas a
250 m b.n.t. (TRAP5), mientras que, hacia el sur en TRAP6, alcanza aproximadamente los
100 m b.n.t. Bajo esta unidad se aprecia valores de resistividad medios-altos (tonalidades
azules a verdes), que indican el paso a materiales más gruesos y/o con mayor grado de
compactación.
78
Entre las mediciones TRAP5-TRAP1-TRAP2, se aprecia un aumento progresivo en los valores
de resistividad mayores a 500 Ωm (tonalidades amarillas a rojas) a una profundidad de 400
m b.n.t. lo cual indicaría la aparición del basamento cristalino, que hacia su techo podría
tener un grado variable de fracturamiento, meteorización y/o alteración, el cual disminuiría
progresivamente en profundidad. Por otro lado, hacia el sur en las mediciones TRAP2-
TRAP3-TRAP4-TRAP6 se identifica un mayor espesor de la capa de resistividad media-alta
(tonalidad verdosa), con la aparición del basamento alterado y/o fracturado bajo los 1.000
m b.n.t. en TRAP6.
79
Figura 4.5:Modelo de resistividad eléctrica y sección interpretada de los puntos MT TRAP1 a TRAP6.
80
o Maravilla (MAR1-3)
En el sector de Maravilla se realizaron un total de 3 mediciones MT cuya información

Tabla 4.2: Ubicación de estaciones geofísicas MT realizadas en el sector de Maravilla (WGS84 – 18S).

Maravilla MAR1 755.483 5.969.362
Maravilla MAR2 755.328 5.969.340
Maravilla MAR3 755.383 5.968.344
Figura 4.6: Ubicación de las estaciones de magnetotelúrica en el sector de Maravilla.
81
El perfil representado en la Figura 4.7, con orientación N-S, al igual que el del sector
Trapiche, indica una primera capa conductiva con valores de resistividad menores a 50 Ωm
(tonalidades lilas) posiblemente asociados a la Fm. Mininco. El espesor de esta unidad se
aprecia relativamente homogéneo con una potencia de 100 m en promedio y con una
creciente extensión hacia el límite sur del perfil (MAR3) donde alcanza profundidades
cercanas a 200 m b.n.t. Bajo esta unidad se aprecia un progresivo aumento en los valores
de resistividad, donde las zonas con tonalidades azules a verdes indicarían el paso a
materiales más gruesos y/o con mayor grado de compactación.
Los valores de resistividad mayores a 500 Ωm (tonalidades amarillas a rojas) se aprecian

a una profundidad cercana a 200 m b.n.t. en el extremo norte del perfil (MAR2), y hacia el
sur se profundiza y empieza a distinguirse bajo los 350 m b.n.t. en MAR3. Esto indicaría la
aparición del basamento cristalino, que hacia su techo podría tener un grado variable de
fracturamiento, meteorización y/o alteración, el cual disminuiría progresivamente en
profundidad. El motivo de su aparición a una menor profundidad que en el sector de
Trapiche se correlaciona con una mayor cercanía a la unidad geomorfológica de la Cordillera
de la Costa.
82
Figura 4.7: Modelo de resistividad eléctrica y sección interpretada de los puntos MT MAR2 – MAR3.
83
o Cato (CATO1-6)
En el sector de Cato se realizaron un total de 6 mediciones MT cuya información geográfica

se presenta en la Tabla 4.3 y en la Figura 4.8.
Tabla 4.3: Ubicación de estaciones geofísicas MT realizadas en el sector de Cato (WGS84 – 18S).

Cato CATO1 782.318 5.950.842
Cato CATO2 782.318 5.949.292
Cato CATO3 783.284 5.949.760
Cato CATO4 783.559 5.950.633
Cato CATO5 785.302 5.954.806
Cato CATO6 785.171 5.954.948
Figura 4.8: Ubicación de las estaciones de magnetotelúrica en el sector de Cato.
84
Las mediciones CATO5 – CATO6 se realizaron cercanas al Río Ñuble con una orientación
NO-SE permitieron generar el perfil de la Figura 4.9, el cual presenta una primera capa
superficial de unos 50 m de espesor con resistividades medias-bajas (tonalidades celestes-
verdosas) asociados a los depósitos fluvio-aluviales que se encuentran rellenando los
cauces de los ríos y que presentarían en estos primeros metros una composición
mayoritariamente gruesa como gravas y ripios. Bajo esta capa y hasta los 100 m de
profundidad, se aprecia un descenso en los valores de resistividades (tonalidades lilas), lo
cual se asocia a la aparición de la Formación Mininco (o posiblemente Fm. La Montaña), la
cual presentaría un predominio de material fino como arcillas y/o limos en este tramo.
Subyaciendo esta capa conductiva, se aprecia un aumento en la resistividad hasta los 200
– 250 m de profundidad que podría corresponder a un miembro distinto de la Formación
Mininco con mayor grado de compactación y contenido de materiales gruesos como bolones
y/o gravas. Desde los 200 – 250 m hasta al menos los 800 m de profundidad las
resistividades vuelven a decrecer hasta valores menores a 30 Ωm, lo cual podría
corresponder a un Miembro inferior de la Formación Mininco con alto contenido arcillo-
limoso o bien a una nueva unidad altamente conductora que no se ha reconocido
anteriormente.
85
Figura 4.9: Modelo de resistividad eléctrica y sección interpretada de los puntos MT CATO5 – CATO6.
86
o Cachapoal (CACHA1-2)
En el sector de Cachapoal se realizaron un total de 4 mediciones MT cuya información

geográfica se presenta en la Tabla 4.4 y Figura 4.10.
Tabla 4.4: Ubicación de estaciones geofísicas MT realizadas en el sector de Cachapoal (WGS84- 18S).

Cachapoal CACHA1 791.676 5.958.305
Figura 4.10: Ubicación de las estaciones de magnetotelúrica en el sector de Cachapoal.
87
El perfil O-E generado a partir de las mediciones CACHA1-CACHA2 en la Figura 4.11, se
ubica hacia el oeste de las mediciones de Cato y cercano al río Ñuble. De manera similar al
perfil descrito anteriormente, se aprecia una primera capa superficial de resistividad media
y con un espesor cercano a los 50 m correlacionado con los depósitos fluvio-aluviales que
rellenan el cauce del río. Bajo esta capa, y de manera similar al perfil de Cato, aparecería
la Formación Mininco con un estrato superior de baja resistividad que alcanza los 100 m de
profundidad, seguido de un aumento en la resistividad asociado a un tramo con mayores
contenidos de materiales granulométricamente más gruesos que alcanzarían profundidades
cercanas a 500 m.
La capa subyacente más profunda podría corresponder a un Miembro inferior de la

Formación Mininco con alto contenido arcillo-limoso o bien a una nueva unidad altamente
conductora que no se ha reconocido anteriormente. Al comparar las profundidades de los
cambios en los valores de resistividad obtenidos en este perfil con los de la localidad de
Cato, se aprecia una profundización de esta última capa a medida que nos desplazamos
hacia el este.
88
Figura 4.11: Modelo de resistividad eléctrica y sección interpretada de los puntos MT CACHA1 –
CACHA2.
89
o Camino Cachapoal (CCACHA1-12)
En el sector de Camino Cachapoal se realizaron un total de 12 mediciones MT cuya

información geográfica se presenta en la Tabla 4.5 y en la Figura 4.12.
Tabla 4.5 Ubicación de estaciones geofísicas MT realizadas en el sector de Camino Cachapoal

(WGS84 – 18S).

Camino Cachapoal CCACHA1 775.308 5.963.595
Las mediciones realizadas en Camino Cachapoal se encuentran al noroeste de la localidad

de Cachapoal. El perfil representado en la Figura 4.13, indica una primera unidad de
resistividad media-alta (tonalidades celestes-verdosas) en sus primeros 60-70 m, asociado
a material disgregado que rellena el valle con altos contenidos de material grueso como
bolones, ripio y gravas. Subyaciendo a esta primera capa, se identifica una progresiva
disminución en los valores de la resistividad (tonalidades azules a moradas), indicando al
igual que los perfiles de Cato y de Cachapoal, la aparición de la Formación Mininco y de su
estrato compuesto mayoritariamente por material fino como arcillas y/o limos.
90
Figura 4.12: Ubicación de las estaciones de magnetotelúrica en el sector Camino Cachapoal.
91
Figura 4.13: Modelo de resistividad eléctrica y sección interpretada de los puntos MT CCACHA12 –
CCACHA11.
92
o Ninquihue (QUI1-4)
En el sector de Ninquihue se realizaron un total de 4 mediciones MT cuya información

Tabla 4.6: Ubicación de estaciones geofísicas MT realizadas en el sector de Ninquihue (WGS84 –

18S).

Ninquihue QUI1 768.056 5.958.622
Ninquihue QUI2 768.135 5.958.639
Ninquihue QUI3 768.850 5.958.609
Ninquihue QUI4 768.770 5.958.612
Figura 4.14: Ubicación de las estaciones de magnetotelúrica en el sector de Ninquihue.
93
El perfil generado a partir de las mediciones QUI3-QUI4 en la Figura 4.15, se ubica al
suroeste de la ciudad de San Carlos, en la unidad geomorfológica del Valle Central. Los
resultados revelan una potente secuencia conductora con resistividades menores a 50 Ωm
(tonalidades moradas) hasta los 110-130 m b.n.t., asociadas a depósitos que se encuentran
rellenando el valle y que presentarían altos porcentajes de material fino como arenas,
arcillas y/o limos. Bajo esta unidad las resistividades decrecen progresivamente
(tonalidades lilas), lo cual, podría corresponder a un estrato de la Formación Mininco con
mayores contenidos de arcillas y/o limos en variables condiciones de humedad y
compactación.
94
Figura 4.15:Modelo de resistividad eléctrica y sección interpretada de los puntos MT QUI3 – QUI4.
95
o El Ala (ALA 1-4)
En el sector de El Ala se encuentra al suroeste de las mediciones anteriores en las cercanías

a la Cordillera de la Costa, en esta localidad se realizaron un total de 4 mediciones MT cuya
información geográfica se presenta en la Tabla 4.7 y la Figura 4.16.
Tabla 4.7: Ubicación de estaciones geofísicas MT realizadas en el sector de El Ala (WGS84 – 18S).

El Ala ALA1 747.001 5.948.238
El Ala ALA2 747.062 5.948.286
El Ala ALA3 749.266 5.948.366
El Ala ALA4 746.457 5.948.963
Figura 4.16: Ubicación de las estaciones de magnetotelúrica en el sector El Ala.
96
Con la medición ALA3 se generó un perfil de resistividad vs profundidad representado en la
Figura 4.17, el cual, presenta una primera capa superficial de 60 m de espesor con
resistividades medias-bajas (tonalidades celestes-verdosas) asociados a los depósitos
fluvio-aluviales que se encuentran rellenando el cauce del río Ñuble y que presentarían una
composición mayoritariamente gruesa con gravas y ripios. Bajo los 60 m y de manera
progresiva se identifica un aumento en las resistividades (tonalidades amarillas-naranjas),
que estarían indicando la aparición de un basamento meteorizado, alterado y fracturado
hasta profundidades cercanas a los 100 m, y a profundidades mayores a 190 m se esperaría
la aparición del basamento cristalino inalterado y fresco. Cabe mencionar que esta
medición, al estar en el lecho del río, los valores de resistividad de la roca fracturada pueden
ser más bajos debido a la infiltración de agua en esta zona.
Figura 4.17: Modelo de resistividad eléctrica y sección interpretada del punto MT ALA3.
97
o Coihueco (COI 1-4)
El sector de Coihueco se emplaza en el extremo este de la cuenca, cercana a la Cordillera

Principal, donde se realizaron un total de 4 mediciones MT cuya información geográfica se
presenta en la Tabla 4.8 y en la Figura 4.18.
Tabla 4.8: Ubicación de estaciones geofísicas MT realizadas en el sector de Coihueco (WGS84 – 18S).

Coihueco COI1 784.385 5.936.968
Coihueco COI2 784.309 5.936.975
Coihueco COI3 784.221 5.936.504
Coihueco COI4 784.150 5.936.510
Figura 4.18: Ubicación de las estaciones de magnetotelúrica en el sector de Coihueco.
98
Las mediciones COI4 – COI3 permitieron genera un perfil con orientación O-E (Figura 4.19),
el cual presenta una primera capa superficial de unos 60 a 70 m de espesor con
resistividades medias-altas (tonalidades verdosas) asociados a los depósitos de la
Formación ¿Montaña/Malleco? que se compone por una mezcla de conglomerados,
conglomerados arenosos, areniscas, intercalaciones de limos y depósitos piroclásticos. Bajo
esta capa y hasta al menos 300 m de profundidad se identifica una unidad de menor
resistividad (tonalidad azul), la cual podría corresponder a condiciones de humedad y/o
saturación con un aumento en la proporción de material fino como limos y/o arcillas lo cual
estaría generando este nivel de mayor conductividad en la zona.
99
Figura 4.19: Modelo de resistividad eléctrica y sección interpretada de los puntos MT COI4 – COI3.
100
❖ Sector Centro
o Viña Riberas del Chillán (VRC 1-2)
En el sector de Viña Riberas del Chillán se realizaron un total de 2 mediciones MT cuya

información geográfica se presenta en la Tabla 4.9 y Figura 4.20.
Tabla 4.9: Ubicación de estaciones geofísicas MT realizadas en el sector de Viña Ribera (WGS84 –
18S).

Viña Ribera VRC1 738.841 5.943.831
Viña Ribera VRC2 738.873 5.943.784
Figura 4.20: Ubicación de las estaciones de magnetotelúrica en el sector de Viña Ribera.
101
El perfil representado en la Figura 4.21, con orientación NO-SE, se encuentra aledaño a la
cordillera de la costa y al igual que los modelos previos descritos en esta zona, indica una
primera capa conductiva con valores de resistividad menores a 50 Ωm (tonalidades lilas)
que alcanza los 10 m de espesor y estaría asociado a material disgregado que rellena el
valle con altos contenidos de arcillas y/o limos. Bajo esta unidad se aprecia un progresivo
aumento en los valores de resistividad, donde, las tonalidades azules a verdes indicarían el
paso a materiales más gruesos y/o con mayor grado de compactación.

a una profundidad cercana a 30 m b.n.t. en el extremo noroeste del perfil (VRC1), mientras
que hacia el sureste se aprecia a 50 m b.n.t. aproximadamente en VRC2. Estas unidades
se asocian con la aparición del basamento con un grado de fracturamiento y/o alteración
hacia el techo, y una disminución de estas en profundidad (sobre los 100 m).
102
Figura 4.21: Modelo de resistividad eléctrica y sección interpretada de los puntos MT VRC1 – VRC2.
103
o Bellavista (RC 1-3)
En el sector de Bellavista se realizaron un total de 3 mediciones MT cuya información

Tabla 4.10: Ubicación de estaciones geofísicas MT realizadas en el sector de Bellavista (WGS84 –

18S).

Bellavista RC1 742.077 5.943.318
Bellavista RC2 742.076 5.943.253
Bellavista RC3 742.019 5.943.197
Figura 4.22: Ubicación de las mediciones de magnetotelúrica en el sector de Bellavista.
104
Hacia el este del perfil anterior se realizaron las mediciones MT RC1 – RC2 – RC3, con
orientación NE-SO, con las que se generó el perfil representado en la Figura 4.23. De
manera análoga al modelo anterior, se aprecia una primera capa conductiva con valores de
resistividad menores a 50 Ωm (tonalidades lilas) con un espesor mayor, que alcanza los 60
m aproximadamente, y que estaría asociado a material disgregado que rellena el valle con
altos contenidos de arcillas y/o limos. Bajo esta unidad se aprecia un progresivo aumento
en los valores de resistividad, donde, las tonalidades azules a verdes indicarían el paso a
materiales más gruesos y/o con mayor grado de compactación.

a una profundidad cercana a 100 - 120 m b.n.t, lo que estaría revelando la aparición del
basamento en condiciones de fracturamiento y/o alteración en la parte superficial del
cuerpo ígneo y en estado fresco a profundidades mayores a 200 m. La mayor profundidad
de aparición del basamento en este sector tiene coherencia debido al mayor
distanciamiento de la Cordillera de la Costa en comparación al sector de Viña Riberas del
Chillán.
105
Figura 4.23: Modelo de resistividad eléctrica y sección interpretada de los puntos MT RC1 – RC2 –
RC3.
106
o Liucura (LIU 1-11)
El sector de Liucura se emplaza en la misma unidad geomorfológica que las mediciones

anteriores, pero en dirección al suroeste, donde se realizaron un total de 11 mediciones MT
cuya información geográfica se presenta en la Tabla 4.11 y en la Figura 4.24.
Tabla 4.11: Ubicación de estaciones geofísicas MT realizadas en el sector de Liucura (WGS84 – 18S).

Liucura LIU1 733.374 5.918.300
Liucura LIU2 733.369 5.918.371
Liucura LIU3 735.159 5.918.988
Liucura LIU4 735.087 5.919.005
Liucura LIU5 735.411 5.917.201
Liucura LIU6 735.400 5.917.143
Liucura LIU7 733.459 5.915.979
Liucura LIU8 733.422 5.915.972
Liucura LIU9 731.574 5.918.734
Liucura LIU10 731.505 5.918.705
Liucura LIU11 731.816 5.915.558
Figura 4.24: Ubicación de las estaciones de magnetotelúrica en el sector de Liucura.
107
El perfil representado en la Figura 4.25, se encuentra orientado SO-NE y cruza a los ríos
Itata y Diguillin en los tramos LIU11 – LIU8 y LIU8 – LIU6 respectivamente. De manera
general, el modelo presenta una primera capa de resistividades medias altas (tonalidades
verdes) seguidas de un progresivo aumento en las resistividades que indicarían la aparición
de roca en un estado variable de alteración, meteorización y/o fractuamiento.
El tramo entre LIU11 – LIU8 presenta un espesor de la primera capa en promedio de 200
m, asociado a material de los depósitos fluvio-aluviales que se encuentran rellenando los
cauces de los ríos que composicionalmente presentan en su mayoría granulometrías
gruesas como gravas y ripios. Bajo esta unidad el drástico incremento de los valores de
resistividad podría indicar la aparición del basamento en condiciones de alteración y
fracturamiento (tonalidades amarillas - naranjas) y se esperaría roca en estado fresco bajo
los 700 m aproximadamente.
Por otro lado, la sección LIU8 – LIU6 presenta un mayor espesor de la cubierta
sedimentaria, inclusive se detectan tramos de media-baja resistividad (tonalidades
celestes) que se pueden asociar a mayores contenidos de arenas, arcillas y/o limos en
comparación a la porción descrita anteriormente. El incremento de las resistividades se
aprecia cercano a los 500 m b.n.t. donde podría aparecer la roca en condiciones de
alteración y/o fracturamiento, que gradualmente cambiaría a roca inalterada a mayores
profundidades (tonalidades rojas).
108
Figura 4.25: Modelo de resistividad eléctrica y sección interpretada de los puntos MT LIU11 – LIU8 –
LIU6.
109
o Pinto (PIN 1-8)
En el sector de Pinto se realizaron un total de 8 mediciones MT cuya información geográfica

se presenta en la Tabla 4.12 y Figura 4.26.
Tabla 4.12: Ubicación de estaciones geofísicas MT realizadas en el sector de Pinto (WGS84 – 18S).

Pinto PIN1 776.504 5.925.073
Pinto PIN2 776.499 5.925.145
Pinto PIN3 775.589 5.928.829
Pinto PIN4 775.582 5.928.885
Pinto PIN5 779.213 5.925.119
Pinto PIN6 779.255 5.925.065
Pinto PIN7 778.531 5.927.940
Pinto PIN8 778.577 5.928.003
Figura 4.26: Ubicación de las estaciones de magnetotelúrica en el sector de Pinto.
110
El perfil representado en la Figura 4.27, se asemeja en gran medida a los perfiles de Cato
y Cachapoal, donde se tiene una primera capa superficial con espesor menor a 10 m y baja
resistividad, asociadas a material de relleno compuesto principalmente por arcillas y/o
limos, seguido de un aumento en los valores de resistividad (tonalidades verdes) hasta los
50 m b.n.t. aproximadamente, que correspondería a la Formación Montaña/Malleco?
mayoritariamente compuesta por bloques, gravas y arena. Bajo esta capa, aparecería la
Formación Mininco con un nivel superior de baja resistividad que alcanzaría los 120 m de
profundidad, seguido de un aumento en la resistividad asociado a un tramo con mayores
contenidos de materiales granulométricamente más gruesos hasta aproximadamente los
400 m.
La capa subyacente más profunda podría corresponder a un Miembro inferior de la

Formación Mininco con alto contenido arcillo-limoso o bien a una nueva unidad altamente
conductora que no se ha reconocido anteriormente.
111
Figura 4.27: Modelo de resistividad eléctrica y sección interpretada de los puntos MT PIN2 – PIN1.
112
o Cartago (CAR1)
En el sector de Cartago se realizó una medición MT cuya información geográfica se presenta

en la Tabla 4.13 y en la Figura 4.28.
Tabla 4.13: Ubicación de estaciones geofísicas MT realizadas en el sector de Cartago (WGS84 – 18S).

Cartago CART1 776.504 5.925.073
Figura 4.28: Ubicación de la medición de magnetotelúrica en el sector de Cartago.
El perfil representado en la Figura 4.29, identifica una potente capa de baja resistividad, la
cual probablemente presente altos contenidos de material fino como arcillas y/o limo que
113
se puede asociar a la Formación Montaña/Malleco? hasta los 600 m aproximadamente. Bajo
esta unidad se aprecia un cambio y aumento de los valores de resistividad (tonalidades
azules a verdes), lo cual es posible asociar con la aparición de la Formación Cura-Mallín
compuesta por una secuencia de rocas piroclásticas, sedimentitas clásticas, tobas y
brechas.
114
Figura 4.29: Modelo de resistividad eléctrica y sección interpretada de la estación de MT CAR1.
115
o Atacalco (ATA 1-8)
En el sector de Atacalco se realizaron un total de 8 mediciones MT cuya información

Tabla 4.14: Ubicación de estaciones geofísicas MT realizadas en el sector de Atacalco (WGS84 –

18S).

Atacalco ATA1 787.092 5.919.268
Atacalco ATA2 787.165 5.919.268
Atacalco ATA3 788.715 5.919.063
Atacalco ATA4 788.704 5.919.131
Atacalco ATA5 788.164 5.917.406
Atacalco ATA6 788.175 5.917.474
Atacalco ATA7 787.483 5.917.600
Atacalco ATA8 787.426 5.917.628
Figura 4.30: Ubicación de las estaciones de magnetotelúrica en el sector de Atacalco.
116
El perfil representado en la Figura 4.31, identifica una primera capa superficial con espesor
menor a 10 m y resistividad bajas, asociadas a material de relleno compuesto
principalmente por arcillas y/o limos, seguido de un aumento en los valores de resistividad
(tonalidades verdes) hasta los 40 m b.n.t. aproximadamente, que corresponderá a la
Formación Montaña/Malleco? mayoritariamente compuesta por bolones, gravas y arena.
Bajo esta capa, la baja de resistividad correspondería a la Formación Mininco de baja
resistividad (tonalidades moradas) que alcanzaría los 180 m de profundidad y que estaría
compuesta por altos contenidos de material fino como arcillas y/o limos.
La capa más profunda se caracteriza por mayores valores de resistividad eléctrica

(tonalidades celestes a verdes), lo cual se podría asociar con la presencia de la Formación
Cura-Mallín, compuesta por una secuencia de rocas piroclásticas, sedimentitas clásticas,
tobas y brechas.
117
Figura 4.31: Modelo de resistividad eléctrica y sección interpretada de los puntos MT ATA8 – ATA7.
118
o San Antonio (SAN 1-7)
En el sector de San Antonio se realizaron un total de 7 mediciones MT cuya información

Tabla 4.15: Ubicación de estaciones geofísicas MT realizadas en el sector de San Antonio (WGS84 –
18S).

San Antonio SAN1 767.653 5.918.962
San Antonio SAN2 767.620 5.918.887
San Antonio SAN3 770.044 5.916.304
San Antonio SAN4 770.093 5.916.291
San Antonio SAN5 766.036 5.916.023
San Antonio SAN6 769.143 5.918.666
San Antonio SAN7 769.072 5.918.688
Figura 4.32: Ubicación de las estaciones de magnetotelúrica en el sector de San Antonio.
119
El perfil representado en la Figura 4.33, presenta condiciones similares a los perfiles de
Pinto, donde se tiene una primera capa superficial con espesor menor a 10 m y bajas
resistividad bajas, asociadas a material de relleno compuesto principalmente por arcillas
y/o limos, seguido de un aumento en los valores de resistividad (tonalidades verdosas)
hasta aproximadamente los 50 m b.n.t., que corresponderá a la Formación
Montaña/Malleco? mayoritariamente compuesta por bolones, gravas y arena. Bajo esta
capa, aparecería la Formación Mininco con un estrato superior de baja resistividad que
alcanzaría los 150 m de profundidad, seguido de un leve aumento en la resistividad
asociado a un tramo con menor porcentaje de materiales finos hasta aproximadamente los
400 m.
La capa más profunda podría corresponder a un Miembro inferior de la Formación Mininco

con alto contenido arcillo-limoso o bien a una nueva unidad altamente conductora que no
se ha reconocido anteriormente.
120
Figura 4.33: Modelo de resistividad eléctrica y sección interpretada de los puntos MT SAN7 – SAN6.
121
o Trehualemu (TRE 1-7)
En el sector de Trehualemu se realizaron un total de 7 mediciones MT cuya información

Tabla 4.16: Ubicación de estaciones geofísicas MT realizadas en el sector de Trehualemu (WGS84 –

18S).

Trehualemu TRE1 777.079 5.908.677
Trehualemu TRE2 777.124 5.908.740
Trehualemu TRE3 777.939 5.910.754
Trehualemu TRE4 778.008 5.910.734
Trehualemu TRE5 777.775 5.910.451
Trehualemu TRE6 779.239 5.909.657
Trehualemu TRE7 779.181 5.909.696
Figura 4.34: Ubicación de las estaciones de magnetotelúrica en el sector de Trehualemu.
122
El perfil representado en la Figura 4.35, presenta una primera capa superficial con espesor
menor a 10 m y resistividad bajas, asociadas a material de relleno compuesto
(tonalidades verdosas) hasta aproximadamente los 80 m b.n.t., que corresponderá a la
Formación Montaña/Malleco? mayoritariamente compuesta por bolones, gravas y arena.
Bajo esta capa, se detectaría la aparición de la Formación Mininco asociado a la baja de
resistividad (tonalidades azules a moradas) que alcanzaría como mínimo los 200 m de
profundidad.
123
Figura 4.35:Modelo de resistividad eléctrica y sección interpretada de los puntos MT TRE3 – TRE4.
124
❖ Sector Sur
o Yungay (YUN 1-4)
En el sector de Yungay se realizaron un total de 4 mediciones MT cuya información

Tabla 4.17: Ubicación de estaciones geofísicas MT realizadas en el sector de Yungay (WGS84 – 18S).

Yungay YUN1 783.042 5.888.464
Yungay YUN2 783.085 5.888.500
Yungay YUN3 782.343 5.891.780
Yungay YUN4 782.328 5.891.709
Figura 4.36: Ubicación de las estaciones de magnetotelúrica en el sector de Yungay.
125
El perfil representado en la Figura 4.37, presenta una primera capa superficial con espesor
de 10 m y baja resistividad, asociadas a material de relleno compuesto principalmente por
arcillas y/o limos, seguido de un aumento en los valores de resistividad (tonalidades azules-
celestes) hasta aproximadamente los 30 m b.n.t., que podrían corresponder tanto a la
Formación Montaña/Malleco? como a la misma unidad de relleno, pero con un aumento en
la concentración de material grueso como gravas y arenas. Esta unidad podría encontrarse
en condiciones de saturación y/o humedad. Bajo esta capa, se detecta un aumento en los
valores de resistividad, que posiblemente indicarían un mayor grado de compactación de
los materiales o bien el paso a la Formación Cura-Mallín.
126
Figura 4.37: Modelo de resistividad eléctrica y sección interpretada de los puntos MT YUN1 – YUN2.
127
o Huepil (HUE 1-2)
En el sector de Huepil se realizaron dos mediciones MT cuya información geográfica se

presenta en la Tabla 4.18 y en la Figura 4.38.
Tabla 4.18: Ubicación de estaciones geofísicas MT realizadas en el sector de Huepil (WGS84 – 18S).

Huepil HUE1 784.035 5.877.364
Huepil HUE2 784.070 5.877.322
Figura 4.38: Ubicación de las estaciones de magnetotelúrica en el sector de Huepil.
128
El perfil representado en la Figura 4.39, identifica una primera capa superficial con un
espesor cercano a 10 m y bajas resistividades, asociadas a material de relleno compuesto
(tonalidades verdes) hasta los 40 m b.n.t. aproximadamente, que correspondería a la
Formación Montaña/Malleco? mayoritariamente compuesta por bloque, gravas y arena. En
la capa subyacente las resistividades aumentan (tonalidades amarillas), lo cual se asociaría
a materiales más compactos o a la aparición de la Formación Cura-Mallín compuesta por
una secuencia de rocas piroclásticas, sedimentitas clásticas, tobas y brechas.
129
Figura 4.39: Modelo de resistividad y secciones interpretadas de los puntos MT HUE1 – HUE2.
130
4.2. Transiente Electromagnético
4.2.1. Modelos de resistividad 1D
Los modelos de inversión 1D generados a partir de las estaciones TEM realizadas se

presentan en el ANEXO C. La Tabla 4.19, muestra la ubicación de los mejores modelos cuya
calidad de información y utilidad a la hora de complementar la información gravimétrica
permitieran interpretar de mejor manera la profundidad del basamento y el espesor
sedimentario mínimo de la cuenca.
Tabla 4.19: Ubicación de estaciones geofísicas TEM de mejor calidad para la interpretación geológica
y geofísica en la cuenca Itata (WGS84 – 18S).
Ubicación Punto TEM Este UTM (m) Norte UTM (m)

Extremo este (Cordillera Principal) TEM5 781.053 5.876.840
Oeste de Trapiche TEM9 745.997 5.974.396
Sur de Maravilla TEM17 755.222 5.964.965
Sur de El Ala TEM20 745.978 5.942.573
Entre Bulnes y Cabrero TEM26 742.825 5.928.323
Trapiche TEM45 761.688 5.977.460
Trapiche TEM50 761.373 5.976.778
El Ala TEM57 749.421 5.947.723
El Ala TEM60 746.669 5.947.491
Las mediciones realizadas en el extremo este, cercanas a la Cordillera Principal

corresponden a los TEM5-TEM7-TEM82-TEM83 y TEM85 (Figura 4.40), e indican de modo
general, estratos de resistividades medias que rondan los 100 Ωm y que corresponderían a
la Formación Montaña/Malleco?, algunos de estos modelos presentan intercalaciones
métricas a decimétricas, lo cual se podría asociar a diferencias porcentuales del material
fino como arcillas y/o limos que posean. Todas estas mediciones en profundidad detectan
131
en última instancia una bajada de resistividad que correspondería a la aparición de la
Formación Mininco con una composición predominante de materiales finos.
Las modelos de las estaciones de Trapiche (TEM45-TEM50) y al sur de Maravilla (TEM17),

representados en la Figura 4.41, presentan un predomino de resistividades bajas a lo largo
de todo el perfil que estaría asociado a la Formación Mininco con altos contenidos de
material fino.
Tal como se describió en los perfiles MT de estas zonas, el espesor de esta capa puede
llegar a alcanzar más de 400 m, por lo que la medición TEM9 se realizó al oeste de la
localidad de Trapiche cercana a la Cordillera de la Costa. Los resultados de esta estación
muestran un alza en los valores de resistividad superando los 1.000 Ωm a una profundidad
de 60 m b.n.t. lo que posiblemente este indicando la aparición del basamento en condiciones
de variable meteorización, alteración y/o fracturamiento.
La Figura 4.42 presenta los modelos de inversión 1D obtenidos de las mediciones TEM57,
TEM60 y TEM20 realizados en la localidad de El Ala y hacia el sur de esta. Los tres modelos
muestran valores de resistividad medios-bajos en sus primeros metros que podrían
asociarse tanto a depósitos fluvio-aluviales como a la Formación Mininco en las zonas de
menor resistividad. Subyaciendo a esta unidad se aprecia un incremento en los valores de
resistividad que se puede atribuir a la aparición del basamento en condiciones variables de
alteración. Esto se correlaciona de buena forma con las mediciones MT desarrolladas en el
mismo sector y que también dan cuenta de esta alza de resistividad. La profundidad a la
cual se detecta el aumento de la resistividad es variable, observándose a 30, 50 y 70 m
aproximadamente.
132
Figura 4.40: Modelos de inversión 1D de las mediciones que se encuentran en el extremo este
cercanos a la Cordillera Principal
133
Figura 4.41: Izquierda: Modelos de inversión 1D de las mediciones que se encuentran Sur de
Maravilla y en Trapiche. Derecha: Modelo de inversión 1D al Oeste de Trapiche.
134
Figura 4.42: Modelos de inversión 1D de las mediciones que se encuentran en El Ala y al Sur de El
Ala.
Por último, las mediciones realizadas entre las ciudades de Bulnes y Cabrero corresponden
a TEM26-TEM29-TEM30-TEM35 y TEM34, las cuales se emplazan en dirección norte-sur
(Figura 4.43). Los modelos de inversión 1D de esta zona muestran de manera general una
capa de resistividad media a baja en sus primeros metros, que se asociaría a material que
rellena los valles con variables contenidos de material grueso como gravas y arenas. Al ir
aumentando la profundidad, los modelos muestran un progresivo descenso en los valores
de resistividad lo cual correspondería a la aparición de la Formación Mininco con elevados
contenidos de material fino como arcillas y/o limos. Esta unidad se extendería por más de
150 m de profundidad correlacionándose de buena forma con las captaciones profundas
analizadas en el sector.
135
Figura 4.43: Modelos de inversión 1D de las mediciones que se encuentran en la franja entre Bulnes
y Cabrero.
136
4.3. Gravimetría
En esta sección se presentan los resultados obtenidos para el modelamiento forward 2D de

los datos gravimétricos.
Los puntos de control de cada uno de los modelos corresponden a información integrada
de mediciones realizadas en este estudio, estratigrafía de pozos e información de estudios
previos. El total de los puntos de control utilizados se presenta en la Figura 4.44.
La nomenclatura utilizada para los puntos de control obtenidos en este trabajo se detalla a
continuación:
• LX-n = Cruce con perfil L número X, en punto de gravimetría n.

• DX-n = Cruce con perfil D número X, en punto de gravimetría n.
• MT-x = Punto X de magnetotelúrica.
• TEMx = Punto X de TEM.
• EST = Estimación de profundidad a partir de estudios anteriores, estratigrafía de
pozos cercanos y espesores estimado en otros perfiles gravimétricos.
Los puntos de control utilizados de trabajos previos en la zona de estudio se abreviaron

como se detalla a continuación:
• -DN= DGA, 2018

• -A = Aguirre et al., 2022
Para reportar la estratigrafía de pozos utilizados en la modelación (ANEXO E) se consideró

la siguiente abreviación X0000, donde X puede ser:
• N = ND-0801-
• R = ND-0882-
• M = ND-1602-
• Q = ND-1601-
• W = ND-1603-
137
Los modelos presentados en la sección de resultados poseen la siguiente estructura:
• En el panel superior, la curva en negro representa la Anomalía de Bouguer Completa

medida, la curva en azul es la Anomalía de Bouguer Completa sintética obtenida a
partir de los modelos desarrollados, mientras que la línea magenta corresponde al
campo de tendencia regional. Los puntos en magenta corresponden a los puntos de
quiebre de la tendencia regional definidos durante la modelación.
• En el panel inferior, la curva en rojo corresponde a la topografía del perfil
gravimétrico a partir de mediciones de GPS diferencial, y la curva en negro
corresponde a la base modelada de la cuenca sedimentaria. Superpuesto a esto se
reportan en verde sondajes y modelos de TEM y MT utilizados como espesores
mínimos durante la modelación. Además, en morado se reportan sondajes y
modelos de TEM y MT que detectan la presencia de basamento según la
interpretación realizada.
138
Figura 4.44: Localización de puntos de control utilizados para modelamiento gravimétrico 2D, obtenidos de este estudio, capaciones subterráneas y estudios
previos. Fuente: Elaboración propia.
139
4.3.1. Perfiles E-W (L1-L7)
• Perfil L1
La sección L1 (Figura 4.45) corresponde a un perfil de aproximadamente 90 km de largo

que cruza las comunas de San Nicolás, San Carlos y San Fabián, atravesando en su
extensión las localidades de La Maravilla, San Carlos, Mutupin, El Sauce, Tres Esquinas,
Cachapoal y San Fabián de Alico. La cuenca posee un espesor máximo de 2.080 m, el que
se encontraría en las proximidades de la ciudad de San Carlos. El espesor medio de la
sección es de 585 m.
• Perfil L2

que cruza las comunas de Ninhue, San Nicolás, San Carlos, Chillán, Coihueco y San Fabián
atravesando en su extensión las localidades de Ninhue, Lajuelas, Santa Juana, San Nicolás,
El Manzano, Villa Illinois, Cato y San Fabián de Alico. La cuenca posee un espesor máximo
de 1.525 m, el que se encontraría formando dos depocentros, uno en las cercanías de Villa
Illinois, y el otro a 13 km al oeste de Cato. El espesor medio de la sección es de 710 m.
Debido a la imposibilidad de encontrar afloramientos de roca en su extremo oriental

comparte el amarre terminal con el perfil L1.
• Perfil L3

que cruza las comunas de Portezuelo, Chillán, Chillán Viejo y Pinto, atravesando en su
extensión las localidades de Confluencia, Quinchamalí, Colliguay, Rucapequén, Nebuco,
Quilmo y Chaul. La cuenca posee un espesor máximo de 1.110 m, el que se encontraría en
Quilmo. El espesor medio de la sección es de 300 m.
140
Figura 4.45: Modelamiento e integración de información de otras metodologías del perfil L1.
141
Figura 4.46: Modelamiento e integración de información de otras metodologías de perfil L2.
142
143
• Perfil L4

que cruza las comunas de Quillón, Bulnes, San Ignacio y El Carmen, atravesando en su
extensión las localidades de Quillón, Bulnes, San Javier, Quiriquina, San Antonio y La
Capilla. La cuenca posee un espesor máximo de 1.120 m, el que se encontraría en Quilmo.
El espesor medio de la sección es de 300 m.
Al no encontrar una profundidad de basamento en el sector Bulnes, lugar donde el campo

regional no estaba bien definido, se utilizó un espesor sedimentario inferido de 200 m, el
cual se encuentra acotado por el espesor sedimentario reportado por Aguirre et al. (2022)
y expedientes de pozos cercanos.
• Perfil L5

que cruza las comunas de Quillón, Bulnes, San Ignacio y El Carmen, atravesando en su
extensión las localidades de Canchillas, Calle Larga, Santa Clara, Pueblo Seco, El Carmen
y Loma Blanca. La cuenca posee un espesor máximo de 950 m, el que se encontraría en
Pueblo Seco. El espesor medio de la sección es de 375 m. Debido a la imposibilidad de
encontrar información respecto a la profundidad de basamento en el sector de Canchillas,
lugar donde el campo regional no se encontraba correctamente definido, este tramo no
pudo ser modelado, y se muestra difuminado en la figura.
• Perfil L6
La sección L6 (Figura 4.50) corresponde a un perfil de unos 75 km que cruza las comunas
de Quillón y Pemuco, atravesando las localidades de General Cruz, San Pedro, Pemuco y
Cartago. La cuenca posee un espesor máximo de 1.470 m, el que se encontraría en las
cercanías de Cartago. El espesor medio de la sección es de 625 m.
144
145
146
147
• Perfil L7
La sección L7 (Figura 4.51) corresponde a un perfil de unos 75 km que cruza las comunas
de Yumbel, Cabrero y Yungay, atravesando las localidades de Cabrero, Colicheo, Yungay.
La cuenca posee un espesor máximo de 800 m, el que se encontraría 10 km al oeste de
Yungay. El espesor medio de la sección es de 290 m.
El campo regional de este perfil se encontró mal definido en los sectores de Cabrero y Los
Mayos (Yungay), debido a la inexistencia de información de profundidad de basamento. Por
esto, se consideraron espesores estimados de 80 y 160 m respectivamente, teniendo en
consideración los espesores mínimos definidos por los expedientes de captaciones
circundantes.
148
149
4.3.2. Perfiles auxiliares (D1-D14)
• Perfil D1
La sección D1 (Figura 4.52) corresponde a un perfil de unos 50 km que cruza las comunas
de San Carlos y Coihueco, atravesando las localidades de Quinquegua, Quillahue, Mutupin,
Las Arboledas y Cato. La cuenca posee un espesor máximo de 1.660 m, el que se
encontraría en las proximidades de Quinquehua. El espesor medio de la sección es de 900
m.
• Perfil D2
de San Carlos, San Nicolás, Chillán, Chillán Viejo, Bulnes, San Ignacio, El Carmen, Pemuco
y Yungay, atravesando las localidades de Ranchillos (San Carlos), El Manzano, Chillán,
Chillán Viejo, Montenegro, Canterrana, La Hermita y Ranchillos (Yungay). La cuenca posee
un espesor máximo de 1.275 m, el que se encontraría en las proximidades de Ranchillos
(San Carlos). El espesor medio de la sección es de 840 m.
• Perfil D3
de San Carlos, Chillán, Chillán Viejo, San Ignacio, El Carmen y Pemuco, atravesando las
localidades de Quinquegua, Ninquihue, Chillán, San Ignacio, El Carmen y Santa Filomena.
La cuenca posee un espesor máximo de 1.790 m, el que se encontraría al norte de
Ninquihue. El espesor medio de la sección es de 1.030 m.
• Perfil D4
de San Carlos, Chillán, Coihueco, Pinto y San Ignacio, atravesando las localidades de
Mutupin, Monte Blanco, Talquipen, Las Mariposas, San Ignacio. La cuenca posee un espesor
máximo de 1.525 m, el que se encontraría al 8 km al norte de Talquipen. El espesor medio
de la sección es de 1.045 m.
150
Figura 4.52: Modelamiento e integración de información de otras metodologías de perfil D1.
151
152
153
154
• Perfil D5
de San Carlos, Coihueco, Pinto, San Ignacio, El Carmen, Pemuco y Yungay, atravesando
las localidades de Cachapoal, Cato, Niblinto, Coleal, Ciruelito, La Capilla, Camarico y El
Saltillo. La cuenca posee un espesor máximo de 1.035 m, el que se encontraría en Cato. El
espesor medio de la sección es de 485 m.
155
156
4.3.3. Mapas de espesor sedimentario y profundidad de basamento
En base a los 24 perfiles gravimétricos 2D modelados, se integraron los espesores

sedimentarios y profundidad de basamento obtenidos para la ubicación de cada estación
gravimétrica, con lo que se obtuvo una visión general en planta de la geometría de la
cuenca que se presenta en las Figura 4.57 y Figura 4.58.
En estas es posible notar la existencia de un depocentro ubicado entre las ciudades de San
Carlos y Coihueco, que posee una orientación NNO-SSE. El depocentro principal, ubicado
al oeste de San Carlos, en la localidad de Santa Margarita alcanzaría espesores superiores
a 1.600 m. Este se extiende hacia el SE, encontrándose un espesor sedimentario >1.400
m al norte de Talquipen. Además, al sur de San Carlos se observa un depocentro secundario
entre las localidades de Chillán y Quiriquina que posee espesores sedimentarios por sobre
los 1.000 m.
Por otra parte, se observa una disminución del espesor sedimentario de la cuenca de Itata
hacia el sur, resolviéndose espesores sedimentarios de 800 m en las proximidades de
Pemuco.
Una diferencia importante con respecto a todos los estudios anteriores es que en este
trabajo se observa la existencia de un depocentro aislado en el sector de Cartago que
alcanza espesores sedimentarios por sobre los 1.000 m. Esta estructura es observada
claramente en la anomalía de Bouguer completa, pero debido a su posición distal con
respecto al Valle Longitudinal, es posible que la densidad utilizada para la modelación no
sea representativa de los sedimentos existentes en la zona, lo que podría resultar en una
diferencia en el espesor sedimentario modelado.
Es importante destacar que las secciones presentadas en este apartado corresponden a los
modelos obtenidos a partir de la utilización de parámetros de modelamiento específicos
detallados en la sección 3.1.3. Si bien, es probable que los valores efectivos de espesor
sedimentario y profundidad de basamento puedan variar en función de la elección de los
parámetros de modelación, los rasgos generales de la geometría de la cuenca, que se
derivan de la geometría de los datos gravimétricos obtenidos, se mantienen. Esto es, la
existencia de una cuenca sedimentaría con sus mayores profundidades en el sector norte
de la zona de estudio, y que se vuelve más somera hacia el sur.
157
Figura 4.57: Visión general en planta del espesor sedimentario.
158
Figura 4.58: Visión general en planta de la profundidad del basamento.
159
4.3.4. Comparación con estudios anteriores
• DGA (2011)
El estudio efectuado por DGA (2011) obtiene una geometría de cuenca (Figura 4.59) similar
a la resuelta por este estudio. Las mayores diferencias serían que el depocentro ubicado en
Coihueco no tiene continuidad hacia el depocentro de San Carlos, y la anomalía resuelta en
este estudio se extiende aún más hacia el oeste, llegando a la localidad de Las Mariposas.
Además, al sur de Pemuco reporta el inicio de una profundización de la cuenca que en el
presente estudio no se logra observar.
• DGA (2012)
El estudio DGA (2012) se enfoca en la zona sur de la cuenca de Itata (Figura 4.60), además
de incluir la región del Biobío. En este trabajo se observa la existencia de un depocentro al
oeste de la ciudad de San Ignacio que coincide con la estructura propuesta en el presente
estudio, aunque la profundidad sugerida es mayor (entre 1.200 y 1.500 m) a la resuelta
en la presente publicación (entre 1000 y 1.200). Por otro lado, en la zona de Pemuco DGA
(2012) exhibe un espesor máximo de 300 m, manifestándose como la región de menor
espesor sedimentario dentro de la Depresión Intermedia al sur de la Cuenca de Itata. Esta
estructura difiere con la obtenida en el presente trabajo, donde en Pemuco se resuelven
espesores sedimentarios de entre 800 y 1.000 m. Finalmente, DGA (2012) es capaz de
generar una estimación de profundidades en el área de Tucapel, zona que no pudo ser
cubierta por el presente trabajo. En esta zona, exhiben la presencia de un depocentro
aislado ente las localidades de Huepil y Tucapel, el cual tendría espesores sedimentarios
por sobre los 1.400 m en Huepil. Esto dista fuertemente de los 400 a 600 m obtenidos para
la misma localidad en este estudio.
160
Figura 4.59: Visión en planta de espesor sedimentario estimado por DGA (2011). Imagen de elaboración propia, a partir
de datos de DGA (2011).
161
Figura 4.60: Visión en planta de espesor sedimentario estimado por DGA (2012). Imagen modificada de DGA (2012).
162
• DGA (2015)
El mapa de espesores sedimentarios expuestos por el estudio DGA (2015) (Figura 4.61)
presenta una alta rugosidad que no se condice con las tendencias exhibidas por el mapa de
anomalía de Bouguer completa residual (Figura 4.62) obtenido en el mismo estudio. Un
ejemplo de esto es que la anomalía de Bouguer residual presenta anomalías negativas al
norte de la zona de estudio, sin embargo, esto no se condice con el modelo de espesor
sedimentario obtenido en el sector. Por esta razón este estudio no será considerado en el
análisis efectuado en este trabajo.
163
Figura 4.61: Modelo de inversión de relleno sedimentario obtenido por DGA (2015). Imagen obtenida directamente de DGA (2015).
164
Figura 4.62: Anomalía de Bouguer completa residual obtenida por DGA (2015). Obtenida directamente de DGA (2015).
165
• Aguirre et al., 2022
Por otro lado, Aguirre et al. (2022) determina la existencia un depocentro más prominente
desde Chillán hacia el sur (Figura 4.63), que alcanza espesores sobre los 1.700 m y que no
se encontraría conectado hacia el norte con el depocentro ubicado en la ciudad de San
Carlos. Además, el depocentro de Coihueco se encuentra ligeramente más al sur,
ubicándose sobre la localidad de Talquipen.
En general, los resultados obtenidos en este trabajo se correlacionan de buena manera con
las profundidades y formas presentadas por Aguirre et al. (2022). Además, dada la cantidad
de mediciones realizadas ha permitido obtener un mapa de mejor resolución y mayor
extensión de la forma de la cuenca.
166
Figura 4.63: Modelo de profundidad del basamento estimado por Aguirre et al., 2022. Imagen obtenida directamente de Aguirre et al., 2022.
167
• Análisis conjunto de DGA (2011), DGA (2012) y el presente estudio
Considerando las diferencias de los espesores sedimentarios obtenidos en este trabajo en

comparación con DGA (2011) y DGA (2012), al sur de la zona de estudio, específicamente
desde Pemuco hacia el sur, y el sector oriental entorno a Tucapel y Huepil, se realizó un
análisis conjunto de las distintas bases de datos gravimétricas. Para esto, se realizó una
nivelación de los datos utilizando puntos de control, que corresponden a aquellos puntos
gravimétricos que se traslapan con mediciones del presente trabajo. Posteriormente, se
calculó la diferencia de la magnitud de gravedad absoluta observada en los puntos de
control y se consideró el promedio de ellas como el valor de desplazamiento de las bases
de datos. En el caso de DGA (2011), la base de datos fue desplazada alrededor de 80 mGal,
mientras que DGA (2012) fue desplazada alrededor de 1 mGal. Sin embargo, para esta
última se utilizó un único punto de traslape, ubicado en los alrededores de Cabrero, lo cual
se considera una cantidad deficiente para hacer una correcta nivelación.
La base de datos conjunta de gravedad absoluta fue reprocesada hasta la obtención de la

anomalía de Bouguer completa. Luego, considerando los puntos de control de profundidad
de basamento cero utilizados en este trabajo, y sumando los puntos de control de los otros
estudios, se hizo el cálculo de la tendencia regional para obtener la anomalía de Bouguer
completa residual conjunta (Figura 4.64). En esta figura se puede observar la presencia de
una anomalía gravimétrica negativa entre las localidades de Campanario, Cholguán y Los
Molinos, lo cual corrobora la profundización de la cuenca hacia el sur de la zona de estudio,
que se observaba en DGA 2011. En base a esto, y considerando las limitantes del análisis
realizado por la falta de puntos comunes entre las bases de datos, se recomienda a futuro,
realizar una campaña gravimétrica específica para poder hacer una correcta nivelación de
los datos existentes, y consecuentemente dar paso a la posibilidad de una modelación en
que se puedan integrar todas las mediciones gravimétricas.
168
Figura 4.64: Anomalía de Bouguer Completa residual de la base de datos gravimétrica conjunta del
presente trabajo, DGA (2011) y DGA (2012). Círculos negros: presente estudio; Triángulos blancos:
DGA (2011); Triángulos azules: DGA (2012).
169
5. CONCLUSIONES Y RECOMENDACIONES
- Respecto a los métodos geofísicos utilizados, se corroboró su buena integración. En

particular los TEM y MT permitieron generar amarres respecto a la profundidad del
basamento en lugares donde no existía roca aflorando para medir gravedad, de esta
forma los perfiles de gravimetría tomaron más confiabilidad respecto a su
modelación.
- Para los métodos electromagnéticos utilizados (MT y TEM), fue clave mantener
distancia de las fuentes de ruido como: poblados, gaseoductos y líneas de alta y
media tensión.
- Respecto a la información generada con TEM, esta fue clave en la cordillera de la
costa para conocer la profundidad del basamento rocoso, así como también para
conocer los valores de resistividad reales en superficie y poder discriminar la
presencia de static shift en las curvas de MT.
- Por otro parte, los resultados de MT a nivel de depresión intermedia, muestran en
general una zona conductiva en profundidad. Esta zona se puede interpretar como
el miembro inferior de Fm. Mininco, así como también podría corresponder a alguna
unidad geológica no descrita, por ejemplo, algunas areniscas marinas, o parte del
basamento metamórfico muy alterado y/o fracturado.
- En la zona precordillerana de Los Andes, fue posible reconocer con las mediciones
MT, unidades que pueden pertenecer a Fm. La Montaña (Malleco?), a Fm. Mininco y
a la Fm. Cura-Mallín.
- Respecto al modelado de la cuenca a nivel regional, se percibe que esta unidad
conductora (de baja resistividad eléctrica) en profundidad, podría representar un
sello hidrogeológico, debido a la presencia de materiales compactados y/o de
granulometría muy fina (limos y/o arcillas).
- El modelo de espesor sedimentario obtenido a partir de la gravimetría sugiere la
existencia de un depocentro ubicado entre las ciudades de San Carlos y Coihueco
que alcanzaría un espesor máximo de 2.080 m a 3.000 m al oeste de San Carlos.
La cuenca hacia el sur se vuelve menos profunda, alcanzando un espesor de 800 m
aproximadamente, en las proximidades de Pemuco.
170
− En base a la información de espesor sedimentario obtenida a partir de los modelos
de este trabajo y la información expuesta por estudios anteriores, se recomienda
realizar campañas de medición gravimétrica que se extiendan hacia el sur de la zona
de estudio, abarcando las localidades de Tucapel y Cholguán, con la finalidad de
definir de mejor manera el inicio de la profundización de la cuenca hacia el sur
observada por DGA (2011) y DGA (2012). Por otra parte, se recomienda extender
el perfil gravimétrico D12 entre las localidades de Culenar y San Carlos, para definir
de mejor manera la geometría del depocentro en esa zona.
− Considerando la existencia de diversos estudios independientes a la fecha en la
cuenca de Itata, se recomienda a futuro realizar un reprocesamiento de todos los
datos existentes, que permita unir todas las bases de datos y hacer una modelación
conjunta que logre integrar la información en un solo modelo regional. Esto teniendo
en cuanta las limitaciones técnicas y supuestos efectuados durante el desarrollo de
cada uno de los trabajos realizados previamente. Para conseguir mejores resultados,
también se recomienda en esta fase, realizar mediciones adicionales que permitan
anclar y unificar las distintas bases de datos. Además, se sugiere la realización de
una inversión gravimétrica tridimensional que resuelva de mejor manera tanto la
geometría de la cuenca como la de estructuras más profundas.
− Adicionalmente, y ya que los modelos de espesor sedimentarios obtenidos a partir
de una modelación 2D pueden variar en función de la elección de los parámetros de
modelación, (contraste de densidad entre basamento y sedimento escogido), se
recomienda realizar una campaña de sondajes con una distribución que sea capaz
de muestrear efectivamente las variaciones latitudinales y longitudinales de las
densidades de los rellenos sedimentarios de la cuenca, para precisar los parámetros
de modelación elegidos en cada sector.
- Para futuros estudios geofísicos a realizar en la cuenca Itata, se recomienda realizar
mayor densidad de mediciones Magnetotelúricas (MT) al sureste del sector El Ala,
con el fin de detectar la profundización del basamento hacia el valle. Además, para
complementar los sectores del valle donde la cantidad de información geofísica es
mínima, se recomienda realizar mediciones Magnetotelúricas (MT) de larga duración
(>24 hrs) en los sectores entre Pemuco y Yungay (donde se debe considerar pedir
permisos forestales), al oeste de Cartago (donde se detectó un depocentro en este
estudio), en el sector de La Hermita y al oeste de la localidad de Coihueco. Esto,
con el fin de obtener información de calidad dentro de los primeros 5 km de
profundidad.
171
- En base a los modelos gravimétricos obtenidos y las faltas de puntos de control en
algunos sectores para definir de buena forma el campo regional, se recomienda
realizar estudios geofísicos a futuro que permitan definir la profundidad de
basamento en sectores específicos como Santa Clara (perfil L4), General Cruz (perfil
L5), Cabrero (perfil L7) y Los Mayos (perfil L7), para constreñir de mejor manera la
modelación de los perfiles involucrados. Considerando que estas zonas poseen alto
nivel de ruido electromagnético, una alternativa a los métodos EM es la realización
de estudios sísmicos.
- Por otra parte, considerando el depocentro observado en este estudio en el sector
de Cartago, que no se había observado previamente en estudios anteriores, se
recomienda la medición de estaciones MT de banda ancha sobre la zona profunda
del depocentro para constreñir de mejor manera la profundidad de esta subcuenca.
- Finalmente, también se recomienda complementar los resultados obtenidos por los
modelos gravimétricos en la zona central de la cuenca, con estaciones de
magnetotelúrica banda ancha, que permitan observar la profundidad de basamento
en base a otras metodologías.
172
6. REFERENCIAS
Aguirre, L., Hervé, F., Godoy, E. 1972. Distribution of metamorphic facies in Chile. An
outline. Kristalinikum 9: 7-19.
Aguirre, I., Maringue, J., Santibáñez, I., Yáñez, G. 2022. El rol de la exploración geofísica
en acuíferos profundos en ambientes semiurbanos y rurales en cuencas andinas de ante
arco: caso de estudio en acuífero del río Ñuble, valle central de Chile. Andean geology
49(1): 18-54.
ABEM. 2021. WalkTEM 2 transient electromagnetics (TEM). Abridged Manual. 20 p.
Ajithabh, K., Patro, P. 2023. SigMT: An open-source Python package for magnetotelluric
data processing. Computers & Geosciences 171 (105270).
Anavalón, A. 2018. Evolución tectonosedimentaria de los depósitos Plio-Pleistocenos del

Valle Central de Chile, entre los 36° 50' S y 38° 50' S. Memoria para optar al Título de
Geólogo. Concepción. Universidad de Concepción. Departamento Ciencias de la Tierra. 133
pp.
Blakely, R. 1995. Potential Theory in Gravity & Magnetic Application. Cambridge. Cambridge
University Press. 441 p.
Caldwell, T., Bibby, H., Brown, C. 2004. The magnetotelluric phase tensor. Geophysical
Journal International 158: 457-469.
Chave, A. y Jones, A. 2012. The magnetotelluric method. Theory and practice. Cambridge.
Cambridge University Press. 552 p.
Dentith, M. y Mudge, S. T. 2014. Geophysics for the mineral exploration geoscientist.

Cambridge. Cambridge University Press. 454 p.
DGA, 2009. Estudio Mapa Hidrogeológico Escala 1:1.000.000.
DGA, 2012. Estudio Hidrogeológico Cuenca Biobío. Tomo I, Informe Final y Planos. S. I. T.
N° 297. Santiago.
DGA, 2011. Estudio hidrogeológico Cuencas Bío Bío e Itata. Tomo I. S.I.T. N° 258,
Santiago.
173
DGA, 2013. Modelación Hidrogeológica Cuenca Itata-Bajo. Región del Biobío. S.I.T. N° 320,
Santiago.
DGA, 2014a. Levantamiento hidrogeológico en cuencas pluviales costeras en la Región del

Libertador Bernardo O'Higgins, del Maule y del Biobío, Etapa 2. S.I.T. N° 346, Santiago.
DGA, 2014b. Estimación Preliminar de las Recargas de Aguas Subterránea y Determinación

de los Sectores Hidrogeológicos de Aprovechamiento Común en la Cuenca de las Regiones
del Maule, Biobío, La Araucanía, Los Ríos y Los Lagos. SDT N° 367, Santiago.
DGA, 2015. Análisis y aplicación metodología para la delimitación y sectorización de

acuíferos en la provincia de Ñuble. VIII Región. S.I.T. N° 377, Santiago.
DGA, 2018. Prospección geofísica sector rural déficit hídrico Valle del Itata. S.I.T. N° 433,
Concepción.
DGA, 2018. Sectorización hidrogeológica del acuífero del Río Laja, Región del Biobío. S.I.T.
N° 430, Santiago.
DGA, 2023. Análisis y caracterización del Acuífero de la cuenca del Rio Itata, Estudio
geofísico Itata Bajo. S.I.T N° 532, Chillan.
Instituto de Investigaciones Geológicas, 1981. Mapas geológicos preliminares de Chile -

Hoja Los Ángeles – Angol. Santiago. Escala 1:250.000.
Instituto de Investigaciones Geológicas, 1981. Mapas geológicos preliminares de Chile -

Hoja Concepción – Chillán. Santiago. Escala 1:250.000.
Kane, M. F. 1962. A comprehensive system of terrain corrections using a digital computer.

Geophysics 27(4): 455-462.
Krieger, L., Peacock, J. 2014. MTpy: A Python toolbox for magnetotellurics. Computers &
Geosciences 72: 167-175.
Longman, I. M. (1959). Formulas for computing the tidal accelerations due to the moon
and the sun. Journal of Geophysical Research 64(12): 2351-2355.
Lowrie, W. y Fichtner, A. 2020. Fundamentals of geophysics. Cambridge. Cambridge

University Press. 390 p.
174
Metronix. 2017. ADU-08e. Installation and operating manual. 168 p. Braunschweig,
Alemania.
Metronix. 2002. PROCMT. User’s guide. 100 p. Braunschweig, Alemania.
Nagy, D. 1966. The prism method for terrain corrections using digital computers. Pure and
applied geophysics 63: 31-39.
Platz, A., Weckmann, U. 2019. An automated new pre-selection tool for noisy
magnetotelluric data using the Mahalanobis distance and magnetic field constraints.
Geophysical Journal International 218: 1853-1872.
Rodi, W., Mackie, R. 2001. Nonlimear conjugate gradients algorithm for 2-D magnetotelluric
inversion. Geophysics 66 (1): 174-187.
SERNAGEOMIN. 2003. Mapa Geológico de Chile: Versión Digital. Servicio Nacional de

Geología y Minería, Publicación Geologica Digital, No. 4 (CD-ROM, versión 1.0). Santiago.
escala 1:1.000.000.
Shuey, R. T., Pasquale, A. S. 1973. End corrections in magnetic profile interpretation.

Geophysics 38(3): 507-512.
Ward, S., Hohmann, G. 1987. Electromagnetic theory for geophysical applications. En:
Nabighian, M. (Ed.), Electromagnetic methods in applied geophysics. Volume 1, Theory (p.
131-311). Society of Exploration Geophysicists.
Won, I. J., Bevis, M. 1987. Computing the gravitational and magnetic anomalies due to a
polygon: Algorithms and Fortran subroutines. Geophysics 52(2): 232-238.
Zhang, P., Roberts. R., Pedersen, L. 1987. Magnetotelluric strike rules. Geophysics 52 (3):
267-278.
175

GEO6069v1

Cargado por

Información del documentohacer clic para expandir la información del documentoAcuífero rio Itata 2024 DGA - 1

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

GEO6069v1

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

GEO6069v1

Cargado por

Copyright:

Formatos disponibles

GOBIERNO DE CHILE

MINISTERIO DE OBRAS PÚBLICAS

ANÁLISIS Y CARACTERIZACIÓN DEL

Chillán, marzo de 2024

Director Regional de Aguas de Ñuble

Geo Resource LTDA.

Geofísico Felipe Reyes Córdova (MSc)

Geofísico Andrés Stuardo Espinoza

Geofísica Daniela Montecinos Cuadros (MSc)

Geofísico Daniel González Parada

Geólogo Ian Espinoza Gutiérrez

Lic. en Geología Javiera Venegas González

Lic. en Geología Cristian Roncagliolo Belmar

Estudio financiado por el Gobierno Regional de Ñuble

1.4. CONCLUSIONES DEL ANÁLISIS DE ANTECEDENTES 18

2.2. UNIDADES HIDROGEOLÓGICAS 22

2.3. CAPTACIONES PROFUNDAS EN LA CUENCA DEL ITATA 24

3.3. TRANSIENTE ELECTROMAGNÉTICO 62

4.2. TRANSIENTE ELECTROMAGNÉTICO 131

Para esto, se realizaron un total de 535 mediciones de gravimetría, 100 mediciones de

La información obtenida con TEM en la Cordillera de la Costa, permitió reconocer la

Los modelos de MT, ubicados principalmente en la Depresión Intermedia muestran en

Finalmente, el modelo de espesor sedimentario obtenido a partir de la gravimetría sugiere

Estos resultados son consistentes con algunos modelos de espesor sedimentario

La Dirección General de Aguas (DGA) es el organismo estatal encargado de gestionar,

El sostenido crecimiento económico y desarrollo social de las últimas décadas ha generado

La escasez hídrica se ha agudizado en los últimos años. En efecto, se observa una

La Dirección General de Aguas actualmente posee una Red de Medición de Aguas

El desarrollo de modelos de simulación que permitan planificar y gestionar los recursos

En la actualidad, para iniciar el proceso de búsqueda y generación de respuestas y

Sin embargo, para el correcto desarrollo e implementación de dicha herramienta y/o

El presente informe hidrogeológico desarrollado por GEO RESOURCE LTDA. a solicitud de

La Región de Ñuble comprende una superficie de 13.104,41 km 2, de los cuales 11.326,37

• Recopilar y analizar antecedentes geológicos y geofísicos de la zona de estudio.

La Fm. Mininco es definida en este documento como sedimentitas lagunares y fluviales

• Mapas geológicos preliminares de Chile, Hoja Concepción – Chillán. Escala

La unidades geológicas y sistemas estructurales presentes en este estudio permiten

• Mapa Hidrogeológico Nacional (DGA, 1989)

Mapa hidrogeológico a nivel nacional, donde se identifica de forma general la importancia

Para el área de estudio se identificaron principalmente unidades de alto potencial

Hacia los extremos este y oeste, la presencia de la Cordillera Principal y Cordillera de la

• Mapa Geológico de Chile: Versión Digital (SERNAGEOMIN, 2003)

Mapa geológico de escala 1:1.000.000 que proporciona una representación conceptual de

Adicionalmente cuenta con información de los principales dominios morfoestructurales a

Para la zona de estudio se destaca la presencia de secuencias sedimentarias y volcánicas

• Estudio hidrogeológico cuencas del Biobío e Itata (DGA, 2011)

El objetivo general del estudio consistió en el levantamiento de información hidrogeológica

En el estudio se presenta una caracterización y un modelo conceptual a partir del análisis

Dentro de los resultados de interés, se destaca la campaña geofísica realizada donde se

• Estudio hidrogeológico cuenca del Biobío (DGA, 2012)

El objetivo general del estudio fue el levantamiento de la información hidrogeológica y, a

Se recopiló la información de la geología y geomorfología del área de interés, estructurando

Adicionalmente el estudio describe las unidades hidrogeológicas presentes,

A partir del relevamiento de información hidrogeológica, DGA elabora un modelo conceptual

El principal objetivo del estudio corresponde a desarrollar una modelación conceptual y

La recopilación de antecedentes y levantamiento de información realizada en el estudio

En la caracterización hidrológica de detalle de la cuenca, se incluye la descripción de la red

La campaña geofísica realizada utilizo el método de sondeos eléctricos verticales (SEV). La