[发明专利]一种反演方法、系统、装置及存储介质

权利要求书

1.一种反演方法，其中，该方法包括：

利用基于Fatti近似的叠前稀疏层反演获取目标区尖脉冲化的纵波阻抗反射率反演结果；

基于所述纵波阻抗反射率反演结果，确定非零反射率的垂向位置向量；

基于所述非零反射率的垂向位置向量，构建分界面约束矩阵；

以构建的分界面约束矩阵作为约束进行基于Russell近似的最小二乘反演获取流体因子反演结果、剪切模量反演结果和密度反演结果中的一种或两种的组合。

2.根据权利要求1所述的反演方法，其中，所述利用基于Fatti近似的叠前稀疏层反演获取目标区尖脉冲化的纵波阻抗反射率反演结果包括：

基于Fatti近似方程和角度子波褶积矩阵确定第一AVO正演矩阵；

基于所述第一AVO正演矩阵并引入反射系数奇偶分解理论，构建关于弹性参数反射率对应的奇偶分量系数的第一地震数据拟合差项；

基于反射系数奇偶分解理论，构建关于弹性参数反射率对应的奇偶分量系数的稀疏层约束项；

基于所述第一地震数据拟合差项、所述稀疏层约束项构建基于Fatti近似的叠前稀疏层反演目标函数即基于Fatti近似的关于弹性参数反射率对应的奇偶分量系数的叠前反演目标函数；

基于构建的所述目标函数进行反演得到纵波阻抗反射率对应的奇偶分量系数；

基于所述纵波阻抗反射率对应的奇偶分量系数，得到所述纵波阻抗反射率数据体即为所述目标区尖脉冲化的纵波阻抗反射率反演结果。

3.根据权利要求2所述的反演方法，其中，所述第一AVO正演矩阵为：

其中，G_Fatti为基于Fatti近似方程和角度子波褶积矩阵确定的第一AVO正演矩阵；W(θ₁),W(θ₂),…W(θ_K)为角度子波褶积矩阵，K指K个所述角度叠加地震数据；A(θ_j)、B(θ_j)和C(θ_j)均为对角矩阵，其中j＝1,2,…K；θ_j为第j个入射角，其中j＝1,2,…K；其中，

4.根据权利要求2或3所述的反演方法，其中，所述第一地震数据拟合差项为：||d-G_FattiΓx||²；

其中，G_Fatti为由Fatti近似和子波褶积矩阵确定的第一AVO正演矩阵；d为角度叠加地震数据；Γ为由奇偶分解矩阵确定的块化矩阵；x为弹性参数反射率对应的奇偶分量系数；其中，

K指K个所述角度叠加地震数据；θ_j为第j个入射角，其中j＝1,2,…K；

其中D为奇偶分解矩阵；其中，

D＝[D_e,D_o]；其中，

其中，x_I为纵波阻抗反射率对应的奇偶分量系数，x_J为横波阻抗反射率对应的奇偶分量系数，x_ρ为密度反射率对应的奇偶分量系数；其中，

其中，x_{I_e}为纵波阻抗反射率对应的偶分量系数，x_{I_o}为纵波阻抗反射率对应的奇分量系数；

其中，x_{J_e}为横波阻抗反射率对应的偶分量系数，x_{J_o}为横波阻抗反射率对应的奇分量系数；

其中，x_{ρ_e}为密度反射率对应的偶分量系数，x_{ρ_o}为密度反射率对应的奇分量系数。