[发明专利]一种地面微地震噪音压制方法

权利要求书

1.一种地面微地震噪音压制方法，其特征在于，包括：

对输入的微地震数据进行时差校正得到校正数据；

对所述校正数据的顺序随机打乱得到待处理数据；

将所述待处理数据进行时窗划分和频带划分得到划分数据；

根据阈值对所述划分数据进行噪音衰减处理得到噪音衰减后的数据；

基于所述随机打乱的方式将所述噪音衰减后的数据的顺序进行恢复处理，得到最终的微地震数据。

2.根据权利要求1所述的地面微地震噪音压制方法，其特征在于，对输入的微地震数据进行时差校正得到校正数据，包括：

基于初始速度和初始震源点位置对输入的微地震数据进行时差校正得到校正数据。

3.根据权利要求1所述的地面微地震噪音压制方法，其特征在于，对所述校正数据的顺序随机打乱得到待处理数据，包括：

利用洗牌算法对所述校正数据的顺序随机打乱得到待处理数据。

4.根据权利要求1所述的地面微地震噪音压制方法，其特征在于，将所述待处理数据进行时窗划分和频带划分得到划分数据，包括：

将所述待处理数据先进行时窗划分得到多个时窗的时域数据；

针对每个时窗，将所述时域数据变换为频域数据；

针对每个时窗，将所述频域数据进行频带划分得到多个频带的划分数据。

5.根据权利要求4所述的地面微地震噪音压制方法，其特征在于，根据阈值对所述划分数据进行噪音衰减处理得到噪音衰减后的数据，包括：

针对每个时窗的每个频带，计算划分数据在每个道的平均能量；

针对每个时窗的每个频带，计算划分数据在每个道的能量阈值；

如果所述道的平均能量大于所述道的能量阈值，则将所述道的相邻至少一道的平均能量的平均值替换所述道的平均能量或将所述道的平均能量替换为零。

6.根据权利要求5所述的地面微地震噪音压制方法，其特征在于，所述平均能量和所述能量阈值根据绝对振幅或均方根振幅计算。

7.根据权利要求4所述的地面微地震噪音压制方法，其特征在于，所述道的能量阈值根据公式T_tfk＝Med_tfk+s_t·M_tfk计算，其中，T_ftk为第t个时窗第k道第f个频带的能量阈值，Med_tfk为第t个时窗第k道第f个频带的中值能量，s_t为系数，M_ftk为第t个时窗第k道第f个频带的中值方差；

根据公式M_tfk＝median_k[|E_tfk-median(E_tfk)|]计算M_ftk，其中，median_k为第t个时窗第k道第f个频带的中值算子，E_ftk为第t个时窗第k道第f个频带的平均能量。

8.根据权利要求4所述的地面微地震噪音压制方法，其特征在于，相邻的两个时窗之间有重合部分。

9.根据权利要求8所述的地面微地震噪音压制方法，其特征在于，相邻的两个时窗的重合部分占两个时窗的比例均小于50％。

10.根据权利要求4所述的地面微地震噪音压制方法，其特征在于，对划分后得到的多个时窗采用Hanning函数进行镶边。