[发明专利]一种基于HTI介质的各向异性单参数反演方法

说明书

技术领域

本发明属于地球物理反演领域，具体涉及一种基于HTI介质的各向异性单参数反演方法，适用于HTI(横向各向异性介质)。

背景技术

随着地质勘探的不断发展，对裂缝性油气藏的研究也越来越广泛，在裂缝预测中对裂缝参数的定量预测主要通过各向异性参数得到。通过叠前方位道集获取各向异性参数目前主要分为三种方法：

1.通过分析NMO速度的方位变化计算裂缝参数，当速度分析精度比较高时，分析NMO速度的方位变化就能够做方位各向异性分析，但在实际问题中，很难确定速度的方位各向异性；

2.通过叠前方位道集中提取剩余方位时差计算，叠前道集中剩余时差比较容易拾取，该方法易于实现，但是由于拾取剩余时差时存在较多误差，无法区分产生剩余方位时差的原因是否是各向异性引起的，导致该方法计算结果准确性较低；

3.利用纵波反射振幅随偏移距(入射角)和方位变化的信息，以Ruger(1998)方程为基础，通过推导近似等转换来化简Ruger方程，计算出与各向异性参数相关的综合参数，这些方法大多相对复杂，假设条件多并且计算效率不高。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种基于HTI介质的各向异性单参数反演方法，基于简化的Ruger方程，通过叠后反演的方法解决叠前各向异性参数反演问题。

本发明是通过以下技术方案实现的：

一种基于HTI介质的各向异性单参数反演方法，包括：

(1)选取任意三个不同方位角对应的叠后地震数据体，通过稀疏脉冲的反演方法计算得到对应该三个不同方位角的反射系数

(2)分析叠前方位道集获得入射角信息，所述入射角信息包括最大入射角θ₂、最小入射角θ₁和间隔，然后通过公式计算出B的值；

(3)利用方位角信息反射系数和参数B，计算得到各向异性参数Δε^(V)，完成反演。

所述方法进一步包括：

(4)将步骤(3)反演得到的Δε^(V)与井旁道或模型数据的Δε^(V)值进行对比，分析反演得到的Δε^(V)值与井旁道或模型数据的Δε^(V)值之间的误差。

所述步骤(2)中，0≤θ₁＜θ₂≤30。

所述步骤(3)是利用下面的公式计算得到各向异性参数Δε^(V)：

与现有技术相比，本发明的有益效果是：

(1)目前叠前反演受叠前道集的影响，稳定性不高，本方法将叠前道集在入射角域进行叠加，将叠加后的具有方位信息的叠后方位体作为反演的输入，使反演方法具有叠后地震反演稳定性高、信噪比高的优点；

(2)目前叠前各向异性参数反演的理论要求输入数据体的个数尽量多，反演方法则更加稳定，而实际应用过程中由于计算效率及数据品质的限制难以得到多个输入数据体，本反演方法从理论上只需要3个叠前方位体数据，在实际应用过程中易于实现；

(3)目前由于观测系统的影响，绝大多数叠前角度域道集存在，近、远道能量弱，中入射角道集能量强的缺陷，本反演方法从理论上实现了叠前角度域道集的数据筛选功能，可直接选取地震数据品质高的部分进行叠加作为反演输入数据；

(4)反演参数少，计算方法简洁。

附图说明

图1本方法的步骤框图

图2理论模型图

图3正演得到的0°方位的叠后数据体

图4正演得到的30°方位的叠后数据体

图5正演得到的60°方位的叠后数据体

图6反演得到的0°方位叠后数据体对应的反射系数

图7反演得到的30°方位叠后数据体对应的反射系数

图8反演得到的60°方位叠后数据体对应的反射系数

图9计算得到的各向异性参数Δε^(V)值

图10模型中各向异性参数Δε^(V)值。

具体实施方式

下面结合附图对本发明作进一步详细描述：

直接进行叠前各向异性参数反演时，Ruger方程内同时含有纵波速度、横波速度、密度和Δδ^(V)，Δε^(V)和Δγ六个未知变量，这就使方程求解时引入多种近似关系且解法复杂，稳定性较低。