一种用谱反演提高地震分辨率方法技术

本发明专利技术公开了一种用谱反演提高地震分辨率方法，包括：获取地震资料并载入地震数据；从地震数据中提取出地震小波数据，并变换地震小波数据到频域W(f)，将地震数据变换到频域S(f)；通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数；通过谱分解对目标函数进行分析。本发明专利技术通过利用谱反演技术能从原地震记录数据体中计算出反射系数，将反射系数与压缩后的地震子波进行卷积，可以得到分辨率提高了的地震记录剖面图，可将其应用于储层定位，地震属性分析等方面，利用分辨率提高了的地震记录剖面更容易识别小尺度的地质目标。地震记录剖面更容易识别小尺度的地质目标。地震记录剖面更容易识别小尺度的地质目标。

【技术实现步骤摘要】
一种用谱反演提高地震分辨率方法


[0001]本专利技术涉及油气勘探
，具体为一种用谱反演提高地震分辨率方法。

技术介绍

[0002]石油、天然气一直以来都是关系到国家经济发展和国民生活便利的重要能源。随着国民经济的快速发展，对油气的需求也逐年大幅度提升，然而由于前几十年对石油、天然气的大规模开采，较易勘探和开采的常规地理区域的油气几乎消耗殆尽。油气勘探的目标开始瞄准非常规地质构造，例如海储层、小断裂等小尺度地质月标，这些区域的油气开发成为重要研究对象。
[0003]然而在对薄储层的勘探开发中，地震记录剖面成像的精度问题日益凸显，因为通过地震记录剖面图定位海层对后续的井位设计及油藏定位起着至关重要的作用，地质统计并预测技术等，这些技术都在一定程度上获得了一些实际地震资料处理效果，但由于绝大多数处理方法都依赖于原始的地震记录数据，必然受到原地震记录分辨率的影响，处理结果始终不够理想。尤其在于小尺度地质目标识别过程中，若单纯依赖于原始地震数据，则更加难以识别，处理结果偏差更大。

技术实现思路

[0004]本专利技术旨在至少在一定程度上解决上述技术中的技术问题之一。为此，本专利技术的一个目的在于提出一种用谱反演提高地震分辨率方法，通过利用谱反演技术能从原地震记录数据体中计算出反射系数，将反射系数与压缩后的地震子波进行卷积，可以得到分辨率提高了的地震记录剖面图，可将其应用于储层定位，地震属性分析等方面，利用分辨率提高了的地震记录剖面更容易识别小尺度的地质目标。
[0005]本专利技术的第二个目的在于提出一种用谱反演提高地震分辨率装置。
[0006]本专利技术的第三个目的在于提出一种用谱反演提高地震分辨率设备。
[0007]本专利技术的第四个目的在于提出一种计算机可读存储介质。
[0008]为实现上述目的，本专利技术提供如下技术方案：一种用谱反演提高地震分辨率方法，包括：
[0009]获取地震资料并载入地震数据；
[0010]从地震数据中提取出地震小波数据，并变换地震小波数据到频域W(f)，将地震数据变换到频域S(f)；
[0011]通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数；
[0012]通过谱分解对目标函数进行分析。
[0013]作为优选，所述从地震数据中提取出地震小波数据，包括从地震数据中通过统计法提取地震子波，并通过多道统计法提取最小相位子波。
[0014]作为优选，所述最小相位子波的获取步骤为：
[0015]设立子波频谱域表达式：W(ω)＝|W(ω)|e
jπω
；
[0016]由子波谱频域表达式可得子波振幅谱从而可以得到子波对数谱，则有子波对数频谱的频域表达式为:
[0017][0018]由InW(ω)求解地震子波的相位谱φ(ω)公式如下:
[0019][0020]通过傅里叶反变换计算时域的地震子波w(t)，由W(ω)及队φ(ω)可得:
[0021][0022]作为优选，所述通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数中，所述目标函数如下：
[0023][0024][0025]其中，w(t,f)为地震子波的振幅谱，s(t,f)为地震记录的振幅谱，r0(t)为反射系数序列的奇分量，r
e
(t)为反射系数序列的偶分量，f
L
为低截频，f
H
为高截频，T
i
为第i层与底N
‑
i+1层之间的间隔，a
e
和a0都是根据信噪比高低来确定其比值的加权系数。
[0026]作为优选，所述通过谱分解对目标函数进行分析中，对谱反演算法进行验证，预设薄互层反射系数模型，选择雷克子波进行合成地震记录进行分析，通过原始反射系数模型与反演结果进行对比，来验证谱反演的正确性及合理性，以及是否能够识别薄层。
[0027]作为优选，所述地震子波通过傅立叶变换将时间域信号转换到频率域中，利用共轲梯度法求取最优解。
[0028]本专利技术还提供一种用谱反演提高地震分辨率装置，包括：
[0029]数据载入模块，其用于获取地震资料并载入地震数据；
[0030]数据转换模块，其用于从地震数据中提取出地震小波数据，并变换地震小波数据到频域W(f)，将地震数据变换到频域S(f)；
[0031]数据分析模块，其用于通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数；
[0032]数据分析模块，还用于通过谱分解对目标函数进行分析。
[0033]本专利技术还提供一种用谱反演提高地震分辨率设备，所述用谱反演提高地震分辨率设备为实体设备，所述用谱反演提高地震分辨率设备包括：
[0034]处理器、存储器，所述处理器、存储器与处理器进行通信连接；
[0035]所述存储器用于储存至少一个所述处理器执行的可执行指令，所述处理器用于执行所述可执行指令以实现如上述的用谱反演提高地震分辨率方法。
[0036]本专利技术还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序，所述计算机程序被处理器执行时实现如上述的用谱反演提高地震分辨率方法。
[0037]与现有技术相比，本专利技术的有益效果是：
[0038]本专利技术通过利用谱反演技术能从原地震记录数据体中计算出反射系数，将反射系数与压缩后的地震子波进行卷积，可以得到分辨率提高了的地震记录剖面图，可将其应用于储层定位，地震属性分析等方面，利用分辨率提高了的地震记录剖面更容易识别小尺度的地质目标。
附图说明
[0039]图1为本专利技术实施例提供的一种用谱反演提高地震分辨率方法的主流程图；
[0040]图2为本专利技术实施例提供的一种用谱反演提高地震分辨率装置的结构框图。
具体实施方式
[0041]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0042]本实施方式的方法的执行主体为终端，所述终端可以为手机、平板电脑、掌上电脑PDA、笔记本或台式机等设备，当然，还可以为其他具有相似功能的设备，本实施方式不加以限制。
[0043]请参阅图1，本专利技术提供一种用谱反演提高地震分辨率，所述方法应用于地震信号谱反演，包括：
[0044]S101、获取地震资料并载入地震数据；
[0045]S102、从地震数据中提取出地震小波数据，并变换地震小波数据到频域W(f)，将地震数据变换到频域S(f)；
[0046]S103、通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数；
[0047]S104、通过谱分解对目标函数进行分析。
[0048]其中，地震信号谱反演技术有它的理论优势，在理想情况下，即地震子波已知，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用谱反演提高地震分辨率方法，其特征在于，包括：获取地震资料并载入地震数据；从地震数据中提取出地震小波数据，并变换地震小波数据到频域W(f)，将地震数据变换到频域S(f)；通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数；通过谱分解对目标函数进行分析。2.根据权利要求1所述的用谱反演提高地震分辨率方法，其特征在于，所述从地震数据中提取出地震小波数据，包括从地震数据中通过统计法提取地震子波，并通过多道统计法提取最小相位子波。3.根据权利要求1所述的用谱反演提高地震分辨率方法，其特征在于，所述最小相位子波的获取步骤为：设立子波频谱域表达式：W(ω)＝|W(ω)|e
jπω
；由子波谱频域表达式可得子波振幅谱从而可以得到子波对数谱，则有子波对数频谱的频域表达式为：由In|W(ω)求解地震子波的相位谱φ(ω)公式如下:通过傅里叶反变换计算时域的地震子波w(t)，由W(ω)及队φ(ω)可得:4.根据权利要求1所述的用谱反演提高地震分辨率方法，其特征在于，所述通过反射系数频域表达式对频域进行奇偶分解，并获取目标函数中，所述目标函数如下：频域表达式对频域进行奇偶分解，并获取目标函数中，所述目标函数如下：其中，w(t,f)为地震子波的振幅谱，s(t,f)为地震记录的振幅谱，r0(t)为反射系数序列的奇分量，r
e
(t)为反射系数序列的偶分量，f
L
为低截频，f
H
为高截频，T
i
为第i层与底N
‑

【专利技术属性】
技术研发人员：于建群，聂可可，
申请(专利权)人：吉奥科思北京能源科技有限公司，
类型：发明
国别省市：