UNIVERSIDAD NACIONAL

AUTNOMA DE MXICO

FACULTAD DE INGENIERA

IMPLEMENTACIN DEL CDIGO ZONDRES2D

PARA LA MODELACIN DIRECTA E INVERSA
DE DATOS DE TOMOGRAFA DE RESISTIVIDAD
ELCTRICA 2D

TESIS PROFESIONAL
QUE PARA OBTENER EL TTULO DE:
INGENIERO GEOFSICO
P R E S E N T A:

ERIKA PAMELA DAZ ANAYA

DIRECTOR: M.C. Javier Urbieta Gmez

CIUDAD UNIVERSITARIA 2010

 Agradecimientos
A mi mama:

Por todo el cario, comprensin y cuidados que me has tenido, pero sobre todo por tu
eterno apoyo en todo lo que he decidido. Siempre has luchado por darme lo mejor en
muchos sentidos y por eso te admiro. Gracias por toda la felicidad que me has
brindado

A mi papa:

Por todos los conocimientos, consejos y enseanzas que me has compartido. Gracias
a ti tuve el coraje de hacer muchas cosas.

A mis hermanas:

Por siempre estar ah cuando las necesito, por llenarme de consejos y lecciones, por
compartir muchos momentos felices y por ser como unas madres para m.

A mis amigos:

Sol, Jos Luis, Hctor, Dani, Pao, Lilo, Leo, Gema, Iza En fin gracias a todos por
compartir diversin, alegras y preocupaciones, pero sobre todo por ofrecerme su
amistad, gracias a ustedes tendr siempre buenos recuerdos de esta etapa de mi
vida. Doy un agradecimiento especial al Muchacho, gracias a ti este trabajo fue un
poco ms sencillo, sin tu apoyo y ayuda no hubiera salido de muchos aprietos.

A mi tutor:

Gracias Javier por confiar en m para elaborar este proyecto, por todos los consejos,
conocimientos y apoyo que me brindaste y sobre todo gracias por toda la paciencia.

Al cuerpo de Sinodales:

Gracias por todo el apoyo y sugerencias para la realizacin de este trabajo.

A Octavio:

Gracias por compartir estos ltimos aos conmigo llenando mi vida de felicidad y
amor, por estar siempre ah apoyndome y alentndome para seguir adelante y ser
una mejor persona. Te amo!!
 NDICE Pgina

Captulo 1. FUNDAMENTOS TEORICOS DEL MTODO DE

CORRIENTE CONTINUA 1

1.1. Introduccin 1
1.2. Fundamentos de la Teora de Resistividad 1
1.2.1. Conductividad Elctrica 4
1.2.2. Concepto de Resistividad Aparente 5
1.3. Propiedades elctricas de los Materiales Terrestres 6
1.4. Clasificacin de las Tcnicas de Prospeccin Elctrica
de Corriente Continua. 10
1.4.1. Calicata Elctrica 10
1.4.2. Sondeo Elctrico Vertical (SEV) 11
[Link]. Efectos laterales en el SEV y
ambigedades en su interpretacin. 12
[Link]. Aplicaciones 12
1.4.3. Descripcin de la tcnica de Tomografa de
Resistividad Elctrica 2D 12
1.5. Dispositivos Elctricos utilizados en la tcnica de
Tomografa de Resistividad Elctrica 2D 14
1.5.1. Dispositivo Wenner-Schlumberger 14
1.5.2. Dispositivo Wenner 15
1.5.3. Dispositivo Polo-Dipolo 15
1.5.4. Dispositivo Dipolo-Dipolo 16
1.5.5. Dispositivo Polo-Polo 16
1.6. Equipo requerido en la tcnica de Tomografa de
Resistividad Elctrica 2D 17
1.7. Planificacin de un estudio de campo 18
1.8. Clculo de la Profundidad de Investigacin 19
1.9. Adquisin de datos de resistividad aparente 20
1.10. Procedimiento para la inversin de los datos 21
1.11. Ventajas y limitaciones del Mtodo 22

Captulo 2. ELPROBLEMA DE LA INVERSIN DE DATOS DE

RESISTIVIDAD APARENTE 24

2.1. Introduccin 24
2.2. El Problema Directo 24
2.3. Descripcin de un modelo 2D 24
2.4. Teora de Sensitividad 25
2.5. El problema inverso 27
2.6. Teora Bsica de Inversin 28
2.7. Planteamiento del Problema no lineal 28
2.8. Mtodos de discretizacin para un modelo 2D 34

Captulo 3. CDIGO DE INVERSIN DE DATOS DE

TOMOGRAFA DE RESISTIVIDAD ELCTRICA ZONDRES2D 37

3.1. Introduccin 37
3.2. Programa de Inversin ZONDRES2D 37
3.3. Caractersticas de ProfileR 39
3.3.1. Especificaciones del archivo 39
3.3.2. Detalles de [Link] 40
3.4. Implementacin del ProfilerV2 42
3.4.1. Creacin de un archivo tipo [Link] con ayuda
de ProfilerV2 43
3.5. Funcionamiento de ZONDRES2D 49

Captulo 4. APLICACIN EN MODELOS SINTTICOS 55

4.1. Introduccin 55
4.2. Modelado directo 2D 55
4.3. Ventajas y desventajas de las diferentes modalidades
de electrodos 56
4.3.1. Modalidad Wenner 56
4.3.2. Modalidad Wenner-Schlumberger 57
4.3.3. Modalidad Dipolo-Dipolo 58
4.3.4. Modalidad Polo-Polo 59
4.4.5. Modalidad Polo-Dipolo 60
4.4. Descripcin de los Modelos Sintticos 61
4.4.1. Modelo de un bloque rectangular 61
4.4.2. Modelo de dos bloques horizontales
rectangulares 63
4.4.3. Modelo de cuatro bloques horizontales
rectangulares 65
4.4.4. Modelo de dos capas y un dique 67
4.4.5. Modelo de una falla y un bloque horizontal 69

Captulo 5. APLICACIN EN DATOS DE CAMPO 72

5.1. Introduccin 72
5.2. Instrumentacin 72
5.3. Procesamiento de los datos 72
5.4. Delegacin lvaro Obregn 73
5.4.1. Marco Geolgico 74
5.4.2. Ubicacin Geogrfica 74
5.4.3. Perfiles de Imgenes Elctricas 75
 [Link]. Calle Pedro Aguirre (Colonia 75
Ampliacin Presiedentes)
[Link]. Calle Segunda Cerrada de Minas 76
(Colonia Arvide) 77
5.5. Coeneo, Michoacn 77
5.5.1. Marco Geolgico 78
5.5.2. Ubicacin Geogrfica 79
5.5.3. Descripcin del levantamiento 79
5.5.4. Resultados 81
5.6. Ixtln de los Hervores, Michoacn 81
5.6.1. Marco Geolgico 82
5.6.2. Ubicacin Geogrfica 82
5.6.3. Descripcin del levantamiento 84
5.6.4. Resultados 84
5.7. El Pahual, Espinal de Veracruz 84
5.7.1. Marco Geolgico 85
5.7.2. Ubicacin Geogrfica 86
5.7.3. Descripcin del levantamiento 86
5.7.4. Resultados 88
5.8. Valle Etla, Oaxaca 88
5.8.1. Marco Geolgico 88
5.8.2. Ubicacin Geogrfica 89
5.8.3. Descripcin del levantamiento 89
5.8.4. Resultados
92
Captulo 6. CONCLUSIONES Y RECOMENDACIONES
94
BIBLIOGRAFA
97
ANEXOS
 RESUMEN

En este trabajo se implementa una herramienta para el procesado e interpretacin

de datos elctricos de resistividad bidimensionales sintticos y reales con la ayuda
del programa de inversin, ZONDRES2D (Zond Geophysical Software) y del
cdigo de modelado directo e inverso, ProfileR versin 2.5 (Andrew Binley,
Universidad de Lancaster). Para la parte de modelado con el programa ProfileR,
se desarroll un cdigo en Matlab que genera archivos de datos de resistividad
con el formato apropiado que utiliza el ProfileR. Con ayuda de este cdigo se
generaron los diferentes modelos sintticos, cada uno de estos modelos se
construyeron utilizando cinco de las configuraciones de electrodos ms utilizadas:
Polo-Polo, Polo-Dipolo, Dipolo-Dipolo, Wenner y Wenner-Schlumberger. Con
respecto al modelado inverso, las secciones de resistividad 2D sintticas y reales
fueron configuradas para su interpretacin final con el programa de inversin
ZONDRES2D.

La aplicacin a datos reales con la tcnica de Tomografa Elctrica de Resistividad

2D se enfoc principalmente a estudios ambientales e ingenieriles someros en
cinco zonas diferentes: la primera, corresponde a la delegacin lvaro Obregn
(Ciudad de Mxico) con el propsito de mapear posibles cavidades; la segunda
pertenece a Coeneo, Michoacn; el estudio de esta zona se hizo con el fin de
delimitar reas con diferente geologa; la tercera; a Ixtln de los Hervores,
Michoacn, esta zona se estudi con el objetivo de caracterizar zonas con
potencial geotrmico; la cuarta; a El Pahualt (Municipio de Espinal, Veracruz) y la
quinta, a Valle Etla, Oaxaca. Estas dos ltimas zonas de estudio se realizaron con
la finalidad de detectar posibles zonas o formaciones susceptibles de explotarse
como abastecimiento de agua.

La adquisicin de los datos de campo se llev a cabo con dos diferentes equipos:
uno de la marca Scintrex, receptor IPR-10A y transmisor IPC-7 2.5 KW y el
segundo, con el sistema STING and SWIFT R1-IP (Advanced Geosciences Inc.).

Los datos de resistividad aparente obtenidos para las cinco zonas mencionadas
anteriormente, se procesaron para su interpretacin final mediante el programa de
inversin ZONDRES2D. Adems, para poder apreciar las bondades, alcances y
efectividad del ZONRES2D en cuanto a la resolucin de las imgenes elctricas
2D producidas con este programa, se comparan estas mismas imgenes con las
generadas por el programa de inversin RES2DINV el cual es un programa
comercial, comnmente utilizado en el procesamiento de datos elctricos y
reconocido por la comunidad Geofsica mundial. Los resultados de las secciones
de resistividad 2D procesadas con el programa ZONRES2D son aceptables.
 INTRODUCCIN

La Tomografa Elctrica de Resistividad es una tcnica de exploracin geofsica

que genera imgenes de la distribucin de la resistividad de una superficie o
volumen, basada en la medicin de numerosos potenciales elctricos que
aparecen debido a la aplicacin de una corriente estacionaria dentro del mismo
volumen. Las imgenes creadas son superficies bidimensionales o volmenes
tridimensionales que reflejan contrastes de resistividad existentes dentro del medio
los cuales se pueden asociar a estructuras o anomalas de inters geolgico,
ambiental o geotcnico.
A principios de 1980, la aproximacin para la inversin de datos de resistividad 2D
fue desarrollada, e.g. los trabajos de Inman (1975) y Tripp et al. (1984). El artculo
de Loke y Baker (1996) fue la base del programa de inversin de resistividad
llamado RES2DINV, comercialmente disponible, el cual ha sido ampliamente
probado en la prctica. A partir de este programa surgieron muchos otros basados
en principios similares, tales como: EarthImager 2D (Advanced Geosciences Inc.),
ResistivityImager 2D (Geotomographie), SensInv2D (Geotomographie), Eletom 32
(Geo&soft international), Resix 2DI (Interpex) y TS2DIP (Zonge Engineering
Research Organization), entre otros. Todos ellos son programas comerciales
utilizados en los trabajos de resistividad elctrica que se realizan por todo el
mundo, pero requieren llave fsica tipo USB y cuestan bastante caros.
Cada vez se requieren ms trabajos de Geofsica empleando esta tcnica, por lo
tanto, se tiene que hacer uso de algn programa de inversin 2D para procesar los
datos.
Por lo anterior, el objetivo principal del presente trabajo es implementar las
herramientas para procesar e interpretar datos elctricos de resistividad
bidimensionales sintticos y reales con el cdigo para modelado directo e
inversin de datos Profiler versin 2.5, desarrollado por Andrew Binley de la
Universidad de Lancaster; el desarrollo del cdigo Profiler V2 de captura de datos,
el cual crea archivos con el formato que utiliza el Profiler para realizar el modelado
directo e inverso y por ltimo, con el programa de inversin ZONDRES2D, escrito
por Alexander Kaminsky de la compaa ZOND GEOPHYSICAL SOFTWARE. El
programa de inversin ZONDRES2D se prob en todas las plataformas de
Windows funcionando adecuadamente y adems no requiere de llave fsica, lo que
lo hace una excelente opcin para utilizarse en programas acadmicos.
En el primer captulo se describen los conceptos fundamentales de los mtodos de
corriente continua, como definiciones tericas, objetivos, las diferentes tcnicas de
corriente continua, las configuraciones de electrodos para la adquisin de datos
de Tomografa de Resistividad Elctrica y las ventajas y limitaciones del mtodo.
En el segundo captulo se explican los principios bsicos de la teora de la
modelacin para la inversin de datos de resistividad. El anlisis cuantitativo de
los datos obtenidos en el campo se realiza resolviendo el problema directo y/o el
problema inverso. El problema directo predice resultados cuando la geologa y la
fsica del medio son conocidas, a partir de la distribucin de las resistividades del
subsuelo se calcula la resistividad aparente en superficie para una configuracin
en particular. De esta forma es posible comparar los datos de campo con modelos
conocidos y as verificar cuantitativamente si estos datos provienen de una
distribucin de resistividades similar a la del modelo propuesto. El problema
inverso construye un modelo razonable del subsuelo a partir de los datos
obtenidos en el campo. Es decir, a partir de los valores de resistividad aparente
medidos en superficie se obtiene la distribucin de resistividades del subsuelo.
En el captulo 3 se describe el funcionamiento y utilidades del programa para
modelado e inversin de datos Profiler. Tambin se describe el programa que se
implement en este trabajo para el modelado directo e inverso, Profiler V2, el cual
se diseo con la finalidad de crear archivos que sean compatibles con el programa
Profiler y hacer ms amigable el uso del mismo. De igual forma se describen las
caractersticas y funcionamiento del programa de inversin ZONDRES2D.
En el cuarto captulo se analiza el comportamiento de cuatro modelos sintticos
que asemejan anomalas de resistividad en el subsuelo con la finalidad de hacer
una recreacin de estructuras geolgicas y estudiar el efecto que tendran con un
estudio de Tomografa de Resistividad Elctrica. Para lograr esto, cada uno de los
modelos se calculo con cinco arreglos de electrodos diferentes y sus respectivas
configuraciones en superficie para la obtencin de la Tomografa. Adems, se
realiz un anlisis de sensitividad para cada uno de estos arreglos con sus
respectivas configuraciones con el propsito de comprobar cul de estas
configuraciones es la ms adecuada basndonos en su respuesta de sensitividad.
Para una fcil comparacin de las imgenes obtenidas con los diferentes arreglos
de electrodos se utiliz el programa de inversin, ZONDRES2D.
En el captulo cinco se presentan estudios de aplicacin a datos reales con la
tcnica de Tomografa de Resistividad Elctrica 2D. Los datos de resistividad
aparente se procesaron con los programas de inversin ZONDRES2D y
RES2DINV. Posteriormente, se realiz la interpretacin final de las imgenes
elctricas 2D generadas con ambos programas.
Finalmente, se presentan las conclusiones derivadas de este trabajo, as como
algunas recomendaciones para estudios futuros.
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

CAPTULO 1. FUNDAMENTOS TERICOS DE LOS MTODOS DE

CORRIENTE CONTINUA.

1.1. INTRODUCCIN

En este captulo se presentan las bases tericas de los mtodos elctricos con la
finalidad de determinar la distribucin de la resistividad del subsuelo haciendo
mediciones sobre la superficie. A partir de estas mediciones, se puede estimar la
resistividad verdadera del subsuelo, dado que la resistividad est relacionada con
diversos parmetros de inters geolgico tales como tipo de roca y grado de
saturacin de agua. Tambin, en este captulo se describen los mtodos elctricos
de exploracin. Estos mtodos investigan potenciales elctricos naturales y
propiedades elctricas intrnsecas de los materiales como son: conductividad,
permeabilidad y permitividad.

El mtodo de resistividad consiste en introducir una corriente elctrica continua en

el terreno por un par de electrodos y posteriormente medir la diferencia de
potencial entre otro par de electrodos. La relacin de estas dos magnitudes es lo
que se conoce como resistividad aparente.

Este mtodo puede aplicarse a la deteccin de reservas geotrmicas, localizacin

de restos arqueolgicos y exploracin de agua subterrnea. Su aplicacin es
menor en la exploracin minera y petrolera. Destaca el bajo precio de los equipos
empleados y su fcil implementacin, si bien requiere de un procedimiento
laborioso (desplazamiento de los electrodos y cables).

1.2. FUNDAMENTOS DE LA TEORIA DE RESISTIVIDAD

El comportamiento fsico de las rocas depende de su contenido de agua y modo

de agregacin de sus minerales, forma, volumen, porosidad, temperatura y
presin.

La Resistividad es una medida de la dificultad de la corriente elctrica a su paso

en un material determinado. La resistencia de un conductor es independiente de la
corriente y del voltaje.

Las medidas de resistividad elctrica del subsuelo son habituales en las

prospecciones geofsicas. Su finalidad es detectar y localizar cuerpos y estructuras
geolgicas basndose en su contraste resistivo. La magnitud de esta medida
depende, entre otras variables, de la distribucin de resistividades de las
estructuras del subsuelo, de las distancias entre los electrodos y de la corriente
inyectada.

El mtodo de Resistividad considera un cuerpo semi-infinito con una resistividad

uniforme, en donde se colocan cuatro electrodos, dos de corriente y dos de
potencial. El gradiente de potencial es medido por un voltmetro o potencimetro.
1
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Este mtodo es una aplicacin prctica de la Ley de Ohm.

Figura 1.1. La Figura muestra el principio de medicin de la resistividad del suelo: se inyecta una
corriente I entre el par de electrodos AB y se mide el potencial V entre el par de electrodos MN.

Las ecuaciones matemticas que explican el fenmeno fsico, son las ecuaciones
de Maxwell (Ec.1.1), que rigen el comportamiento de los campos elctricos y
magnticos

=

(1.1)

= +

Para los regmenes estacionarios no existen cambios temporales por lo tanto las
ecuaciones anteriores (Ec.1.1) se expresan de la siguiente manera:

r =0 (1.2)
= (1.3)

En los mtodos elctricos de corriente continua la ecuacin que es de inters es la

ecuacin 1.2. Si se considera una fuente puntual, es decir un electrodo de
corriente, las ecuaciones que rigen al campo elctrico son:

r =0 (1.4)

2
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

r = I(r) (1.5)

Suponiendo

= () (1.6)

Sustituyendo la ecuacin 1.6 en 1.5, tenemos

2 = () (1.7)

En una regin donde no existe fuente el potencial satisface

2 = 0 (1.8)

Que es la ecuacin de Laplace donde tambin la divergencia de es cero.

La ecuacin de Laplace ser vlida en todo el semiespacio conductor, excepto en

los electrodos.

= 2 2 (1.9)

= (1.10)

Sustituyendo ecuacin 1.9 en ecuacin 1.10, e integrando se obtiene

= 2 . (1.11)

Por lo tanto la diferencia de potencial entre los puntos M y N est indicada por
1 1 1 1
= 2 + (1.12)

Despejando la resistividad

= (1.13)

Donde

2
= 1 1 1 1 (1.14)
+

Esta ltima ecuacin se le conoce como el factor geomtrico del potencial medido
entre MN, y slo depende de la posicin relativa de los cuatro electrodos.

3
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.2.1. Conductividad Elctrica

La conductividad elctrica es la capacidad de algunos materiales para transmitir

corriente elctrica cuando es aplicado un voltaje. Su inverso es la resistividad
elctrica.

En rocas y minerales la conductividad elctrica puede darse de las siguientes

formas:

1. Electrnica. Es la forma normal de flujo de corriente en materiales que

poseen electrones libres, como en el caso de los metales.

2. Electroltica. Cuando el agua presente en los poros tiene alguna sal en

solucin, los iones constituyentes se separan, por lo que son libres para
moverse independientemente en la solucin.

3. Conduccin Dielctrica. Se lleva a cabo en conductores pobres o aislados

bajo la presencia de un campo elctrico variable.

La mayora de las rocas no son buenas conductoras de corriente elctrica, son

prcticamente aislantes debido a que los minerales formadores de roca tambin lo
son, sin embargo, la presencia de poros y fracturas tiene como resultado que los
valores de conductividad se incrementen. Generalmente la conduccin elctrica de
una roca se lleva a cabo por procesos electrolticos debido a que poros y fracturas
contienen algn porcentaje de agua. La conductividad depender de la cantidad
de agua presente as como de su salinidad y distribucin en la roca almacn.
Cuando en la conduccin electroltica se aplica un campo elctrico los iones
positivos (cationes) son atrados hacia el polo negativo, mientras que los iones
negativos (aniones) se ven atrados hacia el polo positivo. Por lo tanto, en una
roca la conductividad depender de la permeabilidad, lo que incluye el nmero de
poros, dimensin e interconexiones, el grado de humedad, la concentracin de
electrolitos disueltos, as como la temperatura y fase en que se encuentra el agua
de poro.

Existe una relacin indirecta entre la conductividad y la edad geolgica de una

roca debido a que ambas estn relacionadas con la porosidad (Tabla 1, Keller y
Frischknecht, 1966). Las rocas ms antiguas son por lo general menos porosas ya
que se encuentran a una mayor profundidad, lo que ocasiona que por presin los
espacios porosos sean cerrados.

Sin embargo, pueden darse excepciones por factores como la presencia de fallas
en las rocas que normalmente no son porosas o altas concentraciones de
minerales conductores, lo que en ambos casos podra dar un incremento en la
conductividad elctrica.

4
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Tabla 1. Relacin de resistividades elctricas de acuerdo con la litologa y edad geolgica

(Keller y Frischknecht, 1966).

Edad geolgica Arena marina, Arenas Rocas volcnicas Granito, gabro Caliza, doloma,
pizarra, terrestres, como basaltos, ahhidrita, sal.
grauwaca argilita, arcosa riolitas, tobas
Cenozoico 1-10 15-50 10-200 500-2000 50-5000
Mesozoico 5-20 25-100 20-500 500-2000 100-10,000
Carbonfero 10-40 50-300 50-1000 1000-5000 200-100,000
Paleozoico 40- 200 100-500 100-2000 1000-5000 10,000-100,000
Precmbrico 100- 2000 300-5000 200-5000 5000-20,000 10,000-100,000
1
Valores de resistividad en ohm-m.

1.2.2. Concepto de Resistividad Aparente

Un semiespacio homogneo de resistividad rara vez se da en la naturaleza ya

que el subsuelo no consta de capas homogneas y que la homogeneidad lateral
es difcil de esperar, por lo tanto en la prospeccin geofsica se trabaja bajo el
concepto de resistividad aparente a .

Definimos a la resistividad aparente en un semiespacio homogneo con

heterogeneidades donde las lecturas cambiaran de una posicin a otra obteniendo
una diferente a la del semiespacio homogneo. Esta variable no es el promedio
de los valores de obtenidos del subsuelo, sino el valor de medicin individual
que se asocia al punto inmediatamente debajo del centro del arreglo de electrodos
y es funcin del espaciamiento de los electrodos y de la distribucin real de la .
Es decir, la resistividad de un terreno imaginario, homogneo e istropo,
equivalente al terreno real heterogneo.

Las rocas gneas y metamrficas tpicamente tienen valores altos de resistividad.

La resistividad de estas rocas depende mucho del grado de fracturacin que
poseen, y el porcentaje de agua que rellena las fracturas del terreno. Las rocas
sedimentarias comnmente son ms porosas y tienen un alto contenido de agua,
lo que normalmente hace disminuir los valores de resistividad.

Los suelos mojados y el agua fresca del terreno tienen an valores ms bajos de
resistividad. Los suelos arcillosos normalmente tienen valores de resistividad ms
bajos que el suelo arenoso. Sin embargo, se destaca que existe una superposicin
en los valores de resistividad de las diferentes clases de rocas y suelos. Esto es
debido a que la resistividad de un muestreo particular del suelo o roca depende de
un nmero de factores tales como la porosidad, el grado de saturacin de agua y
la concentracin de sales disueltas.

Los metales, tal como hierro, tienen valores de resistividad sumamente bajos.
Los qumicos, tal como sodio y cloruro de potasio, pueden reducir mucho el valor
de resistividad del agua del terreno a menos de 1 ohm-m con moderadas
concentraciones. El efecto de electrolitos, tal como el cido actico, es

5
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

comparativamente menor. Los compuestos de hidrocarburos, tpicamente tienen

valores muy altos de resistividad.

1.3. PROPIEDADES ELECTRICAS DE LOS MATERIALES TERRESTRES

Los flujos de corriente en materiales terrestres a profundidades someras se

producen por medio de dos efectos principales. La conduccin electrnica y la
conduccin electroltica. En la conduccin electrnica, el flujo de corriente es por
medio de electrones libres, tal como los metales. En la conduccin electroltica, el
flujo de corriente es a travs del movimiento de iones en el agua subterrnea. En
los estudios de ingeniera y ambientales, la conduccin electroltica es
probablemente el mecanismo ms comn. La conduccin electrnica es
importante cuando existen minerales conductivos, tales como los sulfuros
metlicos y el mineral de grafito.

Metales (Au, Ag, Cu, Sn, grafito): enlace metlico con gran movilidad de los
electrones de valencia, en consecuencia los electrones tienen gran capacidad de
desplazamiento y la resistividad es muy baja (~ 10-8-10-7 ohm-m). Sin embargo la
conduccin metlica no es importante a una escala macroscpica dado que la
proporcin de metales nativos es muy baja, prevaleciendo la resistividad de las
impurezas que en general es mucho ms alta.

Semiconductores (Bornita, magnetita, pirrotita, etc.): No existe la misma

facilidad para el movimiento de los electrones, siendo necesaria la presencia de
campos elctricos fuertes para producir la conduccin. La conductividad de los
semi-conductores aumenta con la temperatura (intrnseca) y con las impurezas.

Dielctricos (Micas, feldespatos, cuarzo, calcita, silicatos en general): En

estos compuestos los electrones estn fuertemente ligados a los tomos mediante
enlaces covalentes o enlaces inicos (electrolitos slidos). La mayor parte de los
minerales constituyentes de la corteza terrestre son dielctricos con resistividades
del orden de 107 ohm-m.

Electrolitos lquidos (soluciones salinas): En las soluciones acuosas los

portadores de la corriente elctrica son los cationes y aniones. Al aplicar un campo
elctrico los iones presentes en las soluciones acuosas tienden a desplazarse en
la direccin de las cargas opuestas neutralizando la polarizacin del medio. Se
llama movilidad de un ion a la velocidad (constante) que adquiere este en un
campo elctrico aplicado. La resistividad de un electrolito es la resistencia elctrica
presentada por un volumen de electrolitos de forma cbica y un metro de lado al
establecer una diferencia de potencial de 1 volt entre dos caras opuestas.

En resumen se aprecia que las propiedades elctricas de los constituyentes de las

rocas, parametrizadas en la resistividad elctrica, presentan una variabilidad de
ms de 10 rdenes de magnitud (~ 10-7 107 ohm-m). Si bien la mayor parte de
los constituyentes de las rocas son dielctricos, lo cual estara implicando que la

6
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

conduccin elctrica debera estar ausente en las estructuras corticales (medio

aislante). Afortunadamente todas las rocas presentan cierto grado de porosidad,
volumen que generalmente est ocupado total o parcialmente por electrolitos que
condicionan una conduccin inica de resistividad muy variable, que dependen
fundamentalmente de las caractersticas de las soluciones acuosas y el tipo de
porosidad y grado de interconexin entre poros.

La resistividad de minerales, rocas comunes, qumicos, aguas, materiales de suelo

y contaminantes (Keller y Frischknecht, 1966; Daniels y Alberty, 1966; Telford et
al. 1990) se muestran en la Tabla 2.

La resistividad de los minerales es muy alta en comparacin con la resistividad de

las rocas que los contienen, esto es debido a la influencia de la porosidad,
humedad y la resistividad de los fluidos que rellenan los poros, en la resistividad
de la roca. En la Tabla 2 se muestra este ejemplo, en donde los rangos de valores
comparativos de la resistividad de los minerales de cuarzo y calcita, y de rocas
que los contienen: arenas y calizas respectivamente. El cuarzo tiene un rango de
resistividad entre 1012 y 1014 ohm-m, las arenas estn compuestas de cuarzo y
su rango de resistividad puede varia de 5 a 103 ohm-m, dependiendo de la
resistividad del agua que la satura. Algo parecido sucede con la calcita, la cual
tiene un valor de resistividad del orden de magnitud entre 6 y 9 veces mayor que
la roca caliza.

Las rocas gneas y metamrficas tpicamente tienen valores grandes de

resistividad. La resistividad de estas rocas es altamente dependiente del grado y el
porcentaje de fracturas rellenas con agua subterrnea. De esta manera una roca
de este tipo puede tener un rango de resistividad, de aproximadamente 1000
hasta 10 millones de ohm-m, dependiendo s esta hmeda o seca.

Esta caracterstica es conveniente en la deteccin de zonas de fractura y otras

caractersticas climticas, tales como estudios de ingeniera y agua subterrnea.
Las rocas sedimentarias, las cuales son ms porosas y tienen alto contenido de
agua, normalmente tienen valores bajos de resistividad comparadas a las rocas
gneas y metamrficas. Los rangos de valores de resistividad de las rocas
sedimentarias es de aproximadamente 10 a 10000 ohm-m, pero la mayora de
ellas se encuentra por debajo de 1000 ohm-m. Los valores de resistividad son
dependientes en gran medida de la porosidad y el contenido de sales en el agua.
Los sedimentos no consolidados tienen valores de resistividad bajos, muy
parecidos a las rocas sedimentarias, con un rango de valores de
aproximadamente 10 a un mximo de 1000 ohm-m. Los valores de resistividad
van a depender de la porosidad (considerando que todos los poros estn
saturados) tal es el caso del contenido de arcillas. El suelo arcilloso tiene un valor
de resistividad ms bajo que un suelo arenoso. Sin embargo, existen valores de
resistividad de las diferentes clases de rocas y suelos se que son iguales o casi
iguales. Esto se debe a que la resistividad de una muestra de roca o de suelo en

7
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

particular depende de factores tales como la porosidad, el grado de saturacin de

agua y la concentracin de sales disueltas.

La resistividad del agua subterrnea vara de 1 hasta 100 ohm-m o ms en medios

terrestres, dependiendo de la concentracin de sales disueltas. La baja
resistividad del agua de mar (aproximadamente 0.2 ohm-m) es debido al
relativamente alto contenido de sal. Esto hace al mtodo de resistividad una
tcnica ideal para el mapeo de la interfase del agua de sal y del agua fresca en las
zonas costeras.

Una ecuacin muy simple de la relacin que existe entre la resistividad de una
roca porosa y el factor de saturacin del fluido es la Ley de Archie. Esta ley es
aplicable para cierto tipo de rocas y sedimentos, particularmente aquellos que
tienen un bajo contenido de arcilla. La conduccin elctrica se asume que es a
travs de los fluidos que rellenan los poros de la roca. La ley de Archie est dada
por
= (1.15)

Donde es la resistividad de la roca, w es la resistividad del fluido, es la

fraccin de la roca rellena con el fluido, mientras a y m son dos parmetros
empricos (Keller y Frischknecht, 1966). Para la mayora de las rocas, a es
aproximadamente 1 mientras m es aproximadamente 2. Para los sedimentos con
un contenido de arcilla significativo, se han propuesto otras ecuaciones ms
complejas (Olivar et al., 1990).

Tambin se muestran las resistividades de diversos tipos de minerales. Sulfuros

metlicos (tales como la pirrotita, galena y pirita) tienen tpicamente valores bajos
de resistividad de menos de 1 ohm-m. El valor de resistividad de un cuerpo
mineral en particular puede diferir grandemente de la resistividad de los cristales
individuales. Otros factores, tales como un cuerpo mineral encontrado en forma
natural (masivo o diseminado) tienen un efecto significativo. Por ejemplo, la pizarra
tiene un valor de resistividad bajo, similar a los sulfuros metlicos, los cuales dan
lugar a problemas en los estudios mineros. La mayora de los xidos, como la
hematita, no tienen un valor de resistividad significativamente bajo. Una de las
excepciones es la magnetita.

Los valores de resistividad de varios contaminantes industriales tambin se

pueden ver en la Tabla 2. Los metales, como el hierro, tienen valores de
resistividad extremadamente bajos. Los productos qumicos que son electrolitos
fuertes, tales como el cloruro de potasio y el cloruro de sodio, reducen en gran
medida la resistividad del agua subterrnea a menos de 1 ohm-m incluso en
concentraciones relativamente bajas. El efecto de los electrolitos dbiles, como el
cido actico, es comparativamente ms pequeo. Los hidrocarbonos, como el
xyleno (6.998x1016 ohm-m), tpicamente tienen valores de resistividad muy altos.
Sin embargo en la prctica el porcentaje de hidrocarbonos en la roca o suelo suele
ser bastante pequeo, y podra no tener un efecto significativo en la resistividad.

8
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Tabla 2. Valores de resistividad de rocas, minerales, aguas, sedimentos qumicos,

minerales y materiales (Tabla tomada de Urbieta, Tesis Maestra, 2009)

Rocas gneas y metamrficas Resistividad (m) Rocas sedimentarias Resistividad (m)

Granito 5x103-106 Arenisca 8-4x103

Granito alterado 1500-10000 Lutita 20-2x103
Granito muy alterado 100-600 Caliza 100-5x103
Basalto 10-1.3x107 Caliza compacta 1000-5000
Ande sita 4.5x104-1.7x102 Caliza blanda 100-300
Diorita 105 Caliza fracturada 500-1000
Gabro 103-106 Limonita 50-107
Sienita 102-106 Dolomita 3.5x102-5x103
Peridotita 3x103-6.5x103 Margas secas 102
Lavas 102-5x104 Margas hmedas 10
Pizarra 6x102-4x107
Mrmol 102-2.5x108
Cuarcita 102-2x108
Aguas, sedimentos y suelos Resistividad (m) Qumicos, minerales y Resistividad ( m)
materiales

Arenas 30-500 Hierro 9.74x10-8

Arena fina 103 0.01M cloruro de potasio 0.708
Arena seca 103-104 0.01M cloruro de sodio 0.843
Arena silcea 200-3000 0.01M cido actico 6.13
Arenas con agua dulce 50-500 Xileno 6.998x1016
Arenas con agua salada 0.5-5 Aluminio 2.65x10-8
Arena Petrolfera 4-2.2x102 Mercurio 95.8x10-8
Arcilla 1-100 Madera 108-1011
Arcilla seca 30 Vidrio 1010-1014
Arcilla hmeda 10 mbar 5x1014
Limos 10-200 Plata 1.47x10-8
Aluvin 10-800 Cobre 1.70x10-8
Gravas 50-103 Oro 2.22x10-8
Gravas hmedas 5x102 Acero 20.0x10-8
Gravas saturadas 102 Platino 11.0x10-8
Conglomerado 2x103-104 Tungsteno 5.25x10-8
Roca salada 3x10-1013 Nquel 6.40x10-8
Aguas dulces superficiales 10-103 Estao 1.5x10-8
Aguas salobres superficiales 2-10 Plomo 22.0x10-8
Agua subterrnea 10-100 Cromo 13.2x10-8
(no contaminada) Calcopirita 1.2x10-5-0.3
Agua marina 0.2 Pirita 2.9x10-5-1.5
Agua de lagos salados 0.1-1 Pirrotita 6.5x10-6-5x10-2
Cinabrio 2x107
Galena 3x10-5-3x102
Cobalto 3.5x10-4-10-1
Bauxita 2x102-6x103
Cuprita 10-3-300
Cromita 1-106
Hematita 3.5x10-3-107
Magnetita 5x10-3-5.7x103
Cuarzo 7.5x1017
Azufre 1015
Carbn 3.5x10-5
Anhidrita 109
Calcita 2x1012
Fluorita 8x1013
Yeso 1000
Grafito 60.00x10-8
Silicio 6.40x102

9
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.4. CLASIFICACIN DE LAS TCNICAS DE PROSPECCIN ELECTRICA

DE CORRIENTE CONTINUA

Existen infinidad de configuraciones de campo para los mtodos elctricos. La

diferencia entre las distintas configuraciones reside en el tipo de transmisor y
receptor de corriente utilizados en el campo, y en cmo se siten en el campo
estos transmisores o receptores. Los mtodos elctricos tienen la ventaja de ser
mucho ms sencillos en el procesado y la interpretacin pero su mayor desventaja
es que cualquier material aislante o muy conductivo puede enmascarar los
materiales que lo rodean.

Los valores de resistividad tienen un gran rango comparado a otras cantidades

fsicas obtenidas por otros mtodos Geofsicos. La resistividad de las rocas y el
suelo en un rea de estudio pueden variar por varios rdenes de magnitud. En
comparacin, la densidad usada en investigaciones gravimtricas usualmente
cambia por debajo de un factor de 2. Esto hace que los mtodos geoelctricos
sean una tcnica geofsica muy verstil.

Las tcnicas de prospeccin elctrica de corriente continua generalmente se

dividen en:

1. Calicata Elctrica
2. Sondeo Elctrico Vertical (SEV).
a. Schlumberger
b. Wenner
3. Tomografa Elctrica
a. Wenner
b. Wenner-Schlumberger
c. Dipolo-Dipolo
d. Polo-Polo
e. Polo-Dipolo

1.4.1. Calicata Elctrica

Son medidas de potencial mediante las cuales se determinan resistividades del

terreno hasta una profundidad fija. Es un mtodo de investigacin horizontal a
profundidad constante para detectar estructuras que se presentan como
heterogeneidades laterales de resistividad.

Las distancias entre los electrodos de corriente y de potencial son constantes, de

manera que el arreglo de electrodos se mueve sin cambiar su forma en la
superficie del terreno. La interpretacin de este mtodo es cualitativa.
Se aplica en el estudio de las variaciones laterales de resistividad en el subsuelo,
pero no puede detectar cambios verticales de resistividad.

10
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.4.2. Sondeo Elctrico Vertical (SEV)

La finalidad del sondeo elctrico vertical (SEV) es averiguar la distribucin vertical

en profundidad de las resistividades aparentes bajo el punto sondeado a partir de
medidas de la diferencia de potencial en la superficie. Se utiliza sobre todo para
detectar y establecer los lmites de capas horizontales de suelo estratificado.

La profundidad de penetracin de la corriente elctrica depende de la separacin

de los electrodos de corriente AB. Si la distancia entre los electrodos AB aumenta,
la corriente circula a mayor profundidad pero u densidad disminuye. Para un
medio istropo y homogneo, el 50% de la corriente circula por encima de la
profundidad AB/2 y el 70.6% por encima de una profundidad AB (Orellana, 1982).
La mayor eficacia del mtodo corresponde al caso en los que los sondeos se
efectan sobre un terreno compuesto por capas lateralmente homogneas en lo
que respecta a la resistividad y limitada por planos paralelos a la superficie del
terreno, medio estratificado. Es tolerablemente vlida para estratos inclinados
hasta unos 30.

Durante mucho tiempo, en prospeccin elctrica, la profundidad de investigacin

ha sido considerada sinnimo de la profundidad de penetracin de la corriente. Sin
embargo, el efecto de una capa en los potenciales o campos observados en
superficie no depende nicamente de la densidad de corriente que la atraviesa.

La profundidad de investigacin caracterstica es la profundidad a la que una capa

delgada de terreno (paralela a la superficie) contribuye con participacin mxima a
la seal total medida en la superficie del terreno. La profundidad de investigacin
viene determinada por la posicin de los electrodos de corriente y de potencial, y
no slo por la penetracin o distribucin de la corriente. Es decir, si se
intercambian entre s las posiciones de los electrodos de potencial con los de
corriente, la distribucin de las lneas de corriente cambia.

Sin embargo, la resistividad aparente y por tanto la profundidad de investigacin

no cambia. Definiendo L como la distancia entre los dos electrodos extremos (sin
considerar los situados en el infinito), los mismos autores determinan la
profundidad de investigacin de diversos dispositivos elctricos en un suelo
homogneo, siendo para el dispositivo polo-polo de 0,35L, para Schlumberger de
0,125L y para Wenner de 0,11L. Barker (1989) la define la profundidad de
investigacin efectiva como la mitad de la seal medida en la superficie debida a
la porcin de suelo superior a esa profundidad y la otra mitad de la seal debida a
la porcin de suelo inferior y muestra la mayor utilidad de sta sobre la utilizada
por Roy y Apparao (1971). Las profundidades de investigacin efectiva para los
dispositivos Wenner, Schlumberger y doble dipolo son respectivamente de 0,17L,
0,19L y 0,25L (para este ltimo la profundidad de investigacin caracterstica es de
0,195L), es decir ligeramente mayores que utilizando la definicin de profundidad
de investigacin caracterstica (Roy y Apparao, 1971).

11
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

[Link]. Efectos laterales en el SEV y ambigedades en su interpretacin

Si el dispositivo elctrico est prximo a un contacto vertical, las lneas de

corriente sern distorsionadas por lo que VM N se ver afectado por el otro
medio, tanto ms cuanto mayor sea la separacin de los electrodos AB. Por lo
tanto, la medida de la resistividad aparente en un SEV est influida por la
distribucin de resistividades en un cierto volumen de terreno. Esto implica que
para distancias AB grandes no se sabr si la resistividad aparente es debida a
cambios de estructuras en la profundidad o a las heterogeneidades laterales por
contraste de resistividades (Orellana, 1982).
Puede ocurrir que las curvas de resistividad aparente para dos casos diferentes de
SEV sean idnticas si la relacin entre profundidad a la que se encuentra n
estrato y su resistividad permanece constante, lo que provoca una ambigedad en
la deduccin del grosor de la capa y su resistividad.

[Link]. Aplicaciones

El SEV es aplicable cuando el objetivo tiene una posicin horizontal y una

extensin mayor que su profundidad. Tal es el caso del estudio de capas
tectnicas, bsqueda de agua subterrnea, infiltracin de agua marina en zonas
costeras, estratificacin y todo tipo de situaciones geolgicas donde el modelo 1D
se aproxime al modelo verdadero. El SEV no es adecuado para contactos
verticales tales como, fallas, diques, cavernas y discordancias.

1.4.3 Descripcin de la tcnica de Tomografa de Resistividad Elctrica 2D

La Tomografa de Resistividad Elctrica 2D es un mtodo de resistividad multi-

electrodo (Figura 1.2), basado en el modelado bidimensional de valores de
resistividad del subsuelo mediante el empleo de tcnicas numricas (elementos
finitos o diferencias finitas).

12
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Figura 1.2. Esquema de la Tomografa de Resistividad Elctrica 2D

La aparicin de la Tomografa de Resistividad Elctrica ha revolucionado de

manera contundente con respecto a los mtodos de resistividad convencionales,
estas tcnicas aunque se han estado utilizando durante varias dcadas en
estudios de filtracin de agua, su limitada resolucin 2D les confera en general un
papel secundario frente a otras tcnicas (i.e. Potencial espontneo).

Este avance se debe fundamentalmente a tres razones:

En primer lugar porque estos mtodos convencionales solo utilizan 4

electrodos, siendo necesario para cada medida variar manualmente sus
posiciones en el terreno (proceso relativamente lento y pesado). En cambio
el mtodo de Tomografa de Resistividad Elctrica es una tcnica multi
electrodos, en el cual todo el proceso de adquisicin de datos est
totalmente automatizado. Esto permite realizar un gran nmero de
mediciones de resistividad, tanto en profundidad como lateralmente, en un
espacio de tiempo relativamente rpido (del orden de 300 medidas en una
hora y media), obteniendo con esto modelos 2D de gran resolucin.

En segundo lugar porque al emplear diferentes tcnicas numricas

permitir procesar eficazmente todo este gran volumen de informacin.

En tercer lugar despus de procesar este gran volumen de informacin se

obtiene una seccin o imagen elctrica 2D en falso color para identificar con
mayor facilidad las diferentes anomalas o zonas de inters. Estas
imgenes elctricas son las que se utilizan para la interpretacin final.

13
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.5. DISPOSITIVOS ELECTRODICOS UTILIZADOS EN LA TCNICA DE

TOMOGRAFA DE RESISTIVIDAD ELCTRICA 2D.

Un conjunto de cuatro electrodos recibe el nombre de dispositivo electrdico. La

medicin del potencial elctrico en un semiespacio homogneo se hace mediante
un dispositivo electrdico o un arreglo geomtrico de electrodos. Consta de cuatro
electrodos, A y B los electrodos de corriente, positivo y negativo respectivamente,
M y N los electrodos de potencial, entre los cuales se mide la diferencia de
potencial del campo elctrico creado por A y B.

Existen diversos tipos de dispositivos, lineales, simtricos, asimtricos,

ecuatoriales. En principio, los cuatro electrodos pueden adoptar cualquier
disposicin geomtrica sobre el plano que representa la superficie del terreno.

La mejor eleccin de configuracin para realizar un estudio de campo depende del

tipo de estructura a ser mapeada (forma, dimensiones), profundidad, sensibilidad
del aparato y nivel de ruido de fondo de la zona de estudio.

A continuacin se describen las caractersticas de los cinco tipos de dispositivos

empleados para la adquisicin de datos con la tcnica de Tomografa de
Resistividad Elctrica

1.5.1. Dispositivo Wenner-Schlumberger

Se trata de una composicin simtrica de los electrodos AMNB dispuestos en

lnea, donde la distancia de los electrodos detectores MN es mucho menor que la
de los inyectores AB.

Figura 1.3. Esquema de configuracin electrdica del dispositivo Schlumberger.

= ( + 1) (1.16)

14
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.5.2. Dispositivo Wenner

En este dispositivo los electrodos AMNB son equidistantes entre si, a una
distancia a. En la figura se observa esquemticamente la disposicin de los
electrodos.

Figura 1.4. Esquema de la configuracin electrdica del dispositivo Wenner.

= 2 (1.17)

1.5.3. Dispositivo Polo-Dipolo

En este dispositivo un electrodo de corriente se lleva a una gran distancia

(tericamente en el infinito) de los otros tres (AMN).

Figura 1.5. Esquema de la configuracin electrdica del dispositivo Polo- Dipolo.

= 2 + 1 (1.18)

15
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.5.4. Dispositivo Dipolo-Dipolo

En este dispositivo los electrodos se disponen sobre una lnea en el orden ABMN
formando as un doble dipolo.

Figura 1.6. Esquema de la configuracin electrdica del dispositivo Dipolo- Dipolo.

= + 1 + 2 (1.19)

1.5.5. Dispositivo Polo-Polo

En este dispositivo los electrodos B y N se llevan a una gran distancia

(tericamente en el infinito) de los otros dos A y M.

Figura 1.7. Esquema de la configuracin electrdica del dispositivo Polo- Polo.

= 2 (1.20)

16
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.6 EQUIPO REQUERIDO EN LA TECNICA DE TOMOGRAFIA DE

RESISTIVIDAD ELECTRICA 2D.

El equipo (Figura 1.8) que se utiliza para la ejecucin del mtodo consta de:

Electrodos: barras metlicas a modo de electrodos cuyo nmero variar en

funcin del problema.
Cable: de gran longitud, con conexiones cada cierto intervalo a fin de poder
conectar los electrodos.
Conectores: pequeos cables de cobre con pinzas en los extremos. Son
los que se encargan de conectar los electrodos al cable.
Computadora porttil: Con el programa de adquisicin instalado se procede
a determinar y exportar las variables de trabajo las cuales son: tipo de
dispositivo electrdico, nmero de electrodos, espaciado entre ellos,
nmero de niveles de profundidad, nmero de medidas a realizar, etctera,
y al terminar la adquisicin de datos, se importan, procesan, visualizan e
interpretan los resultados con el programa de inversin 2D.
Resistivimetro: Se encarga de ejecutar de forma automtica toda la
secuencia de medidas predeterminada, verifica el buen estado de las
conexiones, almacena digitalmente todos los valores de resistividad.
Una batera: es la fuente de alimentacin de todo el sistema.

Figura 1.8. Equipo de campo requerido para adquisicin de datos de Resistividad.

17
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.7 PLANIFICACION DE UN ESTUDIO DE CAMPO.

Los aspectos ms importantes que se definen en la fase de planificacin de la

campaa de campo son:

Delimitar el o los objetivos (resolucin de mapeo requerida).

El nmero de perfiles a realizar y su ubicacin.
Profundidad mxima de investigacin.
El dispositivo electrdico, longitud total de los perfiles, nmero de
electrodos, separacin entre electrodos y niveles de profundidad.

La Tomografa de Resistividad Elctrica se basa en la obtencin de perfiles 2D de

resistividad del subsuelo, por lo en el caso de analizar posibles filtraciones en
alguna estructura lineal (i.e. presa), bastar simplemente con disponer del perfil (o
perfiles) lo ms paralelo posible a la traza de la estructura.

Para investigar la variacin de resistividad en dos dimensiones se puede utilizar

cualquier configuracin de cuatro electrodos y las mediciones pueden ser hechas
en cualquier secuencia.

En cuanto al nmero total de electrodos a emplear ste es muy variable,

regularmente un nmero cercano a 50 electrodos es bastante frecuente. Estos se
disponen de forma equidistante oscilando de 1 a 10 m para problemas ingenieriles
(Griffiths et. al., 1990), si bien para problemas de filtraciones de agua un
espaciado entre 1 y 2 metros es lo ms aconsejable.

La eleccin del mejor dispositivo para realizar un trabajo de campo depende de

diversos factores, tales como: el tipo de estructura geolgica a estudiar, la
sensitividad del dispositivo, la profundidad de investigacin, la intensidad de la
seal y el nivel de ruido ambiental.

Dos aspectos muy importantes a tener en cuenta, son la profundidad mxima de

investigacin y el rea de definicin del mtodo, aspectos ambos estrechamente
ligados.

Al utilizar la tcnica de Tomografa de Resistividad Elctrica el campo de

interaccin disminuye con la profundidad. En consecuencia, a la hora de planificar
el trabajo de campo se debe tener presente esta prdida de la intensidad del
campo elctrico en los laterales del perfil.

En cuanto a la profundidad mxima de investigacin, sta depende esencialmente

de la longitud del perfil (a mayor longitud mayor profundidad), del dispositivo
electrdico utilizado el cual tambin influye en cierta medida y del nmero de
niveles.

18
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

1.8 CLCULO DE LA PROFUNDIDAD DE INVESTIGACIN.

La profundidad a la que se sita espacialmente cada medida de resistividad (punto

de atribucin) es un tema de cierta ambigedad, dado que sta no solo depende
de la disposicin de los electrodos (ms espaciados mayor profundidad), sino
tambin de la distribucin de las resistividades del medio.

La Tabla 3 muestra los factores para calcular la profundidad media de

investigacin (ze) de las diferentes modalidades. L es la longitud total del arreglo.
Note que existen valores idnticos para las modalidades Wenner-Schlumberger y
polo-dipolo. El factor geomtrico esta dado para un espaciamiento entre
electrodos a de un metro. Las figuras de los arreglos se pueden observar en el
punto 1.5 de este captulo.

La profundidad media de investigacin nos da una idea de la profundidad a la cual

se puede penetrar en el subsuelo con una modalidad en particular. Los valores de
la profundidad media son determinados integrando la funcin de sensitividad con
la profundidad (Edwards, 1977).

Esta profundidad media no depende de la resistividad aparente medida o de la

resistividad del modelo de Tierra homogneo. Se debe notar que las
profundidades son estrictamente validas para un modelo de Tierra homogneo,
pero son probablemente lo suficientemente buenas para planear estudios de
campo. Por ejemplo, si existe un contraste de resistividad grande cerca de la
superficie, la profundidad de investigacin podra ser un poco diferente.

Para determinar la profundidad mxima mapeada por un estudio en particular, se

multiplica el mximo espaciamiento entre electrodos a, o la mxima longitud del
arreglo L, por el factor apropiado de nivel de profundidad dado en la Tabla 3
(Loke, 2004). Por ejemplo, si el mximo espaciamiento entre electrodos a usado
por el arreglo Wenner es de 100 metros (o la L mxima es de 300 metros),
entonces la profundidad mxima mapeada es alrededor de 51 metros. Para las
modalidades dipolo-dipolo, polo-dipolo y Wenner-Schlumberger, otro factor que
adems debe tomarse en consideracin es el factor n. Para las modalidades con
cuatro electrodos finitos (tales como dipolo-dipolo, Wenner, y Wenner-
Schlumberger), es probablemente ms sencillo utilizar la longitud L. Esta L no
es la longitud total del tendido es la L mxima situada para un nivel de
profundidad en particular, es decir, si se requiere calcular la profundidad mxima
para el arreglo dipolo-dipolo para una separacin entre electrodos a = 10 metros y
una n =6, la L mxima ser L = a(BA)+na(AM)+a(MN) = 10+6(10)+10 = 80
metros, esto multiplicado por el factor = 0.216 para n = 6 (Tabla 3) resulta una
profundidad mxima de penetracin aproximada de 17 metros.

19
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Tabla 3. Profundidad media de investigacin (ze/a y ze/L) para las diferentes modalidades.
El factor geomtrico est calculado para el valor de a igual a 1 metro (Loke, 2004).

Tipo de Modalidad z e /a z e/L Factor Geomtrico

Wenner 0.519 0.173 6.2832
Dipolo-dipolo n=1 0.416 0.139 18.850
n=2 0.697 0.174 75.398
n=3 0.962 0.192 188.50
n=4 1.220 0.203 376.99
n=5 1.476 0.211 659.73
n=6 1.730 0.216 1055.6
n=7 1.983 0.220 1583.4
n=8 2.236 0.224 2261.9
Wenner-Schlumberger n=1 0.519 0.173 6.2832
n=2 0.925 0.186 18.850
n=3 1.318 0.189 37.699
n=4 1.706 0.190 62.832
n=5 2.093 0.190 94.248
n=6 2.478 0.191 131.95
n=7 2.863 0.191 175.93
n=8 3.247 0.191 226.19
n=9 3.632 0.191 282.74
n=10 4.015 0.191 345.58
Polo-dipolo n=1 0.519 12.566
n=2 0.925 37.699
n=3 1.318 75.398
n=4 1.706 125.66
n=5 2.093 188.50
n=6 2.478 263.89
n=7 2.863 351.86
n=8 3.247 452.39
Polo-polo 0.867 6.28319

1.9 ADQUISICIN DE DATOS DE RESISTIVIDAD APARENTE.

En primer lugar se procede a colocar los electrodos en la superficie de forma

equiespaciada, intentando siempre que formen una lnea recta. Se debe definir
cul es el origen y final del perfil a fin de no tener problemas en la fase de
interpretacin. Se intentar adems que la topografa del terreno sea lo ms plana
posible para evitar tener fenmenos anmalos en las medidas de resistividad.

Los electrodos se clavarn en el terreno lo suficiente como para garantizar un

buen contacto electrodo-suelo (resistencia de contacto baja). Una vez colocados,
se conectan al cable por medio de los conectores caimn-caiman.

20
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Ahora, a travs de la computadora porttil o en el propio resistivimetro se

configuran todas las variables para construir el perfil 2D: tipo de dispositivo,
nmero de electrodos, abertura entre electrodos, nmero de niveles.

Concluida la configuracin, si se realiz a travs de la computadora, se transfieren

estas especificaciones al resistivimetro para que ste realice automticamente
toda la secuencia de mediciones prefijada, para obtener una seccin 2D de
resistividades aparentes.

Sin embargo, antes de iniciar las mediciones, mediante el resistivimetro, se

procede a la verificacin de que todas las conexiones (electrodos-cable) funcionan
correctamente, y que la resistencia de contacto entre electrodo-suelo es
suficientemente baja menores a 100 ohm-m para garantizar buenos resultados. En
el caso de tener valores superiores generalmente del rango entre 500 y 1000 ohm-
m basta con clavar un poco ms los electrodos, o humedecer el terreno para
mejorar la conductividad.

Una vez que est todo dispuesto, se le indica al resistivimetro que inicie la
secuencia de medidas, el proceso tardar dependiendo del nmero de medidas a
realizar. En general obtener el orden de 500 medidas se necesita
aproximadamente una hora y media.

Por ltimo, terminado todo el proceso de captura de datos, la informacin es

almacenada digitalmente dentro del resistivimetro y posteriormente se exporta a la
computadora porttil con el fin de proceder a su procesado e interpretacin con el
programa de inversin 2D correspondiente.

1.10 PROCEDIMIENTO PARA LA INVERSION DE DATOS.

El objetivo de llevar a cabo la inversin, es obtener un modelo 2D de resistividades

reales a partir de la seccin de resistividades aparentes 2D obtenida. Un modelo
es una representacin matemtica idealizada de una seccin de la tierra, los
parmetros del modelo son cantidades fsicas que se obtienen de los datos
observados.

La construccin de la seccin 2D de resistividades aparentes, involucra el total de

datos de resistividad transversales como una seccin a profundidad donde cada
valor de resistividad aparente es graficado como si fuera la resistividad verdadera
de un punto en una profundidad proporcional al espaciamiento de los electrodos.
La configuracin preliminar de los datos proporciona una imagen aproximada de la
distribucin de la resistividad.

El conjunto de datos de resistividad aparente en forma de una seccin se emplean

para generar imgenes elctricas. Las mediciones en la superficie de las
21
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

variaciones de la resistividad aparente generan un modelo de la distribucin

espacial de la resistividad aparente del subsuelo. Todo mtodo de inversin
procura determinar un modelo idealizado del subsuelo que reproduzca los
parmetros medidos de manera adecuada.

Estos algoritmos de inversin se basan en realizar todo un proceso iterativo con el

fin de obtener un modelo de resistividades reales, tal que el modelo de
resistividades aparentes generado, sea lo ms parecida a la seccin de
resistividades aparentes 2D medida en el campo.

1.11. VENTAJAS Y LIMITACIONES DE LA TECNICA.

Esta tcnica presenta las siguientes ventajas:

Debido a que el proceso de captura de los datos est totalmente

automatizado, se pueden realizar un gran nmero de medidas, de manera
que se obtienen secciones elctricas 2D de excepcional resolucin tanto en
profundidad como lateralmente. Esta visualizacin bidimensional de alta
resolucin, con el mtodo de Tomografa de Resistividad Elctrica permite
lo siguiente:

Detectar zonas susceptibles de presentar focos de filtracin

(contaminantes, fugas en presas, infiltracin agua marina, etc.), as como
caracterizar geomtricamente el problema, deteccin de fracturas,
cavidades, tneles, galeras y estudios arqueolgicos.

Abordar de forma eficaz problemas con elevada complejidad estructural e

incorporar datos topogrficos dentro del programa de inversin para realizar
una correccin por topografa.

1. Ubicar con gran eficiencia otro tipo de procedimientos (i.e. sondeos o

piezmetros), dado que marca con mayor exactitud el lugar y profundidad a
la que se encuentra la anomala.

El procesamiento de los datos se puede realizar en el propio campo,

mediante la ayuda de una computadora porttil en un tiempo muy reducido.
Este procesamiento, casi en tiempo real, le confiere al mtodo una gran
versatilidad a la hora de resolver posibles problemas surgidos durante toda
la logstica y adquisicin de datos en campo, o decidir la ejecucin de
nuevos perfiles que inicialmente no estaban proyectados.

Se puede controlar la profundidad mxima de estudio.

Gran variedad de sistemas de resistividad multi-electrodos y programas de

inversin 2D.

22
 Fundamentos Tericos de los Mtodos de Corriente Continua Captulo 1

Las principales limitaciones de este mtodo son:

A diferencia de otros mtodos, la Tomografa de Resistividad Elctrica

requiere de un equipo y programas de inversin 2D relativamente caros.

Como pasa en todos los mtodos de resistividad elctrica, la presencia de

materiales muy resistivos en superficie y a profundidad imposibilita un buen
funcionamiento del mtodo, siendo necesario la utilizacin de otras
tcnicas, en este caso las electromagnticas, en los que no se precisa un
contacto fsico con el suelo.

Se limita a estudios muy someros entre los 50 y los 300 metros.

23
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

CAPITULO 2. El PROBLEMA DE LA INVERSION DE DATOS DE

RESISTIVIDAD APARENTE

2.1. INTRODUCCION

El anlisis cuantitativo de los datos obtenidos en el campo se realiza resolviendo

el problema directo y/o el problema inverso. El problema directo predice resultados
cuando la geologa y la fsica del medio son conocidas, a partir de la distribucin
de las resistividades del subsuelo se calcula la resistividad aparente en superficie
para una configuracin en particular. De esta forma es posible comparar los datos
de campo con modelos conocidos y as verificar cuantitativamente si estos datos
provienen de una distribucin de resistividades similar a la del modelo propuesto.
Por su parte, el problema inverso construye un modelo razonable del subsuelo a
partir de los datos obtenidos en el campo. Es decir, a partir de los valores de
resistividad aparente medidos en superficie se obtiene la distribucin de
resistividades del subsuelo.

2.2. EL PROBLEMA DIRECTO.

Antes de encarar el problema de la inversin de los datos de una prospeccin

geofsica es necesario comprender cmo se resuelve el problema directo. Este
consiste en hallar una expresin para la diferencia de potencial entre cualquier par
de puntos en el espacio como una funcin de la distribucin de resistividades del
medio, de la configuracin utilizada y de la corriente de entrada. A partir de los
valores de potencial calculados se obtienen los de resistividad aparente.

La expresin para los potenciales se obtiene a partir de la aplicacin conjunta de

la ley de Ohm, del principio de conservacin de la carga y de las ecuaciones de
Maxwell para campos estacionarios. En general no es posible despejar la solucin
exacta del problema directo, por lo tanto se recurre a herramientas numricas para
resolverlo, obtenindose una solucin aproximada. La deduccin de las
ecuaciones que definen el problema directo y su resolucin son complejas. Slo
se mencionarn algunos aspectos y consideraciones generales que deben tenerse
en cuenta a la hora de implementarlo.

Una vez encontradas las ecuaciones que definen el problema, antes de resolverlo
se debe asumir a priori algn tipo de distribucin de resistividades del subsuelo
particular. El subsuelo puede modelarse en forma unidimensional (1D),
bidimensional (2D) o tridimensional (3D). Cada dimensin constituye una direccin
a lo largo de la cual la resistividad del suelo puede variar.

2.3. DESCRIPCION DE UN MODELO 2D.

En el modelo bidimensional se consideran variaciones en los valores de

resistividad del subsuelo tanto en profundidad como en una de las direcciones
laterales. Para poder calcular el potencial elctrico en la superficie, dada una
24
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

distribucin de resistividades como la mencionada, el problema se aborda

bsicamente aplicando dos mtodos: el de diferencias finitas (Loke, 1995) y el de
elementos finitos (Coggon, 1971). Cada uno de estos mtodos encara la
resolucin del problema directo de modos distintos y utilizando principios fsicos y
herramientas de clculo numrico diferentes, que no se explicarn dada su
complejidad. Adems requieren modelos de distribuciones de resistividad
bidimensionales diferentes, los que mostraremos a continuacin.

El modelado 2D por el mtodo de diferencias finitas (Mufti, 1976; Loke, 1995)

requiere la divisin o discretizacin del subsuelo en bloques rectangulares,
utilizando una rejilla rectangular (Figura 2.1). En sta, cada bloque posee una
resistividad constante.

El mtodo de diferencias finitas determina el potencial en los nodos de la rejilla

rectangular y por lo tanto se pueden obtener los valores de resistividad aparente a
lo largo del perfil determinado por la rejilla. Destacando que no se pueden obtener
valores de resistividad aparente a lo largo de una direccin diferente a la definida
por la rejilla. Si la rejilla es lo suficientemente fina la distancia entre i y j es
pequea comparada con la longitud del perfil, se pueden modelar estructuras
geolgicas complejas. Para la tcnica de elementos finitos tambin se discretiza el
subsuelo en elementos de resistividad constante pero stos no tienen por qu ser
rectngulos. Se obtienen los potenciales en los nodos de la rejilla y a partir de
stos los valores de resistividad aparente a lo largo del perfil.

Figura 2.1 Subespacio dividido en celdas rectangulares de resistividad constante.

25
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

2.4. TEORA DE SENSITIVIDAD

Un campo elctrico primario se define como aquel a ser medido en la presencia de

un medio conductor uniforme o. Los campos secundarios resultan de
inhomogeneidades empotrados en este medio conductivo. Si se perturba la
conductividad de un pequeo volumen de material, se puede calcular la
sensitividad de la medida del cambio del potencial debida a esa perturbacin.
Matemticamente, una pequea perturbacin de la conductividad por unidad de
perturbacin es conocida como la derivada de Frechet. Por lo tanto, el cambio en
la respuesta como resultado de la presencia de un pequeo cuerpo de resistividad
particular es, en efecto, la sensitividad de una seccin del subsuelo medida debida
a una perturbacin en la conductividad para una localizacin fija de una fuente
puntual en la superficie o en pozo. Las grficas de sensitividad son derivadas a
partir de los clculos de las derivadas de Frechet.

Para ilustrar la sensitividad consideremos la ecuacin de potencial para un

semiespacio inhomogneo
2 2

= 4 2 3 , , (2.1)
2 2 2 2 2 3
+ + [()2 +()2 + ]2

Definiendo la resistividad como:

= (2.2)

Donde:

1 < < +1 , < < +1

= (2.3)
0

El problema de sensitividad, se obtiene de la derivada de Frechet de () con

respecto a como:

( ) 1 2 2
= 4 2 =1 =1 2 2 2 3/2 (2.4)
[ + + ]3/2 [ + + ]

El valor de sensitividad es una medida de la cantidad de informacin acerca de la

resistividad de un modelo de bloques contenido en el conjunto de datos medidos.
El valor ms alto de sensitividad, es el valor ms confiable del modelo de
resistividad. En general, los bloques cerca de la superficie a menudo tienen
valores altos de resistividad debido a que la funcin de sensitividad tiene valores
muy grandes cerca de los electrodos. Los bloques en los lados y en el fondo
tambin tienen valores de sensitividad altos debido a que muchos de estos

26
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

bloques son de gran tamao y adems se extienden a los extremos de la rejilla de

diferencia-finita o elemento-finito.

2.5. EL PROBLEMA INVERSO

Opuestamente a lo que sucede en el problema directo, en el problema inverso los

valores de resistividad aparente que se obtuvieron en las distintas mediciones de
campo o sintticas constituyen datos y lo que se desea obtener es la distribucin
de resistividades del subsuelo que se generaron con dichos valores. En otras
palabras, la aplicacin del problema inverso permite obtener un modelo de la
distribucin de resistividades del subsuelo a partir de los valores de resistividad
aparente medidos en el campo. De esta manera, las ecuaciones que definen al
problema inverso son las mismas que las del problema directo. Este es un
problema auto consistente que se debe resolver en forma iterativa. La resolucin
del problema inverso requiere la eleccin a priori del tipo de modelo del subsuelo
con el que se va a trabajar, 1D, 2D o 3D segn resulte conveniente. Debido a que
los valores de resistividad aparente tienen imprecisiones y/o ambigedades
inherentes al proceso de medicin, el modelo obtenido no debera reproducir en
forma exacta los valores de resistividades aparentes medidas sino dentro de un
determinado rango de error. La diferencia entre los valores de resistividad
aparente con el modelo calculado y los obtenidos en el campo se denomina
desajuste o discrepancia. Para que el modelo obtenido sea realista, el desajuste
debe ser del orden de los errores de medicin. Adoptar en la inversin un
desajuste muy pequeo, puede resultar en estructuras sin sentido fsico que son
una manifestacin del ruido de la medicin. Por otro lado, si el desajuste es muy
grande se puede perder informacin de la estructura interna del subsuelo.

El proceso de inversin de los datos a grandes rasgos es el siguiente:

1) Se parte de un modelo de resistividades inicial 1D, 2D o 3D.

2) Se calculan los valores de resistividad aparente predichos por el
mismo.
3) Se comparan estos valores con los medidos en el campo, es decir se
calcula la discrepancia
4) Si la discrepancia es mayor que el esperado, se modifica el modelo
inicial y se repite el proceso.
5) Si la discrepancia es el esperado, quiere decir que se resolvi el
problema.

El algoritmo de resolucin del problema inverso, por lo expuesto hasta aqu, debe
incluir un algoritmo para la resolucin del problema directo, para llevar a cabo el
punto 2).

Los problemas inversos no tiene solucin nica sino que hay infinitos modelos que
pueden reproducir los datos en forma adecuada. Debido a esto se le indica al
algoritmo de inversin que seleccione un modelo particular de todos los posibles,

27
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

es decir, que el modelo posea caractersticas especficas. De esta forma se logra

una gran flexibilidad para generar diversos modelos y es posible incorporar
informacin adicional acerca de la distribucin de resistividades del medio y
generar un modelo acorde con esa informacin.

2.6. TEORA BASICA DE INVERSION

La teora de inversin en general, est organizada en un conjunto de tcnicas

matemticas y estadsticas (clculo, algebra matricial, estimacin e inferencia
estadstica, etc.), que tienen la finalidad de reproducir informacin de las
propiedades fsicas de un sistema (e.g. la Tierra) partiendo de los datos
observados en este sistema.

En la inversin geofsica se trata de encontrar un modelo con una respuesta que

sea similar a los valores reales medidos. El modelo es una representacin
matemtica idealizada de una seccin de la Tierra. El modelo tiene un conjunto de
parmetros que son las cantidades fsicas que se pretende estimar a partir de los
datos observados. La respuesta del modelo son los datos sintticos que se
calculan a partir de las relaciones matemticas las cuales definen el modelo de un
conjunto de parmetros dados. Todos los mtodos de inversin esencialmente
tratan de determinar un modelo para el subsuelo cuya respuesta se asemeje con
los datos medidos sujetos a ciertas restricciones.

2.7. PLANTEAMIENTO DEL PROBLEMA NO LINEAL.

El conjunto de datos observados se puede escribir como un vector columna y

dado por:

= 1 , 2 . (2.5)

donde m es el nmero de medidas, yj puede ser la resistividad aparente medida o

su logaritmo en base e; el uso del logaritmo permite reducir el rango de valores
con los que se va a trabajar, lo cual es deseable debido a que los valores tpicos
de resistividad oscilan de menos 1 ohm-m a ms de 1000 ohm-m.

Los parmetros del modelo se representan por el siguiente vector columna q

= (1 , 2 , . . , ) (2.6)

donde qj representa las resistividades reales de los n bloques en los que se divide
la malla, es decir, el nmero de los parmetros del modelo.

La respuesta del modelo f se puede escribir en una forma similar

= (1 , 2 , . . . ) (2.7)

28
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

La respuesta del modelo contiene los valores de resistividad aparente calculados

(o sus logaritmos) a partir de los parmetros iniciales q.

Para problemas de resistividad, es comn utilizar los logaritmos de los valores de

resistividad aparente para los datos observados y la respuesta del modelo, y los
logaritmos de los valores del modelo para los parmetros del modelo.

A partir de aqu se plantea que la solucin al problema directo (calculo de

resistividades aparentes a partir de los valores de resistividad real del modelo) se
expresa como yi=Fi[q] donde Fi es una funcin directa no lineal que acta sobre
los parmetros q del modelo discretizado para producir una respuesta del modelo.

La diferencia entre los datos observados (resistividades aparentes de las

mediciones) y la respuesta del modelo (resistividades calculadas a partir de los
parmetros estimados del modelo), esta dado por el vector de discrepancia g que
est definido por

= (2.8)

Adems es posible considerar que cada dato yi tiene una varianza j.

Para un conjunto de variables aleatorias independientes X1,,XN, la probabilidad

P(X1 = x1,,XN = xN) es igual al producto =1 = . Entonces si se supone
que los errores de las mediciones yi son independientes y estn distribuidos de
manera gaussiana, la probabilidad conjunta de los mismos est dada por:

1 2
=1 2 (2.9)

Los valores de parmetros que maximizan esta probabilidad son los estimadores
de mxima certidumbre del modelo. Se puede probar que maximizar esta
expresin es equivalente a minimizar la cantidad denominada chi cuadrado:

2
2 = =1 (2.10)

Por lo que entonces el ajuste por mnimos cuadrados resulta ser un estimador de
mxima certidumbre (Press et al. 1997).

Es por esto que la cantidad definida en la ecuacin puede utilizarse como una
manera de estimar la bondad de un ajuste. En forma matricial, esta ecuacin se
expresa como:
2
2 = [] (2.11)

donde W es una matriz diagonal de peso de M x M que se define como:

29
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

= diag {1/1,1/2, ,1/} (2.12)

Para un modelo inicial m1, los primeros dos trminos del desarrollo de la serie de
Taylor dan la siguiente aproximacin:

+ = + + (2.13)

donde es el resto, y Ji es la matriz jacobiana de M x N de derivadas parciales.

Los elementos de J estn dados por:

= (2.14)

que representa el cambio en la i-sima respuesta del modelo F[qi] debido al

cambio en el j-simo parmetro del modelo. La diferencia

= (2.15)

es una pequea perturbacin alrededor del modelo inicial. F se puede aproximar

quitando el trmino .

Para reducir la discrepancia entre resistividad aparente calculada y medida, se

utiliza la ecuacin de Gauss-Newton para determinar el cambio en los parmetros
del modelo de tal manera que debe reducir la suma cuadrtica de los errores
(Lines y Treitel, 1984)

= (2.16)

Despus de calcular el vector de cambio del parmetro, se obtiene un nuevo

modelo mediante:

= + (2.17)

En la prctica, la ecuacin (2.16) de mnimos cuadrados rara vez se utiliza en la

inversin geofsica. En algunas situaciones el producto de la matriz JT J puede ser
singular, y por lo tanto, la ecuacin de mnimos cuadrados no tiene solucin para
q. Otro problema comn es que el producto de la matriz JT J es casi singular.
Esto puede ocurrir si se introduce un modelo inicial muy diferente al modelo que
se busca ajustar. Podra ocurrir entonces que el residuo de un ajuste por mnimos
cuadrados creciera en lugar de decrecer (Pelton et al. 1978). Para evitar esas
divergencias, resulta necesario imponer restricciones en q para asegurar la
estabilidad del proceso iterativo.

El cambio del vector de parmetros calculado (Ec.2.16) puede tener componentes

demasiado grandes, de tal modo que el nuevo modelo q2 (Ec.2.17) podra tener
valores que no son reales. Un mtodo comn para evitar este problema se logra
30
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

con la modificacin de Marquardt-Levenberg (Lines y Treitel, 1984) en la ecuacin

de Gauss que est dada por

( + ) = (2.18)

donde I es la matriz identidad. El factor es conocido como el Marquardt o factor

de amortiguamiento, y a este mtodo tambin se le conoce como el mtodo de
regresin (Inman, 1975). El factor de amortiguamiento efectivamente condiciona el
rango de valores, cuando se cambian componentes que el vector de parmetros
q puede tomar. Si bien el mtodo de Gauss Newton en la ecuacin (2.16) intenta
minimizar la discrepancia en la suma cuadrtica, la modificacin del mtodo de
Marquardt-Levenberg tambin minimiza la magnitud del vector discrepancia y el
cambio del vector de parmetros. Este mtodo ha sido ampliamente utilizado en la
inversin de datos de sondeo de resistividad, en donde el modelo consiste de un
pequeo nmero de capas. Sin embargo, cuando se tiene un nmero de
parmetros del modelo muy grande, tales como en modelos de inversin 2D y 3D
los cuales consisten de un gran nmero de celdas pequeas, el modelo que
produce este mtodo puede tener una distribucin de resistividad errtica en
zonas de resistividad con altos y bajos espurios (Constable et al. 1987).

Otra opcin es incorporar un suavizado en la ecuacin de mnimos cuadrados, con

lo cual se minimiza la rugosidad del modelo. Para una estructura de dos
dimensiones donde y es la direccin del perfil, una forma de medir la rugosidad, la
inversa de la suavidad, puede darse con la siguiente expresin

= + (2.19)

donde y es una matriz de rugosidad que deriva lateralmente los parmetros del
modelo de los adyacentes. z es una matriz de rugosidad que deriva en la
direccin vertical los parmetros del modelo de los rectngulos contiguos. Ambas
son aproximaciones al operador gradiente por el mtodo de diferencias finitas (de
Groot-Hedlin, 2000). Supongamos que se tiene una malla que consiste de un total
de m elementos, con p elementos en la direccin horizontal, cada uno de un ancho
h, y l elementos en la direccin vertical, con anchos i, i = 1l (se esta suponiendo
una malla en donde la escala horizontal es constante, pero la escala vertical vara
para compensar la prdida del poder de resolucin). Si se enumeran los
elementos de izquierda a derecha comenzando en la parte superior izquierda, la
matriz de rugosidad z de M M es:

1 0 0 1 0
0 1 0 0 1
= (2.20)
1
0

31
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

donde 0 es una matriz de ceros de p M. Hay p 1 ceros entre las entradas en

las columnas de z. Por lo tanto z acta de manera tal de derivar los parmetros
del modelo entre los bloques adyacentes verticales. La matriz de rugosidad
horizontal y de M M est dada por:

1 0
2
= (2.21)

0

donde es la matriz de rugosidad para la capa i , cuya expresin es:

/ /
/ / 0
= .. (2.22)
/ /
0 0

Por tanto, yi deriva los parmetros del modelo entre bloques adyacentes en la
capa i, asignando un peso a las diferencias segn la relacin entre el ancho y el
alto de los bloques de la grilla de discretizacin.

En este caso el problema se formula como la minimizacin de la siguiente funcin:

2 2 2
= + + (2.23)

donde es un multiplicador de Lagrange. La funcin U es minimizada en puntos

donde el gradiente con respecto al modelo es cero. Las funciones de los datos son
no lineales, por lo que la funcin U debe ser linealizada y resuelta iterativamente.

El resultado de minimizar la ecuacin (2.23) produce un sistema de ecuaciones

lineales:

+ = (2.24)

donde es un factor de amortiguamiento (multiplicador de Lagrange), y C una

matriz que es una combinacin lineal (peso) de productos de matrices de
suavizado con sus traspuestas:

= + (2.25)

Este mtodo alternativo minimiza variaciones en los parmetros del modelo

(Modelo de mnimos cuadrados con suavizado), utiliza la norma l2. Por lo que
tiende a producir un modelo con variaciones suaves de la resistividad, este modelo

32
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

solo es aceptable si realmente la resistividad del subsuelo vara de manera

gradual.

La ecuacin de Gauss Newton de mnimos cuadrados se modifica de tal manera

que se reducen al mnimo las variaciones espaciales de los parmetros del
modelo (igualmente, el modelo de valores de resistividad cambia de una manera
suavizada o gradual). La formulacin usual de mnimos cuadrados aplica un
suavizado slo sobre el vector q y no directamente sobre los valores de
resistividad. Se puede formular un mtodo (Modelo de mnimos cuadrados con
suavizado de valores de resistividad) donde tambin los valores de resistividad
cambien suavemente. Este mtodo de mnimos cuadrados con condicionamiento
suavizado (Ellis y Oldenburg, 1994a) tiene la siguiente forma:

+ = (2.26)

El mtodo convencional de mnimos cuadrados que minimiza el cuadrado de la

diferencia entre los datos calculados y las mediciones en general da mejores
resultados si las mediciones tienen fluctuaciones al azar o gaussianas. En el caso
que el proceso sea gaussiano el valor esperado de 2 es N, la cantidad total de
datos. Esto significa que el valor esperado para la cantidad definida como:

1
= 2 (2.27)

y denominada error cuadrtico medio, es 1.

Sin embargo, si el ruido no tiene estas caractersticas aleatorias, proviene de

problemas en el equipo de medicin o de equivocaciones, el estimador de mxima
certidumbre se obtiene minimizando la desviacin absoluta media entre valores
medidos y calculados (Press et al. 1997), y suele dar resultados ms
satisfactorios. Al sustituir el cuadrado por el valor absoluto (es decir, este modelo
alternativo utiliza la norma l1), la probabilidad asociada con este estimador decae
ms lentamente al aumentar las discrepancias. Este modelo es entonces menos
sensible a pequeas separaciones respecto de los supuestos del ajuste, y se lo
denomina robusto; tiende a producir regiones homogneas y con lmites abruptos,
por lo que conviene aplicarlo cuando previamente se esperan estas caractersticas
en el subsuelo.

Como se dijo anteriormente la ecuacin (2.26) procura reducir a un mnimo el

cuadrado de los cambios espaciales, o rugosidad, del modelo de valores de
resistividad. Es de hecho un mtodo de optimizacin con condicionamiento
suavizado con una norma l2. Esto tiende a producir un modelo con una suave
variacin de los valores de resistividad. Esta aproximacin es aceptable si la
resistividad del subsuelo vara de forma paulatina y gradual. En algunos casos, la
geologa del subsuelo se compone de una serie de regiones que son casi
homogneas internamente pero con fuertes lmites entre las diferentes regiones.

33
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

Para tales casos, la formula de inversin en la ecuacin (2.26) se puede modificar

de tal manera que minimiza los cambios en los valores de resistividad del modelo
(Claerbout y Muir, 1973). Esto algunas veces da resultados significativamente
mejores. Tcnicamente esto se conoce como un mtodo de optimizacin con
condicionamiento suavizado con norma l1, o ms comnmente conocida como
mtodo de inversin de bloques (Modelo robusto). Se pueden utilizar un gran
nmero de tcnicas para tal modificacin. Un mtodo simple para aplicar una
norma l1 basada en la optimizacin del mtodo de mnimos cuadrados es la
formulacin iterativa del mtodo de mnimos cuadrados repesados (Wolke y
Schwetlick, 1988). La optimizacin en la ecuacin (2.26) se ha modificado para
obtener un modelo robusto

+ = (2.28)

Donde:

= + (2.29)

donde Rd y Rm son las matrices de peso que introducen diferentes elementos en la

discrepancia y en los vectores del modelo rugoso dndoles el mismo peso en el
proceso de inversin.

La ecuacin (2.28) proporciona un mtodo general que puede ser objeto de

modificaciones si se incluye informacin conocida sobre la geologa del subsuelo.
Por ejemplo, si se sabe que las variaciones de la resistividad del subsuelo son
comnmente confinadas a zonas limitadas, los valores del factor de
amortiguamiento pueden modificarse (Ellis y Oldenburg, 1994a) lo cual permite
cambios ms grandes en estas zonas.

2.8. METODOS DE DISCRETIZACION PARA UN MODELO 2D

En la seccin anterior, se vio que el mtodo de mnimos cuadrados se utiliza para

calcular ciertas caractersticas fsicas del subsuelo, los "parmetros del modelo",
partiendo de las mediciones de resistividad aparente. Los "parmetros del modelo"
se establecen dividiendo en secciones o cubos las diferentes regiones del
subsuelo. La Figura 2.2 muestra las diferentes posibilidades que se pueden
emplear.

El mtodo ms comnmente utilizado en modelos 2D y 3D es una interpretacin

basada nicamente en generar un modelo que divide la seccin de inters en
celdas rectangulares y asigna a cada una de estas regiones valores de
resistividad. Las posiciones de las celdas son fijas y slo se permiten variar a las
que contienen valores de resistividad durante el proceso de inversin. Los
parmetros del modelo son las resistividades de cada una de las celdas. Para el
ejemplo mostrado en la Figura 2.2, los parmetros del modelo son setenta y dos
celdas de valores de resistividad 1 a 72.

34
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

Una aproximacin radicalmente diferente es una frontera basada en el mtodo de

inversin. Este mtodo subdivide el subsuelo en diferentes regiones. La
resistividad se supone que es homognea dentro de cada regin. A la resistividad
se le permite cambiar de manera arbitraria a travs de las fronteras, lo cual es til
en las zonas con una repentina transicin en la geologa. La resistividad de cada
regin y las profundidades a las fronteras son modificadas por el mtodo de
optimizacin de mnimos cuadrados con el fin de que los valores de resistividad
aparente calculados coincidan con los valores observados. Los "parmetros del
modelo" para el ejemplo que se muestra en la Figura 2.2, son dos valores de
resistividad (1 a 2) y de las profundidades en cinco puntos (z1 a z5) resultando
a lo largo de la frontera un total de siete parmetros. Si bien este mtodo funciona
bien para datos sintticos, para conjuntos de datos de campo grandes, este
esquema numrico puede dar lugar a resultados inestables, con altas oscilaciones
en las fronteras (Olayinka y Yaramanci, 2000). Su mayor limitacin es,
probablemente, la hiptesis de una resistividad constante dentro de cada regin.
En particular, los cambios laterales de resistividad cerca de la superficie tienen un
efecto muy grande sobre los valores de resistividad aparente medidos. Este
modelo no toma en cuenta esos cambios laterales, son a menudo errneamente
modelados como cambios en las profundidades de las fronteras.

Esfuerzos recientes han estado combinando los modelos basados en celdas y de

frontera en los mtodos de inversin (Smith et al. 1999). Uno de esos mtodos es
el de inversin lateralmente condicionado (Auken y Christiansen, 2004). En este
mtodo, los cambios laterales (pero no verticales) son permitidos en cada regin
(Figura 2.2c), as como las transiciones abruptas a travs de las fronteras. Los
"parmetros del modelo" para el ejemplo de la Figura 2.2 son: veinticuatro valores
de resistividad (1 a 24) y trece puntos en profundidad (z1 a z13) a lo largo de la
frontera, con lo cual se obtiene un total de treinta y siete parmetros. La
informacin de otras fuentes, tales como pozo o datos ssmicos, se pueden utilizar
para proporcionar una estimacin inicial de la profundidad a la frontera. Una
situacin comn es cuando la informacin detallada est disponible en un solo
pozo. En este caso, la frontera inicial se fija a una profundidad constante. El
mtodo de inversin ajusta las profundidades a travs de una serie de puntos
situados a lo largo de la frontera durante el proceso de inversin. Posteriormente
se aplica un condicionamiento suavizado para reducir al mnimo los cambios en
las profundidades entre puntos adyacentes en la misma frontera (Smith et al.
1999). Este mtodo funciona bien, particularmente cuando el subsuelo se
compone de varias zonas sedimentarias.

Otra generalizacin de este concepto es que permite cambios tanto verticales

como laterales dentro de cada regin (como en un modelo basado nicamente en
celdas) y tambin cambios en las formas a travs de las fronteras (Figura 2.2d). El
modelo que se muestra en la Figura 2.2d tiene setenta y dos valores de
resistividad y cinco valores de profundidad, lo que da un total de setenta y siete
parmetros del modelo. Este tipo de discretizacin es particularmente til cuando

35
 El problema de la Inversin de Datos de Resistividad Aparente Captulo 2

cerca de la superficie existen inhomogeneidades que se producen dentro de la

capa superior las cuales a su vez producen un gran efecto en la medicin de
valores de resistividad aparente. El modelo de celdas utilizadas hasta el momento
tienen formas rectangulares (Figuras 2.2a, 2.2c y 2.2d). Esto en parte debido al
uso del mtodo de diferencias finitas con el cual se calculan los valores de
resistividad aparente del modelo. Una ligera desventaja es que la frontera es
aproximada por una serie de pasos rectangulares. La Figura 2.2e muestra una
posible variacin cuando se utiliza el mtodo de elemento finito con celdas
trapezoidales en donde los bordes de las celdas adyacentes a la frontera son
ajustadas con la finalidad de obtener la forma verdadera de la frontera.

Figura 2.2. Los diferentes modelos para el subsuelo utilizados en la interpretacin de datos de
estudios de imgenes elctricas 2D. (a) Modelo compuesto nicamente por celdas. (b) Modelo
compuesto nicamente por fronteras. (c) Modelo condicionado lateralmente. (d) Modelo combinado
por celdas y fronteras compuesto por celdas rectangulares y (e) Modelo con fronteras compuesto
por celdas trapezoidales.

La finalidad de la inversin de los datos provistos por una TRE es hallar un

modelo, una representacin matemtica idealizada, de la distribucin de
resistividades de una seccin del subsuelo.

36
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

CAPTULO 3. CDIGO DE INVERSIN DE DATOS DE TOMOGRAFA DE

RESISTIVIDAD ELCTRICA ZONDRES2D

3.1. INTRODUCCIN

Despus de un estudio sobre el terreno, las mediciones de resistencia

normalmente se reducen a valores de resistividad aparente. Prcticamente todos
los sistemas comerciales de electrodos mltiples vienen con el programa de
computadora para realizar esta conversin.

En este captulo se describen las caractersticas de tres programas que son utiles
para el procesamiento e interpretacin de datos de Tomografa de resistividad
elctrica 2D. El programa ProfileR (Andrew Binley, Universidad de Lancaster)
realiza tanto modelado directo, como la inversin de datos de resistividad
bidimensionales, sin embargo no tiene forma de visualizarlos grficamente, por lo
que requiere de una interfaz grafica como puede ser el programa Surfer, o, como
se emplea en este trabajo, el programa ZONDRES2D. ProfileR requiere un
formato especial de archivos ([Link]) para poder realizar el modelado, por ello
se implemento el programa ProfilerV2 hecho en Matlab para crear este tipo de
archivos de una manera un poco ms sencilla y as poder manipularlos mediante
ProfileR. Por otra parte tambin se describen las caractersticas del programa
ZONDRES2D (ZOND GEOPHYSICAL SOFTWARE) el cual aparte de servir como
interfaz grafica para los archivos generados por ProfileR, est diseado para el
modelado directo 2D y la inversin de perfiles datos de resistividad.

3.2. PROGRAMA DE INVERSIN ZONDRES2D.

Este algoritmo genera un modelo de las variaciones laterales y verticales de la

distribucin de la resistividad del subsuelo por medio de perfiles. El algoritmo
utiliza cuatro mtodos de inversin por mnimos cuadrados: Smoothness
constrained, Occam, Marquardt y una combinacin de los ltimos dos.

El mtodo de inversin Smoothness constrained es en realidad el mtodo de

inversin ptimo para el anlisis de problemas en donde tengamos variaciones
graduales de resistividad en el terreno, como por ejemplo en filtraciones de agua o
vertidos de contaminantes. Adems este mtodo, basado en el mtodo de
mnimos cuadrados de Gauss-Newton, presenta la gran ventaja de que no
requiere de un modelo inicial, a partir del cual empezar el proceso iterativo.

El algoritmo Occam o de modelo suave no requiere un modelo inicial de capas

sino simplemente la resistividad de un semiespacio homogneo. Esta tcnica
constrie el proceso de inversin minimizando la variacin vertical de la
resistividad en un gran nmero de capas delgadas a la vez que simultneamente
minimiza el error de ajuste entre las respuestas observadas y calculadas. Aunque
este mtodo es prcticamente inmune al sesgo del intrprete, hay un nmero de

37
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

parmetros que tienen influencia en el modelo final, tal como el nmero de capas,
las profundidades a las interfases mnima y mxima, y la suavidad del modelo
final.

El mtodo de inversin Marquardt interpola entre los Algoritmos de Gauss-Newton

y el mtodo de descenso de gradiente. Este algoritmo es ms robusto que el
Algoritmo de Gauss-Newton, lo que significa que en muchos casos se encuentra
una solucin, incluso si se inicia muy lejos del mnimo final. Por otra parte, por
buen comportamiento funciones y parmetros razonables de partida, el algoritmo
de Marquardt tiende a ser un poco ms lento que el Algoritmo de Gauss-Newton.

El proceso de inversin minimiza la diferencia entre los valores de resistividad

aparente calculados y medidos. El suavizado reduce la cantidad de ruido en el
modelo de resistividad. El algoritmo determina la resistividad de los bloques
rectangulares que son los parmetros del modelo y producen una seccin de
resistividades aparentes que corresponde a los datos obtenidos. El nmero de
bloques rectangulares coincide con el mismo nmero de datos obtenidos con la
seccin de resistividades aparentes y adems son arreglados de manera similar.
El mtodo Newton con regularizacin se aplica en la inversin lo cual mejora la
estabilidad de la solucin y permite una distribucin ms uniforme de la
resistividad en el modelo inverso. La diferencia entre los valores medidos y
modelados est dada por el error cuadrtico medio (RMS). Sin embargo, el
modelo con el error RMS ms bajo posible, no siempre es el ms apropiado,
algunas veces, puede mostrar variaciones no realistas en los valores de
resistividad del modelo y adems, no siempre ser el mejor modelo desde el
punto de vista geolgico. En general la aproximacin ms prudente es elegir el
modelo para una iteracin, despus del cual el RMS no cambia significativamente.
Esto usualmente ocurre entre la tercera y quinta iteracin. El mtodo convencional
de mnimos cuadrados minimiza el cuadrado entre los datos calculados y las
mediciones, y tiende a darle gran importancia a datos con una discrepancia muy
grande. Esto lo hace sensitivo a datos errneos, por problemas en las mediciones.

El programa ZONRES2D tiene un nmero de parmetros que pueden ser

modificados por el usuario para la obtencin de resultados que estn ms cerca
de la geologa conocida. Tiene variedad de visualizacin de datos de perfil, editor
de datos y un sistema para analizar sensitividad y la resolucin de la matriz de
medicin tambin estn disponibles. Datos de campo o modelados, as como la
diferencia entre ellos y pesos de las lecturas se puede mostrar en forma de
grficos o secciones de resistividades aparentes. En el editor de electrodos, el
usuario puede leer la informacin sobre los parmetros de la matriz, editar lecturas
o establecer su peso. El programa permite establecer el peso de acuerdo a la
seal de valor o sensibilidad de las lecturas de los parmetros de la seccin. Los
elementos del modelo tambin pueden ser fijos si no afectan a los resultados de
los modelos a seguir. La funcin de anlisis de sensibilidad permite al usuario
estimar el efecto de un elemento determinado modelo en adelante los resultados
de simulaciones.

38
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

El programa reconoce una variedad de formatos de archivo para la inversin de

datos, que permite al usuario procesar los datos incluidos en la muestra con
diferentes tipos de matrices de forma simultnea. Una caracterstica importante del
programa ZONDRES2D es que reconoce un formato particular de datos de tipo
ProfileR. Las caractersticas de este tipo de archivos se describen en el siguiente
punto de este captulo. ZONDRES2D permite salvar los datos invertidos en
formato ascii para poder ser ledos por el programa de interpolacin SURFER con
la finalidad de obtener las imgenes elctricas finales con una mejor interpolacin
que los dos programas anteriores. Esto es til para estudiar la resolucin del
modelo que se puede obtener a travs de distintas estructuras mediante el uso de
diferentes arreglos y configuraciones de electrodos.

3.3. CARACTERSTICAS DE PROFILER

ProfileR est diseado para obtener imgenes elctricas 2D por medio de la

construccin de un perfil de resistividades con un arreglo de electrodos definido en
superficie. El cdigo es una solucin inversa de la distribucin de resistividad 2D
que se basa en una funcin objetivo, el cual se regulariza combinndose con los
mnimos cuadrados ponderados (una solucin de tipo Occam).

El perfil de entrada requiere muy pocos datos para la generacin de la malla de

elementos finitos. Los requisitos mnimos que el usuario debe especificar son el
nmero de electrodos y las coordenadas de los mismos. ProfileR puede
considerar la variacin en la topografa si el usuario especifica la posicin de los
electrodos intermedios y la elevacin. La separacin de los electrodos no tiene
que ser necesariamente uniforme.

La malla generada se compone de una regin frontal y una regin de fondo. La

regin frontal es el rea investigada en el estudio y la variacin de resistividades
en esta rea ser la salida del programa. La regin de fondo se emplea para tener
en cuenta las condiciones de los bordes.

La malla se genera con dos elementos finitos entre los electrodos en la horizontal
de la regin frontal. En la regin de fondo se extiende de izquierda a derecha de la
matriz de electrodos utilizando un aumento exponencial elementos. En los
elementos verticales aumenta el tamao con la profundidad.

La regin est parametrizada en trminos de bloques de resistividad mediante la

agrupacin de parches de elementos. En la regin frontal un bloque de parmetros
se define como un bloque de 2 por 2 de los elementos.

3.3.1. Especificaciones del archivo.

ProfileR solo requiere de un archivo de entrada: profiler. In (Ver caractersticas en

la seccin 3.2.2.)

39
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

ProfileR genera una serie de archivos de salida:

[Link] que contiene el registro principal de la ejecucin.

[Link] que contiene las coordenadas vertical y horizontal de la regin
frontal. La primera lnea muestra el nmero de coordenadas en el archivo.
Le siguen dos columnas con la posicin horizontal y vertical,
respectivamente. Este archivo puede ser usado para superponer la frontera
en Surfer, por ejemplo.
[Link] que contiene las coordenadas de los electrodos. Este puede
ser til para superponer en una imagen la resistencia final.

Si selecciona una solucin inversa entonces obtendr dos archivos adicionales:

[Link] que contendr el resultado de resistividad a la solucin inversa.

El formato de salida tiene cuatro columnas: la columna 1 es la posicin
horizontal, la columna 2 es la elevacin, la columna 3 es la resistividad (en
m), y la columna 4 es la resistividad en log10.
[Link] el contiene cinco columnas. En la primera columna estn los
datos normalizados. La segunda columna contiene los datos observados
como resistividad aparente. La tercera columna contiene las resistividades
aparentes equivalentes para el modelo calculado. La cuarta columna
muestra el peso original de los datos. La quinta columna es el peso final de
los datos (un * puede indicar si algn peso fue modificado).

Si selecciona una solucin directa obtendr el siguiente archivo:

[Link] el cual contiene siete columnas. La primera columna contiene el

nmero de mediciones. Las columnas 2 a 5 contienen los electrodos
empleados para esta medicin. La columna 6 contiene la resistencia y la
columna 7 la resistividad aparente. NOTA: el clculo de la resistividad
aparente asume que el suelo en superficie es horizontal y que z=0
corresponde a el suelo. Si este no es el caso ignore los valores de
resistividad aparente. Tambin en la versin actual el modelo directo es
calculado para una resistividad uniforme de 100 m.

3.3.2. Detalles de [Link]

Lnea 1: Encabezado (No ms de 80 caracteres).

Lnea 2: numero de electrodos (numero entero), profundidad de estudio (nmero

real), resistividad inicial (nmero real), tipo de datos (nmero entero), numero de
polos (nmero real).

Donde el nmero de electrodos puede ser mximo de 128, la profundidad de

estudio es la mxima profundidad (en metros); la resistividad inicial es la
resistividad uniforme con la que comenzara para la solucin inversa en (m); el

40
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

tipo de datos es 0 si los datos utilizados en la inversin no estn transformados a

log, 1 si los datos estn transformados a log y 2 cuando es un modelado directo; el
nmero de polos es el numero de polos situados en el infinito. Para un arreglo
polo-dipolo el nmero de polos ser igual a 1, para un arreglo polo-polo el nmero
de polos ser 2, para cualquier otro tipo de arreglo el nmero de polos cera cero.

Lnea 3: nmero de posiciones de electrodos (nmero entero).

Donde el nmero de posiciones de electrodos es el nmero de posiciones en las

cuales se indica la posicin horizontal y la elevacin de cada electrodo. Cada
posicin al inicio y fin de electrodos debe ser declarada por lo que el nmero de
electrodos debe ser mayor o igual a 2.

Lnea 4: electrodo (nmero entero), electrodo_x (nmero real), electrodo_z

(nmero real).

Donde electrodo es el nmero de electrodo, electrodo_x es la posicin horizontal

(en metros) y electrodo_z es la elevacin (en metros), siendo positivo hacia arriba
del nivel del suelo. Esta lnea se repite para el nmero de posiciones de los
electrodos.

Lnea 5: Error de datos (nmero real).

Donde el error de datos es el tpico % de error. Este valor debe ser mayor que
cero y sus valores regularmente se encuentran entre 2.0 y 5.0. NOTA: Omitir la
lnea 5 si se selecciona un modelado directo.

Lnea 6: nmero de mediciones (nmero entero).

Este nmero es el nmero total de mediciones que se realizaron para el estudio

(mximo 1500).

Lnea 7: n (numero entero), elec1 (nmero entero), elec2 (nmero entero), elec3
(nmero entero), elec4 (nmero entero), R (nmero real).

Donde n es el nmero de medicin, elec1 es el numero de electrodo (no

coordenada) del electrodo de potencial positivo; elec2 es el numero de electrodo
de potencial negativo; elec3 es el numero de electrodo de corriente positivo y
elec4 es el numero de electrodo de corriente negativo, y R es la resistencia en ().

Si el nmero de polos es igual a 1 se omite la columna elec4. Si el nmero de

polos es 2 entonces las columnas de elec2 y elec4 deben ser omitidas.

Nota: R solo se requiere cuando se desea un modelo inverso. Si se desea un

modelo directo solo se requieren los nmeros de electrodos de corriente y
potencial.

41
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

La lnea 7 se repite para todo el nmero de mediciones.

3.4. Implementacin de ProfilerV2.

Como vimos el programa ProfileR nos permite realizar modelado directo e inverso
de datos de TRE, para ello requiere archivos con un formato especfico
([Link]), por ello se implemento un cdigo (ProfilerV2) realizado en Matlab
(versin 8.0) para el pre diseo este tipo de archivos, para poder as generar
archivos para poder modelarlos mediante ProfileR y debido a que dicho programa
no permite visualizar los resultados, estos puedan ser visualizados mediante el
programa de inversin ZONDRES2D o por cualquier otro programa que reconozca
este formato. La finalidad de Profiler V2 es facilitar al usuario la captura manual de
todos los datos requeridos por ProfileR, por lo que tambin ahorra tiempo durante
el modelado.

ProfilerV2 requiere como datos de entrada, algunas caractersticas de la

adquisicin como son nmero total de electrodos que se utilizaron, el
espaciamiento entre ellos, si el levantamiento fue con o sin topografa etc., como
se muestra en la Figura 3.1.

Figura 3.1. Cartula del programa ProfilerV2. Se puede apreciar todos los datos que dicho
programa necesita para trabajar.

El programa permite crear archivos tanto para un modelado directo o inverso,

tambin se ajusta si el levantamiento fue adquirido con o sin topografa.

42
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

En el caso de elegir la opcin de modelado directo los campos Tipo de datos para
la inversin y Error se desactivaran debido a que estos datos slo son
necesarios para crear un archivo de modelado inverso. De igual forma si elige la
opcin de un levantamiento sin topografa el campo Posiciones de electrodos se
desactivara.

3.4.1. Creacin de un archivo tipo [Link] con ayuda de ProfilerV2

1. En Matlab abrir el archivo ProfilerV2.m y aparecer la ventana como se

muestra en la Figura 3.1.
2. Seleccionar el tipo de modelado (directo o inverso), el tipo de levantamiento
(con o sin topografa) y el tipo de arreglo (Wenner, Wenner-Schlumberger,
Polo-Polo, Dipolo-Dipolo o Polo-Dipolo).
3. Llenar los campo de acuerdo con la siguiente informacin:

- # de electrodos: corresponde al nmero total de electrodos que se

emplearon para realizar el levantamiento.
- Posiciones de electrodos: Como se menciono antes, este campo
estar activo si se selecciona la opcin de tipo de levantamiento con
topografa. Como se puede ver en la figura 3.1 en este campo esta una
opcin de Cargar, al dar click aparecer una pantalla como se muestra
en la Figura 3.2., aqu podr seleccionar un archivo de formato txt o xls,
el cual debe contener 3 columnas, la primera indicar el numero de
electrodo, la segunda la posicin en x de ese electrodo y la tercera la
posicin en z. (Figura 3.3)

Figura 3.2. Despliegue de pantalla para poder buscar el archivo que contiene las posiciones de los
electrodos.

Cabe mencionar que no se deben especificar las posiciones de todos los

electrodos sino solo las del primer y ltimo electrodo y de aquellos donde su cota
difiera de cero, como se puede observar en la Figura 3.3.

43
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

Figura 3.3. Formato del archivo con las posiciones de los electrodos.

- Espaciamiento: Es el espaciamiento que existe entre electrodos.

- n: es el nmero de nivel al cual se desea profundizar con el estudio.
Este dato es til para calcular la profundidad media de investigacin de
acuerdo a lo mencionado en la seccin [Link].
- Resistividad inicial: esta es la resistividad que se propone para el
medio. De entrada se propone una de 100 ohm-metro pero el usuario
puede proponer la que mejor crea conveniente.
- # de mediciones: es el nmero total de lecturas que se tomaron
durante el estudio.
- Tipo de datos para la inversin: como mencione esta opcin solo
estar disponible si se desea un tipo de modelado inverso. Aqu el
usuario solo tendr que llenar con un 0 o un 1 segn sea el caso. 0 si
sus datos de resistividad aparente NO estn convertidos en base
logartmica, y 1 si dichos datos lo estn.
- Datos: Al igual que en el campo de Posiciones de electrodos en este
campo se debe cargar un archivo de tipo txt o xls, el cual contendr solo
una columna con los datos de resistividad que se obtuvieron en el
levantamiento (Figura 3.4). Hay que mencionar que el nmero de datos
debe ser igual que el nmero de mediciones reportado. Estos datos
como se menciona en el punto anterior pueden estar o no en base
logartmica.

44
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

Figura 3.4. Formato del archivo que contiene los datos de resistividad del levantamiento. La figura
muestra datos en escala logartmica.

- Error: El usuario deber llenar este campo con el error que considere
apropiado para valorar la inversin. Como ya se menciono este
generalmente tiene valores de entre 2 y 5.
- Archivo de Salida: esta opcin contempla dos campos: Ruta, que ser
la ubicacin en el procesador a donde desea que se guarde el archivo
que se generar; y un campo en blanco donde podr poner un ttulo
para el encabezado de dicho archivo, este ttulo ira aunado a la
modalidad empleada (Ver Figuras 3.5, 3.6 y 3.7)
- Por ltimo se debe dar click en el botn Crear y se generara el
archivo con llamado [Link]. Lo podr localizar en la ruta de acceso
que se sealo en el punto anterior y al abrirlo podr observar un archivo
como el que se muestra en los ejemplos de las figuras 3.5, 3.6 y 3.7.

45
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

Figura 3.5. Ejemplo de un archivo generado por el programa Profiler V2. Este ejemplo muestra el
formato de un modelado directo, hecho con un levantamiento sin topografa, y se utiliz un arreglo
Polo-dipolo para la adquisicin de los datos.

46
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

Figura 3.6. Ejemplo de un archivo generado por el programa Profiler V2. Este ejemplo muestra el
formato de un modelado inverso, hecho con un levantamiento sin topografa, y se utiliz un arreglo
Wenner para la adquisin de los datos.

47
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

Figura 3.7. Ejemplo de un archivo generado por el programa Profiler V2. Este ejemplo muestra el
formato de un modelado inverso, hecho con un levantamiento con topografa, y se utiliz un arreglo
Wenner para la adquisicin de los datos.

Una vez generado el archivo [Link] puede abrir el ejecutable de ProfileR,

para poder generar los modelos directos o invertir los datos, se generaran ciertos
archivos segn sea el caso y despus podr visualizar la inversin mediante el
programa ZONDRES2D (Ver anexo).

48
 Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

3.5. FUNCIONAMIENTO DE ZONDRES2D

Para trabajar con el programa ZONDRES2D se necesita estar en una carpeta

aislado de los archivos de datos y del ProfileR.

A continuacin se describe el procedimiento general para emplear el programa:

1. Abrir ZONDRES2D desde la ubicacin donde se encuentre. (Figura 3.8)

Figura 3.8. Ubicacin del ejecutable del programa ZONDRES2D.

2. Abrir archivo de datos. Como el programa reconoce varios formatos de

datos se puede especificar el tipo para que busque el archivo de interes.
Como se puede observar se pueden abrir los archivos generados por
ProfileR. (Figura. 3.9)

Figura 3.9. Menu para abrir archivos de datos. Se muestran los diferentes formatos que reconoce
ZONDRES2D.

49
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

3. Despues de seleccionar el archivo aparecer una pantalla con las

caracteristicas para formar la malla de datos, se pueden modificar o dejar
los que aparecen por default y luego seleccionar Apply. (Figura. 3.10).

Figura 3.10. Menu para el constructor de malla. Se muestran los diferentes parmetros que
se pueden modificar para construir la malla.

4. El programa mostrar las pseudoseccin de los datos, la pseudoseccin

para el modelo (que en un principio ser una restividad uniforme) y el
modelo en si. (Figura 3.11).

Figura 3.11. Despligue de los datos. (a) Muestra la pseudoseccin de resistividad aparente
de los datos observados. (b) Muestra la pseudoseccin con los datos de resistividad
aparente calculados. (c) Muestra la seccin en bloques de resistividad.

50
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

5. En Options seleccionar Program setup y ah se puede elegir las

diferentes opciones para crear un modelo directo o hacer una inversin de
datos y todos los paramtros que se pueden modificar para obtener
resultados ms afines a nuestro objetivo. Una vez modificados los
parmetros de inters seleccionar Apply (Figura 3.12).

Figura 3.12. Men para cambiar los parametro o metodos de inversion y modelado directo.

6. Para inversin de click en el icono de inversin. (Figura 3.13)

Figura 3.13. Se muestra cual es el icono para iniciar el proceso de inversin.

51
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

7. Usted ver una imagen que aparece junto con la pseudoseccin para ese
modelo. (Figura 3.14)

Figura 3.14. Resultados de la inversin. (a) Muestra la pseudoseccin de resistividad

aparente de los datos observados. (b) Muestra la pseudoseccin con los datos de
resistividad aparente calculados. (c) Muestra la seccin en bloques de resistividad con los
resultados de la inversin.

8. Selecione Options/Model/Contour Section y podrs ver un mapa de

contornos. (Figura 3.15)

Figura 3.15. En la imagen inferior de la figura se muestra la diferencia al seleccionar un

mapa de contornos en contraste con la seccin en bloques.

52
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

9. De doble click en la barra de colores y seleccione setup y podra modificar

la escala de colores, mostrar las isolineas, modificar el nmero, etcetera, si
as lo deseas. (Figura 3.16)

Figura 3.16. Se muestra el menu de la escala de colores, en el cual se pueden modificar

algunos otros parametros.

10. Podr obtener una seccin e interpretarla (Figura 3.17)

Figura 3.17. Se muestra la seccin que se obtiene despues de modicar todos los
parametros necesarios para que el modelo concuerde con nuestro objetivo.

53
Cdigo de Inversin de Datos de Tomografa de Resistividad Elctrica ZONDRES2D Captulo 3

11. Por ultimo en File/Save puede guardar la inversin de los datos en

diferentes formatos para que si lo desea los pueda abrir con otros
programas, guardar los datos de la inversin para abrirlos posteriormente o
simplemente guardar como imagen la seccin obtenida. (Figura. 3.18)

Figura 3.18. Se muestran los diferentes formatos en los que se pueden guardar los datos
obtenido en la inversin para as poderlos abrir en algun otro programa.

54
 Aplicacin en Modelos Sintticos Captulo 4

CAPTULO 4. APLICACION EN MODELOS SINTTICOS

4.1. INTRODUCCIN

En este captulo se presentan las ventajas y desventajas de emplear los diferentes

tipos de arreglos de electrodos con un anlisis de sensitividad. Tambin se
estudiar el poder resolutivo del programa ZONDRES2D con base a diversos
modelos sintticos que asemejan estructuras geolgicas conocidas. Este anlisis
se hace con cinco arreglos diferentes de electrodos, el propsito es hacer una
recreacin de las estructuras geolgicas reales en subsuelo y estudiar el efecto
que tendran con un estudio de resistividad. Esto permitir evaluar la respuesta de
las diferentes modalidades de electrodos para la obtencin de la tomografa, lo
cual adems es recomendable hacer antes de llevar a cabo cualquier estudio de
resistividad. Antes de llevar a cabo un estudio sobre el terreno se debe tener
alguna informacin acerca de la forma y el tamao de los objetivos previstos. Al
tratar con diferentes arreglos en la pantalla digital de la computadora, se puede
evitar el uso de arreglos no adecuados para la deteccin de estructuras de inters.
Tambin se puede tener una idea del espaciamiento adecuado entre los
electrodos adyacentes a utilizar, y la mxima separacin de electrodos necesaria.

4.2. MODELADO DIRECTO 2D

En el problema de modelado directo, se especfica la distribucin de resistividades

en el subsuelo, siendo el objetivo calcular los valores de resistividad aparente por
medio del programa ZONDRES2D. La subrutina de modelado directo, es tambin
parte integral de cualquier programa de inversin, ya que es necesario calcular los
valores de resistividad aparente en cada paso iterativo. Existen tres mtodos
principales para calcular los valores de resistividad aparente para un modelo
especfico. Ellos son (i) los mtodos analticos, (ii) mtodos de elementos en la
frontera y (iii) los mtodos de diferencias finitas y de elementos finitos. Los
mtodos analticos son probablemente los mtodos ms exactos, pero estn
restringidos a geometras relativamente simples (como una esfera o cilindro). Los
mtodos de elementos en las fronteras son ms flexibles, pero el nmero de
regiones con diferentes valores de resistividad que se permite es un tanto limitada
(normalmente menos de 10). En estudios de ingeniera y ambientales el subsuelo
por lo general tiene una distribucin de resistividad arbitraria, de modo que los
mtodos de diferencia finita y elementos finitos son por lo general la nica eleccin
viable. Estos mtodos subdividen el subsuelo en miles de celdas a las cuales se
les asignan diferentes valores de resistividad. Sin embargo, los mtodos analticos
y elementos de frontera usualmente son mtodos independientes que se utilizan
para comprobar la exactitud de los mtodos de diferencia finita y de elemento
finito. Los modelos sintticos que se ilustran en este captulo se generaron
mediante los programas PROFILER V2 y ProfileR. En este ltimo programa se
puede elegir el mtodo de diferencias finitas o el de elementos finitos. El subsuelo
se divide en una rejilla con un gran nmero de pequeas celdas rectangulares
(Figura 2.1) y posteriormente se puede especificar el valor de resistividad de cada

55
 Aplicacin en Modelos Sintticos Captulo 4

una de las celdas. El mtodo de diferencias finitas se basa en un mtodo descrito

por Dey y Morrison (1979), pero con una modificacin introducida por Loke (1994)
para corregir una inconsistencia menor por el mtodo de discretizacin de zona. El
mtodo de elementos finitos utiliza elementos triangulares estndar de primer
orden (Silvester y Ferrari, 1990).

El programa requiere que se suministren los valores de resistividad de los bloques

rectangulares entre las lneas de la malla (y otra informacin). En general se usa
un modelo de resistividades razonable, es decir, que el modelo no est
demasiado pequeo o grande. El programa genera una pseudoseccin de
resistividades aparentes con todos sus datos que son los que posteriormente se
utilizaron con el programa de inversin ZONDRES2D que nos genera la seccin
de resistividades verdaderas. Aunque nuestro principal inters es la inversin de
datos de campo, el programa de modelado directo ProfileR es til tambin, sobre
todo en la etapa de planificacin de estudios superficiales. Con este programa se
pueden calcular los patrones de sensitividad para los diferentes arreglos con la
finalidad de tener una idea de la respuesta que dar cierto tipo de arreglo para una
estructura en particular (por ejemplo una zona de fractura vertical).

4.3. VENTAJAS Y DESVENTAJAS DE LAS DIFERENTES MODALIDADES DE

ELECTRODOS.

Este anlisis basado en el programa ZONDRES2D demuestra que de acuerdo a

las diferentes modalidades de electrodos la tomografa es sensitiva a diferentes
grados sobre una misma estructura. La eleccin de la mejor modalidad para un
estudio de campo depende del tipo de estructura que se pretende encontrar, la
sensitividad del resistivimetro y el nivel de ruido ambiental. En la prctica, las
modalidades que se emplean frecuentemente para estudios de imgenes
elctricas 2D son: Wenner, Wenner-Schlumberger, dipolo-dipolo, polo-polo y polo-
dipolo. En cuanto a las caractersticas que deben considerarse para cada una de
las diferentes modalidades son: 1) La sensitividad de la modalidad a los cambios
vertical y horizontal de la resistividad del subsuelo, 2) La profundidad de
investigacin, 3) La cobertura horizontal de los datos y 4) La intensidad de la
seal.

4.3.1. Modalidad Wenner.

Esta es una modalidad robusta, la cual fue popularizada por el trabajo pionero
hecho por un grupo de investigadores de la Universidad de Birmingham (Griffiths y
Turnbull, 1985; Griffiths, et al.,1990). En la Figura 4.1 se puede observar que la
grfica de sensitividad para esta modalidad tiene contornos casi horizontales bajo
el centro del arreglo, por ello la modalidad Wenner es relativamente sensitiva a
cambios verticales en la resistividad del subsuelo bajo el centro del arreglo. Sin
embargo, es menos sensitivo a cambios horizontales a la resistividad del
subsuelo. En general, el arreglo Wenner es bueno para resolver cambios
verticales, pero relativamente pobre en la deteccin de cambios horizontales. La
intensidad de su seal es inversamente proporcional al factor geomtrico utilizado
56
 Aplicacin en Modelos Sintticos Captulo 4

para calcular el valor de resistividad aparente para la modalidad (Figura 1.3). Para
la modalidad Wenner, el factor geomtrico es 2a, el cual es muy pequeo
comparado con las otras modalidades. En cuanto a otras modalidades comunes,
la modalidad Wenner tiene una intensidad de seal demasiado fuerte. Esto puede
ser un factor importante si el estudio es hecho en reas con alto ruido ambiental.
Una desventaja de esta modalidad para estudios 2D es la cobertura horizontal
relativamente pobre a medida que el espaciamiento entre electrodos es mayor.
Esto podra ser un problema si se utiliza un sistema con un nmero pequeo de
electrodos.

Figura 4.1. Grfica de la funcin de sensitividad de la modalidad Wenner

4.3.2. Modalidad Wenner-Schlumberger

Este es un hbrido entre las modalidades Wenner y Schlumberger (Pazdirek y

Blaha, 1996). Esta modalidad puede utilizarse sobre un sistema con un arreglo de
electrodos de espaciamiento constante (Figura1.2). Ntese que el factor n para
esta modalidad es el radio de la distancia entre los electrodos C1-P1 (o C2-P2) al
espaciamiento entre el par de electrodos de potencial P1-P2, los cuales tienen una
distancia fija a. En la Figura 4.2 se puede apreciar que la grfica de sensitividad
para la modalidad Wenner-Schlumberger es ligeramente diferente a la de la
modalidad Wenner con una ligera curvatura vertical bajo el centro del arreglo, y
valores de sensitividad ligeramente bajos en las regiones comprendidas entre los
electrodos C1 y P1 (tambin C2 y P2). Existe una concentracin ligeramente
grande de valores de alta sensitividad bajo los electrodos P1-P2. Esto nos indica
que esta modalidad es moderadamente sensitiva tanto a estructuras geolgicas
horizontales como verticales. En reas donde existen ambos tipos de estructuras
geolgicas, esta modalidad podra ser un buen compromiso entre las modalidades
Wenner y dipolo-dipolo. La intensidad de la seal para esta modalidad es mucho
ms pequea que la de la modalidad Wenner, pero ms alta que la de la
modalidad dipolo-dipolo. La modalidad Wenner-Schlumberger tiene una cobertura
horizontal ligeramente mejor comparada con la modalidad Wenner, pero mucho
ms reducida que la que se obtiene con la modalidad Dipolo-Dipolo.

57
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.2. Grfica de la funcin de sensitividad de la modalidad Wenner-Schlumberger

4.3.3. Modalidad Dipolo- Dipolo.

El arreglo de electrodos de esta modalidad se puede observar en la Figura 1.5. En

la Figura 4.3 se muestra la grfica de la funcin de sensitividad la cual muestra
que los valores de sensitividad mas grandes estn localizados entre el par dipolar
C1-C2, as como tambin entre el par P1-P2. Esto indica que esta modalidad es
ms sensitiva a cambios de resistividad entre los electrodos de cada par dipolar.
Ntese que la sensitividad en los patrones de los contornos es casi vertical. De
esta manera, se puede decir que la modalidad dipolo-dipolo es muy sensitiva a
cambios horizontales de resistividad pero relativamente insensitivo a cambios
verticales de resistividad. Lo que indica que es muy bueno para mapear
estructuras verticales, pero relativamente pobre mapeando estructuras
horizontales tales como sills o estratos sedimentarios. Una posible desventaja de
esta modalidad es que la intensidad de la seal es muy pequea para valores
grandes del factor n. El voltaje es inversamente proporcional al cubo del factor n,
esto nos indica que para la misma corriente, el voltaje medido por el resistivimetro
decae aproximadamente 200 veces cuando incrementa el factor n de 1 a 6.
Para emplear esta modalidad de una manera ms efectiva, el resistivimetro debe
tener una alta sensitividad y muy buen rechazo de ruido y adems, debe existir un
buen contacto entre los electrodos y el terreno. Esta modalidad ha sido una de las
ms utilizadas para detectar estructuras tales como cavidades, donde la buena
resolucin horizontal de esta modalidad tiene una ventaja mayor. Para la
modalidad dipolo-dipolo, las regiones con valores de sensitividad altos estn
concentradas bajo el par de electrodos C1-C2 y bajo el par de electrodos P1-P2.

58
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.3. Grfica de la funcin de sensitividad de la modalidad Dipolo-Dipolo.

4.3.4. Modalidad Polo-Polo

En la prctica una modalidad polo-polo ideal, con un solo electrodo de corriente y

uno de potencial no existe. Para aproximarnos a una modalidad polo-polo, el
segundo electrodo de corriente y de potencial (C2 yP2) deben llevarse a una
distancia la cual debe ser 20 veces la separacin mxima entre los electrodos C1
y P1 empleados en el estudio. El efecto del electrodo C2 (y similarmente para el
P2) es aproximadamente proporcional al radio de la distancia C1-P1 a la distancia
C2-P1. Si los efectos de los electrodos C2 y P2 no son tomados en cuenta, la
distancia de estos electrodos a la lnea de estudio debe ser al menos 20 veces el
espaciamiento empleado ms grande entre C1-P1. As se asegura que el error
sea menor al 5%. Otra desventaja de esta modalidad es que a causa de la
distancia tan grande que existe entre los electrodos P1 y P2, esto puede producir
una gran cantidad de ruido telrico el cual puede degradar severamente la calidad
de las mediciones. De esta manera esta modalidad se usa principalmente para
estudios donde el espaciamiento entre electrodos sea relativamente pequeo
(menos de 10 metros). Es popular en algunas aplicaciones como estudios
arqueolgicos, para los cuales se requiere un espaciamiento entre electrodos
pequeo. Esta modalidad tiene la cobertura horizontal ms extensa y la mayor
profundidad de investigacin en comparacin con las modalidades anteriores, sin
embargo, tiene la resolucin ms pobre, la cual se ve reflejada por el gran
espaciamiento entre los contornos en la grafica de la funcin de sensitividad
(Figura 4.4)

59
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.4. Grfica de la funcin de sensitividad de la modalidad Polo-Polo.

4.3.5. Modalidad Polo-Dipolo

Esta modalidad tambin tiene una cobertura horizontal relativamente buena, pero
tiene una intensidad de seal relativamente alta en comparacin con la modalidad
dipolo-dipolo y no es tan sensitivo al ruido telrico como la modalidad polo-polo. A
diferencia de las otras modalidades, la modalidad polo-dipolo es una modalidad
asimtrica (Figura 1.4) y sobre estructuras simtricas las anomalas de resistividad
aparente en la pseudoseccin son asimtricas. En algunas situaciones, la
asimetra en los valores de la resistividad aparente medidas podran influenciar en
el modelo obtenido despus de aplicar la inversin.

Esta modalidad tambin requiere de un electrodo remoto, el electrodo C2, el cual

debe estar localizado lo suficientemente lejos de la lnea de estudio. Para esta
modalidad, el efecto del electrodo C2 es aproximadamente al cuadrado del radio
de la distancia C1-P1 a la distancia C2-P1. De esta manera la modalidad polo-
dipolo es menos afectada por el electrodo remoto C2 comparada con la modalidad
polo-polo. Si la distancia del electrodo C2 es ms de 5 veces la distancia ms
grande C1-P1 usada, el error causado por descuidar el efecto del electrodo C2 es
menor al 5%. Debido a su buena cobertura horizontal, esta es una modalidad
atractiva para sistemas multielectrodos con un nmero relativamente pequeo de
nodos. La intensidad de la seal es baja comparada con las modalidades Wenner
y Wenner-Schlumberger (Figura 4.5), pero alto con respecto a la modalidad dipolo-
dipolo. Esta intensidad de la seal decrece con el cuadrado del factor n. Mientras
este efecto no es tan severo como en la modalidad dipolo-dipolo, usualmente no
es recomendable usar valores del factor n mayores de 8 a 10. Ms all de esto, el
espaciamiento a entre el par dipolar P1-P2 se incrementa para obtener una
intensidad de seal ms fuerte.

60
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.5. Grfica de la funcin de sensitividad de la modalidad Polo-Dipolo.

4.4. DESCRIPCIN DE LOS MODELOS SNTETICOS.

Los modelos utilizados semejan diversos tipos de estructuras sencillas que se

pueden encontrar en el subsuelo. Para todos los modelos sintticos se utilizaron 6
niveles y 30 electrodos con un espaciamiento entre electrodos de 3 metros,
cubriendo una distancia horizontal de 87 metros. El espaciamiento nos reflejar
profundidades medias aproximadamente de 3.2 metros para la modalidad Wenner,
7.4 metros en el caso de Wenner-Schlumberger, 5.2 metros para la modalidad
Dipolo-Dipolo, 2.6 metros para la modalidad Polo-Polo y 7.4 metros para la
modalidad Polo-Dipolo. (Ver Tabla 3, Captulo 1).

Todos los modelos estn contaminados con un 2% de ruido aleatorio, por lo que
para la inversin de las pseudosecciones de resistividad aparente se utiliz un
factor de amortiguamiento de 0.1, el cual resulto ser el ms ptimo para la
inversin.

4.4.1. Modelo de un bloque rectangular.

Este modelo consta de un bloque rectangular con un valor de resistividad de 500

Ohm-m y la resistividad del medio es de 100 Ohm-m. (Figura 4.6.)

Figura 4.6. Modelo de un bloque rectangular horizontal

61
 Aplicacin en Modelos Sintticos Captulo 4

En las figuras 4.7 y 4.8 se presentan las secciones de resistividad verdadera, sin
agregar ruido aleatorio, obtenidas con las modalidades Wenner, Wenner-
Schlumberger, Dipolo-Dipolo, Polo-Polo y Polo-Dipolo, de un bloque rectangular.

a) Modalidad Wenner

b) Modalidad Wenner-Schlumberger

c) Modalidad Dipolo-Dipolo

d) Modalidad Polo-Polo

Figura 4.7 Secciones de resistividad obtenidas de un modelo sinttico de un bloque rectangular

para las modalidades a) Wenner, b) Wenner-Schlumberger c) Dipolo-Dipolo y d) Polo-Polo.

62
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.8. Seccin de resistividad obtenida de un modelo sinttico de un bloque rectangular para
la modalidad Polo-Dipolo.

Como se puede observar en las imgenes elctricas, las cinco modalidades dan
una buena resolucin del cuerpo, los modelos definen con claridad una zona de
alta resistividad que coincide con la ubicacin del bloque rectangular. En todos los
casos se empleo un factor de amortiguamiento de 0.1, el cual nos ayudo a definir
bien los cuerpos.

Por otra parte se aprecia que aunque todas las imgenes elctricas muestran un
cuerpo de alta resistividad, la mayora de los arreglos definen mejor la estructura,
sin embargo para el caso de la modalidad Dipolo-Dipolo la estructura se ve como
si fuera vertical, una posible causa se puede deber a que como se observ en el
punto anterior de este captulo, dicha modalidad no es recomendable para mapear
estructuras horizontales.

4.4.2. Modelo de dos bloques horizontales rectangulares.

Este modelo consiste de dos bloques rectangulares horizontales con valores de

resistividad de 10 Ohm-m (izquierda) y 500 Ohm-m (derecha) en un medio de 100
ohm-metro. (Figura 4.9)

Figura 4.9. Modelo de dos bloques horizontales rectangulares.

63
 Aplicacin en Modelos Sintticos Captulo 4

En las figuras 4.10 y 4.11 se muestran las secciones de resistividad verdadera, sin
agregar ruido aleatorio, obtenidas con las modalidades Wenner, Wenner-
Schlumberger, Dipolo-Dipolo, Polo-Polo y Polo-Dipolo.

a) Modalidad Wenner

b) Modalidad Wenner-Schlumberger

c) Modalidad Dipolo-Dipolo

d) Modalidad Polo-Polo

Figura 4.10. Secciones de resistividad obtenidas de un modelo sinttico de dos bloques

rectangulares horizontales para las modalidades a) Wenner, b) Wenner-Schlumberger, c) Dipolo-
Dipolo y d) Polo-Polo.

64
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.11. Seccin de resistividad obtenida de un modelo sinttico de dos bloques rectangulares
horizontales para la modalidad Polo-Dipolo.

Como se puede observar en las imgenes elctricas de todas las modalidades

aparecen los dos cuerpos separados correctamente. Los modelos muestran una
zona de baja resistividad que coincide con el bloque rectangular ubicado a la
izquierda de la seccin y otra zona de alta resistividad que coincide con el bloque
ubicado a la derecha de la seccin. Al igual que en los modelos anteriores se
empleo un factor de amortiguamiento de 0.1, el cual como se puede observar,
ayudo a delinear claramente ambos cuerpos en la imagen elctrica.

Este modelo se gener para evaluar la respuesta de cada modalidad a los

cambios horizontales de resistividad. Al igual que en el ejemplo anterior la
modalidad Dipolo-Dipolo es la ms pobre en cuanto a definir los cuerpos, y
tambin podemos observar que las que mejor los definen son las modalidades
Wenner y Wenner-Schlumberger.

4.4.3. Modelo de cuatro bloques horizontales rectangulares.

Este modelo consiste de cuatro bloques rectangulares horizontales con

resistividades de 100 Ohm-m en un medio de 10 Ohm-m. Tres de dichos bloques
estn situados a la misma profundidad y uno de ellos (el tercero de izquierda a
derecha) est situado a mayor profundidad. (Figura 4.12)

Figura 4.12. Modelo de cuatro bloques horizontales rectangulares.

65
 Aplicacin en Modelos Sintticos Captulo 4

En las figuras 4.13 y 4.14 se pueden observar las secciones de resistividad

verdadera, sin agregar ruido aleatorio, obtenidas con las modalidades Wenner,
Wenner-Schlumberger, Dipolo-Dipolo, Polo-Polo y Polo-Dipolo.

a) Modalidad Wenner

b) Modalidad Wenner-Schlumberger

c) Modalidad Dipolo-Dipolo

d) Modalidad Polo-Polo

Figura 4.13. Secciones de resistividad obtenidas de un modelo sinttico de cuatro bloques

rectangulares horizontales para las modalidades a) Wenner, b) Wenner-Schlumberger, c) Dipolo-
Dipolo y d) Polo-Polo.

66
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.14. Seccin de resistividad obtenida de un modelo sinttico de cuatro bloques

rectangulares horizontales para la modalidad Polo-Dipolo.

Como se puede observar en las imgenes elctricas en todas las modalidades

aparecen los cuatro cuerpos separados correctamente, casi todos los cuerpos se
aprecian con la misma resolucin y en todas las imgenes se aprecia la zona de
baja resistividad y las cuatro zonas de alta resistividad que representan los
cuerpos.

Al igual que en el ejemplo anterior este modelo se diseo con la finalidad de

observar la respuesta de cada modalidad a cambios laterales de resistividad. En
este caso podemos observar que la modalidad que define mejor los cuerpos es la
modalidad Polo-Polo, esto puede ser debido a que como vimos en la seccin 4.3
de este captulo es la modalidad que posee la cobertura horizontal ms extensa en
comparacin con las otras modalidades.

4.4.4. Modelo de dos capas y un dique.

Este modelo consiste de una capa horizontal que asemeja un suelo con un valor
de resistividad de 20 Ohm-m, un dique vertical con resistividad de 1000 Ohm-m,
localizados en un medio de 100 Ohm-m. (Figura 4.15)

Figura 4.15. Modelo de dos capas y un dique.

67
 Aplicacin en Modelos Sintticos Captulo 4

En las figuras 4.16 y 4.17 se pueden observar las secciones de resistividad

verdadera, sin agregar ruido aleatorio, obtenidas con las modalidades Wenner,
Wenner-Schlumberger, Dipolo-Dipolo, Polo-Polo y Polo-Dipolo.

a) Modalidad Wenner

b) Modalidad Wenner-Schlumberger

c) Modalidad Dipolo-Dipolo

d) Modalidad Polo-Polo

Figura 4.16. Secciones de resistividad obtenidas de un modelo sinttico de dos capas y un dique
para las modalidades a) Wenner, b) Wenner-Schlumberger, c) Dipolo-Dipolo y d) Polo-Polo.

68
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.17. Seccin de resistividad obtenida de un modelo sinttico de dos capas y un dique para
la modalidad Polo-Dipolo.

En este caso, este modelo se diseo para observar la respuesta de los arreglos a
cambios verticales de resistividad. Podemos observar en las imgenes elctricas
un cuerpo definido por una zona de alta resistividad (dique) y una capa
representada por una zona de baja resistividad, sin embargo podemos ver que no
todas las modalidades definen bien la estructura y que a diferencia de los
ejemplos anteriores, en este caso la modalidad Dipolo-Dipolo parece ser una de
las que mejor definen la estructura, seguida de las modalidades Wenner-
Schlumberger y Polo-Dipolo.

4.4.5. Modelo de una falla y un bloque horizontal.

Este modelo consiste de dos bloques rectangulares horizontales de diferentes

dimensiones con valores de resistividad de 30 Ohm-m (izquierda) y 10 Ohm-m
(derecha), situados en un medio de 100 Ohm-m. (Figura 4.18)

Figura 4.18. Modelo de una falla y un bloque horizontal.

Las figuras 4.19 y 4.20 muestran las secciones de resistividad verdadera, sin
agregar ruido aleatorio, obtenidas con las modalidades Wenner, Wenner-
Schlumberger, Dipolo-Dipolo, Polo-Polo y Polo-Dipolo.

69
 Aplicacin en Modelos Sintticos Captulo 4

a) Modalidad Wenner

b) Modalidad Wenner-Schlumberger

c) Modalidad Dipolo-Dipolo

d) Modalidad Polo-Polo

Figura 4.19 Secciones de resistividad obtenidas de un modelo sinttico de una falla y un bloque
horizontal para las modalidades a) Wenner, b) Wenner-Schlumberger, c) Dipolo-Dipolo y d) Polo-
Polo.

70
 Aplicacin en Modelos Sintticos Captulo 4

Figura 4.20. Seccin de resistividad obtenida de un modelo sinttico de dos capas y un dique para
la modalidad Polo-Dipolo.

Se puede observar en las imgenes elctricas de todas las modalidades que

aparecen los dos cuerpos separados correctamente representados por zonas de
baja resistividad y tambin se refleja que uno de los cuerpos es mucho ms
grande en dimensiones (falla) que el otro (bloque). Podemos observar que las
modalidades que definen las formas de manera ms correcta son Dipolo-Dipolo,
Wenner-Schlumberger y Polo-Dipolo.

En general, se puede apreciar en los diferentes modelos propuestos, en cualquier

modalidad, que los cuerpos son determinados, sin embargo, depende de cada
modalidad que se puedan definir mejor ciertas estructuras.

71
 Aplicacin en Datos de Campo Captulo 5

CAPTULO 5. APLICACIN EN DATOS DE CAMPO

5.1. INTRODUCCIN.

Actualmente, uno de los campos ms importantes en los cuales se aplica la

prospeccin elctrica es el estudio de problemas ambientales, aunque tambin se
utiliza para la caracterizacin y localizacin de estructuras anmalas como pueden
ser cavidades o antiguos tneles y delineacin de sistemas de acuferos. En este
captulo se aplica la tcnica de Tomografa de Resistividad Elctrica 2D en cinco
zonas diferentes: delegacin lvaro Obregn (Ciudad de Mxico); Coeneo e Ixtln
de los Hervores (Michoacn), el Pahualt (Espinal, Veracruz) y por ltimo Valle Etla,
Oaxaca. Todos los datos obtenidos se procesaron mediante el programa de
inversin ZONDRES2D, para poder apreciar su resolucin, determinar sus
bondades y efectividad, se comparan estas mismas imgenes con las generadas
con el programa RES2DINV que es un programa probado, reconocido y
recomendado para procesar e interpretar datos de resistividad.

5.2. INSTRUMENTACIN

A continuacin se presentan las herramientas necesarias que utilizaron para la

realizacin de cada estudio.

1) 20 a 30 barras metlicas como electrodos.

2) 5 marros.

3) Computadora porttil con software, ZONDRES2D.

4) Cable multi-electrodos

5) Receptor IPR-10A y transmisor IPC-7 2.5 Kw marca SCINTREX, para la

adquisicin de los datos de la Delegacin lvaro Obregn y para obtener los datos
de SEV en la zona de Pahual y Valle Etla y Resistivimetro automatizado Sting R1-
IP, en la adquisicin de los datos de Delegacin lvaro Obregn, Coeneo, Ixtln y
Pahual.

5.3. PROCESAMIENTO DE LOS DATOS

Los datos obtenidos en campo se trataron de manera independiente con dos

programas de inversin, RES2DINV y ZONDRES2D. Los mtodos de inversin
que se aplicaron a todos los perfiles en este trabajo fueron con el mtodo de
mnimos cuadrados con condicionamiento suavizado, debido a que este mtodo
minimiza las variaciones en los parmetros del modelo y en combinacin con un
modelo robusto produce imgenes elctricas de buena resolucin, este mtodo se
describi en el captulo 2.

72
 Aplicacin en Datos de Campo Captulo 5

A todas las secciones de resistividad aparente, los procesos iterativos de inversin

que se aplicaron fueron iguales, para el caso del programa RES2DINV se
realizaron 3 iteraciones y 10 para el caso de ZONDRES2D.

5.4. DELEGACIN LVARO OBREGN.

La Ciudad de Mxico est localizada en el Cinturn Volcnico Trans- Mexicano

(CVTM), en el centro de Mxico. El CVTM es una provincia Plio-Cuaternaria calco-
alcalina que atraviesa Mxico de oeste a este. La regin comprende la mayora del
vulcanismo histrico y actual de Mxico e incluye estratovolcanes andesticos-
dasticos, de vulcanismo rioltico y mayores centros riolticos. La ciudad de Mxico
se localiza en una cuenca a una altura de 2200 metros s.n.m. y aproximadamente
a 300 kilmetros al este de la Trinchera de Amrica Central. Las Sierras
Volcnicas que rodean la cuenca son: la Sierra Nevada al este, la Sierra de
Chichinautzin al sur, la Sierra de las Cruces al oeste y la Sierra de Guadalupe al
norte. La naturaleza de este valle intramontaa afecta directamente la calidad del
aire, el suministro de agua y el desarrollo urbano. Los diferentes procesos
humanos definen el desarrollo urbano e industrial y tambin la poblacin influye en
la calidad de vida. Por lo tanto, administrar los servicios para los habitantes de
esta ciudad es un gran reto (Campos et al. 1197).

Adems la insuficiencia de casas habitacin para la poblacin, un problema que

se ha ido incrementando durante los pasados 40 aos. A finales de los aos
cuarenta y cincuenta, y a causa del crecimiento anrquico de las areas urbanas,
la mayora del material utilizado con propsitos de construccin se obtuvo por
medio de minas excavadas en las sierras del oeste. Dichos materiales fueron
empleados para la construccin de casas, departamentos y edificios del gobierno
en el centro de la Ciudad de Mxico. Cuando esta industria decay a principios de
los aos sesenta, la mayora de estas minas, terminaron en cavidades que fueron
abandonadas o rellenadas con desperdicios. Debido a su tamao, posteriormente
estas cavidades fueron utilizadas por familias enteras como refugios, o tenan
otros usos por la gente del vecindario. Las colonias de esta delegacin, se han ido
asentando sobre zonas minadas en los pasados 20 aos. Actualmente, la
localizacin de algunas de estas estructuras es desconocida. Han ocurrido
muchos accidentes desde entonces; las cavernas se han colapsado y hogares
completos literalmente han sido tragados por el terreno.

La universidad Nacional Autnoma de Mxico y la delegacin lvaro Obregn

firmaron en 1995 un programa de investigacin para localizar zonas de alto riesgo,
en particular zonas minadas sobre las cuales existen muchos asentamientos. En
este trabajo se presentan resultados en un rea con mediciones de resistividad.

73
 Aplicacin en Datos de Campo Captulo 5

5.4.1. Marco Geolgico.

La Ciudad de Mxico fue fundada originalmente sobre una pequea isla en un

lago rodeado por Sierras Volcnicas. Este se localiz en la porcin central del
CVTM (Pasquare et al., 1987). El rea estudiada se encuentra al suroeste del
Valle de Mxico (Figura 5.1). El relieve es abrupto y compuesto de una serie de
colinas y barrancas profundas. Las rocas superficiales son depsitos de materiales
gneos arrojados por volcanes ubicados dentro de la Sierra de las Cruces. Estos
eventos geolgicos produjeron avalanchas de material piroclstico, tobas y
brechas que se depositaron en las faldas de las Sierras. Las rocas volcnicas que
conforman el relieve del rea de estudio son de la era Cenozoica, y es conocida
como la formacin Tarango. Este Horizonte geolgico es una secuencia
compuesta de losas de arena, piedras pmez interestratificadas con arcillas y
calizas. La piedra pmez se produjo de violentas explosiones volcnicas, algunos
fragmentos fueron arrastrados por el viento a grandes distancias del crter y
depositados en las zonas bajas. La columna estratigrfica de la formacin Tarango
se divide en dos unidades. La Unidad A es la ms joven y es encontrada en las
pendientes de las barrancas y consiste de fragmentos de rocas arrastrados por
corrientes de agua. Estos deslizamientos fueron producidos por la saturacin del
agua de lluvia, movindolos hacia abajo a grandes velocidades por efectos de
gravedad. Estos tipos de rocas estn bien compactadas y hospedadas en una
matriz de arena con esquinas angulosas. La Unidad B es una secuencia
piroclstica, un producto de erupciones volcnicas. La piedra pmez, travertino,
arena y material aluvial son los tipos de rocas que forman esta unidad. Dichas
rocas fueron empleadas para elaborar ladrillos ligeros. Esta unidad forma la
corona o partes ms altas de las colinas de esta rea. La mayora de la industria
minera explot las secciones ms bajas en las arenas pumiciticas. Las minas
principales estn localizadas aqu, y la posicin estratigrfica define calidad y
dureza de los materiales. Las cavidades encontradas en esta unidad tienen un
rango de profundidad de 5, 10 y 25 metros. El tamao de las cmaras principales
es variable, del orden de 2 a 3 metros y hasta 30 metros. Las longitudes de los
tneles tienen cientos de metros.

5.4.2. Ubicacin Geogrfica.

La delegacin lvaro Obregn se localiza al poniente del Distrito Federal (Figura

5.1), extendindose de noreste a sureste. Limita al norte con la delegacin Miguel
Hidalgo, al oriente con Benito Jurez y Coyoacn, al sur con Magdalena
Contreras, Tlalpan y con el municipio de Jalatlaco en el Estado de Mxico y al
poniente con Cuajimalpa. En esta delegacin se encuentra el acceso poniente de

74
 Aplicacin en Datos de Campo Captulo 5

la ciudad, sus vialidades regionales carretera federal y autopista, constituyen la

entrada de mercanca y poblacin de los Estados de Mxico y Michoacn.

Figura 5.1 Mapa de ubicacin de la delegacin lvaro Obregn, Ciudad de Mxico.

Geogrficamente la delegacin lvaro Obregn est situada entre los paralelos 19; 14 N y 19;
25'S y los meridianos 99; 10'E W y 99; 20'O W.

5.4.3. Perfiles de Imgenes Elctricas.

Se realiz un perfil geoelctrico 2D con la modalidad Wenner-Schlumberger, en la

colonia Ampliacin Presidentes y otro con la modalidad dipolo-dipolo en la colonia
Arvide, localizadas en la delegacin lvaro Obregn, Cuidad de Mxico. En el
rea de estudio existen varios asentamientos humanos situados sobre zonas
minadas, la cual se puede considerar como zona de alto riesgo.

[Link]. Calle Pedro Aguirre (Colonia ampliacin presidentes).

Este es un ejemplo interesante de un estudio bidimensional con la modalidad

Wenner- Schlumberger con la finalidad del mapeo de posibles cavidades. Para la
construccin de la seccin de resistividades aparentes se emplearon 30 electrodos
con una separacin de 5 metros y posteriormente para cubrir una distancia total
del perfil de 225 metros se hicieron 4 corrimientos de 20 metros cada uno, con 11
niveles, llegando a una profundidad de investigacin aproximada de 21 metros. El
total de datos de resistividad fue de 363.

75
 Aplicacin en Datos de Campo Captulo 5

Resultados

En la Figura 5.2 se puede observar que en ambas secciones de resistividad se

define una posible cavidad, delimitada por una zona de alta resistividad de forma
aproximadamente circular, bajo el centro de la lnea de estudio aproximadamente
a 120 metros en la direccin horizontal y a una profundidad de 15 metros
aproximadamente. Tambin se puede observar una zona anmala de baja
resistividad en la parte superior izquierda a una distancia de 60 metros
aproximadamente. El medio encajonante tiene un valor de resistividad de 60 ohm-
m.

Figura 5.2. Modelos de inversin de una cavidad encontrada en el perfil realizado en la calle Pedro
Aguirre, Colonia Ampliacin Presidentes, Delegacin lvaro Obregn. (a) Modelo de Inversin
realizado con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa
ZONDRES2D.

[Link]. Calle Segunda Cerrada de Minas (Colonia Arvide)

Para la construccin de este perfil se utiliz el arreglo dipolo-dipolo con 28

electrodos equiespaciados 5 metros, cubriendo una distancia horizontal de 135
metros y alcanzando una profundidad aproximada de 30 metros. El total de
lecturas de resistividad adquiridas fue de 213.

Resultados

En ambos modelos de inversin (Figura 5.3) se pueden apreciar tres zonas

anmalas que se pueden asociar a posibles cavidades, una de estas zonas
anmalas representada por una zona de alta resistividad bajo el centro de la lnea
de estudio aproximadamente a unos 90 metros en la direccin horizontal y a una

76
 Aplicacin en Datos de Campo Captulo 5

profundidad aproximada de 10 metros, resulto ser una cavidad vaca. Las otras
dos cavidades estn definidas por dos zonas de resistividad media (alrededor de
200 ohm-m) aproximadamente a 20 y 50 metros en la direccin horizontal de la
lnea de estudio y profundidades de 5 y 10 metros respectivamente. Cabe
mencionar que estas cavidades se muestran con otro valor de resistividad debido
a que se encuentran rellenas.

Figura 5.3. Modelos de inversin de cavidades encontradas en el perfil realizado en la calle

Segunda Cerrada de Minas de la Delegacin lvaro Obregn. (a) Modelo de Inversin realizado
con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa ZONDRES2D.

5.5. Coeneo, Michoacn.

El presente estudio fue realizado en un terreno que se encuentran en las

cercanas del Observatorio de Centelleo Interplanetario del Instituto de Geofsica
de la UNAM (MEXART). Lo anterior con el fin de evaluar la respuesta del terreno y
hacer una correlacin con los resultados arrojados por la Tomografa de
Resistividad Elctrica 2D. En el terreno de coordenadas UTM 217448 en X y
2192493 en Y.

5.5.1. Marco Geolgico

En la zona de estudio, as como en la mayor parte del estado de Michoacn,

convergen tres eventos que se distinguen claramente en tiempo y espacio: el
desarrollo de un arco insular del cretcico inferior desarrollado sobre una
secuencia de esquistos y pizarras pertenecientes al terreno Guerrero; un
vulcanismo terciario del eoceno-oligoceno que identifica la Sierra Madre
Occidental y el vulcanismo ms reciente del eje Neovolcnico. De acuerdo al
77
 Aplicacin en Datos de Campo Captulo 5

informe general complementario a la cartografa geolgica de minera y geoqumica

de Morelia (E14-1), fisiogrficamente la carta Morelia se ubica en la parte sur
central del Eje Neovolcnico y la zona norte de la Sierra Madre del Sur
observndose el traslape del primero sobre la segunda. En la superficie que cubre
la regin convergen tres unidades tectonoestratigrficas con caractersticas
propias, denominadas Terreno Guerrero, Sierra Madre Occidental y Eje
Neovolcnico. La porcin del terreno Guerrero que aflora se subdivide en:
Teteloapan y Arcelia-Huetamo. La relacin entres ambos terrenos es de tipo
tectnico, observndose como estructura principal una cabalgadura que se orienta
sensiblemente N-S.

5.5.2. Ubicacin geogrfica

El municipio de Coeneo se localiza al norte del Estado de Michoacn, en las

coordenadas 1948' de latitud norte y 10135' de longitud oeste, a una altura de
2,040 metros sobre el nivel del mar, cubre una superficie de 393.16 kms2 y
representa el 0.66 por ciento de la superficie del Estado. Limita al norte con los
municipios de Huaniqueo y Jimnez, al este con Morelia (figura 1), al sur con
Quiroga y Erongarcuaro y al oeste con Zacapu. Su distancia a la capital del
Estado es de 80 km.

Figura 5.4. Localizacin Geogrfica del municipio de Coeneo, Michoacn. Se muestra la ubicacin
geogrfica del Observatorio MEXART y la ubicacin del terreno donde se realizo el estudio.

78
 Aplicacin en Datos de Campo Captulo 5

5.5.3. Descripcin del levantamiento

En total se llevaron a cabo nueve perfiles 2D utilizando la modalidad Wenner-

Schlumberger. Para obtener las secciones de resistividad aparente se emplearon
28 electrodos con una separacin de 5 metros, cubriendo una longitud total de 135
metros y un total de 168 lecturas. En el presente trabajo se utilizaron 3 de las 9
secciones de resistividad aparente ms representativas.

Figura 5.5. Localizacin de los perfiles bidimensionales.

5.5.4. Resultados

Como se puede observar en las tres secciones, se diferencian tres unidades. La

unidad A, con los valores ms bajos de resistividad, se puede interpretar como el
suelo saturado con agua, el cual se compone de arcillas plsticas que sirven para
cultivo de maz. La unidad B se asocia a basaltos sin fracturar pero saturados con
agua por lo consiguiente representan valores de resistividad menores a los ms
comunes para este tipo de roca. Por ltimo se tiene la unidad C la cual se puede
interpretar como basaltos fracturados infiltrados con agua.

79
 Aplicacin en Datos de Campo Captulo 5

Figura 5.6. Modelos de inversin del perfil 1 realizado en Coeneo Michoacn. (a) Modelo de
Inversin realizado con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa
ZONDRES2D.

Figura 5.7. Modelos de inversin del perfil 2 realizado en Coeneo Michoacn. (a) Modelo de
Inversin realizado con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa
ZONDRES2D.

80
 Aplicacin en Datos de Campo Captulo 5

Figura 5.8. Modelos de inversin del perfil 3 Realizado en Coeneo Michoacn. (a) Modelo de
Inversin realizado con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa
ZONDRES2D.

5.6. Ixtln de los Hervores, Michoacn.

El sitio de estudio se localiza en un campo geotrmico el cual se encuentra al Este

del lago de Chapala, a una altitud de 1535 m.s.n.m Con un rea de alrededor de 2
km cubierta principalmente de manantiales de aguas termales a una temperatura
promedio de 90, lagunas de lodo, fumarolas y emanaciones de vapor. Se
perforaron 3 pozos de gran profundidad, uno de ellos se descontrol y arroja una
columna de vapor a una altura de 30 m aproximadamente, la cual dur 20 aos
emanando de manera natural. Actualmente tiene un uso recreativo simulando un
giser, el cual a veces se activa de manera natural mientras que otras tantas tiene
que ser auxiliado por un compresor de aire.

En un principio, los estudios en el rea de Ixtln de los Hervores por parte de CFE
tenan como objetivo un posible generador de energa por medio del campo
geotrmico de baja entalpa, sin embargo las temperaturas alcanzadas por los
pozos no fueron suficientes (250) para cumplir el objetivo.

5.6.1. Marco Geolgico.

El rea de estudio se localiza en la subprovincia fisiogrfica de Zonas de Fosas

Tectnicas, la cual se caracteriza por el alineamiento de sus cuencas, que para su
caso particular son fosas tectnicas rellenadas en direccin este-oeste, como es el
caso de los lagos de Chapala, Cuitzeo, Valle de Ixtln y de Zamora, entre otros.
As mismo, su constante actividad volcnica se refleja por innumerables aparatos
volcnicos y rocas extrusivas, conformadas esencialmente por rocas baslticas y
andesticas.

81
 Aplicacin en Datos de Campo Captulo 5

Las rocas que afloran en el rea estudiada son de origen extrusivo; parte de ellas,
principalmente las de menor cohesin y las formadas por minerales lbiles, han
sido fragmentadas, transportadas y depositadas en las depresiones estructurales
cercanas a sus fuentes, dando lugar a la existencia de depsitos clsticos de
diferentes espesores granulometra.

5.6.2. Ubicacin Geogrfica.

Se localiza al Noroeste del Estado de Michoacn, en las coordenadas 2010 de

latitud Norte y 10224 de longitud Oeste, a una altura de 1,530 metros sobre el
nivel del mar.
Limita al Norte con Tanhuato, al Este con Ecuandureo, al Sur con Chavinda y
Zamora, y al Oeste con Pajacuarn y Vista Hermosa. Su distancia a la capital del
Estado es de 174 kms.

Figura 5.9. Localizacin del rea de estudio.

5.6.3. Descripcin del levantamiento.

Se realizaron 7 perfiles de resistividad aparente con el arreglo Wenner-

Schlumberger, con un espaciamiento entre electrodos de 12 metros, cubriendo
una longitud de 324 m y 12 niveles alcanzando una profundidad aproximada de 60
metros. Tres de los tendidos (T1, T2 y T3) se encuentran dentro del balneario en
la parte posterior (Figura 5.10), en esta zona del terreno se localizan varios
rellenos de rocas volcnicas muy poco consolidadas, as mismo, el terreno est

82
 Aplicacin en Datos de Campo Captulo 5

segmentado por un par de riachuelos de 1.50 m de anchura y muy poca

profundidad, entre 15 y 30 cm, los tendidos se realizaron en estas condiciones,
atravesando las dos estructuras mencionadas. Cabe mencionar que el tercer
tendido no corresponde con la secuencia debida, ya que varias edificaciones y
estructuras geolgicas (giser) impedan el libre trnsito.

Los tendidos T4 y T5 (Figura 5.10) se realizaron en la parte frontal del balneario,

sobre la calle. En sta, se encontraban dos tipos de suelo, uno de los cuales
corresponde a un relleno de rocas volcnicas poco consolidadas y el otro a un tipo
de suelo ms consolidado y de mayor suavidad que se encontraba anegado con
poca profundidad, a un lado se hallaban rastros de huellas de desecacin, lo cual
caus problemas con el contacto electrodo-suelo.

En la Figura 5.10, los perfiles T6 y T7 se realizaron a un lado de la carretera Ixtln-

Zamora, en la cual haba un mismo suelo a lo largo de los dos tendidos, ste era
de origen lacustre y muy consolidado, se tuvieron problemas para establecer los
electrodos y con el contacto.

Para fines prcticos solo se procesaron los datos obtenidos en el TE1 por ser uno
de los ms representativos.

Figura 5.10. Localizacin de los sondeos.

83
 Aplicacin en Datos de Campo Captulo 5

5.6.4. Resultados

En los modelos de inversin de la Figura 5.11, se aprecian dos unidades; superior,

que aproximadamente vara entre 1 y 10 m, con resistividades de 3 a 10 ohm-m, y
la inferior, cuyo espesor no es posible definir, las variaciones de resistividades son
bajas, sin embargo con presencia de alto contraste, el cual se puede distinguir por
una zona de muy baja resistividad, la cual se asocia a las zonas geotrmicas que
se encuentran en el rea de estudio.

Figura 5.11. Modelos de inversin de uno de los perfiles realizados en Ixtln de los Hervores,
Michoacn. (a) Modelo de Inversin realizado con el programa RES2DINV. (b) Modelo de Inversin
realizado con el programa ZONDRES2D.

5.7. El Pahual, Espinal de Veracruz.

Este estudio fue realizado con el objetivo comn de la bsqueda de posibles

zonas o formaciones susceptibles de explotarse como abastecimiento de agua. Se
realiz un estudio cuya finalidad era localizar formaciones acuferas en el
subsuelo, tanto someras como profundas, que pudieran ser explotadas como
puntos de captacin de agua para abastecimiento.

5.7.1. Marco Geolgico

La zona de El Pahual est ubicada dentro de las formaciones Guayabal,

Chapopote-Tantoyuca, Palma Real, Meson, Escolin y Tuxpan. Todas contienen
litologas muy semejantes que consisten en areniscas, lutitas y limolitas, y por su
medio ambiente de depsito sus contactos son discordantes.
84
 Aplicacin en Datos de Campo Captulo 5

Figura 5.12. Carta Geolgica- Minera que contiene informacin de la localidad de El Pahual.

5.7.2. Ubicacin Geogrfica.

Se encuentra ubicado en la zona central del Estado, en las coordenadas 2015 de

latitud Norte y 9724 de longitud Oeste, a una altura de 100 metros sobre el nivel
del mar. Limita al Norte con Coatzintla; al Noreste con Papantla; al Oeste con
Coyutla; al Suroeste con Coxquihui; al Sur con Zozocolco y Estado de Puebla. Su
distancia aproximada al Noroeste de la capital del Estado es de 95 Km.

Figura 5.13. Ubicacin geogrfica de la localidad de El Pahual.

85
 Aplicacin en Datos de Campo Captulo 5

5.7.3. Descripcin del levantamiento.

Se realiz un perfil geoelctrico 2D mediante la tcnica denominada Tomografa

de Resistividad Elctrica. La adquisicin se llevo a cabo con un dispositivo Wenner
formado por 20 electrodos espaciados una distancia de 10 m, consiguiendo una
profundidad aproximada de 32 metros.

5.7.4. Resultados

Como se puede observar en la Figura 5.14, los materiales analizados en este perfil
presentan una tendencia a la disposicin horizontal y se encuentran agrupados en
tres capas principales cuyas caractersticas son, de techo a piso, las siguientes:

Una capa superficial, indicada con el nmero 1, de extensin generalizada a lo

largo de todo el perfil, atribuible a la unidad arcillosa suprayacente, con un espesor
que, de manera puntual, podra superar los 2 m.

Una segunda capa, indicada con el nmero 2, atribuible a la zona saturada

(acufero), compuesta de arenas no consolidadas, con un espesor cercano y algo
superior a 25 m. El posible nivel fretico del acufero se encuentra
aproximadamente a 5 o 6 metros de profundidad.

Una tercera capa, indicada con el nmero 3, la capa inferior, que correspondera a
las arcillas infrayacentes las cuales tiene su techo a una profundidad aproximada
de 32 metros.

Figura 5.14. Modelos de inversin del perfil realizado en Pahual. (a) Modelo de Inversin realizado
con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa ZONDRES2D.

86
 Aplicacin en Datos de Campo Captulo 5

Como se puede observar en la imagen elctrica (Figura 5.14), los lmites entre las
capas no son bien definidos, por lo que se realizo un Sondeo Elctrico Vertical
para poder tener una informacin ms certera de ello.

El dispositivo utilizado para la adquisicin de los datos fue el Schlumberger con

una abertura AB/2 mnima de 1 metro y mxima de 62.5 metros alcanzando una
profundidad aproximada de 30 metros como puede apreciarse en la tabla.

El procesamiento de los datos se realiz con el software IPI2WIN.

Figura 5.15. Perfil del Sondeo Elctrico Vertical realizado en la zona de Pahual.

Tabla 4. Valores para cada una de las capas del perfil del SEV.

N Resistividad Espesor en Profundidad en

en Ohm-m m. m.
1 9.015 1.301 1.301
2 3.677 1.522 2.823
3 6.846 26 28.82
4 3.955

Como se puede ver la primera capa tiene una resistividad de aproximadamente 9

Ohm-m y un espesor de 1.3 metros, la segunda capa que es donde se encuentra
el posible acufero esta a una profundidad de 2.8 metros, tiene un espesor
aproximado de 1.5 metros y una resistividad de 3.7 Ohm-m. La tercera capa tiene
una resistividad de 6.8 Ohm-m y un espesor 26 metros. El SEV reporta una capa
ms que la tomografa de una resistividad de alrededor de 4 Ohm-m.

87
 Aplicacin en Datos de Campo Captulo 5

5.8. Valle Etla, Oaxaca

En Mxico, alrededor de 67% del agua total suministrada para las ciudades
principales se obtiene del agua subterrnea. La sobrexplotacin, contaminacin
por la intrusin de aguas residuales salinas, domsticas e industriales est
amenazando los niveles de los acuferos. El acufero de Etla es un ejemplo de
este problema, porque existen una gran cantidad de pozos de extraccin y al
mismo tiempo los residuos domsticos que van a los drenes naturales. Este
acufero produce el 80% del agua total suministrada a la ciudad de Oaxaca y a la
irrigacin agrcola. El valle est situado en la porcin central del estado de
Oaxaca, cerca de la ciudad de Oaxaca, la Capital del Estado. El rea ha sufrido un
incremento en las necesidades de agua en los ltimos veinte aos. La extraccin
del agua se ha obtenido a travs de 190 pozos, aproximadamente, durante los
ltimos 19 aos. En el presente, las fuentes principales de contaminacin se
encuentran en las aguas residuales, y los fertilizantes domsticos e industriales.

5.8.1. Marco Geolgico

El estado de Oaxaca es atravesado por dos importantes sistemas de montaa: la

sierra Madre del Sur, la cual se extiende del sur de Jalisco al istmo de
Tehuantepec, y la sierra Madre Oriental. Ambos convergen hacia el este del
estado formando varios valles con acuferos potenciales. La ciudad de Oaxaca
est situada en la ensambladura de tres valles: el Etla al noroeste, el Tlacolula-
Mitla al sureste y el Zimatln al sur (Figura 5.16). El ro de Atoyac fluye en
direccin Sur, con los valles de Etla y Zimatln, y l es ensamblados por el Mitla
Ro en la ciudad de Oaxaca. El valle de Etla est limitado por la Sierra de Jurez
al este y la Sierra de Oaxaca al oeste (Figura 5.16). El valle de Etla esta
principalmente compuesto por depsitos lacutre-fluviales y en menor proporcin
por rocas volcnicas y depsitos piroclsticos de Edades Mioceno-Cuaternarios.
El sistema fretico del acufero del valle est compuesto por sedimentos de arenas
y gravas del Terciario-Cuaternarios, con un espesor que se extiende de 20 a 50 m
(Flores-Mrquez y otros, 2001). La base de este acufero es una secuencia
consolidada de arcilla, con un espesor variable alcanzando cerca de 36 M.
5.8.2. Ubicacin Geogrfica

El Valle de Etla est localizado a 1703 N y a 9644 W, tiene una altitud de 1700
metros sobre el nivel del mar (Figura 5.16)

88
 Aplicacin en Datos de Campo Captulo 5

Figura 5.16. Localizacin del Valle de Etla.

5.8.3. Descripcin del Levantamiento

Se realiz un perfil geoelctrico 2D mediante la tcnica denominada Tomografa

de Resistividad Elctrica. La adquisicin se llevo a cabo con un dispositivo
Wenner- Schlumberger formado por 21 electrodos espaciados una distancia de 10
m, cubriendo una longitud total de 200 metros, consiguiendo una profundidad
aproximada de 32 metros con 8 niveles y un total de 87 datos.

5.8.4. Resultados

En la Figura 5.17 se puede observar que los materiales analizados en este perfil
presentan una tendencia a la disposicin horizontal y se encuentran agrupados en
tres capas principales.

89
 Aplicacin en Datos de Campo Captulo 5

Una capa superficial, de extensin generalizada a lo largo de todo el perfil,

atribuible a areniscas, con un espesor que, de manera puntual, podra superar los
5 m.

Una segunda capa, atribuible a la zona saturada (acufero), compuesta de arenas

no consolidadas, con un espesor cercano y algo superior a 15 m. El posible nivel
fretico del acufero se encuentra aproximadamente a 10 metros de profundidad.

Una tercera capa, la capa inferior, que correspondera a areniscas infrayacentes

las cuales tienen su techo a una profundidad aproximada de 30 metros.

Figura 5.17. Modelos de inversin del perfil realizado en Valle Etla. (a) Modelo de Inversin
realizado con el programa RES2DINV. (b) Modelo de Inversin realizado con el programa
ZONDRES2D.

En la Figura 5.17 se puede observar que en la imagen elctrica no se alcanzan a

definir bien los horizontes de las capas por lo que se realizo un Sondeo Elctrico
Vertical para poder tener una informacin ms certera de ello, al igual que en el
ejemplo anterior (El Pahual) y los datos que se obtuvieron son los siguientes
(Figura 5.18).

90
 Aplicacin en Datos de Campo Captulo 5

Figura 5.18. Perfil del Sondeo Elctrico Vertical realizado en la zona del Valle de Etla.

Tabla 5. Valores para cada una de las capas del perfil del SEV.

N Resistividad Espesor en Profundidad en

en Ohm-m m. m.
1 77.84 0.5775 0.5775
2 98.48 2.282 2.859
3 8.65 23.38 26.24
4 87.07

Como se puede ver la primera capa tiene una resistividad de aproximadamente 78

Ohm-m y un espesor de 0.6 metros. La segunda capa tiene una resistividad de 6.8
Ohm-m y un espesor 26 metros. La tercera capa que es donde se encuentra el
posible acufero esta a una profundidad de 26 metros, tiene un espesor
aproximado de 23 metros y una resistividad de 8.65 Ohm-m. El SEV reporta una
capa ms que la tomografa de una resistividad de alrededor de 87 Ohm-m.

Como se puedo observar en los dos ltimos ejemplos, los programas tienen
ciertas limitantes al definir los lmites o espesores en un medio estratificado, por
ello se recomienda realizar un estudio con la tcnica de SEV para complementar
la informacin obtenida con la tcnica de Tomografa.

91
 Conclusiones y Recomendaciones Captulo 6

CAPTULO 6. CONCLUSIONES Y RECOMENDACIONES

El mtodo de Tomografa de Resistividad Elctrica resuelve de manera aceptable

estructuras geolgicas complejas con respecto a otras tcnicas de corriente continua.

Los problemas numricos involucrados para resolver un problema inverso caen dentro de
tres categoras: El clculo de la respuesta del modelo (problema directo), el clculo de las
sensitividades y la solucin de un problema de minimizacin a gran escala expresado
como un sistema de ecuaciones lineales.

Los resultados obtenidos de los modelos sintticos muestran la utilidad tcnica para inferir
la presencia de cuerpos en el subsuelo. Esto nos muestra las ventajas y desventajas de
cada modalidad para poder as emplear la ms adecuada al planear un estudio. Las
modalidades Dipolo-Dipolo y Wenner-Schlumberger son menos sensitivas al ruido
aleatorio a profundidad y lateralmente, la modalidad Wenner es mucho ms sensitiva al
ruido aleatorio.

La implementacin del programa Profiler V2 para realizar el modelado directo e inverso de

datos de Tomografa de Resistividad Elctrica 2D resulto ser una herramienta de apoyo
muy verstil.

El programa ZONRES2D til en el procesamiento de datos de Tomografa de Resistividad

Elctrica 2D, es un programa multiusuario, una vez que se tiene registrado presenta una
gran ventaja sobre otras herramientas de procesamiento, debido a que se puede utilizar
en un nmero indefinido de computadoras y trabaja sobre cualquier sistema operativo. No
requiere de instalacin especial en una computadora, se puede trabajar con l inclusive
desde una USB, lo cual lo hace muy porttil. Esto presenta una gran ventaja para
alumnos que requieran de una herramienta de procesamiento de datos de resistividad.

En comparacin con el programa RES2DINV que es un programa comercial probado,

reconocido, recomendado y utilizado en varios pases, el programa ZONDRES2D
resuelve bien las estructuras e incluso en algunas ocasiones obtiene imgenes elctricas
de mejor calidad, porque puede realizar la inversin con un nmero mayor de iteraciones
a diferencia del RES2DINV. Adems, el programa RES2DINV tiene un alto costo y es un
programa monousuario, es decir, trabaja en una sola computadora, necesita una llave
fsica para poder activar todas sus funciones, sin esta llave se puede ocupar solamente
como demo y esto limita a realizar una inversin de datos con tres iteraciones
nicamente.

Las imgenes elctricas 2D obtenidas en todas las zonas estudiadas arrojan buenos
resultados, coinciden con la geologa y ubicacin de las estructuras, aunque tiene ciertas
limitaciones al definir modelos de capas, como se pudo observar en el caso de El Pahual
y Valle Etla.

El algoritmo de inversin que define mejor las imgenes elctricas 2D es el de

suavizamiento condicionado.

92
 Conclusiones y Recomendaciones Captulo 6

RECOMENDACIONES

En la planificacin de un estudio con Tomografa de Resistividad Elctrica, se deben

tomar en cuenta los siguientes aspectos:

Se recomienda emplear un resistivimetro que sea lo suficientemente sensitivo y que los

electrodos tengan buen contacto con el terreno de estudio.

Se recomienda emplear la modalidad Wenner si el rea de estudio es un rea con ruido,

el estudio tiene un tiempo limitado para realizarse y se requiere una buena resolucin
vertical.

Si se requiere una buena resolucin horizontal, una cobertura de datos importante y existe
un buen contacto con el terreno, se recomienda usar la modalidad Dipolo-Dipolo.

En caso de no estar seguros de la estructura a encontrar y se necesita una buena

resolucin horizontal y vertical, se recomienda emplear la modalidad Wenner-
Schlumberger.

Si se tiene un sistema con un nmero limitado de electrodos, la modalidad polo-dipolo

podra ser una opcin viable.

La modalidad Polo-Polo se recomienda para estudios con espaciamientos pequeos entre

electrodos y una buena cobertura horizontal.

Debido a la limitacin del programa ZONDRES2D para definir un modelo de capas, se

recomienda realizar un estudio con la tcnica de Sondeo Elctrico Vertical (SEV) para
corroborar la informacin obtenida en el estudio de tomografa.

Al realizar la interpretacin de las imgenes elctricas se recomienda tener en cuenta la

profundidad media de estudio del arreglo empleado, debido a que el programa
ZONDRES2D hace una extensin de la profundidad y esto puede generar errores en la
interpretacin.

Se recomienda el uso de los programas Profiler V2 y ZONDRES2D como herramientas

para el modelado y procesado de datos adquiridos en prcticas de campo, o en las
asignaturas de Prospeccin Elctrica y Procesamiento de datos Geofsicos donde se
requieren este tipo de herramientas.

Los programas ProfileR y ZONDRES2D se pueden obtener en la siguiente direccin:

[Link]

93
 Bibliografa

BIBLIOGRAFA

lvarez, E.D., Arrubarena, M.M., et al., (2010). Reporte de Prctica Profesional,

Ixtln de los Hervores, Michoacn. Facultad de Ingeniera. UNAM. Mxico, D.F.

Baker,R.D., (1992). A simple algorithm for electrical imaging of the subsurface.

First Break, v. 10 - 2, p. 53-62

Constable, S.C., Parker, R.L. and Constable, C.G., (1987). Occams inversion: A
practical algorithm for generating smooth models from electromagnetic sounding
data. Geophysics, 52, 289-300.

Caballero, J.G., Del Bosque, E.L., et al., (2009). Reporte Final de Prcticas
Profesionales (Coeneo, Michoacn). Facultad de Ingeniera. UNAM. Mxico, D.F.

De Groot-Hedlin, C. and Constable, S., (1990). Occam's inversion to generate

smooth, twodimensional models form magnetotelluric data. Geophysics, 55, 1613-
624.

Dey, A. and Morrison, H.F. (1979). Resistivity modeling for arbitrary shaped two-
dimensional structures. Geophysical Prospecting 27, 1020-1036.

Edwards, L. S., (1977). A modified pseudosection for resistivity and induced

polarization. Geophysics, 42, 1020-1036.

Ellis, R.G. and Oldenburg, D.W., (1994). Applied geophysical inversion.

Geophysical Journal International, 116, 5-11.

Flores, E.L., Chvez, R.E., Martnez, R.G., et al. (2001). Geophysical

Characterization of the Etla Valley aquifer, Oaxaca, Mxico. Instituto de Geofsica,
UNAM. Mxico, D.F.

Keller G.V., Frischknecht F.C., (1966). Electrical Method in Geophysical

Prospecting. Pergamon Press Inc. Oxford

Koefoed O., (1979). Geosounding Principles 1: Resistivity sounding

measurements. Elsevier Science Publishing Company, Amsterdam.

Kufman, A.A., (1992). Geophysical field theory and method, part A. Academic
Press.

Lines L.R. and Treitel S. (1984). Tutorial: A review of least-squares inversion and
its application to geophysical problems. Geophysical Prospecting, 32, 159-186.

Li Y. and Oldenburg D.W. (1992). Approximate inverse mappings in DC resistivity

problems. Geophysical Journal International 109, 343-362.

94
 Bibliografa

Loke, M.H., (1994). The inversion of two-dimensional resistivity data. Unpubl.

PhD thesis, Un. Of Birmingham.

Loke M.H., Baker R.D., (1995). Least-squares deconvolution of apparent resistivity

pseudosection. Geophysics, v. 60, p. 1682-1690

Loke M.H., Baker R.D., (1996). Rapid least-squares inversion of apparent

resistivity pseudosections by quasi- Newton method. Geophysical Prospecting, v.
44, p. 131- 52

Locke, M. H., 1995. RES2DMOD version 2.20: A 2D resistivity forward modeling.

Lopez H., Loke H., Fanton O., Rub E. (2004). Tcnicas prcticas para
investigacin de resistividad en dos y tres dimensiones (tomografa elctrica 2D y
3D).

Lpez Ramos E., (1985). Geologa de Mxico. CONCACYT. Tomo II, tercera
edicin, Mxico, D.F.

Olayinka, A.I. and Yaramanci, U., (2000). Use of block inversion in the 2-D
interpretation of apparent resistivity data and its comparison with smooth
inversion. Journal of Applied Geophysics, 45, 63-82.

Oldenburg, D. W. and Li, Y., (1994). Inversion of induced polarization data.

Geophysics, 59, 13271341.

Orellana E., (1972). Prospeccin Geoelctrica en corriente contina. Paraninfo,

segunda edicin, Madrid, Espaa, p. 24 72

Sasaki, Y., (1992). Resolution of resistivity tomography inferred from numerical

simulation. Geophysical Prospecting, 40, 453-464.

Servicio Geolgico Mexicano (1996), Carta Geolgico Minera de Poza Rica,

Veracruz, F14-12, Esc. 1:250000.

Sumanovac F., Dominkovic A.S., (2007). Determination of Resolution limits of

Electrical Tomography on the Block Model in a homogenous environment by
means of electrical modeling. Rudarsko Geolosko Naftni Zbornik, v. 19, p. 47-56

Telford, W.M., Geldart, L.P. and Sheriff, R.E., (1990). Applied Geophysics
(second edition). Cambridge University Press.

Tong, L. and Yang, C., (1990), Incorporation of topography into two-dimensional

resistivity inversion. Geophysics, 55, 354-361.

Tutorial RES2DINV ver. 3.55. Geotomo Software, enero 2006.

95
 Bibliografa

Urbieta Gmez J., (2001). La Tomografa elctrica en los mtodos de corriente

contina. Tesis de Licenciatura, Facultad de Ingeniera, Universidad Nacional
Autnoma de Mxico.

Urbieta Gmez J., (2009). Evaluacin de las tcnicas de adquisicin de datos 2D

y 3D con tomografa de resistividad elctrica en pozo. Tesis de Maestra.
Posgrado. UNAM.

Urbieta, G.J. (2008). Estimacin De La Distribucin Espacial De Resistividades

Del Subsuelo Para Localizar Un Posible Acufero Por Tomografa Elctrica En El
Pahual, Municipio De Espinal, Veracruz. Departamento de Geofsica. Facultad
Ingeniera, UNAM. Mxico, D.F.

[Link]
=1G1GGLQ_ESMX365&sa=2

[Link]

[Link]
%A9ctrica%20como%20Herramienta%20de%20Diagnostico%20en%20la%20Aus
cultaci%C3%B3n%20de%20Presas%20de%20Materiales%[Link]

[Link]

96
 Anexos

ANEXO A

Algunos modelos sintticos generados con el software RES2DINV.

En las Figuras que se muestran en este anexo se muestran las imgenes

elctricas generadas con el software RES2DIV de los modelos sintticos
presentados en el Captulo 4 del presente trabajo con la finalidad de comparar el
poder resolutivo del programa ZONDRES2D. Se utilizaron 3 de las 5 secciones de
resistividad aparente correspondientes a las diferentes modalidades por ser las
ms representativas

Como se puede observar en las Figuras (A1-A5) los cuerpos son determinados en
los diferentes modelos propuestos, en cualquier modalidad, al igual que lo hace el
programa ZONDRES2D (Capitulo 4), aunque cada programa defina caractersticas
diferentes (por muy pequeas que sean).

Figura A1. Secciones de resistividad obtenidas de un modelo sinttico de un bloque rectangular

para las modalidades a) Wenner, b) Wenner-Schlumberger y c) Dipolo-Dipolo.

97
 Anexos

Figura A2. Secciones de resistividad obtenidas de un modelo sinttico de dos bloques horizontales
para las modalidades a) Wenner, b) Wenner-Schlumberger y c) Dipolo-Dipolo.

Figura A3. Secciones de resistividad obtenidas de un modelo sinttico de cuatro bloques

horizontales para las modalidades a) Wenner, b) Wenner-Schlumberger y c) Dipolo-Dipolo.

98
 Anexos

Figura A4. Secciones de resistividad obtenidas de un modelo sinttico de dos capas y un dique
para las modalidades a) Wenner, b) Wenner-Schlumberger y c) Dipolo-Dipolo.

Figura A5. Secciones de resistividad obtenidas de un modelo sinttico de una falla y un bloque
rectangular para las modalidades a) Wenner, b) Wenner-Schlumberger y c) Dipolo-Dipolo.

99
 Anexos

ANEXO B

Inversin de los datos de campo generados con otros algoritmos de

inversin en el software ZONDRES2D.

A continuacin se muestran las imgenes elctricas generadas con otros tipos de

algoritmos de inversin para algunos de los datos de campo presentados en el
Captulo 5 del presente trabajo. Para fines prcticos se utilizaron los datos de
campo de la Calle Segunda Cerrada de Minas, lvaro Obregn, los datos de
Coeneo, Michoacn y los datos de El Pahual, Espinal de Veracruz, por
considerarse los ms representativos.

Se puede observar que todos los algoritmos de inversin definen las estructuras
en cada seccin, separando los cuerpos por su alta o baja resistividad, sin
embargo no todos determinan bien la forma o la resistividad de los mismos, por
ello, se ocup el algoritmo de Suavizamiento Condicionado para procesar los
datos en el captulo 5, por considerarse el que mejor define las estructuras.

Figura B1. Modelo de inversin generado con el algoritmo de Marquardt del perfil realizado en la
calle Segunda Cerrada de Minas en la Delegacin lvaro Obregn.

Figura B2. Modelo de inversin generado con el algoritmo de Occam del perfil realizado en la calle
Segunda Cerrada de Minas en la Delegacin lvaro Obregn.

100
 Anexos

Figura B3. Modelo de inversin generado con el algoritmo de Occam- Marquardt del perfil realizado
en la calle Segunda Cerrada de Minas en la Delegacin lvaro Obregn.

Figura B4. Modelo de inversin generado con el algoritmo de Marquardt del perfil 1 realizado en
Coeneo Michoacn.

Figura B5. Modelo de inversin generado con el algoritmo de Occam del perfil 1 realizado en
Coeneo Michoacn.

Figura B6. Modelo de inversin generado con el algoritmo de Occam-Marquardt del perfil 1
realizado en Coeneo Michoacn.

101
 Anexos

Figura B7. Modelo de inversin generado con el algoritmo de Marquardt del perfil realizado en El
Pahual, Espinal de Veracruz.

Figura B8. Modelo de inversin generado con el algoritmo de Occam del perfil realizado en El
Pahual, Espinal de Veracruz.

Figura B9. Modelo de inversin generado con el algoritmo de Occam-Marquardt del perfil realizado
en El Pahual, Espinal de Veracruz.

102