[发明专利]基于预测滤波法和纯相移法的可控震源谐波压制方法

说明书

技术领域

本发明涉及地震资料处理领域，特别是涉及到一种基于预测滤波法和纯相移法的可控震源谐波压制方法。

背景技术

利用传统的谐波压制方法进行谐波压制，效果并不理想，原因在于实际资料处理中谐波干扰包括两部分：一是来自本炮内部自身的干扰，二是来自其它炮(特别是下一炮的初至)的干扰。地面力信号滤波是假定输入的炮记录不存在炮内谐波干扰，所以本炮自身谐波干扰的存在势必影响到方法的应用效果。为此我们发明了一种新的基于预测滤波法和纯相移法的可控震源谐波压制方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种既能较好地压制邻炮干扰噪音，又可部分消除炮内谐波干扰的基于预测滤波法和纯相移法的可控震源谐波压制方法。

本发明的目的可通过如下技术措施来实现：基于预测滤波法和纯相移法的可控震源谐波压制方法，该基于预测滤波法和纯相移法的可控震源谐波压制方法包括：步骤1，输入纯相移处理后的互相关地震记录和相应的地面力信号；步骤2，应用纯相移法对本炮谐波干扰进行压制；步骤3，用预测滤波方法对临炮谐波干扰进行压制：步骤4，获得整个连续的压制谐波干扰后的可控震源扫描数据。

本发明的目的还可通过如下技术措施来实现：

在步骤2中，应用纯相移法分离地面力信号g(t)得到扫描信号s(t)和谐波分量h(t)，并进行傅里叶变换得到频率域的地面力信号G(w)、扫描信号S(w)、谐波分量H(w)，求S(w)的共轭为S^*(w)。

在步骤3中，对于线性升频扫描信号，定义时间域的各阶次谐波信号如下式所示：

s_k(t)＝a_k·sin[2πkf(t)t]

式中：t为扫描时间；k为一正整数，表示谐波阶次；a_k(t)为各次谐波的幅值，f(t)表示扫描频率：

f(t)=fL+fH-fL2Tt]]>

其中：f_L为起始频率，f_H为终止频率，T为信号的扫描长度；

可以得出：谐波信号的频率是基波信号频率的整数倍，谐波畸变信号是基波信号和各次谐波分量的线性叠加，则谐波畸变信号可表达为：

sD(t)=Σk=1Msk(t)=Σk=1Mak(t)sin(2πkft)]]>

式中：s_D(t)表示谐波畸变信号；M表示谐波的最高阶次；

基于褶积模型，即第i炮某道检波器所接收到的连续振动记录是该炮地面力信号与该炮炮点和检波点组成的共中心点处反射系数的褶积：