[发明专利]一种地震数据处理方法及装置

说明书

本申请提供一种地震数据处理方法及装置。所述方法包括：基于获取的地震角道集数据，按照入射角的大小将所述地震角道集数据划分为不同角度范围道集数据；根据建立的反射系数函数，将所述不同角度范围道集数据分别进行各向异性参数反演，获得不同角度范围对应的纵波速度，所述反射系数函数包括：不同入射角、纵波速度、地层密度与地震波的反射之间的关系；根据所述反射系数函数、所述不同角度范围对应的纵波速度，获得所述地震角道集数据对应的各向异性参数。利用本申请中各实施例，实现了页岩各向异性参数反演，获得了准确且高分辨率的各向异性参数，为后续地震资料的处理提供了准确的数据基础。

技术领域

本申请属于地球物理勘探技术领域，尤其涉及一种地震数据处理方法及装置。

背景技术

随着对地下介质认识的不断深入，人们意识到简单的将地球介质假设为各向同性介质不能满足现有的勘探开发需求。地震反演技术也开始在各向异性介质中发展起来，但是主要的工作集中在对各向异性介质的裂缝预测中。

在页岩储层发育的地区，通常呈现较强的横向各向异性性质，对于VTI(VerticalTransverse Isotropy，横向各向同性介质)介质各向异性参数的定量反演技术还存在着一定的难度。现有技术中，对于VTI介质的各向异性反演方法研究还不够充分，反演方法存在不稳定性、横波数据信噪比低、纵横波匹配不准确、待求参数敏感性差异大、正则化方法多样、反演结果多解性强等，导致实际叠前道集的弹性参数和各向异性参数反演困难。因此，业内亟需一种能够实现页岩各向异性参数反演，进一步提高地震数据处理准确性的实施方案。

发明内容

本申请目的在于提供一种地震数据处理方法及装置，实现了页岩各向异性参数反演，获得了准确且高分辨率的各向异性参数，为后续地震资料的处理提供了准确的数据基础。

一方面本申请提供了一种地震数据处理方法，包括：

基于获取的地震角道集数据，按照入射角的大小将所述地震角道集数据划分为不同角度范围道集数据；

根据建立的反射系数函数，将所述不同角度范围道集数据分别进行各向异性参数反演，获得不同角度范围对应的纵波速度，所述反射系数函数包括：不同入射角、纵波速度、地层密度与地震波的反射之间的关系；

根据所述反射系数函数、所述不同角度范围对应的纵波速度，获得所述地震角道集数据对应的各向异性参数。

进一步地，所述方法的另一个实施例中，所述反射系数函数包括：

上式中，表示反射系数，θ表示入射角与透射角的平均角度，V_P(θ)表示入射角为θ时的纵波速度，表示入射角为θ时的纵波速度的平均值，V_S0表示垂直方向横波速度，表示垂直方向横波速度的平均值，V_P0表示垂直方向纵波速度，表示垂直方向纵波速度的平均值，ρ表示地层密度，ρ₁表示上层介质的地层密度，ρ₂表示下层介质的地层密度，表示地层密度平均值，δ、ε表示所述各向异性参数，δ₁、ε₁表示上层介质的各向异性参数，δ₂、ε₂表示下层介质的各向异性参数，V_P1(θ)表示入射角为θ时上层介质的纵波速度，V_P2(θ)表示入射角为θ时下层介质的纵波速度，V_P01表示上层介质的垂直纵波速度，V_P02表示下层介质的垂直纵波速度。

进一步地，所述方法的另一个实施例中，所述根据建立的反射系数函数，将所述不同角度范围道集数据分别进行各向异性参数反演，包括：

根据所述反射系数函数，构建所述不同角度范围道集数据对应的纵波速度函数；