[发明专利]一种砂岩储层识别方法

权利要求书

1.一种砂岩储层识别方法，其特征在于，包括以下步骤：

步骤S1、提取敏感识别因子，所述敏感识别因子为纵波速度和横波速度的乘积与密度的n次方的比值，其中，n＝1、2或4；

步骤S2、利用所述敏感识别因子识别砂岩储层。

2.根据权利要求1所述的砂岩储层识别方法，其特征在于，所述敏感识别因子为：SF＝v_pv_s/ρ⁴，其中，SF为所述敏感识别因子，v_p为所述纵波速度，v_s为所述横波速度，ρ为所述密度。

3.根据权利要求1所述的砂岩储层识别方法，其特征在于，提取敏感识别因子，具体为：

获取包含所述敏感识别因子的AVO近似式，以所述AVO近似式为核函数直接反演得到所述敏感识别因子。

4.根据权利要求3所述的砂岩储层识别方法，其特征在于，获取包含所述敏感识别因子的AVO近似式，具体为：

获取Aki-Richards近似式；建立不同的地球物理量与Aki-Richards参数之间的关系，得到AVO近似矩阵；对所述敏感识别因子进行微分计算得到纵波速度反射系数、横波速度反射系数以及密度反射系数与敏感识别因子反射系数之间的反射系数关系式；结合所述AVO近似矩阵以及所述反射系数关系式，得到关于所述敏感识别因子、剪切模量以及密度的AVO近似系数矩阵；对所述AVO近似系数矩阵求逆，并代入所述Aki-Richards近似式，得到所述AVO近似式。

5.根据权利要求3所述的砂岩储层识别方法，其特征在于，以所述AVO近似式为核函数直接反演得到所述敏感识别因子，具体为：

计算不同的多个入射角度下各层界面的反射系数矩阵；对所述反射系数矩阵中的参数进行去相关处理；结合所述反射系数矩阵，通过反演获取反射系数目标函数；对所述反射系数目标函数进行求解得到所述敏感识别因子。

6.根据权利要求5所述的砂岩储层识别方法，其特征在于，计算不同的多个入射角度下各层界面的反射系数矩阵，具体为：

计算不同的多个入射角度下各层界面的初始反射系数矩阵；将所述初始反射系数矩阵简化为分块矩阵，得到所述反射系数矩阵。

7.根据权利要求5所述的砂岩储层识别方法，其特征在于，对所述反射系数矩阵中的参数进行去相关处理，具体为：

获取所述敏感识别因子、剪切模量以及密度之间的协方差矩阵；对所述协方差矩阵进行奇异值分解，进而对所述协方差矩阵的特征向量进行逆扩展；结合逆扩展后的矩阵对所述反射系数矩阵进行变换，得到去相关处理后的反射系数矩阵。

8.根据权利要求5所述的砂岩储层识别方法，其特征在于，通过反演获取反射系数目标函数，具体为：

利用贝叶斯公式反演获取初步目标函数，对所述反射系数目标函数进行去相关处理，在所述初步目标函数中加入约束条件得到所述反射系数目标函数。

9.根据权利要求5-8任一所述的砂岩储层识别方法，其特征在于，对所述反射系数目标函数进行求解得到所述敏感识别因子，具体为：

采用迭代重加权最小二乘法对所述反射系数目标函数进行求解，得到所述敏感识别因子。