[发明专利]基于倾角约束的多道反褶积处理方法

权利要求书

1.基于倾角约束的多道反褶积处理方法，其特征在于：包括以下步骤：

S1：根据原始地震数据d(t,x)计算地层倾角θ(t,x)；

S2：根据地层倾角θ(t,x)计算反射系数r(t,x)沿倾角方向的导数

S3：对反射系数r(t,x)沿倾角方向的导数进行L2范数约束，构建具有倾角约束的多道反褶积目标函数J(r)；

S4：求解反褶积目标函数J(r)，得到反射系数迭代求解公式。

2.根据权利要求1所述的基于倾角约束的多道反褶积处理方法，其特征在于：步骤S1具体如下：

S10：计算原始地震数据d(t,x)沿时间方向的偏导数

S11:计算原始地震数据d(t,x)沿空间方向的偏导数

S12：按照公式：计算地层倾角θ(t,x)。

3.根据权利要求1所述的基于倾角约束的多道反褶积处理方法，其特征在于：步骤S2具体如下：

S20:根据地层倾角θ(t,x)计算反射系数r(t,x)沿时间方向的偏导数在倾角方向的投影

S21：根据地层倾角θ(t,x)计算反射系数r(t,x)沿空间方向的偏导数在倾角方向的投影

S22：根据反射系数r(t,x)沿时间方向的偏导数在倾角方向的投影和沿空间方向的偏导数在倾角方向的投影计算反射系数r(t,x)沿倾角方向的导数

4.根据权利要求3所述的基于倾角约束的多道反褶积处理方法，其特征在于：步骤S20中，反射系数r(t,x)沿时间方向的偏导数在倾角方向的投影步骤S21中，反射系数r(t,x)沿空间方向的偏导数在倾角方向的投影步骤S22中，反射系数r(t,x)沿倾角方向的导数

5.根据权利要求1所述的基于倾角约束的多道反褶积处理方法，其特征在于：所述反褶积目标函数其中，λ为纵向的调节因子，μ为倾角方向的调节因子。

6.根据权利要求5所述的基于倾角约束的多道反褶积处理方法，其特征在于：λ的取值范围为0.01～0.5，μ的取值范围为0.01～1。

7.根据权利要求5所述的基于倾角约束的多道反褶积处理方法，其特征在于：对反褶积目标函数的求解过程如下：先将反射系数r(t,x)的初始值设置为0，然后利用迭代重加权算法求解目标函数J(r),得到反射系数迭代求解公式：其中,r^k+1为经过第(k+1)步迭代后的离散反射系数序列，W为多道地震子波矩阵，为迭代更新算子，D_θ为倾角约束算子。