一种基于基追踪方法的叠后声波阻抗反演方法与流程

技术特征：

1.一种基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述方法包括：

S1，输入叠后地震数据，由叠后地震数据进行构造解释；

S2，输入测井数据，提取或给定地震子波，并对叠后地震数据进行标定，同时计算出过井波阻抗：

S3，以构造解释结果作为约束，对S2得到的过井波阻抗进行内插外推，获得初始波阻抗体，进而获得初始反射系数序列；

S4，对于每一道地震数据，由初始反射系数及子波构建目标函数；

S5，用基追踪方法求解目标函数，获得反演后的反射系数；

S6，由反演后的反射系数以及标定结果计算最终的波阻抗。

2.根据权利要求1所述的基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述S2是这样实现的：

输入测井数据中的速度和密度，利用合成记录对叠后地震数据进行层位标定和子波提取；

所述合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果，即：F(t)＝S(t)*R(t)，初始合成地震记录是由反射系数R(t)与标准地震子波S(t)进行褶积得到；利用初始合成地震记录进行初始标定，在初始标定的基础上通过井旁地震道和速度密度测井曲线联合提取子波，子波提取、合成记录制作及标定是一个迭代的过程，经过多次迭代，即能够得到合适的子波和高精度的合成记录；子波会在S4中用到，其中过井波阻抗Z_i是基于测井数据中的速度和密度利用Z_i＝V_i·ρ_i计算出来的，其中V_i是速度，ρ_i是密度，i表示样点序号。

3.根据权利要求2所述的基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述S3是这样实现的：

利用所述过井波阻抗值，用插值方法获得该层位上其他点的波阻抗值，然后用获得初始反射系数序列。

4.根据权利要求3所述的基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述S4中的目标函数如下：

在约束||d^obs-W·R||_p≤γ下求解min||R||_q，其中d^obs是观测数据，W为子波，R为反射系数模型，γ为任意小的一个数。

5.根据权利要求4所述的基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述S5是这样实现的：

取p＝2，q＝0，变成二次规划问题，然后采用拉格朗日乘子来解此问题。

6.根据权利要求5所述的基于基追踪方法的叠后声波阻抗反演方法，其特征在于：所述S6是这样实现的：

得到反射系数以后，用获得声波阻抗。