[发明专利]获取全地层地震分频剖面的方法及装置

说明书

技术领域

本发明涉及地震记录的处理领域，更具体地讲，涉及一种获取全地层地 震分频剖面的方法及装置。

背景技术

分频地震层序分析技术是地震沉积学的一项重要技术，而分频地震层序 分析是在根据地震记录获取的地震分频剖面的基础上进行的。因此，如何获 取具有全地层信息的地震分频剖面日益成为研究的热点问题。

目前，获取地震分频剖面的方法主要有频率滤波法和时频谱分解法。然 而，利用频率滤波法和时频谱分解法获取地震分频剖面时，需要从地震记录 中选择一个频率范围，从而获取所选择频率范围的地震分频剖面。因此，获 取的地震分频剖面中会丢失很多不在所选择频率范围内的同相轴携带的地层 信息。并且，根据获取的地震分频剖面无法得知该频率范围对应的地震分频 剖面与其他频率范围对应的地震分频剖面间的联系。

因此，需要一种获取地震分频剖面的方法，来得到具有全地层信息的地 震分频剖面。

发明内容

本发明的目的在于提供一种获取全地层地震分频剖面的方法及装置，从 而可获得具有全地层信息的地震分频剖面。

本发明的一方面提供一种获取全地层地震分频剖面的方法，包括：(a) 根据地震记录获取地震记录的双谱；(b)根据地震记录的双谱获取地震子波； (c)根据地震子波获取地层的稀疏反射系数；(d)根据地层的稀疏反射系数 和预定主频的俞氏子波获取预定分辨率的全地层地震分频剖面。

可选地，步骤(b)包括：(b1)根据地震记录的双谱的振幅谱获取地震 子波的振幅谱；(b2)根据地震记录的双谱的相位谱获取地震子波的相位谱； (b3)将地震子波的振幅谱和地震子波的相位谱进行反傅里叶变换，得到时 间域的地震子波。

可选地，步骤(b1)包括：由下面的等式获取地震子波的振幅谱： 其中，为由地震子波的振幅谱的自然对数组成的向量，为 由地震记录的双谱的振幅谱的自然对数组成的向量，A为根据地震记录的双 谱的振幅谱的自然对数获得的系数矩阵。

可选地，向量的元素是双谱的两个地震子波的频率在预定频率范围内 的采样在地震记录的双谱中所对应的振幅谱的自然对数。

可选地，所述预定频率范围为[N₁HZ，N₂HZ]，采样间隔为MHZ，N₁为 小于N₂的正数，M为任意正数，所述两个地震子波的频率在所述预定频率范 围内的采样分别被表示为ω₁和ω₂，其中，ω₁与ω₂的部分或所有组合在地震记 录的双谱中所对应的振幅谱的自然对数作为向量的元素。