Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
MODELACIÓN AMBIENTAL
Hidrología e Hidráulica de Aguas
Subterráneas
Luis A. Camacho B.
Relaciones suelo - agua
1. Sólida
Fases: 2. Líquida
3. Gaseosa
Volumen de control: volumen sólidos + volumen agua +
volumen aire y vapor de agua
Vg Aire Wg 0
Vv
Vw Agua Ww
V W
Vs Suelo Ws
Contenido volumétrico de humedad: Vw V
Volumen de vacíos: Vv V Vs Vw Vg
Porosidad del suelo: n Vv V (n : 0.12 0.55)
Relación de vacíos: e Vv V s
Grado de saturación: S Vw Vv
Vg Aire Wg 0
Vv
Vw Agua Ww
V W
Vs Suelo Ws
Vw V Vv V Vs Vw Vg n Vv V S Vw Vv e Vv V s
nS
Contenido de humedad en peso: w Ww Ws
Peso unitario (densidad): Gravedades específicas:
Total: T W V Masa: Gm T 0
Sólidos: s Ws Vs Sólidos: Gw w 0
Agua: w Ww Vw Agua: G s 0
Seco: seco Ws V 0 densidad de agua a 4C
Vg Aire Wg 0
Vv
Vw Agua Ww
V W
Vs Suelo Ws
Ejemplo 1
Un núcleo cilíndrico de 30 cm de largo y 12 cm de diámetro
tomado de un suelo de masa 6691 g requirió 509 cm3 de agua para
llevarlo a saturación. El núcleo fue secado después en un horno a
105 °C por 24 horas, y después fue pesado nuevamente y su
masa determinada en 6114 g. Calcule la porosidad de la muestra
de suelo, su relación de vacíos y el contenido inicial de humedad.
FACTORES QUE AFECTAN EL MOVIMIENTO DEL
AGUA
1. Textura: distribución de tamaños de partículas del suelo
2. Estructura: ordenamiento y agregación de las partículas
definen el tamaño de poros
1. Poros capilares
2. Poros no capilares
3. Clase de minerales de arcilla
4. Cantidad de materia orgánica
Triángulo Arcilla: < 0.002 mm
Textural Limo: 0.002 – 0.05 mm
Arena: 0.05 – 2 mm
Arena muy fina: 0.05 – 0.1 mm
Arcilloso Arena fina: 0.1 – 0.25 mm
Arena media: 0.25 – 0.5 mm
Arcillo
Arcillo limoso Arena gruesa: 0.5 – 1 mm
arenoso
Franco Franco
Arena muy gruesa: 1 – 2 mm
Franco arcillo- arcilloso arcillo-limoso
arenoso
Franco Franco
Franco limoso
arenoso Limo
Arenoso
Arena (%)
LEY DE DARCY
flujo Saturado
h
Q
Q A
K
x
h
Q= K A
x
h=h2- h1
hp1
hp2
h1
h2
z1 x=x2- x1
h
Q= K x A z2
Datum
z=0
LEY DE DARCY (1856)
Velocidad de Darcy o descarga específica
q Ki Gradiente hidráulico
Conductividad hidráulica
Q qAt
Area total
Q KiAt
Porosidad
v Q / Av q At / Av q / n
Area de vacíos
Velocidad promedio del flujo en el medio
poroso
Ley de Darcy
Ejemplo 2
1. En una cierta área el agua subterránea descarga a un canal. El
suelo tiene una conductividad hidráulica de K = 1.0 cm/s. y
una porosidad de 0.2. El flujo de agua subterráneo es
prácticamente horizontal y el gradiente hidráulico, de
magnitud 0.01, va en una dirección, vista en planta, de 45°
con respecto al alineamiento del canal. Se introduce un
trazador conservativo en el agua subterránea a una distancia
perpendicular al canal de 6m. Si la dispersión y difusión se
asumen despreciables, estime el tiempo que demora el
trazador en aparecer en el canal.
LEY DE DARCY
Aplicabilidad: solamente para condiciones de flujo laminar
qd10 qd10
Diámetro tamiz
que retiene 90%
R de material o pasa 10%
Número de Reynolds
Viscosidad cinemática
Q
q Descarga específica
At
Estrictamente aplicable:R 1
Aproximadamente aplicable: 1 R 10
No válida: R 10 (flujo turbulento)
Heterogeneidad y anisotropía
K = f (propiedades del fluido: peso específico y viscosidad
propiedades del medio poroso: textura, distribución de
granos, etc.)
Anisotropía:
Sea Kx = conductividad hidráulica para el flujo horizontal
Kz = conductividad hidráulica en la dirección vertical
Si Kx = Kz acuífero isotrópico en la conductividad
hidráulica
Lo contrario anisotrópico Los más frecuentes en la naturaleza
Heterogeneidad y anisotropía
La mayoría de formaciones geológicas no son uniformes en
composición y las propiedades del medio poroso varían
espacialmente y con la dirección. Considere la conductividad
hidráulica K:
Homogénea: independiente de la posición en la formación geológica
K(x, y, z) = c (contanste)
Heteregénea: dependiente de la posición K ( x, y, z ) c
Tipos de heterogeneidad: - Estratos geológicos por
capas
- Discontínua (fallas)
Isotrópica: Independiente de la dirección en la formación geológica
Kx Ky Kz
Anisotrópica: Dependiente de la dirección en la formación geológica
Kx Ky Kz
Heterogeneidad y anisotropía
Isotropía transversal: donde la conductividad hidráulica tiene las
componentes horizontal y vertical iguales
Kx Ky Kz
Homogéneo/isotrópico Homogéneo/anisotrópico
Heterogéneo/isotrópico Heterogéneo/anisotrópico
Realidad: Modelación:
Heterogeneidad y anisotropía acuífero homogéneo equivalente
ó acuífero isotrópico equivalente
Usualmente acuíferos
homogéneos anisotrópicos
Acuífero muy extenso
en planta
d1 Arena K1
d2 Rocas fracturadas K2
d3 Arena K3
d4 Arcilla K4
d1 Arena K1
Acuífero muy extenso
d2 Rocas fracturadas K2
en planta
d3 Arena K3
d4 Arcilla K4
N
K d N
d
i i
Kx i 1
N
i
Kz
d
i 1
N
di
i
i 1
i 1 K i
Controlada por la capa de Controlada por la capa de
mayor conductividad menor conductividad
Usualmente: Kx Kz
Heterogeneidad y anisotropía
Ejemplo 3
La superficie piezométrica de un acuífero confinado esta
1.5m bajo la tabla de agua de un acuífero inconfinado
superior. Los dos acuíferos están separados por un
acuitardo de 0.5 m de separación. La tabla de agua en el
acuífero inconfinado está 20m sobre la base del acuitardo
y está en equilibrio con el flujo vertical debido a
infiltración (recarga) desde la superficie. La
conductividad hidráulica K para el acuífero inconfinado
es 08m/día y K para el acuitardo es 0.1 m/día. Cuál es el
flujo o caudal vertical del acuífero inconfinado al
acuífero confinado? Cuál es la velocidad del flujo si la
porosidad es 0.4?
Aplicaciones de la Ley de Darcy
1. Flujo permanente en un acuífero confinado
Considere un acuífero isotrópico y homogéneo de
conductividad hidráulica K, de longitud L, espesor b el
cual está limitado por dos embalses con cabezas constantes
hL y hR respectivamente. Cuál es el caudal que fluye a
través del acuífero confinado?
Cabeza piezométrica
Río
Canal
Acuífero
hL= 3 m Flujo confinado h R= 2 m
Aplicaciones de la Ley de Darcy
1. Flujo permanente en un acuífero confinado
El caudal permanente que fluye por unidad de longitud a
través del acuífero es Q y el caudal específico es q = Q/b =
constante.
De la Ley de Darcy, Cabeza piezométrica
dh
q K
dx
Integrando de x = 0 a x Flujo Acuífero
h q x q confinado
dh dx h h L x
h hL K x 0 K
En particular h = hR para
x = L y por lo tanto,
Donde Kd se conoce como T la
h hR h hR transmisividad del acuífero
q L K y Q L Kd
L L
Aplicaciones de la Ley de Darcy
2. Suposiciones de Dupuit - Forchheimer para flujo inconfinado
Para flujo saturado en un acuífero inconfinado la tabla de agua es
una línea de corriente. Esto implica que el acuífero limitado por dos
embalses de diferente nivel en adición al flujo horizontal tiene que
ocurrir flujo vertical. La solución al problema es complicada siendo
un problema de superficie libre con condiciones no lineales. Dupuit
(1863) desarrolló una solución simple extendida luego por
Forchheimer (1886). Su solución asume que:
a) La velocidad del flujo es proporcional a la tangente del gradiente
hidráulico = pendiente de la tabla de agua.
b) El flujo es horizontal y uniforme en todas partes de la sección
vertical (donde Q = caudal total por unidad de ancho),
dh
Q Kh
dx
Aplicaciones de la Ley de Darcy
2. Suposiciones de Dupuit - Forchheimer para flujo
inconfinado
hL
h
dh
Q Kh Q
hR
dx x
L
h Q x h 2 hL2 Q
Integrando entre x = 0 y x resulta en, hdh dx x
h hL K x 0 2 K
La tabla de agua tiene forma parabólica. En la frontera x = L se
asume que h = hR lo cual da un caudal por unidad de ancho del
acuífero inconfinado de,
Q
K 2
2L
hL hR2
Aplicaciones de la Ley de Darcy – Ejemplo 3
Ejemplo 4
Se utiliza un embalse formado mediante un terraplén de tierra de 50 m de
ancho para contener aguas residuales de una explotación minera. El nivel de
agua al interior del terraplén se mantiene a 7m sobre una base impermeable
y el nivel en un canal de drenaje perimetral se mantiene a 1.5 m sobre la
misma base impermeable. Si la conductividad hidráulica del terraplén de
tierra es de 1.4 x 10-6 ms-1, estime el caudal de estado estable o la tasa de
infiltración del agua residual a través del terraplén, por metro de longitud.
Estime la velocidad del flujo y el tiempo de residencia hidráulica del agua
residual a través del terraplén.
Determinación de la Conductividad Hidráulica en
forma experimental en el laboratorio
Tarea individual (sin entrega): Investigar sobre permeámetros
de cabeza (o carga) constante y permeámetros de cabeza o
carga variable. Cómo se determina K experimentalmente?
ho
h
L SUELO Area = a SUELO L
Area = A
Area Secc. = A
