[发明专利]页岩裂缝型储层稳定反演预测方法及系统

权利要求书

1.一种页岩裂缝型储层稳定反演预测方法，其特征在于，包括：

将含一组垂直或近似垂直裂缝的页岩储层等效为正交各向异性介质，获取页岩裂缝型储层的宽方位地震数据；

基于OA介质三参数的PP波反射系数方程，推导得到不同观测方位的弹性阻抗方程组；

使用贝叶斯反演框架，求取后验概率密度函数的最大后验概率估计，同时加入平滑背景参数模型约束，建立目标函数，对不同观测方位的弹性阻抗体分别进行反演；

基于反演得到的弹性阻抗，以及不同观测方位的弹性阻抗方程组，计算得到纵波阻抗参数以及含各向异性信息的组合参数。

2.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，所述页岩裂缝型储层的宽方位地震数据，包括：地震资料的方位角和入射角。

3.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，基于OA介质三参数的PP波反射系数方程，推导得到不同观测方位的弹性阻抗方程组，具体为：

定义；

其中，表示OA介质的PP波反射系数，表示下覆地层的弹性阻抗，表示上伏地层的弹性阻抗，表示上下两套地层弹性阻抗的均值；

将基于OA介质三参数的PP波反射系数方程代入上式，依次经过两侧取积分、标准化处理以及线性化处理后，得到弹性阻抗方程；

使用不同的方位角和入射角，构建不同观测方位的弹性阻抗方程组。

4.如权利要求3所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，所述的弹性阻抗方程，具体为：

其中，A为纵波阻抗参数，B和C为含各向异性信息的组合参数；A₀、B₀、C₀分别是组合参数A、B、C的平均值；M用来使的尺度与声阻抗一致，并保证在时，；为入射角， 为方位角；表示地层的纵波速度、表示地层的密度。

5.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，使用贝叶斯反演框架，求取后验概率密度函数的最大后验概率估计，具体为：

对地震记录的褶积模型进行去相关处理；

使用贝叶斯反演框架，将柯西分布作为先验概率密度函数，将高斯分布作为似然函数，得到后验概率密度函数；

将先验概率密度函数与全概率定理结合，求取后验概率密度函数的最大后验概率估计。

6.如权利要求5所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，所述平滑背景参数模型约束具体为：将测井上的低频趋势补充到地震所缺失的低频信息中。

7.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，所述建立目标函数，具体为：

其中，N代表采样点个数， 为入射角为，方位角为的弹性阻抗反射系数；和分别代表实际和合成地震记录； 为测井计算的弹性阻抗梯度；和分别代表p阶和q阶范数；常数和分别是权重系数。

8.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，所述对不同观测方位的弹性阻抗体分别进行反演，具体为：

目标函数的最小值即为最优弹性阻抗反射系数估计，进而求得最终的弹性阻抗：

其中，t₀代表初始时刻，t_i代表第i时刻；表示初始时刻的弹性阻抗，表示时刻的PP波反射系数。

9.如权利要求1所述的一种页岩裂缝型储层稳定反演预测方法，其特征在于，将计算得到的纵波阻抗参数以及含各向异性信息的组合参数，用于对页岩储层开展裂缝预测和油气检测。

10.一种页岩裂缝型储层稳定反演预测系统，其特征在于，包括：

数据获取模块，用于将含一组垂直或近似垂直裂缝的页岩储层等效为正交各向异性介质，获取页岩裂缝型储层的宽方位地震数据；

弹性阻抗方程组构建模块，基于OA介质三参数的PP波反射系数方程，推导得到不同观测方位的弹性阻抗方程组；

弹性阻抗反演模块，用于使用贝叶斯反演框架，求取后验概率密度函数的最大后验概率估计，同时加入平滑背景参数模型约束，建立目标函数，对不同观测方位的弹性阻抗分别进行反演；

参数求解模块，用于基于反演得到的弹性阻抗，以及不同观测方位的弹性阻抗方程组，计算得到纵波阻抗参数以及含各向异性信息的组合参数。