1.一种方法，所述方法包括：
在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应，并且所述储层包括多个反射面，其中每个反射面具有反射系数；
对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的共聚焦点(cfp)算子来确定所述多个反射面的反射系数；
比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化；和
显示所述反射系数的所述变化。
2.根据权利要求1所述的方法，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。
3.根据权利要求1所述的方法，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和所述多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子，
f(i+1)(z0，zm)＝f(i)(z0，zm)δf(i)(z0，zm)
和
其中δt(i)是基于在第i次迭代中的差异时移(dts)道集所确定的剩余传播时间。
4.根据权利要求3所述的方法，所述方法还包括，对于所述特定反射面zm和所述特定采集表面位置z0，基于在所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。
5.根据权利要求1所述的方法，所述方法还包括，对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。
6.根据权利要求1所述的方法，其中所述多个采集平面位置包括震源位置和地震检波器位置。
7.根据权利要求1所述的方法，其中所述储层包含烃液体。
8.一种系统，所述系统包括：
计算机存储器；和
一个或多个硬件处理器，所述一个或多个硬件处理器可互操作地与所述计算机存储器连接并且被配置成进行包括以下各项的操作：
在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应，并且所述储层包括多个反射面，其中每个反射面具有反射系数；
对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的共聚焦点(cfp)算子来确定所述多个反射面的反射系数；
比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化；和
显示所述反射系数的所述变化。
9.根据权利要求8所述的系统，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。
10.根据权利要求8所述的系统，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和所述多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子，
f(i+1)(z0，zm)＝f(i)(z0，zm)δf(i)(z0，zm)
其中δt(i)是基于在第i次迭代中的差异时移(dts)道集所确定的剩余传播时间。
11.根据权利要求10所述的系统，其中所述操作还包括，对于所述特定反射面zm和所述特定采集表面位置z0，基于在所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。
12.根据权利要求8所述的系统，其中所述操作还包括，对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。
13.根据权利要求8所述的系统，其中所述多个采集平面位置包括震源位置和地震检波器位置。
14.根据权利要求8所述的系统，其中所述储层包含烃液体。
15.一种非暂时性的计算机可读介质，所述非暂时性的计算机可读介质存储一个或多个由计算机系统可执行以进行包括以下各项的操作的指令：
在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应，并且所述储层包括多个反射面，其中每个反射面具有反射系数；
对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的共聚焦点(cfp)算子来确定所述多个反射面的反射系数；
比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化；和
显示所述反射系数的所述变化。
16.根据权利要求15所述的非暂时性的计算机可读介质，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。
17.根据权利要求15所述的非暂时性的计算机可读介质，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和所述多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子，
f(i+1)(z0，zm)＝f(i)(z0，zm)δf(i)(z0，zm)
其中δt(i)是基于在第i次迭代中的差异时移(dts)道集所确定的剩余传播时间。
18.根据权利要求17所述的非暂时性的计算机可读介质，其中所述操作还包括，对于所述特定反射面zm和所述特定采集表面位置z0，基于在所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。
19.根据权利要求15所述的非暂时性的计算机可读介质，其中所述操作还包括，对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。
20.根据权利要求15所述的非暂时性的计算机可读介质，其中所述多个采集平面位置包括震源位置和地震检波器位置。