Finite element modelling of viscosity-dominated hydraulic fractures

• Published in: Journal of petroleum science & engineering 2012-06, Vol.88-89, p.136-144
• Main Author: Chen, Zuorong
• Format: Article
• Language: English
• Subjects: cohesive zone model; Computational fluid dynamics; Crack propagation; Finite element method; Fluid flow; Fracture mechanics; hydraulic fracture; Hydraulic fracturing; Hydraulics; Mathematical analysis; Mathematical models; mesh transition; petroleum geomechanics
• ISSN: 0920-4105; 1873-4715
• DOI: 10.1016/j.petrol.2011.12.021

Hydraulic fracturing is a highly effective technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to developing more efficient application of this technology, and also presents a significant technical challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on existing pore pressure cohesive finite elements has been established to simulate the propagation of a viscosity-dominated hydraulic fracture in an infinite, impermeable elastic medium. Selected results of the finite element modelling and comparisons with analytical solutions are presented for viscosity-dominated plane strain and penny-shaped hydraulic fractures, respectively. Some important issues such as mesh transition and far-field boundary approximation in the cohesive finite element model have been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by fluid viscosity demonstrates the capability of the cohesive zone finite element model in simulating the hydraulic fracture growth. ► Cohesive finite element modelling of hydraulic fractures. ► Viscosity-dominated plane strain and penny-shaped hydraulic fractures. ► Excellent agreement between the finite element results and analytical solutions. ► Mesh transition and node sharing provides a higher solution accuracy and efficiency. ► Determine the far-field boundary deformation for fractures in an infinite domain.