一种用于求取裂缝各向异性梯度的方法
【技术领域】
[0001] 本发明涉及地球物理技术领域，具体地说，涉及一种求取地层裂缝的各向异性梯 度的方法。
【背景技术】
[0002] 在裂缝性油气藏研究中，需要在裂缝预测过程中获得定量的裂缝参数。裂缝参数 常选用可W很好的反应裂缝发育强度和方向的各向异性梯度。由于各向异性参数的数量级 远远小于地震数据的数量级，因此，多数研究选择将Ruger Η??介质纵波近似反射系数方程 简化，将叠前方位数据作为该方程的输入数据来对该方程反演来得到各向异性梯度。将叠 前方位数据作为输入数据虽然保留着野外采集时地震波的振幅随炮检距变化的规律，但地 震波数据中有效反射信号信噪比较低，反演结果的准确性和可信度受到叠前资料的直接影 响。
【发明内容】
[0003] 为解决上述问题，本发明提供了一种W叠后地震方位数据为基础获取裂缝各向异 性梯度的方法，用W避免叠前地震资料对各向异性梯度结果的直接影响。
[0004] 根据本发明的一个实施例，提供了一种用于求取裂缝各向异性梯度的方法，包括 W下步骤：
[0005] 步骤一、基于采集的地震数据获取叠前方位道集的角度信息；
[0006] 步骤二、基于所述角度信息对所述叠前方位道集进行叠加处理W获取对应的叠后 方位地震数据体；
[0007] 步骤H、对所述叠后方位地震数据体进行稀疏脉冲反演W获取对应的反射系数 体；
[0008] 步骤四、基于简化的Ruger Η??介质纵波近似反射系数方程对所述反射系数体和 所述角度信息进行处理W获取各向异性参数和相关变量；
[0009] 步骤五、基于所述各向异性参数和所述相关变量计算各向异性梯度Δ Γ。
[0010] 根据本发明的一个实施例，所述各向异性参数包括HTI介质上下两层纵波变异系
数差值Δ δ W，所述相关变量为 其中，_为Η??介质上下两层纵波速度平均值， 化 房为ΗΤΙ介质上下两层横波速度平均值，Δ Υ为ΗΤΙ介质上下两层横波各向异性差值。
[0011] 根据本发明的一个实施例，所述角度信息包括叠前方位道集的方位角信息和入射 角信息。
[0012] 根据本发明的一个实施例，在步骤二中，通过将方位角对应的叠前方位道集按照 相同的入射角方向进行叠加 W获取该方位角对应的叠后方位地震数据体。
[0013] 根据本发明的一个实施例，在步骤Η中，对所述叠后方位地震数据体进行稀疏脉 冲反演w获取对应的反射系数体的过程包括：
[0014] 基于最大似然反權积建立所述叠后方位地震数据体的最小目标函数并通过迭代 W求取反射系数；
[0015] 基于所述反射系数和初始阻抗模型，采用递推算法得到初始波阻抗模型；
[0016] 基于测井约束条件对所述初始波阻抗模型和所述反射系数进行调整W获取对应 的反射系数体。
[0017] 根据本发明的一个实施例，所述简化的Ruger Η??介质纵波近似反射系数方程 为：
[0018]
[0019] 其中，Θ为入射角，^为方位角，巧为HTI介质上下两层纵波速度平均值，梦为 HTI介质上下两层横波速度平均值，户为HTI介质上下两层密度平均值，Δ α为HTI介质上 下两层纵波速度差值，A目为ΗΤΙ介质上下两层横波速度差值，Δ Ρ为ΗΤΙ介质上下两层 密度差值，Δ δ W为Η??介质上下两层纵波变异系数差值，Δ Y为HTI介质上下两层横波 各向异性差值。
[0020] 根据本发明的一个实施例，在步骤四中，进一步包括：
[0021] 基于所述入射角信息获取叠前方位道集的最大入射角Θ 2和最小入射角Θ 1，并设 定累加的间隔，将所述简化的Ruger ΗΤΙ介质纵波近似反射系数方程按入射角方向累加 W 得到：
[0022]
[0026] 基于按入射角方向累加后的简化的Ruger HTI介质纵波近似反射系数方程，任选 Η个不同方位角及对应的反射系数体进行处理W得到所述HTI介质上下两层纵波变异系 数差值Δ δ W和所述相关变量
[0027] 根据本发明的一个实施例，在步骤五中，基于Ruger各向异性梯度计算公式对所 述HTI介质上下两层纵波变异系数差值Δ δ W和所述相关变量
进行处理W得到 各向异性梯度Δ Γ，所述Ruger各向异性梯度计算公式为：
[0028]
[0029] 本发明带来了 W下有益效果：
[0030] 本发明W叠后方位数据为基础，避免了叠前地震资料对反演计算结果的直接影 响，使得反演方法稳定性高、信噪比高。同时，本发明计算精度不受输入数据体的个数影响， 在实际应用过程中易于实现。本发明直接选取地震数据品质高的部分进行叠加作为反演输 入数据，避免了叠前角度域道集存在的能量分布不均的问题。
[0031] 本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变 得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利 要求书W及附图中所特别指出的结构来实现和获得。
【附图说明】
[0032] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要的附图做简单的介绍：
[0033] 图1是根据本发明的一个实施例的方法流程图；
[0034] 图2是根据本发明的一个实施例的正演得到的10°、45°、90°方位的叠前数据；
[0035] 图3是根据本发明的一个实施例的反演得到的10°、45°、90°方位对应的反射 系数拟及
[0036] 图4是根据本发明的一个实施例的计算得到的各向异性梯度Δ Γ。
【具体实施方式】
[0037] W下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用 技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据W实施。需要说明 的是，只要不构成冲突，本发明中的各个实施例W及各实施例中的各个特征可W相互结合， 所形成的技术方案均在本发明的保护范围之内。
[0038] 另外，在附图的流程图示出的步骤可W在诸如一组计算机可执行指令的计算机系 统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可不同于此处 的顺序执行所示出或描述的步骤。
[0039] 如图1所示为根据本发明的一个实施例的方法流程图。W下参考图1来对本发明 所述的方法进行详细说明。
[0040] 在步骤S110中，基于采集的地震数据获取叠前方位道集的角度信息。
[0041] 在该步骤中，在进行叠加处理之前，基于采集到的地震数据获取地震数据的叠前 方位道集，并同时获取对应叠前方位道集的方位角信息和叠前方位道集的入射角信息。
[0042] 在步骤S120中，对叠前方位道集进行叠加处理W获取对应的叠后方位地震数据 体。在该步骤中，针对每个叠前方位道集按相同的入射角方向进行叠加处理W获取对应的 叠后方位地震数据体。
[0043] 在步骤S130中，对叠后方位地震数据体进行稀疏脉冲反演W获取对应的反射系 数体。现有的稀疏脉冲反演过程可分为W下几个步骤。
[0044] 首先，采用最大似然反權积进行反射系数的反演。最大似然反權积对地层的假设 认为；地层的反射系数是由较大的反射界面的反射和具有高斯背景的小反射叠加组合而 成，导出一个最小目标函数：
[0045]
[004引其中，R2和妒分别为反射系数的均方值和噪音的均方值，r似和η似表示第K个 采样点的反射系数和噪音，Μ表示反射层数，Ν表示采样总数，λ 1表示给定反射系数的似然 值。通过对该最小目标函数进行多次迭代来求取反射系数。
[0047] 然后，根据最大似然反權积计算得到的反射系数序列，结合初始阻抗模型，采用递 推算法，反演得到初始的波阻抗模型：
[0048]
(2)
[004引其中，Z(i)为第i层的波阻抗值，R(i)为第i层的反射系数。
[0050] 最后，结合测井约束条件进行波阻抗反演。约束稀疏脉冲反演对每一道依据目标 函数对计算出的初始波阻抗进行调整，包括对反射系数的调整。目标优化函数为：
[0051] F = Lp (r) + λ zLq (s-d)+a iLi Δ Z (3)
[0052] 其中，r为反射系数序列，ΔΖ为阻抗趋势的差序列，d为地震道序列，s为合成地 震道序列，λ2为残差权重因子，a为趋势权重因子，p、q为L模因子。具体的，右式第一项 反映了反射系数的绝对值和，第二项反映了合成声波记录与原始地震数据的差值，第Η项 为趋势约束项。
[0053] 通过稀疏脉冲反演计算得到分别对应不同方位角的反射系数体。
[0054] 在步骤S140中，基于简化的Ruger Η??介质纵波近似反射系数方程对所述反射系 数体和所述角度信息进行处理W获取各向异性参数和相关变量。其中，各向异性参数包括 Η??介质上下两层纵波变异系数差值Δ δ W，相