[发明专利]一种页岩储层脆性和各向异性参数预测方法及系统

权利要求书

1.一种页岩储层脆性和各向异性参数预测方法，其特征在于，包括：

采集工区内的测井数据，并基于采集的测井数据确定地震子波；

构建含三属性参数的VTI介质反射系数方程，确定三属性参数；

基于测井数据和三属性参数获取先验分布函数、初始三属性参数模型；

根据初始三属性参数模型、地震子波确定反演残差；

基于先验分布函数和似然函数构建最大后验概率意义下的反演目标函数；

根据反演目标函数和反演残差确定最优三属性参数模型；

依据最优三属性参数模型，预测页岩储层脆性、各向异性参数；

其中，三属性参数包括：第一属性参数第二属性参数第三属性参数其中，E为杨氏模量，ρ为密度，σ为泊松比，δ和ε为弱各向异性参数，g为背景弹性特征。

2.如权利要求1所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述构建含三属性参数的VTI介质反射系数方程的步骤为：

利用弹性波逆散射理论、波恩近似及稳相法，推导出基于纵横波模量的VTI介质反射系数方程；

利用纵横波模量与杨氏模量和泊松比间的定量关系，将基于纵横波模量的VTI介质反射系数方程等价改写为基于岩石力学参数的VTI介质反射系数方程；

通过对基于岩石力学参数的VTI介质反射系数方程进行整合化简，获得含三属性参数的VTI介质反射系数方程。

3.如权利要求1所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述确定地震子波的步骤为：

采集工区内的地震叠前道集、测井数据、实际井旁角度域地震道集；

基于地震叠前道集和测井数据，采取统计方法提取依赖于入射角度的地震子波；

将测井数据为输入模型，利用VTI介质反射系数方程正演模拟角度域的PP波道集；

将PP波道集与实际井旁角度域地震道集对比，得到振幅缩放因子；

将振幅缩放因子应用于所述的地震子波。

4.如权利要求3所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述先验分布函数的获取过程具体为：基于测井数据和三属性参数，统计三属性参数的先验信息，形成先验分布函数。

5.如权利要求4所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述三属性参数的先验信息包括三属性参数的均值和三属性参数的协方差矩阵。

6.如权利要求1所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述初始三属性参数模型是通对测井数据进行平滑处理并沿解释层位外推得到的。

7.如权利要求1所述的一种页岩储层脆性和各向异性参数预测方法，其特征在于，所述根据初始三属性参数模型、地震子波确定反演残差的具体步骤为：

获取实际角道域地震道集；

将所述初始三属性参数模型与敏感度矩阵做乘积得到反射系数矩阵，并与地震子波矩阵进行褶积，得到角度域的角道集；

将实际角道域地震道集与所述角度域的角道集作差，得到反演残差。

8.一种页岩储层脆性和各向异性参数预测系统，其特征在于，包括：

三属性参数确定模块，用于构建含三属性参数的VTI介质反射系数方程，确定三属性参数；

地震子波确定模块，用于基于构建的含三属性参数的VTI介质反射系数方程确定地震子波；

先验分布函数和初始三属性参数模型获取模块，用于基于三属性参数获取先验分布函数、初始三属性参数模型；

反演残差确定模块，用于根据初始三属性参数模型、地震子波确定反演残差；

反演目标函数构建模块，用于基于先验分布函数和似然函数构建最大后验概率意义下的反演目标函数；

最优三属性参数模型确定模块，用于根据反演目标函数和反演残差确定最优三属性参数模型；

预测模块，用于依据最优三属性参数模型，预测页岩储层脆性、各向异性参数；

其中，三属性参数包括：第一属性参数第二属性参数第三属性参数其中，E为杨氏模量，ρ为密度，σ为泊松比，δ和ε为弱各向异性参数，g为背景弹性特征。