Tutorial 4 Mt3d
Tutorial 4 Mt3d
TUTORIAL 4
Visual MODFLOW y MT3D
Modelo de flujo y transporte utilizando el módulo MT3D
En este trabajo se estudia un acuífero de 2050 x 1050 m² de superficie. El acuífero tiene una
conductividad (K) de 20 m/día, un espesor saturado inicial (b) de 50 m, velocidad promedio
lineal (vx ) de 0,5 m/día, porosidad total (n) de 0,2 y un gradiente hidráulico (iH ) de 0,5%. Este
gradiente hidráulico puede ser simulado con dos condiciones de nivel constante (50 m en el
lado izquierdo y 40 m en el lado derecho). El almacenamiento específico (Sy )es de 0,15.
Objetivos:
• Analizar el efecto sobre una pluma de contaminación en el acuífero para una fuente de
de descarga constante o instantánea, conservativa o no conservativa.
• Analizar el efecto de parámetros como la dispersión, porosidad, etc. Por otro lado
analizar el efecto sobre la pluma contaminante, provocado por un pozo de extracción o
por una barrera rocosa.
Se despliega la ventana “Create new model”. En ella se nombra el archivo Transporte, que
automáticamente toma la extensión .vmf.
Paso 3: Posteriormente aparece una ventana que permite definir el dominio del acuífero, la
discretización horizontal y vertical. Además se definen las unidades a usar en el modelo de
simulación.
En la sección “Model Domain” se define el dominio del modelo. En este trabajo se define un
acuífero cuadrado de 2050 x 1050 m². Ingrese los parámetros que pide este item:
En la sección “Units” se ingresan las unidades en que se desea trabajar. Para este caso
ingrese:
Length à meters
Time à days
Conductivity à m/day
Pumping Rate à m³/day
Recharge à mm/year
Concentration à mg/L
Entonces se despliega en la ventana principal el acuífero a estudiar con la grilla que se definió
en el paso 3.
Aparece una ventana que pregunta si se desea guardar las propiedades existentes: “Save
property data before existing?”
Se despliega una nueva ventana con las propiedades hidrogeológicas del acuífero. “Default
Property Values”. En esta se deben ingresar los siguientes valores:
Kx [m/d] à 20
Ky [m/d] à 20
Kz [m/d] à 0
Ss [1/m] à 0.003
Sy [ - ] à 0.15
Eff. Por. [ - ] à 0.15
Tot. Por. [ - ] à 0.2
Para definir el flujo natural en el acuífero, se deben imponer dos condiciones de borde de nivel
constante con diferencia de niveles.Para ello se realiza el siguiente procedimiento.
Paso 9: Hacer Click en el botón “Ok” en la ventana “Save property data before existing?”
Paso 10: Ir a “Assign” y seleccionar “Line”. Ahora con el botón izquierdo del mouse se
selecciona el punto de inicio de la condición de borde de nivel constante (Start Pt.) y se arrastra
hasta el punto donde termina la condición de borde (End Pt.). Una vez encontrado el punto de
término se hace click con el botón derecho del mouse y se despliega la ventana “Assign
Constant Head”.
Paso 12: Repetir nuevamente el paso 10 para insertar la otra condición de borde, en el otro
extremo del modelo.
Para una descarga continua de contaminante en un punto específico dentro del modelo se
deben seguir los siguientes pasos.
Paso 15: Ir al menú “File” y hacer click en “Setup” y seleccionar “Numeric Engines”.
Para ingresar las condiciones iniciales del modelo de trasporte, se deben seguir los siguientes
pasos.
Paso 18: Ir al menú y hacer Click en el botón “Properties” y seleccionar la condición “Initial
Concentration”.
Para ingresar los parámetros de dispersión se debe seguir los siguientes pasos:
Horiz/Long [ ] à 1
Horiz/Vert [ ] à 0
Diff Coeff [m²/day] à 0
Paso 28: Ir al menú y hacer Click en el botón “Boundaries” y seleccionar la condición “Constant
Concentration”.
Fuente de Descarga
Continua
Para insertar el pozo de observación que se ubica a 250 m del punto de la descarga se deben
seguir los siguientes pasos.
Paso 30: Ir al menú y hacer Click en el botón “Wells” y seleccionar la condición “Conc.
Observation Wells”.
Paso 31: Hacer Click en el botón “Ok” en la ventana “Save property data before existing?”
Para agregar el pozo se debe ir a la opción “Add Obs.” y luego hacer click en la posición en que
se ubica el pozo. En la pantalla se despliega la ventana “New Well”. En ella se agregan las
características del pozo, tales como su nombre, las coordenadas, etc.
X [m] à 625
Y [m] à 525
En la sección “Elevation Points” se debe poner:
Elevation [m] à 25
Aparece una nueva pantalla, en que se definen las condiciones de ejecución del programa.
Aparece la ventana “Ouput & Time Step Control” e ingresar los siguientes valores:
Simulation time à 2000 - Max # transport Steps à 2000 - Specify max stepsize à 100
Paso 40: Seleccionar las opciones “MODFLOW” y “MT3D/RT3D”. La primera opción permite
ejecutar el modelo de flujo y la segunda opción ejecuta el modelo de transporte.
Paso 41: Hacer Click en el botón “Translate Run”, para que el programa comience la ejecución.
En esta etapa se debe encontrar los resultados propuestos tales como los contornos
piezométricos, descensos y otros. En este caso los resultados de interés corresponden a la
parte de transporte del contaminante y se analizará el estado de avance de la pluma de
contaminación.
Paso 45: Para ver los distintos estados de ubicación de la pluma se debe hacer click en el
botón next (Time).
Paso 46: Hacer click en el botón “Conc. Vs Time” y activar en la parte derecha la opción
“W1.observacion” y “Conc001” en Species. Luego seleccionar Apply.
Paso 47: Para exportar los resultados en forma numérica, se debe hacer Click
en el botón “Export Data” en la parte superior del gráfico.
Paso 48: Para obtener estos valores en archivo de extensión dat o txt, es necesario hacer click
con el botón “Copy to File”. Para copiar información a planilla Excel se selecciona el botón
“Copy to Clipboard”.
• ¿Cuál es el efecto variar la tasa de dispersividad Horiz/Long en el paso 27? Pruebe con
valores de:
Horiz/Long à 0.01
Horiz/Long à 10
1900 2000 0
Avance de la pluma para un período de 2000 días.
Paso 53: En la sección inferior, seleccionar el menú “Species Params” e ingresar los siguientes
valores:
Paso 55: Seleccionar la opción “Re-initialize the transport Variant and save the new settings”.
Si se desea estudiar el efecto de variar la tasa de decaimiento, se den realizar los siguientes
pasos.
Paso 58: Ir al menú y hacer Click en el botón “Properties” y seleccionar la condición “Species
Parameters”.
Se despliega la ventana “Conc001 Parameters Database” y cambiar los valores deseados para
K_mobile. Notar que si K_mobile es cero, entonces el contaminante se comporta como
compuesto conservativo.