Jjjjjjonf - INFORME SISMICOjjjjjjjjjjjj
Jjjjjjonf - INFORME SISMICOjjjjjjjjjjjj
Jjjjjjonf - INFORME SISMICOjjjjjjjjjjjj
INGENIERO GEOFÍSICO
AREQUIPA – PERU
2018
UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA
FACULTAD DE GEOLOGÍA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERÍA GEOFÍSICA
INGENIERO GEOFÍSICO
AREQUIPA – PERU
2018
DEDICATORIA
i
AGRADECIMIENTO
Quiero agradecer a la empresa Sotelo & Asociados S.A.C. por otorgarme el permiso para utilizar los
datos de campo adquiridos para el “Proyecto Centro de Convenciones Cusco Word, Urubamba,
Cusco” y la realización de este estudio.
A mis compañeros de trabajo quienes conforman la empresa Sotelo & Asociados S.A.C. los cuales
son como mi segunda familia, en especial a la Sra. Rocio Giraldo, Mariano Ignacio, Cristhian Alhuay,
Martin Ruiz, Brayam Velázque, Julio Panca, Waldir Mamani, Marco Pérez, Arnold Quispe, por
ofrecerme su amistad, por compartir experiencias y conocimientos en el trabajo.
A la Dra. Yenny Giraldo por su hospitalidad brindada en el inicio de mi carrera, por sus consejos, por
su paciencia, por su confianza y amistad.
Al Ing. Hebert Sotelo por su profesionalidad, sus sabios consejos, por compartir sus experiencias y
enseñanzas en el campo de la geotecnia, por su amistad incondicional y calidad de persona.
A mi alma mater, Universidad Nacional de San Agustín y mis maestros de la escuela de Ingeniería
Geofísica, a quienes recuerdo con mucho cariño y respeto, Dr. Armando Minaya, Dr. Orlando
Macedo, Mg. Ing. Edgar Gonzales, Mg., e Mg. Ing. Sebastián Zúñiga y al Mg. Ing. Javier Ticona.
A mi asesor Dr. Jorge Soto Vásquez Por sus observaciones y el tiempo dedicado para la revisión de
esta tesis, con el fin de enriquecerla con sus conocimientos.
A mis compañeros de la Promoción XXX “EMIL WIECHERT”, por todos los momentos de alegría que
compartimos en clases durante nuestra formación profesional.
A mis padres, por sus consejos, enseñanzas, sobre todo por su eterno apoyo en todo lo que he
decidido, siempre han luchado para darme lo mejor y por eso los admiro, gracias por toda la felicidad
que me han brindado.
A mi gran amigo Julio Chozo Amaya, por la calidad de persona, por las aventuras y desafíos que
supimos vencer durante los proyectos realizados.
Por último, quiero agradecer a las personas quienes fueron parte de mi formación profesional, gracias
por el aprendizaje que cada uno de ustedes aporto a mi persona y que de una u otra manera llegaron
a conocerme.
ii
RESUMEN
El objetivo general de este estudio fue evaluar las condiciones y características geofísicas y
geotécnicas del suelo que cubre el área de estudio donde se proyecta la construcción de estructuras
importantes.
El primer horizonte podría asociarse a un suelo poco rígido, para el segundo horizonte se puede
deducir que se trata de un material rígido para el tercer horizonte se podría deducir que se trata de un
material medianamente rígido asociado posiblemente gravas con presencia de arcillas, arenas y
limos. Por otro lado, según los estudios de Resistividad de Tomografía Eléctrica. tenemos tres zonas
clasificadas según su resistividad que varían desde 85 a 340 ohm*m. las cuales nos hacen inferir que
posiblemente haya presencia de humedad a profundidad y hacen denotar los horizontes que resaltan
con los estudios de MASW. Adicionalmente se pudo obtener parámetros seudo dinámicos del
subsuelo con la ayuda de la refracción sísmica y MASW.
The study area is located in the area called Patashuaylla - 5th stop in the district and province of
Urubamba department of Cusco where the establishment of the "CENTRO DE CONVENCIÓN
CUSCO WORD" was planned.
The general objective of this study was to evaluate the geophysical and geotechnical conditions and
characteristics of the soil that covers the study area where the construction of important structures is
planned. The field work was carried out in two stages. The first stage consisted in the acquisition of 8
seismic refraction lines, each of which had a length of 46m., two lines of 46m joined, lines overlapped
on three occasions to obtain only 3 lines of 6 for processing. There were also 8 MASW, located in the
same lines of seismic refraction. On the other hand, to complement and verify the water infiltration or
saturation of the subsoil, two geoelectrical resistivity tomography stretches of lengths of 90m. were
performed. In the second stage the topographic survey with drone was done to have a broader view of
the work area serving as a great help for interpretation in the results of the studies.
In the present work of geotechnical investigation the methodology of work for the different applied
methods is described as well as the interpretation of the same giving interesting results, it was
possible to identify 3 main horizons where it can be associated to quaternary materials with
differences in their composition. The first horizon could be associated with a slightly rigid ground, for
the second horizon it can be deduced that it is a rigid material for the third horizon, it could be deduced
that it is a moderately rigid material associated with gravels, sands and slime On the other hand,
according to the studies of Resistivity of Electrical Tomography. we have three zones classified
according to their resistivity that vary from 85 to 340 ohm * m. which make us infer that there may be
moisture present at depth and denote the horizons that stand out with the studies of MASW.
Additionally, pseudo-dynamic parameters of the subsoil could be obtained with the help of seismic
refraction and MASW.
iv
ÍNDICE
DEDICATORIA .......................................................................................................................................... I
AGRADECIMIENTO ................................................................................................................................. II
RESUMEN ............................................................................................................................................... III
ABSTRACT ............................................................................................................................................ IV
LISTA DE CUADROS ............................................................................................................................ XI
LISTA DE FIGURAS ............................................................................................................................ XIII
LISTA DE FOTOS ............................................................................................................................... XVII
CAPITULO I ............................................................................................................................................. 1
INTRODUCCIÓN ...................................................................................................................................... 1
1.1 Ubicación .................................................................................................................................... 2
1.2 Accesibilidad .............................................................................................................................. 3
1.3 Clima y vegetación ..................................................................................................................... 3
1.4 Justificación del problema ........................................................................................................ 4
1.5 Formulación del problema ......................................................................................................... 4
1.5.1 Definición del problema................................................................................................ 4
1.5.2 Formulación del problema ........................................................................................... 4
1.6 Alcances y Limitaciones ............................................................................................................ 4
1.6.1 Alcances ........................................................................................................................ 4
1.6.2 Limitaciones .................................................................................................................. 4
1.7 Variables e Indicadores.............................................................................................................. 5
1.7.1 Independientes .............................................................................................................. 5
1.7.2 Dependientes................................................................................................................. 5
1.7.3 Indicadores .................................................................................................................... 5
1.8 Objetivos ..................................................................................................................................... 5
1.8.1 Objetivos Generales...................................................................................................... 5
1.8.2 Objetivos Específicos ................................................................................................... 5
1.9 Hipótesis ..................................................................................................................................... 6
1.10 Antecedentes de la investigación ............................................................................................. 6
CAPITULO II ............................................................................................................................................ 7
GEOLOGÍA DEL ÁREA DE ESTUDIO .................................................................................................... 7
2.1 Geomorfología ............................................................................................................................ 8
v
2.1.1 Geomorfología regional ................................................................................................ 8
2.1.1.1 Cordillera Oriental ............................................................................................ 8
2.1.1.2 Valle del Vilcanota Urubamba.......................................................................... 8
2.1.2 Unidades Geomorfológicas principales ...................................................................... 9
2.1.2.1 Cerros escarpados .......................................................................................... 9
2.1.2.2 Cárcavas.......................................................................................................... 9
2.1.2.3 Conoide de depósito de flujos de barro I ......................................................... 9
2.1.2.4 Conoide de depósitos de flujo de barro II ........................................................ 9
2.1.2.5 Terrazas fluviales. ............................................................................................ 9
2.1.2.6 Grietas y cavernas de disolución. .................................................................. 10
2.2 Geología local ........................................................................................................................... 10
2.2.1 Estratigrafía ................................................................................................................. 10
2.2.1.1 Formación maras (Ki-ma) .............................................................................. 10
2.2.1.2 Depósitos Fluviales (Q-f) ............................................................................... 11
2.2.1.3 Depósitos coluviales (Q-co) ........................................................................... 11
2.2.1.4 Depósitos aluviales (Q-al) .............................................................................. 12
2.3 Geomorfología local ................................................................................................................. 13
2.4 Drenaje ...................................................................................................................................... 13
2.4.1 Rio Vilcanota ............................................................................................................... 13
CAPITULO III.......................................................................................................................................... 15
FUNDAMENTOS TEÓRICOS ................................................................................................................ 15
3.1 Ondas sísmicas ........................................................................................................................ 16
3.1.1 Parámetros para la descripción de ondas ................................................................ 16
3.1.2 Ondas de cuerpo ......................................................................................................... 17
3.1.2.1 Ondas primarias o de compresión ................................................................. 17
3.1.2.2 Ondas secundarías o de corte ....................................................................... 18
3.1.3 Ondas superficiales .................................................................................................... 19
3.1.3.1 Ondas Love ................................................................................................... 19
3.1.3.2 Ondas Rayleigh ............................................................................................. 20
3.2 Refracción sísmica ................................................................................................................... 20
3.2.1 Ley de Snell ................................................................................................................. 21
3.2.2 Ley de Huygens ........................................................................................................... 22
3.2.3 Ley de Fermat .............................................................................................................. 22
vi
3.2.4 Conceptos básicos de metodología de cálculo en refracción sísmica .................. 23
3.2.4.1 Curvas tiempo-distancia(dromocronas) ......................................................... 23
3.2.4.2 Método de tiempos de intercepto ................................................................... 23
3.2.4.3 Método de distancia crítica ............................................................................ 28
3.2.4.4 Capas buzantes ............................................................................................. 29
3.2.5 Fuente de energía ....................................................................................................... 31
3.2.6 Sismograma................................................................................................................. 33
3.2.7 Ventajas y limitaciones en la aplicación del método ............................................... 34
3.3 Análisis Multicanal de Ondas Superficiales (MASW) ............................................................ 36
3.3.1 Ondas Rayleigh ........................................................................................................... 36
3.3.2 Dispersión de ondas superficiales ............................................................................ 38
3.3.3 Velocidad de fase ........................................................................................................ 39
3.3.3.1 Calculo de velocidad de fase ......................................................................... 39
3.3.4 Curva de dispersión .................................................................................................... 41
3.3.4.1 Construcción de la curva de dispersión ......................................................... 43
3.3.5 Modelo inicial .............................................................................................................. 44
3.3.6 Inversión de datos ...................................................................................................... 46
3.3.7 Profundidad de investigación .................................................................................... 47
3.3.8 Ventajas y limitaciones del método ........................................................................... 48
3.4 Clasificación de suelos en base a las velocidades de corte (Vs) ......................................... 49
3.5 Constantes elásticos dinámicos del suelo............................................................................. 51
3.5.1 Constantes elásticas dinámicas con velocidades de propagación. ....................... 51
3.5.2 Relación de Poisson ................................................................................................... 52
3.5.3 Módulo de corte o Rigidez.......................................................................................... 52
3.5.4 Módulo de elasticidad o Young ................................................................................. 52
3.5.5 Modulo volumétrico o Bulk ........................................................................................ 53
3.6 Tomografía geoeléctrica .......................................................................................................... 53
3.6.1 Ley de ohm y resistividad .......................................................................................... 54
3.6.2 Conductividad eléctrica .............................................................................................. 54
3.6.3 Resistividad aparente ................................................................................................. 55
3.6.4 Propiedades eléctricas de los materiales terrestres ................................................ 56
3.6.5 Dispositivos Electródicos .......................................................................................... 58
3.6.5.1 Dispositivo Wenner – Schlumberger .............................................................. 59
vii
3.6.5.2 Dispositivo Wenner ........................................................................................ 60
3.6.5.3 Dispositivo Polo- Dipolo ................................................................................. 60
3.6.5.4 Dispositivo Dipolo – Dipolo ............................................................................ 60
3.6.5.5 Dispositivo Polo – Polo .................................................................................. 61
3.6.6 Características de los dispositivos Electródico ....................................................... 61
3.6.7 Pseudosección ............................................................................................................ 62
3.6.8 Profundidad de investigación .................................................................................... 63
3.6.9 Inversión de datos ...................................................................................................... 66
3.6.10 Ventajas y limitación del método ............................................................................... 68
CAPITULO IV ......................................................................................................................................... 69
INSTRUMENTACIÓN GEOFÍSICA Y ADQUISICIÓN DE DATOS ......................................................... 69
4.1 Refracción Sísmica – Análisis multicanal de ondas superficiales (MASW) ........................ 70
4.1.1 Equipo de adquisición de datos geosísmicos. ......................................................... 70
4.1.1.1 Sismógrafo geode .......................................................................................... 70
4.1.1.2 Accesorios de sismógrafo geode ................................................................... 71
4.1.1.3 Geófono ......................................................................................................... 72
4.1.1.4 Cable sísmico o cable conector de geófono .................................................. 73
4.1.1.5 Cable Trigger ................................................................................................. 73
4.1.1.6 Fuente de energía ......................................................................................... 74
4.1.2 Parámetros de adquisición ........................................................................................ 74
4.1.3 Descripción de trabajos de campo – método Refracción sísmica.......................... 75
4.1.3.1 Reconocimiento de campo ............................................................................ 76
4.1.3.2 Tendido de línea de refracción sísmica ......................................................... 77
4.1.3.3 Adquisición de datos de refracción sísmica ................................................... 78
4.1.3.4 procesamiento de datos de refracción sísmica .............................................. 79
4.1.4 Descripción de trabajos de campo – método MASW ............................................... 82
4.1.4.1 Tendido de línea sísmica para el método MASW .......................................... 83
4.1.4.2 Adquisición de datos del método MASW ....................................................... 83
4.1.4.3 Procesamiento de datos del método MASW ................................................. 84
4.2 Tomografía Geoeléctrica.......................................................................................................... 87
4.2.1 Equipo de adquisición de datos geoeléctricos ........................................................ 87
4.2.1.1 Resistivímetro ................................................................................................ 87
4.2.1.2 Electrodos de corriente .................................................................................. 88
viii
4.2.1.3 Cables............................................................................................................ 88
4.2.1.4 Otros accesorios ............................................................................................ 89
4.2.2 Parámetros de adquisición ........................................................................................ 89
4.2.3 Descripción de trabajos de campo ............................................................................ 90
4.2.3.1 Tendido de línea ............................................................................................ 90
4.2.3.2 Adquisición de datos ...................................................................................... 92
4.2.3.3 procesamiento de datos geoeléctricos........................................................... 93
4.3 Levantamiento fotométrico con Drone ................................................................................... 96
4.3.1 Equipo de levantamiento fotográfico ........................................................................ 96
4.3.1.1 Drone - Phantom 4 pro .................................................................................. 96
4.3.1.2 Características principales ............................................................................. 98
4.3.2 Descripción de trabajos de campo ............................................................................ 98
4.3.2.1 Plan de vuelo ................................................................................................. 98
4.3.2.2 Apoyo topográfico .......................................................................................... 99
4.3.2.3 Toma fotográfica ............................................................................................ 99
4.3.2.4 Procesamiento de imágenes fotográficas .................................................... 100
CAPITULO V ........................................................................................................................................ 102
ANÁLISIS E INTERPRETACIÓN DE RESULTADOS.......................................................................... 102
5.1 Levantamiento fotométrico de drone.................................................................................... 103
5.2 Perfiles geosísmico de refracción sísmica .......................................................................... 104
5.2.1 Perfil geosísmico RS – 01......................................................................................... 104
5.2.2 Perfil geosísmico RS – 02......................................................................................... 104
5.2.3 Perfil geosísmico RS – 03......................................................................................... 105
5.2.4 Perfil geosísmico RS – 04......................................................................................... 105
5.2.5 Perfil geosísmico RS – 05......................................................................................... 106
5.3 Perfiles unidimensionales de MASW .................................................................................... 107
5.3.1 Perfil unidimensional MASW – 01 ............................................................................ 107
5.3.2 Perfil unidimensional MASW – 02 ............................................................................ 108
5.3.3 Perfil unidimensional MASW – 03 ............................................................................ 108
5.3.4 Perfil unidimensional MASW – 04 ............................................................................ 109
5.3.5 Perfil unidimensional MASW – 05 ............................................................................ 109
5.3.6 Perfil unidimensional MASW – 06 ............................................................................ 110
5.3.7 Perfil unidimensional MASW – 07 ............................................................................ 110
ix
5.3.8 Perfil unidimensional MASW – 08 ............................................................................ 110
5.4 Parámetros pseudo – dinámicos .......................................................................................... 111
5.5 Perfiles geoeléctricos de Tomografía Geoeléctrica - TG .................................................... 115
5.5.1 Perfil Geoeléctrico TG – 01....................................................................................... 116
5.5.2 Perfil Geoeléctrico TG – 02....................................................................................... 116
5.6 Correlación de datos geofísicos ........................................................................................... 117
CONCLUSIONES ................................................................................................................................. 120
RECOMENDACIONES ......................................................................................................................... 121
BIBLIOGRAFÍA .................................................................................................................................... 122
ANEXOS
ANEXO A: EXPLORACIONES DE CAMPO
ANEXO A1: REFRACCIÓN SÍSMICA
ANEXO A2: MASW - ANÁLISIS MULTICANAL DE ONDAS SUPERFICIALES
ANEXO A3: TG - TOMOGRAFÍA GEOELÉCTRICA
ANEXO A4: PARÁMETROS PSEUDO-DINÁMICOS
x
LISTA DE CUADROS
xi
Cuadro 22. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m,
realizados en las líneas LS 02 – MASW 07................................................................... 114
Cuadro 23. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 24 m,
realizados en las líneas LS 05 – MASW 08................................................................... 115
xii
LISTA DE FIGURAS
Figura 1. Sistema hidrológico del río Urubamba (Hydrology of the Urubamba River, by Francisco
Olivera, Ph.D.; Center for Research in Water Resources; University of Texas at Austin;
January 1996). ..................................................................................................................... 14
Figura 2. Trayecto río Urubamba (http://riosdelplaneta.com) .............................................................. 14
Figura 3. Desplazamiento de ondas P a través de un medio solido. ................................................... 16
Figura 4. Desplazamiento de ondas P a través de un medio solido (Keary et al.2002). ..................... 17
Figura 5. Propagación de ondas P (Alsadi, H, N,2017). ...................................................................... 17
Figura 6. Desplazamiento de ondas S a través de un medio solido (Keary et al.2002). ..................... 18
Figura 7. Propagación de ondas SV (Alsadi, H, N,2017). ................................................................... 19
Figura 8. Propagación de ondas SH (Alsadi, H, N,2017). ................................................................... 19
Figura 9. Propagación de ondas Love (Alsadi, H, N,2017). ................................................................ 20
Figura 10. Propagación de ondas Rayleigh (Alsadi, H, N,2017). ........................................................ 20
Figura 11. Camino de rayo sísmico en un medio (Alsadi, H, N,2017). ................................................ 21
Figura 12. Un rayo sísmico impacta el límite de dos medios con diferentes velocidades (v2 > v1) y
se divide en rayos sísmicos refractados y reflectados (Stefan J.,2010) ............................. 22
Figura 13. Propagación de onda plana de acuerdo al Principio de Huygens. a) través de un medio
homogéneo donde la velocidad es constante. b) medio no homogéneo, de velocidad
que está cambiando uniformemente a través de la dirección de propagación (Alsadi, H,
N,2017). ............................................................................................................................. 22
Figura 14. Descripción simple de dos capas, límites paralelos y la correspondiente curva de
tiempo-distancia (Bruce B. Redpath,1973). ....................................................................... 24
Figura 15. Esquema de caso de múltiples capas y Curva de tiempo-distancia. (Bruce B.
Redpath,1973). .................................................................................................................. 26
Figura 16. Ejemplo de interfaz de inmersión y conceptos de disparo inverso y velocidad aparente
(Bruce B. Redpath,1973). .................................................................................................. 30
Figura 17. Se muestra a la izquierda, fuente de energía - caída de peso y a la derecha registro
sísmico (http://www.geofisicaexploraciones.com). ............................................................. 32
Figura 18. Se muestra a la izquierda, fuente de energía – cartucho de 600 gramos de Pyroblast y
a la derecha registro sísmico (http://www.geofisicaexploraciones.com). ........................... 32
Figura 19. Esquema de fuente de energía para refracción sísmica,
(https://es.slideshare.net/JDaviSP/9-ensayo-de-sismica-de-refracción). ........................... 33
Figura 20. Sismograma de registro de ondas P – refracción sísmica (fuente propia) ......................... 34
xiii
Figura 21. Se muestra esquema de configuración de campo del método MASW (Park et al. 2004). . 36
Figura 22. Se muestra en la parte superior la propagación de ondas Rayleigh (retrógrado) y en la
parte inferior, movimiento de partículas elípticas de la propagación de ondas Rayleigh
(Van Der Hilst 2004). ......................................................................................................... 37
Figura 23.Dispersión de ondas de superficie en un medio heterogéneo (SAGEEP, 2003 & Lopes,
I. et al., 2008). .................................................................................................................... 38
Figura 24.Registro sísmico con su respectivo modelo inicial de velocidad de onda de acuerdo a
las ondas superficiales según las características del subsuelo (SAGEEP, 2003). ............ 39
Figura 25.Cálculo de la velocidad de la fase, esquema de métodos de transformación (fuente:
Koichi Hayashi,2003). ........................................................................................................ 41
Figura 26. Trasformadas en dominio de frecuencia – MASW (Koichi Hayashi,2003). ........................ 41
Figura 27. A la izquierda registro de forma de onda, incluyendo modo alto (Koichi Hayashi,2003),
a la derecha una ilustración de cómo se observaría en un modelo 3d (Ásta
Ólafsdóttir,2014) ................................................................................................................ 42
Figura 28. Comparación de curvas de dispersión y su respectivo modelo unidimensional para
materiales con distintas características (Koichi Hayashi,2003).......................................... 43
Figura 29. Modelo inicial basado en 1/3 transformación de longitud de onda (Koichi ayashi,2003). .. 46
Figura 30. Modelo final o perfil unidimensional de ondas de corte – MASW (Koichi Hayashi,2003). .. 47
Figura 31. Máxima profundidad de investigación - MASW. ................................................................. 48
Figura 32. Efecto de onda P en el medio de propagación................................................................... 51
Figura 33. Efecto de onda S en el medio de propagación................................................................... 52
Figura 34. Recopilación de datos en 2D mediante el sistema de resistividad de múltiples
electrodos (Loke et al., 2013)............................................................................................. 53
Figura 35. Diagrama de resistividad (Loke et al., 2013). ..................................................................... 54
Figura 36. Valores de resistividad representativos para rocas y minerales (Palacky, 1987). .............. 57
Figura 37. Dispositivo Tetraelectródico para la medida de la resistividad del suelo. (Moller et al.,
1995).................................................................................................................................. 58
Figura 38. Esquema Electródico del dispositivo Wenner – Schlumberger .......................................... 59
Figura 39. Esquema Electródico del dispositivo Wenner. ................................................................... 60
Figura 40. Esquema Electródico del dispositivo Polo - dipolo. ............................................................ 60
Figura 41. Esquema Electródico del dispositivo Dipolo - dipolo. ......................................................... 61
Figura 42. Esquema Electródico del dispositivo Polo – polo ............................................................... 61
Figura 43. Pseudosección de un dispositivo Electródico Dipolo – dipolo (Fuente propia)................... 63
xiv
Figura 44.Profundidades medias de investigación sobre un medio homogéneo de resistividad
para los dispositivos obtenidas para una distancia interelectródica de 5 metros y
distintos niveles de medición (Edwards 1977). Valores idénticos para los dispositivos
Polo-Dipolo y Wenner–Schlumberger. ............................................................................... 65
Figura 45. Proceso de inversión de datos (Martínez Pagan, 2006). ................................................... 66
Figura 46. Corrección de topografía para inversión de datos-TGE (Phil Sirles, 2016). ....................... 67
Figura 47. Sismógrafo geode (Geometrics Inc., 2009). ....................................................................... 70
Figura 48. Caja NIB (Geometrics Inc., 2009). ..................................................................................... 71
Figura 49. Cable conector de interfaz (Geometrics Inc., 2009). .......................................................... 71
Figura 50. Cable sísmico (Geometrics Inc., 2009). ............................................................................. 73
Figura 51. Armado de línea sísmica (Fuente propia) .......................................................................... 78
Figura 52. Se muestra la distribución de puntos de disparo, el primero se ejecuta antes del 1er
geófono y a una distancia igual a la separación entre geófonos, y el segundo disparo
entre el 6to y 7mo geófono. ............................................................................................... 78
Figura 53. El tercer punto de disparo se realizó entre el 12vo y 13vo geófono. .................................. 78
Figura 54. El 4to punto de disparo se da entre los geófonos 18 y 19, y El ultimo de, se ejecuta
después del 24vo geófono a una distancia igual a la separación entre geófonos. ............ 79
Figura 55.Se muestra el ultimo de punto de disparo, será el primer punto de disparo de la
segunda línea sísmica. ...................................................................................................... 79
Figura 56. Registros de ondas compresionales - LS01. ...................................................................... 80
Figura 57. La figura de la izquierda, muestra el picado de la llegada de la onda de compresión, a
la derecha los respectivos puntos que forma la dromocrónica, la cual muestra los
cambios de velocidades en el terreno de la zona de estudio. ............................................ 80
Figura 58. Perfil de refracción sísmica (Fuente propia). ...................................................................... 81
Figura 59. Perfil de refracción sísmica con proceso de inversión (tomografía sísmica), (Fuente
propia)................................................................................................................................ 81
Figura 60. Distribución de puntos de disparo del ensayo MASW........................................................ 83
Figura 61. a) Tiempo de llegada de la onda de corte para diferentes distancias. b) Curva de
dispersión: velocidad de fase vs modo fundamental (Fuente propia). ............................... 85
Figura 62. Modulo fundamental de la curva de dispersión- MASW (Fuente propia) ........................... 86
Figura 63. Perfil unidimensional de ondas de corte (Vs)- MASW (Fuente propia) .............................. 86
Figura 64. Syscal Pro Swicht 48(http://www.iris-instruments.com/syscal-prosw.html). ....................... 88
Figura 65. Cable INDECO #16 ............................................................................................................ 89
xv
Figura 66. Esquema de tendido de línea de Tomografía geoeléctrica. ............................................... 91
Figura 67. Esquema de adquisición mediante tomografía a eléctrica. ................................................ 93
Figura 68. Esquema de adquisición mediante Tomografía Geoeléctrica (Modificada de Chávez et
al., 2011). ........................................................................................................................... 93
Figura 69. Evaluación de data cruda, Pseudosección - TG 01 (Fuente propia). ................................. 94
Figura 70. Esquema del proceso de inversión - TG 01 (Fuente propia). ............................................. 94
Figura 71. Características generales de inversión, programa Res2DInv - TG 01 (Fuente propia). ..... 95
Figura 72. Perfil geoeléctrico, Res2DInv - TG 01 (Fuente propia). ..................................................... 95
Figura 73. Phantom 4 pro (https://www.dji.com/phantom-4). .............................................................. 96
Figura 74. Diagrama de drone Phantom 4 pro (Manual Phantom 4 pro, 2017)................................... 97
Figura 75. Diagrama de control remoto Phantom 4 pro (Manual Phantom 4 pro, 2017). .................... 97
Figura 76. Diagrama de control remoto Phantom 4 pro (Manual Phantom 4 pro, 2017). .................. 100
Figura 77.Modelo 3d de nube densa de la zona de estudio a partir de fotometría (Agisoft
Photoscan). ...................................................................................................................... 101
Figura 78. A) reconstrucción de imágenes utilizadas para modelo inicial, B) modelo de elevación
digital (DEM) (color azul: cota menor, color verde a amarillo: cota intermedia, color
rojo: cota mayor), C) modelo orto mosaico (Agisoft Photoscan). ..................................... 103
Figura 79. Velocidades VP para líneas de refracción sísmica. ......................................................... 106
Figura 80. Sección de Tomografía Geoeléctrica – TG01 .................................................................. 116
Figura 81. Sección de Tomografía Geoeléctrica – TG02 .................................................................. 117
Figura 82. Sección de Refracción sísmica 03 y su correlación con MASW03 y MASW06................ 118
Figura 83. Sección de Tomografía geoeléctrica 01 y su correlación con MASW03 .......................... 119
xvi
LISTA DE FOTOS
Foto 1. Se observan las dos cárcavas principales que han erosionado los depósitos de flujo de
barro (Guadalupe, et al,2002). ................................................................................................ 12
Foto 2. Se observa las siguientes geoformas: cerros escarpados, cárcava principal, conoide de
depósito de flujos de barro, terraza fluvial (Guadalupe, et al,2002). ....................................... 12
Foto 3. Observamos cárcavas y asentamiento diferencial en la Formación Geológica Maras,
conformado por yeso y lutitas violáceas (Guadalupe, et al,2002). .......................................... 12
Foto 4. Se observa el estado de la explotación de Cantera (Yanahuara) y los materiales de
gravas y arenas (Guadalupe, et al,2002). ............................................................................... 12
Foto 5. Accesorios de sismógrafo geode (Fuente propia). .................................................................. 72
Foto 6. Panorama fotográfico del área de estudio (Fuente propia) ..................................................... 76
Foto 7. Tendido de línea sísmica (Fuente propia) ............................................................................... 77
Foto 8. Proceso de cambio de shot en el área de trabajo- MASW01 (Fuente propia)......................... 83
Foto 9. tendido de cables de línea de Tomografía geoeléctrica (Fuente propia). ................................ 91
Foto 10. Instalación de cable en el Resistivímetro (Fuente propia). .................................................... 92
Foto 11. Adquisición de datos – TG01 (Fuente propia). ...................................................................... 92
Foto 12. Toma de fotos desde la ubicación de manejo de drone (Fuente propia). .............................. 99
xvii
CAPITULO I
INTRODUCCIÓN
La exploración geotécnica del subsuelo se ve limitada en el medio profesional a los costos de esta,
por lo cual y dado su gran peso dentro de los estudios geotécnicos en muchas ocasiones se deja muy
incompleta.
Estos métodos geofísicos, integrados a las técnicas tradicionales, pueden contribuir significativamente
en la optimización de los análisis ingenieriles de un proyecto por la gran cantidad de informaciones
que pueden proporcionar en forma no invasiva, la rápida adquisición de los datos, posibilidad de
trabajar en diferentes tipos de ambientes sin generar ruido, destrucción o impactos al medio urbano,
industrial o ambiental.
En pocos proyectos se utilizan las técnicas de prospección geofísica, debido, en gran medida, al
desconocimiento de esta valiosa herramienta.
El trabajo está conformado por cinco capítulos, en el primer capítulo se describe la ubicación,
accesibilidad, clima y vegetación, en el segundo capítulo se describe la geología que conforma la
zona de estudio, en el tercer capítulo se describe los fundamentos teóricos de cada método realizado
en este trabajo, en el cuarto capítulo se describe la instrumentación geofísica y procesado de datos,
en el quinto capítulo se describe el análisis e interpretación de datos. y por último daremos a conocer
las conclusiones y recomendaciones.
1
1.1 Ubicación
A 803508 8530325
B 803639 8530325
C 803639 8530187
D 803508 8530187
2
CANAL
N INGRESO
CANAL
CANAL
LINEA DE
MEDIA TENSION
AREA DE
ESTUDIO
LINEA DE
MEDIA TENSION
El acceso a la ciudad de Urubamba, se realiza mediante la carretera principal asfaltada Cusco- Pisac-
Calca-Urubamba; con un recorrido aproximado de 65 Km y un tiempo aproximado de viaje en
camioneta de una hora 15 minutos. También se puede acceder, mediante la carretera principal
asfaltada Cusco-Chinchero-Urubamba; con un recorrido aproximado de 60 Km y un tiempo de viaje
de 1 hora.
El acceso a los diferentes barrios, urbanizaciones, asociaciones, áreas de expansión urbana, etc.,
que conforman la ciudad de Urubamba y alrededores, son mediante calles, avenidas, vías asfaltadas
y caminos de herraduras.
El área de estudio, está marcado por dos estaciones: una de estiaje entre los meses de mayo y
noviembre, y otra pluviosa entre los meses de diciembre a abril.
En la región quechua (2300 y 3500 m.s.n.m.), el clima dominante es templado, con notable diferencia
de temperatura entre el día y la noche. La temperatura media anual fluctúa entre 11 y 16 ºC; las
máximas entre 22 y 29 ºC y las mínimas entre 7 y -4 ºC durante el invierno, es decir, de mayo a
agosto. Las lluvias caen con regularidad durante el verano (diciembre a marzo). La vegetación típica
está conformada por: aliso, maíz, calabaza, caigua, tomate, trigo, árboles frutales, ciruelo, almendro,
peral, manzano, membrillo, durazno, etc.
En la región Suni (3500 y 4000 m.s.n.m.), el clima es seco y frío. La temperatura media anual fluctúa
entre 7 y 10 ºC, con máximas superiores a 20 ºC y mínimas invernales de -1 a -6 ºC (mayo-agosto).
La precipitación promedio es de 800 mm por año. La vegetación está compuesta por plantas
silvestres, como el quinual, quishuar, sauco, cantuta, motuy, carhuacasha, wiñayhuayna, suni, la
papa, año, quinua, cañihua, achis, tarwi, haba, oca y olluco.
En la región Puna (4000 y 4800 m.s.n.m.), el clima es frígido, con una temperatura media anual
superior a 0 ºC e inferior a 7 ºC. La precipitación fluctúa entre 400 y 100 mm al año. La vegetación
está compuesta de pajonales, ocsha, ichu, berro, totora, llacho, los bofedales, arbustos de culli,
árboles como la titánica, junco y cunco. Los productos alimenticios son la papa, cebada y maca.
3
1.4 Justificación del problema
la geotecnia clásica permite obtener parámetros del suelo hasta una profundidad somera sin embargo
por ser proyectos con mayor inversión se necesita mayor profundidad de investigación para un mejor
análisis sin alterar el ambiente asumiendo presupuestos que favorecen a la empresa constructora
asimismo haciendo un mejor diseño y eficiente para su construcción siendo estos estudios un
importante complemento en el mundo de la geotecnia.
1.6.1 Alcances
1.6.2 Limitaciones
• Las limitaciones del estudio fue la falta de realización de ensayos directos hasta la
actualidad para una mejor correlación geofísica.
1.7.1 Independientes
• Filtraciones de agua.
1.7.2 Dependientes
• Refracción sísmica.
1.7.3 Indicadores
1.8 Objetivos
5
• Obtener la variación de las velocidades de ondas de corte S (Vs) a través de los
diversos materiales.
1.9 Hipótesis
Los antecedentes de estudios con exploraciones geotécnicas son escasas, se realizaron estudios a
gran escala cubriendo grandes zonas en la provincia de Urubamba, estos estudios fueron cubiertas
con calicatas, ensayo PDC, y extracción de muestras de suelos fuera de la cobertura de nuestra área
de trabajo, pero es necesaria contar con información por la geología y geomorfología similar del lugar,
por consiguiente, tenemos los siguientes proyectos como referencia:
6
CAPITULO II
2.
7
2.1 Geomorfología
Las rocas que afloran en la Cordillera Oriental de la zona de estudio son metamórficas
del Paleozoico inferior, rocas vulcano-sedimentarias del Grupo Mitú (Permo-Triásico), e
Considerado como un valle interandino que tiene una dirección preferencial NOO-SEE.
El valle en la zona de estudio se halla a una altitud de 2630 a 2670 m.s.n.m. Las
vertientes son moderadamente suaves a empinadas. Sin embargo, en la vertiente
Suroeste se aprecian varios deslizamientos.
8
Se observa también que en el paso del río Vilcanota por la ciudad de Urubamba, solo
parte de la margen derecha del río se encuentra encauzada por muros de contención,
razón por la cual en la época de lluvias la erosión e inundación de sus márgenes se
incrementan. También en su margen izquierda se observan numerosos problemas
geodinámicos como deslizamientos, cárcavas y derrumbes.
Es una unidad que se encuentra al sur del área caracterizada por tener pendientes
abruptas entre 70–80%. En esta geoforma se dan los procesos de captación del agua
de lluvia, por lo que cumple la función de una cuenca de recepción del agua que fluirá
por las cárcavas.
2.1.2.2 Cárcavas
Son geoformas originadas por la erosión hídrica, por allí discurre el agua en épocas
lluviosas y está conformada por dos cárcavas principales en la parte baja que se
alimenta de una red de drenaje que la llamamos cárcavas secundarias. Éstas se
encuentran en número de 10 con anchos variables desde 0,50–7 m. Las cárcavas
principales alcanzan los 20 m, se han formado por erosión de los depósitos de flujos de
barro.
Son depósitos formados por erosión del flujo de barro antiguo I, esta erosión formó las
cárcavas y esos materiales fueron depositados al pie de éstas a manera de pequeños
conos.
9
• Terraza 1. Es la terraza más antigua formada por el río Vilcanota, se
encuentra ubicada en la parte alta, al margen derecho del río Vilcanota.
Actualmente, se encuentra una cantera en explotación.
2.2.1 Estratigrafía
En la zona de estudio afloran unidades litológicas, que van desde el Paleozoico inferior hasta el
Cuaternario y se distribuyen tal como se presenta en el Mapa N°02, de acuerdo a la
descripción siguiente:
10
yesos, lutitas y escasas calizas, se estima que tiene un grosor de 100 y 200 metros y
se le asigna una edad Albiano Medio (cretáceo inferior).
Son originados por la acción erosiva y de sedimentación del río Urubamba durante el
cuaternario. gran parte está cubierta por los depósitos de flujos de barro y en el sector
noreste se encuentran los depósitos fluviales, existe también una capa de suelo
orgánico negruzco de 20-40 centímetros de espesor y en otras zonas está cubierta por
los flujos de barro.
Son depósitos originados por la descomposición in situ de las rocas y por la gravedad.
Estos depósitos se reconocen en las laderas de pendiente fuerte de los cerros
circundantes a la ciudad de Urubamba, así como en las laderas de las quebradas
Ajok´asa y Mantacnainyoc.
11
2.2.1.4 Depósitos aluviales (Q-al)
Los conos aluviales están conformados por grandes bloques de granitos, cuarcitas,
rocas volcánicas, etc., envueltos por una matriz gravo-areno-arcillosa. Generalmente
tienen mala selección y estratificación y su permeabilidad es media a alta. Su
comportamiento mecánico es aceptable a bueno.
Foto 1. Se observan las dos cárcavas Foto 2. Se observa las siguientes geoformas:
principales que han erosionado los depósitos cerros escarpados, cárcava principal, conoide
de flujo de barro (Guadalupe, et al,2002). de depósito de flujos de barro, terraza fluvial
(Guadalupe, et al,2002).
12
2.3 Geomorfología local
A nivel local la ciudad de Urubamba se encuentra mayormente sobre depósitos aluviales y depósitos
fluviales. El piso de los valles Ajok´asa y Mantacnainyoc están compuestos por depósitos fluviales
desarrollados sobre depósitos aluviales antiguos.
El piso de los valles Ajok´asa y Mantacnainyoc están cortados por los ríos del mismo nombre, donde
los procesos erosivos se acentúan por el afloramiento de depósitos sedimentarios recientes.
Las laderas de los cerros están conformadas principalmente por depósitos coluviales o por
afloramientos de rocas sedimentarias, volcánicas, e intrusivas en muchos casos muy fracturados que
debido a las pendientes muy empinadas y otros factores forman conos aluviales o deslizamientos
principalmente.
2.4 Drenaje
La zona de estudio está situada en la margen derecha del río Vilcanota y a una altitud media de 2875
m.s.n.m., formando parte del denominado Valle Sagrado de los Incas a continuación se muestra los
principales afluentes hídricos.
Las aguas del Vilcanota en este tramo presentan un régimen torrencial, sobre todo en
temporadas de lluvia y de acuerdo a los mapas de peligros los desbordes probablemente
alcanzarían áreas urbanizadas en margen derecha del sector denominado la Florida. Aguas
abajo en el sector del centro vacacional no se han levantado edificaciones. El ancho del cauce
varía entre 80 y 95 metros.
Cuadro 2. Caudales máximos del río Vilcanota en m3/seg para diferentes períodos de retorno (INDECI,2005).
13
URUBAMBA
Figura 1. Sistema hidrológico del río Urubamba (Hydrology of the Urubamba River, by Francisco Olivera, Ph.D.; Center for
Research in Water Resources; University of Texas at Austin; January 1996).
14
N
AREA DE
ESTUDIO
UBICACION
FUNDAMENTOS TEÓRICOS
En este capítulo se describe los principales fundamentos teóricos de cada método para entender los
principios físico matemáticos, y a su vez hacer la correcta aplicación en investigaciones de suelos en
general.
3.
15
3.1 Ondas sísmicas
Las ondas se clasifican en dos categorías: viajeras y estacionarias. En las primeras hay propagación
de energía mientras que en las otras la energía asociada a la onda permanece confinada entre dos
fronteras (p. ej. Gettys, 1991).
Ondas sísmicas
Los parámetros que se usan para describir una onda son: la frecuencia, f = 1/T, y la frecuencia
angular, ω = 2π /T, donde T es el periodo; y el número de onda, k = 2π / λ, donde λ es la
longitud de la onda.
16
3.1.2 Ondas de cuerpo
Las ondas de cuerpo son ondas que pueden viajar a través de un Medio material elástico en
cualquier dirección. Como se mueven, las olas del mar pueden sufrir cambios en su nivel de
energía y en su geometría de trayectoria de viaje, sujeto a las propiedades físicas del medio.
4
k+ µ
3
𝑉𝑝 =
ρ
Las ondas P son las ondas elásticas más rápidas, 1.7 veces más rápida que las ondas
S y, por lo tanto, son detectados por primero en aplicaciones tales como pruebas
ultrasónicas y reflexión sísmica o refracción sísmica (Miller y Pursey 1955).
17
3.1.2.2 Ondas secundarías o de corte
µ
𝑉𝑠 =
ρ
La relación de la velocidad de la onda p y la onda s viene dada por: (Keary et al. 2002).
𝑉1 2(1 − 𝑣)
=
𝑉2 (1 − 2𝑣)
Por otra parte, cuando una onda S se propaga en un plano vertical, se conoce como
onda SV (Ver Figura 7) Cuando el plano de polarización es horizontal, es llamada onda
SH (Ver Figura 8).
18
Figura 7. Propagación de ondas SV (Alsadi, H, N,2017).
Como su nombre lo indica, las ondas superficiales son ondas que se mueven sobre la
superficie de la tierra.
Las principales características comunes son sus amplitudes relativamente grandes (alto
contenido de energía) y bajas frecuencias en comparación con las ondas de cuerpo. Además
de eso, ellos pueden moverse con una velocidad que es generalmente más lenta que las ondas
corporales y se mueven en el mismo medio.
Fue descubierto, en 1911, por un geofísico inglés llamado Augustus Edward Hough
Love (1863–1940). Se desarrolla solo en los casos en que un medio semi-infinito
elástico sólido está cubierto por una capa horizontal de baja velocidad, al igual que en
el modo de vibración de onda SH, el movimiento de partículas es transversal y se limita
al plano horizontal (Ver Figura 9), las ondas Love viajan por múltiples reflexiones entre
la parte superior y Límite inferior de los planos de la capa superficial.
19
Figura 9. Propagación de ondas Love (Alsadi, H, N,2017).
Es un método geofísico de gran utilidad, las primeras aplicaciones fueron principalmente para la
determinación de la profundidad del lecho de roca en el subsuelo desde 1959 aproximadamente, el
método también se ha utilizado con éxito para la ubicación de zonas de debilidad, tales como zonas
de corte y fallas (Sjögren et al., 1979).
20
Consiste en la medición de la variación de la velocidad de propagación de ondas sísmicas, las cuales
son producidas de forma artificial (impacto de comba con un soporte sobre superficie o utilización de
explosivos).
Esta energía producida por el impacto mecánico o explosivos genera un frente de onda la cual es
detectada por unos sensores llamados geófonos colocados a cada cierta distancia desde la fuente de
energía, por lo tanto, la distancia y el tiempo de llegada de las ondas permite calcular la velocidad de
propagación (ondas refractadas) de ondas primarias (ondas p).
Estas velocidades sísmicas nos permiten estimar la condición del material en el subsuelo o inferir el
tipo de suelo.
La ley de Snell se usa para calcular los ángulos de refracción y de reflexión de los rayos
sísmicos que inciden en un medio, es básicamente como tirar una bola de billar contra una
banda, rebota con determinado ángulo, según cómo se haya lanzado.
Entonces, supongamos dos medios diferentes, caracterizados por un índice de refracción (V1 y
V2). Este índice es dependiente del medio, ósea que es característico para cada medio (Ver
Figura 11).
Cuando un rayo sísmico se encuentra con un medio diferente al medio por el cual se está
propagando pasan 2 cosas. Si se topa con un medio por el que puede pasar, parte del rayo se
desvía (refracción) y parte de él rebota (reflexión) (Ver Figura 12).
21
Figura 12. Un rayo sísmico impacta el límite de dos medios con diferentes velocidades (v2 > v1) y se
divide en rayos sísmicos refractados y reflectados (Stefan J.,2010)
El Principio de Huygens establece que cada punto en un frente de onda actúa como la fuente
de una nueva onda que en un medio homogéneo genera un frente de onda esférico
secundario, cuya envolvente define la posición de una onda generada más tarde. El modelo de
propagación de ondas de Huygens requiere que los frentes de onda secundarios estén activos
solo en los puntos donde la envoltura toca sus superficies. La energía de la onda se está
extendiendo desde los puntos de fuente primarios en todas las direcciones, pero sus
interacciones mutuas hacen que la perturbación resultante cero en todas partes excepto en los
puntos donde tocan el sobre común, (Ver Figura 13).
Figura 13. Propagación de onda plana de acuerdo al Principio de Huygens. a) través de un medio
homogéneo donde la velocidad es constante. b) medio no homogéneo, de velocidad que está
cambiando uniformemente a través de la dirección de propagación (Alsadi, H, N,2017).
En 1650, Pierre Fermat descubrió un notable principio que puede expresarse en los términos
siguientes: “Un rayo luminoso que viaja de un punto a otro seguirá una trayectoria tal que,
22
comparada con trayectorias vecinas, el tiempo empleado en recorrerla es mínimo, máximo o
invariable (es decir, estacionario).”
Para poder dar una mejor interpretación en cuanto al cálculo de resultados se necesita el
conocimiento de las siguientes definiciones.
Es la representación gráfica de los tiempos empleados por las ondas en llegar a los
receptores a partir de la generación de ondas (punto de disparo), estos puntos de
disparo se dan fuera o dentro del tendido sísmico según el objetivo que se requiera.
Una vez construidas las curvas, se procede a la identificación de las secciones que
pertenecen a un mismo refractor, por medio del trazo de líneas que son ajustadas a los
puntos que pertenecen a cada curva t-x.
Así, cada línea recta define las diferentes velocidades como consecuencia de los
distintos estratos presentes en el subsuelo investigado. Los quiebres que se obtengan
al dibujar las líneas rectas, se usarán para determinar las profundidades de cada
estrato.
23
Las velocidades de un modelo simple de dos estratos, con la curva t-x; a continuación,
se calcula el tiempo de llegada del impulso refractado en un detector. Considerando la
ruta de viaje ABCD (Ver Figura 14):
𝑍1
𝐴𝐵 = 𝐶𝐷 =
cos 𝛼
𝐵𝐶 = 𝑋 − 2𝑍1 tan 𝛼
Figura 14. Descripción simple de dos capas, límites paralelos y la correspondiente curva de tiempo-
distancia (Bruce B. Redpath,1973).
𝐴𝐵 + 𝐶𝐷 𝐵𝐶
𝑇= +
𝑉1 𝑉2
24
2𝑍1 𝑋 − 2𝑍1 tan 𝛼
= +
𝑉1 cos 𝛼 𝑉2
1 sen 𝛼 𝑋
= 2𝑍1 − +
𝑉1 cos 𝛼 𝑉2 cos 𝛼 𝑉2
𝑉2 − 𝑉1 sen 𝛼 𝑋
= 2𝑍1 +
𝑉1 𝑉2 cos 𝛼 𝑉2
𝑉1
sen 𝛼 = (1)
𝑉2
1
𝑇 = 2𝑍1 𝑉1 sen 𝛼 − sen 𝛼 + 𝑋
𝑉1 𝑉2 cos 𝛼 𝑉2
1 − sen2 𝛼 𝑋
= 𝑍1 𝑉1 +
𝑉1 𝑉2 sen 𝛼 cos 𝛼 𝑉2
2𝑍1 cos2 𝛼 𝑋
= +
𝑉2 sen 𝛼 cos 𝛼 𝑉2
Y sustituyendo por :
2𝑍1 cos 𝛼 𝑋
𝑇= +
𝑉1 𝑉2
𝑇1 𝑉1
𝑍1 =
2 cos 𝛼
25
𝑉1 2
𝑉1
Por trigonometría cos 𝛼 = 1 − sen2 𝛼 ; sen 𝛼 = ; cos 𝛼 = 1−
𝑉2 𝑉2
𝑇1 𝑉1
𝑍1 =
𝑉2 2 − 𝑉1 2
2
𝑉2
𝑇1 𝑉1 𝑉2
𝑍1 =
2 𝑉2 2 − 𝑉1 2
Figura 15. Esquema de caso de múltiples capas y Curva de tiempo-distancia. (Bruce B. Redpath,1973).
𝑇1 𝑉1
𝑍1 =
2 cos 𝛼
26
𝑇𝑖2 𝑉1
𝑍1 =
𝑉1
2 cos sen−1
𝑉2
𝑉1
cos sen−1
𝑉3
𝑇𝑖3 − 𝑇𝑖2 𝑉 𝑉2
cos sen−1 𝑉1
2
𝑍2 =
𝑉2
2 cos sen−1
𝑉3
𝑉 𝑉
cos sen−1 𝑉1 2𝑍2 cos sen−1 𝑉2
4 4
𝑇𝑖4 − 𝑇𝑖2 𝑉 − 𝑉2 𝑉3
cos sen−1 𝑉1
2
𝑍3 =
𝑉3
2 cos sen−1
𝑉4
∆𝑇2 𝑉2
𝑍2 =
𝑉2 (2)
2 cos sen−1
𝑉3
∆𝑇3 𝑉3 (3)
𝑍3 =
𝑉3
2 cos sen−1
𝑉4
Donde ΔT2 y ΔT3; Son indicados en la Figura 15. Las ecuaciones (2) y (3) darán los
espesores son mayores que los reales, y se sugiere que los espesores se calculen
inicialmente de ambas maneras para saber si el error es significativo en una situación
particular.
Cuando el subsuelo tiene interfaces múltiples, lo que se conoce como el caso de capas
múltiples, siempre horizontales, la obtención de la ecuación de tiempo de viaje se
complica más, pero en definitiva tiene una forma similar a la de dos o tres capas.
27
Nuevamente la primera derivada de “t” respecto de “x” da la pendiente de la curva:
dt/dx = 1/Vn La construcción de las dromocronas desde los registros de campo
requiere de un proceso matemático conocido como regresión lineal. Como los puntos
no están alineados, debe encontrarse la recta que mejor se ajusta a ellos. Esto porque
que el suelo donde se “clava” cada geófono generalmente no es siempre el mismo, ni
tampoco el anclaje ni la fijación.
Si partimos de la ecuación:
2
𝑉1
2𝑍 1 −
𝑋 𝑉2
𝑇= +
𝑉2 𝑉1 𝑉2
Observamos que esta es la ecuación de una línea recta, pero con pendiente 1/V1 y
cuya ordenada en el origen se conoce como tiempo de intercepción (ti).
28
𝑋
𝑡=
𝑉1
𝑋𝐶 𝑋𝐶 2𝑍 𝑉2 2 − 𝑉1 2
= +
𝑉1 𝑉2 𝑉1 𝑉2
Y finalmente se obtiene:
𝑋𝐶 𝑉2 −𝑉1
𝑍=
2 𝑉2 +𝑉1
29
Figura 16. Ejemplo de interfaz de inmersión y conceptos de disparo inverso y velocidad aparente (Bruce
B. Redpath,1973).
𝑉1
𝑉2𝐷 =
sen 𝛼 + 𝛾
30
Donde γ es el ángulo de inmersión de la interfaz en relación con la superficie y α es el
ángulo de incidencia crítico. De manera similar, la velocidad aparente, V2U, observada
para el disparo en la dirección de actualización se da por:
𝑉1
𝑉2𝑈 =
sen 𝛼 − 𝛾
𝑉1
𝛼 + 𝛾 = sen−1
𝑉2𝐷
𝑉1
𝛼 − 𝛾 = sen−1
𝑉2𝑈
1 𝑉1 𝑉1
𝛾= sen−1 − sen−1
2 𝑉2𝐷 𝑉2𝑈
2𝑉2𝑈 𝑉2𝐷
𝑉2 = cos 𝛾
𝑉2𝑈 +𝑉2𝐷
Considerando como una parte importante en la aplicación del método de refracción sísmica,
pues la fuente de energía determinara la calidad de datos al momento de su procesamiento.
Entonces para tener una buena adquisición de datos como primera condición debemos tener
en cuenta que la cantidad de energía sea suficiente para alcanzar el objetivo propuesto, es
decir que la señal recibida sea lo suficientemente amplia y clara. Esto condiciona además la
frecuencia de la señal. Podemos resumirlos en tres grupos:
La caída de pesos es la más usada para investigación del subsuelo a poca profundidad,
digamos hasta los 50 metros. Se genera con una masa de 20 a 25 lb. tipo martillo, con la cual
se golpea una placa metálica afirmada al suelo para mayor efectividad de la transferencia de
31
energía. Cuando se necesita investigar profundidades mayores, digamos de 100 o 200 metros,
se utilizan masas de hasta 500 kg. suspendidas en un aparejo o en la caja de un camión o
camioneta, que se dejan caer pesadamente al suelo para generar la perturbación.
Figura 17. Se muestra a la izquierda, fuente de energía - caída de peso y a la derecha registro sísmico
(http://www.geofisicaexploraciones.com).
Los explosivos como la dinamita o los compuestos químicos como el Anfo (nitrato de amonio
con combustible) son los más usados para profundidades mayores, pero requieren que se los
coloque en perforaciones de unos metros (hasta 20 m) por seguridad y para que la mayor parte
de la energía se propague hacia adentro y no se pierda en la superficie.
Figura 18. Se muestra a la izquierda, fuente de energía – cartucho de 600 gramos de Pyroblast y a la derecha
registro sísmico (http://www.geofisicaexploraciones.com).
Los cañones de aire se utilizan cuando se realiza prospección sísmica en el agua, para
estudiar el fondo de océanos o lagos. Estos consisten de un gran pistón que recibe un golpe de
presión de aire que lo hace “saltar” repentinamente, generando una gran burbuja que se
propaga hasta tocar el fondo donde genera a vez la perturbación elástica deseada.
32
Debe mencionarse que cualquiera de las fuentes de energía descriptas, genera una
perturbación en el aire que también alcanza los sensores, entorpecen los registros y complica
la interpretación.
3.2.6 Sismograma
A cada punto de recepción le corresponde una traza. Las trazas poco agitadas antes del
impulso, señalan la llegada de las ondas directas y refractadas por una perturbación notable.
En ocasiones las llegadas no producen más que una pequeña agitación en las trazas, la cual
es fácilmente observable. El cálculo de los tiempos se realiza con ayuda de la traza que da el
instante inicial y un sistema de líneas normales a los bordes del sismograma.
El viento es un importante agente perturbador, tanto porque azota directamente a los geófonos
como porque sacude la vegetación; por este motivo se entierran los geófonos, lo que permite
mejorar el acoplamiento con el terreno. Si el sismógrafo lleva incorporada una pantalla que
permita seguir las oscilaciones de los receptores, el operador puede elegir un momento de
calma relativa para aplicar un impulso en el punto de disparo y obtener un registro más claro.
La utilización de la mayor cantidad de geófonos tiene por objeto reforzar las llegadas en fase,
es decir, la de las ondas sísmicas. Por otro lado, se pretende debilitar las llegadas que no están
en fase, es decir, las de las ondas parásitas. (Ver Figura 20).
33
Primeros arribos
(ondas P)
• Nos ayuda a determinar los parámetros pseudo - dinámicos con la ayuda del
método MASW.
• Este método es de los más baratos por metro de profundidad y por sondeo a su
vez es de mayor rapidez de ejecución
• Los datos que aportan reflejan las propiedades promedio de un volumen del
subsuelo, es decir, no son puntuales.
34
Dentro de las limitaciones del presente método, se describe a continuación los casos más
particulares:
• El clima es otro factor muy importante para tener en cuenta, ya que si se adquiere
datos con la presencia de lluvias esta emite un ruido por la caída de gotas de
agua sobre los geófonos ocasionando un mal registro de ondas p.
• Otro factor es el ruido del viento y la fuerza que este genere de tal forma que
pueda mover cables de conexión y geófonos.
En las últimas décadas, el uso de este método sísmico está aumentando para la caracterización
geotécnica debido a su rápida ejecución, simplicidad de implementación y rentabilidad. El método de
(MASW) es uno de los métodos sísmicos más interesantes.
Este método consiste en el análisis de dispersión de las ondas sísmicas superficiales (ondas
Rayleigh) para obtener la velocidad de fase a distintas frecuencias. Luego, a través de un proceso de
inversión no lineal, se busca un perfil unidimensional de velocidades de ondas de corte con respecto
a la profundidad.
Figura 21. Se muestra esquema de configuración de campo del método MASW (Park et al. 2004).
Las ondas superficiales son inherentemente dispersivas, lo que significa que la amplitud de la
onda superficial disminuye con la profundidad y la distancia desde la fuente. Dadas las
características dispersivas, se entiende que las ondas superficiales viajan exclusivamente
dentro de los suelos cercanos a la superficie. La velocidad de las ondas de Rayleigh es
36
comparable a las velocidades de las ondas de corte. En una formación rocosa con una relación
de Poisson de aproximadamente 0.25, la velocidad de la onda de Rayleigh es
aproximadamente el 92 por ciento de la velocidad de la onda de corte. En materiales con
relaciones de 0.4 a 0.5, el porcentaje aumenta a 94 a 95.5 por ciento, respectivamente.
Como regla general, se supone que la velocidad de las ondas Rayleigh es aproximadamente
del 90 al 92 por ciento de la onda de corte respectiva (Ivanov, Park y Xia 2009, Parasnis 1997).
La estimación de la velocidad de la onda de corte se puede hacer dentro de un margen de error
del diez por ciento utilizando estas suposiciones.
Las ondas Rayleigh solo existen si pueden viajar a lo largo de una superficie libre, en este caso
la interfaz tierra-aire. Se originan a partir de la interferencia de la componente vertical de las
ondas de corte (Sv) y las ondas P y tienen un movimiento de partículas elípticas. Debido a que
las ondas Rayleigh solo viajan a lo largo de superficies libres. Por otra parte, cuando se golpea
la superficie de un depósito de suelos, se generan ondas elásticas tanto de cuerpo (compresión
y corte) como de superficie (Love y principalmente tipo Rayleigh). El 67% de la energía
aplicada se propaga como ondas Rayleigh, el 27% como ondas de corte y cerca del 7% como
ondas de compresión. (Woods, 1968).
Figura 22. Se muestra en la parte superior la propagación de ondas Rayleigh (retrógrado) y en la parte inferior,
movimiento de partículas elípticas de la propagación de ondas Rayleigh (Van Der Hilst 2004).
37
3.3.2 Dispersión de ondas superficiales
𝜆=𝐶𝑓
Donde “c” es la velocidad de fase.
Es claro que en un medio estratificado; las longitudes de onda más grandes penetran a
profundidades mayores para un modo dado y su velocidad será afectada por las propiedades
del material a mayor profundidad, dado que las velocidades de las ondas superficiales son
representativas de la rigidez del material.
Figura 23. Dispersión de ondas de superficie en un medio heterogéneo (SAGEEP, 2003 & Lopes, I. et al., 2008).
Mediante el uso de ondas superficiales con un rango de longitudes de onda, es posible evaluar
las propiedades de los materiales a través de una gama de profundidades y obtener la curva de
dispersión.
38
Figura 24. Registro sísmico con su respectivo modelo inicial de velocidad de onda de acuerdo a las ondas
superficiales según las características del subsuelo (SAGEEP, 2003).
La velocidad de fase es la rapidez en la cual un punto de una fase de ondícula constante viaja
en dirección normal a la superficie de la onda en un medio homogéneo. Esta es la rapidez de
una onda plana que viaja en la misma dirección de tal medio. En un medio anisotrópico, la
velocidad de fase está sujeta a la dispersión angular, es decir, la rapidez varía con el ángulo.
En un medio atenuante, la velocidad de fase está sujeta a la dispersión de la frecuencia, es
decir, la rapidez varía con la frecuencia. (Winterstein, 1990).
39
El método es mucho mejor que la correlación cruzada, de modo que nos permite
calcular el directorio de velocidad de fase a partir de datos de forma de onda multicanal
(más de tres trazas) y separar el modo fundamental de la curva de velocidad de fase
de modos superiores y Ondas del cuerpo visualmente (Park et al. 1999), también
propusieron una transformación de forma de onda, denominada Análisis de ondas de
superficie multicanal (MASW, por sus siglas en inglés), que puede calcular la velocidad
de fase directamente desde una recopilación de disparo común multicanal similar al
método de McMechan. McMechan y Yedlin (1981) calculan la velocidad aparente (p)
en primer lugar y la transforman en el dominio de la frecuencia secundaria.
40
4. El valor absoluto se calcula y se representa en la velocidad de fase (c)
en función del dominio de la frecuencia.
Figura 25. Cálculo de la velocidad de la fase, esquema de métodos de transformación (fuente: Koichi
Hayashi,2003).
Las curvas de dispersión son un conjunto de curvas que representan la propagación de los
modos de onda que se encuentran en una geometría específica (Becker et.al,2017).
41
Por medio de las curvas de dispersión cada tipo particular de estructura laminar tiene su
conjunto asociado con curvas de dispersión, vale la pena mencionar que también son una
herramienta para comprender la propagación de la onda en la geometría estudiada. Estas
curvas son como un mapa que contiene todos los modos de onda posibles que pueden
propagarse en una onda guiada específica para cada combinación de frecuencia y número de
onda.
Figura 27. A la izquierda registro de forma de onda, incluyendo modo alto (Koichi Hayashi,2003), a la derecha una
ilustración de cómo se observaría en un modelo 3d (Ásta Ólafsdóttir,2014)
La generación de una curva de dispersión es uno de los pasos más críticos para generar
finalmente un perfil preciso de velocidad de onda de corte. Las curvas de dispersión
generalmente se muestran como velocidad de fase en función de la frecuencia, pero también
podrían presentarse en diferentes dominios como longitud de onda en función de la frecuencia
y velocidad de grupo en función de la frecuencia. Estas curvas permiten observar la dispersión
asociada a cada modo de onda.
42
función fase variará linealmente con la distancia. Cuanto mayor sea el número de receptores
utilizados por unidad de longitud en un ensayo, más alta es la resolución espacial de los
resultados, mejorando así la capacidad de detectar cambios en las propiedades horizontales
del medio.
Figura 28. Comparación de curvas de dispersión y su respectivo modelo unidimensional para materiales con
distintas características (Koichi Hayashi,2003).
Determinar una curva de dispersión confiable es el paso más importante del análisis y
ejecución del método MASW. El método asume que los registros están compuestos en
la mayoría por el modo fundamental de ondas superficiales e ignora los modos
superiores. Como los registros obtenidos contienen modos superiores y en algunos
casos con una energía no despreciable.
43
Transformado la función f(x,t) cuyo dominio es el desplazamiento - tiempo a una
función F(x,w) con dominio del desplazamiento – frecuencia.
El espectro de fase P(x,w) conserva todas las propiedades dispersivas de las ondas
generadas y el espectro de amplitud conserva las propiedades de atenuación de las
ondas. De esta manera se puede expresar F (x,w) también como:
Luego la función F (c,w) tendrá máximos a lo largo del dominio de cw, lo cual permitirá
graficar la curva de dispersión.
44
como 1/3 de la longitud de onda y se representa en el gráfico de velocidad-profundidad con su
velocidad de fase (Ver Figura 29(b)). Es bien sabido que la velocidad de fase de la onda
Rayleigh refleja aproximadamente la velocidad de onda S media hasta la profundidad de 1/2 a
1/4 de longitud de onda. Este procedimiento se denomina “transformación de longitud de onda”;
en términos de profundidad aparente (DA) y velocidad de onda de Rayleigh (VR).
45
Figura 29. Modelo inicial basado en 1/3 transformación de longitud de onda (Koichi Hayashi,2003).
46
modelo de capa dada y se compara con la curva de dispersión experimental obtenida
previamente (Ryden et al., 2006).
Figura 30. Modelo final o perfil unidimensional de ondas de corte – MASW (Koichi Hayashi,2003).
Choon B. Park y Mario Carnevale, 2010, en la publicación “Optimum MASW Survey - Revisit
after a Decade of Use”, indica que la profundidad de investigación (Zmax), adoptado como un
criterio empírico, es L/2. Y El análisis más preciso de las velocidades de fase puede lograrse
solo para longitudes de ondas de hasta una longitud extendida (λ ≤ L). Para un desplazamiento
mayor de la fuente (por ejemplo, X1 ≥ L) mejora la energía para ondas de superficie de longitud
de onda más larga y aumenta λmax con una extensión de receptor dada. Por otra parte, un
pequeño desplazamiento de la fuente (por ejemplo, X1 ≤ L) no sólo puede reducir λmax, sino
también dar lugar a velocidades de fase subestimadas. El grado de subestimación, sin
embargo, parece altamente dependiente del sitio y a veces insignificante.
Se debe evitar longitudes de arreglo muy extensas. Ya que las ondas superficiales generadas
por las fuentes sísmicas más comunes se atenúan por debajo del nivel de ruido al final de un
arreglo excesivamente largo, haciendo que la señal registrada por los receptores más lejanos
sea muy ruidosa impidiendo su utilización en el análisis (Park, 1999).
47
Figura 31. Máxima profundidad de investigación - MASW.
48
En cuanto a las limitaciones de este método tenemos:
• Sólo permite diferenciar las capas del subsuelo que presenten mayor contraste de
velocidad y mayor espesor. En el caso de suelos con capas intermedias de menor
velocidad el método si las puede visualizar, siempre y cuando estas capas sean
de regular espesor, no va a poder detectar lentes de material muy delgados.
Debido a la necesidad de conocer la naturaleza del subsuelo, durante las últimas décadas se han
venido desarrollando diversos métodos de exploración entre ellos los métodos de MASW y MAM, los
cuales nos permiten clasificar los suelos de acuerdo al código IBC 2012 y a la norma E.030 de Diseño
Sismorresistente 2016 en base a sus velocidades de ondas de corte, mediante el capítulo 2.3.1
perfiles de suelo, tabla Nro. 2 (Clasificación de los perfiles de suelo), donde identifican el tipo de perfil
del suelo de acuerdo a la velocidad promedio de ondas de corte de los 30 primeros metros.
49
para determinar la velocidad promedio de propagación de ondas de corte se determina con la
siguiente formula:
𝐻
VS30 =
ℎ𝑖
𝑉𝑠𝑖
Donde:
a continuación, se muestra los cuadros donde determinan la clasificación sísmica según las ondas de
corte en los primeros 30 metros de profundidad.
Perfil Descripción
S0 >1500 m/s - - Roca dura
S1 500 m/s a 1500 m/s > 50 > 100 kPa Roca o suelos muy rígidos
S2 180 m/s a 500 m/s 15 a 50 50 kPa a 100 kPa Suelos intermedios
S3 < 180 m/s < 15 25 kPa a 50 kPa Suelos blandos
S4 Clasificación basada en el EMS Condiciones excepcionales
50
3.5 Constantes elásticos dinámicos del suelo
El paso de las ondas “fuerza” la roca, induciendo una deformación proporcional al esfuerzo
según lo estableció Hooke, está dado por: ε = Eσ. La onda “P” provoca un cambio de volumen,
pero no de forma en respuesta al esfuerzo de compresión y tensión; mientras que la onda “S”
produce un cambio de forma, pero no de volumen (Ver Figura 32).
La propagación de las ondas en un medio geológico está afectada por las propiedades físicas
de las rocas; su velocidad es función directa de la elasticidad y densidad de la roca.
Las ondas “P” están afectadas por la “incompresibilidad” del material geológico, así como de su
rigidez.
Las ondas “S” producen un cambio de forma, pero no de volumen (Ver Figura 33), y solamente
dependen de la rigidez.
La velocidad a la cual cada tipo de onda viaja es función de la “elasticidad” del material
geológico, y representan un sistema sinergístico. Por ello, en función de las velocidades de
propagación –compresión y corte, la elasticidad del material se puede calcular; y si se añade el
conocimiento de los pesos volumétricos, es factible evaluar las constantes elásticas dinámicas.
deformación longitudinal.
La velocidad de propagación de las ondas sísmicas está relacionada con las propiedades
elásticas del medio de propagación; entonces, la Relación de Poisson se da con la siguiente
expresión matemática:
2
Vp
− 2
= s2
V
Vp
2. − 2
Vs
Es la relación entre esfuerzo unitario y la deformación longitudinal unitaria, esto es, para la
compresión y tensión respectivamente y se define de la siguiente forma:
Gd = .Vs
2
E d = 2.Gd (1 + )
52
3.5.5 Modulo volumétrico o Bulk
Ed
Kd =
3.(1 − 2 )
Con las mediciones adquiridas se construye una sección en dos dimensiones (2D) que muestra una
primera aproximación de los cambios en el subsuelo. Posteriormente, se aplica un algoritmo de
inversión para obtener la distribución real de resistividades o imagen eléctrica. Imagen que será un
resultado interpretable desde un punto de vista físico y geológico, y que dará información sobre las
características físicas del subsuelo.
Este método geofísico fue desarrollado para la investigación de áreas de geología compleja (Griffiths,
1993). Su aplicación es diversa, como en la exploración de recursos naturales y actualmente en
numerosas aplicaciones ambientales, hidrogeológicas, arqueológicas y geotécnicas.
También se destaca por ser una herramienta de carácter no destructivo, con un bajo costo, es de fácil
implementación, adecuado por optimización de tiempo, combinando adecuadamente resolución
lateral y profundidad de investigación.
Figura 34. Recopilación de datos en 2D mediante el sistema de resistividad de múltiples electrodos (Loke et al.,
2013).
53
3.6.1 Ley de ohm y resistividad
ρ = RS/L
V = IR
Esta relación se representa mejor imaginando la corriente que pasa a través de un cable
delgado. La aplicación expuesta de la Ley de Ohm ha hecho de esta relación un concepto
clave en el estudio de la teoría eléctrica (Gibson y George 2003). Las unidades para potencial
eléctrico, corriente y resistencia son voltios, amperios y ohmios, respectivamente.
54
En este caso la corriente se da por circulación de iones, y la conductividad es función de la
cantidad de agua y de sales ionizadas disueltas en ella. Concretamente, dependerá de los
siguientes factores:
ρa =𝐾 Δ𝑉𝐼
Para el caso de un medio homogéneo la expresión nos arrojará el valor de la resistividad real.
Es decir, si el terreno fuera homogéneo, la resistividad que se obtiene a partir de la expresión,
tendría que ser constante e independiente del espaciamiento que existe entre los electrodos
(Parasnis, 1997).
55
3.6.4 Propiedades eléctricas de los materiales terrestres
La resistividad es una de las magnitudes físicas con mayor amplitud de variación para diversos
materiales. Además, su valor depende de diversos factores como la temperatura, humedad o
presión. Estrictamente hablando todos los cuerpos son eléctricamente conductores dado que
permiten, en mayor o menor medida, el paso de portadores de cargas eléctricas.
La corriente eléctrica puede ser propagada a través de los minerales y las rocas, a través de
tres formas: conducción electrónica (óhmica), conducción electrolítica y conducción dieléctrica.
- Mientras que la conducción dieléctrica se lleva a cabo en conductores pobres o aislantes, los
cuales tienen muy pocas o más bien, ninguna carga libre.
Cuadro 5 Valores de resistividad representativos para diferentes materiales (modificada de Gasulla, 1999).
Resistividad (Ωm)
Material
Mínimo Máximo
Grava 50(agua saturada) > 10^4 (en seco)
Arena 50(agua saturada) > 10^4 (en seco)
Limo 20 50
Marga 30 100
Arcilla (húmeda) 5 30
Arcilla (seca) - >1000
Lodos 15 25
Arenisca <50(húmedo) > 10^5 (en compacto)
Caliza <100(húmedo) > 10^5 (en compacto)
Esquisto <50(húmedo) > 10^5 (en compacto)
Toba 2x10^3 (húmedo) 10^5 (mojada)
<100(degradado,
Rocas ígneas y metamórficas > 10^6 (en compacto)
húmedo)
Roca de sal 30(húmedo) > 10^6 (en compacto)
Residuos domésticos e industriales <1 >1000 (de plástico)
Agua natural 10 300
56
Resistividad (Ωm)
Material
Mínimo Máximo
Agua de mar (35%. NaCI) 0.25 -
Agua salada (salmuera)<0.15 3x10^3
Agua superficial (rocas ígneas) 0.1 100
Agua superficial(sedimentos) 10 10^6
El agua pura es muy poco conductora a causa de su muy reducida disociación. La resistividad
del agua destilada es de unos 105 Ωm por lo que puede considerarse como aislante. Las
aguas que se encuentran en la naturaleza presentan, sin embargo, conductividad apreciable,
pues siempre tienen disuelta alguna sal, generalmente NaCl. Así las aguas de lagos y arroyos
de alta montaña varían entre 103 Ωm y 3,103 Ωm., las aguas subterráneas tienen
resistividades de 1 a 20 Ωm, y las aguas marinas tienen una resistividad de unos 0,2 Ωm.
Figura 36. Valores de resistividad representativos para rocas y minerales (Palacky, 1987).
57
La resistividad de las rocas puede variar en margen amplísimo en función del contenido en
agua, de la salinidad de esta y del modo de distribución de los poros. (Ver Figura 36) presenta
un gráfico de los márgenes de variación más comunes en algunas rocas y minerales. La
fisuración, impregnación en agua salada, etc., pueden extender estos límites.
Figura 37. Dispositivo Tetraelectródico para la medida de la resistividad del suelo. (Moller et al., 1995).
De acuerdo con Parasnis (1997), la resistividad aparente obtenida a partir de una medición con
un dispositivo Electródico, no es un promedio de las resistividades obtenidas del subsuelo, sino
que cada valor de resistividad es una medición individual que está asociada al punto atribuido
al centro del dispositivo (dependiendo de la geometría de éste) y está en función del
espaciamiento de los electrodos y la distribución real de resistividad.
Para el caso más general, en donde los 4 electrodos están dispuestos de modo cualesquiera,
en línea recta y sobre una superficie plana, la diferencia de potencial Δ𝑉 se podrá obtener
mediante la siguiente expresión.
en donde AM, BM, AN y BN son las distancias entre cada par de electrodos. Por lo tanto, si
también tenemos que:
58
Entonces con base en las expresiones anteriores, podemos definir a la resistividad aparente de
un medio de la siguiente manera:
59
3.6.5.2 Dispositivo Wenner
En este dispositivo los electrodos AMNB son equidistantes entre sí, a una distancia a.
En la figura se observa esquemáticamente la disposición de los electrodos.
En este dispositivo los electrodos se disponen sobre una línea en el orden ABMN
formando así un doble dipolo.
60
Figura 41. Esquema Electródico del dispositivo Dipolo - dipolo.
Para cada tipo de dispositivo electródico existen ventajas y desventajas, que depende de la
profundidad de investigación, sensibilidad del dispositivo, cobertura de investigación e
intensidad de la señal, al establecer el panorama de trabajo se debe tomar en cuenta estos
criterios, algunos tomados en el Cuadro 6.
61
Cuadro 6. Ventajas y desventajas en diferentes tipos de dispositivos electródicos (Summer, 1976; Koefoed, 1979;
Ward, 1990; Parasnis, 1996).
Relación señal/
Dispositivo Electródico Ventajas Desventajas
ruido
Arreglo simétrico.
Resolución vertical
Minimiza el acoplamiento
(estructuras horizontales)
inductivo entre electrodos de
pobre.
corriente y potencial. Pobre (intensidad
Dipolo-dipolo Sensibilidad al ruido geológico.
La profundidad de penetración de la señal débil)
La resistencia de contacto de
no depende de la geometría,
los electrodos móviles es más
sino del contraste de
crítica que el polo-polo.
resistividades.
3.6.7 Pseudosección
El conjunto de datos de resistividad aparente (ρa) son utilizados para la construcción de una
pseudosección de (ρa) (Hallof, 1957) los cuales muestran una representación aproximada de la
variación de la resistividad en el subsuelo la cual involucra datos de ρ transversales como una
sección a profundidad, donde cada valor de (ρa) es graficado como si fuera la ρ verdadera de
un punto inmediatamente bajo el centro del arreglo de electrodos en una propiedad
proporcional al espaciamiento de los electrodos a (Edwards, 1977). Para obtener una
62
distribución confiable de la distribución de resistividades en el subsuelo, es necesario aplicar
una inversión a la pseudosección de resistividades aparentes.
La profundidad de investigación puede tener una variedad de definiciones físicas, pero todas
estas se derivan de la necesidad de relacionar una medición hecha en la superficie con alguna
profundidad particular, además de que se puedan optimizar los parámetros del estudio para
una mejor identificación del objetivo (Barker, 1989).
El Cuadro 7, muestra los factores para calcular la profundidad media de investigación (Ze) de
las diferentes modalidades. L es la longitud total del arreglo. Note que existen valores idénticos
para las modalidades Wenner-Schlumberger y polo-dipolo. El factor geométrico esta dado para
un espaciamiento entre electrodos “a” de un metro.
63
Para determinar la profundidad máxima mapeada por un estudio en particular, se multiplica el
máximo espaciamiento entre electrodos “a”, o la máxima longitud del arreglo “L”, por el factor
apropiado de nivel de profundidad dado en el Cuadro 7 (Loke, 2004).
Por ejemplo:
Si el máximo espaciamiento entre electrodos “a” usado por el arreglo Wenner es de 100 metros
(o la L máxima es de 300 metros), entonces la profundidad máxima mapeada es alrededor de
51 metros.
Por ejemplo:
Si se requiere calcular la profundidad máxima para el arreglo dipolo-dipolo para una separación
entre electrodos a = 5 metros y una n =7, la “L” máxima será:
Esto multiplicado por el factor = 0.216 para n = 7 (Ver Cuadro 7) resulta una profundidad
máxima de penetración aproximada de 10 metros aproximadamente.
Cuadro 7. Profundidad media de investigación (Ze/a y Ze/L) para las diferentes modalidades. El factor geométrico
está calculado para el valor de “a” igual a 1 metro (Loke, 2004).
Dispositivo Factor
nivel Z e/a Z e/L
Electródico Geométrico
Wenner - 0.519 0.173 6.2832
n=1 0.416 0.139 18.85
n=2 0.697 0.174 75.398
n=3 0.962 0.192 188.5
n=4 1.22 0.203 376.99
Dipolo-dipolo
n=5 1.476 0.211 659.73
n=6 1.73 0.216 1055.6
n=7 1.983 0.22 1583.4
n=8 2.236 0.224 2261.9
64
Dispositivo Factor
nivel Z e/a Z e/L
Electródico Geométrico
n=3 1.318 0.189 37.699
n=4 1.706 0.19 62.832
n=5 2.093 0.19 94.248
n=6 2.478 0.191 131.95
n=7 2.863 0.191 175.93
n=8 3.247 0.191 226.19
n=9 3.632 0.191 282.74
n=10 4.015 0.191 345.58
n=1 0.519 - 12.566
n=2 0.925 - 37.699
n=3 1.318 - 75.398
n=4 1.706 - 125.66
Polo-dipolo
n=5 2.093 - 188.5
n=6 2.478 - 263.89
n=7 2.863 - 351.86
n=8 3.247 - 452.39
Polo-polo - 0.867 - 6.28319
Por otra parte, también tenemos un gráfico con la relación de la profundidad de investigación
en base a los niveles de medición para cada dispositivo Electródico.
Figura 44. Profundidades medias de investigación sobre un medio homogéneo de resistividad para los
dispositivos obtenidas para una distancia interelectródica de 5 metros y distintos niveles de medición
(Edwards 1977). Valores idénticos para los dispositivos Polo-Dipolo y Wenner–Schlumberger.
65
3.6.9 Inversión de datos
los valores aparentes se pueden usar de manera efectiva para estimar las condiciones reales
del subsuelo. Las secciones aparentes de resistividad proporcionan los datos necesarios para
completar un proceso de inversión iterativo.
El procedimiento de generación del perfil real mediante su inversión, parte de los datos de
resistividad aparente medidos, representados en forma de pseudo-perfil. A continuación, se
genera un modelo hipotético de resistividades verdaderas del subsuelo, y resolviendo lo que se
define como “cálculo del problema directo”, se llega al modelo de resistividades aparentes que
se derivaría del mismo. Estas resistividades aparentes se comparan con las realmente
medidas, y se calcula el error cometido. A través de este error se modifica el modelo hipotético
real de resistividades verdaderas, y se repite el proceso anterior. De este modo, tras una serie
de iteraciones, se consigue un modelo de resistividades verdaderas del subsuelo que da
explicación a las resistividades aparentes medidas.
66
especialista ya que es necesario configurar el programa de cálculo de acuerdo a las
características del objeto de estudio, el dispositivo y la experiencia del mismo.
Figura 46. Corrección de topografía para inversión de datos-TGE (Phil Sirles, 2016).
67
3.6.10 Ventajas y limitación del método
Se pueden realizar un gran número de medidas, dando como resultado una alta resolución
tanto longitudinalmente como a profundidad (perfil 2D).
Este método nos brinda muchas ventajas las cuales serán descritas a continuación:
• Los resultados del estudio previo a la realización del informe se pueden estimar en
campo con rapidez (computadora portátil). Este procesamiento, casi en tiempo
real, le confiere al método una gran versatilidad a la hora de resolver posibles
problemas surgidos durante toda la logística y adquisición de datos en campo,
ayudando a tomar decisiones en cuanto a la realización de más líneas como
complemente para la obtención de mejores resultados y discretizar mejor el
subsuelo.
68
CAPITULO IV
En este capítulo se describe la instrumentación geofísica que se utilizó para la adquisición de datos
para cada método aplicado en este proyecto, así mismo se describe la metodología de campo y el
procesado de datos para su realización.
4.
69
4.1 Refracción Sísmica – Análisis multicanal de ondas superficiales (MASW)
70
4.1.1.2 Accesorios de sismógrafo geode
Cable
digital
Conexión de
cable conector
de interfaz
71
software tiene que estar registrado, se envía el código de acceso del
usuario, y Geometrics proporciona un número de registro.
Cable conector
geode – caja
NIB
Caja NIB
4.1.1.3 Geófono
Para levantamientos de poca profundidad, las señales sísmicas tienden a registrar una
frecuencia mucho más alta y un geófono con una frecuencia natural de alrededor de
4.5 - 14 Hz es un buen arreglo.
Cable de geófono:
72
4.1.1.4 Cable sísmico o cable conector de geófono
Es el cable que nos permite la conexión de todos los geófonos y esta a su vez va
conectado con el sismógrafo
Este cable es el que transporta las señales eléctricas transformadas por los geófonos.
Además, tiene puntos de conexión eléctrica donde se conecta cada receptor a
intervalos uniformes. El espaciamiento entre geófonos varía desde 1 m o mayores,
todo en función del nivel de detalle definido por los objetivos de la investigación.
73
4.1.1.6 Fuente de energía
Las fuentes de energía en el presente proyecto provienen del uso de una comba de
25lb, esta herramienta es sincronizada con el sismógrafo mediante un “cable
interruptor” de sensibilidad al impacto, unido al mango de la comba y conectado
directamente al conector de activación en el sismógrafo. Puede usarse un cable de
extensión, para permitir que la comba este localizada lejos del sismógrafo.
Una placa de aluminio o acero de 15cm x 15cm se coloca en el terreno para recibir el
impacto de la comba. Esta “placa de percusión” proporciona un acoplamiento más
eficiente y una activación más precisa.
Esta fuente de energía es muy común para levantamientos poco profundos. Es muy
popular debido a su poco peso, portátil, de bajo costo, repetible y seguro. Sus únicas
serias limitaciones son su poco grado de profundidad, aunque en este respecto, se
compara favorablemente con un numero de aditamentos menos portátiles y más
costosos. Los sismógrafos de mejoramiento de señal o de apilamiento fueron
desarrollados para usarse con una comba.
Los parámetros utilizados para el método de refracción sísmica y MASW son casi semejantes
en cuanto al tendido de línea sin embargo se diferencian la longitud de registro y el intervalo de
muestreo.
Cuadro 8. Parámetros de adquisición de datos sísmicos para MASW y Refracción Sísmica (Geometrics Inc,
2009).
Parámetro Descripción
Arreglo lineal
74
Parámetro Descripción
MASW:
0.5 o 1 milisegundos (ms)
Intervalo de muestreo
Refracción sísmica:
0.125 – 0.5 milisegundos (ms)
MASW:
1 a 2 segundos.
Longitud de registro
Refracción sísmica:
0.250 – 1 segundos.
Según sea necesario para mejorar la calidad de los datos, espere a
Stacking
momentos callados para registrar.
75
Cuadro 9. Coordenadas de Ubicación de líneas de refracción sísmica.
El área de trabajo constituye parte de un material cuaternario por los cual se podría
hacer una distribución uniforme cubriendo toda el área de trabajo.
76
4.1.3.2 Tendido de línea de refracción sísmica
Para el tendido de línea sísmica (armado de arreglo sísmico) se hizo el uso de una
cinta métrica de 100m, tendiendo de forma lineal cubriendo una longitud de 46m para
algunas líneas y para otras una longitud de 94 m, seguido se empezó a la realización
del plantado de geófonos(14Hz) perpendicularmente a la superficie de terreno con
intervalos de espacio de 2m entre cada geófono, de tal manera que queden bien
acoplados al suelo como si fuera parte del ello.
Una vez ya plantado los geófonos se procedió al tendido de cable sísmico respetando
los intervalos de 2m para cada punto de conexión, al mismo tiempo se hace el tendido
respectivo del cable trigger.
posteriormente se hizo la conexión de los geófonos con el cable sísmico, esta parte del
procedimiento es importante ya que uno de los problemas que se presenta al momento
de hacer la adquisición de datos se genera por la mala conexión entre estos dos
componentes, por eso es necesario limpiar bien los conectores de los geófonos o las
aberturas de conexión del cable sísmico , ya que por razones del mismo ámbito laboral
estos están cubiertos de tierra como consecuencia genera un mal registro con respecto
a ese canal.
77
Figura 51. Armado de línea sísmica (Fuente propia)
Una vez culminado el armado de la línea sísmica se verifico las condiciones de cada
geófono en su conexión verificando si estaba listo para el registro de ondas P.
La fuente para generar las ondas sísmicas fue el impacto de un martillo de 25 libras (o
11.3 kg) sobre una placa metálica.
1° Shot 2° Shot
Figura 52. Se muestra la distribución de puntos de disparo, el primero se ejecuta antes del 1er geófono
y a una distancia igual a la separación entre geófonos, y el segundo disparo entre el 6to y 7mo
geófono.
3° Shot
Figura 53. El tercer punto de disparo se realizó entre el 12vo y 13vo geófono.
78
4° Shot 5° Shot
Figura 54. El 4to punto de disparo se da entre los geófonos 18 y 19, y El ultimo de, se ejecuta después
del 24vo geófono a una distancia igual a la separación entre geófonos.
Durante la adquisición de datos se observó las tendencias en los datos ya que esto
ayuda a hacer selecciones fiables de los tiempos de llegada, de no hacer un buen
registro se mejora la señal golpeando varias veces con la comba. Esto es útil en áreas
que son sísmicamente ruidoso y en áreas con condiciones geológicas complejas.
Por otra parte, para la realización del traslape de línea sísmica se hizo la siguiente
configuración del tendido de cable en cuanto a la adquisición de datos.
Figura 55. Se muestra el ultimo de punto de disparo, será el primer punto de disparo de la segunda línea
sísmica.
Esto nos permite la unión de dos líneas para tener una mayor longitud de estudio en el
área con la profundidad de investigación constante y con la misma resolución.
Para realizar procesamiento de las líneas de refracción sísmica se usó como ayuda el
paquete de procesamiento “SeisImager”.
Este paquete está conformado por cuatro programas los cuales son Pickwin (para
análisis de primeros arribos de ondas), Plotrefa (para análisis de refracción sísmica),
GeoPlot (para visualización de data sísmica), y Surface Wave Análysis Wisard.
El primer programa permite realizar las lecturas de los sismogramas a partir de eso el
procesamiento consiste en seleccionar los primeros arribos de las ondas “P”,
79
Figura 56. Registros de ondas compresionales - LS01.
80
Figura 58. Perfil de refracción sísmica (Fuente propia).
Como último paso se hace una técnica de inversión de cuadrícula para determinar la
velocidad de bloques de 2 dimensiones individuales (píxeles) dentro de un perfil en
lugar de modelar velocidades como capas. Como resultado, la tomografía por fracción
puede, en algunos casos, modelar con mayor precisión y proporcionar una mejor
resolución de la estructura de velocidad compleja del subsuelo. Una limitación de la
fracción sísmica es la incapacidad de referirse a la existencia de ciertas capas,
denominadas capas ocultas o zonas ciegas. Esto se debe a un contraste de velocidad
insuficiente del espesor de la capa.
Figura 59. Perfil de refracción sísmica con proceso de inversión (tomografía sísmica), (Fuente propia).
81
Cuadro 10. Velocidad de la onda de compresión Vp. Según la Norma ASTM – D 5777
Descripción Vp (m/s)
El trabajo de campo consistió en la adquisición de 8 líneas de MASW, las cuales cada uno
tenía una longitud de 57.5m con intervalos de distancia de 2.5m.
Coordenadas UTM
Tipo de
N° MASW WGS 84
onda
Este Norte
1 MASW 01 S 803557 8530203
2 MASW 02 S 803557 8530231
3 MASW 03 S 803563 8530228
4 MASW 04 S 803571 8530227
5 MASW 05 S 803576 8530277
6 MASW 06 S 803570 8530278
7 MASW 07 S 803562 8530278
8 MASW 08 S 803566 8530254
82
4.1.4.1 Tendido de línea sísmica para el método MASW
Para la medición de los tiempos de viaje de las ondas sísmicas de bajas frecuencias y
gran amplitud, se empleó un arreglo de 24 geófonos de 4.5 Hz. de frecuencia, (Sin
embargo, según Park et al 1999, el ensayo MASW también puede ser eficaz con solo
doce canales de grabación conectados a geófonos de baja frecuencia menor a 10 Hz).
La fuente para generar las ondas sísmicas es el mismo del método de MASW (impacto
de un martillo de 25 libras (ó 11.3 kg) sobre una placa metálica). Se realizaron cuatro
puntos de impacto ubicados a 0.2L y 0.4L del extremo de la línea, tal y como se
muestra en la siguiente figura:
Geófonos
Esta transformación elimina todo el ruido cultural ambiental, así como el ruido
generado por la fuente tal como ondas dispersadas de objetos enterrados (cimientos
del edificio, alcantarillas, rocas, etc.). La imagen de dispersión muestra la relación entre
la velocidad de fase y la frecuencia de las ondas propagadas horizontal y directamente
desde el punto de impacto a la línea receptora.
La curva de dispersión extraída se utiliza finalmente como una referencia para volver a
calcular la variación Vs con profundidad por debajo del área de investigación. Este
retrocálculo se denomina inversión y el proceso puede ser automatizado con
restricciones razonables (Xia et al., 1999).
- Imágenes de dispersión.
84
Para el procesamiento se utilizó los programas Surface Wave Analysis Wizard y Wave
Eq (Surface Wave Analysis).
El primer programa (Surface Wave Análisis Wizard) nos muestra el gráfico de distancia
versus tiempo que luego mediante las transformadas de Fourier se obtiene las curvas
de velocidad de fase versus frecuencia en donde se puede ver la tendencia de la onda
de fase que define la velocidad de la onda S que nos permite obtener la curva de
dispersión en el modo fundamental generada en el estudio.
a) b)
Figura 61. a) Tiempo de llegada de la onda de corte para diferentes distancias. b) Curva de dispersión:
velocidad de fase vs modo fundamental (Fuente propia).
85
Figura 62. Modulo fundamental de la curva de dispersión- MASW (Fuente propia)
Figura 63. Perfil unidimensional de ondas de corte (Vs)- MASW (Fuente propia)
86
4.2 Tomografía Geoeléctrica
El método involucra contar con equipos y herramientas para: generar la emisión medida de
corriente continua en el terreno y en la medición con un equipo procesador de la diferencia de
potencial existente en dos puntos fijos, según un esquema geométrico determinado de cuatro
electrodos.
4.2.1.1 Resistivímetro
Características Descripción
Trasmisor (Rx)
Max. Voltaje. 800V en modo switch.
Max. Voltaje. 1000V in modo manual.
Max. Intensidad de corriente. 2500mA.
87
Características Descripción
44800 lecturas, hasta 8 horas de forma de onda
Almacenamiento completa, almacenadas en la memoria de estado
sólido.
Dimensiones
4.2.1.3 Cables
Este cable tiene la finalidad de conducir la corriente eléctrica inyectada por el sistema
transmisor del Resistivímetro Syscal Pro hacia el subsuelo mediante la configuración
realizada en su arreglo Electródico con los respectivos electrodos.
88
El cable utilizado en este estudio es un cable #16 fabricada por la empresa INDECO.
- 4 combas de 8lb.
- Batería de 12v
- Navaja de campo.
- 2 pelacables.
- Cuaderno de campo.
Cuadro 13. Parámetros de adquisición del método de tomografía geoeléctrica (SYSCAL Pro – User’s manual)
Parámetro Descripción
Configuración Dipolo – Dipolo
medición Rho & Ip
89
Parámetro Descripción
Quality factor 3
Voltaje min. /max. (Tx) 500/800v
Numero de Mediciones por Lectura 2-3
Display options Average
voltaje Unsigned
Duración del pulso de inyección 2s
El dispositivo Electródico empleado en esta ocasión fue el dispositivo Dipolo - dipolo debido es
muy sensible a cambios de resistividad horizontales.
Primeramente, se hizo el tendido de una cinta métrica para poder hacer la medición de
los intervalos de cada electrodo.
90
Figura 66. Esquema de tendido de línea de Tomografía geoeléctrica.
Los electrodos fueron ubicados perpendicular al terreno, para esto se hace el uso de
las combas para clavar sobre el suelo de tal forma que queden firmes, lo suficiente
como para garantizar un buen contacto electrodo-suelo.
Una vez instalado todos los cables con sus respectivos electrodos se hace la
instalación del Resistivímetro tomando una posición estratégica para visualizar la línea.
91
Foto 10. Instalación de cable en el Resistivímetro (Fuente propia).
Ahora, a través del Resistivímetro configuraremos todas las variables del dispositivo:
número de electrodos utilizados, distancia entre ellos, dispositivo Electródico de
medida, así como el número de medidas que queremos realizar. (Ver Cuadro 13).
92
permanentes durante las lecturas, de tal forma que la adquisición de datos sea buena
en campo ya que es lo más importante en el estudio.
Una vez verificado todos estos aspectos durante la adquisición de datos, se guarda la
lectura.
Figura 68. Esquema de adquisición mediante Tomografía Geoeléctrica (Modificada de Chávez et al.,
2011).
Una vez terminado el registro de todo el perfil, se ingresaron los datos al programa
Oasis Montaj para el control de calidad de datos de resistividades aparentes (Ver
Figura 69).
93
Figura 69. Evaluación de data cruda, Pseudosección - TG 01 (Fuente propia).
En general los cambios de resistividad aparente eran suaves dentro de un mismo nivel
horizontal en la pseudo sección. Esto podría entenderse como una afirmación o
confirmación de las anomalías observadas en la otra línea geoeléctrica - TG02.
94
De manera iterativa se van ajustando los valores de resistividad de cada bloque hasta
conseguir un error mínimo o alcanzar un nivel máximo de iteraciones o lograr una
pequeña diferencia entre errores para iteraciones sucesivas (Ver Figura 71).
El programa utiliza una distribución fija de celdas y en cada iteración cambia los
valores de resistividad solamente. La rutina de inversión está basada en el método de
mínimos cuadrados con suavizado forzado (Sasaki 1992), (Ver Figura 72).
por otra parte, una vez hecho la inversión se obtiene el perfil geoeléctrico real se podría
decir, para tener una mejor visualización e interpretación, se exporto la data del perfil
en un formato “XYZ” al programa Oasis Montaj para su presentación final.
95
4.3 Levantamiento fotométrico con Drone
En la actualidad la tecnología de drones está transformando muchos campos en los que la fotografía
aérea es una alternativa para medición y procesamiento, en este caso mucho más económico en
tiempo y costo para obtener información topográfica, mediante imágenes y aplicando técnicas
fotogramétricas.
El drone o denominado vehículo aéreo no tripulado (UAV), es un equipo móvil aéreo autónomo
que posee características técnicas particulares para realizar trabajos de todo tipo, se pueden
utilizar tanto en el exterior como en el interior de algunos espacios, se le pueden incorporar
como carga útil diferentes tipos de sensores ópticos: cámaras de video, cámaras fotográficas
de espectro visible, infrarrojo, etc., e incluso otro tipo de sensores. Puede desarrollar diferentes
alturas de vuelo (desde 1 m. hasta 500 m.) en vuelos totalmente automáticos programados y
controlados de modo remoto o de forma manual mediante un mando de radiocontrol.
96
con mayor claridad y resolución y menos ruido, la capacidad de doble frecuencia del
control remoto aumenta la eficacia y la estabilidad de la transmisión de video HD.
El Phantom 4 introduce una serie de características que se supone lo hacen más fácil
de volar.
Figura 74. Diagrama de drone Phantom 4 pro (Manual Phantom 4 pro, 2017).
Figura 75. Diagrama de control remoto Phantom 4 pro (Manual Phantom 4 pro, 2017).
97
4.3.1.2 Características principales
En este caso el vuelo de drone tuvo el fin de visualizar toda el área de trabajo, así
como generar la topografía, pero por factor tiempo no se pudo utilizar alguna imagen
georreferenciada para la delimitación de la zona de vuelo.
Las imágenes se obtienen según el plan de vuelo planificado, entonces se registra una
imagen en cada toma de foto manual permitiéndonos obtener fotografías individuales
para la posterior obtención de orto mosaicos, curvas de nivel, nubes de puntos entre
otros productos.
Foto 12. Toma de fotos desde la ubicación de manejo de drone (Fuente propia).
99
4.3.2.4 Procesamiento de imágenes fotográficas
Una vez tomadas las fotografías con sus parámetros se realiza el cálculo de los
parámetros de orientación externa de cada una de las fotografías.
Para este trabajo, el procesamiento de las imágenes se realizó con el programa Agisoft
Photoscan que es un software destinado a crear modelos 3D de alta calidad a partir de
imágenes tomadas desde puntos de vista de cámaras conocidos, basado en la
tecnología de reconstrucción 3D Multivisión.
Figura 76. Diagrama de control remoto Phantom 4 pro (Manual Phantom 4 pro, 2017).
Por último, el Photoscan parametriza la superficie del modelo asignándole a cada parte
de la malla un trozo de la imagen, creando así un mapa de textura (Ver Figura 77).
100
Figura 77. Modelo 3d de nube densa de la zona de estudio a partir de fotometría (Agisoft Photoscan).
101
CAPITULO V
En este capítulo se describe los resultados obtenidos a partir de las exploraciones geofísicas
aplicados a la geotecnia propiamente dicha debido a la investigación somera, los métodos aplicados
en este estudio nos ayudaran a discretizar mejor los distintos estratos en el subsuelo y la existencia
de un posible nivel freático.
5.
102
5.1 Levantamiento fotométrico de drone
Como apoyo al presente estudio nos fue muy útil contar con este dispositivo de vuelo, ya que nos dio
resultados interesantes pudiendo tener un amplio panorama del lugar de estudio a base de las
fotografías aéreas, por otro lado también se pudo procesar los datos con el objetivo de tener un las
curvas de nivel topográfico, lo cual nos sería útil para la presentación de un modelo seudo 3d de
acuerdo a las secciones de tomografía geoeléctricos, los resultados del levantamiento fotométrico se
presentan a continuación:
Figura 78. A) reconstrucción de imágenes utilizadas para modelo inicial, B) modelo de elevación digital (DEM)
(color azul: cota menor, color verde a amarillo: cota intermedia, color rojo: cota mayor), C) modelo orto
mosaico (Agisoft Photoscan).
Como resultado final, estos modelos nos sirvieron para exportar la data al programa civil 3D, donde
se realizaron los planos respectivos, a continuación, se muestra el plano de ubicación de las
exploraciones geofísicas realizadas en la localidad de Patashuaylla, Provincia de Urubamba,
departamento del Cusco.
103
N
CANAL
N
INGRESO
CANAL
CANAL
LINEA DE
MEDIA TENSION AREA DE
ESTUDIO
MASW-7 UBICACION
MASW-6
MASW-5
LEYENDA
MASW-1
ERT-01
LS-01
MASW-8
MASW-3
MASW-4
COORDENADAS MASW
ESTE NORTE
LS-02
LINEAS
MEDIO/MASW
LS-03
LS-04
ERT-01
ERT-02
LS-01
MASW-1
UNIVERSIDAD NACIONAL DE SAN AGUSTÍN
FACULTAD DE GEOLOGÍA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERÍA GEOFÍSICA
5.2 Perfiles geosísmicos de refracción sísmica
El área donde se realizaron las líneas sísmicas se constituye generalmente de depósitos cuaternarios
(depósitos aluviales), según la geología del lugar, a causa de eso se pudieron obtener los resultados
parecidos y congruentes para todas las líneas.
Esta línea sísmica tiene una longitud de 46 metros se realizó con intervalos de 2m para ganar
resolución ya que el estudio tiene como objetivo identificar los materiales yacentes en el lugar
donde serán el soporte de estructuras importantes.
Tenemos la presencia de dos estratos según las dromocronas en general, el primer estrato
correspondería a un depósito aluvial constituido de material orgánico por ser terreno de cultivos
con presencia de arcillas y gravas las cuales presentan un grado de compactación mala a
media, con velocidades de compresión que varían desde 327 – 588m/s, por otra parte, este
estrato estaría abarcando desde una profundidad de 3.0 a 4.5 m aproximadamente.
Para el segundo estrato se estima la presencia del mismo material aluvial, pero con
características distintas a la primera en cuanto al grado de densidad, sin embargo, se induciría
la presencia de gravas con presencia de arenas con un grado de densidad media a densa, sus
velocidades de compresión que varían desde 588 – 712m/s. y abarca desde una profundidad
promedio desde los 3.5m hasta los 16m aproximadamente según la refracción sísmica.
Esta línea está conformada por dos líneas sísmicas, su adquisición se hizo con traslape de
líneas, abarcando una longitud de 94m. Ya que la idea era abarcar mayor área de investigación
manteniendo la misma resolución del estudio.
El primer estrato estaría conformado por un material aluvial con presencia de material de
cultivo con velocidades de compresión que varían desde 380 a 614m/s, este estrato tendría
una potencia en cuanto a su espesor que varía desde los 3.5 hasta los 5.5m aproximadamente.
Seguido se estima la presencia de gravas con presencia de arenas posiblemente, este material
seria medianamente denso a denso con velocidades de ondas primarias que varían desde 614
104
a 720m/s. donde su espesor estaría variando desde los 3.5m hasta los 16 metros
posiblemente.
Está constituida de dos líneas de refracción sísmica haciendo una longitud total de 94 metros,
la ubicación de esta línea se encuentra junto a la línea de tomografía geoeléctrica 01.
El primer estrato estaría conformado por un material aluvial poco denso, con velocidades de
compresión que varían desde 337 a 592m/s, este estrato tendría una que varía desde los 2.5
hasta los 5.0m aproximadamente.
Está constituida de dos líneas de refracción sísmica haciendo una longitud total de 94 metros,
la ubicación de esta línea se encuentra junto a la línea de tomografía geoeléctrica 02.
El primer estrato estaría conformado por un material aluvial poco denso, con velocidades de
compresión que varían desde 368 a 629m/s, este estrato tendría un espesor que varía desde
los 2.5 hasta los 5.0m aproximadamente.
105
El perfil correspondiente a la línea de RS04 se presenta en el ANEXO A1 - PERFIL 04,
considerando la dromocrona del mismo, así como también el modelo inicial y el modelo final
con su inversión.
Está línea tiene una longitud de 46 metros, con fin de corroborar la data y tener más
información en el área de trabajo su ubicación fue estratégicamente de forma diagonal entre la
línea sísmica 02 y 03.
El primer estrato estaría conformado por un material aluvial poco denso con velocidades de
compresión que varían desde 385 a 617m/s, este estrato tendría una potencia que varía desde
los 2.0 hasta los 4.0m aproximadamente.
Seguido se estima la presencia de material aluvial, este material seria medianamente denso a
denso con velocidades de ondas primarias que varían desde 617 a 715m/s. donde su espesor
estaría variando desde los 2.0m hasta los 16 metros posiblemente.
La interpretación global de las líneas de refracción sísmica se dio a partir de los siguientes intervalos
en cuanto a sus velocidades, tomando velocidades patrón analizados con los resultados entre sí,
dándonos resultados congruentes, Ver Anexo A1.
106
5.3 Perfiles unidimensionales de MASW
En el Anexo A.2 se presenta el perfil unidimensional de todos los ensayos MASW donde se
muestran: la curva de dispersión, el perfil sísmico unidimensional y la tabla de velocidades de corte
promedio para los primeros 30 m de investigación (Vs30).
Por otro lado, se hizo la utilización del Cuadro 3. (Clasificación sísmica del suelo (código internacional
IBC-2012) y Cuadro 4. (Clasificación del suelo Norma técnica E-030,2016), Para clasificar
sísmicamente el suelo definiendo el tipo y nombre de suelo de acuerdo a la velocidad promedio de los
30m más superficiales (Vs 30), tal como se muestra en el siguiente cuadro.
Los estudios de MASW, se obtuvieron con el armado de líneas de 57.5m de longitud, pudiendo
llegar a los 30m de profundidad en promedio.
Se identificaron dos estratos principales dentro de los cuales su velocidad de corte varia de
forma moderada, a continuación, se analiza el perfil unidimensional.
Para este perfil unidimensional tenemos desde superficie la presencia de un material aluvial
poco rígido, su velocidad de corte varía desde 206 a 222m/s hasta una profundidad
aproximada de 4.8m, subyacente a este material se estima la presencia de material aluvial
medianamente rígido a rígido con velocidades de corte que abarcan desde los 236 a 280m/s.
107
Este estrato llega hasta una profundidad aproximada de 11.0 a 30m, a su vez se observa que
la velocidad aumenta a profundidad.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma notable, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 221 a 242m/s hasta una profundidad aproximada de 4.8m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 322 a 355m/s. Este estrato llega hasta una profundidad aproximada de 12m, seguido
tenemos un material medianamente rígido con velocidades de corte que varían desde 239 a
345m/s con una profundidad que abarca desde los 12 a 30m aproximadamente.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma notable, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 242 a 236m/s hasta una profundidad aproximada de 3.7m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 291 a 316m/s. Este estrato llega hasta una profundidad aproximada de 10m, seguido
tenemos un material medianamente rígido a rígido con velocidades de corte que varían desde
284 a 354m/s con una profundidad que abarca desde los 10 a 30m aproximadamente.
108
El perfil unidimensional correspondiente al MASW03 se presenta en el ANEXO A2 - PERFIL
03, donde se considera la curva de dispersión, así como su clasificación de suelos según la
norma IBC (2012) y la norma técnica peruana E-030.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma notable, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 205 a 253m/s hasta una profundidad aproximada de 3.7m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 303 a 319m/s. Este estrato llega hasta una profundidad aproximada de 10m, seguido
tenemos un material medianamente rígido con velocidades de corte que varían desde 277 a
312m/s con una profundidad que abarca desde los 10 a 30m aproximadamente.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma notable, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 237 a 262m/s hasta una profundidad aproximada de 2.6m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 308 a 318m/s. Este estrato llega hasta una profundidad aproximada de 8.8m,
seguido tenemos un material medianamente rígido con velocidades de corte que varían desde
272 a 321m/s con una profundidad que abarca desde los 8.8 a 30m aproximadamente.
109
5.3.6 Perfil unidimensional MASW – 06
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma moderada, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 224 a 252m/s hasta una profundidad aproximada de 2.6m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 298 a 305m/s. Este estrato llega hasta una profundidad aproximada de 7m, seguido
tenemos un material medianamente rígido con velocidades de corte que varían desde 293 a
337m/s con una profundidad que abarca desde los 7 a 30m aproximadamente.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma moderada, a continuación, se analiza el perfil unidimensional.
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 193 a 257m/s hasta una profundidad aproximada de 2.6m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 294 a 294m/s. Este estrato llega hasta una profundidad aproximada de 7.3m,
seguido tenemos un material medianamente rígido con velocidades de corte que varían desde
287 a 305m/s con una profundidad que abarca desde los 7.3 a 30m aproximadamente.
Se identificaron tres estratos principales dentro de los cuales su velocidad de corte varia de
forma moderada, a continuación, se analiza el perfil unidimensional.
110
Se tiene desde superficie la presencia de un material aluvial poco rígido, su velocidad de corte
varía desde 212 a 287m/s hasta una profundidad aproximada de 4.8m, subyacente a este
material se estima la presencia de material aluvial rígido con velocidades de corte que abarcan
desde los 298 a 316m/s. Este estrato llega hasta una profundidad aproximada de 13.5m,
seguido tenemos un material medianamente rígido con velocidades de corte que varían desde
297 a 328m/s con una profundidad que abarca desde los 13.5 a 30m aproximadamente.
En base a los resultados de las velocidades se ondas primarias Vp y ondas de corte Vs, obtenidos de
los estudios geofísicos, se han obtenido los siguientes valores del módulo de elasticidad del suelo
para cada estrato:
Cuadro 16. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 01 & MASW 01.
111
Cuadro 17. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 02 & MASW 02.
Relación
Módulo de Módulo de Módulo
Profundidad Densidad de
Vp (m/s) Vs (m/s) Corte Gd Young Ed Volumétrico
(m) (Tn/m3) poisson
(Kg/cm2) (Kg/cm2) Kd (Kg/cm2)
µ
1.00 401 221 1.8 0.23 897.1 2208.6 1368.4
1.67 451 238 1.8 0.26 1036.0 2607.8 1800.0
2.63 460 240 1.8 0.27 1053.6 2666.0 1892.7
3.68 463 208 1.8 0.33 794.6 2115.8 2090.3
4.83 463 242 1.8 0.26 1075.7 2716.7 1908.8
6.05 463 323 1.8 0.02 1910.3 3896.5 1352.6
7.37 641 356 1.8 0.23 2321.3 5696.6 3478.2
8.77 711 337 1.8 0.31 2079.8 5453.6 4811.7
10.26 728 285 1.8 0.38 1486.7 4090.9 5493.5
11.84 730 249 1.8 0.41 1138.8 3205.4 5766.8
13.51 730 239 1.8 0.42 1049.2 2969.3 5826.5
15.26 730 260 1.8 0.40 1236.9 3460.4 5701.4
Cuadro 18. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16m, realizados en las líneas
LS 03 & MASW 03.
112
Cuadro 19. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 04 &MASW 04.
Cuadro 20. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 04 – MASW 05.
113
Cuadro 21. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 03 – MASW 06.
Cuadro 22. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 16 m, realizados en las líneas
LS 02 – MASW 07.
114
Cuadro 23. Parámetros pseudo - dinámicos, recomendados desde una profundidad de 1 a 24 m, realizados en las líneas
LS 05 – MASW 08.
Se realizó el método de tomografía geoeléctrica analizando los objetivos, dentro de ellos se destaca
principalmente la posible existencia de filtraciones de agua en el subsuelo o la existencia de un nivel
freático, en consecuencia, se ha distribuido estratégicamente dos líneas de tomografía geoeléctrica
por el dispositivo de arreglo Dipolo - dipolo cubriendo el área de estudio.
Las anomalías de color azul-celeste representan valores de resistividad bajas que oscilan entre 85 –
170 ohm*m, que se manifiestan en profundidad, asociados a la presencia de materiales con contenido
de humedad.
Las anomalías de color naranja-rojo representan valores altos de resistividad que oscilan entre 235 –
340 ohm*m., que podrían estar asociados a suelos resistivos.
115
5.5.1 Perfil Geoeléctrico TG – 01
Se pudieron identificar principalmente tres estratos según las resistividades registradas, esta
línea tiene una longitud de 90m como se muestra en la pseudosección, estos estratos
identificados se extienden a lo largo del perfil denotando con claridad las anomalías principales.
Se pudieron identificar principalmente tres estratos según las resistividades registradas, esta
línea tiene una longitud de 90m como se muestra en la pseudosección, estos estratos
identificados se extienden a lo largo del perfil denotando con claridad las anomalías principales.
116
El primer estrato estaría conformado de un material medianamente resistivo posiblemente
estaría asociado a un material aluvial con presencia de arcillas y orgánicos, con resistividades
que varían, desde 173 – 235 ohm*m, con profundidades que varían desde los 2 a 5m de
profundidad, por otra parte.
Los resultados de los distintos estudios realizados para este proyecto nos va permitir sacar
conclusiones particulares para cada método, sin embargo, la correlación entre los métodos aplicados
nos va a permitir una interpretación más relevante, ya que cada método nos brinda características
particulares para cada material, entonces analizando los resultados se llega a conclusiones
interesantes.
117
teniendo resultados congruentes. Ver Lámina 01, Pseudo – 3D Tomografía Geoeléctrica y Lamina
02, Pseudo – 3D de Refracción Sísmica.
Estos depósitos son utilizados como canteras en algunos lugares cerca del lugar de estudio, entonces
tomando como referencia esta aclaración, nos permite afirmar que el material que se obtiene a través
del estudio de refracción sísmica y MASW, serian gravas con contenido de arcillas y limos en la zona
de estudio, ahora el grado de densidad y rigidez lo determina los materiales adjuntos a estas, en todo
caso podría repercutir la humedad a profundidad conservada por los materiales fino como arcillas y
limos, la cual hace que el material sea más blando.
Se saca esta conclusión básicamente a partir del estudio de MASW y Tomografía geoeléctrica ya que
la primera capa se comporta como un material blando a medianamente rígido, seguido se estima la
presencia de un material rígido, lo cual se muestra claramente en el perfil unidimensional de los
MASW y a su vez se denota claramente como un material resistivo en los perfiles de Tomografía
geoeléctrica, por otra parte los perfiles de Refracción sísmica nos ayudan a determinar que el material
que se obtiene hasta la segunda capa u horizonte no tiene un cambio brusco de velocidades de
ondas P como para estimar la presencia de un material totalmente diferente al material aluvial.
Entonces hecha esta aclaración se tendría material aluvial hasta profundidades más allá de las
investigaciones probablemente.
118
TG01 TG02
120
RS02 RS03
RS04
RS05
RS01
121
Sin embargo, los perfiles de Refracción sísmica no nos permiten estimar este tercer estrato debido a
sus principios fundamentales lo cual nos indica que V1>V2>V3, (Ver Figura 82). Tomando esa
definición podemos afirmar que en este caso no se cumple, por ende, nos muestra solo los rayos
refractores en dos capas y la tercera capa lo incluye en la segunda capa.
Por otra parte, los estudios de Tomografía geoeléctrica denotan este tercer estrato con resistividades
medianamente bajas, es por eso que nos induce a la estimación del material que lo conforma, debido
a esta característica nos atrevemos a afirmar que se podría asociar a materiales como gravas con
presencia de arcillas y limos ya que este material conserva la humedad y a su vez se comporta como
un material medianamente blando cuando contiene humedad (Ver Figura 83).
119
N
AREA DE
ESTUDIO
MASW-06 MASW-05
2860
2860
MASW-03 MASW-04
H-01 H-01
ALTURA (m.s.n.m)
ALTURA (m.s.n.m)
2850
2850
H-02
H-02
H-03 H-03
2840
2840
2830
2830
0 10 20 30 40 50 60 70 80 90 -10 0 10 20 30 40 50 60 70 80 90
DISTANCIA (m) DISTANCIA (m)
UBICACION
LEYENDA
SECCIÓN GEOSÍSMICA 01 SECCIÓN GEOSÍSMICA 02
MASW-01
2865 MASW-07
ALTURA (m.s.n.m)
2855
H-01 MASW-02
2860
2850
ALTURA (m.s.n.m)
H-01
H-02 2855
2845
H-03 2850 H-02
2840 MASW
2845 H-03
2835
0 5 10 15 20 25 30 35 40 45 2840
DISTANCIA (m)
2835
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95
DISTANCIA (m)
ESTE NORTE ESTE NORTE
LINEAS
INICIO FINAL
H-01 MASW
2855 H-01 2855 MEDIO/MASW
1. Según las dromocronas obtenidas de los ensayos de refracción sísmica se pudieron obtener
principalmente dos capas representativas, El primer horizonte con velocidades de compresión que
abarcan de 327 a 417 m/s con espesores que no superan los 6.0 m., se trataría de un material
poco denso. Subyaciendo se estima la presencia de un material medianamente densas a densas
con una velocidad de ondas de compresión que llega hasta 764m/s.
2. Según los MASW Tenemos desde superficie la presencia de un material blando a medianamente
blando, su velocidad de ondas de corte varía de 193 a 287m/s hasta una profundidad aproximada
de 4.8m. como máximo, subyaciendo se estima la presencia de un material rígido, con velocidades
de ondas corte que varían desde 236 a 355 m/s. su profundidad varía desde 2.6 – 13.5m.
subyaciendo se estima la presencia de un material medianamente rígido con velocidades de
ondas corte que varían desde 229 a 235 m/s. su profundidad varía desde 7 – 30m en promedio
3. Según los estudios de tomografía geoeléctrica se distinguen tres horizontes. El primer Horizonte
estaría asociada a gravas con presencia de materiales finos, además de eso posiblemente con
presencia de humedad mínima, con resistividades que varían entre 170 – 235 ohm*m
aproximadamente. El segundo material estaría asociado a la presencia de material aluvial denso
permeable. Su resistividad varía desde 235 – 340 ohm*m. y por último el tercer horizonte estaría
representado por gravas con presencia de arcillas y limos posiblemente con presencia de
humedad, cuyas resistividades varían entre 85 - 170 ohm*m aproximadamente.
4. Los materiales que conforman el subsuelo en el área donde se realizó la prospección geofísica se
clasificó de acuerdo a lo propuesto por la International Building Code (IBC, 2009), los cuales se
clasifican como suelo tipo D (suelo rígido) con velocidad promedio de ondas de corte en los
primeros 30 m de profundidad de (Vs30), a su vez también se clasifica de acuerdo a la norma
Peruana E- 030, indicando que pertenece a un suelo tipo S2 (Suelo intermedio).
5. Se puede afirmar que el material que será soporte de estas estructuras importantes es
regularmente bueno a partir del segundo estrato lo cual se presenta desde una profundidad que
varía entre 2.6 a 6m aproximadamente.
6. Los parámetros pseudo - dinámicos de los estratos característicos del subsuelo se calcularon
sobre la base de los resultados obtenidos de los estudios geofísicos y valores de densidad
referenciales, los valores de estos parámetros se presentan en Anexo A4.
120
RECOMENDACIONES
1. Los métodos Geofísicos son estudios indirectos, que tienen un grado de aproximación aceptable,
sin embargo, se requiere confirmar los resultados estimados, con cualquier método directo.
2. Se recomienda la aplicación del método de tomografía geoeléctrica en la parte inferior del área de
trabajo paralela a la línea de refracción 01 y transversales a las líneas de refracción 02,03,04.
121
BIBLIOGRAFÍA
• Ásta Ólafsdóttir Elín, 2014, “Multichannel Analysis of Surface Waves: Methods for Dispersion
Analysis of Surface wave data”, University of Iceland.
• Barker, R.D., 1981, “The offset system of electrical resistivity sounding and its use with a
multicore cable. Geophysical Prospecting”, 29, 128–43.
• Bruce B. Redpath, 1973, “Seismic Refraction Exploration for Engineering Site Investigations”,
California.
• Carlos Huamán Egoávil, Roberth Aguilar Chuquimuni, 2000, “Implementación del Equipo de
Adquisición de Datos de Prospección Sísmica, Programa Científico PC – CISMID”.
• Chávez R., Tejero A., Delgado C., Cifuentes G. y Hernández E., 2011, Poster, L and Corner
arrays for 3D electric resistivity tomography: An alternative for geophysical surveys in urban
zones. American Geophysical Union.
• Choon B. Park, Mario Carnevale, 2010, “Optimum MASW Survey - Revisit after a Decade of
Use”.
• Choon B. Park, Miller Richard, Xia Jianghai, Nils Ryden, and Peter Ulriksen, 2003, "The
MASW Method - What and Where it is." EAGE 65th Conference and Exhibition. Stavanger.
122
• Choon B. Park, Miller, R., Xia Jianghai, J.,1999, “Multichannel analysis of surface waves”.
Geophysics, 64(3), 800 - 808.
• Coral Crespo Erick, 2015, “Generación de Modelos de Elevación Digital a Partir de Imágenes
Obtenidas Desde Vehículos Aéreos no Tripulados (UAVs)”, Pontifica Universidad Javeriana,
Bogotá -Colombia.
• Díaz Anaya Erika, 2010, “Implementación del Código Zondres2D para la Modelación Directa e
Inversa de Datos de Tomografía de Resistividad Eléctrica 2D”, Tesis de grado, Universidad
Nacional Autónoma de México, México.
• Dr. Ing. Alva hurtado Jorge, 2013, “Aplicación de la Geofísica en la Ingeniería Geotécnica”, XXI
CONEIC, Huancayo.
• Enrique Guadalupe G., Manuel Zea A., Iris Villafuerte S., Dante Flores O., “Estudio Geológico
– Geotécnico para el Relleno Sanitario de Machu Picchu y Pueblos Aledaños”.
• Estrada Luis A., 2012, “Prospección Geofísica”, Catedra de Geofísica, Universidad Nacional de
Tucumán.
• García Carrasco Pedro, 2013, “Avance en Técnicas Geofísicas para la Caracterización del
Subsuelo Mediante Innovación y el Uso de Herramientas de Gestión de Información Espacial”,
Tesis Doctoral, Universidad de Salamanca.
• Geometrics. Inc.,2009, ES-3000, GeodeTM y StrataVisorTM NZ/NZC, Manual del Operador P/N
28519-01 Rev K.
• Graff, K.,1991, “Wave motion in elastic solids”. Dover Publications Inc.: New York.
123
• Hamid N. Alsadi, 2017, “Seismic Hydrocarbon Exploration”, 2D and 3D Techniques, Seismic
waves.
• Hassan Asem Ahmed, 2014, “Electrical Resistivity Method for Water Content Characterisation of
Unsaturated Clay Soil”, Durham University.
• Hayashi, K, 2003, “Data Acquisition and Analysis of Active and Passive Surface Wave Methods”.
Short Course - SAGEEP 2003.
• Hayashi, K., 2006, “Surface Wave Data Acquisition and Field Methods”, Short Course
Geometrics Inc.
• Hayashi, K., 2008, “Development of the Surface-wave Methods and Its Application to Site
Investigations”, Thesis or Dissertatio, Kyoto University.
• INDECI – PNUD PER / 02/ 051, 2005, “Estudio de Mapa de Peligros de la Ciudad de Urubamba”.
• INDECI – PNUD PER / 02/ 051, 2005, “Estudio de Mapa de Peligros de la Ciudad de
Ollantaytambo”.
• International Building Code (IBC) por International Code Council., Inc., 2012.
• IRIS instruments, 2003, SYSCAL Pro Standard & Switch (48 - 72 or 96) Version, User´s Manual.
• Javier, M., G., 2008, “Levantamiento Sísmico de Refracción Somera y Levantamiento Geológico
en el Área del Casino, al Norte del Poblado de Urumaco, Proyecto de Grado para optar el Título
de Ingeniero Geofísico”, Universidad Simón Bolívar, Venezuela.
• José Manuel López Avalos, 2008, “Aplicación del método de refracción sísmica para la
determinación de velocidades de ondas P”, Tesis de Grado.
• Joshua Lee Hubbard, 2009, “Use of Electrical Resistivity and Multichannel Analysis of Surface
Wave Geophysical Tomography in Geotechnical Site Characterization of Dam”, The University of
Texas at Arlington.
124
• Keary, P., Brooks, M. and Hill, I., (2002), “An Introduction to Geophysical Exploration (3rd ed.)”.
Blackwell Science Ltd: Malden, MA.
• Loke M.H., Baker R.D., 1996. “Rapid least-squares inversion of apparent resistivity
pseudosections by quasi- Newton method”, Geophysical Prospecting.
• Loke, M. H., Chambers, J. E., Rucker, D. F., Kuras, O. and Wilkinson, P. B., 2013, Recent
developments in the direct-current geoelectrical imaging method, Journal of Applied Geophysics.
• Loke, M.H., 1994, “The inversion of two-dimensional resistivity data”. Unpubl. PhD thesis, Un. Of
Birmingham.
• Lopes, I., Santos, J.A., and Almeida, I.M.d. 2008, “O método das ondas superficiais: aquisição,
processamento e inversão”, in Revista Geotecnia. p. 79-109.
• Lopez H., Loke H., Fanton O., Rubí E., 2004, “Técnicas prácticas para investigación de
resistividad en dos y tres dimensiones (tomografía eléctrica 2D y 3D)”.
• Manilla Aceves Alfonso, 2003, “Geofísica Aplicada en los Proyectos básicos de Ingeniería
Civil”, Publicación Técnica No. 229, México.
• Manilla Aceves, A. A., 2003,” Geofísica Aplicada en los Proyectos básicos de Ingeniería Civil”,
Publicación Técnica No. 229 Sanfandila, Qro, 2003.
125
• Moller, I., Christensen, N.B. and B. H. Jacobsen, 1995, Fast approximate 2D interpretation of
resistivity profile data, in proceedings of the 1st Meeting Environmental and Engineering
Geophysics, Torino 468 – 471.
• Orellana E., 1982, “Prospección Geoeléctrica en corriente continua”. Paraninfo, Madrid, España.
• Parasnis, D.S., 1996, “Principles of Applied Geophysics (Fifth Edition)”, Chapman & Hall,
London, 456 pp.
• PhD. Francisco Olivera, 1996, “Hydrology of the Urubamba River”, Center for Research in
Water Resources, University of Texas at Austin.
• Phil sirles, 2016, “Overview of Geophysical Investigation Methods (for use on Earthen
Embankments)”, course code: E0274, National Dam Safety Program Technical Seminar.
• Rodríguez Manuel Arlande, 2001, “Geofísica Aplicada a la Obra Civil. Método Geoeléctrico y
Sísmica de Refracción, Casos Prácticos”. Geoconsult Ingenieros Consultores, Valencia, España.
• Saavedra Morales Oscar, 2013, “Análisis de Ondas Sísmicas en Entorno MATLAB”, Tesis de
grado, Instituto Politécnico Nacional, México
• Suarez Gómez Luis, 2017, “Optimización de Parámetros para la Adquisición de Datos Medidos
in Situ por Medio del Método de Análisis Multicanal de Ondas Superficiales (MASW)”, Tesis de
Grado, Universidad Técnica Federico Santa María.
• Telford, W.M., Geldart, L.P., Sheriff, R.E. and Keys, D.A., 1990, “Applied Geophysics (Second
Edition)”, Cambridge University Press, Cambridge, 770 pp.
• Van Der Hilst, Robert, 2004, "Chapter 4: Seismology." MIT Open Course Ware.
http://ocw.mit.edu/NR/rdonlyres/Earth--Atmospheric--and-Planetary-Sciences/12-201Fall-
2004/4813C5AE-9A5F-4440-8A49-39559A713D67/0/ch4.pdf.
126
• Zonge, Ken, Jeff Wynn, and Scott Urquatt., 2005, "Resistivity, Induced Polarization, and
Complex Resistivity." In Near-Surface Geophysics, by Dwain K. Butler, 265-300. Tulsa: Society of
Engineering Geophysicist.
WEBGRAFÍA
• https://www.bing.com/maps.
• http://www.geometrics.com/applications/frequency-asked-questions/seismic-refraction/Original
Articles.
• http://www.iris-instruments.com/syscal-pro.html.
• http://riosdelplaneta.com.
• http://www.geofisicaexploraciones.com.
• https://es.wikipedia.org/wiki/Onda#/media/File:Wave.png.
• https://es.slideshare.net/JDaviSP/9-ensayo-de-sismica-de-refracción.
127
ANEXO A
EXPLORACIONES DE CAMPO
ANEXO A1
REFRACCIÓN SÍSMICA
ANEXO A2
MASW