4. 微震事件检测与到时拾取

微震监测系统每天会产生大量连续波形，其中真正属于岩体破裂的事件只占很小一部分。要把这些事件从背景噪声中挑出来，并准确标定 P 波、S 波的到达时刻，是后续震源定位和能量计算的前提。到时拾取的误差会直接传递到定位结果：单个台站 5 ms 的拾取偏差，在台网几何不理想时可能造成数十米甚至更大的定位漂移。因此，事件检测与到时拾取既是"第一道工序"，也是精度链条上最敏感的一环。

一、从连续波形到事件候选：检测与拾取的分工

工程上通常把这一环节拆成两个层次。第一层是事件检测，回答"某段时间窗内是否有事件"；第二层是到时拾取，回答"事件波形的初动在哪个采样点"。检测追求召回率，宁可多报也不能漏报；拾取追求精度，需要在检测给出的时间窗内细化到采样点级。两者常串联使用：先用短时窗/长时窗能量比（STA/LTA）在连续记录上扫描出触发段，再在触发段内用 AIC 准则或人工方法确定初动点。

以常见的 500 Hz 采样微震台站为例，一个采样点间隔为 2 ms。若要求 P 波拾取误差控制在 ±4 ms 以内，实际上就是要求拾取结果落在真实初动前后各 2 个采样点之内，这对自动算法和复核流程都提出了明确要求。

二、STA/LTA：最常用的触发检测器

STA/LTA 的核心思想是比较短时窗与长时窗的平均能量（或包络幅值）。事件到来时短时窗能量迅速上升，而长时窗因包含大量背景噪声上升缓慢，二者比值出现明显峰值。其递推形式为：

STA[i] = STA[i-1] + (x[i]^2 - x[i-Ns]^2) / Ns
LTA[i] = LTA[i-1] + (x[i]^2 - x[i-Nl]^2) / Nl
R[i]   = STA[i] / LTA[i]
触发条件：R[i] > R_on 且持续 T_dur；解除条件：R[i] < R_off

参数选择直接决定检测性能。短时窗 Ns 应覆盖 P 波主周期的一部分，长时窗 Nl 通常取短时窗的 5～10 倍。一个可用的起点是：采样率 500 Hz 时取 Ns = 0.1 s（50 点）、Nl = 1.0 s（500 点），触发阈值 R_on = 3.0～4.0，解除阈值 R_off = 1.5，最短持续时间 T_dur = 0.05 s。若矿区机械振动、凿岩或爆破干扰强，可适当提高 R_on 至 5 以上，但会牺牲小震级事件的检出率，需要结合台网噪声水平标定。

STA/LTA 的优点是计算量小、可实时运行、对突变型初动敏感；缺点是对低信噪比事件和渐变型初动不敏感，且阈值对背景噪声变化敏感，单一固定阈值在昼夜噪声差异大的矿区容易误触发或漏触发。实践中常采用分频段计算（如 10～100 Hz 带通后再算比值）或自适应阈值来改善。

三、AIC 准则：在触发段内精确定位初动

赤池信息准则（AIC）把波形看作两段不同统计特性的信号：初动前为噪声段，初动后为信号段。对每个候选点 k 分别拟合两段 AR 模型，取使 AIC 值最小的 k 作为初动点。工程上更常用其简化形式（Maeda 型）：

AIC(k) = k * log(var(x[1:k])) + (N-k-1) * log(var(x[k+1:N]))

其中 N 为分析窗长度，var 为方差。计算时在 STA/LTA 触发点前后各取一段窗口（例如触发点前 0.2 s、后 0.5 s），逐点计算 AIC(k)，取最小值对应的采样点即为 P 波初动。AIC 对信噪比较好的事件效果显著，拾取精度可达 1～2 个采样点；但对低信噪比、初动不清晰的波形，AIC 曲线可能出现多个局部极小，需要结合偏振分析或人工判断。

实际处理中常采用"STA/LTA 粗定位 + AIC 精拾取"的组合：STA/LTA 给出事件大致起点，AIC 在该点邻域内细化。S 波拾取则更困难，可在 P 波拾取后，在 P 波之后的时间窗内对水平分量做同样处理，或利用 S/P 到时差经验范围约束搜索区间。

四、人工复核：自动拾取的必要闭环

自动拾取不可能完全可靠，人工复核是保证到时目录质量的关键环节。复核不是简单"看一眼波形"，而应形成有依据的判断流程。建议按以下顺序检查：

1. 波形形态检查：确认初动方向清晰、P 波起跳点可辨；若初动被噪声淹没，标记为低置信度。
2. 偏振分析辅助：P 波初动在三分量上应表现为质点运动方向与射线方向一致，可用水平/垂直分量能量比或偏振度辅助判断初动点。
3. 多台一致性检查：同一事件在各台站的 P 波拾取应满足走时曲线关系。以初步定位结果为参考，计算各台站拾取残差，残差超过阈值（如 5 ms）的台站重点复核。
4. S 波合理性检查：S 波到时必须晚于 P 波，且 S/P 时差应随震中距增大而增大，与经验关系或理论走时一致。
5. 质量标记：为每个到时赋予质量等级（如 A/B/C），A 级用于精定位，C 级仅作参考或不参与定位。

五、参数标定与效果评估

自动拾取系统上线前，应选取一段包含已知事件（如爆破或已人工确认的微震事件）的数据做标定。以 20～30 个高信噪比事件为参考集，统计 STA/LTA 的检出率、误报率，以及 AIC 拾取与人工拾取的偏差分布。若 P 波拾取偏差的均方根超过 4 ms，或漏检率超过 10%，应调整短时窗长度、触发阈值或前置滤波频带后重新评估。

运行阶段还需持续监控：每日统计触发数量、拾取残差分布和低置信度事件比例。当背景噪声因施工或季节变化而升高时，固定阈值可能失效，此时应引入自适应阈值或按台站分别标定，避免检测性能悄然下降而未被察觉。

综上，STA/LTA 负责"找得到"，AIC 负责"拾得准"，人工复核负责"信得过"。三者构成从连续波形到高质量到时目录的完整链条，其输出直接决定后续震源定位与能量分析的可靠性。掌握各方法的适用条件、关键参数和复核判据，比单纯追求算法复杂度更能提升实际监测效果。