[发明专利]一种提高地震数据分辨率的方法、装置及存储介质

权利要求书

1.一种提高地震数据分辨率的方法，其特征在于，包括：

获取当前地震道的时空变子波振幅谱；

根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数；

根据所述频变压缩系数计算所述当前地震道的时空变宽频子波振幅谱；

基于所述时空变子波振幅谱和所述时空变宽频子波振幅谱，计算所述当前地震道对应的时空变提频算子；

利用所述时空变提频算子对所述当前地震道进行提频处理。

2.根据权利要求1所述的提高地震数据分辨率的方法，其特征在于，所述获取当前地震道的时空变子波振幅谱，包括：

从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道；

根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱。

3.根据权利要求2所述的提高地震数据分辨率的方法，其特征在于，所述从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道，包括：

利用空间滑动窗按照设定的搜索半径值，从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道。

4.根据权利要求2所述的提高地震数据分辨率的方法，其特征在于，所述根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱，包括：

利用短时傅里叶变换计算各周边地震道的时频谱，并基于所述时频谱计算各周边地震道的平均时频振幅谱；

根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱。

5.根据权利要求4所述的提高地震数据分辨率的方法，其特征在于，所述根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱，包括：

对所述平均时频振幅谱进行时频域对数谱模拟，得到所述当前地震道的时空变子波振幅谱。

6.根据权利要求1所述的提高地震数据分辨率的方法，其特征在于，所述根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数，包括：

基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数。

7.根据权利要求6所述的提高地震数据分辨率的方法，其特征在于，所述基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数，包括：

按照如下表达式确定所述当前地震道对应的频变压缩系数a(t,f)，

其中，f为当前地震道的频率，f_p(t)为当前地震道第t个采样点的时空变子波振幅谱的峰值频率，f_Nyq为尼奎斯特频率，f_g为取值范围为1～2之间的频率常数，a_max为频变子波最大压缩系数。

8.根据权利要求1所述的提高地震数据分辨率的方法，其特征在于，所述利用所述时空变提频算子对所述当前地震道进行提频处理，包括：

根据所述时空变提频算子对所述当前地震道进行时频域提频处理，并获得提频处理后的时频谱；

对所述时频谱在时间方向上进行求和，再对求和结果进行反傅里叶变换，并将反傅里叶变换的实部作为提频处理的结果。

9.一种提高地震数据分辨率的装置，其特征在于，包括：

时空变子波振幅谱获取模块，用于获取当前地震道的时空变子波振幅谱；

压缩系数构建模块，用于根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数；

宽频子波振幅谱获取模块，用于根据所述频变压缩系数计算所述当前地震道的时空变宽频子波振幅谱；

提频算子计算模块，用于基于所述时空变子波振幅谱和所述时空变宽频子波振幅谱，计算所述当前地震道对应的时空变提频算子；

提频处理模块，用于利用所述时空变提频算子对所述当前地震道进行提频处理。