[发明专利]一种用于压制机械因素谐波干扰的方法

权利要求书

1.一种用于压制机械因素谐波干扰的方法，其特征在于，所述方法包括如下步骤：

步骤1：获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线；

步骤2：根据第i炮的地面力信号，计算第i炮中各次谐波分量，并计算得到第i炮的总谐波，其中i为大于等于2的正整数；

所述根据第i炮的地面力信号，计算第i炮中各次谐波分量的具体过程可描述为：

将第i炮的地面力信号转换为频率域表示，得到频率域第i炮的地面力信号g_i(ω)，其中，频率域第i炮的地面力信号g_i(ω)满足：

g_i(ω)满足：

g_i(ω)为频率域第i炮的地面力信号，F_i(ω)为频率域第i炮的基波信号，h_ik(ω)为频率域第i炮的k次谐波，i≧2，k≧2且均为正整数；其中，i为炮的次序，k为第i炮中的总谐波数，m为第i炮中谐波的次序；

将得到的频率域第i炮的k次谐波与频率做互相关，即可得到第i炮的各次谐波分量；

根据第i炮的地面力信号以及从声波测井曲线中提取的滑动扫描第i炮振动记录，计算得到地下反射系数；

步骤3：将步骤2所得第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录做互相关处理，得到谐波畸变信号产生的谐波干扰HARM_i；

步骤4：根据步骤2所得地下反射系数以及步骤3所得谐波畸变信号产生的谐波干扰HARM_i，计算得到第i炮产生的总谐波干扰；

步骤5：从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度t_l以及滑动时间t_h，确定并消除第i炮对其它临近炮的谐波干扰。

2.根据权利要求1所述的一种压制机械因素谐波干扰的方法，其特征在于，所述方法还包括：

步骤6：重复循环步骤2-步骤5，分别逐项确定并消除第i-1炮、第i-2炮、…、第1炮对其它临近炮的谐波干扰。