Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Bmfcis 218 e

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 118

Universidad Austral de Chile

Facultad de Ciencias de la Ingeniera


Escuela de Ingeniera Civil Electrnica

ESTUDIO DE LA ROBUSTEZ DE LA ESTIMACIN


ESPECTRAL DE LA TCNICA H/V DE NAKAMURA
PARA ESTUDIOS DE CARACTERIZACIN
DINMICA DE SUELOS
Tesis para optar al Ttulo de:
Ingeniero Civil Electrnico
Profesor Patrocinante:
Sr. Alfredo Illanes Manrquez
Ingeniero Ejecucin Electrnica
Master en Seales y Telecomunicaciones
Doctor en tratamiento de Seales y
Telecomunicaciones

VCTOR HUMBERTO SANDOVAL GMEZ


MARCO JAVIER VERA CATALN
VALDIVIA CHILE
2013

"Un viaje de mil kilometros debe empezar por un paso."


Lao-Ts

Agradecimientos
Vctor Sandoval Gmez.
En primer lugar agradecer a mis profesores guas en esta aventura Dr. Alfredo Illanes, Dr. Galo Valdebenito y Dr. Pedro Rey, ya que con su generosa e importante
ayuda cerrare un ciclo de mi vida. En el camino me ensearon muchos conocimientos tericos pero sin duda los que mas recordare son los humanos.
Mis mas grandes agradecimientos tambin a mis compaeros de universidad, profesores, amigos, familiares y tantos otros por su aliento constante en esta difcil labor,
y en especial a mis padres Humberto Sandoval Oyarzo y Olga Gmez Aichele quienes
con su esfuerzo y trabajo constante lograron hacer realidad este tan anhelado sueo.

Marco Vera Cataln.


Quisiera agradecer a mis profesores guas, Dr. Alfredo Illanes, Dr. Galo Valdebenito
y Dr. Pedro Rey, por entregar sus conocimientos y apoyo en cada momento del
desarrollo de este proyecto.
A mis padres por entregarme los valores necesarios para transitar por la vida, y en
especial a mi madre Lucy Cataln, que con su esfuerzo y apoyo a lo largo de estos
aos, es posible este logro en mi vida.
A mis compaeros y amigos de la universidad por hacer del pasar por esta etapa, un
recuerdo que permanecer en mi memoria.
A mis amigos de la vida por su constante apoyo y porque siempre han estado ah y
espero que siempre lo esten, adems por alegrarme el da a da.
A mi familia por su constante apoyo y preocupacin en los momentos ms dificiles de
este ciclo de mi vida, en especial a mis tios Fredy, Silvia, Ruben, jeanette, Guillermo
y Maria,.
Y a todos y cada uno de los que siempre han estado presentes apoyando y ayudando
de uno u otra forma en todos los momentos de este transitar por el camino de la
universidad.
Hay tantos a quien nombrar, solo me queda decirles GRACIAS A TODOS!. . .

ndice general
Agradecimientos

Resumen

1. Captulo I INTRODUCCIN
1.1. Antecedentes . . . . . . . . . .
1.2. Planteamiento del problema . .
1.2.1. Problematica General . .
1.2.2. Problematica Especfica
1.3. Objetivos . . . . . . . . . . . .
1.3.1. Objetivo General . . . .
1.3.2. Objetivos Especficos . .

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

2. Captulo II: FUNDAMENTOS Y CONCEPTOS


2.1. Riesgo ssmico . . . . . . . . . . . . . . . . . . . .
2.1.1. Peligrosidad ssmica . . . . . . . . . . . .
2.1.2. Vulnerabilidad ssmica . . . . . . . . . . .
2.1.3. Caracterizacin Dinmica de Suelos . . . .
2.2. Mtodo de Nakamura . . . . . . . . . . . . . . . .
2.2.1. Estimacin del espectro de Nakamura . . .
2.2.2. Resultados de la estimacin de espectros .
2.3. PDS en microvibraciones y comparacin con otras

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.

.
.
.
.
.
.
.

. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
reas de geologa.

3. Captulo III CONCEPTOS TERICOS DEL PROYECTO


3.1. Tcnicas tradicionales de estimacin del espectro de Nakamura.
3.1.1. Filtros . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.1.2. Ventaneo . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.1.3. Suavizado . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2. Tcnicas Modernas de Estimacin de Espectros . . . . . . . . .
3.2.1. Anlisis en Tiempo-Frecuencia . . . . . . . . . . . . . . .
4. Captulo IV MATERIALES Y MTODOS
4.1. Base de datos de las seales de ruido ambiental
4.2. Software Utilizados . . . . . . . . . . . . . . . .
4.2.1. MatLab . . . . . . . . . . . . . . . . . .
4.2.2. Geopsy . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.

8
8
11
11
12
13
13
14

.
.
.
.
.
.
.
.

15
15
15
16
18
23
25
27
31

.
.
.
.
.
.

34
34
37
40
45
47
48

.
.
.
.

54
54
55
55
55

ndice general
4.3. Protocolo de orden y procesamiento para cada seal de la base de datos 56
5. Captulo V RESULTADOS FINALES
5.1. Primera Situacin Observada . . . . . . . . . . . . . . . . . . . . .
5.2. Segunda Situacin Observada . . . . . . . . . . . . . . . . . . . . .
5.3. Tercera Situacin Observada . . . . . . . . . . . . . . . . . . . . . .
5.4. Estudio estadstico comparativo entre los resultados de Obras Civiles
y Geopsy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.5. Situaciones Observadas con Tcnicas Modernas . . . . . . . . . . .
5.5.1. Registro correspondiente a Pitrufqun 10 . . . . . . . . . . .
5.5.2. Registro correspondiente a Purn 5 . . . . . . . . . . . . . .
5.5.3. Registro correspondiente a Mininco 6 . . . . . . . . . . . . .
5.5.4. Registro correspondiente a Capitn Pastene 2 . . . . . . . .
5.6. Eliminacin de ruido transientes . . . . . . . . . . . . . . . . . . . .

57
. 63
. 64
. 66
.
.
.
.
.
.
.

67
72
73
75
78
80
84

6. Captulo VI Conclusiones

88

Bibliografa

91

A. Planes de trabajo y sectorizaciones de las zonas de estudios


A.1. Plan de trabajo y sectorizacin de la ciudad de Valdivia . .
A.1.1. Descripcin de equipos y parmetros de adquisicin
A.1.2. Descripcin de la campaa de campo . . . . . . . .
A.2. Protocolo de procesamiento de seales . . . . . . . . . . .
A.3. Sectorizacin de la localidad de Capitn Pastene . . . . . .
A.4. Sectorizacin de la localidad de Mininco . . . . . . . . . .
A.5. Sectorizacin de la localidad de Lumaco . . . . . . . . . .
A.6. Sectorizacin de la localidad de Purn . . . . . . . . . . .
A.7. Sectorizacin de la localidad de Pitrufqun . . . . . . . . .
B. Protocolo de Procesamiento

. .
.
. .
. .
. .
. .
. .
. .
. .

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.

98
98
99
101
102
103
103
104
104
105
106

C. Aclaracin de Conceptos
112
C.1. Licuefaccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112

ii

ndice de figuras
2.1.
2.2.
2.3.
2.4.
2.5.
2.6.
2.7.
2.8.
2.9.

Diferencias en el grado de Riesgo Sismico, dependiendo del lugar. .


Amplificacin de ondas sismicas segun el tipo de suelo. . . . . . . .
Mapa Isoperiodos Ciudad de Valdivia . . . . . . . . . . . . . . . . .
Ejemplo de relacin espectral horizontal y vertical en un emplazamiento en roca . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Caso 1 peak claro . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Caso 2: Peak poco claro de baja frecuencia . . . . . . . . . . . . . .
Caso 3: Dos peaks (f1>f0) . . . . . . . . . . . . . . . . . . . . . . .
Caso 4: peak ancho o mltiples peaks . . . . . . . . . . . . . . . . .
Caso5: peaks agudos y de origen industrial . . . . . . . . . . . . . .

3.1. Diagrama de bloque estimacin de espectro de Nakamura. . . . . . .


3.2. Seal de entrada de tres componentes EO, NS, y V. . . . . . . . .
3.3. Componentes de una seal de microvibracin con solape de ventanas
de 10 % y 30 %. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.4. Ventana deslizante, proceso de calculo del espectro FFT. . . . . . .
3.5. Espectro FFT promediado de cada componente de una seal. . . . .
3.6. Razn (H/V ) promediada. . . . . . . . . . . . . . . . . . . . . . . .
3.7. Filtro Butterworth de distintos ordenes . . . . . . . . . . . . . . . .
3.8. Operacin ventaneo . . . . . . . . . . . . . . . . . . . . . . . . . .
3.9. Representacin ventana Rectangular . . . . . . . . . . . . . . . . .
3.10. Espectro de la Ventana Rectangular . . . . . . . . . . . . . . . . . .
3.11. Filtro caracterstico de la ventana Rectangular . . . . . . . . . . .
3.12. Representacin Ventana de Hanning o Hann . . . . . . . . . . . .
3.13. Filtro caracterstico de la ventana de Hanning . . . . . . . . . . . .
3.14. Comparacin entre tcnica H/V estndar (a) y H/V incorporando la
funcin de suavizado de Konno y Ohmachi (b) . . . . . . . . . . . .
3.15. Espectro FFT y Suavizado del mismo. . . . . . . . . . . . . . . . .
3.16. Explicacin grfica de la STFT . . . . . . . . . . . . . . . . . . . .
3.17. STFT Aplicada a un registro de microvibracin . . . . . . . . . . .

. 18
. 21
. 23
.
.
.
.
.
.

25
28
28
29
29
30

. 35
. 35
.
.
.
.
.
.
.
.
.
.
.

36
36
36
37
40
42
43
43
44
44
45

.
.
.
.

46
47
49
50

4.1. Etapas Generales del procesamiento de datos . . . . . . . . . . . . . . 56


5.1. Distribucin porcentual de todas las seales procesadas en nuestro
estudio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62

ndice de figuras
5.2. Figura donde se aprecia que dos espectros de Geopsy estiman una frecuencia natural, mientras que las otras dos pruebas restantes estiman
otra frecuencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.3. Porcentaje de seales con la primera situacin . . . . . . . . . . . . .
5.4. Espectro de una prueba realizada en Geopsy no sigue la forma de los
otros espectros realizados con otros conjuntos de parmetros . . . . .
5.5. Porcentaje de seales que poseen un espectro diferente a las dems
pruebas de Geopsy . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.6. Truncamiento de los espectros de Geopsy en comparacin al peak ms
definido de eoc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.7. Porcentaje de seales con peak truncado en comparacin al peak mas
definido de eoc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.8. Histograma 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.9. Histograma 1 detallado (a). . . . . . . . . . . . . . . . . . . . . . . .
5.10. Histograma 1 detallado (b). . . . . . . . . . . . . . . . . . . . . . . .
5.11. Histograma 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.12. Histograma 2 detallado (a) . . . . . . . . . . . . . . . . . . . . . . .
5.13. Histograma 2 detallado (b) . . . . . . . . . . . . . . . . . . . . . . .
5.14. Histograma 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.15. Histograma 3 detallado (a) . . . . . . . . . . . . . . . . . . . . . . .
5.16. Histograma 3 detallado (b) . . . . . . . . . . . . . . . . . . . . . . .
5.17. Seal Pitrufquen 10 por componentes y en el dominio temporal . . . .
5.18. Imagen de los espectros de Pitrufqun 10 calculados con eoc y Geopsy
5.19. Modelamiento AR variante en el tiempo de Pitrufqun 10 . . . . . . .
5.20. mesh del espectro obtenido con modelamiento AR variante en el tiempo
5.21. Registro Purn 5 por componentes y en el dominio temporal . . . . .
5.22. Anlisis espectral de Purn 5 tanto por Geopsy como por eoc . . . . .
5.23. Modelamiento AR variante en el tiempo de Purn 5 . . . . . . . . . .
5.24. Nakamura AR variante en el tiempo del espectro correspondiente al
registro Purn 5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.25. Registro Mininco 6 por componentes y en el dominio temporal . . . .
5.26. Anlisis espectral de Mininco 6 con Geopsy y eoc . . . . . . . . . . .
5.27. Modelado AR variante en el tiempo del registro Mininco 6 . . . . . .
5.28. Capitn Pastene 2 por componentes y en el dominio temporal . . . .
5.29. Tiempo-Frecuencia de la componente E-O de Capitn Pastene 2 . . .
5.30. Tiempo- Frecuencia de la componente N-S de Capitn Pastene 2 . . .
5.31. Tiempo-Frecuencia de la componente Vertical de Capitn Pastene 2 .
5.32. Anlisis espectral de Pastene 2 con Geopsy y eoc . . . . . . . . . . .
5.33. Modelamiento AR variante en el tiempo de Capitn Pastene 2 . . . .
5.34. Seal en con sus tres componentes en el dominio temporal . . . . . .
5.35. Anlisis espectral sin anti triggering . . . . . . . . . . . . . . . . . . .
5.36. Estimaciones Frecuenciales para la segunda prueba con anti triggering
5.37. Comparacin de los (H/V)2 en las dos pruebas realizadas . . . . . . . .

63
64
65
66
66
67
68
68
69
69
70
70
71
71
72
73
74
74
75
76
76
77
78
79
79
80
81
81
82
82
83
83
84
85
86
87

ndice de figuras
A.1.
A.2.
A.3.
A.4.
A.5.
A.6.
A.7.
A.8.

Sectorizacin de la zona de estudio . . . . . . . .


Equipo sismolgico de medicin Geosig GBV 316
Equipo sismolgico de medicin Tromino ENGY .
Plan de trabajo de Capitn Pastene . . . . . . . .
Plan de trabajo de Mininco . . . . . . . . . . . . .
Plan de trabajo de Lumaco . . . . . . . . . . . . .
Plan de trabajo Purn . . . . . . . . . . . . . . .
Plan de trabajo de Pitrufqun . . . . . . . . . . .

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

98
100
100
103
103
104
104
105

ndice de cuadros
4.1. Frecuencias de muestreo . . . . . . . . . . . . . . . . . . . . . . . . . 54
5.1. Nomenclatura de las tablas que poseen el detalle de los registros de
microfizacin ssmicos . . . . . . . . . . . . . . . . . . . . . . . . .
5.2. Seales Procesadas de Pitrufqun . . . . . . . . . . . . . . . . . . .
5.3. Seales Procesadas de Capitn Pastene . . . . . . . . . . . . . . .
5.4. Seales Procesadas de Purn . . . . . . . . . . . . . . . . . . . . .
5.5. Seales Procesadas de Lumaco . . . . . . . . . . . . . . . . . . . .
5.6. Seales Procesadas de Mininco . . . . . . . . . . . . . . . . . . . .
5.7. Seales Procesadas de Valdivia . . . . . . . . . . . . . . . . . . . .
5.8. Seales Procesadas de Valdivia . . . . . . . . . . . . . . . . . . . .
5.9. Distribucin de seales procesadas en cada localidad . . . . . . . .
5.10. Lista de seales donde la frecuencia estimada por Geopsy coincide
con otro conjunto de parmetros de procesado . . . . . . . . . . . .
5.11. Resumen de seales que poseen un espectro diferente a las otras pruebas de Geopsy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.12. Lista de seales que presentan la condicin de peak truncado en comparacin al peak mas definido de eoc . . . . . . . . . . . . . . . . .
5.13. Frecuencias estimadas para la primera prueba, en donde no se utiliz
anti triggering . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.14. Frecuencias estimadas para la segunda prueba, en donde se utiliz anti
triggering . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.15. Corrimiento frecuencial en valor absoluto entre la prueba uno y la
prueba dos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.

57
58
59
59
59
59
60
61
62

. 64
. 65
. 67
. 85
. 86
. 86

A.1. Nmero de puntos de medicin por unidad de superficie en cada sector 99

Resumen
El riesgo ssmico es un concepto que considera la posibilidad de que se produzcan prdidas de vidas humanas y prdidas econmicas debidas a la accin de un
terremoto. Tiene dos componentes principales: la peligrosidad ssmica y la vulnerabilidad ssmica. La peligrosidad ssmica hace referencia a la probabilidad de que
una determinada accin ssmica se produzca sobre un punto del territorio con una
cierta extensin, intensidad y duracin. La vulnerabilidad ssmica hace referencia al
impacto de la accin ssmica sobre la poblacin, edificios, servicios, instalaciones,
etc. De este modo, un mismo sismo puede ser muy destructor si se produce en una
zona habitada y con construcciones de mala calidad, o bien causar pocos daos si
se produce en una zona desierta o con construcciones de buena calidad [ALI].
Los estudios de peligrosidad ssmica, informan de la importancia de la vibracin del
movimiento del suelo esperable en un lugar, debido a la accin ssmica. Combinando
este conocimiento con la evaluacin de la vulnerabilidad de los elementos expuestos,
se puede evaluar el riesgo ssmico, es decir caracterizar el grado de dao que puedan
sufrir los edificios y las infraestructuras debido a la accin ssmica. Las metodologas
utilizadas para el anlisis del riesgo ssmico ocupan una herramienta muy importante, como lo es el anlisis dinmico de suelos [LESA], ya que a travs de l podemos
conocer la respuesta o comportamiento de los suelos frente a las ondas ssmicas; esta
respuesta dinmica est ntimamente ligada con la frecuencia natural de vibracin
del suelo. El anlisis se realiza en un lugar atendiendo a las condiciones geolgicas
locales y su propsito es, entregar mayor informacin local sobre un terreno, ya que
las condiciones geolgicas a muy poca distancia puede causar importantes daos en
infraestructura entre otros.
La evaluacin de la repuesta local de sitio puede ser tratada por mtodos tericos
y empricos. Los primeros permiten un anlisis atendiendo a mltiples factores del
problema, aunque requieren una descripcin detallada de las capas del suelo que
conforman el estrato y de las fuentes que originan los sismos, las cuales en algunos
casos pueden ser muy complejas y por lo tanto, se ve dificultada la realizacin de un
modelo matemtico coherente. Los mtodos empricos se basan en registros de datos
ssmicos en el sitio de estudio, de tal manera que se puedan encontrar amplificaciones
relativas y perodos dominantes de oscilacin. Al contrario del mtodo anterior, no
se requiere una descripcin detallada de las capas del suelo, pero se necesita tener
una gran cantidad de registros o informacin ssmica, lo cual es posible en regiones
donde se tenga una alta sismicidad y se cuente con redes de registro ssmico. En el

caso de no contar con mucha informacin ssmica, se debe recurrir al uso de otros
mtodos. Entre las alternativas para obtener datos ssmicos de una regin en particular se encuentran las microtrepidaciones; estas son vibraciones de baja amplitud
del terreno y estn conformadas por ondas internas y superficiales [Ani]. El uso de
registros de microtrepidaciones para la evaluacin de la respuesta dinmica de suelos
es muy atractivo, debido a su bajo costo y a la facilidad y rapidez para obtener mucha informacin [LESA]. Para el anlisis de microtrepidaciones se han propuesto tres
tcnicas que son: Interpretacin directa de la Transformada de Fourier, tambin llamada Amplitudes Espectrales; clculo de los espectros relativos de amplificacin de
las estaciones en suelo blando versus una estacin en suelo firme o tcnica de Kagami
y clculo de las relaciones espectrales entre la componente horizontal y vertical del
movimiento en un mismo sitio o tcnica de Nakamura. En esta ltima tcnica est
enfocado el estudio que se presenta a continuacin, aunque se debe tener claro que
no se presentarn los resultados de algn estudio de Microzonificacin ssmica de un
lugar en particular, ms bien se trata del estudio de la robustez de la estimacin de
la frecuencia natural, al comparar los resultados obtenidos al utilizar dos softwares
de anlisis de microzonificacin que utilizan herramientas de tratamiento frecuencial
tradicionales. En nuestro trabajo se investigarn los resultados observables cuando
se varan los parmetros del procesado digital de las seales, como duracin de las
ventana, solape entre estas, eliminacin de transientes, adems de utilizar otras tcnicas ms recientes de estimacin espectral, aplicadas en otras reas cientficas, en
este caso consisti en el modelado paramtrico autorregresivo, comparando las principales frecuencias obtenidas con este mtodo, con las calculadas con metodologas
tradicionales.
Como se mencion anteriormente este estudio se basa en las estimaciones frecuenciales de la relacin espectral (H/V ), la cual requiere el clculo espectral de todos
los componentes grabados, que es slo una estimacin del espectro real. El valor
estimado puede variar con las condiciones de grabacin o con el mtodo utilizado
para la estimacin espectral. Incluso el espectro puede variar cuando los parmetros
de un algoritmo de estimacin espectral se modifican. El propsito principal de este
trabajo es realizar un anlisis de la robustez del clculo de la razn espectral (H/V )
utilizando mtodos de estimacin espectral variables en el tiempo. En una primera etapa, para cada grabacin de ruido ambiental, se calculan varios coeficientes
espectrales (H/V ), con dos diferentes programas informticos para el anlisis de microzonificacin. Cada relacin espectral se calcula usando un conjunto diferente de
parmetros del algoritmo y las frecuencias estimadas fundamentales de cada espectro se comparan. En un segundo paso se realiza un anlisis de estimacin espectral
variable en el tiempo para cada grabacin utilizando un modelo paramtrico autorregresivo. Las frecuencias principales obtenidas con este mtodo se comparan con
las calculadas mediante las metodologas tradicionales. Los resultados muestran que
en la mayora de las grabaciones de la tcnica tradicional se realiza una estimacin
correcta de la frecuencia fundamental. Sin embargo, en algunos casos, la frecuencia
obtenida no corresponde a la fundamental, sino a la frecuencia de los intervalos de

seal no estacionarias. Adems, en algunos otros casos espectros calculados con diferentes algoritmos o el mismo algoritmo utilizando diferentes parmetros puede dar
lugar a diferentes frecuencias fundamentales.

1. Captulo I INTRODUCCIN
1.1.

Antecedentes

Gran parte de los problemas que afectan a la infraestructura y a las personas, luego
de ocurrir un evento ssmico en alguna regin, es debido al desconocimiento de
nuestro entorno y de las condiciones del suelo en donde realizamos nuestra vida, es
ms, muchas veces hemos pasado una vida entera pensando que todo es esttico,
hasta que la naturaleza nos muestra otra realidad, y en esta falta de informacin se
encuentran temas como estudio de riesgo ssmico, mapas de isoperiodos, condiciones
geolgicas locales de nuestra localidad, entre otras. Y, es que, en el comn de la gente,
la pregunta Qu es el riesgo ssmico?, no ha tomado la relevancia que debera tener,
aun cuando est nos puede entregar informacin valiosa sobre reas geogrficas con
mayor o menor riesgo, proporcionar ayuda en el desarrollo de planes de respuesta
ante una emergencia ssmica, ayudar a la estimacin dentro de las primeras horas
luego de la ocurrencia de un terremoto del impacto que se podra haber causado. Las
condiciones o efectos locales de sitio son determinantes e influyentes en la hora de
analizar los daos originados por un sismo, ya que de unos kilmetros de distancia
de un lugar a otro, los cambios en los estratos de la tierra pueden ser diametralmente
opuestos, pasando desde un terreno muy blando a otro muy rgido, por ejemplo, esto
trae como consecuencia que la frecuencia con que vibra el terreno es distinta a cortas
distancias, o que puede haber una amplificacin de las ondas ssmicas debido a los
efectos locales del sitio afectado [ALI, SAM+ ].
La respuesta local del sitio puede ser evaluada de manera terica o emprica [Pil07,
LESA]. En ambos casos se hace compleja la tarea, ya que en los mtodos tericos
requieren una descripcin detallada de las capas que conforman el suelo y adems
las fuentes que generan los sismos, informacin que a veces es difcil de conseguir.
En los mtodos empricos se requiere un registro ssmico histrico del terreno a
analizar [Ani], informacin que si no est disponible es muy difcil de poder suplir
o se demora tiempo en recolectar a travs de redes de monitoreo. Existen otras
alternativas para poder obtener datos de forma emprica y aqu encontramos las microvibraciones o ruido ssmico ambiental, las cuales aprovechan el ruido ambiental
y/o artificial, correspondiente a vibraciones producidas por actividades humanas,
vibraciones causadas por fuentes superficiales locales como las industrias y el trfico vehicular, ruido natural de baja frecuencia generado por mareas y vientos. Las

1.1 Antecedentes
ventajas de aprovechar el ruido ambiental han seducido a muchos investigadores,
quienes las han utlilizados en muchas partes del mundo para realizar mapas de riesgo ssmico. Este tipo de mediciones no requiere que las redes de monitoreo estn
pendientes de algn evento ssmico durante das o meses, tampoco se necesita un
nivel de inversin cuantiosa a diferencia de las otras opciones.
Todos estos mtodos tratan de entregar una estimacin de la frecuencia natural de
oscilacin del terreno, pero las consideraciones tericas son distintas existiendo tres
tipos de tcnicas en la estimacin de la frecuencia fundamental para las microtrepidaciones o microvibraciones ssmicas. Estas tcnicas son: amplitudes espectrales,
tcnica de Kagami o espectros relativos a una estacin de referencia, tcnica de Nakamura o espectros horizontales entre verticales [LESA]. Precisamente esta ltima
es donde nos centraremos en el presente trabajo, estudiando la estimacin de los
espectros y analizando los cambios observables cuando al procesado de seales se
cambian parmetros que reiteradamente son utilizados en el rea de estudios de
microzonificacin ssmica. Es casi imposible poder continuar sin mencionar los resultados generados a partir del proyecto liderado por la unin europea SESAME,
el cual pretende dar directrices generales acerca del procesamiento y adquisicin
de seales de ruido ambiental. En trminos generales no es fcil encontrar muchos
estudios que se dediquen a investigar que sucede cuando los parmetros de procesamiento espectral son modificados, casi todos los toman de SESAME [Pro04], pero
no se preguntan si cuando se seleccionan otros parmetros, que perfectamente son
posibles de aplicar, se obtienen frecuencias distintas, as como tambin si con otras
tcnicas de estimacin espectral la determinacin de frecuencias cambia. Pero dentro
de lo poco que se encontr tenemos el trabajo del profesor Raul Ritta, el cual en
su tesis de magister estudio el tema, para un estudio de microfizacin ssmica de la
ciudad de Mayaguez, Puerto Rico [JEA12], en este trabajo se pueden leer interesantes pruebas no solo de estimacin de espectros, adems de comparacin con otras
tcnicas alternativas a la de Nakamura.
Otra lnea que se sigui en el estudio que se realiz se preocup de realizar una
comparacin sobre la aplicacin del procesamiento de seales. En el cual en el rea
de microvibraciones an se utilizan las tcnicas de estimacin de espectros basadas
en la Transformada de Fourier calculada con la clsica Transformada Rpida de
Fourier ms conocida como FFT. Esto se contrapone al avance en otras disciplinas
relacionadas como la sismologa y la vulcanologa, es de especial cuidado este hecho
ya que se podra estar afectando los estudios de riesgo ssmico que se realizan hoy
en da en diversas ciudades del mundo utilizando la tcnica de Nakamura. Todo esto
sin nombrar a los avances en otras reas donde el anlisis espectral est bastante
avanzado como por ejemplo la defensa civil, las telecomunicaciones y la ingeniera
biomdica. Estas reas de aplicacin utilizan tcnicas avanzadas de anlisis espectral
principalmente para anlisis de datos de radares, sonares, redes de sensores inalmbricos, codificacin de datos, bioseales [Sem04, CJ07, SL05] (electrocardiograma,
electroencefalograma, electromiograma, etc), entre otras muchas aplicaciones ms.

1.1 Antecedentes
Dentro de la forma de estimacin de la frecuencia natural tambin se debe tener en
cuenta la forma en que se realiza, ya que no es de una forma usual dentro del abanico
de tcnicas utilizadas en el tratamiento frecuencial. Esto se debe a la informacin
que tratan de extraer desde los registros de ruido ambiental, adems de tratar de
enfrentar el problema de los ruidos impulsionales, los cuales son de corta duracin
pero que tienen una gran magnitud. Estos factores se consideran en el procesado de
seales, pero no existe un estudio acabado si la labor se realiza de manera exitosa o es
realizado con una baja asertividad, actualmente lo que se realiza es hacer mediciones
de ruido ambiental que van desde los 15 a los 20 minutos, luego la persona que
procesa la informacin debe a mano sacar estos ruidos y dejar la seal con la menor
cantidad de estos eventos para proceder a estimar el espectro. Todos estos pasos
hacen lento el proceso de estimacin, adems de carecer de un criterio ampliamente
estudiado en cuanto a la estimacin frecuencial, y nico, ya que es el usuario el
que debe seleccionar con los criterios que el considere como ptimos. Por otro lado
no existen estudios que hablen del error en la estimacin espectral cuando la seal
es analizada con y sin transientes, los nicos intentos por tratar de automatizar
este proceso son los llamados algoritmos anti-trigger los cuales solo se encargan de
escoger un umbral (una vez ms definido arbitrariamente por un usuario o segn
alguna pauta general de procesado) que trata de sacar las parte ms conflictivas de
las seales para estimar el espectro.
Las nuevas tcnicas ocupadas en otras reas que apuntan al anlisis frecuencial de
seales no estacionarias o que presentan dinmicas transientes se basan en Transformada Wavelet y en modelado paramtrico. La primera tcnica es muy til para el
anlisis de transientes, mientras que la segunda para la estimacin espectral variante
en el tiempo (para seales no estacionarias) con una baja cantidad de datos. Estas
tcnicas modernas de anlisis espectral permiten una mayor robustez en la estimacin del espectro utilizando una baja cantidad de observaciones del proceso. Ms
an las tcnicas propuestas son generalmente inmunes a ruidos transientes, lo que
genera espectros no dependientes de las no estacionalidades presentes en los datos.
Con los antecedentes vistos anteriormente, la principal problemtica que pretende
abordar este trabajo, es lo que respecta a la estimacin de los espectros que son
utilizados en el mtodo de la razn espectral H/V de Nakamura, no as cuestionar
el mtodo, sino analizar que tan robusta es la estimacin de los espectros y por ende
la determinacin de la frecuencia natural de un sitio. Para ello se realizaron dos
etapas de anlisis, en la primera se compararon las estimaciones realizadas por dos
software distintos, que utilizan tcnicas tradicionales de estimacin de espectros. En
la segunda etapa se utilizo un mtodo de estimacin espectral variable en el tiempo,
en este caso el modelado paramtrico autorregresivo, estos resultados se compararon
con los mtodos tradicionales.

10

1.2 Planteamiento del problema

1.2.

Planteamiento del problema

1.2.1.

Problematica General

El anlisis de la respuesta del sitio, es una parte fundamental de la evaluacin del


riesgo ssmico en reas propensas a terremotos [Pro04]. La naturaleza y la distribucin del dao ssmico estn muy influenciadas por la respuesta del suelo frente a
cargas cclicas. Dicha respuesta est controlada en gran medida por las propiedades
del suelo. Por ejemplo, la susceptibilidad de un depsito de suelo para amplificar
el movimiento depende de las propiedades dinmicas, siendo el mdulo de corte y
el amortiguamiento principales parmetros para modelar la respuesta ssmica del
suelo. Estos parmetros se relacionan con otras propiedades del suelo, de carcter
ms fsico, como por ejemplo el ndice de poros, la densidad relativa y la velocidad
de las ondas de corte [Rod05]. Los mtodos utilizados para caracterizar la respuesta del suelo ante un evento ssmico pueden ser clasificados en mtodos numricos y
mtodos experimentales. Los mtodos de simulacin numrica son modelos matemticos utilizados en sitios donde no existen los requerimientos necesarios para utilizar
mtodos experimentales [Alf]. Entre los mtodos tericos y numricos se encuentran
los relacionados con la resolucin analtica de las ecuaciones que describen el problema, y los mtodos de simulacin terica y numrica. Constituyen la denominada
aproximacin terica. En la aplicacin de estos mtodos se parte de un modelo
que se aproxime a la situacin real. Los mtodos experimentales tienen por objetivo
determinar una funcin de transferencia de los depsitos superficiales, a partir de registros de sismos reales o mediciones de ruidos ambientales. Entre los ms aceptados
en las ltimas dcadas, se encuentran los basados en microvibraciones ambientales
[Val12]. Los microtremores, tambin conocidos como ruidos ambientales ssmicos,
son vibraciones aleatorias inducidas en las masas de suelo y roca por fuentes de
origen naturales y artificiales [Con].
Las tcnicas experimentales son herramientas muy tiles, puesto que en algunos
casos encontrar un emplazamiento que sirva de referencia no es fcil. Entre ellas cabe
destacar el mtodo de las razones espectrales (H/V ) o HVSR. Este consiste en hallar
las razones espectrales suavizadas entre las componentes horizontales y verticales
correspondientes a un mismo emplazamiento. Otro muy parecido a ste ltimo es
el conocido por mtodo de Nakamura. La diferencia entre estas dos tcnicas radica
en el tipo de registro utilizado: en el primer caso son los registros de aceleracin,
velocidad o desplazamientos sobre los que se hallan las razones, sin embargo, la
tcnica de Nakamura se aplica sobre registros de ruido ambiental [Rub].
Nakamura propuso en el ao 1989 [Nak89] que el cociente entre la componente
horizontal y la componente vertical en el dominio frecuencial era un buen estimador
de la respuesta de sitio. Esta tcnica se aplic por primera vez a microtemblores de
zonas urbanas de Japn.

11

1.2 Planteamiento del problema


Nakamura supone que la razn (H/V ) o razn espectral corresponde a la funcin de
transferencia de los niveles superficiales sometidos a movimientos horizontales. Basa
su trabajo en la hiptesis de que un microtemblor est compuesto por un cuerpo
principal de ondas y por ondas superficiales RayLeigh, inducidas por fuentes artificiales (trenes, trfico, etc). Estas ondas superficiales se consideran perturbaciones
que es necesario eliminar.
En la aplicacin de las tcnicas de clculo de razones espectrales existen algunas
dificultades. La eleccin del emplazamiento de referencia en funcin de los requisitos
que debe cumplir mencionados anteriormente, es determinante en los resultados
finales. A las tcnicas de suavizado tambin se les debe prestar atencin, ya que
una incorrecta aplicacin de stas puede ocultar rasgos importantes de las razones
[yEC08]. Tambin existe gran controversia respecto a la ventana del registro elegida;
si sta no es adecuada, puede que las amplitudes estimadas se contaminen con fases
que estn fuera de los objetivos del trabajo [JEA12].

1.2.2.

Problematica Especfica

Dentro de los estudios de caracterizacin de suelos, se pueden encontrar diversos


algoritmos que son utilizados por los fabricantes de los sensores industriales para
encontrar el periodo fundamental de vibracin en los terrenos, mediante el empleo
del mtodo de la razn espectral (H/V ) en superficie, conocida comnmente como
mtodo de Nakamura. Los datos entregados por estos sistemas de medicin son ampliamente utilizados por gelogos e ingenieros para la elaboracin de mapas ssmicos.
Sin embargo estos algoritmos vienen protegidos dentro de un software propietario y
los usuarios poco conocen el funcionamiento de estos algoritmos. Es decir los usuarios utilizan estos algoritmos como cajas negras teniendo acceso nicamente a la
modificacin de parmetros pero sin tener un conocimiento acabado de la influencia
de stos en el espectro finalmente calculado.
El mtodo de Nakamura requiere como entrada el clculo de espectros por componentes del registro de microvibracin. En general los algoritmos de clculo de la
razn espectral poseen los siguientes pasos [JEA12, eA08]:
Segmentacin de cada componente del registro en intervalos (con o sin solape).
Clculo de la FFT en cada ventana de cada componente.
Suavizado de la FFT en cada ventana de cada componente (se usa normalmente filtros suavizantes como triangulares o el de Konno Ohmachi).
Clculo de la razn espectral (H/V ) a partir de el espectro de las tres componentes en cada ventana

12

1.3 Objetivos
Promediacin de los espectro (H/V )
Hasta hoy en da no se ha realizado ningn estudio acerca de la robustez de la
metodologa para calcular la razn espectral ni el efecto que la modificacin de
parmetros podra provocar en el clculo del espectro final, y por ende en el clculo
de la frecuencia natural del suelo. Es sabido en el mbito de procesamiento de
seales que el clculo de espectros por medio de FFT promediado no es robusto
cuando se aplica a seales no estacionarias como las de microvibraciones. Adems la
no robustez o confiabilidad de las estimaciones espectrales que requiere Nakamura
pueden afectar el clculo de la frecuencia natural y por ende afectar el resultado de
los estudios y mapeos de caracterizacin de riesgo ssmico.
Entre las principales directrices, para la adquisicin de datos y estimacin de espectros se encuentran las entregadas por el proyecto de la unin europea SESAME
(Site effects assessment using ambient excitations) [Pro04], en el cual participaron
ms de 80 investigadores de 14 centros geofsicos, y que pretende controlar mejor
los fundamentos tericos y metodolgicos de la utilizacin del ruido ambiental ssmico, para estimar de forma vlida el efecto de sitio. A partir de lo anterior se han
realizados muchos estudios sobre microzonificacin ssmica en distintas ciudades del
mundo [yEC08], En los cuales los autores comentan la eleccin de parmetros, para
tratar las seales de acuerdo a los lineamientos generales de SESAME. Sin embargo
existen muy pocos estudios que hablen sobre los cambios que se observan en las
estimaciones de los periodos, cuando parmetros como longitud de ventana, solape,
eliminacin de transitorios; son variados por otros que no formen parte de ese estndar de valores, que trae por defecto el software que se gener como resultado de
SESAME.

1.3.

Objetivos

1.3.1.

Objetivo General

El propsito principal de este trabajo consiste en el anlisis de la robustez


de la estimacin de los espectros, utilizados en el mtodo de la razn espectral (H/V ) de Nakamura, con el fin de comparar la robustez en la estimacin
de la frecuencia natural de los suelos, obtenida con tcnicas tradicionales de
procesamiento de seales y la obtenida con tcnicas modernas.

13

1.3 Objetivos

1.3.2.

Objetivos Especficos

Estudiar el mtodo (H/V ) para determinar el riesgo ssmico.


Interpretar los resultados entregados por un software de estimacin espectral
utilizando la metodologa de Nakamura.
Utilizar tcnicas avanzadas de estimacin de espectros para comprobar la estimacin de los espectros de Nakamura.
Realizar un estudio de sensibilidad de parmetros utilizados en la estimacin
de espectros del mtodo Nakamura.
Concluir si la frecuencia natural fue bien o mal determinada en los registros
analizados.

14

2. Captulo II: FUNDAMENTOS Y


CONCEPTOS
2.1.

Riesgo ssmico

Riesgo ssmico se define como el dao potencial que sobre la poblacin y sus bienes,
la infraestructura, el ambiente y la economa pblica y privada, pueda causarse por
la ocurrencia de un evento ssmico. Los aspectos que contempla el estudio del riesgo
ssmico corresponden a factores como la probabilidad de que se produzca un evento
ssmico, efectos locales debido a amplificacin de las ondas ssmicas, fragilidad de las
edificaciones, alta densidad de habitantes o bienes que puedan resultar con daos
[Rsi].
La Peligrosidad y la Vulnerabilidad Ssmica sern definidos con mayor profundidad
en la Subseccin 2.1.1 y la Subseccin 2.1.2.

2.1.1.

Peligrosidad ssmica

La peligrosidad ssmica es la probabilidad de que ocurra un fenmeno fsico como


consecuencia de un terremoto, como pueden ser el movimiento mismo del terreno,
as como la licuefaccin, los deslizamientos de tierra, inundaciones, ruptura de fallas,
etc.
El tamao y localizacin de estos efectos dependen de diversos factores, principalmente de las caractersticas geolgicas y geotcnicas del lugar, pero indudablemente
de las caractersticas del terremoto [Raf07].
Por tal motivo, el primer paso en la evaluacin de la peligrosidad ssmica es caracterizar las zonas sismotectnicas, para posteriormente entender mejor las caractersticas de los terremotos. Generalmente, en su evaluacin se utilizan mtodos o modelos
probabilsticos simplificados de clculo basados en el establecimiento de leyes estadsticas para definir el comportamiento ssmico de una zona, las fuentes ssmicas y la
atenuacin del movimiento del suelo, expresando los resultados en forma de probabilidad de ocurrencia de los distintos tamaos de los terremotos, la probabilidad de

15

2.1 Riesgo ssmico


excedencia de distintos niveles de intensidad del movimiento o a los valores mximos
de aceleracin esperados en un lugar y en un intervalo de tiempo determinado. Sin
embargo, estos modelos involucran una gran cantidad de incertidumbres lo que lleva
inevitablemente a ser calculados a partir de la extrapolacin de datos, a la adaptacin de estudios de otras regiones para que estos modelos sean completamente
funcionales y en muchos casos a la simplificacin de los mismos [Som00].

2.1.2.

Vulnerabilidad ssmica

La otra parte importante para evaluar el riesgo ssmico de una regin es la evaluacin
de la vulnerabilidad ssmica de las estructuras construidas por el ser humano. La
experiencia ha enseado, a travs de los terremotos pasados que existen estructuras
de una misma tipologa que pueden sufrir un mayor grado de dao debido a un
terremoto a pesar de localizarse en el mismo sitio [Yep96, oSR89, DGA88]. Esto se
debe a que existen estructuras con una calidad estructural mejor que otras, o en
otras palabras, su vulnerabilidad es menor. Por lo tanto, se puede llegar a plantear
que la vulnerabilidad ssmica de una estructura o grupo de estructuras, es la calidad
estructural o capacidad de sus elementos estructurales para resistir un terremoto.
El que una estructura sea ms o menos vulnerable ante un terremoto de determinadas caractersticas, es una propiedad intrnseca de cada estructura, es decir, es
independiente de la peligrosidad ssmica del sitio de emplazamiento, por lo tanto,
una estructura puede ser vulnerable pero no estar en riesgo, a menos que se encuentre
en un sitio con una cierta peligrosidad ssmica. Puede observarse, desde este punto
de vista que los estudios de vulnerabilidad ssmica se pueden aplicar a cualquier
obra de ingeniera civil, como son edificaciones, presas, carreteras, puentes, taludes,
depsitos, centrales nucleares y, en general, a toda obra en la que se requiera conocer su comportamiento ante un posible terremoto y las consecuencias que puedan
producir.
Una de las primeras tareas que hay que plantearse en la evaluacin de la vulnerabilidad ssmica de las estructuras, es el alcance que tendrn los estudios de Riesgo
Ssmico. Este paso es importante, ya que la vulnerabilidad ssmica se puede evaluar
de una manera muy especfica realizando un estudio muy exhaustivo de las estructuras o muy general como puede ser un estudio a nivel urbano [MH]. El conocimiento
del comportamiento de una estructura es generalmente complejo y dependiente de
un sinnmero de parmetros que en su mayora son difciles de obtener. Algunos de
estos parmetros incluyen, la obtencin de las caractersticas del movimiento ssmico del lugar, la resistencia de los materiales con los que est construida, la calidad
de la construccin, la interaccin con los elementos no estructurales, el contenido
existente en la estructura en el momento del terremoto, entre otros muchos.
Sin embargo, la realizacin de estudios a nivel urbano, se puede reducir al conocimiento de algunos parmetros bsicos para poder clasificar la estructura, o, en otras

16

2.1 Riesgo ssmico


palabras al conocimiento de su calidad estructural. Normalmente, al plantear la realizacin de un estudio de riesgo ssmico lleva implcito la realizacin del estudio de
grandes reas para lo cual, los estudios a nivel urbano son los ms factibles. Por lo
tanto, es oportuno aclarar que estos estudios generalmente estn englobados dentro
de un marco estadstico, por lo que casi todas las estimaciones de la vulnerabilidad
tienen niveles significativos de incertidumbres asociados [Som00]. Gran parte de esto
se debe a que las estimaciones se realizan sin considerar estudios detallados en la
construccin, condicin y comportamiento de la estructura. A menudo, la estimacin de la vulnerabilidad se realiza basada en observaciones visuales sin referencia
en clculos del comportamiento estructural.
Un aspecto importante en la evaluacin de la vulnerabilidad ssmica de las estructuras es definir el dao que pueda sufrir una estructura debido a un terremoto. En
realidad, el dao va asociado a la vulnerabilidad, ya que una estructura es ms o
menos vulnerable dependiendo del dao que pueda sufrir ante un terremoto. La palabra dao se utiliza ampliamente para describir distintos fenmenos que puedan
ocasionar los movimientos ssmicos a las estructuras, refirindose principalmente el
deterioro fsico [R.03].
Una vez revisados los conceptos de Peligrosidad Ssmica y Vulnerabilidad Ssmica,
se puede observar que el Riesgo Ssmico depende directamente de ellos, es decir,
para que exista verdaderamente riesgo ssmico en un lugar, ambos conceptos deben
producirse y existir respectivamente, es decir, el riesgo ssmico evala y cuantifica
las consecuencias sociales y econmicas potenciales provocadas por un terremoto,
como resultado de la falla de las estructuras cuya capacidad resistente fue excedida.
El riesgo ssmico depende fuertemente de la cantidad y tipo de asentamientos humanos del lugar, tal y como se puede observar en la Figura 2.1. Para ejemplificar
grficamente esto, en la parte derecha de la figura observamos que, aunque el peligro potencial ssmico es muy alto en un determinado lugar, el riesgo ssmico es
pequeo porque es una regin muy deshabitada. En cambio, si observamos el lado
izquierdo de la figura el peligro ssmico no es tan grande porque los grandes sismos
no suelen ser tan frecuentes, pero la cantidad de personas que viven all y el tipo de
construccin, hacen que el riesgo ssmico sea muy grande; esto es claro al comparar
el nmero de vctimas, estructuras daadas o perdidas monetarias, para uno y otro
lugar.

17

2.1 Riesgo ssmico

Figura 2.1.: Diferencias en el grado de Riesgo Sismico, dependiendo del lugar.

2.1.3.

Caracterizacin Dinmica de Suelos

2.1.3.1.

Microzonificacin Ssmica y efectos de sitio

La microzonificacin ssmica consiste en dividir una regin en porciones ms pequeas, cada una de ellas, con un comportamiento similar durante un sismo, de
forma que puedan definirse recomendaciones precisas para el diseo y construccin
de edificios sismo resistentes. Para cada una de estas zonas, adems de especificarse
la fuerza sismica posible, se deben identificar los tipos de fenomenos que se puede
desencadenar a raiz del sismo, como son deslizamientos, amplificacion exagerada del
movimiento o la posibilidad de licuacion del suelo. Las leyes de atenuacin reflejan
la propagacin de las ondas ssmicas en la roca de la corteza. El movimiento en
la superficie del terreno en un sitio dado puede diferir totalmente del que se tiene
en la roca base, debido a alteraciones de las ondas causadas por efectos geolgicos,
topogrficos y de rigidez del suelo [Oba]. A estas alteraciones, llamadas en trminos generales efectos locales o efectos de sitio, se les ha dado mayor importancia
y se han considerado cada vez ms en los aos recientes y han conducido a la necesidad de realizar estudios de microzonificacin en las reas donde se encuentran
los asentamientos humanos, para detectar aquellas zonas que presentan problemas
especiales.
La microzonificacin ssmica es un procedimiento para la estimacin del riesgo ssmico a partir del movimiento del suelo y algunos fenmenos relacionados, tomando
en cuenta los efectos de las condiciones locales del sitio. El subsuelo y las condiciones topogrficas pueden amplificar o reducir el umbral de aceleracin del suelo

18

2.1 Riesgo ssmico


en un sitio con respecto a lo que podra esperarse para un sustrato consolidado en
ese mismo punto [Efr03]. Estos efectos locales de sitio son la base de lo que sera
un mapa de microzonificacin ssmica, donde el rea de estudio queda dividida en
sectores de diferente riesgo potencial, tambin llamados zonas de riesgo ssmico, que
se etiquetan, de acuerdo con su peligrosidad, en orden ascendente. Con esto se obtiene la informacin necesaria para deducir los efectos de un terremoto sobre las
edificaciones, en una ciudad. La informacin ms importante de este proceso es la
obstencion de la frecuencia natural de vibracion del suelo de cada una de estas zonas,
informacion que para nuestro caso correspondia a un estudio de microzonificacion
en diversos sectores de la ciudad de Valdivia, de los cuales se tomaron algunos de
estos registros para ser procesados.
La microzonificacin ssmica se divide en cuatro fases para su estudio [Pil07]. La
primera fase consiste en definir las zonas sismogeneticas (fuentes de terremotos con
caractersticas ssmicas y tectnicas homogneas), a partir del estudio de la sismicidad.
La segunda fase se refiere al desarrollo de diversos mapas de la zona, como: geolgico,
morfolgico, geotcnico de daos por sismos y otros; al combinar todos estos mapas
se obtienen las microzonas.
La tercera fase es el monitoreo ssmico utilizando diferentes tipos de registros, para
determinar movimientos fuertes, dbiles y la vibracin ambiental en las microzonas
estimadas en la segunda fase. Con esto se pueden conocer las distintas respuestas
ssmicas que tendrn las diferentes microzonas (efectos de sitio).
Una de las tcnicas utilizadas frecuentemente para estimar los efectos de sitio, en
regiones de moderada a alta sismicidad, es la razn espectral de registros de temblores, obtenidos de forma simultnea en una estacin en suelo blando y otra como
referencia en suelo ms firme. El problema de este tipo de mediciones, es que toma
mucho tiempo de espera, ya que los temblores fuertes, mayores a 7 grados, ocurren
con menor frecuencia que los temblores dbiles. Una alternativa es utilizar razones
espectrales de temblores dbiles. Esta tcnica tiene la limitacin, en zonas urbanas,
donde debido al alto nivel de ruido ambiental, principalmente por el trfico vehicular, provoca que la obstencin de registros tiles se realice en horas donde el ruido
disminuye, lo cual aumenta el tiempo de trabajo de campo.
Otra tcnica, la cual se ha hecho muy popular en los ltimos tiempos, es la estimacin de la razn espectral entre las componentes horizontales y la vertical de un
mismo registro, propuesta por Nakamura para interpretar mediciones de vibracin
ambiental. Los resultados de estudios de mediciones realizadas por Lermo y Chaves
Garcia [LCG94], sugieren que, si el efecto de sitio es causado por una geologa simple,
entonces Nakamura proporciona una estimacin gruesa de la frecuencia y la amplitud del primer modo de resonancia (ondas S). Una de las ventajas ms importantes
de esta tcnica es que elimina la necesidad de tener una estacin de referencia.

19

2.1 Riesgo ssmico


Para esta estimacin de la frecuencia y amplitud que realiza Nakamura, es necesario el uso del procesamiento digital de seales, utilizando diferentes tipos de filtros
y tcnicas de estimacin de espectros, es a esto ltimo hacia donde se enfoca este
trabajo, realizando un anlisis mas critico a la tcnica de estimacin actual, proponiendo una tcnica mas avanzada y as obtener una mejor estimacin del efecto de
sitio.
El uso de vibracin ambiental para la estimacin del efecto de sitio, fue introducido
en Japn alrededor de los aos 50. Algunos investigadores han cuestionado la utilizacin de vibraciones ambientales [KO98, CB94], aun as las mediciones de estas
vibraciones presentan un gran atractivo para la caracterizacin de la respuesta de
sitio, tanto por ser sencilla de operar y por su bajo costo, como tambin por la
rapidez con que permiten obtener resultados. Existen tres tcnicas, que son las ms
utilizadas para analizar los registros de vibracin ambiental: la determinacin de la
frecuencia mxima, a partir de la densidad espectral de potencia; el clculo de la
razn espectral y la tcnica de Nakamura. Estas tres tcnicas fueron analizadas por
Lermo y Chavez Garcia [LCG94], quienes realizaron comparaciones de resultados
de vibraciones ambientales, con la razn espectral estndar de registros de temblores. Su conclusin es que la vibracin, cuando se analiza con Nakamura, permite la
determinacin, con buena precisin del periodo dominante de sedimentos sujetos a
amplificacin dinmica en un intervalo de frecuencias entre 0.3 y 5 Hz, junto con
una estimacin preliminar del nivel de amplificacin.
La cuarta fase consiste en la clasificacin de las microzonas por orden de peligrosidad y seleccionar puntos idneos para la construccin, as como establecer algunas
recomendaciones para el uso de estos mapas.
La microzonificacin ssmica posee la finalidad de evaluar los riesgos ssmicos minuciosamente para lograr una aceptable planificacin regional y urbanstica y para
ponderar los niveles de las fuerzas o movimientos para la construccin de estructuras
antissmicas [Efr03].
Se ha observado que las intensidades de los sismos varan considerablemente a distancias muy cortas [MH]. Esto lleva a pensar que el factor esencial para la evaluacin
del dao de las estructuras son las condiciones del subsuelo local, es decir, los efectos
de sitio.

2.1.3.2.

Efectos de sitio

A la transformacin que sufren las ondas ssmicas desde la roca base, hasta la superficie, se le conoce como efecto de sitio y es diferente e independiente del efecto
de interaccin suelo-estructura, ya que se presenta est o no una estructura en la
superficie.

20

2.1 Riesgo ssmico


Las condiciones geolgicas y de suelos subyacentes a la superficie del terreno, pueden
influir de distintas maneras en la respuesta de un lugar a un sismo. En ciertos casos,
los desplazamientos permanentes del terreno pueden causar la deformacin de las
estructuras que se apoyan en ellos.
Las amplitudes y formas de las ondas ssmicas dependen del mecanismo focal y
de la cantidad de energa liberada, as como de la trayectoria y de las condiciones
locales. El mecanismo focal define la manera en que la energa es irradiada en el
espacio y el tiempo [Pil07]. La distribucin de la intensidad ssmica est ligada con
las caractersticas de la ruptura. Por otra parte, las ondas ssmicas, una vez emitidas
por la fuente, dependen fundamentalmente de las propiedades de los medios en que
se propagan y de las irregularidades que encuentra en su trayecto. Esto es aplicable
al describir los efectos de las condiciones locales. La Figura 2.2 muestra el cambio
de la seal al variar el tipo de suelo por el cual se propaga la onda. La amplificacin
de la frecuencia es ms elevada cuando se trata de terrenos blandos, mientras que
en terrenos ms rocosos la amplificacin disminuye.

Figura 2.2.: Amplificacin de ondas sismicas segun el tipo de suelo.

El problema fundamental que debe ser resuelto para la estimacin del riesgo ssmico
es evaluar la respuesta dinmica especfica de los depsitos superficiales y es comnmente referido a la estimacin del factor de amplificacin de dichos depsitos. Los
objetivos bsicos en el estudio de la respuesta de sitio son [Pil07]:
- Obtencin del Perodo Fundamental.
- Estimacin de la amplificacin del movimiento del terreno.
- Evaluar el espectro de respuesta para el diseo y la seguridad de las estructuras.
- Evaluacin del potencial de licuefaccin, (definicin del concepto en Seccin C.1).

21

2.1 Riesgo ssmico


La respuesta del sitio puede variar dependiendo de las condiciones del subsuelo y
las caractersticas geolgicas locales. Algunos parmetros importantes, que influencian esta respuesta, son la extensin horizontal y la profundidad de los depsitos
superficiales que suprayacen al basamento, topografa del basamento, sedimentos
depositados y fallas presentes en las zonas [FKG03].

2.1.3.3.

Mapa de isoperiodos

Como parte importante de una microzonificacin se contempla la elaboracin de


mapas donde, mediante curvas que unen valores iguales, se exprese la distribucin
de periodos dominantes del terreno. Comnmente, los mapas de isoperiodos guardan
una cierta correlacin con aquellos de geologa superficial, mostrando periodos cortos
(p. ej. 0.1 0.2 segundos) en terrenos firmes, mientras que en rellenos aluviales
llegan a observarse periodos de 3, 4 5 segundos [Pro04, yEC08, JEA12, VPCB05,
FAM+ 10].
Este tipo de mapas puede integrarse colectando, mediante un acelergrafo o sismgrafo porttil, seales de vibracin ambiental en tantos sitios como sea posible y
calculando sus espectros de amplitud de Fourier para obtener los periodos dominantes.
Aunque un mapa de periodos configurado con el solo anlisis de microtemblores
puede proporcionar una idea ms o menos aproximada de la respuesta dinmica,
se recomienda especialmente que, al menos en aquellos sitios donde se ubicaron
acelergrafos o sismgrafos, se verifiquen los periodos derivados del anlisis de la
vibracin ambiental con el propsito de minimizar las incertidumbres en anlisis
posteriores [UGU].
El conocer los periodos dominantes del terreno, as como su distribucin, resulta
de gran utilidad para estimar la coincidencia de stos con los periodos naturales
de los edificios. Cuando ambos periodos son iguales o muy parecidos, se presenta
un fenmeno de amplificacin del movimiento de la estructura que generalmente se
denomina como resonancia, el cual implica desplazamientos y deformaciones mucho ms grandes, en comparacin con otros edificios donde no hay coincidencia de
periodos y, en consecuencia, una mayor probabilidad de dao.
En la Figura 2.3, se observa un mapa de isoperiodos correspondiente la ciudad de
Valdivia. En el mapa se representan con distintos colores, una distribucin de 6
intervalos, de las mediciones realizadas [Val12].

22

2.2 Mtodo de Nakamura

Figura 2.3.: Mapa Isoperiodos Ciudad de Valdivia

2.2.

Mtodo de Nakamura

Nakamura propuso en el ao 1989 [Nak89] que el cociente entre la componente


horizontal y la componente vertical en el dominio frecuencial era un buen estimador
del efecto de sitio. Esta tcnica se aplic por primera vez a microtemblores de zonas
urbanas de Japn. Este mtodo ha sido ampliamente difundido en los ltimos aos
y utilizado por la fiabilidad de sus resultados, acompaada de rapidez y economa
en su ejecucin.
Nakamura supone que la razn (H/V ) o razn espectral corresponde a la funcin de
transferencia de los niveles superficiales sometidos a movimientos horizontales. Basa
su trabajo en la hiptesis de que un microtemblor est compuesto por un cuerpo
principal de ondas y por ondas superficiales RayLeigh, inducidas por fuentes artificiales (trenes, trfico, etc). Se considera que las componentes horizontal y vertical
del registro son similares en la base rocosa.
Dichos registros son amplificados por las capas de suelo blando que yacen sobre el
sustrato rocoso. Se considera adems que las componentes horizontales se amplifican
debido a multirreflexiones de la onda S, en tanto que los microterremotos verticales
lo hacen por multirreflexiones de la onda P. La velocidad de propagacin de la
onda P es generalmente mayor que 100 m/s, y los movimientos con frecuencias de
aproximadamente 10 Hz o menos pueden no amplificarse, debido a la multirreflexin
en las capas superficiales de varias decenas de metros de espesor.

23

2.2 Mtodo de Nakamura


Por otra parte, el efecto de la onda Rayleigh aparece resaltado en el movimiento
vertical, su efecto puede conocerse usando la componente vertical, haciendo la relacin entre la superficie y el subestrato. Normalmente el efecto de la onda Rayleigh
es cercano a uno, cuando la relacin anterior es aproximadamente uno.
En la medida en que la relacin aumenta, el efecto de las ondas Rayleigh puede
volverse crtico. La propuesta de Nakamura pretende eliminar este efecto.
Para lograrlo, Nakamura realiz los siguientes postulados:
Las ondas RayLeigh slo afectan a la parte superior del nivel y en cambio no
afectan a la base del nivel de suelo.
El efecto de las ondas RayLeigh es visible principalmente en el movimiento vertical,
y por tanto se puede definir la razn.
ES =

SV S
SV B

(2.1)

Donde SV S y SV B corresponden al espectro vertical de amplitudes de Fourier en la


superficie S y en la base B del nivel. Esta razn toma el valor de uno si no existe
presencia de ondas RayLeigh y menor a uno si se incrementa la presencia de estas
ondas superficiales.
El efecto de la onda RayLeigh es igual en superficie para el movimiento vertical
y horizontal.
Nakamura propuso una funcin de transferencia modificada corregida para ondas
RayLeigh. Una expresin de la funcin de transferencia clsica es:
ST =

SHS
SHB

(2.2)

Donde SHS y SHB corresponden al espectro horizontal de amplitud de Fourier en la


superficie S y la base B del nivel del suelo. La funcin de transferencia modificada
propuesta por Nakamura tiene la siguiente expresin:
SM =

ST
SHS SV B
=
ES
SHB SV S

(2.3)

Donde SM es la funcin de transferencia modificada, ST es la funcin de transferencia clsica segn ec. 2.2., ES es la razn definida en la ec. 2.3, y SV B y SV S
son los espectros de amplitudes de Fourier en la base y la superficie del estrato,
respectivamente.
Esta ecuacin se simplific a partir de los resultados en tres ciudades Japonesas que
para el rango frecuencial de 0.1 a 20 Hz establecan que en basamentos rocosos la

24

2.2 Mtodo de Nakamura


propagacin es la misma en todas las direcciones, como es posible observar en la
figura.

Figura 2.4.: Ejemplo de relacin espectral horizontal y vertical en un emplazamiento en roca


Luego se tiene:
SHB
=1
SV B

(2.4)

Por tanto, la funcin de transferencia modificada de Nakamura se expresa como:


SM =

SHS
SV S

(2.5)

Es decir, la funcin de transferencia se puede obtener con una nica medida en


la superficie. El fundamento terico de esta tcnica se ha cuestionado [CB94] pero
estudios experimentales muestran que, utilizando esta tcnica, para suelos blandos
se muestra un pico claro que est correlacionado con la frecuencia de resonancia
fundamental.

2.2.1.

Estimacin del espectro de Nakamura

Hasta hoy en da no se ha realizado ningn estudio acerca de la robustez de la


metodologa para calcular la razn espectral ni el efecto que la modificacin de
parmetros podra provocar en el clculo del espectro final y por ende en el clculo
de la frecuencia natural del suelo.
El mtodo de Nakamura requiere como entrada el clculo de espectros por componentes del registro de microvibracin. En general los algoritmos de clculo de la
razn espectral poseen los siguientes pasos [JEA12, eA08, SLM10]:

25

2.2 Mtodo de Nakamura


Filtrado pasa-banda de cada una de las componentes para atenuar ruidos de
baja y alta frecuencia.
Separacin en intervalos de tiempo (ventanas) de cada una de las seales componentes (NS, EO, y V). No hay uniformidad en el largo del intervalo de tiempo
ni tampoco en si existe solape entre las distintas ventanas. Generalmente se
encuentran ventanas de entre 20 segundos y 1 minuto con un solape que va
entre 0 % y 10 %.
Clculo de la FFT dentro de cada una de las ventanas.
Suavizado de la FFT dentro de cada una de las ventanas. Normalmente en esta
etapa se difiere en los filtros digitales que se utilizan. El filtro ms utilizado es
el llamado de Konno-Ohmachi [KO98]. Sin embargo tambin se han utilizado
filtros triangulares y a promedio mvil.
Clculo de la razn H/V dentro de cada ventana.
Promediacin de todas las ventanas de las razones (H/V ).
Estos pasos requieren la eleccin de bastantes parmetros tales como:
Longitud de la ventana.
Solape de las ventanas.
Filtro de procesamiento.
Filtro de suavizado.
Opcionalmente en algunos casos se extraen de la seal los transientes con el fin de
obtener una mejor estimacin del espectro por ventanas. Esto se realiza tanto de
forma manual como automtica. Para ello se comprueba que cada ventana no est
fuertemente contaminada por vibraciones aleatorias de gran energa comparadas con
el ruido de fondo, que se est midiendo en ninguna de sus tres componentes. Esto es
as, ya que la nica explicacin para obtener aumentos importantes de la vibracin
medida es que alguna fuente importante y transitoria relativamente cercana haya
generado principalmente una onda directa [JEA12].

2.2.1.1.

Desventajas en la estimacin del espectro de Nakamura

Es sabido en el mbito de procesamiento de seales que el clculo de espectros


por medio de FFT promediado no es robusto cuando se aplica a seales no
estacionarias [Kay99, SM05] como las de microvibraciones.

26

2.2 Mtodo de Nakamura


La mayor limitacin del mtodo (H/V ), es que presenta vacos tericos [KO98],
lo que explicara el hecho de que en ciertos casos el mtodo funcione [LCG94]
y en otros no.
Investigaciones tericas y experimentales [CB94] han demostrado que el mtodo nicamente permite la determinacin del perodo predominante.
An no se ha encontrado una relacin numrica entre el perodo obtenido con
microtemblores y el obtenido con movimientos ssmicos fuertes, como si ocurre
en el caso de movimientos ssmicos dbiles y fuertes.
No existe un estudio acabado sobre distintas tcnicas para estimar los espectros
de Nakamura.

2.2.2.

Resultados de la estimacin de espectros

Como se menciono ms arriba, no existen muchos estudios sobre las tcnicas para
estimar los espectros de Nakamura, desde el punto de vista del tratamiento frecuencial y el manejo de los parmetros del procesamiento de las seales. Precisamente
la parte fundamental de nuestro estudio es la estimacin del espectro de Nakamura y la interpretacin de los resultados obtenidos, tomamos como referencias, dos
estudios importantes que se han realizado, y que intentan estandarizar o ser una
gua del uso del mtodo. El primero es un estudio realizado por la unin europea
bajo el nombre de proyecto Sesame [Pro04], en el cual se realizan diversas pruebas
entregando los parmetros que presentan mejores resultados. El segundo estudio corresponde al trabajo realizado por el ingeniero Raul Ritta [JEA12] para su tesis de
maestra, y un posterior paper donde analiza la determinacin del periodo fundamental con vibracin ambiental usando Nakamura. Estos dos estudios nos brindan
algunas referencias enfocadas al anlisis frecuencial que podemos usar en nuestro estudio, como lo son algunos parmetros, filtros e interpretaciones de formas de peaks
que se pueden presentar.

2.2.2.1.

Estudio SESAME

El poyecto de la unin europea SESAME [Pro04], tiene por objetivo entregar algunas
de las principales directrices, para la adquisicin de datos y estimacin de espectros y
pretende controlar mejor los fundamentos tericos y metodolgicos de la utilizacin
del ruido ambiental, para estimar de forma vlida el efecto de sitio. Este estudio fija
algunos parmetros como longitud de ventana, solape, eliminacin de transitorios,
los cuales son sueguidos en los diferentes estudios realizados en ciudades de todo el
mundo [Alf].

27

2.2 Mtodo de Nakamura


En este estudio, se presenta una clasificacin de cinco tipos principales de peaks
[Pro04], con los cuales el realizador de un estudio de microzonificacin, se puede encontrar, segn los criterios mencionados anteriormente. A continuacin se muestran
las figuras de cada tipo de peak, con una pequea resea sobre su interpretacin.
Caso 1: Peak Claro

Figura 2.5.: Caso 1 peak claro

Aqu es evidente el predominio de un peak de mxima amplitud, SESAME da como interpretacin para este caso que todos los criterios utilizados en su estudio se
cumplen y la frecuencia fundamental de este sitio puede ser estimada de manera
fiable.
Caso 2: Peak poco claro de baja frecuencia

Figura 2.6.: Caso 2: Peak poco claro de baja frecuencia

28

2.2 Mtodo de Nakamura


Para este caso, se presenta un peak de baja frecuencia y baja amplitud, las directrices
entregadas por SESAME, recomiendan que cuando se den estos casos se deben
realizar ms pruebas en el sitio, o tomar en cuenta una serie de pasos que se detallan
en uno de los captulos de la publicacin del estudio.
Caso 3: Dos peaks (f1>f0)

Figura 2.7.: Caso 3: Dos peaks (f1>f0)


Aqu se presentan dos peaks de mxima amplitud, uno de ellos con menor frecuencia
que el otro. La interpretacin que le da SESAME a este caso es que no se puede
afirmar con exactitud que la frecuencia natural es el peak de baja frecuencia, realizar
un suavizado con una banda ms estrecha podra satisfacer los criterios y dar un
mejor resultado.
Caso 4: Peak ancho o mltiples peaks

Figura 2.8.: Caso 4: peak ancho o mltiples peaks

29

2.2 Mtodo de Nakamura


Para este caso ms de dos peak cumplen con el criterio de amplitud de SESAME,
la disponibilidad de informacin geolgica, y geofsica del terreno, puede permitir
discriminar cual de los peaks es el de la frecuencia fundamental, la explicacin de
la presencia de los dems peaks puede deberse a las caractersticas del sitio de
medicin, como por ejemplo el borde de un valle. Comparar este registro con uno
cercano puede servir para identificar la frecuencia natural.
Caso 5: Peaks agudos y de origen industrial

Figura 2.9.: Caso5: peaks agudos y de origen industrial

Para este caso se observa la presencia de un peak estrecho, el cual es caracterstico


por presentarse en sitios donde existe la influencia de maquinas industriales, para
corroborar la presencia de este peak, se pueden variar los parmetros de suavizado,
la respuesta a esto, observada por SESAME, es que el peak se estrecha mas y mas,
y a la vez su amplitud crece. Otra caracterstica para confirmar la presencia de este
peak, es que se presenta en los espectros de Fourier de las tres componentes de la
seal.

2.2.2.2.

Estudio Ral Ritta en Ciudad Mayagez

En el trabajo realizado por el Sr Ritta [JEA12], se estudio la aplicabilidad de la


tcnica de los cocientes espectrales (H/V ) a registros de microvibracin ambiental
para la determinacin de la frecuencia fundamental del suelo. Para realizar dicho
estudio se tomaron 134 registros en diferentes sitios de la ciudad. Para el procesado
de las seales, se desarrollo un algoritmo en Matlab, los pasos para el procesado son:
1. remocin de la media del registro.

30

2.3 PDS en microvibraciones y comparacin con otras reas de geologa.


2. seleccin automtica de ventanas basado en parmetros definidos por el usuario.
3. multiplicacin de cada ventana por una ventana de Hanning.
4. clculo de la transformada de Fourier y suavizado de la misma para cada
ventana.
5. clculo de la media de los espectros de Fourier de todas las ventanas y finalmente.
6. clculo de los cocientes espectrales NS/V, EO/V y Hc/V (donde Hc son las
componentes horizontales combinadas).
El algoritmo desarrollado, incluy un modulo para conservar las partes ms estacionarias de ruido ambiental y evitar los transitorios asociados a fuentes puntuales
(vehculos, maquinarias, etc.). Este algoritmo es conocido como anti-triggering. Para
demostrar la estabilidad del mtodo, se tomaron registros de 1 da de duracin, para
luego procesar distintas horas. Como resultado a este procedimiento se tiene que la
frecuencia fundamental fue la misma, independientemente de la hora del da en que
se toma la medicin. Lo cual demuestra la estabilidad del mtodo y que la frecuencia
fundamental es una propiedad del terreno y no depende del tipo de excitacin al
cual se someta.
Como conclusiones del estudio del Sr Ritta, este menciona que el mtodo, entrega
muy buenos resultados en comparacin a otras tcnicas, adems es una herramienta
simple y econmica y se puede aplicar en cualquier momento, sin ser necesaria la
espera de la ocurrencia de un terremoto.

2.3.

PDS en microvibraciones y comparacin con


otras reas de geologa.

Desde hace aos la geologa ha necesitado del procesamiento digital de seales (PDS)
para analizar y extraer informaciones adicionales de los distintos datos y registros
con los que se trabaja [MGC99, Tra10]. Las reas pioneras dentro de geologa que
hicieron uso del procesamiento de seales fueron la sismologa y la vulcanologa
[TK03, HO10, Sch94, Les09]. En ambas reas la motivacin principal para hacer
uso del PDS fue la estimacin de espectros de frecuencias y el uso de filtros digitales
para acondicionar las seales [eA07a, ZTC03, JWH02, BS99, JNS03, LGM02].
El rea del estudio de riesgo ssmico por medio de la caracterizacin dinmica de
suelos no escapa al uso de PDS. El estudio de microvibraciones para el clculo

31

2.3 PDS en microvibraciones y comparacin con otras reas de geologa.


de la razn espectral (H/V ) de Nakamura hace uso de filtros digitales y tcnicas
de estimacin de espectro desde los tiempos en que esta tcnica fue propuesta a
principios de la dcada de los 90 [Her08, Nak08]. Sin embargo las tcnicas de PDS
utilizadas en reas como sismologa y vulcanologa estn bastante ms avanzadas
que en el rea de microvibraciones, lo que nos lleva a preguntarnos el porqu de
este hecho y en que afecta los estudios de riesgo ssmico que hoy en da se realizan
utilizando la tcnica de Nakamura.
Tal como en las reas de investigacin antes mencionadas el rea de caracterizacin
dinmica de suelos necesita del uso de filtros digitales para acondicionar y preprocesar la seal antes de calcular la razn espectral (H/V ). Los filtros utilizados en
el rea de microvibraciones estn acordes con los avances del PDS en dicha rea y no
se diferencian mayormente de los filtros utilizados en sismologa y vulcanologa. Sin
embargo el anlisis frecuencial que se les realiza a los registros de microvibraciones
requieren de una discusin ms detallada y de un anlisis crtico al mtodo de
estimacin de espectros que hoy en da se utiliza en dicha rea de investigacin.
Es sabido que para obtener la razn espectral (H/V ), en primer lugar es necesario
calcular los espectros de cada una de las componentes norte-sur (NS), este-oeste
(EO) y vertical (V) [JEA12, eA08]. Este es el rol del anlisis frecuencial, el cual corresponde al proceso de determinar contenidos frecuenciales en una seal. El anlisis
espectral pude pensarse como un prisma matemtico que descompone una seal
en sus frecuencias constituyentes tal como un prisma descompone luz en sus colores
(longitudes de onda) constituyentes [Sem04]. Un proceso fsico puede durar toda la
eternidad. Sin embargo solo una parte del proceso puede ser captado. Esto quiere
decir que cualquier intento de calcular las componentes frecuenciales de una seal
solo ser una estimacin del espectro real. El problema de estimar el espectro de
una seal se resume a cmo encontrar el verdadero espectro de esta seal sabiendo
que se dispone de un nmero de datos finito y en tiempo discreto, ms an sabiendo
que los datos estn perturbados por ruido [Kay99, SLM10].
En el caso de la tcnica de la razn espectral (H/V ) de Nakamura se ha optado a lo
largo de los aos de investigacin por una tcnica de estimacin del espectro basada
en la Transformada de Fourier calculada con la clsica Transformada Rpida de
Fourier ms conocida como FFT. En general no est del todo claro el procedimiento
detallado por el cual se calcula el espectro (H/V ) utilizando FFT. Normalmente los
ingenieros e investigadores utilizan softwares, muchas veces propietarios, que son
unas verdaderas cajas negras para los usuarios.
Hoy en da el anlisis espectral est bastante avanzado con respecto a lo que se realiza en el rea de microvibraciones. Las ltimas tres dcadas han estado marcadas por
la gran cantidad de tcnicas y algoritmos nuevos que se han propuesto en anlisis
espectral. Estas nuevas tcnicas apuntan al anlisis frecuencial de seales no estacionarias o que presentan dinmicas transientes. Entre estas nuevas tcnicas modernas
de anlisis espectral se destacan las basadas en Transformada Wavelet [Mal08] y

32

2.3 PDS en microvibraciones y comparacin con otras reas de geologa.


en modelado paramtrico [Kit10, MKBB05]. La primera tcnica es muy til para el
anlisis de transientes, mientras que la segunda para la estimacin espectral variante
en el tiempo (para seales no estacionarias) con una baja cantidad de datos. Estas
tcnicas modernas de anlisis espectral permiten una mayor robustez en la estimacin del espectro utilizando una baja cantidad de observaciones del proceso. Ms
an las tcnicas propuestas son generalmente inmunes a ruidos transientes, lo que
genera espectros no dependientes de las no estacionaridades presentes en los datos.
Lo que se propone para contribuir al avance de esta rea y de lo cual esta tesis es
una parte del proyecto, es desarrollar una nueva tcnica de estimacin del espectro
(H/V ) que tenga las tres caractersticas principales siguientes:
Garantice una mayor confiabilidad en la estimacin de la frecuencia natural
utilizando tcnicas de anlisis espectral que sean robustas frente a dinmicas
transientes presentes en los registros. Esto se realizar sin necesidad de suprimir los transientes de la seal, sino que utilizando tcnicas robustas variantes
en el tiempo inmunes a dinmicas transientes.
Permita la estimacin de la frecuencia natural en un registro de menor duracin
a los que deben realizarse hoy en da. Esto mejorara significativamente la
cantidad de tiempo utilizado en realizar los registros de microvibraciones y
analizarlos.
Entregue un espectro variante en el tiempo (espectro (H/V ) en tiempo y en
frecuencia) de fcil anlisis visual que permita al usuario tomar la decisin
acerca de la frecuencia natural observando el registro completo tanto en tiempo
como en frecuencia. Esto ser bastante til sobre todo en registros de suelos
de naturaleza compleja, donde ms de una frecuencia natural puede estar
presente.
Para realizar esto se propone primeramente una estimacin espectral mediante modelado paramtrico autoregresivo (AR) [Kit10, MKBB05] utilizando parmetros variantes en el tiempo con el fin de que el modelo se adapte a las no estacionalidades
de la seal. Estos parmetros sern luego estimados incorporando estadsticas robustas de Huber [Ill03, HA95] con el fin de que el modelo sea inmune a las dinmicas
transientes, especficamente a los ruidos impulsionales. Finalmente una vez el modelo caracterizado y sus parmetros estimados se proceder a estimar la funcin
de transferencia del sistema cuyos polos entregarn las dinmicas ms significativas
presentes en las seales. Con el fin de encontrar el polo directamente relacionado
con la frecuencia natural del suelo se utilizarn mtodos de seguimiento de polos
[Kit10, eA07b].

33

3. Captulo III CONCEPTOS


TERICOS DEL PROYECTO
En este capitulo se vern las tcnicas tradicionales para la estimacin de espectros,
utilizadas en el mtodo de Nakamura, adems de algunas de las tcnicas modernas
utilizadas en el estudio realizado.

3.1.

Tcnicas tradicionales de estimacin del


espectro de Nakamura.

Desde hace algn tiempo el anlisis espectral espacial y temporal se ha introducido


activamente en la investigacin geofsica. El mtodo tradicional, es el anlisis espectral de la serie de tiempo por la transformada de Fourier, qu proporciona las
caractersticas de una seal variable temporal en el espacio de frecuencia. Mientras
que los mtodos adaptativos (no lineales) son una herramienta poderosa del anlisis
espectral para la deteccin de periodicidades ocultas en la serie de tiempo y ms
cuando las series de tiempo registradas son cortas. Tanto los modelos lineales como
los no lineales utilizan la transformada de Fourier para pasar del espacio de tiempo
al espacio de frecuencia y viceversa [Eli08].
El anlisis por medio de la transformada de Fourier provee informacin de las caractersticas de una seal variable temporal. Este proporciona solo su frecuencia
contenida sobre el registro de longitud entera y toda la informacin relacionada a la
fase para eventos individuales, comprendiendo la serie de tiempo perdida en la transformacin del espacio temporal al espacio de frecuencias. As es posible identificar
la magnitud y frecuencia de varios eventos incluidos en la serie de tiempo.
En el mtodo de Nakamura, se utiliza esta tcnica tradicional basada en la transformada de Fourier, para estimar los espectros. Como ya se menciono en Subseccin 2.2.1
para estimar los espectros utilizados en el mtodo de Nakamura, es necesario emplear
algunas tcnicas tradicionales.
A modo introductorio, se recordaran de forma breve los pasos para la estimacin del
espectro de Nakamura, los cuales son:

34

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


Filtrado pasa-banda de cada una de las componentes para atenuar ruidos de
baja y alta frecuencia.
Separacin en intervalos de tiempo (ventanas) de cada una de las seales componentes (NS, EO, y V).
Clculo de la FFT dentro de cada una de las ventanas.
Suavizado de la FFT dentro de cada una de las ventanas.
Clculo de la razn H/V dentro de cada ventana.
Promediacin de todas las ventanas de las razones (H/V ).
A continuacin se muestran los pasos anteriores en un diagrama de bloques, ademas
de algunas imagenes de los procesos ms destacables.

Figura 3.1.: Diagrama de bloque estimacin de espectro de Nakamura.

Figura 3.2.: Seal de entrada de tres componentes EO, NS, y V.

La Figura 3.2, muestra el tipo de seal de microvibracin con el cual se trabaj,


la cual posee tres componentes (E-O, N-S, y Vertical). Esta seal primeramente es
filtrada con un filtro pasa-banda, cuyo objetivo es eliminar los ruidos de baja y alta
frecuencia.

35

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.

Figura 3.3.: Componentes de una seal de microvibracin con solape de ventanas


de 10 % y 30 %.
Luego del filtrado, la seal es separada en intervalos de tiempo o ventanas, como se
observa en la Figura 3.3.

Figura 3.4.: Ventana deslizante, proceso de calculo del espectro FFT.


Los registros con los que se cuentan tienen una duracin de 15 minutos. En la
Figura 3.4, se muestra el proceso de calculo de la FFT a traves de la seal, se
realizan 22 FFT sucesivas, con un largo de 40 seg. de seal cada una. Para este caso
no hay solapamiento de ventanas.

Figura 3.5.: Espectro FFT promediado de cada componente de una seal.

36

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.

Figura 3.3.: Componentes de una seal de microvibracin con solape de ventanas


de 10 % y 30 %.
Luego del filtrado, la seal es separada en intervalos de tiempo o ventanas, como se
observa en la Figura 3.3.

Figura 3.4.: Ventana deslizante, proceso de calculo del espectro FFT.


Los registros con los que se cuentan tienen una duracin de 15 minutos. En la
Figura 3.4, se muestra el proceso de calculo de la FFT a traves de la seal, se
realizan 22 FFT sucesivas, con un largo de 40 seg. de seal cada una. Para este caso
no hay solapamiento de ventanas.

Figura 3.5.: Espectro FFT promediado de cada componente de una seal.

36

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


En la Figura 3.5, se aprecian los espectros promediados de la FFT de cada una de
las componentes de la seal.

Figura 3.6.: Razn (H/V ) promediada.


La Figura 3.6, muestra el final del proceso, en donde se obtiene el espectro de la razn
(H/V ) promedio de cada una de las ventanas. En la figura se aprecian en distintos
colores, los espectro (H/V) de cada ventana, la lnea negra continua corresponde
al espectro promedio, mientras que las lneas negras punteadas, son la desviacin
estndar.
A continuacin se explicara la teora de los procesos necesarios para la estimacin
de los espectros de Nakamura.

3.1.1.

Filtros

Un filtro es un elemento que transmite una seal de tal manera que su salida es el
resultado de convolver la seal de entrada con el impulso de respuesta de la funcin
h(t) del filtro. En el dominio de la frecuencia este corresponde a una multiplicacin
del espectro de la frecuencia de la seal por la respuesta de la funcin del filtro
H(f ). El filtro se caracteriza por su pulso de respuesta en el dominio del tiempo
y por su respuesta en frecuencia en el dominio de la frecuencia [Mig10]. Ambas
caracterizaciones contienen la misma informacin sobre el filtro y son relacionadas
va la transformada de Fourier:

H(f ) = F [h(t)]

(3.1)

37

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


En palabras mas simples, un filtro es un elemento que discrimina una determinada
frecuencia o gama de frecuencias de una seal que pasa a travs de l, pudiendo
modificar tanto su amplitud como su fase.
Se puede llegar a expresar matemticamente la funcin de transferencia en forma de
fraccin mediante las transformaciones en frecuencia adecuadas. Se dice que los valores que hacen nulo el numerador son los ceros y los que hacen nulo el denominador
son polos [Mor09].

H(f ) =

numerador(f )
denominador(f )

(3.2)

El nmero de polos y ceros indica el orden del filtro y su valor determina las caractersticas del filtro, como su respuesta en frecuencia y su estabilidad.

3.1.1.1.

Filtro lineal

Un filtro lineal es aquel filtro que aplica un operador lineal a una seal variable
en el tiempo. Son usados ampliamente en procesamiento de seales. Una de sus
aplicaciones ms frecuentes es la eliminacin de frecuencias no deseadas de una
determinada seal de entrada o, al contrario, discriminar una determinada frecuencia
de las dems [Mor09].
Los filtros lineales se pueden clasificar segn su funcin de transferencia en dos
grupos que se describen a continuacin:
Respuesta en amplitud
Los fitros lineales pueden dividirse en dos clases: filtros de respuesta infinita (IIR,
del ingls infinite impulse response) y filtros de respuesta finita (FIR, del ingls finite
impulse response) [Smi97]:
Los filtros FIR (que slo puede ser implementados en tiempo discreto) pueden
ser descritos como una suma ponderada de entradas con un determinado retardo. Para estos filtros, si la entrada en un determinado instante es cero, la
salida ser cero a partir de un instante posterior a los retardos inducidos por
el filtro. De este modo, solo existir respuesta por un tiempo finito.
Los filtros IIR, por el contrario, pueden presentar salida an cuando la entrada
sea cero, si las condiciones iniciales son distintas de cero. La energa del filtro
decaer con el tiempo, pero no llegar a ser nula. Por tanto, la respuesta al
impulso se extiende infinitamente.

38

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


Respuesta en frecuencia
Hay varios tipos de filtros lineales en lo que respecta a su respuesta en frecuencia:
Filtro paso bajo: permite el paso de frecuencias bajas.
Filtro paso alto: permite el paso de frecuencias alto.
Filtro paso banda: permite el paso de un rango intermedio de frecuencias.
Filtro banda eliminada: bloquea el paso de un rango intermedio de frecuencias.
Filtro paso todo: permite el paso de todas las frecuencias, pudiendo modificar
su fase.

3.1.1.2.

Filtro de Butterworth

El filtro de Butterworth es uno de los filtros electrnicos ms bsicos, diseado para


producir la respuesta ms plana que sea posible hasta la frecuencia de corte. En
otras palabras, la salida se mantiene constante casi hasta la frecuencia de corte,
luego disminuye a razn de 20n dB por dcada, donde n es el nmero de polos
del filtro [MGSJV10]. El espectro de potencia (o magnitud cuadrtica) de un filtro
(pasa-bajos) de Butterworth se puede expresar como:

| H() |2 =

1
2N
1+
c


(3.3)

donde c es la frecuencia de corte y N es el orden del filtro. Esta funcin tiene 2N


polos, en tanto que | H() | tiene N polos. Este filtro se usa como anti-alias y se
caracteriza por tener una banda de paso muy plana. Esto es as porque las primeras
2N 1 derivadas de espectro de potencia en = 0 son nulas.
El filtro Butterworth ms bsico es el tpico filtro pasa bajo de primer orden, el cual
puede ser modificado a un filtro pasa alto o aadir en serie otros formando un filtro
pasa banda o elimina banda y filtros de mayores rdenes.
Segn lo mencionado antes, la respuesta en frecuencia del filtro es extremadamente
plana (con mnimas ondulaciones) en la banda pasante. Visto en un diagrama de
Bode con escala logartmica, la respuesta decae linealmente desde la frecuencia de
corte hacia menos infinito. Para un filtro de primer orden son -20 dB por dcada.

39

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


El filtro de Butterworth es el nico filtro que mantiene su forma para rdenes mayores (slo con una cada de ms pendiente a partir de la frecuencia de corte). Este tipo
de filtros necesita un mayor orden para los mismos requerimientos en comparacin
con otros, como los de Chebyshev o el elptico.

Figura 3.7.: Filtro Butterworth de distintos ordenes

Este filtro se utiliz como filtro de procesamiento, para la seal filtrada utilizada en
el software Matlab, el filtro fue un pasa-banda de orden 3 con frecuencias de corte de
0,05 a 10 Hz. En tanto, en el estudio de sensibilidad de parmetros con el software
Geopsy, se utilizaron rdenes 2 y 3, y frecuencias de corte de 0,1 a 10 Hz 0,1 a 20
Hz, y 0,1 a 15 Hz.

3.1.2.

Ventaneo

Si el espectro de una seal es dependiente del tiempo, es necesario usar elementos


suficientemente cortos de ella (con la certeza de que el espectro es constante sobre
cada segmento), para analizar el espectro. Tomar un segmento funcin del tiempo se
conoce como trabajar con ventanas. Este es uno de los procesos necesarios al aplicar
el mtodo de Nakamura, en donde la seal es seccionada en ventanas, las cuales se
pueden solapar o no, para luego calcular la FFT de cada una de estas ventanas. a
continuacin se describe el marco terico del proceso de ventaneo.
Las ventanas son funciones matemticas usadas con frecuencia en el anlisis y el
procesamiento de seales para evitar las discontinuidades al principio y al final de
los bloques analizados. En procesamiento de seales, una ventana se utiliza cuando
nos interesa una seal de longitud voluntariamente limitada. En efecto, una seal

40

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


real tiene que ser de tiempo finito; adems, un clculo slo es posible a partir de un
nmero finito de puntos. Para observar una seal en un tiempo finito, se multiplica
por una funcin ventana [Mig10]. As, cuando multiplicamos una seal S(t) por una
ventana, obtendremos nicamente los T primeros segundos de la seal. En vez de
estudiar la seal S(t), se estudia la seal truncada. La utilizacin de una ventana
cambia el espectro en frecuencia de la seal.
En microvibraciones, es conocido que para obtener un buen espectro, es necesario
contar con un registro de larga duracin, as de esta forma se tendr la suficiente
cantidad de datos para separar las seales en una cantidad suficiente de ventanas,
para as obtener un buen suavizado del espectro.
Se divide la seal en fragmentos cortos (short time) delimitados en el tiempo por una
funcin ventana. Una ventana no es nada ms que un tipo especfico de envolvente
que se aplica para un anlisis espectral. Para el caso de Nakamura la duracin de
la ventana est normalmente comprendida entre 20 s y 1 minuto, y los segmentos a
veces se solapan, variando entre 0 % y 10 %. A travs del anlisis espectral individual
de cada segmento enventanado, se obtiene una secuencia de medidas (de espectros)
que constituyen el espectro variable a lo largo del tiempo o sonograma [Ign09].

3.1.2.1.

Operacin del Ventaneo

Un segmento de longitud T de la seal es:

= S(t)(t)
S(t)

(3.4)

es, por el teorema de convolucin:


Y la transformada de Fourier de S(t)

S()
= S() ()

(3.5)

Donde S() y () son las transformadas de Fourier de S(t) y (t).


En la Figura 3.8 se muestra grficamente la operacin ventaneo.

41

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.

Figura 3.8.: Operacin ventaneo

3.1.2.2.

Tipos de Ventanas utilizadas

Existen distintos tipos de ventana que permiten obtener distintos resultados en el


dominio de las frecuencias. a continuacin se mostrarn los dos tipos de ventanas
que fueron utilizados en el anlisis realizado en el presente trabajo [Mig10].
En general, las ventanas poseen tres propiedades que son las que interesan a la hora
de realizar la eleccin de cual utilizar [JP07]:
Ancho del lbulo principal
Nivel de lbulos laterales
Pendiente de cada
Por lo general las ventanas con ancho del lbulo principal mayor tienen menores
niveles de lbulos laterales.
Los cuatro tipos de ventanas ms importantes son la ventana Rectangular, la ventana
Hanning, la ventana Hamming y la ventana Blackman. a continuacin se mostraran
los dos tipos de ventanas que fueron utilizados en el anlisis realizado en el presente
trabajo, las cuales son la Rectangular y la de Hanning.

42

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


3.1.2.3.

Ventana Rectangular

Una ventana Rectangular es aquella que posee un valor de 1 en todo el intervalo de


la ventana, y de 0 para cualquier otro valor [JP07, Mig10].

V (n) = 1

(3.6)

Figura 3.9.: Representacin ventana Rectangular


En la Figura 3.9 (izquierda) podemos observar la ventana Rectangular continua, la
misma que posee una amplitud constante, por lo que al utilizar una ventana de
este tipo la seal no se ver afectada. En la Figura 3.9 (derecha) encontramos el
equivalente de esta seal para seales muestreadas, es decir la ventana Rectangular
discreta.
El espectro de la ventana Rectangular es:

Figura 3.10.: Espectro de la Ventana Rectangular

43

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


El filtro caracterstico dado por la integral de la transformada de Fourier de la
ventana Rectangular se muestra en la Figura 3.11. El filtro tiene un lbulo que es
dos veces la anchura del espaciado4f .

Figura 3.11.: Filtro caracterstico de la ventana Rectangular

3.1.2.4.

Ventana Hanning

La ventana de Hanning se define a travs de la funcin [JP07, Mig10]:

1 1
2n
V (n) = Cos
2 2
N

N = 1, ...n

(3.7)

Para valores fuera del rango N = 1, ...n, se tiene una amplitud igual a cero.

Figura 3.12.: Representacin Ventana de Hanning o Hann

En la Figura 3.12 (izquierda) podemos observar la ventana de Hanning continua, la


misma que posee caractersticas especiales con respecto a la ventana Rectangular,
ya que atena la seal en los bordes de la misma. En la Figura 3.12 (derecha)

44

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


encontramos la ventana de Hanning discreta, la que ha sido muestreada a con una
factor de 32 muestras.
Sus filtros caractersticos se muestran en la Figura 3.13. El lbulo es 44f , doble
del ancho de la ventana Rectangular. El primer lbulo es mucho ms atenuado, y la
cada de la tasa es mucho ms rpida que la ponderacin Rectangular.

Figura 3.13.: Filtro caracterstico de la ventana de Hanning

Este ltimo tipo de ventana fue el ms utilizado en este trabajo, variando sus duraciones entre 20, 40, y 60 segundos, con solapes entre 0 y 50 %.

3.1.3.

Suavizado

Durante el calculo del espectro (H/V ), los espectros de Fourier pueden ser suavizados,
con el objetivo de aclarar el aspecto global de las curvas. Para el clculo del espectro
total, el espectro de Fourier de cada ventana de tiempo seleccionada se alisa y
entonces la curva promediada se calcula. Para el clculo de (H/V ), en primer lugar
los espectros horizontales (NS y EW) de Fourier, se combinan y luego se aplica el
suavizado en el espectro de Fourier horizontal resultante y en el espectro de Fourier
vertical [Geoa]. Por un lado el suavizado del espectro permite que no hayan valores
iguales a cero en la FFT calculada, ya que valores iguales a cero en los espectros
de cada una de las componentes significaran valores que tienden a infinito en el
espectro H/V, debido a que para calcular el espectro (H/V ) se necesita realizar una
operacin matemtica de divisin de espectros. Por otro lado el promediado de los
espectros (H/V ) sirve para dar una menor importancia a los espectros de segmentos
de seal que estn muy afectados por ruidos transientes (como pasos de personas,
movimientos de autos, etc.).
Los mtodos de suavizado generalmente utilizados son:

45

3.1 Tcnicas tradicionales de estimacin del espectro de Nakamura.


Suavizado de Konno y Ohmachi [Geob, KO98]: la funcin de Konno y Ohmachi
es de tipo trigonomtrico-logartmico y toma en cuenta la frecuencia central
de mxima amplitud aplicando un coeficiente exponencial. Esta suavizacin
utiliza un ancho de banda constante en escala logartmica y es fuertemente
recomendado ya que esta funcin de suavizado preserva el nmero diferente
de puntos de baja y alta frecuencia. Este suavizado es controlado por una
constante de suavizado que varia entre 0 a 100. una constante de 0 da un
suavizado muy fuerte, mientras que una constante de 100, da un suavizado muy
suave. El suavizado de Konno y Ohmachi, se define por la siguiente funcin:


sin b log

f 4
fc


f
b log
fc

(3.8)

Donde:
f = f recuencia
fc = f recuencia central, la que realiza el suavizado
b = coef iciente de ancho de banda
La funcin de suavizado de Konno y Ohmachi es muy recomendable para el anlisis de frecuencia, ya que garantiza un nmero constante de puntos de baja y alta
frecuencia [Geob].

Figura 3.14.: Comparacin entre tcnica H/V estndar (a) y H/V incorporando la
funcin de suavizado de Konno y Ohmachi (b)
Constante de Suavizado: esta funcin de suavizado tiene una forma triangular
centrada en la frecuencia actual y su anchura es igual al "ancho de banda"
dado por el usuario. Este ancho de banda vara de 0 (suavizado suave) a 100
(suavizado fuerte).

46

3.2 Tcnicas Modernas de Estimacin de Espectros


Suavizado Proporcional: esta funcin de suavizado tiene una forma triangular
y su anchura depende de la frecuencia actual. La anchura media es definido
por porcentaje*frecuencia. El valor de "porcentaje" no puede ser mayor o igual
a 100 % (suavizado fuerte).
Estos mtodos de suavizado descritos, son los que dispone el toolbox del Software Geopsy [Geoa], siendo el ms conocido y mayormente utilizado el de Konno y
Ohmachi.
4

x 10

1.5

0.5

10

15

Figura 3.15.: Espectro FFT y Suavizado del mismo.

La Figura 3.15 muestra el espectro de Fourier de una de las componentes de un


registro de microvibracin. La lnea negra representa el suavizado realizado con el
software Geopsy, donde se utilizo un suavizado de Konno-Ohmachi con constante
de suavizado de 40.

3.2.

Tcnicas Modernas de Estimacin de Espectros

La mayora de las iniciativas que se han desarrollado a lo largo de los aos, para la
estimacin espectral en estudios de suelos, se basan en la transformada de Fourier
calculada con la clsica Transformada Rpida de Fourier ms conocida como FFT.
Hoy en da el anlisis espectral est bastante avanzado en comparacin con lo que
se realiza en microvibraciones, ya que en los ltimas tres dcadas se han desarrollado nuevos algoritmos en otras reas como sonares, radares, bioseales [Sem04],
telecomunicaciones [Pro01], ingeniera biomedica [CJ07, SL05, NA09], adems de
reas similares a las de microvibraciones, como lo son sismologa y vulcanologa

47

3.2 Tcnicas Modernas de Estimacin de Espectros


[Seb10, JE00]. En estas ltimas destaca el uso de la wavelet, con la cual se han obtenido grandes avances y aplicaciones en los anlisis espectrales de estas dos grandes
reas. Estas nuevas tcnicas apuntan al anlisis frecuencial de seales no estacionarias o que presentan dinmicas transientes [JCD]. A continuacin se muestran las
principales tcnicas modernas de procesado de seales utilizadas para estimar espectros de procesos de caractersticas no estacionarias y que sern utilizadas durante
este trabajo.

3.2.1.

Anlisis en Tiempo-Frecuencia

La mayora de las seales no son estacionarias. La amplitud y el espectro de una


seal pueden variar en el espacio temporal. Un ejemplo es la emisin acstica de un
barco, con variaciones debidas a los cambios de velocidad y el encendido y apagado
de la maquinaria de a bordo. Las seales de radar con modulacin frecuencial tienen
un tiempo dependiente del espectro. Al hablar, las frecuencias tambin son distintas
en los diferentes segmentos de la palabra hablada. Una completa caracterizacin de
las seales no estacionarias en el dominio de la frecuencia debe adems incluir el
aspecto temporal, con el que se obtiene el anlisis tiempo-frecuencia de una seal
[Mig10].

3.2.1.1.

Transformada corta de Fourier (STFT)

Como ya es sabido, la transformada de Fourier constituye una herramienta mediante la cual se puede obtener informacin sobre como est distribuida la energa de
una seal a travs de sus distintas componentes de frecuencia, es decir, se pueden
conocer todas las componentes de frecuencia existentes en la seal y sus respectivos
aportes energticos. Todo lo anterior se puede resumir diciendo que la transformada
de Fourier tiene una perfecta resolucin en frecuencia lo que la hace una herramienta muy til para el anlisis de seales estacionarias. Sin embargo, ella no puede ser
aplicada con el objeto de obtener informacin precisa de cundo o dnde las diferentes componentes de frecuencia se encuentran en la seal como es el caso de seales
quasi - estacionarias o no estacionarias cuyo contenido espectral vara con el tiempo.
En otras palabras, la transformada de Fourier posee una resolucin en tiempo nula
[Mig10].
Con la intensin de resolver el problema de resolucin en tiempo de la transformada
de Fourier, Denis Gabor [Gab46] adapt la transformada utilizando un procedimiento llamado ventanamiento. Este procedimiento consiste en dividir una seal x(t) en
pequeos segmentos a travs del tiempo de tal manera que podamos asumir que para cada segmento la seal es estacionaria y as calcular la Transformada de Fourier
clsica para cada porcin de la seal.

48

3.2 Tcnicas Modernas de Estimacin de Espectros


La transformada corta de Fourier (Short-Time Fourier Transform: STFT) es una
adaptacin de la TF introducida por Gabor para incluir la dependencia temporal
mediante la utilizacin de ventanas centradas en distintos instantes temporales de
la seal sobre las que se calcula la TF [F.10]. Las ventanas, w(t), son de duracin
constante para todas las frecuencias y se supone que la seal es estacionaria en las
mismas. La STFT se define como:

ST F Tx (t, f ) =

[x( )w(t )] ejf d

(3.9)

donde w(t ) es la funcion de ventaneo y es la variable que corre la ventana a


traves de la seal x(t).
La eleccin de la duracin y el ancho de la ventana constituyen parmetros de gran
importancia ya que a travs de estos podemos establecer el grado de resolucin
tanto de tiempo como de frecuencia que deseemos. Si la ventana es muy angosta
se analizar una porcin muy pequea de la seal lo que permite tener una buena
resolucin en tiempo pero una mala resolucin en frecuencia ya que se conocer slo
una mnima fraccin del espectro total existente en la seal. Por otro lado, si la
ventana en muy ancha se tendr una buena resolucin en frecuencia pero una mala
resolucin en tiempo. Por lo tanto un defecto de la STFT es que no puede entregar
una buena resolucin tanto en tiempo como en frecuencia y que estos parmetros
de la ventana son fijos.
En consecuencia, existe un compromiso entre buena resolucin en el tiempo o buena
resolucin en frecuencia. Para obtener una seal estacionaria, se elige una ventana
lo suficientemente estrecha que permita ser considerada as. Cuanto ms estrecha
sea la ventana mejor resolucin en el tiempo se obtendr y por lo tanto una mejor
representacin de estacionalidad y peor resolucin en frecuencia. Por lo tanto el
problema consiste en la eleccin de una ventana para el anlisis, dependiendo de la
aplicacin [Mig10].

Figura 3.16.: Explicacin grfica de la STFT

49

3.2 Tcnicas Modernas de Estimacin de Espectros


En la Figura 3.16 [LAS04] se muestra un ejemplo de una funcin ventana de tipo
Gaussiana. La primera funcin (roja) tiene la ventana localizada en t = t10 , la
segunda (azul) en t = t20 y la tercera (verde) en t = t30 . Estas ventanas corresponden
a tres transformadas de Fourier en tiempos distintos. Por lo tanto, se obtendr una
buena representacin tiempo-frecuencia de la seal de entrada.
Adems hay que hacer constar que la funcin ventana est caracterizada por ser
real, simtrica (la transformada de Fourier clsica es simtrica y, la STFT, al ser una
modificacin de esta, tambin lo es) y de duracin finita (esto es un inconveniente ya
que, al aplicarse la ventana tan slo a una parte de la seal, se causa una disminucin
de la resolucin en frecuencia).
Con esto se puede realizar un espectrograma, el cual es una representacin en tres
dimensiones (tiempo, frecuencia y ampitud) de la STFT de la seal. Con el espectrograma adems de conocer las componentes de frecuencia de la seal, tambin se
conoce su localizacin en el tiempo. Para el anlisis realizado, en la Figura 3.17 se
observa lo obtenido al aplicar STFT a los registros de microvibraciones. Se muestra
el espectrograma que agrupa los espectros en cada ventana, mostrando las principales frecuencias presentes en una de las componentes de un registro de microvibracin.
Destacadas en recuadros se aprecian las componentes frecuenciales que permanecen
constantes en el tiempo, siendo posibles "candidatas" a frecuencia natural luego de
realizar el cociente (H/V ).
STFT
0

200

Tiempo seg.

400

600

800

1000

1200
5

10
15
Frecuencia Hz

20

25

Figura 3.17.: STFT Aplicada a un registro de microvibracin

50

3.2 Tcnicas Modernas de Estimacin de Espectros


3.2.1.2.

Modelado Paramtrico (AR)

Los mtodos paramtricos extrapolan los valores de una autocorrelacin para retardos de n mayores o iguales a N . Esta extrapolacin es posible si se tiene informacin
sobre la generacin de los datos. As, el primer paso a seguir es elegir un modelo
apropiado para el proceso. Esta eleccin puede estar basada en conocimiento a priori
de cmo ha sido generado el proceso, o en resultados experimentales que indiquen
que un modelo en particular funciona bien [EPo]. Algunos modelos paramtricos
para la estimacin espectral son:
Autorregresivo (AR) con funcin de transferencia todo-polos.
De media mvil (MA) con funcin de transferencia todo-ceros.
Autorregresivo de media mvil (ARMA) con funcin de transferencia con polos
y ceros.
Los modelos autorregresivos pueden describirse como aqullos en los que una variable
o conjunto de variables se explican, al menos en parte, en funcin de los valores
pasados de esa misma variable o conjunto de variables.
Los modelos autorregresivos se abrevian con la palabra AR tras la que se indica
el orden del modelo: AR(1), AR(2),....etc. El orden del modelo expresa el nmero
de observaciones retrasadas de la serie temporal analizada que intervienen en la
ecuacin [dAM]. As, por ejemplo, una serie Y (t) sigue un modelo autorregresivo de
primer orden o AR(1), si es generada por:

Y (t) = c + 1 Yt1 + (t)

(3.10)

Donde c y 1 < < 1 son cosntantes a dterminar y (t) es un proceso de riudo


blanco con varianza 2 .
Este tipo de modelo explica el valor tomado en la serie por la variable en cada
momento de tiempo(t), (t), como debido a la influencia:
1. Del valor tomado por la variable en el momento de tiempo inmediatamente
anterior (t 1), Yt1 . El coeficiente 1 representa el grado de influencia del
valor anterior de la variable sobre el valor considerado. Este coeficiente, al ser
la serie estacionaria, tomar valores entre -1 y 1.
2. Del valor del ruido, 1 , elemento aleatorio, en el momento (t) en que la variable
toma su valor.

51

3.2 Tcnicas Modernas de Estimacin de Espectros


La expresin genrica de un modelo autorregresivo, no ya de un AR(1) sino de un
AR(p) sera la siguiente:

Y (t) = c + 1 Yt1 + 2 Yt2 + ... + p Ytp + (t)

(3.11)

En este modelo general, el valor de la variable en cada momento (t) viene influido
por los valores tomados por la variable en los p momentos anteriores (t 1, t 2,
t 3, ..., t p) y el ruido que se produce en ese momento (t). El grado de influencia
de cada uno de los valores anteriores de la variable viene dado por los parmetros 1 ,
2 , 3 , ..., p . La condicin de estacionaridad impone a cada uno de estos coeficientes
tomar valores entre -1 y 1.
El modelo AR es el ms empleado. Por un lado, es adecuado para representar espectros con peak estrechos (resonancias). Por otro, produce ecuaciones lineales muy
simples para los parmetros AR. El proceso autorregresivo puede obtenerse filtrando
ruido blanco de varianza unidad con un filtro todo-polos.
La Ecuacin 3.10, se puede expresar de la siguiente forma:

Y (t) =

p
X

(k)Y (t k) + (t) para todo t

(3.12)

K=1

Este sistema usa las p salidas ms recientes y la entrada actual para generar recursivamente la prxima salida. Estos modelos adquieren especial significancia en
procesamiento de seales e imgenes debido a que poseen varias propiedades importantes.
En el caso anterior los coeficientes (K) del modelo AR no dependen del tiempo y
por ende tanto los polos como la funcin de transferencia resultantes son fijos para
durante todo el registro. Esto se presta bien cuando el proceso es invariante en el
tiempo. Sin embargo en este estudio se trabaj con seales que varan en el tiempo.
Es por esta razn que en este trabajo los coeficientes AR del modelo son variantes en
el tiempo. Esto quiere decir que se van actualizando a medida que el tiempo avanza
lo que genera una funcin de transferencia y polos variantes en el tiempo.
En el modelo autorregresivo no estacional general AR(p), la cantidad de valores
anteriores de los cuales recibe influencia o depende el valor de la variable en cada
momento viene dado por p. Para un modelo autorregresivo concreto, si cada valor
solo depende del anterior, en el modelo general tenemos p = 1, si de los dos anteriores, p = 2, y as sucesivamente. A partir del modelo general, para abarcar a todos los

52

3.2 Tcnicas Modernas de Estimacin de Espectros


modelos particulares, se obtienen estos modelos concretos. En la prctica, el orden
de p no suele ser muy alto, debido al debilitamiento de la influencia con el aumento
del tiempo, o los momentos, entre valores de la variable.
Los modelos autorregresivos no estacionales pueden interpretarse como una regresin en la cual el valor de la variable se explica a travs del valor tomado por ella
misma en momentos anteriores recientes. De aqu el nombre de auto-regresivos
no-estacionales dado a los modelos AR(p).
Este modelo, se acompaa de una serie de restricciones conectadas con importantes
hiptesis analticas [dAM]:
El proceso no debe ser anticipante (hiptesis de recursividad temporal); lo que
quiere decir que los valores de una variable en un momento t no dependern
de los que esta misma tome en t+j.
El proceso ha de ser invertible, el supone que la correlacin entre una variable
y su pasado va reducindose a medida que nos alejamos ms en el tiempo
(proceso ergdico).
La magnitud de los coeficientes est limitada en valor absoluto: as, por ejemplo, en el caso de un AR(1), el coeficiente autorregresivo de un proceso estocstico estacionario ha de ser inferior a 1 en valor absoluto; en el caso de un
Ar(2), es la suma de los dos coeficientes la que no puede exceder la unidad.
Estas restricciones expresadas en los coeficientes conectan con las propiedades
de estacionariedad del proceso analizado o, dicho de otro modo: slo los modelos cuyos coeficientes respetan una serie de condiciones (que dependen del
orden p del modelo) representan procesos estocsticos estacionarios y, por
tanto, tienen utilidad analtica.
Este modelo se caracteriza por tener un buen comportamiento en dinmicas de larga
duracin, no as en seales con muchas dinmicas trancientes, es aqu donde aparece
el uso de las Wavelet, la cual proporciona un anlisis de resolucin mltiple (en
tiempo y frecuencia), a diferentes escalas y resoluciones, solucionando as el problema
con los transientes. Para el caso de nuestro estudio no se utilizaron Wavelet, por lo
cual no sern analizadas con mayor profundidad.

53

4. Captulo IV MATERIALES Y
MTODOS
4.1.

Base de datos de las seales de ruido ambiental

La informacin que se utiliz en esta investigacin proviene de estudios de microfizacin ssmica que fueron realizados en el sur de Chile, en donde se aplicaron la
tcnica de la razn espectral H/V para la determinacin de la frecuencia natural
de oscilacin del suelo. El estudio fue realizado por el grupo de tesistas de riesgo
sismico del Instituto de Obras Civiles de la Universidad Austral de Chile, con el
objetivo de obtener periodos fundamentales de vibracin y una estimacin de las
amplificaciones locales en los terrenos. Adems de generar mapas de isoperiodos y
efectos de sitio en las zonas en donde se llev a cabo. Esto permiti construir una
gran base de datos con registros de microvibraciones ssmicas, la cual fue tomada
para nuestro trabajo con el fin de comparar las frecuencias naturales encontradas
por Obras Civiles con otras tcnicas de estimacin espectral.
Los registros de ruido ambiental pertenecientes al estudio de Obras Civiles fueron
tomados con frecuencias de muestreo de 50 Hz y 100 Hz con gefonos Tromino y
Geosig (Subseccin A.1.1), en las localidades de Purn, Valdivia, Pitrufquen, Lumaco, Capitn Pastene y Mininco. En la Tabla 4.1 se aprecia un resumen con la
totalidad de seales en cada sector con su respectiva frecuencia de muestreo y el
equipo con el cual fue tomado el registro.
Sector

N de seales muestreadas a 50 Hz

N de seales muestreadas a 100 Hz

Equipo de adquisicin

Pitrufquen

14

Geosig

Capitn Pastene

Geosig

Mininco

Geosig

Purn

Geosig

Valdivia

60

Geosig
Geosig y Tromino

Lumaco

Total

42

60

Cuadro 4.1.: Frecuencias de muestreo

54

4.2 Software Utilizados


Como se puede ver en la Tabla 4.1 las seales muestreadas a 100 Hz fueron tomadas
en la ciudad de Valdivia, mientras que en el resto de las localidades se adquirieron
seales cuya frecuencia de muestreo es de 50 Hz. Los registros que poseen un mayor
nivel de orden y estructuracin en lo que se refiere a protocolo de adquisicin de
datos son las que fueron tomadas en Valdivia esto se puede revisar con mayor detalle
en la Seccin A.1. Todo este conjunto de datos de microfizacin ssmica forman la
base de datos en la cual se trabaj para realizar este estudio.

4.2.

Software Utilizados

Los software utilizados para el procesamiento de las seales fueron principalmente


dos, MatLab y Geopsy a continuacin se entrega una resea de estos programas

4.2.1.

MatLab

Es una herramienta de software matemtico que ofrece un entorno de desarrollo integrado (IDE) con un lenguaje de programacin propio (lenguaje M). Est disponible
para las plataformas Unix, Windows y Mac OS X.
Entre sus prestaciones bsicas se hallan: la manipulacin de matrices, la representacin de datos y funciones, la implementacin de algoritmos, la creacin de interfaces
de usuario (GUI) y la comunicacin con programas en otros lenguajes y con otros
dispositivos hardware. El paquete MATLAB dispone de dos herramientas adicionales que expanden sus prestaciones, a saber, Simulink (plataforma de simulacin
multidominio) y GUIDE (editor de interfaces de usuario - GUI). Adems, se pueden
ampliar las capacidades de MATLAB con las cajas de herramientas (toolboxes); y
las de Simulink con los paquetes de bloques (blocksets).
Es un software muy usado en universidades y centros de investigacin y desarrollo.
En los ltimos aos ha aumentado el nmero de prestaciones, como la de programar
directamente procesadores digitales de seal [lW].

4.2.2.

Geopsy

El equipo Geopsy es el encargado del desarrollo, distribucin y mantenimiento del


software de cdigo abierto para la investigacin y aplicacin geofsica. Nacido durante el proyecto europeo SESAME, ha proporcionado las herramientas para el procesamiento de vibraciones ambientales con caracterizacin de sitio. Progresivamente,

55

4.3 Protocolo de orden y procesamiento para cada seal de la base de datos


se incluyen las tcnicas ms convencionales, con el fin de ofrecer una alta calidad,
integral y una plataforma libre para la interpretacin de los experimentos geofsicos.
Esta agrupacin cuenta con una pgina web [SESa], en la cual se encuentra una
seccin con toda la documentacin necesaria, para apoyar a los usuarios del software,
dicha documentacin se estructura en torno a tutoriales y manuales de referencia.
Los primeros estn destinados a describir paso a paso algunas de las tareas comunes
de procesamiento de vibracin ambiental. La segunda proporciona explicaciones ms
detalladas de las caractersticas del software [SESb].

4.3.

Protocolo de orden y procesamiento para cada


seal de la base de datos

Para comenzar el trabajo de procesamiento de los datos es necesario la conversion


de los archivos entregados por Obras Civiles, estos estan con extensiones nativas de
los equipos de mediciones de ruido ambiental, las nuevas extensiones deben estar en
archivos de texto, todo esto es realizado para que sea mas facil el trabajo posterior,
los detalles de las extensiones y la forma de nombrar a los nuevos archivos de datos
se puden revisar en Apndice B.
El segundo paso dentro del protocolo de trabajo es con el software Geopsy, en donde
se cargan los archivos creados en la etapa anterior para realizar cuatro pruebas
Apndice B en las cuales se varian los parametros de procesado, y se obtienen ademas
las tres componentes de la seal (E-O, N-S, V) con una configuracion de parametros
de proesado, resultando al termino de esta estapa un total de siete archivos.
El tercer bloque del protocolo se realiza en MatLab, es aqu donde se cargan todos
los archivos anteriores para ir creando las figuras, esto se realiza para poder hacer
comparaciones entre los distintos espectros creados en las etapas anteriores y el
estudio de Obras Civiles, ademas se calculan espectros con tecnicas modernas de
estimacion frecuencial, para luego realizar figuras con los comandos mesh, imagesc.
Finalmente se crea un archivo .txt con las observaciones que se vieron durante el
procesado de la seal, esto lo podemos ver en la Figura 4.1.

Figura 4.1.: Etapas Generales del procesamiento de datos

56

5. Captulo V RESULTADOS
FINALES
A continuacin se presentaran los resultados finales de nuestro estudio, la forma de
exponer estos sera a travs de tablas, en primer lugar presentaremos la totalidad
de las seales que fueron procesadas (aproximadamente cien seales), los resultados
estn divididos por sectores a los cuales pertenecen los registros. Ms adelante se
detallaran los principales problemas que se observaron durante el procesado dando
un ejemplo grfico y luego el listado de seales que cumple con esas caractersticas.
Ahora se mostrara la primera tabla en la cual se detalla la nomenclatura con que se
trabajar de aqu en adelante.
Nomenclatura

Descripcin

Ri C. Pastene

Corresponde a Registro numero i de Capitn Pastene

Ri Purn

Corresponde a Registro numero i de Purn

Ri Lumaco

Corresponde a Registro numero i de Lumaco

Ri Pitrufqun

Corresponde a Registro numero i de Pitrufqun

Ri Mininco

Corresponde a Registro numero i de Mininco

Ri Valdivia

Corresponde a Registro numero i de Valdivia

Frec. Natural de EOC. (Hz)

Frec. Natural determinada por el Estudio de Obras Civiles

Frecuencia Muestreo (Hz)

Frecuencia a la cual fue muestreada la seal

Frec. Natural Geopsy. (Hz)

Frecuencia que se determino con el programa Geopsy

Frec. Natural AR (Orden 100 y 500)

Frecuencia natural determinada con modelamiento Auto Regresivo

R5 NL Valdivia

Corresponde a Registro i Nueva Lectura de Valdivia

Sin Determinar

Corresponde a seal con espectro muy complejo para determinar frecuencia

eoc

Equipo de Obras Civiles

Cuadro 5.1.: Nomenclatura de las tablas que poseen el detalle de los registros de
microfizacin ssmicos

A continuacin se proceder a mostrar el detalle de las seales procesadas, de acuerdo


al sector y segn las consideraciones anteriores, en la Tabla 5.7, Tabla 5.8, Tabla 5.3,
Tabla 5.2, Tabla 5.4, Tabla 5.5 y Tabla 5.6.
La primera columna donde aparece la palabra registro, corresponde a los nombres de todas las muestras tomadas en la ciudad de Valdivia

57

Captulo V RESULTADOS FINALES


En la segunda columna encontramos la frecuencia de muestreo con la cual fue
adquirido el registro
En la tercera columna tenemos la frecuencia natural estimada por el equipo
de Obras Civiles, esta frecuencia estaba en la base de datos y fue estimada con
tcnicas usadas normalmente en el rea de microfizacin ssmica
En la cuarta columna encontramos la frecuencia natural estimada por nuestro
estudio, ocupando el software Geopsy
En la ultima columna tenemos la frecuencia natural de vibracin estimada por
nuestro estudio con modelado AR, en esta columna se aprecia que en algunas
seales aparecen dos frecuencias principales por eso existen como frecuencias
determinadas dos posibilidades. Mientras que en otras seales las frecuencias
aparecen separadas por un guion, lo cual indica que es una ancho de banda
frecuencial y que se encuentra entre las frecuencias que aparecen sealadas en
el cuadro correspondiente. La otra posibilidad es que la naturaleza de la seal
sea muy compleja, razn por la cual es muy difcil determinar una frecuencia,
en tales casos el cuadro dice sin determinar. Para mayor detalle sobre los
parmetros de procesamiento se pueden revisar en el Apndice B.
Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R1 Pitrufqun

50

2,087

2,07

2,1

R2 Pitrufqun

50

3,149

3,06

3,02

R3 Pitrufqun

50

2,722

3,03

3,02

R4 Pitrufqun

50

2,368

2,5

2,65 y 3,41

R5 Pitrufqun

50

3,479

3,3

3,37

R6 Pitrufqun

50

2,588

2,63

2,62

R7 Pitrufqun

50

3,101

2,9

R8 Pitrufqun

50

2,576

2,641

2,85

R9 Pitrufqun

50

2,637

2,372

2,75

R10 Pitrufqun

50

2,612

3,35

3,3

R11 Pitrufqun

50

1,965

2,12

2,16

R12 Pitrufqun

50

1,282

1,21

1,64 y 2,39

R13 Pitrufqun

50

2,502

2,233

2,91

R14 Pitrufqun

50

2,832

2,54

2,35-3,4

Cuadro 5.2.: Seales Procesadas de Pitrufqun

58

Captulo V RESULTADOS FINALES


Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR (Hz)

R1 C. Pastene

50

2,747

2,45

2,71

R2 C. Pastene

50

2,393

3,71

R3 C. Pastene

50

2,673

2,349

2,56

R4 C. Pastene

50

2,087

2,259

2,35

R5 C. Pastene

50

2,466

2,12

2,65

Cuadro 5.3.: Seales Procesadas de Capitn Pastene

Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R1 Purn

50

3,43

3,7

3,71

R2 Purn

50

3,308

3,4

3,54

R3 Purn

50

5,115

4,5

2,67

R4 Purn

50

3,247

3,1

5,21

R5 Purn

50

5,103

2,7

Entre 2,67-3,65 y 5,5

Cuadro 5.4.: Seales Procesadas de Purn

Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R1 Lumaco

50

3,418

3,6

3,72

R2 Lumaco

50

5,603

4,69

4,62 y 7,22

R3 Lumaco

50

8,435

7,06 y 8,09

Entre 5,33 y 14,66

R4 Lumaco

50

4,077

3,7

Entre 3,19 y 4,75

R5 Lumaco

50

10,205

10,05

Entre 6,94 y 13,34

Cuadro 5.5.: Seales Procesadas de Lumaco

Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R1 Mininco

50

1,257

1,232

2,45

R2 Mininco

50

1,05

1,12 y 1,356

1,54;2,4;3,85

R3 Mininco

50

1,221

1,048

3,5

R4 Mininco

50

1,221

1,101

1,3-1,6

R5 Mininco

50

0,989

1,116

1,12-1,5

R6 Mininco

50

11,328

0,857 y 11,68

1,07 y 10,18

R7 Mininco

50

0,952

0,95

1,09-1,5

Cuadro 5.6.: Seales Procesadas de Mininco

59

Captulo V RESULTADOS FINALES


Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R58 Valdivia

100

1,221

1,273

1,2-1,9

R59 Valdivia

100

1,062

1,097

1,2-1,8

R60 Valdivia

100

0,964

1,24-2,34;9,7

R201 Valdivia

100

1,147

1,157

0,98-2,2

R202 Valdivia

100

0,989

1,123

1,3-2

R203 Valdivia

100

1,306

1,224 y 1,382

1,4-2,3

R204 Valdivia

100

1,379

1,457

1,31-2,44

R206 Valdivia

100

0,952

0,9 y 0,953

0,88-1,7

R208 Valdivia

100

0,952

0,974

0,94-1,9

R209 Valdivia

100

0,757

0,769

0,98-1,8

R210 Valdivia

100

1,318

0,726 y 1,284

1,18-2

R211 Valdivia

100

0,793

0,799

0,8-1,8

R213 Valdivia

100

0,928

0,91

1,25-3

R214 Valdivia

100

1,013

0,95-2,2

R215 Valdivia

100

1,88

1,068

1,4-2,5

R222 Valdivia

100

0,769

0,778

0,96-1,72

R701 Valdivia

100

0,732

0,77

Entre 6,48 y 10,25

R702 Valdivia

100

0,732

0,74

Entre 3,3-5,44 y 6,1-7,31

R703 Valdivia

100

0,793

0,8

4,27

R704 Valdivia

100

0,928

0,91

5,65

R705 Valdivia

100

0,72

0,75

Entre 5,1 y 8,69

R706 Valdivia

100

0,72

0,7

0,74 y 1,3

R707 Valdivia

100

0,781

0,8

1,78

R708 Valdivia

100

0,732

0,74

Sin Determinar

R709 Valdivia

100

1,318

1,41

2,85

R710 Valdivia

100

0,732

0,71

5,79

R711 Valdivia

100

0,745

0,71

5,82

R712 Valdivia

100

0,818

0,81

5,62 y 5,82

R1 Valdivia

50

0,928

0,902

1,6-2

R2 Valdivia

50

0,94

0,924

0,96-1,8

R3 Valdivia

50

0,964

0,89

1,07-1,7

R4 Valdivia

50

0,769

0,77

1,05-1,95

R5 NL Valdivia

100

0,708

0,703

0,5-1,4

R6 Valdivia

50

0,696

0,935

1,6-2

Cuadro 5.7.: Seales Procesadas de Valdivia

60

Captulo V RESULTADOS FINALES


Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R7 Valdivia

100

1,331

1,007

1,78

R8 Valdivia

100

1,135

1,056

Sin Determinar

R9 Valdivia

100

1.001

0,988

1,3 y 1,81

R11 Valdivia

100

1,135

1,13

1,64 y 2,19

R12 Valdivia

100

1,16

1,2

1,6-2,26

R16 Valdivia

100

1,05

1,123

1,48-2,27

R18 Valdivia

100

1,16

1,13

Entre 1,18 y 1,94

R20 Valdivia

100

1,135

1,15

1,68 y 3,65

R22 Valdivia

100

1,587

1,71

Entre 1,5-2,23 y 3,78

R24 Valdivia

100

1,184

1,05

2,54

R26 Valdivia

100

0,952

0,98

1,71

R27 Valdivia

100

0,903

1,04

1,64

R28 Valdivia

100

1,196

1,198

1,06-1,77

R29 Valdivia

100

6,03

5,86

5,84

R31 Valdivia

100

0,671

0,733

0,92

R32 Valdivia

100

0,757

0,748

1,24-1,76

R33 Valdivia

100

0,781

0,759

1,2-1,8

R34 Valdivia

100

0,769

0,782

1,34-1,85

R35 Valdivia

100

0,781

0,808 y 0,984

1,36-3,23

R36 Valdivia

100

0,781

0,748

1,11-1,71

R37 Valdivia

100

0,781

0,785

0,99-1,6

R38 Valdivia

100

0,818

0,868

1,16-2,24

R41 Valdivia

50

1,184

1,093

1,4-2,4

R42 Valdivia

100

1,086

1,022

1,5-3

R43 Valdivia

100

0,879

1,018

1,5-2,5

R45 Valdivia

100

1,245

1,38

1,43 y 2,23

R48 Valdivia

100

1,062

0,98

Entre 1,3 y 2,06

R50 Valdivia

100

1,147

1,4

Entre 1,5 y 2,19

R51 Valdivia

100

1,257

1,22

Sin Determinar

R52 Valdivia

100

1,135

1,24

2,78; 4,68 y 5,75

R55 Valdivia

100

0,83

0,83

1,5

R56 Valdivia

100

0,842

0,823

1,2-1,6

Cuadro 5.8.: Seales Procesadas de Valdivia

A modo de presentacin en la Tabla 5.9 podemos encontrar la totalidad de seales


procesadas de acuerdo a cada sector y en la Figura 5.1 podemos observar los porcentajes que ocupan cada uno de los sectores dentro del total de seales procesadas
de la base de datos.

61

Captulo V RESULTADOS FINALES


Sector
Cantidad de seales
Valdivia
66
Capitn Pastene
5
Pitrufqun
14
Purn
5
Lumaco
5
Mininco
7
Total
102
Cuadro 5.9.: Distribucin de seales procesadas en cada localidad

Figura 5.1.: Distribucin porcentual de todas las seales procesadas en nuestro estudio

En las siguientes secciones se mostraran las situaciones particulares observadas en


las distintas experiencias realizadas en este trabajo. A modo de clarificar los ttulos
de cada seccin se presenta una breve descripcin de cada uno.
Primera situacin observada: corresponde a una revisin de los 4 experimentos realizados con Geopsy. Aqu se agrupan los casos en que dos experimentos muestran una frecuencia, mientras que los otros dos muestran una
frecuencia distinta.
Segunda situacin observada: corresponde a la situacin en donde uno
de los 4 experimentos realizados en Geopsy, no coincide en frecuencia con los
dems ni con la frecuencia de eoc.
Tercera situacin observada: corresponde a la situacin en la cual los 4
experimentos poseen un peak truncado en relacin al peak obtenido por eoc.

62

5.1 Primera Situacin Observada

5.1.

Primera Situacin Observada

En base a las observaciones realizadas en nuestro estudio, durante la etapa de procesamiento con el software Geopsy. Y utilizando las cuatro configuraciones de parmetros definidas para este propsito. Se encontr que en ocasiones, dos conjuntos de
parmetros coincidan en la estimacin de la misma frecuencia natural, mientras que
los dos restantes, coincidan en la estimacin de una frecuencia natural distinta. La
siguiente situacin la podemos observar en la Figura 5.2, donde las pruebas (H/V)2
y (H/V)3 muestran una frecuencia natural estimada de 1,36 Hz, mientras que (H/V)1
y (H/V)4 muestran una frecuencia natural estimada de 1,11 Hz.
4
eoc
H/V 1
H/V 2
H/V 3
H/V 4

3.5

1,05 Hz
3

1,37 Hz

2.5

1.5

1,36 Hz
1

1,11 Hz
0.5

0.5

1.5

2.5

3.5

4.5

Figura 5.2.: Figura donde se aprecia que dos espectros de Geopsy estiman una
frecuencia natural, mientras que las otras dos pruebas restantes estiman otra frecuencia

En la Tabla 5.10 observamos el listado completo de seales con esta condicin.

63

5.2 Segunda Situacin Observada


Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R2 Mininco

50

1,05

1,12 y 1,356

1,54;2,4;3,85

R3 Mininco

50

1,221

1,048

3,5

R4 Mininco

50

1,221

1,101

1,3-1,6

R5 Mininco

50

0,989

1,116

1,12-1,5

R12 Pitrufqun

50

1,282

1,21

1,64 y 2,39

R18 Valdivia

100

1,16

1,13

Entre 1,18 y 1,94

R35 Valdivia

100

0,781

0,808 y 0,984

1,36-3,23

R48 Valdivia

100

1,062

0,98

Entre 1,3 y 2,06

R50 Valdivia

100

1,147

1,4

Entre 1,5 y 2,19

R51 Valdivia

100

1,257

1,22

Sin Determinar

R52 Valdivia

100

1,135

1,24

2,78; 4,68 y 5,75

R203 Valdivia

100

1,306

1,224 y 1,382

1,4-2,3

R209 Valdivia

100

0,757

0,769

0,98-1,8

Cuadro 5.10.: Lista de seales donde la frecuencia estimada por Geopsy coincide
con otro conjunto de parmetros de procesado

Que como podemos ver en la Figura 5.3 un 13 % del total de las seales, presentaron
la situacin descrita.

Figura 5.3.: Porcentaje de seales con la primera situacin

5.2.

Segunda Situacin Observada

Aqu agruparemos un segundo tipo de situacin observada, durante el trabajo de


procesado. Ya que se pudo ver en algunas ocasiones que una prueba de Geopsy no se
asemejaba a ninguna de las otras tres, ni tampoco a el estudio realizado por Obras
Civiles, esta situacin queda graficada en la Figura 5.4. Donde el espectro de (H/V)3
por ejemplo no sigue la forma que describen (H/V)1 , (H/V)2 , y (H/V)4 .

64

5.2 Segunda Situacin Observada


4
eoc3
H/V1
H/V2
H/V3
H/V4

5,32 Hz

3.5

1,81 Hz
8,37 Hz

2.5

1.5
4,61 Hz
1

0.5

10

15

Figura 5.4.: Espectro de una prueba realizada en Geopsy no sigue la forma de los
otros espectros realizados con otros conjuntos de parmetros

En la Tabla 5.11 vemos en detalle cuales son los registros que poseen esta condicin, luego de haberlas procesado con Geopsy y en la Figura 5.5 podemos ver que
aproximadamente un 3 % del total de las seales presentaron esta situacin.
Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R3 Lumaco

50

8,435

7,06 y 8,09

Entre 5,33 y 14,66

R6 Valdivia

50

0,696

0,935

1,6-2

R7 Valdivia

100

1,331

1,007

1,78

Cuadro 5.11.: Resumen de seales que poseen un espectro diferente a las otras
pruebas de Geopsy

65

5.3 Tercera Situacin Observada

Figura 5.5.: Porcentaje de seales que poseen un espectro diferente a las dems
pruebas de Geopsy

5.3.

Tercera Situacin Observada

Aqu veremos todas las seales, cuyos espectros realizados con los cuatro conjuntos
de parmetros en Geopsy muestran un truncamiento del peak principal. En comparacin al peak mas fino y definido que muestra eoc, la situacin queda plasmada
en la Figura 5.6 correspondiente al registro Lumaco 4. Donde el cuadro negro nos
muestra la diferencia entre ambos espectros.
6
eoc
H/V 1
H/V 2
H/V 3
H/V 4

2,39 Hz

5
3,05 Hz
4

2,98 Hz

3,52 Hz

10

15

Figura 5.6.: Truncamiento de los espectros de Geopsy en comparacin al peak ms


definido de eoc
En la Tabla 5.12 encontramos en detalle todas las seales que presentan esta situa-

66

5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy
cin y en la Figura 5.7 podemos ver grficamente que aproximadamente un 10 % del
total de las seales procesadas durante el estudio, presentaron esta condicin.
Registro

Frecuencia Muestreo (Hz)

Frec. Natural de EOC. (Hz)

Frec. Natural Geopsy. (Hz)

Frec. Natural AR. (Hz)

R2 C. Pastene

50

2,393

3,71

R5 C. Pastene

50

2,466

2,12

2,65

R2 Lumaco

50

5,603

4,69

4,62 y 7,22

R4 Lumaco

50

4,077

3,7

Entre 3,19 y 4,75

R4 Pitrufqun

50

2,368

2,5

2,65 y 3,41

R9 Pitrufqun

50

2,637

2,372

2,75

R10 Pitrufqun

50

2,612

3,35

3,3

R13 Pitrufqun

50

2,502

2,233

2,91

R3 Valdivia

50

0,964

0,89

1,07-1,7

R8 Valdivia

100

1,135

1,056

Sin Determinar

Cuadro 5.12.: Lista de seales que presentan la condicin de peak truncado en


comparacin al peak mas definido de eoc

Figura 5.7.: Porcentaje de seales con peak truncado en comparacin al peak mas
definido de eoc

5.4.

Estudio estadstico comparativo entre los


resultados de Obras Civiles y Geopsy

Para verificar el porcentaje de similitud entre los estudios realizados por Obras
Civiles y nuestro estudio ocupando el software Geopsy se realizaron tres histogramas, observando las frecuencias que se muestran en Tabla 5.7, Tabla 5.8, Tabla 5.3,
Tabla 5.2, Tabla 5.4, Tabla 5.5 y Tabla 5.6, realizando para ello la diferencia entre
las frecuencias de eoc y Geopsy.

67

5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy

45
40
35
30
25
20
15
10
5
0
0

0.5

1.5

2.5

Figura 5.8.: Histograma 1

Debido a que en algunos casos se presentaban dos frecuencias observadas, se realizaron tres histogramas, la Figura 5.8 muestra la comparacin con la frecuencia ms
baja de los registros con dos frecuencias, los resultados fueron que 93 seales poseen
un error menor a 0,5 Hz, lo cual representa un 91,2 % de las seales analizadas.
Otras 6 se encuentran con un error de entre 0,5 y 1 Hz, es decir un 5,8 % del total,
2 seales con error de 1 a 3 Hz, en porcentaje un 1,96 % y 1 seal que esta fuera de
rango con una diferencia de 10,5 Hz.
50
45
40
35
30
25
20
15
10
5
0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Figura 5.9.: Histograma 1 detallado (a).

68

5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy

3.5

2.5

1.5

0.5

10

11

Figura 5.10.: Histograma 1 detallado (b).


Para tener una mejor apreciacin de los datos obtenidos en los histogramas, se realizo
un anlisis en intervalos ms pequeos. Considerando que un error menor a 0,2 Hz
no es relevante, entre 0,2 y 0,5 Hz sera aceptable y un error mayor a 0,5 Hz, ya
es algo que habra que revisar detalladamente. En la Figura 5.9 y la Figura 5.10 se
aprecia que: 81 seales presentan error bajo los 0,2 Hz, correspondiente a un 79,41 %,
12 seales poseen error entre 0,2 y 0,5 Hz, un 11,76 % y 9 seales presentan error
mayor a 0,5 Hz, correspondiente a un 8,82 % de un total de 102 seales procesadas.

35

30

25

20

15

10

0
0

0.5

1.5

2.5

Figura 5.11.: Histograma 2


En el segundo Histograma, Figura 5.11, se utilizaron las frecuencias ms altas de los
registros con frecuencias dobles, los resultados observados aqu fueron, un total de

69

5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy
95 seales con error menor a 0,5 Hz, un 93,14 % del total, 5 seales entre 0,5 y 1 Hz
de error, un 4,9 % y una sola seal con error entre 1 y 3 Hz, un 0,98 % del total.
35

30

25

20

15

10

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Figura 5.12.: Histograma 2 detallado (a)

3.5

2.5

1.5

0.5

0
0.5

1.5

2.5

Figura 5.13.: Histograma 2 detallado (b)

En mayor detalle, como se muestra en la Figura 5.12 y la Figura 5.13, si se divide en


intervalos ms pequeos, se observan 79 seales con error menor a 0,2 Hz, 77,45 %,
17 seales con error entre 0,2 y 0,5 Hz, 16,66 % y 6 seales con error mayor a 0,5
Hz, 5,88 %.

70

5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy

60

50

40

30

20

10

0
0

10

Figura 5.14.: Histograma 3

El tercer histograma, Figura 5.14, se realiz, observando los espectros de las seales
y eligiendo el peak con mayor amplitud, lo cual arrojo los siguientes resultados, 94
seales con error menor a 0,5 Hz, lo cual representa un 92,1 % del total de las seales,
6 seales con errores entre 0,5 y 1 Hz, 5,88 % del total, y 2 seales con error sobre
1 Hz, un 1,96 % del total de las seales.
70

60

50

40

30

20

10

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Figura 5.15.: Histograma 3 detallado (a)

71

5.5 Situaciones Observadas con Tcnicas Modernas

3.5

2.5

1.5

0.5

10

11

Figura 5.16.: Histograma 3 detallado (b)

Para mayor detalle, dividido en intervalos mas pequeos, como se muestra en la


Figura 5.15 y la Figura 5.16, se tienen 80 seales con error menor a 0,2 Hz, correspondiente a un 78,43 %, 14 seales con errores entre 0,2 y 0,5 Hz, correspondiente
a un 13,73 %, y 8 seales con error mayor a 0,5 Hz, correspondiente al 7,84 % del
total de las seales procesadas.
Finalmente con estos resultados podemos decir que cerca de un 90 % de las seales
posee una alta aceptabilidad de resultados, al encontrar una diferencia menor a 0,5
Hz, presentndose algunos caso particulares debidos a complejidades en el anlisis
frecuencial, con intervencin de ruidos impulsionales demasiado grandes en el registro, que provocan alteraciones en el espectro, produciendo errores a la hora de elegir
la frecuencia natural de un registro. En resumen ms detallado y promediando los
resultados de los 3 histogramas, tenemos que alrededor de un 78,43 % de las seales
poseen un error menor a 0,2 Hz, un 13,88 % presenta un error entre 0,2 y 0,5 Hz y
un 7,51 % presenta el error ms considerable, sobre los 0,5 Hz.

5.5.

Situaciones Observadas con Tcnicas Modernas

En esta parte del anlisis, se utilizaron tcnicas modernas de estimacin espectral,


como STFT y modelamiento auto regresivo variante en el tiempo. Esto con el propsito de tener otros puntos de vista sobre el espectro resultante de las seales. Y
tambin observar las diferencias obtenidas en la estimacin frecuencial, entre las
tcnicas utilizadas normalmente en el rea de microfizacin ssmica y estas nuevas
tcnicas. A continuacin se presentan una serie de registros interesantes, los cuales
fueron tratados con las tcnicas anteriormente mencionadas.

72

5.5 Situaciones Observadas con Tcnicas Modernas

5.5.1.

Registro correspondiente a Pitrufqun 10

El anlisis frecuencial de este registro nos muestra un interesante caso, el cual tiene
como foco de observacin las diferencias obtenidas en la estimacin de la frecuencia
natural obtenidas por el equipo de Obras Civiles y nuestro estudio. Y tambin, de
ver como las dos metodologas hacen una estimacin distinta de los peak de menor
amplitud que muestran como resultado luego del tratamiento frecuencial.
El registro Pitrufquen 10 en el dominio temporal lo encontramos en la Figura 5.17.
Donde se aprecian desde los 0 hasta los 400 segundos una actividad agitada en
relacin, a lo que sucede en el resto del registro.
500
EO
0

500

200

400

600

800

1000

1200

1400

500
NS
0

500

200

400

600

800

1000

1200

1400

2000
V
0
2000
4000

200

400

600

800

1000

1200

1400

Figura 5.17.: Seal Pitrufquen 10 por componentes y en el dominio temporal

Si nos dedicamos a ver el espectro de la seal, tal y como lo muestra la Figura 5.18.
Donde los dos cuadros negros nos indican los peaks predominantes en l, y encerrados
en estos, estn las cuatro pruebas realizadas en Geopsy, que arrojan tres frecuencias
estimadas como resultado, mientras que eoc solo muestra dos frecuencias estimadas.
El cuadro negro de la derecha en la Figura 5.18. Presenta una clara similitud entre
lo que muestra eoc y Geopsy, no destacndose mayores problemas. En cambio en
el cuadro negro de la izquierda de la misma figura, el espectro de eoc solo muestra
un peak, mientras Geopsy muestra dos peaks, y que adems esta situacin cambia
totalmente la eleccin de la posible frecuencia natural de este registro. Las razones
por las cuales se provocara esta diferencia podran ser los ruidos que presenta la
seal en el tiempo, y que cuya presencia excit una frecuencia que los algoritmos
trataron de forma diferente. Ya que mientras uno atenu la frecuencia que esta a los
3,3 Hz el otro la conservo y la entrego como la frecuencia natural estimada.

73

5.5 Situaciones Observadas con Tcnicas Modernas


7
eoc
H/V1
H/V2
H/V3
H/V4

2,6 Hz

2,7 Hz

5,4 Hz

3,3 Hz
1

10

15

Figura 5.18.: Imagen de los espectros de Pitrufqun 10 calculados con eoc y Geopsy

Tambin se compararon resultados con modelamiento AR variante en el tiempo, tal


y como nos muestra la Figura 5.19. En donde se pueden observar con facilidad tres
frecuencias estimadas, las cuales estn demarcadas por tres lineas amarillas y con
sus respectivas frecuencias en el lado derecho de la imagen. Los resultados de esta
prueba se asemejan en mayor medida a lo que nos entrega Geopsy. Ya que la banda
de 3,46 Hz esta mucho mas marcada y definida que la de 2,8 Hz.
0

200

400

2,8 Hz

600

3,46 Hz

800

5,23 Hz
1000

10

15

Figura 5.19.: Modelamiento AR variante en el tiempo de Pitrufqun 10

Gracias a la funcin mesh de MatLab, la cual es una representacin grfica de


superficie en tres dimensiones de una matriz, podemos apreciar, lo que sucede en el
espectro estimado con modelamiento AR variante en el tiempo. En la Figura 5.20

74

5.5 Situaciones Observadas con Tcnicas Modernas


podemos observar que existen tres bandas de frecuencias las cuales estn marcadas
por lineas rojas y numeradas desde izquierda a derecha, en donde tambin existe un
peak muy grande a los 3,02 Hz el cual puede estar siendo excitado por algn ruido
presente en la seal.

Figura 5.20.: mesh del espectro obtenido con modelamiento AR variante en el tiempo

5.5.2.

Registro correspondiente a Purn 5

El registro 5 tomado en la localidad de Purn, nos entrega una muy buena apreciacin sobre la estimacin de la frecuencia natural, cuando existen dos bandas
frecuenciales importantes y muy marcadas en el espectro. De las cuales ninguna es
el resultado de una excitacin provocada por algn ruido particular, y que sin embargo nos hacen preguntarnos como los algoritmos hacen la eleccin de la frecuencia
natural dependiendo del nivel y tipo de ruido presente en ese momento. A continuacin procederemos a revisar en la Figura 5.21 la seal en sus tres componentes y en
el dominio temporal.

75

5.5 Situaciones Observadas con Tcnicas Modernas


2000
E0
0

2000

200

400

600

800

1000

1200

1400

2000
NS
1000
0
1000

200

400

600

800

1000

1200

1400

1000
V
0

1000

200

400

600

800

1000

1200

1400

Figura 5.21.: Registro Purn 5 por componentes y en el dominio temporal

Como podemos corroborar en la Figura 5.22 tanto las pruebas de Geopsy como la
de eoc nos indican la presencia de dos bandas frecuenciales muy marcadas en el
espectro, no destacndose mayores peak a parte de los que se muestran en la figura,
por tal motivo se verifica la existencia de dos frecuencias igualmente importantes y
que predominan en el espectro.
5.5
eoc
H/V1
H/V2
H/V3
H/V4

5,11 Hz

5
3,11 Hz

4.5
4
3.5
3
2.5
2

2,72 Hz

5,08 Hz

1.5
1
0.5

10

15

Figura 5.22.: Anlisis espectral de Purn 5 tanto por Geopsy como por eoc

Para tener una apreciacin distinta del espectro, ocupamos el modelado AR variante
en el tiempo, el cual nos entrega el espectro de la Figura 5.23. Como nos podemos dar
cuenta aqu tambin se aprecian las dos bandas frecuenciales, y como se mencion
antes, ambas bandas son de importancia en el anlisis frecuencial, siendo las dos
posibles candidatas a ser frecuencias naturales estimadas.

76

5.5 Situaciones Observadas con Tcnicas Modernas


2,35Hz3,96Hz 4,36Hz6,73Hz

1000

800

600

400

200

Banda 1

Banda 2

10

15

Figura 5.23.: Modelamiento AR variante en el tiempo de Purn 5

Si observamos el espectro en tres dimensiones, tal y como nos muestra la Figura 5.24
podemos ver con claridad, que la frecuencia de 3,88 Hz, que tambin esta presente
en la Figura 5.23, no aparece en el espectro de la Figura 5.22, esto debido a que los
algoritmos de ambos mtodos desecharon el ruido que origina esta frecuencia.
Lo que sucedera en este caso es que la diferencia entre ambas bandas, sera mayor
o menor dependiendo del ruido presente en ese instante dado, lo cual podra hacer variar la eleccin del peak mas alto por parte de los algoritmos y por ende la
estimacin de la frecuencia natural.

77

5.5 Situaciones Observadas con Tcnicas Modernas

Figura 5.24.: Nakamura AR variante en el tiempo del espectro correspondiente al


registro Purn 5

5.5.3.

Registro correspondiente a Mininco 6

El registro correspondiente a Mininco 6, nos entrega una idea de el trabajo realizado


por los algoritmos en la eleccin de la frecuencia natural, cuando en el suelo existe un
comportamiento en baja y en alta frecuencia. Este registro es de inters en nuestros
resultados porque en los mltiples estudios de Microzonificacin ssmica siempre se
escoge una nica frecuencia natural como resultado, cuando perfectamente es posible
tambin la existencia de mas de una frecuencia natural presente.
Solo que separada a varios Hz de distancia una de la otra, pero igualmente importante que la frecuencia natural escogida. Razn por la cual los cientficos del rea
pueden no estar considerando informacin importante del espectro en un suelo, por
omisin de los softwares que utilizan.
La seal la podemos encontrar en el dominio temporal en la Figura 5.25. Donde
estn sus tres componentes Este-0este, Norte-Sur y Vertical.

78

5.5 Situaciones Observadas con Tcnicas Modernas


5000
EO
0

5000

200

400

600

800

1000

1200

1400

5000
NS
0

5000

200

400

600

800

1000

1200

1400

5000
V
0

5000

200

400

600

800

1000

1200

1400

Figura 5.25.: Registro Mininco 6 por componentes y en el dominio temporal

Para tener una apreciacin de lo que sucede en el espectro de la seal, en la


Figura 5.26 podemos observar de color azul el espectro determinado por el equipo
de Obras Civiles y las cuatro pruebas realizadas en Geopsy con su correspondiente color. Existe un comportamiento frecuencial en baja y alta frecuencia como se
menciono antes, pero el de alta frecuencia es mas disperso.
6
eoc
H/V 1
H/V 2
H/V 3
H/V 4

11,35 Hz

0,84 Hz
0

10

15

Figura 5.26.: Anlisis espectral de Mininco 6 con Geopsy y eoc

Para tener una idea de que es lo que sucede con un anlisis AR variante en el
tiempo, en la Figura 5.27. Encontramos el espectro en donde tambin se puede ver
el comportamiento frecuencial dual que posee el registro, en donde se aprecian tres
bandas frecuenciales, con sus respectivas frecuencias al margen de la imagen.

79

5.5 Situaciones Observadas con Tcnicas Modernas


2,45Hz

4,93Hz5,47Hz

7,59Hz12,44Hz

200

400

600

800

1000

1200

5,73Hz7,08Hz

10

15

Figura 5.27.: Modelado AR variante en el tiempo del registro Mininco 6

Lo que podemos inferir de esta experimentacin, es la falta de atencin que realizan


los algoritmos a los peak que son omitidos, durante el procesado. Ya que solo consideran una frecuencia como resultado, cuando es muy poco probable que en todos
los casos solo exista una frecuencia fundamental.

5.5.4.

Registro correspondiente a Capitn Pastene 2

El anlisis de Capitn Pastene 2, nos lleva a un atractivo caso sobre las repercusiones
que provocara la excitacin de una frecuencia en particular, a causa de un gran ruido
presente en la seal temporal. Y cmo los algoritmos que normalmente se utilizan
en el rea de microfizacin ssmica, enfrentan esta situacin, para poder entregar
una frecuencia natural estimada.
La seal por componentes y en el dominio temporal se presenta en la Figura 5.28.

80

5.5 Situaciones Observadas con Tcnicas Modernas


100
EO
50
0
50

200

400

600

800

1000

1200

1400

100
NS
0

100

200

400

600

800

1000

1200

1400

100
V
0

100

200

400

600

800

1000

1200

1400

Figura 5.28.: Capitn Pastene 2 por componentes y en el dominio temporal

Para ver el espectro por componentes de la seal, ocuparemos un anlisis de tiempofrecuencia, en cada una de ellas. Comenzando por la componente E-O. En donde
podemos ver que existe una frecuencia de 3,69 Hz y un marcado ruido entre los 693
y los 1187 segundos.
Spectogyule orden 100 este
0

200

3,69 Hz

400

600

800

1000

1200

10

15

Figura 5.29.: Tiempo-Frecuencia de la componente E-O de Capitn Pastene 2

Ahora procederemos a revisar la seal en su componente Norte-Sur en la Figura 5.30.


Donde se aprecia nuevamente la frecuencia de 3,69 Hz, mientras que un ruido entre
689 y los 1187 segundos.

81

5.5 Situaciones Observadas con Tcnicas Modernas


0

200

3,69 Hz

400

600

800

1000

1200

10

15

Figura 5.30.: Tiempo- Frecuencia de la componente N-S de Capitn Pastene 2

Y por ultimo tenemos la Figura 5.31 en donde encontramos la componente vertical


de la seal con un anlisis de tiempo-frecuencia. Aqu la frecuencia mas significativa
corresponde a unos 3,71 Hz, y nuevamente un ruido presente entre los 694 y los 1187
segundos.
0

200

3,71 Hz
400

600

800

1000

1200

10

15

Figura 5.31.: Tiempo-Frecuencia de la componente Vertical de Capitn Pastene 2

Para ver el espectro resultante de las tres componentes e identificar cuales son las
frecuencias mas importantes dentro del espectro, observamos la Figura 5.32. Donde
apreciamos que existen al menos cuatro frecuencias claras en el espectro, estas fueron
obtenidas del estudio de Obras Civiles y las cuatro configuraciones de parmetros
de Geopsy.

82

5.5 Situaciones Observadas con Tcnicas Modernas


6

2,39 Hz

eoc
H/V 1
H/V 2
H/V 3
H/V 4

3,05 Hz

2,98 Hz
1

3,52 Hz
0

10

15

Figura 5.32.: Anlisis espectral de Pastene 2 con Geopsy y eoc


Para ver las posibles similitudes o diferencias con lo que nos entrega Geopsy y eoc, se
realiz una ultima prueba, la cual se observa en la Figura 5.33. Donde se entrega el
espectro de Nakamura ocupando modelado AR variante en el tiempo. Es aqu donde
nuevamente se aprecia la frecuencia de 3,69 Hz, algo sumamente extrao ya que
en las figuras Figura 5.29, Figura 5.30 y Figura 5.31. Tambin esta presente, pero
supuestamente al hacer el cociente (H/V) de Nakamura debera haberse eliminado del
espectro resultante. Y por su puesto tambin tenemos en la Figura 5.33 una banda
frecuencial que va desde los 2,19 Hz hasta los 3,36 Hz.
2,19Hz3,36Hz

1000

800

3,69 Hz

600

400

200

10

15

Figura 5.33.: Modelamiento AR variante en el tiempo de Capitn Pastene 2


Ahora si comparamos la Figura 5.32 y la Figura 5.33 nos podemos dar cuenta que
existe un gran ruido que esconde a la frecuencia de 3,69 Hz presente en Nakamura

83

5.6 Eliminacin de ruido transientes


AR y a los 3,05 Hz presente en la Figura 5.32. Y que provocara finalmente que los
algoritmos entreguen como respuesta los 2,39 Hz.

5.6.

Eliminacin de ruido transientes

Aunque nuestro estudio no apunta directamente a determinar que es lo que sucede


en la estimacin espectral, cuando en una misma seal se procesan los datos con
y sin la presencia de un gran ruido transiente. Sin embargo, esta prueba fue parte
de las experimentaciones realizadas durante el procesamiento de las seales, con
el objetivo de tener una idea de lo que sucede con la estimacin de la frecuencia
natural, bajo esta situacin. Para ello se utiliz un registro proveniente de la ciudad
de Valdivia, el cual lleva por nombre Valdivia 709, y cuya frecuencia de muestro
es de 100 Hz. Esta prueba se realiz principalmente en Geopsy. Ocupando cuatro
conjuntos de configuraciones de parmetros para el procesamiento de la seal.
Como se mencion anteriormente existen cuatro conjuntos de parmetros de procesamiento en Geopsy, los cuales se pueden encontrar en el Apndice B. Lo nico
diferente que se presenta en esta seal fue que la primera prueba se realizo con
el ruido transiente presente. Tal como se observa en la Figura 5.34 el ruido esta
encerrado con un cuadro rojo, los resultados luego de las pruebas se muestran en
la Figura 5.35. Donde (H/V)1 , (H/V)2 , (H/V)3 y (H/V)4 representan el espectro resultante con el primer conjunto de parmetros y as sucesivamente con el resto de los
cocientes.
5000
EO
0

5000

200

400

600

800

1000

1200

1400

1600

1800

2000

5000
NS
0

5000

200

400

600

800

1000

1200

1400

1600

1800

2000

5000
Vertical
0

5000

200

400

600

800

1000

1200

1400

1600

1800

2000

Figura 5.34.: Seal en con sus tres componentes en el dominio temporal

84

5.6 Eliminacin de ruido transientes


12

1,31 Hz

0,75 Hz

eoc709
H/V1
H/V2
H/V3
H/V4

10

1,41 Hz

0,72 Hz
0

0.5

1.5

2.5

3.5

4.5

Figura 5.35.: Anlisis espectral sin anti triggering

Las frecuencias encontradas las podemos resumir en la Tabla 5.13, donde en la primera columna encontramos el nombre del conjunto de parmetros de procesado
((H/V)1 , (H/V)2 , (H/V)3 y (H/V)4 ), y en la segunda columna tenemos la frecuencia
estimada para cada conjunto de parmetros, expresada en Hz.
Conjunto de Parmetros Frecuencia Estimada (Hz)
(H/V)1
1,389
H
( /V)2
1,401
(H/V)3
1,408
H
( /V)4
1,386
Cuadro 5.13.: Frecuencias estimadas para la primera prueba, en donde no se utiliz anti triggering

Ahora realizaremos la segunda prueba, la cual consiste en estimar nuevamente la


frecuencia natural, ocupando los mismos (H/V)1 , (H/V)2 , (H/V)3 y (H/V)4 pero esta
vez eliminaremos el ruido transiente presente en la seal, los resultados los podemos observar en la Figura 5.36. Donde claramente encontramos que la frecuencia
estimada ya no corresponde a las encontradas anteriormente. Esto a su vez depende
del conjunto de parmetros que sean escogidos, as como tambin que parte de los
ruidos son extrados en la seal.

85

5.6 Eliminacin de ruido transientes


12
eoc709
H/V1
H/V2
H/V3
H/V4

1,31 Hz

0,74 Hz
10

0,72 Hz

1,40 Hz

0.5

1.5

2.5

3.5

4.5

Figura 5.36.: Estimaciones Frecuenciales para la segunda prueba con anti triggering
En la Tabla 5.14 se pueden ver los resultados de las frecuencias encontradas en
(H/V)1 , (H/V)2 , (H/V)3 y (H/V)4
Conjunto de Parmetros Frecuencia Estimada (Hz)
(H/V)1
1,397
H
( /V)2
0,726
H
( /V)3
1,401
(H/V)4
1,397
Cuadro 5.14.: Frecuencias estimadas para la segunda prueba, en donde se utiliz
anti triggering

En la Tabla 5.15 podemos observar la diferencia frecuencial entre las dos pruebas
en valor absoluto.
Conjunto de Parmetros Diferencia de Frecuencia en valor absoluto (Hz)
(H/V)1
0,008
H
( /V)2
0,675
(H/V)3
0,007
H
( /V)4
0,011
Cuadro 5.15.: Corrimiento frecuencial en valor absoluto entre la prueba uno y la
prueba dos

86

5.6 Eliminacin de ruido transientes


Como se puede ver la frecuencia natural con mayor cambio es cuando se procesan
los datos con la configuracin de parmetros que esta definida en (H/V)2 , en donde
el corrimiento frecuencial es de aproximadamente de 0,675 Hz, y el conjunto de
parmetros donde se da esta situacin es:
Ventana: 40 segundos
Solape: 10 %
Filtro preprocesamiento: Butterworth orden 2, 0.1Hz-20Hz
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
Como se observa en la Tabla 5.15 la eliminacin del ruido transiente altero la estimacin frecuencial, pero esto a su vez se atena o se acrecienta de acuerdo al
conjunto de parmetros de procesado que se utilicen. La diferencia en la estimacin
del espectro para (H/V)2 con las dos pruebas se pueden ver en la Figura 5.37 donde
la linea roja es el espectro obtenido de la seal sin anti triggering, y la linea negra
representa a la prueba (H/V)2 con anti triggering.
8
H/V2 Sin anti triggering
H/V2 Con anti triggering
7

1,04 Hz

0,72 Hz
3

0.5

1.5

2.5

Figura 5.37.: Comparacin de los (H/V)2 en las dos pruebas realizadas

87

6. Captulo VI Conclusiones
El estudio realizado tuvo por objetivo analizar la robustez de la estimacin de los
espectros utilizados en la razn espectral H/V de Nakamura. Para ello se tomaron
alrededor de 100 registros de microvibraciones obtenidas por el equipo de Obras Civiles, correspondientes a localidades del sur de Chile, con ellos se procedi a realizar
el respectivo procesado de seales, utilizando tanto tcnicas tradicionales, as como
tcnicas ms avanzadas de procesamiento de seales. Se utilizo el software de procesamiento de seales sismolgicas geopsy para obtener los espectros de Nakamura,
y de esta forma compararlos con los espectros obtenidos por el equipo de Obras
Civiles. Adems de esto, se procesaron las seales en Matlab, utilizando tcnicas
modernas para obtener los espectros y realizar comparaciones entre los diferentes
mtodos. Es as que hemos llegado a las siguientes conclusiones:
Al efectuar una comparacin entre las frecuencias naturales obtenidas por eoc
y las determinadas con el Software Geopsy y realizar un analisis de dispersin
con Histogramas, se llego a la conclusin de que de 102 seales, el 90 % posee
un dispersin menor a 0,5 Hz, por lo que se puede decir que la estimacin de
los espectros de Nakamura es robusta, y en general existe una clara convergencia. Pero en ciertos casos existen excepciones, que en valor porcentual es
alrededor de un 7,51 %, estas excepciones pueden deberse a diversas causas:
que se tomaron parametros demasiado sencibles a las dinamicas de los suelos,
la forma en que se tomaron los registros, si los registros se tomaron de dia o
de noche, el tipo de dinamica del suelo, fuentes de ruido externas, etc.
Al realizar el estudio de sensibilidad de parmetros, la cual se realizo con el
software de estimacin de espectros de Nakamura, Geopsy, y el cual consista
en realizar 4 experimentos variando los parmetros de ventaneo y filtro de
procesamiento, se observaron en algunos casos variaciones dentro de un mismo
registro entre los 4 experimentos de Geopsy, por lo que existe en cierto grado
un problema de sensibilidad de parmetros.
En los registros en donde se da el caso de peaks claros, las frecuencias naturales seleccionada por el equipo de Obras Civiles y la obtenida al utilizar
Nakamura con modelo AR, coinciden. Existe una tendencia a que los registros
con frecuencias naturales mayores a 1,2 Hz coincidan en mayor cantidad de
veces con lo obtenido por AR, la mayor cantidad de coincidencias se da en las
seales muestreadas a 50 Hz, esto puede deberse a que el algoritmo no posee
un comportamiento optimo con las seales muestreadas a 100 Hz.

88

Captulo VI Conclusiones
En las seales de 100 Hz existen problemas para ver las frecuencias que se
encuentran bajo 1 Hz. Intentando mejorar esto, se utilizo la seal filtrada
de Geopsy, se aumento el orden del modelo AR, y se disminuy el solape
del algoritmo de Matlab en 100 muestras. En estas seales la mayora de los
resultados obtenidos con AR, muestran frecuencias ms altas que las obtenidas
por el equipo de Obras Civiles.
Al observar la comparacin de los registros de 50 y 100 Hz de los registros de
Valdivia, se muestran grandes diferencias entre las frecuencias obtenidas por
eoc y las obtenidas con el mtodo AR, esto puede ser posible debido a que
los registros de Valdivia, poseen una baja frecuencia natural, por lo cual seria
un problema con las bajas frecuencias y no con el muestreo de las seales.
Aqu se pueden dar dos casos: el primero seria que el mtodo utilizado por
eoc funciona mal a bajas frecuencias, o segundo, el mtodo AR funciona mal
en bajas frecuencias, este es el fundamento que da pie a la segunda parte de
esta tesis, ya que se deber revisar los resultados con mas detalles y ver si
mejorando el mtodo AR, se presentan o no estas dificultades.
Existen registros en donde algunos peak elegidos como frecuencias naturales,
corresponden a ruidos, los cuales ocultan la real frecuencia, esto es posible
observar con mayor claridad al aplicar AR y observar mesh de las seales.
Esto nos indica que el utilizar tcnicas ms avanzadas contribuye a obtener
mejores anlisis de cada registro y de esta forma mejorar los resultados.
Uno de los ltimos anlisis realizados y que an no se ha profundizado, es
observar la variacin, al eliminar los transientes de las seales, los cuales producen error en la obtencin de la frecuencia natural. Para este anlisis, se
tomaron algunos registros con dos o ms peaks cercanos entre si y se seleccionaron manualmente los trazos mas estacionarios. En los casos observados,
se aprecia una clara mejora de los espectros, haciendo que los peaks sean ms
claros y varen sus frecuencias. La desventaja de esto es que toma mucho tiempo seleccionar cada trazo de la seal, lo cual no es beneficioso a la hora de
realizar un estudio de Microzonificacin.
Finalmente podemos concluir que al utilizar tcnicas tradicionales se obtienen estimaciones robustas de los espectros, pero se presentan ciertos casos
en los cuales existen algunas particularidades, como ruidos impulsionales que
modifican la frecuencia natural escogida. Al utilizar tcnicas modernas de tratamiento de seales, en este caso el modelo AR, para la estimacin de los
espectros de Nakamura, se obtienen resultados en los que es ms fcil detectar
dichas particularidades, ya que esta tcnica es ms inmune a los ruidos que
ocultan las verdaderas frecuencias naturales de cada registro. Adems al aplicar el algoritmo implementado para este estudio, se obtienen diversas formas
de visualizar los resultados, y as realizar un mejor anlisis y tomar una mejor
decisin. Cabe sealar que an queda mucho trabajo por realizar para tener

89

Captulo VI Conclusiones
un mtodo con la robustez deseada. Este estudio aporta a dicho trabajo al
demostrar que con la aplicacin de mtodos de tiempo-frecuencia, se obtienen
resultados interesantes para un posterior anlisis ms profundo.

90

Bibliografa
[Alf]

J. J. Egozcue y A. Ugalde Alfaro. Determinacin de caractersticas


dinmicas del suelo a partir de microtemblores. Universidad Politcnica
Catalua. Departamento de Ingeniera del Terreno y Cartogrfica.

[ALI]

ALI. Universidad de Alicante, Riesgo Ssmico: Estimacin e incertidumbre de daos ssmicos y prdidas para entornos urbanos,
http://riesgosismico.es/contenido.php@riesgo.html.

[Ani]

Ojeda C. Anibal. Evaluacin de efectos de sitio para la microzonificacin smica de Pereira.

[BS99]

T. Bartosch and D. Seidl. Spectrogram analysis of selected tremor signals using short-time Fourier transform and continuous wavelet transform. Annali di Geofisica, Vol. 42, No. 3, pp. 497-506, 1999.

[CB94]

Lachet C. and P. Y. Bard. Numerical and theoretical investigations


on the possibilities and limitations of the Nakamura is technique, Jour.
Phys. Earth, 42, 377-397., 1994.

[CJ07]

Lara Castro and Omar Jos. Nuevas metodologas no invasivas de diagnosis de defectos incipientes en rodamientos de bola, Cap. 2. Universidad Carlos III de Madrid, Depto de Ingeniera Mecnica, 2007.

[Con]

Sheila Alodia Yauri Condo. Microtremores. Escuela Profesional de Ingeniera Geofsica Universidad Nacional de San Agustn de Arequipa.

[dAM]

Rafael de Arce and Ramn Maha. MODELOS ARIMA: Definiciones


bsicas, Dpto. Economa Aplicada U.D.I. Econometra e Informtica.

[DGA88]

Benedetti D., Benzoni G., and Parisi M. A. Seismic vulnerability and


risk evaluation for old urban nuclei, Earthquake Engineering and Structural Dynamics, 16, 183-201., 1988.

[eA07a]

Emmanuel Chaljub et Al. Spectral Element Analysis in Seismology.


Advances in Geophysics, Vol. 48, pp. 365-419, 2007.

91

Bibliografa
[eA07b]

P.R. Scalassara et Al. Autoregressive decomposition and pole tracking


applied to vocal fold nodule signals. Pattern Recognition Letters, Vol.
28, pp. 1360-1367, 2007.

[eA08]

S. Bonnefoy Claudet et Al. Site effect evaluation in the basin of Santiago


de Chile using ambient noise measurements Geophys. J. Int., 2008.

[Efr03]

Grajales Vargas Efran. Tesis Maestra en Ciencias: Comportamiento


de los Edificios con Dispositivos Disipadores de Energa, Universidad
Autonoma de Nuevo Len, Mexico, 2003.

[Eli08]

Lara Aguayo Adriana Elizabeth. Mtodos lineales y no lineales del anlisis espectral de los procesos y campos aleatorios hidrometeorolgicos.
Departamento de Fsica, CUCEI, Universidad de Guadalajara, 2008.

[EPo]

Anlisis de seales, Capitulo 2, Biblioteca virtual, Escuela Politcnica


Nacional, Quito, Ecuador.

[F.10]

Guerrero Martinez Juan F. Ingeniera Biomdica, Curso 2010-2011.


Escola Tcnica Superior D Enginyeria, Departament d Enginyeria Electrnica, 2010.

[FAM+ 10]

Leyton F., Seplveda S. A., Astroza M., Rebolledo S., Gonzlez L.,
Ruiz S., Foncea C., Herrera M., and Lavado J. ZONIFICACION SSMICA DE LA CUENCA DE SANTIAGO, CHILE, Congreso chileno
de sismologa e ingeniera antissmica. X jornadas. Valdivia- Santiago,
Chile, 2010.

[FKG03]

D. Fah, F. Kind, and D. Giardini. Inversion of local S-wave velocity


structures from average H/V ratios, and their use for the estimation
of site-effects, Journal of Seismology, Volumen 7, Issue 4, pp 449-467.,
2003.

[Gab46]

D. Gabor. Theory of Communication. Journal of Institute for Electrical


Engineering 93, part III (26)., 1946.

[Geoa]

Geopsy. H/V Toolbox, http://www.geopsy.org/wiki/index.php/Geopsy:_H/V_and_Spectum_Toolboxes:_Processing_Tab.

[Geob]

Geopsy.
Filtro Konno y Ohmachi, http://www.geopsy.org/wiki
/index.php/-Smoothing_details.

[Geoc]

Geosig. http://www.geosig.com/.

[HA95]

Panbong Ha and Souguil Ann. Robust time-varying parametric modelling of voiced speech. Signal Processing, Vol. 42, pp. 311-317, 1995.

92

Bibliografa
[Her08]

M. Herak. Model HVSR A Matlabs tool to model horizontal-to-vertical


spectral ratio of ambient noise. Computers and Geosciences, Vol. 34,
pp. 1514-1526, 2008.

[HO10]

Jens Havskov and Lars Ottemoller. Routine Data Processing in Earthquake Seismology. Springer, 2010.

[Ign09]

Begoa Moya Ignacio. Anlisis Tiempo-Frecuencia de la Seal de Vibracin de un Cambiador de Tomas en Carga, 2009.

[Ill03]

A. Illanes. Estimation parametrique robuste pour l analyse de la variabilit cardiaque. Tesis de Master. Laboratoire I3S, CNRS. France,
2003.

[JCD]

Ritta Ral J., Massa Julio C., and Chiappero Guillermo D. Identificacin del perodo fundamental de sitio a partir de registros de vibracin
ambiental. Universidad Nacional de Crdoba, Departamento de estructuras.

[JE00]

Ibez J.M. and Carmona E. Sismicidad Volcnica. Instituto Andaluz


de Geofisica. Universidad de Granada, 2000.

[JEA12]

Ritta Ral J., Surez Luis E., and Pando Miguel A. Determinacin del
Perodo Fundamental del Suelo usando vibracin ambiental y el cociente
espectral H/V, Mecnica Computacional Vol XXXI, pgs. 1399-1419
(artculo completo). Argentina, 13-16 Noviembre, 2012.

[JNS03]

Philippe Jousset, Jurgen Neuberg, and Susan Sturton. Modelling the


time-dependent frequency content of low frequency volcanic earthquakes.
Journal of Volcanology and Geothermal Research, Vol. 128, pp. 201223, 2003.

[JP07]

Cceres Juan-Pablo. Anlisis Espectral 1: Transformada Corta de Fourier y Ventanas, 2007.

[JWH02]

C. Ji, D. J. Wald, and D. V. Helmberger. Source Description of the


1999 Hector Mine, California, Earthquake, Part I: Wavelet Domain
Inversion Theory and Resolution Analysis. Bulletin of the Seismological
Society of America, Vol. 92, No. 4, pp 1192-1207, 2002.

[Kay99]

S.M. Kay. Modern Spectral Estimation: Theory and Application. Prentice Hall, 1999.

[Kit10]

G. Kitagawa. Introduction to time series modeling. CRC Press, 2010.

[KO98]

K. Konno and T. Ohmachi. Ground-Motion Characteristics Estimated from Spectral Ratio between Horizontal and Vertical Components of
Microtremor. Bull. Seism. Soc. Am., 88, 228-241., 1998.

93

Bibliografa
[LAS04]

Arana Lpez-Abad and M. Soledad. Interfaz grfica para el clculo


de la dinmica longitudinal y prestaciones de un vehculo mediante la
herramienta guide de Matlab, Proyecto Fin de Carrera de la Escuela
Politcnica Superior. Universidad Carlos III de Madrid, Madrid, 2004.

[LCG94]

J. Lermo and F. J. Chavez-Garcia. Are microtremors useful in site


response evaluation?, Bull. Seism. Soc. Am., 83,1350-1364., 1994.

[Les09]

P. Lesage. Interactive Matlab software for the analysis of seismic volcanic signals. Computers and Geosciences, Vol. 35, pp. 2137-2144, 2009.

[LESA]

Bermudez Maria L, Franco Luis E., Martnez Samuel, and Ojeda


Anibal. Uso de las microtrepidaciones para la evaluacin de la respuesta dinmica de los suelos.

[LGM02]

Philippe Lesage, Francois Glangeaud, and Jerome Mars. Applications


of autoregressive models and time-frequency analysis to the study of volcanic tremor and long-period events. Journal of Volcanology and Geothermal Research, Vol. 114, pp. 391-417, 2002.

[Lop]

Clara Eugenia Giraldo Lopera. Licuacin del suelo. Mecinica de Fluidos y Recursos Hidriulicos, Escuela de Ingenieria de Antioquia.

[lW]

Enciclopedia libre Wikipedia. Historial de versiones de MatLab.


http://es.wikipedia.org/wiki/MATLAB.

[Mal08]

S. Mallat. A Wavelet Tour of Signal Processing. Academic Press, 2008.

[MGC99]

Jean Luc Mari, Francois Glangeaud, and Francoise Coppens. Signal


Processing for Geologists and Geophysicists. Institut Francais Du Petrole Publications, Editions Technip Hardcover, 1999.

[MGSJV10] M. Martnez, L. Gmez, A. J. Serrano, and J. Gmez J. Vila. Filtros Digitales. Escola Tcnica Superior d Enginyeria, Departament d
Enginyeria Electrnica, Universitat d Valencia, 2010.
[MH]

U. Mena Hernandez. Evaluacin del riesgo ssmico en zonas urbanas.


Tesis Universidad Politcnica Catalua, Barcelona.

[Mig10]

Zamorano Castao Miguel. ANALISIS DE SEALES MEDIANTE


STFT Y WAVELET. APLICACIN A DEFECTOLOGA EN RODAMIENTOS. Universidad Carlos III de Madrid, Legans, 2010.

[MKBB05] O. Meste, B. Khaddoumi, G. Blain, and S. Bermon. Time-Varying


Analysis Methods and Models for the Respiratory and Cardiac System
Coupling in Graded Exercise. IEEE Transactions on Biomedical Engineering, Vol. 52, N0. 11, 2005.

94

Bibliografa
[Mor09]

Asuncin Moreno. Seales y Sistemas I, Filtros. Departamento de Teora de la Seal y Comunicaciones Universitat Politcnica de Catalunya
(UPC), 2009.

[NA09]

A. Nait-Ali. Advanced Biosignal Processing. Springer, 2009.

[Nak89]

Y. S Nakamura. Method for dynamic Characteristics estimation of


Subsurface using microtremor on the Ground Surface. Geotechnical eng.
And Disaster Prevention Laboratory, QR of RTRI, Vol. 30, No. 1, 1989.

[Nak08]

Y. Nakamura. On the H/V spectrum. The 14th World Conference on


Earthquake Engineering October 12-17. Beijing, China, 2008.

[Oba]

Tupak Obando. Microzonificacin ssmica. Doctorado, y Master en


Geologa, y Gestin Ambiental por la Universidad Internacional de Andalucia UNIA.

[oSR89]

EERI Committee on Seismic Risk. The basics of seismic risk analysis,


Earthquake Spectra, 675-702., 1989.

[Pil07]

Cuadra Monreal Pilar. Tesis: Aplicacin de Tecncas de Vibraciones


Ambientales: Anlisis de Microtremores y vibraciones naturales, para la
caracterizacin de sitio. Universidad Simn Bolvar, Sartenejas, Abril,
2007.

[Pro01]

John G. Proakis. Digital communications. McGraw-Hill, 2001.

[Pro04]

Sesame European Research Project. Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation, 2004.

[R.03]

Bonnett R. Tesis Doctoral: Vulnerabilidad y riesgo Ssmico de edificios.


Aplicacin a entornos urbanos en zonas de amenaza alta y moderada.
ETS de Ingeniera de Caminos, Canales y Puertos de Barcelona, Universidad Politcnica Catalua. Barcelona, 2003.

[Raf07]

Caballero Guerrero lvaro Rafael. Tesis de Maestra: Determinacin


de la Vulnerabilidad ssmica por medio del Mtodo del ndice de vulnerabilidad en las estructuras ubicadas en el centro Histrico de la ciudad
de Sincelejo. Facultad de Ingeniera, Fundacin Universidad del Norte,
Sincelejo, 15 de julio, 2007.

[RL10]

S. Ruiz and F. Leyton. Frecuencia fundamental del suelo observado con


acelerogramas del terremoto de Valparaso 1985 y microvibraciones. X
Congreso de Sismologa e Ingeniera Antissmica. Santiago, Chile, 2010.

95

Bibliografa
[Rod05]

M. Rodrguez. Caracterizacin de la Respuesta Ssmica de los Suelos. Aplicacin a la ciudad de Barcelona. Tesis Ingeniera de Caminos,
Canales y Puertos. Barcelona. Universidad Politcnica Catalua, 2005.

[Rsi]

Riesgo ssmico. http://riesgosismico.es/contenido.php@riesgo.html.

[Rub]

Sonia lvarez Rubio. El Efecto Local sobre el Movimiento Ssmico del Suelo: Fenomenologa y Resultados Recientes, Capitulo 2,
E.U.I.T.Industrial. Universidad Politcnica de Madrid.

[SAM+ ]

Figueras Sara, Macau Albert, Belvaux Myriam, Peix Miquel, Benjumea Beatriz, Gabas Anna, and Susagna Teresa Guola Xavier.
Caracterizacin de efectos ssmicos
locales en la ciudad de Lorca. http://hal-brgm.archivesouvertes.fr/docs/00/70/42/71/PDF/Fisica_de_la_tierra_suelos_SF.pdf.

[Sch94]

Frank Scherbaum. Basic Concepts in Digital Signal Processing for Seismologists. Lecture Notes in Earth Sciences. Springer, 1994.

[Seb10]

Araujo Sebastan. Un Modelo de fuente ssmica para el tremor volcnico basado en la turbulencia elstica del magma andestico y el lodo
volcanico. Universidad Politecnica Saleciana, Ecuador, 2010.

[Sem04]

J.L. Semlow. Biosignal and Biomedical Image Processing. Marcel Dekker Inc., 2004.

[SESa]

SESAME. Pgina web de GEOPSY. http://www.geopsy.org.

[SESb]

SESAME.
Descripcin del programa computacional GEOPSY.
http://www.geopsy.org/wiki/index.php/Main_Page.

[SL05]

L. Sommo and P. Laguna. Bioelectrical signal processing in cardiac and


neurological applications. Elsevier Academic Press, 2005.

[SLM10]

V. Soto, F. Leyton, and J. Maripangui. Consideraciones sobre la tcnica


de H/V en microvibraciones y resultados preliminares de mediciones
en Curic y Talca, Regin del Maule. Congreso Chileno de sismologa
e ingeniera antissmica. Valdivia-Santiago, Chile, 2010.

[SM05]

P. Stoica and R. Moses. Spectral Analysis of Signals . Prentice Hall,


2005.

[Smi97]

S.W. Smith. The Scientist & Engineers Guide to Digital Signal Processing. California Technical Pub, 1997.

[Som00]

P. Sommerville. Seismic Hazard Evaluation. State of the art. Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, Paper No.2833., 2000.

96

Bibliografa
[TK03]

Tetsuo Takanami and Genshiro Kitagawa. Methods and Applications


of Signal Processing in Seismic Network Operations . Lecture Notes in
Earth Sciences . Springer, 2003.

[Tra10]

Martn H. Trauth. Matlab Recipes for Earth Sciences. Springer, 2010.

[Tro]

Tromino. http://www.tromino.eu/.

[UGU]

Gua bsica para la elaboracin de atlas estatales y municipales de peligro y riesgo. Fenmenos geolgicos, Captulo III. Ciencias de la salud
ambiental. Universidad de Guadalajara.

[Val12]

Galo Valdebenito. Base de datos de estudios de microzonificacin. Instituto de Obras Civiles, Universidad Austral de Chile, 2012.

[VPCB05]

R. Verdugo, C. Pastn, J. Campos, and F. Bonilla. SO DE LA RAZON ESPECTRAL H/V EN SUPERFICIE PARA LA CARACTERIZACIN SSMICA DE SANTIAGO", Congreso Chileno de Sismologa
e Ingeniera Antissmica IX Jornadas. Concepcin, Chile, 2005.

[yEC08]

Molina Sergio y Elas Cesar. Microzonacin Ssmica del rea Metropolitana de San Salvador (AMSS). Noruega, 2008.

[Yep96]

F. Yepez. Metodologa para la evaluacin de la vulnerabilidad y riesgo


ssmico de estructuras aplicando tcnicas de simulacin. Tesis Doctoral,
Escuela Tcnica Superior de Ingenieros de Caminos, Canales y Puertos. Universidad Politcnica Catalua, Barcelona, 1996.

[ZTC03]

H. Zhang, C. Thurber, and C.Rowe. Automatic P-Wave Arrival


Detection and Picking with Multiscale Wavelet Analysis for SingleComponent Recordings. Bulletin of the Seismological Society of America, Vol. 93, No. 5, pp. 1904-1912, 2003.

97

A. Planes de trabajo y
sectorizaciones de las zonas de
estudios
A.1.

Plan de trabajo y sectorizacin de la ciudad de


Valdivia

A continuacin se explicarn los detalles tcnicos y protocolos de trabajo que se


usaron durante la confeccin de la base de datos, esta informacin fue proporcionada
por personas del Instituto de Obras Civiles de la Universidad Austral de Chile.
El rea de estudio abarc casi en su totalidad el sector urbano de la ciudad de Valdivia. En trminos de superficie, la zona de estudio contempla un total aproximado
de 25.66 km2. A su vez, se ha subdividido la zona de estudio en 7 sectores, esto con
el fin de establecer una mejor planificacin en la campaa de campo y una nomenclatura que facilite la ubicacin de los puntos de medicin en la zona. Teniendo esto
presente, la Figura A.1 muestra la sectorizacin de la zona de estudio.

Figura A.1.: Sectorizacin de la zona de estudio

98

A.1 Plan de trabajo y sectorizacin de la ciudad de Valdivia


La cantidad de puntos de medicin de microvibraciones ambientales en cada sector,
la superficie de los mismos, y la densidad de puntos por km2 se muestra en la
siguiente tabla:
Sector Superficie (km2 ) N Puntos N Puntos/km2
1
1.06
93
88
2
2.9
28
10
3
3.1
22
7
4
6.4
33
5
5
6.6
34
5
6
2.4
18
8
7
3.2
30
9
Total
25.66
258
10
Cuadro A.1.: Nmero de puntos de medicin por unidad de superficie en cada
sector

Cabe destacar que el espaciamiento entre la gran mayora de los puntos de medicin
resulta acorde a las recomendaciones para estudios de microzonificacin propuestos
por SESAME [Pro04], en donde se sugiere adoptar un espaciamiento inicial de 500
m entre dos puntos adyacentes de registro. No obstante, en el caso del Sector 1,
correspondiente a la zona cntrica de la ciudad, el espaciamiento entre puntos fue
mucho menor, por lo cual es el sector de estudio que posee la mayor densidad de
puntos por km2 (88 Puntos/km2). Se adopt este criterio en el plan de trabajo
en vista de que evidentemente en este sector se concentra la mayor cantidad de
actividades de la ciudad, lo cual enfocado desde el punto de vista del riesgo ssmico,
determinar una alta exposicin en la zona acorde a la mayor densidad poblacional
que esta posee.
Con respecto a la nomenclatura de los archivos generados, cada punto de medicin
es referenciado por un nmero de tres dgitos, donde el primer nmero representa
el sector y los dos nmeros siguientes corresponden a la numeracin definida para
dicho punto dentro del sector al cual pertenece. (Ej: Punto N 718 representa el
punto nmero 18 del sector 7).

A.1.1.

Descripcin de equipos y parmetros de adquisicin

Para la adquisicin de datos de microvibraciones ambientales se requieren equipos


triaxiales que presenten curvas de respuesta plana en velocidad (gefonos). En este
estudio se emplearon equipos GeoSIG [Geoc] y equipos Tromino ENGY [Tro].

99

A.1 Plan de trabajo y sectorizacin de la ciudad de Valdivia


Los equipos GeoSIG GVB-316 ( Figura A.2 ), de origen suizo, registran ondas de
frecuencia comprendidas entre 4.5 HZ y 50 Hz. Poseen una respuesta confiable de
hasta 1 Hz, y se han utilizado exitosamente en otros estudios de microzonificacin
[RL10]. Con estos equipos se realizaron mediciones de microvibraciones de 15 a 20
minutos de duracin, cada una a una tasa de muestreo de 50 Hz, con el fin de obtener
longitudes de registro adecuadas para caracterizar frecuencialmente el suelo.

Figura A.2.: Equipo sismolgico de medicin Geosig GBV 316


Los equipos Micromed Tromino ENGY (Figura A.3), de origen italiano, poseen sensores capaces de registrar ondas de frecuencia comprendidas entre 0.1 Hz y 256 Hz.
Con estos equipos se realizaron mediciones de microvibraciones de 16 minutos cada
una, a una tasa de muestreo de 128 Hz.

Figura A.3.: Equipo sismolgico de medicin Tromino ENGY


Cabe destacar que para la toma de datos en terreno, los equipos fueron adecuadamente orientados segn norte magntico, georreferenciados y configurados para la
captura de datos usando los transductores de velocidad en alta ganancia, como es
caracterstico para mediciones de microvibraciones ambientales.

100

A.1 Plan de trabajo y sectorizacin de la ciudad de Valdivia

A.1.2.

Descripcin de la campaa de campo

El procedimiento para la adquisicin de datos en terreno se puede resumir como


sigue:
Una vez ubicado el lugar donde se desea efectuar la medicin, se verifican las condiciones ptimas para la correcta instalacin del equipo. Se debe tener en cuenta
minimizar al mximo los factores que puedan generar error en las mediciones, como
la presencia de cables de alta tensin, excesiva cercana con trfico vehicular y/o
circulacin de peatones, evidencia de tuberas subterrneas, condiciones climticas,
entre otros.
Verificado el punto anterior, se busca una superficie plana de contacto entre el suelo y
el equipo, la cual debe ser preferentemente contra terreno natural, en un suelo donde
se aprecie cierto grado de compactacin. Tambin se pueden efectuar las mediciones
sobre losas de hormign o asfalto. Nunca se debe instalar el instrumento sobre pasto,
suelos blandos, sueltos o con un elevado grado de humedad.
Se orienta el equipo en direccin al norte magntico, para luego proceder a su fijacin
y nivelacin en el terreno, asegurndose que ninguno de los extremos del aparato
quede separado del suelo.
Utilizando un aparato GPS, se registran las coordenadas geogrficas del punto. En el
caso de los equipos GeoSIG, los parmetros de adquisicin de datos son configurados
en un ordenador porttil conectado al equipo. Luego de definir la frecuencia de
muestreo de la grabacin, el inicio de la adquisicin de datos se hace configurando
la opcin Time Trigger, ingresando la fecha, la hora y el tiempo de duracin de la
grabacin.
Los equipos Tromino son configurados manualmente de forma directa en el equipo,
en donde se fija el tiempo de grabacin y se activan los tres canales velocimtricos
para la adquisicin de datos en alta ganancia. El inicio de adquisicin de datos tiene
un desfase de 15 segundos luego de la configuracin manual del equipo.
Una vez comenzado el proceso de grabacin, se debe evitar al mximo el trnsito
de personas, vehculos y animales por las inmediaciones de los equipos, con el fin de
garantizar el menor ruido local posible en la adquisicin de datos.
Se registra en la planilla de campo informacin referente a la fecha y hora de adquisicin, georreferenciacin, superficie de contacto entre el suelo y el equipo, fuentes de
ruido y algn otro tipo de informacin pertinente. Finalizado el tiempo de grabacin,
se apaga el equipo y se realiza el traslado al punto siguiente.

101

A.2 Protocolo de procesamiento de seales

A.2.

Protocolo de procesamiento de seales

Para dar un orden y un procedimiento claro al trabajo que se realiz sobre la base
de datos, se cre un protocolo de trabajo con el objetivo de dejar numeradas y
clasificadas las pruebas realizadas. Este protocolo especfica paso a paso lo que se
debe ir realizando con cada uno de los registros que dispone la base de datos, es
as como se comienza por cargar los datos que fueron procesados por el Instituto
de Obras Civiles, ya que esto sirve para ms adelante realizar comparaciones con
nuestro estudio, la informacin que entrega el equipo del Instituto de Obras Civiles
est bajo la forma de un archivo Excel del cual se extraen los datos y se pasan
a un archivo de texto. Luego se traspasan los datos de microvibraciones desde los
archivos que entrega como salida el gefono a un formato de texto ms estndar,
y cada seal es identificada con el nombre de la componente de la cual proviene
(Norte-Sur, Este-Oeste, Vertical). Luego estos datos son utilizados para crear un
archivo con el nombre de la ciudad y acompaado del nmero de registro.
Una vez realizado el trabajo de conversin de archivos se pasa a un procesado con
Geopsy, se debe cargar el archivo que contiene las tres componentes de la seal
microvibratoria, el cual por nombre lleva el sector y el nmero de registro donde
fue tomada. A continuacin se debe indicar la frecuencia de muestreo y tambin se
introduce el orden en el cual estn las componentes de la seal de microvibraciones.
Luego se trabaja cada seal microvibratoria con cuatro configuraciones de parmetros de procesado distintos, adems se obtiene el espectro de cada una de las tres
componentes con un conjunto de datos definidos tambin en el protocolo.
En una tercera etapa se defini el trabajo que se realiza en MatLab ocupando todos
los archivos anteriores, aqu tambin se calcula la transformada de Fourier de cada
componente y se compara con la obtenida con Geopsy, se estiman espectros con
otras tcnicas de procesado de seales, se guardan imgenes que son obtenidas con
el comando imagesc de Matlab, todo con un estricto orden, de acuerdo al nmero
de registro y al tipo de procesado que se realiz en ella.
Finalmente se estandarizaron las conclusiones que se extraeran de cada seal procesada, colocando nfasis en ver posibles diferencias entre distintas tcnicas de procesado, todo el detalle del protocolo de trabajo se muestra a continuacin, donde se
explica paso a paso cuales fueron los nombres y pruebas que se llevaron a cabo

102

A.3 Sectorizacin de la localidad de Capitn Pastene

A.3.

Sectorizacin de la localidad de Capitn


Pastene

Figura A.4.: Plan de trabajo de Capitn Pastene

A.4.

Sectorizacin de la localidad de Mininco

Figura A.5.: Plan de trabajo de Mininco

103

A.5 Sectorizacin de la localidad de Lumaco

A.5.

Sectorizacin de la localidad de Lumaco

Figura A.6.: Plan de trabajo de Lumaco

A.6.

Sectorizacin de la localidad de Purn

Figura A.7.: Plan de trabajo Purn

104

A.7 Sectorizacin de la localidad de Pitrufqun

A.7.

Sectorizacin de la localidad de Pitrufqun

Figura A.8.: Plan de trabajo de Pitrufqun

105

B. Protocolo de Procesamiento
Lo primero que se debe realizar al momento de comenzar con el trabajo de procesamiento es extraer todos los datos desde los archivos originales, esta labor se realiza
porque se debe trabajar con formatos de texto mas estndar, para un posterior
trabajo en MatLab.
Archivo Excel: Cada punto geogrfico registrado en la base de datos posee un
archivo Excel que contiene el estudio completo realizado por el equipo de obras
civiles. El SRSS corresponde al grfico de H/V con respecto a la frecuencia
en Hz. Si se pincha esta curva se destacan automticamente las columnas que
forman el grfico de H/V. En la mayora de los casos la columna E corresponde
a la frecuencia y la columna F al H/V.
TCO_20110420_144100_TCO___SE__: archivo que contiene la componente Este-Oeste. La primera columna corresponde al tiempo y la segunda
columna a la amplitud. Se debe tener atencin al hecho de que no todos los
registros estn muestreados a la misma frecuencia de muestreo.
TCO_20110420_144100_TCO___SN__: archivo que contiene la componente Norte-Sur y que posee similares consideraciones que el archivo que contiene la componente Este-Oeste.
TCO_20110420_144100_TCO___SZ__: archivo que contiene la componente Vertical y que posee similares consideraciones que el archivo que contiene
la componente Este-Oeste.
Con el fin de hacer la carga ms fcil en Matlab es mejor pasar cada uno de
estos tres archivos a sei.txt, sni.txt, szi.txt. Donde el valor i corresponde al
nmero de registro del lugar muestreado.
Seales a tener en cuenta para ser cargadas con el fin de realizar el anlisis por
registro (i corresponder al nmero de registro):
eoci.txt: Estudio H/V realizado por el equipo de obras civiles (el cual fue
extrado desde el archivo Excel en la etapa anterior). Este corresponde a un
archivo de dos columnas donde la primera representa el vector frecuencia y la
segunda representa el vector amplitudes.

106

Protocolo de Procesamiento
Valdiviai.txt: Archivo donde cada columna contiene los registros de cada una
de las componentes. La primera columna corresponde al registro E-O, la segunda al registro N-S y la tercera al registro V.
Guardar un archivo llamado filtrogeop, en el caso de las seales cuya frecuencia
de muestreo es de 100 Hz, este filtro es un Butterworth orden 3, entre 0.05-10
Hz antes de realizar las pruebas geophv1_i.hv, geophv2_i.hv, geophv3_i.hv,
geophv4_i.hv.
geophv1_i.hv: Primer estudio H/V1 realizado en Geopsy en base a los siguientes parmetros:
Ventana: 20 segundos
Solape: 0 %
Filtro preprocesamiento: Butterworth orden 3, 0.1Hz-10Hz
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
geophv2_i.hv: Segundo estudio H/V2 realizado en Geopsy en base a los siguientes parmetros:
Ventana: 40 segundos
Solape: 10 %
Filtro preprocesamiento: Butterworth orden 2, 0.1Hz-20Hz
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
geophv3_i.hv: Tercer estudio H/V3 realizado en Geopsy en base a los siguientes parmetros:
Ventana: 60 segundos
Solape: 50 %
Filtro preprocesamiento: Butterworth orden 3, 0.1-15Hz

107

Protocolo de Procesamiento
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
geophv4_i.hv: Cuarto estudio H/V4 realizado en Geopsy en base a los siguientes parmetros:
Ventana: 20 segundos
Solape: 0 %
Filtro preprocesamiento: Cosine taper width 0.1
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
Pei.spec: Archivo que contiene el espectro de la componente E-O calculado
con Geopsy utilizando los parmetros que se describen ms abajo.
Pni.spec: Archivo que contiene el espectro de la componente N-S calculado
con Geopsy utilizando los parmetros que se describen ms abajo.
Pzi.spec: Archivo que contiene el espectro de la componente V calculado con
Geopsy utilizando los parmetros que se describen ms abajo.
Parmetros utilizados para calcular el espectro por componentes utilizando Geopsy:
Ventana: 20 segundos
Solape: 0 %
Filtro preprocesamiento: Butterworth orden 3, 0.1Hz-15Hz
Filtro suavizado: Konno Ohmachi, constante = 40, use cosine taper 5 %
Horizontal components: square averaged
Output: 0.01-15Hz con 4000 muestras y en modo lineal
Parmetros utilizados por equipo de obras civiles para clculo de espectro H/V:
Ventana de anlisis: 20 segundos

108

Protocolo de Procesamiento
Solape: 0 %
Filtrado de pre-procesamiento: principalmente Butterworth y muy pocas veces
Chebychev, pasabanda [0.05 10] Hz. A veces frecuencia de corte superior igual
a 20 Hz.
Filtrado de suavizado del espectro: triangular.
Rango de frecuencias naturales encontradas en Valdivia: 0.7 Hz - 10Hz
Anlisis a realizar sobre cada registro:
1. Cargar la seal eoci.txt, valdiviai.txt, Pei.spec, Pni.spec, Pzi.spec, geophv1_i.hv,
geophv2_i.hv, geophv3_i.hv y geophv4_i.hv.
2. Graficar en una misma figura con distintos colores y con leyenda eoci, geophv1_i,
geophv2_i, geophv3_i y geophv4_i. La abscisa debe estar entre 0.01 y 15 Hz.
Esta figura se debe guardar como comparacin_oea_geop_i.fig (donde i corresponde al nmero del registro).
3. Filtrar cada componente temporal con un filtro pasabandas ideal con frecuencias de corte 0.2Hz-50Hz.
4. Calcular espectros por componente utilizando el clsico anlisis de Fourier por
FFT. Estas seales sern llamadas Fe para el espectro E-O, Fn para el espectro
N-S y Fz para el espectro vertical.
5. Graficar para cada componente (tres figuras distintas) Pei vs Fe, Pni vs Fn y
Pzi vs Fz, con dos colores distintos y con leyenda. La abscisa debe estar entre
0.01 y 15 Hz.

Estas figuras deben guardarse como comparacin_espectros_fft_geop_e_i.fig,


comparacin_espectros_fft_geop_n_i.fig y comparacin_espectros_fft_geop_z_i.fig.
Donde e, n, y z corresponden a las componentes E-O, N-S y V respectivamente y la variable i corresponde al nmero de registro.
6. Calcular las STFT (funcin spectog) por componentes utilizando los siguientes
parmetros (estas seales sern llamadas tfe, tfn y tfz para las componentes
N-S, E-O y V respectivamente):
Caso en que la frecuencia de muestreo es de 50 Hz: nfft=1000, noverlap=950.
Caso en que la frecuencia de muestreo es de 100 Hz: nfft=2000, noverlap=1800.

109

Protocolo de Procesamiento
7. Graficar cada STFT por componentes utilizando imagesc: 8. Calcular las AR variante en el tiempo (funcin spectogyule en MatLab) por
componentes utilizando los siguientes parmetros:
Caso en que la frecuencia de muestreo es de 50 Hz: nfft=1000, noverlap=950.
Caso en que la frecuencia de muestreo es de 100 Hz: nfft=2000, noverlap=1800.
Se trabajarn con orden 50 y orden 100 en el caso de las seales cuya
frecuencia de muestreo sea de 50 Hz y con orden 500, cuando la frecuencia
de muestreo sea de 100 Hz.
en caso de orden 50 las variables de salida son are50, arn50 y arz50
respectivamente para cada componente).
en caso de orden 100 las variables de salida son are100, arn100 y
arz100 respectivamente para cada componente).
en caso de orden 500 las variables de salida son are500, arn500 y
arz500 respectivamente para cada componente).
9. Graficar cada espectro de AR variante en el tiempo por componentes utilizando
imagesc.
10. Calcular Nakamura AR (con spectogyulenakamura.m) utilizando los siguientes
parmetros:
a) Nakamura AR es un calculo que se realiza sobre las componentes del
registro de ruido ambiental, las cuales son sei.txt, sni.txt, sz.txt, en primer
lugar se les da un filtraje que va a depender de la frecuencia de muestro,
luego se calcula el espectro en cada ventana y finalmente entrega como
salida sp, f y t, en donde sp es el vector que contiene los datos del espectro
calculado, f es la frecuencia y t es el tiempo.
Caso en que la frecuencia de muestreo es de 50 Hz: nfft=1000, noverlap=950.
Caso en que la frecuencia de muestreo es de 100 Hz: nfft=2000, noverlap=1800.
Se trabajarn con orden 50 y orden 100 en seales cuya frecuencia de
muestreo sea de 50 Hz, y con orden 500 en las seales cuya frecuencia de
muestreo sea de 100Hz.

110

Protocolo de Procesamiento
en caso de orden 50 la variable de salida es sp50.
en caso de orden 100 la variable de salida es sp100.
en caso de orden 500 la variable de salida es sp500.
11. Graficar Nakamura AR utilizando el comando de MatLab imagesc
Conclusiones a extraer de cada experimentacin:
Posibles diferencias entre los espectros H/V de Figura comparacin_oea_geop_i.fig
Posibles diferencias entre espectros por componente calculados por Geopsy y
espectros por componente calculados por FFT. La pregunta a responderse es:
es realmente el espectro calculado por Geopsy un suavizado del espectro de
Fourier?
Verificar si espectros por componente de Geopsy coinciden con los espectros
entregados por spectog y con los espectros entregados por spectogyule.
Comparar los espectros de la figura comparacin_oea_geop_i.fig con los entregados por el Nakamura AR. Coincide la frecuencia natural observada?
Qu componentes importantes frecuenciales aparecen en uno y otro espectro?
Correlacionar resultados obtenidos con nivel de incertidumbre de espectro H/V
obtenido tanto por Geopsy como por el equipo de obras civiles.

111

C. Aclaracin de Conceptos
C.1.

Licuefaccin

La licuefaccin de suelo describe el comportamiento de suelos que, estando sujetos a


la accin de una fuerza externa (carga), en ciertas circunstancias pasan de un estado
slido a un estado lquido, o adquieren la consistencia de un lquido pesado. Es un
tipo de corrimiento, provocado por la inestabilidad de un talud. Es uno de los fenmenos ms dramticos y destructivos y, adems, ms polmicos y peor explicados
que pueden ser inducidos en depsitos por acciones ssmicas.
Debido a las enormes magnitudes de los efectos destructivos de la licuacin durante
el terremoto de Niigata, Japn en 1964, los ingenieros se vieron en la obligacin de
brindarle toda su atencin. Desde entonces se ha avanzado a pasos agigantados para
entender el fenmeno, sus consecuencias, analizando y evaluando el potencial de
licuacin de un lugar y desarrollando nuevas tecnologas que mitiguen las acciones
destructoras de un terremoto [Lop].
El conocimiento del proceso y sus efectos se ha basado en tres diferentes principios:
Observaciones de campo durante y despus de los terremotos.
Experimentos en el laboratorio en muestras de suelos saturados, y en modelos de
fundaciones y estructuras.
Estudios tericos.
Algunos de los efectos de licuacin del suelo son catastrficos, tales como la falla de
grandes pendientes o presas, el desplome de edificios y puentes, el colapso parcial
o total de muros de contencin. Otros sin embargo son un poco menos dramticos
como largas deformaciones de la superficie terrestre, el asentamiento y consecuente
inundacin de grandes reas, etc. An estos efectos laterales en muchos terremotos
causan daos de grandes proporciones a carreteras, vas frreas, tuberas y edificios.

112

También podría gustarte