Bmfcis 218 e
Bmfcis 218 e
Bmfcis 218 e
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.
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 . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . .
reas de geologa.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
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
. .
.
. .
. .
. .
. .
. .
. .
. .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
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.
. 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
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.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
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
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.
1.2.1.
Problematica General
11
1.2.2.
Problematica Especfica
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
13
1.3 Objetivos
1.3.2.
Objetivos Especficos
14
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
15
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
17
2.1.3.
2.1.3.1.
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
19
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
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.3.3.
Mapa de isoperiodos
22
2.2.
Mtodo de Nakamura
23
SV S
SV B
(2.1)
SHS
SHB
(2.2)
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.4)
SHS
SV S
(2.5)
2.2.1.
25
2.2.1.1.
26
2.2.2.
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
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
28
29
2.2.2.2.
30
2.3.
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
32
33
3.1.
34
35
36
36
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
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.1.2.
Filtro de Butterworth
| H() |2 =
1
2N
1+
c
(3.3)
39
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
40
3.1.2.1.
= S(t)(t)
S(t)
(3.4)
S()
= S() ()
(3.5)
41
3.1.2.2.
42
Ventana Rectangular
V (n) = 1
(3.6)
43
3.1.2.4.
Ventana Hanning
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.
44
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
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
x 10
1.5
0.5
10
15
3.2.
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.1.
Anlisis en Tiempo-Frecuencia
3.2.1.1.
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
ST F Tx (t, f ) =
(3.9)
49
200
Tiempo seg.
400
600
800
1000
1200
5
10
15
Frecuencia Hz
20
25
50
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:
(3.10)
51
(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
(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
53
4. Captulo IV MATERIALES Y
MTODOS
4.1.
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
Equipo de adquisicin
Pitrufquen
14
Geosig
Capitn Pastene
Geosig
Mininco
Geosig
Purn
Geosig
Valdivia
60
Geosig
Geosig y Tromino
Lumaco
Total
42
60
54
4.2.
Software Utilizados
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
55
4.3.
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
Ri Purn
Ri Lumaco
Ri Pitrufqun
Ri Mininco
Ri Valdivia
R5 NL Valdivia
Sin Determinar
eoc
Cuadro 5.1.: Nomenclatura de las tablas que poseen el detalle de los registros de
microfizacin ssmicos
57
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
58
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
Registro
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
Registro
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
R4 Lumaco
50
4,077
3,7
R5 Lumaco
50
10,205
10,05
Registro
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
59
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
R702 Valdivia
100
0,732
0,74
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
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
60
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
R20 Valdivia
100
1,135
1,15
1,68 y 3,65
R22 Valdivia
100
1,587
1,71
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
R50 Valdivia
100
1,147
1,4
R51 Valdivia
100
1,257
1,22
Sin Determinar
R52 Valdivia
100
1,135
1,24
R55 Valdivia
100
0,83
0,83
1,5
R56 Valdivia
100
0,842
0,823
1,2-1,6
61
Figura 5.1.: Distribucin porcentual de todas las seales procesadas en nuestro estudio
62
5.1.
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
63
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
R35 Valdivia
100
0,781
0,808 y 0,984
1,36-3,23
R48 Valdivia
100
1,062
0,98
R50 Valdivia
100
1,147
1,4
R51 Valdivia
100
1,257
1,22
Sin Determinar
R52 Valdivia
100
1,135
1,24
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.
5.2.
64
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
R3 Lumaco
50
8,435
7,06 y 8,09
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
Figura 5.5.: Porcentaje de seales que poseen un espectro diferente a las dems
pruebas de Geopsy
5.3.
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
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
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
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
Figura 5.7.: Porcentaje de seales con peak truncado en comparacin al peak mas
definido de eoc
5.4.
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
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
68
5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy
3.5
2.5
1.5
0.5
10
11
35
30
25
20
15
10
0
0
0.5
1.5
2.5
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
3.5
2.5
1.5
0.5
0
0.5
1.5
2.5
70
5.4 Estudio estadstico comparativo entre los resultados de Obras Civiles y Geopsy
60
50
40
30
20
10
0
0
10
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
71
3.5
2.5
1.5
0.5
10
11
5.5.
72
5.5.1.
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
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
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
200
400
2,8 Hz
600
3,46 Hz
800
5,23 Hz
1000
10
15
74
Figura 5.20.: mesh del espectro obtenido con modelamiento AR variante en el tiempo
5.5.2.
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
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
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
1000
800
600
400
200
Banda 1
Banda 2
10
15
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.3.
78
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
11,35 Hz
0,84 Hz
0
10
15
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
4,93Hz5,47Hz
7,59Hz12,44Hz
200
400
600
800
1000
1200
5,73Hz7,08Hz
10
15
5.5.4.
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
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
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
81
200
3,69 Hz
400
600
800
1000
1200
10
15
200
3,71 Hz
400
600
800
1000
1200
10
15
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
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
1000
800
3,69 Hz
600
400
200
10
15
83
5.6.
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
84
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
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
85
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
1,04 Hz
0,72 Hz
3
0.5
1.5
2.5
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]
[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]
[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]
[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]
[DGA88]
[eA07a]
91
Bibliografa
[eA07b]
[eA08]
[Efr03]
[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]
[F.10]
[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]
[Gab46]
[Geoa]
[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]
[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]
[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]
[JP07]
[JWH02]
[Kay99]
S.M. Kay. Modern Spectral Estimation: Theory and Application. Prentice Hall, 1999.
[Kit10]
[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]
[LCG94]
[Les09]
P. Lesage. Interactive Matlab software for the analysis of seismic volcanic signals. Computers and Geosciences, Vol. 35, pp. 2137-2144, 2009.
[LESA]
[LGM02]
[Lop]
Clara Eugenia Giraldo Lopera. Licuacin del suelo. Mecinica de Fluidos y Recursos Hidriulicos, Escuela de Ingenieria de Antioquia.
[lW]
[Mal08]
[MGC99]
[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]
[Mig10]
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]
[Nak89]
[Nak08]
[Oba]
[oSR89]
[Pil07]
[Pro01]
[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]
[Raf07]
[RL10]
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]
[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]
[SESb]
SESAME.
Descripcin del programa computacional GEOPSY.
http://www.geopsy.org/wiki/index.php/Main_Page.
[SL05]
[SLM10]
[SM05]
[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]
[Tra10]
[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]
[ZTC03]
97
A. Planes de trabajo y
sectorizaciones de las zonas de
estudios
A.1.
98
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.
99
100
A.1.2.
101
A.2.
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.
A.4.
103
A.5.
A.6.
104
A.7.
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.
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
112