SeisComP3 文档

scautoloc

定位地震事件

描述

scautoloc 用于准实时自动定位地震事件。通常作为守护进程运行,读取拾取和振幅并实时处理。可离线使用。如果定位符合密集条件,scautoloc 将通报其他应用提供初动。

定位过程

scautoloc 的定位过程包含以下几个步骤:

拾取准备

scautoloc 中每个新进的拾取都需要伴随一个特定的振幅集。因 SeisComP3 数据模块中振幅和拾取都是独立对象,振幅将作为属性添加至相应的拾取。

拾取过滤

过滤每个新进的拾取,检测是否过期,关联振幅是否显示完全。如果台站处理拾取十分频繁,将被认为是假信号,或是递增SNR阈值。

关联

最先尝试将新进拾取与已知初动进行关联。如果关联成功,就可回避集结处理。在确定的环境下,拾取会被关联且不愿成核。

集结

如果直接关联失败,scautoloc 尝试从未关联的预接收的拾取中选择新的初动。这种处理叫作“集结”。scautoloc 在空域和时域上执行栅格检索,需要大量的 CPU 和 RAM 资源。

初动细化

一个通过关联或集结方式生成或更新的初动,仍可能被错误地解释为 P波。scautoloc 试图改良这些基于 SNR 和振幅而拾取的初动。 如果初动包含足够数量的到波,假定能很好定位,scautoloc 则会试图将拾取作为次相位,如 pP。 这种次相位是 弱关联的,不用作定位。

初动过滤

这个过程涉及到最终新建或更新的初动的一致性检查。在此过程中,初动不再被更改。

在集结和关联的过程中,就像初动的细化和过滤,可靠的启发式条件被应用于比较并发初动的 质量。 scautoloc 使用的振幅类型为 snrmb,分别对应于 SNR振幅 和(绝对)mb振幅

栅格文件

栅格配置文件由栅格点(6栏)组成,每个栅格点一行:

-10.00 105.00 20.0 5.0 180.0 8

依次为:栅格点坐标(纬度、经度、深度)、直径、最大台站距离、最小拾取数。上述行设置的栅格点中心位于(南纬 10°,东经 105°),深度为 20km,直径为 5°。最大台站距离为 180°,至少 8 个拾取才能确定初动并定位。

台站配置文件

台站配置文件的每行包含台网代码、台站代码、使用标志(0或1)和最大集结距离。台网和台站代码可以使用通配符*

* * 1 90
GE * 1 180
GE HLG 1 10
TE RGN 0 10

上例,所有台网的所有台站默认创建 90° 范围内的事件。GE 台站可创建任意范围的事件,除了 HLG 台站限制在 10° 范围。TE RGN 台站将被忽略。

日志

scautoloc 产生两种日志文件:包含处理和定位历史的普通日志文件;可选的拾取日志文件。通过关键字 NEWUPDOUT 可查看相关日志:

grep '\(NEW\---UPD\---OUT\)' ~/.seiscomp3/log/scautoloc.log

测试模式

scautoloc 连接消息服务器,接收拾取和振幅,但不发送结果至服务器。日志文件为普通。该模式用于测试新的参数设置。

离线模式

scautoloc 通常作为守护进程在后台运行,不断读取拾取和振幅并实时处理。scautoloc 也可离线操作,用于调试。通过命令行参数--offline设置autoloc.offline来激活。离线模式下,scautoloc 不会连接消息服务器或数据库。实际上,scatuoloc 从拾取文件中读取,例如:

2008-09-25 00:20:16.6 SK LIKS EH __ 4.6 196.953 1.1 A [id]
2008-09-25 00:20:33.5 SJ BEO BH __ 3.0 479.042 0.9 A [id]
2008-09-25 00:21:00.1 CX MNMCX BH __ 21.0 407.358 0.7 A [id]
2008-09-25 00:21:02.7 CX HMBCX BH __ 14.7 495.533 0.5 A [id]
2008-09-24 20:53:59.9 IA KLI BH __ 3.2 143.752 0.6 A [id]
2008-09-25 00:21:04.5 CX PSGCX BH __ 7.1 258.407 0.6 A [id]
2008-09-25 00:21:09.5 CX PB01 BH __ 10.1 139.058 0.6 A [id]
2008-09-25 00:21:24.0 NU ACON SH __ 4.9 152.910 0.6 A [id]
2008-09-25 00:22:09.0 CX PB04 BH __ 9.0 305.960 0.6 A [id]
2008-09-25 00:19:13.1 GE BKNI BH __ 3.3 100.523 0.5 A [id]
2008-09-25 00:23:47.6 RO IAS BH __ 3.1 206.656 0.3 A [id]
2008-09-25 00:09:12.8 GE JAGI BH __ 31.9 1015.304 0.8 A [id]
2008-09-25 00:25:10.7 SJ BEO BH __ 3.4 546.364 1.1 A [id]

这里[id]是占位符,对应实际拾取ID。

Configuration

etc/defaults/global.cfg
etc/defaults/scautoloc.cfg
etc/global.cfg
etc/scautoloc.cfg
~/.seiscomp3/global.cfg
~/.seiscomp3/scautoloc.cfg