CN107290219B - 储气库盖层岩石的综合脆性评价方法 - Google Patents
储气库盖层岩石的综合脆性评价方法 Download PDFInfo
- Publication number
- CN107290219B CN107290219B CN201710639252.2A CN201710639252A CN107290219B CN 107290219 B CN107290219 B CN 107290219B CN 201710639252 A CN201710639252 A CN 201710639252A CN 107290219 B CN107290219 B CN 107290219B
- Authority
- CN
- China
- Prior art keywords
- brittleness
- rock
- rock sample
- evaluation method
- evaluation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/40—Investigating hardness or rebound hardness
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种储气库盖层岩石的综合脆性评价方法,属岩石脆性评价技术领域。本发明首先,分别运用基于应力‑应变曲线的脆性评价方法、基于硬度的脆性评价方法和基于矿物组成的脆性评价方法对所取岩样的脆性进行评价,然后,运用基于优化理论的权重反分析方法求得三种脆性评价方法的权重值,最后,对待评价储气库盖层岩石的脆性进行综合评价。本发明有效解决了以往单一的脆性评价方法的局限性,实用性强,评价可靠性高,而且,本发明所取岩样一芯多用、实验成本低。
Description
技术领域
本发明涉及岩石脆性评价技术领域,具体涉及的是一种储气库盖层岩石的综合脆性评价方法。
背景技术
地下储气库在国家战略投资、事故应急、协调季节用气高峰低谷期以及优化供气系统等方面有着保护国家能源安全、维护社会稳定、缓解高峰油气资源分配压力、调控市场油气价格等作用,世界各国正大力推进储气库的建设。
岩石脆性是指岩石在发生很小变形即破坏的性质。储气库盖层岩石脆性评价对盖层力学性质研究、储气库最大允许压力确定和盖层封闭能力评价有着重要意义。储气库注气后,储层发生膨胀变形,致使盖层上抬;采气后,储层下沉,盖层随之下沉变形。若盖层脆性较大,变形过程中会产生裂缝,导致盖层渗透率提高、封闭能力变差,从而造成油气泄漏、污染环境、发生爆炸引起人员伤亡等问题。因此,有必要开展储气库盖层岩石脆性评价方法研究。迄今为止,表征盖层岩石的脆性定义和度量还没有统一的标准。现有的评价方法多从硬度或坚固性、强度、矿物组成及应力应变曲线等角度进行,存在着试验成本较高不易实现、适用性不强、缺乏理论依据、理论依据缺乏验证、评价方法过于单一等方面的局限。
发明内容
本发明的目的是针对上述现有技术存在的问题,提供一种实验成本低、实用性强、对岩石脆性评价可靠的储气库盖层岩石的综合脆性评价方法。
为实现上述目的,本发明采用如下技术方案:
储气库盖层岩石的综合脆性评价方法,包括以下步骤:
1)对待评价储气库盖层岩石进行钻井取样并将所取岩样打磨加工分成三类岩样,即完好岩样Pm、残缺岩样Ir和破坏岩样Dn;
2)取步骤1)中的完好岩样Pm进行三轴压缩试验以获取完好岩样Pm的应力-应变曲线,然后,运用基于应力-应变曲线的脆性评价方法评价所取岩样的脆性;
3)取步骤1)中的残缺岩样Ir进行硬度测试以获取残缺岩样Ir的硬度,然后,运用基于硬度的脆性评价方法评价所取岩样的脆性;
4)分别取步骤2)和步骤3)中试验后完好岩样Pm的破坏样和残缺岩样Ir的破坏样以及步骤1)中的破坏岩样Dn,并分别将其研磨成粉末后进行物相分析,然后,运用基于矿物组成的脆性评价方法获取脆性指数并进行所取岩样的脆性评价;
5)运用基于优化理论的权重反分析方法计算步骤2)、步骤3)和步骤4)中三种脆性评价方法的权重值,综合评价待评价储气库盖层岩石的脆性。
进一步地,所述步骤2)中基于应力-应变曲线的脆性评价方法的评价公式为:
式中,B'1为峰后应力降的相对大小,B'2为峰后应力降的绝对速率,B'3为屈服影响系数,τp为峰值强度,τr为残余强度、kac(AC)为屈服阶段起点与残余阶段起点两点连线的斜率,εq1为屈服阶段起点应变、εq2为屈服阶段终点应变、εr为残余应变。
进一步地,所述步骤3)中基于硬度的脆性评价方法的评价公式为:
Br1=S11.2
式中,S11.2为冲击试验后统计得到的粒径小于11.2mm的碎屑百分比;
或为:
Br2=qσc
式中,q为冲击试验后统计得到的粒径小于0.60mm的碎屑百分比,σc为抗压强度,取值为步骤2)中完好岩样Pm三轴压缩试验围压为0MPa时的值。
进一步地,所述步骤4)中基于矿物组成的脆性评价方法的评价公式为:
式中,V石英、V方解石、V黏土、V其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的矿物体积,YM石英、YM方解石、YM黏土、YM其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的弹性模量,PR石英、PR方解石、PR黏土、PR其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的泊松比;
或为:
式中,ai为各种脆性矿物相对石英的脆性系数,其中,ai=Ei/μi,i为矿物种类,Mi为地层中矿物的体积百分含量,Ei为第i种脆性矿物的弹性模量、μi为第i种脆性矿物的泊松比。
本发明具有以下有益效果:
本发明分别运用基于应力-应变曲线的脆性评价方法、基于硬度的脆性评价方法和基于矿物组成的脆性评价方法对所取岩样的脆性进行评价,然后,运用基于优化理论的权重反分析方法求得三种脆性评价方法的权重值,最后,对待评价储气库盖层岩石的脆性进行综合评价,有效解决了以往单一的脆性评价方法的局限性,实用性强,评价可靠性高,而且,本发明所取岩样一芯多用、实验成本低。
附图说明
图1为本发明储气库盖层岩石的综合脆性评价方法的流程图;
图2为本发明储气库盖层岩石的综合脆性评价方法中不同围压的简化岩样应力-应变曲线。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,储气库盖层岩石的综合脆性评价方法,包括以下步骤:
1)对待评价储气库盖层岩石进行钻井取样并将所取岩样打磨加工分成三类岩样,即完好岩样Pm、残缺岩样Ir和破坏岩样Dn;
2)取步骤1)中的完好岩样Pm进行三轴压缩试验以获取完好岩样Pm的应力-应变曲线,然后,运用基于应力-应变曲线的脆性评价方法评价所取岩样的脆性;
3)取步骤1)中的残缺岩样Ir进行硬度测试以获取残缺岩样Ir的硬度,然后,运用基于硬度的脆性评价方法评价所取岩样的脆性;
4)分别取步骤2)和步骤3)中试验后完好岩样Pm的破坏样和残缺岩样Ir的破坏样以及步骤1)中的破坏岩样Dn,并分别将其研磨成粉末后进行物相分析,然后,运用基于矿物组成的脆性评价方法获取脆性指数并进行所取岩样的脆性评价;
5)运用基于优化理论的权重反分析方法计算步骤2)、步骤3)和步骤4)中三种脆性评价方法的权重值,综合评价待评价储气库盖层岩石的脆性。
其中,步骤2)中完好岩样Pm进行三轴压缩试验时,将完好岩样Pm加工成尺寸为直径×高=50mm×100mm的标准圆柱体,试验围压取值范围为0MPa~模拟地层条件下的围压值,每隔5MPa取一个围压值,且围压值数量不得少于5个,试验温度为地表温度+地温梯度×地层埋深,地温梯度取值为25℃/km,各岩样在不同围压下的应力-应变曲线的简化结果如图2所示。
由图2可以看出,岩样在三轴试验获取的应力-应变曲线可简化为由弹性阶段、屈服阶段、软化阶段和残余阶段4个阶段组成。随着围压的升高,峰值强度越来越高、屈服阶段越来越长。这是因为三轴试验时的围压对岩样而言相当于套箍,起到约束的作用,且围压越高、套箍约束作用越强,岩样塑性变形能力越好、脆性越弱。当围压较低时(围压为P1、P2),约束作用较弱,此时,岩样的应力-应变曲线没有表现出屈服阶段,但随着围压的升高,峰值强度越来越大。当围压进一步升高至P3、P4、P5时,随着围压的升高,套箍约束作用越来越强,峰值强度进一步增大,从P3开始出现屈服阶段,且随着围压的升高屈服阶段越来越长,当围压足够大(对应P5)时,岩样的应力-应变曲线屈服阶段足够长,没有出现软化阶段和残余阶段。
将图2中试验数据代入基于应力-应变曲线的脆性评价方法的评价公式,可得不同围压下完好岩样Pm的脆性指数,其中,基于应力-应变曲线的脆性评价方法的评价公式为:
式中,B'1为峰后应力降的相对大小,B'2为峰后应力降的绝对速率,B'3为屈服影响系数,τp为峰值强度,τr为残余强度、kac(AC)为屈服阶段起点与残余阶段起点两点连线的斜率,εq1为屈服阶段起点应变、εq2为屈服阶段终点应变、εr为残余应变。
步骤3)中残缺岩样Ir的硬度测试采用冲击试验,试验后将粒径小于11.2mm及0.60mm的碎屑百分比进行统计,并分别记为S11.2、q,然后,将试验统计结果代入基于硬度的脆性评价方法的评价公式,可得残缺岩样Ir的脆性指数,其中,基于硬度的脆性评价方法的评价公式为:
Br1=S11.2
式中,S11.2为冲击试验后统计得到的粒径小于11.2mm的碎屑百分比;
或为:
Br2=qσc
式中,q为冲击试验后统计得到的粒径小于0.60mm的碎屑百分比,σc为抗压强度,取值为步骤2)中完好岩样Pm三轴压缩试验围压为0MPa时的值。
步骤4)中采用X射线衍射仪技术(X-ray diffraction,XRD)进行岩样的物相分析,试验时分别将步骤2)和步骤3)中试验后完好岩样Pm的破坏样和残缺岩样Ir的破坏样以及步骤1)中的破坏岩样Dn研磨至粒度在1~10微米左右,然后,将测试得到的各岩样的各矿物含量代入到基于矿物组成的脆性评价方法的评价公式,可得所取各岩样的脆性指数,其中,基于矿物组成的脆性评价方法的评价公式为:
式中,V石英、V方解石、V黏土、V其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的矿物体积,YM石英、YM方解石、YM黏土、YM其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的弹性模量,PR石英、PR方解石、PR黏土、PR其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的泊松比;
或为:
式中,ai为各种脆性矿物相对石英的脆性系数,其中,ai=Ei/μi,i为矿物种类,Mi为地层中矿物的体积百分含量,Ei为第i种脆性矿物的弹性模量、μi为第i种脆性矿物的泊松比。
步骤5)中基于优化理论的权重反分析方法将直接求权重问题转化为通过实际工程样本反求各影响因素权重问题,不需先知经验,求解过程完全由实际样本驱动,具有很强的可操作性,其基本原理及求解步骤如下:
假设评价等级向量L=(l1,l2,…,lm-1,lm),等级个数为m,评价指标向量I=(i1,i2,…,in-1,in),指标个数为n,评价指标的权重向量为W=(w1,w2,…,wn-1,wn),且评价指标和评价等级之间的综合评价算法函数关系为f,有:
L=f(I,W)
上式的含义为:由评价指标及评价指标的权重根据一定的综合评价算法函数f算出最终的评价等级,其中,综合评价算法有模糊综合评价法、灰色综合评价法。
将上式中“由评价指标、评价指标的权重及综合评价算法函数f求解评价等级”称为正问题,把“已知评价等级求解评价指标的权重及综合评价算法函数f”称为反问题,基于优化理论的权重反分析方法就是在已知评价等级、评价指标及综合评价算法函数f的前提下,依据一定的优化算法来求解评价指标的权重,其中,优化算法可选用粒子群算法、混沌优化算法,依据一定的优化算法来求解评价指标的权重的具体过程为:
若已知样本集S=(s1,s2,…,sk-1,sk),样本由评价指标值和评价等级组成,样本集S的个数k不应小于所求的权重个数n,则由样本集S求解各评价指标权重的数学优化问题表达式为:
上述基于优化理论的权重反分析方法计算步骤2)、步骤3)和步骤4)中三种脆性评价方法的权重值的实际过程是在确定综合评价算法函数f和优化算法后,在MATLAB平台上编制相应的权重反分析程序WeightCal.m即可求解权重向量。
其他未详细描述的部分均属于现有技术。
Claims (2)
1.储气库盖层岩石的综合脆性评价方法,其特征在于:包括以下步骤:
1)对待评价储气库盖层岩石进行钻井取样并将所取岩样打磨加工分成三类岩样,即完好岩样Pm、残缺岩样Ir和破坏岩样Dn;
2)取步骤1)中的完好岩样Pm进行三轴压缩试验以获取完好岩样Pm的应力-应变曲线,然后,运用基于应力-应变曲线的脆性评价方法评价所取岩样的脆性;
3)取步骤1)中的残缺岩样Ir进行硬度测试以获取残缺岩样Ir的硬度,然后,运用基于硬度的脆性评价方法评价所取岩样的脆性;
4)分别取步骤2)和步骤3)中试验后完好岩样Pm的破坏样和残缺岩样Ir的破坏样以及步骤1)中的破坏岩样Dn,并分别将其研磨成粉末后进行物相分析,然后,运用基于矿物组成的脆性评价方法获取脆性指数并进行所取岩样的脆性评价;
5)运用基于优化理论的权重反分析方法计算步骤2)、步骤3)和步骤4)中三种脆性评价方法的权重值,综合评价待评价储气库盖层岩石的脆性;
所述步骤2)中基于应力-应变曲线的脆性评价方法的评价公式为:
式中,B′1为峰后应力降的相对大小,B'2为峰后应力降的绝对速率,B'3为屈服影响系数,τp为峰值强度,τr为残余强度、kac(AC)为屈服阶段起点与残余阶段起点两点连线的斜率,εq1为屈服阶段起点应变、εq2为屈服阶段终点应变、εr为残余应变;
所述步骤3)中基于硬度的脆性评价方法的评价公式为:
Br1=S11.2
式中,S11.2为冲击试验后统计得到的粒径小于11.2mm的碎屑百分比;
或为:
Br2=qσc
式中,q为冲击试验后统计得到的粒径小于0.60mm的碎屑百分比,σc为抗压强度,取值为步骤2)中完好岩样Pm三轴压缩试验围压为0MPa时的值;
所述步骤4)中基于矿物组成的脆性评价方法的评价公式为:
式中,V石英、V方解石、V黏土、V其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的矿物体积,YM石英、YM方解石、YM黏土、YM其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的弹性模量,PR石英、PR方解石、PR黏土、PR其他脆性矿物依次为石英、方解石、黏土及其他脆性矿物的泊松比;
或为:
式中,ai为各种脆性矿物相对石英的脆性系数,其中,ai=Ei/μi,i为矿物种类,Mi为地层中矿物的体积百分含量,Ei为第i种脆性矿物的弹性模量、μi为第i种脆性矿物的泊松比。
2.根据权利要求1所述的储气库盖层岩石的综合脆性评价方法,其特征在于:所述步骤5)中基于优化理论的权重反分析方法是以综合评判算法和优化算法为基础进行的;所述综合评价算法为模糊综合评价法,或为灰色综合评价法;所述优化算法为粒子群算法,或为混沌优化算法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2017105973487 | 2017-07-20 | ||
CN201710597348 | 2017-07-20 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107290219A CN107290219A (zh) | 2017-10-24 |
CN107290219B true CN107290219B (zh) | 2020-04-07 |
Family
ID=60103916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710639252.2A Active CN107290219B (zh) | 2017-07-20 | 2017-07-31 | 储气库盖层岩石的综合脆性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107290219B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112179770A (zh) * | 2020-09-29 | 2021-01-05 | 西南石油大学 | 基于岩屑微纳米压痕实验的页岩单轴抗压强度评价方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109187916B (zh) * | 2018-07-11 | 2021-08-31 | 中国石油天然气股份有限公司 | 一种判别地层脆性的连续深度处理方法及装置 |
CN109612840B (zh) * | 2018-12-27 | 2021-03-09 | 东北大学 | 用于获得脆性岩石峰后曲线和残留强度的实验装置及方法 |
CN110501758B (zh) * | 2019-07-30 | 2021-06-18 | 东北大学 | 一种砂砾岩储层纵向连续脆性指数预测方法 |
CN111271055B (zh) * | 2020-02-26 | 2021-10-08 | 中国石油大学(北京) | 页岩的脆性指数确定方法、装置和设备 |
US12050297B2 (en) | 2020-09-11 | 2024-07-30 | Saudi Arabian Oil Company | Method and system for determining energy-based brittleness |
CN112253101B (zh) * | 2020-10-09 | 2023-02-03 | 中国石油大学(北京) | 油气资源勘探方法、装置、设备及计算机可读存储介质 |
CN112362475B (zh) * | 2020-11-09 | 2022-05-10 | 重庆科技学院 | 油藏型储气库圈闭密封极限压力的评价方法 |
CN112986247B (zh) * | 2021-04-29 | 2021-08-03 | 中国科学院地质与地球物理研究所 | 基于矿物含量和分布的岩石脆性分析方法、系统及设备 |
CN115795916B (zh) * | 2023-01-09 | 2023-05-05 | 北京科技大学 | 一种岩石脆性评价方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104406849B (zh) * | 2014-11-21 | 2017-01-04 | 中国石油天然气股份有限公司 | 一种储层岩石脆性的预测方法及装置 |
CN104677734A (zh) * | 2015-02-12 | 2015-06-03 | 中国石油大学(华东) | 一种多尺度信息融合的致密地层脆性指数测量方法 |
CN106872260B (zh) * | 2017-03-09 | 2019-11-29 | 成都理工大学 | 一种岩石脆性指数的获取方法及岩石的脆性评价方法 |
-
2017
- 2017-07-31 CN CN201710639252.2A patent/CN107290219B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112179770A (zh) * | 2020-09-29 | 2021-01-05 | 西南石油大学 | 基于岩屑微纳米压痕实验的页岩单轴抗压强度评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107290219A (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107290219B (zh) | 储气库盖层岩石的综合脆性评价方法 | |
CN111238931B (zh) | 基于能量演化的页岩脆性指数评价方法 | |
CN105866835A (zh) | 一种基于地应力分布的断层三维封闭性定量评价方法 | |
CN109324345B (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN108827774A (zh) | 煤岩脆性评价方法 | |
CN110715859B (zh) | 一种基于弹塑性变形和破裂强度的脆性指数评价方法 | |
CN110057675A (zh) | 一种土体受扰动后的e-logp’曲线的确定方法 | |
CN111077027A (zh) | 一种高应力下钙质砂的剪切强度确定方法 | |
Luo et al. | Dynamic mechanical properties and constitutive model of shale with different bedding under triaxial impact test | |
Tang et al. | Study on the microscopic fracture process and acoustic emission of shale based on digital image | |
Wei et al. | Investigation and prediction of thermal cracking and permeability enhancement in ultra-low permeability rocks by in-situ thermal stimulation | |
Gao et al. | Mechanical characterization of uniaxial compression associated with lamination angles in shale. | |
Xia et al. | Dynamic tensile test of coal, shale and sandstone using split Hopkinson pressure bar: a tool for blast and impact assessment | |
Dou et al. | Study on the Effect of High‐Temperature Heat Treatment on the Microscopic Pore Structure and Mechanical Properties of Tight Sandstone | |
Gao et al. | New mathematical models for calculating the proppant embedment and conductivity | |
Ren et al. | Hydraulic fracture initiation and propagation in deep coalbed methane reservoirs considering weak plane: CT scan testing | |
Cai et al. | Triaxial Compressive Failure Characteristics and Constitutive Model Study of Jurassic‐Cretaceous Weakly Cemented Sandstone | |
Meng et al. | Experimental study on the mechanical behavior of different lithologic shales in Qingshankou group under in-situ geological conditions | |
Xu et al. | Experimental research on brittleness and rockburst proneness of three kinds of hard rocks under uniaxial compression | |
Tan et al. | A classification and description method for carbonate reservoir mechanical properties based on the Hoek-brown criterion | |
Zhang et al. | Study on impact tendency of coal and rock mass based on different stress paths | |
CN112881184A (zh) | 基于全程应力应变的砂砾岩脆性指数确定方法及系统 | |
Liu et al. | Mechanical Behavior and Impact Resistance Tendency of Coal Mass with Different Water Content | |
Guo et al. | The influence of fracture strain energy on the burst tendency of coal seams and field application | |
Wen et al. | Experimental study on the evolution of rock strength and sealing properties under multi-cycles alternating loading |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |