CN110068860B - 一种时移地震数据处理方法、装置及系统 - Google Patents

Abstract

本说明书实施例公开了一种时移地震数据处理方法、装置及系统，所述方法包括获取目标工区的待匹配地震数据以及参考地震数据；根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。利用本说明书各个实施例，可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

Description

一种时移地震数据处理方法、装置及系统

技术领域

本发明涉及地球物理勘探技术领域，特别地，涉及一种时移地震数据处理方法、装置及系统。

背景技术

时移地震勘探方法是油气田生产过程中，在油气藏开发的不同时期重复进行三维地震勘探的方法，不同时间的地震响应随时间的变化可以表征储藏内部油气的变化规律。同时，通过不同时间地震数据的差异分析，也可以描述油气藏内部物性参数的变化。

但是在实际应用中，时移地震数据往往是在不同时期进行间隔性采集并处理的，各时期采集的地震数据所对应的噪声、物理环境变数、近地表影响、采集仪器和处理参数等具有一定的差异性，使得油藏内部物性参数分析结果受到较大影响。针对三维叠后地震数据时差、振幅及相位进行分别归一化处理，可以在一定程度上消除不同时期数据处理及采集所造成的影响。但是该方法是对地震数据的整体处理，无法对局部参数进行单独调整，归一化处理精度有限。现有的匹配追踪方法虽然能够弥补上述归一化方法的不足，但对油气藏流体流动的前缘及储层内部物性变化的描述精度及准确性不够。

因此，本技术领域亟需一种更为准确的时移地震数据处理方法，以提高油藏内部物性参数分析的准确性。

发明内容

本说明书实施例的目的在于提供一种时移地震数据处理方法、装置及系统，可以提高油藏内部物性参数分析的准确性。

本说明书提供一种时移地震数据处理方法、装置及系统是包括如下方式实现的：

一种时移地震数据处理方法，包括：

获取目标工区的待匹配地震数据以及参考地震数据；

根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

本说明书提供的所述方法的另一个实施例中，所述对所述待匹配地震数据及参考地震数据进行归一化匹配处理，包括：

基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

根据所述Huber范数误差函数的最小化计算结果确定Huber范数的匹配滤波算子；

利用所述匹配滤波算子对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

本说明书提供的所述方法的另一个实施例中，所述对所述待匹配地震数据及参考地震数据进行归一化匹配处理，包括：

根据下述公式确定基于Huber范数的匹配滤波目标函数Obj的最小值，获得Huber范数匹配滤波算子A：

式中，m表示待匹配地震数据，d表示参考地震数据，ε为门槛值，N为离散数据点的个数；

根据上述公式确定的Huber范数匹配滤波算子A对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

本说明书提供的所述方法的另一个实施例中，所述门槛值采用下述方式确定：

其中，d_region为采样点的预设邻域范围内的待匹配地震数据，m_region为采样点的预设邻域范围内的参考地震数据，||·||_p为d_region-m_region的p范数，α∈(0,1)为振幅差异系数。

本说明书提供的所述方法的另一个实施例中，所述方法还包括：

判断所述匹配滤波后的地震数据相对所述参考地震数据的差异数据体是否满足预设储层开发认识精度条件；

如果判断结果为否，则调整所述振幅差异系数重新计算所述匹配滤波算子，直至所述差异数据体满足预设储层开发认识精度条件。

本说明书提供的所述方法的另一个实施例中，所述对所述待匹配地震数据及参考地震数据进行归一化匹配处理之前，还包括：

对所述待匹配地震数据及参考地震数据进行规则化处理；

相应的，所述对所述待匹配地震数据、参考地震数据进行归一化匹配处理包括根据Huber范数匹配滤波方法对规则化处理后的待匹配地震数据及参考地震数据进行归一化处理。

另一方面，本说明书实施例还提供一种时移地震数据处理装置，所述装置包括：

数据获取模块，用于获取目标工区的待匹配地震数据以及参考地震数据；

匹配滤波处理模块，用于根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

本说明书提供的所述装置的另一个实施例中，所述匹配滤波处理模块包括：

误差函数构建单元，用于基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

滤波算子确定单元，用于根据所述Huber范数误差函数的最小化计算结果，确定Huber范数的匹配滤波算子；

匹配滤波处理单元，用于利用所述匹配滤波算子对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

另一方面，本说明书实施例还提供一种时移地震数据处理设备，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

获取目标工区的待匹配地震数据以及参考地震数据；

根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

另一方面，本说明书实施例还提供一种时移地震数据处理系统，所述系统包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个实施例所述方法的步骤。

本说明书一个或多个实施例提供的一种时移地震数据处理方法、装置及系统，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

附图说明

为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本说明书提供的一种时移地震数据处理方法实施例的流程示意图；

图2为本说明书提供的另一种时移地震数据处理方法实施例的流程示意图；

图3为本说明书提供的一个实施例中的匹配滤波处理的流程示意图；

图4为本说明书提供的一种时移地震数据处理装置实施例的模块结构示意图；

图5为根据本说明书的一个示例性实施例的服务器的示意结构图。

具体实施方式

为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。

获取时移地震数据，对时移地震成像结果相减后，油气藏的静态性质，如构造、岩性性质等被消去，可以获得如流体饱和度、压力、温度等油气藏动态流体性质的直接成像结果。因此，在油气藏生产中以时间延迟的形式进行重复地震勘探，可以实现对油气藏生产引起的油气藏内部物性参数(如流体饱和度、压力和温度等)变化的描述，并可以追踪到流体流动的前缘，实现对油气藏进行动态监测与管理。

但实际应用中，时移地震数据通常是间隔性采集、处理的，两次采集很难保证各项因素完全一致。如地下水位的变化会造成地表条件的不一致，环境的变化会造成环境噪声的不一致，震源类型、激发位置或放炮方式的不同会造成能量分布的不一致，采集仪器型号的不同会造成不同的仪器噪声与不同的频谱特征，观测系统的差别会导致两个数据体难以比较等等。所有这些因素不一致都会造成反演结果之间的差异可能无实际物理意义。

相应的，本说明书实施例提供了一种时移地震数据处理方法，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

图1是本说明书提供的所述一种时移地震数据处理方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。

具体的一个实施例如图1所示，本说明书提供的时移地震数据处理方法的一个实施例中，所述方法可以包括：

S102：获取目标工区的待匹配地震数据以及参考地震数据。

时移地震数据可以包括在不同时间对所述目标工区进行重复性地震观测所获取的地震勘探数据。所述待匹配地震数据及所述参考地震数据可以为待分析时间段内不同时期采集的地震数据，以用于对相应时期的油气藏内部物性参数变化进行分析。

通常时移地震数据采集和处理年份不同，由于技术的发展，新的地震不可能与原有的地震采用同样的采集、处理参数。相应的，本说明书的一个实施例中，在获取待匹配地震数据及参考地震数据后，可以先对时移地震数据进行规则化处理。所述规则化处理可以包括对不同时间获得的地震数据的地震道位置、道数、采样间隔、采样点数等进行统一化处理。

在进行地震数据间的匹配滤波之前，可以先对地震数据进行规则化处理。如可以保证不同数据中地震道位置及道数相同，每个地震道重采样到相同的采样间隔及采样点数。同时，还可以设定各地震数据的处理流程一致，保证数据不受地震数据处理及数据范围的影响。通过对地震数据进行规则化处理，可以有效降低因地震数据采集和前期处理的差异性对时移地震数据分析结果准确性的影响。

S104：根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

所述Huber范数匹配滤波方法可以包括对小于等于门槛值的数据利用L2范数进行处理、对大于门槛值的数据利用L1范数进行处理。利用基于Huber范数的匹配滤波方法对待匹配地震数据及参考地震数据进行归一化处理，可以实现对差异性较小的数据进行平滑处理同时，还可以有效保留差异性较大的数据。从而在实现对储层内部物性变化刻画的稳定性的同时，进一步有效提高对油气藏流体流动的前缘及储层内部物性变化描述的精度。

图2是本说明书提供的另一种时移地震数据处理方法实施例流程示意图。如图2所示，本说明书的一个实施例中，可以利用下述方式对进行归一化处理：

S1042：基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

S1044：根据所述Huber范数误差函数的最小化计算结果，确定Huber范数的匹配滤波算子；

S1046：利用所述匹配滤波算子对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

可以根据Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数，然后，可以通过求解该Huber范数误差函数的极小值来确定Huber范数的匹配滤波算子。然后，可以根据计算获得的匹配滤波算子对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

然后，可以将匹配滤波后的地震数据相对所述参考地震数据的差异数据体，作为所述目标工区的时移地震差异数据。然后，对目标工区的时移地震差异数据进行分析，可以确定出待匹配地震数据与参考地震数据所对应的两个不同时期的油气藏内部物性参数变化情况。

一些实施方式中，在构建Huber范数误差函数时，可以利用参考地震数据与待匹配地震数据的差值属性数据体，也可以采用均方根等属性数据体，这里不做限定。相应的，在获得匹配滤波算子之后，可以根据误差函数构建时对应的数据体属性，计算匹配滤波后的地震数据相对所述参考地震数据的差异数据体，作为所述目标工区的时移地震差异数据。

本说明书的一个或者多个实施例中，对时移地震数据的匹配滤波处理可以表示为线性系统d＝Am，其中，m待处理时移地震户数，A为匹配滤波算子，d为参考地震数据。一些实施方式中，基于Huber范数的匹配滤波目标函数可以表示为：

其中，

式中，

为Huber函数，ε为门槛值，N为离散数据点的个数。通过求解上述目标函数Obj的最小化问题，可以得到匹配滤波算子A。

然后，可以根据得到的匹配滤波算子A对待匹配地震数据进行处理，获得匹配滤波后的地震数据Am。相应的，可以将匹配滤波后的地震数据相对所述参考地震数据的差值数据体，作为所述目标工区的时移地震差异数据。

所述门槛值可以根据需要预先设定以及调整。优选的，本说明书的另一个实施例中，可以根据采样点周围一定邻域范围的地震数据来确定该采样点对应的门槛值，以使得不同的采样点对应的门槛值同该采样点周围的地震数据相适应。所述采样点周围一定邻域范围的地震数据可以为以采样点为中心的一定邻域范围的地震数据振幅属性值的均值、均方根值等，邻域范围可以根据实际需要确定，这里不做限定。

地震数据的振幅往往存在较强的空变形，例如浅层地震数据真实异常差异变化门槛值可能和深层非真实数据差异表现为同等级别的振幅差异，使得利用单一门槛值约束的匹配滤波计算方法的计算准确性受限。如单一门槛值设置过高，待匹配地震数据与参考地震数据的振幅差异性较小的地层的储层内部物性参数无法有效刻画；而门槛值设置的过低，则待匹配地震数据与参考地震数据的振幅差异性较大的地层的匹配滤波后获得差异性数据体噪声太大，从而严重影响了最终结果的准确性。

相应的，本说明书一个或者多个实施例中通过进一步考虑不同地层位置处待匹配地震数据与参考地震数据的实际差异性，在滤波处理过程中，通过结合各实际地震数据来自适应的调整门槛值ε，可以使得各采样点的滤波处理与其所处的地层所对应的地震数据特征相适应，从而进一步提高最终处理结果的准确性。

一些实施方式中，所述门槛值还可以采用下述方式确定：

其中，d_region为采样点的预设邻域范围内的待匹配地震数据，m_region为采样点的预设邻域范围内的参考地震数据，||·||_p为d_region-m_region的p范数，α∈(0,1)为振幅差异系数。

本说明书的另一个实施例中，所述方法还可以包括：

判断所述匹配滤波后的地震数据相对所述参考地震数据的差异数据体是否满足预设储层开发认识精度条件；

如果判断结果为否，则调整所述振幅差异系数重新计算所述匹配滤波算子，直至所述差异数据体满足预设储层开发认识精度条件。

图3表示本说明书一个实施例中提供的匹配滤波处理流程示意图。如图3所示，在获得差异数据体后，可以先分析差异数据体是否符合已有储层开发的认识，如是否符合现有生产井的生产信息、对现有油藏流体前缘认识等。如果不符合相应的认识精度条件，则可以根据认识分析结果调整匹配滤波器参数，重新计算匹配滤波算子。反复进行这一过程，最终使得差异数据体满足已有储层开发认识，完成数据归一化计算。

一些实施方式中，可以调整所述门槛值的大小，重新计算获得匹配滤波算子。本说明书的一个实施例中，对于上述根据一定阈值范围的地震数据自适应的确定的门槛值，可以通过调整振幅差异系数α，来重新计算获得匹配滤波算子。

上述实施例提供的方案，通过利用Huber匹配滤波方法对时移地震数据进行匹配滤波处理，可以实现对差异性较小的数据进行平滑处理同时，有效保留差异性较大的数据。从而在过滤噪声的同时，进一步精确有效的刻画出储层内部物性的变化特征。进一步的，在匹配滤波处理过程中，可以根据实际地层的地震数据特征自适应的调整门槛值的大小，使得门槛值的设置更符合地层的地震数据变化特征，从而进一步提高各储层内物性参数刻画的精确性。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

本说明书一个或多个实施例提供的一种时移地震数据处理方法，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，可以同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以油藏内部物性参数分析结果的准确性。

基于上述所述的时移地震数据处理方法，本说明书一个或多个实施例还提供一种时移地震数据处理装置。所述的装置可以包括使用了本说明书实施例所述方法的系统、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图4表示说明书提供的一种时移地震数据处理装置实施例的模块结构示意图，如图4所示，所述装置可以包括：

数据获取模块202，可以用于获取目标工区的待匹配地震数据以及参考地震数据；

匹配滤波处理模块204，可以用于根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

本说明书的另一个实施例中，所述匹配滤波处理模块204可以包括：

误差函数构建单元，可以用于基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

滤波算子确定单元，可以用于根据所述Huber范数误差函数的最小化计算结果，确定Huber范数的匹配滤波算子；

匹配滤波处理单元，可以用于利用所述匹配滤波算子对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

本说明书一个或多个实施例提供的一种时移地震数据处理装置，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种时移地震数据处理设备，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

获取目标工区的待匹配地震数据以及参考地震数据；

根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据。

需要说明的，上述所述的设备根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。

本说明书实施例所提供的方法实施例可以在移动终端、计算机终端、服务器或者类似的运算装置中执行。以运行在服务器上为例，图5是应用本说明书实施例的布尔电路编码服务器的硬件结构框图。如图5所示，服务器10可以包括一个或多个(图中仅示出一个)处理器20(处理器20可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)、用于存储数据的存储器30、以及用于通信功能的传输模块40。本邻域普通技术人员可以理解，图5所示的结构仅为示意，其并不对上述电子装置的结构造成限定。例如，服务器10还可包括比图5中所示更多或者更少的组件，例如还可以包括其他的处理硬件，如数据库或多级缓存、GPU，或者具有与图5所示不同的配置。

存储器30可用于存储应用软件的软件程序以及模块，如本发明实施例中的搜索方法对应的程序指令/模块，处理器20通过运行存储在存储器30内的软件程序以及模块，从而执行各种功能应用以及数据处理。存储器30可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中，存储器30可进一步包括相对于处理器20远程设置的存储器，这些远程存储器可以通过网络连接至计算机终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

传输模块40用于经由一个网络接收或者发送数据。上述的网络具体实例可包括计算机终端的通信供应商提供的无线网络。在一个实例中，传输模块40包括一个网络适配器(Network Interface Controller，NIC)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，传输模块40可以为射频(Radio Frequency，RF)模块，其用于通过无线方式与互联网进行通讯。

所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。

上述实施例所述的一种时移地震数据处理设备，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

本说明书还提供一种时移地震数据处理系统，所述系统可以为单独的时移地震数据处理系统，也可以应用在多种地震勘探数据处理系统中。所述的系统可以为单独的服务器，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的服务器集群、系统(包括分布式系统)、软件(应用)、实际操作装置、逻辑门电路装置、量子计算机等并结合必要的实施硬件的终端装置。所述时移地震数据处理系统可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个或者多个实施例中所述方法的步骤。

需要说明的，上述所述的系统根据方法或者装置实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

上述实施例所述的一种时移地震数据处理系统，可以利用基于Huber范数的匹配滤波方法对时移地震数据中不同数据体进行归一化匹配处理。基于Huber范数的匹配滤波方法，可以结合L1范数和L2范数匹配算法的特点，同时保证对差异数据边界的刻画精度以及对储层内部物性变化刻画的稳定性，从而可以提高利用时移地震数据进行油藏内部物性参数分析的准确性。

需要说明的是，本说明书上述所述的装置或者系统根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类、存储介质+程序实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本说明书实施例并不局限于必须是符合标准数据模型/模板或本说明书实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书的可选实施方案范围之内。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。

本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

以上所述仅为本说明书的实施例而已，并不用于限制本说明书。对于本领域技术人员来说，本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书的权利要求范围之内。

Claims (6)

1.一种时移地震数据处理方法，其特征在于，包括：

获取目标工区的待匹配地震数据以及参考地震数据；

根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据，包括：

根据下述公式确定基于Huber范数的匹配滤波目标函数Obj的最小值，获得Huber范数匹配滤波算子A：

式中，m表示待匹配地震数据，d表示参考地震数据，ε为门槛值，N为离散数据点的个数；其中，所述门槛值采用下述方式确定：

其中，d_region为采样点的预设邻域范围内的待匹配地震数据，m_region为采样点的预设邻域范围内的参考地震数据，||·||_p为d_region-m_region的p范数，α∈(0,1)为振幅差异系数；

利用所述匹配滤波算子A对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：

判断所述匹配滤波后的地震数据相对所述参考地震数据的差异数据体是否满足预设储层开发认识精度条件；

如果判断结果为否，则调整所述振幅差异系数重新计算所述匹配滤波算子，直至所述差异数据体满足预设储层开发认识精度条件。

3.根据权利要求1所述的方法，其特征在于，所述对所述待匹配地震数据及参考地震数据进行归一化匹配处理之前，还包括：

对所述待匹配地震数据及参考地震数据进行规则化处理；

相应的，所述对所述待匹配地震数据、参考地震数据进行归一化匹配处理包括根据Huber范数匹配滤波方法对规则化处理后的待匹配地震数据及参考地震数据进行归一化处理。

4.一种时移地震数据处理装置，其特征在于，所述装置包括：

数据获取模块，用于获取目标工区的待匹配地震数据以及参考地震数据；

匹配滤波处理模块，用于根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据；

其中，所述匹配滤波处理模块包括：

误差函数构建单元，用于基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

滤波算子确定单元，用于根据下述公式确定基于Huber范数的匹配滤波目标函数Obj的最小值，获得Huber范数匹配滤波算子A：

式中，m表示待匹配地震数据，d表示参考地震数据，ε为门槛值，N为离散数据点的个数；其中，所述门槛值采用下述方式确定：

其中，d_region为采样点的预设邻域范围内的待匹配地震数据，m_region为采样点的预设邻域范围内的参考地震数据，||·||_p为d_region-m_region的p范数，α∈(0,1)为振幅差异系数；

匹配滤波处理单元，用于利用所述匹配滤波算子A对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

5.一种时移地震数据处理设备，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

获取目标工区的待匹配地震数据以及参考地震数据；

根据Huber范数匹配滤波方法对所述待匹配地震数据及参考地震数据进行归一化匹配处理，获得所述目标工区的时移地震差异数据，包括：

基于Huber范数匹配滤波方法构建所述待匹配地震数据相对所述参考地震数据的Huber范数误差函数；

根据下述公式确定基于Huber范数的匹配滤波目标函数Obj的最小值，获得Huber范数匹配滤波算子A：

式中，m表示待匹配地震数据，d表示参考地震数据，ε为门槛值，N为离散数据点的个数；其中，所述门槛值采用下述方式确定：

其中，d_region为采样点的预设邻域范围内的待匹配地震数据，m_region为采样点的预设邻域范围内的参考地震数据，||·||_p为d_region-m_region的p范数，α∈(0,1)为振幅差异系数；

利用所述匹配滤波算子A对所述待匹配地震数据进行处理，获得匹配滤波后的地震数据。

6.一种时移地震数据处理系统，其特征在于，所述系统包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1-3中任意一项所述方法的步骤。

Images