[发明专利]一种等长连续时变地震子波获取方法与装置

说明书

本发明提供了一种等长连续时变地震子波获取方法与装置，包括：输入地震数据并按预设参数从地震数据中提取地震波形序列，按预设参数从地震波形中提取地震子波并组成时变地震子波序列，从时变地震子波序列中依次抽取两个相邻的地震子波并组成地震子波对序列，依次对所有地震子波对进行变形处理并形成对应的渐变地震子波序列，构建并输出全时域连续时变地震子波序列。本发明专利提出的一种等长连续时变地震子波获取方法和装置，利用二维图像变形算法对地震子波进行插值，形成连续时变子波序列，既实现了地震成像空间中的地震子波是时变的，又保证了地震子波的变化是连续的，且不至于引起计算资源的大量消耗。

技术领域

本发明涉及地震数据处理领域，具体而言，涉及一种等长连续时变地震子波获取方法与装置。

背景技术

针对中国发明专利/发明专利(申请号：CN201710501413.1)所提出的一种一种地震子波提取方法，该方法可以包括：基于井资料与地震资料，获得多个井的反射系数序列与对应的井旁地震道，其中多个井的反射系数序列以反射系数矩阵表示，对应的井旁地震道以井旁地震道矩阵表示；基于反射系数矩阵与井旁地震道矩阵，建立关于地震子波矩阵的字典学习方程；基于反射系数矩阵、井旁地震道矩阵与地震子波矩阵，建立目标函数，计算使目标函数的值最小的地震子波矩阵，作为最终的地震子波矩阵；对最终的地震子波矩阵进行提取与处理，获得平均子波，可以同时对多个井进行精确高效地地震子波提取。

现有方法在实现子波整形预先获得地震成像空间中的每一个采样点的地震子波，将这些地震子波改造为标准波形，但是现有的地震子波在获取的过程中一方面是由于地层的反射、透射和散射导致地震子波互相叠加和干涉，单一地震子波基本上是不存在的，即地震子波常常是未知的，需要采用一定的数学算法从地震记录中近似地对地震子波进行估算，另一方面，地震子波的估算主要依赖于反演算法，而反演算法不可避免地导致运算量巨大，对地震成像空间中的每一个采样点处的地震子波进行估算是不现实的，在工程实践中，常常假设地震子波是非时变的，一个地震道仅估算一个子波，随着计算机运算能力的提高，很多地球物理学者尝试采用有限数量的时变子波进行地震处理，由于子波距离过大，这些子波的形状可能会差别很大，会导致处理后的地震响应的连续性发生突变，针对地震处理中地震子波非时变或有限时变假设所造成的地震资料分辨率降低和地震响应连续性变差等影响，提出一种等长连续时变地震子波获取方法与装置。

发明内容

本发明的目的在于：针对目前存在的背景技术提出的问题，为了实现利用二维图像变形算法对地震子波进行插值，形成连续时变子波序列，既实现了地震成像空间中的地震子波是时变的，又保证了地震子波的变化是连续的，且不至于引起计算资源的大量消耗上述发明目的，本发明提供了以下技术方案：一种等长连续时变地震子波获取方法与装置，包含以下步骤：

S101输入地震数据并按预设参数从地震数据中提取地震波形序列；

S102按预设参数从地震波形中提取地震子波并组成时变地震子波序列；

S103从时变地震子波序列中依次抽取两个相邻的地震子波并组成地震子波对序列；

S104依次对所有地震子波对进行变形处理并形成对应的渐变地震子波序列；

S105构建并输出全时域连续时变地震子波序列。