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

权利要求书

1.一种地震资料处理中阙值的求取方法，其特征在于，所述求取方法包括以下步骤：

获取地震资料预定时窗内的频谱曲线；

通过等式(2)计算所述频谱曲线上不同频率的叠加频谱，然后通过等式(3)对所述叠加频谱进行求导，得到叠加频谱变化率曲线；

选取所述叠加频谱变化率曲线上叠加频谱变化率最大值所对应的频率作为阙值，

其中，

所述等式(2)为：E_w＝∑X(w) w＝w₁,w₂,…,w_max，

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

所述等式(3)为：

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

2.根据权利要求1所述的地震资料处理中阙值的求取方法，其特征在于，所述获取地震资料预定时窗内的频谱曲线的步骤包括：

在地震资料上选取一个时窗；

对所选时窗内的地震资料进行傅里叶变换，以将所选时窗内的地震资料由时间域变换到频率域，得到频谱曲线。

3.根据权利要求2所述的地震资料处理中阙值的求取方法，其特征在于，通过等式(1)对所选时窗内的地震资料进行傅里叶变换，

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

4.一种地震资料处理中有效信号频率的求取方法，其特征在于，所述求取方法包括以下步骤：

获取地震资料预定时窗内的频谱曲线；

将所述频谱曲线上最大振幅对应的频率确定为主频，并根据所述主频将所述频谱曲线划分为第一段频谱曲线和第二段频谱曲线，其中，第一段频谱曲线的频率小于所述主频，第二段频谱曲线的频率大于所述主频；

通过等式(4)计算所述第一段频谱曲线上不同频率的叠加频谱E′_w，然后通过等式(6)对所述第一段频谱曲线上不同频率的叠加频谱进行求导，得到第一段叠加频谱变化率曲线；

选取所述第一段叠加频谱变化率曲线上叠加频谱变化率最大值所对应的频率作为第一阙值；

通过等式(5)计算所述第二段频谱曲线上不同频率的叠加频谱E″_w，然后通过等式(7)对所述第二段频谱曲线上不同频率的叠加频谱进行求导，得到第二段叠加频谱变化率曲线；

选取所述第二段叠加频谱变化率曲线上叠加频谱变化率最大值所对应的频率作为第二阙值；

将所述第一阙值和所述第二阙值之间的频率确定为有效信号频率，

其中，

等式(4)为E′_w＝∑X(w) w＝w′₁,w′₂,…,w′_max，

在等式(4)中，E′_w为所述第一段频谱曲线上不同频率的叠加频谱，w表示频率，w′₁,w′₂,…,w′_max为所述第一段频谱曲线上的频率，w′_max表示所述第一段频谱曲线上频率的最大值，

所述等式(5)为E″_w＝∑X(w) w＝w″₁,w″₂,…,w″_max，

在等式(5)中，E″_w为所述第二段频谱曲线上不同频率的叠加频谱，w表示频率，w″₁,w″₂,…,w″_max为所述第二段频谱曲线上的频率，w″_max表示所述第二段频谱曲线上频率的最大值，

所述等式(6)为

在等式(6)中，K′_w为所述第一段频谱曲线上不同频率的叠加频谱的一次偏导数，E′_w为所述第一段频谱曲线上不同频率的叠加频谱，w表示频率，w′₁,w′₂,…,w′_max为所述第一段频谱曲线上的频率，w′_max表示所述第一段频谱曲线上频率的最大值，

所述等式(7)为

在等式(7)中，K″_w为所述第二段频谱曲线上不同频率的叠加频谱的一次偏导数，E″_w为所述第二段频谱曲线上不同频率的叠加频谱，w表示频率，w″₁,w″₂,…,w″_max为所述第二段频谱曲线上的频率，w″_max表示所述第二段频谱曲线上频率的最大值。