¿Es Posible Aplicar las Redes Neuronales Artificiales en el Pronóstico Sísmico?
Análisis Preliminar Aplicando la Topología Radial. Caso: México
Can Artificial Neural Networks be Applied in Seismic Predicition? Preliminary
Analysis Applying Radial Topology. Case: Mexico
Cinthya Mota-Hernández1, Luis Esquivel-Rodríguez1 and Rafael Alvarado-Corona2
1Universidad
2SEP,

Autónoma del Estado de México, México. www.uaemex.mx

SARACS Research Group, www.saracs.com.mx, Ciudad de México, México. E-mail: ralvcor@gmail.com

Resumen--- Sismos tectónicos de elevada magnitud pueden
causar pérdidas considerables en términos de vidas humanas,
económicos y de infraestructura, entre otros. De acuerdo con
una evaluación publicada por el Servicio Geológico de los
Estados Unidos, 30 es el número de sismos que han impactado
en gran medida a México desde finales del siglo XIX al
actual.. En base a datos del Servicio Sismológico Nacional, en
el periodo comprendido entre el 1 de Enero del 2006 y el 1 de
Mayo del 2013 se han presentado 5,826 sismos mayores a 4.0
grados de magnitud en la escala de Richter (28.54% del total
de sismos registrados en el territorio nacional), siendo la placa
del Pacífico y la placa de Cocos las de mayor importancia. La
investigación describe el desarrollo de una Red Neuronal
Artificial (RNA) basada en la topología radial que busca
generar una predicción con un margen de error menor al
20% que informe sobre la probabilidad de un sismo futuro,
una de las preguntas principales del enfoque propuesto es:
¿Es posible aplicar las redes neuronales artificiales en el
pronóstico sísmico? Se puede argumentar que la investigación
tiene el potencial de aportar en el pronóstico sísmico, más
investigación es necesaria para consolidar datos y contribuir a
mitigar el impacto causado por tales eventos al vincularse con
la sociedad.
Palabras clave--- Análisis, México,
neuronales artificiales, Sismicidad.

Pronóstico,

Redes

Abstract--- Tectonic earthquakes of high magnitude can cause
considerable losses in terms of human lives, economic and
infrastructure, among others. According to an evaluation
published by the U.S. Geological Survey, 30 is the number of
earthquakes which have greatly impacted Mexico from the
end of the XIX century to this one. Based upon data from the
National Seismological Service, on the period between
January 1, 2006 and May 1, 2013 there have occurred 5,826
earthquakes which magnitude has been greater than 4.0
degrees on the Richter magnitude scale (25.54% of the total of
earthquakes registered on the national territory), being the
Pacific Plate and the Cocos Plate the most important ones.
This document describes the development of an Artificial
Neural Network (ANN) based on the radial topology which
seeks to generate a prediction with an error margin lower

than 20% which can inform about the probability of a future
earthquake one of the main questions is: can artificial neural
networks be applied in seismic forecasting? It can be argued
that research has the potential to bring in the forecast seismic,
more research is needed to consolidate data and help mitigate
the impact caused by such events linked with society.
Keywords--- Analysis, Mexico, Neural Artificial Networks,
Seismicity.

I. INTRODUCCION
Los desastres naturales pueden causar alto impacto
en sociedades vulnerables, existe investigación que se ha
realizado en varias facetas [9], por otro lado, para el
monitoreo y alerta de actividad sísmica en México se
cuenta a partir de 1991 con el Sistema de Alerta Sísmica
(SAS) [1], que cubre las zonas aledañas a Guerrero y
cuenta con un total de 12 estaciones sensoras, además del
Sistema de Alerta Sísmica de la Ciudad de Oaxaca
(SASO), desarrollado en el 2003 y con un total de 36
estaciones sensoras. Estos sistemas han emitido 78 alertas
tempranas de más de 2350 sismos [2], hecho que ha
permitido contemplar la posibilidad de extender la
cobertura de los mismos a los estados comprendidos entre
Jalisco y Chiapas (la mayor parte del territorio se ubica
sobre la placa de Cocos). De cualquier manera, ambos
sistemas reaccionan en cuanto la actividad sísmica es
detectada, sin embargo, se evaluó la posibilidad de emplear
una RNA capaz de buscar un patrón en los registros de
sismos desde los 4.0 grados de magnitud a partir del año
1998 hasta el 1 de Mayo del 2013, la RNA se basó en una
topología radial. Haciendo uso de una muestra de poco más
de 5,000 sismos, misma que fue dividida para integrar los
conjuntos de validación y producción, se llevó a cabo el
entrenamiento de la RNA mediante el aprendizaje
supervisado con el objetivo de obtener a sus salidas la
magnitud y otros datos de un posible sismo futuro. El
presente artículo explica el desarrollo de dicha red con base

en una metodología de análisis para sistemas suaves [3],
además de los resultados obtenidos del entrenamiento.
Existen diversas investigaciones en el área sísmica con
enfoques que enriquecen y aportan de manera integral al
análisis del fenómeno [10].

El diseño elegido para estructura de la red, considerando el
tratamiento de datos para la investigación, se ilustra en la Fig.1.

II. RECOPILACIÓN Y TRATAMIENTO DE DATOS
El primer proceso llevado a cabo fue la adquisición y
procesamiento de los datos de tal manera que el entrenamiento
contara con un número aceptable de muestras, cada una de ellas
completa y en una escala coherente general. Los datos fueron
extraídos del Servicio Sismológico Nacional [4] y se sitúan entre
el 2 de Marzo del año 2006 y el 1 de Mayo del año 2013,
integrando una muestra de un total de 5,798 registros sísmicos,
población que fue segmentada de la siguiente forma:
TABLA I
Segmentación de datos

No.
muestras
% Total

de

Fig. 1 Estructura de la red
Fuente: Elaboración propia

Conjunto de
entrenamiento
3001

Conjunto de
validación
1551

Conjunto de
producción
1246

51.76%

26.75%

21.49%

Tras la formación de los grupos se llevó a cabo un análisis
sobre la importancia para el caso de los campos que conformaban
las muestras extraídas, optando por el empleo de las siguientes
variables:
TABLA II
Variables de entrada y salida
Entradas
Fecha
Hora
Minuto
Zona

Salidas
Latitud
Longitud
Profundidad
Magnitud

El proceso de conversión de datos tanto en tipo como en
escala se llevó a cabo en todos los campos que representarían las
variables de entrada para la red, y en todos los conjuntos
obtenidos en el paso anterior.

III. DESARROLLO DE LA RED NEURONAL ARTIFICIAL
Partiendo del análisis del problema y con base en el
planteamiento de las variables analizadas, se procedió a realizar el
modelado de las medidas de desempeño de la red a entrenar, así
como el análisis funcional de la misma. Para el entrenamiento se
planteó el uso de aprendizaje supervisado, con una terminación
basada en el error, para propósitos de exactitud se establece que el
error fuera causa de terminación del entrenamiento cuando se
alcanzara la condición de sobreentrenamiento; el error máximo
planteado para las salidas de la red fue del 20%, cualquier salida
por encima de ese rango se consideraría como razón para
clasificar la red como no óptima para su uso.

A priori, era deseable establecer el modelo de red que fuera
más óptimo para la solución del problema planteado, tanto en el
tiempo requerido como en el error mínimo obtenido a la salida de
la red. Este proceso se basó en los aspectos siguientes:
A. Modelos de Redes elegidos.
1.
2.
3.
4.
5.

Perceptrón multicapa, modelo conocido por ser
apropiado para una gran cantidad de problemas [5].
Recurrente en series de tiempo [6].
Recurrente generalizada, similar a la anterior y
comúnmente usada en predicción [7].
Red de base radial con perceptrón multicapa [8].
Radial, de propósito general.

B. Número de elementos de procesamiento (EP).
Las redes preliminares entrenadas poseían 7 capas ocultas,
cuyas cantidades de EP son descritos en la Tabla III.
TABLA III
EP por capa
Capa 1
Capa 2
Capa 3
Capa 4
Capa 5
Capa 6
Capa 7

No. de EP
20
60
100
150
100
60
20

C. Total de Muestras a usar.
En todas las redes preliminares se hizo uso de una población
de 200 muestras para el entrenamiento.
D. Reglas de Propagación y Funciones de Activación.

Se hizo uso de las reglas de propagación momentum para el
perceptrón multicapa entrenado (ésta regla es propia de los
perceptrones multicapa) y Quickprop (que intenta hacer uso de un
modelo cuadrático de la superficie de error para acelerar la
convergencia).
E. Criterio de convergencia.
Se estableció como criterio de terminación el mismo a
utilizar en la red final, es decir, terminación basada en el error
mínimo antes de alcanzar el sobreentrenamiento. Se estableció
como error mínimo el 10%. Al finalizar el entrenamiento de las
redes antes mencionadas, se obtuvieron los siguientes resultados a
fin de seleccionar el modelo de red más óptimo para su aplicación
se llevaron a cabo las comparativas siguientes:

Fig. 4 Comparación de conjuntos de producción
Fuente: Elaboración propia

La comparación entre los conjuntos de entrenamiento y
validación demostró que la red 4 (red de base radial con
perceptrón multicapa) era la más adecuada para resolver el
problema. La red final se basó por consiguiente, en una topología
radial con perceptrón multicapa con los siguientes datos:
TABLA V
Detalles de la red final
No. de EP

Fig. 2 Comparación de conjuntos de validación
Fuente: Elaboración propia

Capa de
entrada
Capa 1

100

Capa 2

100

Capa 3

200

Capa 4

400

Capa 5

200

Capa 6

100

Capa 7

50

50

Función de
transferencia
Gaussiana

Regla de
propagación
Quickprop

Tangente
hiperbólica
Tangente
hiperbólica
Tangente
hiperbólica
Tangente
hiperbólica
Tangente
hiperbólica
Tangente
hiperbólica
Tangente
hiperbólica

Quickprop
Quickprop
Quickprop
Quickprop
Quickprop
Quickprop
Quickprop

Obteniendo, a sus salidas lo siguiente:
TABLA VI
Salidas de la red final

Fig. 3 Comparación de conjuntos de entrenamiento
Fuente: Elaboración propia

Dato
Ciclos requeridos
Tiempo requerido
Conjunto de producción
Error promedio Latitud
Error promedio Longitud
Error promedio Magnitud
Conjunto de validación
Error promedio Latitud
Error promedio Longitud
Error promedio Magnitud

Red final
4411
6:24:22
2.40%
5.14%
0.59%
2.70%
5.32%
0.53%

Conjunto de entrenamiento
Error promedio Latitud
Error promedio Longitud
Error promedio Magnitud
MSE
NMSE
% Error

concluyentes.
6.56%
8.32%
0.52%
0.161018021136
1.066836403099
12.519098354399

AGRADECIMIENTOS.- Los autores desean agradecer a todas y
cada una de las instituciones, organizaciones y personas que han
brindado apoyo en la realización del presente proyecto de
investigación.

REFERENCIAS
La comparación generada entre los valores actuales y los
generados por la red neuronal, se muestran en la Fig. 5

Fig. 5 Gráfica comparativa de los valores de las entradas y salidas
de la red
Fuente: Elaboración propia

[1]

J. M. Espinosa-Aranda, A. Cuéllar, A. García, G. Ibarrola, R.
Islas, S. Maldonado and F. H. Rodríguez, “Evolution of the
Mexican Seismic Alert System (SASMEX) (in Earthquake
early warning )”, Seismological Research Letters, 80(5):694706, Sep. 1999.

[2]

J.M. Espinosa-Aranda, A. Cuéllar, F.H. Rodríguez, B.
Frontana, G. Ibarrola, R. Islas, A. García, “The seismic alert
system of Mexico (SASMEX): Progress and its current
applications”, Soil Dynamics and Earthquake Engineering,
Vol. 31(2): 154-62, Feb. 2011.

[3]

Checkland, P., and Scholes, J, La metodología de sistemas
suaves en acción, México: Ed. Noriega Editores, 1994.

[4]

Servicio Sismológico Nacional (2013), México,
Recuperado de http://www.ssn.unam.mx/.

[5]

A.G. Parlos, K.T. Chong, and A.F. Atiya, "Application of the
recurrent multiplayer perception in modeling complex process
dynamics", IEEE Trans. on Neural Networks: Special issue on
Recurrent Dynamic Neural Networks, Vol. 5(1): 255-66, Jan.
1994.

[6]

Hui Liu, Jun-An Lu, Jinhu Lü, David J. Hill, “Structure

IV. DISCUSIÓN
En el proceso que conlleva a la elección del modelo de RNA
más adecuado para la resolución de un problema, es necesario
evaluar aspectos como el tipo de problema a resolver, el error
mínimo esperado a la salida, el número de muestras y el
tratamiento-conversión de los datos contenidos. Pese a que la
elección del número de neuronas además del número y
composición de las capas ocultas es un factor limitante para
ciertos equipos informáticos en cuanto a los recursos con que se
cuenta, su elección correcta es clave para alcanzar el criterio de
convergencia establecido. El proceso de elección del mejor
modelo de RNA implica prueba y error utilizando distintos
modelos y/o variaciones en su aplicación. Las redes neuronales
artificiales tienen potencial aplicación en el ámbito sísmico.

identification

of uncertain

networks with time delay”,

general

complex

dynamical

Automatica, Vol. 45(8):1799-

1807, Aug. 2009.
[7]

Panakkat, A. and Adeli, H., “Recurrent Neural Network for
Approximate Earthquake Time and Location Prediction Using
Multiple Seismicity Indicators”, Computer-Aided Civil and
Infrastructure Engineering, Vol. 24(4): 280–92, May 2009.

[8]

Jayawardena A.W., et al., “Comparison of Multilayer
Perceptron and Radial Basis Function networks as tools for
flood forecasting”. in Destructive Water: Water-Caused
Natural Disasters, Their Abatement and Control , Leavesley
GH, et al. editors. Anaheim, CA: IAHS Publication
Number 239, pp. 107–13.

[9]

ARISTIZABAL, EDWIN FABIAN GARCIA, CARLOS
ALBERTO RIVEROS JEREZ, and MANUEL ALONSO
BUILES
BRAND.
"INFLUENCE
OF RAINFALL
INTENSITY ON INFILTRATION AND DEFORMATION
OF UNSATURATED SOIL SLOPES INFLUENCIA DE LA
INTENSIDAD DE LA LLUVIA EN LA INFILTRACION Y
DEFORMACION DE
TALUDES
PARCIALMENTE
SATURADOS." Dyna 170 (2011): 117.

[10]

Alvarado-Corona, Rafael, and Jaime Santos-Reyes. "Applying
MORT to the Analysis of the Haiti’s Earthquake." DISASTER
ADVANCES 5, no. 4 (2012): 102-109.

V. CONCLUSIÓNES
Un análisis de la aplicación de Redes Neuronales artificiales
aplicadas al ámbito sísmico fue presentado. En conclusión, se
obtuvo un error mínimo ubicado por debajo del 10% en las tres
variables de salida, pese a que una de las salidas planteadas fue
descartada, siendo el modelo de topología radial con perceptrón
multicapa el de mayor eficiencia para la resolución del problema.
La variable que en todos los casos presentó el menor error fue la
magnitud sísmica; no obstante haber excluido variables de salida
referentes a tiempo, se plantea llevar a cabo mayor investigación
sólida que invite a la reflexión y pueda brindar resultados

