2.5. EL EXPERIMENTO DE DARCY
2.5.1. FORMA EXTENDIDA DE LA LEY DE DARCY
2.5.2. EJEMPLOS DEL CÁLCULO DE LA VELOCIDAD DE FLUJO DE AGUAS
SUBTERRÁNEAS EN DOS DIMENSIONES
2.5.3. CONCEPTOS ADICIONALES DE FLUIDOS POTENCIALES
2.5.1. FORMA EXTENDIDA DE LA LEY DE DARCY



El experimento de Darcy relaciona la descarga total Q del gradiente
de carga, el cual es Dh/l.
El factor de proporcionalidad es KA, donde K es la conductividad
hidráulica y A es el área de la columna por donde el agua fluye.
La descarga especifica o flujo volumétrico, es el fluido que pasa a
través del medio porosos en una sección de área A perpendicular al
flujo por unidad de tiempo:
q
Q
A

Esta medida puede ser separada por dl, por un cambio en dh:
q 
Q
A
K
dh
dl
q  K
dh
dl


La ecuación anterior es una relación constitutiva y es equivalente al
momento en la ecuación de balance.
La relación clásica entre la ecuación de momento de balance para un
fluido en medios porosos fue dada por Hubbert en 1954 y a partir de
la derivación la ley de Darcy a través de la ecuación de Neiver-Stock,
se llega a la expresión de flujo de un fluido en un material saturado:
q    Nd

2
Donde el termino h es igual a
  g   h
h 
h
x
i
h
y
j
h
z
k
Figura 2.26. Definición esquemática del concepto de vector unitario

Hubbert definió el producto de Nd2 como la permeabilidad del
medio (k), la cual depende de la geometría de los granos y
representa la conductividad hidráulica K como:
K   Nd

2
  g 
kg

Si el medio es isotrópico se dice que K es un escalar, pero si el
medio es anisotrópico, entonces se tiene que adaptar el concepto de
la K a una matriz que contenga diferentes valores en diferentes
puntos, los valores de la matriz se le llama tensor
q  K h
 K xx

K  0

 0
0
K yy
0
0 

0

K zz 
q x 
 K xx
 

q y     0
q 
 0
 z
0
K yy
0
0  h x 


0 h y 

K zz    h  z 
Figura 2.27. Ilustración de los efectos de la litología en la conductividad
hidráulica cuando los ejes coordenados estas orientados igual a los estratos
 K xx
q x 

 
q


 y
 K yx
q 
K
 z
 zx
K xy
K yy
K zy
K xz    h  x 


K yz    h  y 
K zz    h  z 
Figura 2.28. Ilustración de los efectos de la litología en la conductividad
hidráulica cuando los ejes coordenados no estas orientados igual a los estratos


Las dimensiones del vector de la descarga especifica q es [L/T],
esta no es una medida de la velocidad del agua, sino es una medida
del volumen que pasa a través de una superficie de área A en
tiempo DT dividido por el área A y DT.
Puesto que la porosidad es la relación de espacios vacíos y espacio
total, el área del agua puede ser expresado como el área total
multiplicado por la porosidad, así obtenemos la velocidad para una
partícula de agua en un medio poroso la cual se expresa con la
relación:
v 
q

Índice
2.5.2. EJEMPLOS DEL CÁLCULO DE LA VELOCIDAD DE
FLUJO DE AGUAS
DIMENSIONES
SUBTERRÁNEAS
EN
DOS
Figura 2.30. Experimento de flujo en espacio de dos dimensiones. El agua subterránea
se mueve de derecha a izquierda a través de la caja llenada con arena. El nivel del agua
en la arena esta denotada por la elevación en los manómetros.
Tenemos
h ( x )  a  bx
h3  a  b 0  a
h 1  a  b  L 1 3
Resolvemos para obtener a y b
a  h3
b   h 3  h 1  L 1 3
Los sustituimos en la primera ecuación, tenemos que:
h  h 3  h 1  h 3   x L 1 3
De la ecuación anterior la localización de h=h2 y resolvemos para la
posición x2 tenemos que:
x 2  h 2  h 3  h 1  h 3   L 1 3
X2
Si consideramos el problema como sistema de ecuaciones lineales
tenemos que:
h ( x 1 , y 1 )  a  bx 1  cy 1
h ( x 2 , y 2 )  a  bx 2  cy 2
h ( x 3 , y 3 )  a  bx 3  cy 3
Para obtener los coeficientes b y c tenemos:
b
h
x
1
1
 h 2  y 2  y 3   h 2  h 3  y 1  y 2 
 x 2  y 2  y 3    x 2  x 3  y 1  y 2 
c
h
y
1
1
 h 2  x 2  x 3   h 2  h 3  x 1  x 2 
 y 2  x 2  x 3    y 2  y 3  x 1  x 2 
Habiendo obtenido los coeficientes, podemos determinar del gradiente:
El gradiente de aguas subterráneas
esta dado por:
h 
Figura 2.30. Nodal arreglo para el uso aproximaciones algebraicas
en el cálculo de gradientes
h
x
 b
h
x
i
h
y
h
y
j
c
Figura 2.31. Ejemplo del problema mostrando el cálculo de la constante
de la línea de nivel estático y la resultante del vector velocidad
Si tenemos
h
x
b
1
h
1000
y
c
1
1000
Para el triangulo inferior donde K=4ft/día:
q x   K xx 
h
x
 0 . 004
q y   K yy 
Si la =0.25, la velocidad es: v 
v 1  ( 0 . 016 , 0 . 016 ) ft día
h
x
 0 . 004
q

v 1  . 0226 ft día


Para el triangulo inferior donde K=2ft/día:
q x   K xx 
h
x
 0 . 002
q y   K yy 
Si la =0.25, la velocidad es: v 
v 2  ( 0 . 008 , 0 . 008 ) ft día

h
x
 0 . 002
q

v 2  . 0113 ft día

Índice
2.5.3. CONCEPTOS ADICIONALES DE FLUIDOS
POTENCIALES

Revisando el concepto de fluido potencial, el cual esta definido por
la ecuación:
h  z 0 dz  Patm
z

p
d
g   
Se remplaza la variable de presión por la variable p, y se introduce
en la ley de Darcy:

d 
z
p

q   K   h   K    h  z 0 dz  Patm
g    


Para la evaluación posterior de la expresión es necesario introducir
una relación matemática que describe como diferenciar una integral,
la cual se le conoce como la regla de Leibnitz:
 a ( x ) f ( x ,  ) d   a ( x )  f ( x ,  ) d   f ( x , b ( x ))   b  f ( x , a ( x ))   a
b(x)

b(x)
Utilizando la regla de Leibnitz en la ecuación q, es igual a
 z(x)
d 
p(x)

 K    z 0 ( x ) dz  Patm ( x )
g    


  K  z 0 ( x )  1dz  1 z  1 z 0
z(x)

 P(x)
 P atm 
d
P
  K   Patm ( x ) 



 g    g P   g Patm  


Si las condiciones iniciales z0=0 y Patm=0 entonces:

P 
q   K   z 

 g P  

q 
K
 g P 
  g  P  z   P 
Índice
Descargar

2.5.1. FORMA EXTENDIDA DE LA LEY DE DARCY