4. 微地震数据预处理实战
井中微地震监测使用三分量检波器阵列,事件信号通常是高频、短持续、弱能量的P波和S波初至,而背景中混杂着抽油机机械振动、井筒流体流动、电缆晃动、工频干扰以及随机噪声。原始记录的动态范围可能跨越数个数量级,若直接进入事件检测与定位,弱事件会被噪声淹没,强噪声又可能被误判为事件。因此预处理的目标不是"让剖面好看",而是在保真前提下提高信噪比、恢复相对振幅关系、剔除不可用道、统一时间基准,为后续偏振分析、初至拾取和定位提供可靠输入。
下面以一套典型的井中三分量阵列记录为例,说明预处理流程的组织方式。假设采样率 dt = 0.5 ms(2000 Hz),单级三分量,阵列共 12 级,记录长度 6 s,数据格式为 SEG-Y,道头中已写入检波器深度、方位角与采样间隔。示例中的参数与数值均为教学假设,实际处理需依据现场噪声特征调整。
坏道剔除:先识别,再标记
坏道在井中阵列中很常见,典型表现有三类:一是"死道",整道只有直流偏置或极低幅度的随机噪声;二是"异常强噪道",例如某级检波器耦合不良,持续出现大幅值尖脉冲;三是"极性反转道",三分量接反或电缆接头问题导致同一事件的波形与相邻道反相。识别时不要只看单道波形,应结合阵列统计量:计算每道在事件时窗外的RMS,与全阵列RMS中位数比较,偏离超过3倍即列为可疑;再对相邻道做互相关,若某道与上下相邻道的相关系数显著为负,则判为极性反转。
处理上建议标记而非物理删除:在道头中写入坏道标志位(如SEG-Y的trace_identification_code或自定义字段),后续叠加、偏振分析、定位时按标志位屏蔽。这样既保留了原始数据可追溯性,也避免因误判造成不可逆的信息损失。若某级三分量中仅一个分量损坏,可只屏蔽该分量,另两个分量仍可用于约束。
去噪与滤波:先压制窄带干扰,再做带通
工频干扰(50 Hz及其谐波)在井中记录中往往表现为整道持续的等幅正弦,会严重干扰低频段偏振分析。处理时先用陷波滤波压制50 Hz、100 Hz、150 Hz,陷波器Q值取30~50,带宽不宜过宽,以免在陷波点附近造成明显相位畸变。对于抽油机等周期性机械振动,若其基频稳定,也可用同样的窄带陷波或自适应噪声抵消处理。
随机噪声则用带通滤波配合f-x反褶积。示例中事件主频约150~400 Hz,带通设为80~600 Hz,采用零相位滤波器(正反两次滤波)以避免初至时刻产生时移——这一点对后续走时定位至关重要。滤波后必须做频谱对比:检查事件频带内能量是否保留、初至波形是否被"抹圆"。若发现初至变缓或幅度明显下降,说明高通截止频率设得过高,应回调。
增益恢复:区分"显示增益"与"真振幅补偿"
球面扩散使远离震源的检波器接收到的振幅按传播距离衰减,若不做补偿,深部与浅部道的能量差异会掩盖弱事件。真振幅补偿按几何扩散规律对每道乘以与传播距离相关的补偿因子,补偿后同一事件在不同级的振幅比应接近理论值,这是后续震源机制或能量估算的基础。而AGC(自动增益控制)会强行拉平各时窗幅度,破坏相对振幅关系,只能用于显示,不能用于定量反演。实操中建议保留两套数据:一套真振幅补偿用于分析,一套AGC用于快速浏览。
时间校正:把"钟差"和"延迟"分开处理
时间校正包含两个层面。第一是绝对时间基准:井中记录与地面GPS/井口时钟之间可能存在固定偏移或漂移,需用已知的校准信号(如井口敲击、电缆测试脉冲)对齐,确保事件绝对时刻可用于与压裂施工曲线对照。第二是相对延迟补偿:检波器深度、电缆长度、采样起始延迟都会引入道间时差,需按道头中的深度与采样延迟逐道校正,使同一事件的初至走时只反映传播路径差异。校正后应检查同一事件在相邻道的初至时差是否随深度单调变化,若出现跳变,说明某道时间校正有误。
完成上述五步后,应形成一份可核查的预处理记录:坏道清单及判据、滤波与陷波参数、增益补偿方式、时间校正量。这份记录既是质量控制的依据,也是后续定位结果出现异常时回溯问题的第一手材料。