1.本发明涉及非常规油气勘探与开发领域，涉及地层岩石的脆性指数计算方法，尤其是一种致密砂岩储层岩石脆性指数计算方法。
背景技术：
2.岩石的脆性指数是致密砂岩油气体积压裂设计中应考虑的重要因素之一。一般认为，脆性指数与岩石矿物成分、岩石力学性质等密切相关。岩石脆性指数越高，越容易破裂形成复杂裂缝。
3.当前，国内外岩石脆性评价方法主要可分为三大类：基于弹性应变的评价方法，基于岩石弹性参数的评价方法和基于岩石矿物组成的评价方法。对于非常规储集层，采用不同脆性评价方法计算得到的脆性指数之间差异往往较大，现有脆性指数评价方法可靠性不强，准确性有待进一步提高。
4.基于弹性应变的脆性评价方法是实验室研究岩石脆性指数的主要手段之一。应力-应变实验根据应力—应变实验分析，当应力由某一初始弹性态加载到峰值强度后，将发生突变而迅速跌落至残余强度面上（如附图1）。通过岩石应力跌落的这一特征，可以将脆性指数b1表示为b1=((b-a)-(c-b))/(b-a)。然而，在实际操作中，经常会遇到缺少残余强度或残余强度不明显的问题。对于该类型应力—应变曲线（如附图2），基于弹性应变的评价方法无法准确计算脆性指数。
技术实现要素：
5.本发明的目的在于提高致密砂岩储层岩石脆性指数计算精度与可靠性，有效解决基于弹性应变的脆性指数评价方法中残余强度面不明显的问题，提出了一种设计合理，通用性强，且计算结果准确可靠的致密砂岩储层岩石脆性指数计算方法。
6.本发明一种致密砂岩储层岩石脆性指数计算方法，包括以下步骤：1）制备应力-应变实验所需致密砂岩试样；2）开展应力-应变实验直至岩样破坏；3）将破坏后的岩样进行处理，制备铸体薄片；4）分析铸体薄片中裂缝分布特征，计算分数维度；5）分析分数维度数据，估算岩石脆性指数，估算模型为脆性指数b=(d-d
min
)/(d
max-d
min
)，其中d
max
为通过步骤4）计算得到的最大分数维度，d
min
为通过步骤4）计算得到的最小分数维度，d为某一试样的分数维度；6）对步骤5中估算的岩石脆性指数进行校正，得到致密砂岩储层岩石脆性指数计算模型b0=(d-d
0min
)/(d
0max-d
0min
)，其中d
0min
为校正后的最小分数维度，d
0max
为校正后的最大分数维度；d
0min
和d
0max
的计算如下式：f(d
0min
，d
0max
)= ，
其中1＜d
0min
＜d
0max
＜2，bi为由铸体薄片分数维度计算得到的第i块岩样脆性指数，b1i为由基于弹性应变计算得到的第i块岩样脆性指数，当f(d
0min
，d
0max
)取最小值为目标函数，得到校正后的d
0min
和d
0max
。
7.优选地，步骤4）中分数维度计算方法：采用光学显微镜观察铸体薄片并拍摄图片，对图片进行处理，得到裂缝分布图，根据裂缝分布图计算分数维度。
8.优选地，步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，选定裂缝长度作为裂缝特征，得到边长为r的正方形网格内裂缝总长度l，在双对数坐标中拟合lnl同lnr的关系，线性拟合后得到 ；其中c为比例常数，d为以裂缝长度为特征计算得到的裂缝分数维度。
9.优选地，步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，选定裂缝数量作为裂缝特征，得到边长为r的正方形网格内裂缝数量n，在双对数坐标中拟合lnl同lnn的关系，线性拟合后得到lnl=c+dlnn；其中c为比例常数，d为以裂缝数量为特征计算得到的裂缝分数维度。
10.优选地，步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，选定裂缝面积作为裂缝特征，得到边长为r的正方形网格内裂缝面积s，在双对数坐标中拟合lnl同lns的关系，线性拟合后得到lnl=c+dlns；其中c为比例常数，d为以裂缝面积为特征计算得到的裂缝分数维度。
11.或者优选地，步骤3）中采用蓝色荧光树脂对破坏后的岩样进行染色，标记裂缝初始形态，切割，打磨，封装，最后形成用于光学显微镜下观察的铸体薄片。
12.优选地，步骤2）中应力-应变实验为单轴压缩实验或三轴压缩实验，将试样装入橡胶套，放入压力机承压板中心，以0.5mpa/s的加载速率对试样加载，直至试样破坏。
13.优选地，步骤2）中应力-应变实验为拉伸实验或剪切实验或切削实验。
14.本发明有效的提高了致密砂岩储层岩石脆性指数计算精度与可靠性，有效解决基于弹性应变的脆性指数评价方法中残余强度面不明显的问题。
15.当残余强度面不明显时，基于弹性应变的脆性评价方法失效，而本发明基于分数维度的计算方法依然可以定量评价脆性。
附图说明
16.图1为致密砂岩存在残余强度的应力-应变曲线。
17.图2为致密砂岩不存在残余强度的应力-应变曲线。
18.图3为破坏后的岩样。
19.图4为本发明对破坏后岩样进行处理得到铸体薄片。
20.图5为本发明采用网格覆盖法计算裂缝分数维度示意图。
21.图6为本发明的脆性指数计算结果与基于弹性应变的脆性指数计算结果的对比图。
具体实施方式
22.本发明一种致密砂岩储层岩石脆性指数计算方法，包括以下步骤：1）制备应力-应变实验所需致密砂岩试样；
2）开展应力-应变实验直至岩样破坏；3）将破坏后的岩样进行处理，制备铸体薄片；4）分析铸体薄片中裂缝分布特征，计算分数维度；5）分析分数维度数据，估算岩石脆性指数，估算模型为脆性指数b=(d-d
min
)/(d
max-d
min
)，其中d
max
为通过步骤4）计算得到的最大分数维度，d
min
为通过步骤4）计算得到的最小分数维度，d为某一试样的分数维度；6）对步骤5中估算的岩石脆性指数进行校正，得到致密砂岩储层岩石脆性指数计算模型b0=(d-d
0min
)/(d
0max-d
0min
)，其中d
0min
为校正后的最小分数维度，d
0max
为校正后的最大分数维度；d
0min
和d
0max
的计算如下式：f(d
0min
，d
0max
)=，其中1＜d
0min
＜d
0max
＜2，bi为由铸体薄片分数维度计算得到的第i块岩样脆性指数，b1i为由基于弹性应变计算得到的第i块岩样脆性指数，当f(d
0min
，d
0max
)取最小值为目标函数，得到校正后的d
0min
和d
0max
。
23.步骤4）中分数维度计算方法：采用光学显微镜观察铸体薄片并拍摄图片，对图片进行处理，得到裂缝分布图，根据裂缝分布图计算分数维度。
24.步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，选定裂缝长度作为裂缝特征，得到边长为r的正方形网格内裂缝总长度l，在双对数坐标中拟合lnl同lnr的关系，线性拟合后得到；其中c为比例常数，d为以裂缝长度为特征计算得到的裂缝分数维度。
25.步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，也选定裂缝数量作为裂缝特征，得到边长为r的正方形网格内裂缝数量n，在双对数坐标中拟合lnl同lnn的关系，线性拟合后得到lnl=c+dlnn；其中c为比例常数，d为以裂缝数量为特征计算得到的裂缝分数维度。
26.步骤4）中采用网格覆盖法计算裂缝分数维度，具体为以铸体薄片中心为基准点，划正方形网格，还可以选定裂缝面积作为裂缝特征，得到边长为r的正方形网格内裂缝面积s，在双对数坐标中拟合lnl同lns的关系，线性拟合后得到lnl=c+dlns；其中c为比例常数，d为以裂缝面积为特征计算得到的裂缝分数维度。
27.步骤3）中采用蓝色荧光树脂对破坏后的岩样进行染色，标记裂缝初始形态，切割，打磨，封装，最后形成用于光学显微镜下观察的铸体薄片。
28.步骤2）中应力-应变实验为单轴压缩实验或三轴压缩实验，将试样装入橡胶套，放入压力机承压板中心，以0.5mpa/s的加载速率对试样加载，直至试样破坏。
29.步骤2）中应力-应变实验也可以为拉伸实验或剪切实验或切削实验。
30.拉伸实验：选取圆柱形岩石试样，然后在万能材料试验机上进行单向拉伸，直至拉断。
31.剪切实验：选取圆柱形岩石试样，然后在万能材料试验机上施加剪切力，直至剪断。
32.切削实验：切削实验在岩石切削实验台上进行。实验时，试样用岩石夹紧装置固定，调节切削齿高度，使切削齿高度低于试样顶部，启动电机使试样向切削齿方向运动，直
至试样破坏。
实施例
33.从岩心库选取致密砂岩地层特定深度岩心，采用直径为25mm的套筒钻头钻取小岩心，切割小岩心端面，得到直径为25mm，高为50mm的圆柱形试样；以单轴压缩实验的方式进行岩样破坏，将试样装入橡胶套，放入压力机承压板中心，以0.5mpa/s的加载速率对试样加载，直至试样破坏；采用蓝色荧光树脂对破坏后的岩样进行染色，标记裂缝初始形态，切割，打磨，封装，最后形成用于光学显微镜下观察的铸体薄片试样，如图4所示，在光学显微镜下，由压缩实验产生的裂缝呈现为蓝色。
34.分析铸体薄片中裂缝分布特征，采用网格覆盖法，以铸体薄片中心为基准点，划正方形网格，以裂缝长度作为裂缝特征，得到边长为r的正方形网格内裂缝总长度l如图5所示，在双对数坐标中拟合lnl同lnr的关系，线性拟合后得到；其中c为比例常数，d为以裂缝长度为特征计算得到的裂缝分数维度。最终得到分数维度数据集（1.22，1.29，1.34，1.34，1.35，1.36，1.41，1.42，1.43，1.65）。
35.统计上步骤计算得到的最大分数维度d
max
（1.65）和最小分度维度d
min
（1.22），采用公式b=(d-1.22)/(1.65-1.22)估算脆性指数；估算得到脆性指数数据集（0， 0.16，0.27，0.29，0.31，0.32，0.45，0.46，0.50，1）。
36.利用基于弹性应变计算得到的脆性指数数据集（0.2, 0.28, 0.32, 0.32, 0.41, 0.41, 0.41, 0.45, 0.51, 0.76）对上步骤中通过裂缝分数维度估算的岩石脆性指数进行校正，得到基于裂缝分数维度的致密砂岩储层岩石脆性指数计算模型b0=(d-d
0min
)/(d
0max-d
0min
)，其中d
0min
为校正后的最小分数维度，d
0max
为校正后的最大分数维度；d
0min
和d
0max
的计算如下式：f(d
0min
，d
0max
)=，其中1＜d
0min
＜d
0max
＜2，bi为由铸体薄片分数维度计算得到的第i块岩样脆性指数，b1i为由基于弹性应变计算得到的第i块岩样脆性指数，当f(d
0min
，d
0max
)取最小值为目标函数，得到校正后的d
0min
和d
0max
；本实施例中通过逐步增大d
max
值与逐步减小d
min
值结合的方式，求解使f(d
0min
，d
0max
)取最小值的最大分数维度d
max
与最小分数维度d
min
，此时最大分数维度d
max
为校正后的最大分数维度d
0max
，最小分数维度d
min
为校正后的最小分数维度d
0min
，求解过程中应保证1＜d
min
＜d
max
＜2。
37.校正过程matlab程序：clear, close all, clcd=[1.22,1.29,1.34,1.34,1.35,1.36,1.41,1.42,1.43,1.65];b=[0.2, 0.28, 0.32, 0.32, 0.41, 0.41, 0.41, 0.45, 0.51, 0.76];f=zeros(23,36);m=0;for mina=1.22:-0.01:1
ꢀꢀꢀꢀ
m=m+1;
ꢀꢀꢀꢀ
n=0;
ꢀꢀꢀꢀ
for maxa=1.65:0.01:2
ꢀꢀꢀꢀꢀꢀꢀꢀ
n=n+1;
ꢀꢀꢀꢀꢀꢀꢀꢀ
err=0;
ꢀꢀꢀꢀꢀꢀꢀꢀ
for i=1:10
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
err=err+((((d(i)-mina))/(maxa-mina))-b(i))^2;
ꢀꢀꢀꢀꢀꢀꢀꢀ
end
ꢀꢀꢀꢀꢀꢀꢀꢀ
f(m,n)=err;
ꢀꢀꢀꢀ
endend[m, n]=find(f==min(min(f)));min_real=1.22-(m-1)*0.01max_real=1.65+(n-1)*0.01) 本实施例校正后得到的最大分数维度d
0max
为1.83，最小分度维数d
0min
为1.07，最终确定基于分形维的脆性指数计算方法b=(d-1.07)/(1.83-1.07)；此外还可以采用基于岩石弹性参数或基于岩石矿物计算得到的脆性指数对模型进行校正。
[0038]
图6为本发明脆性指数计算结果与基于弹性应变的脆性指数计算结果的对比。可以看到，采用本发明计算得到的脆性指数与常规方法即基于弹性应变计算脆性指数的计算结果基本一致。更重要的是，采用本发明方法，可以在残余强度缺失的情况下定量计算脆性指数，为砂岩储层可压性评价提供了一种新手段。