一种消除伪频谱的频谱分析方法技术

本发明专利技术属于数字信号处理技术领域，尤其涉及地震勘探资料处理中消除伪频谱的分析方法。一种消除伪频谱的频谱分析方法，所述频谱分析方法包括剔除伪频谱步骤，所述伪频谱为数字信号处理中正常频谱以外存在的与正常频谱呈周期关系的频谱；所述剔除伪频谱包括对数字信号进行希尔伯特变换，然后进行傅里叶变换，且将经过希尔伯特变换的结果数据设置为傅里叶变换时输入数据的虚部，经过傅里叶变换后输出频谱，所述伪频谱数字信号被剔除。本发明专利技术输出的频谱图上“伪频谱”消失，极大地改善了信号的频谱分析效果。具有突出的效果。

【技术实现步骤摘要】

本专利技术属于数字信号处理
，尤其涉及的地震勘探资料处理中消除伪频谱的分析方法。
技术介绍
频谱分析是现代信号处理与分析中的一个重要手段，它是将时域信号变换至频域加以分析的方法称为频谱分析，其目的是把复杂的时间历程波形经过傅里叶变换分解为若干单一的谐波分量来研究，从而获得信号的频率结构以及各谐波和相位信息。对信号进行频谱分析可以获得更多有用信息，如求得动态信号中的各个频率成分和频率分布范围，求出各个频率成分的幅值分布和能量分布，从而得到主要幅度和能量分布的频率值。然而在对信号进行傅里叶变换后，往往存在“伪频谱”，类似于数字滤波中的“伪门”，即除了正常频谱外，还存在与正常频谱呈周期关系的频谱，这对信号的频谱分析非常不利，尤其是当真实频谱与“伪频谱”有重叠频率时会导致求得的频谱错误。本专利技术所描述的对常规傅里叶变换进行优化，将其产生的“伪频谱”进行消除，其剔除和消除方法尚未被公开。
技术实现思路
本专利技术介绍了一种消除“伪频谱”的方法，该方法利用优化的傅里叶变换，将常规傅里叶频谱分析中产生的“伪频谱”进行去除。本专利技术的目的是为频谱分析提供了一种消除“伪频谱”的方法。常规傅里叶变换频谱的特点如下，众所周知，傅里叶变换具有周期性的特点，如公式（1）F(m)=ΔtΣn=0N-1f(n)e-i2πmn/N---(1)]]>N是时间域的抽样点数，也是计算出频谱的频率抽样个数，由连续傅里叶变换过渡到离散傅里叶变换时，应用了公式（2）ΔtΔf=1N---(2)]]>公式（2）是完成一对DFT（Discrete Fourier Transform）的条件，否则就不能进行傅里叶正反变换的对应计算。可以计算出N就是傅里叶变换的频率抽样点周期，于是，（1）式可以写为：F(m+N)=ΔtΣn=0N-1f(n)e-i2π(m+N)n/N]]>(3)=ΔtΣn=0N-1f(n)e-i2πmn/Nei2πn]]>由于e-i2πn＝1，故（3）式可写为：F(m+N)=ΔtΣn=0N-1f(n)e-i2πmn/N]]>(4)=F(m)]]>式（4）表示F(m)确实是以N为频率抽样点数的周期，它表示计算F(m)时，如果给定的f(n)是N个值，那么只要计算N个F(m)就行了，在多计算就重复了。如图（1）所示，此时N为44。本专利技术的技术思路及技术实现由于傅里叶变换存在频率抽样点数的周期性，这对样点数较少且采样率较高的数据频谱分析带了干扰，导致数据的真实频谱与它的周期频谱（本专利称为“伪频谱”）交织重叠在一起，目前尚无对策。首先对信号进行希尔伯特变换，然后将结果数据放入傅里叶变换时输入数据的虚部，然后进行傅里叶变换，计算后的频谱图上“伪频谱”消失，极大地改善了信号的频谱分析效果，尤其是频谱与其周期频谱交织重叠在一起时。具体见流程图，如图（2）所示。设实函数道x(t)的希尔伯特变换为h(t)，则希尔伯特变换表达式为：h(t)=1π∫-∞+∞x(τ)t-τdτ---(5)]]>在常规傅里叶变换的程序设置中将输入时的虚部填零，本专利技术创新性地将原数据希尔伯特变换后的数据输入傅里叶变换时的虚部，由于x(t)与h(t)在希尔伯特域正交，故其频率谱中的“伪频谱”逼近于零。具体的，本专利技术的包括，一种消除伪频谱的频谱分析方法，所述频谱分析方法包括剔除伪频谱步骤，所述伪频谱为数字信号处理中正常频谱以外存在的与正常频谱呈周期关系的频谱；所述剔除伪频谱包括对数字信号进行希尔伯特变换，然后进行傅里叶变换，且将经过希尔伯特变换的结果数据设置为傅里叶变换时输入数据的虚部，经过傅里叶变换后输出频谱，所述伪频谱数字信号被剔除。对所述数字信号x(t)进行希尔伯特变换得到h(t)，表达式为：h(t)=1π∫-∞+∞x(τ)t-τdτ---(1)]]>其中，t为时间，τ为时间延迟，+∞为正无穷，-∞为负无穷。所述傅里叶变换采用公式：F(ω)=∫-∞+∞f(t)eiωtdt---(2)]]>其中t为时间，ω为圆频率，i为虚数，+∞为正无穷，-∞为负无穷。采用欧拉公式e±ix＝cosx±isinx简化（2）式，故（2）式可写为：F(ω)=∫-∞+∞f(t)(cos(ωt)-isin(ωt))dt---(3)]]>（3）式虚部为sin(ωt)，令sin(ωt)＝h(t)。本专利技术应用在地震信号频谱分析的方法中，对于含噪的地震数字信号道集进行剔除为频谱的过程包括：步骤（1）输入含噪地震数据信号道集；（2）对所述含噪地震数据信号道集进行希尔伯特变换:含噪地震数字信号x(t)进行希尔伯特变换得到h(t)，其中，t为时间，τ为时间延迟，+∞为正无穷，-∞为负无穷。（3）进行傅里叶变换：完成步骤（2）后，进行傅里叶变换，且设h(t)为傅里叶变换中的虚部；所述傅里叶变换的公式为：所述傅里叶变换采用公式：F(ω)=∫-∞+∞f(t)eiωtdt---(2)]]>其中t为时间，ω为圆频率，i为虚数，+∞为正无穷，-∞为负无穷。采用欧拉公式e±ix＝cosx±isinx简化（2）式，故（2）式可写为：F(ω)=∫-∞+∞f(t)(cos(ωt)-isin(ωt))dt---(3)]]>（3）式虚部为sin(ωt)，令sin(ωt)＝h(t)。（4）输出频谱：将经过傅里叶变换后的含噪地震数据信号频谱图输出；（5）频谱分析：就是傅里叶变换后的数列，然后呈图，是频率-振幅谱，横轴是频率，纵轴是振幅。本专利技术经过对数字信号进行希尔伯特变换和傅里叶变换，且在操作中创造性的将傅里叶变换中的虚部设置为希尔伯特变换后的结果，由于x(t)与h(t)在希尔伯特域正交，故其频率谱中的“伪频谱”逼近于零。使计算后的频谱图上“伪频谱”消失，极大地改善了信号的频谱分析效果。具有突出的效果。附图说明图1为常规傅里叶变换频谱图图本文档来自技高网...

【技术保护点】
一种消除伪频谱的频谱分析方法，其特征在于，所述频谱分析方法包括剔除伪频谱步骤，所述伪频谱为数字信号处理中正常频谱以外存在的与正常频谱呈周期关系的频谱；所述剔除伪频谱包括对数字信号进行希尔伯特变换，然后进行傅里叶变换，且将经过希尔伯特变换的结果数据设置为傅里叶变换时输入数据的虚部，经过傅里叶变换后输出频谱，所述伪频谱数字信号被剔除。

【技术特征摘要】
1.一种消除伪频谱的频谱分析方法，其特征在于，所述频谱分析方法包括剔除伪
频谱步骤，所述伪频谱为数字信号处理中正常频谱以外存在的与正常频谱呈周
期关系的频谱；所述剔除伪频谱包括对数字信号进行希尔伯特变换，然后进行
傅里叶变换，且将经过希尔伯特变换的结果数据设置为傅里叶变换时输入数据
的虚部，经过傅里叶变换后输出频谱，所述伪频谱数字信号被剔除。
2.根据权利要求1所述的一种消除伪频谱的频谱分析方法，其特征在于，对所
述数字信号x(t)进行希尔伯特变换得到h(t)，表达式为：
h(t)=1π∫-∞+∞x(τ)t-τdτ---(1)]]>其中，t为时间，τ为时间延迟，+∞为正无穷，-∞为负无穷；
所述傅里叶变换采用公式：
F(ω)=∫-∞+∞f(t)eiωtdt---(2)]]>其中t为时间，ω为圆频率，i为虚数，+∞为正无穷，-∞为负无穷；
采用欧拉公式e±ix＝cosx±isinx简化（2）式，故（2）式可写为：
F(ω)=∫-∞+∞f(t)...

【专利技术属性】
技术研发人员：刘志成，谢金娥，方伍宝，贾春梅，宋林，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11