[发明专利]一种角度弹性参数岩性识别方法和装置

权利要求书

1.一种角度弹性参数岩性识别方法，其特征在于，该角度弹性参数岩性识别方法包括以下步骤：

在常规转换横波弹性阻抗方程的基础上，得到与纵波弹性阻抗同量级的转换横波弹性阻抗方程；

根据所述转换横波弹性阻抗方程以及纵波弹性阻抗方程，得到转换横波弹性阻抗伪测井曲线和纵波弹性阻抗伪测井曲线；

利用转换横波弹性阻抗伪井曲线替代横波阻抗测井曲线、纵波弹性阻抗伪井曲线替代纵波阻抗测井曲线，代入常规弹性参数计算公式，得到某一角度范围内弹性参数伪测井曲线；

利用所述弹性参数伪测井曲线结合岩性测井曲线进行敏感性对比，选取敏感性最强的角度弹性参数对应的弹性参数为目标角度弹性参数，对应的角度为最佳角度；

输入常规多波角度部分叠加地震道集；

对地震道集进行层位进行解释；

采用常规约束稀疏脉冲反演算法进行转换横波弹性阻抗反演，得到带有角度信息的转换横波弹性阻抗数据体、纵波弹性阻抗数据体；

根据测井曲线目标角度弹性参数及最佳角度分析结果，代入转换横波弹性阻抗、纵波弹性阻抗反演数据体到相应弹性参数计算公式中，得到角度弹性参数数据体；

根据角度弹性参数数据体岩性取值范围，进行岩性识别，完成所述角度弹性参数岩性的识别；

与纵波弹性阻抗同量级的转换横波弹性阻抗方程按照如下步骤获得：

引入目的层平均横波速度、密度，对常规转换横波弹性阻抗方程进行标准化，得到标准化后的转换横波弹性阻抗方程：

其中：

θ表示入射角，单位为度；

V_p为纵波速度，单位为米/秒；

V_s为横波速度，单位为米/秒；

ρ为密度，单位为千克/立方米；

V_s0为目的层段横波速度，单位为米/秒；

ρ₀为密度的平均值，单位为千克/立方米。

2.根据权利要求1所述的角度弹性参数岩性识别方法，其特征在于，所述常规多波角度部分叠加地震道集包括纵波角度部分叠加地震道集、转换横波角度部分叠加地震道集，其中，要求叠加道集角度与所述最佳角度相同。

3.根据权利要求1所述的角度弹性参数岩性识别方法，其特征在于，转换横波弹性阻抗数据体SEI(θ)，与常规纵波弹性阻抗数据体EI(θ)，通过组合得到角度弹性参数数据体：P_ang(θ)＝f(SEI(θ),EI(θ))，其中，θ表示入射角，单位为度。

4.根据权利要求3所述的角度弹性参数岩性识别方法，其特征在于，P_ang(θ)＝f(SEI(θ),EI(θ))中，常规弹性参数P的构建关系表述为纵波阻抗I_p与横波阻抗I_s的函数P＝f(I_p,I_s)。

5.根据权利要求1所述的角度弹性参数岩性识别方法，其特征在于，最佳角度按照如下步骤确定：

利用纵波、横波、密度测井曲线，根据转换横波弹性阻抗方程以及常规纵波弹性阻抗方程，分别采用不同角度计算角度弹性参数伪井曲线，结合岩性测井曲线进行敏感性对比，选取敏感性最强的角度弹性参数对应的角度为最佳角度。