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

摘要

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

主权项

1.一种各向异性介质地震反射特征方程建立方法，采用各向异性介质反射系数近似公式如下：R_PP(θ,φ)＝P+(G_iso+Γ_xcos²φ+Γ_ysin²φ)sin²θ   (1)，公式(1)中：Γ_x和Γ_y为不同对称轴方向的各向异性梯度，Z＝ρα，代表纵波阻抗，θ为入射角，φ为方位角，α、β分别为各向同性背景部分的纵、横波速度，ρ为密度，Δ代表不同的界面之间的差异，‑代表两层的平均值，是各向同性背景部分的纵波反射系数，弱各向异性参数δ及γ的下标x、y及z代表不同对称轴平面内的各向异性参数定义；其特征在于，所述各向异性介质地震反射特征方程建立方法包括以下步骤：第一步：对公式(1)中的G_iso乘以sin²φ+cos²φ，即得：公式(2)式整理得到：假设即得：令即得：公式(5)整理可得：公式(6)即为：其中，Α＝ρα，第二步：公式(7)写为：公式(8)中，m₁,m₂,m₃为模型参数，定义定义a₁(θ,φ),a₂(θ,φ),a₃(θ,φ)为模型参数权重系数，定义公式(8)改写成矩阵形式：R＝Α(θ,φ)m   (9)，地震道集褶积模型用矩阵表示为：S＝WR   (10)，由公式(9)和(10)得到：S＝Gm   (11)，公式(11)中，S为地震数据矩阵，W为子波矩阵，G为子波与系数矩阵的联合矩阵，在贝叶斯反演框架下构建目标函数时，先假设待反演参数先验服从柯西分布，似然函数服从高斯分布，然后采用平滑背景模型参数构建参数去相关协方差矩阵，实现待反演参数去相关并提高横向连续性，同时将平滑背景模型约束加入目标函数，最终得到目标泛函为：公式(12)中，M为模型参数个数，其中：分别为三个待反演参数约束系数，和分别为噪声分布和模型参数的方差，G'为去相关后子波系数矩阵，m'为待反演去相关后模型参数反射系数的转置矩阵，表示积分公式