[发明专利]一种基于双因素的自动识别与去除工业干扰的方法

权利要求书

1.一种基于双因素的自动识别与去除工业干扰的方法，其特征在于：该自动识别与去除工业干扰的处理方法包括以下步骤:

1）获取常规的地震原始记录；

2）选中地震原始记录中的一道数据，对其进行频率域滤波，频带范围选择工业干扰对应的频率范围；

频率域滤波的具体方法是：对选中的一道数据进行傅里叶变换；将傅里叶变换的结果乘以滤波器的频率响应，得到频率域滤波输出；对频率域滤波输出进行傅里叶反变换，得到滤波结果；

其流程为：

滤波器的频率响应的计算公式为：

式中：x(t)是一道地震原始数据，t为时间，单位是ms；X(ω)是x(t)的傅立叶变换；ω为频率，单位是Hz；H(ω)为滤波器的频率响应，Y(ω)为频率域滤波输出，y(t)是滤波后得到的干扰数据；

3）分别计算地震原始数据x(t)和滤波得到的干扰数据y(t)的能量值，并做比值，具体计算如下：

a=A(y(t))A(x(t))]]>

式中：A(x(t))是x(t)的能量，A(y(t))是y(t)的能量值，t为时间，单位是ms；

4）对滤波得到的干扰数据y(t)求其1s时间内过零点的个数b；

5）如果同时满足a≥0.001且90≤b≤110，则判定该道数据中存在工业干扰，否则判定为不存在工业干扰；

6）将正余弦加权逼近法应用于用上述方法判定的存在干扰的地震记录；判定为不存在的干扰的地震记录则不作任何处理；

正余弦加权逼近法的具体计算公式为：

w₁=c₁,w₂=c₂,u=d;e_i=0,z_i=0;

s_1i=cos(2×π×t_i×50),s_2i=sin(2×π×t_i×50);

e_i=w₁×s_1i+w₂×s_2i;z_i=x_i-e_i;

w₁=w₁+u×z_i×s_1i,w₂=w₂+u×z_i×s_2i;

(i=0,KN-1).

式中，x_i(i=0,KN-1)，为原始数据，N为采样点数；w₁,w₂,u分别为加权值，并给定一个初始值；s₁,s₂分别为给定的频率为50Hz的初始正余弦函数；t_i(i=0,KN-1)为离散的时间值；e_i(i=0,KN-1)为对干扰的估计；z_i(i=0,KN-1)为去除干扰后的地震记录；

7）重复步骤2）到6），直到地震原始记录中所有道的数据都被处理后停止；

8）采用常规方法输出压制干扰后的地震资料。

2.根据权利要求1所述的基于双因素的自动识别与去除工业干扰的方法，其特征在于：所述步骤2）中，所述工业干扰对应的频率范围为48～52Hz，进一步优选49～51Hz。

3.根据权利要求1所述的基于双因素的自动识别与去除工业干扰的方法，其特征在于：所述步骤3）中能量的求取方法是通过求时间序列的范数，然后开平方得到的。