Flujo de aguas subterráneas

Juan Camilo Tabares Prada

Escuela ambiental

Universidad de Antioquia

Resumen
El siguiente trabajo tiene como objetivo alcanzar una primera inducción acerca de los flujos de aguas subterráneas y su
relación con la mecánica de fluidos. Primeramente, se hará una recopilación acerca de cuáles son los modelos utilizados para
representar estos flujos y las ecuaciones utilizadas para su construcción, trataremos los métodos más utilizados en el campo de
la ingeniería, además conoceremos en que actividades se usan mas estos tipos de modelación y cuáles son sus aportes para la
ingeniería y el desarrollo de las comunidades.

Palabras claves: Aguas subterráneas, Ley de Darcy, Hidrogeología, MODFLOW.

1. Introducción sitio, identificando los parámetros mas importantes

con el fin de que nuestro modelo se aproxime lo mas
Como primer acercamiento al campo de las aguas posible a la realidad, MODFLOW es un modelo que
subterráneas es importante conocer los hechos responde bien a la hora de la recopilación de datos,
físicos que dieron a evidenciar el comportamiento de manejo de los parámetros y la ejecución de las
los fluidos a partir de medios porosos. Henry Darcy ecuaciones matemáticas que permiten analizar los
ingeniero de caminos Frances fue el que descubrió flujos en tres dimensiones.
mediante un experimento con un aparato que hoy se
le conoce como “permeámetro”, demostró como un
fluido pasa por un material poroso de un extremo a
otro cambiado su carga potencial hidráulica. Esto 2. Marco conceptual
finalmente se tradujo en la ley de Darcy donde es Respondiendo a la pregunta.
posible determinar el caudal conociendo el área de la ¿Consultar cuáles son los modelos utilizados para
sección del fluido, el gradiente hidráulico y una representar estos fenómenos y cuáles son las
constante k que es conocida como la conductividad ecuaciones y suposiciones que consideran?
hidráulica, definida como la capacidad de un
material a dejar pasar un fluido a través de el sin Ecuación de Bernoulli
alterar su estructura. Por otro lado la ecuación de
Bernoulli aplicada a las aguas subterráneas, es decir, v2 P (1)
considerando la velocidad de esta ecuación como + + Zg=constante
2 p
despreciable y despejando la posición, podemos
decir que el agua se mueve de puntos de mayor
potencial hidráulico a zonas de menor potencial Dividiendo por la gravedad tenemos y denotando
hidráulico, así entendiendo un poco los sistemas de γ= pg, se tiene la ecuación de Bernoulli para el flujo
acuíferos que son encargados de almacenar y de aguas subterráneas.
transmitir grandes cantidades de agua. La
modelación de aguas subterráneas no es mas que la v2 P (2)
consolidación de la información hidrológica del h= + +Z
2g γ
 del medio, la cual es una característica de cada
donde h es, la altura piezométrica (carga hidráulica)
material.
y lo denotamos por
h=he + hp + hv. La ley de Darcy se basa en las siguientes hipótesis
que condicionan la validez de su ley
donde he = z( altura de elevación), hp = P (altura de
γ - Medio continuo, es decir que los poros
vacíos estén intercomunicados.
presión), hv= v2 .
(altura de velocidad ) - Medio isotrópico.
2g
- Medio homogéneo.
Despreciando la velocidad, debido a que los flujos - Flujo del agua en régimen laminar.
en acuíferos son lentos, se tiene que
(Núñez, 2018)
P (3) Ecuación de continuidad.
h= + Z
γ Flujo estacionario (Régimen permanente)
Donde: Salidas – Entradas = 0
h conocido como altura piezométrica que describe el
- Nivel piezométrico del acuífero es
comportamiento del agua, el agua se mueve desde
puntos de mayor a menor potencial hidráulico. constante.
(Núñez, 2018). - En la ecuación no aparece el tiempo
(t).

Flujo transitorio (Régimen no permanente)

- Nivel piezométrico del acuífero
variable
- Pozo bombeando
- En las ecuaciones aparece el tiempo.
(Nuñez,2018)
Ley de Darcy

Darcy encontró que la cantidad de agua fluye a Modelo para la representación de flujos
través de un medio poroso por unidad de tiempo, es subterráneos
decir el caudal o descarga Q, es proporcional a la
sección transversal A y a la perdida de carga MODFLOW
(diferencia entre las alturas del agua de los El programa al que se hace referencia es el modelo
piezómetros Δh= h1-h2, e inversamente
proporcional a la longitud del material poroso hidrogeológico con mayor utilización en el mundo
(Sánchez, 2016). científico. Ofrece una simulación versátil, ya sea con
flujo estacionario o transitorio, en escenarios
Δh (4)
Q=K . A . saturados confinados, no confinados o ambos.
L
Resuelve las ecuaciones diferenciales parciales
Donde k es la constante de proporcionalidad,
usando el método de diferencias finitas, donde cada
llamada conductividad hidráulica o permeabilidad
celda de Mod-Flow es un volumen unitario. En
general la ecuación de flujo de aguas subterráneas -Primeramente, el recurso hídrico es fundamental
que emplea Mod-Flow es la siguiente (Rodríguez. para a vida, es necesario para realizar diferentes
2010): actividades como la ganadería, agricultura y
d dh d dh d dh dh abastecimiento de agua potable para consumo
( Kx )+ ( Ky )+ (Kz )=Ss .
dx dx dy dy dz dz dt humano. Además, tiene un valor económico, social y
ambiental en todos los usos que se le destina.
(5)
-El conocimiento de acuíferos con fines de
La ecuación es conocida como “Ecuación general planificación y determinación de criterios a la hora
del flujo de agua subterránea” de determinar concesiones de aguas subterráneas
para la población. (Hernandez,2018)
Características importantes - Monitoreo de cuencas con presencia de
-Representa la salida de agua subterránea en (Plano contaminantes
X, Plano Y, Plano Z) En un medio heterogéneo y - caracterización de flujos de agua a partir de su
anisotrópico. salinidad para diferentes actividades. (Fagundo, 2016)
-Admite representar flujo en estado permanente y
transitorio en acuíferos del tipo no confinados,
confinados y modelos confinantes.
-Admite simular los efectos de pozos, recargas, ríos, Resolviendo a la pregunta
drenes y evapotranspiración. -¿Como los temas desarrollados durante el curso se
-Muestra una distribución por módulos en el relacionan con estos modelos?
programa y en la documentación que facilita una Como pudimos ver en la primera pregunta, el modelo
de Mod-Flow se basa en el principio de lay de
mejor comprensión del modelo.
continuidad, la ecuación del balance de masas y la ley
´-Usa diferencias finitas. (Yagual, 2016) de Darcy, teoría vista durante el curso.
A continuación, se mostrara la demostración de la
Fortalezas del programa
ecuación general del flujo de agua subterránea, ecuación
-El esquema numérico de resolución por diferencias utilizada para la modelación Mod-Flow.(Núñez, 2018)
finitas, entrega soluciones aceptables a los
La ecuación general del balance de masas involucra
problemas modelados, a través de una metodología consideraciones de entrada y salida de flujos y cambios
en el almacenamiento, es decir.
numéricamente sencilla.
-Permite su aplicación a una gran variedad de Entradas- Salidas = Variación de almacenamiento

condiciones del mundo real.

Analizando un cubo de continuo, con un volumen de
-La representación gráfica de la modelación, facilita control denotado por la línea color azul.

el ingreso de datos y la visualización de resultados.

(Yagual, 2016)

Resolviendo a la pregunta
¿Cuáles son las aplicaciones desde el punto de vista
de la ingeniería?
 1 d (∆ z )
α=
∆ z dP
(12)

d(∆z) = α∆zdP (13)

Variación de la porosidad del terreno con la presión

d (m) = (1 − m) (d (∆z))/ ∆z (14)

Relacionando estas ecuaciones obtenemos.

dρ = βρ(ρgd(h)) (15)
La masa de fluido que entra por una cara es igual a la
densidad por el caudal especifico por el área de paso, es lo cual muestra la relación entre las variaciones de la
decir: densidad y las variaciones de h.
M =¿ ρqA
d(∆z) = α∆z(ρgd(h)) (16)
(6)
La masa de agua en el volumen de control esta dada por
lo cual muestra la relación entre las variaciones de
el producto del volumen por la densidad. ∆z y las variaciones de h.
d [m] = (1 − m)αρgd(h) (17)
M =¿ ρV (7)
Sin embargo, el volumen del material solido dentro del lo cual muestra la relación entre las variaciones de la
volumen de control, las entradas y las salidas para las porosidad con la altura(o potencial piezométrico) h.
caras del volumen de control esta dado por.
Estas ´últimas tres expresiones (15), (16) y (17)
∆x∆y∆z − m∆x∆y∆z ; (1 − m)∆x∆y∆z (8) relacionan las variaciones de presión, densidad y
porosidad con las variaciones del potencial
∂ ∂ ∂ (8.1) piezométrico h. Esto nos permitirá relacionar los
−( ( ρqx)+ ( ρqy)+ (ρqz))¿ flujos y almacenamientos
∂x ∂y ∂z
Ahora retornamos a la expresión (9), que determina
Derivando tenemos que la variación temporal de la la variación temporal de la masa M respecto al
masa M respecto al tiempo seria. tiempo, la cual lo vamos a derivar aplicando la
derivada de un producto y a la vez sustituiremos las
∂M ∂ (9) ecuaciones (15), (16) y (17).
= ¿
∂t ∂t
Resolviendo tenemos finalmente

∂M ∂h
Parámetros que pueden provocar variaciones de la =(αρg+ mβρg)ρ ∆ x ∆ y ∆ z
masa contenida en el volumen de control. ∂t ∂t
Tenemos lo siguiente:
La compresibilidad del agua. (18)
1 dρ
β= Reemplazando la ecuación (8.1) y (18) en el balance
ρ dP
general de masas.
(10)
Entradas– Salidas = Variación de
dρ = βρdP (11) almacenamiento

La compresibilidad del terreno.

Simplificando tenemos que primeros ha sido conocerla la diferencia que
∂ qx ∂ qy ∂ qz ∂h produciría en el futuro la aplicación de los
−( + + )=ρg(α +mβ)
∂x ∂y ∂z ∂t condicionantes establecidos en el plan de
(19) explotación de la masa de agua subterránea de
acuífero. El ultimo escenario trata de representar el
La ley de Darcy nos define el flujo especifico “q” de posible efecto que el cambio climático tendría sobre
acuerdo a las tres direcciones de X,Y,Z. el sistema.
Datos de interés: el acuífero cuenta con 98 unidades
dh ; dh ; dh
q x=−Kx q y =−Ky q x=−Kx de demanda agrícola (1101 pozos) y 92 unidades de
dx dy dy demanda urbana.

Llamamos ∂ h , el cual es el
Ss= ρg(α + mβ)
∂t
“coeficiente de almacenamiento especifico”. Conclusiones al utilizar el modelo MODFLOW.
Reemplazando en la ecuación (19) obtenemos
finalmente “la ecuación general del flujo de agua  El funcionamiento del sistema en un
subterránea.” periodo de 35 años esta regido por el
aumento progresivo del volumen extraído,
d dh d dh d dh dh debido al cambio de usos experimentado y
( Kx )+ ( Ky )+ (Kz )=Ss .
dx dx dy dy dz dz dt a la expansión de las superficies regadas, lo
(20) que provoca una disminución del flujo que
el sistema drena hacia los ríos y un
descenso de las alturas piezométricas en
algunas zonas. Asimismo, a lo largo de este
periodo las precipitaciones sufren una
Casos de estudio ligera disminución con el tiempo. Se
1. confirma que la masa de agua subterránea
El siguiente caso de estudio hace referencia a una de la zona se encuentra en mal estado
modelación matemática de un acuífero localizado en cuantitativo.
la provincia de valencia, con el fin de analizar el  En la relación con la interacción rio-
estudio del plan de explotación de la masa de agua, acuífero, se observa que el volumen que el
que últimamente venía mostrando unas tendencias sistema drena al rio Magro disminuye
no sostenibles en la cuenca, como la disminución en progresivamente, lo que puede constituir
los aportes al embalse de la zona y de las aguas una de las posibles causas de la
subterráneas que son utilizadas para el riego. Es por disminución de los niveles del embalse. Sin
esto que se considera el plan hidrológico de la embargo, la disminución de la precipitación
demarcación hidrográfica del Juca, que considera que se ha producido también es una causa
que el agua de la dicha zona es declarada en mal incluso mas influyente de la reducción de
estado cuantitativo. El plan hidrológico considera la las aportaciones al rio.
elaboración de un modelo matemático (Con código  En cuanto a los escenarios simulados, el
MODFLOW) que permita un mejor conocimiento de horizonte temporal considerado ha sido el
la información que ahora se tiene. (Fernandez,2017) año hidrológico 2030/2031. La aplicación
Por último, con el modelo calibrado y con el fin de del Plan de Explotación de la masa de agua
predecir el comportamiento del sistema ante subterránea Requena-Utiel podría constituir
diferentes situaciones o hipótesis, se han simulado una solución para los usuarios de la zona.
tres diferentes escenarios. La finalidad de los dos No obstante, el estado del sistema
 empeoraría notablemente, alcanzando unas acuífero acuitardo y descendentes desde el
pérdidas de 35 hm3 anuales. Dicha acuífero libre.
situación se podría agravar bajo los efectos  Se identifico, para el acuífero confinado,
del cambio climático. (Fernandez,2017) flujos desde la frontera con la unidad
metamórfica del paleozoico lo cual va en
2. dirección a confirmar la recarga lateral
El siguiente caso hace referencia a la construcción desde la zonas altas y alejadas de la zona de
de una herramienta de simulación, los modelos estudio.
numéricos en hidrogeología ofrecen un camino para
avanzar en el entendimiento de sistemas de 3.
acuíferos. En esta ocasión se hará la modelación El siguiente caso de referencia habla sobre el
numérica con el propósito de ayudar, a interpretar la modelamiento del flujo subterráneo de la
información y la dinámica de flujo subterráneo del microcuenca Macari-Melgar Perú en donde se quiere
Bajo Cauca Antioqueño. conocer básicamente el funcionamiento de esta al ser
El sistema de acuíferos del Bajo Cauca Antioqueño, utilizada como fuente para abastecimiento de
esta conformado por tres unidades hidrogeológicas: actividades agrícolas, ganaderas y crecimiento
Un acuífero libre, un acuitardo y un acuífero poblacional que de por si estos factores contribuyen
confinado. a la intensificación del estrés en los recursos
El acuífero libre está formado por un conjunto hídricos.
estratigráfico que consta de una delgada capa de La microcuenca Macari-Melgar presenta zonas
suelo, depósitos aluviales recientes y el saprofito semiáridas, disminución de cauces, lugares de
poco consolidado de las rocas sedimentarias del empinamiento, poca vegetación y uso
terciario del Miembro superior de la formación indiscriminado del recurso hídrico, con
consecuencias de disminución de fuentes de aguas a
cerrito; con espesores que varían entre 10 y 90
las comunidades aledañas en el distrito de Macari.
metros , su extensión abarca toda la planicie, registra Ocurriendo desabastecimiento en meses de estiaje.
valores de conductividad hidráulica entre 1 y 5 [9]
m/día. El acuitardo muestra una conductividad Para la modelación es imprescindible el uso de las
hidráulica de 0.1m/día los espesores varían entre 100 variables como son el DEM, ríos y pozos
y cero metros en sitios donde la unidad desaparece al Los resultados muestran una recarga de 1.4971m3/s
norte cerca de los afloramientos del acuífero de infiltración en el subsuelo en sus diferentes capas.
confinado. El acuífero confinado, miembro inferior En la descarga tenemos el drenaje por los ríos que
de Cerrito esta formado por intercalaciones de conforma la microcuenca con un total de
areniscas conglomeráticas, sus espesores varían 0.9474m3/s, y descarga por evapotranspiración que
entre 10 y mas de 100 metros; es unidad poco es la absorción de la del subsuelo 0.5498m3/s.
explorada y explotada, podría constituir una En conclusión, se observó mediante el balance
importante reserva de agua subterránea para la hídrico la distribución del agua en el subsuelo y la
región. (Betancur,2009) gran mayoría se descarga por drenaje en los ríos que
conforma la microcuenca.
Conclusiones al utilizar el modelo matemático. El balance hídrico nos muestra que las actividades
humanas, generan un impacto negativo al recurso
 En relación con el acuífero acuitardo y hidrogeológico. Con el pasar del tiempo este recurso
confinado, si bien se modelaron como puede llegar a disminuirse drásticamente por que las
medios homogéneos, asignándoles valores autoridades locales y la población no toman en
de conductividad teóricos. Se observo la conciencia este aspecto que es indispensable para la
ocurrencia de flujos verticales y sostenibilidad del recurso agua y todo el ecosistema
ascendentes desde el acuífero confinado al de la microcuenca Macari-Melgar. (Ahuanlla,2020)
 Referencias
[1] BR. Joel Núñez Mejía(2018).Aplicación del método
de diferencias finitas para el flujo de un contaminante
en un acuífero confinado homogéneo e isotrópico.
Universidad Nacional Pedro Ruiz Gallo. (Perú)
[2] Javier Sánchez San Román. (2016) Ley de Darcy.
Conductividad hidráulica. Universidad Salamanca
(España)
[3] Miranda Rodríguez y Francisco Ariel. Modelación
del comportamiento hidrodinámico del agua
subterránea de la zona comprendida entre Carbó,
Pesqueira y Zamora, Sonora. Universidad de sonora
[4] Omar Yagual Muñoz y Teresa María Reyna (2016).
Modelación de flujo subterráneo. Aplicación en la
cuenca del milagro y análisis de pozo en subcuenca del
Rio Suquía. Universidad Nacional de Córdoba.
(Ecuador).
[5] Sebastián Hernando Hernández Dueñas (2018)
3. Conclusiones. Aplicación de visual MODFLOW y sistemas de flujo de
La modelación de aguas subterráneas es muy Toth, como soporte al modelo hidrogeológico
importante en la actualidad debido a que nos permite conceptual de la secretaria de ambiente, en el bloque
evaluar el recurso hídrico actual y realizar proyecciones centro de la ciudad Bogotá. Universidad de la Salle
de almacenamiento con el fin de poder tomar decisiones (Bogotá).
en bienestar de la comunidad y del ambiente. Como ya [6] Juan R. Fagundo, Margarita M. Alconada, J. Joel
sabemos, el agua es el recurso mas importante del Carillo Rivera y Patricia González Hernández. (2014)
planeta, es la base de la vida y por lo tanto debemos caracterización de los flujos de agua subterránea a partir
hacer un uso racional de ella. El estudio de acuíferos, los de su salinidad. Universidad la Habana (Cuba).
contaminantes y los métodos de extracción, son carácter [7] Irati Fernández Urizar (2017). Modelación
de estudio para los ingenieros de hoy en día. La Matemática del flujo de agua subterránea del acuífero
modelación hidrológica es una herramienta eficaz si de requena- utiel. Estudio del plan de explotación de la
queremos realizar una aproximación a los masa de agua. (Valencia)
comportamientos de los flujos, sin embargo, la labor del [8] Teresita Betancur y Carlos Palacio (2009). La
ingeniero esta en poder interpretar esta información y modelación numérica como herramienta para la
transmitirla de tal manera, que la sociedad este al exploración hidrogeológica y construcción de modelos
pendiente, ya que en la actualidad tenemos algunos conceptuales. Universidad de Antioquia (Colombia).
problemas con el uso de las aguas subterráneas como [9] José Luis Ahuanlla Huaccha y Renny Daniel Diaz
por ejemplo el desabastecimiento, la contaminación, la Aguilar (2020) Aplicación de MODFLOW 6 con Model
extracción excesiva, entre otras, es por eso que el Mouse 4 para el modelamiento del flujo subterráneo de
ingeniero debe promover el buen uso para que se haga la microcuenca Macari-Melgar. Universidad Peruana
una correcto aprovechamiento de tal forma de no alterar Unión (Perú).
significativamente el ciclo hidrológico. El trabajo en
conjunto con la población y los entes encargados del
manejo territorial y de los recursos naturales es de vital
importancia pues se podría analizar por población como
es el estado de las aguas subterráneas para su
aprovechamiento y poder promover actividades que
generen un avance sostenible para las regiones, pues es
evidente que, en algunos casos, se hace un uso indebido
de las aguas subterráneas en zonas donde el fluido
escasea y puede ser utilizado para cosas mas
importantes.