[发明专利]基于叠前地震数据的地层孔隙压力预测方法、装置及电子设备

权利要求书

1.一种基于叠前地震数据的地层孔隙压力预测方法，其特征在于，所述方法包括：

基于待反演弹性参数建立所述待反演弹性参数的反演模型；

获取叠前地震数据，根据所述叠前地震数据及所述反演模型以构建所述反演模型的目标函数，并对所述目标函数进行偏导计算以得到所述待反演弹性参数的改变量；

对所述改变量进行迭代以更新所述改变量，得到所述待反演弹性参数的反演结果，其中，所述待反演弹性参数包括横波速度、纵波速度以及密度数据；

根据反演后得到的密度数据计算上覆地层压力；

根据反演后得到的横波速度、纵波速度及密度数据计算垂直有效应力；

根据上覆地层压力、垂直有效应力及地层孔隙压力之间的关系，并基于得到的所述上覆地层压力和所述垂直有效应力计算地层孔隙压力。

2.根据权利要求1所述的基于叠前地震数据的地层孔隙压力预测方法，其特征在于，所述基于待反演弹性参数建立所述待反演弹性参数的反演模型的步骤，包括：

利用精确Zoeppritz方程基于待反演弹性参数以建立如下的待反演弹性参数的反演模型：

其中V_p1，V_p2，V_s1，V_s2，ρ₁，ρ₂分别表示上下层的纵波速度、横波速度和密度；α，α'分别为纵波的入射角和透射角；β，β'分别为横波的入射角和透射角；R_PP，R_PS为纵波和横波的反射系数；T_PP，T_PS分别为纵波和横波的透射系数。

3.根据权利要求1所述的基于叠前地震数据的地层孔隙压力预测方法，其特征在于，获取叠前地震数据，根据所述叠前地震数据及所述反演模型以构建所述反演模型的目标函数，并对所述目标函数进行偏导计算以得到所述待反演弹性参数的改变量的步骤，包括：

根据所述叠前地震数据及所述反演模型中的待反演弹性参数的模型数据之间的差值以及所述待反演弹性参数建立如下的目标函数：

设定所述目标函数的一阶偏导数为零，以得到如下所示的所述待反演弹性参数的改变量：

Δm_k＝(J(m_k)^TΣ^-1J(m_k)+μ_kI+λ_kL^TL)^-1[g_k-λ_kL^TLm_k]

其中，m＝[V_p1,V_p2,…,V_pN,V_s1,V_s2,…,V_sN,ρ₁,ρ₂,...,ρ_N]为待反演弹性参数，分别代表了纵波速度、横波速度和密度数据；f_k为叠前地震数据和待反演弹性参数的模型数据之间的差值；J_k为Jacobian矩阵，L为一阶导数平滑算子，μ_k和λ_k为与叠前地震数据和模型数据之间误差的平方和||f_k||²有关的参数。

4.根据权利要求1所述的基于叠前地震数据的地层孔隙压力预测方法，其特征在于，所述根据反演后得到的密度数据计算上覆地层压力的步骤，包括：

将反演后得到的密度数据按预设深度间隔数据进行积分叠加以得到上覆地层压力，如下：

其中，P_c为上覆地层压力，ρ(z)为深度Z处的地层密度，g为重力加速度。