[发明专利]一种微震监测数据谐波干扰压制方法及装置

权利要求书

1.一种微震监测数据谐波干扰压制方法，其特征在于，该方法包括：

采集微地震数据，将每个检波点记录的地震数据分选在一起作为一个共检波点道集；

将所述共检波点道集中的每一记录道作FFT变换，获取频域地震数据X(x,f)；其中，x表示采样点所在记录道的道号，f表示采样点对应的频率；

利用所述频域地震数据X(x,f)获得振幅谱A(x,f)；

对所述振幅谱A(x,f)的每个频率采样进行多道中值滤波，获得第一次中值滤波后的振幅谱A₁(x,f)；

对第一次中值滤波后的振幅谱A₁(x,f)的每一记录道进行多点中值滤波，获得第二次中值滤波后的振幅谱A₂(x,f)；

利用所述第一次中值滤波后的振幅谱A₁(x,f)和所述第二次中值滤波后的振幅谱A₂(x,f)获得滤波因子；

利用所述滤波因子对所述频域地震数据进行滤波处理，并将滤波处理后的结果进行逆FFT变换，得到谐波干扰压制后的时间域微地震数据。

2.如权利要求1所述的方法，其特征在于，所述滤波因子为：

F(x,f)＝A₂(x,f)/A₁(x,f)

其中，F(x,f)为滤波因子。

3.如权利要求1所述的方法，其特征在于，所述振幅谱为：

A(x,f)＝sqrt(X_r²(x,f)+X_i²(x,f))

其中，X_r(x,f)、X_i(x,f)分别为频域地震数据X(x,f)的实部和虚部。

4.如权利要求1所述的方法，其特征在于，所述采集微地震数据的步骤包括：

通过地面多条排列或井中多级检波器采集单分量或三分量微地震数据；

将连续记录的微地震数据按固定长度分割为多个文件存储。

5.一种微震监测数据谐波干扰压制装置，其特征在于，该装置包括：

共检波点道集获取单元，用于采集微地震数据，将每个检波点记录的地震数据分选在一起作为一个共检波点道集；

频域地震数据获取单元，用于将所述共检波点道集中的每一记录道作FFT变换，获取频域地震数据X(x,f)；其中，x表示采样点所在记录道的道号，f表示采样点对应的频率；

振幅谱获取单元，用于利用所述频域地震数据X(x,f)获得振幅谱A(x,f)；

第一中值滤波单元，用于对所述振幅谱A(x,f)的每个频率采样进行多道中值滤波，获得第一次中值滤波后的振幅谱A₁(x,f)；

第二中值滤波单元，用于对第一次中值滤波后的振幅谱A₁(x,f)的每一记录道进行多点中值滤波，获得第二次中值滤波后的振幅谱A₂(x,f)；

滤波因子获取单元，用于利用所述第一次中值滤波后的振幅谱A₁(x,f)和所述第二次中值滤波后的振幅谱A₂(x,f)获得滤波因子；

时间域微地震数据获取单元，用于利用所述滤波因子对所述频域地震数据进行滤波处理，并将滤波处理后的结果进行逆FFT变换，得到谐波干扰压制后的时间域微地震数据。

6.如权利要求5所述的装置，其特征在于，所述滤波因子获取单元得到的滤波因子为：

F(x,f)＝A₂(x,f)/A₁(x,f)

其中，F(x,f)为滤波因子。

7.如权利要求5所述的装置，其特征在于，所述振幅谱获取单元得到的振幅谱为：

A(x,f)＝sqrt(X_r²(x,f)+X_i²(x,f))

其中，X_r(x,f)、X_i(x,f)分别为频域地震数据X(x,f)的实部和虚部。

8.如权利要求5所述的装置，其特征在于，所述共检波点道集获取单元包括：

采集模块，用于通过地面多条排列或井中多级检波器采集单分量或三分量微地震数据；

存储模块，用于将连续记录的微地震数据按固定长度分割为多个文件存储。