SeisComP3 文档

过滤器语法

SeisComP3 支持基于字符串的过滤器定义。

过滤器定义支持内建过滤器链(操作符>>->)以及数学运算符(+-*/^| |)。

示例

A(1,2)>>(B(3,4)*2+C(5,6,7))>>D(8)

这里,A、B、C 和 D 是参数不同的过滤器。 如果过滤一个采样,需经历以下阶段:

s = 采样

  1. 使用 A 过滤 ssa = A(1,2)(s)
  2. 使用 B 过滤 sasb = B(3,4)(sa)
  3. sb = sb * 2
  4. 使用 C 过滤 sasc = C(5,6,7)(sa)
  5. sbsc 相加: sbc = sb + sc
  6. 使用 D 过滤 sbcsf = D(8)(sbc)

sf = 最终采样

scautopick 默认使用的过滤器:

RMHP(10) >> ITAPER(30) >> BW(4,0.7,2) >> STALTA(2,80)

过滤器首先移除偏移。 然后应用 30 秒的 ITAPER。 然后使用 4 序巴特沃斯带通(频率 0.7~2Hz)过滤。 最后使用 2 秒短时窗和 80秒长时窗的 STA/LTA 过滤。

过滤器

下列过滤器函数有效。

AVG(timespan)

计算前面的采样的平均值

  • timespan -- 时间范围

BW(order, lo-freq, hi-freq)

巴特沃斯带通过滤器(BW)

  • order -- 过滤器序号
  • lo-freq -- 低频
  • hi-freq -- 高频

BW_LP(order, hi-freq)

巴特沃斯低通过滤器

  • order -- 过滤器序号
  • hi-freq -- 频率

BW_HP(order, lo-freq)

巴特沃斯高通过滤器

  • order -- 过滤器序号
  • lo-freq -- 频率

BW_HLP(order, lo-freq, hi-freq)

巴特沃斯高低通过滤器

  • order -- 过滤器序号
  • lo-freq -- 低频
  • hi-freq -- 高频

DIFF()

分化过滤器

INT([a = 0])

同化过滤器

权重计算基于参数 a

  • a -- 系数

ITAPER(timespan)

单侧余弦锥

  • timespan -- 时间范围

RMHP(timespan)

高通过滤器

  • timespan -- 时间范围

SM5([type = 1])

5 秒地震计模拟

  • type -- 数据类型
    • 0 -- 偏移
    • 1 -- 速度
    • 2 -- 加速度

STALTA(sta, lta)

短时均值与长时均值比率(STA/LTA)过滤器

  • sta -- 短时窗
  • lta -- 长时窗

WA([type = 1])

伍德安德森地震计模拟

  • type -- 数据类型
    • 0 -- 偏移
    • 1 -- 速度
    • 2 -- 加速度

WWSSN_LP([type = 1])

世界标准地震台网长周地震计模拟

  • type -- 数据类型
    • 0 -- 偏移
    • 1 -- 速度
    • 2 -- 加速度

WWSSN_SP([type = 1])

世界标准地震台网短周地震计模拟

  • type -- 数据类型
    • 0 -- 偏移
    • 1 -- 速度
    • 2 -- 加速度