[发明专利]地震面波自动识别及压制方法

权利要求书

1.地震面波自动识别及压制方法，其特征在于，该地震面波自动识别及压制方法包括：

步骤1，根据原始地震单炮记录，利用修正S变换对不同偏移距的地震单炮记录进行时频分析；

步骤2，在时间—频率—偏移距域谱内进行面波分布时间、分布频率、分布范围、分布能量的同步自动识别；

步骤3，在修正S变换时频谱中提取单频段地震单炮记录；

步骤4，对单频段地震单炮记录进行径向道变换；

步骤5，根据面波的频散特性和能量变化特征，在径向道变换域内对不同的单频段地震单炮记录进行变频率的动态滤波；

步骤6，将所有单频段地震单炮记录中提取的面波噪音进行径向道逆变换。

2.根据权利要求1所述的地震面波自动识别及压制方法，其特征在于，在步骤1中，根据原始地震单炮记录，利用修正S变换对不同偏移距的地震单炮记录进行时频分析，地震单炮记录的时间—频率—偏移距域谱计算公式如下：

其中：s(t,offset)为导入的不同偏移距原始地震单炮记录，t为时间，单位ms，offset为偏移距，单位m；ST(α,fre,offset)是s(t,offset)的时间—频率—偏移距域谱，α为窗函数的时间位置，单位ms，fre为频率，单位Hz；参数β和λ是修正S变换的调节参数。

3.根据权利要求1所述的地震面波自动识别及压制方法，其特征在于，在步骤2中，面波分布时间、分布频率、分布范围、分布能量的四同步自动识别公式如下：

其中：R(t,fre,w,eng)为四同步识别结果，t为时间，单位ms，fre为频率，单位Hz，w为范围，单位m，eng为能量，无量纲单位；当能量差值|ST(t_i,fre,w)|-|ST(t_i-1,fre,w)|满足识别条件时，t₁＝t_i，i为采样点数，t₂为地震单炮记录长度，单位ms；当能量差值|ST(t,fre_j,w)|-|ST(t,fre_j-1,w)|满足识别条件时，fre₁为起始频率值，fre₂＝fre_j，j为时频谱频率采样间隔，单位Hz；当能量差值|ST(t,fre,w_k)|-|ST(t,fre,w_k-1)|满足识别条件时，w₁为起始偏移距，w₂＝w_k，k为偏移距间隔，单位m；eng₁和eng₂为分布能量识别范围，是根据地震单炮记录中面波和有效信号的能量计算得到的。

4.根据权利要求3所述的地震面波自动识别及压制方法，其特征在于，在步骤2中，将地震单炮记录分解成含有面波数据和不含有面波数据，只针对面波分布时间、分布频率、分布范围、分布能量内的数据进行噪音压制，从而保护有效信号不受损失。

5.根据权利要求1所述的地震面波自动识别及压制方法，其特征在于，在步骤3中，充分考虑面波的频散特性和能量变化特征，面波在不同频率段的视速度和能量是不同的，在修正S变换时频谱中提取单频段地震单炮记录，即将地震单炮记录分解成N个单频段地震单炮记录s_j(t,offset)，j＝1,2,…,N，N为面波分布最大频率值，t为时间，单位ms，offset为偏移距，单位m。

6.根据权利要求1所述的地震面波自动识别及压制方法，其特征在于，在步骤4中，对单频段地震单炮记录进行径向道变换，单频段地震单炮记录的径向道变换公式如下：

其中：RT[]为径向道变换，s_j(t,offset)为单频段地震单炮记录，t为时间，单位ms，offset为偏移距，单位m；为径向道变换后数据，t₀为起始时间，单位ms，offset₀为起始偏移距，单位m。