一种基于地震资料的地层孔隙压力预测方法与流程

技术特征：

1.一种基于地震资料的地层孔隙压力预测方法，其特征在于：所述方法基于深度域测井数据Well(z)和叠后时间域地震数据Seis(x,y,t)实现对地层孔隙压力P_pore(x,y,t)的预测，其中z代表深度域，单位为m或者ft，x,y代表每个地震道对应的空间位置，是x,y坐标或者是纵横测线号，t代表时间域，单位为s或者ms。

2.根据权利要求1所述的基于地震资料的地层孔隙压力预测方法，其特征在于：所述方法包括以下步骤：

(1)基于叠后时间域地震数据Seis(x,y,t)，通过叠后波阻抗反演，获取纵波阻抗信息Ip(x,y,t)；

(2)基于深度域测井数据Well(z)，通过时深转换和数据分析，获取基于测井资料的地层有效应力P_diff(t)和纵波阻抗Ip(t)的关系P_diff(t)＝g(Ip(t))；

(3)基于步骤(1)获得的纵波阻抗Ip(x,y,t)和步骤(2)获得的关系式P_diff(t)＝g(Ip(t))，计算地层有效应力P_diff(x,y,t)；

(4)基于步骤(1)获得的纵波阻抗Ip(x,y,t)，确定上覆地层压力P₀(x,y,t)；

(5)基于步骤(3)获得的地层有效应力P_diff(x,y,t)和步骤(4)获得的上覆地层压力P₀(x,y,t)，计算得到地层孔隙压力P_pore(x,y,t)。

3.根据权利要求2所述的基于地震资料的地层孔隙压力预测方法，其特征在于：所述步骤(2)是这样实现的：

深度域测井数据Well(z)要求包含如下测井曲线，①从声波时差数据获取的纵波速度vp(z)；②测井解释的岩石密度ρ(z)；③测井解释的孔隙压力p_pore(z)；

通过如下关系式将深度域测井数据转换到与地震匹配的时间域数据

$t={\∫}_0^z\frac{2\mathrm{dz}}{\mathrm{vp}\left(z\right)}$

经过时深转换后，获得如下时间域的测井数据：①时间域的纵波速度曲线vp(t)；②时间域的密度曲线ρ(t)；③时间域的孔隙压力p_pore(t)；

通过速度和密度的乘积计算得到基于测井资料的纵波阻抗：

Ip(t)＝ρ(t)*vp(t)

基于测井资料的上覆地层压力p₀(t)实际上等效于上覆地层的重力，通过如下计算公式计算得到：

$p_0\left(t\right)={\∫}_0^z\mathrm{\ρ}\left(z\right)\mathrm{gdz}={\∫}_0^t\mathrm{\ρ}\left(t\right)g\frac{\mathrm{vp}\left(t\right)}{2}\mathrm{dt}$

用所述上覆地层压力p₀(t)同测井解释得到的孔隙压力p_pore(t)通过减法运算得到基于测井资料的地层有效应力p_diff(t)：

p_diff(t)＝p₀(t)-p_pore(t)

将纵波阻抗Ip(t)表示成地层有效应力p_diff(t)的函数：

Ip(t)＝f(P_diff(t))＝Ip₀+cP_diff(t)^b

对方程进行重新推导，将p_diff(t)表示成Ip(t)的函数：

P_diff(t)＝g(Ip(t))＝((Ip(t)-Ip₀)/c)^1/b

这里，Ip₀、c和b是常数，需要根据实际测井数据经过交会图分析和线性回归的方法获得。

4.根据权利要求3所述的基于地震资料的地层孔隙压力预测方法，其特征在于：所述步骤(3)是这样实现的：

计算得到与波阻抗数据范围大小一致的地层有效应力数据体：

P_diff(x,y,t)＝g(Ip(x,y,t))＝((Ip(x,y,t)-Ip₀)/c)^1/b。

5.根据权利要求4所述的基于地震资料的地层孔隙压力预测方法，其特征在于：所述步骤(4)是这样实现的：

用纵波阻抗的积分形式表示上覆地层压力：

$P_0(x,y,t)={\∫}_0^t\frac{1}{2}\mathrm{Ip}(x,y,t)\mathrm{\ξdt}$

这里ξ是一个常数，代表重力加速度。

6.根据权利要求5所述的基于地震资料的地层孔隙压力预测方法，其特征在于：所述步骤(5)是这样实现的：

利用上覆地层压力P₀(x,y,t)和地层有效应力P_diff(x,y,t)之间的差获得地层孔隙压力：

P_pore(x,y,t)＝P₀(x,y,t)-P_diff(x,y,t)。