Libro-Geofisica-Aplicada-Prot PDF
Libro-Geofisica-Aplicada-Prot PDF
Libro-Geofisica-Aplicada-Prot PDF
MÉTODO GEOELÉCTRICO
Antonio Carrasco Cortez
Ingeniero geólogo
Oruro - Bolivia
2020
Geofísica Aplicada Antonio Carrasco Cortez
Tengamos el valor de decir la verdad el cadalso la decapitación, puede que esto purifique la vida
que olvidamos vivir por seguir estrellas fugaces que nunca podemos alcanzar.
2
Geofísica Aplicada Antonio Carrasco Cortez
PROLOGO
Pretendo que este texto sirva de manual práctico a los métodos de la prospección
geofísica, se encuentra escrito principalmente para ingenieros geólogos, agrónomos e Ing.
civiles espero sin embargo que otras ramas de la ingeniería puedan considerarlo útil.
Muchas partes del libro tienen como base casos prácticos tratados breve y sucintamente,
que se trata deliberadamente en forma esquemática, aunque podría continuarse hasta
alcanzar un grado de detalle relativamente elevado y complejo.
He de mencionar finalmente para no cansarlo más, que el libro hoy escrito pueda inspirar
a otros a comprender y aplicar uno de los múltiples lenguajes de la geología aplicada.
3
Geofísica Aplicada Antonio Carrasco Cortez
INDICE GENERAL
Prologo ..................................................................................................................... 3
Capitulo VII Métodos de interpretación del sondeo eléctrico asistidos por ordenador .... 93
Capitulo XII Sección vertical de resistividad aparente y sección eléctrica .................. 157
4
Geofísica Aplicada Antonio Carrasco Cortez
Capítulo I
El hombre no existe más que por sus obras; aquel que nada ha hecho, pasa por la vida
como un fantasma.
No se puede empezar a hablar de los métodos geoelectricos sin antes realizar un análisis
básico de las corrientes eléctricas que circulan por el cuerpo solido de la tierra, la atmosfera
y los océanos que la envuelven, todas ellas constituyen un complejo sistema de cargas
móviles relacionadas entre sí.
5
Geofísica Aplicada Antonio Carrasco Cortez
La tierra solida parece tener un exceso de cargas negativas que provocan una influencia
constante de iones con carga positiva del aire, aunque la magnitud de estas cargas negativas
fluctúa con el tiempo y varia de un lugar a otro se mantiene de un modo continuo y de una
manera aún desconocida.
El grado de ionización del aire varía notablemente con la altura a la cual se realice la
medida, hora del día y la latitud, cerca de la superficie de la tierra la ionización es reducida.
Todo campo magnético es generado por un campo eléctrico o viceversa, ambos son
componentes de los campos electromagnéticos, el párrafo siguiente abordara brevemente el
concepto de campo eléctrico de la tierra, enfocando desde el campo magnético terrestre.
Se puede decir que dos terceras partes del campo magnético terrestre, es de origen externo
(es decir se deben a la acción de los astros y el sol), el tercio restante se debe a las
corrientes que circulan por el globo terrestre, si intentamos calcular la distribución de la
corriente en el interior de la tierra, que podría producir la componente interna del campo
magnético de seguro uno se tropezaría con la dificultad de hallar un modelo de tierra que
satisfaga todas las observaciones realizadas.
Las corrientes en la superficie terrestre giran alrededor de la tierra siguiendo al sol, se cree
que esto es una parte del sistema de corriente distribuida, (la corteza superficial al
interactuar con la radiación solar y por el efecto fotoeléctrico genera diferenciales de
potencial eléctrico en la misma) cuando menos por varios cientos de kilómetros de
profundidad e inducidas por el sistema de corriente diurno atmosférico, no se ha podido
demostrar de manera evidente la existencia en la superficie de la tierra de corrientes
verticales, solo se tiene certeza de la existencia de las desviaciones de corrientes
horizontales provocadas por variaciones horizontales en la litología de los materiales de la
corteza estos se reflejan físicamente en variaciones de conductividad en las rocas.
6
Geofísica Aplicada Antonio Carrasco Cortez
TENCIONES LOCALES
Las tenciones locales se deben a diversas causas, las reacciones químicas entre los niveles
de sub saturación y saturación de agua conjuntamente con los minerales formadores de
rocas y sedimentos, dan lugar a modelos característicos alrededor de ciertos cuerpos
mineralizados.
7
Geofísica Aplicada Antonio Carrasco Cortez
El progreso de la física a fines del siglo XIX y principios del XX, dio un impulso al
desarrollo de la geofísica aplicada en general.
Fig. 2 En 1912 Conrad Schlumberger, utilizando un equipo muy básico, elaboro el primer
mapa de curvas de potencial en su finca cerca de Caen, Normandía.
En 1912, utilizando un equipo muy básico, elaboro el primer mapa de curvas de potencial
en su finca cerca de Caen, en Normandía, antes de ampliar sus encuestas para las minas de
8
Geofísica Aplicada Antonio Carrasco Cortez
En los años 1919 y 1922 los científicos Suecos N. Lundberg y K. Sundberg, iniciaron con
sus trabajos la exploración eléctrica con campos variables y en particular los métodos
vasados en la observación de las líneas equipotenciales del campo eléctrico y la intensidad
del campo magnético.
9
Geofísica Aplicada Antonio Carrasco Cortez
En el año 1924 el geofísico Ruso A.A. Petrovskii realiza los primeros trabajos de
prospección eléctrica, mediante el método del campo eléctrico natural dándole al mismo
una metodología de interpretación cuantitativa hasta entonces desconocida.
En los años 1928 y 1929 la prospección eléctrica empezó, por primera vez con éxito a
utilizar en la exploración e investigación de estructuras petrolíferas y gasíferas (como
muestra los trabajos de M. N. Berdichepskii, en el compendio de geofísica aplicada de la
U.S.S.R.), en los años siguientes el volumen de estos trabajos crese de modo considerable,
bueno se podría decir que esta etapa de la prospección geoelectrica solo sucedió en la parte
oriental del mundo, ya que en la parte oxidental, el método eléctrico todavía permanecía y
permanece relegado en el proceso de exploración de hidrocarburos.
El francés Raimond Maillet, introduce a finales de los años treinta, la teoría de Darr
Zarrouk, demostrando así con todo rigor qua cada curva de sondeo eléctrico, no refleja, en
general el efecto independiente de espesores y resistividades, sino lo que influye
decisivamente en la morfología de dichas curvas, son el producto y cociente de espesores y
resistividades de cada capa.
En 1946 el ruso Zeleznak demuestra la gran utilidad del método de polarización inducida
en la localización de minerales metálicos y desde entonces el método se usa ampliamente
en la industria minera.
En 1953 el Francés L. Cagniart publica por primera vez un artículo sobre el método
magnetotelúrico.
10
Geofísica Aplicada Antonio Carrasco Cortez
Bueno a partir de los cincuenta desde mi punto de vista no se dieron más descubrimientos
en el campo de la geoelectrica, sin embargo, se perfeccionaron y mejoraron los principios
de los métodos descubiertos hasta entonces.
11
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo II
(Jaime Sáenz)
La geofísica aplicada por el método eléctrico tiene amplia aplicación para el estudio de las
estructuras geológicas, los yacimientos y criaderos minerales de la parte superior de la
corteza terrestre.
12
Geofísica Aplicada Antonio Carrasco Cortez
Sondeo electromagnético.
Calicata electromagnética.
Georradar (GPR.)
-Método de las ondas radiofónicas.
13
Geofísica Aplicada Antonio Carrasco Cortez
Entre todos estos métodos los más empleados y difundidos son los de resistividad y la
polarización inducida.
- Las reservas indicadas son reconocidas por labores mineras, sondeos mecánicos muy
separados o métodos geofísicos. (Yacimientos minerales de E. W. Petrascheck).
- Con el uso de sistemas geofísicos es posible determinar muchas de las propiedades físicas
de las rocas del subsuelo, estos datos si se interpretan con propiedad pueden arrojar mucha
luz sobre la litología de las rocas, y si se dispone de suficientes datos, de la forma de los
cuerpos de la roca. (Geología estructural de N. Billings).
14
Geofísica Aplicada Antonio Carrasco Cortez
V
Resistencia = = ΔV/I=R (ohm)
A
A = Intensidad de corriente.
15
Geofísica Aplicada Antonio Carrasco Cortez
Estos conceptos ayudaran a discernir mejor los encontrados en el acápite resistividad de las
rocas además de ser de importancia en el caso de realizarse una prospección eléctrica en
una zona geotérmica.
cuanto más grande es la resistividad tanto mayor es el campo eléctrico que se necesita para
generar atreves del material una densidad de corriente determinada las unidades de la
resistividad son: (v/m) / (A/(m*m)) = (V*m) /A donde el diferencial de potencial (Δv)
dividido entre la intensidad de corriente (A) según la ley de Ohm es igual a la resistencia
cuya unidad de medida es el ohmio (ohm) por tanto las unidades de medida de la
resistividad también pueden ser expresadas en “ohm*m”.
16
Geofísica Aplicada Antonio Carrasco Cortez
Donde la densidad de corriente (J) es igual a la intensidad de corriente dividida entre el área
de la sección (S) por tanto:
J=I/S
Y el campo eléctrico es igual al diferencial de potencial (ΔV) dividido entre la longitud (L):
E=ΔV/L
De donde “ΔV” puede ser medido fácilmente atreves del voltímetro simbolizado en el
grafico 3, con la letra “V” y la intensidad de corriente respectivamente por el amperímetro
simbolizado por la letra “A”, donde “S” es el área en metros cuadrados y “L” es la longitud
del cuerpo cilíndrico en metros.
17
Geofísica Aplicada Antonio Carrasco Cortez
Para entender de mejor manera el concepto de resistividad esta será explicada desde el
concepto de conductividad.
Metálica
Electrolítica
CONDUCTIVIDAD METALICA
Esta se caracteriza por el transporte de la corriente por la misma matriz de la roca es decir
que la corriente pasa por la masa rocosa este tipo de conductividad, se da en algunas menas
metálicas y esquistos grafiticos.
CONDUCTIVIDAD ELECTROLITICA
Este tipo de conductividad que presentan la mayoría de las rocas el primer gráfico de la fig.
6, muestra una barra de cobre por el cual atraviesa la corriente, de un modo análogo es
como actúa la conductividad metálica.
El segundo gráfico de la fig. 6 muestra el experimento del paso de la corriente atreves del
agua con sal (en si el agua sin sales disueltas es no conductora de corriente eléctrica) de un
modo análogo se desarrolla el fenómeno de conductividad electrolítica.
18
Geofísica Aplicada Antonio Carrasco Cortez
La resistividad de las rocas que tienen conductividad electrolítica tiene como función
decreciente a los siguientes factores:
- La porosidad total comunicante (se refiere a que, si los poros de roca se encuentran
comunicados entre sí, la conductividad será buena y si los poros no están comunicados
sucederá lo contrario).
- Resistividad del agua de imvicion (o resistividad del agua que se encuentra en los poros
de la roca).
- La resistividad del agua de imvicion disminuye cuando aumenta la cantidad de sales que
lleva disuelta, también depende de la naturaleza y composición de estas sales.
La porosidad (P) de una roca o suelo es la expresión del porcentaje de espacios vacíos
respecto al volumen total.
P = (100%*vp) / V
Donde (vp) es el volumen requerido de agua para llenar los poros y (V) el volumen total.
En la zona saturada de la formación sedimentaria el agua llena todos los intersticios, por lo
que la porosidad es una medida directa del agua contenida en la unidad de volumen.
Sin embargo, si realizamos un poso y lo bombeamos no podremos extraer toda esa agua en
razón de fuerzas moleculares y de tención superficial.
Er = (100*vr) /V
Ef = (100*vf) /V
Y como: V = vr + vf
Resulta: P = Er + Ef
19
Geofísica Aplicada Antonio Carrasco Cortez
Porosidad eficaz (Pz). - Es el volumen de poros conectados entre sí por unidad de volumen
total estos poros son los encargados de la conductividad electrolítica.
Pz = (100%*vpc) / V
Porosidad total = Er + Pz
Er: Porosidad de retención específica.
La porosidad eficaz varía entre 0,01 (1 %) Y 0,4 (40%).
Ejemplo: Una formación litológica compuesto por gravas y arenas que tiene 100 Km2 de
superficie (Ap) y un espesor saturado (h) de 50 m puede almacenar el siguiente volumen de
agua.
Volumen de agua Vol=superficie del cuerpo (Ap)* espesor del cuerpo (h)*Pz
20
Geofísica Aplicada Antonio Carrasco Cortez
21
Geofísica Aplicada Antonio Carrasco Cortez
Se nota que el mismo corresponde a dos campos de existencia como muestra la tabla de la
fig.8, entre los porcentajes probabilísticos se escoge al mayor y se procede a correlacionar
los porcentajes escogidos por medio de las filas con el correspondiente nombre de mineral
o roca.
22
Geofísica Aplicada Antonio Carrasco Cortez
TABLAS DE RESISTIVIDAD
23
Geofísica Aplicada Antonio Carrasco Cortez
24
Geofísica Aplicada Antonio Carrasco Cortez
Bueno para entender mejor el problema lo abordaremos desde el siguiente punto de vista, si
sobre un medio homogéneo eléctricamente se introduce una corriente continua de
intensidad “I” atreves del electrodo puntual “A” como muestra la figura 12, (b)
considerando el electrodo “B” en el infinito (en la práctica muy lejos del electrodo A).
Grafico (a)
Grafico (b)
26
Geofísica Aplicada Antonio Carrasco Cortez
Entonces esta fluye uniformemente en forma radial en todas las direcciones las líneas de
corriente en este caso coincidirán con los radios que parten del centro de los electrodos que
inyectan corriente al subsuelo y las superficies equipotenciales que son perpendiculares a
las líneas de corriente forman un sistema de semiesferas concéntricas en “A” (ver fig.13).
El signo negativo se debe a que la corriente fluye de los puntos de potencial elevados a los
puntos de menor potencial, la magnitud “2πr 2 ” es el área del casquete esférico, la
expresión nos muestra como al aumentar “r” el potencial disminuye, integrando la
expresión anterior se tiene que:
r
V =∫0 Ra Idr/(2πr 2 )
r
V= ( Ra I /2π )∫0 dr/(r 2 ) = Ra*I /2πr………..… (6)
27
Geofísica Aplicada Antonio Carrasco Cortez
Si el electrodo “B” (-) que cierra el circuito se encuentra a una distancia finita de “A” (+)
figura 12 (a) el potencial en el punto N será:
Y en el punto M será:
Ra= (2*π/ ((1/AM) -(1/BM) -(1/AN) + (1/BN))) *(ΔV /I) =Ra (ohm*m)
Dónde:
28
Geofísica Aplicada Antonio Carrasco Cortez
Salvo que el libro está destinado a geólogos, le comunicarle que las cosas no son tan fáciles
lo mejor será seguir leyendo el texto, este se encuentra especialmente endulzado para el
apetito siempre antipático hacia las ciencias exactas que tiene el buen geólogo a
continuación, algunos problemas que evaluaran el capítulo:
1.-Un alambre metálico de resistencia R, se divide en tres partes iguales y se conecta una
pegada a la otra para formar un alambre nuevo de longitud igual a una tercera parte de la
longitud original cual es la resistencia de este nuevo alambre
Respuesta: R/9.
3.- Calcule la resistividad aparente (Ra) del suelo donde se pretende instalar una planta
telefónica tomando en cuenta los siguientes datos de campo:
A M es = a en metros
M N es = a en metros
N A es = a en metros
Fig-14
Respuesta: Ra = (ΔV/I) *2 π a.
29
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo III
Cuando pones la proa visionaria hacia una estrella y llevas el resorte misterioso de un
ideal, solo vives por esa partícula de ensueño que se sobrepone a lo real, cuando se te
anuda la garganta al recordar la cicuta impuesta a Sócrates, la cruz izada para Jesucristo,
inclinándote con igual veneración ante los creadores de verdad o de belleza.
(Jose Ingenieros)
Fig.15 a Fig. 15 b
30
Geofísica Aplicada Antonio Carrasco Cortez
Los esquemas de las figuras 15 y 16 pertenecen a los primeros circuitos que muestran la
arquitectura básica de los equipos de prospección eléctrica.
TRANSMISOR
La labor del transmisor es de inyectar corriente continua o directa al sub suelo, claro que la
corriente deberá de ser la suficiente para que en otro punto de la superficie del terreno se
tenga la recepción de la corriente inyectada en forma de diferencial de potencial.
31
Geofísica Aplicada Antonio Carrasco Cortez
fig. 18 Circuito transmisor con elevador de voltaje alterno y generador de baja frecuencia
(fuente de diseño propia).
32
Geofísica Aplicada Antonio Carrasco Cortez
- Generador de corriente:
- Elevador de voltaje:
- Conmutador selector
Es simplemente una llave tipo interruptor que selecciona los voltajes otorgados por el
inversor para que un único voltaje seleccionado ingrese a la etapa de rectificación de
corriente, por lo general este conmutador debe de soportar intensidades de corriente que
fluctúan entre los: 0.1 o 5 amperios.
- Rectificador de corriente
- Filtrado de onda
Esta etapa se encuentra compuesta por un banco de filtros capacitores electrolíticos según
los requerimientos del circuito, el capacitor es un dispositivo eléctrico que acumula energía
en forma de campo eléctrico, cuando la tención varia a lo largo de un ciclo eléctrico, la
energía puede ser acumulada durante algún tiempo y luego ser devuelta el mismo, cumple
la función de convertir la onda eléctrica cuadrada proveniente del banco de diodos en la
33
Geofísica Aplicada Antonio Carrasco Cortez
etapa de rectificación de corriente, en una función eléctrica recta continua, que es el tipo de
corriente requerida para la prospección eléctrica.
- Reóstato
Es una resistencia variable capaz de soportar la intensidad de corriente que circula, por el
suelo mediante los electrodos de corriente, este dispositivo denominado reóstato se
encarga de regular la intensidad de corriente de salida del transmisor aumentando o
disminuyendo la resistencia eléctrica entre los electrodos de corriente (es muy común en
la geología del altiplano boliviano, que en el terreno la primera capa tenga una alta
conductividad debido al alto contenido de arcillas del suelo, con aperturas cortas de
electrodos de corriente, la intensidad de corriente en los mismos suele sobre pasar la
capacidad de flujo de corriente de los diodos rectificadores, para evitar esto se debe de
conectar una resistencia variable en serie antes de la conexión del amperímetro que registra
la intensidad de corriente suministrada a los electrodos de corriente A y B, también
conectada en serie para poder regular la intensidad de corriente suministrada al suelo).
-El amperímetro se conectan tal y como muestra el esquema de la fig. 19 este registra la
intensidad de salida del transmisor.
Este cumple la función de variar la frecuencia de la onda cuadrada entregada por el banco
de capacitores, para así evitar la interferencia con el ruido eléctrico o SP, sus principales
características son:
RESEPTOR
Este voltaje opuesto puede ser regulado mediante los reóstatos del esquema de la fig. 19,
hasta que la batería de bajo voltaje solo genere la corriente equivalente al ruido una vez
logrado este proceso, solo se registrara el voltaje correspondiente a los inyectados por el
transmisor, claro que el proceso de compensación de ruido mediante los reóstatos se la debe
de realizar mientras el transmisor este apagado y para cada posición de M N en el terreno.
LOS OHMIMETROS
Se deben de disponer de un par de ohmiometros, uno para ser conectados a las salidas del
transmisor mientras el mismo se encuentre apagado, para medir la resistencia del suelo y
así verificar que los electrodos de transmisión estén bien conectados a tierra.
Un segundo ohmímetro deberá de estar conectado a las entradas del receptor también para
verificar que los electrodos de recepción o de potencial estén bien conectados al suelo, el
35
Geofísica Aplicada Antonio Carrasco Cortez
uso de los mismos es de suma importancia puesto que controlaran la calidad de la conexión
del equipo al suelo como muestra el grafico 21.
Una vez conectados los ohmimetros, en paralelo a la línea AB o MN, los valores
registrados en el mismo no deberán exceder en rango los 100 M ohm. Para aperturas cortas
de AB y MN, es decir que la resistencia del suelo nunca deberá de tender al infinito.
36
Geofísica Aplicada Antonio Carrasco Cortez
c) ELECTRODOS
Electrodos metálicos. - Pueden ser utilizados tanto para conectar el transmisor a tierra como
para conectar el receptor a tierra están conformados por barras de acero inoxidable o por
barras de bronce o cobre endurado, no es recomendable usar este tipo de electrodos como
receptores porque en los mismos puede llegar a producirse un diferencial de potencial
variable debido a la acción electrolítica lo cual distorsiona o enmascara el potencial real
medido.
Electrodos impolarizables.- Estos solo pueden ser usados como electrodos de recepción,los
electrodos impolarizables más sencillos consisten en una varilla de cobre sumergidos en
una solución saturada de sulfato de cobre colocada en una vasija porosa, el contacto se
realiza enterrando la parte porosa de la vasija, como los electrodos tienen la misma
concentración no se produce ningún potencial adicional en estos y las medidas tomadas
solo corresponden al diferencial de potencial producido durante la lectura de ruido o la
lectura de voltaje tomado por el receptor.
37
Geofísica Aplicada Antonio Carrasco Cortez
El cable que se emplea en la prospección eléctrica, no tiene nada de especial solo requiere
que sea liviano y que se a resistente además de tener una buena envoltura aislante se emplea
generalmente cable multifibras de cobre aislado.
Respecto a los carretes estos son fáciles de construir o se pueden usar cualquier carrete de
cable que se tenga a disposición al igual que los cables estos deben de ser resistentes y
livianos para poder ser movilizados en el campo.
e) IMPLEMENTOS
- Agua con sal común al 10% la solución debe de ser vertida en el contacto entre
los electrodos y el suelo para mejorar la conducción eléctrica (recuerde el siguiente
consejo no importa el tamaño de los electrodos sino la calidad de la conexión).
38
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 26 Circuito transmisor y receptor con corrector de ruido (fuente de diseño propia).
La desventaja de estos instrumentos yace en que si los cables son largos los efectos de
inducción electromagnética pueden ser molestos e incluso podrían malograr el proceso de
toma de datos.
39
Geofísica Aplicada Antonio Carrasco Cortez
A continuación, los siguientes problemas que reforzaran los conceptos del capitulo
1.- Cinco pilas iguales que despiden un voltaje de 2 voltios individualmente y 0.6 ohmios
de resistencia interna, hallar la intensidad de corriente que suministran a una resistencia
exterior de carga de 17 ohmios cuando las pilas se conectan en a) serie b) paralelo
b) 0.117 amperios
2.-Un amperímetro tiene una resistencia interna de 0.006 ohmios y cada división de su
escala de medida equivale a un amperio con que resistencia R conectada en paralelo se
debe shuntar para que cada división de su escala corresponda a una corriente de 5 amperios
de intensidad.
3.-Un condensador cuya capacidad es iguala 250 uF. Se conecta a la red de corriente alterna
determinar su resistencia siendo la frecuencia de 50 200 y 400 Hz.
40
Geofísica Aplicada Antonio Carrasco Cortez
Capítulo IV
CONFIGURACION de ELECTRODOS
EN EL METODO ELÉCTRICO
La ciencia tiene más sinceridad cuando habla de su impotencia que de sus triunfos.
(Primo Arrieta)
Se dice que los arreglos son lineales, por el hecho que los electrodos son colocados a lo
largo de una línea imaginaria trazada en la zona de estudio y simétricos a razón de que a
partir del punto de sondeo se divide la recta imaginaria en dos partes iguales donde los
electrodos deberán de ser colocados en ambas partes de la recta en las mismas posiciones
como en una imagen especular.
41
Geofísica Aplicada Antonio Carrasco Cortez
DISPOSITIVO SCHLUMBERGER
En el cual los electrodos de potencial ocupan la parte central y la distancia entre los
electrodos de potencial o recepción MN es mucho menor a los de inyección de corriente
AB como muestra la fig. 28, teóricamente la distancia MN debe de ser 1/5 de la distancia
AB, pero en la práctica se obtienen muy buenos resultados con MN igual o mayor a 1/3 de
la distancia AB.
K (m) = π ((AM *AN) /MN); el factor geométrico adquiere este valor si el arreglo es
simétrico, dónde: AM = es la distancia en metros en el terreno entre los electrodos A y M
DISPOSITIVO WENNER
K (metros) = 2πa
Este arreglo puede ser de suma utilidad cuando los valores de voltaje en recepción son muy
bajos, el punto medio entre M y N es el punto de sondeo.
Este difiere a los anteriores debido a que los electrodos de corriente se encuentran
separados de los electrodos de recepción por la distancia en metros denominada “na”
siendo “a” la distancia mutua que existe entre pares de electrodos que conforman los
dipolos, el factor geométrico está dado por:
La medida de recepción del potencial varía con el cuadrado de la distancia de los electrodos
de corriente a medida que la distancia entre el centro de los electrodos de potencial MN y el
electrodo de corriente A va creciendo, se debe también aumentar la apertura de los
electrodos de potencial MN a ¼ a 1/5 de la distancia “a” según el esquema de la fig. -32, el
factor geométrico se calcula a partir de la fórmula: K = (2*π/ ((1/AM) - (1/AN)), donde se
toma a BM=BN=0
El factor geométrico está dado por: K (en metros) =2 π n a (n+1), el punto de medida del
sondeo es el punto medio entre N y A.
DISPOSITIVO LEE
44
Geofísica Aplicada Antonio Carrasco Cortez
En la práctica estos arreglos son poco usados por el hecho que la representación de los
datos correspondientes a estos es relativamente más compleja que la de los arreglos
lineales.
45
Geofísica Aplicada Antonio Carrasco Cortez
Los conceptos de sensibilidad nos ayudaran a elegir la metodología óptima para el registro
de datos en los métodos eléctricos de prospección a continuación, se describe la
sensibilidad y las características de uso de los principales arreglos electrodicos:
ARREGLO WENNER
En la figura 36, se muestra la configuración Wenner, donde los contornos graficados para el
diagrama de sensibilidad son sub horizontales abajo del centro del arreglo electrodico, lo
que nos indica que este dispositivo es sensible a cambios de resistividad vertical del
subsuelo no así para aquellos cambios de resistividad horizontales en general el dispositivo
Wenner es bueno resolviendo cambios verticales de resistividad (estructuras horizontales)
pero relativamente pobre en detectar los cambios horizontales (estructuras verticales
estrechas por ejemplo fallas, filones y otros de buzamiento sub vertical y de poca
potencia).
46
Geofísica Aplicada Antonio Carrasco Cortez
47
Geofísica Aplicada Antonio Carrasco Cortez
este conjunto es muy sensible a cambios de resistividad entre cada par de dipolos, note que
la sensibilidad en el modelo de contorno es casi vertical.
Así el dispositivo dipolo-dipolo es muy sensible a cambios de resistividad horizontal de
modo que es útil en estructuras verticales tales como muros enterrados, cavidades, diques y
plumas de contaminación, pero relativamente pobre en estructuras horizontales como las
capas sedimentarias.
La profundidad media de investigación de este conjunto también depende del factor “n”, así
como del factor “a” en general, este conjunto posee una profundidad más somera de
investigación comparado con el dispositivo Wenner, sin embargo, para sondeos en 2
dimensiones (2D) este conjunto tiene una cobertura mejor de datos horizontales que el de
Wenner.
Una posible desventaja de este conjunto es que la señal es muy pequeña para valores
grandes del factor “n”; el voltaje es inversamente proporcional al cubo del factor “n”, lo
que significa que para la misma corriente, el voltaje medido por el receptor del
resistivímetro disminuye más de 200 veces cuando “n” aumenta desde 1 a 6 para usar este
conjunto el resistivímetro debería tener comparativamente una alta sensibilidad y un muy
buen rechazo de ruido en el circuito receptor además de asegurar un muy buen contacto
entre los electrodos MN y el terreno.
ARREGLO SCHLUMBERGER
48
Geofísica Aplicada Antonio Carrasco Cortez
La configuración polo dipolo tiene una buena cobertura horizontal y también una fortaleza
de señal significativamente más alta comparada con el dispositivo dipolo-dipolo y las otras
configuraciones comunes, el conjunto polo dipolo es un dispositivo asimétrico cuando este
se aplica sobre estructuras simétricas las anomalías de resistividad aparente en la
pseudosección, son asimétricas en algunas situaciones, la asimetría en los valores de
resistividad aparente podría influir en el modelo obtenido después de la interpretación de
datos un método para eliminar el efecto de esta asimetría sería repetir las medicines con un
arreglo de electrodos inverso por combinaciones entre las medidas realizadas en forma
“Directa e "Inversa", cualquier predisposición en el modelo debido a la naturaleza
asimétrica de este conjunto se eliminaría.
El conjunto Polo-dipolo requiere también un electrodo remoto (el electrodo B) que debe
implantarse suficientemente lejos de la línea del sondeo, la fortaleza de señal es
comparativamente menor a los dispositivos Wenner y Schlumberger, pero más alta que en
el dispositivo Dipolo-dipolo.
49
Geofísica Aplicada Antonio Carrasco Cortez
A continuación, los siguientes problemas que reforzaran los conceptos del capítulo.
1.-Se realizó la medida de la resistividad del suelo usando el arreglo electrodico tipo
wenner calcular las profundidades aproximada de investigación si las aperturas de los
electrodos AB son las siguientes 3m, 5m, 8m, 10m, 30m, 1000m
2.-Dedusca la fórmula del factor geométrico tipo Lee a partir de la siguiente formula base:
3.- Se desea realizar el reconocimiento de una estructura sub vertical (filon de sulfuros)
sobre una topografía sub vertical, mediante el método eléctrico que arreglo electrodico es el
más recomendado.
50
Geofísica Aplicada Antonio Carrasco Cortez
Capítulo V
Una unidad geoelectrica es aquella capa o estrato del subsuelo identificada por su
resistividad y espesor, un corte geoelectrico es la sección formada por diferentes unidades
geoelectricas de acuerdo con su disposición en el terreno.
La realización de una serie de sondeos eléctricos a lo largo de un perfil que permite obtener
una sección del subsuelo a lo largo de este correlacionando datos geológicos y geofísicos
simultáneamente de los resultados obtenidos individualmente por cada sondeo eléctrico.
Analicemos el caso de un substrato de dos capas como muestran los gráficos 41 y 42,
ambos análisis intentan explicar el principio básico del sondeo eléctrico desde dos puntos
51
Geofísica Aplicada Antonio Carrasco Cortez
distintos de vista, espero que el lector asimile el principio físico del sondeo eléctrico
vertical (sev) ya sea partir del primer caso o el segundo.
mientras que los electrodos de resepcion MN permanesen fijos a la distacia MN/2 = 1.5
metros.
52
Geofísica Aplicada Antonio Carrasco Cortez
corriente AB, lo suficiente para que la corriente fluya atraves de la vase de la capa 1 y
atrves de la capa 2 de modo que la resistividad se aproxime a la resistividad de la segunda
capa, este proseso se logra en la pocicion AB 2.
Es fácil ver que la resistividad aparente seria, virtualmente la misma que la resistividad de
la primera capa cuando la separación entre los electrodos AB es muy pequeña en
comparación con el espesor de la capa superior, por el hecho que en este caso sería muy
poca la corriente que penetraría en el sustrato, para espaciamientos muy grandes en
comparación con el espesor de la primera capa la resistividad aparente se aproxima a la
resistividad de la segunda capa, porque la parte de la corriente confinada a la capa
superficial o primera capa se hace despreciable.
53
Geofísica Aplicada Antonio Carrasco Cortez
ADQUISICIÓN DE DATOS
En el arreglo Schlumberger las distancias de AB/2 (serie típica de valores de AB/2 =1, 1.5,
2, 3, 6, etc.) cuando la medida de potencial es muy pequeña con relación a la escala en el
milivoltímetro del receptor se incrementa la separación entre los electrodos MN, en este
caso se debe de realizar la labor denominada de empalme, que será descrita a continuación:
Una vez tomada la lectura con una posición de MN-1 para cambiar a la siguiente posición
de MN-2, no se deben de mover los electrodos de corriente debido a que se debe de tomar
la resistividad a las dos posiciones de “MN”, MN-1 y MN-2 con una misma apertura de
electrodos “AB” en teoría la lectura de resistividad deberá ser la misma o algo próxima a
este proceso se le denomina empalme de curva.
54
Geofísica Aplicada Antonio Carrasco Cortez
La aplicación del sondeo eléctrico vertical SEV., exige que las formaciones o capas que se
desean investigar presenten suficiente contraste de resistividad, que haga posible su
diferenciación eléctrica, además que estas formaciones o capas deben tener un buzamiento
sub horizontal no mayor a 3 o 5 grados, esto se debe a que la teoría de interpretación del
sondeo eléctrico solo se aplica a terrenos con estratos sub horizontales.
Otra condición necesaria para el correcto uso del método es que la superficie del terreno sea
relativamente llana, aunque esta condición es menos estricta para los SEV profundos hasta
ahora no han podido establecer procedimientos físico matemáticos para corregir la
55
Geofísica Aplicada Antonio Carrasco Cortez
influencia de los accidentes topográficos a causa de las complicaciones teóricas que entraña
este problema.
Este efecto en general no muy notable puede ponerse fastidioso en las regiones montañosas,
como muestra el grafico 44.
Para la compresión total de este importante acápite me resulta más didáctico presentarla en
forma de un ejemplo:
56
Geofísica Aplicada Antonio Carrasco Cortez
Bueno vale la pena aclarar el rango y el multiplicador son datos que se obtienen en el
receptor del equipo, el producto de ambos nos da la lectura de potencial en milivoltios, la
intensidad de corriente es medida directamente en el transmisor de baja frecuencia.
K= π ((S*S/a) -(a/4))
Ra = (ΔV /I) *K
Posteriormente los datos de resistividad aparente deben de ser ploteados en una hoja de
escala bilogaritmica donde en el eje horizontal se encuentre la separación de electrodos
AB/2 y en el eje vertical el valor correspondiente de resistividad aparente como muestra la
figura 47 la misma corresponde a los datos de resistividad calculados anteriormente.
58
Geofísica Aplicada Antonio Carrasco Cortez
Noten como se forma una curva relativamente suave y sin cambios abruptos de pendiente o
de puntos dispersos el caso de existir puntos dispersos realizados por malas medidas.
Fig-47 Puntos ploteados datos en corriente de baja frecuencia, hoja bilogaritmica para
sondeos tipo Schlumberger.
De la misma manera que el caso anterior, se ejemplificara este paso: nótese que la
metodología de toma de datos, es la misma que en corriente de baja frecuencia, salvo que se
compensa o de no ser así se registra el parámetro denominado ruido geoelectrico, este
parámetro es un voltaje telúrico, que debe de ser registrado por los electrodos de potencial
antes de inyectar corriente por el transmisor en el suelo generalmente son lecturas muy
bajas de voltaje, que varían en función a la posición de los electrodos en el proceso de
medición.
59
Geofísica Aplicada Antonio Carrasco Cortez
60
Geofísica Aplicada Antonio Carrasco Cortez
Para hallar la resistencia “R” en ohmios se divide el diferencial de potencial “ΔVF” en mili
voltios entre la intensidad de corriente “I” en miliamperios.
Para hallar “K” en metros o factor geométrico se usa la siguiente formula tomada del
capítulo V:
K= π ((S*S/a) -(a/4))
a =Distancia MN en (m)
Ra = (ΔVF /I) *K
61
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 50 Puntos ploteados datos en corriente continua hoja bilogaritmica para sondeos tipo
Schlumberger.
Los datos tomados en el receptor son principalmente dos el ruido geoelectrico y la lectura
de potencial producido por el transmisor, como ya se explicó la diferencia de los mismos se
denomina lectura filtrada, el valor de la lectura filtrada por recomendación debe ser mayor
a 10 milivoltios.
Tenga cuidado al tomar la lectura de ruido geoelectrico puesto que este varía con el tiempo,
si existiese una tormenta eléctrica o las corrientes telúricas son tan intensas que no permiten
tener un valor de ruido constante es mejor que abandone el trabajo y vuelva otro día.
En muchas ocasiones una mala lectura se debe a la mala conexión de los electrodos con el
suelo, lo que se recomienda para solucionar esto es clavar muy bien los electrodos al suelo
62
Geofísica Aplicada Antonio Carrasco Cortez
y si es posible moje un poco la zona donde se implanto el electrodo con una solución de
agua salina o sulfato de cobre en solución al 10%.
Verifique que los cables estén bien conectados y su aislado se encuentre en buen estado.
Por otra parte, la utilización de coordenadas logarítmicas tiene una gran ventaja para la
interpretación cuantitativa por medio de curvas patrones, puesto que una curva para un
determinado modelo puede ser empleado para otros modelos cuyas relaciones de espesor y
resistividad son iguales.
2 MAXIMA RESISTIVIDAD
£= √
MINIMA RESISTIVIDAD
En otras palabras si se extienden los electrodos de corriente en una dirección por ejemplo
norte sur y se toma el valor de resistividad a la apertura de AB/2 igual a 5m y se calcula la
resistividad y nos da el valor de 100 ohm*m, luego en el mismo centro del cual se midió el
63
Geofísica Aplicada Antonio Carrasco Cortez
anterior AB/2 se extiende la misma apertura de electrodos de 5 metros pero esta vez en
dirección este oeste y se toma la medida de resistividad y la misma resulta ser de 90
ohm*m entonces se calcula el coeficiente de anisotropía.
Esto quiere decir en la práctica que todos los sondeos eléctricos deben de realizarse en una
misma dirección (extender los electrodos AB en una dirección prederminada, en función a
la situación geológica y el valor de la anisotropía) puesto de no ser así no se podrán
correlacionar los mismos por el efecto de anisotropía.
Ya adelantamos que los datos de resistividad deben de ser ploteados en una hoja
bilogaritmica donde el eje horizontal represente la apertura de electrodos y el eje vertical la
resistividad medida.
i. Las curvas de campo pueden ser suaves y continuas, si estas han sido obtenidas en
un terreno plano y el sub suelo está constituido por estructuras tabulares y
lateralmente homogéneas.
ii. Cuando existen variaciones laterales o variaciones topográficas estas se manifiestan
en las curvas de resistividad aparente como datos dispersos, aunque la dispersión
puede también ser producida por cambios de resistividad debido a la presencia de
unidades geoelectricas de reducido espesor, las variaciones topográficas pueden ser
observadas durante la adquisición de datos y eventualmente las variaciones laterales
cuando estas son superficiales lo cual permite determinar si es conveniente suavizar
o no la curva de resistividad aparente.
-La forma de las curvas es concordante con la resistividad verdadera de las unidades
geoelectricas del suelo, es decir que cuando la resistividad de las unidades es mayor o
menor en profundidad esta se manifiesta con el ascenso o descenso de la curva de
resistividad aparente.
65
Geofísica Aplicada Antonio Carrasco Cortez
El grafico 53, muestra los modelos de dos capas estos son relativamente sencillos a
continuación se muestran los gráficos correspondientes a las variantes de las curvas de 3
capas (fig. 54).
El grafico 54, se toma como ejemplo la curva de tipo H, este se explica de la siguiente
forma: cuando la capa más cercana a la superficie tiene la resistividad aparente mayor a la
segunda capa o capa infrayasente, pero el sustrato infrayasente a las dos primeras capas es
66
Geofísica Aplicada Antonio Carrasco Cortez
de mayor resistividad que la segunda capa, esta tercera, puede ser de mayor o menor
resistividad que la primera, pero su resistividad siempre será mayor a la de la segunda capa.
K = (Ra1<Ra2>Ra3)
H = (Ra1>Ra2<Ra3)
Q = (Ra1>Ra2>Ra3)
A = (Ra1<Ra2<Ra3)
FUNDAMENTOS DE LA INTERPRETACIÓN
Hasta que el problema de la interpretación de los sondeos eléctricos (sev), fue resuelto por
Hummel quien elaboro el modelo de comportamiento eléctrico de un suelo estratificado
asta para 3 capas basándose en la analogía óptica de las imágenes especulares, a este
método se le denomino método de las imágenes.
Los primeros métodos de interpretación cuantitativa del sev., fueron métodos gráficos entre
los que se destacan el método de superposición de Tagg y el del punto auxiliar, con el
advenimiento de los ordenadores electrónicos se desarrollaron numerosos métodos de
naturaleza numérica como los de aproximación y los iterativos, los mismos desarrollados
por el D.Koefoed, vale la pena recalcar que dichos métodos modernos se basan en los
estudios de Maillet y los parámetros de Darr Zarrouk.
67
Geofísica Aplicada Antonio Carrasco Cortez
(Si el lector se ha pasado el trabajo de leer el texto notara que la expresión fue deducida en
el capítulo II del presente texto).
Δv= Z (I*Ra1/2π) * 1 / r1
Donde “r1“es la distancia del observador al punto de reflexión del potencial en la interface
ente las capas y Z es un coeficiente de reflexión el potencial total en el primer medio es, por
consiguiente
Un observador en el segundo medio, vera solo la luz de la fuente verdadera, pero atreves de
la separación solo pasa una parte (1-Z), puesto que la fracción “Z” ha sufrido reflexión, el
potencial en el medio dos puede por tanto escribirse:
Ahora analicemos el caso del potencial en el contacto donde el cual en todos los puntos
para cualquier (r) = distancia desde la fuente hasta el punto el potencial hade ser el mismo
Tanto la expresión (1) como la (2) son iguales al potencial en el punto (p) punto en la
interface entre los dos medios de resistividad Ra1 y Ra2
68
Geofísica Aplicada Antonio Carrasco Cortez
Ra1+Ra1*Z-Ra2+Ra2*Z=0
Ra1*Z+Ra2*Z = Ra2-Ra1
Z= (Ra2-Ra1) / (Ra2+Ra1)
En las primeras etapas de desarrollo del SEV, se pensaba que las formas de las curvas de
resistividad y sus puntos característicos como máximos, mínimos, puntos de inflexión, etc.
Dependían directa y separadamente de las resistividades y espesores de las distintas capas
componentes del corte geoelectrico sobre el que se efectúa el sondeo eléctrico.
Ahora bien, los trabajos de los geofísicos de la escuela francesa en especial los de Raymond
Maillet demostraron con todo rigor que cada curva de SEV no refleja en general el efecto
independiente de espesores y resistividades sino lo que influye decisivamente sobre la
morfología de dichas curvas, son los productos y cocientes de espesor y resistividad de
69
Geofísica Aplicada Antonio Carrasco Cortez
cada capa a estas magnitudes combinadas Maillet asigno el nombre de parámetros de Darr
zarrouk y sus definiciones son las que siguientes:
Es el espesor E de una capa dividida por su resistividad se representa por la letra S entonces
para el manto o capa “i” como índice numérico de manto o capa será:
Si = Ei / Rai
Esta magnitud representa la conductancia de un prisma recto cortado en la capa, cuya altura
es igual al espesor de esta, con base cuadrada de un metro de lado y de la misma
resistividad de la capa.
Ti= Ei * Rai
Una de las propiedades más importantes de estos parámetros es que son aditivos es decir
que la S o la T de un conjunto de capas es igual a la suma de las S o las T, respectivamente,
de todas las capas.
El análisis de los parámetros de Darr zarrouk desde mi punto de vista será más fructífero,
con un ejemplo imaginemos una serie de unidades geoelectricas isótropas y homogéneas de
resistividad Ra y espesor E dada la naturaleza vectorial de la corriente, esta fluirá en
diferentes direcciones si esta corriente se descompone en dos direcciones una paralela a la
estratificación y otra perpendicular, para la segunda componente las unidades se llegan a
T=T1+T2+T3+T4+…………Tn
Tn=∑𝑛1 𝐸𝑖 ∗ 𝑅𝑎𝑖
Del mismo modo para los componentes que fluyen paralelos a la estratificación, las
diferentes unidades se comportan como conductores en paralelo y la conductancia total será
igual a la suma de conductancias de cada una de las unidades, la conductancia total será:
S=S1+S2+S3+S4+……………Sn
Sn=∑𝑛1 𝐸𝑖/𝑅𝑎𝑖
PRINCIPIOS DE EQUIVALENCIA Y DE SUPRECION
La equivalencia en S tiene lugar en capas cuya resistividad mucho menor que las capas
contiguas inferior y superior, cuyo parámetro S es mayor que el del conjunto de capas
supreyasentes.
Ejemplo:
La segunda capa de un corte geoelectrico definido por los espesores E1=10, E2=5 y las
resistividades Ra1=1, Ra=20, Ra3=1 es equivalente a la del corte E1=10 E2=2.5 y Ra1=1,
Ra2=40, Ra3=1 pues en ambos casos T2=100 y S2<<S1.
71
Geofísica Aplicada Antonio Carrasco Cortez
No se debe tomar de forma determinante el análisis del acápite que prosigue, puesto que en
el sentido práctico no suele usarse salvo se requiera analizar alguna variante en el método
de interpretación.
𝛻*E=0
𝛻*B=J…. (1)
E=Ra*J
El potencial asociado a un campo que se rige por la ecuación (1), satisface la ecuación de
Laplace:
𝛻 ∗ 𝛻*V=0
1
(𝜕 2 𝑉/𝜕 𝑟^2) + (( ) 𝜕 𝑉/𝜕 𝑟) + (𝜕 2 𝑉/𝜕 𝑧^2) + ((1/(𝑟𝑟)𝜕 2 𝑉/𝜕 𝑝^2) = 0
𝑟
72
Geofísica Aplicada Antonio Carrasco Cortez
Dónde:
Cualquier combinación lineal de la función (3) también es una solución de la expresión (2)
por lo tanto una forma más general de expresar la solución particular, es tomar en cuenta
todos los posibles valores de G:
∞
V=∫0 (Ж(𝐺 ) ∗ 𝑒 −𝐺𝑧 + ¥(𝐺 ) ∗ 𝑒 +𝐺𝑧 )𝐽𝑜(𝐺𝑟) ∗ 𝑑𝐺 ..(4)
El potencial V en cualquier punto del terreno es una función de (r,z) la expresión (4)
representa el potencial respecto a una referencia muy lejana ,teóricamente en el infinito
Donde Ж y ¥ son funciones que serán sustituidas por nuevas funciones según la siguiente
deducción que considera a un terreno homogéneo
73
Geofísica Aplicada Antonio Carrasco Cortez
u=√𝑟 2 + 𝑧 2
J= I/ (2πU)
El campo eléctrico en cualquier punto del terreno se obtiene por la ley puntual de ohm
E=Ra*J
El potencial en la superficie del terreno se obtiene haciendo z=0, en la ecuación (5) esta
puede escribirse en forma integral utilizando la función de Bessel de orden cero mediante la
integral de Lipschitz:
∞
∫0 𝑒 −𝐺𝑧 ∗ 𝐽𝑜(𝐺𝑟)𝑑𝐺 =1/(u)… (6)
74
Geofísica Aplicada Antonio Carrasco Cortez
∞
V= ((Ra*I)/ (2π)) ∫0 (𝛷(𝐺 ) ∗ 𝑒 −𝐺𝑧 + 𝛸(𝐺 ) ∗ 𝑒 𝐺𝑧 + 𝑒 −𝐺𝑧 )𝐽𝑜(𝐺𝑟) ∗ 𝑑𝐺 … (8)
Donde Ra es la resistividad del estrato más superficial, estrato1, en la ecuación (8) las
funciones Ж y ¥ han sido sustituidas por nuevas funciones 𝛷 y 𝛸, si el terreno es
homogéneo (8) debe de reducirse a la expresión (7) lo que implica que 𝛷 y 𝛸 sean igual a
cero si por el contrario el terreno tiene más de un estrato, la solución (8) corresponde al
potencial en cada estrato del terreno, en general las funciones 𝛷 y 𝛸 son diferentes para
cada estrato .es decir cada estrato individualmente tendrá una distribución de potencial Vi
perteneciente individualmente a cada uno la misma es dependiente de las funciones 𝛷i
y 𝛸i respectivas ,la solución (8) debe de satisfacer las condiciones de contorno impuestas
por la fuente y los limites ,entre los diferentes estratos del terreno estas condiciones son las
siguientes :
. En cada límite entre estratos Z=espesor hi, el potencial debe de ser continuo Vi=V
i+1, lo que se traduce en:
i. En cada limite la componente vertical del vector densidad de corriente debe de ser
continuo, después de ejecutar la operación gradiente sobre la solución para el
potencial se obtiene:
Rai+1 (1+( 𝛷𝑖(𝐺)) ∗ 𝑒 −𝐺ℎ𝑖 - 𝛸𝑖(𝐺) ∗ 𝑒 𝐺ℎ𝑖 ) = Rai (1+( 𝛷𝑖 + 1(𝐺)) ∗ 𝑒 –𝐺ℎ𝑖 –
𝛸𝑖 + 1(𝐺) ∗ 𝑒 𝐺ℎ𝑖
ii. En la superficie del terreno, z=0, i=1, el componente vertical del campo eléctrico i
del vector densidad de corriente deberán de anularse excepto en el punto de
inyección de corriente ya que se considera al aire como un medio aislante casi
perfecto:
Ez=Jz=0... (9)
∞ ∞
∫0 𝐺 ∗ 𝐽𝑜(𝐺𝑟)𝑑𝐺 +∫0 (𝛷 (𝐺 ) − 𝛸(𝐺 ))𝐽𝑜(𝐺𝑟) ∗ 𝑑𝐺 =0
75
Geofísica Aplicada Antonio Carrasco Cortez
𝛷(𝐺) = 𝛸(𝐺)
Cuando z tiende al infinito el potencial deberá de tender a cero, esto implica que la función
𝛸(𝐺) = 0
Particularizando las condiciones de contorno para un suelo de dos estratos se obtiene las
siguientes relaciones:
Donde el subíndice i =1
RAI+1 (1+( 𝛷𝑖(𝐺)) ∗ 𝑒 −𝐺ℎ𝑖 - 𝛸𝑖(𝐺 ) ∗ 𝑒 𝐺ℎ𝑖 ) = RAI (1+( 𝛷𝑖 + 1(𝐺)) ∗ 𝑒 –𝐺ℎ𝑖 –
𝛸𝑖 + 1(𝐺 ) ∗ 𝑒 𝐺ℎ𝑖 …. (12)
Donde el subíndice i =1
𝛷𝑖(𝐺)
Z= (Ra2-Ra1) / (Ra2+Ra1)
El potencial sobre la superficie del terreno respecto a una referencia muy lejana se
determina mediante la ecuación (10) sustituyendo Φi(G) se obtiene:
76
Geofísica Aplicada Antonio Carrasco Cortez
∞
V= ((Ra1*I)/ (2π)) ( ∫0 ((1 + Z ∗ e −2Gh )/ (1 − Z ∗ e −2Gh )) Jo(Gr)dG)… (14)
El factor que involucra a Z y la función exponencial pueden expresarse como una sumatoria
infinita:
De la ecuación (17) podemos obtener el valor de R=VNM/I medido para una disposición i
de los electrodos AB, MN si substituimos el valor de R en la ecuación
Ro = (V /I) *K
77
Geofísica Aplicada Antonio Carrasco Cortez
Capítulo VI
La conciencia no es más que un espejo opaco que refleja las opiniones del medio en que
uno vive.
Los métodos de interpretación gráficos del sondeo eléctrico vertical (SEV), en la actualidad
casi han sido reemplazados por completo por los métodos de interpretación por ordenador
(varios autores de textos en prospección geofísica incluso los considera obsoletos y
anticuados), cabe la pena mencionar que no todos podemos exceder a un software de
ordenador, en el caso mío los primeros trabajos que realice fueron interpretados por el
método del punto auxiliar, donde los resultados pueden llegar a ser tan precisos como los
resultados de un programa de ordenador.
Este método propuesto por Tagg. (1934) es aplicable para interpretar hasta 2 capas, se
considera un método directo ya que la solución se obtiene directamente de un cálculo
sencillo con los datos experimentales, el método se basa que para una determinada distancia
78
Geofísica Aplicada Antonio Carrasco Cortez
79
Geofísica Aplicada Antonio Carrasco Cortez
El primer término de esta ecuación es la expresión que nos da como resultado el valor de la
resistividad aparente cuando se emplea el dispositivo Wenner, por lo tanto, se puede
escribir la expresión de la siguiente manera:
Ra=R1 (1+4F)
Para determinar la profundidad h a la cual se encuentra una interfaz, se realiza una serie de
medidas de la resistividad aparente aumentando la distancia inter electrodica después de
cada observación como cuando el valor de la separación electrodica “a” (en el caso del
dispositivo Wenner) tiende a cero, el valor de Ra, tiende a R1 es posible determinar el valor
de R1 utilizando aperturas interelectrodicas muy pequeñas, para mayores separaciones las
únicas aperturas son h y k ya que Ra se obtiene de los datos observados en cada estación
(estación=puntos de medida o posiciones de A,B) y R1 se ha determinado previamente, por
consiguiente se deduce que al efectuar dos mediciones de Ra correspondientes a dos
distancias interelectrodicas distintas Ra1 y Ra2 sustituyendo los valores en la ecuación (1)
Y resolviendo el sistema de ecuaciones que resulta puede determinarse los valores de h y k
estas son las únicas incógnitas de la función “F”, esta tiene una expresión complicada, no
resulta nada conveniente resolver el sistema por métodos algebraicos, por lo que adopta una
solución grafica para ello se utiliza una serie de curvas patrón calculadas a partir de la
ecuación (1) que también puede ser escrita de la siguiente forma:
2h 2 2h 2 4h 2
Ra=R1 (1+4((k (√(1 + ( a ) ))- (k (√(4 + ( a ) )))+ ((k^2 (√(1 + ( a ) ))- (k^2
4h 2 6h 2
(√(4 + ( ) )))+ ((k^3 (√(1 + ( ) ))- (k^3
a a
6h 2
(√(4 + ( a ) )))+…………………………….)
Tagg ha construido una serie de curvas basadas en la ecuación (1) como muestra el grafico
61.
80
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 61 Curvas patrón para 2 capas según Tagg (Introducción a la Prospección Geofísica
(1984) – Milton B. Dobrin).
81
Geofísica Aplicada Antonio Carrasco Cortez
Una serie para R2/R1 mayor a 1 y otra para R2/R1 menor a 1, cada línea horizontal
corresponde a un valor de Ra/R1, la misma cruzara de una a diez a las curvas que
representan decimas de los valores de “k”.
Al aplicar estas curvas a los valores obtenidos en el campo, se supone que los valores de
“Ra” para valores de apertura de “a” (dispositivo Wenner) que se aproximen a cero es R1 y
dividiendo todas las lecturas de Ra entre R1, la relación Ra/R1 versus los valores de a
correspondientes, darán como resultado una gráfica que se asemeje a una recta que cortara
la serie apropiada de curvas patrón.
En cada intersección con una curva, que corresponda a una k diferente, el valor de h/a es
convertido en “k”, multiplicándolo sencillamente por el valor de la apertura “a” después los
valores de “k” son comparados con los de “h” en otra gráfica, y todos los puntos
correspondientes a una separación de electrodos dada son unidos por una línea.
EL MÉTODO DE SUPERPOSICIÓN
82
Geofísica Aplicada Antonio Carrasco Cortez
83
Geofísica Aplicada Antonio Carrasco Cortez
El método del punto auxiliar permite determinar el corte geoelectrico, haciendo uso de
curvas teóricas de dos o tres capas o bien combinando ambas y cuatro ábacos auxiliares uno
para cada tipo de cortes de tres unidades.
Este método consiste en reducir la curva de resistividad aparente de varias unidades a otra
de menor número, reemplazando dos unidades a un estrato equivalente descomponiendo de
esta manera la curva de campo en tantas curvas de dos unidades, cuantas veces amerite el
caso.
La curva graficada se superpone sobre el ábaco de dos capas hasta que simultáneamente
manteniendo los dos ejes paralelos (ejes de la escala bilogaritmica del papel
semitransparente y del Abaco patrón de 2 capas) hasta que la curva de campo coincida con
alguna curva teórica patrón, o de no conseguir esto se pueda interpolar entre dos curvas
teóricas.
Una vez conseguido esto, se proyecta sobre el papel transparente, de la curva de campo la
cruz (origen de coordenadas del Abaco patrón).
Luego se procede a anotar el valor numérico ubicado en el eje vertical del Abaco, que
coincida o concuase con la parte final de la curva teórica (cola de la curva que tiende a una
recta) elegida por comparación para la interpretación, a este valor numérico le
denominaremos “W” para facilitar la explicación.
84
Geofísica Aplicada Antonio Carrasco Cortez
Para hallar el valor de espesor de la primera capa se debe interpolar el punto cruz de la hoja
bilogaritmica de datos de campo, hacia el eje horizontal que corresponde a los valores de
aperturas AB/2 en metros, donde se obtendrá un valor numérico de espesor en metros que
corresponde al espesor real de la primera capa.
85
Geofísica Aplicada Antonio Carrasco Cortez
Para la interpretación de tres capas, se realiza el artificio de sustituir las dos primeras capas
por una capa ficticia e equivalente a las dos capas sustituidas, la capa ficticia tendrá una
resistividad y un espesor equivalentes que serán dados por las abscisas y ordenadas de la
segunda cruz a continuación se explica el proceso de interpretación con un poco más de
detalle:
Se superpone la parte izquierda de la curva de campo, con las curvas teórica hasta encontrar
la que mejor ajuste del Abaco correspondiente a las curvas patrón de dos capas y se
proyecta la primera cruz (origen de coordenadas del Abaco patrón), ojo solo la parte
izquierda de la curva de campo debe coincidir con la curva patrón, no toda la curva de
campo.
86
Geofísica Aplicada Antonio Carrasco Cortez
Luego se procede a anotar el valor numérico ubicado en el eje vertical del Abaco, que
coincida o concuase con la parte final de la curva teórica (cola de la curva que tiende a una
recta) elegida por comparación para la interpretación, a este valor numérico le
denominaremos “W”.
Para hallar el valor de resistividad de la primera capa se procede a interpolar el punto cruz-
1 de la hoja bilogaritmica de datos de campo, hacia el eje vertical, correspondiente al eje de
resistividades, donde se obtendrá un valor numérico de resistividad que corresponde a la
resistividad verdadera de la primera capa.
Para hallar el valor de espesor de la primera capa se debe interpolar el punto cruz-1 de la
hoja bilogaritmica de datos de campo, hacia el eje horizontal que corresponde a los valores
de aperturas AB/2, donde se obtendrá un valor numérico de espesor en metros que
corresponde a la potencia real, de la primera capa.
Posteriormente se coloca la primera cruz (la marcada en la hoja de campo), sobre el origen
de coordenadas de los gráficos auxiliares, (estos son de cuatro tipos: TIPO H, TIPO Q,
TIPO K, TIPO A, el uso de cada tipo depende de la morfología de la curva de campo es
decir si la curva de campo es de tipo H, entonces se adoptará el grafico auxiliar tipo H para
realizar la interpretación.
Después de poner la primera cruz sobre el origen de coordenadas de los gráficos auxiliares,
teniendo la precaución de tener siempre los ejes paralelos se proyecta, sobre la hoja
biligaritmica de campo la curva de los gráficos auxiliares, que parta de la primera cruz y
que corresponde al valor numérico “W”.
87
Geofísica Aplicada Antonio Carrasco Cortez
88
Geofísica Aplicada Antonio Carrasco Cortez
Para obtener el valor de espesor de la segunda capa se multiplica el valor numérico de “P”
por el espesor de la primera capa.
Como generar curvas teóricas para su uso en el método de superposición o el método del
punto auxiliar.
89
Geofísica Aplicada Antonio Carrasco Cortez
Paso 2.- Una vez ejecutado el programa elija la pestaña opciones de sistema (Options
sistem) de la barra superior de tareas, se desplegará una ventana en la misma elija la opción
iniciar un nuevo modelo (Begin a new model), como muestra el grafico 71.
90
Geofísica Aplicada Antonio Carrasco Cortez
Paso 3.- Luego en la misma ventana iniciar un nuevo modelo en la opción rango de
evaluación de AB/2 (range of AB/2 values), llene el espacio denominado mínimo con el
valor de AB/2 mínimo con la cual se tomaron los datos que se desean interpretar, y el
espacio denominado máximo con el valor de AB/2 máximo con la cual se ejecutó el sondeo
a interpretar, finalmente llene el espacio denominado puntos con el número de datos de
resistividad aparente tomados en el campo, finalmente presione aceptar (ok) como
muestra el grafico - 73.
Paso 4.- Posteriormente se desplegará la ventana de trabajo del programa, luego ingrese a
la pestaña opciones de pantalla (display options) y elija la opción tamaña de hoja de trabajo
(set the display size), y marque la opción 100% y presione aceptar (ok), como muestra el
grafico 74.
Inmediata mente aparecerá en su pantalla la curva teórica que usted género como muestra
el grafico 75.
5.- Finalmente utilice la pestaña opciones de sistema (system options) para realizar la
impresión de la curva teórica generada.
92
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo VII
(Agustín Aspiazu)
En los años 1950 a 1960 se publicaron muchos artículos atacando a los métodos gráficos de
interpretación en geoeléctrica, los procedimientos correctos, sin embargo, no estaban muy
adelantados por aquel entonces pues, aunque Stefanescu (1930) había resuelto el problema
del comportamiento eléctrico para medios estratificados, su cálculo numérico era lento y
tedioso que asía que este sea difícil de aplicar.
La aparición de los ordenadores produjo cambios enormes en esta situación, por una parte,
se pudieron calcular problemas directos de diversos métodos primero para casos sencillos y
luego para otros cada vez más complejos, también fue aumentando la velocidad de cálculo
de modo que muchos se realizan ahora en tiempos muy cortos.
93
Geofísica Aplicada Antonio Carrasco Cortez
94
Geofísica Aplicada Antonio Carrasco Cortez
-Luego se abre una ventana que corresponde al proceso de instalación del programa como
muestra el grafico 77, simplemente se debe presionar siguiente y aceptar a todas las
opciones que indique la ventana, luego el programa se instalara automáticamente, los
requerimientos mínimos de instalación son: un ordenador Pentium III con 512 Megabits de
memoria RAM y sistema operativo Windows 2000.
95
Geofísica Aplicada Antonio Carrasco Cortez
El registro de datos fue realizado con un equipo de corriente continua, sin corrector de
ruido, en el mismo se llega a calcular, en el programa Excel, hasta la resistividad aparente
(siguiendo la metodología indicada en el capítulo V el sondeo eléctrico; subtitulo registro
de datos en corriente continua).
Paso 1.- Abra el programa Ipi2win ingresando a la barra de menú todos los programas.
Paso 2.-Ingrese al menú archivo y elija la opción crear como muestra el grafico 79.
Fig. 79
Paso 3.-Luego aparecerá una ventana de nombre nuevo punto de sev (new ves point), elija
el nombre del arreglo electrodico usado en el campo para su registro de datos, en la opción
que tiene por defecto al arreglo Sclumberger, como muestra el grafico 80:
Fig. 80
Paso 4.-En la misma ventana new ves point, presione la opción simbolizada con rho (p) y
seleccione la opción mostrar puntos (show numbers).
- Luego ingrese los datos de apertura de electrodos AB/2 debajo de la columna AB/2 siendo
el caso de haber elegido el arreglo Sclumberger.
- Ingrese los datos de apertura de MN bajo la columna denominada MN, tanto los datos de
AB/2 y MN deberán estar en la unidad de metros.
- Ingrese el dato de resistividad aparente calculada bajo la columna RO_a como muestra el
grafico 81:
96
Geofísica Aplicada Antonio Carrasco Cortez
- Luego presione la opción “ok” y se abrirá una ventana donde usted podrá elegir donde su
guardara el archivo que está creando, elija un lugar donde guardar y nombre, para
posteriormente cerrar la ventana.
97
Geofísica Aplicada Antonio Carrasco Cortez
- Entonces se abrirá una ventana denominada opciones, elija la pestaña nuevo modelo esta
le proporcionará espacios para ingresar los siguientes datos: el error mínimo de ajuste de
curva, el error máximo de ajuste y el número de capas para el cual será generada la curva
teórica en función a las características geológicas y a la morfología de la curva de
resistividad aparente se deben eligen estos parámetros según criterio de interpretación.
Paso 6.-Ingrese de nuevo a la opción punto sev y presione la opción nuevo modelo,
inmediatamente el programa encontrara la curva de menor error de ajuste para la curva de
campo con las características indicadas en la ventana de opciones, como muestra la fig. 84
Fig. 84
98
Geofísica Aplicada Antonio Carrasco Cortez
Terminando así el análisis del sondeo eléctrico vertical con la siguiente ventana (fig-85)
Fig. 85
- Las líneas en gradiente color azul representan los parámetros de la curva teórica.
- La curva de color rojo que se superpone a los puntos de campo representa a la curva
teórica.
99
Geofísica Aplicada Antonio Carrasco Cortez
Los valores numéricos obtenidos en la tabla superior son los datos finales de la
interpretación por ordenador, los mismos al igual que los métodos gráficos, corresponden a
datos de resistividad verdadera y espesores de capas.
A primera apreciación pareciese que interpretar un sondeo eléctrico por ordenador fuese
sencillo, salvo que uno debe de tomar en cuenta que los valores de resistividad y los
espesores hallados, deben de concuasar con la geología conceptual de la zona de no ser así
se tendrá que volver iterar el modelo por el ordenador hasta que esto suceda, por volver a
iterar el modelo me refiero a seleccionar el número de capas de la curva teórica, los errores
máximo y mínimo tal vez suavizar la curva de campo manualmente o por medio del
ordenador, ajustar la curva teórica mediante la interfaz gráfica, hasta lograr un resultado
que encaje perfectamente con el modelo conceptual geológico de la zona de no ser así la
interpretación seria incorrecta.
Bueno no quiero asustar a nadie, pero llegar a la interpretación correcta y apropiada a una
determinada situación litológica estructural, a beses puede llegar a tardar horas iterando
modelos que satisfagan al modelo geológico conceptual.
6.- La opción crea un modelo automático de curva teórica para los puntos de campo.
100
Geofísica Aplicada Antonio Carrasco Cortez
7.-Esta opción crea un nuevo modelo de curva teórica ajustada a los parámetros de
inversión previamente definidos en la ventana parámetros de inversión.
14.-Ingreso al bloc de notas, por el mismo se insertan los datos para la creación de perfiles
de resistividad, o correlación de puntos sev.
17.-Es una pestaña que despliega, opciones para modificar manualmente por medio del
ratón los parámetros de la curva teórica iterada.
101
Geofísica Aplicada Antonio Carrasco Cortez
Fig-87
El modelado directo se realiza comparando los datos de campo con la curva obtenida para
un modelo teórico dado asistido por un ordenador y empleando un algoritmo adecuado para
este fin.
Este proceso se realiza hasta que la curva del modelo teórico, coincida con la curva de
campo, ósea que prácticamente se emplea el mismo procedimiento que con el método
grafico de superposición, solamente que en este caso la comparación se la realiza en la
pantalla del monitor del computador.
102
Geofísica Aplicada Antonio Carrasco Cortez
La inversión de los datos se puede realizar mediante los métodos de capa simple y el de
multicapas, el de capa simple es el método tradicional y se obtienen resultados similares a
los de interpretación gráfica, es decir el sub suelo dividido en pocas unidades geoelectricas
que corresponden a otros tantos estratos geológicos.
103
Geofísica Aplicada Antonio Carrasco Cortez
Paso 1.- Se debe de proceder al creado, de la planilla de datos mediante el sub programa bloc
de notas, tomando en cuenta el siguiente formato:
Nombre xxxxxx
Comentario xxxxx
Punto xxxxxx
AB/21 Ra1
. .
. .
. .
. .
AB/2 n Ra n
104
Geofísica Aplicada Antonio Carrasco Cortez
Fig-90
El archivo creado deberá de ser guardado en el directorio que uno desee, con el cuidado de
guardarlo con la extensión .dat para que así el programa SEV-1.3 pueda reconocerlo.
Y elija la carpeta denominada ves. exe como muestra el grafico siguiente (fig-91).
Fig-91
Paso 3.-Una vez ejecutado el programa elija la pestaña opciones de sistema (Options
sistem) de la barra superior de tareas, se desplegará una ventana en la misma elija la opción
de ventana iniciar un nuevo modelo (Begin a new model), como muestra la figura-92.
105
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 92
Fig. 93
Simultáneamente se abre una ventana para asesar al directorio o discos de la PC, se debe
buscar en los mismos el archivo creado con los datos del sondeo eléctrico, es decir el
archivo con la extensión .dat, como muestra el grafico 94.
Fig. 94
Presione la opción abrir, inmediatamente ingresara a la ventana campo de datos (field data),
verifique en la misma que los datos sean los correctos, y presione aceptar (ok) como
muestra el grafico 95.
106
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 95
Paso 4.- en la ventana de trabajo del programa ara su aparición una sub ventana
denominada interpretación automática de 2 y 3 capas, elija en la opción número de capas
(n layer) e ingrese el número de capas adecuado para la interpretación en función a la
morfología de la curva de resistividad aparente, luego presione aceptar (ok) como muestra
el grafico 96
Fig. 96
Paso 5.-El programa ajustara una curva teórica a la curva de campo, luego ingrese a la
pestaña opciones de pantalla (Display options) y elija la opción juego de colores del
esquema (Set de color scheme) y elija la opción color en la ventana desplegada y
finalmente presione aceptar (ok), como muestra el grafico 97.
Fig. 97
107
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 99
En el mismo es posible modificar la curva teórica, para que la misma logre un mejor ajuste
con los datos de campo.
1.-Contiene las opciones iniciar nuevo modelo, guardar imprimir, salir del programa.
2.-Se encuentra dentro de la pestaña editar modelo, y sirve para añadir capas al modelo.
5.- Se encuentra dentro de la pestaña editar modelo, y sirve para eliminar capas al modelo.
3.- Se encuentra dentro de la pestaña editar modelo, y sirve para invertir modelos de curva
de resistividad.
4.- Se encuentra dentro de la pestaña editar modelo, y sirve para abrir la ventana de
interpretación automática.
8.-Abre la ventana donde se guarda los datos del fabricante del programa.
Al igual que con el programa ipi2win la interpretación con el programa SEV 1.30, se
encuentra condicionada a ajustar un modelo de curva cuyos parámetros cumplan con un
modelo geológico determinado, este programa no interpreta datos, por el método multicapa.
El programa sev 1.30 es libre o gratuito se encuentran en el buscador de internet, solo tenga
cuidado que el programa descargado no sea una versión de prueba o su algoritmo base este
alterado, compruebe el correcto funcionamiento del programa antes de usarlo (el presente
texto cuenta con un cd de respaldo en su versión física.
109
Geofísica Aplicada Antonio Carrasco Cortez
110
Geofísica Aplicada Antonio Carrasco Cortez
-Luego se abre una ventana que corresponde al proceso de instalación del programa como
muestra, se debe presionar siguiente y aceptar a todas las opciones que indique la ventana,
luego el programa se ejecutara automáticamente.
Fig. 102 Planilla de datos y cálculos para obtener la resistividad aparente, de un sondeo
eléctrico vertical con la configuración Sclumberguer, sin empalme de curva.
El registro de datos fue realizado con un equipo de corriente continua, sin corrector de
ruido, en el mismo se llega a calcular, en el programa Excel, hasta la resistividad aparente
(siguiendo la metodología indicada en el capítulo V el sondeo eléctrico; subtitulo registro
de datos en corriente continua).
111
Geofísica Aplicada Antonio Carrasco Cortez
Paso 0.- Abra el programa ingresando a la barra de menú todos los programas.
Paso 1.- Ingrese al menú file y elija la opción VES note pad como muestra el grafico 103.
Fig. 103
Paso 2.- Luego aparecerá una ventana de nombre VES notepad, elija el nombre del arreglo
electrodico usado en el campo para su registro de datos, en la opción que tiene por defecto
al arreglo Schlumberger, como muestra el grafico 104:
Fig. 104
Paso 3.- Ingrese los datos de apertura de electrodos AB/2 debajo de la columna AB/2
siendo el caso de haber elegido el arreglo Sclumberger.
- Ingrese los datos de apertura de MN bajo la columna denominada MN, tanto los datos de
AB/2 y MN deberán estar en la unidad de metros.
112
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 105
Fig. 106
113
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 107
Paso 5.- Ingrese de nuevo a la ventana de trabajo del programa y presione la opción “run
inversión”, inmediatamente el programa encontrara la curva de menor error de ajuste para
114
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 109
-terminando así el análisis del sondeo eléctrico vertical con la siguiente ventana (fig-110)
Fig. 110
- En la ventana Model editor, los puntos representan los datos tomados en el campo.
115
Geofísica Aplicada Antonio Carrasco Cortez
- En la ventana Model editor, las líneas en gradiente color rojo representan los parámetros
de la curva teórica.
- En la ventana” Model editor”, la curva de color rojo que se superpone a los puntos de
campo representa a la curva teórica.
Los valores numéricos obtenidos en la tabla superior son los datos finales de la
interpretación por ordenador, los mismos al igual que los métodos gráficos, corresponden a
datos de resistividad verdadera y espesores de capas.
Uno debe de tomar en cuenta que los valores de resistividad y los espesores hallados, deben
de concuasar con la geología conceptual de la zona de no ser así se tendrá que volver iterar
el modelo por el ordenador, hasta que esto suceda, por volver a iterar el modelo me refiero
a seleccionar el número de capas de la curva teórica, los errores máximo y mínimo tal vez
suavizar la curva de campo manualmente o por medio del ordenador, ajustar la curva
teórica en el ordenador, hasta lograr un resultado que encaje perfectamente con el modelo
conceptual geológico de la zona, de no ser así la interpretación seria incorrecta.
File; ves notepad. -Abre la ventana crear un nuevo punto de sondeo eléctrico
116
Geofísica Aplicada Antonio Carrasco Cortez
. -La opción crea un modelo automático de curva teórica para los puntos de campo.
. -Esta opción crea un nuevo modelo de curva teórica ajustada a los parámetros de
inversión previamente definidos en la ventana “Program setub”.
Options; set / line cordenates. -Opción ventana de trabajo para la inserción de datos
topográficos en la creación de perfiles por correlación de puntos de sondeo eléctrico.
117
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo VIII
Por esta razón el hombre penetra con los ojos de la apreciación y el juicio dentro de las
montañas y puede ver casi exactamente los lugares donde existe mena y su cantidad, dirige
la excavación hacia ellos, pues de otra manera se iría al azar
V. Birinnguccio
La primera depende de las circunstancias del problema, si se busca una capa de gravas
dentro de arcillas, se produce un contraste de resistividad, que da lugar a la aplicación del
método eléctrico, si el objeto es una capa «dura» y sana bajo un recubrimiento de roca
alterada y meteorizada, se origina un contraste en la velocidad sísmica, lo que hace
adecuado el método de refracción sísmica, si se buscan restos de cerámica y objetos
antiguos, se tiene un contraste de susceptibilidades magnéticas que nos lleva a la
magnetometría, etc.
Lo fundamental es la existencia de un contraste de alguna propiedad física entre el objeto
buscado y el medio circundante, elegido el método es necesario determinar la distribución y
densidad de las mediciones de campo.
Si esta planificación es correcta, los datos de campo reflejarán la presencia y posición del
objeto buscado.
La interpretación físico matemático es la etapa más delicada, que debe efectuarse por
geofísicos experimentados y con buen conocimiento de las leyes físicas que intervienen en
el método empleado.
118
Geofísica Aplicada Antonio Carrasco Cortez
Los datos de campo se representan normalmente mediante un gráfico, como la curva del
sondeo eléctrico vertical que, aunque parece relativamente simple contiene mucha
información, el problema de la interpretación consiste principalmente en partir de estos
gráficos sencillos para llegar a deducir un modelo de subsuelo, que en muchos casos puede
llegar a ser compleja ya que la naturaleza suele ser más complicada de lo que uno puede
llegar a visar.
Para efectuar una interpretación correcta es necesario conocer con el mayor detalle posible
el problema que se trata de resolver, las leyes físicas implicadas y la variación de la
respuesta del objeto buscada según sea su disposición en el subsuelo, también es preciso
coordinar toda la información geofísica obtenida para evitar discordancias, considerando
también la posibilidad de su existencia lo que puede ocurrir dada la naturaleza geológica de
cada caso en específico.
El impacto de los ordenadores sobre la prospección geofísica muestra que entre los
programas disponibles los hay de «interpretación automática», cuyo mal uso puede llevar a
grandes errores.
Además, no debe usarse ningún programa que no tenga «contraprueba», es decir, que una
vez hallado un modelo solución, no calcule los datos de campo correspondientes y
compruebe su identidad con los realmente observados.
119
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo IX
CALICATAS ELÉCTRICAS
120
Geofísica Aplicada Antonio Carrasco Cortez
Los electrodos de corriente se mantienen separados por una distancia moderada, las de
potencial se encuentran situadas en el centro de aquellas a poca distancia entre sí, el
conjunto entero se mueve a lo largo de un perfil como en el método de Wenner, desde cada
estación a la siguiente avanzando todos los electrodos la misma distancia a la vez.
Tanto el dispositivo wenner como el Schlumberger se caracterizan por presentar una pobre
sencibilidad a los cambios de resistividad en el sentido lateral, ademas de no presentar
grandes ventajas en la ejecucion practica de las calicatas (vease capitulo IV del texto).
121
Geofísica Aplicada Antonio Carrasco Cortez
Este arreglo es muy utilizado en la prospección por calicatas, por el hecho que permite
modificar rápidamente la profundidad de penetración de la calicata, también es un
dispositivo sumamente sensible a cambios de resistividad lateral.
Luego se interpolan los valores de resistividad obtenidos en cada nodo o punto de medida
para formar curvas de iso resistividad y así construir la pseudoseccion eléctrica.
122
Geofísica Aplicada Antonio Carrasco Cortez
123
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 116 Esquema del cateo eléctrico cuando los electrodos de corriente A y C están
conectados al transmisor (arreglo polo dipolo).
124
Geofísica Aplicada Antonio Carrasco Cortez
125
Geofísica Aplicada Antonio Carrasco Cortez
Con frecuencia la línea de inyección de corriente puede ser elegida entre amplios límites y
en la elección intervienen las razones de índole práctico siguientes:
- Si se trata de un dispositivo móvil, cuanto menor sea su longitud más fácil resulta los
desplazamientos entre estaciones (caso dispositivo Wenner, Schlumberger, etc.)
126
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 119 Variación del perfil de resistividad aparente, en función a la apertura de los
electrodos de corriente según Kunetz.
-Cilindro horizontal. - Sea un cilindro horizontal de resistividad nula o infinita que yace
en un medio de resistividad =Ro y si “r” es el radio del cilindro y p la profundidad de su
127
Geofísica Aplicada Antonio Carrasco Cortez
eje, la anomalía que produce en una calicata realizada sobre la vertical de su eje es
aproximadamente igual a:
Se observa que un cilindro cuyo techo yace a una profundidad igual a su diámetro produce
una anomalía del 22% del valor de fondo.
128
Geofísica Aplicada Antonio Carrasco Cortez
En los gráficos 121 y 122, se han reproducido curvas de calicata obtenidas sobre una falla
vertical que pone en contacto dos medios de diferente resistividad realizados sobre el medio
conductor y resistivo, en estos casos se observa en los perfiles que la resistividad aparente
no alcanza el valor máximo de resistividad hasta que el ultimo electrodo de corriente no
haiga cruzado el contacto litológico.
-Capas delgadas:
129
Geofísica Aplicada Antonio Carrasco Cortez
Este método es análogo a los métodos gráficos de interpretación del sondeo electrico,
puesto que está basado en la construcción de perfiles teóricos, que se encuentran en función
a valores de resistividad, el proceso consiste en realizar el ajuste del perfil de campo al
perfil calculado, posteriormente se asumen los parámetros de resistividad (con los cuales
fue construido el perfil teórico), como parámetros de resistividad real de la calicata
eléctrica.
Fórmulas para la construcción del perfil teórico en la configuración polo dipolo (dispositivo
AMN o BNM)
130
Geofísica Aplicada Antonio Carrasco Cortez
x=distancia en metros entre el electrodo de corriente (A) y el punto medio entre los
electrodos de potencial
-Los valores de Z12, X, d y Ro deberán ser iterados hasta que se logre la creación de un
perfil que ajuste al perfil de campo.
-La expresión (1) es exacta para valores de “d” mayor o igual a “X”, la expresión (2) para
“d” menor igual a “X” y mayor igual a cero, las expresiones que define el valor de “Ra”
correspondientes al caso en que todo el dispositivo se encuentra sobre la superficie del
medio de resistividad igual a “Ra” se puede hallar por la formula (1) sustituyendo en ella:
Ro1 por Ro2 y Z12 = (Ro2-Ro1) / (Ro2+Ro1) por Z21 =(Ro1-Ro2)/ (Ro1+Ro2).
Se calcula dos valores de resistividad “Ra”, uno para la configuración AMN y otro para la
configuración BNM a partir de las siguientes formulas
Se calcula por simple media aritmética el valor de Ra AMNB Donde Ra AMN y Ra BNM
son los valores de resistividad calculados a partir de las expresiones 1,2 y 3
131
Geofísica Aplicada Antonio Carrasco Cortez
Las gráficas teóricas son construidas con el eje de resistividades a escala logarítmica, y el
eje de progresivas en el perfil a escala aritmética como muestra el grafico 125, por tanto,
para su comparación la curva de campo también será graficada con la escala logaritmica de
resistividades, y a la misma escala horizontal de la gráfica generada.
132
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo X
o POTENCIAL INDUCIDO
No todos aceptan que un Dios haya hecho al hombre; lo que nadie niega es que el hombre haya
hecho un millar de dioses.
(Franz Tamayo)
La polarización como fenómeno físico puede ser originado debido a dos causas principales:
la polarización electrónica y la polarización de membrana a continuación, serán descritas
con un poco de más detalle:
133
Geofísica Aplicada Antonio Carrasco Cortez
Si la zona cargada positivamente persiste, repele a suficiente distancia otros iones positivos,
y actúa como una membrana impenetrable que impide el movimiento de las mismas atreves
de las capilares, cuando se hace pasar una corriente eléctrica atreves de las arcillas los iones
positivos se desplazan y al interrumpirse la corriente las cargas positivas vuelven a su
posición inicial, el proceso de redistribución se manifiesta como una tención decreciente en
el contacto de electrones con la arcilla ,el cual se traduce en el fenómeno de polarización
inducida.
134
Geofísica Aplicada Antonio Carrasco Cortez
La adquisición de datos se realiza una vez establecidos los parámetros concernientes a los
perfiles de medida (número de perfiles extensión de los mismos y correspondiente
dirección, la distancia entre puntos de medida o estaciones, arreglo electrodico, etc. deberá
de ser elegida con mucho cuidado en base a toda la información disponible para este fin
(mapas geológicos, topografía, otros trabajos de geofísica en la zona, etc.).
Ósea se emplea un sistema simétrico cuya distancia mínima entre los electrodos BM es
igual a la separación de los dipolos AB y MN, la distancia BM (na) se incrementa
proporcionalmente a la distancia AB denominada “a”
Se debe emplear una serie de electrodos para facilitar el trabajo de campo, con lo cual en la
práctica lo que se desplaza son los cables de conexión y no así los electrodos.
También suelen emplearse, aunque con menos frecuencia los arreglos Polo polo, Polo
dipolo y Wenner.
Como ya se recalcó la técnica consiste en inyectar corriente continua al sub suelo a través
de los electrodos A y B, durante un cierto intervalo de tiempo (tiempo suficiente para que el
potencial primario induzca el potencial secundario) y después medir la tención por medio
de los electrodos de potencial MN durante un tiempo “t” antes que el potencial secundario
tienda a cero (fig-130).
PI = polarización aparente.
PI = Vt(v) / Vo (mv)
Otra cantidad que caracteriza el fenómeno de polarización inducida es el área que encierra
la curva de potencial residual y su respectiva asíntota que tiende a cero, sin embargo, en la
práctica solo se mide el área entre un intervalo de tiempo t1 y t2 por integración
automática, la cantidad medida se denomina cargavilidad aparente o simplemente
cargavilidad y la designan por la letra “m”.
t1
m=∫to Vt dt/Vo
138
Geofísica Aplicada Antonio Carrasco Cortez
Se puede calcular la intensidad de corriente en el proceso de carga del modo siguiente sea
“ñ” el valor mínimo admisible de Vt, entonces para el valor medio dado de μ (polarización
aparente), la diferencia de potencial del campo polarizarte es:
Vo= (Ra*I) /K
Reemplazando se tiene:
μ *Vo * 100 ≥ ñ
μ *((Ra*I) / K)*100 ≥ ñ
I ≥ (ñ*K) / (100*Ra* μ)
I = (ñ*K)/ (100*Ra* μ)
139
Geofísica Aplicada Antonio Carrasco Cortez
Por ejemplo: una línea de perfil eléctrico de a=200 metros y n=5, se toman medidas con un
instrumento que toma como valor mínimo de Vt = 9 milivoltios, tomando en cuenta que el
valor medio de resistividad de zona de estudio es de 200 ohm*m y el valor de la
polarización aparente regional es del 2%, entonces se calcula la intensidad de corriente
requerida en la prospección:
I=71.2513 miliamperios
Claro que estos valores solo son un parámetro puesto que en la realidad se procede a elevar
el voltaje de salida del transmisor hasta obtener una respuesta de potencial inducido, mayor
a 0.5 milivoltios, tomando la lectura del mismo al tiempo mínimo de medida (por ejemplo
0.3 segundos), una vez conseguido esto la intensidad de corriente que registre el transmisor
será la adecuado para la prospección procurar nunca hacer trabajar al equipo a la máxima
capacidad de salida del transmisor.
-Tiempo de medida “t” después del corte de corriente para la medición de Vt.: Para el
cálculo del mismo, se procede a tomar lecturas de Vt a distintos tiempos (por ejemplo, a 0.5
segundos el valor de Vt es 10 milivoltios, a 1segundo Vt es igual a 5 milivoltios, a 1.5
segundos el valor de Vt es de 2 milivoltios, a 2 segundos el valor de Vt es de 1 milivoltio y
así sucesivamente hasta que se forme la asíntota de la curva de tención residual).
140
Geofísica Aplicada Antonio Carrasco Cortez
De acuerdo a la figura 132, que representa la curva de tención residual se deberán tomar
como tiempos de medida para el registro de Vt a los valores de “t” que no formen la
asíntota de la curva es decir a los valores de 0.5 segundos, 1 y 1.5 segundos, estos son los
tiempos óptimos de medida de Vt para la prospección.
141
Geofísica Aplicada Antonio Carrasco Cortez
I RUIDO Vo I K R Ra Vt
n A B M N (A ) SP(mv) LECTURA (m V) (m A) ( m) (ohm*m) (ohm*m) 1 - 1,5 s Vt μ%
1 0 10 20 30 0,35 37,1 156 118,9 350 188,5 0,34 64,03462 39,9 2,8 2,35
2 0 10 30 40 0,37 38,9 75 36,1 370 753,98 0,1 73,56419 40.1 1.2 3.32
3 0 10 40 50 0,37 37,3 53,2 15,9 370 1885 0,04 81,00212 38,5 1,2 7,55
4 0 10 50 60 0,37 38,9 48,3 9,4 370 3769,9 0,03 95,77610 39,8 0,9 9,57
5 0 10 60 70 0,39 37,2 45 7,8 390 6597,3 0,02 131,9468 39,8 2,6 33,3
6 0 10 70 80 0,3 39,7 42,4 2,7 300 10556 0,01 95,00174 41 1.3 48.1
7 0 10 80 90 0,4 37,9 42,2 4,3 400 15834 0,01 170,211455 39,2 1,3 30,2
8 0 10 90 100 0,39 34,5 36,3 1,8 390 22619 0.01 104,397519 35,21 0,71 39,4
9 0 10 100 110 0,4 32 33,7 1,7 400 31102 0.01 132,182483 32,5 0,5 29,4
10 0 10 110 120 0,44 38,9 37,4 1,5 440 41469 0.01 141,37164 39.5 0.6 40
11 0 10 120 130 0,44 48,4 49 0,6 440 53910 0.01 73,5132528 48,6 0,2 33,3
12 0 10 130 140 0,44 48,4 49 0,6 440 68612 0.009 93,5623217 48,6 0,2 33.3
13 0 10 140 150 0,45 39,5 39,8 0,3 450 85765 0.009 57,1769744 39,7 0,2 66,7
Fig. 134 Fragmento de planilla de dados de polarización aparente con medida simultanea de
resistividad aparente, las medidas se realizaron con el arreglo electrodico Dipolo dipolo.
Vt = SP –Vt1
Los datos bajo las columnas ABMN representan la posición de los electrodos en un perfil
lineal, realizados con el arreglo dipolo dipolo.
142
Geofísica Aplicada Antonio Carrasco Cortez
Ra = (Vo / I) *K
Otro parámetro que se emplea en el dominio de las frecuencias es el factor metálico (FM)
que se obtiene dividiendo el efecto de frecuencia por la resistividad aparente que
corresponde a la frecuencia menor o corriente continua y multiplicando esto por un factor
constante arbitrario 2* π *(10^5) con fin que el valor resultante no sea demasiado
pequeño, a continuación la expresión que identifica al FM:
143
Geofísica Aplicada Antonio Carrasco Cortez
144
Geofísica Aplicada Antonio Carrasco Cortez
intersección de dos líneas imaginarias trazadas desde el centro de cada dipolo a 45 grados
respecto de la horizontal en dirección del dipolo opuesto como muestra la figura 136.
L= (a/2) + (na/2)
Dónde: a = distancia en metros entre los electrodos AB = distancia entre electrodos MN,
Una vez obtenida la profundidad aparente y la posición del punto de medida se procede a
construir la pseudoseccion, que se elabora interpolando los valores de los puntos de
medida, al final de dicho proceso se obtienen 2 pseudosecciones una de resistividad
aparente y otra de polarización aparente, para otros arreglos electrodicos (Slumberguer,
Wenner) el punto de medida es asignado al centro de los electrodos de potencial y la
profundidad aparente debe de ser calculada por los procedimientos ya explicados en el
capítulo IX.
145
Geofísica Aplicada Antonio Carrasco Cortez
146
Geofísica Aplicada Antonio Carrasco Cortez
Los mapas de isolineas de PI pueden ser elaborados a partir de los perfiles de y las pseudo
secciones de PI, mediante la simple correlación por interpolación del valor de polarización
y su respectiva ubicación en el espacio (coordenadas N, E) estos, permiten definir mejor la
anomalía y la forma del cuerpo y además clasificar las anomalías de acuerdo a su
intensidad.
La configuración de las pseudosecciones, si bien nos pueden proporcionar una idea acerca
de la distribución vertical de los parámetros de PI., medidos no nos dan una imagen cabal
de estas, debido a que normalmente las zonas anómalas presentan un ensanchamiento, para
mayor separación de electrodos, aunque la fuente generadora de la anomalía tenga paredes
verticales, de ahí que la interpretación cualitativa no siempre resulta ser la apropiada, si no
se la realiza cuidadosamente y se cuenta con información geológica adicional de la zona
prospectada.
147
Geofísica Aplicada Antonio Carrasco Cortez
Goriachko-1968.
148
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo XI
Por respeto a la libertad, la ley permite que el pobre entable la lucha por la vida sin más
armas que su debilidad y su hambre.
El principio del método de carga se basa en conectar uno de los polos (electrodo de
corriente A) de la fuente de corriente eléctrica (transmisor de corriente continua),
directamente al cuerpo estudiado (afloramientos de menas producto de la erosión,
intercepción de menas por taladros, etc.) y un segundo electrodo de corriente (B) fuera de
los límites de la zona de estudio, a una distancia suficiente mente grande para que la
influencia de su campo pueda ser despreciada.
149
Geofísica Aplicada Antonio Carrasco Cortez
Se debe contar con un equipo de resistividad o uno de polarización inducida, para emplear
el transmisor según el siguiente procedimiento:
150
Geofísica Aplicada Antonio Carrasco Cortez
Los datos obtenidos manteniendo fijo el electrodo M y usando algunas estaciones base
secundarias, deben ser corregidas por enlace es decir relacionarlas a la estación base esta
corrección consiste en sumar o restar, el valor obtenido en las estaciones base secundarias
con relación a la estación principal, a las lecturas de las estaciones medidas con relación a
las estaciones secundarias, por ejemplo:
En la figura 141 las lecturas obtenidas con relación a las estaciones secundarias se
encuentran entre paréntesis en la parte inferior de los perfiles y las corregidas en la parte
superior.
151
Geofísica Aplicada Antonio Carrasco Cortez
Los datos obtenidos son referidos a un nivel base o nivel local cero, que es el potencial
promedio medido sobre sectores donde no existen perturbaciones apreciables.
REPRESENTACION DE DATOS
2000
8019820
1950
1900
8019800 1850
1800
1750
8019780 1700
1650
1600
8019760
1550
1500
8019740 1450
1400
1350
8019720 1300
1250
1200
8019700 1150
1100
1050
8019680
1000
716060 716080 716100 716120 716140 716160 716180 716200 716220 716240
152
Geofísica Aplicada Antonio Carrasco Cortez
Posteriormente se interpola los puntos de medida y se unen los puntos de igual voltaje,
como resultado se obtiene un mapa donde es posible observar la ocurrencia de anomalías,
de existir las mismas se deben de realizar cortes en la misma para evaluar la magnitud de
estas, generalmente el método es empleado para la correlación o verificar la continuidad de
estructuras mineralizadas.
Bueno al parecer la interpretación del método de carga es netamente cualitativa, pero sus
resultados son de mucha utilidad para la correlación de estructuras geológicas.
Se inyecta corriente (se debe tomar la lectura de corriente introducida por los electrodos
AB) en un cuerpo conductor (una mena mineralizada) de fuerte buzamiento, la carga
inyectada es producida por corriente continua el campo eléctrico se estudia midiendo el
gradiente de potencial a lo largo de perfiles situados en la superficie del terreno,
perpendicularmente a la dirección del cuerpo estudiado.
153
Geofísica Aplicada Antonio Carrasco Cortez
En la figura 145, Se explica el siguiente procedimiento: para la primera medida del perfil,
el electrodo M debe de estar en la posición 1 y el electrodo N en la posición 2 del perfil
entonces se toma la medida de diferencial de potencial entre uno y dos, luego se traslada
el electrodo de potencial M a la posición 2 y el electrodo de potencial N a la posición 3 y se
toma la correspondiente lectura de potencial, así sucesivamente, hasta barrer con todo el
perfil, uno debe de tener cuidado de no cambiar la polaridad de los electrodos de corriente
hasta que haiga terminado la prospección, también verifique que el electrodo M esté
conectado al borne positivo del mili voltímetro o instrumento de medida de potencial, tenga
mucho cuidado con los signos del voltaje resepcionado, anote los mismos seguido de la
medida realizada.
154
Geofísica Aplicada Antonio Carrasco Cortez
calcula el valor medio del gradiente de potencial (¥) en cada intervalo de puntos de
observación.
REPRESENTACION DE DATOS
155
Geofísica Aplicada Antonio Carrasco Cortez
Siendo “u” la distancia en metros entre el perfil dado y el que pasa por el punto de carga, el
factor “u” compensa la disminución del gradiente de potencial al alejarse del punto de
carga.
Fig-148 La distancia “u” en metros entre el perfil dado y el que pasa por el punto de carga,
método de gradiente de potencial.
156
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo XII
Y SECCIÓN ELÉCTRICA
¿Qué es la vida? Un frenesí. ¿Qué es la vida? Una ilusión, una sombra, una ficción; y el
mayor bien es pequeño; que toda la vida es sueño, y los sueños, sueños son.
Estas secciones se construyen en base a una serie de sondeos eléctricos (los mismos pueden
ejecutarse con cualquier arreglo electrodico) ejecutados sobre una línea geofísica, para su
representación se considera un sistema de coordenadas (x, y, z), en las abscisas (eje x) se
ubican los diferentes sondeos eléctricos y en las ordenadas (eje y) los valores de apertura
de separación de electrodos AB/2 donde se designa a cada punto de medida, el valor de
resistividad aparente correspondiente (el valor de resistividad aparente o verdadera según
sea el caso del tipo de sección de resistividad a construir, se estaría representando en el eje
z, como en un plano topográfico salvo el caso de que en vez de elevaciones se representan
datos de resistividad) posteriormente se configuran las líneas de ISO resistividades,
interpolando los valores de resistividad aparente de cada sondeo eléctrico, obteniéndose
como resultado la sección de eléctrica.(fig. 151)
157
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 150 Sección de resistividad aparente (Ricardo Morales, René Gonzales O. Congreso
nacional de geología Oruro -Bolivia -1979)
158
Geofísica Aplicada Antonio Carrasco Cortez
Estas secciones solo pueden dar información cualitativa, puesto que han sido construidas en
base a la interpolación de datos de resistividad aparente que es un parámetro cualitativo y a
la separación de electrodos que no tiene relación directa con la profundidad investigada,
salvo aquella que, a mayor separación de electrodos de corriente, mayor la profundidad de
investigación.
La pseudosección es una aproximación del subsuelo, ya que cada valor es afectado por la
resistividad aparente que incluye una porción de suelo sobre y debajo del punto de
medición.
Para introducir los datos de una sección de múltiples puntos de sondeo de resistividad es
necesaria la creación de un archivo de datos, el programa más usado para la creación de
este archivo de datos es el block de notas de Windows.
En la tercera línea se debe escribir los datos topográficos del sondeo respecto a su
ubicación ingresando primero la coordenada este simbolizada por la letra “x”,
posteriormente la coordenada norte simbolizada por la letra “y” finalmente la elevación
simbolizada por la letra “z”, como muestra el texto siguiente XYZ: 10 0 0
En la cuarta fila especificar los datos a ingresar del sondeo eléctrico de la siguiente manera
escriba primero “ab/2” seguido un espacio vacío” mn” y con un último espacio vacío
escriba “Rk” que será el denominador de la columna de resistividad aparente calculada,
como muestra el siguiente texto: ab/2 mn rk.
Posteriormente inserte los datos del sondeo eléctrico en forma de columnas debajo de cada
denominativo realizado en la cuarta fila, de la siguiente manera:
159
Geofísica Aplicada Antonio Carrasco Cortez
Por último, se debe de repetir el procedimiento debajo de la última columna llenada con
datos del sondeo eléctrico, desde el llenado del número de sondeo como muestra la fig. 152,
este último paso debe realizarse para el total de sondeos realizados para la construcción de
la sección de resistividad.
Al concluir el llenado de datos de los sondeos se debe cerrar el programa block de notas y
guardar el archivo con la extensión punto ves “(. ves)”.
160
Geofísica Aplicada Antonio Carrasco Cortez
Paso 1.- Abrir el archivo con extensión. ves, creado en el bloc de notas desde la pestaña
FILE y luego hacer clic en OPEN DATA.
Se desplegará una ventana llamada “abrir” en el mismo buscar el archivo creado en el block
de notas y hacer clic en abrir, como muestra la figura 153.
Fig. 153
161
Geofísica Aplicada Antonio Carrasco Cortez
Paso 2.- en la parte superior del primer grafico que es la sección de resistividad aparente se
encuentran ubicados los sondeos eléctricos realizados, simbolizados por un triángulo con
una esquina dirigida hacia abajo, si se desea ver los puntos de sondeo en vista de planta
para correlacionar con los datos del perfil de debe entrar a la pestaña “Options” y hacer clic
en Set lines / coordinates (figura 155).
Fig. 155
Una vez ubicados los puntos de sondeo eléctrico hacer click, en la barra del cursor de
selección para elegir el sondeo a interpretar al ser seleccionado el sondeo, el triángulo que
señala su ubicación se tiznara de rojo (fig. 157), alado del cursor de selección se encuentra
la opción de configuración de escala de valores para la sección de resistividad.
162
Geofísica Aplicada Antonio Carrasco Cortez
fig. 157
Paso 3.- Una vez seleccionado el sondeo eléctrico este debe de ser interpretado como se
describe en el capítulo VII, (fig. 158) todos los sondeos de la sección de resistividad, serán
interpretados de la misma forma como resultado del ajuste de curva de resistividad se ira
conformando automáticamente la sección de resistividad y espesores verdaderos (figura
159).
- Los puntos de la ventana Graphics representan los datos tomados en el campo para el
sondeo seleccionado a interpretar.
- En la ventana Model editor, las líneas en gradiente color rojo representan los parámetros
de la curva teórica.
163
Geofísica Aplicada Antonio Carrasco Cortez
- En la ventana Model editor, la curva de color rojo que se superpone a los puntos de campo
representa a la curva teórica generada por el programa.
Los valores numéricos obtenidos individualmente para cada sondeo de la ventana “Model
Table” son los datos finales de la interpretación por ordenador, los mismos corresponden a
datos de resistividad verdadera y espesores de capas.
Fig. 159 Grafico superior sección de resistividad aparente y el grafico inferior sección de
resistividad y espesores verdaderos.
La sección de resistividad aparente debe ser interpretada de forma cualitativa los cambios
de contraste en la resistividad, además de la morfología de las curvas de isoresistividad dan
a conocer la probable estructura litológica en subsuelo.
164
Geofísica Aplicada Antonio Carrasco Cortez
Existen varios programas que pueden generan secciones verticales de resistividad entro otro
de los usuales se encuentra el programa IPI 2WIN, la metodología de interpretación y
manejo de estas secciones es total mente análogo de un programa a otro.
Fig. 160 Grafico superior sección de resistividad aparente y el grafico inferior sección de
resistividad y espesores verdaderos generados con el programa IPI2Win.
Regresión de procesos gaussianos o Kriging, procede del nombre del geólogo Sudafricano
D. G. Krige, cuyos trabajos en la predicción de reservas de oro, realizados en la década del
cincuenta son considerarse como pioneros en los métodos de interpolación espacial.
Kriging encierra un conjunto de métodos de predicción espacial que se fundamentan en la
minimización del error cuadrático medio de predicción.
Teóricamente sea Z(X) sea Z1(X) el valor que toma la variable en X y sea Z (X+h) el valor
que toma la variable en el punto a una distancia vectorial “h” del punto anterior, la función
del variograma Ύ (h) para una distancia vectorial “h” viene dada por:
C= Meseta lugar donde la varianza aumenta con la distancia alcanzando un valor más o
menos constante llamado umbral o meseta.
166
Geofísica Aplicada Antonio Carrasco Cortez
Caso cuando la malla de datos es regular: Si se denomina Xi a los puntos donde se registra
una determinada magnitud física (resistividad, polarización aparente, cargabilidad, etc.) y Z
(Xi) es la magnitud de la muestra ubicada en Xi y tomando en cuenta que la distancia entre
muestras es una distancia “a” utilizando la siguiente ecuación:
Donde “k” es el que indica el múltiplo de “a” el valor de “k” varia de la unidad hasta la
longitud total de la malla en uno u otro sentido, en la cual esta forme un último par entre
puntos de medida, y “n” es el número de muestras disponibles para formar pares a una
distancia (k*a).
El siguiente ejemplo: extraído del texto de Isobel Clark, (Geoestadistica práctica) ayudara
en gran manera a comprender el proceso de construcción de un variograma en 2D con
malla regular.
Considere los valores numéricos de los vértices dela figura 162 como datos que representen
una magnitud física:
+ (30-32)² ] (2 36)
γ*(100)= 1.46(% )²
168
Geofísica Aplicada Antonio Carrasco Cortez
+ (29-32)² ] (2 33)
γ*(200) = 3.30(% )²
169
Geofísica Aplicada Antonio Carrasco Cortez
Que interviene en el cálculo Ύ (h) para un vector “h” dado es frecuentemente igual a 0 o a
1……, entonces el variograma experimental tiene un aspecto errático y resulta casi
imposible interpretarlo, para volverlo más robusto se permite alguna tolerancia de cálculo
sobre la distancia y las direcciones:
2
Ύ (h) = (1/2(N h)) ∑((Z(X) − Z(X + h)) )
Nh
170
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 166 Región de tolerancia T (h) alrededor del vector “h” (caso bidimensional).
El ancho de bandas sirve para limitar la separación del cono de tolerancia a una extensión
máxima.
171
Geofísica Aplicada Antonio Carrasco Cortez
De esta forma cada intervalo define ahora un cuadrilátero curvilíneo del grafico 167, el
cuadrilátero de vértices ABCD donde el semivariograma se calcula tomando todas aquellas
parejas de puntos cuyo vector de separación (Xi-Yi) caiga dentro del mencionado
cuadrilátero, puesto que Ύ (h) es una función de la distancia vectorial “h” y de la dirección
π π
θ, conviene en general calcular Ύ (h) en varias direcciones (por ejemplo θ = 0 , 4 , 2 , etc).
No es posible dar un criterio general para la elección adecuada del tamaño de intervalo L y
del ángulo de tolerancia Дφ se debería simplemente adoptar valores que aseguren un
número de puntos “M” suficientemente grande.
En general cuanto más grande sea L y Дφ, mas suavizado es la curva función
semibariograma.
173
Geofísica Aplicada Antonio Carrasco Cortez
174
Geofísica Aplicada Antonio Carrasco Cortez
175
Geofísica Aplicada Antonio Carrasco Cortez
-Mapas de error.
176
Geofísica Aplicada Antonio Carrasco Cortez
177
Geofísica Aplicada Antonio Carrasco Cortez
100 110 120 130 1 67 340,46 273,4601 426 0,64193 121 188,49552
100 110 130 140 2 5 89,75 84,75002 426 0,19894 150 753,98208
100 110 140 150 3 10 66,5 56,50002 426 0,13263 250 1884,9552
100 110 150 160 4 11 38,346 27,34601 426 0,06419 242 3769,9104
178
Geofísica Aplicada Antonio Carrasco Cortez
100 110 160 170 5 11 25,916 14,916 426 0,03501 231 6597,3432
100 110 170 180 6 12 22,0489 10,04893 426 0,02359 249 10555,749
100 110 180 190 7 20 26,484 6,48405 426 0,01522 241 15833,624
110 120 130 140 1 6 140,751 134,7512 200 0,67376 127 188,49552
110 120 140 150 2 3 46,7676 43,76762 200 0,21884 165 753,98208
110 120 150 160 3 6 27,2207 21,22066 200 0,1061 200 1884,9552
110 120 160 170 4 5 16,7775 11,77747 200 0,05889 222 3769,9104
110 120 170 180 5 2 8,12368 6,123677 200 0,03062 202 6597,3432
110 120 180 190 6 7 11,8125 4,812543 200 0,02406 254 10555,749
110 120 190 200 7 8 11,0568 3,056786 200 0,01528 242 15833,624
120 130 140 150 1 9 199,986 190,986 300 0,63662 120 188,49552
120 130 150 160 2 10 72,8662 62,86622 300 0,20955 158 753,98208
120 130 160 170 3 12 43,831 31,831 300 0,1061 200 1884,9552
120 130 170 180 4 15 34,019 19,01902 300 0,0634 239 3769,9104
120 130 180 190 5 16 28,3686 12,36862 300 0,04123 272 6597,3432
120 130 190 200 6 14 20,5367 6,536722 300 0,02179 230 10555,749
120 130 200 210 7 18 21,8273 3,827298 300 0,01276 202 15833,624
130 140 150 160 1 12 384,423 372,4226 360 1,03451 195 188,49552
130 140 160 170 2 13 88,9169 75,91692 360 0,21088 159 753,98208
130 140 170 180 3 11 52,062 41,06198 360 0,11406 215 1884,9552
130 140 180 190 4 10 32,9183 22,91832 360 0,06366 240 3769,9104
130 140 190 200 5 10 21,0226 11,02262 360 0,03062 202 6597,3432
140 150 160 170 1 9 454,634 445,6339 400 1,11408 210 188,49552
140 150 170 180 2 6 101,493 95,49299 400 0,23873 180 753,98208
140 150 220 230 7 2 7,05254 5,052539 400 0,01263 200 15833,624
150 160 170 180 1 2 316,066 314,0658 320 0,98146 185 188,49552
150 160 190 200 3 3 33,388 30,38799 320 0,09496 179 1884,9552
150 160 200 210 4 4 19,5335 15,53353 320 0,04854 183 3769,9104
150 160 210 220 5 5 15,671 10,67096 320 0,03335 220 6597,3432
150 160 220 230 6 6 12,9725 6,972504 320 0,02179 230 10555,749
150 160 230 240 7 7 11,8504 4,850437 320 0,01516 240 15833,624
160 170 180 190 1 10 376,693 366,6931 320 1,14592 216 188,49552
179
Geofísica Aplicada Antonio Carrasco Cortez
160 170 190 200 2 20 73,0517 53,05166 320 0,16579 125 753,98208
160 170 200 210 3 15 50,6507 35,65071 320 0,11141 210 1884,9552
160 170 210 220 4 12 33,2207 21,22066 320 0,06631 250 3769,9104
160 170 220 230 5 13 25,9022 12,90216 320 0,04032 266 6597,3432
160 170 230 240 6 20 27,4272 7,427232 320 0,02321 245 10555,749
160 170 240 250 7 12 17,2142 5,21422 320 0,01629 258 15833,624
170 180 190 200 1 22 159,934 137,9343 200 0,68967 130 188,49552
170 180 200 210 2 32 82,3991 50,39908 200 0,252 190 753,98208
170 180 210 220 3 6 31,0404 25,04038 200 0,1252 236 1884,9552
170 180 220 230 4 4 16,8385 12,8385 200 0,06419 242 3769,9104
170 180 230 240 5 5 12,0028 7,002819 200 0,03501 231 6597,3432
170 180 240 250 6 6 11,4188 5,418848 200 0,02709 286 10555,749
170 180 250 260 7 8 11,4105 3,410464 200 0,01705 270 15833,624
180 190 200 210 1 9 366,701 357,7008 465 0,76925 145 188,49552
180 190 210 220 2 2 113,011 111,0106 465 0,23873 180 753,98208
180 190 220 230 3 5 54,338 49,33804 465 0,1061 200 1884,9552
180 190 230 240 4 4 37,3032 33,30318 465 0,07162 270 3769,9104
180 190 240 250 5 12 30,9599 18,9599 465 0,04077 269 6597,3432
180 190 250 260 6 12 23,7618 11,76184 465 0,02529 267 10555,749
190 200 210 220 1 15 145,507 130,5071 205 0,63662 120 188,49552
190 200 220 230 2 14 67,8342 53,83417 205 0,26261 198 753,98208
190 200 230 240 3 12 33,7512 21,75118 205 0,1061 200 1884,9552
190 200 240 250 4 13 25,0719 12,0719 205 0,05889 222 3769,9104
190 200 250 260 5 14 20,8361 6,836085 205 0,03335 220 6597,3432
200 210 220 230 1 55 303,282 248,2818 360 0,68967 130 188,49552
200 210 230 240 2 66 166,268 100,2676 360 0,27852 210 753,98208
200 210 240 250 3 7 49,0169 42,01691 360 0,11671 220 1884,9552
200 210 250 260 4 5 25,913 20,91296 360 0,05809 219 3769,9104
210 220 230 240 1 6 185,049 179,0493 250 0,7162 135 188,49552
210 220 240 250 2 2 68,3146 66,31457 250 0,26526 200 753,98208
210 220 250 260 3 5 33,5153 28,51527 250 0,11406 215 1884,9552
180
Geofísica Aplicada Antonio Carrasco Cortez
La posición aparente de medida (p0) definida por: xd1, yd1 con el arreglo Dipolo dipolo,
transformado a sistema cartesiano es X1, Y1 que es la proyección resultante de “p0” en los
ejes de las abscisas y las ordenadas del sistema de coordenadas relativas adoptado.
Dónde: a = distancia en metros entre los electrodos AB = distancia entre electrodos MN, y
Yd2 es igual a la cota terreno, na = distancia en metros entre los electrodos BM.
progresiva Elevación
0 118
10 117
20 116
30 114
40 112
50 108
60 105
70 103
181
Geofísica Aplicada Antonio Carrasco Cortez
80 100
90 98
100 98
110 97
120 95
130 92
140 87
150 83
160 79
170 76
180 74
190 76
195 79
200 80
210 83
220 86
230 90
240 95
250 100
Tabla 2 Topografía tab1.
no X Y Z Res(OHM*m)
1 27 99 120
2 30,6 93 150
3 34 87 200
4 38 81 202
5 41 74,6 229
6 44 68 232
7 47,5 62 287
8 36,5 96,7 126
9 40 90,5 155
10 43,7 84,3 205
11 47 78 272
12 50 72 268
13 53,5 65,7 215
14 57 59,4 279
15 46 94 128
182
Geofísica Aplicada Antonio Carrasco Cortez
183
Geofísica Aplicada Antonio Carrasco Cortez
184
Geofísica Aplicada Antonio Carrasco Cortez
185
Geofísica Aplicada Antonio Carrasco Cortez
Una vez generado el archivo de datos base, ingresar a la pestaña “Grid” del programa y
elegir la opción “variogram”, hacer click en “New variogram”, seleccionar los datos base
en formato Exel de tab1g, luego en la ventana “new variograma” verificar la correlación
correcta para los datos XYZ con las posiciones en el panel de la sección eléctrica y el
respectivo valor de resistividad para Z, presione la opción ok, simultáneamente se
desplegara el variograma generado por el programa, asiendo clik en el grafico del
variograma se despliega la ventana de modelado, elegir la pestaña “experimental” donde es
posible realizar el ajuste de datos para definir la dirección, tolerancia angular, etc. del
variograma, en la pestaña “Model” se ajusta y calibra el modelo teórico de variograma para
los datos experimentales.
186
Geofísica Aplicada Antonio Carrasco Cortez
nuevamente los datos del archivo Exel de tab 1g, se desplegara la ventana Grid data, de esta
en la sub opción “Gridding method” elegir “Kriging”.
Luego hacer click en la pestaña “advanced options”, en esta nueva ventana seleccionar
“Remove” para eliminar el variograma generado por defecto por el programa, luego hacer
click en “get variograma” automáticamente asesara a los datos del modelado de variograma
teórico realizado, luego presione la opción ok, volverá a la ventana “Grid data” en esta
seleccione la opción ok, finalmente se desplegará una ventana con un reporte estadístico del
ajuste de modelo.
Para poder construir la sección eléctrica seleccionar la pestaña “Map” y elegir la opción
“New”, hay se desplegara todas las opciones de generado de isolineas en 2D y 3D, hacer
ckik en la opción que más convenga y se asesara al buscador que de forma automática le
mostrara el archivo generado por modelado y ajuste de variograma, este archivo es de tipo
common grid files o GRD surfer grid, elija la opción abrir y se desplegara simultáneamente
el panel de isocontornos en 2d o 3d o cual sea la opción de generado seleccionado.
Primer protocolo de ajuste de variograma generado con una tolerancia angular de 45 grados
en dirección X, y construcción la pseudoseccion por kriging.
187
Geofísica Aplicada Antonio Carrasco Cortez
Column C: res
Direction: 0.0 Tolerance: 45.0
8000
7000
6000
5000
Variogram
4000
3000
2000
1000
0
0 10 20 30 40 50 60 70 80
Lag Distance
Fig. 179 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “tab1g”, aplicando el modelo esférico con efecto de pepita, los parámetros
de modelo son: rango=65 y meseta=7800.
188
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 180 Pseudoseccion eléctrica de los datos “tab1g”, calculada mediante kriging usando el
variograma con tolerancia angular de 45 grados en dirección X modelo esférico con efecto
de pepita, los parámetros de modelo son: rango=65 y meseta=7800.
189
Geofísica Aplicada Antonio Carrasco Cortez
190
Geofísica Aplicada Antonio Carrasco Cortez
Column C: res
Direction: 90.0 Tolerance: 60.0
16000
14000
12000
10000
Variogram
8000
6000
4000
2000
0
0 10 20 30 40 50 60 70 80
Lag Distance
Fig. 183 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “dat1g”, aplicando el modelo esférico sin efecto de pepita, los parámetros de
modelo son: rango=68 y meseta=15000.
191
Geofísica Aplicada Antonio Carrasco Cortez
192
Geofísica Aplicada Antonio Carrasco Cortez
Se realizó el variograma con una tolerancia angular de 60 grados en dirección del eje Y, a
razón de resaltar las estructuras geológicas verticales (fallas diaclasas, etc.).
Tercer protocolo, ajustar el variograma con una tolerancia angular de 90 grados esto quiere
decir que la ventana, toma en cuenta el total de los datos en los 4 cuadrantes para el cálculo
o estimación de un punto determinado, y construya la pseudoseccion por kriging.
Fig. 186 Parámetros del variograma teórico, datos pseudoseccion “tab 1g”.
193
Geofísica Aplicada Antonio Carrasco Cortez
Column C: res
Direction: 0.0 Tolerance: 90.0
11000
10000
9000
8000
7000
Variogram
6000
5000
4000
3000
2000
1000
0
0 10 20 30 40 50 60 70 80
Lag Distance
Fig. 187 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “tab1g”, aplicando el modelo efecto de hoyo (con efecto de pepita), los
parámetros de modelo son: rango=11 y meseta =8200.
194
Geofísica Aplicada Antonio Carrasco Cortez
195
Geofísica Aplicada Antonio Carrasco Cortez
196
Geofísica Aplicada Antonio Carrasco Cortez
Capitulo XIII
TOMOGRAFÍA ELÉCTRICA
EN 2 DIMENSIONES
(Franz Tamayo)
Se estudia a lo largo de líneas en una dirección determinada, donde se inyecta una corriente
eléctrica y mediante un dispositivo receptor se observa la respuesta del subsuelo, siendo el
resultado, datos de resistividades aparentes.
Se caracteriza por ser una técnica de resistividad tetraelectrodica, cuya configuración varía
dependiendo del objetivo de estudio.
197
Geofísica Aplicada Antonio Carrasco Cortez
Con las mediciones adquiridas se construye una sección en dos dimensiones que muestra
una primera aproximación de los cambios en el subsuelo, posteriormente se aplica un
algoritmo de inversión para obtener la distribución real de resistividades o imagen eléctrica,
que será el resultado interpretable, que dará información sobre las características físicas y
geológicas de subsuelo.
Este método geofísico fue desarrollado para la investigación de áreas de geología compleja
su aplicación es diversa, como en la exploración de recursos naturales, aplicaciones
ambientales, hidrogeológicas, arqueológicas y geotécnicas.
198
Geofísica Aplicada Antonio Carrasco Cortez
Calcular la profundidad L:
L= (a/2) + (na/2)
Dónde: a = distancia en metros entre los electrodos AB = distancia entre electrodos MN,
Una vez obtenida la profundidad aparente y la posición del punto de medida se procede a
construir la pseudoseccion que se elabora interpolando los valores de los puntos de medida,
al final de dicho proceso se obtienen una imagen de resistividad aparente, para otros
arreglos electrodicos (Slumberguer, Wenner) el punto de medida es asignado al centro de
los electrodos de potencial y la profundidad aparente debe de ser calculada por los
procedimientos ya explicados en el capítulo IX del presente texto.
No existe gran variación acerca del registro de datos en el método tomografico respecto a
los otros métodos de resistividad, algunos autores de bibliografía geofísica, que no han
palpado el método a profundidad tienden a complicarse, expresando que es necesario
arreglos multielectrodicos y un sistema automático de registro de datos, úsese un sistema
automático multielectrodico para el registro de datos o úsese un sistema tetraelectrodico
manual para la adquisición de datos el resultado será el mismo yo en lo personal para cuidar
la calidad y fidelidad de los datos recomiendo el uso de un sistema tetraelectrodico de
cambio de posición en el perfil tomografico de forma manual, ya que los equipos de
registro automático tienden a cometer errores y detectar falsas anomalías en la TRE.
Fig. 191 Esquema registro de datos para TRE con arreglo Wenner.
a) Mediciones del nivel 1; b) Mediciones del nivel 2, los niveles son profundidades
aparentes de registro para la construcción de la pseudeseccion.
199
Geofísica Aplicada Antonio Carrasco Cortez
Se colocan cada uno de los electrodos a lo largo de la línea de estudio a una distancia “a”
de separación entre ellos; para realizar las medidas en el primer nivel (n=1) se utilizan
cuatro electrodos: dos de corriente y dos de potencial y las posiciones de los electrodos
serán descritas de forma correlativa en el perfil separadas por la distancia denominada “a”
en el arreglo Wenner, donde los electrodos en posición 1 y 4 son los de corriente y los
electrodos en posición 2 y 3 son los de potencial (figura 191a), para la medida del siguiente
punto al mismo nivel, se utilizan las posiciones de los electrodos 2, 3, 4 y 5 donde los
electrodos 2 y 5 son los de corriente y los electrodos 3 y 4 son los de potencial y así
sucesivamente hasta culminar las mediciones del nivel 1.
Para realizar las mediciones a lo largo de la línea de estudio en el nivel 2 (n=2), se utilizan
los electrodos en posición 1, 3, 5 y 7 donde los electrodos de corriente son 1 y 7 y los de
potencial P1 y P2 son los electrodos 3 y 5, nótese que en este caso la separación entre C1-
P1, P1-P2 y P2-C1 es “na” en este caso 2a (figura 192.a), se repite el mismo procedimiento
hasta llegar al “n” deseado, en los arreglos tomograficos se usa la nomenclatura de C1 y C2
para los electrodos de corriente “AB”, y respectivamente P1 y P2 para los electrodos de
potencial “MN”, es decir las tomografías son simplemente sondeos eléctricos sucesivos a lo
largo de un perfil.
Fig. 192 Esquema registro de datos TRE para Schlumberger a) niveles de la tomografía
tipo Schlumberger. b) Rastreo horizontal y vertical respecto de la posición de los
electrodos.
200
Geofísica Aplicada Antonio Carrasco Cortez
Para realizar una tomografía eléctrica con arreglo Dipolo-dipolo se colocan todos los
electrodos necesarios equidistantes a lo largo de la línea de estudio y se procede a realizar
las mediciones; existen dos formas de realizar las mediciones de tomografía eléctrica con
este dispositivo, una es por niveles tal y como se ejemplificó para el dispositivo Wenner y
la otra es “profundizando”.
201
Geofísica Aplicada Antonio Carrasco Cortez
202
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 195 Tomógrafos 2d con adquisición de datos por ordenador, construidos por el autor
del presente texto.
Entradas analogicas
Entradas
A DC resolucion.......................................16 bits
Rangos
203
Geofísica Aplicada Antonio Carrasco Cortez
Common-mode rejection
Device on
Device off
Anti-aliasing filter...................................None
AI+ or AI – to AGND.............................±16 V
Analog Output
Rango
204
Geofísica Aplicada Antonio Carrasco Cortez
TRANSMISOR con salida de voltaje variable en AC y DC, con una potencia de 3000 W.
puede cubrir un terreno de hasta 1 Km y 200 m de profundidad aproximadamente en una
prueba de medición.
Alimentado por una fuente de 220 VAC con capacidad de suministro de 50, 150, 250, 300,
400, 500, 700, 900, 1200 V AC/DC (en AC señal cuadrada a frecuencia variable) con una
potencia máxima de 3000 W.
ACCESORIOS
La finalidad de la inversión de los datos provistos por una Tomografía Eléctrica Resistiva
es hallar un modelo, una representación matemática idealizada, de la distribución de
resistividades de una sección del subsuelo consistente con las mediciones superficiales de
tensión sobre la superficie.
205
Geofísica Aplicada Antonio Carrasco Cortez
A estos valores de la propiedad física que se quiere estimar de las mediciones los
denominaremos parámetros del modelo.
En este caso se considerará que los valores de resistividades aparentes que se calculan a
partir de la posición de los electrodos, la corriente inyectada y la diferencia de potencial
registrada superficialmente son las mediciones.
La conexión matemática entre los parámetros del modelo y la respuesta del modelo estará
dada por el método de elementos finitos o de diferencias finitas.
En la actualidad el método de elemento finitos y el método de las diferencias finitas son los
más usados para realizar el modelado directo como subrutina dentro del procesos de
inversión pero de estos dos es el método de los elementos finitos el que se ha destacado por
ser el que mejor resuelve la inversión de los datos cuando los mismos fueron tomados sobre
una zona con desniveles topográficos, ya que en este método de cálculo las condiciones de
contorno quedan incorporadas como integrales en la funcional que se requiere minimizar
206
Geofísica Aplicada Antonio Carrasco Cortez
Del modelo idealizado se pueden obtener respuestas cuando a través del mismo se simula la
propagación de un campo de onda eléctrica, estas son datos sintéticos que pueden ser
calculados a través de relaciones matemáticas, para así definir el modelo y sus variables
que se denominaran “parámetros calculados del modelo”.
Sobre una sección real del subsuelo donde se propaga un campo electrico, se obtienen a
través del conjunto transmisor receptor (equipo resistivimentro) una serie de mediciones a
las que se llamarán “parámetros medidos”, estos parámetros físicos son calculados
mediante ecuaciones que explican el fenómeno de propagación.
Todo método de inversión procura determinar un “modelo idealizado” de subsuelo que esté
de acuerdo con los parámetros medidos y que esté sujeto a ciertas restricciones de ajuste y
calibración.
La relación matemática que determina la diferencia entre los parámetros medidos y los
calculados viene dada por la aproximación a mínimos cuadrados.
207
Geofísica Aplicada Antonio Carrasco Cortez
Donde “m” es el número de medidas. Cada dato “yi” tiene una varianza “σj”.
Estos datos surgen de aplicar una funcional sobre los parámetros del modelo “q”, que se
resuelven a través de métodos numéricos entre los cuales está el método de elementos
finitos ideal para resolver el problema directo aun cuando existen condiciones de frontera
como la topografía.
La diferencia entre los datos observados y los datos de las respuestas modeladas se da por
un vector de discrepancia “g” que se define como:
g = y-f (4)
El modelo inicial “q” es un modificado de tal forma que la suma de los errores al cuadrado
“E” de la diferencia entre la respuesta modelada “f” y los datos observados “y” sea que esta
tienda al mínimo.
208
Geofísica Aplicada Antonio Carrasco Cortez
(5)
A su vez la raíz del error cuadrático medio (RMS) donde “N” es el número total de
mediciones, también representa la bondad del ajuste para alguno de los modelos de
resistividades reales.
Para reducir el valor del error, se usa la ecuación de Gauss-Newton para determinar el
cambio en el modelo de parámetros que se reducirán a la suma de los errores al cuadrado.
(6)
(7)
Tal que Jij el cambio en el i-ésimo modelo de respuestas debido al cambio en el j-ésimo
modelo de parámetros. Una vez obtenido el vector de cambios en los parámetros es posible
obtener un nuevo modelo dado por:
(8)
209
Geofísica Aplicada Antonio Carrasco Cortez
Este proceso se repite con cada iteración, hasta que el error cuadrático varíe acercándose al
mínimo posible de la iteración anterior.
Para realizar los cálculos de resistividad aparente se necesita una malla de elementos finitos
que al menos debe incluir los límites definidos por la malla de regularización, pues debe
haber al menos un nodo en cada borde, para mayor precisión puede incluir más nodos, en
especial cerca de la superficie y en los límites inferiores.
(1)
El medio infinito se restringe a una región finita, que representa una porción de una sección
vertical de la tierra, introduciendo una frontera.
Figura 198: Elemento arbitrario P del modelo discretizado con sus cuatro puntos vecinos E,
N, W y S
El objetivo es hallar la función vi,j para los elementos con las condiciones:
2. ∂v/∂z = 0 en la superficie
211
Geofísica Aplicada Antonio Carrasco Cortez
3. qi,j = 0 excepto en los dos lugares donde se ubican los electrodos de corriente.
A continuación se obtendrá una aproximación por diferencias finitas del potencial vi,j en P
en términos de los potenciales en los elementos vecinos.
(1)
(2)
Por tanto:
(3)
(4)
212
Geofísica Aplicada Antonio Carrasco Cortez
(5)
(6)
(7)
Dónde:
(8)
Los coeficientes α dependen del punto de la grilla considerado, resolver la ecuación (7)
para cada uno de los N puntos de la grilla es equivalente a resolver el sistema.
213
Geofísica Aplicada Antonio Carrasco Cortez
Rv = b (9)
LA TOPOGRAFÍA EN LA INVERSIÓN
En adquisición de datos geoeléctricos sobre áreas con cambios significativos en la inversión
de la superficie del terreno, los efectos de la topografía deben ser tomados en cuenta cuando
se lleva a cabo una inversión del conjunto de datos adquiridos.
El uso de los factores de corrección para un modelo de subsuelo homogéneo no da
resultados lo suficientemente precisos si hay grandes variaciones de la resistividad cerca de
la superficie, en ambientes de alta resistividad (Ralt.), cubiertos de mantos de baja
resistividad (Rbaj.) ocurre una asociación geológica heterogénea eléctricamente de gran
contraste de resistividad entre la zona Rbaj. y la zona Ralt., aunque de por si las rocas de
grano fino representan un problema para los métodos eléctricos, cuando los mismos están
en la superficie y combinados con materiales altamente resistivos dejando un contraste de
relación mayor a 10:1 no es recomendable el uso de factores de corrección porque los
resultados obtenidos serán poco precisos sin embargo en ambientes que no presentan las
condiciones antes mencionadas es posible que el método arroje resultados confiables.
No solo existe la posibilidad de corregir datos de resistividad aparente que presenten
modificación por efecto de la topografía, con el método de modelado de subsuelo
homogéneo, también se puede utilizar el método de modelado de inversión con topografía,
que incorpora directamente los datos de altimetría del terreno dentro y la subrutina de
inversión, el método de inversión con topografía se ha convertido en la actualidad en el más
utilizado por estar incluido en varios software de inversión comerciales.
Para resolver el problema en presencia de topografía se utiliza en la resolución del
problema directo el método de cálculo es conocido como “Elemento Finito”.
214
Geofísica Aplicada Antonio Carrasco Cortez
Esto sucede en el caso en que las condiciones de contorno incluyan derivadas y formas
irregulares en los límites de las regiones de interés.
De acuerdo a lo expuesto en los párrafos anteriores, en el método de diferencias finitas cada
condición de contorno que incluye una derivada es aproximada por un cociente de
diferencias en los puntos de la grilla de datos en la sección de resistividad.
Esto resulta ser un inconveniente en el caso en que los contornos son irregulares, dicha
desventaja se elimina en el método de elementos finitos pues las condiciones de contorno
quedan incorporadas como integrales en la funcional que se quiere minimizar.
De esta manera el procedimiento de construcción del método resulta ser independiente de
las condiciones de borde del problema particular.
Algunos perfiles del subsuelo incluyen topografía, en este caso el problema numérico
directo se resolverá por medio de elementos finitos, en los casos en que la topografía sea
directamente plana se utilizará el método de diferencias finitas expuesto en la sección
anterior.
A continuación, se describe el método para el caso general en que se tiene una ecuación
diferencial de la siguiente forma:
(1)
(2)
Donde θ1 y θ2 son los ángulos directores de la normal a la superficie en el punto (x, y), bajo
las condiciones de que p, q, r y f sean continuas en D ∪ S, p y q tienen primeras derivadas
parciales continuas, g1 y g2 son continuas en S2 y que p (x, y) > 0, q (x, y) >0, r (x, y) ≤ 0 y
g1(x, y) > 0, se tiene que la solución a la ecuación anterior se obtiene unívocamente si se
minimiza la función:
215
Geofísica Aplicada Antonio Carrasco Cortez
(3)
φ (x, y) = a + bx + cy (4)
Fig. 200: Región debe ser dividida en un número finito de secciones, o elementos de
forma regular, tales como rectángulos o triángulos.
Por simplicidad se supondrá que D está dividida en elementos triangulares, los vértices del
triángulo se denominan nodos, el método busca una aproximación de la forma:
216
Geofísica Aplicada Antonio Carrasco Cortez
(6)
Algunas de estas constantes γn+1, γn+2, . . ., γn, son utilizadas para que se satisfaga la
condición de contorno Φ (x, y) = g (x, y) en S1.
Las constantes restantes γ1, γ2,. . ., γn son para minimizar la funcional “I”.
(7)
La condición de mínimo para la funcional “I” como función de las variables γ1, γ2, . . . , γn
es:
(8)
217
Geofísica Aplicada Antonio Carrasco Cortez
(9)
(10)
Ac = b, (11)
(12)
Para i = 1, 2, . . ., n y j = 1, 2, . . ., m y
(13)
218
Geofísica Aplicada Antonio Carrasco Cortez
(14)
(15)
Para regiones sin topografía, todas las interfaces entre las distintas regiones son verticales u
horizontales.
En el caso en que exista topografía, las coordenadas verticales varían para ajustarse al
relieve.
Ku = s (16)
219
Geofísica Aplicada Antonio Carrasco Cortez
(17)
Luego se calculan los valores de resistividad aparente para el modelo propuesto, esto se
realiza a través del método de diferencias finitas o de elementos finitos.
Para el paso inicial en el ejemplo dado, donde la tierra es homogénea, la resistividad real y
la aparente coinciden, es decir yj = m0, a partir de lo cual es sencillo calcular el vector
discrepancia g = d − y0 = d −m0.
220
Geofísica Aplicada Antonio Carrasco Cortez
221
Geofísica Aplicada Antonio Carrasco Cortez
En caso de tener topografía estos datos se agregan después de la última línea con
información de resistividad y no se colocan los 6 ceros sino hasta haber incluido dentro del
archivo de datos, los valores de topografía el orden en que se almacenan los datos se
explica a continuación:
La última línea con información de resistividad, supongamos que es la línea N.
Fila N+1 Etiqueta, se coloca 2 si se tienen datos topográficos sino se coloca 0.
Fila N+2 Esta línea corresponde al número entero K del total datos topográficos adquiridos.
Fila N+3 son dos valores, el primero es la posición x respecto del primer electrodo y el
segundo es el valor de la altura en ese punto, así se ordenan todos los valores de altimetría
en la próximas K líneas.
222
Geofísica Aplicada Antonio Carrasco Cortez
Luego de la última fila con información topográfica, se coloca una fila con un número
entero, que corresponde a la posición en los datos topográficos del primer electrodo.
Finalmente se agregan 6 ceros.
Para generar el archivo tipo dat., se utilizarán los datos de la tabla 1 “tab1” del capítulo XII,
que muestra el cálculo de resistividad aparente con arreglo Dipolo dipolo (arreglo
tomografico) y sus respectivas cotas topográficas enmarcadas en la tabla 2 (datos
topográficos tab1), el grafico 201, muestra los datos de resistividad de la planilla tab1 en
formato texto tipo .dat
85 10 2 156
95 10 2 152
105 10 2 157
115 10 2 150
125 10 2 160
135 10 2 158
145 10 2 159
155 10 2 180
165 10 2 26
175 10 2 125
185 10 2 190
195 10 2 180
205 10 2 198
215 10 2 210
225 10 2 200
10 10 3 200
20 10 3 205
30 10 3 200
40 10 3 199
50 10 3 206
60 10 3 208
70 10 3 200
80 10 3 204
90 10 3 200
100 10 3 200
110 10 3 200
120 10 3 250
130 10 3 200
140 10 3 200
150 10 3 215
160 10 3 65
170 10 3 179
180 10 3 210
190 10 3 236
200 10 3 200
210 10 3 200
220 10 3 220
230 10 3 215
15 10 4 202
25 10 4 272
35 10 4 287
45 10 4 201
55 10 4 207
65 10 4 217
75 10 4 263
85 10 4 255
95 10 4 212
105 10 4 257
115 10 4 242
224
Geofísica Aplicada Antonio Carrasco Cortez
125 10 4 222
135 10 4 239
145 10 4 240
155 10 4 64
165 10 4 183
175 10 4 250
185 10 4 242
195 10 4 270
205 10 4 222
215 10 4 219
225 10 4 220
20 10 5 229
30 10 5 268
40 10 5 246
50 10 5 213
60 10 5 236
70 10 5 229
80 10 5 276
90 10 5 265
100 10 5 219
110 10 5 269
120 10 5 231
130 10 5 202
140 10 5 16
150 10 5 220
160 10 5 266
170 10 5 231
180 10 5 269
190 10 5 220
200 10 5 220
210 10 5 220
220 10 5 220
25 10 6 232
35 10 6 215
45 10 6 233
55 10 6 219
65 10 6 286
75 10 6 263
85 10 6 289
95 10 6 233
105 10 6 245
115 10 6 231
125 10 6 254
135 10 6 230
145 10 6 55
155 10 6 44
165 10 6 238
175 10 6 245
185 10 6 286
225
Geofísica Aplicada Antonio Carrasco Cortez
195 10 6 267
205 10 6 260
215 10 6 260
30 10 7 287
40 10 7 279
50 10 7 256
60 10 7 239
70 10 7 256
80 10 7 255
90 10 7 272
100 10 7 248
110 10 7 267
120 10 7 263
130 10 7 241
140 10 7 242
150 10 7 202
160 10 7 10
170 10 7 200
180 10 7 240
190 10 7 258
200 10 7 270
210 10 7 270
2
27
0 118
10 117
20 116
30 114
40 112
50 108
60 105
70 103
80 100
90 98
100 98
110 97
120 95
130 92
140 87
150 83
160 79
170 76
180 74
190 76
195 79
200 80
210 83
220 86
230 90
226
Geofísica Aplicada Antonio Carrasco Cortez
240 95
250 100
1
0
0
0
0
0
0
Fig. 201 Datos en texto de la tomografía”tab1” en el formato .dat (datos de la tabla 1 “tab1”
del capítulo XII)
Paso 2. - Ingresar al icono de acceso al programa Zond Res 2d, (figura 202) en la ventana
principal de trabajo abrir el archivo en formato .dat, generado a partir de los datos de la
tabla” tab1” (fig. 201), desde el programa asiendo click en la pestaña file luego elija la
opción “Open file” y agá “Enter” en la opción como muestra la fig. 203.
227
Geofísica Aplicada Antonio Carrasco Cortez
Emergerán algunas ventanas que indique las características de los datos, agá aceptar en
estas ventanas, tomando atención de lo que estas muestren y corregir el archivo .dat, de ser
necesario.
Finalmente emergerá una ventana “Mesh constructor” donde se visualiza toda la
información paramétrica de los datos revisar detenidamente estos parámetros y hacer clik
en “Apply” (figura 204).
Pasó 3.- Se desplegará de forma simultánea la ventana de trabajo del programa, (figura 205)
donde:
El primer grafico muestra los datos de campo en forma de sección de resistividad en 2d.
El segundo grafico muestra la primera aproximación de cálculo de modelo de los datos de
resistividad aparente.
El tercer grafico muestra la topografía del perfil de datos y el diagrama en bloques de la
sección de resistividad.
228
Geofísica Aplicada Antonio Carrasco Cortez
Paso 4.- Hacer Ckik en la pestaña “Start inversión” (figura 206), el programa
iterara sucesivamente modelos hasta lograr el mínimo error de ajuste entre los datos de
campo y los datos calculados por el programa (fig. 207).
229
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 207 Ventanas de trabajo mostrando el modelo calculado por iteración múltiple por el
programa.
Paso 5.- Hacer Click en la opción “Layered mode” y la ventana de trabajo mostrara los
resultados de la inversión en forma de capas de resistividad tomando en cuenta en el
modelo los datos topográficos, luego hacer Click en la pestaña “Layered” y seleccionar la
opción Transparent (figura 208).
230
Geofísica Aplicada Antonio Carrasco Cortez
Paso 6.- Para finalizar hacer Click en la pestaña “Options” y elegir la opción “3D sections
plot” (figura 210), se desplegará la ventana “3D sections viewer” hacer Click en la pestaña
“3D view” (fig. 211).
231
Geofísica Aplicada Antonio Carrasco Cortez
232
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 213 Perfil tomografico del terreno (datos “tab1” del capítulo XII) (Resultado del
proceso de inversión tomando en cuenta la topografía en el programa Zond Res2D).
.
El programa Zond res2D no es gratuito, se encuentran en el buscador de internet, solo tenga
cuidado que el programa descargado, sea una versión de prueba funcional limitado
generalmente por solo el tiempo de ejecución, lo cual es una gran ventaja debido a que deja
que el programa se ejecute con todas sus funciones y todo su potencial por durante 30
minutos permitiendo el guardado del archivo hasta el momento del cierre del programa así
es posible volver a ejecutar la interpretación y continuar con la misma, se recomienda
comprobar el correcto funcionamiento del programa antes de usarlo.
Zond res2D puede encontrarse a la venta, este cuenta con una llave física con puerto USB,
si compra el programa cerciórese que esta llave este en correcto estado de funcionamiento,
no es necesario adquirir un tomógrafo para poder tener acceso a este programa puesto que
el software es de uso genérico puede trabajar con varias marcas de equipos y equipos de
fabricación local.
233
Geofísica Aplicada Antonio Carrasco Cortez
Para generar el archivo tipo dat., se utilizarán los datos de la tabla 1 “tab1” (del capítulo
XII), que muestra el cálculo de resistividad aparente con arreglo Dipolo dipolo (arreglo
tomografico) y sus respectivas cotas topográficas enmarcadas en la tabla 2 (datos
topográficos tab1).
234
Geofísica Aplicada Antonio Carrasco Cortez
3 Paso. - Se escogen los parámetros de inversión con los que se desea trabajar, para invertir
los datos, se escogió el método de inversión que incluye un suavizado del modelo de
resistividad, opción que se escoge cuando los datos son muy ruidosos para obtener un
mejor resultado.
235
Geofísica Aplicada Antonio Carrasco Cortez
Una vez que los datos hayan sido cargados al programa ingresar a la pestaña inversión y
hacer click con el ratón en la opción “Include smootihing of model resistivity”, se
desplegara una ventana emergente denominada “Smoothness Constrain”, hacer click en la
opción “yes” y luego en “ok” como muestra la figura 216.
4 Paso. - Se escogió el método para resolver el problema directo, en este caso se escogió el
método de los elementos finitos para la inversión de los datos.
236
Geofísica Aplicada Antonio Carrasco Cortez
Para elegir la opción del método de inversión de datos se debe de ingresar a la opción
“Change settings”, en este se desplegará una ventana, elegir la opción Use “Finite element
method” y hacer clik en la misma inmediatamente se desplegará una ventana emergente
denominada “Type of forward modelling method” (figura 218).
El método de inversión deberá ser elegido según las condiciones de frontera del modelo, en
términos prácticos, el método de elementos finitos se utiliza cuando las condiciones de
topografía del terreno presentan fuertes desniveles respecto de la ubicación de los
electrodos a lo largo del perfil tomografico y el método de diferencias finitas puede ser
237
Geofísica Aplicada Antonio Carrasco Cortez
usado como método de inversión satisfactorio en el caso que la topografía del terreno sea
suave y presente bajos desniveles respecto de la ubicación de los electrodos en el perfil.
Para iniciar la inversión por elementos finitos del modelo hacer click en la pestaña
inversión y se desplegará una ventana en esta elegir la opción “Least squares inversión”, se
desplegará, otra ventana la cual permite navegar en el ordenador para elegir la ubicación
del archivo de inversión, escoger a preferencia la ubicación de este archivo y hacer click en
guardar.
Durante el proceso de inversión aparece una ventana de dialogo que pregunta si se quiere
realizar otra iteración, la respuesta a la misma queda a juicio del interpretador.
En este caso se realizó un máximo de 17 iteraciones, cabe destacar que con cada iteración
se recalcula el modelo idealizado y por lo tanto los parámetros calculados con el fin de
aproximar estos últimos a los parámetros medidos.
Luego de finalizado el proceso de inversión se despliega el resultado del mismo, el cual
consiste en tres gráficos, el gráfico de los parámetros medidos, el gráfico de los parámetros
calculados y el gráfico del modelo de resistividades reales encontrado, también se puede
desplegar el gráfico de la diferencia entre los parámetros medidos y los calculados.
238
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 220 Resultado del proceso de inversión por elementos finitos en el programa
RES2DINV.
Para poder visualizar todas las ventanas de resultados de la inversión hacer click en la
opción Display emergerá una ventana con una única opción llamada “show inversión
results” hacer click en esta como muestra la figura 221.
239
Geofísica Aplicada Antonio Carrasco Cortez
Fig-223 Perfil tomografico correspondiente a los datos tab1 (Resultado del proceso de
inversión por elementos finitos tomando en cuenta la topografía en el programa
RES2DINV).
Con los resultados de la inversión se puede dar paso a la interpretación cuantitativa y
análisis de la tomografía de los datos tab1.
El resultado del análisis de la tomografía tab1, tomando en cuenta que la geología base de
la zona corresponde a deposición de material volcánico y erosión fluvio lacustre:
240
Geofísica Aplicada Antonio Carrasco Cortez
La resistividad comprendida entre 60 y 100 (ohm*m) está compuesta por arena y limos de
grano fino (zona 1 perfil tomografía tab1) su espesor promedio es de 3.5 metros.
La resistividad comprendida entre 200 y 390 (ohm*m) se encuentra compuesta por material
arenoso con clastos de grava (los mismos se encuentran concentrados principalmente en el
lecho del cauce de rio (zona 2) el espesor de esta capa es de aproximadamente 5 metros.
La resistividad comprendida entre 160 y 212 (ohm*m) está compuesta por material alterado
de la roca madre ocasionado por infiltración de agua subterránea en fracturas de
consolidación (zona 3) su espesor promedio es de 1.5 metros.
Res2dinv puede encontrarse a la venta este cuenta con una llave física de puerto USB, si
compra el programa cerciórese que esta llave este en correcto estado de funcionamiento, no
es necesario adquirir un tomógrafo para poder tener acceso a este programa puesto que el
241
Geofísica Aplicada Antonio Carrasco Cortez
software Res2DInv es de uso genérico puede trabajar con varias marcas de equipos y
también con equipos de fabricación local.
242
Geofísica Aplicada Antonio Carrasco Cortez
resultados de una tomografía ya interpretada su manejo es sencillo y casi intuitivo para los
que ya han aprendió a manejar los programas de interpretación del T.E.R.
Paso 1.- Abra el programa e ingresar a la pestaña “file”, se desplegará una cascada de
opciones hacer Click en “Open” (fig. 213).
Fig. 213
Simultáneamente se abre una ventana para asesar al directorio o discos de la PC, se debe
explorar en este el archivo creado con los datos de la tomografia, como muestra el grafico
214.
Fig. 214
243
Geofísica Aplicada Antonio Carrasco Cortez
Paso 3.- La ventana de trabajo “Table” en la pestaña “Data” muestra los datos de la
tomografia, donde N= número de datos, PsZ= profundidad teórica de registro de datos, X =
poción en el progresivo horizontal del perfil del dato de resistividad y App Res el dato de
resistividad registrado en la posición X, Psz.
244
Geofísica Aplicada Antonio Carrasco Cortez
El espacio de trabajo que muestra la pseudosección eléctrica (fig. 219) de los datos de la
tomografia, la cual se encuentra enlazada mediante interfaz gráfica con la ventana Table,
fig. 219
Fig. 220
245
Geofísica Aplicada Antonio Carrasco Cortez
Paso 5 la opción Filtering fiel data, permite realizar el tratado básico estadististico
de los datos consiguiendo con esto el filtrado por ruido eléctrico y errores de lectura (fig.
221).
246
Geofísica Aplicada Antonio Carrasco Cortez
ESTUDIO DE CASOS
EN GEOFÍSICA APLICADA
247
Geofísica Aplicada Antonio Carrasco Cortez
EVALUACION GEOFISICA
Nueva Jayoma (Labor geofísica sev 1)
Antonio Carrasco Cortez
ORURO - 2010
DESCRIPCION DEL METODO, UBICACIÓN ZONA DE ESTUDIO Y OBJETIVOS
Un sondeo eléctrico vertical (SEV) permite obtener información 1D del terreno mediante la
aplicación de un pulso de corriente como estímulo y el simultáneo registro de la diferencia
de potencial generada por el terreno a modo de respuesta. Utilizando la Ecuación de
Laplace es posible deducir el perfil estratigráfico a partir de una SEV de datos "Estímulo /
Respuesta" obtenido en terreno.
El ADD entrega información 2D del terreno y es muy utilizado en minería para detectar
zonas con alta concentración de minerales metálicos. La fase de terreno genera un mapa
"Aparente" del subsuelo, el que debe ser convertido en un mapa "Real" por intermedio de la
Ecuación de Laplace
Se cuenta que los "buenos" Geofísicos pueden deducir mentalmente el mapa real con sólo
echarle un vistazo al mapa aparen.
248
Geofísica Aplicada Antonio Carrasco Cortez
N de
sev. zona o población Este Norte elevación
Nueva Jayoma
1 771234 7979865 3200
Fig-2 Ubicación puntos sev.
OBJETIVOS
OBJETIVOS GENERALES:
OBJETIVO ESPECIFICO:
ADQUISICION DE DATOS
Procedimiento de Adquisición de Datos.
1. Escoger el punto de sondeo, es decir, se recomienda dependiendo del caso, que sea un
punto “centrado”, que exista el espacio adecuado a ambos lados de éste punto para ubicar
los electrodos tan lejos como sea necesario.
2. Ubicar el Equipo, cables, etc. en el punto central.
3. Conectar el amperímetro en serie y el voltaje en paralelo a los cables de corriente (cables
del carrete) y de potencial (cables cortos), respectivamente.
4. Con la huincha, medir una distancia AN inicial a ambos lados del centro y enterrar los
electrodos (barras de cobre) a esa distancia.
5. Con la huincha, medir una distancia AN inicial a ambos lados del centro y enterrar los
electrodos (barras de cobre) a esa distancia.
6. Tratar de que los electrodos y el centro queden alineados.
7. Conectar los cables de corriente a los electrodos de los extremos y los cables de potencial
a los electrodos centrales.
8. Cerrar el circuito de corriente en serie, uniendo un cable de corriente al amperímetro
externo y éste a su vez conectarlo a una conexión de corriente del equipo, mientras que el
otro cable de corriente se conecta directamente al equipo SEV.
9. Conectar los cables de potencial al voltímetro externo o a las conexiones de potencial del
SEV. También es posible conectar los cables a ambos voltímetros con el cuidado de que
sean en paralelo.
10. Prender el equipo SEV, el amperímetro y voltímetro.
11. Dejar el amperímetro para que mida en A y el voltímetro en mV, aunque esto dependerá
de cada experiencia.
249
Geofísica Aplicada Antonio Carrasco Cortez
250
Geofísica Aplicada Antonio Carrasco Cortez
251
Geofísica Aplicada Antonio Carrasco Cortez
IPI2Win está diseñado para la interpretación 1D de las curvas de sondeo eléctrico vertical a
lo largo de un perfil. El programa es un interpretador que resuelve problemas geológicos y
logra el ajuste de las curvas de sondeos teórica y calculada. Su enfoque a la obtención de un
resultado geológico en la característica que distingue a IPI2Win de otros programas
populares de inversión automática.
Especial atención está puesta en la interpretación iterativa flexible y cómoda para el
usuario, la solución dado por el criterio del interpretador es posible elegir, entre un juego de
soluciones equivalentes, la que mayor se ajuste tanto a los datos geofísicos (e.j. la que
produzca un menor error de ajuste) como a los datos geológicos (ej. sección geoeléctrica en
correspondencia a la información geológica). Mediante la Comparación de varias
definiciones de la estructura geológica a lo largo del perfil y no del resultado de un sondeo
aislado se determinada la mejor aproximación en IPI2Win. Esta aproximación provee la
oportunidad de utilizar la información geológica a priori, extrayendo la mayor cantidad de
datos posible en situaciones geológicas complejas.
252
Geofísica Aplicada Antonio Carrasco Cortez
253
Geofísica Aplicada Antonio Carrasco Cortez
De la gráfica:
En el caso son 4 capas siendo la capa Nª1 la más superficial, con un error de ajuste del 12.3
%
CONCLUSIONES:
254
Geofísica Aplicada Antonio Carrasco Cortez
El mismo consta de las siguientes capas sedimentarias que serán descritas desde superficie
o cota del terreno, hasta el tope del sondeo eléctrico:
La tercera capa se encuentra compuesta por limo arcillas esta capa se encuentra sub
saturada de agua, el espesor de esta capa es de 7.7 metros, debido a los valores de
permeabilidad y porosidad efectiva se estima que el caudal proporcionado por los
materiales de esta capa es bajo, por los valores de resistividad que fueron comparadas en
una tabla de resistividades y dureza relativa del agua, se concluye que el agua presenta una
salinidad baja.
255
Geofísica Aplicada Antonio Carrasco Cortez
256
Geofísica Aplicada Antonio Carrasco Cortez
EVALUACION GEOFÍSICA
HIDROGEOLOGICA
CONSTRUCCIÓN DE POZOS PARA GANADO CAMÉLIDO
MUNICIPIO BELÉN DE ANDAMARCA
Antonio Carrasco Cortez
Oruro – 2015
Es una prueba que se usa para conocer la distribución de resistividad del suelo, dadas las
propiedades eléctricas de los materiales que constituyen del subsuelo, éste manifiesta un
determinado comportamiento ante el paso de corriente eléctrica. Dicho comportamiento se
manifiesta en líneas de campo o flujo eléctrico constante para materiales homogéneos, por
ello cuando los valores y características de estos campos varían, indican un cambio de
materiales o una discontinuidad dentro de un mismo depósito. Así, los suelos gruesos o
rocosos se caracterizan por presentar una alta resistividad, mientras que las zonas arcillosas
o con altos contenidos de humedad y sales corresponden con anomalías de baja resistividad.
257
Geofísica Aplicada Antonio Carrasco Cortez
OBJETIVOS
OBJETIVOS GENERALES:
OBJETIVO ESPECIFICO:
258
Geofísica Aplicada Antonio Carrasco Cortez
259
Geofísica Aplicada Antonio Carrasco Cortez
SEV 1
260
Geofísica Aplicada Antonio Carrasco Cortez
SEV 2
SEV 3
SEV 4
261
Geofísica Aplicada Antonio Carrasco Cortez
De las gráficas:
262
Geofísica Aplicada Antonio Carrasco Cortez
INTERPRETACIÓN SEV 1
Cuadro interpretación geológica y perfil geoelectrico 1d, las capas sedimentarias son
descritas desde la cota terreno hasta el tope del sondeo eléctrico.
263
Geofísica Aplicada Antonio Carrasco Cortez
264
Geofísica Aplicada Antonio Carrasco Cortez
De las gráficas:
265
Geofísica Aplicada Antonio Carrasco Cortez
INTERPRETACIÓN SEV 2
266
Geofísica Aplicada Antonio Carrasco Cortez
267
Geofísica Aplicada Antonio Carrasco Cortez
De las gráficas:
268
Geofísica Aplicada Antonio Carrasco Cortez
INTERPRETACIÓN SEV 3
269
Geofísica Aplicada Antonio Carrasco Cortez
270
Geofísica Aplicada Antonio Carrasco Cortez
De las gráficas:
271
Geofísica Aplicada Antonio Carrasco Cortez
INTERPRETACIÓN SEV 4
272
Geofísica Aplicada Antonio Carrasco Cortez
273
Geofísica Aplicada Antonio Carrasco Cortez
CONCLUSIONES Y RECOMENDACIONES
Los Perfiles geoelectricos 1d, presentan el paleorelieve de depósitos de facie fluvio lacustre
de edad cuaternaria.
RECOMENDACIONES
Coordenadas
SALINIDAD
CAUDAL PROFUNDIDAD Recomendado
APARENTE
APARENTE PERFORACIÓN para
basada en la
(l/s) (m) perforación
Población Este Norte N sev conductividad
JACHA
PAMPA 618723 7914871 1 1.5 a 2 45m si Baja a media
QALAZAYA 615938 7917096 2 0.7 a 1 45m si Baja a media
RIO BARRAS 623154 7916070 3 2a4 48m si baja
WAIÑA
JAWIRA 626327 7920115 4 1 a 1.5 38m si baja
274
Geofísica Aplicada Antonio Carrasco Cortez
El ADD entrega información 2D del terreno y es muy utilizado en minería para detectar
zonas con alta concentración de minerales metálicos. La fase de terreno genera un mapa
aparente del subsuelo, el que debe ser convertido en un mapa "Real" por intermedio de la
Ecuación de Laplace.
275
Geofísica Aplicada Antonio Carrasco Cortez
Zona o
población Norte Este Elevación msnm.
NUEVA
Inicio perfil ESPERANZA 1
eléctrico 8205354.00 m S 295328.00 m E 160
NUEVA
Fin perfil ESPERANZA 2
eléctrico 8205394.00 m S 295441.00 m E 160
Fig-2 Ubicación labor geofísica
OBJETIVOS
OBJETIVOS GENERALES:
OBJETIVO ESPECIFICO:
276
Geofísica Aplicada Antonio Carrasco Cortez
9. Conectar los cables de potencial al voltímetro externo o a las conexiones de potencial del
SEV. También es posible conectar los cables a ambos voltímetros con el cuidado de que
sean en paralelo.
10. Prender el equipo, el amperímetro y voltímetro.
11. Dejar el amperímetro para que mida en A y el voltímetro en mV, aunque esto dependerá
de cada experiencia.
12. El equipo está listo para comenzar la medición.
14. Tomar nota de los valores de Corriente y Voltaje dados por los instrumentos. Tratar de
que los valores sean tomados al mismo tiempo.
15. Mover el switch de polaridad hacia el centro para dejar de inyectar corriente.
16. Mover el switch de polaridad en sentido contrario (hacia abajo) (Ciclo Reverso) y
anotar nuevamente los valores de Corriente y Voltaje.
17. Mover nuevamente el switch hacia el centro y apagar el equipo.
18. Desconectar los electrodos de corriente para poder moverlos.
19. Desplazar los electrodos de corriente unos metros más hacia los extremos, es decir,
variar la distancia electródica y mantener constante la separación entre los electrodos de
potencial.
20. Conectar nuevamente los electrodos de corriente y realizar una nueva medición,
repitiendo los pasos desde el 10 al 18.
21. Realizar nuevas medidas con distintos valores de AB hasta que la diferencia de
potencial entre el ciclo directo y el reverso sea pequeña (< 10 mV)
22. Aumentar la distancia MN, dejando fija la última distancia AB para poder comparar la
medición ya que se estaría midiendo el mismo punto en profundidad.
277
Geofísica Aplicada Antonio Carrasco Cortez
278
Geofísica Aplicada Antonio Carrasco Cortez
279
Geofísica Aplicada Antonio Carrasco Cortez
CONCLUSIONES:
El mismo consta de las siguientes capas sedimentarias que serán descritas desde superficie
o cota del terreno, hasta el tope del perfil eléctrico:
La segunda capa se encuentra compuesta por limo y limo arenas, esta capa se encuentra sub
saturada de agua en su base, debido a los valores de permeabilidad y porosidad efectiva se
estima que el caudal proporcionado por los materiales de esta capa es alto, por los valores
de resistividad que fueron comparadas en una tabla de resistividades y dureza relativa del
agua, se concluye que el agua presenta una salinidad baja.
281
Geofísica Aplicada Antonio Carrasco Cortez
EVALUACIÓN GEOFÍSICA
VENERO SORATA
(Provincia Larecaja rio Consata)
Antonio Carrasco Cortez
Oruro – 2018
El ADD entrega información 2D del terreno y es muy utilizado en minería y geotecnia, para
determinar las características litológicas estructurales en el subsuelo. La fase de terreno
genera un perfil "Aparente" del subsuelo, el que debe ser convertido en un perfil "Real" por
intermedio de la Ecuación de elementos y diferencias finitas combinado al método de
iteraciones sucesivas.
282
Geofísica Aplicada Antonio Carrasco Cortez
OBJETIVOS GENERALES
Se tiene como principal objetivo el de determinar por métodos indirectos en este caso
geofísicos el relieve isopaco y las características litológicas, de del sector Venero Sorata.
OBJETIVO ESPECIFICO:
ADQUISICIÓN DE DATOS
283
Geofísica Aplicada Antonio Carrasco Cortez
Se emplea un sistema simétrico cuya distancia mínima entre los electrodos BM es igual a la
separación de los dipolos AB y MN, la distancia BM (na) se incrementa proporcionalmente
a la distancia AB denominada “a”
284
Geofísica Aplicada Antonio Carrasco Cortez
El primer grafico muestra los datos de campo, el segundo muestra el modelo iterado y
ajustado por el método de elementos finitos y finalmente el tercer grafico muestra el
modelo recalculado sin tomar en cuenta la topografía.
285
Geofísica Aplicada Antonio Carrasco Cortez
PERFIL ELECTRICO 1
286
Geofísica Aplicada Antonio Carrasco Cortez
El primer grafico muestra los datos de campo, el segundo muestra el modelo iterado y
ajustado por el método de elementos finitos y finalmente el tercer grafico muestra el
modelo recalculado sin tomar en cuenta la topografía.
287
Geofísica Aplicada Antonio Carrasco Cortez
PERFIL ELECTRICO 2
288
Geofísica Aplicada Antonio Carrasco Cortez
El primer grafico muestra los datos de campo, el segundo muestra el modelo iterado y
ajustado por el método de elementos finitos y finalmente el tercer grafico muestra el
modelo recalculado sin tomar en cuenta la topografía.
289
Geofísica Aplicada Antonio Carrasco Cortez
PERFIL ELECTRICO 3
290
Geofísica Aplicada Antonio Carrasco Cortez
El primer grafico muestra los datos de campo, el segundo muestra el modelo iterado y
ajustado por el método de elementos finitos y finalmente el tercer grafico muestra el
modelo recalculado sin tomar en cuenta la topografía.
PERFIL ELECTRICO 4
292
Geofísica Aplicada Antonio Carrasco Cortez
El primer grafico muestra los datos de campo, el segundo muestra el modelo iterado y
ajustado por el método de elementos finitos y finalmente el tercer grafico muestra el
modelo recalculado sin tomar en cuenta la topografía.
293
Geofísica Aplicada Antonio Carrasco Cortez
PERFIL ELECTRICO 5
294
Geofísica Aplicada Antonio Carrasco Cortez
CONCLUSIONES:
Como principal conclusión se tienen los 5 perfiles de paleorelieve Venero Sorata, que
delimita el perfil isopaco, de la zona de interés, en este se determina la profundidad y el
espesor de las distintas capas litológicas, que componen el mismo, y su variación en
función a su posición en el corte adema que el mismo pone en manifiesto las principales
estructuras geológicas, detectadas por la tomografía eléctrica.
Perfil de paleo relieve Venero Sorata. - Este presenta el paleo relieve de deposición
sedimentaria de facie fluvial, se encuentra compuesto por las siguientes capas enumeradas
desde superficie:
Capa 2.- Sedimento cuaternario fluvial y coluvial principalmente limo arena y gravilla
dispersa
Los datos topográficos para la construcción de los 5 perfiles de paleo relieve fueron
tomados de un modelo de elevación satelital.
ESTRUCTURAS SEDIMENTARIAS
Múltiples paleocanales y lechos de erosión en la roca madre, tomando en cuenta que el
estudio se realizó para la ubicación de un yacimiento tipo venero se recomienda tomar
énfasis en la deposición sedimentaria de la capa 1 Sedimento cuaternario fluvial
principalmente arena y grava, y las zonas de erosión máxima en el lecho de roca madre,
puesto que estos pueden ser zonas de acumulación de minerales de alto peso específico.
295
Geofísica Aplicada Antonio Carrasco Cortez
296
Geofísica Aplicada Antonio Carrasco Cortez
297
Geofísica Aplicada Antonio Carrasco Cortez
298
Geofísica Aplicada Antonio Carrasco Cortez
BIBLIOGRAFÍA
- A. S. Eve y D.A. Keis. (1956). Geofísica Aplicada. España: Paraninfo.
- Daniela josefina Ravelo Matos. (2000). Tesis de grado inversión de datos e interpretación
de datos de resistividad aparente en 2D adquiridos en relieve topográfico. Venezuela:
Universidad Simón Bolívar
- Sabrina Yesica Fazzito. (2006). Tesis de grado Método de procesamiento de datos para
determinar resistividad dieléctrica del subsuelo. Argentina: Universidad de Buenos Aires.
- Xavier Emeri. (2000). Geoestadistica lineal. Chile: Claus Von Plate y Cia.
299
Geofísica Aplicada Antonio Carrasco Cortez
- Jose Rafael Lajara Vizcaino. (2007). Lab View entorno grafico de programación. España:
Alfaomega.
- Porfirio Díaz Machicao. (1966). Antología Prosa y verso de Bolivia. Bolivia: Los amigos
del libro.
300
Geofísica Aplicada Antonio Carrasco Cortez
INDICE DE FIGURAS
Imagen portada (Georgius Agricola (1556). De re Metallica 1st ed. Francia: Hieronymus
Froben & Nicolaus Episcopius) ............................................................................................ 1
Imagen contraportada (Georgius Agricola (1556). De re Metallica 1st ed. Francia:
Hieronymus Froben & Nicolaus Episcopius) ....................................................................... 2
Imagen capítulo I. (Georgius Agricola. (1556). De re Metallica 1st ed. Francia:
Hieronymus Froben & Nicolaus Episcopius) ....................................................................... 1
Fig. 1 Esquema generación de tenciones locales debidas a un cuerpo mineral. (Antonio
Carrasco Cortez. (2019)) ...................................................................................................... 7
Fig. 2 En 1912 Conrad Schlumberger, utilizando un equipo muy básico, elaboro el primer
mapa de curvas de potencial en su finca cerca de Caen, Normandía. (Geophysical
Research. 1950) .................................................................................................................... 8
Fig. 3 Primer registro eléctrico de pozo realizado por Schlumberger 1920-1929.
(Geophysical Research. 1950) ............................................................................................. 9
Imagen capítulo II. (F. Sears, M. Zemansky, H. Young. (1986). Física Universitaria.
México: Addison Wesley Iberoamericana). ..................................................................... 12
Fig. 4 Esquema ley de ohm. (Antonio Carrasco Cortez. (2019)) ...................................... 15
Fig. 5 Esquema acerca del principio de resistividad. (Antonio Carrasco Cortez. (2019)) 17
Fig. 6 Tipos de conductividad. (Antonio Carrasco Cortez. (2019)) .................................. 18
Fig. 7 Esquemas volumen de almacenamiento de agua. (Antonio Carrasco Cortez. (2019))
.............................................................................................................................................. 20
Fig. 9 Tabla de resistividades por campos de transición de valores. (D. S. Parasnis y Philips
Kearey. (1971). Geofísica Minera. España: Paraninfo) .................................................... 23
Fig. 12 Campo producido por un electrodo puntual. (Antonio Carrasco Cortez. (2019)) 26
301
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 16 Esquema equipo de resistividad según A. S. EVE y D.A. KEIS Geofísica aplicada
(1956) (A. S. EVE y D.A. KEIS Geofísica aplicada (1956)) .............................................. 31
fig. 18 Circuito transmisor con elevador de voltaje alterno y generador de baja frecuencia
(fuente de diseño propia). (Antonio Carrasco Cortez. (2019)). ............................................ 32
Fig-20 Esquema receptor con corrector de ruido (Antonio Carrasco Cortez. (2019)). ........ 35
Fig-22. Imagen resistivimetros construidos por el autor del presente texto (Antonio
Carrasco Cortez. (2019)). ..................................................................................................... 36
Fig-25. Esquema carretes de fácil construcción (Antonio Carrasco Cortez. (2019)) .......... 38
Fig. 26 Circuito transmisor y receptor con corrector de ruido (fuente de diseño propia)
(Antonio Carrasco Cortez. (2019)) ...................................................................................... 39
Imagen capítulo IV. (F. Sears, M. Zemansky, H. Young. (1986). Física Universitaria.
México: Addison Wesley Iberoamericana). ..................................................................... 41
Fig. 28 Arreglo electrodico Schlumberger (vista de perfil) (Antonio Carrasco Cortez.
(2019)). ................................................................................................................................. 42
Fig. 29. Arreglo electrodico Wenner. (Antonio Carrasco Cortez. (2019)) .......................... 42
Fig. 30 Arreglo electrodico dipolo dipolo (Antonio Carrasco Cortez. (2019)) .................... 43
Fig. 31 Arreglo electrodico medio Wenner (Antonio Carrasco Cortez. (2019)) ................. 43
302
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 35 Arreglos polares (vista de planta). (Antonio Carrasco Cortez. (2019)). .................. 45
Fig-36 Grafico de tipo de sensibilidad arreglo Wenner. (Antonio Carrasco Cortez. (2019))
.............................................................................................................................................. 46
Fig-37 Grafico de tipo de sensibilidad arreglo dipolo dipolo. (Antonio Carrasco Cortez.
(2019)). ................................................................................................................................. 47
Fig-38 Grafico de tipo de sensibilidad arreglo schlumberger. (Antonio Carrasco Cortez.
(2019)). ................................................................................................................................. 48
Fig-39 Grafico de tipo de sensibilidad arreglo polo dipolo (Antonio Carrasco Cortez.
(2019)) .................................................................................................................................. 49
Fig-40 Tabla profundidad aparente y la sensibilidad de arreglos electrodicos. (Antonio
Carrasco Cortez. (2019)) ...................................................................................................... 50
Fig. 42 esquemas principios del sev. (Antonio Carrasco Cortez. (2019)) ........................... 53
Fig. 49 Cálculo de la resistividad aparente a partir de los datos de la planilla dcx1. (Antonio
Carrasco Cortez. (2019)) ...................................................................................................... 60
Fig. 50 Puntos ploteados datos en corriente continua hoja bilogaritmica para sondeos tipo
Schlumberger. (Antonio Carrasco Cortez. (2019)) .............................................................. 62
303
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 54 tipos de curvas de resistividad de 3 capas. (Antonio Carrasco Cortez. (2019)) ....... 66
Fig-56 Esquema método de reflexión de capas. (Antonio Carrasco Cortez. (2019)) ........... 69
Fig-57 Esquema Parámetros de Darr zarrouk. (Antonio Carrasco Cortez. (2019)) ............. 70
Fig-59 comportamiento eléctrico del sub suelo (Antonio Carrasco Cortez. (2019)) ........... 74
Fig. 61 Curvas patrón para 2 capas según Tagg (Introducción a la Prospección Geofísica
(1984) – Milton B. Dobrin). ................................................................................................. 81
Fig. 65 Grafico donde solo concuasa la parte izquierda de la curva. (Antonio Carrasco
Cortez. (2019)) ..................................................................................................................... 86
Fig-66 Esquema morfología de tipos de curva de campo. (Antonio Carrasco Cortez. (2019))
.............................................................................................................................................. 87
Fig-67 Grafico auxiliar, esquema aclarativo de su uso. (Antonio Carrasco Cortez. (2019)) 88
Fig-70 Ventana de ingreso al programa ves-1.30 (Antonio Carrasco Cortez. (2019)) ........ 90
304
Geofísica Aplicada Antonio Carrasco Cortez
Fig-72 Ventana iniciando un nuevo modelo. (Antonio Carrasco Cortez. (2019)) ............... 91
Fig. 73 Ventana nuevo modelo de sev. (Antonio Carrasco Cortez. (2019)) ........................ 91
Fig. 74 Ventana tamaño de hoja de trabajo. (Antonio Carrasco Cortez. (2019)) ................. 92
Fig. 81 Ventana de creación de archivo de datos del sondeo eléctrico. (Antonio Carrasco
Cortez. (2019)) ..................................................................................................................... 97
Fig. 86 Barra de iconos del programa Ipi2win. (Antonio Carrasco Cortez. (2019)) .......... 100
Fig. 88 Método de capa simple (Antonio Carrasco Cortez. (2019)) .................................. 103
Fig-101 Icono del instalador. (Antonio Carrasco Cortez. (2019)) ..................................... 111
Fig. 102 Planilla de datos y cálculos para obtener la resistividad aparente, de un sondeo
eléctrico vertical con la configuración Sclumberguer, sin empalme de curva. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 111
Fig-103 (Antonio Carrasco Cortez. (2019)) ....................................................................... 112
Fig-104 (Antonio Carrasco Cortez. (2019)) ....................................................................... 112
Fig-105 (Antonio Carrasco Cortez. (2019)) ....................................................................... 113
Fig-106 (Antonio Carrasco Cortez. (2019)) ....................................................................... 113
Fig-107 (Antonio Carrasco Cortez. (2019)) ....................................................................... 114
Fig-108 Grafico de la ventana de configuración (Setub). (Antonio Carrasco Cortez. (2019))
............................................................................................................................................ 114
Fig. 112 Esquema de perfil de resistividad con el arreglo Wenner. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 221
Fig. 113 Esquema de perfil de resistividad con el arreglo Schlumberger. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 121
Fig-114 Calicata Dipolo dipolo a profundidad variable (Antonio Carrasco Cortez. (2019))
............................................................................................................................................ 122
Fig. 115 Esquema cateo eléctrico combinado. (Antonio Carrasco Cortez. (2019)) ........... 124
Fig. 116 Esquema del cateo eléctrico cuando los electrodos de corriente A y C están
conectados al transmisor (arreglo polo dipolo). (Antonio Carrasco Cortez. (2019)) ......... 124
Fig. 117 Esquema del cateo eléctrico combinado (Antonio Carrasco Cortez. (2019)) ..... 125
Fig. 118 Tabla empírica para el cálculo de profundidad de penetración (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 126
Fig. 119 Variación del perfil de resistividad aparente, en función a la apertura de los
electrodos de corriente según Kunetz. (Kunetz Geophisics (1950)) ................................. 127
Fig-120 Anomalía producida por un cilindro horizontal. (Antonio Carrasco Cortez. (2019))
............................................................................................................................................ 128
Fig. 121 Anomalía por contacto vertical. (Antonio Carrasco Cortez. (2019)) ................... 128
306
Geofísica Aplicada Antonio Carrasco Cortez
Fig-122 Anomalía por contacto sub vertical. (Antonio Carrasco Cortez. (2019)) ............. 129
Fig-123 Anomalía producida por una capa delgada sub vertical (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 129
Fig-124 Anomalía producida por capas delgadas verticales según Kunetz. (Kunetz
geophisics (1950)) .............................................................................................................. 130
Fig.-125 Esquema perfil teórico de resistividad. (Antonio Carrasco Cortez. (2019))........ 132
Fig. 126 Esquema polarización de electrodos en una mena. (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 134
Fig. 130 Esquema potencial inducido. (Antonio Carrasco Cortez. (2019)) ..................... 137
Fig. 131 Curva carga y descarga o tención residual. (Antonio Carrasco Cortez. (2019)) .. 138
Fig. 132 Curva de tención residual. (Antonio Carrasco Cortez. (2019))............................ 140
Fig. 133 Valor de la polarización aparente de fondo. (Antonio Carrasco Cortez. (2019)) 141
Fig. 134 Fragmento de planilla de dados de polarización aparente con medida simultanea de
resistividad aparente, las medidas se realizaron con el arreglo electrodico Dipolo dipolo.
(Antonio Carrasco Cortez. (2019)) ..................................................................................... 142
307
Geofísica Aplicada Antonio Carrasco Cortez
Fig-139 Esquema del equipo para mediciones de polarización inducida según (Goriachko
(1968) ................................................................................................................................. 148
Imagen capítulo XI. (F. Sears, M. Zemansky, H. Young. (1986). Física Universitaria.
México: Addison Wesley Iberoamericana) ..................................................................... 149
Fig-140 Esquema respecto del Método de carga o Mise a la masse (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 150
Fig. 142 Esquemas mapa de iso voltajes. (Antonio Carrasco Cortez. (2019)) ................... 152
Fig-143 Esquema corte en el mapa de iso voltajes. (Antonio Carrasco Cortez. (2019)) .. 153
Fig.- 144 Esquema correlación mediante perfiles en el método de carga. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 153
Fig. 145 Esquemas registro de datos técnica del gradiente de potencial. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 154
Fig-148 La distancia “u” en metros entre el perfil dado y el que pasa por el punto de carga,
método de gradiente de potencial. (Antonio Carrasco Cortez. (2019)) .............................. 156
Fig. 150 Sección de resistividad aparente (Ricardo Morales, René Gonzales O. Congreso
nacional de geología Oruro -Bolivia -1979) ....................................................................... 158
308
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 154 Ventana de trabajo del programa ZOND- IP -1D. (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 161
fig. 156 Ventana line Settings. (Antonio Carrasco Cortez. (2019)) ................................... 162
Fig. 159 Grafico superior sección de resistividad aparente y el grafico inferior sección de
resistividad y espesores verdaderos. (Antonio Carrasco Cortez. (2019)).......................... 164
Fig. 160 Grafico superior sección de resistividad aparente y el grafico inferior sección de
resistividad y espesores verdaderos generados con el programa IPI2Win. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 165
Fig-161 Esquema partes del variograma. (Antonio Carrasco Cortez. (2019)) ................ 166
Fig. -162 Malla regular cada nodo muestra el valor de Z (Xi) (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 167
Fig. 166 Región de tolerancia T (h) alrededor del vector “h” (caso bidimensional). (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 171
Fig. 168 Variograma experimental tolerancia angular de 90. (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 172
Fig. 169 Modelo Efecto Pepita Puro. (Antonio Carrasco Cortez. (2019)) ......................... 173
Fig. 170 Modelo Esférico (Antonio Carrasco Cortez. (2019)) ........................................... 173
309
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 171 Modelo Exponencial (Antonio Carrasco Cortez. (2019)) .................................... 174
Fig. 172 Modelo Gaussiano (Antonio Carrasco Cortez. (2019)) ..................................... 174
Fig. 174 Modelo Seno Cardinal (Antonio Carrasco Cortez. (2019)) ................................. 175
Fig. 175 Modelo Potencial (Antonio Carrasco Cortez. (2019)) ......................................... 176
Fig. 177 Dirección y tolerancia angular datos pseudoseccion “tab1g”. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 187
Fig. 179 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “tab1g”, aplicando el modelo esférico con efecto de pepita, los parámetros
de modelo son: rango=65 y meseta=7800. (Antonio Carrasco Cortez. (2019)) ................. 188
Fig. 180 Pseudoseccion eléctrica de los datos “tab1g”, calculada mediante kriging usando el
variograma con tolerancia angular de 45 grados en dirección X modelo esférico con efecto
de pepita, los parámetros de modelo son: rango=65 y meseta=7800. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 189
Fig. 181 Dirección y tolerancia angular datos pseudoseccion “tab1g”. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 190
Fig. 182 Parámetros del variograma teórico, datos pseudoseccion “tab1g”. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 190
Fig. 183 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “dat1g”, aplicando el modelo esférico sin efecto de pepita, los parámetros de
modelo son: rango=68 y meseta=15000. (Antonio Carrasco Cortez. (2019)) ................... 191
Fig. 185 Dirección y tolerancia angular datos pseudoseccion “tab1g”. (Antonio Carrasco
Cortez. (2019)) ................................................................................................................... 193
310
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 186 Parámetros del variograma teórico, datos pseudoseccion “tab 1g”. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 193
Fig. 187 Variograma teórico ajustado a los datos del variograma experimental de la
pseudoseccion “tab1g”, aplicando el modelo efecto de hoyo (con efecto de pepita), los
parámetros de modelo son: rango=11 y meseta =8200. (Antonio Carrasco Cortez. (2019))
............................................................................................................................................ 194
Fig. 191 Esquema registro de datos para TRE con arreglo Wenner.a) Mediciones del nivel
1; b) Mediciones del nivel 2, los niveles son profundidades aparentes de registro para la
construcción de la pseudeseccion. (Antonio Carrasco Cortez. (2019)) .............................. 199
Fig. 192 Esquema registro de datos TRE para Schlumberger a) niveles de la tomografía
tipo Schlumberger. b) Rastreo horizontal y vertical respecto de la posición de los
electrodos. (Antonio Carrasco Cortez. (2019)) .................................................................. 200
Fig. 193: Modelo 1 de adquisición de datos con el dispositivo Dipolo dipolo. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 201
Fig. 194: Modelo 2 de adquisición de datos con el dispositivo Dipolo dipolo. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 202
Fig. 195 Tomógrafos 2d con adquisición de datos por ordenador, construidos por el autor
del presente texto. (Antonio Carrasco Cortez. (2019)) ....................................................... 203
311
Geofísica Aplicada Antonio Carrasco Cortez
Figura 198: Elemento arbitrario P del modelo discretizado con sus cuatro puntos vecinos E,
N, W y S (Antonio Carrasco Cortez. (2019)) ..................................................................... 211
Figura 199: El modelo de regularización está definido por la malla de regularización, que
posee un valor determinado de resistividad en cada bloque, el problema directo requiere
una grilla de elementos finitos que al menos incluya los bordes de la malla de
regularización, pero que para mayor precisión se pueden incluir más nodos. (Antonio
Carrasco Cortez. (2019)) .................................................................................................... 211
Fig. 200: Región debe ser dividida en un número finito de secciones, o elementos de forma
regular, tales como rectángulos o triángulos. (Antonio Carrasco Cortez. (2019)) ............. 216
Fig. 201 Datos en texto de la tomografía”tab1” en el formato .dat (datos de la tabla 1 “tab1”
del capítulo XII) (Antonio Carrasco Cortez. (2019)) ......................................................... 227
Fig. 202 Icono de acceso al programa Zond Res 2d. (Antonio Carrasco Cortez. (2019)) . 227
Fig. 203 Opción “Open file”. (Antonio Carrasco Cortez. (2019)) ..................................... 227
Fig. 204 Ventana Mesh constructor. (Antonio Carrasco Cortez. (2019)) .......................... 228
Fig. 205 Ventana de trabajo del programa Zond res 2d. (Antonio Carrasco Cortez. (2019))
............................................................................................................................................ 229
Fig. 206 Pestaña “Start inversión”. (Antonio Carrasco Cortez. (2019)) ........................... 229
Fig. 207 Ventanas de trabajo mostrando el modelo calculado por iteración múltiple por el
programa. (Antonio Carrasco Cortez. (2019)).................................................................... 230
Fig. 209 Gráfico del modelo de resistividades reales. (Antonio Carrasco Cortez. (2019)) 231
Figura 210 Pestaña “Options” y opción “3D sections plot”. (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 231
Fig. 211 Ventana “3D sections viewer”. (Antonio Carrasco Cortez. (2019)) .................... 232
Fig. 212 Ventana Opción “3D view”. (Antonio Carrasco Cortez. (2019)) ........................ 232
Fig. 213 Perfil tomografico del terreno (datos “tab1” del capítulo XII) (Resultado del
proceso de inversión tomando en cuenta la topografía en el programa Zond Res2D).
(Antonio Carrasco Cortez. (2019)) .................................................................................... 233
Fig. 216 Ventana emergente Smoothness Constrain. (Antonio Carrasco Cortez. (2019)) . 236
Fig-217 Ventana opción método de inversión de datos en el programa RES2DINV
(Antonio Carrasco Cortez. (2019)) ..................................................................................... 237
Fig. 218 Ventana emergente Type of forward modelling method. (Antonio Carrasco Cortez.
(2019)) ................................................................................................................................ 237
312
Geofísica Aplicada Antonio Carrasco Cortez
Fig. 219 Ventana de inicio de inversión de modelo. (Antonio Carrasco Cortez. (2019)) .. 238
Fig. 220 Resultado del proceso de inversión por elementos finitos en el programa
RES2DINV. (Antonio Carrasco Cortez. (2019))................................................................ 239
Fig. 221 Ventana resultados de la inversión. (Antonio Carrasco Cortez. (2019)) ............. 239
Fig. 215 Ventana principal de trabajo X2IPI (Antonio Carrasco Cortez. (2019)) ............. 244
Fig. 216 Ventana Read. (Antonio Carrasco Cortez. (2019)) .............................................. 244
Fig. 217 Ventana Table, pestaña Data (Antonio Carrasco Cortez. (2019))........................ 244
Fig. 218 Ventana Table, pestaña Electrode (Antonio Carrasco Cortez. (2019)) ................ 245
Fig. 221 Ventana de trabajo Filtering fiel data (Antonio Carrasco Cortez. (2019))........... 246
EVALUACION GEOFISICA
Fig-1 Esquema teórico de la ejecución de un sev. (Antonio Carrasco Cortez. (2019)) ..... 247
Fig-2 Ubicación puntos sev. (Antonio Carrasco Cortez. (2019)) ....................................... 249
fig. 3- Imagen adquisición de datos (Antonio Carrasco Cortez. (2019)) ........................... 250
313
Geofísica Aplicada Antonio Carrasco Cortez
EVALUACION GEOFÍSICA
HIDROGEOLOGICA
Fig-1 Esquema teórico de la ejecución de un sev. (Antonio Carrasco Cortez. (2019)) ..... 257
Fig-2 Ubicación puntos sev. (Antonio Carrasco Cortez. (2019)) ....................................... 258
fig. 3- Imagen adquisición de datos (Antonio Carrasco Cortez. (2019)) ........................... 259
Fig-1 Esquema teórico de la ejecución de un sev. (Antonio Carrasco Cortez. (2019)) ..... 275
Fig-2 Ubicación puntos sev. (Antonio Carrasco Cortez. (2019)) ....................................... 276
fig. 3- imagen adquisición de datos (Antonio Carrasco Cortez. (2019)) ............................ 277
EVALUACIÓN GEOFÍSICA
VENERO SORATA
Fig-1 Esquema adquisición de datos. (Antonio Carrasco Cortez. (2019)) ......................... 283
fig. 2- imagen adquisición de datos (Antonio Carrasco Cortez. (2019)) ............................ 284
314
Geofísica Aplicada Antonio Carrasco Cortez
ÍNDICE DE TABLAS
315