5/18/21

Fundamentos de
Ingeniería Ambiental

Diana Ita Nagy

18 de mayo 2021

Aguas subterráneas

Características generales, Ley de Darcy,

contaminación, recuperación y conos
de depresión y curvas de captura

1
 5/18/21

Acuífero y Manantial

Acuífero Manantial
• Capa geológica saturada • Resulta cuando un
suficientemente acuífero se llena hasta el
permeable para permitir punto en que el agua se
el flujo de agua a su través desborda a la superficie
de la tierra.
• Se asienta sobre una base • Se encuentran
o lecho de confinamiento principalmente en las
• Estos dos términos no calizas y dolomitas
están definidos de forma • Mecanismo de disolución
precisa y tienen a menudo y fractura forman espacios
un significado relativo para el flujo de agua

Nivel Freático vs Nivel Piezométrico

En ambos casos corresponde al nivel superior de un
acuífero
Nivel Freático Nivel piezométrico
• La presión del agua es • La presión del agua es
igual a la presión mayor que la atmosférica
atmosférica • Puede alcanzar una altura
• Acuíferos libres mayor al nivel en que se
encuentra confinada
• Acuíferos confinados

*Nivel hidrostático: Es la distancia vertical entre un nivel de referencia

(normalmente el nivel del mar) y el nivel freático.

2
 5/18/21

Subsidencia del suelo por bombeo de

agua subterránea
• Disminución de la elevación de la superficie terrestre
por cambios subterráneos.
• Causas:
o Bombeo de agua, petróleo y gas de las reservas
subterráneas.
o Disolución de acuíferos de caliza.
o Colapso de minas subterráneas.
o Drenaje de los suelos orgánicos.
o Humedecimiento de suelos secos.

3
 5/18/21

Ganancia vs. Pérdida de Corriente

Los intercambios de agua entre los depósitos de agua
superficial y los acuíferos son importantes
Ganancia de corriente Pérdida de corriente
• Los ríos usualmente • Es posible que las
empiezan como pequeños corrientes entreguen agua
arroyos y aumentan en al subsuelo en algunos
caudal a medida que puntos
fluyen hacia el mar • Es la “recarga” de agua al
• El agua que ellos ganan acuífero
frecuentemente viene del
agua subterránea

Movimientos del Agua Subterránea

Fuente: [Link]

4
 5/18/21

Porosidad

La cantidad de agua que puede ser almacenada

en un acuífero depende de la porosidad de la
roca o del suelo:

!"#$%&' (& )*+í"-

Porosidad (η)= !"#$%&' .".*# (-ó#1("- 2 )*+í"-)

Porosidad Efectiva
• No toda el agua que se encuentra en el acuífero
puede ser aprovechada.

• El volumen aprovechable (recuperable) en un

acuífero no confinado se llama rendimiento
específico o porosidad efectiva.

• En acuíferos confinados, se usa el coeficiente de

almacenamiento, ya que el rendimiento está
afectado por la presión desarrollada al remover el
agua.

10

5
 5/18/21

Valores de Porosidad y Rendimiento

Específico

Rendimiento de un acuífero (m3) = Rendimiento específico (%) × Volumen total del acuífero

Material Porosidad (%) Rendimiento específico

(%)
Arcilla 45 3
Arena 34 25
Grava 25 22
Gravilla y arena 20 16
Arenisca 15 8
Caliza, pizarra 5 2
Cuarcita, granito 1 0.5

11

Ejemplo

Para un acuífero de arena, ¿qué volumen de

agua puede ser almacenado en una columna
saturada si tiene una sección transversal con
área de 1 m2 y profundidad de 2 metros?
¿Cuánta agua podría ser extraída?
1. Volumen del material: 1 m2 x 2 m = 2 m3
2. Volumen de agua: Porosidad x Vol de material
= 0.34 x 2m3 = 0.68 m3
3. Cantidad que puede extraerse:
Rendimiento: Rendimiento específico x Volumen de material
= 0.25 x 2m3 = 0.5 m3

12

6
 5/18/21

Nivel hidrostático y líneas equipotenciales

• Nivel hidrostático (h):

o Es la distancia vertical entre un nivel de referencia
(normalmente el nivel del mar) y el nivel freático.
Dimensión: metros de agua
• Líneas equipotenciales:
o Líneas que unen los puntos de misma altura
hidrostática.
• Líneas de corriente:
o Dirección del flujo de agua subterránea
(perpendiculares a las líneas equipotenciales).

13

Gradiente hidráulico

• Acuífero no confinado: pendiente del nivel

freático medida en dirección al nivel máximo.
• El agua en el acuífero se va a mover en la
dirección del mayor gradiente hidráulico.

!"#$%& '( )%*(+ ,%'-&./á/%!&

Gradiente hidráulico =
'%./")!%" ,&-%1&)/"+

,!2," ',
Gradiente hidráulico = 3
= '3

14

7
 5/18/21

Gradiente hidráulico
El agua en el acuífero se va a mover en la
dirección del mayor gradiente hidráulico.

Nivel hidrostático y gradiente en un acuífero libre. El flujo de agua subterránea va de

izquierda a derecha en el plano de la página, y el gradiente es Δh/L.
15

Gradiente hidráulico
Red de flujo bidimensional:
Consistente en líneas de corriente y líneas equipotenciales. Si
dos pozos coinciden en la misma línea de corriente, el
gradiente es justamente Δh/L.

16

8
 5/18/21

¿Cómo calcular el gradiente

hidráulico con tres pozos?

17

¿Cómo calcular el gradiente

hidráulico con tres pozos?

La altura de columna de agua (h) en

cualquier punto de esta línea es constante.
Flujo de agua es perpendicular a esta línea.

18

9
 5/18/21

Ejemplo
¿Gradiente hidráulico?
• 2 pozos a 200m de distancia este-oeste.
• Pozo oeste: h = 30.2 m
• Pozo este: h = 30.0 m
• Tercer pozo 100m al sur del pozo este: h = 30.1 m
200 m 30.1
30.2 30.0
La distancia entre la equipotencial y
el pozo este se determina
L geométricamente:
100 m
L = 100/√2.
Línea equipotencial
!".$%!"." &
30.1 Gradiente = !"" = 0.00141
&
#

19

Ley de Darcy
!"
𝑄 =𝐾∗𝐴∗
!#

Donde:
Q= caudal (m3/día)
K= conductividad hidráulica o coeficiente de
permeabilidad (m/día)
A= área de sección transversal (m2)
#$
= gradiente hidráulico
#%

Se asume:
Velocidad = constante (en toda la sección de análisis)
Flujo laminar

20

10
 5/18/21

La Ley de Darcy
• Es la ecuación básica que describe el flujo de agua
subterránea en un lecho poroso
• Limitaciones
o No se cumple en zonas turbulentas (cerca al pozo de extracción)
o No válido en arcillas coloidales (materiales finos)
o Solo en medios saturados
o Fuerzas inerciales despreciables
o Dependiente de la temperatura
• La constante de proporcionalidad 𝐾 puede variar para un
mismo suelo dependiendo de la diferente ordenación de las
partículas (anisotropía) o forma. También por
concentraciones de limo o arcilla. Es aproximada

21

Ley de Darcy
Caracterización de acuíferos, según su conductividad hidráulica

• Acuíferos homogéneos à misma conductividad hidráulica (K)

en todos sus puntos.

• Acuíferos heterogéneos à diferente conductividad (K)

• Acuíferos isotrópicos à misma conductividad hidráulica (K)

en cualquier dirección.

• Acuíferos anisotrópicos à conductividad depende de la

dirección. Común: K más alta en dirección horizontal que
vertical.

22

11
 5/18/21

Velocidad de Darcy 𝑣
𝑑ℎ
𝑄 = 𝐾𝐴
𝑑𝐿
𝑄 𝑑ℎ
𝑣= =𝐾
𝐴 𝑑𝐿
- Q: Caudal, m3/día
- K: Conductividad hidráulica, m/día
- A: Sección de flujo, m2
- v: Velocidad de Darcy, m/día
- dh/dL: Gradiente hidráulico

No es una velocidad «real». Asume que toda la sección A es capaz de dejar fluir el agua a su
través. Puesto que gran parte de la sección está formada por sólidos, el área real a
través de la cual fluye el agua es mucho menor y como resultado la velocidad real es mayor.

23

Velocidad Real promedio 𝑣′

(de filtrado)
Como el volumen de control no está
totalmente vacío:
𝑄 = 𝐴 A 𝑣 = 𝐴′ A 𝑣′
Q
𝐴
𝑣' = 𝑣
Partículas 𝐴′
sólidas
𝐴A𝐿 𝑉
-L- 𝑣' = 𝑣= 𝑣
𝐴′ A 𝐿 𝑉′
𝑉′
Sea: 𝜂=
𝑉
• 𝐴: Sección transversal total
𝑣
• 𝐴’: Sección transversal “vacía” 𝑣' =
𝜂
• 𝑉 = 𝐿𝐴: Volumen total
• 𝑉′ = 𝐿𝐴′: Volumen “vacío” 𝑣 𝐾 𝑑ℎ
𝑣' = =
• 𝜂 = 𝑉 ' /𝑉: porosidad 𝜂 𝜂 𝑑𝐿

24

12
 5/18/21

Conductividad hidráulica

Material Conductividad (m/día)

Arcilla 0.0004
Arena 41
Grava 4100
Grava y arena 410
Arenisca 4.1
Caliza 0.041
Cuarcita 0.0004
Linsley et al., 1992

25

Ejemplo
Se tiene un acuífero confinado de 20 metros de espesor con dos pozos
de monitoreo espaciados 500 metros en la dirección del flujo. La
diferencia del nivel freático del pozo es de 2,0 metros. La
conductividad hidráulica es de 50 metros/día.

Estimar el tiempo (t) que necesita el agua del primer al segundo pozo,
la velocidad de Darcy y la velocidad promedio lineal. La porosidad del
acuífero es de 35%.

()
Velocidad de Darcy ʋ = K = 50 m/día × 0.004 = 0.20 m/día
(*

ʋ´ = (0.20 m/día)/0.35 = 0.57 m/día

t = 500m/0.57m/d = 877 d

26

13
 5/18/21

Contaminación de aguas
subterráneas

27

Contaminación del agua subterránea

Las vías por las que distintas sustancias llegan a los acuíferos
contaminando las aguas subterráneas son muy diversas:
• Infiltración de sustancias depositadas en superficie o de la lluvia a
través de ellas.
• Filtración de sustancias almacenadas bajo tierra, o disolución de
ellas por el agua subterránea.
• Filtración desde un río influente.
• Derrames accidentales de depósitos superficiales o enterrados.
• Desde la superficie, a través de captaciones abandonadas o mal
construidas.
• Inyecciones en pozos.

Una vez contaminada => muy difícil de restaurar

28

14
 5/18/21

Contaminación del agua subterránea

29

Contaminación del agua subterránea

30

15
 5/18/21

Contaminación del agua subterránea

Líquidos de fase no acuosa (NAPL) = no se disuelven bien en el
agua.
àRecuperación del acuífero difícil o imposible
2 Grupos:
• Más densos que el agua (DNAPL): Tricloroetileno (TCE),
tetracloroetileno, policlorobifenilos (PCB), algunos
pesticidas, hidrocarburso aromaticos policiclicos (PAH),…
• Menos densos que el agua (LNAPL): Hidrocarburos como
petróleo, gasolina, benceno,….

31

Contaminación del agua subterránea

32

16
 5/18/21

Transporte de contaminantes en aguas

subterráneas
En el agua subterránea los contaminantes se
pueden mover por:

• Advección.
• Dispersión.
• Difusión.

33

Transporte de contaminantes
Advección
Los solutos son
transportados por el
movimiento de la masa
de agua.

Movimiento en la misma
dirección y velocidad
que el flujo de agua.

à Velocidad de flujo real

34

17
 5/18/21

Transporte de contaminantes
Difusión
El transporte difusivo es debido al movimiento de iones disueltos y moléculas debido
a la existencia de gradientes de concentración (de mayor a menor concentración).

Movimiento aleatorio tanto del soluto como del solvente

En ausencia de movimiento de agua subterránea, el contaminante se difundirá en
todas direcciones.

Movimiento lento (comparado con la advección)

Un aporte instantáneo (pulso) en un

flujo crea una pluma que se va
difundiendo según se mueve
gradiente abajo:
(a) en una dimensión
(b) en dos dimensiones

35

Transporte de contaminantes
Dispersión
Los diferentes tipos, tamaños y orientación de los poros dan
lugar a velocidades que difieren de unos puntos a otros.

Los solutos tienden también a moverse con diferente

velocidad. Esta desviación de la trayectoria ideal se denomina
dispersión hidráulica.

36

18
 5/18/21

Transporte de contaminantes
Difusión y Dispersión
Cuando la dispersión se produce debido a la difusión molecular y a la
no homogeneidad de velocidades microscópicas, se denomina
dispersión hidrodinámica (o mecánica).
Dispersión
• Caminos irregulares
• El agua podría pasar más rápido o
lento

+
• Trayecto irregular

Difusión molecular
• Movimiento del contaminante en
el agua
• Debido a una diferencia de

=
concentraciones

Dispersión hidrodinámica
• Ambos vuelven difusos los bordes
de la pluma
• Combinación de los dos
mecanismos

37

Transporte de contaminantes
Sorción
A medida que los contaminantes se movilizan por el
acuífero, estos son adsorbidos por los sólidos o
absorbidos:

• Absorción à el contaminante pasa a formar

parte de los sólidos, hay una transformación.
• Adsorción à el contaminante se adhiere a la
superficie de los sólidos.

38

19
 5/18/21

Transporte de contaminantes
Retardo
• La razón del total de contaminante en una unidad de
volumen que es disuelto en el agua subterránea es
llamado factor de retardo.

• El factor de retardo de los iones de cloruro es de 1 (todos

son disueltos en el agua).
• Un factor de retardo de 5:
o 20% del contaminante es disuelto en el agua
o 80% sorbido por los sólidos del acuífero.

Velocidad del contaminante = velocidad lineal/factor de retardo

39

Transporte de contaminantes
Retardo

• En estado estacionario, es la razón de la velocidad

promedio del agua subterránea 𝑣’ a la velocidad
del material sorbido, 𝑣𝑠:
𝑣′
𝑅=
𝑣$
• Por ejemplo, un factor de retardo de 5 indica que
el agua subterránea viaja cinco veces más rápido
en un periodo determinado de lo que lo hace el
contaminante

40

20
 5/18/21

Transporte de contaminantes
Retardo

Contaminantes por un acuífero

• Unos son absorbidos (entran) a

las partículas a lo largo de su
camino

• Otros son adsorbidos (se adhieren

a la superficie de las partículas)

41

Transporte de contaminantes
Ejemplo
Consideramos 3 pozos de observación, cada uno situado
en un vértice de un triangulo equilátero. La distancia
entre cualquier par de pozos es de 300 m. Si el borde
frontal de una pluma es perfectamente recto y la
velocidad lineal es de 5 m/día, señale cuanto tiempo
tardaría en viajar del pozo 1 después de haber llegado
(simultáneamente) a los pozos 2 y 3. Asuma un factor de
retardo de 2.
L = (300sin60°) m

t = (300sin60°) m / (5.00 m/día) = 51.96 d

Con factor de retardo:

tcontaminante = 51,96 días * 2 = 103.92 días

42

21
 5/18/21

Conos de depresión y curvas de

zona de captura

43

Cono de Depresión
• Para proveer al pozo de agua debe haber un
gradiente hacia el pozo para proveerlo de agua
• Es más, a mayor rapidez de bombeo más
pronunciado es el gradiente en la vecindad del
pozo
• Al bombear agua de un pozo, el nivel freático de
un acuífero libre, o el piezométrico en uno
cautivo, forma un cono de depresión, en las
inmediaciones de dicho pozo

44

22
 5/18/21

Cono de Depresión (Acuífero Libre)

Tasa de bombeo
Ley de Darcy para volumen de control cilíndrico
(acuífero libre):
𝑑ℎ$ 𝑑ℎ$
𝑄 = 𝐾𝐴 = 𝐾2𝜋𝑟ℎ$
𝑑𝑟 𝑑𝑟
Integrando:
%" 𝑑𝑟 &#,"
/ 𝑄 = / 𝐾2𝜋𝑟ℎ𝑑ℎ$
%! 𝑟 &#,!
𝑟' ( − ℎ(
𝑄𝑙𝑛 = 𝜋𝐾 ℎ$,' $,(
𝑟(
( − ℎ(
𝜋𝐾 ℎ$,' $,(
𝑄=
𝑙𝑛 𝑟'⁄𝑟(

• K = Conductividad hidráulica.
• (dhc/dr) = gradiente hidráulico de un punto del borde
del cono de depresión al pozo.
• hc,1 y hc,2 = Profundidad del acuífero para puntos
dentro del cono de depresión (pozos de
observación).
• r = Distancia de los pozos de observación del pozo
*Relación entre la tasa de bombeo y el descenso del principal.
nivel freático en el área del cono de depresión.
45

Cono de Depresión (Acuífero Confinado)

Procedimiento similar. Altura del cilindro de radio r, a través del cual fluye el agua para
alcanzar el pozo, es una constante igual al espesor del acuífero, B.
El cono de depresión aparece ahora en la superficie piezométrica.

𝑑ℎ+
𝑄 = 𝐾2𝜋𝑟𝐵
𝑑𝑟

𝐾2𝜋𝐵 ℎ+,$ − ℎ+,-

𝑄=
𝑙𝑛 𝑟$⁄𝑟-

hc = Profundidad del cilindro = la distancia entre el límite inferior del

acuífero y la superficie piezométrica.

46

23
 5/18/21

Curvas de Zona de Captura

• Como un pozo de extracción bombea agua, la superficie
hidrostática del acuífero (la superficie piezométrica en caso
de que sea cautivo) en su vecindad experimenta una
depresión, por lo que crea un gradiente hidráulico que
arrastra la pluma hacia el pozo
• Con pozos adecuadamente situados, las aguas subterráneas
contaminadas pueden ser extraídas del acuífero para ser
descontaminadas y utilizadas o devueltas al acuífero
• La zona de captura debe ser del mismo ancho o mayor que
el penacho para que el pozo pueda limpiar eficientemente
el agua contaminada.

47

Curvas de Zona de Captura

La zona de captura es la región en la cual todas las
líneas de flujo convergen en el pozo de extracción.

48

24
 5/18/21

Curvas de Zona de Captura

Javandel y Tsang (1986)
obtienen una relación entre las
coordenadas x e y de la zona de
captura para un pozo único de
extracción en el origen de
coordenadas.

𝑄 𝑄 𝑦
𝑦=± − 𝑡𝑎𝑛%$
2𝐵𝑣 2𝜋𝐵𝑣 𝑥
Donde:
𝑦 B = espesor del acuífero (m).
𝑡𝑎𝑛𝜙 =
𝑥 ʋ = Velocidad de Darcy (m/día).
Q = Tasa de bombeo (m3/día).
𝑄 𝜙
𝑦= 1−
2𝐵𝑣 𝜋

49

Curvas de Zona de Captura

Superposición de la pluma en una curva de zona de captura para un solo pozo
de extracción.

50

25
 5/18/21

Curvas de Zona de Captura

𝑄 𝜙
𝑦= 1−
2𝐵𝑣 𝜋

Esta ecuación nos y = Q/2Bʋ

ayuda predecir datos y = Q/4Bʋ

importantes de la zona
de captura:

En caso que: x = infinito

Entonces: ɸ = 0 e y = Q/2Bʋ
à anchura total de la zona de captura = Q/Bʋ

En caso que: ɸ = π/2

Entonces: x = 0 e y = Q/4Bʋ
à anchura de la zona de captura medida sobre el eje x = Q/2Bʋ

51

Curvas de Zona de Captura

Ejemplo
Supongamos un derrame de 0.10 m3 de TCE se distribuye
homogéneamente por un acuífero de 10m de espesor,
formando una pluma rectangular de 250 m de ancho. La
velocidad de Darcy es de 10-6 m/s. ¿Cuál sería a tasa de
bombeo de un pozo situado sobre el borde frontal de la
pluma, para capturar toda la pluma?

Q = 2Bʋ × 250

Q = 2 × 10m × 10-6 m/s × 250 m

=250m
Q = 0.005 m3/s

52

26
 5/18/21

Curvas de Zona de Captura

• Un pozo no siempre es la mejor solución.
• Si el único pozo está muy lejos à mejor varios pozos más
cerca a la cabeza de la pluma.
• Bombeando al mismo ritmo, alineados a lo largo del eje y.

53

Curvas de Zona de Captura

Una ecuación general para la mitad positiva de la
curva de captura para n pozos óptimamente
espaciado, situados regularmente a lo largo del eje
y, es:
(
𝑄 1
𝑦= 𝑛 − 1 𝜙%
2𝐵𝑣 𝜋
%&'

ɸi es el ángulo que forma la recta que une el pozo i-

ésimo con un punto de la curva de captura, con la
horizontal.

54

27