基于TV约束和Toeplitz矩阵分解的波阻抗反演

摘要

利用波阻抗剖面的非高斯分布特点以及地震子波褶积矩阵的Toeplitz结构,对波阻抗剖面进行全变分(TV)约束,可以在压制随机噪声的同时保持剖面的不连续性,对地震子波褶积矩阵进行Toeplitz稀疏矩阵分解得到地震子波的稀疏表达.地震资料的低频损失导致无法反演出波阻抗的低频背景,故将测井或解释层位信息通过最小二乘法建立约束条件.建立的约束最优化目标函数可以同时反演子波和波阻抗.文中将地震剖面整体处理,反演结果比常规的逐道反演具有更高的精度、横向连续性和抗噪性.