[发明专利]一种快速高精度叠后地震阻抗反演方法

权利要求书

1.一种快速高精度叠后地震阻抗反演方法，包括：

通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型；

其特征在于，还包括：

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

其中，所述W为子波矩阵；所述C是指积分矩阵；所述D为反射系数分解算子；所述ξ_low是指纵波阻抗低频模型；所述s为地震记录；所述L为低频滤波矩阵；所述λ为正则化调节参数，用来控制待反演参数的稀疏度；所述α为权重参数；所述m表示奇偶分量对应的系数；

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

2.根据权利要求1所述的方法，其特征在于，所述通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型包括：

通过对所述反射系数序列进行奇偶分解，得到分解后的反射系数序列：r＝Dm；

通过将所述分解后的反射系数序列r＝Dm代入所述地震褶积模型中，该地震褶积模型的表达式为s＝Wr+n，得到反射系数序列分解后的地震褶积模型，该反射系数序列分解后的地震褶积模型的表达式为：s＝WDm+n；

其中，所述r是指反射系数序列；所述D是指反射系数分解算子；所述m是指奇偶分量对应的系数；所述s是指地震记录；所述W是指子波矩阵；所述n为地震噪声序列。

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

根据测井资料，建立纵波阻抗低频模型；

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

4.根据权利要求3所述的方法，其特征在于，所述利用L_1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数包括：

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

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

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

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

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

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