非常规油气储层岩石脆性定量评价方法

本发明涉及非常规洁净环保能源勘探开发，特别是一种非常规油气储层岩石脆性定量评价方法。

背景技术：

1、随着社会的高速发展，能源需求量日益增常，非常规油气资源潜力巨大，引起了国内外的广泛关注。非常规油气储层的低孔、低渗特点决定了其在自然状态下的开采有一定的难度，实现非常规油气储层的经济有效开发，依赖于压裂增产改造措施。近年来，随着水力压裂技术的不断更新，已成为此类油气藏改造的主流技术，而脆性评价是压裂改造和储层评价的关键参数。

2、由于脆性的影响因素众多，目前对于脆性的定义尚未统一，国内外学者根据自己的研究内容提出了多种脆性评价方法。对于脆性评价，本质上包含的信息越多，评价越准确，也就要求需要将多种因素进行综合考虑分析，然而这些因素存在相互影响的关系，简单的将多种因素组合起来，会导致信息的重复，因此探索出一种更加合理的脆性评价方法，值得去研究。

3、正交试验设计(orthogonal experimental design)是一种利用标准化的正交表来安排试验方案，研究多因素多水平的一种设计方法，根据正交性从全面试验中挑选出部分有代表性的点进行试验，这些点具备了"均匀分散，齐整可比"的特点，从而实现减少实验次数的目的。对正交实验数据的评价分析，探究各因素对实验结果的影响程度，判断出各因素对试验结果影响的主次顺序，从而达到筛选主控因素的目的。

4、网络层次分析法(anp)是美国匹兹堡大学的t·l·saaty教授提出的在层次分析法(analytic hierarchy process,简称ahp)的基础上进行改进的一种新的决策方法，anp的网络层次结构相对于ahp递阶层次结构更为复杂，既存在递阶层次结构，又存在内部循环相互支配的层次结构，而且层次结构内部还存在依赖性和反馈性，网络层次分析法能够更好的表征因素间的耦合关系。脆性评价方法包括基于岩石所表现出的力学响应参数建立的弹性参数法、抗压抗拉强度法、应力应变特征法，基于岩石矿物组成的矿物组分法。传统的单一因素的脆性评价方法，并不能完整的表征岩石的脆性特点，综合考虑多因素的传统脆性评价方法，无法厘清成岩过程及加载过程中的脆性主控因素，没有解耦各因素之间的交互耦合影响作用，不能准确的表征岩石的脆性。

技术实现思路

1、为解决现有技术中存在的问题，本发明的目的是提供一种非常规油气储层岩石脆性定量评价方法，本发明综合考虑成岩过程及加载全过程各阶段脆性影响因素的作用，遴选出各阶段主控因素，同时考虑了各阶段主控因素之间的交互耦合关系，实现了主控因素互相影响的解耦及推算了脆性特征贡献的权重，建立了一种新型脆性定量评价方法，也为非常规油气储层岩石油气的高效勘探开发提供理论基础。

2、为实现上述目的，本发明采用的技术方案是：一种非常规油气储层岩石脆性定量评价方法，包括以下步骤：

3、s1、综合考虑非常规油气储层岩石油气储层勘探及开发过程中脆性特征影响因素及响应特征，将脆性特征影响过程分为成岩过程，加载弹性过程及加载损伤过程，初步遴选各阶段脆性影响成岩参数，弹性参数及破裂参数，建立一级指标层；

4、s2、基于成岩过程及加载全过程脆性影响成岩参数，弹性参数及破裂参数特点，基于实验数据结合正交实验设计算法及极差分析对影响脆性的各成岩参数进行重要程度排序，筛选并确定脆性特征成岩主控参数；基于加载过程应力应变曲线变化规律确定弹性阶段主控参数及损伤过程破裂主控参数；

5、s3、将各阶段主控因素设立为二级指标层，并与一级指标层形成非独立递阶层次结构的网络层次结构模型；明确各阶段主控因素之间的耦合关系，采用网络层次分析算法，基于专家评分及实践经验，采用九分法构造优势度判断矩阵，经过一致性检验后构造超矩阵，求解加权超矩阵解耦各主控因素间的相互影响，并进一步推算各阶段主控参数对脆性的贡献权重；

6、s4、输入实验数据，计算各阶段脆性主控参数及贡献权重大小，结合成岩及加载全过程交互影响主控因素解耦赋权的脆性定量评价模型，得到最终脆性指数。

7、作为本发明的进一步改进，在步骤s1中，所述加载包括拉应力加载和压应力加载；所述成岩参数包括矿物组分类别及含量、成岩温度、成岩压力和胶结程度，其中所述矿物组分类别包括石英、方解石、长石、黏土、黄铁矿和toc；所述弹性参数包括杨氏模量，泊松比、阻抗、刚量和纵横波速度；所述破裂参数包括抗压强度、抗拉强度、内摩擦角和残余强度。

8、作为本发明的进一步改进，在步骤s2中，所述实验数据包括xrd矿物组分全岩分析数据、toc含量数据和成岩压力数据；所述正交实验设包括：选择原则是依据正交表均匀分散，齐整可比的特点，选用较小的正交表，通过部分试验分析了解全面试验的情况，然后找到最优水平组合，并基于极差分析确定因素的主次影响顺序，实现筛选成岩过程脆性影响主控因素的目的；所述加载过程应力应变曲线变化规律包括：加载弹性过程及加载损伤破裂过程，以及应力增大过程中的径向应变及轴向应变响应规律，抗压及抗拉破裂强度响应规律。

9、作为本发明的进一步改进，在步骤s3中，非独立递阶层次结构表示各指标以及指标集之间存在相互耦合影响且层次递阶；各阶段主控因素之间的耦合关系表示成岩、弹性及破裂某一阶段的主控因素对脆性的影响又受另外两阶段主控因素影响和支配，存在相互耦合影响作用的情况。

10、作为本发明的进一步改进，在步骤s3中，采用九分法构造优势度判断矩阵中，九分法为采用数字1-9来表示两个元素的重要程度，1表示同等重要，数字越大重要性越强。

11、作为本发明的进一步改进，在步骤s3中，所述一致性检验为判断优势度判断矩阵是否满足一致性的要求，具体如下：

12、

13、

14、

15、其中，i指矩阵的第i个特征根，为第i个特征根的特征矢量，λmax为最大特征根，sbw,i为优势度判断矩阵第i列数值之和，n为优势度判断矩阵的阶数；c.i.为一致性指标，ri为随机一致性指标；c.r.为一致性比率；c.i.值越大，表明优势度判断矩阵偏离完全一致性的程度越大；c.i.值越小(越接近于0)，表明优势度判断矩阵的一致性越好；当n＝1或2时，c.i.为0，此时优势度判断矩阵完全一致；当n≥3时，再继续求解c.r.，若c.r.≤0.1，则优势度判断矩阵满足一致性要求。

16、作为本发明的进一步改进，在步骤s3中，所述加权超矩阵的建立以一级指标下的二级指标优势度比较，确定优势度判断矩阵，通过特征根法求得排序向量，如果特征向量经检验符合一致性判断，将其改写成矩阵形式，得到局部权重向量矩阵，将局部矩阵中的列向量和行向量进行重要度排序，建立超矩阵w：

17、

18、超矩阵列向量是采用特征向量法得到的归一化特征向量值，其中wiji＝1,2,..n.,；j＝1,2,...,n的列向量表示指控因素i对主控因素j中参数影响程度的排序向量；

19、将超矩阵中每一个子块wij都经过归一化处理得到一个归一化的排序向量：aj＝[a1j,a2j,...,anj]t，由此得出加权矩阵a：

20、

21、加权矩阵a与子块wij相乘得到加权超矩阵

22、

23、作为本发明的进一步改进，步骤s4中，所述脆性定量评价模型为根据网络层次分析法求解极限加权超矩阵获得的脆性主控因素权重值与归一化的实验值，综合评价表示如下：

24、b＝∑aibi

25、其中，b为最终的脆性综合指数，ai为第i项的指标权重，bi为第i项指标值。

26、本发明基于正交试验设计筛选成岩参数中的脆性主控因素，基于加载过程应力应变曲线变化规律确定弹性阶段主控参数及损伤过程破裂主控参数，基于网络层次分析法准确全面地表征并解耦各阶段脆性主控因素之间的交互影响关系，进一步推算各阶段主控因素的贡献权重构建新型脆性定量评价模型，从而使脆性预测结果更加准确。

27、本发明的有益效果是：

28、1、本发明综合考虑了非常规油气储层成岩过程及加载全过程各阶段脆性影响因素的作用，有机分类了成岩过程脆性影响参数，加载弹性过程脆性影响参数，加载损伤破裂过程脆性影响参数，采用正交实验设计方法遴选了成岩过程脆性影响主控参数，减少实验次数的基础上能够定量表征各因素的相关性和各因素对实验结果的影响程度，实现此过程筛选脆性主控参数的目的；基于加载过程应力应变曲线变化规律确定弹性阶段主控参数及损伤过程破裂主控参数，使脆性评价结果更加可靠；

29、2、本发明建立了非常规油气储层岩石脆性评价的非独立递阶层次的网络层次结构模型，基于网络层层次分析算法框架，结合九分法判断矩阵及加权超矩阵求解，全面完整的表征并解耦了各阶段脆性主控因素间的交互影响关系，降低主控因素影响的重复性，使脆性评价结果更加可靠。

30、3、本发明提出的基于不同阶段脆性主控因素解耦的脆性定量评价方法，求解了各主控因素对脆性响应的贡献权重，并在赋权过程中结合了主观经验，可以针对不同地区作出相应调整，结合实验数据，能够大大提高脆性预测的准确性，同时也为脆性评价方法提供了新的研究方向。