基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法

作     者：汪玲玲 高静怀 赵谦 徐宗本 姜秀娣 

作者机构：西安交通大学数学与统计学院西安710049 海洋石油勘探国家工程实验室西安710049 西安交通大学电子与信息工程学院西安710049 中海石油研究中心北京100029 

出 版 物：《地球物理学报》 (Chinese Journal of Geophysics)

年 卷 期：2017年第60卷第2期

页      面：639-654页

核心收录：

学科分类：081801[工学-矿产普查与勘探] 081802[工学-地球探测与信息技术] 08[工学] 0818[工学-地质资源与地质工程] 

基　　金：国家自然科学基金项目(41404107 41390454 91330204) 中国博士后科学基金面上项目(2013M532028)共同资助 

主　　题：波阻抗 矩阵分解 稀疏 横向连续性 递推反演 

摘      要：针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果.