一种基于改进匹配追踪算法的时频分析方法技术

技术编号:11385371 阅读:66 留言:0更新日期:2015-05-01 12:33
本发明专利技术涉及一种基于匹配追踪算法的时频分析方法,其步骤包括:构建地震信号的复分析信号,确定出到时、延时处的瞬时相位、最佳主频;计算延时位置的主频、相位,作为雷克子波参数,确定原子形态;计算子波能量归一化条件下的子波瞬时振幅;原子振幅可以按给定一个百分比进行匹配,确定出原子到时、原子主频、原子相位、原子振幅这四个关键参数后,从地震道中减去与这个与地震信号最大相关的匹配子波,得到残余信号,对残余信号在一定阈值范围内递推迭代,完成原始信号的分解;将所有匹配子波的时频谱叠加起来得到地震信号的时频谱。本发明专利技术在缩短匹配追踪时频分析的计算效率的同时,保证了原子主频匹配的准确度,解决了地震分频体剖面同相轴不连续的问题。

【技术实现步骤摘要】

本专利技术涉及勘探地球物理领域,具体涉及石油地球物理勘探领域的一种基于改进匹配追踪算法的时频分析方法。 
技术介绍
地震数据匹配追踪时频分析方法在时频原子的匹配上有很多技术方案,最开始技术方案:投影能量最大为原则的自适应搜索方法,这种方式缺点在于,对每一个采样点都要进行多个参数进行扫描,每次迭代过程中都要进行大量的内积运算,计算量过于庞大,计算耗时巨大。如何降解匹配追踪算法的计算复杂度,提高匹配追踪算法计算效率显得尤为关键。其他国内、外的专家学者也在这个方面做出很多尝试,将遗传优化算法引入原子匹配中或是把信号的瞬时特征引入到匹配追踪的算法中,缩短了耗时,但是效率提高的同时也降低了原子匹配的正确性以及地震分频剖面同相轴的连续性。 
技术实现思路
针对现有技术的不足,本专利技术的目的在于:提供一种基于改进匹配追踪算法的时频分析方法。 本专利技术采用的技术方案是:基于改进匹配追踪算法的时频分析方法,其特征在于,该基于改进匹配追踪算法的时频分析方法包括下列步骤: 步骤1:输入一道地震数据; 步骤2:构建地震信号的复分析信号,确定出地震信号的瞬时振幅和瞬时相位,根据地震信号波峰、波谷位置确定出原子到时,并求出此到时处地震信号的瞬时相位; 求取地震信号的瞬时频率,以原子到时处一个波长范围内的所有瞬时频率分别作为原子主频,构建不同的匹配原子,确定出与地震信号最匹配的原子,对应的主频即为到时处的真实原子主频; 步骤3:根据已确定的原子到时、主频、相位,作为雷克子波参数,确定原子波形; 步骤4:计算子波能量归一化条件下的子波瞬时振幅;按用户给定的瞬时振幅百分比参数,确定匹配原子振幅。 步骤5:在原子到时、原子主频、原子相位、原子振幅这四个关键参数确定后,从地震道中减去由这四个参数确定的匹配子波,得到残余信号;若残余信号大于给定阈值,则返回步骤2重新迭代,否则,停止递归,即认为完成原始信号的分解; 步骤6:计算信号二维谱,将所有匹配子波的时频谱叠加起来得到地震信号的时频谱,完成该道地震信号的匹配追踪时频分析; 步骤7:重复上述步骤,完成所有道地震数据的匹配追踪时频分析。 进一步,在步骤2中,所述构建地震信号的复分析信号采用下式来完成 S(t)=s(t)+jHT(s(t))      公式(9)。 所述地震信号的瞬时振幅的确定采用公式(10)来完成; A(t)=|S(t)|=s2(t)+HT2(s(t))]]>公式(10)。 地震信号的瞬时相位的确定采用公式(11)来完成, θ(t)=arctanHT(s(t))s(t)]]>公式(11)。 进一步,在步骤2中,所述匹配是指: 给定原子库D={wk,k=1,2,...K本文档来自技高网...

【技术保护点】
一种基于改进匹配追踪算法的时频分析方法,其特征在于,该基于改进匹配追踪算法的时频分析方法包括下列步骤: 步骤1:输入一道地震数据; 步骤2:构建地震信号的复分析信号,确定出地震信号的瞬时振幅和瞬时相位,根据地震信号波峰、波谷位置确定出原子到时,并求出此到时处地震信号的瞬时相位; 求取地震信号的瞬时频率,以原子到时处一个波长范围内的所有瞬时频率分别作为原子主频,构建不同的匹配原子,确定出与地震信号最匹配的原子,对应的主频即为到时处的真实原子主频; 步骤3:根据已确定的原子到时、主频、相位,作为雷克子波参数,确定原子波形; 步骤4:计算子波能量归一化条件下的子波瞬时振幅;按用户给定的瞬时振幅百分比参数,确定匹配原子振幅; 步骤5:在原子到时、原子主频、原子相位、原子振幅这四个关键参数确定后,从地震道中减去由这四个参数确定的匹配子波,得到残余信号;若残余信号大于给定阈值,则返回步骤2重新迭代,否则,停止递归,即认为完成原始信号的分解; 步骤6:计算信号二维谱,将所有匹配子波的时频谱叠加起来得到地震信号的时频谱,完成该道地震信号的匹配追踪时频分析; 步骤7:重复上述步骤,完成所有道地震数据的匹配追踪时频分析。...

【技术特征摘要】
1.一种基于改进匹配追踪算法的时频分析方法,其特征在于,该基于改进匹配追踪算法的时频分析方法包括下列步骤: 
步骤1:输入一道地震数据; 
步骤2:构建地震信号的复分析信号,确定出地震信号的瞬时振幅和瞬时相位,根据地震信号波峰、波谷位置确定出原子到时,并求出此到时处地震信号的瞬时相位; 
求取地震信号的瞬时频率,以原子到时处一个波长范围内的所有瞬时频率分别作为原子主频,构建不同的匹配原子,确定出与地震信号最匹配的原子,对应的主频即为到时处的真实原子主频; 
步骤3:根据已确定的原子到时、主频、相位,作为雷克子波参数,确定原子波形; 
步骤4:计算子波能量归一化条件下的子波瞬时振幅;按用户给定的瞬时振幅百分比参数,确定匹配原子振幅; 
步骤5:在原子到时、原子主频、原子相位、原子振幅这四个关键参数确定后,从地震道中减去由这四个参数确定的匹配子波,得到残余信号;若残余...

【专利技术属性】
技术研发人员:韩文功冯德永李红梅王兴谋刘力辉宫红波陆蓉王红刘海宁孙兴刚
申请(专利权)人:中国石油化工股份有限公司中国石油化工股份有限公司胜利油田分公司物探研究院
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1