A continuum model for coupled stress and fluid flow in discrete fracture networks

Quan Gan, Derek Elsworth

Research output: Contribution to journalArticle

30 Citations (Scopus)
4 Downloads (Pure)

Abstract

We present a model coupling stress and fluid flow in a discontinuous fractured mass represented as a continuum by coupling the continuum simulator TF_FLAC3D with cell-by-cell discontinuum laws for deformation and flow. Both equivalent medium crack stiffness and permeability tensor approaches are employed to characterize pre-existing discrete fractures. The advantage of this approach is that it allows the creation of fracture networks within the reservoir without any dependence on fracture geometry or gridding. The model is validated against thermal depletion around a single stressed fracture embedded within an infinite porous medium that cuts multiple grid blocks. Comparison of the evolution of aperture against the results from other simulators confirms the veracity of the incorporated constitutive model, accommodating stress-dependent aperture under different stress states, including normal closure, shear dilation, and for fracture walls out of contact under tensile loading. An induced thermal unloading effect is apparent under cold injection that yields a larger aperture and permeability than during conditions of isothermal injection. The model is applied to a discrete fracture network to follow the evolution of fracture permeability due to the influence of stress state (mean and deviatoric) and fracture orientation. Normal closure of the fracture system is the dominant mechanism where the mean stress is augmented at constant stress obliquity ratio of 0.65—resulting in a reduction in permeability. Conversely, for varied stress obliquity (0.65–2) shear deformation is the principal mechanism resulting in an increase in permeability. Fractures aligned sub-parallel to the major principal stress are near-critically stressed and have the greatest propensity to slip, dilate and increase permeability. Those normal to direction of the principal stress are compacted and reduce the permeability. These mechanisms increase the anisotropy of permeability in the rock mass. Furthermore, as the network becomes progressively more sparse, the loss of connectivity results in a reduction in permeability with zones of elevated pressure locked close to the injector—with the potential for elevated pressures and elevated levels of induced seismicity.
Original languageEnglish
Pages (from-to)43-61
Number of pages19
JournalGeomechanics and Geophysics for Geo-Energy and Geo-Resources
Volume2
Issue number1
Early online date5 Jan 2016
DOIs
Publication statusPublished - Mar 2016

Fingerprint

fracture network
fluid flow
permeability
obliquity
simulator
fracture orientation
induced seismicity
fracture geometry
dilation
unloading
porous medium
stiffness
connectivity
crack
anisotropy

Keywords

  • coupled simulation
  • discrete fracture network
  • geothermal reservoir
  • fracture permeability
  • stress-dependent aperture

Cite this

A continuum model for coupled stress and fluid flow in discrete fracture networks. / Gan, Quan; Elsworth, Derek.

In: Geomechanics and Geophysics for Geo-Energy and Geo-Resources, Vol. 2, No. 1, 03.2016, p. 43-61.

Research output: Contribution to journalArticle

@article{b3a17de4dc4a4e9d93cdd883143b7b9c,
title = "A continuum model for coupled stress and fluid flow in discrete fracture networks",
abstract = "We present a model coupling stress and fluid flow in a discontinuous fractured mass represented as a continuum by coupling the continuum simulator TF_FLAC3D with cell-by-cell discontinuum laws for deformation and flow. Both equivalent medium crack stiffness and permeability tensor approaches are employed to characterize pre-existing discrete fractures. The advantage of this approach is that it allows the creation of fracture networks within the reservoir without any dependence on fracture geometry or gridding. The model is validated against thermal depletion around a single stressed fracture embedded within an infinite porous medium that cuts multiple grid blocks. Comparison of the evolution of aperture against the results from other simulators confirms the veracity of the incorporated constitutive model, accommodating stress-dependent aperture under different stress states, including normal closure, shear dilation, and for fracture walls out of contact under tensile loading. An induced thermal unloading effect is apparent under cold injection that yields a larger aperture and permeability than during conditions of isothermal injection. The model is applied to a discrete fracture network to follow the evolution of fracture permeability due to the influence of stress state (mean and deviatoric) and fracture orientation. Normal closure of the fracture system is the dominant mechanism where the mean stress is augmented at constant stress obliquity ratio of 0.65—resulting in a reduction in permeability. Conversely, for varied stress obliquity (0.65–2) shear deformation is the principal mechanism resulting in an increase in permeability. Fractures aligned sub-parallel to the major principal stress are near-critically stressed and have the greatest propensity to slip, dilate and increase permeability. Those normal to direction of the principal stress are compacted and reduce the permeability. These mechanisms increase the anisotropy of permeability in the rock mass. Furthermore, as the network becomes progressively more sparse, the loss of connectivity results in a reduction in permeability with zones of elevated pressure locked close to the injector—with the potential for elevated pressures and elevated levels of induced seismicity.",
keywords = "coupled simulation, discrete fracture network, geothermal reservoir, fracture permeability, stress-dependent aperture",
author = "Quan Gan and Derek Elsworth",
note = "Acknowledgments This work is a partial result of support from the US Department of Energy under project DOE-DE-343 EE0002761. This support is gratefully acknowledged. We also acknowledge the data from the University of Oklahoma and University of Texas at Austin.",
year = "2016",
month = "3",
doi = "10.1007/s40948-015-0020-0",
language = "English",
volume = "2",
pages = "43--61",
journal = "Geomechanics and Geophysics for Geo-Energy and Geo-Resources",
issn = "2363-8419",
publisher = "Springer International Publishing AG",
number = "1",

}

TY - JOUR

T1 - A continuum model for coupled stress and fluid flow in discrete fracture networks

AU - Gan, Quan

AU - Elsworth, Derek

N1 - Acknowledgments This work is a partial result of support from the US Department of Energy under project DOE-DE-343 EE0002761. This support is gratefully acknowledged. We also acknowledge the data from the University of Oklahoma and University of Texas at Austin.

PY - 2016/3

Y1 - 2016/3

N2 - We present a model coupling stress and fluid flow in a discontinuous fractured mass represented as a continuum by coupling the continuum simulator TF_FLAC3D with cell-by-cell discontinuum laws for deformation and flow. Both equivalent medium crack stiffness and permeability tensor approaches are employed to characterize pre-existing discrete fractures. The advantage of this approach is that it allows the creation of fracture networks within the reservoir without any dependence on fracture geometry or gridding. The model is validated against thermal depletion around a single stressed fracture embedded within an infinite porous medium that cuts multiple grid blocks. Comparison of the evolution of aperture against the results from other simulators confirms the veracity of the incorporated constitutive model, accommodating stress-dependent aperture under different stress states, including normal closure, shear dilation, and for fracture walls out of contact under tensile loading. An induced thermal unloading effect is apparent under cold injection that yields a larger aperture and permeability than during conditions of isothermal injection. The model is applied to a discrete fracture network to follow the evolution of fracture permeability due to the influence of stress state (mean and deviatoric) and fracture orientation. Normal closure of the fracture system is the dominant mechanism where the mean stress is augmented at constant stress obliquity ratio of 0.65—resulting in a reduction in permeability. Conversely, for varied stress obliquity (0.65–2) shear deformation is the principal mechanism resulting in an increase in permeability. Fractures aligned sub-parallel to the major principal stress are near-critically stressed and have the greatest propensity to slip, dilate and increase permeability. Those normal to direction of the principal stress are compacted and reduce the permeability. These mechanisms increase the anisotropy of permeability in the rock mass. Furthermore, as the network becomes progressively more sparse, the loss of connectivity results in a reduction in permeability with zones of elevated pressure locked close to the injector—with the potential for elevated pressures and elevated levels of induced seismicity.

AB - We present a model coupling stress and fluid flow in a discontinuous fractured mass represented as a continuum by coupling the continuum simulator TF_FLAC3D with cell-by-cell discontinuum laws for deformation and flow. Both equivalent medium crack stiffness and permeability tensor approaches are employed to characterize pre-existing discrete fractures. The advantage of this approach is that it allows the creation of fracture networks within the reservoir without any dependence on fracture geometry or gridding. The model is validated against thermal depletion around a single stressed fracture embedded within an infinite porous medium that cuts multiple grid blocks. Comparison of the evolution of aperture against the results from other simulators confirms the veracity of the incorporated constitutive model, accommodating stress-dependent aperture under different stress states, including normal closure, shear dilation, and for fracture walls out of contact under tensile loading. An induced thermal unloading effect is apparent under cold injection that yields a larger aperture and permeability than during conditions of isothermal injection. The model is applied to a discrete fracture network to follow the evolution of fracture permeability due to the influence of stress state (mean and deviatoric) and fracture orientation. Normal closure of the fracture system is the dominant mechanism where the mean stress is augmented at constant stress obliquity ratio of 0.65—resulting in a reduction in permeability. Conversely, for varied stress obliquity (0.65–2) shear deformation is the principal mechanism resulting in an increase in permeability. Fractures aligned sub-parallel to the major principal stress are near-critically stressed and have the greatest propensity to slip, dilate and increase permeability. Those normal to direction of the principal stress are compacted and reduce the permeability. These mechanisms increase the anisotropy of permeability in the rock mass. Furthermore, as the network becomes progressively more sparse, the loss of connectivity results in a reduction in permeability with zones of elevated pressure locked close to the injector—with the potential for elevated pressures and elevated levels of induced seismicity.

KW - coupled simulation

KW - discrete fracture network

KW - geothermal reservoir

KW - fracture permeability

KW - stress-dependent aperture

U2 - 10.1007/s40948-015-0020-0

DO - 10.1007/s40948-015-0020-0

M3 - Article

VL - 2

SP - 43

EP - 61

JO - Geomechanics and Geophysics for Geo-Energy and Geo-Resources

JF - Geomechanics and Geophysics for Geo-Energy and Geo-Resources

SN - 2363-8419

IS - 1

ER -