Modelación Numérica de Flujo

Subterráneo

M Sc Ing Alex Armas Blancas

 Contenido
1. Introducción a los métodos numéricos
2. Conceptos básicos de flujo y transporte
3. Principios de procesos de modelación
4. Desarrollo de modelo conceptual y construcción de modelo numérico
5. Ejemplos de modelación de flujo (software VISUAL MODFLOW)
6. Evaluación de resultados
7. Ejemplos y tutoriales en el Modelo en diferencias finitas.
 Algunas interrogantes
• ¿Que es un modelo?
• ¿Cuál es el propósito de un modelo?
• ¿Que tipos de modelos existen?

Modelación Hidrogeológica - Enfoque de sistemas

Entradas, I(t) Salidas, Q(t)
Ecuación Transformación
Q(t) =  I(t)
Modelo: Es una representación a escala de un sistema o proceso real.

Modelo conceptual: Hipótesis de como funciona o

se comporta un sistema o proceso.

Modelo matemático: Son abstracciones que

representan los procesos en forma de
ecuaciones, las propiedades físicas se
representan como constantes o coeficientes de
las ecuaciones y las medidas del estado o el
potencial del sistema se hace variable.
Un modelo matemático es cualquier conjunto de ecuaciones o estructuras
matemáticas, completo y consistente, que es elaborado para corresponder a alguna
otra entidad (Física, biológica, Psicológica, conceptual, u otro modelo
matemático).
El modelo matemático se utiliza para:
a) Obtener respuestas sobre lo que sucederá en el mundo físico.
b) Influir en las experimentaciones u observaciones posteriores.
c) Promover el progreso y la comprensión conceptual.
d) Auxiliar a la axiomatización de la situación física.
 Modelo Conceptual y Matemático
• Modelo Conceptual
• Idealización e implicancia de las
condiciones iniciales
• Representación descriptiva del
sistema de agua subterránea
que incorpora la interpretación
de las condiciones geológicas y
hidrológicas
• Que procesos son importantes
para el modelo?
• Cuales son las fronteras?
• Que datos necesitan ser
colectados y cuales están
disponibles?.
 Conceptos básicos
• Identificación del Problema Identificación del problema y
• Sistema y elementos a ser modelados descripción
• Relación e interacción entre ellos
• Grado de precisión Datos de
conceptualización
• Conceptualización y Desarollo entrada
• Representación matemática
• Modelo tipo
• Método numérico desarrollo
• Código de computo
• Condiciones iniciales
• Condiciones de frontera Calibración
• Calibración
• Estimación de Parámetros
• Comparación de modelados con simulados Validación y
• Ajuste de parámetros análisis de Sensibilidad
• Análisis de incertidumbre
• Validación
• Usar un grupo de datos independiente al de calibración Aplicación, escenarios
• Comparar resultados del modelo con observados
Presentación
resultados
 Modelo Conceptual y Matemático
Modelo Matemático
• Representación matemática de los procesos hidrológicos,
hidrología subterránea y del transporte de contaminantes
• Utiliza las ecuaciones fundamentales del flujo y
conservación de masa
• Se basa en observaciones reales o extrapolaciones
Tipos de Modelos

Modelo Conceptual

Estático, descripción y cuantificación Dinámico, predictivo, soluciona

general del medio ecuaciones y cuantifica los procesos
 Métodos de Solución
• Analítico
• Uso de formulas simples
• Solución exacta en el punto de Cálculo
• Modelo homogéneo
• Se basa en observaciones reales o extrapolaciones
• Numérico
• Discretizacion espacial y temporal
• Solución aproximada
• Admite heterogeneidades
• Transforman las Ecuaciones en Derivadas Parciales (PDFs) que
gobiernan el flujo en Ecuaciones Diferenciales Ordinarias (ODEs) o
ecuaciones algebraicas para su solución
El Ciclo Hidrológico
 Principales Procesos de un Modelo
• La evapotranspiración
• Escorrentía superficial
• Flujo sub superficial
• Infiltración
• Flujo en zona no saturada
• Las aguas subterráneas (flujo saturado)
• Flujo en río
 Requerimientos de Datos
• Datos de precipitación (el mayor ingreso al modelo)
• Datos de suelo
• Cobertura vegetal (para la ET)
• Serie de caudales históricos (para calibración)
• Datos de evaporación
• Datos de Temperatura (max, min, promedio)
• Modelo de Elevación Digital (DEM)- Red de canales
(transito de río)
• Datos Nivel Freático (flujo en zona saturada)
 Modelación Numérica de Aguas Subterráneas

Modelación numérica: Comparación de

software
Visual Groundwater
Ventajas/desventajas GMS ModelMuse
MODFLOW Vistas
Gratuito No No No Sí
Fácil de aprender a utilizar Sí Sí No No
Capacidades de
Alta Alta Baja Mediana
visualización
Soporte técnico Medio Medio Alto Medio
Compatibilidad con
versiones actuales de Media Alto Alto Medio
MODFLOW

Frecuencia actualizaciones Baja Alta Alta Alta

Integración GIS Media Alta Alta Alta
Compatibilidad con otros
Media Alta Alta Media
códigos: MT3D, MODPATH
 Aspectos en la Calibración de Modelos
• Los modelos son aproximaciones de la realidad; no pueden
representar con precisión los sistemas naturales
• No hay una sola, estadística aceptado o prueba que determina si
un modelo es válido.
• Se requieren ambas comparaciones gráficas y pruebas
estadísticas en la calibración del modelo y validación
• Los modelos no se puede esperar que sean más precisos que los
errores (intervalos de confianza) en la entrada y los datos
observados
 Modelos Numéricos:
Diferencias Finitas y Elementos Finitos
Differencias Finitas Elementos Finitos Diferencias Finitas con
Volúmen de Control

Knight Piésold and Co. , 2010

Gomes et al (2007) Yanacocha FEFLOW Model EIA Modflow USG, AQUAVEO’s images (GMS)
 Comparación de técnicas numéricas de
Modelación
Método Ventajas Limitaciones
Diferencias Finitas ▪ Conserva la masa Red Retangular resulta en discretización
▪ Intuitivo, simple, conocido y comprobado, facil de innecesaria fuera de la zona de interes Dificil
programar adaptación a sistema muy complejos
▪ Más programas disponibles con capas deformes
▪ Open source (Modflow) Pozos extraen agua de toda la celda
Capas, lineas y conulas continuas
Tensor de K sin derivadas cruzadas
Elementos Finitos ▪ Discretización mas eficiente en zonas de interes solamente No tan intuitivo, teoría más complejaConserva
▪ Pozos pueden ser representados como nudos discretos masa solo globalmente
▪ Bien adaptado a geología compleja Errores locales en el balance de aguas pueden
▪ Tensor de K totalmente representado ocurrir (raros en FEFLOW)
▪ Capas no tienen que ser continuas Solución mas compleja, discretización mas
elaborada, requiere software Especializados

▪ Discretización eficiente Poco uso en aguas subterráneas

Diferencias Finitas con ▪ Bien adaptado a geología compleja (adopción?)
Volumen de Control ▪ Tensor de K totalmente representado Solución mas compleja, discretización mas
▪ Capas no tienen que ser continuas elaborada, requiere software
especializados
R
Modelación numérica con Visual
MODFLOW
• Objetivo general:
• Construir un modelo numérico de flujo en régimen
permanente
• Objetivos específicos:
• Comprender el modelo conceptual y replicarlo en un
modelo numérico
• Definir el área de modelación
• Definir la geometría de las capas
• Definir las condiciones de borde
• Asignar parámetros hidráulicos
• Ejecutar modelo y analizar los resultados del modelo
 Versiones de Modflow: Características Principales
Malla Manejo de Solucionadores
Celdas secas
MODFLOW 88 / i,j,k, fijo/continuo Rewetting SIP,SSOR,PCG2
MODFLOW 96
MODFLOW-2000 i,j,k, fijo/continuo Rewetting SIP,SSOR,PCG2, AMG,GMG, PCGN
MODFLOW-2005 i,j,k, fijo/continuo Rewetting SIP,SSOR,PCG2,GMG
i,j,k, fijo/continuo; modelos Rewetting SIP,SSOR,PCG2,GMG
MODFLOW-LGR
anidados padre-hijo
i,j,k, fijo/continuo Pseudo soil GMRES
MODFLOW-NWT (2005) χMD
function
Malla desestructurada, Pseudo soil GMRES
MODFLOW USG (2013) poligonos variables function χMD

i,j,k, fijo/continuo,modelos UZF package PCGN

MODFLOW-OWHM anidados padre-hijo (1DRichards) GMRES
χMD

Obs: VMOD incluye otros solucionadores: WHS y SAMG, descritos más adelante
 Modflow: Procesos y Paquetes
▪ El Código de Modlow esta organizado en Procesos y paquetes
▪ Proceso básico: GWP (Groundwater Flow Process) e incluye datos
básicos de tipos de celdas (BAS6), malla
(DIS), multiplicadores y definición de zonas
y parámetros para el proceso de
sensibilidad.
▪ Dentro de este proceso, hay los varios paquetes
▪ Groundwater Flow packages:
diferentes métodos de
construcción de DF incluyendo
parámetros del acuifero
(BCF,LPF,UPW etc)
▪ Paquetes de condiciones de borde (Rio, Dren, strem, GHB, …)
▪ Solvers: solucionadores de matriz: PCG, SIP, SOR …
▪ Control de Salidas (Output Control):
determina que pasos de tiempo
grabar y en que forma (binario, asci,
etc)

Fuente: Modflow online

guide:http://water.usgs.gov/ogw/modfl
ow-owhm/Guide/index.html
 Modflow: Procesos y Paquetes
▪ Otros Procesos:
▪ Proceso de Observación: información de
observaciones para fines de calibración
▪ Procesos especializados:
▪ SWR (formulación de flujo de aguas superficiales 1D /2D,
basado en la ecuación de Saint Venant con aprox. De ola
difusa.
▪ FARM: proceso que incluye la optimización tecnico-
economica del uso del agua en areas agricolas. Incorpora
aguas superficiales/subterráneas
Fuente: Modflow online guide:http://water.usgs.gov/ogw/modflow-
owhm/Guide/index.html
 Principales características que debe cumplir un modelo de flujo subterráneo:
• Debe ser sencillo.
• De fácil aplicación.
• Tiempos de simulación razonables.
• Entendible para el resto de la comunidad científica.
• Debe representar el modelo conceptual establecido.
 Elementosfinitos: Diferenciasfinitas:

• Los cálculos se hacen en los • Los cálculos se hacen para el centro

• nodos de los elementos de la celda.
Variedad de forma de los • Cada celda son prismas rectangulares.
elementos
• Se pueden tener nodos
intermedios
 Elementosfinitos Diferenciasfinitas:
Ventajas Ventajas
• Discretización más versatil • Conceptualización simple
• Refinamiento de la malla únicamente en • Mas facil de usar
zonas de interes
• Más software disponible
• Las condiciones de contorno se
• Conservación del balance de las masas
representan de forma discreta más
ajustadas a la realidad • Balances de agua locales más faciles de
hacer
• La malla desestructurada permite una
mejor adaptación a las morfologías • Testado mundialmente tanto para flujo
complejas como para transporte
 Elementosfinitos Diferenciasfinitas:
Desventajas Desventajas
• Teoría más compleja • Morfologías limitadas

• Velocidades discontinuas en los • Refinamienmto de la malla ineficiente

elementos • Las condiciones de contorno ocupan la
• La conservación de masa local no esta celda completa
garantizada • La representación de las geometrías
• Difícil de entender a la comunidad no escalonadas es más complicada
especializada • Secado de celdas/rehumectación
inestable
 Diferenciasfinitas Block-centered Finite Difference Grid
Inflow
Mas celdas y más puntos de cálculo

Well
Representation

Elementosfinitos
Finite Element Grid
Menos puntos de cálculo (con una
representación más eficiente)

Well
Representation

River
 Diferenciasfinitas
Las condiciones de controno
se discretizan de forma
escalonada

Elementosfinitos
Mejores ajustes para
representar las condiciones de
contorno
 • Formas irregulares
17,080 celdas/capa • Distribución de los pozos
4,500 nodos/capa
• Numero de puntos de calculo
 Land
a River
Rehumedecer las celdas en Surface

MODFLOW genera
inestabilidades. River

En elementos finitos es posible b

deformar la malla para adaptarse X X X X X X X X X
X X X
X X
la posición de la superficie X X
X X
piezometrica (Free and movible)
Active Inactive
Constant head X Cell which receive
recharge

z
h=z

x
1. Recopilación y análisis de la información.
2. Desarrollo del modelo conceptual
3. Elección del código de simulación
4. Diseño del modelo
a) Dominio del modelo
b) Estructura (Número de capas)
c) Propiedades
d) Condiciones de contorno
5. Calibración en régimen permanente
6. Calibración en régimen transitorio
7. Análisis de sensibilidad Simulaciones de
futuro
Régimen estacionario: Estado de equilibrio del sistema hidrogeológico,
dónde el volumen total de entradas al sistema son iguales al volumen total de salidas.

Entradas = Salidas
Régimen natural: Sistema hidrogeológico que no recibe ninguna influencia
de origen antrópico que modifique el volumen de entradas y/o salidas.
Entradas = Salidas
Régimen influenciado: Conjunto de presiones antrópicas que se ejercen
sobre el sistema hidrogeológico y que provocan una variación en el volumen de
entradas y salidas.

Entradas ≠ Salidas
MODFLOW es un modelo matemático de flujo cuasi-tridimensional que resuelve
la ecuación de flujo mediante el método de diferencias finitas en estado Estacionario
/transitorio.

▪ Desarrollado por McDonald & Harbaugh (USGS) en 1983

▪ Actualizado posteriormente: MODFLOW-88, MODFLOW-96, MODFLOW-2000, MF2005
▪ Calcula el potencial hidráulico en el centro de cada una de las celdas de un malla
cuadrangular y el flujo a través de las caras de los prismas o celdas, resolviendo la
ecuación de flujo

▪ Modeliza el flujo saturado (NO el flujo no saturado)

▪ Medios poros o medios de porosidad equivalente
▪ Asume que la temperatura y la densidad son uniformes
 ∆zk
hi,j,k qi,j+½,k hi,j+1,k

∆yi

∆xj+1/2
Visual MODFLOW es un Interface Gráfico de Usuario (GUI)
que puede ser usado para desarrollar modelos en :

1.MODFLOW modelos de flujo subterráneo

2.MT3D/RT3D modelos de transporte (modflow)
3.SEAWAT modelos de doble densidad
4.WinPEST Módulo de autocalibración.
 Construcción del Modelo con Visual
Modflow
1. Eligir la version del Modflow
2. Definición de la Malla Input
3. Entrada de las Propiedades del Medio Poroso
4. Elección y Entrada de Condiciones de Borde
5. Importar Observaciones
6. Ejecutar la simulación
▪ Eligir el modo de tiempo: permanente o transitorio
▪ Selecionar Solucionadores Run
▪ Discretización temporal
▪ Opciones de rewetting
▪ Otras opciones: condición inicial, anisotropía, recarga, ET
7. Evaluar los resultados (planta, corte, 3D, gráficos de calibración, hidrogramas…) Output
 Construcción del Modelo con Modflow:
Discretización Espacial
▪ Tipos de Discretización en Modflow: (1)

▪ DF Clásica:
▪ Lineas, columnas y capas continuas
▪ Discretización local ineficiente, pero facil de implementar
▪ Local Grid Refinement (LGR)
▪ Lineas, columnas y capas continuas (2)

▪ Modelos anidados padre-hijo, iteración entre padre y hijos

▪ Discretización local eficiente, facil de implementar, muchos archivos
▪ Malla Desestructuradas: Modflow USG
▪ Poligonos variados, discretización local eficiente, simulaciones rapidas
▪ Capas pueden ser discontinues y subdiscretizadas individualmente
▪ Aun en desarrollo (3)
 Mallas Rectangulares Estructuradas

Principios Básicos:
2. Mallas homogeneas son las mejores:
▪ Respetan supuestos de DF
▪ Pero suelen generar muchas celdas
3. Refinar la malla:
▪ Areas con gradientes altos
▪ Alrededor de pozos, minas, etc
▪ Areas con muchas
heterogeneidades
▪ Areas de acuñamientos y capas
truncadas/con cambios topograficos
importantes
▪ Areas com plumas de
contaminación

4. Refinado debe ser gradual, con celdas no

más que 1.5 a 2 veces las dimensiones de
su celda vecina Celdas inactivas
 Mallas Rectangulares Estructuradas vista
en corte
Topografia de las capas:
▪ Anisotropía no-orthogonal
▪ El tensor de K en Modflow no tiene las
derivadas cruzadas y se assume
paralelo a las caras del bloque.

(c) Celdas inactivas (d)

▪ Esto puede generar errores numéricos
en casos como (b) al lado
▪ Dos soluciónes para dicha solución se
presentan abajo:
▪ (c) Grid no deformado (respeta DF; más
celdas, menor error) Celdas
▪ (d) Grid Deformado (en efecto deforma inactivas
la K)
Modificado de Anderson et al, 2014
 Mallas Rectangulares Estructuradas vista
en corte
Topografia de las capas:
▪ Pendientes muy
altos:problemas
▪ (celdas
secas/inestabilidades)

▪ Evitar:
▪ Celdas desconectadas
▪ Overlap de ~75% es
ideal
▪ Refinar ayuda!
 Efectos de la malla

razón de Δy Δy Mejor discretización

aspecto Δx Δx
Flujo de trabajo avanzado: Leapfrog

1. Importar litologias de
pozos

2. Realizar Modelo de
Facies sedimentaria

3. Asignar valores a
cada facie
Source: Schlumberger Water Services
File:Abrir y cerrar problemas.
Importar ficheros modflow.

Input: Entrada de datos del

problema.

Run: Procesar el problema

Output: Resultados del problema

Setup: Selección de código

modflow

Help:Ayuda
Crear un proyecto
nuevo
Abrir un proyecto
existente
MENU
INPUT
Construcción del Modelo con Modflow:
Condiciones de Borde Principales
▪ Condiciones de Borde de Modflow incluyen:
▪ Tipo 1: Carga constante (especificada):
▪ El nivel es especificado en el nudo
▪ Tipo 2: Flujo especificado:
▪ Caudales se especifican en la celda directamente.
▪ Tipo 3: Flujo dependiente de la Carga (Nivel):
▪ El modelo calcula un caudal que es basado en la diferencia de nivel,
multiplicado por una conductancia (Resistencia hidráulica)
MENU
INPUT
MENU
INPUT
 Condiciones de Borde: Rio (RIV)
▪ Condición tipo 3: Hcelda
▪ Modflow calcula el flujo que
Ganador
depende de la diferencia de Hrio
nivel entre el rio y el acuifero Hrio RBot
Perdedor

(celda) k’,b’ Colgado

▪ Q=C (Hrio-Hacuifero) => =
perdedor
C
▪ Q=C (Hacuifero - Hrio) =>
ganador
▪ Q=C (HrioRBot) => rios W
colgados C = k’/b’ (Areario)
L
▪ Conductancia usada para Areario = WL
representar deposición de
finos en el lecho del rio:
 Condiciones de Borde: Lago (LAK3)

▪ Condición tipo 3
▪ Fuente/Sumidero
▪ Cambio de flujo con el modelo
igual que condición tipo Río
▪ Permite añadir al lago:
▪ Precipitación ET, Escurrentía,
bombeo, descargas al lago
▪ Flujos entrantes/salientes de
rios (paquete SFR)
▪ Permite lagos desconectados
▪ Usado para representar minas en
fase de cierre

Anderson, Woessner and Hunt (2014)

Lago en VMOD
Lago en VMOD
• T= 300 d
• Hlago=14.987m
Ejemplo de uso del paquete MNW

• Pozo
abandonado

Konikow and Hornberger, 2006

MENU
INPUT
MENU
INPUT
MENU RUN
MENU RUN
MENU OUTPUT
 Simulaciones
▪ Régimen de la simulación:
▪ Estacionario (Permanente): H=Constante
▪ ‘Quasi-estacionario’ en realidad
▪ Transitorio: H=f(t)
▪ H depende del tiempo
▪ Requiere definir:
▪ Condición Inicial
▪ Variación temporal de condiciones de borde
▪ Discretización temporal:
▪ Periodos de Estress (cambios de H o Q en bordes/pozos)
▪ Pasos de tiempo (Δt)
▪ Control de salidas (Output Control):
▪ Que pasos de tiempo grabar durante la simulación
 Condición Inicial
Q

▪ Regimen Estacionario: Ho = nivel pre-bombeo

(condición inicial)
▪ h no cambia con el tiempo (h=cte.)
▪ Si Q=0 en la figura, h =cte. s1 h1 h2
▪ Si Q≤Recarga, despues de un cierto s Hn = nivel estabilizado com
tiempo => nivel se estabiliza, h=hn bombeo (estacionario)
2
▪ Regimen Transitorio:
▪ h=f(t); se requiere una condición
inicial para iniciar los cálculos
▪ ho = Condición de contorno en el
tiempo
▪ Simulaciones empiezan a partir de
ho :
▪ s1= ho – h1
▪ s2= ho – h2 …
 Condición Inicial
▪ Caso especial Transitorio:
Monitoring Well KW001a

1.5

▪ Equilibrio dinámico 1

▪ Condicion inicial requiere 0.5

Timd (days)
simulacion transitoria por algun 0

tiempo, hasta que se estabelezca -0.5

la condicion dinamica (dias hasta -1

meses) -1.5

▪ Especialmente complejo en 0 0.5 1 1.5

Hydraulic Head (m)
2 2.5 3

modelos de intrusion salida (varios KW001a obs KW001a calc

años, debido a efectos de
densidad).
 Discretización Temporal
Periodos de Estrés (P.E.): Periodo de estrés 1 Periodo de estrés 1

▪ Cualquier cambio en H o Q en
condiciones de borde en el modelo

pozo
Q
= 1 PE
▪ Cada Periodo de Estres debe ser
subsividido en pasos de tiempo
(Δt) , para calcular los niveles Tiempo
correctamente
▪ EL número de pasos de tiempo y
su duración influencian los
resultados!
▪ Default: 10 pasos de tiempo por PE
 Discretización Temporal
Entrada de PE/ Δt en GWVistas:
Pasos de Tiempo:
▪ Modflow requiere:
▪ No. del Periodo de Estres
▪ Duración del P.E.
▪ No. de Pasos de tiempo
▪ Multiplicador del paso de timepo
(TSMULT)
▪ Algunas GUI requieren entradas en PE,
no en tiempos (dias): GWV por
ejemplo

▪ En Visual Modflow se entra Start time

y Stop Time en cada C.Borde y la GUI
prepara PE y Δt automaticamente
Seleccionar modo transitorio en Setup
Discretización temporal
Entrar Condición inicial
Control de salidas: Eligir que grabar
Activación de pozos en VMOD

Activado

Desactivado
Simulación Transitoria: Salidas
Simulación Transitoria: Salidas –
3D
Simulación transitoria: abatimientos (2D)

280 dias 1023 dias

OJO: Necesitas tener una buena condición inicial para que los
CAPA 1
abatimientos estén correctos!
Simulación transitoria: abatimientos (2D)

280 dias 1023 dias

CAPA 3
Simulación transitoria: abatimientos (3D)
 Condiciones de Borde (Límites)
▪ La solución de la ecuación de Flujo require la
definición de condiciones de borde
▪ Condiciones de borde defininen la interacción
entre el modelo y areas externas
▪ Se require condiciones de borde en todos los
limites del modelo.
▪ Si no se especifica nada, el modelo assume
que no hay intercamio de agua entre el
modelo y las areas externas (no Flow o flujo
nulo).
▪ Hay 3 tipos principales de condiciones de
borde.
Predicciones de Abatimientos/Flujos

SCENARIO 1 - 450 L/s - 3 PUMPING WELLS Scenario 1- 450 L/s, 3 Wells

100%

90%
1.4
80%
1.2 70%

Pe rcent Reduction on Flow

1.0 60%
Drawdown (m) 0.8 50%

0.6 40%

30%
0.4
20%
0.2
10%
0.0
0%
0 20 40 60 80 40 60
0 20 80 100
time (years) time (years)
Aplicaciones Prácticas
Los modelo son herramientas prácticas de manejo de acuíferos y nos periten
realizar simulaciones predictivas. Algunos ejemplos de aplicaciones prácticas:
▪ Aplicaciones de Flujo
▪ Explotación sostenible de acuíferos: diseño de campo de pozos,
impactos hidrolóicos, interferencia entre pozos, estimación de
caudales sostenibles etc.
▪ Drenaje de Minas/excavaciones: caudales requeridos, ubicaion y
numero de pozos, drenes,
zanjas de bombeo pantallas impermeables etc.
▪ Impactos hidrologicos: reducción del flujo de base a areas sensibles
▪ Aplicaciones de Transporte
▪ Evaluación de las tasas de migración y dimensiones de las plumas de
contaminación
▪ Diseno de Sistemas de remediación
▪ Evaluación del riesgo a la salud humana
Para Qué Sirven las Zonas de Captura?

▪ Definición de Perímetros de Protección de Pozos (PPP)

▪ De donde viene el agua del pozo?
▪ Cuanto tiempo demora desde la recarga hasta el pozo?
▪ Diseño de sistemas de Remediación
(Descontaminación)
• Cual es la mejor manera de limpiar la contaminación
(cobertura, tasas de extracción, tiempo de tránsito)?
• Sistemas de Bombeo/Tratamiento
• Barreras impermeables
• Barreras permeables (‘funnel and gate’)
• Drenes de extracción
• Sistemas de remediación in situ
• Biodegradación natural o estimulada (inyección/extracción)
• Sistema de Extracción de Vapores del Suelo
Ejemplos: Definición de zonas de captura y
protección con rastreo de partículas

Fuente: Cleary, Princeton remediation course

 ¿2D ó 3D?

• Zonas de captura
mas chicas en 3D,
debido a la
componente 2D 3D
vertical del flujo 1 Layer 8
Layers

Source:
USGS
2004-
5038
2D v/s 3D
• Importante para la remediación de areas

contaminadas!

 600m???

Q = 2-3
gpm/pozo

Source: Princeton Groundwater

Ejemplos: simulación de zonas de captura

Zona
estacionaria
Zona 10 años

Zonas de captura basadas en el

tiempo de tránsito

Source: Frind et al (2009)

Drenaje de Minas

Source: Schlumberger Water Services, 2006

Despresurización de taludes
▪ Bingham Canyon Pit - Utah
▪ Modelo en sección transversal
▪ Modelo de corte con espesor
unitario (1m o 1 pie)
▪ Asume flujo transversal nulo, lo
que es una simplificación no
siempre realista.
▪ Muy util para evaluar los efectos
de heterogeneidades en detalles
▪ FEFLOW – Elementos Finitos
(SEEP/W tambien es una Buena
opción)
▪ Usa elementos unidimensionales
para similar drenes radials
▪ Muy dificil de hacerlo con
Modflow
▪ Notar el efecto de los los tuneles
Source: Schlumberger Water Services
Simulación de Presas y embalses

Source: Hernandez y Gomes, 2003

Simulación de Presas y embalses

Presa de
tierra
Dren
Barrera
Impermeable
 Modflow: Control de la Convergencia
• Depende del solucionador utilizado
• Los más modernos poseen los seguientes
controles más importantes: H inicial
▪ Número máximo de iteraciones (externas):
▪ Default: 50 iteraciones Iteración:
Calcula Nuevo nivel
▪ Si Modflow no llega al criterio de convergencia de
nivel hasta el no. max. de iteraciones, no converge Compara nivel antiguo
▪ Criterio de convergencia de nivel con nivel nuevo

▪ Compara niveles en todas las celdas entre las dos

ultimas iteraciones Δh<Criterio
Sí
Modelo
convergenci
Converg e
▪ Si la diferencia es inferior al criterio de a
(0.01m)
convergencia de niveles: modelo converge
No
▪ Si no, intentara una nueva iteración ¿
No. Mx Modelo
▪ Criterio de convergencia de residuo No Iteracione
s
Sí NO
▪ en algunos solucionadores (PCG, NWT, USG) se Fue Converge
alcanzado
refiere a diferencias de caudales en el modelo ?

(m3/d)
▪ En otros solucionadores es una manera de
terminar las iteraciones internas
 Modflow: Control de la Convergencia
H inicial
▪ Solucionadores tipo Gradiente Conjugado Pre-
condicionado: Iteraciones internas y externas: Iteración:
Calcula Nuevo nivel
▪ Internas: iteración de niveles solamente
▪ Externas: actualización de la matrix de conductancia
Compara nivel antiguo
▪ Dumping Factor: con nivel nuevo
▪ Disminuye la agresividad del solucionador:
▪ 1.0: sin dumping;
Δh<Criterio
▪ 0.6 es el limite más bajo en general convergenci Sí Modelo
▪ Valores demasiado bajos suelen fuerzar la a
(0.01m)
Converge
convergencia, pero potencialmente con soluciones
equivocadas ¿
No
▪ Uso del dumping aumenta el tiempo de la solución, No. Mx
Modelo
No Iteracione
pero disminuye la generación de oscilaciones y de Sí
s NO
Fue
celdas secas alcanzado Converg e
?
 Modflow: Solucionadores
• SIP – Strongly Implicit Procedure: muy estable, pero lento y require alto número de iteraciones
• SSOR – Successive Slice Over Relaxation: inestable en 3D, recomendado para modelos de corte 2D
• PCG2 – Pre-Conditioned Conjugate Gradient: estable, robusto, no muy agresivo, buena opción para minas
• PCGN – PCG- with Non Linear control: similar al PCG2 pero mejor para problemas no-lineales
• WHS – PCG-BISTAB: estable, robusto, más rapido que el PCG, pero agresivo. Usar con dumping (<1.0; no >0.6)
• GMG – Geometric Multi Grid: estable, robusto, rápido, indicado para modelos con muchas celdas. Mejor opción si paquete
LPF es usado.
• SAMG – System Algeabric Multi Grid: estable, robusto, rápido para modelos con muchas celdas y alto contraste de K
• GMRES – PCG method for irregular matrices (Modflow NWT/USG, OWHM): similar a los PCG, muy estable
• χMD – PCG BISTAB for irregular matrices (Modflow NWT/USG, OWHM): similar al WHS, tambien muy estable y rápido
• DE4 – Direct solver: ocupa bastante memoria pero brinda una solución exacta en casos lineales. Para casos no lineales,
tiene que iterar también
 Convergencia
▪ Modelo no Converge, ¿que hacer?
▪ Malla:
▪ Sin contrastes grandes de tamaño de celda (no más que 2 veces la celda vecina).
▪ Use herramienta de suavización (‘Smoothing’)
▪ Capas con celdas conextadas ~ 75% de overlap
▪ Parametros:
▪ Evite altos contrastes de K
▪ use linea de celdas con valores intermediarios entre altaK y baja K
▪ ej: use linea de K= E-4 entre zonas de K= E-7 y K=E-2
▪ Intente con campos interpolados si se puede
▪ Use solucionador GMG o SAMG
▪ Celdas secas:
▪ Verifique si oscilación no es causada por pozos con pequeña rejilla en la capa 1. Busque por ‘WET’ en archivo .LST
(VMOD) o .out (GWV)
▪ Saque la rejilla de la capa 1 si esto pasa
▪ Aumente el umbral seco-húmedo (‘Wet-Dry Threashold’) en opciones de rewetting
▪ Cambie a MODFLOW NWT (en Set-up)
Recomendaciones
▪ Usar PCG y WHs en casos communes
▪ WHS más rapido, pero más agresivo; usar dumping ~ 0.7 o parecido
▪ PCG: el valor default es en unidades de segundos (m3/s), si se usa unidades
de dia, aumentar criterio de residuo para 864 o 8.64
▪ PCG: suele ser sensible al número de iteraciones internas; en caso de
problemas, aumente las iteraciones internas bastante (~100 o más)
▪ Usar GMG en casos más complejos, o con uso del paquete LPF
▪ Usar SAMG en casos de modelos muy grandes y con alta contraste de K
▪ Modflow NWT o USG: requieren los solucionadores GMRES o χMD
▪ Tienen dumping dinámico, aplicado cuando necesario. Usar en caso de
oscilaciones