一种确定地质界面的方法
【技术领域】
[0001] 本申请涉及地球物理勘探技术领域，特别涉及一种确定地质界面的方法。
【背景技术】
[0002] 目前，地震勘探法是在陆地和海洋勘探石油和天然气的重要手段。在地震勘探的 过程中，为了建立各个地层的速度模型，通常需要在地震剖面中确定不同地层之间的地质 界面。在确定地质界面后，当对地层速度模型中的速度进行更新后，一般需要重新确定地震 剖面，从而需要重新确定地质界面，以获取新的地质界面，确保地质界面的地质意义。
[0003] 现有技术中，一般利用射线追踪法重新确定地质界面，具体如下：
[0004] 首先，选取速度更新前地震剖面中地质界面上的一点，沿着该点的法线方向发射 射线，记录射线的出射时间t及出射角Θ ;然后，在出射点处以角度Θ反方向发射射线，采 用更新后的速度在时间t内前进的距离作为新的地质界面上的一点；然后，依次选取速度 更新前地震剖面中地质界面上的每个点，重复以上的步骤；最后，对得到的新的地质界面上 的点进行拟合，从而得到新的地质界面。
[0005] 上述现有技术中重新确定地质界面的方法，当地质界面较为复杂时，由于难以准 确的确定射线的出射时间t及出射角Θ，从而难以准确地拟合出与速度更新后的地震剖面 相一致的地质界面。
【发明内容】
[0006] 本申请实施例的目的是提供一种确定地质界面的方法，以在对地层速度模型中的 速度进行更新后，较准确地拟合出与速度更新后的地震剖面相一致的地质界面。
[0007] 为解决上述技术问题，本申请实施例提供一种确定地质界面的方法是这样实现 的：
[0008] 一种确定地质界面的方法，包括：
[0009] 1)以第一点为中心点，从第一地震剖面内获取水平方向长度为Lx、竖直方向长度 为L y的第一图像块，所述第一点为从第一地震剖面内的地质界面上选取的点；
[0010] 2)以第二地震剖面内的第二点为中心点，从第二地震剖面内获取水平方向长度为 Lx、竖直方向长度为My的第二图像块，所述第二点在第二剖面内的位置与所述第一点在第 一剖面内的位置相同，所述第二地震剖面为第一地震剖面内地层的速度模型更新后的地震 剖面，其中，M y> Ly;
[0011] 3)从第二图像块内获取My-Ly+1个水平方向长度为L x、竖直方向长度为Ly的子图 像块，从My-Ly+1个子图像块中选取子图像块，获取选取的子图像块的中心点与第二点之间 的位置关系，并建立第一点与该位置关系的对应关系；
[0012] 4)依据预设规则选取第一地震剖面内地质界面上的点，针对选取的每个点，重复 步骤1)至步骤3);
[0013] 5)根据所述位置关系和所述对应关系，确定第二地震剖面内的地质界面。
[0014] 由以上本申请实施例提供的技术方案可见，本申请实施例通过对比速度更新前后 地震剖面之间图像块的相似程度，来确定速度更新前后地震剖面内地质界面上点的位置关 系，并依据所述位置关系计算速度更新后地震剖面内地质界面上的点，最后通过曲线拟合 得到速度更新后地震剖面内的地质界面。与现有技术相比，本申请实施例的方法不需要准 确的确定射线的出射时间t及出射角Θ，因而可以高精度的拟合出与速度更新后的地震剖 面相一致的地质界面。
【附图说明】
[0015] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可以根据这些附图获得其他的附图。
[0016] 图1为本申请实施例一个地震剖面图；
[0017] 图2为本申请实施例一种确定地质界面方法的流程图；
[0018] 图3为本申请实施例一种根据位置关系和对应关系确定地震剖面内地质界面方 法的流程图。
【具体实施方式】
[0019] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护 的范围。
[0020] -般地，在地震剖面内，各个地层之间的分界面称为地震剖面内的地质界面。在地 震勘探的过程中，在通过地震数据得到地震剖面后，可以在地震剖面内确定不同地层之间 的地质界面。附图1示出了某个地震剖面，在该地震剖面内，5条黑色的线条表示5个地质 界面。线条上的小黑点为地质界面上的点。该5个地质界面可以将该地震剖面分为6个地 层。对于每个地层，可以建立该地层的速度模型。
[0021] 在对地层建立速度模型的过程中，可能需要对速度模型进行多次调整，以对速度 模型中的速度进行更新。当通过某个地震剖面建立速度模型后，若需要对速度模型中的速 度进行更新，相应地，该地震剖面也需要进行更新。对于速度更新后的地震剖面，需要重新 确定该地震剖面内的地质界面，以确保地质界面的地质意义。下面介绍本申请确定地质界 面方法的一个实施例。如图2所示，该实施例包括：
[0022] S201 :以第一点为中心点，从第一地震剖面内获取水平方向长度为Lx、竖直方向长 度为Ly的第一图像块，所述第一点为从第一地震剖面内的地质界面上选取的点。
[0023] 所述第一地震剖面一般为速度更新前的地震剖面。
[0024] -般地，一个地震剖面内可以有多个地质界面。在步骤S201中，可以从第一地震 剖面内选取一个地质界面，然后从该地质界面上选取一点，将选取的点记为第一点。
[0025] 进一步地，对于第一地震剖面，可以以该地震剖面的左上角为原点，以水平向右的 方向为X轴的正方向，以竖直向下的方向为Y轴的正方向，建立第一直角坐标系。所述第 一点在第一直角坐标系内的坐标值可以记为（Px，Py)。以第一点为中心，可以从第一地震剖 面内得获取大小为L xXLy的第一图像块。其中，L#第一图像块沿X轴方向的长度，LyS 第一图像块沿Y轴方向的长度。所述1^与Ly均为正整数。例如，第一点的坐标值可以为 (20, 100)。Lx的值可以为8。Ly的值可以为24。
[0026] S202 :以第二点为中心点，从第二地震剖面内获取水平方向长度为Lx、竖直方向长 度为My的第二图像块，所述第二点为第二地震剖面内的点，并且所述第二点在第二剖面内 的位置与所述第一点在第一剖面内的位置相同，所述第二地震剖面为所述第一地震剖面内 地层的速度模型更新后的地震剖面，其中，M y> Ly。
[0027] 所述第二地震剖面一般为速度更新后的地震剖面。根据在地震勘探过程中获取的 地震数据，可以通过现有的偏移方法得到速度更新后的地震剖面，即第二地震剖面。
[0028] 进一步地，对于第二地震剖面，可以以该地震剖面的左上角为原点，以水平向右的 方向为X轴的正方向，以竖直向下的方向为Y轴的正方向，建立第二直角坐标系。在第二直 角坐标系内，可以选取坐标为（P x，Py)的点，将该点记为第二点。以第二点为中心，可以从第 二地震剖面内获取大小为L xXMy的第二图像块。其中，LxS第二图像块沿X轴方向的长度， M#第二图像块沿Y轴方向的长度，My> Ly。所述LAMy均为正整数。例如，若第一点在 第一直角坐标系的坐标值为（20, 100)，Lx的值为8,Ly的值为24。那么，第二点在第二直角 坐标系的坐标值可以为（20, 100)，My的值可以为50。
[0029] S203 :从第二图像块内获取My-Ly+1个水平方向长度为Lx、竖直方向长度为L y的子 图像块，从My-Ly+1个子图像块中选取子图像块，获取选取的子图像块的中心点与第二点之 间的位置关系，并建立第一点与该位置关系的对应关系。
[0030] 具体地，从第二图像块内获取乂-Ly+l个水平方向长度为Lx、竖直方向长度为1^的 子图像块，分别将每个子图像块与第一图像块进行对比，以从M y-Ly+1个子图像块中选取子 图像块，获取选取的子图像块的中心点与第二点之间的位置关系，并建立第一点与该位置 关系的对应关系。所述中心点一般为第二地震剖面地质界面上的点。所述位置关系可以为 位置差值，例如可以为中心点与第二点在第二直角坐标系中的坐标差值。
[0031] 可选地，沿竖直方向从第二图像块内获取My_Ly+l个水平方向长度为L x、竖直方向 长度为Ly的子图像块，相邻子图像块的上边沿或下边沿之间的间隔为1。例如，若M y的值 为50, Lx的值为8, L y的值为24。那么，可以从第二图像块内获取27个水平方向长度为8、 竖直方向长度为24的子图像块。
[0032] 可选地，可以获取My_Ly+l个子图像块。针对每个子图像块，计算该子图像块与第 一图像块的相似度。根据每个子图像块的相似度，从M y-Ly+1个子图像块中获取与第一图像 块最相似的子图像块，将获取的子图像块记为第三子图像块。
[0033] 具体地，可以通过公式￡(4-珲)2计算子图像块与第一图像块的相似度。其中，N /=1 为子图像块与第一图像块的像素点的个数；由于子图像块与第一图像块的大小相同，所以 子图像块与第一图像块的像素点的个数相同。i为像素点。化