使用辐射方向图校正重新标定虚拟源的制作方法

优先权要求

本申请要求于2017年1月27日提交的美国专利申请no.62/451,297的优先权，其全部内容通过引用合并于此。

背景技术：

虚拟源(vs)重新标定是一种基于干涉的采集和处理方法，可以克服时移地震图像应用中的近地表复杂性。该技术将下行波能量与上行能量互相关，以将地表源记录重新标定到埋藏式接收器位置。通过将这两个路径的波场互相关，vs记录可以帮助解决在地表源和埋藏式接收器之间发生的时移变化。但是，存在重复性问题，尤其是在具有复杂的近地表结构的近场土地应用中。

技术实现要素：

本公开内容描述了如下方法和系统，包括用于虚拟源(vs)重新标定的计算机实现的方法、计算机程序产品以及计算机系统。

在一个实现方式中，将接收到的炮集分类到一个公共接收器集。根据合成数据计算出地震波场直接到达的目标三维(3d)幅度谱。针对每个接收器计算场数据中的地震波场直接到达的3d幅度谱。从场数据的3d幅度谱到目标响应计算出匹配滤波器，并将其应用于与公共接收器集分离的下行地震波场，以生成滤波后的下行地震波场。对滤波后的下行地震波场进行时间相关的平滑处理，以生成平滑的下行地震波场。计算与公共接收器集分离的上行地震波场和平滑的下行地震波场之间的互相关。

可以使用计算机实现的方法来实现前述的实现方式；存储用于执行计算机实施的方法的计算机可读指令的非暂时性计算机可读介质；以及计算机实现的系统，该系统包括与硬件处理器可互操作地耦接的计算机存储器，所述硬件处理器被配置为执行该计算机实现的方法/存储在该非暂时性计算机可读介质上的指令。

在本说明书中描述的主题可以实现在特定实施方案中，以便实现以下优点中的一个或多个。首先，描述了下行近似过程，该过程将时间门控和地表源阵列渐缩技术与时间相关平滑以及查找/应用针对幅度谱的匹配滤波器的技术相结合，从而近似出理想的p波下行波场，以及其多维辐射方向图。描述了迭代步骤，这些步骤用于查找匹配的滤波器并将其应用于3d空间中的下行直接到达的理想幅度谱以及与时间相关的平滑处理。通过在与上行反射相关之前、将下行直接p波的最佳近似作为目标，获得与较干净的下行p波场相关联的正确图像的显著增强，同时抑制与串扰相关的伪像。其次，所描述的方法解决了冲突的可重复性问题，例如下行直接到达内的波干扰，与上行反射到达的频率含量不匹配以及场数据中存在的时移勘测内部和之间的变化的地表源谱。第三，与单独的时间门控和阵列渐缩相比，所描述的方法解决了土地监测员勘测中不可重复的问题。第四，所描述的方法通过减少近地表变化的影响来提高可重复性。第五，如果适用埋藏式传感器，则所描述的方法可以提高二维(2d)/3d地震图像的质量。其它优点对于本领域普通技术人员将是明显的。

本说明书的主题的一个或多个实施方式的细节在附图和具体实施方式中阐述。通过说明书、附图和权利要求书，所述主题的其它特征、方面和优点将变得明显。

附图说明

本专利或申请文件包含至少一个彩色附图。含有彩色附图的本专利申请公开的拷贝将在提出请求并支付必要的费用后由专利和商标局提供。

图1a和图1b是示出根据本公开的实施方式的、用于四维地震储层监控的重新标定虚拟源(vs)的方法的框图，该方法改进了图像质量和可重复性。

图2a是示出根据本公开的实施方式的汉明窗的使用的数据图。

图2b是根据本公开的实施方式的表示余弦时间锥形的数据图。

图3是示出根据本公开的实施方式的特定直接到达地震迹线的时间门控的数据图。

图4a至图4e是示出根据本公开的实施方式的、使用图1a和图1b中描述的方法处理典型源阵列产生的地震数据及其频率-波数功率谱的四个阶段的的图像。

图5a和图5b是分别示出根据本公开的实施方式的示例理想信号和噪声信号的数据图。

图6a至图6d是示出根据本公开的实施方式绘制的地震勘测数据的公共深度点(cdp)堆栈的图像。

图7a至图7d是示出了根据本公开的实施方式的、绘制的从1120(毫秒(ms))至1300ms的十三个样本地震勘测中每一个的共同深度点(cdp)迹线的图像。

图8是示出根据本公开的实施方式的、与每种方法相关联的归一化均方根(nrms)值的直方图的图。

图9是示出根据本公开的实施方式的、通过计算相邻的cdp之间的nrms(nrmsc)来比较下行直接p波近似之间的图像质量的图。

图10是根据本公开的实施方式的、用于提供与如本公开内容中描述的算法、方法、功能、处理、流程和过程相关联的计算功能的示例性计算机系统的框图。

在各个附图中，相同的附图标记和标记表示相同的元件。

具体实施方式

以下详细描述描述了虚拟源(vs)重新标定，并且被呈现为使得本领域技术人员能够在一个或多个具体实施方式的上下文中做出和使用所公开的主题。对公开的实施方式的各种修改、改变和排列可以实现并且对本领域技术人员而言将显而易见，并且在不背离本公开的范围的情况下，所定义的一般原理可适用于其他实施方式和应用。因此，本公开并非意在限于所描述的或示出的实施方式，而应赋予与所公开的原理和特征一致的最宽范围。

vs是一种基于干涉的采集和处理方法，可以克服时移(time-lapse)地震图像应用中的近地表结构复杂性。vs方法将下行地震波能量与上行(反射的)地震波能量互相关，以将地表源记录重新标定(redatum)到埋藏式接收器位置。通过将这两个路径的波场互相关，vs记录可以帮助解决在地表地震源和埋藏式地震接收器之间发生的时移变化。换句话说，vs重新标定用于简化记录的地震波场，并消除与位于地震源和地震接收器之间的异质性相关的失真。但是，存在重复性问题，尤其是在具有复杂的近地表结构的近场土地应用中。

vs重新标定的另一目标是通过校正时移噪声来改善时移勘测的可重复性(例如近地表昼夜和季节周期以及采集几何体和炮耦合的微小变化)。在较小的时间窗口内使用总的早期到达地震波场来估计下行直接p波已被证实是vs发展的必不可少的部分。作为时间窗的一种替代方法，也可以通过将vs记录与下行早期到达地震波场的估计幅度辐射方向图去卷积，或者通过在满足记录相位的情况下用所需的方向图替换该方向图，来改善vs数据。

迄今为止，大多数场应用都专注于在深海环境或水平井中使用地震接收器采集的vs数据。关于陆上vs场数据的报道甚少，因为在这些地区，埋藏式地震接收器靠近活跃(近场)地震源。利用vs重新标定，近场地震源可能会导致下行的早期到达地震波场与上行的地震反射之间的带宽不匹配，从而导致输出vs记录出现高频偏差。在下行/上行波窗口中发生的多种波类型之间的过度干扰也会使vs质量下降。

描述了一种vs重新标定方法，该方法包括匹配滤波器的应用和迭代估计，用于在三维(3d)空间中的下行直接地震波能量到达的理想幅度谱以及在互相关之前应用时间相关的平滑处理。时频波数滤波是基于下行、直接的p波估计和近场近似进行的，包括源阵列空间渐缩和相位不变近似。幅度谱替换提高了vs重新标定的可重复性和图像质量。所描述的方法使用干涉测量法的概念将利用地表上的地震源和地表下的地震接收器采集的3d地震数据集转换为看起来好像源和接收器都位于复杂的近地表结构下方的数据集。

所描述的技术获得了显著改善的vs响应，并且可以克服前面提到的出现的可重复性问题。特别是在具有复杂的近地表结构的近场土地应用中。所描述的方法还保留了已记录的下行早期到达的地震波场的原始相位，但控制了vs的辐射方向图，使其紧密代表理想的下行直接p波的辐射方向图。

使用从场数据中提取的可重复性指标可以量化性能改进。与以前的方法(例如vs时间门控技术)相比，结果表明，新描述的匹配滤波器可以显著提高从所得的时移vs地震数据产生的可重复性和图像质量。

在较高的层次上，vs方法将地表地震源重新标定到埋藏式地震传感器位置，并旨在消除或减轻近地表复杂性的影响。由所有参与地震源贡献的弹性vs互相关双向波场如下面的等式(1)中所述：

其中×表示时间互相关，ra、rb和rs分别表示两个地震接收器在空间位置a和b处的空间坐标，以及震源位置。当ra被视为vs时，所得的v(rb|ra)是将在接收器rb中记录的干涉数据。等式(1)右侧总和中的左列表示由浅埋式地震传感器在接近偏移处接收的下行直接波场，而右列表示上行反射的地震波场。dp、up,、dm、um和us是接收的分别与直接p波到达、倍数波和其他剪切波有关的地震波场。

在等式(1)，只有下行的直接p波dp(ra|rs；t)与上行up(rb|rs；t)的第一相关才形成正确的地震图像。其余的相关性会生成称为“串扰”的伪像。更一般地，在属于同一波模式的ra和rb处的地震波场的相关性将产生正确的事件，而其他分组可能由于相位不正确而不会。

所描述的方法引入了两种重叠的方法，该方法能够减轻与串扰相关的地震图像伪像，同时增强下行地震波场窗口内的p波，从而提高vs重新标定的可重复性和图像质量。所描述的方法包括迭代步骤，该迭代步骤用于查找匹配滤波器并将其应用于3d中下行直接到达的理想幅度谱，以及时间相关的平滑处理。为了恢复vs的理想辐射方向图，所描述的方法实现了具有时间相关的平滑处理的匹配滤波器。通过在与上行地震反射相关之前、将下行直接p波的最佳近似作为目标，可以增强与较干净的下行p波场相关联的正确图像，同时抑制与串扰相关联的伪影。所描述的与辐射方向图校正和相关幅度控制互相关的方法可潜在地用于信号处理，例如电气工程、医学成像和测井降噪。

在较高层次上，两种重叠方法用于：

a)如图1a所示，根据公共接收器集计算下行地震波场直接到达的3d频率-波数谱。通过替换近地表层来估计直接p波的理想3d频率-波数谱。迭代求解匹配滤波器，以最小化(仅)在计算的和理想的直接到达之间的3d幅度谱的失配函数(如图1b所示)。该速度用于计算理想的直接p波谱，该谱在每次迭代中都用残差更新。

b)如图1b所示，应用在下行地震波场窗口内获得的3d匹配滤波器。原始地震迹线的时间门控捕获与直接地震波场到达相关联的地震信号。对地震炮集的输入源阵列进行空间相关的滤波，目的是有利地权衡偏移范围内的输入源迹线的贡献，以使远场近似有效且固定相的干涉法假设近似成立。时间相关的平滑处理在窗口门控内执行。每次迭代后都可以执行质量控制。

方法b)的滤波器是通过根据等式(2)在vs位置定义以下二维(2d)高斯加权函数中心来构造的：

其中，r是源阵列孔径，r是源埋藏式接收器对之间的偏移，而g(r)是输出高斯权重。

类似地，时间相关的滤波器可以表示为从直接到达小波的峰值开始的汉明窗。参考图2a，图2a是示出根据本公开的实施方式的汉明窗的使用的数据图200a。在图2a中，示出了没有时间锥形(ttaper)202以及具有时间锥形204的直接地震到达。

如图2b所示，图2b是根据本公开的实施方式的表示余弦时间锥形200b的数据图。ttaper1(毫秒(ms))206是汉明窗渐缩至直接到达的开始时间(缺省为0ms)。ttaper2(ms)208是汉明窗渐缩至直接到达的结束时间(缺省为150ms)。

根据等式(3)：

t是时间门控长度，t代表从小波峰开始计数的时间样本。

在方法a)中，针对下行直接地震波场到达的3d幅度谱，迭代地查找并应用理想匹配滤波器可被视为频率-波数(f-k)滤波，从而实现与所需的下行p波场相关联的理想辐射方向图，同时保持相位不变。具体而言，可以将理想匹配滤波器作为一个最小化问题来实现，如等式(4)所示：

其中d(kx,ky,ω)和dp(kx,ky,ω|f)分别是输入迹线和所需的p波场的3d幅度谱(波数kx、波数ky和频率ω)，f(kx,ky,ω)是对应的匹配滤波器，其旨在使用最小二乘反演过程中的最小能量准则来求解。所需的p波幅度谱dp(kx,ky,ω|f)基于初始均质速度模型，该模型使用与每次迭代中获得的匹配滤波器相关的残留误差进行更新。n和i是总迭代次数和迭代索引。

第一次迭代后，可以使用等式(5)表示的正则化反演形式迭代地更新滤波器：

fi+1＝fi+(d^td+εi)-¹d^t(dp-dfi)(5),

其中ε是正则化因子。需要迭代，因为dp也会作为f的函数进行更新。

通过选择接近偏移的首次到达来建立p波的初始均质速度模型，从而相应地计算幅度谱dp(kx,ky,ω|f)。常数均质p速度模型被迭代更新，其与匹配的滤波器更新耦合。使用相同的最小化成本函数(即等式(4))来实现速度，但是会受到速度的影响。

建立初始模型后，可以使用ricker小波的3d解析格林函数来迭代地更新速度，如等式(6)所示：

其中p是恒定的速度模型，δd是数据残差，r是偏移量，t是到达时间。ω表示ricker小波的中心频率。

为了应用从方法a)获得的匹配滤波器，数据被重新组织到公共炮(shot)域中-以每个源-接收器对的形式编写，如等式(7)所示：

其中v(ra|rb)是应用匹配滤波器后的输出迹线，当ra被视为虚拟源时记录在接收器rb中，sm是时间相关的平滑算子，并且u(rb|rs；ω)表示频域中的上行波场迹线。

自适应更新的匹配滤波器仅在与上行(反射)波场互相关之前应用于下行(直接p)波场。然后，在互相关之后，匹配的滤波器免于对反射波场携带的四维(4d)信号产生负面影响。

在地震接收器埋在浅地中的情况下(例如，比一个接收器波场短的30米(m)-150m)，典型的基于波方程或基于射线的速度更新方法不适用于如此小的规模。相反，均质速度模型可用于幅度谱校正，而不是使用相位。此外，使用解析形式求解3d格林函数比使用数值求解器更有效地进行计算。

图1a和图1b是示出根据本公开的实施方式的、用于4d地震储层监控的重新标定虚拟源(vs)的方法100的框图，该方法改进了图像质量和可重复性。为了清楚说明，下面的描述在总体上描述了在本说明书中的其他附图的上下文中的方法。然而，应理解，该方法可以例如是情况由适当的任意合适的系统、环境、软件和硬件、或系统、环境、软件和硬件的组合来执行。在一些实施方式中，该方法的各个步骤可以并行、组合、循环或以任意顺序运行。

在较高的层次上，图1a根据实际公共接收器集计算下行直接到达的3d频率-波数谱。通过替换近地表层来估计直接p波的理想3d频率波数谱。迭代地求解匹配滤波器，以最小化(仅)在计算出的和理想的直接地震波场到达之间的3d幅度谱的失配函数。在每次迭代中都会使用残留误差来更新用于计算理想直接p波谱的速度。注意，等式(4)代表整个方法。

在102处，接收可用的炮集。在典型的实施方式中，以反映标准地震仓(bin)格网的数据结构(例如，地震仓格网数据交换格式或其他格式的数据阵列)接收可用的炮集数据。方法100自102进行到104。

在104处，将接收到的炮集分类到公共接收器集中。在典型的实施方式中，任何公用排序算法都可以用于根据各种数据标准(例如，按时间、地理位置或其他数据)进行分类。方法100从104前进到106和116。

在一些实施方式中，可以利用并行类型的处理功能(例如，多线程)。例如，每个处理分支(即，以106和116开始)可以通过适当的功能来管理，以确保在特定处理、操作或功能之前(例如，方法100的118处)完成每个分支。

在106处，选择接近偏移的第一到达以提供初始恒定的p波速度模型。在典型的实施方式中：1)在很小的偏移量内(例如，<20m)选择直接的首次到达；2)自动选择攻击(onset)到达；3)将每个地震炮到埋藏式接收器之间的距离除以选定的到达时间，以沿每个炮-接收器波径生成恒定的速度值；以及4)来自3)的平均速度值用于输出被视为初始恒定的p波速度模型的统一均质速度。方法100自106进行到108。

在108处，使用等式(6)将恒定的p波速度模型从地表地震源迭代地更新到埋藏式地震接收器。注意，在第一次通过所描述的方法时，不更新恒定的p波速度模型。方法100自108进行到110。

在110，使用3d傅里叶变换将恒定的p波速度模型变换为频率-波数域。方法100自110进行到112。

在112处，根据合成数据计算直接到达的目标3d幅度谱。在典型的实施方式中，通过使用标准有限差分求解器以106的恒定p波速度模型对声波方程进行数值求解来生成合成数据。

为了在实施方式中生成合成数据集，可以通过下式给出提供矢量场粒子速度的示例声波方程：

其中是拉普拉斯算子，u是粒子速度，t是时间，c是声速。方法100自112进行到118。

在118处，使用等式(5)计算场数据的3d幅度谱到目标响应的匹配滤波器。在一些实施方式中，将计算出的匹配滤波器存储在计算机化的数据存储库中。方法100从118前进到图1b中的120。

在120，确定是否应当重新计算用于更新匹配滤波器的迭代过程。通常，确定是基于公差值和场数据的3d幅度谱(请参阅116)与合成数据的目标3d幅度谱(112的等式(8))的计算差值(失配)之间的比较。在典型情况下，公差值设置为目标幅度谱的平均值的大约1％。如果确定计算出的失配值小于或等于设定的公差值，则方法100前进至108以重新计算匹配的滤波器。场数据的3d幅度谱(来自116)和计算的差(116和112之间)返回到108。然而，如果确定计算出的失配值大于设定的公差值，则方法100进行到图1b中的122。

转到另一个处理分支，在114处，使用3d傅里叶变换将公共接收器集转换成频率-波数域。方法100自114进行到116。

在116处，使用等式(5)(例如，标准傅里叶变换)针对每个埋藏式接收器(vs)计算场数据中直接到达的3d幅度谱。如前所述，方法100从116进行到118。

在较高的层次，图1b在下行波窗口内应用获得的3d匹配滤波器。在窗口门控内执行时间相关的平滑处理。每次迭代后执行质量控制。

在122，接收图1a(来自104)的公共接收器集。方法100自122进行到124。

在124处，执行将公共接收器集分离成上行/下行地震波场的操作。请注意，如果地震检波器和水听器同时安装，则水听器只会记录标量压力响应，而不会区分上下波场。地震检波器记录埋藏式接收器位置的矢量位移，该位移对于上下波场是不同的。这些分量的有效组合(例如，自适应求和或减法)可以输出上行和下行波场数据。如果场中可以安装地震检波器和水听器，则需要上行/下行分离，以满足vs的请求前数据要求。否则，直接的早期到达被视为下行波，而目标反射数据则被视为上行波场。方法100从124前进到126和138。

在一些实施方式中，可以利用并行类型的处理功能(例如，多线程)。例如，每个处理分支(即，以126和138开始)可以由适当的功能来管理，以确保在特定处理、操作或功能之前(例如，在142处)完成每个分支。

在126处，在一些实施方式中，将分离的下行波场数据简单地传递以用于进一步处理(例如，到128)。在其他实施方式中，可以在使用之前(例如，归一化)在进一步处理之前对分离的下行波场数据进行预处理。

方法100自126进行到128。

在128处，对每个地震迹线进行时间门控。参考图3，图3是示出根据本公开的实施方式的特定直接到达地震迹线的时间门控的数据图300。在图3中，相对于直接到达迹线302，tda(ms)304是直接到达地震迹线302的迹线长度，门控306是直接到达的门控长度(缺省50ms)，并且fb时间308是地震波的开始时间。在数据图300的底部轴上指示直接到达迹线302的幅度310(图示单位为0.001m)。方法100自128进行到130。

在130处，使用等式(2)执行空间相关孔径滤波。方法100自130进行到132。

在132，从计算机化的数据存储库中加载根据图1a的迭代计算的匹配滤波器。方法100自132进行到134。

在134处，将加载的匹配滤波器应用于下行波场。方法100自134进行到136。

在136处，使用等式(3)对滤波后的下行波场执行时间相关的平滑处理。注意，可以在回到132的循环中执行质量控制(qc)和参数测试/选择。通常，qc是人为干预步骤，用于在每次匹配滤波器迭代之后检查场数据和合成数据的幅度谱之间的上述失配计算的收敛性。如果不收敛，则需要根据图1a重新计算匹配滤波器。在人员操作者确定需要重新计算匹配滤波器的情况下，人员操作者可以将处理引导回图1a的处理，以重新计算/更新匹配滤波器，该匹配滤波器可以在132处被重新加载以和方法100一起使用。在一些实施方式中，可以通过自动化过程或机器学习/人工智能过程来执行qc功能。方法100自136进行到142。

在138处，在一些实施方式中，将分离的上行波场数据简单地传递以用于进一步处理(例如，到140)。在其他实施方式中，可以在使用之前(例如，归一化)在进一步处理之前对分离的上行波场数据进行预处理。方法100自138进行到140。

在140处，在每个上行地震波场执行地滚波去除。在炮集(shotgather)中，地滚波的特点是晚到、高幅度和低频事件，这些事件定义了陡峭的三角形中心区域，该区域遮盖了反射到达。地滚波事件的视在速度较小(smallapparentvelocity)(或等效地，倾角较大)，可以通过例如执行2d傅立叶变换进行隔离和去除。在一个示例中，地滚波可位于扇形区域中。通过将扇形区域中的傅立叶变换值清零，然后进行傅立叶逆变换，可以去除地滚波。方法100使用来自136的平滑数据从140前进到142。

在142，计算每个时间相关的平滑上行地震波场和时间相关的平滑下行地震波场之间的互相关，使用等式(7)在每个反射和方向到达迹线之间执行互相关。请注意，从124、138、140和142执行循环遍历来自一个炮的所有接收器。方法100自142进行到144。

在144处，vs炮集数据被堆叠。请注意，从122、124、138、140、142和144执行循环遍历来自一个地震勘测的所有炮集。方法100自144进行到146。

在146处，开始堆叠的vs炮集的输出。例如，可以将堆叠的vs炮输出到视觉显示设备进行可视化，或者输出到计算机化数据存储库。在146之后，方法100停止。

图4a至图4e是根据本公开的实施方式的数据图400a-400e，其分别示出了使用图1a和图1b中描述的方法100来处理典型震源阵列生成的地震数据(2d–7.5m采样时的孔径30m)及其f-k功率谱(vs的辐射方向图)的四个阶段。参考图5a和图5b，图5a和图5b是根据本公开的实施方式的数据图500a和500b，其分别示出了示例理想信号和噪声信号。如从图5a中可以看出，信号迹线502a由平滑曲线表示。在图5b中，地震信号502b包含各种噪声(例如，噪声504b和506b)。

图4a是根据本公开的实施方式的以相应的时间(ms)/迹线和频率(hz)/波数x(1/m)示出输入的早期到达总地震波场的数据图400a。例如，并参考图5b，信号迹线502b包含各种噪声(例如，噪声504a和506a)。作为参考，时间/迹线数据图中的理想信号迹线应类似于图5a的信号迹线502a以及频率/波数数据图应类似于对称的菱形形状(如图4e中所示)。

图4b示出了根据本公开的实施方式的在与图4a相关联的操作之后的门控(窗口化)和空间相关的滤波之后的数据图400b。与从图4a的比较可以看出，在门控和空间相关的滤波操作之后，所有信号迹线中的总噪声已降低。

图4c示出了根据本公开的实施方式的在与图4b相关联的操作之后查找并应用匹配滤波器之后的数据图400c。从与图4b的比较可以看出，在查找并应用匹配滤波器之后，所有信号迹线中的总噪声已经进一步降低。

图4d示出了根据本公开的实施方式的在与图4c相关联的操作之后，将时间相关的平滑加上匹配滤波器应用于加窗和空间渐缩(taper)数据之后的数据图400d。从与图4c的比较可以看出，在应用了时间相关的平滑处理加上匹配滤波器之后，所有信号迹线中的总噪声都得到了进一步降低。匹配的滤波器对接收器集的频谱带宽进行归一化，并抑制尾波。时间相关的平滑处理进一步抑制了噪声并稳定了小波特征。

图4e示出了根据本公开的实施方式的在与图4d相关联的操作之后，在数值上求解用于真实的下行p波的解析解之后的数据图400e。从与图4d的比较可以看出，在数值上求解真正的下行p波的解析解之后，所有信号迹线中的总噪声已经进一步减小到无噪声(或基本上没有噪声)。在这里，通过计算波形及其f-k频谱，通过使用解析解用均质层替换覆盖层，可以量化估计的p波近似于真实的p波的质量。如图4e中所示，下行波场在时间-偏移域以及频率-波数域中都逐渐接近真实的p波。

图6a-6d是根据本公开的实施方式的图像600a-600d，其示出了地震勘测数据的绘制的公共深度点(cdp)堆叠。对于示例数据，在沙漠环境中进行了多次(这里是十三次)2d地震勘测。这些勘测是在注入二氧化碳之前获得的，并用于可重复性研究。输出vs记录是从堆叠在公共偏移量内的相关后合集(水平轴)获得的。cdp堆叠阶段包括对平面标定的较小静态校正、正常移出、逐迹线幅度平衡和静音。显示的时间间隔(垂直轴)为0ms至1900ms，每100ms显示一条定时线。图6a至图6d中的圆圈602表示每个空间位置的目标储层。

图6a示出了根据本公开的实施方式的非vs控制部分的数据图600a。图6a示出了在多个反射器上的良好信号连续性。

图6b示出了根据本公开的实施方式的在将传统vs与图6a的数据一起使用之后的数据图600b。使用传统vs可以看到增加的噪声。

图6c示出了根据本公开的实施方式的在使用利用图6a的数据的具有匹配滤波器的vs的数据图600c。注意，降低了总噪声，但是响应类似于图6a的响应。

图6d示出了根据本公开的实施方式的使用具有匹配滤波器加上具有图6a的数据的时间相关的平滑处理的vs的数据图600d。图6d示出了与数据图600b或600c的信噪比相比的改进的信噪比、最佳的连续性和最强的信号。

图7a至图7d是根据本公开的实施方式的图像700a-700d，其示出了针对从1120ms至1300ms的十三个样本地震勘测中的每一个绘制的公共深度点(cdp)迹线。显示的时间间隔(垂直轴)为0ms至1900ms，每100ms显示一条定时线。对于示例数据，在非vs控制区域中的勘测s1-s6和s7-s13(水平轴)之间，特别是在指示的目标储层(例如，图7a至图7d中的圆圈702)处，示出了重复性问题。圆圈702指示从前述的十三次地震勘测收集的目标储层的波形。包括具有匹配滤波器的vs的处理流程的应用似乎在目标储层水平具有改进的可重复性。

图7a示出了根据本公开的实施方式的非vs的数据图700a。注意在勘测s1-s6和s7-s13之间的所示目标储层(圆圈702)中示出的重复性问题。

图7b示出了根据本公开的实施方式的在将传统vs与图7a的数据一起使用之后的数据图700b。使用传统vs可以看到增加的噪声。

图7c示出了根据本公开的实施方式的在将具有匹配滤波器的vs与具有图7a的数据一起使用之后的数据图700c。注意，降低了总噪声，但是响应类似于图7a的响应。

图7d示出了在使用具有匹配滤波器加上时间相关的平滑处理的vs以及图7a的数据之后的数据图700d。图7d示出了图7a至图7d中最佳的可重复性。

可以使用在所示储层周围短窗口内计算的归一化均方根(nrms)对图4a至图4e、图6a-图6d以及图7a-图7d的观察结果进行量化。没有指定基线，因为在所有示例勘测之间都计算了nrms。考虑总共13次勘测，这表示每个cdp上有68个nrms组合。

图8是示出根据本公开的实施方式的与每种方法相关联的nrms值的直方图的曲线图800。非vs控制部分802(蓝线)示出了在分开勘测s1-s6和s7-s13的十七个月中的双峰nrms分布(40％变化)。通过传统vs804(红线)也观察到了类似的分布，尽管比非vs窄。幅度谱替换(asr)vs806(黑线)(具有匹配滤波器的vs)和具有时间权重vs808(浅蓝线)的asr(具有匹配滤波器和时间相关的平滑处理的vs)基本上将nrms以29％降低到单个峰值(也就是说，具有匹配滤波器和时间相关的平滑处理的vs在所有四个数据集中具有最小的nrms值，表明改进的可重复性)。图例中显示了分布的中间值。这些结果支持在图4a至图4e以及图6a至图6d中得出的观察结果，并示出了具有匹配滤波器的vs如何显著减少由多个月的时间间隔分开的地震勘测之间的不可重复性。

图9是示出根据本公开的实施方式的通过计算相邻cdp之间的nrms(nrmsc)在下行直接p波近似之间的图像质量进行比较的曲线图900。该新的连续性度量(nrmsc)很有用，因为其具有与nrms相同的输入时间窗口和相同的单位，并且旨在测量图像反射器的连续性。如图所示，具有匹配滤波器加上时间相关的平滑处理的vs具有最小的nrmsc值。与其他三个分布(非vs控制部分902(蓝线)、asrvs906(黑线)(具有匹配滤波器的vs)以及具有时间权重vs908(浅蓝线)的asr)相比，加窗的传统vs904(红线)的图像质量较差。这与在图6a和图6d中所作的观察一致。图像连续性改善了在曲线图900中向左移动的可能性。

图10是示出根据本公开的实施方式的用于提供与所描述的算法、方法、功能、处理、流程和过程相关联的计算功能的计算机实现的系统1000的示例的框图。在所示的实施方式中，系统1000包括计算机1002和网络1030。

示出的计算机1002旨在包括任意计算设备，例如服务器、台式计算机、膝上型/笔记本计算机、无线数据端口、智能电话、个人数字助理(pda)、平板计算机、这些设备内的一个或多个处理器、另一计算设备、或计算设备的组合(包括计算设备的物理或虚拟实例，或计算设备的物理或虚拟实例的组合)。另外，计算机1002可以包括可以接受用户信息的输入设备(例如键区、键盘、触摸屏、另一输入设备或输入设备的组合)、以及输出设备(该输出设备在图形型用户界面(ui)(或gui)或其他ui上传送与计算机1002的操作相关联的信息，包括数字数据、视觉、音频、其他类型的信息或这些类型的信息的组合)。

计算机1002可以用作用于执行本公开中描述的主题的客户端、网络组件、服务器、数据库或其他持久性存储器的分布式计算系统中的角色、其他角色、或多个角色的组合。示出的计算机1002可通信地与网络1030耦接。在一些实施方式中，计算机1002的一个或多个组件可以被配置为在包括基于云计算、本地、全局、或其他环境在内的环境(或者环境的组合)中操作。

从较高的层次来看，计算机1002是可操作为接收、发送、处理、存储或管理与所描述的主题相关联的数据和信息的电子计算设备。根据一些实施方式，计算机1002还可以包括或可通信地耦接到服务器，包括应用服务器、电子邮件服务器、web服务器、缓存服务器、流传输数据服务器或其他服务器、或服务器的组合。

计算机1002可以通过网络1030(例如，从在另一台计算机1002上执行的客户端软件应用)接收请求，并且通过使用软件应用或软件应用的组合处理接收到的请求来响应接收到的请求。另外，还可以从内部用户(例如，从命令控制台或通过其他内部访问方法)、外部或第三方、或其它实体、个人、系统或计算机向计算机1002发送请求。

计算机1002的每个组件可以使用系统总线1003进行通信。在一些实施方式中，计算机1002的任意或所有组件(包括硬件、软件或硬件和软件的组合)可以使用应用编程接口(api)1012、服务层1013、或api1012和服务层1013的组合，通过系统总线1003进行接口交互。api1012可以包括针对例程、数据结构和对象类的规范。api1012可以是独立于或依赖于计算机语言，并且指的是完整的接口、单个功能或甚至是一组api。服务层1013向计算机1002或可通信地耦接到计算机1002的其他组件(无论是否被示出)提供软件服务。计算机1002的功能可以对于使用该服务层1013的所有服务消费者是可访问的。软件服务(例如由服务层1013提供的软件服务)通过定义的接口提供可重用的、定义的功能。例如，接口可以是以java、c++、另一计算语言编写的软件、或者可以是以可扩展标记语言(xml)格式、另一格式或格式组合提供数据的计算语言的组合编写的软件。虽然被示为计算机1002的集成组件，但是备选实施方式可以将api1012和/或服务层1013示为作为相对于计算机1002的其他组件或可通信地耦接到计算机1002的其他组件(无论是否被示出)独立的组件。此外，在不脱离本公开的范围的情况下，api1012和/或服务层1013的任意或所有部分可以被实现为另一软件模块、企业应用或硬件模块的子模块或副模块。

计算机1002包括接口1004。虽然被示为单个接口1004，但是可以根据计算机1002的特定需要、期望或特定实现而使用两个或更多个接口1004。计算机1002使用接口1004与在分布式环境中通信地链接到网络1030的另一计算系统(无论是否示出)通信。通常，接口1004可操作为与网络1030通信并且包括以软件、硬件或软件和硬件的组合编码的逻辑。更具体地，接口1004可以包括支持与通信相关联的一个或多个通信协议的软件，使得网络1030或接口1004的硬件可操作为在所示出的计算机1002内部和外部传送物理信号。

计算机1002包括处理器1005。虽然被示为单个处理器1005，但是可以根据计算机1002的特定需要、期望或特定实现而使用两个或更多个处理器1005。通常，处理器1005执行指令并操纵数据以执行计算机1002的操作以及如本公开中所描述的任何算法、方法、功能、处理、流程和过程。

计算机1002还包括数据库1006(该数据库1006可以保存计算机1002的数据)、通信地链接到网络1030的另一个组件(无论是否示出)、或者计算机1002和另一个组件的组合。例如，数据库1006可以是存储与本公开一致的数据的内部存储器、常规或其他类型的数据库。在一些实施方式中，根据计算机1002的特定需要、期望或特定实现和所描述的功能，数据库1006可以是两个或更多个不同数据库类型(例如，混合的内部存储器和常规数据库)的组合。尽管被示出为单个数据库1006，根据计算机1002的特定需要、期望或特定实现和所描述的功能，可以使用相似或不同类型的两个或更多个数据库。虽然数据库1006被示出为计算机1002的集成组件，但是在备选实施方式中，数据库1006可以在计算机1002的外部。如所示出的，数据库1006保存原始地震数据1016(例如，炮集)和处理后的地震图像1018(例如，vs炮集)。

计算机1002还包括存储器1007(该存储器1007可以保存计算机1002的数据)、通信地链接到网络1030的另一个组件或多个组件(无论是否示出)、或者计算机1002和另一个组件的组合。存储器1007可以存储与本公开一致的任意数据。在一些实施方式中，根据计算机1002的特定需要、期望或特定实现和所描述的功能，存储器1007可以是两个或更多个不同类型的存储器的组合(例如，半导体和磁存储器的组合)。尽管被示出为单个存储器1007，根据计算机1002的特定需要、期望或特定实现和所描述的功能，可以使用相似或不同类型的两个或更多个存储器1007。虽然存储器1007被示为计算机1002的集成组件，但是在备选实施方式中，存储器1007可以在计算机1002的外部。

应用1008是根据计算机1002的特定需要、期望或特定实现提供功能(尤其是关于本公开中描述的功能)的算法软件引擎。例如，应用1008可以用作一个或多个组件、模块、应用等。此外，尽管被示为单个应用1008，但是应用1008可以被实现为计算机1002上的多个应用1008。另外，虽然被示出为与计算机1002集成，但是在备选实施方式中，应用1008可以在计算机1002的外部。

计算机1002还可以包括电源1014。电源1014可以包括可以被配置为用户或非用户可更换的可再充电或不可再充电电池。在一些实施方式中，电源1014可以包括电力转换或管理电路(包括再充电、备用或其他电力管理功能)。在一些实施方式中，电源1014可以包括电源插头，以允许计算机1002插入墙上插座或另一电源以例如为计算机1002供电或为可再充电电池充电。

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

对于所描述的方法，计算资源取决于要处理的数据量。例如，如果接收到的地震数据小于100gb，则当前的台式或便携式计算机就足以执行所述的处理。但是，如果接收到的地震数据大于1tb，则优选具有高性能计算系统/群集的超级计算来增强计算性能。

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

例如，在第一实施方式中，一种计算机实现的方法，包括：将接收到的炮集分类为公共接收器集；根据合成数据计算地震波场直接到达的目标三维(3d)幅度谱；针对每个接收器计算场数据中地震波场直接到达的3d幅度谱；计算从场数据的3d幅度谱到目标响应的匹配滤波器；将计算出的匹配滤波器应用于与公共接收器集分离的下行地震波场，以产生滤波后的下行地震波场；对滤波后的下行地震波场进行时间相关的平滑处理，以产生平滑的下行地震波场；以及计算从公共接收器集分离的上行地震波场与平滑下行地震波场之间的互相关。

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

可与以下任何特征组合的第一特征，还包括：选择接近偏移的第一地震波场直接到达作为初始恒定的p波速度模型；将恒定的p波速度模型从地表源更新到埋藏式接收器；以及执行3d傅立叶变换，将更新后的恒定的p波速度模型转换到频率-波数域。

第二特征，可以与先前或之后的特征中的任何一个组合，还包括：确定公差值和场数据的3d幅度谱与合成数据的目标3d幅度谱的计算差之间的比较小于或等于公差值；以及重新计算从场数据的3d幅度谱到目标响应的匹配滤波器。

可与前述或后续特征中的任何一个组合的第三特征，还包括：将公共接收器集分离为下行地震波场和上行地震波场。

可以与之前或之后的任何特征组合的第四特征，还包括：对每个下行地震波场进行时间门控；以及对每个时间门控的下行地震波场执行空间相关的孔径滤波。

第五特征，可以与先前或之后的特征中的任何一个组合，还包括对上行的地震波场进行地滚波去除。

第六特征，可以与先前或之后的特征中的任何一个组合，进一步包括：堆叠虚拟源炮集；以及开始输出虚拟源炮集。

在第二实施方式中，一种存储一个或多个指令的非暂时性计算机可读介质，所述指令可由计算机系统执行以：将接收到的炮集分类为公共接收器集；根据合成数据计算地震波场直接到达的目标三维(3d)幅度谱；针对每个接收器计算场数据中地震波场直接到达的3d幅度谱；计算从场数据的3d幅度谱到目标响应的匹配滤波器；将计算出的匹配滤波器应用于与公共接收器集分离的下行地震波场，以产生滤波后的下行地震波场；对滤波后的下行地震波场进行时间相关的平滑处理，以生成平滑的下行地震波场；以及计算从公共接收器集分离的上行地震波场与平滑下行地震波场之间的互相关。

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

可与以下任何特征组合的第一特征，还包括一个或多个可由计算机系统执行的指令以：选择接近偏移的第一地震波场直接到达作为初始恒定的p波速度模型；将恒定的p波速度模型从地表源更新到埋藏式接收器；和执行3d傅里叶变换，将更新后的恒定的p波速度模型转换到频率波数域。

可与先前或随后的特征中的任何一个组合的第二特征，还包括可由计算机系统执行以执行以下操作的一个或多个指令：选择接近偏移的第一地震波场直接到达作为初始恒定的p波速度模型；将恒定的p波速度模型从地表源更新到埋藏式接收器；以及执行3d傅里叶变换，将更新后的恒定的p波速度模型转换到频率波数域。

第三特征，可以与先前或以下特征中的任何一个组合，还包括一个或多个可由计算机系统执行的指令，以：确定公差值和场数据的3d幅度谱与合成数据的目标3d幅度谱的计算差之间的比较小于或等于公差值；以及重新计算从场数据的3d幅度谱到目标响应的匹配滤波器。

可以与先前或随后的特征中的任何一个组合的第四特征，还包括一个或多个可由计算机系统执行以将公共接收器集分离为下行地震波场和上行地震波场的指令。

可以与先前或以下特征中的任何一个组合的第五特征，还包括一个或多个可由计算机系统执行以执行以下指令的指令：对每个下行地震波场进行时间门控；以及对每个时间门控的下行地震波场执行空间相关的孔径滤波。

第六特征，可以与先前或随后的特征中的任何一个组合，进一步包括一个或多个可由计算机系统执行的指令，以：堆叠虚拟炮集；以及开始输出虚拟炮集。

在第三实施方式中，计算机实现的方法包括：计算机存储器；以及硬件处理器，其与计算机存储器可互操作地耦合，并被配置为：将接收到的炮集分类为公共接收器集；根据合成数据计算地震波场直接到达的目标三维(3d)幅度谱；针对每个接收器计算场数据中地震波场直接到达的3d幅度谱；计算从场数据的3d幅度谱到目标响应的匹配滤波器；将计算出的匹配滤波器应用于与公共接收器集分离的下行地震波场，以产生滤波后的下行地震波场；对滤波后的下行地震波场进行时间相关的平滑处理，以生成平滑的下行地震波场；以及计算从公共接收器集分离的上行地震波场与平滑下行地震波场之间的互相关。

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

第一项功能可与以下任何功能组合，进一步配置为：选择接近偏移的第一地震波场直接到达作为初始恒定的p波速度模型；将恒定的p波速度模型从地表源更新到埋藏式接收器；以及执行3d傅立叶变换，将更新后的恒定的p波速度模型转换到频率波数域。

第二个功能可以与之前或之后的任何功能组合，进一步配置为：确定公差值和场数据的3d幅度谱与合成数据的目标3d幅度谱的计算差之间的比较小于或等于公差值；以及重新计算从场数据的3d幅度谱到目标响应的匹配滤波器。

可与之前或之后的特征中的任何一个组合的第三特征还被配置为将公共接收器集分离为下行地震波场和上行地震波场。

可以与先前或随后的任何特征组合的第四特征，进一步配置为：对每个下行地震波场进行时间门控；以及对每个时间门控的下行地震波场执行空间相关的孔径滤波。

可以与先前或随后的特征中的任何一个组合的第五特征还被配置为在上行的地震波场上执行地滚波去除。

第六项功能可以与之前或之后的任何功能组合，进一步配置为：堆叠虚拟炮集；和开始输出虚拟炮集。

在一些实施方式中，所描述的方法可被配置为向计算机实现的控制器、数据库或其他计算机实现的系统发送消息、指令或其他通信，以动态地启动对另一计算机实现的系统的控制，控制或使得另一计算机实现的系统执行计算机实现的功能/操作或其他功能/操作。例如，可以传输基于数据的操作、操作、输出或与gui的交互，以使与计算机、数据库、网络或其他基于计算机的系统关联的操作执行存储效率、数据检索或与本公开一致的其他操作。在另一示例中，与任何所示的gui进行交互可以自动地导致从gui传输的一个或多个指令触发对数据的请求、数据的存储、数据的分析或与本公开一致的其他操作。

在一些情况下，传输的指令可以导致关于有形的真实世界的计算设备或其他设备的控制、操作、修改、增强或其他操作。例如，所描述的gui可以发送请求以减慢或加速计算机数据库磁/光盘驱动器，关闭/激活计算系统，使网络接口设备禁用、节制或增加通过网络连接允许的数据带宽，或发出可听见/可视警报(例如，机械警报/发光设备)作为关于与和所描述的方法相关的(或者与和所描述的方法相关的计算系统交互的)计算系统的结果、行为、确定或分析的通知。

在一些实施方式中，所描述的方法的输出可以用于动态地影响、指导、控制、影响或管理与碳氢化合物生产、分析和回收有关的有形的、真实世界的设备，或者用于与本公开一致的其他目的。例如，可以将从正在进行的钻井作业中接收的实时数据合并到使用所述方法进行的分析中。所产生的2d/3d地震图像(包括实时数据)的提高的质量可用于多种目的。例如，根据所描述的方法的生成结果，可以修改井眼轨迹，可以增加或减少钻孔速度，可以停止钻孔，可以激活/停用警报(例如视觉、听觉或语音警报)，可以影响(例如，停止、重新启动、加速或减少)精炼或泵送作业。其他示例可以包括在检测到地下障碍物时(例如，通过视觉、听觉或声音警报)提醒地质导向和定向钻井人员。在一些实施方式中，可以将所描述的方法学集成为动态计算机实现的控制系统的一部分，以控制、影响或与和本公开内容一致的任何与碳氢化合物有关的或其他有形的、真实世界的设备一起使用。

在本说明书中描述的主题和功能操作的实施可以实现在下述形式中：数字电子电路、有形体现的计算机软件或固件、计算机硬件，包括在本说明书中公开的结构及其结构等同物、或它们中的一个或多个的组合。描述的主题的软件实施可以被实现为在有形非暂时计算机可读介质上编码的一个或多个计算机程序，即计算机程序指令的一个或多个模块，所述程序用于由计算机或计算机实现的系统执行或者控制计算机或计算机实现的系统的操作。备选地或另外地，程序指令可以在人工生成的传播信号(例如，机器生成的电、光或电磁信号)上编码，所述信号被生成以对信息进行编码以传输到接收器装置，以供计算机或计算机实现的系统执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或计算机存储介质的组合。配置一个或多个计算机意味着一个或多个计算机安装了硬件、固件或软件(或硬件、固件和软件的组合)，以便当软件由一个或多个计算机执行时，执行特定的计算操作。

术语“实时”、“实时(快速)(rft)”、“接近实时(nrt)”、“准实时”或类似术语(如本领域的普通技术人员所理解的)意味着动作和响应在时间上接近，使得个人感知动作和响应基本上同时发生。例如，在个人做出了访问数据的动作之后对数据显示的响应的时间差(或用于启动显示)可以小于1毫秒(ms)、小于1秒(s)或小于5秒。尽管所请求的数据不需要被即时显示(或启动以显示)，但是考虑到所描述的计算系统的处理限制和例如收集、精确测量、分析、处理、存储或传输所需的时间，在没有任何有意的延迟的情况下显示(或启动以显示)该数据。

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

可以以任何形式的编程语言来编写计算机程序(也可以称作或描述为程序、软件、软件应用、单元、模块、软件模块、脚本代码或另一组件)，所述编程语言包括：编译或解释语言、或者声明或程序语言，并且可以以任何形式来部署计算机程序，包括部署为例如单独的程序或者用于计算环境的模块、组件或子例程。计算机程序可以(但不是必须)与文件系统中的文件相对应。程序可以存储在保持其它程序或数据(例如，存储在标记语言文档中的一个或多个脚本)的文件的一部分中、存储在专用于所讨论的程序的单个文件中、或者存储在多个协同文件中(例如，存储一个或多个模块、子程序或代码的一部分的文件)。计算机程序可以被部署为在一个计算机上或者在位于一个站点或分布在多个站点并且通过通信网络互连的多个计算机上执行。

尽管各图中所示的程序的部分被示为使用各种对象、方法或其他处理实现所描述的特征和功能的单独组件(例如单元或模块)，但是视情况，程序可以替代地包括多个子单元、子模块、第三方服务、组件、库和其他组件。相反，各种组件的特征和功能可以视情况组合成单个组件。可以统计地、动态地或者统计地且动态地确定用于进行计算确定的阈值。

所描述的方法、过程或逻辑流程表示与本公开一致的功能的一个或多个示例，并且不旨在将本公开限制为所描述或示出的实施方式，而是被赋予与所描述的原理和特征一致的最宽范围。所描述的方法、处理或逻辑流可以由一个或多个可编程计算机来执行，所述一个或多个可编程计算机执行一个或多个计算机程序以通过操作输入数据并且生成输出数据来执行功能。方法、处理或逻辑流也可以由专用逻辑电路(例如cpu、fpga或asic)来执行，并且计算机也可以实现为专用逻辑电路(例如cpu、fpga或asic)。

用于执行计算机程序的计算机可以基于通用或专用微处理器、这两者或其它类型的cpu。通常，cpu将从存储器接收指令和数据并写入存储器。计算机的必不可少的元件是用于执行指令的cpu和用于存储指令和数据的一个或更多个存储器设备。通常，计算机还将包括用于存储数据的一个和或更多个大容量存储设备(例如，磁盘、磁光盘或光盘)，或可操作耦接以便从所述一个或更多个大容量存储设备接收或向其发送数据。然而，计算机不需要具有这些设备。此外，可以将计算机嵌入到另一个设备中，例如移动电话，个人数字助理(pda)，移动音频或视频播放器，游戏控制台，全球定位系统(gps)接收器或便携式存储器存储设备。

用于存储计算机程序指令和数据的非暂时性计算机可读介质可以包括所有形式的永久性/非永久性或易失性/非易失性存储器、介质和存储器设备，例如包括半导体存储器设备，例如随机(ram)、只读存储器(rom)、相变存储器(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、可擦除可编程只读存储器(eprom)、电可擦除可编程只读存储器(eeprom)和闪存设备；磁设备，例如磁带、盒式磁带、磁带盒、内部/可移动盘；磁光盘；以及光学存储器设备，例如数字通用/视频光盘(dvd)、紧凑盘(cd)-rom、dvd+/-r、dvd-ram、dvd-rom、高分辨率/密度(hd)-dvd和蓝光/蓝光盘(bd)以及其他光学存储器技术。存储器可以存储各种对象或数据，包括：高速缓存器、类(class)、框架、应用、模块、备份数据、工作、网页、网页模板、数据结构、数据库表格、存储动态信息的知识库、或者包括任意参数、变量、算法、指令、规则、约束、引用在内的任意其它适当的信息。此外，存储器还可以包含其他适当的数据，例如日志、策略、安全或访问数据或报告文件。处理器和存储器可以由专用逻辑电路来补充或者并入到专用逻辑电路中。

为了提供与用户的交互，本说明书中描述的主题可以实现在计算机上，该计算机具有用于向用户显示信息的显示设备(例如，阴极射线管(crt)、液晶显示器(lcd)、发光二极管(led)或等离子监视器)和用户可以向计算机提供输入的键盘和指点设备(例如，鼠标、轨迹球或轨迹板)。还可以使用触摸屏(诸如具有压敏性的平板计算机表面，使用电容或电感测的多点触摸屏或其他类型的触摸屏)向计算机提供输入。可以使用其他类型的设备与用户交互。例如，向用户提供的反馈可以是任意形式的感官反馈(例如，视觉、听觉、触觉或反馈类型的组合)。可以以任意形式(包括声音，语音或触觉输入)来接收来自用户的输入。此外，计算机可以通过向用户使用的客户端计算设备发送文档或者从该设备接收文档来与用户交互(例如，通过响应于从用户的移动计算设备上的web浏览器接收到的请求而向所述web浏览器发送网页，来与用户交互)。

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

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

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

尽管本说明书包含许多具体实施细节，然而这些细节不应被解释为对要求保护的范围或任何发明构思的范围构成限制，而是用于说明特定于具体发明构思的具体实施例的特征。在单个实施方式中，还可以组合实现本说明书中在独立实施方式的上下文中描述的特定特征。相反的，单个实施方式的上下文描述的不同特征也可在多个实施方式中各自实现，或以任意子组合来实现。此外，虽然前述特征可以被描述为在某些组合中起作用并且甚至最初如此要求保护，但是来自所要求保护的组合的一个或多个特征在一些情况下可以从组合中删除，并且所要求保护的组合可以针对子组合或子组合的变体。

已经描述了本主题的特定实施方式。对于本领域技术人员显而易见的是，所描述的实施方式的其它实现，改变和置换在所附权利要求的范围内。尽管在附图和权利要求中以特定顺序描述了操作，这不应被理解为：为了实现期望的结果，要求按所示的特定次序或按顺序次序来执行这些操作，或者要求执行所有图示的操作(一些操作可以看作是可选的)。在某些情况下，多任务或并行处理(或者多任务和并行处理的组合)可以是优选地并且视情况来执行。

此外，在前述的实现中的各种系统模块和组件的分离或集成不应被理解为在所有实施方式中要求这样的分离或集成，并且应该理解的是，所描述的程序组件和系统一般可以一起集成在单个软件产品中或封装为多个软件产品。

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

此外，任何要求保护的实施方式被认为适用于至少一种计算机实施的方法；存储用于执行计算机实施的方法的计算机可读指令的非暂时性计算机可读介质；以及计算机系统，该系统包括与硬件处理器可互操作地耦接的计算机存储器，所述硬件处理器被配置为执行计算机实现的方法或存储在非暂时性计算机可读介质上的指令。