[发明专利]一种高精度的纵横波阻抗反演方法

权利要求书

1.一种高精度的纵横波阻抗反演方法，其特征在于，包括：

通过对纵、横波阻抗和密度的地震反射系数近似关系进行优化推导处理，得到优化后的纵、横波阻抗反射系数近似关系，其中，所述优化后的纵、横波阻抗反射系数近似关系为：

利用地震褶积模型和所述优化后的纵、横波阻抗反射系数近似关系，得到反射系数分解后的地震褶积模型；

利用所述反射系数分解后的地震褶积模型及拉格朗日乘子法，构建纵、横波阻抗低频模型同步约束的反演目标函数；

利用L_1-2最小化算法对所述纵、横波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数，其包括：

通过对所述反演目标函数进行化简处理，得到L_1-2最小化反演目标函数；

利用增广拉格朗日乘子法对所述L_1-2最小化反演目标函数进行约束，得到最终目标函数

利用交替方法乘子法ADMM对所述最终目标函数进行求解，得到奇偶分量系数；

利用所述奇偶分量系数，得到纵、横波阻抗；

其中，θ为入射角，b为实际工区测井资料拟合的参数，分别为t时刻对应的纵波阻抗和横波阻抗，α为权重参数。

2.根据权利要求1所述的方法，其特征在于，所述利用地震褶积模型和所述优化后的纵、横波阻抗反射系数近似关系，得到反射系数分解后的地震褶积模型包括：

通过对所述优化后的纵、横波阻抗反射系数近似关系中的反射系数进行奇偶分解，得到反射系数分解后的纵横、波阻抗反射系数近似关系；

根据地震褶积模型和所述反射系数分解后的纵、横波阻抗反射系数近似关系，得到反射系数分解后的地震褶积模型。

3.根据权利要求2所述的方法，其特征在于，所述反射系数分解后的纵、横波阻抗反射系数近似关系包括：

r＝Dm；

其中，所述r是指纵横波阻抗反射系数；所述D是指反射系数分解算子；所述m是指奇偶分量系数。

4.根据权利要求3所述的方法，其特征在于，所述根据地震褶积模型和所述反射系数分解后的纵、横波阻抗反射系数近似关系，得到反射系数分解后的地震褶积模型包括：

通过将所述反射系数分解后的纵、横波阻抗反射系数近似关系r＝Dm代入所述地震褶积模型s＝WCr中，得到反射系数分解后的地震褶积模型：s＝WCDm；

其中，所述s是指多角度叠前地震记录；所述W是指w(θ)组成的多角度子波矩阵；所述C是指由C_p(θ)、C_s(θ)组成的系数矩阵。

5.根据权利要求4所述的方法，其特征在于，所述纵、横波阻抗同步约束的反演目标函数包括：

其中，所述β₁是指纵波阻抗低频约束的权重系数；所述β₂是指横波阻抗低频约束的权重系数；所述所述所述以及所述B均是指积分算子；所述I_{p_low}是指纵波阻抗低频模型；所述I_{s_low}是指横波阻抗低频模型；所述L_LP为低通滤波矩阵。

6.根据权利要求5所述的方法，其特征在于，所述利用所述反射系数分解后的地震褶积模型及拉格朗日乘子法，构建纵、横波阻抗低频模型同步约束的反演目标函数包括：

根据测井资料，分别建立纵波低频模型和横波低频模型；

根据所述拉格朗日乘子法、所述纵波低频模型以及所述横波低频模型对基追踪标准公式进行同步约束，得到所述反演目标函数。

7.根据权利要求1所述的方法，其特征在于，所述利用所述奇偶分量系数，得到纵、横波阻抗包括：

根据所述奇偶分量系数，得到纵、横波反射系数；

通过对所述纵、横波反射系数进行积分处理，得到纵、横波阻抗。