UNIVERSIDAD DE LOS ANDES

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA CIVIL Y AMBIENTAL
BOGOTÁ D.C
2015

Actualización del modelo numérico del

acuífero Morroa utilizando Visual
Modflow Flex

PROYECTO DE GRADO EN INGENIERÍA CIVIL

Sergio Eduardo López Ramírez

Asesor: Carlos Molano

Co-asesor: Luis Alejandro Camacho
Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Tabla de contenido
Resumen ............................................................................................................................................. iv
1 Descripción General .................................................................................................................... 1
Localización geográfica .................................................................................................................. 1
Población ......................................................................................................................................... 2
Hidrografía ...................................................................................................................................... 3
Geomorfología ................................................................................................................................ 3
Obejtivos de modelación ................................................................................................................. 3
2 Modelo Hidrogeológico Conceptual ........................................................................................... 4
Geología .......................................................................................................................................... 4
Hidrología ....................................................................................................................................... 6
Antecedentes ................................................................................................................................... 6
3 Modelo Numérico ....................................................................................................................... 8
Antecedentes ................................................................................................................................... 8
Software utilizado ........................................................................................................................... 8
Área de modelación......................................................................................................................... 8
Datos de entrada .............................................................................................................................. 8
Construcción de superficies y capas ............................................................................................ 8
Grilla ......................................................................................................................................... 10
Pozos ......................................................................................................................................... 12
Propiedades hidráulicas ............................................................................................................. 14
Condiciones de frontera ............................................................................................................ 15
Simulación en Estado Estacionario ............................................................................................... 18
Simulación en Estado Transitorio ................................................................................................. 20
4 Conclusiones ............................................................................................................................. 26
5 Agradecimientos........................................................................................................................ 27
6 Referencias ................................................................................................................................ 28

Sergio Eduardo López Ramírez i

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Índice de Ilustraciones
Ilustración 1. Ubicación general del acuífero Morroa ......................................................................... 1
Ilustración 2 Área de modelación........................................................................................................ 9
Ilustración 3. Vista en planta y en perfil de las capas generadas Visual MODFLOW ..................... 10
Ilustración 4. Recorte parte norte de la grilla .................................................................................... 11
Ilustración 5. Recorte de la configuración vertical de la grilla fila 305 ............................................ 11
Ilustración 6. Recorte conductividad hidráulica fila 305 .................................................................. 15
Ilustración 7. Ríos asignados al modelo. De norte a sur: Arroyo Mancomoján, Arroyo Grande de
Corozal y Rio Camítico ..................................................................................................................... 16
Ilustración 8. Histograma de valores de recarga ............................................................................... 16
Ilustración 9. Zonas de recarga modelo numérico ............................................................................ 17
Ilustración 10. Valores de cabeza GBH Oriente ............................................................................... 18
Ilustración 11. Valores de cabeza GBH Sur ...................................................................................... 18
Ilustración 12. Calibración estado estacionario ................................................................................. 19
Ilustración 13. Tabla de agua estado estacionario ............................................................................. 20
Ilustración 14. Tabla de agua estado estacionario con pozos en funcionamiento ............................. 23
Ilustración 15. Vista 3D Tabla de agua costado oriental ................................................................... 23
Ilustración 16. Tabla de agua zona de mayor abatimiento ................................................................ 24
Ilustración 17. Vista tridimensional tabla de agua de la zona de mayor abatimiento ....................... 24

Sergio Eduardo López Ramírez ii

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Índice de Tablas
Tabla 1. Cobertura Total de Acueducto (Superintendencia de Servicios Públicos Domiciliarios, 2011)
............................................................................................................................................................. 2
Tabla 2. Espesor de las principales formaciones relacionadas al acuífero de Morroa (Abreu & Diaz,
2004) ................................................................................................................................................... 7
Tabla 3. Pozos de Bombeo (CARSUCRE) ....................................................................................... 13
Tabla 4. Resumen propiedades hidráulicas por capas ....................................................................... 14
Tabla 5. Valores de entrada condición de frontera ríos. Fuente: Juicio de expertos ......................... 15
Tabla 6. Proyectos de perforación próximos a entrar en operación .................................................. 21
Tabla 7. Caudales de bombeo ........................................................................................................... 22

Sergio Eduardo López Ramírez iii

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Resumen
El acuífero de Morroa es la principal fuente de agua para casi medio millón de personas en el
departamento de Sucre. Dada su gran importancia, es indispensable contar con herramientas que
ayuden a la toma de decisiones sobre su manejo y aprovechamiento. Después de la realización de
varios modelos hidrogeológicos conceptuales (Rodriguez, 1993) (Abreu & Diaz, 2004) (Pérez,
Obregón, & García, 2009), el desarrollo de un modelo numérico es un paso necesario con el fin de
establecer guías para la regulación y el mejor aprovechamiento de uno de los principales recursos de
los sucreños.
La necesidad de mejorar el servicio de abastecimiento de agua potable (principalmente para los
Municipios de Sincelejo y Carmen de Bolívar) trajo consigo la construcción de pozos de una gran
magnitud desde el año 2006; los últimos perforados tienen alrededor de 800 m de profundidad y
extraen caudales de cerca de 100 l/s. Debido al impacto que estos nuevos proyectos de perforación
pueden tener en el acuífero y en los pequeños pozos que sirven para el abastecimiento de comunidades
rurales, es importante modelar los impactos de estas iniciativas de desarrollo.
En éste trabajo se realiza la modelación numérica de 40 de los principales pozos de extracción
localizados en el acuífero de Morroa. El principal objetivo del actual modelo es determinar los efectos
producidos por los pozos de gran magnitud desarrollados a partir de 2006 y de 4 pozos más de
características similares próximos a ponerse en operación.
Los resultados del modelo muestran que los pozos de grandes profundidades y caudales tendrán un
impacto significativo sobre todo el acuífero siendo mayor principalmente en el campo de pozos de
Corozal.

Sergio Eduardo López Ramírez iv

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

1 Descripción General
El acuífero Morroa es la fuente de agua más asequible para los Municipios de Sincelejo, Corozal,
Ovejas, Morroa, Los Palmitos, San Juan de Betulia y Sampués en el departamento de Sucre. La
sobreexplotación del acuífero es tal vez el problema más grave que enfrente este recurso (Ingeominas,
2002). (Herrera, Vargas, & Taupin, 2006). Debido a su origen geológico, este acuífero posee una gran
heterogeneidad e anisotropía en sus propiedades hidráulicas .Los recursos dinámicos del acuífero
Morroa se calculan en 160 Mm3/año (5m3/s) basados en unas reservas del orden de 1’661.026 Mm3
(Ingeominas, 2002)

Localización geográfica
El acuífero Morroa se ubica en el departamento de Sucre, zona noroccidental del territorio
Colombiano. Se extiende en dirección suroeste-noreste en una franja de aproximadamente unos 5 km
de ancho que va desde Sampués en límites con el departamento de Córdoba, hasta Ovejas en cercanías
del departamento de Bolívar.

Ilustración 1. Ubicación general del acuífero Morroa

Sergio Eduardo López Ramírez 1

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Población
El departamento de Sucre tiene una población aproximada de 843.000 habitantes; se estima que el
total de la población de los 7 municipios que se abastecen actualmente del acuífero Morroa es de
438.000 habitantes (DANE, 2014). Lo anterior significa que aproximadamente el 52% de la población
de éste departamento tiene como fuente de agua potable este acuífero.
En cuanto a la cobertura de acueducto, la Tabla 1 muestra el porcentaje de cubrimiento para cada uno
de los 7 municipios que en la actualidad se abastecen del acuífero Morroa. Nótese que 5 de estos
municipios tienen una cobertura inferior al 80% para el año 2011. Este escenario muestra la necesidad
de realizar nuevas desarrollos en el acuífero con el fin de aumentar significativamente las tasas de
cobertura en estos municipios.
Tabla 1. Cobertura Total de Acueducto (Superintendencia de Servicios Públicos Domiciliarios, 2011)

Municipio Cobertura Total de Acueducto (%)

Corozal 90.10
Los Palmitos 76.50
Morroa 74.80
Ovejas 79.80
Sampués 74.30
San Juan de Betulia 69.60
Sincelejo 79.90

En este mismo sentido, la Gráfica 1 muestra la necesidad de mejorar sustancialmente las condiciones
de suministro de agua; nótese que el indicador de Necesidades Básicas Insatisfechas de los municipios
del acuífero se encuentra cerca de un 10 % por encima del promedio nacional.

Necesidades básicas insatisfechas

80

70 Necesidades
Básicas
60 Insatisfechas
- NBI
50
Promedio 7
(%)

40 municipios
30

20 Promedio
Colombia
10

Municipio
Gráfica 1. Necesidades básicas insatisfechas municipios acuífero Morroa (DANE, 2010)

Sergio Eduardo López Ramírez 2

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Hidrografía
La zona de estudio se ubica en la margen oriental de los Montes de María donde los arroyos son
temporales y se forman de manera impetuosa después de aguaceros fuertes. En épocas de sequía el
agua escasea por lo que se hace necesaria la construcción de sistemas de almacenamiento temporal
del agua conocidos como jagüeyes. El principal río que se encuentra en la zona es el arroyo grande
de corozal que recoge las aguas residuales de varios municipios, apartándose así de su carácter
intermitente.

Geomorfología
Sucre se encuentra ubicado dentro de la llanura caribe que cubre la parte nor-occidental del territorio
colombiano. Dentro de este departamento se pueden distinguir 4 subregiones principales: la llanura
costera, la zona montañosa o de colinas, las sabanas y la depresión inundable del bajo Cauca y San
Jorge. El acuífero Morroa se encuentra ubicado en la zona de sabanas. Ésta tiene un área de unos
2097 km2, está formada por una serie lomas pequeñas con alturas entre 10 y 20 metros por encima
del nivel circundante. Corresponde a la parte oriental de los Montes de María donde a manera de
declive van hacia las áreas bajas de los ríos San Jorge, Cauca y Magdalena.

Obejtivos de modelación
El principal objetivo es el montaje de un modelo numérico del acuífero Morroa haciendo uso de la
herramienta computacional MODFLOW. El desarrollo del modelo se hará mediante el uso de la
herramienta comercial Visual Modflow Flex, la cual hace uso del sistema de modelación de
diferencias finitas MODFLOW-2005 desarrollado por el servicio geológico de la Unión Americana
(USGS por su siglas en inglés). A partir de la construcción del modelo se espera:
 Obtención de un modelo numérico a escala regional del acuífero de la formación Morroa
 Obtención de un modelo que represente de manera satisfactoria las condiciones del acuífero
de tal manera que pueda serle útil a la autoridad ambiental
 Evaluar la afectación generada por los pozos de gran magnitud construidos desde 2006.
 Evaluar la afectación generada por los próximos pozos a entrar en operación en la zona.

Sergio Eduardo López Ramírez 3

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

2 Modelo Hidrogeológico Conceptual

Geología
El siguiente mapa muestra el resumen de la información geológica de la zona de estudio

Mapa 1 Mapa geológico del acuífero Morroa (CARSUCRE, 2014)

Las formaciones que representan el mayor interés para el proyecto de modelación numérica son: la
Formación Morroa, la Formación Sincelejo (Inferior y Superior) y la Formación Betulia. La primera
en mención, constituye la unidad principal del acuífero; la segunda infrayace a la formación Morroa

Sergio Eduardo López Ramírez 4

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

y la tercera suprayace a la primera en mención. A continuación se presenta de manera sucinta la

descripción geológica de estas tres formaciones de acuerdo al trabajo realizado por el INGEOMINAS
en 1992 (Arias & Fuquen, 1992)
Formación Sincelejo Inferior
Secuencia de areniscas friables de grano variable, lentejones erráticos de conglomerados con
predominancia de cuarzo blanco, liditas y rocas ígneas, lentes de arcilla y limolitas de acuerdo a la
descripción dada por Kassem en 1967 (Cáceres, Cuacalón, & Kassem, 1967). Su espesor es variable
considerando un promedio de 200 m. Esta formación aflora en los alrededores del municipio de
Sincelejo descansando discordantemente sobre arcillolitas rojas ferruginosas de la formación El
Carmen. A ésta formación se le asigna una edad del Plioceno.
Formación Sincelejo Superior
Es una secuencia compuesta por arcillolitas con intercalaciones de arenisca y lentes de conglomerado.
En las canteras situadas al norte de Ovejas está en contacto fallado con la formación El Carmen,
mientras que en los alrededores de la ciudad de Sincelejo descansa, al parecer en conformidad, sobre
el Sincelejo Inferior. Su espesor es de aproximadamente 200 m y su edad es del Plioceno.
Formación Morroa
Es una secuencia de areniscas, conglomerados y arcillas con nódulos calcáreos, la cual presenta
rápidos cambios de facies tanto laterales como verticales. De acuerdo a lo observado en el mapa
geológico, esta formación se extiende desde el Norte de Sampués hasta el Oriente de Ovejas teniendo
una dirección mayoritariamente de N25°E. El buzamiento de la formación es variable siendo menor
en el tope que en la base, el cual está entre 5 y 25° aproximadamente. La formación Morroa descansa
con un ángulo de cerca de 20° de inclinación sobre el Sincelejo Superior. El contacto entre estas dos
formaciones es discordante. Por su parte, el contacto con la Formación Betulia es aparentemente
concordante. El espesor de la formación es de aproximadamente 1170 m y se le asigna una edad
Plioceno- Pleistoceno.
La formación Morroa es de origen fluvial caracterizada por la presencia de conglomerados y areniscas
de grano grueso en capas lenticulares con estratificación cruzada y rápidos cambios faciales tanto
verticales como horizontales. Lo anterior, impide la presencia de niveles guía continuo que faciliten
la correlación de los diferentes estratos con los lugares donde están expuestos.
Formación Betulia
Formación de carácter lacustre, está constituida por una serie de arcillolitas abigarradas con niveles
arenosos. Aflora en la parte oriental del área de estudio y se encuentra suprayaciendo la Formación
Morroa en aparente en concordancia. Puede alcanzar hasta1000 m de espesor en su parte más
desarrollada (Abreu & Diaz, 2004).

Sergio Eduardo López Ramírez 5

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Hidrología
Precipitación
La precipitación en la zona está entre los 1100 y los 1500 mm/año de acuerdo al mapa realizado por
INGEOMINAS para 1992 (Rodriguez, 1993). Se distingue un periodo seco de diciembre a marzo y
un periodo húmedo con una distribución uniforme de abril a noviembre. Tres cuartas partes de la
precipitación total se tiene para los meses de mayo a octubre.
Recarga
Los valores de recarga para la zona han sido estimados entre 34 y 117 mm/año (CARSUCRE, 2005),
registros que son similares a los encontrados para 1992 en el estudio hidrogeológico de
INGEOMINAS (Rodriguez, 1993) los cuales oscilan entre 59.3 y 179 mm/año. Por una parte los
procesos de meteorización y erosión de las rocas superficiales aumentan las arcillas en este nivel; lo
anterior lleva a una continua disminución de la recarga del sistema. De otra parte, mediante el análisis
de la situación local, los bajos valores de recarga pueden estar asociados en parte a la compactación
del suelo producida por la prácticas agrícolas y ganaderas inadecuadas (CARSUCRE, 2005).

Antecedentes
Los últimos modelos hidrogeológicos conceptuales desarrollados para la zona corresponden a los
trabajos de Abreu y Díaz en 2004 (Abreu & Diaz, 2004)y al trabajo publicado por Pérez, Obregón y
García en 2009 (Pérez, Obregón, & García, 2009). Aunque el trabajo más reciente propone un
concepto mucho más riguroso en cuanto a la distribución de las propiedades hidráulicas de acuerdo a
la caracterización geoestadística del terreno, el trabajo de Abreu y Díaz resulta en una aproximación
más práctica para el alcance de este proyecto de grado.
En el primer trabajo señalado, se hizo una recopilación de información de líneas sísmicas, registros
de pozos petroleros, imágenes satelitales y trabajo de campo para proponer un nuevo modelo
hidrogeológico de la formación Morroa. Los aspectos más importantes de éste trabajo resaltan la
diferenciación de la formación Morroa en tres unidades acuíferas principales: Morroa, Sincé y
Albania llamadas así por la determinación de los autores.
La unidad Albania es descrita como un paquete grueso de areniscas de grano grueso a
conglomeráticas con intercalaciones de delgados lentes de arcilla. Por su parte, el acuífero de Sincé
se distingue como un paquete grueso de conglomerados y areniscas conglomeráticas con intercalación
de algunos lentes arenosos.
En cuanto a la geometría de los acuíferos, Abreu y Díaz coinciden con las afirmaciones anteriores en
el hecho de que el acuífero de Morroa trata de un gran monoclinal con buzamiento y acuñamineto
hacia el este. En este punto es importante destacar que la interpretación de líneas sísmicas realizada
en este trabajo ayudó a la consolidación del modelo geológico hacia el oriente del afloramiento de la
formación. Los autores sugieren que el límite oriental del acuífero puede coincidir con la falla
Romeral-Bucaramanga.
De otro lado, Abreu y Díaz destacan que los acuíferos de Sincé y Albania se corresponden a un
sinclinal bastante amplio cuya forma obedece al peso por la columna de sedimentos. Debido a la falta
de información sobre la rata de perforación de los pozos petroleros, solamente se tiene un referente
del orden de magnitud del espesor de cada una de las formaciones analizadas en este trabajo. Dicha
información es presentada en la Tabla 2.

Sergio Eduardo López Ramírez 6

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Finalmente, en cuanto a las propiedades hidráulicas, Abreu y Díaz realizaron la interpretación de18
pruebas de bombeo mediante métodos manuales y haciendo uso del software Aquifer Test. .Lo
anterior les permitió observar que los valores de conductividad hidráulica obtenidos correspondían a
un acuífero confinado. Solamente una de las pruebas de bombeo tenía la suficiente información como
para determinar el valor del coeficiente de almacenamiento; determinado este valor se corroboró que
se trata de acuífero confinado.
Tabla 2. Espesor de las principales formaciones relacionadas al acuífero de Morroa (Abreu & Diaz, 2004)

Formación Espesor
(m)
Sincelejo 1x102
Morroa 1x103
Albania 3x102
Sincé 2x102

Sergio Eduardo López Ramírez 7

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

3 Modelo Numérico
Antecedentes
El primer trabajo de modelación numérica de este acuífero fue realizado por Santisteban en 1989
(Santisteban, 1989) quién diseño un programa para la solución de las ecuaciones de flujo en dos
dimensiones a través de un esquema de malla centrada.
Después, en el año 2000, Buitrago y Donado (Buitrago & Donado, 2000) construyeron un modelo
para el acuífero Morroa haciendo uso del programa Visual MODFLOW en la versión 2.7.1. Este
trabajó consto tanto del desarrollo del modelo hidrogeológico conceptual como del modelo numérico.
El modelo desarrollado fue de carácter regional cubriendo un área de 1471 km2.
Las características principales de este modelo señalan la zonificación del modelo hidrogeológico a
partir del estudio del INGEOMINAS realizado en 1993 (Rodriguez, 1993) De esta manera se
adoptaron valores de las propiedades hidráulicas y de la recarga por zonas. En este estudio no se
incluyó la existencia de capas impermeables en medio de alguna de las zonas del proyecto. En cuanto
a las principales características de la grilla, se utilizaron celdas de 90x80 m y para el modelamiento
tanto del estado estacionario como el transitorio se tomó como año base 1993.

Software utilizado
Para el desarrollo de este modelo numérico se utilizó el software Visual MODFLOW Flex versión
2014.2 desarrollado por Schlumberger Water Services. El anterior software permite el ingreso de un
modelo conceptual al programa a partir del cual se generan las celdas con sus respectivas propiedades
hidráulicas. Lo anterior significa que a partir de la creación de la grilla y de las zonas con sus
respectivas propiedades dentro de un modelo conceptual, el programa realiza la asignación de las
propiedades a cada una de las celdas del modelo.

Área de modelación
El modelo desarrollado es de carácter regional debido a que cubre los dos campos de pozos principales
del acuífero: el de Corozal y el de Ovejas. A su vez, cobija los nuevos pozos desarrollados para
suministrar agua a Carmen de Bolívar, los cuales, están localizados en Los Palmitos. De igual manera,
dentro del área de modelación del proyecto se encuentran los pozos de gran profundidad desarrollados
en Corozal a partir del 2006. El límite oriental del modelo está a 2.5 km del pozo más al este de la
modelación, mientras que el occidental coincide con el límite en al afloramiento de la formación
Morroa. El área tiene una extensión de 870 km2 teniendo un ancho de 15 km y un largo de 58 km. El
área de modelación se puede observar en la Ilustración 2.

Datos de entrada

Construcción de superficies y capas

Para la construcción de la superficies del modelo numérico se procedió a dibujar perfiles transversales
con puntos de interpolación de acuerdo al Mapa Hidrogeológico y los perfiles desarrollados en 1992
por el INGEOMINAS (Rodriguez, 1993). A su vez, para tener en cuenta los últimos trabajos
desarrollado en la descripción del modelo hidrogeológico (Abreu & Diaz, 2004), se disminuyó el
buzamiento de las capas a partir del lugar en que el espesor las zonas más permeables alcanzaba los
1000 m de espesor. Los perfiles se realizaron para puntos a cada 1500m.

Sergio Eduardo López Ramírez 8

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 2 Área de modelación

Una vez se definidos los perfiles se procedió a dibujar las superficies que conforman los límites de
cada una de las capas del modelo. Se utilizó la herramienta kirging incluida en el paquete de Visual
Modflow Flex para el dibujo de éstas superficies. En algunas capas fue necesaria la creación de puntos
adicionales con el fin de emular con mayor precisión las formas dadas por el mapa hidrogeológico de
la zona. El límite que coincide con la superficie del terreno se construyó a partir del modelo de
elevación digital ASTER desarrollado por el Ministerio de Economía Industria y Comercio de Japón

Sergio Eduardo López Ramírez 9

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

(METI por sus siglas en inglés) y la NASA. En la Ilustración 3 se puede observar una vista en planta
y un perfil de las capas generadas mediante el método anteriormente mencionado.

Ilustración 3. Vista en planta y en perfil de las capas generadas Visual MODFLOW

El modelo se construyó a partir del dibujo de 4 capas permeables dentro de la formación Morroa
separadas por 3capas arcillosas. A su vez la Formación Morroa está limitada por capas impermeable
en su base y en su tope: la formación Sincelejo Superior y la formación Betulia respectivamente. Es
decir, en total se construyeron 9 capas: 5 impermeables representando las formaciones que infrayacen
y suprayacen la formación Morroa y las intercalaciones de arcillolitas en medio de la formación; y 4
capas permeables representando los principales niveles permeables descritos en el modelo
hidrogeológico conceptual. Debido a que el programa Visual MODFLOW Flex permite la inclusión
de valores a partir de las zonas, la construcción de superficies y el cruzamiento de las mismas
conformaron capas que se pinchan y por ende los valores de conductividad hidráulica son
independientes de la forma en cómo se construya la grilla. Es decir, no existe una completa
independencia entre capas impermeables tal como lo sugiere Rodríguez para 1993.
Los cortes realizados por el INGEOMINAS sirvieron para determinar el espesor de las capas
impermeables en la zona sur del acuífero. Para definir el espesor de las mismas en la parte norte, se
utilizó la información de registros eléctricos y diseños de los Pozos A y B ubicados en Los Palmitos.
Se asumió que los afloramientos arcillosos de la formación Morroa en la zona de Ovejas corresponden
a la continuación de las capas que se encuentran en Corozal. Debido a la propia construcción de las
superficies, las zonas de baja permeabilidad no aislaron completamente las capas permeables.
Grilla

Debido a que se trata de un modelo de escala regional se utilizaron grillas de un tamaño 250x250m
en las zonas exteriores del área de modelación, donde se encuentra la menor cantidad de pozos. Con
el fin de mejorar la precisión del modelo en las zonas de explotación más importantes, se hizo un
refinamiento de la grilla para cubrir la mayor parte de los pozos del campo de Corozal; en dicho lugar,
la grilla es de 62.5x62.5m. La grilla se giró 335° en el sentido opuesto a las manecillas del reloj con
el fin de que coincidiera con el rumbo de 25° NE predominante en el acuífero. En la Ilustración 4 se

Sergio Eduardo López Ramírez 10

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

muestra un recorte de la parte norte de la grilla creada donde se puede observar la rotación, las celdas
inactivas (color azul) y el refinamiento de la grilla.

Ilustración 4. Recorte parte norte de la grilla

En la vertical, se utilizaron 10 capas del tipo “deformed”; lo anterior significa que las capas se
construyen deformadas de acuerdo al espesor asignado en el modelo conceptual. El espesor mínimo
de las capas es de 2 m, el cual, se establece con el fin de mejorar la estabilidad numéroca del modelo.
Un ejemplo de la configuración vertical del acuífero se puede observar en la Ilustración 5. En total
para 10 capas, 583 filas y 169 columnas se obtuvieron 985.279 celdas para la modelación regional
del acuífero.

Ilustración 5. Recorte de la configuración vertical de la grilla fila 305

Sergio Eduardo López Ramírez 11

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Pozos

Pozos de bombeo
Con la información suministrada por CARSUCRE se ingresaron al modelo 40 pozos de bombeo. La
Tabla 3 muestra el resumen de la ubicación y la profundidad de los mismos. De los pozos
mencionados en la Tabla 3, se tienen valores de medición de caudal con precisión a partir de los años
2001 y 2002, momento en que entró en funcionamiento el Proyecto de Protección Integral de Aguas
Subterráneas (PPIAS). A partir de la restricción anterior, los valores de caudal para fechas anteriores
a éstos registros se asumieron como los valores de caudal del primer registro de medición.
Niveles piezométricos
Se obtuvo información histórica de niveles para 38 de los 40 pozos de bombeo mencionados en la
Tabla 3 (no se obtuvo información para los pozos 44-II-D-PP-13 y Pozo B). Al igual que lo sucedido
con los registros de caudales, la información de niveles se presenta en detalle a partir de los años 2001
y 2002.

Sergio Eduardo López Ramírez 12

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Tabla 3. Pozos de Bombeo (CARSUCRE)

Código X Y Z Profundidad
44-II-D-PP-07 876678 1546720 214 255
44-II-D-PP-08 876524 1547175 219 208
44-II-D-PP-13 876134 1547310 210 194
44-IV-C-PP-01 864789 1523252 178 250
44-IV-C-PP-06(31) 864222 1522713 171 359
44-IV-C-PP-07(33) 864900 1521400 140 333
44-IV-C-PP-08(34) 864330 1520560 159 385
44-IV-D-PP-16(16) 866410 1523510 165 146
44-IV-D-PP-24(24) 867670 1522400 144 233
44-IV-D-PP-25(25) 866800 1524100 164 207
44-IV-D-PP-30 866598 1525810 184 391
44-IV-D-PP-31(32) 867703 1523516 157 359
44-IV-D-PP-34 866640 1521160 136 205
44-IV-D-PP-35(36) 865746 1520291 160 367
44-IV-D-PP-36 868120 1524880 177 235
44-IV-D-PP-38(02) 866848 1523429 151 404
44-IV-D-PP-39 (Concesión) 865579 1525441 170 312
44-IV-D-PP-40(42) 868104 1523010 148 359
44-IV-D-PP-41(44) 867850 1521330 150 398
44-IV-D-PP-42 (03) 869080 1523263 178 401
44-IV-D-PP-43 (43B) 870073 1523217 197 468
44-IV-D-PP-44 (45) 865664 1522222 142 505
44-IV-D-PP-45 (Nuevo Rincón) 868263 1526963 120 279
44-IV-D-PP-46 871591 1525158 160 800
44-IV-D-PP-47 872666 1522898 134 867
44-IV-D-PP-48 869956 1528106 160 239
52-II-A-PP-01 857060 1512470 140 42
52-II-A-PP-11(37) 864787 1518406 164 367
52-II-A-PP-14 (38) 861400 1515320 172 385
52-II-A-PP-15(41) 863880 1519150 188 378
52-II-A-PP-16(40) 862650 1515420 168 280
52-II-A-PP-22 862900 1514500 165 348
52-II-B-PP-01(35) 865622 1519231 155 354
52-II-B-PP-02(49) 866414 1515964 148 730
52-II-C-PP-18 (Concesión) 857140 1506230 153 145
52-II-C-PP-21 855210 1504950 141 145
52-II-C-PP-22 855170 1509500 148 104
52-II-C-PP-24 (Carretero) 855504 1506930 152 41
Pozo A 876337 1539054 210 778
Pozo B 877653 1538184 229 472

Sergio Eduardo López Ramírez 13

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Propiedades hidráulicas

La obtención de ésta información requiere tanto de una adecuada interpretación de las pruebas de
bombeo, el diseño del pozo y la estratigrafía en el lugar de perforación. La información utilizada para
la asignación de propiedades a cada una de las 9 zonas del modelo conceptual se resume en la Tabla
4.
Tabla 4. Resumen propiedades hidráulicas por capas

Zona T (m²/d) Espesor Ke, xy (m/d) Kz (m/d) Ss Sy

1 - 1000 5.00E-04 5.00E-05 1.00E-07 0.08
2 20-215 300 1.00E-01 1.00E-02 6.70E-04 0.15
3 - 30 5.00E-04 5.00E-05 1.00E-10 0.08
4 8-235 333 1.50E-01 1.50E-02 3.30E-04 0.13
5 - 40 5.00E-04 5.00E-05 1.00E-10 0.08
6 85 180 1.00E-01 1.00E-02 1.00E-04 0.11
7 - 35 5.00E-04 5.00E-05 1.00E-10 0.08
8 100 200 1.00E-01 1.00E-02 1.00E-04 0.1
9 - 60 5.00E-04 5.00E-05 1.00E-10 0.08

Para la definición de las propiedades hidráulicas se contó con el apoyo del geólogo Alcides Hugget.
Se utilizaron los datos de las pruebas de bombeo de los pozos 44-IV-D-PP-14, 44-IV-D-PP-17, 44-
IV-D-PP-18, 44-IV-D-PP-24, 44-IV-D-PP-02, 44-IV-D-PP-05, 44-IV-D-PP-16, 44-IV-D-PP-22 y
44-IV-C-PP-26 para la asignación de las propiedades hidráulicas a cada una de las capas del modelo.
Para las capas productoras los valores de conductividad hidráulica se asignaron mediante el análisis
de las pruebas de bombeo y la identificación, a partir de los registros eléctricos, del espesor
equivalente de cada capa para cada uno de los pozos. Para las capas arcillosas del modelo se asumió
una conductividad de 5E-4 m/d. Para simular en el modelo las condiciones de anisotropía se asumió
que la conductividad hidráulica vertical equivale al 10% de la horizontal.
En cuanto a la capacidad de almacenamiento de las capas del modelo, para las aquellas productoras
se asumió como el promedio aritmético hallado para los diferentes pozos que captan el agua de cada
capa. Para las aquellas arcillosas se tomó el valor de 1E-10.La Ilustración 6 muestra un corte
representativo de la conductividad hidráulica asignada a las celdas del modelo.

Sergio Eduardo López Ramírez 14

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 6. Recorte conductividad hidráulica fila 305

Condiciones de frontera

Se asignaron condiciones de frontera de tres tipos para la definición del modelo: se ingresaron tres
ríos, se incluyó la recarga y fronteras del tipo III o generales (GBH por sus siglas en inglés).
Ríos
Se utilizaron tres condiciones de frontera de tipo rio en el modelo: el Rio Camítico, el Arroyo Grande
de Corozal y Arroyo Mancomoján. Debido a la falta de información suficiente de estos ríos se
asumieron los valores de entrada para el programa de acuerdo a lo mostrado en la siguiente tabla.
Tabla 5. Valores de entrada condición de frontera ríos. Fuente: Juicio de expertos

Espesor del lecho del Rio (m) 3

K (m/d) 2
Ancho (m) 4

Se obtuvo información de los registros de niveles para el Arroyo Grande de Corozal disponible en el
IDEAM. De acuerdo a los niveles registrados para este rio se asumió una profundidad media de 10
cm para todas las condiciones de este tipo. La Ilustración 7 muestra la ubicación de estos ríos dentro
del área de modelación.

Sergio Eduardo López Ramírez 15

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 7. Ríos asignados al modelo. De norte a sur: Arroyo Mancomoján, Arroyo Grande de Corozal y Rio
Camítico

Recarga
Se utilizaron las zonas definidas por Rodríguez (Rodriguez, 1993) para la asignación inicial de la
recarga en el modelo, tal como se muestra en la Ilustración 9. De acuerdo a dicho esquema, la Zona
1 corresponde a los lugares donde no hay recarga. Inicialmente se supuso una recarga de 40 mm/año
para todas las zonas (zona 2 a zona 6). Para la mejora de la calibración y de la modelación numérica,
se buscó aumentar de manera gradual los valores de recarga a medida que se pasaba de las zonas
impermeable a las permeables. El valor de recarga máximo utilizado, de acuerdo a la información
recopilada, fue de 115 mm/año; el mínimo utilizado dentro de las zonas de recarga fue de 0 mm/año.
El promedio de los valores de recarga utilizado fue de 24.6 mm/año y la mediana es de 18 mm/año.
El histograma de recarga se muestra en la Ilustración 8

Histograma (Recarga)
0.35

0.3
Frecuencia relativa

0.25

0.2

0.15

0.1

0.05

0
0 20 40 60 80 100 120 140
Recarga

Ilustración 8. Histograma de valores de recarga

Sergio Eduardo López Ramírez 16

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Tal como se puede observar en la anterior tabla, la mayor parte de los valores de recarga se encuentran
entre 0 y 36 mm/año. Después de esto vienen dos grupos de recarga importantes: el primero (48-60
mm/año) es un aumento de recarga que se hizo al sur de la zona VI; el segundo (96-108 mm/año)
corresponde a los valores de recarga de la zona V.

Ilustración 9. Zonas de recarga modelo numérico

Condiciones de frontera de tipo general

Se utilizaron dos fronteras de este tipo para emular la influencia del agua que se encuentra almacenada
a las afueras del área representada en el modelo numérico. La primera se ubicó en la zona oriental del
modelo con valor de cabeza variable entre los 80 m y los 160 m; siendo éste valor más alto en la zona
norte del acuífero y disminuyendo de manera gradual hacia el sur. Con esta frontera se busca
representar el almacenamiento que existe en la zona oriental del modelo sin tener la necesidad de
extenderse hasta esa área.
La segunda frontera de este tipo se ubicó en el límite sur del modelo numérico con el fin de emular
el volumen de agua que pueda existir en el acuífero hacia las inmediaciones del departamento de
Córdoba. El valor de cabeza constante para esta frontera fue de 128m, el cual se obtuvo con el fin de
conseguir una mejor calibración del modelo en estado estacionario. Para ambas condiciones de
frontera se tomó una conductividad hidráulica 0.0503 m/d, la cual corresponde a la conductividad
promedio de las capas ingresadas en el sistema. La frontera oriental se supuso a 5000m del modelo
mientras que la sur a 1000 m. En la Ilustración 10 y en la Ilustración 11 se muestran de manera
esquemática los valores de cabeza obtenidos a lo largo de las filas y columnas del modelo
respectivamente

Sergio Eduardo López Ramírez 17

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

GBH Oriente
180
160
140
120
Cabeza (m)

100
80
60
40
20
0
0 100 200 300 400 500 600
Fila

Ilustración 10. Valores de cabeza GBH Oriente

GBH Sur
128.5

128

127.5
Cabeza (m)

127

126.5

126

125.5

125
0 20 40 60 80 100 120 140 160 180
Columna

Ilustración 11. Valores de cabeza GBH Sur

Simulación en Estado Estacionario

Para realizar la calibración en estado estacionario se supusieron los niveles estáticos de 38 de los 40
pozos para el tiempo de inicio de la explotación del acuífero (aproximadamente para el año 1955).
Los primeros registros del acuífero indican que el nivel del agua se encontraba aproximadamente a 5
m de la superficie (Insfopal, 1981). Dado esto, se construyeron los niveles observados como la resta
de 5 m al nivel de superficie para cada pozo. A través de un proceso de prueba y error, se ajustaron
los valores de las condiciones de frontera con el fin de lograr lo mejor calibración de la curva de
niveles cales calculados contra niveles observados.

Sergio Eduardo López Ramírez 18

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Configuración de las corridas

Para el estado estacionario -obtención de la solución de la ecuación de flujo en el medio poroso
haciendo igual a cero el término de almacenamiento- se tuvo en cuenta los siguientes aspectos. Para
la solución de las ecuaciones flujo se utilizó el método GMG (Geometric Multigrid Solver) debido a
su rapidez en la solución del sistema. El número máximo de iteraciones externas fue de 500 mientras
que el de iteraciones internas fue de 8; los valores tanto del criterio de convergencia de cambio de
cabezas así como del criterio de convergencia residual fueron de 0.01. Usualmente la solución del
sistema en estado estacionario tardaba alrededor de 2 minutos teniendo como valor inicial de cabezas
la superficie del terreno (nótese que este valor no tiene ninguna influencia en el resultado final del
cálculo más allá del tiempo de simulación).
Se utilizó el paquete BCF6 para la asignación de transmisividad de las capas, siendo las capas de la
1 a la 10 del tipo 3. Lo anterior significa que la transmisividad podría ser calculada como si se tratara
de un acuífero confinado o uno inconfinado, dependiendo del valor de la cabeza en la celda. Para
todas las capas la conductancia entre celdas adyacentes fue calculada a partir de uso de la media
armónica. No se hizo uso del remojo para la solución del modelo ya que no fue posible lograr la
convergencia una vez se habilitaba esta función. A continuación, se muestran los resultados obtenidos
en la calibración en estado estacionario. El error medio de la calibración resultó ser de -0.99 m y el
error estándar de 1.18m. El máximo error en calibración se dio para el pozo 52-II-A-PP-01 con un
valor de + 11m.

Ilustración 12. Calibración estado estacionario

Sergio Eduardo López Ramírez 19

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 13. Tabla de agua estado estacionario

En la Ilustración 13 se puede observar la tabla de agua obtenida como resultado de la simulación del
estado estacionario. De acuerdo al valor de las cabezas se puede ver como la dirección del agua toma
una orientación principal en sentido noroeste –sureste. Los niveles más altos para el año en que se
inicia la explotación del acuífero se encuentran al norte del mismo; los niveles de la tabla de agua de
esta manera emulan las condiciones naturales del terreno.

Simulación en Estado Transitorio

Se recopiló información de los nuevos proyectos de perforación a desarrollarse en la zona; las
principales características de éstos pozos se resumen en la Tabla 6.

Sergio Eduardo López Ramírez 20

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Tabla 6. Proyectos de perforación próximos a entrar en operación

Código del Pozo X Y Z Profundidad Caudal (m3/d)

PP-48-ADESA 873376 1528102 195 850 6480
PP-51-ADESA 875421 1529786 184 825 6480
PP-50-ADESA 874764 1525992 195 800 6480
PP-55-ADESA 867379 1595039 171 800 6480

A pesar de los numerosos intentos en tratar de realizar la simulación en estado transitorio, ésta no fue
posible de hacer por Problemas encontrados principalmente con el uso del software Visual Modflow
Flex, los cuales impidieron la ejecución apropiada del modelo en este estado. Varios de los errores
encontrados sobre el funcionamiento del software fueron resueltos gracias al soporte técnico de la
empresa comercializadora del mismo. Sin embargo, al momento de la terminación del tiempo de
presentación del proyecto de grado, un error relacionado con la forma en que el modelo corría en
estado transitorio no había sido corregido. Lo anterior, supone que el software utilizado, bajo las
condiciones actuales, no permite la obtención de una herramienta adecuada para la toma de decisiones
por parte de la autoridad ambiental. Es importante aclarar que la corrección de dicho error se
encuentra fuera del alcance del presente proyecto de grado.

Como alternativa a la corrida del modelo numérico en estado transitorio se realizó la corrida del
mismo en estado estacionario haciendo ingreso de los caudales de explotación esperados para el año
2035. Para dicho año se asumió la entrada en operación de los pozos anteriormente mencionados en
la Tabla 6, junto con los últimos caudales de explotación reportados para los pozos que ya se
encuentran en funcionamiento en la región. El caudal total utilizado para dicha modelación es de 994
L/s tal como se muestra en la Tabla 7. Nótese que dicho caudal no pudo ser calibrado debido a la
imposibilidad de hacer la corrida en estado transitorio. En la Ilustración 14 y en la Ilustración 15 se
presentan los resultados obtenidos bajo esta simulación; tales valores no pueden ser vistos de manera
absoluta. Lo más importante a tener en cuenta sobre la modelación realizada son los órdenes de
magnitud obtenidos dentro de la simulación.

A partir de la comparación de las imágenes obtenidas se puede observar que el descenso promedio
esperado puede estar cerca de 100 m. El 93 % de la tabla de agua para la simulación con el efecto de
los pozos se encuentra entre -344 y 101 m, mientras que para el estado estacionario el 96% de las
cabezas se encontraban entre 115 y 207m. De manera general se puede observar que se genera un
gran abatimiento en la zona oriental del campo de pozos de Corozal; justamente en el lugar donde se
ubicarán los principales pozos de explotación en un futuro inmediato.

La zona que más se afectará, en la que más se disminuirá la altura de su tabla de agua, se encuentra
cercana al campo de pozos de Corozal. Esta área se ve principalmente influida por los siguientes
pozos: 44-IV-PP-47, 44-IV-D-PP-46 y PP-50-ADESA. En esta región la tabla de agua señalada
alcanzaría incluso una altura de 600m por debajo del nivel medio del mar. Los abatimientos de esta
zona se pueden evidenciar tanto en la Ilustración 16 como en la Ilustración 17. Nótese también que
esta área tiene un aumento significativo de cabezas en el costado sur del pozo 44-IV-PP-47, donde se
encuentra justamente el Arroyo Grande de Corozal. En la Ilustración 15 se puede ver como al sur del
Arroyo Grande de Corozal se alcanza otro gran abatimiento de la tabla de agua atribuido
principalmente al pozo 52-II-B-PP-02(49).

Sergio Eduardo López Ramírez 21

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Tabla 7. Caudales de bombeo

Caudal
Pozo
(m3/d)
44-II-D-PP-07 125
44-II-D-PP-08 85
44-II-D-PP-13 20
44-IV-C-PP-01 137
44-IV-C-PP-06(31) 89
44-IV-C-PP-07(33) 387
44-IV-C-PP-08(34) 70
44-IV-D-PP-16(16) 299
44-IV-D-PP-24(24) 315
44-IV-D-PP-25(25) 158
44-IV-D-PP-30 72
44-IV-D-PP-31(32) 97
44-IV-D-PP-34 138
44-IV-D-PP-35(36) 87
44-IV-D-PP-36 31
44-IV-D-PP-38(02) 135
44-IV-D-PP-39 (concesión) 74
44-IV-D-PP-40(42) 73
44-IV-D-PP-41(44) 110
44-IV-D-PP-42 (03) 114
44-IV-D-PP-43 (43B) 67
44-IV-D-PP-44 (45) 49
44-IV-D-PP-45 (Nuevo Rincón) 25
44-IV-D-PP-46 158
44-IV-D-PP-47 824
44-IV-D-PP-48 50
52-II-A-PP-01 29
52-II-A-PP-11(37) 98
52-II-A-PP-14 (38) 72
52-II-A-PP-15(41) 63
52-II-A-PP-16(40) 129
52-II-A-PP-22 32
52-II-B-PP-01(35) 142
52-II-B-PP-02(49) 1047
52-II-C-PP-18 (concesión) 124
52-II-C-PP-21 39
52-II-C-PP-22 8
52-II-C-PP-24 (Carretero) 13
POZO_A 102
POZO_B 135
PP-48-ADESA 6480
PP-51-ADESA 6480
PP-50-ADESA 6480
PP-55-ADESA 6480

Sergio Eduardo López Ramírez 22

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 14. Tabla de agua estado estacionario con pozos en funcionamiento

Ilustración 15. Vista 3D Tabla de agua costado oriental

Sergio Eduardo López Ramírez 23

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

Ilustración 16. Tabla de agua zona de mayor abatimiento

Ilustración 17. Vista tridimensional tabla de agua de la zona de mayor abatimiento

Sergio Eduardo López Ramírez 24

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

De otro lado, se encuentra la parte norte del acuífero, donde los abatimientos esperados son los menos
notorios, es decir, la tabla de agua no se muestra afectada en las mismas proporciones de la parte
oriental del campo de pozos de Corozal. Sin embargo, los descensos generales de la zona norte son
de 100 a 200 m, muy similares al promedio general del acuífero. Los principales pozos de grandes
caudales en ésta región, el Pozo A y el Pozo B, no afectan de manera dramática los niveles de agua.
Al igual que lo sucedido en la parte norte del acuífero, en la parte sur no se observan abatimientos
de la misma magnitud como los que se encuentran alrededor de la zona de Corozal. Para esta zona
el abatimiento más significativo se tiene alrededor del pozo 52-II-C-PP-24 (Carretero).

Sergio Eduardo López Ramírez 25

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

4 Conclusiones
 Debido a la heterogeneidad y anisotropía del acuífero, la calibración del modelo en estado
estacionario ha sido compleja.

 El modelo numérico propuesto tiene una calibración en estado estacionario aceptable,

representando de forma adecuada los niveles registrados para el año 1955, los cuáles, de
acuerdo a la información recopilada eran saltantes o muy cercanos a la superficie.

 A partir de la construcción de un modelo en estado estacionario con la entrada en

funcionamiento de los pozos de bombeo se representaron de forma aproximada los efectos
esperados por la entrada en funcionamiento de los nuevos pozos de gran magnitud y caudal.

 La principal afectación a la tabla de agua se encuentra al oriente del campo de pozos de

Corozal, donde los pozos 44-IV-PP-47, 44-IV-D-PP-46 y PP-50-ADESA (pozos de gran
profundidad y caudal) tienen influencia.

 El pozo 52-II-B-PP-02(49) tiene una afectación importante en el costado sur del Arroyo
Grande de Corozal y así mismo sobre los demás ubicados en el campo pozos de Corozal

 De los pozos de gran magnitud, aquellos ubicados en la parte norte del acuífero (Pozo A y
Pozo B) son los que tienen el menor impacto. Sin embargo, los abatimientos en esta zona no
dejan de ser significativos pues según la modelación se ubican en un rango de 100 a 200m.

 En la parte sur se presentan abatimientos sustanciales pero no de la misma magnitud que los
modelados para la parte norte del acuífero.

 Una vez se realicen las correcciones del software por parte de la empresa distribuidora del
mismo, la corrida del modelo en estado transitorio puede resultar satisfactoria. Lo anterior,
hará muy probable que el modelo sirva como herramienta para la toma de decisiones de la
autoridad ambiental en Sucre.

Sergio Eduardo López Ramírez 26

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

5 Agradecimientos
A los profesores Carlos Molano y Luis Alejandro Camacho por su apoyo en el desarrollo de este
trabajo de grado
Al departamento de Ingeniería Civil y Ambiental de la Universidad de los Andes por su apoyo en la
gestión de los diferentes trámites necesarios para la gestión de obtención de la información.
A la corporación autónoma de Sucre, CARSUCRE por su colaboración en la adquisición de
información de niveles, diseños de pozos, registros eléctricos y trabajos de investigación realizados
en el acuífero de Morroa. En especial, un agradecimiento al departamento de aguas liderado por el
Ingeniero Geólogo Héctor Mario Herrera Parra, quien en compañía de un excelente y amable grupo
de profesionales brindó toda su colaboración para la recopilación de información sobre el acuífero de
Morroa.
A todas las personas de Hidrogeocol por su colaboración sobre el manejo del software y su asesoría
en los puntos clave del modelo.
A la empresa Aguas de la Sabana (ADESA) por su ayuda sobre la información de registros y pruebas
de bombeo de los últimos pozos desarrollados en la zona.

Sergio Eduardo López Ramírez 27

Universidad de los Andes
Departamento de Ingeniería Civil y Ambiental
Actualización del modelo numérico del Acuífero Morroa utilizando Visual MODFLOW Flex

6 Referencias
Abreu & Diaz. (2004). Propuesta de un modelo hidrogeológíco conceptual para el acuífero de la
región de Sincelejo, Corozal, Morroa, Chingalé, Betulia y Sampúes (Departamentos de Sucre
y Cordoba, Colombia). Tesis para optar al titulo de Geólogo. Bogotá, Colombia: Universidad
Nacional de Colombia.
Arias & Fuquen. (1992). Descripción Geológica de la Formación Morroa. 1992: INGEOMINAS.
Buitrago, J., & Donado, L. (2000). Evalución de las Condiciones de Explotación de la Zona de
Recarga del Acuífero Morroa . Depratamentos de Sucre y Córdoba, Colombia: Proyecto de
Grado para optar al titulo de Ingeniero Civil. Facultad de Ingeniería. Universidad Nacional
de Colombia.
Cáceres, .., Cuacalón, I., & Kassem, T. (1967). Geología del cuadrángulo E8 Sincelejo. Bogotá:
Servicio Geológico Nacional.
CARSUCRE. (2005). Proyecto de Proteccción Integral de Aguas Subterráneas PPIAS. Sincelejo:
CARSUCRE.
CARSUCRE. (2014). Mapa Geológico del Acuífero Morroa. Sincelejo, Sucre.
DANE. (2010). SIGOT. Obtenido de SIGOT: http://sigotn.igac.gov.co/sigotn/default.aspx
DANE. (2014). SIGOT. Obtenido de SIGOT: http://sigotn.igac.gov.co/sigotn/default.aspx
Herrera, H., Vargas, M., & Taupin, J. (2006). Estudio Hidrogeológico con Énfasis en
Hidrogeoquímica de Acuifero Morroa (Colombia). Estudios de Hidrología Isotópica en
América Latina, 29-46.
Ingeominas. (2002). Mapa Hidrogeológico de Sucre. Mapa de Hidroquímica. Bogotá D.C, Colombia:
Ministerio de Minas y Energía.
Ingeominas. (2002). Memoria Técnica del Mapa de Aguas Subterráneas del Departamento de Sucre
en Escala 1:250.000. Bogotá: Ministerio de Minas y Energía.
Insfopal. (1981). Estudio Hidrogeológico del Flanco Nororiental de la Serranía de San Jacinto y de
la Zona Litoral del Golfo de Morrosquillo.
Pérez, A., Obregón, N., & García, O. (2009). Análisis de Metodologías Geoestadísticas en la
Modelación del Acuífero Morroa. Revista Facultad de Ingeniería Universidad de Antioquia,
63-76.
Rodriguez, G. (1993). Evaluación Hidrogeológica del Acuifero Morroa. Santa Fé de Bogotá:
INGEOMINAS.
Santisteban, H. (1989). Evaluación del potencial del acuífero de la formación Morroa. Tesis de Grado
para optar al título de Magister en Recursos Hidráulicos e Ingeniería Ambiental. Bogotá
D.C: Universidad de los Andes.
Superintendencia de Servicios Públicos Domiciliarios. (2011). SIGOT. Obtenido de SIGOT:
http://sigotn.igac.gov.co/sigotn/default.aspx

Sergio Eduardo López Ramírez 28