一种提高地震数据分辨率的方法、装置及存储介质制造方法及图纸

本发明专利技术公开了一种提高地震数据分辨率的方法、装置及存储介质，方法包括：获取当前地震道的时空变子波振幅谱；根据时空变子波振幅谱构建当前地震道对应的频变压缩系数；根据频变压缩系数计算当前地震道的时空变宽频子波振幅谱；基于时空变子波振幅谱和时空变宽频子波振幅谱，计算当前地震道对应的时空变提频算子；利用时空变提频算子对当前地震道进行提频处理。本发明专利技术方案提取的基于时间和空间的频变子波，更加适用于复杂构造工区的需要；基于时空变子波振幅谱和时空变宽频子波振幅谱得到的时空变提频算子可以根据地震数据的局部频谱特征，自适应确定出相对应的时空变提频算子，即使是应用于复杂构造工区的地震数据，也可以有效提高其分辨率。可以有效提高其分辨率。可以有效提高其分辨率。

【技术实现步骤摘要】
一种提高地震数据分辨率的方法、装置及存储介质


[0001]本专利技术涉及油气探勘及数据处理
，尤其涉及一种提高地震数据分辨率的方法、装置及存储介质。

技术介绍

[0002]对于地球物理勘探而言，高分辨率的地震数据有助于分辨更薄的地层，以及更精细地描述小断块、小断层以及油气藏边缘储层的特点。而且，如何提高地震数据的分辨率一直是地球物理勘探的难点。
[0003]现有技术中提出了一种通过设计频变子波的压缩系数以提高地震数据分辨率的方法，通过对叠后地震数据进行提频处理，可以达到较好的提高地震数据分辨率的效果。但是，现有技术中的方案针对全工区提取的是单一子波，根据地震数据的频谱特征可知，当工区的地下构造非常复杂时，地震数据的频谱特征在时空域的差异也较大，也就是说，现有技术中仅提取单一子波的方式无法有效描述复杂构造工区的地震频谱特征，从而造成提频处理后的地震数据在一些区域的垂向分辨率得到有效提升，而另外一些区域并没有得到明显改善。更严重的是，现有技术中仅提取单一子波的方式还会造成地震数据的振幅、能量的时空相对关系发生较大改变，严重影响地震数据提高分辨率后的保幅性。

技术实现思路

[0004]本专利技术要解决的技术问题是：提供一种适用于复杂构造工区的提高地震数据分辨率的方案，使得不仅可以提高地震数据的分辨率，还可以使得提高分辨率后的地震数据具有较好的保幅性。
[0005]为解决上述技术问题，本专利技术提供了一种提高地震数据分辨率的方法，包括：
[0006]获取当前地震道的时空变子波振幅谱；
[0007]根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数；
[0008]根据所述频变压缩系数计算所述当前地震道的时空变宽频子波振幅谱；
[0009]基于所述时空变子波振幅谱和所述时空变宽频子波振幅谱，计算所述当前地震道对应的时空变提频算子；
[0010]利用所述时空变提频算子对所述当前地震道进行提频处理。
[0011]可选地，所述获取当前地震道的时空变子波振幅谱，包括：
[0012]从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道；
[0013]根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱。
[0014]可选地，所述从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道，包括：
[0015]利用空间滑动窗按照设定的搜索半径值，从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道。
[0016]可选地，所述根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变
子波振幅谱，包括：
[0017]利用短时傅里叶变换计算各周边地震道的时频谱，并基于所述时频谱计算各周边地震道的平均时频振幅谱；
[0018]根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱。
[0019]可选地，所述根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱，包括：
[0020]对所述平均时频振幅谱进行时频域对数谱模拟，得到所述当前地震道的时空变子波振幅谱。
[0021]可选地，所述根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数，包括：
[0022]基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数。
[0023]可选地，所述基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数，包括：
[0024]按照如下表达式确定所述当前地震道对应的频变压缩系数a(t,f)，
[0025][0026]其中，f为当前地震道的频率，f
p
(t)为当前地震道第t个采样点的时空变子波振幅谱的峰值频率，f
Nyq
为尼奎斯特频率，f
g
为取值范围为1～2之间的频率常数，a
max
为频变子波最大压缩系数。
[0027]可选地，所述利用所述时空变提频算子对所述当前地震道进行提频处理，包括：
[0028]根据所述时空变提频算子对所述当前地震道进行时频域提频处理，并获得提频处理后的时频谱；
[0029]对所述时频谱在时间方向上进行求和，再对求和结果进行反傅里叶变换，并将反傅里叶变换的实部作为提频处理的结果。
[0030]为解决上述技术问题，本专利技术提供了一种提高地震数据分辨率的装置，包括：
[0031]时空变子波振幅谱获取模块，用于获取当前地震道的时空变子波振幅谱；
[0032]压缩系数构建模块，用于根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数；
[0033]宽频子波振幅谱获取模块，用于根据所述频变压缩系数计算所述当前地震道的时空变宽频子波振幅谱；
[0034]提频算子计算模块，用于基于所述时空变子波振幅谱和所述时空变宽频子波振幅谱，计算所述当前地震道对应的时空变提频算子；
[0035]提频处理模块，用于利用所述时空变提频算子对所述当前地震道进行提频处理。
[0036]可选地，所述时空变子波振幅谱获取模块包括周边地震道确定单元和时空变子波振幅谱获取单元；其中，
[0037]所述周边地震道确定单元，用于从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道；
[0038]所述时空变子波振幅谱获取单元，用于根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱。
[0039]可选地，所述周边地震道确定单元，具体用于利用空间滑动窗按照设定的搜索半径值，从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道。
[0040]可选地，所述时空变子波振幅谱获取单元包括平均时频振幅谱计算子单元和时空变子波振幅谱获取子单元；其中，
[0041]所述平均时频振幅谱计算子单元，利用短时傅里叶变换计算各周边地震道的时频谱，并基于所述时频谱计算各周边地震道的平均时频振幅谱；
[0042]所述时空变子波振幅谱获取子单元，用于根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱。
[0043]可选地，所述时空变子波振幅谱获取子单元，具体用于对所述平均时频振幅谱进行时频域对数谱模拟，得到所述当前地震道的时空变子波振幅谱。
[0044]可选地，所述压缩系数构建模块，用于基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数。
[0045]可选地，所述压缩系数构建模块，具体用于按照如下表达式确定所述当前地震道对应的频变压缩系数a(t,f)，
[0046][0047]其中，f为当前地震道的频率，f
p
(t)为当前地震道第t个采样点的时空变子波振幅谱的峰值频率，f
Nyq
为尼奎斯特频率，f
g
为取值范围为1～2之间的频率常数，a
max
为频变子波最大压缩系数。
[0048]可选地，所述提频处本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种提高地震数据分辨率的方法，其特征在于，包括：获取当前地震道的时空变子波振幅谱；根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数；根据所述频变压缩系数计算所述当前地震道的时空变宽频子波振幅谱；基于所述时空变子波振幅谱和所述时空变宽频子波振幅谱，计算所述当前地震道对应的时空变提频算子；利用所述时空变提频算子对所述当前地震道进行提频处理。2.根据权利要求1所述的提高地震数据分辨率的方法，其特征在于，所述获取当前地震道的时空变子波振幅谱，包括：从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道；根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱。3.根据权利要求2所述的提高地震数据分辨率的方法，其特征在于，所述从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道，包括：利用空间滑动窗按照设定的搜索半径值，从三维目标工区的地震数据中确定出当前地震道对应的各周边地震道。4.根据权利要求2所述的提高地震数据分辨率的方法，其特征在于，所述根据各周边地震道的平均时频振幅谱获取所述当前地震道的时空变子波振幅谱，包括：利用短时傅里叶变换计算各周边地震道的时频谱，并基于所述时频谱计算各周边地震道的平均时频振幅谱；根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱。5.根据权利要求4所述的提高地震数据分辨率的方法，其特征在于，所述根据所述平均振幅谱获取所述当前地震道的时空变子波振幅谱，包括：对所述平均时频振幅谱进行时频域对数谱模拟，得到所述当前地震道的时空变子波振幅谱。6.根据权利要求1所述的提高地震数据分辨率的方法，其特征在于，所述根据所述时空变子波振幅谱构建所述当前地震道对应的频变压缩系数，包括：基于所述时空变子波振幅谱中的峰值频率，确定所述当前地震道对应的频变压缩系数。7.根据权利要求6所述的提高地震数据分辨率的...

【专利技术属性】
技术研发人员：张丰麒，刘俊州，刘兰锋，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：