CN104267429B - 确定地层压力的方法及装置 - Google Patents

Abstract

本发明公开了一种确定地层压力的方法及装置，其中方法包括：对纵波PP+转换波PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；将纵波阻抗体、横波阻抗体和密度体转换为深度域；根据密度体计算上覆地层压力；根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；根据上覆地层压力和骨架有效应力计算地层压力。采用本发明可以提高地层压力地震预测的精度。

Description

确定地层压力的方法及装置

技术领域

本发明涉及石油地球物理勘探技术领域，尤其涉及确定地层压力的方法及装置。

背景技术

地层压力是石油及天然气勘探、开发中非常重要的一个参数。地层压力预测是油气勘探和开发面临的重大地质问题。在钻探过程中，精准的地层压力预测对保证钻井安全、提高钻探效率和降低钻井成本具有极其重要的作用。近年来发现，地层压力预测对保护储层尤其是天然气储层也具有特别重要的意义。

四十多年来，石油行业逐步形成了一套钻前地层压力预测、监测的基本理论，从基本原理上大体可分为两大类：一类利用各种参数随深度的变化在正常段内建立起压实趋势线，然后根据实测值偏离趋势线的程度来估算地层压力；另一类则不需要直接建立正常压实趋势线，而是建立测量值与地层压力间的经验关系，以判定和估算地层压力。

对于第一类方法，准确的压力预测关键在于建立真实准确的正常压实趋势线，尤其是对于埋深较大地层。但只能通过经验来建立的正常压实趋势线决定了第一类方法的主观不确定性，故此类方法的适用前提是提供足够多的资料建立正常压实趋势线。这类方法主要包括：伊顿法(具体可以参见：Eaton B A.The equation for geopressure prediction from welllogs.SPE，1975，5544)，以及等效深度公式法等。

第二类方法不依赖于正常压实趋势线，通过建立地震地层压力计算模型来预测地层压力，故应用起来比较方便，尤其适用于勘探初期的区块。现有比较流行的地层压力计算模型有如下几种：菲利普恩法(具体可以参见：Fillippone W R.On the prediction of abnormallypressured sedimentary rocks from seismic data.OTC，1979，3662；Fillippone W R.Estimationof formation parameters and the prediction of overpressures from seismic data.SEG，1982，0502)，以及刘震法(具体可以参见：刘震，张万选，张厚福等；辽西凹陷北洼下第三系异常地层压力分析；石油学报，1993，14(1)：14-24)等。

现行地震地层压力计算模型均最终归结到与纵波速度有关的公式，容易造成地层压力地震预测的精度不够、分辨率不高等问题。

发明内容

本发明实施例提供一种确定地层压力的方法，用以提高地层压力地震预测的精度，该方法包括：

对纵波PP+转换波PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

将纵波阻抗体、横波阻抗体和密度体转换为深度域；

根据密度体计算上覆地层压力；

根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

根据上覆地层压力和骨架有效应力计算地层压力；

所述根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力，包括：

获取多口全波列测井的实测压力值；

计算相应的上覆地层压力曲线；

通过最小二乘法拟合公式与H的二阶多项式关系；其中，为单位厚度压缩量，ΔH为厚度压缩量，H为厚度；

利用测井上求取的与H二阶多项式关系，联合纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力，所述修正的骨架有效应力计算公式如下：

$P_e=\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H};$

其中，P_e为骨架有效应力，P_imp ^D为深度域纵波阻抗体，S_imp ^D为深度域横波阻抗体，ρ^D为深度域密度体。

一个实施例中，所述对PP+PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体，包括：

获取保幅处理的n个PP叠前道集的分偏移距部分叠加数据和n个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域弹性参数体：纵波阻抗体、横波阻抗体和密度体，其中，n≥3。

一个实施例中，所述将纵波阻抗体、横波阻抗体和密度体转换为深度域，包括：

获取叠加速度谱资料和地表基准面资料；

利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

一个实施例中，所述深度域是以地表为基准面的埋深域。

一个实施例中，所述根据密度体计算上覆地层压力，包括：通过对密度体的深度积分来计算上覆地层压力。

本发明实施例还提供一种确定地层压力的装置，用以提高地层压力地震预测的精度，该装置包括：

联合反演模块，用于对PP+PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

深度域模块，用于将纵波阻抗体、横波阻抗体和密度体转换为深度域；

上覆地层压力模块，用于根据密度体计算上覆地层压力；

骨架有效应力模块，用于根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

地层压力模块，用于根据上覆地层压力和骨架有效应力计算地层压力；

所述骨架有效应力模块包括：

获取单元，用于获取多口全波列测井的实测压力值；

计算单元，用于计算相应的上覆地层压力曲线；

拟合单元，用于通过最小二乘法拟合公式与H的二阶多项式关系；其中，为单位厚度压缩量，ΔH为厚度压缩量，H为厚度；

骨架有效应力单元，用于利用测井上求取的与H二阶多项式关系，联合纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力，所述修正的骨架有效应力计算公式如下：

$P_e=\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H};$

其中，P_e为骨架有效应力，P_imp ^D为深度域纵波阻抗体，S_imp ^D为深度域横波阻抗体，ρ^D为深度域密度体。

一个实施例中，所述联合反演模块具体用于获取保幅处理的n个PP叠前道集的分偏移距部分叠加数据和n个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域弹性参数体：纵波阻抗体、横波阻抗体和密度体，其中，n≥3。

一个实施例中，所述深度域模块包括：

资料单元，用于获取叠加速度谱资料和地表基准面资料；

模型单元，用于利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

转换单元，用于以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

一个实施例中，所述深度域是以地表为基准面的埋深域。

一个实施例中，所述上覆地层压力模块具体用于通过对密度体的深度积分来计算上覆地层压力。

在本发明实施例提供的技术方案中，通过叠前PP+PS联合同时反演得到高精度的纵波阻抗体、横波阻抗体和密度体，并创造性地转换至深度域，为下一步压力预测奠定了基础，随后根据密度体计算上覆地层压力，根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力，最终求得地层压力，克服了以往只利用纵波速度来预测地层压力的缺陷，引入了高精度的横波阻抗和密度弹性参数体，提高了地层压力地震预测的精度，提高了对薄储层压力预测的分辨率，特别适用于致密气田薄层的地震压力预测。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本发明实施例中确定地层压力的方法实施流程示意图；

图2为本发明实施例中转换为深度域的实施流程示意图；

图3为本发明实施例中获取方法实施流程示意图；

图4为本发明实施例中确定地层压力的装置结构示意图；

图5为本发明实施例中深度域模块的示意图；

图6为本发明实施例中骨架有效应力模块的示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

发明人发现，现行地震地层压力计算模型均最终归结到与纵波速度有关的公式，容易造成地层压力地震预测的精度不够、分辨率不高等问题，是由于现行地震地层压力计算模型不考虑横波速度与密度等的影响，仅涉及纵波速度项而舍弃横波速度与密度项。基于此，在本发明实施例提供的技术方案中，考虑针对致密气田钻前地层压力地震预测问题，采用刘震法地震地层压力计算模型，并做改进，然后在叠前多波多分量地震资料基础上通过叠前多波联合反演获得可靠的高分辨率的弹性参数信息(包括纵波阻抗、横波阻抗和密度体)，并最终形成一套适用于致密气田的基于叠前多波参数的高精度钻前地层压力地震预测方法。下面对实施方式进行说明。

图1为本发明实施例中确定地层压力的方法实施流程示意图，如图1所示，可以包括：

步骤101：对PP(纵波)+PS(转换波)地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

具体实施中，可以获取保幅处理的n(n≥3)个PP叠前道集的分偏移距部分叠加数据和n(n≥3)个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域P_imp ^T(纵波阻抗体)、S_imp ^T(横波阻抗体)和ρ^T(密度体)等弹性参数体。

步骤102：将纵波阻抗体、横波阻抗体和密度体转换为深度域；

对步骤101中所得弹性参数体时深转换，得到深度域P_imp ^D(纵波阻抗体)、S_imp ^D(横波阻抗体)和ρ^D(密度体)等弹性参数体；

图2为本发明实施例中转换为深度域的实施流程示意图，如图2所示，可以包括：

步骤201：获取叠加速度谱资料、地表基准面资料；

步骤202：利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

实施中，也即：首先利用多口已知的测井进行精细井震标定，产生速度模型，用其来标定利用速度谱建立的速度模型；

步骤203：以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

具体的，实施例中的时深转换的含义与常规的时深转换不同：常规时深转换一般以海平面为最终基准面，得到的是海拔深度；而实施例中时深转换是以地表为基准面，得到埋深深度数据体。时间域数据体转换为埋深深度体为后续的上覆岩层压力计算提供了基础。

步骤103：根据密度体计算上覆地层压力；

也即，根据密度体计算上覆地层压力P_ov。

具体实施中，可以对密度体ρ^D进行深度积分，得到上覆岩层压力体P_ov：

$P_{ov}={\∫}_0^{H_0}{\mathrm{\ρ}}^D\left(H\right)gdh---\left(1\right)$

其中H₀为从地表算起的观察点深度。

步骤104：根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

也即：根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力P_e。

具体实施中，刘震等提出的地震地层压力计算模型Ⅰ(SFPM-Ⅰ)中骨架有效应力计算公式为(具体可以参见：刘震，张万选等，辽西凹陷北洼下第三系异常地层压力分析，石油学报，1993，14(1)：14-23)：

$P_e=\frac{1+\mathrm{\σ}}{3-3\mathrm{\σ}}\·\mathrm{\ρ}\·{V_p}^2\·\frac{\mathrm{\Δ}H}{H}---\left(2\right)$

其中σ为泊松比，ρ为密度，V_p为纵波速度，为单位厚度压缩量，ΔH为厚度压缩量，H为厚度。

刘震把此公式中σ取为0.12的常数值，并最终归结为纵波速度V_p的一元二次方程，为此在此模型的基础上，本发明实施例将σ保留，并将其变形为：

$P_e=\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H}---\left(3\right)$

利用此变形后的公式进行骨架有效应力P_e计算。

步骤105：根据上覆地层压力和骨架有效应力计算地层压力。

也即：根据上覆地层压力和骨架有效应力计算地层压力P_p。

具体实施中，计算地层压力计算公式可以为：

$P_p=P_{ov}-\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H}---\left(4\right)$

其中，为一个受深度、岩性和沉积埋藏史等多种因素有关的变量。

图3为本发明实施例中获取方法实施流程示意图，如图3所示，可以包括：

步骤301：获取多口全波列测井的实测压力值；

步骤302：计算井上相应的上覆地层压力曲线；

步骤303：通过最小二乘法拟合与H的二阶多项式关系。

具体实施中，将井上已知的P_p、P_imp ^D、S_imp ^D、ρ^D和计算得到的P_ov代入到步骤105中的公式中(4)即可以求得该点出的值，然后再利用该地区的实测地层压力资料和求得的值，通过最小二乘法即可拟合与H的二阶多项式关系为：

$\frac{\mathrm{\Δ}H}{H}={\mathrm{aH}}^2+bH+c---\left(5\right)$

其中，a、b、c为常数。

步骤304：利用测井上求取的与H二阶多项式关系，联合求取的纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力。

具体实施中，地层压力计算公式可以为：

相应的地层压力系数计算公式为：

具体实施时，还可以利用实测钻井压力资料与计算的地层压力结果分析与标定，对计算结果进行修正。

具体实施时，还可以绘制地层压力/压力系数剖面图。

下面结合鄂尔多斯盆地的实例说明本发明实施例提供的技术方案的运用：致密气田叠前多波参数法进行钻前地层压力地震预测。

鄂尔多斯盆地SLG地区为典型的低孔、低渗致密气田，主力气层为上古生界二叠系的盒8、山1砂岩层，为强非均质性的薄储层，该区古生界地层压力系数一般在0.76-0.95之间，属于异常低压。将本发明实施例提供的技术方案应用于鄂尔多斯盆地SLG地区，其具体的实现步骤可以包括：

1、叠前PP+PS联合同时反演

选定多波测线S077053进行叠前PP+PS联合同时反演，选取各7个变偏移距叠加的PP和PS道集的部分叠加数据，在PP叠前反演获得准确的纵横波速度比的基础上，将此结果作为模型，将PS数据从PS域转换到PP时间域，再用标志反射层对转换波数据进行层位压缩，获得最佳的纵波、转换波匹配剖面；然后分别输入PP和PS部分叠加地震数据体和对应的角度子波，给出不同数据(纵波阻抗、密度、纵横波速度比)的纵向变化趋势以及横向上的约束范围，通过质量控制优选出合适的一组参数。按照求解Knott—Zoeppritz方程对PP波、PS波进行联合同时反演，最终获得更加稳定的高分辨率纵波阻抗、横波阻抗和密度数据体：S077053叠前PP+PS联合同时反演得到的纵波阻抗、横波阻抗和密度数据体，其中CDP7488SU69井处，盒8段储层纵波阻抗中低值(11182g/cc*m/s)、横波阻抗中高值(6145.6g/cc*m/s)、密度中低值(2.46g/cc)。

2、弹性参数体转至埋深域

首先利用过线测井SU69、SU74、SU165等进行精细井震标定，产生速度模型，用其来标定利用速度谱建立的速度模型，然后以地表为基准面，将纵波阻抗、横波阻抗和密度体转换至埋深域，获得S077053埋深域纵波阻抗、横波阻抗和密度数据体。

3、多井拟合与H的二阶多项式关系

利用本区内全波列测井SHEN15、MI35、SU151、SU48、SU163、SU69、SD51-42目的层盒8的地层压力资料和已知或者计算的弹性参数P_p、P_imp ^D、S_imp ^D、ρ^D、P_ov，求得这7口井点处的值，然后再利用实测地层压力资料和7口井的值，通过最小二乘法拟合，并利用实测压力与计算的地层压力结果分析与标定，对计算结果进行修正，最终得到本地区与H的二阶多项式关系为：

$\frac{\mathrm{\Δ}H}{H}=-8.957745\×{10}^{-14}{H}^2+5.7128389325\×{10}^{-10}H-1.1755738181174\×{10}^{-7}---\left(8\right)$

4、求取地层压力/压力系数剖面

联合公式(6)、(7)和公式(8)，即可求取地层压力/压力系数剖面，其中D取经验参数-27MPa。可以得到，S077053地震预测的地层压力及压力系数剖面中，CDP7488处预测的地层压力及压力系数分别为29.77MPa和0.92；对应的SU69井处实测地层压力及压力系数分别为29.52MPa和0.91，井点处与预测的吻合较好，预测结果较为理想，而且预测的二维剖面分辨率较高，提升了薄储层压力预测的实际效果，本发明实施例提供的技术方案特别适用于致密气田薄层的地震压力预测。

基于同一发明构思，本发明实施例中还提供了一种确定地层压力的装置，由于该装置解决问题的原理与确定地层压力的方法相似，因此该装置的实施可以参见方法的实施，重复之处不再赘述。

图4为本发明实施例中确定地层压力的装置结构示意图，如图4所示，装置中可以包括：

联合反演模块401，用于对PP+PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

深度域模块402，用于将纵波阻抗体、横波阻抗体和密度体转换为深度域；

上覆地层压力模块403，用于根据密度体计算上覆地层压力；

骨架有效应力模块404，用于根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

地层压力模块405，用于根据上覆地层压力和骨架有效应力计算地层压力。

实施中，联合反演模块具体可以用于获取保幅处理的n个PP叠前道集的分偏移距部分叠加数据和n个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域弹性参数体：纵波阻抗体、横波阻抗体和密度体，其中，n≥3。

图5为本发明实施例中深度域模块的示意图，如图5所示，深度域模块可以包括：

资料单元501，用于获取叠加速度谱资料和地表基准面资料；

模型单元502，用于利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

转换单元503，用于以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

实施中，所述深度域是以地表为基准面的埋深域。

实施中，上覆地层压力模块具体可以用于通过对密度体的深度积分来计算上覆地层压力。

图6为本发明实施例中骨架有效应力模块的示意图，如图6所示，骨架有效应力模块可以包括：

获取单元601，用于获取多口全波列测井的实测压力值；

计算单元602，用于计算相应的上覆地层压力曲线；

拟合单元603，用于通过最小二乘法拟合公式与H的二阶多项式关系；

骨架有效应力单元604，用于利用测井上求取的与H二阶多项式关系，联合纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力。

由上述实施例可见，在本发明实施例提供的技术方案中，通过叠前PP+PS联合同时反演得到高精度的纵波阻抗体、横波阻抗体和密度体，并创造性地转换至深度域，为下一步压力预测奠定了基础，随后根据密度体计算上覆地层压力，根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力，最终求得地层压力，克服了以往只利用纵波速度来预测地层压力的缺陷，引入了高精度的横波阻抗和密度弹性参数体，提高了地层压力地震预测的精度，提高了对薄储层压力预测的分辨率，特别适用于致密气田薄层的地震压力预测。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种确定地层压力的方法，其特征在于，包括：

对纵波PP+转换波PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

将纵波阻抗体、横波阻抗体和密度体转换为深度域；

根据密度体计算上覆地层压力；

根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

根据上覆地层压力和骨架有效应力计算地层压力；

根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力，包括：

获取多口全波列测井的实测压力值；

计算相应的上覆地层压力曲线；

通过最小二乘法拟合公式与H的二阶多项式关系；其中，为单位厚度压缩量，ΔH为厚度压缩量，H为厚度；

利用测井上求取的与H二阶多项式关系，联合纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力，所述修正的骨架有效应力计算公式如下：

$P_e=\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H};$

其中，P_e为骨架有效应力，P_imp ^D为深度域纵波阻抗体，S_imp ^D为深度域横波阻抗体，ρ^D为深度域密度体。

2.如权利要求1所述的方法，其特征在于，所述对PP+PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体，包括：

获取保幅处理的n个PP叠前道集的分偏移距部分叠加数据和n个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域弹性参数体：纵波阻抗体、横波阻抗体和密度体，其中，n≥3。

3.如权利要求1所述的方法，其特征在于，所述将纵波阻抗体、横波阻抗体和密度体转换为深度域，包括：

获取叠加速度谱资料和地表基准面资料；

利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

4.如权利要求1或3所述的方法，其特征在于，所述深度域是以地表为基准面的埋深域。

5.如权利要求1所述的方法，其特征在于，所述根据密度体计算上覆地层压力，包括：通过对密度体的深度积分来计算上覆地层压力。

6.一种确定地层压力的装置，其特征在于，包括：

联合反演模块，用于对PP+PS地震数据进行保幅处理获取多个部分叠加数据，并联合测井资料，通过叠前PP+PS联合同时反演方法求取纵波阻抗体、横波阻抗体和密度体；

深度域模块，用于将纵波阻抗体、横波阻抗体和密度体转换为深度域；

上覆地层压力模块，用于根据密度体计算上覆地层压力；

骨架有效应力模块，用于根据纵波阻抗体、横波阻抗体和密度体计算骨架有效应力；

地层压力模块，用于根据上覆地层压力和骨架有效应力计算地层压力；

所述骨架有效应力模块包括：

获取单元，用于获取多口全波列测井的实测压力值；

计算单元，用于计算相应的上覆地层压力曲线；

拟合单元，用于通过最小二乘法拟合公式与H的二阶多项式关系；其中，为单位厚度压缩量，ΔH为厚度压缩量，H为厚度；

骨架有效应力单元，用于利用测井上求取的与H二阶多项式关系，联合纵波阻抗体、横波阻抗体和密度体，利用修正的骨架有效应力计算公式计算骨架有效应力，所述修正的骨架有效应力计算公式如下：

$P_e=\frac{\left({P_{imp}^2}^D-\frac{4}{3}{S_{imp}^2}^D\right)}{{\mathrm{\ρ}}^D}\·\frac{\mathrm{\Δ}H}{H};$

其中，P_e为骨架有效应力，P_imp ^D为深度域纵波阻抗体，S_imp ^D为深度域横波阻抗体，ρ^D为深度域密度体。

7.如权利要求6所述的装置，其特征在于，所述联合反演模块具体用于获取保幅处理的n个PP叠前道集的分偏移距部分叠加数据和n个PS叠前道集的分偏移距部分叠加数据，进行PP和PS数据对齐后，进行叠前PP+PS联合同时反演，得到时间域弹性参数体：纵波阻抗体、横波阻抗体和密度体，其中，n≥3。

8.如权利要求6所述的装置，其特征在于，所述深度域模块包括：

资料单元，用于获取叠加速度谱资料和地表基准面资料；

模型单元，用于利用标定后的井速度模型校正用速度谱建立的速度模型，产生最终速度模型；

转换单元，用于以地表为基准面，将纵波阻抗体、横波阻抗体和密度体转换为埋深域。

9.如权利要求6或8所述的装置，其特征在于，所述深度域是以地表为基准面的埋深域。

10.如权利要求6所述的装置，其特征在于，所述上覆地层压力模块具体用于通过对密度体的深度积分来计算上覆地层压力。