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

摘要

本发明公开了一种基于双因素的自动识别与去除工业干扰的方法，获取地震资料后,首先对每一道数据做工业干扰频带范围的频率域滤波，分别计算滤波结果的能量A(y(t))和原始数据的能量A(x(t))，并用A(y(t))除以A(x(t))。如果这个比值大于0.001且滤波结果在1s时间内过零点的次数为90～110时，则判定该道存在50Hz工业干扰，否则判定为不存在。对检测出存在干扰的道应用正余弦加权逼近法进行处理，得到去除干扰后的地震数据。本发明是在原始炮数据的时间域对工业干扰进行检测与压制，无需变换到检波点域或频率域，节省了大量的时间和空间；仅对存在干扰的记录做处理，能够有效提高地震资料信噪比，是一种高效率的保真处理技术。

主权项

一种基于双因素的自动识别与去除工业干扰的方法，其特征在于：该自动识别与去除工业干扰的处理方法包括以下步骤: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）将正余弦加权逼近法应用于用上述方法判定的存在干扰的地震记录；判定为不存在的干扰的地震记录则不作任何处理；正余弦加权逼近法的具体计算公式为：w1=c1,w2=c2,u=d;ei=0,zi=0;s1i=cos(2×π×ti×50),s2i=sin(2×π×ti×50);ei=w1×s1i+w2×s2i;zi=xi‑ei;w1=w1+u×zi×s1i,w2=w2+u×zi×s2i;(i=0,KN‑1).式中，xi(i=0,KN‑1)，为原始数据，N为采样点数；w1,w2,u分别为加权值，并给定一个初始值；s1,s2分别为给定的频率为50Hz的初始正余弦函数；ti(i=0,KN‑1)为离散的时间值；ei(i=0,KN‑1)为对干扰的估计；zi(i=0,KN‑1)为去除干扰后的地震记录；7）重复步骤2）到6），直到地震原始记录中所有道的数据都被处理后停止；8）采用常规方法输出压制干扰后的地震资料。