Actualización de los Parámetros

Sismológicos en la Evaluación del Peligro

Sísmico en el Perú

Dr. Jorge Alva Hurtadp

Con el fin de proveer una mejor información,

I Resumen

El presente trabajo presenta los

resultados de la evaluación de la
actividad sísmica en cada una de las
fuentes sismogénicas propuestas para
el Perú, habiéndose calculado la
relación de recurrencia de la actividad
sísmica ocurrida en el país en el
periodo comprendido entre los años
1963 y 2003.

Se ha realizado una revisión de la

actividad sísmica del pasado,
reafirmando la definición de las fuentes
sismogénicas, teniendo en cuenta las
características tectónicas de la región,
donde la probabilidad de ocurrencia de
sismos de distintas magnitudes es
homogénea en toda la fuente. Todo esto, con el los parámetros sismológicos en la evaluación
propósito de efectuar la actualización de los del Peligro Sísmico en el Perú, hasta ahora
parámetros estadísticos de recurrencia para planteados, deben ser reevaluados y
cada una de las fuentes, utilizando la magnitud actualizados.
Ms y la magnitud de momento Mw.
En el presente estudio se ha realizado la
I Introducción actualización de los parámetros de recurrencia
sísmica utilizando información sismológica
Teniendo en cuenta el alto nivel de riesgo registrada hasta el año 2003 y se han
sísmico a que está expuesto el Perú y empleado diversas leyes de atenuación de
reconociendo que el entorno sísmico y los aceleraciones de sismos tanto de subducción
métodos de análisis cambian continuamente, como continentales para evaluar los resultados
es necesario proveer datos adecuados obtenidos con estos nuevos parámetros. Los
destinados a estimar la peligrosidad sísmica. resultados son presentados en forma de curvas
de peligro sísmico.
DR. JORGE ALVA HURTADO En este artículo se describe la metodología de
cálculo de los parámetros de recurrencia y un
Ingeniero Civil egresado de la breve ejemplo de referencia, seleccionado para
Universidad Nacional de Ingeniería, distinguir las diferencias de los resultados con
CE MSCE PhD Doctorado en los EEUU cada una de las leyes de atenuación utilizadas.
Investigador del CISMID
Profesor Principal de Antegrado y I Evaluación de Fuentes Sismogénicas
Jefe de la Sección de Postgrado de
la FIC-UNI Una fuente sismogénica es aquella línea, zona
Consultor en Ingeniería Geotécnica o volumen geográfico que tiene similitudes
Actualmente se desempeña como geológicas, geofísicas y sísmicas tales, que se
Presidente del Capitulo de Ingeniería pueda considerar que posee un potencial
Civil del Consejo Departamental de sísmico homogéneo, es decir, en la que el
Lima del Colegio de Ingenieros del Perú. proceso de generación y recurrencia de sismos
es espacial y temporalmente
similar. La definición de las fuentes
sismogénicas y su demarcación es
de gran importancia para estudios
cuantitativos de Peligro Sísmico.

En el presente trabajo se ha
realizado una revisión de la
actividad sísmica del pasado,
reafirmando la definición de las
fuentes sismogénicas realizada
por Castillo (1993). La
determinación de estas fuentes
sismogénicas se basó en el mapa
de distribución de epicentros, así
como en las características
tectónicas de nuestro país,
agrupándose en Fuentes de
Subducción y Fuentes
Continentales. suficientes datos para modelar las fallas activas
como fuentes lineales en este tipo de análisis.
La Figura N° 1 muestra las 20 fuentes
Tabla Nº 01. coordenadas Geográficas de las Fuentes de
Subducción Superficiales y de las Fuentes Continentales sismogénicas definidas. Las Fuentes 1, 2, 3, 4
y 5 están ubicadas a lo largo de la costa y
FUENTES COORDENADAS GEOGRAFICAS (º) representan la sismicidad superficial en la zona
-80.29 +02.00 -78.32 +02.00 de Benioff (0 - 70 Km). Las Fuentes 13, 14, 15,
FUENTE 1 -81.39 -00.97 -79.65 -01.21 16, 17, 18 y 19 representan la sismicidad
-81.52 -02.39 -80.19 -02.50
intermedia (71 a 300 Km). La Fuente 20
-82.00 -03.39 -80.17 -03.45
FUENTE 2 -82.00 -06.83 -80.67 -05.42 representa la sismicidad profunda (500 a 700
-81.17 -09.00 -79.27 -07.90 Km). Las Fuentes Continentales 6, 7, 8, 9, 10,
-81.17 -09.00 -79.27 -07.90 11 y 12 están relacionadas, con la actividad
FUENTE 3
-77.00 -14.80 -75.84 -13.87
-77.00 -14.80 -75.84 -13.87 sísmica superficial andina. La Tabla N° 1
FUENTE 4 -74.16 -17.87 -73.00 -16.53 presenta las coordenadas geográficas de las
-74.16 -17.87 -73.00 -16.53 fuentes sismogénicas de subducción
FUENTE 5 -71.85 -19.87 -69.21 -19.00 superficiales y continentales y la Tabla N° 2
-71.85 -22.00 -69.21 -22.00
-77.50 +01.58 -76.92 +01.19
presenta las coordenadas geográficas para las
-79.83 -01.65 -78.90 -02.53 fuentes de subducción intermedias y profundas
FUENTE 6 -79.96 -02.46 -78.97 -03.43 (Castillo, 1993)
-80.92 -02.96 -80.79 -03.44
-78.28 -08.20 -77.86 -08.07
FUENTE 7 -77.21 -10.47 -76.83 -10.23 I Evaluación de la Recurrencia Sísmica
-75.84 -13.87 -74.76 -13.13
Para evaluar la variación del tamaño de los
FUENTE 8 -73.00 -16.53 -71.41 -14.67 eventos sísmicos que cada fuente sísmica
-73.00 -16.53 -71.41 -14.67 pueda generar es necesario determinar la
FUENTE 9 -69.71 -18.67 -68.12 -16.13 recurrencia sísmica de la fuente. La recurrencia
-76.92 +01.19 -76.50 +01.00 sísmica representa el número de eventos
FUENTE 10 -78.90 -02.53 -77.35 -02.40 mayores o iguales a alguna magnitud dentro de
-79.10 -05.20 -77.00 -04.77
-79.10 -05.20 -75.10 -04.33 la fuente y está descrita por la pendiente de la
FUENTE 11 -76.34 -10.67 -74.17 -09.33 relación de recurrencia de Gutenberg y Richter
-74.76 -13.13 -72.48 -11.40
(1954) b, la razón media anual de actividad
-74.76 -13.13 -72.48 -11.40 sísmica, la magnitud mínima y la magnitud
FUENTE 12 -68.12 -16.13 -67.76 -13.80
máxima.

Las Fuentes de Subducción modelan la

interacción de las placas Sudamericana y de Para cuantificar la relación de recurrencia de la
Nazca. Las Fuentes Continentales están actividad sísmica de la zona en estudio se
relacionadas con la actividad sísmica utilizó la expresión propuesta originalmente por
superficial andina. Se han presentado las Ishimoto-Ida en 1939 y posteriormente
fuentes como áreas, ya que no existen adecuada por Richter (1958).
 una recta en función de la densidad de datos
Log N = a - bM que existen en determinada zona. Los datos
utilizados para el cálculo de a y b se
donde: encuentran dentro de los valores de Magnitud
Mínima de Homogeneidad y Magnitud Máxima.
N = número de sismos de magnitud M o Con la finalidad de hallar la Magnitud Mínima
mayor por unidad de tiempo. de Homogeneidad se hace uso del diagrama
a,b = parámetros que dependen de la región de Gutenberg y Richter, llamada también curva
y son constantes determinadas de una regre- de recurrencia sísmica, en dicha curva se
sión no-lineal del catálogo de sismicidad local. observa un punto de inflexión cuyo valor en las
abscisas corresponde a la Magnitud Mínima de
La expresión anterior también se puede Homogeneidad. A partir de esta Magnitud
describir como: mínima se presenta un comportamiento lineal,
de esta relación lineal se obtienen los valores
Donde: N = Γ0 eβM de a y b.
La tasa es la tasa media anual de ocurrencia
Donde: de eventos mayores o iguales que la magnitud
mínima de homogeneidad.
Γ0 = 10a es el número de sismos por unidad
Tabla Nº 02. Coordenadas Geográficas de las Fuentes de Subducción
de tiempo con M > 0. Intermedias y Profundas
Β = b x ln 10
FUENTES COORDENADAS GEOGRAFICAS (º)
-78.73 +02.00 -76.00 +01.82
Para la obtención de los parámetros a y b es FUENTE 13 -81.00 -00.67 -79.59 -02.55
necesario que la base de datos disponible esté -81.00 -03.07 -79.20 -03.07
completa y lo más actualizada posible, por tal -81.00 -03.07 -79.20 -03.07
-81.93 -05.73 -78.60 -04.00
motivo se utilizó datos sismológicos desde el FUENTE 14
-79.80 -08.13 -77.17 -06.53
año 1963 al año 2003. La información -79.80 -08.13 -77.17 -06.53
sismológica fue obtenida del catálogo sísmico FUENTE 15 -76.38 -14.30 -73.86 -12.46
del Proyecto SISRA (1985) hasta el año 1992 FUENTE 16
-76.38 -14.30 -73.86 -12.46
-73.28 -16.87 -71.21 -14.40
con los datos verificados publicados por el ISC
-73.28 -16.87 -71.21 -14.40
(Internacional Seismological Centre) y FUENTE 17 -70.86 -18.80 -68.93 -15.73
actualizados hasta el año 2003 por el IGP -70.38 -22.00 -67.98 -22.00
(2003). -79.59 -02.55 -77.50 -00.73
FUENTE 18 -78.60 -04.00 -75.51 -02.06
-77.17 -06.53 -75.27 -05.33
Los parámetros estadísticos de recurrencia -77.17 -06.53 -75.27 -05.33
para cada una de las fuentes sismogénicas se FUENTE 19
-73.86 -12.46 -72.03 -11.13
han calculado utilizando la magnitud Ms y la FUENTE 20
-72.31 -06.67 -71.00 -06.33
-71.14 -11.30 -69.69 -10.93
magnitud de momento Mw, debido a que las
nuevas leyes de atenuación utilizadas están
expresadas en magnitud de momento, y se Para determinar las profundidades
requiere uniformizar la entrada de datos para la representativas de los hipocentros en las zonas
integración de la amenaza sísmica. sismogénicas se realizó un trabajo estadístico
Se utilizó la siguiente relación entre las del cálculo de frecuencias de sismos versus
magnitudes mb y Ms (Castillo, 1993). profundidad. Las Tablas N° 3 y N° 4 presentan
los parámetros de recurrencia en base a
mb = 3.30 + 0.40 Ms magnitudes Ms y Mw, utilizados para el caso de
considerar las fuentes sismogénicas de la
Para la relación entre Ms y Mw, se consideró la Figura N° 1.
obtenida por Bolaños y Monroy (2004). La
relación propuesta se presenta en tres rangos: I Leyes de Atenuación

Mw = 0.740 Ms + 1.742 Ms 6; Las leyes de Atenuación son expresiones que

Mw = 0.683 Ms + 2.039 6 < Ms <8; relacionan magnitud, posición relativa fuente-
Mw = 1.093 Ms 0.593 Ms 8 sitio e intensidad sísmica, éstas son necesarias
para evaluar los efectos que produce cada una
En el análisis estadístico de los parámetros de de las fuentes en el sitio de interés, en términos
recurrencia se utilizó el método de mínimos de intensidad sísmica.
cuadrados. Este método ajusta los valores a
 - Dahle, A.; Climent, A;
TABLA Nº 03 Parámetros sísmicos calculados en base a Magnitudes Ms
Taylor, W; Bungum, H.;
Ms Santos, P.; Ciudad Real, M.;
FUENTE PROF.(Km)
a b Mmin Mmax BETA TASA Lindholm, C.; Strauch, W y
F1 2.8532 0.3468 3.0 8.1 0.80 1.62 40 Segura, F. (1995).
F2 4.4932 0.5955 4.5 7.9 1.37 1.63 40 - Schmidt, V; Dahle, A. y
F3 4.1264 0.4836 3.5 8.0 1.11 6.79 40
F4 4.2187 0.5393 4.0 8.2 1.24 2.88 40
Bungum, H. (1997)
F5 4.1756 0.5254 4.0 8.2 1.21 2.97 40 - Youngs, R.R; Chiou,S-J.;
F6 2.7701 0.3822 2.8 7.4 0.88 1.25 40 Silva, W.J. y Humphrey, J.R.
F7 2.2523 0.4252 3.5 7.4 0.98 0.15 40
F8 3.2145 0.5414 2.8 7.0 1.25 1.25 40 (1997)
F9 3.1971 0.4536 3.5 7.5 1.04 1.02 40 - Sadigh, K.; Chang, C.-Y.;
F10 3.2445 0.4265 3.8 7.3 0.98 1.05 40 Egan, J.A.; Makdisi, M;
F11 3.7534 0.4362 3.0 7.1 1.00 6.96 40
F12 2.9039 0.4447 3.0 7.1 1.02 0.93 45
Youngs, R.R. (1997)
F13 3.0047 0.4711 3.0 6.9 1.08 0.98 125 -Casaverde y Vargas (1980)
F14 3.2452 0.4747 3.5 6.5 1.09 0.96 130
F15 3.4352 0.4352 3.8 7.2 1.00 1.51 130 Las cuatro primeras fórmulas
F16 4.4482 0.6110 4.0 7.2 1.41 2.52 115
F17 5.5512 0.6915 4.8 7.5 1.59 4.26 130 fueron generadas con la
F18 3.5942 0.4026 3.5 7.5 0.93 3.83 155 base de datos de
F19 4.3398 0.5671 4.3 7.0 1.31 1.99 160
3.0434 0.3660 4.5 7.5 0.84 0.62 580
movimientos del oeste de los
F20
EE.UU., donde sistemas de
fallas transcurrentes y de
desplazamiento de rumbo son las
TABLA Nº 04. Parámetros sísmicos Calculados en base a Magnitudes Mw fuentes sismogenéticas. Sin embargo
Ms
FUENTE a b Mmin Mmax BETA TASA PROF.(Km) corresponden a marcos tectónicos
F1 3.7217 0.4795 4.0 8.3 1.10 1.59 40
diferentes a los del Perú.
F2 6.0894 0.8392 5.1 7.4 1.93 1.61 40 Campbell (1981) desarrolló leyes de
F3 5.4177 0.6828 4.3 8.2 1.57 7.57 40 atenuación para fuentes continentales.
F4 5.6176 0.7528 4.7 8.4 1.73 3.00 40
F5 5.5403 0.7337 4.7 8.4 1.69 3.09 40
Las relaciones de atenuación de Dahle
F6 3.7363 0.5305 3.8 7.1 1.22 1.39 40 et al. (1995) para Centroamérica y de
F7 3.2533 0.5746 4.3 7.1 1.32 0.15 40 Schmidt et al. (1997) para Costa Rica,
F8 4.4890 0.7316 3.8 6.8 1.68 1.28 40 fueron establecidas para marcos
F9 4.3551 0.6309 4.3 7.2 1.45 1.10 40
tectónicos donde un proceso de
F10 4.2485 0.5763 4.5 7.0 1.33 1.13 40
F11 4.8948 0.6126 4.0 6.9 1.41 6.96 40 subducción es el fenómeno dominante.
F12 3.9508 0.6010 4.0 6.9 1.38 0.88 45
F13 4.1138 0.6367 4.0 6.8 1.47 0.92 125 Youngs, Chiou, Silva y Humphrey
F14 4.3627 0.6415 4.3 6.5 1.48 1.01 130
F15 4.5287 0.6015 4.5 7.0 1.39 1.66 130 (1997), han desarrollado relaciones de
F16 5.8866 0.8257 4.7 7.0 1.90 2.53 115 atenuación para la aceleración máxima
F17 7.4704 0.9855 5.3 7.2 2.27 4.42 130 del suelo y para la aceleración
F18 4.6431 0.5638 4.3 7.2 1.30 4.14 155 espectral de respuesta, para dos tipos
F19 5.7901 0.7880 4.9 6.8 1.81 2.12 160 de terremotos que ocurren en la zona
F20 4.0533 0.5213 5.1 7.2 1.20 0.62 580
de subducción: terremotos de interfase
Para evaluar la variabilidad de resultados al ser
y terremotos intraplaca que tengan una
integrados con diferentes leyes de atenuación y
magnitud de momento M=5 o mayor, y para
debido a que no existe una relación de
distancias de 10 a 500 kilómetros. Estas
atenuación generada con los movimientos
relaciones fueron desarrolladas por medio de
fuertes del terreno registrados en el Perú
un análisis de regresión utilizando una base de
mediante las técnicas estadísticas aceptadas
datos de 174 terremotos de las características
en el .campo profesional para derivar dichas
anotadas que han ocurrido en diversas partes
expresiones matemáticas, se utilizaron diversas
del mundo, tales como Alaska, Chile, Cascadia,
relaciones de atenuación y se compararon los
Japón, México, Perú y Salomon Islands. Esta
resultados que de ahí se derivan. Entre las más
ley se aplica a sismos de subducción.
conocidas y utilizadas se encuentran las
siguientes:
Las ecuaciones dadas por Sadigh, Chang,
- Idriss, I. M. (1993)
Egan, Makdisi y Youngs (1997), han
- Sadigh, K. et al. (1993)
desarrollado relaciones de atenuación para la
- Boore, D.; Joyner, W. y Fumal, T. (1993)
aceleración máxima del suelo y para la
- McGuire(1974)
aceleración espectral de respuesta, utilizando
- Campbell,K.W(1981)
una base de datos de terremotos de California.
Esta ley de atenuación se aplica a sismos
continentales.
La ley propuesta por Casaverde y Vargas
(1980) está basada en los registros de los
acelerógrafos de las componentes
horizontales de diez sismos peruanos
registrado en Lima y alrededores, siendo
evidente que existe escasez de datos de
registros de aceleraciones en el Perú.

I Determinación del Peligro Sísmico

El Peligro Sísmico puede calcularse

considerando la suma de los efectos de la
totalidad de las fuentes sísmicas y la
distancia entre cada fuente y el sitio donde
se encuentra la estructura. Idriss, 1993 se aplicó a las 20 fuentes; JB93
significa que la relación de atenuación de
Para la determinación del Peligro Sísmico se Boore, Joyner y Fumal, 1993 se aplicó a todas
han utilizado los programa de cómputo: RISK y las fuentes; S93 significa que la relación de
EZRISK desarrollados por R. McGuire (1976, atenuación de Sadigh et al., 1993 se aplicó a
1995), y el programa MRIESGO desarrollado todas las fuentes; S97 significa que la fórmula
por Bolaños y Monroy (2004). La diversidad de de Schmidt, 1997 para subducción se utilizó en
programas de cómputo utilizados se debe a su las fuentes de subducción y la fórmula de
capacidad para incluir las leyes de atenuación Schmidt, 1997 para eventos someros se aplicó
propuestas. a las fuentes continentales; YO97 + CA81
significa que la ley de atenuación de Youngs et
Estos programas se utilizaron para evaluar las al. (1997), se aplicó a las fuentes de
probabilidades de excedencia correspondientes subducción y la fórmula de Campbell (1981), se
a determinados niveles de aceleración pico. El aplicó a las fuentes continentales; YO97+SA97
período de retorno de la aceleración pico es significa que la ley de atenuación de Youngs et
simplemente el inverso de la probabilidad anual al. (1997), se aplicó a las fuentes de
de excedencia. Se utilizaron las relaciones de subducción y la fórmula de Sadigh et al. (1997),
atenuación arriba mencionadas para evaluar la se aplicó a las fuentes continentales.
variabilidad que pueden tener los resultados al CV80+MG74 significa que la ley de atenuación
ser integrados con diferentes relaciones de de Casaverde y Vargas, 1980 se aplicó a las
fuentes de
TABLA Nº 5. aceleraciones Máximas esperadas (% G) obtenidas con diferentes Leyes de Atenuación subducción y la
propuesta por
PERIODO DE RETORNO
McGuire, 1974 se
LEY DE ATENUACION aplicó a las fuentes
100 200 500 1 000 2 000 3 000 5 000 10 000
continentales. Para
Idriss, I. M. (1993) I 93 0.05 0.07 0.10 0.13 0.15 0.16 0.19 0.21 estas últimas leyes
Sadigh, K. et al. (1993) S 93 0.04 0.05 0.07 0.08 0.10 0.11 0.14 0.15
se utilizaron los
valores de
Boore, Joyner y Fumal (1993, 4): JB 93 0.23 0.29 0.38 0.46 0.53 0.57 0.66 0.76 atenuación y
Dahle, A. et al. (1995) DA 95 0.10 0.13 0.15 0.18 0.20 0.22 0.25 0.29 recurrencia sísmica
en función de Ms.
Schmidt, et al. (1997): S 97 0.13 0.14 0.18 0.20 0.23 0.25 0.28 0.32
Para las otras
atenuación. leyes se utilizaron las magnitudes de momento
Mw.
En la Tabla N°5 y la Figura N°2 se muestran los Se puede observar la amplia dispersión de los
resultados obtenidos para la Coordenada resultados. Es evidente que algunas de las
Longitud -71.33°, Latitud -14.90°. El símbolo leyes de atenuación utilizadas no son
DA95 significa que la relación de atenuación de aplicables al caso en estudio. Por ejemplo I93,
Dahle et al., 1995 se aplicó a todas las fuentes; S93 y JB93 fueron desarrolladas para la costa
I93 significa que la relación de atenuación de oeste de EEUU y producen en nuestro caso la
mayor variabilidad de resultados. Las otras
leyes son para sismos de subducción o han
sido utilizadas en combinación con sismos
continentales, pudiéndose observar que el
rango de variación de resultados es mucho
menor.

I Conclusiones y Recomendaciones

En el presente trabajo se ha realizado una

revisión de la actividad sísmica del pasado,
reafirmando la definición de las fuentes
sismogénicas realizada por Castillo (1993).
Para la actualización de los parámetros de
recurrencia sísmica se utilizó información
sismológica obtenida del catálogo sísmico del
Proyecto SISRA (1985) hasta el año 1992 con
los datos verificados publicados por el ISC
(Internacional Seismological Centre) y
actualizados
hasta el año 2003 por el IGP. En el estudio
probabilístico de peligro sísmico se han
considerado las fuentes sismogénicas como
áreas. Se han utilizado los catálogos sísmicos y
determinado recurrencias sísmicas con
magnitudes Ms y Mw actualizados al año 2003.

Se han utilizado diversas fórmulas de

atenuación publicadas en la literatura técnica
para evaluar los movimientos fuertes del
terreno en el sitio generados por los eventos
sísmicos ocurridos en las fuentes. Entres éstas
se incluyen: Casaverde y Vargas (1980), R.;
McGuire (1974); Idriss, I.M. (1993); Sadigh, K.
et al. (1993); Boore, Joyner y Fumal (1993,4);
Campbell, K.W. (1981); Dahle, et al. (1995);
Schmidt, et al. (1997); Youngs, et al. (1997);
Sadigh, et al. (1997).

Se puede observar la amplia variabilidad de los

resultados. Es evidente que algunas de las
leyes de atenuación utilizadas no son
aplicables al Perú. Por ejemplo las leyes de
atenuación I93, S93 y JB93 fueron
desarrolladas para la costa oeste de EEUU y
producen en nuestro caso la mayor variabilidad
de resultados. Las otras leyes son para sismos
de subducción o han sido utilizadas en
combinación con sismos continentales,
pudiéndose observar que el rango de variación
de resultados es mucho menor.

Reproducimos el artículo con autorización del

autor.
En: Ingeniería Civil, Nº 03. nov. 2005.