基于时频二次谱提高地震数据分辨率的方法及装置制造方法及图纸

本发明专利技术公开了一种基于时频二次谱提高地震数据分辨率的方法及装置，其涉及地球物理勘探技术领域，所述基于时频二次谱提高地震数据分辨率的方法包括以下步骤：对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱；对所述时频二次谱进行滤波进而获得时变子波振幅谱；基于时变子波振幅谱得到提高分辨率算子；基于所述提高分辨率算子与所述时频谱得到提高分辨率后的输出结果。本发明专利技术中基于时频二次谱提高地震数据分辨率的方法能够有效提高层析反演结果的精度和可靠性,进而服务于薄储层油藏等的后续处理解释工作。

Method and device for improving seismic data resolution based on time frequency two spectrum

The invention discloses a method based on time-frequency two spectral method and apparatus for improving the resolution of seismic data, which relates to the technical field of geophysical exploration, two time-frequency spectrum method to improve the resolution of seismic data includes the following steps: Based on the seismic records of the non steady state of S transform spectrum, then the the spectrum of two-dimensional Fourier transform time-frequency spectrum obtained two times; the two time frequency spectrum filtering and obtain time-varying wavelet amplitude spectrum; based on time-varying wavelet amplitude spectrum and improve the resolution of the operator; improve the resolution operator and the spectrum output results based on the improved resolution. The method for improving the resolution of seismic data based on the time-frequency two spectrum can effectively improve the accuracy and the reliability of the tomographic inversion result, and then serve the subsequent processing and interpretation of the thin oil reservoir.

【技术实现步骤摘要】

本专利技术涉及地球物理勘探
，特别涉及一种基于时频二次谱提高地震数据分辨率的方法及装置。
技术介绍
随着勘探技术的进步，勘探目标也越来越复杂，简单油气藏勘探发展技术已经非常成熟，寻找岩性油气藏、断块油气藏、薄层油气藏等复杂地质目标体已经成为勘探开发的主要研究对象，这就需要尽可能的消除地震子波的影响，提高剖面的分辨率。目前时变子波估计大都建立在稳态子波提取的基础上，在时域上通过小时窗滑动或插值处理来实现。但是，实际地震子波是随时间变化的，地震子波的稳态假设与实际情况并不相符，以此为基础的方法往往不能达到满足分辨率的要求，因此迫切需要与实际情况相符的、在非稳态情况下能够提高分辨率的方法和技术。
技术实现思路
为了克服现有技术的上述缺陷，本专利技术实施例中提供了一种基于时频二次谱提高地震数据分辨率的方法及装置，其能够提高层析反演结果的精度和可靠性。本专利技术实施例的具体技术方案是：一种基于时频二次谱提高地震数据分辨率的方法，其包括以下步骤：对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱；对所述时频二次谱进行滤波进而获得时变子波振幅谱；基于时变子波振幅谱得到提高分辨率算子；基于所述提高分辨率算子与所述时频谱得到提高分辨率后的输出结果。优选地，在所述对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱的步骤中，具体为：对非平稳的地震记录d(t,x)进行S变换得到时频谱A(τ,f,x)，对所述时频谱沿τ和f方向进行二维傅里叶变换得到时频二次谱A(2)(ω,f(2),x)，计算公式如下： A ( 2 ) ( ω , f ( 2 ) , x ) = ∫ - ∞ ∞ ∫ - ∞ ∞ A ( τ , f , x ) e - j 2 π ω τ e - j 2 πff ( 2 ) d τ d f ]]>其中，x表示空间位置，τ表示时间变量，f表示频率变量，ω表示对应于τ的傅里叶变换，f(2)表示对应于f的傅里叶变换，j表示虚数单位。优选地，在所述对所述时频二次谱进行滤波进而获得时变子波振幅谱的步骤中，具体为：将二维低通滤波器H(2)(ω,f(2),x)与所述时频二次谱A(2)(ω,f(2),x)相乘进行二维傅里叶反变换以得到所述时变子波振幅谱，计算公式如下： W ( τ , f , x ) = ∫ - ∞ ∞ ∫ - ∞ ∞ H ( 2 ) ( ω , f ( 2 ) , x ) A ( 2 ) ( ω , f ( 2 ) , x ) e j 2 π ω τ e j 2 πff ( 2 ) dωdf ( 2 ) ]]>其中，f(2)表示对应于f的傅里叶变换，x表示空间位置，ω表示对应于τ的傅里叶变换，τ表示时间变量，j表示，f表示频率变量，W(τ,f,x)表示时变子波振幅谱。优选地，在所述基于时变子波振幅谱得到提高分辨率算子的步骤中，具体为令AMAX(x)＝max{W(τ,f,x)本文档来自技高网...

【技术保护点】
一种基于时频二次谱提高地震数据分辨率的方法，其特征在于，其包括以下步骤：对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱；对所述时频二次谱进行滤波进而获得时变子波振幅谱；基于时变子波振幅谱得到提高分辨率算子；基于所述提高分辨率算子与所述时频谱得到提高分辨率后的输出结果。

【技术特征摘要】
1.一种基于时频二次谱提高地震数据分辨率的方法，其特征在于，其包括以下步骤：对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱；对所述时频二次谱进行滤波进而获得时变子波振幅谱；基于时变子波振幅谱得到提高分辨率算子；基于所述提高分辨率算子与所述时频谱得到提高分辨率后的输出结果。2.根据权利要求1所述的基于时频二次谱提高地震数据分辨率的方法，其特征在于，在所述对非稳态的地震记录进行S变换得到时频谱，再对所述时频谱进行二维傅里叶变换得到时频二次谱的步骤中，具体为：对非平稳的地震记录d(t,x)进行S变换得到时频谱A(τ,f,x)，对所述时频谱沿τ和f方向进行二维傅里叶变换得到时频二次谱A(2)(ω,f(2),x)，计算公式如下： A ( 2 ) ( ω , f ( 2 ) , x ) = ∫ - ∞ ∞ ∫ - ∞ ∞ A ( τ , f , x ) e - j 2 π ω τ e - j 2 πff ( 2 ) d τ d f ]]>其中，x表示空间位置，τ表示时间变量，f表示频率变量，ω表示对应于τ的傅里叶变换，f(2)表示对应于f的傅里叶变换，j表示虚数单位。3.根据权利要求1所述的基于时频二次谱提高地震数据分辨率的方法，其特征在于，在所述对所述时频二次谱进行滤波进而获得时变子波...

【专利技术属性】
技术研发人员：董烈乾，王德营，汪长辉，张慕刚，杨敬磊，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11