[发明专利]一种地震数据拓频处理方法、装置、储存介质及电子设备

说明书

本发明公开了一种地震数据拓频处理方法，包括：将地震数据经过薄层匹配追踪处理，得到初始反射系数模型；将所述初始反射系数模型进行稀疏反演处理，得到精准反射系数模型；将所述地震数据进行拓展频率处理，同时得到高频子波数据；将所述高频子波数据与所述反射系数模型褶积得到地震数据拓频结果。该方法稳定性高，得到的反射系数更加稀疏可靠；新剖面频谱拓展不失自然，在处理过程中不损失数据振幅特征。

技术领域

本发明涉及地址勘测技术领域，尤其涉及一种地震数据拓频处理方法、装置、储存介质及电子设备。

背景技术

随着勘探程度的日益加深，对于薄层、精细构造等识别能力的要求日益提高。地震数据是有效描述地下空间介质的手段。但是在地震数据采集过程中，由于受到激发条件和地层吸收的影响，得到的地震主频往往较低；而且在处理过程中，基于褶积模型的反射波处理方法受到其分辨率的影响，分辨能力仅能够达到1/4波长。

随着高分辨处理技术的发展，目前已经发展了测井约束反演、拓频处理、反褶积、小波分析等等方法，还有时频域联合反演反射系数方法、匹配追踪滤波方法、基于盲源反褶积等等技术。

但是这些方法由于未充分利用反射系数的稀疏性和全频段性，对于复杂地区的拓频能力十分有限，且在大多地区难以做到保幅保真处理。提高地震数据分辨率的同时，得到保幅保真真高分辨率的地震数据是亟待解决的问题。

发明内容

(一)发明目的

本发明的目的是提供一种地震数据拓频处理方法、装置、储存介质及电子设备以解决现有技术分辨率低且不保真的问题。

(二)技术方案

为解决上述问题，本发明的第一方面提供了一种地震数据拓频处理方法，包括：

将地震数据经过薄层匹配追踪处理，得到初始反射系数模型；

将所述初始反射系数模型进行稀疏反演处理，得到精准反射系数模型；

将所述地震数据进行拓展频率处理，同时得到高频子波数据；

将所述高频子波数据与所述反射系数模型褶积得到地震数据拓频结果。

进一步地，所述将地震数据经过薄层匹配追踪处理，得到初始反射系数模型具体为：

假设两个单反射层位于时间t_a和t_b，反射系数为a和b，其中a,b均为时间t的时间序列，w(t_a)和w(t_b)分别代表两个基函数，该基函数为小波基或科沃莱特基函数；

从目标方程求解a和b，式中d(t)为地震数据：

obj＝argmin||d(t)-a(t)w(t_a)-b(t)w(t_b)||₂；

argmin||·||₂代表求解L2范数；

能够得到解

得到初始反射系数模型r(t)＝a(t)+b(t)，其中，t为时间。

进一步地，所述将所述初始反射系数模型进行基于L0范数的稀疏反演处理，得到精准反射系数模型具体为：

目标方程为：

obj＝argmin||R(t)||₀ s.t.d(t)＝R(t)w(t)；

d(t)为输入，求解L0范数，得到精准反射系数模型R(t)