5.2 - REDES SÍSMICAS -...

25
5.2 - REDES SÍSMICAS

Transcript of 5.2 - REDES SÍSMICAS -...

Page 1: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

5.2 - REDES SÍSMICAS

Page 2: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

5 - REDES SISMICAS

5.1 - Redes de cobertura mundial, regional y local+Arreglos5.2 - Redes sísmicas de México5.3 - Boletines de sismicidad5.4 - Catálogos de sismicidad y su utilidad5.5 - Criterios de selección de instrumentación sísmica5.6 - Estudios de niveles de ruido

Page 3: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Tamaño o apertura: máxima distancia entre 2 estaciones

Diferencia red - arreglo: técnicas usadas en procesado (red ↔ arreglo ↔ estación). “una red puede utilizarse como un arreglo o un arreglo como una red “.

ARREGLOS: REDES SÍSMICAS CON ALTA COHERENCIA ENTRE SEÑALES

ARCES (N Noruega):4 anillos concéntricos 150, 325, 700 y 1500 m

(3) (5) (7) (9)24 SP Z + BB centro

Page 4: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

MÉTODOS DE ARREGLOS

•Método establecidos. Por datos nuevos es muy “de moda”.

•Desarrollado para hacer detecciones y localizaciones muy precisas de explosiones nucleares.

• Unos nuevos aplicaciones son estudios de temblores grandes y localización de tremor no-volcánicos.

Page 5: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Ventajas:- detección señales muy débiles, mal identificables 1 sitio (detección nuclear)- mejora localización epicentral: ‘orientación’ del array como antena (prueba ≠ vapp y Φ) - identificación de fases (≠ vpropagación)

Condiciones:- r >> 20·Øarreglo aproximación ondas planas- r > 10·λ “ “ “- r suficientemente pequeña comportamiento no puntual- alta precisión medida t relativos entre estaciones - uso técnicas específicas:

- filtrado en velocidad o formación de haz (i: beam forming)- apilado (i: stacking)- análisis f-k

supresión ruido y mejora ratio s/n

MÉTODOS DE ARREGLOS: PARA AMPLIFICAR LA SEÑAL Y DETECTAR LA DIRECCIÓN

NORSAR (Noruega):60 km Ø42 sitios

7 subarreglos: 6 SP Z, ~ 3 km Ø

Page 6: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Yellowknife (Canadá):SP (azul y rojo) + BB (verde)

↑↑ volumen datos tratamiento automático - detección (filtros STA/LTA, formación de haz) - procesado - atributos de la señal (f-k, tllegada, T, amplitud, polarización)

Page 7: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Conceptos básicos:

rj: vector posición sitio j respecto sitio referencia

Φ: azimut (realmente back-azimut)Θ: dirección de propagación del frente de ondas resp. N Θ = Φ ± 180º

i: ángulo incidencia (≤ 90º)vapp: velocidad aparente con que frente ‘barre’ arreglo

[vc, ∞) -↑,→

f(i, vc)

s = 1/vapp: lentitud (constante para un rayo específico)

s/km local o regional s/º telesísmico –parám. rayo: 1/(vapp·p), p = 6371π/180º ≈ 111.19 km/º-

k = ω·s = ω/vapp = 2π/λ número de onda (km-1)

τj: tiempo de retraso de sitio j respecto sitio de referencia

>0: llegada antes a j que a sitio referencia

Page 8: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

τ2 = t2-t1 = L/vc retraso sitio 2 resp. a 1vapp = d/(t2-t1) = vc/sen ii) Si diferencias altitud entre sitios << no corrección (sup. vapp,z = ∞ y sz = 0)Para sitio j (xj,yj):

ii) Si diferencias altitud entre sitios imp.(p.ej., pozos)

corrección:

El tiempo de retraso entre estaciones depende de la velocidad aparente y el back azimut

Page 9: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Apilado (stacking)

si no existe gran atenuación local entre sitios arreglo distancia entre sitios próximos sftem. pequeña registros en fase sftem. grande ruido no común en cada sitio

ruido más incoherente que la señal mejora ratio s/n por suma (apilado)

señal observada: w(t) = S(t) + n(t)

obteniendo el registro suma o haz (beam) de las trazas de los M sitios del arreglo:

Page 10: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Apilado (stacking)

Suponiendo n(t) distribución normal de amplitudes valor medio 0 varianza = σ2 para todos los sitios

Al sumar las trazas la varianza del ruido queda σs2 = M·σ2 σs

= √M · σ

ruido apilado α √M señal apilada α M ganancia o mejora en la ganancia G2 = M

vapp = 10 km/s, Φ = 158º, evento en Grecia fase PcP sismo mar Tierreno, 9.6º

Page 11: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Una de las mayores dificultades: detección de señales diferentes del ruido de fondo

uso algoritmo STA/LTA

Page 12: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Función de transferencia de un arreglo:

describe sensibilidad y resolución arreglo = f (sfase observada, kfase observada, geometría arreglo)

arreglo óptimo para detectar señales con lentitud s0 (señales con otra s: supresión parcial)

Influencia de los distintos parámetros del arreglo:apertura: define resolución para k pequeños + apertura menor k medible (mayor λ)

λ máxima analizable ~ apertura

nº sitios arreglo (M): controla habilidad arreglo para suprimir energía que cruza arreglo al mismo tiempo que señal y con diferente s (≡ filtro en k)

distancia entre sitios: define mayor k que puede resolverse (menor λ)

menor distancia menor k para una velocidad dadageometría: define dependencia de los puntos anteriores con azimut

Page 13: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

ej. geometría muy diferente función transferencia muy distinta distinta resolución a la s (lentitud) de un frente

Yellowknife: gran apertura alta resolución para medir vapp

pobre resolución azimutal

ARCES: pequeña apertura incapacidad para resolver ondas con peq. difs en k resolución azimutal perfecta

Page 14: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Análisis f-k

dominio f

empleado para estimar s de una fase

red de 51x51 puntos equiespaciados de s

entre -0.4 y 0.4 s/km

para cada punto se evalúa potencia del haz:

máxima potencia define s de la fase

Page 15: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

ARREGLOS Y ESTRUCTURA

Arreglos 2D (lineares) o 3D

más útiles donde los cambios en estructura son más grande en una dimensión que otra

usados en experimentos de reflexión y refracción

MASE

Iglesias et al 2010

Page 16: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

TEMBLOR DE SUMATRA 2004

• Empezar con una localización de prueba

• Alinear los registros

• Estimar la coherencia entre sismogramas para la localización de prueba para una ventana de tiempo pequeña (aquí 30 segundos).

• Repetir para todas las localizaciones en un mapa.

• Desafortunadamente funciona mejor para temblores muy grandes.

• Difícil interpretar (qué significan los colores?)

• En el caso de Sumatra 2004, los colores fueron interpretados como deslizamiento

Ishii et al 2007

Page 17: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

TEMBLOR DE SUMATRA 2004

• Empezar con una localización de prueba

• Alinear los registros

• Estimar la coherencia entre sismogramas para la localización de prueba para una ventana de tiempo pequeña (aquí 30 segundos).

• Repetir para todas las localizaciones en un mapa.

• Desafortunadamente funciona mejor para temblores muy grandes.

• Difícil interpretar (qué significan los colores?)

• En el caso de Sumatra 2004, los colores fueron interpretados como deslizamiento

Ishii et al 2007

Page 18: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

TEMBLOR DE SUMATRA 2004

• Empezar con una localización de prueba

• Alinear los registros

• Estimar la coherencia entre sismogramas para la localización de prueba para una ventana de tiempo pequeña (aquí 30 segundos).

• Repetir para todas las localizaciones en un mapa.

• Desafortunadamente funciona mejor para temblores muy grandes.

• Difícil interpretar (qué significan los colores?)

• En el caso de Sumatra 2004, los colores fueron interpretados como deslizamiento

Ishii et al 2007

Page 19: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

TEMBLOR DE SUMATRA 2004

• Empezar con una localización de prueba

• Alinear los registros

• Estimar la coherencia entre sismogramas para la localización de prueba para una ventana de tiempo pequeña (aquí 30 segundos).

• Repetir para todas las localizaciones en un mapa.

• Desafortunadamente funciona mejor para temblores muy grandes.

• Difícil interpretar (qué significan los colores?)

• En el caso de Sumatra 2004, los colores fueron interpretados como deslizamiento

Ishii et al 2007

Page 20: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Ishii, Harvard

TEMBLOR DE JAPON 2011

Page 21: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Meng y Ampuero, Caltech

Page 22: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

Ji, UCSB

SISMOGRAMAS, USANDO MULTIPLE

VENTANAS

Page 23: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

• Los circulos y rectangulos son las localizaciones de los altas frecuencias, uno usando el arreglo de los estaciones en Europa y el otro con los de USArray

• El superficie de colores es el deslizamiento en derivado por datos de GPS.

on

May

25,

201

1w

ww

.sci

ence

mag

.org

Dow

nloa

ded

from

Simons et al 2011

Page 24: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

ESTUDIOS DE TREMOR NO-VOLCANICO

• T

“The million dollar question”: El tremor no-volcanico es causado por (a) deslizamiento en el interface (b) por movimiento de fluidos en la corteza arriba de el interface (c) otroPara contestar es escencial saber las localizaciones de los tremores.

Kostoglodov

Page 25: 5.2 - REDES SÍSMICAS - usuarios.geofisica.unam.mxusuarios.geofisica.unam.mx/vala/cursos/Registros_Sismicos_2013-1... · -uso técnicas específicas: ... c) s = 1/v app: ... app,z

G-GAPParte del proyecto son studios de tremor no-volcanico en la zona de subducción de GuerreroUna localización muy precisa es esencial para entender el proceso físico que causa los tremores

Kostoglodov