[发明专利]基于谱分解提高叠后地震资料分辨率的方法

说明书

技术领域

本发明涉及石油地震勘探，属于地震勘探资料处理领域，更具体地讲，涉及一种基于谱分解提高叠后地震资料分辨率的方法。

背景技术

随着勘探、开发的发展，地震勘探已经从简单的构造勘探转向复杂地表、复杂构造、裂缝、薄储层和老油田剩余油的勘探。薄储层和老油田剩余油的勘探和地震成像分辨率有直接的关系。目前，地震子波的主频大约在30Hz左右，对于25米以下的薄层就很难区分开它的顶部反射和底部反射同相轴。但是工业的目标地质体的厚度一般在10米甚至更小。地震同相轴追踪不足以对储层的顶底界面作精细地刻画。谱模拟等常规反褶积虽然是实际地震资料处理中应用最为广泛的提高地震资料分辨率的方法，但也很难以分辨10米以下的目标。因此，亟需开发一种能够提高地震资料分辨率、改善地震资料分辨薄互层的能力的方法。

发明内容

针对现有技术中存在的不足，本发明的目的之一在于解决上述现有技术中存在的一个或多个问题。例如，本发明的目的之一在于提供一种基于谱分解提高叠后分辨率方法，使得提高分辨率后的地震剖面能分辨分辨率更高的薄层。

为了实现上述目的，本发明提供了一种基于谱分解提高叠后地震资料分辨率的方法。所述方法包括以下步骤：a、获得叠后地震记录；b、在叠后地震记录的二次谱中分离提取出地震子波的二次谱，并构成地震记录的零相位子波；c、根据反射系数奇偶分解理论，利用谱分解技术建立基于多层模型的地震记录反射系数谱分解目标函数；d、利用甚快模拟退火非线性优化算法对所述地震记录反射系数谱分解目标函数进行全局寻优，反演出地震记录反射系数；e、提高所述地震记录反射系数的偶分量比重，并将提高偶分量比重后的地震记录反射系数与所述地震记录的零相位子波进行褶积，形成提高分辨率后的地震记录。

根据本发明基于谱分解提高叠后地震资料分辨率的方法的一个实施例，所述步骤b包括：

B1、根据忽略噪声的地震记录褶积模型，得到下面的等式(1)：

A_x(w)＝A_w(w)×A_r(w)  (1)

在等式(1)中，A_x(w)、A_w(w)、A_r(w)分别为地震记录的振幅谱、地震子波的振幅谱和反射系数序列的振幅谱，×表示乘积运算；

B2、对等式(1)两边进行傅立叶(傅里叶)正变换并取模，得到下面的等式(2)：

Ax2(w)=Aw2(w)*Ar2(w)---(2)]]>

在等式(2)中，及分别为地震记录的二次谱、地震子波的二次谱以及反射系数序列的二次谱，*表示褶积运算；

B3、利用低通滤波器对地震记录的二次谱进行低通滤波以剔除反射系数序列的二次谱，保留地震子波的二次谱；

B4、根据步骤B3提取出的地震子波的二次谱，构成地震记录的零相位子波。

根据本发明基于谱分解提高叠后地震资料分辨率的方法的一个实施例，

所述步骤c包括：

C1、按下面的等式(3)～(5)对反射系数序列r(t)进行奇偶分解，

r(t)＝r_o(t)+r_e(t)  (3)

r_o(t)＝[r(t)-r(-t)]/2  (4)

r_e(t)＝[r(t)+r(-t)]/2  (5)

在等式(3)～(5)中，r_o(t)和r_e(t)分别为反射系数序列的奇分量和偶分量；