[发明专利]地震资料处理方法及其阙值和有效信号频率的求取方法

说明书

本发明涉及一种地震资料处理方法及其阙值和有效信号频率的求取方法。所述阙值求取方法包括获取地震资料预定时窗内的频谱曲线；计算频谱曲线上不同频率的叠加频谱，然后对叠加频谱进行求导，得到叠加频谱变化率曲线；选取叠加频谱变化率曲线上叠加频谱变化率最大值所对应的频率作为阈值。所述有效信号频率求取方法包括根据所述阙值确定出有效频带范围。所述地震资料处理方法采用所述地震资料处理中有效信号频率的求取方法求取有效信号频率。本发明能够定量的得到阙值大小，使得求取的阙值更加准确有效；并且能够有效地将噪声与有效信号区分开来，求取的有效频率更为准确、可靠，从而保证了地震资料处理与解释的准确性。

技术领域

本发明属于石油地震勘探地震资料处理领域，具体地说，涉及一种地震资料处理中阙值的求取方法、地震资料处理中有效信号频率的求取方法以及地震资料处理方法。

背景技术

当前，在石油勘探领域，噪声是影响地震资料品质的一个重要因素，噪声越多、越强，对有效信号的影响越大，叠加产生的地质剖面也就越差。在整个地震资料处理的过程中，如果能有效地压制噪声，保护有效信号，处理结果就能较为真实的还原地质形态面貌，这对于地震解释分析是有利的。但是由于野外采集到的地震资料中，噪声通常会与有效信号混杂在一起，因此需要选取一个阙值将噪声与有效信号区分开来，这对于地震处理的过程质量控制也是十分有必要的。

而在阙值的求取过程中，通常是根据经验人为地进行噪声与有效信号的划分，从而确定阙值大小。这种判定方法很模糊、很主观，并且阙值的取值缺乏定量的判定和严密的计算，因此，在这种情况下，处理效果往往会大打折扣。

因此，需要一种地震资料处理中阙值和有效频率的求取方法来对地质勘探过程中产生的地震数据进行有效地处理和准确的解释。

发明内容

本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的之一是解决在地震资料处理过程中确定阙值大小时的随机性与不准确性，缺乏定量判定和严密计算的技术问题。

总体来讲，本发明通过对时窗内数据进行傅里叶变换得到频谱曲线，频谱曲线进行叠加求和后，进而计算叠加频谱曲线的变化率的方式求取阙值。

为了实现上述目的，本发明一方面提供了一种地震资料处理中阙值的求取方法，所述方法可以包括以下步骤：获取地震资料预定时窗内的频谱曲线；计算所述频谱曲线上不同频率的叠加频谱，然后对所述叠加频谱进行求导，得到叠加频谱变化率曲线；选取所述叠加频谱变化率曲线上叠加频谱变化率最大值所对应的频率作为阈值。

在本发明的地震资料处理中阙值的求取方法的一个示例性实施例中，所述获取地震资料预定时窗内的频谱曲线的步骤可以包括：在地震资料上选取一个时窗；对所选时窗内的地震资料进行傅里叶变换，以将所选时窗内的地震资料由时间域变换到频率域，得到频谱曲线。其中，可以通过等式(1)对所选时窗内的地震资料进行傅里叶变换：

在等式(1)中，x(t)为所选时窗的地震资料，t表示时间，单位为秒，w表示频率，单位为Hz，

在本发明的地震资料处理中阙值的求取方法的一个示例性实施例中，可以通过等式(2)计算所述频谱曲线上不同频率的叠加频谱：

E_w＝∑X(w)w＝w₁,w₂,…,w_max 等式(2)

在等式(2)中，E_w为频谱曲线上不同频率的叠加频谱，w表示频率，w₁,w₂,…,w_max表示频谱曲线上的不同频率，w_max表示频谱曲线上频率的最大值。

在本发明的地震资料处理中阙值的求取方法的一个示例性实施例中，可以通过等式(3)对所述叠加频谱进行求导，得到叠加频谱变化率曲线：