转换波井震标定方法及装置的制造方法
【技术领域】
[0001]本申请涉及地震勘探技术领域，尤其是涉及一种转换波井震标定方法及装置。
【背景技术】
[0002]井震标定是地震解释中的一个关键的、重要的步骤，是连接测井、地震与地质信息的桥梁，其结果准确与否直接决定着地震解释和反演的准确性。以纵波为例，井震标定的过程一般包括以下几个步骤:首先基于井曲线计算反射系数；其次将地震子波与反射系数进行褶积产生合成记录；然后将合成记录与地震道进行时移偏差计算；最后更新时深关系曲线。
[0003]而转换波资料的信噪比通常比纵波低，且地层对转换波能量的吸收比纵波强，使得转换波的能量较小、频率较低。在转换波地震数据和转换波合成记录中存在的这些不利因素，使得转换波井震标定比纵波井震标定难度大得多。与纵波类似，转换波井震标定中也需要对合成记录进行压缩和拉伸。目前常用的标定软件提供了一个交互界面以便于进行测井和地震之间的人工标定，即允许用户通过交互进行时移、拉伸和压缩合成地震记录的方式来匹配地震道。当用户手动进行时移偏差计算合成记录和地震道时，标定软件会自动更新时深曲线。然而，这个过程是费时费力的，标定结果完全取决于用户的技能和经验，难免会有较大的人为误差，因而其标定结果的可靠性较低，从而不利于后续地震解释和反演结果的可信度和准确性。
【发明内容】
[0004]本申请实施例的目的在于提供一种转换波井震标定方法及装置，以提高转换波井震标定的可靠性。
[0005]为达到上述目的，一方面，本申请实施例提供了一种转换波井震标定方法，包括以下步骤:
[0006]获取井旁地震道的转换波合成记录；
[0007]确定所述转换波合成记录的整体时移量和相位旋转量，并分别基于所述整体时移量和所述相位旋转量对所述转换波合成记录进行整体时移和相位旋转，获得初优化转换波合成记录；
[0008]基于预设算法对所述初优化转换波合成记录和所述井旁地震道进行时移偏差计算，获得所述初优化转换波合成记录中每个样点相对于所述井旁地震道中对应样点的时移量;
[0009]基于每个样点的时移量对所述初优化转换波合成记录进行时移调整，获得再优化转换波合成记录；
[0010]基于所述再优化转换波合成记录获取时深曲线。
[0011]另一方面，本申请实施例还提供了一种转换波井震标定装置，包括:
[0012]第一合成记录获取模块，用于获取井旁地震道的转换波合成记录；
[0013]第二合成记录获取模块，用于确定所述转换波合成记录的整体时移量和相位旋转量，并分别基于所述整体时移量和所述相位旋转量对所述转换波合成记录进行整体时移和相位旋转，获得初优化转换波合成记录；
[0014]样点时移量获取模块，用于基于预设算法对所述初优化转换波合成记录和所述井旁地震道进行时移偏差计算，获得所述初优化转换波合成记录中每个样点相对于所述井旁地震道中对应样点的时移量；
[0015]第三合成记录获取模块，用于基于每个样点的时移量对所述初优化转换波合成记录进行时移调整，获得再优化转换波合成记录；
[0016]时深曲线获取模块，用于基于所述再优化转换波合成记录获取时深曲线。
[0017]本申请实施例在基于整体时移量和相位旋转量对转换波合成记录进行整体时移和相位旋转，获得初优化转换波合成记录后，又用预设算法来匹配转换波合成记录和井旁地震道，该匹配过程能替换了现有井震标定中的手工拉伸和压缩，从而将人为因素导致的井震标定误差减小到最低，因而可极大地提高转换波井震标定的可靠性。从而有利于后续地震解释和反演结果的可信度和准确性。因此，本申请实施例在地震资料解释和储层预测方面具有广阔的应用前景。
【附图说明】
[0018]此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的一部分，并不构成对本申请实施例的限定。在附图中:
[0019]图1为本申请实施例的转换波井震标定方法的流程图；
[0020]图2为本申请实施例中转换波叠加数据示意图；
[0021]图3为本申请实施例中选取的雷克子波的波形示意图；
[0022]图4为本申请实施例的测井数据示意图；
[0023]图5为本申请实施例中合成记录的相位旋转20度时的示意图；
[0024]图6为本申请实施例中转换波数据与合成记录的对比示意图；
[0025]图7为采用本申请实施例更新前后的时深关系曲线示意图(虚线为更新前，实线为更新后)；
[0026]图8为本申请实施例的转换波井震标定装置的结构框图；
[0027]图9为本申请实施例中第二合成记录获取模块的结构框图。
【具体实施方式】
[0028]为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于解释本申请实施例，但并不作为对本申请实施例的限定。
[0029]下面结合附图，对本申请实施例的【具体实施方式】作进一步的详细说明。
[0030]参考图1所示，本申请实施例的转换波井震标定方法包括以下步骤:
[0031]S1、获取井旁地震道的转换波合成记录。具体获取步骤如下:
[0032]首先，在野外利用纵波震源激发地震波并利用检波器记录地震波，按照常规地震资料处理流程对采集的转换波数据进行相对振幅保持的高保真处理，形成转换波叠加数据，如图2所示；
[0033]其次，通过对井旁地震道进行振幅谱分析，例如图2中的道号为332的地震道(以下简称道)获得振幅谱的中心频谱为25Hz。
[0034]然后，以25Hz中心频谱作为子波的峰值频率，产生25Hz的雷克子波，如图3所示；
[0035]再次，根据测井数据(例如图4所示)计算设定转换波入射角(例如20度)的转换波反射系数；
[0036]最后，将转换波反射系数和雷克子波进行褶积可得到转换波合成记录。
[0037]S2、确定所述转换波合成记录的整体时移量和相位旋转量，并分别基于所述整体时移量和所述相位旋转量对所述转换波合成记录进行整体时移和相位旋转，获得初优化转换波合成记录。所述转换波合成记录的整体时移量和相位旋转量是通过对所述井旁地震道和所述转换波合成记录进行相关分析得到的。具体的，
[0038]本申请实施例中，所述转换波合成记录的整体时移量具体可通过以下步骤得到:
[0039]将转换波合成记录按照预设的时移顺序(例如由上至下)和时移步进量进行整体时移，并计算各次整体时移后的转换波合成记录与所述井旁地震道的互相关系数；
[0040]选取其中互相关系数最大时的时移增量作为所述转换波合成记录的整体时移量。
[0041]本申请实施例中，所述转换波合成记录的相位旋转量通过具体以下步骤得到:
[0042]将基于所述整体时移量整体时移后的转换波合成记录的相位，按照预设的旋转顺序和旋转步进量(例如I度)进行旋转，并计算各次旋转后的转换波合成记录与所述井旁地震道的互相关系数；
[0043]选取其中互相关系数最大时的相位增量作为所述转换波合成记录的相位旋转量。一般的，转换波合成记录的相位的变化范围为I?360度，例如经过计算分析，当转换波合成记录的相位旋转20度时，转换波合成记录与井旁地震道的互相关系数最大，这时对转换波合成记录进行20度的常相位旋转(例如图5所示)。
[0044]本申请实施例中，上述整体时移和相位旋转可先后进行，比如先计算出整体时移量并以该整体时移量为依据进行整体时移，然后计算相位旋转量并以该相位旋转量为依据进行相位旋转。当然，在本申请的另一实施例中也可以采用其他顺序。
[0045]S3、基于预设算法对所述初优化转换波合成记录和所述井旁地震道进行时移偏差计算，获得所述初优化转换波合成记录中每个样点相对于所述井旁地震道中对应样点的时移量。所述预设算法可以采用动态时间规整算法(DTW，Dynamic Time Warping)。本申请实施例中DTW算法的一个具体实现可以是:设定一个时移量范围，依此参数DTW计算初优化转换波合成记录和井旁地震道这两个时间序列之间的二维对齐误差，进而得到二维累积距离，在二维累积距离中反向追踪寻找最小路径，沿此路径可得到每个样点的时移量。
[0046]S4、基于每个样点的时移量对所述初优化转换波合成记录进行时移调整，获得再优化转换波合成记录。
[0047]S5、基于所述再优化转换波合成记录获取时深曲线。结合图7所示，更新前，转换波合成记录和井旁地震道之间的相关系数小，由此得到的时深关系曲线可靠性低；更新后，转换波合成记录和井旁地震道之间的相关系数大，由此得到的时深关系曲线可靠性高。
[0048]此外，本申请实施例中，为了更好地将数据进行显示和对比，同一数据可进行重复显示。例如图6中所示的第16?20道和第21?25道数据；从第一道开始，每5道是同一数据重复了 5次。其中第I?5道是转换波井旁地震道，第6?10道是转换波合成记录，其与第1-5道的相关系数为0.57 ;第11?15道是进行相位旋转后的转换波合成记录，其与第1-5道的相关系数为0.60 ;第16-20道是对第11?15道进行简单动态调整后的合成记录，其与第1-5道的相关系数为0.86 ;第21?25道是对第11?15道进行优化动态调整后的合成记录，其与第I?5道的相关系数为0.92。
[0049]本申请实施例在进行初步的调整后，用动态时间规整算法来匹配转换波合成记录和井旁地震道，该匹配过程能替换了现有井震标定中的手工拉伸和压缩，从而将人为因素导致的井震标定误差减小到最低，因而可极大地提高转换波井震标定的可靠性。从而有利于后续地震解释和反演结果的可信度和准确性。因此，本申请实施例在地震资料解释和储层预测方面具有广阔的应用前景。
[0050]结合图8所示，与上述转换波井震标定方法实施例对应，本申请实施例的转换波井震标定装置包括:
[0051]第一合成记录获取模块81，用于获取井旁地震道的转换波合成记