基于L 1/2 正则化理论的地震稀疏反褶积

地震反褶积是一种重要的压缩地震子波、提高薄层纵向分辨率的地震数据处理方法。在层状地层的假设下, 反射系数可视作稀疏的脉冲序列, 所以地震反褶积可以描述为一个稀疏求解问题, L1正则化被广泛用于解决稀疏问题, 但近年来一些文献证明L1正则化的稀疏表达能力不是最优的。针对这一问题, 基于快速发展的L1/2正则化理论, 提出将L1/2正则化作为反射系数的稀疏约束进行地震反褶积处理, 并使用其特定的阈值迭代算法进行求解, 对单道模型的测试证实了该方法对正则化参数和噪声有较好的适应能力。简单二维模型和Marmousi2模型数据的测试结果表明, 基于该方法的反演结果能较好地拟合反射系数振幅, 并且对噪声干扰的鲁棒性更强, 能够更好地保护弱反射系数。实际数据应用结果表明, 该方法能有效消除子波影响, 较好地分辨出薄层结构和透镜体结构, 为地震数据高分辨处理提供了有力工具。