使用共聚焦点法的连续地震储层监测的制作方法

优先权要求

本申请要求于2017年3月16日提交的美国专利申请号62/472,407的优先权，将其整个内容通过引用结合于此。

本公开内容涉及地震(矿震，seismic)数据处理。

背景

在一个时间期内的连续储层监测涉及地震深度成像方法，其可以产生快速且准确的图像以与应对所采集的大量数据集。在用于连续储层监测的地震成像中的关键步骤是在其中不能进行直接观察的地球固体内部内的波场的估算。除非在成像前高分辨地震速度模型是可用的，否则基于沿着表面检波器(接收器，receiver)的开放边界记录的地震数据的标准估算通常不足以解释能量在复杂的地下如何传播，或者无法准确预测地表下的多次散射波(多重的)。在一些情况中，对于整个上覆层(overburden)使用严苛的深度定向速度估算方法以产生足够准确的速度模型。

概述

本公开内容描述了使用共聚焦点(cfp)法进行连续地震储层监测的方法和系统，包括计算机实施的方法、计算机程序产品和计算机系统。

在一个实施方式中，在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应。所述储层包括多个反射面(reflector)，其中每个反射面具有反射系数。对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的共聚焦点(cfp)算子来确定所述多个反射面的反射系数。比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化。显示所述反射系数的所述变化。

前述实施方式可以使用计算机实施的方法、非暂时性的计算机可读介质和计算机实施的系统来实施，所述非暂时性的计算机可读介质存储计算机可读指令以进行所述计算机实施的方法，所述计算机实施的系统包括与硬件处理器可互操作地(能共同操作地，interoperably)连接的计算机存储器，所述硬件处理器被配置成进行所述计算机实施的方法/存储在所述非暂时性的计算机可读介质上的指令。

本公开内容中所描述的主题通过高效且有效地处理大量地震数据集以快速且准确地产生地震图像而能够实现四维连续的地震储层监测。所描述的方法基于稳健的迭代算子更新过程，同时不派生详尽的速度模型并且不处理整个采集的三维数据集。所产生的地震图像和所确定的在一个时间期内反射系数的变化可以用于有效的油气勘探，如确定用于油井的钻井参数。其他优点对于本领域普通技术人员来说将会是明显的。

本说明书的主题的一种或多种实施方式的细节在附图和描述中提出。所述主题的其他特征、方面和优点根据该说明书、附图和权利要求书将变得明显。

附图说明

图1是根据一些实施方式的使用共聚焦点(cfp)法进行连续地震储层监测的一种示例方法的流程图。

图2是根据一些实施方式的用于迭代更新一个反射层的cfp聚焦算子的一种示例方法的流程图。

图3a示出了根据一些实施方式的使用一个目标反射面的初始聚焦算子的一组差异时移(dts)道集(gather)。

图3b示出了根据一些实施方式的使用一个目标反射面的更新的聚焦算子的一组dts道集。

图4a示出了根据一些实施方式的对于一个目标反射面的初始聚焦算子。

图4b示出了根据一些实施方式的对于一个目标反射面的在一次迭代后的更新的聚焦算子。

图4c示出了根据一些实施方式的在对于一个目标反射面的初始聚焦算子和更新的聚焦算子之间的差。

图5a示出了根据一些实施方式的在目标水平面的反射系数的归一化均方根(normalizedrootmeansquare，nrms)振幅。

图5b示出了根据一些实施方式的勘探(勘察，survey)2和勘探17之间的nrms振幅差。

图6是示出了根据一些实施方式的一种示例计算机系统的框图，所述示例计算机系统用于提供如本公开内容中所描述的与所述算法、方法、函数、过程、流程和程序相关的计算功能。

各个附图中相同的附图标号和标记指示相同的要素。

详述

以下详述描述了使用共聚焦点(cfp)方法的连续地震储层监测，并且被提供用于使得本领域技术人员能够完成和使用在一个或多个具体实施方式的情况下所公开的主题。可以进行所公开实施方式的多种改变、变化和置换，并且其对本领域技术人员来说是非常明显的。在不背离本公开内容的范围的情况下，所公开的实施方式中限定的一般原理可以适用于其他实施方式和应用。因此，本公开内容并不打算限于所描述或所示出的实施方式，而是旨在符合与所公开的原理和特征一致的最宽范围。

可以使用在一个时间期内收集的地震数据来产生地震图像以监测储层的变化。例如，可以在一个时间期内每月采集一次储层的地震反射数据。可以使用所采集的多个地震数据集来产生多个三维(3d)地震图像以确定储层在所述一个时间期内的变化，例如，烃液体的组成或位置的变化。在高水平，所描述的方法对连续3d地震数据集使用cfp法以基于cfp道集和差异时移(dts)道集迭代更新cfp算子，同时不派生详尽的速度模型。

图1是根据一些实施方式的使用cfp法进行连续地震储层监测的一种示例方法100的流程图。为了表述清楚，以下描述总体上描述了在本公开内容中其他附图的情况下的方法100。例如，方法100可以通过图6中所述的计算机系统来进行。然而，将会理解的是，方法100可以例如视情况通过任何合适的系统、环境、软件和硬件，或者系统、环境、软件和硬件的组合来进行。在一些实施方式中，方法100的多个步骤可以并行地、组合地、循环地或以任何顺序进行。

方法100在框102处开始，其中鉴别储层。例如，储层可以包括在地球内的一个或多个地表下层，并且储层的地理界线可以基于来自用户的输入来鉴别。储层可以包括多个反射面，并且每个反射面可以具有反射系数。例如，储层可以是长5千米(km)、宽5km且深5km的3d区域。如果每个反射面建模为长5米(m)、宽5m且深5m的立方体，则储层可以由10⁹个反射面表示。反射面的位置可以由立方体的中心或立方体中的其他点表示。在一些实施方式中，可以将储层分为多个立方体，立方体角可以表示储层中的反射面。例如，如果将所述区域分为5m乘5m乘5m的立方体，则储层可以由1001³个反射面表示。也可以使用其他方法来确定表示储层的反射面的组。

在框104处，在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应。例如，储层的地震数据可以在18个月内每月采集一次，并且可以采集总计18组地震数据。在地震数据采集期间，可以将许多个检波器(例如，地震检波器或水用检波器)定位在地表之上或之下。震源可以将地震波发送到地球中，并且检波器可以记录被地球内的各个地表下层反射的波。震源可以例如被卡车拖曳，并且在不同位置处产生地震波。例如，震源可以在第一位置处射出第一发射(炮，shot)以使检波器记录反射的波，并且震源移动到第二位置以射出第二发射。在与单个发射对应的各检波器处记录的数据称为道(trace)。例如，如果震源在100,000个不同位置处射出发射并且存在1,000个检波器，则所得地震数据集可以具有10⁸个道。

在框106处，对于每个地震数据集，可以通过迭代更新与反射面和多个采集表面位置相关的cfp聚焦算子来确定储层中的反射面的反射系数。采集表面位置可以包括震源位置和地震检波器位置。在先前的100,000个发射位置和1,000个检波器的储层实例中，存在总计101,000个采集表面位置。

迭代更新聚焦算子可以使用如等式(1)和(2)所示的基于cfp法的单向格林函数(one-waygreen’sfunction)进行数据驱动，

f⁽ⁱ⁺¹⁾(z0，zm)＝f⁽ⁱ⁾(z0，zm)δf⁽ⁱ⁾(z0，zm)，i＝0，1，...(1)

其中cfp算子f⁽ⁱ⁾(z0，zm)描述在第i次迭代时在采集表面位置z0和目标反射面zm之间的传播时间或走时(途径时间，traveltime)，δf⁽ⁱ⁾(z0，zm)将采集表面位置z0和目标反射面zm之间的走时更新δt⁽ⁱ⁾/2，并且ω为角频率。注意，等式(1)和(2)以频域表示，其中算子f⁽ⁱ⁾(z0，zm)为狄拉克δ函数(diracdeltafunction)的傅里叶变换，所述狄拉克δ函数在与采集表面位置z0和目标反射面zm之间的走时对应的时间处具有无穷的值而在其他位置为零。迭代更新操作(算法，operation)也可以通过以下式以时域表示

f⁽ⁱ⁺¹⁾(z0，zm)＝f⁽ⁱ⁾(z0，zm)*δf⁽ⁱ⁾(z0，zm)，i＝0，1，...，(3)

其中*表示卷积，δ为狄拉克δ函数，并且t为初始格林函数的走时。在一些情况下，如之后将会讨论的，等式(2)或(4)中的δt⁽ⁱ⁾为与差异时移(dts)道集中的时间零点的拾取时间偏差(pickedtimedeviation)。

在一些实施方式中，初始cfp算子f⁽⁰⁾(z0，zm)可以为在z0和zm之间的走时的任何粗略估算。在一个典型的实施方式中，在目标水平面处定义的正常时差(nmo)速度可以用于得出初始算子。注意，在之后将会解释的cfp道集中，仍可以存在除来自目标反射面的反射事件以外的反射事件。等走时原理确认，如果算子的运动学是正确的，则在每个cfp道集中的目标反射面的走时将会与相应聚焦算子的时间反演(time-reverse)相同。根据等走时原理，之后将会解释的dts道集可以通过将聚焦算子逐道与其相应cfp道集关联的时间产生。基于对dts道集的自动化拾取，使用等式(1)或(3)迭代更新在第i次迭代的聚焦算子，直到对于目标反射面的每个地表下格点都满足等走时的原理(即，dts道集中的目标事件变平)。等式(2)或(4)中的δt⁽ⁱ⁾为与dts道集中的时间零点的拾取时间偏差。目标导向的叠加图像可以通过叠加来自最终聚焦算子的dts道集来得出。

图2是根据一些实施方式的用于迭代更新一个反射面的cfp聚焦算子的一种示例方法200的流程图。换言之，对于每个反射面进行方法200，并且在先前的10⁹个反射面的储层实例中，可以将方法200进行10⁹次。为了清晰的表述，以下的描述总体上描述了在本公开内容中其他附图的情况下的方法200。例如，方法200可以通过图6所述的计算机系统来进行。然而，将理解的是，方法200可以视情况通过系统、环境、软件、硬件、或者系统、环境、软件和硬件的组合来进行。在一些实施方式中，方法200的多个步骤可以并行地、组合地、循环地或以任何顺序进行。

方法200始于框202，其中迭代计数i初始化为零，并且选择目标反射面zm。

在框204处，对于每个采集表面位置，初始化描述在采集表面位置z0和目标反射面zm之间的走时的算子f⁽⁰⁾(z0，zm)。在先前的101,000个采集表面位置的储层实例中，将对101,000对的z0和zm确定初始算子f⁽⁰⁾(z0，zm)。在一些实施方式中，在目标水平面处定义的正常时差(nmo)速度可以用于得出初始算子。例如，由于目标反射面位置和采集表面位置是已知的，可以基于nmo速度以及在目标反射面与采集表面位置之间的距离得出初始走时。在一些实施方式中，f⁽⁰⁾(z0，zm)可以在频域中等于或者在时域中等于δ(t+tnmo)，其中tnmo为基于nmo速度估算的初始走时。

在框206处，使用f⁽ⁱ⁾(z0，zm)产生cfp道集。在初始迭代中，使用f⁽⁰⁾(z0，zm)产生cfp道集。在一些实施方式中，基于采集表面位置对地震数据进行分类。例如，基于检波器位置对地震数据进行分类。在先前的储层实例中，可以按以下顺序对地震数据中的10⁸个道进行分类：在第一检波器处接收的100,000个道，在第二检波器处接收的100,000个道，等等。来自第一检波器的100,000个道中的每个都与相应的算子f⁽ⁱ⁾(z0，zm)相关联，其中z0为与特定道对应的震源位置。通过与算子f⁽ⁱ⁾(z0，zm)相关联，每个道都将会基于算子中的走时而时移(或校正)。在将相应算子f⁽ⁱ⁾(z0，zm)应用于在第一检波器处的100,000个道之后，所得的100,000个道将加起来成为一个cfp道。类似地，在将相应算子f⁽ⁱ⁾(z0，zm)应用于在第二检波器处的100,000个道之后，所得的道加起来成为第二个cfp道。在将算子应用于在1,000个检波器处的道之后，产生1,000个cfp道。在一些实施方式中，多个cfp道可以形成一个cfp道集，例如，所述1000个cfp道形成一个cfp道集。在一些实施方式中，一个cfp道为一个cfp道集。

在框208处，可以通过将cfp道集与算子f⁽ⁱ⁾(z0，zm)关联来产生dts道集。每个cfp道集都可以与相应的算子f⁽ⁱ⁾(z0，zm)关联，其中z0为与特定cfp道集对应的检波器位置。在先前的储层实例中，可以基于1,000个cfp道集产生1,000个dts道集。

在框210处，可以基于dts道集确定双向剩余走时δt⁽ⁱ⁾。例如，可以由dts道集鉴别反射事件，并且可以基于反射事件与时间零点的时间偏差来确定双向剩余走时δt⁽ⁱ⁾。对于先前的储层实例，图3a(之后将会更多讨论)示出了与1,000个索引为1至1000的检波器对应的1,000个dts道集。红线302示出了反射事件。对于每个dts道集，可以确定δt⁽ⁱ⁾，其等于在红线302与时间零点之间的时间差，并且对于来自图3a的1,000个检波器位置，可以确定1000个δt⁽ⁱ⁾。在一些实施方式中，可以对由dts道集得到的双向剩余走时进行插值以产生所有采集表面位置的整组剩余走时。例如，可以对由图3a得到的1000个δt⁽ⁱ⁾进行插值以产生对于101,000个采集表面位置的101,000个δt⁽ⁱ⁾。可以使用插值法如德洛内三角剖分法(delaunaytriangulationmethod)。在一些情况下，也可以基于震源位置对地震数据进行分类，并且可以对分类的数据进行步骤206-210以产生与100,000个源位置对应的100,000个δt⁽ⁱ⁾。

在框212处，确定来自dts道集的δt⁽ⁱ⁾是否为零或是接近于零(例如，在与零相差的预定阈值内)。换言之，框202确定dts道集是否包括在时间零点的平坦反射事件。在时间零点的平坦反射事件表明已经得到最终聚焦算子并且可以停止迭代，其中方法200继续进行框218。如果dts道集不包括在时间零点的平坦事件，则方法200继续进行框214以迭代地更新算子。

在框214处，对于每个采集表面位置，可以基于等式(3)和(4)(或等式(1)和(2))确定算子f⁽ⁱ⁺¹⁾(z0，zm)。例如，使用与1,000个检波器对应的1000个δt⁽ⁱ⁾更新对于1,000个检波器位置的cfp算子，并且使用与100,000个震源位置对应的100,000个δt⁽ⁱ⁾更新对于100,000个源位置的cfp算子。

在框216处，迭代计数i增加一，并且方法200返回到框206并将更新的算子应用于通过检波器记录的原始地震数据以产生cfp和dts道集。

在框218处，在已经得到最终聚焦算子之后(在dts道集显示出平坦反射事件时)，可以通过将来自最终算子的dts道集加起来来产生dts叠加。在先前的储层实例中，将1000个dts道集加起来以产生一个dts叠加。

在框220处，可以将目标反射面zm的反射系数确定为dts叠加在时间零点处的值。

在将方法200应用于各个反射面zm之后，可以由地震数据集产生地震图像，其中地震图像包括储层中的所有反射面的反射系数。在一些实施方式中，所产生的地震图像可以为3d图像。

返回到图1，在框108处，可以比较与不同地震数据集对应的反射面的反射系数(或地震图像)以确定在所述时间期内的反射系数的变化。例如，每月采集的地震数据可以产生地震图像。可以比较在18个月期间的18个地震图像以监测储层的变化，如储层中的烃液体的组成或位置变化。

图3a-5b示出了使用所描述的方法由在将co2注入到储层中的项目中采集的时移现场地震数据集而产生高效且准确的3d目标导向地震图像。在所述项目中，每个月以密集的发射分布(在x和y两个方向上10m间隔)和在地表下70m深度的约1000个埋入检波器进行完整的3d地震勘探。埋入检波器以减轻近地表复杂性的影响并且提高勘探之间的可重复性。采用互易性来产生cfp和dts道集，由所述道集得出聚焦算子。将在co2注入之前采集的基线勘探用于使用cfp和dts道集产生聚焦算子。将所描述的方法应用于时移数据集揭示了在储层水平处的变化，其与co2注入史相符。

图3a示出了根据一些实施方式的使用对于一个目标反射面的初始聚焦算子的一组dts道集300a。在图3a中，水平轴表示dts道集索引，并且垂直轴表示以秒计的时间。如先前所讨论的，该组的dts道集300a通过使用由nmo速度得出的初始聚焦算子f⁽⁰⁾(z0，zm)(在图4a中示出并且之后将讨论)包括约1,000个dts道集。dts道集的组300a包括由红线302示出的自动拾取的反射事件。由于所述项目中的最大偏移为3km，所以来自初始聚焦算子的dts道集表现为包括在时间零点处的几乎平坦的事件，而剩余时间变化仍可以由反射事件302计算并且用于更新聚焦算子。

图3b示出了根据一些实施方式的使用对于一个目标反射面的更新的聚焦算子的一组dts道集300b。在图3b中，水平轴表示dts道集索引，并且垂直轴表示以秒计的时间。该组的dts道集300b通过使用在更新初始聚焦算子达一次迭代后的更新的算子f⁽¹⁾(z0，zm)(在图4b中示出并且之后将讨论)包括约1,000个dts道集。该组的dts道集300b包括在时间零点处为平坦事件的反射事件304。

图4a示出了根据一些实施方式的对于一个目标反射面的初始聚焦算子400a。图4a示出了在目标反射面和所有采集表面位置之间的初始聚焦算子。水平轴和垂直轴分别表示在采集表面的x方向和y方向上的网格索引。绿点402表示目标反射面的位置。彩条408描绘在彩色光谱的颜色和以秒计的连续走时范围之间的映射。例如，红色区域如406表示走时为约1.4秒的聚焦算子，而深蓝色区域如404表示走时为约0.6秒的聚焦算子。类似地，图4b示出了根据一些实施方式的对于一个目标反射面在一次迭代后的更新的聚焦算子400b。图4c示出了根据一些实施方式的在对于一个目标反射面的初始聚焦算子400a和更新的聚焦算子400b之间的差400c。彩条410描绘在彩色光谱的颜色和以秒计的连续走时范围之间的映射。例如，红色区域如414表示大约为零的走时差，而深蓝色区域如416表示大约150毫秒的走时差。蓝点412表示埋入检波器的位置。由于图4c中的蓝点所描绘的埋入检波器的稀疏性，所以使用德洛内三角剖分法作为插值和外推工具来产生所有采集表面位置的更新的剩余走时。

图5a示出了根据一些实施方式的在目标水平面的反射系数的归一化均方根(nrms)振幅500a。在图5a中，水平轴和垂直轴分别表示在目标水平面的x方向和y方向上的网格索引。图5a示出了来自在第二个月采集的勘探2的反射系数的nrms振幅。黑色方块502、504、506和508表示四个co2注入位置。彩条510描绘在彩色光谱的颜色和连续nrms振幅值范围之间的映射。例如，红色区域如514表示值接近于一的nrms振幅，而深蓝色区域如512表示值接近于零的nrms振幅。假设不同勘探之间的上覆层速度相同，因此，对不同的勘探应用由勘探2(其因为在勘探2和后续勘探中得到的可重复性较好而用作基线勘探，这可能是由于表面的部分固结)得出的相同聚焦算子以产生地震图像。

图5b示出了根据一些实施方式的在勘探2和勘探17之间的nrms振幅差500b。勘探17(未示出)在勘探2之后15个月时进行。彩条516描绘在彩色光谱的颜色和连续nrms振幅差范围之间的映射。例如，橙色区域如518表示值为10的nrms振幅差，而淡蓝色区域如520表示值为约-7的nrms振幅。红色圆圈522中的异常被认为与co2注入有关，其与使用储层模拟器产生的预测一致，但是更好地界定了co2羽流(plume)的边界。

图6是根据一些实施方式的示例计算机系统600的框图，所述示例计算机系统用于提供如本公开内容中所描述的与所描述的算法、方法、函数、过程、流程和程序相关的计算功能。所示出的计算机602旨在涵盖任何计算装置，如服务器、台式计算机、膝上型/笔记本电脑、无线数据端口、智能电话、个人数字助理(pda)、平板式计算装置或在这些装置内的一个或多个处理器，包括计算装置的物理或虚拟实例(或两者)。此外，计算机602可以包括：包括输入装置(如可以接受用户信息的键座、键盘或触控屏)和传输与计算机602的操作相关的信息(例如，在图形用户界面(gui)上传输数字数据、视觉或音频信息(或信息的组合))的输出装置的计算机。

计算机602可以起到作为客户端、网络组件、服务器、数据库的作用或这些作用的组合以进行在本公开内容中所描述的主题。所示出的计算机602与网络630可通信地连接。在一些实施方式中，计算机602的一个或多个组件可以被配置成在包括基于云计算的、本地的、全球的环境或这些环境的组合的环境内运行。

在高水平，计算机602是可操作用于接收、传输、处理、存储或管理与所描述主题相关的数据和信息的电子计算装置。根据一些实施方式，计算机602还可以包括应用服务器、电子邮件服务器、网络服务器、缓存服务器、流数据服务器或这些服务器的组合，或者与它们可通信地连接。

计算机602可以通过网络630从客户端应用(例如，在另一计算机602上执行的)接收请求，并且通过使用一个或多个适当的软件应用处理所接收的请求来响应所接收的请求。另外，请求也可以从内部用户(例如，从命令控制台)、外部或第三方、其他自动化应用以及任何其他适当的实体、个体、系统或计算机发送到计算机602。

计算机602的每个组件都可以使用系统总线603通信。在一些实施方式中，计算机602的任何或全部组件(硬件或软件两者，或者硬件和软件的组合)可以通过系统总线603使用应用程序化接口(api)612或服务层613(或者api612和服务层613的组合)彼此或者与界面604(或者两者的组合)相互联系。api612可以包括例行程序、数据结构和对象类别的规格。api612可以是独立于或依赖于计算机语言的，并且是指完整界面、单个函数或甚至一组api。服务层613向计算机602或与计算机602可通信地连接的其他组件(无论是否示出)提供软件服务。计算机602的功能可以是对于使用该服务层的所有服务消费者可获得的。软件服务如由服务层613提供的那些通过限定的界面提供可重复使用的、限定的功能。例如，界面可以是用java、c++或者以可扩展标记语言(xml)格式或格式的组合提供数据的计算语言的组合写成的软件。尽管作为计算机602的集成组件示出，但是备选的实施方式可以示出关于计算机602的其他组件或与计算机602可通信地连接的其他组件(无论是否示出)为独立组件的api612或服务层613。此外，在不背离本公开的范围的情况下，api612或服务层613的任何或全部部分可以实施为另一软件模块的子模块或次模块、企业应用或硬件模块。

计算机602包括界面604。尽管在图6中作为单个界面604示出，但是根据具体的需求、要求或计算机602的具体实施方式，可以使用两个以上界面604。界面604被计算机602用于与在分布式环境中连接至网络630的其他系统(无论是否示出)通信。通常，界面604包括以软件或硬件(或者软件和硬件的组合)编码的逻辑，并且可操作用于与网络630通信。更具体地，界面604可以包括支持与通信相关的一个或多个通信协议的软件，以使得网络630或界面的硬件可操作用于在所示的计算机602之内和外部通信物理信号。

计算机602包括处理器605。尽管在图6中作为单个处理器605示出，但是根据具体的需求、要求或计算机602的具体实施方式，可以使用两个以上处理器。通常，处理器605执行指令并且操纵数据以进行计算机602的操作和本公开内容中所描述的任何算法、方法、函数、过程、流程和程序。

计算机602还包括数据库606，其可以为计算机602或可以连接至网络630的其他组件(无论是否示出)(或两者的组合)保存数据。例如，数据库606可以为存储符合本公开内容的数据的内存数据库或传统数据库。在一些实施方式中，根据具体的需求、要求或计算机602的具体实施方式以及所描述的功能，数据库606可以为两种以上不同的数据库类型的组合(例如，混合型内存和传统数据库)。尽管在图6中作为单个数据库606示出，但是根据具体的需求、要求或计算机602的具体实施方式以及所描述的功能，可以使用两个以上数据库(例如相同的或多种类型的组合)。尽管数据库606作为计算机602的整体组件示出，但是在备选实施方式中，数据库606可以在计算机602外部。例如，数据库606可以保存连续地震储层监测的地震数据集。

计算机602还包括存储器607，其可以为计算机602或可以连接至网络630的其他组件(无论是否示出)(或两者的组合)保存数据。例如，存储器607可以为存储符合本公开的数据的随机存取存储器(ram)、只读存储器(rom)、光学存储器、磁存储器等。在一些实施方式中，根据具体的需求、要求或计算机602的具体实施方式以及所描述的功能，存储器607可以为两种以上不同存储器类型的组合(例如，ram和磁存储的组合)。尽管在图6中作为单个存储器607示出，但是根据具体的需求、要求或计算机602的具体实施方式以及所描述的功能，可以使用两个以上存储器607(例如相同的或多种类型的组合)。尽管存储器607作为计算机602的整体组件示出，但是在备选实施方式中，存储器607可以在计算机602外部。

应用608为算法软件引擎，其根据具体的需求、要求或计算机602的具体实施方式、特别是相对于本公开内容中所描述的功能提供功能。例如，应用608可以作为一个或多个组件、模块或应用。此外，尽管作为单个应用608示出，但是应用608可以实施为在计算机602上的多个应用608。另外，尽管与计算机602一体地示出，但是在备选实施方式中，应用608可以在计算机602外部。

可以存在任意数量的与包含计算机602的计算机系统相关联或者在其外部的计算机602，各个计算机602通过网络630通信。此外，在不背离本公开的范围的情况下，术语“客户端”、“用户”和其他适当的术语可以视情况可互换地使用。此外，本公开考虑了多个用户可以使用一个计算机602，或者一个用户可以使用多个计算机602。

所描述主题的实施方式可以单独地或组合地包括一个或多个特征。

例如，在第一个实施方式中，一种方法包括在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应。所述储层包括多个反射面，其中每个反射面具有反射系数。对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子来确定所述多个反射面的反射系数。比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化。显示所述反射系数的所述变化。

前述和其他所述实施方式可以各自任选地包括以下特征中的一个或多个：

可与以下特征中的任一个组合的第一特征，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。

可与前述或以下特征中的任一个组合的第二特征，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子，

f⁽ⁱ⁺¹⁾(z0，zm)＝f⁽ⁱ⁾(z0，zm)δf⁽ⁱ⁾(z0，zm)

和

其中δt⁽ⁱ⁾是基于在第i次迭代中的差异时移(dts)道集所确定的剩余传播时间。

可与前述或以下特征中的任一个组合的第三特征，所述方法还包括：对于所述特定反射面zm和所述特定采集表面位置z0，基于在所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。

可与前述或以下特征中的任一个组合的第四特征，所述方法还包括：对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。

可与前述或以下特征中的任一个组合的第五特征，其中所述多个采集平面位置包括震源位置和地震检波器位置。

可与前述或以下特征中的任一个组合的第六特征，其中所述储层包含烃液体。

在第二实施方式中，一种系统包括：计算机存储器，和与所述计算机存储器可互操作地连接的一个或多个硬件处理器。所述一个或多个硬件处理器被配置成进行包括以下各项的操作：在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应。所述储层包括多个反射面，其中每个反射面具有反射系数。对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子来确定所述多个反射面的反射系数。比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化。显示所述反射系数的所述变化。

前述和其他所述实施方式可以各自任选地包括以下特征中的一个或多个：

可与以下特征中的任一个组合的第一特征，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。

可与前述或以下特征中的任一个组合的第二特征，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子

f⁽ⁱ⁺¹⁾(z0，zm)＝f⁽ⁱ⁾(z0，zm)δf⁽ⁱ⁾(z0，zm)

和

其中δt⁽ⁱ⁾是基于在第i次迭代时的差异时移(dts)道集所确定的剩余传播时间。

可与前述或以下特征中的任一个组合的第三特征，其中所述操作还包括：对于所述特定反射面zm和所述特定采集表面位置z0，基于所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。

可与前述或以下特征中的任一个组合的第四特征，其中所述操作还包括：对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。

可与前述或以下特征中的任一个组合的第五特征，其中所述多个采集平面位置包括震源位置和地震检波器位置。

可与前述或以下特征中的任一个组合的第六特征，其中所述储层包含烃液体。

在第三实施方式中，非暂时性的计算机可读介质存储能够由计算机系统执行的用于进行包括以下各项的操作的一个或多个指令：在一个时间期内接收与储层相关的多个地震数据集，其中所述多个地震数据集与在所述时间期内的不同时间接收的地震数据对应。所述储层包括多个反射面，其中每个反射面具有反射系数。对于所述多个地震数据集中的每一个，通过迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子来确定所述多个反射面的反射系数。比较与不同地震数据集对应的反射系数以确定在所述时间期内所述反射系数的变化。显示所述反射系数的所述变化。

前述和其他所述实施方式可以各自任选地包括以下特征中的一个或多个：

可与以下特征中的任一个组合的第一特征，其中每个cfp算子表示在特定反射面和特定采集表面位置之间的传播时间。

可与前述或以下特征中的任一个组合的第二特征，其中对于特定反射面zm和特定采集表面位置z0，迭代更新与所述多个反射面和多个采集表面位置相关的cfp算子包括使用以下式计算在第(i+1)次迭代时更新的cfp算子，

f⁽ⁱ⁺¹⁾(z0，zm)＝f⁽ⁱ⁾(z0，zm)δf⁽ⁱ⁾(z0，zm)

和

其中δt⁽ⁱ⁾是基于在第i次迭代是的差异时移(dts)道集所确定的剩余传播时间。

可与前述或以下特征中的任一个组合的第三特征，其中所述操作还包括：对于所述特定反射面zm和所述特定采集表面位置z0，基于所述特定反射面zm和特定采集表面位置z0之间的传播时间的初始估值来确定初始cfp算子。

可与前述或以下特征中的任一个组合的第四特征，其中所述操作还包括：对于所述多个地震数据集中的每一个，基于所确定的所述多个反射面的反射系数来产生所述储层的三维地震图像。

可与前述或以下特征中的任一个组合的第五特征，其中所述多个采集平面位置包括震源位置和地震检波器位置。

本说明书中所描述的主题和函数运算的实施方式可以以数字电子电路、有形地体现的计算机软件或固件、计算机硬件(包括本说明书中所公开的结构及其结构等同体)或它们中的一个或多个的组合来实施。本说明书中所描述的主题的实施方式可以实施为一个或多个计算机程序，即在有形的、非暂时性的计算机可读的计算机存储介质上编码的计算机程序指令的一个或多个模块，其由数据处理设备执行或者控制数据处理设备的操作。备选地或另外地，程序指令可以在人工生成的传播信号例如机器生成的电信号、光信号或电磁信号中/上编码，所述机器生成的电信号、光信号或电磁信号产生用于编码用于传输至合适接收器设备以由数据处理设备执行的信息。计算机存储介质可以为机器可读的存储装置、机器可读的存储基底、随机或串行存取存储装置或计算机存储介质的组合。

术语“实时(real-time)”、“即时(realtime)”、“同时(realtime)”、“(快速)实时(rft)”、“近实时(nrt)”、“准实时”或类似术语(如本领域普通技术人员理解的)意指动作或响应是时间接近的，使得个体感觉所述动作和响应基本上同时发生。例如，在个体访问数据的动作后显示(或开始显示)数据的响应的时间差可以小于1ms，小于1s，或小于5s。尽管不必即刻显示(或开始显示)所请求的数据，但是考虑到所描述的计算系统的处理限制和例如采集、准确测量、分析、处理、存储或传输所述数据所需的时间，在没有任何有意延迟的情况下显示(或开始显示)数据。

术语“数据处理设备”、“计算机”或“电子计算装置”(或如本领域普通技术人员所理解的等同体)是指数据处理硬件并且涵盖所有种类的用于处理数据的设备、装置和机器，包括例如可编程处理器、计算机或者多个处理器或计算机。所述设备也可以是或者还可以包括专用逻辑电路，例如中央处理器(cpu)、fpga(现场可编程门阵列)或asic(专用集成电路)。在一些实施方式中，数据处理设备或专用逻辑电路(或者数据处理设备或专用逻辑电路的组合)可以基于硬件或基于软件的(或者是基于硬件和基于软件两者的组合)。所述设备可以任选地包括产生用于计算机程序的执行环境的代码，例如，构成处理器固件、协议栈、数据库管理系统、操作系统或执行环境的组合的代码。本公开考虑在有或没有传统操作系统的情况下使用数据处理设备，所述传统操作系统例如为linux、unix、windows、macos、android、ios或操作系统的组合。

也可以称为或者描述为程序、软件、软件应用、模块、软件模块、脚本或代码的计算机程序可以任何形式的编程语言写成，包括编译语言或解释语言或者说明性语言或过程语言。计算机程序可以以任何形式部署，包括作为独立程序、作为模块、组件、子例程或适合在计算环境中使用的其他单元。计算机程序可以但不必与文件系统中的文件对应。程序可以存储在文件的一部分中，所述文件保存其他程序或数据，例如，存储在标记语言文档中、在专用于所讨论程序的单个文件中或在多个协同文件(例如，存储一个或多个模块、子程序或部分代码的文件)中的一个或多个脚本。计算机程序可以部署为在一个计算机上或在位于一个地点或者分布在多个地点并通过通信网络相互连接的多个计算机上执行。尽管在多个图中示出的程序的部分显示为通过多个目的、方法或其他过程实施多个特征和功能的个体模块，但是相反地，程序可以视情况包括多个子模块、第三方服务、组件、库等。相反地，可以视情况将多个组件的特征和功能合并为单个组件。用于进行计算确定的阈值可以是静态地、动态到或同时静态地和动态地确定的。

本说明书中所描述的方法、过程或逻辑流程可以由一个或多个可编程计算机进行，所述一个或多个可编程计算机执行一个或多个计算机程序以通过对输入数据进行运算并且产生输出来发挥功能。所述方法、过程或逻辑流程也可以由专用逻辑电路，例如，cpu、fpga或asic进行，并且设备也可以实施为专用逻辑电路，例如，cpu、fpga或asic。

适用于执行计算机程序的计算机可以基于通用或专用的微处理器。通常，cpu将从只读存储器(rom)或随机存取存储器(ram)或两者接收指令和数据。计算机的基本元件为用于进行或执行指令的cpu和一个或多个用于存储指令和数据的存储装置。通常，计算机将还包括一个或多个用于存储数据的大容量存储装置(例如，磁盘、磁光盘或光盘)，或者操作连接为从所述大容量存储装置接收数据或将数据传输至所述大容量存储装置，或两者。然而，计算机不必具有这样的装置。此外，计算机可以嵌入另一装置(例如，举例来说，移动电话、个人数字助理(pda)、移动音频或视频播放器、游戏控制器、全球定位系统(gps)接收器或便携式存储装置例如通用串行总线(usb)闪存驱动器)中。

适用于存储计算机程序指令和数据的计算机可读介质(暂时性或非暂时性的，视情况而定)包括所有形式的非易失性存储器、介质和存储装置，包括例如半导体存储装置(例如，可擦除编程只读存储器(eprom)、电可擦除编程只读存储器(eeprom)和闪存装置)、磁盘(例如，内置硬盘或可移动盘)、磁光盘和光学存储装置(例如，cd-rom、dvd+/-r、dvd-ram和dvd-rom盘)。存储器可以存储多种对象或数据，包括高速缓存、类、框架、应用、备份数据、工作、网页、网页模板、数据库表、存储动态信息的储存库和任何其他适当的信息，包括任何参数、变量、算法、指令、规则、约束或其参考。另外，存储器可以包括任何其他适当的数据，如日志、策略、安全性或访问数据或者报告文件。处理器和存储器可以通过专用逻辑电路进行补充，或者并入专用逻辑电路中。

为了提供与用户的交互，可以在具有用于向用户显示信息的显示装置(例如crt(阴极射线管)、lcd(液晶显示器)、led(发光二极管)或等离子体监视器)以及键盘和点击装置(例如，鼠标、轨迹球或轨迹板)的计算机上实施本说明书中所描述的主题的实施方式，用户可以通过所述键盘和点击装置向计算机提供输入。也可以使用触控屏如具有压力敏感性的平板计算机表面或者利用电容感应或电感应的多点触控屏向计算机提供输入。也可以使用其他种类的装置来提供与用户的交互；例如，向用户提供的反馈可以为任何形式的感觉反馈，例如，视觉反馈、听觉反馈或触觉反馈；并且来自用户的输入可以以任何形式接收，包括声学输入、语音输入或触觉输入。另外，计算机可以通过将文档发送至用户所使用的装置和从用户所使用的装置接收文档与用户交互；例如，通过响应于从网络浏览器接收的请求将网页发送至在用户的客户端装置上的网络浏览器。

术语“图形用户界面”或“gui”可以以单数形式或复数形式使用以描述一个或多个图形用户界面以及特定图形用户界面的显示中的每一个。因此，gui可以表示任何图形用户界面，包括但不限于网络浏览器、触控屏或命令行界面(cli)，其处理信息并且高效地向用户提供信息结果。通常，gui可以包括多个用户界面(ui)要素，其中的一些或全部与网络浏览器相关，如交互域、下拉列表和按钮。这些和其他ui要素可以与网络浏览器的功能有关或者表现网络浏览器的功能。

可以在计算系统中实施本说明书中所描述的主题的实施方式，所述计算系统包括后端组件(例如，数据服务器)，或者包括中间组件(例如，应用服务器)，或者包括前端段组件(例如，具有图形用户界面或网络浏览器的客户端计算机，用户可以通过其与本说明书中所描述的主题的实施方式交互)，或者一个或多个这样的后端、中间或前端组件的任意组合。系统的组件可以通过有线或无线数字数据通信(或数据通信的组合)的任何形式或介质例如通信网络相互连接。通信网络的实例包括局域网(lan)、无线电接入网(ran)、城域网(man)、广域网(wan)、全球微波接入互操作性(wimax)、使用例如802.11a/b/g/n或802.20(或802.11x和802.20的组合或其他符合本公开的协议)的无线局域网(wlan)、互联网的全部或一部分、或者通信网络的组合。网络可以在网络地址之间通信数据，例如，互联网协议(ip)数据包、帧中继框架、异步传输模式(atm)信元、语音或视频。

计算系统可以包括客户端和服务器。客户端和服务器通常远离彼此，并且一般通过通信网络交互。客户端和服务器的关系由在相应计算机上运行并且彼此具有客户端-服务器关系的计算机程序产生。

尽管本说明书含有许多具体的实施方式细节，但是这些不应被解释为对可以要求保护的范围的限制，而应解释为可以特定于具体构思的具体实施方式的特征的描述。在本说明中描述的在不同实施方式情况下的某些特征也可以组合地在单个实施方式中实施。相反地，在单个实施方式情况下描述的多个特征也可以在多个实施方式中单独地实施，或以任何合适的子组合实施。此外，尽管前述特征可以描述为以特定组合起作用并且甚至最初这样要求，但是，在一些情况下，来自所要求的组合的一个或多个特征可以从组合中省去，并且所要求的组合可以涉及子组合或子组合的变化方案。

已经描述了本主题的具体实施方式。对本领域技术人员来说明显的，所述实施方式的其他实施方式、变化方案和置换方案在所附权利要求的范围内。尽管操作在附图中进行描绘或者要求以特定顺序，但是这不应被理解为要求这类操作以所示的特定顺序或以先后顺序进行，或者进行所有示出的操作(一些操作可以认为是任选的)以实现期望的结果。在某些情形中，多任务处理或并行处理(或多任务处理和并行处理的组合)可以是有利的，并且视情况进行。

此外，前述实施方式中的多个系统模块和组件的分离或集成不应理解为在所有实施方式中都要求这样的分离或集成。应理解，所述程序组件和系统通常可以在单个软件产品中集成到一起或者封装为多个软件产品。

因此，前述示例性实施方式不限定或约束本公开。在不脱离本公开的精神和范围的情况下，其他改变、替换和变化也是可能的。

此外，任何要求保护的实施方式被认为是至少适用于计算机实施的方法、非暂时性的计算机可读介质和计算机系统来实施，所述非暂时性的计算机可读介质存储用于进行所述计算机实施的方法的计算机可读指令，所述计算机系统包括与硬件处理器能够互相操作地连接的计算机存储器，所述硬件处理器被配置成进行所述计算机实施的方法或存储在所述非暂时性的计算机可读介质上的指令。