BASES DEL FLUJO DE AGUA SUBTERRANEA

Conductividad hidráulica :Es una variable de cuan fácil fluye el agua por un material.
Depende de permeabilidad intrínseca, densidad y la gravedad ,donde el agua más
cálida presenta mayor k debido a que presenta una viscosidad dinámica más pequeña.

La porosidad efectiva n=Vp/Vt

Acuíferos y acuitardos

Capas con baja permeabilidad son llamados acuitardos ,capas con fugas o capas
semiconfinada .Una verdadera impermeabilidad NO existe .

Con respecto a los acuíferos .Existen 2 tipos :

Acuíferos no confinados : Donde el dominio saturado es el contorno del topo de la napa

freática .

Acuíferos confinados :Donde se separa de la superficie por un aquitardo .El flujo se

denomina confinado cuando la cabeza está por encima de la parte superior del acuífero
.Ojito , el head en acuíferos confinados puede caer por debajo de la parte superior del
acuífero lo que da lugar a un flujo no confinado

El flujo en acuitardos puede a menudo aproximarse como vertical. Cuánto mayor es el

espesor de un acuitardo más difícil es que el agua fluya a través de él. Por lo que la
proporción del espesor y la conductividad hidráulica vertical es llamada resistencia
hidráulica y es una medida de impermeabilidad de esta capa.

En acuíferos un mayor espesor indica que el agua puede transmitirse más. La habilidad
de transmitir agua en dirección horizontal es llamada transmisividad.

Almacenamiento y flujo transitorio

Cuando no hay cambios con el tiempo ,el volumen almacenado permanece constante
en el tiempo ,esto implica que entrada y salida de flujo de agua son considerados
iguales. El coeficiente de almacenamiento es igual al volumen que se puede extraer
que se obtiene por 1m3 de un material .

El almacenamiento en un acuífero no confinado incrementa cuando la napa freática

aumenta. En pocas palabras es igual a la porosidad eficaz.

∆V =SA ∆h , donde S es el coeficiente de almacenamiento de acuífero no confinado. A

es área y ∆h es la altura que aumenta por el aumento de napa freática.

El almacenamiento de un acuífero no confinado es referido como coeficiente de

almacenamiento freático o campo específico y es representado por Sp .

El almacenamiento en un acuífero confinado es mas complicado ,debido a que todos

los poros están enteramente llenados con agua. Sin embargo, se puede almacenar
agua a través de la expansión del espacio de poros y aunque solo en pequeña
medida .Cuando el head aumenta por ∆h ,el volumen ∆V pero el coeficiente de
almacenamiento es mucho mas pequeño que en acuíferos no confinados.
En la mayoría de los acuíferos no consolidados la habilidad de expandir es
significativamente mas alta que la habilidad de compresión del agua. El coeficiente de
almacenamiento de acuíferos confinados depende del espesor del acuífero .

S=SsH Donde Ss es el almacenamiento especifico de un acuífero (almacenamiento por

unidad de espesor) .Se observa que es adimensional debido a que hay relación de
volumen de agua /volumen roca

Se observa que cuando se descarga agua en acuífero no confinado el nivel freático

también disminuye, en acuífero confinado la superficie piezométrica disminuye.

Es muy importante debido a que el almacenamiento junto a la calidad son principales

factores de decisión de explotar dicho acuífero

Ley de Darcy para modelamiento de agua subterránea

Descarga por unidad de área es q=Q/Aq = k *∆h /L se denomina a la descarga

especifica y lo mal refieren a velocidad de Darcy ya que no es una velocidad sino
descarga por unidad de área.

El flujo de agua subterránea puede estar en alguna dirección y es expresado en

términos de descarga específica en el vector x,y,z con componentes en cada de las 3
direcciones.

En algunos acuíferos la conductividad hidráulica

es anisotrópica y la formula puede reescribirse de la siguiente manera

La velocidad promedio puede obtenerse dividiendo q/n ,donde q es la descarga

específica y n la porosidad.

Dupuit y Forchheimer

Una simplificación de la dimensionalidad del agua subterránea para facilitar su análisis

y modelado.

En escala regional la mayoría de los acuíferos son relativamente delgados en

comparación con su extensa área y por lo tanto referidos como acuíferos poco
profundos. El vector de descarga es el caudal total del acuífero integrado sobre el
espesor saturado H del acuífero. En otras palabras es la descarga por ancho unitario del
acuífero

Flujo unidimensional constante con transmisividad constante

El vector descarga Q x es la descarga

por unidad de espesor (L /T)
2

La recarga es la descarga por unidad de área (L/T)

N es recarga

Donde T es transmisividad

A y B son constantes a determinar de condiciones de contorno

Condiciones de frontera : Tipo Dirichlet (valor h conocido) ,tipo Neumann (flujo
conocido),mixto

Con respecto a la ecuación que abarca A y B el tipo de frontera depende del problema.

Se puede observar que es una función para observar como se comporta h.

El parámetro B sale cuando x=0 y A ,cuando x=L

Flujo entre dos ríos

Donde la recarga es 0 ,las condiciones de contorno me brindan información sobre h 0 y

h1

h|x=0 =h0 h|x=L=hL

Se observa que QxT es constante

Codigo en Python:

Se aplica las formula para hallar el h

para cada 100 metros y vemos que el
Qx es constante.

Indica la descarga específica en cada altitud desde el

punto 0 a L
Para lo ultimo se usa streamplot .Donde se tiene primero que colocar una malla x e y ;
luego componentes de flujo.

Recarga areal entre dos rios

La conductividad en acuífero no confinado es posible aproxime a constante e igual a kH

mientras que la elevación de nivel freático es pequeña en comparación con el espesor
total.

Un acuífero no confinado esta limitado a la izquierda y derecha por dos rios paralelos
que penetran totalmente el [Link] rios están a distancia L y en contacto directo
con el acuífero.

h|x=0 =h0 h|x=L=hL

T es aprox constante

De esta manera obtenemos la función de h y de Q x

Pero observamos como el head se expresa de una manera diferente ,así como el Q x

Ejemplo 3 : Considera el caso que h0=h1 =10m ,el head entre los dos rios es 9
m .Computariza la recarga para este caso usando los parámetros anteriores .Graficar h
vs x

Cuando hay N ,cambia

todo.
Descarga específica vertical qz . Relacionamos ∂ qx / ∂ x + ∂ qz / ∂ z =0

Entonces se hace derivada y solo nos quedamos con

Donde una derivada crea una

función

De descarga específica vertical.

Área de recarga entre un limite impermeable y un río

Se observa que el head esta en caída ,así como que el flujo Qx inicia en 0 y
aumenta ,no hay Qx negativo.

Al siguiente ejemplo se le añade una capa de muy baja permeabilidad entre el acuífero
y el río

Descarga específica en la capa de baja

permeabilidad.

c es la resistencia hidraulica y es una dimension de tiempo

C es conduntancia del lecho del

río y es H/c

Para este problema el flujo en el rio es conocido y es igual a la recarga total NL

Nos dan datos extras de la capa de baja permeabilidad hR y c

Fujo a través de dos zonas de diferentes transmisividades

Datos : k,H y T son

diferentes .

No hay recarga debido a

que hay capa confinada.

Cond. de continuidad :interfaces entre diferentes capas de acuíferos o entre acuíferos

y no acuíferos

Por lo que si no hay N se elimina y nos quedamos con esta función para CADA zona

Por lo que hay 2 condiciones de limite y dos condiciones de continuidad

Observe que h y Q se
calcula con todos los
parámetros menos con hLo
que es el head en la
transición de T0 y T1

Se observa que el Q es constante debido a que no hay N es decir la capa esta

confinada. Asimismo es posible calcular la Transmisividad efectiva que en sí es la media
armónica de las T0 y T1 y el valor es mas bajo que la media aritmética de ambas
transmisividades.
CAPITULO 2
FLUJO SEMICONFINADO UNIDIMENSIONAL ESTACIONARIO

Se basa en el modelo donde el acuífero no esta cubierto por el tope por una capa
confinada sino una capa que presenta baja permeabilidad ,por lo que permite
infiltración.

Vemos en el grafico dos modelos .Uno cuando el h* sobrepasa el leaky layer y otro
donde el h* esta en el leaky layer . El termino h* hace referencia a la cabeza hidrauica
fija o nivel piezométrico constante.

Como no cambia el nivel significa que es la referencia o el nivel de presión del agua por
encima de la capa semiconfinante es decir h* no depende del flujo del agua dentro del
acuífero.

c es una medida de que tan difícil es el flujo vertical a través de la capa

seminconfinante ,un mayor c ,mas difícil es el flujo vertical por dicha capa.

T es la transmisividad del acuifero subyacente y sus valores se basan en dicho acuífero

KyH

Leaky factor es una medida que combina la resistencia al flujo vertical en capa
confinante como la transmisividad del acuifero [Link] que da una idea de que
tan fácilmente fluye el agua a través del sistema que considera la capa semiconfinante
como el acuífero subyacente.

Formulas a usar :

Flujo de un canal a un área drenada (polder)

onsiderar un flujo de un canal a un [Link] acuifero es semiconfinado ,la cabeza en el
polder es h* .Se usa la ecuación de Helmholtz.

Donde se demuestra que B es 0 y A es h0-h*

Donde labda es leaky layer y es 200 ,200*5 =1000 metros .Se observa como el head
desciende en forma parabólica y también la descarga

Flujo entre un lago y una área drenada

Se observa un lago y un área [Link]

Están separados del acuífero por una

Capa [Link] la resistencia en

c0 es la resistencia bajo el lago y c1 resistencia

bajo el agua [Link] lago esta separado por

un dique de insignificante ancho.

Aquí dicho dique es x=0 .La transmisividad del acuífero es constante.

Flujo para un rio en un sistema de dos acuíferos

Muchos acuíferos están limitados en el tope por un leaky layer con otro acuífero por
[Link] para un flujo estacionario en tal sistema de dos acuíferos puede ser
escrito como la suma de flujos confinados y la solución para flujo semiconfinado.

Primera aplicación : Flujo a un rio en un sistema de dos aquíferos .Recarga es

insignificante .El rio es estrecho y en contacto superior con el acuífero superior .La
condición de contorno es h0 x=0=hr
Una cantidad flujo σ (L2/T) ,la descarga σ es obtenida de continuidad de flujo.

Por lo que habrá 4 funciones dos para acuifero 0 y 2 para acuifero 1 ,asimismo el Q
tendrá 4 funciones

Se observa que en función de heads ,hay 4 funciones ,el punto de truncamiento es

cuando Q=0

AREA DE RECHARGA ENTRE DOS RIOS EN UN SISTEMA DE DOS ACUÍFEROS

Flujo entre dos rios con área de recarga N en un

El rio atraviesa el acuifero superior pero no el

Acuifero inferior .Ahora ambos acuíferos se dividen

De una leaky layer .El origen es en el medio de los

Dos rios y la distancia entre los rios es 2L

Ten en cuenta que el acuifero inferior esta

delimitado por capas impermeables lo que da
a entender la derivada es 0
Cuando se hace un único modelo usando la formula se sección 1.2 varia el grafico de
head y de Q
Flujo no confinado unidimensional constante con espesor saturado
variable