[发明专利]一种基于混合高阶分数阶ATpV稀疏正则化的地震反演方法

摘要

本发明公开了一种基于混合高阶分数阶ATpV稀疏正则化的地震反演方法，属于地球物理反演及油气储层预测领域，解决现有的采用一阶差分求全变分导致的问题。本发明对地震记录进行预处理后提取层位信息，提取层位信息后基于测井数据获取待反演参数的初始模型；基于待反演参数的初始模型和预处理后的地震记录计算初始波阻抗对数，再构建基于混合高阶分数阶ATpV稀疏正则化的正演模型；基于地震记录、地震子波，结合交替方向乘子方法和正演模型对初始波阻抗对数进行更新，得到更新后的波阻抗对数和用于更新下一个波阻抗对数的参数；判断更新前后的值是否满足条件，若是，更新下一个波阻抗对数，若否，获得待反演参数。本发明用于地震反演。

主权项

1.一种基于混合高阶分数阶ATpV稀疏正则化的地震反演方法，其特征在于，包括如下步骤：步骤1：获取地震数据，地震数据包括地震记录S、地震子波w和测井数据，对地震记录S进行预处理后提取层位信息，提取层位信息后基于测井数据获取待反演参数的初始模型；步骤2：基于待反演参数的初始模型和预处理后的地震记录S₀计算初始波阻抗对数L₀，再基于初始波阻抗对数L₀计算分数阶差分矩阵构建基于混合高阶分数阶ATpV稀疏正则化的正演模型；步骤3：基于预处理后的地震记录S₀、地震子波w，结合交替方向乘子方法和正演模型对L₀进行更新，得到更新后的波阻抗对数Lⁱ⁺¹和用于更新下一个波阻抗对数的参数；步骤4：判断更新前后的值是否满足||Lⁱ⁺¹‑Lⁱ||₂/||Lⁱ||₂＞tol，若是，i＝i+1，再转到步骤3进行下一个波阻抗对数的更新，若否，则根据待反演参数与波阻抗对数的关系获得待反演参数。