[发明专利]三维页岩储层孔隙压力预测方法、装置和电子设备

权利要求书

1.一种三维页岩储层孔隙压力预测方法，其特征在于，包括：

获取待预测区域的密度测井曲线、纵波速度测井曲线、页岩储层实测孔隙压力值和三维地震偏移数据体；

根据所述密度测井曲线和所述纵波速度测井曲线，确定所述待预测区域的页岩储层的上覆地层压力和纵波模量；

根据所述上覆地层压力和所述页岩储层实测孔隙压力值，确定页岩储层等效应力；

将所述页岩储层实测孔隙压力值作为约束条件，通过多元线性回归的方式拟合所述纵波速度测井曲线、所述页岩储层等效应力和所述纵波模量，得到拟合参数；

根据所述密度测井曲线和所述纵波速度测井曲线，利用所述三维地震偏移数据体，通过测井约束波阻抗反演得到所述待预测区域的波阻抗数据体、密度数据体和纵波速度数据体；

根据所述波阻抗数据体、所述密度数据体和所述纵波速度数据体，确定储层上覆地层压力数据体和纵波模量数据体；

根据所述拟合参数、所述纵波速度数据体和所述纵波模量数据体，确定页岩储层等效应力数据体；

根据所述储层上覆地层压力数据体和所述页岩储层等效应力数据体，确定页岩储层孔隙压力；

将所述页岩储层实测孔隙压力值作为约束条件，通过多元线性回归的方式拟合所述纵波速度测井曲线、所述页岩储层等效应力和所述纵波模量，得到拟合参数的步骤，包括：

建立方程V_P0(z)＝a₀+a₁S_eff0(z)+a₂M₀(z)；其中，V_P0(z)为所述纵波速度测井曲线，a₀、a₁和a₂为拟合参数；S_eff0(z)为所述页岩储层等效应力、M₀(z)为所述纵波模量；

通过多元线性回归的方式求解所述方程，得到所述拟合参数a₀、a₁和a₂的值。

2.根据权利要求1所述的方法，其特征在于，根据所述密度测井曲线和所述纵波速度测井曲线，确定所述待预测区域的页岩储层的上覆地层压力和纵波模量的步骤，包括：

计算所述待预测区域的页岩储层的上覆地层压力其中，i为样点序号；z_i为第i个样点对应的埋深；g为重力加速度；△h为相邻样点的间距；ρ₀(j)为第j个样点处的密度测井值；

计算纵波模量其中，ρ₀(z_i)为z_i埋深处的密度测井值；V_P0(z_i)为z_i埋深处的纵波速度测井值。

3.根据权利要求1所述的方法，其特征在于，根据所述上覆地层压力和所述页岩储层实测孔隙压力值，确定页岩储层等效应力的步骤，包括：

计算页岩储层等效应力S_eff0(z_i)＝S_V0(z_i)-P₀(z_i)；其中，S_V0(z_i)为所述上覆地层压力；P₀(z_i)为所述页岩储层实测孔隙压力值；z_i为第i个样点对应的埋深。