[发明专利]基于深度波数谱的瞬变波包提取方法、存储介质与终端

说明书

本公开提供了一种基于深度波数谱的瞬变波包提取方法、存储介质与终端，通过对深度域地震信号进行积分变换，得到深度波数谱；将每一深度点的波数谱变换到同态域，得到对应的波数复赛谱，经低通滤波后反变换为波数谱；确定波数谱中每一深度点的波数谱最大值，根据该最大值对应的波数，计算对应深度的瞬变波包；合成深度域地震记录，并与观测地震记录对比，若差别超过设定值则返回积分变换，改变积分变换的调节因子，重新得到深度波数谱，直到最终计算的差值满足设定值。

技术领域

本公开涉及一种基于深度波数谱的瞬变波包提取方法、存储介质与终端。

背景技术

本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。

随着叠前深度偏移技术的发展，深度域地震资料越来越多。与时间域资料相比，深度域地震资料具有构造真实直观的特点，可以综合利用测井资料得到更精确的反映地下地层情况的资料。在深度地震资料反演中，影响反演过程的至关重要一步是深度域合成地震记录，进行深度域层位标定。深度域层位标定是联系深度域地震资料和测井资料的桥梁，其准确与否取决于深度域合成地震记录的好坏，而精确的深度域合成地震记录依赖于准确的地震瞬变波包的提取。

目前，地震波包提取方法一般是先通过时间域地震资料提取地震子波，然后通过时深转换及深度域重新采样，将时间域的子波转换为深度域的波包，然而由于时间域和深度域地震的非线性特征，用这种方法提取到的波包用于深度域合成记录误差较大。同时，由于地层的吸收衰减作用，地震波包的波数在每个深度点是不一样的，随深度变化而变化。在深度域用同一个波包进行合成地震记录不符合实际情况，因此合成地震记录与观测地震记录相似度不高，影响了深度域反演的准确性。

林伯香在《利用变换深度域速度函数制作深度域合成地震记录》一文中指出，深度域波包是地层速度的函数，其传播过程不满足线性时不变条件，必须对深度域速度函数进行适当的变换，使其满足褶积方法计算合成地震记录的条件。此后，很多学者通过不同的方法变换速度函数，在深度域提取波包并合成地震记录。例如张杰在《基于Gibbs抽样方法的深度域地震子波提取》一文中提出基于褶积理论，首先对层速度和层厚度用一个常速度进行变换，在变换后的常速度下，利用抽样方法进行波包提取；陈学华等人在《深度域地震子波的优化提取方法》中提出一种深度域速度变换与脊回归算法相结合的深度域地震子波提取方法。此类深度域波包提取方法都依赖地层速度模型，速度变换方法不同则波包提取结果不同，在实际应用中还存在诸多问题。

发明内容

本公开为了解决上述问题，提出了一种基于深度波数谱的瞬变波包提取方法、存储介质与终端，本公开利用深度域地震记录的深度波数谱来提取瞬变波包，实现了深度域瞬变波包的精确提取。利用此瞬变波包通过瞬变褶积模型合成深度地震记录，与实际地震记录匹配良好。

根据一些实施例，本公开采用如下技术方案：

一种基于深度波数谱的瞬变波包提取方法，包括：

对深度域地震信号进行积分变换，得到深度波数谱；

将每一深度点的波数谱变换到同态域，得到对应的波数复赛谱，经低通滤波后反变换为波数谱；

确定波数谱中每一深度点的波数谱最大值，根据该最大值对应的波数，计算对应深度的瞬变波包；

合成深度域地震记录，并与观测地震记录对比，若差别超过设定值则返回积分变换，改变积分变换的调节因子，重新得到深度波数谱，直到最终计算的差值满足设定值。

作为进一步的限定，对深度域地震信号进行积分变换的具体过程为：

式中S(γ,ω)表示深度波数谱；γ表示目标深度点；ω表示圆波数；s(h)表示深度域地震信号；h表示对应的深度范围；W(h,ω)是深度域高斯窗函数：