[发明专利]一种可控震源地震勘探中消除谐波干扰的方法

摘要

本发明是可控震源地震勘探中消除谐波干扰的方法，首先提取震源力信号，也可以通过重锤和平板的质量以及它们的加速度来求取震源力信号；然后利用震源力信号计算各次谐波分量与基波信号振幅谱的比值，求取预测滤波算子；利用预测滤波算子对地震道记录进行预测滤波处理，估算出每一道的谐波分量；最后从未相关的地震道记录中减去预测出来的谐波分量，达到消除谐波干扰的目的。本发明还可以对相关后道记录进行处理，既可以在时间域进行处理，也可以在频域率进行处理；既可以有选择的处理近道数据，也可以对整个炮记录进行处理，灵活方便，处理效率非常高，能够满足目前野外采集的要求。

主权项

一种可控震源地震勘探中消除谐波干扰的方法，特点是包括如下步骤：1)采集记录可控震源力信号；2)通过震源力信号计算各次谐波信号与基波信号频谱比值，求取预测滤波算子；3)对可控震源地震道数据进行预处理；4)对共炮点道集地震数据第一道做快速傅里叶变换，得到该道的频率域数据；5)利用预测滤波算子对第一道数据进行谐波估测，求取该道的谐波分量： H = R a 1 H 1 · HPO = R a 1 H 1 · Σ m = 2 M ρ m H m H 1 - - - ( 2 ) 式中：H为预测的谐波分量；R为地层的反射系数序列在频率域的表现形式；a1为基波信号的权系数；H1为基波信号的频谱值；HPO为预测滤波算子，用于对道记录数据进行谐波估计；m为谐波阶次；M为不能忽略的谐波的最高阶次；ρm为m次谐波信号与基波信号频谱的比值；Hm为第m次谐波的频谱值；6)从该道的频率域数据中减去预测出来的谐波分量，得到消除谐波后的结果；7)将去掉谐波分量的频率域道数据，进行傅里叶反变换，变到时间域并保存；8)对共炮点道集的其余道依次重复步骤3)至7)，直到完成当前炮的所有道记录；9)对下一个共炮点道集地震数据重复步骤2)至8)，直到完成整个可控震源数据压制谐波干扰的处理。