Tesis Defendida Por: Ahiezer Alarcón Barradas
Tesis Defendida Por: Ahiezer Alarcón Barradas
Tesis Defendida Por: Ahiezer Alarcón Barradas
Dr. Marco Antonio Pérez Flores Dr. José Manuel Romo Jones
Miembro del Comité Miembro del Comité
24 de Septiembre de 2010.
CENTRO DE INVESTIGACIÓN CIENTÍFICA Y DE EDUCACIÓN SUPERIOR
DE ENSENADA
TESIS
que para cubrir parcialmente los requisitos necesarios para obtener el grado de
MAESTRO EN CIENCIAS
Presenta:
AHIEZER ALARCÓN BARRADAS
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.
A mis padres, Nieves Barradas Ramirez y Antonio Alarcón Alfonseca, por todo su
apoyo, fuerza y dedicación...
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 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),
Página
Resumen español……………………………………...……...…………... I
Resumen ingles…………………………………………………...……….. II
Dedicatorias………………………………………………………..……...... III
Agradecimientos…………………………………………………..……….. IV
Contenido…………………………………………………………..……….. V
Lista de Figuras…………………………………………………….…..….. VII
Lista de Tablas……………………………………………………….……... XII
Página
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
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
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
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
Figura Página
43 Modelo estático de capas subhorizontales de PE residual,
equivalencia de circuito eléctrico........................................... 82
Tabla Página
Introducción
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.
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).
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.
da
nsena
E
rroyo
A
b)
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
I.5 Antecedentes
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
II.1 Teoría
I a −I c
E d = RT ln C 1 /C 2 (1)
nF
E= p (3)
4
donde e,r y m son la contante dieléctrica (F/m), la resistividad ( Ohms-m)
II.2 Metodología
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
II.3 Interpretación
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.
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 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.
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
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:
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 2scos z sen 2xi cos 2z sen xi −2s cos z sen
V xx x i ; z , , s , q= − (7)
4s 2 [ x i 2s2 z2 ]q x 2i z2 q [xi −2s2 z 2 ]q
V xx 0=
4s2 [
−2kz sen 4s2 z 2 q−z 2q
z 2q 4s 2 z 2 q ] (8)
V xx s z 2q 4sz 2q
= 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 4sz 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 2z 2q ] (10)
V xx sV xx −s
d 2= (11)
V xx 0
y por el lado derecho resulta:
donde:
d 2= f 2 s ; q , z (13)
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
V 0 z 2q−1
k= , sen≠0 (14)
sen
−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.
- Perfil
b)
1 2
zona de interacción
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).
1 2
zona de interacción
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
0.01 M A.
Químicos
Tabla II. Resistividades de las rocas más comunes, algunos minerales y químicos
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
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.
∇ 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).
E=−∇ V (18)
J = E (19)
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)
∂V
E=− (22)
∂r
2 r 2 J =I (24)
sustituyendo E=c/r2 y despejando c queda:
I I
c= = (25)
2 2
I
V= (26)
2 r
V=
[
I 1 1 1 1
− − =
I
2 r1 r2 r3 r 4 K ] (27)
V
= K
I (28)
43
V
a =K (29)
I
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
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 n1 a V
a = n n1 aI
I
V
a = n n1n2 a
I V
a = n n1n2 a
I
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
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).
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.
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)
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)
1 2 1 2
F x = ∥Y − Ax∥ ∥ Dx∥ (33)
2 2
III.4.1 Línea 1
a)
W E
2 3
b)
2
3
1
c)
Línea 5
2 3
III.4.3 Línea 3
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
III.4.5 Línea 5
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
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)
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 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
-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
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)
Lixiviado
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.
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
Perfil L3
Error Absoluto= 17.20%
W E
170
160
PE
150
140
130 q=0.5 = cilindro vertical
120
110
Basamento
100
Conclusiones y recomendaciones
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.
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 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.
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. 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.
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.
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.
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.
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.
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.
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.
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.
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.
Referencias electrónicas:
National Map Seamless Server del U.S. Geological Survey, consultado Octubre
2009:
http://seamless.usgs.gov/
91
Apéndices
Apéndice I
d2=(entr(s1,2)+entr(s2,2))/entr(s3,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
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
Estimación media:
jo
Flu rente s
fe do
pre ixivia Estimación mínima:
l
de
100
Apéndice V
Memoria Fotográfica.