一种调整用于叠前偏移成像的地震数据能量一致性的方法
【技术领域】
[0001] 本发明涉及地震勘探资料处理技术，具体地，涉及一种调整用于叠前偏移成像的 地震数据能量一致性的方法。
【背景技术】
[0002] 目前，随着油气勘探的深入和数字计算技术的飞速发展，叠前偏移已取代叠后偏 移，成为解决复杂地质条件精细成像的常规处理技术。特别是基于射线理论的Kichhoff?积 分法叠前偏移方法，在目前工业生产中应用最为广泛。
[0003] Kichhoff积分法偏移的原理是沿椭球面对地层振幅求和，相邻面元能量不均通常 会产生画弧现象，影响地震资料的成像精度。叠前偏移对道集数据的能量一致性要求相对 较高，为了消除地震波在传播过程中球面扩散、吸收衰减、震源能量差异、检波器耦合差异 以及不同时间采集的区块间的能量差异对反射波振幅的影响，现有的处理流程中通常采用 几何扩散补偿、地表一致性振幅补偿和剩余振幅补偿来消除能量差异，尽量提高振幅保真 度，保持地震波振幅能量的相对强弱关系。
[0004] 为更好的改善能量不均对叠前时间偏移成像的影响，目前在工业化生产处理中， 正逐步得到重视。吴清岭等提出基于覆盖次数的叠前振幅归一化处理，改善振幅一致性处 理对连片处理的影响；张学涛等为了解决在地震资料连片处理中，能量不均严重影响叠前 时间偏移成像质量的问题，采用了 CMP域逐点自动能量调整法；张伟等在假设共偏移距道 集覆盖次数均匀情况下，分析了能量一致性对偏移成像影响，研究认为如果处理过程中振 幅补偿不均就会造成分辨率横向变化。综合采用几何扩散补偿、地表一致性振幅补偿和剩 余振幅补偿来消除能量差异，利用监控手段保证补偿效果，可以有效提高偏移成像质量。
【发明内容】
[0005] 目前，地震资料的处理中通常仍是采用补偿前后的道集、纯波叠加剖面以及均方 根振幅属性的监控，判断能量补偿的效果。
[0006] 但是，在实现本发明的过程中，本发明人发现：为有效提高偏移成像质量，目前的 这些监控方法对于叠后偏移是能满足的，但对于叠前偏移仍是不足的。原因在于叠前偏移 的工业化算法设计，其在实现过程中没有考虑观测系统设计和切除等原因造成的成像点间 的覆盖次数差异以及炮检距的分布的均匀性带来的能量差异。为利用叠前偏移方法更好的 提高成像质量，一方面要充分利用质量监控做好偏移前振幅补偿；另一方面还要从叠前偏 移的工业化算法着手，弥补算法的不足，在偏移前做好道集能量调整，以解决考虑观测系统 设计和切除等原因造成的成像点间的覆盖次数差异以及炮检距的分布的均匀性带来的能 量差异。
[0007] 为此，非常需要一种调整用于叠前偏移成像的地震数据能量一致性的方法。
[0008] 为了实现上述目的，本发明实施例提供一种调整用于叠前偏移成像的地震数据能 量一致性的方法，包括：
[0009] 步骤I，采集地震勘探数据；
[0010] 步骤2,对所述地震勘探数据进行预处理，并分选成共中心点道集；
[0011] 步骤3,将所述共中心点道集分组；
[0012] 步骤4,针对每一组都进行如下处理：
[0013] 步骤41，按照如下公式计算考虑覆盖次数的有效叠加地震数据：
[0014]
[0015] 步骤42,按照如下公式计算不考虑覆盖次数的有效叠加地震数据：
[0016]
[0017] 步骤43,按照如下公式计算能量调整因子：
[0018] opz (t) = szl(t)/sz2(t)
[0019] 步骤44,按照如下公式将所述能量调整因子应用于该组的所有地震数据：
[0020] A-；(〇 = .Yy (/)〇/;(/)
[0021] 其中，z表示该组在所述共中心点道集中的编号；
[0022] szl (t)表示考虑覆盖次数的有效叠加地震数据；
[0023] t表示时间；
[0024] N表示覆盖次数；
[0025] j表示该组的地震数据的编号；
[0026] X]⑴表示该组的地震数据；
[0027] sz2 (t)表示不考虑覆盖次数的有效叠加地震数据；
[0028] opz (t)表示能量调整因子；
[0029] 表示应用能量因子调整后的地震数据；
[0030] 步骤5,得到能量调整后的、用于叠前偏移成像的地震数据。
[0031] 借助于上述技术方案，本发明利用考虑覆盖次数的有效叠加除以不考虑覆盖次数 的有效叠加计算叠前偏移前能量调整因子，不仅考虑了覆盖次数的影响，更考虑到了实际 叠加数据的有效值，是纯数据驱动，能量调整更稳定，容易理解和操作。应用本发明进行叠 前偏移前能量调整，更加符合叠前偏移的工业化算法，可获得更好的地震勘探成像质量。
【附图说明】
[0032] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述 中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些 实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附 图获得其他的附图。
[0033] 图1是本发明示例性方法的具体流程图；
[0034] 图2是本发明实施例一的具体流程图；
[0035] 图3是未应用本实施例进行的叠前时间偏移成像结果；
[0036] 图4是应用本实施例进行的叠前时间偏移成像结果。
【具体实施方式】
[0037] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0038] 本领域技术技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、 方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软 件（包括固件、驻留软件、微代码等），或者硬件和软件结合的形式。
[0039] 根据本发明的实施方式，提出了一种调整用于叠前偏移成像的地震数据能量一致 性的方法。
[0040] 在本文中，需要理解的是，附图中的任何元素数量均用于示例而非限制，以及任何 命名都仅用于区分，而不具有任何限制含义。
[0041] 示例件方法
[0042] 本示例性方法提供一种调整用于叠前偏移成像的地震数据能量一致性的方法，如 图1所示，该方法包括：
[0043] 步骤S11，采集地震勘探数据。
[0044] 步骤S12,对地震勘探数据进行预处理，并分选成共中心点道集。
[0045] 可选地，该步骤可以对地震勘探数据进行包括静校正，去噪，球面扩散补偿、地层 吸收衰减、地表一致性等能量补偿，以及子波一致性校正等预处理。
[0046] 步骤S13,将共中心点道集分组。
[0047] 可选地，该步骤可以选择按炮检距分组，也可以选择炮检距和方位角分组，也可以 按炮检距向量片技术（Offset Vector Tile，简称0VT)进行分组。
[0048] 分组后，每组在该共中心点道集中都有唯一的编号。
[0049] 步骤S14,针对每一组都进行如下处理：
[0050] 步骤S141，按照如下公式计算考虑覆盖次数的有效叠加地震数据：
[0051]
(公式 1)
[0052] 步骤S142,按照如下公式计算不考虑覆盖次数的有效叠加地震数据：
[0053]
(公式 2)
[0054] 步骤S143,按照如下公式计算能量调整因子：
[0055] opz(t) = szl(t)/sz2(t) (公式 3)
[0056] 步骤S144,按照如下公式将能量调整因子应用于该组的所有地震数据：
[0057] X11 (t) = X1(I