MODELACIÓN DEL TRANSPORTE REACTIVO EN ACUÍFEROS MEDIANTE MODFLOW 6 Y FLOPY: DESARROLLO DE UN CÓDIGO PROPIO EN PYTHON
Se presenta la modelación de flujo y transporte reactivo en acuíferos usando MODFLOW 6 (USGS) junto con el enfoque de componentes (de Simoni et al. 2005). El trabajo se centra en el desarrollo de un código propio en Python con la biblioteca FLOPY y en FORTRAN para incorporar las reacciones y la especiación (interfaz REMIX). Es validado usando un caso de referencia unidimensional de disolución de yeso y su solución analítica.
El modelo, en régimen transitorio, acopló en cada paso temporal MODFLOW 6 y la interfaz REMIX, mientras FLOPY gestionó la ejecución y el posprocesamiento en Python. Los perfiles de nivel piezométrico, velocidad y concentración obtenidos coincidieron satisfactoriamente con la solución analítica, superando obstáculos relacionados con cambios en la formulación matemática y en los formatos de datos de las nuevas versiones de las herramientas. Realizado en el marco del proyecto SECIHTI-CF-2023-G-904 “Soluciones basadas en la naturaleza a la contaminación del agua subterránea por especies de nitrógeno”.
de Simoni, M., Carrera, J., Sánchez-Vila, X. y Guadagnini, A. (2005). A procedure for the solution of multicomponent reactive transport problems. Water Resources Research, 41(11), 1–16. https://doi.org/10.1029/2005WR004056