一种多分量地震资料面波压制方法

摘要

本发明公开了一种多分量地震资料面波压制方法，该方法包括：获取在多个时间采样点采样得到的单个多分量矢量地震道信号Sk(t)；对Sk(t)进行时频双域分解，以得到时频谱信号WgSk(t，a)；对时频谱信号WgSk(t，a)进行瞬时极化分析，其中该瞬时极化分析包括计算出针对各个时间采样点的瞬时极化椭圆率P(t，a)；基于视速度大小，将WgSk(t，a)分为非面波区信号和面波区信号，并对面波区信号进行面波信号压制处理，以得到压制面波信号后的时频谱信号W′gSk(t，a)；对W′gSk(t，a)进行反变换，并重构以得到压制面波后的矢量地震道信号S′k(t)；其中，面波信号压制处理包括：统计出在各个频率下的非面波区信号的均方根能量EnonGR(a)；针对每个频率下的面波区信号，以该频率下的均方根能量EnonGR(a)为能量阈值将面波区信号分为强能量信号和弱能量信号后进行分类约束极化滤波。

主权项

一种多分量地震资料面波压制方法，该方法包括：获取在多个时间采样点采样得到的单个多分量矢量地震道信号Sk(t)，k为正整数；对所述Sk(t)进行时频双域分解，以得到时频谱信号WgSk(t，a)，t表示时间，a表示频率；对所述时频谱信号WgSk(t，a)进行瞬时极化分析，其中该瞬时极化分析包括计算出针对各个时间采样点的瞬时极化椭圆率P(t，a)；基于视速度大小，将所述WgSk(t，a)分为非面波区信号和面波区信号并对所述面波区信号进行面波信号压制处理，以得到压制面波信号后的时频谱信号W′gSk(t，a)；对所述W′gSk(t，a)进行反变换并重构以得到压制面波后的矢量地震道信号S′k(t)；其中，所述面波信号压制处理包括：统计出在各个频率下的所述非面波区信号的均方根能量EnonGR(a)；针对每个频率下的所述面波区信号，以该频率下的所述均方根能量EnonGR(a)为能量阈值将所述面波区信号分为强能量信号和弱能量信号后进行分类约束极化滤波。