EJERCICIO PRACTICO DE MODFLOW
Problema
Se considera el siguiente aluvial de la figura. Se supone que el río es muy caudaloso y
que está bien conectado con el acuífero. Se han medido las siguientes medidas de
piezometría. Hay dos pozos de bombeo. En uno de ellos se ha medido el nivel.
Suponiendo régimen casi-estacionario para el flujo y acuífero confinado, se pide:
a) Plantear las condiciones de contorno del acuífero (para resolver el flujo)
b) Calibrar la Transmisividad del acuífero (se supone homogéneo) en base a los datos
medidos.
c) Dibujar la red de flujo con la Transmisividad calibrada.
d) Realizar el ejercicio cambiando el régimen a transitorio de un año suponiendo que el
caudal extraído de los pozos haya cambiado con los siguientes datos:
1. Periodo: 90 días Q=2l/s Q=5l/s
2. Periodo: 90 días Q=4l/s Q=10l/s
3. Periodo: 90 días Q=1l/s Q=3l/s
4. Periodo: 90 días Pozos inactivos
1
a) Plantear las condiciones de contorno del acuífero (para resolver el flujo)
A partir de las siguientes posibilidades de condiciones de contorno:
- Nivel preescrito
- Caudal preescrito
- Condición mixta
Se elije una condición de contorno de nivel preescrito, fijándose el nivel en la parte
superior e inferior del acuífero (marcado en azul en la figura 1. Siendo esta condición
de contorno la mas adecuada para este caso ya que no se dispone de suficiente
información para las otras dos.
En base a la piezometría del modelo conceptual, véase figura 2. manualmente se fija
el nivel en cada celda del modelo numérico correspondiente con el nivel del modelo
conceptual. De esta forma queda el nivel preescrito en la parte superior e inferior del
acuífero.
En el modelo numérico se le debe introducir debidamente situados los niveles de los
10 piezómetros, las dos extracciones del bombeo de los dos pozos Q 1=2 l/s y Q2=5 l/s,
como también la interpolación de los niveles del río, ver figura 1. Todo ello necesario
para una buena obtención de la dirección del flujo.
2
Figura 1. Modelo numérico, donde se observa el nivel preescrito en la parte superior e inferior
marcado en azul, la situación de los piezómetros, de los pozos y los niveles del río interpolados
b) Calibrar la Transmisividad del acuífero (se supone homogéneo) en base a los datos
medidos.
Para calibrar la Transmisividad se hace manualmente, cambiando el valor de la misma
hasta lograr que los niveles (h) queden bien ajustados en la grafica de niveles
calculados y niveles medidos.
Teniendo en cuenta que no se dispone de la geología ni de las cotas de las capas que
conforman el aluvial para tener una idea de rango de valores de Transmisividad. Se
determina la Transmisividad solamente con el comportamiento del flujo a partir de las
líneas equipotenciales representadas por el modelo numérico y con su calibración.
Se observa que para valores superiores de 1000 m2/d el flujo no queda afectado con la
extracción de los pozos que están bombeando, en cambio para valores inferiores a
1000 m2/d el flujo si se ve afectado por estas extracciones. La calibración en la gráfica
de niveles calculados y niveles medidos es parecida para valores superiores de 1000
m2/d y para valores de Transmisividad inferiores a 500 m2/d los niveles medidos y los
niveles calculados por el programa no quedan bien ajustados.
Para este determinado modelo se asume una Transmsividad homogénea en todo el
aluvial y se obtiene un valor de T= 500 m2/d, ver figuras 3 y 4. Teniendo en cuenta que
el bombeo de los pozos afecta al flujo, véase figura 2.
3
Figura 3. Gráfica de niveles calculados y niveles medidos, se observan los diferentes niveles de
los piezómetros que aproximadamente quedan ajustados en la recta, para T=500 m 2/d.
Figura 4. Se observa el modelo numérico una vez calibrada la Transmisividad
c) Dibujar la red de flujo con la Transmisividad calibrada.
Una vez se han introducido las condiciones de contorno y se ha calibrado la
transmisividad se procede a dibujar la red de flujo, formada por las líneas isopiezas y
las líneas de flujo, ver figura 5.
Se observa que a pesar de que algunos niveles no quedan bien ajustados en la gráfica
de calibración del modelo numérico, las líneas de flujo se ajustan perfectamente a las
del modelo conceptual, véase figura 2.
4
Figura 5. Red de flujo del modelo numérico, se observan las isopiezas en color azul y las líneas
de flujo de color rojo.
Para que el modelo de flujo sea correcto es necesario comprobar que el balance de
las entradas y salidas de caudal sea 0, ver figura 6.
Figura 6. Balance de flujo del modelo numérico, se observa que el balance es igual a 0.
d) Realizar el ejercicio cambiando el régimen a transitorio de un año suponiendo que
el caudal extraído de los pozos haya cambiado con los siguientes datos:
1. Periodo: 90 días Q=2l/s Q=5l/s
2. Periodo: 90 días Q=1l/s Q=3l/s
3. Periodo: 90 días Q=4l/s Q=10l/s
5
4. Periodo: 90 días Pozos inactivos
En este caso el régimen del acuífero es transitorio es necesario introducir un nuevo
dato el almacenamiento especifico Ss. En este caso es Ss=0.0001m
Estos serían los resultados obtenidos Figura 7, Figura 8, Figura 9 y Figura 10.
