Tfe Men Aba

Abatimiento del nivel piezométrico en acuíferos confinados
con pozos de bombeo. Comparación entre la solución

analítica y numérica mediante Modflow.
Alumno: PABLO JOSÉ MÉNDEZ CÁNOVAS

Director: IVAN ALHAMA MANTECA
AGRADECIMIENTOS
Este Trabajo Final de Grado de Ingeniería Civil, se centra en la rama de geotecnia y

suelos. Este área despertó un cierto interés en mi gracias al profesor de esta asignatura y
director de este trabajo: Iván Alhama Manteca, quien hace que la geotecnia sea mucho
más interesante.
Quería aprovechar estas palabras para agradecerle su ayuda, y todo el interés que ha
mostrado durante la realización del mismo, donde ha sido un excelente tutor
prestándose a ayudarme en la resolución de cualquier problema que me surgiera.
Así mismo, me gustaría tener en consideración a mis padres José Salvador y

Encarnación, a mi hermana Belén, a mi novia María, y a todos aquellos amigos que
siempre han estado ahí para apoyarme y animarme, tanto en los buenos como en los
malos momentos, y no solo durante la realización de este trabajo sino durante toda la
carrera.
Ya que este es un trabajo que pone fin a mis estudios de Ingeniería Civil, también quería
agradecer al resto de profesores, su dedicación y todo lo que me han enseñado durante
estos años.
Gracias a todos.
ÍNDICE
1 INTRODUCCIÓN ................................................................................................... 1
2 OBJETIVOS Y ALCANCE ...................................................................................... 2
3 FUNDAMENTOS TEÓRICOS ................................................................................. 3
3.1 El agua en el medio poroso. Definición de acuífero y parámetros fundamentales....... 3
3.2 Piezometría ...................................................................................................... 7
3.3 La ecuación de Darcy ...................................................................................... 11
3.4 La ecuación general de flujo: Ecuación de Laplace y ecuación de Poisson.............. 13
3.5 Resolución analítica para el diseño de abatimientos mediante ensayos de bombeo... 17
3.5.1 Ensayos de bombeo en régimen permanente. Deducción de las ecuaciones...... 18
3.5.2 Formulación analítica en el diseño de abatimientos de acuíferos confinados.

Flujo proveniente de una fuente circular; caso de pozo simple y pozos múltiples............ 21
3.6 Fuerzas de filtración. El fenómeno del sifonamiento ............................................ 29
4 EL PROGRAMA MODFLOW ................................................................................ 35
4.1 Descripción general ......................................................................................... 35
4.2 Interfaz de comunicación con modflow: Visual Modflow Classic.......................... 37
5 APLICACIONES ................................................................................................... 39
5.1 Solución analítica ............................................................................................ 40
5.1.1 Comprobación de sifonamiento.................................................................. 43
5.2 Solución numérica........................................................................................... 46
5.2.1 Creación del modelo (‘Input’) .................................................................... 49
5.2.2 Ejecución del modelo (‘Run’) .................................................................... 56
5.2.3 Resultados obtenidos (‘Output’) ................................................................. 58
5.3 Comparación entre la solución analítica y numérica y discusión de resultados ........ 67
6 CONCLUSIONES ................................................................................................. 71
7 REFERENCIAS..................................................................................................... 72
Abatimiento del nivel piezométrico en acuíferos confinados mediante pozos Autor: Pablo José Méndez Cánovas
de bombeo. Comparación entre solución analítica y numérica mediante Modflow Grado en Ingeniería Civil
1 INTRODUCCIÓN
En este trabajo se aborda un problema de abatimiento del nivel piezométrico de una

excavación en un acuífero confinado mediante el empleo de métodos numéricos. La
estructura del documento describe, en primer lugar, los fundamentos teóricos del flujo de
agua en medios porosos, los distintos parámetros hidrogeológicos que se deben tener en
consideración y la explicación de la solución analítica del problema a resolver.
Posteriormente se hace una descripción del software utilizado para la solución numérica:
Visual Modflow Classic. Finalmente se plantea un problema y se calcula la solución
analítica y numérica, comentando y comparando los resultados obtenidos.
En este trabajo se afronta un campo fundamental en ingeniería del terreno, la presencia de

agua en el suelo, desde el punto de vista de los métodos numéricos. La tecnología de los
softwares de métodos numéricos ha alcanzado una gran precisión y exactitud, lo que hace
que se hayan convertido hoy día en una herramienta que facilita la labor del experto en
geotecnia. En concreto, el empleo de estos métodos numéricos para la resolución de
problemas de flujo en medios porosos, es en ocasiones exigido en los pliegos de
condiciones técnicas de obras.
1
2 OBJETIVOS Y ALCANCE
El objetivo principal que tiene la realización de este proyecto será la familiarización con el
manejo del software Visual Modflow Classic. Este programa permite la resolución
numérica de problemas de flujo de agua subterránea, o más exactamente, la evolución de
los niveles debido a condiciones especificadas: bombeos, ríos, drenes, infiltración de
precipitaciones, etc…
Además, como objetivo secundario, se adquirirán las técnicas, ya sean analíticas o

numéricas, para el diseño de un abatimiento del nivel piezométrico en acuíferos confinados
mediante pozos de bombeo y para estudiar el fenómeno del sifonamiento, algo que resulta
especialmente interesante en abatimientos para excavaciones. La comparación de los
resultados analíticos y numéricos, aportará una visión profunda del problema y de las
técnicas de resolución.
Con todo esto, parece necesario hacer una revisión de los aspectos teóricos del flujo de
agua en medios porosos, para poder ponerse en disposición de controlar y entender los
parámetros que se introducen en el input e interpretar los resultados del output.
Pese a la simplificación del escenario planteado, se adquirirán los conocimientos básicos

necesarios para abordar problemas más complejos, como pueda ser el de acuíferos
multicapa
2
3 FUNDAMENTOS TEÓRICOS
3.1 El agua en el medio poroso. Definición de acuífero y parámetros fundamentales
Un acuífero es una formación geológica que almacena agua y que es capaz de transmitirla
de manera que puede ser aprovechada como recurso. El agua de superficie se infiltra a
través del suelo hasta encontrarse con una capa impermeable; entonces se acumula y llena
los poros y fisuras de las rocas, dando lugar a una zona de saturación. El agua subterránea
puede brotar de modo natural (manantiales o fuentes) o ser extraída mediante pozos u otros
elementos de captación. Las características fundamentales de un acuífero son la baja
velocidad en el movimiento del agua subterránea, los grandes volúmenes de reservas y el
gran tiempo de renovación del agua en el sistema. Según el comportamiento hidráulico del
agua, los acuíferos se pueden clasificar en:
- Acuífero libre: es un acuífero no limitado en la parte superior por un techo impermeable,

de manera que existe un nivel freático a una cierta profundidad. Cuando una perforación
alcanza este nivel aparece un cierto volumen de agua libre. Estos acuíferos pueden
recargarse desde la superficie mediante una fracción de la pluviometría o excedentes de
riego.
- Acuífero confinado: es un acuífero limitado en la parte superior por un techo

impermeable, de manera que el agua que el agua se encuentra a presión superior a la
atmosférica en todos los puntos. Cuando una perforación llega al techo del acuífero el agua
sube por el tubo hasta que se estabiliza a una profundidad correspondiente a su nivel
piezométrico. La recarga de estos acuíferos se produce lateralmente.
- Acuífero semiconfinado: es un acuífero de tipo confinado que se diferencia de este en que

se encuentra limitado en la parte superior por una formación geológica muy poco
permeable pero que es capaz de producir recarga al acuífero.
3
Figura 3.1. Tipos de acuíferos

Fuente: http://www.artinaid.com/2013/04/que-es-un-acuifero/
La capacidad de almacenar y transmitir el agua por parte de un acuífero debe ser

cuantificada para evaluar las formaciones desde el punto de vista hidrológico. Por ello, se
asocian a las formaciones cuatro parámetros básicos en cuanto a su comportamiento
respecto al agua que puedan contener y transmitir. Estos cuatro parámetros
hidrogeológicos característicos de las formaciones geológicas son los siguientes:
- Porosidad: Es la relación entre el volumen de huecos y el volumen total de una roca. Es

un parámetro adimensional y depende únicamente de la constitución de la roca o suelo, es
decir, de su textura característica, sin que intervenga la forma geométrica ni la potencia de
la formación o su mecanismo de funcionamiento hidráulico en la naturaleza. Según el tipo
de formación, los poros pueden ser debidos a espacios intergranulares, en las formaciones
detríticas, o a grietas y fisuras, en el caso de rocas fisuradas o cársticas; el concepto de
porosidad puede asociarse a unas y otras. La porosidad de una formación puede estar
ligada únicamente a la textura de la misma o depender además de las características del
fluido que se mueve en su interior. La primera de ellas es la porosidad total n, referida al
volumen total de poros, independientemente de que el fluido tenga o no la capacidad de
circular entre ellos
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑙𝑜𝑠 𝑝𝑜𝑟𝑜𝑠

𝑛=
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙
4
- La porosidad eficaz ne: también denominada porosidad cinemática, referida al volumen

de poros conectados por los que el transporte de fluido es posible, y queda ligada no solo a
la textura de la formación, sino también a las características del fluido:
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑝𝑜𝑟𝑜𝑠 𝑐𝑜𝑛𝑒𝑐𝑡𝑎𝑑𝑜𝑠

𝑛𝑒 =
𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑡𝑜𝑡𝑎𝑙
- Coeficiente de almacenamiento: Este coeficiente representa la capacidad para liberar agua

de un acuífero. Se define como el volumen de agua que es capaz de liberar un prisma de
base unitaria y de altura la del acuífero, cuando desciende 1 m el nivel piezométrico. Es, en
consecuencia, un parámetro físico referido al volumen cedido por unidad de volumen de
acuífero y, al igual que la porosidad, es adimensional. La forma en que las formaciones
geológicas permeables liberan agua es distinta según su situación en la naturaleza y el
estado de sus niveles piezométricos. Así, cabe distinguir formaciones o acuíferos en
funcionamiento confinado o acuíferos en funcionamiento libre.
- Permeabilidad: La permeabilidad es el parámetro que permite evaluar la capacidad de

transmitir agua de una formación en función de la textura de la misma, sin relacionarla con
su estructura o forma geométrica. Dentro de este concepto pueden definirse otros dos:
La permeabilidad ligada tanto a las características texturales del medio físico como del
fluido que transmite, denominada conductividad hidráulica o representada por el
coeficiente de permeabilidad k, y la denominada permeabilidad intrínseca K, que es sólo
dependiente de las características internas del medio poroso.
La permeabilidad efectiva k se define como el caudal que es capaz de atravesar una sección
unitaria del acuífero, normal al flujo, bajo un gradiente piezométrico unitario. Depende, en
consecuencia, tanto de las características del medio físico (acuífero), como de las del fluido
que lo atraviesa (agua). Por razones de su definición (caudal por unidad de superficie), las
dimensiones de la permeabilidad efectiva pueden establecerse:
𝐿3 𝑇¯1
[𝑘] = = 𝐿𝑇¯¹
𝐿2
5
Se usan cm/s y m/día como unidades más habituales (m/día suele utilizarse en
hidrogeología con carácter prácticamente general), y cm/s en ingeniería. La permeabilidad
intrínseca y la permeabilidad efectiva se relacionan únicamente por los parámetros capaces
de definir las características del fluido.
𝛾 𝑔
𝑘=𝐾 𝜇
O bien 𝑘=𝐾 𝜂
donde: k = permeabilidad efectiva [𝐿𝑇¯1 ]

K = permeabilidad intrínseca
𝛾 = peso específico del agua [𝑀𝐿¯2 𝑇¯2 ]
𝜇 = viscosidad dinámica del fluido [𝑀𝐿¯1 𝑇¯1 ]
𝜂 = viscosidad cinemática [𝐿²𝑇¯1 ]
g = aceleración de la gravedad [𝐿𝑇¯2 ]
y en consecuencia:
𝜇 𝐿𝑇¯1 ∗ 𝑀𝐿¯¹𝑇¯¹
[𝐾] = �𝑘 � = = 𝐿²
𝛾 𝑀𝐿¯²𝑇¯²
No ha sido posible encontrar una relación capaz de definir, con carácter generalista, la
permeabilidad intrínseca en función del medio, a pesar de los esfuerzos de diversos
investigadores. En cualquier caso, no influye solamente el tamaño de las partículas que
conforma el medio, sino también su forma y superficie, no admitiendo generalización estas
dos últimas características. Diversos autores de esta temática han establecido las siguientes
relaciones:
K = c * d² Hazen
K = c * m * d² Slichter
𝑚−0.13
K = λ * ((1−𝑚)⅓)² * d² Terzaghi
La variable d es el denominado diámetro eficaz, que, para el caso más general es el d10 de
la curva granulométrica de la muestra; es decir, el tamaño de abertura de tamiz que permite
6
el paso del 10 % en peso de la muestra de material granular. Los coeficientes c, m, n y λ

dependen de las características geométricas y de fricción de los granos, pudiendo por ello
admitirse que para cualquier caso:
K = cte * d² y por tanto: K = cte * d10²
Para un mismo fluido a una misma temperatura. La permeabilidad, al contrario que ocurría
con la porosidad, sí depende del tamaño del grano; tal y como se ha visto, el valor d10
condiciona el valor de la permeabilidad.
En definitiva, la conductividad hidráulica o permeabilidad efectiva k nos da una idea de la

relación existente entre el caudal por unidad de sección y el gradiente hidráulico. A mayor
conductividad hidráulica mayor caudal, o menor gradiente hidráulico se precisa para un
mismo caudal. Guarda relación con otros parámetros que intervienen en otros procesos
físicos, como la conductividad eléctrica (inverso de la resistividad) y depende de las
características del medio poroso y del fluido. Mientras que la permeabilidad intrínseca nos
da una idea de las características del terreno, y al medirse en unidades de superficie y tener
valores tan bajos, podemos considerarla como “una medida de la superficie del poro”.
- Transmisividad, T: es el parámetro definido para evaluar la capacidad de transmitir agua

de los acuíferos, teniendo en cuenta no sólo la textura del acuífero y las características del
fluido, sino también las características estructurales o geométricas. Se define como el
producto de la conductividad hidráulica k y el espesor del acuífero b:
T=k·b
Acuíferos muy permeables pero de bajo espesor pueden, a pesar de sus excelentes
características texturales, no ser aptos para una eficaz transmisión del agua debido a su
bajo espesor y, consecuentemente, baja transmisividad.
3.2 Piezometría
7
El movimiento del agua en las formaciones geológicas está condicionado por la capacidad
de realizar trabajo de sus partículas elementales, es decir, por su energía. La energía
depende, de una parte, del campo de fuerzas a que está sometido el fluido y del agente
sensible a dicho campo y, de otra, del tipo de energía en sí o forma en que el cuerpo tenga
almacenada su capacidad de desarrollar trabajo.
Los tres tipos de energía que pueden incidir sobre el agua que se mueve en el seno de una
formación permeable, asumiendo que no van a existir cambios de calor, son las siguientes:
Energía potencial: Eh = m*g*z
Energía cinética: Ec = ½*m*v²
Energía de presión: Ep = p*V
Donde g representa la intensidad del campo gravitatorio y m la masa o agente sensible al

campo gravitatorio, z es la altura respecto a una cota de referencia, v es la velocidad del
fluido (agua), p es la presión a que está sometido el fluido y V es el volumen. Por el
principio de conservación de la energía, ésta va a permanecer constante a lo largo de una
línea de corriente, y por tanto:
E = Eh + Ec + Ep = cte (Caso de los fluidos ideales y cuando el régimen es estacionario)
O bien: m*g*z + ½*m*v² + p*V = cte
𝑣² 𝑝
Sacando m*g factor común: m*g * (z + 2∗𝑔 + 𝜌∗𝑔
) = cte
Siendo ρ la densidad del agua.

Ante un mismo fluido y un mismo campo gravitatorio, el terrestre, se tiene:
𝑣² 𝑝
(z + 2∗𝑔 + 𝜌∗𝑔
) = cte
Que no es más que la expresión del teorema de Bernoulli y en la que todos sus términos
tienen dimensiones de longitud:
8
[z] = L
𝑣² 𝐿²𝑇¯²
[2∗𝑔] = 𝐿𝑇²
=L
𝑝 𝑀𝐿𝑇¯²𝐿¯²
[𝜌∗𝑔] = 𝑀𝐿¯³𝐿𝑇¯² = L
En consecuencia, a los términos de la ecuación se les denomina alturas:
z → Altura de posición
𝑣²
2∗𝑔
→ Altura de velocidad
𝑝
𝜌∗𝑔
→ Altura de presión
El término v²/2*g, o altura de velocidad, en el caso del movimiento del agua en el seno de
un medio permeable, puede ser despreciado frente a la altura de presión y la altura de
posición, dado lo bajo del valor de la velocidad del agua en dichos medios y, por tanto, la
capacidad para realizar trabajo por parte del agua en un punto A de un acuífero, vendrá
dada por su nivel piezométrico o suma de las alturas de posición y presión:
𝑃𝑎
ha = za + 𝛾
Siendo 𝛾 el peso específico del agua.
Al nivel piezométrico se le suele denominar también altura piezométrica, cota piezométrica

o carga hidráulica. El nivel piezométrico se mantiene constante en la vertical de un
almacenamiento de agua libre o en la vertical de un medio permeable saturado cuando éste
es homogéneo e isótropo.
9
Figura 3.2. Nivel piezométrico. Fuente: Ingeniería Geológica. González de Vallejo, 2002
Tal y como se ve en la Figura 3.2, el punto A está a presión atmosférica y por tanto su
nivel piezométrico, ha, es igual a la altura de posición za. El punto B tiene un nivel
piezométrico, hb, que se quiere comparar con ha.
Sea Pb la presión en el punto B y zb su altura de posición:
ha = za
𝑃𝑏 (𝑍𝑎−𝑍𝑏) ∗ 𝛾
hb = zb + 𝛾
= zb + 𝛾
= za = ha
Es decir, el nivel piezométrico es el mismo en los puntos A y B.
Figura 3.3. Teorema de Bernoulli

Fuente: http://fisica24im10semestre2014.blogspot.com.es/2014/04/actividad-12-principio-de-
bernoulli-y.html
10
3.3 La ecuación de Darcy
La ley de Darcy fue enunciada por su autor en 1856, después de diversos trabajos y
experimentaciones. Establece que el caudal Q que es capaz de atravesar un medio
permeable, es proporcional al área de paso del medio permeable A, normal al flujo, y al
gradiente de niveles piezométricos entre la entrada y la salida del flujo en el medio
permeable, i. La constante de proporcionalidad es la permeabilidad del medio, en la que
quedan incluidas las características del fluido, es decir, es la permeabilidad efectiva,
conductividad de Darcy o conductividad hidráulica. En consecuencia y con carácter global:
Q = -k * A * i
Siendo: Q = Caudal de paso [L³T¯¹]

k = Permeabilidad de Darcy [LT¯¹]
A = Área de paso [L²]
i = Gradiente piezométrico
El gradiente de un escalar es un vector, definido por su módulo, dirección y sentido. El

módulo es la derivada direccional de la función escalar, la dirección es la normal a las
curvas de nivel de la función escalar y el sentido es el del aumento de valor de la función.
Al moverse el flujo de zonas de mayor a menor nivel piezométrico, que es en este caso la
función escalar, el flujo y el gradiente tienen direcciones distintas, lo que justifica el signo
negativo de la ecuación de Darcy. Llamando h al nivel piezométrico y x al espacio en el
sentido del flujo, la ecuación puede expresarse en forma diferencial, considerando la
variación de niveles piezométricos a lo largo del recorrido del flujo por el medio poroso:
𝑑ℎ
Q = -k * A * 𝑑𝑥
En el caso de la Figura 3.4, el experimento mantiene constantes los niveles h1 y h2, y al ser
constantes k y A para todo el medio permeable se establece la forma más sencilla de Darcy:
ℎ1 − ℎ2
𝑄 =𝑘∗𝐴∗
𝐿
11
Figura 3.4. Ley de Darcy

Fuente: http://geolabs.com.mx/blog/trabajo/propiedades-hidraulicas-del-suelo/
La ley de Darcy es únicamente válida para el caso en que el régimen sea laminar. Para
evaluar si el régimen es realmente laminar hay que acudir al denominado número de
Reynolds, que representa la relación entre las fuerzas de inercia y viscosas en el seno de un
fluido:
𝜌∗𝑣∗𝑑
Re = 𝜇
donde: v = Velocidad del fluido

d = Tamaño medio de las partículas, se considera el d10 en medios
granulares y 2e en medios fisurados, siendo e el ancho medio de las fisuras.
ρ = Densidad del fluido
𝜇 = Viscosidad dinámica
El número de Reynolds es adimensional.
(𝑀𝐿¯³)(𝐿)
[Re] = (𝑀𝐿𝑇¯²)(𝑇)(𝐿¯²) → Adimensional
De cualquier modo, para su cálculo habrá que utilizar unidades homogéneas; por ejemplo,
v en cm/s, d en cm, ρ en g/cm³ y 𝜇 en dina * seg/cm². Un valor del número de Reynolds
comprendido entre 1 y 10 denota un régimen laminar en el que puede considerarse
12
aplicable la ley de Darcy, aunque es particularmente fiable para valores inferiores a 4. El

régimen se manifiesta totalmente turbulento para valores de 60 a 180, y en los valores
intermedios puede asumirse como cuasi-turbulento. En cualquier caso, con esta última
gama de parámetros, el flujo no cumpliría la ley de Darcy.
3.4 La ecuación general de flujo: Ecuación de Laplace y ecuación de Poisson
Ecuación de Laplace: La ecuación de Laplace es la primera de las ecuaciones diferenciales

en derivadas parciales que gobiernan el flujo en medios porosos (permeables en general,
cuando el medio pueda ser asumido como poroso). Se refiere al flujo en condiciones de
régimen permanente y en ausencia de fuentes y sumideros.
Esta ecuación, como todas las ecuaciones del flujo en medios porosos, se obtiene aplicando
conjuntamente la ecuación de la continuidad y la ley de Darcy. El flujo se estudia
aceptando, de una parte, que la masa entrante debe ser igual a la saliente, y que ambas,
dentro del medio poroso, están condicionadas por los planteamientos de Darcy. Así, se
tiene:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧

+ + = 0 (continuidad)
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕ℎ 𝜕ℎ 𝜕ℎ
𝑣𝑥 = − 𝑘𝑥𝑥 ∗ 𝜕𝑥 ; 𝑣𝑦 = − 𝑘𝑦𝑦 ∗ 𝜕𝑦 ; 𝑣𝑧 = − 𝑘𝑧𝑧 ∗ 𝜕𝑧
sustituyendo:
𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ
�−𝑘𝑥𝑥 ∗ � + �−𝑘𝑦𝑦 ∗ � + �−𝑘𝑧𝑧 ∗ � = 0
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧
Si se considera el medio homogéneo (la permeabilidad es la misma en cualquier punto del

espacio) e isótropo (la permeabilidad es la misma en cualquier dirección del espacio), se
tiene que:
𝑘𝑥𝑥 = 𝑘𝑦𝑦 = 𝑘𝑧𝑧 = 𝑘
y:
𝜕 2ℎ 𝜕 2ℎ 𝜕 2ℎ
+ + =0
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2
13
Que es la ecuación de Laplace (régimen permanente, medio homogéneo e isótropo y

ausencia de fuentes y sumideros). Hay dos tipos característicos de condiciones de
contorno:
- Condiciones Dirichlet: nivel piezométrico conocido en frontera.
- Condiciones Neuman: flujo conocido en frontera.
La resolución de la ecuación solo es posible por métodos analíticos en sistemas físicos

sencillos, como es el caso de flujo hacia pozos y algunos problemas de flujo entre zanjas.
Para casos reales más complejos debe acudir se a métodos aproximados de resolución
numérica, o también, en algunos otros casos, a métodos gráficos de resolución, como redes
de flujo.
Ecuación de Poisson: Laplace no integra en sus planteamientos la existencia de fuentes y

sumideros. Sin embargo, y aún en problemas de régimen permanente, es preciso incluir
bombeos y recargas (fuentes y sumideros) asociados a problemas reales. La ecuación de
Poisson de flujo en medios porosos en régimen permanente, y con posibilidad de integrar
bombeos y recargas, puede desarrollarse con la ayuda del esquema de la Figura 3.5.
Figura 3.5. Esquema para la deducción de la ecuación de Poisson

Fuente: Ingeniería Geológica. González de Vallejo, 2002
Se considera una recarga R(x, y) por unidad de área y unidad de tiempo. Aplicando
condiciones de flujo estacionario:
14
𝜕𝑣𝑥 𝜕𝑣𝑦
∆𝑥(𝑏∆𝑦) = ∆𝑦(𝑏∆𝑥) = 𝑅(𝑥, 𝑦)∆𝑥∆𝑦
𝜕𝑥 𝜕𝑦
Aplicando Darcy con hipótesis de homogeneidad e isotropía
𝑘𝑥𝑥 = 𝑘𝑦𝑦 = 𝑘𝑧𝑧 = 𝑘
y siendo
𝑇 = 𝑘∗𝑏
se obtiene la ecuación de Poisson:
𝜕 2ℎ 𝜕 2ℎ 𝑅(𝑥, 𝑦)
+ = −
𝜕𝑥 2 𝜕𝑦 2 𝑇
En régimen permanente, presencia de fuentes y sumideros y permeabilidad homogénea e

isótropa. En esta ecuación se tiene:
- Si R(x, y) = 0, es coincidente con la ecuación de Laplace
- Si R(x, y) es positiva = Recarga
- Si R(x, y) es negativa = Bombeo
Todo lo dicho para la resolución de la ecuación de Laplace es válido para la de Poisson.

Sin embargo, esta última es más realista y completa y se ajusta mejor a la mayoría de los
problemas reales que suelen presentarse. Establece de modo simplista que la suma de
entradas y salidas de agua en un recinto acotado debe ser igual al bombeo o recaiga que
incida sobre tal recinto, siempre en condiciones de régimen permanente (no variación de
los niveles piezométricos en el recinto) y de medio homogéneo e isótropo en cuanto a la
permeabilidad.
Ecuación del flujo en régimen transitorio: En el régimen permanente se obtienen
soluciones de equilibrio para un fenómeno concreto. No se trata de cómo ni cuándo se
15
alcanza ese equilibrio, únicamente se calcula la solución de equilibrio. Si un sistema está

en equilibrio, para cambiar esta situación se requieren ciertas acciones sobre el sistema,
tales como bombeos o recaigas.
El sistema evoluciona hasta alcanzar de nuevo la situación de equilibrio correspondiente;

mientras esta nueva situación es alcanzada, se suceden un conjunto de situaciones a lo
largo del tiempo, que es lo que se denomina régimen transitorio. Los niveles van variando
y en consecuencia se va almacenando o desalmacenando agua en el recinto. En este caso,
la ecuación de la continuidad, en régimen transitorio, se aplica al recinto según:
Volumen saliente = Volumen entrante – Volumen desalmacenado

(En un intervalo de tiempo ∆t)
Figura 3.6. Esquema para la deducción de la ecuación general del flujo transitorio
En el intervalo ∆t, el nivel varía ∆h, como indica la Figura 3.6. El volumen de roca vaciado
vendrá dado por ∆h * ∆x * ∆y, y el volumen de agua contenido en este volumen de roca
vendrá dado por S * ∆h * ∆x * ∆y. El volumen desalmacenado en el tiempo ∆t será S *
∆h * ∆x * ∆y/ ∆t. Introduciendo este término en la ecuación de la continuidad en presencia
de fuentes y sumideros, es decir, incluyendo el término del volumen desalmacenado en la
unidad de tiempo en el término de recargas, se tiene:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕ℎ
∆𝑥(𝑏∆𝑦) + ∆𝑦(𝑏∆𝑥) = 𝑅(𝑥, 𝑦, 𝑡)∆𝑥∆𝑦 − 𝑆 (∆𝑥∆𝑦)
𝜕𝑥 𝜕𝑦 𝜕𝑡
Aplicando Darcy generalizado como en anteriores ecuaciones:
16
𝜕ℎ 𝜕ℎ 𝜕ℎ
𝑣𝑥 = − 𝑘𝑥𝑥 ∗ 𝜕𝑥 ; 𝑣𝑦 = − 𝑘𝑦𝑦 ∗ 𝜕𝑦 ; 𝑣𝑧 = − 𝑘𝑧𝑧 ∗ 𝜕𝑧
Asumiendo una vez más condiciones de homogeneidad e isotropía (kxx = kyy = kzz = k),
discretización cuadrada (∆x = ∆y) y T = k*b, se obtiene:
𝜕 2 ℎ 𝜕 2ℎ 𝑆 𝜕ℎ 𝑅(𝑥, 𝑦, 𝑡)
2
+ 2
= −
𝜕𝑥 𝜕𝑦 𝑇 𝜕𝑡 𝑇
Ecuación general del flujo en régimen transitorio, presencia de fuentes y sumideros y
permeabilidad homogénea e isótropa. Los términos de esta ecuación tienen el siguiente
significado:
𝜕2 ℎ 𝜕2 ℎ
𝜕𝑥 2
+ 𝜕𝑦 2
∶ suma de caudales entrantes por caras laterales
𝑆 𝜕ℎ
𝑇 𝜕𝑡
∶ volúmenes desalmacenados en la unidad de tiempo
𝑅(𝑥,𝑦,𝑡)
𝑇
∶ bombeos o recargas en la unidad de tiempo
En definitiva, esta ecuación indica que el balance de caudales entrantes y salientes en un

recinto por razones de gradientes piezométricos, más las recargas externas al sistema,
como bombeos y recargas, debe ser igual a los caudales de llenado o vaciado del recinto.
𝑆 𝜕ℎ
Como se puede observar, si 𝑇 𝜕𝑡
= 0 , se obtiene la ecuación de Poisson de régimen
𝑆 𝜕ℎ 𝑅(𝑥,𝑦,𝑡)
permanente con fuentes y sumideros. Si 𝑇 𝜕𝑡
=0y 𝑇
= 0 , se obtiene la ecuación de
Laplace.
La ecuación general del flujo y su resolución para cada caso concreto con sus adecuadas
condiciones de contorno es de uso general en los problemas de hidrogeología e ingeniería
geológica.
3.5 Resolución analítica para el diseño de abatimientos mediante ensayos de

bombeo
Los ensayos de bombeo son, sin lugar a dudas, el método más completo y fiable para el
cálculo de los parámetros hidrogeológicos de un acuífero. Consisten en bombear un pozo,
17
en principio a caudal constante, y analizar el descenso de los niveles piezométricos, tanto

en el pozo de bombeo como en pozos o piezómetros circundantes.
Existen dos tipos básicos de métodos: ensayos de bombeo en régimen permanente y
ensayos de bombeo en régimen transitorio o variable. En los primeros los niveles
piezométricos no varían, y se interpretan los descensos producidos en el entorno del pozo
como consecuencia de un bombeo constante en el mismo; en régimen variable lo que se
analiza es la evolución de los niveles en el pozo de bombeo y pozos de observación a lo
largo de la prueba.
3.5.1 Ensayos de bombeo en régimen permanente. Deducción de las ecuaciones
La ecuación general del flujo bidimensional bajo condiciones de acuífero confinado y

régimen permanente viene dada por:
𝜕 2ℎ 𝜕 2ℎ
+ =0
𝜕𝑥 2 𝜕𝑦 2
siendo h el nivel piezométrico. Ello implica que el acuífero pueda ser considerado
homogéneo, isótropo e infinito, junto a la no existencia de recargas verticales.
Para el caso de flujo radial hacia un punto donde se sitúa el pozo de bombeo, puede
transformarse la ecuación a coordenadas polares, tomando como eje de coordenadas el eje
del pozo, que será un punto al proyectarlo sobre el plano xy. La transformación se realiza
según:
𝑥 = 𝑟 cos 𝛼
𝑦 = 𝑟 𝑠𝑒𝑛 𝛼
Anulando los términos en α, se tiene la ecuación:
𝜕 2 ℎ 1 𝜕ℎ
+ =0
𝜕𝑟 2 𝑟 𝜕𝑟
18
que al depender solo de r puede escribirse de la forma:
1 𝑑 𝑑ℎ
�𝑟 � = 0
𝑟 𝑑𝑟 𝑑𝑟
por lo tanto:
𝑑ℎ
�𝑟 � = 𝑐𝑡𝑒
𝑑𝑟
Para calcular la constante, solo hay que suponer que el caudal que sale por el pozo debe ser
igual al que atraviesa una superficie cilíndrica de altura constante b y radio genérico r. En
este caso, según Darcy:
𝑑ℎ
𝑄 = 2𝜋𝑟𝑏𝑘
𝑑𝑟
luego:
𝑑ℎ 𝑄 𝑄
𝑐𝑡𝑒 = 𝑟 = =
𝑑𝑟 2𝜋𝑏𝑘 2𝜋𝑇
y por tanto:
𝑄 𝑑𝑟
𝑑ℎ =
2𝜋𝑇 𝑟
Si se realiza la integración entre el radio genérico r y el radio de influencia R, distancia a la

cual el descenso es nulo, se tiene que el nivel variará entre el valor h a distancia genérica y
el valor ℎ0 representativo de niveles en reposo a escala regional:
ℎ0 𝑅
𝑄 𝑑𝑟
� 𝑑ℎ = �
ℎ 2𝜋𝑇 𝑟 𝑟
obteniéndose la solución de Thiem (1906):

𝑄 𝑅
ℎ0 − ℎ = 𝑙𝑛
2𝜋𝑇 𝑟
donde:
ℎ0 = nivel piezométrico inicial
h = nivel piezométrico a la distancia r
19
Q = caudal de bombeo
T = transmisividad
R = radio de influencia
r = distancia al eje del pozo
Todas estas fórmulas parten de la aceptación de que la transmisividad es constante en

cualquier punto del espacio; es decir, asumen un medio homogéneo e isótropo. A la
relación de Thiem también puede llegarse por métodos más simples, sin acudir a la
ecuación general. A continuación se incluye este último desarrollo matemático por ser
bastante representativo e intuitivo.
Según la Figura 3.7, si se considera un acuífero cautivo o confinado en el que un pozo

perforado en el mismo bombea un caudal constante Q, y existen dos pozos de observación
(nº 1 y nº 2) estando el nivel estabilizado en todo el cono de bombeo, se puede establecer
que el caudal que sale por el pozo es igual al que atraviesa una superficie cilíndrica ideal
de radio genérico r y altura el espesor del acuífero b.
Figura 3.7. Esquema para la deducción de la fórmula de Thiem

Las depresiones en el pozo de bombeo, y en los piezómetros nº 1 y nº 2 serán 𝑑𝑝 , 𝑑1 y 𝑑2 ,
y los niveles ℎ𝑝 , ℎ1 y ℎ2 , respectivamente. Los piezómetros se encuentran a las distancias
𝑟1 y 𝑟2 del pozo de bombeo. A la distancia genérica r, la depresión y el nivel h.
20
Según la ley de Darcy, puede establecerse que el caudal que atraviesa la superficie
cilíndrica porosa será igual al producto de la permeabilidad del acuífero por el área de paso
y por el gradiente hidráulico que existe entre uno y otro lado de la superficie cilíndrica
considerada. Es decir:
𝑄 = 𝑘𝐴𝑖
donde A = 2πrb es el área de paso, i = dh/dr es el gradiente
𝑑ℎ
𝑄 = 𝑘2𝜋𝑟𝑏 𝑘𝑏 = 𝑇
𝑑𝑟
𝑑ℎ 𝑄 𝑑𝑟
𝑄 = 2𝜋𝑟𝑇 𝑑ℎ = ·
𝑑𝑟 2𝜋𝑇 𝑟
Para integrar esta última expresión, se tienen que definir los límites de integración. Cuando
h varía entre un nivel genérico h y el nivel inicial del acuífero ho' variará entre un radio
genérico r, en el que el nivel vale h, y una distancia R (radio de influencia) a la cual el
descenso es nulo y, por tanto, h vale ℎ0 :
ℎ0 𝑅
𝑄 𝑑𝑟 𝑄 𝑅
� 𝑑ℎ = � y ℎ0 − ℎ = 𝑙𝑛
ℎ 2𝜋𝑇 𝑟 𝑟 2𝜋𝑇 𝑟
o bien:
𝑄 𝑅
𝑑= 𝑙𝑛
2𝜋𝑇 𝑟
3.5.2 Formulación analítica en el diseño de abatimientos de acuíferos confinados. Flujo

proveniente de una fuente circular; caso de pozo simple y pozos múltiples
FLUJO RADIAL EN ACUIFEROS CONFINADOS

Si consideramos una capa de espesor constante D, de arena libremente drenante cargada
con agua y confinada arriba y abajo esencialmente por estratos impermeables, como se
21
muestra en la Figura 3.9. Si se instala un pozo cilíndrico de manera que penetre por
completo en la arena, el agua comenzará a llenar el pozo hasta una altura determinada por
el nivel piezométrico de equilibrio en la arena. Al bombear el agua del pozo su nivel
desciende, la cabeza total en la arena situada inmediatamente al lado del pozo se reduce y
por tanto el agua fluye a través de la arena de la zona de cabeza total alta hasta la zona
alrededor del pozo, y en consecuencia hacia el interior del mismo. Cuando el caudal en el
pozo es constante e igual a su caudal de bombeo, se establece la condición de flujo
estacionario y el nivel piezométrico en la arena desciende como lo muestra en la Figura
3.9. El flujo en cualquier punto de la arena se produce solo en el plano horizontal, siempre
y cuando el nivel piezométrico no descienda por debajo del nivel superior de la arena, las
fronteras del flujo superior e inferior están definidas por las superficies de la arena por
arriba y por abajo. Estas son las condiciones para tener un flojo confinado y el estrato de
arena se denomina acuífero confinado.
Figura 3.9. Flujo en un acuífero confinado proveniente de una fuente circular: Pozo simple
Fuente: Mecánica de Suelos. Peter L. Berry & David Reid, 1993
FLUJO PROVENIENTE DE UNA FUENTE CIRCULAR

Aunque muchos depósitos sedimentarios tienen una permeabilidad mucho más alta en la
dirección horizontal que en la vertical, las variaciones de permeabilidad en el plano
horizontal son a menudo mucho menos pronunciadas. Si se supone que la arena es
isotrópica en el plano horizontal, entonces el flujo tendrá una simetría radial alrededor del
pozo y el abatimiento del nivel piezométrico será el mismo en todos los puntos situados a
22
igual distancia del pozo. Además cuando la distancia al pozo aumenta, el abatimiento del
nivel piezométrico disminuye, hasta una distancia radial re a la cual el pozo no tiene
influencia significativa en el nivel piezométrico de la arena. Esta distancia re es el radio de
influencia del pozo. Por tanto puede considerarse que el flujo proviene de una fuente
circular localizada en el radio de influencia del pozo.
Como el flujo en cualquier punto de la arena se produce en el plano horizontal x, y para el

cual se supone kx = ky, la ecuación de flujo estacionario se reduce a una ecuación de
Laplace en dos dimensiones:
𝜕²ℎ 𝜕²ℎ
+ =0
𝜕𝑥² 𝜕𝑦²
Para el flujo radial es más conveniente trabajar en coordenadas polares r y ϴ que en la

Figura 3.9 se relacionan con x e y mediante:
𝑦
𝑟 = �𝑥 2 − 𝑦 2 𝛳 = 𝑎𝑟𝑐𝑡𝑎𝑛 �𝑥 �
Al transformar la ecuación de Laplace a coordenadas polares se obtiene:
𝜕²ℎ 1 𝜕ℎ 1 𝜕²ℎ
+ + =0
𝜕𝑟² 𝑟 𝜕𝑟 𝑟² 𝜕𝛳²
Para un flujo simétrico alrededor del eje z, h es la misma en todos los puntos situados a
igual distancia del pozo. Entonces, para cualquier valor de r, 𝜕ℎ/𝜕𝛳= 0, y la ecuación se
reduce a una diferencial ordinaria:
𝑑²ℎ 1 𝑑ℎ 1𝑑 𝑑ℎ
+ =0→ �𝑟 � = 0
𝑑𝑟² 𝑟 𝑑𝑟 𝑟 𝑑𝑟 𝑑𝑟
23
Se observa ahora que h es función solo de la variable r y el problema se reduce

matemáticamente a un caso unidimensional. Al resolver esta última ecuación pueden
obtenerse soluciones para los valores de la cabeza total en el dominio del flujo y para el
caudal hacia el pozo.
Pozo simple: Si en la Figura 3.9 Se toma la base de la arena como nivel de referencia, la
cabeza total h en cualquier radio r está dada por la altura del nivel piezométrico desde la
base. Si se designa el radio del pozo como 𝑟𝑤 , la cabeza total en la arena adyacente al pozo
como ℎ𝑤 y la cabeza total en el radio de influencia como ℎ𝑒 , las condiciones de frontera
son:
ℎ = ℎ𝑤 𝑒𝑛 𝑟 = 𝑟𝑤 y ℎ = ℎ𝑒 𝑒𝑛 𝑟 = 𝑟𝑒
Si se integra una vez la ecuación descrita justo antes:
𝑑ℎ
𝑟 = 𝑐1
𝑑𝑟
Y si se integra por segunda vez se obtiene:
ℎ = 𝑐1 ln 𝑟 + 𝑐2
donde 𝑐1 y 𝑐2 son constantes de integración. Al sustituir las condiciones de frontera en la

ecuación anterior y resolver el sistema de ecuaciones simultáneas se obtiene:
ℎ𝑒 −ℎ𝑤 ℎ𝑒 −ℎ𝑤
𝑐1 = 𝑐2 = ln 𝑟𝑒
ln 𝑟𝑒 /𝑟𝑤 ln 𝑟𝑒 /𝑟𝑤
Si se sustituyen las constantes en la segunda integral se obtiene:
ℎ𝑒 − ℎ𝑤 𝑟𝑒
ℎ = ℎ𝑒 − ln
ln 𝑟𝑒 /𝑟𝑤 𝑟
24
El principio de continuidad requiere que el caudal estacionario Q en el pozo sea igual al

flujo a través de una sección cilíndrica de radio r en la arena. Por tanto:
𝑑ℎ 𝑑ℎ
𝑄 = 𝑘𝑖𝐴 = 𝑘 2𝜋D = 2𝜋𝑘Dr
𝑑𝑟 𝑑𝑟
𝑑ℎ
La ecuación establece que 𝑟 𝑑𝑟 = 𝑐1 y la ecuación define una expresión para 𝑐1, al
remplazar se obtiene:
2𝜋𝑘D(ℎ𝑒 − ℎ𝑤 )
𝑄=
ln 𝑟𝑒 /𝑟𝑤
Al reordenar la ecuación y sustituirla en la anterior de h, se obtiene la siguiente expresión

para la cabeza total h en función de Q para cualquier radio r:
𝑄 𝑟𝑒
ℎ = ℎ𝑒 − ln
2𝜋𝑘D 𝑟
Estas ecuaciones dependen del radio de influencia 𝑟𝑒 . Sin embargo, una determinación
precisa de 𝑟𝑒 es difícil debido a las dificultades inherentes a la estimación de la distancia a
la cual el pozo deja de tener influencia. Además, el radio de influencia del pozo no es
contante, un incremento en el abatimiento produce un aumento de 𝑟𝑒 y en cuanto mayor sea
la permeabilidad de la arena mayor será el radio de influencia del pozo. Afortunadamente
la influencia de 𝑟𝑒 está en términos de logaritmo natural y, por tanto, la solución no es
sensible a pequeñas inexactitudes. Si diésemos valores típicos a 𝑟𝑒 dentro de un rango de
valores de Q que se esperan para una arena, observaríamos que al duplicar el valor de 𝑟𝑒
los cambios de Q son inferiores al 10%. Por lo que las inexactitudes no serán
significativas. En la siguiente tabla podemos ver un ejemplo:
25
Figura 3.10. Efecto de valores de 𝑟𝑒 dentro del rango de valores de Q que se esperan en arena
Debido a la incertidumbre asociada con el valor de 𝑟𝑒 y al hecho de que el nivel del agua en
el pozo será inferior a ℎ𝑤 en un valor equivalente a la pérdida de cabeza hidráulica en el
pozo, algunas veces los niveles piezométricos en la arena se registran directamente al
instalar pozos de observación a distancias adecuadas del pozo de bombeo. En la Figura 3.9,
por ejemplo, los pozos de observación a las distancias 𝑟1 y 𝑟2 registran niveles
piezométricos correspondientes a ℎ1 y ℎ2 , respectivamente. Al sustituir estos valores en la
ecuación ℎ = 𝑐1 ln 𝑟 + 𝑐2 , y resolverla como se hizo antes, se obtienen las siguientes
expresiones alternas para el caudal Q en el pozo y para la cabeza total h en cualquier radio
r:
2𝜋𝑘D(ℎ2 − ℎ1 ) 𝑄 𝑟1
𝑄= 𝑦 ℎ = ℎ1 − ln
ln 𝑟𝑒 /𝑟𝑤 2𝜋𝑘D 𝑟
Pozos múltiples: Las excavaciones en arcillas situadas sobre acuíferos de arena o grava
tienden a fallar por levantamiento del fondo. Para asegurar su estabilidad a menudo se
necesita reducir el nivel piezométrico en el acuífero. Esto puede obtenerse mediante un
sistema de pozos de alivio que penetren completamente en el acuífero, los cuales se
instalan alrededor del perímetro de la excavación propuesta, como se observa en la Figura
3.11. Los cálculos de estabilidad definen el mínimo abatimiento de nivel piezométrico
necesario en el área de la excavación y para verificar que el sistema de pozos propuesto es
capaz de realizar dicho abatimiento se necesita una expresión para la cabeza total en
cualquier punto P de la arena que resulta de bombear en varios pozos simultáneamente.
Ahora para el flujo hacia un pozo simple situado en la posición 1, la cabeza total en P se
𝑄 𝑟
obtiene de la ecuación ℎ = ℎ𝑒 − 2𝜋𝑘D ln 𝑟𝑒, de la manera siguiente:
26
𝑄1 𝑟𝑒1
ℎ𝑝1 = ℎ𝑒 − ln
2𝜋𝑘D 𝑟1
Y el abatimiento en P está dado por:
𝑄1 𝑟𝑒1
ℎ𝑒 − ℎ𝑝1 = ln
2𝜋𝑘D 𝑟1
Donde ℎ𝑝1 =cabeza total en P debida al bombeo únicamente en el pozo 1.

𝑄1= caudal estacionario en el pozo 1
𝑟1= distancia de P a la línea central del pozo 1
𝑟𝑒1=radio de influencia del pozo 1
𝑄 𝑟
Como la ecuación ℎ = ℎ𝑒 − 2𝜋𝑘D ln 𝑟𝑒, es una solución de una ecuación diferencial lineal
(una ecuación de la Place), para calcular el abatimiento total en P debido a varios pozos
que actúan juntos puede aplicarse el principio de superposición mediante la suma
algebraica de los efectos separados de cada pozo individual. Por tanto, si hay n pozos el
abatimiento total en P está dado por:
1 𝑟𝑒1 𝑟𝑒2 𝑟𝑒3 𝑟𝑒𝑛

ℎ𝑒 − ℎ𝑝 = �𝑄1 ln + 𝑄2 ln + 𝑄3 ln + ⋯ 𝑄𝑛 ln �
2𝜋𝑘D 𝑟1 𝑟2 𝑟3 𝑟𝑛
Donde los subíndices 2, 3, … , n corresponden a los pozos 2, 3 , … , n, respectivamente.

Entonces:
𝑛
1 𝑟𝑒𝑗
ℎ𝑒 − ℎ𝑝 = � Q𝑗 ln
2𝜋𝑘D 𝑟𝑗
𝑗=1
En la práctica, es usual tomar medidas para que las características del abatimiento asociado
con cada pozo sean similares y entonces puede tomarse el mismo radio de influencia 𝑟𝑒
para cada pozo. En todo caso, las pequeñas inexactitudes en el valor de 𝑟𝑒 no son
significativas. Para alcanzar un abatimiento similar en cada pozo, éstos pueden localizarse
de tal manera que cada uno intercepte aproximadamente la misma cantidad de filtración y
que pueda bombearse más o menos la misma cantidad de agua desde cada pozo. Como el
flujo proveniente de una fuente circular, para el correcto posicionamiento de los pozos solo
será necesario que estén espaciados a igual distancia alrededor de la excavación. Cuando se
27
bombea desde pozos profundos es común utilizar bombas sumergibles localizadas bajo el
nivel del agua en los pozos. Si se utiliza el mismo tipo de bomba en todo el sistema y cada
bomba funciona a la misma velocidad bajo las mismas condiciones de carga, entonces el
caudal de cada bomba será aproximadamente el mismo. Para estas condiciones puede
suponerse que el caudal Q en cada pozo es el mismo e igual a la capacidad de operación de
las bombas. Por consiguiente, al tomar Q y 𝑟𝑒 constantes para cada pozo, la ecuación queda
reducida a:
𝑛
𝑄 𝑟𝑒
ℎ𝑒 − ℎ𝑝 = � ln
2𝜋𝑘D 𝑟𝑗
𝑗=1
Donde para la cabeza total en el punto P se obtiene la expresión siguiente:
𝑛
𝑄 𝑟𝑒
ℎ𝑝 = ℎ𝑒 − � ln
2𝜋𝑘D 𝑟𝑗
𝑗=1
Q es la capacidad de operación de las bombas y 𝑟𝑗 es la distancia de P a la línea central del

pozo j, donde j = 1, 2, 3, …, n.
Figura 3.11. Flujo en un acuífero confinado proveniente de una fuente circular: Pozo múltiple.
28
3.6 Fuerzas de filtración. El fenómeno del sifonamiento
Resulta intuitivo comprender que el agua en movimiento a través del terreno induce sobre
éste una fuerza que tiende a arrastrarlo. De hecho, como se ha expuesto anteriormente, para
que exista flujo es necesaria una diferencia de altura piezométrica, de manera que el agua
circule desde puntos de mayor carga ℎ𝑎 hacia puntos de menor carga ℎ𝑏 . La diferencia
𝛥ℎ = ℎ𝑎 − ℎ𝑏 representa el trabajo o energía gastados en resistir la fuerza de arrastre del
agua al moverse a través de los pequeños poros del suelo.
Así, si las fuerzas que resisten al flujo son menores que la fuerza erosiva de éste, las
partículas de suelo podrán ser arrastradas por el agua. Este fenómeno puede originar serios
problemas en diversas aplicaciones geotécnicas como se muestra en la Figura 3.12.
Figura 3.12. a) Tubificación de una presa de materiales por erosión interna.

b) Tubifiación de una excavación por erosión interna
Las fuerzas resistentes a la erosión dependen de la cohesión del suelo, de su granulometría,

compacidad y densidad, siendo las arenas finas uniformes y flojas los suelos más
susceptibles de sufrir el arrastre del agua. Por su parte la fuerza de arrastre del agua
depende del gradiente hidráulico (i=Δh/l).
Si se produce una concentración de flujo (y gradiente) suficiente en las proximidades de la

superficie de salida del flujo, las primeras partículas de suelo pueden ser arrastradas. A
29
continuación, al mantenerse aproximadamente la diferencia de altura piezométrica (Δh)

pero disminuir el recorrido de la filtración (L) por la pérdida de suelo, se incrementa el
gradiente y también la fuerza erosiva. De esta forma, el mecanismo de erosión puede
progresar hacia el interior del terreno hasta, en caso extremo, conducir a la ruina de la obra.
Desde un punto de vista teórico, la forma más sencilla y clásica de establecer las
condiciones de equilibrio del suelo frente a las fuerzas de filtración es la mostrada en la
Figura 3.13. En ella se representan tres situaciones posibles de flujo vertical en un sencillo
dispositivo llamado permeámetro de carga constante.
Se trata de un recipiente en donde se encuentra una masa de suelo de altura (L), que se
supondrá constituido por una arena fina, confinado entre dos rejillas. Por encima del suelo
existe una lámina de agua libre, cuya cota se mantiene constante en todo momento (punto
D). Por debajo del suelo el permeámetro se conecta a un conducto que termina en otro
recipiente anexo, siempre lleno de agua (hasta el punto A), pero que puede moverse hacia
arriba o hacia abajo. Finalmente, desde el interior de la masa de suelo se pueden disponer
unos piezómetros abiertos (P1, P2 y P3). Tomando como plano de referencia (z=0) la cota
inferior de la masa de suelo, resulta sencillo controlar en todo momento las alturas
piezométricas. Basta para ello con colocar una regleta graduada como la mostrada en la
Figura 3.13 y recordar que la cota de agua leída en la regleta será igual a la altura
piezométrica (ℎ = 𝑧 + 𝑢/ϒ𝑤 ).
30
Figura 3.13. Condiciones de equilibrio en presencia de filtración

En el caso de la Figura 3.13a) la cota de agua del recipiente anexo (A) se hace coincidir
con la cota de agua en el permeámetro (D), de manera que resulta inmediato comprobar
que las condiciones son hidrostáticas, sin flujo de agua:
ℎ𝐴 = ℎ𝐵 = ℎ𝐶 = ℎ𝐷 = ℎ𝑃1,𝑃2,𝑃3
Lo que puede comprobarse fácilmente si más que observar que los niveles de agua en el
permeámetro, el recipiente anexo y en los piezómetros alcanzan la misma cota. Las
presiones intersticiales en los puntos extremos de suelo resultarán:
31
𝑢𝐴
ℎ𝐴 = 𝑧𝐴 +
= 𝑧𝐴
ϒ𝑤
𝑢𝐶 𝑢𝐶
ℎ𝐴 = ℎ𝐶 = 𝑧𝐶 + =𝐿+ → 𝑢𝐶 = 𝛥𝐿 · ϒ𝑤
ϒ𝑤 ϒ𝑤
𝐿 + 𝛥𝐿 � 𝑢𝐵 𝑢𝐵
ℎ𝐴 = ℎ𝐵 = 𝑧𝐵 + =0+ → 𝑢𝐵 = (𝐿 + 𝛥𝐿) · ϒ𝑤
ϒ𝑤 ϒ𝑤
las presiones totales verticales:
𝜎𝑣𝐶 = 𝛥𝐿ϒ𝑤
𝜎𝑣𝐵 = 𝛥𝐿ϒ𝑤 +Lϒ𝑠𝑎𝑡
y las tensiones efectivas verticales:
𝜎′𝑣𝐶 = 𝜎𝑣𝐶 − 𝑢𝐶 = 0
𝜎′𝑣𝐵 = 𝜎𝑣𝐵 − 𝑢𝐵 = 𝐿(ϒ𝑠𝑎𝑡 − ϒ𝑤 )
La Figura 3.13b) muestra una situación en la que el recipiente anexo se encuentra a una
altura 𝛥ℎ por encima de la superficie libre del permeámetro, lo que establece una
diferencia de carga o altura piezométrica. Asumiendo que no se producen pérdidas de
carga en los recorridos sin suelo (tramos AB y CD), se tendrá:
𝑢
Punto B: ℎ𝐵 = ℎ𝐴 = 𝑧𝐴 + ϒ 𝐴 = 𝑧𝐴 = 𝐿 + 𝛥𝐿 · 𝛥ℎ
𝑤
𝑢𝐷
Punto C: ℎ𝐶 = ℎ𝐷 = 𝑧𝐷 + ϒ = 𝑧𝐷 = 𝐿 + 𝛥𝐿
𝑤
De esta forma, la diferencia de altura piezométrica es precisamente igual a 𝛥ℎ, el flujo

resulta ascendente en la masa de suelo (ℎ𝐵 > ℎ𝑐 ) y el gradiente hidráulico producido será:
𝑖 = 𝛥ℎ/𝐿.
Otra forma sencilla de comprobar las condiciones de flujo es observar directamente los
piezómetros. En la Figura 2.23b) se aprecia claramente que la cota del piezómetro P1 es
mayor que en el P2, y ésta a su vez mayor que en el P3, lo que indica que ℎ𝑃1 > ℎ𝑃2 > ℎ𝑃3
y que el flujo es ascendente. El gradiente se puede también determinar de forma directa,
32
sin más que leer las alturas del agua de cada piezómetro en la regleta, que se recuerda son
directamente las alturas piezométricas, y dividir por los recorridos del agua entre los
piezómetros, que no son más que las diferencias en la altura geométrica:
ℎ𝑝1 − ℎ𝑝2 ℎ𝑝2 − ℎ𝑝3 ℎ𝑝1 − ℎ𝑝3

𝑖= = =
𝑧2 − 𝑧1 𝑧3 − 𝑧2 𝑧3 − 𝑧1
Con respecto a las presiones intersticiales en los puntos extremos de la masa de suelo se
tendrá:
𝑢𝐵 𝑢𝐵
ℎ𝐵 = ℎ𝐴 = 𝐿 + 𝛥𝐿 + 𝛥ℎ = 𝑧𝐵 + = 0+ → 𝑢𝐵 = (𝐿 + 𝛥𝐿 + 𝛥ℎ) · ϒ𝑤
ϒ𝑤 ϒ𝑤
𝑢𝐶 𝑢𝐶
ℎ𝐶 = ℎ𝐷 = 𝐿 + 𝛥𝐿 = 𝑧𝐶 + =𝐿+ → 𝑢𝐶 = 𝛥𝐿 · ϒ𝑤
ϒ𝑤 ϒ𝑤
Como se deduce de las relaciones anteriores y de las columnas piezométricas en la

situación de flujo ascendente, las presiones intersticiales en la masa de suelo resultan
superiores a las de la condición hidrostática. Habida cuenta de que las tensiones totales
verticales no han sufrido variación alguna (se conserva la misma altura de suelo saturado
sobre cada punto y la misma lámina de agua CD), las tensiones efectivas verticales habrán
disminuido. Así, en el punto B se tendrá:
𝜎′𝑣𝐵 = 𝜎𝑣𝐵 − 𝑢𝐵 = (𝛥𝐿 · ϒ𝑤 + 𝐿 · ϒ𝑠𝑎𝑡 ) − (𝐿 + 𝛥𝐿 + 𝛥ℎ) · ϒ𝑤

𝜎′𝑣𝐵 = 𝐿(ϒ𝑠𝑎𝑡 − ϒ𝑤 ) − 𝛥ℎ · ϒ𝑤
La expresión anterior sugiere que si se aumenta lo suficiente la diferencia de carga 𝛥ℎ se

podrían llegar a anular las tensiones efectivas del suelo, situación que se conoce como
sifonamiento. En estas condiciones un suelo sin cohesión pierde completamente su
resistencia al corte y pasa a comportarse como un fluido. Un ejemplo típico de este caso
son las arenas movedizas. La expresión anterior se puede formular en función del gradiente
hidráulico (𝑖 = 𝛥ℎ/𝐿):
𝜎′𝑣𝐵 = 𝐿(ϒ𝑠𝑎𝑡 − ϒ𝑤 ) − 𝑖 · 𝐿 · ϒ𝑤 = 𝐿 · (ϒ𝑠𝑎𝑡 − ϒ𝑤 − 𝑖 · ϒ𝑤 ),
33
de forma que el sifonamiento se alcanzaría para un gradiente determinado 𝑖𝐶 , llamado

gradiente crítico:
ϒ𝑠𝑎𝑡 − ϒ𝑤
ϒ𝑠𝑎𝑡 − ϒ𝑤 − 𝑖 · ϒ𝑤 = 0 → 𝑖𝐶 =
ϒ𝑤
Si se tiene en cuenta que un orden de magnitud habitual para el peso específico saturado de
un suelo es ϒ𝑠𝑎𝑡 = 20 𝑘𝑁/𝑚3 y que el peso específico del agua es próxima a ϒ𝑤 =10 𝑘𝑁/
𝑚3 , el gradiente crítico suele encontrarse en torno a 𝑖𝐶 = 1.
En la Figura 3.13c) muestra una tercera alternativa de flujo, en la que la lámina de agua del
recipiente anexo se encuentra por debajo de la lámina de agua del permeámetro. En estas
circunstancias la diferencia de altura piezométrica 𝛥ℎ, pero en esta ocasión el flujo resulta
descendente en la masa de suelo (ℎ𝐶 > ℎ𝐵 ) con gradiente hidráulico (𝑖 = 𝛥ℎ/𝐿).
De nuevo, observando los piezómetros se aprecia directamente que la cota de agua en el
piezómetro P3 es mayor que en el P2, y ésta a su vez mayor que en el P1, lo que indica que
ℎ𝑃3 > ℎ𝑃2 > ℎ𝑃1 y que el flujo es descendente. El gradiente se puede también determinar
de forma directa a partir de los piezómetros.
Con la relación de las leyes de tensiones en los puntos extremos de la masa de suelo, con
respecto a las presiones intersticiales, se tendrá:
𝑢𝐵 𝑢𝐵
ℎ𝐵 = ℎ𝐴 = 𝐿 + 𝛥𝐿 − 𝛥ℎ = 𝑧𝐵 + = 0+ → 𝑢𝐵 = (𝐿 + 𝛥𝐿 − 𝛥ℎ) · ϒ𝑤
ϒ𝑤 ϒ𝑤
𝑢𝐶 𝑢𝐶
ℎ𝐶 = ℎ𝐷 = 𝐿 + 𝛥𝐿 = 𝑧𝐶 + =𝐿+ → 𝑢𝐶 = 𝛥𝐿 · ϒ𝑤
ϒ𝑤 ϒ𝑤
Las presiones intersticiales en la masa de suelo resultan inferiores a las de la condición

hidrostática y por tanto las tensiones efectivas verticales habrán aumentado. Así en el
punto B se tendrá:
𝜎′𝑣𝐵 = 𝜎𝑣𝐵 − 𝑢𝐵 = (𝛥𝐿 · ϒ𝑤 + 𝐿 · ϒ𝑠𝑎𝑡 ) − (𝐿 + 𝛥𝐿 − 𝛥ℎ) · ϒ𝑤

𝜎′𝑣𝐵 = 𝐿(ϒ𝑠𝑎𝑡 − ϒ𝑤 ) − 𝛥ℎ · ϒ𝑤
34
4 EL PROGRAMA MODFLOW
4.1 Descripción general
Modflow es un programa de modelación tridimensional para flujo de agua subterránea que

utiliza una expresión de diferencias finitas de aproximación. Fue desarrollado por el
Servicio Geológico de los Estados Unidos (USGS), bajo el sistema de Fortran 88 y
documentado en los años 1995 – 1996. El programa fue generado por un proceso
matemático y matrices complejas bajo la ley de Darcy y sus resultados finales presentaban
una calidad muy pobre. Posteriormente la Universidad de Waterloo logró codificar un
módulo para una mejor calidad visual en los resultados el cual le dieron el nombre de
Visual Modflow. Hoy en día se han estructurado diferentes módulos codificados para otras
aplicaciones particulares y otras nuevas compatibilidades.
Modflow puede simular un fluido en estado transitorio y estacionario, combinando varias

tipos de capas, así como tipo de acuíferos confinados y semiconfinados, también puede
calcular el movimiento de una partícula, es decir modelos de transportes de contaminantes.
El modelo puede requerir de varios parámetros, depende del cual sea su aplicación requiere
datos como: pozos de extracción o de inyección, área de recarga, flujos y drenes,
conductividad hidráulica, conductividad hidráulica vertical, topografía, flujos de entrada y
salida, coeficiente de almacenamiento, condiciones de fronteras, coordenadas para capas
horizontales, geometría del acuífero, etc.
Los modelos son, por definición, aproximaciones y simplificaciones de un sistema natural.

Las suposiciones y las soluciones de compromiso debidas a las limitaciones en la
capacidad de cómputo, la precisión limitada del proceso de solución, la escasez de fondos
y la insuficiencia de los datos de entrada, pueden combinar los resultados y las
predicciones.
35
Las limitaciones en Modflow se basan en que los datos de entrada afectan la precisión y
aplicabilidad de los modelos. Una importante información pudiera ser la arquitectura del
modelo: el tamaño de la malla y el número de capas. Los primeros modelos matemáticos
estaban limitados por la capacidad de cómputo y la poca disponibilidad de información
hidrogeológica; el tamaño de la malla era grueso y los datos eran limitados. Hoy en día, el
tamaño de las celdas puede ser mucho más pequeño debido a que los sistemas de cómputo
se encuentran muy avanzados tanto en su capacidad como en memoria.
Sin embargo, la información hidrogeológica no ha aumentado a la misma velocidad que la

capacidad de simular la presencia y el movimiento del agua subterránea utilizando celdas
muy pequeñas. Actualmente el tamaño de la celda está limitado más bien por el objetivo
del modelo y los datos hidrogeológicos, y no por el tiempo de cómputo.
Los objetivos de un modelo matemático tomando en cuenta la conservación del acuífero

son principalmente tres:
1) Implementar un modelo de simulación hidrodinámica del acuífero para hacer

balances.
2) Elaborar escenarios de manejo del acuífero.
3) Hacer recomendaciones para la administración del agua subterránea, operando
los escenarios de incremento paulatino de extracción.
Modflow resuelve las ecuaciones diferenciales parciales usando el método de diferencias

finitas, donde cada celda de Modflow es un volumen unitario.
En general la ecuación de flujo de aguas subterráneas que emplea Modflow es la siguiente:
𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕 𝜕ℎ 𝜕ℎ
�𝐾𝑥𝑥 � + �𝐾𝑦𝑦 � + �𝐾𝑧𝑧 � + 𝑊 = 𝑆𝑆
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑡
donde: Kxx, Kyy y Kzz = Valores de la conductividad hidráulica para cada eje (L/T).
h = Es la pérdida de carga hidráulica.
36
W = Flujo volumétrico por unidad de volumen representada como el

suministro o descarga de agua, donde los valores negativos indican extracción
de agua y los positivos inyección de agua (𝐿¯¹).
Ss = Almacenamiento específico del medio poroso (𝐿¯¹).
t = Tiempo (T).
Diferencias finitas: La forma de la diferencial parcial por diferencias finitas en un espacio

discretizado del dominio del acuífero representado por filas, columnas y capas es:
𝑚 𝑚 𝑚 𝑚 𝑚 𝑚
𝐶𝑅𝑖,𝑗−1,𝑘 �ℎ𝑖,𝑗−1,𝑘 − ℎ𝑖,𝑗,𝑘 � + 𝐶𝑅𝑖,𝑗+1,𝑘 �ℎ𝑖,𝑗+1,𝑘 − ℎ𝑖,𝑗,𝑘 � + 𝐶𝐶𝑖−1,𝑗,𝑘 �ℎ𝑖−1,𝑗,𝑘 − ℎ𝑖,𝑗,𝑘 �
2 2 2
𝑚 𝑚 𝑚 𝑚
+ 𝐶𝐶𝑖+1,𝑗,𝑘 �ℎ𝑖+1,𝑗,𝑘 − ℎ𝑖,𝑗,𝑘 � + 𝐶𝑉𝑖,𝑗,𝑘−1 �ℎ𝑖,𝑗,𝑘−1 − ℎ𝑖,𝑗,𝑘 �
2 2
𝑚 𝑚 𝑚
+ 𝐶𝑉𝑖,𝑗,𝑘+1,𝑘 �ℎ𝑖,𝑗,𝑘+1,𝑘 − ℎ𝑖,𝑗,𝑘 � + 𝑃𝑖,𝑗,𝑘 ℎ𝑖,𝑗,𝑘 + 𝑄𝑖,𝑗,𝑘
2
𝑚 𝑚−1
ℎ𝑖,𝑗,𝑘 − ℎ𝑖,𝑗,𝑘
= 𝑆𝑆𝑖,𝑗,𝑘 (∆𝑟𝑗 ∆𝑐𝑖 ∆𝑣𝑘 ) 𝑚
𝑡 − 𝑡 𝑚−1
donde: hi,j,km = Pérdida de carga hidráulica en la celda i,j,k al paso del tiempo m.
CV, CR y CC = Conductancia hidráulica, o un pedazo de conductancias
entre los nodos i,j,k y un nodo vecino *
Pi,j,k es la suma de los coeficientes de la pérdida de carga de las fuentes y de
las descargas.
Qi,j,k = Suma de las constantes de los términos de las fuentes y las descargas,
cuando Qi,j,k< 0.0 es el flujo del sistema de agua subterránea (como el
bombeo) y Qi,j,k> 0.0 es el flujo en superficie (como la inyección).
tm = en el paso del tiempo m.
4.2 Interfaz de comunicación con modflow: Visual Modflow Classic
Para crear un modelo, lo primero que se debe hacer es crear una carpeta donde vayamos a
guardar los archivos que se generan de cada modelización. El programa crea unos 40
37
archivos de los cuales va haciendo uso conforme los va necesitando según modelizaciones
del mismo escenario. Si alguna vez queremos crear nuevos escenarios a partir de uno
existente, debemos crear una copia de la carpeta y trabajar con ella cada vez que se desee.
De esta forma nos aseguramos de que siempre tendremos un archivo base al que recurrir
(contenido en la carpeta inicial) para poder recuperar las simulaciones y los datos iniciales.
La metodología usada por el programa para trabajar con el modelo creado, consiste en
generar archivos a partir de unos valores de ‘input’ (datos de inicio como: geometría,
parámetros físicos, parámetros químicos, tiempos…), de ‘run’ (ejecución) y de ‘output’
(resultados de la modelización).
Más adelante, en el apartado 5.2, se explicará cómo funciona cada una de estas partes
mediante la realización numérica del ejemplo propuesto.
38
5 APLICACIONES
El problema a resolver consiste en una excavación hasta 6m de profundidad en un depósito

de arcilla de 8m de espesor situado sobre 11m de arena gravosa, como lo muestra la Figura
5.1a. La excavación tendrá un área a nivel del terreno de 84mx64m. El estudio del
subsuelo reveló que el nivel del terreno, y los cálculos de estabilidad correspondientes
indicaron que éste debe reducirse 7m (hasta los -5,5m, si contamos desde la cota 0) en el
área de la excavación para tener la seguridad adecuada para que no falle por levantamiento
de la base.
Se realizó un ensayo de bombeo en un pozo de 0,5m de diámetro que penetró

completamente en la arena en el centro de la excavación propuesta y registró descensos del
nivel piezométrico de 1,2m y 0,43m en pozos de observación situados a 12m y 150m del
pozo de bombeo, respectivamente, cuando el caudal estacionario en el pozo era de
11,4 x 10−2 𝑚3/s.
A partir de estos datos tendremos que determinar la disposición adecuada de pozos de

alivio de 0,5m de diámetro que penetran completamente en la arena para alcanzar la
reducción necesaria del nivel piezométrico. Suponer que las bombas sumergibles pueden
utilizarse con una capacidad de operación de 6,5 𝑥 10−2 𝑚3 /𝑠, y que el radio de influencia
del sistema de pozos es de 1400m.
Figura 5.1a. Perfil de la excavación sobre acuífero confinado (Abatimiento con pozo múltiple)
39
Figura 5.1b. Planta de la excavación sobre acuífero confinado (Abatimiento con pozo múltiple)
5.1 Solución analítica
El primer problema que se nos plantea es obtener la permeabilidad, k, que tiene la capa de
arenas. Si se toma la base de la arena como nivel de referencia, ℎ𝑒 = 20,5𝑚. De este
modo, sabemos que para 𝑄 = 11,4 𝑥 10−2 𝑚3 /𝑠, los resultados del ensayo de bombeo son:
En 𝑟1 = 12𝑚, abatimiento = 1,20m, ℎ1 = 20,5 − 1,2 = 19,3𝑚

En 𝑟2 = 150𝑚, abatimiento = 0,43m, ℎ2 = 20,5 − 0,43 = 20,07𝑚
Y así se puede hallar el valor de k:
𝑄 · ln 𝑟2 /𝑟1 11,4 𝑥 10−2 𝑥 ln 150/12

𝑘= = = 0,541𝑥 10−2 𝑚/𝑠
2𝜋D(ℎ2 − ℎ1 ) 2 𝑥 𝜋 𝑥11 𝑥 (20,07 − 19,30)
Una vez que tenemos este dato, se puede continuar con el cálculo del abatimiento. Si los
pozos se instalan a 2m del borde del talud, el sistema encerrará un área rectangular de 88m
x 68m. Para determinar el número de pozos necesario, dicha área rectangular puede
representarse como un pozo circular equivalente de radio 𝑅𝑤 dado por:
88 𝑥 68
𝑅𝑤 = � = 43,6𝑚
𝜋
Para un abatimiento de 7m, el caudal en este pozo equivalente está dado por:
40
2𝜋𝑘D(ℎ𝑒 − ℎ𝑤 ) 2 𝑥 𝜋 𝑥 0,541 𝑥 10−2 𝑥 11 𝑥 7

𝑄= = = 75,4 𝑥 10−2 𝑚3 /𝑠
ln 𝑟𝑒 /𝑅𝑤 ln 1400/43,6
Esta es una estimación de la cantidad de agua que deberá evacuarse de la arena; si se tiene
en cuenta que la capacidad de las bombas es de 6,5 𝑥 10−2 𝑚3 /𝑠 se necesitarán alrededor
de 12 bombas. Como el flujo proveniente de una fuente circular, entonces los 12 pozos
pueden localizarse a distancias iguales alrededor del perímetro del área de la excavación,
como podemos ver en la Figura 5,1b. Ahora puede verificarse que la disposición propuesta
es adecuada calculando el abatimiento en el centro de la excavación propuesta y la cabeza
total en la arena en los pozos mediante la ecuación:
𝑛
𝑄 𝑟𝑒
ℎ𝑝 = ℎ𝑒 − � ln
2𝜋𝑘D 𝑟𝑗
𝑗=1
Donde Q es la capacidad de operación de las bombas, en este caso de 6,5 𝑥 10−2 𝑚3 /𝑠 y 𝑟𝑗

es la distancia desde P hasta el centro del pozo j, que se obtiene a partir de la Figura 5,1b.
En la siguiente tabla tenemos los cálculos para el centro de la excavación y para el pozo 3.
Número de Centro excavación Pozo 3

pozo rj (m) ln re/rj rj (m) ln re/rj
1 44 3,46 38,5 3,59
2 51,1 3,31 19,7 4,26
3 42,8 3,49 0,25 8,63
4 34 3,72 26 3,99
5 42,8 3,49 52 3,29
6 51,1 3,31 70,5 2,99
7 44 3,46 77,8 2,89
8 51,1 3,31 92,2 2,72
9 42,8 3,49 85,6 2,79
10 34 3,72 72,8 2,96
11 42,8 3,49 68 3,02
12 51,1 3,31 62,6 3,11
Σ 41,55 Σ 44,25
Tabla 5.1. Cálculos ln 𝑟𝑒 /𝑟𝑗 en el centro de la excavación y pozo 3

Por tanto, en el centro de la excavación propuesta se tiene:
41
6,5 𝑥 10−2 𝑥 41,56

ℎ𝑝 = 20,5 − = 13,28𝑚
2 𝑥 𝜋 𝑥 0,541 𝑥 10−2 𝑥 11
Entonces el abatimiento = 7,22m. (es favorable ya que el requerido era 7m).

En el pozo 3:
6,5 𝑥 10−2 𝑥 44,24
ℎ𝑝 = 20,5 − = 12,81𝑚
2 𝑥 𝜋 𝑥 0,541 𝑥 10−2 𝑥 11
Por simetría del sistema de pozos ésta es también la cabeza total en la arena en los pozos 5,
9 y 11. De manera similar, puede demostrarse que la cabeza total en los pozos restantes
está en el rango de 12,72m a 12,94m. Se observa entonces que si los pozos están
localizados correctamente la cabeza total es casi la misma. La cabeza total en este caso es
más o menos 12,8m, la cual es satisfactoria y proporciona una profundidad adecuada para
la operación correcta de las bombas sumergibles.
Estas unidades están medidas tomando como cota 0, el nivel inferior de la capa de gravas.
Si tomásemos como cota 0 el nivel superior de la capa de arcillas, el nivel piezométrico en
el centro de la excavación sería:
ℎ𝑝 = 19 − 13,28 = 5,72𝑚
Y el nivel en el pozo 3:
ℎ𝑝 = 19 − 12,81 = 6,19𝑚
42
5.1.1 Comprobación de sifonamiento
Para probar la utilidad del abatimiento, se va a explicar con en nuestro ejemplo como se
produciría sifonamiento al realizar la excavación, si no fuese porque se abate el perfil
piezométrico hasta la altura indicada.
Si consideramos que el peso específico de la capa de arcillas es ϒ𝑎𝑟𝑐𝑖𝑙𝑙𝑎 = 17 𝑡𝑛/𝑚2 , y el

peso del agua ϒ𝑎𝑔𝑢𝑎 = 9,81 𝑡𝑛/𝑚2 , vamos a calcular la tensión total antes de producirse
la excavación, si se produjese la excavación sin abatimiento, y si se hace la excavación
después del abatimiento.
- Antes de producirse la excavación (caso 1), Figura 5.2a:
El nivel piezométrico del acuífero se encuentra en 1,5m, como se muestra en la

Calculamos la tensión efectiva en el punto A:
ϒ𝐴′ = ϒ 𝑇 − 𝑢𝐴 = 136 − 93,2 = 42,8 𝑡𝑛/𝑚

ϒ 𝑇 = 𝐻𝑎𝑟𝑐𝑖𝑙𝑙𝑎 · ϒ𝑎𝑟𝑐𝑖𝑙𝑙𝑎 = 8 · 17 = 136 𝑡𝑛/𝑚
𝑢𝐴 = (1,5 + 𝐻𝑎𝑟𝑐𝑖𝑙𝑙𝑎 ) · ϒ𝑎𝑔𝑢𝑎 = (1,5 + 8) · 9,81 = 93,2 𝑡𝑛/𝑚
𝑡𝑛
Como ϒ𝐴′ = 42,8 𝑚 > 0, no se produce sifonamiento.
Figura 5.2a. Perfil del terreno (caso 1)
43
- Si se produjese la excavación sin abatir el nivel piezométrico (caso 2), Figura 5.2b:
El nivel piezométrico del acuífero continúa en 1,5m. Calculamos la tensión efectiva en el

punto A:
′
ϒ𝐴2 = ϒ 𝑇 − 𝑢𝐴 = 34 − 93,2 = −59,2 𝑡𝑛/𝑚
ϒ 𝑇2 = 𝐻𝑎𝑟𝑐𝑖𝑙𝑙𝑎 · ϒ𝑎𝑟𝑐𝑖𝑙𝑙𝑎 = 2 · 17 = 34 𝑡𝑛/𝑚
𝑢𝐴2 = 𝑢𝐴 = 93,2 𝑡𝑛/𝑚
𝑡𝑛
Como ϒ𝐴′ = −59,2 𝑚 < 0, se produciría sifonamiento, y el terreno colapsaría
Figura 5.2b. Perfil del terreno (caso 2)
- Si se hace la excavación después de abatir el nivel piezométrico (caso 3), Figura 5.2c:
En este caso el nivel piezométrico del acuífero ya ha sido abatido y se encuentra a -

5,5m. Calculamos la nueva tensión efectiva en el punto A:
′
ϒ𝐴3 = ϒ 𝑇 − 𝑢𝐴 = 34 − 24,525 = 9,475 𝑡𝑛/𝑚
ϒ 𝑇3 = 𝐻𝑎𝑟𝑐𝑖𝑙𝑙𝑎 · ϒ𝑎𝑟𝑐𝑖𝑙𝑙𝑎 = 2 · 17 = 34 𝑡𝑛/𝑚
𝑢𝐴3 = 𝐻𝑎𝑟𝑐𝑖𝑙𝑙𝑎 · ϒ𝑎𝑔𝑢𝑎 = 2,5 · 9,81 = 24,525 𝑡𝑛/𝑚
𝑡𝑛
Como ϒ𝐴′ = 9,475 𝑚 < 0, no se produce fallo por sifonamiento.
44
Figura 5.2c. Perfil del terreno (caso 3)
Podemos concluir a la vista de estos resultados, que el abatimiento del nivel piezométrico
es clave en esta excavación, ya que de no hacerlo se produciría un colapso por fallo del
terreno debido al fenómeno de sifonamiento.
45
5.2 Solución numérica
Para comenzar, se crea el archivo Mediante la opción ‘New’ que se encuentra en la pestaña
‘File’, como se muestra en la Figura 5.3, y donde se abrirá una nueva ventana en la que se
asignará un nombre a nuestro archivo y seleccionará la carpeta donde se quiera guardar.
Figura 5.3. Creación del archivo
Una vez que tenemos el nuevo archivo, se debe introducir las condiciones del modelo de
trabajo. Aparecen las siguientes ventanas:
• ‘Project Outline’: Es la primera ventana que aparece al crear el archivo. Se dejará las
diferentes variables como aparecen en la Figura 5.4. El tipo de flujo saturado con
densidad constante, ya que consideramos nuestro acuífero saturado y homogéneo.
Figura 5.4. Project outline
46
• ‘Flow option’: En las opciones de tiempo, se asignará al bombeo un flujo transitorio

‘Transient flow’ de manera que se represente el estudio del flujo a lo largo del tiempo,
y no sólo respecto al resultado final ‘Steady-State Flow’. Además al elegir el tiempo
que tardará en llegar al régimen estacionario se cogerá suficientemente elevado,
concretamente de 3650 días (10 años), para que de esta manera estemos seguros de
que hemos llegado a dicho régimen.
Con esto, se podría introducir la permeabilidad del estrato de arena pero se hará más
adelante, ya que tenemos que diferenciar dos capas con permeabilidades distintas.
Figura 5.5. Flow Option.
• ‘Model domain’: Las dimensiones de la malla se van a diseñar según la distribución

de pozos en el ejercicio. Al tener un radio de influencia de los pozos de 1400 metros, y
una excavación de dimensiones de 88*68 (contando que se pone los pozos a dos
47
metros del perímetro) se otorgará a nuestra malla una extensión de 4000m x 4000m,
con 80 celdas para cada dirección, de manera que queden espaciadas cada 50m.
En cuanto a la profundidad se asigna dos celdas ‘Layers’, debido a que es un acuífero

confinado y tiene dos capas con permeabilidades distintas, de 0 a -8m la capa de
arcillas y de -8 a -19m la capa de arenas, obteniendo un espesor total de -19m
(considero el origen 0 en superficie).
Figura 5.6. Model Domain
48
5.2.1 Creación del modelo (‘Input’)
Una vez introducida la forma de trabajo, se va a ajustar e introducir los datos del ejercicio,
de manera que se pueda obtener los resultados de forma numérica. Para ello, se clica en la
pestaña ‘Input’.
En la paleta de herramientas de la izquierda, hay diferentes opciones con las que ir

editando los datos de Inicio de nuestro modelo. Los comandos ‘View column’, ‘View Row’
y ‘View Layer’ sirven para ir cambiando las vistas del terreno viendo su planta, alzado y
perfil según se necesiten. Justo debajo se puede abrir un desplegable con las distintas
variables de Inicio que cabría la posibilidad de cambiar en nuestro modelo como aparece
en la Figura 5.7.
Figura 5.7. Paleta de herramientas.
49
Para el modelo propuesto se usarán los siguientes parámetros:
• Mallado (‘Grid’): El mallado ya se había configurado previamente, pero en esta

opción se va a refinar (dividir en 2) aquellas celdas en las que se precise de mayor
detalle, con el fin de aproximarse con mayor rigor a la resolución analítica. En el
ejercicio la zona más interesante e influenciada es la zona perimetral de la excavación
donde se sitúan los pozos. Por ello se va a refinar un perímetro de 1000x1000 para
estar seguros de coger toda la zona.
Figura 5.8. Mallado
• Condiciones iniciales de piezometría (‘Initial heads’): Se sitúa el nivel piezométrico

inicial en la capa inferior a la cota 1,5m, como se indica en el enunciado del modelo.
Para ello se selecciona la opción ‘Assing’ y con ‘Window’, y se marca toda la malla en
planta de la capa inferior. En el desplegable que aparece se añade una nueva
condición, donde se marca los 1,5m correspondientes. Si se hace click en ‘Database’
se ven todas las piezometrías, Figura 5,9.
50
Figura 5.9. Piezometrías iniciales (‘Initial heads’)
• Conductividad (‘Conductivity’): En este apartado hay que darle a cada capa su

conductividad hidráulica. En el caso de la capa de arenas su conductividad será
K = 0,541 x 10−2 en los ejes x e y, mientras que en el eje z será 10 veces más
pequeña. Para el caso de la capa de arcillas que es prácticamente impermeable, se
escoge un valor de un rango de valores de arcillas de baja plasticidad (obtenido del
libro de Gonzalez Vallejo 2002) situado entre 10−5 a 10−8 . Por tanto, se elige para su
conductividad K = 10−7 al igual que antes en los ejes x e y, mientras que en el z, es
10 veces menor.
Figura 5.10. Parámetros de conductividad
51
• Condición de contorno para establecer gradiente hidráulico (‘Constant head’):

Es la representación de llegada al pozo de agua proveniente del acuífero, en este caso
vamos a establecer que la primera capa de arcillas el agua se encuentra a la cota 0, y en la
de arenas que el agua se encuentra a la cota 1,5m a partir de aquella zona donde no se ve
afectada por el radio de influencia de los pozos, que es 1400m, quedando la malla como se
muestra en la Figura 5.11 (derecha).
Figura 5.11. Condiciones de contorno para establecer gradiente hidráulico
• Pozos de extracción (‘Pump Wells’): Al igual que en la solución analítica, hay que
situar 12 pozos de extracción situados a 2m del perímetro de la excavación, como se
indica en la Figura 5.12.
Figura 5.12. Posición de los pozos de extracción
52
Éstos quedarán inscritos en un rectángulo de 88 x 68m. Cada uno de estos pozos

tendrá una capacidad de bombeo de 6.5 𝑥 10−2 m3 /s, o lo que es lo mismo 5616 m3 /
día, el signo de este valor será negativo ya que se trata de una extracción y no de una
recarga de agua. En este apartado se le dirá al programa la profundidad a la que actúan
los pozos, en este caso se escoge de -10m hasta los -19m, para asegurar que solo se
coge de la zona permeable. También se especificará el día en el que comienza y acaba
la extracción, para el cual se va a escoger que empieza a los 365 días (1 año) y
finalizaremos a los 3650 días (10 años). En la opción de ‘Database’ podemos ver las
características de cada uno, Figura 5.13.
Figura 5.13. Pozos de extracción
53
• Pozos de observación: Los pozos de observación se ponen con el fin de representar

donde se encuentra el nivel piezométrico y nos pueden servir de comparación con la
solución analítica. En este caso se van a escoger 5:
P-A: en el centro de la excavación

P-B: en la posición del pozo 3 de extracción.
P-C: en la posición del pozo 6 de extracción.
P-D: en la posición del pozo 9 de extracción.
P-E: en la posición del pozo 12 de extracción.
Todos tienen la rejilla para la medición a una profundidad de 15 metros.
Figura 5.14. Pozos de observación
54
Figura 5.15. Distribución final de los pozos en planta
Figura 5.16. Distribución final de los pozos en alzado (centro excavación)
55
5.2.2 Ejecución del modelo (‘Run’)
Este apartado, sirve para escoger como queremos obtener la solución (cantidad y tipo de
datos a obtener de cada periodo). Para empezar, en la pestaña ‘Modflow2000’, se cogerá la
opción ‘Time steps’, donde se divide el tiempo total de cada periodo en intervalos. En este
caso, se eligen 10 intervalos para el primer periodo (de los 365 primeros días) donde no se
bombea caudal y 20 intervalos para el segundo periodo (del día 366 al 3650) en el que son
más años y sí bombeamos caudal, por lo que necesitaremos más datos de salida. La opción
de ‘Multiplier’ sirve para que los datos salgan de manera logarítmica dando una mayor
frecuencia de datos al principio, algo que nos conviene ya que al comienzo del bombeo es
donde mayor cambio se produce en la piezometría.
Figura 5.17. Time steps options
El siguiente paso también en la pestaña ‘Modflow2000’ se cogerá la opción ‘Output

control’ donde tenemos que seleccionar el tipo de datos que queramos que se represente.
En este caso, se escoge tanto el nivel piezométrico como el abatimiento de todos los
periodos de tiempo en los que lo hemos dividido.
56
Figura 5.18. Resultados a representar
Por último, al clicar en la pestaña ‘Run’ saldrá una ventana donde se seleccionará los
‘Engines to run’ de ‘Modflow 2000’ (para representar lo anteriormente dicho) y Modpath
(para representar el movimiento de las partículas), y se ordena la opción ‘Translate&Run’
para que el programa ejecute los datos introducidos en ‘Input’
Figura 5.19. Selección de ‘Engines’ to run
57
5.2.3 Resultados obtenidos (‘Output’)

Este último apartado es el que nos muestra los resultados obtenidos del modelo creado. El
programa muestra datos de los diferentes periodos de tiempo que se le pidieron en la fase
anterior ‘Run’. Para visualizar los distintos periodos en los que se ha dividido el bombeo,
debemos ir a la opción time, pudiendo pasar al siguiente o al anterior directamente con las
opciones ‘Next’ y ‘Previous’ o abriendo el desplegable completo (como se muestra en la
Figura 5.20) si se hace click en ‘Time’.
Figura 5.20. Vista en planta de toda la superficie (capa arenas)
Como son demasiados datos los que nos da el programa, solamente mostraremos los
resultados finales del estudio del bombeo (a los 3650 días).
58
Figura 5.21. Vista en planta de toda la superficie (capa arenas)
Figura 5.22. Vista en planta de toda la superficie (capa arcillas)
59
Figura 5.23. Vista en planta de la zona de excavación (capa arenas)
Figura 5.24. Vista en planta de la zona de excavavión (capa arcillas)
Figura 5.25. Vista de perfil de toda la superficie (centro excavación, pozo obs. A)
60
Figura 5.26. Vista de perfil de la zona de excavación (centro excavación, pozo obs. A)
Figura 5.27. Vista de perfil de la zona de excavación (pozo de observación B)
Figura 5.28. Movimiento de las partículas (planta)
61
Figura 5.29. Movimiento de las partículas (perfil)
Otra de las funcionalidades que encontramos en la parte de output son los gráficos que
ofrece. Concretamente se va a utilizar el gráfico que relaciona la piezometría ‘Head’ (m)
en función del tiempo, ‘Time’ (s). Para ello iremos a la opción ‘Graphs’, y en ‘Time
Series’ seleccionamos ‘Head’. Aparece un gráfico como el de la Figura 5.30, donde
podemos elegir de que pozos queremos que aparezcan las mediciones. En la Figura 5.30 se
han seleccionado los 5 pozos de observación que habían puesto anteriormente.
Figura 5.30. Gráfica Nivel piezométrico-Tiempo
Además si ampliamos el eje x de la gráfica (tiempo), justo en el tramo donde comienza el

bombeo, Figura 5.31, podemos observar cuantos días se tarda en realizar el abatimiento
que se busca. Se puede apreciar que serán unos 17 o 18 días, desde el comienzo del
bombeo.
62
Figura 5.31. Gráfica Nivel piezométrico-Tiempo (Ampliación eje x)
Continuando con las opciones que se ofrecen de ‘Output’, merece la pena comentar los
resultados que ofrece su funcionalidad 3D. Pulsando la tecla F2 o en la barra de
herramientas inferior la opción 3D, aparecerá una nueva ventana donde podemos trabajar
con las distintas opciones que ofrece.
En la parte superior izquierda de la pantalla están todas las opciones que podemos
modificar, como se muestra en la Figura 5.32 izquierda. Al hacer click en ellas, aparece en
la parte inferior izquierda de la pantalla los parámetros que podemos modificar de dicha
opción, pudiendo hacer que esta aparezca o no en pantalla, o modificando los colores con
los que aparecen entre otras cosas, al igual que aparece en la Figura 5.32 derecha.
Figura 5.32. Opciones y parámetros que pueden modificarse en el modelo 3D
63
Además, como se muestra en la Figura 5.33, se puede modificar la vista haciendo uso de
las pestañas de la parte inferior de la pantalla, pudiendo rotar ‘Rotate’ o mover ‘Shift’ la
posición en cualquiera de los 3 ejes, cambiar la posición de incidencia de la luz ‘light
position’ o modificar el periodo de tiempo ‘Time’ al que pertenece la imagen en pantalla,
según los periodos en que fue dividida anteriormente en la fase de ‘Run’.
Figura 5.33. Pestañas utilizadas para rotar la vista
Utilizando estas herramientas que se han explicado, estamos en disposición para sacar
diferentes capturas de nuestro modelo. Se han elegido las siguientes vistas para aclarar un
poco más, los resultados ‘Output’ obtenidos.
Figura 5.34. Vista en planta del modelo 3D
64
Figura 5.35. Vista en planta ampliada del modelo 3D
Figura 5.36. Vista en perspectiva del modelo 3D
Figura 5.37. Vista en perspectiva ampliada del modelo 3D
65
Figura 5.38. Vista de una sección del modelo 3D
Figura 5.39. Vista de una sección ampliada del modelo 3D
Figura 5.40. Vista en perspectiva de la trayectoria de las partículas del modelo 3D
66
5.3 Comparación entre la solución analítica y numérica y discusión de resultados
Si recordamos la solución que se obtuvo analíticamente, el nivel piezométrico en el centro

de la excavación se abatía 5,72 metros, mientras que en el extremo de la excavación
correspondiente al pozo de extracción 3 (o pozo de observación A), el nivel piezométrico
resultaba 6,19m. Se tendrán en cuenta estos 2 valores ya que son los que hallamos de
manera analítica y podemos comparar.
En el caso del resultado numérico obtenido con Modflow a los 3650 días, es decir cuando
se ha alcanzado el estado estacionario, es de 5,725m en el pozo de observación colocado
en el punto A (centro de la excavación). En cuanto al pozo de observación situado en el
pozo B (pozo de extracción 3) es 5,73m.
Si observamos el resto de pozos de observación que habíamos colocado en nuestro modelo

numérico (Figura 5.15), vemos que tienen valores comprendidos entre los 5,73m y los
6,12m (Figura 5.30) de abatimiento (alcanzando un valor muy parecido al que resultaba en
el modelo analítico para el pozo de observación B). Por lo tanto, se puede decir que se
cumple la condición que nos pedía el enunciado de que debía abatirse al menos 5,5 metros.
La diferencia de piezometría que surge en los valores del pozo de extracción 3 (pozo de
observación B), se deben a que los resultados calculados numéricamente son más precisos
que los analíticos, pudiendo surgir estas diferencias.
Para ver esta comparación de piezometrías visualmente, se han realizado las gráficas de la
Figuras 5.41 y 5.42. En ellas se ve la evolución de la piezometría en el tiempo de los dos
pozos de observación por ambos métodos. En el caso del pozo del centro de la excavación
(Pozo A), se puede ver que la solución analítica (‘P-A observed’), que corresponde a la
línea recta de color rojo, es muy similar a la solución obtenida con Modflow (‘P-A
interpolated’). Sí nos fijamos en el pozo situado en el pozo de extracción 3 (Pozo B), esta
diferencia entre la solución analítica (‘P-B observed’) y la solución numérica (‘P-B
interpolated’) es un tanto mayor como ya hemos comentado anteriormente.
67
Figura 5.41. Gráfica comparativa de la piezometría con respecto al tiempo, calculada de manera
analítica y numérica
Figura 5.42. Gráfica ampliada de la comparativa de la piezometría con respecto al tiempo,

calculada de manera analítica y numérica
68
A continuación se planteará una posible alternativa al caso de estudio que puede resultar de
gran utilidad para cálculos de casos reales. Se trata de ver qué ocurriría si, por ejemplo,
debido a la disponibilidad que haya de bombas en el mercado, éstas tuviesen distinta
capacidad de la que se utilizó para el cálculo, o simplemente porque son las que tiene la
empresa para la extracción.
Si las bombas que se encontrasen disponibles tuviesen una capacidad de extracción de

8 𝑥 10−2 m3 /s (6912 m3 /día), sería mayor a la que calculamos en la solución analítica de
6.5 𝑥 10−2 m3 /s (5616 m3 /día). Con este caudal de extracción mayor, el abatimiento del
nivel piezométrico que se obtiene es de 7,39m en el caso del pozo de observación A
(centro de la excavación) y 7,4m en el pozo de observación B (situado en el pozo de
extracción 3), como se muestra en la Figura 5.43, donde vemos la comparativa entre estos
valores numéricos, con respecto a la restricción de 5,5m.
Figura 5.43. Gráfica ampliada de la comparativa de la piezometría con respecto al tiempo con el
nuevo caudal de extracción de 8 𝑥 10−2 𝑚3 /𝑠 (6912 𝑚3 /𝑑í𝑎)
Figura 5.44. Vista de perfil de la zona de excavación (centro excavación, pozo obs. A), con el
nuevo caudal de extracción de 8 𝑥 10−2 𝑚3 /𝑠 (6912 𝑚3 /𝑑í𝑎)
69
Como solamente es necesario abatirlo hasta los 5,5m para que cumpla con el factor de
seguridad necesario para que el terreno no colapse, y nuestro abatimiento conseguido con
este caudal es mayor, podríamos ahorrarnos la construcción de algunos pozos. Sí
eliminamos los pozos de extracción 1 y 7, obtenemos un abatimiento en el centro de la
excavación y en el pozo de extracción 3, de 5,87m y 5,88m respectivamente (como se
muestra en la Figura 5.44), por lo que cumpliríamos igualmente con la restricción, y
además se ahorraría la colocación de dos pozos y sus respectivas bombas.
Figura 5.45. Gráfica ampliada de la comparativa de la piezometría con respecto al tiempo con el
nuevo caudal de extracción de 8 𝑥 10−2 𝑚3 /𝑠 (6912 𝑚3 /𝑑í𝑎), y eliminando 2 pozos
Figura 5.46. Vista de perfil de la zona de excavación (centro excavación, pozo obs. A) con el nuevo
caudal de extracción de 8 𝑥 10−2 𝑚3 /𝑠 (6912 𝑚3 /𝑑í𝑎), y eliminando 2 pozos
70
6 CONCLUSIONES
El caso de estudio que se ha propuesto consiste en realizar un abatimiento del nivel

piezométrico de un acuífero confinado que se encuentra bajo un depósito de arcillas de 8
metros de espesor y k = 10-7m/s. El acuífero está en una capa inmediatamente inferior de
11m de espesor y de k = 0,541 x 10-2m/s. La excavación tiene unas dimensiones de
84x64m. Mediante la aplicación del método analítico, se obtiene una solución que consiste
en una red de 12 pozos completamente penetrantes dispuestos perimetralmente, con una
caudal de bombeo de 6,5 x 10-2 m3/s cada uno. El descenso en el centro de la excavación es
de 5,72m y en el registro del pozo 3, de 6,19m, cumpliendo con la condición de que el
descenso debe ser superior a 5,5m.
Por otro lado, se ha creado un modelo con Visual Modflow, teniendo en cuenta las
condiciones de contorno planteadas en el problema (dimensiones de la excavación,
permeabilidad, piezometría inicial…) y la solución analítica proporcionada por el método
(disposición de pozos y caudal de extracción de estos). Los resultados obtenidos con el
programa han sido satisfactorios, puesto que proporciona valores de 5,725m en el pozo de
observación A (situado en el centro de la excavación) y 5,73m en el pozo de observación B
(situado en el pozo de extracción 3).
Como se puede comprobar, las diferencias entre ambos métodos son muy leves y, en
cualquier caso, ambas suponen que el nivel piezométrico se sitúe por debajo de los 5,5m de
profundidad,
Además, se ha comprobado que no se produce la condición de sifonamiento, por la cuál el

terreno colapsaría.
Para ampliar información, se han discutido los resultados obtenidos proponiendo otra
solución numérica, para el caso de que las bombas de extracción disponibles en el mercado
fuesen de otro caudal diferente al que se obtuvo por solución analítica. En este nuevo caso,
el caudal es de 8 x 10−2 m3 /s (6912 m3 /día). El abatimiento obtenido con este caudal es
7,39m en el pozo de observación A (situado en el centro de la excavación) y 7,4m en el
pozo de observación B (situado en el pozo de extracción 3). Como es muy superior al
71
necesario de 5,5m que ya incluye un factor de seguridad, se opta por suprimir los pozos de
extracción 1 y 7, que habíamos colocado en la solución inicial. El abatimiento del perfil
piezométrico obtenido con esta nueva disposición de pozos es de 5,87m en el pozo de
observación A (centro de excavación) y 5,88m en el pozo de observación B (pozo de
extracción 3), cumpliendo con la condición de 5,5m que evita el sifonamiento en la
excavación propuesta en este caso.
7 REFERENCIAS
72
- GONZÁLEZ DE VALLEJO, Luis Ignacio. Ingeniería Geológica. Mercedes

Ferrer Gijón, Luis Ortuño Abad, y Carlos Oteo Mazo. Pearson Educación,
Madrid, 2002. ISBN: 84-205-3104-9.
- BERRY, Peter L. & REID, David. Mecánica de suelos. Editado por Martha
Edna Suárez R. McGraw Hill: Santafé de Bogotá, 1993. Traducido de la
primera edición en inglés de An introduction to soil mechanics. ISBN: 958-600-
172-5
- WATERLOO HIDROGEOLOGIC, Inc. Visual Modflow 2011.1 User’s

Manual.
73

Tfe Men Aba

Cargado por

Copyright:

Formatos disponibles

Tfe Men Aba

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tfe Men Aba

Cargado por

Copyright:

Formatos disponibles

Abatimiento del nivel piezométrico en acuíferos confinados

con pozos de bombeo. Comparación entre la solución

Alumno: PABLO JOSÉ MÉNDEZ CÁNOVAS

Este Trabajo Final de Grado de Ingeniería Civil, se centra en la rama de geotecnia y

Así mismo, me gustaría tener en consideración a mis padres José Salvador y

2 OBJETIVOS Y ALCANCE ...................................................................................... 2

3 FUNDAMENTOS TEÓRICOS ................................................................................. 3

3.1 El agua en el medio poroso. Definición de acuífero y parámetros fundamentales....... 3

3.2 Piezometría ...................................................................................................... 7

3.3 La ecuación de Darcy ...................................................................................... 11

3.4 La ecuación general de flujo: Ecuación de Laplace y ecuación de Poisson.............. 13

3.5 Resolución analítica para el diseño de abatimientos mediante ensayos de bombeo... 17

3.5.1 Ensayos de bombeo en régimen permanente. Deducción de las ecuaciones...... 18

3.5.2 Formulación analítica en el diseño de abatimientos de acuíferos confinados.

3.6 Fuerzas de filtración. El fenómeno del sifonamiento ............................................ 29

4 EL PROGRAMA MODFLOW ................................................................................ 35

4.1 Descripción general ......................................................................................... 35

4.2 Interfaz de comunicación con modflow: Visual Modflow Classic.......................... 37

5.1 Solución analítica ............................................................................................ 40

5.1.1 Comprobación de sifonamiento.................................................................. 43

5.2 Solución numérica........................................................................................... 46

5.2.1 Creación del modelo (‘Input’) .................................................................... 49

5.2.2 Ejecución del modelo (‘Run’) .................................................................... 56

5.2.3 Resultados obtenidos (‘Output’) ................................................................. 58

5.3 Comparación entre la solución analítica y numérica y discusión de resultados ........ 67

En este trabajo se aborda un problema de abatimiento del nivel piezométrico de una

En este trabajo se afronta un campo fundamental en ingeniería del terreno, la presencia de

Además, como objetivo secundario, se adquirirán las técnicas, ya sean analíticas o

Pese a la simplificación del escenario planteado, se adquirirán los conocimientos básicos

3.1 El agua en el medio poroso. Definición de acuífero y parámetros fundamentales

- Acuífero libre: es un acuífero no limitado en la parte superior por un techo impermeable,

- Acuífero confinado: es un acuífero limitado en la parte superior por un techo

- Acuífero semiconfinado: es un acuífero de tipo confinado que se diferencia de este en que

Figura 3.1. Tipos de acuíferos

La capacidad de almacenar y transmitir el agua por parte de un acuífero debe ser

- Porosidad: Es la relación entre el volumen de huecos y el volumen total de una roca. Es

𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑙𝑜𝑠 𝑝𝑜𝑟𝑜𝑠

- La porosidad eficaz ne: también denominada porosidad cinemática, referida al volumen

𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑝𝑜𝑟𝑜𝑠 𝑐𝑜𝑛𝑒𝑐𝑡𝑎𝑑𝑜𝑠

- Coeficiente de almacenamiento: Este coeficiente representa la capacidad para liberar agua

- Permeabilidad: La permeabilidad es el parámetro que permite evaluar la capacidad de

donde: k = permeabilidad efectiva [𝐿𝑇¯1 ]

el paso del 10 % en peso de la muestra de material granular. Los coeficientes c, m, n y λ

K = cte * d² y por tanto: K = cte * d10²

En definitiva, la conductividad hidráulica o permeabilidad efectiva k nos da una idea de la

- Transmisividad, T: es el parámetro definido para evaluar la capacidad de transmitir agua

Energía cinética: Ec = ½*m*v²

Energía de presión: Ep = p*V

Donde g representa la intensidad del campo gravitatorio y m la masa o agente sensible al

E = Eh + Ec + Ep = cte (Caso de los fluidos ideales y cuando el régimen es estacionario)

O bien: m*g*z + ½*m*v² + p*V = cte

Siendo ρ la densidad del agua.

En consecuencia, a los términos de la ecuación se les denomina alturas:

Siendo 𝛾 el peso específico del agua.

Al nivel piezométrico se le suele denominar también altura piezométrica, cota piezométrica

Es decir, el nivel piezométrico es el mismo en los puntos A y B.

Figura 3.3. Teorema de Bernoulli

3.3 La ecuación de Darcy

Energía cinética: Ec = ½mv²

O bien: mgz + ½mv² + p*V = cte