CN110208086A - 脉动压裂岩石累积损伤求解方法 - Google Patents
脉动压裂岩石累积损伤求解方法 Download PDFInfo
- Publication number
- CN110208086A CN110208086A CN201910488507.9A CN201910488507A CN110208086A CN 110208086 A CN110208086 A CN 110208086A CN 201910488507 A CN201910488507 A CN 201910488507A CN 110208086 A CN110208086 A CN 110208086A
- Authority
- CN
- China
- Prior art keywords
- pressure break
- rock
- stress
- damage
- under
- 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.)
- Granted
Links
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
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
-
- 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/0001—Type of application of the stress
- G01N2203/0005—Repeated or cyclic
-
- 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/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- 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/0058—Kind of property studied
- G01N2203/0076—Hardness, compressibility or resistance to crushing
- G01N2203/0085—Compressibility
-
- 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/0076—Hardness, compressibility or resistance to crushing
- G01N2203/0087—Resistance to crushing
-
- 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
-
- 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/025—Geometry of the test
- G01N2203/0256—Triaxial, i.e. the forces being applied along three normal axes of the specimen
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及的是脉动压裂岩石累积损伤求解方法,它包括步骤一、获取脉动压裂岩石应力‑应变滞回环曲线,进而获取一定脉动频率和应力水平下脉动压裂岩石应力‑应变滞回环曲线的耗散能和弹性能;步骤二、建立描述脉动压裂过程中岩石损伤演化计算模型;步骤三、将步骤一获取的目的压裂区块一定脉动频率和应力水平下岩石应力‑应变滞回环曲线的耗散能和弹性能,代入步骤二中的计算模型,获取岩石累积损伤变量;步骤四、重复步骤一至三,得到目的压裂区块不同脉动频率和应力水平条件下的岩石累积损伤变量;步骤五、分析不同脉动频率和应力水平条件下目的压裂区块岩石损伤演化规律。本发明可为目的压裂区块实际脉动压裂施工设计提供技术依据。
Description
技术领域:
本发明涉及的是脉动压裂岩石损伤演化的表征方法,具体涉及的是脉动压裂岩石累积损伤求解方法。
背景技术:
脉动压裂是一项非常规油气资源开发的新型技术,该技术的核心是施加周期性变化的脉动压力。在扰动应力的反复加-卸载作用下,使岩石发生疲劳损伤和破坏,一方面导致岩石强度降低,另一方面当扰动应力达到岩石相应强度极限后使岩石发生破裂形成裂缝,从而增加储层改造体积,提高油气产量。
岩石可视为一种非均质的多相复合结构,而且存在大量的天然缺陷,且这些缺陷的分布完全是随机的。岩石的破坏过程是非常复杂的,如果只是单纯地用经典弹塑性力学或断裂力学的方法来描述,将难以获得理想的结果,因此,有必要将损伤理论引入到对岩石的研究中。
损伤,即在外界荷载作用下,由细观结构缺陷萌生、扩展等不可逆变化引起的材料或结构宏观力学性能的劣化。在损伤力学范畴内,材料内部的孔隙、裂隙等天然缺陷的结构演化可反映出材料力学性能的劣化,由此将损伤与材料力学性能联系起来。材料的损伤是客观存在的,并不是一种独立的物理性质,其结构演化同时伴随着能量变化,往往通过损伤变量来描述材料损伤的优劣性。
损伤变量是岩石力学特性和结构特征的中间桥梁,也是损伤力学中有待解决的核心问题,其定义应具有涉及参数少且能够充分体现材料损伤特征的特点,既要考虑岩石实际特点,又要具有实际适用性和可行性。
目前对于岩石损伤变量的定义主要从宏观和微观两个方面:
宏观上基于弹性模量定义损伤变量,该方法实际可操作性强,但条件过于理想化,岩石中富含天然缺陷,实际获取的无损弹性模量与真实弹性模量数值仍有一定差异。
基于微缺陷面积定义损伤变量,尽管从微观角度定义的损伤变量物理意义明确,能够较为真实地反映岩石材料性能,但这些微观参数很难测量,不利于推广使用。
以往岩石损伤演化模型中累积损伤变量的表达方式多种多样,但在表征岩石损伤演化规律时都存在不同程度的不适用性,因此,需要提出适用于脉动压裂岩石累积损伤求解新方法。
发明内容:
本发明的目的是提供脉动压裂岩石累积损伤求解方法,这种脉动压裂岩石累积损伤求解方法用于解决现有技术中表征岩石损伤演化规律的方法时都存在不同程度的不适用性问题。
本发明解决其技术问题所采用的技术方案是:这种脉动压裂岩石累积损伤求解方法:
步骤一、开展目的压裂区块岩芯周期循环载荷压缩的破坏实验,通过开展一定脉动频率和应力水平下周期循环载荷压缩破坏实验,获取脉动压裂岩石应力-应变滞回环曲线,进而获取一定脉动频率和应力水平下脉动压裂岩石应力-应变滞回环曲线的耗散能和弹性能;
步骤二、根据脉动压裂岩石应力-应变曲线滞回环所反映的能量演化特征,建立描述脉动压裂过程中岩石损伤演化计算模型:
式中:D为累积损伤变量;Di为第i次脉动循环下的单次损伤变量;Udi为第i次脉动循环下的耗散能;Uei为第i次脉动循环下的弹性能;Ui为第i次脉动循环外界输入总能;
步骤三、将步骤一获取的目的压裂区块一定脉动频率和应力水平下岩石应力-应变滞回环曲线的耗散能和弹性能,代入步骤二中的脉动压裂过程中岩石损伤演化计算模型,获取岩石累积损伤变量;
步骤四、重复步骤一至步骤三,得到目的压裂区块不同脉动频率和应力水平条件下的岩石累积损伤变量;
步骤五、依据目的压裂区块岩石累积损伤变量计算结果,分析不同脉动频率和应力水平条件下目的压裂区块岩石损伤演化规律,为目的压裂区块实际脉动压裂施工设计提供技术依据。
本发明具有以下有益效果:
1、通过本发明提供的脉动压裂岩石累积损伤求解方法,可以更准确评价脉动压裂岩石损伤;
2、通过本发明提供的脉动压裂岩石累积损伤求解方法进行求解时,只需要开展周期循环载荷压缩破坏实验,简单易行,易于实现;
3、通过本发明提供的脉动压裂岩石累积损伤求解方法,能够为储层脉动压裂选取施工频率和应力水平而获取最大损伤提供依据,有利于改善压裂效果;
4、本发明提供的脉动压裂岩石累积损伤变量要根据周期循环载荷压缩破坏实验结果进行计算,应用于不同区块时,通过开展不同区块岩芯周期循环载荷压缩破坏实验来重新计算,能够确保该方法适用于不同区块,实用性强。
附图说明:
图1为周期循环载荷下岩石应力-应变滞回环曲线示意图;
图2为周期循环载荷下单个应力-应变滞回环曲线的耗散能和弹性能示意图。
具体实施方式
下面结合附图对本发明做进一步的说明:
这种脉动压裂岩石累积损伤求解方法如下:
步骤1、开展目标压裂井层的岩石现场取芯,制备周期循环载荷压缩破坏实验标准岩样,φ25*50mm或φ50*100mm岩芯,开展一定脉动频率和应力水平下周期循环载荷压缩破坏实验;
通过周期循环载荷压缩破坏实验获取应力-应变滞回环曲线(如图1所示)和每个滞回环下的耗散能和弹性能(如图2所示),图2为单个滞回环下的耗散能和弹性能示意图,其中OA为加载曲线,AB为卸载曲线,面积OAB的大小即为该滞回环下的耗散能,面积BAC的大小即为该滞回环下的弹性能,实验可直接获取各部分面积大小;
步骤2、根据脉动压裂岩石应力-应变曲线滞回环所反映的能量演化特征,建立全新的描述脉动压裂过程中岩石损伤演化的计算模型:
式中:D为累积损伤变量;Di为第i次脉动循环下的单次损伤变量;Udi为第i次脉动循环下的耗散能;Uei为第i次脉动循环下的弹性能;Ui为第i次脉动循环外界输入总能;
步骤3、将步骤一获取的目的压裂区块一定脉动频率和应力水平下岩石应力-应变滞回环曲线的耗散能和弹性能,代入步骤二中的脉动压裂过程中岩石损伤演化计算模型,获取岩石累积损伤变量;
步骤4、重复步骤1至步骤3,可计算得到目的压裂区块不同脉动频率和应力水平条件下的岩石累积损伤变量;
步骤5、依据目的压裂区块岩石累积损伤变量计算结果,分析不同脉动频率和应力水平条件下目的压裂区块岩石损伤演化规律,为实际脉动压裂施工设计提供技术依据。
Claims (1)
1.一种脉动压裂岩石累积损伤求解方法,其特征在于:
步骤一、开展目的压裂区块岩芯周期循环载荷压缩的破坏实验,通过开展一定脉动频率和应力水平下周期循环载荷压缩破坏实验,获取脉动压裂岩石应力-应变滞回环曲线,进而获取一定脉动频率和应力水平下脉动压裂岩石应力-应变滞回环曲线的耗散能和弹性能;
步骤二、根据脉动压裂岩石应力-应变曲线滞回环所反映的能量演化特征,建立描述脉动压裂过程中岩石损伤演化计算模型:
式中:D为累积损伤变量;Di为第i次脉动循环下的单次损伤变量;Udi为第i次脉动循环下的耗散能;Uei为第i次脉动循环下的弹性能;Ui为第i次脉动循环外界输入总能;
步骤三、将步骤一获取的目的压裂区块一定脉动频率和应力水平下岩石应力-应变滞回环曲线的耗散能和弹性能,代入步骤二中的脉动压裂过程中岩石损伤演化计算模型,获取岩石累积损伤变量;
步骤四、重复步骤一至步骤三,得到目的压裂区块不同脉动频率和应力水平条件下的岩石累积损伤变量;
步骤五、依据目的压裂区块岩石累积损伤变量计算结果,分析不同脉动频率和应力水平条件下目的压裂区块岩石损伤演化规律,为目的压裂区块实际脉动压裂施工设计提供技术依据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910488507.9A CN110208086B (zh) | 2019-06-05 | 2019-06-05 | 脉动压裂岩石累积损伤求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910488507.9A CN110208086B (zh) | 2019-06-05 | 2019-06-05 | 脉动压裂岩石累积损伤求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110208086A true CN110208086A (zh) | 2019-09-06 |
CN110208086B CN110208086B (zh) | 2021-07-09 |
Family
ID=67791155
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910488507.9A Active CN110208086B (zh) | 2019-06-05 | 2019-06-05 | 脉动压裂岩石累积损伤求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110208086B (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101967860A (zh) * | 2010-10-19 | 2011-02-09 | 哈尔滨工业大学 | 最不利设计地震动记录的确定方法 |
WO2013105995A2 (en) * | 2011-02-25 | 2013-07-18 | Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College | Fatigue monitoring for composite materials |
CN103344514A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法 |
CN104020259A (zh) * | 2014-06-17 | 2014-09-03 | 哈尔滨工程大学 | 阻尼材料的损耗因子与能量及温度耦合关系的测试装置及测试方法 |
CN104833536A (zh) * | 2014-02-12 | 2015-08-12 | 大连理工大学 | 一种基于非线性累积损伤理论的结构疲劳寿命计算方法 |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
CN105866373A (zh) * | 2016-03-25 | 2016-08-17 | 安徽理工大学 | 一种煤与瓦斯动力灾害的新型判别方法 |
CN106202913A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 时间相关的蠕变疲劳损伤评定方法 |
CN107103140A (zh) * | 2017-04-28 | 2017-08-29 | 电子科技大学 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
US20170296705A1 (en) * | 2015-12-03 | 2017-10-19 | The American University In Cairo | Titanium-based functional nano-architectures for drug eluting stents |
CN107505204A (zh) * | 2017-07-12 | 2017-12-22 | 河海大学 | 一种基于最小耗能原理建立岩石损伤本构模型的方法 |
CN108827774A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤岩脆性评价方法 |
CN108829993A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤层脉动水力压裂振幅和频率设计方法 |
CN108920739A (zh) * | 2018-04-27 | 2018-11-30 | 天津大学 | 一种考虑损伤累积效应的材料本构模型数值分析方法 |
CN109253928A (zh) * | 2018-10-22 | 2019-01-22 | 山东科技大学 | 一种冲击地压启动能量阈值确定方法及冲击地压预测方法 |
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
-
2019
- 2019-06-05 CN CN201910488507.9A patent/CN110208086B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101967860A (zh) * | 2010-10-19 | 2011-02-09 | 哈尔滨工业大学 | 最不利设计地震动记录的确定方法 |
WO2013105995A2 (en) * | 2011-02-25 | 2013-07-18 | Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College | Fatigue monitoring for composite materials |
CN103344514A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法 |
CN104833536A (zh) * | 2014-02-12 | 2015-08-12 | 大连理工大学 | 一种基于非线性累积损伤理论的结构疲劳寿命计算方法 |
CN104020259A (zh) * | 2014-06-17 | 2014-09-03 | 哈尔滨工程大学 | 阻尼材料的损耗因子与能量及温度耦合关系的测试装置及测试方法 |
CN105158084A (zh) * | 2015-09-15 | 2015-12-16 | 华东理工大学 | 一种材料的蠕变-疲劳寿命预测方法 |
US20170296705A1 (en) * | 2015-12-03 | 2017-10-19 | The American University In Cairo | Titanium-based functional nano-architectures for drug eluting stents |
CN105866373A (zh) * | 2016-03-25 | 2016-08-17 | 安徽理工大学 | 一种煤与瓦斯动力灾害的新型判别方法 |
CN106202913A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 时间相关的蠕变疲劳损伤评定方法 |
CN107103140A (zh) * | 2017-04-28 | 2017-08-29 | 电子科技大学 | 一种基于双线性累积损伤的时变疲劳可靠性分析方法 |
CN107505204A (zh) * | 2017-07-12 | 2017-12-22 | 河海大学 | 一种基于最小耗能原理建立岩石损伤本构模型的方法 |
CN108920739A (zh) * | 2018-04-27 | 2018-11-30 | 天津大学 | 一种考虑损伤累积效应的材料本构模型数值分析方法 |
CN108827774A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤岩脆性评价方法 |
CN108829993A (zh) * | 2018-06-23 | 2018-11-16 | 东北石油大学 | 煤层脉动水力压裂振幅和频率设计方法 |
CN109253928A (zh) * | 2018-10-22 | 2019-01-22 | 山东科技大学 | 一种冲击地压启动能量阈值确定方法及冲击地压预测方法 |
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
Non-Patent Citations (5)
Title |
---|
I. EINAV等: "《Coupled damage and plasticity models derived from energy and dissipation potentials》", 《INTERNATIONAL JOURNAL OF SOLIDS AND STRUCTURES》 * |
YUWEI LI等: "《Model of T-Type Fracture in Coal Fracturing and Analysis of Influence Factors of Fracture Morphology》", 《ENERGIES 》 * |
宋卫东等: "《三轴加-卸载下分层充填体能耗及损伤特性》", 《中国矿业大学学报》 * |
李玉伟等: "《双重孔隙渗流作用井眼围岩损伤力学模型研究》", 《数学的实践与认识》 * |
王桂林等: "《单轴压缩下非贯通节理岩体损伤破坏能量演化机制研究》", 《岩土工程学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110208086B (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104807702A (zh) | 一种煤岩体水力压裂起裂行为的实验测试方法 | |
CN104234433A (zh) | 轮辐式环向悬臂索承网格结构张拉成形方法及张拉装置 | |
Han et al. | Seismic performance of concrete-encased column base for hexagonal concrete-filled steel tube: numerical study | |
CN110208087A (zh) | 脉动压裂循环载荷下岩石强度计算方法 | |
CN110208086A (zh) | 脉动压裂岩石累积损伤求解方法 | |
CN107368630A (zh) | 一种获取耦合损耗因子的数值方法 | |
CN110823729B (zh) | 一种库区涉水厚层危岩体溃屈失稳的评价方法 | |
CN109855991A (zh) | 一种预测复合材料疲劳寿命的方法 | |
Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
CN110455650A (zh) | 一种快速确定预制裂缝混凝土梁疲劳寿命的方法 | |
CN114293964A (zh) | 一种室内评价空间应力场主动利用效果的系统及方法 | |
CN109142005B (zh) | 一种沥青混合料水稳定性评价方法 | |
Xia et al. | Ultimate strength and post ultimate strength behaviors of hull plates under extreme longitudinal cyclic load | |
Wu et al. | A novel mechanical metamaterial with dramatically different elastic parameters in tension and compression | |
CN112287445B (zh) | 一种相邻大直径桩水平动力响应分析方法及系统 | |
CN113392505A (zh) | 一种高温-水冷却循环作用后岩石统计损伤计算方法 | |
Wei et al. | Cyclic plastic behavior of dovetail under fretting load | |
CN104127173A (zh) | 一种基于动态脉搏波特征参数的病理信息提取方法 | |
CN110954419A (zh) | 一种预置裂隙水力疲劳压裂设计方法 | |
WO2022135166A1 (zh) | 确定核电厂安全壳预应力钢束张拉顺序的方法及系统 | |
Green et al. | Computation of number of equivalent strain cycles: A theoretical framework | |
Qing et al. | A constitutive model considering post-liquefaction deformation based on the logarithmic skeleton curve | |
Wati | Komunikasi Interpersonal dan Motivasi Terhadap Produktivitas Kerja Bidan Puskesmas Jatinegara | |
Han et al. | Experimental research on seismic behavior of unbonded precast reinforced concrete frame | |
CN106484952A (zh) | 一种混响场随机面压载荷模型的等效技术 |
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 |