339

Rev. Inst. investig. Fac. minas metal. cienc. geogr. vol 24 n° 48, 2021: 329 - 340 ISSN-L:1561-0888
https://doi.org/10.15381/iigeo.v24i48.21785

Relación entre la cargabilidad en el dominio del tiempo y la

cargabilidad espectral en un estudio de IP realizado por una
empresa de exploración el sector de Checa Alta, Valle de
Tambo Arequipa, 2015

Relation between time-domain chargability and spectral chargability in an IP

survey conducted by an exploration company in the high sector Checa, Valley
of Tambo Arequipa, 2015

Ricardo Fabian Peláez Mc Evoy1, Jesús Alberto Torres Guerra2

Recibido: 15/12/2020 - Aprobado: 26/10/2021 – Publicado: 23/12/2021

RESUMEN
Los estudios de geofísica orientados a la caracterización juegan un papel muy importante para proporcionar un mejor entendimiento de la estructura
interna de la tierra y específicamente de una nueva zona de estudio y como consecuencia poder mapear sus propiedades, sus características y
los efectos que generan los procesos de alteración geológica. Esto finalmente generará en una mejor interpretación geológica - geofísica en la
zona nueva zona de estudio. En este trabajo de investigación se está utilizando los datos de cargabilidad registrados en el dominio del tiempo;
primeramente, se utiliza los valores registrados a través de las 20 ventanas muestreadas en el tiempo y esta curva de decaimiento obtenida es
la que se utiliza para determinar la función de transferencia que expresa el comportamiento del subsuelo al aplicar una corriente y observar su
cargabilidad. Seguidamente se utiliza esta función de transferencia y se aplica la transformada de Fourier para descomponer esta nueva señal en
sus valores espectrales. Estos valores espectrales extraídos de la función de transferencia contienen valores de amplitud, Fase y frecuencia angular.
Con estas nuevas variables se ha generado 3 cubos adicionales y esto ha hecho posible tener más elementos y variables que ayuden en una mejor
caracterización e interpretación de la zona de estudio. Lo importante de este trabajo de investigación es la forma como se utiliza los datos registrados
en el tiempo y a través de una transformada de Fourier se cambia al dominio espectral que poseen características sensibles para cada tipo de roca o
suelo. Esto ha permitido poder entender cómo se comporta un medio al aplicarle un campo eléctrico y mediante estos resultados podemos otorgarles
un valor agregado y adicional a los datos de cargabilidad registrados en el tiempo. Todo lo señalado anteriormente se ha conseguido utilizando una
serie de algoritmos que se han desarrollado en MATLAB para finalmente obtener las variables frecuencia angular, amplitud y fase. Finalmente, estos
resultados se observan con mejor claridad con la variable amplitud donde se puede observar una diferenciación entre un tipo de roca masiva, una
roca fracturada, zonas de alteración y zonas de baja compactación. Con respecto a la variable fase se puede confirmar una diferenciación de bloques
con un buen alineamiento de la falla confirmando el rumbo y azimut de esta.
Palabras claves: Cargabilidad en el dominio del tiempo; cargabilidad espectral; transformada de Fourier.

ABSTRACT
Geophysical studies oriented to characterization play a very important role in providing a better understanding of the internal structure of the earth
and specifically of a new study area and consequently being able to map its properties, its characteristics and the effects generated by the processes
of geological alteration. This will finally generate a better geological - geophysical interpretation in the new study zone. In this research work the
chargeability data registered in the time domain is being used, firstly the values registered through the 20 windows sampled in time are used and
this decay curve obtained is the one used to determine the function transfer that expresses the behavior of the subsoil when applying a current
and observing its chargeability. This transfer function is then used, and the Fourier transform is applied to decompose this new signal into its
spectral values. These spectral values extracted from the transfer function contain values of amplitude, phase, and angular frequency. With these
new variables, 3 additional cubes have been generated and this has made it possible to have more elements and variables that help in a better
characterization and interpretation of the study area. The important thing about this research work is the way in which the data recorded over time is
used and through a Fourier transform it is changed to the spectral domain that have sensitive characteristics for each type of rock or soil. This has
allowed us to understand how a medium behaves when an electric field is applied to it and through these results, we can give an added and additional
1 Bachiller e Ingeniero Geofísico de la Universidad Nacional de San Agustín de Arequipa, Perú.
Autor para correspondencia: ricpelaez_6363@hotmail.com - ORCID: https://orcid.org/0000-0002-2724-7375
2 Universidad Nacional Mayor de San Marcos, Facultad de Ingeniería Geológica, Minera, Metalúrgica y Geográfica, Unidad de Posgrado, Lima, Perú. Docente.
E-mail: jtorresgu@unmsm.edu.pe - ORCID: https://orcid.org/0000-0001-8186-5249

© Los autores. Este artículo es publicado por la Revista del Instituto de investigación de la Facultad de minas, metalurgia y ciencias
geográficas de la Universidad Nacional Mayor de San Marcos. Este es un artículo de acceso abierto, distribuido bajo los términos de la
licencia Creative Commons Atribución 4.0 Internacional (CC BY 4.0) [https://creativecommons.org/licenses/by/4.0/deed.es] que permite
el uso, distribución y reproducción en cualquier medio, siempre que la obra original sea debidamente citada de su fuente original. Para
mas información, por favor póngase en contacto con iigeo@unmsm.edu.pe
Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

value to the chargeability data recorded over time. Everything mentioned above has been achieved using a series of algorithms that have been
developed in MATLAB to finally obtain the variables angular frequency, amplitude, and phase. Finally, these results are observed with better clarity
with the variable amplitude where a differentiation can be observed between a type of massive rock, a fractured rock, alteration zones and zones of
low compaction. Regarding the phase variable, a differentiation of blocks with a good alignment of the fault can be confirmed, confirming the bearing
and azimuth of the fault.
Keywords: Time domain chargability; spectral chargability; Fourier transform.

I. INTRODUCCIÓN relación de la anomalía de cargabilidad con algún aspecto

de mineralización que es de interés en los estudios de
El método de polarización inducida (Seigel et al., 2007) ha exploración.
madurado en los últimos 50 años desde una herramienta
que busca anomalías relativamente simples hasta utilizarlo Hay una diferencia entre los datos registrados en
como un método analítico cuantitativo. En todos los el dominio espectral aplicando una fuente de corriente
estudios realizados con el método de polarización inducida alterna donde se genera frecuencias controladas y los datos
el objetivo principal ha sido la detección y búsqueda de registrados en el tiempo de manera convencional utilizando
metales base como los sulfuros y en otros casos como
una fuente de corriente pulsante, que es nuestro caso. En
un método para determinar zonas de arcillas. En esta
este segundo caso se está encontrando una respuesta
oportunidad el trabajo de investigación está orientado a
espectral intrínseca y natural que es propia del terreno y
demostrar el valor adicional que contiene las curvas de
que está relacionada a la cargabilidad registrada.
decaimiento obtenidas en campo que reflejan los valores de
cargabilidad y mediante un método de transformación del
Mientras que las mediciones espectrales de IP
dominio del tiempo al dominio de la frecuencia, generar
generadas con una fuente alterna son aquellas que también
nuevas variables que permitan caracterizar una nueva
se registran en el dominio del tiempo, y son afectadas por
zona de estudio. Este tema de investigación nos permite
frecuencia que corresponden a efectos electromagnéticos
demostrar que las señales registradas en el tiempo contienen
información muy importante que ayudan y generan valor (EM) que podemos considerarlos como ruido y es necesario
agregado a los datos de cargabilidad. corregirlos y filtrarlos para que no altere los valores
primarios, propios de la zona en estudio.
Esta relación es muy importante porque la cargabilidad
y su respuesta espectral descompuesta en tres elementos, 1.1. Ubicación del estudio de IP
tales como la frecuencia, amplitud y fase permite ingresar El área de investigación se encuentra ubicada en el
al dominio espectral y esto abre la posibilidad de poder departamento de Arequipa, provincia de Islay. En la figura
generar varios productos que ayudaran en la interpretación 1 se muestra la ubicación del estudio y en la figura 2 se
geológica y relacionarlas con algún aspecto estructural, muestra las líneas de IP. El tiempo de viaje de la ciudad
algún cambio textural en las rocas y finalmente encontrar la de Arequipa a la zona de estudio tomo 2 horas y media a

Figura 1. mapa de ubicación de la zona de estudio.

Relación entre la cargabilidad en el dominio del tiempo y la cargabilidad espectral en un estudio de IP realizado por una empresa de
330 exploración el sector de Checa Alta, Valle de Tambo Arequipa, 2015
 Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

Figura 2. mapa de ubicación de las líneas de IP.

través de la vía de la Panamericana sur con dirección a la Seguidamente encontramos al Grupo Cabanillas de
ciudad de Moquegua. El centro de operaciones se ubicó manera discordante sobre el gneis Mesoproterozoico,
en el pueblo de Chucarapi, lugar desde donde se operó este grupo está constituido por una secuencia de
diariamente y el tiempo de viaje fue de 1 hora y media a la areniscas cuarzosas que están intercaladas con arcillitas
zona de trabajo. Las coordenadas se expresan en el sistema silicificadas del Devónico. Inmediatamente tenemos a la
UTM, Datum WGS 84 del hemisferio Sur, los vértices del Formación Chocolate del Jurásico inferior constituida por
área de estudio son las siguientes: lavas andesíticas y dacíticas con algunos horizontes de
conglomerados, sus afloramientos se ubican en Punta de
231944.7, 8114521.6 Bombón y Clemesí.

234057.4, 8114507.4 Dentro de los volcánicos eruptivos, el volcánico

Chocolate se considera uno de ellos.
234064.4, 8113798.5
La Formación Guaneros del Jurásico superior está en
231944.7, 8113805.5 discordancia al Volcánico Chocolate y está compuesta por
lavas andesíticas y dacíticas porfiríticas, aglomerados y
1.2. Contexto geológico de la zona de estudio areniscas, su origen eruptivo está ubicada más al Oeste en
El complejo basal de la costa (Bellido Bravo & Narvaez, una zona ahora cubierta por el océano Pacifico.
1960) ubicado en la costa sur del Perú, entre los poblados
Un nuevo ciclo volcánico se generó en un periodo
de San Juan de Marcona, Mollendo e Ilo forma la
entre el Cretácico superior e inicios del Terciario y esto
Cordillera de la Costa. Este tipo de rocas también son originó los extensos depósitos del Grupo Toquepala afloran
conocidas como macizo de Arequipa (Cobbing & Pitcher, desde el límite del rio Tambo hacia el Sureste, conformando
1972). El complejo basal de la costa está constituido por las zonas bajas y medias del lado occidental andino. Este
rocas del Meso-proterozoico y Paleozoico, tales como depósito muestra varias discordancias entre sus varias
gneis, esquistos, dioritas, granitos y migmatitas. Este fases, que indican interrupciones en sus ciclos, seguidos
basamento rocoso ha sufrido una deformación continua por varios procesos erosivos, que están constituidos por
y una cristalización del granito por lo que se considera brechas de lava y flujos piroclásticos que muestran una
como un ortogneis y ha sufrido inyecciones de granitos composición andesítica y dacíticas.
rojos que corresponden a una edad entre los 2000 a 600
Ma. (Cobbing E. , 1985). En la Formación Moquegua se observa la
sedimentación que se inició en el Oligoceno temprano,
INGEMMET ha realizado dataciones radiométricas encima de un terreno irregular establecido por sus rocas
en muestras de zircones determinando edades de 1862.6 volcánicas, la Formación Camaná formada de areniscas
Millones de años (+-8.1) para un gneis diorítico y de 996.4 y arcillitas de color crema y abigarradas, del Oligoceno
Millones de años (+-5.3) para rocas graníticas, asignándoles medio, la Formación Millo constituida por conglomerados
una edad que corresponde al Mesoproterozoico al complejo polimícticos con un origen diverso de las gravas
basal de la costa. Estas rocas están relacionadas a depósitos redondeadas, areniscas conglomeráticas, tobas y cenizas
de feldespatos en forma de diques. volcánica, las cuales recubren una importante área de

Peláez & Torres 331

Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

estudio y por finalmente viene la disposición de depósitos en estudio. Este diseño no experimental, se clasifica en
cuaternarios que se encuentran en las partes altas. una modalidad transversal que tienen lugar cuando se
recolectan los datos en un solo momento. El presente
Durante el Mesozoico inferior se intruyó la súper estudio se encuadra en esta investigación transversal, dado
unidad Punta Coles constituida principalmente por dioritas que se obtienen las variables espectrales y se analiza su
y granodioritas de 190 Ma. (Pitcher & Bussell, 1985); está relación después de aplicar una transformación matemática
súper unidad es el componente más antiguo del Batolito utilizando una función de transferencia como respuesta del
de la Costa, la cual intruyó a la Formación Chocolate terreno.
(Jenks, 1948), en el Cretáceo tardío y el Terciario medio se
emplazó la super unidad Ilo compuesta de diorita, diorita La metodología aplicada en el presente trabajo de
cuarcífera, granodiorita hornblendica de 103 Ma. (Pitcher investigación estuvo compuesta por las siguientes fases:
& Bussell, 1985), la cual intruyó a los intrusivos de la súper
unidad Punta Coles y en el Cretáceo Superior se emplazó • Concepción de la idea,
la súper unidad Yarabamba compuesta por intrusivos de
composición diorítica, monzonítica y granodiorítica con • Descripción del problema,
hornblenda.
• Marco teórico,
1.3. Objetivos de la investigación • Alcance de la investigación,
• Determinar la relación entre la cargabilidad en el
dominio del tiempo y su respuesta espectral en un • Definición y formulación de las hipótesis,
estudio de Polarización Inducida realizado en el
sector de Checa Alta, Valle de Tambo, Arequipa. • Diseño de la investigación,

• Determinar la relación que existe entre la • Selección de muestras,

cargabilidad y la respuesta espectral de la • Recolección de datos,
frecuencia, la fase y la amplitud en un estudio de
polarización inducida. • Análisis de los datos y
• Resultados de la investigación.
II. MÉTODOS
2.1. Polarización inducida
El método aplicado para obtener los resultados en este El método de polarización inducida consiste primero en
estudio consistió en utilizar la transformada de Fourier que inyectar corriente y generar una polarización a diferentes
es una herramienta de gran valor para analizar las funciones niveles en profundidad en el subsuelo. Hay medios en
no periódicas. Por lo tanto, la curva de decaimiento no la cual debido a la presencia de minerales tales como
considerada como una función no periódica y que ha los sulfuros generan una cargabilidad. Seguidamente se
sido registrada en el dominio del tiempo se transformó hace mediciones del decaimiento de este sobre voltaje en
al dominio de la frecuencia utilizando la transformada de el tiempo después de la cesación del pulso de corriente
Fourier mediante una función de transferencia obtenida de
(Método del dominio del tiempo) (Sheriff, 2002).
los datos registrados.
Hay dos mecanismos por los cuales el comportamiento
Esta transformada producirá valores de amplitud,
natural del medio tiene un comportamiento tipo
frecuencia y valores de ángulo de fase. Intuitivamente,
condensador. El primero se denomina polarización de
uno espera una relación uno a uno entre la cargabilidad
membrana o también conocido como electrolítico y el
en el dominio del tiempo y el ángulo de fase en el
segundo como polarización por grano o de electrodos. Este
dominio de la frecuencia, pero a pesar del hecho de que
último puede estar asociado con minerales potencialmente
esta relación es reconocida hay pocos estudios realizados
económicos, mientras que el primero no lo es.
al respecto. Se observa diferencias en la respuesta de la
polarización inducida espectral para diferentes tipos de 2.2. Polarización de membrana
rocas, y particularmente para casos de rocas mineralizadas
que exhiben una dispersión en la frecuencia de manera La polarización de membrana se genera por la impedancia
significativa (Sogade et al., 2006). de la interfase entre dos medios y es el resultado para ciertas
regiones que poseen conducción normal en una trayectoria
Este estudio de investigación tiene un carácter de poro, una capa de difusión es muy gruesa para bloquear
descriptivo ya que permite obtener las características y dicho poro, de manera que se realiza una transferencia
propiedades del fenómeno de la Cargabilidad observada. de iones de cierto tamaño y una polaridad, reduciendo el
El tipo de diseño que corresponde a este estudio es no movimiento de las cargas generando un almacenamiento
experimental, por lo tanto, no se manipula la variable que crea una “membrana” polarizada. En general este tipo
Cargabilidad, solo se registra esta variable de manera de polarización ocurre en suelos que contienen arcillas,
indirecta como una respuesta del terreno después de haber debido a que poseen una capacidad de intercambio iónico
inyectado una corriente para generar un campo eléctrico alto que bloquea parcialmente la trayectoria de la solución.
y después registrar la diferencia de potencial. En otras Por lo tanto, la polarización de membrana es un efecto de
palabras, se interpreta la realidad de una situación ya la polarización inducida que se observa comúnmente en
existente no provocada sino como una respuesta del terreno estudios de mapeos de contaminación (Sogade, 2006).

Relación entre la cargabilidad en el dominio del tiempo y la cargabilidad espectral en un estudio de IP realizado por una empresa de
332 exploración el sector de Checa Alta, Valle de Tambo Arequipa, 2015
 Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

2.3. Polarización de granos • Potencia del transmisor 5000W (Walcer

Cuando los iones fluyen a través del fluido de los poros en Geophysics Ltda.)
respuesta a un potencial aplicado, ellos pueden encontrar
granos electrónicamente conductores, por ejemplo, los • Corriente máxima: 10Amp
sulfuros. Además, ellos pueden formar una barrera con la • Dominio de Tiempo: 1, 2, 4 y 8 segundos.
cual los iones pueden interactuar eléctricamente. Para que
la corriente fluya, los iones en un lado de la barrera deben • Frecuencias de salida 0.125 Hertz
intercambiar electrones con el conductor grano, que a su
vez intercambia electrones con iones en el otro lado de la • Voltajes de salida: 200, 300, 500, 700, 1000,
barrera. El circuito eléctrico localmente consiste en una 1500, 2000, 2500 y 3000 Voltios
combinación de conducción iónica y electrónica e implica
una reacción química entre el mineral y la solución (el Este trabajo de investigación se ha divido en 4 partes.
fluido en los poros) (Dentith & Mudge, 2015).
• Determinación de la Función de transferencia
2.4. Cargabilidad
• Aplicación de la Transformada de Fourier a los
Este valor de cargabilidad se observa cuando una
datos
corriente estacionaria se inyecta en el subsuelo mediante
dos electrodos y luego esta corriente se interrumpe. Esta • Cálculo de las variables de Frecuencia-Amplitud
tensión no retornará a cero en forma instantánea, sino que y Frecuencia-Fase
decaerá lentamente en el tiempo, lo que indica que esta
carga ha sido almacenada en el medio donde se inyecto • Establecer la relación entre los valores de la
esta corriente. cargabilidad y la cargabilidad espectral.

Esta Cargabilidad Integral se define como el área bajo 2.6. Determinación del modelo exponencial
la curva de descarga donde V (t) y Vm tienen las mismas La determinación del modelo exponencial es el punto de
unidades y la Cargabilidad es expresada en milisegundos partida y tomando en consideración los datos registrados
(Telford & Geldart, L., & Sheriff, 1931). en las 20 ventanas en el dominio del tiempo se transforma
estos valores como logaritmos en base 10. Seguidamente
𝟏𝟏 𝐭𝐭𝟐𝟐 se utiliza la metodología por mínimos cuadrados para
𝐌𝐌 = ∫ 𝐕𝐕(𝐭𝐭)𝐝𝐝𝐝𝐝 (1) obtener los coeficientes a0 y a1 y determinar la función
𝐕𝐕𝐕𝐕 𝐭𝐭𝐭𝐭
exponencial que expresa el comportamiento de la curva de
2.5. Equipos utilizados decaimiento en cada punto a lo largo de cada perfil.
Instrumentación GDD INC
Una vez obtenidos los valores de los coeficientes a1 y
• Receptor digital GDD de 32 canales. aplicando el inverso del logaritmo en base 10 al valor de a2,
se obtendrá los dos coeficientes que representan la función
• Modelo GRx8-32 exponencial final que es la que representa la siguiente curva
de decaimiento. (Figura 3).
• Ventanas programables
Y = 89.183 X-0.629 (2)
• Modos disponible utilizado Cole-Cole R2 = 0.9845 (3)

Figura 3. Curva de decaimiento

Peláez & Torres 333

Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

La primera ecuación corresponde a la ecuación ecuación exponencial obtenida expresada en la variable

de potencias, mientras que la segunda corresponde al tiempo. Por lo tanto, podemos proponer lo siguiente:
coeficiente de correlación o de Pearson que presenta un
valor de 0.9845, es decir es próximo a 1. Este valor expresa 𝝏𝝏𝝏𝝏
= 𝒃𝒃 𝑿𝑿𝒂𝒂 (6)
la relación entre las dos variables. Este valor vario entre 𝝏𝝏𝝏𝝏
0 y 1, mientras este valor esté más cerca de 1 habrá una Expresando el valor de x en valor del tiempo
mayor similitud entre las dos ecuaciones. En la mayoría de
las comparaciones de los datos de las curvas este valor ha 𝝏𝝏𝝏𝝏
sido muy próximo a 1. = 𝒃𝒃 𝒕𝒕𝒂𝒂 (7)
𝝏𝝏𝝏𝝏
2.7. Regresión polinomial Aplicando integración en ambos lados de la ecuación
La Regresión Polinomial se considera un caso especial de 𝒕𝒕𝒕𝒕
la Regresión Lineal, en la cual se extiende el modelo lineal 𝑴𝑴 = 𝒃𝒃 ∫ 𝒕𝒕𝒂𝒂 𝒅𝒅𝒅𝒅 (8)
al agregar predictores adicionales que se obtienen cuando 𝒕𝒕𝒕𝒕
se elevan cada uno de los predictores originales a una Esta ecuación expresa la relación entre la cargabilidad
potencia. Este método extiende la Regresión Lineal a una y la ecuación exponencial obtenida mediante el
relación que no es lineal entre las variables dependientes procesamiento de los datos registrados en campo que
e independientes reemplazando un modelo lineal por una determina los valores de cargabilidad. Esta nueva relación
función polinomial. Por lo tanto, el método de regresión es la que define la nueva curva de decaimiento expresada
polinomial es el más adecuado para ajustar un tipo de datos en tiempo y a través de ella podemos obtener los tiempos
que exhiben un comportamiento curvilíneo. correspondientes a cada punto en profundidad para nuestra
caracterización de la zona de estudio.
Esto permite de una forma directa realizar el
modelamiento de las curvas sin tener que proponer modelos 2.10. Transformada continua de Fourier
complicados no lineales. Se realizan los ajustes lineales Una transformada de Fourier es una operación que
dentro de un espacio de mayor dimensión construido con convierte funciones en el dominio del tiempo al domino de
estas funciones básicas, este modelo tiene la facilidad de la frecuencia. Esta es la fórmula que se ha utilizado y define
adaptarse a una serie de opciones de datos de manera más el principio de esta transformada (Quinquis, 2013).
amplia.
∞
En nuestro caso se planteó obtener un polinomio de 𝓧𝓧(𝐟𝐟) = ∫ 𝐱𝐱(𝐭𝐭) 𝐞𝐞−𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣 𝐝𝐝𝐝𝐝 (9)
segundo grado o también llamado cuadrático. −∞
t : Tiempo
𝒚𝒚 = 𝒂𝒂𝟎𝟎 + 𝒂𝒂𝟏𝟏 𝒙𝒙𝒊𝒊 + 𝒂𝒂𝟐𝟐 𝒙𝒙𝟐𝟐𝒊𝒊 + 𝒆𝒆 (4)
f = Frecuencia en Hertz
2.8. Determinación de la función de transferencia
La función de transferencia es una función matemática x(t) : Funcion de Transferencia
lineal o no lineal que utiliza los modelos matemáticos
e-j2πftdt : Funcion Kernel
para caracterizar la zona de estudio y permite representar
el comportamiento dinámico y estacionario de cualquier X(f) : Espectro en funcion de la frecuencia
sistema.
Una transformada de Fourier inversa (IFT) se convierte
Una vez que se han determinados los coeficientes del del dominio de la frecuencia al dominio del tiempo. Ambas
polinomio cuadrático a través de la regresión polinomial son consistentes (Quinquis, 2013).
para cada punto del perfil se pudo definir esta función como
∞
la función de transferencia y a esta función se aplicará la
transformada de Fourier. 𝐱𝐱(𝐭𝐭) = ∫ 𝓧𝓧(𝐟𝐟) 𝐞𝐞−𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣𝐣 𝐝𝐝𝐝𝐝 (10)
−∞
𝒚𝒚 = 𝒂𝒂𝟎𝟎 + 𝒂𝒂𝟏𝟏 𝒙𝒙𝒊𝒊 + 𝒂𝒂𝟐𝟐 𝒙𝒙𝟐𝟐𝒊𝒊 (5)
III. RESULTADOS
Al aplicar la transformada de Fourier a la función de
Los estudios de Polarización Inducida se muestran
transferencia, los nuevos datos extraídos de este proceso
sobre la variable de cargabilidad. En la mayoría de los
generaran las nuevas variables de frecuencia vs Amplitud y
resultados obtenidos se obtiene una relación directa entre la
frecuencia vs Fase. concentración de sulfuros y los valores de cargabilidad. El
tipo y grado de mineralización no es fijo y se han propuesto
2.9. Relación de la cargabilidad con la ecuación
algunas tablas para mostrar estas relaciones de porcentaje
exponencial
(Telford & Geldart, L., & Sheriff, 1931).
Adicionalmente a los resultados obtenidos anteriormente,
se ha propuesto la siguiente ecuación integral, en la cual Los valores de cargabilidad encima de 20
se iguala dos funciones obtenidas. La primera es el valor milisegundos muestran valores de 0.5 % y valores mayores
de la cargabilidad M obtenida por integración de la curva a 30 milisegundos valores de 1% en sulfuros totales. No se
de decaimiento y la segunda ecuación corresponde a la ha definido una relación directa entre los valores obtenidos

Relación entre la cargabilidad en el dominio del tiempo y la cargabilidad espectral en un estudio de IP realizado por una empresa de
334 exploración el sector de Checa Alta, Valle de Tambo Arequipa, 2015
 Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

de magnetita y los valores de la cargabilidad. Esto se debe Los valores de fase dentro del cubo 3D, muestran
a que la magnetita es un mineral de origen magmático cierta diferenciación y coinciden con la anomalía de la
formado por diferenciación magmática o metamorfismo cargabilidad por lo que puede utilizarse esta variable para
de contacto. Mientras que los sulfuros como la pirita y la poder identificar zonas de mineralización de sulfuros.
calcopirita son de origen hidrotermal y metamórfico. Se
puede encontrar una asociación entre ambos, pero no es El resultado en el procesado e interpretación del
una relación directa (Aguilef Carmona, 2012). cubo 3D ha tenido como objetivo principal poder observar
en un mismo dominio ya normalizado los cambios de
Se realizó un control de calidad a los datos de resistividad, cargabilidad y las variables espectrales.
cargabilidad registrados en campo porque se detectó Finalmente, seleccionar las zonas más favorables para el
que había un acople de ruido de carácter telúrico que emplazamiento de sondeos y la identificación de zonas de
enmascaraba a las curvas de decaimiento. Por lo tanto,
alteración y zonas de fracturas y fallas.
se generó un algoritmo (Filtro de baja frecuencia) para
eliminar este ruido y limpiar la curva de decaimiento, este Estos estudios mostraron resultados positivos, lo que
proceso se aplicó antes de la inversión de los datos.
permitió inferir la presencia de un cuerpo en profundidad
En el procesamiento de datos e inversión se generó con alto grado de cargabilidad y baja resistividad. Esta
un cubo 3D con los datos de cargabilidad y resistividad. respuesta física es propia de los pórfidos profundos con alto
El objetivo prioritario fue identificar de acuerdo con la grado de mineralización dentro de un cuerpo porfirítico,
anomalía de cargabilidad algún depósito de sulfuros, lo que hace de este proyecto un buen candidato para una
asociado con una zona de alteración, este reconocimiento campaña de perforación y logeo.
fue efectuado con el fin de realizar recomendaciones en
la siguiente campaña de perforación. Este cubo 3D ha Con los datos registrados se ha generado hasta 8 cubos
permitido obtener una mejor idea de la zona de estudio y 3D con las siguientes variables que fueron obtenidas de las
poder tomar decisiones en la perforación. líneas registradas, utilizando las curvas de decaimiento.

Mediante este estudio de IP se ha determinado la • Cubo 3D con la variable Fase Normal (ver Figura 4)
existencia de una anomalía de cargabilidad en profundidad
que está asociada a una mineralización de sulfuros, • Cubo 3D con la variable Fase Logaritmo (ver Figura 5)
destacando que la zona con resistividades bajas está
relacionada a zonas de infiltración con fuerte fracturamiento • Cubo 3D con la variable Amplitud Logaritmo (ver
y alto grado de alteración que origina una disminución en Figura 6)
sus valores de resistividad, mientras que valores altos en
resistividad corresponderían a rocas de mayor dureza. • Cubo 3D con la variable Amplitud Natural (ver
Figura 7)
Los resultados de este proceso muestran que las
variables espectrales tales como la amplitud, fase y • Cubo 3D con la variable Cargabilidad (ver Figura 8)
frecuencia ofrecen un valor adicional a los datos de
cargabilidad. Las variaciones de las amplitudes generadas • Cubo 3D con la variable Resistividad Logaritmo (ver
por la transformada nos muestran que los valores de Figura 9)
amplitud definen muy claramente el límite de los dos
bloques por donde pasa una falla de carácter regional. • Cubo 3D con la variable Cargabilidad filtrada (ver
Además, se definen dentro del primer bloque tres pequeños Figura 10)
bloques y que coinciden con dos fracturas transversales a
la falla principal. • Cubo 3D con la variable Tiempo (ver Figura 11)

Figura 4. Cubo 3D con variable Fase Normal

Peláez & Torres 335

Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

Figura 5. Cubo 3D con variable Fase Logaritmo

Figura 6. Cubo 3D con variable Amplitud

Figura 7. Cubo 3D con variable Amplitud Logaritmo

Figura 8. Cubo 3D con variable Cargabilidad

Relación entre la cargabilidad en el dominio del tiempo y la cargabilidad espectral en un estudio de IP realizado por una empresa de
336 exploración el sector de Checa Alta, Valle de Tambo Arequipa, 2015
 Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

Figura 9. Cubo 3D con variable Fase Logaritmo

Figura 10. Cubo 3D con variable Amplitud

Figura 11. Cubo 3D con variable Amplitud Logaritmo

IV. DISCUSIÓN variable espectral. Esto se confirmó con la variable fase

y esto ayudará en un futuro utilizar este método para
A partir de los cubos 3D obtenidos, aceptamos la hipótesis
identificar y caracterizar zonas exploratorias nuevas.
general y especifica que está orientada a encontrar una
relación entre el valor de la cargabilidad y su valor en el
Estos resultados guardan relación con los que sostienen
dominio espectral, en este caso un valor de las variables de
(Martin et al., 2020) donde recientes investigaciones han
amplitud, fase y frecuencia. Mediante estos resultados se
pudieron determinar e identificar zonas de mineralización demostrado que mediciones en el dominio de la frecuencia
en profundidad en la zona de estudio. De los resultados (FDIP), tienen una dependencia de la frecuencia con la
obtenidos, nuestro objetivo fue correlacionar la anomalía conductividad/resistividad compleja a través de mediciones
de cargabilidad encontrada en profundidad con alguna recopiladas en el dominio del tiempo (TDIP).

Peláez & Torres 337

Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

Al comparar la similitud entre las mediciones en • Después de aplicar el filtro digital a las curvas de
el dominio de tiempo, frecuencia y fase como se usa decaimiento, estas mostraron un comportamiento
en el método de polarización inducida para estudios normal sin incluir valores negativos a estas
de prospección geofísica, primero se deben definir los curvas.
parámetros y establecer la base para la comparación. Las
condiciones para la equivalencia matemática siguen luego • Se obtuvo la función de transferencia utilizando
en una secuencia natural. Las mediciones de laboratorio se los datos registrados y a través de la cual se ha
realizan en una variedad de muestras de rocas para indicar podido extraer el contenido espectral a los datos
el tipo de comparación y correlación que se espera de las de cargabilidad.
mediciones de campo (ZONGE et al., 1972).
• Las anomalías de resistividad están controladas
Los resultados indican que las mediciones de fase y por fallas que presentan una pendiente de 60 y
dominio de frecuencia son equivalentes, pero se necesitan 80 grados aproximadamente. Estos eventos son
dos mediciones de dominio de frecuencia para producir de moderados a débiles y estarían relacionados
la misma cantidad de información de polarización que a procesos de fallamiento, así como con la
una medición monofásica. Las mediciones en el dominio presencia de cuerpos resistivos, fracturados y con
del tiempo y la frecuencia serán similares, pero nunca control estructural, como se muestra en la sección
equivalentes debido a las diferencias básicas en sus de resistividad.
respectivas definiciones de parámetros (ZONGE et al.,
1972). Esto está acorde con los que se ha encontrado en • Las anomalías de resistividad mostraron valores
este estudio. moderados a bajos, exponiendo mayores valores
a mayores profundidades.

V. CONCLUSIONES • En el cubo que se ha generado, el contenido de

sulfuros totales es un factor de primer orden
• De los datos de cargabilidad registrados en el que explica las variaciones en la intensidad del
dominio del tiempo se ha extraídos datos que efecto de la polarización inducida (cargabilidad
corresponden al dominio espectral y cuyas eléctrica), obteniendo que la relación entre la
variables son la frecuencia angular, la amplitud concentración de sulfuros y su cargabilidad
y la fase eléctrica presenta en general una proporcionalidad
directa, pero no muestra una relación lineal. Por
• Con estas nuevas variables hemos generado 3
cubos adicionales y esto ha hecho posible mostrar lo tanto, se puede establecer, que en este caso
una mejor caracterización de la zona de estudio y los valores de cargabilidad fueron de un valor
estos resultados se ha utilizado en la parte de la moderado a bajo.
interpretación geofísica.

• El uso de una mono frecuencia de 0.125 Hertz VI. AGRADECIMIENTOS

utilizada por el transmisor en la inyección de Agradecimientos a la empresa SGA GEOFISICA por
corriente a lo largo del registro de los datos de IP haber proporcionado los datos del estudio de polarización
no ha limitado la obtención de nuevas frecuencias inducida. Asimismo, a todas aquellas personas que han
angulares que reflejan la respuesta del terreno en apoyado en la revisión, evaluación y chequeo de este
superficie como en profundidad. en superficie
trabajo de investigación.
como en profundidad.

• Mediante este estudio fue posible confirmar la

presencia de la falla regional Tambo- Locumba, VII. REFERENCIAS
la cual tiene una dirección Noroeste - Sur Este. Aguilef Carmona, S. G. (2012). Analisis geológico de estudios
geofísicos de Polarización Inducida (IP) en sistemas IOCG
• Se obtuvo información sobre las propiedades y pórfido cuprifero del Norte de Chile. In Universidad de
físicas de la zona de estudio, en base a lo cual, Chile. Universidad de Chile.
se sustentó el desarrollo de un modelo geológico-
geofísico que facilitará la interpretación de este Bellido Bravo, E., & Narvaez, S. (1960). Geología del
estudio en la siguiente fase. cuadrángulo de Atico (Hoja 33-o)-[Boletín A 2]. Repositorio
Institucional INGEMMET, 2. https://repositorio.ingemmet.
• Los valores de cargabilidad obtenidos en el cubo gob.pe/handle/20.500.12544/32
3D varían entre 6 y 12 milivoltios/voltio y esto
Dentith, M. C., & Mudge, S. (2015). Geophysics for the mineral
correspondería a un 0.4 % de sulfuros, lo cual exploration geoscientist. Choice Reviews Online, 52(06),
coincide con una zona en la que se presenta un 52-3156-52–3156. https://doi.org/10.5860/choice.186586
dique en profundidad.
Martin, T., Günther, T., Orozco, A. F., & Dahlin, T. (2020).
• Se realizo un buen control de calidad a los datos Evaluation of spectral induced polarization field
registrados y se aplicó un filtro digital a los datos measurements in time and frequency domain. Journal
para eliminar el ruido telúrico que se metió en la of Applied Geophysics, 180. https://doi.org/10.1016/j.
señal primaria y la distorsiono. jappgeo.2020.104141

Relación entre la cargabilidad en el dominio del tiempo y la cargabilidad espectral en un estudio de IP realizado por una empresa de
338 exploración el sector de Checa Alta, Valle de Tambo Arequipa, 2015
 Rev. Inst. investig. Fac. minas metal. cienc. geogr. 24(48) 2021

Quinquis, A. (2013). Digital signal processing using Telford, W., & Geldart, L., & Sheriff, R. (1931). Applied
MATLAB®. In Digital Signal Processing Using geophysics. In Nature (Vol. 127, Issue 3212, pp. 783–785).
MATLAB® (3rd ed.). Global Engineering: https://doi. https://doi.org/10.1038/127783a0
org/10.1002/9780470610992
ZONGE, K. L., SAUCK, W. A., & SUMNER, J. S. (1972).
Seigel, H., Nabighian, M., Parasnis, D. S., & Vozoff, K. (2007). COMPARISON OF TIME, FREQUENCY, AND PHASE
The early history of the induced polarization method. MEASUREMENTS IN INDUCED POLARIZATION.
Leading Edge (Tulsa, OK), 26(3), 312–321. https://doi. Geophysical Prospecting, 20(3), 626–648. https://doi.
org/10.1190/1.2715054 org/10.1111/j.1365-2478.1972.tb00658.x

Sheriff, R. E. (2002). Encyclopedic Dictionary of Applied

Geophysics. In Encyclopedic Dictionary of Applied
Geophysics. https://doi.org/10.1190/1.9781560802969

Sogade, J. A., Scira-Scappuzzo, F., Vichabian, Y., Shi, W., Rodi,

W., Lesmes, D. P., & Morgan, F. D. (2006). Induced-
polarization detection and mapping of contaminant plumes.
Geophysics, 71(3), 71. https://doi.org/10.1190/1.2196873

Peláez & Torres 339