Flujo Estacionario Unidimensional.

19

4 FLUJO ESTACIONARIO UNIDIMENSIONAL

4.1 Acuíferos Confinados.

La representación esquemática de un acuífero confinado se muestra en la Figura 4.1.

Figura 4.1 Acuífero Confinado.

De la ley de Darcy se obtiene la siguiente ecuación diferencial simple:

Q dh
 q  K (4.1)
A dx

la solución general es:

Q q
h x  Co  x  Co (4.2)
KA K

que representa una línea recta con pendiente – Q/ (KA)= –q/ K

en donde:

Co: Constante de integración, que se obtiene resolviendo las condiciones de frontera (o de

borde)

Utilizando la ecuación anterior con condiciones de borde x=0, h=ho se tiene que Co=ho.

Reemplazando:

Q Q q
h  ho  x  ho  x  ho  x (4.3)
KA Kwb K
20 Flujo y Contaminación de Aguas Subterráneas.

Siendo ésta la ecuación de una línea recta,

en donde:

w: Ancho del acuífero

b: Espesor

T: Kb (Transmisividad)

Transmisividad

Usualmente en acuíferos confinados el espesor b es constante y si el el medio es homogéneo K

es constante, por lo tanto el producto Kb también es constante conocido con el nombre de
transmisividad T. De manera que se define como la propiedad de un acuífero para transmitir el
agua en todo su espesor saturado; una definición más técnica a partir de la ley de Darcy,
transmisividad es el caudal (Q) que pasa a través de todo el espesor (b) del acuífero, en un
ancho unitario (w=1), bajo un gradiente unitario a las condiciones ambientales.

T=Kb

dh
Q  K A (4.4)
dx

dh
si w  1 y 1 (4.5)
dx

A  b *1 (4.6)

Por lo tanto la ecuación (4.5) queda:

Q  Kb  T , si w  1 y dh  1 (4.7)

dx

Los términos entre corchetes, indican que los valores numéricos son iguales, pero claramente
las unidades son diferentes.

Finalmente incluyendo la transmisividad, la ecuación (4.3) queda:

Q q'
h  ho  x  ho  x (4.8)
Tw T

Ejemplo 4.1
Flujo Estacionario Unidimensional. 21

4.2 Acuíferos Libres

La representación esquemática de un acuífero libre se muestra en la Figura 4.2.

Figura 4.2 Acuífero Libre.

Ley de Darcy:

dh
Q  K A (4.9)
dx

A  wh (4.10)

dh
Q   Kwh (4.11)
dx

reorganizando:

Q dx   Kwh dh (4.12)

Con la condición de contorno de la Figura 4.2. x=0, h=ho:

x h

 Q dx    Kw h dh
0 ho
(4.13)

Qx 
KW
2

 ho 2  h 2  (4.14)

2Q
h 2  ho 2  x (4.15)
Kw

La ecuación (4.15) se conoce como Ecuación de Dupuit; en su solución se tienen que considerar
las hipótesis de Dupuit ( y Forcheimmer), que se pueden resumir en la siguiente forma:
22 Flujo y Contaminación de Aguas Subterráneas.

La velocidad se puede tomar proporcional a la tangente del gradiente hidráulico. El flujo es

horizontal y uniforme en cualquier sección (o las líneas equipotenciales verticales).

Por lo tanto la Ecuación de Dupuit, al no tomar en cuenta las componentes verticales de la

velocidad arrojará valores del espesor saturado mayores que el real.

Estas hipótesis limitan la solución, pero son necesarias para resolver la ecuación de Laplace,
como se mostrará más adelante; además son tenidas en cuenta en muchos desarrollos
posteriores.

Es importante anotar que cuando se trabaja en acuíferos libres las cabezas se deben tomar con
respecto al fondo del acuífero, no así en acuíferos confinados en los cuales se puede definir un
datum en cualquier sitio.

Si se considera ahora que el acuífero libre tiene recarga ó infiltración (I), como se muestra en la
Figura 4.3., se pueden plantear las siguientes expresiones:

(a) Ley de Darcy:

dh
Q   Kh (4.16)
dx

considerando una franja de ancho unitario. La infiltración por elemento dx es:

dQ  Idx (4.17)

Figura 4.3 Acuífero Libre con Recarga.

(b) Continuidad:

Aplicando el principio de continuidad entre el flujo que entra al acuífero por efecto de la
recarga y el que se mueve a lo largo de él se igualan las dos expresiones

d  dh 
 Kh    I (4.18)
dh  dx  Ecuación de Boussinesq
Flujo Estacionario Unidimensional. 23

d  dh 
T   I (4.19)
dh  dx 

Si en la expresión anterior no se considera la infiltración se puede llegar a la ecuación de

Laplace de la siguiente forma:

I 0 (4.20)

d 2h
0 (4.21)
dx 2

2h  0 (4.22)
Ecuación de Laplace

Retomando el caso en el que se presenta infiltración en un medio homogéneo ( k constante ), la

Ecuación de Boussinesq obtenida se puede reescribir así:

d  dh  I
h    (4.23)
dx  dx  K

Linealizando:

d  dh 2  2I
   
dx  dx
(4.24)
 K

d  dh 2  2I
    dx
dx  dx
(4.25)
 K

Integrando:

dh 2 2I
  x  C1 (4.26)
dx K

2I
 dh 
K 
x dx  C1  dx  C 2
2
(4.27)

I 2
h2   x  C1 x  C2 (4.28)
K

Si se considera el caso particular en el cual se tienen dos cauces separados una distancia L,
como se muestra en la Figura 4.4, y se toma el origen de coordenadas x en el canal con altura
h= h1, se pueden plantear las siguientes condiciones de contorno:
24 Flujo y Contaminación de Aguas Subterráneas.

Figura 4.4 Acuífero libre con infiltración separado por dos cauces

x  0 h  h1 (4.29)

x  L h  h2 (4.30)

Lo que implica que:

C 2  h1
2
(4.31)

I 
h2  h1    L2
2 2

C1  K (4.32)
L

o sea

I 2  h2  h1 I L 
2 2

h2   x   x  h1
2
 (4.33)
K  L K 

Esta ecuación representa la curva de nivel de tabla de agua entre los dos canales.

Para encontrar el punto donde se tiene la mayor altura de tabla de agua (divisoria de aguas) se
tiene que:

dh
0 (4.34)
dx

luego se reemplaza en la ecuación (4.33) para obtener x.

Flujo Estacionario Unidimensional. 25

4.3 Acuíferos Semiconfinados

La representación grafica del acuífero semiconfinado está en la Figura 4.5, en donde se definen
las siguientes variables:

ho: Nivel de la tabla de agua por encima del estrato semipermeable

h: Cabeza del acuífero semiconfinado.

K’: Conductividad hidráulica vertical de la capa semipermeable

K: Conductividad hidráulica del acuífero semiconfinado.

T: Transmisividad del acuífero semiconfinado.

b’: Espesor de la capa semipermeable.

b: Espesor de la capa semiconfinada.

c = b’ / K’: Resistencia hidráulica [T]. Su significado físico consiste en que si se divide por
la porosidad (c/n), se obtiene el tiempo que una gota de agua tarda en pasar el estrato
semipermeable cuando está sometida a un gradiente unitario.

B2 = (Kb)c = T c: Factor de goteo [L]

Aplicando la ley de Darcy al acuífero semiconfinado:

dh
Q  T (4.35)
dx

dQ d 2h
 T 2 (4.36)
dx dx

Figura 4.5 Acuífero Semiconfinado

26 Flujo y Contaminación de Aguas Subterráneas.

Aplicando también la ley de Darcy al flujo vertical en la capa semiconfinante a través de la

sección dx se puede plantear:

h  ho
dQ   K '  1 * dx  (4.37)
b'

dQ h  ho

dx b' (4.38)
K'

dQ ho  h
 (4.39)
dx c

De nuevo se aplica el principio de continuidad, ya que el caudal que pasa por el estrato
semipermeable debe ser igual al caudal que se mueve en el estrato semiconfinado, si no se
tienen otras fuentes de recarga.

Por lo tanto igualando:

d 2 h ho  h
T  (4.40)
dx 2 c

d 2 h h  ho
 (4.41)
dx 2 B2

d 2 h  ho  h  ho 
 (4.42)
dx 2 B2

La solución a esta ecuación es de la forma:

x x

h  ho  C1 e B  C 2 e B (4.43)

los valores de C1 y C2, constantes de integración, se obtienen con las condiciones de borde
propias de cada situación particular.

En algunos textos se emplea el valor 1/c [T-1], que se conoce con el nombre de conductancia
hidráulica.

4.4 FLUJO RECTILINEO EN ACUIFEROS ESTRATIFICADOS

Una manera de considerar la heterogeneidad en los medios porosos consiste en idealizarlos

como una serie de capas paralelas y tomar el valor de la Conductividad Hidráulica de cada una
de ellas diferente. Esta esquematización es muy realista ya que las rocas de origen sedimentario,
que se forman por depositación de capas (sedimentos), constituyen los mejores acuíferos.
Flujo Estacionario Unidimensional. 27

El problema ahora consiste en determinar una Conductividad Hidráulica equivalente para toda
la estratigrafía, lo que lleva a considerar dos situaciones: Cuando el flujo es paralelo a la
estratigrafía y cuando este es perpendicular a ella.

4.4.1 Flujo Paralelo a la Estratificación

La Figura 4.6 presenta las condiciones a analizar.

Aplicando la Ley de Darcy al caudal que fluye por cada estrato:

h
Q1   K 1 b1 (4.44)
x

h
Q2   K 2 b2 (4.45)
x

Figura 4.6 Flujo paralelo a la estratificación

para el n’esimo estrato:

h
Q n  K n bn (4.46)
x

El caudal total se puede expresar de acuerdo a la Ley de Darcy en términos de una K

equivalente:

h
Qt   K eq
x
 bi (4.47)

a su vez será igual a la suma de cada uno de los diferentes caudales que pasan a través de cada
estrato:
28 Flujo y Contaminación de Aguas Subterráneas.

Q t   Qi (4.48)

reemplazando:

h  h 
K eq
x
 bi    K i
 x 
bi  (4.49)

K eq 
 K b 
i i
(4.50)
 b i

Teq  Ti (4.51)

Nótese que la conductividad hidráulica Keq es un promedio ponderado.

4.4.2 Flujo Perpendicular a la Estratificación

La Figura 4.7 presenta las condiciones a analizar.

Aplicando la Ley de Darcy a cada estrato:

h1
Q1   K 1 A (4.52)
b1

h2
Q2   K 2 A (4.53)
b2

para el n’esimo estrato:

h2
Q2   K 2 A (4.54)
b2

En general para un estrato i se puede plantear que:

hn
Qn   K n A (4.55)
bn

Tomando la columna como un solo medio poroso con un valor de K equivalente resulta que:

h
Qt   K eq A (4.56)
b
Flujo Estacionario Unidimensional. 29

Qt b
h   (4.57)
K eq A

Figura 4.7 Flujo perpendicular a la estratificación

Para estas condiciones es válido que:

h   hi (4.58)

b   bi (4.59)

Q t  Q1  Q2    Qn (4.60)

por lo tanto se puede replantear lo escrito como:

 Qi bi  Qt b
  K   (4.61)
 i A  K eq A

Qt  bi  Qt  b
  K   (4.62)
A  i  A K eq
30 Flujo y Contaminación de Aguas Subterráneas.

Keq 
 b 
i

 bi  (4.63)
  K 
 i

La forma de la ecuación (4.63) corresponde a un promedio armónico, que en modelos

numéricos de flujo tiene también bastante aplicación.