• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

岩石射孔开裂的初步数值模拟

林英松 姚劲松 刘莹 乔继延 丁雁生

引用本文:
Citation:

岩石射孔开裂的初步数值模拟

    作者简介: 林英松(1964- ),女,博士,教授,linyingsong@upc.edu.cn;
  • 中图分类号: O389

Preliminary numerical simulation of rock perforation cracking

  • CLC number: O389

  • 摘要: 岩石射孔作业后孔眼周围裂缝分布规律对后续压裂有不可忽视的影响。选取射孔围岩的横切面为研究对象,将三维射孔侵彻过程简化为二维扩孔过程。考虑岩石细观非均匀性,设细观强度参数服从韦布尔分布。应用拉伸破坏准则和Mohr-Coulomb压剪破坏准则,并用模量折减法处理单元开裂,从而用FEPG软件实现了有限元数值模拟。模拟结果表明:射孔后的岩石可根据裂缝产生原因及分布由内而外划分为四个区域:压剪破坏区、拉伸破坏集中区、拉伸破坏扩展区和未破坏区。分析了不同射孔弹规格及围压条件下裂纹分布变化规律。与室内模拟实验结果进行对比分析,初步验证了模型的有效性。
  • 图 1  几何模型示意图

    Figure 1.  Schematic diagram of the geometric model

    图 2  数值模拟结果图

    Figure 2.  Numerical simulation result

    图 3  岩石破坏分区图

    Figure 3.  The distribution of rock damage zones

    图 4  数值模拟裂纹数分布图

    Figure 4.  Crack number distribution of numerical simulation

    图 5  数值模拟不同缝宽裂纹数分布图

    Figure 5.  Simulated number distribution of different width crack

    图 6  侵彻孔直径7 mm围压2 MPa时岩石切片破坏形态

    Figure 6.  The rock damage form with 7 mm perforating charge under 2 MPa confining pressure

    图 7  不同侵彻孔直径下岩石切片裂纹分布图

    Figure 7.  Distribution of rock crack number with perforating charge of different diameters

    图 8  侵彻孔径9 mm围压为0时岩石切片裂缝形态

    Figure 8.  The rock damage form with 9 mm charge under no confining pressure

    图 9  不同围压下岩石切片裂纹分布图

    Figure 9.  The distribution of rock crack number under different confining pressures

    图 10  实验装置示意图

    Figure 10.  Schematic diagram of experimental device

    图 11  水泥试样安装图

    Figure 11.  Diagram of cement sample installation

    图 12  射孔后水泥靶材破坏形态

    Figure 12.  Damaged form of cement sample surface after perforation

    图 13  靶材剖面裂缝数目分布图

    Figure 13.  The distribution of crack number on cement samples

    图 14  X1-B面不同缝宽裂缝数目分布图

    Figure 14.  The distribution of crack with different width on X1-B surface

    图 15  数值模拟与物理实验靶材破坏形态对比图

    Figure 15.  Comparison of damage forms between physical experiment and numerical simulation

    图 16  数值模拟与物理实验裂缝数目分布图

    Figure 16.  Crack number distribution by physical experiment and numerical simulation

    表 1  模拟计算参数

    Table 1.  Parameters of simulation calculation

    参数
    名称
    弹性模
    量/GPa
    抗拉强
    度/MPa
    内摩擦
    角/(°)
    内聚力/
    MPa
    泊松比侵彻孔
    直径/mm
    参数大小401017200.311
    下载: 导出CSV

    表 2  天然页岩与水泥试样数据对比表

    Table 2.  Comparison of performance parameters between natural shale and cement sample

    试样种类渗透率/(10−6 μm2)孔隙度/%
    天然页岩4.04.1
    水泥试样6.73.3
    下载: 导出CSV

    表 3  水泥试样尺寸及编号表

    Table 3.  Size and number of cement samples

    围压/
    MPa
    试样直径/
    mm
    试样长度/
    mm
    试样
    编号
    上端面
    编号
    下端面
    编号
    30200200X1X1-AX1-B
    300X2X2-AX2-B
    500X3X3-AX3-B
    20200200Y1Y1-AY1-B
    300Y2Y2-AY2-B
    500Y3Y3-AY3-B
    下载: 导出CSV
  • [1] IEA. World energy outlook 2011 [R]. France: International Energy Agency, 2011.
    [2] 邹才能, 张国生, 杨智, 等. 非常规油气概念、特征、潜力及技术——兼论非常规油气地质学 [J]. 石油勘探与开发, 2013, 40(4): 385–399. DOI: 10.11698/PED.2013.04.01.
    ZOU Caineng, ZHANG Guosheng, YANG Zhi, et al. Geological concepts, characteristics, resource potential and key techniques of unconventional hydrocarbon: On unconventional petroleum geology [J]. Petroleum Exploration and Development, 2013, 40(4): 385–399. DOI: 10.11698/PED.2013.04.01.
    [3] 刘合, 王峰, 王毓才, 等. 现代油气井射孔技术发展现状与展望 [J]. 石油勘探与开发, 2014, 41(6): 731–737. DOI: 10.11698/PED.2014.06.13.
    LIU He, WANG Feng, WANG Yucai, et al. Oil well perforation technology: status and prospects [J]. Petroleum Exploration and Development, 2014, 41(6): 731–737. DOI: 10.11698/PED.2014.06.13.
    [4] 胡柳青, 李夕兵, 龚声武. 冲击载荷作用下裂纹动态响应的数值模拟 [J]. 爆炸与冲击, 2006, 26(3): 214–221. DOI: 10.11883/1001-1455(2006)03-0214-08.
    HU Liuqing, LI Xibing, GONG Shengwu. Simulation on dynamic response of crack subjected to impact loading [J]. Explosion and Shock Waves, 2006, 26(3): 214–221. DOI: 10.11883/1001-1455(2006)03-0214-08.
    [5] 朱秀星, 张林, 薛世峰, 等. 砂岩储层射孔压实伤害评价 [J]. 中国石油大学学报(自然科学版), 2014, 38(1): 137–142. DOI: 10.3969/j.issn.1673-5005.2014.01.021.
    ZHU Xiuxing, ZHANG Lin, XUE Shifeng, et al. Evaluation of perforating damage for sandstone [J]. Journal of China University of Petroleum, 2014, 38(1): 137–142. DOI: 10.3969/j.issn.1673-5005.2014.01.021.
    [6] 薛世峰, 逄铭玉, 朱秀星, 等. 砂岩储层射孔压实带孔隙度与渗透率损伤研究 [J]. 岩土力学, 2015, 36(6): 1529–1536. DOI: 10.16285/j.rsm.2015.06.002.
    XUE Shifeng, PANG Mingyu, ZHU Xiuxing, et al. Study of porosity and permeability damage of perforation compaction zone in sandstone reservoir [J]. Rock and Soil Mechanics, 2015, 36(6): 1529–1536. DOI: 10.16285/j.rsm.2015.06.002.
    [7] 单清林, 金衍, 王亚军, 等. 螺旋射孔多孔眼起裂裂缝形态有限元模拟 [J]. 中国海上油气, 2017, 29(4): 123–130. DOI: 10.11935/j.issn.1673-1506.2017.04.016.
    SHAN Qinglin, JIN Yan, WANG Yajun, et al. Finite element simulation of fracture initiation from multiple perforation channels if spiral pattern [J]. China Offshore Oil and Gas, 2017, 29(4): 123–130. DOI: 10.11935/j.issn.1673-1506.2017.04.016.
    [8] 王成, 王万军, 宁建国. 聚能装药对混凝土靶板的侵彻研究 [J]. 力学学报, 2015, 47(4): 672–886. DOI: 10.6052/0459-1879-14-336.
    WANG Cheng, WANG Wanjun, DING Jianguo. Investigation on shaped charge penetrating into concrete targets [J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(4): 672–886. DOI: 10.6052/0459-1879-14-336.
    [9] NABIPOUR A, SARMADIVALEH M. A DEM study on perforation induced damaged zones and penetration length in sandstone reservoirs [C] // 44th US Rock Mechanics Symposium and 5th US-Canada Rock Mechanics Symposium. American Rock Mechanics Association, 2010.
    [10] 向旭, 赵世华, 周伏虎. 石油射孔弹金属射流穿孔机理及金属粉末选用研究 [J]. 测井技术, 2000, 24(6): 448–449. DOI: 10.3969/j.issn.1004-1338.2000.06.012.
    XIANG Xu, ZHAO Shihua, ZHOU Fuhu. Perforation mechanism of metal perforating boreholes and selection of metal powders [J]. Well Logging Technology, 2000, 24(6): 448–449. DOI: 10.3969/j.issn.1004-1338.2000.06.012.
    [11] 范裕如. 串联外套式增效射孔技术研究[D]. 太原: 中北大学, 2013: 65.
    [12] 谭延栋. 聚能喷流射孔[M]. 北京: 石油工业出版社, 1958.
    [13] 徐芝纶. 弹性力学上册[M]. 北京: 高等教育出版社, 2006.
    [14] 林英松, 蒋金宝, 朱天玉, 等. 爆炸载荷对水泥试样损伤破坏规律研究 [J]. 中国石油大学学报(自然科学版), 2006, 30(3): 55–58. DOI: 10.3321/j.issn:1000-5870.2006.03.012.
    LIN Yingsong, JIANG Jinbao, ZHU Tianyu, et al. Research of cement sample's damage and fracture by exploding load [J]. Journal of China University of Petroleum (Natural Science Edition), 2006, 30(3): 55–58. DOI: 10.3321/j.issn:1000-5870.2006.03.012.
    [15] 荣峰. 非均匀脆性介质损伤演化的多尺度数值模拟[D]. 北京: 中国科学院力学研究所, 2006.
    [16] Recommended practices for evaluation of well perforators [S]. Washington: API Publishing Services, 2006.
    [17] 林英松, 刘莹, 魏晓菲, 等. 射孔作业对孔道周围岩石损伤规律实验研究 [J]. 实验室研究与探索, 2017, 36(12): 47–51. DOI: 10.3969/j.issn.1006-7167.2017.12.013.
    LIN Yingsong, LIU Ying, WEI Xiaofei, et al. Experimental study on failure laws of rock around perforating tunnel caused by perforation [J]. Research and Exploration in Laboratory, 2017, 36(12): 47–51. DOI: 10.3969/j.issn.1006-7167.2017.12.013.
  • [1] 曹丽娜韩秀清张永光曹宇欣 . 射孔弹爆轰波间干扰的有限元分析. 爆炸与冲击, 2011, 31(2): 220-224. doi: 10.11883/1001-1455(2011)02-0220-05
    [2] 张永康李玉龙汪海青 . 典型梁-缘结构鸟撞破坏的有限元分析. 爆炸与冲击, 2008, 28(3): 236-242. doi: 10.11883/1001-1455(2008)03-0236-07
    [3] 朱秀云潘蓉林皋李亮 . 基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析. 爆炸与冲击, 2015, 35(2): 222-228. doi: 10.11883/1001-1455(2015)02-0222-07
    [4] 李志强刘晓明赵勇刚赵隆茂 . 微爆索线性切割航空有机玻璃的实验研究和有限元分析. 爆炸与冲击, 2007, 27(5): 385-389. doi: 10.11883/1001-1455(2007)05-0385-05
    [5] 任波王乘 . 并行冲击/侵彻有限元数值模拟技术. 爆炸与冲击, 2005, 25(3): 260-264. doi: 10.11883/1001-1455(2005)03-0260-05
    [6] 张晓天贾光辉黄海 . 基于节点分离Lagrange有限元方法的超高速碰撞碎片云数值模拟. 爆炸与冲击, 2010, 30(5): 499-504. doi: 10.11883/1001-1455(2010)05-0499-06
    [7] 范志庚陈常青万强 . 泡沫铝率相关性能的有限元模拟. 爆炸与冲击, 2014, 34(6): 742-747. doi: 10.11883/1001-1455(2014)06-0742-06
    [8] 崔春生马铁华祖静 . 复合射孔压力加载下铜柱测压计测试实验和数值分析. 爆炸与冲击, 2012, 32(4): 439-444. doi: 10.11883/1001-1455(2012)04-0439-06
    [9] 韩小平张元冲沈亚鹏张泰华赵壮华 . 冲击载荷下炸药装药动态响应的有限元分析及热点形成机理的数值模拟. 爆炸与冲击, 1997, 17(2): 143-152.
    [10] 谢宗蕻卞文杰昂海松乔新 . 峰窝夹芯结构雷达罩鸟撞有限元分析与模拟. 爆炸与冲击, 1999, 19(3): 235-242.
    [11] 陈二云马大为乐贵高赵改平朱孙科 . 环形激波绕射、反射和聚焦流场的间断有限元模拟. 爆炸与冲击, 2009, 29(2): 177-181. doi: 10.11883/1001-1455(2009)02-0177-05
    [12] 王肖钧肖绍平胡秀章 . 通量修正输运法在动态有限差分和有限元中的应用. 爆炸与冲击, 1998, 18(3): 198-202.
    [13] 李海涛罗伟姜雨省张俊松 . 复合射孔爆燃气体压裂裂缝起裂扩展研究. 爆炸与冲击, 2014, 34(3): 307-314. doi: 10.11883/1001-1455(2014)03-0307-08
    [14] 谭多望陈森华刘文翰 . 有限元方法在爆炸动力学中的应用. 爆炸与冲击, 1997, 17(4): 340-346.
    [15] 孙振生张世英 . RKDG有限元方法计算流体与刚体耦合. 爆炸与冲击, 2008, 28(1): 80-85. doi: 10.11883/1001-1455(2008)01-0080-06
    [16] 王肖钧胡秀章李永池 . 高速碰撞中的有限元方法及其应用. 爆炸与冲击, 1993, 13(4): 296-304.
    [17] 宋顺成曹法和田时雨刘心德陈兆生 . 高速弹丸撞击地质材料的有限元计算. 爆炸与冲击, 1994, 14(2): 161-168.
    [18] 李妍妍郑志军虞吉林王长峰 . 闭孔泡沫金属变形模式的有限元分析. 爆炸与冲击, 2014, 34(4): 464-470. doi: 10.11883/1001-1455(2014)04-0464-07
    [19] 卫丰张光 . 药筒发射应力和抽壳力的有限元分析. 爆炸与冲击, 2001, 21(1): 72-75.
    [20] 刘晶波姚玲 . 有限元畸变单元离散模型中的SH波动. 爆炸与冲击, 1993, 13(1): 9-17.
  • 加载中
图(16)表(3)
计量
  • 文章访问数:  322
  • HTML全文浏览量:  319
  • PDF下载量:  26
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-04-12
  • 录用日期:  2018-08-31
  • 刊出日期:  2019-07-01

岩石射孔开裂的初步数值模拟

    作者简介:林英松(1964- ),女,博士,教授,linyingsong@upc.edu.cn
  • 1. 中国石油大学(华东)石油工程学院,山东 青岛 266555
  • 2. 中国科学院力学研究所,北京 100081

摘要: 岩石射孔作业后孔眼周围裂缝分布规律对后续压裂有不可忽视的影响。选取射孔围岩的横切面为研究对象,将三维射孔侵彻过程简化为二维扩孔过程。考虑岩石细观非均匀性,设细观强度参数服从韦布尔分布。应用拉伸破坏准则和Mohr-Coulomb压剪破坏准则,并用模量折减法处理单元开裂,从而用FEPG软件实现了有限元数值模拟。模拟结果表明:射孔后的岩石可根据裂缝产生原因及分布由内而外划分为四个区域:压剪破坏区、拉伸破坏集中区、拉伸破坏扩展区和未破坏区。分析了不同射孔弹规格及围压条件下裂纹分布变化规律。与室内模拟实验结果进行对比分析,初步验证了模型的有效性。

English Abstract

  • 近年来,页岩气的开发日益成为各国关注的焦点[1]。而页岩气开采的难点在于储层具有低渗透率和低孔隙度的属性,为产生工业油气流通常需要采取水力压裂等增产措施[2]。压裂作业通常是在射孔作业基础上进行的,因此射孔孔眼处岩石的损伤,对压裂作业有不可忽视的影响[3]

    关于射孔周围岩石损伤情况,仅有少数学者研究了对砂岩靶射孔后孔道周围的损伤量、应力分布及渗透率变化等,胡柳青等[4]对不同冲击载荷下裂纹响应进行数值模拟,得到了一系列随时间变化的动态应力场及动态应变场。朱秀星等[5]、薛世峰等[6]利用LS-DYNA软件模拟了射孔后砂岩骨架应力、塑性应变等参数变化。单清林等[7]建立了螺旋射孔损伤模型以研究岩体破坏后弹性参数、强度及渗透率变化规律。王成等[8]模拟了聚能射孔弹各参数变化时对漏斗坑直径及深度、侵彻孔直径及穿深等的影响。Nabipour等[9]对射孔后的砂岩进行研究,发现岩石原始孔隙度越大,射孔损伤范围越大。而针对射孔周围岩石的微裂纹分布规律的研究较为缺乏。

    本文采用FEPG有限元程序生成系统,根据射孔弹作用机理,考虑射孔作业时岩石上的载荷分布特征,建立模型进行计算。考虑岩石非均质性,将材料细观强度等参数设置为空间内随机分布,使其在载荷作用下能产生随机裂纹。模拟计算后,将结果与物理实验结果进行对比,以验证力学模型有效性,并分析不同条件下射孔孔道周围岩石裂纹形态分布规律和宽度分布规律。

    • 射孔弹金属射流侵彻岩石过程主要分3个阶段[10]

      第一阶段,射孔开始。射流以6~7 km/s的速度撞击岩石,该速度高于岩石中的声速,射流碰撞靶材时,在撞击点处产生强冲击波并向靶材内部传入,造成波后的岩石损伤,形成失效波。与此同时,高速射流在岩石表面砸出一个漏斗形坑,尺度比射流直径大几倍,仅占总孔深的很小一部分[11-12]。冲击波、失效波与开坑,构成射流侵彻的第一个阶段。

      第二阶段,后续射流继续在接触面产生极高压力,约为200 GPa,岩石产生塑性变形,侵彻孔持续加深,直径为射流几倍。此阶段持续时间长,对穿深贡献大,因与波动效应无关,称为准定常阶段[10]

      第三阶段为终止阶段,该阶段射流速度低,碰撞点压力大幅下降,逐渐趋近岩石强度,孔深停止增加。射流在孔底堆积,并在孔底前方产生压实区,远处产生开裂区。

      上述3个阶段是通常的说法,针对的是孔深和孔径。射孔周围岩石的拉伸裂纹,作为压裂问题的初始条件,常常不被注意,却是本文关注的问题。

    • 射孔过程较为复杂,建模前需简化。

      对于射流侵彻过程,第二阶段持续时间最长,射孔深度最深,产生裂纹的区域最大,因此选取此阶段产生的侵彻孔邻域内压实和远处的开裂作为研究对象。划分网格时单元尺度应与裂纹宽度保持一致,若进行三维模拟,单元和节点数量太大。为减少计算量,选择射孔围岩某一切片作为数值模拟对象,从而将侵彻过程简化为二维问题。由于二维情形射流开孔时中心点处应变率无穷大,较难模拟,初始时刻在岩石切片中心预制一个内径极小的孔,将开孔简化为扩孔。最后,为保证模拟过程接近实际,需考虑岩石非均匀性。本模型中,为保留裂缝形态和方向的随机性,不宜使用预制裂缝法,而将岩石视为细观非均匀材料,岩石的力学参数在空间中随机分布,较弱点会率先发生破坏,从而随机产生裂纹。

      综上所述,做出如下假设:

      (1)岩石切片无沿射孔轴向的位移,只有径向和周向运动;

      (2)忽略冲击波及失效波对岩石损伤的贡献;

      (3)在岩石切片中心有一个小孔,按照给定径向速度脉冲扩孔,使孔周围岩石破坏;

      (4)岩石为细观非均匀弹脆性材料,其弹性模量、抗剪强度和抗拉强度等在空间随机分布,服从韦布尔分布。

      综上所述,将三维射孔过程简化为平面应变问题;忽略冲击波效应的运动方程退化为平衡方程;细观非均匀弹脆性材料的宏观力学行为不限于弹脆性。

    • 为便于描写随机裂纹,使用直角坐标系。细观弹性本构关系为[13]

      $\left( {\begin{array}{*{20}{c}} {{\sigma _{xx}}} \\ {{\sigma _{yy}}} \\ {{\sigma _{xy}}} \end{array}} \right) = \frac{E}{{1 - {\mu ^2}}}\left( {\begin{array}{*{20}{c}} 1&\mu &0 \\ \mu &1&0 \\ 0&0&{\displaystyle\frac{{1 - \mu }}{2}} \end{array}} \right) \cdot \left( {\begin{array}{*{20}{c}} {{\varepsilon _{xx}}} \\ {{\varepsilon _{yy}}} \\ {{\varepsilon _{xy}}} \end{array}} \right)$

      式中:E为弹性模量,μ为泊松比。将式(1)简记为:

      $ {{\sigma }} = {{D}} \cdot {{\varepsilon }} $

      式中:σ为应力张量,D为弹性矩阵,ε为应变张量。

      FEPG (finite element program generator)平台适用许多偏微分方程描写的问题,但需要给出偏微分方程的弱形式。针对该问题,运用虚位移原理,由平衡方程得到虚功表达式:

      $\int_S {\left[ {\left( {\frac{{\partial {\sigma _{xx}}}}{{\partial x}} + \frac{{\partial {\sigma _{xy}}}}{{\partial y}} + {f_x}} \right){\text{δ}}u + \left( {\frac{{\partial {\sigma _{x{\rm{y}}}}}}{{\partial x}} + \frac{{\partial {\sigma _{yy}}}}{{\partial y}} + {f_y}} \right){\text{δ}}v} \right]{\rm{d}}S = 0} $

      式中:δu、δv分别为两个方向上的虚位移,fxfy分别为两个方向的质量力分量,S为岩石切片面积。利用分部积分公式将上式化为:

      $\int_S {\left( {{\sigma _{xx}}{\text{δ}}{\varepsilon _{xx}} + {\sigma _{yy}}{\text{δ}}{\varepsilon _{yy}} + {\sigma _{xy}}{\text{δ}}{\varepsilon _{xy}}} \right){\rm{d}}S} = \int_S {\left( {{f_x}{\text{δ}}u + {f_y}{\text{δ}}v} \right){\rm{d}}S + \int_\Gamma {\left( {{T_x}{\text{δ}}u + {T_y}{\text{δ}}v} \right)} } {\rm{d}}\varGamma $

      式中:Γ为岩石切片的边界,Tx、Ty分别为边界上两个方向的作用力。式(4)改写成向量形式:

      $\int_S {{\text{δ}}{{{\varepsilon }}^{\rm{T}}} \cdot {{\sigma }}{\rm{d}}S = \int_S {{\text{δ}}{{{u}}^{\rm{T}}} \cdot {{f}}{\rm{d}}S} + \int_\Gamma {{\text{δ}}{{{u}}^{\rm{T}}} \cdot {{T}}{\rm{d}}\varGamma } } $

      式中:εTuT分别为应力边界上的应变与位移,fT分别为质量力矩阵和边界力矩阵。代入本构方程式(2),得:

      $\int_S {{\text{δ}}{{{\varepsilon }}^{\rm{T}}} \cdot {{D}} \cdot {{\varepsilon }}{\rm{d}}S = } \int_S {{\text{δ}}{{{u}}^{\rm{T}}} \cdot {{f}}{\rm{d}}S} + \int_\Gamma {{\text{δ}}{{{u}}^{\rm{T}}} \cdot {{T}}{\rm{d}}\varGamma } $

      由外边界上的边界条件,可得:

      $\int_S {{\text{δ}}{{{\varepsilon }}^{{\rm{(o)}}}} \cdot {{D}} \cdot {{\varepsilon }}{\rm{d}}S = } \int_S {{\text{δ}}{{{u}}^{{\rm{(o)}}}} \cdot {{f}}{\rm{d}}S} + \int_{{\Gamma _o}} {{\text{δ}}{{{u}}^{{\rm{(o)}}}} \cdot {{{T}}^{{\rm{(o)}}}}_{}{\rm{d}}{\varGamma ^{{\rm{(o)}}}}_{}} $

      式中:上标中的o表示外边界。这就是本问题的弱形式。本模型中预制孔内径为1.5 mm,内边界给定位移$ u^{\rm{(o)}}_i$。内边界Γi的虚位移为零,使式(7)中内边界项为零。边界条件表示为:

      $\left\{ {\begin{array}{*{20}{l}} {{{T}} = {{{T}}^{{\rm{(o)}}}}\quad\;\;{\text{在}}\;{\varGamma ^{{\rm{(o)}}}}{\text{上}}}\\ {{{u}} = {{{u}}_{{i}}}\quad\,\quad{\text{在}}\;{\varGamma _i}{\text{上}}} \end{array}} \right.$

      将随侵彻过程变化的扩孔速度简化为两个阶段,在起始时间步位移量很大,后续时间步位移量很小,如:

      $ {u_r}\left( t \right) = \left\{ {\begin{array}{*{20}{l}} {0.1\Delta {u_r}}\\ {{{0.9\Delta {u_r}} / {999}}} \end{array}} \right.\begin{array}{*{20}{l}} {t = 1}\\ {1 \text{<} t \text{≤} 1\;000} \end{array} $

      式中:urt)为第t步时内边界位移量,Δur为整个扩孔过程中内边界总位移量。式(2)、(7)和(8)~(9)等,就是用FEPG软件平台生成本问题的有限元程序所需要的数学表达式。

    • 岩石的细观压剪开裂遵循Mohr-Coulomb破坏准则,细观拉伸开裂遵循拉伸破坏准则。

      岩石一旦开裂,基于连续介质建立的偏微分方程便不再成立,因此进行有限元计算时,需对裂纹进行处理。本模型中使用“模量折减法”进行处理。当单元达到破坏条件时,仍保持材料连续,折减单元模量以扩大其变形,从而反映宏观断裂。

      经过分析[14],孔眼附近岩石将发生压剪破坏。单元一旦发生剪切破坏,该单元的剪切弹性模量衰减到原值的若干分之一。孔眼远处岩石发生拉伸破坏,缝宽大于前者,因此破坏后单元的拉伸弹性模量衰减量大于前者的衰减量。模量衰减若干倍,即相同应力下单元变形扩大若干倍。模量的衰减量在计算中根据拟合图像确定。

      为了能够使模拟结果产生近似真实的随机裂纹,令弹性模量、剪切强度和拉伸强度等细观参数在空间中随机分布,分布模型选择韦布尔分布。

      为方便计算,使用一系列满足均匀分布的伪随机数通过变化得到满足韦布尔分布的随机数列A[15]

      $A = {\eta _{\rm{A}}}\exp \left\{ {\frac{{\ln \left[ { - \ln \left( {1 - F} \right)} \right]}}{{{m_{\rm{A}}}}}} \right\}$

      式中:F$[0,1)$区间内均匀分布的伪随机数列,ηAmA是关于韦布尔分布的两个参数。则岩石的弹性模量为:

      $E = {E_0}{{\rm{e}}^{{{\ln \left( { - \ln \left( {1 - F} \right)} \right)} / {{E_{\rm m}}}}}}$

      式中:$E$为满足韦布尔分布的随机弹性模量序列;E0为初始弹性模量,决定E的整体大小;Em决定E的分散程度。

      同理设置岩石抗拉强度与抗剪强度在切片中随机分布。

    • 拉伸破坏开裂后,需对缝宽进行计算。由于模型中岩石切片只受径向压力作用,因此环向应变为第一主应变,缝宽计算公式为:

      $h = c{\varepsilon _1} = c\left[ {{\varepsilon ^2} - \left( {{\varepsilon _{xx}} + {\varepsilon _{yy}}} \right)\varepsilon - \frac{1}{4}\varepsilon _{xy}^2 + {\varepsilon _{xx}}{\varepsilon _{yy}}} \right]$

      式中:h为裂缝宽度,ε1为第一主应变,c为缝宽系数。

    • 侵彻岩石切片网格划分如图1所示。

      图  1  几何模型示意图

      Figure 1.  Schematic diagram of the geometric model

    • 根据现场岩心资料及射孔资料,设置模型参数如表1。除此之外,为方便计算,韦布尔分布均匀度参数mA均取3,围压取1 MPa。

      参数
      名称
      弹性模
      量/GPa
      抗拉强
      度/MPa
      内摩擦
      角/(°)
      内聚力/
      MPa
      泊松比侵彻孔
      直径/mm
      参数大小401017200.311

      表 1  模拟计算参数

      Table 1.  Parameters of simulation calculation

    • 数值模拟计算结果如图2所示。

      图  2  数值模拟结果图

      Figure 2.  Numerical simulation result

    • 根据破坏形态,由内而外可将切片分为四块区域:首先是近孔眼处的压剪破坏区,该区域的破坏几乎全为压剪破坏,分析认为其原因是射孔时中心孔瞬间产生大变形,该区域首先达到压剪破坏条件,产生环带状压剪破坏区。第二层为拉伸破坏集中区,该区域单元径向压力不足以产生压剪破坏,但能导致足以产生拉伸破坏的环向拉力。第三层为拉伸裂纹扩展区,该区域主要是初期大变形后,在准静态载荷下稍远处单元逐渐拉伸破坏而形成的。最外层为未破坏区,随着距孔眼距离逐渐增大,单元应力逐渐降低,不足以产生破坏。切片破坏区域划分如图3所示。

      图  3  岩石破坏分区图

      Figure 3.  The distribution of rock damage zones

      记切片上半径与侵彻孔直径之比r/d0为无量纲半径,统计距射孔中心不同无量纲半径处裂纹数量和缝宽,做裂纹数量分布图和缝宽分布图如图45所示。

      图  4  数值模拟裂纹数分布图

      Figure 4.  Crack number distribution of numerical simulation

      图  5  数值模拟不同缝宽裂纹数分布图

      Figure 5.  Simulated number distribution of different width crack

    • 使用不同型号射孔弹对岩石造成的破坏不同。由于本模型中,中心孔边界条件为应变边界条件,故不同射孔弹作用载荷不同,反映在该模型中影响的是侵彻孔直径。保持其他参数不变(围压为2 MPa),分别设置侵彻孔直径d0为7、9、11、13、15 mm,分析裂纹分布的变化规律。图6为其中侵彻孔直径为7 mm时切片破坏形态:

      图  6  侵彻孔直径7 mm围压2 MPa时岩石切片破坏形态

      Figure 6.  The rock damage form with 7 mm perforating charge under 2 MPa confining pressure

      对比图2,分析可知随侵彻孔直径增大,切片压剪破坏区、拉伸破坏集中区、拉伸破坏扩展区都有不同程度的增大。其中拉伸破坏扩展区增大迅速,即对载荷大小最敏感,侵彻孔直径为13 mm和15 mm时拉伸破坏扩展区已扩大到切片边缘。裂纹数目分布如图7所示。

      图  7  不同侵彻孔直径下岩石切片裂纹分布图

      Figure 7.  Distribution of rock crack number with perforating charge of different diameters

      记压剪破坏区、拉伸破坏集中区和拉伸破坏扩展区外边界半径分别为r1r2r3,如图3所示。根据模拟结果,归纳侵彻孔直径对岩石切片裂纹数目分布在围压不变时的影响如下:

      (1)不同侵彻孔直径d0条件下,裂纹数目分布形态类似。压剪破坏区边界无量纲半径r1/d0≈1.5;无边界影响时拉伸破坏扩展区无量纲半径r3/d0<15。

      (2)裂纹数目峰值半径落在C1r/d0C2范围,称C2C1为峰值带宽。侵彻孔孔径越大,裂纹数目越多,峰值带宽越大。切片半径记作rB。当rB/d0>20时,峰值带宽C2C1很小;当rB/d0≈18时,峰值带宽上升到称C2C1≈2;当rB/d0≈15时,边界开始出现裂纹;当rB/d0≈14时,边界出现大量裂纹,峰值上界C2迅速增长。

      (3)随侵彻孔直径增大,峰值前裂纹数目增加速度及峰值后裂纹数目衰减速度加快。

      (4)侵彻孔直径较大,rB/d0接近15和小于15时,峰值后裂纹数目受边界效应的影响出现波动。

    • 保持其他参数不变(侵彻孔直径为11 mm),分别将模拟围压设为0、0.5、1.0、1.5、2.0 MPa,分析裂纹变化规律。图8为其中零围压下岩石切片破坏情况。

      图  8  侵彻孔径9 mm围压为0时岩石切片裂缝形态

      Figure 8.  The rock damage form with 9 mm charge under no confining pressure

      对比图2,总结模拟结果可得,随着围压增大,压剪破坏区范围基本不变,拉伸破坏集中区和拉伸破坏扩展区范围有不同程度的减小,其中前者减小缓慢,后者对围压变化较敏感。统计5种围压下裂缝数目分布如图9所示。

      图  9  不同围压下岩石切片裂纹分布图

      Figure 9.  The distribution of rock crack number under different confining pressures

      归纳围压对岩石切片裂纹分布规律影响如下:

      (1)围压变化基本不影响压剪破坏区和拉伸破坏集中区,峰值前裂纹数目的增长几乎相同;

      (2)围压减小时,裂纹数目的峰值有所增大,峰值带宽C2C1逐渐变大,拉伸破坏扩展区无量纲半径逐渐增长,直至达到切片边缘;

      (3)围压越小,裂纹数目越多,且峰值后裂纹数目衰减得越慢;

      (4)围压较小时,峰值后裂纹数目受边界效应的影响出现波动。

    • 因天然岩心力学性质难以重复,采用水泥石代替天然岩石进行模拟实验。为保证实验可靠性,根据天然岩心的渗透率、孔隙度等参数选择性能相近的水泥型号。经对比选择胜维G级水泥(水灰比0.44)作为试样,两者参数对比如表2

      试样种类渗透率/(10−6 μm2)孔隙度/%
      天然页岩4.04.1
      水泥试样6.73.3

      表 2  天然页岩与水泥试样数据对比表

      Table 2.  Comparison of performance parameters between natural shale and cement sample

      经实验测定,水泥试样抗拉强度为1.7 MPa,抗压强度为36 MPa。

      使用胜利油田测井公司高压射孔设备,按照API RP 19B标准[16]流程进行实验。实验装置如图10所示。由于实验设备尺寸受限,水泥试样直径最大为200 mm,为尽量避免边界效应,试样直径应尽可能大[17],因此试样直径设为200 mm。为观察横剖面裂纹分布,采用预制剖面法,将3段分别长200、300和500 mm的试样堆叠放置,如图11所示。

      图  10  实验装置示意图

      Figure 10.  Schematic diagram of experimental device

      图  11  水泥试样安装图

      Figure 11.  Diagram of cement sample installation

      此外,为方便后续工作,将试样预制剖面编号,各剖面编号名称见表3

      围压/
      MPa
      试样直径/
      mm
      试样长度/
      mm
      试样
      编号
      上端面
      编号
      下端面
      编号
      30200200X1X1-AX1-B
      300X2X2-AX2-B
      500X3X3-AX3-B
      20200200Y1Y1-AY1-B
      300Y2Y2-AY2-B
      500Y3Y3-AY3-B

      表 3  水泥试样尺寸及编号表

      Table 3.  Size and number of cement samples

    • 实验后[17]发现X1、Y1试样被射穿,X2、Y2仅侵彻一半,未被完全射穿。X1、Y1下端面裂纹分布如图12所示。

      图  12  射孔后水泥靶材破坏形态

      Figure 12.  Damaged form of cement sample surface after perforation

      用水润湿试样表面以方便肉眼观察宏观裂纹。使用便携式数码显微镜(最高放大倍数为150),结合Nano Measurer软件观察微观裂纹。以射孔孔道中心为圆心,圆心至靶材边缘最近处为半径,将试样半径16等分,分别以孔道中心为圆心做圆,孔道边缘处标为0。统计以孔眼中心,与各同心圆相交径向裂纹数量及X1-B面裂缝缝宽。统计结果如图1314所示。

      图  13  靶材剖面裂缝数目分布图

      Figure 13.  The distribution of crack number on cement samples

      图  14  X1-B面不同缝宽裂缝数目分布图

      Figure 14.  The distribution of crack with different width on X1-B surface

    • 物理模型与数值模拟岩石剖面裂纹形态对比及裂纹分布如图1516所示。

      图  15  数值模拟与物理实验靶材破坏形态对比图

      Figure 15.  Comparison of damage forms between physical experiment and numerical simulation

      图  16  数值模拟与物理实验裂缝数目分布图

      Figure 16.  Crack number distribution by physical experiment and numerical simulation

      对比数值模型与物理实验参数设置,发现前者围压相较于后者小一个数量级时,裂缝数目分布才较相似。初步分析知,径向裂纹始于拉伸破坏,它达到一定长度后的扩展取决于断裂机制,后者在较小拉伸应力作用下可以继续延伸。本文的力学模型尚未引入断裂机制,因此导致这一偏差。

      分析可知物理实验与数值模拟裂纹分布较相似,中心一圈裂纹较少,外围出现大量裂纹并呈辐射状扩展,后逐渐衰减。由于尺寸有限,物理实验靶材剖面不存在未破坏区,且靶材受反射波影响,边界效应明显,剖面中存在几条的环向裂纹,而模拟计算时忽略爆炸冲击波且尺寸较大,没有此类问题。统计两结果裂纹数目作图16。可以看出了两结果中裂缝数目分布规律趋势类似。

    • (1)实验结果初步验证了所建立的射流侵彻岩石二维简化力学模型的有效性,并提示应当增加断裂机制。

      (2)根据数值模拟计算结果,以破坏类型及裂缝形态数量为依据,可将靶材切片由内而外分为4个区域:压剪破坏区、拉伸破坏集中区、拉伸破坏扩展区和未破坏区。

      (3)随着侵彻孔孔径的增大及围压的降低,压剪破坏区、拉伸破坏集中区及拉伸破坏扩展区范围均有不同程度的增大,其中拉伸破坏扩展区最为敏感,增长最快;且侵彻孔孔径改变时,裂纹数目无量纲半径上的分布趋势大致相同。

参考文献 (17)

目录

    /

    返回文章
    返回