Evaluación Vulnerabilidad Aguasubterranea Botadero
Evaluación Vulnerabilidad Aguasubterranea Botadero
Evaluación Vulnerabilidad Aguasubterranea Botadero
PRESENTADO POR:
GONZALO DE JESUS MENA LOPEZ
JUAN ORLANDO MONTES GARCIA
RECTOR :
SECRETARIO GENERAL :
DECANO :
SECRETARIO :
DIRECTOR :
INGENIERO CIVIL
Título :
Presentado por :
Docentes Directores :
Docentes Directores:
iv
A nuestras familias, por el apoyo y comprensión recibida de parte de ellos en todo este
tiempo.
A nuestros asesores: Ing. Edgar Alfredo Gavidia Paredes e Ingra. Susan Elizabeth
Campos de Orellana; por su guía, trabajo, dedicación y consejos recibidos a lo largo de
la ejecución de este trabajo.
Y así a todas las personas no profesionales y profesionales, a las cuales se les consulto o
pidió ayuda, la cual recibimos sin ningún interés.
A mi Madre: Silvia López, quien siempre me ha apoyado con su gran amor y que sus
palabras y oraciones siempre me han respaldado.
A mi Padre: Francisco Mena, quien es un ejemplo de lucha, por su apoyo y por estar
cerca de mí en la distancia en esta etapa de mi vida.
De manera especial a Fina Rubio y a mí hermano Walter Mena, por sus palabras de
ánimo y apoyo recibido en esta ultima etapa de mi vida.
Y de una manera muy especial a cada profesor(a) que he tenido a lo largo de toda mi
formación académica, principalmente a Don Osmin Elías (QEPD), quien fue un gran
amigo y a quien recordare siempre.
A todas estas personas y otras que formaron parte de la historia de mi vida y a todas
aquellas quienes cayeron o siguen luchando por lo que creen…
A mis Padres: Orlando Montes Flores y Blanca Lidia García de Montes, por haberme
apoyado hasta hoy en día, a pesar de las dificultades enfrentadas en ocasiones en la
universidad siempre logre sobreponerme a ellas con su ayuda.
A mi Novia: Dina Lily Merino Ruiz, con tu amor y alegría siempre me alentaste a lo
largo del presente trabajo siendo tu ayuda de gran importancia.
A mi compañero de tesis: Gonzalo, que a pesar de que un inicio teníamos muchas dudas
respecto a este tema por los pocos conocimientos adquiridos de este en nuestra carrera
logramos finalizarlo con los objetivos propuestos.
A mis grandes amigos del Club Fifa: Adán, Aldana, Caleb, Giovanni, Miguel, Roberto y
Víctor, son las personas con las que paso mis ratos de ocio y los aprecio mucho a todos.
A mis amigas, Esmeralda Azahar, Paoly Aguiluz y Mónica Mata, siempre recibí buenos
consejos y palabras de aliento de su parte por lo que estimo mucho su amistad.
Y a todos aquellos que en el transcurso de mi vida han ido formando la persona que soy
hoy.
viii
2.3.1. Introducción a los SIG ...................................................................................... 70
2.3.2. Aplicación de los SIG en el manejo de los recursos hídricos
subterráneos ...................................................................................................... 77
2.3.3. Uso del ILWIS para la elaboración de mapas de vulnerabilidad de
acuíferos ............................................................................................................ 78
ix
4.4. Propuestas de medidas de seguimiento medioambientales y técnicas de
remediación de aguas subterráneas ................................................................. 152
x
INDICE DE FIGURAS
xi
Figura 3.11 Perfil geoeléctrico de un sector donde se ubica el antiguo botadero de
Mariona ................................................................................................................. 100
Figura 3.12 Ventana para importar datos en ILWIS 3.7 ........................................................... 106
Figura 3.13 Ventana de configuración para creación de un mapa de puntos a partir de
una tabla en ILWIS 3.7 .......................................................................................... 107
Figura 3.14 Mapa de puntos del perímetro de antiguo botadero de Mariona .......................... 108
Figura 3.15 Ventana para creación de mapa de polígonos a partir de mapa de puntos
en ILWIS 3.7 .......................................................................................................... 109
Figura 3.16 Mapa de polígonos de perímetro de antiguo botadero de Mariona ...................... 109
Figura 3.17 Ventana para conversión de formato vectorial a raster en ILWIS 3.7 .................. 110
Figura 3.18 Ventana para creación de georeferenciación en ILWIS 3.7 .................................. 110
Figura 3.19 Ventana para conversión de mapa de polígonos a mapa de segmentos en
ILWIS 3.7 ............................................................................................................... 111
Figura 3.20 Ventana para unir tablas en ILWIS 3.7 ................................................................. 113
Figura 3.21 Ventana para interpolación por método Kringing en ILWIS 3.7........................... 114
Figura 3.22 Ventana para intersección de mapas raster en ILWIS 3.7 .................................... 115
Figura 3.23 Ventana para creación de mapas de atributos raster en ILWIS 3.7 ...................... 116
Figura 3.24 Ventana para creación de dominio en ILWIS 3.7 .................................................. 118
Figura 3.25 Ventana de edición de dominio en ILWIS 3.7 ........................................................ 119
Figura 3.26 Ventana para aplicar dominio a mapa raster en ILWIS 3.7 .................................. 120
Figura 3.27 Mapa del parámetro D de GOD de la zona donde se localiza el antiguo
botadero de Mariona ............................................................................................. 121
Figura 3.28 Mapa de vulnerabilidad de las aguas subterráneas a la contaminación de
la zona donde se ubica el antiguo botadero de Mariona. Resultado de la
aplicación del método GOD ................................................................................. 123
Figura 3.29 Mapa de recarga potencial de la cuenca del rio Acelhuate .................................. 125
Figura 3.30 Mapa del parámetro D de DRASTIC de la zona donde se localiza el
antiguo botadero de Mariona................................................................................ 130
Figura 3.31 Ventana para crear georeferencia de una imagen en ILWIS 3.7 .......................... 131
Figura 3.32 Ventana de editor de georeferenciaciones en ILWIS 3.7 ....................................... 133
Figura 3.33 Ventana para creación de mapa de segmentos en ILWIS 3.7 ................................ 134
Figura 3.34 Ventana para poligonizar segmentos en ILWIS 3.7 ............................................... 136
xii
Figura 3.35 Mapa del parámetro R de DRASTIC de la zona donde se localiza el
antiguo botadero de Mariona................................................................................ 138
Figura 3.36 Mapa de parámetro S de DRASTIC de la zona donde se localiza el
antiguo botadero de Mariona................................................................................ 141
Figura 3.37 Mapa de parámetro T de DRASTIC de la zona donde se localiza el antiguo
botadero de Mariona ............................................................................................. 144
Figura 3.38 Mapa de vulnerabilidad de las aguas subterráneas a la contaminación de
la zona donde se ubica el antiguo botadero de Mariona. Resultado de la
aplicación del método DRASTIC.......................................................................... 146
Figura 4.1 Distribución de la vulnerabilidad de las aguas subterráneas a la
contaminación de la zona donde se ubica el antiguo botadero de Mariona.
Método GOD. ........................................................................................................ 148
Figura 4.2 Distribución de la vulnerabilidad de las aguas subterráneas a la
contaminación de la zona donde se ubica el antiguo botadero de Mariona.
Método DRASTIC. ................................................................................................. 149
Figura 4.3 Impacto de la expectativa de limpieza en los costes y duración de un sistema
convencional de bombeo y tratamiento. ................................................................ 163
Figura 4.4 Los líquidos de fase no acuosa no se disuelven bien en el agua subterránea. ........ 164
Figura 4.5 Un único pozo de extracción situado en x=0, y=0, en un acuífero con flujo a
lo largo del eje x. ................................................................................................... 166
Figura 4.6 Curva de zona de captura para un pozo único de extracción situado en el
origen, en un acuífero con velocidad de flujo v, espesor B y tasa de
bombeo Q. ............................................................................................................. 168
Figura 4.7 Superposición de la pluma en una curva de zona de captura para un solo
pozo de extracción ................................................................................................. 169
Figura 4.8 Ubicación de un único pozo de extracción por delante de una pluma de
contaminante ......................................................................................................... 170
Figura 4.9 Curvas de captura para pozos óptimamente espaciados a lo largo del eje y,
bombeando a una tasa Q: (a) dos pozos; (b) tres pozos. ...................................... 171
Figura 4.10 Ubicación de dos pozos de extracción en el frente de una pluma de
contaminante ......................................................................................................... 173
Figura 4.11 Manipulación del gradiente hidráulico con múltiples pozos. ................................ 175
xiii
Figura 4.12 Técnica de remediación de aguas subterráneas por Biorrecuperación in
situ. ........................................................................................................................ 178
Figura 4.13 Barrera reactiva permeable para el tratamiento de una pluma
contaminante en aguas subterráneas .................................................................... 179
xiv
INDICE DE TABLAS
xv
Tabla 4.2 Presentación de medida de seguimiento medioambiental 1 ...................................... 154
Tabla 4.3 Presentación de medida de seguimiento medioambiental 2 ...................................... 155
Tabla 4.4 Presentación de medida de seguimiento medioambiental 3 ...................................... 156
Tabla 4.5 Presentación de medida de seguimiento medioambiental 4 ...................................... 157
Tabla 4.6 Presentación de medida de seguimiento medioambiental 5 ...................................... 158
Tabla 4.7 Presentación de medida de seguimiento medioambiental 6 ...................................... 159
Tabla 4.8 Presentación de medida de seguimiento medioambiental 7 ...................................... 160
Tabla 4.9 Presentación de medida de seguimiento medioambiental 8 ...................................... 161
Tabla 4.10 Facilidad relativa de limpieza de acuíferos contaminados en función de la
química del contaminante y la hidrogeología del lugar .......................................... 165
xvi
INTRODUCCION
xvii
Luego en el segundo capitulo se presenta la teoría respecto a lo que son los acuíferos
como: Su definición, clasificación, características principales y contaminación de estos;
además de la definición e importancia de conocer la vulnerabilidad de las aguas
subterráneas en el manejo de recursos hídricos, presentando para ello los métodos GOD
y DRASTIC principalmente, ya que han sido los que se han aplicado. Finalizando este
con la teoría de los SIG y su aplicación en el manejo de recursos hídricos, enfocado en el
proceso de la elaboración de un mapa que muestre la vulnerabilidad de las aguas
subterráneas a la contaminación, en donde se habla sobre el uso del ILWIS para obtener
dicho mapa.
xvii
geoeléctrico resultado de un estudio previo. Posteriormente se muestra la información
obtenida sobre los parámetros del método GOD, luego se describe el procesamiento
realizado en el software ILWIS 3.7 para obtener el mapa de vulnerabilidad de las aguas
subterráneas a la contaminación de toda la zona de estudio; así mismo se presenta para el
método DRASTIC.
Finalizando con el quinto capitulo donde se mencionan las conclusiones a las se llego
luego de conocer y analizar los resultados obtenidos y las recomendaciones que se
proponen.
xviii
CAPITULO 1. Consideraciones Generales
1.1. Antecedentes
Mediante los parámetros solicitados en cada uno de los métodos y el uso del álgebra de
mapas se logra evaluar la vulnerabilidad de un acuífero a la contaminación. Las
metodologías antes mencionadas ya han sido utilizadas en trabajos de graduación
anteriormente realizados como:
2
acrónimo en ingles). Este estudio contó con el apoyo de instituciones como la
Administración Nacional de Acueductos y Alcantarillados (ANDA) y
Fortalecimiento Institucional de ANDA para Investigación de Aguas Subterráneas
(FIAS).
3
Además de los trabajos de graduación mencionados anteriormente que contienen
información acerca de los métodos GOD y DRASTIC, se cuenta con el trabajo de
graduación denominado “Evaluación Geoeléctrica del Antiguo Botadero de Residuos
Sólidos Urbanos de Mariona”, el cual fue realizado por alumnos de la Facultad de
Ciencias Naturales y Matemáticas de la UES en el 2009, que contiene entre otras
temáticas los sondeos geofísicos, inventario de pozos cercanos que muestran el nivel
freático y otros tipos de información acerca de nuestra zona de estudio.
4
contaminación por causa del ex botadero, debido a que además, éste no posee una
membrana impermeable en la parte inferior, convirtiéndolo en un foco de contaminación
para los recursos hídricos de la zona.
1.3. Objetivos
• Aplicar los métodos GOD y DRASTIC para evaluar la vulnerabilidad de las aguas
subterráneas en el sector donde se localiza el antiguo botadero de desechos sólidos
de Mariona.
5
• Proponer algunas medidas de seguimiento medioambientales que vayan enfocadas a
minimizar el daño causado por este antiguo botadero y algunas técnicas de
remediación de aguas subterráneas.
1.4. Alcances
Con la aplicación de los métodos GOD y DRASTIC y haciendo uso del programa
ILWIS 3.7, el cual es un software SIG gratuito, se obtendrán los mapas de
vulnerabilidad, uno por cada método y así determinar las zonas más sensibles al peligro
de contaminación por este antiguo botadero; además se compararan los resultados
obtenidos por los dos métodos utilizados para su posterior análisis y se propondrán
algunas medidas de seguimiento medioambientales que vayan enfocadas a minimizar el
daño causado por la presencia de este ex botadero y técnicas de remediación de aguas
subterráneas.
6
1.5. Limitaciones
• Este trabajo se hará en gran parte con información proporcionada por estudios ya
realizados y con la colaboración de instituciones públicas; por lo que nos
tendremos que condicionar a los trámites exigidos para que se nos proporcione la
información necesaria en la aplicación de los métodos a utilizar.
1.6. Justificación
En términos generales, las aguas subterráneas son de mejor calidad que las aguas
superficiales, las primeras representan un recurso importante para el desarrollo socio-
económico nacional; convirtiéndose en la principal fuente de abastecimiento de agua
potable en nuestro país. Su explotación, se realiza a través de la perforación de pozos y
mediante la captación de manantiales.
En los alrededores del antiguo botadero de Mariona, el vital líquido es usado por
asentamientos humanos y por industrias que lo utilizan en sus procesos y que a la vez
pueden ser posibles causantes de contaminación, como lo son las empresas: Coca Cola,
Jugos Jumex e Ingenio El Ángel (Rivera y Aparicio, 2009).
Considerando que el deterioro de la calidad de las aguas subterráneas suele ser lento y
zonal, permite que los problemas de alteración de la calidad se detecten transcurrido
7
bastante tiempo desde el inicio de la actividad contaminante, producto de las diversas
actividades que potencialmente generan una carga contaminante importante al subsuelo
y ocasionan un riesgo de contaminación a las aguas subterráneas.
8
CAPITULO 2. Generalidades
2.1. Acuíferos
2.1.1. Introducción
El agua procedente de la lluvia fluye hacia los ríos y arroyos, y vuelve a la atmósfera por
evaporación o transpiración o se filtra hacia el subsuelo para formar parte de las aguas
subterráneas. Cuando el agua se filtra a través de las grietas o los poros del terreno y las
rocas, pasa a través de la llamada zona insaturada, que se caracteriza por la presencia
tanto de aire como de agua en los espacios intersticiales del suelo. El agua vadosa, de la
zona insaturada o de aireación no está disponible para el aprovechamiento humano, es
decir, no se puede bombear, si bien las plantas ciertamente absorben el agua del suelo
que se encuentra cerca de la superficie. En la zona saturada los espacios intersticiales
están completamente rellenos de agua. El agua de esta zona se llama agua subterránea, y
su techo es el nivel freático. Hay una zona de transición entre las dos zonas mencionadas
anteriormente que es la franja de capilaridad, en donde el agua sube por las pequeñas
grietas como resultado de la atracción entre el agua y las superficies rocosas. La Figura
2.1 ilustra esta distribución (Gilbert M. Masters, 2008).
El hidrogeólogo debe tener una clara concepción del proceso de la naturaleza, que se
conoce con el nombre de ciclo hidrológico y que ya Meinzer (1942), consideraba el
concepto central de la hidrología (Custodio & Llamas, 2001). En la Figura 2.2 se
presenta un esquema simplificado del ciclo hidrológico, cuya comprensión no requiere
grandes explicaciones.
El agua subterránea constituye una porción del agua que circula a través del
ciclo hidrológico; esta proviene de la infiltración o percolación a través de la
superficie de la tierra, después de lo cual se transporta lentamente a grandes
distancias hasta que nuevamente regresa a la superficie por acción del flujo natural, de
las plantas o el hombre.
10
del total de agua dulce y 98.4% del agua dulce no congelada en la hidrosfera (Davis,
2005, p. 218), por lo que representa una fuente importante para el suministro de agua en
el mundo.
Las descargas del agua subterránea ocurren cuando el agua emerge; cuando fluye hacia
el mar o a la superficie en forma de fuentes, ríos y lagos. Dentro del ciclo
hidrológico el agua subterránea muy cercana a la superficie puede evaporarse en la
tierra y por la transpiración de la cobertura vegetal.
Podemos decir también que el agua subterránea es el agua que circula en la zona
saturada de los acuíferos y que es posible de captar por obras de ingeniería (pozos,
drenes, etc.). Un acuífero es una capa geológica saturada suficientemente permeable para
permitir el flujo de agua a su través. Se sienta sobre una base o lecho de confinamiento,
que es una capa relativamente impermeable que limita el movimiento del agua; cómo
podemos ver el concepto de aguas subterráneas esta estrechamente relacionado con el de
acuífero, temática que se abordara más adelante.
11
¾ Concepto de hidrogeología.
La hidrogeología es una rama de las ciencias geológicas, que estudia las aguas
subterráneas en lo relacionado con su circulación, sus condicionamientos geológicos y
su captación; así su definición dice “La hidrogeología es la ciencia que estudia el origen
y la formación de las aguas subterráneas, las formas de yacimiento, su difusión,
movimiento, régimen y reservas, su interacción con los suelos y rocas, su estado
(líquido, sólido y gaseoso) y propiedades (físicas, químicas, bacteriológicas y
radiactivas); así como las condiciones que determinan las medidas de su
aprovechamiento, regulación y evacuación” (Mijailov, 1985, p. 285).
3. El estudio de los procesos que rigen los movimientos de las aguas subterráneas en el
interior de las rocas y de los sedimentos.
12
¾ Flujo en medios porosos: Ley de Darcy
Experiencia de Darcy.
En 1856, en la ciudad francesa de Dijon, el Ingeniero Henry Darcy fue encargado del
estudio de la red de abastecimiento a la ciudad. Parece que también debía diseñar filtros
de arena para purificar el agua, así que se interesó por los factores que influían en el
flujo del agua a través de los materiales arenosos, y presentó el resultado de sus trabajos
como un apéndice a su informe de la red de distribución. Ese pequeño apéndice ha sido
la base de todos los estudios físico-matemáticos posteriores sobre el flujo del agua
subterránea.
Donde:
Q = Caudal
Δh = Diferencia de Potencial entre A y B
Δl = Distancia entre A y B
Gradiente hidráulico = Δh /Δl
13
Básicamente un permeámetro es un recipiente de sección constante por el que se hace
circular agua conectando a uno de sus extremos un depósito elevado de nivel constante.
En el otro extremo se regula el caudal de salida mediante un grifo que en cada
experimento mantiene el caudal también constante. Finalmente, se mide la altura de la
columna de agua en varios puntos (como mínimo en dos, como en la Figura 2.4).
i = Δh / L (2.1)
14
En laboratorio, el permeámetro se sitúa verticalmente y con el flujo ascendente para
facilitar la evacuación del aire contenido inicialmente en el material poroso. Es decir;
variando el caudal con un grifo y/o moviendo el depósito elevado, los niveles del agua
en los tubos varían. Podemos probar también con permeámetros de distintos diámetros y
midiendo la altura de la columna de agua en puntos más o menos próximos. Pues bien:
cambiando todas las variables, siempre que utilicemos la misma arena, entonces se
cumple con la relación de la Ecuación 2.2.
Δh (2.2)
Q = K .Seccion.
Δl
Si utilizamos otra arena (más gruesa o fina, o mezcla de gruesa y fina, etc.) y jugando de
nuevo con todas las variables, se vuelve a cumplir la ecuación anterior, pero la constante
de proporcionalidad lineal es otra distinta. Darcy concluyó, por tanto, que esa constante
era propia y característica de cada arena. Esta constante se llamó permeabilidad (K)
aunque actualmente se denomina conductividad hidráulica.
⎛ dh ⎞ (2.3)
q = −K ⎜ ⎟
⎝ dl ⎠
Donde:
q = Q / sección (es decir: caudal que circula por m2 de sección)
K = Conductividad Hidráulica
15
dh/dl= Gradiente hidráulico expresado en incrementos infinitesimales (el signo menos se
debe a que el caudal es una magnitud vectorial, cuya dirección es hacia los Δh
decrecientes; es decir, que Δh o dh es negativo y, por tanto, el caudal será
positivo).
Sabemos que en cualquier conducto por el que circula un fluido se cumple que:
Caudal
Velocidad Darcy =
Seccion total (2.5)
La parte de la sección total por la que puede circular el agua es la porosidad eficaz; si
una arena tiene una porosidad del 10% (0.10), el agua estaría circulando por el 10% de
la sección total del tubo. Y para que el mismo caudal circule por una sección 10 veces
menor, su velocidad debería ser 10 veces mayor. Por tanto, se cumplirá que:
Velocidad Darcy
Velocidad lineal media = (2.6)
ηe
Donde:
ηe = porosidad eficaz
16
Denominamos velocidad lineal media, y no velocidad real, al resultado de la Formula
2.6, debido a lo siguiente: esa fórmula refleja correctamente la velocidad real de las
partículas en una sección cualquiera del medio poroso, por ejemplo, en la mostrada en la
Figura 2.6, pero no es exacta para calcular con ella el tiempo de recorrido entre dos
puntos.
En la Figura 2.7 se muestra un tubo de longitud L1 lleno de arena por el que se hace
circular agua. Calculamos la velocidad lineal media mediante las Formulas 2.5 y 2.6, y
con esa velocidad evaluamos el tiempo de recorrido a lo largo del tubo de dicha figura
(tiempo=L1 /velocidad). Si después medimos experimentalmente ese tiempo de recorrido
añadiendo un colorante al agua, obtendríamos un tiempo ligeramente superior, ya que la
distancia recorrida ha sido mayor: no L1 sino L2 (que es desconocida).
17
Si llamamos velocidad real a la registrada a lo largo de un recorrido a través de un medio
poroso, ésta sería igual a:
• Acuífero
• Acuícludo
• Acuitardo
• Acuífugo
18
Acuícludo: Del latín claudere que significa encerrar o cerrar, se define como aquella
formación geológica que conteniendo agua en su interior, incluso hasta la saturación no
la transmite y por lo tanto no es posible su explotación. Dentro de este grupo podemos
incluir los cienos y legamos (arcillas por lo general) de origen deltaico y/o estuario. Que
a pesar de poseer enormes cantidades de agua (superiores al 50% en volumen) no son
hidrogeológicamente aptas para construcción de captaciones de aguas subterráneas
(Custodio & Llamas, 2001, p: 260).
Acuitardo: Del latín tardare que significa retardar, hace referencia a la existencia de
numerosas formaciones geológicas que, conteniendo apreciables cantidades de agua la
transmiten muy lentamente por lo que tampoco son aptos para el emplazamiento de
captaciones, pero sin embargo, bajo condiciones especiales permiten una recarga vertical
de otros ciertos casos. Por ejemplo un nivel de arcillas limosas o arenosas, puede
comportarse como un acuitardo, si esta dispuesto encima o debajo de un acuífero más
importante, al cual puede recargar, o incluso recibir agua del mismo (Custodio &
Llamas, 2001, p: 260).
Acuifugo: Del latín fugere que significa huir, son aquellas formaciones geológicas que
no contienen agua ni la pueden transmitir, como por ejemplo, un macizo granítico no
alterado, o unas rocas metamórficas sin apenas meteorización ni fracturación (Custodio
& Llamas, 2001, p: 260).
Los acuíferos los podemos clasificar por la presión hidrostática del agua encerrada en los
mismos, y dentro de esta clasificación podemos señalar:
19
Acuíferos libres no confinados o freáticos: Son aquellos en los cuales existe una
superficie libre del agua encerrada en ellos, que está en contacto directo con el aire y por
lo tanto a presión atmosférica (Custodio & Llamas, 2001, p: 264). La superficie libre se
conoce como superficie freática y el nivel a que ella se eleva, respecto a otro de
referencia, nivel freático. Está formado en general por un estrato permeable parcialmente
saturado de agua que yace sobre otro estrato impermeable o relativamente impermeable.
La ilustración de la Figura 2.8 muestra un acuífero libre, un pozo perforado hasta la zona
saturada de tal acuífero, tendrá agua a la presión atmosférica en el nivel freático. El agua
subterránea contenida en un acuífero cautivo, que es aquel encerrado entre dos capas de
confinamiento; puede estar bajo presión, de modo que en un pozo perforado en dicho
acuífero; el agua alcanzará una altura mayor del nivel a que se encuentra confinada, en
cuyo caso se denomina pozo artesiano. El nivel que alcanza el agua en un pozo artesiano
se denomina: nivel piezométrico. En algunas circunstancias el agua puede alcanzar la
superficie del pozo e incluso superarla sin necesidad de bombeo, debido a la presión del
acuífero confinado. En tal caso se llama pozo artesiano surgente. En la Figura 2.8
también se muestra una capa impermeable en medio de la zona no saturada, sobre la
20
masa de agua principal. El agua que se filtra hacia abajo es atrapada por esta capa y crea
un acuífero suspendido.
No todas las formaciones geológicas, o rocas en general, poseen, la misma facilidad para
transmitir y proporcionar agua en cantidades apreciables económicamente; es por ello
que los acuíferos tienen ciertas características que los diferencian de otras formaciones
geológicas. Estas características y propiedades permiten definir y en algunos casos,
predecir, el funcionamiento o respuesta del acuífero frente a unas determinadas acciones
exteriores (Custodio & Llamas, 2001).
Los intersticios de la roca, son como los caminos por los que pasa el flujo de agua
subterránea y se clasifican en:
Poros: son los intersticios intergranulares que hay entre los granos de los sedimentos
clásticos consolidados y no consolidados o en las tobas volcánicas sueltas (Figura 2.9A).
21
Fisuras: son fracturas o grietas en las rocas (Figura 2.9B) y Cavidades: presentes en las
rocas carbonatadas y en los tubos de lava de las rocas volcánicas (Figura 2.9C).
Como es evidente, no todas las formaciones geológicas poseen las propiedades que
caracterizan a un acuífero, es decir, que sus intersticios o espacios huecos, sean capaces
de almacenar el agua y de servirle de conductos a la vez para transmitirla, permitiendo
su aprovechamiento para satisfacer las necesidades del hombre. Desde este punto de
vista, el concepto de acuífero es en cierto modo relativo, ya que una formación
geológica que sea capaz, por ejemplo, de producir 1 l/s, no se considerará como
acuífero, en un lugar donde haya otras en que sea posible captar 50 l/s o más; pero en
una región árida, donde no haya otras posibilidades, a pesar del ínfimo caudal, esa
formación que brinda 1 l/s, para resolver las necesidades de una pequeña población,
constituye sin duda un acuífero, aunque pobre.
Las formaciones geológicas que constituyen acuíferos pueden incluirse en uno de estos
dos grupos:
22
Los acuíferos formados por depósitos no consolidados, están constituidos por materiales
sueltos, fundamentalmente, arenas, gravas o mezclas de ambas de origen geológico muy
diverso. Según Todd (1959) el 90% de los acuíferos en explotación en el mundo, caen
bajo esta categoría.
Porosidad
23
presión que se produce cuando el agua va siendo extraída, así para calcular el
rendimiento es necesario introducir otro término: el coeficiente de almacenamiento
(Masters y Ela, 2008).
Los poros pueden ser el resultado de procesos muy diferentes; unos (porosidad primaria)
se originan al mismo tiempo que las rocas en los que se encuentran, por ejemplo los que
existen entre los diversos fragmentos que constituyen una roca detrítica (arena de una
playa o de una duna, depósitos aluviales de un río, etc.); otros (porosidad secundaria) se
originan con posterioridad a la formación de la roca y a menudo por procesos
completamente diferentes a los que dieron origen a aquella, por ejemplo, fisuras o
diaclasas en una roca fracturada, o lo que quizás es mejor conocido, la red de conductos
y cavidades producida por disolución de rocas solubles (simas, sumideros o cuevas).
En la porosidad influyen varios factores, entre los que se puede distinguir los siguientes:
24
define como el volumen de agua contenido en una roca, liberada por la acción de la
gravedad y que se denomina Val (volumen de agua libre). La porosidad efectiva, esta
expresada por la Ecuación 2.9.
V al (2.9)
ηe = × 100
Vtot
Donde:
ηe = Porosidad efectiva
Val = Volumen de agua libre
Vtot= Volumen total
h (2.10)
Q = kA
I
La Ecuación 2.10 se cita aquí sólo para definir el concepto de permeabilidad y obtener
sus unidades: Despejando el valor de k en la fórmula anterior se comprueba que las
unidades de k son las de una velocidad (longitud/tiempo). En el Sistema Internacional
25
serían m/seg., pero para manejar números más cómodos, por tradición se continúa
utilizando m/día. En Geotecnia y otras ramas de ingeniería se utiliza el cm/ seg.
Los factores que determinan la permeabilidad pueden ser intrínsecos y extrínsecos. Los
intrínsecos son los propios del acuífero y dependen del tamaño de los poros. Los factores
extrínsecos son los que dependen del fluido y son fundamentalmente su viscosidad y su
peso específico. La viscosidad de un fluido es la medida de la fuerza resistente, por
unidad de área y por unidad de gradiente de velocidad transversal a la dirección del
movimiento del fluido. Tanto la viscosidad como el peso específico dependen de la
temperatura, por lo que en casos especiales debe tenerse en cuenta ésta, sobre todo por la
notable influencia de la viscosidad en la permeabilidad.
γ g (2.11)
kD = k= k
μ v
Y también:
v (2.12)
k= KD
g
Donde:
γ, peso específico absoluto del fluido (unidades de fuerza/longitud3)
μ, viscosidad dinámica del fluido (unidades de (fuerzaxtiempo)/longitud2)
g, aceleración de la gravedad (unidades de longitud/tiempo2)
v, viscosidad cinemática del fluido (unidades de longitud2/tiempo)
26
Por otra parte, el autor Pérez Franco (1978) ha demostrado que al considerar el flujo en
medios porosos en su forma más general no lineal, es necesario tomar en consideración
un nuevo parámetro adimensional característico de cada medio, o sea una nueva
propiedad intrínseca, que por analogía con el flujo en tuberías denomina rugosidad
equivalente, y que se representa por el símbolo “C”. Esta propiedad está relacionada con
la conductividad hidráulica turbulenta (KT) por la Fórmula 2.13.
1
⎛ 12 ⎞2
⎜ gk ⎟
KT = ⎜ ⎟ (2.13)
⎜ C ⎟
⎝ ⎠
Transmisividad
Concepto introducido por Theis en 1935, se define como el caudal que se filtra a través
de una franja vertical de terreno, de ancho unitario y de altura igual a la del manto
permeable saturado bajo un gradiente unitario a una temperatura fija determinada; es
decir que es el producto del espesor saturado del acuífero y la conductividad hidráulica
(Custodio y Llamas, 2001). Este coeficiente es representando por la letra T y puede
determinarse mediante ensayos de bombeo. Cuando el coeficiente es introducido a la
ecuación de Darcy (Ecuación 2.2), el flujo a través de cualquier sección transversal
vertical de acuífero viene expresado por (Johnson División OUP Inc, 1975):
Q = TIW (2.14)
27
Donde I es el gradiente hidráulico y W es el espesor del acuífero a través del cual tiene
lugar el flujo.
Si observamos la Figura 2.10 intuimos que los dos estratos acuíferos deben proporcionar
el mismo caudal: uno tiene la mitad de permeabilidad, pero el doble de espesor que el
otro.
Efectivamente, el parámetro que nos indica la facilidad del agua para circular
horizontalmente por una formación geológica, es una combinación de la permeabilidad
y del espesor, siendo para este ejemplo la transmisividad en ambos casos de 150 m2/día
Coeficiente de almacenamiento
Hay dos tipos de acuíferos desde el punto de vista de la presión hidrostática del agua
encerrada en los mismos: el acuífero freático, que tiene el techo a la presión atmosférica,
y el acuífero confinado, que tiene en su techo una capa confinante. Si disminuye la
presión hidrostática de un acuífero confinado, por ejemplo por la extracción de agua
subterránea, aumenta el peso a soportar por el acuífero y la presión que resulta, expulsa
cierta cantidad de agua del acuífero. Al mismo tiempo, el descenso de la presión produce
una pequeña expansión y consiguiente liberación de agua. La capacidad de producir
agua de un acuífero cautivo se expone mediante el coeficiente de almacenamiento.
28
El coeficiente de almacenamiento se define como el volumen de agua que puede ser
liberado por un prisma vertical del acuífero de sección igual a la unidad y altura igual a
la del acuífero saturado si se produce un descenso unidad del nivel piezométrico o de
carga hidráulica; (Custodio & Llamas, 2001) y se expresa como se muestra en la
Ecuación 2.15.
Podemos encontrar tanto características químicas como físicas del agua, cada una de
estas es importante en el estudio de las aguas subterráneas, ya que nos permite establecer
si un agua es adecuada para consumo humano, o para qué tipo de usos es que puede
llegar a ser beneficiosa en diversas actividades antrópicas. A continuación mencionamos
de forma resumida lo más importante referente a la química y física del agua
subterránea, iniciando con las principales sustancias disueltas que se encuentran
presentes en el agua subterránea.
29
Los iones menores más importantes que suelen estar en concentraciones entre 0.01 y 10
ppm son: anión nitrato (NO3-), ión carbonato (CO3=), catión potasio (K+), catión ferroso
(Fe++), ión nitrito (NO2-), ión fluoruro (F-), ión amonio (NH4+) y Ion estroncio (Sr++).
También hay que mencionar otros que se encuentran en concentraciones entre 0.0001 y
0.1 ppm, y los cuales son: ión bromuro (Br-), ión sulfuro (S=), catión fosfato (PO4-3), ión
borato diácido (BO3H2-), ión nitrito (NO2-), ión hidróxido (OH-), ión yoduro (I-), ión
férrico (Fe+++), ión manganeso (Mn++), ión amonio (NH4+), ión hidruro (H-), ión
alumínico (Al+++); entre otros.
Características físicas
Dentro de este tipo de características podemos mencionar (Custodio & Llamas, 2001):
30
• Conductividad y resistividad: Conductividad es la capacidad de un agua para
conducir electricidad y resistividad, es la reciproca de la conductividad. La
conductividad se puede medir por medio de un conductivímetro o
resistivímetro portátil.
31
Características químicas y fisicoquímicas del agua
Las características a mencionar en este apartado son (Custodio & Llamas, 2001):
32
• Demanda química de oxigeno o materia orgánica: Medida de la cantidad de
oxigeno necesario para eliminar la materia orgánica contenida en un agua
mediante procesos biológicos aerobios. Comúnmente se refiere a 5 días. Es una
medida importante de contaminación y debe referirse a un cierto tiempo.
En los últimos años ha tomado gran relevancia el estudio de las aguas subterráneas, ya
que estas disponen la mayor parte del recurso hídrico dulce que puede ser usado en
diferentes actividades por el ser humano; pero actualmente ya no sólo se trata de cómo
encontrar agua, este vital liquido se encuentra afectado por las diferentes actividades
humanas que han producido que la calidad de este ya no sea la misma que en el pasado,
y es por medio del estudio de acuíferos y estudio de las diferentes fuentes de
contaminación de éstos, que podemos llegar a entender en qué forma evitar que las
aguas subterráneas sigan siendo contaminadas o por lo menos predecir la evolución del
problema.
33
durante las últimas decenas de miles años; pero también podemos encontrar en menor
medida carbonatos de sodio y magnesio, que dan dureza al agua y se depositan en forma
de sarro en las superficies metálicas, losas, vidrios, etc.
9 Aguas residuales
9 Actividades agrícolas
9 Ganadería
9 Otras
Existen diferentes tipos de procesos contaminantes que pueden afectar las aguas
subterráneas, y entre ellas podemos mencionar (Figura 2.11):
34
Entre las actividades que causan contaminación puntual encontramos a:
infiltran en el terreno.
9 Explotación excesiva de los acuíferos que facilita el que las aguas salinas
invadan la zona de aguas dulces, por desplazamiento de la interfase entre los dos
tipos de aguas.
35
El tiempo con que se produce la contaminación obedece a varios factores, los cuales se
describen a continuación:
a) Dilución
La importancia de este proceso viene ligada a la forma de producirse la
contaminación; el agua contaminada es dispersada en el acuífero, lo cual supone
un cierto grado de mezcla de los contaminantes con el agua, por lo que el
ambiente físico y la estructura, mineralogía y características hidráulicas de los
acuíferos juegan un papel muy importante en la resistencia a la contaminación.
b) Filtración mecánica
En materiales granulares finos pueden retenerse partículas en suspensión y
microorganismos. En arenas gruesas, gravas o acuíferos fracturados, es casi
inexistente. Además, la capacidad de cualquier filtro es limitada.
c) Precipitación
Pueden precipitar diversas sustancias según los cambios en el pH o el potencial
redox del medio, al igual que la filtración, puede disminuir fuertemente la
permeabilidad del medio.
f) Desintegración o descomposición
Los elementos radiactivos y algunos productos químicos, como los pesticidas,
desaparecen con el paso del tiempo, se habla de desintegración radiactiva o de
descomposición, en el caso de los pesticidas. En ambos, el factor crucial es la
vida media (tiempo en el que la concentración de la sustancia se reduce a la
mitad). Esto será efectivo en elementos cuya vida media sea corta en relación con
el tiempo de tránsito del agua en el acuífero. La vida media de algunos pesticidas
es de unos pocos días, mientras que la de ciertos elementos radiactivos es de
miles de años.
La principal solución correctiva que debería de tomarse para la gestión de los residuos
sólidos, es la creación de un Relleno Sanitario que cumpla con todas las normas técnicas
necesarias; donde debe de haber un manejo integral de los residuos sólidos, desde su
generación hasta su disposición final; y también la población debe tomar conciencia que
es muy importante el reciclaje, para no desperdiciar productos que necesitan de un gran
tiempo para su biodegradación.
38
• Transmisión de organismos patógenos de animales infectados al hombre, por
contacto con el suelo, alimentos, agua y por la crianza de animales alimentados
en estas condiciones.
39
• Riesgos a la salud de los trabajadores del botadero.
Las poblaciones expuestas a los agentes físicos, químicos y biológicos de los residuos
sólidos municipales son los trabajadores formales e informales que manipulan residuos;
la población no servida; la que vive cerca de los sitios de tratamiento y disposición de
los residuos sólidos; la población de segregadores (personas que se dedican a la
segregación de la basura) y sus familias; y la población en general, a través de la
contaminación de los cuerpos de agua superficiales y subterráneos, del consumo de
carne de animales criados en basurales, y de la exposición a residuos peligrosos. Los
principales factores que contribuyen a esta situación son la poca atención de las
autoridades relacionadas con el sector y la deficiente calidad de servicios prestados
(Acurio, 2007).
Vrba y Zaporzec (1994) definen a la vulnerabilidad como una propiedad intrínseca del
sistema de agua subterránea que depende de la sensibilidad del mismo a los impactos
humanos y naturales.
40
Foster e Hirata (1991) dicen que la vulnerabilidad del acuífero a la contaminación,
representa su sensibilidad para ser adversamente afectado por una carga contaminante
impuesta.
Zektser (2000) dice que el concepto de vulnerabilidad del agua subterránea se basa en
el supuesto de que el medio ambiente geológico puede proporcionar un cierto grado de
protección a las aguas subterráneas naturales de un impacto antropogénico, la
contaminación de la atmósfera particularmente, el suelo y la escorrentía superficial.
41
Pettyjohgn et al (1991) dice que la susceptibilidad del acuífero se relaciona con la
posibilidad de contaminación, es decir, los acuíferos que tienen un alto grado de
vulnerabilidad son los que están en zonas de alta densidad poblacional.
Cuando tenemos un recurso finito que puede llegar a ser afectado por diferentes razones,
resulta necesario realizar una planificación de cómo se puede usar tal recurso y de qué
forma se puede proteger éste también.
42
El agua es un vital recurso que se ve afectado por la demanda de la población para
diversas actividades antropogenias, a la vez que se ve afectado por ciertas actividades
contaminantes que limitan el uso futuro del mismo.
Es en este contexto que surge la denominada Gestión Integral del Manejo de los
Recursos Hídricos, en nuestro caso priorizando las aguas subterráneas, sabemos que en
el país es la mayor fuente de abastecimiento de agua potable a la población; también es
conocido que los acuíferos tienen cierto nivel de autodepuración, dependiendo del tipo
de roca que los resguarde. Es muy favorable la abundancia de arcillas y de materia
orgánica, en cambio en los depósitos aluviales o las zonas kársticas la purificación del
agua es mucho más difícil y estos tipos de acuíferos son mucho más sensibles a la
contaminación; pero a pesar de la capacidad de purificarse es necesario tomar las
medidas cautelarias para no tener que recurrir a procesos muy difíciles y de alto costo
económico, por lo que cada vez existe mayor conciencia sobre la necesidad de
desarrollar una gestión integral segura de las aguas subterráneas. Como resultado de esta
situación, actualmente se generaliza el uso de técnicas para cartografiar la vulnerabilidad
de acuíferos a la contaminación, como una herramienta de gran utilidad para
compatibilizar la carga contaminante generada por actividades antrópicas.
43
• El conocimiento de la vulnerabilidad de los acuíferos es clave en la formulación
de políticas de gestión del agua subterránea y protección de la misma; así como en el
desarrollo de normativa sobre vertidos y desechos sólidos.
44
• Sistemas Paramétricos: Utilizan parámetros convenientemente seleccionados para
representar la vulnerabilidad, a los cuales les asignan distintos rangos e
interacciones para generar un índice de vulnerabilidad absoluto o relativo
Por otra parte podemos evaluar la vulnerabilidad de las aguas subterráneas por medio de
modelos de simulación, los que son muy difíciles de realizar en países en vías de
desarrollo, debido a la gran cantidad de información necesaria. Sin embargo, este
método es muy apropiado para analizar situaciones puntuales o con poca extensión
espacial. Los métodos estadísticos o de monitoreo no permiten determinar la
vulnerabilidad, sino que se utilizan para realizar estudios probabilísticas, que evalúen la
posible dispersión del contaminante. Debido a lo anterior, resulta conveniente el uso de
los métodos paramétricos, para los cuales existen experiencias a nivel tanto nacional
como internacional que avalan su correcto funcionamiento.
45
Cada uno de los parámetros tiene un valor de asignación dependiendo de sus
características, el valor de vulnerabilidad es el resultado de la multiplicación de estos
parámetros, podemos entenderlo mejor mediante la siguiente fórmula (Cerros, 2007):
V = GxOxD (2.16)
La estimación del índice de vulnerabilidad GOD involucra una serie de etapas concretas
(Foster et al., 2002), siendo éstas las siguientes:
46
1. Identificar el grado de confinamiento hidráulico del acuífero y asignarle un
índice a este parámetro en una escala de 0.0 a 1.0.
47
Figura 2.12 Método GOD para la evaluación de la vulnerabilidad a la contaminación de acuíferos.
Fuente: Gloria Cerros. Un enfoque de gestión ambiental del recurso hídrico desde la perspectiva de la
vulnerabilidad de acuíferos. 2007. p:69
48
depósitos diferentes se debe seleccionar la litología predominante o limitante para
caracterizar al sustrato suprayacente (Alvarado, 2007).
A continuación se explicara por separado en qué consiste cada uno de los parámetros a
determinar en la metodología GOD para evaluar la vulnerabilidad de las aguas
subterráneas a la contaminación.
• Ninguno
• Surgente
• Confinado
• Semiconfinado
• No confinado (cubierto)
• No confinado
49
Litología de la cobertura (O)
Este es igual que el parámetro I del método DRASTIC, y se refiere al tipo de litología de
la zona no saturada del acuífero, y al grado de consolidación de los materiales que lo
constituyen. Estos valores varían entre 0.4 y 1.0; y para la obtención de este parametro
es necesario contar con información hidrogeológica.
Para este índice se toma en cuenta la profundidad a la cual se encuentra el nivel del agua
subterránea en los acuíferos no confinados y al techo de éste cuando es confinado, el
índice varia de 0.6 a 1.0 (Alvarado, 2007).
J.M. Vías et al (1999) manifiesta que este método fue desarrollado por Aller en 1987, es
utilizado para cartografiar la vulnerabilidad a la contaminación en acuíferos, forma parte
de los métodos nominados de sistemas paramétricos, se aplica independientemente del
contaminante potencial y permite trabajar con unidades espaciales homogéneas. En este
método se consideran 7 variables cuyas iniciales (en ingles) dan nombre al método:
50
C: Hydraulic Conductivity of the aquifer (Conductividad Hidráulica)
A cada una de estas variables se les asigna un puntaje (los parámetros con subíndice R
en la Ecuación 2.17 indican el puntaje en cada variable), índices de 1 a 10, los cuales
reflejan su importancia en la vulnerabilidad; también se les asigna un peso o
ponderación (los parámetros con subíndice W en la Ecuación 2.17 indican el peso o
ponderación en cada variable), de acuerdo a la influencia respecto a la vulnerabilidad.
Para el peso ponderado se emplea índices entre 1 y 5, adoptando generalmente el mayor
de 5 para la profundidad del agua (D) y la litología de la sección subsaturada (I) y el
menor de 1 para la topografía (T); posteriormente ambos índices se multiplican y luego
se suman los 7 resultados, esto lo podemos visualizar de mejor manera en la siguiente
ecuación:
V = DR DW + RR RW + AR AW + S R SW + TRTW + I R I W + CR CW (2.17)
51
Foster & Skinner (1995) argumentan que el índice de vulnerabilidad obtenido a través
del método DRASTIC, es el resultado de la interacción de varios parámetros, cuyas
ponderaciones son cuestionables y en ciertos casos no independientes, pudiendo dar
origen a situaciones en que parámetros relevantes, tales como la litología del suelo,
quedan ocultos por otros parámetros menos importantes, tales como la movilidad del
contaminante en la zona saturada.
Este método se recomienda sea aplicado en zonas de estudio superiores a los 42.7 Has,
lo que permite una evaluación del sistema de flujo regional más que local; también se
puede usar en sitios para disposición de desechos sólidos, porque permite evaluar la
susceptibilidad de un sitio a la acción de una carga contaminante permanente; también
para la evaluación de zonas de recarga y para creación de zonas de protección de
acuíferos, porque permite discriminar áreas homogéneas con el mismo nivel de
vulnerabilidad (Belmonte, 2000).
A continuación se explica en qué consiste las distintas variables del método DRASTIC,
de qué forma se puede obtener información necesaria sobre ellas, y sobre los rangos y
respectivos puntajes que se les pueden asignar a cada una de las variables.
La información a emplear en este parámetro debe estar en metros, y puede haber sido
obtenida por medio de pozos perforados y excavados observados en visitas de campo a
la zona en estudio, también mediante informes previos, donde puede ser necesaria una
verificación (Escobar, Garcia, & Guerrero, 2002).
52
La medición de este parámetro puede dar resultados erróneos cuando el acuífero que se
estudia es semiconfinado o confinado, ya que el nivel del agua medido no corresponde a
la distancia desde el suelo hasta el acuífero; por lo que es mejor trabajar solo con
información de pozos excavados (Escobar, Garcia, & Guerrero, 2002).
Tabla 2.4 Puntajes de parámetro de Profundidad al nivel freatico en método DRASTIC (Dr)
Rango de Profundidades (m) Valor
<1.5 10
1.5-5.0 9
5.0-10.0 7
10.0-20.0 5
20.0-30.0 2
>30.0 1
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 43
Recarga (R)
La recarga es el volumen de agua por unidad de área que ingresa al acuífero durante un
período de tiempo, para lo cual debe atravesar la superficie y la zona no saturada. En su
camino el agua arrastra consigo muchas sustancias, siendo así el principal vehículo para
la lixiviación y el transporte de contaminantes líquidos y sólidos (Cerros, 2007).
En la Tabla 2.5 se describe los puntajes a asignar dependiendo del rango del rango en
que se encuentre los valores obtenidos de la recarga neta.
53
Tabla 2.5 Puntajes de Recarga Neta en método DRASTIC (Rr)
Rango (m/día) Valor
0 – 0.00014 1
0.00014 - 0.00028 3
0.00028 – 0.00049 6
0.00049 – 0.00070 8
> 0.00070 9
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 43
Material acuífero(A)
En la Tabla 2.6 se presenta los puntajes a asignar dependiendo del tipo de la litología
presente en el acuífero en estudio.
54
Tabla 2.6 Puntajes para Tipo de Acuífero en método DRASTIC (Ar)
Tipo de Acuífero Rango Valor Típico
Arcillas, limos, margas 1-3 2
Rocas ígneas, metamórficas 2-5 3
Rocas ígneas, metamórficas 3-5 4
alteradas
Alternancia de areniscas, arcillas 5-9 6
y calizas
Areniscas masivas 4-9 6
Calizas masivas 4-9 6
Arenas, gravas y conglomerados 4-9 8
Volcánicas 2-10 9
Calizas carstificadas 9-10 10
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 43
Tipo de suelo(S)
El presente parámetro al igual que el medio acuífero, los materiales que conforman esta
zona influyen en las propiedades del contaminante, atenuándolo, con lo que la
vulnerabilidad a partir de este parámetro depende del tipo de materiales que se
encuentran en ella (Escobar, Garcia, & Guerrero, 2002).
En la Tabla 2.7 se describe los puntajes que se deben usar dependiendo del tipo de suelo
que exista en la zona en estudio.
55
Tabla 2.7 Puntajes para Tipo de Suelo en método DRASTIC (Sr)
Tipo de Suelo Valor
Arcilla no expansiva agregada 1
Suelo orgánico 2
Marga arcillosa 3
Marga limosa 4
Marga 5
Marga arenosa 6
Arcilla expansiva y/o agregada 7
Turba 8
Arena 9
Grava 10
Delgada o ausente 10
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 43
El agua que llega a la superficie del suelo, puede escurrir o infiltrarse. La influencia de
la topografía consiste en que permite que un contaminante pueda infiltrarse a través de la
zona no saturada (Cerros, 2007).
En la Tabla 2.8 se muestra los puntajes a usar dependiendo del rango en que se
encuentre la pendiente en la zona en estudio.
56
Impacto de la zona no saturada según su litología (I)
El tiempo que demora el agua en recorrer el camino hasta la zona saturada no sólo
depende de la longitud, sino también del tipo de material con que se encuentra. Las
características del espesor, composición litológica, propiedades químicas de las rocas
determinan las condiciones de atenuación del material que se encuentra por debajo del
horizonte del suelo y por encima de la capa freática. La biodegradación, la
neutralización, la fricción, reacción química, volatización y dispersión son fenómenos
que pueden ocurrir y contribuyen al proceso de atenuación en la zona vadosa
(Hernández, 2006).
Tabla 2.9 Puntajes para parámetro Impacto de la Zona Vadosa en método DRASTIC (Ir)
Formación de la zona vadosa Valor
Arcilla, limo. Marga 1
Esquistos, pizarras 3
Calizas 6
Areniscas 6
Alternancias: calizas, areniscas y arcillas 6
Arenas y gravas arcillosas 6
Rocas: ígneas metamórficas 4
Arenas, gravas 8
Volcánica 9
Calizas carstificadas 10
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 44
57
Conductividad Hidráulica (C)
Este parámetro debe ser medido en metros sobre día, y la información es obtenida de
informes de pozos perforados o excavados, además de complementarse de estudios
hidrogeológicos previos (Escobar, Garcia, & Guerrero, 2002).
En la Tabla 2.10 se presenta los puntajes a usar dependiendo del rango en que se
encuentre la conductividad hidraulica determinada en el acuífero en estudio.
Tabla 2.10 Puntajes utilizados para la Conductividad Hidráulica en método DRASTIC (Cr)
Conductividad (m/día) Valor
0-4 1
4 - 12 2
12-28 4
28-41 6
41-82 8
> 82 10
Fuente: Mario Escobar, et al. Metodología para la evaluación de la vulnerabilidad de las aguas
subterráneas aplicada en el Valle de Zapotitán. 2002. p: 45
En Tabla 2.11 se da una descripción general sobre otros métodos que se pueden llegar a
usar para evaluar la vulnerabilidad de las aguas subterráneas.
58
Tabla 2.11 Resumen de métodos para evaluar la vulnerabilidad de las aguas subterráneas
METODO DESCRIPCION
SINTACS Auge (2006) nos explica que SINTACS es un método derivado de DRASTIC,
desarrollado por Civita et al en 1990, para adecuarlo a las características
hidrogeológicas de Italia, el acrónimo (en italiano) SINTACS comprende:
S: Profundidad del agua
I : Infiltración
N: Sección subsaturada
T: Tipo de suelo
A: Características hidrogeológicas del acuífero
C: Conductividad hidráulica
S: pendiente topográfica
AVI: Acrónimo de aquifer vulnerability index (índice de vulnerabilidad de acuifero),
desarrollado por Van Stempvoort et al en 1992, para el mapeo de la vulnerabilidad del
agua subterránea de la provincia de Prairie en Canadá. Cuantifica la vulnerabilidad a
través de un parámetro llamado resistencia hidráulica (C) función del espesor de la
zona subsaturada (d) y la permeabilidad vertical de los componentes de la misma (k).
En la siguiente formula podemos visualizar lo dicho anteriormente: [Auge, 2006]
di
C = ∑
i = estratos ki
(2.18)
EKv Es un método creado para ser usado en acuíferos libres, y cuyos parámetros a
determinar son la profundidad de la superficie freática (E) y la permeabilidad vertical
en la zona subsaturada (Kv), parámetros que también son considerados en el método
AVI; se les asigna índices que van de 1 a 5, posteriormente ambos valores se suman y
nos brindan un índice final, con extremos de 2 y 10, lo escrito lo podemos visualizar en
la siguiente fórmula: [Auge, 2006]
EKV = E + K V (2.19)
EPIK Método paramétrico desarrollado por Doerfliger y Zwahlen en 1997, para ser usado
específicamente en acuíferos kársticos, el acrónimo significa: Zona de fisuras y grietas
en terreno kárstico (Epikarst, E), Cubierta de protección (Protective cover, P),
Condiciones de infiltración (Infiltracion conditions, I), Desarrollo de la red kársticas
(Karst network development, K); el método incluye 4 factores de ponderación (α, β, γ
y δ), aplicables a cada parámetro, y para valorar su peso relativo en el cálculo del
índice de vulnerabilidad intrínseca, podemos representarlo en la siguiente fórmula:
[Auge, 2004]
VULNERABILIDADi = (α .Ei) + ( β .Pi) + (γ .Ii) + (δ .Ki) (2.20)
ΔhT’ Auge (2004), nos explica que este método fue desarrollado para acuíferos
semiconfinados, la vulnerabilidad de estos esta controlada por propiedades físicas y
geométricas del acuitardo que conforma su techo y también por la diferencia de
potencial hidráulico que guarda con el libre sobrepuesto; donde las variables a
determinar serán: Δh: diferencia de potencial hidráulico, T’: transmisividad vertical del
acuífero. Considerando las dos variables mencionadas (Δh y T') se establecen tres
grados de vulnerabilidad, determinados en primer lugar por el gradiente vertical de
potenciales hidráulicos y secundariamente por la transmisividad del acuitardo
59
Continuación Tabla 2.11…
GODE Auge et al (2004) desarrolla una modificación al método GOD, para usarlo en
acuíferos confinados o semiconfinados; la diferencia entre GODE y GOD es que
además de incluir los parámetros tipo de acuífero (G), litología de la cobertura (O),
profundidad del acuífero (D); el primero agrega una nueva variable a determinar, la
cual es el espesor de la zona saturada correspondiente al acuífero libre sobrepuesto y el
espesor del acuitardo (E).
GODS Desarrollado por Foster et al, en 1996, el presente método utiliza los mismos
parámetros que el índice GOD que son tipo de acuífero (G), litología de la cobertura
(O), profundidad del acuífero (D), surge debido al consenso de algunos investigadores
en que al método GOD le faltaba incluir el efecto que tiene el suelo (S), y se representa
mediante la siguiente formula: [Ríos y Vélez, 2008]
VULNERABILIDADGODS = G * O * D * S (2.21)
PATHS Método creado por Valcarce y Rodríguez en 2004, cuyo índice se obtiene a partir del
análisis de factores como la Precipitación (Pp), profundidad del nivel freático (P),
litología del acuífero (A), tectónica (T), pendiente topográfica (PT) y la resistencia
hidráulica del suelo (HS). Cada uno de los parámetros es ponderado por un factor W,
según su importancia. La valoración del índice se hace según la siguiente ecuación:
[Ríos y Vélez, 2008]
VULNERABILIDADPATHS = WP + WPp + WA + WT + WPT + WHS (2.22)
BGR Creado por Hölting et al en 1995, el método estima la protección de los materiales de
la zona no saturada al acuífero, los parámetros que son necesarios determinar en la
evaluación de la protección (PT) son: la capacidad de campo efectiva del suelo (S),
recarga natural del acuífero (W), litología (R), espesor de cada una de las capas de la
zona no saturada (T), presencia de acuíferos colgados (Q), condiciones de presión (HP).
A cada uno de estos parámetros se les asocia un puntaje, los cuales ponderados y
sumados entre si, dan como resultado un valor que se traduce en vulnerabilidad[Ríos y
Vélez, 2008]
⎛ N ⎞
PT = W * S + W ⎜ ∑ R S * Ti ⎟ + Q + H P (2.23)
⎝ i =1 ⎠
9 Información de pozos
60
9 Modelo de elevación digital (MED)
9 Hidrogeología
9 Información de suelos
Perforación de pozos
El método más seguro para conocer las características de las formaciones que yacen por
debajo de la superficie terrestre del terreno, es el de perforar a través de éstas,
obteniendo de este modo muestras geológicas mientras se perfora y llevando un registro
litológico del agujero (Johnson División OUP Inc, 1975).
Las muestras que usualmente se obtienen por otros sistemas, son los cortes que resultan
de la acción del barreno, aunque no son completamente representativas de las formación
que se ha penetrado, y por lo general bastante confiables en las exploraciones de aguas
subterráneas (Johnson División OUP Inc, 1975, p. 180).
61
El método convencional de rotación resulta a menudo más rápido y más económico
cuando se perfora en formaciones no consolidadas, que no contienen guijarros o cantos
rodados. La ventaja del método de rotación crece conforme se perfora más
profundamente en tales materiales. Una de las desventajas más generalizadas en el
método de rotación, es la que en un agujero de prueba perforado por el método, el nivel
estático del agua por lo general no se puede medir a menos que se instale un ademe y
extraiga gran parte del lodo de perforación (Johnson División OUP Inc, 1975, p. 181).
Pozos de observación
62
b) El valor aproximado de su permeabilidad o transmisividad y de su coeficiente de
almacenamiento
Inventario de pozos
Métodos geofísicos
Dentro de los métodos eléctricos el que más utiliza es el método eléctrico resistivo y la
característica física utilizada es la resistividad o resistencia específica de los diferentes
63
terrenos. Si en un terreno introducimos una corriente eléctrica, ésta se propaga de forma
tal que podemos asignar al terreno una resistencia, una capacidad y una inductancia. Si
la corriente eléctrica es continua tendremos la resistencia como único parámetro
(Custodio & Llamas, 2001).
Los instrumentos portátiles sismográficos utilizan golpes sobre la superficie del terreno,
para producir las ondas de choque. El ensayo se efectúa colocando el detector en
contacto con el terreno, conectando un cable eléctrico con el contador de tiempo y dando
golpes con un mazo a varias distancias del detector. Se comienza en un punto a 3 m del
detector y se prosigue con los puntos siguientes de impacto, a intervalos de 3 m
(Johnson División OUP Inc, 1975, p.201).
Estudios hidrogeológicos.
Mapas topográficos
Los mapas topográficos nos permiten determinar las curvas de nivel de la zona en
estudio, es por medio de esta y la interpolación en un SIG que se logra determinar el
MED. La topografía ilustra la distribución de zonas montañosas y llanuras. Permite
apreciar la extensión de las distintas cuencas de recepción de cada uno de los cursos de
aguas superficiales que existen en la misma, la distribución de las pendientes del terreno
etc.
Mapas geológicos
Los mapas topográficos, mapas geológicos y las fotografías aéreas; constituyen una de
los instrumentos básicos para el estudio de las aguas subterráneas. Esta fuente de
información es de la mas básica y es por medio de esta que se puede determinar la
localización y movimiento de las aguas subterráneas (Custodio & Llamas, 2001).
65
Las fotografías aéreas constituyen un excelente instrumento para la cartografía geológica
pues, pone en relieve aspectos difíciles de ver sobre el propio terreno (Custodio &
Llamas, 2001).
Mapas de suelos
66
Estudios climatológicos
Pruebas de infiltración
Visitas de campo
Las visitas de campo son importantes para la verificación de datos obtenidos existentes y
generar un criterio propio sobre la problemática encontrada; en esta se puede delimitar la
zona de estudio mediante la ubicación de coordenadas geográficas que delimiten el
perímetro de donde se realizara la investigación.
67
2.2.5. Uso de los mapas de vulnerabilidad
Entre los años 1960 – 1970, se desarrolló la construcción de los mapas de vulnerabilidad
en diferentes países de Europa, fundamentalmente en Francia, Checoslovaquia y
Alemania. En los Estados Unidos, Walker construyó un mapa de contaminación
potencial de acuíferos en el estado de Illinois, (Walker, 1969), el cual es considerado el
primer mapa de vulnerabilidad producido en ese país.
Para realizar estos mapas se utiliza un Sistema de Información Geográfica (SIG), que es
una herramienta diseñada para consultar, manipular y aplicar los datos de: información
geográfica, cartografía digital, bases de datos de atributos y estadísticas, de una manera
simultánea y automatizada (López y Granados, 2000).
Los mapas de vulnerabilidad deben ser vistos como una de las herramientas principales
para el manejo del medio ambiente y pueden ser usados con tres propósitos
fundamentales (Baez, 2001):
68
a. En procesos de toma de decisión concerniente a la protección y manejo de los
recursos hídricos subterráneos,
En general hay que tomar en cuenta que los mapas de vulnerabilidad se ven limitados
por la falta de datos representativos debido a malas técnicas de recopilación de datos,
descripción inadecuada del sistema, o por falta de metodología generalizada; es por ello
que es conveniente no usar demasiados datos para evitar una mala interpretación
(Escobar, Garcia, & Guerrero, 2002).
Según Vrba y Zaporozec (1994), las tendencias futuras en la producción y uso de los
mapas de vulnerabilidad son las siguientes:
b. Definir las clases de vulnerabilidad con mayor precisión cuantitativa sobre la base
del incremento del conocimiento del transporte de contaminantes.
69
e. Incremento de la producción de mapas de vulnerabilidad a escalas 1:10000 a
1:25000.
Martínez et al (1999) define los SIG como una colección de hardware, software, datos
geográficos y personal integrado para la captura eficiente, almacenamiento,
actualización, manipulación y análisis de información geográfica (Fígura 2.13).
70
Figura 2.13Definición grafica de Sistemas de Información Geográfica
Fuente: Curso de Sistemas de Información Geográfica con aplicaciones al Manejo de Recursos Naturales
71
9 Manipulación de mapas con operaciones algebraicas (por ejemplo: sumas,
multiplicaciones, restas, y división entre mapas).
9 Operaciones de geoprocesamiento.
Los elementos geográficos básicos usados por los SIG se describen en la Tabla 2.12.
72
¾ Modelos de representación de los SIG
SIG Raster
Los SIG raster se centran en las propiedades del espacio más que en la precisión de la
localización. Dividen el espacio en celdas regulares donde cada una de ellas representa
un único valor. Los datos raster se componen de filas y columnas de celdas, cada celda
almacena un valor único, los datos raster pueden ser imágenes (imágenes raster), con un
valor de color en cada celda (o píxel). Otros valores registrados para cada celda pueden
ser un valor discreto, como el uso del suelo, valores continuos, como temperaturas, o un
valor nulo si no se dispone de datos.
9 Utilizan celdas a las cuales se les asigna el valor del atributo (otros nombres:
sistema cuadrícula o matricial).
SIG Vectorial
9 Utiliza coordenadas x, y para definir puntos, los puntos forman líneas y las líneas
cerradas, polígonos.
74
Figura 2.14 Tipos de representación de la información en los SIG
Fuente: http://www.geotecnologias.com/Documentos/GIS.pdf
¾ Programas SIG
75
Tabla 2.13 Características de algunos programas SIG
Modelo de
Fecha de Sistemas Licencia de
Nombre datos Descripción
creación operativos programa
usado
Programa que trabaja junto a
ESRI (Enviromental Systems
Research Institute); este
organiza, visualiza,
Programa no
ERDAS 1978 Windows Raster manipula,mide e integra
libre
cualquier tipo de imágenes
geográficas e información
geoespacial en ambientes de
2D y 3D.
Distribuido por ITC Enschede
(International Institute for
Geo-Information Science and
Earth Observation) Permite
entrar, manipular, analizar y
presentar datos geográficos
Raster (mapas). A partir de estos
ILWIS 1985 Windows Libre
Vectorial se puede generar
información espacial,
modelos espaciales y evaluar
diferentes procesos de la
superficie de la tierra. Es un
software libre a partir de su
versión 3.4.
Desarrollado por ESRI, este
programa permite administrar,
organizar, crear y
previsualizar tanto datos
geográficos como
alfanuméricos; trabaja en
Raster Programa no
ArcGis 1999 Windows ambiente 2D y 3D. Se
Vectorial libre
distribuye comercialmente
bajo tres niveles de licencias
que son, en orden creciente de
funcionalidades (y coste):
ArcView, ArcEditor y
ArcInfo.
Programa basado en 350
módulos y herramientas que
GRASS Raster ejecutan tareas concretas,
1982 Linux Libre
GIS Vectorial puede visualizar y manipular
datos vectoriales y también
raster.
Desarrollado por Bentley
Systems, combina potencia de
Micro Raster Programa no
1985 Windows diseño para el Diseño
Station Vectorial libre
Asistidos por computador
(CAD) 2D y 3D profesional.
76
Continuación Tabla 2.13…
Nombre Fecha de Sistemas Modelo de Licencia de Descripción
creación operativos datos programa
usado
Desarrollado por Golden
Software, este programa
Raster Programa no permite digitalizar, organizar,
Surfer 1983 Windows
Vectorial libre visualizar y manipular
información geoespacial en
ambiente 2D y 3D
El papel de los SIG en la gestión del manejo del recurso hídrico subterráneo es como
instrumento, esto se refiere a las herramientas de gestión como la evaluación y la
información. Incluye además las competencias y destrezas de las instituciones (Cerros,
2007).
77
Los SIG en la gestión de recursos hídricos nos permite elaborar mapas de
vulnerabilidad, evaluar riesgo, realizar diagnósticos de determinadas zonas previo a la
formulación de políticas locales tendientes al desarrollo sustentable.
En ILWIS se manejan dos tipos de datos, éstos son (Escobar, García, & Guerrero, 2002):
78
En Figura 2.15 se muestra el ambiente ILWIS, además de mostrar de forma general
todas las herramientas que se pueden usar al conocer ILWIS.
79
En el área llamada “catálogo” que se muestra en la Figura 2.15 hay una serie de iconos,
que representan toda la información que se crea en ILWIS, dicha información es creada
con ayuda de la barra de operaciones, y a continuación en la Tabla 2.14 se presenta cada
uno.
Georeferencia: Define la relación entre la rejilla en un mapa raster y las coordenadas del
mapa
Mapas raster
Mapas de segmentos
Mapas de puntos
80
Todas las herramientas mencionadas anteriormente son vitales en la elaboración de
mapas de vulnerabilidad, pero necesita de lo que llamamos “algebra de mapas”, que son
operaciones de cálculo que podemos realizar mediante un software SIG para obtener un
nuevo mapa, en nuestro caso nos interesa obtener un mapa de vulnerabilidad para el
método GOD y DRASTIC. Al realizar operaciones de algebra de mapas en ILWIS, es
importante saber que una expresión usualmente contiene operadores y funciones para
especificar el cálculo a realizarse. Los nombres de mapas y las constantes que se usan en
una fórmula son llamadas operandos. Cuando la expresión es ejecutada, el programa
realiza el cálculo en cada píxel, empezando por el primer píxel en la primera fila, y
continuando hasta el último píxel de la última línea (ITC, 2001).
En la Tabla 2.15 se explican las funciones básicas que son necesarias en la elaboración
de mapas de vulnerabilidad.
Una de las ventajas de ILWIS es que considera que no todos los datos que se manejan
son estáticos, ya que los parámetros pueden sufrir cambios intrínsecos como es el caso
de la precipitación media y la recarga neta, también puede haber por adición de nueva
información, o mejoras en la información original. La modificación de cualquier dato
perteneciente a una tabla de ser ajustada con el resto de la información a partir de la
ejecución de la opción “Update all columns” ubicada en el menú principal de la tabla.
81
Una vez modificada la información de las tablas es necesario actualizar los diferentes
mapas relacionados a estas. Para ello se debe desplegar el cuadro de propiedades de tales
mapas y ejecutar la opción de actualización (make up-to-date) ubicada en las
propiedades de dependencia del mapa (Dependency) en la cual son listados todos los
mapas relacionadas a éste.
82
Figura 2.16 Ejemplo de un mapa de vulnerabilidad a la contaminación de un acuífero, resultado de
la aplicación del método GOD y uso de un SIG.
Fuente: Gloria Cerros. Un enfoque de gestión ambiental del recurso hídrico desde la perspectiva de la
vulnerabilidad de acuíferos. Trabajo de Graduación UCA. El Salvador. 2007.
83
CAPITULO 3. Aplicación de los métodos GOD y DRASTIC
para la evaluación de la vulnerabilidad de las
aguas subterráneas en la zona de estudio
Tabla 3.1 Coordenadas y elevaciones de los puntos del perímetro de la zona de estudio
Punto Latitud Longitud Elevación
M‐1 297613.00 477431.00 463.00
M‐2 297643.00 477406.00 470.00
M‐3 297670.00 477392.00 470.00
M‐4 297767.00 477281.00 478.00
M‐5 297753.00 477145.00 487.00
M‐6 297738.00 476810.00 497.00
M‐7 297717.00 476768.00 498.00
M‐8 297720.00 476668.00 496.00
M‐9 297714.00 476622.00 480.00
M‐10 297676.00 476632.00 480.00
M‐11 297645.00 476647.00 480.00
M‐12 297613.00 476670.00 480.00
M‐13 297595.00 476819.00 472.00
M‐14 297541.00 477012.00 468.00
M‐15 297590.00 477212.00 465.00
M‐16 297601.00 477320.00 463.00
M‐17 297594.00 477376.00 473.00
M-4
M-5
M-6
M-8 M-7
M-9
M-10 M-3
Antiguo Botadero de Mariona M-2
M-11
M-1
M-12 M-16
M-13 M-15 M-17
M-14
85
Figura 3.3 Ubicación de la zona de estudio entre los municipios de Nejapa y Apopa
Se desconoce las razones para la selección de este sitio para que funcionara el botadero;
no hubo un análisis técnico, ni mucho menos una Evaluación del Impacto Ambiental
(EIA), en el cual se debía conocer el grado de vulnerabilidad de las aguas subterráneas a
la contaminación para determinar si el lugar era el adecuado o no. La decisión fue
tomada por la administración del entonces Alcalde de San Salvador, Dr. Armando
Calderón Sol.
87
¾ Fuentes de Contaminación
Otras fuentes de contaminación que se ubican en la periferia del antiguo botadero son el
establecimiento de compañías industriales tales como: Coca Cola, Jugos Jumex, Ingenio
88
El Ángel y la construcción de proyectos habitacionales que gastan y/o contaminan los
recursos hídricos de la zona (Rivera y Aparicio, 2009).
Fue hasta el año 2001 que la administración de la Alcaldía de San Salvador comienza
con el cierre técnico del antiguo botadero de Mariona conjuntamente con la OPAMSS y
la empresa MIDES, y comenzar así el uso del relleno sanitario administrado por la
empresa MIDES en el municipio de Nejapa (Rivera y Aparicio, 2009).
Las actividades que se realizaron en el cierre técnico del botadero fueron las siguientes
(ver Figura 3.4):
• Recubrimiento de desechos.
• Siembra de grama.
• Ventilación de biogás.
89
90
Figura 3.4 Obras que se realizaron en cierre técnico del botadero de Mariona
Fuente: MIDES, presentaciones de “Cierre técnico de botadero de Mariona”, s.f.
De las visitas realizadas a la zona de estudio se observó durante los recorrido, que las
mayorías de medidas tomadas en el cierre técnico han desaparecido, como el cerco
perimetral, la ventilación del biogás, se encuentra la presencia de basura en una parte del
antiguo botadero; específicamente cerca de los asentamientos humanos que se ubican a
las orillas del mismo, presencia de lloraderos; asimismo se vio flujo de agua
contaminada en la misma zona, posiblemente de aguas residuales producidas por la
misma gente, además de flujo de lixiviados. El terreno es ocupado actualmente para
pastar ganado, lo cual no es recomendable para la salud de las personas que consumen
alimentos derivados de animales que se alimentan en estas condiciones.
91 ubica este antiguo botadero. En la Figura 3.5 se muestra una secuencia de fotografías de
la evolución que se ha dado en la zona de estudio.
91
A continuación se muestra información de diversas características de la zona, como: La
geología, hidrogeología, relieve, hidrología, climatología, vegetación e información del
perfil geoeléctrico resultado del sondeo realizado por Rivera y Aparicio (2009).
¾ Geología
92
efusivas ácidas (s3’b). Es en parte contemporáneo con los miembros s2, s3’b, s5’a, s5’c,
nombrado informalmente tierra blanca por Williams & Meyer-Abich (1953: 5-6), Weyl
(1955:81) y no incluye las tobas ácidas subyacentes (miembro c1), equivalente a la Serie
Superior de los Estratos de San Salvador (Dürr & Klinge; 1960a, 1960b).
Stirton & Gealey (1949: 1737) mencionan la secuencia en las Llanuras Volcánicas del
Área de Valle Central Interior y la sección esta constituida por tobas pomezosas y
brechas. Tanto los autores como Williams & Meyer Abich (1953: 5-6, 1955: 33) asumen
que este material fue producto de la actividad del Volcán de San Salvador, guardando
estrecha relación con los eventos volcano-tectónicos del lago de Ilopango.
Posteriormente Meyer-Abich (1956: 76) reubica el evento en la depresión de Ilopango,
durante el tercer y último hundimiento.
93
El espesor es superior a los 50 m en las cercanías del Lago de Ilopango (Schmidt-
Thomé, 1975: 216) y se reduce a 4 m en la Ciudad de San Salvador; los fragmentos de
pómez alcanzan tamaños de 25 a 30 cm (Baxter, 1979: 12), y xenolitos son frecuentes.
94
¾ Hidrogeología
95
Figura 3.7 Ubicación de la zona de estudio en el mapa hidrogeológico de la cuenca del Río
Acelhuate
Fuente: Laura Gil. Mapa y estudio de recarga acuífera del área metropolitana de San Salvador. 2007
¾ Relieve
Apopa y Nejapa son municipios situados a una altitud de 200 a 500 msnm (ver Figura
3.8), y se encuentran entre la Cadena Costera y la Cordillera Fronteriza caracterizadas
por una topografía plana comparada con las fuertes pendientes de los sistemas
orográficos que la limitan. Este ex - botadero se sitúa a una altura promedio de 480
msnm (Rivera y Aparicio, 2009).
96
Figura 3.8 Mapa de elevaciones de los municipios de Apopa – Nejapa
Fuente: Francisco Rivera, Osmany Aparicio. Evaluación geoeléctrica del antiguo botadero de residuos
sólidos urbanos de Mariona. 2009.
¾ Hidrología
97
desembocadura, recibe el aporte de varios efluentes, entre ellos: La quebrada
Montserrat, Arenal de Mejicanos, quebrada Tutunichapa -Tomayate, y los ríos El Ángel,
Ilohuapa, Las Cañas, Guazapa y Tasajera (Rivera y Aparicio, 2009).
¾ Climatología
La zona de estudio presenta una precipitación promedio anual que varía entre 1400 –
2400 mm, un rango de temperatura que varía desde 10° C hasta 28° C. (Rivera y
Aparicio, 2009).
98
¾ Vegetación
El uso de suelo para el cual estaba destinada esta zona era para pastos y granos básicos,
pero además este antiguo botadero está muy cercano a áreas urbanas (ver Figura 3.10).
¾ Perfil Geoeléctrico
99
Figura 3.11 Perfil geoeléctrico de un sector donde se ubica el antiguo botadero de Mariona
Arriba: Muestra la ubicación de la realización de los Sondeos Eléctricos Verticales. Abajo: Perfil
elaborado con la información de los Sondeos Eléctricos Verticales
Fuente: Francisco Rivera, Osmany Aparicio. Evaluación geoeléctrico del antiguo botadero de residuos
sólidos urbanos de Mariona. 2009.
100
• La primera capa, se encuentra en un rango de 30 Ω - m a 90 Ω - m, la cual se
extiende a lo largo de todo el perfil con un espesor de aproximadamente de 25 m.
de tierra blanca (Piroclastitas ácidas y eplicastitas volcánicas subordinadas).
101
3.2. Elaboración del mapa de vulnerabilidad de las aguas
subterráneas con la aplicación del método GOD
Por ser relativamente una zona pequeña concluimos que en toda ella se presenta el
mismo tipo de acuífero.
De acuerdo a la Figura 3.11 que muestra el perfil geoeléctrico de la zona del antiguo
botadero y del mapa geológico del país, se pudo determinar que la litología de la
cobertura de la zona de estudio está compuesta por piroclástitas ácidas y eplicastitas
volcánicas subordinadas, conocidas comúnmente como tierra blanca.
102
que se tomó esta de una tabla adaptada de la litología de las diferentes coberturas de San
Salvador (Tabla 3.2), dando como resultado la ponderación de este parámetro 0.80;
dicha tabla tiene origen en el trabajo de graduación de tesis de postgrado en Geología de
Arévalo Herrera en 2005, el cual se denomina “Evaluación de la Vulnerabilidad
Intrínseca del Acuífero Opico-Quezaltepeque, a través de datos geofísicos, pozos y
análisis Hidrogeoquímicos”.
Esta misma litología de la cobertura se presenta en toda la zona de estudio; por lo cual
toda ella tendrá el mismo índice de evaluación.
Para la determinación de este parámetro conocido por “D” que en inglés es “depth” que
significa la profundidad al nivel freático desde la superficie del terreno, se necesitaba de
información sobre niveles freáticos y la altimetría de la zona, por lo que se contó con la
información de los pozos de alrededor de la zona de estudio, tomada del inventario (ver
103
Tablas 3.3 y 3.4, para visualizar localización de puntos de agua dirigirse a Anexo 2) que
aparece en la tesis “Evaluación geoeléctrica del antiguo botadero de residuos sólidos
urbanos de Mariona”, además del archivo ArcView Shape (.shp) con curvas de nivel a
cada diez metros, dando como resultado el mapa que muestra el nivel freático en cada
sector de este lugar, ponderando para ello según sea su profundidad, tal como aparece en
la Figura 2.12 del capítulo 2.
Al usar el software ILWIS se necesita crear una carpeta que contenga cada uno de los
archivos que se originarán para la creación del mapa de vulnerabilidad, esta carpeta se
creó desde el explorar de Windows y la llamamos “Antiguo Botadero de Mariona”,
luego abrimos ILWIS y como si estuviéramos en el explorador de Windows buscamos el
104
sitio donde creamos esa carpeta para comenzar a crear el mapa de vulnerabilidad.
También creamos otra carpeta, denominada “Información base” para guardar mapas de
formato ArcView Shape y tablas que será necesario importar a ILWIS durante el
proceso de creación de mapas de vulnerabilidad.
3.1. En la ventana principal del programa fuimos a la barra de menús, luego a “File”,
“Import”.
3.2. Se localizó las herramientas para importar tablas en “Table”, y se eligió el
formato “Comma delimeted” (ver Figura 3.12).
3.3. Se buscó el archivo Excel previamente hecho en “input”.
105
3.4. Nos dirigimos a “Options” (esta herramienta aparece cuando ya se ha elegido el
archivo a importar) donde se vió inicialmente la información que contenía
nuestro archivo para tener claro qué tipo de dominio debíamos de asignarle; en
nuestro caso teníamos coordenadas y elevaciones las cuales son asignadas como
dominios “value”; aunque por defecto el software identifica el tipo de dominio
que debe ser la información, dentro de esta ventana se asigno nombres a las
columnas dando doble clic en la columna que se deseaba asignar nombre en
“Column Name”; no habiendo otros cambios se dio clic en finalizar, lo que nos
regresaba a la ventana para importar datos.
3.5. Se asignó el nombre a la tabla, colocando el nombre en “Output”, en nuestro caso
la llamamos “perimetro” luego sólo se dió “OK” para que fuese importada la
tabla a ILWIS.
106
4.1. En la ventana principal del programa se dió clic derecho sobre la tabla que
llamamos “perimetro” y se eligió la opción “Table Operations”, “Table to Point
Map” (ver Figura 3.13).
4.2. Elegimos las coordenadas “X” y “Y” de los puntos que deseábamos ingresar.
4.3. En “Coordinate System” seleccionamos “Unknow”, ya que esta opción acomodó
automáticamente el mapa dependiendo de las coordenadas de la tabla.
4.4. Dejamos activada la opción “Use Table Domain” porque no fue necesario usar
las elevaciones obtenidas por medio del GPS.
4.5. Por último en la barra de “Output Point Map” dimos el nombre de nuestro mapa
de puntos en este caso lo llamamos “perimetro”.
4.6. Procedimos a presionar “Show”, para visualizar el mapa que estábamos creando
(ver Figura 3.14), en caso de no querer verlo sólo se presionaba “Define”.
Figura 3.13 Ventana de configuración para creación de un mapa de puntos a partir de una tabla en
ILWIS 3.7
107
Figura 3.14 Mapa de puntos del perímetro de antiguo botadero de Mariona
108
Figura 3.15 Ventana para creación de mapa de polígonos a partir de mapa de puntos en ILWIS 3.7
6.1. En la ventana principal del programa dimos clic derecho al mapa de polígonos
creado llamado “perimetro”, elegimos la opción “Polygon to Raster” (ver Figura
3.17).
6.2. En “Output Raster Map” llamamos al nuevo mapa “perimetro”.
6.3. Creamos una georeferencia, esto puede hacerse antes de entrar a esta ventana o
dentro de la misma, seleccionando el icono para creación de georeferencia que
indica la Figura 3.17, esto abrirá la ventana que se muestra en Figura 3.18
respecto a la creación de referencias geográficas para formato raster, este archivo
109
nos será de mucha importancia para la creación de posteriores mapas; en
“GeoReference Name” llamamos nuestra georeferencia “perimetro”; en
“Coordinate System” seleccionamos “Unknow” y dejamos activado “Center of
Corner pixels”, esto nos sirvió para que las coordenadas de las esquinas de
nuestro mapa vectorial se colocaran en el centro del píxel de las esquinas del
mapa raster, luego solo dimos “OK” para regresar a la ventana de la Figura 3.17.
6.4. Para finalizar solo elegimos en “GeoReference” la georeferencia creada y dimos
un clic en “Show”.
Figura 3.17 Ventana para conversión de formato vectorial a raster en ILWIS 3.7
110
7. Creación de mapa de segmentos de perímetro, éste lo obtuvimos de la siguiente
manera:
7.1. En la ventana principal del programa ILWIS dimos clic derecho en el mapa de
polígonos “perimetro” y nos dirigimos a “Vectorize”, “Polygon to Segment”.
7.2. Encontramos dos opciones que dejaremos por defecto; la opción “Mask”, que
permite buscar polígonos con un nombre específico. Si se tiene activado “*”
buscará todos los polígonos que existan, y la opción “Single Name” es usada
cuando queremos imponer un nombre único a todos los polígonos dentro de un
mapa de polígonos seleccionado (ver Figura 3.19).
7.3. Asignamos el nombre al mapa de segmentos en “Output Segment Map” en
nuestro caso lo llamamos “perimetro”.
7.4. Clic en “Show” para ver el mapa de segmentos.
Figura 3.19 Ventana para conversión de mapa de polígonos a mapa de segmentos en ILWIS 3.7
111
A continuación se describe el proceso realizado para crear los mapas de cada uno de los
parámetros del método GOD, cuando así se requiera, para generar mediante la
multiplicación de éstos, el mapa que nos muestre la vulnerabilidad de las aguas
subterráneas existente en la zona del antiguo botadero de Mariona.
Para este parámetro primero explicaremos como se generó el mapa de nivel freático,
luego el MED y por último el mapa de profundidad de las aguas subterráneas.
112
2. Importación de tabla de datos de pozos y manantiales; se siguió el proceso
mencionado en la delimitación de perímetro numeral 3, solo que se asignó el nombre
de “NF” a la columna que tiene el dato de nivel freático.
3.1. Clic derecho a la tabla de “manantiales” (también pudo haber sido la tabla
“pozos”) dentro de la ventana principal y nos dirigimos a “Table Operation”,
“Glue Tables” (ver Figura 3.20).
3.2. En “Number of Input Tables” seleccionamos “2”.
3.3. En “2nd Table” elegimos la tabla “pozos”.
3.4. En “Output Table” escribimos el nombre de la nueva tabla, la cual llamamos en
nuestro caso “manatiales y pozos”.
3.5. Activamos la casilla “Vertical”.
3.6. Para finalizar clic en “Show”.
113
5. Interpolación de mapa puntos de nivel freático; para este procedimiento hicimos uso
del método “Kriging”, ya que es de los métodos más utilizados en interpolación y
creación de curvas de contorno; y seguimos los siguientes pasos:
5.1. Clic derecho a nuestro mapa de puntos en la ventana principal y nos dirigimos a
“Interpolation”, “Kriging…” (ver Figura 3.21).
5.2. No cambiamos las opciones por defecto, pero sí seleccionamos en
“GeoReference” la georeferencia llamada “perímetro”.
5.3. En “Output Raster Map” llamamos al mapa “nivel_freatico”, para terminar, clic
en “Show”.
Figura 3.21 Ventana para interpolación por método Kringing en ILWIS 3.7
114
6. Intersección de mapa raster de nivel freático con perímetro; teniendo el mapa raster
“nivel_freatico” procedimos a interceptar está información con nuestro perímetro,
para esto seguimos los siguientes pasos:
7. Creación de mapa de atributos; para tener el mapa de nivel freático deseado, debimos
crear un mapa de atributos, ya que al revisar la información observamos que nos daba
115
combinada la información del perímetro y los niveles freáticos, los pasos a seguir
fueron los siguientes:
7.1. Clic derecho dentro de la ventana principal al mapa “cross_nf” luego “Raster
Operations”, “Attribute Map” (ver Figura 3.23).
7.2. Seccionamos en “Table” la tabla “cross_nf”.
7.3. En “Attribute” elegimos la columna “nivel_freatico”.
7.4. En “Output Raster Map” llamamos al mapa “botadero_nf”.
7.5. Por último dimos clic en “Show” y obtuvimos nuestro mapa de nivel freático.
Figura 3.23 Ventana para creación de mapas de atributos raster en ILWIS 3.7
116
9. Intersección de curvas de nivel con perímetro, debido al tiempo que involucraba
trabajar con todo el archivo curvas de nivel a cada 10 metros se decidió hacer un
procedimiento que conllevará menos tiempo de análisis, por lo que el archivo
“curvasdenivel” fue intersectado con el perímetro de nuestra zona en estudio,
haciendo uso del procedimiento siguiente:
11. Interpolación de mapa de segmentos, con el mapa de segmentos ordenado con sus
respectivas elevaciones nos dispusimos a interpolar el mapa generado anteriormente,
con el procedimiento expuesto en el parámetro D de GOD numeral 5, esto creó un
archivo que abarcaba todo lo georefenciado y no sólo el perímetro.
13. Creación de mapa raster de atributos, la columna que usamos para generar este
mapa fue “elevaciones” y lo llamamos “botadero_elevaciones”.
14. Creación de mapa de profundidad de agua subterránea; como lo que nos interesaba
era la profundidad del agua y no su nivel freático, tuvimos que realizar la diferencia
117
entre la topografía de la zona (MED de nuestra zona de estudio) y el nivel freático,
por lo que se aplicó en la barra de comandos de ILWIS la Ecuación 3.1.
15.5. Dimos un clic en “OK” y nos dirigió a la ventana de edición de dominio (ver
Figura 3.25)
118
15.6. Agregamos los intervalos de profundidad en “Edit”, “Add Item”; una vez ahí en
“Upper Bound” seleccionamos el límite superior de nuestro intervalo y en
“Name” el intervalo para identificarlo dentro del mapa.
119
Figura 3.26 Ventana para aplicar dominio a mapa raster en ILWIS 3.7
17. Creación de mapa de puntajes; el último paso para la creación del mapa del
parámetro “D” fue la aplicación de la Ecuación 3.2 en la barra de comandos de
ILWIS.
D = iff (botadero_ D ="0 − 5",0.9, iff (botadero_ D ="5 − 20",0.8, iff (botadero_ D ="20− 50,0.7,0.6)))3.2
Al dar “enter” después de colocar la ecuación, se abrió una ventana donde se mostraba el
tipo de dominio que contendría el mapa y otras características, como no nos interesaba
cambiar los valores por defecto, solo dimos un clic en “OK” para ver nuestro mapa.
En la Figura 3.27 se muestra el mapa del parámetro “D” del método “GOD” con sus
respectivos puntajes.
120
121
121
121
• Aplicación de método GOD
122
123
123
3.3. Elaboración del mapa de vulnerabilidad de las aguas
subterráneas con la aplicación del método DRASTIC
Este parámetro es igual al parámetro “D” del método GOD que ya ha sido mencionado
anteriormente, y el origen de la información es la misma, cambiando solamente el
procesamiento de datos por la diferente manera de evaluar los rangos de profundidades
al nivel freático el método DRASTIC, lo cual se hace en base a la Tabla 2.4.
Esta información fue tomada de un informe elaborado por la OPAMSS en el año 2007
denominado “Mapa y Estudio de Recarga Acuífera del Área Metropolitana de San
Salvador, Zona Prioritaria No 1: Faldas del Volcán de San Salvador – Cuenca del Río
Acelhuate”, la información de recarga neta de la cuenca de Río Acelhuate se muestra en
la Figura 3.29.
Para aplicarlo a nuestro caso primero fue necesario referenciar geográficamente el mapa
de recarga de Rio Acelhuate y digitalizarlo. Se estableció el valor de este parámetro de
acuerdo a la Tabla 2.5, en la cual dependiendo de en qué rango esté el valor de recarga,
así será asignado el puntaje, quedando como se muestra en la Tabla 3.6.
124
Figura 3.29 Mapa de recarga potencial de la cuenca del rio Acelhuate
Fuente: Beatriz Gil. Mapa y estudio de Recarga Acuifera del Area Metropolitana de San Salvador, zona
prioritaria No 1: Faldas del Volcan de San Salvador-Cuenca del Rio Acelhuate. 2007
125
piroclástitas ácidas y eplicastitas volcánicas subordinadas, las cuales se conocen son de
origen volcánico, se le asignó un valor de 9.0 de calificación, tomando en cuenta lo que
dice la Tabla 2.6, siendo el valor típico que recomienda la metodología para este tipo de
medio acuífero.
La información con la que se contó para elaborar el mapa de este parámetro es el archivo
ArcView Shape de curvas de nivel a cada diez metros, el cual consta de un mapa de
126
segmentos y una tabla donde se indica información sobre cada curva, como su longitud y
elevación. Su evaluación se hizo en base a la Tabla 2.8.
3.3.1.6. Impacto del tipo de Zona Vadosa (I = Impact of the Vadose Zone Media)
Este parámetro es similar al parámetro “O” del método GOD, que se refiere al estrato de
la zona vadosa, la cual está formada por piroclastitas ácidas y eplicastitas volcánicas
subordinadas. Esta información se conoció en base al sondeo geoeléctrico realizado en el
trabajo de graduación “Evaluación Geoeléctrica del Antiguo Botadero de Residuos
Sólidos Urbanos de Mariona”, este perfil geoeléctrico se muestra en la Figura 3.11.
De la Tabla 2.9 se estableció una nota de “9.0”, ya que este tipo de litología es volcánica
y está presente en toda la zona de estudio.
127
Tabla 3.8 Determinación de la conductividad hidráulica de la zona de estudio
Datos Fórmula Conductividad (K)
El dato de conductividad hidráulica será el mismo para toda la zona de estudio y que
según la Tabla 2.10 se le asigna un valor de puntuación de 8.
Es de mencionar que estos datos de los cuales se dedujo la conductividad hidráulica son
aproximados, con la salvedad que para asignarle una calificación de 8, el valor de este
podría estar en el rango de 41 a 82 m/día, tal como lo establece la Tabla 2.10.
A continuación se describe el proceso realizado para crear los mapas de cada una de las
variables del método DRASTIC, partiendo de que ya se tienen creado los mapas de
puntos, segmentos, polígono y raster de la zona de estudio, como parte del
procesamiento de la información para el método GOD.
128
Dr = iff (botadero_ Dr"0 − 15",10, iff (botadero_ Dr = 1.5 − 5",9, iff (botadero_ Dr ="5 − 10",7,
iff (botadero_ Dr ="10 − 20",5, iff (botadero_ Dr ="30 − 20",2,1))))) 3.4
En la Figura 3.30 se muestra el mapa del parámetro D del método DRASTIC con sus
respectivos puntajes.
129
130
130
• Parámetro “R” de DRASTIC
Para la elaboración del mapa de este parámetro se aplicaron los siguientes pasos:
Figura 3.31 Ventana para crear georeferencia de una imagen en ILWIS 3.7
131
2.4. En “GeoReference Name” llamamos a la nueva georeferencia “recarga potencial”
2.5. Elegimos la opción “GeoRef Tiefpoints”, ya que consiste en un método de
referenciar mediante la ubicación de puntos dentro de la imagen, asignándoles
posteriormente sus coordenadas; para las otras alternativas fue necesario contar
con el MED (GeoRef Direct Linear), usado para ortofotografías (GeoRef Ortho
Photo) y para fotografías (GeoRef Parallel Projective).
2.6. En Coordinate System se eligió “unkown”, después “OK”, entonces el programa
nos dirigió a la ventana para editar georeferenciaciones como se muestra en la
Figura 3.32.
2.7. En el editor de georeferencias seleccionamos el método “Conformal”, este
permitió referenciar las imágenes de forma debida con solo tres puntos, para otros
métodos se necesita mayor información de la que proporciona el mapa con que
contamos.
2.8. Elegimos el icono que se indica en Figura 3.32 y dimos un clic en uno de los
puntos que conocíamos las coordenadas dentro de la imagen, cuando dimos ese
clic se desplegó una ventana donde apareció la información de la columna y fila
donde se encontraba el punto introducido y al lado colocamos la coordenada X y
Y que correspondía a dicho punto, esto se hizo en los otros dos puntos que se
tenían.
2.9. Para finalizar salimos del editor buscando en el menú “File”, “Exit Editor” y la
imagen quedó referenciada.
132
Figura 3.32 Ventana de editor de georeferenciaciones en ILWIS 3.7
3. Creación de dominio con rangos de recarga neta; se creó un dominio “Class” con los
diferentes intervalos que nos indicaba el mapa de recarga para tenerlo disponible
cuando se fuera digitalizando, a este lo llamamos “recarga neta”.
4. Digitalización del mapa de recarga neta; siendo tantos los polígonos que deberían
digitarse se tomo la decisión de eligir aquellos que intersectaran con el antiguo
botadero, y realizamos el siguiente proceso:
133
4.2. Nos dirigimos a “File”, “Create”, “Segment Map” (ver Figura 3.33)
4.3. En “Map Name” se nombró al mapa de segmentos como “segment_recarga”
4.4. No cambiamos las coordenadas que se indican en la ventana
4.5. En dominio elegimos “recarga_neta” y después “OK”. Se abrió la ventana de
edición de mapas de segmentos en la que visualizamos el mapa.
4.6. Para dibujar los segmentos se usó “Edit”, “User mode”
4.7. Apareció una ventana para elegir el código del segmento que deseabamos digitar,
se pudo elegir de la lista que ingresamos cuando se hizo el dominio. Estos
segmentos se podían cambiar posteriormente, pero se recomienda ir trabajando de
forma ordenada para evitar confusiones posteriores
4.8. Definimos el contorno de los polígonos, clickeando por cada cambio de dirección
para introducir líneas o también de modo continuo manteniendo el clic izquierdo y
desplazándonos para ir dibujando los contornos de los polígonos, cuando ya se
terminó de dibujar el segmento sólo se dió un doble clic. Es recomendable que en
dado caso nos equivocamos a la hora de digitalizar, no se interrumpa el dibujo de
segmentos, pues se debe proseguir y posteriormente se puede modificar el error en
el menú “Edit”, “Move Points”.
134
4.9. Con los segmentos realizados debimos hacer una serie de comprobaciones para
evitar errores a la hora de crear los polígonos, estas comprobaciones se
encontraban dentro del menú “File”, “Check Segment” y se describen en Tabla
3.9; si teníamos algún error el programa nos indicaba donde estaba ese error.
135
4.14. Activamos la casilla “Edit” para verificar que nuestro mapa se generará de forma
correcta, y por último “OK”.
4.15. En la ventana de edición de mapa de polígonos, verificamos que todo estuviera
bien y salimos del editor, “file”, “Exit Editor”.
7. Creación de mapa de puntajes; los últimos pasos a realizar para la obtención del mapa
final fue la aplicación de la Ecuación 3.5 en la barra de comandos de ILWIS:
136
Rr = iff (botadero _ rec arg a ="145 − 340" ,8,9) 3.5
En la Figura 3.35 se muestra el mapa del parámetro R del método DRASTIC con sus
respectivos puntajes.
137
138
138
• Parámetro “A” de DRASTIC
139
Tabla 3.10 Clases de suelos según mapa pedológico de El Salvador
CLASE TIPO DE SUELO
1 ALUVIALES
2 GRUMOSOLES
3 ANDISOLES
4 LATOSOLES ARCILLO ROJIZOS
5 REGOSOLES Y HALOMORFICOS
6 LITOSOLES
7 LATOSOLES ARCILLOSOS
ACIDOS
8 AREA URBANA
4. Intersección de mapa de tipos de suelos con perímetro; se uso la función “Cross” para
realizar tal intersección, a los mapas y tabla que se generaron los llamamos
“cross_suelos” y como ya ha sido explicado anteriormente, creamos un mapa de
atributos para obtener un mapa con los códigos que deseamos; elegimos dentro de la
ventana de mapas de atributos la tabla “cross_suelos” y el atributo “suelos” y
llamamos a este mapa “botadero_suelo”.
En la Figura 3.36 se muestra el mapa del parámetro S del método DRASTIC con sus
respectivos puntajes
140
141
141
• Parámetro “T” de DRASTIC
142
6. Creación de mapa de puntajes; con el mapa de pendientes finalmente terminado
procedimos a asignar los puntajes respectivos a cada intervalo haciendo uso de la
Ecuación 3.10.
Tr = iff (botadero _ pendientes ="0 − 2" ,10, iff (botadero _ pendientes ="2 − 6" ,9,
3.10
iff (botadero _ pendientes ="6 − 12" ,5, iff (botadero _ pendientes ="12 − 18,3,1))))
En la Figura 3.37 se muestra el mapa del parámetro T del método DRASTIC con sus
respectivos puntajes.
No se hará procesamiento para crear el mapa, ya que se tiene un solo tipo de material,
por lo cual este parámetro tendría una sola calificación en toda la zona de estudio.
Este parámetro tuvo un solo valor y por ello una sola calificación en toda la zona de
estudio, por lo que también no se presenta su mapa correspondiente.
143
144
144
• Aplicación de método DRASTIC
El resumen de los valores a sustituir o los mapas a utilizar se presentan en la Tabla 3.11.
DRASTIC = 5 * Dr + 4 * Rr + 3 * 9 + 2 * Sr + Tr + 5 * 9 + 3 * 8 3.11
Así como se realizó antes en GOD, usamos la función “slicing” para una mejor
interpretación de los resultados, el domino a crear se llamó “rangos_DRASTIC” y se
utilizó de base la Tabla 2.5. El resultado de ésto fue el mapa final que se llamó
“Vulnerabilidad_DRASTIC”.
146
CAPITULO 4. Análisis de resultados y presentación de
medidas de seguimiento medioambientales y
técnicas de remediación de aguas
subterráneas
En la Figura 3.28 del capítulo 3 puede observarse el resultado final de la aplicación del
método GOD, el cual es un mapa donde se obtuvieron valores de vulnerabilidad que
varían de 0.34 a 0.38, encontrándose dentro del rango de “0.3 – 0.5” que propone la
metodología, determinándose que el área de estudio tiene vulnerabilidad “moderada”;
misma que se encuentra entre la vulnerabilidad baja y alta; que indica que estas aguas
subterráneas tienen una sensibilidad media ante la contaminación, por lo que los
esfuerzos por remediar el daño hecho al acuífero también serian de un nivel intermedio.
En consecuencia, este grado de vulnerabilidad demuestra que se debe tener cuidado con
las actividades antropogénicas que se deban de realizar en la zona, ya que la ubicación
de estas podrían ser una fuente potencial de contaminación en el área. Para haber
establecido este botadero se debió conocer que vulnerabilidad presentaban las aguas
subterráneas del lugar; aunque para establecer un botadero ningún lugar es adecuado,
mucho menos aquellos que presentan vulnerabilidades que van de moderadas a
extremas. Lamentablemente en este sitio no se tomaron las medidas necesarias para
evitar el flujo de contaminantes en el subsuelo, ya que ha sido comprobada previamente
la existencia de una pluma de contaminación, compuesta por los lixiviados que se
producen a causa de este antiguo botadero dirigida al municipio de Apopa (Rivera y
Aparicio, 2009); esto tampoco quiere decir que sea el único municipio afectado.
Para analizar los resultados obtenidos, debemos comprender que el método GOD está
destinado a aplicarse a grandes regiones, zonas que no poseen mucha información
relacionada a las condiciones de un acuífero, por lo que se espera definiciones menos
claras en cuanto a la evaluación de la vulnerabilidad que al usar métodos como
DRASTIC o SINTACS. Además, al tener condiciones tan homogéneas por la
información obtenida y no ser una extensión muy grande la zona estudiada, era de
esperarse que no hubiera mucha variación entre los valores de vulnerabilidad en toda el
área, siendo el único parámetro que presentaba un poco de variación el de la profundidad
al nivel freático (D), y al ser tan amplios los intervalos para asignar el índice de
evaluación; no se esperaba también demasiada diferencia en los resultados, por lo que se
obtiene un mapa con una sola clasificación de vulnerabilidad.
100
80
Area afectada (%)
60
40
20
0
0)
)
)
)
3)
.7
.5
.1
1.
0.
0
0
0
‐
5‐
3‐
‐
0‐
0.7
.1
0.
0.
0.
(0
(
(
(
(
a
ta
a
aja
ja
ad
Al
Ba
m
er
b
tre
uy
od
Ex
M
Vulnerabilidad
Figura 4.1 Distribución de la vulnerabilidad de las aguas subterráneas a la contaminación de la
zona donde se ubica el antiguo botadero de Mariona. Método GOD.
148
4.2. Método DRASTIC
Con el método DRASTIC se obtuvo valores para la vulnerabilidad que varían entre 144
a 191, por lo que la zona de estudio presenta tres tipos de vulnerabilidad, tal como se
muestra en la Figura 3.38, las cuales son: Moderada, alta y muy alta; en donde la
primera tiene 31.9 por ciento del área total de la zona, la segunda 52.0 por ciento y la
tercera 16.1 por ciento respectivamente, estas distribuciones se muestran en la Tabla 4.1
y Figura 4.2 con sus respectivos rangos de calificación correspondientes y color que los
caracteriza.
Tabla 4.1 Distribución de la vulnerabilidad de las aguas subterráneas de la zona donde se ubica el
antiguo botadero de Mariona. Método DRASTIC
Vulnerabilidad (rango) Área afectada (%)
Despreciable (<100) 0.0
Muy baja (101-119) 0.0
Baja (120-139) 0.0
Moderada (140-159) 31.9
Alta (160-179) 52.0
Muy alta (180-199) 16.1
Extrema (>200) 0.0
149
Para el análisis de los resultados debemos tomar en cuenta que el método posee varias
limitaciones, algunas asociadas a factores externos como la calidad de la información y
otras asociadas a las bases con que fue planteado el método. Entre los factores externos
podemos mencionar que al preparar el mapa de vulnerabilidad de las aguas subterráneas
se tuvo que concluir que teníamos un área de estudio homogénea, esto a partir de un
perfil geoeléctrico, ya que no se contaba con columnas litológicas que demostraran lo
contrario, se consideró que por ser una zona relativamente pequeña, no variaría mucho
su litología; además se tuvo que asignar un único valor de conductividad hidráulica
porque no se logró obtener un inventario de pozos con este tipo de información, ni
tampoco existían estudios que relacionaran la conductividad o la transmisividad del
acuífero con datos obtenidos de sondeos eléctricos, esto ya se ha hecho en estudios de
algunas cuencas de Latinoamérica. También se tuvo que trabajar con curvas de nivel a
cada 10 metros, ya que se quería la información topográfica de la zona antes de que se
ubicara el botadero en la zona.
Entre las limitaciones que tiene el método por como fue formulado podemos mencionar
el hecho de que algunos autores critican la reiteración en el alcance de algunos
parámetros como “R” y “C” ya que ambos están vinculados a la renovación del agua en
el acuífero y la poca incidencia que tiene el parámetro “S”.
Los resultados obtenidos indican que es necesaria una mayor protección en el lado sur-
oeste y este del área estudiada, ya que posee vulnerabilidad de alta a muy alta, y de las
visitas realizadas se verificó la presencia de basura a nivel superficial en el lado sur-
oeste.
En el lado norte del antiguo botadero es menor la vulnerabilidad comparada con el lado
sur, ya que tenemos presencia de vulnerabilidad moderada a alta. No existe presencia
significativa de basura, pero se observan lloraderos superficiales en este sector.
150
En general más de la mitad del área de la zona de estudio se encuentra en un rango alto
de vulnerabilidad, lo que es un indicativo de que en la zona no se debió ubicar un
botadero a cielo abierto y que lamentablemente no se tomaron las medidas adecuadas
para evitar la contaminación de las aguas subterráneas de este sector.
Los resultados obtenidos por los métodos GOD y DRASTIC coinciden únicamente en
un 31.9 por ciento, el cual es el porcentaje de vulnerabilidad “moderada” que presenta la
zona aplicando el método DRASTIC, mientas que con la aplicación de GOD se obtuvo
que todo la zona presenta este mismo tipo de vulnerabilidad.
La diferencia entre los métodos GOD y DRASTIC son hasta cierta manera obvias,
primeramente el número de parámetros que considera cada método para evaluar la
vulnerabilidad de las aguas subterráneas de una determinada zona son diferentes, el
método GOD considera tres tipos de parámetros; mientras el DRASTIC considera siete.
Dos parámetros se utilizan en común por los métodos, que son: la profundidad al nivel
freático (“D” en ambos métodos) y la cobertura del acuífero (“O” para GOD e “I” para
DRASTIC), el método DRASTIC no toma en cuenta el grado de confinamiento del
acuífero(G); pero sí toma en cuenta la recarga neta (R), el material que constituye el
acuífero (A), tipo de suelo (S), la topografía (T) y la conductividad hidráulica (C); en
donde el otro método no los considera, por lo que podemos decir que DRASTIC es más
representativo y que en lo económico requiere de un presupuesto más alto que GOD,
esto debido a la cantidad de información que se debe de utilizar para aplicarlo. Por lo
cual su uso se ve limitado muchas veces.
151
el método GOD y Tabla 2.4 para el método DRASTIC, ambas en el capítulo 2, en donde
se puede observar que DRASTIC asigna una mayor calificación que GOD dependiendo
del rango evaluado, y es por ello también que se tiene menos diferencia en el resultado
obtenido por el método GOD. En el caso del parámetro “I” de DRASTIC, la calificación
obtenida de la Tabla 2.9 del capítulo 2 es mayor a la propuesta en la Tabla 3.2 del
capítulo 3. Además los parámetros “D” e “I” son los de mayor trascendencia para
DRASTIC, ya que son a los que se les asigna un mayor peso o factor de importancia
(ver Tabla 3.11), por lo que se puede decir; que es una de las causas por las que se tienen
valores más altos de vulnerabilidad en DRASTIC que en GOD.
Cabe mencionar que para la realización de este estudio se trabajó en su mayoría con
información de estudios realizados anteriormente y que no se contó con el inventario de
pozos de la zona y de sus alrededores que incluyera información de la litología del lugar,
por lo que se usó un perfil geoeléctrico de la zona. Además se trabajó con las
condiciones originales del lugar, con el objetivo de conocer el grado de sensibilidad de
las aguas subterráneas de esta zona a la contaminación; tal como se muestra en los
mapas de vulnerabilidad, resultados de la aplicación del método GOD y DRASTIC.
152
Dentro de las fallas que se pudo observar están las que se presentan en las siguientes
Tablas, desde la 4.2 a la 4.9; en donde se detalla su ubicación en la zona, a qué afecta, la
posible solución y recomendaciones de ser necesario; también se muestran fotografías
que ilustran su estado, las cuales se tomaron en el primer semestre de este año 2010.
Las soluciones en relación a las fallas detectadas que se plantean son relativamente
sencillas, no requieren de tecnologías caras y económicamente pudiesen ser factibles
para quien debería ser el responsable de las medidas de seguimiento medioambientales
en donde funcionó este botadero a cielo abierto.
153
Tabla 4.2 Presentación de medida de seguimiento medioambiental 1
Foto(s):
N° 1/8
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
Posible solución Reconstrucción
inmediata de muro
perimetral, que impida
el acceso a esta zona.
Recomendaciones Se debe vigilar regularmente, ya que el cerco anterior fue robado
por partes y tampoco se le dio el mantenimiento adecuado.
154
Tabla 4.3 Presentación de medida de seguimiento medioambiental 2
Foto(s):
N° 2/8
Falla detectada Presencia de lloraderos
superficiales y cercanos
a la población.
Ubicación Sector norte de antiguo
botadero.
Fecha de visita 22 de febrero de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
155
Tabla 4.4 Presentación de medida de seguimiento medioambiental 3
Foto(s):
N° 3/8
Falla detectada Daños en el sistema de
evacuación de gases
(Chimeneas).
Ubicación Distribuidas en toda la
zona
Fecha de visita 22 de febrero de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
156
Tabla 4.5 Presentación de medida de seguimiento medioambiental 4
Foto(s):
N° 4/8
Falla detectada Incursión de ganado.
Ubicación Toda la zona de este
antiguo botadero.
Fecha de visita 22 de febrero de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
Posible solución Restringir completamente
el ingreso de animales y
personas al lugar; para
eso ya se planteó la
construcción de un muro
perimetral.
Recomendaciones Los productos derivados de estos animales alimentados en estas
condiciones son dañinos para la salud humana, por ello se debe
inventariar el ganado y hacer estudios de los derivados de estos.
157
Tabla 4.6 Presentación de medida de seguimiento medioambiental 5
Foto(s):
N° 5/8
Falla detectada Obstrucción de canaletas de
desalojo de aguas lluvias y
lixiviados.
Ubicación Sector sureste de botadero
Fecha de visita 28 de Mayo de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros
158
Tabla 4.7 Presentación de medida de seguimiento medioambiental 6
Foto(s):
N° 6/8
Falla detectada Perdida de material de
cobertura por procesos
de erosión.
Ubicación Varios sectores.
Fecha de visita 28 de Mayo de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros
159
Tabla 4.8 Presentación de medida de seguimiento medioambiental 7
Foto(s):
N° 7/8
Falla detectada Promontorios de basura.
Ubicación Principalmente en
sectores de zona sur‐
oeste.
Fecha de visita 28 de Mayo de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
160
Tabla 4.9 Presentación de medida de seguimiento medioambiental 8
Foto(s):
N° 8/8
Falla detectada Pozo de visita sin reja.
Ubicación Sector Noroeste de
botadero.
Fecha de visita 22 de febrero de 2010
Afecta a:
Pobladores de la zona 9
Aguas Subterráneas 9
Otros 9
Recomendaciones Fijar la reja de una forma segura para que no sea quitada de este
lugar.
161
Respecto a las técnicas que se mencionan a continuación; éstas van enfocadas a la
remediación de aguas subterráneas debido a la contaminación generada por este antiguo
botadero, pero su uso depende en gran medida del factor económico y de la tecnología
con que se disponga, ya que la mayoría de ellas tienen un alto costo económico y
requiere de diversa información, tal como se detallará en cada una de ellas.
La limpieza del agua subterránea puede hacerse después de que el agua ha sido extraída
(tratamiento ex situ) o en el acuífero (in situ). En esta sección discutiremos primero los
problemas planteados al usar el método tradicional ex situ de bombeo y tratamiento.
Finalmente, métodos in situ que, cuando funcionan, disminuyen enormemente los costes
de recuperación de las aguas.
162
observar una relación hipotética de estos datos para un sistema convencional. El ejemplo
corresponde a una pluma de 190 millones de litros con una concentración de 1 mg/l de
TCE (Tricloroetileno C2HCl3), usando los costes estimados por la EPA y considerando
una tasa de descuento del 4%. Según N. R. Council (1994) para eliminar el 80% del
contaminante se calcula que hacen falta unos 15 años y unos 2.8 millones de dólares
(incluyendo capital, operación y costes de mantenimiento). Para limpiar el 99.99%, el
sistema costaría 6 millones de dólares y se tardarían 84 años. Obviamente hay muchas
suposiciones e incertidumbres en un cálculo de este tipo, pero las implicaciones están
claras.
Figura 4.3 Impacto de la expectativa de limpieza en los costes y duración de un sistema convencional
de bombeo y tratamiento.
Los puntos superiores corresponden a una eliminación del TCE del 99 y del 99,99%. Datos para un
volumen de pluma de 190 millones de litros, una tasa de bombeo de 380 l/min, un factor de retardo de 4,8
(TCE), y una tasa de descuento del 4%.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 271
Una hipótesis implícita en la Figura 4.3 es que la zona puede ser recuperada usando
técnicas de bombeo y tratamiento. Un estudio reciente del Comité de Alternativas de
Limpieza de Aguas Subterráneas (Committee on Ground Water Cleanup Alternatives)
trabajando para el National Research Council, ha evaluado 77 lugares contaminados con
sistemas de bombeo y tratamiento en operación (N. R. Council, 1994) y ha encontrado
sólo unas pocas circunstancias en las cuales estos sistemas conseguirán probablemente
163
su objetivo en un lapso y a un coste razonables. Ha establecido una clasificación de
lugares basada en la química del contaminante y su hidrogeología que intenta ordenarlos
por su facilidad relativa de limpieza. Se han encuadrado en cuatro categorías. Los
lugares de categoría 1 son los más fáciles de limpiar y los de categoría 4 los más
difíciles. Los de categoría 1 se caracterizan por ser acuíferos homogéneos, de una sola
capa, contaminados por sustancias que se mueven y disuelven fácilmente y que se
eliminan hasta cierto punto de una forma natural por degradación o volatilización. Las
técnicas de bombeo y tratamiento deberían ser capaces de restaurar la calidad de agua
potable en estos pocos lugares. En el otro extremo, los puntos de categoría 4 cuentan con
una compleja hidrogeología de materiales heterogéneos y múltiples capas, con
intrincados recorridos de los contaminantes muy difíciles de predecir. Los contaminantes
de la categoría 4 son NAPL (líquidos de fase no acuosa), con DNAPL (NAPL densos) y
LNAPL que son NAPL ligeros (ver Figura 4.4 y Anexo 3).
Figura 4.4 Los líquidos de fase no acuosa no se disuelven bien en el agua subterránea.
(a) Los LNAPL flotan sobre ella; (b) Los DNAPL forman piscinas.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 269
Master y Ela (2008) señalan que en los lugares de categoría 4 no es probable alcanzar el
objetivo de calidad de agua potable. De los 77 sitios estudiados por el Comité, sólo dos
fueron clasificados como categoría 1, mientras que en la categoría 4 se situaron 42. La
siguiente Tabla 4.10, muestra las características de cada grupo en esas cuatro categorías.
164
Tabla 4.10 Facilidad relativa de limpieza de acuíferos contaminados en función de la química del
contaminante y la hidrogeología del lugar
Química del contaminante
Altamente
Móvil, se
Hidrogeología absorbido, Altamente Fase Fase
disuelve (se Móvil, se
se disuelve absorbido, separada separada
degrada/ disuelve
(se degrada se disuelve LNAPL DNAPL
volatiliza)
volatiliza)
Homogénea,
1-2 2 2-3 2-3 1
capa simple
Homogénea,
1* 1-2 2 2-3 2-3 3
múltiples capas
Heterogénea,
2* 2 3 3 4
capa simple
Heterogénea,
2 2 3 3 3 4
múltiples capas
Fracturada 3* 3 3 4* 4
(*) Facilidad relativa de limpieza, donde l es la más fácil y 4 la más difícil.
Fuente: National Research Council. 1994.
Para nuestro caso la hidrogeología que presenta la zona es homogénea, con una sola
capa y según la Tabla 4.10 tendría una facilidad de limpieza que va de 1 a 3,
dependiendo de la química del contaminante, de los cuales no se tiene información
actualizada, por lo que no podríamos asignar un valor específico de facilidad de
limpieza.
165
La Figura 4.5 muestra el efecto de un pozo de extracción situado en una zona con un
flujo estable y uniforme de agua subterránea, paralelo al eje x y en la dirección negativa
de éste. Cuando se extrae el agua las líneas de corriente se curvan hacia el pozo como se
observa en la Figura 4.5. La envolvente de estas líneas de corriente que convergen en el
pozo se denomina curva de zona de captura. El agua comprendida dentro de esta zona
será extraída; la que cae fuera, no lo será. Las líneas de flujo exteriores a la zona de
captura pueden curvarse hacia el pozo pero el flujo asociado al gradiente hidráulico
natural es suficientemente fuerte como para que el agua pase de largo por él.
Figura 4.5 Un único pozo de extracción situado en x=0, y=0, en un acuífero con flujo a lo largo del
eje x.
La zona de captura es la región en la cual todas las líneas de flujo convergen en el pozo de extracción.
Dibujado para Q/Bn — 2000.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 260
Javandel y Tsang (1986) desarrollaron la técnica de las curvas de zona de captura como
ayuda para diseñar campos de pozos de extracción para limpieza de acuíferos. Su
análisis se basa en un acuífero ideal (es decir, homogéneo, isótropo, uniforme en sección
e infinito en anchura). También, sea libre o cautivo, con unos descensos de nivel en
pozos, insignificantes en comparación con el espesor total del acuífero. Suponen que los
pozos de extracción se adentran en el total del acuífero y extraen uniformemente de él a
cualquier nivel. Estas son hipótesis muy restrictivas que es improbable que sean
166
satisfechas en condiciones reales; a pesar de ello los análisis resultantes proporcionan
una buena perspectiva de los factores principales que se obtienen con modelos más
realistas pero más complejos.
Q Q y (4.1)
y=± − tan −1
2 Bv 2πBv x
Donde:
B = Espesor del acuífero (m).
v = Velocidad de Darcy (m/día).
Q = Tasa de bombeo del pozo (m3/día).
La Ecuación (4.1) puede reescribirse en términos del ángulo (en radianes) del radio
vector al punto de la zona de captura considerado, como se muestra en la Figura 4.6.
Esto es:
y (4.2)
tan φ =
x
167
Figura 4.6 Curva de zona de captura para un pozo único de extracción situado en el origen, en un
acuífero con velocidad de flujo v, espesor B y tasa de bombeo Q.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 244
168
Figura 4.7 Superposición de la pluma en una curva de zona de captura para un solo pozo de
extracción
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 262
169
Así la pluma de 80 m de anchura cabrá perfectamente en la zona de captura si el pozo
está situado a cierta distancia gradiente abajo respecto del borde frontal. Usando la
Figura 4.8 como guía, podemos determinar la distancia x que debe separar la pluma del
pozo. De (4.3) con y = 40 m:
Q ⎛ φ⎞ ⎛ φ⎞
y= ⎜1 − ⎟ = 40 = 50⎜1 − ⎟
2Bv ⎝ π ⎠ ⎝ π⎠
Así, el ángulo en radianes en el punto donde la pluma toca la zona de captura es:
φ = 0.2π rad
Y, de la Figura 4.8:
y 40
x= = = 55m
tan φ tan (0.2π )
El pozo de extracción debería perforarse 55 m por delante de la pluma y directamente
en línea con ella.
Figura 4.8 Ubicación de un único pozo de extracción por delante de una pluma de contaminante
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 263
170
Una mejor solución implica más pozos de extracción colocados próximos a la cabeza de
la pluma. Javandel y Tsang han deducido curvas de zona de captura para un conjunto de
n pozos óptimamente emplazados bombeando al mismo ritmo Q, y alineados a lo largo
del eje y. Se define optimización como el máximo espaciado entre pozos que no
permiten el paso de flujo a través de la línea que los une. La distancia de separación
entre dos pozos se ha calculado en Q/πBv. Si los pozos están más separados que esto,
algo de flujo puede pasar entre ellos y no ser capturado. Con este espaciado óptimo, los
dos pozos capturarán una pluma de anchura Q/Bv medida sobre el eje y, y de anchura 2
Q/Bv suficientemente lejos de los pozos, como se muestra en la Figura 4.9(a).
Parámetros análogos para el caso de tres pozos óptimamente espaciados se dan en la
Figura 4.9 (b).
Figura 4.9 Curvas de captura para pozos óptimamente espaciados a lo largo del eje y, bombeando a
una tasa Q: (a) dos pozos; (b) tres pozos.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 264
Una ecuación general para la mitad positiva de la curva de captura para n pozos
óptimamente espaciados, situados regularmente a lo largo del eje y, es:
Q ⎛ 1 ⎞ (4.4)
∑ φ⎟
n
y= ⎜n −
2 Bv ⎝ π i =1 i
⎠
171
Donde es el ángulo que forma la recta que une el pozo i-ésimo con un punto de la
curva de captura, con la horizontal.
Para demostrar el uso de estas curvas, repetiremos el Ejemplo 4.2, pero esta vez con dos
pozos.
Solución.
a) La anchura de la pluma sobre el eje y (y el borde de ataque de ésta) es de 80 m,
así a la vista de la Figura 4.9.
Q Q
= = 80m
Bv 20mx2 x10 −6 m
s
3
Q = 0.0032 m (cada uno)
s
De la Figura 4.9, el espaciado óptimo entre dos pozos viene dado por
172
3
Q 0.0032 m
Espaciado optimo = = s = 25.5m
πBv πx 20mx2 x10 −6 m
s
Figura 4.10 Ubicación de dos pozos de extracción en el frente de una pluma de contaminante
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 265
480.000m 3
t= = 2.4 años
m2 s h dias
0.0064 x3600 x 24 x365
s h dia año
173
Las plumas de contaminación reales y los acuíferos son mucho más complejos y diseñar
un campo de pozos es una tarea mucho más difícil que la que se ha presentado aquí, por
lo que se requiere de una amplia información para su aplicación.
Los pozos de alrededor del antiguo botadero están siendo contaminados por los
lixiviados que éste genera, pero existe un modo de proteger a estos y es combinando
pozos de extracción con pozos de inyección. Los pozos de extracción se usan para
rebajar el nivel freático (o superficie piezométrica), creando un gradiente hidráulico que
arrastra la pluma a los pozos. Los pozos de inyección elevan el nivel freático y empujan
la pluma. Por medio de un cuidadoso diseño del emplazamiento y ritmo de bombeo de
estos pozos, podemos manipular el gradiente hidráulico de tal modo que las plumas
puedan mantenerse apartadas de los pozos de agua potable y dirigirlas hacia los de
extracción. Una vez extraída, el agua contaminada puede ser tratada y reinyectada de
nuevo al acuífero, utilizada o liberada en el sistema hídrico superficial.
La Figura 4.11 muestra dos estrategias para controlar una pluma contaminante que
amenaza un pozo cercano, aquí llamado pozo de producción. En la Figura 4.11 (a) un
pozo de inyección empuja la pluma lejos del de producción. Pozos adicionales serían
perforados sobre la pluma misma para extraerla. En la Figura 4.11 (b), se usa diferente
método. El pozo de inyección se usa para empujar la pluma hacia el de extracción y así
ésta puede ser evacuada más rápidamente del acuífero.
174
Figura 4.11 Manipulación del gradiente hidráulico con múltiples pozos.
(a) Los pozos de inyección y extracción empujan y arrastran la pluma lejos del pozo de producción; (b) El
pozo de inyección empuja la pluma hacia el de extracción.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 266
175
impermeable de tierra y bentonita. La barrera de lodo mantiene la pluma confinada
mientras se aplican otras medidas de recuperación para limpiar el acuífero.
Intentar conseguir la limpieza del agua subterránea hasta el estándar de agua potable
parece un objetivo vano para los acuíferos más contaminados. Esta eficacia limitada ha
motivado el desarrollo de técnicas alternativas para potenciar o sustituir los sistemas
convencionales. Siguen ahora unas descripciones generales de las técnicas más
ampliamente usadas.
• Biorrecuperación in situ
Las técnicas de biorrecuperación de acuíferos pueden ser una alternativa viable a aplicar
en un sitio que este generando lixiviados, porque existen numerosos microorganismos
que biotransforman sustancias inorgánicas, como los metales pesados; es por ello que se
presenta esta técnica como una opción para remediar las aguas subterráneas
contaminadas de este sector del antiguo botadero y de sus alrededores. Estos sistemas in
situ se han diseñado para degradar la polución del subsuelo sin necesidad de extraer los
contaminantes y someterlos a un sistema de tratamiento en superficie. Tratando los
contaminantes in situ se elimina la necesidad de excavación, tratamiento en superficie y
transporte a un vertedero, y el riesgo de exposición de los operarios a los productos
químicos peligrosos se reduce en gran medida. Hay otras ventajas también, entre ellas el
potencial para degradar compuestos que son sorbidos por los materiales del subsuelo.
176
ser tan simple como un túnel de infiltración que permita percolar agua rica en oxígeno y
nutrientes hasta la zona de contaminación. Cuando se encuentra más profunda, es
necesario un sistema de inyección como el que se muestra en la Figura 4.12 que haga
llegar la preparación acuosa a la zona deseada. Se puede suministrar el oxígeno
necesario para la descomposición aerobia haciendo burbujear aire u oxígeno en el agua
inyectada, o añadiendo peróxido de hidrógeno a ésta. Los nutrientes añadidos suelen ser
nitrógeno y fósforo cuando los contaminantes son hidrocarburos. Se pueden añadir
pozos de extracción para controlar el flujo del agua inyectada.
177
Figura 4.12 Técnica de remediación de aguas subterráneas por Biorrecuperación in situ.
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 275
Los disolventes clorados pueden también degradarse de forma anaerobia. En este caso
las bacterias biodegradantes son otras y entre los nutrientes requeridos se cuentan
donadores de electrones, como el metanol, la glucosa o el acetato. Algunos productos
intermedios de esta biodegradación, como el cloruro de vinilo, son peligrosos en sí
mismos, por lo que hay que actuar con precaución. Es posible llegar con esta
degradación a etano y eteno simples, pero los procesos anaerobios transcurren mucho
más lentamente que los correspondientes aerobios, por lo que esta degradación puede
llevar años.
Las barreras reactivas permeables (PRB), son estructuras subterráneas que sirven para
tratar aguas subterráneas contaminadas en situ, producidas por el avance de la pluma de
contaminación generada por los lixiviados de los desechos peligrosos enterrados en
zonas como en nuestro caso, tratándose de un antiguo botadero y conociendo que esta
técnica es aplicada en sitios contaminados con diversas sustancias, entre las que se
cuentan: metales pesados (cromo, níquel, plomo, cobre, manganeso, uranio, arsénico,
178
selenio, vanadio, etc), Compuestos Orgánicos Volátiles (que incluyen solventes clorados
como el Tetracloroeteno (PCE), Tricloroeteno (TCE) y Cloruro de vinilo) e
hidrocarburos (BTEX); se propone evaluar y aplicar esta técnica de ser factible para
remediar las aguas subterráneas de este sector del antiguo botadero y de sus alrededores.
Desde mediados de los 90’s, las barreras reactivas permeables han suscitado más
atención y han recabado más estudio que ninguna otra técnica para descontaminar el
subsuelo (LJ.S. EPA, 1998). Como se muestra en la Figura 4.13, la PRB es una barrera
permeable de material especialmente escogido para que pasen las aguas contaminadas.
Figura 4.13 Barrera reactiva permeable para el tratamiento de una pluma contaminante en aguas
subterráneas
Fuente: Masters y Ela. Introducción a la Ingeniería Medioambiental. 2008. p: 275
En el tránsito a su través, los contaminantes son eliminados por uno o más procesos,
como la degradación abiótica redox, biotransformación, adsorción y precipitación.
Puesto que la PRB es sólo efectiva si la pluma lo atraviesa, la conductividad hidráulica
de la PRB debe ser mayor que la del acuífero en torno a ella. Si esto no sucede, o si la
PRB se obstruye por la filtración de partículas, precipitación de sólidos o concreciones
de productos de reacción sólidos, el agua se desviará y la rodeará.
179
El material más común usado es el hierro de valencia cero en limaduras (ZVI), que se
coloca típicamente en una barrera de 2 a 6 pies de espesor. El hierro en el estado de
oxidación cero (Fe(0)) se oxida rápidamente a ion ferroso (Fe(II)), y los electrones
perdidos pueden reducir destructivamente muchos de los contaminantes orgánicos
halogenados (ver Anexo 3). Usando tricloroetileno, como ejemplo, la descloración
reductiva completa puede seguir una reacción global tal como:
De forma similar a la Ecuación 4.7, el puede ser usado para tratar ciertos iones
inorgánicos sensibles-redox como el cromo desde el estado Cr(VI) al menos móvil
Cr(III).
Para eliminar otros contaminantes se pueden emplear otros materiales solos o junto con
él. Del mismo modo se pueden incorporar a las barreras reactivas nutrientes destinados
al crecimiento microbiano y con ello conseguir a la vez una biorrecuperación y la
eliminación de algunos contaminantes tanto orgánicos como inorgánicos. Los microbios
estimulados por las PRB no solamente se han usado para degradar directamente un
contaminante sino también para producir reactivos, como el sulfuro, que, o bien
reacciona con él o lo hace precipitar. Si se usa cal en la barrera el pH puede subir
suficientemente cómo para que ciertos metales precipiten. En otras aplicaciones se ha
utilizado un sorbente, el hidróxido férrico, para adsorber contaminantes como el arsénico
o el plomo (U.S. EPA, 1998).
Las barreras reactivas son atractivas porque suelen tener unos costes relativamente
bajos, permiten limpiar ciertos contaminantes difíciles de tratar por otros medios y no
necesitan mantenimiento una vez instaladas. Sin embargo sus prestaciones a largo plazo
180
no siempre han sido completamente satisfactorias. Las PRB no han funcionado o lo han
hecho por debajo de las expectativas por muchas razones: Perdida de permeabilidad,
pérdida de hierro u otros reactivos, cinética de reacción menor de la esperada,
obstrucción del medio o inactivación por biomasa generada por los microorganismos, y
mala conexión hidráulica entre el acuífero y la zona reactiva permeable.
181
CAPITULO 5. Conclusiones y recomendaciones
5.1. Conclusiones
• Dadas las características de cada método, para nuestro país el más factible es el
método GOD, ya que su aplicación requiere de menor cantidad de información y
sus resultados son confiables.
183
5.2. Recomendaciones
• Para ubicar cualquier tipo de proyecto, se debe realizar un estudio que determine
las areas de protección de pozos existentes en la zona, con el fin de proteger el
recurso hídrico subterráneo.
184
• Dado que con el presente estudio se ha logrado conocer la sensibilidad que
muestra la zona a la contaminación de las aguas subterráneas, se recomienda
realizar un estudio que determine el riesgo generado por los contaminantes
generados por la presencia de este antiguo botadero en este lugar.
185
BIBLIOGRAFIA
Beatriz Gil. (2007) Mapa y estudio de Recarga Acuífera del Área Metropolitana
de San Salvador, zona prioritaria No 1: Faldas del Volcán de San Salvador-Cuenca del
Río Acelhuate. El Salvador.
Escobar, M., Garcia, J., & Guerrero, T. (2002). Metodología para la evaluación
de la vulnerabilidad de las aguas subterráneas aplicada en el Valle de Zapotitán.
Antiguo Cuscatlan, El Salvador.
187
ITC. (2007). proyecto UNESCO-RAPCA. Recuperado 22 de Junio 2010, de
http://www.itc.nl/external/unesco-rapca/Pagina%20Casos%20de%20estudio.html
Johnson División OUP Inc. (1975). El agua subterránea y los pozos. Minnesota.
188
ANEXOS
ANEXO 1
190
ANEXO 2
191
ANEXO 3
Tabla A3. Los veinticinco contaminantes más frecuentemente detectados en el agua subterránea de
lugares con residuos peligrosos
192
ANEXO 4
Cuando los LNAPL entran en la zona insaturada, algunos pueden disolverse en agua en
los poros, otros pueden evaporarse y acabar mezclados con el aire de los espacios
porosos, y otros pueden ser adsorbidos por las partículas del suelo. Cuando alcanzan el
193
nivel freático no se disuelven bien y se difunden hasta formar una capa de contaminante
flotando sobre la superficie de la zona saturada, como se representaba en la Figura 4.4,
en la Figura 4.4(a). Los DNAPL también se hunden, se esparcen, adsorben y evaporan
en la zona insaturada, pero cuando alcanzan el nivel freático siguen hundiéndose hasta
que alcanzan una capa de material relativamente impermeable, donde forman estanques
que pueden desbordarse y seguir su curso hasta la siguiente capa impermeable, como se
ve en la Figura 4.4(b). Partes de los DNAPL forman pequeños glóbulos que acaban
atrapados en los espacios porosos, haciendo virtualmente imposible sacarlos mediante el
bombeo. Como quedan alojados en pequeñas fisuras y se acumulan en estanques al
fondo de los acuíferos se van disolviendo lentamente al pasar el flujo de agua
subterránea. Su solubilidad es tan baja que pueden transcurrir siglos antes de ser
completamente arrastrados. Los intentos de eliminar los DNAPL extrayendo agua son
completamente vanos.
194
ANEXO 5
GLOSARIO
Aguas subterráneas: Es el agua que circula en la zona saturada de los acuíferos y que
es posible de captar por obras de ingeniería.
Aguas superficiales: Son las aguas que discurren por la superficie de la tierra.
Agua vadosa: Agua somera que se encuentra entre la capa freática y la superficie.
Capilaridad: Es una propiedad de los líquidos que depende de su tensión superficial que
le confiere la capacidad de subir o bajar por un tubo capilar.
Desechos sólidos: Es cualquier material sólido acumulado como basura. Por ejemplo,
ropa, objetos plásticos, metálicos o cerámicos, latas de bebidas, llantas, poliestireno,
poliuretano, materiales de construcción, etc.
195
Espacios intersticiales: Espacios presentes entre las partículas que conforman el suelo.
Lixiviado: Es el líquido que percola a través de los residuos sólidos y que acarrea
materiales disueltos o suspendidos. La infiltración de agua de lluvia es el principal
generador de lixiviados en los rellenos sanitarios y en los botaderos.
196
Mapas de isopiezas: Refleja la forma de la superficie freática.
Nivel piezométrico: Punto concreto de una superficie virtual que alcanzaría el agua si se
hicieran infinitas perforaciones en un acuífero.
Ordenamiento Territorial: Es una normativa, con fuerza de ley, que regula el uso del
territorio, definiendo los usos posibles para las diversas áreas en que se ha dividido el
territorio, ya sea: el país como un todo, o una división administrativa del mismo.
Partes por millón (ppm): Es una unidad de medida de concentración de una solución.
Se refiere a la cantidad de mg (miligramos) que hay en un kg (kilogramo) de disolución.
197
Pluviometría: Estudio y tratamiento de los datos de precipitación que se obtienen en los
pluviómetros ubicados a lo largo y ancho del territorio.
Precipitación: Caída de agua sólida o líquida por la condensación del vapor sobre la
superficie terrestre.
198
Vulnerabilidad de las aguas subterráneas: Naticional Research Council (1993) define
vulnerabilidad de las aguas subterráneas como la tendencia o la probabilidad de que los
contaminantes puedan llegar a una posición específica en el sistema de aguas
subterráneas después de la introducción de contaminantes en algún lugar por encima del
acuífero superior.
Zona insaturada: También llamada zona vadosa o zona de aireación, esta comprendida
entre la superficie del terreno y el nivel freático del agua en los acuíferos y en ella los
poros no están completamente ocupados por agua.
Zona saturada: Zona del subsuelo en la que todos los huecos están ocupados por agua,
por debajo del nivel freático.
199