[发明专利]基于多子波的地震数据分频方法

权利要求书

1.基于多子波的地震数据分频方法，其特征在于，包括：

S1：生成地震子波序列；

S2：多子波地震道数据分解；

S3：选定参与分频计算的子波或子波范围；

S4：计算选定参与分频计算的子波的频谱和振幅包络；

S5：计算选定参与分频计算的子波的给定分频频率的分频；

S6：对选定参与分频计算的子波的给定分频频率的分频进行叠加处理，所有选定参与分频计算的子波所对应的地震道数据的给定分频频率的分频，等于单个选定参与分频计算的子波的给定分频频率的分频的叠加，即：

其中，f_p为给定分频频率，N是选定参与分频计算的子波个数。

2.根据权利要求1所述的基于多子波的地震数据分频方法，其特征在于，所述步骤S5中，对于选定参与分频计算的子波w_i(t-t_0i)，计算其频谱F_i(f)和振幅包络E_i(t-t_0i)的乘积，得到该子波在时频平面的能量分布，即：

V_i(t,f)＝F_i(f)*E_i(t-t_0i) (3)

对于给定分频频率f_p，依公式(3)计算得到选定参与分频计算子波w_i(t-t₀)的分频为：

V_i(t,f_p)＝F_i(f_p)*E_i(t-t_0i) (4)

其中，*表示两边函数相乘。

3.根据权利要求2所述的基于多子波的地震数据分频方法，其特征在于，所述步骤S3中，选择全部子波参与分频计算，得到相应的原始地震道的分频结果。

4.根据权利要求3所述的基于多子波的地震数据分频方法，其特征在于，所述步骤S3中，选择部分子波参与分频计算，得到相应的部分子波重构的地震道的分频结果。

5.根据权利要求2或4所述的基于多子波的地震数据分频方法，其特征在于，所述步骤S2中，基于生成的地震子波序列，对每个地震道进行分解，得到与每个地震道对应的一组l¹范数优化的，稀疏的，精确重构出原始地震道的地震子波。