Eh Cobriza Modelamiento
Eh Cobriza Modelamiento
Eh Cobriza Modelamiento
Elaborado por:
Abril – 2019
ÍNDICE GENERAL
ÍNDICE DE CUADROS
Cuadro 3.1-1 Medición de niveles freáticos en piezómetros.................................................. 19
Cuadro 3.3-1 Pruebas Permeabilidad .................................................................................... 22
Cuadro 4.1-1 Conductividad hidráulica y unidades hidrogeológicas ...................................... 27
Cuadro 4.1-2 Parámetros hidráulicos históricos zona del proyecto ....................................... 27
Cuadro 4.1-3 Coeficiente de almacenamiento para la zona del proyecto .............................. 28
Cuadro 4.1-4 Leyenda hidrogeológica para la zona del proyecto .......................................... 31
Cuadro 5.3-1 Estimación de la Recarga para las Metodologías Utilizadas ............................ 37
Cuadro 5.3-2 Estimación de la conductancia......................................................................... 38
Cuadro 5.3-3 Tipo de condición de frontera........................................................................... 39
ÍNDICE DE ANEXOS
LISTA DE MAPAS
1.1 Introducción
1.2 Ubicación
El acceso principal para llegar al proyecto, el acceso es por carretera y se realiza siguiendo la ruta:
Lima – La Oroya; La Oroya - Huancayo, Huancayo – Pampas, y de Pampas a Cobriza siendo esta
ruta la más rápida. Sin embargo, existe una ruta alterna desde Lima, que pasa por Ayacucho, no
obstante, que esta ruta esta mejor conservada, pero es mucho más larga.
En el siguiente Cuadro 1.3-1 se muestran las características de la ruta de acceso hacia la zona
del proyecto.
1.4 Objetivos
2. INFORMACIÓN DISPONIBLE
Consiste en la identificación y registro de los cuerpos de agua superficial, como; ríos y quebradas,
alrededor de los componentes mineros y ambientales, habiendo considerado el tipo y la condición
de las fuentes superficiales, entre estos:
Las unidades hidrográficas pertenecen al río Mantaro que drenan sus aguas por la margen
izquierda conjuntamente con otras principales fuentes de alimentación, con aportes principalmente
de la precipitación.
La Qda. Huaribamba, se origina en la cota de 4600 con el nombre de Qda. Chontahuayjo,
uniéndose aguas abajo por la margen derecha con otra quebrada. recibiendo el nombre de Qda.
Huaribambaca. Por la margen izquierda recibe las aguas de las quebradas Boticayoc,
Jauramayoc, Yuracyacu, Amapato, Huerta Amaru y otras hasta desembocar al río Mantaro.
La Qda. Lindero, es una zona que se ubica en la margen izquierda del río Mantaro y abarca varias
quebradas pequeñas, entre ellas la principal de nombre Lindero. Esta zona se ubica también en
la cuenca hidrográfica del río Mantaro, perteneciendo a la cuenca oriental del Océano Atlántico.
La Qda. Lucumayo, también se ubica en la margen izquierda del río Mantaro y que derivan sus
aguas a la cuenca oriental del Océano Atlántico.
El mapeo hidrogeológico se realizó con el fin de evaluar la condición actual de las fuentes
subterráneas y las características hidráulicas de las unidades hidrogeológicas, procediendo con el
monitoreo de fuentes, la que consistió en la identificación y caracterización de cuerpos de agua
subterránea con descarga natural y artificial, los que descargan directamente a superficie y las que
son emanantes de procedencia subterránea (manantiales); filtraciones, sondajes, piezómetros,
canales, y cuerpos de agua en superficie; captaciones, canales, etc., los que discurren en el área
de estudio y que son directa influencia para la determinación del aporte de recarga sobre la
hidrogeología del área de estudio.
Las pruebas hidráulicas se realizaron en cada perforación ejecutada, eligiendo el tipo de prueba
según las condiciones de la roca atravesada, aplicando el método Lefranc carga constante o carga
variable para suelos y rocas fracturadas (Ver Anexo 4), los cuales se describen a continuación.
Este método consiste en introducir agua al sondeo por el interior de la tubería de perforación hasta
la boca de la tubería, manteniendo el nivel del agua a ese nivel sin que rebalse, luego se registra
el consumo de agua cada minuto por medio de un caudalímetro, durante un periodo de 10 minutos
cuando el consumo es alto, ó entre 20 y 30 minutos en tramos de bajo consumo. Este ensayo se
Consiste en inyectar agua al sondeo por dentro de la tubería de perforación hasta alcanzar la boca
de dicha tubería, creando una sobrecarga hidráulica por encima del nivel freático del pozo, luego
se procede a realizar las medidas de los niveles de descenso con una sonda eléctrica manual a
intervalos de tiempos determinados hasta alcanzar el nivel de agua natural del pozo. Este tipo de
ensayo se realiza en rocas fracturadas a poco fracturadas por debajo del nivel freático y en tramos
donde existe poco consumo de agua.
LF CC 71 72 1.00E-02 Grava
LF CC 80.25 81.25 1.00E-02 Grava
LF CC 91.5 72.5 1.00E-02 Grava
LF CC 99 100.25 1.00E-02 Grava
LFCC: Ensayo de permeabilidad Lefranc carga constante
Fuente: Geoservice, 2014
Elaboración: Propia, 2019
LF NC 24 25 2.70E-03 Roca
LF NC 29 30 1.30E-03 Roca
PA-1
LF NC 34 35 1.50E-03 Roca
LF NC 39 40 3.20E-03 Roca
LF NC 14 15 4.40E-03 Grava
LF NC 19 20 6.19E-04 Grava
PA-2
LF NC 24 25 1.04E-03 Grava
LF NC 29 30 1.40E-03 Grava
LF NC 4 5 5.60E-03 Arena
LF NC 9 10 3.10E-03 Grava
LF NC 14 15 1.04E-03 Grava
LF NC 19 20 6.18E-04 Arena
LF NC 24 25 1.80E-03 Arena
PA-3
LF NC 29 30 1.28E-03 Grava
LF NC 34 35 2.68E-03 Roca
LF NC 39 40 3.58E-04 Roca
LF NC 44 45 1.43E-04 Roca
LF NC 49 50 1.80E-04 Roca
LF NC 4 5 9.62E-03 SP-SM
LF NC 9 10 3.29E-03 SP-SM
PA-4
LF NC 14 15 1.70E-03 SP
LF NC 19 20 2.99E-04 SP
LF NC 24 25 4.90E-03 GP
LF NC 29 30 1.36E-03 Roca
LF NC 4 5 4.90E-03 Grava
LF NC 9 10 2.13E-03 Grava
PA-5
LF NC 34 35 1.27E-03 Grava
LF NC 39 40 1.40E-03 Arena
LFNC: Ensayo de permeabilidad Lefranc Nivel Constante
LFNV: Ensayo de permeabilidad Lefranc Nivel Variable
LU: Ensayo de permeabilidad Lugeon
Fuente: Geoservice, 2014
Elaboración: Propia, 2019
Se realizaron en total 09 aforos puntuales que permiten conocer el uso de las aguas que la mina
cobriza utiliza con fines mineros y poblacional en julio 2018 (ver Anexo 5). Cada uno de los puntos
de monitoreo fue registrado con coordenadas UTM mediante el empleo de GPS portátil (ver Mapa
N° 03, Anexo 5). También se tomó en cuenta el inventario realizado en año 2014 como puntos de
registró los parámetros de campo: conductividad eléctrica (CE), pH, temperatura, oxígeno disuelto
y STD.
Aguas de drenaje de mina Nivel Punto de Monitoreo de Efluente Nivel Aguas que se derivan directamente por la margen
1 AFO-1 20/07/2018 5:40:00 p.m. 566441 8610772 1922 0.204 203.7
0 (en tuberia de descarga) Cero 807 (Drenaje de Minas) derecha del río Mantaro.
Reservorio Almacenamiento Reservorio Co Armado donde se El manantial se encuentra aguas arriba en la
2 - 21/07/2018 10:25:00 p.m. 563661 8611942 2560 deposita el agua del Mantial Huerta - - Minero
Manantial Huerta Amaru margen derecha de la Qda. Huerta Amaru.
Amaru El aforo se realizao en la bocatoma que se ubica
Aforo Qda. Huaribamba (en
3 AFO-2 21/07/2018 11:05:00 a.m. 563290 8611316 2542 Punto de la Bocatoma Qda. Huaribamba 0.767 767.0 en la margen derecha del río Huaribamba, en el
bocatoma)
canal de aducción de forma semitrapezoidal.
Aforo Captación Huaribamba
4 AFO-3 21/07/2018 11:15:00 p.m. 563290 8611316 2542 Captación Agua Fresca a Mina 0.136 136.4 Minero El aforo se realizó en la ventana de captación
(En ventana de captación)
Aforo Manatial Huerta Amaru Aforo del Manantial Huerta Amaru en la Aforo de manantial se efectuó en el reservorio que
5 AFO-4 21/07/2018 11:50:00 a.m. 565660 8610342 2450 0.0046 4.6 Minero
(en reservorio Nivel 51) Estación Nivel 51 se ubica en el nivel 51.
Se aforó en un canal de derivación tipo
Aforo de Efluente minero (en trapezoidal que se ubica en la Planta de
6 AFO-5 21/07/2018 12:00:00 a.m. 566460 8610288 2026 Aforo en Planta de Tratamiento en canal 0.1538 153.8
Planta de Tratamiento) Tratamiento y se deriva por la margen izquierda
del río Mantaro.
Se aforó en una tuberia de descarga de PVC que
Aforo de Efluente poblacional
7 AFO-6 21/07/2018 06:10:00 p.m. 568413 8607888 2204 Planta de Tratamiento/ Tuberia 0.003 3.2 se ubica en la Planta de Tratamiento y se deriva
(en Planta de Tratamiento)
por la margen izquierda del río Mantaro.
Aforo Mantial Larian (en Captación en Reservorio derivado por Poblacional
8 AFO-7 22/07/2018 10:15:00 a.m. 565649 8608824 2619 0.0036 3.6 Se aforó en el reservorio de almacenamiento.
reservorio de almacenamiento) Tuberia Parco
Aforo Mantial Soccoshuayjo (En Captación en Reservorio derivado por
9 AFO-8 22/07/2018 10:20:00 p.m. 565649 8608808 2615 0.0031 3.1 Sin Uso Se aforó en el reservorio de almacenamiento.
reservorio de almacenamiento) Tuberia
Las aguas que vienen del manantial Lucumayo se
Aforo Lucumayo (en Tanque de Captación en Reservorio derivado por
10 AFO-9 22/07/2018 11:30:00 a.m. 568107 8607416 2278 0.0097 9.7 Poblacional aforó en uno de los tres tanques de
Almacenamiento) Tuberia
almacenamiento.
1 ML-01 8 608 197 565 630 0.20 Manantial Lindero (parte alta del Club Parco)
3 T-1 8 607 559 569 041 1.20 Alumbramiento de agua termal, baño Coris
4 Machahuayco 8 607 669 567 248 1.20 Alumbramiento de agua sector Machahuayco
Fuente: Geoservice, 2014
Elaboración: Propia, 2019
Nivel de agua
Código Este Norte Cota msnm. Profundidad (m)
(m)
2.6.1 Geomorfología
Esta unidad geomorfológica presenta un conjunto de colinas con un relieve truncado por la
actividad fluvioglaciar del sistema montañoso de cumbres, este conjunto de colinas circundantes
con relieves moderados se le denomina “Superficie Puna”.
En el aspecto morfodinámico, en el área no se destacan procesos relevantes que puedan tener
un notable impacto en la morfología del ámbito de emplazamiento minero.
El presente proyecto se encuentra regionalmente emplazado en la unidad geomorfológica
denominada Cordillera Oriental (Flanco Oeste). Esta unidad está dividida de manera natural por
un profundo valle labrado por el río Mantaro, de tal manera que se le describe como dos bloques
separados, denominados: Cordillera de Razuhuilca (margen derecha) con su ladera de Marcas-
Chincho y Cordillera de Yanaorco (margen izquierda) con su ladera de cobriza.
Fuente: INGEMMET
Las unidades regionales que contienen al área de estudio están enmarcadas entre las quebradas
Huaribambilla y Ramadapampa, estas han sido identificadas por el INGEMMET y están en el
Boletín A-72 – cuadrángulo 26-ñ y son las siguientes.
Proyecto
3. HIDROGEOLOGÍA
En la zona de estudio no existe un control del reservorio acuíferos de este tipo, como son los
macizos rocosos; el cumple un papel esencial en las discontinuidades (diaclasas, fallas, juntas y
planos de estratificación); al ser planos de debilidad del macizo que facilitan el almacenamiento y
la circulación del agua, actúan como drenes facilitando y acelerando el movimiento de las aguas
las cuales generan caídas de carga hidráulica bruscas.
En época de lluvia la infiltración a través de los materiales del cuaternario es mayor lo que genera
la recarga al acuífero del proyecto de Ampliación de Operaciones Mineras de la UM Cobriza.
Se revisaron los niveles de agua subterráneas medidos en los cuatro (04) manantiales y un (01)
piezómetros que se encuentran en la zona de estudio del Proyecto, los cuales constituyeron parte
del programa de investigación hidrogeológica, ver Cuadro 3.1-1 y el Mapa 3.
Este Norte Ts_Id Nivel Freático Radio Cota Sup Cota Inf Código
El relieve topográfico juega un papel importante en la recarga del acuífero, las áreas de las
microcuencas para el Proyecto de Ampliación de Operaciones Mineras de la UM Cobriza, permite
la recepción de las precipitaciones que ocurren en el área a implementar las operaciones y los
componentes del proyecto de profundización de los Nv 200 y el Nv 100, estas precipitaciones
anuales se encargan de alimentar anualmente al sistema acuífero mediante la infiltración del agua,
a través de los suelos, fracturas y fallas, aunque es deficitario durante los meses de mayo, junio,
julio, agosto y octubre, sin embargo en los demás meses se incrementa, otro proceso de recarga
es producido por los ríos de las quebradas: Microcuencas Huaribamba, Lindero y Lucumayo
presentes en estas microcuencas, a diferencia de las precipitaciones el área de recarga de
infiltración solo es por donde se encuentra el flujo de agua siendo menor a la infiltración que se da
por las precipitaciones, esta recarga se da en los meses de estiaje. Es debido a este proceso de
recarga natural, que existe una napa freática que en general fluctúa, el acuífero aumenta su nivel
freático por las precipitaciones producidas en época de lluvia y va disminuyendo lentamente en
época de estiaje.
Una manifestación importante del proceso de descarga del agua subterránea, es la existencia de
algunas zonas de humedad como manantiales que en forma local indican la proximidad del nivel
freático a partir de la cota 2 200 y es la descarga del mismo.
La principal fuente de alimentación o de recarga del acuífero es:
• Agua proveniente de la infiltración producto de las precipitaciones.
• Agua proveniente de la infiltración de los riachuelos de las quebradas Huaribamba, Lindero
y Lucumayo.
La zona de estudio del Proyecto la recarga y descarga del acuífero son las mismas quebradas y
las labores existentes a la actualidad en la zona de estudio, tal y como se observa en el modelo
conceptual construido para estos fines, ver Mapa 6.
La descarga de las aguas subterráneas del sistema acuífero, se produce de forma natural hacia
la base del cauce del rio Mantaro. El mecanismo de descarga del acuífero no es distinguible a
simple vista ya que el primer estrato del suelo no permite una descarga sobre el talud porque estos
percolan directamente hacia el río, los cuales pasan a formar parte de caudal base de los ríos
adyacentes.
A su vez la existencia de las labores subterráneas en el área de estudio, evidencian la recarga
directa del primer estrato, debido a la gran permeabilidad y conformación de suelos semi
permeables, estos se encuentran en la zona de estudio y estacional en la base del talud del cauce
de las quebradas, donde existen descarga difusas, permitiendo el humedecimiento superficial
estacional, asociados a materiales de la zona de llanuras aluviales conformado por arenas finas y
limos de llanura de inundación de las quebradas adyacentes al área de estudio.
A partir de los datos obtenidos del modelo se han podido interpretar un mapa de flujo subterráneo
para las aguas subterráneas que se muestran en el Mapa 8, donde se observa que la dirección se
concentra en el cauce principal de la quebrada Huaribamba.
Para la interpretación se han considerado los ensayos de permeabilidad puesto que como se había
mencionado en el Cuadro 3.4-1 los resultados para cada tramo y prueba realizada y sus
resultados, el cual determina que estos parámetros son bajos en permeabilidad por debajo de los
primeros 30 metros de profundidad y en función a los datos geológicos e investigaciones
hidrogeológicas aledañas del proyecto, se puede determinar que las variables hidrogeológicas e
investigar y su dinámica potencial del sistema acuífero es medio acuífero - acuifugo.
El acuífero está limitado verticalmente por la roca sana (no fracturada) de la formación del grupo
Tarma, pero en otras se encuentra debajo del granito de cobriza (ambos formados por andesitas).
En época de lluvia la infiltración a través de los materiales del cuaternario (de escaso o nulo
espesor) es mayor lo que genera la recarga del acuífero – acuitardo, y el movimiento del agua es
complejo, pero siempre es por la red de grietas interconectadas y a favor de diferencias de carga
hidráulica o presiones piezométricas y adyacentes a las quebradas.
Definimos permeabilidad como la capacidad de un cuerpo (en términos particulares, un suelo) para
permitir el paso de un fluido (en términos particulares, el agua) sin que dicho tránsito altere la
estructura interna del cuerpo. Dicha propiedad se determina objetivamente mediante la imposición
de un gradiente hidráulico en una sección del cuerpo, y a lo largo de una trayectoria determinada.
El concepto permeabilidad puede recibir también las acepciones de conductividad o transmisividad
hidráulica, dependiendo del contexto en el cual sea empleado.
La permeabilidad se cuantifica en base al coeficiente de permeabilidad, definido como la velocidad
de traslación del agua en el seno del terreno y para un gradiente unitario. El coeficiente de
permeabilidad puede ser expresado según la siguiente función:
K=Q/IA
Donde:
- k: coeficiente de permeabilidad o conductividad hidráulica [m/s]
- Q: caudal [m3/s]
- I: gradiente [m/m]
Slide 1
Slide 3
Slide 2
Con el objeto de conocer en el área de estudio la forma de la superficie piezométrica del manto
estudiado, así como seguir la evolución de las capas acuíferas subterráneas y calcular las
fluctuaciones de flujo subterráneo, se elaboró el plano de hidroisohipsas de cada nivel o capa
(Slide) del nivel freático (ver Mapa 8 y Anexos 6.7), que se ha basado en las medidas de la
profundidad del nivel del agua en reposo del piezómetro y de los manantiales existentes y
monitoreados, cuyo análisis permitió calcular la altura de la superficie piezométrica en metros
sobre el nivel del mar y determinar el gradiente hidráulico, trazar las líneas de corriente, determinar
la dirección y sentido del movimiento de las aguas subterráneas.
En toda el área de estudio el nivel piezométrico se encuentra en equilibrio con la presión
atmosférica dando lugar a un acuífero de tipo confinado. La superficie piezométrica se sitúa
aproximadamente entre los 3 900 msnm y los 1820 msnm en el área de estudio.
Se ha estimado que la recarga se produce principalmente en la zona este del área de estudio, la
que presenta una topografía más elevada y en donde ocurren las mayores precipitaciones del área
(de 300 a 950 mm al año).
4.1.1 Permeabilidad
La transmisividad (T) se refiere a la facilidad del acuífero para transmitir agua a través de su
espesor saturado (b), y está directamente relacionada con la conductividad hidráulica (K), que es
la facilidad de un fluido para moverse a través de los espacios interconectados de los materiales,
y depende tanto del fluido como de los medios que atraviesa. En el Cuadro 4.2-1 muestra valores
de conductividad hidráulica teóricos que podrían ser característicos de las unidades
hidrogeológicas de Struckmeier & Margat (1995) (capítulo 5.2) y la clasificación del tipo de acuífero
según Custodio & Llamas (2001), en función de la permeabilidad.
104 UH-1
103 Muy alta UH-1
UH 1
102
UH-1 y UH-2 Acuíferos buenos a moderados
Alta
UH-1 y UH-2
101
UH-1, UH-2, UH-3
Media
UH-2 y UH-3
1
UH-3
10-1 Baja UH-3
Acuíferos pobres
UH-4 y UH-5
10-2 a 10-4
UH-5
Muy baja
10-4 a 10-9 UH- 5 No acuífero
Fuente: Custodio & Llamas, 2001 y Struckmeier & Margat, 1995
Elaboración: Propia, 2019
Código Este Norte Nivel Freático Cota Sup Cota Inf Kx (m/s) Ky (m/s) Kz (m/s)
Sector 1 565 670 8 608 815 36.15 2 697 2 640 0.0001 0.0001 1.00E-05
Sector 2 565 630 8 608 197 4.6 2 918 2 913 1.12E-05 1.12E-05 1.12E-06
Sector 3 565 575 8 608 977 2.5 2 618 2 616 3.68E-05 3.68E-05 3.68E-06
Sector 4 569 041 8 607 559 0.1 1 962 1 962 1.50E-05 1.50E-05 1.50E-06
Sector 5 567 248 8 607 669 0.5 2 630 2 630 1.06E-07 1.06E-07 1.06E-08
Elaboración: Propia, 2019
UH-1
UH-2
UH-3
La construcción y validación de un modelo numérico implica seguir las pautas de un protocolo que
incluye la selección y verificación del código de modelamiento, el diseño del modelo, la calibración,
análisis de sensibilidad y predicción finalmente. Para el caso de estudio se ha seguido el protocolo
propuesto por Anderson y Woessner, 1992, los pasos se describen a continuación y han sido
esquematizados en la Figura 5.1-1.
1. Establecer el propósito del modelo. El objetivo es determinar que ecuación gobernante es
conveniente utilizar y el código que se seleccionará para resolverá la ecuación.
2. Desarrollar un modelo conceptual del sistema. Se identifican las unidades hidroestratigráficas
y las fronteras del sistema. El propósito de la construcción del modelo conceptual es simplificar
el problema y organizar los datos para que el sistema pueda analizarse fácilmente. Se
recopilan datos de campo, incluyendo información necesaria para realizar un balance de
aguas subterráneas y, para asignar valores a los parámetros del modelo y a sus esfuerzos.
3. Seleccionar la ecuación gobernante y el código computacional. El código es el programa de
computación que contiene un algoritmo para resolver numéricamente el modelo matemático.
La verificación de la ecuación gobernante demuestra que la aproximación describe el proceso
físico que ocurre en el medio poroso y/o fracturado.
4. Diseño del modelo. El modelo conceptual se pone de manera adecuada para el modelado
numérico. Este paso incluye el diseño de la malla, la selección de los pasos de tiempo,
estableciendo condiciones iniciales y de frontera y preliminarmente se seleccionan los valores
para los parámetros del acuífero y capas hidrológicas.
5. Calibración. El propósito de la calibración es comprobar que el modelo pueda reproducir el
flujo y el nivel freático medido en el terreno. Para ello se ajustan los parámetros del modelo
hasta minimizar las diferencias entre las variables calculadas y las observadas, en este caso
caudales de flujo y altura piezométrica.
6. Verificación del modelo. El propósito de la verificación del modelo es establecer una mayor
confianza en el modelo utilizando el conjunto de valores de los parámetros y capas calibrados
para reproducir un segundo conjunto de datos de campo.
7. Predicción. Cuantifica la respuesta del sistema para eventos futuros. El modelo se corre con
los valores de los parámetros y capas calibrados, excepto para las capas en las que se espera
que cambien en el futuro. Para realizar la simulación son necesarias las estimaciones futuras
de las capas. La incertidumbre en una simulación predictiva surge de la incertidumbre en el
modelo calibrado y la incapacidad para estimar los valores exactos de la magnitud y las capas
futuras.
8. Presentación del diseño del modelo y resultados. La presentación clara del diseño del modelo
y los resultados esenciales para una comunicación eficaz de los esfuerzos del modelado.
Se ha realizado la equivalencia del macizo fracturado aun medio poroso continuo que se
fundamenta en el concepto de un medio poroso equivalente al fracturado (Anderson & Woessner,
1990) ampliamente utilizado con buenos resultados.
A continuación, se muestra la ecuación diferencial que define el flujo de agua subterránea, esta se
obtiene a partir de la Ley de Darcy y la ley de conservación de la masa:
𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ
(𝐾𝑋𝑋 ) + (𝐾𝑦𝑦 ) + (𝐾𝑧𝑧 ) − 𝑊 = 𝑆𝑠
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡
Dónde:
H : altura potencial o carga hidráulica
Kxx,Kyy,Kzz : valores de conductividad hidráulica a lo largo de los ejes coordenados x, y, z ejes,
que se supone están paralelos a los ejes principales de la conductividad hidráulica
(L / T);
W : razón de flujo volumétrico por unidad de volumen, y representa a las fuentes y/o
sumideros de agua (T-1);
Ss : almacenamiento específico del material (L-1), y
T : tiempo (T).
Además, se han incluido elementos discretos que representan los componentes del proyecto como
las labores de profundización del Nv 200 y el Nv 100 sobre todo. Por tanto, se trata de un tipo de
modelación en modo permanente que superpone un medio continuo con la presencia de
elementos discretos que conducen las aguas subterráneas a modo de vías preferenciales.
El código de modelamiento es el programa que contiene los algoritmos para resolver mediante
métodos numéricos la ecuación del flujo.
Existen dos grandes grupos de métodos numéricos usados en la resolución de la ecuación del
flujo: el método de los elementos finitos (MEF). La selección del método numérico condiciona en
gran medida la geometría de los elementos (triángulos) que forman la malla. Se podría decir que
en general los MEF conllevan mallas libres de los MEF es totalmente arbitraria. De todas formas,
la selección del método numérico es subjetiva y cada uno de los métodos tiene ventajas y
desventajas. Los MEF tienen una interpretación física (balance de masa de los elementos), una
entrada de datos cómoda y están tradicionalmente muy bien documentados. La ventaja es los
MEF se pueden ajustar perfectamente al contorno, pero su interpretación física es difícil debido a
una mayor complejidad del desarrollo matemático.
Se ha elegido el software Feflow 7.1 que está basado en la técnica de los elementos finitos. Las
ventajas principales de este método y que se adaptan al caso de este estudio hidrogeológico para
el Proyecto de Ampliación de Operaciones Mineras de la UM Cobriza y son las siguientes:
• Mejor representación de estratos con gran buzamiento y anisotropía.
• Permite un refinado local para una mayor precisión sin necesidad de refinar las condiciones
de contorno como ocurre en el caso de los elementos finitos.
• Permite la conservación del balance masas.
• Realiza un balance de aguas de manera fácil.
Las condiciones de frontera o contorno aplicado al modelo numérico, a fin de representar las
condiciones naturales del medio hidrogeológico, incluyeron: carga hidráulica de flujo específico;
carga dependiente, (ver Anexo 6 y el Mapa 6). Estas condiciones de contorno se ilustran en el
Figura 5.3-2 la que se describe a continuación:
Carga específica: (Hydraulic Head BC) Los bordes de carga constante de borde tipo 1 fueron
asignados sobre la superficie superior del modelo a lo largo de los cursos de agua permanentes,
de acuerdo al modelo de elevación digital del terreno. Los límites representaron las entradas y
salidas del sistema, habiendo asignado también esta condición de borde en la parte superior del
modelo.
Para los casos puntuales se le asignó el tipo de borde superficies de filtración (Seepage face)
asignado a las labores mineras para representar el cono de depresión, así como los caudales de
filtración que se producirían al profundizar los niveles proyectados.
Sin flujo: se asignaron condiciones de borde tipo no flujo, para representar divisiones inferidas de
las aguas subterráneas y líneas de flujo en todo el perímetro del modelo.
Se define como el volumen de agua que puede ser liberado por un prisma vertical del acuífero, de
sección igual a la unidad y altura la del espesor saturado, si se produce un descenso unidad del
nivel piezométrico.
En el caso que nos ocupa el régimen del acuífero es de tipo semi confinado y por tanto el software
utilizará el valor indicado en la porosidad efectiva. Como se indicó anteriormente no se ha calibrado
al realizarse solo simulaciones en régimen estacionario.
La condición de borde superior del modelo numérico se corresponde con una recarga distribuida,
este es el único mecanismo de recarga natural al sistema de agua subterránea que se corresponde
a la infiltración de las precipitaciones, para ello se ha tenido en consideración las características
meteorológicas del estudio hidrologico, considerando la precipitación y una vez analizado el
balance hídrico se ha considerado una tasa anual de infiltración de 50 mm y una recarga al acuífero
de la zona de estudio de 142.5 mm/año.
A pesar de que las quebradas Huaribamba, Lindero y Lucumayo, se considera, como aluviales
asociados a ellos pueden estar parcialmente saturados y pueden estar cediendo agua al sistema
regional. Por esta razón han sido considerados en el modelo, aunque no son río totalmente
Donde:
K: permeabilidad en la zona saturada del río L: longitud del río
W: ancho del río
M: espesor saturado del río
En em Cuadro 5.3-2 se muestran las variables y resultados obtenidos de conductancia en cada
uno de los ríos del dominio del modelo.
K L W M Conductancia
Quebradas
m/d m m m m2/d
Feflow permitió realizar tareas de modelado en la zona de estudio para cada capa (Slide) y
determinar las características de las capas del movimiento del flujo subterráneo (a través de
Slide 1 Slide 2
• La geología de mina se caracteriza por su gran exposición de rocas con calizas grises, limo
arcillitas laminadas gris oscuras del grupo Tarma.
• Se realizó una campaña de monitoreo; la 1ra corresponde a la época estiaje (inventario Julio
de 2018) registró 4 puntos de monitoreo superficiales y 01 fuentes subterráneas.
• Se observa que los principales aportantes de agua subterránea en el proyecto son de las
quebradas Huaribamba, Lindero y Lucumayo, la que evidencia que la gran mayoría son aguas
neutras de similar composición química.
• Se diferenciaron 3 unidades hidrogeológicas; con predominancia de Depósitos Cuaternarios
(Q-co/al) cubierto por suelos fluvioglaciarios, Se trata de un sistema hidrogeológico no
complejo, dentro de un medio de rocas de baja a media conductividad hidráulica, localmente
intruida por un sistema de vetas ligado a las anomalías estructurales, que puntualmente
conducen agua subterránea atribuida a su moderada conductividad del Grupo Tarma (Cs-t) y
se el Granito de cobriza (PsTr-g-c) el cual se presenta desarrolada como unidades
hidrogeológicas de Acuifugo
• Los resultados de conductividad hidráulica de han sectorizado en el sector 1 se tiene 1e-3 a
1e-5 en el sector 2 se presenta una conductividad de 1.12e-5 a 1.12 e-6 y el sector 4 se tiene
3.68e-5 a 3.68e-6, en el sector 4 de 1.50e-5 a e-6 y el sector 5 de 1.06e-7 a e-08.
• Se tiene un nivel freático; una de superficie ligado a suelos cuaternarios y flujos superficiales
temporales y otro nivel profundo dentro de un medio “semiconfinado a confinado” de baja
presión, debido a la escasa conductividad hidráulica del macizo rocoso, en la cual se
presentan las actividades de laboreo actual la que ejerce “baja a moderada” carga hidráulica.
El mapa piezométrico muestra isopiezas y direcciones de flujo hacia las labores más
profundas, esta recibe descarga de flujo subterráneo localizado.
• Se ha elaborado un modelo matemático de flujo en régimen estacionario o en equilibrio como
punto de partida para el modelo en régimen transitorio operacional.
• El drenaje que producen las canteras en el sistema hidrogeológico supone un 1.3 % del
balance de agua oscilando entre 0.6 y 8 L/s. El drenaje mayor se genera en el sector de la
cantera Cristina Norte donde la profundización final es mayor alcanzando una cota mínima de
225 msnm.
• No se observan cambios significativos en el comportamiento hidrodinámico de los ríos
Mantaro, con una recarga del acuífero del 142.5 mm/año.
• Se produce un descenso piezométrico muy localizado con un radio de influencia sobre el
laboreo existente, este descenso puede alcanzar hasta un máximo de 80 m en su parte más
profunda.