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

摘要

本发明涉及地震资料保幅去噪技术领域，特别是涉及一种用于压制机械因素谐波干扰的方法，该压制方法可靠性高、易于实现，可以有效的压制机械因素产生的谐波干扰，且不会对与谐波干扰频率相近的有效信号产生损害。该压制方法包括如下步骤：获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线；计算第i炮中各次谐波分量、第i炮的总谐波、地下反射系数；计算谐波畸变信号产生的谐波干扰；计算得到第i炮产生的总谐波干扰；确定并消除第i炮对其它临近炮的谐波干扰。

主权项

1.一种用于压制机械因素谐波干扰的方法，其特征在于，所述方法包括如下步骤：步骤1：获取可控震源原始单炮记录、单炮地面力信号以及声波测井曲线；步骤2：根据第i炮的地面力信号，计算第i炮中各次谐波分量，并计算得到第i炮的总谐波，其中i为大于等于2的正整数；所述根据第i炮的地面力信号，计算第i炮中各次谐波分量的具体过程可描述为：将第i炮的地面力信号转换为频率域表示，得到频率域第i炮的地面力信号gi(ω)，其中，频率域第i炮的地面力信号gi(ω)满足：g_i(ω)满足：gi(ω)为频率域第i炮的地面力信号，Fi(ω)为频率域第i炮的基波信号，hik(ω)为频率域第i炮的k次谐波，i≧2，k≧2且均为正整数；其中，i为炮的次序，k为第i炮中的总谐波数，m为第i炮中谐波的次序；将得到的频率域第i炮的k次谐波与频率做互相关，即可得到第i炮的各次谐波分量；根据第i炮的地面力信号以及从声波测井曲线中提取的滑动扫描第i炮振动记录，计算得到地下反射系数；步骤3：将步骤2所得第i炮的总谐波与从声波测井曲线中提取的滑动扫描第i炮振动记录做互相关处理，得到谐波畸变信号产生的谐波干扰HARMi；步骤4：根据步骤2所得地下反射系数以及步骤3所得谐波畸变信号产生的谐波干扰HARMi，计算得到第i炮产生的总谐波干扰；步骤5：从声波测井曲线中提取的滑动扫描第i炮振动记录中确定出扫描信号长度T、记录长度tl以及滑动时间th，确定并消除第i炮对其它临近炮的谐波干扰。