关键词:
弹性阻抗
叠前反演
微分超拉普拉斯
贝叶斯框架
分辨率
摘要:
岩石物理模型的引入,使得我们可以充分利用叠前地震资料进行反演从而得到更加全面的地下介质信息。然而以反射系数方程为基础的弹性阻抗方程中包含的需要反演的参数的增多,不同参数对反射系数贡献量的不同等问题也增加了反演的难度。因此,为了保证多参数反演的稳定性,两步法反演得到了广泛的应用。即首先进行弹性阻抗反演以得到不同入射角、方位角下的弹性阻抗信息,随后以弹性阻抗作为输入进行下一步的多参数反演。弹性阻抗作为波阻抗反演技术的拓展已经十分成熟,但仍然存在解的非唯一性、分辨率不足和稀疏性较差等问题。为此,本文将信号与图像处理领域中被广泛应用的超拉普拉斯约束引入到反演过程当中。该方法能够避免0范数求导问题,同时能够根据数据本身特征选取最优p值,使反演结果获得较高分辨率的同时保证一定的稳定性。并且引入线性化反射系数公式,以求直接通过地震数据反演得到弹性阻抗的反演结果。首先,对线性化弹性阻抗近似公式进行推导并分析其精度。随后,基于贝叶斯反演理论框架,推导了引入微分超拉普拉斯块约束后的弹性阻抗反演目标函数,并选用迭代重加权最小二乘法求解目标函数来获得弹性阻抗反演结果。最后,将新方法结果与传统方法获得的结果进行比较,进而验证方法的准确性以及稳定性。相关系数以及弹性阻抗曲线的对比均有效地证明了所提出方法的反演结果优于传统L1范数和L2范数约束的反演结果。将新方法应用于合成数据以及实际数据的反演当中,证明了其有效性和可行性,可以广泛应用于弹性阻抗的反演流程当中,为后续多参数的反演提供了更加可靠全面的信息。