Search

In reservoir simulation, it is important to understand the mechanical behaviour of fractured rocks and the effect of shear and tensile displacements of fractures on their aperture. Tensile opening directly enhances the fracture aperture, whereas shear of a preexisting rough-walled fracture creates aperture changes dependent on the local stress state. Since fracture dilatation increases reservoir permeability, both processes must be included in a realistic and consistent manner into the mechanical reservoir simulation model. Here, we use the extended finite volume method (XFVM) to conduct flow and geomechanics simulations. In XFVM, fractures are embedded in a poroelastic matrix and are modelled with discontinuous basis functions. On each fracture segment the tractions and compressive forces are calculated, and one extra degree of freedom is added for both the shear and tensile displacement. In this particular XFVM implementation we assume that linear elasticity and steady state fluid pressure adequately constrain the effective stress. In this paper, shear dilation is not calculated a posteriori, but it enters the equations such that aperture changes directly affect the stress state. This is accomplished by adding shear dilation to the displacement gradients and therefore ascertains a consistent representation in the stress-strain relations and force balances. We illustrate and discuss the influence of this extra term in two simple test cases and in a realistic layer-restricted two-dimensional fracture network subjected to plausible in situ stress and pore pressure conditions. [ABSTRACT FROM AUTHOR]

We present a fracture-only reservoir simulator for multiphase flow: the fracture geometry is modeled explicitly, while fluid movement between fracture and matrix is accommodated using empirical transfer functions. This is a hybrid between discrete fracture discrete matrix modeling where both the fracture and matrix are gridded and dual-porosity or dual-permeability simulation where both fracture and matrix continua are upscaled. The advantage of this approach is that the complex fracture geometry that controls the main flow paths is retained. The use of transfer functions, however, simplifies meshing and makes the simulation method considerably more efficient than discrete fracture discrete matrix models. The transfer functions accommodate capillary- and gravity-mediated flow between fracture and matrix and have been shown to be accurate for simple fracture geometries, capturing both the early- and late-time average behavior. We verify our simulator by comparing its predictions with simulation results where the fracture and matrix are explicitly modeled. We then show the utility of the approach by simulating multiphase flow in a geologically realistic fracture network. Waterflooding runs reveal the fraction of the fracture–matrix interface area that is infiltrated by water so that matrix imbibition can occur. The evolving fraction of the fracture–matrix interface area turns out to be an important characteristic of any particular fracture system to be used as a scaling parameter for capillary driven fracture–matrix transfer. [ABSTRACT FROM AUTHOR]