用于对钻井泥岩地层压力进行预测的方法与流程

[0001]
本发明属于油气勘探技术领域，具体涉及一种用于对钻井泥岩地层压力进行预测的方法。


背景技术：

[0002]
在含油气盆地中，如果地下特定深度处地层压力超过该深度静水压力，则称之为超压。超压与油气成藏、储量评价密切相关，同时决定了安全钻井窗口，因而准确预测地层压力对油气运移研究和钻井工程实施都十分重要。目前常用的地层压力预测方法都是根据泥岩孔隙度变化特征进行的，并且所有的压力预测均需要建立描述孔隙度随深度的变化的趋势，即“正常压实曲线”。泥岩压实过程中，随着埋深增加，孔隙流体能够顺利排出，孔隙度随之不断减小，地层压力保持静水压力，该过程称之为正常压实。利用能够反映泥岩孔隙度变化的测井响应开展压实状态分析已经成为泥岩压实研究及地层压力预测的主要手段。“正常压实曲线”定义了正常压力条件下泥岩孔隙度或其测井响应与深度的关系，与之相对应的孔隙度/测井响应与垂直有效应力的关系称之为加载曲线。建立可靠的正常压实曲线和加载曲线是地层压力预测及超压成因判识等研究的基础。
[0003]
研究人员对泥岩机械压实作用机理及岩石物理响应开展了大量研究，提出利用athy公式来表征泥岩孔隙度随深度变化的定量关系，并且将这种关系应用于地层压力预测及超压成因机制判识。传统上对于泥岩压实的认识仅仅基于单一的机械压实作用。然而，这种基于单一的机械压实曲线进行地层压力预测的方法在应用于中深层地层压力预测及超压成因机制分析时存在严重的误差。随着埋深增大、地层温度升高，中深层泥岩岩石物理性质及其测井响应特征相对于浅层泥岩岩石物理性质及其测井响应特征存在显著的差异，进而导致中深层泥岩压实的孔隙度-有效应力关系与浅层泥岩的机械压实作用传统认识呈现较大差异。浅层泥岩机械压实作用获得的泥岩正常压实曲线不能外推到中深层泥岩。因此，忽略中深层泥岩压实对泥岩正常压实曲线的影响，必然导致中深层地层压力预测存在严重的错误，给油气勘探带来了巨大的困难。
[0004]
随着研究的深入，已经认识到不同埋藏深度的泥岩具有不同的孔隙度演化趋势，因而泥岩的正常压实曲线不应该被认为是单一的趋势。前人尝试利用蒙脱石和伊利石的压实趋势作为机械压实阶段和化学压实阶段的正常压实曲线。然而实际盆地中泥岩矿物组分并非全是蒙脱石或伊利石，不同矿物组成导致的泥岩压实趋势线差异可以很明显。zhang j c，yin s x.real-time pore pressure detection：indicators and improved methods[j].geofluids，2017，9：1-12通过假设蒙脱石-伊利石转化过程用线性关系来描述，建立了蒙脱石正常压实曲线、蒙-伊过渡阶段正常压实曲线和伊利石正常压实曲线并基于这些正常压实曲线进行压力预测的方法。尽管该方法将泥岩压实分为三个阶段，但其准确性却依赖于盆地内复杂的蒙-伊转化反应的经验认识，并且该方法仅将蒙-伊过渡阶段压实曲线认为是一种简单的线性关系。然而，实际盆地中，泥岩的黏土矿物组成不仅是单一的蒙脱石或伊利石，蒙-伊转化过程也并非简单的线性关系可以解决，不同矿物组分导致泥岩的正常压
实趋势线差异很大，简单的用两个端元组分来预测孔隙压力很可能带来较大的误差。更重要的是，该方法提出伊利石正常压实曲线通过对正常压力段内声波时差-深度数据直接拟合而获得，但是实际上对于伊利石压实仍然保持正常压力的泥岩并不多，进而影响压力预测结果的准确性，并且伊利石正常压实曲线往往很难直接观测获得，因此该方法难以在实际工作中推广应用。
[0005]
因此，迫切需要提出一种新的适用于泥岩不同压实阶段的地层压力预测方法。


技术实现要素：

[0006]
针对现有技术中存在的技术问题，本发明的目的是提供一种用于对钻井泥岩地层压力进行预测的方法，该方法可以更精确地预测泥岩地层压力。
[0007]
为了达到上述目的，本发明提供一种用于对钻井泥岩地层压力进行预测的方法，所述方法包括以下步骤：获取钻井地层的自然伽马测井曲线、声波测井曲线、密度测井曲线和中子测井曲线；利用自然伽马测井曲线获得钻井地层的泥质含量；利用中子测井曲线和密度测井曲线获得钻井地层的中子孔隙度与密度孔隙度之差值；根据钻井地层的泥质含量以及中子孔隙度与密度孔隙度之差值确定泥岩段；利用声波测井曲线和密度测井曲线获得所述泥岩段的声波时差随密度的变化趋势线；根据所述泥岩段的声波时差随密度的变化趋势线将所述泥岩段划分为机械压实阶段、过渡阶段和化学压实阶段；通过对所述机械压实阶段的声波测井曲线进行指数函数拟合，以声波时差与深度关系建立机械压实阶段正常压实曲线；利用所述化学压实阶段的声波测井曲线，通过实测地层压力约束的压实参数反演方法以修正的声波时差与深度关系建立化学压实阶段正常压实曲线；通过分别利用机械压实阶段正常压实曲线和化学压实阶段正常压实曲线获得过渡阶段起点和终点处的声波时差并且基于这两个点的声波时差和地层深度进行计算，以修正的声波时差与深度关系建立过渡阶段正常压实曲线；和利用获得的机械压实阶段正常压实曲线、化学压实阶段正常压实曲线和过渡阶段正常压实曲线对钻井泥岩地层压力进行预测。
[0008]
根据本发明，进一步地，利用获得的机械压实阶段正常压实曲线、化学压实阶段正常压实曲线和过渡阶段正常压实曲线对钻井泥岩地层压力进行预测的步骤包括：将所述机械压实阶段正常压实曲线、所述化学压实阶段正常压实曲线和所述过渡阶段正常压实曲线直接用于对钻井泥岩地层压力进行预测；或者将所述机械压实阶段正常压实曲线、所述化学压实阶段正常压实曲线和所述过渡阶段正常压实曲线分别转换为相应的加载曲线，然后将所述加载曲线用于对钻井泥岩地层压力进行预测。
[0009]
根据本发明，进一步地，在利用自然伽马测井曲线获得钻井地层的泥质含量的步骤中通过以下方程进行计算获得所述钻井地层的泥质含量：
[0010][0011][0012]
式中，i
gr
为泥质含量指数；gr是地层自然伽马值；gr
min
和gr
max
分别是纯砂岩和纯泥岩自然伽马值；v
sh
为泥质含量；并且gcur为经验参数。
[0013]
根据本发明，进一步地，根据所述泥岩段的声波时差随密度的变化趋势线将所述泥岩段划分为机械压实阶段、过渡阶段和化学压实阶段的步骤包括：将其声波时差随密度
的变化趋势线斜率绝对值最大的泥岩段部分确定为过渡阶段；将过渡阶段起点以上的泥岩段部分确定为机械压实阶段；和将过渡阶段终点以下的泥岩段部分确定为化学压实阶段。
[0014]
根据本发明，进一步地，所述机械压实阶段正常压实曲线由以下关系式表示：
[0015][0016]
式中，ac为机械压实阶段泥岩声波时差；z为地层深度；ac
0m
为机械压实阶段泥岩地表声波时差；并且c
m
为机械压实阶段泥岩压实系数。
[0017]
根据本发明，进一步地，所述化学压实阶段正常压实曲线由以下关系式表示：
[0018][0019]
式中，ac为化学压实阶段泥岩声波时差；z为地层深度；z
c
为化学压实阶段起点处的地层深度；ac
c
为化学压实阶段起点处的泥岩声波时差；并且c
c
为化学压实阶段泥岩压实系数。
[0020]
根据本发明，进一步地，所述过渡阶段正常压实曲线由以下关系式表示：
[0021][0022]
式中，ac为过渡阶段泥岩声波时差；z为地层深度；z
m
为过渡阶段起点处的地层深度；ac
t
为过渡阶段起点处的泥岩声波时差；并且c
t
为过渡阶段泥岩压实系数。
[0023]
根据本发明，进一步地，通过以下方程将所述机械压实阶段正常压实曲线、所述化学压实阶段正常压实曲线和所述过渡阶段正常压实曲线分别转换为相应的加载曲线：
[0024][0025]
式中，σ
v
为垂直有效应力；z为地层深度；δz为密度测井采样间隔；ρ
ri
为采样间隔内岩石密度；g为重力加速度；n为采样间隔数目；并且ρ
f
为地层水的密度。
[0026]
根据本发明，进一步地，所述实测地层压力约束的压实参数反演方法包括：
[0027]
首先对c
c
赋予初值，基于平衡深度法由以下方程估算地层压力p：
[0028][0029]
式中，p为地层压力，z为地层深度，ρ
r
为岩石体积密度，ρ
w
为孔隙水密度，g为重力加速度，c
c
为化学压实阶段泥岩压实系数，ac
0c
为化学压实阶段外推地表声波时差，并且ac为化学压实阶段泥岩声波时差；
[0030]
将估算的地层压力p与实测地层压力进行比较；和
[0031]
根据估算的地层压力p与实测地层压力是否近似来确定是否重新对c
c
赋予初值并且估算地层压力p直至估算的地层压力p与实测地层压力之间的误差最小，以获得该初值作为c
c
。
[0032]
根据本发明，进一步地，从所述自然伽马测井曲线、声波测井曲线、密度测井曲线和中子测井曲线中剔除受到井筒扩径或坍塌影响的中浅层测井数据，以及井径大于钻头直径20％以上的扩径层段测井数据。
[0033]
有益效果
[0034]
根据本发明的方法通过建立更符合泥岩实际压实的不同压实阶段正常压实曲线，可以更精确地预测泥岩地层压力。此外，本方法还可以进一步利用所得到的不同压实阶段正常压实曲线获得相应的加载曲线并将其用于超压成因分析，从而提高超压成因分析的准
确性。本方法还有助于指导泥岩物性预测、地球物理解释、构造演化恢复。
附图说明
[0035]
图1为说明根据本发明的实测地层压力约束的压实参数反演方法的示意图。
[0036]
图2显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井地层的自然伽马测井曲线(gr)、声波测井曲线(ac)、密度测井曲线(den)、中子测井曲线(cnl)、井径测井曲线(cal)、钻头直径测井曲线(bit)、泥质含量(v
sh
)-深度曲线和中子孔隙度与密度孔隙度之差-深度曲线。
[0037]
图3显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井地层的泥质含量v
sh
与中子孔隙度-密度孔隙度之差值交会图。
[0038]
图4显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井泥岩段的声波时差-深度图和密度-深度图。
[0039]
图5显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井泥岩段的声波时差-密度交会图。
[0040]
图6显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井泥岩不同压实阶段正常压实曲线以及根据本发明实施例的方法预测的地层压力和基于单一机械压实曲线预测的地层压力与实测地层压力的比较。
[0041]
图7显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井泥岩的利用不同正常压实曲线的地层压力预测的误差分析。
[0042]
图8显示了根据本发明的一个实施例的渤海湾盆地东营凹陷ly1井泥岩不同压实阶段加载曲线及超压成因判识。
具体实施方式
[0043]
根据本发明的方法包括以下步骤：获取钻井地层的自然伽马测井曲线、声波测井曲线、密度测井曲线和中子测井曲线；利用自然伽马测井曲线获得测井地层的泥质含量；利用中子测井曲线和密度测井曲线获得钻井地层的中子孔隙度与密度孔隙度之差值；根据钻井地层的泥质含量以及中子孔隙度与密度孔隙度之差值确定泥岩段；利用声波测井曲线和密度测井曲线获得所述泥岩段的声波时差随密度的变化趋势线；根据所述泥岩段的声波时差随密度的变化趋势线将所述泥岩段划分为机械压实阶段、过渡阶段和化学压实阶段；通过对所述机械压实阶段的声波测井曲线进行指数函数拟合，以声波时差与深度关系建立机械压实阶段正常压实曲线；利用所述化学压实阶段的声波测井曲线，通过实测地层压力约束的压实参数反演方法以修正的声波时差与深度关系建立化学压实阶段正常压实曲线；通过分别利用机械压实阶段正常压实曲线和化学压实阶段正常压实曲线获得过渡阶段起点和终点处的声波时差并且基于这两个点的声波时差和地层深度进行计算，以修正的声波时差与深度关系建立过渡阶段正常压实曲线；和利用获得的机械压实阶段正常压实曲线、化学压实阶段正常压实曲线和过渡阶段正常压实曲线对钻井泥岩地层压力进行预测。
[0044]
在本方法中，自然伽马测井曲线是指钻井地层的自然伽马值随深度变化的曲线，通过自然伽马测井而获得；声波测井曲线是指钻井地层的声波时差随深度变化的曲线，通过声波测井而获得；密度测井曲线是指钻井地层的体积密度随深度变化的曲线，通过密度
测井而获得；而中子测井曲线是指钻井地层的中子值随深度变化的曲线，通过中子测井而获得。自然伽马测井、声波测井、密度测井和中子测井方法在本领域中是已知的。在本方法中，还包括获取井径(cal)及钻头直径(bit)等测井曲线。在获得这些测井曲线之后，为了保证测井数据的可靠性，可以对这些测井曲线进行处理。首先剔除受到井筒扩径或坍塌影响的中浅层测井数据，其次，通过对比井径与钻头直径，剔除井径大于钻头直径20％以上的扩径层段测井数值，还可以剔除受到断层、裂缝、有机质含量等因素导致的异常测井数值。
[0045]
在利用自然伽马测井曲线获得钻井地层的泥质含量的步骤中，通过以下方程进行计算获得所述钻井地层的泥质含量：
[0046][0047][0048]
式中，i
gr
为泥质含量指数；gr是地层自然伽马值；gr
min
和gr
max
分别是纯砂岩和纯泥岩自然伽马值；v
sh
为泥质含量；并且gcur为经验参数。在本领域中，gcur通常如下确定：新地层(古近系以后)取值为2；老地层(古近系以前)取值为3.7。
[0049]
在本方法中，可以利用中子测井曲线和密度测井曲线获得钻井地层的中子孔隙度与密度孔隙度之差值。具体地，可以利用中子测井曲线和密度测井曲线分别计算中子孔隙度和密度孔隙度
[0050][0051][0052]
式(3)中，为中子孔隙度，％；φ
ma
为骨架中子值，取值为-0.05；φ
f
为孔隙流体中子值，取值为1.00；φ为地层中子值，％。式(4)中，为密度孔隙度，％；ρ
ma
为骨架密度，取值为2.7g/cm3；ρ
f
为孔隙流体密度，取值为1.0g/cm3；ρ为地层体积密度，g/cm3。
[0053]
在如上获得中子孔隙度与密度孔隙度之后，计算两者之差值即可得到中子孔隙度与密度孔隙度之差值然后，可以根据钻井地层的泥质含量v
sh
以及中子孔隙度与密度孔隙度之差值确定泥岩段。具体地，可以绘制v
sh
与中子孔隙度-密度孔隙度差值交会图。在一个或多个实施例中，可以根据该交会图选择v
sh
＞75％及作为纯泥岩的阈值，从而确定泥岩段。
[0054]
在确定泥岩段之后，可以利用声波测井曲线和密度测井曲线获得所述泥岩段的声波时差随密度的变化趋势线。可以选取厚度大于5～10m的纯泥岩段，在此基础上统计泥岩段的声波时差和体积密度测井的平均值。可以使用泥岩段的声波时差和体积密度测井的平均值来绘制泥岩声波时差-密度交会图。泥岩声波时差-密度交会图可以非常直观地反映声波时差随密度的变化关系。可以根据泥岩段的声波时差随密度的变化趋势线将所述泥岩段划分为机械压实阶段、过渡阶段和化学压实阶段。根据本发明，泥岩在浅层、低温条件下压实过程表现为声波时差减小、密度增大的特征，遵循机械压实趋势；而后泥岩压实趋势从机械压实曲线逐渐过渡到化学压实曲线，在此过程中黏土矿物快速转化导致泥岩密度快速增大现象，对应于泥岩压实的过渡阶段，随后沿着化学压实趋势继续发生压实作用。因此，可以依据声波时差随密度的变化趋势线的斜率变化来确定泥岩部分所属压实阶段，具体地，
可以首先将其声波时差随密度的变化趋势线斜率绝对值最大的泥岩段部分确定为过渡阶段，然后将过渡阶段起点以上的泥岩段部分确定为机械压实阶段，并且将过渡阶段终点以下的泥岩段部分确定为化学压实阶段。在一个或多个实施例中，泥岩的压实过程可以为划分如下三个阶段：机械压实阶段(＜2000m)、过渡阶段(2000～2500m)和化学压实阶段(＞2500m)。
[0055]
根据本发明，处于过渡阶段起点(z
m
，亦即，机械压实阶段终点)以上的泥岩压实过程定义为泥岩的机械压实阶段。机械压实阶段正常压实曲线可以表示为声波时差(ac)-深度(z)的指数函数：
[0056][0057]
式中，ac为机械压实阶段泥岩声波时差，μs/ft；z为地层深度，m；ac
0m
为机械压实阶段泥岩地表声波时差，μs/ft；并且c
m
为机械压实阶段泥岩压实系数。
[0058]
因此，可以通过对所述机械压实阶段的声波测井曲线进行指数函数拟合，以声波时差与深度关系建立机械压实阶段正常压实曲线。
[0059]
根据本发明，处于过渡阶段终点(z
c
，亦即，化学压实阶段起点)以下的泥岩压实过程定义为泥岩化学压实阶段。化学压实阶段泥岩正常压实曲线用修正的声波时差(ac)-深度(z)关系表述：
[0060][0061]
式中，ac为化学压实阶段泥岩声波时差；z为地层深度；z
c
为化学压实阶段起点处的地层深度；ac
c
为化学压实阶段起点处的泥岩声波时差；并且c
c
为化学压实阶段泥岩压实系数。
[0062]
基于本发明，化学压实阶段泥岩往往发育超压致使声波时差-深度关系偏离正常压实曲线，如图1所示。基于本发明，为了建立符合实际压实的化学压实阶段正常压实曲线，特别采用了实测地层压力约束的压实参数反演方法。根据不同矿物组分的机械压实实验，不同矿物组分在未压实状态下原始孔隙度近似相等，因而化学压实与机械压实泥岩应该具有相同的地表声波时差，即：ac
0c
＝ac
0m
，因此可以将方程(6)转化为：
[0063][0064]
式中，ac为化学压实阶段泥岩声波时差；z为地层深度；ac
0c
为化学压实阶段外推地表声波时差；并且c
c
为化学压实阶段泥岩压实系数。
[0065]
在一个或多个实施方案中，实测地层压力约束的压实参数反演方法包括：
[0066]
首先对c
c
赋予初值c
ci
(c
ci
≥c
m
)，基于平衡深度法由以下方程估算地层压力p：
[0067][0068]
式中，p为地层压力，mpa；z为地层深度，m；ρ
r
为岩石体积密度，g/cm3；ρ
w
为孔隙流体密度，1.0g/cm3；g为重力加速度；c
c
为化学压实阶段泥岩压实系数；ac
0c
为化学压实阶段外推地表声波时差，μs/ft；并且ac为化学压实阶段泥岩声波时差，μs/ft；
[0069]
将估算的地层压力p与实测地层压力进行比较；和
[0070]
根据估算的地层压力p与实测地层压力是否近似来确定是否重新对c
c
赋予初值并且估算地层压力p直至估算的地层压力p与实测地层压力之间的误差最小，以获得该初值作为c
c
。
[0071]
根据本发明，如果估算的地层压力p与实测地层压力近似，则将此时的初值作为c
c
；如果估算的地层压力p与实测地层压力不近似，则重新对c
c
赋予初值并且估算地层压力p直至估算的地层压力p与实测地层压力之间的误差最小，以获得该初值作为c
c
根据本发明，估算的地层压力p与实测地层压力近似的情况包括估算的地层压力p与实测地层压力相等的情况以及估算的地层压力p与实测地层压力之间的误差最小的情况。具体地，使用实测地层压力(或泥浆压力)作为约束，如果预测地层压力等于实测地层压力(或泥浆压力)，则将该c
ci
值作为c
c
。如果预测地层压力小于实测地层压力(或泥浆压力)，则以一定的梯度(δc
i
)不断增大c
c
值，重新预测地层压力，对比预测地层压力与实测地层压力(或泥浆压力)，当两者误差达到最小时，此时对应的cc值即可用于建立化学压实阶段正常压实曲线。根据本发明，实测地层压力(或泥浆压力)可以利用钻井过程中dst、mdt或rft测试(或泥浆密度)获得。
[0072]
根据本发明，介于z
m
～z
c
深度区间的泥岩处于泥岩压实的过渡阶段。过渡阶段内泥岩的正常压实曲线用修正的声波时差(ac)-深度(z)关系表述：
[0073][0074]
式中，ac为过渡阶段泥岩声波时差，μs/ft；z为地层深度，m；z
m
为过渡阶段起点处的地层深度，m；ac
t
为过渡阶段起点处的泥岩声波时差，μs/ff；并且c
t
为过渡阶段泥岩压实系数。
[0075]
过渡阶段起点处声波时差ac
t
等于机械压实阶段终点处(z
m
)的声波时差(ac
m
)，可以从机械压实曲线计算获得，同时过渡阶段终点处声波时差等于化学压实阶段起点处(z
c
)的声波时差(ac
c
)，可以通过从化学压实曲线计算获得。因此，利用机械压实曲线和化学压实曲线进行计算，可以获得两点(z
m
，ac
m
)和(z
c
，ac
c
)。基于该两点进行计算，可以获得过渡阶段泥岩压实系数，从而建立过渡阶段正常压实曲线。
[0076]
因此，以上述方式获得了耦合泥岩不同压实阶段的综合正常压实曲线：
[0077][0078]
同一地区泥岩压实作用具有区域相似性，因而在建立的综合正常压实曲线基础上，如下利用平衡深度法可以实现不同压实阶段的地层压力预测：
[0079][0080]
式中，p为地层压力，mpa；z为地层深度，m；z
i
为不同压实阶段起点处的地层深度，m；ρ
r
为岩石体积密度，g/cm3；ρ
w
为孔隙流体密度，1.0g/cm3；g为重力加速度；c
i
为不同压实阶段泥岩压实系数；ac
i
为不同压实阶段起点处的声波时差，μs/ft；并且ac为不同压实阶段泥岩声波时差，μs/ft；其中，对于机械压实阶段(z＞z
m
)：z
i
＝0，c
i
＝c
m
，ac
i
＝ac
0m
；对于过渡阶段(z
m
＜z＜z
c
)：z
i
＝z
m
，c
i
＝c
t
，ac
i
＝ac
t
；对于化学压实阶段(z＞z
c
)：z
i
＝z
c
，c
i
＝c
c
，ac
i
＝ac
c
。
[0081]
在其他一些实施方案中，可以将所述机械压实阶段正常压实曲线、所述化学压实阶段正常压实曲线和所述过渡阶段正常压实曲线分别转换为相应的加载曲线，即，将不同压实阶段的声波时差-深度关系转化为声波时差-垂直有效应力关系。然后，将所述加载曲
线用于对钻井泥岩地层压力进行预测。加载曲线除了可以用于孔隙压力预测之外，还可以用于超压成因分析。
[0082]
以下对将不同压实阶段正常压实曲线转换为相应的加载曲线进行说明。
[0083]
根据有效应力基本原理，正常压力条件下地层承受的垂直有效应力等于相同深度静岩压力(s
v
)与静水压力(p
h
)的差值：
[0084]
σ
v
＝s
v-p
h
ꢀꢀꢀ
(12)
[0085]
式中，σ
v
为垂直有效应力，mpa；s
v
为静岩压力，mpa；p
h
为静水压力，mpa。
[0086]
静岩压力为上覆岩石施加的压力，特定深度的静岩压力可以根据上覆岩石的体积密度计算获得：
[0087]
s
v
＝ρ
r
·
g
·
z
ꢀꢀꢀ
(13)
[0088]
式中，s
v
为静岩压力，mpa；z为深度，m；ρ
r
为地层体积密度，g/cm3；g为重力加速度。
[0089]
通常密度测井以特定的深度间隔记录井筒周围的地层密度，对于每一个采样间隔内，使用以下公式计算每一测井间隔深度内静岩应力(sv
i
)：
[0090]
s
vi
＝ρ
ri
·
g
·
δz
ꢀꢀꢀ
(14)
[0091]
式中，sv
i
为每一测井间隔深度内静岩应力，mpa；δz为密度测井采样间隔，m；为采样深度区间内岩石密度，g/cm3；g为重力加速度。
[0092]
为了获得深度z处的垂直应力，需要将每个单独的垂直应力增量求和：
[0093][0094]
式中，s
v
为静岩压力，mpa；n为采样间隔数目；δz为密度测井采样间隔，m；为采样深度区间内岩石密度，g/cm3；g为重力加速度。
[0095]
特定深度的静水压力(p
h
)为该深度以上连通水柱的压力，可以表示为：
[0096]
p
h
＝ρ
f
·
g
·
z
ꢀꢀꢀ
(16)
[0097]
式中，p
h
为静水压力，mpa；为孔隙水密度，g/cm3；g为重力加速度；z为深度，m。
[0098]
进而，特定深度处的垂直有效应力为：
[0099][0100]
式中，σ
v
为垂直有效应力，mpa；n为采样间隔数目；为采样深度区间内岩石密度，g/cm3；g为重力加速度；δz为密度测井采样间隔，m；为孔隙水密度，g/cm3；g为重力加速度；z为深度，m。
[0101]
根据上述公式，可以获得不同压实阶段正常压实曲线中相应深度处的垂直有效应力，进而将声波时差-深度关系转化为声波时差-有效应力关系，从而获得机械压实阶段、过渡阶段和化学压实阶段的加载曲线。
[0102]
为了使本发明的目的、技术方案及优点更加清楚明白，以下将通过实施例来具体说明本发明。应当懂得，这些实施例仅用于说明的目的，它们并不以任何方式对本发明构成限制。
[0103]
实施例
[0104]
以渤海湾盆地东营凹陷ly1井为例。通过自然伽马测井、声波测井、密度测井和中子测井方法获取钻井地层的自然伽马测井曲线、声波测井曲线、密度测井曲线和中子测井
曲线。图2显示了渤海湾盆地东营凹陷ly1井的自然伽马测井曲线(gr)、声波测井曲线(ac)、密度测井曲线(den)和中子测井曲线(cnl)。
[0105]
对上述测井曲线进行处理。首先，剔除中浅部由于地层压实和成岩胶结程度低导致井筒扩径或坍塌引起的波动较大的测井数值。进一步，通过对比井径(cal)与钻头直径(bit)测井曲线(如图2中的cal/bit)，确定了井径大于钻头直径20％以上的扩径层段，这些测井数据也被排除在外。
[0106]
利用图2所示的自然伽马测井曲线(gr)，通过公式(1)和(2)计算获得ly1井地层的泥质含量v
sh
(参见图2)。利用图2所示的中子测井曲线(cnl)和密度测井曲线(den)，通过公式(3)和(4)计算获得钻井地层的中子孔隙度与密度孔隙度之差值(参见图2)。绘制v
sh
与中子孔隙度-密度孔隙度之差值交会图(参见图3)。选择v
sh
＞75％及作为纯泥岩的阈值，确定泥岩段。
[0107]
为降低围岩对目标泥岩段测井响应的影响，选取厚度一般大于5～10m的纯泥岩段，统计该泥岩段的声波测井响应和密度测井响应的平均值，绘制声波时差-深度图和密度-深度图(见图4)。
[0108]
利用图4所示的声波时差-深度图和密度-深度图编制泥岩段声波时差-密度交会图(见图5)。根据声波时差-密度交会图，可以确定埋深小于2000m(具体为1500-2000m)的泥岩段处于机械压实阶段；埋深介于2000-2500m泥岩段密度表现出快速增大的趋势，处于过渡阶段；而埋深超过2500m(具体为2500-4000m)的泥岩段处于化学压实阶段。
[0109]
对埋深小于2000m的机械压实阶段浅泥岩声波时差-深度曲线进行指数函数拟合，从而确定机械压实阶段正常压实曲线为：
[0110]
ac＝189.2098
·
e-0.0001985
·
z
ꢀꢀꢀ
(18)
[0111]
以实测地层压力作为约束，对cc赋予初值c
ci
＝0.0002，利用方程(8)估算的地层压力小于实测地层压力，以一定的梯度(δc
i
＝0.0000002)不断调整压实系数c
c
，利用方程(8)估算不同c
c
条件下的地层压力，对比估算地层压力与实测地层压力，当估算地层压力与实测地层压力误差最小时，确定此时的c
ci
＝0.0002036作为化学压实阶段压实系数c
c
。因此泥岩化学压实阶段正常压实曲线为：
[0112]
ac＝189.2098
·
e-0.00002036
·
z
ꢀꢀꢀ
(19)
[0113]
2000-2500m之间为泥岩压实的过渡阶段，此时利用机械压实阶段正常压实曲线(18)和化学压实阶段正常压实曲线(19)计算出过渡阶段起点为z
m
＝2000m，ac
m
＝127.1701μs/ft；终点为z
c
＝2500m，ac
c
＝104.9658μs/ft。基于这两个点(z
m
，ac
m
)和(z
c
，ac
c
)进行计算，获得ac
t
＝273.9894μs/ff，c
t
＝0.00003838。所以可以得到过渡阶段正常压实曲线为：
[0114]
ac＝273.9894
·
e-0.00003838
·
z
ꢀꢀꢀ
(20)
[0115]
进而，获得泥岩压实的综合正常压实曲线：
[0116][0117]
所建立的综合正常压实曲线如图6左图所示。以建立的综合正常压实曲线为依据，利用方程(11)所示的平衡深度法可以预测不同压实阶段地层压力，预测地层压力与实测地层压力达到良好的匹配关系，而基于单一机械压实曲线预测的地层压力与实测地层压力偏
差较大(见图6右图)。
[0118]
利用背景技术中描述的现有技术(即，利用单一机械压实正常压实曲线对钻井泥岩地层压力进行预测的方法以及zhang j c，yin s x.real-time pore pressure detection：indicators and improved methods[j].geofluids，2017，9：1-12所述的预测方法以及本发明的方法分别对渤海湾盆地东营凹陷ly1井进行压力预测，统计不同深度点实测地层压力处对应的地层压力预测值，如图7所示。对比实测地层压力与预测地层压力可以发现，通过利用单一机械压实正常压实曲线对钻井泥岩地层压力进行预测的方法预测的地层压力普遍小于实测地层压力，绝对误差超过10mpa，而利用zhang j c，yin s x.real-time pore pressure detection：indicators and improved methods[j].geofluids，2017，9：1-12所述的预测方法预测的压力则大于实测地层压力，绝对误差也可达10mpa。而通过本发明方法预测的地层压力与实测地层压力具有远远更好的吻合关系，几乎保持一致，充分证明本发明提出的方法能够更精确地预测地层压力。
[0119]
接着，根据图2所示的密度测井曲线，获得岩石体积密度随深度变化规律，通过方程(15)对体积密度进行积分获得静岩压力与深度关系：
[0120]
s
v
＝0.0099
·
z
1.0962
ꢀꢀꢀ
(22)
[0121]
通过方程(17)以静岩压力与相同深度静水压力之差获得特定深度处垂直有效应力：
[0122]
σ
v
＝0.0028
·
z
1.1792
ꢀꢀꢀ
(23)
[0123]
将泥岩不同压实阶段正常压实曲线(21)中深度计算为垂直有效应力，进而获得如下所示的泥岩机械压实阶段、过渡阶段及化学压实阶段的加载曲线(参见图8)：
[0124][0125]
同一地区泥岩压实作用具有区域相似性，在计算获得泥岩不同压实阶段对应的加载曲线的基础上，可以应用于地层压力准确预测，也可以用于超压成因判识。
[0126]
在声波时差-垂直有效应力交会图(图8)中，压力系数介于1.0-1.2之间的常压数据点基本落在过渡阶段加载曲线和化学压实阶段加载曲线上；压力系数介于1.2-1.4之间的超压数据点落在化学压实阶段加载曲线上，表明超压成因主要为泥岩化学压实作用；而压力系数大于1.4时的超压数据点则落在化学压实阶段加载曲线之外，反映了卸载成因机制的贡献。
[0127]
上述实施例仅例示性的说明了本发明，而非用于限制本发明。熟知本领域的技术人员应当理解，在不偏离本发明的精神和范围的情况下，对本发明实施例所作的任何更改和变化均落在本发明的范围内。且本发明的保护范围应由所附的权利要求确定。