Diferencia entre revisiones de «Nivel piezométrico G5»

Revisión del 13:25 14 may 2014

Definimos nivel piezométrico como la altura que alcanzaría el agua al realizar un sondeo en un punto de un acuífero confinado. Este valor depende de la presión a la que esté el propio acuifero.

Si construimos sobre el acuífero confinado un pozo circular de radio [math]\rho _{0}[/math], el nivel piezométrico varía. Para que el problema sea más sencillo utilizaremos coordenadas cilíndricas. Para poder conocer la variación del nivel piezométrico nos apoyaremos en la ecuación de la conservación de la masa y la ley de Darcy:

[math] S ·\frac{ \partial h }{ \partial t } + div q = 0[/math]

[math] q = - k ·\nabla h [/math]

La ley de Darcy establece que "el flujo de agua q a través de un medio poroso es proporcional a la diferencia de presión, que a su vez se puede escribir en términos del gradiente del nivel piezométrico en cada punto". La constante K se deduce experimentalmente para cada material y se conoce como la conductividad hidráulica o permeabilidad. La constante S en la ley de conservación de la masa se conoce como almacenamiento específico y se interpreta como la cantidad de agua que libera el acuífero al descender el nivel piezométrico en una unidad, por unidad de volumen. Combinando las ecuaciones de conservación de la masa con la ley de Darcy, obtenemos la ecuación:

[math] \frac{ \partial h }{ \partial t } - D · \Delta h = 0, \quad \rho \gt \rho _{0} \quad θ\in (0,2\pi ) \quad t\gt0 \quad (1) [/math]

Donde [math] D= \frac{k}{s}[/math] se conoce como difusividad hidráulica.

Pozo

1 Obtención del Laplaciano y ecuación diferencial en polares

El Laplaciano es la divergencia del gradiente. Al aplicarle ambos operadores a nuestra función, nos da que:

[math]\Delta h(\rho,\theta)[/math] = [math] (\frac{\partial ^2 h}{\partial \rho^2}+ \frac{1}{\rho}·\frac{\partial h}{\partial \rho}+\frac{\partial^2 h}{\partial \theta^2})= 0[/math]

y aplicándolo a la fórmula (1) nos quedaría:

[math] \frac{ \partial h }{ \partial t } - D·(\frac{\partial ^2 h}{\partial \rho^2}+ \frac{1}{\rho}·\frac{\partial h}{\partial \rho}+\frac{\partial^2 h}{\partial \theta^2})= 0[/math], [math]\rho \gt[/math] [math]\rho _{0} [/math]

Finalmente, al decirnos que h debe depender solamente de [math]\rho[/math],es decir h=h([math]\rho[/math]), se deduce que [math]h _{θ}[/math]=0 y por lo tanto [math]h _{θθ}[/math]=0, por lo que nuestra ecuación final será:

[math]\frac{\partial h}{\partial t}- D·(\frac{\partial ^2 h}{\partial \rho^2}+ \frac{1}{\rho}·\frac{\partial h}{\partial \rho})=0[/math]

2 Sistema completo de ecuaciones

Las condiciones de frontera serán:

[math]h(\rho _{0},t) = h _{\rho0} \quad \quad h(20,t) = h _{0}[/math]

Para que el problema tenga una única solución, nos falta una condición inicial del tipo [math]h(\rho,0) = h _{i}(\rho)[/math]. El problema quedará:

\[\left\{\begin{matrix}\ \frac{\partial h}{\partial t}- D·(\frac{\partial ^2 h}{\partial \rho^2}+ \frac{1}{\rho}·\frac{\partial h}{\partial \rho})=0 , & \\ h(\rho _{0},t)=h _{\rho} \quad h(20,t)=h _{0}, & h(20,t)=h _{0} \\ h(\rho,0)=h _{i}(\rho) , & \end{matrix}\right.\]

3 Resolución del problema por diferencias finitas y método del trapecio

4 Dibujo del comportamiento del nivel piezométrico

5 Método de Euler

6 Nivel piezométrico en intervalos de tiempo grandes

7 Estado estacionario

[math]\frac{\partial h}{\partial t}=0 [/math] [math]D·(\frac{\partial ^2 h}{\partial \rho^2}+ \frac{1}{\rho}·\frac{\partial h}{\partial \rho})=0[/math]

8 Capacidad de recuperación en acuífero

9 Método de Fourier