[发明专利]一种基于岩性-物性分类的碎屑岩地层压力计算方法

权利要求书

1.一种基于岩性-物性分类的碎屑岩地层压力计算方法，其特征在于，所述的方法包括如下步骤：

S1，对碎屑岩地层进行岩性和物性分类，划分为泥岩层和不同物性条件的砂岩层；

S2，根据测井解释渗透率数据，将砂岩层划分为极低渗低渗砂岩层和中高渗砂岩层；

S3，建立电缆测压-深度图，对测压数据进行分段线性拟合；拟合直线近似平行于上覆压力时，砂岩为孤立型砂岩层；拟合直线近似平行于静水压力时，砂岩为连通型砂岩层；结合步骤S2的结果，将中高渗砂岩层分为孤立型和连通型；

S4，根据文献调研，计算泥岩层的地层压力；输入测井声波时差的数据，采用等效深度法、伊顿法、Bowers法、经验法计算泥岩地层孔隙压力；

S5，计算极低渗低渗砂岩层和孤立型中高渗砂岩层的地层压力；输入测井声波时差的数据，采用上下相邻泥岩声波时差的平均值或线性内插结果，采用的地层压力计算方法与步骤S4中泥岩层相同；

S6，对连通型中高渗砂岩采用压力侧向传递法计算地层压力，包括开展地震资料构造解释，追踪砂岩顶底界面，编制砂岩顶底界面的深度构造图，在深度构造图中读取砂岩的最大埋深和最小埋深；

S7，不同岩性、物性地层计算地层压力按深度合并；同一深度存在多个计算压力数据时，连通型中高渗砂岩层计算压力优先级最高，采用连通型中高渗砂岩层计算压力，从而得到全部地层的地层压力。

2.根据权利要求1所述的一种基于岩性-物性分类的碎屑岩地层压力计算方法，其特征在于，所述步骤S4、步骤S5、步骤S6中采用的数据还包括声波速度、电阻率和VSP速度；步骤S5、步骤S6中所述声波速度、电阻率和VSP速度分别采用相邻泥岩层的声波时差、电阻率、地震层速度值的均值。

3.根据权利要求2所述的一种基于岩性-物性分类的碎屑岩地层压力计算方法，其特征在于，所述伊顿法的计算公式为：

式中，

PP为孔隙压力梯度，g/cm³；

OBG为上覆压力梯度，g/cm³；

PP_N为静水压力梯度，g/cm³；

DT_N为正常压实地层声波时差，us/ft；

DT_o为实测地层声波时差，us/ft；

x为伊顿指数；

所述等效深度法的计算公式为：

PP＝PP_E+(OBG-OBG_E)

式中，

PP为孔隙压力，PSI或Pa；

PP_E为等效深度点的孔隙压力，PSI或Pa；

OBG为上覆压力，PSI或Pa；

OBG_E为等效深度点的上覆压力，PSI或Pa；

所述Bowers法的计算公式为：

式中，

PP为孔隙压力梯度，g/cm³；

OBG为上覆压力梯度，g/cm³；

V₀为海底沉积物速度，m/s；

V为实测地层速度，m/s；

A和B为经验常数；

Depth为垂深，m；

所述经验法的计算公式为：

PP＝F(V-V_N)

式中，

PP为孔隙压力梯度，g/cm³；

V_N为正常压实地层速度，m/s；

V为实测地层速度，m/s；

F为速度与测量压力梯度间的经验函数。

4.根据权利要求3所述的一种基于岩性-物性分类的碎屑岩地层压力计算方法，其特征在于，所述步骤S6中连通型中高渗砂岩压力计算公式为：

假设连通型中高渗砂岩的深度中值位置为H，

H＝(H_max+H_min)/2时砂岩孔隙压力与泥岩孔隙压力相等，

即P_SandH＝P_ShaleH，

采用压力侧向传递法计算连通型中高渗砂岩的地层孔隙压力，砂岩最小埋深H_min位置的孔隙压力:

P_SandHmin＝P_ShaleH-Kρg(H_max-H_min)/2，

砂岩最大埋深Hmax位置的孔隙压力:

P_SandHmax＝P_ShaleH+Kρg(H_max-H_min)/2，

式中，ρ为砂岩流体密度，g/cm³；g为重力加速度，9.81m/s²；P为地层孔隙压力，MPa；H为深度，m；K为单位转换系数，1×10^-3。