一种压制地震资料深层多次波的方法技术

本发明专利技术提供了一种压制地震资料深层多次波的方法，属于地震数据处理技术领域。考虑到层间多次波在频率和速度方面与一次有效波重叠，难以识别，本发明专利技术根据地层反射系数确定若干多次波算子，通过多次波算子对地震资料中的多次波进行压制，通过比较压制后得到的一次有效波的能量，确定出对应能量最小的一次有效波，进而确定出一次有效波对应的压制后地震资料。采用本发明专利技术，能够提高地震资料的保真比，满足油气勘探需求。足油气勘探需求。足油气勘探需求。

【技术实现步骤摘要】
一种压制地震资料深层多次波的方法


[0001]本专利技术涉及一种压制地震资料深层多次波的方法，属于地震数据处理
，尤其涉及地震资料处理时保真提高信噪比处理阶段中压制多次波的方法。

技术介绍

[0002]复杂地质构造的深层勘探开发展示了深层油气藏勘探的良好前景，然而在深层地质勘探过程中，由于上覆地层存在强反射界面，造成深层多次波发育，严重影响地震资料的信噪比和保真性，制约着油气勘探的进程。
[0003]现有技术中，为解决地震资料处理中深层多次波的问题，通常利用多次波与一次有效波在频率、速度等方面的差异进行压制，例如，公开号为CN112925023A的中国专利文件，公开了一种全波场反演地震数据多次波压制方法，该方法基于τ
‑
p变换，利用双平面波域对地震数据的压缩特性，通过多次迭代反演进行多次波压制和一次波估计。然而由于层间多次波在频率和速度方面与一次有效波重叠，难以识别，导致对深层多次波的压制效果有限，难以满足当前精细勘探的需求。

技术实现思路

[0004]本专利技术的目的在于提供一种压制地震资料深层多次波的方法，用于解决对地震资料处理中对深层多次波压制时难以对层间多次波进行识别和压制的问题。
[0005]为了实现上述目的，本专利技术提供了一种压制地震资料深层多次波的方法，包括如下步骤：S1、对待压制地震资料进行τ
‑
p变换，得到τ
‑
p域地震资料；
[0006]S2、通过预先得到的若干多次波算子，分别对τ
‑r/>p域地震资料中的层间多次波进行压制，得到对应各多次波算子的一次有效波地震资料；
[0007]S3、计算各一次有效波地震资料的能量并进行比较，确定出与能量最小的一次有效波地震资料对应的多次波算子作为最优多次波算子；
[0008]S4、将通过最优多次波算子对τ
‑
p域地震资料中的层间多次波进行压制得到的一次有效波地震资料进行反τ
‑
p变换，得到多次波压制后地震资料。
[0009]进一步地，在上述方法中，步骤S2中所述的多次波算子根据地层反射系数进行确定，将地层反射系数与待压制地震资料做不同预测步长的褶积运算得到不同多次波算子。
[0010]进一步地，在上述方法中，通过如下公式得到一次有效波地震资料：
[0011]P
‑
A(ω)P2+A2(ω)P3+
…
＝P0[0012]式中，P为τ
‑
p域地震资料的总波场响应，A(ω)为多次波算子，P0为一次有效波地震资料的一次有效波响应。
[0013]进一步地，在上述方法中，步骤S3中，通过对如下公式进行求解得到一次有效波地震资料能量：
[0014][0015]式中，为第i炮的能量，i为炮数，表示第i炮总波场响应，t为地震波在介质中的传播时间，x
r
表示检波点位置，表示第i个炮点的位置，表示滤波算子，*表示褶积，n为迭代次数，为多次波响应。
[0016]进一步地，在上述方法中，还通过沿着P轴改变预测步长的方法对τ
‑
p域地震资料进行多次波压制。
[0017]进一步地，在上述方法中，通过求解如下公式得到所述预测步长Δt(p)：
[0018][0019]式中，为n阶多次波旅行时，为n
‑
1阶多次波旅行时，z为介质深度，Δt(0)为地面接收点到地下反射点的垂直旅行时，p为地震资料的水平波慢度，v为多次波在介质中的传播速度。
附图说明
[0020]图1为本专利技术方法实施例中压制地震资料深层多次波的方法的流程框图；
[0021]图2为本专利技术方法实施例中不同的多次波算子A(ω进行多次波压制时的对比图；
[0022]图3
‑
a为本专利技术方法实施例中压制前的地震资料；
[0023]图3
‑
b为本专利技术方法实施例中的多次波模型；
[0024]图3
‑
c为本专利技术方法实施例中采用多次波模型进行多次波压制后的地震资料；
[0025]图4为本专利技术方法实施例中采用常规压制方法和本专利技术方法后道集的效果对比图；
[0026]图5为采用常规方法与采用本专利技术进行多次波压制后的地震资料的偏移剖面对比图；
[0027]图6为本专利技术方法实施例中采用常规压制方法和本专利技术方法进行多次波压制后的地震记录对比示意图。
具体实施方式
[0028]为了使本专利技术的目的、技术方案及优点更加清楚明了，以下结合附图及实施例，对本专利技术进行进一步详细说明。
[0029]方法实施例：
[0030]本专利技术的一种压制地震资料深层多次波的方法，如图1所示，包括如下步骤：
[0031]S1、对处于时空域(T
‑
x域)的地震资料进行动校正，并在动校正后进行τ
‑
p变换，从而得到τ
‑
p域地震资料。
[0032]S2、τ
‑
p域地震资料在每一个P道上的长周期全程多次波具有较好的周期性，因此，通过沿着P轴改变预测步长的方法，对τ
‑
p域地震资料中的长周期全程多次波进行压制。
[0033]预测步长Δt(p)的推导公式如下所示：
[0034][0035]式中，为n阶多次波旅行时，为n
‑
1阶多次波旅行时，z为介质深度，
Δt(0)为地面接收点到地下反射点的垂直旅行时，p为地震资料的水平波慢度，p＝sinθ/v，v为多次波在介质中的传播速度，θ为波传播方向与竖直方向的夹角。
[0036]作为其他实施方式，预测步长还可根据经验确定，预测步长Δt(p)取值为4ms。
[0037]在第1个P道上预测步长为Δt(p)，在第2个P道上预测步长为2Δt(p)，第3个P道上预测步长为3Δt(p)，
…
，第k个P道上预测步长为kΔt(p)。在每个P道上根据对应的预测步长来预测长周期全程多次波，进而对τ
‑
p域地震资料中的长周期全程多次波进行压制。
[0038]S3、采用数据驱动的自适应短周期全程多次波及层间多次波的压制方式，通过预先得到的若干多次波算子，以及如下公式，分别对步骤S2中对短周期全程多次波进行压制后的地震资料进行层间多次波压制，得到与各多次波算子对应的一次有效波地震资料。
[0039]P
‑
A(ω)P2+A2(ω)P3+
…
＝P0[0040]式中，P为τ
‑
p域地震资料的总波场响应，A(ω)为多次波算子，P0为一次有效波地震资料的一次有效波响应。
[0041]多次波算子通过如下方法得到：根据测井资料得到地层反射系数，将地层反射系数与地震资料做不同预测步长的褶积运算，从而得到不同多次波算子A(ω)。多次波算子的数值过小会损坏地震资料中的一次有效波，数值过大会导致本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种压制地震资料深层多次波的方法，其特征在于，包括如下步骤：S1、对待压制地震资料进行τ
‑
p变换，得到τ
‑
p域地震资料；S2、通过预先得到的若干多次波算子，分别对τ
‑
p域地震资料中的层间多次波进行压制，得到对应各多次波算子的一次有效波地震资料；S3、计算各一次有效波地震资料的能量并进行比较，确定出与能量最小的一次有效波地震资料对应的多次波算子作为最优多次波算子；S4、将通过最优多次波算子对τ
‑
p域地震资料中的层间多次波进行压制得到的一次有效波地震资料进行反τ
‑
p变换，得到多次波压制后地震资料。2.根据权利要求1所述的压制地震资料深层多次波的方法，其特征在于，步骤S2中所述的多次波算子根据地层反射系数进行确定，将地层反射系数与待压制地震资料做不同预测步长的褶积运算得到不同多次波算子。3.根据权利要求2所述的压制地震资料深层多次波的方法，其特征在于，通过如下公式得到一次有效波地震资料：P
‑
A(ω)P2+A2(ω)...

【专利技术属性】
技术研发人员：杨兴圣，张高成，赵伟，李俨，尹安，
申请(专利权)人：中国石油化工股份有限公司河南油田分公司勘探开发研究院，
类型：发明
国别省市：