CN106154318B - 一种地震频谱重建方法 - Google Patents

Abstract

本发明提供了一种地震频谱重建方法，涉及地球物理勘探领域。本方法包括：(1)获取地震信号x(t)，对所述地震信号x(t)进行傅里叶变换，获得该地震信号的振幅谱X(f)；(2)选取子波频谱公式，拾取子波频谱的峰值频率m的初始值；(3)计算子波频谱的峰值振幅a；(4)根据最大匹配投影原理对子波频谱的峰值振幅a和子波频谱的峰值频率m进行局部优化得到局部最优的子波频谱的峰值振幅和子波频谱的峰值频率；(5)将子波频谱的峰值振幅a和子波频谱的峰值频率m带入子波频谱公式获得子波的振幅谱；利用局部最优的子波频谱的峰值振幅a和子波频谱的峰值频率m对地震信号的振幅谱X(f)进行重建，获得重建后的地震频谱以及残差。

Description

一种地震频谱重建方法

技术领域

本发明涉及地球物理勘探领域，具体涉及一种地震频谱重建方法。

背景技术

地震频谱成像技术已成为一项储层特色解释技术，已在地震属性分析方面发挥了重要的作用。频谱分解的核心是信号的时频分析方法，时频分析方法一般分为线性时频分析方法、双线性时频分析方法和参数化时频分析方法等，所有的时频分析方法都有其自身固有的优点和缺陷。利用时频分析方法获得的频率道集，在每一个频率上计算的所谓的“单频信息”都是一定频带范围内信息的叠合，而且反应的是所有地震响应在局部频带范围内信息的叠合；同时，在每一个时间深度上，瞬时谱展示的是整个频率范围内所有的频率成分，没有提供给用户感兴趣的成分，不能获得每一个地震响应所对应的频率成分。

分频属性能够更好地刻画储集层，但是解释人员往往趋向于通过带通滤波拾取一定频带的成分。然而地震子波的主瓣的震荡是有限的，转换到频率域后它对应着一个特定频带范围，解释人员很难从频谱中拾取某一个地震响应的频率范围。如果能够获得每一个频率对应的地震响应，那么就能够对感兴趣的特定储层进行更细致的研究。2010年，Tomasso提出了“频谱重构”的概念，通过正演模拟来估计地震信号的频率成分。不同于传统的频谱分解方法，频谱重构是利用模型的频谱来重建信号的频谱，从地震信号中提取出有意义的分量，而不是分解该信号。频谱重构能够只显示某一个特定的频谱成分，这个特定的成分与某一个特定的层相联系，而这个特定的层可能就是用户感兴趣的目的层。然而，Tomasso采取手工拾取频谱成分和振幅，这种做法不准确，而且严重的依赖于个人的经验，而且针对大数据体应用时也受限制。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种地震频谱重建方法，采用最大匹配投影原理自动对频谱进行重建，重建过程不包含方程求解，稳定性较好，且重建结果不受人为因素影响。

本发明是通过以下技术方案实现的：

一种地震频谱重建方法，包括：

(1)获取地震信号x(t)，对所述地震信号x(t)进行傅里叶变换，获得该地震信号的振幅谱X(f)；

(2)选取子波频谱公式，拾取子波频谱的峰值频率m的初始值；

(3)计算子波频谱的峰值振幅a；

(4)根据最大匹配投影原理对子波频谱的峰值振幅a和子波频谱的峰值频率m进行局部优化得到局部最优的子波频谱的峰值振幅和子波频谱的峰值频率；

(5)将子波频谱的峰值振幅a和子波频谱的峰值频率m带入子波频谱公式获得子波的振幅谱；利用局部最优的子波频谱的峰值振幅a和子波频谱的峰值频率m对地震信号的振幅谱X(f)进行重建，获得重建后的地震频谱以及残差；

(6)判断所述残差是否满足给定条件，如果是，则转入步骤(7)，如果否，则返回步骤(2)；

(7)结束。

所述步骤(2)中子波频谱公式如下：

其中a为该子波频谱的峰值振幅、m为该频谱的峰值频率，f为频率。

所述步骤(3)是利用下式计算子波频谱的峰值振幅a：

其中为子波振幅谱集合，<rⁱX,R_i>表示rⁱX和R_i的内积，为第i次重建振幅谱的标准化；如果i＝1，r^i-1X＝r⁰X即初始振幅谱X；如果i＞1，则r^i-1X为第i-1次重建后的残差；argmax表示取最大投影。

所述步骤(4)是这样实现的：

根据公式(3)估算参数a_i：给定参数a_i的范围，在该范围内取多值作为a_i，将各值与峰值频率m的初始值m_i一起带入公式(1)获得多个子波的振幅谱，这些子波的振幅谱集合即为在该子波振幅谱集合中确定满足公式(3)的一个子波振幅谱R_i(m_i,f)，该振幅谱对应的振幅即为估算的参数a_i；其中参数a_i的范围为n*Peak:Peak，其中n∈[0.1,0.9]，Peak为rⁱX的峰值振幅的值；

在获得初始的估算的参数a_i和m_i后，在局部范围内搜索参数集ξ，局部范围可表示为[ξ-Δξ,ξ+Δξ]，其中Δa为振幅间隔，Δm为峰值频率间隔；在局部范围内，根据公式(3)确定最优的a_i和m_i，获得最优的参数集ξ_i，具体实现方法如下：在局部范围内取多个值作为a_i和m_i，形成多组重建参数ξ_i＝{a_i,m_i}，将每一组参数ξ_i带入公式(1)获得子波的振幅谱，然后寻找一组满足公式(3)参数ξ_i，即为第i次重建的参数ξ_i＝{a_i,m_i}，即得到局部最优的子波频谱的峰值振幅和子波频谱的峰值频率。

所述步骤(5)中的利用局部最优的子波频谱的峰值振幅和子波频谱的峰值频率对地震信号的振幅谱X(f)进行重建，获得残差是这样实现的：

重建后的振幅谱的形式为公式(2)所定义的形式：

其中n为重建次数、a_i为与第i次重建的频谱R_i对应的局部最优的子波频谱的峰值振幅、m_i为与第i次重建的频谱R_i对应的局部最优的子波频谱的峰值频率，rⁿX为第n次重建后的残差，且r⁰X＝X，标记i为1到n的正整数，为第i次重建的子波振幅谱，对于每一次重建，都有

所述步骤(6)中的给定条件为：

所述残差rⁿ⁺¹X的能量应小于等于原始振幅谱能量r⁰X的5％，r⁰X＝X。

与现有技术相比，本发明提供的地震频谱重建方法采用最大匹配投影原理自动对频谱进行重建，子波频谱的振幅和峰值频率可自适应的调节，获得的重建频谱的和更能接近原始频谱；并且重建过程不包含方程求解，稳定性更好；重建过程不受人为因素影响。

附图说明

图1是根据本发明的一个实施方式的步骤框图。

具体实施方式

下面结合附图对本发明作进一步详细描述。

如图1所示，本发明地震频谱重建方法包括以下步骤：

(1)选取子波频谱形式，该子波频谱形式是由公式(1)定义的Ricker子波振幅谱形式：

其中a为该子波频谱的峰值振幅、m为该频谱的峰值频率，f为频率(公式R(m,f)是计算在频率为f时的振幅值)，因此Ricker子波的振幅谱由参数集ξ＝{a,m}(如果选取不同的子波，将采用不同的参数集，Ricker子波只是由{a,m}两个参数确定，如果采用其他的子波，将使用该子波所对应的几个参数来形成参数集即可)来描述；选取子波振幅谱形式可以从本领域相关的基础知识中获取，不限于Ricker子波振幅谱，比如常用的Morlet子波振幅谱也可作为子波振幅谱形式，如果采用其他的地震子波，那么子波的振幅谱形式也会变化，不过计算的整个思路是一样的。对于其他子波振幅谱形式就是其对应的公式。

(2)获取地震信号x(t)，对所述地震信号x(t)进行傅里叶变换，获得该地震信号的振幅谱X(f)；

对地震信号进行傅里叶变换，求取振幅谱都是公知的常用方法，其变换和计算方法在公知教科书中都有详细说明。

(3)对振幅谱X(f)进行多次重建，直到满足重建退出条件(从第一次重建开始，如果第一次重建的结果即满足了退出条件，那么就退出，如果不满足条件，则迭代进行第二次重建，依此类推，直到满足重建条件即停止计算。)，重建后的振幅谱的形式为公式(2)所定义的形式：

其中n为重建次数、a_i为与第i次重建的频谱R_i对应的振幅、m_i为与第i次重建的频谱R_i对应的峰值频率，rⁿX为第n次重建后的残差，且r⁰X＝X，标记i为1到n的正整数，为第i次重建的子波振幅谱，对于每一次重建，都有

根据最大匹配投影原理计算第i次重建时的参数集ξ_i＝{a_i,m_i}；计算第i次重建参数ξ_i＝{a_i,m_i}时包括：拾取振幅谱r^i-1X的峰值振幅对应的频率作为初始m_i，根据公式(3)估算参数a_i：

其中为子波振幅谱集合，<rⁱX,R_i>表示riX和R_i的内积，为第i次重建振幅谱的标准化；如果i＝1，r^i-1X＝r⁰X即初始振幅谱X；如果i＞1，则r^i-1X为第i-1次重建后的残差；argmax表示取最大投影。

具体的，根据公式(3)估算参数a_i包括：给定参数a_i的范围，在该范围内取多值作为a_i，将各值与初始m_i一起带入公式(1)获得多个子波的振幅谱，这些小波的振幅谱集合即为在该子波振幅谱集中确定满足公式(3)的一个子波振幅谱R_i(m_i,f)，该振幅谱对应的振幅即为估算的参数a_i。其中参数a_i的范围为n*Peak:Peak，其中n∈[0.1,0.9]，Peak为rⁱX的峰值振幅的值。

为了更精确的得到第i次的重建参数ξ_i＝{a_i,m_i}，在获得初始的估算的参数a_i和m_i后，还应在局部范围内搜索参数集ξ，局部范围可表示为[ξ-Δξ,ξ+Δξ]，其中Δa为振幅间隔，Δm为峰值频率间隔；在局部范围内，根据公式(3)确定最优的a_i和m_i，获得最优的参数集ξ_i。具体实现方法为：在局部范围内取多个值作为a_i和m_i，形成多组重建参数ξ_i＝{a_i,m_i}，将每一组参数ξ_i带入公式(1)获得子波的振幅谱，然后寻找一组满足公式(3)参数ξ_i，即为第i次重建的参数ξ_i＝{a_i,m_i}(即最优的参数集ξ_i，在范围[ξ-Δξ,ξ+Δξ]内生成多个参数ξi＝{a_i,m_i}，通过公式(3)判断最好的一组作为搜寻结果)。

其所述重建退出条件可以为所述残差rⁿ⁺¹X的能量应小于原始振幅谱能量(即r0X＝X)的5％，甚至更小，具体应视实际情况而定。每一次迭代的信号振幅谱的能量即为按频率对各能量(表达式就是公式(1)，公式(1)表示的是一个频率f下的能量，信号振幅谱的能量包含了多个频率，比如如果信号采样间隔为2ms，那么最大频率就是250HZ，同时如果频率域采样间隔为1HZ，那么就有251个公式(1)，公式(1)中的f就是0-250，将计算得到的251个值相加，再乘以频率域采样间隔，即得到信号振幅谱的能量。)进行积分，即将每个频率下的振幅相加乘以频率域采样间隔即可。

(4)将所述计算出的参数a_i和m_i(就是步骤(3)最后得到的重建的参数ξ_i＝{a_i,m_i})带入公式(1)得到第i次重建出的雷克子波振幅谱。

本发明提供的地震频谱重建方法是多次迭代的过程。对于第一次迭代重建，有即此时将X(f)作为初始振幅谱，根据上述方法确定ξ₁＝{a₁,m₁}，由此得到第一次迭代重建的子波振幅谱R₁(m₁,f)。由于则有在第二次迭代中，将r¹X作为初始振幅谱，根据上述方法确定ξ₂＝{a₂,m₂}，由此获得第二次重建的子波振幅谱R₂(m₂,f)。在第三次迭代中，将r²X作为初始振幅谱。以此类推，在经过N次迭代之后就能得到第N次的残差rⁿX，当残差满足退出条件后，便不再对残差进行重建，即地震频谱重建结束。重建结束后，残差rⁿX应被视为噪音数据对其进行处理，如直接滤掉。

总之，本发明提供的地震频谱重建方法采用最大匹配投影原理自动对频谱进行重建，子波频谱的振幅和峰值频率可自适应的调节，获得的重建频谱的和更能接近原始频谱；并且重建过程不包含方程求解，稳定性更好；重建过程不受人为因素影响，可方便的应用于较大的地震数据，且计算效率较高。

上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。

Claims (3)

1.一种地震频谱重建方法，其特征在于：所述方法包含以下步骤：

(1)获取地震信号x(t)，对所述地震信号x(t)进行傅里叶变换，获得该地震信号的振幅谱X(f)；

(2)选取子波频谱公式，拾取子波频谱的峰值频率m的初始值；

(3)计算子波频谱的峰值振幅a；

(4)根据最大匹配投影原理对子波频谱的峰值振幅a和子波频谱的峰值频率m进行局部优化得到局部最优的子波频谱的峰值振幅和子波频谱的峰值频率；

(5)将子波频谱的峰值振幅a和子波频谱的峰值频率m带入子波频谱公式获得子波的振幅谱；利用局部最优的子波频谱的峰值振幅a和子波频谱的峰值频率m对地震信号的振幅谱X(f)进行重建，获得重建后的地震频谱以及残差；

(6)判断所述残差是否满足给定条件，如果是，则转入步骤(7)，如果否，则返回步骤(2)；

(7)结束,

所述步骤(2)中子波频谱公式如下：

其中a为该子波频谱的峰值振幅、m为该频谱的峰值频率，f为频率；

所述步骤(3)是利用下式计算子波频谱的峰值振幅a：

其中为子波振幅谱集合，<rⁱX,R_i>表示rⁱX和R_i的内积，为第i次重建振幅谱的标准化；如果i＝1，r^i-1X＝r⁰X即初始振幅谱X；如果i＞1，则r^i-1X为第i-1次重建后的残差；argmax表示取最大投影，

所述步骤(4)是这样实现的：

根据公式(3)估算参数a_i,a_i为与第i次重建的频谱R_i对应的局部最优的子波频谱的峰值振幅,标记i为1到n的正整数：给定参数a_i的范围，在该范围内取多值作为a_i，将各值与峰值频率m的初始值m_i一起带入公式(1)获得多个子波的振幅谱，这些子波的振幅谱集合即为在该子波振幅谱集合中确定满足公式(3)的一个子波振幅谱R_i(m_i,f)，该振幅谱对应的振幅即为估算的参数a_i；其中参数a_i的范围为j*Peak～Peak，其中j∈[0.1,0.9]，Peak为rⁱX的峰值振幅的值；

在获得初始的估算的参数a_i和m_i后，在局部范围内搜索参数集ξ，局部范围可表示为[ξ-Δξ,ξ+Δξ]，其中Δa为振幅间隔，Δm为峰值频率间隔；在局部范围内，根据公式(3)确定最优的a_i和m_i，获得最优的参数集ξ_i，具体实现方法如下：在局部范围内取多个值作为a_i和m_i，形成多组重建参数ξ_i＝{a_i,m_i}，将每一组参数ξ_i带入公式(1)获得子波的振幅谱，然后寻找一组满足公式(3)参数ξ_i，即为第i次重建的参数ξ_i＝{a_i,m_i}，即得到局部最优的子波频谱的峰值振幅和子波频谱的峰值频率。

2.根据权利要求1所述的地震频谱重建方法，其特征在于：所述步骤(5)中的利用局部最优的子波频谱的峰值振幅和子波频谱的峰值频率对地震信号的振幅谱X(f)进行重建，获得残差是这样实现的：

重建后的振幅谱的形式为公式(2)所定义的形式：

其中n为重建次数、a_i为与第i次重建的频谱R_i对应的局部最优的子波频谱的峰值振幅、m_i为与第i次重建的频谱R_i对应的局部最优的子波频谱的峰值频率，rⁿX为第n次重建后的残差，且r⁰X＝X，标记i为1到n的正整数，为第i次重建的子波振幅谱，对于每一次重建，都有

3.根据权利要求2所述的地震频谱重建方法，其特征在于：所述步骤(6)中的给定条件为：

所述残差rⁿ⁺¹X的能量应小于等于原始振幅谱能量r⁰X的5％，r⁰X＝X。