一种地震波形匹配方法及装置与流程

本申请涉及地震数据处理与分析技术领域，特别涉及一种地震波形匹配方法及装置。

背景技术：

基于地震数据确定地层的构造形态是地震勘探的主要目的之一，这一过程通常被称为地震构造解释。在地震构造解释过程中，例如，层位解释或断层解释，都需要对地震数据中不同地震道对应的地震波形进行匹配。地震波形匹配是指对两个地震波形之间的特征点的对应关系进行分析，以提取地震波形的特征点匹配关系序列。

现有的地震波形匹配方法的主要过程是：计算地震数据中目标地震波形上每一个目标采样点分别与参考地震波形上指定参考采样点之间的距离，即以采样点为中心的相同时窗窗长内的波形之间的欧式距离；将距离最小的目标采样点作为指定参考采样点对应的匹配点；重复上述步骤，从目标地震波形中搜寻参考地震波形上每一个参考采样点对应的匹配点。针对信噪比较低的地震数据，采用现有的方法进行地震波形匹配，匹配过程可能容易受到噪声干扰，出现匹配错误的情况。

技术实现要素：

本申请实施例的目的是提供一种地震波形匹配方法及装置，以提高地震波形匹配的准确度。

为解决上述技术问题，本申请实施例提供一种地震波形匹配方法及装置是这样实现的：

一种地震波形匹配方法，提供有地震数据，其中，所述地震数据包括：参考地震道数据和目标地震道数据；所述方法包括：

从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同；

从所述参考地震波形序列中确定参考地震稀疏脉冲序列，以及从所述目标地震波形序列中确定目标地震稀疏脉冲序列；

确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，建立动态时间弯曲模板；所述动态时间弯曲模板用于表示确定波形匹配路径过程中所采用的搜索模式；

基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列。

优选方案中，所述从所述参考地震波形序列中确定参考地震稀疏脉冲序列，包括：

基于预设拟合尺度集中指定拟合尺度，采用最小二乘拟合方法确定指定拟合多项式，并分别对所述指定拟合多项式进行一阶求导和二阶求导，得到与所述指定拟合多项式对应的一阶多项式导数和二阶多项式导数；

将所述一阶多项式导数与所述参考地震波形序列的卷积结果作为所述参考地震波形序列的一阶导数，以及将所述二阶多项式导数与所述参考地震波序列的卷积结果作为所述参考地震波形序列的二阶导数；

基于所述参考地震波形序列的一阶导数和二阶导数，分别确定所述参考地震波形序列中极值点的位置和极值点类型；根据所述参考地震波形序列中极值点的位置和极值点类型，确定与所述指定拟合尺度对应的参考地震稀疏脉冲序列；其中，所述参考地震稀疏脉冲序列中参考脉冲点为所述参考地震波形序列中极值点；

按照预设拟合尺度集中拟合尺度由小到大的顺序，遍历确定与当前拟合尺度对应的当前参考地震稀疏脉冲序列，直至当前参考地震稀疏脉冲序列中的参考脉冲点个数与前一个参考地震稀疏脉冲序列中的参考脉冲点个数之间的差值小于预设差异阈值为止，并将当前参考地震稀疏脉冲序列作为最终的参考地震稀疏脉冲序列。

优选方案中，所述从所述目标地震波形序列中确定目标地震稀疏脉冲序列，包括：

基于预设拟合尺度集中指定拟合尺度，采用最小二乘拟合方法确定指定拟合多项式，并分别对所述指定拟合多项式进行一阶求导和二阶求导，得到与所述指定拟合多项式对应的一阶多项式导数和二阶多项式导数；

将所述一阶多项式导数与所述目标地震波形序列的卷积结果作为所述目标地震波形序列的一阶导数，以及将所述二阶多项式导数与所述目标地震波序列的卷积结果作为所述目标地震波形序列的二阶导数；

基于所述目标地震波形序列的一阶导数和二阶导数，分别确定所述目标地震波形序列中极值点的位置和极值点类型；根据所述目标地震波形序列中极值点的位置和极值点类型，确定与所述指定拟合尺度对应的目标地震稀疏脉冲序列；其中，所述目标地震稀疏脉冲序列中目标脉冲点为所述目标地震波形序列中极值点；

按照预设拟合尺度集中拟合尺度由小到大的顺序，遍历确定与当前拟合尺度对应的当前目标地震稀疏脉冲序列，直至当前目标地震稀疏脉冲序列中的目标脉冲点个数与前一个目标地震稀疏脉冲序列中的目标脉冲点个数之间的差值小于预设差异阈值为止，并将当前目标地震稀疏脉冲序列作为最终的目标地震稀疏脉冲序列。

优选方案中，所述确定波形距离类型和波形匹配步长集，包括：

获取所述地震数据的信噪比，基于预设判别条件和所述信噪比，确定所述波形距离类型和最大波形匹配步长；所述最大波形匹配步长的取值为大于1的整数；

将取值为1至最大波形匹配步长的整数的波形匹配步长所构成的集合作为所述波形匹配步长集。

优选方案中，所述预设判别条件为：

当所述信噪比大于预设信噪比阈值时，将所述波形距离类型设置为平均绝对误差方法，确定所述最大波形匹配步长的取值范围为2～4；或者，

当所述信噪比小于或等于预设信噪比阈值时，将所述波形距离类型设置为高阶绝对误差方法，确定所述最大波形匹配步长的取值范围为5至预设采样点数量。

优选方案中，所述基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列，包括：

基于所述波形距离类型，分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，并存储于波形距离矩阵；所述波形距离矩阵中行号与所述参考地震稀疏脉冲序列的参考脉冲点序列号相对应，所述波形距离矩阵中列号与所述目标地震稀疏脉冲序列的目标脉冲点序列号相对应，所述波形距离矩阵中的任一元素的值为所述波形距离；

在所述波形距离矩阵中确定正向起始元素；

基于所述动态时间弯曲模板，从所述波形距离矩阵中确定所述正向起始元素指向的目标元素；

遍历从所述波形距离矩阵中确定当前目标元素所指向的下一个目标元素；

将各个所述目标元素对应的行号和列号作为所述目标地震波形匹配路径序列中各个匹配路径中参考脉冲点序列号和目标脉冲点序列号。

优选方案中，所述基于所述波形距离类型，分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，包括：

针对所述参考地震稀疏脉冲序列中指定参考脉冲点和所述目标地震稀疏脉冲序列中指定目标脉冲点，基于所述波形距离类型，计算以所述指定参考脉冲点为中心的预设时窗内的参考脉冲点集合，与以所述指定目标脉冲点为中心的预设时窗内的目标脉冲点集合之间的波形距离。

优选方案中，所述在所述波形距离矩阵中确定正向起始元素，包括：

将所述波形距离矩阵中第一行的各个元素中波形距离最小的元素作为所述正向起始元素。

优选方案中，所述基于所述动态时间弯曲模板，从所述波形距离矩阵中确定所述正向起始元素所指向的目标元素，包括：

基于所述动态弯曲模板，将所述正向起始元素所指向的所述波形距离矩阵中的元素中波形距离最小的元素作为所述目标元素。

优选方案中，基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列，还包括：

基于所述波形距离类型，分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，并将所述波形距离存储于波形距离矩阵；所述波形距离矩阵中行号与所述参考地震稀疏脉冲序列的参考脉冲点序列号相对应，所述波形距离矩阵中列号与所述目标地震稀疏脉冲序列的目标脉冲点序列号相对应，所述波形距离矩阵中的任一元素的值为所述波形距离；

对所述波形距离矩阵中第一行和第一列的各个元素中指定元素的值作为预先设置的累积距离矩阵中与所述指定元素相同行号和相同列号对应的元素的值；其中，所述累积距离矩阵与所述波形距离矩阵的行数和列数相同；

基于所述动态时间弯曲模板和所述波形距离矩阵，确定所述累积距离矩阵中除第一行和第一列的元素以外的其他元素的值；

基于所述动态时间弯曲模板和所述累积距离矩阵，确定所述目标地震波形匹配路径序列。

优选方案中，所述基于所述动态时间弯曲模板和所述累积距离矩阵，确定所述目标地震波形匹配路径序列，包括：

在所述累积距离矩阵中确定反向起始元素；

基于所述动态时间弯曲模板，从所述累积距离矩阵中确定所述反向起始元素所指向的目标元素；

遍历从所述累积距离矩阵中确定当前目标元素所指向的下一个目标元素；

将各个所述目标元素对应的行号和列号作为所述目标地震波形匹配路径序列中各个匹配路径中参考脉冲点序列号和目标脉冲点序列号。

一种地震波形匹配装置，所述装置提供地震数据，其中，所述地震数据包括：参考地震道数据和目标地震道数据；所述装置包括：地震波形序列获取模块、地震脉冲序列确定模块、动态时间弯曲模板建立模块和匹配处理模块；其中，

所述地震波形序列获取模块，用于从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同；

所述地震脉冲序列确定模块，用于分别从所述参考地震波形序列和所述目标地震波形序列确定参考地震稀疏脉冲序列和目标地震稀疏脉冲序列；

所述动态时间弯曲模板建立模块，用于确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，建立动态时间弯曲模板；所述动态时间弯曲模板用于表示确定波形匹配路径过程中所采用的搜索模式；

所述匹配处理模块，用于基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列。

本申请实施例了一种地震波形匹配方法及装置，可以从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同；可以分别从所述参考地震波形序列和所述目标地震波形序列确定参考地震稀疏脉冲序列和目标地震稀疏脉冲序列；确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，可以建立动态时间弯曲模板；所述动态时间弯曲模板用于表示确定波形匹配路径过程中所采用的搜索模式；基于所述波形距离类型和所述动态时间弯曲模板，可以对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，可以得到目标地震波形匹配路径序列。通过所确定的波形匹配步长集建立的动态时间弯曲模板，进行匹配处理，可以跳过脉冲序列中噪声脉冲，从而可以提高地震波形匹配的准确度。

附图说明

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

图1是本申请一种地震波形匹配方法实施例的流程图；

图2是本申请实施例中参考地震波形序列和目标地震波形序列的示意图；

图3是本申请实施例中参考地震稀疏脉冲序列和目标地震稀疏脉冲序列的示意图；

图4是本申请实施例中正向动态时间弯曲模板的示意图；

图5是本申请实施例中反向动态时间弯曲模板的示意图；

图6是本申请实施例中波形距离矩阵的示意图；

图7是本申请实施例中累积距离矩阵的示意图；

图8是本申请实施例中地震波形匹配路径的示意图；

图9是本申请实施例中地震波形匹配路径序列的示意图；

图10是本申请地震波形匹配装置实施例的组成结构图。

具体实施方式

本申请实施例提供一种地震波形匹配方法及装置。

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

本申请实施例提供一种地震波形匹配方法。所述方法可以提供有地震数据，其中，所述地震数据可以包括：参考地震道数据和目标地震道数据。

在一种实施方式中，所述地震数据的采样时间间隔可以为1毫秒。所述参考地震道数据可以包括：参考采样点，以及与所述参考采样点对应的采样时间和振幅。所述目标地震道数据可以包括：目标采样点，以及与所述目标采样点对应的采样时间和振幅。

图1是本申请一种地震波形匹配方法实施例的流程图。所述地震波形匹配方法，包括以下步骤。

步骤s101：从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同。

具体地，可以从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列。所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同。其中，所述参考地震波形序列和所述目标地震波形序列均可以为一维数组。数组中的元素与采样点一一对应。

例如，图2是本申请实施例中参考地震波形序列和目标地震波形序列的示意图。图2中波形a表示参考地震波形序列，波形b表示目标地震波形序列。纵向坐标表示采样时间，单位为毫秒。

步骤s102：从所述参考地震波形序列中确定参考地震稀疏脉冲序列，以及从所述目标地震波形序列中确定目标地震稀疏脉冲序列。

在一种实施方式中，从所述参考地震波形序列中确定参考地震稀疏脉冲序列，具体可以包括，基于预设拟合尺度集中指定拟合尺度，采用最小二乘拟合方法确定指定拟合多项式，并分别对所述指定拟合多项式进行一阶求导和二阶求导，得到与所述指定拟合多项式对应的一阶多项式导数和二阶多项式导数。可以将所述一阶多项式导数与所述参考地震波形序列的卷积结果作为所述参考地震波形序列的一阶导数，以及将所述二阶多项式导数与所述参考地震波序列的卷积结果作为所述参考地震波形序列的二阶导数。基于所述参考地震波形序列的一阶导数和二阶导数，可以分别确定所述参考地震波形序列中各个极值点的位置和极值点类型。根据所述参考地震波形序列中各个极值点的位置和极值点类型，确定与所述指定拟合尺度对应的参考地震稀疏脉冲序列。其中，所述参考地震稀疏脉冲序列中参考脉冲点为所述参考地震波形序列中极值点。可以按照预设拟合尺度集中拟合尺度由小到大的顺序，遍历确定与当前拟合尺度对应的当前参考地震稀疏脉冲序列，直至当前参考地震稀疏脉冲序列中的参考脉冲点个数与前一个参考地震稀疏脉冲序列中的参考脉冲点个数之间的差值小于预设差异阈值为止，并将当前参考地震稀疏脉冲序列作为最终的参考地震稀疏脉冲序列。

在一种实施方式中，从所述目标地震波形序列中确定目标地震稀疏脉冲序列，具体可以包括，基于预设拟合尺度集中指定拟合尺度，采用最小二乘拟合方法确定指定拟合多项式，并分别对所述指定拟合多项式进行一阶求导和二阶求导，得到与所述指定拟合多项式对应的一阶多项式导数和二阶多项式导数。可以将所述一阶多项式导数与所述目标地震波形序列的卷积结果作为所述目标地震波形序列的一阶导数，以及将所述二阶多项式导数与所述目标地震波序列的卷积结果作为所述目标地震波形序列的二阶导数。基于所述目标地震波形序列的一阶导数和二阶导数，分别确定所述目标地震波形序列中各个极值点的位置和极值点类型。根据所述目标地震波形序列中各个极值点的位置和极值点类型，确定与所述指定拟合尺度对应的目标地震稀疏脉冲序列。其中，所述目标地震稀疏脉冲序列中目标脉冲点为所述目标地震波形序列中极值点。可以按照预设拟合尺度集中拟合尺度由小到大的顺序，遍历确定与当前拟合尺度对应的当前目标地震稀疏脉冲序列，直至当前目标地震稀疏脉冲序列中的目标脉冲点个数与前一个目标地震稀疏脉冲序列中的目标脉冲点个数之间的差值小于预设差异阈值为止，并将当前目标地震稀疏脉冲序列作为最终的目标地震稀疏脉冲序列。

例如，图3是本申请实施例中从图2的参考地震波形序列和目标地震波形序列中分别确定的参考地震稀疏脉冲序列和目标地震稀疏脉冲序列的示意图。图3中脉冲a表示参考地震稀疏脉冲序列，脉冲b表示目标地震稀疏脉冲序列。图3中的纵坐标表示采样时间，单位为毫秒。按照采样时间由小到大，可以对所述参考地震稀疏脉冲序列中参考脉冲点排序，以及对目标地震稀疏脉冲序列中目标脉冲点排序并编号。如此，所述参考地震稀疏脉冲序列可以由脉冲a中参考脉冲点s0、s1、s2、...、s30构成。所述目标地震稀疏脉冲序列可以由脉冲b中目标脉冲点t0、t1、t2、...、t29构成。

步骤s103：确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，建立动态时间弯曲模板；所述动态时间弯曲模板用于表示确定波形匹配路径过程中所采用的搜索模式。

在一种实施方式中，可以获取所述地震数据的信噪比，基于预设判别条件和所述信噪比，可以确定所述波形距离类型和最大波形匹配步长。所述最大波形匹配步长的取值可以为大于1的整数。所述波形距离类型可以表示两个波形之间的波形距离的计算方法。可以将取值为1至最大波形匹配步长的整数的波形匹配步长所构成的集合作为所述波形匹配步长集。

在一种实施方式中，所述预设判别条件可以包括：当所述信噪比大于预设信噪比阈值时，可以将所述波形距离类型设置为平均绝对误差方法，可以确定所述最大波形匹配步长的取值范围为2～4；或者，当所述信噪比小于或等于预设信噪比阈值时，可以将所述波形距离类型设置为高阶绝对误差方法，可以确定所述最大波形匹配步长的取值范围为5至预设采样点数量。所述预设信噪比阈值的取值范围可以为50～60分贝。

在一种实施方式中，基于所述波形匹配步长集，建立动态时间弯曲模板，具体可以包括，基于所述波形匹配步长集，可以确定多种搜索模式和所述搜索模式的数量。可以将各种所述搜索模式所构成的模板作为所述动态时间弯曲模板。其中，所述动态时间弯曲模板可以包括：正向动态时间弯曲模板和反向动态时间弯曲模板。例如，图4是本申请实施例中正向动态时间弯曲模板的示意图。图5是本申请实施例中反向动态时间弯曲模板的示意图。所述波形匹配步长集可以包括值为1、2、3、4、5的波形匹配步长。如图4所示，所述正向动态时间弯曲模板中包括的搜索模式的数量为9，各种搜索模式可以分别为：

g(i-1,j-1)→g(i,j)；

g(i-2,j-1)→g(i,j)；

g(i-3,j-1)→g(i,j)；

g(i-4,j-1)→g(i,j)；

g(i-5,j-1)→g(i,j)；

g(i-1,j-2)→g(i,j)；

g(i-1,j-3)→g(i,j)；

g(i-1,j-4)→g(i,j)；

g(i-1,j-5)→g(i,j)。

如图5所示，所述反向动态时间弯曲模板中包括的搜索模式的数量为9，各种搜索模式可以分别为：

g(i,j)→g(i-1,j-1)；

g(i,j)→g(i-2,j-1)；

g(i,j)→g(i-3,j-1)；

g(i,j)→g(i-4,j-1)；

g(i,j)→g(i-5,j-1)；

g(i,j)→g(i-1,j-2)；

g(i,j)→g(i-1,j-3)；

g(i,j)→g(i-1,j-4)；

g(i,j)→g(i-1,j-5)。

步骤s104：基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列。

具体地，基于所述波形距离类型，可以分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，并存储于波形距离矩阵；所述波形距离矩阵中行号与所述参考地震稀疏脉冲序列的参考脉冲点序列号相对应，所述波形距离矩阵中列号与所述目标地震稀疏脉冲序列的目标脉冲点序列号相对应，所述波形距离矩阵中的任一元素的值为所述波形距离。可以在所述波形距离矩阵中确定正向起始元素。基于所述动态时间弯曲模板，可以从所述波形距离矩阵中确定所述正向起始元素所指向的目标元素。遍历从所述波形距离矩阵中确定当前目标元素所指向的下一个目标元素。可以将各个所述目标元素对应的行号和列号作为所述目标地震波形匹配路径序列中各个匹配路径中参考脉冲点序列号和目标脉冲点序列号。

在一种实施方式中，基于所述波形距离类型，分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，具体可以包括，针对所述参考地震稀疏脉冲序列中指定参考脉冲点和所述目标地震稀疏脉冲序列中指定目标脉冲点，基于所述波形距离类型，可以计算以所述指定参考脉冲点为中心的预设时窗内的参考脉冲点集合，与以所述指定目标脉冲点为中心的预设时窗内的目标脉冲点集合之间的波形距离。

在一种实施方式中，在所述波形距离矩阵中确定正向起始元素，具体可以包括，可以将所述波形距离矩阵中第一行的各个元素中波形距离最小的元素作为所述正向起始元素。

在另一种实施方式中，基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列，具体还可以包括：基于所述波形距离类型，可以分别计算所述参考地震稀疏脉冲序列中各个参考脉冲点与所述目标地震稀疏脉冲序列中各个目标脉冲点之间的波形距离，并将所述波形距离存储于波形距离矩阵。所述波形距离矩阵中行号与所述参考地震稀疏脉冲序列的参考脉冲点序列号相对应，所述波形距离矩阵中列号与所述目标地震稀疏脉冲序列的目标脉冲点序列号相对应，所述波形距离矩阵中的任一元素的值为所述波形距离。可以对所述波形距离矩阵中第一行和第一列的各个元素中指定元素的值作为预先设置的累积距离矩阵中与所述指定元素相同行号和相同列号对应的元素的值。其中，所述累积距离矩阵与所述波形距离矩阵的行数和列数相同。基于所述动态时间弯曲模板和所述波形距离矩阵，可以确定所述累积距离矩阵中除第一行和第一列的元素以外的其他元素的值。基于所述动态时间弯曲模板和所述累积距离矩阵，可以确定所述目标地震波形匹配路径序列。

在一种实施方式中，基于所述动态时间弯曲模板和所述累积距离矩阵，确定所述目标地震波形匹配路径序列，具体可以包括，可以在所述累积距离矩阵中确定反向起始元素。基于所述动态时间弯曲模板，可以从所述累积距离矩阵中确定所述反向起始元素所指向的目标元素。可以遍历从所述累积距离矩阵中确定当前目标元素所指向的下一个目标元素。可以将各个所述目标元素对应的行号和列号作为所述目标地震波形匹配路径序列中各个匹配路径中参考脉冲点序列号和目标脉冲点序列号。

在一种实施方式中，基于所述动态时间弯曲模板和所述波形距离矩阵，确定所述累积距离矩阵中除第一行和第一列的元素以外的其他元素的值，具体可以包括，可以采用下述公式计算所述累积距离矩阵中除第一行和第一列的元素以外的其他元素中一元素的值：

d(i,j)＝d(i,j)+min{d(m,n)}

其中，d(i,j)表示所述累积距离矩阵中第i行和第j列对应的元素的值，d(i,j)表示所述波形距离矩阵中第i行和第j列对应的元素的值，d(m,n)表示所述累积距离矩阵中第m行和第n列对应的元素的值，其中，i＞1，j＞1，当m＝i-1,i-2,...,i-k时，n＝j-1；当m＝i-1时，n＝j-1,j-2,...,j-k。

在一种实施方式中，在所述累积距离矩阵中确定反向起始元素，具体可以包括，可以将所述累积距离矩阵中最后一行的各个元素中波形距离最小的元素作为所述反向起始元素。

例如，图6是本申请实施例中基于图3得到的波形距离矩阵的示意图。图7是本申请实施例中基于图3和图6得到的累积距离矩阵的示意图。图6和图7中的行号和列号分别表示所述参考地震稀疏脉冲序列的参考脉冲点序列号，以及所述目标地震稀疏脉冲序列的目标脉冲点序列号。图6和图7中的灰度值分别表示波形距离和累积距离。图8是本申请实施例中基于图7得到的地震波形匹配路径的示意图。图9是本申请实施例中基于图8得到的地震波形匹配关系和地震波形匹配路径序列的示意图。图9中(a)和(b)分别为基于图8得到的地震波形匹配关系和地震波形匹配路径序列的示意图。图9中(a)中的横坐标和纵坐标分别为振幅和采样时间。其中，采样时间的单位为毫秒(ms)。图9中(a)的各个参考脉冲点与各个目标脉冲点之间连线表示所述地震波形匹配关系，即连线两端的参考脉冲点与目标脉冲点是互相匹配的。在图9中(b)的地震波形匹配路径序列中，参考列中元素的值分别与参考采样点排序后的编号一一对应。根据所述地震波形匹配关系，可以将与参考脉冲点相匹配的目标脉冲点排序后的编号作为目标列中元素的值。其中，互相匹配的参考脉冲点的编号和目标脉冲点的编号分别作为在参考列和目标列的同一行处元素的值。如参考脉冲点s0与目标脉冲点t0匹配，参考列和目标列的第一行处元素的分别为0和0。某些参考脉冲点没有相互匹配的目标脉冲点，这时目标列中与这些参考脉冲点的编号对应的行处的元素的值为-1。如参考脉冲点s6、s9和s10没有相互匹配的目标脉冲点，目标列中与这些参考脉冲点的编号对应的行处的元素的值为-1。从图9中，可以看出，目标地震稀疏脉冲序列中的目标脉冲点t1、t2、t8没有对应匹配的参考采样点，表明在采用本申请实施例的地震波形匹配方法进行匹配处理时，可以跳过如目标脉冲点t1、t2、t8这样的噪声，从而可以提高波形匹配的准确度。

所述地震波行匹配方法实施例，可以从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同；可以分别从所述参考地震波形序列和所述目标地震波形序列确定参考地震稀疏脉冲序列和目标地震稀疏脉冲序列；确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，可以建立动态时间弯曲模板；所述动态时间弯曲模板用于表示确定波形匹配路径过程中所采用的搜索模式；基于所述波形距离类型和所述动态时间弯曲模板，可以对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，可以得到目标地震波形匹配路径序列。通过所确定的波形匹配步长集建立的动态时间弯曲模板，进行匹配处理，可以跳过脉冲序列中噪声脉冲，从而可以提高地震波形匹配的准确度。

图10是本申请地震波形匹配装置实施例的组成结构图。所述地震波形匹配装置可以提供地震数据，其中，所述地震数据包括：参考地震道数据和目标地震道数据。如图10所示，所述地震波形匹配装置可以包括：地震波形序列获取模块100、地震脉冲序列确定模块200、动态时间弯曲模板建立模块300和匹配处理模块400。

所述地震波形序列获取模块100，可以用于从所述参考地震道数据中获取参考地震波形序列，以及从所述目标地震道数据中获取目标地震波形序列；所述参考地震波形序列与所述目标地震波形序列分别对应的采样时间段相同。

所述地震脉冲序列确定模块200，可以用于分别从所述参考地震波形序列和所述目标地震波形序列确定参考地震稀疏脉冲序列和目标地震稀疏脉冲序列。

所述动态时间弯曲模板建立模块300，可以用于确定波形距离类型和波形匹配步长集；基于所述波形匹配步长集，建立动态时间弯曲模板；所述动态时间弯曲模板可以表示确定波形匹配路径过程中所采用的搜索模式。

所述匹配处理模块400，可以用于基于所述波形距离类型和所述动态时间弯曲模板，对所述参考地震稀疏脉冲序列和目标地震稀疏脉冲序列进行匹配处理，得到目标地震波形匹配路径序列。

所述地震波形匹配实施例与所述地震波形匹配方法实施例相对应，可以实现所述地震波形匹配方法实施例，并取得方法实施例的技术效果。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(programmablelogicdevice,pld)(例如现场可编程门阵列(fieldprogrammablegatearray，fpga))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片pld上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logiccompiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(hardwaredescriptionlanguage，hdl)，而hdl也并非仅有一种，而是有许多种，如abel(advancedbooleanexpressionlanguage)、ahdl(alterahardwaredescriptionlanguage)、confluence、cupl(cornelluniversityprogramminglanguage)、hdcal、jhdl(javahardwaredescriptionlanguage)、lava、lola、myhdl、palasm、rhdl(rubyhardwaredescriptionlanguage)等，目前最普遍使用的是vhdl(very-high-speedintegratedcircuithardwaredescriptionlanguage)与verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

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

上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitorymedia)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络pc、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

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

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。