一种基于双重校正的泥页岩地层孔隙压力的计算方法与流程

1.本发明属于石油钻井技术领域，具体涉及一种基于双重校正的泥页岩地层孔隙压力计算方法。


背景技术：

2.储层地质力学指在储层开发中遇到的各种地质力学问题，一般涉及岩石力学、孔隙压力、地应力等参数，其中孔隙压力对于井壁稳定性、压裂改造等有重要的影响。在储层孔隙压力预测方面，主要可分为三类，基于欠压实理论的等效深度法和eaton法、基于有效应力理论的bowers和有效应力法、综合解释法，这些预测方法在地层孔隙压力预测方面都有一定的适应性。随着能源需求的增加，非常规页岩油气储层受到越来越多的重视。现有生产表明页岩油气储层中存在异常高压地层，其形成机制除了与生烃作用有关外，还与页岩储层岩石中富含有机质有关。随着岩石中有机质含量增大，岩石的声波时差将增大，基于压实理论获取的正常压实趋势曲线时，实测地层声波时差将偏离正常压实趋势曲线，造成地层压力预测值偏高的现象。另外，在页岩油气储层钻井过程中，钻井液侵入地层将发生较明显的水化作用，水化作用也将造成岩石声波时差增大，而裸眼测井获得地层的声波数据，都是水化后地层的声波数据。
3.因此，页岩油气储层钻井过程中，会受到有机质和水化作用的双重影响，常规的孔隙度计算方法无法直接应用。


技术实现要素：

4.本发明的目的是提供一种基于双重校正的泥页岩地层孔隙压力计算方法，通过对水化作用和有机质的先后校正，计算的地层孔隙压力大小更加接近原始地层状态，提高了预测精度，相比于常规方法更符合泥页岩地质特征，算出的孔隙压力更为合理。
5.本发明采取的技术方案是：
6.一种基于双重校正的泥页岩地层孔隙压力计算方法，包括以下步骤：
7.步骤一：准备研究工区的实验岩心资料、地质资料、现场钻井用钻井液、测井资料、钻完井报告；
8.步骤二：对获取的研究工区井下实验样品进行岩心描述，且钻取标准岩样数60
‑
80块，对岩样烘干后，测量干燥岩样的密度、电阻率和超声波波速，获取原岩岩样的密度、电阻率和声波时差(分别表示为den0、r0、δt0，指未发生浸泡水化的值)；
9.步骤三：采用现场钻井用钻井液对步骤二的干燥岩样进行浸泡实验(浸泡过程中岩样与钻井液将发生水化作用进而改变岩石的物理性质)，对浸泡后的岩样进行密度、电阻率、声波测量，获取浸泡后岩样的密度、电阻率和声波时差(分别表示为den、r、δt，发现浸泡后岩样密度增大、电阻率减小、声波时差增大)；
10.结合步骤二，获得密度、电阻率和声波时差的改变量(δden、δr、δ(δt)1)；以电阻率改变量为自变量，分别以密度改变量、声波时差改变量为因变量进行拟合回归(分析发
现密度改变量、声波时差改变量与电阻率改变量存在较好的关系)，见式(1)和式(2)所示，
11.另外，电阻率改变量又能通过双侧向电阻率测井资料获得，见式(3)所示(钻井过程中，钻井液会侵入地层，在井壁附近形成侵入带，侵入带地层的岩石物理性质会发生改变，现场的双侧向电阻率测井可反映出钻井液对侵入带地层和原状地层电阻率的差异，其中深侧向电阻率r
d
反映原状地层电阻率(即r0＝r
d
)，浅侧向电阻率r
s
反映侵入带地层电阻率，进而可获得钻井液侵入后地层电阻率改变量。通过双侧向电阻率测井值，可将本技术的室内实验结果与现场实际测量结果建立关系)，
12.δ(δt)1＝f1(δr)＝aδr b
ꢀꢀꢀ
(1)
13.δden＝f2(δr)＝dln(δr)
‑
e
ꢀꢀꢀ
(2)
14.δr＝r
d
‑
r
s
ꢀꢀꢀ
(3)
15.式中：δ(δt)1、δr、δden分别为岩样因为浸泡水化而导致的声波时差改变量、电阻率改变量、密度改变量；r
d
为深侧向电阻率值；r
s
为浅侧向电阻率值；a、b、d、e为待定系数；
16.步骤四：基于现场测井资料中的声波时差、密度、深侧向电阻率、浅侧向电阻率资料，获得经水化校正后的原状地层声波时差、密度值为：
17.δt0＝δt
‑
δ(δt)1ꢀꢀꢀ
(4)
18.den0＝den
‑
δden
ꢀꢀꢀ
(5)
19.式中：δt为测井资料中的地层声波时差(水化后的)，us/m；δt0为水化校正后的地层声波时差(原状地层的)，us/m；den为测井资料中的地层密度，g/cm3；den0为水化校正后的地层密度(原状地层的)，g/cm3；
20.步骤五：通过实验同步对步骤二的干燥岩样进行有机碳含量测试，获得每个岩样的有机碳含量，从而得到不同含量有机碳岩样的声波时差与有机碳含量间的关系式，
21.δt(toc)＝fln(ω(toc)) p
ꢀꢀꢀ
(6)
22.式中：δt(toc)为不同含量有机碳的岩样的声波时差，us/m；ω(toc)为岩样的有机碳含量，％；f、p为待定系数；
23.对式(6)求取当岩样有机碳含量趋于0时极限，获取岩样声波时差的极限值δt
lim
：
[0024][0025]
步骤六：以岩样声波时差极限值作为基准值，获取岩样声波时差与该极限值的差值，定义该差值为岩样声波时差的校正量，并获取该校正量与有机碳含量间的关系式，见式(7)；
[0026]
δ(δt)2＝δt(toc)
‑
δt
lim
＝fln(ω(toc)) p
‑
δt
lim
ꢀꢀꢀ
(7)
[0027]
式中：δ(δt)2为(原状地层)基于有机碳含量的岩样声波时差校正量；
[0028]
步骤七：基于水化校正后的地层声波时差，结合测井资料中的有机碳含量信息，对富含有机质地层声波时差进行有机质校正，最终，经过水化校正和有机质校正后的地层声波时差为
[0029]
δt
(水化 toc)校正
＝δt0‑
δ(δt)2ꢀꢀꢀ
(8)
[0030]
式中：δt
(水化 toc)校正
为经过双重校正(水化校正和有机碳校正)后的真实声波时差，us/m；
[0031]
步骤八：利用测井资料中自然伽马信息，计算地层的泥质含量，利用测井资料中井径信息和钻完井报告中钻头尺寸信息，计算地层的扩径率，进而，选择泥质含量大于85％、井径扩径率小于10％泥岩段地层，基于压实理论，构建研究区块的基于校正后声波时差的压实趋势曲线，其方程如式(9)，
[0032]
h＝hln(δt
(水化 toc)校正
) j
ꢀꢀꢀ
(9)
[0033]
式中：h为地层深度，m；h、j为待定系数；
[0034]
步骤九：结合钻完井报告的中途测试或静压测试获取的地层压力实测资料，利用式(10)计算求出不同地层压力测试点下的压实校正系数c，发现压实校正系数c值不是固定值，是一个与深度有关的函数，以此建立了压实校正系数c与深度间的关系，见式(11)，
[0035][0036]
式中：p
p
为地层压力，mpa；σ
v
为上覆岩层压力，mpa；p
h
为地层静水压力，mpa；δt
n
为同一深度正常压实趋势线上的声波时差，us/m；c为压实校正系数，
[0037]
c＝mh
n
ꢀꢀꢀ
(11)
[0038]
式中：h为地层深度，m；m、n为待定系数；
[0039]
其中，
[0040][0041]
式中：h0为起始测井点深度，m；为未测井段深度h0点的密度，g/cm3；g为重力加速度，m/s2；den
水化前
为水化校正后的密度测井值，g/cm3。
[0042]
本发明的有益效果：
[0043]
本发明首先通过室内岩样进行密度、电阻率和声波时差水化校正，再分别建立密度和声波时差与电阻率的拟合关系；之后对原状地层岩心进行有机质校正，结合上一步的水化校正，实现水化和有机质的双重校正；在此基础上，建立地层孔隙压力与压实系数、上覆岩层压力、地层静水压力等参数的函数关系，最终得到进行了两种类型校正的孔隙压力计算方法。
[0044]
在有机质(有机碳)校正中，通过求极限的方法，既能准确获得有机质含量对测井资料影响的阈值，又能解决传统方法以某个较小值近似代替产生的精度差的问题。
[0045]
经过上述校正，计算的地层孔隙压力大小更加接近原始地层状态，提高了预测精度，这是因为消除了因为钻井过程中钻井液浸泡水化而对测井声波时差、密度的影响，考虑了时间因素的基础上修正了测井曲线值；除此之外，进行测井资料的有机质含量校正，相比于常规方法更符合泥页岩地质特征，算出的孔隙压力更为合理。
附图说明
[0046]
图1岩样声波时差与极限值的校正量与toc含量间的关系；
[0047]
图2研究区地层校正后压实趋势线；
[0048]
图3压实校正系数与深度间的关系。
具体实施方式
[0049]
实施例
[0050]
一种基于双重校正的泥页岩地层孔隙压力计算方法，分别对地层声波时差水化校正、地层声波时差的有机质含量校正，再进行地层压实趋势方程的构建，得到不同深度点的压实校正系数，从而构建压实校正系数的计算方法，得到地层孔隙压力的计算方程；
[0051]
具体包括以下步骤：
[0052]
步骤一：准备研究工区的实验岩心资料、地质资料、现场钻井用钻井液、测井资料、钻完井报告；
[0053]
步骤二：对获取的研究工区井下实验样品进行岩心描述，且钻取标准岩样数60块，对岩样烘干后，测量干燥岩样的密度、电阻率和超声波波速，获取原岩岩样的密度、电阻率和声波时差(分别表示为den0、r0、δt0，指未发生浸泡水化的值)；
[0054]
步骤三：采用现场钻井用钻井液对步骤二的干燥岩样进行浸泡实验(浸泡过程中岩样与钻井液将发生水化作用进而改变岩石的物理性质)，对浸泡后的岩样进行密度、电阻率、声波测量，获取浸泡后岩样的密度、电阻率和声波时差(分别表示为den、r、δt，发现浸泡后岩样密度增大、电阻率减小、声波时差增大)；
[0055]
结合步骤二，获得密度、电阻率和声波时差的改变量(δden、δr、δ(δt)1)；以电阻率改变量为自变量，分别以密度改变量、声波时差改变量为因变量进行拟合回归(分析发现密度改变量、声波时差改变量与电阻率改变量存在较好的关系)，见式(1)和式(2)所示，
[0056]
另外，电阻率改变量又能通过双侧向电阻率测井资料获得，见式(3)所示(钻井过程中，钻井液会侵入地层，在井壁附近形成侵入带，侵入带地层的岩石物理性质会发生改变，现场的双侧向电阻率测井可反映出钻井液对侵入带地层和原状地层电阻率的差异，其中深侧向电阻率r
d
反映原状地层电阻率(即r0＝r
d
)，浅侧向电阻率r
s
反映侵入带地层电阻率，进而可获得钻井液侵入后地层电阻率改变量。通过双侧向电阻率测井值，可将本技术的室内实验结果与现场实际测量结果建立关系)，
[0057]
δ(δt)1＝f1(δr)＝1.5896δr
‑
1.0593
ꢀꢀꢀ
(1)
[0058]
δden＝f2(δr)＝0.0663ln(δr)
‑
0.1692
ꢀꢀꢀ
(2)
[0059]
δr＝r
d
‑
r
s
ꢀꢀꢀ
(3)
[0060]
式中：δ(δt)1、δr、δden分别为岩样因为浸泡水化而导致的声波时差改变量、电阻率改变量、密度改变量；r
d
为深侧向电阻率值；r
s
为浅侧向电阻率值；
[0061]
步骤四：基于现场测井资料中的声波时差、密度、深侧向电阻率、浅侧向电阻率资料，获得经水化校正后的原状地层声波时差、密度值为(将步骤三的公式(1)～(3)代入公式(4)和(5)即可)：
[0062]
δt0＝δt
‑
δ(δt)1ꢀꢀꢀ
(4)
[0063]
den0＝den
‑
δden
ꢀꢀꢀ
(5)
[0064]
式中：δt为测井资料中的地层声波时差(水化后的)，us/m；δt0为水化校正后的地层声波时差(原状地层的)，us/m；den为测井资料中的地层密度，g/cm3；den0为水化校正后的地层密度(原状地层的)，g/cm3；
[0065]
步骤五：通过实验同步对步骤二的干燥岩样进行有机碳含量测试，获得每个岩样的有机碳含量，从而得到不同含量有机碳岩样的声波时差与有机碳含量间的关系式，
[0066]
δt(toc)＝58.951ln(ω(toc)) 233.98
ꢀꢀꢀ
(6)
[0067]
式中：δt(toc)为不同含量有机碳的岩样的声波时差，us/m；ω(toc)为岩样的有机碳含量，％；
[0068]
对式(6)求取当岩样有机碳含量趋于0时极限，获取岩样声波时差的极限值δt
lim
：
[0069][0070]
步骤六：以岩样声波时差极限值作为基准值，获取岩样声波时差与该极限值的差值，定义该差值为岩样声波时差的校正量，并获取该校正量与有机碳含量间的关系式，见式(7)，如图1所示；
[0071][0072]
式中：δ(δt)2为(原状地层)基于有机碳含量的岩样声波时差校正量；
[0073]
步骤七：基于水化校正后的地层声波时差，结合测井资料中的有机碳含量信息，对富含有机质地层声波时差进行有机质校正，有机质校正后的地层声波时差为
[0074]
δt
(水化 toc)校正
＝δt0‑
δ(δt)2ꢀꢀꢀ
(8)
[0075]
式中：δt
(水化 toc)校正
为经过双重校正(水化校正和有机碳校正)后的真实声波时差，us/m；需要注意的是当地层中不含有有机质时，可不必进行地层声波时差的有机质校正，其声波时差为水化校正后的地层声波时差。
[0076]
步骤八：利用测井资料中自然伽马信息，计算地层的泥质含量，利用测井资料中井径信息和钻完井报告中钻头尺寸信息，计算地层的扩径率，进而，选择泥质含量大于85％、井径扩径率小于10％泥岩段地层，基于压实理论，构建研究区块的基于校正后声波时差的压实趋势曲线，如图2所示，其方程如式(9)，
[0077]
h＝
‑
6430.3ln(δt
(水化 toc)校正
) 31240
ꢀꢀꢀ
(9)
[0078]
式中：h为地层深度，m；
[0079]
步骤九：结合钻完井报告的中途测试或静压测试获取的地层压力实测资料，利用式(10)计算求出不同地层压力测试点下的压实校正系数c，发现压实校正系数c值不是固定值，是一个与深度有关的函数，以此建立了压实校正系数c与深度间的关系，如图3所示，见式(11)，
[0080][0081]
式中：p
p
为地层压力，mpa；σ
v
为上覆岩层压力，mpa；p
h
为地层静水压力，mpa；δt
n
为同一深度正常压实趋势线上的声波时差，us/m；c为压实校正系数，
[0082]
c＝5
×
10
‑
13
h
3.5363
ꢀꢀꢀ
(11)
[0083]
式中：h为地层深度，m，
[0084]
其中，
[0085]
[0086]
式中：h0为起始测井点深度，m；为未测井段深度h0点的密度，g/cm3；g为重力加速度，m/s2；den
水化前
为水化校正后的密度测井值，g/cm3。
[0087]
联合公式(8)～公式(12)得到计算孔隙压力的计算式，如表1所示。经过声波时差校正和有机质校正后，地层压力预测结果明显要好于校正前，说明该方法利于提高地层压力的预测精度。
[0088]
表1研究地层校正前后预测压力和实测压力对比
[0089]