逆时偏移与波形反演源波场重建问题中的检查点方法

摘要

对于叠前逆时偏移成像和全波形反演来说,为了满足互相关成像条件以及计算模型校正梯度方向时对于震源波场和检波点波场延拓方向一致性的要求,需要对其中一方进行波场重建.通常情况下这一步骤针对震源波场进行,称之为'源波场重建'.对于源波场重建问题目前主要的解决方法有两种:根据地震波场正演传播算法推导出相应的波场反传计算公式，再配合‘存储/重载’吸收边界位置一定层厚的波场信息的方法实现，本文中将其简称为‘显式反传公式法’或‘显式公式法’。利用Symes教授提出的检查点方法(checkpointing)实现。经过作者的研究、实践，这种方法算法实现简单，适用性强，并且重建波场与原始正传波场之间理论上应该不存在任何量级的误差。作者目前已利用检查点方法实现了在时间离散方面采用四阶精度LDDRK方案的谱元法地震叠前逆时深度偏移成像算法，并将其应用于了北京地区陆上爆炸震源的实际反射地震数据处理，取得了良好的效果。