[发明专利]提取叠前校正道集地震子波的方法

说明书

技术领域

本发明总体来说涉及地震勘探领域。更具体地讲，涉及一种提取叠前校正道集地震子波的方法。

背景技术

地震勘探技术逐渐向高精度、高分辨率方向发展，提取精确的地震子波是地震信号反褶积、波阻抗反演的基础，因此，如何提取精确的地震子波日益成为研究的热点问题。

目前，提取地震子波的方法大都基于自相关函数法或统计性自相关法，然而，自相关函数法需要提出最小相位地震子波的假设，且该方法对于相位是盲目的，不能提取到地震子波的相位信息，因此在此种假设前提下得到的处理结果是不可靠的；此外，统计性自相关地震子波提取方法不计算相位谱，默认为最小相位或者人为假设一个常相位，但是，该方法抗噪声干扰的能力十分有限，只有在高斯白噪声中才能得到相对较为可靠的结果。

综上所述，现有的提取地震子波的方法并不能满足人们对地震子波精度和抗噪能力方面的需求。

发明内容

本发明的目的在于提供一种提取叠前校正道集地震子波的方法，以克服现有技术中地震子波精度低和抗噪能力差的问题。

根据本发明示例性实施例的一方面，提供一种提取叠前校正道集地震子波的方法，包括：(A)对采集的地震记录和测井数据分别进行预处理，以获取地震叠前道集数据和测井道集数据；(B)对地震叠前道集数据和测井道集数据分别进行非线性滤波来获取待校正地震叠前道集数据和待处理测井道集数据；(C)利用混合算法确定每道待校正地震叠前道集数据的相位校正角，其中，所述混合算法为最大方差模法和解析法；(D)利用所述相位校正角分别对相应道的待校正地震叠前道集数据进行相位校正，以获取待处理地震叠前道集数据；(E)基于待处理地震叠前道集数据和待处理测井道集数据，利用各道待处理地震叠前道集数据的预定匹配滤波算子来确定待处理地震叠前道集数据的最终匹配滤波算子；(F)将最终匹配滤波算子与待处理地震叠前道集数据进行匹配滤波，以获取待处理地震叠前道集数据和待处理测井道集数据的匹配道集数据；(G)基于匹配道集数据和待处理测井道集数据来确定地震叠前道集数据的地震子波频谱，并利用所述地震子波频谱获取最优地震子波。

可选地，步骤(C)可包括：(C1)利用待校正地震叠前道集数据来确定待校正地震叠前道集数据的初始模型道；(C2)基于所述初始模型道，利用最大方差模法来获取待校正地震叠前道集数据的参考模型道，其中，所述参考模型道是指相位为零时的初始模型道；(C3)针对待校正地震叠前道集数据中的每道数据，基于参考模型道，利用解析法来确定所述每道数据的相位校正角。

可选地，步骤(E)可包括：(E1)将每道待处理地震叠前道集数据分别与相应的预定匹配滤波算子进行匹配滤波，使得匹配滤波后的每道待处理地震叠前道集数据的实际输出趋近于期望输出的相应道的待处理测井道集数据，并确定匹配滤波后每道数据的实际输出与期望输出的误差；(E2)利用每道数据的误差来获取每道数据的误差能量，基于所述误差能量应用最小二乘法来确定每道数据的匹配滤波算子；(E3)将需要进行处理的各道数据的匹配滤波算子进行平均处理，以得到待处理地震叠前道集数据的最终匹配滤波算子。

可选地，可根据各道待处理叠前道集数据的相关性和信噪比来确定需要进行处理的匹配滤波算子。

可选地，步骤(G)可包括：(G1)利用所述匹配道集数据和待处理测井道集数据来确定地震叠前道集数据的地震子波频谱；(G2)将所述地震子波频谱进行傅里叶反变换，得到时间域的地震叠前道集地震子波；(G3)将时间域的地震叠前道集地震子波进行加窗平滑处理来获取最优地震子波。

可选地，步骤(G1)可包括：(G11)利用所述匹配道集数据和待处理测井道集数据确定地震叠前道集数据的初始子波频谱和谱一致函数；(G12)利用所述谱一致函数对所述初始子波频谱进行无偏项替代转换，来获取地震子波频谱。

在根据本发明示例性实施例的提取叠前校正道集地震子波的方法中，能够有效提高地震子波的精度和抗噪能力，并且不需要最小相位的假设，提高了地震子波的可靠性，为后续的反褶积和反演处理奠定了有利基础。

附图说明

通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中：

图1示出根据本发明示例性实施例的提取叠前校正道集地震子波的方法的流程图；

图2示出根据本发明示例性实施例的非线性滤波前后的地震叠前道集数据的示意图；

图3示出根据本发明示例性实施例的确定待校正地震叠前道集数据的相位校正角的步骤的流程图；