UNIVERSIDAD NACIONAL DE SAN AGUSTIN

MAESTRIA EN CIENCIAS DE LA TIERRA, MECION GEOTECNIA

GEOFISICA APLICADA
NOMBRE: Wilbet Mamani Hañari DNI 43695735
mahawilgl@gmail.com

RAZONES ESPECTRALES H/V

1 INTRODUCCION

En el curso de geofísica aplicada del post-grado se tiene una fase teórica durante el cual se dan a conocer conceptos sobre métodos
sísmicos con las cuales se aplican a la geotecnia como métodos indirectos. En la cual se pone en práctica el método de RAZON
ESPECTRAL H/V, con datos sísmicos para la zona del Mato Groso – Brasil.

En esta práctica se considera el uso de registros de ruido superficiales, que corresponden a mediciones de velocidad en tiempo
continuo en todas las dimensiones del espacio, las dos componentes horizontes Norte-Sur (N-S) y Este-Oeste (E-O) y la componente
vertical (Z).

Como lo muestra la Figura 1.1 en la vista aérea de Mato Groso las estaciones abarcan un tramo en estudio. Los sismógrafos están
ubicados sobre. La distancia promedio entre estaciones aledañas es de 5 Km abarcando 198 Km a lo largo de la vía, barriendo una
superficie de 617 Km2 aproximadamente.

Figura 1.1 ubicación de los puntos de registros de ruido sísmico.

 2 METODOLOGIA

2.1 INSTRUMENTOS QUE PUEDEN SER UTILIZADOS

Existen diferentes tipos de estaciones sísmicas: sismómetros de periodo corto, sismómetros de banda ancha, acelografos, etc. En esta
memoria se trabajó con sismómetros de banda ancha, los cuales presentan una respuesta instrumental lineal para el rango de
frecuencias que interesa estudiar en esta memoria.

Cada estación utilizada en esta memoria está compuesta por 4 equipos: un digitalizador, una unidad registradora de datos, sensores y
una antena GPS.

Cada sismómetro de cada estación cuenta con un digitalizador Quanterra Q330 (Figura 2.2). Este tiene 3 ó 6 canales de alta
resolución y es de baja potencia. El Q330 puede escribir datos en cuatro sistemas receptores diferentes al mismo tiempo. En los
controles de sensores, se tiene una función automática de recentrado de masas y las funciones de calibración.

Los sismómetros son Trillium Compacto 120s (Figura 2.2), que corresponden a un sismómetro de banda ancha con un geófono
resistente. Este instrumento incorpora un diseño de un sensor de retroalimentación de fuerza triaxial simétrica con una respuesta plana
a la velocidad de 120 segundos.

Cada sismógrafo tiene una frecuencia de muestreo de 100 [muestras/s] registrando en tiempo continuo, desde su instalación, la
velocidad del sitio en las 2 componentes horizontales Norte-Sur (N-S) y Este-Oeste (E-O) y la vertical (Z). El registro se entrega en
cuentas y la constante para obtener m/s tiene un valor de 5.017 x 108[cuentas/m/s].

Figura 2.2 Digitalizador Quanterra 330 a la izquierda y Sismómetro Trillium Compacto 120s a la derecha

2.2 INSTALACIÓN DE ESTACIONES

Al momento de elegir el lugar donde se situarán las estaciones se debe buscar que ésta permanezca lo más alejada de cualquier
perturbación, debido a la gran sensibilidad de los instrumentos. La Figura 2.3 se muestra el sitio de una de las estaciones en donde se
puede ver la forma de posicionar los sensores in situ, además de la caja con las herramientas necesarias para obtener las señales.

Figura 2.3 Instalación de la Estación una excavación superficial tipo para dejar los sensores en una ubicación permanente, en
(b) se presenta los sensores en la excavación y en (c) se observa la excavación completamente cubierta. Por ultimo en la (d),
se puede ver la caja con el digitalizador y demás instrumentos que transforman las vibraciones del suelo en los registros.

Páginá 2
 2.3 RECOPILACIÓN DE LOS REGISTROS

Existen dos formas de obtener los datos de las estaciones: estas son copiando los datos del digitalizador o recibirlos mediante señal on
line, estas señales son ingresadas aun base de almacenamiento de donde se copian los registros para ser usados. Una vez obtenidos
los datos se deben transformar de formato Antilope a formato SAC. Se eligió este formato por ser un tipo legible por Geopsy programa
con el cual se analizaran las señales.

2.4 REGISTROS DEL ARREGLO

Luego de tener los registros en formato SAC, se procede a filtrar la señal entre 0.1 y 20 [Hz], para evitar la influencia de periodos
largos. Además las señales son corregidas por tendencia y promedio.

En las Figura 2.4 se tiene un ejemplo de registro correspondiente ente las 14:18:27 pm y las 14:31:07, del día 10 de Junio del 2017, en
donde se muestran entradas las 3 componentes, Vertical (Z), Este-Oeste (E-O) y Norte-Sur (N-S) sin filtrar ni removida la tendencia y
filtrado y con la tendencia removida respectivamente. Dependiendo del análisis que se quiera realizar, el registro se puede cortar para
la ventana temporal necesaria,

Figura 2.4 Registros del día 10 de junio del 2017 de la Estación 141827 sin Filtrar ni Corregir (arriba) y con ventanas temporal
necesaria como filtro (abajo).

Páginá 3
 3 RAZÓN ESPECTRAL SUPERFICIAL ENTRE LAS COMPONENTES HORIZONTAL Y VERTICAL (H/V)
USANDO EL SOFTWARE GEOPSY®.

3.1 MARCO TEORICO

Al reconocer las ondas Rayleigh correspondidas con un movimiento elíptico retrogrado, se aprecia en todas las componentes, por lo
que la amplificación de su efecto (H y V) puede ser cuantificado con la razón entre los espectros (A) de los movimientos de la superficie
(S) y del basamento (B) del estrato siguiente (Nakamura, 1989) se tiene:

Luego la razón H/V

Dado que el basamento corresponde generalmente a roca dura, ésta se comporta como suelo duro por lo que el cociente
que refiere al basamento tiende a 1.

Por lo que

3.2 PROCEDIMIENTO PARA OBTENER EL H/V

Geopsy® es un programa que permite calcular las razones espectrales de las señales ingresadas de una manera bastante simple para
el usuario, sin embargo hay que definir algunos inputs para poder analizar correctamente los registro y así obtener resultados
razonables, dado que las razones espectrales se calculan a partir de los espectros, se deben tener las mismas consideraciones para
su cálculo.

Uno de los parámetros principales es determinar las ventanas temporales de la señal que serán utilizadas y el tamaño de estas. Una
ventana se define como la multiplicación de una parte de la señal por una función suave y nula en los extremos, y así evitar aparición
de altas frecuencias producidas por el corte abrupto de la señal al inicio y al final. En cada ventana se calculan los espectros y luego las

Páginá 4
 razones espectrales, luego se calcula una curva promedio y un error asociado. Se debe tener cuidado que el largo de la ventana debe
ser acorde al rango de frecuencias que se busca interpretar, en este caso la frecuencia mínima mostrada es de 0.2 [hz] lo que
representa un periodo máximo mostrado de 5 [seg] por lo que se eligió un largo de ventana Lw de 60 [seg].

El criterio de selección de ventanas de ruido blanco se hace en definición los términos STA y LTA. STA es la amplitud promedio de la
señal en un período corto de tiempo, tsta, y LTA es la amplitud promedio de la señal en un período más largo de tiempo, tlta. El proceso
de selección de ventanas se basa en la comparación de la razón STA/LTA con valores umbrales máximo, Smax, y mínimo, Smin,
definidos de antemano, al igual que un porcentaje de traslape (Pasten, 2007).

Para seleccionar una ventana de vibraciones naturales estacionarias se toma una parte de la señal de largo tlta, en este caso igual al
largo de la ventana Lw, y se calcula LTA. Posteriormente una ventana móvil de largo tsta recorre el intervalo definido por tlta de
principio a fin y va calculando instantáneamente STA. Si la razón STA/LTA se enmarca dentro de los valores umbrales, Smax y Smin,
durante todo el recorrido de la ventana móvil, esa parte de la señal será considerada como una ventana apropiada para el cálculo de
razones espectrales de la segunda etapa (Pasten, 2007).

En la Tabla 3.1 se pueden observar los criterios ingresados, en la Figura 3.1 las ventanas seleccionadas para un registro típico y la
Figura 3.2 muestra el resultado de la razón espectral.

Tabla 3.1 Criterios de selección de ventanas ingresados a Geopsy

Parámetro Símbolo Valor

Largo de ventanas LV 60 [s]
Tiempo de la ventana corta STA tSTA 1 [s]
Tiempo de la ventana Larga LTA tLTA 60 [s]
Mínimo valor del cociente STA/LTA (STA/LTA)min 0.2
Máximo valor del cociente STA/LTA (STA/LTA)max 2.5
Coeficiente de ancho de banda b 40

Figura 3.1 Ventanas seleccionadas para el registro. Cada color representa una ventana de 60 s.

Páginá 5
 Figura 3.2 Razón espectral resultante de la señal mostrada en la Figura 3.1 proporcionada por Geopsy, la curva continua
representa el promedio, mientras de las segmentadas la dispersión producto de calcular la razón espectral para todas las
ventanas seleccionadas.

4 RESULTADOS.

La vibración dinamica de un suelo producto de un sismo presenta con una amplificación de las ondas sísmicas caracterizada por las
propiedades del estrato por donde se propagan las ondas (Kanai, 1956, 1957). El argumento para encontrar el periodo fundamental
usando el método de Nakamura es la consideración de que las fuentes generadoras de ruido o vibraciones naturales y sismos
producen una solicitación dinámica aleatoria, por lo cual un depósito de suelo tiende a vibrar mayormente de acuerdo a su periodo
fundamental, que se evidenciada en los peak de las curvas de los H/V que muestra que dicha frecuencia tiene un factor de
amplificación mayor (Lachet et al,, 1996).

Las diferencias entre las formas y amplitudes de los registros depende del tipo de suelo, por lo que es de esperar que un suelo menos
rígido muestre un H/V con un peak y un suelo más rígido tenga una razón espectral plana para todo el rango de frecuencias estudiado
(Bonnefoy-Claudet et al, 2009).

En la Figura 4.2 se representa las frecuencias predominantes para cada punto de medición.

4.1 CLASIFICACIÓN DE LOS H/V

Dadas las formas de las curvas de la razón espectral H/V obtenidas para los sitios de las estaciones, se proponen cuatro grupos (ver
Figura 4.1):

 Planos, aquellos con amplitudes menores a 2 en donde no es posible identificar un peak en la razón espectral H/V
 Con Peak Somero, aquellos en donde la razón espectral H/V es mayormente plana, pero es posible identificar un peak
predominante de amplitud cercana a 2
 Con Peak Marcado, aquellos en donde se identifica uno o más un peak predominantes en donde alguno alcanza amplitud
superior a 2, pero inferior a 3
 Con Peak muy Marcado, aquellos en donde se identifica uno o más un peaks predominantes en donde alguno alcanza una
amplitud superior a 3

Páginá 6
 Espectro Pláno Espectro con
Peák Somero

Espectro con
Espectro con Peák
Peák Márcádo
muy Márcádo

Figura 4.1 Razones Espectrales. La curva de color negro es el H/V promedio.

4.2 PERIODOS FUNDAMENTALES

Como se expuso en la introducción del capítulo, relacionando el valor espectral máximo de la curva H/V se puede obtener la frecuencia
natural de oscilación y con el recíproco de ella, el periodo fundamental del sitio. Es importante recalcar que para las curvas planas no
se puede distinguir una frecuencia fundamental por lo que para este tipo no se calculó un periodo natural.

En las Tablas 4.1 a la 4.3 se muestran los resultados obtenidos para los distintos tipos de razones espectrales.
En la figura 4.3 se observa el tipo de razón espectral para cada punto de medición.

Páginá 7
 Tabla 4.1 Periodos Naturales H/V con Peak Somero

Código E N f (Hz) H/V max Tn (S)

134450 475196 8611242 6.75 1.52 0.15
174841 463897 8639125 1.34 1.53 0.75
170732 438681 8701186 2 1.61 0.50
134611 428155 8718265 2.23 1.65 0.45
183520 442186 8691724 1.01 1.68 0.99
160703 478210 8591865 2.08 1.68 0.48
115135 447736 8677796 1.84 1.71 0.54
154922 464151 8648600 4.39 1.74 0.23
122825 449499 8673099 0.61 1.78 1.64
142225 474326 8606398 1.08 1.78 0.93
153231 476914 8596689 0.68 1.81 1.47
134034 429790 8740711 1.58 1.83 0.63
130959 473400 8615910 0.91 1.87 1.10
130108 452917 8669434 1.42 1.88 0.70
112527 467852 8629823 1.68 1.91 0.60
123303 471571 8620588 1.68 1.91 0.60
125637 426932 8723109 9.74 1.92 0.10

Tabla 4.2 Periodos Naturales H/V con Peak Marcado

Código E N f (Hz) H/V max Tn (S)

144828 436715 8748428 3.46 2.01 0.29
152344 441202 8750794 5.78 2.01 0.17
133323 456392 8665773 3.99 2.25 0.25
110850 454898 8724702 2.63 2.28 0.38
122040 428271 8727917 0.62 2.4 1.61
130124 432982 8736794 1.95 2.57 0.51
182453 465890 8634465 0.65 2.6 1.54
145602 475596 8601525 0.72 2.66 1.39
141827 432779 8744650 3.02 2.68 0.33
140508 459856 8662121 1.25 2.82 0.80

Tabla 4.3 Periodos Naturales H/V con Peak Muy Marcado

Código E N f (Hz) H/V max Tn (S)

122143 431073 8732118 1.09 2.99 0.92
164536 479748 8587172 0.74 3.04 1.35
171419 462382 8643925 1.54 6.86 0.65

Páginá 8
 Figura 4.2 Razones Espectrales. Con el Rango de Frecuencias.

Páginá 9
 Figura 4.3 con los tipos de H/V obtenidos. Los marcadores indican la posición de las estaciones y el color indica el tipo de
razón espectral. Las zonas coloreadas indican la unidad geológica descrita por el instituto geológico de Brasil.

Páginá 10
 5 DISCUSIÓN Y CONCLUSIONES

A partir de la posición de las estaciones en el mapa (figura 4.3), se puede apreciar una clara correlación entre los depósitos más rígidos
con razones espectrales sin un peak marcado lo que muestra que entre esto suelos y el basamento hay un menor contraste. En otras
palabras, al realizar la razón espectral H/V en la gravas, se obtiene un espectro plano o con un peak somero.

Las estaciones que muestran un peak muy marcado y que alcanzan a ubicarse dentro del mapa están claramente enmarcadas en
depósitos de suelos más finos y por tanto de menor rigidez por lo que es de esperar que en estos sitios se tenga una mayor
impedancia de las ondas sísmicas y que esta propagación incorpore el marcado periodo natural de vibración del sitio.

Las curvas resultantes de la razón espectrales siguen algunos patrones por lo que se propone una clasificación usando la forma de las
curvas de ellas. Por ejemplo las razones planas, se mantienen en un vecindario bastante pequeño de amplitudes espectrales H/V igual
a 1, no hay razones planas en torno a otro valor, lo que implica que la componente horizontales, no sólo es similar en forma sino
también en tamaño respecto de las componentes superficiales, lo que se condice con que en este tipo de razones no exista un periodo
predominante. Asimismo las razones espectrales que muestran un peak somero, permanecen predominantemente en la unidad, pero
con máximos visibles, que casi no superan H/V igual a 2.

6 BIBLIOGRAFIA

Nakamura,Y. “A method for dynámic chárácteristics estimátion of subsurfáce using microtremor on the
ground surfáce” Quárterly report of ráilwáy technicál reseárch institute (RTRI), 30(1), 25-33. Jápán.
1989.

Pardo, M. (2012). Apuntes del curso Dinámicá de Suelos. Depártámento de Ingenieríá Civil.
Universidád de Chile.

Woolery, E W and Street, R: “3D neár-surfáce soil response from H/V ámbient-noise rátios”. Soil
Dynámics ánd Eárthquáke Engineering, 22 (9), 865-876. 2002

Páginá 11