[发明专利]宽频时空变地震子波提取方法

说明书

本发明提供一种宽频时空变地震子波提取方法，该宽频时空变地震子波提取方法包括：步骤1，获取地震数据，并对地震数据进行预处理；步骤2，从被预处理后的地震数据中获取地震解释层位，并在地震解释层位中截取层位时窗，层位时窗表示目标区域；步骤3，在目标区域中，划定一个或多个子波空变区域；步骤4，在一个或多个子波空变区域中，划定一个或多个时变滑动延拓时窗；步骤5，将所有子波空变区域中，将所有时变滑动延拓时窗得到的估算地震子波进行组合，得到宽频时空变地震子波。该宽频时空变地震子波提取方法采用高效的算法，同时提供时变、空变地震子波，相比于常规同态域子波提取算法稳定性和精度更高，为高精度地震反演提供支撑。

技术领域

本发明涉及地球物理勘探技术领域，特别是涉及到一种宽频时空变地震子波提取方法。

背景技术

地震子波估计是地震资料反演、正演和反褶积等处理的必要环节，也是地震资料处理领域的重要研究课题之一，而提取子波的准确性是制约处理精度的关键因素之一。

现阶段，国内外地震子波提取方法主要是基于褶积模型，分为确定性和统计性两大类。确定性子波提取方法是指在有测井数据的情况下，利用测井数据和井旁地震道数据，基于褶积模型采用数学方法估算出地震子波，是反演研究中最为常用的子波提取方法之一，该方法的优点是无需对反射系数的分布做出某种假设，但由于需要测井数据的参与，子波提取的精度与测井曲线准确性以及测井地震标定准确性密切相关。统计性子波提取方法的特点是不需要测井曲线的参与，仅应用地震道数据通过统计方法来估算子波，通常需要对地震数据和反射系数分布做出一定的假设，同态域子波提取方法就是一种高阶统计性子波提取方法。

传统的地震子波提取方法基于子波时不变的假设，然而由于地下介质对地震波能量的吸收衰减作用，实际地震资料中的地震子波具有时变特征；同时，地下地质特征的变化以及采集处理的差异也使地震子波具有空变特征。时空变子波的准确提取是实现高精度地震资料反演的必要条件，也是目前迫切需要解决的重要问题。

发明人发现，现有的地震子波提取技术无法有效解决时空变子波的高效、高精度提取问题，因此亟待开发一种宽频时空变地震子波提取方法。

在申请号：CN202010013147.X的中国专利申请中，涉及到一种基于压缩感知的高斯频率域提高薄互层分辨率的方法，包括：进行薄互层正演模拟，分析其地震响应特性；利用傅里叶变换把地震资料从时间域变换到频率域；引入对信号具有去噪和平滑作用的高斯函数，把频率域地震资料变换到高斯频率域；对原始地震资料利用多道复赛谱技术提取地震子波；利用提取的地震子波和高斯频率域地震资料进行基于压缩感知的稀疏反演，得到稀疏反射系数；基于提取的地震子波进行拓频处理得到宽频地震子波，把宽频地震子波与稀疏反射系数进行褶积，获得高分辨率地震资料。

在申请号：CN201810855372.0的中国专利申请中，涉及到一种基于线性扫描的提高薄互层识别能力的方法，包括：通过对地震资料进行频谱和波形特征分析，得到地震资料的优势频带范围和相位信息；利用复赛谱提取每道地震资料的子波，调整所获得的子波频带和相位；利用线性扫描信号的宽频信息,滑动时窗选择最佳线性扫描子波；对于单道记录利用提取子波与最佳线性扫描子波进行子波整形反褶积，得到该目标道的整形算子；利用新的复合整形算子与该目标地震记录道褶积，获得新的高分辨率地震记录。

在申请号：CN201810465628.7的中国专利申请中，涉及到一种智能化时变盲反褶积宽频处理方法及装置，其中，该方法包括以下步骤：从衰减地震数据中提取出初始地震子波；随机产生包含多个地层品质因子Q值的Q值群；根据所述初始地震子波和所述Q值群建立时变子波矩阵群；根据所述衰减地震数据和所述时变子波矩阵群，获得多个反射系数；建立适应度函数；根据所述适应度函数和所述多个反射系数，获得多个适应值，其中，一个适应值对应一个地层品质因子Q值；根据所述多个适应值确定最优地层品质因子Q值。

以上现有技术均与本发明有较大区别，未能解决我们想要解决的技术问题，为此我们发明了一种新的宽频时空变地震子波提取方法。

发明内容