CN110726608B - 一种基于应力-应变曲线能量演化的页岩脆性评价方法 - Google Patents
一种基于应力-应变曲线能量演化的页岩脆性评价方法 Download PDFInfo
- Publication number
- CN110726608B CN110726608B CN201911087573.1A CN201911087573A CN110726608B CN 110726608 B CN110726608 B CN 110726608B CN 201911087573 A CN201911087573 A CN 201911087573A CN 110726608 B CN110726608 B CN 110726608B
- Authority
- CN
- China
- Prior art keywords
- energy
- brittleness
- stress
- peak
- shale
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 239000011435 rock Substances 0.000 claims abstract description 40
- 238000013210 evaluation model Methods 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims description 8
- 230000006378 damage Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 4
- 230000000977 initiatory effect Effects 0.000 abstract description 4
- 238000012512 characterization method Methods 0.000 abstract description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 10
- 239000007789 gas Substances 0.000 description 5
- 239000003345 natural gas Substances 0.000 description 5
- 238000005336 cracking Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Pathology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Medicinal Chemistry (AREA)
- Fluid Mechanics (AREA)
- Food Science & Technology (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于应力‑应变曲线能量演化的页岩脆性评价方法,包括以下步骤:基于应力‑应变曲线能量演化,建立页岩脆性评价模型,根据所述页岩脆性评价模型,评价页岩的脆性。本发明克服了现有技术中只考虑了峰后应力‑应变状态对于岩体脆性的影响,未考虑峰前特性对岩体脆性的表征的缺陷;克服了现有技术中引入起裂应力参数,但由于起裂应力与起裂应变取值方法较复杂且取值精度低,导致工程应用难度大的缺陷。本发明参数取值方便,脆性指数计算结果与实验现象吻合,能够准确地反映岩石的脆性特征,体现不同岩样间脆性差异,为水力压裂提供重要的技术支撑。
Description
技术领域
本发明涉及页岩脆性评价技术领域,特别涉及一种基于应力-应变曲线能量演化的页岩脆性评价方法。
背景技术
2018年我国天然气进口量约1254亿立方米,成为全球第一大天然气进口国,对外依存度升至45.3%。然而事实上,我国的页岩气储量高居全球第一位,可采资源量达21.8万亿立方米,开发潜力巨大。因此非常规资源尤其是页岩气资源的开发价值不言而喻。
页岩储层具有低孔低渗的特点,开发难度大,因此必须对其进行增产改造。水力压裂是一种对页岩储层改造以提高采收率的有效方法,主要手段是通过向地层注入超过其吸液能力的高粘压裂液,从而在储层中压出裂缝,形成复杂的裂缝网络。而脆性已被证实与复杂缝网的形成密切相关。岩石脆性越大,破裂后形成的裂缝越多,形成有效复杂裂缝网络的面积也越大;并且裂缝发育会促进页岩储层中游离态天然气体积的增加和吸附态天然气的解吸附,有助于提高气藏富集程度,增加产气量。因此,准确评价岩石脆性,对提高页岩气采收率至关重要。
现有技术中,周辉,孟凡震,张传庆等.基于应力–应变曲线的岩石脆性特征定量评价方法[J].岩石力学与工程学报,2014,33(06):1114-1122中公开了一种脆性评价模型,但该模型只考虑了峰后应力-应变状态对于岩体脆性的影响,没有考虑峰前特性对岩体脆性的表征,因此相对于整个应力-应变过程并不全面;陈国庆,赵聪,魏涛.基于全应力–应变曲线及起裂应力的岩石脆性特征评价方法[J].岩石力学与工程学报,2018,37(01):51-59中公开了一种脆性评价模型,该模型引入了起裂应力参数,用峰后应力跌落速率表征峰后脆性指数,用峰前起裂点至峰值点之间的应力增长速率表征峰前脆性指数,然而实际应用发现,其中起裂应力与起裂应变取值方法较复杂且取值精度低,将直接导致该模型计算的脆性评价结果误差较大,工程应用难度大。
发明内容
针对上述问题,本发明旨在提供一种能够准确评价岩石脆的一种基于应力-应变曲线能量演化的页岩脆性评价方法。
本发明的技术方案如下:
一种基于应力-应变曲线能量演化的页岩脆性评价方法,包括以下步骤:
基于应力-应变曲线能量演化,建立页岩脆性评价模型,所述页岩脆性评价模型为:
式中:
Bi为脆性评价指标,无量纲;
Bi1为基于应力-应变曲线能量演化的能量指数,无量纲;
Bi2为考虑峰后跌落应力及峰后跌落应变大小的峰后脆性指标,无量纲;
E为杨氏弹性模量,MPa;
M为峰后弹性模量,MPa;
M=dσ/dε (2)
式中:
d为微分符号;
σ为应力,MPa;
ε为应变,无量纲;
根据所述页岩脆性评价模型,评价页岩的脆性。
进一步地,所述能量指数的计算方法为:
式中:
dwr为峰值后的破裂能,MPa;
dwe为岩石破坏过程中被消耗的弹性能,MPa。
进一步地,所述弹性能的计算方法为:
式中:
σB为峰值应力,MPa;
σC为残余应力,MPa。
进一步地,对于I类岩石,所述破裂能为所述弹性能与附加能之和,所述附加能为峰后破坏过程中的弹性能不足以产生破坏而需要的额外能量;而对于II类岩石,所述破裂能为所述弹性能与释放能之差,所述释放能为峰后破坏过程中释放的多余能量;所述附加能或所述释放能的计算表达式如下:
因M=dσ/dε,故I类岩石中M<0、dwa<0,故II类岩石中M>0、dwa>0,因此所述破裂能的计算式为:
式(5)、式(6)中:dwa为附加能或释放能,MPa。
进一步地,所述峰后脆性指标的计算方法为:
式中:
εC为残余应变,无量纲;
εB为峰值应变,无量纲。
与现有技术相比,本发明具有如下优点:
本发明建立的模型克服了现有技术中只考虑了峰后应力-应变状态对于岩体脆性的影响,未考虑峰前特性对岩体脆性的表征的缺陷;克服了现有技术中引入起裂应力参数,但由于起裂应力与起裂应变取值方法较复杂且取值精度低,导致工程应用难度大的缺陷。本发明参数取值方便,脆性指数计算结果与实验现象吻合,能够准确地反映岩石的脆性特征,体现不同岩样间脆性差异,为水力压裂提供重要的技术支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为第I类和第II类岩石破坏过程的应力-应变曲线能量演化;
图2为一个具体实施例的页岩脆性评价模型计算结果示意图;
图3为当峰值应力作为单因素变化时,峰值应力与脆性指数关系示意图;
图4为当峰值应变作为单因素变化时,峰值应变与脆性指数关系示意图;
图5为当残余应力作为单因素变化时,残余应力与脆性指数关系示意图;
图6为当残余应变作为单因素变化时,残余应变与脆性指数关系示意图;
图7为当弹性模量作为单因素变化时,弹性模量与脆性指数关系示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
一种基于应力-应变曲线能量演化的页岩脆性评价方法,包括以下步骤:
基于如图1所示的第I类和第II类岩石破坏过程的应力-应变曲线能量演化,建立页岩脆性评价模型,根据所述页岩脆性评价模型,评价页岩的脆性。
所述页岩脆性评价模型为:
式中:
Bi为脆性评价指标,无量纲;
Bi1为基于应力-应变曲线能量演化的能量指数,无量纲;
Bi2为考虑峰后跌落应力及峰后跌落应变大小的峰后脆性指标,无量纲;
E为杨氏弹性模量,MPa;
M为峰后弹性模量,MPa;
M=dσ/dε (2)
式中:
d为微分符号;
σ为应力,MPa;
ε为应变,无量纲;
所述能量指数的计算方法为:
式中:
dwr为峰值后的破裂能,MPa;
dwe为岩石破坏过程中被消耗的弹性能,MPa。
所述弹性能的计算方法为:
式中:
σB为峰值应力,MPa;
σC为残余应力,MPa。
对于I类岩石,所述破裂能为所述弹性能与附加能之和,所述附加能为峰后破坏过程中的弹性能不足以产生破坏而需要的额外能量;而对于II类岩石,所述破裂能为所述弹性能与释放能之差,所述释放能为峰后破坏过程中释放的多余能量;所述附加能或所述释放能的计算表达式如下:
因M=dσ/dε,故I类岩石中M<0、dwa<0,故II类岩石中M>0、dwa>0,因此所述破裂能的计算式为:
式(5)、式(6)中:dwa为附加能或释放能,MPa。
所述峰后脆性指标的计算方法为:
式中:
εC为残余应变,无量纲;
εB为峰值应变,无量纲。
在一个具体的实施例中,以官东页岩区块为例,验证所述页岩脆性评价模型的准确性。所述页岩脆性评价模型输入参数及计算结果如表1所示:
表1模型输入参数及计算结果表
根据表1可绘制如图2所示的Bi指数与岩石破裂模式关系图,根据图2可以看出,本发明所述的页岩脆性评价模型符合实验现象,能够很好地反映脆性指数与岩样破坏模式间的关系,即岩样脆性指数越大,岩样越偏向于劈裂破坏,脆性越大。
在另一个具体的实施例中,还对本发明的所述页岩脆性评价模型进行了参数敏感性分析,采用单一参数变化的方法探讨各参数对所述页岩脆性评价模型的脆性评价结果的影响。
具体的,当峰值应力作为单一参数变化时,峰值应力与脆性指数关系如图3所示,由图3可知,峰值应力与脆性指数呈明显的线性关系,随着峰值应力增大,脆性指数也随之增大,峰值应力每增加10MPa,脆性指数相应增加632左右,分析原因是因为峰值应力增大将使得峰后指数这一参数增大。
具体的,当峰值应变作为单一参数变化时,峰值应变与脆性指数关系如图4所示,由图4可知,随着峰值应变增加,脆性指数也随之增加,二者呈正相关关系,随着峰值应变增大,脆性指数也随之增大,并且增大的速率越来越快,分析其原因是因为峰值应力增大将使得峰后应变减小,从而使得峰后指数这一参数增大。
具体的,当残余应力作为单一参数变化时,残余应力与脆性指数关系如图5所示,由图5可知,残余应力与脆性指数呈负相关关系,随着残余应力增大,脆性指数先随之线性减小,到后来减小速率加快,分析其原因是因为残余应力增大使得峰后应力降幅小,从而使得峰后指数这一参数减小。
具体的,当残余应变作为单一参数变化时,残余应变与脆性指数关系如图6所示,由图6可知,残余应变与脆性指数呈负相关关系,随着残余应变增大,脆性指数随之减小,并且减小速率逐渐变小,分析原因是因为残余应变增大将导致峰后应变增大,而这将使得峰后指数这一参数减小。
具体的,当弹性模量作为单一参数变化时,弹性模量与脆性指数关系如图7所示,由图7可知,弹性模量与脆性指数呈明显的线性关系,脆性指数随弹性模量增大而增大,弹性模量每增加5000MPa,脆性指数增加7.7左右,分析其原因是因为弹性模量增大将导致能量指数增大。
综上可知,在单因素变量的情况下,页岩脆性与峰值应力、峰值应变、弹性模量成正比,与残余应力、残余应变成反比。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087573.1A CN110726608B (zh) | 2019-11-08 | 2019-11-08 | 一种基于应力-应变曲线能量演化的页岩脆性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087573.1A CN110726608B (zh) | 2019-11-08 | 2019-11-08 | 一种基于应力-应变曲线能量演化的页岩脆性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110726608A CN110726608A (zh) | 2020-01-24 |
CN110726608B true CN110726608B (zh) | 2022-04-05 |
Family
ID=69225053
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911087573.1A Active CN110726608B (zh) | 2019-11-08 | 2019-11-08 | 一种基于应力-应变曲线能量演化的页岩脆性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110726608B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111504779B (zh) * | 2020-04-07 | 2021-04-16 | 武汉大学 | 利用脆性指数确定岩石软化曲线的方法及装置 |
US12050297B2 (en) | 2020-09-11 | 2024-07-30 | Saudi Arabian Oil Company | Method and system for determining energy-based brittleness |
CN112525672B (zh) * | 2021-01-21 | 2021-05-28 | 东北石油大学 | 基于裂纹扩展能量演化的页岩储层脆性测试装置及方法 |
CN113051727B (zh) * | 2021-03-13 | 2022-12-23 | 西华大学 | 基于岩石峰前起裂及峰后应力特征的脆性评价方法及系统 |
CN115165567B (zh) * | 2022-07-08 | 2024-06-07 | 重庆交通大学 | 考虑围压及方位角的层状岩石脆性评估方法 |
CN115795916B (zh) * | 2023-01-09 | 2023-05-05 | 北京科技大学 | 一种岩石脆性评价方法及装置 |
CN117932951B (zh) * | 2024-01-30 | 2024-07-16 | 同济大学 | 基于测井数据的枯竭砂岩型储气库注采岩层可压性评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361211A (zh) * | 2014-10-24 | 2015-02-18 | 中冶长天国际工程有限责任公司 | 岩石统计损伤本构模型的构建和应用方法 |
CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105114047A (zh) * | 2015-07-22 | 2015-12-02 | 西南石油大学 | 一种利用测井资料优选体积压裂射孔层位的方法 |
CN105156103B (zh) * | 2015-09-29 | 2018-06-05 | 西南石油大学 | 一种岩屑-岩心-井眼-储层多尺度的页岩储层三维可压裂性评价方法 |
US10846447B2 (en) * | 2016-04-29 | 2020-11-24 | Exxonmobil Upstream Research Company | Method and system for stacking fracture prediction |
CN106908322B (zh) * | 2017-02-23 | 2019-04-23 | 成都理工大学 | 一种基于全应力-应变曲线的岩石脆性指数评价方法 |
CN106872260B (zh) * | 2017-03-09 | 2019-11-29 | 成都理工大学 | 一种岩石脆性指数的获取方法及岩石的脆性评价方法 |
CN108593436B (zh) * | 2018-05-11 | 2020-08-11 | 北京石油化工学院 | 一种基于应力应变曲线评价致密储层可压性的方法 |
CN110006738B (zh) * | 2019-02-28 | 2020-08-18 | 中国石油大学(北京) | 一种基于应力应变曲线和划痕测试的岩石脆性评价方法 |
CN110146374B (zh) * | 2019-03-01 | 2020-07-03 | 中国石油大学(北京) | 脆性指数的测定方法和装置 |
-
2019
- 2019-11-08 CN CN201911087573.1A patent/CN110726608B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361211A (zh) * | 2014-10-24 | 2015-02-18 | 中冶长天国际工程有限责任公司 | 岩石统计损伤本构模型的构建和应用方法 |
CN104775810A (zh) * | 2015-03-03 | 2015-07-15 | 西南石油大学 | 一种页岩气储层可压性评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110726608A (zh) | 2020-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110726608B (zh) | 一种基于应力-应变曲线能量演化的页岩脆性评价方法 | |
CN110296892B (zh) | 基于能量分析的岩石损伤演化过程中特征应力确定方法 | |
CN109632510B (zh) | 一种预测水化损伤页岩强度的方法 | |
Wu et al. | Long‐term strength determination and nonlinear creep damage constitutive model of salt rock based on multistage creep test: implications for underground natural gas storage in salt cavern | |
CN105626025A (zh) | 一种页岩储层压裂的可压性评价方法 | |
CN108827774A (zh) | 煤岩脆性评价方法 | |
US10928282B2 (en) | Method and device for determining elasticity of cement stone utilized in well cementing of oil-gas well | |
CN102128741A (zh) | 硬脆性岩石的三轴流变试验工艺与方法 | |
Petrakov et al. | The effect of fluid saturation on the elastic-plastic properties of oil reservoir rocks | |
CN104268377A (zh) | 一种基于煤岩工业组分的脆性指数确定方法 | |
Liu et al. | Experimental study of the seepage characteristics of loaded coal under true triaxial conditions | |
CN105628486A (zh) | 测量页岩岩石力学性质的方法及装置 | |
Jiang et al. | Mechanical properties investigation and damage constitutive models of red sandstone subjected to freeze-thaw cycles | |
LIU et al. | Experimental study on frost heaving pressure evolution of rock ice cracks under freezing-thawing cycles | |
Wang et al. | Evaluation of bedding effect on the bursting liability of coal and coal-rock combination under different bedding dip angles | |
CN112814669A (zh) | 一种页岩油藏全生命周期采收率预测方法和系统 | |
CN103837418A (zh) | 一种测定破裂后岩石三轴流变特性的加载路径方法 | |
Qian et al. | Experimental study on crack evolution behavior and constitutive model of granite based on the deviatoric stress to peak strength ratio | |
CN111189711B (zh) | 一种基于能量耗散的岩石屈服强度确定方法 | |
CN112505298A (zh) | 一种致密油气储产性能参数原位测试方法 | |
Chen et al. | Creep Properties of Mudstone Interlayer in Bedded Salt Rock Energy Storage Based on Multistage Creep Test: A Case Study of Huai’an Salt Mine, Jiangsu Province | |
Li et al. | Experimental study on stress sensitivity of high‐temperature and high‐pressure sandstone gas reservoirs in Yingqiong Basin | |
Xu et al. | Experimental investigation into compressive behaviour and preconsolidation pressure of structured loess at different moisture contents | |
Hou et al. | A strain rate dependent constitutive model for the Lower Silurian Longmaxi Formation shale in the fuling gas field of the Sichuan Basin, China | |
BU et al. | Experimental study on the influence of temperature and damage degree on the permeability of sandstone |
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 |