Sismo Impulsivo Carlos Fraud
Sismo Impulsivo Carlos Fraud
Sismo Impulsivo Carlos Fraud
852
1 INTRODUCCION
Los movimientos sismicos cercanos a la falla sismognica, han sido tema de mucha
investigacin en aos recientes, se puede decir que las vibraciones del suelo proximas a la
ruptura que genera un terremoto pueden ser caracterizadas por movimientos impulsivos de
corta duracin, que exponen a las estructuras un alto input de energa. (Krawinkler et al.
,2005). La zona comprendida dentro de una estrecha franja a cada lado la traza de la falla se
conoce como zona cercana o prxima a fallas (near-fault, near-field o near-source). Dentro de
esta zona, el movimiento del suelo est fuertemente influenciado por el mecanismo y la
direccin de propagacin de la ruptura respecto del sitio como as tambin por posibles
desplazamientos permanentes del suelo como resultado del deslizamiento relativo de ambas
caras de la falla. (Sasani, M. And Bertero, 2000)
Estos movimientos con pulsos de velocidad intensos pueden afectar severamente el
desempeo ssmico de estructuras (Anderson y Bertero, 1987; Hall et al, 1995;. Alavi y
Krawinkler, 2000). El perodo del pulso en relacin con el perodo fundamental de la
estructura tambin afecta en gran medida el rendimiento estructural (Anderson y Bertero,
1987). Los registros de terremotos obtenidos en estaciones ubicadas en las proximidades de la
falla sismognica muestran significativas diferencias con aquellos registrados en estaciones
que se encuentran a cierta distancia.
La presencia de pulsos de aceleracin y velocidad en los registros sismicos le confiere alto
poder destructivo. A la fecha distintos investigadores han propuestos formas diversas para
evaluar si un registro es de caracter impulsivo o no. En el presente trabajo se describe el criterio
mayormente aceptado internacionalmente para catalogar a un acelerograma como impulsivo, el cual
utiliza la transformada wavelet a la vez que se presenta un nuevo mtodo de bajo costo
computacional y eficaz para identificar registros ssmicos .
2 PULSO DE VELOCIDAD
El termino pulso ha sido usado como referencia de la aceleracin, velocidad y
desplazamientos de grandes movimientos del suelo debido a terremotos. Para este estudio los
pulsos son considerados en un registro como ondas de baja frecuencia, con temprana
ocurrencia en la historia del tiempo de la velocidad y con una gran amplitud.
Varias investigadores han desarrollado modelos analticos detallados, describiendo las
formas de pulsos de velocidad resultante de directividad (Fu and Menun 2004; Makris and
Black 2004; Mavroeidis and A. Papageorgiou 2003; Rodrguez-Marek and Bray 2004).
Ninguno de estos modelos, sin embargo, son capaces de determinar si un movimiento del
suelo arbitrario contiene un pulso, sino que slo se utilizan para caracterizar los pulsos que
han sido previamente identificados por el juicio del usuario (Baker, 2007).
Por lo general, los usuarios clasifican los registros de forma manual usando su mejor juicio
visual, pero esto da lugar a clasificaciones que varan de un autor a otro (Fu and Menun 2004;
Mavroeidis and A. Papageorgiou 2003; Somerville 2003; Akkar et al. 2005; Cox and Ashford
2002).
3 CRITERIO DE BAKER PARA CLASIFICAR REGISTROS IMPULSIVOS
Un enfoque del proceso de seales basado en ondas wavelet propuestas por J. Baker 2007,
2008; es usado para identificar empricamente pulsos en registros ssmicos. El procedimiento
consiste en "aislar" o "extraer" los pulsos contenidos en un registro de velocidades utilizado
anlisis de ondas wavelet.
Una comprensin bsica puede ser lograda a travs de la comparacin con el anlisis de
Copyright 2013 Asociacin Argentina de Mecnica Computacional http://www.amcaonline.org.ar
853
Fourier. El anlisis de Fourier representa una seal que utiliza una combinacin lineal de
ondas sinusoidales, cada una representa una seal de longitud infinita y una sola frecuencia.
En contraste, el anlisis wavelet se descompone en una seal de ondas que se localizan en el
tiempo y que representan un rango estrecho de frecuencias. Para seales no estacionarias tales
como movimientos ssmicos, puede ser ventajoso para representar la seal como una suma de
pequeas ondas en lugar de una suma de ondas sinusoidales estacionarias. (J. Baker, 2007).
En el anlisis de Fourier las ondas seno y coseno sirven como funcin base y cada onda es
muy precisa en cuanto a la frecuencia que representa, pero no precisa al todo en el intervalo de
tiempo que representa (por ser de longitud infinita). En contraste, la funcin base wavelet, es
ms precisa en el tiempo y el rango de frecuencia que representa. La localizacin del tiempo
es particularmente ventajosa con estudios de fenmenos de corta duracin tal como pulsos de
velocidad. (S. Shahi and J.Baker, 2010).
Una funcin bsica de ondas wavelet definida en el tiempo:
s.l (t ) =
1 t l
s s
(1)
C s ,l =
f (t ) s.l (t )dt =
f (t )
1 t l
dt
s s
(2)
El mtodo para el anlisis de las seales digitales, tales como la historia de la velocidad del
suelo, utiliza la Ec. 2 en forma discretizada, tomando una suma en lugar de una integral.
La utilidad de este procesamiento de seales es que las funciones bsicas de la
transformada de wavelet son similares a la forma de los pulsos de velocidad, causado por
directividad. El pulso de velocidad se muestra sobre la transformada wavelet, como un
coeficiente importante de la onda que tiene asociado una escala y ubicacin con el pulso, as
como se extrae el pulso del movimiento del suelo.
Para cada registro ssmico analizado, se extrae el mayor pulso de velocidad, usando la
descomposicin de wavelet. Ello se muestra en la figura 1, en la cual se observa claramente
como captura el pulso de la velocidad del registro original. El registro del pulso extrado se le
sustrae al registro original obteniendo el registro residual. Para registros no impulsivos el
registro residual es prcticamente igual al movimiento del registro ssmico original.
854
Figura 1: Ilustracin del procedimiento de utilizado para extraer la parte de pulso para el terremoto de 1979
Imperial Valley estacin El Centro Array #5. (J. Baker 2008)
1
1+ e
(3)
La primera variable es la relacin del cociente de la velocidad mxima del suelo (PGV) del
registro residual y del registro original, la segunda variable es el cociente de la energa del
registro residual y la energa del registro original, donde la energa se considera como
cuadrado de la velocidad acumulada de la seal.
Los valores varan entre 0 y 1; para registros con puntajes por encima de 0,85 se consideran
que proporcionan un fuerte indicador de que el movimiento del suelo es tipo pulso y por
debajo de 0,15 se consideran que el movimiento es del tipo no pulso, y entre 0.85 y 0.15 son
clasificados como ambiguos.
A su vez considera que los efectos de directividad son de inters primordial, lo cual
incorpora un criterio adicional a fin de identificar los pulsos que llegan temprano en la
evolucin temporal de velocidad del suelo. Los pulsos tardos se identifican mediante el
clculo de la velocidad cuadrtica acumulada del registro original y del registro del pulso
extrado.
Copyright 2013 Asociacin Argentina de Mecnica Computacional http://www.amcaonline.org.ar
855
CSV (t ) = V 2 (u )du
(4)
Figura 2: Se muestra el tiempo donde se produce el 10 % de CSV del registro del pulso extrado y el tiempo del
20% del CSV del registro original. Registro de Imperial Valley 1979 estacin FF Overpass Meloland CE.
(J. Baker 2007)
Se calcula el CSV para el registro original y el registro del pulso extrado. Luego se
identifican los valores de los tiempos donde se produce el 10% del CSV del registro del pulso
extrado y el tiempo donde se produce el 20% del CSV del registro original. Para que el
registro se considere tipo pulso el tiempo del CSV del registro del pulso extrado debe ser
menor al tiempo del CSV del registro original. De esta manera se asegura que el pulso se
inicia antes de un valor significativo de la velocidad cuadrtica acumulada del registro
original. En la figura 2 se ilustra una muestra del procedimiento descripto.
4 METODO PROPUESTO PARA IDENTIFICAR REGISTROS IMPULSISVOS
El presente trabajo propone un mtodo cuantitativo para identificar registros impulsivos, de
forma tal, que sus resultados reflejen buena correlacin con la observacin manual o visual
por parte del investigador. El mtodo se concentra en el anlisis del trazo de la historia de
tiempo de velocidades obtenidas por integracin de las aceleraciones medidas por los
instrumentos ssmicos.
El mtodo se orienta a identificar el carcter impulsivo del registro de velocidad, pero es
posible aplicarlo tambin a los registros de aceleracin o desplazamiento. El mismo est
basado en la determinacin de un nuevo parmetro para caracterizar el registro de la velocidad
de cada una de las componentes horizontales.
856
l i =1 xi + yi
n
(5)
Ld v = i =1 t i + vi
n
(6)
Ip =
Ld v
PGV
(7)
Siendo PGV , el valor mximo absoluto del registro de velocidad en unidades de cm/s. La
definicin captura de forma sencilla y eficaz el aspecto impulsivo que se visualiza
manualmente en un registro.
Copyright 2013 Asociacin Argentina de Mecnica Computacional http://www.amcaonline.org.ar
857
5 REGISTROS SSMICOS
A los fines de evaluar la efectividad del ndice de Impulsividad se seleccion un set de 17
registros de aceleracin, los cuales se integraron en el tiempo a fin de obtener la historia del
tiempo de velocidades.
Los registros ssmicos utilizados corresponden a terremotos destructivos con magnitudes
momento comprendidas entre 5.7 y 8.0. La seleccin es un subconjunto de los 398 registros
ssmicos utilizados por J. Baker a fin de poder contrastar los resultados y fueron obtenidos de
la base de datos internacional PEER NGA (Pacific Earthquake Engineering Research Center Next Generation Attenuation). La tabla 1, muestra los datos sismolgicos de los registros
seleccionados.
Los registros corresponden a la componente horizontal, dicha componente es normal a la
traza de la falla, de los cuales cinco fueron generados por fallas tipo inversa oblicua, seis a
fallas de transcurrente, cinco a fallas tipo Inversa, y uno a fallas normal.
Los datos de aceleracin del suelo en funcin del tiempo (acelerogramas), tienen definido
un paso de tiempo t constate, por lo cual cada incremento de la aceleracin queda definido
por un intervalo de tiempo t en segundos.
Los pasos de tiempo t de los registros utilizados varan de 0.005; 0.003; 0.01; 0.02
segundos, dependiendo de la calibracin del acelermetro registrador, los equipos modernos
comnmente registran 200 muestras por segundo.
NOMBRE
DE TERREMOTO
ESTACIN
DIST.
PLANO
AO
FALLA
(Km)
DIST.
EPIC.
(Km)
MW
TIPO DE
FALLA
Coyote Lake
1979
20.67
23.91
5.74
Desgarre
Coyote Lake
1979
20.67
23.91
5.74
Desgarre
Chi-Chi, Taiwan
TCU117
1999
25.44
48.33
7.62
Inversa oblic.
Chi-Chi, Taiwan
WGK
1999
9.96
31.96
7.62
Inversa oblic.
Imperial Valley-06
Chihuahua
1979
7.29
18.88
6.53
Desgarre
San Fernando
Lake Hughes #9
1971
22.57
23.10
6.61
Inversa
Irpinia, Italy-01
Bovino
1980
46.25
51.61
6.90
Normal
1987
20.79
15.29
5.99
Inversa oblic.
Cape Mendocino
Cape Mendocino
1992
6.96
10.36
7.01
Inversa
Northridge-01
1994
5.92
20.27
6.69
Inversa
Loma Prieta
1989
9.31
27.05
6.93
Inversa oblic.
Imperial Valley-06
El Centro Array #4
1979
7.1
27.1
6.5
Desgarre
Erzican, Turkey
Erzincan
1992
4.4
9.0
6.7
Desgarre
Northridge-01
1994
5.5
21.6
6.7
Inversa
Chi-Chi, Taiwan
TCU068
1999
0.3
47.9
7.6
Inversa oblic.
Imperial Valley-06
EC Meloland Overpass FF
1979
0.1
19.4
6.5
Desgarre
Chi-Chi, Taiwan-03
TCU076
1999
14.7
20.8
6.2
Inversa
858
ESTACIN
PGV
(s)
(cm/s)
Ld v
Ip
Coyote Lake
0.005
5.87
177.9
30.3
Coyote Lake
0.005
5.47
173.2
31.6
Chi-Chi, Taiwan
TCU117
0.005
57.81
1522.3
26.3
Chi-Chi, Taiwan
WGK
0.005
67.21
1857.1
27.6
Imperial Valley-06
Chihuahua
0.01
30.44
1167.0
38.3
San Fernando
Lake Hughes #9
0.01
4.55
251.3
55.3
Irpinia, Italy-01
Bovino
0.003
2.76
245.0
88.9
0.02
36.51
498.9
13.7
Cape Mendocino
Cape Mendocino
0.02
58.40
980.7
16.8
Northridge-01
0.02
119.77
1641.3
13.7
Loma Prieta
0.005
71.47
989.3
13.8
Imperial Valley-06
El Centro Array #4
0.005
77.9
762
9.8
Erzican, Turkey
Erzincan
0.005
95.4
865
9.1
Northridge-01
0.01
87.7
756
8.6
Chi-Chi, Taiwan
TCU068
0.005
191.1
1822
9.5
Imperial Valley-06
EC Meloland Overpass FF
0.005
114.9
816
7.1
Chi-Chi, Taiwan-03
TCU076
0.005
59.3
601
10.1
Ip
Ip < 12
12 < Ip < 20
20 < Ip < 30
Ip > 30
859
Clasificacin
Fuertemente Impulsivo
Impulsivo
Moderadamente Impulsivo
No Impulsivo
860
a)
100
b)
PGV=36.51
100
Pulso Indicador=1
Pulso Tardio=V
Pulso Indicador=0.32
Pulso Tardio=F
50
cm/s
cm/s
50
PGV=58.52
-50
-50
Ip=13.7
-100
10
15
20
Tiempo
25
30
Ip=16.8
35
-100
40
10
15
20
Tiempo
c)
100
30
100
Pulso Indicador=0.99
Pulso Tardio=V
40
PGV=71.47
Pulso Indicador=0.45
Pulso Tardio=F
cm/s
50
-50
-50
Ip=13.7
-100
35
d)
PGV=119.77
50
cm/s
25
10
15
20
Tiempo
25
30
Ip=13.8
35
-100
40
10
15
20
Tiempo
25
30
35
40
Figura 4: Registros de la historia de la Velocidad en funcin del tiempo. a) Whittier Narrows, estacin Downey
Birchdale; b) Cape Mendocino, estacin Cape Mendocino; c) Northridge, estacin Newhall - Fire Sta;
d) Loma Prieta estacin Saratoga - W Valley Coll.
a)
10
b)
PGV=2.76
10
Pulso Indicador=0
Pulso Tardio=F
Pulso Indicador=0.88
Pulso Tardio=F
cm/s
cm/s
PGV=5.47
-5
-5
Ip=88.9
-10
10
15
20
Tiempo
25
30
Ip=31.6
35
-10
40
10
15
c)
10
25
10
Pulso Indicador=0
Pulso Tardio=F
35
40
PGV=5.87
Pulso Indicador=0.90
Pulso Tardio=F
cm/s
-5
-5
Ip=55.3
-10
30
d)
PGV=4.55
cm/s
20
Tiempo
10
15
20
Tiempo
25
30
Ip=30.3
35
40
-10
10
15
20
Tiempo
25
30
35
40
Figura 5: Registros de la historia de la Velocidad en funcin del tiempo. a) Irpinia, Italy estacin Bovino;
b) Coyote Lake estacin SJB Overpass, Bent 5 g.l.; c) San Fernando estacin Lake Hughes #9; d) Coyote Lake
estacin SJB Overpass, Bent 3 g.l.
861
a)
b)
80
80
PGV=30.44
Pulso Indicador=0.02
Pulso Tardio=F
Pulso Indicador=0.98
Pulso Tardio=V
60
40
40
20
20
cm/s
cm/s
60
PGV=57.81
-20
-20
-40
-40
-60
-60
Ip=38.3
-80
10
20
30
40
50
Tiempo
60
70
Ip=26.3
80
-80
90
10
20
30
40
50
Tiempo
60
70
80
90
c)
80
PGV=67.21
Pulso Indicador=1
Pulso Tardio=V
60
40
cm/s
20
0
-20
-40
-60
Ip=27.6
-80
10
20
30
40
50
Tiempo
60
70
80
90
Figura 6: Registros de la historia de la Velocidad en funcin del tiempo. a) Imperial Valley-06, estacin
Chihuahua; b) Chi-Chi Taiwan, estacin TCU117; c) Chi-Chi Taiwan, estacin WGK.
a)
b)
200
200
PGV=77.90
Pulso Indicador=1
Pulso Tardio=F
Pulso Indicador=1
Pulso Tardio=F
100
cm/s
cm/s
100
PGV=95.40
0
-100
0
-100
Ip=9.9
-200
10
20
30
40
50
Tiempo
c)
60
70
80
Ip=9.1
90
-200
100
200
10
20
30
40
50
Tiempo
d)
60
70
Pulso Indicador=1
Pulso Tardio=F
PGV=191.10
0
-100
10
20
30
40
50
Tiempo
e)
60
70
80
Ip=9.5
90
-200
100
200
10
20
30
40
50
Tiempo
f)
60
70
80
90
100
200
PGV=114.90
Pulso Indicador=1
Pulso Tardio=F
PGV=59.30
Pulso Indicador=1
Pulso Tardio=F
100
cm/s
100
cm/s
100
-100
Ip=8.6
0
-100
0
-100
Ip=7.1
-200
90
Pulso Indicador=1
Pulso Tardio=F
100
cm/s
100
-200
80
200
PGV=87.70
cm/s
10
20
30
40
50
Tiempo
60
70
80
Ip=10.1
90
100
-200
10
20
30
40
50
Tiempo
60
70
80
90
100
Figura 7: Registros de la historia de la Velocidad en funcin del tiempo. a) Imperial Valley-06, estacin El
Centro Array #4; b) Erzican Turkey, estacin Erzincan; c) Northridge estacin Newhall - W Pico Canyon Rd; d)
Chi-Chi Taiwan, estacin TCU068; e) Imperial Valley-06, estacin EC Meloland Overpass FF; f) Chi-Chi
Taiwan, estacin TCU076.
862
NOMBRE
DE
TERREMOTO
ESTACIN
PGV
(cm/s)
ANLISIS J. BAKER
Indicador
de Pulso
Pulso
Tardo
PULSO
Ip
Coyote Lake
5.87
0.90
FALSO
NO
30.3
Coyote Lake
5.47
0.88
FALSO
NO
31.6
Chi-Chi, Taiwan
TCU117
57.81
0.98
VERDADERO
NO
26.3
Chi-Chi, Taiwan
WGK
67.21
1.00
VERDADERO
NO
27.6
Imperial Valley-06
Chihuahua
30.44
0.02
FALSO
NO
38.3
San Fernando
Lake Hughes #9
4.55
0.00
FALSO
NO
55.3
Irpinia, Italy-01
Bovino
2.76
0.00
FALSO
NO
88.9
1.00
VERDADERO
NO
13.7
0.32
FALSO
NO
16.8
0.99
VERDADERO
NO
13.7
0.45
FALSO
NO
13.8
1.00
FALSO
SI
9.8
Whittier Narrows-01
Cape Mendocino
Northridge-01
Loma Prieta
Imperial Valley-06
Downey - Birchdale
Cape Mendocino
Newhall - Fire Sta
Saratoga - W Valley Coll.
El Centro Array #4
36.51
58.52
119.77
71.47
77.9
Erzican, Turkey
Erzincan
95.4
1.00
FALSO
SI
9.1
Northridge-01
87.7
1.00
FALSO
SI
8.6
Chi-Chi, Taiwan
TCU068
191.1
1.00
FALSO
SI
9.5
Imperial Valley-06
EC Meloland Overpass FF
114.9
1.00
FALSO
SI
7.1
1.00
FALSO
SI
10.1
Chi-Chi, Taiwan-03
TCU076
59.3
Tabla 4: Contrastacin de los valores de impulsividad entre el criterio de J. Baker y el Mtodo propuesto.
8 CONCLUSIONES
El mtodo propuesto a travs del ndice de Impulsividad Ip se muestra como sencillo y
eficaz para evaluar el carcter impulsivo de un registro ssmico a travs de la historia de
velocidades horizontales del suelo.
El rango de valores que puede tomar Ip permite una clasificacin del grado de
impulsividad a fin de considerar en los anlisis distintos niveles. Tambin permite identificar
registros No impulsivos.
El procedimiento propuesto presenta una buena correlacin con el criterio de J. Baker, en
los casos que ambos criterios difieren, se observa que el mtodo propuesto se corresponde con
la observacin visual.
REFERENCIAS
Akkar, S., Yazgan, U., and Glkan, P., Drift estimates in frame buildings subjected to
nearfault ground motions, Journal of Structural Engineering, ASCE 131(5), 10141024,
2005.
Alavi, B., and Krawinkler, H, Consideration of near-fault ground motion effects in seismic
design. Proceedings, 12th World Conf. on Earthquake Engineering., Auckland, New
Zealand, 2000.
Anderson, J. C., and Bertero, V. V., Uncertainties in establishing design arthquakes. J.
Struct. Eng., 113_8_, 17091724, 1987.
863
Baker, Jack W., Quantitative Classification of Near-Fault Ground Motions Using Wavelet
Analysis, Bulletin of the Seismological Society of America, Vol. 97, No. 5, pp. 14861501,
2007.
Baker, Jack W, Identification of near-fault velocity pulses and prediction of resulting response
spectra. Geotechnical Earthquake Engineering and Soil Dynamics IV, Sacramento,
California, 10, 2008.
Bray, J. and Rodriguez-Marek A., Characterization of forward-directivity ground motions in
the near-fault region. Soil Dynamics and Earthquake Engineering, 24, 815-828. Elsevier
Ltd, 2004.
Cox, K. E., and Ashford, S. A., Characterization of large velocity pulses for laboratory testing.
Pacific Earthquake Engineering Research Center, University of California at Berkeley,
Berkeley, California, 60, 2002.
FEMA, Quantification of Building Seismic Performance Factors, FEMA P-695 Report,
prepared by Applied Technology Council for Federal Emergency Management Agency,
Washington, D.C, 2009.
Fu, Q., and Menun, C., "Seismic-environment-based simulation of near-fault ground
motions." Proceedings, 13th World Conference on Earthquake Engineering, Vancouver,
Canada, 15, 2004.
Hall, J.F., Heaton, T.H., Halling, M.W. and Wald, D.J., Near-source ground motion and its
effects on flexible buildings. Earthquake Spectra 11:4, 569-605,1995.
Krawinkler, H., Alavi B., and Zareian F.(2005), Impact Of Near-Fault Pulses On Engineering
Design, Directions In Strong Motion Instrumentation, 83106.Springer.Printed In The
Netherlands, 2005.
Mavroeidis, G. and Papageorgiou, A., A mathematical representation of near-fault ground
motions.Bulletin of the Seismological Society of America 93:3, 1099-1131, 2003.
Makris, N., and Black, C. J., Evaluation of peak ground velocity as a good intensity
measure for near-source ground motions. J. Eng. Mech., 130_9_, 10321044, 2004.
PEER NGA (Pacific Earthquake Engineering Research Center).http://peer.berkeley.edu/smcat
Sasani, M. and Bertero, V., Importance of severe pulse-type ground motions in perfomancebased engineering: historical an critical review. 12th World Conference on Earthquake on
Engineering, New Zeland, paper 1302, 2000.
Shahi S. and Baker J.W., Signal Processing And Probabilistic Seismic Hazard Analysis Tools
For Characterizing The Impact Of Near-Fault Directivity, Joint Conference Proceedings
7th International Conference On Urban Earthquake Engineering (7cuee) & 5th
International Conference On Earthquake Engineering (5icee), Tokyo Institute Of
Technology, Tokyo, Japan, 2010.
Sofware MatLab. www.manualespdf.es/manual-matlab-6-1
Somerville, P. G., Magnitude scaling of the near fault rupture directivity pulse, Phys. Earth
Planet. Interiors 137, no. 1, 12, 2003.
AGRADECIMIENTOS
Programa de Becas Doctorales para Docentes UTN. Subsecretara de Posgrado. Secretara de
Ciencia y Tecnologa. Rectorado. UTN.
Universidad Tecnolgica Nacional Facultad Regional Mendoza.
CeReDeTeC. Centro Regional de Desarrollos Tecnolgicos para la Construccin, Sismologa
e Ingeniera Ssmica.