CAPITULO IV

FUENTES SISMICAS Y LEYES DE

ATENUACIÓN PARA EL PERU

4.1 Introducción

Para analizar probabilísticamente el peligro sísmico es necesario

conocer la contribución que cada fuente sísmica pueda realizar al sitio en
estudio. Para analizar la contribución de cada fuente es necesario distinguir
además de su geometría, la magnitud mínima y máxima que cada fuente
pueda desarrollar, seleccionar una profundidad representativa de los
eventos sísmicos dentro de la fuente y reconocer el tipo de mecanismo que
gobierna el movimiento del suelo dentro de la fuente para asignarle un
modelo de atenuación.

En este capitulo se analizan los métodos más usados para describir los
parámetros que se utilizan en un análisis de peligro sísmico. Se desarrollan
los parámetros para una fuente sísmica y se describen los modelos de
atenuación elegidos para representar el movimiento del suelo y respuesta
estructural. Los valores del movimiento del suelo y respuesta estructural
obtenidos de los modelos de atenuación se comparan luego con algunos
valores medidos en el suelo firme de Lima.

69
4.2 Fuentes sismogénicas
Las fuentes sismogénicas utilizadas para representar la sismicidad del País
corresponden a las propuestas por Castillo y Alva (1993). Castillo y Alva
definieron 20 fuentes sismogénicas basados en la distribución espacial de la
sismicidad y en los aspectos neotectónicos del Perú. Las fuentes fueron
agrupadas en fuentes de subducción y corteza.

La interacción de las placas de Nazca y Sudamericana están representadas

por las fuentes de subducción superficial (0-70 km), fuentes con
profundidad intermedia (71-300 km) y fuentes que representan la
sismicidad profunda (500 - 700 km).

La subducción superficial esta representada por las fuentes 1 a la 5

ubicadas a lo largo de la costa. Las fuentes de corteza están representadas
por las fuentes 6 a la 12 (Ver figura 4.1). Las fuentes con profundidad
intermedia están representadas por las fuentes 13 a la 19. La sismicidad
profunda la representa la fuente 20 (Ver figura 4.2).

Todas las fuentes sísmicas fueron modeladas como fuentes volumétricas,

debido a la insuficiencia de datos para modelar fallas como fuentes lineales.
(Castillo y Alva 1993).

70
 -82 -80 -78 -76 -74 -72 -70 -68
+2 +2

F1 F6
0 0

-2 -2

-4 F10 -4

-6 -6

F2 F11
-8 F7 -8

-10 -10
F3
-12 -12
F12
-14 -14

F8
-16
F4 -16
F9

-18 -18
FUENTES SUPERFICIALES
H < 70 km
Castillo-Alva (1993)
-20 Catálogo SISRA (1963-2000) Mag. F5 -20
4.0
km 5.0
100 200 300 400 6.0
7.0
-22 -22
-82 -80 -78 -76 -74 -72 -70 -68

Figura 4.1. Ubicación de las fuentes para representar la sismicidad superficial (h <
70 km) para la zona de subducción (F1 a F5) y para la zona de corteza (F6 a F12).

71
 -82 -80 -78 -76 -74 -72 -70 -68
+2 +2

F13
0 0

-2 -2

F18
-4 -4

F14
-6 -6

-8 -8
F15 F20
-10 -10

F19
-12 -12

-14 -14
F16
-16 -16
F17

-18 -18
F. INTERMEDIAS Y PROFUNDAS
71 < H < 700 km
Castillo - Alva (1993)
-20 Catálogo SISRA (1963-2000) Mag. -20
4.0
km
5.0
100 200 300 400 6.0
7.0
-22 -22
-82 -80 -78 -76 -74 -72 -70 -68
Figura 4.2. Ubicación de las fuentes para representar la sismicidad intermedia
(70< h < 300 km) de la zona de subducción (F13 a F19) y sismicidad profunda
(500< h < 700km) (F20).

72
Sin embargo, la variación en la distribución espacial de sismos presentada
en los últimos 10 años, principalmente en la costa Peruana, hace necesaria
la modificación de las fuentes de subducción F3 y F4 y de la fuente de
corteza F8 propuestas por Castillo y Alva (1993). Esta variación se evidencia
de manera más clara en los cortes transversales C4, C5, C6 y C7 que se
encuentran en el Anexo 1.

La figura 4.3 muestra uno de los cortes realizado sobre la ciudad de Lima
(corte C6) con los eventos proyectados 100 km a cada lado del corte, donde
se ha graficado la proyección en profundidad de la fuente F3 propuesta por
Castillo y Alva (1993) y la modificación propuesta en este trabajo sobre la
línea de corte.

Figura 4.3. Distribución de sismos en profundidad donde se muestra la

proyección de la fuente F3 propuesta por Castillo y Alva (1993) y la proyección de
la fuente F3 propuesta en este trabajo. Nótese el incremento de sismos en la zona
de acoplamiento entre la placa de nazca y Sudamericana.

En la figura se evidencia un grupo de sismos que sigue la pendiente de la

zona de subducción, en el límite de interacción entre las placas, que no es
considerado en la fuente F3 propuesta por Castillo. Este incremento de

73
sismos hace necesario la modificación de la fuente F3. La modificación de
las fuentes F4 y F8 se hizo con el mismo criterio.

La figura 4.4 muestra la ubicación de las nuevas fuentes superficiales

propuestas en este trabajo. Las fuentes que representan la sismicidad
intermedia y profunda se indican en la figura 4.2. Las coordenadas
geométricas que definen la forma y ubicación de las fuentes sísmicas se
detallan en la Tabla 4.1.

74
 -82 -80 -78 -76 -74 -72 -70 -68
+2 +2

0 0

F1
-2 -2
F6
-4 F10 -4

-6 -6
F11
F2 F7
-8 -8

F3
-10 -10

-12 -12

F12
-14 -14

F4 F8
-16
F9 -16

F5
-18 -18
FUENTES SUPERFICIALES
H < 70 km
Monroy - Bolaños (2004)
Mag.
-20 Catálogo SISRA (1963-2000) -20
4.0
km
5.0
100 200 300 400 6.0
7.0
-22 -22
-82 -80 -78 -76 -74 -72 -70 -68

Figura 4.4. Ubicación de las fuentes para representar la sismicidad superficial (h <
70 km) de la zona de subducción (F1 a la F5) y para la zona de corteza (F6 a la
F12) propuestas en este trabajo.

75
 Tabla 4.1.- Coordenadas geográficas de las fuentes de análisis
COORDENADAS GEOGRAFICAS
FUENTE LONGITUD LATITUD LONGITUD LATITUD
1 -80,29 2,00 -78,32 2,00
-81,39 -0,97 -79,65 -1,21
-81,52 -2,39 -80,19 -2,50
2 -82,00 -3,39 -80,17 -3,45
-82,00 -6,83 -80,67 -5,42
-81,17 -9,00 -79,27 -7,90
3 -81,17 -9,00 -79,27 -7,90
-77,00 -14,80 -75,54 -13,755
4 -77,00 -14,80 -75,54 -13,755
-74,16 -17,87 -73,00 -16,53
5 -74,16 -17,87 -73,00 -16,53
-71,85 -19,87 -69,21 -19,00
-71,85 -22,00 -69,21 -22,00
6 -77,50 1,58 -76,92 1,19
-79,83 -1,65 -78,90 -2,53
-79,96 -2,46 -78,97 -3,43
-80,92 -2,96 -80,79 -3,44
7 -78,28 -8,20 -77,86 -8,07
-77,21 -10,47 -76,83 -10,23
8 -75,54 -13,755 -74,76 -13,13
-73,00 -16,53 -71,41 -14,67
9 -73,00 -16,53 -71,41 -14,67
-69,71 -18,67 -68,12 -16,13
10 -76,92 1,19 -76,50 1,00
-78,90 -2,53 -77,35 -2,40
-79,10 -5,20 -77,00 -4,77
11 -79,10 -5,20 -75,10 -4,33
-76,34 -10,67 -74,17 -9,33
-74,76 -13,13 -72,48 -11,40
12 -74,76 -13,13 -72,48 -11,40
-68,12 -16,13 -67,76 -13,80
13 -78,73 2,00 -76,00 1,82
-81,00 -0,67 -79,59 -2,55
-81,00 -3,07 -79,20 -3,07

76
 FUENTE COORDENADAS GEOGRAFICAS
LONGITUD LATITUD LONGITUD LATITUD
-81,00 -3,07 -79,20 -3,07
14 -81,93 -5,73 -78,60 -4,00
-79,80 -8,13 -77,17 -6,53

15 -79,80 -8,13 -77,17 -6,53

-76,38 -14,30 -73,86 -12,46
16 -76,38 -14,30 -73,86 -12,46
-73,28 -16,87 -71,21 -14,40
17 -73,28 -16,87 -71,21 -14,40
-70,86 -18,80 -68,93 -15,73
-70,38 -22,00 -67,98 -22,00
18 -79,59 -2,55 -77,50 -0,73
-78,60 -4,00 -75,51 -2,06
-77,17 -6,53 -75,27 -5,33
19 -77,17 -6,53 -75,27 -5,33
-73,86 -12,46 -72,03 -11,13
20 -72,31 -6,67 -71,00 -6,33
-71,14 -11,30 -69,69 -10,93

4.3 Evaluación del tamaño de los eventos sísmicos

Para evaluar la variación en el tamaño de los eventos sísmicos que cada
fuente sísmica pueda generar es necesario conocer la recurrencia sísmica de
la fuente. La recurrencia sísmica representa el número de eventos mayores
o iguales a alguna magnitud dentro de la fuente y está descrita por la
pendiente de la relación de recurrencia de Gutenberg y Ritcher (b), la razón
media anual de actividad sísmica (), la magnitud mínima y la magnitud
máxima.

Las ecuaciones para definir la recurrencia sísmica fueron desarrolladas en el

capitulo II. A continuación se presentan algunas ecuaciones por
conveniencia. La relación de recurrencia de Gutenberg y Ritcher es la
siguiente:
Log Nm = a – bm

77
De la que se puede obtener el número de sismos Nm con la ecuación
(2.11):
exp  (m  mmin )  exp  (mmax  mmin )
Nm  
1  exp  (mmax  mmin )

Donde:
-m
=e min

 = b*ln(10)
m = magnitud de interés
Mmax= magnitud máxima
Mmin = magnitud mínima.

4.3.1 Determinación del valor de b

El valor de b puede ser determinar por dos métodos: el método de los
mínimos cuadrados o e método de máxima verosimilitud. El método de los
mínimos cuadrados ajusta los valores de la muestra de sismos a una recta
en función de la cantidad de datos, por lo que si existen escasos datos con
valores grandes la influencia de estos en la suma de los cuadrados será
desproporcionada.

El método de los mínimos cuadrados debe ser evitado en los casos en que
uno desee ajustar datos de una distribución acumulativa, como en el caso
de la relación de Gutenberg-Ritcher, debido a que lleva a estimaciones
inestables de los valores de la muestra (Vere-Jones, Yosihiko 2003).

El método de máxima verosimilitud ajusta la recta al valor medio más

probable de la muestra, en función a valores de magnitud mínima y
máxima. Además, debido a que su formulación está en base a una función
de densidad de probabilidades su uso representa mejor la muestra que el
método de mínimos cuadrados. El valor de b por el método de máxima
verosimilitud (Utsu 2003) es igual a:
log e
b= (4.1)
( M  M min)
Donde M es el valor medio de la muestra de magnitudes.

78
4.3.2 Determinación de la magnitud mínima y razón media anual
La magnitud mínima de la fuente sísmica corresponde al punto de inflexión
de la curva de recurrencia sísmica y se obtiene calculando el punto de
máxima curvatura (Wiemer and Katsumata,1999; Wiemer and Wyss, 2000).
El punto de máxima curvatura es definido como el punto que corresponde al
valor máximo del número de sismos no acumulativos versus la magnitud.
La figura 4.6 muestra el valor correspondiente al punto de máxima
curvatura encontrado en la fuente 3, donde la magnitud mínima hallada fue
de 3.40.

Figura 4.5.- Curva de recurrencia sísmica de la fuente F3 con la ubicación del

punto de inflexión que representa a la magnitud mínima (Mmin).

79
Figura 4.6.- Número de sismos para cada magnitud encontrados en la fuente F3.

La correcta estimación de la magnitud mínima es de mucha importancia en

los análisis de sismicidad y peligro sísmico debido a que esta relacionada
directamente con el número de sismos esperados en la fuente. La magnitud
mínima influencia significativamente en la determinación de las
aceleraciones y ordenadas espectrales máximas esperadas.

La razón media anual () de excedencia es definida como el número de

sismos anuales que pueden ocurrir mayores o iguales a la magnitud mínima
y es igual a =exp (-mmin).
Donde:
 = ln10*a;  = ln10*b;
Mmin = Magnitud mínima.

La recta de ajuste que representa la recurrencia sísmica para la fuente F3

se encuentra graficada en la figura 4.7. En la figura se muestra el valor de
“b” obtenido y su dispersión basada en el número de eventos y las
magnitudes de la fuente. Además se muestra el valor de “a” para obtener el
número de sismos mayores o iguales a la magnitud mínima.

La figura 4.8(a) muestra la distribución acumulativa que tendrían los sismos

basados en la recta de recurrencia hallada con el método de máxima

80
verosimilitud y los sismos acumulativos propios de la fuente F3. En esta se
puede observar que la distribución acumulativa basada en la recta de ajuste
encontrada calcula bastante bien el número de sismo para la fuente.

La figura 4.8(b) muestra la distribución acumulativa que tendrían los sismos

basados en la recta de recurrencia hallada con el método de mínimos
cuadrados y los sismos acumulativos propios de la fuente F3. Las
variaciones son evidentes.

Figura 4.7.- Recta de ajuste (Nm=a-bm) que representará la sismicidad aportada

por la fuente F3 al peligro sísmico del Perú.

81
 500
Nm-Mw
450

400
Nm-Mw(Ajuste)
Num. Acumulativo Sismos

350 Beta=1.027
Razón=12.75
300 Mmin=3.4

250

200

150

100

50

0
3.00 4.00 5.00 6.00 7.00 8.00 9.00
Mw

Figura 4.8(a).- Número de sismos mayores o iguales a diferentes magnitudes

obtenidos del ajuste de la figura 4.7 vs. el número de sismos de la fuente F3.

800
Nm-Mw
700

600
Num. Acumulativo Sismos

Nm-Mw(Ajuste)
Beta=1.51
500
Razón=20.59
Mmin=3.4
400

300

200

100

0
3.00 4.00 5.00 6.00 7.00 8.00 9.00
Mw
Figura 4.8(b).- Número de sismos mayores o iguales a diferentes magnitudes
utilizando el método de mínimos cuadrados vs. el número de sismos de la fuente
F3.

4.3.3 Determinación de la magnitud máxima

La magnitud máxima (Mmax) es definida como el evento más grande que
pueda ocurrir en la fuente sísmica (McGuire 1976). El criterio para
determinarla esta basado en que el evento más grande ocurrido en el
pasado en la fuente sísmica también puede ocurrir en el futuro.

82
La determinación de la magnitud máxima se obtuvo recopilando información
sísmica histórica e instrumental ocurrida en las diferentes fuentes sísmicas
utilizadas en este trabajo.

La información sísmica histórica estuvo basada en el trabajo de Dorbath et

al. (1990). La información instrumental se obtuvo del catálogo del proyecto
SISRA (IGP 2001). Las magnitudes máximas utilizadas en este trabajo se
encuentran en la tabla 4.2.

4.4 Profundidades representativas de las fuentes sismogénicas

Debido a que la caracterización de las fuentes fue hecha de manera
volumétrica es necesario determinar una profundidad que represente a los
eventos dentro de la fuente. Las profundidades representativas fueron
obtenidas de histogramas de la frecuencia del número de sismos para cada
profundidad de los eventos presentados en el catálogo actualizado del
proyecto SISRA. Los histogramas se encuentran en el anexo 1.

Los parámetros sismológicos encontrados para las 20 fuentes sísmicas que

servirán para el análisis del peligro sísmico se resumen en la tabla 4.2
detallada a continuación.

83
 Tabla 4.2. Parámetros sismológicos de las 20 fuentes.

FUENTE Mmin Mmax BETA RAZON PROFUNDIDAD (Km)

(Mw) (Mw) () SISMOS ()
1 4,50 7,33 1,248 0,939 35,00
2 4,00 8,00 1,142 4,050 30,00
3 3,40 8,70 1,027 12,75 30,00-60,00
4 3,60 8,80 1,098 6,832 30,00-60,00
5 3,60 8,80 0,928 4,797 30,00-60,00
6 3,60 7,16 1,091 2,215 35,00
7 4,30 7,90 1,950 0,288 30,00-50,00
8 4,30 6,80 2,190 0,552 30,00
9 3,80 7,20 1,087 1,143 45,00
10 3,80 7,33 1,032 2,087 35,00
11 4,50 8,15 1,881 3,940 35,00-60,00
12 4,10 6,90 1,395 0,808 40,00
13 3,80 7,00 1,637 1,259 90,00
14 3,80 7,10 1,506 2,435 90,00
15 3,60 8,20 1,289 5,319 100,00
16 3,60 7,90 1,059 6,249 120,00
17 4,50 8,40 1,628 10,48 90,00-110,00-160,00
18 3,30 7,33 1,038 9,850 110,00-170,00
19 3,40 7,00 1,013 8,240 120,00-180,00
20 4,30 7,10 0,764 0,987 610,00

4.5 Leyes de atenuación del movimiento del suelo

La estimación del movimiento del suelo depende de la magnitud y la
distancia desde el origen del sismo hasta el sitio. Esta estimación se obtiene
de las relaciones o leyes de atenuación. La elección de la ley de atenuación
constituye un paso muy importante en la determinación del peligro sísmico
del lugar.

Actualmente es posible desarrollar modelos de atenuación del movimiento

del suelo para sismos con características regionales y tectónicas diferentes.
Por ejemplo sismos que se desarrollan en la corteza en regiones
tectónicamente activas, sismos en la zona de subducción de interfase y
sismos que ocurren en la placa de subducción descendente.

84
El Perú, caracterizado por movimientos sísmicos que ocurren debido al
proceso de subducción de la placa de Nazca, cuenta únicamente con 14
registros corregidos propios de esta característica tectónica.

Debido a esta limitada cantidad de registros de aceleraciones, es difícil

desarrollar una relación de atenuación para estimar parámetros como la
aceleración del suelo y ordenadas espectrales para los diferentes tipos de
sismos en el Perú (sismos de interfase e intraplaca y sismos de corteza).
Además, debido a que los datos estarían muy esparcidos, se presentarían
valores de incertidumbre muy grandes en la extrapolación utilizada para
desarrollar el modelo de atenuación.

En el Perú, sin embargo, existen leyes de atenuación para el movimiento del

suelo propuestas por Casaverde (1980), Huaco (1980) y Ruiz (1999). Estas
relaciones lamentablemente no pudieron utilizarse en este trabajo debido a
los escasos datos utilizados en su estimación y en la medida del tamaño del
sismo en que están basadas (generalmente mb o Ms).

En la actualidad es posible utilizar sismos que hayan ocurrido en regiones

diferentes a las de análisis, siempre y cuando obedezcan al mismo proceso
tectónico de generación de sismos.

Youngs et al. (1997) desarrollaron leyes de atenuación para zonas de

subducción de sismos de interfase e intraplaca usando datos de sismos
registrados en Alaska, Chile, Cascadia, Japon, Mexico, Perú (14 registros) y
las Islas Salomón para distancias entre 10 y 500 km, teniendo en cuenta
además las características del sitio.

Youngs et al. (1997) definieron las características del sitio clasificándolas en

tres grupos: roca, suelo duro poco profundo y suelo profundo. Se consideró
eventos en roca a todos aquellos con velocidades de corte cercanos a los
750 m/seg, eventos en suelo profundo aquellos con distancias a la roca
mayores de 20 m y con velocidades de corte entre 180 y 360 m/seg y
eventos en suelo poco profundo aquellos donde la profundidad del suelo es
menor a los 20 m.

85
Youngs utilizo la magnitud momento Mw (Hanks and Kanamori 1979) para
la medida del evento. La localización epicentral, profundidad, magnitud y
mecanismo focal fueron obtenidos de publicaciones especiales o del
Harvard centroid moment tensor solutions. Youngs mostró que las
aceleraciones máximas obtenidas de sismos de subducción atenúan más
lentamente que aquellas producidas por sismos de corteza, y que los sismos
de intraplaca producen aceleraciones mayores que los de interfase para la
misma magnitud y distancia. En la figura 4.9 se muestran algunas curvas
de atenuación estimadas por Youngs para eventos con distintas
magnitudes.

Figura 4.9.- Valores de aceleración estimados por Youngs para sismos de

interfase, los números en paréntesis corresponden a la profundidad promedio de los
eventos ( Youngs et al. 1997).

86
Youngs propuso la siguiente ley de atenuación en roca:

ln(A) =0.2418+1.414M-2.552ln(R+1.7818 e0.554M)+0.00607H+0.3846Zt (4.2)

Desviación estándar del ln(A) = 1.45-0.1M (4.3)
Donde:
A = aceleración del suelo (g).
M = magnitud momento (Mw).
R = distancia más cercana a la ruptura (km).
H = profundidad (km).
Zt = 0 para interfase, 1 para intraplaca.

Las relaciones de atenuación para las aceleraciones en roca propuesta por

Youngs para diferentes magnitudes se encuentran graficadas en la figura
4.10. En la misma figura se muestran algunos sismos registrados en el
Perú. Los sismos peruanos, representados por triángulos, corresponden a
eventos con Mw entre 6,0 y 6,9 y H entre 50 y 58 km.

La tabla 4.3 muestra las aceleraciones obtenidas con la relación de

atenuación propuesta por Youngs y la propuesta por Casaverde y las
aceleraciones de las dos componentes horizontales registradas en diferentes
sismos peruanos en la estación del Instituto Geofísico del Perú (IGP). Se
eligieron únicamente los eventos registrados en la estación del Instituto
Geofísico del Perú (Parque de la reserva) debido a que las velocidades de
corte de su suelo están entre los 525 m/seg y 700 m/seg (Repetto et al.
1980), correspondientes aproximadamente con la clasificación hecha por
Youngs en roca. El valor promedio (P.50) corresponde al valor obtenido
directamente de la ley de atenuación, mientras que el valor del 84% perc.
es el valor del promedio más la desviación estándar (P.84).

En la tabla 4.3 y en la figura 4.10 se observa que los valores de aceleración

registrados en el IGP se encuentran dentro del rango de valores obtenidos
con la relación de atenuación de Youngs. Esto permite dar validez a la ley
de atenuación elegida.

87
 YOUNGS-INTERFASE - ROCA - PROF=50 km

1.000
Mw=8

Mw=7

Mw=6

Mw=5

Mw 6-6.9

Mw 7-7.9
Aceleración máxima (g)

0.100

0.010
10 100 1000
Distancia (Km)

Figura 4.10.- Curvas de atenuación propuestas por Youngs vs aceleraciones

registradas en Lima para eventos Mw entre 6.0 y 7.9.

88
Tabla 4.3. Aceleraciones registradas en diferentes eventos por el IGP, y estimadas
usando la ley de Youngs et al. (1997) y Casaverde (1980)
[a(cm/seg2)=68.7*(exp(0.8*Ms))*(R+25)-1].

Evento Mw Ms R H IGP Youngs Youngs Casaverde Casaverde

(km) (km) (g) (P.50) (P.84) (P.50) (P.84)

(g)
(g) (g) (g)

31/01/51 6,18 6,0 116 50 0,046 0,022 0,043 0,060 0,108

0,061

17/10/66 8,10 7,5 165 24 0,180 0,055 0,094 0,148 0,2669

0,269

31/05/70 7,90 7,8 260 56 0,097 0,029 0,050 0,126 0,226

0,104

29/11/71 5,81 5,5 138 54 0,087 0,011 0,023 0,034 0,062

0,053

5/01/74 6,55 6,6 123 98 0,072 0,055 0,106 0,093 0,167

0,066

3/10/74 8,10 7,6 74 13 0,192 0,119 0,205 0,309 0,555

0,179

9/11/74 7,00 7,2 75 15 0,069 0,065 0,121 0,222 0,3989

0,046

Para representar la atenuación de la aceleración de los sismos de corteza se

utilizó la relación propuesta por Sadigh et al. (1997). Esta relación esta
basada principalmente en sismos de la Costa Oeste de los Estados Unidos y
en datos obtenidos de los sismos de Gazli (Rusia, 1976) y Tabas (Irán,
1978). La medida del tamaño del evento fue caracterizada por la magnitud
momento, Mw y la distancia definida como la menor distancia a la ruptura.

La ley de atenuación propuesta por Sadigh en roca es la siguiente:

ln(A)=-0.7488+1.2M-2.52ln(R + exp(1.55579 + 0.3M)) (4.4)
Desviación estándar del ln (A) = 1.39-0.14M; 0.38 para M  7.21 (4.5)

89
 Donde:
A = aceleración del suelo en g.
M = magnitud momento (Mw).
R = distancia más cercana a la ruptura en km.

Las relaciones de atenuación para las aceleraciones en roca propuesta por

Sadigh et al. 1997 para diferentes magnitudes se encuentran graficadas en
la figura 4.11.

1.000
Mw=8
Mw=7
Mw=6
Mw=5

0.100
Aceleración Máxima (g)

0.010

0.001
10 100 1000

Distancia (km)

Figura 4.11.- Curvas de atenuación propuestas por Sadigh k., Chang C., Egan J.,
Makdisi F., Youngs R. (1997) para sismos de corteza.

90
4.6 Leyes de atenuación para ordenadas espectrales
Para estimar la respuesta estructural que un sismo de una magnitud y
distancia específica pueda inducir a diferentes edificaciones se utilizaron
modelos de atenuación para ordenadas espectrales en función de periodos
de vibración estructural específicos. Los modelos de atenuación usados
corresponden a los propuestos por Youngs et al. (1997) para sismos de
subducción de interfase e intraplaca y el modelo de Sadigh et al. (1997)
para los sismos de corteza.

La metodología utilizada por ambos es similar y fue descrita en el acápite

2.3.3 resumida a continuación por conveniencia. El proceso consiste en
procesar los registros de aceleraciones de los eventos sísmicos más fuertes
contenidos en la base de datos y obtener relaciones de amplificación del
movimiento del suelo con respecto a la respuesta estructural (Sa/Amax). La
relación de atenuación de ordenadas espectrales es luego obtenida de
combinar los coeficientes de amplificación espectral con la relación de
atenuación apropiada de la aceleración del suelo.

Las relaciones de atenuación propuestas por ambos investigadores

corresponden a un amortiguamiento de 5% y a periodos estructurales de
0,07; 0,075; 0,1; 0,2; 0,3; 0,4; 0,5; 0,75; 1,0; 1,5; 2,0 y 3,0 segundos.

La ecuación de atenuación para ordenadas espectrales propuesta por

Youngs en roca es:
ln(Sa)=0.2418+1.414M+C1+C2(10-M)^3+C3ln(R+1.7818 e0.554M)+
0.00607H+0.3846Zt (4.6)
Desviación estándar de ln(Sa) = C4+C5M (4.7)
Donde:
Sa = ordenada espectral (g); M = magnitud momento (Mw).
R = distancia más cercana a la ruptura (km).
H = profundidad (km); Zt = 0 para interfase, 1 para intraplaca.
Los coeficientes se detallan a continuación:

91
 Tabla 4.4. Coeficientes para la ley de atenuación de ordenadas espectrales
en roca propuesta por Youngs et al. (1997).
Periodo C1 C2 C3 C4 C5
(seg)
0,075 1,275 0,0000 -2,707 1,45 -0,1
0,100 1,188 -0,0011 -2,655 1,45 -0,1
0,200 0,722 -0,0027 -2,528 1,45 -0,1
0,300 0,246 -0,0036 -2,454 1,45 -0,1
0,400 -0,115 -0,0043 -2,401 1,45 -0,1
0,500 -0,400 -0,0048 -2,360 1,45 -0,1
0,750 -1,149 -0,0057 -2,286 1,45 -0,1
1,000 -1,736 -0,0064 -2,234 1,45 -0,1
1,500 -2,634 -0,0073 -2,160 1,50 -0,1
2,000 -3,328 -0,0080 -2,107 1,55 -0,1
3,000 -4,511 -0,0089 -2,033 1,65 -0,1

La ecuación de atenuación para ordenadas espectrales propuesta por

Sadigh en roca es:
ln(Sa)=C1+C2M+C3(8.5M)^2.5+C4ln(R+exp(C5+C6M))+C7ln(R+2) (4.8)
Donde:
Sa = aceleración del suelo (g).
M = magnitud momento (Mw).
R = distancia más cercana a la ruptura (km).

Los coeficientes de la ley de atenuación se describen en la tabla 4.5 y 4.6.

Tabla 4.5. Coeficientes para la ley de atenuación de ordenadas espectrales en roca

propuestos por Sadigh et al. (1997), para eventos con M  6.5.
Periodo C1 C2 C3 C4 C5 C6 C7
(seg)
0,07 0,1320 1,2 0,0072 -2,5536 1,55579 0,3 -0,0989
0,10 0,3300 1,2 0,0072 -2,5776 1,55579 0,3 -0,0492
0,20 0,1836 1,2 -0,0048 -2,4960 1,55579 0,3 0,0000
0,30 -0,0684 1,2 -0,0204 -2,4336 1,55579 0,3 0,0000
0,40 -0,3576 1,2 -0,0336 -2,3880 1,55579 0,3 0,0000
0,50 -0,7056 1,2 -0,0480 -2,3340 1,55579 0,3 0,0000
0,75 -1,4496 1,2 -0,0600 -2,2380 1,55579 0,3 0,0000
1,00 -2,0460 1,2 -0,0660 -2,1600 1,55579 0,3 0,0000
1,50 -2,8880 1,2 -0,0780 -2,0700 1,55579 0,3 0,0000
2,00 -3,5340 1,2 -0,0840 -2,0040 1,55579 0,3 0,0000
3,00 -4,4400 1,2 -0,0960 -1,9320 1,55579 0,3 0,0000

92
Tabla 4.6. Coeficientes para la ley de atenuación de ordenadas espectrales en roca
propuestos por Sadigh et al. (1997), para eventos con M > 6.5.

Periodo C1 C2 C3 C4 C5 C6 C7
(seg)
0,07 -0,6480 1,32 0,0072 -2,5536 -0,5814 0,6288 -0,0984
0,10 -0,4500 1,32 0,0072 -2,5776 -0,5814 0,6288 -0,0492
0,20 -0,5964 1,32 -0,0048 -2,4960 -0,5814 0,6288 0,0000
0,30 -0,8484 1,32 -0,0204 -2,4330 -0,5814 0,6288 0,0000
0,40 -1,1376 1,32 -0,0336 -2,3880 -0,5814 0,6288 0,0000
0,50 -1,4856 1,32 -0,0480 -2,3340 -0,5814 0,6288 0,0000
0,75 -2,2296 1,32 -0,0600 -2,2380 -0,5814 0,6288 0,0000
1,00 -2,8260 1,32 -0,0660 -2,1600 -0,5814 0,6288 0,0000
1,50 -3,6684 1,32 -0,0780 -2,0700 -0,5814 0,6288 0,0000
2,00 -4,3140 1,32 -0,0840 -2,0040 -0,5814 0,6288 0,0000
3,00 -5,2200 1,32 -0,0960 -1,9320 -0,5814 0,6288 0,0000

Las desviaciones estándar del logaritmo natural de las ordenadas

espectrales estimadas para cada periodo son:

Tabla 4.7
Periodo Sigma[ln(y)]
0,07 1,40-0,14M; 0,39 para M  7,21
0,10 1,41-0,14M; 0,40 para M  7,21
0,20 1,43-0,14M; 0,42 para M  7,21
0,30 1,45-0,14M; 0,44 para M  7,21
0,40 1,48-0,14M; 0,47 para M  7,21
0,50 1,50-0,14M; 0,49 para M  7,21
0,75 1,52-0,14M; 0,51 para M  7,21
 1,00 1,53-0,14M; 0,52 para M  7,21

Con la finalidad de conocer la respuesta estructural que estima la relación

propuesta por Youngs y la respuesta estructural para un sismo peruano, se
construyeron espectros basados en las características de sismos ocurridos
en el pasado. Los eventos elegidos corresponden a los sismos de 1974
ocurridos en Lima, debido a que están muy bien documentados. Los
parámetros de estos eventos se encuentran detallados en la tabla 4.3.

Las figuras 4.12 y 4.13 muestran la respuesta estructural basada en la

ecuación 4.6 para el valor medio P(.50) y el valor medio más la desviación
estándar P(.84) correspondientes a los sismos de 1974. Asimismo, se

93
muestran los espectros de respuesta para las dos componentes horizontales
registradas de los sismos.

En las graficas se observa que los espectros propuestos por Youngs

aproximan bastante bien la respuesta estructural que tuvieron las
edificaciones durante los sismos de 1974. En ambos casos los espectros que
incluyen la variabilidad de la respuesta estructural (P(.84)) se presentan
como un límite superior a los espectros de respuesta de los sismos de 1974.

Se hace necesario remarcar que el espectro que se utiliza para el diseño es

el espectro suavizado del espectro de respuesta, por lo que la comparación
ante los sismos de 1974 debe hacerse con este criterio.

Esta comparación da validez para poder calcular ahora espectros que

tengan en cuenta la posibilidad de ocurrencia de sismos con diferentes
magnitudes y distancias e incluir en su determinación la incertidumbre en la
respuesta estructural. Los espectros calculados se encuentran en el capitulo
VI.

94
 Respuesta Espectral Sismo 3/10/74
Mw=8.1, H=13km, Rrup=74km
IGP-C7410N08
0.70

IGP-C7410N82
0.60

0.50 YOUNGS - P(.50)

Mw=8.1,Rrup=74km
Sa (g)

0.40 YOUNGS -P(.84)

Mw=8.1,Rrup=74km
0.30

0.20

0.10

0.00
0.0 0.5 1.0 1.5 2.0 2.5 3.0

Periodo (seg.)

Figura 4.12.- Espectros de respuesta para el sismo de Lima del 3/10/74 vs. el
espectro obtenido con el modelo de Youngs para el valor medio P(.50) y el valor
medio más la desviación estándar P(.84).

Respuesta Espectral Sismo 9/11/74

Mw=7.0, H=15km, Rrup=75km
0.30
IGP-C7411N08

0.25
IGP-C7411N82

0.20 YOUNGS - P(.50)

Sa (g)

Mw=7.0,Rrup=75km
0.15 YOUNGS -P(.84)
Mw=7.0,Rrup=75km
0.10

0.05

0.00
0.0 0.5 1.0 1.5 2.0 2.5 3.0

Periodo (seg.)

Figura 4.13.- Espectros de respuesta para el sismo de Lima del 9/11/74 vs. el
espectro obtenido con el modelo de Youngs para el valor medio P(.50) y el valor
medio más la desviación estándar P(.84).

95
 CAPITULO V

DISTRIBUCION PROBABILISTICA DE
ACELERACIONES EN EL PERU

5.1 Introducción

La amplitud de las aceleraciones máximas esperadas en un lugar

específico esta en función de las fuentes que representan la sismicidad de la
zona y la ley de atenuación usada para representar el movimiento en cada
fuente.

Las fuentes que definen la sismicidad en el Perú han ido variando conforme
se obtenían más datos. Es debido a esta constante variación en la
distribución espacial de sismos que ha sido necesario revisar y actualizar las
últimas fuentes definidas para el Perú.

Las leyes de atenuación, por otro lado, son actualizadas cada cierto periodo
de tiempo utilizando métodos cada vez más sofisticados e incluyendo
nuevas variables en su estimación. En este trabajo se revisaron distintas
leyes de atenuación y se eligieron modelos diferentes a los utilizados en
estudios anteriores de peligro en el Perú. Entre las características más
importantes de la ley de atenuación usada esta la capacidad de distinguir
eventos de subducción de interfase e intraplaca.

La información de las fuentes y leyes de atenuación puede ser combinada

con la ayuda de programas de computadora. Hoy en día existen diversos
programas para el cálculo de las aceleraciones. Sin embargo, en este
trabajo se ha preferido desarrollar un nuevo programa que calcule la
aceleración del suelo y la respuesta estructural de una manera didáctica y
simple.

96
En este capitulo se presenta trabajos anteriores de peligro sísmico y se
recalcula la distribución de aceleraciones en el Perú. Se describe además las
características del programa de cómputo desarrollado.

5.2 Trabajos anteriores

El peligro sísmico para el Perú fue estudiado principalmente en los trabajos
de Casaverde y Vargas (1980), Sharma y Candia-Gallegos (1991) y Castillo
y Alva (1993). Cada uno estudio características diferentes en cuanto a
fuentes sísmicas, leyes de atenuación y probabilidades de excedencia.

El trabajo de Casaverde y Vargas (1980) fue uno de los primeros estudios

de análisis de peligro sísmico desarrollado en el país. Ellos desarrollaron una
ley de atenuación basada en 10 registros sísmicos obtenidos en la ciudad de
Lima y plantearon las primeras fuentes sísmicas y calcularon la distribución
de aceleraciones. Con esto contribuyeron a la zonificación sísmica del País.
La distribución de aceleraciones, sin embargo, fue hecha para diferentes
probabilidades de excedencia. Debido a esto sus valores no pueden
compararse con los obtenidos en este trabajo.

Sharma y Candia-Gallegos (1991) obtuvieron valores de aceleración para

un periodo de retorno de 475 años en roca para el Perú. Su trabajo se basó
en cinco fuentes sísmicas y en la ley de atenuación de Krinitzsky et al.
(1988). Los valores de aceleración máxima fueron para las ciudades de
Arequipa (0,52g), Lima (0,44g), Cuzco (0,43g), Huancayo (0,39g), Chiclayo
(0,36g) y Piura (0.43g).

Castillo y Alva (1993) obtuvieron valores de aceleraciones para un 10% de

probabilidad de excedencia en suelo firme para periodos de exposición de
50 y 100 años utilizando una base de datos correspondiente al periodo
entre 1963 y 1990, con 3892 eventos sísmicos. En su trabajo definieron 20
fuentes sismogénicas y utilizaron la ley de atenuación propuesta por
Casaverde y Vargas (1980) para representar el movimiento del suelo.

Los valores más altos obtenidos por Castillo y Alva corresponden a 0.50g y
0.60g en 50 y 100 años respectivamente para las ciudades de Tumbes,

97
Piura, Ica, Tacna y el Norte de Chile. Los valores más bajos fueron desde
0,06g y 0,08g ubicados en el departamento de Loreto. Asimismo la ciudad
de Lima presentó aceleraciones entre 0,44g y 0,42g para un periodo de
exposición de 50 años.

5.3 Distribución de aceleraciones en el Perú

El cálculo de la aceleración debe ser hecha para algún grado de peligro o
periodo de retorno. La selección del grado de peligro o periodo de retorno
depende de consideraciones sociales, económicas y políticas envueltas en
su estimación. Para muchos análisis de peligro sísmico a nivel mundial, el
movimiento sísmico de diseño para estructuras comunes ha sido establecido
para un 10% de probabilidad de ser excedido en un periodo de exposición
de 50 años (periodo de retorno de 475 años).

La distribución de aceleraciones realizada en este trabajo corresponde a

periodos de exposición de 50 y 100 años y un 10 % de probabilidad de
excedencia. Los valores hallados se grafican en mapas, con curvas de igual
aceleración.

El calculo de aceleraciones en este trabajo esta basado en 20 fuentes

sísmicas propuestas por Castillo y Alva (1993). Sin embargo, debido a la
variación en la distribución de sismos principalmente en la costa central se
cambio la geometría de las fuentes F3, F4 y F8. Los parámetros utilizados
en el análisis de peligro sísmico de estas fuentes se muestran en la Tabla
5.1.

Tabla 5.1.- Parámetros sismológicos de las fuentes actualizadas en este trabajo.

FUENTE Mmin Mmax BETA RAZON PROFUNDIDAD (Km)

(Mw) (Mw) () SISMOS
3 3,40 8,70 1,027 12,750 30,00-60,00
4 3,60 8,80 1,098 6,832 30,00-60,00
8 4,30 6,80 2,190 0,552 30,00

98
El movimiento del suelo que se produce debido a eventos que ocurren en la
zona de subducción de interfase e intraplaca fue caracterizado utilizando la
ley de atenuación de Youngs et al. (1997); mientras que para representar el
movimiento del suelo que se produce debido a sismos de corteza se uso la
ley de atenuación propuesta por Sadigh et al. (1997).

El calculo de aceleraciones se hizo para una cuadricula que se construyó

dividiendo los ejes correspondientes a las longitudes en 0,1° y los ejes
correspondientes a las latitudes en 0,5°; desde la coordenada geográfica
0,0° hasta la coordenada -19,0° de latitud sur, y desde la coordenada -
82,0° hasta la coordenada -68,0° de longitud oeste.

Los valores obtenidos para cada punto de la cuadricula fueron luego

interpolados a las aceleraciones mostradas en las figuras 5.1 y 5.2. Los
valores de aceleración como se observa van disminuyendo de Oeste a Este.
Los valores de aceleración más altos se distribuyen paralelos a la línea de
costa. Esto se debe a que la mayor cantidad de sismos se producen en la
zona donde empieza la subducción, en la fosa Peruano – Chilena. En la zona
de Arequipa, entre los 15,5° y los 16,5° de latitud sur aproximadamente, se
observa un cambio en la pendiente de las curvas, que coincide con el
cambio en la subducción de la placa de Nazca debido a una contorsión de
esta.

Los resultados en este trabajo fueron comparados con el trabajo de Castillo

y Alva. Se encontraron diferencias en la amplitud de las aceleraciones,
debido a que en este trabajo se utilizo una base de datos de sismicidad
actualizada en 10 años, diferente medida del tamaño del evento (Castillo y
Alva utilizaron el máximo entre Mb y Ms) y diferente ley de atenuación.

La principal variación en cuanto a amplitud de aceleraciones se ha

encontrado en la zona norte del País, en los departamentos de Piura y
Tumbes, donde los valores han disminuido en un 32% debido al poco
incremento en la sismicidad de la zona. Esto se evidencia en la disminución
de la razón anual de sismos presentada en la fuente F2 que disminuyó de
4,58 sismos de magnitud Mw por año presentada para el periodo entre

99
1963 y 1990 a 4,05 sismos de magnitud Mw por año presentada para el
periodo entre 1963 y 2000. La tabla 5.2 muestra los valores obtenidos en
este trabajo y los valores obtenidos por Castillo y Alva para algunas
ciudades del País.

Tabla 5.2.- Aceleraciones presentadas en diferentes ciudades del Perú.

Lugar Este Trabajo Alva Castillo (1993)

(g) (g)
Lima (-77,00;-12,00) 0,40 0,43
Ica (-75,48;-14,02) 0,41 0,42
Arequipa (-71,60;-16,30) 0,36 0,38
Tacna (-70,30;-18,10) 0,39 0,41
Piura (-80,69;-6,00) 0,36 0,47

100
 -82 -80 -78 -76 -74 -72 -70 -68
+1 +1

0 0

ECUADOR COLOMBIA
-2 -2
0.
36
g

-4 0.
22 -4
g

0.
0.

04
0.
20

g
06
0.
g

g
08
-6 -6
0.

g
18 g
0.
20
g

0.
10
g
0.

0.
12
36

0.

g
g

14
g
-8 -8
BRASIL
0.
16
0.

g
18

0.45 g 0.
g

14
0.

g
20
0.

g
22
0.

0.40 g
g
0.

24

-10 -10
0..302g 0 .3.38 0g
0 .3 4g 0 .4

26
0 .3

g
0

28 g

g
g

0.34 g
OC

0
6g g

0.28 g
EA
NO

-12 -12
0.

0.22 g
12

0.
g
0.
PA

0 16
14

0. .18 g
g

0.16 g 0. 20 gg
CI

0 2
0 .2 2g
0 .26 4g
FI

0. .28g g
CO

30
0.10 g 0.
32 g
g
-14 0 0.1 -14
0. .14 2g
0 . 6g

0.04 g 1 g
0 0.1 6g
34
0. .38 g
0 0 2g

0. .20 8g
3 g
0. .4 44gg

g
0 44
4

22 g
00..

0.02 g 0g
0 .2
0. .264g
28 g
g
-16 BOLIVIA -16
PONTIFICIA UNIVERSIDAD CATOLICA DEL PERU
0 .30 2g

ESCUELA DE GRADUADOS
0.3

MAPA DE ORDENADAS ESPECTRALES

g

Periodo estructural (Tn) : 0.00 seg

0.34

-18 Probabilidad de excedencia : 10% -18

0.

g
36

Periodo de exposición : 50 años 0.

38
g

0 . 0 . 40 g
MANUEL MONROY, ANA BOLAÑOS - 2004 42 g
g
km
CHILE
-20 100 200 300 400 -20
-82 -80 -78 -76 -74 -72 -70 -68

Figura 5.1.- Distribución de aceleraciones en el Perú para un periodo de exposición

de 50 años y 10% de probabilidad de excedencia.

101
 -82 -80 -78 -76 -74 -72 -70 -68
+1 +1

0 0

ECUADOR COLOMBIA
-2 -2
0. 38 g
0. .40 2g
36 g
0 .4 4g

g
0 .4
0

-4 0. -4
0.
26
08
g
g

0.
0.

10
24

g
0.

0.
2
0g

12
g
-6 -6
0. .24
22 g
0

0.
g

1
0.

4g

0.1
44

6
g

g
0.1
-8 8g -8
BRASIL

0.52 g
0.1
8g
0.48 g
-10 -10

0.16g
0 .2
0g

0.42 g
0.2 2g
OC

0.36 g
EA

0.

0.
2
4g

0. 26g
NO

-12 0.30 g
28
0.30 g
-12
0 .3 g
PA

2g
0.34g
0.24 g
CI

0.36 g
F

0.38g
IC

0.18 g 0.4
O

0g
-14 -14
0. 44 6g
0. .448gg

0.12 g
42 g
0 . 50 g
0 . 2

g
0 .5

0.0.2.246g 0.30
0

20 2g g .2 g 0 . 3 4g
0 .2 0

0.08 g
0

g 8g 3 2g

-16 -16
PONTIFICIA UNIVERSIDAD CATOLICA DEL PERU
0

BOLIVIA
0.368g 0.40g
.
0 .3

ESCUELA DE GRADUADOS
g

MAPA DE ORDENADAS ESPECTRALES

Periodo estructural (Tn) : 0.00 seg
0 .4

-18 Probabilidda de excedencia : 10% -18

2 g 4g g
0. 0. .48 0g

Periodo de exposición : 100 años

4 46 g
0 .5
0

MANUEL MONROY, ANA BOLAÑOS - 2004

km

100 200 300 400

CHILE
-20 -20
-82 -80 -78 -76 -74 -72 -70 -68

Figura 5.2.- Distribución de aceleraciones en el Perú para un periodo de exposición

de 100 años y 10% de probabilidad de excedencia.

102
5.4 Programa de cómputo para el cálculo de aceleraciones
Actualmente existen distintos algoritmos de computadora (McGuire 1976,
1978; Bender and Perkins 1987), que incluyen rutinas para el cálculo del
peligro sísmico en las que las fuentes sísmicas pueden caracterizarse como
fuentes lineales o fuentes volumétricas.

En la base de estos algoritmos y en la metodología planteada en los

capítulos II, III y IV se ha desarrollado un programa para el cálculo del
peligro sísmico. El programa se llama MRiesgo y sirve para calcular tanto la
aceleración del suelo como ordenadas espectrales (Ver Anexo II). El
programa incluye una interfase gráfica para el almacenamiento de datos y
presentación de resultados. El programa MRiesgo fue escrito en el entorno
Matlab debido a la gran librería de funciones matemáticas y estadísticas que
presenta.

El programa MRiesgo grafica los resultados en función de probabilidades de

excedencia versus un parámetro del movimiento del suelo (aceleración
máxima esperada o respuesta estructural) para un solo sitio o distintos
sitios. El programa crea además un archivo de texto con los datos de
entrada proporcionados por el usuario y los resultados más importantes
obtenidos como el número de sismos anuales aportados por cada fuente a
cada aceleración deseada y la probabilidad de excedencia de cada
aceleración para distintos periodos de exposición o vida útil de las
estructuras.

La sismicidad es caracterizada en el programa por medio de fuentes

volumétricas. El usuario debe introducir las coordenadas geográficas de los
bordes de cada fuente en el sentido de las agujas del reloj. Luego de
definida la geometría de cada fuente, el programa muestra un cuadro de
diálogo para ingresar los principales parámetros sismológicos de la fuente.
Los parámetros son: un coeficiente entre 0 y 1, magnitud mínima, magnitud
máxima, beta, razón de sismos, profundidad y tipo de atenuación. El
coeficiente se utiliza para modificar la sismicidad de la fuente. El coeficiente
será cero si no quiere incluirse la sismicidad de alguna fuente y uno para
incluirla.

103