一种页岩水力压裂过程中页岩压裂脆性指数的预测方法

1.本发明涉及一种页岩水力压裂过程中页岩压裂脆性指数的预测方法，属于石油天然气开采技术领域。


背景技术：

2.页岩气作为非常规油气资源的重要组成部分，由于其资源量的巨大以及能源的清洁特性得到了人们的广泛重视与开发。因其储层压裂效果与脆性有着很紧密的联系，于是利用脆性的方法来评价储层裂缝的复杂程度是目前常用的手段，正确知晓目的层段的脆性大小有助于更好的节约成本和提高产量。
3.伴随着各国对页岩气的高度开发，脆性指数的建立种类也越来越多。javie等人(2007)认为页岩的脆性取决于脆性矿物的占比，主要取决于石英含量的多少来判断地层的脆性；richman等人(2008)根据某一特定页岩区块进行试验得出杨氏模量和脆性、泊松比与脆性的散点图，发现页岩脆性随着岩石泊松比的减小而增大以及杨氏模量的增大而增大；其余主要的还有对地层同时利用矿物法与弹性参数法，并对二者的预测结果与实际试验出的情况做拟合，以求预测结果与真实值更加接近。
4.而实际情况在使用这些方法时，会出现以下问题：
5.(1)在不同地区的页岩，其主要脆性物质不一定是石英，而各种矿物所占脆性大小也很难确定，常规的矿物法很难给出一个精准的目的层段页岩脆性预测。
6.(2)rickman所提出来的弹性参数法，考虑的因素太少，仅有泊松比和弹性模量，然而实际情况却是随着页岩层段的深度增加，围压的影响远超弹性模量和泊松比，并且不同的页岩区域表现出来的情况也不一样，很难找到这些力学参数之间的规律。
7.(3)采用拟合法的前提是需要一个有一定表征地层脆性的方法做支撑，在矿物法和弹性参数法都偏离实际情况时，拟合起来的效果也不见得好，并且拟合法本身就具有一定的概率性，具有风险。
8.总的来说前面所提方法的使用都会有地域限制，这些方法的使用在不同区域将会有很大的差异性，不具有普遍适用性。


技术实现要素：

9.本发明的目的是针对现有的方法大多不能有效地针对不同区域计算出准确有效的脆性指数，所提供的一种页岩水力压裂过程中页岩压裂脆性指数的预测方法。
10.本发明解决上述技术问题所提供的技术方案是：一种页岩水力压裂过程中页岩压裂脆性指数的预测方法，包括：
11.分别测定各组岩芯的各矿物组分的百分含量；
12.利用三轴岩石力学测试系统分别测定各组岩芯的力学性质；
13.根据各组岩芯的各矿物组分的百分含量、各组岩芯的力学性质确定各个矿物与力学性质之间的相关性；
14.根据各个矿物与力学性质之间的相关性确定各个矿物的标准量；
15.根据各组岩芯的各矿物组分的百分含量确定各组岩芯的新脆性矿物百分含量；
16.根据各组岩芯的新脆性矿物百分含量、每个矿物的标准量确定各组岩芯的相差量；
17.根据各组岩芯的相差量确定各组岩芯的矿物脆性指数；
18.根据各组岩芯的矿物脆性指数确定各组岩芯的页岩脆性指数。
19.进一步的技术方案是，所述各矿物组分的百分含量m
j
包括石英百分含量m1、钾长石百分含量m2、斜长石百分含量m3、黄铁矿百分含量m4、方解石百分含量m5、白云石百分含量m6以及其他矿物组分百分含量m7；
20.其中m1+m2+m3+m4+m5+m6＝100％。
21.进一步的技术方案是，根据各组岩芯的各矿物组分的百分含量、各组岩芯的杨氏模量和泊松比确定各个矿物与力学性质之间的相关性包括：
22.根据各组岩芯的各矿物组分的百分含量、各组岩芯的杨氏模量和泊松比确定各个矿物与杨氏模量和泊松比相关性的图；
23.再在各个矿物与杨氏模量和泊松比相关性的图上的数据进行拟合，获得各个矿物的杨氏模量相关性拟合直线、泊松比相关性拟合直线；
24.根据各个矿物的杨氏模量相关性拟合直线、泊松比相关性拟合直线确定各个矿物与力学性质之间的相关性。
25.进一步的技术方案是，各个矿物与力学性质之间的相关性的计算公式如下：
26.θ
j
＝θ
aj
+θ
bj
27.θ
aj
＝|tan
‑1(w
1j
)|
28.θ
bj
＝|tan
‑1(w
2j
)|
29.式中：θ
j
为各个矿物与力学性质之间的相关性；θ
aj
为各个矿物与杨氏模量之间的相关性；θ
bj
为各个矿物与泊松比之间的相关性；w
1j
为各个矿物的杨氏模量相关性拟合直线斜率；w
2j
为各个矿物的泊松比相关性拟合直线斜率。
30.进一步的技术方案是，各个矿物的标准量的计算公式如下：
[0031][0032]
式中：θ
j
为各个矿物与力学性质之间的相关性；b
j
为各个矿物的标准量；j为矿物组分的排位数。
[0033]
进一步的技术方案是，各组岩芯的新脆性矿物百分含量的计算公式如下：
[0034][0035]
式中：a
j
为各组岩芯的新脆性矿物百分含量；m
j
为各矿物组分的百分含量。
[0036]
进一步的技术方案是，所述相差量的计算公式为：
[0037]
[0038]
式中：g为相差量；a
j
为各组岩芯的新脆性矿物百分含量；b
j
为各个矿物的标准量。
[0039]
进一步的技术方案是，所述矿物脆性指数的计算公式如下：
[0040][0041][0042]
式中：g为相差量；b为矿物脆性指数；g
max
为最大相差量；b
j
为各个矿物的标准量。
[0043]
进一步的技术方案是，所述页岩脆性指数的计算公式为：
[0044]
b
i
＝b
·
(1
‑
m7)
[0045]
式中：b
i
为页岩脆性指数；b为矿物脆性指数；m7为其他矿物组分百分含量。
[0046]
本发明具有以下有益效果：本发明的方法利用矿物力学非均质性和不同矿物对页岩力学性质的不同贡献度来将页岩的矿物特征和力学特征做到平等化和区域化，得到适用于不同区域的脆性指数计算方法，这种方法将会为目的储层提供较为准确的脆性指数计算方法，解决了现有技术方法不具有普遍适用性的问题。
附图说明
[0047]
图1是矿物含量与力学性质之间的相关性；
[0048]
图2是实施例1中石英矿物对杨氏模量和泊松比的影响图；
[0049]
图3是实施例1中长石矿物对杨氏模量和泊松比的影响图；
[0050]
图4是实施例1中碳酸盐矿物对杨氏模量和泊松比的影响图。
具体实施方式
[0051]
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0052]
本发明的一种页岩水力压裂过程中页岩压裂脆性指数的预测方法，包括以下步骤：
[0053]
步骤s1、分别测定各组岩芯的各矿物组分的百分含量；六种脆性矿物组分石英、钾长石、斜长石、黄铁矿、方解石、白云石的百分含量分别记为：m1、m1、m2、m3、m5、m6；其他矿物组分含量为m7；其他矿物指的是包括粘土在内的所有对矿物脆性有抑制作用的矿物；m1+m2+m3+m4+m5+m6＝100％；
[0054]
步骤s2、利用三轴岩石力学测试系统分别测定各组岩芯的力学性质；
[0055]
步骤s3、根据各组岩芯的各矿物组分的百分含量、各组岩芯的力学性质确定各个矿物与力学性质之间的相关性；
[0056]
其中测定的力学性质主要包括杨氏模量和泊松比与矿物百分含量之间的关系，由rickman脆性指数可知，脆性大小与杨氏模量成正比，与泊松比成反比，于是可分别对两者进行数据拟合，以两条线相交处的角度来表征每种矿物对力学性质的影响程度如何，确定方式如图1所示；
[0057]
由θ
j
来表征相关性的大小，左边的斜线表示矿物含量与杨氏模量之间的相关性，
右边的斜线表示矿物含量与泊松比之间的相关性，θ
j
越大，则此矿物的力学性质影响程度越大；
[0058]
θ
j
＝θ
aj
+θ
bj
[0059]
θ
aj
＝|tan
‑1(w
1j
)|
[0060]
θ
bj
＝|tan
‑1(w
2j
)|
[0061]
式中：θ
j
为各个矿物与力学性质之间的相关性；θ
aj
为各个矿物与杨氏模量之间的相关性；θ
bj
为各个矿物与泊松比之间的相关性；w
1j
为各个矿物的杨氏模量相关性拟合直线斜率；w
2j
为各个矿物的泊松比相关性拟合直线斜率；
[0062]
步骤s4、根据各个矿物与力学性质之间的相关性确定各个矿物的标准量；
[0063]
标准量即为在此区域下每种矿物的含量刚好使得矿物的脆性指数达到最大时各个矿物的百分含量，此时，页岩的非均质性最强；页岩中具有脆性的物质为石英、钾长石、斜长石、黄铁矿、方解石和白云石这六种；
[0064]
各个矿物的标准量的计算公式如下：
[0065][0066]
式中：θ
j
为各个矿物与力学性质之间的相关性；b
j
为各个矿物的标准量；
[0067]
其中，石英b1、钾长石b2、斜长石b3、方解石b4、白云石b5、黄铁矿b6；
[0068]
步骤s5、根据各组岩芯的各矿物组分的百分含量确定各组岩芯的新脆性矿物百分含量；
[0069]
剔除非脆性矿物m7的影响，以六种脆性矿物的百分含量和为分母，各个矿物的百分含量分别作为分子，重新计算百分含量为；
[0070][0071]
式中：a
j
为各组岩芯的新脆性矿物百分含量；m
j
为各矿物组分的百分含量；
[0072]
其中，石英a1、钾长石a2、斜长石a3、方解石a4、白云石a5、黄铁矿a6；
[0073]
步骤s6、根据各组岩芯的新脆性矿物百分含量、每个矿物的标准量确定各组岩芯的相差量；
[0074]
相差量即为脆性矿物百分含量与矿物的标准量之差的绝对值之和，计算如下：
[0075][0076]
式中：g为相差量；a
j
为各组岩芯的新脆性矿物百分含量；b
j
为各个矿物的标准量；
[0077]
步骤s7、根据各组岩芯的相差量确定各组岩芯的矿物脆性指数；
[0078]
相差量g越小，矿物的脆性越强，最小值为0；相差量g越大，矿物的脆性越弱。于是，脆性矿物的脆性指数与相差量成反比，令脆性指数b在0
‑
1之间变化，此时需要求得相差量的最大值，计算过程如下：
[0079]
令b
min
＝(b
j
)
min
，那么改进为那么脆性指数b的表达式为：
[0080]
式中：g为相差量；b为矿物脆性指数；g
max
为最大相差量；b
j
为各个矿物的标准量；
[0081]
步骤s8、根据各组岩芯的矿物脆性指数确定各组岩芯的页岩脆性指数；
[0082]
其他矿物组分主要为高岭石、绿泥石等粘土矿物成分，这些矿物组分的存在会影响页岩总的脆性大小，对脆性起到抑制作用。于是，针对于此，将其他矿物组分的m7作为无效含量，类比为页岩矿物中由m7的区域为无效区域，要将这部分区域的影响考虑进去，于是页岩脆性指数的表达式为：
[0083]
b
i
＝b
·
(1
‑
m7)
[0084]
式中：b
i
为页岩脆性指数；b为矿物脆性指数；m7为其他矿物组分百分含量。
[0085]
实施例1
[0086]
将本发明的页岩压裂脆性指数的预测方法应用于具体页岩样品的测试，具体步骤如下：
[0087]
(1)取永页2井岩芯42块进行岩石矿物分析，编号为1
‑
42，采用x射线衍射仪测定各种矿物的组分含量，具体数据如表1所示。
[0088]
表1页岩矿物百分含量
[0089]
[0090]
[0091][0092]
(2)对42块岩芯利用三轴岩石力学测试系统测定其对应的杨氏模量和泊松比，具体数据如表2所示。
[0093]
表2三轴力学实验实验结果
[0094]
[0095]
[0096][0097]
于是可以结合表1和表2的数据做出各个矿物与杨氏模量和泊松比相关性的图，为了方便计算，将钾长石与斜长石统称为长石一类，剩下的矿物除石英外统称为碳酸盐矿物，如图2
‑
4所示。
[0098]
于是，由图中数据可以得出，
[0099]
θ1＝tan
‑1(75.28)+|tan
‑1(
‑
0.0012)|＝89.24
°
+0.07
°
＝89.31
°
[0100]
θ2＝|tan
‑1(
‑
0.0022)|+|tan
‑1(
‑
55.359)|＝0.13
°
+88.79
°
＝88.92
°
[0101]
θ3＝|tan
‑1(
‑
0.0022)|+|tan
‑1(
‑
55.359)|＝0.13
°
+88.79
°
＝88.92
°
[0102]
θ4＝tan
‑1(119.83)＝89.52
°
[0103]
θ5＝tan
‑1(119.83)＝89.52
°
[0104]
θ6＝tan
‑1(119.83)＝89.52
°
[0105]
(3)确定标准量，其中于是，
[0106][0107][0108][0109]
(4)确定新的脆性矿物百分含量，如表3所示。
[0110]
表3脆性矿物百分含量
[0111]
[0112]
[0113][0114]
(5)计算相差量g，其中b
min
＝0.166，那么g
max
＝2
×
(0.167+0.166+3
×
0.167)＝1.668，各组岩芯的相差量g由表4所示。
[0115]
表4相差量g
[0116]
序号相差量g(％)序号相差量g(％)152.32259.3249.32371.1364.12444.9448.52541.7
573.82639.3662.42763.3750.92836.98502918.5959.530112.71048.231100.81148.73266.31252.13361.31372.13450.414593527.71549.63636.91651.63755.217503899.21839.63991.21964.640372081.44145.82176.74244
[0117]
(6)确定矿物矿物脆性表达式b，于是每组的脆性指数由表5所示。
[0118]
表5矿物脆性指数b
[0119][0120][0121]
(7)考虑其他矿物对页岩脆性指数的负影响，计算公式为：b
i
＝b
·
(1
‑
m7)，表示结果如表6所示。
[0122]
表6页岩脆性指数b
i
[0123][0124][0125]
由表6可以看出，最大脆性指数为序号18岩芯，最小脆性指数为序号13和序号30岩芯，于是可以得出18号岩芯对应区域的脆性最强，13号岩芯和30号岩芯对应区域的脆性最弱。
[0126]
综上所述，本发明发明了一种用于页岩水力压裂过程中页岩脆性指数的预测方法。该方法步骤清晰，且使用起来能够根据各个区域的不同而相应的变化形式，具有广泛的适用性。
[0127]
以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。