en la ciudadela inca de machu picchu: cusco-perú - Ingemmet
en la ciudadela inca de machu picchu: cusco-perú - Ingemmet
en la ciudadela inca de machu picchu: cusco-perú - Ingemmet
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
XIII Congreso Peruano <strong>de</strong> Geología. Resúm<strong>en</strong>es Ext<strong>en</strong>didos<br />
Sociedad Geológica <strong>de</strong>l Perú<br />
APLICACIÓN DE LA GEOFÍSICA (RESISTIVIDAD 2D Y GEORADAR-GPR)<br />
EN LA CIUDADELA INCA DE MACHU PICCHU: CUSCO-PERÚ<br />
Walter Pari 1 , Víctor Carlotto 1,2 , José Cárd<strong>en</strong>as 2 , Martín Oviedo 2 , Mart<strong>en</strong> Douma 3<br />
Mel Best 3 & Peter Bobrowsky 3<br />
1 INGEMMET, Av. Canadá 1470, San Borja, Lima-Perú wpari@ingemmet.gob.pe<br />
2 UNSAAC Av. <strong>de</strong> <strong>la</strong> Cultura, s/n Cusco-Perú<br />
3 Servicio Geológico <strong>de</strong> Canadá<br />
INTRODUCCIÓN<br />
En vista <strong>de</strong> <strong>la</strong> hipótesis p<strong>la</strong>nteada <strong>de</strong> <strong>la</strong> exist<strong>en</strong>cia <strong>de</strong> <strong>de</strong>slizami<strong>en</strong>tos profundos <strong>en</strong> <strong>la</strong> ciuda<strong>de</strong><strong>la</strong> Inca <strong>de</strong><br />
Machu Picchu era necesario realizar investigaciones geofísicas sistemáticas <strong>en</strong> el sitio con el objeto <strong>de</strong><br />
obt<strong>en</strong>er información geológica y geotécnica <strong>de</strong>l subsuelo, características <strong>de</strong> los sedim<strong>en</strong>tos<br />
superficiales, <strong>la</strong> topografía <strong>de</strong>l substrato rocoso, así como <strong>la</strong> pres<strong>en</strong>cia <strong>de</strong> fal<strong>la</strong>s y p<strong>la</strong>nos <strong>de</strong><br />
<strong>de</strong>slizami<strong>en</strong>to. En el marco <strong>de</strong>l Proyecto Multinacional Andino: Geoci<strong>en</strong>cias para <strong>la</strong>s Comunida<strong>de</strong>s<br />
Andinas (MAP:GAC) con <strong>la</strong> asist<strong>en</strong>cia <strong>de</strong>l Servicio Geológico <strong>de</strong> Canadá, el INGEMMET y <strong>la</strong><br />
Universidad <strong>de</strong>l Cusco, pudieron realizar <strong>la</strong>s investigaciones geofísicas cuyos resultados contribuy<strong>en</strong> a<br />
<strong>la</strong> evaluación <strong>de</strong> los peligros, como un aporte para <strong>la</strong> conservación <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong> Inca.<br />
Los trabajos <strong>de</strong> campo han consistido <strong>en</strong> prospección geofísica aplicando, <strong>en</strong> su primera fase los<br />
métodos <strong>de</strong> Resistividad Eléctrica y Electromagnétismo EM-31/EM-34 (Mayo <strong>de</strong>l 2004), y <strong>en</strong> una<br />
segunda fase aplicando los Métodos <strong>de</strong> Georadar <strong>de</strong> p<strong>en</strong>etración-GPR y Resistividad 2D (Junio-Julio,<br />
2005). El método electromagnético no es tratado aquí por no haber t<strong>en</strong>ido resultados satisfactorios <strong>en</strong><br />
Machu Picchu.<br />
UBICACIÓN - PRESENTACIÓN<br />
La Ciuda<strong>de</strong><strong>la</strong> Inca <strong>de</strong> Machu Picchu se <strong>en</strong>cu<strong>en</strong>tra ubicada a 112.5 Km al noroeste <strong>de</strong> <strong>la</strong> ciudad <strong>de</strong>l<br />
Cusco, a una altura <strong>de</strong> 2,350 msnm (Fig. 1).<br />
Se localiza <strong>en</strong> pl<strong>en</strong>a Cordillera Ori<strong>en</strong>tal, don<strong>de</strong> el río<br />
Urubamba forma el cañón <strong>de</strong>l Urubamba. Descubierta a<br />
inicios <strong>de</strong>l siglo pasado, el sitio recibe<br />
aproximadam<strong>en</strong>te un millón <strong>de</strong> turistas por año. En<br />
Machu Picchu y sus alre<strong>de</strong>dores se pres<strong>en</strong>tan<br />
frecu<strong>en</strong>tem<strong>en</strong>te <strong>de</strong>slizami<strong>en</strong>tos superficiales, caída <strong>de</strong><br />
rocas y aluviones, produci<strong>en</strong>do daños materiales y<br />
pérdidas humanas como ocurrió <strong>en</strong> Aguas Cali<strong>en</strong>tes el<br />
2004.<br />
MÉTODO DE RESISTIVIDAD-WENNER 2D<br />
Este método se caracteriza por el estudio <strong>de</strong> <strong>la</strong>s<br />
variaciones <strong>de</strong> los parámetros físicos <strong>de</strong> <strong>la</strong>s rocas o <strong>de</strong><br />
los suelos y <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong> factores que afectan <strong>la</strong><br />
resistividad <strong>de</strong> los materiales, como <strong>la</strong> porosidad, agua<br />
<strong>en</strong> los poros, conductividad <strong>en</strong> los granos minerales,<br />
grado <strong>de</strong> compactación, etc. La resistividad y<br />
Fig. 1. Mapa <strong>de</strong> ubicación <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong> conductividad <strong>de</strong> los materiales terrestres se basa <strong>en</strong> <strong>la</strong><br />
Inca <strong>de</strong> Machu Picchu<br />
aplicación <strong>de</strong> <strong>la</strong> Ley <strong>de</strong> Ohm (V = IR) don<strong>de</strong> V = voltaje<br />
(Voltios), I = corri<strong>en</strong>te (Amperios) y R = resist<strong>en</strong>cia<br />
(Ohmios). El método dispone <strong>de</strong> difer<strong>en</strong>tes configuraciones y para el pres<strong>en</strong>te estudio se ha aplicado <strong>la</strong><br />
50
XIII Congreso Peruano <strong>de</strong> Geología. Resúm<strong>en</strong>es Ext<strong>en</strong>didos<br />
Sociedad Geológica <strong>de</strong>l Perú<br />
configuración W<strong>en</strong>ner con 48 electrodos con espaciami<strong>en</strong>tos <strong>de</strong> 5 m, que ofrece mejor resolución <strong>en</strong><br />
profundidad pero que su p<strong>en</strong>etración es somera (50 m). Colecta datos <strong>de</strong> resistividad <strong>en</strong> dos<br />
dim<strong>en</strong>siones (2D) usando un sistema <strong>de</strong> múltiples electrodos.<br />
Se realizaron 8 líneas <strong>en</strong> <strong>la</strong> carretera <strong>de</strong> acceso Hiram Bingham y una línea <strong>en</strong> <strong>la</strong> P<strong>la</strong>za Principal <strong>de</strong> <strong>la</strong><br />
Ciuda<strong>de</strong><strong>la</strong>. El equipo geofísico empleado fue el Sistema Syscal Iris con 48 electrodos <strong>de</strong> pot<strong>en</strong>cial. El<br />
procesami<strong>en</strong>to y los resultados obt<strong>en</strong>idos se han <strong>de</strong>sarrol<strong>la</strong>do <strong>en</strong> base al Programa <strong>de</strong> Mo<strong>de</strong><strong>la</strong>mi<strong>en</strong>to <strong>de</strong><br />
Resistividad <strong>en</strong> 2 D. (Inversión <strong>de</strong> resistividad, RES2DINV V.4).<br />
RESULTADOS<br />
El background da valores <strong>de</strong> resistividad mayores a los 3000 ohm-m que son típicos <strong>de</strong> rocas<br />
graníticas, y correspon<strong>de</strong> a los aflorami<strong>en</strong>tos <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong>. Los datos <strong>de</strong> resistividad <strong>de</strong> los ocho<br />
tramos superiores <strong>de</strong> <strong>la</strong> carretera (línea 0 a línea 7) nos indican, <strong>en</strong> color amarillo (Fig. 2), valores<br />
bajos por comparación al background, <strong>en</strong>tre 700 y 1200 ohm-m. Estas zonas <strong>de</strong> baja resistividad se<br />
interpretan como rocas int<strong>en</strong>sam<strong>en</strong>te fracturadas con pres<strong>en</strong>cia <strong>de</strong> flujos <strong>de</strong> agua y coincid<strong>en</strong> con una<br />
fal<strong>la</strong> que atraviesa <strong>la</strong> zona.<br />
Fig. 2. Perfiles <strong>de</strong> resistividad y t<strong>en</strong>d<strong>en</strong>cias <strong>en</strong> <strong>la</strong> carretera Hiram Bingham<br />
En <strong>la</strong> ciuda<strong>de</strong><strong>la</strong> <strong>inca</strong> se realizó un perfil transversal a <strong>la</strong> P<strong>la</strong>za Principal don<strong>de</strong> zonas puntuales <strong>de</strong> baja<br />
resistividad (color azul) se asocian a <strong>la</strong>s terrazas superiores al oeste y este, probablem<strong>en</strong>te <strong>de</strong>bido a<br />
suelos <strong>de</strong> mayor espesor para <strong>la</strong> agricultura. Así mismo se observa baja resistividad <strong>en</strong> el área bajo <strong>la</strong>s<br />
construcciones como una capa <strong>de</strong>lgada quizás <strong>de</strong> arcil<strong>la</strong>s o suelos que agregaron los <strong>inca</strong>s para hacer<br />
una bu<strong>en</strong>a superficie <strong>de</strong> construcción para los edificios. La zona <strong>de</strong> baja resistividad al <strong>la</strong>do este <strong>de</strong> <strong>la</strong><br />
P<strong>la</strong>za Principal pue<strong>de</strong> ser re<strong>la</strong>cionada a <strong>la</strong>s fracturas saturadas con agua. La alta resistividad mostrada<br />
<strong>en</strong> <strong>la</strong> base <strong>de</strong> <strong>la</strong> p<strong>la</strong>za está re<strong>la</strong>cionada con <strong>la</strong> parte superior <strong>de</strong>l substrato rocoso <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong>. En<br />
los extremos este y oeste <strong>de</strong>l perfil, se pued<strong>en</strong> dibujar p<strong>la</strong>nos <strong>de</strong> <strong>de</strong>slizami<strong>en</strong>tos superficiales que<br />
corroboran <strong>la</strong> cartografía geológica. Bajo <strong>la</strong> P<strong>la</strong>za Principal no se aprecia p<strong>la</strong>nos <strong>de</strong> fal<strong>la</strong> ni <strong>de</strong><br />
<strong>de</strong>slizami<strong>en</strong>tos.<br />
Fig. 3. Perfil geoeléctrico <strong>en</strong> <strong>la</strong> P<strong>la</strong>za Principal<br />
51
XIII Congreso Peruano <strong>de</strong> Geología. Resúm<strong>en</strong>es Ext<strong>en</strong>didos<br />
Sociedad Geológica <strong>de</strong>l Perú<br />
MÉTODO GEORADAR DE PENETRACIÓN-GPR (Ground P<strong>en</strong>etrating Radar)<br />
El levantami<strong>en</strong>to georadar fue realizado empleando el Sistema Pulse EKKO 100 fabricado <strong>en</strong> Canadá.<br />
Incluye una unidad <strong>de</strong> control, un procesador y dos ant<strong>en</strong>as (trasmisor y un receptor) <strong>de</strong> alta frecu<strong>en</strong>cia<br />
(50, 100 y 200 Mhz). El procedimi<strong>en</strong>to está principalm<strong>en</strong>te re<strong>la</strong>cionado al método <strong>de</strong> reflexión<br />
sísmica, el cual está constituido por un trasmisor (Tx) que emite una señal hacia el subsuelo. El<br />
regreso <strong>de</strong> <strong>la</strong> onda electromagnética es <strong>de</strong>tectado y registrado por el receptor (Rx).<br />
El GPR funciona g<strong>en</strong>eralm<strong>en</strong>te con una separación fija <strong>de</strong>l trasmisor-receptor. La adquisición continua<br />
<strong>de</strong> los datos se ejecuta movi<strong>en</strong>do <strong>la</strong>s ant<strong>en</strong>as <strong>de</strong>l GPR continuam<strong>en</strong>te a lo <strong>la</strong>rgo <strong>de</strong>l perfil. La<br />
información recolectada se pres<strong>en</strong>ta <strong>en</strong> tiempo real y <strong>de</strong> manera continua <strong>en</strong> <strong>la</strong> pantal<strong>la</strong> <strong>de</strong>l<br />
computador, pres<strong>en</strong>tándose como una imag<strong>en</strong> <strong>de</strong>l material subyac<strong>en</strong>te a lo <strong>la</strong>rgo <strong>de</strong> <strong>la</strong> línea <strong>de</strong><br />
investigación. Las reflexiones ocurr<strong>en</strong> don<strong>de</strong> hay un cambio <strong>en</strong> <strong>la</strong> constante dieléctrica y <strong>la</strong><br />
conductividad <strong>de</strong>l medio (Fig. 4a). Las secciones <strong>de</strong> los datos ploteados son semejantes a <strong>la</strong>s secciones<br />
sísmicas, con tiempo <strong>en</strong> nano-segundos <strong>en</strong> el eje vertical y <strong>la</strong> posición <strong>en</strong> el eje horizontal (Fig. 4b). La<br />
profundidad <strong>de</strong> investigación varía <strong>de</strong>s<strong>de</strong> m<strong>en</strong>os <strong>de</strong> un metro hasta aproximadam<strong>en</strong>te 30 a 40 m.<br />
Fig. 4. Perfil <strong>de</strong> Reflexión <strong>de</strong>l GPR<br />
Las líneas <strong>de</strong>l Georadar-GPR muestran reflectores electromagnéticos que correspond<strong>en</strong> a patrones <strong>de</strong><br />
frecu<strong>en</strong>cia bi<strong>en</strong> <strong>de</strong>finidos, los que indican c<strong>la</strong>ram<strong>en</strong>te cambios litológicos causados <strong>la</strong> permitividad<br />
re<strong>la</strong>tiva y conductividad eléctrica <strong>de</strong>l medio.<br />
RESULTADOS<br />
Uno <strong>de</strong> los b<strong>la</strong>ncos <strong>de</strong>l GPR fue el área <strong>de</strong> <strong>la</strong> P<strong>la</strong>za Principal que separa el Intihuatana y el sector<br />
urbano <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong>. El objetivo principal fue obt<strong>en</strong>er <strong>la</strong> configuración <strong>de</strong>l subsuelo y <strong>de</strong>l substrato<br />
rocoso bajo <strong>la</strong> p<strong>la</strong>za principal, consecu<strong>en</strong>tem<strong>en</strong>te calcu<strong>la</strong>r el espesor <strong>de</strong> rell<strong>en</strong>o y <strong>la</strong> topografía <strong>de</strong>l<br />
substrato. Otro <strong>de</strong> los objetivos fue <strong>de</strong>terminar <strong>la</strong> pres<strong>en</strong>cia o <strong>la</strong> aus<strong>en</strong>cia <strong>de</strong> p<strong>la</strong>nos <strong>de</strong> fal<strong>la</strong> o<br />
<strong>de</strong>slizami<strong>en</strong>to <strong>de</strong> dirección norte-sur.<br />
Tres facies o unida<strong>de</strong>s principales fueron id<strong>en</strong>tificadas <strong>en</strong> <strong>la</strong>s secciones georadar. La primera facies se<br />
sitúa al techo <strong>de</strong> <strong>la</strong>s secciones, <strong>en</strong>tre 0.5 y 1 m. (Fig. 5). Muestra una fuerte reflexión que pue<strong>de</strong><br />
correspon<strong>de</strong>r al suelo agríco<strong>la</strong> ya que los <strong>inca</strong>s normalm<strong>en</strong>te t<strong>en</strong>ían suelos <strong>de</strong> 0.5 m. (Wright &<br />
Val<strong>en</strong>cia, 1999).<br />
La segunda unidad parece ext<strong>en</strong><strong>de</strong>rse hasta aproximadam<strong>en</strong>te 4 m <strong>de</strong> profundidad y repres<strong>en</strong>ta una<br />
unidad <strong>de</strong> subsuelo compuesto principalm<strong>en</strong>te por bloques, fragm<strong>en</strong>tos <strong>de</strong> rocas, y otros productos<br />
posiblem<strong>en</strong>te acarreados <strong>de</strong> <strong>la</strong>s canteras por los <strong>inca</strong>s (rell<strong>en</strong>o). La falta <strong>de</strong> coher<strong>en</strong>cia, <strong>en</strong> muchas <strong>de</strong><br />
<strong>la</strong>s reflexiones, repres<strong>en</strong>ta <strong>la</strong> dispersión <strong>de</strong> <strong>la</strong>s ondas <strong>de</strong> radar ocasionado por los escombros. Los<br />
perfiles muestran que <strong>la</strong> p<strong>la</strong>za principal es una <strong>de</strong>presión que habría sido rell<strong>en</strong>ada durante <strong>la</strong><br />
construcción <strong>de</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong>.<br />
Infrayaci<strong>en</strong>do <strong>la</strong>s facies <strong>de</strong> suelos e indicando los límites <strong>de</strong> capacidad <strong>de</strong> p<strong>en</strong>etración <strong>de</strong>l radar, se<br />
ti<strong>en</strong>e una tercera unidad pobrem<strong>en</strong>te organizada u ocasionalm<strong>en</strong>te <strong>de</strong> reflexiones hiperbólicas<br />
<strong>en</strong>trecruzadas. Este tipo <strong>de</strong> apari<strong>en</strong>cia es consist<strong>en</strong>te ya sea con el substrato rocoso o con una litología,<br />
52
XIII Congreso Peruano <strong>de</strong> Geología. Resúm<strong>en</strong>es Ext<strong>en</strong>didos<br />
Sociedad Geológica <strong>de</strong>l Perú<br />
<strong>en</strong> <strong>la</strong> cual gran<strong>de</strong>s bloques son <strong>la</strong> mayoría <strong>de</strong> compon<strong>en</strong>tes. Des<strong>de</strong> el punto <strong>de</strong> vista geológico<br />
repres<strong>en</strong>ta <strong>la</strong> transición substrato-caos granítico.<br />
En los perfiles georadar no se observó ninguna fal<strong>la</strong> principal o p<strong>la</strong>no <strong>de</strong> <strong>de</strong>slizami<strong>en</strong>to a lo <strong>la</strong>rgo <strong>de</strong><br />
cualquiera <strong>de</strong> <strong>la</strong>s líneas <strong>de</strong> estudio, aunque hay indicaciones <strong>de</strong> fracturami<strong>en</strong>to rocoso <strong>de</strong> m<strong>en</strong>or<br />
importancia con infiltraciones <strong>de</strong> agua tanto <strong>en</strong> <strong>la</strong> carretera Hiram Bingham, así como <strong>en</strong> <strong>la</strong> Ciuda<strong>de</strong><strong>la</strong><br />
Inca.<br />
Fig. 5. Ubicación <strong>de</strong> los perfiles <strong>de</strong> georadar y perfiles interpretados<br />
CONCLUSIONES<br />
De los tres métodos empleados el electromagnético no fue efectivo para el mapeo <strong>en</strong> Machu Picchu.<br />
La resistividad fue efectiva para <strong>la</strong>s rocas graníticas, pues es capaz <strong>de</strong> mapear variaciones d<strong>en</strong>tro <strong>de</strong><br />
rocas altam<strong>en</strong>te resistivas. El georadar <strong>de</strong> p<strong>en</strong>etración GPR fue efectivo para <strong>de</strong>terminar el espesor <strong>de</strong>l<br />
rell<strong>en</strong>o, <strong>la</strong> litología y <strong>la</strong> topografía <strong>de</strong>l sustrato.<br />
No se ha observado ninguna fal<strong>la</strong> principal o p<strong>la</strong>nos <strong>de</strong> <strong>de</strong>slizami<strong>en</strong>to a lo <strong>la</strong>rgo <strong>de</strong> cualquiera <strong>de</strong> <strong>la</strong>s<br />
líneas o perfiles <strong>de</strong> estudio, aunque hay indicaciones <strong>de</strong> fracturas <strong>de</strong> m<strong>en</strong>or importancia rell<strong>en</strong>adas con<br />
agua <strong>en</strong> <strong>la</strong> carretera Hiram Bingham y cerca <strong>de</strong> <strong>la</strong> p<strong>la</strong>za principal. El relieve <strong>de</strong>l substrato bajo <strong>la</strong> p<strong>la</strong>za<br />
es significativo y <strong>la</strong> parte superior parece haber sido rell<strong>en</strong>ada por fragm<strong>en</strong>tos y escombros. A<strong>de</strong>más el<br />
georadar permite interpretar bloques gran<strong>de</strong>s <strong>de</strong> granito <strong>en</strong> <strong>la</strong> parte superior <strong>de</strong>l substrato que<br />
correspon<strong>de</strong>ría al caos granítico.<br />
Las especu<strong>la</strong>ciones por otras investigaciones que un p<strong>la</strong>no <strong>de</strong> <strong>de</strong>slizami<strong>en</strong>to profundo <strong>de</strong> dirección<br />
norte-sur podría pres<strong>en</strong>tarse bajo <strong>la</strong> P<strong>la</strong>za Principal no están sust<strong>en</strong>tadas por el radar <strong>de</strong> p<strong>en</strong>etración ni<br />
por <strong>la</strong> resistividad.<br />
REFERENCIAS<br />
Carlotto, V. Cárd<strong>en</strong>as, J. Romero, D. Valdivia, W. Tintaya, D. (1999). Geología <strong>de</strong> los Cuadrángulos <strong>de</strong><br />
Quil<strong>la</strong>bamba y Machu Picchu. Boletín Nº 127, Serie “A”: Carta Geológica Nacional. INGEMMET, 317 p.<br />
Mucho, R. Carlotto, V. Pari, W. Oviedo, M. Douma, M. Best, M. Bobrowsky, P.(2005). The application of<br />
Ground p<strong>en</strong>etrating Radar (GPR) at Machu Picchu, Peru (C101-1) En Landsli<strong>de</strong>s Risk Analysis and<br />
Sustainable Disaster Managem<strong>en</strong>t I Asamblea G<strong>en</strong>eral Consorcio Internacional <strong>de</strong> Deslizami<strong>en</strong>tos-ICL-<br />
Washington , 2005: 55-59.<br />
Wright KR, Val<strong>en</strong>cia, A. (1999) Anci<strong>en</strong>t Machu Picchu drain age <strong>en</strong>gineering. J Irrig Drain E-ASCE<br />
125:360-369.<br />
53