本发明涉及生成感兴趣的地下区域的图像的方法和系统。一般说来，本发明的一个实施例包括建立利用地震数据的边界条件和包括来自地球模型中的源位置的激发的初始条件。然后，源波场通过地球模型向前传播达到最大时间，并且在时间上稀疏的多个检验点上被保存，并且每个时间步长上的源波场的相应边界值被保存。此外，利用多个可用检验点和每个时间步长上的保存边界值从最大时间开始通过地球模型向后传播源波场。从最大时间开始通过地球模型同时向后传播接收器波场。成像条件在所选时间步长上被应用于向后传播的源波场和接收器波场两者，并且将那些波场被用于生成与地下区域有关的图像。