一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法与流程

技术特征：

1.一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法，其特征在于包括：

采用弯曲射线追踪法将经过叠前偏移的共成像点道集从偏移距域转换到射线参数域；对所有的共成像点道集转换后，通过抽取常数射线参数的数据和部分数据叠加构建共射线参数剖面；对共射线参数剖面迭代反演得到高精度地震反射系数序列；在射线参数域中采用广义线性反演算法对反射系数道进行反演求取射线波阻抗；采用能够反映岩性参数随频率变化特征的射线阻抗，结合对其它多种弹性阻抗参数的综合分析，联合探区储层物性特征和地质特征，开展基于统计岩石物理的岩性定量解释方法，实现不整合圈闭储层的岩性预测。

2.根据权利要求1所述的基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法，其特征在于，

(1)将地震采集得到并经过偏移成像的CIP道集通过弯曲射线追踪法，从原本的偏移距域转换到射线参数域,根据下面公式通过迭代计算得到共成像点道集中各个反射点的水平方向坐标x_i(i＝0,1,…,2k)

其中u_i为水平慢度，遵从u_i＝u_2k-i对称，其中i＝0,1,…,2k；s_i是点x_i和x_i-1 之间的射线线段长度；

利用sinc方程对射线参数域的波形进行空间插值重构，得到对应射线参数域的共成像点道集；

(2)对所有的共成像点道集转换后，通过抽取常数射线参数的数据和部分数据叠加构建共射线参数剖面；

(3)通过优化Lp范数和柯西约束两种稀疏反演方法将稀疏约束与地震反射系数分布相结合，通过迭代反演得到高精度地震反射系数序列；

(4)在射线参数域中，将反射系数剖面上的相对振幅值转换成实际的反射系数值后采用广义线性反演算法对反射系数道进行反演求取射线波阻抗；

将从CRP剖面反演得到的反射系数剖面用于反演，求取射线波阻抗，沿着一条具有常数射线参数值的射线路径，其反射系数可以表示为

r(t_i；p)＝f(RI_i(p),RI_i+1(p))， (2)

将该表达式进行Taylor级数展开，并表示为矩阵―向量的形式：

r(RI)≈r(RI₀)+G·(RI-RI₀), (3)

式中RI是待求解的波阻抗剖面，RI₀是波阻抗剖面的初始模型，G是反射系数对波阻抗的偏微分矩阵，求解此一阶Taylor展开式得到模型修正量δRI＝RI-RI₀；

通过线性迭代方式反演最终得到对应不同射线参数的射线波阻抗；

(5)利用反演得到的对应不同射线参数的对岩性敏感的射线波阻抗，结合其它类型阻抗，基于统计岩石物理方法，定量地震属性与储层描述之间的关系，预测岩性。