煤体黏聚裂纹本构方程研究及其在压裂工程中的应用

发布时间：2020-10-14 02:15

　　 煤层气作为一种重要的清洁能源,近年来得到了广泛的关注,水力压裂是提高煤层气采收率的主要技术手段,煤层中水力压裂裂纹扩展行为将直接影响到煤层气的开采效果,因而需要对煤的断裂行为进行深入研究。线弹性断裂力学作为一种十分成功的断裂理论框架,已被广泛地应用于表征固体材料中的裂纹扩展行为。对于线弹性岩石断裂力学来说,岩石一般被简化为线弹性脆性材料,相对于固体材料试件尺寸,其裂纹尖端前断裂过程区范围很小,可以被忽略。但另一方面,煤的破坏形式通常表现为准脆性破坏,即其应力峰值后存在明显的应变软化区。对于这种准脆性材料,其断裂过程区尺寸范围相对较大,且会对材料的断裂行为产生很大的影响。因此,线弹性断裂理论已不再适用于对煤体中裂纹扩展行为的研究。而黏聚力模型被证明是一种有效的理论工具,能够描述准脆性材料断裂过程区中的非线性断裂行为。在该理论模型中,固体材料裂纹尖端前断裂过程区被简化为一条闭合的裂纹或闭合的裂纹面(分别对应二维及三维情况),其中断裂过程区内的非线性断裂行为通过黏聚力与上下裂纹面相对位移之间的本构关系进行表征。在本研究中,通过物理实验建立了不同煤阶煤的I型及I/II混合型黏聚裂纹本构方程。同时,将所建立的黏聚裂纹本构方程引入到煤体压裂裂纹扩展数值计算模型中,对煤体压裂裂纹扩展进行数值模拟。此外,对不同煤阶煤体进行了物理压裂实验,并将压裂实验结果与基于黏聚力模型的压裂数值模拟结果进行了对比。主要研究内容及结论如下:(1)通过对煤体圆盘形试件紧凑拉伸DC(T)实验确立了不同煤阶煤的I型黏聚裂纹本构方程。对于煤的I型黏聚裂纹紧凑拉伸DC(T)实验,随着煤试件煤阶的升高,其初始刚度及峰值载荷逐渐升高,最大张开位移逐渐减小,实验峰后软化阶段载荷与裂纹尖端张开位移曲线趋于线性变化,且破坏形式逐渐趋于脆性破坏;同时,随着煤化程度的提高,煤DC(T)试件的平均I型断裂能逐渐降低,断裂能实验结果变异系数值不断增加。对于较低阶煤试件来说,煤试件中I型裂纹扩展路径更加曲折,且断裂面粗糙度系数数值相对较大。更为重要的是,由黏聚裂纹应力场及位移场关系推导得到的I型黏聚裂纹本构关系的一般形式Karihaloo多项式本构方程,对于五种不同煤阶煤软化曲线的拟合度最高,且能够对其进行统一表征,因而被确定为不同煤体的I型黏聚裂纹本构方程。此外,通过对不同煤阶煤进行的I型单边缺口梁三点弯曲实验及与之相对应的数值模拟,验证了通过上述实验所建立的Karihaloo多项式黏聚裂纹本构方程对描述煤体中I型裂纹扩展行为的适用性。(2)通过对煤DC(T)试件的紧凑拉伸实验与煤PTS试件剪切贯穿实验建立了不同煤阶煤基于PPR势能函数的I/II混合型黏聚裂纹本构方程。对于煤II型黏聚裂纹的PTS试件剪切贯穿实验,随着煤试件煤阶的升高,其初始刚度及峰值载荷逐渐升高,且最大裂纹切向位移逐渐减小,同时II型黏聚裂纹断裂能逐渐降低。此外,对不同煤阶煤进行了I/II混合型单边切口梁三点弯曲实验及与之对应的引入上述I/II混合型黏聚裂纹本构方程的数值模拟,验证了通过上述实验所建立的基于PPR势能函数的I/II混合型黏聚裂纹本构方程对表征煤体中I/II混合型裂纹扩展行为的适用性。(3)对不同煤阶煤,其中包括弱粘煤、肥煤及无烟煤,进行了水力压裂物理实验。同时,对不同煤阶煤进行了液态CO_2及超临界态CO_2无水压裂实验。实验结果表明:煤试件煤阶越高,其水力压裂实验中临界起裂压力值越大;同时,起裂时间逐渐缩短。与水力压裂相比,煤体无水压裂的起裂压力值有明显降低;其中超临界态CO_2压裂的起裂压力最低,相较于水力压裂的起裂压力值,肥煤、无烟煤及泥岩的超临界态CO_2压裂的起裂压力分别降低了30.42%、33.95%及35.68%。经过无水压裂的煤岩试件,其裂纹数量明显增多,对于超临界态CO_2压裂后的煤岩试件,其中形成了相互交错的裂隙网络。(4)基于黏聚力模型对不同煤阶煤进行了水力压裂数值模拟,其中黏聚力模型中的黏聚裂纹本构关系采用本研究中通过实验所建立的煤体黏聚裂纹本构方程,模拟结果表明,采用黏聚力模型建立的数值模拟结果与水力压裂实验结果相符合,而基于线弹性断裂理论进行的煤体水力压裂数值模拟结果与水力压裂实验结果存在较大偏差。这说明黏聚力模型相比传统线弹性断裂理论更适合于研究煤体水力压裂裂纹扩展。此外,通过零厚度黏聚型单元方法实现了不同流体压裂煤层的多裂纹扩展数值模拟。
【学位单位】：太原理工大学
【学位级别】：博士
【学位年份】：2019
【中图分类】：TD712.6
【部分图文】：

并进行了现场试验研究。经过近 70 年理论上的不断完新，水压致裂技术已成功广泛地应用于油气工程[6,7],及非常规天力压裂增产煤层气工艺过程如图 1-1 所示，该技术是通过地面高一定速度连续地注入井筒，使得井底附近产生高压，当该注入的应力与储气层自身抗拉强度时，储气煤层发生破裂。随着压裂液生的裂缝持续扩展，最终在煤层井筒底部一定范围内形成裂缝网入混有支撑剂的携砂液，高压携砂液将继续促进煤层中裂缝扩展致裂的裂缝内。当携砂液泵送完毕后，通过注入低黏度破胶液，井筒返排到地面，仅流下支撑剂用于支撑储层中裂缝。水力压裂高了储层渗透率，煤层中游离的甲烷气体通过裂缝通道经井筒排，从而达到增产目的。

太原理工大学博士研究生学位论文(2) 剪切型(II 型)：裂纹受到外部平行于裂纹面的剪切应力作用，使得裂纹上下表纹面相对滑开，形成面内剪切；(3) 撕开型(III 型)：裂纹受到外部平行于裂纹面及裂纹前缘的剪切应力的共同作用裂纹上下表面相对错开，裂纹上下表面位移间断。

图 1-3 断裂行为的分类(Ba ant, 1985[147])Fig. 1-3 Types of fracture process behaviorDugdale[156]首先对裂纹前断裂过程区内裂纹扩展过程进行了研究，其对预制有穿透裂纹的金属薄板进行拉伸试验，观察到裂纹尖端前发生塑性变形形成一扁平带状区(即为断裂过程区)，Dugdale 认为在该区域内存在一种垂直于预制裂纹面的应力，只有克服该应力裂纹方可扩展，该应力被认为是材料屈服应力。同时期学者 Barenblatt[157]在对脆性材料断裂过程的研究中，提出断裂过程区内存在黏聚力限制裂纹扩展，该力来源于原子间吸引力, 吸引力的大小与裂缝张开位移有关。Hillerborg 等[158]将裂纹尖端前断裂过程区简化为一条虚拟裂纹，其上作用的黏聚力使其闭合，该黏聚力与虚拟裂纹张开度相关，裂纹扩展需克服此黏聚力，Hillerborg 将此模型成功应用于混凝土中裂纹扩展分析，由此黏聚力模型(cohesive zone mode, CZM)得以建立。黏聚力模型是将断裂过程区简化为
【相似文献】

相关期刊论文 前10条

1 沈怀荣;;弹粘塑性材料的大变形内时本构方程[J];强度与环境;1987年06期

2 Manabu Gotoh,范德鑫;一类有角点效应的塑性本构方程 Ⅰ.一般理论[J];力学进展;1988年01期

3 张培源,杨绪灿;多相固体的本构等效问题[J];重庆大学学报(自然科学版);1989年03期

4 胡强,童忠钫;非线性粘弹性材料的微分型本构方程[J];浙江大学学报(自然科学版);1989年04期

5 刘林娜,范镜泓;包含相变的非等温内时本构方程及其在2Cr13钢淬火应力分析中的应用[J];力学学报;1989年S1期

6 鲍忠兴,顾惠成,张芝婷;层压酚醛玻璃钢在高压下压缩率的研究[J];科学通报;1989年07期

7 郭兴明;非线性弹性动力学解的存在性[J];应用数学和力学;1997年03期

8 王春敏;织物本构方程的研究[J];纺织基础科学学报;1994年03期

9 阳代军;宁林新;林建国;DEAN T A;;连续冷却过程中的相变本构方程[J];机械工程学报;2009年02期

10 王少林，阮雪榆，俞新陆，陈森灿，胡宗式;金属高温塑性本构方程的研究[J];上海交通大学学报;1996年08期

相关博士学位论文 前10条

1 杨健锋;煤体黏聚裂纹本构方程研究及其在压裂工程中的应用[D];太原理工大学;2019年

2 刘健;油气井水泥石力学行为本构方程与完整性评价模型研究[D];西南石油大学;2013年

3 郭彪;铁基材料粉末锻造及致密化成形技术研究[D];西南交通大学;2012年

4 任述光;基于弹塑性本构理论的岩石相变折曲研究[D];天津大学;2008年

5 汤泽军;镁合金管材热态内压成形失稳行为研究[D];哈尔滨工业大学;2010年

6 李哲;硬质合金刀具切削高强度钢力热特性及粘结破损机理研究[D];哈尔滨理工大学;2013年

7 朱大雷;火灾下钢筋混凝土结构热弹塑性分析[D];哈尔滨工程大学;2011年

8 杨佳;聚乙烯及其共混物熔体在挤出和拉伸过程中流变行为和机理的研究[D];华南理工大学;2012年

9 徐爽;有筋扩张网免拆模板现浇外保温墙体研究[D];东北大学;2016年

10 魏国兵;镁锂合金挤压过程中的组织演变及数值模拟研究[D];重庆大学;2015年

相关硕士学位论文 前10条

1 李肖;AZ31镁合金板材在450℃变形中考虑晶粒长大及法向各向异性的本构方程[D];辽宁科技大学;2018年

2 冯启隆;张量函数表示的热弹塑性本构方程应用研究[D];太原科技大学;2018年

3 崔超;B_4C/Al复合材料本构方程及轧制特性研究[D];沈阳理工大学;2018年

4 于海富;橡胶材料本构方程的研究[D];北京化工大学;2017年

5 高小平;盐岩力学特性时温效应实验研究及其本构方程[D];中国科学院研究生院（武汉岩土力学研究所）;2005年

6 张文超;计及相变的QP980钢的动态本构方程与断裂研究[D];宁波大学;2017年

7 胡海烽;粉末高能率有限元模型压制条件实验修正与本构方程研究[D];宁波大学;2017年

8 李风;基于电磁胀环的紫铜高速率本构方程研究[D];武汉理工大学;2008年

9 王海任;各向同性与横观各向同性材料热弹塑性本构理论研究[D];太原科技大学;2016年

10 刘猛;形状记忆聚合物本构方程的构建及实验验证[D];浙江理工大学;2016年

本文编号：2840028