INETER - Geofísica - Sismología
Detección y Procesamiento de Sismos
con SEISLOG y SEISAN

(en la Central Sísmica de INETER)

SEISLOG, SEISAN, SEISNET -
son sistemas de software elaborados por la
Universidad de Bergen, Noruega

El sistema SEISLOG (en las máquinas QNX, QNCBIG y QNXBB)

(
descripción amplia de SEISLOG, ver http://www.cepredenac.org/proyects/resis/seisman.html)

Funciona bajo el sistema operativo QNX, que es similar al UNIX
(por lo tanto existen los comandos de UNIX, como ls, cp, rm, ftp, telnet, …)

Tareas

- Digitalización de los canales analógicos
- Grabación continua de la información, en disco duro
- Filtración de los datos
- Detección de eventos con el método STA/LTA, usando diferentes criterios
- Determinación preliminar de la llegada de ondas P
- Almacenamiento de eventos en disco duro

El acceso a la información es posible usando las funciones del sistema operativo, vía red local o línea telefónica (modem). El SEISLOG tiene también posibilidades de procesamiento básico de los datos, p.ej. un programa interactivo para picar fases y un sencillo programa de localización .

Programas importantes del monitoreo de las funciones de SEISLOG:
(en la PC del QNX, QNXBB o QNXBB)

MON - monitorea la adquisición de datos, detección y funcionamiento del sistema de tiempo
  entrar : mon

DPL   - monitorear la señal de un solo canal
 entrar: dpl, número del canal, 1

PAR - Programa para cambiar parámetros de SEISLOG
            uso de este programa solo por personas autorizadas!
   Con este programa se determina cuales de los canales se usan para la detección,
   como se hace la detección, cuales de los canales canales se archivan, etc…

Copiar datos del QNX, QNXBIG, QNXBB a la SUN

Manualmente:

Las máquinas están conectadas vía red local.

Se usa el programa ftp: por ejemplo con la máquina qnx :

ftp -i qnx - conectarse con el qnx, en algunos casos pide el programa un password
pwd saber en que directorio se trabaja
ls ver el contenido del directorio
bi             copiar en modo binario
get nombre_del_archivo
mget *_1 copiar todos los archivos con extensión _1 a la máquina local

Nota: De la misma manera se copian datos de cualquier máquina en el INTERNET:
ftp -i qnxbb.unah.hon o ftp -i número_de_la_máquina
ls - lista de archivos
bi - cambiar a modo binario
mget nombre _del_archivo

Registrar manualmente eventos en la base de datos (después de copiarlos a la SUN)

dirf nombre_archivos     p.ej:dirf 2000*_1,   hacer lista de eventos, con extensión _1
qnxsei cambiar el formato de los datos, del formato del QNX al SEISAN
dirf nombre_evento(s) lista de eventos en nuevo formato
mulplt
0
1                 ver los sismogramas decidir sobre integración en base de datos

Procesar eventos

eev añofechahora p.ej. eev 2000011619 , el evento de 2000/01/16 19horas

Proceso automático para copiar y registrar eventos

La SUN chequea cada minuto si el QNX tiene un nuevo evento registrado.
Si eso es el caso - el sistema
- copia el sismograma digital al directorio WOR/seisnet,
- registra el evento en la base de datos,
- pica fases,
- localiza el sismo,
- determina la magnitud Mc,
- actualiza el mapa de sismos de Nicaragua y de América Central,
- actualiza la lista de eventos de los últimos 10 días,
- copia mapas y listas al servidor de INTERNET.

La localización y la magnitud automática no son muy confiables. Por eso, el sismólogo de turno tiene que revisar el evento inmediatamente!


Procesamiento de los registros sísmicos, en la SUN (con SEISAN)

1. Interpretar evento

local
- frecuencias altas 3..10 Hz, es decir períodos de ondas P de entre 0.3 y 0.1s
- los sismogramas en las estaciones se ven muy diferentes entre si.
- Eventos fuertes se ven muy bien en los acelerogramas MGAAHZ, ..N, ..E

regional
- sismogramas muy largos, frecuencias bajas, 0.3.. 10 Hz, ver en la banda ancha!
- Los registros en las estaciones se ven bastante parecidos.
- Salen bien después de filtrar con x (0.1 .. 1 Hz)
- Se ve poco o nada en los acelerógrafos.

lejos
sismogramas muy largos, registros de las estaciones se ven muy parecidos
La frecuencia es muy baja, es decir períodos de 1 s o más bajo.
Filtrar con x (0.1 ..1 Hz) o z (0.01 .. 0.1 Hz)
Se ven muy bién en la estación de banda ancha.
No se ve nada en los acelerógrafos

volcánico
Se ve bién solo en una estación, la que está instalada en el volcán.
Mucho más pequeño en estaciones más lejanas.

2. Picar fases, con EEV

P - no hay problemas.. Cuidado: picar en ventana de pocos segundos, no en ventana
del evento completo! Para garantizar exactitud.
S - buscar en las componentes horizontales !
eventos normales en MGAN, MGAE, PYTN, TISN, TISE
para eventos fuertes en MGAAHN, MGAAHE (acelerógrafo Managua)
y MGABBN, MGABBN (banda ancha Managua)
para eventos fuertes lejanos en MGABB (banda ancha)
No inventar S donde no se ve nada !
Es suficiente picar una S.

3. Magnitud, con EEV

para eventos locales MC en la coda de todas las estaciones,
y/o ML, especialmente si la coda está cortada !!
para eventos fuertes locales ML !!! en MGAAHN, MGAAHE (acelerógrafo)
y/o MGABBN, MGABBE (banda ancha)
para eventos fuertes lejanos
   MB de la fase P en MGABBZ,
   MS de las ondas superficiales MGABBZ, MGABBN, MGABBE
  se puede probar determinar MB en MGAZ, MGAN, MGAE

4. Localizar, con EEV

eev yymmddhhss (base_de_datos)
L - el programa localiza el evento y presenta epicentro, magnitud, nombre de área ...
- ver los errores en el tiempo de las estaciones.
- ver sismogramas de las estaciones con errores grandes, corregir
- para eventos locales los errores deben estar debajo de 1 segundo.
- Para eventos regionales se puede permitir errores más grandes, inclusive algunos segundos.

5. Hacer mapa del evento
Usar el programa epimap… y/o los comandos:

comando función   archivo del mapa
mga                        Mapa de Managua y alrededores, con epimap   epimap.plt
nic                          Mapa de Nicaragua, con epimap    epimap.plt
cam                        Mapa de Centroamérica, con epimap    epiamp.plt
nica                        mapa bonito de Nicaragua, con     nica.ps
came                      mapa bonito de Centroamérica, con GMT     nica.ps

Los mapas se imprimen con
lp archivo_del_mapa

6. Editar e Imprimir el comunicado con comando comunic

chequear el texto de comunicado!


Comandos del EEV (los comandos se dan en mayuscula o minuscula)

E      - editar el archivo de lecturas,
T      - ver el archivo de lectura (S-file), en pantalla
TT    - ver solo cabecera del archivo de lectura (S-file), en pantalla
P       - plotear evento en pantalla, procesarlo (picar fases, determinar amplitudes, etc..)
Q      - salir del EEV
L       - localizar evento, calcula también magnitud (con el programa Hypocenter)
Lxx   - localiza esta lectura y la con el número xx, juntos, como un solo evento
LL    - localiza esta lectura juntos con la siguiente como un solo evento
xx     - pasar al evento número xx
Dxx  - pasar al primer evento del día xx
Axx  - integra lectura xx en el archivo de la lectura actual, el evento original se mantiene
AA   - integra lectura actual a la próxima lectura, lectura original se mantiene
D      - deletear esta lectura, el programa pide confirmación
B      - va para la lectura anterior
C      - copiar esta lectura a otra base de datos o en archivo del directorio corriente
R      - renombrar el tipo del evento (local, regional, distante)
Z      - picar fases automáticamente, fases ya existentes se quedan pero una marca
F      - calcular una solución del plano focal
FC    - acumular datos para una solución del plano focal de varios eventos
I       - arranque del programa PITSA
Jyymm BAS   - salta al año yy y e mes mm de la base de datos BAS
Eyymm           - terminar la sesión de trabajo del EEV con el año yy y el mes mm
O     - ejecuta un comando del sistema operativo (p.ej. ols ejecuta ls)
Sxxxx - buscar los próximos 2 eventos que se encuentran
            en un período de tiempo de xxxx segundos
MAC  - entrar información macrosísmica
SYNT - preparar parámetros para modelaje sintético
BOUCH - correr el programa de Bouchon para calcular sismogramas sintéticos
BOUSEI - hacer un archivo en el formato de SEISAN usando el resultado de BOUCH
IASP      - generar tiempos de arribo sintéticos usando el modelo de la IASPEI
HERRMANN - correr el programa de Herrmann para modelar sismogramas sintéticos
HERSEI   -  hacer un archivo en el formato de SEISAN
                    usando el resultado de HERRMANN
WKBJ      -  hacer sismogramas sintéticos con el programa WKBJ
INPUTEPI -  entrar datos del hipocentro y tiempo de origen en el archivo de lectura
INPUTONE - integrar una línea del tipo uno (hipocentro) adicional en el archivo de lectura
RMSDEP     - calcular y plotear RMS en función de la profundidad
UPDATE     - completar el archivo de lectura con hipocentro, magnitud



Comandos en MULPLT

f  (un canal) :   seleciona el canal siguiente
f  (multi canal) :   selecciona el próximo evento
q   - Quit quitar programa
r    - Plot re-plotear evento
z,x,v,b,n,m - filtrar canal(es)
w    - simular sismograma obtenido con sismógrafo Wood-Anderson,
         da movimiento del suelo en Nanometros (nm) en el sismograma siguiente.
         sirve para determinar magnitud ML (Richter original)
o     - otro, seleccionar otros canales
a     - Determinar amplitud. Posicionar cursor al máximo de la oscilación y presionar A.
          Después poner cursor en el minimo de la misma oscilación y prsionar A.
          La amplitud se graba con la próxima fases picada. Se presentan amplitud y período.
          Sirve para determinar magnitudes. Primero filtrar (W para ML, J para Mb)
         después A para determinar amplitud y opicar una fase. Con esta fase se graba la
         amplitud y con la localización se determina la magnitud correspondiente
t     - Toggle Cambiar entre modo un-canal y multi-canal
c     - determinar el final de la coda (para magnitud de caoda, Mc)
d     - Deletear fase. Posicionarcursor cerca de la fase y prosionar D
d     - en el modo multi canal d borra el evento, el archivo del sismograma digital!
s     - Spectrum, Espectro (funciona solo en modo un-canal): computar un espectro.
         Presionar S, seleccionar ventana con el cursor,
         entrar información que pide el programa, sale espectro con el espectro del ruido
         (ruido se toma del sismograma antes de la fase)

S    - lo mismo como s, sin espectro del ruido!
g    - Ground sintetizar el movimiento del suelo,
         solo si existe archivo con característica del sismómetro
h    - Azimut solo en modo un canal: hacer un analisis de 3 componentes
                 (solo si existen 3 comp.).
                  Determina azimut y velocidad aparente. Estos parámetros se
                  integran en el archivo de lectura (S-file) con la próxima fase que se pica.
y     - seleccionar una fase sintética que aparece en pantalla.
                Poner cursor cerca donde podría encontrarse la fase.
                 Presionar y, el programa selecciona la fase más cercana.
e, i   - poner fases E y I
j      - Mb determinar la magnitud Mb (para eventos lejanos).
                Genera un sismograma sintético donde se pueden picar las amplitudes
                (ver arriba, amplitudes) de la fase P!
p     - Registrar Integrar evento en base de datos.
               Programa pide información sobre tipo de evento.
               Se copia archivo de onda al directorio WAV
>     - Print, solo en modo multi canal : imprimir pantalla
*     - Scale escala fija de los sismograms, se selecciona amplitud máxima.
u     - Rotat Rotar componentes. Se usa el azimut determinado por la localización.
            Si se determinó el azimut con H y aplica U,
          el programa pide el azimut al usuario.
-     -   Distance, plotear los canales en orden de la distancia
q    - Quit, Quitar el programa



preparado Wilfried Strauch, wil.gf@ineter.gob.ni
actualizado: 22 diciembre 2000