用于确定断层展布格局的地震资料解释方法和设备的制造方法
【技术领域】
[0001]本发明涉及地震资料解释领域，具体地，涉及一种用于确定断层展布格局的地震资料解释方法和设备。
【背景技术】
[0002]我国大陆经历多次板块构造旋回的叠加，地质构造复杂，尤其是断层构造发育。在能源、工程勘探中，首先需要对断层构造展布进行确定，以便科学确定能源开发、工程布置方案。地震资料解释是用于确定构造展布的一种重要方法。目前，常规的断层展布地震解释方法主要是:首先对目的层反射波进行层位追踪，然后通过手工解释地震剖面，画出断面形状，层位与断面的交点进行平面投影，通过多个交点圈出断层，获得断层展布格局。这种方法需要较大的人工解释工作，而且属于一种二维的地震解释方法，只能对较大的构造有较好效果，工作效率较低。为了提高断层构造展布格局的解释效率，一般解释方法是通过地震属性提取水平切片或者沿层切片。该方法需要两个大步骤，一是地震属性的计算，二是目的层地震反射波追踪。地震属性体(例如，方差体、地震振幅等参数)的计算往往涉及较大的计算量，对于普通的计算机，则需要较长的时间，在高性能的机器才能获得较短的计算时间。同时，该属性体的计算结果主要是帮助人工识别断点，在进行目的层地震反射波追踪中，主要是人工解释为主，因此在实际应用中存在一定的局限性。
【发明内容】
[0003]本发明的目的是提供一种能够快速确定断层展布格局的地震资料解释方法和设备。
[0004]为了实现上述目的，本发明提供一种用于确定断层展布格局的地震资料解释方法。该方法包括:步骤S1:对目的层反射波进行全三维层位自动追踪，获得原始层位信息；步骤S2:对所述原始层位信息进行平滑，获得平滑层位信息；步骤S3:基于所述原始层位信息与所述平滑层位信息之间的差，得出反射波时间差信息；以及步骤S4:确定所述反射波时间差信息的各个波峰极值点和各个波谷极值点，并对所述各个波峰极值点进行拟合以及对所述各个波谷极值点进行拟合，得出断层交线。
[0005]优选地，采用互相关法对所述目的层反射波进行所述全三维层位自动追踪，获得所述原始层位信息。
[0006]优选地,所述步骤SI包括:步骤SlOl:在一地震剖面上,确定所述目的层反射波在基准地震道上的位置，该位置作为种子点；步骤S102:设定一时窗；步骤S103:在所述地震剖面上的目标地震道上，按照设定的采样间隔移动所述时窗，并确定对应的相关系数，其中，该相关系数表示在所述时窗未移动时所述基准地震道上被包括在该时窗内的波形与在所述时窗移动之后所述目标地震道上被包括在该时窗内的波形之间的相关性；步骤S104:判断所述相关系数是否满足设定的阈值条件，在所述相关系数不满足所述阈值条件的情况下，返回步骤S103继续进行，直到所确定出的相关系数满足所述阈值条件为止；在所述相关系数满足所述阈值条件的情况下，进行步骤S105:基于所述采样间隔和与该相关系数对应的移动次数，确定所述种子点在所述目标地震道上的移动量，以确定所述目的层反射波在所述目标地震道上的位置；步骤S106:判断是否已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置，若是，进行步骤S108:基于所述目的层反射波在所述地震剖面上的各个地震道上的位置，得出所述原始层位信息；若不是，进行步骤S107:更换目标地震道，并复位所述时窗，之后返回步骤S103继续进行，直到确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置为止。
[0007]优选地，所述步骤SI包括:步骤SlOl:在一地震剖面上，确定目的层反射波在基准地震道上的位置，该位置作为种子点；步骤S102:设定一时窗；步骤S103:在所述地震剖面上的目标地震道上，按照设定的采样间隔移动所述时窗，并确定对应的相关系数，其中，该相关系数表示在所述时窗未移动时所述基准地震道上被包括在该时窗内的波形与在所述时窗移动之后所述目标地震道上被包括在该时窗内的波形之间的相关性；步骤S104:判断所述时窗在所述目标地震道上的移动次数是否已达到最大移动次数，其中，该最大移动次数是基于设定的最大倾角时差和所述采样间隔来确定的，在所述移动次数未达到所述最大移动次数的情况下，返回步骤S103继续进行，直到所述移动次数达到所述最大移动次数为止；在所述移动次数达到所述最大移动次数的情况下，进行步骤S105:基于所述采样间隔和与最大相关系数对应的移动次数，确定所述种子点在所述目标地震道上的移动量，以确定所述目的层反射波在所述目标地震道上的位置；步骤S106:判断是否已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置，若是，进行步骤S108:基于所述目的层反射波在所述地震剖面上的各个地震道上的位置，得出所述原始层位信息；若不是，进行步骤S107:更换目标地震道，并复位所述时窗，之后返回步骤S103继续进行，直到确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置为止。
[0008]优选地，该方法还包括:在所述移动次数达到所述最大移动次数的情况下，先进行步骤S109:判断所述最大相关系数是否满足设定的阈值条件，在所述最大相关系数满足所述阈值条件的情况下，再进行步骤S105。
[0009]优选地，采用最小二乘平滑法对所述原始层位信息进行平滑，获得所述平滑层位信息。
[0010]优选地，采用以下方法之一对所述各个波峰极值点进行拟合以及对所述各个波谷极值点进行拟合:最小二乘拟合法或蚂蚁算法。
[0011]本发明还提供一种用于确定断层展布格局的地震资料解释设备。该设备包括:用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置；用于对所述原始层位信息进行平滑，获得平滑层位信息的装置；用于基于所述原始层位信息与所述平滑层位信息之间的差，得出反射波时间差信息的装置；以及用于确定所述反射波时间差信息的各个波峰极值点和各个波谷极值点，并对所述各个波峰极值点进行拟合以及对所述各个波谷极值点进行拟合，得出断层交线的装置。
[0012]优选地，对所述目的层反射波进行所述全三维层位自动追踪，获得所述原始层位信息是采用互相关法进行的。
[0013]优选地，所述用于对目的层反射波进行全三维层位自动追踪，获得原始层位信息的装置包括:用于在一地震剖面上，确定所述目的层反射波在基准地震道上的位置的第一装置，该位置作为种子点；用于设定一时窗的第二装置；用于在所述地震剖面上的目标地震道上，按照设定的采样间隔移动所述时窗，并确定对应的相关系数的第三装置，其中，该相关系数表示在所述时窗未移动时所述基准地震道上被包括在该时窗内的波形与在所述时窗移动之后所述目标地震道上被包括在该时窗内的波形之间的相关性；用于判断所述相关系数是否满足设定的阈值条件的第四装置；用于在所述相关系数不满足所述阈值条件的情况下，重复运行所述第三装置，直到所确定出的相关系数满足所述阈值条件为止的第五装置；用于在所述相关系数满足所述阈值条件的情况下，基于所述采样间隔和与该相关系数对应的移动次数，确定所述种子点在所述目标地震道上的移动量，以确定所述目的层反射波在所述目标地震道上的位置的第六装置；用于判断是否已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的第七装置；用于在尚未确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，更换目标地震道，并复位所述时窗，之后重复运行所述第三装置，直到确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置为止的第八装置；以及用于在已确定出所述目的层反射波在所述地震剖面上的所有地震道上的位置的情况下，基于所述目的层反射波在所述地震剖面上的各个地震道上的位置，得出所述原始层位信息的第九装置。
[0014]优选地，所述用