U2.8b Aguas Subterraneas RevA PDF

UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERIA CIVIL

1
Departamento Académico de Hidráulica e Hidrología
RECURSOS HIDRAULICOS HH333

RECURSOS HIDRAULICOS HH333 G,I,J,K
CIO
Secciones: G,I, J, K
Profesores: Walter Cabrera Cabrera / Marisa Silva Dávila / Carlos
Iparraguirre Ortiz / Ada Arancibia Samaniego
Unidad 2.8b:
UNI - FIC - DAHH
AGUAS SUBTERRÁNEAS
2
CIO
AGUAS SUBTERRÁNEAS
INTERDEPENDENCIA DEL AGUA SUPERFICIAL Y SUBTERRÁNEA
UNI - FIC - DAHH
IDENTIFICACIÓN DEL ÁREA DE RECARGA CON HIDROLOGÍA ISOTÓPICA

INTERDEPENDENCIA
3 CONCEPTO CLAVE El agua de la precipitación y deshielo
que se infiltra en el terreno, fluye en
dirección predominantemente vertical
hasta que llega a la napa acuífera.
CIO
Luego fluye a lo largo de un gradiente

casi horizontal, la mayoría para
descargar naturalmente en algún lugar
aguas abajo, donde el acuífero
intersecta lagos, corrientes y ríos en
zonas vecinas o alejadas.
Por lo tanto, el agua subterránea constituye agua superficial prestada.

UNI - FIC - DAHH
Es solamente una cuestión de tiempo antes de que la recarga natural se convierta en

descarga natural, es decir, en evapotranspiración o agua superficial, en algún lugar aguas
abajo.
Fuente: Ponce, 2010
INTERDEPENDENCIA
4 CONCEPTO CLAVE Si el agua subterránea se extrae y se consume
(no es regresada al acuífero), resulta en un
decremento en la descarga natural, es decir, en
una disminución del flujo de base que alimenta
las aguas superficiales.
CIO
Es importante evaluar el impacto en el flujo

base por su importancia ambiental, y por ello la
recarga y la descarga.
Pero la relación entre el lugar de recarga y el

lugar de la descarga es difícil de determinar. La
circulación del agua subterránea puede ser (1)
UNI - FIC - DAHH
local, (2) subregional o (3) regional.
Estudios geológicos y de calidad de agua ayudan a inferir el origen de las agua

subterráneas (el área de recarga), pero también técnicas especiales como la isotópica.
Fuente: Ponce, 2010
INTERDEPENDENCIA
5 EVOLUCIÓN DE LA CALIDAD DEL AGUA SUBTERRÁNEA
CIO
Fuente: H. Jegat, 2015

UNI - FIC - DAHH
En el agua de recorrido local predominan los carbonatos, luego progresivamente

combinaciones de carbonatos y sulfatos, sulfatos, combinaciones de sulfatos y cloruros, y
en las aguas de mayor recorrido y tiempo de permanencia o flujos regionales, los
cloruros, a la vez que aumenta la salinidad.
HIDROLOGÍA ISOTÓPICA
6 CONCEPTO
Aprovechar un elemento presente en la
masa de agua para estudiar su
comportamiento durante las distintas
CIO
fases del ciclo hidrológico.

Trazadores Naturales:
Técnicas Isotópicas
Isótopos Ambientales
UNI - FIC - DAHH
 Protio ( 1 H, H ) 99,98%, estable

 Deuterio ( 2 H, D ) 0,016 %, estable
 Tritio ( 3 H, T ) 5.10-6, radioactivo
Fuente: H. Jegat, 2015
7 CONCEPTO
Isótopos Ambientales
 Oxígeno O ( 16 O ), 99,76%, estable
Oxígeno 17 ( 17 O ) 0,04%, estable
CIO
 Oxigeno 18 ( 18 O ) 0,20%, estable
Isótopos Estables del Agua
 1H2 18O
UNI - FIC - DAHH
 1H2H 16O
Fuente: Jegat, H.,2015
8 CONTENIDO ISOTÓPICO DE LA LLUVIA
La lluvia en cada lugar de la
tierra tiene un particular
composición isotópica.
CIO
Varios factores, de la
formación de nubes
(temperatura, latitud) y de
precipitación (vientos,
altitud) y el efecto del
mayor peso isotópico,
concentra los mayores
contenidos más cerca de los
UNI - FIC - DAHH
bordes continentales y
decrece hacia el interior de
los continentes. Composición en 2H y 18O de la lluvia
d = variación del contenido respecto al océano estándar
Fuente: Soler, A, 2015
CIO
El contenido isotópico d18O y dD

en toda la tierra se agrupan en
UNI - FIC - DAHH
una estrecha franja lineal

Composición ponderada anual en d18O de la lluvia

Los procesos
hidrológicos
afectan la
CIO
composición
isotópica del agua
UNI - FIC - DAHH

11 APLICACIONES EN HIDROGEOLOGÍA
 Edad
CIO
 Origen
 Tiempo de residencia
 Mezclas
 Recarga
UNI - FIC - DAHH
 Hidrodinámica
 Vías Preferenciales

12 APLICACIONES EN HIDROGEOLOGÍA - EJEMPLO
Definir el modelo
conceptual del
CIO
acuífero del Valle

Quíbor (Venezuela):
Establecer si la
recarga para el
bombeo agrícola
proviene de la
UNI - FIC - DAHH
quebrada Atarigua o
de la parte alta del
río Turbio.

13 APLICACIONES EN HIDROGEOLOGÍA - EJEMPLO
Tres campañas:
Sep 2001: 11 muestras
May 2002: 37 muestras
Oct 2002: 58 muestras

CIO
Análisis de D y 18O
Interpretación cualitativa
y cuantitativa
Resultado:
UNI - FIC - DAHH
Asociación del agua de los

pozos profundos con la
lluvia de la parte alta
(Río Turbio)
14
DETERIORO DE LAS AGUAS

CIO
SUBTERRÁNEAS
INTRUSIÓN SALINA
UNI - FIC - DAHH
SOBREXPLOTACIÓN DE ACUÍFEROS
CONTAMINACIÓN DE AGUAS SUBTERRÁNEAS
15
UNI - FIC - DAHH CIO
INTRUSIÓN SALINA
16
La sobrexplotación en zonas litorales favorece el avance del frente del agua salada sobre
el agua dulce, extrayendo los pozos agua salada y salinizando el subsuelo.
La zona afectada prácticamente no se recupera ni puede explotarse agua dulce.

SOBREEXPLOTACIÓN
17 MÁS PROPIAMENTE ES
CUANDO SE SUPERA LA
EXPLOTACIÓN
SOSTENIBLE, QUE ES
DEL ORDEN DE UN
CIO
PORCENTAJE DE LA
RECARGA (10% - 70%,
PREFERIBLE TECHO DE
40%) O ALREDEDOR DE
2% DE LA
PRECIPITACIÓN EN EL
UNI - FIC - DAHH
ÁREA DE RECARGA.
18
CONTAMINACIÓN DE AGUAS SUBTERRÁNEAS
19 FUENTES PUNTUALES Y DIFUSAS DE CONTAMINACIÓN
PLUMA CONTAMINANTE
20 MOVILIDAD PREDOMINANTE DE CONTAMINANTES
CIO
Principalmente los contaminantes se

UNI - FIC - DAHH
movilizan según el mecanismo de

movilización del agua subterránea Desde los ríos desconectados la
(efecto advectivo) contaminación progresa más lenta por la
ralentización del flujo en medio no saturado
PLUMA CONTAMINANTE
21 CONCENTRACIÓN DE CONTAMINANTES
La concentración de solutos va disminuyendo desde la fuente con el recorrido del flujo, por
causa de varios fenómenos físico químicos. Se debe comparar con los límites máximos
admisibles. Cuando estos se superan el agua se considera contaminada.
22
TRANSPORTE DE CONTAMINANTES
PROCESOS DE TRANSPORTE
23
REQUERIMIENTOS DEL MODELO
 PROCESOS FÍSICOS
 ADVECCIÓN
 DISPERSIÓN
CIO
 SORCIÓN (ADSORCIÓN-DESORCIÓN)
Fronteras de
Transporte
 PROCESOS QUÍMICO-TÉRMICOS Propiedades

de los
 DEGRADACIÓN POR REACCIÓN. Materiales
UNI - FIC - DAHH
Parámetros
de Reacción
FUENTE: GUIGUER, N., 2017.

ECUACIÓN DE TRANSPORTE
24
𝜕 𝜕 𝜕𝐶 𝑞𝑠 𝐶 𝜕𝐶 C = Concentración
− 𝑣𝑖 𝐶 + 𝐷𝑖𝑗 ± 𝐶𝑠 − 𝜆 𝐶 + 𝜌𝑏 =𝑅 [M L-3]
𝜕𝑥𝑖 𝜕𝑥𝑖 𝜕𝑥𝑗 𝜃 𝜃 𝜕𝑡 mgr / L
rb = densidad del suelo

CIO
Advección Dispersión Sumidero/Fuente Reacciones Retardo  = porosidad

R = Factor de retardo
Dij = tensor de difusión
UNI - FIC - DAHH
 La Advección se debe al flujo

 La dispersión y el retardo se deben a las propiedades de los materiales
 Las reacciones son dependientes de las especies
FUENTES Y SUMIDEROS qs
 Cs
25 EJEMPLO DE FUENTES 
Fuente debajo del
Se requiere nivel freático Fuente en
Fuente sobre el Recarga de
saber: el nivel
nivel freático Agua Limpia
freático
CIO
• Tipo de
Frontera
(sobre o
debajo del
nivel freático)
• Localización
UNI - FIC - DAHH
• Historia en el
tiempo
• Concentración

ADVECCIÓN

26 − (viC )
xi
 Transporte conjunto con el flujo
 Fuente dominante del transporte para la mayoría de los sistemas acuíferos
CIO
t Fuente
“flujo conectado”
c x = v t
Pluma
x
UNI - FIC - DAHH
t=0 t=t1
fuente de “pulso” fuente “constante”

DISPERSIÓN
27
  C 
Dij

xi  x j 
Recarga
Fuente
CIO
D = DM + D * acuífero Pluma
Dos procesos
ZONA DE DISPERSIÓN
 Dispersión Mecánica (DM)
por fricción, tortuosidad.
UNI - FIC - DAHH
 Difusión Molecular (D*)

por gradientes de
concentración, tortuosidad

DISPERSIÓN
28   C 
Dij
 La dispersión expande la pluma 
xi  x j 
 No hay movimiento lento del centro de masa
ADVECCIÓN + DISPERSIÓN
CIO
Fuente A+D
c x = v t
Pluma
x
UNI - FIC - DAHH
t=0 t=t1
Fuente de“pulso” Fuente “constante”

ADSORCIÓN Y DESORCIÓN
29
ADSORCIÓN
Ciertos contaminantes atacan por si mismos los granos de suelo
(partición), y se quedan impregnando las partículas sólidas.
CIO
FLUJO
UNI - FIC - DAHH
Nota: ADSORCIÓN es retención de solutos en superficies.

ABSORCIÓN es ingreso de solutos al agua desde mezclas gaseosas.
ADSORCIÓN Y DESORCIÓN
30
DESORCIÓN
Cuando las concentraciones del contaminante en el agua subterránea
caen, las partículas adsorbidas son disueltas de nuevo en el agua.
CIO
FLUJO
UNI - FIC - DAHH

ADSORCIÓN C
R
31 t
la adsorción causa “retardo” en relación al movimiento del agua
subterránea porque parte de los solutos se quedaron impregnados en las
superficies sólidas.
CIO
A + D + ADSORCIÓN
Fuente
c x < v t/R
Pluma
UNI - FIC - DAHH
x
t=0 t=t1
fuente de “pulso” Fuente “constante”
DEGRADACIÓN 𝐶
−𝜆 𝐶 + 𝜌𝑏
32 𝜃
 Es una consecuencia de la reacción, por la remoción de unas especies
químicas y la producción de otras.
 Se degradan solutos disueltos y solutos absorbidos.
CIO
A + D + A + DEGRADACIÓN
Fuente
c x < v t/R
Pluma
UNI - FIC - DAHH
x
t=0 t=t1
fuente de “pulso” Fuente “constante”
33
CIO
MANEJO DE ACUÍFEROS
RECARGA ARTIFICIAL DE ACUÍFEROS
UNI - FIC - DAHH
EXPLOTACIÓN CONJUNTA AGUAS SUPERFICIALES Y

SUBTERRÁNEAS
34 RECARGA INDUCIDA
CIO
• Diques y empozamientos aumentan el área de recarga y favorecen la

UNI - FIC - DAHH
extensión de la infiltración y de las tasas de infiltración.

• Se usa el acuífero como un embalse.
• Se tiene la ventaja adicional la capacidad de depuración de aguas
del acuífero
35 BOMBEO PARA RECARGA ARTIFICIAL
Se alterna bombeo de recarga al acuífero, con bombeo

de extracción del acuífero.
EXPLOTACIÓN CONJUNTA DE AGUAS SUPERFICIALES Y
36
SUBTERRÁNEAS
CIO
Para satisfacer mejor a la demanda en cantidad y calidad se combinan funciones de captación,

UNI - FIC - DAHH
almacenamiento y depuración. Las aguas subterráneas afianzan la garantía por su propiedad de

flujo sostenido.
Por ejemplo en función de la disponibilidad hídrica en un embalse, una determinada demanda
puede ser atendida a partir del agua superficial contenida en el mismo, conjuntamente con aguas
subterráneas o exclusivamente con estas últimas. Este sistema se puede completar con recarga
artificial que se almacene en el acuífero el agua circulante por el río.
37
PRINCIPIOS
VULNERABILIDAD
VULNERABILIDAD
38 DEFINICIONES
•Vulnerabilidad intrínseca
Propiedad referida exclusivamente al medio (tipo de acuífero y cobertura,
permeabilidad, profundidad, recarga, etc.), sin tener en cuenta la incidencia de las
CIO
sustancias contaminantes.
•Vulnerabilidad específica
Además del comportamiento del medio, trascendencia al tipo y carga del contaminante
•Riesgo
(Foster, 1987) definen al riesgo como el peligro de deterioro en la calidad de un
acuífero, por la existencia real o potencial de sustancias contaminantes en su entorno.
UNI - FIC - DAHH
La vulnerabilidad específica incluye parcialmente el concepto de riesgo, toda vez que

se refiere al peligro de deterioro en relación a sustancias contaminantes específicas.
http://tierra.rediris.es/hidrored/ebvulnerabilidad.html
VULNERABILIDAD
39 VULNERABILIDAD Y RIESGO
Con el objeto de clarificar, se cita un ejemplo sencillo:
 a) Un depósito con puertas sin llave ni candado, ubicado en la zona
rural, es muy vulnerable debido al fácil acceso; sin embargo el riesgo
CIO
de hurto o robo es bajo, porque en la región no hay vándalos.

 b) Otro depósito, ubicado en la zona urbana y provisto de fuertes
cerraduras, candados y ventanas enrejadas, es poco vulnerable pero de
alto riesgo, debido a la presencia de delincuentes en la región.
UNI - FIC - DAHH
40
INDICE DRASTIC
VULNERABILIDAD
41
INDICE DRASTIC
VULNERABILIDAD
42
INDICE DRASTIC
VULNERABILIDAD
43
INDICE GOD
VULNERABILIDAD
VULNERABILIDAD
44
PLANIFICACIÓN - ZONIFICACIÓN DE PROTECCIÓN
DE ACUÍFEROS
Estudios de vulnerabilidad y de transporte de hipotéticos contaminantes sirven para definir los perímetros de
protección para salvaguardar la calidad y cantidad de las aguas subterráneas que se extraen de las captaciones;
ante el riesgo potencial de la actividad antrópica en los alrededores, restringiendo de forma gradual las actividades
o instalaciones susceptibles de contaminar las aguas subterráneas, o que puedan afectar al caudal destinado al
abastecimiento.
45
MODELACIÓN NUMÉRICA
MODELACIÓN CONCEPTUAL
MODELACIÓN
ECUACION DE FLUJO SUBTERRÁNEO
46
ECUACIÓN GENERAL
La ecuación general es una ecuación diferencial, obtenida en base a los
principios de conservación de masa y de cantidad de movimiento:
  h    h    h  h
CIO
Kx  +  K y  +  K z  = SS +W
x  x  y  y  z  z  t
Incluyendo factores de reducción y corrección, en función de la humedad,

se expande para aplicar para medios no saturados:
  h    h    h  h S w
UNI - FIC - DAHH
 x rw 
K k + 
 K k 
 +  z rw 
K k = S S +  W
x  x  y  y  z  z  t t
y rw w S
Como toda ecuación diferencial, deben especificarse las condiciones

iniciales y de borde. Hay tres tipos de condiciones, I, II y III.
CONDICIONES INICIALES Y DE BORDE
47
CONDICIÓN TIPO I (Dirichlet) – POTENCIAL IMPUESTO
h( xb , y b , z b , t ) = hb (t )
 Condición sobre la variable h.

CIO
 h es el mismo en el contorno:
línea o superficie
equipotencial.
 Suele estar asociado a
contactos entre el acuífero y
una masa de agua superficial:
UNI - FIC - DAHH
 lagos
 mares
 ríos
48
CONDICIÓN TIPO II (Newman) – FLUJO INCONDICIONAL IMPUESTO
h
q n ( xb , y b , z b , t ) = qb (t ) q n = − K . = qb (t )
n
 Condición relacionada con la
derivada, con el flujo q.
CIO
 Flujo de agua definido que sale

del acuífero o penetra en el a
través del borde.
 Puede ser nulo en el caso del
contacto entre el acuífero y una
unidad impermeable, o un borde
alejado del objeto principal del
modelo.
UNI - FIC - DAHH
 Las divisorias de agua también se

ajustan a este tipo de condición
de contorno.
CONDICIÓN TIPO III (Cauchy) – FLUJO CONDICIONADO AL POTENCIAL
49
q n ( xb , y b , z b , t ) = C b (hb (t ) − h( xb , y b , z b , t ) ) Cb = K b / Db Conductancia
 Condición relacionada con la

derivada, con el flujo q.
CIO
 Se aplica a las entradas y salidas de

agua del acuífero a través de capas
semiconfinantes que lo separan de
otra fuente de recarga externa.
 El flujo que sale del acuífero o
penetra en él depende de la
diferencia del potencial entre el H1(t) = Carga Hidráulica en el reservorio
H2 (en AB, t) = Carga hidráulica por calcular por el modelo
acuífero y la fuente externa, de la
UNI - FIC - DAHH
conductividad hidráulica vertical del Flujo condicionado a través AB

acuitardo o capa semiconfinante, de
su extensión superficial y de su 𝑯𝟏 −𝑯𝟐 𝑲′
espesor. 𝒒= 𝑲′ 𝑨 𝑪= ′
𝒃′ 𝒃
DESARROLLO DEL MODELO
TIPOS DE MODELOS
50 Primero se desarrolla el modelo conceptual y luego este se implementa en el modelo numérico.
El modelo numérico tiene dos fases, no transitorio y transitorio (el último requiere datos en el tiempo).
Modelo Conceptual Modelo Matemático

Estático, descriptivo, cuantifica el medio en Estático y Dinámico, predictivo, soluciona

forma general. ecuaciones y cuantifica los procesos de
Define la dinámica general del medio físico. transporte de agua y de contaminantes
Define las condiciones de borde 50
Fuente: Gomes et al, 1997
MODELOS CONCEPTUALES
Un modelo conceptual es una representación estática del sistema físico:
51
Identifica los elementos principales de las condiciones hidrogeológicas:
▪ Geología:
▪ Geometría, tope, base y espesor de las capas.
▪ Estructuras (i.e. fracturas, fallas, pliegues, dirección/buzamiento, intrusiones, diques, sills, zonas
CIO
de sizalle, etc)
▪ Litologías y su distribución en el espacio (capas, acuñamientos, paleo-canales, abanicos aluviales,
de acuerdo al ambiente de deposición …)
▪ Hidrogeología/Climatología/Hidrología
▪ Acuíferos y acuitardos, su geometría y condiciones de borde (ríos, lagos, tajos, túneles, etc).
▪ Parámetros hidráulicos (K, Sy, Ss, nef, nt, densidad)
▪ Fuentes y sumideros (pozos de extracción/inyección; drenes, etc)
▪ Mecanismos y tasas de recarga y descarga.
UNI - FIC - DAHH
▪ Contaminación
▪ Condiciones de borde: fuentes de contaminación, carga contaminante, ubicación, duración,
intensidad.
▪ Parámetros de migración/atenuación de los contaminantes (partición, dispersión, difusión,
sorción, biodegradación transformaciones.
Fuente: Guiguer, N., 201751
MODELOS CONCEPTUALES
52
Identifica los elementos principales de las condiciones hidrogeológicas:
PARÁMETROS DE ENTRADA
CIO
 Geometría de las capas

 Propiedades hidráulicas (K, Ss, Sy, espesores, nef, etc.)
 Recarga/Descarga (tasas, distribución especial y temporal)
 Condiciones iniciales y de borde
 ríos, lagos, tajos, túneles, etc.
 Fuentes y sumideros (pozos de inyección, extracción, ubicación, tasas, variación temporal)
UNI - FIC - DAHH
 Fuentes (tasa, ubicación, histórico) y propiedades de los contaminantes

 Niveles piezometricos/concentraciones iniciales

MODELO NUMÉRICO
UTILIDAD
53
 MEJORAR EL CONOCIMIENTO hidrogeológico (3-D!)
 AUXILIAR EN EL MANEJO de las aguas subterráneas
 Determinar caudales de bombeo sostenibles
 Calcular zonas de protección de aguas subterráneas y pozos de captación
 Estimar caudales de drenaje de minas/instañaciones (tajo, túneles, bocaminas, drenes, etc)
CIO
 PREDECIR IMPACTOS hidrológicos

 Impacto del drenaje de la mina en un rio, humedal, campos agrícolas vecinos a la mina
 Influencia de una represa en el flujo regional del acuífero
 ¿Qué pasa si los sistemas de drenaje fallan?
 ¿Hasta dónde una pluma de contaminación migraría sin medidas de remediación?
 DISEÑAR Y OPTIMIZAR sistemas de bombeo/remediación
 Drenaje de minas, control de presión de poros.
 Control de plumas de contaminación por bombeo, remediación in situ, etc.
UNI - FIC - DAHH
 ANALIZAR LOS RIESGOS AMBIENTALES

 Determinar direcciones mas probables de transporte de contaminantes.
 Predecir concentraciones de contaminantes en un cuerpo receptor para estimar el riesgo a la
salud humana.
53
MODELOS NUMÉRICOS
DIFERENCIAS FINITAS Y ELEMENTOS FINITOS
54
Diferencias Finitas Elementos Finitos Diferencias Finitas con

Volúmen de Control
Gomes et al (2007) FEFLOW Model EIA, KP and Co., 2010 Modflow USG, AQUAVEO’s images (GMS)
54
MODELOS NUMÉRICOS
COMPARACIÓN
55
Método Ventajas Limitaciones
Diferencias Finitas ▪ Conserva la masa ▪ Red Retangular resulta en
▪ Intuitivo, simple, conocido y comprobado, fácil de discretización innecesaria fuera de
programar la zona de interes
▪ Más programas disponibles ▪ Dificil adaptación a sistema muy
CIO
▪ Hay versión open source (Modflow) complejos con capas deformes

▪ Pozos extraen agua de toda la celda
▪ Capas, lineas y conulas continuas
▪ Tensor de K sin derivadas cruzadas
Elementos Finitos ▪ Discretización mas eficiente en zonas de interes ▪ No tan intuitivo, teoría más
solamente compleja
▪ Pozos pueden ser representados como nudos ▪ Conserva masa solo globalmente.
discretos Errores locales en el balance de
▪ Bien adaptado a geología compleja aguas pueden ocurrir (raros en
▪ Tensor de K totalmente representado FEFLOW)
▪ Capas no tienen que ser continuas ▪ Solución mas compleja,
discretización mas elaborada,
UNI - FIC - DAHH
requiere software especializados

Diferencias Finitas ▪ Discretización eficiente ▪ Poco uso en aguas subterráneas
▪ Bien adaptado a geología compleja (adopción?)
con Volumen de ▪ Tensor de K totalmente representado ▪ Solución mas compleja,
▪ Capas no tienen que ser continuas discretización mas elaborada,
Control requiere software especializados
55
MODELOS NUMÉRICOS
MATRICES
56
La mayoría de modelos numéricos se reducen a ecuaciones lineales que
despejan h o h y q.
Ecuación matricial del MODFLOW para un modelo con 3 lineas, 4 columnas y 2 capas! Fuente:
Manual del MODFLOW (McDonald and Harbaugh, 1988) 56
MODELO NUMÉRICO
RESULTADOS
57
 Niveles freáticos (hidroisohipsas) y
piezométricos, isoprofundidades
(equidistancia de la superficie a la
napa).
CIO
 Concentraciones en cada celda.

 Redes de flujo en secciones y en 3D.
 Caudales y velocidades.
 Balance de agua y de contaminantes.
 Distribución de presiones (para diseño
de taludes y túneles).
 Series temporales o hidrógrafas y
UNI - FIC - DAHH
curvas de paso de contaminantes en

cualquier punto, superficie o volumen.
Hidroisohipsas en Acuífero de Lima

57
MODELO NUMÉRICO
RESULTADOS
58
 Trayectoria de partículas lanzadas en
el campo de flujo
58
MODELO CONCEPTUAL
EJEMPLO
59 Condiciones de borde
TIPO I : Niveles del río, océano (0 m), divisoria (60,5 m),
sección del río (101,5 m), pozos y piezómetros. DIVISORIA CON
TIPO II: Divisoria (Q=0), roca madre (Q=0), ACUIFERO VECINO
Recarga = Ineficiencia de riego = RIEGO – ETR
Recarga = 265 mm/año, efecto Qpozos incluido.
CIO
RIEGO AFLUENTE POR EL

ETReal = ETPotencial LECHO ALUVIAL
( lo calculará el modelo
en base a sección, K y
cota nivel freático)
DESCARGA
AL MAR Llueve poquísimo, no
hay lluvia efectiva
UNI - FIC - DAHH
RIO ZAÑA
Pozos
59
Fuente: Mendoza, M., 2009
MODELO CONCEPTUAL
EJEMPLO
60
ESTRATO ACUÍFERO
CIO
CAPA 1
ESTRATO ARCILLOSO
CAPA 2
Dz (m) BORDE IMPERMEABLE
R (Ohm.m)
Profundidad (m)
UNI - FIC - DAHH
BASAMENTO ROCOSO
IMPERMEABLE
SECCIÓN TRANSVERSAL DEL VALLE CON BASE EN INVESTIGACIONES GEOFÍSICAS

MODELO CONCEPTUAL
EJEMPLO
PARÁMETROS DE PRUEBAS DE POZO
61 Trasmisividad
Espesor (m2/sec) Cond.
Coef. de
Prueba Pozo Saturado hidráulica
Almacenamietno S
(m) Descenso Recuperación K (cm/s)
1 15_029 26.59 1.10E-02 4.14E-02

2 15_031 12.70 8.30E-02 4.40E-02 4.76E-01
3 16_007 24.25 3.00E-02 4.60E-02 1.53E-01
4 16_017 14.72 4.00E-02 1.90E-02 1.87E-01
CIO
5 16_032 24.70 1.20E-02 1.50E-02 5.43E-02

16_051 31.16 1.20E-02 3.00E-02 6.09E-02
6
16_052 (P) 31.00 7.00E-02 6.00E-02 2.09E-01
16_065 (P) 34.70 1.50E-02 4.32E-02
7
16_066 34.70 4.50E-02 2.00E-02 8.65E-02
16_075 32.20 6.80E-02 2.11E-01 0.0120
8
16_068 (P) 22.45 2.90E-02 1.29E-01
9 16_094 28.30 2.00E-02 2.20E-02 7.41E-02 0.0900
10 07_098 30.00 2.60E-02 8.67E-02 0.0110
11 16_092 32.75 2.90E-03 8.85E-03 0.0002
07_095 30.64 1.83E-02 5.96E-02 0.1259
12
07_096 (P) 25.62 5.68E-02 2.22E-01 0.0003
13 07_132 49.59 3.02E-03 6.10E-03 0.1294
07_138 26.79 1.57E-02 5.87E-02 0.0014
14
07_137 (P) 8.47 2.59E-02 3.06E-01 0.0040
UNI - FIC - DAHH
CONDUCTIVIDAD HIDRAULICA INICIAL DE LAS CAPAS (el modelo numérico luego las calibra)
Capa Kx (cm/s) Ky (cm/s) Kz (cm/s)

1 8.35E-02 8.35E-02 8.35E-03
2 5.00E-08 5.00E-08 5.00E-08
MODELO NUMÉRICO
EJEMPLO – MODELO EN DIFERENCIAL FINITAS MODFLOW-SURFACT
62
62
MODELO NUMÉRICO
EJEMPLO
CONTROL DE RESULTADOS
63
Cargas Hidráulicas Calculadas vs. Frecuencia de Residuales

Observadas
100 30
CIO
Frecuencia
25
Cargas Hidraúlicas Calculadas (msnm)
90
80 20
70 15
y = 1.0109x
R² = 0.9974 10
60
50 5
40
0
-8 -6 -4 -2 0 2 4 6 y8 mayor...
30
Residuales (m)
20
10
UNI - FIC - DAHH
0
0 20 40 60 80 100
Cargas Hidraúlicas Observadas (msnm)

Tiene controles similares a los de los
modelos precipitación escorrentía

MODELO NUMÉRICO
EJEMPLO
RESULTADOS
64
Khx Khy Kvz
La permeabilidad se Zona (cm/seg) (cm/seg) (cm/seg) Kh/Kv
distribuye en 6 zonas. 1
2
9.73E-02
7.03E-05
9.73E-02
7.03E-05
4.63E-03
6.88E-06
21
10.2
3 6.90E-02 6.90E-02 3.47E-03 19.9
4 1.90E-03 1.90E-03 8.10E-05 23.4
CIO
5 1.30E-01 1.30E-01 4.63E-03 28.1

6 7.08E-03 7.08E-03 2.38E-03 3
5
4
1 3 2
UNI - FIC - DAHH
CAPA 1 CAPA 2 Fuente: Mendoza, M., 2009

MODELO NUMÉRICO
EJEMPLO
RESULTADOS
NAPA FREÁTICA Y DIRECCIÓN DEL FLUJO SUBTERRÁNEO
65
SECCIONES Y LÍNEAS EQUIPOTENCIALES

MODELO NUMÉRICO
EJEMPLO
RESULTADOS
APLICACIÓN PREDICTIVA
66
Hipótesis:
Caudal de explotación de 23 L/s en los pozos
del distrito de Cayalti y de 5.8 L/s en los distritos
Restantes (138.4 L/s en total)
CIO
CONCLUSIONES
La margen izquierda del río Zaña (parte baja) presenta

excelentes condiciones para el aprovechamiento.
UNI - FIC - DAHH
Las altas tasas de recarga y la alta conductividad

hidráulica de la primera capa hacen que los impactos
negativos en la napa freática sean mínimos.

67
PARACAS
GRACIAS

U2.8b Aguas Subterraneas RevA PDF

Cargado por

Copyright:

Formatos disponibles

U2.8b Aguas Subterraneas RevA PDF

Cargado por

Información del documento

Título original

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

U2.8b Aguas Subterraneas RevA PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERIA CIVIL

RECURSOS HIDRAULICOS HH333

IDENTIFICACIÓN DEL ÁREA DE RECARGA CON HIDROLOGÍA ISOTÓPICA

Luego fluye a lo largo de un gradiente

Por lo tanto, el agua subterránea constituye agua superficial prestada.

Es solamente una cuestión de tiempo antes de que la recarga natural se convierta en

Es importante evaluar el impacto en el flujo

Pero la relación entre el lugar de recarga y el

local, (2) subregional o (3) regional.

Estudios geológicos y de calidad de agua ayudan a inferir el origen de las agua

Fuente: H. Jegat, 2015

En el agua de recorrido local predominan los carbonatos, luego progresivamente

fases del ciclo hidrológico.

 Protio ( 1 H, H ) 99,98%, estable

 Oxigeno 18 ( 18 O ) 0,20%, estable

Isótopos Estables del Agua

El contenido isotópico d18O y dD

una estrecha franja lineal

Fuente: Soler, A, 2015

Fuente: Soler, A, 2015

Fuente: Jegat, H.,2015

acuífero del Valle

Fuente: Jegat, H.,2015

Oct 2002: 58 muestras

Asociación del agua de los

DETERIORO DE LAS AGUAS

La zona afectada prácticamente no se recupera ni puede explotarse agua dulce.

Principalmente los contaminantes se

movilizan según el mecanismo de

 PROCESOS QUÍMICO-TÉRMICOS Propiedades

FUENTE: GUIGUER, N., 2017.

rb = densidad del suelo

Advección Dispersión Sumidero/Fuente Reacciones Retardo  = porosidad

 La Advección se debe al flujo

FUENTE: GUIGUER, N., 2017.

FUENTE: GUIGUER, N., 2017.

 Difusión Molecular (D*)

FUENTE: GUIGUER, N., 2017.

FUENTE: GUIGUER, N., 2017.

Nota: ADSORCIÓN es retención de solutos en superficies.

FUENTE: GUIGUER, N., 2017.

EXPLOTACIÓN CONJUNTA AGUAS SUPERFICIALES Y

• Diques y empozamientos aumentan el área de recarga y favorecen la

extensión de la infiltración y de las tasas de infiltración.

Se alterna bombeo de recarga al acuífero, con bombeo

Para satisfacer mejor a la demanda en cantidad y calidad se combinan funciones de captación,

almacenamiento y depuración. Las aguas subterráneas afianzan la garantía por su propiedad de

La vulnerabilidad específica incluye parcialmente el concepto de riesgo, toda vez que

de hurto o robo es bajo, porque en la región no hay vándalos.

Incluyendo factores de reducción y corrección, en función de la humedad,

Como toda ecuación diferencial, deben especificarse las condiciones

 Condición sobre la variable h.

 Flujo de agua definido que sale

 Las divisorias de agua también se

 Condición relacionada con la

 Se aplica a las entradas y salidas de

conductividad hidráulica vertical del Flujo condicionado a través AB

Modelo Conceptual Modelo Matemático