高分辨率时变谱模拟方法、装置、电子设备及介质制造方法及图纸

本申请公开了一种高分辨率时变谱模拟方法、装置、电子设备及介质。该方法可以包括：针对地震数据进行Gabor变换，获得时频振幅谱与时频对数谱；针对时频对数谱进行非线性平滑处理，获得均值滤波后的时频对数谱；获得时频对数谱时间方向上的衰减趋势，计算时变子波时频对数谱；根据时变子波时频对数谱，计算时变子波振幅谱；计算时变的低频端与高频端的频带拓宽量，确定时变子波期望输出；根据时变子波振幅谱与时变子波期望输出，计算高分辨率算子；计算高分辨率时变谱结果。本发明专利技术依据数据信噪比时变特征设计合理的高分辨率算子，实现浅中深层数据分辨率提升合理，参数易于调控的目的，为后续反演、解释等提供高分辨率的数据基础。础。础。

【技术实现步骤摘要】
高分辨率时变谱模拟方法、装置、电子设备及介质


[0001]本专利技术涉及油汽地球物理
，更具体地，涉及一种高分辨率时变谱模拟方法、装置、电子设备及介质。

技术介绍

[0002]地震波在地下介质中传播时，受介质的非完全弹性因素的影响，地震波能量存在吸收衰减，并且高频能量衰减更加严重，导致地震子波是时变的，地震数据的分辨率逐渐降低。实际观测到的地震数据存在噪声，要想在噪声环境中合理的提高地震资料的分辨率，进行相对准确的时变子波估计和合理的提高分辨率算子设计是非常必要的。
[0003]目前研究人员主要从时变子波谱估计和反褶积算子设计两方面来提高地震分辨率：
[0004](1)时变子波振幅谱估计：地震数据是时变的，进行时变振幅谱估计需要将地震数据变换到时频域，将频域谱模拟思想应用到时频域之中，基于Gabor变换、S变换和广义S变换都可以进行时变子波振幅谱估计。
[0005](2)提高分辨率算子设计：常用的时变提高分辨率算子设计一般是在时频域中将瞬时子波的频谱展平，使得瞬时子波的期望输出为单位脉冲函数。<br/>[0006]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种高分辨率时变谱模拟方法，其特征在于，包括：针对地震数据进行Gabor变换，获得时频振幅谱与时频对数谱；针对所述时频对数谱进行非线性平滑处理，获得均值滤波后的时频对数谱；根据所述均值滤波后的时频对数谱，获得时频对数谱时间方向上的衰减趋势，进而计算时变子波时频对数谱；根据所述时变子波时频对数谱，计算时变子波振幅谱；计算时变的低频端与高频端的频带拓宽量，确定时变子波期望输出；根据所述时变子波振幅谱与所述时变子波期望输出，计算高分辨率算子；根据所述时频振幅谱与所述高分辨率算子，计算高分辨率时变谱结果。2.根据权利要求1所述的高分辨率时变谱模拟方法，其中，针对所述时频对数谱进行非线性平滑处理，获得均值滤波后的时频对数谱包括：针对所述时频对数谱的每一个频率进行如下步骤，获得所述均值滤波后的时频对数谱：针对所述时频对数谱开时窗，将所述时窗内数据由小到大排序；删除数据总数中设定比例的极小值和极大值部分，计算剩余部分的均值作为时窗中心位置处的滤波输出；滑动时窗得到当前频率的所述时频对数谱的alpha
‑
trim均值滤波输出。3.根据权利要求1所述的高分辨率时变谱模拟方法，其中，针对所述时频对数谱时间方向上的衰减趋势，沿频率方向通过公式(1)进行多项式拟合，获得所述时变子波时频对数谱：其中，Q
tf
(t
i
,f)为时变子波时频对数谱，i＝1,2,
…
,M，M为每道数据的采样点个数，p
f
为多项式阶数，f为频率，c
k
为拟合系数。4.根据权利要求1所述的高分辨率时变谱模拟方法，其中，通过公式(2)计算所述时变子波振幅谱：w(t,f)＝exp(Q
tf
(t,f))
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中，w(t,f)为时变子波振幅谱，Q
tf
(t
i
,f)为时变子波时频对数谱。5.根据权利要求1所述的高分辨率时变谱模拟方法，其中，计算时变的低频端与高频端的频带拓宽量，确定时变子波期...

【专利技术属性】
技术研发人员：刘炯，刘喜武，霍志周，张丰麒，刘宇巍，张金强，刘倩，
申请(专利权)人：中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：