INSTITUTO POTOSINO DE INVESTIGACIÓN

CIENTÍFICA Y TECNOLÓGICA, A.C.

POSGRADO EN GEOCIENCIAS APLICADAS

Estudio de la estructura somera en el municipio
de Cerro de San Pedro por medio del uso de
“Título degeofísicos
métodos la tesis”

Tesis que presenta

Jorge Luis Flores Salgado

Para obtener el grado de
Maestro en Geociencias Aplicadas

Director de la Tesis:

Dr. Vsevolod Yutsis

San Luis Potosí, S.L.P., Julio de

2019

i
 Constancia de aprobación de la tesis

La tesis "Estudio de la estructura somera en el municipio de Cerro de San

Pedro por medio del uso de métodos geofísicos" presentada para obtener el
Grado de Maestro en Geociencias Aplicadas, fue elaborada por Jorge Luis Flores
Salgado y aprobada el tres de julio del dos mil diecinueve por los suscritos,
designados por el Colegio de Profesores de la División de Geociencias Aplicadas
del Instituto Potosino de Investigación Científica y Tecnológica, A.C.

Dr. Vsevolod Yutsis

Director de la tesis

Dr. José Alfredo Ramos Leal

Miembro del Comité Tutoral

.·~

ii
 Créditos Institucionales

Esta tesis fue elaborada en el Laboratorio de Geofísica de la División de Geociencias

Aplicadas del Instituto Potosino de Investigación Científica y Tecnológica, A.C., bajo la
dirección del [Link] Yutsis.

Durante la realización del trabajo el autor recibió una beca académica del Consejo
Nacional de Ciencia y Tecnología N° Registro 590573 y del Instituto Potosino de
Investigación Científica y Tecnológica, A. C.

ii
 Instituto Potosino de Investigación
Científica y Tecnológica/ A.C.
Acta de Examen de
1PICYT
El Secretario Académico del Instituto Potosino de Investigación Científica y Tecnológica, A.C.,
certifica que en el Acta 039 del Libro Primero de Actas de Exámenes de Grado del Programa de
Maestría en Geociencias Aplicadas está asentado lo siguiente:

En la ciudad de San Luis Potosí a los 3 días del mes de julio del año 2019, se reunió a las 10:00
horas en las instalaciones del Instituto Potosino de Investigación Científica y Tecnológica, A.C., el
Jurado integrado por:

Dr. Vsevolod Yutsis Presidente IPICYT

Dr. José Alfredo Ramos Leal Secretario IPICYT
Dr. Ornar Delgado Rodríguez Sinodal IPICYT
a fin de efectuar el examen, que para obtener el Grado de:

MAESTRO EN GEOCIENCIAS APLICADAS

sustentó el C.
Jorge Luis Flores Salgado
sobre la Tesis intitulada:

Estudio de la estructura somera en el municipio de Cerro de San Pedro por medio del uso de
métodos geoflsicos
que se desarrolló bajo la dirección de

Dr. Vsevolod Yutsis

El Jurado, después de deliberar, determinó

APROBARLO

Dándose por terminado el acto a las 11: 15 horas, procediendo a la firma del Acta los integrantes
del Jurado. Dando fe el Secretario Académico del Instituto.
A petición del interesado y para los fines que al mismo convengan, se extiende el presente
documento en la ciudad de San Luis Potosí, S.L.P., México, a los 3 días del mes de julio de 2019.

o~,~
Secretario Académic~--==----===~~.,.
,,,.~~;:s.;,
 Agradecimientos
Al Instituto de Investigación Científica y Tecnológica (IPICyT) por todo el apoyo
refrendando a lo largo de este tiempo.
Al Consejo Nacional de Ciencia y Tecnología (CONACyT) por la beca y facilidades
otorgadas durante este periodo.
A mi director de tesis el Dr. Yutsis, por toda la paciencia que mostró y nunca me dejó solo
en este proceso, sin su apoyo no habría podido hacer nada. “¡Último esfuerzo Doc!”
A mis sinodales el Doc. José Alfredo que siempre veía lo positivo en trabajo mostrado,
también al Doc. Omar por siempre buscar hasta el más mínimo detalles. ¡Gracias!
A mis padres que siempre estuvieron al pendiente de mis estudios, ellos llevan el crédito
de cada pequeño logro que puedo culminar. Sin ustedes no sería nada en esta vida.
A mi ahora esposa Rosario (Chayo) que estuvo ahí día con día para apoyarme
incondicionalmente a lo largo de todo este proceso, que la veces que creí claudicar ahí
estuvo para hacerme fuerte. Gracias por todos los momentos que compartimos, sin duda
alguna una de las mejores de nuestras vidas.
A mis fieles compañeros: Daniel Velásquez, César Muñoz, Jesús Mota, Karla Pérez, en
general a toda la pecera. Esas idas a campo créanme que jamás las voy a olvidar, toda
una travesía ir hacer levantamientos geofísicos en SLP.
A la Familia Islas Cruz por todo el apoyo brindado, por acogerme y hacerme parte de su
familia.
A todos mis amigos de FCT que llevo en el corazón (Rafael Guerra, Diana Cruz, Jorge
Castro … Y toda la Familia Sable)
Por último, quiero agradecer a Dios por siempre darme fuerzas para poder continuar en
cualquier circunstancia.

iv
Contenido

Constancia de aprobación de la tesis ............................................................................................... i

Créditos Institucionales ...................................................................................................................ii
Agradecimientos ............................................................................................................................ iv
Contenido .......................................................................................................................................v
Lista de Figuras .............................................................................................................................. ix
Lista de Tablas ............................................................................................................................... xi
Resumen....................................................................................................................................... xii
Abstract ....................................................................................................................................... xiii
1 Generalidades ............................................................................................................................. 1
1.1 Introducción ......................................................................................................................... 1
1.2 Objetivos .............................................................................................................................. 2
1.2.1 General.......................................................................................................................... 2
1.2.2 Específicos ..................................................................................................................... 2
1.3 Delimitación del problema ................................................................................................... 2
1.4 Hipótesis .............................................................................................................................. 3
1.5 Metodología......................................................................................................................... 3
1.6 Localización .......................................................................................................................... 4
1.7 Ubicación fisiográfica ........................................................................................................... 5
1.8 Antecedentes ....................................................................................................................... 6
2 Marco Geológico ......................................................................................................................... 8
2.1 Geología Regional................................................................................................................. 8
2.1.2 Plataforma Valles San Luis Potosí ....................................................................................... 9
2.1.3. Formación la Peña .......................................................................................................... 11
2.1.4. Formación Cuesta del Cura ............................................................................................. 11
2.1.5. Formación Indidura ........................................................................................................ 11
2.1.6. Formación Cárdenas ....................................................................................................... 11

v
 2.1.7. Formación Casita Blanca ................................................................................................. 11
2.1.8. Ignimbrita Santa María ................................................................................................... 12
2.2. Geología y estratigrafía local ............................................................................................. 12
3 Marco Teórico ........................................................................................................................... 14
3.1 Métodos de Prospección Geofísica una breve introducción. ............................................... 14
3.2 Método magnético ............................................................................................................. 14
3.2.1 El campo magnético terrestre ...................................................................................... 15
3.2.2 Componentes vectoriales del campo magnético .......................................................... 16
3.2..3 Fuentes principales del campo magnético terrestre .................................................... 19
3.2.4 Variaciones del campo geomagnético .......................................................................... 21
3.3 Método gravimétrico.......................................................................................................... 22
3.3.1. Introducción ............................................................................................................... 22
3.3.2 Medición de la Gravedad ............................................................................................. 24
3.3.3 El geoide ...................................................................................................................... 24
3.3.4 Densidad de las Rocas .................................................................................................. 26
3.3.5 Método de Nettleton ................................................................................................... 28
3.3.6 Reducción de datos Gravimétricos ............................................................................... 29
3.3.7 Corrección por deriva .................................................................................................. 30
3.3.8 Corrección por mareas................................................................................................. 31
3.3.9 Corrección por latitud .................................................................................................. 31
3.3.10 Corrección por aire libre ............................................................................................ 33
3.3.11 Corrección de Bouguer .............................................................................................. 34
3.3.12 Corrección por terreno (topográfica) ......................................................................... 35
3..3.13Corrección isostática.................................................................................................. 36
3.3.14 Anomalía de Aire Libre .................................................................................................. 37
3.3.15 Anomalía de Bouguer ................................................................................................ 37
3.4 Método electromagnético .................................................................................................. 38
3.4.1 GPR ............................................................................................................................. 38
3.4.2 Fundamentos Teóricos del GPR ................................................................................... 39
3.4.3 Ecuaciones de Maxwell ................................................................................................ 39
3.4.4 Ecuación del GPR ......................................................................................................... 40
3.4.5 Parámetros Electromagnéticos .................................................................................... 41
3.4.6 Permitividad Eléctrica .................................................................................................. 42

vi
 3.4.7 Velocidad de Propagación y Longitud de Onda............................................................. 42
3.4.9 Velocidad de Propagación y Longitud de Onda............................................................. 42
3.4.10 Procesamiento de Información de GPR ...................................................................... 44
3.4.11 Filtros ........................................................................................................................ 44
3.4.12 Filtros Horizontales .................................................................................................... 44
3.4.13 Filtros Verticales ........................................................................................................ 44
3.4.14 Deconvolución ........................................................................................................... 45
3.4.15 Migración .................................................................................................................. 45
3.4.16 Transformada Rápida de Fourier (TDF) ....................................................................... 46
3.4.17 Transformada Hilbert................................................................................................. 46
4 Adquisición de Datos ................................................................................................................. 47
4.1 Levantamiento Gravimétrico y Magnético .......................................................................... 47
4.1.1 Establecimiento de estación base .................................................................................... 48
4.2 Procedimiento.................................................................................................................... 49
4.3 Instrumentación ................................................................................................................. 49
4.3.2 Gravimétrico................................................................................................................ 49
4.3.2 Magnetometría ........................................................................................................... 50
4.4 Levantamiento topográfico ................................................................................................ 52
4.5 Levantamiento con radar de penetración ........................................................................... 54
4.5.1 Procedimiento ................................................................................................................. 55
4.5.2 Instrumentación .............................................................................................................. 55
5 Resultados ................................................................................................................................ 57
5.1 Levantamiento magnético .................................................................................................. 57
5.2 Reducción al polo ............................................................................................................... 58
5.3 Dominios magnéticos DM................................................................................................... 59
5.3.1 DM I ............................................................................................................................ 59
5.3.2 DM II ........................................................................................................................... 59
5.3.3. DM III ......................................................................................................................... 60
5.4 Primera derivada en Z ........................................................................................................ 60
5.5 Señal analítica .................................................................................................................... 62
5.5.1 DMSA I ........................................................................................................................ 65
5.5.2 DMSA II ....................................................................................................................... 65
5.5.3 DMSA III ...................................................................................................................... 65

vii
 5.6 Análisis Gravimétrico .......................................................................................................... 65
5.6.1 Análisis Regional .......................................................................................................... 66
5.6.2 Anomalías residuales ................................................................................................... 68
5.7 Resultados GPR .................................................................................................................. 72
5.7.1 Ubicación perfiles GPR ................................................................................................. 73
5.7.3 Interpretación general de Radargramas ....................................................................... 74
5.7.4 Calculo de profundidades a partir de velocidades de propagación ............................... 75
5.8 Modelo de interpretación................................................................................................... 80
5.8.1 Modelo de interpretación con base en radargramas 1 ..................................................... 81
5.8.2 Modelo de interpretación con base en radargramas 2 ..................................................... 82
5.8.3 Modelo de interpretación con base en radargramas 3 ..................................................... 83
5.9 Conclusiones .......................................................................................................................... 84
Referencias .................................................................................................................................. 85
ANEXO A ...................................................................................................................................... 88
Anexo B ....................................................................................................................................... 89
Anexo C ....................................................................................................................................... 90

viii
Lista de Figuras

Figura 1. Localización de la zona de estudio Cerro de San Pedro, SLP ............................... 4

Figura 2. Ubicación fisiográfica de la zona de estudio (Tomado de Nieto-Samaniego,1999)
...................................................................................................................................................... 5
Figura 3. Modelo de elevación digital de la República Mexicana indicando ubicación de la
Mesa Central, los colores naranjas señalan las cotas superiores a 2000 metros y en color
rojo las menores a 2000 metros. (Tomado de Nieto Samaniego, 1999) ............................... 6
Figura 4. Ubicación de la Plataforma Valles-San Luis Potosí. ................................................ 9
Figura 5, Columna estratigráfica local Cerro de San Pedro. ................................................. 10
Figura 6. Mapa geológico Regional del municipio de Cerro de San Pedro, SLP. Tomado y
modificado de (INEGI, 1998)................................................................................................... 13
Figura 7. Campos potenciales de la Tierra. ............................................................................ 15
Figura 8. Generación del Campo magnético terrestre. .......................................................... 16
Figura 9. Inclinación magnética, representa el ángulo entre una línea imaginaria orientada
de acuerdo con la posición de los polos magnéticos y la horizontal (superficie terrestre).
(b) La declinación magnética de la Tierra es el ángulo entre el norte geográfico y las
líneas del campo magnético. ................................................................................................... 17
Figura 10. Componentes del campo magnético total. (Tomado y modificado de Lillie,
1999) .......................................................................................................................................... 17
Figura 11. Representación esquemática de la expansión del fondo oceánico. .................. 19
Figura 12. Diferencias entre Geoide y elipsoide. ................................................................... 25
Figura 13. Relación geométrica entre Elipsoide y esferoide. ................................................ 26
Figura 14. EL geoide como referencia en una superficie equipotencial (Telford et al.,1990).
.................................................................................................................................................... 26
Figura 15. Estimación de denisdad utilizando el método de Nettelton, (Telford et al., 1990)
.................................................................................................................................................... 28
Figura 16. Corrección por deriva a lo largo de un perfil donde se realiza la corrección a
partir de una estación base. ..................................................................................................... 30
Figura 17. La variación de la gravedad debido a la posición con respecto a la latitud. ...... 32
Figura 18. Correcciones por altitud (Tomado y modificado de Cantos,1974) ..................... 33
Figura 19. Corrección de exceso de masa bajo el punto B................................................... 35
Figura 20. Corrección por terreno. (Burger, 1992) ................................................................. 35
Figura 21. Esquema de hipótesis de Pratt (A) y Airy (B). ...................................................... 36
Figura 22.- Levantamiento microgravimetría. ......................................................................... 48
Figura 23. Gravímetro modelo CG-5 (a) vista lateral (b) vista de planta. ............................ 50
Figura 24. Magnetómetro modelo GSM-19. ........................................................................... 51
Figura 25. Levantamiento magnetométrico terrestre en Cerro de San Pedro, SLP. En la
figura A se puede observar cableado eléctrico que debe ser tomado en cuenta a la hora
del procesamiento debido a que provoca ruido de la señal. En la figura B se observa una
zona libre de cableado donde existe poca influencia de ruido. ............................................ 52
Figura 26. Levantamiento topográfico utilizando GPS diferencial. ....................................... 53
Figura 27. Esquema del funcionamiento del GPR. ................................................................ 54
Figura 28. Levantamiento líneas de GPR usando una antena de 100Hz. ........................... 56

ix
Figura 29. Mapa del Campo magnético Residual de Cerro de San Pedro,SLP.................. 57
Figura 30. Mapa Campo magnético reducido al polo Cerro de San Pedro, SLP. ............... 58
Figura 31. Mapa del Campo magnético reducido al polo Cerro de San Pedro, SLP .......... 59
Figura 32. Primera derivada en “Z” Cerro de San Pedro, SLP ............................................. 61
Figura 33. Principales lineamientos marcados con base en el mapa de primera derivada
en”Z”. .......................................................................................................................................... 62
Figura [Link] señal analítica Cerro de San Pedro, SLP .................................................. 63
Figura 35. Dominios magnéticos con base en el mapa de señal analítica de Cerro de San
Pedro, SLP. ............................................................................................................................... 64
Figura 36. Figura A Mapa de elevación topográfica. Figura B Mapa de anomalía de
Bouger completa. ...................................................................................................................... 66
Figura 37. Contraste entre el mapa de Anomalía de Bouger Completa (A) y la anomalía
Regional calculada (B).............................................................................................................. 68
Figura 38. Figura relación entre el MDE (A) y la anomalía Residual (B). ............................ 69
Figura 39. Interpretación del mapa de anomalías residuales (B) con respecto al MDE (A.)
.................................................................................................................................................... 70
Figura 40. Señal analítica Cerro de San Pedro, SLP ............................................................ 71
Figura 41. Gradiente Horizontal, Cerro de San Pedro, SLP ................................................. 72
Figura 42. Distribuciones líneas de radar de penetración a lo largo de la cabecera
municipal de Cerro de San Pedro, SLP .................................................................................. 73
Figura 43. Línea 1 interpretada de Radar de penetración .................................................... 77
Figura 44. Línea 2 interpretada de radar de penetración ...................................................... 77
Figura 45. Línea 4 interpretada de radar de penetración. ..................................................... 78
Figura 46. Línea 5 interpretada de radar de penetración. ..................................................... 79
Figura 47. Modelo geológico-geofísico I con orientación W-E .............................................. 81
Figura 48. Modelo geológico-geofísico II. ............................................................................... 82
Figura 49. Modelo geológico-geofísico III ............................................................................... 83

x
Lista de Tablas

Tabla 1. Susceptibilidad magnética de rocas más comunes. ............................................... 22

Tabla 2. Rango de densidades Rocas y minerales más comunes (Telford er al.,1990) .... 27
Tabla 3. Velocidades de onda materiales y rocas. ................................................................ 43
Tabla 4. Geolocalización inicio y fin de puntos de líneas de radar de penetración ............ 74

xi
Resumen

Los métodos geofísicos son una herramienta importante en la exploración de la

estructura interna de la Tierra, estos estudios pueden abarcar enormes
profundidades o hasta los primeros metros del subsuelo, esto depende el enfoque
y el objetivo de la investigación.
Dentro de esta investigación se concentró el análisis de la estructura somera del
Municipio de Cerro de San Pedro, SLP utilizando los métodos geofísicos potenciales
gravimétrico y magnetométrico y el radar de penetración terrestre GPR (por sus
siglas en inglés.
Estas herramientas permiten generar un enfoque de las condiciones geológicas
presentes en los primeros metros de profundidad de la zona.
La integración de métodos geofísicos permite realizar un análisis detallado cuyos
resultados se presentan de forma gráfica como tres modelos geológicos-geofísicos
2D. Estos modelos indican que el área de estudio es una zona en gran medida
fracturada.
Dentro del análisis se encontró que dentro de la zona estudiada no se hallan
estructuras que se correlacionen con túneles como se había planteado.

Palabras clave: Métodos geofísicos, gravimetría, magnetometría, Radar de

penetración terrestre (GPR), métodos potenciales, modelo geológico-geofísico.

xii
Abstract

Geophysical methods are a very important tool for the exploration of internal
structure of the earth, these studies can cover large depths or just the first meters of
the ground, it depends on the objective or the kind of investigation.

This investigation is major focus on the analysis of the shallow structure of Cerro de
San Pedro city hall of San Luis Potosi applying the gravimetric and magnetometric
potential geophysical methods and the ground penetrating radar (GPR) method.

These tools allow to create an approach about the geological condition of the first
meters of depth of the study zone.

Geophysical method integration enables a detailed analysis whose results are

presented graphically as three 2D geological-geophysical models. These models
point out that the study area is a largely fractured zone.

In the analysis, it was found that within the studied area there are no structures that
correlate with tunnels as had been proposed.

Keywords: Geophysical methods, magnetometric prospection, gravimetric

prospection, ground penetrating radar (GPR), potential geophysical methods,
geological-geophysical model.

xiii
1 Generalidades

1.1 Introducción

La geofísica como herramienta de exploración es una herramienta poderosa para la

interpretación de estructuras del subsuelo que dependiendo el tipo de método y el
objetivo a perseguir se manipula el alcance de medición, que puede ir desde los
primeros metros del subsuelo hasta mediciones a escalas kilométricas.
Dentro de los métodos de gran alcance podemos mencionar a los métodos
potenciales, en particular el método gravimétrico y magnético. Ambos miden
distintas componentes físicas, y estas se pueden interpretar después según la
litología presente o estructuras geológicas.
Por otro lado, el método de radar de penetración es un método limitado en cuanto
penetración y alcance. El cual consiste en la emisión de ondas emitidas por una
fuente y estas viajan a través del subsuelo para su posterior recepción dando
información de las condiciones físicas del subsuelo.
Mediante una integración de los métodos descritos con anterioridad se puede llegar
a una correlación de ambas escalas y brindar una interpretación mucho más
detallada y completa.
La importancia de la presencia de la geofísica de exploración y el estudio de las
geociencias en asentamientos es sumamente necesaria. Estas investigaciones
aportan información estructural que puede ser utilizada para definir las carencias y
la vulnerabilidad a la que se es perceptible, así como proporcionar explicaciones en
caso de acontecimientos extraños correspondientes a la estructura sobre la cual
yace el poblado.
La investigación que aquí se presenta fue realizada con base en mediciones
magnéticas, gravimétricas y de georradar, relacionadas con el campo de atracción
gravitacional y el campo magnético terrestre, así como mediciones
electromagnéticas realizadas por el método de radar de penetración terrestre,
Ground penetrating radar (GPR, por sus siglas en ingles), cuyos datos fueron
integrados para su posterior interpretación.

1
1.2 Objetivos

1.2.1 General

Verificar cualitativamente las condiciones estructurales del subsuelo del

municipio de Cerro de San Pedro en relación con la posible presencia de túneles
dentro de la zona de estudio.

1.2.2 Específicos

Realizar un estudio mediante el uso y la implementación de métodos geofísicos a lo

largo de la cabecera municipal de Cerro de San Pedro, SLP donde se puedan
describir las propiedades y condiciones de la estructura somera de esta locación,
interpretadas a través de un modelo geológico-geofísico 2D.

1.3 Delimitación del problema

Los conocimientos de las condiciones estructurales del subsuelo establecen los

principios necesarios para la construcción o la remediación de edificaciones
dañadas por algún acontecimiento relacionado debido a estas condiciones del
terreno manifestadas de manera posterior a la construcción. El municipio de Cerro
de San Pedro, SLP es fácilmente reconocido por sus antecedentes mineros que
datan desde la época de la colonia, algunos recuentos históricos señalan que dentro
de este periodo de explotación y desarrollo se construyeron además de
edificaciones, como lo son los templos de Cerro de San Pedro y San Nicolás de
Tolentino, por debajo de ellos se instauraron una serie de túneles que se conectan
entre ellos y algunos otros puntos que fueron de gran interés para los pobladores
de esa época.
Agrietamientos, daños en la estructura y eventos relacionados con micro-sismicidad
han estado presentes de forma más recurrente en el municipio de Cerro de San
Pedro, es por ello que se plantea el estudio de la estructura somera de esta zona

2
para verificar la existencia de estos túneles o descartar esta posibilidad mediante la
aplicación de métodos geofísicos.
Cabe mencionar que la explotación minera de manera prolongada puede ocasionar
daños en edificaciones y en el subsuelo.
Con base en un modelo geológico-geofísico se busca estudiar la estructura somera
de CSP para encontrar posibles causas de las afecciones que sufren las
edificaciones de la zona.

1.4 Hipótesis

La liberación de energía producida durante la explotación minera, constantes

excavaciones pueden provocar desestabilidad entre la superficie del terreno y el
subsuelo, esto se traduce en afecciones materiales al mismo terreno y edificaciones
presentes en él.
Las respuestas del radar de penetración pueden servir como interpretación de
presencia de túneles o vacíos presentes en el subsuelo a un nivel somero del
subsuelo, que es la posible ventana de localización de estos túneles.

1.5 Metodología

Este trabajo fue realizado en cuatro etapas:

Se recopiló información geológica y geofísica del área de estudio de carácter
regional, conjunto a ella se revisaron antecedentes de estudios previos en la
zona, esta recopilación de información sirvió como punto de partida del
estudio, herramienta para entender mejor los procesos geológicos que ha
pasado.

Se llevó a cabo una visita de la zona de estudio con la finalidad de decidir la

logística de adquisición de datos, posteriormente se realizó la primera
campaña de levantamiento de datos correspondiente a magnetometría y
adquisición de datos de radar de penetración, en una segunda campaña se
obtuvieron los datos gravimétricos.

En esta tercera etapa se realizó el procesado de los datos gravimétricos, de

GPR y magnéticos, se realizó un análisis del tipo de respuesta obtenido
relacionando este con el tipo de estructuras geológicas presentes, para así
marcar los lineamientos principales observados. Llegar a estas imágenes
fue el resultado de aplicar una serie de filtros para limpieza de ruido de los
datos adquiridos.

3
 La cuarta etapa fue el análisis cuantitativo y la integración de métodos
para la generación de modelos, para esto se ubicaron una serie de
perfiles para obtener modelos geológicos-geofísicos 2D locales, con los
cuales se pudieron dar resultados más concretos acerca de la estructura
somera del municipio de Cerro de San Pedro, SLP.

1.6 Localización

La Zona de estudio del presente trabajo se encuentra dentro del municipio de Cerro
de San Pedro (CSP), ubicado en la parte SW del Estado de San Luis Potosí entre
los paralelos 22° 16’ y 22° 08’, los meridianos 100° 42’ y 100° 54’ de longitud oeste,
con una altura promedio entre 2040 metros sobre el nivel del mar.
El estudio se concentra dentro de la cabecera municipal del municipio de Cerro de
San Pedro, realizándose mediciones a lo largo de cada calle transitable.

Figura 1. Localización de la zona de estudio Cerro de San Pedro, SLP

4
1.7 Ubicación fisiográfica

Fisiográficamente el área de estudio se encuentra localizada dentro de la Mesa

Central, la cual se encuentra delimitada en la parte Este por El sistema de fallas
Taxco-San Miguel de Allende dispuesta de Norte a Sur y de NW a SE por el Sistema
San Luis-Tepehuanes. En la parte Norte colinda con la provincia de la Sierra Madre
Oriental, en la parte Oeste por la Sierra Madre Occidental y en el lado Sur con el eje
Neovolcánico.

Figura 2. Ubicación fisiográfica de la zona de estudio (Tomado de Nieto-Samaniego,1999)

La Mesa Central (MC) se localiza en la parte centro-norte de México más de la mitad

de su superficie se encuentra por encima de la cota 2 000 msnm y las elevaciones
topográficas en su interior son moderadas, la mayoría forman desniveles inferiores
a los 600 m, descrita como una zona elevada de amplias llanuras y sierras
dispersas, constituido principalmente por una cubierta de rocas volcánicas
cenozoicas (Hérnandez- Sánchez Mejorada, 1968)

La MC fue propuesta como provincia fisiográfica por Raisz en 1959. Dentro de la

MC se pueden reconocer dos regiones. La región sur, donde se localizan las cotas
más elevadas, se encuentra por encima de los 2 000 msnm. Esa región de la MC

5
es una región montañosa, cubierta en su mayor parte por rocas volcánicas
cenozoicas (Nieto-Samaniego et al., 1999) (Nieto Samaniego & Alanís-Álvarez,
2005). En contraste, la región norte, que es la de mayor tamaño, muestra un estado
de erosión más avanzado con grandes cuencas continentales rellenas de
sedimentos aluviales y lacustres.

Figura 3. Modelo de elevación digital de la República Mexicana indicando ubicación de la Mesa Central, los colores naranjas
señalan las cotas superiores a 2000 metros y en color rojo las menores a 2000 metros. (Tomado de Nieto Samaniego,
1999)

1.8 Antecedentes

Algunos de los estudios realizados dentro de la zona de estudio comprenden

estudios de tipo geológico-minero, además como los de impacto ambiental y social,
estos últimos debido a la disputa de la operación de la minera dentro del municipio
de Cerro de San Pedro.
Se encontró que los depósitos metalíferos de tipo skarn son especialmente
abundantes en los límites de la Mesa Central, es ahí donde las condiciones son
esenciales para su formación (J.J, 1978) (Rubin, 1988).
La mineralización presente de tipo skarn es formado por la intrusión de un cuerpo
monzodiorítico en las calizas cretácicas de la Formación Cuesta del Cura; así como,
en las Formaciones Soyatal y Tamaulipas (Coremi, 1996) (Petersen, 2001).

6
En la microcuenca Cerro San Pedro, a lo largo de más de 400 años han existido
diferentes periodos de actividad minera, los principales elementos de interés
económico del yacimiento mineral han sido Au, Ag, Zn y Pb (Caddey, 1993).

Dentro de un estudio de lixiviados y modelación hidrogeoquímica se encontró que

las unidades geológicas existentes en la microcuenca incluyen rocas sedimentarias
del Mesozoico: i) Formación La Peña, ii) Formación Cuesta del Cura, que
constituyen el complejo basal o basamento (Aranda, 1989) y una cubierta
Cenozoica que incluye: i) Formación Cenicera, ii) Pórfido San Pedro, iii) Andesita
Casita Blanca, iv) Ignimbrita Santa María, v) Latita Portezuelo, vi) Ignimbrita
Panalillo (Labarthe y Tristán, 1978).
El pórfido fue emplazado dentro de una zona de dilatación formado con la
intersección de la parte oeste de empujes e inversos esfuerzos axiales planares
(Winterbourner, 1999).
Trabajo más reciente indica que el pórfido de San Pedro fue emplazado durante
deslizamientos inversos a través de posteriores fallas compresionales (Keith, 2004).

7
2 Marco Geológico

2.1 Geología Regional

Cerro de San Pedro (CSP) se encuentra localizado dentro del napa de Parras * el
cual corresponde a una sección de la parte Este del cinturón plegado de la Sierra
Madre Oriental del noreste de México y se extiende al occidente a lo largo del límite
fisiográfico con la Sierra Madre Occidental (SMOc) dentro de una secuencia de
calizas plegadas de edad cretácica y al oriente limita con la Mesa central, esta
secuencia es cortada por un pórfido diorítico de edad terciario temprano, debido al
contacto entre las calizas y el pórfido diorítico, estas presentan cierto grado de
mineralización.

8
2.1.2 Plataforma Valles San Luis Potosí

La Plataforma Valles-San Luis Potosí (PVSL) es una unidad paleogeográfica de

gran extensión que abraca geográficamente los estados de Tamaulipas, San Luis
Potosí, Hidalgo, Querétaro de edad mesozoica.
La PLVS consiste en un segmento de rocas precámbricas, paleozoicas y triásicas,
plegadas y falladas sobre la cual se depositaron una secuencia delgada de
sedimentos marinos del Jurásico Superior y Cretácico.

Figura 4. Ubicación de la Plataforma Valles-San Luis Potosí.

9
Figura 5, Columna estratigráfica local Cerro de San Pedro.

10
2.1.3. Formación la Peña

La Formación la Peña (Aptiano Superior), de cuenca se encuentra conformada por

intercalaciones de calizas de tipo mudstone de estratos delgados a medios,
intercalados con laminaciones y capas delgadas de lutitas y limolita con nódulos y
lentes de pedernal negro, mientras que la Formación la Peña, de talud se
caracteriza por la intercalación de capas brechosas con fragmentos de fósiles y
pedernal, en matriz calcárea, además esta presenta mayor cantidad de capas de
terrígenos, muestra estructuras “slump”, nódulos, lentes y capas de pedernal.

2.1.4. Formación Cuesta del Cura

La Formación Cuesta del Cura (Albiano) está constituida por sedimentos de cuenca
que se denotan por ser una interestratificación de capas delgadas de calizas
micríticas y bandas de pedernal en una secuencia rítmica, en donde los sedimentos
de talud difieren a los de cuenca por la presencia de horizontes de capas brechosas
donde los fragmentos pertenecen a corales, rudistas, y de otra procedencia distinta
a las del arrecife.
2.1.5. Formación Indidura

La Formación Indidura (Turoniano), descansa de forma transicional sobre la

Formación Cuesta del Cura y se distingue por ser una secuencia de capas de
calizas arcillosas de espesor delgado a medio de color gris claro a ligeramente
rojizas intercaladas con limolita del mismo color.
2.1.6. Formación Cárdenas

Sobre yace de la Formación Indidura, se encuentra conformada por una secuencia

de estratos medianos de lutitas, areniscas y limolitas, la cual solo se puede ver en
el fono de los arroyos y está cubierta por carpeta de conglomerados.

2.1.7. Formación Casita Blanca

La Formación Casita Blanca (Eoceno), con sedimentos de tipo lacustre y

conglomerados polimícticos que se depositaron debido a un evento volcánico que
ocasionó el emplazamiento de derrames de lava, y en menor proporción flujos
piroclásticos.

11
2.1.8. Ignimbrita Santa María

La Ignimbrita Santa María, se caracteriza por la presencia de fragmentos alargado

de pómez bien soldados.

2.2. Geología y estratigrafía local

La sección estratigráfica local está compuesta por calizas de edades Cretácicas,

Formación La peña (Aptiano) y Cuesta del cura (Albiano) (Bello-Barradas, 1963),
(Labarthe, 1982) representan depósitos de aguas profundas y de cuenca
depositados a lo largo del margen occidental de la Plataforma Valles- San Luis
(Carrillo, 1971), (Aranda-Gómez, 2005) que han sido intrusionadas por el pórfido
diorítico de San Pedro, un cuerpo ígneo intrusivo monzodiorítico de edad del
Cretácico Tardío al Terciario Temprano emplazado por una falla de escala regional
de tipo inverso.

La litología está comprendida por capas de calizas tipo mudstone que varían de
espesores delgados a medianos, y pueden ser separados en dos tipos de facies. La
primera, unidades de cuenca y de talud. La Formación la Peña, de cuenca, se
encuentra conformada por una intercalación de calizas micríticas de estratos
delgados a medios, intercalados con laminaciones y capas delgadas de lutita y
limolita con nódulos y lentes de pedernal negro, y la de talud se distingue por la
intercalación de capas brechosas con fragmentos de fósiles, y pedernal, en matriz
calcárea.

12
Figura 6. Mapa geológico Regional del municipio de Cerro de San Pedro, SLP. Tomado y modificado de (INEGI, 1998).

13
3 Marco Teórico

3.1 Métodos de Prospección Geofísica una breve introducción.

La geofísica comprende el estudio de la tierra en su composición, esto gracias a

mediciones de tipo físico que normalmente son llevados a cabo en la superficie.
Los métodos geofísicos se basan en la medición de las propiedades físicas de los
elementos que componen el subsuelo (rocas, sedimentos, agua, cavidades,
objetos, etc.) y se pueden clasificar en dos tipos: pasivos y activos.

Pasivos, son aquellos que se basan en la detección de las variaciones de las

propiedades naturales del terreno, sin introducir señales artificiales. Ejemplos de
estos podrían ser: Gravimetría, Potencial espontáneo, Magnetometría, Geoquímica,
Radiometría, Geobotánica y Geobiología, entre otros.

Activos, son aquellos que estudian el comportamiento de los materiales del

subsuelo mediante la creación de señales artificiales y viendo cómo estas señales
interaccionan con dichos materiales. Ejemplos de estos podrían ser: Polarización
inducida, Métodos sísmicos (sísmica de refracción, sísmica de reflexión), métodos
eléctricos (SEV, tomografía, métodos geoeléctricos), métodos electromagnéticos,
georradar, entre otros

3.2 Método magnético

El método magnético es una herramienta de exploración geofísica que consiste en

realizar mediciones de las variaciones en el campo magnético total terrestre (F) y
en base en la respuesta del medio a través de anomalías magnéticas poder ser
interpretadas y correlacionadas según el tipo de unidades geológicas a las que
estas correspondan.

14
3.2.1 El campo magnético terrestre

El campo magnético de la Tierra, también conocido como campo geomagnético se

comporta de manera semejante al campo generado por una barra magnética, es
como si tuviera un enorme imán de barra en su centro (núcleo), debido a que el
campo magnético generado por la tierra se comporta como un gran imán este
expone líneas de flujo que penetran el polo norte magnético y convergen por el polo
sur magnético.

Figura 7. Campos potenciales de la Tierra.

No obstante, las temperaturas dentro del interior de la Tierra son demasiado

elevadas, lo que impide la magnetización de manera permanente de cualquier tipo
de material que pueda entrar en contacto (Fe, Ni, Ir, Os, etc.) este fenómeno de
anulación de la magnetización es explicado con el entendimiento de la temperatura
de Curie. Precisamente porque la temperatura de Curie del hierro es de +/- 700°C
aproximadamente, y sabemos que las temperaturas del núcleo son mucho mayores
a esta, en consecuencia, el material no presentará propiedades magnéticas.

Debido a la explicación anterior sabemos que la generación del campo

geomagnético por esa vía es imposible, por ello se adoptó una nueva hipótesis que
describe de una manera más coherente el proceso físico, el magnetismo de la Tierra
debe producirse por el efecto del flujo de corrientes eléctricas constantes dentro del
interior de la Tierra, en consecuencia, del movimiento de corrientes convectivas
dentro del núcleo y la rotación exterior del núcleo líquido (efecto de coriolisis).

15
Figura 8. Generación del Campo magnético terrestre.

3.2.2 Componentes vectoriales del campo magnético

El campo magnético esta descrito como un vector que posee magnitud y sentido.
La orientación observada en una brújula exhibe la dirección del campo magnético
terrestre, las componentes del vector del campo magnético se pueden adquirir
mediante una medición en un punto P.

La inclinación magnética (i) es el ángulo entre la línea imaginaria que se encuentra

orientada hacia el campo magnético total (F) y la horizontal o la superficie terrestre
(Figura 9a). Los valores de las líneas verticales en el polo Sur (geográfico y
magnético), donde salen las líneas del campo (hacia abajo) pueden variar con
valores de -90°, mientras que en el Polo Norte donde entran las líneas del campo
(hacia arriba) el valor es de +90°. En el Ecuador, donde las líneas del campo son
paralelas a la superficie, el ángulo que representa la inclinación magnética es de 0°.

Otro componente se encuentra representado por el ángulo que se origina por el

Meridiano Magnético respecto a la dirección del Meridiano geográfico, a este ángulo
se le denomina declinación magnética (Figura 9b). Se estableció que las
declinaciones magnéticas posicionadas al Oeste del meridiano geográfico serán
negativas (D-) y las que se encuentren posicionadas al Este serán positivas.

16
Figura 9. Inclinación magnética, representa el ángulo entre una línea imaginaria orientada de acuerdo con la posición de
los polos magnéticos y la horizontal (superficie terrestre). (b) La declinación magnética de la Tierra es el ángulo entre el
norte geográfico y las líneas del campo magnético.

Partiendo de un Punto P Los elementos del campo magnético son los siguientes:

Figura 10. Componentes del campo magnético total. (Tomado y modificado de Lillie, 1999)

17
La magnitud (F) que corresponde al vector del campo magnético total o también
llamado intensidad total del campo, se explica por la siguiente ecuación:

𝐹 = √𝐹𝐻 2 + 𝐹𝑉 2 (Ec 3.1) Y √𝐹𝑁 2 + 𝐹𝐸 2 + 𝐹𝑁 2 (Ec 3.2)

La inclinación y declinación magnéticas quedan definidas en las ecuaciones 1.1 y

1.2
𝐹
𝑖 = tan−1 𝐹𝑉 (Ec 3.3)
𝐻

𝐹
𝐷 = tan−1 𝐹𝑉 (Ec 3.4)
𝐻

El modelo de tipo axial genera una respuesta a la discusión de un planeta con un

campo magnético, esto gracias a que las ecuaciones pueden ser desarrolladas para
describir la intensidad (i) y la dirección del campo magnético. Las magnitudes
vectoriales FH y Campo total FV, sus componentes horizontales y verticales
respectivamente pueden ser calculadas, así como el campo Total (F) (Butler, 1992).

𝑀⋅cos 𝜙
𝐹𝐻 = (Ec 3.5)
𝑅3

2𝑀⋅sen 𝜙
𝐹𝑉 = 𝑅3
(Ec 3.6)

En donde:
R= Radio Terrestre
M= Intensidad del campo total en el Ecuador.
𝞥= Latitud magnética con una inclinación verdadera de 11°
La inclinación magnética para un dipolo axila también varia de forma sistemática
con la latitud:
𝑡𝑎𝑛 𝑖 = 2𝑡𝑎𝑛𝜙 (Ec 3.7)

18
 1
A diferencia con el método gravimétrico donde la razón de decaimiento es , la
𝑅2
1
intensidad del campo magnético decae por u factor de .
𝑅3

3.2..3 Fuentes principales del campo magnético terrestre

La principal fuente generadora del campo magnético terrestre o geomagnético,

puede ser encontrada en el núcleo terrestre, caracterizado como un dipolo, con la
dirección del dipolo algo inclinada con respecto al eje de rotación (aproximadamente
11), una cantidad importante de componentes no dipolares procedentes de la misma
tierra hacen una contribución importante como influencia de los complejos sistemas
de corriente ionosféricas y magnetosfericas.

El campo llamado campo principal es aquel generado por el núcleo y tiene una
variación temporal denominada variación secular, la variación secular actual se
caracteriza por una disminución en la magnitud del momento dipolar que es
aproximadamente de 0.5% por año por una rotación hacia el Oeste del Vector
dipolar de 0.05° de longitud por año.
Considerando este campo en la escala de tiempo geológico (millones de años, Ma)
se sabe que ha experimentado numerosas inversiones de polaridad, cuya
naturaleza todavía no ha podido explicarse claramente, pero que se han podido
registrar hasta una antigüedad de 170 millones de años. En la Figura 11, están
representadas esquemáticamente las inversiones de los últimos 4.5 Ma.

Figura 11. Representación esquemática de la expansión del fondo oceánico.

En la capa superior fría y sólida de la tierra, la corteza, se encuentran rocas con

propiedades magnéticas diversas. El campo principal actúa sobre ellas, generando

19
magnetización que es básicamente de dos tipos: inducida, que está determinada
por el campo magnético actual, o remanente, que guarda información del campo
presente en la época en que la roca se formó. Estos dos tipos de magnetización se
deben a magnetización espontánea, una propiedad compleja de los minerales
ferromagnéticos en la corteza terrestre (Blakely, 1996).

La magnetización espontánea depende de la temperatura. A medida que el material

se calienta, el espaciamiento entre momentos magnéticos vecinos a nivel atómico
aumenta, hasta que se alcanza un punto en que la magnetización espontánea se
anula. La temperatura a la cual esto ocurre se denomina temperatura de Curie.
Entonces, tanto la magnetización inducida como la remanente desaparecen a
temperaturas mayores que la de Curie. Aun cuando algunos efectos
paramagnéticos y diamagnéticos persisten a estas temperaturas, para un estudio
de la corteza terrestre se puede considerar que las rocas con temperatura mayores
que la de Curie son esencialmente no magnéticas.

En el espacio exterior, el campo creado por la Tierra interactúa con el flujo de

partículas cargadas que constituyen el viento solar y resulta severamente
distorsionado, a esta cavidad "cerrada" se le denomina magnetósfera.

Dentro de la magnetósfera existen complejos sistemas de corrientes, que dan lugar

a variados fenómenos electromagnéticos, lo cuales están directamente
relacionados con la actividad solar, cuya variación temporal típica se expande en un
espectro bastante amplio: desde los segundos (micro pulsaciones), hasta períodos
típicos de los ciclos solares. En la región comprendida entre los 50 y 1500 km de
distancia a la superficie terrestre se encuentra la ionósfera, formada por partículas
ionizadas como resultado de la acción de la radiación ultravioleta proveniente del
sol sobre los constituyentes atmosféricos. Las corrientes debidas a estas partículas
crean a su vez componentes magnéticas importantes. El estudio del campo externo
es de fundamental importancia para el conocimiento del espacio exterior a la tierra
(necesario para las comunicaciones, la navegación espacial, los satélites 81
artificiales, etc.), y también por la influencia que los eventos electromagnéticos
tienen en la vida terrestre. Las erupciones solares provocan las llamadas tormentas
magnéticas, que pueden afectar gasoductos, líneas de transmisión, etc.

El campo electromagnético externo induce campos eléctricos y magnéticos en la

tierra, dada su conductividad. Estos campos afectan tanto a la corteza como al
manto, y son también empleados para inferir las propiedades eléctricas y
magnéticas. Con lo anterior se puede concluir que cada medición del campo
magnético realizada cerca de la superficie de la tierra nos da como resultado la
suma de 4 principales fuentes de distintos orígenes:

1. El núcleo, la magnitud del campo magnético dipolar es de aproximadamente

50000 nT en los polos y 25000 nT en el ecuador (1 nT = 10-9 Tesla; 1 Tesla = 104
Gauss).
2. Las rocas magnetizadas, cuya intensidad del campo magnético puede llegar a
los 1000 nT.
20
3. Las corrientes externas, que pueden dar valores del campo magnético entre unos
pocos nT, en los días solares tranquilos, hasta alrededor de 1000 nT en las
tormentas magnéticas. Estos valores en levantamientos de campo pueden ser
consideradas como ruido e incluso hacer que los datos levantados sean inservibles.
4. Las corrientes inducidas en la corteza y el manto, cuyos valores del campo
magnético dependen de la conductividad de las rocas.

3.2.4 Variaciones del campo geomagnético

Las variaciones en el campo magnético terrestre pueden verse afectadas por

diferentes factores, y composición de los materiales, estas variaciones sufridas
corresponden en el dominio del tiempo y el espacio de los diferentes tipos de
materiales tanto internos como externos presentes. Las variaciones del tipo
externas que genera variaciones en el campo magnético terrestre son producto de
la interacción del campo magnético global con el campo magnético asociado con el
viento solar.

Los principales efectos que son de especial interés durante los levantamientos de
magnetometría son enumerados a continuación (Nabighian, 1988) (Grauch, 1966).

• La variación diurna, debe su origen a la acción del viento solar sobre las corrientes
ionosféricas, además de provocar el achatamiento del campo magnético.
terrestre en la parte de la Tierra que está expuesta al Sol. Esta variación es de
aproximadamente 60 nT, en periodos de 24 horas (Telford, 1990).
• La interacción del viento solar con el campo magnético terrestre causa la formación
de un anillo de cargas en el ecuador que produce variaciones de ~10 nT en latitudes
cercanas a él.
• Transiciones del campo magnético terrestre en los polos, con permanencia de
algunas horas, por la interacción de partículas cargadas y las líneas de flujo del
campo terrestre que concurren durante las tormentas magnéticas (~20 nT).

21
Tabla 1. Susceptibilidad magnética de rocas más comunes.

Tipos de Rocas Susceptibilidad magnética X 10-3 (SI)

Sedimentarias Rango Promedio
Dolomita 0-0.9 0.1
Caliza 0-3 0.3
Arenisca 0-20 0.4
Lutita 0.01-15 0.9
Metamórficas
Anfibolita 0.7
Gneiss 0.1-25
Cuarcita 4
Serpentinita 3-17
Ígneas
Granito 0-50 2.5
Riolita 0.2-35
Gabro 1-90 70
Basalto 0.2-35 70
Diorita 1-90 85
Peridotita 0.6-120 150
Andesita 90-200 160

3.3 Método gravimétrico

3.3.1. Introducción

Dentro de este apartado se darán a conocer las principales bases teóricas del
método gravimétrico; así como, sus fundamentos más importantes como método de
exploración geofísico no destructivo.
La gravedad es una de las principales fuerzas presentes en todo el universo tanto
a una escala macroscópica y microscópica, cuando se habla de la fuerza de
atracción de la gravedad relacionamos las principales componentes presentes, que
corresponden a la masa y la distancia que separe a dichas masas de objetos, donde
la cantidad de masa de dichos objetos estará relacionado con el producto de sus
masas, esto quiere decir que entre mayor masa tenga un objeto este tendrá mucha
mayor fuerza de atracción o capacidad para atraer a otros objetos y será
inversamente proporcional a la distancia que separe a los objetos. La Ley de
Gravitación Universal formulada por Isaac Newton describe perfectamente la
relación anterior entre masa y distancia de separación, la cual enuncia que Ley
de gravitación universal formulada por Isaac Newton, afirma que la fuerza de
atracción que experimentan dos cuerpos dotados de masa es directamente

22
proporcional al producto de sus masas e inversamente proporcional al cuadrado de
la distancia que los separa.

El método gravimétrico es uno de los métodos geofísicos más empleados, el cual

consiste en medir las variaciones del campo gravitacional de la Tierra. Estas
variaciones son generadas por de las diferencias de densidad entre las rocas del
subsuelo. El método está basado en la ley de Gravitación Universal decretada por
Newton, la cual establece que la fuerza de atracción, F, entre dos partículas o
cuerpos con masas m1 y m2, separadas por una distancia r, es representada por la
siguiente expresión (Burger, Exploration Geophysics of the Shallow Subsurfac,
1996).
𝑚1 𝑚2
𝐹=𝐺 (Ec 3.8)
𝑟2

Donde, G es la constante de gravitación universal con un valor de 6.673 x 10-11

(NM2/Kg2), ésta fue determinada por primera vez en 1798 por Lord Cavendish.
En donde se asume que la Tierra es esférica, por lo tanto, la fuerza ejercida por ésta
sobre un cuerpo esférico, de masa m en reposo, sobre su superficie es (Burger,
Introduction to Applied Geophysics, 2006)(Burguer et al., 2006):
𝑚𝑀
𝐹 = 𝐺 𝑅2 (Ec 3.9)

Donde, M es la masa de la Tierra y R su radio. Esto también supone que la densidad

varía únicamente con la variación de la distancia al centro de la Tierra.

En el caso de la fuerza de aceleración, ésta también es dada por la segunda Ley de

Newton, la cual establece que ésta es directamente proporcional a la masa (m) y a
la aceleración (a):
𝐹 = 𝑚𝑎 (3.10)

cuando la aceleración “a” es causada únicamente por el efecto de atracción

gravitacional de la Tierra, ésta es la constante de aceleración gravitacional “g”
(Burguer et al., 2006)”, por lo que para deducirla se tiene que:

𝑚𝑀
𝐹 = 𝑚𝑔 = 𝐺 (3.11)
𝑅2

Entonces la Aceleración gravitacional queda definida por:

𝐺𝑀
𝑔 = 𝑅2 (3.12)

23
La magnitud de la aceleración es igual a la gravedad en la tierra (g), y es
directamente proporcional a la masa de la Tierra (M) e inversamente proporcional
al cuadrado del radio de la Tierra R (ecuación 3.15). En teoría, la gravedad sobre la
superficie de la Tierra es de 9.8 m/s2, pero en realidad varía de un lugar a otro debido
a que la tierra no tiene la forma de una esfera perfecta y tiene una topografía
irregular; además, las masas de cada cuerpo son variables especialmente cerca de
la superficie, ya que a profundidad tienen un comportamiento más constante.

Las dimensiones de g quedan establecidas como 𝑚⁄𝑠 2 dentro del Sistema

Internacional de unidades (SI) y para el sistema CGS queda como 𝑐𝑚⁄𝑠 2

Dentro de la geofísica las variaciones monitoreadas de la gravedad (g) son

cantidades muy pequeñas, es por eso que se utiliza el Gal, la unidad de aceleración
de gravedad, el cual es proporcional a 1 𝑐𝑚⁄𝑠 2 . La unidad gravimétrica es más
común encontrarla como mili Gal, donde: 1 miliGal es igual a 0.001 Gal

3.3.2 Medición de la Gravedad

Considerando el achatamiento de la Tierra, esta causa un aumento en la gravedad

del ecuador al polo. Así podemos dividir las mediciones de gravedad en dos clases:
Absolutas y relativas, el primer tipo de medición corresponde a la determinación
absoluta de la gravedad en cualquier punto o lugar de la Tierra; la segunda
corresponde en registrar el cambio de gravedad de un lugar con respecto a otro.
Como se mencionó con anterioridad en la exploración gravimétrica es necesario
obtener mediciones de gran precisión para los cambios de gravedad causados por
estructuras del subsuelo, esto se obtiene mediante el uso de un instrumento
electrónico capaz de percibir señales del orden de 0.01 mGal, la principal
importancia del uso de este tipo de instrumentos es la adquisición de valores de
gravedad más acertados para la determinación de anomalías gravimétricas para su
posterior interpretación y correlación geológica.

3.3.3 El geoide

24
Es muy bien sabido que el planeta tierra al igual que otros planetas no presentan
una forma esférica, sino que debido a efectos de la gravitación y de la fuerza
centrífuga producida al rotar sobre su mismo eje genera un achatamiento polar y
ensanchamiento ecuatorial; a esta formo volumétrica se le conoce como esferoide.

Figura 12. Diferencias entre Geoide y elipsoide.

Considerando que si se tomara en cuenta la corteza terrestre no se podría hablar

de que la Tierra podría ser representada como un geoide, solo tomamos a
consideración el nivel medio de las mareas, si se pudiera considerar un geoide, un
geoide es una superficie de nivel que representa al nivel medio del mar, la cual se
prolonga por debajo de los continentes y cubre a la Tierra en su totalidad. Puede
ser imaginada como la superficie del mar en condiciones ideales de quietud y es en
todo punto perpendicular a la línea de plomada o dirección de la gravedad.

El geoide es un modelo físico que busca representar la verdadera forma de la Tierra

calculándola como una superficie del campo de gravedad con potencial constante y
es utilizada como referencia para determinar la elevación del terreno. Las imágenes
incluidas en esta página muestran la ubicación del geoide con respecto de la
superficie topográfica y el elipsoide geodésico de referencia.

25
Figura 13. Relación geométrica entre Elipsoide y esferoide.

Figura 14. EL geoide como referencia en una superficie equipotencial (Telford et al.,1990).

3.3.4 Densidad de las Rocas

Para poder entender e interpretar de una forma más coherente las anomalías
gravimétricas se necesita de información sobre la densidad de las rocas que se
encuentran en el área de estudio. Factores como el intemperismo y la erosión de
las rocas puede afectar en el valor de la densidad de las rocas que afloran en
nuestra zona de estudio, perdiendo competencia y haciéndolas menos densas que
el resto de las rocas aflorando, información de estas puede ser obtenido por
estudios previos como sondeos en pozos o registros radiométricos, si no se cuenta
con aluna información extra de este tipo podemos recurrir al uso de tablas
estandarizadas para los valores más típicos de densidad de rocas y minerales,
como la que se presenta a continuación:

26
Tabla 2. Rango de densidades Rocas y minerales más comunes (Telford er al.,1990)

Rocas Rango de Minerales Rango de

Comunes Densidades densidades
aproximadas
Arena seca 1.4-1.65 Esfalerita 3.8-4.2
Serpentinita 2.5-2.6

Arena 1.95-2.05 Galena 7.3-7.7

húmeda
Gneises 2.65-2.75
Calcopirita 4.1-4.3
Carbón 1.2-1.5

Granito 2.5-2.7
Cromita 4.5-4.8
Yeso 2.1-2.4

Dolerita 2.7-3.1 Pirrotita 4.4-4.7

Sal 2.6-2.7

Basalto 2.7-3.1 Hematita 5.0-5.2

Caliza 2.6-2.7
Pirita 4.9-5.2
Gabro 2.7-3.3

Cuarcita 2.6-2.7
Magnetita 5.1-5.3
Peridotita 3.1-3.4

27
3.3.5 Método de Nettleton

El método de Nettleton es un método empleado para la estimación de la densidad

cerca de la superficie, este cálculo se obtiene mediante la realización de mediciones
gravimétricas en determinados perfiles, áreas con diferentes elevaciones
topográficas, las cuales son correlacionados con la densidad (ρ), el valor promedio
de densidad será el que produzca la mínima correlación (positiva o negativa) entre
la topografía y el mapa de anomalía de Buguer.

Este método proporciona información sobre las densidades únicamente a

profundidades relativamente someras, y puede ser fácilmente usado sólo cuando la
litología cercana a la superficie es homogénea. Ya que en muchas áreas las
características topográficas son producto de afloramientos con litología
heterogénea, y el valor de densidad determinado puede ser incorrecto (Milsom,
2003).

Figura 15. Estimación de denisdad utilizando el método de Nettelton, (Telford et al., 1990)

28
3.3.6 Reducción de datos Gravimétricos

Las suposiciones teóricas que indican que si la Tierra tuviera una simetría esférica
perfecta, que careciera de su movimiento de rotación para anular el efecto de
achatamiento y ensanchamiento en los polos y el Ecuador respectivamente,
además que se tratara de un material completamente homogéneo podríamos hablar
de que la atracción gravitatoria terrestre sería igual en cualquier punto sobre esta,
sabemos que esto es solo existe en un modelo teórico perfecto, por lo tanto
trasladamos esta información para encontrar y estudiar las variaciones que generan
estos cambios y traducirlos en interpretaciones del subsuelo.
Durante la adquisición de datos gravimétricos, estos no arrojan una medición exacta
del valor de gravedad, sino que corresponde a una lectura que ha sido multiplicada
por un factor de calibración Un factor importante es la posición de observación sobre
la superficie de la Tierra:

El valor absoluto de gravedad en el ecuador es 978 Gals, mientras que en los polos
dicho valor es de 983.2 Gals (Burguer et al., 2006). Dicha lectura se le conoce como
gravedad observada (Gobs) A partir de ella se obtienen anomalías de la fuerza de
gravedad sobre la superficie del elipsoide y la magnitud de estas depende de los
cambios en las densidades de las rocas y de los excedentes de altura sobre el
geoide (Mironov, 1977).
Antes de poder interpretar los resultados de un estudio gravimétrico, se debe
corregir todas las variaciones del campo gravitacional de la Tierra que no provengan
de las diferencias de densidad de las rocas (Kearey, 2002), debido a que algunas
de estas variaciones son producidas por fuentes que no son de interés geológico
directo (Milsom,2003).

El proceso inicia con la medición del valor de la G obs y el cálculo a partir de la

formula internacional de gravedad, con referencia a una estación base a partir de la
cual se comienza a hacer los cálculos que se enumeran enseguida:
Corrección por deriva
Corrección por Mareas
Corrección por Latitud
Corrección por Elevación
Corrección por Aire Libre
Corrección de Bouguer
Corrección por Terreno
Corrección isostática

29
3.3.7 Corrección por deriva

Eliminar el efecto de deriva es el primer paso para comenzar con el análisis de datos
gravimétricos. Estos cambios corresponden a variaciones en intervalos de tiempo y
están relacionadas por deriva del mismo instrumento, esto debido a los cambios
elásticos que presenta el resorte, generalmente unidireccionales. Las variaciones
pueden ser observadas cuando se repite una misma medición en un mismo punto,
se recomienda al momento de hacer un levantamiento registrar la primera y última
lectura en una estación designada como estación base para corregir los datos y
evitar esa variación en todos nuestros datos. La corrección por deriva en el tiempo
t es d, la cual es sustraída de los valores observados (Kearey et al., 2002).
La deriva instrumental está a menudo relacionada con la temperatura, y es poco
probable que sea lineal entre dos lecturas consecutivas tomadas en un ambiente
frío al inicio y al final del día, cuando las temperaturas fueron mayores al mediodía.

Para corregir manualmente por la deriva, usando la suposición lineal, las lecturas
primero son corregidas por el efecto de mareas, después se obtiene el coeficiente
de deriva, el cual resulta de la diferencia entre las lecturas inicial y final realizadas
en la base local, este resultado es entonces sustraído de cada lectura a su vez. El
signo del coeficiente de deriva es dictado por el requerimiento después de la
corrección, todos los valores relativos a la deriva deben ser cero. Después de la
corrección por deriva, los valores absolutos de gravedad de cada estación
gravimétrica son obtenidos por la suma de la gravedad absoluta de la base local
con los valores corregidos por deriva.

Figura 16. Corrección por deriva a lo largo de un perfil donde se realiza la corrección a partir de una estación base.

30
3.3.8 Corrección por mareas

Las mareas terrestres provocan variaciones en la gravedad de un poco más de 0.3

mGal, estas variaciones son generadas debido a la interacción del agua de los
océanos con la atracción generada por la de los cuerpos de la Luna y el Sol, también
depende la latitud donde nos encontremos y la hora a la que se realicen las
mediciones.
El Gravímetro CG-5 Scintrex realiza esta corrección en automático, será el
instrumento que se utilice durante nuestras campañas de medición.

3.3.9 Corrección por latitud

La Corrección por Latitud (CLat ), es realiza debido a la variación de la aceleración

de la gravedad, que depende de la latitud, ya que la aceleración es mayor conforme
nos acercamos a los Polos, y es menor en el Ecuador. Como el radio polar
(6357Km) es 21 Km más pequeño que el radio ecuatorial (6378Km), los puntos de
los polos se encuentran más cercanos al centro de masa de la Tierra, entonces el
valor de la gravedad en los polos es más grande ( aproximadamente 0.7%) que en
el ecuador. Y también es afectada, pero en menor proporción, por la fuerza
centrífuga por efecto de rotación, que es menor en los Polos y mayor en el Ecuador,
y tiene un sentido contrario a la fuerza de gravedad.

El principal objetivo de aplicar esta corrección a los datos es corregir la posición en

elipsoide para que sea lo más exacto posible, para ello se han desarrollado fórmulas
para realizarlo por medio de la latitud que medimos en campo.
Existe una ecuación desarrollada para el cálculo de la aceleración teórica (gΦ)
usando la latitud geográfica (Φ), relativa al nivel del mar (g0). También se toman en
cuenta α в β, constantes que dependen de la forma del esferoide y de la velocidad
de rotación de la Tierra.

Ec (3.13)

31
Figura 17. La variación de la gravedad debido a la posición con respecto a la latitud.

En 1930, la International Union of Geodesy and Geophysics adoptó la Fórmula de

Gravedad Internacional, la cual se convirtió en el estándar para trabajar. Sin
embargo, con el tiempo se han tratado de mejorar los cálculos y usando el poder de
las computadoras y el mejor conocimiento de los parámetros de la tierra el resultado
fue la fórmula conocida como Geodetic Referent System 1967(GRS67) y en 1974
la IAG mostró una fórmula más actualizada, la cual se utilizó para el procesamiento
de nuestros datos.
Fórmulas utilizadas para calcular la gravedad normal usando la latitud medida en
cada punto:

32
3.3.10 Corrección por aire libre

Esta corrección se refiere a la intensidad de la gravedad que varía de acuerdo con

la distancia del punto de medición con respecto al centro de la Tierra, es decir, cómo
varía la fuerza de gravedad con la elevación (z) de la estación base a la estación
medida. Para ello, se considera la elevación de la estación base como referencia, y
se le resta a la obtenida en la estación en la cual se quiera hacer la corrección,
obteniendo así la altura h de separación entre ellas (Fig. 18). Esta corrección se
obtiene de la siguiente fórmula:

m1 ∂g m 2g
g=G ; ∂r = −2G r3 = − = −0.3086mgal/m (Ec 3.14)
r2 r

CAL = (−0.3086∗ Δh) (Ec 3.15)

La corrección es positiva si el punto de medicion esta sobre el nivel de la estacion

base, de lo contrario sera negativa.

Figura 18. Correcciones por altitud (Tomado y modificado de Cantos,1974)

33
3.3.11 Corrección de Bouguer

Es similar a la corrección por elevación, sólo que toma en cuenta el material entre
la estación base y la estación medida, aproximando todas las masas arriba o abajo
del plano de referencia con una base homogénea de extensión lateral infinita y de
espesor igual a la altura h de la estación respecto al plano de referencia (estación
base). Esto quiere decir que se calcula hallando el efecto gravitatorio del terreno
comprendido entre el nivel de la estación (plano de la cota h) y el nivel de referencia,
suponiendo que la estación se ha hecho en un plano topográfico horizontal, el efecto
gravitatorio sobre la unidad de masa es
(CB = −2πkδh = −0.419δΔh0.3086∗ Δh) (Ec 3.16)

Donde es la densidad de la corteza terrestre (en promedio es de 2.67 gr/cm3) y K

es la constante de la gravitación universal. La corrección es negativa si la estación
medida está por encima del nivel de referencia, ya que la densidad del cuerpo
aumenta la atracción, entre mayor sea la altura mayor será.

La corrección de Bouguer calcula el empuje extragravitacional ejercido por un

bloque de rocas con espesores de h metros y una densidad principal ρ dada en
kg/m3 o g/cm3 (Fig. 19), el cual encontramos en las mediciones de gravedad (gobs),
siendo calculada por la siguiente cantidad 0.0419ρh mGal/m.

Este valor debe ser calculado a partir de los puntos medidos sobre el nivel del mar.
Se utilizan valores promedios de densidades de la roca de 2.67 g/cm3 para rocas
cristalinas y 2.3 g/cm3 para rocas sedimentarias, dependiendo qué tipo de roca
exista en el subsuelo. En áreas marinas la corrección es ligeramente diferente ya
que la densidad del agua marina es reemplazada por las densidades de los
diferentes tipos de rocas que existen en el área de estudio, como se muestra en el
siguiente cuadro:

𝛿𝑔 = 2𝜋𝐺𝑃ℎ = 𝛽𝑝ℎ (Ec 3.17)

Donde
Β= 2πG= 0.4192 Mg-1 ó 0.0419mGal7m
G= 6.667 x 10-8 m3Mg-1s-2

34
Figura 19. Corrección de exceso de masa bajo el punto B.

3.3.12 Corrección por terreno (topográfica)

Se considera el efecto de atracción de masa dado por la topografía (valles y

montañas) de alrededor de la estación, tanto las situadas por encima de la estación
como las depresiones situadas debajo de la estación, pero como este material en
realidad no existe, debemos sumar su atracción para compensar. Por lo tanto, la
corrección topográfica se suma siempre, ya se trate de un valle o de una montaña,
ya que éstos generan cierta atracción a la estación tomada.

Figura 20. Corrección por terreno. (Burger, 1992)

35
𝛿𝐸 = 𝛿𝑔𝐹 − 𝛿𝑔𝐵 (Ec 3.18)
𝛿𝑔𝐹 = 0.3086ℎ y 𝛿𝑔𝐵 = 0.0419𝑝ℎ (Ec 3.19)
𝛿𝐸 = (0.3086 − 0.0419)𝑚𝐺𝑎𝑙/𝑚 (Ec 3.20)
Donde p es el promedio de la densidad de la roca en el estudio.
3..3.13Corrección isostática

Las anomalías de Bouguer son casi siempre negativas en zonas montañosas y

altiplanicies y positivas sobre los mares profundos. Además, la magnitud de la
anomalía de Bouguer depende de la altitud del continente o de la profundidad de
las aguas. Tales observaciones sugieren que las correcciones de Bouguer son
excesivamente grandes. El exceso de materiales entre la superficie y el geoide en
los continentes y el déficit de material que implica la presencia de agua en vez de
rocas en los mares, no tienen efectos gravitacionales tan grandes como los
deducidos a partir de sus masas y dimensiones.
Las montañas no son simplemente masas yacentes sobre la superficie de la Tierra,
sino que se originaron a causa de las temperaturas anormalmente elevadas del
interior de la Tierra, que determinaron la dilatación de los materiales y por
consiguiente una disminución de su densidad.

Figura 21. Esquema de hipótesis de Pratt (A) y Airy (B).

36
3.3.14 Anomalía de Aire Libre

La anomalía de aire libre (Free Air) (Δ𝑔AL), se calcula sustrayendo la corrección

por latitud (Gravedad Teórica) de la gravedad observada y se adiciona la corrección
por elevación (Aire Libre).

𝛥𝑔𝐴𝐿 = 𝑔𝑜𝑏𝑠 −𝑔𝑡𝑒𝑜 +𝐶𝐴𝐿 [𝑚𝐺𝑎𝑙] (Ec 3.21)

Donde
gobs= Gravedad observada
gteo=Gravedad teórica
CAL=Corrección por aire libre

3.3.15 Anomalía de Bouguer

El resultado de la reducción de datos gravimétricos es la Anomalía de Bouguer, la

cual correlacionamos solo con las variaciones de densidades laterales de la corteza
terrestre, información que es de especial interés para aplicaciones geofísicas y
geológicas.
La anomalía de Bouguer es la diferencia entre el valor de la gravedad observada
(gobs), ajustado por la suma algebraica de todas las correcciones necesarias (Σcorr)
y restándole la medida hecha en la base (gbase). La variación de la anomalía de
Bouguer nos debe reflejar la variación lateral de la densidad tales como una alta
densidad cambiando a bajas densidades lo que debe de producirnos un incremento
positivo en la anomalía. Por el contrario, si tenemos una variación de bajas
densidades hasta llegar a altas densidades nos da como resultado una anomalía
negativa.

𝛥𝑔𝐵𝑠 = 𝛥𝑔𝐴𝐿 − 𝐶𝐵𝑆 [mGal] (Ec 3.22)

Donde:
𝛥𝑔𝐴𝐿 = 𝐴𝑛𝑜𝑚𝑎𝑙𝑖𝑎 𝑑𝑒 𝑎𝑖𝑟𝑒 𝑙𝑖𝑏𝑟𝑒
𝛥𝐶𝐵𝑆 = 𝐶𝑜𝑟𝑟𝑒𝑐𝑐𝑖ó𝑛 𝐵𝑜𝑢𝑔𝑒𝑟

Para el cálculo de la anomalía de Bouger simple (𝛥𝑔𝐵𝑠), además de tener en cuenta

la elevación de cada estación gravimétrica se consideró la densidad del material,
p=2.67 gr/cm3.

37
Por su parte, la anomalía de Bouger completa (𝛥𝑔𝐵𝐶 ) compensa el excedente o
déficit de masa no contemplado por la corrección de Bouger simple.

Δg BC = Δg Bs + CT + CC (Ec 3.23)

Donde:
𝛥𝑔𝐵𝑠 = 𝐴𝑛𝑜𝑚𝑙𝑎𝑖𝑎 𝐵𝑜𝑢𝑔𝑒𝑟
𝐶𝑇 = 𝐶𝑜𝑟𝑟𝑒𝑐𝑐𝑖ó𝑛 𝑝𝑜𝑟 𝑇𝑒𝑟𝑟𝑒𝑛𝑜
𝐶𝐶 = 𝐶𝑜𝑟𝑟𝑒𝑐𝑐𝑖ó𝑛 𝑝𝑜𝑟 𝑐𝑢𝑟𝑣𝑎𝑡𝑢𝑟𝑎

3.4 Método electromagnético

Este capítulo abraca el entendimiento de las bases teóricas que conforman la teoría
electromagnética, la propagación de estas ondas y su interacción con medios
físicos; así como, también el método electromagnético como herramienta de
exploración geofísica.
El Georradar (Ground Penetrating Radar, GPR)
El radar de penetración terrestre o GPR por sus siglas en ingles es una técnica
geofísica de carácter no destructiva y no invasiva para la exploración del subsuelo,
para su caracterización y monitoreo (Olohoeft, 2006).

3.4.1 GPR

El radar de penetración terrestre o mejor conocido por sus siglas en inglés, ground
penetrating radar (GPR) es una técnica de exploración geofísica que emplea un

38
rango de señales electromagnéticas diseñadas principalmente para la localización
de objetos o estructuras enterradas a una profundidad somera con respecto de la
superficie terrestre.

El rango de aplicaciones del método de GPR es muy amplio y crece de gran manera
gracias al incremento de la adquisición tecnología hoy en día, esto ayuda a la
optimización en reconstrucción de señales; así como, en el diseño de softwares más
sofisticados.

Sus aplicaciones son numerosas, algunos ejemplos son: caracterización de

elementos de construcción (cimentaciones, hormigones, asfaltos, etcétera),
prospección minera, cartografía, estratigrafía, arqueología, paleontología, tectónica
y sismología, contaminación, detección de objetos enterrados, entre muchas otras.

Esta versatilidad del GPR se basa en que ofrece un registro de alta resolución, es
un método rápido y la utilización de intervalos de disparo pequeños permite
considerarlo como un método casi continuo (Pueyo, 2016) (Zamora, 2005).

El geo radar puede proporcionar información sobre la localización exacta (posición

y profundidad) del objeto detectado. Esta técnica, presenta como ventaja, además,
su mayor poder de resolución y la posibilidad de detectar objetos no magnéticos
(Busquet, 1996) (Cukavac, 2008) La respuesta del medio a las señales emitidas,
captada en tiempos determinados, se denomina traza, la cual contiene las
características electromagnéticas del medio atravesado. La sucesiva acumulación
de estas trazas genera lo que se ha denominado Radargramas, la cual, para ser
visualizada, se suele presentar en escalas de colores, generando las imágenes.

3.4.2 Fundamentos Teóricos del GPR

La teoría clásica del electromagnetismo conceptualiza el estudio de la relación entre

fenómenos eléctricos y magnéticos generados por cargas eléctricas en movimiento
o reposo; así como, la interacción entre cada uno de sus campos con cargas y
corrientes.
Estructura de la Teoría electromagnética:
La teoría electromagnética (EM) está conformada por una serie de leyes
fundamentales que rigen el comportamiento de los campos eléctricos y magnéticos
inducidos por cargas eléctricas estáticas y móviles.

3.4.3 Ecuaciones de Maxwell

39
Las ecuaciones de Maxwell son la base teórica en la que se sustenta el fenómeno,
en conjunto con las denominadas ecuaciones constitutivas, que relacionan la
intensidad del campo eléctrico y magnético con el desplazamiento eléctrico y la
inducción magnética, respectivamente, permitiendo correlacionar las interacciones
entre las ondas electromagnéticas y los medios en los cuales se propagan. A
continuación, se presentan las ecuaciones de Maxwell.

∇𝐷 =𝜌𝑓 (Ec 3.24) ∇𝐵 =0 (Ec

3.26)

∇𝐸 =−𝜕𝐵 /𝜕𝑡+𝑀 (Ec 3.25) ∇𝐻 =−𝜕𝐷 /𝜕𝑡+𝐽 (Ec 3.27)

Donde E y H representan respectivamente el campo eléctrico y magnético; D y B el

vector desplazamiento eléctrico y vector inducción magnética; ρ la densidad de
carga libre y J la densidad de corriente. Este sistema de ecuaciones en conjunto
con las ecuaciones de continuidad, las cuales se presentan a continuación, permite
el entendimiento de las interacciones entre las ondas electromagnéticas y los
medios en los cuales se propagan.
⃗ = 𝜀𝐸⃗ (Ec 3.28)
𝐷
⃗ = 𝐵⃗ (Ec 3.29)
𝐻 𝜇
𝐽 = 𝜎𝐸⃗ (Ec 3.30)

Este sistema de ecuaciones contiene los tres parámetros que caracterizan un medio
electromagnéticamente: ε, σ, μ. La conductividad σ representa la medida de la
capacidad de un material de conducir corriente eléctrica. La permitividad dieléctrica
ε brinda la medida de la capacidad de polarización de un material en presencia de
un campo eléctrico. Finalmente, la medida de la capacidad de un medio para atraer
los campos magnéticos viene descrita por la permeabilidad magnética μ.

3.4.4 Ecuación del GPR

La ecuación del georadar para equipos monoestáticos (antena en función emisor-

receptor), presentada en la literatura está determinada por la relación (Q) entre la

40
potencia recibida, PR, y la potencia emitida, PE. La Ecuación 2.8, representa esta
relación.

(Ec 3.31)
Donde:
𝜆: Longitud de onda en aire.
G: Ganancia de la antena.
S: Superficie efectiva del reflector.
𝛼: Constante de atenuación.
r: Distancia al reflector.

De la Ecuación 2.8, se puede obtener la relación entre amplitudes de onda (QAmp)

recibida, AR, y emitida, AE, de la siguiente manera (Identification of cavities by
extraction of characteristic parameter from ground probing radar reflection data,
1997):

(Ec 3.32)
3.4.5 Parámetros Electromagnéticos

Los parámetros electromagnéticos (conductividad, permitividad dieléctrica y

permeabilidad), son los que definen un medio al paso de una onda
electromagnética.
La conductividad de un medio proporciona la medida de la respuesta de las cargas
libres en él, en presencia de un campo eléctrico externo, siendo el factor de
proporcionalidad entre el campo libre aplicado y la densidad de volumen de corriente
debido al movimiento de estas cargas libres. Es decir, proporciona la medida de la
capacidad de un material de conducir corriente eléctrica y es descrita según la ley
de Ohm, cuya formulación fue presentada en la Ecuación 2.7.
Las unidades de medida de la conductividad, en el S.I., son omhs/metro o el
equivalente de Siemens/metro. La conductividad de un medio es la inversa de su

41
resistividad. En general se puede distinguir entre materiales conductores,
semiconductores y aislantes.

3.4.6 Permitividad Eléctrica

La permitividad es una constante que da medida de la capacidad de polarización de

un material en presencia de un campo eléctrico. Proporciona un valor de la
respuesta estática del material cuando está en presencia de un campo eléctrico
externo. La constante dieléctrica del vacío es 1. Para la mayoría de los materiales
se pueden encontrar en el subsuelo al realizar prospección electromagnética, los
valores de permitividad relativa se encontrarán entre 1 (aire) y 81, siendo está ultima
la constante dieléctrica del agua a 20° de temperatura (Vega, 2002).

3.4.7 Velocidad de Propagación y Longitud de Onda

Las ondas electromagnéticas viajan en un medio a una velocidad que depende de

la constante dieléctrica relativa del subsuelo, de la velocidad de la luz en el vacío, y
de la permeabilidad magnética relativa. La velocidad de la onda electromagnética
(Vm) en un cierto material viene dada por: La conductividad de un medio
proporciona la medida de la respuesta de las cargas libres en él, en presencia de
un campo eléctrico externo, siendo el factor de proporcionalidad entre el campo libre
aplicado y la densidad de volumen de corriente debido al movimiento de estas
cargas libres. Es decir, proporciona la medida de la capacidad de un material de
conducir corriente eléctrica y es descrita según la ley de Ohm, cuya formulación fue
presentada en la Ecuación 2.7.
Las unidades de medida de la conductividad, en el S.I., son omhs/metro o el
equivalente de Siemens/metro. La conductividad de un medio es la inversa de su
resistividad. En general se puede distinguir entre materiales conductores,
semiconductores y aislantes.

3.4.9 Velocidad de Propagación y Longitud de Onda

Las ondas electromagnéticas viajan en un medio a una velocidad que depende de

la constante dieléctrica relativa del subsuelo, de la velocidad de la luz en el vacío, y
de la permeabilidad magnética relativa. La velocidad de la onda electromagnética
(Vm) en un cierto material viene dada por:

(Ec 3.33)

42
donde P es el factor de pérdida, de manera que P = σ / ωε y σ es la conductividad,
ω = 2πf ( f es la frecuencia), y ε la constante dieléctrica absoluta. En un material con
un bajo factor de pérdida tenemos P → 0; entonces la velocidad de la onda
electromagnética es dada por:

(Ec 3.34)

Tabla 3. Velocidades de onda materiales y rocas.

Roca/Material Velocidad
m/ns
Aire 0.3
Hielo 0..16-0.17
Suelo Seco 0.15
Arena Seco 0.15
Granito 0.13
Caliza 0.12
Concreto 0.08-0.12
Pavimento 0.10
Arenisca 0.09
Arcilla 0.06
Agua dulce 0.033
Agua salada 0.033

La profundidad de penetración (D) puede ser determinada por la Ecuación 2.12, a

partir de la velocidad en el medio Vm. Finalmente, el tiempo doble de propagación
puede ser determinado mediante representación gráfica de las señales del GPR.
Con estas premisas la profundidad resulta ser (Tavera, 2006)

(Ec 3.35)

43
3.4.10 Procesamiento de Información de GPR

Los registros obtenidos con las prospecciones realizadas con georadar pueden ser
interpretados, en algunas ocasiones, de manera directa. Sin embargo, dada la
inmensa variabilidad de los medios y sus constituyentes (objetos enterrados tales
como tuberías, entre otros), en la mayoría de las ocasiones, se hace necesario tratar
los datos, con el objetivo de mejorar la relación señal/ruido o para resaltar algún
evento. Se debe mencionar que los mismos deben ser empleados con especial
cuidado, ya que, aunque puedan contribuir a mejorar el registro, una incorrecta o
excesiva aplicación, puede conllevar a la pérdida de información importante o llevar
a interpretación y conclusiones incorrectas. A continuación, se presentan algunos
de los más relevantes procedimientos aplicados para este fin.

3.4.11 Filtros

El término filtrado, hace referencia a una variedad amplia de técnicas cuyo objetivo
es modificar selectivamente la amplitud o la frecuencia de las señales recibidas.
Este tipo de procedimiento podría clasificarse en dos grupos: filtros horizontales y
verticales.

3.4.12 Filtros Horizontales

Consiste en un tratamiento que se aplica horizontalmente en el registro, en dirección

espacial. El objetivo perseguido con la aplicación de estos filtros es la eliminación
del ruido de fondo de un radargrama. Este ruido que puede ser de alta frecuencia
(se registra en pocas trazas seguidas), generado al deslizar la antena por un terreno
irregular.
Esta metodología es una de las más aplicadas en el tratamiento de datos en
registros de georadar. El hablar de contenidos frecuenciales, para este tipo de de
procesado, está referido a la longitud horizontal de las anomalías. Una anomalía de
gran longitud de onda (y por ello de baja frecuencia) quedará registrada en un
elevado número de trazas, mientras que una anomalía de pequeña longitud de onda
(alta frecuencia) lo será en pocas trazas. Este tipo de tratamiento puede realizarse
con dos tipos de filtros: los denominados filtros de respuesta finita (FIR) y los
denominados filtros de respuesta infinita (IIR).

3.4.13 Filtros Verticales

44
Son tratamientos que se aplican en la dirección vertical o temporal de los registros,
sobre cada una de las trazas. Pueden utilizarse para realizar correcciones de línea
base o para eliminar parte del ruido electrónico inherente a cada traza. Consisten
en multiplicar cada señal por un operador que vale la unidad en la banda de
frecuencias que interesa conservar, y que tiende a cero en la banda de frecuencias
que no interesa. Cuando se aplica estos filtros es conveniente realizar estudios de
las frecuencias antes de la aplicación y después de la misma, para comprobar que
se ha realizado correctamente y para evitar perder información que podría ser
relevante en el estudio realizado. Pueden ser también de pasa altas o pasa bajas;
los filtros de paso alto, impiden que se registren bajas frecuencias, normalmente
aquellas que se encuentran por debajo de 10 MHz, las cuales son consideradas
como ruidos del sistema. Por lo tanto, tienen como objetivo eliminar el ruido de baja
frecuencia, que depende en gran parte de la antena que se utiliza, del cableado y
de su longitud. Por otra parte, de manera opuesta los filtros pasan bajas, se utilizan
para eliminar frecuencias más elevadas que las esperadas. Con la aplicación de
estos dos filtros queda demarcado el rango de frecuencia del estudio.

3.4.14 Deconvolución

Se trata de un procedimiento que tiene como objetivo la mejora de la resolución

vertical en los radargramas. Es una técnica que se ha aplicado a la prospección con
georadar y que es proveniente de los estudios de sísmica. Se aplica porque se
considera que el medio por el que se produce la propagación de la señal actúa de
filtro sobre la misma, de manera que la traza grabada en el radargrama es en
realidad una convolución del pulso inicial emitido. El objetivo que persigue es
deconvolucionar la onda para eliminar los efectos del medio sobre el pulso inicial.
De esta manera, se persigue transformar el radargrama de campo en un registro
donde las señales representan únicamente la reflexión primaria.

3.4.15 Migración

Este procedimiento se realiza para trasladar a su posición real los efectos de la

inclinación de las reflexiones registradas y colapsar las difracciones que producen.
Para utilizar este tratamiento se necesita estimar la velocidad promedio de
propagación por el medio, pero puede conocerse a partir de las hipérbolas
generadas en las difracciones. Consiste básicamente en el cálculo de la ecuación
de cada una de las hipérbolas, desplazando elreflector al vértice de la hipérbola.
Como la ecuación de la hipérbola depende de la velocidad de propagación del
medio, un error en la determinación de este parámetro puede ocasionar que no
llegue a colapsar totalmente dicha hipérbola (velocidad estimada inferior a la real) o
bien que se produzca el efecto conocido como sobremigrado (la velocidad estimada
es superior a la real).

45
3.4.16 Transformada Rápida de Fourier (TDF)

Este tratamiento de datos, conocido ampliamente, consiste en la realización de un

análisis frecuencial de los registros, lo que permite el estudio de la evolución de los
parámetros de la señal en función de la frecuencia. A menudo se utiliza para obtener
los espectros de frecuencia de las trazas, con el fin de seleccionar mejor los filtros
frecuenciales que se pueden aplicar, o para comprobar los efectos que sobre las
frecuencias ha tenido la aplicación de alguno de los mencionados filtros.

3.4.17 Transformada Hilbert

Este procedimiento matemático expresa la relación entre la parte real y la parte

imaginaria de las señales, con el que pueden obtenerse envolventes de las
amplitudes en los registros, la fase y la frecuencia instantáneas. Estas últimas
pueden utilizarse como indicadores del tipo de filtrado que sufre la señal durante su
propagación por el medio, proporcionando información del comportamiento del
medio frente a la señal electromagnética aplicada. La obtención de la envolvente de
una traza proporciona información acerca de la amplitud instantánea de la señal y,
por tanto, acerca de la energía implicada en cada evento.

46
4 Adquisición de Datos

De acuerdo con nuestro objetivo dentro del área de estudio consiste en encontrar
vacíos en el subsuelo posteriormente estos serían relacionados con posibles
túneles dentro del pueblo de Cerro de San Pedro, para ello se hizo uso de métodos
geofísicos de baja penetración en el caso del método del radar de penetración
terrestre y de métodos potenciales (gravimetría y magnetometría) esto debido a que
la ventana de ubicación de un túnel es comprendida de forma somera con respecto
a otros estudios de mayor interés de penetración.

4.1 Levantamiento Gravimétrico y Magnético

Los métodos potenciales nos dan información desde cientos de km hasta los
primeros metros del subsuelo esto dependen según sea nuestro objeto de estudio,
por ello. Durante un estudio gravimétrico el espaciamiento entre los puntos de
medición puede variar desde unos cuantos metros, si el caso requiere detalle, como
en estudios arqueológicos, geotécnicos o mineros; hasta varios kilómetros, como
en el estudio de reconocimientos regionales (Milsom, 2003).
Cuando se trata de estudios donde el espaciamiento entre cada estación se reduce
a tan solo unos cuantos metros es considerado como microgravimetría (Fig. 22).
La microgravimetría es una técnica que está basada en la medición de diminutas
variaciones que existen debido a la acción gravitatoria de la Tierra correspondientes
a cambios de distintas densidades de los materiales que constituyen el subsuelo.
Este método propone hacer uso de la variación de la gravedad debido a las
alteraciones de la densidad de las rocas en función de su contenido.
Dentro de las principales aplicaciones de esta técnica se encuentra la detección de
vacios, debido a que esto supone un considerable déficit de masa en el subsuelo y
genera considerables contrastes de densidad con respecto al área circundante al
terreno. Contrastes de alta densidad de igual manera pueden ser registrados y estos
son asociados a altas concentraciones minerales o zonas de alta permeabilidad
p.e., las de baja concentración son relacionadas a zonas de fractura o zonas de
alteración p.e.

47
Figura 22.- Levantamiento microgravimetría.
Para este caso de estudio que fue realizado dentro del municipio de Cerro de San
Pedro, SLP el espaciamiento entre estación de medición fue de 5 metros se
recolectó información de la medición de gravedad, magnética, así como la
información topográfica de cada punto. Fueron integrados un total de 214 puntos, la
distribución de estos puntos fue un tato uniforme ya que estos fueron tomados a lo
largo de carreteras y las calles principales con acceso para este lugar.

4.1.1 Establecimiento de estación base

Para el cálculo de valores de gravedad absolutos se requiere cuando menos de una
estación base que se encuentre cerca de la zona de estudio, la estación base por
lo general es de fácil acceso y además se conoce el valor absoluto de la gravedad,
para los casos en los cuales no se cuente con el acceso de una estación base, se
puede establecer una, midiendo la diferencia en gravedad entre la estación base y
la base local.

48
4.2 Procedimiento

Se realizó la adquisición de datos en una primera campaña para el levantamiento

magnético el cual fue realizado el 17-11-2017, en donde se estableció una base
local a las afueras de la entrada principal al pueblo. Para el levantamiento
gravimétrico y topográfico fue realizado en tres campañas llevadas a cabo 6-12-
2016, 28-01-2017 y [Link] lecturas gravimétricas se obtuvieron a partir de
la diferencia de aceleración gravitacional entre la estación base de gravedad
absoluta ubicada en Morales, SLP. A partir de ahí se estableció una base de
gravedad absoluta para el área de estudio. De forma simultánea para el
levantamiento gravimétrico se realizó el estudio topográfico de cada estación
mediante el uso de un GPS diferencial.

4.3 Instrumentación

4.3.2 Gravimétrico

Para el levantamiento gravimétrico se implementó el uso de gravímetro SCINTREX

Autograv, modelo CG-5 (Fig. 23), el cual tiene una resolución estándar de 1
microGal con una desviación estándar inferior a 5 microGals. El elemento de
detección del Autograv se basa en un sistema elástico de cuarzo fundido con
anulación electrostática (SCINTREX, 2006). La fuerza gravitacional sobre la masa
de prueba es balanceada por un resorte y una fuerza de restauración electrostática
relativamente pequeña. La posición de la masa, la cual es detectada por un
transductor de desplazamiento capacitivo, se altera con el cambio en gravedad. Un
circuito de retroalimentación automática aplica voltaje DC a las placas del capacitor
produciendo una fuerza electrostática sobre la masa que lo lleva de nuevo a una
posición estable. El voltaje de retroalimentación, que es una medida del valor
relativo de gravedad en el lugar de la lectura, se convierte a una señal digital y luego
se transmiten al sistema de adquisición de datos del instrumento para el
procesamiento, visualización y almacenamiento.

49
Figura 23. Gravímetro modelo CG-5 (a) vista lateral (b) vista de planta.

Para la operación del gravímetro, puede ser de forma remota mediante el uso de un
control o de forma manual, para este estudio se manipuló mediante la primera
opción ya que de esta manera se reduce el ruido al momento del aparato estar
midiendo. Previo a la medición, es nivelado sobre el tripee. Para la medición el
gravímetro es programado para la duración de la medición, el gravímetro puede
hacer mediciones por segundo, debido a que la zona en la que se trabajó es
considerada de baja influencia de ruido este se programó para realizar la lectura de
gravedad en 60 segundos por estación.

Dentro de los datos de campo sólo se utilizaron aquellos en los que el error no
superaba 0.002 mGales de error, si por alguna razón las mediciones se veían
afectadas por ruido y este excedía el parámetro anterior la lectura de la estación era
repetida. Para cada campaña llevada a cabo se realizó una lectura de apertura en
la estación de base y una de cierre, esto para facilitar la corrección por deriva
instrumental.

4.3.2 Magnetometría

Para el levantamiento magnético se utilizó un magnetómetro Gem system

modelo GSM-19 (Fig. 24) el cual opera mediante los siguientes pasos:
a) Polarización: Una fuerte corriente FR es pasada atreves del sensor
creando un flujo rico de protones en el sensor.

50
 b) Desvió: Un pequeño pulso es desvía al protón magnetizado dentro del
plano de precisión.
c) Pausa: La pausa permite que los transitorios eléctricos mueran,
dejando un protón que se descompone lentamente.
d) Cuenta: La frecuencia de precisión de campo magnético por encima
del nivel de ruido.
e) Almacenamiento: Los resultados se almacenan en la memoria junto
con la fecha, la hora y las coordenadas de la lectura, cuando es usado
como base solo el tiempo total y la lectura del campo son
almacenados.

Figura 24. Magnetómetro modelo GSM-19.

Previo al levantamiento se ubicó la estación base para realizar correcciones por

deriva instrumental del aparato para la mejora de la calidad de los datos. El
magnetómetro fue operado en modo móvil, en esta opción nos permite tomar la
medición de forma manual, la duración por estación fue programada para 30
segundos por estación mientras que en la estación base se programó para que
realizara una lectura cada 120 segundos.
El levantamiento magnético al igual que el gravimétrico fueron realizados a lo largo
de calles y vías de acceso del pueblo de Cerro de San Pedro, un tanto irregular,
pero se trató de cubrir la mayor área posible del área de estudio (Fig. 25).
Los datos originales utilizados fueron aquellos que no superaban el 5% de error, se
debe tomar en cuenta elementos como autos, cableado eléctrico es considerado
como ruido, las lecturas en lugares cercanos a estos objetos fueron registrados en
campo para a la hora de procesar tomarlos en cuenta.

51
Figura 25. Levantamiento magnetométrico terrestre en Cerro de San Pedro, SLP. En la figura A se puede observar cableado
eléctrico que debe ser tomado en cuenta a la hora del procesamiento debido a que provoca ruido de la señal. En la figura B
se observa una zona libre de cableado donde existe poca influencia de ruido.

4.4 Levantamiento topográfico

52
El levantamiento topográfico fue realizado mediante el uso de un GPS de tipo
diferencial de la marca TRIMBLE, se designó una estación de monitoreo para la
antena tipo base (Fig. 26) y mediante la otra antena móvil se realizó el levantamiento
de cada estación para recaudar la información de su posición (componente x, y, z).
Fue programado usando el modo Estático rápido.

Ubicación base monitoreo.

Debido a que no se contaba con un banco de nivel cercano para nuestra base de
monitoreo se designó una, y para poder conocer su posición se procesaron los datos
de rinex de 4 estaciones de monitoreo cercanas (Aguascalientes, Zacatecas,
Guanajuato y Querétaro) mediante este procesado se georreferenció nuestro punto
de monitoreo y atreves del vectorizar cada estación en la que se tomaron lecturas
de gravedad y magnética.

Figura 26. Levantamiento topográfico utilizando GPS diferencial.

53
4.5 Levantamiento con radar de penetración

El radar de penetración terrestre o mejor conocido por sus siglas en inglés, ground
penetrating radar (GPR) es una técnica de exploración geofísica que emplea un
rango de señales electromagnéticas diseñadas principalmente para la localización
de objetos o estructuras enterradas a una profundidad somera con respecto de la
superficie terrestre (Fig.27).

El rango de aplicaciones del método de GPR es muy amplio y crece de gran manera
gracias al incremento de la adquisición tecnología hoy en día, esto ayuda a la
optimización en reconstrucción de señales, así como en el diseño de softwares más
sofisticados.

Figura 27. Esquema del funcionamiento del GPR.

A continuación, se enlistan algunas de las principales aplicaciones del método:

54
• Investigaciones Arqueológicas
• Análisis de puentes
• Inspección de Yacimientos Petroleros
• Evaluación de estados de Construcción
• Investigación de suelos contaminados
• detección de minas enterradas
• evaluación de hormigón armado
• investigaciones forenses
• investigaciones geofísicas
• Detección de tuberías y cableado

4.5.1 Procedimiento

La adquisición de líneas de Georadar se realizó a lo largo de las calles y vías de

acceso al igual que el levantamiento de métodos potenciales, en un principio se
recogieron datos con una antena de 800 Hz, al momento de visualizar los datos se
observó que la penetración no fue la necesaria para diferenciar elementos posibles
presentes. Por lo que se optó por usar una antena de 100 Hz con la cual se
consiguió un alcance de +/- 20 metros (Fig. 28).

4.5.2 Instrumentación

Se utilizó el radar de penetración de la marca MALA, para la localización de estas

líneas se tomaron puntos de inicio y fin de cada línea capturada empleando el GPS
diferencial Trimble.

55
Figura 28. Levantamiento líneas de GPR usando una antena de 100Hz.

El principio físico del radar de penetración está basado en la emisión de una señal
de impulsos electromagnéticos para la detección y localización de elementos,
formaciones o anomalías presentes en el subsuelo a nivel somero empleando el
principio de reflexión de ondas electromagnéticas propagadas en un medio. (de
forma análoga se compara con el método sísmico)
Las emisiones de impulsos electromagnéticos son de muy corta duración (entre 0.6
ns y 10 ns) en la banda de frecuencias UHF y VHF que varían de 100Mhz y 1.5
GHz.
Las lecturas tomadas por el GPR se efectúan esquemáticamente siguiendo el
siguiente orden:
1- Mediante el uso de la unidad central de mando se programa el pulso
electromagnético generado por medio de la antena emisora.
2- La onda electromagnética se propaga por el medio, hasta chocar con
elementos de diferentes propiedades electromagnéticas respecto al material
circundante.
3- Una parte de la señal de la onda se refleja hacia la superficie, mientras que
parte de la energía sigue propagándose, generándose múltiples reflexiones.
4- La onda reflejada es almacenada en la unidad central para su posterior
tratamiento e interpretación.

56
5 Resultados
5.1 Levantamiento magnético
Después de las correcciones correspondientes aplicados a los datos obtenidos en
campo, se obtuvo el campo magnético total (CMT). Posteriormente a estos datos se
les aplicó la resta del campo magnético producido por la tierra (IGRF), al resultado
de esta operación se le conoce como Campo Residual (Fig. 29). Directamente con
esta información se pueden correlacionar con estructuras o información según el
tipo de material que esté presente en el área de estudio, teniendo así una primera
referencia de la zona de estudio, en la que destacan valores altos los cuales se
pueden apreciar que se encuentran dispuestos en diagonal desde el SW al NE. Se
puede observar una anomalía de carácter negativa de forma alargada localizada en
el segundo cuadrante dispuesta de norte a sur, lo cual se puede atribuir a una
dipolaridad entre esta estructura y la anterior donde los valores oscilan entre los -
80nT y 140nT indicados en el mapa.

Figura 29. Mapa del Campo magnético Residual de Cerro de San Pedro,SLP

57
5.2 Reducción al polo

Al mapa obtenido del CMR (Fig. 30), se le aplicó la reducción al polo, dicho proceso
sitúa en su posición verdadera a las anomalías presentes, esto debido a que las
fuentes difieren de la posición de las anomalías. Al aplicar el filtro las anomalías
magnéticas se sitúan sobre las fuentes que las causan, simulando del traslado de
la zona de trabajo al polo Norte en donde la inclinación del campo magnético es
vertical.

Figura 30. Mapa Campo magnético reducido al polo Cerro de San Pedro, SLP.

58
5.3 Dominios magnéticos DM

5.3.1 DM I
Este dominio magnético es caracterizado por presentar los valores más bajos que
van desde los -80 hasta -20 nT, presenta una forma alargada de Norte a Sur,
localizado en la parte SW del mapa (parte superior izquierda) y es atribuido a la
zona de remoción producto de la minera, como también a material aluvial o poco
consolidado.

Figura 31. Mapa del Campo magnético reducido al polo Cerro de San Pedro, SLP

5.3.2 DM II
Dominio magnético presente en la parte central del mapa, de gran extensión y con
valores que oscilan entre los -30 nT a 60 nT. Este dominio es asociado a la
presencia de rocas de tipo sedimentario, Calizas y Lutitas.

59
5.3.3. DM III
Este dominio magnético está formado por dos altos magnéticos ubicados al SW y
SE del mapa (inferior derecha e izquierda respectivamente) con valores que varían
entre los 80 nT y 120 nT, son caracterizados por ser los valores más altos y son
atribuidos a la presencia de rocas de tipo volcánicas, Dioritas y riolitas.

5.4 Primera derivada en Z

A partir del mapa de Reducción al polo se le aplicó el filtro de primera derivada en

Z (1aDz) lo cual permite identificar lineamientos magnéticos(Fig. 32). Un lineamiento
constituye una flexión o rasgo físico identificable en la superficie terrestre o en su
caso en el subsuelo, a través de los patrones de anomalías geofísicas, el cual, como
su nombre lo indica, tiene una forma lineal, rectilínea o suavemente curvilínea
(O’Learly et al., 1976). Los lineamientos se correlacionan a zonas de alto gradiente
que pueden estar asociados a fallas geológicas, fracturas y/o contactos litológicos.

60
Figura 32. Primera derivada en “Z” Cerro de San Pedro, SLP

Se identificaron los lineamientos presentes para el mapa de primera derivada en Z

donde se pueden identificar preferencialmente la dirección de las estructuras
presentes, las cuales se ubican mayormente con la dirección 15 a 25 NE (Fig. 33).
Algunas de las anomalías son cruzadas por algún lineamiento lo cual puede ser
inferido por una falla o fracturas, todos estos lineamientos fueron de carácter
superficial.

61
Figura 33. Principales lineamientos marcados con base en el mapa de primera derivada en”Z”.

5.5 Señal analítica

Al mapa de reducción al polo se le aplicó la señal analítica la cual permite resaltar

los bordes de las anomalías y así delimitar las posibles estructuras a las cuales
estas están asociadas (Fig. 34).

62
Figura [Link] señal analítica Cerro de San Pedro, SLP

Se marcaron Tres dominios a partir de la señal analítica en donde estos dominios

delimitan las estructuras presentes de la zona de estudio debido a la relación
topográfica y geológica presente en la zona (Fig. 35).

63
Figura 35. Dominios magnéticos con base en el mapa de señal analítica de Cerro de San Pedro, SLP.

64
5.5.1 DMSA I
Este dominio ubicado en la zona central de nuestro mapa con una marca de tipo “U”
es asociado a la depresión existente en la zona de estudio la cual corresponde a un
caudal.

5.5.2 DMSA II
Corresponde a la parte de valores medios y es asociado a rocas sedimentarias
como calizas y Lutitas además corresponde a la zona de media altura de la zona de
estudio.

5.5.3 DMSA III

Este dominio es caracterizado por valores altos magnéticos y al cual le
corresponden las zonas altas presentes a la topografía del área de estudio.
Se efectuó la aplicación del filtro derivada en la componente “x” (Fig. 36), donde se
pudo observar una serie de lineamientos magnéticos que contrastan y se pueden
relacionar de manera directa con las fracturas presentes a un mismo nivel de
profundidad como se encontró más adelante con la información proporcionada por
los radargramas usando una antena de 100 MHz.

5.6 Análisis Gravimétrico

65
5.6.1 Análisis Regional

Para realizar el análisis regional se utilizó el modelo de elevación digital obtenido

dentro del levantamiento topográfico descrito anteriormente (Fig. 37a), se
compararon las anomalías gravimétricas obtenidas con las respuestas topográficas,
llegando a la siguiente relación:
Los mínimos gravimétricos que oscilan los 220mGal mostrados con el número uno
corresponde a los sitios donde el nivel topográfico es menor dentro de la zona de
estudio, mientras que los valores de nuestras anomalías gravimétricas más altas se
les atribuye los puntos de mayor elevación de esta zona y en la cual predominan
valores de aproximadamente 230 mGals (Fig. 37b).

Figura 36. Figura A Mapa de elevación topográfica. Figura B Mapa de anomalía de Bouger completa.

A la anomalía de Bouger se le aplicó la separación de sus dos componentes,

anomalía residual y regional (Reynolds, 2007), esta primera se utiliza para analizar
de forma somera las fuentes las cuales pueden estar relacionadas (Fig. 38a), y la
segunda corresponde a las anomalías de mayor influencia atribuidas a cuerpos a
mayor profundidad (Fig. 38b).

66
Se obtuvo la anomalía regional y residual aplicando un filtro Gaussiano de regresión
polinomial de tercer grado.

Se puede observar en la siguiente imagen donde se hace una comparativa entre la

Anomalía de Bouger (AB) y la anomalía Regional (AR) calculada, donde se puede
inferir que la mayor influencia de las anomalías presentes son debido a estructuras
de gran profundidad, además se puede observar como el mapa de la AB y AR son
sumamente parecidas, la interpretación a la cual se puede llegar es que las
anomalías positivas corresponden a los puntos más elevado así como los puntos
de menor elevación están relacionados con las anomalías de carácter negativo.

También se puede hacer una interpretación donde se observa que en la parte centro
de nuestro mapa existen una tendencia homogénea de anomalías intermedias, lo
cual se atribuye litológicamente a rocas sedimentarías presentes en la zona que
corresponden a calizas y lutitas. En la parte SW del mapa (esquina inferior derecha)
se observa un bajo gravimétrico el cual topográficamente correlaciona con el mde,
aunado a esto esta anomalía negativa es formada por la presencia de material
cuaternario (aluvión y coluvión).

Se observa en la parte norte y cortando de NE a SW la presencia de anomalías

positivas que encajan de forma acertada con el mde, se le atribuye a esta respuesta
como rocas volcánicas como riolita o latita.

67
Figura 37. Contraste entre el mapa de Anomalía de Bouger Completa (A) y la anomalía Regional calculada (B).

5.6.2 Anomalías residuales

Las anomalías residuales son atribuidas a elementos a que se encuentran a una

escala somera respectivamente con respecto a las anomalías regionales, debido a
que el estudio se remonta a una escala de menor penetración se hará un mayor
enfoque a estas últimas.

Comparando el mapa de Anomalía residual (Ares) (Fig. 39a) y el modelo de

elevación (mde) (Fig. 39b) se corrobora aquellas mismas estructuras que
corresponden a puntos elevado y aquellos lugares de mayor depresión, pero se
puede distinguir también aquellas anomalías que no corresponden con el mde, es
donde se trata de atribuir el tipo de estructura o elemento litológico al cual
corresponde.

68
Figura 38. Figura relación entre el MDE (A) y la anomalía Residual (B).

Se pueden distinguir en el mapa de Ares una anomalía de tipo dipolar ubicada en la

parte centro de nuestra zona de estudio, indicadas con dos círculos enumerados
como 1 y 2 (Fig. 40a), estas anomalías se pueden observar que encajan de manera
correcto con el mde, tanto la anomalía positiva con valores promedio de 0.60 mGal
que corresponde a una zona de alta elevación como la anomalía negativa con
valores intermedios de -0.74 mGal (Fig. 40b) correspondientemente.

Se ilustran dos pentágonos enumerados como 3 y 4 respectivamente, los cuales

corresponden a anomalías presentes que no corresponden con el mde, para la
anomalía número 3 se puede describir como una anomalía de tipo positiva por lo
cual se esperaría tener un punto alto dentro del mde y se encuentra lo contrario, a
este valor positivo se puede correlacionar con la presencia de rocas de tipo
volcánicas como lo son las dioritas o directamente al pórfido que intrusiona a las
rocas calizas presentes dentro de nuestra zona de estudio, es generada una
anomalía positiva debido a la diferencia contrastante de la densidad relativa de este
tipo de rocas volcánicas.

69
Mientras que en la anomalía representada con el número cuatro se observa que
está localizada en una zona de alto relieve local como lo manifiesta nuestro mde,
pero esta presenta una anomalía negativa, podemos correlacionar esta anomalía
con la geología presente de nuestra zona, y se puede llegar a inferir como la
presencia de rocas de tipo sedimentarias que corresponden a Calizas y Lutitas, las
cuales se encuentran en alternancia dentro de nuestra zona de estudio.

Figura 39. Interpretación del mapa de anomalías residuales (B) con respecto al MDE (A.)

Al mapa de AB se le aplicó la señal analítica, esta para definir de una manera más
clara las estructuras presentes dentro de la zona de estudio (Fig. 41).

70
Figura 40. Señal analítica Cerro de San Pedro, SLP

Al mapa de AB se le aplicó el gradiente horizontal, esta para definir las estructuras

presentes dentro de la zona de estudio (Fig. 42).

71
Figura 41. Gradiente Horizontal, Cerro de San Pedro, SLP

5.7 Resultados GPR

A las líneas levantadas durante la campaña (Fig. 43; Tabla 4), se les aplico una
serie de filtros para eliminar el ruido presente, y también para realzar elementos que
correspondan a estructuras u objetos enterrados en el subsuelo, las imágenes
resultantes arrojan como resultado el contraste de evidente debido a reflectores
superficiales.

A continuación, se presentan las imágenes de los seis perfiles tomados en la zona

de estudio, así como una interpretación de lo que se puede distinguir principalmente
en ellas.

72
5.7.1 Ubicación perfiles GPR

Figura 42. Distribuciones líneas de radar de penetración a lo largo de la cabecera municipal de Cerro de San Pedro, SLP

73
Tabla 4. Geolocalización inicio y fin de puntos de líneas de radar de penetración

Número Coordenadas punto inicio de línea (A) Coordenadas punto fin de línea (B)
de Línea
GPR
1 Latitud: Longitud: Latitud: Longitud:
22.2156 86N -100.801059 W 22.215982N -100.800206 W

2 Latitud: Longitud: Latitud: Longitud:

22.216896N -100.8000060W 22.2171121N -100.800477 W

3 Latitud: Longitud: Latitud: Longitud:

22.217053N -100.800045 W 22.217999N -100.8000168W

4 Latitud: Longitud: Latitud: Longitud:

22.216868N -100.800054 W 22.216856N -100.7999891W

5 Latitud: Longitud: Latitud: Longitud:

22.216302N -100.800591W 22.218057N -100.800350 W

6 Latitud: Longitud: Latitud: Longitud:

22.217823N -100.799963 W 22.217534N -100.79866184 W

5.7.3 Interpretación general de Radargramas

Las principales estructuras que se encontraron y lograron interpretar fueron,

fracturas, tuberías o materiales enterrados, tres tipos de materiales rocosos (Figs.
44, 45, 46 y 47).
En los primeros arribos de nuestros Radargramas se puede distinguir un horizonte
de alrededor de poco más de un metro de espesor, donde todos los reflectores
muestran las mismas propiedades físicas, la baja ganancia expuesta, además como
se trata de los primeros metros del subsuelo, podemos asociarlo con el aluvión
presente en la zona de estudio de manera tendencial se aparece este horizonte en
nuestros seis registros, caracterizado por estas mismas propiedades. (Figs. 44).
Se encuentra siguiendo con la observación de reflectores, con otros dos diferentes
que muestran condiciones en los que se delimitan dos diferentes tipos de roca, el
primero marcado como un reflector delgado y muy bien consolidado, para el
segundo se caracteriza por mostrar un espesor de mayor grado y se observa
menormente consolidado, estos se correlacionan geológicamente como calizas y
lutitas respectivamente, se registra la alternancia de estos dos tipos de rocas,
geológicamente podemos ver la intercalación de estas rocas (Figs. 45).

74
A lo largo de los seis Radargramas se observa el fracturamiento de las dos unidades
rocosas presentes, el grado de fracturamiento presente en los Radargramas es de
alto grado (Figs. 46).
Por ser un asentamiento activo la influencia de redes de drenaje y eléctrico juegan
un papel importante dentro de la interpretación de nuestros Radargramas, es por
eso por lo que se interpretan estas estructuras que se presentan en la mayoría de
nuestros levantamientos, las cuales se observan en los registros como hipérbolas,
la mayoría se identifican en una ventana de 3 a 5 metros, estas alineadas de forma
transversal con respecto a la dirección de levantamiento de nuestras líneas. (Figs.
47)
Se pueden observar hipérbolas que no están asociadas a tuberías, esto debido a
que la ventana en la cual se encuentran es superior a la de 4 metros de profundidad,
se puede inferir que se trata de otro tipo de objetos enterrados, los cuales se
desconocen.
Dentro de este análisis no se encontró la presencia de las representaciones típicas
de vacíos, los cuales se pueden atribuir a túneles, ni la presencia de material de
relleno en el dado caso de que hayan existido estos, colapsado y al final rellenado
con material rocos presente, por este método deja de forma clara la idea de la
inexistencia de túneles dentro de nuestra área de estudio, de forma detallada se da
una interpretación y las condiciones de la estructura somera del área de Cerro de
San Pedro.
5.7.4 Calculo de profundidades a partir de velocidades de propagación

La velocidad utilizada para el cálculo de las profundidades de las anomalías

presentes fue de 0.12 m/ns tomada como un promedio dentro de zona de estudio,
además de que es un valor de velocidad promedio para la roca caliza que varía

75
 entre los 0.9 m/ns y los 15 m/ns.

Los recuadros en color Los recuadros en Los círculos en color

rojo indican hipérbolas verde y en azul indican rojo corresponden a
que son interpretadas dos tipos de reflectores fracturas, muchas de
como la presencia de con propiedades ellas presentes a lo
tuberías dentro de la distintas, lo cual es largo de todos lo
zona. interpretado como dos perfiles levantados.
tipos de materiales Algo muy típico en la

76
distintos. zona.
 Figura 44. Línea 2 interpretada de radar de penetración
Figura 43. Línea 1 interpretada de Radar de penetración

Dentro de este perfil

igual podemos La tendencia de Se distingue un cambio
distinguir la presencia fracturas dentro de dentro de la transición
de tuberías a lo largo este perfil es denotada del primer metro de
de nuestra sección. mediante los círculos e nuestra sección, lo
color rojo, podemos cual es atribuido al
observar que se trata cambio de material, los
de una zona afecta por reflectores engrosan a

77
este tipo de medida que aumenta
la profundidad.
 Se puede distinguir dentro del
radargrama al igual que los anteriores
una serie de fracturas presentes dentro
del rango de los primeros 5 metros de
profundidad. Además, se identifican
tuberías presentes a los 6 y 8 metros
que cruzan de manera transversal a la
línea.

Figura 45. Línea 4 interpretada de radar de penetración.

78
 Dentro del siguiente radargrama se
puede identificar una gran cantidad de
fracturas presentes a lo largo de la
sección a una profundidad aproximada
de los 4 a 5 metros. Así como la
presencia de dos anomalías referentes
a objetos enterrados a una profundidad
de 6 metros aproximadamente.

Figura 46. Línea 5 interpretada de radar de penetración.

79
5.8 Modelo de interpretación

Con base a los resultados obtenidos de los métodos geofísicos empleados se

generaron tres modelos de la estructura somera del área estudiada, en ella se
infieren las principales estructuras geológicas que pueden presentarse (Figs. 48, 49
y 50). Los modelos geológico-geofísicos presentes se ajustaron a una profundidad
de 20 metros siendo estos de carácter somero, puesto que fue la mayor penetración
que se logró alcanzar utilizando el radar de penetración en donde se descartó la
presencia de túneles en la zona.

80
 Figura 47. Modelo geológico-geofísico I con orientación W -E
Modelo geológico 1, este perfil tiene una longitud de 75 metros de largo
orientado de Oeste a Este, se caracteriza por presentar una serie de
fracturas a lo largo de todo el perfil con longitudes registradas hasta los
5.8.1 Modelo de interpretación con base en radargramas 1

10 metros de profundidad aproximadamente. Se aprecia con cruces de

color rojo la presencia de tuberías que cortan de forma transversal en
esta zona situadas alrededor de los 3 y 4 metros de profundidad.

81
 Figura 48. Modelo geológico-geofísico II.
Modelo geológico 2, este perfil tiene una longitud de 120
metros de largo orientado de Norte a Sur, se caracteriza por
presentar una serie de fracturas a lo largo de todo el perfil con
longitudes registradas hasta los 14 metros de profundidad
aproximadamente. Se aprecia con cruces de color rojo la
presencia de tuberías que cortan de forma transversal en esta
zona situadas alrededor de los 3y 5 metros de profundidad,
5.8.2 Modelo de interpretación con base en radargramas 2

se aprecia una posible tubería u objeto enterrado a los 10

metros de profundidad.

82
 Figura 49. Modelo geológico-geofísico III
Modelo geológico 3, este perfil tiene una longitud de 195 metros de largo orientado
de Sur a Norte, se caracteriza por presentar una serie de fracturas a lo largo de todo
el perfil con longitudes registradas hasta los 8 metros de profundidad
aproximadamente. Se aprecia con cruces de color rojo la presencia de tuberías que
cortan de forma transversal en esta zona situadas alrededor de los 5 y 7 metros de
profundidad. Se presenta como el perfil con más densidad de fracturas en
comparativa con los demás.
5.8.3 Modelo de interpretación con base en radargramas 3

83
5.9 Conclusiones

En base al análisis y respuestas geológico-geofísico de la estructura somera de la

cabecera municipal de Cerro de San Pedro se descarta la presencia de túneles.
Las respuestas gravimétricas y magnéticas muestran una relación directa con el
modelo de elevación y la geología presente en la zona de estudio.
La interpretación de los Radargramas genera una imagen clara de las condiciones
del subsuelo de los primeros 20 metros del municipio de Cerro de San Pedro.
El fracturamiento presente en la zona se puede atribuir a la presencia de actividad
minera presente durante un periodo prolongado de tiempo en la que construcciones
como el templo “San Pedro” pueden llegar a ser afectadas con daños en la
infraestructura, cabe señalar que es sólo una posible respuesta a esta condición en
donde aunado a ello los daños se puede conceder a la evolución tectónica de Cerro
de San Pedro o a una combinatoria de ambos eventos.
El uso del Radar de penetración es una herramienta de suma importancia para la
caracterización del subsuelo dentro de un intervalo somero como fue empleado en
esta investigación, tuberías y objetos enterrados son fáciles de reconocer
empleando esta técnica de exploración geofísica.
Se descarta la posible localización de túneles dentro de esta zona, no se percibe
mediante la implementación de dichos métodos geofísicos la ubicación o presencia
de este tipo de estructuras por lo menos dentro de una ventana de los +/- 20 metros
de profundidad.

84
Referencias

Aranda-Gómez. (2005). El volcanismo tipo intraplaca del Cenozoico tardío en el

centro y norte de México . Boletín de la sociedad geológica mexicana, 187-
225.
Bello-Barradas, D. C. (1963). contrasting styles of laramide folding across thewest-
central margin of the cretaceous valles OF SAN LUIS POTOSI. Revista
Mexicana de Ciencias Geológicas,, 97-111.
Blakely. (1996). Potential Theory in Gravity & Magnetic. Melburne: Cambridge
University Press.
Burger, H. (1996). Exploration Geophysics of the Shallow Subsurfac. United States:
N.P.
Burger, H. (2006). Introduction to Applied Geophysics. [Link]: N.A.
Busquet. (1996). Geofísica ambiental ténicas no destructivas para el reconocimiento
de zonas contaminadas por vertidos. Acta Geologica Hispanica, 73-82.
Butler, R. (1992). Paleomagnetism: Magnetic Domains to Geologic Terranes.
Portland, Oregon : Blackwell Scientific Publications.
Caddey. (1993). Unfolding the Geology of the West. USA: The geological society of
america .
Carrillo. (1971). La plataforma Valles- San Luis Potosí. Boletin de la asosciación
mexicana de geologos petroleros, 1-15.
Coremi. (1996). Monografía geológicominera del estado de San Luis Potosí.
Secretaria de Energía, Minas e Industria Paraestatal, Consejo de Recursos
Minerales, 217.
Cukavac. (2008). Reconstruction of Buried Objects by Implementation of Ground
Penentrating Radar Technique. Proceeding of the International Conference,
333-338.
Grauch. (1966). High speed gravity and magnetic calculations of uniform cylindrical
bodies of arbitrary cross-section and finite length. Pure and Applied
Geophysics PAGEOPH, 48-53.

85
Hérnandez- Sánchez Mejorada, L. R. (1968). Carta Geológica de la República
Mexicana escala 1:2 000 000. México, Instituto de geología: INEGI.
Identification of cavities by extraction of characteristic parameter from ground
probing radar reflection data. (1997). Tesis Doctoral, Christian-Albrechts
Universität, 142-144.
INEGI. (1998). Carta geologica minera F14-4. Cd. de México : Servicio geológico
mexicano .
J.J, A. G. (1978). Metamorphism, mineral zoning, and paragenesis in the San Martín
mine, Zacatecas, Mexico. Boletin de la sociedad geologica mexicana , 223-
229.
Kearey, P. (2002). An introduction to geophysical exploration. [Link]: Blackwell
Sicience.
Keith. (2004). Geology of the Cerro San Pedro, Precious and Base Metals Mining
District San Luis Potosí, Mexico. Society of Economic Geologists, 220.
Labarthe. (1982). ronología y distribución espacial del vulcanismo en el Campo
Volcánico de San Luis Potosí. BOLETÍN DE LA SOCIEDAD GEOLÓGICA
MEXICANA, 287-303.
Milsom. (2003). Field Geophysics. [Link]: Wiley.
Mironov. (1977). Curso de prospección gravimpetrica . REVERTE.
Nabighian. (1988). Electromagnetic theory for geophysical applications. Sydeny:
Cambridge university press.
Nieto Samaniego, A., & Alanís-Álvarez. (2005). La mesa central de México:
estratigrafía, estructura,y evolución tecotonica cenozoica. Boletin de la
sociedad geologica mexicana, 285-318.
Olohoeft. (2006). Ground-Penetrating Radar (GPR) Mapping as a Method for
Planning Excavation Strategies. Society for American Archaeology, 78-84.
Petersen, M. (2001). Geology of the Cerro San Pedro porphyry-related gold-silver
deposit, San Luis Potosí, Mexico. THE SOCIETY OF ECONOMIC
GEOLOGISTS, 157-162.
Pueyo. (2016). Ground penetrating radar evaluation of the internal structure of fluvial
tufa deposits (Devanos-A ´ navieja system, NE Spain). Geophysical Journal
International, 557-573.
Rubin, J. K. (1988). Mineralogy and geochemistry of the San Martín skarn deposits.
Economic Geology, 1782-1792.

86
Tavera. (2006). Métodos no-destructivos para la mejora del rendimiento de una red
hidráulica: casos de estudio con el georadar. Congreso Latinoamericano de
Hidráulica.
Telford. (1990). Applied geophysics. Melbourne: Press Syndicate of the University
of Cambridge.
Vega. (2002). “Radar de Subsuelo. Evaluación para Aplicaciones en Arqueología y
en Patrimonio Histórico-Artistico. Tesis Docotroal,Universidad Politécnica de
Cataluña,España., 101-105.
Winterbourner. (1999). Alteration and mineralization of the San Pedro porphyry,
Cerro San Pedro, San Luis Potosí, Mexico. The geological society of america
, 208.
Zamora. (2005). Aplicación de la Prospección Geofísica por Georradar (GPR) a la
Delimitación y Caracterización de Estructuras Arqueológicas. Geogaceta,
131-134.

87
ANEXO A
Línea Fecha Hora medición Campo Estación
Magnético
medido

1 17/11/2016 [Link] 41978.61 1

1 17/11/2016 [Link] 41973.34 2
1 17/11/2016 [Link] 41958.01 3
1 17/11/2016 [Link] 41927.1 4
1 17/11/2016 [Link] 41943.03 5
1 17/11/2016 [Link] 41966.61 6
1 17/11/2016 [Link] 42006.95 7
1 17/11/2016 [Link] 42100.97 8
1 17/11/2016 [Link] 42012.07 9
1 17/11/2016 [Link] 41984.63 10
1 17/11/2016 [Link] 41985.4 11
1 17/11/2016 [Link] 41979.56 12
1 17/11/2016 [Link] 42003.12 13
1 17/11/2016 [Link] 42022.97 14
1 17/11/2016 [Link] 42011.26 15
1 17/11/2016 [Link] 42004.25 16
1 17/11/2016 [Link] 41973.07 17
1 17/11/2016 [Link] 41915.89 18
1 17/11/2016 [Link] 41845.8 19
1 17/11/2016 [Link] 41930.33 20

88
Anexo B

Línea Estación Elevación Gravedad Desviación

medida estándar

0 0 26.9041 3180.516 0.017

1 0 27.6365 3184.965 0.023
1 1 26.9041 3184.989 0.024
1 2 26.9041 3184.975 0.017
1 3 26.9041 3184.891 0.029
1 4 27.1482 3184.962 0.023
1 5 27.1482 3184.955 0.015
1 6 27.3924 3184.922 0.016
1 7 27.6365 3184.934 0.013
1 8 27.6365 3184.919 0.016
1 9 27.3924 3184.961 0.013
1 10 26.9041 3184.98 0.018
1 11 26.66 3184.943 0.017
1 12 26.66 3184.951 0.023
1 13 26.9041 3185.015 0.014
1 14 26.4158 3185.034 0.014
1 15 26.1717 3185.038 0.014
1 16 26.1717 3185.132 0.016
1 17 26.1717 3185.129 0.018
1 18 26.66 3185.099 0.019
1 19 27.3924 3185.162 0.012
1 20 28.3689 3185.223 0.019

89
Anexo C

Punto X Y Z
1 314456.837 2457960.9 2035.751
2 314457.54 2457965.86 2035.905
3 314458.158 2457970.82 2035.905
4 314458.74 2457975.78 2035.976
5 314461.059 2457980.56 2037.896
6 314460.167 2457985.68 2036.362
7 314460.609 2457990.62 2036.524
8 314461.252 2457995.55 2037.396
9 314459.914 2458000.35 2037.855
10 314458.596 2458005.14 2038.103
11 314457.171 2458009.92 2038.398
12 314456.113 2458014.75 2038.994
13 314455.231 2458019.61 2039.788
14 314454.855 2458024.55 2040.358
15 314454.08 2458029.51 2040.438
16 314453.48 2458034.4 2040.991
17 314452.834 2458039.24 2041.72
18 314452.041 2458044.11 2042.354
19 314452.017 2458047.65 2042.986
20 314448.986 2458049.4 2040.424

90