本发明涉及地球物理勘探技术领域，尤其涉及一种数据处理方法和装置。
背景技术：
在地震数据处理过程中，利用地震数据中面波进行近地表结构反演已经得到广泛应用。目前，在利用面波进行近地表结构反演中通常需要建立地层模型，地层模型建立的好坏经常直接影响到反演运算效率和反演结果。因此，需要建立高质量的地层模型。
现有技术中，地层模型的建立的通常可以利用面波不同频率波长，以及不同频率波长与深度的之间的转换系数(即半波长系数)进行半波长法近地表结构反演，得到相应控制点的近地表结构；然后，利用所述目标面波频散点近地表结构进行建模，继而可以基于确定的地层模型确定近地表结构反演结果数据。但目前在利用面波进行近地表结构反演中，常常基于经验采用统一半波长系数，与实际情况存在差异，可靠性差，影响后续建立的地层模型的精度，进而无法得到准确的得到近地表结构反演数据。
技术实现要素：
本申请的目的是提供一种数据处理方法和装置，可以获取准确的获取半波长系数，进而可以提高地层模型，和近地表结构反演数据的精度。
本申请提供的数据处理方法和装置是这样实现的：
一种数据处理方法，所述方法包括：
基于预设半波长系数范围设置多个目标面波频散点第一预设数量的半波长系数；
利用所述目标面波频散点的层速度、频率和半波长系数计算所述目标面波频散点的所述半波长系数所对应的表层厚度；
基于预设控制点表层速度和表层厚度按照预设深度内插间隔进行表层速度扩展处理，得到第二预设数量的控制点的表层数据，所述控制点的表层数据包括一一对应的表层速度和表层厚度；
计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数；
将最高的相关系数所对应的半波长系数作为所述目标面波频散点在相应表层厚度的半波长系数。
一种数据处理装置，所述装置包括：
系数设置模块，用于基于预设半波长系数范围设置多个目标面波频散点第一预设数量的半波长系数；
第一计算模块，用于利用所述目标面波频散点的层速度、频率和半波长系数计算所述目标面波频散点的所述半波长系数所对应的表层厚度；
表层速度扩展处理模块，用于基于预设控制点表层速度和表层厚度按照预设深度内插间隔进行表层速度扩展处理，得到第二预设数量的控制点的表层数据，所述控制点的表层数据包括一一对应的表层速度和表层厚度；
第二计算模块，用于计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数；
半波长系数确定模块，用于将最高的相关系数所对应的半波长系数作为所述目标面波频散点在相应表层厚度的半波长系数。
本申请提供的技术方案可以通过计算设置的第一预设数量的半波长系数所对应的目标面波频散点在同一表层厚度所对应的层速度与实际地层所对应的控制点的表层速度之间的相关性，来确定出最符合地层实际情况的半波长系数。与现有技术相比，利用本申请提供的技术方案可以准确的确定出符合地层实际情况的半波长系数，进而可以提高后续地层模型，和近地表结构反演数据的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的数据处理方法的一种实施例的流程图；
图2是本申请提供的基于预设半波长系数范围设置目标面波频散点第一预设数量的半波长系数一种实施例的流程示意图；
图3是现有技术中实际测量的微测井以及现有固定半波长系数的表层厚度与表层速度曲线的对比图；
图4是现有技术中实际测量的微测井以及本申请确定的半波长系数所对应的表层厚度与表层速度曲线的对比图；
图5是本申请提供的数据处理装置的一种实施例中的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
以下首先介绍本申请一种数据处理方法的一种实施例。图1是本申请提供的数据处理方法的一种实施例的流程图，本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的如图1所示，所述方法可以包括：
S110：基于预设半波长系数范围设置多个目标面波频散点第一预设数量的半波长系数。
本申请实施例中，具体的，所述多个目标面波频散点可以包括面波所对应的频散点。如图2所示，图2是本申请提供的基于预设半波长系数范围设置目标面波频散点第一预设数量的半波长系数一种实施例的流程示意图，具体的，可以包括：
S210：利用所述预设半波长系数范围和所述目标面波频散点的数量确定半波长系数最小值，半波长增量和半波长系数最大值；
在实际应用中，面波波长对深度转换系数介于(0，1)，通常在[0.5，0.7]，也称之为面波的半波长系数。在一个具体的实施例中，假设所述预设半波长系数范围为[0.5，0.7]，相应的，所述半波长系数最小值为0.5，所述半波长系数最大值为0.7，假设目标面波频散点为3个(在实际应用中一般目标面波频散点的数量较多，这里仅仅以目标面波频散点为3个作为示例。)，相应的，所述半波长增增量可以为0.1。
S220：基于所述半波长系数最小值，所述半波长增量和所述半波长系数最大值确定多个目标面波频散点的第一半波长系数组。
结合上述3个目标面波频散点示例，可以得到3个目标面波频散点的第一半波长系数组为(0.5，0.6，0.7)。
S230：从所述第一半波长系数组中选取一个半波长系数，作为替换系数。
结合上述3个目标面波频散点示例，假设选取0.5为替换系数。
S240：以所述替换系数替换所述第一半波长系数组中与所述替换系数不同的半波长系数，得到所述多个目标面波频散点的第一数量的第二半波长系数组；
结合上述3个目标面波频散点示例，以所述替换系数替换所述多个目标面波频散点中任一目标面波频散点所对应的半波长系数可以得到2个第二半波长系数组(0.5，0.5，0.7)和(0.5，0.6，0.5)。
S250：以所述替换系数替换所述第二半波长系数组中与所述替换系数不同的半波长系数，得到所述多个目标面波频散点第二数量的第三半波长系数组。
结合上述3个目标面波频散点示例，第二数量的第三半波长系数组可以为一个半波长系数组(0.5，0.5，0.5)。
S260：重复替换半波长系数的操作至得到的半波长系数组中的半波长系数都为所述替换系数。
270：将当前获得的半波长系数组作为待替换半波长系数组；
S280：选取所述第一半波长系数组中新的半波长系数作为替换系数，所述新的半波长系数为第一半波长系数组中未作为替换系数的半波长系数；
S290：重复上述替换处理的步骤至预设收敛条件。
其中，所述预设收敛条件包括选取所述第一半波长系数组中的全部半波长系数作为替换系数，且最后一个替换系数替换得到半波长系数组中的全部半波长系数为所述最后替换系数。
S300：将所述第一半波长系数组和达到所述预设收敛条件之前每次替换得到的半波长系数组中的半波长系数作为目标面波频散点第一预设数量的半波长系数。
结合上述3个目标面波频散点示例，可以选取0.6为替换系数，以上述的4个半波长系数组(0.5，0.6，0.7)、(0.5，0.5，0.7)和(0.5，0.6，0.5)以及(0.5，0.5，0.5)为待替换半波长系数组，进而可以得到7个第二半波长系数组可以为(0.6，0.6，0.7)、(0.5，0.6，0.6)、(0.6，0.5，0.7)、(0.5，0.5，0.6)、(0.6，0.6，0.5)、(0.5，0.6，0.6)和(0.6，0.5，0.5)；进一步的，以该7个第二半波长系数组为待替换半波长系数组可以得到一个第三半波长系数组(0.6，0.6，0.6)。
进一步的，可以选取0.7为替换系数，以上述12个半波长系数组作为待替换半波长系数组，可以得到12个第二半波长系数组为(0.7，0.6，0.7)、(0.5，0.7，0.7)、(0.7，0.5，0.7)、(0.7，0.6，0.5)、(0.5，0.7，0.5)、(0.7，0.5，0.5)、(0.5，0.7，0.5)、(0.6，0.7，0.7)、(0.7，0.6，0.6)、(0.5，0.7，0.6)、(0.7，0.6，0.5)和(0.6，0.7，0.5)。进一步的，以该12个第二半波长系数组为待替换半波长系数组可以得到3个第三半波长系数组为(0.7，0.7，0.5)、(0.7，0.7，0.7)(0.7，0.7，0.6)。
有上述可见，上述3个目标面波频散点可以对应27个半波长系数组。上述设置半波长系数的过程中可以全面的覆盖多种不同的半波长系数，进而可以保证得到的目标面波频散点的半波长系数中涵盖了符合地层实际情况的半波长系数。
S120：利用所述目标面波频散点的层速度、频率和半波长系数计算所述目标面波频散点的所述半波长系数所对应的表层厚度。
本申请实施例中，利用所述目标面波频散点的层速度、频率和半波长系数计算所述目标面波频散点的所述半波长系数所对应的表层厚度可以包括采用下述公式计算：
Hij＝Cij*Vi/fi
其中，Hij表示第i个目标面波频散点的第j个半波长系数所对应的表层厚度；Cij表示第i个目标面波频散点的第j个半波长系数；Vi表示第i个目标面波频散点的层速度；fi表示第i个目标面波频散点的频率。
这里由于同一目标面波频散点对应j个(即多个)半波长系数，相应的，同一目标面波频散点会对应不同的表层厚度，相应的同一厚度会对应多个目标面波频散点。
另外，在实际应用中，在面波频散分析中获得的面波相速度为可以作为平均相速度，将平均相速度转换为平均纵波速度后，可以根据平均纵波速度和层速度之间的换算关系，计算获得得到层速度。
S130：基于预设控制点表层速度和表层厚度按照预设深度内插间隔进行表层速度扩展处理，得到第二预设数量的控制点的表层数据，所述控制点的表层数据包括一一对应的表层速度和表层厚度。
在实际的勘探采集过冲中，常常会选取一些控制点来采集数据，本申请实施例中可以选取几个控制点作为预设控制点，所述预设控制点表层速度和表层厚度可以直接在采集数据中获取，然后，基于预设控制点表层速度和表层厚度按照预设深度内插间隔进行表层速度扩展处理，得到第二预设数量的控制点的表层数据，所述控制点的表层数据包括一一对应的表层速度和表层厚度。
具体的，这里的预设深度内插间隔可以结合实际应用中计算得到的目标面波频散点的半波长系数所对应的表层厚度进行设置，使得到第二预设数量的控制点的表层数据中的表层厚度与目标面波频散点的半波长系数所对应的表层厚度向对应。
另外，这个预设控制点一般选取多个，而多个预设控制点的表层厚度和表层速度按照预设深度内插间隔进行表层速度扩展处理，得到的控制点的表层数据之间会存在差异，相应的，同一控制点在同一表层厚度会对应多个表层速度。
S140：计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数。
本申请实施例中，可以计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数。具体的，可以包括采用下述公式计算：
其中，rhi表示表层厚度为h时，多组模板表层数据中表层速度与第i个目标面波频散点目标面波频散点的层速度之间的相关系数；xhik表示表层厚度为h时，第i个目标面波频散点的第k个半波长系数所对应的层速度；k表示表层厚度为h时，第i个目标面波频散点所对应的第k个表层速度；n表示表层厚度为h时，第i个目标面波频散点所对应的层速度数量；yhk表层厚度为h时，第k个控制点的层速度；k表示表层厚度为h时，第k个控制点所对应的表层速度；n表示表层厚度为h时，第k个控制点所对应的表层速度数量；当表层速度相同时，第i个目标面波频散点所对应的层速度数量与第k个控制点所对应的表层速度数量相同。
S150：将最高的相关系数所对应的半波长系数作为所述目标面波频散点在相应表层厚度的半波长系数。
在实际应用中，深度模型越接近于实际地层，其与实际地层相关性越好，反之，其与实际地层相关性越差。因此，当已知实际采集得到的控制点的转换得到的地层速度和厚度信息时，通过分析不同半波长系数对应的层速度、表层深度与实际地层控制点所对应的表层速度、表层深度之间的相关性，可以确定符合实际地层情况的半波长系数。因此，本申请实施例中将最高的相关系数所对应的半波长系数作为所述目标面波频散点在相应表层厚度的半波长系数。
如图3所示，图3是现有技术中实际测量的微测井以及现有固定半波长系数的表层厚度与表层速度曲线的对比图。图中310是现有固定半波长系数的表层厚度与表层速度曲线，320是实际测量的微测井的表层厚度与表层速度曲线。从图中可见，两者存在较大差异。
如图4所示，图4是现有技术中实际测量的微测井以及本申请确定的半波长系数所对应的表层厚度与表层速度曲线的对比图。图中410是本申请确定的半波长系数所对应的表层厚度与表层速度曲线，420是实际测量的微测井的表层厚度与表层速度曲线。从图4中可见，两者之间的差异较小。
在另一些实施例中，由于目标面波频散点的数量较大，本申请实施例中，还可以对目标面波频散点进行分段后依次进行计算，继而可以减少计算量，相应的，在计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数之前，所述方法还可以包括：
对所述多个目标面波频散点进行分段处理，第四预设数量的分段目标面波频散点；
相应的，所述计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数包括计算同一表层厚度所对应的控制点的表层速度与分段目标面波频散点的层速度之间的相关系数。
在另一些实施例中，所述方法还包括：
基于所述目标面波频散点在相应表层厚度的半波长系数进行半波长法近地表结构反演，得到所述目标面波频散点近地表结构；
利用所述目标面波频散点近地表结构进行建模，得到地层模型；
基于所述地层模型确定近地表结构反演数据。
由以上本申请一种数据处理方法的实施例可见，本申请通过计算设置的第一预设数量的半波长系数所对应的目标面波频散点在同一表层厚度所对应的层速度与实际地层所对应的控制点的表层速度之间的相关性，来确定出最符合地层实际情况的半波长系数。与现有技术相比，利用本申请提供的技术方案可以准确的确定出符合地层实际情况的半波长系数，进而可以提高后续地层模型，和近地表结构反演数据的精度。
本申请另一方面还提供一种数据处理装置，图5是本申请提供的数据处理装置的一种实施例中的结构示意图；如图5所示，所述装置500可以包括：
系数设置模块510，可以用于基于预设半波长系数范围设置多个目标面波频散点第一预设数量的半波长系数；
第一计算模块520，可以用于利用所述目标面波频散点的层速度、频率和半波长系数计算所述目标面波频散点的所述半波长系数所对应的表层厚度；
表层速度扩展处理模块530，可以用于基于预设控制点表层速度和表层厚度按照预设深度内插间隔进行表层速度扩展处理，得到第二预设数量的控制点的表层数据，所述控制点的表层数据包括一一对应的表层速度和表层厚度；
第二计算模块540，可以用于计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数；
半波长系数确定模块550，可以用于将最高的相关系数所对应的半波长系数作为所述目标面波频散点在相应表层厚度的半波长系数。
另一实施例中，所述系数设置模块510可以包括：
第一系数设置单元，可以用于利用所述预设半波长系数范围和所述目标面波频散点的数量确定半波长系数最小值，半波长增量和半波长系数最大值；
第一半波长系数组确定单元，可以用于基于所述半波长系数最小值，所述半波长增量和所述半波长系数最大值确定多个目标面波频散点的第一半波长系数组；
第一替换系数选取单元，可以用于从所述第一半波长系数组中选取一个半波长系数，作为替换系数；
第一待替换系数组确定单元，可以用于将所述第一半波长系数组作为待替换系数组；
第一替换单元，用于以所述替换系数替换所述待替换半波长系数组中任一与所述替换系数不同的半波长系数，得到所述多个目标面波频散点的第一数量的第二半波长系数组；
第二替换单元，可以用于以所述替换系数替换所述第一数量的第二半波长系数组中与所述替换系数不同的半波长系数，得到所述多个目标面波频散点的第二数量的第三半波长系数组；
第一数据处理单元，可以用于重复替换半波长系数的操作至得到的半波长系数组中的半波长系数都为所述替换系数；
第二待替换系数组确定单元，可以用于将当前获得的半波长系数组作为待替换半波长系数组；
第二替换系数选取单元，可以用于选取所述第一半波长系数组中新的半波长系数作为替换系数，所述新的半波长系数为第一半波长系数组中未作为替换系数的半波长系数；
第二数据处理单元，可以用于重复上述替换处理的步骤至预设收敛条件；
半波长系数确定单元，可以用于将所述第一半波长系数组和达到所述预设收敛条件之前每次替换得到的半波长系数组中的半波长系数作为目标面波频散点第一预设数量的半波长系数；
其中，所述预设收敛条件包括选取所述第一半波长系数组中的全部半波长系数作为替换系数，且最后一个替换系数替换得到半波长系数组中的全部半波长系数为所述最后替换系数。
另一实施例中，所述第一计算模块520中采用的计算公式可以包括下述计算公式：
Hij＝Cij*Vi/fi
其中，所述H表示第i个目标面波频散点的第j个半波长系数所对应的表层厚度；所述Cij表示第i个目标面波频散点的第j个半波长系数；Vi表示第i个目标面波频散点的层速度；fi表示第i个目标面波频散点的频率。
另一实施例中，所述第二计算模块540中采用的计算公式可以包括下述计算公式：
其中，rhi表示表层厚度为h时，多组模板表层数据中表层速度与第i个目标面波频散点目标面波频散点的层速度之间的相关系数；xhik表示表层厚度为h时，第i个目标面波频散点的第k个半波长系数所对应的层速度；k表示表层厚度为h时，第i个目标面波频散点所对应的第k个表层速度；n表示表层厚度为h时，第i个目标面波频散点所对应的层速度数量；yhk表层厚度为h时，第k个控制点的层速度；k表示表层厚度为h时，第k个控制点所对应的表层速度；n表示表层厚度为h时，第k个控制点所对应的表层速度数量；当表层速度相同时，第i个目标面波频散点所对应的层速度数量与第k个控制点所对应的表层速度数量相同。
另一实施例中，所述装置500还可以包括：
分段处理模块，可以用于在计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数之前，对所述多个目标面波频散点进行分段处理，第四预设数量的分段目标面波频散点；
相应的，所述第二计算模块540计算同一表层厚度所对应的控制点的表层速度与目标面波频散点的层速度之间的相关系数包括计算同一表层厚度所对应的控制点的表层速度与分段目标面波频散点的层速度之间的相关系数。
另一实施例中，所述装置500还可以包括：
半波长法反演模块，可以用于基于所述目标面波频散点在相应表层厚度的半波长系数进行半波长法近地表结构反演，得到所述目标面波频散点近地表结构；
地层模型建立模块，可以用于利用所述目标面波频散点近地表结构进行建模，得到地层模型；
近地表结构反演数据确定模块，可以用于基于所述地层模型确定近地表结构反演数据。
由以上本申请一种数据处理方法或装置的实施例可见，本申请通过计算设置的第一预设数量的半波长系数所对应的目标面波频散点在同一表层厚度所对应的层速度与实际地层所对应的控制点的表层速度之间的相关性，来确定出最符合地层实际情况的半波长系数。与现有技术相比，利用本申请提供的技术方案可以准确的确定出符合地层实际情况的半波长系数，进而可以提高后续地层模型，和近地表结构反演数据的精度。
本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。