PRÁCTICAS

MÓDULO 4

Luis Miguel Millán Galante

Octubre 2020
EJERCICIO 1
Módulo 4. Aplicación de modelos de simulación (práctico)

Práctica 1. Construcción de un modelo numérico de flujo sencillo

Suponemos un sector de acuífero isótropo, con bordes norte y sur impermeables y este y oeste de piezo-
metría conocida, tal y como se muestra en el esquema siguiente:

Se asume el siguiente escenario:

• En este acuífero hay dos pozos, cuyas coordenadas locales y caudales de explotación en continuo son
los siguientes:

P1. X: 425 / Y: 360 m / Q:100 l/s

P2. X: 610 / Y: 180 m / Q: 75 l/s

• La permeabilidad del acuífero es de 5 m/d.

• La superficie topográfica es plana, con cota 300 msnm y el acuífero tiene un espesor continuo de 100
msnm.

1. Construye el modelo numérico según la situación descrita, teniendo en cuenta un tamaño de celda de
50x50m y obtén los siguientes resultados:

a. Piezometría con pozos en explotación (isopiezas cada 1 metro)

b. Balance hídrico con pozos en explotación

c. Piezometría con pozos inactivos (isopiezas cada 1 metro)

d. Balance hídrico con pozos inactivos

3
 Time
Time == 365
365
IN
Time
Time == 365
365
OUT
IN
1000000200000030000004000000

OUT
Volume [m^3]

4000000
0

STORAGECONSTANT HEAD TOTAL

3000000
Volume [m^3]
2000000
1000000
0

STORAGE CONSTANT HEAD TOTAL

Mass Balance: MODFLOW

Date: 10/27/2020 EJERCICIO 1 - PREGUNTA 2
 SIN POZOS
OUT ENTRADAS SALIDAS
Volume [m^3] STORAGE 0 0
STORAGE 0 CONSTANT HEAD 4725029 4724683
CONSTANT HEAD 4724683 TOTAL 4725029 4724683
TOTAL 4724683

IN
Volume [m^3]
STORAGE 0
CONSTANT HEAD 4725029
TOTAL 4725029
 Time
Time == 365
365
IN
Time
Time == 365
365
OUT
IN
OUT
2000000 4000000 6000000
Volume [m^3]

6000000
0

STORAGE
CONSTANT HEAD
WELLS TOTAL
Volume [m^3]
4000000
2000000
0

STORAGE CONSTANT HEAD WELLS TOTAL

Mass Balance: MODFLOW

Date: 10/27/2020 EJERCICIO 1 - PREGUNTA 4
 CON POZOS
OUT ENTRADAS SALIDAS
Volume [m^3] STORAGE 0 0
STORAGE 0 CONSTANT HEAD 7372550 1853517
CONSTANT HEAD 1853517 WELLS 0 5518800
WELLS 5518800 TOTAL 7372550 7372317
TOTAL 7372317

IN
Volume [m^3]
STORAGE 0
CONSTANT HEAD 7372550
WELLS 0
TOTAL 7372550
EJERCICIO 2
Simulación numérica de aguas subterráneas con Visual Modflow

Práctica 2. Acuífero aluvial ligado a un río

En un acuífero aluvial, en cuyos bordes norte y sur aflora el basamento impermeable, con una anchura
media de 500 metros, discurre un río por su parte central.

El río tiene 5 metros de ancho, un calado de 30 centímetros y la cota de la lámina de agua coincide con la
cota topográfica del terreno.

La permeabilidad vertical del lecho de sedimentos es de 1E-3 m/d, y tiene un espesor medio de 10 centí-
metros.

Puesto que lo que se quiere estudiar es un tramo del acuífero aluvial se ha tomado una sección del mismo
donde se conoce su cota piezométrica, siendo ésta en la sección oeste de 46 mnsm y en la sección este de 44
msnm.

Se ha realizado un levantamiento topográfico del terreno y se han obtenido las siguientes cotas topográ-
ficas:

4
Módulo 4. Aplicación de modelos de simulación (práctico)

Según la información litológica de la zona y los ensayos de permeabilidad realizados, el acuífero se com-
pone de 3 capas:

• Capa superior, de 15 m de espesor y 1 m/d de permeabilidad en la horizontal y 0.1 m/d en la vertical

• Capa intermedia, de 5 m de espesor, 0.001 m/m de permeabilidad en la horizontal y en la vertical.

• Capa inferior, de espesor variable pero la cota de la base se encuentra a 20 msnm.

El acuífero tiene una recarga de 20 mm/año.

2. Ejecuta un modelo en régimen estacionario para un periodo de simulación de 365, con celdas cua-
dradas de 25 metros de lado y obtén los siguientes resultados:

a. Describir el régimen hidrodinámico (red de flujo y relación río-acuífero)

b. Balance hídrico

3. Ejecuta un modelo en régimen estacionario, suponiendo que hay 4 pozos en explotación, para un
periodo de simulación de 365, con celdas cuadradas de 25 metros de lado y obtén los siguientes re-
sultados:

a. Describir el régimen hidrodinámico (red de flujo y relación río-acuífero)

b. Balance hídrico

5
RIO PERDEDOR
 SIN POZOS
OUT ENTRADAS SALIDAS
Time = 365 STORAGE 0 0
Volume [m^3] CONSTANT HEAD 9250 3592
STORAGE 0 RIVER LEAKS 0 5526
CONSTANT HEAD 9250 TOTAL 9250 9118
RIVER LEAKS 0
TOTAL 9250

IN
Time = 365
Volume [m^3]
STORAGE 0
CONSTANT HEAD 3592
RIVER LEAKS 5526
TOTAL 9118
 Time
Time == 365
365
IN
Time
Time == 365
365
OUT
IN
OUT
8000
Volume [m^3]
4000 6000
2000

8000
0

STORAGE
CONSTANTRIVER
HEAD LEAKAGE
TOTAL
6000
Volume [m^3]
4000
2000
0

STORAGE CONSTANT HEAD RIVER LEAKAGE TOTAL

Mass Balance: MODFLOW

Date: 10/27/2020 EJERCICIO 2 - PREGUNTA 2
EQUILIBRIO RÍO PERDEDOR RÍO GANADOR EQUILIBRIO
 CON POZOS
OUT ENTRADAS SALIDAS
Volume [m^3] STORAGE 0 0
STORAGE 0 CONSTANT HEAD 0 66592
CONSTANT HEAD 0 WELLS 73000 0
WELLS 73000 RIVER LEAKS 0 5526
RIVER LEAKS 0 TOTAL 73000 72118
TOTAL 73000

IN
Volume [m^3]
STORAGE 0
CONSTANT HEAD 66592
WELLS 0
RIVER LEAKS 5526
TOTAL 72118
 Time
Time == 365
365
IN
Time
Time == 365
365
OUT
IN
OUT
60000
Volume [m^3]
20000 40000

60000
0

STORAGE
CONSTANT HEAD
WELLS
RIVER LEAKAGE
TOTAL
Volume [m^3]
40000
20000
0

STORAGE CONSTANT HEAD WELLS RIVER LEAKAGE TOTAL

Mass Balance: MODFLOW

Date: 10/27/2020 EJERCICIO 2 - PREGUNTA 4
EJERCICIO 3
Simulación numérica de aguas subterráneas con Visual Modflow

Práctica 3. Simulación del abatimiento del nivel freático en

excavaciones
Se pretende realizar una excavación para un edificio con planta cuadrada y para la cual es necesario depri-
mir el nivel freático por debajo de la cota de excavación.

La cota del terreno actualmente está a 20 mnsm, y la base del acuífero a 0 msnm y la excavación a realizar
tendrá una profundidad máxima de 10 msnm (se facilita archivo .grd con la topografía y profundidad excava-
ción).

Se supone un nivel constante en todos los bordes del área a simular de 18 mnsm, lo cual es indicativo de
que habrá un flujo constante hacia el interior de la excavación.

La permeabilidad del acuífero es de 15 m/d en la horizontal y de 1.5 m/d en la vertical.

El esquema definido es el siguiente:

• Sección transversal por zona central:

6
Módulo 4. Aplicación de modelos de simulación (práctico)

Se solicita:

Con un mallado de 10 m de lado y un área de 200 x 200 metros:

1. Definir red de pozos de bombeo que logren deprimir el nivel freático para un periodo de simulación
de 365 días (número, ubicación y caudal de bombeo).

2. Documentar con figuras el cono de depresión obtenido tanto en planta como en sección transversal.
Si es posible, también en 3D.

7
well ID x y flow (m3/s)
1 50 110 -11.5
1(10) 70 70 -11.5
1(11) 70 150 -11.5
1(12) 110 150 -11.5
1(2) 90 150 -11.5
1(3) 130 110 -11.5
1(4) 90 70 -11.5
1(5) 50 90 -11.5
1(6) 130 90 -11.5
1(7) 50 130 -11.5
1(8) 131 130 -11.5
1(9) 110 70 -11.5
EJERCICIO 4
Simulación numérica de aguas subterráneas con Visual Modflow

Práctica 4. Simulación flujo en una presa

Se pretende definir la red de flujo subterráneo bajo una presa, así como determinar el caudal de paso sub-
terráneo que se produce de una parte a otra.

Crear 1 fila y 20 columnas, con tamaño de celda de 10 m, 10 capas en la vertical, con base 0 mnsm y cota
superior 100 m, salvo para la parte central cuya cota es de 150 msnm, tal y como se muestra en las figuras
siguientes:

• Vista en planta del mallado:

• Vista en perfil del mallado:

La altura piezométrica en la parte izquierda de la presa es constante, siendo de 150 msnm y en la parte
derecha de 100 msnm.

1. Introducir valores de permeabilidad de Kxy=1 m/d Kz=0.1 m/d y para Kxy=0.01 m/d y obtener:

a. Red de flujo bajo la presa

b. ¿Cuánto agua pasa de un lado a otro de la presa?

2. Introducir valores de permeabilidad de Kxy=0.01 m/d Kz=0.001 m/d

a. ¿Hay alguna variación respecto a los resultados anteriores?

b. En caso afirmativo, descríbelos.

8
APARTADO 1.1

Keq = 0.316 m/d

APARTADO 1.2

Keq = 0. 032 m/d

APARTADO 2

Keq = 0. 003 m/d

 EJERCICIO 4

APARTADO 1.1 APARTADO 1.2 APARTADO 2

kx (m/d) 1 0.01 0.01

ky (m/d) 1 0.01 0.01

kz (m/d) 0.1 0.1 0.001

keq (m/d) 0.316 0.032 0.003

i (m/m) 2.5 2.5 2.5

v (m/d) 0.7906 0.0791 0.0079

A (m2) 1000 1000 1000

Q (m3/d) 790.569 79.057 7.906

Se aprecia que el Caudal que pasa bajo la presa varía de acuerdo con el coeficiente de permeabilidad esquivalente; lógico,
puesto que si el gradiente y el área son los mismos en los tres casos, el caudal será directamente proporcional a la
permeabilidad