Indexed by CSTPCD
Scopus

Petroleum Science Bulletin ›› 2025, Vol. 10 ›› Issue (4): 719-735. doi: 10.3969/j.issn.2096-1693.2025.02.019

Previous Articles     Next Articles

FEM numerical simulation for hydraulic fracture propagation in shale reservoirs influenced by weak bedding planes

WANG Fei1,2(), LIU Wei1,2,*(), DENG Jingen1,2, LI Donggang3, TAN Yawen1,2, FENG Yongcun1,2   

  1. 1 College of Petroleum Engineering, China University of Petroleum, Beijing 102249, China
    2 State Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum, Beijing 102249, China
    3 Gas Production Company, Daqing Oilfield Company Limited, Daqing 163000, China
  • Received:2024-12-02 Revised:2025-02-17 Online:2025-08-15 Published:2025-08-05
  • Contact: LIU Wei E-mail:wangfei_cup@126.com;liuwei@cup.edu.cn

基于有限元方法的层理弱面对页岩水力裂缝扩展影响规律

王菲1,2(), 刘伟1,2,*(), 邓金根1,2, 李东刚3, 谭雅文1,2, 冯永存1,2   

  1. 1 中国石油大学(北京)石油工程学院,北京 102249
    2 中国石油大学(北京)油气资源与工程全国重点实验室,北京 102249
    3 大庆油田有限责任公司采气分公司,大庆 163000
  • 通讯作者: 刘伟 E-mail:wangfei_cup@126.com;liuwei@cup.edu.cn
  • 作者简介:王菲(1994年—),在读博士研究生,主要从事水力压裂室内实验、裂缝形态数值预测与控制技术等方面的研究,wangfei_cup@126.com
  • 基金资助:
    国家自然科学基金面上项目(52074313)

Abstract:

The Linxing gas field was selected as the research object, where weak bedding planes represent typical features and significantly influence hydraulic fracture propagation. This study provides valuable insights on hydraulic fracture propagation in bedded shale formations and offers guidance for optimizing fracturing techniques. The characteristics of shale featuring bedding planes were examined by utilizing rock mechanics experiments and direct shear tests. Considering the cementation strength and friction properties of bedding planes, a computational subroutine was developed to characterize the contact behavior for bedding planes. A 3D Finite Element Method-Cohesive Zone Model (FEM-CZM) has been established for multi-field coupling analysis of stress-damage-fluid flow, specifically incorporating bedding planes. This model incorporates a contact constitutive relationship that accounts for both friction and cementation strength of the bedding. A comprehensive and systematic quantitative analysis is conducted to investigate the influence of various factors on bedding shear slip and the propagation of hydraulic fractures. These factors include the initial opening of bedding fractures, friction coefficient, cementation strength, number of bedding planes surrounding the wellbore, and fracturing operation parameters. The results indicate that the presence of weakly bonded bedding planes leads to complex fracture propagation patterns involving both tensile and shear fractures. Bedding plane apertures serve as preferential flow pathways for fracturing fluid, significantly inhibiting fracture propagation. When the bedding aperture increases to 300 μm, the fractures are unable to cross bedding planes, which limits the fracture scale. Compared to the bonding strength, the bedding friction coefficient plays a more dominant role in determining whether fractures penetrate. Higher friction coefficients facilitate the penetration, regardless of whether the bedding planes are bonded or not. The penetration probability increases exponentially as the friction coefficient rises. In contrast, with lower friction coefficients and weakly bonded bedding planes, fractures are intercepted, while higher cementation strength allows for effective penetration. Furthermore, with a rise in the number of bedding planes, the shear fractures along these beddings expands considerably, which results in a more intricate fracture pattern. The shear failure of multiple bedding planes restricts the development of tensile-dominated fractures, which reduces the efficiency of reservoir stimulation. Optimizing fracturing fluid injection, increasing high-viscosity fracturing fluid volumes, and raising injection rates can enhance vertical fracture propagation and improve the stimulated reservoir area. Further validation of the influence of bedding planes on fracture propagation is provided by analyzing distributed temperature sensing (DTS) profiles, as well as the post-fracturing performance in Linxing.

Key words: bedding planes, 3D finite element model, hydraulic fracture propagation, frictional slip, stimulation reservoirs area

摘要:

层理弱面是页岩储层水力裂缝扩展的关键控制因素。本文以临兴气区页岩为研究对象,通过室内岩石力学实验明确该区域页岩基质和层理力学特征;建立了含层理弱面的“应力—损伤—渗流”多场耦合三维FEM-CZM模型(Finite Element Method-Cohesive Zone Model),引入考虑层理界面摩擦及胶结强度的层理接触本构,系统分析层理初始水力开度、摩擦系数、胶结强度、井周层理数目以及压裂施工参数对层理剪切滑移和裂缝扩展的影响。研究结果表明,层理干扰使水力裂缝扩展机制复杂,张拉和剪切裂缝共存。相比胶结强度,层理摩擦系数对裂缝穿层行为起主导作用;高摩擦系数促使裂缝抵抗层理干扰而更早穿层,且层理剪切缝范围减小。低摩擦系数及弱胶结强度时缝高扩展受阻;层理胶结增强会削弱层理滑移,促进裂缝穿层。层理初始水力开度达到300 μm时,压裂液沿层理渗流会导致裂缝无法穿层和改造体积受限。层理缝水力开度越大、摩擦系数越小、胶结强度越低,层理越容易剪切滑移。井周弱层理数目增多导致剪切缝面积增加,但张拉裂缝沟通面积减小。通过优化压裂液泵注程序,增加高黏压裂液用量,提高压裂液泵注排量,可以促进裂缝纵向延伸,提升储层动用程度。现场压裂井的井温监测进一步验证了层理弱面对压裂效果的影响。

关键词: 层理型页岩, 三维有限元, 水力裂缝扩展, 层理剪切滑移, 压裂改造面积

CLC Number: