142 Capítulo 7: Filtración

C A P Í T U L O 7
Filtración

7.1 Introducción
En el capítulo 6 se consideraron algunos casos simples para los que se requiere la aplicación
directa de la ley de Darcy para calcular el flujo de agua a través del suelo. En muchos casos el
flujo de agua a través del suelo no es sólo en una dirección y no es uniforme sobre toda el área
perpendicular al flujo. En tales casos el flujo de agua subterránea se calcula generalmente por
el uso de gráficos que se refieren como redes de flujo. El concepto de la red de flujo se basa en
la ecuación de continuidad de Laplace, que rige la condición de flujo constante para un punto
dado en la masa de suelo. En este capítulo vamos a deducir la ecuación de continuidad de La-
place y aplicarla a la elaboración de redes de flujo.

7.2 Ecuación de continuidad de Laplace

Para deducir la ecuación diferencial de continuidad de Laplace, consideremos una sola fila de
pilotes que han sido clavados en una capa de suelo permeable, como se muestra en la figura
7.1a. Se supone que la hilera de pilotes es impermeable. El flujo de agua en estado estacionario
del lado aguas arriba hacia el lado aguas abajo a través de la capa permeable es un flujo de dos
dimensiones. Para el flujo en un punto A, consideramos un bloque elemental de suelo. El bloque
tiene dimensiones dx, dy y dz (la longitud dy es perpendicular al plano del papel); se muestra
en una escala ampliada en la figura 7.1b. Sean vx y vz las componentes de la velocidad de des-
carga en las direcciones horizontal y vertical, respectivamente. El caudal de agua en el bloque
elemental en la dirección horizontal es igual a vx dz dy, y en la dirección vertical es vz dx dy. Los
tipos de flujo de salida desde el bloque en las direcciones horizontal y vertical son

avx
0vx
dxb dz dy
0x
y
0vz
avz dzb dx dy
0z

142
 http://estudiantesingcivil.blogspot.mx/
7.2 Ecuación de continuidad de Laplace 143

Pilote

H1 h
H2

dz A

dx

Capa
impermeable

(a)

(y +z
yz
z )
dz dx dy

dy

yx dz dy (y +
x
yx
x )
dx dz dy

dz

dx
yz dx dy

(b)

Figura 7.1 (a) Pilotes en una sola fila clavados en una capa permeable; (b) flujo en A

respectivamente. Suponiendo que el agua es incompresible y que no se produce ningún cambio

en el volumen en la masa de suelo, sabemos que el caudal total de entrada debe ser igual al
caudal total de salida. Por lo tanto,

0vz
c a vx dx b dz dy a vz dz b dx dy d
0vx
[vx dz dy vz dx dy] 0
0x 0z
144 Capítulo 7: Filtración

o
0vx 0vz
0 (7.1)
0x 0z

Con la ley de Darcy, las velocidades de descarga pueden expresarse como

kx a b
0h
vx kxix (7.2)
0x
y

kz a b
0h
vz kziz (7.3)
0z

donde kx y kz son las conductividades hidráulicas en las direcciones horizontal y vertical, res-
pectivamente.
De las ecuaciones (7.1), (7.2) y (7.3), podemos escribir

0 2h 0 2h
kx kz 0 (7.4)
0x2 0z2

Si el suelo es isotrópico con respecto a la conductividad hidráulica, es decir kx 5 kz, la

ecuación de continuidad anterior para flujo en dos dimensiones se simplifica a

0 2h 0 2h
0 (7.5)
0x2 0z2

7.3 Redes de flujo

La ecuación de continuidad [ecuación(7.5)] en un medio isotrópico representa dos familias de
curvas ortogonales: las líneas de flujo y las líneas equipotenciales. Una línea de flujo es una
línea a lo largo de la cual una partícula de agua se desplazará desde el lado aguas arriba hacia
el lado aguas abajo en un medio de suelo permeable. Una línea equipotencial es una línea a lo
largo de la cual el potencial de carga en todos los puntos es igual. Por lo tanto, si se colocan
piezómetros en diferentes puntos a lo largo de una línea equipotencial, el nivel de agua subirá
a la misma elevación en todos ellos. La figura 7.2a muestra la definición de flujo y líneas equi-
potenciales para el flujo de la capa de suelo permeable alrededor de la fila de pilotes mostradas
en la figura 7.1 (para kx 5 kz 5 k).
A la combinación de un número de líneas de flujo y líneas equipotenciales se le llama red
de flujo. Las redes de flujo se construyen para calcular el flujo de las aguas subterráneas en el
medio. Para completar la construcción gráfica de una red de flujo se debe dibujar el flujo y las
líneas equipotenciales, de tal manera que las líneas equipotenciales intersecten a las líneas de
flujo en ángulo recto y los elementos de flujo formados son cuadrados aproximados.
La figura 7.2b muestra un ejemplo de una red de flujo completa. Otro ejemplo de una red
de flujo en una capa permeable isotrópica se muestra en la figura 7.3. En estas figuras, Nƒ es
el número de canales de flujo en la red y Nd es el número de caídas de potencial (definido más
adelante en este capítulo).
 7.3 Redes de flujo 145

Pilote

H
H

Línea de flujo
kx kz k
Línea
equipotencial

Capa
impermeable

Pilote
Nivel de agua

H Nivel de agua
b a d H e

kx kz k c
Nf
Nd

f g

Capa
impermeable

Figura 7.2 (a) Definición de líneas de flujo y líneas equipotenciales, (b) red de flujo completa
146 Capítulo 7: Filtración

H
H1
H2

Filtro de punta

kx kz k
Nf 5
Nd 9

Figura 7.3 Flujo neto bajo un dique con filtro de punta

Dibujar una red de flujo toma varios ensayos. Mientras se construye la red de flujo, deben
mantenerse las condiciones de contorno en mente. Para la red de flujo mostrada en la figura
7.2b, se aplican las cuatro condiciones de contorno siguientes:
1. Las superficies de aguas arriba y aguas abajo de la capa permeable (líneas ab y de) son
líneas equipotenciales.
2. Debido a que ab y de son líneas equipotenciales, todas las líneas de flujo se intersectan en
ángulo recto.
3. El límite de la capa impermeable, es decir, la línea fg, es una línea de flujo y también lo es
la superficie del pilote impermeable, la línea acd.
4. Las líneas equipotenciales intersectan a las líneas acd y fg en ángulo recto.

7.4 Cálculo de la filtración a partir de una red de flujo

En toda red de flujo, la franja entre dos líneas de flujo adyacentes se llama canal de flujo. La
figura 7.4 muestra un canal de flujo con las líneas equipotenciales formando elementos cuadra-
dos. Sean h1, h2, h3, h4, . . . , hn los niveles piezométricos correspondientes a las líneas equipo-
tenciales. La tasa de filtración a través del canal de flujo por unidad de longitud (perpendicular
a la sección vertical a través de la capa permeable) se puede calcular de la siguiente manera:
debido a que no hay flujo a través de las líneas de flujo,

¢q1 ¢q2 ¢q3 p ¢q (7.6)

h1

h2
∆q
h3
l1 h4
l1

l2
l2 l3
∆q 1 l3

∆q 2 ∆q
∆q 3

Figura 7.4 Filtración a través de un canal de flujo con elementos cuadrados

 7.4 Cálculo de la filtración a partir de una red de flujo 147

A partir de la ley de Darcy, la velocidad de flujo es igual a kiA. Por lo tanto, la ecuación (7.6)
se puede escribir como

ka b l1 ka b l2 ka b l3
h1 h2 h2 h3 h3 h4 p
¢q (7.7)
l1 l2 l3
La ecuación (7.7) muestra que si los elementos de flujo se dibujan como cuadrados aproxima-
dos, entonces la caída en el nivel piezométrico entre dos líneas equipotenciales adyacentes es la
misma. Esto se conoce como caída de potencial. Por lo tanto,

p H
h1 h2 h2 h3 h3 h4 (7.8)
Nd
y
H
¢q k (7.9)
Nd
donde

H 5 diferencia de carga entre el lado de aguas arriba y el de aguas abajo

Nd 5 número de caídas de potencial

En la figura 7.2b los elementos de flujo son cuadrados aproximados. Para cualquier canal de
flujo, H 5 H1 ! H2 y Nd 5 6.
Si el número de canales de flujo en una red de flujo es igual a Nƒ, el caudal total a través
de todos los canales por unidad de longitud puede ser dado por
HNf
q k (7.10)
Nd
Aunque dibujar los elementos cuadrados para una red de flujo es conveniente, no siempre
es necesario. Alternativamente, se puede dibujar una malla rectangular para un canal de flujo,
como se muestra en la figura 7.5, a condición de que las razones de anchura a longitud para
todos los elementos rectangulares en la red de flujo sea la misma. En este caso, la ecuación (7.7)
para el caudal a través del canal puede ser modificada para

ka b b1 ka b b2 ka b b3
h1 h2 h2 h3 h3 h4 p
¢q (7.11)
l1 l2 l3

h1

h2
∆q b1
l1 h3
h4
b2
l2 b3
l3
∆q
b 1 b 2 b 3 ...
= = = =n
l1 l2 l3

Figura 7.5 Filtración a través de un canal de flujo con elementos rectangulares

148 Capítulo 7: Filtración

Si b1/l1 5 b2/l2 5 b3/l3 5 . . . 5 n (es decir, los elementos no son cuadrados), las ecuaciones (7.9)
y (7.10) pueden ser modificadas:

kH a b
n
¢q (7.12)
Nd
o

Nf
q kH a bn (7.13)
Nd

La figura 7.6 muestra una red de flujo de filtraciones en torno a una sola fila de pilotes.
Tenga en cuenta que los canales de flujo 1 y 2 tienen elementos cuadrados. Por lo tanto, el cau-
dal a través de estos dos canales se puede obtener de la ecuación (7.9):

H H H
¢q1 ¢q2 k k 2k
Nd Nd Nd

Sin embargo, el flujo del canal 3 tiene elementos rectangulares. Estos elementos tienen una re-
lación entre anchura y longitud de alrededor de 0.38, por lo que, de la ecuación (7.12), tenemos

kH a b
0.38
¢q3
Nd

Así, la tasa total de la filtración se puede dar como

kH
q ¢q1 ¢q2 ¢q3 2.38
Nd

Nivel de agua

H
Nivel freático

Superficie del suelo

Canal de flujo 1
l
=1
b
Canal de flujo 2
l
=1
b
Canal de flujo 3
l 1
≈
b 0.38
Superficie
impermeable

Figura 7.6 Red de flujo para filtración en torno a una sola fila de pilotes
 7.4 Cálculo de la filtración a partir de una red de flujo 149

Ejemplo 7.1
En la figura 7.7 se muestra una red de flujo para el flujo en torno de una sola fila de pilotes
en una capa de suelo permeable. Tenemos que kx 5 kz 5 k 5 5 " 10!3 cm/s.
a. ¿Qué tan alto (por encima de la superficie del suelo) crecerá el agua si se colocan piezó-
metros en los puntos a, b, c y d?
b. ¿Cuál es la tasa de filtración a través de canal de flujo II por unidad de longitud (perpen-
dicular a la sección mostrada)?
Solución
a. En la figura 7.7 vemos que Nf 5 3 y Nd 5 6. La diferencia de carga entre el lado de
aguas arriba y aguas abajo es 3.33 m, por lo que la pérdida de carga por cada gota es
3.33/6 5 0.555 m. El punto a está situado en la línea equipotencial 1, lo que significa
que la caída de potencial en a es 1 " 0.555 m. El agua en el piezómetro en a subirá a
una elevación de (5 ! 0.555) 5 4.445 m por encima de la superficie del suelo. Del
mismo modo, podemos calcular el resto de los niveles piezométricos:
b 5 (5 ! 2 " 0.555) 5 3.89 m por encima de la superficie del suelo
c 5 (5 ! 5 " 0.555) 5 2.225 m por encima de la superficie del suelo
d 5 (5 ! 5 " 0.555) 5 2.225 m por encima de la superficie del suelo
b. De la ecuación (7.9), tenemos
H
¢q k
Nd
3 5
k 5 10 cm/s 5 10 m/s
¢q (5 10 5)(0.555) 2.775 10 5
m3/s/m

Pilote

5m
Superficie del suelo 1.67 m
0 6
I
c II

III
10 m a

b
1 2 3 4 5 d

Superficie
impermeable

Figura 7.7 Red de flujo para el flujo en torno a una fila de pilotes en una capa de suelo permeable
150 Capítulo 7: Filtración

7.5 Redes de flujo en un suelo anisotrópico

La construcción de una red de flujo descrita hasta el momento y las ecuaciones deducidas (7.10)
y (7.13) para el cálculo de la filtración se han basado en la suposición de que el suelo es isotrópi-
co. Sin embargo, en la naturaleza la mayoría de los suelos presenta algún grado de anisotropía.
Para tener en cuenta la anisotropía del suelo con respecto a la conductividad hidráulica, hay que
modificar la construcción de la red de flujo.
La ecuación diferencial de continuidad para un flujo de dos dimensiones [ecuación (7.4)] es

0 2h 0 2h
kx kz 0
0x2 0z2
Para suelos anisotrópicos, kx Z kz. En este caso, la ecuación representa dos familias de
curvas que no se unen a 90°. Sin embargo, podemos reescribir la ecuación anterior como

0 2h 0 2h
0 (7.14)
(kz /kx)0x2 0z2

Sustituyendo x 2kz kx x, podemos expresar la ecuación (7.14) como

0 2h 0 2h
0 (7.15)
0x¿ 2 0z2

Ahora la ecuación (7.15) está en una forma similar a la de la ecuación (7.5), con x remplazada
por x', que es el nuevo sistema de coordenadas transformado. Para construir la red de flujo,
utilice el siguiente procedimiento:

Paso 1: Adopte una escala vertical (es decir, eje z) para el dibujo de la sección transversal.
Paso 2: Adopte una escala horizontal (es decir, el eje x) de tal manera que escala hori-
zontal 2kz kx escala vertical.
Paso 3: Con escalas adoptadas como en los pasos 1 y 2, trace la sección vertical a través
de la capa permeable paralela a la dirección del flujo.
Paso 4: Dibuje la red de flujo de la capa permeable de la sección obtenida en el paso 3,
con las líneas de flujo intersectando a las líneas equipotenciales en ángulo recto
y los elementos como cuadrados aproximados.

La tasa de filtración por unidad de longitud se puede calcular mediante la modificación

de la ecuación (7.10)

HNf
q 2kx kz (7.16)
Nd

donde H 5 pérdida de carga total

Nf y Nd 5 número de canales de flujo y caídas potenciales, respectivamente (de la red
de flujo elaborada en el paso 4)

Tenga en cuenta que cuando las redes de flujo se dibujan en secciones transformadas (en
suelos anisotrópicos), las líneas de flujo y las líneas equipotenciales son ortogonales. Sin em-
bargo, cuando son redibujadas en una sección alineada, estas líneas no están en ángulos rectos
 7.5 Redes de flujo en un suelo anisotrópico 151

kz 1
=
kx 6
Escala vertical 6m

Escala horizontal
6( 6) 14.7 m

(a)

Escala 6 m

(b)

Figura 7.8 Elemento de flujo en un suelo anisotrópico: (a) en una sección transformada; (b) en una sec-
ción alineada

entre sí. Este hecho se muestra en la figura 7.8. En esta figura, se supone que kx 5 6kz. La figura
7.8a muestra un elemento de flujo en una sección transformada. El elemento de flujo ha sido
rediseñado en una sección alineada en la figura 7.8b.

Ejemplo 7.2
En la figura 7.9a se muestra una sección de la presa. Las conductividades hidráulicas de
la capa permeable en las direcciones vertical y horizontal son 2 " 10!2 mm/s y 4 " 10!2
mm/s, respectivamente. Dibuje una red de flujo y calcule la pérdida de filtración de la presa
en m3/día /m.
 http://estudiantesingcivil.blogspot.mx/
152 Capítulo 7: Filtración

6.1 m

7.6 m (a)

6.1 m

1.0
1.0
0.5

Escala horizontal 7.6 2 10.75 m

Escala vertical 7.6 m

Capa permeable Capa impermeable

(b)

Figura 7.9

Solución
A partir de los datos dados,
kz 5 2 10!2 mm/s 5 1.728 m/día
kx 5 4 10!2 mm/s 5 3.456 m/día
y h 5 6.1 m. Para trazar la red de flujo,

B4
2
2 10
Escala horizontal 2
(escala vertical)
10

22
1
(escala vertical)

Sobre la base de esto, la sección de la presa vuelve a representarse, y la red de flujo se

dibuja como en la figura 7.9b. La tasa de filtración está dada por q 2kx kz H(Nf Nd).
De la figura 7.9b, Nd 5 8 y Nf 5 2.5 (el canal de flujo más inferior tiene una anchura con
longitud de 0.5). Por lo tanto,

q 2(1.728)(3.456) (6.1)(2.5/8) 4.66 m3/día /m

 Problemas 153

7.6 Resumen

1. Para la condición de flujo de dos dimensiones, la ecuación de continuidad de Laplace se

puede dar como [ecuación (7.4)]:

0 2h 0 2h
kx kz 0
0x2 0z2

2. Para la condición isotrópica con respecto a la conductividad hidráulica, kx 5 kz. Así

0 2h 0 2h
0
0x2 0z2

3. Las redes de flujo se pueden dibujar usando la ecuación de continuidad de Laplace.

4. Las líneas de flujo y líneas equipotenciales son dos familias de curvas. Cuando kx 5 kz, las
líneas de flujo y las líneas equipotenciales son dos familias de curvas ortogonales.
5. La tasa de filtración (q) bajo estructuras hidráulicas (kx 5 kz 5 k) se puede escribir como
[ecuación (7.13)]:
Nf
q kH a bn
Nd
donde n 5 razón anchura-longitud de los elementos de flujo.

Problemas
7.1 Consulte la figura 7.10 y utilice estos valores:
H1 5 7 m D 5 3.5 m
H2 5 1.75 m D1 5 7 m

Pilote

H1
H2

D1
k = 6.5 × 10−4 cm/s

Capa
impermeable

Figura 7.10
154 Capítulo 7: Filtración

29.5 m

1.7 m
10 m

Pilote 7.2 m
2.4 m
19.5 m

1.02 × 10 −5 m/s

Capa impermeable

Figura 7.11

Dibuje una red de flujo. Calcule la pérdida de la filtración por metro de longitud del pilote
(en ángulo recto a la sección transversal mostrada).
7.2 Dibuje una red de flujo para una sola fila de pilotes clavados en una capa permeable, como
se muestra en la figura 7.10, dada la siguiente información:

H1 5 5 m D54m
H2 5 0.7 m D1 5 10 m

Calcule la pérdida de la filtración por metro de longitud del pilote (en ángulo recto con
respecto a la sección transversal mostrada).
7.3 Dibuje una red de flujo para el vertedero que se muestra en la figura 7.11. Calcule la tasa
de filtración debajo del vertedero.