一种波形自适应匹配的可控震源谐波消除方法及系统技术方案

本发明专利技术提出了一种波形自适应匹配的可控震源谐波消除方法及系统，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波前后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。本发明专利技术将震源特征信号中的参考扫描信号作为匹配滤波道，设计自适应滤波算子，作用于产生谐波畸变的震源力信号，消除力信号中的谐波干扰。求取原始力信号与滤波后的力信号振幅变化系数，将该系数对滤波后的振幅进行处理，有效保留力信号中有效信号振幅信息，最终达到消除可控震源力信号中谐波干扰的目的。

A Method and System of Vibroseis Harmonic Elimination Based on Waveform Adaptive Matching

The invention provides a method and system for harmonic elimination of vibroseis with adaptive matching of waveforms. The method includes: extracting the reference scanning signal and output signal from the characteristic signal of vibroseis; using the reference scanning signal as the waveform matching channel, adopting the adaptive matching method, matching the force signal with the reference scanning signal channel; and finding out the pre-and post-filtering of the output signal. Amplitude change coefficient, which acts on the output signal after filtering. The reference scanning signal in the source characteristic signal is used as a matched filter channel, and an adaptive filter operator is designed to act on the source force signal which produces harmonic distortion and eliminate the harmonic interference in the force signal. The amplitude variation coefficients of the original force signal and the filtered force signal are obtained. The filtered amplitude is processed by the coefficients. The amplitude information of the effective signal in the force signal is effectively retained, and the harmonic interference in the force signal of the vibroseis can be eliminated finally.

【技术实现步骤摘要】
一种波形自适应匹配的可控震源谐波消除方法及系统
本专利技术属于地震资料处理领域，通过提取可控震源出力信号中有效信号、消除谐波干扰信号，提供一种波形自适应匹配的可控震源谐波消除方法及系统。
技术介绍
在可控震源勘探过程中，由于震源机械装置和震动装置的非线性，以及震板与大地的耦合问题，导致可控震源输出的信号产生畸变。震源平板向地下传输扫描信号的同时，也向地下传输参考信号整数倍的谐波干扰。根据其频率的倍数，将这种谐波干扰依次分为二次谐波、三次谐波及n次谐波。国内外学者提出了许多消除地震数据中高阶谐波，提取有效波信号的方法。Li等(1995)针对线性扫描方式提出了纯相移滤波方法，可有效地压制近场和远场的谐波畸变，它是针对传统的单炮点采集提出的，不适用于可控震源滑动扫描采集方式。针对滑动扫描数据，Fleure(2002)提出了一种在频率域分离谐波并对其压制的技术，Meunier和Bianchi等(2002)提出了一种通过预测谐波，再在每道记录中依次减去谐波的方法。Moerig(2007)提出通过对相关后数据中负时间域的谐波噪声进行统计，来估计地震记录中谐波能量，进而消除谐波畸变的方法。AbdEl-Aal(2010)和ClaudioBagaini(2010)分别在时间域和频率域推导谐波预测算子来预测相关后地震记录中的谐波噪声，但它需要用到相关前数据，无法直接处理相关后数据。CharlesSicking(2009)利用地面力信号和扫描起始时间设计滤波器，将滤波器与基波能量褶积来估计该炮对上一炮记录产生的谐波噪声，该方法可直接用来处理相关后记录。因此，本领域急需一种在有效保留原始地震道振幅变化信息的同时，消除震源信号中的谐波干扰信号的方法。
技术实现思路
可控震源激发过程中，震源平板向地下激发有效信号的同时，有效信号产生畸变生成整数倍频率的谐波干扰。为了解决现有技术中的问题，本专利技术提供一种自适应的波形匹配滤波方法，并对滤波后的振幅进行处理，在有效保留原始地震道振幅变化信息的同时，消除震源信号中的谐波干扰信号。根据本专利技术的一个方面，提供一种波形自适应匹配的可控震源谐波消除方法，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波前后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。进一步地，提取可控震源特征信号包括：提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；分别将参考扫描信号和出力信号的两道地震信号数据提取出来。进一步地，将出力信号与参考扫描信号进行匹配滤波包括：设计自适应滤波器；将自适应滤波器作用于出力信号，得到滤波后的出力信号道。进一步地，设计自适应滤波器包括：采用最小均方LMS算法作为自适应滤波器；对自适应滤波器进行初始化，给定n时刻的抽头输入向量；计算自适应滤波器输出对输入信号的响应；比较输出结果与期望响应产生的估计误差；根据估计误差自动调整滤波器参数。进一步地，求取出力信号滤波后的振幅变化系数包括以下步骤：对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。进一步地，将振幅变化系数作用于滤波后的出力信号道包括：对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；将振幅变化系数和滤波后出力信号振幅谱作乘积，作为新的振幅谱，相位谱保持不变；对新的振幅谱和相位谱作傅里叶逆变换，得到振幅变化作用后的出力信号。根据本专利技术的另一方面，提供一种波形自适应匹配的可控震源谐波消除系统，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波前后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。本专利技术利用震源特征信号中的参考扫描信号，消除出力信号中谐波干扰，从而获得激发的实际震源有效信号。附图说明通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了本专利技术实施例的参考扫描信号的示意图。图2示出了本专利技术实施例的含谐波干扰的出力信号。图3示出了本专利技术实施例的去除谐波的出力信号和参考扫描信号对比图。图4示出了本专利技术实施例的滤波后的出力信号与参考扫描信号的误差值。图5示出了本专利技术实施例的实际出力信号。图6示出了本专利技术实施例的出力信号的时频分析图。图7示出了本专利技术实施例的滤波后的出力信号。图8示出了本专利技术实施例的滤波后出力信号的时频图。图9示出了本专利技术方法的流程图。具体实施方式下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。本专利技术提供了一种基于波形匹配的滤波方法，将震源特征信号中的参考扫描信号作为匹配滤波道，设计自适应滤波算子，作用于产生谐波畸变的震源出力信号，消除出力信号中的谐波干扰。求取原始出力信号与滤波后的出力信号振幅变化系数，将该系数对滤波后的振幅进行处理，有效保留出力信号中有效信号振幅信息，最终达到消除可控震源出力信号中谐波干扰的目的。根据本专利技术的一个实施方式，如图9所示，提供一种波形自适应匹配的可控震源谐波消除方法，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波前后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。优选地，提取可控震源特征信号包括：提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；分别将参考扫描信号和出力信号的两道地震信号数据提取出来。优选地，将出力信号与参考扫描信号进行匹配滤波包括：设计自适应滤波器；具体地，采用最小均方LMS算法作为自适应滤波器，即将均方误差替换为瞬时误差的平方，以瞬时误差信号平方的梯度作为均方误差函数梯度的估计；对自适应滤波器进行初始化，给定n时刻的抽头输入向量；计算自适应滤波器输出对输入信号的响应；比较输出结果与期望响应产生的估计误差；根据估计误差自动调整滤波器参数；最后将自适应滤波器作用于出力信号，得到滤波后的出力信号道。优选地，求取出力信号滤波后的振幅变化系数包括以下步骤：对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。最后，将振幅变化系数作用于滤波后的出力信号道包括：对滤波后本文档来自技高网...

【技术保护点】
1.一种波形自适应匹配的可控震源谐波消除方法，其特征在于，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。

【技术特征摘要】
1.一种波形自适应匹配的可控震源谐波消除方法，其特征在于，该方法包括：提取可控震源特征信号中的参考扫描信号和出力信号；将参考扫描信号作为波形匹配道，采用自适应匹配方法，将出力信号与参考扫描信号道进行匹配滤波；求取出力信号滤波后的振幅变化系数；将振幅变化系数作用于滤波后的出力信号。2.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，提取可控震源特征信号包括：提取可控震源辅助道中参考扫描信号、出力信号、平板加速度信号、重锤加速度信号；分别将参考扫描信号和出力信号的两道地震信号数据提取出来。3.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，将出力信号与参考扫描信号进行匹配滤波包括：设计自适应滤波器；将自适应滤波器作用于出力信号，得到滤波后的出力信号道。4.根据权利要求3所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，设计自适应滤波器包括：采用最小均方LMS算法作为自适应滤波器；对自适应滤波器进行初始化，给定n时刻的抽头输入向量；计算自适应滤波器输出对输入信号的响应；比较输出结果与期望响应产生的估计误差；根据估计误差自动调整滤波器参数。5.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，求取出力信号滤波后的振幅变化系数包括以下步骤：对参考扫描信号作傅里叶变换，获得参考扫描信号的振幅谱；对滤波前出力信号作傅里叶变换，获得出力信号的振幅谱；计算每一个频率点处参考扫描信号振幅值与滤波前出力信号振幅值的比值，将该比值作为出力信号滤波后在此频率点处的振幅变化系数。6.根据权利要求1所述的波形自适应匹配的可控震源谐波消除方法，其特征在于，将振幅变化系数作用于滤波后的出力信号道包括：对滤波后的出力信号作傅里叶变换，获得滤波后的出力信号的振幅谱和相位谱；将振幅变化系数和滤波后出力信号振幅谱作乘积，...

【专利技术属性】
技术研发人员：张洁，殷厚成，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11