1.
ZONIFICACION SISMICA
1.1. Introduccion
El borde occidental de América del Sur es una de las regiones
sísmicamente más activas en el mundo y su actividad más importante esta
asociada al proceso de subducción de la placa de Nazca por debajo de la
Sudamericana. El Perú es parte de esta región y frente a su línea de costa se
genera de manera frecuente sismos de diversas magnitudes a diferentes
niveles de profundidad, siendo los más grandes los que han generado
destrucción importante en los diversos departamentos distribuidos, de norte a
sur, cercanos a la costa. Un segundo tipo de sismicidad es la producida por las
deformaciones corticales, presentes sobre la Cordillera Andina y a lo largo de la
zona denominada Subandina, en donde se producen sismos menores en
magnitud y frecuencia. La ocurrencia continua de sismos en el Perú, como
en cualquier otra región sísmicamente activa, juegan un rol importante en los
diversos procesos sismotectónicos que se producen en el tiempo como la
deformación de rocas, formación de fallas, fenómenos precursores de
grandes sismos y la ruptura sísmica que en general, tiene duraciones del
orden de 60 segundos para un sismo de magnitud Ms igual a 7.0.
1.2.- distribucion espacial de los sismos
1.2.1. - Sismicidad Histórica
La información disponible sobre la sismicidad histórica de Perú
viene del tiempo de la conquista y colonización del país, y se encuentra
esparcida en diferentes obras inéditas, manuscritos, crónicas, narraciones,
informes administrativos de aquellos años .
Gran parte de esta información ha sido recolectada y publicada por Silgado
(1978) y Dorbath et al (1990). Es importante considerar que la calidad de estos
datos depende de la distribución y densidad de la población en las regiones
afectadas en el pasado; por lo tanto, es posible que hayan ocurrido sismos
importantes en áreas no pobladas o próximas a localidades con las
cuales era difícil establecer comunicación. Esto explicaría la ausencia de
información sobre sismos que pudieran haber ocurrido en la Alta Cordillera y
zona Subandina.
En la Figura , se muestra la
localización epicentral de los
sismos históricos ocurridos
entre 1513 y 1959 (Imax >
VII en la escala Mercalli
Modificada), los mismos que
en su mayoría se
distribuyen principalmente
entre la línea de fosa y la
costa, siendo mayor el
número de sismos en la
región Centro y Sur de Perú
debido probablemente a que
estas regiones eran las mas pobladas y donde se constituyeron las ciudades
más importantes después del siglo XVI.
Los sismos de 1586 (IX MM), primer gran sismo para el cual se tiene
documentación histórica; 1687 (VIII MM) y 1746 (X MM) destruyeron casi
completamente a la ciudad de Lima. Este último generó un tsunami con olas de
15-20 metros de altura que inundo totalmente al puerto del Callao. Asimismo,
los sismos de 1604 (IX MM), 1784 (X MM) y 1868 (X MM) destruyeron
parcialmente a las ciudades de Arequipa, Moquegua, Tacna, Puno y Norte de
Chile. Según la información disponible, estos sismos habrían producido daños
moderados en la ciudad de Puno.
TALLER DE OBRA HIFRAULICA - pág. 2
Para este periodo de tiempo (1500-1959), no existe información sobre los
sismos ocurridos en la Alta Cordillera y menos para el departamento de
Puno, pero es posible que el sismo de 1650 haya sido sentido con intensidades
del orden de V (MM) en Ilave con daños menores.
En el mapa ya mostrado se indica el área de estudio y en ella los sismos
históricos que también habrían producido daños en el departamento de Puno y
cuyo detalle.
1.2.2.- Sismicidad Instrumenta
El análisis y evaluación de la distribución espacial de los sismos ocurridos en
Perú durante el periodo instrumental, ha sido realizado por diversos autores
utilizando información telesísmica (Cahill y Isacks, 1992; Tavera y Buforn, 2001
y Bernal, 2002) y regional (Hasegawa y Sacks, 1981; Grange, 1984;
Schneider y Sacks, 1987; Rodríguez y Tavera, 1991) y los resultados
obtenidos han permitido conocer las características de la geometría de las
principales fuentes sismogénicas asociadas a la deformación cortical y al
proceso de subducción de la placa de Nazca bajo la Sudamericana.
TALLER DE OBRA HIFRAULICA - pág. 3
La sismicidad en el borde oriental de la Cordillera Andina es mayoritaria
en las regiones norte y centro de Perú; mientras que, en la región sur el número
de sismos es menor debido posiblemente a que la cordillera presenta mayor
raíz (75 km) y espesor (200 km), imposibilitando el desarrollo de rápidas
deformaciones y la ocurrencia frecuente de sismos. En la zona andina de la
región sur, la ocurrencia de sismos es mínima, aunque en la frontera de Puno
con Arequipa, se nota la presencia de un alineamiento de sismos debido
probablemente a una pequeña falla presente cerca de la localidad de Santa
Lucia y que fuera reactivada con el sismo de Arequipa del 2001 con un sismo
de magnitud ML=5.0 (29 de Junio 2001; 22h 33 min HL) que causó alarma en
las localidades de Santa Lucia y Lampa porque produjeron deslizamientos de
tierra y piedras, además de emanaciones de agua. La intensidad evaluada en
la ciudad Ilave fue de IV (MM), en Juliaca y Ayaviri de III (MM).
TALLER DE OBRA HIDRAULICA
En general, la distribución de sismos en todo el territorio peruano permite
identificar la presencia de diversas fuentes sismogénicas y de estas, las más
importantes son las que producen los sismos de elevada magnitud a niveles
superficiales: proceso de subducción y deformación cortical. El potencial de
estas fuentes sismogénicas puede ser cuantificado con la cantidad de energía
liberada por los sismos ocurridos en dichas fuentes. Así, en la Figura 3 se
puede observar que la principal fuente configura un área paralela a línea de
costa en donde se observa la presencia de acumulación de valores altos de
energía producto de la ocurrencia de grandes sismos debidos al proceso de
contacto de placas (placa de Nazca y Sudamérica) y que fue descrito
anteriormente. Una segunda fuente sismogénica con valores altos de energía
liberada se ubica en el interior del continente y se distribuye paralela a la
Cordillera Andina y tiene su origen en los sismos que ocurren debido a la
deformación cortical puesta en evidencia con la presencia de importantes
sistemas de fallas geológicas. Obsérvese que en los alrededores de la ciudad
de Puno existe una pequeña concentración de valores de energía que indican
la ocurrencia de sismos con magnitudes menores a 5.0Ms.
TALLER DE OBRAS HIDRAULICAS
2.3.- deformacion continenta
En el caso de la deformación cortical, la información sísmica disponible permite
considerar que los periodos de recurrencia de los sismos de magnitud elevada en
el interior del continente, aparentemente son muy grandes, ya que para los últimos
500 años no se cuenta con información confiable que permita confirmar la
ocurrencia de sismos. Esto sugiere que la deformación continental es lenta y poco
uniforme. Los pocos sismos para los cuales se cuenta con información han
permitido observar en superficie rasgos tectónicos que indican el tipo de
deformación producida (fallas geológicas). En general, las zonas de mayor
deformación superficial han sido identificadas y descritas por diversos autores
desde el punto de vista sismológico y tectónico. Recientemente, Tavera y Buforn
(2001) han estudiado la distribución de los esfuerzos que causan la deformación
en el interior del continente a partir de la elaboración de mecanismos focales para
CIUDAD DE ESTUDIO
ZONA DE
ESTUDIO
SISTEMAS DE
sismos con magnitudes MsFALLAS PRESENRES
mayores a 5.8. EN LA REGION DEL SUR
TALLER DE OBRAS HIDRAULICAS
Los resultados indican que la deformación es mayor en la zona Subandina de
la región norte en comparación con otras regiones y con la Alta Cordillera. La
orientación de los esfuerzos causantes de la deformación superficial es
perpendicular a la dirección que sigue la Cordillera Andina, tal como se
observa en la Figura 5. En la región sur, la ausencia de sismos de magnitud
elevada sustenta el desarrollo de un lento proceso de deformación que
conlleva a la ocurrencia de sismos con periodos de retorno de muchas
décadas.
TALLER DE OBRAS HIDRAULICAS
1.4.- zonificacion sismica
El mapa de zonificación sísmica, base primaria para que el Gobierno
del Perú en el año 1997 y modificado en el año 2003 elabore la Norma
Sismorresistente NTE-030-97, En la elaboración de este mapa se ha
considerado las condiciones regionales de sismicidad que engloban a todos
los parámetros sísmicos significativos como son la ubicación de los focos
sísmicos, trayectoria de recorrido de las ondas, distancia epicentral y
energía liberada, todos contenidos y definidos en términos de la aceleración
del suelo producido por los sismos y que son medidos sobre roca o
terreno firme; a lo que se agregan las amplificaciones que sufren las
ondas sísmicas por las condiciones locales.
1.4.1. Parametros sismicos
Calculo del factor de zona (z)
Ubicando la zona de nuestro sitio, Puno
TALLER DE OBRAS HIDRAULICAS
2. ANALISIS DE PELIGRO SISMICO Y ESTIMADO DEL MOVIMIENTO
SISMICO DE DISEÑO
2.1. Introduccion
La ocurrencia de terremotos desastrosos nos reafirma la importancia de los
Análisis de Peligro y Riesgo Sísmicos para estimar las consecuencias de
Estos eventos. Aún cuando se han logrado grandes avances en la predicción
sísmica, el tiempo de ocurrencia, la magnitud o la ubicación de un terremoto no
pueden aún ser predichos con certeza. Por lo tanto un adecuado análisis de
peligro y riesgo sísmico, si bien no podrá eliminar el daño potencial, ayudará a
reducir sus efectos considerablemente al realisar el diseño de la presa.
2.2. PELIGRO SISMICO
El Análisis de Peligro Sísmico involucra la estimación cuantitativa de la
amenaza sísmica en un determinado lugar. Este peligro sísmico puede ser
analizado determinísticamente, asumiendo un determinado terremoto
escenario, o probabilísticamente, considerando explícitamente las
incertidumbres en el tamaño, ubicación y tiempo de ocurrencia de los
terremotos. Por otro lado, el Análisis de Riesgo Sísmico permite estimar los
daños y pérdidas potenciales en una región debido a terremotos futuros; es
decir incluye la evaluación del comportamiento de las estructuras ante la
amenaza sísmica de la zona donde éstas se ubican.
Los pasos básicos en un procedimiento de Análisis de Peligro y Riesgo
Sísmico regional incluyen:
Identificación de las
fuentes sismogénicas
Modelamiento de la ocurrencia de terremotos en estas fuentes
Estimación de la atenuación del movimiento sísmico entre la fuente y la
región
Evaluación de los efectos de locales de sitio como amplificación
sísmica, licuación, deslizamientos y manifestación superficial de la
ruptura de la falla.
Estimación de daños en los inventarios regionals
Estimación de las pérdidas esperadas en la región.
TALLER DE OBRAS HIDRAULICAS
2.2.1 I d e n t i f i c a c i o n y evaluacion de las fuentes sismogenicas
Para evaluar el peligro sísmico en un determinado lugar o región, se
deben identificar todas las posibles fuentes de actividad sísmica y evaluar su
potencial para generar futuros terremotos. La identificación de las fuentes
sismogénicas requiere de un cierto trabajo de detección, observando e
identificando rasgos
a) Evidencias Geológicas
La teoría de la tectónica de placas nos asegura que la ocurrencia de los
terremotos está escrita en los registros geológicos, principalmente en forma
de dislocamientos o desplazamientos de varios estratos. En algunas partes del
mundo, estos registros geológicos son fácilmente accesibles y relativamente
fáciles de interpretar por geólogos especialistas.
b) Actividad de las Fallas
La mera presencia de una falla no indica la probabilidad de ocurrencia de
futuros terremotos. El concepto de actividad de las fallas es importante y ha
sido un tópico de mucha discusión y controversias a través de los años.
c) Indicadores de Magnitud
La evidencia geológica también puede ser usada para estimar la magnitud de
terremotos pasados, mediante correlación de las características de
deformación observadas con las magnitudes de terremotos registrados.
Estudios de terremotos en todo el mundo muestran que las fallas no se
rompen en toda su longitud durante un evento individual.
d) Evidencias Tectónicas
La tectónica de placas y la teoría del rebote elástico nos dicen que los
terremotos ocurren para liberar la energía de deformación que se acumula tras
el movimiento relativo de las placas. La razón de movimiento por lo tanto
estaría relacionada a la razón de acumulación de energía de deformación
y también a la razón de liberación de energía de deformación.
TALLER DE OBRAS HIDRAULICAS
2.3 RIESGO SISMICO
El siguiente estudio de riesgo sísmico es valido para todas las obras del
proyecto, de acuerdo a lo que se realiza para este tipo de estudios y proyectos.
2.3.1 Análisis Estadístico e Información histórica sísmica:
Seguidamente se va a indicar del catalogo sísmico usado para la evaluación
del Riesgo Sísmico de la zona de estudio, una tabla resumen con los eventos
sísmicos que se han presentado cada año los cuales presentaron magnitudes
mayores a 4.50, para lo cual se han considerado datos correspondientes a un
área de radio 400 km. tomando como centro la zona de estudio coordenadas
13.808º Latitud Sur y 73.233º Longitud Oeste.
2.3.2. Periodo de Retorno
Para la evaluación del período de retorno, se realizó un análisis estadístico de
la sismicidad del área del proyecto, para lo cual se toma en consideración lo
siguiente:
- Distribución Probabilística de Poisson
- Probabilidad de Ocurrencia
- Periodo de Retorno
1) Distribución Probabilística De Poisson
La distribución de Poisson, permite determinar si los datos utilizados son los
adecuados para este tipo de estudios. El método consiste en plotear los
números de sismos observados (N) de frecuencia de clases, teniendo 0,1, 2,
…,n sismos por intervalos de tiempo (t) y comparando con la distribución
aleatoria de Poisson (P). Los datos utilizados corresponden a sismos
superficiales y intermedios correspondientes a los años 1973-2008 con
magnitud ≥ 4.0 Mw.
TALLER DE OBRAS HIDRAULICAS
La ecuación que se utiliza es:
Donde:
n = 1,2,3,…. n
t = período (para un año)
V = Tasa de ocurrencia por “t” debe ser un año
Donde V = N° de sismos
N° de Años
K(N = n, t) = k.P(N = n, t)
K= cálculo del número esperado de eventos
k = numero de intervalos de tiempo usado (Nº de sismos)
Para esto se determina el número de eventos que se presenta por cada año
o Ocurrencia de sismos desde 1973-2008, para magnitudes ≥ 4.5 Mw
Cálculo del número esperado de eventos (K) (Poisson) y su comparación con
los datos originales
Frecuencia de Nº de años con Nº total de P K
n clases n eventos eventos
0 0 0 3.79763E-09 1.36715E-07
1 0 0 7.36318E-08 2.65075E-06
2 0 0 7.13819E-07 2.56975E-05
3 0 0 4.61339E-06 0.000166082
4 0 0 2.23621E-05 0.000805036
TALLER DE OBRAS HIDRAULICAS
5 0 0 8.67153E-05 0.003121752
6 0 0 0.000280219 0.010087884
7 0 0 0.000776162 0.027941838
8 0 0 0.001881115 0.067720148
9 0 0 0.004052526 0.145890936
10 0 0 0.007857398 0.282866314
11 1 11 0.013849655 0.498587594
12 2 24 0.022377452 0.805588289
13 2 26 0.033374918 1.201497064
14 2 28 0.046221613 1.663978076
15 2 30 0.059745715 2.150845736
16 1 16 0.072400189 2.606406812
17 4 68 0.082574072 2.972666593
18 1 18 0.088945528 3.202039015
19 4 76 0.090766051 3.267577825
20 2 40 0.087992644 3.167735169
21 2 42 0.081241885 2.924707868
22 2 44 0.07159954 2.577583449
23 4 92 0.060358067 2.172890396
24 3 72 0.048761494 1.755413769
25 0 0 0.037817247 1.361420901
26 2 52 0.028201323 1.015247637
27 0 0 0.020251568 0.729056431
28 1 28 0.014023407 0.504842648
29 0 0 0.009375803 0.337528897
30 0 0 0.006059547 0.218143676
31 1 31 0.003789931 0.136437532
Luego se determina el valor de V: V = 19.39
TALLER DE OBRAS HIDRAULICAS
Cálculo de los valores para la utilización del Método de Max-Likelihood.
Por el método de mínimos cuadrados obtenemos “a” y “b” mediante las
siguientes expresiones:
a = 5.01
b = - 0.80
Donde:
N = Numero de intervalos “P”
Luego la ecuación lineal es: Y = a + bX
Donde:
Y = Numero de sismos
X = Magnitud
2) Probabilidad de ocurrencia
TALLER DE OBRAS HIDRAULICAS
Para esto se tiene que hacer una corrección a la anterior ecuación,
obteniéndose una ecuación para un año, en base a nuestro catálogo de 34
años que es el tiempo t;
Y1 = 10 5.01 - 0.80 x
a1 = a – log(tc)
a1 = 5.01– log(36)
a1 = 3.45
Obteniéndose la siguiente ecuación:
Y1 = 10 3.45 -0.80 x
Ahora podemos usar la ecuación de Poisson para la probabilidad de
ocurrencia:
Pr (M) = 1 – exp(-Y1T)
Cálculo de la Probabilidad de Ocurrencia
Tiempo(añ Probabilidad de ocurrencia %
os) Mw=5 Mw=6 Mw=7 Mw=7.5
T=1 1 24.56043 4.368541 0.705445 0.281441
34 78 75 5
T=2 2 43.08871 8.546241 1.405914 0.562090
8 99 97 91
T=5 5 75.56592 20.01585 3.477813 1.399308
39 8 22 84
T=10 10 94.02975 36.02537 6.834674 2.779037
92 03 59 03
T=15 15 98.54122 48.83044 10.07479 4.139458
68 13 06 55
T=20 20 99.64356 59.07246 13.20222 5.480843
23 75 14 58
T=30 30 99.97871 73.81676 19.13456 8.107565
98 27 71 94
T=40 40 99.99872 83.24937 24.66145 10.66129
TALLER DE OBRAS HIDRAULICAS
95 09 63 07
T=50 50 99.99992 89.28384 29.81060 13.14404
41 7 06 65
T=100 100 100 98.85164 50.73448 24.56043
07 21 34
Probabilidad de Ocurrencia
120
100
Probabilidad en porcentaje
80
60
40
20
0
0 20 40 60 80 100 120
Tiempo en Años
5 Mw 6 Mw 7 Mw 7.5 Mw
3) Período de retorno
El período medio de retorno presenta la siguiente ecuación: θ = 1/y
Periodo de Retorno
Mw Y1 O
4.5 0.707945784 1.4
5.0 0.281838293 3.5
5.5 0.112201845 8.9
6.0 0.044668359 22.4
6.5 0.017782794 56.2
7.0 0.007079458 141.3
7.5 0.002818383 354.8
TALLER DE OBRAS HIDRAULICAS
400.0
Periodo de Retorno
350.0
300.0
Tiempo en años
250.0
200.0
150.0
100.0
50.0
0.0
4 4.5 5 5.5 6 6.5 7 7.5 8
Magnitud (Mw)
2.3.2. Estimación de la Aceleración de la Zona
Definidas las fuentes de sismicidad, se determina la mayor solicitación que
representa el terremoto máximo posible, a considerar en la seguridad de las
obras.
a) Factores de Atenuación
Para sismos originados en zona de subducción, tipo Perú y profundidades
focales superiores a 6 km., según Woorward and Clyde (1982):
PGA = 0.214 ( e 0.5Mw) x ( R + C)-0.85
o PGA = Aceleración de pico horizontal en g.
o MW = Magnitud-momento, equivale a magnitud MS hasta 8.0
o R = Distancia inclinada en Km. hasta el sector del plano de ruptura.
o C= 0.864 e0.463Mw, Fundación en rocas firmes.
Para sismos superficiales originados en la corteza continental: según Campbell
(1981):
PGA = 0.016 (e 0.87M) x ( R + C)-1.09
o PGA = Aceleración de pico horizontal en g.
o M = Magnitud (en general ML magnitud local)
o R= Distancia a la falla en km.
TALLER DE OBRAS HIDRAULICAS
o C = 0.06 e0.7 M
Supone que en la fuente el sismo de magnitud máxima ocurre en cualquier
parte con igual probabilidad. Con la magnitud máxima y distancia mínima, se
halla el valor de aceleración del sismo máximo. Considerando varias fuentes, el
terremoto máximo creíble es el mayor de los valores calculados:
Valores de aceleracion maxima (pga, g)
PRESA
FUENTE M máx. D mín. (KM) PGA (g)
Benioff
M máx. = 8.5 100 0.22
Terremoto maximo creible
Es el valor máximo de solicitación sísmica a considerar en el cálculo
sismorresistente de obras, cuya integridad es imprescindible por la seguridad
de las poblaciones Este sismo ha sido calculado, para un sismo de magnitud
8.5 Mw, cuyo valor resultante de aceleración es de 0.22g.
Terremoto de operación de la presa
Se determina por imperativos económicos, la ocurrencia del terremoto que no
perjudique la operación normal y la rentabilidad de la obra. Se analiza la
utilidad, costo del resguardo antisísmico, vulnerabilidad en relación al nivel de
solicitación, perjuicio socio económico de una interrupción de operación, costo
de reparación, etc.
Se considera el periodo de retorno de 50 años para el terremoto de operación,
calculándose con una magnitud de 6.5, habiendo resultado el valor más alto de
0.10g, que debe aplicarse como sismo de operación de la presa, para los
respectivos cálculos.
TALLER DE OBRAS HIDRAULICAS
Valores de aceleracion maxima (pga, g)
FUENTE M máx. 6.50 PRESA
Periodo de Retorno 50 D mín. PGA (g)
años (KM)
Benioff
M máx. = 6.50 100 0.10
Infromacion general
Descripción de los mapas de Isositas en el Perú.
TALLER DE OBRAS HIDRAULICAS
TALLER DE OBRAS HIDRAULICAS
Zona de
Estudio
Mapa de distribución de máximas intensidades
sísmicas observadas.
TALLER DE OBRAS HIDRAULICAS
ZONA DE ESTUDIO
Placa de Nazca
Placa
Continental
Subducción
Los sismos y las placas tectónicas Nazca y Continental –
Proceso de Subducción.
PARÁMETROS DE DISEÑO
De esta manera, el territorio nacional queda segmentado tal como se muestra
en las figuras siguientes, donde se observa que el proyecto se encuentra
ubicado dentro del área de influencia de la Zona 2, se espera que el factor de
aceleración máxima del terreno (A) sea de 0.34 g, con una probabilidad de 10%
de ser excedida en 50 años, equivalente a un periodo de recurrencia de 475
años.
TALLER DE OBRAS HIDRAULICAS
TALLER DE OBRAS HIDRAULICAS
ACTUALIZACIÓN DEL PROGRAMA DE CÓMPUTO ORIENTADO A
LA DETERMINACIÓN DEL PELIGRO SÍSMICO EN EL PAÍS
Probabilidad anual de excedencia
La primera pestaña de análisis es para las Curvas de “Probabilidad anual de
excedencia VS. Aceleración espectral”, donde predeterminadamente se
analiza para un amortiguamiento de 5% y un periodo de 0.00 (PGA).
Esta gráfica tiene por abscisas a la Probabilidad anual de excedencia que es la
inversa del Tiempo de retorno en años, y a la Aceleración espectral dadas en
unidades gal en las ordenadas.
Para el eje de las abscisas el límite máximo es 1.00, ya que no se requiere
análisis para periodos de retorno menores a 1 año.
TALLER DE OBRAS HIDRAULICAS
Espectro de diseño
En esta sección, la gráfica presentada es el Espectro de Peligro Uniforme,
donde de manera predeterminada aparece la gráfica para un periodo de
retorno de 75 años.
En el eje de las abscisas se tienen los valores de aceleración espectral en
unidades g, y en las ordenadas están los periodos estructurales de diseño.
Los periodos estructurales para los cuales se analizó desde 0.000, siguiendo
0.050, y 0.010, a partir del cual se procede en intervalos de 0.05 s hasta el
TALLER DE OBRAS HIDRAULICAS
periodo 1.0 s, donde el intervalo cambia a 0.1 s hasta llegar al periodo
máximo de análisis de 3.0 s.
TALLER DE OBRAS HIDRAULICAS
