[发明专利]一种各向异性介质地震反射特征方程建立方法

说明书

本发明涉及油气勘探开发领域，具体公开了一种各向异性介质地震反射特征方程建立方法，该方法将现有的采用各向异性介质反射系数近似公式，改造为一个能够进行OA介质各向异性反演的三项PP波近似公式，改进公式在保证原公式精度的基础上，减少了待求参数的个数，能够更好地适用于页岩气储层的岩性预测和流体解释。

技术领域

本发明涉及油气勘探开发领域，具体包含一种各向异性介质地震反射特征方程建立方法。

背景技术

伴随着非常规油气储层研究的发展，尤其是页岩储层研究的深入，适用于垂直裂缝发育的页岩地层等效正交各向异性介质(OA)成为地球物理学家们关注的焦点。当研究发育垂向裂缝的非常规页岩气储层时，仅考虑页岩成层性引起的VTI特征影响是不合理的。然而，Ruger提出的近似公式仅能计算两个平面内的反射系数近似值，无法全方位展示页岩储层的地震反射特征，在实际发育垂直裂缝的页岩储层中应用时具有一定的局限。受δ,ε,γ等多参数的影响，现有反射特征方程开展各向异性反演是一件困难的事，因此需要建立新的有利于反演的方程指导实际勘探，为油气地震勘探提供更加可靠的技术支持。

发明内容

本发明要解决的技术问题是：克服现有关于基于计算OA介质反射系数近似公式的反演方法待求参数多且敏感性差异大、正则化方法多样、反演结果多解性强等缺陷，提出一种新的计算方法，一个能够进行OA介质各向异性反演的三项PP波近似公式，改进公式在保证原公式精度的基础上，减少了待求参数的个数，能够更好地适用于页岩气储层的岩性预测和流体解释。

为了解决上述技术问题，本发明所采用的技术方案是：一种各向异性介质地震反射特征方程建立方法，设采用各向异性介质反射系数近似公式如下：

R_PP(θ,φ)＝P+(G_iso+Γ_xcos²φ+Γ_ysin²φ)sin²θ (1)，

公式(1)中：

Γ_x和Γ_y为不同对称轴方向的各向异性梯度，Z＝ρα，代表纵波阻抗，θ为入射角，φ为方位角，α、β分别为各向同性背景部分的纵、横波速度，ρ为密度，Δ代表不同的界面之间的差异，^-代表两层的平均值，是各向同性背景部分的纵波反射系数，弱各向异性参数ε、δ及γ的下标x、y及z代表不同对称轴平面内的各向异性参数定义；

所述各向异性介质地震反射特征方程建立方法包括以下步骤：

第一步：对公式(1)中的G_iso乘以sin²φ+cos²φ，即得：

公式(2)式整理得到：

假设即得：

令即得：

公式(5)整理可得：

公式(6)即为：

其中，

A＝ρα，

第二步：公式(7)写为：

公式(8)改写成矩阵形式：

R＝Αm (9)，