所述目的层反射波在所述基准地震道上的位置(即，重新确定种子点)、重新设定所述时窗7的大小、重新设定最大倾角时差、重新设定所述采样间隔或重新设定所述阈值条件。也就是说，在时窗7已经移动了最大移动次数之后仍没有对应的相关系数能够满足阈值条件，那么，可以重新确定种子点6和/或调整一些参数。之后，基于新的设定，按照上述同样的方法，重新在地震剖面上的各个目标地震道对种子点6进行自动追踪。
[0054]以上描述了三种示例方式，无论采样哪种方式，均可以实现对目的层反射波的全三维层位自动追踪，并获得原始层位信息。
[0055]图6中示出了根据本发明提供的方法进行层位自动追踪之后得到的原始层位信息的一个示例。如图6所示，曲线17表示含有断层构造的煤层，曲线14为对煤层反射波的原始层位追踪结果，即所述原始层位信息。
[0056]在获得原始层位信息之后，参考图1中的步骤S2，可以对所述原始层位信息进行平滑，获得平滑层位信息。对所述原始层位信息进行平滑的方法有很多，在本发明中，采用最小二乘平滑法对所述原始层位信息进行平滑。最小二乘平滑法对于本领域的技术人员是公知的，因而，本发明在此不进行赘述。以图6的示例为例，原始层位信息14经平滑后，得到平滑层位信息15。从平滑层位信息15可以看出，经过平滑后，在断层位置，反射波更为平滑，从而与原始层位信息14出现了一个时间差。
[0057]之后，步骤S3，基于所述原始层位信息与所述平滑层位信息之间的差，得出反射波时间差信息。在本发明中，可以将平滑层位信息15减去原始层位信息14，得到反射波时间差信息16(反之亦然，即，也可以利用原始层位信息减去平滑层位信息来得出反射波时间差信息)。在反射波时间差信息16上能明显看出断层带表现为一对波峰和波谷，分别对应时间差明显大值和时间差明显小值。需要注意的是，在断层带对应的时间差上，不一定表现为一对完好的波峰和波谷，也可能表现为一对近似的波峰和波谷。
[0058]之后，步骤S4，确定所述反射波时间差信息的各个波峰极值点和各个波谷极值点，并对所述各个波峰极值点进行拟合以及对所述各个波谷极值点进行拟合，得出断层交线。
[0059]可以通过对所述反射波时间差信息求解一阶导数，来提取所述各个波峰极值点和各个波谷极值点。之后，将所述各个波峰极值点和各个波谷极值点映射到以横测线19 (cdp方向)为横坐标、纵测线18 (inline方向)为纵坐标的坐标系中。
[0060]虽然图6中示出的反射波时间差信息16仅具有一对波峰和波谷，但这仅是局部的示意图，在实际中，反射波时间差信息16上会出现多对波峰和波谷，由此可以提取出多个波峰极值点和多个波谷极值点，并将它们映射到所述坐标系中，如图7所示。
[0061]之后，可以利用拟合算法，对各个波峰极值点进行拟合以及对各个波谷极值点进行拟合，从而得出断层交线。在本发明中，可以采用最小二乘拟合法或蚂蚁算法来进行拟合。拟合得出的波峰断层交线如图7中的曲线20所示，拟合得出的波谷断层交线如曲线21所示。需要说明的是，最小二乘拟合法和蚂蚁算法均是本领域技术人员公知的拟合算法，因此，其具体原理在本发明中不进行描述。
[0062]由此，通过本发明提供的地震资料解释方法，能够快速、准确地识别出断层展布格局，并生成断层交线。相较于现有的解释方法，大约能缩减30% -60%的时间，并且能够广泛地应用在各种实际应用中。
[0063]本发明还提供一种用于确定断层展布格局的地震资料解释设备。该设备可以包括:用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置；用于对所述原始层位信息进行平滑，获得平滑层位信息的装置；用于基于所述原始层位信息与所述平滑层位信息之间的差，得出反射波时间差信息的装置；以及用于确定所述反射波时间差信息的各个波峰极值点和各个波谷极值点，并对所述各个波峰极值点进行拟合以及对所述各个波谷极值点进行拟合，得出断层交线的装置。
[0064]其中，对所述目的层反射波进行所述全三维层位自动追踪，获得所述原始层位信息可以采用互相关法来进行。
[0065]在这种情况下，用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置可以包括:用于在一地震剖面上，确定所述目的层反射波在基准地震道上的位置的第一装置，该位置作为种子点；用于设定一时窗的第二装置；用于在所述地震剖面上的目标地震道上，按照设定的采样间隔移动所述时窗，并确定对应的相关系数的第三装置，其中，该相关系数表示在所述时窗未移动时所述基准地震道上被包括在该时窗内的波形与在所述时窗移动之后所述目标地震道上被包括在该时窗内的波形之间的相关性；用于判断所述相关系数是否满足设定的阈值条件的第四装置；用于在所述相关系数不满足所述阈值条件的情况下，重复运行所述第三装置，直到所确定出的相关系数满足所述阈值条件为止的第五装置；用于在所述相关系数满足所述阈值条件的情况下，基于所述采样间隔和与该相关系数对应的移动次数，确定所述种子点在所述目标地震道上的移动量，以确定所述目的层反射波在所述目标地震道上的位置的第六装置；用于判断是否已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的第七装置；用于在尚未确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，更换目标地震道，并复位所述时窗，之后重复运行所述第三装置，直到确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置为止的第八装置；以及用于在已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，基于所述目的层反射波在所述地震剖面上的各个地震道上的位置，得出所述原始层位信息的第九装置。
[0066]可替换地，用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置可以包括:用于在一地震剖面上，确定目的层反射波在基准地震道上的位置的第一装置，该位置作为种子点；用于设定一时窗的第二装置；用于在所述地震剖面上的目标地震道上，按照设定的采样间隔移动所述时窗，并确定对应的相关系数的第三装置，其中，该相关系数表示在所述时窗未移动时所述基准地震道上被包括在该时窗内的波形与在所述时窗移动之后所述目标地震道上被包括在该时窗内的波形之间的相关性；用于判断所述时窗在所述目标地震道上的移动次数是否已达到最大移动次数的第四装置，其中，该最大移动次数是基于设定的最大倾角时差和所述采样间隔来确定的；用于在所述移动次数未达到所述最大移动次数的情况下，重复运行所述第三装置，直到所述移动次数达到所述最大移动次数为止的第五装置；用于在所述移动次数达到所述最大移动次数的情况下，基于所述采样间隔和与最大相关系数对应的移动次数，确定所述种子点在所述目标地震道上的移动量，以确定所述目的层反射波在所述目标地震道上的位置的第六装置；用于判断是否已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的第七装置；用于在尚未确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，更换目标地震道，并复位所述时窗，之后重复运行所述第三装置，直到确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置为止的第八装置；以及用于在已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，基于所述目的层反射波在所述地震剖面上的各个地震道上的位置，得出所述原始层位信息的第九装置。
[0067]在这种替换的实施方式中，用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置还可以包括:用于在所述移动次数达到所述最大移动次数的情况下，判断所述最大相关系数是否满足设定的阈值条件，在所述最大相关系数满足所述阈值条件的情况下，再运行所述第六装置的第十装置。
[0068]此外，对所述原始层位信息进行平滑，获得所述