Tesis Defendida Por: Ahiezer Alarcón Barradas

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 115

TESIS DEFENDIDA POR

Ahiezer Alarcón Barradas


Y APROBADA POR EL SIGUIENTE COMITÉ

Dr. Rogelio Vázquez González


Director del Comité

Dr. Marco Antonio Pérez Flores Dr. José Manuel Romo Jones
Miembro del Comité Miembro del Comité

Dr. Jesús Arturo Martín Barajas Dra. Idania Valdez Vázquez


Miembro del Comité Miembro del Comité

Dr. Antonio González Fernández Dr. David Hilario Covarrubias Rosales


Coordinador del programa de posgrado Director de Estudios de Posgrado
en Ciencias de la Tierra.

24 de Septiembre de 2010.
CENTRO DE INVESTIGACIÓN CIENTÍFICA Y DE EDUCACIÓN SUPERIOR
DE ENSENADA

PROGRAMA DE POSGRADO EN CIENCIAS


EN CIENCIAS DE LA TIERRA

Exploración geofísica en terrenos del Exbasurero de Ensenada, B.C.

TESIS
que para cubrir parcialmente los requisitos necesarios para obtener el grado de
MAESTRO EN CIENCIAS

Presenta:
AHIEZER ALARCÓN BARRADAS

Ensenada, Baja California, México, Septiembre de 2010.


I
Resumen de la tesis de Ahiezer Alarcón Barradas, presentada como requisito
parcial para la obtención del grado de Maestro en Ciencias, en Ciencias de la
Tierra con orientación en Geofísica Aplicada. Ensenada, Baja California.
Septiembre de 2010.

Exploración geofísica en terrenos del Exbasurero de Ensenada, B.C.

Resumen aprobado por:


________________________________
Dr. Rogelio Vázquez González

En este trabajo se evalúa el problema de infiltraciones de lixiviados del


Exbasurero de Ensenada B.C. mediante técnicas geofísicas. Con este estudio se
establece la ubicación espacial de las infiltraciones del sitio, cuantificación de los
residuos que ahí yacen y una primera etapa de evaluación posclausura.
Se empleó el método de Potencial Espontáneo (PE) a manera de inspección
general del sitio y reconocimiento de zonas anómalas. Los datos se aproximaron a
un potencial residual y se interpretaron de forma cualitativa y cuantitativa por
medio de métodos geométricos simples.
A su vez los resultados se corroboraron con métodos más especializados como la
Tomografía Eléctrica Resistiva en dos dimensiones (TER 2D). Se realizaron 9
perfiles y un análisis de la conveniencia de los distintos tipos de arreglos
eléctrodicos. Se obtuvieron modelos de resistividad por medio de los programas
RES2DINV y Multielectrodico RESIST6 CICESE.
Los resultados de ambos métodos muestran una zona de polaridad negativa y alta
conductividad eléctrica hacia la parte norte del sitio, que es correlacionada con
presencia de lixiviados.
Ambos métodos probaron ser herramientas eficaces en la investigación de
infiltraciones acuosas con contraste eléctrico resistivo.

Palabras Clave: Exbasurero de Ensenada B.C., Potencial Espontáneo (PE),


Tomografía Electrica Resistiva 2D (TER 2D), lixiviados, RES2DINV, Multielectrodo
RESIST6 CICESE.
II
Abstract of the thesis presented by Ahiezer Alarcon Barradas as a partial
requirement to obtain the Master degree in Earth Science, with orientation in
Applied Geophysics. Ensenada, Baja California, México. September, 2010.

Exploration geophysics on the Old Dump site Ensenada, B.C.

In this tesis we evaluate the problem of infiltration of leachate from Old Dump
Ensenada's B.C. with geophysical techniques. This research establishes the
spatial location of the infiltration in the site, the quantification of waste, and a first
stage of post-closure assessment.
We used the method of Spontaneous Potential (SP) as an overall site inspection
tool and recognition of anomalous zones. The data were close to a residual
potential, and were interpreted qualitatively and quantitatively by simple
geometrical methods.
At the same time the results were corroborated with more specialized methods
such as the two-dimensional Electrical Resistive Tomography (ERT 2D). Nine
profiles and an analysis of the desirability of different types of electrodes arrays
were perfomed. Resistivity models were obtained through programs RES2DINV
and RESIST6 Multielectrode CICESE.
The results of both methods show a area with negative polarity and high electrical
conductivity towards the northern part of the site, correlated with presence of
leachate.
Both methods proved effective tools in the investigation of aqueous infiltration with
a resistive electric contrast.

Keywords: Old Dump Ensenada's B.C., Spontaneous Potential (SP), Electric


Resistive Tomography 2D (ERT 2D), leachate, RES2DINV, RESIST6
Multielectrodo CICESE.
III
Dedicatorias

A mis padres, Nieves Barradas Ramirez y Antonio Alarcón Alfonseca, por todo su
apoyo, fuerza y dedicación...

A mis hermanos, Misael e Ivonne.


IV
Agradecimientos

Al Dr. Rogelio Vázquez González por su apoyo, dedicación y paciencia en el


desarrollo de esta tesis.

A los integrantes de mi comité: Dr. José M. Romo Jones, Dr. Marco A. Pérez
Flores, Dr. J. Arturo Martín Barajas y la Dra. Idania Valdez Vázquez por sus
importantes y enriquecedores comentarios, muestra del interés de la tesis.

A todo el personal de Ciencias de la Tierra, por su ayuda y especialmente a


Martha E. Barrera López, Guadalupe Zepeda Tapia, Ma. Bábara Uribe Martínez.

Al personal técnico Alejandro Díaz Fernández, C. Enrique Castillo Guerrero, Mario


Vega Aguilar, José de J. Mojarro Bermúdez, Humberto S. Benítez Pérez y Jaime
Calderón González.

A los alumnos del Verano Científico 2009, por el trabajo realizado en el sitio.

A mis cuates: José (el Koopa), Ismael (el Ratón), Armando (el Chat), Jesús (el
Chuy marihuano), Román (el Musolini),

Al Centro de Investigación Científica y de Educación Superior de Ensenada


(CICESE) y al Consejo Nacional de Ciencia y Tecnología (CONACYT) por la
oportunidad de realizar mi Maestría y el apoyo económico otorgado.
V
CONTENIDO

Página
Resumen español……………………………………...……...…………... I
Resumen ingles…………………………………………………...……….. II
Dedicatorias………………………………………………………..……...... III
Agradecimientos…………………………………………………..……….. IV
Contenido…………………………………………………………..……….. V
Lista de Figuras…………………………………………………….…..….. VII
Lista de Tablas……………………………………………………….……... XII

Capítulo I. Introducción …................................................................... 1


I.1 Objetivos..................................................................................... 5
I.2 Área de estudio........................................................................... 5
I.3 Fisiografía y Geología................................................................. 7
I.4 Situación del Exbasurero............................................................ 9
I.5 Antecedentes.............................................................................. 9

Capítulo II. Método de Potencial Espontáneo…………………........... 11


II.1 Teoría......................................................................................... 12
II.2 Metodología............................................................................... 15
II.3 Interpretación............................................................................. 18
II.3.1 Interpretación Cualitativa........................................................... 18
II.3.1.1 Corrección por deriva temporal................................................... 20
II.3.1.2 Corrección por Topografía........................................................... 23
II.3.2 Interpretación Cuantitativa.......................................................... 28
II.4 Resultados.................................................................................. 34

Capítulo III. Método eléctrico de corriente continua.…..................... 37


III.1 Teoría......................................................................................... 39
III.2 Metodología............................................................................... 44
III.2.1 Instrumentación.......................................................................... 53
III.3 Inversión 2D Tomografía Eléctrica Resistiva.............................. 54
III.3.1 Inversión 2D por Método de Elemento Finito............................. 55
III.3.2 Inversión 2D Multielectródica por Método de Ecuación
Integral.................................................................................................... 56
III.4 Resultados e Interpretación....................................................... 59
III.4.1 Línea 1....................................................................................... 59
III.4.2 Línea 2....................................................................................... 61
III.4.3 Línea 3....................................................................................... 61
III.4.4 Línea 4....................................................................................... 63
III.4.5 Línea 5....................................................................................... 63
III.4.6 Línea 6....................................................................................... 66
VI
CONTENIDO (continuación)

Página

III.4.7 Línea 7....................................................................................... 69


III.4.8 Línea 8....................................................................................... 71
III.4.9 Línea 9....................................................................................... 71
III.4.10 Secciones de resistividad 3D..................................................... 74

Capítulo IV. Discusiones ….................................................................. 76

Capitulo V. Conclusiones y recomendaciones................................... 84

Bibliografía……………….………........................................................... 86

Apéndices…........................................................................................... 93
I Datos de campo de Potencial Espontáneo con errores de cierre..... 93
Programa de Modelación-Inversión Matlab de PE residual.............. 94
II Modelos observados y calculados de TER 2D................................. 96
III Análisis químico básico de muestras de agua.................................. 99
IV Cálculo de volumen de residuos del Exbasurero.............................. 101
V Memoria fotográfica........................................................................... 102
VII
LISTA DE FIGURAS

Figura Página
1 Mapa de ubicación de la zona de estudio, Exbasurero de
Ensenada que se localiza a 6.6 km de la ciudad.................... 6

2 a) Perfil escavado de 7m en granito intemperizado que


muestra el suelo residual tipo grus 1. b) Mapa geológico del
área de estudio. (INEGI).......................................................... 8

3 Técnica de campo de gradiente para la medición de PE........ 16

4 Técnica de campo de campo total o de base fija.................... 16

5 Secuencia de medición de PE en campo................................ 18

6 Situación análoga en campo por condiciones redox,


metabolismo de materia orgánica........................................... 20

7 Espectro de frecuencia de la intensidad de campo


magnético y su correspondiente espectro de campo eléctrico
inducido para un modelo uniforme de resistividad de 20
Ohms-m................................................................................... 21

8 Variación de la componente Hp paralela al eje de rotación


de la Tierra en los días de medición y ajuste lineal de la
deriva....................................................................................... 22

9 Corrección por deriva temporal de observaciones de PE....... 23

10 Interacción iónica interfase mineral-solución y variación del


potencial. Doble capa: (1) Capa de Helmholtz. (2) Capa de
Gouy-Chapman....................................................................... 24

11 Distribución general del nivel freático...................................... 25

12 Relación entre la elevación topográfica y la señal de PE........ 26

13 a) Corrección del efecto topográfico de PE. b) Relieve


topográfico y ubicación del arroyo Ensenada.......................... 27
VIII
LISTA DE FIGURAS (continuación)

Figura Página
14 Izquierda. Modelo sintético de PE. Derecha. Inversión del
modelo sintético obtenido por análisis de segunda
derivada................................................................................. 31

15 Mapa de ubicación del perfil de anomalía de PE


residual................................................................................... 32

16 a) Solución del perfil de anomalía de PE en el Exbasurero


debido al efecto de dos formas simples (q 1=0.8≈cilindro
horizontal y q2=0.5=cilindro vertical). b) Ajuste del perfil a
los efectos de las formas y zona de interacción
considerada............................................................................ 33

17 Mapa de contornos de PE residual. Resaltan valores


positivos en la parte suroeste y valores negativos en la
parte norte a partir de los pozos de
respiración.............................................................................. 34

18 Desarrollo de las zonas redox en una pluma contaminante


por lixiviados donde Eh es el potencial de reducción
Eh1<Eh2<Eh3<Eh0............................................................... 35

19 Ajuste de los datos de PE residual observado y el


calculado................................................................................ 36

20 Electrodo puntual de emsión de corriente en un medio


homogéneo de resistividad ρ................................................. 41

21 Determinación de resistividad de un medio, arreglo


tetraeléctrodico. Líneas de corriente (azul) y de potenciales
(negras discontinuas)............................................................. 43

22 Arreglos electródicos utilizados en este trabajo..................... 45

23 Simulación de cobertura horizontal de datos con 28 electros


para arreglo electródico Wenner, Schlumberger y Dipolo-
dipolo, utilizados para este trabajo........................................ 46

24 Niveles y puntos de atribución de la tomografía tipo


Schlumberger en el relevamiento electródico........................ 47
IX
LISTA DE FIGURAS (continuación)

Figura Página
25 Profundidades medias de investigación sobre un medio
homogéneo obtenidas para una distancia eléctrodica de
10m a los diferentes niveles de medición (Edwards, 1977)... 49

26 Función de sensibilidad para arreglo Wenner,


Schlumberger (n=2) y Dipolo-dipolo (n=4) correspondiente
a un modelo homogéneo de resistividad, (tomada de López
Hidalgo et al. 2007)................................................................ 50

27 Decaimiento de la señal para un modelo homogéneo de


100 Ohms-m, I=1.5A y a=10m.............................................. 51

28 Mapa de ubicación de líneas de TER en el Exbasurero de


Ensenada............................................................................... 52

29 Consideración de distribución de electrodos de transmisión


y recepción para el método de Ecuación Integral.................. 58

30 Modelos de la Línea 1. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 60

31 Modelos de la Línea 2. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 62

32 Modelos de la Línea 3. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 64

33 Modelos de la Línea 4. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 65
X
LISTA DE FIGURAS (continuación)

Figura Página
34 Modelos de la Línea 5. a,b) Resultados por el programa
RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 67

35 Modelos de la Línea 6. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 68

36 Modelos de la Línea 7. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 70

37 Modelos de la Línea 8. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 72

38 Modelos de la Línea 9. a,b) Resultados por el programa


RES2DINV para el arreglo Wenner y Schlumberger
respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.................................................................................. 73

39 Cortes de resistividad del terreno del Exbasurero de


Ensenada a -3.5, -8.5, -16, -30, -45 y -62.5 metros de
profundidad............................................................................ 75

40 Relación entre volumen contenido de agua y resistividad


(Greller et al. 2005)................................................................ 78

41 Comparación entre el mapa de contornos de PE residual y


el modelo de resistividad a una profundidad de 10m........... 79

42 Izquierda. Evolución temporal de volumen de lixiviado en


un sitio de disposición (Farquhar, 1989). Derecha.
Tendencia de la concentración de lixiviados con el tiempo a
partir de la clausura................................................................ 80
XI
LISTA DE FIGURAS (continuación)

Figura Página
43 Modelo estático de capas subhorizontales de PE residual,
equivalencia de circuito eléctrico........................................... 82

44 Comparación entre la interpretación cuantitativa de PE


residual y TER........................................................................ 83

45 Modelos observados y calculados de TER 2D. a) Para los


arreglos Wenner y Schlumberger, respectivamente por el
programa RES2DINV. b) Gráfico de ajuste de los datos al
modelo con coeficientes de correlación de Pearson y
Spearman por el programa RESIST6 Multielectrodo
CICESE.................................................................................. 96

46 Localización de muestras obtenidas en el aroyo Ensenada.. 99

47 Topografía actual y anterior del Exbasurero de Ensenada. 102


XII
LISTA DE TABLAS

Tabla Página

I Parámetros obtenidos por descomposición de formas del


36
análisis cuantitativo de PE residual.....................................

II Resistividades de las rocas más comunes, algunos


37
minerales y químicos...........................................................

III Especificaciones técnicas del Super Sting AGI 54

IV Observaciones de campo de PE miliVolts (mV). Circuito


Exterior se refiere al polígono realizado en el Exbasurero
con las respectivas direcciones. Líneas Interiores se 93
refiere a los perfiles al interior este polígono.(E de
C.=Error de cierre)..............................................................

V Ubicación, altitud (m.s.n.m.m), nivel freático estático (nF),


conductividad eléctrica (CE) y sólidos disueltos totales
99
(STD) de las muestras colectadas sobre el Arroyo
Ensenada.............................................................................

VI Concentración de cationes y aniones principales en


partes por millón (ppm) de las muestras colectadas sobre 100
el arroyo Ensenada y del Lixiviado......................................
Capítulo I

Introducción

La falta de mecanismos adecuados para el manejo y control de los residuos


sólidos generados por una ciudad derivan en una gama de problemas y sitios
improvisados que, con el tiempo, se convierten en grandes estructuras con gran
cantidad de sustancias, entre ellas tóxicas peligrosas, que no se degradan con
facilidad y persisten en la naturaleza. Particularmente los residuos sólidos
municipales (RSM) se colocan, desde hace 30 años e incluso actualmente, a cielo
abierto sobre barrancas, cañadas, terrenos baldíos y márgenes de cuerpos de
agua.

La importancia del agua para el desarrollo de los ecosistemas así como para las
actividades productivas del hombre es indiscutible, pero la irresponsable y
negligente contaminación de este recurso justifican la necesaria atención para su
protección y para asegurar el acceso este líquido elemental.

El agua subterránea es una fracción significativa de los recursos hídricos, a nivel


global representa unas veinte veces más que el total de las aguas superficiales de
todos los continentes e islas y en muchas partes del mundo, donde la lluvia es
escasa, constituye la fuente fundamental de abastecimiento, de ahí la importancia
en su cuidado como reserva de agua dulce. El hecho de que el agua subterránea
yace en el subsuelo limita el control de los mecanismos de interacción, por lo que
es difícil su gestión, lo cual la hace susceptible a la contaminación y
sobreexplotación cuando se carece de un plan de aprovechamiento sostenible.
2
Resulta alarmante la inserción a los sistemas hídricos de sustancias tóxicas
(fertilizantes, detergentes, productos químicos, etc.) lo suficientes para causar
efectos dañinos mensurables en la flora o fauna (incluido el hombre) de manera
directa o indirecta, que alteran la calidad y limitan la disponibilidad del agua,
comprometiendo uno de los bienes básicos de conservación.

El problema de la contaminación del agua tiene múltiples perspectivas (social,


ecológica, biológica, geohidrológica, etc.) y se presenta en formas diversas tanto
geológico-naturales como inducidas (físico-químicas y biológicas), con
asociaciones y sinergismos difíciles de preveer, con efectos que derivan en un
abanico de efectos negativos. El aumento desmedido de desechos
antropogénicos han provocado la improvisación de sitios sin el más mínimo
cuidado a lo recursos hídricos subterráneos, comprometiendo su disponibilidad.

Durante la Revolución Industrial, empezaron los problemas de aglomeración de


población y por ende problemas de saneamiento de las ciudades. En 1848, Gran
Bretaña publica la Ley de Salud Pública e inicia el proceso de regulación de
residuos. En la década de 1970 se establecieron normas ambientales para la
eliminación de residuos y bases para la construcción de vertederos con
revestimiento geológico. Fue hasta 1974, cuando se reconoció por primera vez el
problema de los lixiviados generados en los vertederos.

Los lixiviados son líquidos, incluidos todos los componentes en suspensión, que
han percolado o drenado a través de los residuos sólidos. En este contexto los
residuos sólidos pueden ser definidos de manera más amplia y podría incluir por
ejemplo un almacén de sal o suelo contaminados. En este trabajo nos
enfocaremos en los lixiviados procedentes de los vertederos de basura y, en
particular de los residuos sólidos urbanos (RSM).

Los daños potenciales debidos a los lixiviados se producen cuando sustancias


3
tóxicas se concentran, penetran y se mueven a través del subsuelo hasta
interactuar con el sistema hídrico (río, lago o manto acuífero), o bien pueden
permanecer en el subsuelo como amenaza potencial, pudiendo modificar la
calidad del agua y limitar su uso. No obstante, el daño más importante es su
potencial para reducir el contenido de oxígeno disuelto (OD) dentro de cualquier
curso o fuente de agua.

En México es prácticamente reciente la normalización de vertederos municipales


de acuerdo a la legislación NOM-083-ECOL-1994, y la más reciente en vigor, la
NOM-083-SEMARNAT-2003, que establece las especificaciones de selección del
sitio, el diseño, construcción, operación, monitoreo, clausura y obras
complementarias de un sitio de disposición final de residuos sólidos urbanos y de
residuos que requieren manejo especial.

En los planes de clausura de los vertederos, las herramientas geofísicas ofrecen


una fuente de información rápida y económica para la caracterización de estos
sitios. Particularmente se usan en la detección de filtraciones de fluidos
contaminantes subterráneos, ya que proveen datos útiles que permiten evaluar su
distribución y extensión espacial en el sitio. Y es en el caso del Exbasurero de la
cuidad de Ensenada, B.C. donde se aplicaron con este fin. Informes técnicos y
testimonios recabados señalan un tiempo de uso mayor a 40 años y un foco
latente de contaminación (Vázquez González, 2005). El acondicionamiento del
lugar en donde se ubica el Exbasurero corresponde a la época en donde no se
preparaban adecuadamente los sitios de depósito, ni mucho menos se construía
obra civil para el diseño y operación de los residuos (Moreno et al. 2004).

En este sitio, debido a sus características no planeadas, los subproductos que se


generan por la descomposición, degradación y acción bacterial de los materiales
que ahí yacen forman diversos compuestos, entre ellos lixiviados que migran en
forma aleatoria, alterando la composición del suelo-agua de manera acumulativa
4
(Vázquez González, 2007; Noguera Cáceres, 2003).

El presente trabajo forma parte importante del proyecto “Identificación de suelos


contaminados utilizando métodos geofísicos” del Departamento de Geofísica
Aplicada del Centro de Investigación Científica y de Educación Superior de
Ensenda (CICESE), en donde se propone evaluar, mediante la aplicación de los
métodos geofísicos, zonas contaminadas, utilizando métodos eléctricos y
electromagnéticos. En el caso del Exbasurero de Ensenada, B.C., se realizó una
primera evaluación del sitio, mediante dos técnicas geofísicas: el método de
Potencial Espontáneo (PE) y la Tomografía Eléctrica Resistiva en dos dimensiones
(TER 2D).

Esta tesis aborda el problema de la detección de lixiviados en el subsuelo con la


siguiente capitulación: en el capítulo II se describe el método de Potencial
Espontáneo (PE), y se interpretan los datos tomados en el Exbasurero durante el
mes de septiembre del 2009, tanto en forma cualitativa como cuantitativamente.

El capítulo III se trata el problema con Tomografía Eléctrica Resistiva en dos


dimensiones (TER 2D), hecha en el lugar durante el mes de noviembre del 2009,
con configuración electródica Wenner y Schlumberger. Se obtiene una imagen de
la distribución de resistividad del subsuelo, mediante los programas de inversión
RES2DINV (Loke y Barker, 1996), y el programa RESIST 6 Multielectrodo CICESE
(Pérez Flores, 1995; Pérez et al. 2001).

En el capítulo IV, se discuten los resultados de estos dos métodos y se realiza la


integración de la información para la construcción de un modelo único que
evidencia la existencia de infiltraciones de lixiviados.
Finalmente en el capítulo V, se presentan las conclusiones y recomendaciones con
base en los resultados obtenidos.
5
I.1 Objetivos

Objetivo particular: En este trabajo se investiga la distribución y el alcance espacial


de las infiltraciones de lixiviados del antiguo basurero de Ensenada B.C. mediante
métodos eléctricos geofísicos, particularmente Potencial Espontáneo (PE) y
Tomografía Eléctrica Resistiva 2D (TER 2D). Además se hace un análisis de la
conveniencia de cada una de las técnicas empleadas.

Objetivo general: Proporcionar una base técnica sólida de evaluación con estudios
geofísicos posclausura de Exbasurero de Ensenada B.C., como base para una
segunda etapa de estudios geoquímicos dirigidos a desarrollar un plan de
estabilización-remediación del sitio.

I.2 Área de estudio

El Exbasurero de Ensenada, B.C. tiene un área de 6.6 hectáreas; se localiza al


noroeste de esta ciudad, entre los paralelos 31° 52’ a 31° 53’ de latitud norte y los
meridianos 116° 33’ a 116° 34’ de longitud oeste, hacia el norte, a 300 m, colinda
con el arroyo Ensenada. (Figura 1). El acceso se tiene por la carretera federal
número 3 Ensenada-San Felipe. Se encuentra rodeado de 4 colonias que fueron
desarrollándose en los últimos 20 años (Fraccionamiento Mar, Ampliación Gómez
Morín, Popular 89 y 17 de Abril).
6

da
nsena
E
rroyo
A

Figura 1.Mapa de ubicación de la zona de estudio, Exbasurero de Ensenada que se


localiza a 6.6 km al este de la ciudad.
7
I.3 Fisiografía y Geología

De acuerdo a la clasificación de las provincias fisiográficas realizada por el INEGI


(1997), la zona de estudio se encuentra en la subprovincia de la Sierra de Baja
California, a su vez contenida en la Provincia de Baja California. Esta subprovincia
se caracteriza por un basamento formado por rocas batolíticas graníticas,
metamórficas mesozoicas cubiertas por depósitos sedimentarios recientes (Raisz,
1964). La zona urbana está asentada prácticamente sobre conglomerados con
clastos y arena de origen granítico del Cretácico Tardío y Terciario Temprano, y los
arroyos presentan terrazas de origen fluvial del Pleistoceno (CNA, 2002).

El subsuelo de Ensenada consiste en un basamento conformado por rocas


metavolcánicas de la Formación Alisitos (Almazán Vázquez, 1988). Este cuerpo
basal, aflora al oriente de la zona costera, profundizándose conforme se adentra
hacia la playa. La configuración obtenida a partir de métodos potenciales, muestra
que el basamento y los depósitos volcánicos de la Formación Alisitos se
profundizan hacia la zona del mar y se encuentran subyaciendo a la cuenca
sedimentaria de la Formación Rosario en los extremos norte y sur (Gastil et al.
1975); la Formación Rosario es predominantemente arenosa, sin embargo tiende
a mostrarse más conglomerática en los sectores central y meridional, lo cual se
considera como la zona donde la unidad presenta los mayores índices de
permeabilidad.

Cubriendo discordantemente a la Formación Rosario, se encuentra una secuencia


de materiales granulares con buenos índices de permeabilidad, representados por
arenas de grano medio a fino, arcillas y cuerpos de conglomerados, los cuales
fueron depositados en los causes de los arroyos provenientes de las zonas
montañosas, así como en las zonas de talud al pie de las prominencias
topográficas; los espesores de esta secuencia varían de 30 a 80 m (CNA, 2002).
8
Específicamente, el sitio se encuentra en una capa arenosa somera derivada del
intemperismo del gránito Cretásico denominada como grus 1 (Figura 2a).
a)

b)

Figura 2. a) Perfil excavado de 7m en granito intemperizado que muestra el suelo residual


tipo grus1. b) Mapa geológico del área de estudio. (INEGI, 1974).

1 Término usado para describir productos granulares in-situ derivados de la desintegración de rocas cuarzo-
granulares, muy característico del granito. (Migón, 2006).
9
I.4 Situación del Exbasurero

El sitio de estudio funcionó como vertedero municipal alrededor de 40 años, fue


clausurado de manera oficial en el 2004. En el 2009, el XIX Ayuntamiento de
Ensenada propuso la construcción de un cementerio, pero a falta de estudios de
factibilidad éste se derogó (Madrigal, 2009).

Actualmente, sin estudios técnicos y ambientales previos, el XIX Ayuntamiento de


Ensenada pretende nuevamente construir un panteón en el predio. Por otra parte,
la Secretaría de Protección al Ambiente del Estado mantiene un proceso
administrativo en contra del Gobierno Municipal (Cruz Aguirre, 2008) debido a que
se ha ignorado la NOM-083-SEMARNAT-2003, que establece especificaciones
tales como un programa posclausura por un periodo de 20 años, así como de un
programa de monitoreo por la emisión de biogás y generación de lixiviados.

I.5 Antecedentes

Entre los trabajos previos consultados se encuentran:


– Cartas emitidas por INEGI como:
Carta Topográfica 1:50,000; Ensenada H11B12. INEGI.
Carta Geológica 1:50,000; Ensenada H11B12. INEGI.
– Estudio realizado por el CICESE con la finalidad de detectar una posible
capa arcillosa somera en el Exbasurero (Gómez Treviño y Espinoza, 1990).
– Determinación de la disponibilidad de agua en el acuífero de Ensenada,
Edo. de B.C. por la Comisión Nacional del Agua (CNA), Subgerencia de
Evaluación y Modelación Hidrogeológica (CNA, 2002).
– Proyecto para la exploración y aprovechamiento de biogás producto de la
descomposición de las millones de toneladas que alberga, verano de 2006,
10
con una inversión de 50 mil dólares se buscó extraer gas metano del
Exbasurero para la generación de electricidad, hecho por Secretaría de
Ecología del Edo. y aprobado por la Comisión de Cooperación Ecológica
Fronteriza (COCEF), la Agencia de Protección medioambiental
estadounidense (EPA), el Banco de Desarrollo América del Norte y el
Gobierno Federal y Estatal (Cruz Aguirre, 2008).
– Informe de Actividades Primer Semestre 2007, Proyecto EPA-STAG,
Contaminación por lixiviados en Basureros: Caracterización,
Biorremediación y Prevención. CICESE-UAM-Universidad de Arizona
(Vázquez González, 2007).
– Estudios para valorar gases, toxicidad y otros parámetros con el fin de
comprobar la estabilización del lugar, por la Secretaría de Protección al
Ambiente.

Se piensa que existe una cantidad importante de residuos que yacen en el lugar,
dado que Ensenada para 1990 contaba 259 979 habitantes (COESPO, 2009).
Considerando una población constante, sólo a manera de ilustración, durante los
40 años de actividad del basurero y tomando en cuenta que aproximadamente se
generan 0.939 kg/día per capita de residuos domésticos (INE-SEMARNAP, 1999).
Se calcula que aproximadamente se generaron 244 t/día, 89 000 t/año lo que
significa que alrededor de 2.7 millones de toneladas fueron acumuladas en el lugar
a lo largo de su funcionamiento, medianamente concordante con el Apéndice IV.
Capitulo II

Método de Potencial Espontáneo

El método de prospección de Potencial Espontáneo (PE) tuvo su origen en la


búsqueda de recursos minerales, posteriormente se convirtió en una herramienta
básica para el registro geofísico de pozos, consecutivamente su uso se amplió a
aplicaciones en ingeniería civil (Rozycki et al. 2006; Revil et al. 2007),
medioambiental (Naudet et al. 2008; Alatorre Zamora et al. 2008) y geotermia
(Manzella, 2000; Yasukawa et al. 2000). Es una herramienta útil, así como una
técnica auxiliar de varias aplicaciones debido a que es robusto, sencillo y rápido.

El método de PE o potencial eléctrico natural es una técnica pasiva, en el que se


miden diferencias de potencial eléctrico natural en la superficie del suelo, éstas
diferencias son debidas a mecanismos físico-químicos en los que los fluidos y sus
propiedades eléctricas juegan un papel importante (Orellana, 1982). Su principal
límitación es que no se tiene control de las fuentes que lo provocan, por lo tanto
su interpretación se hace un tanto difícil y es generalmente a nivel cualitativo
(Nyquist y Corry, 2002).

Además, el método es sensible a flujos de corriente que ocurren a profundidad del


terreno por procesos de oxido-reducción, razón por la que se uso este método en
el sitio de estudio de esta tesis, a fin de inspección general y para tratar de
identificar donde ocurre la degradación de la materia contenida y zonas de
infiltración de fluidos contaminantes (Naudet et al. 2004).

Se llevo a cabo una campaña de exploración de PE en el Exbasurero durante el


12
mes de septiembre del 2009, cubriendo un polígono de 560 m por 560 m y 3
líneas con espaciamiento electródico de 40 m con 98 estaciones de medición en
total.

Los datos se procesaron con la intensión de obtener el potencial residual e


interpretaron de forma cualitativa. También de manera cuantitativa por medio de
modelación-inversión de métodos fijos geométricos.

II.1 Teoría

El PE tiene como fuentes reconocidas que inducen este fenómeno:

– Potencial de Difusión y potencial de membrana. Es el potencial que está


asociado con diferencias en concentración de especies iónicas en el
subsuelo que provocan potenciales de difusión. Si los aniones y cationes
involucrados tienen diferente movilidad llegan a formarse dos zonas, en
cada una de las cuales predominan los iones de una carga, resultando
entonces diferencias de equilibrio que provocan un potencial eléctrico
(Orellana, 1982). El campo eléctrico resultante dependerá de la rapidez o
lentitud de movilidad iónica para alcanzar la electroneutralidad. En
equilibrio, el potencial de difusión Ed está dado por:

 I a −I c 
E d = RT ln C 1 /C 2  (1)
nF

donde: Ia y Ib son las mobilidades de los aniones y de los cationes


respectivamente, n es la carga eléctrica del ion, R es la constante universal
de los gases, T es la temperatura absoluta, F es la constante de Faraday y
C1 y C2 son las concentraciones de las soluciones que crean el gradiente
de difusión. Este mecanismo puede crear anomalías de decenas de
13
miliVolts (mV), y puede ser fuente de ruido en la mayoría de estudios de PE
(Nyquist y Corry, 2002).

– Potencial de membrana. Cuando se tienen arcillas en contacto con arenas,


se desarrolla un voltaje en el contacto, porque la arcilla es permeable a
cationes como el Na+ pero no lo es para aniones como el Cl -. La difusión del
Na+ de la arena en la arcilla crea un gradiente de potencial que tiende a
conducir los iones de Na+ de nuevo a la arena. El equilibrio entre la difusión
y la migración de iones eléctricamente conducidos viene dado por la
ecuación de Nernst, que es el potencial de difusión sin la diferencia de
movilidad de iones.
RT
Ed= lnC 1 /C 2  (2)
nF

– Potencial de Filtración. Surge cuando el agua o cualquier otro fluido circula


a través de una membrana porosa, se produce, entre los dos lados de ésta,
una diferencia de potencial. Las rocas pueden ser consideradas
membranas cuando poseen, por su porosidad, una red capilar a través de
las cuales pueden filtrarse las aguas subterráneas. Las paredes capilares
pueden absorber los aniones, las cuales atraen cationes, formándose una
capa eléctrica doble. Este fenómeno fue estudiado por primera vez por
Helmholtz en el siglo XIX. Para el flujo en un tubo capilar a través del cual
circula un electrolito el campo eléctrico, E, (V/m) esta dado por:


E= p (3)
4
donde e,r y m son la contante dieléctrica (F/m), la resistividad ( Ohms-m)

y la viscosidad (Pa-s) del electrolito, respectivamente; z es un parámetro


determinado por el material de la pared capilar y el electrolito llamado
potencial de doble capa: p es el gradiente de presión (Pa/m). El campo
eléctrico ocurre en la misma dirección del gradiente de presión, y en
14
sentido contrario al flujo del electrolito.

– Potencial Redox. Es provocado por reacciones geoquímicas de oxido-


reducción, equivalentes a la de una pila electroquímica. Sato y Mooney
(1960) propusieron el modelo para un depósito de sulfuros con una zona
por encima y otra por debajo del nivel freático, con reacciones de oxidación
que toman lugar por encima del cuerpo de sulfuros y reducción hacia la
parte de abajo del cuerpo, mismo que actúa como un cuerpo inerte
conductor. Así, las porciones alta y baja del depósito se comportan como
dos medias pilas galvánicas. La cima actúa como el ánodo con carga
negativa y el fondo como el cátodo con carga positiva. En respuesta a la
carga eléctrica, los iones migran a través de los poros de la roca
circundante del cuerpo mineral.

– Corrientes Telúricas. La radiación ultravioleta del Sol ioniza las móleculas


del aire en una delgada capa de la atmósfera llamada ionósfera, localizada
a una altura de 80 km sobre la superficie de la tierra, las corrientes
eléctricas en la ionosfera se forman sistemáticamente debido a movimiento
de los iones, los cuales son afectados por factores como las mareas diarias
y mensuales, variaciones estacionales, periodos de fluctuación del Sol
(manchas solares). Estas corrientes de iones producen campos magnéticos
que a su vez producen corrientes inducidas en la corteza terrestre llamadas
corrientes telúricas.

Existen otros fenómenos que también producen PE como el termoeléctrico,


presente en áreas volcánicas y geotérmicas en donde se han observado
anomalías considerables, debido a gradientes de temperatura y movimiento de
fluidos subterráneos iónizados, con amplitudes de decenas hasta centenas de mV
y polaridad positiva que reflejan mecanismos de polarización que ocurren a
profundidad (Finizola et al. 2004). El efecto bioeléctrico, por la acción de las raíces
15
de las plantas en la selección de iones y bombeo de agua que pueden ser notados
cuando hay cambios de vegetación, asociados a variaciones en la composición del
suelo o rocas subyacentes (Nyquist y Corry, 2002).

II.2 Metodología

El equipo necesario para la adquisición de PE consiste en un voltímetro de gran


impedancia interna (>10MOhms) y altamente sensible, éste mide la diferencia de
potencial eléctrico entre los dos eléctrodos plantados en tierra. Se deben usar
eléctrodos impolarizables, ya que las simples estacas metálicas distosionan la
señal por las reacciones electroquímicas entre el metal y la humedad del suelo.
Los electrodos se separan a cierta distancia (decena de metros) por medio de un
carrete de cable.

El equipo usado para el proyecto fue un voltímetro METEX M-4640A, con


sensibilidad de 0.01 mV y con una impedancia interna de 20 MOhms y electrodos
impolarizables Cu/CuSO4.

Existen dos formas comunes para la adquisición de datos en campo:

– El método de gradientes, emplea una separación constante entre los


electrodos M y N a lo largo del perfil; se recomienda alternar el orden de los
electrodos una vez tomada la lectura en la próxima estación, es decir, el
electrodo M pasa a ser el electrodo N y N - M, o bien, el electrodo 1 que
está al inicio del primer dipolo pasa a ser el último del siguiente dipolo
mientras el electrodo 2 permanece fijo y es hasta el siguiente dipolo donde
se moverá al final del dipolo 3 (Figura 3); esto para evitar remanentes de
polarización en los electrodos. El potencial total por tanto se obtiene por la
suma de las diferencias del potencial entre cada dipolo.
16

Figura 3: Técnica de campo de gradiente para la medición de PE.

– El método de campo total, en donde un electrodo permanece fijo a un punto


de partida y el otro se mueve en dirección de la línea de prospección
(Figura 4), la implementación en campo de esta técnica es impráctica
porque la longitud del cable aumenta conforme se avanza en el perfil. Su
aplicación principal es en registro de pozos.

Figura 4: Técnica de campo de campo total o de base fija

La técnica de campo usada para la recolección de datos de PE en este estudio fue


la de gradientes con dipolos de 40 m de separación. También se recomienda cavar
de 15 a 20 centímetros de profundidad para colocar el electrodo, asegurar su
17
contacto en el suelo y tiempo para su estabilización (20 segundos) antes de
efectuar la medición, para evitar efectos de potenciales bioeléctricos de superficie.

Un control de calidad para los datos de campo de PE, es la aplicación de la ley de


mallas de circuitos eléctricos o de tensiones de Kirchhoff que enuncia: “En toda
malla, la suma de todas las caídas de tensión es igual a la suma de todas las
subidas de tensión”, es decir, que en todo circuito cerrado la suma de diferencias
de potencial es igual a cero.

Para cumplir con el criterio anterior se realizan polígonos cerrados o bien líneas
de partida y de regreso (Figura 5), el fin es obtener un circuito cerrado que según
la ley de Kirchhoff, la suma algebraica deberá ser igual a cero. Sin embargo, en la
práctica no es así, se obtienen valores no nulos debidos a variaciones temporales
del campo eléctrico, tensiones parásitas y errores de medición que producen
errores de cierre (ecuación 4). Según las normas soviéticas el error de cierre debe
ser menor o igual a 5% (Orellana, 1982). Cabe mencionar que este error de cierre
debe repartirse entre todas las lecturas del polígono. Es importante señalar que
las lecturas deben estar referidas a un electrodo base, mismo que es colocado
fuera del área de estudio para observar si en verdad existen diferencias
significativas de potencial eléctrico en el sitio. Para este trabajo el error de cierre
promedio fue de 4.02% (Apéndice I).
100∣∑  V i∣
i
e= (4)
∑ ∣ V i∣
18

Área de
estudio

Figura 5. Secuencia de medición de PE en


campo

II.3 Interpretación

Descritas la causas que provocan el fenómeno de PE y la imposibilidad de


controlar parámetros como la intensidad, la frecuencia y profundidad de
investigación, el registro de campo obtenido se ve afectado por una gran cantidad
de variedades de PE integradas en un sóla lectura (voltaje), ya que en la
naturaleza se combinan efectos de redox, filtración, difusión, deriva del campo
eléctrico natural (corrientes telúricas) y si es el caso, efecto termoeléctrico,
además de efectos perturbadores de sitio como la corrosión de elementos
químicos, actividad biológica, líneas de transmisión eléctrica, etc., por tanto es
difícil tratar de aislar sólo una componente de interés del resto de la señal.

II.3.1 Interpretación Cualitativa

Los datos se presentan en mapas de configuración de PE en voltaje, ya que


facilita la inspección visual que puede revelar anomalías relacionadas con cuerpos
atractivos para el proyecto. La presentación de los datos en un mapa lleva un
proceso de acondicionamiento de la señal en donde se sustraen componentes no
19
deseados a fin de que el resultado sea el efecto debido a una anomalía.

Para el caso de datos de PE, el objetivo es buscar o al menos aproximar el


potencial residual de la señal debida a condiciones de oxido-reducción a la que
llamaremos potencial residual. Para tal tarea, se deben eliminar componentes
como la deriva del potencial en tiempo a causa de corrientes telúricas terrestres,
electrofiltración y difusión (Yüngül, 1950).

La falta de disponibilidad de información como datos de pozos de observación,


determinación de concentraciones y movilidad iónica de pares de sales
involucradas y cuantificación de potencial de doble capa, etc. es un factor que
limita la separación de aquellos fenómenos no deseados. Para este trabajo se
realizan aproximaciones a los datos para remover efectos como la deriva
temporal y de electrofiltración.

El PE medido sobre sitios contaminados está asociado a efectos que indican


condiciones favorables para la existencia y funcionamiento de procesos similares
al de una pila galvánica, una razón es por el rico contenido en materia orgánica y
otra es que generalmente la parte superior de un acuífero se encuentra en estado
de oxidación y al incorporarse fluidos con poder reductor modifican el ambiente
habitual de redox (Naudet et al. 2004). Según Nyquist y Corry (2002), también
existe correlación entre PE y zonas de metabolismo anaeróbico, en donde el suelo
actuá como una barrera permeable que permite el flujo de iones y desacelera la
mezcla directa entre la zonas anaeróbica y aeróbica (Figura 6).
20

Figura 6. Situación análoga en campo por condiciones redox,


metabolismo de materia orgánica.

II.3.1.1 Corrección por deriva temporal.

El primer paso consiste entonces en eliminar la componente de deriva temporal de


los datos de PE debidas a fluctuaciones del campo magnético terrestre que induce
fluctuaciones de corrientes eléctricas en la corteza y manto de la Tierra (Lowrie,
2007). Los efectos ionosféricos aparecen en el espectro de energía de campo
magnético como picos que representan distintos periodos que van desde
fracciones de segundo hasta años. También es posible observar que a mayor
amplitud de campo magnético, mayor amplitud también de campo eléctrico con
deflexiones muy similares (Figura 7).
21

Figura 7. Espectro de frecuencia de la intensidad


de campo magnético y su correspondiente espectro
de campo eléctrico inducido para un modelo
uniforme de resistividad de 20 Ohms-m.

Para tal objetivo, se tomaron datos de la actividad magnética terrestre durante los
días de medición del Geostationary Operational Environmental Satellite (GOES)
de la National Weather Services “NWS” de la NOAA, EEUU, que proveen
información del pronóstico del tiempo, seguimiento de tormentas severas e
investigación de meteorología. Satélites como el GOES-11 (GOES-oeste),
comúnmente localizado a 135° W sobre el Océano Pacífico y GOES-12 (GOES-
este), localizado a 75° W sobre el río Amazonas, tienen como una de sus
funciones monitorear el campo magnético terrestre mediante magnetómetros tipo
fluxgate de tres ejes.

La variación diurna, causada por la asimetría del campo magnético de la Tierra


22
con respecto a la dirección del Sol, es monitoreada con componentes llamadas Hp
(componente paralelo al eje de rotación de la Tierra) y Ht (Intensidad de campo
total), en donde valores cercanos o menores a cero y el satélite colocado en lado
diurno (mediodía) puede ser debido a una compresión de la magnetósfera. En
cambio, valores cercanos o menores a cero y posicionado de lado nocturno
(media noche) indica fuertes corrientes que se asocian a subtormentas e
intensificación de corrientes de la cola geomagnética de Tierra (Singer et al. 2008).

Figura 8. Variación de la componente Hp paralela al eje de rotación de la Tierra en los


días de medición y ajuste lineal de la deriva

En la Figura 8 es evidente observar en general, un aumento en la señal del campo


magnético en los días donde se tomaron datos de PE (recuadro en amarillo), por
tanto, el campo eléctrico también debió haber registrado un aumento en magnitud.

Considerando que las lecturas tienen efecto de variación temporal del campo
eléctrico (corrientes telúricas), por lo que se debe remover esta deriva a los datos
observados, de igual manera con un ajuste lineal ascendente (Figura 9) como la
encontrada en la variación del Hp (Figura 8). Ésta aproximación justifica una de las
razones de remover el llamado “efecto regional” o de ruido de fondo de los datos
23
de PE.

Deriva

Figura 9. Corrección por deriva temporal de observaciones de PE.

II.3.1.2 Corrección por Topografía.

La relación entre el PE y los cambios de elevación en el terreno, se debe


principalmente a la inducción de potenciales de filtración (Ernstson y Scherer,
1986) causados por la disposición de las capas geológicas que gobiernan el flujo
gravitacional del agua subterránea. Esto tiene relación con el nivel piezométrico, o
bien, con el grosor de la zona no-saturada, la cual tiende a incrementarse
proporcionalmente con la elevación (Finizola et al. 2004). Esta atribución es válida
cuando la roca está saturada por el fluido y éste tiene una conductividad eléctrica
muy alta en comparación con la de la matriz rocosa, donde la conductividad entre
24
los granos (interface sólido-líquido) es despreciable. El potencial de filtración
también se genera incluso en rocas parcialmente saturadas.

El potencial de filtración se da por interacciones iónicas en la interface sólido-


líquido del subsuelo (Figura 10), la descomposición de los minerales de la roca
generan una intensidad de campo eléctrico cuasi-estático el cual depende de la
constante dieléctrica, la viscosidad y la conductividad del fluido. La diferencia de
presión requerida para que suceda el potencial de filtración es obtenida por la
variación en el nivel de saturación del subsuelo.

Figura 10. Interacción iónica interfase mineral-solución y variación del potencial. Doble
capa: (1) Capa de Helmholtz. (2) Capa de Gouy-Chapman.

Se considera que el campo eléctrico es cuasi-estático debido a que el tiempo de


descomposición mineral y el tiempo de residencia del agua subterránea en el
punto muestreado es mucho mayor que el tiempo requerido para la toma de datos
de PE. Por tanto, para que esto sea válido los datos deberán tomarse en época de
estiaje, cuando la circulación es más lenta y se evita la componente vertical de
filtración y el flujo subsuperficial.
25
Se estima que los datos de PE tienen un efecto causado por la topografía, donde
el potencial incrementa cuando la altitud de la topografía disminuye y es más
próxima a la interfase sólido-líquido (Zlotnicki y Nishida, 2003), o bien, a la zona
vadosa y zona saturada, debido a que generalmente la distribución del nivel
freático varía de acuerdo a la forma de la topografía del terreno, aunque el relieve
del nivel freático no es tan abrupto como el relieve topográfico (Figura 11).

Figura 11. Distribución general del nivel freático

Para retirar el efecto topográfico de los datos obtenidos en el Exbasurero se utilizó


un Modelo Digital de Elevación con resolución de un noveno de segundo de arco
producido por el National Map Seamless Server del U.S. Geological Survey.
Posteriormente se obtuvo la elevación puntual de cada una de las estaciones en el
Exbasurero con el programa ArcGIS versión 9.2 y finalmente se gráficaron cada
uno de los perfiles de elevación y su correspondiente perfil de PE. En la Figura 12
es claro observar que conforme aumenta el nivel de elevación topográfico la señal
de PE disminuye.
26

Figura 12. Relación entre la elevación topográfica y la señal de PE.

A los datos PE tomados en el Exbasurero se le sustrajo una superficie plana


(Figura 13), ya que la contribución por efecto topográfico varía linealmente con la
elevación.

La sustracción de un plano a los datos de PE para eliminar el efecto por topografía


es razonable porque, como ya se mencionó, a 300 metros se encuentra el
afluente del arroyo Ensenada y a medida que el nivel topográfico se eleva
aumenta la profundidad al nivel freático. Por otra parte en un pozo cercano,
aguas arriba del arroyo, a una elevación de 153 m.s.n.m.m. el nivel fréatico
estático se encuentra a 4.1 m de profundidad (Apéndice III).
27

Electrodo Base
a)

PE (mV)

b)

Figura 13. a) Corrección del efecto topográfico de PE. b) Relieve topográfico y ubicación
del arroyo Ensenada.
28
II.3.2 Interpretación Cuantitativa

La interpretación cuantitativa de los datos de PE puede ser agrupada en dos


categorías: modelos numéricos en 2D y 3D en donde el subsuelo se discretiza en
celdas con propiedades físicas determinadas. Por otro lado, los métodos
geométricos fijos simples, que buscan la determinación y forma de la estructura
mediante modelos sencillos, cuerpos de forma regular y simple como cilindros,
esferas, barras y hojas con finita o semi-infinita longitud (Srigutomo et al. 2006).

La ventaja de los métodos geométricos fijos sobre los modelos numéricos


discretos es que no requieren información de la variación de densidad de
corriente, resistividad y estimación de la profundidad de la estructura en el terreno,
éstos sólo requieren los datos de PE residual y la interpretación de anomalías
aisladas es rápida.

Para este estudio, se utilizaron cuerpos geométricos simples como la esfera y el


cilindro para la interpretación cuantitativa de anomalías de PE residual con base
en el modelo de Yüngül, (1950) y Bhattacharya y Roy, (1981), en donde la
anomalía de PE puede ser aproximada a una polarización oblicua de la forma
inmersa en un medio homogéneo, mediante el método de análisis de derivadas
(Mohamed et al. 1998; Mohamed et al. 2003; Srigutomo et al. 2006).

Yüngül, (1950), formuló la modificación de la expresión de una anomalía de PE


producida por una estructura simple polarizada y sugirió remover el efecto
topográfico y el efecto regional para así obtener el PE residual.

Una anomalía de PE causada por un objeto simple puede ser representada como:

x i cos  z sen 
V  x i ; z , , q =k (5)
 x 2i z 2 q
29
donde: xi es la coordenada de posición horizontal, z, es la profundidad al centro
del cuerpo, Q es el ángulo de polarización, q es el factor de forma, 0.5, 1.0 y 1.5
para un cilindro vertical, cilindro horizontal y una esfera respectivamente. Cuando
q cercano a cero la estructura es una hoja horizontal; k es el momento dipolar
eléctrico.

Con el fin de reducir el número de variables del que depende el PE en la ec. 5 (k,
z, Q, q) y obtener los parámetros de un perfil de anomalía, se linealiza ésta
ecuación por medio de una análisis de derivadas, donde se toman 4 ventanas
arbitrarias móviles de longitud “s” a lo largo del perfil, se calcula el gradiente
horizontal de PE (dV/dx) en segunda derivada numérica sobre los datos como:

V  xi 2s −2V xi V  x i−2s


V xx  x i = (6)
4s2

donde V(xi) es el valor observado en el perfil. Por otra parte la segunda derivada
analítica de la ecuación 5 queda:

[ ]
2
k  xi 2scos z sen  2xi cos 2z sen   xi −2s cos z sen 
V xx  x i ; z , , s , q= −  (7)
4s 2 [ x i 2s2 z2 ]q x 2i z2  q [xi −2s2 z 2 ]q

sí xi=0 la ecuación anterior resulta:

V xx 0=
4s2 [
−2kz sen   4s2 z 2 q−z 2q
z 2q 4s 2 z 2 q ] (8)

donde Vxx(0) es el valor de la segunda derivada en el origen del cuerpo, que


corresponde con la posición del centro del cuerpo en superficie.
Ahora, tomando la ecuación 7 y 8, obtenemos la ecuación normalizada para x i=s y
xi=-s, que es la longitud de la ventana:

V xx  s z 2q 4sz 2q
= 2 2 q 2q 2
V xx 0 2z  4s z  − z  4s  z  [
4s cot 2z z
2 q
− 2q −
z
4s cot z
16s2 z 2 q ] (9)
30

V xx −s z 2q 4sz 2 q
=
V xx 0 2z  4s2 z 2 q−z 2q [
−4s cot 2z z
2 2 q
 4s  z 
− 2q 
z
4s cot  z
16s 2z 2q ] (10)

si sumamos las ecuaciones 9 y 10 por el lado izquierdo nos queda:

V xx  sV xx −s 
d 2= (11)
V xx 0
y por el lado derecho resulta:

z 2 4s 2z 2 q [ s 2 z 2 q − 9s2 z 2 q ]


f 2  s ; q , z = 2 2 q 2 2 q 2q (12)
 s  z  9s z  [ z − 4s2 z 2 q ]

donde:
d 2= f 2  s ; q , z  (13)

donde 2 subíndice es el orden de la derivada, y s la longitud de la ventana.

Se obtiene así el valor numérico de la segunda derivada de los datos sobre el


perfil de anomalía de PE (ec. 6), para cada valor de s en la ec. 11, que debe ser
igual (ec. 13) a f2(s; q,z) de la ec. 12, función dependiente ahora de sólo dos
variables (q, z) fácil de resolver.

En la Figura 14 se muestra un perfil de anomalía de PE causado por un cilindro


horizontal y la solución por análisis de segunda derivada, para valores de s de 2,
3, 4 y 5. Para cada valor de s obtendremos una línea, a cada el factor geométrico
q corresponde una profundidad z, la solución es la intersección de las líneas para
las 4 ventanas arbitrarias (s).
31
S o lu c ió n 2 a d e r q - z
4 .5

3 .5

p ro f u n d id a d z
3

2 .5

s = 2
2 s = 3
s = 4
s = 5
1 .5
0 .5 0 .6 0 .7 0 .8 0 .9 1 1 .1 1 .2 1 .3 1 .4 1 .5
F a c to r g e o m é t ric o q

Figura 14. Izquierda. Modelo sintético de PE. Derecha. Inversión del modelo sintético
obtenido por análisis de segunda derivada

Obtenidos el factor geométrico q y la profundidad z, el momento dipolar eléctrico


puede ser calculado por la relación de la ecuación 5, si xi=0 para todas las formas
q, despejando k queda:

V 0 z 2q−1
k= , sen≠0 (14)
sen 

donde V(0) es el valor de la anomalía en el origen.


Ahora, si asumimos que V(xi; z,Q,q)=0 de la ecuación 5, resulta:

−z
cot = (15)
x0

x0 es la distancia a la intersección con cero del perfil de la anomalía, para que sea
valida la asunción y es con esta ecuación con la que se obtiene el ángulo de
polarización.

Para determinar el origen de la anomalía en el perfil (xi=0) se traza una línea recta
32
uniendo el máximo y mínimo del perfil, la intersección con la curva con ésta recta
es el punto xi=0.

La manera en que se resolvió el problema para este trabajo fue mediante un


programa implementado en MatLab (Apéndice I) con el siguiente procedimiento:
– Se obtiene el valor de la segunda derivada numérica para cierta ventana s
(20,30,40,50).
– Posteriormente, usando la ecuación 13 y 12, para un conjunto de z
estimadas (1-100) se logran valores para q=0.5.
– Igualmente se adquieren valores de la ecuación 12 para q=1.5, que son los
extremos de las formas geométricas.
– Gráfica ambas curvas y se introducen los valores de arranque para aplicar
el método numérico de optimización por bisecciones.
– Terminado para todas las ventanas se grafican las curvas solución en
donde la intersección o región determina el q y z óptimo del perfil.

La orientación del perfil se muestra en la Figura 15:

- Perfil

Figura 15. Mapa de ubicación del perfil de anomalía de PE residual.


33
En la Figura 16 se muestra el del perfil de anomalía de PE del Exbasurero.
Después de varios ensayos se interpreta que este perfil se debe al efecto de dos
cuerpos, por lo que se obtuvieron los parámetros para dos formas simples
independientes. Posteriormente se combinaron en una zona de interacción de
frontera entre ellos, logrando un ajuste del 90%.
a)

b)

1 2

zona de interacción

Figura 16. a) Solución del perfil de anomalía de PE en el Exbasurero debido al efecto de


dos formas simples (q1=0.8 ≈cilindro horizontal y q2=0.5=cilindro vertical).b) Ajuste del
perfil a los efectos de las formas y zona de interacción considerada .
34
II.4 Resultados

Como resultados de la interpretación cuantitativa, se obtuvo el mapa de contornos


de PE residual (Figura 17) en donde es claro apreciar una zona hacia la parte sur-
oeste con polaridad positiva. Hacia la porción norte, a partir de los pozos de
respiración del Exbasurero, una anomalía con polaridad negativa que tiende a
prolongarse hacia el arroyo Ensenada. Ambas anomalías corresponden con los
extremos del cuerpo del Exbasurero, lo que sugiere que éste es el eje del dipolo,
que a su vez indican reacciones existentes al interior y fuera del mismo.

Figura 17. Mapa de contornos de PE residual. Resaltan valores positivos en la parte


suroeste y valores negativos en la parte norte a partir de los pozos de respiración.
35
Se intuye que la anomalía negativa en la porción norte del sitio es la ruta por
donde se dispersan los lixiviados, debido al cambio abrupto de la topografía
existente, que va de 176 a 132 m.s.n.m.m. en una línea de 200 m, marca un
camino preferente descendente de fluidos hacia el arroyo Ensenada. Por otra
parte también es consistente en comparación con observaciones hechas por
varios autores (Naudet et al. 2004; Nyquist y Corry, 2002; Revil et al. 2007),
quienes advierten claramente encontrar una anomalía de polaridad negativa
asociada a la extensión de una pluma contaminante, por reacciones de oxido-
reducción que ocurren en ésta.

Naudet et al. (2004), propone que debido al extremo poder reductor del lixiviado
dentro del sistema acuífero, la extensión del lixiviado modifica las condiciones
naturales de redox, pues éste es rico en materia orgánica y NH 4+ con gran
capacidad de ceder electrones durante la biodegradación provocado por bacterias.
Así entonces, las condiciones de redox marcan un gradiente de potencial de
reducción que establece el paso preferente de lixiviado en una pluma
contaminante y aumenta al alejarse del sitio (Figura 18).

Figura 18. Desarrollo de las zonas redox en una pluma


contaminante por lixiviados donde Eh es el potencial de
reducción Eh1<Eh2<Eh3<Eh0
36
La anomalía positiva hacia el sur del sitio corresponde con la parte más antigua
del Exbasurero, donde se sospecha existen reacciones anaeróbicas asociadas a
metanogénesis2.

En cuanto a los resultados de la interpretación cuantitativa se muestran en la


siguiente tabla para las dos formas previstas en el perfil de anomalía de PE:

Cuerpo x0 (m) q z (m) θ (grados) k (mV)

1 195 0.8 53.29 -60.7 1518.5

2 380.5 0.5 20.62 40.52 46.98

Tabla I. Parámetros obtenidos por descomposición de formas del análisis cuantitativo de


PE residual.

Con un error de ajuste del 10.16% (Figura 19).

1 2
zona de interacción

Figura 19. Ajuste de los datos de PE residual observado y el calculado.

2 Fase final de la descomposición orgánica donde el ácido acético, hidrógeno, agua, ácido fórmico y
metanol son convertidos por acción microbiana a metano y dióxido de carbono.
Capítulo III

Método eléctrico de corriente continua

Las estudios de resistividad proporcionan una imagen de la distribución de esta


propiedad intrínseca del subsuelo. Para convertir esta imagen en información de
utilidad geológica es necesario conocer cómo varía este parámetro en los
diferentes materiales que constituyen el subsuelo del área bajo estudio. En la
Tabla II se describen brevemente los valores de resistividad de las rocas más
comunes, algunos materiales de suelo y sustancias químicas.

Material Resistividad (Ω m) Material Resistividad (Ω m)


R. Ígneas y Met

Granito 5x103 – 106 Arcilla 1 – 100


Suelos y agua

Basalto 103 – 106 Aluvión 10 – 800

Mármol 102 – 2.5x 108 Agua Fresca 10 – 100

Arenisca 8 – 4 x103 Agua de mar 0.2


R. Sedimentarias

0.01 M A.
Químicos

Esquisto 20 – 2 x103 6.13


Acético

Limolita 50 – 4 x102 Xileno 6.998x1016

Tabla II. Resistividades de las rocas más comunes, algunos minerales y químicos

Como se puede apreciar en la Tabla II, las rocas ígneas y metamórficas


típicamente tienen altos valores de resistividad. La resistividad de estas rocas
38
depende mucho del grado de fracturación y del porcentaje del agua que las sature
(Orellana, 1982). Así, las rocas sedimentarias que comúnmente son más porosas
y tienen alto contenido de agua tienen valores más bajos de resistividad. La
resistividad del agua en los poros de la roca, tiene valores aún más bajos pues
depende de la concentración y tipos del sales disueltas, así como de la
temperatura.

Los valores de resistividad tienen un amplio rango de variación en las rocas de un


área de estudio, comparado con otras propiedades físicas obtenidas por métodos
geofísicos, como la densidad que usualmente cambia por un factor menor a 2 y la
velocidad sísmica que normalmente no varía más de un factor de 10 (López
Hidalgo et al. 2007) hacen que los métodos eléctricos resistivos sean técnicas
geofísicas muy versátiles.

La aplicación de los métodos resistivos particularmente de corriente continua,


Tomografía Eléctrica Resistiva (TER 2D) es diversa y se usa con éxito en:
– Geotécnia (profundidad de basamento, detección de fallas y cavidades).
– Hidrogeológia (profundidad del nivel freático, flujo preferente del agua, etc.).
– Medioambiente (detección de infiltraciones subterráneas, dispersión de
contaminantes, etc.).
– Arqueología (localización de restos de murallas, pavimentos y edificaciones
antiguas).

Para este trabajo se empleo TER 2D con el fin de identificar en el subsuelo del
Exbasurero de Ensenada, el flujo preferente de dispersión de lixiviados, evaluar su
extensión y corroborar los resultados obtenidos mediante Potencial Espontáneo.

Se llevaron a cabo 9 líneas de TER 2D que cubren el área del Exbasurero, con
longitudes que varían de 270 a 380 m, con un interespaciado electródico de 10 m.
39
En su mayoría, se utilizó configuración electródica Wenner y Schlumberger. Sólo
en la primera línea se usó Dipolo-dipolo. Estas líneas se realizaron en la
temporada de secas durante el mes de noviembre del 2009.

III.1 Teoría

El desarrollo de la TER 2D empezó con la aparición de los sistemas multi-


electródicos igualmente espaciados, en donde se tiene control del procedimiento
de adquisición de datos por medio de un computador que se encarga de alternar
el orden de las terminales a lo largo del conjunto electródico. Esto proporciona
gran eficiencia y bajo coste en las mediciones de campo (Sumanovac y
Dominkovic, 2007), por la rapidez y facilidad de toma de datos que permite
abarcar áreas considerables.

Básicamente, tiene el mismo principio físico que los Sondeos Eléctricos Verticales
(SEV) con la diferencia que se obtiene sensibilidad horizontal, porque al permutar
el orden de los electrodos se realiza un barrido que cubre lateralmente el
subsuelo.

Si se inyecta corriente en un punto A en la superficie de un medio homogéneo


conductor, el circuito tenderá a cerrarse, y la corriente regresará por un punto B.
Si en el medio no hay pérdidas ni fuentes de corriente, ésta es igual en el punto A
y en el punto B.

Para explicar el fenómeno partimos de las leyes de Maxwell para campos


electromagnéticos:
∇ x E=−∂ B/∂ t (16)

∇ x H =J ∂ D/∂t (17)
40
donde: E, es el campo eléctrico; B, es la densidad de flujo magnético (inducción
magnética); t, es el tiempo; H, es la intensidad de campo magnético; J, es la
densidad de corriente eléctrica y D, es la densidad de flujo eléctrico
(desplazamiento eléctrico).

Para condiciones estacionarias en donde no hay cambios temporales, las


ecuaciones 16 y 17, son igual a cero respectivamente, se entiende que el campo
eléctrico E y magnético B son conservativos o irrotacionales, sí se deriva E en un
potencial escalar V, se encuentra la siguiente relación entre ambos:

E=−∇ V (18)

a su vez el campo E debe satisfacer la ley de Ohm en su forma reducida para


medios isótropos:

J = E (19)

La ecuación 19 expresa que la densidad de corriente J en un punto, tiene la


misma dirección y sentido que el campo E en ese punto y es proporcional al factor
de σ. Este valor corresponde a la conductividad del subsuelo (inverso de la
resistividad “ρ”).

Considerando que en todos los puntos del semiespacio homogéneo son excitados
por el flujo de corriente se satisface la ecuación de continuidad:
∇⋅J =0 (20)

el medio es homogéneo y estacionario por tanto la ecuación 20 es igual a cero,


combinando la ecuación 18, 19 y 20 y desarrollando la ecuación resultante usando
propiedades vectoriales y dado que la conductividad es uniforme, resulta que el
potencial V cumple con la ecuación de Laplace:
41

∇⋅ E=∇ ⋅∇ V  ∇ 2 V ⇒ ∇ 2 V =0 (21)

que indica que en cada superficie equipotencial del medio encontraremos el


mismo valor de densidad de corriente J, dirigida radialmente respecto del centro,
excepto en los electrodos y en una superficie de discontinuidad resistiva (Figura
20).

Figura 20. Electrodo puntual de emsión de corriente en un medio


homogéneo de resistividad ρ

La solución de la ecuación 21 es del tipo V=c/r, donde c es una constante y r es la


distancia. Por la ecuación 18, para un medio homogéneo:

∂V
E=− (22)
∂r

evaluando el campo la distancia r


(23)
E=−
∂V
∂r
=−
∂ c c
∂r r 
= 2
r
42
para conocer la densidad de corriente J sobre la región, se realiza la integral de
superficie de la ecuación 20, en el volumen de la fuente (superficie semiesférica,
Figura 20) y aplicando el teorema de la divergencia de Gauss, resulta:

2  r 2 J =I (24)
sustituyendo E=c/r2 y despejando c queda:

I I
c= = (25)
2  2 

sustituyendo (25) en la solución de la ec. de Laplace:

I
V= (26)
2 r

esto es sólo aplicable a un punto fuente y un punto receptor, con distancia r.


Para determinar la resisitividad del medio, por tanto, es necesario conocer la
distribución del potencial eléctrico y la corriente suministrada. Para suministrar
corriente a través del espacio conductor comúnmente se usa un par de electrodos
(A,B) y un par de varillas (M,N) que determinan el potencial entre cada una de
ellas y su sustracción define la diferencia de potencial (Figura 21).

Para un arreglo tetraelectródico la diferencia de potencial se define como:

V=
[
I 1 1 1 1
− −  =
I
2  r1 r2 r3 r 4 K ] (27)

donde rn es la distancia entre los electrodos de corriente-potencial (Figura 21), K


contiene las distancias entre los electrodos, llamado “factor geométrico” y depende
de la distribución de los eléctrodos. Si despejamos ρ de la ecuación 27 finalmente
obtenemos:

V
= K
I (28)
43

Figura 21. Determinación de resistividad de un medio, arreglo tetraeléctrodico.


Líneas de corriente (azul) y de potenciales (negras discontinuas).

Hasta ahora se ha hablado de un medio homogéneo, pero en la naturaleza esta


condición rara vez es encontrada, puesto que el subsuelo varía de resistividad en
las tres dimensiones del espacio. Una primera aproximación es suponer que el
medio sólo cambia de resistividad en dirección horizontal, por tanto una medida
tomada no representa la resistividad real y se le llama resistividad aparente (ρa)

V
a =K (29)
I

donde: K, es una constante que depende de la geometría del arreglo electródico;


ΔV, es el voltaje medido; I,intensidad de corriente medida.

Los sondeos eléctricos de TER 2D nos ofrecen un modelo más preciso del
44
subsuelo, ya que consideran que la resistividad cambia tanto en la dirección
vertical como en horizontal a lo largo de la línea de medición y supone que la
resistividad en la dirección perpendicular a la línea de sondeo no varía, presunción
muy razonable cuando existen cuerpos o anomalías elongadas.

III.2 Metodología

Se conoce como “arreglo electródico” al conjunto de electrodos que inyectan


corriente en el subsuelo y registran el potencial generado. Generalmente, los
arreglos suelen ser de cuatro electrodos, aunque también suelen emplearse los
dispositivos de tres o de dos electrodos. Un par de varillas son los responsables
de inyectar corriente al subsuelo y otro par se encarga de registrar el potencial
eléctrico generado por la corriente inyectada.

Existen varias formas de disponer los electrodos en campo, se diferencian entre sí


por la distancia relativa entre ellos y por la posición de los electrodos de corriente
respecto a los de potencial. En la Figura 22 se muestran los diferentes dispositivos
empleados y la expresión de resistividad aparente ρa de cada arreglo.

La TER 2D toma como base estos arreglos y los repite definidamente de acuerdo
al número de eléctrodos y separación máxima electródica en función de la
profundidad y cobertura horizontal deseada.
45

 V V
 a a
a =2
a =2
I I

V
a = n n1 a  V
a = n n1 aI
I

V
a = n n1n2 a
I V
a = n n1n2 a
I

Figura 22. Arreglos electródicos utilizados en este trabajo.

Para este trabajo se usó un arreglo de 28 electrodos espaciados a 10 metros con


configuración electródica Wenner, Schlumberger y una línea de Dipolo-dipolo. Los
tipos de arreglos se programan y simulan en el software de ordenador propio del
equipo AGISSAdmin de Advanced Geosciences, Inc. para posteriormente
cargarlos al equipo y ejecutarlos en campo. En la Figura 23 se muestran las
pseudosecciones simuladas por el programa AGISSadmin para los arreglos
Wenner, Schlumberger y Dipolo-dipolo. Es importante señalar la distinta cantidad
de lecturas tomadas por cada tipo arreglo y la cobertura horizontal de datos sobre
46
el subsuelo en el perfil.

Figura 23. Simulación de cobertura horizontal de datos con


28 electros para arreglo electródico Wenner, Schlumberger y
Dipolo-dipolo, utilizados para este trabajo.

Para el arreglo Wenner se obtienen 117 lecturas con un factor de expansión


máxima de electrodos de a=9. Para Schlumberger se simulan 136 valores con una
relación AB/MN máxima de 11, es decir, que si a=1 entonces na=5 en unidades de
separación de electródica. Para Dipolo-dipolo se adquieren 237 mediciones con
separación máxima de n=11 y una separación máxima de dipolos de 3 unidades
de espaciamiento, es decir, a=3.
47
Conforme se mueve el orden de los electrodos a lo largo del perfil y su
espaciamiento entre ellos se obtienen niveles teóricos de investigación, que por
convención se grafican a 45 grados en el punto central del arreglo tetraelectródico
relativo (Figura 24) y es el sitio de atribución de la medición. Esta técnica permite
adquirir en campo una imagen distorsionada de la distribución de resistividad en el
subsuelo, llamada “pseudosección”, al unir por el método de interpolación curvas
de isoresistividad.

Figura 24. Niveles y puntos de atribución de la tomografía tipo Schlumberger en el


relevamiento electródico.
La selección de este tipo de arreglos electródicos para este trabajo consistió en
considerar los siguientes puntos:

– El tipo de estructura a ser estudiada. Se pretende hacer la detección de


lixiviados y su flujo preferencial de dispersión, que de existir, idealmente
podría asemejarse a una estructura de forma cónica en donde el centro es
el punto donde se originan los lixiviados, o bien, a una forma tubular
alargada, que razonablemente se aproxima a una estructura 2D (Zume et
48
al. 2006).

– La profundidad de investigación. En la Figura 25 se muestra la profundidad


media de investigación para cada uno de ellos. Esto da una idea del nivel
de exploración al que podemos sondear; los valores de profundidad son
determinados al integrar la función de sensibilidad con la profundidad
(Edwards, 1977). Cabe destacar que estas profundidades son válidas para
un modelo homogéneo de tierra, por lo que si hay grandes contrastes de
resistividad, la profundidad real de investigación se ve modificada.
En el caso de dos o más capas de diferente resistividad, teóricamente la
profundidad de penetración cambia porque las líneas de corriente que
atraviesan la frontera entre las diferentes capas requieren condiciones de
continuidad en la componente de densidad de corriente J normal a la
interfase y la componente de campo eléctrico E tangencial al linde.

– La sensibilidad de cada uno de los arreglos. Esta es una función que


básicamente nos muestra el grado en que un cambio en la resistividad de
una sección del subsuelo influirá en la medida del potencial. A valores altos
de sensibilidad mayor es la influencia de la región del subsuelo sobre la
medida (McGillivray y Oldenburg, 1990).
En la Figura 26, se muestra la función de sensibilidad para los arreglos
eléctrodicos Wenner, Schlumberger y Dipolo-dipolo, nótese que los valores
más altos de sensibilidad se encuentran cerca de los electrodos. Por otro
lado, la diferencia en la geometría de los contornos ayuda a determinar la
respuesta de los conjuntos a distintos tipos de estructuras (1D, 2D o 3D).
Para la configuración Wenner los contornos de sensibilidad son casi
horizontales abajo del centro del arreglo, lo que indica que es sensible a
cambios de resistividad vertical del subsuelo, no del mismo modo a
cambios horizontales.
El arreglo Schlumberger muestra los contornos de sensibilidad son
49
ligeramente diferentes al Wenner con contornos más delgados y alargados
en el centro del conjunto, sensibilidad inferior en las regiones entre A y M
(también N y B) y concentración mayor de sensibilidad entre A y B. Esto
significa que el arreglo Schlumberger es moderadamente sensible a
estructuras horizontales y verticales.
Respecto a la configuración Dipolo-dipolo los valores más importantes de
sensibilidad se ubican entre los electrodos de corriente así como también
en los de potencial, por lo que este arreglo es muy sensible a cambios de
resisitividad horizontal, pero relativamente insensible a cambios verticales.

Figura 25. Profundidades medias de investigación sobre un medio homogéneo obtenidas


para una distancia eléctrodica de 10m a los diferentes niveles de medición (Edwards,
1977).

– Relación señal-ruido. Está estrechamente relacionada con el valor de la


constante K de cada dispositivo; a mayor valor de esta constante menor
será la diferencia de potencial medido (ecuación 29) y viceversa para un
mismo modelo de subsuelo, es decir, la señal es inversamente proporcional
al factor geométrico usado para calcular el valor de resistividad aparente.
50

Figura 26. Función de sensibilidad para arreglo Wenner, Schlumberger (n=2) y


Dipolo-dipolo (n=4) correspondiente a un modelo homogéneo de resistividad,
(tomada de López Hidalgo et al. 2007).

En la Figura 27 se muestra una gráfica del decaimiento del voltaje V (señal)


a diferente separación de electrodos para los arreglos Wenner,
Schlumberger y Dipolo-dipolo para un modelo homogéneo.
Es importante señalar que Wenner es el arreglo con menor decaimiento en
la diferencia de potencial, por tanto se deduce que éste tiene mayor
fortaleza en la señal. En cambio, el Dipolo-dipolo el voltaje es inversamente
51
proporcional al cubo del factor “n”, lo que significa que para la misma
corriente, el voltaje medido disminuye más de 200 veces cuando “na”
aumenta de 1 a 6. Esto es una característica considerable si las líneas de
medición son llevadas a cabo en áreas con un ruido de fondo alto. Como
sucede en el sitio del Exbasurero donde PE superficial medido, apreciado
para la TER 2D como ruido, es de casi 100 mV.

Figura 27. Decaimiento de la señal para un modelo homogéneo de 100 Ohms- m, I=1.5A
y a=10m.
Se realizaron 9 líneas de TER 2D en el Exbasurero con la ubicación mostrada en
la Figura 28.
52

Figura 28. Mapa de ubicación de líneas de TER en el Exbasurero de Ensenada.

La configuración eléctrodica usada para este trabajo fue Wenner y Schlumberger


(Carpenter et al. 1991; Bernstone et al. 2000). Se realizó un sólo tendido en
arreglo Dipolo-dipolo en la línea 1 con fines de prueba pero con resultados no
satisfactorios y con pobre contribución a la construcción del modelo subterráneo.
También se realizó Polarización Inducida (PI) con resultados no satisfactorios, esto
debido a fenómenos como el fuerte acoplamiento de los electrodos (varillas de
metal) y el alto ruido de fondo en el sitio, que no permite la relajación armónica de
polarización en el tiempo.
El equipo se colocó en el centro del tendido de la línea, entre los electrodos 12 y
13, con el fin de evitar que la señal decreciera debido a la longitud del cable y para
realizar en las líneas 3 y 6 la técnica denominada “roll-along”, que consiste en
53
ubicar al principio del tendido cierto número de electrodos (12) para después éstos
pasarlos al final detrás de la línea de medición.

III.2.1 Instrumentación

El equipo usado para este trabajo fue un resistivimetro multielectrodo Super Sting
R1/IP de Advanced Geosciences, Inc. de 28 electrodos (Apéndice V), donde el R1
se refiere al uso de un sólo canal para 4 terminales por medida, es decir, 2 de
corriente y 2 de potencial únicamente. IP se refiere a que es capaz de medir
polarización inducida en el dominio del tiempo (cargabilidad).

La alternancia de los polos es automática de acuerdo al tipo de arreglo y a las


separaciones máximas programadas con el software AGISSAdmin. Además, se
puede simular la línea de medición con el software para: estimar el tiempo
requerido de cada línea, comunicarse con el equipo para cargar archivos de
comando (instrucciones de orden de electrodos en el arreglo), descargar los datos
obtenidos en campo hacia el ordenador y por último convertir los datos en otro
formato para distintos programas (2DINV y 3DINV, Surfer para Windows, Resix
Plus y Resix IP2D).

En la Tabla III se muestran las especificaciones técnicas importantes del equipo


utilizado para este trabajo Super Sting AGI.
54

Especificación Técnica Rango

Resolución de medición Max. 30 nV, dependiendo del nivel de ruido.

Corriente de salida 1mA- 2A

Voltaje de salida 400 V, dependiendo de la transmisión de


corriente y resistividad del suelo.

Potencia de salida 200 W

Impedancia de entrada >20 MΩ

Supresión de ruido Mejor que 100dB en f>20Hz

Fuente de alimentación 12V DC externa

Tabla III. Especificaciones técnicas del Super Sting AGI.

Las varillas usadas en campo para la adquisición de datos fueron de fierro de un


metro de longitud, también se usó solución salina para asegurar el buen contacto
entre el suelo y el electrodo.

III.3 Inversión 2D Tomografía Eléctrica Resistiva.

Se debe tomar en cuenta que las observaciones obtenidas en campo son datos
de resistividad aparente (ρa) y no corresponden a la resistividad real del subsuelo.
Si graficamos la pseudosección de estos valores se muestra efectivamente la
distribución de resistividad en el subsuelo pero con una visión deformada de
contornos y alteración en magnitud de valores, esto debido al tipo de arreglo
electródico usado y a la presencia de cuerpos de diferentes resistividades.

Para obtener la distribución de resistividad real de la estructura en estudio, los


55
datos se deben someter a un método de inversión, que consiste en obtener las
resistividades reales a partir de las resistividades aparentes medidas. La ecuación
diferencial que describe el comportamiento de la resistividad para medios
heterogéneos e isotrópicos es la ecuación de Laplace (ec. 21), que ahora es igual
a una constante y se puede resolver al menos con 3 tipos de métodos: Diferencias
Finitas, Elemento Finito y Ecuación Integral.

El método de Diferencias Finitas es rápido pero no contempla la topografía en el


modelo, para este trabajo se trataron los datos con el método de Elemento Finito,
que si contempla corrección topográfica y Ecuación Integral, que no contempla
topografía, pero se realizó una deformación de la malla utilizada a manera de
aproximación, con resultados aceptables.

III.3.1 Inversión 2D por el método de mínimos cuadrados con


suavizamiento constreñido.

Los datos tomados se trataron con el método de optimización Gauss- Newton con
suavizamiento constreñido mediante el software RES2DINV (Loke y Barker, 1996).
Este método es usado para determinar la resistividad de cada prisma (paramétros
del modelo) de un arreglo de bloques rectangulares y minimizar la diferencia entre
la resistividad aparente medida y calculada, de acuerdo a la siguiente ecuación:

 J T J uF  d = J T g (30)

donde: F=fxfxt+fzfzt y a su vez, fx es el filtro de achatamiento horizontal y f z es el filtro


de achatamiento vertical, el efecto de estos filtros es restringir la suavidad de las
resistividades del modelo a un valor casi constante, es decir, permite introducir una
matriz de peso o de ponderación generalmente se consideran los datos con el
mismo valor de peso y es usada especialmente para datos ruidosos; J, es la
56
matriz de derivadas parciales (Jacobiano); u, es el factor de amortiguamiento; d,
vector del modelo de perturbación, o bien, el vector de corrección de los
parámetros del modelo; g, vector de discrepancia que contiene las diferencias
logarítmicas entre los datos observados y calculados de resistividad.

La inversión por mínimos cuadrados consta de tres pasos principales:


– Calcula los valores de resistividad aparente para el modelo
– Calcula la matrix del Jacobiano (J) de derivadas parciales
– Resuelve el sistema de ecuaciones lineales de la ecuación 30
los valores de resistividad aparente y derivadas parciales son calculados usando
el método de Elementos Finitos (o bien Diferencias Finitas).

Al tratar de minimizar el cuadrado de los cambios espaciales, o rugosidad, en el


modelo de resistividad se tiende a producir modelos con variación suave de
resistividad, aproximación que es aceptable si la resistividad del subsuelo varia de
manera sutil y gradual. En algunos casos, el subsuelo geológico consiste en cierto
número de regiones que internamente pueden ser apreciadas como zonas
homogéneas con fronteras a los extremos que pueden o no corresponder con los
contornos del modelo de resistividad (“Principio de Supresión”).

III.3.2 Inversión 2D Multielectródica por Método de Ecuación Integral

Este método considera la relación directa entre la resistividad aparente y


verdadera como funciones ponderadas, con el programa RESIST6 Multielectrodo
CICESE. La Ecuación Integral (ecuación 31) descrita por (Gómez Treviño, 1987),
sugiere que cualquier medición electromagnética puede ser representada como un
promedio ponderado de la distribución de conductividad en el subsuelo
dv (31)
E r , r 0=−∫ G ' r , r ' ⋅E r ' , r 0 
V'
 r ' 
57
donde: G'(r,r') es la función de Green; E(r',r0), es el campo eléctrico en el subsuelo;
r0, es el vector de posición en el origen o punto de referencia; r' es el vector de
posición del punto interno del subsuelo; y r es el vector de posición en el punto de
medición.

A partir de la ecuación 31, se adapta al método de resistividad de corriente


continua (Pérez Flores, 1995; Pérez et al. 2001). La integral caracteriza un semi-
espacio con bajos contrastes de resistividad por aproximación lineal.

K −1
log a r A , r B , r M , r N =
2 V
∫ N  r A , r B , r M , r N , r '  log r ' dv (32)

donde: ρa, es la resistividad aparente y depende de la localización de la fuente


(A,B) y receptor (M,N); rA;B;M;N, es el vector de distancia hacia los electrodos de
inyección (A,B) y medición (M,N); K, es el factor geométrico del arreglo usado; N,
es el kernel que depende de la localización y características de la fuente receptora
y transmisora y del arreglo rectangular de prismas (discretización).

El método de inversión lo que hace es obtener el valor de resistividad de cada


prisma de la discretización del medio, mediante un sistema lineal del tipo y=Ax
donde: y es el conjunto de datos de resistividad aparente observado; A es la matriz
que describe el modelo y x es el vector de dimensión superior que tiene infinitas
combinaciones de valores, solución del sistema indeterminado.

Para encontrar la solución adecuada, que mantenga el compromiso entre el ajuste


a los datos y variación espacial de la solución se utiliza un método de
regularización que minimiza la norma cuadrática de los residuales de la forma:

1 2 1 2
F  x = ∥Y − Ax∥  ∥ Dx∥ (33)
2 2

donde: Y, son las observaciones (ρa); x, son las incógnitas; A, es la función de


58
peso correspondiente; D, es la matriz que contiene las derivadas parciales de
resistividad del subsuelo; β, factor de suavizamiento.

Figura 29. Consideración de distribución de electrodos de transmisión y recepción


para el método de Ecuación Integral.

La consideración de la posición de los electrodos para el método de la Ecuación


Integral, como vectores de distancia hacia los de corriente y potencial desde un
punto arbitrario (Figura 29), hace posible combinar varios tipos de arreglos para
caracterizar un mismo terreno. A dos distintos arreglos electródico, por ejemplo,
Wenner y Schlumberger, al aplicarle esta descomposición vectorial, obtendremos
sólo el log ρa (rA,rB,rM,rN) (ecuación 32) en función de su distancia a los electrodos
para cada observación y así, un conjunto de datos indiferentes a la configuración
usada, que caracterizan el subsuelo. Por tanto un modelo más robusto con las
ventajas de cada configuración en sensibilidad horizontal y vertical (Pérez Flores,
1995).
59
III.4 Resultados e Interpretación

La inversión de los datos primeramente se realizó con el programa RES2DINV


(Loke y Barker, 1996) por el método de mínimos cuadrados con suavizamiento
constreñido, de manera separada, para cada tipo de arreglo electródico.
Posteriormente, de forma múltiple los datos se invirtieron con el programa
RESIST6 Multielectrodo CICESE (Pérez Flores, 1995; Pérez et al. 2001).
Los modelos de resistividad se muestran para cada línea, las dimensiones en el
eje de las abscisas es en metros y en el eje de las ordenadas con respecto al
nivel medio del mar y en el Apéndice II el ajuste de los datos a cada modelo.

III.4.1 Línea 1

Se observa que los modelos obtenidos con la configuración electródica Wenner y


Schlumberger son similares, se reconocen principalmente 3 zonas (Figura 30):
1. Resistividades mayores a 250 Ohms·m hacia el extremo oeste cerca de la
superficie y que se profundiza hacia el este, con una protuberancia en el centro del
modelo. Los primeros electrodos, al este, fueron colocados en terreno natural por
lo que sugiere ser ésta la resistividad del gránito intemperizado ”grus” del sitio,
congruente con Gómez Treviño y Espinoza (1990).
2. Resistividades en la parte superficial que abarca el Exbasurero del orden del 0.3
a 250 Ohms·m se intuye que corresponde efectivamente al relleno de éste.
3. Hacia el extremo este se observa una anomalía que sugiere la presencia de un
cuerpo conductor y que en relación espacial se ubica a unos 50 metros
perpendicular a un respiradero, indicando una posible concentración de fluidos.
60

a)
W E

2 3

b)

2
3
1

c)
Línea 5

2 3

Figura 30.Modelos de la Línea 1. a,b) Resultados por el programa RES2DINV para el


arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6 Multielectrodo
CICESE.
61
III.4.2 Línea 2

La Línea 2 se realizó en modalidad Wenner y Schlumberger con modelos


independientes igualmente semejantes. Predominan principalmente dos zonas
(Figura 31):
1. Hacia el oeste del perfil donde se colocaron los primeros electrodos, no es
posible observar la misma tendencia en relación a la línea anterior, paralela a 80
metros de distancia, ni el basamento en la parte central, sólo hacia esta parte
parece ser más profundo el relleno del Exbasurero (modelo de Ecuación Integral).
Hacia el este, sugiere volver encontrar resistividades que corresponden con el
basamento del sitio.
2. En la longitud de 190 metros sobre el perfil 2, se encuentra un respiradero del
Exbasurero que corresponde con la zona más conductora del modelo,
confirmando la existencia y extensión de fluidos, se tienen pruebas de presencia
de fluidos al interior de los respiraderos, o bien, podría corresponder de zonas
húmedas.

III.4.3 Línea 3

La Línea 3 esta localizada hacia el extremo norte del sitio, corresponde a


relevamiento en “roll-along” de 380 metros de longitud en Wenner y Schlumberger,
con repetición de datos al centro del perfil.
Es de esperarse encontrar el basamento muy cerca de la superficie dado que la
Línea está ubicada hacia el norte a fueras del límite del Exbasurero. Se reconocen
2 zonas principales (Figura 32):
1. Hacia el extremo oeste se observa un cuerpo conductor alrededor de la posición
80 m, lo que sugieren la presencia de suelo con alto contenido de humedad, que
puede corresponder a concentración de fluidos contaminantes.
62

a)
W Respiradero E

b)

c)
Línea 5 Línea 7 Línea 8

Figura 31. Modelos de la Línea 2. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
63
2. Resistividades igual y superiores a 250 Ohms·m en la parte baja del modelo a
partir de los 120 m de altitud, que por control de campo sugieren corresponder al
granito intemperizado grus del sitio que se vuelve más resistivo a profundidad, lo
que es consistente con un basamento menos alterado en comparación con el
grus.

III.4.4 Línea 4

La Línea 4 se localiza en el extremo oeste del Exbasurero, con configuración


electródica Wenner y Schlumberger, con orientación sur-norte. Se reconocen
principalmente dos áreas (Figura 33):
1. Resistividades mayores a 500 Ohms·m, localizadas hacia el sur del perfil 4,
sugieren corresponder al granito sano que aflora en el terreno natural del sitio y
se extiende por todo el basamento del sitio. Es de importancia la flexión en la cota
a los 75 m del resistivo, corresponde con el control de campo de este afloramiento.
2. Dos cuerpos conductores hacia el extremo norte del perfil, de resistividad similar
a la anomalía oeste del perfil 3, nótese que estas líneas son a su vez
perpendiculares, indicando posiblemente concentración de una fase acuosa de
fluidos contaminantes.

III.4.5 Línea 5

La Línea 5 se localiza hacia el este, es paralela a la línea 4, intersecta a las Líneas


1, 2 y 6, abarca al sitio de norte-sur. Igualmente con modalidad electródica Wenner
y Schlumberger. Predominan principalmente 2 zonas (Figura 34):
1. Una zona de alta conductividad hacia el extremo norte y otra hacia el centro del
modelo, concerniente a una probable concentración de humedad y/o fluidos.
2. Hacia la parte extrema sur a profundidad se observan resistividades mayores a
250 Ohms·m, que son adecuadas para el granito a profundidad y grus en
superficie, congruente a su vez con el control de ubicación de los electrodos.
64

a)
W E

b)

c)

Línea 8
Línea 7

Figura 32. Modelos de la Línea 3. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
65

S N
a)

2 2
1

b)

2
2
1

c)

2
2
1

Figura 33. Modelos de la Línea 4. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
66
III.4.6 Línea 6

La Línea 6 se localiza al centro del sitio de estudio, va de oeste a este y es


paralela a la línea 2 a 40 metros al norte, intersecta a las Líneas 5, 7, 8 y 9,
realizada en configuración electródica Wenner y Schlumberger, este último en
modalidad roll-along de 380 m.
Aunque los modelos difieren en la localización del basamento, en la Línea 6 es
importante señalar tres zonas de interés (Figura 35):
1. Hacia el centro del perfil y hacia el este se identifican dos cuerpos de muy baja
resistividad <2 Ohms·m, que por observaciones de campo corresponden a puntos
donde existen los respiraderos del Exbasurero (Figura 28), vistos con mayor
nitidez en el modelo por método de Ecuación Integral, confirman la presencia de
fluidos acuosos contaminantes “lixiviados” al interior del subsuelo del Exbasurero.
2. Hacia el extremo oeste un pequeño cuerpo conductor de la misma magnitud
que los dos cuerpos anteriores, es también observado en la Línea 2. Este
conductor sugiere también la presencia de fluidos de lixiviado al interior del
Exbasurero, que aunque no se obtuvo muestra y es determinado por correlación a
partir de la zona 1.
3. Hacia la parte inferior del modelo, se observan resistividades mayores a 250
Ohms·m que corresponden al granito-grus de la zona.
Al extremo este no es posible identificar el basamento de granitico, pues el modelo
aparece distorsionado con resistividades de 8 Ohms·m. Observaciones de campo
en esta zona, indican que es una parte inestable del bloque del Exbasurero y se
intuye que posiblemente exista una capa de material con alta humedad. Por
características del programa de inversión no es posible identificar la interfase del
basamento al interior del sitio. Esto se debe a que donde existen cuerpos
altamente conductores, éstos tienden a extenderse hacia abajo (Pérez Flores,
1995), debido a que la ecuación de regularización busca la suavidad vertical en
el modelo.
67

a)

1 1
2

b)

1 1
2

c)
Línea 1
Línea 6 Línea 2

1 1

Figura 34. Modelos de la Línea 5. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
68

a)
W E

2
1 1

b)

2 1

c)
Línea 5 Línea 7
Respiraderos Línea 8 Línea 9

2 1 1

Figura 35. Modelos de la Línea 6. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
69
III.4.7 Línea 7

Esta Línea se ubicada cerca de la zona central del sitio, va de sur a norte y pasa a
5 metros de distancia de un pozo de respiración y a 35 metros de otro y fue
realizada con configuración electródica Wenner y Schlumberger de 270 metros.
Aunque no es posible diferenciar altas resistividades, que correspondan con el
basamento, hacia el extremo sur y norte, el control de campo indica que esos
electródos fueron colocados sobre el terreno natural. Se reconocen tres zonas de
interés (Figura 36).
1. Hacia el sur se observa un cuerpo de muy baja resistividad, que como antes se
mencionó, confirma la presencia y extensión de fluidos contaminantes en el
subsuelo del Exbasurero.
2. Cerca de la parte central del perfil, se tiene otro respiradero y su efecto es
observado en el modelo como la presencia de otro cuerpo conductor al interior del
terreno.
3. Hacia el extremo norte es notorio otro cuerpo altamente conductor, que
corresponde con el conductor localizado hacia las partes bajas del Exbasurero.
Aunque no es evidente la conexión entre ésta y la zona 2, sólo en el modelo
Schlumberger, es muy probable que sí exista, debido a que en el nivel 140
m.s.n.m.m. del terreno natural, el perfil se encuentra ya fuera del relleno de
basura, la pendiente es favorable y no existe evidencia de estructura alguna que
impida su paso. Puede corresponder con el camino preferente de fluidos
contaminantes en dirección hacia el arroyo Ensenada, donde la zona 2 es el foco
de lixiviados.
70

a)
S N

1
2
3

b)

2
3

c)

Respiradero Línea 2 Línea 6 Respiradero

Línea 3
1
2

Figura 36. Modelos de la Línea 7. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
71
III.4.8 Línea 8

La Línea 8 se ubica hacia la zona central oeste del sitio, va de sur a norte, fue
realizada en modalidad electródica Wenner y Schlumberger, a 10 metros de un
pozo de respiración del Exbasurero y es paralela a la Línea 7.
En el perfil 8 se reconocen básicamente tres zonas de importacia (Figura 37):
1. Hacia el extremo sur y central del perfil con resistividades mayores a 100
Ohms·m que sugieren corresponder a la localización del basamento de granítico.
2. Hacia la parte superficial sur con resistividades menores a 8 Ohms·m se
observa la presencia de un cuerpo de considerable extensión, que indica la
presencia de fluidos contaminantes lixiviados o bien zonas húmedas.
3. Hacia el extremo norte, indica la presencia de un cuerpo conductor extenso y
sucede la misma característica que en la Línea 7, la conexión entre esta zona y la
zona 2 está evidenciada sólo por el modelo Schlumberger y ligeramente en el
modelo de Ecuación Integral. Se intuye que en realidad existe una conexión entre
estas dos zonas, de forma análoga a la Línea 7.

III.4.9 Línea 9

Esta ubicada hacia el extremo oeste, va de sur a norte con configuración Wenner
y Schlumberger. En esta línea se localizan 3 zonas importantes (Figuras 38):
1. Al extremo sur un cuerpo conductor y otro al centro de alta conductividad con
resistividades menores a 8 Ohms·m, similar a la Línea 8, que indican la presencia
de fluidos contaminantes lixiviados o zonas húmedas.
2. Hacia el centro a profundidad con resistividades mayores se tiene indicio de la
localización del basamento de granítico.
3. Al extremo norte del perfil, coherente a las otras líneas paralelas, se ubica un
cuerpo de alta conductividad que se intuye es la extensión de fluidos
contaminantes provenientes del relleno del Exbasurero hacia las partes bajas
debido a la preferente vía que proporciona la pendiente del sitio.
72

a)
S N

2 3

b)

2
3

c)

Línea 2 Línea 6 Respiradero

Línea 3
2

1 3

Figura 37. Modelos de la Línea 8. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
73

a)
S N

1
1
3
2

b)

1 1
3
2

c) Línea 1 Línea 6

Línea 3

1
1
3
2

Figura 38. Modelos de la Línea 9. a,b) Resultados por el programa RES2DINV para el
arreglo Wenner y Schlumberger, respectivamente. c) Resultado por RESIST6
Multielectrodo CICESE.
74
III.4.10 Secciones de resistividad 3D

Finalmente, se presentan algunos cortes a diferentes profundidades de


investigación de la TER 2D en el Exbasurero con el fin de facilitar la interpretación
en la búsqueda de bajas resistividades que son correlacionadas con senderos
preferentes de infiltración de lixiviados hacia afueras del sitio (Figura 39).

-3.5m -8.5m

-16m -30m
75

-45m -62.5m

Figura 39. Cortes de resistividad del terreno del Exbasurero de Ensenada a -3.5, -8.5, -16,
-30, -45 y -62.5 metros de profundidad.

Como se puede apreciar existen cuerpos de muy baja resistividad hacia la parte
central del Exbasurero que predominan hasta una profundidad de -16 metros,
congruentes también con la ubicación de los respiraderos. Las observaciones de
campo evidencian que al interior de estos tubos de respiración existen fluidos, lo
que confirma que estas bajas resistividades, son atribuibles a concentraciones de
fluidos contaminantes lixiviados. También se observa respuestas de baja
resistividad con extensión alargada hacia el norte que a su vez posiblemente
indican su flujo preferente de migración pendiente abajo.

Hacia el extremo este de sitio sobresale siempre una zona de alta conductividad
que toca una parte interna del Exbasurero, corresponde a una parte donde no se
efectuaron líneas de TER 2D, por que se consideró como una zona no atractiva
dado que se ubica sobre el terreno natural en el borde del sitio. Ésta se muestra
de forma elongada al este, ya que los únicos datos para la interpolación de la
sección son las resistividades del modelo de la Línea 6.
Capítulo IV

Discusiones

Mediante la aplicación del método de PE y su interpretación cualitativa, es posible


identificar hacia la parte norte una anomalía de predominancia negativa (Figura
17) que indica el paso de lixiviados hacia afuera del Exbasurero (Nyquist y Corry,
2002; Naudet et al. 2003; Naudet et al. 2004) consistente a su vez con la
preferente pendiente topográfica (Figura 13).

El método de TER 2D demuestra que la corriente eléctrica es capaz de fluir a


través de la masa de residuos acorde a las observaciones de Carpenter et al.
(1991) y Bernstone et al. (2000).

En general, entre las variables que afectan la resistividad eléctrica del terreno la
principal es el contenido de humedad y la concentración de sales disueltas en los
fluidos que ocupan los poros y/o fisuras en la roca (Orellana, 1982), por tanto, la
resistividad del medio (ρ) se relaciona directamente con la resistividad eléctrica del
fluido (ρw) contenido en sus poros por la ecuación:
= F w (34)

donde: F es el factor de formación.

El factor de formación depende de la geometría de los poros y del grado de


saturación de los mismos. Se relaciona la resistividad ρ, la porosidad Φ y el factor
de formación F a través de la forma general de la ley de Archie:
=w a 
−m
S
−n
(35)
77
donde: m, es el exponente de la cementación de la roca (1.3-2.5); n, es el
exponente de saturación (~2); a, es una constante que depende de la tortuosidad
de los poros (0.5-2.5) y S, es el porcentaje de poros ocupados por agua.

Particularmente en un Exbasurero podemos además citar que la resistividad


medida depende también de la composición de los residuos, de la velocidad de
generación y movilidad de lixiviados, de la generación de gas, la temperatura
interna y de la variación de densidad de compactación.

La resistividad típica de material de residuos a granel es de 15 a 30 Ohms·m si


poseen cierto grado de humedad y de 30 a 70 Ohms·m si se encuentran secos y
con algún índice de compactación (Bernstone et al. 2000). Resistividades
observadas en los modelos de TER 2D del Exbasurero principalmente en la parte
superior.

La composición de los residuos que yacen en el sitio se estima que está


compuesta un 6.9% de vidrio, 5.3% metales, 11.7% plásticos, 19.6% celulósicos y
48.5% orgánicos (SEMARNAP e INE, 1999). De estos sobresale el contenido de
materia orgánica que favorece la proliferación de poblaciones de
microorganismos, que en la degradación de la fracción orgánica da como
resultado una fase acuosa de lixiviado y productos gaseosos conocidos como
biogás.

La velocidad de generación y movilidad de lixiviados está estrechamente


relacionada con la precipitación, evaporación y temperatura (Moreno et al. 2004).
De la climatología obtenida de 1984 a 2004, se tiene que Ensenada presenta una
precipitación promedio anual de 279.6 mm, que sucede en la temporada invernal
(noviembre-febrero), y una temperatura media anual de 17.8 °C.

Lo anterior significa que para un área de 6.6 hectáreas, considerando precipitación


78
constante y sin contar la evaporación, caen aproximadamente 18 450 metros
cúbicos de agua al año, disponibles para la formación de lixiviado y si la infiltración
es completa. Aproximación que es razonable pues no se tienen en cuenta los
fluidos residuales que se generan por la descomposición de la materia orgánica
que compensan parte que se pierde por evaporación en un clima templado. Ahora
bien, durante la lluvia la capacidad de infiltración del suelo del Exbasurero no
permite que el agua se encharque puesto que es altamente poroso por la
heterogeneidad de los material y la ligera compactación.

Por tanto, las bajas resistividades encontradas en el Exbasurero (menores a 15


Ohms·m) pertenecen a cuerpos con concentración de fluidos o a cuerpos con
cierto grado de saturación de humedad que corresponden a la ubicación de
lixiviados y es acorde con los resultados de los experimentos realizados por
Greller et al. (2005) con diferentes grados de compactación (ley de Archie) (Figura
40).

Prim e r grupo de m u e stras


S e gun do grupo de m u e stras
Te rce r grupo de m u e stras
Estado se co
Le y de Archie (a=1, m =2)
Le y de Archie (a=1, m =2.5)
Re sistividad (O h m s-m )

Estado de se cado por


grave dad

Lixiviado

Estado satu rado

Volu m e n con te nido de agu a (%)

Figura 40. Relación entre volumen contenido de agua y


resistividad (Greller et al. 2005).
79
Nótese en la Figura 39 hacia el norte del Exbasurero, una zona predominante de
muy baja resistividad hasta una profundidad de -16 metros que se extiende desde
interior del polígono a partir de los pozos de respiración, que es congruente
también con la área donde se encuentra la parte negativa de Potencial
Espontáneo (Figura 41), por lo que permite establecer que esta es la zona donde
ocurre la dispersión de lixiviados en el Exbasurero hacia el arroyo.

Figura 41. Comparación entre el mapa de contornos de PE residual y el modelo de


resistividad a una profundidad de 10m.
80
La evolución temporal del volumen de lixiviado producido de acuerdo al desarrollo
de disposición de residuos va creciendo en dimensiones, pues depende del
incremento de la superficie y disminuye posterior a la clausura (Figura 42)
(Farquhar, 1989; González Herrera, 2005). Si el Exbasurero fue clausurado de
manera oficial en 2004, entonces se encuentra en plena atenuación de producción
de lixiviados.

Figura 42. Izquierda. Evolución temporal de volumen de lixiviado en un sitio de


disposición (Farquhar, 1989). Derecha. Tendencia de la concentración de lixiviados con el
tiempo a partir de la clausura.

Sin embargo, la concentración de lixiviados en un tiradero, tiende a disminuir con


el tiempo a partir de su clausura. A 7 años del cierre del basurero de Ensenada se
considera que se encuentra en declive de producción de lixiviado debido a que la
degradación biológica y de contaminantes solubles se generan en un periodo
máximo de 5 años, persistiendo los contaminantes poco solubles que son de
mayor afectación al acuífero local.

En cuanto al perfil de interpretación cuantitativa de PE, que es interceptado por la


línea 1 de TER 2D al sur y la línea 3 al norte, se obtuvo un valor de profundidad al
81
centro (z) del primer cuerpo de anomalía de 53.29 metros con un factor
geométrico (q) de 0.8. Por parte de la TER 2D línea 1 (Figura 30) al oeste, a los 24
metros de profundidad (150 m.s.n.m.m.), se encuentra un cuerpo altamente
conductor, pudiendo ser éste la causa de la anomalía de PE. El segundo cuerpo
de anomalía de PE a 20.62 metros de profundidad al centro de éste y 0.5 de factor
geométrico, es interceptado por la línea 3 de TER 2D, donde al oeste a 10 metros
de profundidad (130 m.s.n.m.m.) sobresale un cuerpo de alta conductividad,
posible causa de la anomalía de PE.

La estimación de las profundidades de los dos cuerpos de anomalía de PE, que se


atribuyen por TER 2D son la causa de ésta, no son coherentes y representan
aproximadamente la mitad de lo evaluado en la interpretación cuantitativa de PE.

Sin embargo, una anomalía de PE residual es también interpretada como la


relación entre el grosor de una zona saturada y otra insaturada, dada una alta
razón de resistividad entre los substratos y homogeneidad en la zona insaturada
(Aubert y Atangana, 1996; Gurk y Frank, 2001). Esta aproximación, funciona para
un modelo estático de capas subhorizontales homogéneas similar al
comportamiento de un capacitor con distintos medios dieléctricos (Figura 43),
donde un cambio brusco de resistividad eléctrica entre la interfase de los medios
corresponde a la intensidad de campo eléctrico, por tanto al potencial eléctrico. La
profundidad estimada por la interpretación cuantitativa es, entonces, de interfase
entre aquella relación de las dos zonas y que es más congruente con TER 2D
(Figura 44).

La capa relleno en el Exbasurero contrasta significativamente por más de 10


veces en magnitud resistiva en comparación con el basamento granítico al que se
puede considerar homogéneo e insaturado, aunque su relieve no es horizontal.
Esta aproximación es razonablemente válida porque, además de que existe el
contraste resistivo requerido para el modelo, la capa relleno del Exbasurero se
82
encuentra con cierto grado de saturación de fluidos, acorde a la vez con este
modelo (Figura 43).

Figura 43. Modelo estático de capas subhorizontales de PE residual,


equivalencia de circuito eléctrico.

Por último, el ajuste de los modelos de TER 2D (Apéndice II) tienen en promedio
errores del 13.72%, que se consideran aceptables pues se reconoce que el
Exbasurero es una estructura con variaciones en resistividad en las tres
direcciones. Caracterizarlo con herramientas que tienen por principio teórico
variaciones sólo en dos dimensiones (2D) representa una aproximación válida que
proporciona información valiosa de la estructura del subsuelo en el sitio.

En el Apéndice II es posible observar la bondad del ajuste del modelo en cuanto la


magnitud de resistividad observada en campo y la calculada, valores quedan fuera
del alcance del modelo como las muy bajas y altas resistividades, debido a que:
– El programa de inversión esta diseñado para variaciones pequeñas de
resistividad lateral y vertical, por lo que contrastes fuertes quedan fuera del
alcance del método de inversión.
– Existen de cuerpos resistivos con variación 3D, que superficialmente
pueden corresponder a cuerpos vacíos dentro del relleno del Exbasurero, o
bien cuerpos caóticos de biogás, que se muestran en los datos observados
83
pero no son posibles de invertir por el programa dado su alto valor resistivo
y no continuidad lateral en el perfil.

Perfil L1
Error Absoluto= 6.168%

W PE E

170

160

150
q=0.8
≈cilindro horizontal Basamento
140

130

120
0 50 100 150 200 250

0 2 8 20 40 100 200 300 400 500 1000

Perfil L3
Error Absoluto= 17.20%
W E
170
160
PE
150
140
130 q=0.5 = cilindro vertical
120
110
Basamento
100

0 50 100 150 200 250 300 350

0 2 8 20 40 100 200 300 400 500 1000


Resistividad Ohms m

Figura 44. Comparación entre la interpretación cuantitativa de PE residual y TER.


Capítulo V

Conclusiones y recomendaciones

El método geofísico de PE resultó una técnica simple de fácil implementación en


campo y robusta, conveniente para la caracterización de infiltraciones
subterráneas, como el caso del Exbasurero de Ensenada. Su interpretación es
sencilla y proporciona información de zonas de anomalía eléctrica. Se considera
que es de gran utilidad para una inspección general, ya que auxilia a enfocar
áreas de interés, que después pueden ser cubiertas con métodos más
especializados.

La aproximación hecha para obtener el PE residual, resultó apropiada, ya que


permitió remover efectos no deseados para esta investigación y sirvió de base
para la adquirir parámetros en la interpretación cuantitativa de PE, dado que no se
contó con otro tipo de información para su obtención.

Mediante la interpretación cuantitativa de los datos de PE fue posible obtener, con


información de TER 2D, una buena estimación de la profundidad del basamento.
Que de acuerdo a la base teórica del modelo geométrico simple resultó
discrepante, ya que “z” debería corresponder al centro del cuerpo de anomalía de
PE y fué más congruente con la localización de la interfase entre el relleno del
Exbasurero y el basamento granítico, por lo que se sugiere hacer nuevas pruebas
de este método.

Los modelos de TER 2D obtenidos se consideran robustos, porque se obtuvieron


tanto por el programa RES2DINV, de manera separada para cada arreglo
85
electródico como de forma conjunta, por el método de Ecuación Integral RESIST6
Multielectrodo CICESE. Los modelos presentan grandes similitudes gráficas en al
menos 6 perfiles. Este último método ofrece cierta ventaja, ya que incorpora
características de cada tipo de arreglo en un sólo modelo, a diferencia del primero
que obtiene un modelo para cada configuración electródica, de un mismo terreno.

Las configuraciones electródicas usadas en el sitio, donde existe alto ruido de


fondo (1500 mV por PE), muestran ventaja las medidas hechas con los arreglos
Wenner y Schlumberger, que combinados logran capacidad de sensibilidad
horizontal y una mayor cobertura de datos que caracterizan el subsuelo. La
desventaja del arreglo Dipolo-dipolo, es que está sujeto a un rápido decaimiento
de la señal conforme se aumenta la distancia “na”.

Los resultados obtenidos muestran que existen infiltraciones de lixiviado del


Exbasurero de Ensenada hacia la parte norte y alcanzan unos 200m de extensión
hacia el arroyo Ensenada, localizado hacia las partes bajas topográficas. Los
resultados de PE residual cualitativo concuerdan con los resultados de TER 2D e
indican que la parte negativa de PE y la zona conductora de TER 2D, que va de
los 3.5 hasta los 20m de profundidad, se propagan dos extremidades (una al
oriente y otra al occidente) delimitan el camino preferente de dispersión (Figura
41).

Con este trabajo se establece una primera fase de evaluación del Exbasurero,
ubicación espacial y percance al sistema acuífero local, así como cuantificación de
los residuos que ahí yacen (Apendice IV). Hará falta elaborar:
– Sistema de pozos de monitoreo
– Sistema de captación de lixiviados
– Valoración de técnicas de remediación, con el fin de asegurar la
estabilización del sitio.
86

Bibliografía

Alatorre Zamora, M. A., Cortina, L., Campos Enríquez, J. O. y Birch Francis, S. 2008.
Interpretación de anomalías de potencial espontáneo aplicada a un vertero de desechos
urbanos en el occidente de México. e-Gnosis. 6(3): 1-24 p.

Almazan Vázquez, E. 1988. Marco Paleosedimentario y Geodinámino de la Formación


Alisitos en la Peninsula de la Baja California. Univ. Nal. Autón. México, Inst. de
Geología. Revista, v.7(1): 41-51 p.

Aubert, M. y Atangana Yéné, Q. 1996. Self-Potential Method in Hydrogeological


Exploration of Volcanic Areas. Ground Water. 34(6): 1010-1016 p.

Bernstone, C., Dahlin, T., Ohlsson, T. y Hogland W. 2000. DC-resistivity mapping of


internal landfill structures: two pre-excavation surveys. Environmental Geology. 39
(3-4): 360-371 p.

Bhattacharya, B.B., y Roy, N. 1981. A Note on the Use of a Nomogram for Self-Potential
Anomalies. Geophysical Prospecting. 29: 102-107 p.

Carpenter J.,P., Calkin F., S. y Kaufmann S., R. 1991. Assessing a fractured landfill cover
using electrical resistivity and seismic refraction techniques. Geophysics. 56(11):
1896-1904 p.

Comisión Estatal de Población (COESPO). 2009. Perfil Sociodemográfico de Ensenada.


Gobierno del Estado de Baja California. 02: 1067-1081 p.

Comisión Nacional del Agua (CNA). 2002. Determinación de la Disponibilidad del Agua
en el Acuífero Ensenada, Estado de Baja California, México. Subgerencia de
Evaluación y Modelación Hidrogeológica. 13 pp.

Cruz Aguirre, J. 2008. Vecinos es desacuerdo con el proyecto del ayuntamiento de


Ensenada: De basurero a panteón. Revista ZETA. Semana del 2 al 8 de octubre
(1853): 253 p.

Edwards, L.S. 1977. A modified pseudosection for resistivity and IP. Geophysics. 42 (5):
1020-1036 p.
87
Ernstson, K. y Scherer Ulrich, H. 1986. Self-potential variations with time and their
relation to hydrogeologic and meteorological parameters. Geophysics. 51(10): 1967
- 1972 p.

Farquhar G., J. 1989. Leachate: production and characterization. Canadian Journal of Civil
Engineering. 16: 317-325 p.

Finizola, A., François Lénat, J., Macedo, O. y Ramos, D. 2004. Fluid circulation and
structural discontinuities inside Misti volcano (Peru) inferred from self-potencial
measurements. Journal of Volcanology and Geothermal Research. 135: 343-360 p.
doi: 10.1016/j.jvolgeores.2004.03.009.

Gastil R., Gordon, R.P. Phillips y E.C., Allison. 1975. Reconnaissance geology of the state
of Baja California. The Geological Society of America, Inc. Mem. 140 p.

Gómez Treviño, E. 1987. Nonlinear integral equations for electromagnetic inverse


problems. Geophysics. 52(9): 1297-1302 p.

Gómez Treviño, E. y Espinoza, J. M. 1990. Estudio Geofísico en los Terrenos del Basurero
Municipal de Ensenada, B.C. Div. de Ciencias de la Tierra CICESE. 7 pp.

González Herrera, R. A. 2005. Evaluación del impacto y afectación hidrogeológica del


basurero de Mérida, Yucatán. Tesis Doctoral, UNAM. 290 pp.

Greller, S., Bouyé M., J., Guérin, R., Robain, H. y Skhiri, N. 2005. Electrical Resistivity
Tomography (ETR) applied to moisture measurements in bioreactor: principies, in
situ measurements and results. International Workshop Hydro-Physico-Mechanic of
Landfills. 5pp.

Gurk, M. y Frank, B. 2001. Cave detection using the Self-Potential-Surface (SPS)


technique on a karstic terrain in the Jura mountains (Switzerland). Kolloquium
Elektromagnetische Tiefenforschung. 9 pp

Institituto Nacional de Ecología; Secretaría de Medio Ambiente Recursos Naturales y


Pesca (INE-SEMARNAP). 1999. Reporte de Estado Ambiental y de los Recursos
Naturales en la Frontera Norte de México. INE. 256-285 p.

Loke H.M. y Barker D.R. 1996. Rapid least-square inversion of apparent resistivity
pseudosections by a quasi–Newton method. Geophysical Propecting. 44: 131-152 p.

López Hidalgo, A. A., Loke H. M.,Fanton O., G. y Cara Rubí ,E. 2007. Técnicas prácticas
para la investigación de resistividad en dos y tres dimensiones (tomografía
electrica 2D y 3D). Universidad Nacional de San Juan. 35 pp
88
Lowrie, W. 2007. Fundamentals of Geophysics. Cambridge University Press. New
York, USA. Capítulo 4.3, 256-276 p.

Madrigal, N. 2009. Prensa local El Vigía, Ensenada, B.C.; Será panteón el exbasurero:
Porfirio Vargas.

Manzella, A. 2006. Geophysical Methods in Geothermal Exploration. Italian National


Research Council. Pisa, Italy. 41 pp.

McGillivray R.,P. y Oldenburg W.,D. 1990. Methods for Calculating Fréchet Derivates and
Sensitivities for the Non-Linear Inverse Problem: A Comparative Study.
Geophysical Propecting.38: 499-524 p.

Migón, P. 2006. Granite Landscapes of the World. Oxford University Press Inc. New
York. 53-86 p.

Mohamed Abdelrahman, E., Barker Ammar, A.A., Ismail Hassanein, H. y


Abdelmottaleb Hafez, M. 1998. Derivate analysis of SP anomalies.
Geophysics .63(3): 890-897 p.

Mohamed Abdelrahman, E., M. El-Araby, H., G. Hassaneen, A.R. y A. Hafez, M. 2003.


New methods for shape and depth determinations form SP data. Geophysics. 68(4):
1202 - 1210 p. doi: 10.1190/1.1598112

Moreno, A., Gárfias, J., y Morales, P. 2004. Transporte de Lixiviado Proveniente de un


Basurero en un Acuífero Superficial Areno-Arcilloso, Metepec, Estado de México.
Congress XXXIII, International Association of Hidrogeologists. Zacatecas, Méx 11
-15 Octubre.

Naudet V., Lazzari, M., Perrone, A., Loperte, A., Piscitelli, S. y Lapenna, V. 2008.
Integrated geophysical and geomorphological approach to investigate the snowmelt-
triggered landslide of Bosco Piccolo village (Basilicata, southern Italy). Engineering
Geology. 34 pp. doi: 10.1016/j.enggeo.2008.02.008

Naudet, V., Revil, A. y Bottero Y.,J. 2003. Relationship between self-potential (SP) signal
and redox conditions in contaminated groundwater. Geophysical Research Letters.
30(21): 4 p. doi: 10.1029/2003GL018096

Naudet, V., Revil, A., Rizzo, E., Bottero Y., J. y Bégassat, P. 2004. Groundwater redox
conditions and conductivyty in a contaminant plume from geoelectrical
investigation. Hydrology and Earth System Science. 8(1):8 - 22 p.

Noguera Cáceres, J. F. 2003. Caracterización de zonas contaminadas por métodos


geoquímicos: Área de Gavá-Viladecans (Delta del Llobregat). Tesis de Doctorado
89
en Geología, Unidad de Petrología y Geoquímica. 229 pp.

Nyquist E.J. y Corry E. C. 2002. Self-potential: The ugly duckling of environmental


geophysics. The Leading Edge. 446 - 451 p.

Orellana, E. 1982. Prospección Geoeléctrica en corriente continua. Paraninfo. Madrid


España. 578 pp.

Pérez Flores, M. A. 1995. Inversión Rapida en 2-D de Datos de Resistividad,


Magnetoteluricos y Electromagnéticos de Fuente Controlada a Bajos Números de
Inducción. Tesis de Doctorado, CICESE. 352 pp.

Pérez Flores, M.A., Méndez Delgado, S. y Gómez Treviño, E. 2001. Imaging low-
frecuency and dc electromagnetic fields using a simple linear approximation.
Geophysics. 66(4): 1067 - 1081 p.

Raisz, E. 1964. Landforms of Mexico: Cambridge, Mass. U.S. Navy Geographical


Branch. Mapa, escala aprox. 1:3 000 000.

Revil, A., Titov, K., Doussan, C. y Lapenna, V. 2007. Applied Hydrogeophysics. Springer
Netherland. 383 pp.

Rozycki, A., Ruiz Fonticiella J. Ma., J. de la Cuadra, A.. 2006. Delimitación de las zonas
de flujo de agua con el método del potencial espontáneo en el espandón aguas
arriba de una presa. Ingeniería Civil 144. 121 - 126 p.

Sato, M. y Mooney ,M. H. 1960. The Electrochemical Mechanism od Sulfide Self-


Potentials. Geophysics. XXV(1): 226 - 249 p.

Secretaria del Medio Ambiente, Recursos Naturales y Pesca e Instituto Nacional de


Ecología. 1999. Reporte del Estado Ambiental y de los Recursos Naturales en la
Frontera Norte de México..INE, México D.F.: 349 pp.

Singer, H.J., Matherson, L., Grubb, R., Newman a. y Bouwer S.D. 2008. Monitoring space
weather with the GOES magnetometers. NOAA Space Environment Center. 11 pp.

Srigutomo, Wahyu, Eleonora, A. y Hablri Zen, M. 2006. Quantitative Analysis of Self-


Potential Anomaly: Derivative Analysis, Least-Square Method and Non- Linear
Inversion. Indonesian Journal of Physics. 17(2): 49 – 55 p.

Sumanovac, F. y Dominkovic Alavanja, S. 2007. Determination of resolution limits of


electrical tomography on the block model in a homogenous environment by means
of electrical modelling. Rudarsko-geolosko-naftni zbornik.19: 47 - 56 p.
90
Vázquez González, R. 2005 Informe ambiental antiguo basurero de la ciudad de Ensenada.
CICESE. Ensenada, B.C., reporte interno, 88 pp.

Vázquez González, R. 2007. Informe de actividades primer semestre 2007, Contaminación


por lixiviados en basureros: Caracterización, Bioremediación y Prevención.
CICESE. Ensenada B.C., reporte interno, 15 pp.

Yasukawa, K., Andan, A. Kusuma, S.D. Y Uchida, T. 2000. Self-Potential Survey in the
Mataloko Geothermal Prospect, Flores, Indonesia. Proceedings World Geothermal
Congress, May 28-June 10, Tohoku, Japan: 1985 - 1990 p.

Yüngül, S. 1950. Interpretation of spontaneous polarization anomalies caused by spherical


ore bodies. Geophysics.15: 237 - 246 p.

Zlotnicki, J. y Nishida, Y. 2003. Review on Morphological Insights of Self- Potential


Anomalies on Volcanoes. Surveys in Geophysics. 24: 291 - 338 p.

Zume T., J., Tarhule, A. y Christenson, S. 2006. Subsurface Imaging of an


Abandoned Solid Waste Landfill Site in Norman, Oklahoma. Ground Water
Monitoring & Remediation. 26(2):62 - 69 p.

Referencias electrónicas:

Geostationary Operational Environmental Satellite (GOES) de la National Weather


Services “NWS” de la NOAA, EEUU, consultado Octubre 2009:
http://www.swpc.noaa.gov/rt_plots/mag_3d.html

National Map Seamless Server del U.S. Geological Survey, consultado Octubre
2009:
http://seamless.usgs.gov/
91
Apéndices

Apéndice I

Datos de campo de Potencial Espontáneo con errores de cierre

Polígono Circuito Exterior (mV) Líneas Interiores (mV)


Est Sur Oeste Norte Este Norte Media Sur
1 20.44 10.23 -23.82 -3.86 -14.61 10.44 1.98
2 -7.86 -2.84 23.74 -24.22 5.44 -13.55 -5.52
3 6.65 -6.43 21.63 -19.27 -29.69 26.88 7.62
4 -20.65 -12.65 -2.46 6.52 -18.92 48.97 -6.58
5 -14.83 -1.38 -2.34 5.58 23.02 -46.42 12.45
6 13.32 16.06 7.21 4.45 -8.06 7.35 9.52
7 -7.66 0.5 -20.15 -29.78 -7.78 -18.82 3.85
8 -37.52 2.87 -26.78 42.3 33.76 10.43 -3.95
9 -24.25 -8.64 10.05 1.54 -27.63 7.47 4.52
10 41.35 -2.94 -56.41 15.96 -2.11 28.67 -1.23
11 15.56 0.95 28.31 3.19 -10.45 5.87 3.11
12 -5.65 27.24 44.22 8.87 11.93 -3.33 32.69
13 -13.84 -36.11 -10.52 31.41 -1.23 -8.16 -25.45
14 -2.11 12.26 14.88 -9.14 26.12 22.12 -23.25
Suma Suma Suma Suma E. de C. E de C. E. de C.
-37.05 -0.88 7.56 33.55 3.81 4.58 3.67
E. de C. 4.02
Tabla IV. Observaciones de campo de PE miliVolts (mV). Circuito Exterior se refiere al
polígono realizado en el Exbasurero con las respectivas direcciones. Líneas Interiores se
refiere a los perfiles al interior este polígono.(E de C.=Error de cierre)
92
Programa Matlab de Modelación-Inversión de PE residual.

%Programa de Modelación-Inversión para anomalías de Polarización Inducida


Residual con base en el método de interpretación geométrico fijo simple.
Mohamed et al. (1998), Mohamed et al. (2003) y Srigutomo et al. (2006)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all, close all
datos=load('prueba.dat'); % lee los datos de PE residual
x=datos(:,1);

%Abre ventanas de muestreo para la realizar la segunda derivada numérica


(Vxx)
for s=2:5;
close all
v=[zeros(2*s,1);datos(:,2);zeros(2*s,1)];
% Derivada Vxx
count=(2*s)+1; %acoplamiento de arranque para cualquier ventana
den=4*(s*s);
for i=count:length(datos)+2*s;
vxx(i-2*s)=(v(i+2*s)-2*v(i)+v(i-2*s))/den;
end
% Dibuja los datos y la segunda derivada numérica
plot(x,datos(:,2),x,vxx),grid
entr=[x,vxx'];

% Encuentra el valor de la segunda derivada en la longitud de la


ventana
s1=find (entr == s);
s2=find (entr == -s);
s3=find (entr == 0);

d2=(entr(s1,2)+entr(s2,2))/entr(s3,2)

% Para cada factor q encontramos z optima


q=[0.5:0.2:1.5];

% Propuesta de z iniciales arbitrarios


z=[0:10];
% Cálculo de la segunda derivada analítica para q mínimo
f2=f2s(z,0.5,s);
% Calcula la diferencia entre el valor de la segunda derivada numérica
y analítica.
dif=d2-f2;
% Cálculo de la segunda deriva analítica para q máximo
f2b=f2s(z,1.5,s);
difb=d2-f2b;
% Dibuja las diferencias entre las derivadas con el mínimo y máximo
valor de q, pide introducir valores de arranque para minimizar ésta
diferencia hasta la tolerancia, esto para diferentes ventanas s.
figure(2)
plot(z,dif,z,difb),grid on
title('Indroducir zup y zdown')
93
legend(['s=',num2str(s)])
hold on
zn=ginput(2);
zup=zn(1);
zdw=zn(2);

for k=1:length(q)
zu=zup;
zd=zdw;
tol=0.0001;
diff=d2-f2s(((zu+zd)/2),q(k),s);
while abs(diff)>tol
zop=(zu+zd)/2;
f22up=d2-f2s(zu,q(k),s);
diff=d2-f2s(zop,q(k),s);
if f22up*diff<0
zd=zop;
else
zu=zop;
end
end
zzz(k,s-1)=zop;
end
end

% Dibuja las curvas solución para cada ventana s


figure (3)
plot(q,zzz),grid
title('Solución 2a der q - z'), xlabel('Factor geométrico q'),
ylabel('profundidad z'), legend('s=2','s=3','s=4','s=5',4)
94
Apéndice II

Figura 45. Modelos observados y calculados de TER 2D. a) Para los arreglos
Wenner y Schlumberger, respectivamente por el programa RES2DINV. b) Gráfico
de ajuste de los datos al modelo con coeficientes de correlación de Pearson y
Spearman por el programa RESIST6 Multielectrodo CICESE.

a) Wenner b) Multielectrodo
Pearson =0.6002
Spearman =0.6864

Schlumberger
95
96
97
Apéndice III

Análisis químico básico de muestras de agua. De muestras tomadas sobre el


arroyo Ensenada, mapa de ubicación y de lixiviado del interior del Exbasurero.

Nombre X (m) Y (m) Alt.(m) n.F. (m) C.E.(mS/cm) S.D.T. (mg/L)


P DM 542325 3528107 153 4.1 2.21 1414.4
P DM’ 542299 3528125 153 4.3 1.91 1222.4
N OGlez 542116 3527988 147 0 2.82 1804.8
Zanja 541219 3527369 97.46 4.6 2.54 1625.6
Zanja 2 541235 3527396 97.44 4.6 2.79 1785.6
A Ens 1 541211 3527300 95 0 2.86 1830.4
A Ens 2 540852 3527217 86 0 6.19 4952
Tabla V. Ubicación, altitud en m.s.n.m.m., nivel freático estático (nF), conductividad
eléctrica (CE) y sólidos disueltos totales (SDT) de las muestras colectadas sobre el Arroyo
Ensenada.

Figura 46. Localización de muestras obtenidas en el arroyo Ensenada.


98

Muestra Mg2+ K+ Si+ Na+ Ca2+ F- Cl- S042-


PDM 65.64 2.52 20.39 168.3 76.95 0.23 456.4 135.81
N OGlez 86.1 6.46 18.81 210.6 98.04 0.25 680.11 121.62
Zanja 78.06 6.71 29.05 283.9 58.35 - 734.55 152.48
A Ens 182.1 8.72 18.77 441.2 152.2 - 1481.74 161.77

Lixiviado 28.4 99.74 26.66 1913.8 39.54 4.34 3352.68 55.62


Conductividad: 21.1 mS/cm Temperatura: 38.6 °C pH: 7.46
Tabla VI. Concentración de cationes y aniones principales en partes por millón (ppm) de
las muestras colectadas sobre el arroyo Ensenada y del lixiviado.
99
Apéndice IV
Cálculo de volumen de residuos del Exbasurero
El cálculo se realizó utilizando los datos de las líneas de TER 2D en la localización
del basamento. Se construyó una malla regular de 100 filas en “x” y 92 columnas
en “y” y mediante la técnicas numéricas se encontró el volumen diferencial entre
las mallas con la ayuda de los programa SURFER y VOXLER.

Figura 47. Topografía actual y anterior del Exbasurero de Ensenada.


Actualmente

Residuos sólidos Terreno natural

Anteriormente Estimación máxima:

Estimación media:

jo
Flu rente s
fe do
pre ixivia Estimación mínima:
l
de
100
Apéndice V

Memoria Fotográfica.

La hazaña Los residuos solidos

Potencial Espontáneo (PE) Super Sting R1/PI (TER)


101

Tendido de TER 2D Línea de TER 2D

Toma de datos de TER 2D Toma de muestra de agua

Toma de muestra en un respiradero El personal


del Exbasurero

También podría gustarte