一种地层倾角检测方法
【技术领域】
[0001] 本发明涉及油气勘探技术领域，具体地说，涉及一种地层倾角检测方法。
【背景技术】
[0002] 地震数据处理中深度域偏移剖面的倾角、方位角是三维地下地质体的重要属性。 作为中间结果，它可以用来进行深度域网格层析反演，提高深度域速度模型精度。
[0003] 现有的地层倾角检测方法包括倾斜叠加法、平面波分解法、复地震道分析法、地震 波相干法、梯度结构张量法和给予加权结构张量法等方法。其中，利用倾斜叠加法能够计算 二维地震剖面同相轴的倾角；利用平面波分解法、复地震道分析法和地震波相干法能够获 取得到局部地层倾角；利用梯度结构张量法能够对三维地震数据的梯度向量进行分析，从 而得到局部地层倾角和方位角；利用给予加权结构张量法能够提高倾角估计的一致性。
[0004] 但是，现有的地层倾角检测方法在对低信噪比数据时的地层倾角进行检测时，会 出现检测结果不准确的问题。同时，当地层比较稀疏时，现有的地层倾角检测方法无法检测 缺乏反射波同相轴的地层的地层倾角。
[0005] 基于上述情况，亟需一种准确、可靠地进行地层倾角检测的方法。
【发明内容】
[0006] 为解决上述问题，本发明提供了一种地层倾角检测方法，所述方法包括：
[0007] 地震波形数据获取步骤，根据当前道分析点和各个预设扫描倾角，确定对应于所 述各个预设扫描倾角的临近道分析点，根据所述当前道分析点和临近道分析点，分别获取 当前道的地震波形数据和对应于所述各个预设扫描倾角的临近道的地震波形数据；
[0008] 相似能量谱计算步骤，根据所述当前道的地震波形数据和对应于所述各个预设扫 描倾角的临近道的地震波形数据，确定所述各个预设扫描倾角的相似能量谱；
[0009] 地层倾角确定步骤，从各个相似能量谱中提取最大相似能量谱，并将所述最大相 似能量谱所对应的预设扫描倾角作为地层倾角。
[0010] 根据本发明的一个实施例，在所述地震波形数据获取步骤中，
[0011] 在当前道中，获取以所述当前道分析点为中心的预设深度窗口的地震波形数据， 作为所述当前道的地震波形数据；
[0012] 在各个临近道中，获取以相应临近道分析点为中心的预设深度窗口的地震波形数 据，作为所述临近道的地震波形数据。
[0013] 根据本发明的一个实施例，所述当前道两侧的临近道的数量相等。
[0014] 根据本发明的一个实施例，所述相似能量谱计算步骤包括：
[0015] 根据所述当前道的地震波形数据和对应于所述各个预设扫描倾角的临近道的地 震波形数据，计算对应于所述各个预设扫描倾角的各个临近道的地震波形数据的相关值；
[0016] 根据所述各个临近道的地震波形数据的相关值，计算所述各个预设扫描倾角的相 似能量谱。
[0017] 根据本发明的一个实施例，根据如下公式计算所述各个临近道的地震波形数据的 相关值：
[0019] 其中，Cini表示对应于扫描倾角m的第i个临近道的地震波形数据的相关值，Wd表 示预设深度窗口所包含的分析点的总数，4表示所述当前道的预设深度窗口中第k个分析 点，b ink表示第i个临近道中对应于倾角m的预设深度窗口中的第k个分析点，Dak和Dbink 分别表示分析点ak和binik的地震波形数据。
[0020] 根据本发明的一个实施例，根据如下公式计算所述各个预设扫描倾角的相似能量 谱：
[0022] 其中，Sm表示扫描倾角m的相似能量谱，η表示临近道的总数。
[0023] 根据本发明的一个实施例，在所述地层倾角确定步骤中，还从各个相似能量谱中 提取与所述最大相似能量谱临近的若干相似能量谱，
[0024] 其中，根据所述最大相似能量谱和所述若干临近的能量谱确定地层倾角插值曲 线，并根据所述地层倾角插值曲线确定地层倾角。
[0025] 根据本发明的一个实施例，从各个相似能量谱中提取与所述最大相似能量谱相邻 的两个相似能量谱。
[0026] 根据本发明的一个实施例，在所述地层倾角确定步骤中，
[0027] 计算所述地层倾角插值曲线的一阶导数；
[0028] 将所述一阶导数等于零时所对应的倾角作为地层倾角。
[0029] 根据本发明的一个实施例，根据如下公式计算所述地层倾角：
[0031] 其中，X表示地层倾角，Y1表示最大相似能量谱，X1表示最大相似能量谱对应的预 设扫描倾角，^和y 2表示与所述最大相似能量谱相邻的两个相似能量谱。
[0032] 本发明所提供的地层倾角检测方法通过利用预设窗口长度的深度窗口来其他区 相应的地震波形数据，这样相较于现有的地层倾角检测方法，能够明显减少同相轴稀疏性 的影响，能够得到更为稳定、精度更高的倾角剖面。
[0033] 同时，通过相关属性的同相叠加，使得相较于现有的地层倾角检测方法，本方法对 数据信噪比的依赖性更低，并具有更好的通用性。本方法流程简单、易于实现。同时由于在 进行地层倾角检测时，是针对地震剖面上行的一个具体的点来计算的，与其他点之间没有 联系，因此本方法的计算密度高、并行性好，适合大规模并行计算。
[0034] 此外，本方法还可以根据最大相似能量谱及其取值邻近的相似能量谱来确定地层 倾角插值曲线，并根据该插值曲线来确定地层倾角。这样所得到的检测结果能够更加接近 真实值。而通过地层倾角插值曲线能够得到连续的地层倾角信息，还可以作为输入信息进 行深度域层析速度反演，具有更加广泛的应用价值。
[0035] 本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变 得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利 要求书以及附图中所特别指出的结构来实现和获得。
【附图说明】
[0036] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要的附图做简单的介绍：
[0037] 图1是根据本发明一个实施例的地层倾角检测方法的流程图；
[0038] 图2是根据本发明另一个实施例的地层倾角检测方法的流程图；
[0039] 图3是根据本发明一个实施例的输入的地震剖面图；
[0040] 图4是利用倾斜叠加法检测得到的图3的地层倾角剖面图；
[0041] 图5是根据本发明一个实施例的地层倾角检测方法得到的图3的地层倾角剖面 图。
【具体实施方式】
[0042] 以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用 技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明 的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合， 所形成的技术方案均在本发明的保护范围之内。
[0043] 同时，在以下说明中，出于解释的目的而阐述了许多具体细节，以提供对本发明实 施例的彻底理解。然而，对本领域的技术人员来说显而易见的是，本发明可以不用这里的具 体细节或者所描述的特定方式来实施。
[0044] 另外，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系 统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处 的顺序执行所示出或描述的步骤。
[0045] 实施例一:
[0046] 针对现有的地层倾角检测方法所存在的缺陷，本实施例基于窗口内地震波形的相 关值进行倾斜叠加分析，并根据叠加得到的能量最大值来确定最佳的地层倾角。
[0047] 图1示出了本实施例所提供的地层倾角检测方法的流程图。
[0048] 本实施例首先在地震波形数据获取步骤中获取当前道和对应于各个预设扫描倾 角的临近道的地震波形数据。具体地，如图1所示，首先在步骤SlOl中根据当前道分析点 和各个预设扫描倾角，确定对应于各个预设扫描倾角的临近道分析点。临近道分布在当前 道的两侧且靠近当前道，也可以说是以当前道为中心向两侧延伸预设道数的地震道，优选 地分布在当前道两侧的临近道的数量相等。当然，在本发明的其他实施例中，分布在当前到 两侧的临近道的数量也可以不等，本发明不限于此。
[0049] 当前道分析点是已知的，那么过当前道分析点以预设扫描倾角构建一条直线，而 这条直线与各个临近道的交点便是各个临近道的临近道分析点。对于不同的预设扫描倾 角，各个临近道的临近道分析点是不同的。所以在步骤S102中，根据当前道分析点和各个 临近道分析点，便可以分别获取当前道的地震波形数据和对应于各个预设扫描倾角的临近 道的地震波形数据。
[0050] 本实施例中，对于当前道，获取以当前道分析点为中心的预设深度窗口的地震波 形数据，并将获取的地震波形数据作为当前道的地震波形数据。同样，对于各个临近道，获 取以相应临近道分析点为中心的预设深度窗口的地震波形数据，并将获取的地震波形数据 作为相应临近道的地震波形数据。
[0051] 需要说明的是，在本发明的其他实施例中，当前道和各个临近道的地震波形数据 还可以通过其他合理方式获得，例如对于当前道，将以当前道分析点为起点的预设深度窗 口的地震波形数据作为当前道的地震波形数据，而各个临近道的地震波形数据也可以通过 同样方式获得，但本发明不限于此。
[0052] 本实施例所提供的地层检测方法通过利用预设窗口长度的深度窗口来其他区相 应的地震波形数据，这样相较于现有的地层倾角检测方法，能够明显减少同相轴稀疏性的 影响，从而得到更为稳定、精度更高的倾角剖面。
[0053] 获得当前道和各个临近道的地震波形数据后，本实施例在相似能量谱计算步骤 中，根据当前道的地震波形数据和对应于各个预设扫描倾角的临近道的地震波形数据，来 确定各个预设扫描倾角的相似能量谱。
[0054] 如图1所示，具体地，本实施例在步骤S103中，根据当前道的地震波形数据和对应 于各个预设扫描倾角的临近道的地震波形数据，计算对应于各个预设扫描倾角的各个临近 道的地震波形数据的相关值。本实施例中，根据如下公式计算各个临近道的地震波形数据 的相关值：
[0056] 其中，Cini表示对应于倾角m的第i个临近道的地震波形数据的相关值，Wd表示预 设深度窗口