Comparative verification of discrete and smeared numerical approaches for the simulation of hydraulic fracturing

Research output: Contribution to journalResearch articleContributedpeer-review

Contributors

  • Keita Yoshioka - , Helmholtz Centre for Environmental Research (Author)
  • Francesco Parisio - , Helmholtz Centre for Environmental Research (Author)
  • Dmitri Naumov - , Helmholtz Centre for Environmental Research (Author)
  • Renchao Lu - , Chair of Applied Environmental Systems Analysis, Helmholtz Centre for Environmental Research (Author)
  • Olaf Kolditz - , Chair of Applied Environmental Systems Analysis, Helmholtz Centre for Environmental Research (Author)
  • Thomas Nagel - , Helmholtz Centre for Environmental Research, Freiberg University of Mining and Technology (Author)

Abstract

The numerical treatment of propagating fractures as embedded discontinuities is a challenging task for which an analyst has to select a suitable numerical method from a range of options. Since their inception in the mid-80s, smeared approaches for fracture simulation such as non-local damage, gradient damage or more lately phase-field modelling have steadily gained popularity. One of the appeals of a smeared implicit fracture representation, the ability to handle complex topologies with unknown crack paths in relatively coarse meshes as well as multiple-crack interaction and multiphysics, is a fundamental requirement for the numerical simulation of hydraulic fracturing in complex situations which is technically more difficult to achieve with many other methods. However, in hydraulic fracturing simulations, not only the prediction of the fracture path but also the computation of fracture width and propagation pressure (frac pressure) is crucial for reliable and meaningful applications of the simulation tool; how to determine some of these quantities in smeared representations is not immediately obvious. In this study, two of the most popular smeared approaches of recent, namely non-local damage and phase-field models, and an approach in which the solution space is locally enriched to capture a strong discontinuity combined with a cohesive-zone model are verified against fundamental hydraulic fracture propagation problems in the toughness-dominated regime. The individual theoretical foundations of each approach are discussed and differences in the treatment of physical and numerical properties of the methods when applied to the same physical problems are highlighted through examples.

Details

Original languageEnglish
Article number13
JournalGEM - International Journal on Geomathematics
Volume10
Issue number1
Publication statusPublished - 1 Dec 2019
Peer-reviewedYes

Keywords

Keywords

  • Brittle fracture, Cohesive zone models, GEMex, GeomInt, Hydraulic fracturing, Non-local damage, OpenGeoSys, Phase field method